1
|
Brook R, Symonds C, Shalashilin DV. Full wave function cloning for improving convergence of the multiconfigurational Ehrenfest method: Tests in the zero-temperature spin-boson model regime. J Chem Phys 2024; 161:064102. [PMID: 39120032 DOI: 10.1063/5.0221184] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/30/2024] [Accepted: 07/26/2024] [Indexed: 08/10/2024] Open
Abstract
In this paper, we report a new algorithm for creating an adaptive basis set in the Multiconfigurational Ehrenfest (MCE) method, which is termed Full Cloning (FC), and test it together with the existing Multiple Cloning (MC) using the spin-boson model at zero-temperature as a benchmark. The zero-temperature spin-boson regime is a common hurdle in the development of methods that seek to model quantum dynamics. Two versions of MCE exist. We demonstrate that MC is vital for the convergence of MCE version 2 (MCEv2). The first version (MCEv1) converges much better than MCEv2, but FC improves its convergence in a few cases where it is hard to converge it with the help of a reasonably small size of the basis set.
Collapse
Affiliation(s)
- Ryan Brook
- School of Chemistry, University of Leeds, Leeds LS2 9JT, United Kingdom
| | | | | |
Collapse
|
2
|
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
|
3
|
Makhov DV, Hutton L, Kirrander A, Shalashilin DV. Ultrafast electron diffraction of photoexcited gas-phase cyclobutanone predicted by ab initio multiple cloning simulations. J Chem Phys 2024; 160:164310. [PMID: 38661201 DOI: 10.1063/5.0203683] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/15/2024] [Accepted: 04/02/2024] [Indexed: 04/26/2024] Open
Abstract
We present the result of our calculations of ultrafast electron diffraction (UED) for cyclobutanone excited into the S2 electronic state, which is based on the non-adiabatic dynamics simulations with the Ab Initio Multiple Cloning (AIMC) method with the electronic structure calculated at the SA(3)-CASSCF(12,12)/aug-cc-pVDZ level of theory. The key features in the UED pattern were identified, which can be used to distinguish between the reaction pathways observed in the AIMC dynamics, although there is a significant overlap between representative signals due to the structural similarity of the products. The calculated UED pattern can be compared with the experiment.
Collapse
Affiliation(s)
- Dmitry V Makhov
- School of Chemistry, University of Leeds, Leeds LS2 9JT, United Kingdom
- School of Mathematics, University of Bristol, Fry Building, Woodland Road, Bristol BS8 1UG, United Kingdom
| | - Lewis Hutton
- Physical and Theoretical Chemistry Laboratory, Department of Chemistry, University of Oxford, Oxford OX1 3QZ, United Kingdom
| | - Adam Kirrander
- Physical and Theoretical Chemistry Laboratory, Department of Chemistry, University of Oxford, Oxford OX1 3QZ, United Kingdom
| | | |
Collapse
|
4
|
Klētnieks Ē, Alonso YC, Vaníček JJL. Isotope Effects on the Electronic Spectra of Ammonia from Ab Initio Semiclassical Dynamics. J Phys Chem A 2023; 127:8117-8125. [PMID: 37738161 PMCID: PMC10561269 DOI: 10.1021/acs.jpca.3c04607] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/08/2023] [Revised: 08/22/2023] [Indexed: 09/24/2023]
Abstract
Despite its simplicity, the single-trajectory thawed Gaussian approximation has proven useful for calculating the vibrationally resolved electronic spectra of molecules with weakly anharmonic potential energy surfaces. Here, we show that the thawed Gaussian approximation can capture surprisingly well even more subtle observables, such as the isotope effects in the absorption spectra, and we demonstrate it on the four isotopologues of ammonia (NH3, NDH2, ND2H, and ND3). The differences in their computed spectra are due to the differences in the semiclassical trajectories followed by the four isotopologues, and the isotope effects─narrowing of the transition band and reduction of the peak spacing─are accurately described by this semiclassical method. In contrast, the adiabatic harmonic model shows a double progression instead of the single progression seen in the experimental spectra. The vertical harmonic model correctly shows only a single progression but fails to describe the anharmonic peak spacing. Analysis of the normal-mode activation upon excitation provides insight into the elusiveness of the symmetric stretching progression in the spectra.
Collapse
Affiliation(s)
- Ēriks Klētnieks
- Laboratory of Theoretical Physical
Chemistry, Institut des Sciences et Ingénierie Chimiques, Ecole Polytechnique Fédérale de Lausanne
(EPFL), Lausanne CH-1015, Switzerland
| | - Yannick Calvino Alonso
- Laboratory of Theoretical Physical
Chemistry, Institut des Sciences et Ingénierie Chimiques, Ecole Polytechnique Fédérale de Lausanne
(EPFL), Lausanne CH-1015, Switzerland
| | - Jiří J. L. Vaníček
- Laboratory of Theoretical Physical
Chemistry, Institut des Sciences et Ingénierie Chimiques, Ecole Polytechnique Fédérale de Lausanne
(EPFL), Lausanne CH-1015, Switzerland
| |
Collapse
|
5
|
Freixas VM, Malone W, Li X, Song H, Negrin-Yuvero H, Pérez-Castillo R, White A, Gibson TR, Makhov DV, Shalashilin DV, Zhang Y, Fedik N, Kulichenko M, Messerly R, Mohanam LN, Sharifzadeh S, Bastida A, Mukamel S, Fernandez-Alberti S, Tretiak S. NEXMD v2.0 Software Package for Nonadiabatic Excited State Molecular Dynamics Simulations. J Chem Theory Comput 2023; 19:5356-5368. [PMID: 37506288 DOI: 10.1021/acs.jctc.3c00583] [Citation(s) in RCA: 3] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 07/30/2023]
Abstract
We present NEXMD version 2.0, the second release of the NEXMD (Nonadiabatic EXcited-state Molecular Dynamics) software package. Across a variety of new features, NEXMD v2.0 incorporates new implementations of two hybrid quantum-classical dynamics methods, namely, Ehrenfest dynamics (EHR) and the Ab-Initio Multiple Cloning sampling technique for Multiconfigurational Ehrenfest quantum dynamics (MCE-AIMC or simply AIMC), which are alternative options to the previously implemented trajectory surface hopping (TSH) method. To illustrate these methodologies, we outline a direct comparison of these three hybrid quantum-classical dynamics methods as implemented in the same NEXMD framework, discussing their weaknesses and strengths, using the modeled photodynamics of a polyphenylene ethylene dendrimer building block as a representative example. We also describe the expanded normal-mode analysis and constraints for both the ground and excited states, newly implemented in the NEXMD v2.0 framework, which allow for a deeper analysis of the main vibrational motions involved in vibronic dynamics. Overall, NEXMD v2.0 expands the range of applications of NEXMD to a larger variety of multichromophore organic molecules and photophysical processes involving quantum coherences and persistent couplings between electronic excited states and nuclear velocity.
Collapse
Affiliation(s)
- Victor M Freixas
- Departments of Chemistry and Physics and Astronomy, University of California, Irvine, California 92697-2025, United States
| | - Walter Malone
- Department of Physics, Tuskegee University, Tuskegee, Alabama 36088, United States
| | - Xinyang Li
- Theoretical Division, Center for Nonlinear Studies (CNLS), and Center for Integrated Nanotechnologies (CINT), Los Alamos National Laboratory, Los Alamos, New Mexico 87545, United States
| | - Huajing Song
- Theoretical Division, Center for Nonlinear Studies (CNLS), and Center for Integrated Nanotechnologies (CINT), Los Alamos National Laboratory, Los Alamos, New Mexico 87545, United States
| | - Hassiel Negrin-Yuvero
- Departamento de Ciencia y Tecnologia, Universidad Nacional de Quilmes/CONICET, B1876BXD Bernal, Argentina
| | - Royle Pérez-Castillo
- Departamento de Ciencia y Tecnologia, Universidad Nacional de Quilmes/CONICET, B1876BXD Bernal, Argentina
| | - Alexander White
- Theoretical Division, Center for Nonlinear Studies (CNLS), and Center for Integrated Nanotechnologies (CINT), Los Alamos National Laboratory, Los Alamos, New Mexico 87545, United States
| | - Tammie R Gibson
- Theoretical Division, Center for Nonlinear Studies (CNLS), and Center for Integrated Nanotechnologies (CINT), Los Alamos National Laboratory, Los Alamos, New Mexico 87545, United States
| | - Dmitry V Makhov
- School of Chemistry, University of Leeds, Leeds LS2 9JT, United Kingdom
- School of Mathematics, University of Bristol, Bristol BS8 1TW, United Kingdom
| | | | - Yu Zhang
- Theoretical Division, Center for Nonlinear Studies (CNLS), and Center for Integrated Nanotechnologies (CINT), Los Alamos National Laboratory, Los Alamos, New Mexico 87545, United States
| | - Nikita Fedik
- Theoretical Division, Center for Nonlinear Studies (CNLS), and Center for Integrated Nanotechnologies (CINT), Los Alamos National Laboratory, Los Alamos, New Mexico 87545, United States
| | - Maksim Kulichenko
- Theoretical Division, Center for Nonlinear Studies (CNLS), and Center for Integrated Nanotechnologies (CINT), Los Alamos National Laboratory, Los Alamos, New Mexico 87545, United States
| | - Richard Messerly
- Theoretical Division, Center for Nonlinear Studies (CNLS), and Center for Integrated Nanotechnologies (CINT), Los Alamos National Laboratory, Los Alamos, New Mexico 87545, United States
| | - Luke Nambi Mohanam
- Department of Electrical and Computer Engineering, College of Engineering, Boston University, Boston, Massachusetts 02215, United States
| | - Sahar Sharifzadeh
- Department of Electrical and Computer Engineering, College of Engineering, Boston University, Boston, Massachusetts 02215, United States
| | - Adolfo Bastida
- Departamento de Química Física, Universidad de Murcia, Murcia 30100, Spain
| | - Shaul Mukamel
- Departments of Chemistry and Physics and Astronomy, University of California, Irvine, California 92697-2025, United States
| | | | - Sergei Tretiak
- Theoretical Division, Center for Nonlinear Studies (CNLS), and Center for Integrated Nanotechnologies (CINT), Los Alamos National Laboratory, Los Alamos, New Mexico 87545, United States
| |
Collapse
|
6
|
Negrin-Yuvero H, Freixas VM, Ondarse-Alvarez D, Alfonso-Hernandez L, Rojas-Lorenzo G, Bastida A, Tretiak S, Fernandez-Alberti S. Vibrational Funnels for Energy Transfer in Organic Chromophores. J Phys Chem Lett 2023; 14:4673-4681. [PMID: 37167537 DOI: 10.1021/acs.jpclett.3c00748] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/13/2023]
Abstract
Photoinduced intramolecular energy transfers in multichromophoric molecules involve nonadiabatic vibronic channels that act as energy transfer funnels. They commonly take place through specific directions of motion dictated by the nonadiabatic coupling vectors. Vibrational funnels may support persistent coherences between electronic states and sometimes delineate the presence of minor alternative energy transfer pathways. The ultimate confirmation of their role on the interchromophoric energy transfer can be achieved by performing nonadiabatic excited-state molecular dynamics simulations by selectively freezing the nuclear motions in question. Our results point out this strategy as a useful tool to identify and evaluate the impact of these vibrational funnels on the energy transfer processes and guide the in silico design of materials with tunable properties and enhanced functionalities. Our work encourages applications of this methodology to different chemical and biochemical processes such as reactive scattering and protein conformational changes, to name a few.
Collapse
Affiliation(s)
- Hassiel Negrin-Yuvero
- Departamento de Ciencia y Tecnologia, Universidad Nacional de Quilmes/CONICET, Bernal B1876BXD, Argentina
| | - Victor Manuel Freixas
- Departamento de Ciencia y Tecnologia, Universidad Nacional de Quilmes/CONICET, Bernal B1876BXD, Argentina
| | - Dianelys Ondarse-Alvarez
- Departamento de Ciencia y Tecnologia, Universidad Nacional de Quilmes/CONICET, Bernal B1876BXD, Argentina
| | - Laura Alfonso-Hernandez
- Departamento de Ciencia y Tecnologia, Universidad Nacional de Quilmes/CONICET, Bernal B1876BXD, Argentina
| | - German Rojas-Lorenzo
- Departamento de Física Atómica y Molecular, Instituto Superior de Tecnologías y Ciencias Aplicadas, Universidad de La Habana, La Habana 10400, Cuba
| | - Adolfo Bastida
- Departamento de Química Física, Universidad de Murcia, Murcia 30100, Spain
| | - Sergei Tretiak
- Theoretical Division, Center for Nonlinear Studies (CNLS), and Center for Integrated Nanotechnologies (CINT), Los Alamos National Laboratory, Los Alamos 87545, New Mexico, USA
| | | |
Collapse
|
7
|
Nam Y, Song H, Freixas VM, Keefer D, Fernandez-Alberti S, Lee JY, Garavelli M, Tretiak S, Mukamel S. Monitoring vibronic coherences and molecular aromaticity in photoexcited cyclooctatetraene with an X-ray probe: a simulation study. Chem Sci 2023; 14:2971-2982. [PMID: 36937575 PMCID: PMC10016608 DOI: 10.1039/d2sc04335a] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/03/2022] [Accepted: 02/07/2023] [Indexed: 02/16/2023] Open
Abstract
Understanding conical intersection (CI) dynamics and subsequent conformational changes is key for exploring and controlling photo-reactions in aromatic molecules. Monitoring of their time-resolved dynamics remains a formidable experimental challenge. In this study, we simulate the photoinduced S3 to S1 non-adiabatic dynamics of cyclooctatetraene (COT), involving multiple CIs with relaxation times in good agreement with experiment. We further investigate the possibility to directly probe the CI passages in COT by off-resonant X-ray Raman spectroscopy (TRUECARS) and time-resolved X-ray diffraction (TRXD). We find that these signals sensitively monitor key chemical features during the ultrafast dynamics. First, we distinguish two CIs by using TRUECARS signals with their appearances at different Raman shifts. Second, we demonstrate that TRXD, where X-ray photons scatter off electron densities, can resolve ultrafast changes in the aromaticity of COT. It can further distinguish between planar and non-planar geometries explored during the dynamics, as e.g. two different tetraradical-type CIs. The knowledge gained from these measurements can give unique insight into fundamental chemical properties that dynamically change during non-adiabatic passages.
Collapse
Affiliation(s)
- Yeonsig Nam
- Department of Chemistry, University of California Irvine California 92697-2025 USA
| | - Huajing Song
- Physics and Chemistry of Materials, Theoretical Division, Los Alamos National Laboratory Los Alamos New Mexico 87545 USA
| | - Victor M Freixas
- Departamento de Ciencia y Tecnologia, Universidad Nacional de Quilmes/CONICET B1876BXD Bernal Argentina
| | - Daniel Keefer
- Department of Chemistry, University of California Irvine California 92697-2025 USA
| | | | - Jin Yong Lee
- Department of Chemistry, Sungkyunkwan University Suwon 16419 Korea
| | - Marco Garavelli
- Dipartimento di Chimica Industriale "Toso Montanari,", Universita' degli Studi di Bologna I-40136 Bologna Italy
| | - Sergei Tretiak
- Physics and Chemistry of Materials, Theoretical Division, Los Alamos National Laboratory Los Alamos New Mexico 87545 USA
| | - Shaul Mukamel
- Department of Chemistry, University of California Irvine California 92697-2025 USA
| |
Collapse
|
8
|
Freixas VM, Keefer D, Tretiak S, Fernandez-Alberti S, Mukamel S. Ultrafast coherent photoexcited dynamics in a trimeric dendrimer probed by X-ray stimulated-Raman signals. Chem Sci 2022; 13:6373-6384. [PMID: 35733898 PMCID: PMC9159119 DOI: 10.1039/d2sc00601d] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/28/2022] [Accepted: 04/11/2022] [Indexed: 12/14/2022] Open
Abstract
The photoinduced ultrafast coherent inter-chromophore energy redistribution in a triarylamine trimer is explored using nonadiabatic excited state molecular dynamics followed by simulations of X-ray Raman signals. The nitrogencentered system ensures strong interchromophore interactions and, thus, the presence of coherences. Nevertheless, the multitude of non-deterministic photoinduced pathways during the ultrafast inter-branch migration of the excitation results in random confinement on some branches and, therefore, spatial exciton scrambling and loss of phase information at long times. We show that the vibronic coherence dynamics evolving into the incoherent scrambling mechanism on ultrafast 50 fs timescale, is accurately probed by the TRUECARS X-ray stimulated Raman signal. In combination with previous results, where the technique has revealed long-lived coherences in a rigid heterodimer, the signal is most valuable for detecting ultrafast molecular coherences or their absence. We demonstrate that X-ray Raman spectroscopy is a useful tool in the chemical design of functional molecular building blocks. The photoinduced ultrafast coherent inter-chromophore energy redistribution in a triarylamine trimer is explored using nonadiabatic excited state molecular dynamics followed by simulations of X-ray Raman signals.![]()
Collapse
Affiliation(s)
- Victor M Freixas
- Departamento de Ciencia y Tecnologia, Universidad Nacional de Quilmes/CONICET B1876BXD Bernal Argentina
| | - Daniel Keefer
- Department of Chemistry and Physics and Astronomy, University of California Irvine California 92697-2025 USA
| | - Sergei Tretiak
- Theoretical Division, Los Alamos National Laboratory Los Alamos New Mexico 87545 USA
| | | | - Shaul Mukamel
- Department of Chemistry and Physics and Astronomy, University of California Irvine California 92697-2025 USA
| |
Collapse
|
9
|
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
|
10
|
Esch MP, Levine BG. An accurate, non-empirical method for incorporating decoherence into Ehrenfest dynamics. J Chem Phys 2021; 155:214101. [PMID: 34879667 DOI: 10.1063/5.0070686] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
In mixed quantum-classical nonadiabatic molecular dynamics methods, the anchoring of the electronic wave function to a single nuclear geometry results in both quantitative and qualitative errors in the dynamics. In the context of both Ehrenfest and trajectory surface hopping methods, methods for incorporating decoherence are widely used to eliminate these errors. However, the accuracy of these methods often depends strongly on the parameterization of the decoherence time and/or other related quantities. Here, we present a refinement of the recently introduced collapse to a block (TAB) scheme for incorporating decoherence into Ehrenfest dynamics. The proposed approach incorporates an approximation to the history of the population dynamics and treats the coherence decay as Gaussian, rather than exponential. This method uses parameters that can be obtained from first principles, rather than empirical fitting. Application to one-dimensional models indicates excellent agreement with numerically exact simulations. We also introduce a second refinement to the TAB method: a robust linear least-squares algorithm for determining collapse probabilities.
Collapse
Affiliation(s)
- Michael P Esch
- Department of Chemistry and Institute for Advanced Computational Science, Stony Brook University, Stony Brook, New York 11794, USA
| | - Benjamin G Levine
- Department of Chemistry and Institute for Advanced Computational Science, Stony Brook University, Stony Brook, New York 11794, USA
| |
Collapse
|
11
|
Westermayr J, Marquetand P. Machine Learning for Electronically Excited States of Molecules. Chem Rev 2021; 121:9873-9926. [PMID: 33211478 PMCID: PMC8391943 DOI: 10.1021/acs.chemrev.0c00749] [Citation(s) in RCA: 162] [Impact Index Per Article: 54.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/17/2020] [Indexed: 12/11/2022]
Abstract
Electronically excited states of molecules are at the heart of photochemistry, photophysics, as well as photobiology and also play a role in material science. Their theoretical description requires highly accurate quantum chemical calculations, which are computationally expensive. In this review, we focus on not only how machine learning is employed to speed up such excited-state simulations but also how this branch of artificial intelligence can be used to advance this exciting research field in all its aspects. Discussed applications of machine learning for excited states include excited-state dynamics simulations, static calculations of absorption spectra, as well as many others. In order to put these studies into context, we discuss the promises and pitfalls of the involved machine learning techniques. Since the latter are mostly based on quantum chemistry calculations, we also provide a short introduction into excited-state electronic structure methods and approaches for nonadiabatic dynamics simulations and describe tricks and problems when using them in machine learning for excited states of molecules.
Collapse
Affiliation(s)
- Julia Westermayr
- Institute
of Theoretical Chemistry, Faculty of Chemistry, University of Vienna, Währinger Strasse 17, 1090 Vienna, Austria
| | - Philipp Marquetand
- Institute
of Theoretical Chemistry, Faculty of Chemistry, University of Vienna, Währinger Strasse 17, 1090 Vienna, Austria
- Vienna
Research Platform on Accelerating Photoreaction Discovery, University of Vienna, Währinger Strasse 17, 1090 Vienna, Austria
- Data
Science @ Uni Vienna, University of Vienna, Währinger Strasse 29, 1090 Vienna, Austria
| |
Collapse
|
12
|
Abstract
Electronically excited states of molecules are at the heart of photochemistry, photophysics, as well as photobiology and also play a role in material science. Their theoretical description requires highly accurate quantum chemical calculations, which are computationally expensive. In this review, we focus on not only how machine learning is employed to speed up such excited-state simulations but also how this branch of artificial intelligence can be used to advance this exciting research field in all its aspects. Discussed applications of machine learning for excited states include excited-state dynamics simulations, static calculations of absorption spectra, as well as many others. In order to put these studies into context, we discuss the promises and pitfalls of the involved machine learning techniques. Since the latter are mostly based on quantum chemistry calculations, we also provide a short introduction into excited-state electronic structure methods and approaches for nonadiabatic dynamics simulations and describe tricks and problems when using them in machine learning for excited states of molecules.
Collapse
Affiliation(s)
- Julia Westermayr
- Institute of Theoretical Chemistry, Faculty of Chemistry, University of Vienna, Währinger Strasse 17, 1090 Vienna, Austria
| | - Philipp Marquetand
- Institute of Theoretical Chemistry, Faculty of Chemistry, University of Vienna, Währinger Strasse 17, 1090 Vienna, Austria
- Vienna Research Platform on Accelerating Photoreaction Discovery, University of Vienna, Währinger Strasse 17, 1090 Vienna, Austria
- Data Science @ Uni Vienna, University of Vienna, Währinger Strasse 29, 1090 Vienna, Austria
| |
Collapse
|
13
|
Song H, Freixas VM, Fernandez-Alberti S, White AJ, Zhang Y, Mukamel S, Govind N, Tretiak S. An Ab Initio Multiple Cloning Method for Non-Adiabatic Excited-State Molecular Dynamics in NWChem. J Chem Theory Comput 2021; 17:3629-3643. [PMID: 34014085 DOI: 10.1021/acs.jctc.1c00131] [Citation(s) in RCA: 7] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
The recently developed ab initio multiple cloning (AIMC) approach based on the multiconfigurational Ehrenfest (MCE) method provides a powerful and accurate way of describing the excited-state dynamics of molecular systems. The AIMC method is a controlled approximation to nonadiabatic dynamics with a particular strength in the proper description of decoherence effects because of the branching of vibrational wavepackets at a level crossing. Here, we report a new implementation of the AIMC algorithm in the open source NWChem computational chemistry program. The framework combines linear-response time-dependent density functional theory with Ehrenfest mean-field theory to determine the equations of motion for classical trajectories. The multidimensional wave function is decomposed into a superposition of Gaussian coherent states guided by Ehrenfest trajectories (i.e., MCE approach), which can clone with fully quantum mechanical amplitudes and phases. By using an efficient time-derivative based nonadiabatic coupling approach within the AIMC method, all observables are calculated on-the-fly in the nonadiabatic molecular dynamics process. As a representative example, we apply our implementation to study the ultrafast photoinduced electronic and vibrational energy transfer in a pyridine molecule. The effects of the cloning procedure on electronic and vibrational coherence, relaxation and unidirectional energy transfer are discussed. This new AIMC implementation provides a high-level nonadiabatic molecular dynamics framework for simulating photoexcited dynamics in complex molecular systems and experimentally relevant ultrafast spectroscopic probes, such as nonlinear coherent optical and X-ray signals.
Collapse
Affiliation(s)
- Huajing Song
- Physics and Chemistry of Materials, Theoretical Division, Los Alamos National Laboratory, Los Alamos, New Mexico 87545, United States
| | - Victor M Freixas
- Departamento de Ciencia y Tecnologia, Universidad Nacional de Quilmes/CONICET, B1876BXD, Bernal, Argentina
| | | | - Alexander J White
- Physics and Chemistry of Materials, Theoretical Division, Los Alamos National Laboratory, Los Alamos, New Mexico 87545, United States
| | - Yu Zhang
- Physics and Chemistry of Materials, Theoretical Division, Los Alamos National Laboratory, Los Alamos, New Mexico 87545, United States
| | - Shaul Mukamel
- Departments of Chemistry, Physics, and Astronomy, University of California, Irvine, California 92697, United States
| | - Niranjan Govind
- Physical and Computational Sciences Directorate, Pacific Northwest National Laboratory, Richland, Washington 99352, United States
| | - Sergei Tretiak
- Physics and Chemistry of Materials, Theoretical Division, Los Alamos National Laboratory, Los Alamos, New Mexico 87545, United States.,Center for Integrated Nanotechnologies, Los Alamos National Laboratory, Los Alamos, New Mexico 87545, United States
| |
Collapse
|
14
|
Freixas VM, White AJ, Nelson T, Song H, Makhov DV, Shalashilin D, Fernandez-Alberti S, Tretiak S. Nonadiabatic Excited-State Molecular Dynamics Methodologies: Comparison and Convergence. J Phys Chem Lett 2021; 12:2970-2982. [PMID: 33730495 DOI: 10.1021/acs.jpclett.1c00266] [Citation(s) in RCA: 17] [Impact Index Per Article: 5.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/12/2023]
Abstract
Direct atomistic simulation of nonadiabatic molecular dynamics is a challenging goal that allows important insights into fundamental physical phenomena. A variety of frameworks, ranging from fully quantum treatment of nuclei to semiclassical and mixed quantum-classical approaches, were developed. These algorithms are then coupled to specific electronic structure techniques. Such diversity and lack of standardized implementation make it difficult to compare the performance of different methodologies when treating realistic systems. Here, we compare three popular methods for large chromophores: Ehrenfest, surface hopping, and multiconfigurational Ehrenfest with ab initio multiple cloning (MCE-AIMC). These approaches are implemented in the NEXMD software, which features a common computational chemistry model. The resulting comparisons reveal the method performance for population relaxation and coherent vibronic dynamics. Finally, we study the numerical convergence of MCE-AIMC algorithms by considering the number of trajectories, cloning thresholds, and Gaussian wavepacket width. Our results provide helpful reference data for selecting an optimal methodology for simulating excited-state molecular dynamics.
Collapse
Affiliation(s)
- Victor M Freixas
- Universidad Nacional de Quilmes, Roque Saénz Peña 352, B1876BXD Bernal, Argentina
| | - Alexander J White
- Theoretical Division, Los Alamos National Laboratory, Los Alamos, New Mexico 87545, United States
| | - Tammie Nelson
- Theoretical Division, Los Alamos National Laboratory, Los Alamos, New Mexico 87545, United States
| | - Huajing Song
- Theoretical Division, Los Alamos National Laboratory, Los Alamos, New Mexico 87545, United States
| | - Dmitry V Makhov
- School of Chemistry, University of Leeds, Leeds LS2 9JT, U.K
- School of Mathematics, University of Bristol, Bristol BS8 1TW, U.K
| | | | | | - Sergei Tretiak
- Theoretical Division, Los Alamos National Laboratory, Los Alamos, New Mexico 87545, United States
| |
Collapse
|
15
|
Makhov DV, Shalashilin DV. Simulation of the effect of vibrational pre-excitation on the dynamics of pyrrole photo-dissociation. J Chem Phys 2021; 154:104119. [DOI: 10.1063/5.0040178] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/31/2023] Open
Affiliation(s)
- Dmitry V. Makhov
- School of Chemistry, University of Leeds, Leeds LS2 9JT, United Kingdom
- School of Mathematics, University of Bristol, Bristol BS8 1TW, United Kingdom
| | | |
Collapse
|
16
|
Keefer D, Freixas VM, Song H, Tretiak S, Fernandez-Alberti S, Mukamel S. Monitoring molecular vibronic coherences in a bichromophoric molecule by ultrafast X-ray spectroscopy. Chem Sci 2021; 12:5286-5294. [PMID: 34168779 PMCID: PMC8179640 DOI: 10.1039/d0sc06328b] [Citation(s) in RCA: 14] [Impact Index Per Article: 4.7] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/17/2020] [Accepted: 02/24/2021] [Indexed: 11/21/2022] Open
Abstract
The role of quantum-mechanical coherences in the elementary photophysics of functional optoelectronic molecular materials is currently under active study. Designing and controlling stable coherences arising from concerted vibronic dynamics in organic chromophores is the key for numerous applications. Here, we present fundamental insight into the energy transfer properties of a rigid synthetic heterodimer that has been experimentally engineered to study coherences. Quantum non-adiabatic excited state simulations are used to compute X-ray Raman signals, which are able to sensitively monitor the coherence evolution. Our results verify their vibronic nature, that survives multiple conical intersection passages for several hundred femtoseconds at room temperature. Despite the contributions of highly heterogeneous evolution pathways, the coherences are unambiguously visualized by the experimentally accessible X-ray signals. They offer direct information on the dynamics of electronic and structural degrees of freedom, paving the way for detailed coherence measurements in functional organic materials.
Collapse
Affiliation(s)
- Daniel Keefer
- Departments of Chemistry and Physics and Astronomy, University of California Irvine California 92697-2025 USA
| | - Victor M Freixas
- Departamento de Ciencia y Tecnologia, Universidad Nacional de Quilmes/CONICET B1876BXD Bernal Argentina
| | - Huajing Song
- Theoretical Division, Los Alamos National Laboratory Los Alamos New Mexico 87545 USA
| | - Sergei Tretiak
- Theoretical Division, Los Alamos National Laboratory Los Alamos New Mexico 87545 USA
| | | | - Shaul Mukamel
- Departments of Chemistry and Physics and Astronomy, University of California Irvine California 92697-2025 USA
| |
Collapse
|
17
|
Freixas VM, Nelson T, Ondarse-Alvarez D, Nijjar P, Mikhailovsky A, Zhou C, Fernandez-Alberti S, Bazan GC, Tretiak S. Experimental and theoretical study of energy transfer in a chromophore triad: What makes modeling dynamics successful? J Chem Phys 2020; 153:244114. [PMID: 33380074 DOI: 10.1063/5.0028126] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/31/2022] Open
Abstract
Simulation of electronic dynamics in realistically large molecular systems is a demanding task that has not yet achieved the same level of quantitative prediction already realized for its static counterpart. This is particularly true for processes occurring beyond the Born-Oppenheimer regime. Non-adiabatic molecular dynamics (NAMD) simulations suffer from two convoluted sources of error: numerical algorithms for dynamics and electronic structure calculations. While the former has gained increasing attention, particularly addressing the validity of ad hoc methodologies, the effect of the latter remains relatively unexplored. Indeed, the required accuracy for electronic structure calculations to reach quantitative agreement with experiment in dynamics may be even more strict than that required for static simulations. Here, we address this issue by modeling the electronic energy transfer in a donor-acceptor-donor (D-A-D) molecular light harvesting system using fewest switches surface hopping NAMD simulations. In the studied system, time-resolved experimental measurements deliver complete information on spectra and energy transfer rates. Subsequent modeling shows that the calculated electronic transition energies are "sufficiently good" to reproduce experimental spectra but produce over an order of magnitude error in simulated dynamical rates. We further perform simulations using artificially shifted energy gaps to investigate the complex relationship between transition energies and modeled dynamics to understand factors affecting non-radiative relaxation and energy transfer rates.
Collapse
Affiliation(s)
- Victor M Freixas
- Universidad Nacional de Quilmes/CONICET, Roque Saenz Peña 352, B1876BXD Bernal, Argentina
| | - Tammie Nelson
- Physics and Chemistry of Materials, Theoretical Division, Los Alamos National Laboratory, Los Alamos, New Mexico 87545, USA
| | | | - Parmeet Nijjar
- Department of Chemistry, University of Southern California, Los Angeles, California 90089-1062, USA
| | - Alexander Mikhailovsky
- Department of Chemistry and Biochemistry, Center for Polymers and Organic Solids, University of California Santa Barbara, Santa Barbara, California 93106, USA
| | - Cheng Zhou
- Department of Chemistry, National University of Singapore, Singapore 117543, Singapore
| | | | - Guillermo C Bazan
- Department of Chemistry and Biochemistry, Center for Polymers and Organic Solids, University of California Santa Barbara, Santa Barbara, California 93106, USA
| | - Sergei Tretiak
- Physics and Chemistry of Materials, Theoretical Division, Los Alamos National Laboratory, Los Alamos, New Mexico 87545, USA
| |
Collapse
|
18
|
Esch MP, Levine BG. Decoherence-corrected Ehrenfest molecular dynamics on many electronic states. J Chem Phys 2020; 153:114104. [DOI: 10.1063/5.0022529] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/06/2023] Open
Affiliation(s)
- Michael P. Esch
- Department of Chemistry, Michigan State University, East Lansing, Michigan 48824, USA
| | - Benjamin G. Levine
- Department of Chemistry, Michigan State University, East Lansing, Michigan 48824, USA
| |
Collapse
|
19
|
Song H, Fischer SA, Zhang Y, Cramer CJ, Mukamel S, Govind N, Tretiak S. First Principles Nonadiabatic Excited-State Molecular Dynamics in NWChem. J Chem Theory Comput 2020; 16:6418-6427. [DOI: 10.1021/acs.jctc.0c00295] [Citation(s) in RCA: 9] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/13/2022]
Affiliation(s)
- Huajing Song
- Physics and Chemistry of Materials, Theoretical Division, Los Alamos National Laboratory, Los Alamos, New Mexico, 87545, United States
| | - Sean A. Fischer
- Chemistry Division, U.S. Naval Research Laboratory, Washington, District of Columbia 20375, United States
| | - Yu Zhang
- Physics and Chemistry of Materials, Theoretical Division, Los Alamos National Laboratory, Los Alamos, New Mexico, 87545, United States
| | - Christopher J. Cramer
- Department of Chemistry, Supercomputing Institute and Chemical Theory Center, University of Minnesota, Minneapolis, Minnesota 55455, United States
| | - Shaul Mukamel
- Departments of Chemistry, and physics and astronomy, University of California, Irvine, California 92697, United States
| | - Niranjan Govind
- Physical and Computational Sciences Directorate, Pacific Northwest National Laboratory, Richland, Washington 99352, United States
| | - Sergei Tretiak
- Physics and Chemistry of Materials, Theoretical Division, Los Alamos National Laboratory, Los Alamos, New Mexico, 87545, United States
- Center for Integrated Nanotechnologies, Los Alamos National Laboratory, Los Alamos, New Mexico 87545, United States
| |
Collapse
|
20
|
Fedorov DA, Seritan S, Fales BS, Martínez TJ, Levine BG. PySpawn: Software for Nonadiabatic Quantum Molecular Dynamics. J Chem Theory Comput 2020; 16:5485-5498. [PMID: 32687710 DOI: 10.1021/acs.jctc.0c00575] [Citation(s) in RCA: 9] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/13/2023]
Abstract
The ab initio multiple spawning (AIMS) method enables nonadiabatic quantum molecular dynamics simulations in an arbitrary number of dimensions, with potential energy surfaces provided by electronic structure calculations performed on-the-fly. However, the intricacy of the AIMS algorithm complicates software development, deployment on modern shared computer resources, and postsimulation data analysis. PySpawn is a nonadiabatic molecular dynamics software package that addresses these issues. The program is designed to be easily interfaced with electronic structure software, and an interface to the TeraChem software package is described here. PySpawn introduces a task-based reorganization of the AIMS algorithm, allowing fine-grained restart capability and setting the stage for efficient parallelization in a future release. PySpawn includes a user-friendly and interactive Python analysis module that will enable novice users to painlessly adopt AIMS. As a demonstration of PySpawn's simulation capability and analysis module, we report complete active space self-consistent field-based AIMS simulations of the 1,2-dithienyl-1,2-dicyanoethene molecule, a promising molecular photoswitch.
Collapse
Affiliation(s)
- Dmitry A Fedorov
- Department of Chemistry, Michigan State University, East Lansing, Michigan 48824, United States
| | - Stefan Seritan
- Department of Chemistry and the PULSE Institute, Stanford University, Stanford, California 94305, United States.,SLAC National Accelerator Laboratory, Menlo Park, California 94305, United States
| | - B Scott Fales
- Department of Chemistry and the PULSE Institute, Stanford University, Stanford, California 94305, United States.,SLAC National Accelerator Laboratory, Menlo Park, California 94305, United States
| | - Todd J Martínez
- Department of Chemistry and the PULSE Institute, Stanford University, Stanford, California 94305, United States.,SLAC National Accelerator Laboratory, Menlo Park, California 94305, United States
| | - Benjamin G Levine
- Department of Chemistry, Michigan State University, East Lansing, Michigan 48824, United States
| |
Collapse
|
21
|
Malone W, Nebgen B, White A, Zhang Y, Song H, Bjorgaard JA, Sifain AE, Rodriguez-Hernandez B, Freixas VM, Fernandez-Alberti S, Roitberg AE, Nelson TR, Tretiak S. NEXMD Software Package for Nonadiabatic Excited State Molecular Dynamics Simulations. J Chem Theory Comput 2020; 16:5771-5783. [DOI: 10.1021/acs.jctc.0c00248] [Citation(s) in RCA: 33] [Impact Index Per Article: 8.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/20/2022]
Affiliation(s)
- Walter Malone
- Physics and Chemistry of Materials, Theoretical Division, Los Alamos National Laboratory, Los Alamos, New Mexico 87545, United States
- Center for Nonlinear Studies, Los Alamos National Laboratory, Los Alamos, New Mexico 87545, United States
| | - Benjamin Nebgen
- Physics and Chemistry of Materials, Theoretical Division, Los Alamos National Laboratory, Los Alamos, New Mexico 87545, United States
| | - Alexander White
- Physics and Chemistry of Materials, Theoretical Division, Los Alamos National Laboratory, Los Alamos, New Mexico 87545, United States
| | - Yu Zhang
- Physics and Chemistry of Materials, Theoretical Division, Los Alamos National Laboratory, Los Alamos, New Mexico 87545, United States
| | - Huajing Song
- Physics and Chemistry of Materials, Theoretical Division, Los Alamos National Laboratory, Los Alamos, New Mexico 87545, United States
| | - Josiah A. Bjorgaard
- Physics and Chemistry of Materials, Theoretical Division, Los Alamos National Laboratory, Los Alamos, New Mexico 87545, United States
| | - Andrew E. Sifain
- U.S. Army Research Laboratory, Aberdeen Proving Ground, Aberdeen, Maryland 21005, United States
| | | | - Victor M. Freixas
- Universidad Nacional de Quilmes/CONICET, Roque Saenz Peña 352, B1876BXD Bernal, Argentina
| | | | - Adrian E. Roitberg
- Department of Chemistry, University of Florida, Gainesville, Florida 32611, United States
| | - Tammie R. Nelson
- Physics and Chemistry of Materials, Theoretical Division, Los Alamos National Laboratory, Los Alamos, New Mexico 87545, United States
- Center for Integrated Nanotechnologies, Los Alamos National Laboratory, Los Alamos, New Mexico 87545, United States
| | - Sergei Tretiak
- Physics and Chemistry of Materials, Theoretical Division, Los Alamos National Laboratory, Los Alamos, New Mexico 87545, United States
| |
Collapse
|
22
|
Freixas VM, Tretiak S, Makhov DV, Shalashilin DV, Fernandez-Alberti S. Vibronic Quantum Beating between Electronic Excited States in a Heterodimer. J Phys Chem B 2020; 124:3992-4001. [DOI: 10.1021/acs.jpcb.0c01685] [Citation(s) in RCA: 9] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- V. M. Freixas
- Departamento de Ciencia y Tecnologia, Universidad Nacional de Quilmes/CONICET, B1876BXD Bernal, Argentina
| | - S. Tretiak
- Theoretical Division, Center for Nonlinear Studies (CNLS), and Center for Integrated Nanotechnologies (CINT), Los Alamos National Laboratory, Los Alamos, New Mexico 87545, United States
| | - D. V. Makhov
- School of Chemistry, University of Leeds, Leeds LS2 9JT, U.K
- School of Mathematics, University of Bristol, Bristol BS8 1TW, U.K
| | | | - S. Fernandez-Alberti
- Departamento de Ciencia y Tecnologia, Universidad Nacional de Quilmes/CONICET, B1876BXD Bernal, Argentina
| |
Collapse
|
23
|
Zhang Y, Li L, Tretiak S, Nelson T. Nonadiabatic Excited-State Molecular Dynamics for Open-Shell Systems. J Chem Theory Comput 2020; 16:2053-2064. [DOI: 10.1021/acs.jctc.9b00928] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/07/2023]
Affiliation(s)
- Yu Zhang
- Physics and Chemistry of Materials, Theoretical Division, Los Alamos National Laboratory, Los Alamos, New Mexico 87545, United States
| | - Linqiu Li
- Department of Chemistry, University of Southern California, Los Angeles, California 90089, United States
| | - Sergei Tretiak
- Physics and Chemistry of Materials, Theoretical Division, Los Alamos National Laboratory, Los Alamos, New Mexico 87545, United States
- Center for Integrated Nanotechnologies, Los Alamos National Laboratory, Los Alamos, New Mexico 87545, United States
| | - Tammie Nelson
- Physics and Chemistry of Materials, Theoretical Division, Los Alamos National Laboratory, Los Alamos, New Mexico 87545, United States
| |
Collapse
|
24
|
Nelson TR, White AJ, Bjorgaard JA, Sifain AE, Zhang Y, Nebgen B, Fernandez-Alberti S, Mozyrsky D, Roitberg AE, Tretiak S. Non-adiabatic Excited-State Molecular Dynamics: Theory and Applications for Modeling Photophysics in Extended Molecular Materials. Chem Rev 2020; 120:2215-2287. [PMID: 32040312 DOI: 10.1021/acs.chemrev.9b00447] [Citation(s) in RCA: 198] [Impact Index Per Article: 49.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/22/2022]
Abstract
Optically active molecular materials, such as organic conjugated polymers and biological systems, are characterized by strong coupling between electronic and vibrational degrees of freedom. Typically, simulations must go beyond the Born-Oppenheimer approximation to account for non-adiabatic coupling between excited states. Indeed, non-adiabatic dynamics is commonly associated with exciton dynamics and photophysics involving charge and energy transfer, as well as exciton dissociation and charge recombination. Understanding the photoinduced dynamics in such materials is vital to providing an accurate description of exciton formation, evolution, and decay. This interdisciplinary field has matured significantly over the past decades. Formulation of new theoretical frameworks, development of more efficient and accurate computational algorithms, and evolution of high-performance computer hardware has extended these simulations to very large molecular systems with hundreds of atoms, including numerous studies of organic semiconductors and biomolecules. In this Review, we will describe recent theoretical advances including treatment of electronic decoherence in surface-hopping methods, the role of solvent effects, trivial unavoided crossings, analysis of data based on transition densities, and efficient computational implementations of these numerical methods. We also emphasize newly developed semiclassical approaches, based on the Gaussian approximation, which retain phase and width information to account for significant decoherence and interference effects while maintaining the high efficiency of surface-hopping approaches. The above developments have been employed to successfully describe photophysics in a variety of molecular materials.
Collapse
Affiliation(s)
- Tammie R Nelson
- Theoretical Division , Los Alamos National Laboratory , Los Alamos , New Mexico 87545 , United States
| | - Alexander J White
- Theoretical Division , Los Alamos National Laboratory , Los Alamos , New Mexico 87545 , United States
| | - Josiah A Bjorgaard
- Theoretical Division , Los Alamos National Laboratory , Los Alamos , New Mexico 87545 , United States
| | - Andrew E Sifain
- Theoretical Division , Los Alamos National Laboratory , Los Alamos , New Mexico 87545 , United States.,U.S. Army Research Laboratory , Aberdeen Proving Ground , Maryland 21005 , United States
| | - Yu Zhang
- Theoretical Division , Los Alamos National Laboratory , Los Alamos , New Mexico 87545 , United States
| | - Benjamin Nebgen
- Theoretical Division , Los Alamos National Laboratory , Los Alamos , New Mexico 87545 , United States
| | | | - Dmitry Mozyrsky
- Theoretical Division , Los Alamos National Laboratory , Los Alamos , New Mexico 87545 , United States
| | - Adrian E Roitberg
- Department of Chemistry , University of Florida , Gainesville , Florida 32611 , United States
| | - Sergei Tretiak
- Theoretical Division , Los Alamos National Laboratory , Los Alamos , New Mexico 87545 , United States
| |
Collapse
|
25
|
Titov E, Humeniuk A, Mitrić R. Comparison of moving and fixed basis sets for nonadiabatic quantum dynamics at conical intersections. Chem Phys 2020. [DOI: 10.1016/j.chemphys.2019.110526] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/24/2023]
|
26
|
Fedorov DA, Levine BG. Nonadiabatic Quantum Molecular Dynamics in Dense Manifolds of Electronic States. J Phys Chem Lett 2019; 10:4542-4548. [PMID: 31342748 DOI: 10.1021/acs.jpclett.9b01902] [Citation(s) in RCA: 13] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/10/2023]
Abstract
Most nonadiabatic molecular dynamics methods require the determination of a basis of adiabatic or diabatic electronic states at every time step, but in dense manifolds of electronic states, such approaches become intractable. A notable exception is Ehrenfest molecular dynamics, which can be implemented without explicit determination of such a basis but suffers from unphysical behavior when propagation on a mean-field potential energy surface (PES) does not accurately reflect the true dynamics on multiple electronic states. Here we introduce the multiple cloning for dense manifolds of states (MCDMS) method, a systematically improvable approximation to the multiple cloning method. MCDMS avoids both the mean-field PES problem and the need to compute the full electronic spectrum. This is achieved by reformulating multiple cloning to use a subspace of approximate eigenstates constructed from the time-dependent Ehrenfest electronic wave function. By application to model systems, we show that this approach allows a substantial reduction in the size of the required electronic basis without significant loss in accuracy.
Collapse
Affiliation(s)
- Dmitry A Fedorov
- Department of Chemistry , Michigan State University , East Lansing , Michigan 48824 , United States
| | - Benjamin G Levine
- Department of Chemistry , Michigan State University , East Lansing , Michigan 48824 , United States
| |
Collapse
|
27
|
Freixas VM, Ondarse-Alvarez D, Tretiak S, Makhov DV, Shalashilin DV, Fernandez-Alberti S. Photoinduced non-adiabatic energy transfer pathways in dendrimer building blocks. J Chem Phys 2019; 150:124301. [DOI: 10.1063/1.5086680] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- V. M. Freixas
- Departamento de Ciencia y Tecnologia, Universidad Nacional de Quilmes/CONICET, B1876BXD Bernal, Argentina
| | - D. Ondarse-Alvarez
- Departamento de Ciencia y Tecnologia, Universidad Nacional de Quilmes/CONICET, B1876BXD Bernal, Argentina
| | - S. Tretiak
- Theoretical Division, Center for Nonlinear Studies (CNLS), and Center for Integrated Nanotechnologies (CINT), Los Alamos National Laboratory, Los Alamos, New Mexico 87545, USA
| | - D. V. Makhov
- School of Chemistry, University of Leeds, Leeds LS2 9JT, United Kingdom
- School of Mathematics, University of Bristol, Bristol BS8 1TW, United Kingdom
| | - D. V. Shalashilin
- School of Chemistry, University of Leeds, Leeds LS2 9JT, United Kingdom
| | - S. Fernandez-Alberti
- Departamento de Ciencia y Tecnologia, Universidad Nacional de Quilmes/CONICET, B1876BXD Bernal, Argentina
| |
Collapse
|
28
|
Baskov R, White AJ, Mozyrsky D. Improved Ehrenfest Approach to Model Correlated Electron-Nuclear Dynamics. J Phys Chem Lett 2019; 10:433-440. [PMID: 30621396 DOI: 10.1021/acs.jpclett.8b03061] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/09/2023]
Abstract
Mixed quantum-classical mechanical descriptions are critical to modeling coupled electron-nuclear dynamics, i.e., nonadiabatic molecular dynamics, relevant to photochemical and photophysical processes. We introduce an efficient description of such dynamics in terms of an effective Hamiltonian that not only properly captures electron-nuclear correlation effects but also helps develop an efficient computational method. In particular, we introduce a coupled Gaussian wavepacket parametrization of the nuclear wave function, which generalizes the Ehrenfest approach to account for electron-nuclei correlations. We test this new approach, Ehrenfest-Plus, on a suite of model problems that probe electron-nuclear correlation in nonadiabatic transitions. The high accuracy of our approach, combined with mixed quantum-classical efficiency, opens a path for improved simulation of nonadiabatic molecular dynamics in realistic molecular systems.
Collapse
Affiliation(s)
- Roman Baskov
- Institute of Physics of the National Academy of Sciences of Ukraine , Pr. Nauky 46 , Kyiv-28 MSP 03028 , Ukraine
| | - Alexander J White
- Theoretical Division , Los Alamos National Laboratory , Los Alamos , New Mexico 87545 , United States
| | - Dmitry Mozyrsky
- Theoretical Division , Los Alamos National Laboratory , Los Alamos , New Mexico 87545 , United States
| |
Collapse
|