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.     

Quantum wavepacket ab initio molecular dynamics: an approach for computing dynamically averaged vibrational spectra including critical nuclear quantum effects

We have introduced a computational methodology to study vibrational spectroscopy in clusters inclusive of critical nuclear quantum effects. This approach is based on the recently developed quantum wavepacket ab initio molecular dynamics method that combines quantum wavepacket dynamics with ab initio molecular dynamics. The computational efficiency of the dynamical procedure is drastically improved (by several orders of magnitude) through the utilization of wavelet-based techniques combined with the previously introduced time-dependent deterministic sampling procedure measure to achieve stable, picosecond length, quantum-classical dynamics of electrons and nuclei in clusters. The dynamical information is employed to construct a novel cumulative flux/velocity correlation function, where the wavepacket flux from the quantized particle is combined with classical nuclear velocities to obtain the vibrational density of states. The approach is demonstrated by computing the vibrational density of states of [Cl-H-Cl]-, inclusive of critical quantum nuclear effects, and our results are in good agreement with experiment. A general hierarchical procedure is also provided, based on electronic structure harmonic frequencies, classical ab initio molecular dynamics, computation of nuclear quantum-mechanical eigenstates, and employing quantum wavepacket ab initio dynamics to understand vibrational spectroscopy in hydrogen-bonded clusters that display large degrees of anharmonicities.     

Enzymatic GTP hydrolysis: insights from an ab initio molecular dynamics study

Ab initio methods were used to shed light on fundamental aspects of the enzymatic mechanism of guanosine triphosphate hydrolysis in the Cdc42/Cdc42GAP complex. The calculations focused on the nucleophilic addition of the catalytic water molecule to the gamma-phosphate phosphorus atom. A large model system was required to correctly reproduce the electrostatic properties on the active site. The model turned out to reproduce most of the electrostatic field of the biological complex at the reactants. Our calculations established the H-bond pattern of the catalytic water (WAT), which turned out to interact with Q61 and T35, in the most stable conformation. This ruled out the possibility that the catalytic water transferred its proton directly to the gamma-phosphate. Furthermore, the calculations suggested that the electronic structure of WAT was very different from that in the bulk. Finally, this study showed that during the reaction, WAT transferred a proton to Gln61, consistent with the available X-ray data on a transition-state analogue/enzyme complex(19) and with the decrease of activity in the Q61E mutant.     

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.     

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.     

Vibrational predissociation dynamics of methane-Ar: an ab initio approach

We calculated the cross sections for vibrational predissociation of methane-Ar induced by excitation of the methane nu 3 mode with the aid of an ab initio CH4-Ar potential depending explicitly on the nu 3 and nu 1 normal coordinates of the CH4 monomer. We found that dissociation into CH4 fragments excited in the nu 1 mode, a V-->V' process with very low kinetic energy release, strongly dominates over direct dissociation into Ar and ground state CH4, and is responsible for the line broadening observed experimentally. The (observed and calculated) strong variation of the line widths for the Van der Waals levels excited in combination with the nu 3 mode (giving states of A, F and E symmetry) is related to the opening up of appropriate nu 1 dissociation channels and the occurrence of rotational resonances in the nu 1 continuum in the energy range of the quasi-bound nu 3 levels.     

Li ion diffusion mechanisms in LiFePO4: an ab initio molecular dynamics study

The mechanisms for thermal (self) diffusion of Li ions in fully lithiated LiFePO(4) have been investigated with spin polarized ab initio molecular dynamics calculations. The effect of electron correlation is taken into account with the GGA+U formalism. It was found that Li ion diffusion is not a continuous process but through a series of jumps from one site to another. A dominant process is the hopping between neighboring Li sites around the PO(4) groups, which results in a zigzag pathway along the crystallographic b-axis. This observation is in agreement with a recent neutron diffraction experiment. A second process involves the collaborative movements of the Fe ions leading to the formation of antisite defects and promotes Li diffusion across the Li ion channels. The finding of the second mechanism demonstrates the benefit of ab initio molecular dynamics simulation in sampling diffusion pathways that may not be anticipated.     

On the calculation of ab initio quantum molecular similarities for large systems: Fitting the electron density

A set of procedures for rapid calculation of quantum molecular similarities from ab initio wave functions is discussed. In all cases a density fitting is carried out to eliminate the need of calculating costly four-centered integrals. It is proved that this methodology can be applied to large systems to reproduce exact quantum molecular similarity measures at an extremely low computational cost. © 1994 by John Wiley & Sons, Inc.     

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.     

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.     

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.     

Chemistry at surfaces: from ab initio structures to quantum dynamics

Recent years have witnessed an ever growing interest in theoretically studying chemical processes at surfaces. Apart from the interest in catalysis, electrochemistry, hydrogen economy, green chemistry, atmospheric and interstellar chemistry, theoretical understanding of the molecule–surface chemical bonding and of the microscopic dynamics of adsorption and reaction of adsorbates are of fundamental importance for modeling known processes, understanding new experimental data, predicting new phenomena, controlling reaction pathways. In this work, we review the efforts we have made in the last few years in this exciting field. We first consider the energetics and the structural properties of some adsorbates on metal surfaces, as deduced by converged, first-principles, plane-wave calculations within the slab-supercell approach. These studies comprise water adsorption on Ru(0001), a subject of very intense debate in the past few years, and oxygen adsorption on aluminum, the prototypical example of metal passivation. Next, we address dynamical processes at surfaces with classical and quantum methods. Here the main interest is in hydrogen dynamics on metallic and semi-metallic surfaces, because of its importance for hydrogen storage and interstellar chemistry. Hydrogen sticking is studied with classical and quasi-classical means, with particular emphasis on the relaxation of hot–atoms following dissociative chemisorption. Hot atoms dynamics on metal surfaces is investigated in the reverse, hydrogen recombination process and compared to Eley–Rideal dynamics. Finally, Eley–Rideal, collision-induced desorption, and adsorbate-induced trapping are studied quantum mechanically on a graphite surface, and unexpected quantum effects are observed.     

Water mediated proton transfer in a mesostructured aluminosilicate framework: an ab initio molecular dynamics study

The proton transfer process mediated by water molecules adsorbed in an aluminosilicate framework has been studied using ab initio molecular dynamics simulations. This investigation has been carried out using a quasi-one-dimensional model simulating the mesoporous aluminosilicate channel structures. The effects of both the water loading and temperature of the system have been considered. At low coverage (one water molecule per acid site), the hydroxonium ion (H(3)O)(+) is found to be a transition state, in agreement with earlier studies on zeolites. At a higher water coverage (two water molecules per acid site), the (H(5)O(2))(+) species and the hydrogen bonded "neutral complex" structure are both found to be stable complexes at finite temperatures. The vibrational frequency spectrum is simulated by performing a Fourier transform of the velocity autocorrelation function (VAF), and the peak positions in the VAF are compared with IR measurements and zero-temperature calculations.     

Lithium hydroxide phase transition under high pressure: an ab initio molecular dynamics study.

The high-pressure phase transition in the deuterated lithium hydroxide crystalline state has been studied by Car-Parrinello molecular dynamics simulations, in the constant-pressure, constant-temperature ensemble. The recently developed metadynamics approach has been applied to encourage the system to transform into different phases in an affordable simulation time. A previously not completely characterized high-pressure phase has been obtained. The structural and spectroscopic properties have been studied and compared with the neutron scattering, infrared and Raman measurements. It has been found that the calculated structure differs slightly from the experimental hypothesis, and that the presence of strong hydrogen bonds is the source of the red shift and of the characteristic features of the OD-stretching bands in both IR and Raman spectra.     

Nanosize icosahedral quasicrystal in Mg90Ca10 glass: an ab initio molecular dynamics study

Rapid solidification of Mg(90)Ca(10) from its liquid state is studied by means of an ab initio molecular dynamics technique, and its local structure is investigated by various analyzing methods. The liquid and amorphous states are found to have slightly different short range order even though the perfect and defective icosahedral bonding environments are major bonding elements of both liquid and amorphous states. Perfect icosahedrons with a small frequency exist in the liquid state, more develop during the cooling process and they become the leading building units in the glass state, indicating an icosahedral short range order in Mg(90)Ca(10) glass. Also the linked icosahedrons lead to an icosahedral medium range order. Furthermore, an ordered arrangement of some icosahedrons in the hexagonal symmetry is observed in the glass model, representing a nanoscale icosahedral quasicrystalline phase in Mg(90)Ca(10) glass.     

Scalable fine-grained parallelization of plane-wave-based ab initio molecular dynamics for large supercomputers

Many systems of great importance in material science, chemistry, solid-state physics, and biophysics require forces generated from an electronic structure calculation, as opposed to an empirically derived force law to describe their properties adequately. The use of such forces as input to Newton's equations of motion forms the basis of the ab initio molecular dynamics method, which is able to treat the dynamics of chemical bond-breaking and -forming events. However, a very large number of electronic structure calculations must be performed to compute an ab initio molecular dynamics trajectory, making the efficiency as well as the accuracy of the electronic structure representation critical issues. One efficient and accurate electronic structure method is the generalized gradient approximation to the Kohn-Sham density functional theory implemented using a plane-wave basis set and atomic pseudopotentials. The marriage of the gradient-corrected density functional approach with molecular dynamics, as pioneered by Car and Parrinello (R. Car and M. Parrinello, Phys Rev Lett 1985, 55, 2471), has been demonstrated to be capable of elucidating the atomic scale structure and dynamics underlying many complex systems at finite temperature. However, despite the relative efficiency of this approach, it has not been possible to obtain parallel scaling of the technique beyond several hundred processors on moderately sized systems using standard approaches. Consequently, the time scales that can be accessed and the degree of phase space sampling are severely limited. To take advantage of next generation computer platforms with thousands of processors such as IBM's BlueGene, a novel scalable parallelization strategy for Car-Parrinello molecular dynamics is developed using the concept of processor virtualization as embodied by the Charm++ parallel programming system. Charm++ allows the diverse elements of a Car-Parrinello molecular dynamics calculation to be interleaved with low latency such that unprecedented scaling is achieved. As a benchmark, a system of 32 water molecules, a common system size employed in the study of the aqueous solvation and chemistry of small molecules, is shown to scale on more than 1500 processors, which is impossible to achieve using standard approaches. This degree of parallel scaling is expected to open new opportunities for scientific inquiry.     

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.     

Exciton migration dynamics in a dendritic molecule: quantum master equation approach using ab initio molecular orbital configuration interaction method

We investigate the exciton migration dynamics in a dendritic molecular model composed of pi-conjugation linear-leg units (acetylenes and diacetylene) and a benzene ring (branching point) using the quantum master equation approach with the ab initio molecular orbital (MO) configuration interaction (CI) method. The efficient migration of exciton from short-length linear legs (acetylenes) to long-length linear leg (diacetylene) via a benzene ring is observed. As predicted in previous studies, the exciton (electron and hole) distributions are relatively well localized in each generation segmented by the meta-branching point (meta-substituted benzene ring) though the electron and hole distributions are delocalized and are somewhat spatially different from each other within each generation. It is found that the excitons localized in the generation composed of short linear legs occupy in higher-lying exciton states, while those in the generation composed of long linear legs do in lower-lying ones. These features suggest the decoupling of pi-conjugation at the meta-branching point. On the other hand, the relaxation effect between exciton states is found to be caused by the exciton-phonon coupling, in which the existence of common configurations (electron-hole pairs) in CI wave functions between adjacent exciton states (having primary distributions on short and long linear-leg regions, respectively) is important for the relaxation between their exciton states. This feature indicates the importance of partial penetration of pi-conjugation through the meta-substituted benzene ring in excited states for such exciton migration.     

Hydrogen bonded structure and dynamics of liquid-vapor interface of water-ammonia mixture: an ab initio molecular dynamics study

We have carried out ab initio molecular dynamics simulations of a liquid-vapor interfacial system consisting of a mixture of water and ammonia molecules. We have made a detailed analysis of the structural and dynamical properties of the bulk and interfacial regions of the mixture. Among structural properties, we have looked at the inhomogeneous density profiles of water and ammonia molecules, hydrogen bond distributions, orientational profiles, and also vibrational frequency distributions of bulk and interfacial molecules. It is found that the interfacial molecules show preference for specific orientations so as to form water-ammonia hydrogen bonds at the interface with ammonia as the acceptor. The structure of the system is also investigated in terms of inter-atomic voids present in the system. Among the dynamical properties, we have calculated the diffusion, orientational relaxation, hydrogen bond dynamics, and vibrational spectral diffusion in bulk and interfacial regions. It is found that the diffusion and orientation relaxation of the interfacial molecules are faster than those of the bulk. However, the hydrogen bond lifetimes are longer at the interface which can be correlated with the time scales found from the decay of frequency time correlations.     

Evolution of intermolecular structure and dynamics in supercritical carbon dioxide with pressure: an ab initio molecular dynamics study

The effect of pressure on supercritical carbon dioxide (scCO2) has been characterized by using Car-Parrinello molecular dynamics simulations. Structural and dynamical properties along an isotherm of 318.15 K and at pressures ranging from 190 to 5000 bar have been obtained. Intermolecular pair correlation functions and three-dimensional atomic probability density map calculations indicate that the local environment of a central CO2 molecule becomes more structured with increasing pressure. The closest neighbors are predominantly oriented in a distorted T-shaped geometry while neighbors separated by larger distances are likely oriented in a slipped parallel arrangement. The structure of scCO2 at high densities has been compared with that of crystalline CO2. The probability distributions of intramolecular distances narrow down with increasing pressure. A marginal but non-negligible effect of pressure on the instantaneous intramolecular OCO angle is observed, lending credence to the idea that intermolecular interactions between CO2 molecules in an inhomogeneous near neighbor environment could contribute to the observed instantaneous molecular dipole moment. The extent of deviation from a perfect linear geometry of the carbon dioxide molecule decreases with increasing pressure. Time constants derived from reorientational time correlation functions of the molecular backbone compare well with experimental data. Within the range of thermodynamic conditions explored here, no significant changes are observed in the frequencies of intramolecular vibrational modes. However, a blue shift is observed in the low-frequency cage rattling mode with increasing pressure.