Ion-orbital coupling in Car-Parrinello calculations of hydrogen-bond vibrational dynamics: case study with the NH3-HCl dimer

We have performed Car-Parrinello molecular dynamics (CPMD) calculations of the hydrogen-bonded NH(3)-HCl dimer. Our main aim is to establish how ionic-orbital coupling in CPMD affects the vibrational dynamics in hydrogen-bonded systems by characterizing the dependence of the calculated vibrational frequencies upon the orbital mass in the adiabatic limit of Car-Parrinello calculations. We use the example of the NH(3)-HCl dimer because of interest in its vibrational spectrum, in particular the magnitude of the frequency shift of the H-Cl stretch due to the anharmonic interactions when the hydrogen bond is formed. We find that an orbital mass of about 100 a.u. or smaller is required in order for the ion-orbital coupling to be linear in orbital mass, and the results for which can be accurately extrapolated to the adiabatic limit of zero orbital mass. We argue that this is general for hydrogen-bonded systems, suggesting that typical orbital mass values used in CPMD are too high to accurately describe vibrational dynamics in hydrogen-bonded systems. Our results also show that the usual application of a scaling factor to the CPMD frequencies to correct for the effects of orbital mass is not valid. For the dynamics of the dimer, we find that the H-Cl stretch and the N-H-Cl bend are significantly coupled, suggesting that it is important to include the latter degree of freedom in quantum dynamical calculations. Results from our calculations with deuterium-substitution show that both these degrees of freedom have significant anharmonic interactions. Our calculated frequency for the H-Cl stretch using the Becke-exchange Lee-Yang-Parr correlation functional compares reasonably well with a previous second-order M?ller-Plesset calculation with anharmonic corrections, although it is low compared to the experimental value for the dimer trapped in a neon-matrix.     

On the theory underlying the Car-Parrinello method and the role of the fictitious mass parameter

The theory underlying the Car-Parrinello extended-Lagrangian approach to ab initio molecular dynamics (CPMD) is reviewed and reexamined using "heavy" ice as a test system. It is emphasized that the adiabatic decoupling in CPMD is not a decoupling of electronic orbitals from the ions but only a decoupling of a subset of the orbital vibrational modes from the rest of the necessarily coupled system of orbitals and ions. Recent work [J. Chem. Phys. 116, 14 (2002)] has pointed out that, due to the orbital-ion coupling that remains once adiabatic decoupling has been achieved, a large value of the fictitious mass mu can lead to systematic errors in the computed forces in CPMD. These errors are further investigated in the present work with a focus on those parts of these errors that are not corrected simply by rescaling the masses of the ions. It is suggested that any comparison of the efficiencies of Born-Oppenheimer molecular dynamics (BOMD) and CPMD should be performed at a similar level of accuracy. If accuracy is judged according to the average magnitude of the systematic errors in the computed forces, the efficiency of BOMD compares more favorably to that of CPMD than previous comparisons have suggested.     

Vibrational state dependence of beta and D asymmetry parameters: the case of the highest occupied molecular orbital photoelectron spectrum of methyl-oxirane

The beta angular asymmetry and D dichroic asymmetry parameters of the methyl-oxirane highest occupied molecular orbital (HOMO) band have been experimentally investigated with vibrational resolution using synchrotron radiation. A theoretical calculation of the Franck-Condon factors between vibrational ground state and different ionic vibrational states, in the Born-Oppenheimer harmonic approximation, has been performed in order to gain information on the vibrational states mainly involved in the HOMO photoelectron band. The general good agreement between theoretical and experimental results allows a reliable assignment of the major features. The experimental determination of beta and D shows their dependence on the different final vibrational states. This paper reports, for the first time, experimental evidence of the dependence of the dichroic D parameter on the vibrational excitation of the ion.     

Intermolecular Interactions and Spectroscopic Signatures of the Hydrogen-Bonded System—n-Octanol in Experimental and Theoretical Studies

n-Octanol is the object of experimental and theoretical study of spectroscopic signatures and intermolecular interactions. The FTIR measurements were carried out at 293 K for n-octanol and its deuterated form. Special attention was paid to the vibrational features associated with the O-H stretching and the isotope effect. Density Functional Theory (DFT) in its classical formulations was applied to develop static models describing intermolecular hydrogen bond (HB) and isotope effect in the gas phase and using solvent reaction field reproduced by Polarizable Continuum Model (PCM). The Atoms in Molecules (AIM) theory enabled electronic structure and molecular topology study. The Symmetry-Adapted Perturbation Theory (SAPT) was used for energy decomposition in the dimers of n-octanol. Finally, time-evolution methods, namely classical molecular dynamics (MD) and Car-Parrinello Molecular Dynamics (CPMD) were employed to shed light onto dynamical nature of liquid n-octanol with emphasis put on metric and vibrational features. As a reference, CPMD gas phase results were applied. Nuclear quantum effects were included using Path Integral Molecular Dynamics (PIMD) and a posteriori method by solving vibrational Schrödinger equation. The latter applied procedure allowed to study the deuterium isotope effect.     

Elucidating the vibrational spectra of hydrogen-bonded aggregates in solution: electronic structure calculations with implicit solvent and first-principles molecular dynamics simulations with explicit solvent for 1-hexanol in n-hexane

Fourier transform infrared spectroscopy is a popular method for the experimental investigation of hydrogen-bonded aggregates, but linking spectral information to microscopic information on aggregate size distribution and aggregate architecture is an arduous task. Static electronic structure calculations with an implicit solvent model, Car-Parrinello molecular dynamics (CPMD) using the Becke-Lee-Yang-Parr (BLYP) exchange and correlation energy functionals and classical molecular dynamics simulations for the all-atom version of the optimized parameters for liquid simulations (OPLS-AA) force field were carried out for an ensemble of 1-hexanol aggregates solvated in n-hexane. The initial configurations for these calculations were size-selected from a distribution of aggregates obtained from a large-scale Monte Carlo simulation. The vibrational spectra computed from the static electronic structure calculations for monomers and dimers and from the CPMD simulations for aggregates up to pentamers demonstrate the extent of the contribution of dangling or nondonating hydroxyl groups found in linear and branched aggregates to the "monomeric" peak. Furthermore, the computed spectra show that there is no simple relationship between peak shift and aggregate size nor architecture, but the effect of hydrogen-bond cooperativity is shown to differentiate polymer-like (cooperative) and dimer-like (noncooperative) hydrogen bonds in the vibrational spectrum. In contrast to the static electronic structure calculations and the CPMD simulations, the classical molecular dynamics simulations greatly underestimate the vibrational peak shift due to hydrogen bonding.     

Solvation of Co(III)-cysteinato complexes in water: a DFT-based molecular dynamics study

Structural, dynamical, and vibrational properties of complexes made of metal cobalt(III) coordinated to different amounts of cysteine molecules were investigated with DFT-based Car-Parrinello molecular dynamics (CPMD) simulations in liquid water solution. The systems are composed of Co(III):3Cys and Co(III):2Cys immersed in liquid water which are modeled by about 110 explicit water molecules, thus one of the biggest molecular systems studied with ab initio molecular simulations so far. In such a way, we were able to investigate structural and dynamical properties of a model of a typical metal binding site used by several proteins. Cobalt, mainly a toxicological agent, can replace the natural binding metal and thus modify the biochemical activity. The structure of the surrounding solvent around the metal-ligands complexes is reported in detail, as well as the metal-ligands coordination bonds, using radial distribution functions and electronic analyses with Mayer bond orders. Structures of the Cocysteine complexes are found in very good agreement with EXAFS experimental data, stressing the importance of considering the surrounding solvent in the modeling. A vibrational analysis is also conducted and compared to experiment, which strengthens the reliability of the solvent interactions with the Cocysteine complexes from our molecular dynamics simulations, as well as the dynamics of the systems. From this preliminary analysis, we could suggest a vibrational fingerprint able to distinguish Co(III):2Cys from Co(III):3Cys. Our simulations also show the importance of considering a quantum explicit solvent, as solute-to-solvent proton transfer events have been observed.     

Quantum and dynamical effects of proton donor-acceptor vibrational motion in nonadiabatic proton-coupled electron transfer reactions

This paper presents a general theoretical formulation for proton-coupled electron transfer (PCET) reactions. The solute is represented by a multistate valence bond model, and the active electrons and transferring proton(s) are treated quantum mechanically. This formulation enables the classical or quantum mechanical treatment of the proton donor-acceptor vibrational mode, as well as the dynamical treatment of the proton donor-acceptor mode and the solvent. Nonadiabatic rate expressions are presented for PCET reactions in a number of well-defined limits for both dielectric continuum and molecular representations of the environment. The dynamical rate expressions account for correlations between the fluctuations of the proton donor-acceptor distance and the nonadiabatic PCET coupling. The quantities in the rate expressions can be calculated with a dielectric continuum model or a molecular dynamics simulation of the full system. The significance of the quantum and dynamical effects of the proton donor-acceptor mode is illustrated with applications to model PCET systems.     

Ab initio molecular-dynamics study of liquid formamide

Properties of neat liquid formamide (HCONH2) have been studied by the combination of gradient-corrected density-functional theory, norm-conserving pseudopotentials, and the adaptive finite-element method. The structural and dynamical quantities have been calculated through molecular dynamics simulations under the Born-Oppenheimer approximation. Satisfactory agreement with experimental data was obtained for both intramolecular and intermolecular properties. Our results are also compared with those of the empirical potential functions to clarify their accuracies.     

A density functional theory study of vibrational coupling in the amide I band of beta-sheet models

We report the first molecular orbital/density functional theory (DFT) calculations on the vibrational frequencies involved in the amide I band of completely geometrically optimized models for beta-sheet peptides based upon (up to 16) glycine residues. These calculations use the B3LYP/D95** level of DFT. The primary means of vibrational coupling occurs through H bond, rather than through space, interactions, which is consistent with a previous report on alpha-helical polyalanines and H-bonding chains of both formamides and 4-pyridones. We decoupled the C=O stretching vibrations using selected 14C substitutions to probe the coupling mechanism and to determine "natural" frequencies for individual 14C=Os. The intermolecular H-bonding interactions affect the geometries of the amide groups. Those near the center of H-bonding chains have long C=O bonds. The C=O bond lengths correlate with these "natural" frequencies, The frequencies obtained from the DFT calculations are generally more coupled, and the most intense are more red shifted than those calculated by transition dipole coupling (TDC). TDC inverts the order of the shifted frequencies compared to DFT in several cases.     

Nonstatistical effects in the dissociation of ethyl radical: Finding order in chaos

How does one identify order in complex dynamical systems? A Born-Oppenheimer molecular dynamics simulation of the dissociation of ethyl radical, C(2)H(5), produces an ensemble of classical trajectories which are decomposed in the time-frequency domain using wavelets. A time-dependent scalar metric, the normalized instantaneous orbital complexity, is constructed and shown to correlate not only to the more conventional Lyapunov exponents but also to the dissociation time for an individual trajectory. The analysis of the ensemble of trajectories confirms that the long-lived trajectories are associated with a low degree of ergodicity. While the analysis of molecular dissociation dynamics is the narrow focus of the present work, the method is more general for discovery and identification of ordered regimes within large sets of chaotic data.     

Vibrational echo spectroscopy of aqueous sodium bromide solutions from first principles simulations

A theoretical study of the time-dependent vibrational echo spectroscopy of sodium bromide solutions in deuterated water at two different concentrations of 0.5 and 5.0 M and at temperatures of 300 and 350 K is presented using the method of ab initio molecular dynamics simulations. The instantaneous fluctuations in frequencies of local OD stretch modes are calculated using time-series analysis of the simulated trajectories. The third-order polarization and intensities of three pulse photon-echo are calculated from ab initio simulations. The timescales of vibrational spectral diffusion are determined from the frequency time correlation functions (FTCF) and short-time slope of three pulse photon echo (S3PE) calculated within the second-order cumulant and Condon approximations. It is found that under ambient conditions, the rate of vibrational spectral diffusion becomes slower with increase in ionic concentration. Decay of S3PE calculated for different systems give timescales, which are in close agreement with those of FTCF and also with the results of experimental time-dependent vibrational spectroscopic experiments. © 2019 Wiley Periodicals, Inc.     

Molecular vibrations of methane molecules in the structure I clathrate hydrate from ab initio molecular dynamics simulation

Vibrational frequencies of guest molecules in clathrate hydrates reflect the molecular environment and dynamical behavior of molecules. A detailed understanding of the mechanism for the vibrational frequency changes of the guest molecules in the clathrate hydrate cages is still incomplete. In this study, molecular vibrations of methane molecules in a structure I clathrate hydrate are calculated from ab initio molecular dynamics simulation. The vibrational spectra of methane are computed by Fourier transform of autocorrelation functions, which reveal distinct separation of each vibrational mode. Calculated symmetric and asymmetric stretching vibrational frequencies of methane molecules are lower in the large cages than in the small cages (8 and 16 cm(-1) for symmetric and asymmetric stretching, respectively). These changes are closely linked with the C-H bond length. The vibrational frequencies for the bending and rocking vibrational modes nearly overlap in each of the cages.     

Constrained dynamics and extraction of normal modes from ab initio molecular dynamics: application to ammonia

We present a methodology for extracting phonon data from ab initio Born-Oppenheimer molecular dynamics calculations of molecular crystals. Conventional ab initio phonon methods based on perturbations are difficult to apply to lattice modes because the perturbation energy is dominated by intramolecular modes. We use constrained molecular dynamics to eliminate the effect of bond bends and stretches and then show how trajectories can be used to isolate and define in particular, the eigenvalues and eigenvectors of modes irrespective of their symmetry or wave vector. This is done by k-point and frequency filtering and projection onto plane wave states. The method is applied to crystalline ammonia: the constrained molecular dynamics allows a significant speed-up without affecting structural or vibrational modes. All Gamma point lattice modes are isolated: the frequencies are in agreement with previous studies; however, the mode assignments are different.     

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.     

Car-Parrinello simulation of an O-H stretching envelope and potential of mean force of an intramolecular hydrogen bonded system: application to a Mannich base in solid state and in vacuum

Car-Parrinello molecular dynamics (CPMD) study was performed for an anharmonic system-an intramolecularly hydrogen bonded Mannich-base-type compound, 4,5-dimethyl-2(N,N-dimethylaminemethyl)phenol, to investigate the vibrational spectrum associated with the O-H stretching. Calculations were carried out for the solid state and for an isolated molecule. The classical CPMD simulation was performed and then the proton potential snapshots were extracted from the trajectory. The vibrational Schrodinger equation for the snapshots was solved numerically, and the (O-H) envelope was calculated as a superposition of the 0-->1 transitions. The potential of mean force for the proton stretching mode was calculated from the proton vibrational eigenfunctions and eigenvalues incorporating statistical sampling, nuclear quantum effects, and effects of the environment. Perspectives for application of the presented methodology in the computational support of biocatalysis are given in the study.     

Raman spectroscopic study on solvation of diphenylcyclopropenone and phenol blue in room temperature ionic liquids

We investigated the solvation of several room temperature ionic liquids by Raman spectroscopy using diphenylcyclopropenone (DPCP) and phenol blue (PB) as probe molecules. We estimated acceptor numbers (AN) of room temperature ionic liquids by an empirical equation associated with the Raman band of DPCP assigned as a C=C stretching mode involving a significant C=O stretching character. According to the dependence of AN on cation and anion species, the Lewis acidity of ionic liquids is considered to come mainly from the cation charge. The frequencies and bandwidths of the C=O and C=N stretching modes of phenol blue are found to be close to those in conventional polar solvents such as methanol and dimethyl sulfoxide. The frequencies of these vibrational modes show similar dependence upon the electronic absorption band center as is observed in conventional liquid solvents. However, peculiar behavior was found in the Raman bandwidths and the excitation wavelength dependence of the C=N stretching mode in room temperature ionic liquids. Both the bandwidth of the C=N stretching mode and the extent of the excitation wavelength dependence of the Raman shift of the C=N stretching mode tend to decrease as the absorption band center decreases, in contrast to the case of conventional solvents. This anomaly is discussed in terms of the properties of room temperature ionic liquids.