51
|
Jung KA, Videla PE, Batista VS. Multi-time formulation of Matsubara dynamics. J Chem Phys 2019; 151:034108. [PMID: 31325942 DOI: 10.1063/1.5110427] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Matsubara dynamics has recently emerged as the most general form of a quantum-Boltzmann-conserving classical dynamics theory for the calculation of single-time correlation functions. Here, we present a generalization of Matsubara dynamics for the evaluation of multitime correlation functions. We first show that the Matsubara approximation can also be used to approximate the two-time symmetrized double Kubo transformed correlation function. By a straightforward extension of these ideas to the multitime realm, a multitime Matsubara dynamics approximation can be obtained for the multitime fully symmetrized Kubo transformed correlation function. Although not a practical method, due to the presence of a phase-term, this multitime formulation of Matsubara dynamics represents a benchmark theory for future development of Boltzmann preserving semiclassical approximations to general higher order multitime correlation functions. It also reveals a connection between imaginary time-ordering in the path integral and the classical dynamics of multitime correlation functions.
Collapse
Affiliation(s)
- Kenneth A Jung
- Department of Chemistry, Yale University, P.O. Box 208107, New Haven, Connecticut 06520-8107, USA
| | - Pablo E Videla
- Department of Chemistry, Yale University, P.O. Box 208107, New Haven, Connecticut 06520-8107, USA
| | - Victor S Batista
- Department of Chemistry, Yale University, P.O. Box 208107, New Haven, Connecticut 06520-8107, USA
| |
Collapse
|
52
|
Kaur R, Welsch R. Probing photodissociation dynamics using ring polymer molecular dynamics. J Chem Phys 2019; 150:114105. [PMID: 30901996 DOI: 10.1063/1.5086218] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Abstract
The performance of the ring polymer molecular dynamics (RPMD) approach to simulate typical photodissociation processes is assessed. The correct description of photodissociation requires the calculation of correlation functions or expectation values associated with non-equilibrium initial conditions, which was shown to be possible with RPMD very recently [J. Chem. Phys. 145, 204118 (2016)]. This approach is combined with treatment of the nonadiabatic dynamics employing the ring polymer surface hopping approach (RPSH), which is based on Tully's fewest switches surface hopping (FSSH) approach. The performance is tested using one-dimensional photodissociation models. It is found that RPSH with non-equilibrium initial conditions can well reproduce the time-dependent dissociation probability, and adiabatic and diabatic populations for cases where the crossing point is below and above the Franck-Condon point, respectively, while standard FSSH fails to reproduce the exact quantum dynamics in the latter case. Thus, it is shown that RPSH is an efficient and accurate alternative to standard FSSH, which is one of the most widely employed approaches to study photochemistry.
Collapse
Affiliation(s)
- Rajwant Kaur
- Center for Free-Electron Laser Science, DESY, Notkestrasse 85, 22607 Hamburg, Germany
| | - Ralph Welsch
- Center for Free-Electron Laser Science, DESY, Notkestrasse 85, 22607 Hamburg, Germany
| |
Collapse
|
53
|
Yu Q, Bowman JM. Classical, Thermostated Ring Polymer, and Quantum VSCF/VCI Calculations of IR Spectra of H7O3+ and H9O4+ (Eigen) and Comparison with Experiment. J Phys Chem A 2019; 123:1399-1409. [DOI: 10.1021/acs.jpca.8b11603] [Citation(s) in RCA: 25] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/07/2023]
Affiliation(s)
- Qi Yu
- Cherry L. Emerson Center for Scientific Computation and Department of Chemistry, Emory University, Atlanta, Georgia 30322, United States
| | - Joel M. Bowman
- Cherry L. Emerson Center for Scientific Computation and Department of Chemistry, Emory University, Atlanta, Georgia 30322, United States
| |
Collapse
|
54
|
Vaillant CL, Althorpe SC, Wales DJ. Path Integral Energy Landscapes for Water Clusters. J Chem Theory Comput 2019; 15:33-42. [PMID: 30550261 DOI: 10.1021/acs.jctc.8b00675] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/03/2023]
Abstract
The energy landscapes for a discretized path integral representation of the water dimer, trimer and pentamer are characterized in terms of the localized (classical) and delocalized minima and transition states. The transition states are finite-temperature approximations to the exact instanton path, and they are typically used to calculate the tunneling splittings or reaction rates. The features of the path integral landscape are explored, thus elucidating procedures that could usefully be automated when searching for instantons in larger systems. Our work not only clarifies the role of minima and transition states in path integral calculations but also enables us to analyze the quantum-to-classical transition.
Collapse
Affiliation(s)
- Christophe L Vaillant
- University Chemical Laboratories , Lensfield Road , Cambridge CB2 1EW , United Kingdom
| | - Stuart C Althorpe
- University Chemical Laboratories , Lensfield Road , Cambridge CB2 1EW , United Kingdom
| | - David J Wales
- University Chemical Laboratories , Lensfield Road , Cambridge CB2 1EW , United Kingdom
| |
Collapse
|
55
|
Qu C, Bowman JM. Quantum approaches to vibrational dynamics and spectroscopy: is ease of interpretation sacrificed as rigor increases? Phys Chem Chem Phys 2019; 21:3397-3413. [DOI: 10.1039/c8cp04990d] [Citation(s) in RCA: 25] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/09/2023]
Abstract
The subject of this Perspective is quantum approaches, beyond the harmonic approximation, to vibrational dynamics and IR spectroscopy.
Collapse
Affiliation(s)
- Chen Qu
- Department of Chemistry and Cherry L. Emerson Center for Scientific Computation, Emory University
- Atlanta
- USA
| | - Joel M. Bowman
- Department of Chemistry and Cherry L. Emerson Center for Scientific Computation, Emory University
- Atlanta
- USA
| |
Collapse
|
56
|
Elton DC, Fritz M, Fernández-Serra M. Using a monomer potential energy surface to perform approximate path integral molecular dynamics simulation of ab initio water at near-zero added cost. Phys Chem Chem Phys 2018; 21:409-417. [PMID: 30534683 DOI: 10.1039/c8cp06077k] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/07/2023]
Abstract
It is now established that nuclear quantum motion plays an important role in determining water's hydrogen bonding, structure, and dynamics. Such effects are important to include in density functional theory (DFT) based molecular dynamics simulation of water. The standard way of treating nuclear quantum effects, path integral molecular dynamics (PIMD), multiplies the number of energy/force calculations by the number of beads required. In this work we introduce a method whereby PIMD can be incorporated into a DFT simulation with little extra cost and little loss in accuracy. The method is based on the many body expansion of the energy and has the benefit of including a monomer level correction to the DFT energy. Our method calculates intramolecular forces using the highly accurate monomer potential energy surface developed by Partridge-Schwenke, which is cheap to evaluate. Intermolecular forces and energies are calculated with DFT only once per timestep using the centroid positions. We show how our method may be used in conjunction with a multiple time step algorithm for an additional speedup and how it relates to ring polymer contraction and other schemes that have been introduced recently to speed up PIMD simulations. We show that our method, which we call "monomer PIMD", correctly captures changes in the structure of water found in a full PIMD simulation but at much lower computational cost.
Collapse
Affiliation(s)
- Daniel C Elton
- Department of Physics and Astronomy, Stony Brook University, Stony Brook, New York 11794-3800, USA.
| | | | | |
Collapse
|
57
|
Karsten S, Ivanov SD, Bokarev SI, Kühn O. Simulating vibronic spectra via Matsubara-like dynamics: Coping with the sign problem. J Chem Phys 2018; 149:194103. [DOI: 10.1063/1.5046874] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Sven Karsten
- Institute of Physics, University of Rostock, Albert-Einstein-Str. 23-24, 18059 Rostock, Germany
| | - Sergei D. Ivanov
- Institute of Physics, University of Rostock, Albert-Einstein-Str. 23-24, 18059 Rostock, Germany
| | - Sergey I. Bokarev
- Institute of Physics, University of Rostock, Albert-Einstein-Str. 23-24, 18059 Rostock, Germany
| | - Oliver Kühn
- Institute of Physics, University of Rostock, Albert-Einstein-Str. 23-24, 18059 Rostock, Germany
| |
Collapse
|
58
|
|
59
|
Trenins G, Althorpe SC. Mean-field Matsubara dynamics: Analysis of path-integral curvature effects in rovibrational spectra. J Chem Phys 2018; 149:014102. [DOI: 10.1063/1.5038616] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- George Trenins
- Department of Chemistry, University of Cambridge, Lensfield Road, Cambridge CB2 1EW, United Kingdom
| | - Stuart C. Althorpe
- Department of Chemistry, University of Cambridge, Lensfield Road, Cambridge CB2 1EW, United Kingdom
| |
Collapse
|
60
|
Jung KA, Videla PE, Batista VS. Inclusion of nuclear quantum effects for simulations of nonlinear spectroscopy. J Chem Phys 2018; 148:244105. [PMID: 29960352 DOI: 10.1063/1.5036768] [Citation(s) in RCA: 16] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
The computation and interpretation of nonlinear vibrational spectroscopy is of vital importance for understanding a wide range of dynamical processes in molecular systems. Here, we introduce an approach to evaluate multi-time response functions in terms of multi-time double symmetrized Kubo transformed thermal correlation functions. Furthermore, we introduce a multi-time extension of ring polymer molecular dynamics to evaluate these Kubo transforms. Benchmark calculations show that the approximations are useful for short times even for nonlinear operators, providing a consistent improvement over classical simulations of multi-time correlation functions. The introduced methodology thus provides a practical way of including nuclear quantum effects in multi-time response functions of non-linear optical spectroscopy.
Collapse
Affiliation(s)
- Kenneth A Jung
- Department of Chemistry, Yale University, P.O. Box 208107, New Haven, Connecticut 06520-8107, USA
| | - Pablo E Videla
- Department of Chemistry, Yale University, P.O. Box 208107, New Haven, Connecticut 06520-8107, USA
| | - Victor S Batista
- Department of Chemistry, Yale University, P.O. Box 208107, New Haven, Connecticut 06520-8107, USA
| |
Collapse
|
61
|
Willatt MJ, Ceriotti M, Althorpe SC. Approximating Matsubara dynamics using the planetary model: Tests on liquid water and ice. J Chem Phys 2018; 148:102336. [PMID: 29544307 DOI: 10.1063/1.5004808] [Citation(s) in RCA: 21] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Matsubara dynamics is the quantum-Boltzmann-conserving classical dynamics which remains when real-time coherences are taken out of the exact quantum Liouvillian [T. J. H. Hele et al., J. Chem. Phys. 142, 134103 (2015)]; because of a phase-term, it cannot be used as a practical method without further approximation. Recently, Smith et al. [J. Chem. Phys. 142, 244112 (2015)] developed a "planetary" model dynamics which conserves the Feynman-Kleinert (FK) approximation to the quantum-Boltzmann distribution. Here, we show that for moderately anharmonic potentials, the planetary dynamics gives a good approximation to Matsubara trajectories on the FK potential surface by decoupling the centroid trajectory from the locally harmonic Matsubara fluctuations, which reduce to a single phase-less fluctuation particle (the "planet"). We also show that the FK effective frequency can be approximated by a direct integral over these fluctuations, obviating the need to solve iterative equations. This modification, together with use of thermostatted ring-polymer molecular dynamics, allows us to test the planetary model on water (gas-phase, liquid, and ice) using the q-TIP4P/F potential surface. The "planetary" fluctuations give a poor approximation to the rotational/librational bands in the infrared spectrum, but a good approximation to the bend and stretch bands, where the fluctuation lineshape is found to be motionally narrowed by the vibrations of the centroid.
Collapse
Affiliation(s)
- Michael J Willatt
- Department of Chemistry, University of Cambridge, Lensfield Road, Cambridge CB2 1EW, United Kingdom
| | - Michele Ceriotti
- Laboratory of Computational Science and Modeling, Institute of Materials, École Polytechnique Fédérale de Lausanne, 1015 Lausanne, Switzerland
| | - Stuart C Althorpe
- Department of Chemistry, University of Cambridge, Lensfield Road, Cambridge CB2 1EW, United Kingdom
| |
Collapse
|
62
|
Gaiduk AP, Gustafson J, Gygi F, Galli G. First-Principles Simulations of Liquid Water Using a Dielectric-Dependent Hybrid Functional. J Phys Chem Lett 2018; 9:3068-3073. [PMID: 29768015 DOI: 10.1021/acs.jpclett.8b01017] [Citation(s) in RCA: 19] [Impact Index Per Article: 3.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/08/2023]
Abstract
We carried out first-principles simulations of liquid water under ambient conditions using a dielectric-dependent hybrid functional, where the fraction of exact exchange is set equal to the inverse of the high-frequency dielectric constant of the liquid. We found excellent agreement with experiment for the oxygen-oxygen partial correlation function at the experimental equilibrium density and 311 ± 3 K. Other structural and dynamical properties, such as the diffusion coefficient, molecular dipole moments, and vibrational spectra, are also in good agreement with experiment. Our results, together with previous findings on electronic properties of the liquid with the same functional, show that the dielectric-dependent hybrid functional accurately describes both the structural and electronic properties of liquid water.
Collapse
Affiliation(s)
- Alex P Gaiduk
- Institute for Molecular Engineering , The University of Chicago , Chicago , Illinois 60637 , United States
- Materials Science Division , Argonne National Laboratory , Argonne , Illinois 60439 , United States
| | - Jeffrey Gustafson
- Department of Chemistry , The University of Chicago , Chicago , Illinois 60637 , United States
| | - François Gygi
- Department of Computer Science , University of California , Davis , California 95616 , United States
| | - Giulia Galli
- Institute for Molecular Engineering , The University of Chicago , Chicago , Illinois 60637 , United States
- Materials Science Division , Argonne National Laboratory , Argonne , Illinois 60439 , United States
- Department of Chemistry , The University of Chicago , Chicago , Illinois 60637 , United States
| |
Collapse
|
63
|
Robertson C, Habershon S. Harmonic-phase path-integral approximation of thermal quantum correlation functions. J Chem Phys 2018; 148:102316. [PMID: 29544325 DOI: 10.1063/1.5002189] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
We present an approximation to the thermal symmetric form of the quantum time-correlation function in the standard position path-integral representation. By transforming to a sum-and-difference position representation and then Taylor-expanding the potential energy surface of the system to second order, the resulting expression provides a harmonic weighting function that approximately recovers the contribution of the phase to the time-correlation function. This method is readily implemented in a Monte Carlo sampling scheme and provides exact results for harmonic potentials (for both linear and non-linear operators) and near-quantitative results for anharmonic systems for low temperatures and times that are likely to be relevant to condensed phase experiments. This article focuses on one-dimensional examples to provide insights into convergence and sampling properties, and we also discuss how this approximation method may be extended to many-dimensional systems.
Collapse
Affiliation(s)
- Christopher Robertson
- Department of Chemistry and Centre for Scientific Computing, University Of Warwick, Coventry CV4 7AL, United Kingdom
| | - Scott Habershon
- Department of Chemistry and Centre for Scientific Computing, University Of Warwick, Coventry CV4 7AL, United Kingdom
| |
Collapse
|
64
|
Videla PE, Rossky PJ, Laria D. Isotope effects in aqueous solvation of simple halides. J Chem Phys 2018; 148:102306. [DOI: 10.1063/1.4986231] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/15/2022] Open
Affiliation(s)
- Pablo E. Videla
- Departamento de Química Inorgánica Analítica y Química-Física e INQUIMAE, Facultad de Ciencias Exactas y Naturales, Universidad de Buenos Aires, Ciudad Universitaria, Pabellón II, 1428 Buenos Aires, Argentina
| | - Peter J. Rossky
- Department of Chemistry, Rice University, Houston, Texas 77251-1892, USA
| | - D. Laria
- Departamento de Química Inorgánica Analítica y Química-Física e INQUIMAE, Facultad de Ciencias Exactas y Naturales, Universidad de Buenos Aires, Ciudad Universitaria, Pabellón II, 1428 Buenos Aires, Argentina
- Departamento de Física de la Materia Condensada, Comisión Nacional de Energía Atómica, Avenida Libertador 8250, 1429 Buenos Aires, Argentina
| |
Collapse
|
65
|
Rossi M, Kapil V, Ceriotti M. Fine tuning classical and quantum molecular dynamics using a generalized Langevin equation. J Chem Phys 2018; 148:102301. [DOI: 10.1063/1.4990536] [Citation(s) in RCA: 42] [Impact Index Per Article: 7.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Mariana Rossi
- Laboratory of Computational Science and Modelling, Institute of Materials, Ecole Polytechnique Fédérale de Lausanne, Lausanne, Switzerland
| | - Venkat Kapil
- Laboratory of Computational Science and Modelling, Institute of Materials, Ecole Polytechnique Fédérale de Lausanne, Lausanne, Switzerland
| | - Michele Ceriotti
- Laboratory of Computational Science and Modelling, Institute of Materials, Ecole Polytechnique Fédérale de Lausanne, Lausanne, Switzerland
| |
Collapse
|
66
|
Tao X, Shushkov P, Miller TF. Path-integral isomorphic Hamiltonian for including nuclear quantum effects in non-adiabatic dynamics. J Chem Phys 2018; 148:102327. [DOI: 10.1063/1.5005544] [Citation(s) in RCA: 36] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Xuecheng Tao
- Division of Chemistry and Chemical Engineering, California Institute of Technology, Pasadena, California 91125, USA
| | - Philip Shushkov
- Division of Chemistry and Chemical Engineering, California Institute of Technology, Pasadena, California 91125, USA
| | - Thomas F. Miller
- Division of Chemistry and Chemical Engineering, California Institute of Technology, Pasadena, California 91125, USA
| |
Collapse
|
67
|
Ojha D, Henao A, Kühne TD. Nuclear quantum effects on the vibrational dynamics of liquid water. J Chem Phys 2018; 148:102328. [PMID: 29544291 DOI: 10.1063/1.5005500] [Citation(s) in RCA: 24] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Abstract
Based on quantum-mechanical path-integral molecular dynamics simulations, the impact of nuclear quantum effects on the vibrational and hydrogen bond dynamics in liquid water is investigated. The instantaneous fluctuations in the frequencies of the O-H stretch modes are calculated using the wavelet method of time-series analysis, while the time scales of the vibrational spectral diffusion are determined from frequency-time correlation functions, joint probability distributions, and the slope of three-pulse photon echo. We find that the inclusion of nuclear quantum effects leads not only to a redshift of the vibrational frequency distribution by around 130 cm-1 but also to an acceleration of the vibrational dynamics by as much as 30%. In addition, quantum fluctuations also entail a significantly faster decay of correlation in the initial diffusive regime, which is in agreement with recent vibrational echo experiments.
Collapse
Affiliation(s)
- Deepak Ojha
- Dynamics of Condensed Matter and Center for Sustainable Systems Design, Department of Chemistry, University of Paderborn, Warburger Str. 100, D-33098 Paderborn, Germany
| | - Andrés Henao
- Dynamics of Condensed Matter and Center for Sustainable Systems Design, Department of Chemistry, University of Paderborn, Warburger Str. 100, D-33098 Paderborn, Germany
| | - Thomas D Kühne
- Dynamics of Condensed Matter and Center for Sustainable Systems Design, Department of Chemistry, University of Paderborn, Warburger Str. 100, D-33098 Paderborn, Germany
| |
Collapse
|
68
|
Feng Y, Wang Z, Guo J, Chen J, Wang EG, Jiang Y, Li XZ. The collective and quantum nature of proton transfer in the cyclic water tetramer on NaCl(001). J Chem Phys 2018; 148:102329. [PMID: 29544296 DOI: 10.1063/1.5004737] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Proton tunneling is an elementary process in the dynamics of hydrogen-bonded systems. Collective tunneling is known to exist for a long time. Atomistic investigations of this mechanism in realistic systems, however, are scarce. Using a combination of ab initio theoretical and high-resolution experimental methods, we investigate the role played by the protons on the chirality switching of a water tetramer on NaCl(001). Our scanning tunneling spectroscopies show that partial deuteration of the H2O tetramer with only one D2O leads to a significant suppression of the chirality switching rate at a cryogenic temperature (T), indicating that the chirality switches by tunneling in a concerted manner. Theoretical simulations, in the meantime, support this picture by presenting a much smaller free-energy barrier for the translational collective proton tunneling mode than other chirality switching modes at low T. During this analysis, the virial energy provides a reasonable estimator for the description of the nuclear quantum effects when a traditional thermodynamic integration method cannot be used, which could be employed in future studies of similar problems. Given the high-dimensional nature of realistic systems and the topology of the hydrogen-bonded network, collective proton tunneling may exist more ubiquitously than expected. Systems of this kind can serve as ideal platforms for studies of this mechanism, easily accessible to high-resolution experimental measurements.
Collapse
Affiliation(s)
- Yexin Feng
- International Center for Quantum Materials and School of Physics, Peking University, Beijing 100871, People's Republic of China
| | - Zhichang Wang
- International Center for Quantum Materials and School of Physics, Peking University, Beijing 100871, People's Republic of China
| | - Jing Guo
- International Center for Quantum Materials and School of Physics, Peking University, Beijing 100871, People's Republic of China
| | - Ji Chen
- International Center for Quantum Materials and School of Physics, Peking University, Beijing 100871, People's Republic of China
| | - En-Ge Wang
- International Center for Quantum Materials and School of Physics, Peking University, Beijing 100871, People's Republic of China
| | - Ying Jiang
- International Center for Quantum Materials and School of Physics, Peking University, Beijing 100871, People's Republic of China
| | - Xin-Zheng Li
- International Center for Quantum Materials and School of Physics, Peking University, Beijing 100871, People's Republic of China
| |
Collapse
|
69
|
|
70
|
Affiliation(s)
- Xinzijian Liu
- Beijing National Laboratory For Molecular Sciences, Institute of Theoretical and Computational Chemistry, College of Chemistry and Molecular Engineering, Peking University, Beijing, China
| | - Jian Liu
- Beijing National Laboratory For Molecular Sciences, Institute of Theoretical and Computational Chemistry, College of Chemistry and Molecular Engineering, Peking University, Beijing, China
| |
Collapse
|
71
|
Reddy SK, Moberg DR, Straight SC, Paesani F. Temperature-dependent vibrational spectra and structure of liquid water from classical and quantum simulations with the MB-pol potential energy function. J Chem Phys 2017; 147:244504. [DOI: 10.1063/1.5006480] [Citation(s) in RCA: 65] [Impact Index Per Article: 9.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/03/2023] Open
Affiliation(s)
- Sandeep K. Reddy
- Department of Chemistry and Biochemistry, University of California San Diego, La Jolla, California 92093, USA
| | - Daniel R. Moberg
- Department of Chemistry and Biochemistry, University of California San Diego, La Jolla, California 92093, USA
| | - Shelby C. Straight
- Department of Chemistry and Biochemistry, University of California San Diego, La Jolla, California 92093, USA
| | - Francesco Paesani
- Department of Chemistry and Biochemistry, Materials Science and Engineering, and San Diego Supercomputer Center, University of California San Diego, La Jolla, California 92093, USA
| |
Collapse
|
72
|
Buxton SJ, Habershon S. Accelerated path-integral simulations using ring-polymer interpolation. J Chem Phys 2017; 147:224107. [DOI: 10.1063/1.5006465] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/17/2022] Open
Affiliation(s)
- Samuel J. Buxton
- Department of Chemistry and Centre for Scientific Computing, University of Warwick, Coventry CV4 7AL, United Kingdom
| | - Scott Habershon
- Department of Chemistry and Centre for Scientific Computing, University of Warwick, Coventry CV4 7AL, United Kingdom
| |
Collapse
|
73
|
Folpini G, Siebert T, Woerner M, Abel S, Laage D, Elsaesser T. Water Librations in the Hydration Shell of Phospholipids. J Phys Chem Lett 2017; 8:4492-4497. [PMID: 28858510 DOI: 10.1021/acs.jpclett.7b01942] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/07/2023]
Abstract
The hydrophilic phosphate moiety in the headgroup of phospholipids forms strong hydrogen bonds with water molecules in the first hydration layer. Time-domain terahertz spectroscopy in a range from 100 to 1000 cm-1 reveals the influence of such interactions on rotations of water molecules. We determine librational absorption spectra of water nanopools in phospholipid reverse micelles for a range from w0 = 2 to 16 waters per phospholipid molecule. A pronounced absorption feature with maximum at 830 cm-1 is superimposed on a broad absorption band between 300 and 1000 cm-1. Molecular dynamics simulations of water in the reverse micelles suggest that the feature at 830 cm-1 arises from water molecules forming one or two strong hydrogen bonds with phosphate groups, while the broad component comes from bulk-like environments. This behavior is markedly different from water interacting with less polar surfaces.
Collapse
Affiliation(s)
- Giulia Folpini
- Max-Born-Institut für Nichtlineare Optik und Kurzzeitspektroskopie , D-12489 Berlin, Germany
| | - Torsten Siebert
- Max-Born-Institut für Nichtlineare Optik und Kurzzeitspektroskopie , D-12489 Berlin, Germany
| | - Michael Woerner
- Max-Born-Institut für Nichtlineare Optik und Kurzzeitspektroskopie , D-12489 Berlin, Germany
| | - Stephane Abel
- Institut de Biologie Intégrative de la Cellule (I2BC), Institut Frédéric Joliot, CEA, CNRS, Univ Paris-Sud Université Paris-Saclay, 91405 Gif-Sur-Yvette Cedex, France
| | - Damien Laage
- École Normale Supérieure, PSL Research University, UPMC Univ Paris 06, CNRS, Département de Chimie, PASTEUR, 24 rue Lhomond, 75005 Paris, France
- Sorbonne Universités, UPMC Univ Paris 06, ENS, CNRS, PASTEUR, 75005 Paris, France
| | - Thomas Elsaesser
- Max-Born-Institut für Nichtlineare Optik und Kurzzeitspektroskopie , D-12489 Berlin, Germany
| |
Collapse
|
74
|
Xu ZH, Meuwly M. Vibrational Spectroscopy and Proton Transfer Dynamics in Protonated Oxalate. J Phys Chem A 2017; 121:5389-5398. [PMID: 28597659 DOI: 10.1021/acs.jpca.7b02234] [Citation(s) in RCA: 15] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
The dynamics and infrared spectroscopic signatures of proton transfer in protonated oxalate (p-Oxa) are studied using classical and quantum dynamics. The intermolecular interactions are described by a force field suitable to follow proton transfer. This allows to carry out multiple extended classical molecular dynamics (MD) and ring polymer MD simulations from which the infrared spectrum is determined. Simulations at 600 K sample the quantum mechanical ground state probability distribution and best reproduce the experimentally observed maximum absorption wavelength and part of the line shape. Comparison with the experimentally measured spectrum provides an estimate for the barrier height for proton transfer which can not be determined directly from experiment. A barrier of 4.2 kcal/mol is found to best reproduce the position and width of the infrared absorption of the transferring proton in p-Oxa and also leads to an infrared (IR) spectrum in good agreement with experiment for the deuterated species d-Oxa. A novel means to capture the two resonance forms of oxalate depending on the localization of the excess proton on either CO moiety is found to yield improved results for the spectroscopy in the framework region between 1000 and 2000 cm-1.
Collapse
Affiliation(s)
- Zhen-Hao Xu
- Department of Chemistry, University of Basel , Klingelbergstrasse 80, CH-4056 Basel, Switzerland
| | - Markus Meuwly
- Department of Chemistry, University of Basel , Klingelbergstrasse 80, CH-4056 Basel, Switzerland
| |
Collapse
|
75
|
Marsalek O, Markland TE. Quantum Dynamics and Spectroscopy of Ab Initio Liquid Water: The Interplay of Nuclear and Electronic Quantum Effects. J Phys Chem Lett 2017; 8:1545-1551. [PMID: 28296422 DOI: 10.1021/acs.jpclett.7b00391] [Citation(s) in RCA: 137] [Impact Index Per Article: 19.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/06/2023]
Abstract
Understanding the reactivity and spectroscopy of aqueous solutions at the atomistic level is crucial for the elucidation and design of chemical processes. However, the simulation of these systems requires addressing the formidable challenges of treating the quantum nature of both the electrons and nuclei. Exploiting our recently developed methods that provide acceleration by up to 2 orders of magnitude, we combine path integral simulations with on-the-fly evaluation of the electronic structure at the hybrid density functional theory level to capture the interplay between nuclear quantum effects and the electronic surface. Here we show that this combination provides accurate structure and dynamics, including the full infrared and Raman spectra of liquid water. This allows us to demonstrate and explain the failings of lower-level density functionals for dynamics and vibrational spectroscopy when the nuclei are treated quantum mechanically. These insights thus provide a foundation for the reliable investigation of spectroscopy and reactivity in aqueous environments.
Collapse
Affiliation(s)
- Ondrej Marsalek
- Department of Chemistry, Stanford University , Stanford, California 94305, United States
| | - Thomas E Markland
- Department of Chemistry, Stanford University , Stanford, California 94305, United States
| |
Collapse
|
76
|
Hirshberg B, Sagiv L, Gerber RB. Approximate Quantum Dynamics using Ab Initio Classical Separable Potentials: Spectroscopic Applications. J Chem Theory Comput 2017; 13:982-991. [DOI: 10.1021/acs.jctc.6b01129] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Barak Hirshberg
- Institute
of Chemistry and the Fritz Haber Center for Molecular Dynamics, The Hebrew University, Jerusalem 9190401, Israel
| | - Lior Sagiv
- Institute
of Chemistry and the Fritz Haber Center for Molecular Dynamics, The Hebrew University, Jerusalem 9190401, Israel
| | - R. Benny Gerber
- Institute
of Chemistry and the Fritz Haber Center for Molecular Dynamics, The Hebrew University, Jerusalem 9190401, Israel
- Department
of Chemistry, University of California, Irvine, California 92697, United States
| |
Collapse
|
77
|
Agarwal A, Clementi C, Delle Site L. Path integral-GC-AdResS simulation of a large hydrophobic solute in water: a tool to investigate the interplay between local microscopic structures and quantum delocalization of atoms in space. Phys Chem Chem Phys 2017; 19:13030-13037. [DOI: 10.1039/c7cp01629h] [Citation(s) in RCA: 16] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/26/2023]
Abstract
We perform large scale quantum (path integral) molecular dynamics simulations of a C60 -like molecule in water.
Collapse
|
78
|
An analysis of nonadiabatic ring-polymer molecular dynamics and its application to vibronic spectra. Chem Phys 2017. [DOI: 10.1016/j.chemphys.2016.09.036] [Citation(s) in RCA: 31] [Impact Index Per Article: 4.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/19/2022]
|
79
|
Welsch R, Song K, Shi Q, Althorpe SC, Miller TF. Non-equilibrium dynamics from RPMD and CMD. J Chem Phys 2016; 145:204118. [DOI: 10.1063/1.4967958] [Citation(s) in RCA: 42] [Impact Index Per Article: 5.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/17/2023] Open
Affiliation(s)
- Ralph Welsch
- Division of Chemistry and Chemical Engineering, California Institute of Technology,1200 E. California Blvd., Pasadena, California 91125, USA
| | - Kai Song
- Beijing National Laboratory for Molecular Sciences, State Key Laboratory for Structural Chemistry of Unstable and Stable Species, Institute of Chemistry, Chinese Academy of Sciences, Zhongguancun, Beijing 100190, China
| | - Qiang Shi
- Beijing National Laboratory for Molecular Sciences, State Key Laboratory for Structural Chemistry of Unstable and Stable Species, Institute of Chemistry, Chinese Academy of Sciences, Zhongguancun, Beijing 100190, China
| | - Stuart C. Althorpe
- Department of Chemistry, University of Cambridge, Lensfield Road, Cambridge CB2 1EW, United Kingdom
| | - Thomas F. Miller
- Division of Chemistry and Chemical Engineering, California Institute of Technology,1200 E. California Blvd., Pasadena, California 91125, USA
| |
Collapse
|
80
|
Analysis of the Spectral Characteristics of Pure Moxa Stick Burning by Hyperspectral Imaging and Fourier Transform Infrared Spectroscopy. EVIDENCE-BASED COMPLEMENTARY AND ALTERNATIVE MEDICINE 2016; 2016:1057878. [PMID: 27721889 PMCID: PMC5045997 DOI: 10.1155/2016/1057878] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 03/21/2016] [Revised: 08/02/2016] [Accepted: 08/10/2016] [Indexed: 11/18/2022]
Abstract
The objective of this study was to investigate the spectra characteristics (SC) at wavelengths of 400~1000 nm and 2.5~15.5 μm of pure moxa stick (MS) during its 25-minute burning process using new spectral imaging techniques. Spectral images were collected for the burning pure MS at 5, 10, 15, 20, and 25 min using hyperspectral imaging (HSI) and Fourier transform infrared spectroscopy (FTIR) for the first time. The results showed that, at wavelengths of 400~1000 nm, the spectral range of the cross section of MS burning was 750~980 nm; the peak position was 860 nm. At wavelengths of 2.5~15.5 μm, the spectral range of the cross section of MS burning was 3.0~4.0 μm; the peak position was approximately 3.5 μm. The radiation spectra of MS burning include litter red and amount of infrared (but mainly near infrared) wavelengths. The temperature, blood perfusion, and oxygen saturation increase of Shenshu (BL23) after moxibustion radiation were observed too. According to mechanism of photobiological effects and moxibustion biological effects, it was inferred that moxibustion effects should be linked with moxibustion SC. This study provided new data and means for physical properties of moxibustion research.
Collapse
|
81
|
Litman YE, Videla PE, Rodriguez J, Laria D. Positional Isotope Exchange in HX·(H2O)n (X = F, I) Clusters at Low Temperatures. J Phys Chem A 2016; 120:7213-24. [DOI: 10.1021/acs.jpca.6b06681] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Yair E. Litman
- Departamento
de Química Inorgánica, Analítica y Química-Física
and INQUIMAE-CONICET, Facultad de Ciencias Exactas y Naturales, Universidad de Buenos Aires Ciudad Universitaria, Pabellón II, 1428 Buenos Aires, Argentina
| | - Pablo E. Videla
- Departamento
de Química Inorgánica, Analítica y Química-Física
and INQUIMAE-CONICET, Facultad de Ciencias Exactas y Naturales, Universidad de Buenos Aires Ciudad Universitaria, Pabellón II, 1428 Buenos Aires, Argentina
| | - Javier Rodriguez
- Departamento
de Física de la Materia Condensada, Comisión Nacional de Energía Atómica, Avenida Libertador 8250, 1429 Buenos Aires, Argentina
- ECyT, UNSAM, Martín
de Irigoyen 3100, 1650 San Martín, Pcia. de Buenos Aires, Argentina
| | - Daniel Laria
- Departamento
de Química Inorgánica, Analítica y Química-Física
and INQUIMAE-CONICET, Facultad de Ciencias Exactas y Naturales, Universidad de Buenos Aires Ciudad Universitaria, Pabellón II, 1428 Buenos Aires, Argentina
- Departamento
de Física de la Materia Condensada, Comisión Nacional de Energía Atómica, Avenida Libertador 8250, 1429 Buenos Aires, Argentina
| |
Collapse
|
82
|
Minoshima Y, Seki Y, Takayanagi T, Shiga M. Effects of temperature and isotopic substitution on electron attachment dynamics of guanine–cytosine base pair: Ring-polymer and classical molecular dynamics simulations. Chem Phys 2016. [DOI: 10.1016/j.chemphys.2016.02.019] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/22/2022]
|
83
|
Burnham CJ, English NJ. Study of clathrate hydrates via equilibrium molecular-dynamics simulation employing polarisable and non-polarisable, rigid and flexible water models. J Chem Phys 2016; 144:164503. [PMID: 27131553 DOI: 10.1063/1.4947039] [Citation(s) in RCA: 22] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Christian J. Burnham
- School of Chemical and Bioprocess Engineering, University College Dublin, Belfield, Dublin 4, Ireland
| | - Niall J. English
- School of Chemical and Bioprocess Engineering, University College Dublin, Belfield, Dublin 4, Ireland
| |
Collapse
|
84
|
Ceriotti M, Fang W, Kusalik PG, McKenzie RH, Michaelides A, Morales MA, Markland TE. Nuclear Quantum Effects in Water and Aqueous Systems: Experiment, Theory, and Current Challenges. Chem Rev 2016; 116:7529-50. [DOI: 10.1021/acs.chemrev.5b00674] [Citation(s) in RCA: 339] [Impact Index Per Article: 42.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/15/2022]
Affiliation(s)
- Michele Ceriotti
- Laboratory
of Computational Science and Modeling, Institute of Materials, École Polytechnique Fédérale de Lausanne, 1015 Lausanne, Switzerland
| | - Wei Fang
- Thomas
Young Centre, London Centre for Nanotechnology and Department of Physics
and Astronomy, University College London, London WC1E 6BT, United Kingdom
| | - Peter G. Kusalik
- Department
of Chemistry, University of Calgary, 2500 University Drive NW, Calgary, Alberta T2N 1N4, Canada
| | - Ross H. McKenzie
- School
of Mathematics and Physics, University of Queensland, Brisbane, 4072 Queensland Australia
| | - Angelos Michaelides
- Thomas
Young Centre, London Centre for Nanotechnology and Department of Physics
and Astronomy, University College London, London WC1E 6BT, United Kingdom
| | - Miguel A. Morales
- Lawrence Livermore National Laboratory, Livermore, California 94550, United States
| | - Thomas E. Markland
- Department
of Chemistry, Stanford University, 333 Campus Drive, Stanford, California 94305, United States
| |
Collapse
|
85
|
John C, Spura T, Habershon S, Kühne TD. Quantum ring-polymer contraction method: Including nuclear quantum effects at no additional computational cost in comparison to ab initio molecular dynamics. Phys Rev E 2016; 93:043305. [PMID: 27176426 DOI: 10.1103/physreve.93.043305] [Citation(s) in RCA: 27] [Impact Index Per Article: 3.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/12/2016] [Indexed: 06/05/2023]
Abstract
We present a simple and accurate computational method which facilitates ab initio path-integral molecular dynamics simulations, where the quantum-mechanical nature of the nuclei is explicitly taken into account, at essentially no additional computational cost in comparison to the corresponding calculation using classical nuclei. The predictive power of the proposed quantum ring-polymer contraction method is demonstrated by computing various static and dynamic properties of liquid water at ambient conditions using density functional theory. This development will enable routine inclusion of nuclear quantum effects in ab initio molecular dynamics simulations of condensed-phase systems.
Collapse
Affiliation(s)
- Christopher John
- Dynamics of Condensed Matter, Department of Chemistry, University of Paderborn, Warburger Strasse 100, D-33098 Paderborn, Germany
| | - Thomas Spura
- Dynamics of Condensed Matter, Department of Chemistry, University of Paderborn, Warburger Strasse 100, D-33098 Paderborn, Germany
| | - Scott Habershon
- Department of Chemistry and Centre for Scientific Computing, University of Warwick, Coventry CV4 7AL, United Kingdom
| | - Thomas D Kühne
- Dynamics of Condensed Matter, Department of Chemistry, University of Paderborn, Warburger Strasse 100, D-33098 Paderborn, Germany and Paderborn Center for Parallel Computing and Institute for Lightweight Design, Department of Chemistry, University of Paderborn, Warburger Strasse 100, D-33098 Paderborn, Germany
| |
Collapse
|
86
|
Medders GR, Paesani F. On the interplay of the potential energy and dipole moment surfaces in controlling the infrared activity of liquid water. J Chem Phys 2016; 142:212411. [PMID: 26049431 DOI: 10.1063/1.4916629] [Citation(s) in RCA: 30] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Infrared vibrational spectroscopy is a valuable tool for probing molecular structure and dynamics. However, obtaining an unambiguous molecular-level interpretation of the spectral features is made difficult, in part, due to the complex interplay of the dipole moment with the underlying vibrational structure. Here, we disentangle the contributions of the potential energy surface (PES) and dipole moment surface (DMS) to the infrared spectrum of liquid water by examining three classes of models, ranging in complexity from simple point charge models to accurate representations of the many-body interactions. By decoupling the PES from the DMS in the calculation of the infrared spectra, we demonstrate that the PES, by directly modulating the vibrational structure, primarily controls the width and position of the spectroscopic features. Due to the dependence of the molecular dipole moment on the hydration environment, many-body electrostatic effects result in a ∼100 cm(-1) redshift in the peak of the OH stretch band. Interestingly, while an accurate description of many-body collective motion is required to generate the correct (vibrational) structure of the liquid, the infrared intensity in the OH stretching region appears to be a measure of the local structure due to the dominance of the one-body and short-ranged two-body contributions to the total dipole moment.
Collapse
Affiliation(s)
- Gregory R Medders
- Department of Chemistry and Biochemistry, University of California, San Diego, La Jolla, California 92093, USA
| | - Francesco Paesani
- Department of Chemistry and Biochemistry, University of California, San Diego, La Jolla, California 92093, USA
| |
Collapse
|
87
|
Liu J, Zhang Z. Path integral Liouville dynamics: Applications to infrared spectra of OH, water, ammonia, and methane. J Chem Phys 2016; 144:034307. [DOI: 10.1063/1.4939953] [Citation(s) in RCA: 27] [Impact Index Per Article: 3.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Jian Liu
- Beijing National Laboratory for Molecular Sciences, Institute of Theoretical and Computational Chemistry, College of Chemistry and Molecular Engineering, Peking University, Beijing 100871, China
- State Key Joint Laboratory of Environmental Simulation and Pollution Control, College of Environmental Sciences and Engineering, Peking University, Beijing 100871, China
| | - Zhijun Zhang
- Beijing National Laboratory for Molecular Sciences, Institute of Theoretical and Computational Chemistry, College of Chemistry and Molecular Engineering, Peking University, Beijing 100871, China
| |
Collapse
|
88
|
Watanabe HC, Banno M, Sakurai M. An adaptive quantum mechanics/molecular mechanics method for the infrared spectrum of water: incorporation of the quantum effect between solute and solvent. Phys Chem Chem Phys 2016; 18:7318-33. [DOI: 10.1039/c5cp07136d] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Quantum effects in solute–solvent interactions, such as the many-body effect and the dipole-induced dipole, are known to be critical factors influencing the infrared spectra of species in the liquid phase.
Collapse
Affiliation(s)
- Hiroshi C. Watanabe
- Research Center for Advanced Science and Technology
- The University of Tokyo
- Meguro-ku
- Japan
| | - Misa Banno
- Center for Biological Resources and Informatics
- Tokyo Institute of Technology
- Midori-ku
- Japan
| | - Minoru Sakurai
- Center for Biological Resources and Informatics
- Tokyo Institute of Technology
- Midori-ku
- Japan
| |
Collapse
|
89
|
Baldauf C, Rossi M. Going clean: structure and dynamics of peptides in the gas phase and paths to solvation. JOURNAL OF PHYSICS. CONDENSED MATTER : AN INSTITUTE OF PHYSICS JOURNAL 2015; 27:493002. [PMID: 26598600 DOI: 10.1088/0953-8984/27/49/493002] [Citation(s) in RCA: 22] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/05/2023]
Abstract
The gas phase is an artificial environment for biomolecules that has gained much attention both experimentally and theoretically due to its unique characteristic of providing a clean room environment for the comparison between theory and experiment. In this review we give an overview mainly on first-principles simulations of isolated peptides and the initial steps of their interactions with ions and solvent molecules: a bottom up approach to the complexity of biological environments. We focus on the accuracy of different methods to explore the conformational space, the connections between theory and experiment regarding collision cross section evaluations and (anharmonic) vibrational spectra, and the challenges faced in this field.
Collapse
Affiliation(s)
- Carsten Baldauf
- Fritz Haber Institute, Faradayweg 4-6, 14195 Berlin, Germany
| | | |
Collapse
|
90
|
Wan Q, Galli G. First-Principles Framework to Compute Sum-Frequency Generation Vibrational Spectra of Semiconductors and Insulators. PHYSICAL REVIEW LETTERS 2015; 115:246404. [PMID: 26705645 DOI: 10.1103/physrevlett.115.246404] [Citation(s) in RCA: 20] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 06/19/2015] [Indexed: 05/22/2023]
Abstract
We present a first-principles framework to compute sum-frequency generation (SFG) vibrational spectra of semiconductors and insulators. The method is based on density functional theory and the use of maximally localized Wannier functions to compute the response to electric fields, and it includes the effect of electric field gradients at surfaces. In addition, it includes quadrupole contributions to SFG spectra, thus enabling the verification of the dipole approximation, whose validity determines the surface specificity of SFG spectroscopy. We compute the SFG spectra of ice I_{h} basal surfaces and identify which spectra components are affected by bulk contributions. Our results are in good agreement with experiments at low temperature.
Collapse
Affiliation(s)
- Quan Wan
- Institute for Molecular Engineering, The University of Chicago, Chicago, Illinois 60637, USA
| | - Giulia Galli
- Institute for Molecular Engineering, The University of Chicago, Chicago, Illinois 60637, USA
| |
Collapse
|
91
|
Medders GR, Götz AW, Morales MA, Bajaj P, Paesani F. On the representation of many-body interactions in water. J Chem Phys 2015; 143:104102. [DOI: 10.1063/1.4930194] [Citation(s) in RCA: 78] [Impact Index Per Article: 8.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Gregory R. Medders
- Department of Chemistry and Biochemistry, University of California, San Diego, La Jolla, California 92093, USA
| | - Andreas W. Götz
- San Diego Supercomputer Center, University of California, San Diego, La Jolla, California 92093, USA
| | - Miguel A. Morales
- Lawrence Livermore National Laboratory, 7000 East Avenue, Livermore, California 94550, USA
| | - Pushp Bajaj
- Department of Chemistry and Biochemistry, University of California, San Diego, La Jolla, California 92093, USA
| | - Francesco Paesani
- Department of Chemistry and Biochemistry, University of California, San Diego, La Jolla, California 92093, USA
| |
Collapse
|
92
|
Habershon S. Sampling reactive pathways with random walks in chemical space: Applications to molecular dissociation and catalysis. J Chem Phys 2015; 143:094106. [DOI: 10.1063/1.4929992] [Citation(s) in RCA: 59] [Impact Index Per Article: 6.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
|
93
|
Videla PE, Rossky PJ, Laria D. Isotopic Preferential Solvation of I– in Low-Temperature Water Nanoclusters. J Phys Chem B 2015; 119:11783-90. [DOI: 10.1021/acs.jpcb.5b05561] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Pablo E. Videla
- Departamento
de Química Inorgánica Analítica y Química-Física
e INQUIMAE, Facultad de Ciencias Exactas y Naturales, Universidad de Buenos Aires, Ciudad Universitaria, Pabellón II, 1428 Buenos Aires, Argentina
| | - Peter J. Rossky
- Department
of Chemistry, Rice University, Houston, Texas 77005-1892, United States
| | - Daniel Laria
- Departamento
de Química Inorgánica Analítica y Química-Física
e INQUIMAE, Facultad de Ciencias Exactas y Naturales, Universidad de Buenos Aires, Ciudad Universitaria, Pabellón II, 1428 Buenos Aires, Argentina
- Departamento
de Física de la Materia Condensada, Comisión Nacional de Energía Atómica, Avenida Libertador 8250, 1429 Buenos Aires, Argentina
| |
Collapse
|
94
|
Kessler J, Elgabarty H, Spura T, Karhan K, Partovi-Azar P, Hassanali AA, Kühne TD. Structure and Dynamics of the Instantaneous Water/Vapor Interface Revisited by Path-Integral and Ab Initio Molecular Dynamics Simulations. J Phys Chem B 2015; 119:10079-86. [DOI: 10.1021/acs.jpcb.5b04185] [Citation(s) in RCA: 51] [Impact Index Per Article: 5.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Jan Kessler
- Institute
of Physical Chemistry and Center of Computational Sciences, Johannes Gutenberg University Mainz, Staudinger Weg 9, D-55128 Mainz, Germany
| | - Hossam Elgabarty
- Institute
of Physical Chemistry and Center of Computational Sciences, Johannes Gutenberg University Mainz, Staudinger Weg 9, D-55128 Mainz, Germany
| | - Thomas Spura
- Dynamics
of Condensed Matter, Department of Chemistry, University of Paderborn, Warburger Str. 100, D-33098 Paderborn, Germany
| | - Kristof Karhan
- Dynamics
of Condensed Matter, Department of Chemistry, University of Paderborn, Warburger Str. 100, D-33098 Paderborn, Germany
| | - Pouya Partovi-Azar
- Dynamics
of Condensed Matter, Department of Chemistry, University of Paderborn, Warburger Str. 100, D-33098 Paderborn, Germany
| | - Ali A. Hassanali
- The Abdus Salam International Centre for Theoretical Physics, Strada Costiera 11, I-34151 Trieste, Italy
| | - Thomas D. Kühne
- Dynamics
of Condensed Matter, Department of Chemistry, University of Paderborn, Warburger Str. 100, D-33098 Paderborn, Germany
- Paderborn Center for Parallel Computing and Institute for Lightweight Design with Hybrid Systems, Warburger Str. 100, D-33098 Paderborn, Germany
| |
Collapse
|
95
|
Liu H, Wang Y, Bowman JM. Quantum calculations of the IR spectrum of liquid water using ab initio and model potential and dipole moment surfaces and comparison with experiment. J Chem Phys 2015; 142:194502. [PMID: 26001464 DOI: 10.1063/1.4921045] [Citation(s) in RCA: 42] [Impact Index Per Article: 4.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
The calculation and characterization of the IR spectrum of liquid water have remained a challenge for theory. In this paper, we address this challenge using a combination of ab initio approaches, namely, a quantum treatment of IR spectrum using the ab initio WHBB water potential energy surface and a refined ab initio dipole moment surface. The quantum treatment is based on the embedded local monomer method, in which the three intramolecular modes of each embedded H2O monomer are fully coupled and also coupled singly to each of six intermolecular modes. The new dipole moment surface consists of a previous spectroscopically accurate 1-body dipole moment surface and a newly fitted ab initio intrinsic 2-body dipole moment. A detailed analysis of the new dipole moment surface in terms of the coordinate dependence of the effective atomic charges is done along with tests of it for the water dimer and prism hexamer double-harmonic spectra against direct ab initio calculations. The liquid configurations are taken from previous molecular dynamics calculations of Skinner and co-workers, using the TIP4P plus E3B rigid monomer water potential. The IR spectrum of water at 300 K in the range of 0-4000 cm(-1) is calculated and compared with experiment, using the ab initio WHBB potential and new ab initio dipole moment, the q-TIP4P/F potential, which has a fixed-charged description of the dipole moment, and the TTM3-F potential and dipole moment surfaces. The newly calculated ab initio spectrum is in very good agreement with experiment throughout the above spectral range, both in band positions and intensities. This contrasts to results with the other potentials and dipole moments, especially the fixed-charge q-TIP4P/F model, which gives unrealistic intensities. The calculated ab initio spectrum is analyzed by examining the contribution of various transitions to each band.
Collapse
Affiliation(s)
- Hanchao Liu
- Cherry L. Emerson Center for Scientific Computation and Department of Chemistry, Emory University, Atlanta, Georgia 30322, USA
| | - Yimin Wang
- Cherry L. Emerson Center for Scientific Computation and Department of Chemistry, Emory University, Atlanta, Georgia 30322, USA
| | - Joel M Bowman
- Cherry L. Emerson Center for Scientific Computation and Department of Chemistry, Emory University, Atlanta, Georgia 30322, USA
| |
Collapse
|
96
|
Ni Y, Skinner JL. IR and SFG vibrational spectroscopy of the water bend in the bulk liquid and at the liquid-vapor interface, respectively. J Chem Phys 2015; 143:014502. [DOI: 10.1063/1.4923462] [Citation(s) in RCA: 69] [Impact Index Per Article: 7.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Yicun Ni
- Theoretical Chemistry Institute and Department of Chemistry, University of Wisconsin, Madison, Wisconsin 53706, USA
| | - J. L. Skinner
- Theoretical Chemistry Institute and Department of Chemistry, University of Wisconsin, Madison, Wisconsin 53706, USA
| |
Collapse
|
97
|
Smith KKG, Poulsen JA, Nyman G, Cunsolo A, Rossky PJ. Application of a new ensemble conserving quantum dynamics simulation algorithm to liquid para-hydrogen and ortho-deuterium. J Chem Phys 2015; 142:244113. [DOI: 10.1063/1.4922888] [Citation(s) in RCA: 21] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Kyle K. G. Smith
- Institute for Computational Engineering and Sciences and Department of Chemistry,University of Texas at Austin, Austin, Texas 78712, USA
| | - Jens Aage Poulsen
- Physical Chemistry, Department of Chemistry and Molecular Biology, University of Gothenburg, SE 41296 Gothenburg, Sweden
| | - Gunnar Nyman
- Physical Chemistry, Department of Chemistry and Molecular Biology, University of Gothenburg, SE 41296 Gothenburg, Sweden
| | - Alessandro Cunsolo
- Photon Sciences Directorate, Brookhaven National Laboratory, P.O. Box 5000, Upton, New York 11973, USA
| | - Peter J. Rossky
- Department of Chemistry, Rice University, Houston, Texas 77251, USA
| |
Collapse
|
98
|
Hele TJH, Willatt MJ, Muolo A, Althorpe SC. Boltzmann-conserving classical dynamics in quantum time-correlation functions: “Matsubara dynamics”. J Chem Phys 2015; 142:134103. [DOI: 10.1063/1.4916311] [Citation(s) in RCA: 71] [Impact Index Per Article: 7.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Timothy J. H. Hele
- Department of Chemistry, University of Cambridge, Lensfield Road, Cambridge CB2 1EW, United Kingdom
| | - Michael J. Willatt
- Department of Chemistry, University of Cambridge, Lensfield Road, Cambridge CB2 1EW, United Kingdom
| | - Andrea Muolo
- Department of Chemistry, University of Cambridge, Lensfield Road, Cambridge CB2 1EW, United Kingdom
| | - Stuart C. Althorpe
- Department of Chemistry, University of Cambridge, Lensfield Road, Cambridge CB2 1EW, United Kingdom
| |
Collapse
|
99
|
Honda T, Minoshima Y, Yokoi Y, Takayanagi T, Shiga M. Semiclassical dynamics of electron attachment to guanine–cytosine base pair. Chem Phys Lett 2015. [DOI: 10.1016/j.cplett.2015.03.006] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/25/2022]
|
100
|
Medders GR, Paesani F. Infrared and Raman Spectroscopy of Liquid Water through "First-Principles" Many-Body Molecular Dynamics. J Chem Theory Comput 2015; 11:1145-54. [PMID: 26579763 DOI: 10.1021/ct501131j] [Citation(s) in RCA: 161] [Impact Index Per Article: 17.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/07/2023]
Abstract
Vibrational spectroscopy is a powerful technique to probe the structure and dynamics of water. However, deriving an unambiguous molecular-level interpretation of the experimental spectral features remains a challenge due to the complexity of the underlying hydrogen-bonding network. In this contribution, we present an integrated theoretical and computational framework (named many-body molecular dynamics or MB-MD) that, by systematically removing uncertainties associated with existing approaches, enables a rigorous modeling of vibrational spectra of water from quantum dynamical simulations. Specifically, we extend approaches used to model the many-body expansion of interaction energies to develop many-body representations of the dipole moment and polarizability of water. The combination of these "first-principles" representations with centroid molecular dynamics simulations enables the simulation of infrared and Raman spectra of liquid water under ambient conditions that, without relying on any ad hoc parameters, are in good agreement with the corresponding experimental results. Importantly, since the many-body energy, dipole, and polarizability surfaces employed in the simulations are derived independently from accurate fits to correlated electronic structure data, MB-MD allows for a systematic analysis of the calculated spectra in terms of both electronic and dynamical contributions. The present analysis suggests that, while MB-MD correctly reproduces both the shifts and the shapes of the main spectroscopic features, an improved description of quantum dynamical effects possibly combined with a dissociable water potential may be necessary for a quantitative representation of the OH stretch band.
Collapse
Affiliation(s)
- Gregory R Medders
- Department of Chemistry and Biochemistry, University of California, San Diego , La Jolla, California 92037, United States
| | - Francesco Paesani
- Department of Chemistry and Biochemistry, University of California, San Diego , La Jolla, California 92037, United States
| |
Collapse
|