1
|
Joy C, Mandal B, Bostan D, Dubernet ML, Babikov D. Mixed quantum/classical theory (MQCT) approach to the dynamics of molecule-molecule collisions in complex systems. Faraday Discuss 2024; 251:225-248. [PMID: 38770664 DOI: 10.1039/d3fd00166k] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/22/2024]
Abstract
We developed a general theoretical approach and a user-ready computer code that permit study of the dynamics of collisional energy transfer and ro-vibrational energy exchange in complex molecule-molecule collisions. The method is a mixture of classical and quantum mechanics. The internal ro-vibrational motion of collision partners is treated quantum mechanically using a time-dependent Schrödinger equation that captures many quantum phenomena including state quantization and zero-point energy, propensity and selection rules for state-to-state transitions, quantum symmetry and interference phenomena. A significant numerical speed up is obtained by describing the translational motion of collision partners classically, using the Ehrenfest mean-field trajectory approach. Within this framework a family of approximate methods for collision dynamics is developed. Several benchmark studies for diatomic and triatomic molecules, such as H2O and ND3 collided with He, H2 and D2, show that the results of MQCT are in good agreement with full-quantum calculations in a broad range of energies, especially at high collision energies where they become nearly identical to the full quantum results. Numerical efficiency of the method and massive parallelism of the MQCT code permit us to embrace some of the most complicated collisional systems ever studied, such as C6H6 + He, CH3COOH + He and H2O + H2O. Application of MQCT to the collisions of chiral molecules such as CH3CHCH2O + He, and to molecule-surface collisions is also possible and will be pursued in the future.
Collapse
Affiliation(s)
- Carolin Joy
- Chemistry Department, Wehr Chemistry Building, Marquette University, Milwaukee, Wisconsin 53201-1881, USA.
| | - Bikramaditya Mandal
- Chemistry Department, Wehr Chemistry Building, Marquette University, Milwaukee, Wisconsin 53201-1881, USA.
| | - Dulat Bostan
- Chemistry Department, Wehr Chemistry Building, Marquette University, Milwaukee, Wisconsin 53201-1881, USA.
| | - Marie-Lise Dubernet
- Observatoire de Paris, PSL University, Sorbonne Universite, CNRS, SYRTE, Paris, France
| | - Dmitri Babikov
- Chemistry Department, Wehr Chemistry Building, Marquette University, Milwaukee, Wisconsin 53201-1881, USA.
| |
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
|
Galiana J, Lasorne B. Excitation energy transfer and vibronic relaxation through light-harvesting dendrimer building blocks: A nonadiabatic perspective. J Chem Phys 2024; 160:104104. [PMID: 38465688 DOI: 10.1063/5.0193264] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/21/2023] [Accepted: 02/20/2024] [Indexed: 03/12/2024] Open
Abstract
The light-harvesting excitonic properties of poly(phenylene ethynylene) (PPE) extended dendrimers (tree-like π-conjugated macromolecules) involve a directional cascade of local excitation energy transfer (EET) processes occurring from the "leaves" (shortest branches) to the "trunk" (longest branch), which can be viewed from a vibronic perspective as a sequence of internal conversions occurring among a connected graph of nonadiabatically coupled locally excited electronic states via conical intersections. The smallest PPE building block that is able to exhibit EET, the asymmetrically meta-substituted PPE oligomer with one acetylenic bond on one side and two parallel ones on the other side (hence, 2-ring and 3-ring para-substituted pseudo-fragments), is a prototype and the focus of the present work. From linear-response time-dependent density functional theory electronic-structure calculations of the molecule as regards its first two nonadiabatically coupled, optically active, singlet excited states, we built a (1 + 2)-state-8-dimensional vibronic-coupling Hamiltonian model for running subsequent multiconfiguration time-dependent Hartree wavepacket relaxations and propagations, yielding both steady-state absorption and emission spectra as well as real-time dynamics. The EET process from the shortest branch to the longest one occurs quite efficiently (about 80% quantum yield) within the first 25 fs after light excitation and is mediated vibrationally through acetylenic and quinoidal bond-stretching modes together with a particular role given to the central-ring anti-quinoidal rock-bending mode. Electronic and vibrational energy relaxations, together with redistributions of quantum populations and coherences, are interpreted herein through the lens of a nonadiabatic perspective, showing some interesting segregation among the foremost photoactive degrees of freedom as regards spectroscopy and reactivity.
Collapse
|
4
|
Titov E. Visible Light Induced Exciton Dynamics and trans-to- cis Isomerization in Azobenzene Aggregates: Insights from Surface Hopping/Semiempirical Configuration Interaction Molecular Dynamics Simulations. ACS OMEGA 2024; 9:8520-8532. [PMID: 38405525 PMCID: PMC10882624 DOI: 10.1021/acsomega.3c09900] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 12/11/2023] [Revised: 01/19/2024] [Accepted: 01/24/2024] [Indexed: 02/27/2024]
Abstract
Assemblies of photochromic molecules feature exciton states, which govern photochemical and photophysical processes in multichromophoric systems. Understanding the photoinduced dynamics of the assemblies requires nonadiabatic treatment involving multiple exciton states and numerous nuclear degrees of freedom, thus posing a challenge for simulations. In this work, we address this challenge for aggregates of azobenzene, a prototypical molecular switch, performing on-the-fly surface hopping calculations combined with semiempirical configuration interaction electronic structure and augmented with transition density matrix analysis to characterize exciton evolution. Specifically, we consider excitation of azobenzene tetramers in the nπ* absorption band located in the visible (blue) part of the electromagnetic spectrum, thus extending our recent work on dynamics after ππ* excitation corresponding to the ultraviolet region [Titov, J. Phys. Chem. C2023, 127, 13678-13688]. We find that the nπ* excitons, which are initially strongly localized by ground-state conformational disorder, undergo further (very strong) localization during short-time photodynamics. This excited-state localization process is extremely ultrafast, occurring within the first 10 fs of photodynamics. We observe virtually no exciton transfer of the localized excitons in the nπ* manifold. However, the transfer may occur via secondary pathways involving ππ* states or the ground state. Moreover, we find that the nπ* quantum yields of the trans-to-cis isomerization are reduced in the aggregated state.
Collapse
Affiliation(s)
- Evgenii Titov
- Institute of Chemistry, Theoretical
Chemistry, University of Potsdam, Karl-Liebknecht-Straße 24-25, 14476 Potsdam, Germany
| |
Collapse
|
5
|
Fischermeier D, Turkin A, Selby J, Lambert C, Mitrić R. Simulation of exciton spectra in disordered supramolecular polymers: exciton localization in cisoid indolenine squaraine hexamers. Phys Chem Chem Phys 2023; 26:219-229. [PMID: 38055887 DOI: 10.1039/d3cp04557a] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/08/2023]
Abstract
In order to understand the effects of disorder and defects in oligomers and polymers on the localization of excitons, we investigated the spectral properties of the squaraine B hexamer using long range corrected tight-binding TDDFT (lc-TDDFTB) and Frenkel-exciton model based calculations. Employing classical molecular dynamics, the cisoid indolenine squaraine hexamers helix was propagated in DCM and acetone to obtain ensembles of realistic structures, which naturally exhibit considerable disorder. The trajectories together with several model squaraine systems were studied to show the profound effects of disorder in the superstructure and disorder of the local monomer geometry on optical properties like absorption and exciton localization. We further compared lc-TDDFTB and exciton theory derived spectral data to related experimental data on absorption, exciton transfer and localization in squaraine polymers and oligomers.
Collapse
Affiliation(s)
- David Fischermeier
- Institut für Physikalische und Theoretische Chemie, Julius-Maximilians-Universität Würzburg, Emil-Fischer-Str. 42, 97074 Würzburg, Germany.
| | - Arthur Turkin
- Institut für Organische Chemie, Universitüt Würzburg, Am Hubland, D-97074 Würzburg, Germany.
| | - Joshua Selby
- Institut für Organische Chemie, Universitüt Würzburg, Am Hubland, D-97074 Würzburg, Germany.
| | - Christoph Lambert
- Institut für Organische Chemie, Universitüt Würzburg, Am Hubland, D-97074 Würzburg, Germany.
- Center for Nanosystems Chemistry, Universität Würzburg, Theodor-Boveri-Weg, D-97074 Würzburg, Germany
| | - Roland Mitrić
- Institut für Physikalische und Theoretische Chemie, Julius-Maximilians-Universität Würzburg, Emil-Fischer-Str. 42, 97074 Würzburg, Germany.
| |
Collapse
|
6
|
Mejia-Rodriguez D, Aprà E, Autschbach J, Bauman NP, Bylaska EJ, Govind N, Hammond JR, Kowalski K, Kunitsa A, Panyala A, Peng B, Rehr JJ, Song H, Tretiak S, Valiev M, Vila FD. NWChem: Recent and Ongoing Developments. J Chem Theory Comput 2023; 19:7077-7096. [PMID: 37458314 DOI: 10.1021/acs.jctc.3c00421] [Citation(s) in RCA: 4] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/25/2023]
Abstract
This paper summarizes developments in the NWChem computational chemistry suite since the last major release (NWChem 7.0.0). Specifically, we focus on functionality, along with input blocks, that is accessible in the current stable release (NWChem 7.2.0) and in the "master" development branch, interfaces to quantum computing simulators, interfaces to external libraries, the NWChem github repository, and containerization of NWChem executable images. Some ongoing developments that will be available in the near future are also discussed.
Collapse
Affiliation(s)
- Daniel Mejia-Rodriguez
- Physical Sciences Division, Pacific Northwest National Laboratory, Richland, Washington 99352, United States
| | - Edoardo Aprà
- Environmental Molecular Sciences Laboratory, Pacific Northwest National Laboratory, Richland, Washington 99352, United States
| | - Jochen Autschbach
- Department of Chemistry, University at Buffalo, State University of New York, Buffalo, New York 14260-3000, United States
| | - Nicholas P Bauman
- Physical Sciences Division, Pacific Northwest National Laboratory, Richland, Washington 99352, United States
| | - Eric J Bylaska
- Physical Sciences Division, Pacific Northwest National Laboratory, Richland, Washington 99352, United States
| | - Niranjan Govind
- Physical Sciences Division, Pacific Northwest National Laboratory, Richland, Washington 99352, United States
| | - Jeff R Hammond
- Accelerated Computing, NVIDIA Helsinki Oy, Porkkalankatu 1, 00180 Helsinki, Finland
| | - Karol Kowalski
- Physical Sciences Division, Pacific Northwest National Laboratory, Richland, Washington 99352, United States
| | - Alexander Kunitsa
- Zapata Computing, Inc., 100 Federal Street, Boston, Massachusetts 02110, United States
| | - Ajay Panyala
- Advanced Computing, Mathematics, and Data Division, Pacific Northwest National Laboratory, Richland, Washington 99352, United States
| | - Bo Peng
- Physical Sciences Division, Pacific Northwest National Laboratory, Richland, Washington 99352, United States
| | - John J Rehr
- Department of Physics, University of Washington, Seattle, Washington 98195, United States
| | - Huajing Song
- Physics and Chemistry of Materials, Theoretical Division, 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
- Center for Integrated Nanotechnologies, Los Alamos National Laboratory, Los Alamos, New Mexico 87545, United States
| | - Marat Valiev
- Physical Sciences Division, Pacific Northwest National Laboratory, Richland, Washington 99352, United States
| | - Fernando D Vila
- Department of Physics, University of Washington, Seattle, Washington 98195, United States
| |
Collapse
|
7
|
Gu B. A Discrete-Variable Local Diabatic Representation of Conical Intersection Dynamics. J Chem Theory Comput 2023; 19:6557-6563. [PMID: 37737832 DOI: 10.1021/acs.jctc.3c00560] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 09/23/2023]
Abstract
Conical intersections (CIs) are ubiquitous in polyatomic molecules and are responsible for a wide range of phenomena in photochemistry and photophysics. Modeling the conical intersection dynamics with adiabatic electronic states is hindered by the divergence of the first- and second-order derivative couplings at CIs due to electronic degeneracy. We introduce and implement a novel diabatic representation for exact correlated electron-nuclear wave packet dynamics through conical intersections. It directly employs the adiabatic electronic states but avoids the singular first- and second-order derivative couplings and is robust to different gauge choices of the electronic wave function phases. The reference nuclear geometries defining the adiabatic electronic states are determined by a discrete-variable representation of the nuclear coordinates. The nonadiabatic effects are accounted for by the electronic overlap matrix instead of derivative couplings as in the adiabatic representation. Illustrated by a two-mode conical intersection model, this representation captures all nonadiabatic effects, including electronic transitions, electronic coherence, and geometric phases. Thus, this representation provides a singularity-free framework for modeling ab initio conical intersection wave packet dynamics.
Collapse
Affiliation(s)
- Bing Gu
- Department of Chemistry & Department of Physics, Westlake University, Hangzhou, Zhejiang 310030, China
- Institute of Natural Sciences, Westlake Institute for Advanced Study, Hangzhou, Zhejiang 310024, China
| |
Collapse
|
8
|
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
|
9
|
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
|
10
|
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
|
11
|
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
|
12
|
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
|
13
|
Maskri R, Joubert-Doriol L. The moving crude adiabatic alternative to the adiabatic representation in excited state dynamics. PHILOSOPHICAL TRANSACTIONS. SERIES A, MATHEMATICAL, PHYSICAL, AND ENGINEERING SCIENCES 2022; 380:20200379. [PMID: 35341311 DOI: 10.1098/rsta.2020.0379] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 10/06/2021] [Accepted: 12/09/2021] [Indexed: 06/14/2023]
Abstract
The choice of the electronic representation in on-the-fly quantum dynamics is crucial. The adiabatic representation is appealing since adiabatic states are readily available from quantum chemistry packages. The nuclear wavepackets are then expanded in a basis of Gaussian functions, which follow trajectories to explore the potential energy surfaces and approximate the potential using a local expansion of the adiabatic quantities. Nevertheless, the adiabatic representation is plagued with severe limitations when conical intersections are involved: the diagonal Born-Oppenheimer corrections (DBOCs) are non-integrable, and the geometric phase effect on the nuclear wavepackets cannot be accounted for unless a model is available. To circumvent these difficulties, the moving crude adiabatic (MCA) representation was proposed and successfully tested in low energy dynamics where the wavepacket skirts the conical intersection. We assess the MCA representation in the case of non-adiabatic transitions through conical intersections. First, we show that using a Gaussian basis in the adiabatic representation indeed exhibits the aforementioned difficulties with a special emphasis on the possibility to regularize the DBOC terms. Then, we show that MCA is indeed able to properly model non-adiabatic transitions. Tests are done on linear vibronic coupling models for the bis(methylene) adamantyl cation and the butatriene cation. This article is part of the theme issue 'Chemistry without the Born-Oppenheimer approximation'.
Collapse
Affiliation(s)
- Rosa Maskri
- Univ Gustave Eiffel, Univ Paris Est Creteil, CNRS, UMR 8208, MSME, F-77454 Marne-la-Vallée, France
| | - Loïc Joubert-Doriol
- Univ Gustave Eiffel, Univ Paris Est Creteil, CNRS, UMR 8208, MSME, F-77454 Marne-la-Vallée, France
| |
Collapse
|
14
|
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
|
15
|
Bonardd S, Díaz Díaz D, Leiva A, Saldías C. Chromophoric Dendrimer-Based Materials: An Overview of Holistic-Integrated Molecular Systems for Fluorescence Resonance Energy Transfer (FRET) Phenomenon. Polymers (Basel) 2021; 13:4404. [PMID: 34960954 PMCID: PMC8705239 DOI: 10.3390/polym13244404] [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: 11/23/2021] [Revised: 12/09/2021] [Accepted: 12/10/2021] [Indexed: 12/15/2022] Open
Abstract
Dendrimers (from the Greek dendros → tree; meros → part) are macromolecules with well-defined three-dimensional and tree-like structures. Remarkably, this hyperbranched architecture is one of the most ubiquitous, prolific, and recognizable natural patterns observed in nature. The rational design and the synthesis of highly functionalized architectures have been motivated by the need to mimic synthetic and natural-light-induced energy processes. Dendrimers offer an attractive material scaffold to generate innovative, technological, and functional materials because they provide a high amount of peripherally functional groups and void nanoreservoirs. Therefore, dendrimers emerge as excellent candidates since they can play a highly relevant role as unimolecular reactors at the nanoscale, acting as versatile and sophisticated entities. In particular, they can play a key role in the properties of light-energy harvesting and non-radiative energy transfer, allowing them to function as a whole unit. Remarkably, it is possible to promote the occurrence of the FRET phenomenon to concentrate the absorbed energy in photoactive centers. Finally, we think an in-depth understanding of this mechanism allows for diverse and prolific technological applications, such as imaging, biomedical therapy, and the conversion and storage of light energy, among others.
Collapse
Affiliation(s)
- Sebastián Bonardd
- Departamento de Química Orgánica, Universidad de La Laguna, Avda. Astrofísico Francisco Sánchez S/N, La Laguna, 38206 Tenerife, Spain; (S.B.); (D.D.D.)
- Instituto Universitario de Bio-Orgánica Antonio González, Universidad de La Laguna, Avda. Astrofísico Francisco Sánchez 2, La Laguna, 38206 Tenerife, Spain
| | - David Díaz Díaz
- Departamento de Química Orgánica, Universidad de La Laguna, Avda. Astrofísico Francisco Sánchez S/N, La Laguna, 38206 Tenerife, Spain; (S.B.); (D.D.D.)
- Instituto Universitario de Bio-Orgánica Antonio González, Universidad de La Laguna, Avda. Astrofísico Francisco Sánchez 2, La Laguna, 38206 Tenerife, Spain
- Institutfür Organische Chemie, Universität Regensburg, Universitätsstr. 31, 93053 Regensburg, Germany
| | - Angel Leiva
- Departamento de Química Física, Facultad de Química y de Farmacia, Pontificia Universidad Católica de Chile, Macul, Santiago, CL 7820436, USA;
| | - César Saldías
- Departamento de Química Física, Facultad de Química y de Farmacia, Pontificia Universidad Católica de Chile, Macul, Santiago, CL 7820436, USA;
| |
Collapse
|
16
|
Zhang J, Peng J, Hu D, Lan Z. Investigation of nonadiabatic dynamics in the photolysis of methyl nitrate (CH 3ONO 2) by on-the-fly surface hopping simulation. Phys Chem Chem Phys 2021; 23:25597-25611. [PMID: 34546246 DOI: 10.1039/d1cp03226g] [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
The photolysis mechanism of methyl nitrate (CH3ONO2) was studied using the on-the-fly surface hopping dynamics at the XMS-CASPT2 level. Several critical geometries, including electronic state minima and conical intersections, were obtained, which play essential roles in the nonadiabatic dynamics of CH3ONO2. The ultrafast nonadiabatic decay dynamics to the ground state were simulated, which gives a proper explanation on the broad and structureless absorption spectra of CH3ONO2. The photodissociation channels, including CH3O + NO2, CH3O + NO + O, and others, as well as their branching ratios, were identified. When the dynamics starts from the lowest two electronic states (S1 and S2), the CH3O + NO2 channel is the dominant photolysis pathway, although we observed the minor contributions of other channels. In contrast, when the trajectories start from the third excited state S3, both CH3O + NO2 and CH3O + NO + O channels become important. Here the CH3O-NO2 bond dissociation takes place first, and then for some trajectories, the N-O bond of the NO2 part breaks successively. The quasi-degeneracy of electronic states may exist in the dissociation limits of both CH3O + NO2 and CH3O + NO + O channels. The current work provides valuable information in the understanding of experimental findings of the wavelength-dependent photolysis mechanism of CH3ONO2.
Collapse
Affiliation(s)
- Juanjuan Zhang
- 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 Environment, South China Normal University, Guangzhou 510006, China
| | - Deping Hu
- SCNU Environmental Research Institute, Guangdong Provincial Key Laboratory of Chemical Pollution and Environmental Safety & MOE Key Laboratory of Environmental Theoretical Chemistry, South China Normal University, Guangzhou 510006, China. .,School of Environment, South China Normal University, Guangzhou 510006, China
| | - Zhenggang Lan
- SCNU Environmental Research Institute, Guangdong Provincial Key Laboratory of Chemical Pollution and Environmental Safety & MOE Key Laboratory of Environmental Theoretical Chemistry, South China Normal University, Guangzhou 510006, China. .,School of Environment, South China Normal University, Guangzhou 510006, China
| |
Collapse
|
17
|
Ibele LM, Curchod BFE. Dynamics near a conical intersection-A diabolical compromise for the approximations of ab initio multiple spawning. J Chem Phys 2021; 155:174119. [PMID: 34742188 DOI: 10.1063/5.0071376] [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
Full multiple spawning (FMS) offers an exciting framework for the development of strategies to simulate the excited-state dynamics of molecular systems. FMS proposes to depict the dynamics of nuclear wavepackets by using a growing set of traveling multidimensional Gaussian functions called trajectory basis functions (TBFs). Perhaps the most recognized method emanating from FMS is the so-called ab initio multiple spawning (AIMS). In AIMS, the couplings between TBFs-in principle exact in FMS-are approximated to allow for the on-the-fly evaluation of required electronic-structure quantities. In addition, AIMS proposes to neglect the so-called second-order nonadiabatic couplings and the diagonal Born-Oppenheimer corrections. While AIMS has been applied successfully to simulate the nonadiabatic dynamics of numerous complex molecules, the direct influence of these missing or approximated terms on the nonadiabatic dynamics when approaching and crossing a conical intersection remains unknown to date. It is also unclear how AIMS could incorporate geometric-phase effects in the vicinity of a conical intersection. In this work, we assess the performance of AIMS in describing the nonadiabatic dynamics through a conical intersection for three two-dimensional, two-state systems that mimic the excited-state dynamics of bis(methylene)adamantyl, butatriene cation, and pyrazine. The population traces and nuclear density dynamics are compared with numerically exact quantum dynamics and trajectory surface hopping results. We find that AIMS offers a qualitatively correct description of the dynamics through a conical intersection for the three model systems. However, any attempt at improving the AIMS results by accounting for the originally neglected second-order nonadiabatic contributions appears to be stymied by the hermiticity requirement of the AIMS Hamiltonian and the independent first-generation approximation.
Collapse
Affiliation(s)
- Lea M Ibele
- Department of Chemistry, Durham University, South Road, Durham DH1 3LE, United Kingdom
| | - Basile F E Curchod
- Department of Chemistry, Durham University, South Road, Durham DH1 3LE, United Kingdom
| |
Collapse
|
18
|
Hu D, Peng J, Chen L, Gelin MF, Lan Z. Spectral Fingerprint of Excited-State Energy Transfer in Dendrimers through Polarization-Sensitive Transient-Absorption Pump-Probe Signals: On-the-Fly Nonadiabatic Dynamics Simulations. J Phys Chem Lett 2021; 12:9710-9719. [PMID: 34590858 DOI: 10.1021/acs.jpclett.1c02640] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/13/2023]
Abstract
The time-resolved polarization-sensitive transient-absorption (TA) pump-probe (PP) spectra are simulated using on-the-fly surface-hopping nonadiabatic dynamics and the doorway-window representation of nonlinear spectroscopy. A dendrimer model system composed of two linear phenylene ethynylene units (2-ring and 3-ring) is taken as an example. The ground-state bleach (GSB), stimulated emission (SE), and excited-state absorption (ESA) contributions as well as the total TA PP signals are obtained and carefully analyzed. It is shown that intramolecular excited-state energy transfer from the 2-ring unit to the 3-ring unit can be conveniently identified by employing pump and probe pulses with different polarizations. Our results demonstrate that time-resolved polarization-sensitive TA PP signals provide a powerful tool for the elucidation of excited-state energy-transfer pathways, notably in molecular systems possessing several optically bright nonadiabatically coupled electronic states with different orientations of transition dipole moments.
Collapse
Affiliation(s)
- Deping Hu
- 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 Environment, South China Normal University, Guangzhou 510006, China
| | - Lipeng Chen
- Max Planck Institute for the Physics of Complex Systems, 01187 Dresden, Germany
| | - Maxim F Gelin
- School of Sciences, Hangzhou Dianzi University, Hangzhou 310018, 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
|
19
|
Weight BM, Mandal A, Huo P. Ab initio symmetric quasi-classical approach to investigate molecular Tully models. J Chem Phys 2021; 155:084106. [PMID: 34470343 DOI: 10.1063/5.0061934] [Citation(s) in RCA: 9] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
We perform on-the-fly non-adiabatic molecular dynamics simulations using the symmetrical quasi-classical (SQC) approach with the recently suggested molecular Tully models: ethylene and fulvene. We attempt to provide benchmarks of the SQC methods using both the square and triangle windowing schemes as well as the recently proposed electronic zero-point-energy correction scheme (the so-called γ correction). We use the quasi-diabatic propagation scheme to directly interface the diabatic SQC methods with adiabatic electronic structure calculations. Our results showcase the drastic improvement of the accuracy by using the trajectory-adjusted γ-corrections, which outperform the widely used trajectory surface hopping method with decoherence corrections. These calculations provide useful and non-trivial tests to systematically investigate the numerical performance of various diabatic quantum dynamics approaches, going beyond simple diabatic model systems that have been used as the major workhorse in the quantum dynamics field. At the same time, these available benchmark studies will also likely foster the development of new quantum dynamics approaches based on these techniques.
Collapse
Affiliation(s)
- Braden M Weight
- Department of Physics and Astronomy, University of Rochester, Rochester, New York 14627, USA
| | - Arkajit Mandal
- Department of Chemistry, University of Rochester, Rochester, New York 14627, USA
| | - Pengfei Huo
- Department of Chemistry, University of Rochester, Rochester, New York 14627, USA
| |
Collapse
|
20
|
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: 15] [Impact Index Per Article: 5.0] [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
|
21
|
Brian D, Liu Z, Dunietz BD, Geva E, Sun X. Three-state harmonic models for photoinduced charge transfer. J Chem Phys 2021; 154:174105. [PMID: 34241055 DOI: 10.1063/5.0050289] [Citation(s) in RCA: 9] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/29/2022] Open
Abstract
A widely used strategy for simulating the charge transfer between donor and acceptor electronic states in an all-atom anharmonic condensed-phase system is based on invoking linear response theory to describe the system in terms of an effective spin-boson model Hamiltonian. Extending this strategy to photoinduced charge transfer processes requires also taking into consideration the ground electronic state in addition to the excited donor and acceptor electronic states. In this paper, we revisit the problem of describing such nonequilibrium processes in terms of an effective three-state harmonic model. We do so within the framework of nonequilibrium Fermi's golden rule (NE-FGR) in the context of photoinduced charge transfer in the carotenoid-porphyrin-C60 (CPC60) molecular triad dissolved in explicit tetrahydrofuran (THF). To this end, we consider different ways for obtaining a three-state harmonic model from the equilibrium autocorrelation functions of the donor-acceptor, donor-ground, and acceptor-ground energy gaps, as obtained from all-atom molecular dynamics simulations of the CPC60/THF system. The quantum-mechanically exact time-dependent NE-FGR rate coefficients for two different charge transfer processes in two different triad conformations are then calculated using the effective three-state model Hamiltonians as well as a hierarchy of more approximate expressions that lead to the instantaneous Marcus theory limit. Our results show that the photoinduced charge transfer in CPC60/THF can be described accurately by the effective harmonic three-state models and that nuclear quantum effects are small in this system.
Collapse
Affiliation(s)
- Dominikus Brian
- Division of Arts and Sciences, NYU Shanghai, 1555 Century Avenue, Shanghai 200122, China
| | - Zengkui Liu
- Division of Arts and Sciences, NYU Shanghai, 1555 Century Avenue, Shanghai 200122, China
| | - Barry D Dunietz
- Department of Chemistry and Biochemistry, Kent State University, Kent, Ohio 44242, USA
| | - Eitan Geva
- Department of Chemistry, University of Michigan, Ann Arbor, Michigan 48109, USA
| | - Xiang Sun
- Division of Arts and Sciences, NYU Shanghai, 1555 Century Avenue, Shanghai 200122, China
| |
Collapse
|
22
|
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
|
23
|
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
|
24
|
Shu Y, Zhang L, Sun S, Truhlar DG. Time-Derivative Couplings for Self-Consistent Electronically Nonadiabatic Dynamics. J Chem Theory Comput 2020; 16:4098-4106. [PMID: 32456433 DOI: 10.1021/acs.jctc.0c00409] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Electronically nonadiabatic dynamics methods based on a self-consistent potential, such as semiclassical Ehrenfest and coherent switching with decay of mixing, have a number of advantages but are computationally slower than approximations based on an unaveraged potential because they require evaluation of all components of the nonadiabatic coupling vector. Here we introduce a new approximation to the self-consistent potential that does not have this computational drawback. The new approximation uses time-derivative couplings evaluated by overlap integrals of electronic wave functions to approximate the nonadiabatic coupling terms in the equations of motion. We present a numerical test of the method for ethylene that shows there is little loss of accuracy in the ensemble-averaged results. This new approximation to the self-consistent potential makes direct dynamics calculations with self-consistent potentials more efficient for complex systems and makes them practically affordable for some cases where the cost was previously too high.
Collapse
Affiliation(s)
- Yinan Shu
- Department of Chemistry and Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, United States
| | - Linyao Zhang
- Department of Chemistry and Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, United States.,School of Energy Science and Engineering, Harbin Institute of Technology, Harbin 150001, People's Republic of China
| | - Shaozeng Sun
- School of Energy Science and Engineering, Harbin Institute of Technology, Harbin 150001, People's Republic of China
| | - Donald G Truhlar
- Department of Chemistry and Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, United States
| |
Collapse
|
25
|
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
|
26
|
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
|
27
|
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: 213] [Impact Index Per Article: 53.3] [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
|
28
|
Wohlgemuth M, Mitrić R. Excitation energy transport in DNA modelled by multi-chromophoric field-induced surface hopping. Phys Chem Chem Phys 2020; 22:16536-16551. [DOI: 10.1039/d0cp02255a] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/01/2023]
Abstract
Revealing the extended excited state lifetime due to excitation energy transport in DNA by multi-chromophoric field-induced surface-hopping (McFISH).
Collapse
Affiliation(s)
- Matthias Wohlgemuth
- Institut für Physikalische und Theoretische Chemie
- Julius-Maximilians-Universität Würzburg
- 97074 Würzburg
- Germany
| | - Roland Mitrić
- Institut für Physikalische und Theoretische Chemie
- Julius-Maximilians-Universität Würzburg
- 97074 Würzburg
- Germany
| |
Collapse
|
29
|
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]
|
30
|
Zhou W, Mandal A, Huo P. Quasi-Diabatic Scheme for Nonadiabatic On-the-Fly Simulations. J Phys Chem Lett 2019; 10:7062-7070. [PMID: 31665889 DOI: 10.1021/acs.jpclett.9b02747] [Citation(s) in RCA: 21] [Impact Index Per Article: 4.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/15/2023]
Abstract
We use the quasi-diabatic (QD) propagation scheme to perform on-the-fly nonadiabatic simulations of the photodynamics of ethylene. The QD scheme enables a seamless interface between accurate diabatic-based quantum dynamics approaches and adiabatic electronic structure calculations, explicitly avoiding any efforts to construct global diabatic states or reformulate the diabatic dynamics approach to the adiabatic representation. Using the partial linearized path-integral approach and the symmetrical quasi-classical approach as the diabatic dynamics methods, the QD propagation scheme enables direct nonadiabatic simulation with complete active space self-consistent field on-the-fly electronic structure calculations. The population dynamics obtained from both approaches are in a close agreement with the quantum wavepacket-based method and outperform the widely used trajectory surface-hopping approach. Further analysis of the ethylene photodeactivation pathways demonstrates the correct predictions of competing processes of nonradiative relaxation mechanism through various conical intersections. This work provides the foundation of using accurate diabatic dynamics approaches and on-the-fly adiabatic electronic structure information to perform ab initio nonadiabatic simulation.
Collapse
Affiliation(s)
- Wanghuai Zhou
- Advanced Functional Material and Photoelectric Technology Research Institution, School of Science , Hubei University of Automotive Technology , Shiyan , Hubei 442002 , People's Republic of China
- Department of Chemistry , University of Rochester , 120 Trustee Road , Rochester , New York 14627 , United States
| | - Arkajit Mandal
- Department of Chemistry , University of Rochester , 120 Trustee Road , Rochester , New York 14627 , United States
| | - Pengfei Huo
- Department of Chemistry , University of Rochester , 120 Trustee Road , Rochester , New York 14627 , United States
| |
Collapse
|
31
|
Ye L, Xu C, Gu FL, Zhu C. Functional and Basis Set Dependence for Time-Dependent Density Functional Theory Trajectory Surface Hopping Molecular Dynamics: Cis-Azobenzene Photoisomerization. J Comput Chem 2019; 41:635-645. [PMID: 31743473 DOI: 10.1002/jcc.26116] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/11/2019] [Revised: 10/26/2019] [Accepted: 11/04/2019] [Indexed: 11/08/2022]
Abstract
Within three functionals (TD-B3LYP, TD-BHandHLYP, and TD-CAM-B3LYP) in combination with four basis sets (3-21g, 6-31g, 6-31g(d), and cc-pvdz), global switching (GS) trajectory surface hopping molecular dynamics has been performed for cis-to-trans azobenzene photoisomerization up to the S1 (nπ*) excitation. Although all the combinations show artificial double-cone structure of conical intersection between ground and first excited states, simulated quantum yields and lifetimes are in good agreement with one another; 0.6 (±5%) and 40.5 fs (±10%) by TD-B3LYP, 0.5 (±10%) and 35.5 fs (±4%) by TD-BHandHLYP, and 0.44 (±9%) and 35.2 fs (±10%) by TD-CAM-B3LYP. By analyzing distributions of excited-state population decays, hopping spots, and typical trajectories with performance of 12 functional/basis set combinations, it has been concluded that functional dependence for given basis set is slightly more sensitive than basis set dependence for given functional. The present GS on-the-fly time-dependent density functional theory (TDDFT) trajectory surface hopping simulation can provide practical benchmark guidelines for conical intersection driven excited-state molecular dynamics simulation involving in large complex system within ordinary TDDFT framework. © 2019 Wiley Periodicals, Inc.
Collapse
Affiliation(s)
- Linfeng Ye
- Key Laboratory of Theoretical Chemistry of Environment, Ministry of Education, School of Chemistry and Environment of South China Normal University, Guangzhou, 51006, People's Republic of China
| | - Chao Xu
- Key Laboratory of Theoretical Chemistry of Environment, Ministry of Education, School of Chemistry and Environment of South China Normal University, Guangzhou, 51006, People's Republic of China
| | - Feng Long Gu
- Key Laboratory of Theoretical Chemistry of Environment, Ministry of Education, School of Chemistry and Environment of South China Normal University, Guangzhou, 51006, People's Republic of China
| | - Chaoyuan Zhu
- Key Laboratory of Theoretical Chemistry of Environment, Ministry of Education, School of Chemistry and Environment of South China Normal University, Guangzhou, 51006, People's Republic of China.,Institute of Molecular Science and Department of Applied Chemistry, National Chiao Tung University, Hsinchu, 30010, Taiwan.,Key Center for Emergent Functional Matter Science, National Chiao Tung University, Hsinchu, 30010, Taiwan
| |
Collapse
|
32
|
Ho EKL, Lasorne B. Diabatic pseudofragmentation and nonadiabatic excitation-energy transfer in meta-substituted dendrimer building blocks. COMPUT THEOR CHEM 2019. [DOI: 10.1016/j.comptc.2019.03.013] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/27/2022]
|
33
|
Agostini F, Curchod BFE. Different flavors of nonadiabatic molecular dynamics. WILEY INTERDISCIPLINARY REVIEWS-COMPUTATIONAL MOLECULAR SCIENCE 2019. [DOI: 10.1002/wcms.1417] [Citation(s) in RCA: 42] [Impact Index Per Article: 8.4] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 01/03/2023]
Affiliation(s)
- Federica Agostini
- Laboratoire de Chimie Physique UMR 8000 CNRS/University Paris‐Sud Orsay France
| | | |
Collapse
|
34
|
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
|
35
|
Mandal A, Sandoval C. JS, Shakib FA, Huo P. Quasi-Diabatic Propagation Scheme for Direct Simulation of Proton-Coupled Electron Transfer Reaction. J Phys Chem A 2019; 123:2470-2482. [DOI: 10.1021/acs.jpca.9b00077] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/30/2022]
Affiliation(s)
- Arkajit Mandal
- Department of Chemistry, University of Rochester, 120 Trustee Road, Rochester, New York 14627, United States
| | - Juan S. Sandoval C.
- Department of Chemistry, University of Rochester, 120 Trustee Road, Rochester, New York 14627, United States
| | - Farnaz A. Shakib
- Department of Chemistry, University of Rochester, 120 Trustee Road, Rochester, New York 14627, United States
| | - Pengfei Huo
- Department of Chemistry, University of Rochester, 120 Trustee Road, Rochester, New York 14627, United States
| |
Collapse
|
36
|
Fedorov DA, Levine BG. A discontinuous basis enables numerically exact solution of the Schrödinger equation around conical intersections in the adiabatic representation. J Chem Phys 2019; 150:054102. [PMID: 30736673 DOI: 10.1063/1.5058268] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Solving the vibrational Schrödinger equation in the neighborhood of conical intersections in the adiabatic representation is a challenge. At the intersection point, first- and second-derivative nonadiabatic coupling matrix elements become singular, with the singularity in the second-derivative coupling (diagonal Born-Oppenheimer correction) being non-integrable. These singularities result from discontinuities in the vibronic functions associated with the individual adiabatic states, and our group has recently argued that these divergent matrix elements cancel when discontinuous adiabatic vibronic functions sum to a continuous total nonadiabatic wave function. Here we describe the realization of this concept: a novel scheme for the numerically exact solution of the Schrödinger equation in the adiabatic representation. Our approach is based on a basis containing functions that are discontinuous at the intersection point. We demonstrate that the individual adiabatic nuclear wave functions are themselves discontinuous at the intersection point. This proves that discontinuous basis functions are essential to any tractable method that solves the Schrödinger equation around conical intersections in the adiabatic representation with high numerical precision. We establish that our method provides numerically exact results by comparison to reference calculations performed in the diabatic representation. In addition, we quantify the energetic error associated with constraining the density to be zero at the intersection point, a natural approximation. Prospects for extending the present treatment of a two-dimensional model to systems of higher dimensionality are discussed.
Collapse
Affiliation(s)
- Dmitry A Fedorov
- 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
|
37
|
Floquet Hamiltonian for incorporating electronic excitation by a laser pulse into simulations of non-adiabatic dynamics. Chem Phys 2018. [DOI: 10.1016/j.chemphys.2018.07.048] [Citation(s) in RCA: 13] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/18/2022]
|
38
|
Wang R, Chernyak VY. Dynamical consequences of time-reversal symmetry for systems with odd number of electrons: Conical intersections, semiclassical dynamics, and topology. Chem Phys 2018. [DOI: 10.1016/j.chemphys.2018.09.026] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/26/2022]
|
39
|
Titov E, Humeniuk A, Mitrić R. Exciton localization in excited-state dynamics of a tetracene trimer: a surface hopping LC-TDDFTB study. Phys Chem Chem Phys 2018; 20:25995-26007. [PMID: 30298878 DOI: 10.1039/c8cp05240a] [Citation(s) in RCA: 21] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/09/2023]
Abstract
Excitons in the molecular aggregates of chromophores are key participants in important processes such as photosynthesis or the functioning of organic photovoltaic devices. Therefore, the exploration of exciton dynamics is crucial. Here we report on exciton localization during excited-state dynamics of the recently synthesized tetracene trimer [Liu et al., Org. Lett., 2017, 19, 580]. We employ the surface hopping approach to nonadiabatic molecular dynamics in conjunction with the long-range corrected time-dependent density functional tight binding (LC-TDDFTB) method [Humeniuk and Mitrić, Comput. Phys. Commun., 2017, 221, 174]. Utilizing a set of descriptors based on the transition density matrix, we perform comprehensive analysis of exciton dynamics. The obtained results reveal an ultrafast exciton localization to a single tetracene unit of the trimer during excited-state dynamics, along with exciton transfer between units.
Collapse
Affiliation(s)
- Evgenii Titov
- Institut für Physikalische und Theoretische Chemie, Julius-Maximilians-Universität Würzburg, Emil-Fischer-Straße 42, 97074 Würzburg, Germany.
| | | | | |
Collapse
|
40
|
Sandoval C. JS, Mandal A, Huo P. Symmetric quasi-classical dynamics with quasi-diabatic propagation scheme. J Chem Phys 2018; 149:044115. [DOI: 10.1063/1.5036787] [Citation(s) in RCA: 23] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/25/2022] Open
Affiliation(s)
| | - Arkajit Mandal
- Department of Chemistry, University of Rochester, 120 Trustee Road, Rochester, New York 14627, USA
| | - Pengfei Huo
- Department of Chemistry, University of Rochester, 120 Trustee Road, Rochester, New York 14627, USA
| |
Collapse
|
41
|
Joubert-Doriol L, Izmaylov AF. Variational nonadiabatic dynamics in the moving crude adiabatic representation: Further merging of nuclear dynamics and electronic structure. J Chem Phys 2018; 148:114102. [PMID: 29566517 DOI: 10.1063/1.5020655] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
A new methodology of simulating nonadiabatic dynamics using frozen-width Gaussian wavepackets within the moving crude adiabatic representation with the on-the-fly evaluation of electronic structure is presented. The main feature of the new approach is the elimination of any global or local model representation of electronic potential energy surfaces; instead, the electron-nuclear interaction is treated explicitly using the Gaussian integration. As a result, the new scheme does not introduce any uncontrolled approximations. The employed variational principle ensures the energy conservation and leaves the number of electronic and nuclear basis functions as the only parameter determining the accuracy. To assess performance of the approach, a model with two electronic and two nuclear spacial degrees of freedom containing conical intersections between potential energy surfaces has been considered. Dynamical features associated with nonadiabatic transitions and nontrivial geometric (or Berry) phases were successfully reproduced within a limited basis expansion.
Collapse
Affiliation(s)
- Loïc Joubert-Doriol
- Department of Physical and Environmental Sciences, University of Toronto Scarborough, Toronto, Ontario M1C 1A4, Canada and Chemical Physics Theory Group, Department of Chemistry, University of Toronto, Toronto, Ontario M5S 3H6, Canada
| | - Artur F Izmaylov
- Department of Physical and Environmental Sciences, University of Toronto Scarborough, Toronto, Ontario M1C 1A4, Canada and Chemical Physics Theory Group, Department of Chemistry, University of Toronto, Toronto, Ontario M5S 3H6, Canada
| |
Collapse
|
42
|
Affiliation(s)
- Basile F. E. Curchod
- Department of Chemistry, Durham University, South Road, Durham DH1 3LE, United Kingdom
| | - Todd J. Martínez
- Department of Chemistry and PULSE Institute, Stanford University, Stanford, California 94305, United States
- SLAC National Accelerator Laboratory, Menlo Park, California 94025, United States
| |
Collapse
|
43
|
Freixas VM, Fernandez-Alberti S, Makhov DV, Tretiak S, Shalashilin D. An ab initio multiple cloning approach for the simulation of photoinduced dynamics in conjugated molecules. Phys Chem Chem Phys 2018; 20:17762-17772. [DOI: 10.1039/c8cp02321b] [Citation(s) in RCA: 21] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/29/2022]
Abstract
Multidimensional wave function: a superposition of Gaussian coherent states guided by Ehrenfest trajectories suited to clone and swap their electronic amplitudes.
Collapse
Affiliation(s)
| | | | - Dmitry V. Makhov
- School of Chemistry
- University of Leeds
- Leeds LS2 9JT
- UK
- School of Mathematics
| | - Sergei Tretiak
- Center for Nonlinear Studies (CNLS), and Center for Integrated Nanotechnologies (CINT)
- Los Alamos National Laboratory
- Los Alamos
- USA
| | | |
Collapse
|
44
|
A b initio quantum direct dynamics simulations of ultrafast photochemistry with Multiconfigurational Ehrenfest approach. Chem Phys 2017. [DOI: 10.1016/j.chemphys.2017.04.003] [Citation(s) in RCA: 79] [Impact Index Per Article: 11.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/19/2022]
|
45
|
Izmaylov AF, Joubert-Doriol L. Quantum Nonadiabatic Cloning of Entangled Coherent States. J Phys Chem Lett 2017; 8:1793-1797. [PMID: 28375623 DOI: 10.1021/acs.jpclett.7b00596] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/07/2023]
Abstract
We propose a systematic approach to the basis set extension for nonadiabatic dynamics of entangled combination of nuclear coherent states (CSs) evolving according to the time-dependent variational principle (TDVP). The TDVP provides a rigorous framework for fully quantum nonadiabatic dynamics of closed systems; however, the quality of results strongly depends on available basis functions. Starting with a single nuclear CS replicated vertically on all electronic states, our approach clones this function when replicas of the CS on different electronic states experience increasingly different forces. Created clones move away from each other (decohere), extending the basis set. To determine a moment for cloning, we introduce generalized forces based on derivatives that maximally contribute to a variation of the total quantum action and thus account for entanglement of all basis functions.
Collapse
Affiliation(s)
- Artur F Izmaylov
- Department of Physical and Environmental Sciences, University of Toronto Scarborough , Toronto, Ontario M1C 1A4, Canada
- Chemical Physics Theory Group, Department of Chemistry, University of Toronto , Toronto, Ontario M5S 3H6, Canada
| | - Loïc Joubert-Doriol
- Department of Physical and Environmental Sciences, University of Toronto Scarborough , Toronto, Ontario M1C 1A4, Canada
- Chemical Physics Theory Group, Department of Chemistry, University of Toronto , Toronto, Ontario M5S 3H6, Canada
| |
Collapse
|
46
|
Joubert-Doriol L, Sivasubramanium J, Ryabinkin IG, Izmaylov AF. Topologically Correct Quantum Nonadiabatic Formalism for On-the-Fly Dynamics. J Phys Chem Lett 2017; 8:452-456. [PMID: 28036173 DOI: 10.1021/acs.jpclett.6b02660] [Citation(s) in RCA: 15] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/06/2023]
Abstract
On-the-fly quantum nonadiabatic dynamics for large systems greatly benefits from the adiabatic representation readily available from electronic structure programs. However, conical intersections frequently occurring in this representation introduce nontrivial geometric or Berry phases which require a special treatment for adequate modeling of the nuclear dynamics. We analyze two approaches for nonadiabatic dynamics using the time-dependent variational principle and the adiabatic representation. The first approach employs adiabatic electronic functions with global parametric dependence on the nuclear coordinates. The second approach uses adiabatic electronic functions obtained only at the centers of moving localized nuclear basis functions (e.g., frozen-width Gaussians). Unless a gauge transformation is used to enforce single-valued boundary conditions, the first approach fails to capture the geometric phase. In contrast, the second approach accounts for the geometric phase naturally because of the absence of the global nuclear coordinate dependence in the electronic functions.
Collapse
Affiliation(s)
- Loïc Joubert-Doriol
- Department of Physical and Environmental Sciences, University of Toronto Scarborough , Toronto, Ontario M1C 1A4, Canada
- Chemical Physics Theory Group, Department of Chemistry, University of Toronto , Toronto, Ontario M5S 3H6, Canada
| | - Janakan Sivasubramanium
- Department of Physical and Environmental Sciences, University of Toronto Scarborough , Toronto, Ontario M1C 1A4, Canada
| | - Ilya G Ryabinkin
- Department of Physical and Environmental Sciences, University of Toronto Scarborough , Toronto, Ontario M1C 1A4, Canada
- Chemical Physics Theory Group, Department of Chemistry, University of Toronto , Toronto, Ontario M5S 3H6, Canada
| | - Artur F Izmaylov
- Department of Physical and Environmental Sciences, University of Toronto Scarborough , Toronto, Ontario M1C 1A4, Canada
- Chemical Physics Theory Group, Department of Chemistry, University of Toronto , Toronto, Ontario M5S 3H6, Canada
| |
Collapse
|
47
|
Heaps CW, Mazziotti DA. Accurate non-adiabatic quantum dynamics from pseudospectral sampling of time-dependent Gaussian basis sets. J Chem Phys 2016. [DOI: 10.1063/1.4959872] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Charles W. Heaps
- Department of Chemistry and The James Franck Institute, The University of Chicago, Chicago, Illinois 60637, USA
| | - David A. Mazziotti
- Department of Chemistry and The James Franck Institute, The University of Chicago, Chicago, Illinois 60637, USA
| |
Collapse
|