1
|
Rey M, Carrington T. Using nested tensor train contracted basis functions with group theoretical techniques to compute (ro)-vibrational spectra of molecules with non-Abelian groups. J Chem Phys 2024; 161:044102. [PMID: 39037133 DOI: 10.1063/5.0219434] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/17/2024] [Accepted: 06/30/2024] [Indexed: 07/23/2024] Open
Abstract
In this paper, we use nested tensor-train contractions to compute vibrational and ro-vibrational energy levels of molecules with five and six atoms. At each step, we fully exploit symmetry by using symmetry adapted basis functions obtained from an irreducible tensor method. Contracted basis functions are determined by diagonalizing reduced dimensional Hamiltonian matrices. The size of matrices of eigenvectors, used to account for coupling between groups of coordinates, is reduced by discarding rows and columns. The size of the matrices that must be diagonalized is thus substantially reduced, making it possible to use direct eigensolvers, even for molecules with five and six atoms. The symmetry-adapted contracted vibrational basis functions have been used to compute J = 0 energy levels of the CH3CN (C3v) and J > 0 levels of CH4.
Collapse
Affiliation(s)
- Michaël Rey
- Groupe de Spectrométrie Moléculaire et Atmosphérique, UMR CNRS 7331, BP 1039, F-51687 Reims Cedex 2, France
| | - Tucker Carrington
- Chemistry Department, Queen's University, Kingston, Ontario K7L 3N6, Canada
| |
Collapse
|
2
|
Einsele R, Mitrić R. Nonadiabatic Exciton Dynamics and Energy Gradients in the Framework of FMO-LC-TDDFTB. J Chem Theory Comput 2024. [PMID: 39051619 DOI: 10.1021/acs.jctc.4c00539] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 07/27/2024]
Abstract
We introduce a novel methodology for simulating the excited-state dynamics of extensive molecular aggregates in the framework of the long-range corrected time-dependent density-functional tight-binding fragment molecular orbital method (FMO-LC-TDDFTB) combined with the mean-field Ehrenfest method. The electronic structure of the system is described in a quasi-diabatic basis composed of locally excited and charge-transfer states of all fragments. In order to carry out nonadiabatic molecular dynamics simulations, we derive and implement the excited-state gradients of the locally excited and charge-transfer states. Subsequently, the accuracy of the analytical excited-state gradients is evaluated. The applicability to the simulation of exciton transport in organic semiconductors is illustrated on a large cluster of anthracene molecules. Additionally, nonadiabatic molecular dynamics simulations of a model system of benzothieno-benzothiophene molecules highlight the method's utility in studying charge-transfer dynamics in organic materials. Our new methodology will facilitate the investigation of excitonic transfer in extensive biological systems, nanomaterials, and other complex molecular systems consisting of thousands of atoms.
Collapse
Affiliation(s)
- Richard Einsele
- Institut für Physikalische und Theoretische Chemie, Julius-Maximilians-Universität, Würzburg 97074, Germany
| | - Roland Mitrić
- Institut für Physikalische und Theoretische Chemie, Julius-Maximilians-Universität, Würzburg 97074, Germany
| |
Collapse
|
3
|
Liu S, Peng J, Bao P, Shi Q, Lan Z. Ultrafast Excited-State Energy Transfer in Phenylene Ethynylene Dendrimer: Quantum Dynamics with the Tensor Network Method. J Phys Chem A 2024. [PMID: 39047261 DOI: 10.1021/acs.jpca.4c00322] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 07/27/2024]
Abstract
Photoinduced excited-state energy transfer (EET) processes play an important role in solar energy conversions. Owing to their excellent photoharvesting and exciton-transport properties, phenylene ethynylene (PE) dendrimers display great potential for improving the efficiency of solar cells. In this work, we investigated the intramolecular EET dynamics in a dendrimer composed of two linear PE units (2-ring and 3-ring) using a fully quantum description based on the tensor network method. We first constructed a diabatic model Hamiltonian based on the electronic structure calculations. Using this diabatic vibronic coupling model, we tried to obtain the main features of the EET dynamics in terms of the several diabatic models with different numbers of vibrational modes (from 4 modes to 129 modes) and to explore the corresponding vibronic coupling interactions. The results show that the EET in this PE dendrimer is ultrafast. Four modes of A' symmetry play dominant roles in the dynamics; the remaining 86 modes of A' symmetry can dampen the electronic coherence; and the modes of A″ symmetry do not exhibit significant influence on the EET process. Overall, the first-order intrastate vibronic coupling terms show the dominant role in the EET dynamics, while the second-order intrastate vibronic coupling terms cause damping of the electronic coherence and slow down the overall EET process. This work provides a microscopic understanding of the EET dynamics in PE dendrimers.
Collapse
Affiliation(s)
- Sisi Liu
- SCNU Environmental Research Institute, Guangdong Provincial Key Laboratory of Chemical Pollution and Environmental Safety & MOE Key Laboratory of Environmental Theoretical Chemistry, South China Normal University, Guangzhou 510006, China
- School of Environment, South China Normal University, Guangzhou 510006, China
| | - Jiawei Peng
- SCNU Environmental Research Institute, Guangdong Provincial Key Laboratory of Chemical Pollution and Environmental Safety & MOE Key Laboratory of Environmental Theoretical Chemistry, South China Normal University, Guangzhou 510006, China
- School of Chemistry, South China Normal University, Guangzhou 510006, China
| | - Peng Bao
- 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, Beijing, Zhongguancun 100190, China
- 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, Beijing, Zhongguancun 100190, China
- University of Chinese Academy of Sciences, Beijing 100049, China
| | - Zhenggang Lan
- SCNU Environmental Research Institute, Guangdong Provincial Key Laboratory of Chemical Pollution and Environmental Safety & MOE Key Laboratory of Environmental Theoretical Chemistry, South China Normal University, Guangzhou 510006, China
- School of Environment, South China Normal University, Guangzhou 510006, China
| |
Collapse
|
4
|
Martínez-Fernández L, Green JA, Esposito L, Jouybari MY, Zhang Y, Santoro F, Kohler B, Improta R. The photoactivated dynamics of dGpdC and dCpdG sequences in DNA: a comprehensive quantum mechanical study. Chem Sci 2024; 15:9676-9693. [PMID: 38939156 PMCID: PMC11206432 DOI: 10.1039/d4sc00910j] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/06/2024] [Accepted: 05/04/2024] [Indexed: 06/29/2024] Open
Abstract
Study of alternating DNA GC sequences by different time-resolved spectroscopies has provided fundamental information on the interaction between UV light and DNA, a process of great biological importance. Multiple decay paths have been identified, but their interplay is still poorly understood. Here, we characterize the photophysics of GC-DNA by integrating different computational approaches, to study molecular models including up to 6 bases described at a full quantum mechanical level. Quantum dynamical simulations, exploiting a nonadiabatic linear vibronic coupling (LVC) model, coupled with molecular dynamics sampling of the initial structures of a (GC)5 DNA duplex, provide new insights into the photophysics in the sub-picosecond time-regime. They indicate a substantial population transfer, within 50 fs, from the spectroscopic states towards G → C charge transfer states involving two stacked bases (CTintra), thus explaining the ultrafast disappearance of fluorescence. This picture is consistent with that provided by quantum mechanical geometry optimizations, using time dependent-density functional theory and a polarizable continuum model, which we use to parametrize the LVC model and to map the main excited state deactivation pathways. For the first time, the infrared and excited state absorption signatures of the various states along these pathways are comprehensively mapped. The computational models suggest that the main deactivation pathways, which, according to experiment, lead to ground state recovery on the 10-50 ps time scale, involve CTintra followed by interstrand proton transfer from the neutral G to C-. Our calculations indicate that CTintra is populated to a larger extent and more rapidly in GC than in CG steps and suggest the likely involvement of monomer-like and interstrand charge transfer decay routes for isolated and less stacked CG steps. These findings underscore the importance of the DNA sequence and thermal fluctuations for the dynamics. They will also aid the interpretation of experimental results on other sequences.
Collapse
Affiliation(s)
- Lara Martínez-Fernández
- Departamento de Química Física de Materiales, Instituto de Química Física Blas Cabrera, CSIC 28006 Madrid Spain
| | - James Alexander Green
- Institut für Physikalische Theoretische Chemie, Goethe-Universität Frankfurt am Main Frankfurt am Main Germany
| | - Luciana Esposito
- Istituto di Biostrutture e Bioimmagini-CNR (IBB-CNR) Via De Amicis 95 I-80145 Napoli Italy
| | - Martha Yaghoubi Jouybari
- Department of Chemistry and Biomolecular Sciences, University of Ottawa 10 Marie Curie Ottawa Ontario K1N 6N5 Canada
- National Research Council of Canada 100 Sussex Drive Ottawa Ontario K1A 0R6 Canada
- Istituto di Chimica dei Composti Organometallici (ICCOM-CNR) Area della Ricerca del CNR, Via Moruzzi 1 I-56124 Pisa Italy
| | - Yuyuan Zhang
- Department of Chemistry and Biochemistry, The Ohio State University 100 West 18th Avenue Columbus Ohio 43210 USA
| | - Fabrizio Santoro
- Istituto di Chimica dei Composti Organometallici (ICCOM-CNR) Area della Ricerca del CNR, Via Moruzzi 1 I-56124 Pisa Italy
| | - Bern Kohler
- Department of Chemistry and Biochemistry, The Ohio State University 100 West 18th Avenue Columbus Ohio 43210 USA
| | - Roberto Improta
- Istituto di Biostrutture e Bioimmagini-CNR (IBB-CNR) Via De Amicis 95 I-80145 Napoli Italy
| |
Collapse
|
5
|
Kim CW, Franco I. General framework for quantifying dissipation pathways in open quantum systems. I. Theoretical formulation. J Chem Phys 2024; 160:214111. [PMID: 38833366 DOI: 10.1063/5.0202860] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/07/2024] [Accepted: 05/13/2024] [Indexed: 06/06/2024] Open
Abstract
We present a general and practical theoretical framework to investigate how energy is dissipated in open quantum system dynamics. This is performed by quantifying the contributions of individual bath components to the overall dissipation of the system. The framework is based on the Nakajima-Zwanzig projection operator technique, which allows us to express the rate of energy dissipation into a specific bath degree of freedom by using traces of operator products. The approach captures system-bath interactions to all orders, but is based on second-order perturbation theory on the off-diagonal subsystem's couplings and a Markovian description of the bath. The usefulness of our theory is demonstrated by applying it to various models of open quantum systems involving harmonic oscillators or spin baths and connecting the outcomes to existing results such as our previously reported formula derived for locally coupled harmonic baths [Kim and Franco, J. Chem. Phys. 154, 084109 (2021)]. We also prove that the dissipation calculated by our theory rigorously satisfies thermodynamic principles such as energy conservation and detailed balance. Overall, the strategy can be used to develop the theory and simulation of dissipation pathways to interpret and engineer the dynamics of open quantum systems.
Collapse
Affiliation(s)
- Chang Woo Kim
- Department of Chemistry, Chonnam National University, Gwangju 61186, South Korea
| | - 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
|
6
|
Zhou Y, Lu Y, Zhang Z, Zhang DH. An efficient and universal parallel algorithm for high-dimensional quantum dynamics in poly-atomic reactions. J Chem Phys 2024; 160:202502. [PMID: 38785279 DOI: 10.1063/5.0209245] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/19/2024] [Accepted: 05/06/2024] [Indexed: 05/25/2024] Open
Abstract
This study presents a parallel algorithm for high-dimensional quantum dynamics simulations in poly atomic reactions, integrating distributed- and shared-memory models. The distributions of the wave function and potential energy matrix across message passing interface processes are based on bundled radial and angular dimensions, with implementations featuring either two- or one-sided communication schemes. Using realistic parameters for the H + NH3 reaction, performance assessment reveals linear scalability, exceeding 90% efficiency with up to 600 processors. In addition, owing to the universal and concise structure, the algorithm demonstrates remarkable extensibility to diverse reaction systems, as demonstrated by successes with six-atom and four-atom reactions. This work establishes a robust foundation for high-dimensional dynamics studies, showcasing the algorithm's efficiency, scalability, and adaptability. The algorithm's potential as a valuable tool for unraveling quantum dynamics complexities is underscored, paving the way for future advancements in the field.
Collapse
Affiliation(s)
- Yong Zhou
- Anhui Province Key Laboratory for Control and Applications of Optoelectronic Information Materials, Department of Physics, Anhui Normal University, Wuhu 241000, People's Republic of China
- State Key Laboratory of Molecular Reaction Dynamics and Center for Theoretical and Computational Chemistry, Dalian Institute of Chemical Physics, Chinese Academy of Sciences, Dalian 116023, People's Republic of China
| | - Yunpeng Lu
- School of Chemistry, Chemical Engineering and Biotechnology, Nanyang Technological University, Singapore 637371, Singapore
| | - Zhaojun Zhang
- State Key Laboratory of Molecular Reaction Dynamics and Center for Theoretical and Computational Chemistry, Dalian Institute of Chemical Physics, Chinese Academy of Sciences, Dalian 116023, People's Republic of China
| | - Dong H Zhang
- State Key Laboratory of Molecular Reaction Dynamics and Center for Theoretical and Computational Chemistry, Dalian Institute of Chemical Physics, Chinese Academy of Sciences, Dalian 116023, People's Republic of China
- Hefei National Laboratory, Hefei 230088, People's Republic of China
| |
Collapse
|
7
|
Reddy SR, Coto PB, Thoss M. Intramolecular singlet fission: Quantum dynamical simulations including the effect of the laser field. J Chem Phys 2024; 160:194306. [PMID: 38767260 DOI: 10.1063/5.0209546] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/21/2024] [Accepted: 04/30/2024] [Indexed: 05/22/2024] Open
Abstract
In the previous work [Reddy et al., J. Chem. Phys. 151, 044307 (2019)], we have analyzed the dynamics of the intramolecular singlet fission process in a series of prototypical pentacene-based dimers, where the pentacene monomers are covalently bonded to a phenylene linker in ortho, meta, and para positions. The results obtained were qualitatively consistent with the experimental data available, showing an ultrafast population of the multiexcitonic state that mainly takes place via a mediated (superexchange-like) mechanism involving charge transfer and doubly excited states. Our results also highlighted the instrumental role of molecular vibrations in the process as a sizable population of the multiexcitonic state could only be obtained through vibronic coupling. Here, we extend these studies and investigate the effect of the laser field on the dynamics of intramolecular singlet fission by explicitly including the coupling to the laser field in our model. In this manner, and by selectively tuning the laser field to the different low-lying absorption bands of the systems investigated, we analyze the wavelength dependence of the intramolecular singlet fission process. In addition, we have also analyzed how the nature of the initially photoexcited electronic state (either localized or delocalized) affects its dynamics. Altogether, our results provide new insights into the design of intramolecular singlet fission-active molecules.
Collapse
Affiliation(s)
- S Rajagopala Reddy
- Department of Chemistry, School of Chemical Sciences and Pharmacy, Central University of Rajasthan, NH-8, Bandarsindri, Ajmer, Rajasthan 305817, India
| | - Pedro B Coto
- Materials Physics Center (CFM), Spanish National Research Council (CSIC) and Donostia International Physics Center (DIPC), Paseo Manuel de Lardizabal 5, 20018 Donostia-San Sebastián, Spain
| | - Michael Thoss
- Institute of Physics, Albert-Ludwigs University Freiburg, Hermann-Herder-Str. 3, 79104 Freiburg, Germany
| |
Collapse
|
8
|
García-González F, Otero JC, Ávila Ferrer FJ, Santoro F, Aranda D. Linear Vibronic Coupling Approach for Surface-Enhanced Raman Scattering: Quantifying the Charge-Transfer Enhancement Mechanism. J Chem Theory Comput 2024; 20:3850-3863. [PMID: 38687961 PMCID: PMC11099975 DOI: 10.1021/acs.jctc.4c00061] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/17/2024] [Revised: 04/12/2024] [Accepted: 04/15/2024] [Indexed: 05/02/2024]
Abstract
The outstanding amplification observed in surface-enhanced Raman scattering (SERS) is due to several enhancement mechanisms, and standing out among them are the plasmonic (PL) and charge-transfer (CT) mechanisms. The theoretical estimation of the enhancement factors of the CT mechanism is challenging because the excited-state coupling between bright plasmons and dark CT states must be properly introduced into the model to obtain reliable intensities. In this work, we aim at simulating electrochemical SERS spectra, considering models of pyridine on silver clusters subjected to an external electric field E⃗ that represents the effect of an electrode potential Vel. The method adopts quantum dynamical propagations of nuclear wavepackets on the coupled PL and CT states described with linear vibronic coupling models parametrized for each E⃗ through a fragment-based maximum-overlap diabatization. By presenting results at different values of E⃗, we show that indeed there is a relation between the population transferred to the CT states and the total scattered intensity. The tuning and detuning processes of the CT states with the bright PLs as a function of the electric field are in good agreement with those observed in experiments. Finally, our estimations for the CT enhancement factors predict values in the order of 105 to 106, meaning that when the CT and PL states are both in resonance with the excitation wavelength, the CT and PL enhancements are comparable, and vibrational bands whose intensity is amplified by different mechanisms can be observed together, in agreement with what was measured by typical experiments on silver electrodes.
Collapse
Affiliation(s)
- Francisco García-González
- Andalucía
Tech, Facultad de Ciencias, Departamento de Química Física, Universidad de Málaga, 29071 Málaga, Spain
| | - Juan Carlos Otero
- Andalucía
Tech, Facultad de Ciencias, Departamento de Química Física, Universidad de Málaga, 29071 Málaga, Spain
| | - Francisco J. Ávila Ferrer
- Andalucía
Tech, Facultad de Ciencias, Departamento de Química Física, Universidad de Málaga, 29071 Málaga, Spain
| | - Fabrizio Santoro
- Istituto
di Chimica dei Composti Organometallici (ICCOM-CNR), Area della Ricerca
del CNR, Via Moruzzi 1, I-56124 Pisa, Italy
| | - Daniel Aranda
- Andalucía
Tech, Facultad de Ciencias, Departamento de Química Física, Universidad de Málaga, 29071 Málaga, Spain
| |
Collapse
|
9
|
Hino K, Kurashige Y. Encoding a Many-Body Potential Energy Surface into a Grid-Based Matrix Product Operator. J Chem Theory Comput 2024; 20:3839-3849. [PMID: 38647101 DOI: 10.1021/acs.jctc.4c00046] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 04/25/2024]
Abstract
An efficient algorithm for compressing a given many-body potential energy surface (PES) of molecular systems into a grid-based matrix product operator (MPO) is proposed. The PES is once represented by a full-dimensional or truncated many-body expansion form, which is obtained by ab initio calculations at each grid mesh point, and then all terms in the expansion are compressed and merged into a single MPO while maintaining the bond dimension of the MPO as small as possible. It was shown that the ab initio PES of the H2CO was compressed by more than 2 orders of magnitude in the size of the site operators without loss of accuracy. By the use of grid basis, the tensor rank of the site operators of the MPO is reduced from four to three due to the diagonal nature of the position-dependent operators on grid basis, which significantly reduces the computational cost of the tensor contractions required in the real and imaginary time evolution of the matrix product state (MPS) wave functions with the grid-based MPO (Grid-MPO) Hamiltonian. Similar to other grid-based methods, Grid-MPO is easily applicable to any kinds of potentials of molecular systems, such as analytical empirical model potentials expressed by position operators and ab initio potentials, if the values at the grid points are available. Using the Grid-MPO combined with the MPS, we calculated the time correlation function of the Eigen cation H 3 O + ( H 2 O ) 3 to predict the infrared spectrum and compared with the experimental and the previous theoretical studies. The actual scaling with the size of systems was examined for the multidimensional Henon-Heiles Hamiltonian. It was shown that the method is considerably accelerated by the graphic processing unit (GPU) because the sizes of site operators were kept small and all tensors were able to be stored on the VRAM of a GPU.
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
- FOREST, JST, Honcho 4-1-8, Kawaguchi, Saitama 332-0012, Japan
- CREST, JST, Honcho 4-1-8, Kawaguchi, Saitama 332-0012, Japan
| |
Collapse
|
10
|
Boettcher V, Hartmann R, Beyer K, Strunz WT. Dynamics of a strongly coupled quantum heat engine-Computing bath observables from the hierarchy of pure states. J Chem Phys 2024; 160:094108. [PMID: 38436445 DOI: 10.1063/5.0192075] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/18/2023] [Accepted: 02/07/2024] [Indexed: 03/05/2024] Open
Abstract
We present a fully quantum dynamical treatment of a quantum heat engine and its baths based on the Hierarchy of Pure States (HOPS), an exact and general method for open quantum system dynamics. We show how the change of the bath energy and the interaction energy can be determined within HOPS for arbitrary coupling strength and smooth time dependence of the modulation protocol. The dynamics of all energetic contributions during the operation can be carefully examined both in its initial transient phase and, also later, in its periodic steady state. A quantum Otto engine with a qubit as an inherently nonlinear work medium is studied in a regime where the energy associated with the interaction Hamiltonian plays an important role for the global energy balance and, thus, must not be neglected when calculating its power and efficiency. We confirm that the work required to drive the coupling with the baths sensitively depends on the speed of the modulation protocol. Remarkably, departing from the conventional scheme of well-separated phases by allowing for temporal overlap, we discover that one can even gain energy from the modulation of bath interactions. We visualize these various work contributions using the analog of state change diagrams of thermodynamic cycles. We offer a concise, full presentation of HOPS with its extension to bath observables, as it serves as a universal tool for the numerically exact description of general quantum dynamical (thermodynamic) scenarios far from the weak-coupling limit.
Collapse
Affiliation(s)
- Valentin Boettcher
- Institute of Theoretical Physics, TUD Dresden University of Technology, 01062 Dresden, Germany
- Department of Physics, McGill University, Montréal, Québec H3A 2T8, Canada
| | - Richard Hartmann
- Institute of Theoretical Physics, TUD Dresden University of Technology, 01062 Dresden, Germany
| | - Konstantin Beyer
- Institute of Theoretical Physics, TUD Dresden University of Technology, 01062 Dresden, Germany
- Department of Physics, Stevens Institute of Technology, Hoboken, New Jersey 07030, USA
| | - Walter T Strunz
- Institute of Theoretical Physics, TUD Dresden University of Technology, 01062 Dresden, Germany
| |
Collapse
|
11
|
Hou E, Sun K, Gelin MF, Zhao Y. Finite temperature dynamics of the Holstein-Tavis-Cummings model. J Chem Phys 2024; 160:084116. [PMID: 38421073 DOI: 10.1063/5.0193471] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/22/2023] [Accepted: 02/05/2024] [Indexed: 03/02/2024] Open
Abstract
By employing the numerically accurate multiple Davydov Ansatz (mDA) formalism in combination with the thermo-field dynamics (TFD) representation of quantum mechanics, we systematically explore the influence of three parameters-temperature, photonic-mode detuning, and qubit-phonon coupling-on population dynamics and absorption spectra of the Holstein-Tavis-Cummings (HTC) model. It is found that elevated qubit-phonon couplings and/or temperatures have a similar impact on all dynamic observables: they suppress the amplitudes of Rabi oscillations in photonic populations as well as broaden the peaks and decrease their intensities in the absorption spectra. Our results unequivocally demonstrate that the HTC dynamics is very sensitive to the concerted variation of the three aforementioned parameters, and this finding can be used for fine-tuning polaritonic transport. The developed mDA-TFD methodology can be efficiently applied for modeling, predicting, optimizing, and comprehensively understanding dynamic and spectroscopic responses of actual molecular systems in microcavities.
Collapse
Affiliation(s)
- Erqin Hou
- School of Science, Hangzhou Dianzi University, Hangzhou 310018, China
| | - Kewei Sun
- School of Science, Hangzhou Dianzi University, Hangzhou 310018, China
| | - Maxim F Gelin
- School of Science, Hangzhou Dianzi University, Hangzhou 310018, China
- School of Materials Science and Engineering, Nanyang Technological University, Singapore 639798, Singapore
| | - Yang Zhao
- School of Materials Science and Engineering, Nanyang Technological University, Singapore 639798, Singapore
| |
Collapse
|
12
|
Brey D, Burghardt I. Coherent Transient Localization Mechanism of Interchain Exciton Transport in Regioregular P3HT: A Quantum-Dynamical Study. J Phys Chem Lett 2024; 15:1836-1845. [PMID: 38334949 DOI: 10.1021/acs.jpclett.3c03290] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/10/2024]
Abstract
Transient localization has been proposed as a transport mechanism in organic materials, for both charge carriers and excitons. Here, we characterize a quantum coherent transient localization mechanism using full quantum simulations of an H-aggregated model system representative of regioregular poly(3-hexylthiophene) (rrP3HT). A Frenkel-Holstein Hamiltonian parametrized from first principles is considered, including local high-frequency modes and anharmonic, site-correlated interchain modes. Quantum-dynamical calculations are carried out using the Multi-Layer Multi-Configuration Time-Dependent Hartree (ML-MCTDH) method for a 13-site system with 195 vibrational modes, under periodic boundary conditions. It is shown that temporary localization of exciton polarons alternates with resonant transfer driven by interchain modes. While the transport process is mainly determined by exciton-polarons at the low-energy band edge, persistent coupling with the excitonic manifold is observed, giving rise to a nonadiabatic excitonic flux. This elementary transport mechanism remains preserved for limited static disorder and gives way to Anderson localization when the static disorder becomes dominant.
Collapse
Affiliation(s)
- Dominik Brey
- Institute of Physical and Theoretical Chemistry, Goethe University Frankfurt, Max-von-Laue-Str. 7, 60438 Frankfurt, Germany
| | - Irene Burghardt
- Institute of Physical and Theoretical Chemistry, Goethe University Frankfurt, Max-von-Laue-Str. 7, 60438 Frankfurt, Germany
| |
Collapse
|
13
|
Song Q, Zhang X, Miao Z, Meng Q. Construction of a Mode-Combination Hamiltonian under the Grid-Based Representation for the Quantum Dynamics of OH + HO 2 → O 2 + H 2O. J Chem Theory Comput 2024; 20:597-613. [PMID: 38199964 DOI: 10.1021/acs.jctc.3c01090] [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
In this work, a systematic construction framework on a mode-combination Hamiltonian operator of a typical polyatomic reaction, OH + HO2 → O2 + H2O, is developed. First, a set of Jacobi coordinates are employed to construct the kinetic energy operator (KEO) through the polyspherical approach ( Phys. Rep. 2009, 484, 169). Second, due to the multiconfigurational electronic structure of this system, a non-adiabatic potential energy surface (PES) is constructed where the first singlet and triplet states are involved with spin-orbital coupling. To improve the training database, the training set of random energy data was optimized through a popular iterative optimization approach with extensive trajectories. Here, we propose an automatic trajectory method, instead of the classical trajectory on a crude PES, where the gradients are directly computed by the present ab initio calculations. Third, on the basis of the training set, the potential function is directly constructed in the canonical polyadic decomposition (CPD) form ( J. Chem. Theory Comput. 2021, 17, 2702-2713) which is helpful in propagating the nuclear wave function under the grid-based representation. To do this, the Gaussian process regression (GPR) approach for building the CPD form, called the CPD-GPR method ( J. Phys. Chem. Lett. 2022, 13, 11128-11135) is adopted where we further revise CPD-GPR by introducing the mode-combination (mc) scheme leading to the present CPD-mc-GPR approach. Constructing the full-dimension non-adiabatic Hamiltonian operator with mode combination, as test calculations, the nuclear wave function is propagated to preliminarily compute the reactive probability of OH + HO2 → O2 + H2O where the reactants are prepared in vibrational ground states and in the first triplet electronic state.
Collapse
Affiliation(s)
- Qingfei Song
- Department of Chemistry, Northwestern Polytechnical University, West Youyi Road 127, Xi'an 710072, China
| | - Xingyu Zhang
- Department of Chemistry, Northwestern Polytechnical University, West Youyi Road 127, Xi'an 710072, China
| | - Zekai Miao
- Department of Chemistry, Northwestern Polytechnical University, West Youyi Road 127, Xi'an 710072, China
| | - Qingyong Meng
- Department of Chemistry, Northwestern Polytechnical University, West Youyi Road 127, Xi'an 710072, China
| |
Collapse
|
14
|
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
|
15
|
Jensen AB, Højlund MG, Zoccante A, Madsen NK, Christiansen O. Efficient time-dependent vibrational coupled cluster computations with time-dependent basis sets at the two-mode coupling level: Full and hybrid TDMVCC[2]. J Chem Phys 2023; 159:204106. [PMID: 38010335 DOI: 10.1063/5.0175506] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/07/2023] [Accepted: 11/05/2023] [Indexed: 11/29/2023] Open
Abstract
The computation of the nuclear quantum dynamics of molecules is challenging, requiring both accuracy and efficiency to be applicable to systems of interest. Recently, theories have been developed for employing time-dependent basis functions (denoted modals) with vibrational coupled cluster theory (TDMVCC). The TDMVCC method was introduced along with a pilot implementation, which illustrated good accuracy in benchmark computations. In this paper, we report an efficient implementation of TDMVCC, covering the case where the wave function and Hamiltonian contain up to two-mode couplings. After a careful regrouping of terms, the wave function can be propagated with a cubic computational scaling with respect to the number of degrees of freedom. We discuss the use of a restricted set of active one-mode basis functions for each mode, as well as two interesting limits: (i) the use of a full active basis where the variational modal determination amounts essentially to the variational determination of a time-dependent reference state for the cluster expansion; and (ii) the use of a single function as an active basis for some degrees of freedom. The latter case defines a hybrid TDMVCC/TDH (time-dependent Hartree) approach that can obtain even lower computational scaling. The resulting computational scaling for hybrid and full TDMVCC[2] is illustrated for polyaromatic hydrocarbons with up to 264 modes. Finally, computations on the internal vibrational redistribution of benzoic acid (39 modes) are used to show the faster convergence of TDMVCC/TDH hybrid computations towards TDMVCC compared to simple neglect of some degrees of freedom.
Collapse
Affiliation(s)
| | - Mads Greisen Højlund
- Department of Chemistry, Aarhus University, Langelandsgade 140, DK-8000 Aarhus C, Denmark
| | - Alberto Zoccante
- Dipartimento di Scienze e Innovazione Tecnologica, Università del Piemonte Orientale (UPO), Via T. Michel 11, 15100 Alessandria, Italy
| | - Niels Kristian Madsen
- Department of Chemistry, Aarhus University, Langelandsgade 140, DK-8000 Aarhus C, Denmark
| | - Ove Christiansen
- Department of Chemistry, Aarhus University, Langelandsgade 140, DK-8000 Aarhus C, Denmark
| |
Collapse
|
16
|
Van Haeften A, Ash C, Worth G. Propagating multi-dimensional density operators using the multi-layer-ρ multi-configurational time-dependent Hartree method. J Chem Phys 2023; 159:194114. [PMID: 37982483 DOI: 10.1063/5.0172956] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/19/2023] [Accepted: 10/19/2023] [Indexed: 11/21/2023] Open
Abstract
Solving the Liouville-von-Neumann equation using a density operator provides a more complete picture of dynamical quantum phenomena than by using a wavepacket and solving the Schrödinger equation. As density operators are not restricted to the description of pure states, they can treat both thermalized and open systems. In practice, however, they are rarely used to study molecular systems as the computational resources required are even more prohibitive than those needed for wavepacket dynamics. In this paper, we demonstrate the potential utility of a scheme based on the powerful multi-layer multi-configurational time-dependent Hartree algorithm for propagating multi-dimensional density operators. Studies of two systems using this method are presented at a range of temperatures and including up to 13 degrees of freedom. The first case is single proton transfer in salicylaldimine, while the second is double proton transfer in porphycene. A comparison is also made with the approach of using stochastic wavepackets.
Collapse
Affiliation(s)
- Alice Van Haeften
- Department of Chemistry, University College London, London WC1H 0AJ, UK
| | - Ceridwen Ash
- Department of Chemistry, University College London, London WC1H 0AJ, UK
| | - Graham Worth
- Department of Chemistry, University College London, London WC1H 0AJ, UK
| |
Collapse
|
17
|
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: 1] [Impact Index Per Article: 1.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
|
18
|
Nag P, Babu P M J, Vennapusa SR. Significance of Nonadiabatic Effects on Efficient Triplet Generation in Lumazines. J Phys Chem A 2023; 127:7739-7746. [PMID: 37677159 DOI: 10.1021/acs.jpca.3c04121] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 09/09/2023]
Abstract
The photophysics of lumazines leading to triplet formation and the effect of thionation are explored in the presence of near-degenerate electronic states. Wave packet simulations are performed on model potential energy surfaces to understand the nonadiabatic population transfer among close-lying excited states. Ultrafast population transfer among singlets opens up new intersystem crossing channels from the higher states. An increased spin-orbit coupling strength originating from thionation enhances intersystem crossing and populates the higher triplets first. The rapid internal conversion in the triplet manifold eventually brings the molecules to their respective low-lying long-lived triplet state.
Collapse
Affiliation(s)
- Probal Nag
- School of Chemistry, Indian Institute of Science Education and Research Thiruvananthapuram, Maruthamala PO, Vithura, Thiruvananthapuram, Kerala 695551, India
| | - Janaarthana Babu P M
- School of Chemistry, Indian Institute of Science Education and Research Thiruvananthapuram, Maruthamala PO, Vithura, Thiruvananthapuram, Kerala 695551, India
| | - Sivaranjana Reddy Vennapusa
- School of Chemistry, Indian Institute of Science Education and Research Thiruvananthapuram, Maruthamala PO, Vithura, Thiruvananthapuram, Kerala 695551, India
| |
Collapse
|
19
|
Gera T, Chen L, Eisfeld A, Reimers JR, Taffet EJ, Raccah DIGB. Simulating optical linear absorption for mesoscale molecular aggregates: An adaptive hierarchy of pure states approach. J Chem Phys 2023; 158:2887556. [PMID: 37125709 DOI: 10.1063/5.0141882] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/09/2023] [Accepted: 03/07/2023] [Indexed: 05/02/2023] Open
Abstract
In this paper, we present dyadic adaptive HOPS (DadHOPS), a new method for calculating linear absorption spectra for large molecular aggregates. This method combines the adaptive HOPS (adHOPS) framework, which uses locality to improve computational scaling, with the dyadic HOPS method previously developed to calculate linear and nonlinear spectroscopic signals. To construct a local representation of dyadic HOPS, we introduce an initial state decomposition that reconstructs the linear absorption spectra from a sum over locally excited initial conditions. We demonstrate the sum over initial conditions can be efficiently Monte Carlo sampled and that the corresponding calculations achieve size-invariant [i.e., O(1)] scaling for sufficiently large aggregates while trivially incorporating static disorder in the Hamiltonian. We present calculations on the photosystem I core complex to explore the behavior of the initial state decomposition in complex molecular aggregates as well as proof-of-concept DadHOPS calculations on an artificial molecular aggregate inspired by perylene bis-imide to demonstrate the size-invariance of the method.
Collapse
Affiliation(s)
- Tarun Gera
- Department of Chemistry, Southern Methodist University, P.O. Box, Dallas, Texas 750314, USA
| | - Lipeng Chen
- Max Planck Institute for the Physics of Complex Systems, Nöthnitzer Str. 38, Dresden, Germany
| | - Alexander Eisfeld
- Max Planck Institute for the Physics of Complex Systems, Nöthnitzer Str. 38, Dresden, Germany
| | - Jeffrey R Reimers
- International Centre for Quantum and Molecular Structures and the School of Physics, Shanghai University, 200444 Shanghai, China
- School of Mathematical and Physical Sciences, University of Technology Sydney, Sydney NSW 2007, Australia
| | - Elliot J Taffet
- Department of Chemistry, Southern Methodist University, P.O. Box, Dallas, Texas 750314, USA
| | - Doran I G B Raccah
- Department of Chemistry, Southern Methodist University, P.O. Box, Dallas, Texas 750314, USA
| |
Collapse
|
20
|
Mendive-Tapia D, Meyer HD, Vendrell O. Optimal Mode Combination in the Multiconfiguration Time-Dependent Hartree Method through Multivariate Statistics: Factor Analysis and Hierarchical Clustering. J Chem Theory Comput 2023; 19:1144-1156. [PMID: 36716214 DOI: 10.1021/acs.jctc.2c01089] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/31/2023]
Abstract
The multiconfiguration time-dependent Hartree (MCTDH) method and its multilayer extension (ML-MCTDH) are powerful algorithms for the efficient computation of nuclear quantum dynamics in high-dimensional systems. By providing time-dependent variational orbitals and an optimal choice of layered effective degrees of freedom, one is able to reduce the computational cost to an amenable number of configurations. However, choices related to selecting properly the mode grouping and tensor tree are strongly system dependent and, thus far, subjectively based on intuition and/or experience. Therefore, herein we detail a new protocol based on multivariate statistics─more specifically, factor analysis and hierarchical clustering─for a reliable and convenient guiding in the optimal design of such complex "system-of-systems" tensor-network decompositions. The advantages of employing the new algorithm and its applicability are tested on water and two floppy protonated water clusters with large amplitude motions.
Collapse
Affiliation(s)
- David Mendive-Tapia
- Theoretische Chemie, Universität Heidelberg, Im Neuenheimer Feld 229, D-69120Heidelberg, Germany
| | - Hans-Dieter Meyer
- Theoretische Chemie, Universität Heidelberg, Im Neuenheimer Feld 229, D-69120Heidelberg, Germany
| | - Oriol Vendrell
- Theoretische Chemie, Universität Heidelberg, Im Neuenheimer Feld 229, D-69120Heidelberg, Germany
| |
Collapse
|
21
|
Dutta S, Lode AUJ, Alon OE. Fragmentation and correlations in a rotating Bose-Einstein condensate undergoing breakup. Sci Rep 2023; 13:3343. [PMID: 36849498 PMCID: PMC9971194 DOI: 10.1038/s41598-023-29516-w] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/28/2022] [Accepted: 02/06/2023] [Indexed: 03/01/2023] Open
Abstract
The theoretical investigation of rotating Bose-Einstein condensates has mainly focused on the emergence of quantum vortex states and the condensed properties of such systems. In the present work, we concentrate on other facets by examining the impact of rotation on the ground state of weakly interacting bosons confined in anharmonic potentials computed both at the mean-field level and particularly at the many-body level of theory. For the many-body computations, we employ the well-established many-body method known as the multiconfigurational time-dependent Hartree method for bosons. We present how various degrees of fragmentation can be generated following the breakup of the ground state densities in anharmonic traps without ramping up a potential barrier for strong rotations. The breakup of the densities is found to be associated with the acquisition of angular momentum in the condensate due to the rotation. In addition to fragmentation, the presence of many-body correlations is examined by computing the variances of the many-particle position and momentum operators. For strong rotations, the many-body variances become smaller than their mean-field counterparts, and one even finds a scenario with opposite anisotropies of the mean-field and many-body variances. Further, it is observed that for higher discrete symmetric systems of order k, namely three-fold and four-fold symmetry, breakup to k sub-clouds and emergence of k-fold fragmentation take place. All in all, we provide a thorough many-body investigation of how and which correlations build up when a trapped Bose-Einstein condensate breaks up under rotation.
Collapse
Affiliation(s)
- Sunayana Dutta
- Department of Physics, University of Haifa, 3498838, Haifa, Israel.
- Haifa Research Center for Theoretical Physics and Astrophysics, University of Haifa, 3498838, Haifa, Israel.
| | - Axel U J Lode
- Institute of Physics, Albert-Ludwig University of Freiburg, Hermann-Herder-Strasse 3, 79104, Freiburg, Germany
| | - Ofir E Alon
- Department of Physics, University of Haifa, 3498838, Haifa, Israel
- Haifa Research Center for Theoretical Physics and Astrophysics, University of Haifa, 3498838, Haifa, Israel
| |
Collapse
|
22
|
Lassmann Y, Hollas D, Curchod BFE. Extending the Applicability of the Multiple-Spawning Framework for Nonadiabatic Molecular Dynamics. J Phys Chem Lett 2022; 13:12011-12018. [PMID: 36541684 PMCID: PMC9806853 DOI: 10.1021/acs.jpclett.2c03295] [Citation(s) in RCA: 6] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/31/2022] [Accepted: 12/07/2022] [Indexed: 06/17/2023]
Abstract
Ab initio multiple-spawning (AIMS) describes the nonadiabatic dynamics of molecules by expanding nuclear wave functions in a basis of traveling multidimensional Gaussians called trajectory basis functions (TBFs). New TBFs can be spawned whenever nuclear amplitude is transferred between electronic states due to nonadiabatic transitions. While the adaptive size of the TBF basis grants AIMS its characteristic accuracy in describing nonadiabatic processes, it also leads to a fast and uncontrolled growth of the number of TBFs, penalizing computational efficiency. A different flavor of AIMS, called AIMS with informed stochastic selections (AIMSWISS), has recently been proposed to reduce the number of TBFs dramatically. Herein, we test the performance of AIMSWISS for a series of challenging nonadiabatic processes─photodynamics of two-dimensional model systems, 1,2-dithiane and chromium (0) hexacarbonyl─and show that this method is robust and extends the range of molecular systems that can be simulated within the multiple-spawning framework.
Collapse
|
23
|
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
|
24
|
Asha H, Green JA, Esposito L, Martinez-Fernandez L, Santoro F, Improta R. Effect of the Thermal Fluctuations of the Photophysics of GC and CG DNA Steps: A Computational Dynamical Study. J Phys Chem B 2022; 126:10608-10621. [PMID: 36508709 DOI: 10.1021/acs.jpcb.2c05688] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/14/2022]
Abstract
Here we refine and assess two computational procedures aimed to include the effect of thermal fluctuations on the electronic spectra and the ultrafast excited state dynamics of multichromophore systems, focusing on DNA duplexes. Our approach is based on a fragment diabatization procedure that, from a given Quantum Mechanical (QM) reference method, can provide the parameters (energy and coupling) of the reference diabatic states on the basis of the isolated fragments, either for a purely electronic excitonic Hamiltonian (FrDEx) or a linear vibronic coupling Hamiltonian (FrD-LVC). After having defined the most cost-effective procedure for DNA duplexes on two smaller fragments, FrDEx is used to simulate the absorption and Electronic Circular Dichroism (ECD) spectra of (GC)5 sequences, including the coupling with the Charge Transfer (CT) states, on a number of structures extracted from classical Molecular Dynamics (MD) simulations. The computed spectra are close to the reference TD-DFT calculations and fully consistent with the experimental ones. We then couple MD simulations and FrD-LVC to simulate the interplay between local excitations and CT transitions, both intrastrand and interstrand, in GC and CG steps when included in a oligoGC or in oligoAT DNA sequence. We predict that for both sequences a substantial part of the photoexcited population on G and C decays, within 50-100 fs, to the corresponding intrastrand CT states. This transfer is more effective for GC steps that, on average, are more closely stacked than CG ones.
Collapse
Affiliation(s)
- Haritha Asha
- Consiglio Nazionale delle Ricerche, Istituto di Biostrutture e Bioimmagini (IBB-CNR), Via De Amicis 95,I-80145Napoli, Italy
| | - James A Green
- Consiglio Nazionale delle Ricerche, Istituto di Biostrutture e Bioimmagini (IBB-CNR), Via De Amicis 95,I-80145Napoli, Italy.,Institute of Physical and Theoretical Chemistry, Goethe University Frankfurt, Max-von-Laue-Strasse 7, 60438Frankfurt am Main, Germany
| | - Luciana Esposito
- Consiglio Nazionale delle Ricerche, Istituto di Biostrutture e Bioimmagini (IBB-CNR), Via De Amicis 95,I-80145Napoli, Italy
| | - Lara Martinez-Fernandez
- Departamento de Química, Facultad de Ciencias and Institute for Advanced Research in Chemistry (IADCHEM), Universidad Autònoma de Madrid, Campus de Excelencia UAM-CSIC, 28049Madrid, Spain
| | - Fabrizio Santoro
- Consiglio Nazionale delle Ricerche, Istituto di Chimica dei Composti Organo Metallici (ICCOM-CNR), SS di Pisa, Area della Ricerca, via G. Moruzzi 1, I-56124Pisa, Italy
| | - Roberto Improta
- Consiglio Nazionale delle Ricerche, Istituto di Biostrutture e Bioimmagini (IBB-CNR), Via De Amicis 95,I-80145Napoli, Italy.,DTU Chemistry, Technical University of Denmark, Kemitorvet Bldg 207, DK-2800Kongens Lyngby, Denmark
| |
Collapse
|
25
|
The coupling of the hydrated proton to its first solvation shell. Nat Commun 2022; 13:6170. [PMID: 36257946 PMCID: PMC9579203 DOI: 10.1038/s41467-022-33650-w] [Citation(s) in RCA: 8] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/22/2022] [Accepted: 09/28/2022] [Indexed: 11/30/2022] Open
Abstract
The Zundel (\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$${H}_{5}{O}_{2}^{+}$$\end{document}H5O2+) and Eigen (\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$${H}_{9}{O}_{4}^{+}$$\end{document}H9O4+) cations play an important role as intermediate structures for proton transfer processes in liquid water. In the gas phase they exhibit radically different infrared (IR) spectra. The question arises: is there a least common denominator structure that explains the IR spectra of both, the Zundel and Eigen cations, and hence of the solvated proton? Full dimensional quantum simulations of these protonated cations demonstrate that two dynamical water molecules and an excess proton constitute this fundamental subunit. Embedded in the static environment of the parent Eigen cation, this subunit reproduces the positions and broadenings of its main excess-proton bands. In isolation, its spectrum reverts to the well-known Zundel ion. Hence, the dynamics of this subunit polarized by an environment suffice to explain the spectral signatures and anharmonic couplings of the solvated proton in its first solvation shell. The Zundel [H(H2O)2]+ and Eigen [H(H2O)4]+ cations exhibit radicallly different infrared spectra and are the limiting dynamical structures involved in proton mobility in liquid water. Here, the authors find through quantum dynamics simulations that two polarized water molecules and a proton suffice to explain the key spectroscopic features connected to proton mobility for both species.
Collapse
|
26
|
Green JA, Gómez S, Worth G, Santoro F, Improta R. Solvent Effects on Ultrafast Charge Transfer Population: Insights from the Quantum Dynamics of Guanine-Cytosine in Chloroform. Chemistry 2022; 28:e202201731. [PMID: 35950519 PMCID: PMC9828530 DOI: 10.1002/chem.202201731] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/06/2022] [Indexed: 01/12/2023]
Abstract
We study the ultrafast photoactivated dynamics of the hydrogen bonded dimer Guanine-Cytosine in chloroform solution, focusing on the population of the Guanine→Cytosine charge transfer state (GC-CT), an important elementary process for the photophysics and photochemistry of nucleic acids. We integrate a quantum dynamics propagation scheme, based on a linear vibronic model parameterized through time dependent density functional theory calculations, with four different solvation models, either implicit or explicit. On average, after 50 fs, 30∼40 % of the bright excited state population has been transferred to GC-CT. This process is thus fast and effective, especially when transferring from the Guanine bright excited states, in line with the available experimental studies. Independent of the adopted solvation model, the population of GC-CT is however disfavoured in solution with respect to the gas phase. We show that dynamical solvation effects are responsible for this puzzling result and assess the different chemical-physical effects modulating the population of CT states on the ultrafast time-scale. We also propose some simple analyses to predict how solvent can affect the population transfer between bright and CT states, showing that the effect of the solute/solvent electrostatic interactions on the energy of the CT state can provide a rather reliable indication of its possible population.
Collapse
Affiliation(s)
- James A. Green
- Istituto di Biostrutture e Bioimmagini-CNRVia De Amicis 95I-80145Napoli
| | - Sandra Gómez
- Department of ChemistryUniversity College London20 Gordon StreetLondonWC1H 0AJUnited Kingdom,Departamento de Química FísicaUniversity of SalamancaSalamanca37008Spain
| | - Graham Worth
- Department of ChemistryUniversity College London20 Gordon StreetLondonWC1H 0AJUnited Kingdom
| | - Fabrizio Santoro
- Istituto di Chimica die Composti Organometallici (ICCOM-CNR), Area della Ricerca del CNRVia Moruzzi 1I-56124Pisa
| | - Roberto Improta
- Istituto di Biostrutture e Bioimmagini-CNRVia De Amicis 95I-80145Napoli
| |
Collapse
|
27
|
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
|
28
|
Song Q, Zhang X, Gatti F, Miao Z, Zhang Q, Meng Q. Multilayer Multiconfiguration Time-Dependent Hartree Study on the Mode-/Bond-Specific Quantum Dynamics of Water Dissociation on Cu(111). J Phys Chem A 2022; 126:6047-6058. [PMID: 36054932 DOI: 10.1021/acs.jpca.2c03092] [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/28/2022]
Abstract
In this work, full-dimensional (9D) quantum dynamics calculations on mode-/bond-specific surface scattering of a water molecule on a copper (111) rigid surface are performed through the multilayer multiconfiguration time-dependent Hartree (ML-MCTDH) method. To easily perform the ML-MCTDH calculations on such a triatomic molecule-surface system, we first choose specific Jacobi coordinates as a set of coordinates of water. Next, to efficiently perform the 9D ML-MCTDH wavepacket propagation, the potential energy surface is transferred to a canonical polyadic decomposition form with the aid of a Monte Carlo-based method. Excitation-specific dissociation probabilities of H2O on Cu(111) are computed, and mode-/bond-specific dynamics are demonstrated by comparison with a probability curve computed for a water molecule in the ground state. The dependence of the dissociation probability of the initial state of H2O is studied, and it is found that the excitation-specific dissociation probabilities can be divided into three groups. We find that the vibrationally excited states enhance the dissociation reactivity of H2O, while the rotationally excited states hardly influence it.
Collapse
Affiliation(s)
- Qingfei Song
- Department of Chemistry, Northwestern Polytechnical University, West Youyi Road 127, 710072 Xi'an, China.,Institut des Sciences Moléculaires d'Orsay, CNRS-UMR 8214, Université Paris-Saclay, Bâtiment 520, F-91405 Orsay, France
| | - Xingyu Zhang
- Department of Chemistry, Northwestern Polytechnical University, West Youyi Road 127, 710072 Xi'an, China
| | - Fabien Gatti
- Institut des Sciences Moléculaires d'Orsay, CNRS-UMR 8214, Université Paris-Saclay, Bâtiment 520, F-91405 Orsay, France
| | - Zekai Miao
- Department of Chemistry, Northwestern Polytechnical University, West Youyi Road 127, 710072 Xi'an, China
| | - Qiuyu Zhang
- Department of Chemistry, Northwestern Polytechnical University, West Youyi Road 127, 710072 Xi'an, China
| | - Qingyong Meng
- Department of Chemistry, Northwestern Polytechnical University, West Youyi Road 127, 710072 Xi'an, China
| |
Collapse
|
29
|
Peng WT, Brey D, Giannini S, Dell’Angelo D, Burghardt I, Blumberger J. Exciton Dissociation in a Model Organic Interface: Excitonic State-Based Surface Hopping versus Multiconfigurational Time-Dependent Hartree. J Phys Chem Lett 2022; 13:7105-7112. [PMID: 35900333 PMCID: PMC9376959 DOI: 10.1021/acs.jpclett.2c01928] [Citation(s) in RCA: 6] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/22/2022] [Accepted: 07/19/2022] [Indexed: 05/20/2023]
Abstract
Quantum dynamical simulations are essential for a molecular-level understanding of light-induced processes in optoelectronic materials, but they tend to be computationally demanding. We introduce an efficient mixed quantum-classical nonadiabatic molecular dynamics method termed eXcitonic state-based Surface Hopping (X-SH), which propagates the electronic Schrödinger equation in the space of local excitonic and charge-transfer electronic states, coupled to the thermal motion of the nuclear degrees of freedom. The method is applied to exciton decay in a 1D model of a fullerene-oligothiophene junction, and the results are compared to the ones from a fully quantum dynamical treatment at the level of the Multilayer Multiconfigurational Time-Dependent Hartree (ML-MCTDH) approach. Both methods predict that charge-separated states are formed on the 10-100 fs time scale via multiple "hot-exciton dissociation" pathways. The results demonstrate that X-SH is a promising tool advancing the simulation of photoexcited processes from the molecular to the true nanomaterials scale.
Collapse
Affiliation(s)
- Wei-Tao Peng
- Department
of Physics and Astronomy and Thomas Young Centre, University College London, London WC1E 6BT, United Kingdom
| | - Dominik Brey
- Institute
of Physical and Theoretical Chemistry, Goethe
University Frankfurt, Max-von-Laue-Strasse 7, 60438 Frankfurt am Main, Germany
| | - Samuele Giannini
- Department
of Physics and Astronomy and Thomas Young Centre, University College London, London WC1E 6BT, United Kingdom
| | - David Dell’Angelo
- Department
of Physics and Astronomy and Thomas Young Centre, University College London, London WC1E 6BT, United Kingdom
| | - Irene Burghardt
- Institute
of Physical and Theoretical Chemistry, Goethe
University Frankfurt, Max-von-Laue-Strasse 7, 60438 Frankfurt am Main, Germany
| | - Jochen Blumberger
- Department
of Physics and Astronomy and Thomas Young Centre, University College London, London WC1E 6BT, United Kingdom
| |
Collapse
|
30
|
Mondelo-Martell M, Brey D, Burghardt I. Quantum dynamical study of inter-chain exciton transport in a regioregular P3HT model system at finite temperature: HJ vs. H-aggregate models. J Chem Phys 2022; 157:094108. [DOI: 10.1063/5.0104729] [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 report on quantum dynamical simulations of inter-chain exciton transport in a model of regioregular poly(3-hexylthiophene), rr-P3HT, at finite temperature, using the Multi-Layer Multi-Configuration Time-Dependent Hartree (ML-MCTDH) method for a system of up to 63 electronic states and 180 vibrational modes. A Frenkel Hamiltonian of HJ aggregate type is used, along with a reduced H-aggregate representation; electron-phonon coupling includes local high-frequency modes as well as anharmonic intermolecular modes. The latter are operative in mediating inter-chain transport, by a mechanism of transient localization type. Strikingly, this mechanism is found to be of quantum coherent character and involves non-adiabatic effects. Using periodic boundary conditions, a normal diffusion regime is identified from the exciton mean-squared displacement, apart from early-time transients. Diffusion coefficients are found to be of the order of 3 x 10-3 cm2/s, showing a non-monotonous increase with temperature.
Collapse
Affiliation(s)
- Manel Mondelo-Martell
- Institut für Physikalische u. Theoretische Chemie, Goethe-Universitat Frankfurt am Main Institut fur Physikalische und Theoretische Chemie, Germany
| | | | - Irene Burghardt
- Institute of Physical and Theoretical Chemistry, Goethe University Frankfurt, Germany
| |
Collapse
|
31
|
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
|
32
|
Ten Brink M, Gräber S, Hopjan M, Jansen D, Stolpp J, Heidrich-Meisner F, Blöchl PE. Real-time non-adiabatic dynamics in the one-dimensional Holstein model: Trajectory-based vs exact methods. J Chem Phys 2022; 156:234109. [PMID: 35732530 DOI: 10.1063/5.0092063] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
We benchmark a set of quantum-chemistry methods, including multitrajectory Ehrenfest, fewest-switches surface-hopping, and multiconfigurational-Ehrenfest dynamics, against exact quantum-many-body techniques by studying real-time dynamics in the Holstein model. This is a paradigmatic model in condensed matter theory incorporating a local coupling of electrons to Einstein phonons. For the two-site and three-site Holstein model, we discuss the exact and quantum-chemistry methods in terms of the Born-Huang formalism, covering different initial states, which either start on a single Born-Oppenheimer surface, or with the electron localized to a single site. For extended systems with up to 51 sites, we address both the physics of single Holstein polarons and the dynamics of charge-density waves at finite electron densities. For these extended systems, we compare the quantum-chemistry methods to exact dynamics obtained from time-dependent density matrix renormalization group calculations with local basis optimization (DMRG-LBO). We observe that the multitrajectory Ehrenfest method, in general, only captures the ultrashort time dynamics accurately. In contrast, the surface-hopping method with suitable corrections provides a much better description of the long-time behavior but struggles with the short-time description of coherences between different Born-Oppenheimer states. We show that the multiconfigurational Ehrenfest method yields a significant improvement over the multitrajectory Ehrenfest method and can be converged to the exact results in small systems with moderate computational efforts. We further observe that for extended systems, this convergence is slower with respect to the number of configurations. Our benchmark study demonstrates that DMRG-LBO is a useful tool for assessing the quality of the quantum-chemistry methods.
Collapse
Affiliation(s)
- M Ten Brink
- Institut für Theoretische Physik, Georg-August-Universität Göttingen, Friedrich-Hund-Platz 1, 37077 Göttingen, Germany
| | - S Gräber
- Institut für Theoretische Physik, Georg-August-Universität Göttingen, Friedrich-Hund-Platz 1, 37077 Göttingen, Germany
| | - M Hopjan
- Institut für Theoretische Physik, Georg-August-Universität Göttingen, Friedrich-Hund-Platz 1, 37077 Göttingen, Germany
| | - D Jansen
- Institut für Theoretische Physik, Georg-August-Universität Göttingen, Friedrich-Hund-Platz 1, 37077 Göttingen, Germany
| | - J Stolpp
- Institut für Theoretische Physik, Georg-August-Universität Göttingen, Friedrich-Hund-Platz 1, 37077 Göttingen, Germany
| | - F Heidrich-Meisner
- Institut für Theoretische Physik, Georg-August-Universität Göttingen, Friedrich-Hund-Platz 1, 37077 Göttingen, Germany
| | - P E Blöchl
- Institut für Theoretische Physik, Georg-August-Universität Göttingen, Friedrich-Hund-Platz 1, 37077 Göttingen, Germany
| |
Collapse
|
33
|
Nag P, Isukapalli SVK, Nath A, Vennapusa SR. Revisiting the Dynamics of Triplet Formation in Anthraquinones. J Phys Chem A 2022; 126:3680-3687. [PMID: 35655394 DOI: 10.1021/acs.jpca.2c01980] [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
Triplet formation pathways in 9,10-anthraquinone (AQ) and its hydroxy derivative, 1-hydroxyanthraquinone (HAQ), are studied theoretically. Dynamics simulations on the model singlet-triplet potential energy surfaces within the linear vibronic coupling framework are performed to elucidate possible internal conversion (IC) and intersystem crossing (ISC) pathways in these molecules. An ultrafast IC decay from the "bright" S4 to S1 followed by efficient ISC via S1-T4 and S1-T5 pathways fosters a high triplet quantum yield (ΦT = 0.90) in AQ. In HAQ, a new nonradiative channel of "barrierless" excited-state intramolecular proton transfer (ESIPT) opens up and competes with the IC decay to S1 upon photoexcitation to the "bright" S2. Extremely fast ESIPT on S2 reduces the efficiency of triplet formation via possible ISC pathways involving S1 and S2, resulting in a low ΦT (=0.17).
Collapse
Affiliation(s)
- Probal Nag
- School of Chemistry, Indian Institute of Science Education and Research Thiruvananthapuram, Maruthamala PO, Vithuran, Thiruvananthapuram, Kerala 695551, India
| | - Sai Vamsi Krishna Isukapalli
- School of Chemistry, Indian Institute of Science Education and Research Thiruvananthapuram, Maruthamala PO, Vithuran, Thiruvananthapuram, Kerala 695551, India
| | - Akashdeep Nath
- School of Chemistry, Indian Institute of Science Education and Research Thiruvananthapuram, Maruthamala PO, Vithuran, Thiruvananthapuram, Kerala 695551, India
| | - Sivaranjana Reddy Vennapusa
- School of Chemistry, Indian Institute of Science Education and Research Thiruvananthapuram, Maruthamala PO, Vithuran, Thiruvananthapuram, Kerala 695551, India
| |
Collapse
|
34
|
Lyu N, Soley MB, Batista VS. Tensor-Train Split-Operator KSL (TT-SOKSL) Method for Quantum Dynamics Simulations. J Chem Theory Comput 2022; 18:3327-3346. [PMID: 35649210 DOI: 10.1021/acs.jctc.2c00209] [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/30/2022]
Abstract
Numerically exact simulations of quantum reaction dynamics, including nonadiabatic effects in excited electronic states, are essential to gain fundamental insights into ultrafast chemical reactivity and rigorous interpretations of molecular spectroscopy. Here, we introduce the tensor-train split-operator KSL (TT-SOKSL) method for quantum simulations in tensor-train (TT)/matrix product state (MPS) representations. TT-SOKSL propagates the quantum state as a tensor train using the Trotter expansion of the time-evolution operator, as in the tensor-train split-operator Fourier transform (TT-SOFT) method. However, the exponential operators of the Trotter expansion are applied using a rank-adaptive TT-KSL scheme instead of using the scaling and squaring approach as in TT-SOFT. We demonstrate the accuracy and efficiency of TT-SOKSL as applied to simulations of the photoisomerization of the retinal chromophore in rhodopsin, including nonadiabatic dynamics at a conical intersection of potential energy surfaces. The quantum evolution is described in full dimensionality by a time-dependent wavepacket evolving according to a two-state 25-dimensional model Hamiltonian. We find that TT-SOKSL converges faster than TT-SOFT with respect to the maximally allowed memory requirement of the tensor-train representation and better preserves the norm of the time-evolving state. When compared to the corresponding simulations based on the TT-KSL method, TT-SOKSL has the advantage of avoiding the need to construct the matrix product state Laplacian by exploiting the linear scaling of multidimensional tensor-train Fourier transforms.
Collapse
Affiliation(s)
- Ningyi Lyu
- Department of Chemistry, Yale University, P.O. Box 208107, New Haven, Connecticut 06520-8107, United States
| | - Micheline B Soley
- Department of Chemistry, Yale University, P.O. Box 208107, New Haven, Connecticut 06520-8107, United States.,Yale Quantum Institute, Yale University, P.O. Box 208334, New Haven, Connecticut 06520-8263, United States
| | - Victor S Batista
- Department of Chemistry, Yale University, P.O. Box 208107, New Haven, Connecticut 06520-8107, United States.,Yale Quantum Institute, Yale University, P.O. Box 208334, New Haven, Connecticut 06520-8263, United States
| |
Collapse
|
35
|
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
|
36
|
Fischer EW, Werther M, Bouakline F, Grossmann F, Saalfrank P. Non-Markovian Vibrational Relaxation Dynamics at Surfaces. J Chem Phys 2022; 156:214702. [DOI: 10.1063/5.0092836] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Vibrational dynamics of adsorbates near surfaces plays both an important role for applied surface science and as model lab for studying fundamental problems of open quantum systems. We employ a previously developed model for the relaxation of a D-Si-Si bending mode at a D:Si(100)-(2$\times$1) surface, induced by a ``bath' of more than $2000$ phonon modes [U. Lorenz, P. Saalfrank, Chem. Phys. {\bf 482}, 69 (2017)], to extend previous work along various directions. First, we use a Hierarchical Effective Mode (HEM) model [E.W. Fischer, F. Bouakline, M. Werther, P. Saalfrank, J. Chem. Phys. {\bf 153}, 064704 (2020)] to study relaxation of higher excited vibrational states than hitherto done, by solving a high-dimensional system-bath time-dependent Schr\"odinger equation (TDSE). In the HEM approach, (many) real bath modes are replaced by (much less) effective bath modes. Accordingly, we are able to examine scaling laws for vibrational relaxation lifetimes for a realistic surface science problem. Second, we compare the performance of the multilayer multiconfigurational time-dependent Hartree (ML-MCTDH) approach with the recently developed coherent-state based multi-Davydov D2 {\it ansatz} [N. Zhou, Z. Huang, J. Zhu, V. Chernyak, Y. Zhao, {J. Chem. Phys.} {\bf 143}, 014113 (2015)]. Both approaches work well, with some computational advantages for the latter in the presented context. Third, we apply open-system density matrix theory in comparison with basically ``exact' solutions of the multi-mode TDSEs. Specifically, we use an open-system Liouville-von Neumann (LvN) equation treating vibration-phonon coupling as Markovian dissipation in Lindblad form to quantify effects beyond the Born-Markov approximation.
Collapse
Affiliation(s)
| | | | - Foudhil Bouakline
- Institute of Chemistry, Universität Potsdam Institut für Chemie, Germany
| | - Frank Grossmann
- Institute for Theoretical Physics, Technische Universität Dresden Fachrichtung Physik, Germany
| | - Peter Saalfrank
- Institut für Chemie, Universität Potsdam Institut für Chemie, Germany
| |
Collapse
|
37
|
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
|
38
|
Fregoni J, Garcia-Vidal FJ, Feist J. Theoretical Challenges in Polaritonic Chemistry. ACS PHOTONICS 2022; 9:1096-1107. [PMID: 35480492 PMCID: PMC9026242 DOI: 10.1021/acsphotonics.1c01749] [Citation(s) in RCA: 42] [Impact Index Per Article: 21.0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 11/15/2021] [Revised: 02/01/2022] [Accepted: 02/02/2022] [Indexed: 06/14/2023]
Abstract
Polaritonic chemistry exploits strong light-matter coupling between molecules and confined electromagnetic field modes to enable new chemical reactivities. In systems displaying this functionality, the choice of the cavity determines both the confinement of the electromagnetic field and the number of molecules that are involved in the process. While in wavelength-scale optical cavities the light-matter interaction is ruled by collective effects, plasmonic subwavelength nanocavities allow even single molecules to reach strong coupling. Due to these very distinct situations, a multiscale theoretical toolbox is then required to explore the rich phenomenology of polaritonic chemistry. Within this framework, each component of the system (molecules and electromagnetic modes) needs to be treated in sufficient detail to obtain reliable results. Starting from the very general aspects of light-molecule interactions in typical experimental setups, we underline the basic concepts that should be taken into account when operating in this new area of research. Building on these considerations, we then provide a map of the theoretical tools already available to tackle chemical applications of molecular polaritons at different scales. Throughout the discussion, we draw attention to both the successes and the challenges still ahead in the theoretical description of polaritonic chemistry.
Collapse
|
39
|
Bera A, Nag P, Pandey D, Vennapusa SR. Theoretical approach to modeling the early nonadiabatic events of ESIPT originating from three-state conical intersection in quinophthalone. Photochem Photobiol Sci 2022; 21:1287-1298. [PMID: 35403961 DOI: 10.1007/s43630-022-00220-4] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/22/2021] [Accepted: 03/25/2022] [Indexed: 11/30/2022]
Abstract
We explore the excited-state intramolecular proton transfer process of quinophthalone theoretically. This molecule possesses three low-lying singlet excited states ([Formula: see text] and [Formula: see text]) in a narrow energy gap of less than the N-H stretching frequency. Dynamics simulations show nonadiabatic wavepacket transfer to [Formula: see text] and [Formula: see text] upon initiating the wavepacket on [Formula: see text]. Multiple accessible conical intersections that lie in the Franck-Condon region facilitate the nonadiabatic wavepacket transfer. Nuclear densities associated with the proton transfer promoting vibrations would start accumulating on [Formula: see text] and [Formula: see text] within a few tens of femtoseconds, validating the involvement of these vibrations in the nonadiabatic events that occur before the proton transfer process. Our findings emphasize the necessity of refined kinetic models for assigning the time constants of ultrafast transient spectroscopy measurements due to the simultaneous evolution of nonadiabatic events and proton transfer kinetics in quinophthalone.
Collapse
Affiliation(s)
- Anshuman Bera
- School of Chemistry, Indian Institute of Science Education and Research Thiruvananthapuram, Maruthamala PO, Vithura, Thiruvananthapuram, 695551, India
| | - Probal Nag
- School of Chemistry, Indian Institute of Science Education and Research Thiruvananthapuram, Maruthamala PO, Vithura, Thiruvananthapuram, 695551, India
| | - Diksha Pandey
- School of Chemistry, Indian Institute of Science Education and Research Thiruvananthapuram, Maruthamala PO, Vithura, Thiruvananthapuram, 695551, India
| | - Sivaranjana Reddy Vennapusa
- School of Chemistry, Indian Institute of Science Education and Research Thiruvananthapuram, Maruthamala PO, Vithura, Thiruvananthapuram, 695551, India.
| |
Collapse
|
40
|
Asido M, Hamerla C, Weber R, Horz M, Niraghatam MS, Heckel A, Burghardt I, Wachtveitl J. Ultrafast and efficient energy transfer in a one- and two-photon sensitized rhodamine-BODIPY dyad: a perspective for broadly absorbing photocages. Phys Chem Chem Phys 2022; 24:1795-1802. [PMID: 34985062 DOI: 10.1039/d1cp04528h] [Citation(s) in RCA: 4] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
In view of the demand for photoactivatable probes that operate in the visible (VIS) to near infrared (NIR) region of the spectrum, we designed a bichromophoric system based on a rhodamine fluorophore and a BODIPY photocage. Two-photon excited fluorescence (TPEF) measurements and quantum chemical calculations reveal excellent two-photon properties of the employed rhodamine derivative. Excitation of the rhodamine unit via a one- or two-photon process leads to excitation energy transfer (EET) onto the BODIPY part, which is followed by the liberation of the leaving group. Ultrafast transient absorption spectroscopy provides evidence for a highly efficient EET dynamics on a sub-500 femtosecond scale. Complementary quantum dynamical calculations using the multi-layer multiconfiguration time-dependent Hartree (ML-MCTDH) approach highlight the quantum coherent character of the EET transfer. Photorelease of p-nitroaniline (PNA) was investigated by UV/vis absorption spectroscopy by either excitation of the rhodamine or the BODIPY moiety. Even though a quantitative assessment of the PNA yield could not be achieved for this particular BODIPY cage, the present study provides a design principle for a class of photocages that can be broadly activated between 500 and 900 nm.
Collapse
Affiliation(s)
- Marvin Asido
- Institute of Physical and Theoretical Chemistry, Goethe University Frankfurt, Max-von-Laue Straße 7, 60438 Frankfurt am Main, Germany.
| | - Carsten Hamerla
- Institute of Physical and Theoretical Chemistry, Goethe University Frankfurt, Max-von-Laue Straße 7, 60438 Frankfurt am Main, Germany.
| | - Rebekka Weber
- Institute of Organic Chemistry and Chemical Biology, Goethe University Frankfurt, Max-von-Laue Straße 7, 60438 Frankfurt am Main, Germany.
| | - Maximiliane Horz
- Institute of Physical and Theoretical Chemistry, Goethe University Frankfurt, Max-von-Laue Straße 7, 60438 Frankfurt am Main, Germany.
| | - Madhava Shyam Niraghatam
- Institute of Physical and Theoretical Chemistry, Goethe University Frankfurt, Max-von-Laue Straße 7, 60438 Frankfurt am Main, Germany.
| | - Alexander Heckel
- Institute of Organic Chemistry and Chemical Biology, Goethe University Frankfurt, Max-von-Laue Straße 7, 60438 Frankfurt am Main, Germany.
| | - Irene Burghardt
- Institute of Physical and Theoretical Chemistry, Goethe University Frankfurt, Max-von-Laue Straße 7, 60438 Frankfurt am Main, Germany.
| | - Josef Wachtveitl
- Institute of Physical and Theoretical Chemistry, Goethe University Frankfurt, Max-von-Laue Straße 7, 60438 Frankfurt am Main, Germany.
| |
Collapse
|
41
|
Richings GW, Habershon S. Predicting Molecular Photochemistry Using Machine-Learning-Enhanced Quantum Dynamics Simulations. Acc Chem Res 2022; 55:209-220. [PMID: 34982533 DOI: 10.1021/acs.accounts.1c00665] [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/30/2022]
Abstract
The processes which occur after molecules absorb light underpin an enormous range of fundamental technologies and applications, including photocatalysis to enable new chemical transformations, sunscreens to protect against the harmful effects of UV overexposure, efficient photovoltaics for energy generation from sunlight, and fluorescent probes to image the intricate details of complex biomolecular structures. Reflecting this broad range of applications, an enormously versatile set of experiments are now regularly used to interrogate light-driven chemical dynamics, ranging from the typical ultrafast transient absorption spectroscopy used in many university laboratories to the inspiring central facilities around the world, such as the next-generation of X-ray free-electron lasers.Computer simulations of light-driven molecular and material dynamics are an essential route to analyzing the enormous amount of transient electronic and structural data produced by these experimental sources. However, to date, the direct simulation of molecular photochemistry remains a frontier challenge in computational chemical science, simultaneously demanding the accurate treatment of molecular electronic structure, nuclear dynamics, and the impact of nonadiabatic couplings.To address these important challenges and to enable new computational methods which can be integrated with state-of-the-art experimental capabilities, the past few years have seen a burst of activity in the development of "direct" quantum dynamics methods, merging the machine learning of potential energy surfaces (PESs) and nonadiabatic couplings with accurate quantum propagation schemes such as the multiconfiguration time-dependent Hartree (MCTDH) method. The result of this approach is a new generation of direct quantum dynamics tools in which PESs are generated in tandem with wave function propagation, enabling accurate "on-the-fly" simulations of molecular photochemistry. These simulations offer an alternative route toward gaining quantum dynamics insights, circumventing the challenge of generating ab initio electronic structure data for PES fitting by instead only demanding expensive energy evaluations as and when they are needed.In this Account, we describe the chronological evolution of our own contributions to this field, focusing on describing the algorithmic developments that enable direct MCTDH simulations for complex molecular systems moving on multiple coupled electronic states. Specifically, we highlight active learning strategies for generating PESs during grid-based quantum chemical dynamics simulations, and we discuss the development and impact of novel diabatization schemes to enable direct grid-based simulations of photochemical dynamics; these developments are highlighted in a series of benchmark molecular simulations of systems containing multiple nuclear degrees of freedom moving on multiple coupled electronic states. We hope that the ongoing developments reported here represent a major step forward in tools for modeling excited-state chemistry such as photodissociation, proton and electron transfer, and ultrafast energy dissipation in complex molecular systems.
Collapse
Affiliation(s)
- Gareth W. Richings
- Department of Chemistry, University of Warwick, Coventry, United Kingdom CV4 7AL
| | - Scott Habershon
- Department of Chemistry, University of Warwick, Coventry, United Kingdom CV4 7AL
| |
Collapse
|
42
|
Brey D, Binder R, Martinazzo R, Burghardt I. Signatures of coherent vibronic exciton dynamics and conformational control in two-dimensional electronic spectroscopy of conjugated polymers. Faraday Discuss 2022; 237:148-167. [DOI: 10.1039/d2fd00014h] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Two-dimensional electronic spectroscopy (2DES) signals for homo-oligomer J-aggregates are computed, with a focus on the role of structural change induced by low-frequency torsional modes along with quasi-stationary trapping effects induced...
Collapse
|
43
|
Titov E, Kopp T, Hoche J, Humeniuk A, Mitrić R. (De)localization dynamics of molecular excitons: comparison of mixed quantum–classical and fully quantum treatments. Phys Chem Chem Phys 2022; 24:12136-12148. [DOI: 10.1039/d2cp00586g] [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/21/2022]
Abstract
Molecular excitons play a central role in processes of solar energy conversion, both natural and artificial. It is therefore no wonder that numerous experimental and theoretical investigations in the last...
Collapse
|
44
|
Richings GW, Habershon S. Analyzing Grid-Based Direct Quantum Molecular Dynamics Using Non-Linear Dimensionality Reduction. Molecules 2021; 26:molecules26247418. [PMID: 34946499 PMCID: PMC8708769 DOI: 10.3390/molecules26247418] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/31/2021] [Revised: 11/30/2021] [Accepted: 12/02/2021] [Indexed: 11/18/2022] Open
Abstract
Grid-based schemes for simulating quantum dynamics, such as the multi-configuration time-dependent Hartree (MCTDH) method, provide highly accurate predictions of the coupled nuclear and electronic dynamics in molecular systems. Such approaches provide a multi-dimensional, time-dependent view of the system wavefunction represented on a coordinate grid; in the case of non-adiabatic simulations, additional information about the state populations adds a further layer of complexity. As such, wavepacket motion on potential energy surfaces which couple many nuclear and electronic degrees-of-freedom can be extremely challenging to analyse in order to extract physical insight beyond the usual expectation-value picture. Here, we show that non-linear dimensionality reduction (NLDR) methods, notably diffusion maps, can be adapted to extract information from grid-based wavefunction dynamics simulations, providing insight into key nuclear motions which explain the observed dynamics. This approach is demonstrated for 2-D and 9-D models of proton transfer in salicylaldimine, as well as 8-D and full 12-D simulations of cis-trans isomerization in ethene; these simulations demonstrate how NLDR can provide alternative views of wavefunction dynamics, and also highlight future developments.
Collapse
|
45
|
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
|
46
|
Christopoulou G, Tran T, Worth GA. Direct nonadiabatic quantum dynamics simulations of the photodissociation of phenol. Phys Chem Chem Phys 2021; 23:23684-23695. [PMID: 34642723 DOI: 10.1039/d1cp01843d] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Gaussian wavepacket methods are becoming popular for the investigation of nonadiabatic molecular dynamics. In the present work, a recently developed efficient algorithm for the Direct Dynamics variational Multi-Configurational Gaussian (DD-vMCG) method has been used to describe the multidimensional photodissociation dynamics of phenol including all degrees of freedom. Full-dimensional quantum dynamic calculations including for the first time six electronic states (1ππ, 11ππ*, 11πσ*, 21πσ*, 21ππ*, 31ππ*), along with a comparison to an existing analytical 4-state model for the potential energy surfaces are presented. Including the fifth singlet excited state is shown to have a significant effect on the nonadiabatic photodissociation of phenol to the phenoxyl radical and hydrogen atom. State population and flux analysis from the DD-vMCG simulations of phenol provided further insights into the decay mechanism, confirming the idea of rapid relaxation to the ground state through the 1ππ/11πσ* conical intersection.
Collapse
Affiliation(s)
| | - Thierry Tran
- Department of Chemistry, University College London, London WC1H 0AJ, UK.,Department of Chemistry, Imperial College London, Imperial College London, White City Campus, W12 0BZ London, UK.
| | - Graham A Worth
- Department of Chemistry, University College London, London WC1H 0AJ, UK
| |
Collapse
|
47
|
Lin K, Peng J, Gu FL, Lan Z. Simulation of Open Quantum Dynamics with Bootstrap-Based Long Short-Term Memory Recurrent Neural Network. J Phys Chem Lett 2021; 12:10225-10234. [PMID: 34647736 DOI: 10.1021/acs.jpclett.1c02672] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/13/2023]
Abstract
The recurrent neural network with the long short-term memory cell (LSTM-NN) is employed to simulate the long-time dynamics of open quantum systems. The bootstrap method is applied in the LSTM-NN construction and prediction, which provides a Monte Carlo estimation of a forecasting confidence interval. Within this approach, a large number of LSTM-NNs are constructed by resampling time-series sequences that were obtained from the early stage quantum evolution given by numerically exact multilayer multiconfigurational time-dependent Hartree method. The built LSTM-NN ensemble is used for the reliable propagation of the long-time quantum dynamics, and the simulated result is highly consistent with the exact evolution. The forecasting uncertainty that partially reflects the reliability of the LSTM-NN prediction is also given. This demonstrates the bootstrap-based LSTM-NN approach is a practical and powerful tool to propagate the long-time quantum dynamics of open systems with high accuracy and low computational cost.
Collapse
Affiliation(s)
- Kunni Lin
- Key Laboratory of Theoretical Chemistry of Environment, Ministry of Education, School of Chemistry, South China Normal University, Guangzhou 510006, P. R. China
| | - Jiawei Peng
- Guangdong Provincial Key Laboratory of Chemical Pollution and Environmental Safety and MOE Key Laboratory of Environmental Theoretical Chemistry, SCNU Environmental Research Institute, School of Environment, South China Normal University, Guangzhou 510006, P. R. China
| | - Feng Long Gu
- Key Laboratory of Theoretical Chemistry of Environment, Ministry of Education, School of Chemistry, South China Normal University, Guangzhou 510006, P. R. China
| | - Zhenggang Lan
- Guangdong Provincial Key Laboratory of Chemical Pollution and Environmental Safety and MOE Key Laboratory of Environmental Theoretical Chemistry, SCNU Environmental Research Institute, School of Environment, South China Normal University, Guangzhou 510006, P. R. China
| |
Collapse
|
48
|
Jouybari MY, Green JA, Improta R, Santoro F. The Ultrafast Quantum Dynamics of Photoexcited Adenine-Thymine Basepair Investigated with a Fragment-based Diabatization and a Linear Vibronic Coupling Model. J Phys Chem A 2021; 125:8912-8924. [PMID: 34609880 PMCID: PMC9281421 DOI: 10.1021/acs.jpca.1c08132] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
![]()
In
this contribution we present a quantum dynamical study of the
photoexcited hydrogen bonded base pair adenine–thymine (AT)
in a Watson–Crick arrangement. To that end, we parametrize
Linear Vibronic Coupling (LVC) models with Time-Dependent Density
Functional Theory (TD-DFT) calculations, exploiting a fragment diabatization
scheme (FrD) we have developed to define diabatic states on the basis
of individual chromophores in a multichromophoric system. Wavepacket
propagations were run with the multilayer extension of the Multiconfiguration
Time-Dependent Hartree method. We considered excitations to the three
lowest bright states, a ππ* state of
thymine and two ππ* states (La and Lb) of adenine, and we found that on the 100 fs time
scale the main decay pathways involve intramonomer population transfers
toward nπ* states of the same nucleobase. In AT this transfer
is less effective than in the isolated nucleobases, because hydrogen
bonding destabilizes the nπ* states. The population transfer
to the A → T charge transfer state is negligible, making the
ultrafast (femtosecond) decay through the proton coupled electron
transfer mechanism unlikely, in line with experimental results in
apolar solvents. The excitation energy transfer is also very small.
We carefully compare the predictions of LVC Hamiltonians obtained
with different sets of diabatic states, defined so to match either
local states of the two separated monomers or the base pair adiabatic
states in the Franck–Condon region. To that end we also extend
the flexibility of the FrD-LVC approach, introducing a new strategy
to define fragments diabatic states that account for the effect of
the rest of the multichromohoric system through a Molecular Mechanics
potential.
Collapse
Affiliation(s)
- Martha Yaghoubi Jouybari
- Consiglio Nazionale delle Ricerche, Istituto di Chimica dei Composti Organo Metallici (ICCOM-CNR), SS di Pisa, Area della Ricerca, via G. Moruzzi 1, I-56124 Pisa, Italy
| | - James A Green
- Consiglio Nazionale delle Ricerche, Istituto di Biostrutture e Bioimmagini (IBB-CNR), via Mezzocannone 16, I-80136 Napoli, Italy
| | - Roberto Improta
- Consiglio Nazionale delle Ricerche, Istituto di Biostrutture e Bioimmagini (IBB-CNR), via Mezzocannone 16, I-80136 Napoli, Italy
| | - Fabrizio Santoro
- Consiglio Nazionale delle Ricerche, Istituto di Chimica dei Composti Organo Metallici (ICCOM-CNR), SS di Pisa, Area della Ricerca, via G. Moruzzi 1, I-56124 Pisa, Italy
| |
Collapse
|
49
|
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
|
50
|
Fischer EW, Saalfrank P. A thermofield-based multilayer multiconfigurational time-dependent Hartree approach to non-adiabatic quantum dynamics at finite temperature. J Chem Phys 2021; 155:134109. [PMID: 34624972 DOI: 10.1063/5.0064013] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
We introduce a thermofield-based formulation of the multilayer multiconfigurational time-dependent Hartree (MCTDH) method to study finite temperature effects on non-adiabatic quantum dynamics from a non-stochastic, wave function perspective. Our approach is based on the formal equivalence of bosonic many-body theory at zero temperature with a doubled number of degrees of freedom and the thermal quasi-particle representation of bosonic thermofield dynamics (TFD). This equivalence allows for a transfer of bosonic many-body MCTDH as introduced by Wang and Thoss to the finite temperature framework of thermal quasi-particle TFD. As an application, we study temperature effects on the ultrafast internal conversion dynamics in pyrazine. We show that finite temperature effects can be efficiently accounted for in the construction of multilayer expansions of thermofield states in the framework presented herein. Furthermore, we find our results to agree well with existing studies on the pyrazine model based on the ρMCTDH method.
Collapse
Affiliation(s)
- Eric W Fischer
- Theoretische Chemie, Institut für Chemie, Universität Potsdam, Karl-Liebknecht-Straße 24-25, D-14476 Potsdam-Golm, Germany
| | - Peter Saalfrank
- Theoretische Chemie, Institut für Chemie, Universität Potsdam, Karl-Liebknecht-Straße 24-25, D-14476 Potsdam-Golm, Germany
| |
Collapse
|