Essential dynamics for the study of microstructures in liquids

Essential Dynamics (ED) is a powerful tool for analyzing molecular dynamics (MD) simulations and it is widely adopted for conformational analysis of large molecular systems such as, for example, proteins and nucleic acids. In this study, we extend the use of ED to the study of clusters of arbitrary size constituted by weakly interacting particles, for example, atomic clusters and supramolecular systems. The key feature of the method we present is the identification of the relevant atomic‐molecular clusters to be analyzed by ED for extracting the information of interest. The application of this computational approach allows a straightforward and unbiased conformational study of the local microstructures in liquids, as emerged from semiclassical MD simulations. The good performance of the method is demonstrated by calculating typical observables of liquid water, that is, NMR, NEXAFS O1s, and IR spectra, known to be rather sensitive both to the presence and to the conformational features of hydrogen‐bonded clusters. © 2014 Wiley Periodicals, Inc.

[1]  Bernd G. Pfrommer,et al.  NMR chemical shifts of ice and liquid water: The effects of condensation , 2000 .

[2]  W. Kulig,et al.  A 'clusters-in-liquid' method for calculating infrared spectra identifies the proton-transfer mode in acidic aqueous solutions. , 2013, Nature chemistry.

[3]  Andrea Amadei,et al.  Essential dynamics: foundation and applications , 2012 .

[4]  T. Darden,et al.  Particle mesh Ewald: An N⋅log(N) method for Ewald sums in large systems , 1993 .

[5]  Y. Maréchal.  Infrared spectra of water. I. Effect of temperature and of H/D isotopic dilution , 1991 .

[6]  Peter S. Shenkin,et al.  Cluster analysis of molecular conformations , 1994, J. Comput. Chem..

[7]  A. Becke,et al.  Density-functional exchange-energy approximation with correct asymptotic behavior. , 1988, Physical review. A, General physics.

[8]  Trygve Helgaker,et al.  Ab Initio Methods for the Calculation of NMR Shielding and Indirect Spin-Spin Coupling Constants , 1999 .

[9]  Walter Thiel,et al.  QM/MM calculation of solvent effects on absorption spectra of guanine , 2010, J. Comput. Chem..

[10]  Vincenzo Barone,et al.  Solvent effects on the UV (n --> pi*) and NMR (13C and 17O) spectra of acetone in aqueous solution. An integrated car-parrinello and DFT/PCM approach. , 2005, The journal of physical chemistry. B.

[11]  P. Wernet,et al.  The Structure of the First Coordination Shell in Liquid Water , 2004, Science.

[12]  H. Berendsen,et al.  Interaction Models for Water in Relation to Protein Hydration , 1981 .

[13]  S. Franzen,et al.  Infrared spectra of H(2)16O, H(2)18O and D(2)O in the liquid phase by single-pass attenuated total internal reflection spectroscopy. , 2004, Spectrochimica acta. Part A, Molecular and biomolecular spectroscopy.

[14]  E. Fileti,et al.  Isotropic and anisotropic NMR chemical shifts in liquid water: a sequential QM/MM study , 2007 .

[15]  Walter Kauzmann,et al.  The Structure and Properties of Water , 1969 .

[16]  Hans Ågren,et al.  Calculations of near-edge x-ray-absorption spectra of gas-phase and chemisorbed molecules by means of density-functional and transition-potential theory , 1998 .

[17]  Evert Jan Baerends,et al.  Self-consistent molecular Hartree—Fock—Slater calculations I. The computational procedure , 1973 .

[18]  Collective Coordinates Kinetics and Dynamic Structural Factor in Classical Liquids , 1995 .

[19]  J. G. Snijders,et al.  Towards an order-N DFT method , 1998 .

[20]  H. Berendsen,et al.  Essential dynamics of proteins , 1993, Proteins.

[21]  M. Karplus,et al.  Collective motions in proteins: A covariance analysis of atomic fluctuations in molecular dynamics and normal mode simulations , 1991, Proteins.

[22]  Berk Hess,et al.  LINCS: A linear constraint solver for molecular simulations , 1997 .

[23]  D. Bohm,et al.  COLLECTIVE TREATMENT OF LIQUID HELIUM. , 1967 .

[24]  M. Odelius,et al.  X-Ray Absorption Spectrum of Liquid Water From Molecular Dynamics Simulations: Asymmetric Model , 2006 .

[25]  Luca Frediani,et al.  The Dalton quantum chemistry program system , 2013, Wiley interdisciplinary reviews. Computational molecular science.

[26]  J. Hasted,et al.  The Structure and Properties of Water , 1969 .

[27]  Theoretical modeling and interpretation of X-ray absorption spectra of liquid water. , 2006, Physical chemistry chemical physics : PCCP.

[28]  Arvi Rauk,et al.  On the calculation of multiplet energies by the hartree-fock-slater method , 1977 .

[29]  R. C. Cohen,et al.  Probing the local structure of liquid water by X-ray absorption spectroscopy. , 2006, The journal of physical chemistry. B.

[30]  Parr,et al.  Development of the Colle-Salvetti correlation-energy formula into a functional of the electron density. , 1988, Physical review. B, Condensed matter.

[31]  A. Amadei,et al.  Entropy–energy balance in base catalyzed keto-enol interconversion: A joint theoretical and experimental investigation† , 2011 .

[32]  A. Amadei,et al.  Theoretical characterization of temperature and density dependence of liquid water electronic excitation energy: comparison with recent experimental data. , 2008, The Journal of chemical physics.

[33]  V. Barone,et al.  Accurate density functional calculations of near-edge x-ray and optical absorption spectra of liquid water using nonperiodic boundary conditions: the role of self-interaction and long-range effects. , 2008, Physical review letters.

[34]  Gerrit Groenhof,et al.  GROMACS: Fast, flexible, and free , 2005, J. Comput. Chem..

[35]  B. Halle,et al.  Proton magnetic shielding tensor in liquid water. , 2002, Journal of the American Chemical Society.

[36]  D. Sebastiani,et al.  NMR chemical shifts as a tool to analyze first principles molecular dynamics simulations in condensed phases: the case of liquid water , 2010, Magnetic resonance in chemistry : MRC.

[37]  J. Perdew,et al.  Erratum: Density-functional approximation for the correlation energy of the inhomogeneous electron gas , 1986, Physical review. B, Condensed matter.

[38]  Andrew E. Torda,et al.  Algorithms for clustering molecular dynamics configurations , 1994, J. Comput. Chem..

[39]  Yoshio Inoue,et al.  General parameterization of a reaction field theory combined with the boundary element method , 1994, J. Comput. Chem..

[40]  M. Chiarini,et al.  Inclusion of cybotactic effect in the theoretical modeling of absorption spectra of liquid-state systems with perturbed matrix method and molecular dynamics simulations: the UV–Vis absorption spectrum of para-nitroaniline as a case study , 2014, Theoretical Chemistry Accounts.

[41]  H. Berendsen,et al.  Molecular dynamics with coupling to an external bath , 1984 .

[42]  Jared D. Smith,et al.  Revisiting the total ion yield x-ray absorption spectra of liquid water microjets , 2008, Journal of physics. Condensed matter : an Institute of Physics journal.

[43]  Yi Luo,et al.  LETTER TO THE EDITOR: Spectroscopic probing of local hydrogen-bonding structures in liquid water , 2002 .