1
|
Zhou L, Gao X, Shuai Z. A stochastic Schrödinger equation and matrix product state approach to carrier transport in organic semiconductors with nonlocal electron-phonon interaction. J Chem Phys 2024; 161:084118. [PMID: 39212211 DOI: 10.1063/5.0221143] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/30/2024] [Accepted: 08/12/2024] [Indexed: 09/04/2024] Open
Abstract
Evaluation of the charge transport property of organic semiconductors requires exact quantum dynamics simulation of large systems. We present a numerically nearly exact approach to investigate carrier transport dynamics in organic semiconductors by extending the non-Markovian stochastic Schrödinger equation with complex frequency modes to a forward-backward scheme and by solving it using the matrix product state (MPS) approach. By utilizing the forward-backward formalism for noise generation, the bath correlation function can be effectively treated as a temperature-independent imaginary part, enabling a more accurate decomposition with fewer complex frequency modes. Using this approach, we study the carrier transport and mobility in the one-dimensional Peierls model, where the nonlocal electron-phonon interaction is taken into account. The reliability of this approach was validated by comparing carrier diffusion motion with those obtained from the hierarchical equations of motion method across various parameter regimes of the phonon bath. The efficiency was demonstrated by the modest virtual bond dimensions of MPS and the low scaling of the computational time with the system size.
Collapse
Affiliation(s)
- Liqi Zhou
- MOE Key Laboratory of Organic OptoElectronics and Molecular Engineering, Department of Chemistry, Tsinghua University, Beijing 100084, China
| | - Xing Gao
- School of Materials, Sun Yat-sen University, Shenzhen, Guangdong 518107, China
| | - Zhigang Shuai
- MOE Key Laboratory of Organic OptoElectronics and Molecular Engineering, Department of Chemistry, Tsinghua University, Beijing 100084, China
- School of Science and Engineering, The Chinese University of Hong Kong, Shenzhen, Guangdong 518172, China
| |
Collapse
|
2
|
Takahashi H, Borrelli R. Tensor-Train Format Hierarchical Equations of Motion Formalism: Charge Transfer in Organic Semiconductors via Dissipative Holstein Models. J Chem Theory Comput 2024. [PMID: 39152908 DOI: 10.1021/acs.jctc.4c00711] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 08/19/2024]
Abstract
Hierarchical Equations of Motion (HEOM) in the Tensor-Train (TT) representation is applied to study the charge-transfer dynamics in organic semiconductors (OSCs). The theoretical formulation as well as the basic computational aspects of HEOM-TT are discussed in detail. Charge transfer in OSCs is modeled using dissipative polaronic models that incorporate the effects of both high- and low-frequency molecular vibrations, and it is simulated in a fully quantum and nonperturbative manner, which has not been studied intensively. The capability of treating complex electron-vibrational systems is examined by analyzing and comparing the numerical behavior of the time-dependent variational approach and the time-Alternating Minimal Energy methods and by calculating the current autocorrelation function and diffusivity across various models. Our results indicate that the HEOM-TT framework offers a robust tool for the detailed analysis of complex polaronic systems, suggesting its potential for broader applications.
Collapse
|
3
|
Le Dé B, Jaouadi A, Mangaud E, Chin AW, Desouter-Lecomte M. Managing temperature in open quantum systems strongly coupled with structured environments. J Chem Phys 2024; 160:244102. [PMID: 38913841 DOI: 10.1063/5.0214051] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/15/2024] [Accepted: 06/06/2024] [Indexed: 06/26/2024] Open
Abstract
In non-perturbative non-Markovian open quantum systems, reaching either low temperatures with the hierarchical equations of motion (HEOM) or high temperatures with the Thermalized Time Evolving Density Operator with Orthogonal Polynomials Algorithm (T-TEDOPA) formalism in Hilbert space remains challenging. We compare different ways of modeling the environment. Sampling the Fourier transform of the bath correlation function, also called temperature dependent spectral density, proves to be very effective. T-TEDOPA [Tamascelli et al., Phys. Rev. Lett. 123, 090402 (2019)] uses a linear chain of oscillators with positive and negative frequencies, while HEOM is based on the complex poles of an optimized rational decomposition of the temperature dependent spectral density [Xu et al., Phys. Rev. Lett. 129, 230601 (2022)]. Resorting to the poles of the temperature independent spectral density and of the Bose function separately is an alternative when the problem due to the huge number of Bose poles at low temperatures is circumvented. Two examples illustrate the effectiveness of the HEOM and T-TEDOPA approaches: a benchmark pure dephasing case and a two-bath model simulating the dynamics of excited electronic states coupled through a conical intersection. We show the efficiency of T-TEDOPA to simulate dynamics at a finite temperature by using either continuous spectral densities or only all the intramolecular oscillators of a linear vibronic model calibrated from ab initio data of a phenylene ethynylene dimer.
Collapse
Affiliation(s)
- Brieuc Le Dé
- Institut des Nanosciences de Paris, Sorbonne Université, CNRS, F-75005 Paris, France
| | - Amine Jaouadi
- LyRIDS, ECE Paris, Graduate School of Engineering, Paris F-75015, France
| | - Etienne Mangaud
- MSME, Université Gustave Eiffel, UPEC, CNRS, F-77454 Marne-La-Vallée, France
| | - Alex W Chin
- Institut des Nanosciences de Paris, Sorbonne Université, CNRS, F-75005 Paris, France
| | | |
Collapse
|
4
|
Chen X, Franco I. Bexcitonics: Quasiparticle approach to open quantum dynamics. J Chem Phys 2024; 160:204116. [PMID: 38814013 DOI: 10.1063/5.0198567] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/18/2024] [Accepted: 05/09/2024] [Indexed: 05/31/2024] Open
Abstract
We develop a quasiparticle approach to capture the dynamics of open quantum systems coupled to bosonic thermal baths of arbitrary complexity based on the Hierarchical Equations of Motion (HEOM). This is done by generalizing the HEOM dynamics and mapping it into that of the system in interaction with a few bosonic fictitious quasiparticles that we call bexcitons. Bexcitons arise from a decomposition of the bath correlation function into discrete features. Specifically, bexciton creation and annihilation couple the auxiliary density matrices in the HEOM. The approach provides a systematic strategy to construct exact quantum master equations that include the system-bath coupling to all orders even for non-Markovian environments. Specifically, by introducing different metrics and representations for the bexcitons it is possible to straightforwardly generate different variants of the HEOM, demonstrating that all these variants share a common underlying quasiparticle picture. Bexcitonic properties, while unphysical, offer a coarse-grained view of the correlated system-bath dynamics and its numerical convergence. For instance, we use it to analyze the instability of the HEOM when the bath is composed of underdamped oscillators and show that it leads to the creation of highly excited bexcitons. The bexcitonic picture can also be used to develop more efficient approaches to propagate the HEOM. As an example, we use the particle-like nature of the bexcitons to introduce mode-combination of bexcitons in both number and coordinate representation that uses the multi-configuration time-dependent Hartree to efficiently propagate the HEOM dynamics.
Collapse
Affiliation(s)
- Xinxian Chen
- Department of Chemistry, University of Rochester, Rochester, New York 14627, USA
| | - Ignacio Franco
- Department of Chemistry, University of Rochester, Rochester, New York 14627, USA
- Department of Physics, University of Rochester, Rochester, New York 14627, USA
| |
Collapse
|
5
|
Kang M, Nuomin H, Chowdhury SN, Yuly JL, Sun K, Whitlow J, Valdiviezo J, Zhang Z, Zhang P, Beratan DN, Brown KR. Seeking a quantum advantage with trapped-ion quantum simulations of condensed-phase chemical dynamics. Nat Rev Chem 2024; 8:340-358. [PMID: 38641733 DOI: 10.1038/s41570-024-00595-1] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Accepted: 03/07/2024] [Indexed: 04/21/2024]
Abstract
Simulating the quantum dynamics of molecules in the condensed phase represents a longstanding challenge in chemistry. Trapped-ion quantum systems may serve as a platform for the analog-quantum simulation of chemical dynamics that is beyond the reach of current classical-digital simulation. To identify a 'quantum advantage' for these simulations, performance analysis of both analog-quantum simulation on noisy hardware and classical-digital algorithms is needed. In this Review, we make a comparison between a noisy analog trapped-ion simulator and a few choice classical-digital methods on simulating the dynamics of a model molecular Hamiltonian with linear vibronic coupling. We describe several simple Hamiltonians that are commonly used to model molecular systems, which can be simulated with existing or emerging trapped-ion hardware. These Hamiltonians may serve as stepping stones towards the use of trapped-ion simulators for systems beyond the reach of classical-digital methods. Finally, we identify dynamical regimes in which classical-digital simulations seem to have the weakest performance with respect to analog-quantum simulations. These regimes may provide the lowest hanging fruit to make the most of potential quantum advantages.
Collapse
Affiliation(s)
- Mingyu Kang
- Duke Quantum Center, Duke University, Durham, NC, USA.
- Department of Physics, Duke University, Durham, NC, USA.
| | - Hanggai Nuomin
- Department of Chemistry, Duke University, Durham, NC, USA
| | | | - Jonathon L Yuly
- Lewis-Sigler Institute for Integrative Genomics, Princeton University, Princeton, NJ, USA
| | - Ke Sun
- Duke Quantum Center, Duke University, Durham, NC, USA
- Department of Physics, Duke University, Durham, NC, USA
| | - Jacob Whitlow
- Duke Quantum Center, Duke University, Durham, NC, USA
- Department of Electrical and Computer Engineering, Duke University, Durham, NC, USA
| | - Jesús Valdiviezo
- Kenneth S. Pitzer Theory Center, University of California, Berkeley, CA, USA
- Department of Chemistry, University of California, Berkeley, CA, USA
- Departamento de Ciencias, Sección Química, Pontificia Universidad Católica del Perú, Lima, Peru
| | - Zhendian Zhang
- Department of Chemistry, Duke University, Durham, NC, USA
| | - Peng Zhang
- Department of Chemistry, Duke University, Durham, NC, USA
| | - David N Beratan
- Department of Physics, Duke University, Durham, NC, USA.
- Department of Chemistry, Duke University, Durham, NC, USA.
- Department of Biochemistry, Duke University, Durham, NC, USA.
| | - Kenneth R Brown
- Duke Quantum Center, Duke University, Durham, NC, USA.
- Department of Physics, Duke University, Durham, NC, USA.
- Department of Chemistry, Duke University, Durham, NC, USA.
- Department of Electrical and Computer Engineering, Duke University, Durham, NC, USA.
| |
Collapse
|
6
|
Liu Z, Hu H, Sun X. Multistate Reaction Coordinate Model for Charge and Energy Transfer Dynamics in the Condensed Phase. J Chem Theory Comput 2023; 19:7151-7170. [PMID: 37815937 PMCID: PMC10601487 DOI: 10.1021/acs.jctc.3c00770] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/13/2023] [Indexed: 10/12/2023]
Abstract
Constructing multistate model Hamiltonians from all-atom electronic structure calculations and molecular dynamics simulations is crucial for understanding charge and energy transfer dynamics in complex condensed phases. The most popular two-level system model is the spin-boson Hamiltonian, where the nuclear degrees of freedom are represented as shifted normal modes. Recently, we proposed the general multistate nontrivial extension of the spin-boson model, i.e., the multistate harmonic (MSH) model, which is constructed by extending the spatial dimensions of each nuclear mode so as to satisfy the all-atom reorganization energy restrictions for all pairs of electronic states. In this work, we propose the multistate reaction coordinate (MRC) model with a primary reaction coordinate and secondary bath modes as in the Caldeira-Leggett form but in extended spatial dimensions. The MRC model is proven to be equivalent to the MSH model and offers an intuitive physical picture of the nuclear-electronic feedback in nonadiabatic processes such as the inherent trajectory of the reaction coordinate. The reaction coordinate is represented in extended dimensions, carrying the entire reorganization energies and bilinearly coupled to the secondary bath modes. We demonstrate the MRC model construction for photoinduced charge transfer in an organic photovoltaic caroteniod-porphyrin-C60 molecular triad dissolved in tetrahydrofuran as well as excitation energy transfer in a photosynthetic light-harvesting Fenna-Matthews-Olson complex. The MRC model provides an effective and robust platform for investigating quantum dissipative dynamics in complex condensed-phase systems since it allows a consistent description of realistic spectral density, state-dependent system-bath couplings, and heterogeneous environments due to static disorder in reorganization energies.
Collapse
Affiliation(s)
- Zengkui Liu
- Division
of Arts and Sciences, NYU Shanghai, 567 West Yangsi Road, Shanghai, 200124, China
- NYU-ECNU
Center for Computational Chemistry at NYU Shanghai, 3663 Zhongshan Road North, Shanghai, 200062, China
- Department
of Chemistry, New York University, New York, New York, 10003, United States
| | - Haorui Hu
- Division
of Arts and Sciences, NYU Shanghai, 567 West Yangsi Road, Shanghai, 200124, China
| | - Xiang Sun
- Division
of Arts and Sciences, NYU Shanghai, 567 West Yangsi Road, Shanghai, 200124, China
- NYU-ECNU
Center for Computational Chemistry at NYU Shanghai, 3663 Zhongshan Road North, Shanghai, 200062, China
- Department
of Chemistry, New York University, New York, New York, 10003, United States
- Shanghai
Frontiers Science Center of Artificial Intelligence and Deep Learning, NYU Shanghai, 567 West Yangsi Road, Shanghai, 200124, China
| |
Collapse
|
7
|
Gelin MF, Borrelli R. Thermo-Field Dynamics Approach to Photo-induced Electronic Transitions Driven by Incoherent Thermal Radiation. J Chem Theory Comput 2023; 19:6402-6413. [PMID: 37656914 DOI: 10.1021/acs.jctc.3c00590] [Citation(s) in RCA: 4] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 09/03/2023]
Abstract
The effects of thermal light-matter interaction on the dynamics of photo-induced electronic transitions in molecules are investigated using a novel first principles approach based on the thermo-field dynamics description of both the molecular vibrational modes and of the radiation field. The developed approach permits numerically accurate simulations of quantum dynamics of electronic/excitonic systems coupled to nuclear and photonic baths kept at different temperatures. The baths can be described by arbitrary spectral densities and can have any system-bath coupling strengths. In agreement with the results obtained previously by less rigorous methods, we show that the excitation process obtained by the continuous interaction with the suddenly turned-on thermal radiation field creates a mixed ensemble having a nonnegligible component consisting of a superposition of vibronic eigenstates which can sustain coherent oscillations for relatively long times. The results become especially relevant for the dynamics of electronic transitions upon sunlight excitation. Analytical results based on time-dependent perturbation theory support the numerical simulations and provide a simple interpretation of the time evolution of quantum observables.
Collapse
Affiliation(s)
- Maxim F Gelin
- School of Sciences, Hangzhou Dianzi University, Hangzhou 310018, China
| | - Raffaele Borrelli
- DISAFA, University of Torino, Largo Paolo Braccini 2, Grugliasco I-10095, Italy
| |
Collapse
|
8
|
Chen ZH, Wang Y, Xu RX, Yan Y. Open quantum systems with nonlinear environmental backactions: Extended dissipaton theory vs core-system hierarchy construction. J Chem Phys 2023; 158:074102. [PMID: 36813728 DOI: 10.1063/5.0134700] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/27/2023] Open
Abstract
In this paper, we present a comprehensive account of quantum dissipation theories with the quadratic environment couplings. The theoretical development includes the Brownian solvation mode embedded hierarchical quantum master equations, a core-system hierarchy construction that verifies the extended dissipaton equation of motion (DEOM) formalism [R. X. Xu et al., J. Chem. Phys. 148, 114103 (2018)]. Developed are also the quadratic imaginary-time DEOM for equilibrium and the λ(t)-DEOM for nonequilibrium thermodynamics problems. Both the celebrated Jarzynski equality and Crooks relation are accurately reproduced, which, in turn, confirms the rigorousness of the extended DEOM theories. While the extended DEOM is more numerically efficient, the core-system hierarchy quantum master equation is favorable for "visualizing" the correlated solvation dynamics.
Collapse
Affiliation(s)
- Zi-Hao Chen
- Hefei National Research Center for Physical Sciences at the Microscale and Department of Chemical Physics, University of Science and Technology of China, Hefei, Anhui 230026, China
| | - Yao Wang
- Hefei National Research Center for Physical Sciences at the Microscale and Department of Chemical Physics, University of Science and Technology of China, Hefei, Anhui 230026, China
| | - Rui-Xue Xu
- Hefei National Research Center for Physical Sciences at the Microscale and Department of Chemical Physics, University of Science and Technology of China, Hefei, Anhui 230026, China
| | - YiJing Yan
- Hefei National Research Center for Physical Sciences at the Microscale and Department of Chemical Physics, University of Science and Technology of China, Hefei, Anhui 230026, China
| |
Collapse
|
9
|
Zhang D, Zuo L, Ye L, Chen ZH, Wang Y, Xu RX, Zheng X, Yan Y. Hierarchical equations of motion approach for accurate characterization of spin excitations in quantum impurity systems. J Chem Phys 2023; 158:014106. [PMID: 36610957 DOI: 10.1063/5.0131739] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/23/2022] Open
Abstract
Recent technological advancement in scanning tunneling microscopes has enabled the measurement of spin-field and spin-spin interactions in single atomic or molecular junctions with an unprecedentedly high resolution. Theoretically, although the fermionic hierarchical equations of motion (HEOM) method has been widely applied to investigate the strongly correlated Kondo states in these junctions, the existence of low-energy spin excitations presents new challenges to numerical simulations. These include the quest for a more accurate and efficient decomposition for the non-Markovian memory of low-temperature environments and a more careful handling of errors caused by the truncation of the hierarchy. In this work, we propose several new algorithms, which significantly enhance the performance of the HEOM method, as exemplified by the calculations on systems involving various types of low-energy spin excitations. Being able to characterize both the Kondo effect and spin excitation accurately, the HEOM method offers a sophisticated and versatile theoretical tool, which is valuable for the understanding and even prediction of the fascinating quantum phenomena explored in cutting-edge experiments.
Collapse
Affiliation(s)
- Daochi Zhang
- Hefei National Research Center for Physical Sciences at the Microscale and Synergetic Innovation Center of Quantum Information and Quantum Physics, University of Science and Technology of China, Hefei, Anhui 230026, China
| | - Lijun Zuo
- Hefei National Research Center for Physical Sciences at the Microscale and Synergetic Innovation Center of Quantum Information and Quantum Physics, University of Science and Technology of China, Hefei, Anhui 230026, China
| | - Lyuzhou Ye
- Hefei National Research Center for Physical Sciences at the Microscale and Synergetic Innovation Center of Quantum Information and Quantum Physics, University of Science and Technology of China, Hefei, Anhui 230026, China
| | - Zi-Hao Chen
- Hefei National Research Center for Physical Sciences at the Microscale and Synergetic Innovation Center of Quantum Information and Quantum Physics, University of Science and Technology of China, Hefei, Anhui 230026, China
| | - Yao Wang
- Hefei National Research Center for Physical Sciences at the Microscale and Synergetic Innovation Center of Quantum Information and Quantum Physics, University of Science and Technology of China, Hefei, Anhui 230026, China
| | - Rui-Xue Xu
- Hefei National Research Center for Physical Sciences at the Microscale and Synergetic Innovation Center of Quantum Information and Quantum Physics, University of Science and Technology of China, Hefei, Anhui 230026, China
| | - Xiao Zheng
- Hefei National Research Center for Physical Sciences at the Microscale and Synergetic Innovation Center of Quantum Information and Quantum Physics, University of Science and Technology of China, Hefei, Anhui 230026, China
| | - YiJing Yan
- Hefei National Research Center for Physical Sciences at the Microscale and iChEM, University of Science and Technology of China, Hefei, Anhui 230026, China
| |
Collapse
|
10
|
Lin K, Peng J, Xu C, Gu FL, Lan Z. Automatic Evolution of Machine-Learning-Based Quantum Dynamics with Uncertainty Analysis. J Chem Theory Comput 2022; 18:5837-5855. [DOI: 10.1021/acs.jctc.2c00702] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/06/2022]
Affiliation(s)
- Kunni Lin
- School of Chemistry, South China Normal University, Guangzhou510006, P. R. China
- MOE Key Laboratory of Environmental Theoretical Chemistry, South China Normal University, Guangzhou510006, P. R. China
| | - Jiawei Peng
- School of Chemistry, South China Normal University, Guangzhou510006, P. R. China
- MOE Key Laboratory of Environmental Theoretical Chemistry, South China Normal University, Guangzhou510006, P. R. China
| | - Chao Xu
- MOE Key Laboratory of Environmental Theoretical Chemistry, South China Normal University, Guangzhou510006, P. R. China
- SCNU Environmental Research Institute, Guangdong Provincial Key Laboratory of Chemical Pollution and Environmental Safety, School of Environment, South China Normal University, Guangzhou510006, P. R. China
| | - Feng Long Gu
- MOE Key Laboratory of Environmental Theoretical Chemistry, South China Normal University, Guangzhou510006, P. R. China
- SCNU Environmental Research Institute, Guangdong Provincial Key Laboratory of Chemical Pollution and Environmental Safety, School of Environment, South China Normal University, Guangzhou510006, P. R. China
| | - Zhenggang Lan
- MOE Key Laboratory of Environmental Theoretical Chemistry, South China Normal University, Guangzhou510006, P. R. China
- SCNU Environmental Research Institute, Guangdong Provincial Key Laboratory of Chemical Pollution and Environmental Safety, School of Environment, South China Normal University, Guangzhou510006, P. R. China
| |
Collapse
|
11
|
Mitrić P, Janković V, Vukmirović N, Tanasković D. Spectral Functions of the Holstein Polaron: Exact and Approximate Solutions. PHYSICAL REVIEW LETTERS 2022; 129:096401. [PMID: 36083668 DOI: 10.1103/physrevlett.129.096401] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 01/10/2022] [Revised: 05/02/2022] [Accepted: 08/05/2022] [Indexed: 06/15/2023]
Abstract
It is generally accepted that the dynamical mean field theory gives a good solution of the Holstein model, but only in dimensions greater than two. Here, we show that this theory, which becomes exact in the weak coupling and in the atomic limit, provides an excellent, numerically cheap, approximate solution for the spectral function of the Holstein model in the whole range of parameters, even in one dimension. To establish this, we make a detailed comparison with the spectral functions that we obtain using the newly developed momentum-space numerically exact hierarchical equations of motion method, which yields electronic correlation functions directly in real time. We crosscheck these conclusions with our path integral quantum Monte Carlo and exact diagonalization results, as well as with the available numerically exact results from the literature.
Collapse
Affiliation(s)
- Petar Mitrić
- Institute of Physics Belgrade, University of Belgrade, Pregrevica 118, 11080 Belgrade, Serbia
| | - Veljko Janković
- Institute of Physics Belgrade, University of Belgrade, Pregrevica 118, 11080 Belgrade, Serbia
| | - Nenad Vukmirović
- Institute of Physics Belgrade, University of Belgrade, Pregrevica 118, 11080 Belgrade, Serbia
| | - Darko Tanasković
- Institute of Physics Belgrade, University of Belgrade, Pregrevica 118, 11080 Belgrade, Serbia
| |
Collapse
|
12
|
Cainelli M, Borrelli R, Tanimura Y. Effect of mixed Frenkel and charge transfer states in time-gated fluorescence spectra of perylene bisimides H-aggregates: Hierarchical equations of motion approach. J Chem Phys 2022; 157:084103. [DOI: 10.1063/5.0102000] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
We theoretically investigated the effect of mixed Frenkel (F) and charge transfer (CT) states on the spectral properties of perylene bisimide (PBI) derivatives, focusing on the role of strong electron-phonon interactions. The model consists of a four-level system described by the Holstein Hamiltonian coupled to independent local heat-baths on each site, described by Brownian spectral distribution functions. We employ the reduced hierarchical equations of motion (HEOM) approach to calculate the time evolution of the system and compare it to the pure F exciton cases. We compute the absorption and time-gated fluorescence (TGF) spectra for different exciton transfer integrals and F-CT band gap conditions. The coherence length of excitons ($N_{coh}$) is evaluated employing two different definitions. We observe the presence of an excited hot state peak whose intensity is associated with the delocalization of the excited species and ultrafast dynamics that are solely dependent on the frequency of the local bath. The results indicate that the inclusion of CT states promotes localization of the excitons which is manifested in a decrease in the intensity of the hot state peak and the 0--1 peak, and an increase in the intensity of the 0--0 emission peak in TGF spectrum, leading to a decrease of $N_{coh}$.
Collapse
Affiliation(s)
| | - Raffaele Borrelli
- Department of Agricoltural Science, Università degli Studi di Torino, Italy
| | | |
Collapse
|
13
|
Bose A. Zero-Cost Corrections to Influence Functional Coefficients from Bath Response Functions. J Chem Phys 2022; 157:054107. [DOI: 10.1063/5.0101396] [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/15/2022] Open
Abstract
Recent work has shown that it is possible to circumvent the calculation of the spectral density and directly calculate the coefficients of the discretized influence functionals using data from classical trajectory simulations. However, the accuracy of this procedure depends on the validity of the high temperature approximation. In this work, an alternative derivation based on the Kubo formalism is provided. This enables the calculation of additional correction terms that increases the range of applicability of the procedure to lower temperatures. Because it is based on the Kubo-transformed correlation function, this approach enables the direct use of correlation functions obtained from methods like ring-polymer molecular dynamics and centroid molecular dynamics in determining the influence functional coefficients for subsequent system-solvent simulations. The accuracy of the original procedure and the corrected procedure is investigated across a range of parameters. It is interesting that the correction term comes at zero additional cost. Furthermore, it is possible to improve upon the correction using zero-cost physical intuition and heuristics making the method even more accurate.
Collapse
Affiliation(s)
- Amartya Bose
- Chemistry, Princeton University Department of Chemistry, United States of America
| |
Collapse
|
14
|
Xing T, Li T, Yan Y, Bai S, Shi Q. Application of the imaginary time hierarchical equations of motion method to calculate real time correlation functions. J Chem Phys 2022; 156:244102. [DOI: 10.1063/5.0095790] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/03/2023] Open
Abstract
We investigate the application of the imaginary time hierarchical equations of motion method to calculate real time quantum correlation functions. By starting from the path integral expression for the correlated system–bath equilibrium state, we first derive a new set of equations that decouple the imaginary time propagation and the calculation of auxiliary density operators. The new equations, thus, greatly simplify the calculation of the equilibrium correlated initial state that is subsequently used in the real time propagation to obtain the quantum correlation functions. It is also shown that a periodic decomposition of the bath imaginary time correlation function is no longer necessary in the new equations such that different decomposition schemes can be explored. The applicability of the new method is demonstrated in several numerical examples, including the spin-Boson model, the Holstein model, and the double-well model for proton transfer reaction.
Collapse
Affiliation(s)
- Tao Xing
- 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
| | - Tianchu Li
- 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
| | - Shuming Bai
- 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
|
15
|
Ke Y, Borrelli R, Thoss M. Hierarchical equations of motion approach to hybrid fermionic and bosonic environments: Matrix product state formulation in twin space. J Chem Phys 2022; 156:194102. [DOI: 10.1063/5.0088947] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
We extend the twin-space formulation of the hierarchical equations of motion approach in combination with the matrix product state representation (introduced in J. Chem. Phys. 150, 234102, [2019]) to nonequilibrium scenarios where the open quantum system is coupled to a hybrid fermionic and bosonic environment. The key ideas used in the extension are a reformulation of the hierarchical equations of motion for the auxiliary density matrices into a time-dependent Schrödinger-like equation for an augmented multi-dimensional wave function as well as a tensor decomposition into a product of low-rank matrices. The new approach facilitates accurate simulations of non-equilibrium quantum dynamics in larger and more complex open quantum systems. The performance of the method is demonstrated for a model of a molecular junction exhibiting current-induced mode-selective vibrational excitation.
Collapse
Affiliation(s)
- Yaling Ke
- Institute of Physics, Albert-Ludwigs-Universität Freiburg, Germany
| | - Raffaele Borrelli
- Department of Agricoltural Science, Università degli Studi di Torino, Italy
| | - Michael Thoss
- University of Freiburg Institute of Physics, Germany
| |
Collapse
|
16
|
Ren J, Li W, Jiang T, Wang Y, Shuai Z. Time‐dependent density matrix renormalization group method for quantum dynamics in complex systems. WIRES COMPUTATIONAL MOLECULAR SCIENCE 2022. [DOI: 10.1002/wcms.1614] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 01/28/2023]
Affiliation(s)
- Jiajun Ren
- MOE Key Laboratory of Organic OptoElectronics and Molecular Engineering, Department of Chemistry Tsinghua University Beijing People's Republic of China
| | - Weitang Li
- MOE Key Laboratory of Organic OptoElectronics and Molecular Engineering, Department of Chemistry Tsinghua University Beijing People's Republic of China
| | - Tong Jiang
- MOE Key Laboratory of Organic OptoElectronics and Molecular Engineering, Department of Chemistry Tsinghua University Beijing People's Republic of China
| | - Yuanheng Wang
- MOE Key Laboratory of Organic OptoElectronics and Molecular Engineering, Department of Chemistry Tsinghua University Beijing People's Republic of China
| | - Zhigang Shuai
- MOE Key Laboratory of Organic OptoElectronics and Molecular Engineering, Department of Chemistry Tsinghua University Beijing People's Republic of China
| |
Collapse
|
17
|
Ikeda T, Nakayama A. Collective bath coordinate mapping of "hierarchy" in hierarchical equations of motion. J Chem Phys 2022; 156:104104. [DOI: 10.1063/5.0082936] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
| | - Akira Nakayama
- Department of Chemical System Engineering, University of Tokyo, Japan
| |
Collapse
|
18
|
Li T, Yan Y, Shi Q. A low-temperature quantum Fokker-Planck equation that improvesthe numerical stability of the hierarchical equations of motion for the Brownian oscillator spectral density. J Chem Phys 2022; 156:064107. [DOI: 10.1063/5.0082108] [Citation(s) in RCA: 4] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Affiliation(s)
- Tianchu Li
- Institute of Chemistry Chinese Academy of Sciences, China
| | | | - Qiang Shi
- State Key Laboratory for Structural Chemistry of Unstable and Stable Species, Institute of Chemistry Chinese Academy of Sciences, China
| |
Collapse
|
19
|
Abstract
Surface hopping has seen great success in describing molecular phenomena where electronic excitations tend to be localized, but its application to materials with band-like electronic properties has remained limited. Here, we derive a formulation of fewest-switches surface hopping where both the quantum and classical equations of motion are solved entirely in terms of reciprocal-space coordinates. The resulting method is directly compatible with band structure calculations and allows for the efficient description of band-like phenomena by means of a truncation of the Brillouin zone. Using the Holstein and Peierls models as examples, we demonstrate the formal equivalence between real-space and reciprocal-space surface hopping and assess their accuracy against mean-field mixed quantum-classical dynamics and numerically exact results. Having very similar equations of motion, reciprocal-space surface hopping can be straightforwardly incorporated in existing (real-space) surface hopping implementations.
Collapse
Affiliation(s)
- Alex Krotz
- Department of Chemistry, Northwestern University, 2145 Sheridan Road, Evanston, Illinois 60208, USA
| | - Roel Tempelaar
- Department of Chemistry, Northwestern University, 2145 Sheridan Road, Evanston, Illinois 60208, USA
| |
Collapse
|
20
|
Bose A, Walters PL. A multisite decomposition of the tensor network path integrals. J Chem Phys 2022; 156:024101. [PMID: 35032978 DOI: 10.1063/5.0073234] [Citation(s) in RCA: 12] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Tensor network decompositions of path integrals for simulating open quantum systems have recently been proven to be useful. However, these methods scale exponentially with the system size. This makes it challenging to simulate the non-equilibrium dynamics of extended quantum systems coupled with local dissipative environments. In this work, we extend the tensor network path integral (TNPI) framework to efficiently simulate such extended systems. The Feynman-Vernon influence functional is a popular approach used to account for the effect of environments on the dynamics of the system. In order to facilitate the incorporation of the influence functional into a multisite framework (MS-TNPI), we combine a matrix product state (MPS) decomposition of the reduced density tensor of the system along the sites with a corresponding tensor network representation of the time axis to construct an efficient 2D tensor network. The 2D MS-TNPI network, when contracted, yields the time-dependent reduced density tensor of the extended system as an MPS. The algorithm presented is independent of the system Hamiltonian. We outline an iteration scheme to take the simulation beyond the non-Markovian memory introduced by solvents. Applications to spin chains coupled to local harmonic baths are presented; we consider the Ising, XXZ, and Heisenberg models, demonstrating that the presence of local environments can often dissipate the entanglement between the sites. We discuss three factors causing the system to transition from a coherent oscillatory dynamics to a fully incoherent dynamics. The MS-TNPI method is useful for studying a variety of extended quantum systems coupled with solvents.
Collapse
Affiliation(s)
- Amartya Bose
- Department of Chemistry, Princeton University, Princeton, New Jersey 08544, USA
| | - Peter L Walters
- Department of Chemistry, University of California, Berkeley, California 94720, USA
| |
Collapse
|
21
|
Zhang D, Ding X, Zhang HD, Zheng X, Yan Y. Adiabatic terminator for fermionic hierarchical equations of motion. CHINESE J CHEM PHYS 2021. [DOI: 10.1063/1674-0068/cjcp2110212] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022]
Affiliation(s)
- Daochi Zhang
- Hefei National Laboratory for Physical Sciences at the Microscale & Synergetic Innovation Center of Quantum Information and Quantum Physics & CAS Center for Excellence in Nanoscience, University of Science and Technology of China, Hefei 230026, China
| | - Xu Ding
- Hefei National Laboratory for Physical Sciences at the Microscale & Synergetic Innovation Center of Quantum Information and Quantum Physics & CAS Center for Excellence in Nanoscience, University of Science and Technology of China, Hefei 230026, China
| | - Hou-Dao Zhang
- Hefei National Laboratory for Physical Sciences at the Microscale & Synergetic Innovation Center of Quantum Information and Quantum Physics & CAS Center for Excellence in Nanoscience, University of Science and Technology of China, Hefei 230026, China
| | - Xiao Zheng
- Hefei National Laboratory for Physical Sciences at the Microscale & Synergetic Innovation Center of Quantum Information and Quantum Physics & CAS Center for Excellence in Nanoscience, University of Science and Technology of China, Hefei 230026, China
| | - YiJing Yan
- Hefei National Laboratory for Physical Sciences at the Microscale & iChEM, University of Science and Technology of China, Hefei 230026, China
| |
Collapse
|
22
|
Jang SJ, Burghardt I, Hsu CP, Bardeen CJ. Excitons: Energetics and spatiotemporal dynamics. J Chem Phys 2021; 155:200401. [PMID: 34852498 DOI: 10.1063/5.0075292] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Seogjoo J Jang
- Department of Chemistry and Biochemistry, Queens College, City University of New York, 65-30 Kissena Boulevard, Queens, New York 11367, USA and PhD Programs in Chemistry and Physics, and Initiative for the Theoretical Sciences, Graduate Center, City University of New York, 365 Fifth Avenue, New York, New York 10016, USA
| | - Irene Burghardt
- Institute of Physical and Theoretical Chemistry, Goethe University Frankfurt, Frankfurt, Germany
| | - Chao-Ping Hsu
- Institute of Chemistry, Academia Sinica, Taipei 115, Taiwan and Physics Division, National Center for Theoretical Sciences, Taipei 106, Taiwan
| | - Christopher J Bardeen
- Department of Chemistry, University of California, Riverside, Riverside, California 92521, USA
| |
Collapse
|
23
|
Gelin MF, Velardo A, Borrelli R. Efficient quantum dynamics simulations of complex molecular systems: A unified treatment of dynamic and static disorder. J Chem Phys 2021; 155:134102. [PMID: 34624969 DOI: 10.1063/5.0065896] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/05/2023] Open
Abstract
We present a unified and highly numerically efficient formalism for the simulation of quantum dynamics of complex molecular systems, which takes into account both temperature effects and static disorder. The methodology is based on the thermo-field dynamics formalism, and Gaussian static disorder is included into simulations via auxiliary bosonic operators. This approach, combined with the tensor-train/matrix-product state representation of the thermalized stochastic wave function, is applied to study the effect of dynamic and static disorders in charge-transfer processes in model organic semiconductor chains employing the Su-Schrieffer-Heeger (Holstein-Peierls) model Hamiltonian.
Collapse
Affiliation(s)
- Maxim F Gelin
- School of Sciences, Hangzhou Dianzi University, Hangzhou 310018, China
| | | | | |
Collapse
|
24
|
Hu Z, Brian D, Sun X. Multi-state harmonic models with globally shared bath for nonadiabatic dynamics in the condensed phase. J Chem Phys 2021; 155:124105. [PMID: 34598571 DOI: 10.1063/5.0064763] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/09/2023] Open
Abstract
Model Hamiltonians constructed from quantum chemistry calculations and molecular dynamics simulations are widely used for simulating nonadiabatic dynamics in the condensed phase. The most popular two-state spin-boson model could be built by mapping the all-atom anharmonic Hamiltonian onto a two-level system bilinearly coupled to a harmonic bath using the energy gap time correlation function. However, for more than two states, there lacks a general strategy to construct multi-state harmonic (MSH) models since the energy gaps between different pairs of electronic states are not entirely independent and need to be considered consistently. In this paper, we extend the previously proposed approach for building three-state harmonic models for photoinduced charge transfer to the arbitrary number of electronic states with a globally shared bath and the system-bath couplings are scaled differently according to the reorganization energies between each pair of states. We demonstrate the MSH model construction for an organic photovoltaic carotenoid-porphyrin-C60 molecular triad dissolved in explicit tetrahydrofuran solvent. Nonadiabatic dynamics was simulated using mixed quantum-classical techniques, including the linearized semiclassical and symmetrical quasiclassical dynamics with the mapping Hamiltonians, mean-field Ehrenfest, and mixed quantum-classical Liouville dynamics in two-state, three-state, and four-state harmonic models of the triad system. The MSH models are shown to provide a general and flexible framework for simulating nonadiabatic dynamics in complex systems.
Collapse
Affiliation(s)
- Zhubin Hu
- Division of Arts and Sciences, NYU Shanghai, 1555 Century Avenue, Shanghai 200122, China
| | - Dominikus Brian
- Division of Arts and Sciences, NYU Shanghai, 1555 Century Avenue, Shanghai 200122, China
| | - Xiang Sun
- Division of Arts and Sciences, NYU Shanghai, 1555 Century Avenue, Shanghai 200122, China
| |
Collapse
|
25
|
Gelin MF, Borrelli R. Simulation of Nonlinear Femtosecond Signals at Finite Temperature via a Thermo Field Dynamics-Tensor Train Method: General Theory and Application to Time- and Frequency-Resolved Fluorescence of the Fenna-Matthews-Olson Complex. J Chem Theory Comput 2021; 17:4316-4331. [PMID: 34076412 DOI: 10.1021/acs.jctc.1c00158] [Citation(s) in RCA: 15] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/01/2023]
Abstract
Addressing needs of contemporary nonlinear femtosecond optical spectroscopy, we have developed a fully quantum, numerically accurate wave function-based approach for the calculation of third-order spectroscopic signals of polyatomic molecules and molecular aggregates at finite temperature. The systems are described by multimode nonadiabatic vibronic-coupling Hamiltonians, in which diagonal terms are treated in harmonic approximation, while off-diagonal interstate couplings are assumed to be coordinate independent. The approach is based on the Thermo Field Dynamics (TFD) representation of quantum mechanics and tensor-train (TT) machinery for efficient numerical simulation of quantum evolution of systems with many degrees of freedom. The developed TFD-TT approach is applied to the calculation of time- and frequency-resolved fluorescence spectra of the Fenna-Matthews-Olson (FMO) antenna complex at room temperature taking into account finite time-frequency resolution in fluorescence detection, orientational averaging, and static disorder.
Collapse
Affiliation(s)
- Maxim F Gelin
- School of Sciences, Hangzhou Dianzi University, Hangzhou 310018, China
| | | |
Collapse
|
26
|
Krotz A, Provazza J, Tempelaar R. A reciprocal-space formulation of mixed quantum-classical dynamics. J Chem Phys 2021; 154:224101. [PMID: 34241207 DOI: 10.1063/5.0053177] [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/15/2022] Open
Abstract
We derive a formulation of mixed quantum-classical dynamics for modeling electronic carriers interacting with phonons in reciprocal space. For dispersionless phonons, we start by expressing the real-space classical coordinates in terms of complex variables. Taking these variables as a Fourier series then yields the reciprocal-space coordinates. Evaluating the electron-phonon interaction term through Ehrenfest's theorem, we arrive at a reciprocal-space formalism that is equivalent to mean-field mixed quantum-classical dynamics in real space. This equivalence is numerically verified for the Holstein and Peierls models, for which we find the reciprocal-space Hellmann-Feynman forces to involve momentum-derivative contributions in addition to the position-derivative terms commonly seen in real space. To illustrate the advantage of the reciprocal-space formulation, we present a proof of concept for the inexpensive modeling of low-momentum carriers interacting with phonons using a truncated reciprocal-space basis, which is not possible within a real-space formulation.
Collapse
Affiliation(s)
- Alex Krotz
- Department of Chemistry, Northwestern University, 2145 Sheridan Road, Evanston, Illinois 60208, USA
| | - Justin Provazza
- Department of Chemistry, Northwestern University, 2145 Sheridan Road, Evanston, Illinois 60208, USA
| | - Roel Tempelaar
- Department of Chemistry, Northwestern University, 2145 Sheridan Road, Evanston, Illinois 60208, USA
| |
Collapse
|
27
|
Yan Y, Xu M, Li T, Shi Q. Efficient propagation of the hierarchical equations of motion using the Tucker and hierarchical Tucker tensors. J Chem Phys 2021; 154:194104. [PMID: 34240893 DOI: 10.1063/5.0050720] [Citation(s) in RCA: 21] [Impact Index Per Article: 7.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
We develop new methods to efficiently propagate the hierarchical equations of motion (HEOM) by using the Tucker and hierarchical Tucker (HT) tensors to represent the reduced density operator and auxiliary density operators. We first show that by employing the split operator method, the specific structure of the HEOM allows a simple propagation scheme using the Tucker tensor. When the number of effective modes in the HEOM increases and the Tucker representation becomes intractable, the split operator method is extended to the binary tree structure of the HT representation. It is found that to update the binary tree nodes related to a specific effective mode, we only need to propagate a short matrix product state constructed from these nodes. Numerical results show that by further employing the mode combination technique commonly used in the multi-configuration time-dependent Hartree approaches, the binary tree representation can be applied to study excitation energy transfer dynamics in a fairly large system including over 104 effective modes. The new methods may thus provide a promising tool in simulating quantum dynamics in condensed phases.
Collapse
Affiliation(s)
- 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; University of Chinese Academy of Sciences, Beijing 100049, China; and Physical Science Laboratory, Huairou National Comprehensive Science Center, Beijing 101407, China
| | - 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; University of Chinese Academy of Sciences, Beijing 100049, China; and Physical Science Laboratory, Huairou National Comprehensive Science Center, Beijing 101407, China
| | - Tianchu Li
- 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; University of Chinese Academy of Sciences, Beijing 100049, China; and Physical Science Laboratory, Huairou National Comprehensive Science Center, Beijing 101407, 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; University of Chinese Academy of Sciences, Beijing 100049, China; and Physical Science Laboratory, Huairou National Comprehensive Science Center, Beijing 101407, China
| |
Collapse
|
28
|
Borrelli R, Dolgov S. Expanding the Range of Hierarchical Equations of Motion by Tensor-Train Implementation. J Phys Chem B 2021; 125:5397-5407. [DOI: 10.1021/acs.jpcb.1c02724] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Raffaele Borrelli
- DISAFA, University of Torino, Largo Paolo Braccini 2, Grugliasco 10095, Italy
| | - Sergey Dolgov
- University of Bath, Claverton Down, BA2 7AY Bath, United Kingdom
| |
Collapse
|
29
|
Borrelli R, Gelin MF. Finite temperature quantum dynamics of complex systems: Integrating
thermo‐field
theories and
tensor‐train
methods. WIRES COMPUTATIONAL MOLECULAR SCIENCE 2021. [DOI: 10.1002/wcms.1539] [Citation(s) in RCA: 18] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/18/2022]
Affiliation(s)
| | - Maxim F. Gelin
- School of Sciences Hangzhou Dianzi University Hangzhou China
| |
Collapse
|
30
|
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
|
31
|
Peng J, Xie Y, Hu D, Lan Z. Analysis of bath motion in MM-SQC dynamics via dimensionality reduction approach: Principal component analysis. J Chem Phys 2021; 154:094122. [PMID: 33685149 DOI: 10.1063/5.0039743] [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/14/2022] Open
Abstract
The system-plus-bath model is an important tool to understand the nonadiabatic dynamics of large molecular systems. Understanding the collective motion of a large number of bath modes is essential for revealing their key roles in the overall dynamics. Here, we applied principal component analysis (PCA) to investigate the bath motion in the basis of a large dataset generated from the symmetrical quasi-classical dynamics method based on the Meyer-Miller mapping Hamiltonian nonadiabatic dynamics for the excited-state energy transfer in the Frenkel-exciton model. The PCA method clearly elucidated that two types of bath modes, which either display strong vibronic coupling or have frequencies close to that of the electronic transition, are important to the nonadiabatic dynamics. These observations were fully consistent with the physical insights. The conclusions were based on the PCA of the trajectory data and did not involve significant pre-defined physical knowledge. The results show that the PCA approach, which is one of the simplest unsupervised machine learning dimensionality reduction methods, is a powerful one for analyzing complicated nonadiabatic dynamics in the condensed phase with many degrees of freedom.
Collapse
Affiliation(s)
- Jiawei Peng
- SCNU Environmental Research Institute, Guangdong Provincial Key Laboratory of Chemical Pollution and Environmental Safety & MOE Key Laboratory of Theoretical Chemistry of Environment, South China Normal University, Guangzhou 510006, China
| | - Yu Xie
- SCNU Environmental Research Institute, Guangdong Provincial Key Laboratory of Chemical Pollution and Environmental Safety & MOE Key Laboratory of Theoretical Chemistry of Environment, South China Normal University, Guangzhou 510006, China
| | - Deping Hu
- SCNU Environmental Research Institute, Guangdong Provincial Key Laboratory of Chemical Pollution and Environmental Safety & MOE Key Laboratory of Theoretical Chemistry of Environment, South China Normal University, Guangzhou 510006, China
| | - Zhenggang Lan
- SCNU Environmental Research Institute, Guangdong Provincial Key Laboratory of Chemical Pollution and Environmental Safety & MOE Key Laboratory of Theoretical Chemistry of Environment, South China Normal University, Guangzhou 510006, China
| |
Collapse
|