A combined molecular dynamics+quantum mechanics method for investigation of dynamic effects on local surface structures

A combined molecular dynamics (MD)+quantum mechanics (QM) method for studying processes on ionic surfaces is presented. Through the combination of classical MD and ab initio embedded-cluster calculations, this method allows the modeling of surface processes involving both the structural and dynamic features of the substrate, even for large-scale systems. The embedding approach used to link the information from the MD simulation to the cluster calculation is presented, and rigorous tests have been carried out to ensure the feasibility of the method. The electrostatic potential and electron density resulting from our embedded-cluster model have been compared with periodic slab results, and confirm the satisfying quality of our embedding scheme as well as the importance of applying embedding in our combined MD+QM approach. We show that a highly accurate representation of the Madelung potential becomes a prerequisite when the embedded-cluster approach is applied to temperature-distorted surface snapshots from the MD simulation.     

O2 interaction and reactivity on a model hydroxylated rutile(110) surface

Recently several theoretical studies have examined oxygen adsorption on the clean, reduced TiO2(110) surface. However the photocatalytic behavior of TiO2 and the scavenging ability of oxygen are known to be influenced by the presence of surface hydroxyls. In this paper the chemistry of O2 on the hydroxylated TiO2 surface is investigated by means of first-principles total energy calculations and molecular dynamics (MD) simulations. The MD trajectories show a direct, spontaneous reaction between O2 and the surface hydroxyls, thus supporting the experimental hypothesis that the reaction does not necessarily pass through a chemisorbed O2 state. Following this reaction, the most stable chemisorbed intermediates are found to be peroxide species HO2 and H2O2. Although these intermediates are very stable on the short time scale of MD simulations, the energetics suggests that their further transformation is connected to a new 300 K feature observed in the experimental water temperature programmed desorption (TPD) spectrum. The participation of two less stable intermediate states, involving terminal hydroxyls and/or chemisorbed water plus oxygen adatoms, to the desorption process, is not supported by the total energy calculations. Analysis of the projected density of states, however, suggests the possibility that these intermediates have a role in completing the surface oxidation immediately before desorption.     

Hydrogen-bond dynamics and Fermi resonance in high-pressure methane filled ice

High-pressure, variable temperature infrared spectroscopy and first-principles calculations on the methane filled ice structure (MH-III) at high pressures are used to investigate the vibrational dynamics related to pressure induced modifications in hydrogen bonding. Infrared spectroscopy of isotopically dilute solutions of H(2)O in D(2)O is employed together with first-principles calculations to characterize proton dynamics with the pressure induced shortening of hydrogen bonds. A Fermi resonance is identified and shown to dominate the infrared spectrum in the pressure region between 10 and 30 GPa. Significant differences in the effects of the Fermi resonance observed between 10 and 300 K arise from the double-well potential energy surface of the hydrogen bond and quantum effects associated with the proton dynamics.     

An automated approach for the parameterization of accurate intermolecular force‐fields: Pyridine as a case study

An automated protocol is proposed and validated, which integrates accurate quantum mechanical calculations with classical numerical simulations. Intermolecular force fields, (FF) suitable for molecular dynamics (MD) and Monte Carlo simulations, are parameterized through a novel iterative approach, fully based on quantum mechanical data, which has been automated and coded into the PICKY software, here presented. The whole procedure is tested and validated for pyridine, whose bulk phase, described through MD simulations performed with the specifically parameterized FF, is characterized by computing several of its thermodynamic, structural, and transport properties, comparing them with their experimental counterparts. © 2011 Wiley Periodicals, Inc.     

Finite‐field method with unbiased polarizable continuum model for evaluation of the second hyperpolarizability of an open‐shell singlet molecule in solvents

The static second hyperpolarizability γ of the complexes composed of open‐shell singlet 1,3‐dipole molecule involving a boron atom and a water molecule in aqueous phase are investigated by the finite‐field (FF) method combined with a standard polarized continuum model (PCM) and with a newly proposed unbiased PCM (UBPCM). On the basis of the comparison with the results calculated by the FF method using the full quantum and the quantum‐mechanical/molecular‐mechanical and molecular‐dynamics (QM/MM‐MD) treatments, the present FF‐UBPCM method is demonstrated to remedy the artificial overestimation of the γ caused by standard FF‐PCM calculations and to well reproduce the FF‐QM/MM‐MD and FF‐full‐QM results with much lower costs. © 2013 Wiley Periodicals, Inc.     

Adsorption of atomic and molecular oxygen on the LaMnO3(001) surface: ab initio supercell calculations and thermodynamics

We present and discuss the results of ab initio DFT plane-wave supercell calculations of the atomic and molecular oxygen adsorption and diffusion on the LaMnO(3) (001) surface which serves as a model material for a cathode of solid oxide fuel cells. The dissociative adsorption of O(2) molecules from the gas phase is energetically favorable on surface Mn ions even on a defect-free surface. The surface migration energy for adsorbed O ions is found to be quite high, 2.0 eV. We predict that the adsorbed O atoms could penetrate the electrode first plane when much more mobile surface oxygen vacancies (migration energy of 0.69 eV) approach the O ions strongly bound to the surface Mn ions. The formation of the O vacancy near the O atom adsorbed atop surface Mn ion leads to an increase of the O-Mn binding energy by 0.74 eV whereas the drop of this adsorbed O atom into a vacancy possesses no energy barrier. Ab initio thermodynamics predicts that at typical SOFC operation temperatures (approximately 1200 K) the MnO(2) (001) surface with adsorbed O atoms is the most stable in a very wide range of oxygen gas pressures (above 10(-2) atm).     

氧化铜表面臭氧分解路径及表面氧物种生成机理研究

基于密度泛函理论（DFT）计算研究了O3在完整和具有氧空位的CuO(111)表面吸附的吸附位、吸附结构、吸附能和电子转移情况,比较了O3在完整表面和具有氧空位的表面分解的路径和能垒,分析了氧空位和表面吸附氧的生成机理。结果表明,在完整CuO表面,O3分子通过化学吸附或物理吸附表面结合,吸附能最高为-1.22eV（构型bri(2)）。O3在具有氧空位的CuO表面均为化学吸附,吸附能最高为-2.95eV（构型ovbri(3)）,显著高于完整表面的吸附能。O3吸附后,Cu吸附位的电荷密度减小,O3中的O原子附近的电荷密度显著增强,电荷从CuO表面转移到O3,并形成Cu-O离子键。O3分解后形成了超氧物种,提高了表面的氧化活性。在完整表面,以构型bri(2)为起始构型的路径反应能垒最低,为0.52eV;O2*在完整表面的脱附所需要的最低能量为0.42eV,形成氧空位的O2*脱附能为2.06eV。在具有氧空位的表面,O3分解的反应能垒为0.30eV（构型ovbri(1)）和0.12eV（构型ovbri(3)）,均低于完整表面的反应能垒;分解形成的O2*的最低脱附能也低于完整表面,为0.27eV。可见,氧空位的形成提高了吸附能,降低了反应能垒,使O3分子更容易吸附在CuO表面,并加快了O3的催化分解。     

New force field for calcium binding sites in annexin-membrane complexes

For accurate classical molecular dynamics (MD) simulations of the calcium mediated bound complexes of annexin and membrane we have developed new force-field parameters correctly describing the interaction of the Ca ion with its environment. We have used quantum chemical calculations to investigate the potential energy surface experienced by the Ca ion within the three different binding sites found in domain 1 of annexin V (ANX V/1). Based on these calculations we were able to quantify the charge polarization of atoms within the binding sites, and to determine the geometry and force constants of harmonic restraints between the Ca ion and its coordinating oxygen atoms. Harmonic restraints were introduced to compensate for the deviations between the quantum mechanical potential energy surface and that of the classical force field. Our analysis has shown that using the refined force field for the Ca binding sites enables long-time MD simulations that conserve the initial structure of ANX V/1 significantly better than MD simulations using the standard force field.     

Nonstatistical behavior of reactive scattering in the (18)O+(32)O(2) isotope exchange reaction

The recombination of oxygen atoms with oxygen molecules to form ozone exhibits several strange chemical characteristics, including unusually large differences in formation rate coefficients when different isotopes of oxygen participate. Purely statistical chemical reaction rate theories cannot describe these isotope effects, suggesting that reaction dynamics must play an important role. We investigated the dynamics of the 18O + 32O2 --> O3(*) --> 16O + 34O2 isotope exchange reaction (which proceeds on the same potential energy surface as ozone formation) using crossed atomic and molecular beams at a collision energy of 7.3 kcal mol(-1), providing the first direct experimental evidence that the dissociation of excited ozone exhibits significant nonstatistical behavior. These results are compared with quantum statistical and quasi-classical trajectory calculations in order to gain insight into the potential energy surface and the dynamics of ozone formation.     

Promoting vibrations in human purine nucleoside phosphorylase. A molecular dynamics and hybrid quantum mechanical/molecular mechanical study

Crystallographic studies of human purine nucleoside phosphorylase (hPNP) with several transition-state (TS) analogues in the immucillin family showed an unusual geometric arrangement of the atoms O-5', O-4', and O(P), the nucleophilic phosphate oxygen, lying in a close three-oxygen stack. These observations were corroborated by extensive experimental kinetic isotope effect analysis. We propose that protein-facilitated dynamic modes in hPNP cause this stack, centered on the ribosyl O-4' oxygen, to squeeze together and push electrons toward the purine ring, stabilizing the oxacarbenium character of the TS. As the N-ribosidic bond is cleaved during the reaction, the pK(a) values of N-7 and O-6 increase by the electron density expelled by the oxygen-stack compression toward the purine ring. Increased electron density in the purine ring improves electrostatic interactions with nearby residues and facilitates the abstraction of a proton from a solvent proton or an unidentified general acid, making the purine a better leaving group, and accelerating catalysis. Classical and mixed quantum/classical molecular dynamics (MD) simulations of the Michaelis complex of hPNP with the substrates guanosine and phosphate were performed to assess the existence of protein-promoting vibrations (PPVs). Analogous simulations were performed for the substrates in aqueous solution. In the catalytic site, the O-5', O-4', and O(P) oxygens vibrate at frequencies of ca. 125 and 465 cm(-1), as opposed to 285 cm(-1) in the absence of hPNP. The hybrid quantum mechanical/molecular mechanical method was used to assess whether this enzymatic vibration pushing the oxygens together is coupled to the reaction coordinate, and thus has a direct positive impact on catalysis. The potential energy surface for the phosphorolysis reaction for several snapshots taken from the classical MD simulation showed substantial differences in oxygen compression. Our calculations showed the existence of PPVs coupled to the reaction coordinate, which effect electronic alterations in the active site by pushing the three oxygen centers together in proximity, and accelerate substrate turnover in the phosphorolysis reaction catalyzed by hPNP.     

Optical absorption of methoxy and carboethoxy derivatives of 1,3-diphenyl-1H-pyrazolo[3,4-b]quinoline

Paper deals with experimental investigations and quantum chemical calculations of the optical absorption spectra of methoxy and carboethoxy 1,3-diphenyl derivatives of the pyrazoloquinoline ([PQ]): 6-methoxy-1,3-dyphenil-[PQ], 6-methoxy-1,3-(p-methoxyphenyl)-[PQ], 6-methoxy-1-(p-methoxyphenyl)-[PQ] and 6-carboethoxy-1,3-diphenyl-[PQ]. The quantum chemical calculations are performed by means of the semiempirical quantum chemical methods (AM1 or PM3) applied to: (a) the equilibrium molecular conformation in vacuo (T=0 K); (b) the molecular dynamic (MD) trajectory (T=300 K) which includes the dynamics of a certain molecular fragment (moiety) only (fragmental MD simulations); or (c) the MD trajectory obtained for most general case within the total MD simulations at T=300 K. The results of these calculations are compared with the measured spectra of the optical absorption. The quantum chemical simulations show that the dynamics of the methoxy or carboethoxy groups practically does not influence the absorption spectrum whereas the strongest its modification (300相似文献   

QM/MM calculation of solvent effects on absorption spectra of guanine

Electronic spectra of guanine in the gas phase and in water were studied by quantum mechanical/molecular mechanical (QM/MM) methods. Geometries for the excited‐state calculations were extracted from ground‐state molecular dynamics (MD) simulations using the self‐consistent‐charge density functional tight binding (SCC‐DFTB) method for the QM region and the TIP3P force field for the water environment. Theoretical absorption spectra were generated from excitation energies and oscillator strengths calculated for 50 to 500 MD snapshots of guanine in the gas phase (QM) and in solution (QM/MM). The excited‐state calculations used time‐dependent density functional theory (TDDFT) and the DFT‐based multireference configuration interaction (DFT/MRCI) method of Grimme and Waletzke, in combination with two basis sets. Our investigation covered keto‐N7H and keto‐N9H guanine, with particular focus on solvent effects in the low‐energy spectrum of the keto‐N9H tautomer. When compared with the vertical excitation energies of gas‐phase guanine at the optimized DFT (B3LYP/TZVP) geometry, the maxima in the computed solution spectra are shifted by several tenths of an eV. Three effects contribute: the use of SCC‐DFTB‐based rather than B3LYP‐based geometries in the MD snapshots (red shift of ca. 0.1 eV), explicit inclusion of nuclear motion through the MD snapshots (red shift of ca. 0.1 eV), and intrinsic solvent effects (differences in the absorption maxima in the computed gas‐phase and solution spectra, typically ca. 0.1–0.3 eV). A detailed analysis of the results indicates that the intrinsic solvent effects arise both from solvent‐induced structural changes and from electrostatic solute–solvent interactions, the latter being dominant. © 2009 Wiley Periodicals, Inc. J Comput Chem 2010     

Communication: Quantum polarized fluctuating charge model: a practical method to include ligand polarizability in biomolecular simulations

We present a simple and practical method to include ligand electronic polarization in molecular dynamics (MD) simulation of biomolecular systems. The method involves periodically spawning quantum mechanical (QM) electrostatic potential (ESP) calculations on an extra set of computer processors using molecular coordinate snapshots from a running parallel MD simulation. The QM ESPs are evaluated for the small-molecule ligand in the presence of the electric field induced by the protein, solvent, and ion charges within the MD snapshot. Partial charges on ligand atom centers are fit through the multi-conformer restrained electrostatic potential (RESP) fit method on several successive ESPs. The RESP method was selected since it produces charges consistent with the AMBER/GAFF force-field used in the simulations. The updated charges are introduced back into the running simulation when the next snapshot is saved. The result is a simulation whose ligand partial charges continuously respond in real-time to the short-term mean electrostatic field of the evolving environment without incurring additional wall-clock time. We show that (1) by incorporating the cost of polarization back into the potential energy of the MD simulation, the algorithm conserves energy when run in the microcanonical ensemble and (2) the mean solvation free energies for 15 neutral amino acid side chains calculated with the quantum polarized fluctuating charge method and thermodynamic integration agree better with experiment relative to the Amber fixed charge force-field.     

QuanPol: A full spectrum and seamless QM/MM program

The quantum chemistry polarizable force field program (QuanPol) is implemented to perform combined quantum mechanical and molecular mechanical (QM/MM) calculations with induced dipole polarizable force fields and induced surface charge continuum solvation models. The QM methods include Hartree–Fock method, density functional theory method (DFT), generalized valence bond theory method, multiconfiguration self‐consistent field method, Møller–Plesset perturbation theory method, and time‐dependent DFT method. The induced dipoles of the MM atoms and the induced surface charges of the continuum solvation model are self‐consistently and variationally determined together with the QM wavefunction. The MM force field methods can be user specified, or a standard force field such as MMFF94, Chemistry at Harvard Molecular Mechanics (CHARMM), Assisted Model Building with Energy Refinement (AMBER), and Optimized Potentials for Liquid Simulations‐All Atom (OPLS‐AA). Analytic gradients for all of these methods are implemented so geometry optimization and molecular dynamics (MD) simulation can be performed. MD free energy perturbation and umbrella sampling methods are also implemented. © 2013 Wiley Periodicals, Inc.     

Path-integral molecular dynamics simulations of hydrated hydrogen chloride cluster HCl(H2O)4 on a semiempirical potential energy surface

Path-integral molecular dynamics simulations for the HCl(H₂O)₄ cluster have been performed on the ground-state potential energy surface directly obtained on-the-fly from semiempirical PM3-MAIS molecular orbital calculations. It is found that the HCl(H₂O)₄ cluster has structural rearrangement above the temperature of 300 K showing a liquid-like behavior. Quantum mechanical fluctuation of hydrogen nuclei plays a significant role in structural arrangement processes in this cluster.     

First principle and ReaxFF molecular dynamics investigations of formaldehyde dissociation on Fe(100) surface

Detailed formaldehyde adsorption and dissociation reactions on Fe(100) surface were studied using first principle calculations and molecular dynamics (MD) simulations, and results were compared with available experimental data. The study includes formaldehyde, formyl radical (HCO), and CO adsorption and dissociation energy calculations on the surface, adsorbate vibrational frequency calculations, density of states analysis of clean and adsorbed surfaces, complete potential energy diagram construction from formaldehyde to atomic carbon (C), hydrogen (H), and oxygen (O), simulation of formaldehyde adsorption and dissociation reaction on the surface using reactive force field, ReaxFF MD, and reaction rate calculations of adsorbates using transition state theory (TST). Formaldehyde and HCO were adsorbed most strongly at the hollow (fourfold) site. Adsorption energies ranged from ?22.9 to ?33.9 kcal/mol for formaldehyde, and from ?44.3 to ?66.3 kcal/mol for HCO, depending on adsorption sites and molecular direction. The dissociation energies were investigated for the dissociation paths: formaldehyde → HCO + H, HCO → H + CO, and CO → C + O, and the calculated energies were 11.0, 4.1, and 26.3 kcal/mol, respectively. ReaxFF MD simulation results were compared with experimental surface analysis using high resolution electron energy loss spectrometry (HREELS) and TST based reaction rates. ReaxFF simulation showed less reactivity than HREELS observation at 310 and 523 K. ReaxFF simulation showed more reactivity than the TST based rate for formaldehyde dissociation and less reactivity than TST based rate for HCO dissociation at 523 K. TST‐based rates are consistent with HREELS observation. © 2013 Wiley Periodicals, Inc.     

Mechanically induced generation of highly reactive excited-state oxygen molecules in cluster scattering

Molecular electronic excitation in (O(2))(n) clusters induced by mechanical collisions via the "chemistry with a hammer" is investigated by a combination of molecular dynamics simulations and quantum chemistry calculations. Complete active space self-consistent field augmented with triple-zeta polarizable basis set quantum chemistry calculations of a compressed (O(2))(2) cluster model in various configurations reveal the emergence of possible pathways for the generation of electronically excited singlet O(2) molecules upon cluster compression and vibrational excitation, due to electronic curve-crossing and spin-orbit coupling. Extrapolation of the model (O(2))(2) results to larger clusters suggests a dramatic increase in the population of electronically excited O(2) products, and may account for the recently observed cluster-catalyzed oxidation of silicon surfaces, via singlet oxygen generation induced by cluster impact, followed by surface reaction of highly reactive singlet O(2) molecules. Extensive molecular dynamics simulations of (O(2))(n) clusters colliding onto a hot surface indeed reveal that cluster compression is sufficient under typical experimental conditions for nonadiabatic transitions to occur. This work highlights the importance of nonadiabatic effects in the "chemistry with a hammer."     

Tyrosine absorption spectroscopy: Backbone protonation effects on the side chain electronic properties

The UV–vis spectrum of Tyrosine and its response to different backbone protonation states have been studied by applying the Perturbed Matrix Method (PMM) in conjunction with molecular dynamics (MD) simulations. Herein, we theoretically reproduce the UV–vis absorption spectrum of aqueous solution of Tyrosine in its zwitterionic, anionic and cationic forms, as well as of aqua‐p‐Cresol (i.e., the moiety that constitutes the side chain portion of Tyrosine). To achieve a better accuracy in the MD sampling, the Tyrosine Force Field (FF) parameters were derived de novo via quantum mechanical calculations. The UV–vis absorption spectra are computed considering the occurring electronic transitions in the vertical approximation for each of the chromophore configurations sampled by the classical MD simulations, thus including the effects of the chromophore semiclassical structural fluctuations. Finally, the explicit treatment of the perturbing effect of the embedding environment permits to fully model the inhomogeneous bandwidth of the electronic spectra. Comparison between our theoretical–computational results and experimental data shows that the used model captures the essential features of the spectroscopic process, thus allowing to perform further analysis on the strict relationship between the quantum properties of the chromophore and the different embedding environments. © 2018 Wiley Periodicals, Inc.     

Communication: highly accurate ozone formation potential and implications for kinetics

Atmospheric ozone is formed by the O + O(2) exchange reaction followed by collisional stabilization of the O(3)(?) intermediate. The dynamics of the O + O(2) reaction and to a lesser extent the O(3) stabilization depend sensitively on the underlying potential energy surface, particularly in the asymptotic region. Highly accurate Davidson corrected multi-state multi-reference configuration interaction calculations reported here reveal that the minimal energy path for the formation of O(3) from O + O(2) is a monotonically decaying function of the atom-diatom distance and contains no "reef" feature found in previous ab initio calculations. The absence of a submerged barrier leads to an exchange rate constant with the correct temperature dependence and is in better agreement with experiment, as shown by quantum scattering calculations.