An ab initio quantum mechanical charge field molecular dynamics simulation of a dilute aqueous HCl solution

An ab initio quantum mechanical charge field (QMCF) molecular dynamics simulation has been performed to study the structural and dynamical properties of a dilute aqueous HCl solution. The solute molecule HCl and its surrounding water molecules were treated at Hartree‐Fock level in conjunction with Dunning double‐ζ plus polarization function basis sets. The simulation predicts an average H? Cl bond distance of 1.28 Å, which is in good agreement with the experimental value. The H_HCl···O_w and Cl_HCl···H_w distances of 1.84 and 3.51 Å were found for the first hydration shell. At the hydrogen site of HCl, a single water molecule is the most preferred coordination, whereas an average coordination number of 12 water molecules of the full first shell was observed for the chloride site. The hydrogen bonding at the hydrogen site of HCl is weakened by proton transfer reactions and an associated lability of ligand binding. Two proton transfer processes were observed in the QMCF MD simulation, demonstrating acid dissociation of HCl. A weak structure‐making/breaking effect of HCl in water is recognized from the mean residence times of 2.1 and 0.8 ps for ligands in the neighborhood of Cl and H sites of HCl, respectively. © 2009 Wiley Periodicals, Inc. J Comput Chem, 2010     

Elucidating ultrafast nonradiative decay of photoexcited uracil in aqueous solution by ab initio molecular dynamics

Nonradiative decay of the photoexcited RNA base uracil has been studied in fully explicit aqueous solution using nonadiabatic ab initio molecular dynamics. Detailed comparison of the time-dependent nonadiabatic transition probability with specific molecular vibrational motions provides insight into the mechanism of the ultrafast internal conversion. From a monoexponential fit to the excited state ensemble population, the lifetime of the first electronically excited ππ^* singlet state has been determined to be 359 fs. Additional, reference, nonadiabatic simulations have been carried out in the gas phase, pinpointing the effects of the solvent on the photophysics of uracil. The gas phase excited state lifetime is calculated to be 608 fs, somewhat longer than in solution. In terms of excitation energies and geometrical parameters, the differences between gas phase and aqueous solution are found to be generally small. A notable exception is the excited state out-of-plane torsional motion about the CC double bond, which appears severely damped by the solvent. Moreover, hydrogen bond interactions between the uracil oxygens and the solvent hydrogens are seen to enhance internal conversion.     

Hyperfine interactions in aqueous solution of Cr3+: an ab initio molecular dynamics study

We performed an ab initio molecular dynamics simulation of the paramagnetic transition metal ion Cr³⁺ in aqueous solution. Isotropic hyperfine coupling constants between the electron spin of the chromium ion and nuclear spins of all water molecules have been determined for instantaneous snapshots extracted from the trajectory. The coupling constant of first sphere oxygen, A _iso(¹⁷O_I)=1.9±0.3 MHz, is independent on Cr–O_I distance but increases with the tilt angle for the water molecule approaching 180°. First sphere hydrogen spins have A _iso(¹ H_I)=2.1±0.2 MHz which decreases with increasing tilt angle and shows a Cr–H_I distance dependence. The hyperfine coupling constants for second sphere ¹⁷O is negative and an order of magnitude smaller (−0.20±0.02 MHz) compared to first sphere.     

Time-reversible ab initio molecular dynamics

Time-reversible ab initio molecular dynamics based on a lossless multichannel decomposition for the integration of the electronic degrees of freedom [Phys. Rev. Lett. 97, 123001 (2006)] is explored. The authors present a lossless time-reversible density matrix molecular dynamics scheme. This approach often allows for stable Hartree-Fock simulations using only one single self-consistent field cycle per time step. They also present a generalization, introducing an additional "forcing" term, that in a special case includes a hybrid Lagrangian, i.e., Car-Parrinello-type, method, which can systematically be constrained to the Born-Oppenheimer potential energy surface by using an increasing number of self-consistency cycles in the nuclear force calculations. Furthermore, in analog to the reversible and symplectic leapfrog or velocity Verlet schemes, where not only the position but also the velocity is propagated, the authors propose a Verlet-type density velocity formalism for time-reversible Born-Oppenheimer molecular dynamics.     

Density functional theory for efficient ab initio molecular dynamics simulations in solution

We present a density functional for first-principles molecular dynamics simulations that includes the electrostatic effects of a continuous dielectric medium. It allows for numerical simulations of molecules in solution in a model polar solvent. We propose a smooth dielectric model function to model solvation into water and demonstrate its good numerical properties for total energy calculations and constant energy molecular dynamics.     

Controlling the shape and flexibility of arylamides: a combined ab initio, ab initio molecular dynamics, and classical molecular dynamics study

Using quantum chemistry plus ab initio molecular dynamics and classical molecular dynamics methods, we address the relationship between molecular conformation and the biomedical function of arylamide polymers. Specifically, we have developed new torsional parameters for a class of these polymers and applied them in a study of the interaction between a representative arylamide and one of its biomedical targets, the anticoagulant drug heparin. Our main finding is that the torsional barrier of a C(aromatic)-C(carbonyl) bond increases significantly upon addition of an o-OCH2CH2NH3+ substituent on the benzene ring. Our molecular dynamics studies that are based on the original general AMBER force field (GAFF) and GAFF modified to include our newly developed torsional parameters show that the binding mechanism between the arylamide and heparin is very sensitive to the choice of torsional potentials. Ab initio molecular dynamics simulation of the arylamide independently confirms the degree of flexibility we obtain by classical molecular dynamics when newly developed torsional potentials are used.     

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.     

The Jahn-Teller effect of the TiIII ion in aqueous solution: extended ab initio QM/MM molecular dynamics simulations.

Combined ab initio quantum mechanical/molecular mechanical (QM/MM) molecular dynamics simulations, including only the first and the first and second hydration shells in the QM region, were performed for TiIII in aqueous solution. The hydration structure of TiIII is discussed in terms of radial distribution functions, coordination-number distributions and several angle distributions. Dynamical properties, such as librational and vibrational motions and TiIII-O vibrations, were evaluated. A fast dynamical Jahn-Teller effect of TiIII(aq) was observed in the QM/MM simulations, in particular when the second hydration shell was included into the QM region. The results justify the computational effort required for the inclusion of the second hydration shell into the QM region and show the importance of this effort for obtaining accurate hydration-shell geometries, dynamical properties, and details of the Jahn-Teller effect.     

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.     

The solvation structure of Pb(II) in dilute aqueous solution: an ab initio quantum mechanical/molecular mechanical molecular dynamics approach

Structural properties of the hydrated Pb(II) ion have been investigated by ab initio quantum mechanical/molecular mechanical molecular dynamics simulations at Hartree-Fock quantum mechanical level. The first shell coordination number was found to be nine, and several other structural parameters such as angular distribution functions, radial distribution functions, and tilt- and theta-angle distributions allow the full characterization of the hydration structure of the Pb(II) ion.     

Quantum wavepacket ab initio molecular dynamics for extended systems

In this paper, we introduce a symmetry-adapted quantum nuclear propagation technique that utilizes distributed approximating functionals for quantum wavepacket dynamics in extended condensed-phase systems. The approach is developed with a goal for implementation in quantum-classical methods such as the recently developed quantum wavepacket ab intio molecular dynamics (QWAIMD) to facilitate the study of extended systems. The method has been numerically benchmarked for extended electronic systems as well as protonic conducting systems that benefit from quantum nuclear treatment. Vibrational properties are computed for the case of the protonic systems through use of a novel velocity-flux correlation function. The treatment is found to be numerically accurate and efficient.     

Structure and dynamics of Au+ ion in aqueous solution: ab initio QM/MM MD simulations

Structure and dynamics of hydrated Au(+) have been investigated by means of molecular dynamics simulations based on ab initio quantum mechanical molecular mechanical forces at Hartree-Fock level for the treatment of the first hydration shell. The outer region of the system was described using a newly constructed classical three-body corrected potential. The structure was evaluated in terms of radial and angular distribution functions and coordination number distributions. Water exchange processes between coordination shells and bulk indicate a very labile structure of the first hydration shell whose average coordination number of 4.7 is a mixture of 3-, 4-, 5-, 6-, and 7-coordinated species. Fast water exchange reactions between first and second hydration shell occur, and the second hydration shell is exceptionally large. Therefore, the mean residence time of water molecules in the first hydration shell (5.6 ps/7.5 ps for t*= 0.5 ps/2.0 ps) is shorter than that in the second shell (9.4 ps/21.2 ps for t*= 0.5 ps/2.0 ps), leading to a quite specific picture of a "structure-breaking" effect.     

Accelerating chemical reactions: exploring reactive free-energy surfaces using accelerated ab initio molecular dynamics

A biased potential molecular dynamics simulation approach, accelerated molecular dynamics (AMD), has been implemented in the framework of ab initio molecular dynamics for the study of chemical reactions. Using two examples, the double proton transfer reaction in formic acid dimer and the hypothetical adiabatic ring opening and subsequent rearrangement reactions in methylenecyclopropane, it is demonstrated that ab initio AMD can be readily employed to efficiently explore the reactive potential energy surface, allowing the prediction of chemical reactions and the identification of metastable states. An adaptive variant of the AMD method is developed, which additionally affords an accurate representation of both the free-energy surface and the mechanism associated with the chemical reaction of interest and can also provide an estimate of the reaction rate.     

Like-charge guanidinium pairing from molecular dynamics and ab initio calculations

Pairing of guanidinium moieties in water is explored by molecular dynamics simulations of short arginine-rich peptides and ab initio calculations of a pair of guanidinium ions in water clusters of increasing size. Molecular dynamics simulations show that, in an aqueous environment, the diarginine guanidinium like-charged ion pairing is sterically hindered, whereas in the Arg-Ala-Arg tripeptide, this pairing is significant. This result is supported by the survey of protein structure databases, where it is found that stacked arginine pairs in dipeptide fragments exist solely as being imposed by the protein structure. In contrast, when two arginines are separated by a single amino acid, their guanidinium groups can freely approach each other and they frequently form stacked pairs. Molecular dynamics simulations results are also supported by ab initio calculations, which show stabilization of stacked guanidinium pairs in sufficiently large water clusters.     

Photodissociation of diiodomethane in acetonitrile solution and fragment recombination into iso-diiodomethane studied with ab initio molecular dynamics simulations

Series of hydrates of the most stable glycine-H+/2H2+ in the gas phase are presented at the B3LYP level. The results show that only the amino hydrogens and hydroxyl hydrogens can be monohydrated for the glycine-H+, and the amino hydrogens are preferred. The H6(O4) of glycine-2H2+ is the best site for a water molecule to attach, i.e., the corresponding hydrate is the most stable one among its isomers. Calculations reveal that the binding energies of hydrated hydrogens decrease relative to their counterparts in the isolated glycine-H+/2H2+ complexes and they are positive values and without proton transfer except those of monohydrated glycine-2H2+ complexes with the combination modes of H3O+...(glycine-H+). The complex H3O+...(glycine-H+) is formed by the combination of a H2O molecule and one hydroxyl-site proton of glycine-2H2+, and with the proton transfer to H2O. Here the interaction between the proton of H3O+ and the glycine-H+ mainly depends on an electronic one instead of an initial covalent one of the isolated glycine-2H2+. The generation of the bond between the H3O+ and the glycine-H+ makes the energy of the complex higher than the energy sum of its two separated species (or two reactants of the complex), just like the case of M+...(glycine-H+) bond (M = Li,Na). The observation can explain satisfactorily why the combinations of both a proton and an alkali ion or two alkali ions to a glycine molecule can make the corresponding complex hold reservation energy bond(s), while the combination of two protons and a glycine in our previous work cannot [H. Ai et al., J. Chem. Phys. 117, 7593 (2002)]. For the glycine-2H2+, monohydration at the any site of its amino hydrogens can make the binding strength of any other neighboring proton (hydrogens) stronger relative to its counterpart in the isolated glycine-2H2+. Further hydration, especially at the site of either of hydroxyl hydrogens, would disfavor the reservation energy of the system.     

Electron hydration dynamics in water clusters: A direct ab initio molecular dynamics approach

Electron attachment dynamics of excess electron in water cluster (H2O)n (n = 2 and 3) have been investigated by means of full-dimensional direct ab initio molecular dynamics (MD) method at the MP26-311++G(d,p) level. It was found that the hydrogen bond breaking due to the excess electron is an important process in the first stage of electron capture in water trimer. Time scale of electron localization and hydrogen bond breaking were determined by the direct ab initio MD simulation. The initial process of hydration in water cluster is clearly visualized in the present study. In n = 3, an excess electron is first trapped around the cyclic water trimer with a triangular form, where the excess electron is equivalently distributed on the three water molecules at time zero. After 50 fs, the excess electron is concentrated into two water molecules, while the potential energy of the system decreases by -1.5 kcal/mol from the vertical point. After 100 fs, the excess electron is localized in one of the water molecules and the potential energy decreases by -5.3 kcal/mol, but the triangular form still remained. After that, one of the hydrogen bonds in the triangular form is gradually broken by the excess electron, while the structure becomes linear at 100-300 fs after electron capture. The time scale of hydrogen bond breaking due to the excess electron is calculated to be about 300 fs. Finally, a dipole bound state is formed by the linear form of three water molecules. In the case of n = 2, the dipole bound anion is formed directly. The mechanism of electron hydration dynamics was discussed on the basis of theoretical results.