1
|
Dutta R, Cabral DGA, Lyu N, Vu NP, Wang Y, Allen B, Dan X, Cortiñas RG, Khazaei P, Schäfer M, Albornoz ACCD, Smart SE, Nie S, Devoret MH, Mazziotti DA, Narang P, Wang C, Whitfield JD, Wilson AK, Hendrickson HP, Lidar DA, Pérez-Bernal F, Santos LF, Kais S, Geva E, Batista VS. Simulating Chemistry on Bosonic Quantum Devices. J Chem Theory Comput 2024. [PMID: 39068594 DOI: 10.1021/acs.jctc.4c00544] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 07/30/2024]
Abstract
Bosonic quantum devices offer a novel approach to realize quantum computations, where the quantum two-level system (qubit) is replaced with the quantum (an)harmonic oscillator (qumode) as the fundamental building block of the quantum simulator. The simulation of chemical structure and dynamics can then be achieved by representing or mapping the system Hamiltonians in terms of bosonic operators. In this Perspective, we review recent progress and future potential of using bosonic quantum devices for addressing a wide range of challenging chemical problems, including the calculation of molecular vibronic spectra, the simulation of gas-phase and solution-phase adiabatic and nonadiabatic chemical dynamics, the efficient solution of molecular graph theory problems, and the calculations of electronic structure.
Collapse
Affiliation(s)
- Rishab Dutta
- Department of Chemistry, Yale University, New Haven, Connecticut 06520, United States
| | - Delmar G A Cabral
- Department of Chemistry, Yale University, New Haven, Connecticut 06520, United States
| | - Ningyi Lyu
- Department of Chemistry, Yale University, New Haven, Connecticut 06520, United States
| | - Nam P Vu
- Department of Chemistry, Yale University, New Haven, Connecticut 06520, United States
- Department of Chemistry, Lafayette College, Easton, Pennsylvania 18042, United States
| | - Yuchen Wang
- Department of Chemistry, Department of Physics, and Purdue Quantum Science and Engineering Institute, Purdue University, West Lafayette, Indiana 47907, United States
| | - Brandon Allen
- Department of Chemistry, Yale University, New Haven, Connecticut 06520, United States
| | - Xiaohan Dan
- Department of Chemistry, Yale University, New Haven, Connecticut 06520, United States
| | - Rodrigo G Cortiñas
- Department of Applied Physics and Department of Physics, Yale University, New Haven, Connecticut 06520, United States
- Yale Quantum Institute, Yale University, New Haven, Connecticut 06511, United States
| | - Pouya Khazaei
- Department of Chemistry, University of Michigan, Ann Arbor, Michigan 48109, United States
| | - Max Schäfer
- Department of Applied Physics and Department of Physics, Yale University, New Haven, Connecticut 06520, United States
- Yale Quantum Institute, Yale University, New Haven, Connecticut 06511, United States
| | - Alejandro C C D Albornoz
- Department of Applied Physics and Department of Physics, Yale University, New Haven, Connecticut 06520, United States
- Yale Quantum Institute, Yale University, New Haven, Connecticut 06511, United States
| | - Scott E Smart
- Division of Physical Sciences, College of Letters and Science and Department of Electrical and Computer Engineering, University of California, Los Angeles, Los Angeles, California 90095, United States
| | - Scott Nie
- Division of Physical Sciences, College of Letters and Science and Department of Electrical and Computer Engineering, University of California, Los Angeles, Los Angeles, California 90095, United States
| | - Michel H Devoret
- Department of Applied Physics and Department of Physics, Yale University, New Haven, Connecticut 06520, United States
- Yale Quantum Institute, Yale University, New Haven, Connecticut 06511, United States
| | - David A Mazziotti
- Department of Chemistry and The James Franck Institute, The University of Chicago, Chicago, Illinois 60637, United States
| | - Prineha Narang
- Division of Physical Sciences, College of Letters and Science and Department of Electrical and Computer Engineering, University of California, Los Angeles, Los Angeles, California 90095, United States
| | - Chen Wang
- Department of Physics, University of Massachusetts - Amherst, Amherst, Massachusetts 01003, United States
| | - James D Whitfield
- Department of Physics and Astronomy, Dartmouth College, Hanover, New Hampshire 01003, United States
| | - Angela K Wilson
- Department of Chemistry, Michigan State University, East Lansing, Michigan 48864, United States
| | - Heidi P Hendrickson
- Department of Chemistry, Lafayette College, Easton, Pennsylvania 18042, United States
| | - Daniel A Lidar
- Department of Electrical & Computer Engineering, Department of Chemistry, Department of Physics & Astronomy, and Center for Quantum Information Science & Technology, University of Southern California, Los Angeles, California 90089, United States
| | - Francisco Pérez-Bernal
- Departamento de Ciencias Integradas y Centro de Estudios Avanzados en Física, Matemáticas y Computación, Universidad de Huelva, Huelva 21071, Spain
- Instituto Carlos I de Física Teórica y Computacional, Universidad de Granada, Granada 18071, Spain
| | - Lea F Santos
- Department of Physics, University of Connecticut, Storrs, Connecticut 06269, United States
| | - Sabre Kais
- Department of Chemistry, Department of Physics, and Purdue Quantum Science and Engineering Institute, Purdue University, West Lafayette, Indiana 47907, United States
| | - Eitan Geva
- Department of Chemistry, University of Michigan, Ann Arbor, Michigan 48109, United States
| | - Victor S Batista
- Department of Chemistry, Yale University, New Haven, Connecticut 06520, United States
- Yale Quantum Institute, Yale University, New Haven, Connecticut 06511, United States
| |
Collapse
|
2
|
Sayer T, Montoya-Castillo A. Generalized quantum master equations can improve the accuracy of semiclassical predictions of multitime correlation functions. J Chem Phys 2024; 161:011101. [PMID: 38949578 DOI: 10.1063/5.0219205] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/15/2024] [Accepted: 06/12/2024] [Indexed: 07/02/2024] Open
Abstract
Multitime quantum correlation functions are central objects in physical science, offering a direct link between the experimental observables and the dynamics of an underlying model. While experiments such as 2D spectroscopy and quantum control can now measure such quantities, the accurate simulation of such responses remains computationally expensive and sometimes impossible, depending on the system's complexity. A natural tool to employ is the generalized quantum master equation (GQME), which can offer computational savings by extending reference dynamics at a comparatively trivial cost. However, dynamical methods that can tackle chemical systems with atomistic resolution, such as those in the semiclassical hierarchy, often suffer from poor accuracy, limiting the credence one might lend to their results. By combining work on the accuracy-boosting formulation of semiclassical memory kernels with recent work on the multitime GQME, here we show for the first time that one can exploit a multitime semiclassical GQME to dramatically improve both the accuracy of coarse mean-field Ehrenfest dynamics and obtain orders of magnitude efficiency gains.
Collapse
Affiliation(s)
- Thomas Sayer
- Department of Chemistry, University of Colorado Boulder, Boulder, Colorado 80309, USA
| | | |
Collapse
|
3
|
Sayer T, Montoya-Castillo A. Efficient formulation of multitime generalized quantum master equations: Taming the cost of simulating 2D spectra. J Chem Phys 2024; 160:044108. [PMID: 38270238 DOI: 10.1063/5.0185578] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/31/2023] [Accepted: 12/26/2023] [Indexed: 01/26/2024] Open
Abstract
Modern 4-wave mixing spectroscopies are expensive to obtain experimentally and computationally. In certain cases, the unfavorable scaling of quantum dynamics problems can be improved using a generalized quantum master equation (GQME) approach. However, the inclusion of multiple (light-matter) interactions complicates the equation of motion and leads to seemingly unavoidable cubic scaling in time. In this paper, we present a formulation that greatly simplifies and reduces the computational cost of previous work that extended the GQME framework to treat arbitrary numbers of quantum measurements. Specifically, we remove the time derivatives of quantum correlation functions from the modified Mori-Nakajima-Zwanzig framework by switching to a discrete-convolution implementation inspired by the transfer tensor approach. We then demonstrate the method's capabilities by simulating 2D electronic spectra for the excitation-energy-transfer dimer model. In our method, the resolution of data can be arbitrarily coarsened, especially along the t2 axis, which mirrors how the data are obtained experimentally. Even in a modest case, this demands O(103) fewer data points. We are further able to decompose the spectra into one-, two-, and three-time correlations, showing how and when the system enters a Markovian regime where further measurements are unnecessary to predict future spectra and the scaling becomes quadratic. This offers the ability to generate long-time spectra using only short-time data, enabling access to timescales previously beyond the reach of standard methodologies.
Collapse
Affiliation(s)
- Thomas Sayer
- Department of Chemistry, University of Colorado Boulder, Boulder, Colorado 80309, USA
| | | |
Collapse
|
4
|
Amati G, Mannouch JR, Richardson JO. Detailed balance in mixed quantum-classical mapping approaches. J Chem Phys 2023; 159:214114. [PMID: 38054513 DOI: 10.1063/5.0176291] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/12/2023] [Accepted: 11/07/2023] [Indexed: 12/07/2023] Open
Abstract
The violation of detailed balance poses a serious problem for the majority of current quasiclassical methods for simulating nonadiabatic dynamics. In order to analyze the severity of the problem, we predict the long-time limits of the electronic populations according to various quasiclassical mapping approaches by applying arguments from classical ergodic theory. Our analysis confirms that regions of the mapping space that correspond to negative populations, which most mapping approaches introduce in order to go beyond the Ehrenfest approximation, pose the most serious issue for reproducing the correct thermalization behavior. This is because inverted potentials, which arise from negative electronic populations entering the nuclear force, can result in trajectories unphysically accelerating off to infinity. The recently developed mapping approach to surface hopping (MASH) provides a simple way of avoiding inverted potentials while retaining an accurate description of the dynamics. We prove that MASH, unlike any other quasiclassical approach, is guaranteed to describe the exact thermalization behavior of all quantum-classical systems, confirming it as one of the most promising methods for simulating nonadiabatic dynamics in real condensed-phase systems.
Collapse
Affiliation(s)
- Graziano Amati
- Department of Chemistry and Applied Biosciences, ETH Zürich, 8093 Zürich, Switzerland
| | - Jonathan R Mannouch
- Department of Chemistry and Applied Biosciences, ETH Zürich, 8093 Zürich, Switzerland
| | - Jeremy O Richardson
- Department of Chemistry and Applied Biosciences, ETH Zürich, 8093 Zürich, Switzerland
| |
Collapse
|
5
|
Dominic AJ, Cao S, Montoya-Castillo A, Huang X. Memory Unlocks the Future of Biomolecular Dynamics: Transformative Tools to Uncover Physical Insights Accurately and Efficiently. J Am Chem Soc 2023; 145:9916-9927. [PMID: 37104720 DOI: 10.1021/jacs.3c01095] [Citation(s) in RCA: 4] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 04/29/2023]
Abstract
Conformational changes underpin function and encode complex biomolecular mechanisms. Gaining atomic-level detail of how such changes occur has the potential to reveal these mechanisms and is of critical importance in identifying drug targets, facilitating rational drug design, and enabling bioengineering applications. While the past two decades have brought Markov state model techniques to the point where practitioners can regularly use them to glimpse the long-time dynamics of slow conformations in complex systems, many systems are still beyond their reach. In this Perspective, we discuss how including memory (i.e., non-Markovian effects) can reduce the computational cost to predict the long-time dynamics in these complex systems by orders of magnitude and with greater accuracy and resolution than state-of-the-art Markov state models. We illustrate how memory lies at the heart of successful and promising techniques, ranging from the Fokker-Planck and generalized Langevin equations to deep-learning recurrent neural networks and generalized master equations. We delineate how these techniques work, identify insights that they can offer in biomolecular systems, and discuss their advantages and disadvantages in practical settings. We show how generalized master equations can enable the investigation of, for example, the gate-opening process in RNA polymerase II and demonstrate how our recent advances tame the deleterious influence of statistical underconvergence of the molecular dynamics simulations used to parameterize these techniques. This represents a significant leap forward that will enable our memory-based techniques to interrogate systems that are currently beyond the reach of even the best Markov state models. We conclude by discussing some current challenges and future prospects for how exploiting memory will open the door to many exciting opportunities.
Collapse
Affiliation(s)
- Anthony J Dominic
- Department of Chemistry, University of Colorado Boulder, Boulder, Colorado 80309, USA
| | - Siqin Cao
- Department of Chemistry, Theoretical Chemistry Institute, University of Wisconsin-Madison, Madison, Wisconsin 53706, USA
| | | | - Xuhui Huang
- Department of Chemistry, Theoretical Chemistry Institute, University of Wisconsin-Madison, Madison, Wisconsin 53706, USA
| |
Collapse
|
6
|
Montoya-Castillo A, Markland TE. A derivation of the conditions under which bosonic operators exactly capture fermionic structure and dynamics. J Chem Phys 2023; 158:094112. [PMID: 36889969 DOI: 10.1063/5.0138664] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/15/2023] Open
Abstract
The dynamics of many-body fermionic systems are important in problems ranging from catalytic reactions at electrochemical surfaces to transport through nanojunctions and offer a prime target for quantum computing applications. Here, we derive the set of conditions under which fermionic operators can be exactly replaced by bosonic operators that render the problem amenable to a large toolbox of dynamical methods while still capturing the correct dynamics of n-body operators. Importantly, our analysis offers a simple guide on how one can exploit these simple maps to calculate nonequilibrium and equilibrium single- and multi-time correlation functions essential in describing transport and spectroscopy. We use this to rigorously analyze and delineate the applicability of simple yet effective Cartesian maps that have been shown to correctly capture the correct fermionic dynamics in select models of nanoscopic transport. We illustrate our analytical results with exact simulations of the resonant level model. Our work provides new insights as to when one can leverage the simplicity of bosonic maps to simulate the dynamics of many-electron systems, especially those where an atomistic representation of nuclear interactions becomes essential.
Collapse
Affiliation(s)
| | - Thomas E Markland
- Department of Chemistry, Stanford University, Stanford, California 94305, USA
| |
Collapse
|
7
|
Lyu N, Mulvihill E, Soley MB, Geva E, Batista VS. Tensor-Train Thermo-Field Memory Kernels for Generalized Quantum Master Equations. J Chem Theory Comput 2023; 19:1111-1129. [PMID: 36719350 DOI: 10.1021/acs.jctc.2c00892] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/01/2023]
Abstract
The generalized quantum master equation (GQME) approach provides a rigorous framework for deriving the exact equation of motion for any subset of electronic reduced density matrix elements (e.g., the diagonal elements). In the context of electronic dynamics, the memory kernel and inhomogeneous term of the GQME introduce the implicit coupling to nuclear motion and dynamics of electronic density matrix elements that are projected out (e.g., the off-diagonal elements), allowing for efficient quantum dynamics simulations. Here, we focus on benchmark quantum simulations of electronic dynamics in a spin-boson model system described by various types of GQMEs. Exact memory kernels and inhomogeneous terms are obtained from short-time quantum-mechanically exact tensor-train thermo-field dynamics (TT-TFD) simulations and are compared with those obtained from an approximate linearized semiclassical method, allowing for assessment of the accuracy of these approximate memory kernels and inhomogeneous terms. Moreover, we have analyzed the computational cost of the full and reduced-dimensionality GQMEs. The scaling of the computational cost is dependent on several factors, sometimes with opposite scaling trends. The TT-TFD memory kernels can provide insights on the main sources of inaccuracies of GQME approaches when combined with approximate input methods and pave the road for the development of quantum circuits that implement GQMEs on digital quantum computers.
Collapse
Affiliation(s)
- Ningyi Lyu
- Department of Chemistry, Yale University, New Haven, Connecticut 06520, United States
| | - Ellen Mulvihill
- Department of Chemistry, Yale University, New Haven, Connecticut 06520, United States
| | - Micheline B Soley
- Department of Chemistry, Yale University, New Haven, Connecticut 06520, United States.,Yale Quantum Institute, Yale University, New Haven, Connecticut 06511, United States.,Department of Chemistry, University of Wisconsin-Madison, Madison, Wisconsin 53706, United States
| | - Eitan Geva
- Department of Chemistry, University of Michigan, Ann Arbor, Michigan 48109, United States
| | - Victor S Batista
- Department of Chemistry, Yale University, New Haven, Connecticut 06520, United States.,Yale Quantum Institute, Yale University, New Haven, Connecticut 06511, United States
| |
Collapse
|
8
|
Atsango AO, Montoya-Castillo A, Markland TE. An accurate and efficient Ehrenfest dynamics approach for calculating linear and nonlinear electronic spectra. J Chem Phys 2023; 158:074107. [PMID: 36813724 DOI: 10.1063/5.0138671] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/05/2023] Open
Abstract
Linear and nonlinear electronic spectra provide an important tool to probe the absorption and transfer of electronic energy. Here, we introduce a pure state Ehrenfest approach to obtain accurate linear and nonlinear spectra that is applicable to systems with large numbers of excited states and complex chemical environments. We achieve this by representing the initial conditions as sums of pure states and unfolding multi-time correlation functions into the Schrödinger picture. By doing this, we show that one can obtain significant improvements in accuracy over the previously used projected Ehrenfest approach and that these benefits are particularly pronounced in cases where the initial condition is a coherence between excited states. While such initial conditions do not arise when calculating linear electronic spectra, they play a vital role in capturing multidimensional spectroscopies. We demonstrate the performance of our method by showing that it is able to quantitatively capture the exact linear, 2D electronic spectroscopy, and pump-probe spectra for a Frenkel exciton model in slow bath regimes and is even able to reproduce the main spectral features in fast bath regimes.
Collapse
Affiliation(s)
- Austin O Atsango
- Department of Chemistry, Stanford University, Stanford, California 94305, USA
| | | | - Thomas E Markland
- Department of Chemistry, Stanford University, Stanford, California 94305, USA
| |
Collapse
|
9
|
Sayer T, Montoya-Castillo A. Compact and complete description of non-Markovian dynamics. J Chem Phys 2023; 158:014105. [PMID: 36610963 DOI: 10.1063/5.0132614] [Citation(s) in RCA: 3] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/06/2023] Open
Abstract
Generalized master equations provide a theoretically rigorous framework to capture the dynamics of processes ranging from energy harvesting in plants and photovoltaic devices to qubit decoherence in quantum technologies and even protein folding. At their center is the concept of memory. The explicit time-nonlocal description of memory is both protracted and elaborate. When physical intuition is at a premium, one would desire a more compact, yet complete, description. Here, we demonstrate how and when the time-convolutionless formalism constitutes such a description. In particular, by focusing on the dissipative dynamics of the spin-boson and Frenkel exciton models, we show how to: easily construct the time-local generator from reference reduced dynamics, elucidate the dependence of its existence on the system parameters and the choice of reduced observables, identify the physical origin of its apparent divergences, and offer analysis tools to diagnose their severity and circumvent their deleterious effects. We demonstrate that, when applicable, the time-local approach requires as little information as the more commonly used time-nonlocal scheme, with the important advantages of providing a more compact description, greater algorithmic simplicity, and physical interpretability. We conclude by introducing the discrete-time analog and a straightforward protocol to employ it in cases where the reference dynamics have limited resolution. The insights we present here offer the potential for extending the reach of dynamical methods, reducing both their cost and conceptual complexity.
Collapse
Affiliation(s)
- Thomas Sayer
- Department of Chemistry, University of Colorado Boulder, Boulder, Colorado 80309, USA
| | | |
Collapse
|
10
|
Amati G, Saller MAC, Kelly A, Richardson JO. Quasiclassical approaches to the generalized quantum master equation. J Chem Phys 2022; 157:234103. [PMID: 36550031 DOI: 10.1063/5.0124028] [Citation(s) in RCA: 7] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022] Open
Abstract
The formalism of the generalized quantum master equation (GQME) is an effective tool to simultaneously increase the accuracy and the efficiency of quasiclassical trajectory methods in the simulation of nonadiabatic quantum dynamics. The GQME expresses correlation functions in terms of a non-Markovian equation of motion, involving memory kernels that are typically fast-decaying and can therefore be computed by short-time quasiclassical trajectories. In this paper, we study the approximate solution of the GQME, obtained by calculating the kernels with two methods: Ehrenfest mean-field theory and spin-mapping. We test the approaches on a range of spin-boson models with increasing energy bias between the two electronic levels and place a particular focus on the long-time limits of the populations. We find that the accuracy of the predictions of the GQME depends strongly on the specific technique used to calculate the kernels. In particular, spin-mapping outperforms Ehrenfest for all the systems studied. The problem of unphysical negative electronic populations affecting spin-mapping is resolved by coupling the method with the master equation. Conversely, Ehrenfest in conjunction with the GQME can predict negative populations, despite the fact that the populations calculated from direct dynamics are positive definite.
Collapse
Affiliation(s)
- Graziano Amati
- Laboratory of Physical Chemistry, ETH Zürich, 8093 Zürich, Switzerland
| | | | - Aaron Kelly
- The Hamburg Center for Ultrafast Imaging, Universität Hamburg, Max Planck Institute for the Structure and Dynamics of Matter and Center for Free-Electron Laser Science, 22761 Hamburg, Germany
| | | |
Collapse
|
11
|
Lai Y, Geva E. Electronic Absorption Spectra from Off-Diagonal Quantum Master Equations. J Chem Phys 2022; 157:104115. [DOI: 10.1063/5.0106888] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Abstract
Quantum master equations (QMEs) provide a general framework for describing electronic dynamics within a complex molecular system. Off-diagonal QMEs (OD-QMEs) correspond to a family of QMEs that describe the electronic dynamics in the interaction picture based on treating the off-diagonal coupling terms between electronic states as a small perturbation within the framework of second-order perturbation theory. The fact that OD-QMEs are given in terms of the interaction picture makes it non-trivial to obtain Schrodinger picture electronic coherences from them. A key experimental quantity that relies on the ability to obtain accurate Schrodinger picture electronic coherences is the absorption spectrum. In this paper, we propose using a recently introduced procedure for extracting Schrodinger picture electronic coherences from interaction picture inputs to calculate electronic absorption spectra from electronic dynamics generated by OD-QMEs. The accuracy of the absorption spectra obtained in this way is studied in the context of a biexciton benchmark model, by comparing spectra calculated based on time-local and time-nonlocal OD-QMEs to spectra calculated based on a Redfield-type QME and the non-perturbative and quantum-mechanically exact hierarchical equations of motion (HEOM) method.
Collapse
Affiliation(s)
- Yifan Lai
- Chemistry, University of Michigan, United States of America
| | - Eitan Geva
- Department of Chemistry, University of Michigan Department of Chemistry, United States of America
| |
Collapse
|
12
|
Runeson JE, Lawrence JE, Mannouch JR, Richardson JO. Explaining the Efficiency of Photosynthesis: Quantum Uncertainty or Classical Vibrations? J Phys Chem Lett 2022; 13:3392-3399. [PMID: 35404611 PMCID: PMC9036581 DOI: 10.1021/acs.jpclett.2c00538] [Citation(s) in RCA: 11] [Impact Index Per Article: 5.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/22/2022] [Accepted: 04/01/2022] [Indexed: 06/14/2023]
Abstract
Photosynthetic organisms are known to use a mechanism of vibrationally assisted exciton energy transfer to efficiently harvest energy from light. The importance of quantum effects in this mechanism is a long-standing topic of debate, which has traditionally focused on the role of excitonic coherences. Here, we address another recent claim: that the efficient energy transfer in the Fenna-Matthews-Olson complex relies on nuclear quantum uncertainty and would not function if the vibrations were classical. We present a counter-example to this claim, showing by trajectory-based simulations that a description in terms of quantum electrons and classical nuclei is indeed sufficient to describe the funneling of energy to the reaction center. We analyze and compare these findings to previous classical-nuclear approximations that predicted the absence of an energy funnel and conclude that the key difference and the reason for the discrepancy is the ability of the trajectories to properly account for Newton's third law.
Collapse
|
13
|
Dan X, Xu M, Yan Y, Shi Q. Generalized master equation for charge transport in a molecular junction: Exact memory kernels and their high order expansion. J Chem Phys 2022; 156:134114. [PMID: 35395901 DOI: 10.1063/5.0086663] [Citation(s) in RCA: 1] [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 derive a set of generalized master equations (GMEs) to study charge transport dynamics in molecular junctions using the Nakajima-Zwanzig-Mori projection operator approach. In the new GME, time derivatives of population on each quantum state of the molecule, as well as the tunneling current, are calculated as the convolution of time non-local memory kernels with populations on all system states. The non-Markovian memory kernels are obtained by combining the hierarchical equations of motion (HEOM) method and a previous derived Dyson relation for the exact kernel. A perturbative expansion of these memory kernels is then calculated using the extended HEOM developed in our previous work [M. Xu et al., J. Chem. Phys. 146, 064102 (2017)]. By using the resonant level model and the Anderson impurity model, we study properties of the exact memory kernels and analyze convergence properties of their perturbative expansions with respect to the system-bath coupling strength and the electron-electron repulsive energy. It is found that exact memory kernels calculated from HEOM exhibit short memory times and decay faster than the population and current dynamics. The high order perturbation expansion of the memory kernels can give converged results in certain parameter regimes. The Padé and Landau-Zener resummation schemes are also found to give improved results over low order perturbation theory.
Collapse
Affiliation(s)
- Xiaohan Dan
- 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 and University of Chinese Academy of Sciences, Beijing 100049, China
| | - Meng Xu
- 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 and University of Chinese Academy of Sciences, Beijing 100049, China
| | - Yaming Yan
- 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 and University of Chinese Academy of Sciences, Beijing 100049, 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 and University of Chinese Academy of Sciences, Beijing 100049, China
| |
Collapse
|
14
|
Mulvihill E, Geva E. Simulating the dynamics of electronic observables via reduced-dimensionality generalized quantum master equations. J Chem Phys 2022; 156:044119. [DOI: 10.1063/5.0078040] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/06/2023] Open
Affiliation(s)
- Ellen Mulvihill
- Department of Chemistry, University of Michigan, Ann Arbor, Michigan 48109, USA
| | - Eitan Geva
- Department of Chemistry, University of Michigan, Ann Arbor, Michigan 48109, USA
| |
Collapse
|
15
|
Lai Y, Geva E. On simulating the dynamics of electronic populations and coherences via quantum master equations based on treating off-diagonal electronic coupling terms as a small perturbation. J Chem Phys 2021; 155:204101. [PMID: 34852488 DOI: 10.1063/5.0069313] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Quantum master equations provide a general framework for describing the dynamics of electronic observables within a complex molecular system. One particular family of such equations is based on treating the off-diagonal coupling terms between electronic states as a small perturbation within the framework of second-order perturbation theory. In this paper, we show how different choices of projection operators, as well as whether one starts out with the time-convolution or the time-convolutionless forms of the generalized quantum master equation, give rise to four different types of such off-diagonal quantum master equations (OD-QMEs), namely, time-convolution and time-convolutionless versions of a Pauli-type OD-QME for only the electronic populations and an OD-QME for the full electronic density matrix (including both electronic populations and coherences). The fact that those OD-QMEs are given in terms of the interaction picture makes it non-trivial to obtain Schrödinger picture electronic coherences from them. To address this, we also extend a procedure for extracting Schrödinger picture electronic coherences from interaction picture populations recently introduced by Trushechkin in the context of time-convolutionless Pauli-type OD-QME to the other three types of OD-QMEs. The performance of the aforementioned four types of OD-QMEs is explored in the context of the Garg-Onuchic-Ambegaokar benchmark model for charge transfer in the condensed phase across a relatively wide parameter range. The results show that time-convolution OD-QMEs can be significantly more accurate than their time-convolutionless counterparts, particularly in the case of Pauli-type OD-QMEs, and that rather accurate Schrödinger picture coherences can be obtained from interaction picture electronic inputs.
Collapse
Affiliation(s)
- Yifan Lai
- Department of Chemistry, University of Michigan, Ann Arbor, Michigan 48109, USA
| | - Eitan Geva
- Department of Chemistry, University of Michigan, Ann Arbor, Michigan 48109, USA
| |
Collapse
|
16
|
Mathematical Models with Nonlocal Initial Conditions: An Exemplification from Quantum Mechanics. MATHEMATICAL AND COMPUTATIONAL APPLICATIONS 2021. [DOI: 10.3390/mca26040073] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/19/2022]
Abstract
Nonlocal models are ubiquitous in all branches of science and engineering, with a rapidly expanding range of mathematical and computational applications due to the ability of such models to capture effects and phenomena that traditional models cannot. While spatial nonlocalities have received considerable attention in the research community, the same cannot be said about nonlocality in time, in particular when nonlocal initial conditions are present. This paper aims at filling this gap, providing an overview of the current status of nonlocal models and focusing on the mathematical treatment of such models when nonlocal initial conditions are at the heart of the problem. Specifically, our representative example is given for a nonlocal-in-time problem for the abstract Schrödinger equation. By exploiting the linear nature of nonlocal conditions, we derive an exact representation of the solution operator under assumptions that the spectrum of Hamiltonian is contained in the horizontal strip of the complex plane. The derived representation permits us to establish the necessary and sufficient conditions for the problem’s well-posedness and the existence of its solution under different regularities. Furthermore, we present new sufficient conditions for the existence of the solution that extend the existing results in this field to the case when some nonlocal parameters are unbounded. Two further examples demonstrate the developed methodology and highlight the importance of its computer algebra component in the reduction procedures and parameter estimations for nonlocal models. Finally, a connection of the considered models and developed analysis is discussed in the context of other reduction techniques, concentrating on the most promising from the viewpoint of data-driven modelling environments, and providing directions for further generalizations.
Collapse
|
17
|
Brian D, Sun X. Generalized quantum master equation: A tutorial review and recent advances. CHINESE J CHEM PHYS 2021. [DOI: 10.1063/1674-0068/cjcp2109157] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/03/2023]
Affiliation(s)
- Dominikus Brian
- Division of Arts and Sciences, NYU Shanghai, Shanghai 200122, China
- NYU-ECNU Center for Computational Chemistry at NYU Shanghai, Shanghai 200062, China
- Department of Chemistry, New York University, New York 10003, USA
| | - Xiang Sun
- Division of Arts and Sciences, NYU Shanghai, Shanghai 200122, China
- NYU-ECNU Center for Computational Chemistry at NYU Shanghai, Shanghai 200062, China
- Department of Chemistry, New York University, New York 10003, USA
- State Key Laboratory of Precision Spectroscopy, East China Normal University, Shanghai 200241, China
| |
Collapse
|
18
|
Mulvihill E, Geva E. A Road Map to Various Pathways for Calculating the Memory Kernel of the Generalized Quantum Master Equation. J Phys Chem B 2021; 125:9834-9852. [PMID: 34424700 DOI: 10.1021/acs.jpcb.1c05719] [Citation(s) in RCA: 11] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
The generalized quantum master equation (GQME) provides a powerful framework for simulating electronic energy, charge, and coherence transfer dynamics in molecular systems. Within this framework, the effect of the nuclear degrees of freedom on the time evolution of the electronic reduced density matrix is fully captured by a memory kernel superoperator. However, the actual memory kernel depends on the choice of projection operator and is therefore not unique. Furthermore, calculating the memory kernel can be done in multiple ways that use different forms of projection-free inputs. Although the electronic dynamics is invariant to those choices when quantum-mechanically exact projection-free inputs are used, this is not the case when they are obtained via more feasible semiclassical or mixed quantum-classical approximate methods. Furthermore, the accuracy and numerical stability of the resulting electronic dynamics has been observed to be sensitive to the above-mentioned choices when approximate methods are used to calculate the projection-free inputs. In this article, we provide a systematic road map to 30 possible pathways for calculating the memory kernel and highlight how they are related as well as the ways in which they differ. We also compare the performance of different pathways in the context of the spin-boson benchmark model, with the projection-free inputs obtained via a mapping Hamiltonian linearized semiclassical method. In this case, we find that expressing the memory kernel with an exponential operator where the projection operator precedes the Liouvillian yields the most accurate and most numerically stable results.
Collapse
Affiliation(s)
- Ellen Mulvihill
- Department of Chemistry, University of Michigan, Ann Arbor, Michigan48109, United States
| | - Eitan Geva
- Department of Chemistry, University of Michigan, Ann Arbor, Michigan48109, United States
| |
Collapse
|
19
|
Mulvihill E, Lenn KM, Gao X, Schubert A, Dunietz BD, Geva E. Simulating energy transfer dynamics in the Fenna-Matthews-Olson complex via the modified generalized quantum master equation. J Chem Phys 2021; 154:204109. [PMID: 34241158 DOI: 10.1063/5.0051101] [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/15/2022] Open
Abstract
The generalized quantum master equation (GQME) provides a general and formally exact framework for simulating the reduced dynamics of open quantum systems. The recently introduced modified approach to the GQME (M-GQME) corresponds to a specific implementation of the GQME that is geared toward simulating the dynamics of the electronic reduced density matrix in systems governed by an excitonic Hamiltonian. Such a Hamiltonian, which is often used for describing energy and charge transfer dynamics in complex molecular systems, is given in terms of diabatic electronic states that are coupled to each other and correspond to different nuclear Hamiltonians. Within the M-GQME approach, the effect of the nuclear degrees of freedom on the time evolution of the electronic density matrix is fully captured by a memory kernel superoperator, which can be obtained from short-lived (compared to the time scale of energy/charge transfer) projection-free inputs. In this paper, we test the ability of the M-GQME to predict the energy transfer dynamics within a seven-state benchmark model of the Fenna-Matthews-Olson (FMO) complex, with the short-lived projection-free inputs obtained via the Ehrenfest method. The M-GQME with Ehrenfest-based inputs is shown to yield accurate results across a wide parameter range. It is also found to dramatically outperform the direct application of the Ehrenfest method and to provide better-behaved convergence with respect to memory time in comparison to an alternative implementation of the GQME approach previously applied to the same FMO model.
Collapse
Affiliation(s)
- Ellen Mulvihill
- Department of Chemistry, University of Michigan, Ann Arbor, Michigan 48109, USA
| | - Kristina M Lenn
- Department of Chemistry, University of Michigan, Ann Arbor, Michigan 48109, USA
| | - Xing Gao
- Department of Chemistry, University of Michigan, Ann Arbor, Michigan 48109, USA
| | - Alexander Schubert
- Department of Chemistry, University of Michigan, Ann Arbor, Michigan 48109, USA
| | - Barry D Dunietz
- Department of Chemistry and Biochemistry, Kent State University, Kent, Ohio 44242, USA
| | - Eitan Geva
- Department of Chemistry, University of Michigan, Ann Arbor, Michigan 48109, USA
| |
Collapse
|
20
|
Brian D, Liu Z, Dunietz BD, Geva E, Sun X. Three-state harmonic models for photoinduced charge transfer. J Chem Phys 2021; 154:174105. [PMID: 34241055 DOI: 10.1063/5.0050289] [Citation(s) in RCA: 9] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/29/2022] Open
Abstract
A widely used strategy for simulating the charge transfer between donor and acceptor electronic states in an all-atom anharmonic condensed-phase system is based on invoking linear response theory to describe the system in terms of an effective spin-boson model Hamiltonian. Extending this strategy to photoinduced charge transfer processes requires also taking into consideration the ground electronic state in addition to the excited donor and acceptor electronic states. In this paper, we revisit the problem of describing such nonequilibrium processes in terms of an effective three-state harmonic model. We do so within the framework of nonequilibrium Fermi's golden rule (NE-FGR) in the context of photoinduced charge transfer in the carotenoid-porphyrin-C60 (CPC60) molecular triad dissolved in explicit tetrahydrofuran (THF). To this end, we consider different ways for obtaining a three-state harmonic model from the equilibrium autocorrelation functions of the donor-acceptor, donor-ground, and acceptor-ground energy gaps, as obtained from all-atom molecular dynamics simulations of the CPC60/THF system. The quantum-mechanically exact time-dependent NE-FGR rate coefficients for two different charge transfer processes in two different triad conformations are then calculated using the effective three-state model Hamiltonians as well as a hierarchy of more approximate expressions that lead to the instantaneous Marcus theory limit. Our results show that the photoinduced charge transfer in CPC60/THF can be described accurately by the effective harmonic three-state models and that nuclear quantum effects are small in this system.
Collapse
Affiliation(s)
- Dominikus Brian
- Division of Arts and Sciences, NYU Shanghai, 1555 Century Avenue, Shanghai 200122, China
| | - Zengkui Liu
- Division of Arts and Sciences, NYU Shanghai, 1555 Century Avenue, Shanghai 200122, China
| | - Barry D Dunietz
- Department of Chemistry and Biochemistry, Kent State University, Kent, Ohio 44242, USA
| | - Eitan Geva
- Department of Chemistry, University of Michigan, Ann Arbor, Michigan 48109, USA
| | - Xiang Sun
- Division of Arts and Sciences, NYU Shanghai, 1555 Century Avenue, Shanghai 200122, China
| |
Collapse
|
21
|
Cao S, Montoya-Castillo A, Wang W, Markland TE, Huang X. On the advantages of exploiting memory in Markov state models for biomolecular dynamics. J Chem Phys 2021; 153:014105. [PMID: 32640825 DOI: 10.1063/5.0010787] [Citation(s) in RCA: 26] [Impact Index Per Article: 8.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/24/2023] Open
Abstract
Biomolecular dynamics play an important role in numerous biological processes. Markov State Models (MSMs) provide a powerful approach to study these dynamic processes by predicting long time scale dynamics based on many short molecular dynamics (MD) simulations. In an MSM, protein dynamics are modeled as a kinetic process consisting of a series of Markovian transitions between different conformational states at discrete time intervals (called "lag time"). To achieve this, a master equation must be constructed with a sufficiently long lag time to allow interstate transitions to become truly Markovian. This imposes a major challenge for MSM studies of proteins since the lag time is bound by the length of relatively short MD simulations available to estimate the frequency of transitions. Here, we show how one can employ the generalized master equation formalism to obtain an exact description of protein conformational dynamics both at short and long time scales without the time resolution restrictions imposed by the MSM lag time. Using a simple kinetic model, alanine dipeptide, and WW domain, we demonstrate that it is possible to construct these quasi-Markov State Models (qMSMs) using MD simulations that are 5-10 times shorter than those required by MSMs. These qMSMs only contain a handful of metastable states and, thus, can greatly facilitate the interpretation of mechanisms associated with protein dynamics. A qMSM opens the door to the study of conformational changes of complex biomolecules where a Markovian model with a few states is often difficult to construct due to the limited length of available MD simulations.
Collapse
Affiliation(s)
- Siqin Cao
- Department of Chemistry, The Hong Kong University of Science and Technology, Kowloon, Hong Kong
| | | | - Wei Wang
- Department of Chemistry, The Hong Kong University of Science and Technology, Kowloon, Hong Kong
| | - Thomas E Markland
- Department of Chemistry, Stanford University, Stanford, California 94305, USA
| | - Xuhui Huang
- Department of Chemistry, The Hong Kong University of Science and Technology, Kowloon, Hong Kong
| |
Collapse
|
22
|
Head-Marsden K, Flick J, Ciccarino CJ, Narang P. Quantum Information and Algorithms for Correlated Quantum Matter. Chem Rev 2020; 121:3061-3120. [PMID: 33326218 DOI: 10.1021/acs.chemrev.0c00620] [Citation(s) in RCA: 35] [Impact Index Per Article: 8.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/13/2022]
Abstract
Discoveries in quantum materials, which are characterized by the strongly quantum-mechanical nature of electrons and atoms, have revealed exotic properties that arise from correlations. It is the promise of quantum materials for quantum information science superimposed with the potential of new computational quantum algorithms to discover new quantum materials that inspires this Review. We anticipate that quantum materials to be discovered and developed in the next years will transform the areas of quantum information processing including communication, storage, and computing. Simultaneously, efforts toward developing new quantum algorithmic approaches for quantum simulation and advanced calculation methods for many-body quantum systems enable major advances toward functional quantum materials and their deployment. The advent of quantum computing brings new possibilities for eliminating the exponential complexity that has stymied simulation of correlated quantum systems on high-performance classical computers. Here, we review new algorithms and computational approaches to predict and understand the behavior of correlated quantum matter. The strongly interdisciplinary nature of the topics covered necessitates a common language to integrate ideas from these fields. We aim to provide this common language while weaving together fields across electronic structure theory, quantum electrodynamics, algorithm design, and open quantum systems. Our Review is timely in presenting the state-of-the-art in the field toward algorithms with nonexponential complexity for correlated quantum matter with applications in grand-challenge problems. Looking to the future, at the intersection of quantum information science and algorithms for correlated quantum matter, we envision seminal advances in predicting many-body quantum states and describing excitonic quantum matter and large-scale entangled states, a better understanding of high-temperature superconductivity, and quantifying open quantum system dynamics.
Collapse
Affiliation(s)
- Kade Head-Marsden
- John A. Paulson School of Engineering and Applied Sciences, Harvard University, Cambridge, Massachusetts 02138, United States
| | - Johannes Flick
- Center for Computational Quantum Physics, Flatiron Institute, New York, New York 10010, United States
| | - Christopher J Ciccarino
- John A. Paulson School of Engineering and Applied Sciences, Harvard University, Cambridge, Massachusetts 02138, United States.,Department of Chemistry and Chemical Biology, Harvard University, Cambridge, Massachusetts 02138, United States
| | - Prineha Narang
- John A. Paulson School of Engineering and Applied Sciences, Harvard University, Cambridge, Massachusetts 02138, United States
| |
Collapse
|
23
|
Mannouch JR, Richardson JO. A partially linearized spin-mapping approach for nonadiabatic dynamics. I. Derivation of the theory. J Chem Phys 2020; 153:194109. [DOI: 10.1063/5.0031168] [Citation(s) in RCA: 13] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
24
|
Mannouch JR, Richardson JO. A partially linearized spin-mapping approach for nonadiabatic dynamics. II. Analysis and comparison with related approaches. J Chem Phys 2020; 153:194110. [DOI: 10.1063/5.0031173] [Citation(s) in RCA: 9] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
|
25
|
Kundu S, Makri N. Modular path integral for finite-temperature dynamics of extended systems with intramolecular vibrations. J Chem Phys 2020; 153:044124. [DOI: 10.1063/5.0014838] [Citation(s) in RCA: 12] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/19/2023] Open
Affiliation(s)
- Sohang Kundu
- Department of Chemistry, University of Illinois, Urbana, Illinois 61801, USA
| | - Nancy Makri
- Department of Chemistry, University of Illinois, Urbana, Illinois 61801, USA
- Department of Physics, University of Illinois, Urbana, Illinois 61801, USA
| |
Collapse
|
26
|
Cohen G, Galperin M. Green’s function methods for single molecule junctions. J Chem Phys 2020; 152:090901. [DOI: 10.1063/1.5145210] [Citation(s) in RCA: 25] [Impact Index Per Article: 6.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/12/2023] Open
Affiliation(s)
- Guy Cohen
- The Raymond and Beverley Sackler Center for Computational Molecular and Materials Science, Tel Aviv University, Tel Aviv 69978, Israel
- School of Chemistry, Tel Aviv University, Tel Aviv 69978, Israel
| | - Michael Galperin
- Department of Chemistry and Biochemistry, University of California San Diego, La Jolla, California 92093, USA
| |
Collapse
|
27
|
Runeson JE, Richardson JO. Generalized spin mapping for quantum-classical dynamics. J Chem Phys 2020; 152:084110. [DOI: 10.1063/1.5143412] [Citation(s) in RCA: 27] [Impact Index Per Article: 6.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/24/2023] Open
Affiliation(s)
- Johan E. Runeson
- Laboratory of Physical Chemistry, ETH Zürich, 8093 Zürich, Switzerland
| | | |
Collapse
|
28
|
Chatterjee S, Makri N. Real-Time Path Integral Methods, Quantum Master Equations, and Classical vs Quantum Memory. J Phys Chem B 2019; 123:10470-10482. [DOI: 10.1021/acs.jpcb.9b08429] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Sambarta Chatterjee
- Department of Chemistry, University of Illinois, 505 South Mathews Avenue, Urbana, Illinois 61801, United States
| | - Nancy Makri
- Department of Chemistry, University of Illinois, 505 South Mathews Avenue, Urbana, Illinois 61801, United States
| |
Collapse
|
29
|
Mulvihill E, Gao X, Liu Y, Schubert A, Dunietz BD, Geva E. Combining the mapping Hamiltonian linearized semiclassical approach with the generalized quantum master equation to simulate electronically nonadiabatic molecular dynamics. J Chem Phys 2019; 151:074103. [DOI: 10.1063/1.5110891] [Citation(s) in RCA: 22] [Impact Index Per Article: 4.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Ellen Mulvihill
- Department of Chemistry, University of Michigan, Ann Arbor, Michigan 48109, USA
| | - Xing Gao
- Department of Chemistry, University of Michigan, Ann Arbor, Michigan 48109, USA
| | - Yudan Liu
- Department of Chemistry, University of Michigan, Ann Arbor, Michigan 48109, USA
| | - Alexander Schubert
- Department of Chemistry, University of Michigan, Ann Arbor, Michigan 48109, USA
- Department of Chemistry and Biochemistry, Kent State University, Kent, Ohio 44242, USA
| | - Barry D. Dunietz
- Department of Chemistry and Biochemistry, Kent State University, Kent, Ohio 44242, USA
| | - Eitan Geva
- Department of Chemistry, University of Michigan, Ann Arbor, Michigan 48109, USA
| |
Collapse
|
30
|
Runeson JE, Richardson JO. Spin-mapping approach for nonadiabatic molecular dynamics. J Chem Phys 2019; 151:044119. [DOI: 10.1063/1.5100506] [Citation(s) in RCA: 33] [Impact Index Per Article: 6.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Johan E. Runeson
- Laboratory of Physical Chemistry, ETH Zürich, 8093 Zürich, Switzerland
| | | |
Collapse
|
31
|
Liu J, Segal D. Interplay of Direct and Indirect Charge-Transfer Pathways in Donor–Bridge–Acceptor Systems. J Phys Chem B 2019; 123:6099-6110. [DOI: 10.1021/acs.jpcb.9b04958] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Junjie Liu
- Department of Chemistry and Centre for Quantum Information and Quantum Control, University of Toronto, Toronto, Ontario M5S 3H6, Canada
| | - Dvira Segal
- Department of Chemistry and Centre for Quantum Information and Quantum Control, University of Toronto, Toronto, Ontario M5S 3H6, Canada
| |
Collapse
|
32
|
Pfalzgraff WC, Montoya-Castillo A, Kelly A, Markland TE. Efficient construction of generalized master equation memory kernels for multi-state systems from nonadiabatic quantum-classical dynamics. J Chem Phys 2019; 150:244109. [DOI: 10.1063/1.5095715] [Citation(s) in RCA: 21] [Impact Index Per Article: 4.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- William C. Pfalzgraff
- Department of Chemistry, Stanford University, Stanford, California 94305, USA
- Department of Chemistry, Chatham University, Pittsburgh, Pennsylvania 15232, USA
| | | | - Aaron Kelly
- Department of Chemistry, Dalhousie University, Halifax, Nova Scotia B3H 4R2, Canada
| | - Thomas E. Markland
- Department of Chemistry, Stanford University, Stanford, California 94305, USA
| |
Collapse
|
33
|
Kelly A. Mean field theory of thermal energy transport in molecular junctions. J Chem Phys 2019; 150:204107. [DOI: 10.1063/1.5089885] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Aaron Kelly
- Department of Chemistry, Dalhousie University, Halifax, Nova Scotia B3H 4R2, Canada
| |
Collapse
|
34
|
Polimeno A, Zerbetto M, Abergel D. Stochastic modeling of macromolecules in solution. II. Spectral densities. J Chem Phys 2019; 150:184108. [PMID: 31091922 DOI: 10.1063/1.5077066] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/10/2023] Open
Abstract
In Paper I [Polimeno et al., J. Chem. Phys. 150, 184107 (2019)], we proposed a general approach for interpreting relaxation properties of a macromolecule in solution, derived from an atomistic description. A simple scheme (the semiflexible Brownian, SFB, model) has been defined for the case of limited internal flexibility, but retaining full coupling with external degrees of freedom, inclusion of all of the momenta, and dissipation. Here we discuss the application of the SFB model to the practical evaluation of orientation spectral densities, based on two complementary computational treatments.
Collapse
Affiliation(s)
- Antonino Polimeno
- Dipartimento di Scienze Chimiche, Università degli Studi di Padova, I-35131 Padova, Italy
| | - Mirco Zerbetto
- Dipartimento di Scienze Chimiche, Università degli Studi di Padova, I-35131 Padova, Italy
| | - Daniel Abergel
- Laboratoire des Biomolécules, LBM, Département de Chimie, Ecole Normale Supérieure, PSL University, Sorbonne Université, CNRS, 75005 Paris, France
| |
Collapse
|
35
|
Saller MAC, Kelly A, Richardson JO. On the identity of the identity operator in nonadiabatic linearized semiclassical dynamics. J Chem Phys 2019; 150:071101. [DOI: 10.1063/1.5082596] [Citation(s) in RCA: 37] [Impact Index Per Article: 7.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
| | - Aaron Kelly
- Department of Chemistry, Dalhousie University, Halifax, Nova Scotia B3H 4R2, Canada
| | | |
Collapse
|
36
|
Mulvihill E, Schubert A, Sun X, Dunietz BD, Geva E. A modified approach for simulating electronically nonadiabatic dynamics via the generalized quantum master equation. J Chem Phys 2019; 150:034101. [DOI: 10.1063/1.5055756] [Citation(s) in RCA: 26] [Impact Index Per Article: 5.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Ellen Mulvihill
- Department of Chemistry, University of Michigan, Ann Arbor, Michigan 48109, USA
| | - Alexander Schubert
- Department of Chemistry, University of Michigan, Ann Arbor, Michigan 48109, USA
- Department of Chemistry and Biochemistry, Kent State University, Kent, Ohio 44242, USA
| | - Xiang Sun
- Department of Chemistry, University of Michigan, Ann Arbor, Michigan 48109, USA
| | - Barry D. Dunietz
- Department of Chemistry and Biochemistry, Kent State University, Kent, Ohio 44242, USA
| | - Eitan Geva
- Department of Chemistry, University of Michigan, Ann Arbor, Michigan 48109, USA
| |
Collapse
|
37
|
Kidon L, Wang H, Thoss M, Rabani E. On the memory kernel and the reduced system propagator. J Chem Phys 2018; 149:104105. [DOI: 10.1063/1.5047446] [Citation(s) in RCA: 24] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Lyran Kidon
- Department of Chemistry, University of California, Berkeley, California 94720, USA
- Materials Science Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, USA
| | - Haobin Wang
- Department of Chemistry, University of Colorado Denver, Denver, Colorado 80217-3364, USA
| | - Michael Thoss
- Institute of Physics, University of Freiburg, Hermann-Herder-Strasse 3, 79104 Freiburg, Germany
| | - Eran Rabani
- Department of Chemistry, University of California, Berkeley, California 94720, USA
- Materials Science Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, USA
- The Sackler Center for Computational Molecular and Materials Science, Tel Aviv University, Tel Aviv 69978, Israel
| |
Collapse
|
38
|
Xu M, Yan Y, Liu Y, Shi Q. Convergence of high order memory kernels in the Nakajima-Zwanzig generalized master equation and rate constants: Case study of the spin-boson model. J Chem Phys 2018; 148:164101. [DOI: 10.1063/1.5022761] [Citation(s) in RCA: 12] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Meng Xu
- Beijing National Laboratory for Molecular Sciences, State Key Laboratory for Structural Chemistry of Unstable and Stable Species, CAS Research/Education Center for Excellence in Molecular Sciences, Institute of Chemistry, Chinese Academy of Sciences, Zhongguancun, Beijing 100190, China and University of Chinese Academy of Sciences, Beijing 100049, China
| | - Yaming Yan
- Beijing National Laboratory for Molecular Sciences, State Key Laboratory for Structural Chemistry of Unstable and Stable Species, CAS Research/Education Center for Excellence in Molecular Sciences, Institute of Chemistry, Chinese Academy of Sciences, Zhongguancun, Beijing 100190, China and University of Chinese Academy of Sciences, Beijing 100049, China
| | - Yanying Liu
- Beijing National Laboratory for Molecular Sciences, State Key Laboratory for Structural Chemistry of Unstable and Stable Species, CAS Research/Education Center for Excellence in Molecular Sciences, Institute of Chemistry, Chinese Academy of Sciences, Zhongguancun, Beijing 100190, China and University of Chinese Academy of Sciences, Beijing 100049, China
| | - Qiang Shi
- Beijing National Laboratory for Molecular Sciences, State Key Laboratory for Structural Chemistry of Unstable and Stable Species, CAS Research/Education Center for Excellence in Molecular Sciences, Institute of Chemistry, Chinese Academy of Sciences, Zhongguancun, Beijing 100190, China and University of Chinese Academy of Sciences, Beijing 100049, China
| |
Collapse
|
39
|
Kananenka AA, Sun X, Schubert A, Dunietz BD, Geva E. A comparative study of different methods for calculating electronic transition rates. J Chem Phys 2018; 148:102304. [DOI: 10.1063/1.4989509] [Citation(s) in RCA: 17] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Alexei A. Kananenka
- Department of Chemistry, University of Michigan, Ann Arbor, Michigan 48109, USA
| | - Xiang Sun
- Department of Chemistry, University of Michigan, Ann Arbor, Michigan 48109, USA
| | - Alexander Schubert
- Department of Chemistry, University of Michigan, Ann Arbor, Michigan 48109, USA
- Department of Chemistry and Biochemistry, Kent State University, Kent, Ohio 44242, USA
| | - Barry D. Dunietz
- Department of Chemistry and Biochemistry, Kent State University, Kent, Ohio 44242, USA
| | - Eitan Geva
- Department of Chemistry, University of Michigan, Ann Arbor, Michigan 48109, USA
| |
Collapse
|
40
|
Montoya-Castillo A, Reichman DR. Approximate but accurate quantum dynamics from the Mori formalism. II. Equilibrium time correlation functions. J Chem Phys 2017; 146:084110. [DOI: 10.1063/1.4975388] [Citation(s) in RCA: 18] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
| | - David R. Reichman
- Department of Chemistry, Columbia University, New York, New York 10027, USA
| |
Collapse
|
41
|
Montoya-Castillo A, Reichman DR. Path integral approach to the Wigner representation of canonical density operators for discrete systems coupled to harmonic baths. J Chem Phys 2017; 146:024107. [DOI: 10.1063/1.4973646] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
| | - David R. Reichman
- Department of Chemistry, Columbia University, New York, New York 10027, USA
| |
Collapse
|
42
|
Kelly A, Montoya-Castillo A, Wang L, Markland TE. Generalized quantum master equations in and out of equilibrium: When can one win? J Chem Phys 2016; 144:184105. [DOI: 10.1063/1.4948612] [Citation(s) in RCA: 42] [Impact Index Per Article: 5.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/03/2023] Open
Affiliation(s)
- Aaron Kelly
- Department of Chemistry, Stanford University, Stanford, California 94305, USA
| | | | - Lu Wang
- Department of Chemistry and Chemical Biology, Rutgers University, Piscataway, New Jersey 08854, USA
| | - Thomas E. Markland
- Department of Chemistry, Stanford University, Stanford, California 94305, USA
| |
Collapse
|