Ab initio molecular dynamics of protonated water clusters by integrated multicenter molecular-orbital method

To analyze large-scale cluster systems theoretically, we recently developed an "integrated multicenter molecular-orbital" (IMiC-MO) method. This method calculates the force of an entire system by dividing the system into small regions. We used the method to analyze the effect of cluster size and the process of hydrogen bond network (HBN) growth to form H(+)(H(2)O)(n) (n = 9, 17, and 33) clusters. Our simulations reveal that H(3)O(+) and water molecules in the first solvation shell function take an important role to grow the HBN. In addition, the number of hydrogen donors in each water molecule is strongly related to the shape of the HBN.     

Ab initio molecular dynamics approach to tunneling splitting in polyatomic molecules

An ab initio molecular dynamics approach is combined with the semiclassical tunneling method of Makri and Miller, which is applied to estimations of tunneling splitting in the umbrella inversion of ammonia and the intramolecular hydrogen transfer in malonaldehyde. In the application to malonaldehyde, effects of multidimensionality are examined by assigning quantum zero-point energies only to significant vibrational modes and changing the amount of energy given to other degrees of freedom. The calculated tunneling splitting values are in good agreement with the corresponding experimental values for both molecules.     

Ab initio molecular dynamics study of the keto-enol tautomerism of acetone in solution.

We have studied the keto-enol interconversion of acetone to understand the mechanism of tautomerism relevant to numerous organic and biochemical processes. Applying the ab initio metadynamics method, we simulated the keto-enol isomerism both in the gas phase and in the presence of water. For the gas-phase intramolecular mechanism we show that no other hydrogen-transfer reactions can compete with the simple keto-enol tautomerism. We obtain an intermolecular mechanism and remarkable participation of water when acetone is solvated by neutral water. The simulations reveal that C deprotonation is the kinetic bottleneck of the keto-enol transformation, in agreement with experimental observations. The most interesting finding is the formation of short H-bonded chains of water molecules that provide the route for proton transfer from the carbon to the oxygen atom of acetone. The mechanistic picture that emerged from the present study involves proton migration and emphasizes the importance of active solvent participation in tautomeric interconversion.     

Ab initio static and molecular dynamics study of 4-styrylpyridine.

We report an in‐depth theoretical study of 4‐styrylpyridine in its singlet S₀ ground state. The geometries and the relative stabilities of the trans and cis isomers were investigated within density functional theory (DFT) as well as within Hartree–Fock (HF), second‐order Møller–Plesset (MP2), and coupled cluster (CC) theories. The DFT calculations were performed using the B3LYP and PBE functionals, with basis sets of different qualities, and gave results that are very consistent with each other. The molecular structure is thus predicted to be planar at the energy minimum, which is associated with the trans conformation, and to become markedly twisted at the minimum of higher energy, which is associated with the cis conformation. The results of the calculations performed with the post‐HF methods approach those obtained with the DFT methods, provided that the level of treatment of the electronic correlation is high enough and that sufficiently flexible basis sets are used. Calculations carried out within DFT also allowed the determination of the geometry and the energy of the molecule at the biradicaloid transition state associated with the thermal cis?trans isomerization and at the transition states associated with the enantiomerization of the cis isomer and with the rotations of the pyridinyl and phenyl groups in the trans and cis isomers. Car–Parrinello molecular dynamics simulations were also performed at 50, 150, and 300 K using the PBE functional. The studies allowed us to evidence the highly flexible nature of the molecule in both conformations. In particular, the trans isomer was found to exist mainly in a nonplanar form at finite temperatures, while the rotation of the pyridinyl ring in the cis isomer was incidentally observed to take place within ≈1 ps during the simulation carried out at 150 K on this isomer.     

Ab initio molecular dynamics simulation of a water-hydrogen fluoride equimolar mixture.

Hydrogen fluoride and water can be mixed in any proportion. The resulting solutions have unique acidic properties. In particular, hydrogen fluoride undergoes a weak-to-strong acidity transition with increasing concentration of HF To supplement the knowledge already obtained on dilute or moderately concentrated solutions and gas-phase aggregates, an equimolar mixture is studied here by Car-Parrinello molecular dynamics. The natures of the ions and of the complexes formed in the equimolar liquid were determined. Specifically, H3O+, H5O2+, FH-OH2, and HF2- were spontaneously obtained while only hydronium and fluoride ions pre-exist in the equimolar crystal. The behaviour of the proton in the equimolar liquid was compared with mixtures of other proportions simulated previously in an attempt to relate proton dynamics to acidity. In the same way, the behaviour of HF2- was also examined. In this case, proton localization and transfer appeared to be driven by the fluctuating environment of the solvated ion.     

Modeling electrified metal/water interfaces from ab initio molecular dynamics: Structure and Helmholtz capacitance

The structure and potential distribution of electric double layers (EDLs) are of close relevance to the performance of electrode materials. In the past years, despite tremendous efforts devoted to this topic, an atomistic picture of the EDL is still lacking, let alone understanding on how the EDL structure is related to the dielectric property of interface water. In this article, we briefly review the recent progress in modeling electrified metal/water interfaces using ab initio molecular dynamics (AIMD). The ab initio methods for EDL modeling is firstly summarized, and then we discuss the structures of interface water on metal electrodes at different potential conditions. Moreover, we illustrate the potential-dependent behavior of chemisorbed water on Pt(111) surface and its relationship with the peak of the differential Helmholtz capacitance observed by experiment. At last, we give some perspective for future development in ab initio modeling of electrochemical interfaces.     

Organolithium complex as a gas sensing material for oxides from ab initio calculations and molecular dynamics simulations

Gas sensing study of C₂H₄Li complex toward oxides viz. CO, CO₂, NO, NO₂, SO, and SO₂ gas molecules has been carried out using ab initio method. Different possible configurations of gas molecule adsorption on C₂H₄Li complex are considered. The structural parameters of most stable configuration of gas molecule adsorbed complexes are thoroughly analysed. Electronic properties are studied using total density of states (DOS) plot. Charge transferred between the gas molecule and the substrate is studied using NBO charge analysis. Gas sensing of all the six gas molecules is possible at ambient conditions. Atom centred density matrix propagation (ADMP) molecular dynamics simulations confirmed that all the gas molecules remain adsorbed on C₂H₄Li complex at room temperature during the simulation. This study suggests that the C₂H₄Li complex acts as a novel gas sensing material for CO, CO₂, NO, NO₂, SO, and SO₂ gas molecules at ambient conditions, below room temperature as well as at high pressure.     

Ab initio molecular dynamics simulations study on initial decompositions of β‐HMX at high temperature coupled with high pressures

We performed ab initio molecular dynamics simulations to investigate initial decomposition mechanisms and subsequent chemical processes of β‐HMX (cyclotetramethylene tetranitramine) (octahydro‐1,3,5,7‐tetranitro‐1,3,5,7‐tetrazocine) crystals at high temperature coupled with high pressures. It was found that the initial decomposition step is the simultaneous C–H and N–NO₂ bond cleavage at 3,500 K. When the pressure (1–10 GPa) is applied, the first reaction steps are primarily the C–N and C–H bond fission at 3,500 K. The C–H bond cleavage is a triggering decomposition step of the HMX crystals at 3,500 K coupled with 16 GPa. This indicates that the C–H bonds are much easier to be broken and the hydrogen radicals are much more active. The applied pressures (1–10 GPa) accelerate the decompositions of HMX at 3,500 K. The decomposition pathways and time evolution of the main chemical species demonstrate that the temperature is the foremost factor that affects the decomposition at high pressures (1–10 GPa). However, the decomposition of HMX is dependent on both the temperature (3,500 K) and the pressure (16 GPa). This work will enrich the knowledge of the decompositions of condensed energetic materials under extreme conditions.     

Liquid properties of dimethyl ether from molecular dynamics simulations using Ab Initio force fields

We have used molecular dynamic simulations to study the structural and dynamical properties of liquid dimethyl ether (DME) with a newly constructed ab initio force field in this article. The ab initio potential energy data were calculated at the second order Møller‐Plesset (MP2) perturbation theory with Dunning's correlation consistent basis sets (up to aug‐cc‐pVQZ). We considered 17 configurations of the DME dime for the orientation sampling. The calculated MP2 potential data were used to construct a 3‐site united atom force field model. The simulation results are compared with those using the empirical force field of Jorgensen and Ibrahim (Jorgensen and Ibrahim, J Am Chem Soc 1981, 103, 3976) and with available experimental measurements. We obtain quantitative agreements for the atom‐wise radial distribution functions, the self‐diffusion coefficients, and the shear viscosities over a wide range of experimental conditions. This force field thus provides a suitable starting point to predict liquid properties of DME from first principles intermolecular interactions with no empirical data input a priori. © 2012 Wiley Periodicals, Inc.     

Coupling effects of high temperature and pressure on the decomposition mechanisms of 1,1-diamino-2,2-dinitroehethe crystal: Ab initio molecular dynamics simulations

Recently, Dreger et al. experimentally investigated the phase diagram and decomposition of 1,1-diamino-2,2-dinitroethene (FOX-7) single crystal compressed hydrostatically up to 10 GPa and heated over a range of 293–750 K (J. Phys. Chem. C 2016 , 120, 11092–11098). As a continuation, we performed ab initio molecular dynamic simulations to study the initiation mechanisms and subsequent decomposition of FOX-7 at a temperature of 504 K (initial decomposition temperature) coupled with a pressure of 1–5 GPa, 604 K at 5GPa, and 704 K at 5 GPa. However, our two compressing ways are different: the former is static hydrostatical compression, while our way is dynamic compression. Our results indicate that the initial decomposition mechanism was dependent on the temperature but independent of the pressure. The initial decomposition step is the bimolecular intermolecular hydrogen transfer. The subsequent decomposition of FOX-7 is sensitive to both the temperature and pressure. At 504 K, the decomposition of FOX-7 was accelerated from 1 to 2 GPa and from 3 to 5 GPa but decelerated from 2 to 3 GPa. The temperature exhibits a positive effect on the decomposition. Overall, the temperature and pressure have great cooperative effects on the decomposition of FOX-7. Our study may provide new insight into understanding the initial mechanisms and decomposition reactions of energetic materials at relatively low temperatures coupled with different pressures in atomic detail.     

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.     

Increasing the time step with mass scaling in Born‐Oppenheimer ab initio QM/MM molecular dynamics simulations

Born‐Oppenheimer ab initio QM/MM molecular dynamics simulation with umbrella sampling is a state‐of‐the‐art approach to calculate free energy profiles of chemical reactions in complex systems. To further improve its computational efficiency, a mass‐scaling method with the increased time step in MD simulations has been explored and tested. It is found that by increasing the hydrogen mass to 10 amu, a time step of 3 fs can be employed in ab initio QM/MM MD simulations. In all our three test cases, including two solution reactions and one enzyme reaction, the resulted reaction free energy profiles with 3 fs time step and mass scaling are found to be in excellent agreement with the corresponding simulation results using 1 fs time step and the normal mass. These results indicate that for Born‐Oppenheimer ab initio QM/MM molecular dynamics simulations with umbrella sampling, the mass‐scaling method can significantly reduce its computational cost while has little effect on the calculated free energy profiles. © 2009 Wiley Periodicals, Inc. J Comput Chem, 2009     

Insights into the dynamics of evaporation and proton migration in protonated water clusters from Large‐scale Born–Oppenheimer direct dynamics

Large‐scale on‐the‐fly Born–Oppenheimer molecular dynamics simulations using recent advances in linear scaling electronic structure theory and trajectory integration techniques have been performed for protonated water clusters around the magic number (H₂O)_nH⁺, for n = 20 and 21. Besides demonstrating the feasibility and efficiency of the computational approach, the calculations reveal interesting dynamical details. Elimination of water molecules is found to be fast for both cluster sizes but rather insensitive to the initial geometry. The water molecules released acquire velocities compatible with thermal energies. The proton solvation shell changes between the well‐known Eigen and Zundel motifs and is characterized by specific low‐frequency vibrational modes, which have been quantified. The proton transfer mechanism largely resembles that of bulk water but one interesting variation was observed. © 2012 Wiley Periodicals, Inc.     

Ab initio molecular dynamics study of prebiotic production processes of organic compounds at meteorite impacts on ocean

Recent experiments concerning prebiotic materials syntheses suggest that the iron-bearing meteorite impacts on ocean during Late Heavy Bombardment provided abundant organic compounds associated with biomolecules such as amino acids and nucleobases. However, the molecular mechanism of a series of chemical reactions to produce such compounds is not well understood. In this study, we simulate the shock compression state of a meteorite impact for a model system composed of CO₂, H₂O, and metallic iron slab by ab initio molecular dynamics combined with multiscale shock technique, and clarify possible elementary reaction processes up to production of organic compounds. The reactions included not only pathways similar to the Fischer–Tropsch process known as an important hydrocarbon synthesis in many planetary processes but also those resulting in production of a carboxylic acid. It is also found that bicarbonate ions formed from CO₂ and H₂O participated in some forms in most of these observed elementary reaction processes. These findings would deepen the understanding of the full range of chemical reactions that could occur in the meteorite impact events. © 2018 Wiley Periodicals, Inc.     

The mechanism of formamide hydrolysis in water from ab initio calculations and simulations

The neutral hydrolysis of formamide in water is a suitable reference to quantify the efficiency of proteolytic enzymes. However, experimental data for this reaction has only very recently been obtained and the kinetic constant determined experimentally is significantly higher than that predicted by previous theoretical estimations. In this work, we have investigated in detail the possible mechanisms of this reaction. Several solvent models have been considered that represent a considerable improvement on those used in previous studies. Density functional and ab initio calculations have been carried out on a system which explicitly includes the first solvation shell of the formamide molecule. Its interaction with the bulk has been treated with the aid of a dielectric continuum model. Molecular dynamics simulations at the combined density functional/molecular mechanics level have been carried out in parallel to better understand the structure of the reaction intermediates in aqueous solution. Overall, the most favored mechanism predicted by our study involves two reaction steps. In the first step, the carbonyl group of the formamide molecule is hydrated to form a diol intermediate. The corresponding transition structure involves two water molecules. From this intermediate, a water-assisted proton transfer occurs from one of the hydroxy groups to the amino group. This reaction step may lead either to the formation of a new reaction intermediate with a marked zwitterionic character or to dissociation of the system into ammonia and formic acid. The zwitterionic intermediate dissociates quite easily but its lifetime is not negligible and it could play a role in the hydrolysis of substituted amides or peptides. The predicted pseudo-first-order kinetic constant for the rate-limiting step (the first step) of the hydrolysis reaction at 25 degrees C (3.9x10(-10) s(-1)) is in excellent agreement with experimental data (1.1x10(-10) s(-1)).     

Characterization of dynamics and reactivities of solvated ions by ab initio simulations

Based on a systematic investigation of trajectories of ab initio quantum mechanical/molecular mechanical simulations of numerous cations in water a standardized procedure for the evaluation of mean ligand residence times is proposed. For the characterization of reactivity and structure-breaking/structure-forming properties of the ions a measure is derived from the mean residence times calculated with different time limits. It is shown that ab initio simulations can provide much insight into ultrafast dynamics that are presently not easily accessible by experiment.     

Chemical mechanism of surface‐enhanced raman scattering spectrum of pyridine adsorbed on Ag cluster: Ab initio molecular dynamics approach

The surface‐enhanced Raman scattering (SERS) spectrum of pyridine adsorbed on Ag₂₀ cluster (pyridine‐Ag₂₀) at room temperature is calculated by performing ab initio molecular dynamics simulations in connection with a Fourier transform of the polarizability autocorrelation function to investigate the static chemical enhancement behind the SERS spectrum. The five enhanced vibrational modes of pyridine, namely, υ_6a, υ₁, υ₁₂, υ_9a, and υ_8a, can be assigned and identified by using a new analytical scheme, namely, single‐frequency‐pass filter, which is based on a Fourier transform filtering technique. To understand the factors evoking the enhancement in the SERS spectrum, the dynamic properties of molecular structures and charges for both of the free pyridine and adsorbed pyridine are analyzed. The calculated results indicate that the vibrational amplitudes of adsorbed pyridine are enhanced due to both of the electron transfer from pyridine to Ag₂₀ cluster and the softening of pyridine bond. In addition, the N‐Ag stretching within pyridine‐Ag₂₀ will couple with these five vibrational modes of pyridine. Consequently, the electron transfer between pyridine and Ag₂₀ cluster induced by different molecular vibrational modes prompts the redistribution of electron density of pyridine. These factors collectively cause the noticeable change in polarizability during molecular vibrations and hence result in the enhancement of Raman peaks. © 2013 Wiley Periodicals, Inc.