1
|
Lawrence JE, Lieberherr AZ, Fletcher T, Manolopoulos DE. Fast Quasi-Centroid Molecular Dynamics for Water and Ice. J Phys Chem B 2023; 127:9172-9180. [PMID: 37830934 PMCID: PMC10614180 DOI: 10.1021/acs.jpcb.3c05028] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/26/2023] [Revised: 09/21/2023] [Indexed: 10/14/2023]
Abstract
We describe how the fast quasi-centroid molecular dynamics (f-QCMD) method can be applied to condensed-phase systems by approximating the quasi-centroid potential of mean force as a sum of inter- and intramolecular corrections to the classical interaction potential. The corrections are found by using a regularized iterative Boltzmann inversion procedure to recover the inter- and intramolecular quasi-centroid distribution functions obtained from a path integral molecular dynamics simulation. The resulting methodology is found to give good agreement with a previously published QCMD dipole absorption spectrum for liquid water and satisfactory agreement for ice. It also gives good agreement with spectra from a recent implementation of CMD that uses a precomputed elevated temperature potential of mean force. Modern centroid molecular dynamics methods, therefore, appear to be reaching a consensus regarding the impact of nuclear quantum effects on the vibrational spectra of water and ice.
Collapse
Affiliation(s)
| | - Annina Z. Lieberherr
- Physical
and Theoretical Chemistry Laboratory, Department of Chemistry, University of Oxford, South Parks Road, Oxford OX1 3QZ, United
Kingdom
| | - Theo Fletcher
- Physical
and Theoretical Chemistry Laboratory, Department of Chemistry, University of Oxford, South Parks Road, Oxford OX1 3QZ, United
Kingdom
| | - David E. Manolopoulos
- Physical
and Theoretical Chemistry Laboratory, Department of Chemistry, University of Oxford, South Parks Road, Oxford OX1 3QZ, United
Kingdom
| |
Collapse
|
2
|
Litman Y, Lan J, Nagata Y, Wilkins DM. Fully First-Principles Surface Spectroscopy with Machine Learning. J Phys Chem Lett 2023; 14:8175-8182. [PMID: 37671886 PMCID: PMC10510433 DOI: 10.1021/acs.jpclett.3c01989] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/18/2023] [Accepted: 08/29/2023] [Indexed: 09/07/2023]
Abstract
Our current understanding of the structure and dynamics of aqueous interfaces at the molecular level has grown substantially due to the continuous development of surface-specific spectroscopies, such as vibrational sum-frequency generation (VSFG). As in other vibrational spectroscopies, we must turn to atomistic simulations to extract all of the information encoded in the VSFG spectra. The high computational cost associated with existing methods means that they have limitations in representing systems with complex electronic structure or in achieving statistical convergence. In this work, we combine high-dimensional neural network interatomic potentials and symmetry-adapted Gaussian process regression to overcome these constraints. We show that it is possible to model VSFG signals with fully ab initio accuracy using machine learning and illustrate the versatility of our approach on the water/air interface. Our strategy allows us to identify the main sources of theoretical inaccuracy and establish a clear pathway toward the modeling of surface-sensitive spectroscopy of complex interfaces.
Collapse
Affiliation(s)
- Yair Litman
- Yusuf
Hamied Department of Chemistry, University
of Cambridge, Lensfield Road, Cambridge CB2 1EW, U.K.
- Max
Planck Institute for Polymer Research, Ackermannweg 10, 55128 Mainz, Germany
| | - Jinggang Lan
- Department
of Chemistry, New York University, New York, New York 10003, United States
- Simons
Center for Computational Physical Chemistry at New York University, New York, New York 10003, United States
| | - Yuki Nagata
- Max
Planck Institute for Polymer Research, Ackermannweg 10, 55128 Mainz, Germany
| | - David M. Wilkins
- Centre
for Quantum Materials and Technologies School of Mathematics and Physics, Queen’s University Belfast, Belfast BT7 1NN, Northern Ireland, United Kingdom
| |
Collapse
|
3
|
Prada A, Pós ES, Althorpe SC. Comparison of Matsubara dynamics with exact quantum dynamics for an oscillator coupled to a dissipative bath. J Chem Phys 2023; 158:114106. [PMID: 36948794 DOI: 10.1063/5.0138250] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 03/17/2023] Open
Abstract
We report the first numerical calculations in which converged Matsubara dynamics is compared directly with exact quantum dynamics with no artificial damping of the time-correlation functions (TCFs). The system treated is a Morse oscillator coupled to a harmonic bath. We show that, when the system-bath coupling is sufficiently strong, the Matsubara calculations can be converged by explicitly including up to M = 200 Matsubara modes, with the remaining modes included as a harmonic "tail" correction. The resulting Matsubara TCFs are in near-perfect agreement with the exact quantum TCFs, for non-linear as well as linear operators, at a temperature at which the TCFs are dominated by quantum thermal fluctuations. These results provide compelling evidence that incoherent classical dynamics can arise in the condensed phase at temperatures at which the statistics are dominated by quantum (Boltzmann) effects, as a result of smoothing of imaginary-time Feynman paths. The techniques developed here may also lead to efficient methods for benchmarking system-bath dynamics in the overdamped regime.
Collapse
Affiliation(s)
- Adam Prada
- Yusuf Hamied Department of Chemistry, University of Cambridge, Lensfield Road, Cambridge CB2 1EW, United Kingdom
| | - Eszter S Pós
- Yusuf Hamied Department of Chemistry, University of Cambridge, Lensfield Road, Cambridge CB2 1EW, United Kingdom
| | - Stuart C Althorpe
- Yusuf Hamied Department of Chemistry, University of Cambridge, Lensfield Road, Cambridge CB2 1EW, United Kingdom
| |
Collapse
|
4
|
Tan JA, Takahashi K. Vibrational spectrum of a 1D oscillator: The quantum, the Wigner, and the classical ways. J CHIN CHEM SOC-TAIP 2022. [DOI: 10.1002/jccs.202200398] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/07/2022]
Affiliation(s)
- Jake A. Tan
- Institute of Atomic and Molecular Sciences Academia Sinica Taipei Taiwan
| | - Kaito Takahashi
- Institute of Atomic and Molecular Sciences Academia Sinica Taipei Taiwan
| |
Collapse
|
5
|
Jung KA, Markland TE. 2D spectroscopies from condensed phase dynamics: Accessing third-order response properties from equilibrium multi-time correlation functions. J Chem Phys 2022; 157:094111. [DOI: 10.1063/5.0107087] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
The third-order response lies at the heart of simulating and interpreting nonlinear spectroscopies ranging from two dimensional infrared (2D-IR) to 2D electronic (2D-ES), and 2D sum frequency generation (2D-SFG). The extra time and frequency dimensions in these spectroscopies provides access to rich information on the electronic and vibrational states present, the coupling between them, and the resulting rates at which they exchange energy that are obscured in linear spectroscopy, particularly for condensed phase systems that usually contain many overlapping features. While the exact quantum expression for the third-order response is well established it is incompatible with the methods that are practical for calculating the atomistic dynamics of large condensed phase systems. These methods, which include both classical mechanics and quantum dynamics methods that retain quantum statistical properties while obeying the symmetries of classical dynamics, such as LSC-IVR, Centroid Molecular Dynamics (CMD) and Ring Polymer Molecular Dynamics (RPMD) naturally provide short-time approximations to the multi-time symmetrized Kubo transformed correlation function. Here, we show how the third-order response can be formulated in terms of equilibrium symmetrized Kubo transformed correlation functions. We demonstrate the utility and accuracy of our approach by showing how it can be used to obtain the third-order response of a series of model systems using both classical dynamics and RPMD. In particular, we show that this approach captures features such as anharmonically induced vertical splittings and peak shifts while providing a physically transparent framework for understanding multidimensional spectroscopies.
Collapse
|
6
|
Simulation of Nuclear Quantum Effects in Condensed Matter Systems via Quantum Baths. APPLIED SCIENCES-BASEL 2022. [DOI: 10.3390/app12094756] [Citation(s) in RCA: 5] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/10/2022]
Abstract
This paper reviews methods that aim at simulating nuclear quantum effects (NQEs) using generalized thermal baths. Generalized (or quantum) baths simulate statistical quantum features, and in particular zero-point energy effects, through non-Markovian stochastic dynamics. They make use of generalized Langevin Equations (GLEs), in which the quantum Bose–Einstein energy distribution is enforced by tuning the random and friction forces, while the system degrees of freedom remain classical. Although these baths have been formally justified only for harmonic oscillators, they perform well for several systems, while keeping the cost of the simulations comparable to the classical ones. We review the formal properties and main characteristics of classical and quantum GLEs, in relation with the fluctuation–dissipation theorems. Then, we describe the quantum thermostat and quantum thermal bath, the two generalized baths currently most used, providing several examples of applications for condensed matter systems, including the calculation of vibrational spectra. The most important drawback of these methods, zero-point energy leakage, is discussed in detail with the help of model systems, and a recently proposed scheme to monitor and mitigate or eliminate it—the adaptive quantum thermal bath—is summarised. This approach considerably extends the domain of application of generalized baths, leading, for instance, to the successful simulation of liquid water, where a subtle interplay of NQEs is at play. The paper concludes by overviewing further development opportunities and open challenges of generalized baths.
Collapse
|
7
|
Begušić T, Tao X, Blake GA, Miller TF. Equilibrium-nonequilibrium ring-polymer molecular dynamics for nonlinear spectroscopy. J Chem Phys 2022; 156:131102. [PMID: 35395895 DOI: 10.1063/5.0087156] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/14/2022] Open
Abstract
Two-dimensional Raman and hybrid terahertz-Raman spectroscopic techniques provide invaluable insight into molecular structures and dynamics of condensed-phase systems. However, corroborating experimental results with theory is difficult due to the high computational cost of incorporating quantum-mechanical effects in the simulations. Here, we present the equilibrium-nonequilibrium ring-polymer molecular dynamics (RPMD), a practical computational method that can account for nuclear quantum effects on the two-time response function of nonlinear optical spectroscopy. Unlike a recently developed approach based on the double Kubo transformed (DKT) correlation function, our method is exact in the classical limit, where it reduces to the established equilibrium-nonequilibrium classical molecular dynamics method. Using benchmark model calculations, we demonstrate the advantages of the equilibrium-nonequilibrium RPMD over classical and DKT-based approaches. Importantly, its derivation, which is based on the nonequilibrium RPMD, obviates the need for identifying an appropriate Kubo transformed correlation function and paves the way for applying real-time path-integral techniques to multidimensional spectroscopy.
Collapse
Affiliation(s)
- Tomislav Begušić
- Division of Chemistry and Chemical Engineering, California Institute of Technology, Pasadena, California 91125, USA
| | - Xuecheng Tao
- Division of Chemistry and Chemical Engineering, California Institute of Technology, Pasadena, California 91125, USA
| | - Geoffrey A Blake
- 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
|
8
|
Fletcher T, Zhu A, Lawrence JE, Manolopoulos DE. Fast quasi-centroid molecular dynamics. J Chem Phys 2021; 155:231101. [PMID: 34937347 DOI: 10.1063/5.0076704] [Citation(s) in RCA: 7] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
We describe a fast implementation of the quasi-centroid molecular dynamics (QCMD) method in which the quasi-centroid potential of mean force is approximated as a separable correction to the classical interaction potential. This correction is obtained by first calculating quasi-centroid radial and angular distribution functions in a short path integral molecular dynamics simulation and then using iterative Boltzmann inversion to obtain an effective classical potential that reproduces these distribution functions in a classical NVT simulation. We illustrate this approach with example applications to the vibrational spectra of gas phase molecules, obtaining excellent agreement with QCMD reference calculations for water and ammonia and good agreement with the quantum mechanical vibrational spectrum of methane.
Collapse
Affiliation(s)
- Theo Fletcher
- Physical and Theoretical Chemistry Laboratory, Department of Chemistry, University of Oxford, South Parks Road, Oxford OX1 3QZ, United Kingdom
| | - Andrew Zhu
- Physical and Theoretical Chemistry Laboratory, Department of Chemistry, University of Oxford, South Parks Road, Oxford OX1 3QZ, United Kingdom
| | - Joseph E Lawrence
- Laboratory of Physical Chemistry, ETH Zürich, 8093 Zürich, Switzerland
| | - David E Manolopoulos
- Physical and Theoretical Chemistry Laboratory, Department of Chemistry, University of Oxford, South Parks Road, Oxford OX1 3QZ, United Kingdom
| |
Collapse
|
9
|
Haggard C, Sadhasivam VG, Trenins G, Althorpe SC. Testing the quasicentroid molecular dynamics method on gas-phase ammonia. J Chem Phys 2021; 155:174120. [PMID: 34742190 DOI: 10.1063/5.0068250] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Quasicentroid molecular dynamics (QCMD) is a path-integral method for approximating nuclear quantum effects in dynamics simulations, which has given promising results for gas- and condensed-phase water. In this work, by simulating the infrared spectrum of gas-phase ammonia, we test the feasibility of extending QCMD beyond water. Overall, QCMD works as well for ammonia as for water, reducing or eliminating blue shifts from the classical spectrum without introducing the artificial red shifts or broadening associated with other imaginary-time path-integral methods. However, QCMD gives only a modest improvement over the classical spectrum for the position of the symmetric bend mode, which is highly anharmonic (since it correlates with the inversion pathway). We expect QCMD to have similar problems with large-amplitude degrees of freedom in other molecules but otherwise to work as well as for water.
Collapse
Affiliation(s)
- Christopher Haggard
- Yusuf Hamied Department of Chemistry, University of Cambridge, Lensfield Road, Cambridge CB2 1EW, United Kingdom
| | - Vijay Ganesh Sadhasivam
- Yusuf Hamied Department of Chemistry, University of Cambridge, Lensfield Road, Cambridge CB2 1EW, United Kingdom
| | - George Trenins
- Laboratory of Physical Chemistry, ETH Zürich, 8093 Zürich, Switzerland
| | - Stuart C Althorpe
- Yusuf Hamied Department of Chemistry, University of Cambridge, Lensfield Road, Cambridge CB2 1EW, United Kingdom
| |
Collapse
|
10
|
Plé T, Huppert S, Finocchi F, Depondt P, Bonella S. Anharmonic spectral features via trajectory-based quantum dynamics: A perturbative analysis of the interplay between dynamics and sampling. J Chem Phys 2021; 155:104108. [PMID: 34525824 DOI: 10.1063/5.0056824] [Citation(s) in RCA: 15] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
The performance of different approximate algorithms for computing anharmonic features in vibrational spectra is analyzed and compared on model and more realistic systems that present relevant nuclear quantum effects. The methods considered combine approximate sampling of the quantum thermal distribution with classical time propagation and include Matsubara dynamics, path integral dynamics approaches, linearized initial value representation, and the recently introduced adaptive quantum thermal bath. A perturbative analysis of these different methods enables us to account for the observed numerical performance on prototypes for overtones and combination bands and to draw qualitatively correct trends for the numerical results obtained for Fermi resonances. Our results prove that the unequal performances of these approaches often derive from the method employed to sample initial conditions and not, as usually assumed, from the lack of coherence in the time propagation. Furthermore, as confirmed by the analysis reported in Benson and Althorpe, J. Chem. Phys. 130, 194510 (2021), we demonstrate, both via the perturbative approach and numerically, that path integral dynamics methods fail to reproduce the intensities of these anharmonic features and follow purely classical trends with respect to their temperature behavior. Finally, the remarkably accurate performance of the adaptive quantum thermal bath approach is documented and motivated.
Collapse
Affiliation(s)
- Thomas Plé
- Sorbonne Université, CNRS, Institut des NanoSciences de Paris, INSP, 4 Place Jussieu, F-75005 Paris, France
| | - Simon Huppert
- Sorbonne Université, CNRS, Institut des NanoSciences de Paris, INSP, 4 Place Jussieu, F-75005 Paris, France
| | - Fabio Finocchi
- Sorbonne Université, CNRS, Institut des NanoSciences de Paris, INSP, 4 Place Jussieu, F-75005 Paris, France
| | - Philippe Depondt
- Sorbonne Université, CNRS, Institut des NanoSciences de Paris, INSP, 4 Place Jussieu, F-75005 Paris, France
| | - Sara Bonella
- CECAM Centre Européen de Calcul Atomique et Moléculaire, École Polytechnique Fédérale de Lausanne, Batochimie, Avenue Forel 2, 1015 Lausanne, Switzerland
| |
Collapse
|