Vibrational frequencies in Car-Parrinello molecular dynamics

Second-Order Mass-Weighting Scheme for Atom-Centered Density Matrix Propagation Molecular Dynamics

The atom-centered density matrix propagation (ADMP) method is an extended Lagrangian approach to ab initio molecular dynamics, which includes the density matrix in an orthonormalized atom-centered Gaussian basis as additional, fictitious, electronic degrees of freedom, classically propagated along with the nuclear ones. A high adiabaticity between the nuclear and electronic subsystems is mandatory in order to keep the trajectory close to the Born-Oppenheimer (BO) surface. In this regard, the fictitious electronic mass μ, being a symmetric, nondiagonal matrix in its most general form, represents a free parameter, exploitable to optimize the propagation of the electronic density. Although mass-weighting schemes in ADMP exist, a systematic procedure to define an optimal value of the fictitious masses is not available yet. In this work, in order to rationally evaluate the electronic mass, fictitious electronic normal modes are defined through the diagonalization of the Hessian of the electronic density matrix. If the same frequency is imposed on all such modes (compatible with the chosen integration time step), then the corresponding μ matrix can be calculated and then employed for the following propagation. Analysis of several ADMP test simulations reveals that such Hessian-based mass-weighting approach is able to ensure, together with a 0.1/0.2 fs time steps, a high separation between the (real) nuclear and the (fictitious) electronic frequencies, which determines a high adiabaticity. This high, unprecedented, accuracy in the propagation leads, in turn, to low errors in the estimated nuclear vibrational frequencies, making the ADMP method totally comparable to a fully converged BO molecular dynamics simulation but more computationally efficient. This work, therefore, contributes to a further development of the ADMP ab initio molecular dynamics method, aimed at improving its accuracy through a more rational evaluation of the fictitious electronic mass parameter.

Born-Oppenheimer Molecular Dynamics with a Linear Combination of Atomic Orbitals and Hybrid Functionals for Condensed Matter Simulations Made Possible. Theory and Performance for the Microcanonical and Canonical Ensembles

The implementation of an original Born-Oppenheimer molecular dynamics module is presented, which is able to perform simulations of large and complex condensed phase systems for sufficiently long time scales at the level of density functional theory with hybrid functionals, in the microcanonical (NVE) and canonical (NVT) ensembles. The algorithm is fully integrated in the Crystal code, a program for quantum mechanical simulations of materials, whose peculiarity stems from the use of atom-centered basis functions within a linear combination of atomic orbitals to describe the wave function. The corresponding efficiency in the evaluation of the exact Fock exchange series has led to the implementation of a rich variety of hybrid density functionals at a low computational cost. In addition, the molecular dynamics implementation benefits also from the effective MPI parallelization of the code, suited to exploit high-performance computing resources available on current generation supercomputer architectures. Furthermore, the information contained in the trajectory of the dynamics is extracted through a series of postprocessing algorithms that provide the radial distribution function, the diffusion coefficient and the vibrational density of states. In this work, we present a detailed description of the theoretical framework and the algorithmic implementation, followed by a critical evaluation of the accuracy and parallel performance (e.g., strong and weak scaling) of this approach, when ice and liquid water simulations are performed in the microcanonical and canonical ensembles.

A molecular dynamics study of guest-host hydrogen bonding in alcohol clathrate hydrates

Clathrate hydrates are typically stabilized by suitably sized hydrophobic guest molecules. However, it has been experimentally reported that isomers of amyl-alcohol C5H11OH can be enclosed into the 5(12)6(4) cages in structure II (sII) clathrate hydrates, even though the effective radii of the molecules are larger than the van der Waals radii of the cages. To reveal the mechanism of the anomalous enclathration of hydrophilic molecules, we performed ab initio and classical molecular dynamics simulations (MD) and analyzed the structure and dynamics of a guest-host hydrogen bond for sII 3-methyl-1-butanol and structure H (sH) 2-methyl-2-butanol clathrate hydrates. The simulations clearly showed the formation of guest-host hydrogen bonds and the incorporation of the O-H group of 3-methyl-1-butanol guest molecules into the framework of the sII 5(12)6(4) cages, with the remaining hydrophobic part of the amyl-alcohol molecule well accommodated into the cages. The calculated vibrational spectra of alcohol O-H bonds showed large frequency shifts due to the strong guest-host hydrogen bonding. The 2-methyl-2-butanol guests form strong hydrogen bonds with the cage water molecules in the sH clathrate, but are not incorporated into the water framework. By comparing the structures of the alcohols in the hydrate phases, the effect of the location of O-H groups in the butyl chain of the guest molecules on the crystalline structure of the clathrate hydrates is indicated.

Asymmetric environmental effects on the structure and vibrations of cis-[Pt(NH3)2Cl2] in condensed phases

We report the structural and vibrational properties of anticancer drug cisplatin (cis-[Pt(NH3)2Cl2]) in gas phase, in solid phase, and in aqueous solution using density functional theory (DFT) calculations, quantum mechanical/molecular mechanical (QM/MM) molecular dynamics, and effective normal modes analysis. In contrast with the gas-phase case, asymmetric hydrogen bonding environments are found in both solid phase and aqueous solution. It is shown that the discrepancy of the molecular geometry between previous gas phase calculations and the X-ray crystal structure can be resolved by considering intermolecular hydrogen bonds in the calculations of solid phase. In addition, our simulations in solid phase and aqueous solution reveal that asymmetric environmental effects lead to several spectral features observed in experiments, such as the blue-shift in the N-H stretching region and the frequency splitting of NH3 symmetric deformation modes. Furthermore, a similar decoupling and localization of several vibrational modes of cisplatin is found in solid phase and aqueous solution, in comparison to those of O-H stretching modes of water molecules in liquid water [ J. Phys. Chem. Lett. 2013 , 4 ( 19 ), 3245 - 3250 ].

Systematic improvement of a classical molecular model of water

We report the iAMOEBA ("inexpensive AMOEBA") classical polarizable water model. The iAMOEBA model uses a direct approximation to describe electronic polarizability, in which the induced dipoles are determined directly from the permanent multipole electric fields and do not interact with one another. The direct approximation reduces the computational cost relative to a fully self-consistent polarizable model such as AMOEBA. The model is parameterized using ForceBalance, a systematic optimization method that simultaneously utilizes training data from experimental measurements and high-level ab initio calculations. We show that iAMOEBA is a highly accurate model for water in the solid, liquid, and gas phases, with the ability to fully capture the effects of electronic polarization and predict a comprehensive set of water properties beyond the training data set including the phase diagram. The increased accuracy of iAMOEBA over the fully polarizable AMOEBA model demonstrates ForceBalance as a method that allows the researcher to systematically improve empirical models by efficiently utilizing the available data.

Lagrangian formulation with dissipation of Born-Oppenheimer molecular dynamics using the density-functional tight-binding method

An important element determining the time requirements of Born-Oppenheimer molecular dynamics (BOMD) is the convergence rate of the self-consistent solution of Roothaan equations (SCF). We show here that improved convergence and dynamics stability can be achieved by use of a Lagrangian formalism of BOMD with dissipation (DXL-BOMD). In the DXL-BOMD algorithm, an auxiliary electronic variable (e.g., the electron density or Fock matrix) is propagated and a dissipative force is added in the propagation to maintain the stability of the dynamics. Implementation of the approach in the self-consistent charge density functional tight-binding method makes possible simulations that are several hundred picoseconds in lengths, in contrast to earlier DFT-based BOMD calculations, which have been limited to tens of picoseconds or less. The increase in the simulation time results in a more meaningful evaluation of the DXL-BOMD method. A comparison is made of the number of iterations (and time) required for convergence of the SCF with DXL-BOMD and a standard method (starting with a zero charge guess for all atoms at each step), which gives accurate propagation with reasonable SCF convergence criteria. From tests using NVE simulations of C(2)F(4) and 20 neutral amino acid molecules in the gas phase, it is found that DXL-BOMD can improve SCF convergence by up to a factor of two over the standard method. Corresponding results are obtained in simulations of 32 water molecules in a periodic box. Linear response theory is used to analyze the relationship between the energy drift and the correlation of geometry propagation errors.