Show simple item record

dc.contributor.authorWelsch, Ren
dc.contributor.authorSong, Ken
dc.contributor.authorShi, Qen
dc.contributor.authorAlthorpe, Stuarten
dc.contributor.authorMiller, TFen
dc.date.accessioned2017-02-02T14:03:43Z
dc.date.available2017-02-02T14:03:43Z
dc.date.issued2016-11-28en
dc.identifier.issn0021-9606
dc.identifier.urihttps://www.repository.cam.ac.uk/handle/1810/262237
dc.description.abstractWe investigate the calculation of approximate non-equilibrium quantum time correlation functions (TCFs) using two popular path-integral-based molecular dynamics methods, ring-polymer molecular dynamics (RPMD) and centroid molecular dynamics (CMD). It is shown that for the cases of a sudden vertical excitation and an initial momentum impulse, both RPMD and CMD yield non-equilibrium TCFs for linear operators that are exact for high temperatures, in the $\textit{t}$ = 0 limit, and for harmonic potentials; the subset of these conditions that are preserved for non-equilibrium TCFs of non-linear operators is also discussed. Furthermore, it is shown that for these non-equilibrium initial conditions, both methods retain the connection to Matsubara dynamics that has previously been established for equilibrium initial conditions. Comparison of non-equilibrium TCFs from RPMD and CMD to Matsubara dynamics at short times reveals the orders in time to which the methods agree. Specifically, for the position-autocorrelation function associated with sudden vertical excitation, RPMD and CMD agree with Matsubara dynamics up to O($\textit{t}$$^{4}$) and O($\textit{t}$$^{1}$), respectively; for the position-autocorrelation function associated with an initial momentum impulse, RPMD and CMD agree with Matsubara dynamics up to O($\textit{t}$$^{5}$) and O($\textit{t}$$^{2}$), respectively. Numerical tests using model potentials for a wide range of non-equilibrium initial conditions show that RPMD and CMD yield non-equilibrium TCFs with an accuracy that is comparable to that for equilibrium TCFs. RPMD is also used to investigate excited-state proton transfer in a system-bath model, and it is compared to numerically exact calculations performed using a recently developed version of the Liouville space hierarchical equation of motion approach; again, similar accuracy is observed for non-equilibrium and equilibrium initial conditions.
dc.description.sponsorshipR.W. acknowledges financial support from the Deutsche Forschungsgemeinschaft under Grant No. We 5762/1-1. S.C.A. acknowledges support from the Leverhulme Trust, the Chinese Academy of Sciences visiting professorship for senior international scientists (Grant No. 2013T2J0022) program, and the hospitality of the Miller group at Caltech. Q.S. acknowledges support from the National Natural Science Foundation of China (Grant No. 21290194). T.F.M. acknowledges support from the National Science Foundation (NSF) CAREER Award under Grant No. CHE-1057112 and from the Office of Naval Research (ONR) under Grant No. N00014-16-1-2761.
dc.languageengen
dc.language.isoenen
dc.publisherAmerican Institute of Physics
dc.titleNon-equilibrium dynamics from RPMD and CMDen
dc.typeArticle
prism.issueIdentifier20en
prism.number204118en
prism.publicationDate2016en
prism.publicationNameJournal of Chemical Physicsen
prism.volume145en
dc.identifier.doi10.17863/CAM.7492
dcterms.dateAccepted2016-11-07en
rioxxterms.versionofrecord10.1063/1.4967958en
rioxxterms.versionVoRen
rioxxterms.licenseref.urihttp://www.rioxx.net/licenses/all-rights-reserveden
rioxxterms.licenseref.startdate2016-11-28en
dc.contributor.orcidAlthorpe, Stuart [0000-0003-1288-8070]
dc.identifier.eissn1089-7690
rioxxterms.typeJournal Article/Reviewen
pubs.funder-project-idLeverhulme Trust (RF-2014-091)
cam.issuedOnline2016-11-30en
rioxxterms.freetoread.startdate2017-11-28


Files in this item

Thumbnail
Thumbnail

This item appears in the following Collection(s)

Show simple item record