1
|
Drużbicki K, Gila-Herranz P, Marin-Villa P, Gaboardi M, Armstrong J, Fernandez-Alonso F. Cation Dynamics as Structure Explorer in Hybrid Perovskites-The Case of MAPbI 3. CRYSTAL GROWTH & DESIGN 2024; 24:391-404. [PMID: 38188269 PMCID: PMC10768891 DOI: 10.1021/acs.cgd.3c01112] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 09/19/2023] [Revised: 11/20/2023] [Accepted: 12/01/2023] [Indexed: 01/09/2024]
Abstract
Hybrid organic-inorganic perovskites exhibit remarkable potential as cost-effective and high-efficiency materials for photovoltaic applications. Their exceptional chemical tunability opens further routes for optimizing their optical and electronic properties through structural engineering. Nevertheless, the extraordinary softness of the lattice, stemming from its interconnected organic-inorganic composition, unveils formidable challenges in structural characterization. Here, by focusing on the quintessential methylammonium lead triiodide, MAPbI3, we combine first-principles modeling with high-resolution neutron scattering data to identify the key stationary points on its shallow potential energy landscape. This combined experimental and computational approach enables us to benchmark the performance of a collection of semilocal exchange-correlation functionals and to track the local distortions of the perovskite framework, hallmarked by the inelastic neutron scattering response of the organic cation. By conducting a thorough examination of structural distortions, we introduce the IKUR-PVP-1 structural data set. This data set contains nine mechanically stable structural models, each manifesting a distinct vibrational response. IKUR-PVP-1 constitutes a valuable resource for assessing thermal behavior in the low-temperature perovskite phase. In addition, it paves the way for the development of accurate force fields, enabling a comprehensive understanding of the interplay between the structure and dynamics in MAPbI3 and related hybrid perovskites.
Collapse
Affiliation(s)
- Kacper Drużbicki
- Materials
Physics Center, CSIC-UPV/EHU, Paseo Manuel de Lardizabal 5, Donostia-San Sebastian 20018, Spain
- Polish
Academy of Sciences, Centre of Molecular and Macromolecular Studies, Sienkiewicza 112, Lodz 90-363, Poland
| | - Pablo Gila-Herranz
- Materials
Physics Center, CSIC-UPV/EHU, Paseo Manuel de Lardizabal 5, Donostia-San Sebastian 20018, Spain
| | - Pelayo Marin-Villa
- Materials
Physics Center, CSIC-UPV/EHU, Paseo Manuel de Lardizabal 5, Donostia-San Sebastian 20018, Spain
| | - Mattia Gaboardi
- Materials
Physics Center, CSIC-UPV/EHU, Paseo Manuel de Lardizabal 5, Donostia-San Sebastian 20018, Spain
- C.S.G.I.
& Chemistry Department, University of
Pavia, Viale Taramelli,
16, Pavia 27100, Italy
| | - Jeff Armstrong
- ISIS
Neutron and Muon Facility, Rutherford Appleton
Laboratory, Didcot OX11 0QX, U.K.
| | - Felix Fernandez-Alonso
- Materials
Physics Center, CSIC-UPV/EHU, Paseo Manuel de Lardizabal 5, Donostia-San Sebastian 20018, Spain
- Donostia
International Physics Center (DIPC), Paseo Manuel de Lardizabal 4, Donostia-San
Sebastian 20018, Spain
- IKERBASQUE,
Basque Foundation for Science, Plaza Euskadi 5, Bilbao 48009, Spain
| |
Collapse
|
2
|
Hasegawa T. Nuclear Quantum Dynamics of Three-Dimensional Condensed-Phase Systems by Constant Uncertainty Molecular Dynamics. J Phys Chem Lett 2023; 14:8043-8049. [PMID: 37651720 DOI: 10.1021/acs.jpclett.3c01254] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 09/02/2023]
Abstract
Constant uncertainty molecular dynamics (CUMD), which was developed to include nuclear quantum effects (NQEs) in molecular dynamics (MD) simulations, is extended to three-dimensional condensed-phase systems. Its applicability was verified via CUMD simulations of bulk water and ice Ih with the q-TIP4P/F potential model. The simulated radial distribution functions, proton momentum distributions, and infrared spectra were compared with those of previous studies in which the simulations were carried out on the basis of the imaginary-time path-integral formalism. We found that CUMD can function as an alternative trajectory-based simulation method for including static and dynamic NQEs in MD simulations of three-dimensional condensed-phase systems.
Collapse
Affiliation(s)
- Taisuke Hasegawa
- Center for Computational Sciences, University of Tsukuba, 1-1-1 Tennodai, Tsukuba 305-8577, Japan
| |
Collapse
|
3
|
Videla PE, Batista VS. Matsubara dynamics approximation for generalized multi-time correlation functions. J Chem Phys 2023; 158:2889027. [PMID: 37154285 DOI: 10.1063/5.0146654] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/15/2023] [Accepted: 04/21/2023] [Indexed: 05/10/2023] Open
Abstract
We introduce a semi-classical approximation for calculating generalized multi-time correlation functions based on Matsubara dynamics, a classical dynamics approach that conserves the quantum Boltzmann distribution. This method is exact for the zero time and harmonic limits and reduces to classical dynamics when only one Matsubara mode is considered (i.e., the centroid). Generalized multi-time correlation functions can be expressed as canonical phase-space integrals, involving classically evolved observables coupled through Poisson brackets in a smooth Matsubara space. Numerical tests on a simple potential show that the Matsubara approximation exhibits better agreement with exact results than classical dynamics, providing a bridge between the purely quantum and classical descriptions of multi-time correlation functions. Despite the phase problem that prevents practical applications of Matsubara dynamics, the reported work provides a benchmark theory for the future development of quantum-Boltzmann-preserving semi-classical approximations for studies of chemical dynamics in condensed phase systems.
Collapse
Affiliation(s)
- Pablo E Videla
- Department of Chemistry, Yale University, 225 Prospect Street, New Haven, Connecticut 06520-8107, USA
| | - Victor S Batista
- Department of Chemistry, Yale University, 225 Prospect Street, New Haven, Connecticut 06520-8107, USA
| |
Collapse
|
4
|
Plé T, Mauger N, Adjoua O, Inizan TJ, Lagardère L, Huppert S, Piquemal JP. Routine Molecular Dynamics Simulations Including Nuclear Quantum Effects: From Force Fields to Machine Learning Potentials. J Chem Theory Comput 2023; 19:1432-1445. [PMID: 36856658 DOI: 10.1021/acs.jctc.2c01233] [Citation(s) in RCA: 9] [Impact Index Per Article: 9.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 03/02/2023]
Abstract
We report the implementation of a multi-CPU and multi-GPU massively parallel platform dedicated to the explicit inclusion of nuclear quantum effects (NQEs) in the Tinker-HP molecular dynamics (MD) package. The platform, denoted Quantum-HP, exploits two simulation strategies: the Ring-Polymer Molecular Dynamics (RPMD) that provides exact structural properties at the cost of a MD simulation in an extended space of multiple replicas and the adaptive Quantum Thermal Bath (adQTB) that imposes the quantum distribution of energy on a classical system via a generalized Langevin thermostat and provides computationally affordable and accurate (though approximate) NQEs. We discuss some implementation details, efficient numerical schemes, and parallelization strategies and quickly review the GPU acceleration of our code. Our implementation allows an efficient inclusion of NQEs in MD simulations for very large systems, as demonstrated by scaling tests on water boxes with more than 200,000 atoms (simulated using the AMOEBA polarizable force field). We test the compatibility of the approach with Tinker-HP's recently introduced Deep-HP machine learning potentials module by computing water properties using the DeePMD potential with adQTB thermostatting. Finally, we show that the platform is also compatible with the alchemical free energy estimation capabilities of Tinker-HP and fast enough to perform simulations. Therefore, we study how NQEs affect the hydration free energy of small molecules solvated with the recently developed Q-AMOEBA water force field. Overall, the Quantum-HP platform allows users to perform routine quantum MD simulations of large condensed-phase systems and will help to shed new light on the quantum nature of important interactions in biological matter.
Collapse
Affiliation(s)
- Thomas Plé
- Sorbonne Université, LCT, UMR 7616 CNRS, F-75005 Paris, France
| | - Nastasia Mauger
- Sorbonne Université, LCT, UMR 7616 CNRS, F-75005 Paris, France
| | - Olivier Adjoua
- Sorbonne Université, LCT, UMR 7616 CNRS, F-75005 Paris, France
| | | | - Louis Lagardère
- Sorbonne Université, LCT, UMR 7616 CNRS, F-75005 Paris, France
| | - Simon Huppert
- Institut des Nanosciences de Paris (INSP), CNRS UMR 7588, and Sorbonne Université, F-75005 Paris, France
| | - Jean-Philip Piquemal
- Sorbonne Université, LCT, UMR 7616 CNRS, F-75005 Paris, France.,Institut Universitaire de France, 75005 Paris, France.,Department of Biomedical Engineering, The University of Texas at Austin, Austin, Texas 78712, United States
| |
Collapse
|
5
|
Videla PE, Batista VS. An exact imaginary-time path-integral phase-space formulation of multi-time correlation functions. J Chem Phys 2023; 158:094101. [PMID: 36889944 DOI: 10.1063/5.0137898] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/10/2023] Open
Abstract
An exact representation of quantum mechanics using the language of phase-space variables provides a natural starting point to introduce and develop semiclassical approximations for the calculation of time correlation functions. Here, we introduce an exact path-integral formalism for calculations of multi-time quantum correlation functions as canonical averages over ring-polymer dynamics in imaginary time. The formulation provides a general formalism that exploits the symmetry of path integrals with respect to permutations in imaginary time, expressing correlations as products of imaginary-time-translation-invariant phase-space functions coupled through Poisson bracket operators. The method naturally recovers the classical limit of multi-time correlation functions and provides an interpretation of quantum dynamics in terms of "interfering trajectories" of the ring-polymer in phase space. The introduced phase-space formulation provides a rigorous framework for the future development of quantum dynamics methods that exploit the invariance of imaginary time path integrals to cyclic permutations.
Collapse
Affiliation(s)
- Pablo E Videla
- Department of Chemistry, Yale Energy Sciences Institute and Yale Quantum Institute, Yale University, 225 Prospect Street, New Haven, Connecticut 06520, USA
| | - Victor S Batista
- Department of Chemistry, Yale Energy Sciences Institute and Yale Quantum Institute, Yale University, 225 Prospect Street, New Haven, Connecticut 06520, USA
| |
Collapse
|
6
|
Gardner J, Douglas-Gallardo OA, Stark WG, Westermayr J, Janke SM, Habershon S, Maurer RJ. NQCDynamics.jl: A Julia package for nonadiabatic quantum classical molecular dynamics in the condensed phase. J Chem Phys 2022; 156:174801. [PMID: 35525649 DOI: 10.1063/5.0089436] [Citation(s) in RCA: 5] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/12/2023] Open
Abstract
Accurate and efficient methods to simulate nonadiabatic and quantum nuclear effects in high-dimensional and dissipative systems are crucial for the prediction of chemical dynamics in the condensed phase. To facilitate effective development, code sharing, and uptake of newly developed dynamics methods, it is important that software implementations can be easily accessed and built upon. Using the Julia programming language, we have developed the NQCDynamics.jl package, which provides a framework for established and emerging methods for performing semiclassical and mixed quantum-classical dynamics in the condensed phase. The code provides several interfaces to existing atomistic simulation frameworks, electronic structure codes, and machine learning representations. In addition to the existing methods, the package provides infrastructure for developing and deploying new dynamics methods, which we hope will benefit reproducibility and code sharing in the field of condensed phase quantum dynamics. Herein, we present our code design choices and the specific Julia programming features from which they benefit. We further demonstrate the capabilities of the package on two examples of chemical dynamics in the condensed phase: the population dynamics of the spin-boson model as described by a wide variety of semiclassical and mixed quantum-classical nonadiabatic methods and the reactive scattering of H2 on Ag(111) using the molecular dynamics with electronic friction method. Together, they exemplify the broad scope of the package to study effective model Hamiltonians and realistic atomistic systems.
Collapse
Affiliation(s)
- James Gardner
- Department of Chemistry, University of Warwick, Gibbet Hill Road, Coventry CV4 7AL, United Kingdom
| | - Oscar A Douglas-Gallardo
- Department of Chemistry, University of Warwick, Gibbet Hill Road, Coventry CV4 7AL, United Kingdom
| | - Wojciech G Stark
- Department of Chemistry, University of Warwick, Gibbet Hill Road, Coventry CV4 7AL, United Kingdom
| | - Julia Westermayr
- Department of Chemistry, University of Warwick, Gibbet Hill Road, Coventry CV4 7AL, United Kingdom
| | - Svenja M Janke
- Department of Chemistry, University of Warwick, Gibbet Hill Road, Coventry CV4 7AL, United Kingdom
| | - Scott Habershon
- Department of Chemistry, University of Warwick, Gibbet Hill Road, Coventry CV4 7AL, United Kingdom
| | - Reinhard J Maurer
- Department of Chemistry, University of Warwick, Gibbet Hill Road, Coventry CV4 7AL, United Kingdom
| |
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
|
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
|
9
|
Benson RL, Althorpe SC. On the "Matsubara heating" of overtone intensities and Fermi splittings. J Chem Phys 2021; 155:104107. [PMID: 34525826 DOI: 10.1063/5.0056829] [Citation(s) in RCA: 7] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/19/2022] Open
Abstract
Classical molecular dynamics (MD) and imaginary-time path-integral dynamics methods underestimate the infrared absorption intensities of overtone and combination bands by typically an order of magnitude. Plé et al. [J. Chem. Phys. 155, 2863 (2021)] have shown that this is because such methods fail to describe the coupling of the centroid to the Matsubara dynamics of the fluctuation modes; classical first-order perturbation theory (PT) applied to the Matsubara dynamics is sufficient to recover most of the lost intensity in simple models and gives identical results to quantum (Rayleigh-Schrödinger) PT. Here, we show numerically that the results of this analysis can be used as post-processing correction factors, which can be applied to realistic (classical MD or path-integral dynamics) simulations of infrared spectra. We find that the correction factors recover most of the lost intensity in the overtone and combination bands of gas-phase water and ammonia and much of it for liquid water. We then re-derive and confirm the earlier PT analysis by applying canonical PT to Matsubara dynamics, which has the advantage of avoiding secular terms and gives a simple picture of the perturbed Matsubara dynamics in terms of action-angle variables. Collectively, these variables "Matsubara heat" the amplitudes of the overtone and combination vibrations of the centroid to what they would be in a classical system with the oscillators (of frequency Ωi) held at their quantum effective temperatures [of ℏΩi coth(βℏΩi/2)/2kB]. Numerical calculations show that a similar neglect of "Matsubara heating" causes path-integral methods to underestimate Fermi resonance splittings.
Collapse
Affiliation(s)
- Raz L Benson
- 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
|
10
|
Morresi T, Paulatto L, Vuilleumier R, Casula M. Probing anharmonic phonons by quantum correlators: A path integral approach. J Chem Phys 2021; 154:224108. [PMID: 34241203 DOI: 10.1063/5.0050450] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
We devise an efficient scheme to determine vibrational properties from Path Integral Molecular Dynamics (PIMD) simulations. The method is based on zero-time Kubo-transformed correlation functions and captures the anharmonicity of the potential due to both temperature and quantum effects. Using analytical derivations and numerical calculations on toy-model potentials, we show that two different estimators built upon PIMD correlation functions fully characterize the phonon spectra and the anharmonicity strength. The first estimator is associated with the force-force quantum correlators and, in the weak anharmonic regime, yields reliable zero-point motion frequencies and thermodynamic properties of the quantum system. The second one is instead connected to displacement-displacement correlators and accurately probes the lowest-energy phonon excitations, regardless of the anharmonicity strength of the system. We also prove that the use of generalized eigenvalue equations, in place of the standard normal mode equations, leads to a significant speed-up in the PIMD phonon calculations, both in terms of faster convergence rate and smaller time step bias. Within this framework, using ab initio PIMD simulations, we compute phonon dispersions of diamond and of the high-pressure I41/amd phase of atomic hydrogen. We find that in the latter case, the anharmonicity is stronger than previously estimated and yields a sizeable red-shift in the vibrational spectrum of atomic hydrogen.
Collapse
Affiliation(s)
- T Morresi
- Institut de Minéralogie, de Physique des Matériaux et de Cosmochimie (IMPMC), Sorbonne Université, CNRS UMR 7590, MNHN, 4 Place Jussieu, 75252 Paris, France
| | - L Paulatto
- Institut de Minéralogie, de Physique des Matériaux et de Cosmochimie (IMPMC), Sorbonne Université, CNRS UMR 7590, MNHN, 4 Place Jussieu, 75252 Paris, France
| | - R Vuilleumier
- PASTEUR, Département de chimie, École normale supérieure, PSL University, Sorbonne Université, CNRS, 75005 Paris, France
| | - M Casula
- Institut de Minéralogie, de Physique des Matériaux et de Cosmochimie (IMPMC), Sorbonne Université, CNRS UMR 7590, MNHN, 4 Place Jussieu, 75252 Paris, France
| |
Collapse
|
11
|
Chowdhury SN, Huo P. Non-adiabatic Matsubara dynamics and non-adiabatic ring-polymer molecular dynamics. J Chem Phys 2021; 154:124124. [PMID: 33810665 DOI: 10.1063/5.0042136] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/06/2023] Open
Abstract
We present the non-adiabatic Matsubara dynamics, a general framework for computing the time-correlation function (TCF) of electronically non-adiabatic systems. This new formalism is derived based on the generalized Kubo-transformed TCF using the Wigner representation for both the nuclear degrees of freedom and the electronic mapping variables. By dropping the non-Matsubara nuclear normal modes in the quantum Liouvillian and explicitly integrating these modes out from the expression of the TCF, we derived the non-adiabatic Matsubara dynamics approach. Further making the approximation to drop the imaginary part of the Matsubara Liouvillian and enforce the nuclear momentum integral to be real, we arrived at the non-adiabatic ring-polymer molecular dynamics (NRPMD) approach. We have further justified the capability of NRPMD for simulating the non-equilibrium TCF. This work provides the rigorous theoretical foundation for several recently proposed state-dependent RPMD approaches and offers a general framework for developing new non-adiabatic quantum dynamics methods in the future.
Collapse
Affiliation(s)
- Sutirtha N Chowdhury
- Department of Chemistry, University of Rochester, 120 Trustee Road, Rochester, New York 14627, USA
| | - Pengfei Huo
- Department of Chemistry, University of Rochester, 120 Trustee Road, Rochester, New York 14627, USA
| |
Collapse
|
12
|
Jung KA, Videla PE, Batista VS. Ring-polymer, centroid, and mean-field approximations to multi-time Matsubara dynamics. J Chem Phys 2020; 153:124112. [PMID: 33003707 DOI: 10.1063/5.0021843] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Based on a recently developed generalization of Matsubara dynamics to the multi-time realm, we present a formal derivation of multi-time generalizations of ring-polymer molecular dynamics, thermostatted ring-polymer molecular dynamics (TRPMD), centroid molecular dynamics (CMD), and mean-field Matsubara dynamics. Additionally, we analyze the short-time accuracy of each methodology. We find that for multi-time correlation functions of linear operators, (T)RPMD is accurate up to order t3, while CMD is only correct up to t, indicating a degradation in the accuracy of these methodologies with respect to the single-time counterparts. The present work provides a firm justification for the use of path-integral-based approximations for the calculation of multi-time 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
|
13
|
Abstract
Dynamics at molecule-metal interfaces are a subject of intense current interest and come in many different flavors of experiments: gas-phase scattering, chemisorption, electrochemistry, nanojunction transport, and heterogeneous catalysis, to name a few. These dynamics involve nuclear degrees of freedom entangled with many electronic degrees of freedom (in the metal), and as such there is always the possibility for nonadiabatic phenomena to appear: the nuclei do not necessarily need to move slower than the electrons to break the Born-Oppenheimer (BO) approximation. In this Feature Article, we review a set of dynamical methods developed recently to deal with such nonadiabatic phenomena at a metal surface, methods that serve as alternatives to Tully's independent electron surface hopping (IESH) model. In the weak molecule-metal coupling regime, a classical master equation (CME) can be derived and a simple surface hopping approach is proposed to propagate nuclear and electronic dynamics stochastically. In the strong molecule-metal interaction regime, a Fokker-Planck equation can be derived for the nuclear dynamics, with electronic DoFs incorporated into the overall friction and random force. Lastly, a broadened classical master equation (BCME) can interpolate between the weak and strong molecule-metal interactions. Here, we briefly review these methods and the relevant benchmarking data, showing in particular how the methods can be used to calculate nonequilibrium transport properties. We highlight several open questions and pose several avenues for future study.
Collapse
Affiliation(s)
- Wenjie Dou
- Department of Chemistry , University of California, Berkeley , Berkeley , California 94720 , United States
| | - Joseph E Subotnik
- Department of Chemistry , University of Pennsylvania , Philadelphia , Pennsylvania 19104 , United States
| |
Collapse
|
14
|
Benson RL, Trenins G, Althorpe SC. Which quantum statistics–classical dynamics method is best for water? Faraday Discuss 2020; 221:350-366. [DOI: 10.1039/c9fd00077a] [Citation(s) in RCA: 20] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
We make a systematic comparison of TRPMD, CMD, QCMD and LSC-IVR by calculating the infrared spectrum of water in its gas, liquid and ice phases.
Collapse
Affiliation(s)
- Raz L. Benson
- Department of Chemistry
- University of Cambridge
- Cambridge
- UK
| | - George Trenins
- Department of Chemistry
- University of Cambridge
- Cambridge
- UK
| | | |
Collapse
|
15
|
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
|
16
|
Ekanayake NT, Garashchuk S. Evaluation of the quantum time-correlation functions employing the Hamilton–Jacobi dynamics framework. Theor Chem Acc 2018. [DOI: 10.1007/s00214-018-2398-6] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
|
17
|
Abstract
Electronic friction is a correction to the Born-Oppenheimer approximation, whereby nuclei in motion experience a drag in the presence of a manifold of electronic states. The notion of electronic friction has a long history and has been (re-)discovered in the context of a wide variety of different chemical and physical systems including, but not limited to, surface scattering events, surface reactions or chemisorption, electrochemistry, and conduction through molecular-(or nano-) junctions. Over the years, quite a few different forms of electronic friction have been offered in the literature. In this perspective, we briefly review these developments of electronic friction, highlighting the fact that we can now isolate a single, unifying form for (Markovian) electronic friction. We also focus on the role of electron-electron interactions for understanding frictional effects and offer our thoughts on the strengths and weaknesses of using electronic friction to model dynamics in general.
Collapse
Affiliation(s)
- Wenjie Dou
- Department of Chemistry, University of Pennsylvania, Philadelphia, Pennsylvania 19104, USA
| | - Joseph E Subotnik
- Department of Chemistry, University of Pennsylvania, Philadelphia, Pennsylvania 19104, USA
| |
Collapse
|
18
|
Church MS, Hele TJH, Ezra GS, Ananth N. Nonadiabatic semiclassical dynamics in the mixed quantum-classical initial value representation. J Chem Phys 2018; 148:102326. [DOI: 10.1063/1.5005557] [Citation(s) in RCA: 21] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Matthew S. Church
- Department of Chemistry and Chemical Biology, Cornell University, Ithaca, New York 14853, USA
| | - Timothy J. H. Hele
- Department of Chemistry and Chemical Biology, Cornell University, Ithaca, New York 14853, USA
| | - Gregory S. Ezra
- Department of Chemistry and Chemical Biology, Cornell University, Ithaca, New York 14853, USA
| | - Nandini Ananth
- Department of Chemistry and Chemical Biology, Cornell University, Ithaca, New York 14853, USA
| |
Collapse
|