Lattice dynamics of tetrahedrally bonded boron nitride probed by inelastic X-ray scattering

We report experimental results on the lattice dynamics of zincblende and wurtzite boron nitride obtained by inelastic X-ray scattering (IXS) on polycrystalline samples. The generalized vibrational density of states and the orientation averaged longitudinal acoustic phonon dispersion are determined, and the longitudinal and shear sound speeds as well as the Debye temperature are derived. Our results are compared to ab initio lattice dynamics calculations and available elastic and thermodynamic data.     

Lattice dynamics of solid α-carbon monoxide

The lattice dynamics of α-CO is studied assuming a fully-ordered antiferro P2₁3 structure, using a potential model which includes electrostatic, dispersion and exchange interactions. The coupling of translational and librational motions leads to mixed modes at finite wave vectors and to anomalous acoustic dispersion. Anharmonic frequency shifts and the damping due to phonon—phonon interactions are evaluated for the zero wave vector modes. The single-particle potential wells for the high-frequency librational modes are investigated and suggest that translation—rotation coupling plays a primary role in the molecular reorientations.     

Ab initio lattice dynamics of nonconducting crystals by systematic fragmentation

A systematic method for approximating the ab initio electronic energy of crystal lattices has been improved by the incorporation of long range electrostatic and dispersion interactions. The effect of these long range interactions on the optimization of the crystal structure is reported. The harmonic lattice dynamics have been evaluated to give phonon frequencies and neutron scattering intensities. Exemplary results are reported for diamond, silicon, and α-quartz using Hartree-Fock, M?ller-Plesset perturbation, and coupled-cluster levels of ab initio theory.     

Quantum wave packet ab initio molecular dynamics: an approach to study quantum dynamics in large systems

A methodology to efficiently conduct simultaneous dynamics of electrons and nuclei is presented. The approach involves quantum wave packet dynamics using an accurate banded, sparse and Toeplitz representation for the discrete free propagator, in conjunction with ab initio molecular dynamics treatment of the electronic and classical nuclear degree of freedom. The latter may be achieved either by using atom-centered density-matrix propagation or by using Born-Oppenheimer dynamics. The two components of the methodology, namely, quantum dynamics and ab initio molecular dynamics, are harnessed together using a time-dependent self-consistent field-like coupling procedure. The quantum wave packet dynamics is made computationally robust by using adaptive grids to achieve optimized sampling. One notable feature of the approach is that important quantum dynamical effects including zero-point effects, tunneling, as well as over-barrier reflections are treated accurately. The electronic degrees of freedom are simultaneously handled at accurate levels of density functional theory, including hybrid or gradient corrected approximations. Benchmark calculations are provided for proton transfer systems and the dynamics results are compared with exact calculations to determine the accuracy of the approach.     

Ab initio centroid path integral molecular dynamics: application to vibrational dynamics of diatomic molecular systems

An ab initio centroid molecular dynamics (CMD) method is developed by combining the CMD method with the ab initio molecular orbital method. The ab initio CMD method is applied to vibrational dynamics of diatomic molecules, H2 and HF. For the H2 molecule, the temperature dependence of the peak frequency of the vibrational spectral density is investigated. The results are compared with those obtained by the ab initio classical molecular dynamics method and exact quantum mechanical treatment. It is shown that the vibrational frequency obtained from the ab initio CMD approaches the exact first excitation frequency as the temperature lowers. For the HF molecule, the position autocorrelation function is also analyzed in detail. The present CMD method is shown to well reproduce the exact quantum result for the information on the vibrational properties of the system.     

Extracting effective normal modes from equilibrium dynamics at finite temperature

A general method for obtaining effective normal modes of a molecular system from molecular dynamics simulations is presented. The method is based on a localization criterion for the Fourier transformed velocity time-correlation functions of the effective modes. For a given choice of the localization function used, the method becomes equivalent to the principal mode analysis (PMA) based on covariance matrix diagonalization. On the other hand, a proper choice of the localization function leads to a novel method with a strong analogy with the usual normal mode analysis of equilibrium structures, where the Hessian system at the minimum energy structure is replaced by the thermal averaged Hessian, although the Hessian itself is never actually calculated. This method does not introduce any extra numerical cost during the simulation and bears the same simplicity as PMA itself. It can thus be readily applied to ab initio molecular dynamics simulations. Three such examples are provided here. First we recover effective normal modes of an isolated formaldehyde molecule computed at 20 K in very good agreement with the results of a normal mode analysis performed at its equilibrium structure. We then illustrate the applicability of the method for liquid phase studies. The effective normal modes of a water molecule in liquid water and of a uracil molecule in aqueous solution can be extracted from ab initio molecular dynamics simulations of these two systems at 300 K.     

Multiple time scale Hartree–Fock molecular dynamics

This is the first application of a rigorous, established multiple time-step method to ab initio molecular dynamics. The resulting algorithm is conceptually simple and easy to implement, but very effective. It translates the large mass differences present in ab initio molecular dynamics into substantial savings in computer time while retaining high accuracy. This transforms ab initio molecular dynamics from a desirable but prohibitively expensive possibility into a viable method, at least for short-time phenomena in small systems or for otherwise inaccessibly complicated potential energy surfaces. © 1993 John Wiley & Sons, Inc.     

Nonequilibrium solid phase formation studied by lattice dynamics calculation and ion beam mixing in an immiscible Co-Ag system

For the equilibrium immiscible Co-Ag system, a proven realistic ab initio derived n-body potential is applied to study the nonequilibrium solid phase formation at three chemical stoichiometries of Co/Ag = 1:3, 1:1, and 3:1. To predict the structural stability, the elastic constants and the phonon spectra are calculated at the chosen stoichiometries with a total of eight hypothetical crystalline structures. The calculated results suggest that four compounds, that is, D0(3) CoAg3, B1 CoAg, B2 CoAg, and D0(3) Co3Ag, are unstable, as they all feature negative elastic constants as well as imaginary phonons, and that another four compounds of both fcc-type L1(2) and hcp-type D0(19) structures at chemical stoichiometries of Co/Ag = 1:3 and 3:1, respectively, may elastically be favored and therefore obtainable under some specific conditions. It is also found that all the calculated elastic constants and phonon spectra are coincident within the framework of the elastic theory. Moreover, the calculated elastic constants are in good agreement with those acquired directly from ab initio calculations, lending support to the validity of the ab initio derived n-body Co-Ag potential as well as its resultant elastic constants and the phonon spectra. Interestingly, some of the predicted nonequilibrium solid phases, that is, two hcp-type compounds at chemical stoichiometries of Co/Ag = 1:3 and 3:1, respectively, are indeed obtained in ion beam mixing experiments and their lattice constants determined by diffraction analysis are in good agreement with those from calculations.     

Ab initio molecular dynamics study of formate ion hydration

We perform ab initio molecular dynamics simulations of the aqueous formate ion. The mean number of water molecules in the first solvation shell, or the hydration number, of each formate oxygen is found to be consistent with recent experiments. Our ab initio pair correlation functions, however, differ significantly from many classical force field results and hybrid quantum mechanics/molecular mechanics predictions. They yield roughly one less hydrogen bond between each formate oxygen and water than force field or hybrid methods predict. Both the BLYP and PW91 exchange correlation functionals give qualitatively similar results. The time dependence of the hydration numbers are examined, and Wannier function techniques are used to analyze electronic configurations along the molecular dynamics trajectory.     

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.     

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.     

Modeling the sorption dynamics of NaH using a reactive force field

We have parametrized a reactive force field for NaH, ReaxFF(NaH), against a training set of ab initio derived data. To ascertain that ReaxFF(NaH) is properly parametrized, a comparison between ab initio heats of formation of small representative NaH clusters with ReaxFF(NaH) was done. The results and trend of ReaxFF(NaH) are found to be consistent with ab initio values. Further validation includes comparing the equations of state of condensed phases of Na and NaH as calculated from ab initio and ReaxFF(NaH). There is a good match between the two results, showing that ReaxFF(NaH) is correctly parametrized by the ab initio training set. ReaxFF(NaH) has been used to study the dynamics of hydrogen desorption in NaH particles. We find that ReaxFF(NaH) properly describes the surface molecular hydrogen charge transfer during the abstraction process. Results on heat of desorption versus cluster size shows that there is a strong dependence on the heat of desorption on the particle size, which implies that nanostructuring enhances desorption process. To gain more insight into the structural transformations of NaH during thermal decomposition, we performed a heating run in a molecular dynamics simulation. These runs exhibit a series of drops in potential energy, associated with cluster fragmentation and desorption of molecular hydrogen. This is consistent with experimental evidence that NaH dissociates at its melting point into smaller fragments.     

Ab initio dynamics with wave-packets and density matrices

Efficient methodologies to conduct simultaneous dynamics of electrons and nuclei are discussed. Particularly, attention is directed to a recent development that combines quantum dynamics with ab initio molecular dynamics. The two components of the methodology, namely, quantum dynamics and ab initio molecular dynamics, are harnessed together using a time-dependent self-consistent field-like coupling procedure. An approach to conduct quantum dynamics using an accurate banded, sparse and Toeplitz representation for the discrete free propagator is highlighted with suitable review of other related approaches. One notable feature of the method is that all important quantum dynamical effects including zero-point effects, tunneling as well as over-barrier reflections are accurately treated. Computational methodologies for improved efficiency of the quantum dynamics are also discussed. There exists a number of ways to carry out simultaneous ab initio molecular dynamics (such as Born–Oppenheimer dynamics and extended Lagrangian dynamics, Car–Parrinello dynamics being a prime example of the latter); our prime focus remains on atom-centered density-matrix propagation and Born–Oppenheimer dynamics. The electronic degrees of freedom are handled at accurate levels of density functional theory, using hybrid or gradient corrected approximations. Benchmark calculations are provided for a prototypical proton transfer system. Future generalizations and goals are discussed.     

Towards accurate solvation dynamics of divalent cations in water using the polarizable amoeba force field: From energetics to structure

Molecular dynamics simulations were performed using a modified amoeba force field to determine hydration and dynamical properties of the divalent cations Ca2+ and Mg2+. The extension of amoeba to divalent cations required the introduction of a cation specific parametrization. To accomplish this, the Thole polarization damping model parametrization was modified based on the ab initio polarization energy computed by a constrained space orbital variation energy decomposition scheme. Excellent agreement has been found with condensed phase experimental results using parameters derived from gas phase ab initio calculations. Additionally, we have observed that the coordination of the calcium cation is influenced by the size of the periodic water box, a recurrent issue in first principles molecular dynamics studies.     

First principles calculation of the thermodynamic properties of silicon clusters

The thermodynamic properties of Si clusters are calculated using first principles quantum methods combined with molecular dynamics for simulating the trajectories of clusters. A plane wave basis is used with ab initio pseudo potentials and the local density approximation for determining the electronic energies and forces. Langevin molecular dynamics simulates thermal contact with a constant temperature reservoir. Vibrational spectra, moments of inertia, anharmonic corrections, and free energies are predicted for Si₂ through Si₅. The translational contribution is based on the ideal gas limit. The rotation contribution is approximated using a classical rigid rotator. Vibrational modes are determined from the dynamical matrix in the harmonic approximation. Corrections due to anharmonicity and coupling between rotational and vibrational modes are fit from the molecular dynamics simulations. Received: 17 September 1997 / Accepted: 14 October 1997     

"Pump-probe" atom-centered density matrix propagation studies to gauge anharmonicity and energy repartitioning in atmospheric reactive adducts: case study of the OH + isoprene and OH + butadiene reaction intermediates

Time-resolved "pump-probe" ab initio molecular dynamics studies are constructed to probe the stability of reaction intermediates, the mechanism of energy transfer, and energy repartitioning, for moieties involved during the interaction of volatile organic compunds with hydroxyl radical. These systems are of prime importance in the atmosphere. Specifically, the stability of reaction intermediates of hydroxyl radical adducts to isoprene and butadiene molecules is used as a case study to develop novel computational techniques involving "pump-probe" ab initio molecular dynamics. Starting with the various possible hydroxyl radical adducts to isoprene and butadiene, select vibrational modes of each of the adducts are populated with excess energy to mimic the initial conditions of an experiment. The flow of energy into the remaining modes is then probed by subjecting the excited adducts to ab initio molecular dynamics simulations. It is found that the stability of the adducts arises directly due to the anhormonically driven coupling of the modes to facilitate repartitioning of the excess vibrational energy. This kind of vibrational repartitioning has a critical influence on the energy density.     

Improved hybrid functional for solids: the HSEsol functional

We introduce the hybrid functional HSEsol. It is based on PBEsol, a revised Perdew-Burke-Ernzerhof functional, designed to yield accurate equilibrium properties for solids and their surfaces. We present lattice constants, bulk moduli, atomization energies, heats of formation, and band gaps for extended systems, as well as atomization energies for the molecular G2-1 test set. Compared to HSE, significant improvements are found for lattice constants and atomization energies of solids, but atomization energies of molecules are slightly worse than for HSE. Additionally, we present zero-point anharmonic expansion corrections to the lattice constants and bulk moduli, evaluated from ab initio phonon calculations.     

Hybrid approach for ab initio molecular dynamics simulation combining energy density analysis and short-time Fourier transform: energy transfer spectrogram

We propose a new analysis technique for specifying molecular vibrational modes related with intramolecular and/or intermolecular energy transfer in ab initio molecular dynamics simulation of chemical reaction. The technique combines the short-time Fourier transform method with energy density analysis, which partitions the quantum chemical potential energy in the system into atomic contributions. The image obtained by the combined scheme, termed an energy transfer spectrogram (ETS), enables us to understand the dynamics of energy transfer by time-frequency representation. The time change of the local energy is quite important in chemical reactions. In order to assess the performance of the ETS, its application to the collision reaction between two carbon dioxide molecules is shown.     

The influence of temperature and density functional models in ab initio molecular dynamics simulation of liquid water

The performance of density functional theory methods for the modeling of condensed aqueous systems is hard to predict and validation by ab initio molecular simulation of liquid water is absolutely necessary. In order to assess the reliability of these tests, the effect of temperature on the structure and dynamics of liquid water has been characterized with 16 simulations of 20 ps in the temperature range of 280-380 K. We find a pronounced influence of temperature on the pair correlation functions and on the diffusion constant including nonergodic behavior on the time scale of the simulation in the lower temperature range (which includes ambient temperature). These observations were taken into account in a consistent comparison of a series of density functionals (BLYP, PBE, TPSS, OLYP, HCTH120, HCTH407). All simulations were carried out using an ab initio molecular dynamics approach in which wave functions are represented using Gaussians and the density is expanded in an auxiliary basis of plane waves. Whereas the first three functionals show similar behavior, it is found that the latter three functionals yield more diffusive dynamics and less structure.