SCC-DFTB parameters for simulating hybrid gold-thiolates compounds.

We present a parametrization of a self-consistent charge density functional-based tight-binding scheme (SCC-DFTB) to describe gold-organic hybrid systems by adding new Au-X (X = Au, H, C, S, N, O) parameters to a previous set designed for organic molecules. With the aim of describing gold-thiolates systems within the DFTB framework, the resulting parameters are successively compared with density functional theory (DFT) data for the description of Au bulk, Aun gold clusters (n = 2, 4, 8, 20), and Aun SCH3 (n = 3 and 25) molecular-sized models. The geometrical, energetic, and electronic parameters obtained at the SCC-DFTB level for the small Au3 SCH3 gold-thiolate compound compare very well with DFT results, and prove that the different binding situations of the sulfur atom on gold are correctly described with the current parameters. For a larger gold-thiolate model, Au25 SCH3 , the electronic density of states and the potential energy surfaces resulting from the chemisorption of the molecule on the gold aggregate obtained with the new SCC-DFTB parameters are also in good agreement with DFT results.

Parametrization of the SCC-DFTB Method for Halogens.

Parametrization of the approximative DFT method SCC-DFTB for halogen elements is presented. The new parameter set is intended to describe halogenated organic as well as inorganic molecules, and it is compatible with the established parametrization of SCC-DFTB for carbon, hydrogen, oxygen, and nitrogen. The performance of the parameter set is tested on a representative set of molecules and discussed.

Proton transport in functionalised additives for PEM fuel cells: contributions from atomistic simulations.

The conventional polymer electrolyte membrane (PEM) materials for fuel cell applications strongly rely on temperature and pressure conditions for optimal performance. In order to expand the range of operating conditions of these conventional PEM materials, mesoporous functionalised SiO(2) additives are developed. It has been demonstrated that these additives themselves achieve proton conductivities approaching those of conventional materials. However, the proton conduction mechanisms and especially factors influencing charge carrier mobility under different hydration conditions are not well known and difficult to separate from concentration effects in experiments. This tutorial review highlights contributions of atomistic computer simulations to the basic understanding and eventual design of these materials. Some basic introduction to the theoretical and computational framework is provided to introduce the reader to the field, the techniques are in principle applicable to a wide range of other situations as well. Simulation results are directly compared to experimental data as far as possible.

Kinetic aspects of the thermostatted growth of ice from supercooled water in simulations.

In experiments, the growth rate of ice from supercooled water is seen to increase with the degree of supercooling, that is, the lower the temperature, the faster the crystallization takes place. In molecular dynamics simulations of the freezing process, however, the temperature is usually kept constant by means of a thermostat that artificially removes the heat released during the crystallization by scaling the velocities of the particles. This direct removal of energy from the system replaces a more realistic heat-conduction mechanism and is believed to be responsible for the curious observation that the thermostatted ice growth proceeds fastest near the melting point and more slowly at lower temperatures, which is exactly opposite to the experimental findings [M. A. Carignano, P. B. Shepson, and I. Szleifer, Mol. Phys. 103, 2957 (2005)]. This trend is explained by the diffusion and the reorientation of molecules in the liquid becoming the rate-determining steps for the crystal growth, both of which are slower at low temperatures. Yet, for a different set of simulations, a kinetic behavior analogous to the experimental finding has been reported [H. Nada and Y. Furukawa, J. Crystal Growth 283, 242 (2005)]. To clarify this apparent contradiction, we perform relatively long simulations of the TIP4P/Ice model in an extended range of temperatures. The temperature dependence of the thermostatted ice growth is seen to be more complex than was previously reported: The crystallization process is very slow close to the melting point at 270 K, where the thermodynamic driving force for the phase transition is weak. On lowering the temperature, the growth rate initially increases, but displays a maximum near 260 K. At even lower temperatures, the freezing process slows down again due to the reduced diffusivity in the liquid. The velocity of the thermostatted melting process, in contrast, shows a monotonic increase upon raising the temperature beyond the normal melting point. In this case, the effects of the increasing thermodynamic driving force and the faster diffusion at higher temperatures reinforce each other. In the context of this study, we also report data for the diffusion coefficient as a function of temperature for the water models TIP4P/Ice and TIP4P/2005.

Proton conductivity of SO3 H-functionalized benzene-periodic mesoporous organosilica.

The proton conductivity of benzene-periodic mesoporous silica (PMO) materials functionalized with sulfonic acid groups is investigated using experimental and theoretical techniques. The SO(3) H functionalization of pristine benzene-PMO is realized by three different pathways based on a grafting method in which surface silanol groups and/or benzene rings are used to anchor SO(3) H groups for enhanced proton conductivity. The functionalized material is experimentally characterized using X-ray diffraction, small-angle neutron scattering, and argon adsorption isotherms. After pressing the functionalized benzene-PMOs into pellets, the proton conductivity is deduced from Bode plots of impedance spectra taken in the temperature range of 333-413 K at 100% relative humidity. Using quantum mechanical approaches for selected proton-conduction mechanisms, the free energy barriers for proton transport as well as the local water environment at the surface are calculated. These calculations indicate that different mechanisms from purely bulk water transport are important for the benzene-PMO proton conduction, in agreement with experimental data.

Relativistic parametrization of the self-consistent-charge density-functional tight-binding method. 1. Atomic wave functions and energies.

A detailed treatment of a confined relativistic atom, needed as an initial step for the parametrization of the self-consistent-charge density-functional tight-binding method, is presented and discussed. The required one-component quantities, i.e., orbital energies, orbital wave functions, and Hubbard parameters, are obtained by weighted averaging of the corresponding numbers determined for the atomic spinors. The wave function and density confinement is achieved by introducing the Woods-Saxon potential in the atomic four-component Dirac-Kohn-Sham problem. The effect of the additional confining potential on energy eigenvalues and the shape of atomic wave functions and densities is discussed and numerical examples are presented for the valence spinors of carbon, germanium, and lead.

Treatment of collinear and noncollinear electron spin within an approximate density functional based method.

We report benchmark calculations of the density functional based tight-binding method concerning the magnetic properties of small iron clusters (Fe2 to Fe5) and the Fe13 icosahedron. Energetics and stability with respect to changes of cluster geometry of collinear and noncollinear spin configurations are in good agreement with ab initio results. The inclusion of spin-orbit coupling has been tested for the iron dimer.

Parameter Calibration of Transition-Metal Elements for the Spin-Polarized Self-Consistent-Charge Density-Functional Tight-Binding (DFTB) Method: Sc, Ti, Fe, Co, and Ni.

Recently developed parameters for five first-row transition-metal elements (M = Sc, Ti, Fe, Co, and Ni) in combination with H, C, N, and O as well as the same metal (M-M) for the spin-polarized self-consistent-charge density-functional tight-binding (DFTB) method have been calibrated. To test their performance a couple sets of compounds have been selected to represent a variety of interactions and bonding schemes that occur frequently in transition-metal containing systems. The results show that the DFTB method with the present parameters in most cases reproduces structural properties very well, but the bond energies and the relative energies of different spin states only qualitatively compared to the B3LYP/SDD+6-31G(d) density functional (DFT) results. An application to the ONIOM(DFT:DFTB) indicates that DFTB works well as the low level method for the ONIOM calculation.

Computational studies on polymer adhesion at the surface of gamma-Al2O3. I. The adsorption of adhesive component molecules from the gas phase.

We calculate the minimum energy paths and reaction energies of the adsorption of the epoxide adhesive components diglycidylesterbisphenol A (DGEBA), diethyltriamine (DETA), and the adhesion promoter 3-aminopropylmethoxysilane (AMEO) at two different sites on a model of the native Al2O3 surface, using the nudged elastic band algorithm in conjunction with self-consistent charge-density functional based tight binding. Our results show that the chosen combination of methods is well suited to obtain an overview of the reaction mechanisms and kinetics of the adsorption of organic molecules on inorganic surfaces. The obtained MEP-s show that there is preference for the adsorption of the adhesion promoter, AMEO, over the resin, DGEBA, while the adsorption of the curing agent, DETA, is unfavorable. Our approach also gives an insight into the ranges of the mechanical and electronic influences of the adsorption process on the interface, which neither full ab initio methods nor force field approaches can provide. These results will help to develop a quantum mechanics-molecular mechanics multiscale embedding scheme for more detailed studies of organic/inorganic hybrid interface reactions.