Characterization of dynamics and reactivities of solvated ions byab initio simulations

The stability of [Zn(NH(3))(4)](2+) in water: A quantum mechanical/molecular mechanical molecular dynamics study

To investigate the structural and dynamical properties of the tetraamminezinc(ii) complex (Zn-tetraamine) in aqueous solution, ab initio quantum mechanical/molecular mechanical (QM/MM) molecular dynamics (MD) simulation was performed for 50 ps at the Hartree-Fock (HF) level of theory. A predominant 4-coordinate solvation structure with a maximum probability of the Zn-N distance at approximately 2.1 A was observed, which seems to be involved in the associative mode of water exchange reactions to produce a short-lived, 5-coordinated trigonal bipyramidal structure. Several sets of structural and dynamical parameters such as radial distribution functions (RDF), coordination number distributions (CND), angular distributions (ADF), ligands' mean residence times (MRT) and ion-ligand stretching frequencies have been evaluated in order to get an in depth knowledge of the physical and chemical properties of the tetraamminezinc(ii) complex in aqueous solution. A comparative study of the tetraamminezinc(ii) complex with previously published mono-, di- and triamminezinc(ii) complexes has been also performed, which yielded significant insights into the complex properties as a function of an increasing number of first-shell ammonia ligands.

Structural and dynamical aspects of the unsymmetric hydration of Sb(III): an ab initio quantum mechanical charge field molecular dynamics simulation

An ab initio quantum mechanical charge field molecular dynamics (QMCF MD) simulation was performed to investigate the behavior of the Sb(3+) ion in aqueous solution. The simulation reveals a significant influence of the residual valence shell electron density on the solvation structure and dynamics of Sb(3+). A strong hemidirectional behavior of the ligand binding pattern is observed for the first hydration shell extending up to the second hydration layer. The apparent domain partitioned structural behavior was probed by solvent reorientational kinetics and three-body distribution functions. The three-dimensional hydration space was conveniently segmented such that domains having different properties were properly resolved. The approach afforded a fair isolation of localized solvent structural and dynamical motifs that Sb(3+) seems to induce to a remarkable degree. Most intriguing is the apparent impact of the lone pair electrons on the second hydration shell, which offers insight into the mechanistic aspects of hydrogen bonding networks in water. Such electronic effects observed in the hydration of Sb(3+) can only be studied by applying a suitable quantum mechanical treatment including first and second hydration shell as provided by the QMCF ansatz.

Temperature effects on the structure and dynamics of the Jahn-Teller distorted Cr(2+) ion in aqueous solution: a hybrid QM/MM molecular dynamics simulation

The influences of an elevated temperature on the structure and dynamics of the Jahn-Teller distorted [Cr(H(2)O)(6)](2+) complex have been studied using an ab initio quantum mechanical/molecular mechanical (QM/MM) molecular dynamics simulation, showing that the increased temperature affects the lifetime distortions of the hydrated Cr(2+) ion by decreasing the inversion time to 0.5-2 ps.

Structure and dynamics of the hydrated magnesium ion and of the solvated magnesium carbonates: insights from first principles simulations

We report first principles molecular dynamics simulations based on the density functional theory and the Car-Parrinello method to study the structures and dynamics of the hydrated Mg(2+) ion and of the solvated MgHCO(3)(+) and MgCO(3) complexes in aqueous solution. According to these simulations, the first hydration shell of the hydrated magnesium ion consists of six water molecules, whereas in the solvated magnesium bicarbonate and magnesium carbonate complexes the Mg(2+) is mostly five-coordinated, which indicates that when coordinated to magnesium the HCO(3)(-) and CO(3)(2-) anions reduce its the coordination sphere. Our simulations show that the structures of the most stable monomers of magnesium bi-carbonate and magnesium carbonate in solution are Mg[eta(1)-HCO(3)](H(2)O)(4)(+) and Mg[eta(1)-CO(3)](H(2)O)(4), i.e. the preferred hydration number is four, while the (bi-)carbonate is coordinated to the magnesium in a monodentate mode. The analysis of the exchange processes of the water molecules in the first and second hydration shell of Mg(2+) shows that the HCO(3)(-) or CO(3)(2-) ligands affect the dynamics of the magnesium coordination spheres by making its hydration shell more "labile". Furthermore, molecular dynamics simulations of the non-associated Mg(2+)/Cl(-) pair in water suggest that, despite negligible differences in the coordination spheres of Mg(2+), the chloride anion has a significant influence on the water exchange rates in the second hydration shell of Mg(2+).

Structure and dynamics of the U4+ ion in aqueous solution: an ab initio quantum mechanical charge field molecular dynamics study

The structure and dynamics of the stable four-times positively charged uranium(IV) cation in aqueous solution have been investigated by ab initio quantum mechanical charge field (QMCF) molecular dynamics (MD) simulation at the Hartree-Fock double-zeta quantum mechanical level. The QMCF-MD approach enables investigations with the accuracy of a quantum mechanics/molecular mechanics approach without the need for the construction of solute-solvent potentials. Angular distribution functions; radial distribution functions; coordination numbers of the first, second, and third shell (9, 19, and 44, respectively); coordination number distribution functions; tilt- and Theta-angle distribution functions; as well as local density corrected triangle distribution functions have been employed for the evaluation of the hydrated ion's structure. Special attention was paid to the determination of the geometry of the first hydration layer, and the results were compared to experimental large-angle X-ray scattering and extended X-ray absorption fine structure data. The solvent dynamics around the ion were also investigated using mean ligand residence times and related data and, resulting from the unavailability of any experimental data, were compared to ions with similar properties.

Structure and dynamics of phosphate ion in aqueous solution: an ab initio QMCF MD study

A simulation of phosphate in aqueous solution was carried out employing the new QMCF MD approach which offers the possibility to investigate composite systems with the accuracy of a QMMM method but without the time consuming creation of solute-solvent potential functions. The data of the simulations give a clear picture of the hydration shells of the phosphate anion. The first shell consists of 13 water molecules and each oxygen of the phosphate forms in average three hydrogens bonds to different solvent molecules. Several structural parameters such as radial distribution functions and coordination number distributions allow to fully characterize the embedding of the highly charged phosphate ion in the solvent water. The dynamics of the hydration structure of phosphate are described by mean residence times of the solvent molecules in the first hydration shell and the water exchange rate.

Structural and dynamical properties of hydrogen fluoride in aqueous solution: an ab initio quantum mechanical charge field molecular dynamics simulation

The novel ab initio quantum mechanical charge field (QMCF) molecular dynamics simulation at the Hartree-Fock level has been employed to investigate hydration structure and dynamics of hydrogen fluoride in aqueous solution. The average H-F bond length of 0.93 A obtained from the QMCF MD simulation is in good agreement with the experimental data. The HHF...Ow distance of 1.62 A was evaluated for the first hydration shell, and 2.00 A was observed for the FHF...Hw distance. The stability of hydrogen bonding is more pronounced in the hydrogen site of hydrogen fluoride, with a single water molecule in this part of the first hydration shell. A wide range of coordination numbers between 3 and 9 with an average value of 5.6 was obtained for the fluorine site. The force constants of 819.1 and 5.9 N/m were obtained for the HHF-FHF and HHF...Ow interactions, respectively, proving the stability of the nondissociated form of hydrogen fluoride in aqueous solution. The mean residence times of 2.1 and 2.5 ps were determined for ligand exchange processes in the neighborhood of fluorine and hydrogen atoms of hydrogen fluoride, respectively, indicating a weak structure-making effect of hydrogen fluoride in water. The corresponding H-bond lifetimes attribute this effect to the H atom site of HF.

Al(III) hydration revisited. An ab initio quantum mechanical charge field molecular dynamics study

To assess the novel quantum mechanical charge field (QMCF) molecular dynamics (MD) approach, two simulations of hydrated Al(III) have been carried out, as this system proved to be a well-suited test case for hybrid ab initio/molecular mechanics simulations. Two different population analysis schemes according to Mulliken and Lowdin have been applied to evaluate the atomic charges in the QM region. It is shown that the QMCF MD approach yields a substantially improved description of the system and that, due to the fact that solute-solvent potentials can be renounced, the QMCF MD framework is a more convenient approach to investigate solvated systems compared to conventional ab initio QM/MM MD approaches.

Capabilities of chemical simulation methods in the elucidation of structure and dynamics of solutions

As a result of recent methodological developments in connection with enhanced computational capacity, theoretical methods have become increasingly valuable and reliable tools for the investigation of solutions. Simulation techniques utilizing a quantum mechanical (QM) approach for the treatment of the chemically most relevant region so-called hybrid quantum mechanical/molecular mechanical (QM/MM) simulations have reached a level of accuracy that often equals or may even surpass experimental methods. The latter is true in particular whenever ultrafast (i.e., picosecond) dynamics prevail, such as in labile hydrates or structure-breaking systems. The recent development of an improved QM/MM framework, the quantum mechanical charge field (QMCF) ansatz, enables a broad spectrum of solute systems to be elucidated. As this novel methodology does not require any solute solvent potential functions, the applicability of the QMCF method is straightforward and universal. This advantage is bought, however, at the price of a substantial increase of the QM subregion, and an attendant increase in computational periods to levels of months, and even a year, despite parallelizing high-performance computing (HPC) clusters. Molecular dynamics (MD) simulations of chemical systems showing increasing complexity have been performed, and demonstrate the superiority of the QMCF ansatz over conventional QM/MM schemes. The systems studied include Pd2+, Pt2+, and Hg22+, as well as composite anions such as PO43- and ClO4-.

Hydration Behavior of Alkyl Amines and Their Corresponding Protonated Forms. 1. Ammonia and Methylamine

The structure and dynamics of hydration of ammonia/ammonium and methylamine/methylammonium systems have been studied by Car-Parrinello molecular dynamics simulation. While methylamine interacts weakly with the aqueous environment, the interaction of ammonia is found to be much stronger than expected. Both protonated species show a highly structured first solvation sphere. The solvent exchange mechanisms for all species were also investigated, along with the geometry of the hydration spheres. Comparison of these exchange mechanisms with that published for the ammonium ion shows only a minor difference. Analysis of the respective distribution functions has allowed insight into the thermodynamics of solvation for both systems. The calculated pKa values (9.23/10.65) correspond very closely with the published experimental values of 9.25 and 10.65.

Simulation of Ca2+ and Mg2+ solvation using polarizable atomic multipole potential

The alkaline earth metals calcium and magnesium are critically involved in many biomolecular processes. To understand the hydration thermodynamics of these ions, we have performed molecular dynamics simulations using a polarizable potential. Particle-mesh Ewald for point multipoles has been applied to the calculation of electrostatic interactions. The parameters in this model have been determined from an ab initio quantum mechanical calculation of dimer interactions between ions and water. Two methods for ion solvation free energy calculation, free energy perturbation, and the Bennett acceptance ratio have been compared. Both predict results consistent with other theoretical estimations while the Bennett approach leads to a much smaller statistical error. Based on the Born theory and the ion-oxygen radial distribution functions, we estimate the effective size of the ions in solution, concluding that K(+) > Na(+) congruent with Ca(2+) > Mg(2+). There appears to be much stronger perturbation in water structure, dynamics, and dipole moment around the divalent cations than the monovalent K(+) and Na(+). The average water coordination numbers for Ca(2+) and Mg(2+) are 7.3 and 6, respectively. The lifetime of water molecules in the first solvation shell of Mg(2+) is on the order of hundreds of picoseconds, in contrast to only few picoseconds for Ca(2+), K(+), or Na(+).

Ab initio QM/MM dynamics of H3O+ in water

A molecular dynamics (MD) simulation based on a combined ab initio quantum mechanics/molecular mechanics (QM/MM) method has been performed to investigate the solvation structure and dynamics of H3O+ in water. The QM region is a sphere around the central H3O+ ion, and contains about 6-8 water molecules. It is treated at the Hartree-Fock (HF) level, while the rest of the system is described by means of classical pair potentials. The Eigen complex (H9O4+) is found to be the most prevalent species in the aqueous solution, partly due to the selection scheme of the center of the QM region. The QM/MM results show that the Eigen complex frequently converts back and forth into the Zundel (H5O2+) structure. Besides the three nearest-neighbor water molecules directly hydrogen-bonded to H3O+, other neighbor waters, such as a fourth water molecule which interacts preferentially with the oxygen atom of the hydronium ion, are found occasionally near the ion. Analyses of the water exchange processes and the mean residence times of water molecules in the ion's hydration shell indicate that such next-nearest neighbor water molecules participate in the rearrangement of the hydrogen bond network during fluctuative formation of the Zundel ion and, thus, contribute to the Grotthuss transport of the proton.

Pair interaction potentials with explicit polarization for molecular dynamics simulations of La3+ in bulk water

Pair interaction potentials (IPs) were defined to describe the La(3+)-OH(2) interaction for simulating the La(3+) hydration in aqueous solution. La(3+)-OH(2) IPs are taken from the literature or parametrized essentially to reproduce ab initio calculations at the second-order Moller-Plesset level of theory on La(H(2)O)(8) (3+). The IPs are compared and used with molecular dynamics (MD) including explicit polarization, periodic boundary conditions of La(H(2)O)(216) (3+) boxes, and TIP3P water model modified to include explicit polarization. As expected, explicit polarization is crucial for obtaining both correct La-O distances (r(La-O)) and La(3+) coordination number (CN). Including polarization also modifies hydration structure up to the second hydration shell and decreases the number of water exchanges between the La(3+) first and second hydration shells. r(La-O) ((1))=2.52 A and CN((1))=9.02 are obtained here for our best potential. These values are in good agreement with experimental data. The tested La-O IPs appear to essentially account for the La-O short distance repulsion. As a consequence, we propose that most of the multibody effects are correctly described by the explicit polarization contributions even in the first La(3+) hydration shell. The MD simulation results are slightly improved by adding a-typically negative 1r(6)-slightly attractive contribution to the-typically exponential-repulsive term of the La-O IP. Mean residence times are obtained from MD simulations for a water molecule in the first (1082 ps) and second (7.6 ps) hydration shells of La(3+). The corresponding water exchange is a concerted mechanism: a water molecule leaving La(H(2)O)(9) (3+) in the opposite direction to the incoming water molecule. La(H(2)O)(9) (3+) has a slightly distorded "6+3" tricapped trigonal prism D(3h) structure, and the weakest bonding is in the medium triangle, where water exchanges take place.

Stability of different zinc(II)-diamine complexes in aqueous solution with respect to structure and dynamics: a QM/MM MD study

In the context of our detailed study of the chemical behavior of aquo- and ammine-Zn(II) complexes, ab initio quantum mechanical/molecular mechanical (QM/MM) molecular dynamics (MD) simulations were performed at the Hartree-Fock (HF) level for the zinc(II)-diamine complexes in aqueous solution. The initial structures of cis and trans isomers of the tetraaquodiamminezinc(II) complex were found to transform into the triaquodiamminezinc(II) complex by releasing one water ligand after approximately 6 and approximately 22 ps of simulation time, respectively. The structural and dynamical properties of these three zinc complexes, i.e., cis-[Zn(NH3)2(H2O)4]2+, trans-[Zn(NH3)2(H2O)4]2+, and [Zn(NH3)2(H2O)3]2+, were analyzed in terms of radial distribution functions (RDF), coordination number distributions (CND), angular distribution functions (ADF), tilt and theta angle distributions, ligands' mean residence times (MRTs), and ion-ligand stretching frequencies. One considerably elongated Zn-O bond of 2.43 A was observed in the case of the cis isomer for one of the water ligands located in the trans position to an ammonia ligand. In the trans isomer the average Zn-O bond length was observed to be 2.23 A, while in the triaquodiamminezinc(II) complex two distinct Zn-O bonds, namely 2.12 A for the ligands in the trigonal plane and 2.26 A for axial water molecules, were observed. As both of the octahedral isomers are transformed into the pentacoordinated structure within the picosecond range, they might be regarded as "metastable species or intermediates", while the triaquodiamminezinc(II) complex is the most stable species of the zinc(II)-diamine complex in aqueous solution.

Structure and Dynamics of Sulfate Ion in Aqueous SolutionAn ab initio QMCF MD Simulation and Large Angle X-ray Scattering Study

The hydrated sulfate ion has been characterized in aqueous solution in structural and dynamic aspects using ab initio quantum mechanical charge field (QMCF) molecular dynamics (MD) simulation and large angle X-ray scattering (LAXS) methods. The LAXS data show an average coordination number of the sulfate ion of up to 12 water molecules bound through hydrogen bonding, while the QMCF MD simulation displays a wide range of coordination numbers between 8 and 14 with an average value of approximately 11. The Os...Ow distance cannot be distinguished from the Ow...Ow distance in the LAXS experiment; the weighted mean O...O distance is 2.880(10) A. In the simulation, the Os...Ow and Ow...Ow distances are found to be very similar, namely, 2.86 and 2.84 A, respectively. The S-Os bond and S...Ow distance have been determined by the LAXS experiment as 1.495(6) and 3.61(2) A, respectively, indicating an average nearly tetrahedral S-Os...Ow angle. The approximately 5% deviations of simulation distances (1.47 and 3.82 A) from the experimental ones can probably be ascribed to the neglect of correlation energy in the quantum mechanical method. The mean residence time of water ligands at O atoms, 2.57 ps, is longer than that in pure water, 1.7 ps, characterizing the sulfate ion as a weak structure maker.

Adaptive Partitioning in Combined Quantum Mechanical and Molecular Mechanical Calculations of Potential Energy Functions for Multiscale Simulations

In many applications of multilevel/multiscale methods, an active zone must be modeled by a high-level electronic structure method, while a larger environmental zone can be safely modeled by a lower-level electronic structure method, molecular mechanics, or an analytic potential energy function. In some cases though, the active zone must be redefined as a function of simulation time. Examples include a reactive moiety diffusing through a liquid or solid, a dislocation propagating through a material, or solvent molecules in a second coordination sphere (which is environmental) exchanging with solvent molecules in an active first coordination shell. In this article, we present a procedure for combining the levels smoothly and efficiently in such systems in which atoms or groups of atoms move between high-level and low-level zones. The method dynamically partitions the system into the high-level and low-level zones and, unlike previous algorithms, removes all discontinuities in the potential energy and force whenever atoms or groups of atoms cross boundaries and change zones. The new adaptive partitioning (AP) method is compared to Rode's "hot spot" method and Morokuma's "ONIOM-XS" method that were designed for multilevel molecular dynamics (MD) simulations. MD simulations in the microcanonical ensemble show that the AP method conserves both total energy and momentum, while the ONIOM-XS method fails to conserve total energy and the hot spot method fails to conserve both total energy and momentum. Two versions of the AP method are presented, one scaling as O(2N) and one with linear scaling in N, where N is the number of groups in a buffer zone separating the active high-level zone from the environmental low-level zone. The AP method is also extended to systems with multiple high-level zones to allow, for example, the study of ions and counterions in solution using the multilevel approach.