Ab initio quantum mechanical charge field molecular dynamics simulation (QMCF-MD) of Bi(3+) in water.

The hydration of the Bi(III) ion was determined via an ab initio quantum mechanical charge field molecular dynamics (QMCF-MD) simulation. Ten picosecond sampling was carried out to determine structural and dynamical properties of the Bi(III) ion in aqueous solution. In the first hydration shell, the ion is 9-fold coordinated with a maximum probability of the Bi-O distance at 2.51 Å. In total, 11 exchanges were observed in the first-shell showing associative, dissociative, and interexchange character. As with the dominant existence of 9-fold coordination, the geometry of the Bi(III) ion is in between the tricapped trigonal prism and the capped square antiprism.

A QMCF-MD investigation of the structure and dynamics of Ce4+ in aqueous solution.

A quantum-mechanical charge-field molecular dynamics simulation has been performed for a tetravalent Ce ion in aqueous solution. In this framework, the complete first and second hydration spheres are treated by ab initio quantum mechanics supplemented by an electrostatic embedding technique, making the construction of non-Coulombic solute-solvent potentials unnecessary. During the 10 ps of simulation time, the structural aspects of the solution were analyzed by various methods. Experimental results such as the mean Ce-O bond distance and the predicted first-shell coordination number were compared to the results obtained from the simulation resolving some ambiguities in the literature. The dynamics of the system were characterized by mean ligand residence times and frequency/force constant calculations. Furthermore, Ce-O and Ce-H angular radial distribution plots were employed, yielding deeper insight into the structural and dynamical aspects of the system.

Guanidinium in aqueous solution studied by quantum mechanical charge field-molecular dynamics (QMCF-MD).

Structure and dynamics of guanidinium in aqueous solution were examined via a double zeta HF level Quantum Mechanical Charge Field-Molecular Dynamics (QMCF-MD) simulation, as well as two Molecular Mechanics-Molecular Dynamics (MM-MD) simulations, parametrised via the Amber99 parameter set, employing the side chain of arginine as a template. Coulombic parameters were fitted via Mulliken population data of the QM simulation, as well as via the recommended restrained electrostatic potential fit (RESP). Although guanidinium is one of the most weakly hydrated cations yet characterised, its hydration pattern is quite complex and pronounced in the quantum mechanical simulation. The positive charge is mainly located on the central carbon, resulting in strong solute-oxygen coordination. Hydrogen bonds are mainly donated by the amide hydrogens, but are also accepted via the nitrogens to a very low extent. Detailed analysis of structure and dynamics, comparing the applied QM and MM models, provides evidence that the arginine parametrisation leads to highly different results than the quantum mechanical treatment, and that the RESP parametrisation is too polarised.

Combining a Dissociative Water Model with a Hybrid QM/MM Approach-A Simulation Strategy for the Study of Proton Transfer Reactions in Solution.

An implementation strategy for a dissociative water potential in hybrid QM/MM simulations is outlined. As the knowledge of the time-dependent topology is crucial for the assignment of solvent molecules to the QM or MM subregion, proton transfer events and the associated change of the molecular composition have to be monitored as the simulation is progressing. A simple and effective update criterion is proposed, which was found to be an efficient tool to identify sustained proton transfer reactions. The resulting topology data enable the application of the dissociative solvent model in QM/MM simulations and serve as a reference for the analysis of time-dependent properties such as the proton hopping rate and the diffusion coefficient. For the latter, an interpolation scheme is proposed linking subsequent proton transfer events into a single diffusive entity. Suitable settings of key parameters for the topology update and the interpolation scheme have been determined by analyzing MD trajectories of an excess proton in water and 1 M HCl. The resulting values for the proton hopping rate and the diffusion coefficient are well within the range estimated by EVB models and CPMD approaches. An investigation of the hydrolytic conversion of As(III) to [As(OH)2](+) serves as an exemplary application of the dissociative model in a QM/MM simulation study.

Hydration of trivalent lanthanum revisited - An ab initio QMCF-MD approach.

The previously investigated La3+-hydrate has been re-evaluated by means of the quantum mechanical charge field (QMCF) molecular dynamics (MD) approach. Improved description of the hydration characteristics has been realised by including the full second hydration shell into the quantum mechanically treated region and by introducing the influence of the surrounding bulk via an electrostatic embedding technique. Analytical tools such as the ligand angular radial distribution analysis have been employed to gain deeper insight into the structural features of the hydrate. La3+ simultaneously forms nona- and decahydrates with capped trigonal and quadratic prismatic structure, besides small amounts of an octahydrate.

Computational study of the cerium(III) ion in aqueous environment.

This work comprises the first quantum chemical simulation study of the Ce3+ ion in aqueous environment. The structural and dynamical properties have been investigated by means of the quantum mechanical charge field (QMCF) molecular dynamics (MD) approach and the results, where applicable, have been compared to experimental data. Besides conventional analytical tools, angular radial distribution functions have been employed to gain deeper insight into the structure of the hydrate. The ion-oxygen stretching motion's wavenumber, further characterising the Ce-O bond, is in excellent agreement with experimental results, same as the structural values obtained from the simulation.

Hydrogen bond formation of formamide and N-methylformamide in aqueous solution studied by quantum mechanical charge field-molecular dynamics (QMCF-MD).

The formation of hydrogen-bonds of formamide and "cis"-N-methylformamide in aqueous solution was examined using double zeta HF level Quantum Mechanical Charge Field-Molecular Dynamics (QMCF-MD) simulations. Basic attributes such as structure and dynamics of the solvates and hydrogen-bonds were studied in particular by means of coordination number distributions, mean residence times and radial distribution functions, on which spatial restrictions in the form of planes and cones were applied. Advanced methods of analysis gave detailed information about the sterical environment and the dynamic behavior of strong and weak hydrogen-bonds formed by the residues. The comparison of both molecules over a sampling period of 12 ps provided information on the influence of methylation of the amide function on molecular and hydration properties.

Sulfur dioxide in water: structure and dynamics studied by an ab initio quantum mechanical charge field molecular dynamics simulation.

An ab initio Quantum Mechanical Charge Field Molecular Dynamics Simulation (QMCF MD) was performed to investigate structure and dynamics behavior of hydrated sulfur dioxide (SO(2)) at the Hartree-Fock level of theory employing Dunning DZP basis sets for solute and solvent molecules. The intramolecular structural characteristics of SO(2), such as S═O bond lengths and O═S═O bond angle, are in good agreement with the data available from a number of different experiments. The structural features of the hydrated SO(2) were primarily evaluated in the form of S-O(wat) and O(SO(2))-H(wat) radial distribution functions (RDFs) which gave mean distances of 2.9 and 2.2 Å, respectively. The dynamical behavior characterizes the solute molecule to have structure making properties in aqueous solution or water aerosols, where the hydrated SO(2) can easily get oxidized to form a number of sulfur(VI) species, which are believed to play an important role in the atmospheric processes.

Structure and dynamics of methanol in water: a quantum mechanical charge field molecular dynamics study.

An ab initio quantum mechanical charge field molecular dynamics simulation was carried out for one methanol molecule in water to analyze the structure and dynamics of hydrophobic and hydrophilic groups. It is found that water molecules around the methyl group form a cage-like structure whereas the hydroxyl group acts as both hydrogen bond donor and acceptor, thus forming several hydrogen bonds with water molecules. The dynamic analyses correlate well with the structural data, evaluated by means of radial distribution functions, angular distribution functions, and coordination number distributions. The overall ligand mean residence time, τ identifies the methanol molecule as structure maker. The relative dynamics data of hydrogen bonds between hydroxyl of methanol and water molecules prove the existence of both strong and weak hydrogen bonds. The results obtained from the simulation are in excellent agreement with the experimental results for dilute solution of CH(3)OH in water. The overall hydration shell of methanol consists in average of 18 water molecules out of which three are hydrogen bonded.

Hydration of highly charged ions.

Based on a series of ab initio quantum mechanical charge field molecular dynamics (QMCF MD) simulations, the broad spectrum of structural and dynamical properties of hydrates of trivalent and tetravalent ions is presented, ranging from extreme inertness to immediate hydrolysis. Main group and transition metal ions representative for different parts of the periodic system are treated, as are 2 threefold negatively charged anions. The results show that simple predictions of the properties of the hydrates appear impossible and that an accurate quantum mechanical simulation in cooperation with sophisticated experimental investigations seems the only way to obtain conclusive results.

Structure and dynamics of the Zr(4+) ion in water.

The four-times positively charged zirconium ion in aqueous solution was simulated, using an ab initio quantum mechanical charge field molecular dynamics approach. As no hydrolysis reaction occurred during the simulation time of 10 ps, the target of this study was the evaluation of the structure and dynamics of the monomeric hydrated zirconium(iv) ion. The ion forms three hydration shells. In the first hydration shell the ion is 8-fold coordinated with a maximum probability of the Zr-O distance at 2.25 Å. While no exchanges occurred between the first and second shell, the mean residence time of the water molecules in the second shell is 5.5 ps. A geometry of the first hydration shell in-between a bi-capped trigonal prism and a square antiprism was found and a Zr-O force constant of 188 N m(-1) was evaluated.

Hydrolysis of tetravalent group IV metal ions: an ab initio simulation study.

Simulations of the tetravalent group IV metal ions, Ge(IV), Sn(IV), and Pb(IV), in aqueous solution were performed using ab initio quantum mechanical charge field molecular dynamics (QMCF MD). The process of hydrolysis, which occurred for each of the metal ions, was analyzed in terms of the time evolution of solvent configuration. Several important factors involved in the initiation of proton dissociation from first shell water molecules are discussed in connection to the nature of proton mobility in aqueous solution.

Structure and dynamics of the chromate ion in aqueous solution. An ab initio QMCF-MD simulation.

An ab initio quantum-mechanical charge-field molecular-dynamics (QMCF-MD) simulation of the chromate ion in aqueous solution at ambient temperature was performed to study the structure and dynamics of this ion and its hydration shell. In contrast to conventional quantum-mechanical molecular-mechanics molecular-dynamics (QM/MM-MD) simulations, the QMCF-MD approach offers the possibility of investigating composite systems with the accuracy of a QM/MM method but without the time-consuming construction of solute-solvent potential functions. The data of the simulation give a clear picture of the first hydration shell of the chromate anion, which consists of 14 water molecules. The mean distance between the oxygen atoms of the chromate and the hydrogen atoms of water is 1.82 A. Each chromate oxygen atom is in average coordinated to 2.6 water molecules. The first-shell mean ligand residence time was evalulated as 2.2 ps; the vibrational frequency of the nu(OH) mode was found to be 185 cm(-1). Several structural parameters such as the radial distribution functions, angular distribution functions, and coordination number distributions enable a full characterization of the embedding of the chromate ion in the solvent water. The dynamics of the hydration structure are described by mean residence times of the water molecules in the first hydration shell, distance plots, and velocity autocorrelation functions.

Structure and dynamics of the UO2+ ion in aqueous solution: an ab initio QMCF-MD study.

The quantum mechanical charge field molecular dynamics (QMCF-MD) framework was applied in a simulation of the uranyl(v) ion in aqueous solution. The structure was evaluated on the basis of overall and sectorial radial distribution functions, angular distribution functions, tilt- and Theta-angle distribution functions and coordination number distributions. The cation is strongly coordinated by 4 water ligands at an average distance of 2.51 A, while the oxygen atoms are on average bound to 1.2 water molecules at a distance of 2.9 A. A mean residence time of 2.83 ps was evaluated for the oxygen sites of the uranyl(v) ion. The results are in good agreement with previous experimental and theoretical data on the hydration of similar ions.

A quantum mechanical charge field molecular dynamics study of Fe(2+) and Fe(3+) ions in aqueous solutions.

Ab initio quantum mechanical charge field molecular dynamics (QMCF-MD) simulations have been performed for aqueous solutions of Fe(2+) and Fe(3+) ions at the Hartree-Fock level of theory to describe and compare their structural and dynamical behavior. The structural features of both hydrated ions are characterized by radial distribution functions that give the maximum probability of the ion-O distance for Fe(2+) and Fe(3+) ions at 2.15 and 2.03 A, respectively. The angular distribution functions of both ions prove the octahedral arrangement of six water ligands, whereas the second shells of these ions differ. Both ions show influence on the water molecules beyond the second shells. The structure-forming abilities of both ions are visible from the ligand mean residence times and ion-O stretching frequencies evaluated for both ions. The substantially improved data obtained from these QMCF-MD simulations show better correlation with available experimental results than the conventional quantum mechanics/molecular mechanics molecular dynamics (QM/MM MD) approaches with one hydration shell treated by quantum mechanics.

Hydrolytic conversion of AsO(-3)4 to HAsO(-2)4: a QMCF MD study.

A quantum mechanical charge field molecular dynamics (QMCF MD) study of AsO in water was carried out to gain insight into its conversion from the hydrated anion resulting in OH(-) ions and HAsO, which occurs on the scale of a few hundred femtoseconds. The OH(-) ion undergoes further proton exchange with water molecules, while HAsO is a stable species.

Structural and dynamic aspects of hydration of HAsO4(-2): an ab initio QMCF MD simulation.

An ab initio quantum mechanical charge field simulation has been carried out in order to obtain molecular level insight into the hydration behavior of HAsO4(-2), one of the major biologically active components of As(V) oxoanion in neutral to slightly alkaline aqueous medium. Moreover, a geometrical definition of hydrogen bonding has been used to probe and characterize both solute-solvent and solvent-solvent hydrogen bonding present in the system. The asymmetry of the anion induced by the protonation of one of the oxygens of the arsenate anion causes rather irregular hydration structure. The nonprotonated oxygen atoms preferably form relatively stable hydrogen bonds with two to three water molecules in their vicinity, while the protonated oxygen forms one or two hydrogen bonds, weaker than water-water hydrogen bonds. The two types of As-O distances obtained from the simulation (1.68 and 1.78 A for the protonated and nonprotonated oxygens, respectively) are in good agreement with the experimental data. The two types of As-O stretching frequencies obtained from the simulation (855 and 660 cm(-1) reproduce well the experimental ATR-FTIR results (859 and 680-700 cm(-1)).

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 dependence of structure and dynamics of the hydrated Ca2+ ion according to ab initio quantum mechanical charge field and classical molecular dynamics.

Simulations using ab initio quantum mechanical charge field molecular dynamics (QMCF MD) and classical molecular dynamics using two-body and three-body potentials were performed to investigate the hydration of the Ca(2+) ion at different temperatures. Results from the simulations demonstrate significant effects of temperature on solution dynamics and the corresponding composition and structure of hydrated Ca(2+). Substantial increase in ligand exchange events was observed in going from 273.15 K to 368.15 K, resulting in a redistribution of coordination numbers to lower values. The effect of temperature is also visible in a red-shift of the ion-oxygen stretching frequencies, reflecting weakened ligand binding. Even the moderate increase from ambient to body temperature leads to significant changes in the properties of Ca(2+) in aqueous environment.

Hydrated germanium (II): irregular structural and dynamical properties revealed by a quantum mechanical charge field molecular dynamics study.

Structural and dynamical properties of Ge (II) in aqueous solution have been investigated using the novel ab initio quantum mechanical charge field (QMCF) molecular dynamics (MD) formalism. The first and second hydration shells were treated by ab initio quantum mechanics at restricted Hartree-Fock (RHF) level using the cc-pVDZ-PP basis set for Ge (II) and Dunning double-zeta plus polarization basis sets for O and H. Besides ligand exchange processes and mean ligand residence times to observe dynamics, tilt- and theta-angle distributions along with an advanced structural parameter, namely radial and angular distribution functions (RAD) for different regions were also evaluated. The combined radial and angular distribution depicted through surface plot and contour map is presented to provide a detailed insight into the density distribution of water molecules around the Ge(2+) ion. A strongly distorted hydration structure with two trigonal pyramidal substructures within the first hydration shell is observed, which demonstrates the lone-pair influence and provides a new basis for the interpretation of the catalytic and pharmacological properties of germanium coordination compounds.