1
|
Citty B, Lynd JK, Gera T, Varvelo L, Raccah DIGB. MesoHOPS: Size-invariant scaling calculations of multi-excitation open quantum systems. J Chem Phys 2024; 160:144118. [PMID: 38619062 DOI: 10.1063/5.0197825] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/14/2024] [Accepted: 03/11/2024] [Indexed: 04/16/2024] Open
Abstract
The photoexcitation dynamics of molecular materials on the 10-100 nm length scale depend on complex interactions between electronic and vibrational degrees of freedom, rendering exact calculations difficult or intractable. The adaptive Hierarchy of Pure States (adHOPS) is a formally exact method that leverages the locality imposed by interactions between thermal environments and electronic excitations to achieve size-invariant scaling calculations for single-excitation processes in systems described by a Frenkel-Holstein Hamiltonian. Here, we extend adHOPS to account for arbitrary couplings between thermal environments and vertical excitation energies, enabling formally exact, size-invariant calculations that involve multiple excitations or states with shared thermal environments. In addition, we introduce a low-temperature correction and an effective integration of the noise to reduce the computational expense of including ultrafast vibrational relaxation in Hierarchy of Pure States (HOPS) simulations. We present these advances in the latest version of the open-source MesoHOPS library and use MesoHOPS to characterize charge separation at a one-dimensional organic heterojunction when both the electron and hole are mobile.
Collapse
Affiliation(s)
- Brian Citty
- Department of Chemistry, University of Texas at Austin, Austin, Texas 78712, USA
| | - Jacob K Lynd
- Department of Chemistry, University of Texas at Austin, Austin, Texas 78712, USA
| | - Tarun Gera
- Department of Chemistry, University of Texas at Austin, Austin, Texas 78712, USA
| | - Leonel Varvelo
- Department of Chemistry, Southern Methodist University, PO Box 750314 Dallas, Texas 75205, USA
| | - Doran I G B Raccah
- Department of Chemistry, University of Texas at Austin, Austin, Texas 78712, USA
| |
Collapse
|
2
|
Lorenzoni N, Cho N, Lim J, Tamascelli D, Huelga SF, Plenio MB. Systematic Coarse Graining of Environments for the Nonperturbative Simulation of Open Quantum Systems. PHYSICAL REVIEW LETTERS 2024; 132:100403. [PMID: 38518302 DOI: 10.1103/physrevlett.132.100403] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 10/23/2023] [Revised: 01/22/2024] [Accepted: 02/13/2024] [Indexed: 03/24/2024]
Abstract
Conducting precise electronic-vibrational dynamics simulations of molecular systems poses significant challenges when dealing with realistic environments composed of numerous vibrational modes. Here, we introduce a technique for the construction of effective phonon spectral densities that capture accurately open-system dynamics over a finite time interval of interest. When combined with existing nonperturbative simulation tools, our approach can reduce significantly the computational costs associated with many-body open-system dynamics.
Collapse
Affiliation(s)
- Nicola Lorenzoni
- Institut für Theoretische Physik und IQST, Albert-Einstein-Allee 11, Universität Ulm, D-89081 Ulm, Germany
| | - Namgee Cho
- Institut für Theoretische Physik und IQST, Albert-Einstein-Allee 11, Universität Ulm, D-89081 Ulm, Germany
| | - James Lim
- Institut für Theoretische Physik und IQST, Albert-Einstein-Allee 11, Universität Ulm, D-89081 Ulm, Germany
| | - Dario Tamascelli
- Institut für Theoretische Physik und IQST, Albert-Einstein-Allee 11, Universität Ulm, D-89081 Ulm, Germany
- Dipartimento di Fisica "Aldo Pontremoli," Università degli Studi di Milano, Via Celoria 16, 20133 Milano, Italy
| | - Susana F Huelga
- Institut für Theoretische Physik und IQST, Albert-Einstein-Allee 11, Universität Ulm, D-89081 Ulm, Germany
| | - Martin B Plenio
- Institut für Theoretische Physik und IQST, Albert-Einstein-Allee 11, Universität Ulm, D-89081 Ulm, Germany
| |
Collapse
|
3
|
Dey D, Woodhouse JL, Taylor MP, Fielding HH, Worth GA. On the multiphoton ionisation photoelectron spectra of phenol. Phys Chem Chem Phys 2024; 26:3451-3461. [PMID: 38205824 DOI: 10.1039/d3cp05559k] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/12/2024]
Abstract
The phenol molecule is a prototype for non-adiabatic dynamics and the excited-state photochemistry of biomolecules. In this article, we report a joint theoretical and experimental investigation on the resonance enhanced multiphoton ionisation photoelectron (REMPI) spectra of the two lowest ionisation bands of phenol. The focus is on the theoretical interpretation of the measured spectra using quantum dynamics simulations. These were performed by numerically solving the time-dependent Schrödinger equation using the multi-layer variant of the multiconfiguration time-dependent Hartree algorithm together with a vibronic coupling Hamiltonian model. The ionising laser pulse is modelled explicitly within the ionisation continuum model to simulate experimental femtosecond 1+1 REMPI photoelectron spectra. These measured spectra are sensitive to very short lived electronically excited states, providing a rigorous benchmark for our theoretical methods. The match between experiment and theory allows for an interpretation of the features of the spectra at different wavelengths and shows that there are features due to both 'direct' and 'indirect' ionisation, resulting from non-resonant and resonant excitation by the pump pulse.
Collapse
Affiliation(s)
- Diptesh Dey
- Department of Chemistry, University College London, 20 Gordon Street, London WC1H 0AJ, UK.
- Department of Chemistry, Northwestern University, Evanston, Illinois 60208, USA
| | - Joanne L Woodhouse
- Department of Chemistry, University College London, 20 Gordon Street, London WC1H 0AJ, UK.
- Department of Chemistry, University of Southampton, Southampton, SO17 1BJ, UK
| | - Marcus P Taylor
- Department of Chemistry, University College London, 20 Gordon Street, London WC1H 0AJ, UK.
| | - Helen H Fielding
- Department of Chemistry, University College London, 20 Gordon Street, London WC1H 0AJ, UK.
| | - Graham A Worth
- Department of Chemistry, University College London, 20 Gordon Street, London WC1H 0AJ, UK.
| |
Collapse
|
4
|
Hoppe H, Manthe U. Eigenstate calculation in the state-averaged (multi-layer) multi-configurational time-dependent Hartree approach. J Chem Phys 2024; 160:034104. [PMID: 38230812 DOI: 10.1063/5.0188748] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/23/2023] [Accepted: 12/15/2023] [Indexed: 01/18/2024] Open
Abstract
A new approach for the calculation of eigenstates with the state-averaged (multi-layer) multi-configurational time-dependent Hartree (MCTDH) approach is presented. The approach is inspired by the recent work of Larsson [J. Chem. Phys. 151, 204102 (2019)]. It employs local optimization of the basis sets at each node of the multi-layer MCTDH tree and successive downward and upward sweeps to obtain a globally converged result. At the top node, the Hamiltonian represented in the basis of the single-particle functions (SPFs) of the first layer is diagonalized. Here p wavefunctions corresponding to the p lowest eigenvalues are computed by a block Lanczos approach. At all other nodes, a non-linear operator consisting of the respective mean-field Hamiltonian matrix and a projector onto the space spanned by the respective SPFs is considered. Here, the eigenstate corresponding to the lowest eigenvalue is computed using a short iterative Lanczos scheme. Two different examples are studied to illustrate the new approach: the calculation of the vibrational states of methyl and acetonitrile. The calculations for methyl employ the single-layer MCTDH approach, a general potential energy surface, and the correlation discrete variable representation. A five-layer MCTDH representation and a sum of product-type Hamiltonian are used in the acetonitrile calculations. Very fast convergence and order of magnitude reductions in the numerical effort compared to the previously used block relaxation scheme are found. Furthermore, a detailed comparison with the results of Avila and Carrington [J. Chem. Phys. 134, 054126 (2011)] for acetonitrile highlights the potential problems of convergence tests for high-dimensional systems.
Collapse
Affiliation(s)
- Hannes Hoppe
- Theoretische Chemie, Fakultät für Chemie, Universität Bielefeld, Universitätsstr. 25, D-33615 Bielefeld, Germany
| | - Uwe Manthe
- Theoretische Chemie, Fakultät für Chemie, Universität Bielefeld, Universitätsstr. 25, D-33615 Bielefeld, Germany
| |
Collapse
|
5
|
Lyu N, Miano A, Tsioutsios I, Cortiñas RG, Jung K, Wang Y, Hu Z, Geva E, Kais S, Batista VS. Mapping Molecular Hamiltonians into Hamiltonians of Modular cQED Processors. J Chem Theory Comput 2023; 19:6564-6576. [PMID: 37733472 DOI: 10.1021/acs.jctc.3c00620] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 09/23/2023]
Abstract
We introduce a general method based on the operators of the Dyson-Masleev transformation to map the Hamiltonian of an arbitrary model system into the Hamiltonian of a circuit Quantum Electrodynamics (cQED) processor. Furthermore, we introduce a modular approach to programming a cQED processor with components corresponding to the mapping Hamiltonian. The method is illustrated as applied to quantum dynamics simulations of the Fenna-Matthews-Olson (FMO) complex and the spin-boson model of charge transfer. Beyond applications to molecular Hamiltonians, the mapping provides a general approach to implement any unitary operator in terms of a sequence of unitary transformations corresponding to powers of creation and annihilation operators of a single bosonic mode in a cQED processor.
Collapse
Affiliation(s)
- Ningyi Lyu
- Department of Chemistry, Yale University, New Haven, Connecticut 06520, United States
| | - Alessandro Miano
- Department of Applied Physics, Yale University, New Haven, Connecticut 06520, United States
- Department of Physics, Yale University, New Haven, Connecticut 06520, United States
- Yale Quantum Institute, Yale University, New Haven, Connecticut 06511, United States
| | - Ioannis Tsioutsios
- Department of Applied Physics, Yale University, New Haven, Connecticut 06520, United States
- Department of Physics, Yale University, New Haven, Connecticut 06520, United States
- Yale Quantum Institute, Yale University, New Haven, Connecticut 06511, United States
| | - Rodrigo G Cortiñas
- Department of Applied Physics, Yale University, New Haven, Connecticut 06520, United States
- Department of Physics, Yale University, New Haven, Connecticut 06520, United States
- Yale Quantum Institute, Yale University, New Haven, Connecticut 06511, United States
| | - Kenneth Jung
- Department of Chemistry, Yale University, New Haven, Connecticut 06520, 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
| | - Zixuan Hu
- 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
| | - Sabre Kais
- Department of Chemistry, Department of Physics and Purdue Quantum Science and Engineering Institute, Purdue University, West Lafayette, Indiana 47907, 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
|
6
|
Kim Y, Mitchell Z, Lawrence J, Morozov D, Savikhin S, Slipchenko LV. Predicting Mutation-Induced Changes in the Electronic Properties of Photosynthetic Proteins from First Principles: The Fenna-Matthews-Olson Complex Example. J Phys Chem Lett 2023; 14:7038-7044. [PMID: 37524046 DOI: 10.1021/acs.jpclett.3c01461] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 08/02/2023]
Abstract
Multiscale molecular modeling is utilized to predict optical absorption and circular dichroism spectra of two single-point mutants of the Fenna-Matthews-Olson photosynthetic pigment-protein complex. The modeling approach combines classical molecular dynamics simulations with structural refinement of photosynthetic pigments and calculations of their excited states in a polarizable protein environment. The only experimental input to the modeling protocol is the X-ray structure of the wild-type protein. The first-principles modeling reproduces changes in the experimental optical spectra of the considered mutants, Y16F and Q198V. Interestingly, the Q198V mutation has a negligible effect on the electronic properties of the targeted bacteriochlorophyll a pigment. Instead, the electronic properties of several other pigments respond to this mutation. The molecular modeling demonstrates that a single-point mutation can induce long-range effects on the protein structure, while extensive structural changes near a pigment do not necessarily lead to significant changes in the electronic properties of that pigment.
Collapse
Affiliation(s)
- Yongbin Kim
- Department of Chemistry, Purdue University, 560 Oval Drive, West Lafayette, Indiana 47907, United States
| | - Zach Mitchell
- Department of Physics and Astronomy, Purdue University, 525 Northwestern Avenue, West Lafayette, Indiana 47907, United States
| | - Jack Lawrence
- Department of Chemistry, Purdue University, 560 Oval Drive, West Lafayette, Indiana 47907, United States
| | - Dmitry Morozov
- Nanoscience Center and Department of Chemistry, University of Jyväskylä, P.O. Box 35, 40014 Jyväskylä, Finland
| | - Sergei Savikhin
- Department of Physics and Astronomy, Purdue University, 525 Northwestern Avenue, West Lafayette, Indiana 47907, United States
| | - Lyudmila V Slipchenko
- Department of Chemistry, Purdue University, 560 Oval Drive, West Lafayette, Indiana 47907, United States
| |
Collapse
|
7
|
Kaiser A, Daoud RE, Aquilante F, Kühn O, De Vico L, Bokarev SI. A Multiconfigurational Wave Function Implementation of the Frenkel Exciton Model for Molecular Aggregates. J Chem Theory Comput 2023; 19:2918-2928. [PMID: 37115036 DOI: 10.1021/acs.jctc.3c00185] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 04/29/2023]
Abstract
We present an implementation of the Frenkel exciton model into the OpenMolcas program package enabling calculations of collective electronic excited states of molecular aggregates based on a multiconfigurational wave function description of the individual monomers. The computational protocol avoids using diabatization schemes and, thus, supermolecule calculations. Additionally, the use of the Cholesky decomposition of the two-electron integrals entering pair interactions enhances the efficiency of the computational scheme. The application of the method is exemplified for two test systems, that is, a formaldehyde oxime and a bacteriochlorophyll-like dimer. For the sake of comparison with the dipole approximation, we restrict our considerations to situations where intermonomer exchange can be neglected. The protocol is expected to be beneficial for aggregates composed of molecules with extended π systems, unpaired electrons such as radicals or transition metal centers, where it should outperform widely used methods based on time-dependent density functional theory.
Collapse
Affiliation(s)
- Andy Kaiser
- Institut für Physik, Universität Rostock, Albert-Einstein-Str. 23-24, 18059 Rostock. Germany
| | - Razan E Daoud
- Dipartimento di Biotecnologie, Chimica e Farmacia, Università degli Studi di Siena, Via A. Moro 2, 53100 Siena, Italy
| | - Francesco Aquilante
- Theory and Simulation of Materials (THEOS) and National Centre for Computational Design and Discovery of Novel Materials (MARVEL), École Polytechnique Fédérale de Lausanne (EPFL), CH-1015 Lausanne, Switzerland
| | - Oliver Kühn
- Institut für Physik, Universität Rostock, Albert-Einstein-Str. 23-24, 18059 Rostock. Germany
| | - Luca De Vico
- Dipartimento di Biotecnologie, Chimica e Farmacia, Università degli Studi di Siena, Via A. Moro 2, 53100 Siena, Italy
| | - Sergey I Bokarev
- Institut für Physik, Universität Rostock, Albert-Einstein-Str. 23-24, 18059 Rostock. Germany
- Chemistry Department, Technical University of Munich, Lichtenbergstr. 4, 85748 Garching, Germany
| |
Collapse
|
8
|
Malpathak S, Ananth N. Non-linear correlation functions and zero-point energy flow in mixed quantum-classical semiclassical dynamics. J Chem Phys 2023; 158:104106. [PMID: 36922136 DOI: 10.1063/5.0133222] [Citation(s) in RCA: 3] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/17/2023] Open
Abstract
Mixed quantum classical (MQC)-initial value representation (IVR) is a recently introduced semiclassical framework that allows for selective quantization of the modes of a complex system. In the quantum limit, MQC reproduces the semiclassical Double Herman-Kluk IVR results, accurately capturing nuclear quantum coherences and conserving zero-point energy. However, in the classical limit, although MQC mimics the Husimi-IVR for real-time correlation functions with linear operators, it is significantly less accurate for non-linear correlation functions with errors even at time zero. Here, we identify the origin of this discrepancy in the MQC formulation and propose a modification. We analytically show that the modified MQC approach is exact for all correlation functions at time zero, and in a study of zero-point energy (ZPE) flow, we numerically demonstrate that it correctly obtains the quantum and classical limits as a function of time. Interestingly, although classical-limit MQC simulations show the expected, unphysical ZPE leakage, we find that it is possible to predict and even modify the direction of ZPE flow through selective quantization of the system, with the quantum-limit modes accepting energy but preserving the minimum quantum mechanically required energy.
Collapse
Affiliation(s)
- Shreyas Malpathak
- Department of Chemistry and Chemical Biology, Baker Laboratory, Cornell University, Ithaca, New York 14853, USA
| | - Nandini Ananth
- Department of Chemistry and Chemical Biology, Baker Laboratory, Cornell University, Ithaca, New York 14853, USA
| |
Collapse
|
9
|
Gelin MF, Chen L, Domcke W. Equation-of-Motion Methods for the Calculation of Femtosecond Time-Resolved 4-Wave-Mixing and N-Wave-Mixing Signals. Chem Rev 2022; 122:17339-17396. [PMID: 36278801 DOI: 10.1021/acs.chemrev.2c00329] [Citation(s) in RCA: 17] [Impact Index Per Article: 8.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/25/2023]
Abstract
Femtosecond nonlinear spectroscopy is the main tool for the time-resolved detection of photophysical and photochemical processes. Since most systems of chemical interest are rather complex, theoretical support is indispensable for the extraction of the intrinsic system dynamics from the detected spectroscopic responses. There exist two alternative theoretical formalisms for the calculation of spectroscopic signals, the nonlinear response-function (NRF) approach and the spectroscopic equation-of-motion (EOM) approach. In the NRF formalism, the system-field interaction is assumed to be sufficiently weak and is treated in lowest-order perturbation theory for each laser pulse interacting with the sample. The conceptual alternative to the NRF method is the extraction of the spectroscopic signals from the solutions of quantum mechanical, semiclassical, or quasiclassical EOMs which govern the time evolution of the material system interacting with the radiation field of the laser pulses. The NRF formalism and its applications to a broad range of material systems and spectroscopic signals have been comprehensively reviewed in the literature. This article provides a detailed review of the suite of EOM methods, including applications to 4-wave-mixing and N-wave-mixing signals detected with weak or strong fields. Under certain circumstances, the spectroscopic EOM methods may be more efficient than the NRF method for the computation of various nonlinear spectroscopic signals.
Collapse
Affiliation(s)
- Maxim F Gelin
- School of Science, Hangzhou Dianzi University, Hangzhou 310018, China
| | - Lipeng Chen
- Max-Planck-Institut für Physik komplexer Systeme, Nöthnitzer Strasse 38, D-01187 Dresden, Germany
| | - Wolfgang Domcke
- Department of Chemistry, Technical University of Munich, D-85747 Garching,Germany
| |
Collapse
|
10
|
Hino K, Kurashige Y. Matrix Product State Formulation of the MCTDH Theory in Local Mode Representations for Anharmonic Potentials. J Chem Theory Comput 2022; 18:3347-3356. [PMID: 35606892 DOI: 10.1021/acs.jctc.2c00243] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
The matrix product state formulation of the multiconfiguration time-dependent Hartree theory, MPS-MCTDH, reported previously [Kurashige, J. Chem. Phys. 2018, 19, 194114] is extended to realistic anharmonic potentials with n-mode representations beyond the linear vibronic coupling model. For realistic vibrational potentials, the local mode representation should give a more compact representation of the potentials, i.e., lowering the dimensionality of the entanglements, than the normal coordinates, and the MPS-MCTDH formulation should work more efficiently and maintain the accuracy with a small bond dimension of the MPS ansatz. In fact, it was confirmed that the use of the local coordinates made the interaction matrices diagonal dominant and the number of terms in the n-body expansion of the potentials was significantly reduced. The method was applied to the IR spectrum of the CH2O molecule, the zero-point energies, and the vibrational energy redistribution dynamics of polyenes C2nH2n+2. The results showed that the efficiency of the MPS-MCTDH method is significantly accelerated by the use of local coordinates even if the long-range interactions are included in the potential.
Collapse
Affiliation(s)
- Kentaro Hino
- Department of Chemistry, Graduate School of Science, Kyoto University, Kitashirakawa Oiwake-cho, Sakyo-ku, Kyoto 606-8502, Japan
| | - Yuki Kurashige
- Department of Chemistry, Graduate School of Science, Kyoto University, Kitashirakawa Oiwake-cho, Sakyo-ku, Kyoto 606-8502, Japan
| |
Collapse
|
11
|
Predicting the future of excitation energy transfer in light-harvesting complex with artificial intelligence-based quantum dynamics. Nat Commun 2022; 13:1930. [PMID: 35411054 PMCID: PMC9001686 DOI: 10.1038/s41467-022-29621-w] [Citation(s) in RCA: 13] [Impact Index Per Article: 6.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/22/2021] [Accepted: 03/16/2022] [Indexed: 01/20/2023] Open
Abstract
Exploring excitation energy transfer (EET) in light-harvesting complexes (LHCs) is essential for understanding the natural processes and design of highly-efficient photovoltaic devices. LHCs are open systems, where quantum effects may play a crucial role for almost perfect utilization of solar energy. Simulation of energy transfer with inclusion of quantum effects can be done within the framework of dissipative quantum dynamics (QD), which are computationally expensive. Thus, artificial intelligence (AI) offers itself as a tool for reducing the computational cost. Here we suggest AI-QD approach using AI to directly predict QD as a function of time and other parameters such as temperature, reorganization energy, etc., completely circumventing the need of recursive step-wise dynamics propagation in contrast to the traditional QD and alternative, recursive AI-based QD approaches. Our trajectory-learning AI-QD approach is able to predict the correct asymptotic behavior of QD at infinite time. We demonstrate AI-QD on seven-sites Fenna–Matthews–Olson (FMO) complex. Simulations of energy transfer in light-harvesting complexes are computationally very demanding. Here the authors apply an artificial intelligence quantum dissipative algorithm to study the excited state energy transfer dynamics in a light-harvesting complex.
Collapse
|
12
|
Ellerbrock R, Manthe U. A non-hierarchical correlation discrete variable representation. J Chem Phys 2022; 156:134107. [PMID: 35395891 DOI: 10.1063/5.0088509] [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/15/2022] Open
Abstract
The correlation discrete variable representation (CDVR) facilitates (multi-layer) multi-configurational time-dependent Hartree (MCTDH) calculations with general potentials. It employs a layered grid representation to efficiently evaluate all potential matrix elements appearing in the MCTDH equations of motion. The original CDVR approach and its multi-layer extension show a hierarchical structure: the size of the grids employed at the different layers increases when moving from an upper layer to a lower one. In this work, a non-hierarchical CDVR approach, which uses identically structured quadratures at all layers of the MCTDH wavefunction representation, is introduced. The non-hierarchical CDVR approach crucially reduces the number of grid points required, compared to the hierarchical CDVR, shows superior scaling properties, and yields identical results for all three representations showing the same topology. Numerical tests studying the photodissociation of NOCl and the vibrational states of CH3 demonstrate the accuracy of the non-hierarchical CDVR approach.
Collapse
Affiliation(s)
- Roman Ellerbrock
- Department of Chemistry and The PULSE Institute, Stanford University, Stanford, California 94305, USA
| | - Uwe Manthe
- Theoretische Chemie, Fakultät für Chemie, Universität Bielefeld, Universitätsstr. 25, D-33615 Bielefeld, Germany
| |
Collapse
|
13
|
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
|
14
|
Serwatka T, Roy PN. Ground state of asymmetric tops with DMRG: Water in one dimension. J Chem Phys 2022; 156:044116. [DOI: 10.1063/5.0078770] [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] Open
Affiliation(s)
- Tobias Serwatka
- Department of Chemistry, University of Waterloo, Waterloo, Ontario N2L 3G1, Canada
| | - Pierre-Nicholas Roy
- Department of Chemistry, University of Waterloo, Waterloo, Ontario N2L 3G1, Canada
| |
Collapse
|
15
|
Soley MB, Bergold P, Gorodetsky AA, Batista VS. Functional Tensor-Train Chebyshev Method for Multidimensional Quantum Dynamics Simulations. J Chem Theory Comput 2021; 18:25-36. [PMID: 34898201 DOI: 10.1021/acs.jctc.1c00941] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Methods for efficient simulations of multidimensional quantum dynamics are essential for theoretical studies of chemical systems where quantum effects are important, such as those involving rearrangements of protons or electronic configurations. Here, we introduce the functional tensor-train Chebyshev (FTTC) method for rigorous nuclear quantum dynamics simulations. FTTC is essentially the Chebyshev propagation scheme applied to the initial state represented in a continuous analogue tensor-train format. We demonstrate the capabilities of FTTC as applied to simulations of proton quantum dynamics in a 50-dimensional model of hydrogen-bonded DNA base pairs.
Collapse
Affiliation(s)
- Micheline B Soley
- Yale Quantum Institute, Yale University, P.O. Box 208334, New Haven, Connecticut 06520-8263, United States.,Department of Chemistry, Yale University, P.O. Box 208107, New Haven, Connecticut 06520, United States
| | - Paul Bergold
- Zentrum Mathematik, Technical University of Munich, Boltzmannstr. 3, 85748 Garching, Germany
| | - Alex A Gorodetsky
- Department of Aerospace Engineering, University of Michigan, 1320 Beal Avenue, Ann Arbor, Michigan 48109-2140, United States
| | - Victor S Batista
- Yale Quantum Institute, Yale University, P.O. Box 208334, New Haven, Connecticut 06520-8263, United States.,Department of Chemistry, Yale University, P.O. Box 208107, New Haven, Connecticut 06520, United States.,Energy Sciences Institute, Yale University, P.O. Box 27394, West Haven, Connecticut 06516-7394, United States
| |
Collapse
|
16
|
Cho KH, Rhee YM. Computational elucidations on the role of vibrations in energy transfer processes of photosynthetic complexes. Phys Chem Chem Phys 2021; 23:26623-26639. [PMID: 34842245 DOI: 10.1039/d1cp04615b] [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/21/2022]
Abstract
Coupling between pigment excitations and nuclear movements in photosynthetic complexes is known to modulate the excitation energy transfer (EET) efficiencies. Toward providing microscopic information, researchers often apply simulation techniques and investigate how vibrations are involved in EET processes. Here, reports on such roles of nuclear movements are discussed from a theory perspective. While vibrations naturally present random thermal fluctuations that can affect energy transferring characteristics, they can also be intertwined with exciton structures and create more specific non-adiabatic energy transfer pathways. For reliable simulations, a bath model that accurately mimics a given molecular system is required. Methods for obtaining such a model in combination with quantum chemical electronic structure calculations and molecular dynamics trajectory simulations are discussed. Various quantum dynamics simulation tools that can handle pigment-to-pigment energy transfers together with their vibrational characters are also touched on. Behaviors of molecular vibrations often deviate from ideality, especially when all-atom details are included, which practically forces us to treat them classically. We conclude this perspective by considering some recent reports that suggest that classical descriptions of bath effects with all-atom details may still produce valuable information for analyzing sophisticated contributions by vibrations to EET processes.
Collapse
Affiliation(s)
- Kwang Hyun Cho
- Department of Chemistry, Korea Advanced Institute of Science and Technology (KAIST), Daejeon 34141, Korea.
| | - Young Min Rhee
- Department of Chemistry, Korea Advanced Institute of Science and Technology (KAIST), Daejeon 34141, Korea.
| |
Collapse
|
17
|
Lindoy LP, Kloss B, Reichman DR. Time evolution of ML-MCTDH wavefunctions. II. Application of the projector splitting integrator. J Chem Phys 2021; 155:174109. [PMID: 34742222 DOI: 10.1063/5.0070043] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
The multi-layer multi-configuration time-dependent Hartree (ML-MCTDH) approach can suffer from numerical instabilities whenever the wavefunction is weakly entangled. These instabilities arise from singularities in the equations of motion (EOMs) and necessitate the use of regularization of the EOMs. The Projector Splitting Integrator (PSI) has previously been presented as an approach for evolving ML-MCTDH wavefunctions that is free of singularities. Here, we will discuss the implementation of the multi-layer PSI with a particular focus on how the steps required relate to those required to implement standard ML-MCTDH. We demonstrate the efficiency and stability of the PSI for large ML-MCTDH wavefunctions containing up to hundreds of thousands of nodes by considering a series of spin-boson models with up to 106 bath modes and find that for these problems, the PSI requires roughly 3-4 orders of magnitude fewer Hamiltonian evaluations and 2-3 orders of magnitude fewer Hamiltonian applications than standard ML-MCTDH and 2-3/1-2 orders of magnitude fewer evaluations/applications than approaches that use improved regularization schemes. Finally, we consider a series of significantly more challenging multi-spin-boson models that require much larger numbers of single-particle functions with wavefunctions containing up to ∼1.3×109 parameters to obtain accurate dynamics.
Collapse
Affiliation(s)
- Lachlan P Lindoy
- Department of Chemistry, Columbia University, 3000 Broadway, New York, New York 10027, USA
| | - Benedikt Kloss
- Department of Chemistry, Columbia University, 3000 Broadway, New York, New York 10027, USA
| | - David R Reichman
- Department of Chemistry, Columbia University, 3000 Broadway, New York, New York 10027, USA
| |
Collapse
|
18
|
Lindoy LP, Kloss B, Reichman DR. Time evolution of ML-MCTDH wavefunctions. I. Gauge conditions, basis functions, and singularities. J Chem Phys 2021; 155:174108. [PMID: 34742180 DOI: 10.1063/5.0070042] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/27/2023] Open
Abstract
We derive a family of equations-of-motion (EOMs) for evolving multi-layer multiconfiguration time-dependent Hartree (ML-MCTDH) wavefunctions that, unlike the standard ML-MCTDH EOMs, never require the evaluation of the inverse of singular matrices. All members of this family of EOMs make use of alternative static gauge conditions than those used for standard ML-MCTDH. These alternative conditions result in an expansion of the wavefunction in terms of a set of potentially arbitrary orthonormal functions, rather than in terms of a set of non-orthonormal and potentially linearly dependent functions, as is the case for standard ML-MCTDH. We show that the EOMs used in the projector splitting integrator (PSI) and the invariant EOM approaches are two special cases of this family obtained from different choices for the dynamic gauge condition, with the invariant EOMs making use of a choice that introduces potentially unbounded operators into the EOMs. As a consequence, all arguments for the existence of parallelizable integration schemes for the invariant EOMs can also be applied to the PSI EOMs.
Collapse
Affiliation(s)
- Lachlan P Lindoy
- Department of Chemistry, Columbia University, 3000 Broadway, New York, New York 10027, USA
| | - Benedikt Kloss
- Department of Chemistry, Columbia University, 3000 Broadway, New York, New York 10027, USA
| | - David R Reichman
- Department of Chemistry, Columbia University, 3000 Broadway, New York, New York 10027, USA
| |
Collapse
|
19
|
Weike T, Manthe U. The multi-configurational time-dependent Hartree approach in optimized second quantization: thermal ensembles and statistical sampling. Chem Phys 2021. [DOI: 10.1016/j.chemphys.2021.111413] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/27/2022]
|
20
|
Wang C, Ulusoy IS, Aebersold LE, Wilson AK. Multi-configuration electron-nuclear dynamics: An open-shell approach. J Chem Phys 2021; 155:154103. [PMID: 34686063 DOI: 10.1063/5.0063478] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
The multi-configuration electron-nuclear dynamics for open shell systems with a spin-unrestricted formalism is described. The mean fields are evaluated using second-order reduced density matrices for electronic and nuclear degrees of freedom. Applications to light-element diatomics including equilibrium geometries, electronic energies, dipole moments, and absorption spectra are presented. The von Neumann entropies for different spin states of a LiH molecule are compared.
Collapse
Affiliation(s)
- Cong Wang
- Department of Chemistry, Michigan State University, 578 S Shaw Lane, East Lansing, Michigan 48824-1322, USA
| | - Inga S Ulusoy
- Department of Chemistry, Michigan State University, 578 S Shaw Lane, East Lansing, Michigan 48824-1322, USA
| | - Lucas E Aebersold
- Department of Chemistry, Michigan State University, 578 S Shaw Lane, East Lansing, Michigan 48824-1322, USA
| | - Angela K Wilson
- Department of Chemistry, Michigan State University, 578 S Shaw Lane, East Lansing, Michigan 48824-1322, USA
| |
Collapse
|
21
|
MacDonell RJ, Dickerson CE, Birch CJT, Kumar A, Edmunds CL, Biercuk MJ, Hempel C, Kassal I. Analog quantum simulation of chemical dynamics. Chem Sci 2021; 12:9794-9805. [PMID: 34349953 PMCID: PMC8293981 DOI: 10.1039/d1sc02142g] [Citation(s) in RCA: 10] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/16/2021] [Accepted: 06/14/2021] [Indexed: 11/21/2022] Open
Abstract
Ultrafast chemical reactions are difficult to simulate because they involve entangled, many-body wavefunctions whose computational complexity grows rapidly with molecular size. In photochemistry, the breakdown of the Born-Oppenheimer approximation further complicates the problem by entangling nuclear and electronic degrees of freedom. Here, we show that analog quantum simulators can efficiently simulate molecular dynamics using commonly available bosonic modes to represent molecular vibrations. Our approach can be implemented in any device with a qudit controllably coupled to bosonic oscillators and with quantum hardware resources that scale linearly with molecular size, and offers significant resource savings compared to digital quantum simulation algorithms. Advantages of our approach include a time resolution orders of magnitude better than ultrafast spectroscopy, the ability to simulate large molecules with limited hardware using a Suzuki-Trotter expansion, and the ability to implement realistic system-bath interactions with only one additional interaction per mode. Our approach can be implemented with current technology; e.g., the conical intersection in pyrazine can be simulated using a single trapped ion. Therefore, we expect our method will enable classically intractable chemical dynamics simulations in the near term.
Collapse
Affiliation(s)
- Ryan J MacDonell
- School of Chemistry, University of Sydney NSW 2006 Australia
- University of Sydney Nano Institute, University of Sydney NSW 2006 Australia
| | - Claire E Dickerson
- School of Chemistry, University of Sydney NSW 2006 Australia
- School of Physics, University of Sydney NSW 2006 Australia
- ARC Centre of Excellence for Engineered Quantum Systems, University of Sydney NSW 2006 Australia
- University of Sydney Nano Institute, University of Sydney NSW 2006 Australia
| | - Clare J T Birch
- School of Chemistry, University of Sydney NSW 2006 Australia
- University of Sydney Nano Institute, University of Sydney NSW 2006 Australia
| | - Alok Kumar
- Department of Chemistry, Indian Institute of Technology Bombay Mumbai 400076 India
| | - Claire L Edmunds
- School of Physics, University of Sydney NSW 2006 Australia
- ARC Centre of Excellence for Engineered Quantum Systems, University of Sydney NSW 2006 Australia
| | - Michael J Biercuk
- School of Physics, University of Sydney NSW 2006 Australia
- ARC Centre of Excellence for Engineered Quantum Systems, University of Sydney NSW 2006 Australia
| | - Cornelius Hempel
- School of Physics, University of Sydney NSW 2006 Australia
- ARC Centre of Excellence for Engineered Quantum Systems, University of Sydney NSW 2006 Australia
- University of Sydney Nano Institute, University of Sydney NSW 2006 Australia
| | - Ivan Kassal
- School of Chemistry, University of Sydney NSW 2006 Australia
- University of Sydney Nano Institute, University of Sydney NSW 2006 Australia
| |
Collapse
|
22
|
Varvelo L, Lynd JK, Bennett DIG. Formally exact simulations of mesoscale exciton dynamics in molecular materials. Chem Sci 2021; 12:9704-9711. [PMID: 34349941 PMCID: PMC8293828 DOI: 10.1039/d1sc01448j] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/11/2021] [Accepted: 05/31/2021] [Indexed: 02/04/2023] Open
Abstract
Excited state carriers, such as excitons, can diffuse on the 100 nm to micron length scale in molecular materials but only delocalize over short length scales due to coupling between electronic and vibrational degrees-of-freedom. Here, we leverage the locality of excitons to adaptively solve the hierarchy of pure states equations (HOPS). We demonstrate that our adaptive HOPS (adHOPS) methodology provides a formally exact and size-invariant (i.e., ) scaling algorithm for simulating mesoscale quantum dynamics. Finally, we provide proof-of-principle calculations for exciton diffusion on linear chains containing up to 1000 molecules. The adaptive hierarchy of pure states (adHOPS) algorithm leverages the locality of excitons in molecular materials to perform formally-exact simulations with size-invariant (i.e., ) scaling, enabling efficient simulations of mesoscale exciton dynamics.![]()
Collapse
Affiliation(s)
- Leonel Varvelo
- Department of Chemistry, Southern Methodist University PO Box 750314 Dallas TX USA
| | - Jacob K Lynd
- Department of Chemistry, Southern Methodist University PO Box 750314 Dallas TX USA
| | - Doran I G Bennett
- Department of Chemistry, Southern Methodist University PO Box 750314 Dallas TX USA
| |
Collapse
|
23
|
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
|
24
|
Weike T, Manthe U. Symmetries in the multi-configurational time-dependent Hartree wavefunction representation and propagation. J Chem Phys 2021; 154:194108. [PMID: 34240912 DOI: 10.1063/5.0054105] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
In multi-configurational time-dependent Hartree (MCTDH) approaches, different multi-layered wavefunction representations can be used to represent the same physical wavefunction. Transformations between different equivalent representations of a physical wavefunction that alter the tree structure used in the multi-layer MCTDH wavefunction representation interchange the role of single-particle functions (SPFs) and single-hole functions (SHFs) in the MCTDH formalism. While the physical wavefunction is invariant under these transformations, this invariance does not hold for the standard multi-layer MCTDH equations of motion. Introducing transformed SPFs, which obey normalization conditions typically associated with SHFs, revised equations of motion are derived. These equations do not show the singularities resulting from the inverse single-particle density matrix and are invariant under tree transformations. Based on the revised equations of motion, a new integration scheme is introduced. The scheme combines the advantages of the constant mean-field approach of Beck and Meyer [Z. Phys. D 42, 113 (1997)] and the singularity-free integrator suggested by Lubich [Appl. Math. Res. Express 2015, 311]. Numerical calculations studying the spin boson model in high dimensionality confirm the favorable properties of the new integration scheme.
Collapse
Affiliation(s)
- Thomas Weike
- Theoretische Chemie, Fakultät für Chemie, Universität Bielefeld, Universitätsstr. 25, D-33615 Bielefeld, Germany
| | - Uwe Manthe
- Theoretische Chemie, Fakultät für Chemie, Universität Bielefeld, Universitätsstr. 25, D-33615 Bielefeld, Germany
| |
Collapse
|
25
|
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
|
26
|
Mainali S, Gatti F, Iouchtchenko D, Roy PN, Meyer HD. Comparison of the multi-layer multi-configuration time-dependent Hartree (ML-MCTDH) method and the density matrix renormalization group (DMRG) for ground state properties of linear rotor chains. J Chem Phys 2021; 154:174106. [PMID: 34241072 DOI: 10.1063/5.0047090] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Abstract
We demonstrate the applicability of the Multi-Layer Multi-Configuration Time-Dependent Hartree (ML-MCTDH) method to the problem of computing ground states of one-dimensional chains of linear rotors with dipolar interactions. Specifically, we successfully obtain energies, entanglement entropies, and orientational correlations that are in agreement with the Density Matrix Renormalization Group (DMRG), which has been previously used for this system. We find that the entropies calculated by ML-MCTDH for larger system sizes contain nonmonotonicity, as expected in the vicinity of a second-order quantum phase transition between ordered and disordered rotor states. We observe that this effect remains when all couplings besides nearest-neighbor are omitted from the Hamiltonian, which suggests that it is not sensitive to the rate of decay of the interactions. In contrast to DMRG, which is tailored to the one-dimensional case, ML-MCTDH (as implemented in the Heidelberg MCTDH package) requires more computational time and memory, although the requirements are still within reach of commodity hardware. The numerical convergence and computational demand of two practical implementations of ML-MCTDH and DMRG are presented in detail for various combinations of system parameters.
Collapse
Affiliation(s)
- Samrit Mainali
- Université Paris-Saclay, Institut des Sciences Moléculaires d'Orsay ISMO, UMR CNRS 8214, F-91405 Orsay, France
| | - Fabien Gatti
- Université Paris-Saclay, Institut des Sciences Moléculaires d'Orsay ISMO, UMR CNRS 8214, F-91405 Orsay, France
| | - Dmitri Iouchtchenko
- Department of Chemistry, University of Waterloo, Waterloo, Ontario N2L 3G1, Canada
| | - Pierre-Nicholas Roy
- Department of Chemistry, University of Waterloo, Waterloo, Ontario N2L 3G1, Canada
| | - Hans-Dieter Meyer
- Theoretische Chemie, Physikalisch-Chemisches Institut, Universität Heidelberg, Im Neuenheimer Feld 229, 69120 Heidelberg, Germany
| |
Collapse
|
27
|
Popp W, Brey D, Binder R, Burghardt I. Quantum Dynamics of Exciton Transport and Dissociation in Multichromophoric Systems. Annu Rev Phys Chem 2021; 72:591-616. [PMID: 33636997 DOI: 10.1146/annurev-physchem-090419-040306] [Citation(s) in RCA: 22] [Impact Index Per Article: 7.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/09/2022]
Abstract
Due to the subtle interplay of site-to-site electronic couplings, exciton delocalization, nonadiabatic effects, and vibronic couplings, quantum dynamical studies are needed to elucidate the details of ultrafast photoinduced energy and charge transfer events in organic multichromophoric systems. In this vein, we review an approach that combines first-principles parameterized lattice Hamiltonians with accurate quantum dynamical simulations using advanced multiconfigurational methods. Focusing on the elementary transfer steps in organic functional materials, we address coherent exciton migration and creation of charge transfer excitons in homopolymers, notably representative of the poly(3-hexylthiophene) material, as well as exciton dissociation at polymer:fullerene heterojunctions. We emphasize the role of coherent transfer, trapping effects due to high-frequency phonon modes, and thermal activation due to low-frequency soft modes that drive a diffusive dynamics.
Collapse
Affiliation(s)
- Wjatscheslaw Popp
- Institute for Physical and Theoretical Chemistry, Goethe University Frankfurt, 60438 Frankfurt, Germany;
| | - Dominik Brey
- Institute for Physical and Theoretical Chemistry, Goethe University Frankfurt, 60438 Frankfurt, Germany;
| | - Robert Binder
- Institute for Physical and Theoretical Chemistry, Goethe University Frankfurt, 60438 Frankfurt, Germany;
| | - Irene Burghardt
- Institute for Physical and Theoretical Chemistry, Goethe University Frankfurt, 60438 Frankfurt, Germany;
| |
Collapse
|
28
|
Wang H, Meyer HD. Importance of Appropriately Regularizing the ML-MCTDH Equations of Motion. J Phys Chem A 2021; 125:3077-3087. [DOI: 10.1021/acs.jpca.0c11221] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Haobin Wang
- Department of Chemistry, University of Colorado Denver, Denver, Colorado 80217-3364, United States
| | - Hans-Dieter Meyer
- Physikalisch-Chemisches Institut, Universität Heidelberg, Im Neuenheimer Feld 229, D-69120 Heidelberg, Germany
| |
Collapse
|
29
|
Goswami S, Mundakkathparambil Parameswaran A, Köppel H. A quantum dynamical investigation of the excitation transfer in two doubly hydrogen-bonded molecular dimers. Mol Phys 2020. [DOI: 10.1080/00268976.2020.1762011] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/24/2022]
Affiliation(s)
- Sugata Goswami
- Theoretische Chemie, Physikalisch-Chemisches Institut, Universität Heidelberg, Heidelberg, Germany
| | | | - H. Köppel
- Theoretische Chemie, Physikalisch-Chemisches Institut, Universität Heidelberg, Heidelberg, Germany
| |
Collapse
|
30
|
Duan HG, Nalbach P, Miller RJD, Thorwart M. Intramolecular vibrations enhance the quantum efficiency of excitonic energy transfer. PHOTOSYNTHESIS RESEARCH 2020; 144:137-145. [PMID: 32306173 PMCID: PMC7203599 DOI: 10.1007/s11120-020-00742-x] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 12/26/2019] [Accepted: 03/26/2020] [Indexed: 06/11/2023]
Abstract
We study the impact of underdamped intramolecular vibrational modes on the efficiency of the excitation energy transfer in a dimer in which each state is coupled to its own underdamped vibrational mode and, in addition, to a continuous background of environmental modes. For this, we use the numerically exact hierarchy equation of motion approach. We determine the quantum yield and the transfer time in dependence of the vibronic coupling strength, and in dependence of the damping of the incoherent background. Moreover, we tune the vibrational frequencies out of resonance with the excitonic energy gap. We show that the quantum yield is enhanced by up to 10% when the vibrational frequency of the donor is larger than at the acceptor. The vibronic energy eigenstates of the acceptor acquire then an increased density of states, which leads to a higher occupation probability of the acceptor in thermal equilibrium. We can conclude that an underdamped vibrational mode which is weakly coupled to the dimer fuels a faster transfer of excitation energy, illustrating that long-lived vibrations can, in principle, enhance energy transfer, without involving long-lived electronic coherence.
Collapse
Affiliation(s)
- Hong-Guang Duan
- I. Institut für Theoretische Physik, Universität Hamburg, Jungiusstraße 9, 20355, Hamburg, Germany
- Max Planck Institute for the Structure and Dynamics of Matter, Luruper Chaussee 149, 22761, Hamburg, Germany
- The Hamburg Center for Ultrafast Imaging, Luruper Chaussee 149, 22761, Hamburg, Germany
| | - Peter Nalbach
- Westfälische Hochschule, Münsterstr. 265, 46397, Bocholt, Germany
| | - R J Dwayne Miller
- Max Planck Institute for the Structure and Dynamics of Matter, Luruper Chaussee 149, 22761, Hamburg, Germany
- The Hamburg Center for Ultrafast Imaging, Luruper Chaussee 149, 22761, Hamburg, Germany
- The Departments of Chemistry and Physics, University of Toronto, 80 St. George Street, Toronto, M5S 3H6, Canada
| | - Michael Thorwart
- I. Institut für Theoretische Physik, Universität Hamburg, Jungiusstraße 9, 20355, Hamburg, Germany.
- The Hamburg Center for Ultrafast Imaging, Luruper Chaussee 149, 22761, Hamburg, Germany.
| |
Collapse
|
31
|
|
32
|
Kim Y, Morozov D, Stadnytskyi V, Savikhin S, Slipchenko LV. Predictive First-Principles Modeling of a Photosynthetic Antenna Protein: The Fenna-Matthews-Olson Complex. J Phys Chem Lett 2020; 11:1636-1643. [PMID: 32013435 DOI: 10.1021/acs.jpclett.9b03486] [Citation(s) in RCA: 12] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/10/2023]
Abstract
High efficiency of light harvesting in photosynthetic pigment-protein complexes is governed by evolutionary-perfected protein-assisted tuning of individual pigment properties and interpigment interactions. Due to the large number of spectrally overlapping pigments in a typical photosynthetic complex, experimental methods often fail to unambiguously identify individual chromophore properties. Here, we report a first-principles-based modeling protocol capable of predicting properties of pigments in protein environment to a high precision. The technique was applied to successfully uncover electronic properties of the Fenna-Matthews-Olson (FMO) pigment-protein complex. Each of the three subunits of the FMO complex contains eight strongly coupled bacteriochlorophyll a (BChl a) pigments. The excitonic structure of FMO can be described by an electronic Hamiltonian containing excitation (site) energies of BChl a pigments and electronic couplings between them. Several such Hamiltonians have been developed in the past based on the information from various spectroscopic measurements of FMO; however, fine details of the excitonic structure and energy transfer in FMO, especially assignments of short-lived high-energy sites, remain elusive. Utilizing polarizable embedding quantum mechanics/molecular mechanics with the effective fragment potentials, we computed the electronic Hamiltonian of FMO that is in general agreement with previously reported empirical Hamiltonians and quantitatively reproduces experimental absorption and circular dichroism spectra of the FMO protein. The developed computational protocol is sufficiently simple and can be utilized for predictive modeling of other wild-type and mutated photosynthetic pigment-protein complexes.
Collapse
Affiliation(s)
- Yongbin Kim
- Department of Chemistry, Purdue University, 560 Oval Drive, West Lafayette, Indiana 47907, United States
| | - Dmitry Morozov
- Nanoscience Center and Department of Chemistry, University of Jyväskylä, P.O. Box 35, Jyväskylä 40014, Finland
| | - Valentyn Stadnytskyi
- Department of Physics and Astronomy, Purdue University, 525 Northwestern Avenue, West Lafayette, Indiana 47907, United States
- Laboratory of Chemical Physics, National Institute of Diabetes, Digestion and Kidney Diseases, National Institutes of Health, 5 Memorial Drive, Bethesda, Maryland 20892, United States
| | - Sergei Savikhin
- Department of Physics and Astronomy, Purdue University, 525 Northwestern Avenue, West Lafayette, Indiana 47907, United States
| | - Lyudmila V Slipchenko
- Department of Chemistry, Purdue University, 560 Oval Drive, West Lafayette, Indiana 47907, United States
| |
Collapse
|
33
|
Li W, Ren J, Shuai Z. Numerical assessment for accuracy and GPU acceleration of TD-DMRG time evolution schemes. J Chem Phys 2020; 152:024127. [PMID: 31941314 DOI: 10.1063/1.5135363] [Citation(s) in RCA: 30] [Impact Index Per Article: 7.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/23/2022] Open
Abstract
The time dependent density matrix renormalization group (TD-DMRG) has become one of the cutting edge methods of quantum dynamics for complex systems. In this paper, we comparatively study the accuracy of three time evolution schemes in the TD-DMRG, the global propagation and compression method with the Runge-Kutta algorithm (P&C-RK), the time dependent variational principle based methods with the matrix unfolding algorithm (TDVP-MU), and with the projector-splitting algorithm (TDVP-PS), by performing benchmarks on the exciton dynamics of the Fenna-Matthews-Olson complex. We show that TDVP-MU and TDVP-PS yield the same result when the time step size is converged and they are more accurate than P&C-RK4, while TDVP-PS tolerates a larger time step size than TDVP-MU. We further adopt the graphical processing units to accelerate the heavy tensor contractions in the TD-DMRG, and it is able to speed up the TDVP-MU and TDVP-PS schemes by up to 73 times.
Collapse
Affiliation(s)
- Weitang Li
- MOE Key Laboratory of Organic Optoelectronics and Molecular Engineering, Department of Chemistry, Tsinghua University, Beijing 100084, People's Republic of China
| | - Jiajun Ren
- MOE Key Laboratory of Organic Optoelectronics and Molecular Engineering, Department of Chemistry, Tsinghua University, Beijing 100084, People's Republic of China
| | - Zhigang Shuai
- MOE Key Laboratory of Organic Optoelectronics and Molecular Engineering, Department of Chemistry, Tsinghua University, Beijing 100084, People's Republic of China
| |
Collapse
|
34
|
Zheng J, Peng J, Xie Y, Long Y, Ning X, Lan Z. Study of the exciton dynamics in perylene bisimide (PBI) aggregates with symmetrical quasiclassical dynamics based on the Meyer–Miller mapping Hamiltonian. Phys Chem Chem Phys 2020; 22:18192-18204. [DOI: 10.1039/d0cp00648c] [Citation(s) in RCA: 9] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/16/2022]
Abstract
The exciton dynamics in one-dimensional stacked PBI (Perylene Bisimide) aggregates was studied with SQC-MM dynamics (Symmetrical Quasiclassical Dynamics based on the Meyer–Miller mapping Hamiltonian).
Collapse
Affiliation(s)
- Jie Zheng
- Industrial Research Institute of Nonwovens & Technical Textiles
- Shandong Center for Engineered Nonwovens (SCEN)
- College of Textiles Clothing
- Qingdao University
- Qingdao 266071
| | - 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
| | - Yunze Long
- Industrial Research Institute of Nonwovens & Technical Textiles
- Shandong Center for Engineered Nonwovens (SCEN)
- College of Textiles Clothing
- Qingdao University
- Qingdao 266071
| | - Xin Ning
- Industrial Research Institute of Nonwovens & Technical Textiles
- Shandong Center for Engineered Nonwovens (SCEN)
- College of Textiles Clothing
- Qingdao University
- Qingdao 266071
| | - 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
|
35
|
Larsson HR. Computing vibrational eigenstates with tree tensor network states (TTNS). J Chem Phys 2019; 151:204102. [DOI: 10.1063/1.5130390] [Citation(s) in RCA: 29] [Impact Index Per Article: 5.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/29/2023] Open
Affiliation(s)
- Henrik R. Larsson
- Division of Chemistry and Chemical Engineering, California Institute of Technology, Pasadena, California 91125, USA
| |
Collapse
|
36
|
Chin A, Mangaud E, Chevet V, Atabek O, Desouter-Lecomte M. Visualising the role of non-perturbative environment dynamics in the dissipative generation of coherent electronic motion. Chem Phys 2019. [DOI: 10.1016/j.chemphys.2019.110392] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/30/2022]
|
37
|
Popp W, Polkehn M, Hughes KH, Martinazzo R, Burghardt I. Vibronic coupling models for donor-acceptor aggregates using an effective-mode scheme: Application to mixed Frenkel and charge-transfer excitons in oligothiophene aggregates. J Chem Phys 2019; 150:244114. [DOI: 10.1063/1.5100529] [Citation(s) in RCA: 19] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/14/2022] Open
Affiliation(s)
- Wjatscheslaw Popp
- Institute of Physical and Theoretical Chemistry, Goethe University Frankfurt, Max-von-Laue-Str. 7, 60438 Frankfurt, Germany
| | - Matthias Polkehn
- Institute of Physical and Theoretical Chemistry, Goethe University Frankfurt, Max-von-Laue-Str. 7, 60438 Frankfurt, Germany
| | - Keith H. Hughes
- School of Natural Sciences, Bangor University, Bangor, Gwynedd LL572UW, United Kingdom
| | - Rocco Martinazzo
- Department of Chemistry, Università degli Studi di Milano, v. Golgi 19, 20133 Milano, Italy
| | - Irene Burghardt
- Institute of Physical and Theoretical Chemistry, Goethe University Frankfurt, Max-von-Laue-Str. 7, 60438 Frankfurt, Germany
| |
Collapse
|
38
|
Abstract
A comprehensive approach to modeling open quantum systems consistent with thermodynamics is presented. The theory of open quantum systems is employed to define system bath partitions. The Markovian master equation defines an isothermal partition between the system and bath. Two methods to derive the quantum master equation are described: the weak coupling limit and the repeated collision model. The role of the eigenoperators of the free system dynamics is highlighted, in particular, for driven systems. The thermodynamical relations are pointed out. Models that lead to loss of coherence, i.e., dephasing are described. The implication of the laws of thermodynamics to simulating transport and spectroscopy is described. The indications for self-averaging in large quantum systems and thus its importance in modeling are described. Basic modeling by the surrogate Hamiltonian is described, as well as thermal boundary conditions using the repeated collision model and their use in the stochastic surrogate Hamiltonian. The problem of modeling with explicitly time dependent driving is analyzed. Finally, the use of the stochastic surrogate Hamiltonian for modeling ultrafast spectroscopy and quantum control is reviewed.
Collapse
Affiliation(s)
- Ronnie Kosloff
- The Institute of Chemistry and The Fritz Haber Centre for Theoretical Chemistry, The Hebrew University of Jerusalem, Jerusalem 9190401, Israel
| |
Collapse
|
39
|
Baiardi A, Reiher M. Large-Scale Quantum Dynamics with Matrix Product States. J Chem Theory Comput 2019; 15:3481-3498. [DOI: 10.1021/acs.jctc.9b00301] [Citation(s) in RCA: 48] [Impact Index Per Article: 9.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/26/2022]
Affiliation(s)
- Alberto Baiardi
- Laboratorium für Physikalische Chemie, ETH Zürich, Vladimir-Prelog-Weg 2, 8093 Zürich, Switzerland
| | - Markus Reiher
- Laboratorium für Physikalische Chemie, ETH Zürich, Vladimir-Prelog-Weg 2, 8093 Zürich, Switzerland
| |
Collapse
|
40
|
Kaur R, Welsch R. Probing photodissociation dynamics using ring polymer molecular dynamics. J Chem Phys 2019; 150:114105. [PMID: 30901996 DOI: 10.1063/1.5086218] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Abstract
The performance of the ring polymer molecular dynamics (RPMD) approach to simulate typical photodissociation processes is assessed. The correct description of photodissociation requires the calculation of correlation functions or expectation values associated with non-equilibrium initial conditions, which was shown to be possible with RPMD very recently [J. Chem. Phys. 145, 204118 (2016)]. This approach is combined with treatment of the nonadiabatic dynamics employing the ring polymer surface hopping approach (RPSH), which is based on Tully's fewest switches surface hopping (FSSH) approach. The performance is tested using one-dimensional photodissociation models. It is found that RPSH with non-equilibrium initial conditions can well reproduce the time-dependent dissociation probability, and adiabatic and diabatic populations for cases where the crossing point is below and above the Franck-Condon point, respectively, while standard FSSH fails to reproduce the exact quantum dynamics in the latter case. Thus, it is shown that RPSH is an efficient and accurate alternative to standard FSSH, which is one of the most widely employed approaches to study photochemistry.
Collapse
Affiliation(s)
- Rajwant Kaur
- Center for Free-Electron Laser Science, DESY, Notkestrasse 85, 22607 Hamburg, Germany
| | - Ralph Welsch
- Center for Free-Electron Laser Science, DESY, Notkestrasse 85, 22607 Hamburg, Germany
| |
Collapse
|
41
|
Scheit S, Goswami S, Meyer HD, Köppel H. Fully quantal treatment of nonadiabatic molecular photodynamics: General considerations and application to the benzene cation. COMPUT THEOR CHEM 2019. [DOI: 10.1016/j.comptc.2019.01.011] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/25/2022]
|
42
|
Kim CW, Lee WG, Kim I, Rhee YM. Effect of Underdamped Vibration on Excitation Energy Transfer: Direct Comparison between Two Different Partitioning Schemes. J Phys Chem A 2019; 123:1186-1197. [DOI: 10.1021/acs.jpca.8b10977] [Citation(s) in RCA: 13] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/09/2023]
Affiliation(s)
- Chang Woo Kim
- Department of Chemistry, Korea Advanced Institute of Science and Technology (KAIST), Daejeon 34141, Korea
| | - Weon-Gyu Lee
- Department of Chemistry, Sungkyunkwan University, Suwon 16419, Korea
| | - Inkoo Kim
- Samsung Advanced Institute of Technology, Samsung Electronics, 130 Samsung-ro, Yeongtong-gu, Suwon 16678, Korea
| | - Young Min Rhee
- Department of Chemistry, Korea Advanced Institute of Science and Technology (KAIST), Daejeon 34141, Korea
| |
Collapse
|
43
|
Schulze J, Shibl MF, Al-Marri MJ, Kuhn O. High-dimensional exciton-vibrational wave-packet dynamics in the FMO complex. influence of site-specific spectral densities. EPJ WEB OF CONFERENCES 2019. [DOI: 10.1051/epjconf/201920510010] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
The correlated exciton-vibrational dynamics of the Fenna-Matthews-Olson (FMO) complex is studied using Multi-layer Multi-configuration Time-dependent Hartree (ML-MCTDH) wavepacket propagation. Exciton populations and coherences are shown to be sensitive to the details of the spectral density.
Collapse
|
44
|
Kurashige Y. Matrix product state formulation of the multiconfiguration time-dependent Hartree theory. J Chem Phys 2018; 149:194114. [DOI: 10.1063/1.5051498] [Citation(s) in RCA: 20] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Yuki Kurashige
- Department of Chemistry, Graduate School of Science, Kyoto University, Kitashirakawa Oiwake-cho, Sakyo-ku Kyoto 606-8502, Japan
| |
Collapse
|
45
|
Plötz PA, Megow J, Niehaus T, Kühn O. All-DFTB Approach to the Parametrization of the System-Bath Hamiltonian Describing Exciton-Vibrational Dynamics of Molecular Assemblies. J Chem Theory Comput 2018; 14:5001-5010. [PMID: 30141929 DOI: 10.1021/acs.jctc.8b00493] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/08/2023]
Abstract
Spectral density functions are central to the simulation of complex many body systems. Their determination requires making approximations not only to the dynamics but also to the underlying electronic structure theory. Here, blending different methods bears the danger of an inconsistent description. To solve this issue we propose an all-DFTB approach to determine spectral densities for the description of Frenkel excitons in molecular assemblies. The protocol is illustrated for a model of a PTCDI crystal, which involves the calculation of monomeric excitation energies and Coulomb couplings between monomer transitions, as well as their spectral distributions due to thermal fluctuations of the nuclei. Using dynamically defined normal modes, a mapping onto the standard harmonic oscillator spectral densities is achieved.
Collapse
Affiliation(s)
- Per-Arno Plötz
- Institut für Physik , Universität Rostock , Albert-Einstein-Strasse 23-24 , 18059 Rostock , Germany
| | - Jörg Megow
- Institut für Chemie , Universität Potsdam , Karl-Liebknecht-Strasse 24-25 , 14476 Potsdam , Germany
| | - Thomas Niehaus
- Univ Lyon, Université Claude Bernard Lyon 1, CNRS, Institut Lumière Matière, F-69622 Villeurbanne , France
| | - Oliver Kühn
- Institut für Physik , Universität Rostock , Albert-Einstein-Strasse 23-24 , 18059 Rostock , Germany
| |
Collapse
|
46
|
Affiliation(s)
- Haobin Wang
- Department of Chemistry, University of Colorado Denver, Denver, Colorado 80217-3364, USA and Beijing Computational Science Research Center, No. 10 East Xibeiwang Road, Haidian District, Beijing 100193, China
| | - Hans-Dieter Meyer
- Theoretische Chemie, Physikalisch-Chemisches Institut, Universität Heidelberg, Im Neuenheimer Feld 229, D-69120 Heidelberg, Germany
| |
Collapse
|
47
|
Wang Y, Ke Y, Zhao Y. The hierarchical and perturbative forms of stochastic Schrödinger equations and their applications to carrier dynamics in organic materials. WILEY INTERDISCIPLINARY REVIEWS-COMPUTATIONAL MOLECULAR SCIENCE 2018. [DOI: 10.1002/wcms.1375] [Citation(s) in RCA: 13] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 01/31/2023]
Affiliation(s)
- Yu‐Chen Wang
- State Key Laboratory for Physical Chemistry of Solid Surfaces, Collaborative Innovation Center of Chemistry for Energy Materials, Fujian Provincial Key Lab of Theoretical and Computational Chemistry, and College of Chemistry and Chemical Engineering, Xiamen University Xiamen China
| | - Yaling Ke
- State Key Laboratory for Physical Chemistry of Solid Surfaces, Collaborative Innovation Center of Chemistry for Energy Materials, Fujian Provincial Key Lab of Theoretical and Computational Chemistry, and College of Chemistry and Chemical Engineering, Xiamen University Xiamen China
| | - Yi Zhao
- State Key Laboratory for Physical Chemistry of Solid Surfaces, Collaborative Innovation Center of Chemistry for Energy Materials, Fujian Provincial Key Lab of Theoretical and Computational Chemistry, and College of Chemistry and Chemical Engineering, Xiamen University Xiamen China
| |
Collapse
|
48
|
Wang L, Fujihashi Y, Chen L, Zhao Y. Finite-temperature time-dependent variation with multiple Davydov states. J Chem Phys 2018; 146:124127. [PMID: 28388128 DOI: 10.1063/1.4979017] [Citation(s) in RCA: 42] [Impact Index Per Article: 7.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
The Dirac-Frenkel time-dependent variational approach with Davydov Ansätze is a sophisticated, yet efficient technique to obtain an accurate solution to many-body Schrödinger equations for energy and charge transfer dynamics in molecular aggregates and light-harvesting complexes. We extend this variational approach to finite temperature dynamics of the spin-boson model by adopting a Monte Carlo importance sampling method. In order to demonstrate the applicability of this approach, we compare calculated real-time quantum dynamics of the spin-boson model with that from numerically exact iterative quasiadiabatic propagator path integral (QUAPI) technique. The comparison shows that our variational approach with the single Davydov Ansätze is in excellent agreement with the QUAPI method at high temperatures, while the two differ at low temperatures. Accuracy in dynamics calculations employing a multitude of Davydov trial states is found to improve substantially over the single Davydov Ansatz, especially at low temperatures. At a moderate computational cost, our variational approach with the multiple Davydov Ansatz is shown to provide accurate spin-boson dynamics over a wide range of temperatures and bath spectral densities.
Collapse
Affiliation(s)
- Lu Wang
- Department of Physics, Zhejiang University, Hangzhou 310027, People's Republic of China
| | - Yuta Fujihashi
- Division of Materials Science, Nanyang Technological University, Singapore 639798, Singapore
| | - Lipeng Chen
- Division of Materials Science, Nanyang Technological University, Singapore 639798, Singapore
| | - Yang Zhao
- Division of Materials Science, Nanyang Technological University, Singapore 639798, Singapore
| |
Collapse
|
49
|
Reprint of “The effect of site-specific spectral densities on the high-dimensional exciton-vibrational dynamics in the FMO complex”. Chem Phys 2018. [DOI: 10.1016/j.chemphys.2018.05.015] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/19/2022]
|
50
|
Shi Q, Xu Y, Yan Y, Xu M. Efficient propagation of the hierarchical equations of motion using the matrix product state method. J Chem Phys 2018; 148:174102. [DOI: 10.1063/1.5026753] [Citation(s) in RCA: 41] [Impact Index Per Article: 6.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- 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
| | - Yang 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
| | - 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
| |
Collapse
|