1
|
S Mattos R, Mukherjee S, Barbatti M. Quantum Dynamics from Classical Trajectories. J Chem Theory Comput 2024. [PMID: 39235064 DOI: 10.1021/acs.jctc.4c00783] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 09/06/2024]
Abstract
Nonadiabatic molecular dynamics plays an essential role in exploring the time evolution of molecular systems. Various methods have been developed for this study, with varying accuracy and computational cost. One very successful among them is trajectory surface hopping, which propagates nuclei as classical trajectories using forces from a quantum description of the electrons and incorporates nonadiabatic effects through stochastic state changes during each trajectory propagation. A statistical analysis of an ensemble of the independent trajectories recovers the simulated system's behavior. This approach can give good results, but it is known to overlook nuclear quantum effects, leading to inaccurate predictions. Here, we present quantum dynamics from classical trajectories (QDCT), a new protocol to recover the quantum wavepacket from the classical trajectories generated by surface hopping. In this first QDCT implementation, we apply it to recover results at the multiple spawning level from postprocessing surface hopping precomputed trajectories. With a series of examples, we demonstrate QDCT's potential to improve the accuracy of the dynamics, correct decoherence effects, and diagnose problems or increase confidence in surface hopping results. All that comes at virtually no computational cost since no new electronic calculation is required.
Collapse
Affiliation(s)
- Rafael S Mattos
- Aix Marseille University, CNRS, ICR, 13397 Marseille, France
| | - Saikat Mukherjee
- Aix Marseille University, CNRS, ICR, 13397 Marseille, France
- Faculty of Chemistry, Nicolaus Copernicus University in Torun, 87100 Torun, Poland
| | - Mario Barbatti
- Aix Marseille University, CNRS, ICR, 13397 Marseille, France
- Institut Universitaire de France, 75231 Paris, France
| |
Collapse
|
2
|
Ibele LM, Memhood A, Levine BG, Avagliano D. Ab Initio Multiple Spawning Nonadiabatic Dynamics with Different CASPT2 Flavors: A Fully Open-Source PySpawn/OpenMolcas Interface. J Chem Theory Comput 2024. [PMID: 39228232 DOI: 10.1021/acs.jctc.4c00855] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 09/05/2024]
Abstract
We introduce an interface between PySpawn, a simulation package to run ab initio multiple spawning (AIMS) nonadiabatic dynamics, and OpenMolcas, a software package to perform multiconfigurational perturbations theory (CASPT2) electronic structure calculations. Our interface allows us to exploit all the functionalities of the two codes: the modular and efficient Python implementation of the AIMS algorithm and the extensive analysis tools offered by PySpawn, with the cutting-edge implementation of CASPT2 equations in OpenMolcas, including the recently introduced analytical gradients and different flavors. Both are fully open-source and free of charge, making the following implementation unique in the current plethora of software for nonadiabatic dynamics. This represents an important step toward a wider application of AIMS-based nonadiabatic dynamics combined with high-accuracy excited-state calculations. The importance and the need for such an implementation are demonstrated by application to the ultrafast relaxation of fulvene from S1 to S0, which is drastically affected by the potential energy surface on which the nuclear wavepacket is propagated. Additionally, the decay is influenced by the CASPT2 flavor adopted, posing interesting questions in the choice of one over the other and opening the door to deeper studies on the effect of CASPT2 formulations in nonadiabatic dynamics.
Collapse
Affiliation(s)
- Lea M Ibele
- CNRS, Institut de Chimie Physique UMR8000, Université Paris-Saclay, 91405 Orsay, France
| | - Arshad Memhood
- Department of Chemistry, Institute for Advanced Computational Science, Stony Brook, New York 11794, United States
| | - Benjamin G Levine
- Department of Chemistry, Institute for Advanced Computational Science, Stony Brook, New York 11794, United States
| | - Davide Avagliano
- Chimie ParisTech, CNRS, Institute of Chemistry for Life and Health Sciences (iCLeHS UMR 8060), PSL University, 75005 Paris, France
| |
Collapse
|
3
|
Myers CA, Miyazaki K, Trepl T, Isborn CM, Ananth N. GPU-accelerated on-the-fly nonadiabatic semiclassical dynamics. J Chem Phys 2024; 161:084114. [PMID: 39193942 DOI: 10.1063/5.0223628] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/16/2024] [Accepted: 08/11/2024] [Indexed: 08/29/2024] Open
Abstract
GPU-accelerated on-the-fly nonadiabatic dynamics is enabled by interfacing the linearized semiclassical dynamics approach with the TeraChem electronic structure program. We describe the computational workflow of the "PySCES" code interface, a Python code for semiclassical dynamics with on-the-fly electronic structure, including parallelization over multiple GPU nodes. We showcase the abilities of this code and present timings for two benchmark systems: fulvene solvated in acetonitrile and a charge transfer system in which a photoexcited zinc-phthalocyanine donor transfers charge to a fullerene acceptor through multiple electronic states on an ultrafast timescale. Our implementation paves the way for an efficient semiclassical approach to model the nonadiabatic excited state dynamics of complex molecules, materials, and condensed phase systems.
Collapse
Affiliation(s)
- Christopher A Myers
- Department of Chemistry and Biochemistry, University of California Merced, Merced, California 95343, USA
| | - Ken Miyazaki
- Department of Chemistry, Northwestern University, Evanston, Illinois 60208, USA
| | - Thomas Trepl
- Theoretical Physics IV, University of Bayreuth, 95440 Bayreuth, Germany
| | - Christine M Isborn
- Department of Chemistry and Biochemistry, University of California Merced, Merced, California 95343, USA
| | - Nandini Ananth
- Department of Chemistry and Chemical Biology, Baker Laboratory, Cornell University, Ithaca, New York 14853, USA
| |
Collapse
|
4
|
Curchod BFE, Orr-Ewing AJ. Perspective on Theoretical and Experimental Advances in Atmospheric Photochemistry. J Phys Chem A 2024; 128:6613-6635. [PMID: 39021090 PMCID: PMC11331530 DOI: 10.1021/acs.jpca.4c03481] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/27/2024] [Revised: 07/02/2024] [Accepted: 07/03/2024] [Indexed: 07/20/2024]
Abstract
Research that explores the chemistry of Earth's atmosphere is central to the current understanding of global challenges such as climate change, stratospheric ozone depletion, and poor air quality in urban areas. This research is a synergistic combination of three established domains: earth observation, for example, using satellites, and in situ field measurements; computer modeling of the atmosphere and its chemistry; and laboratory measurements of the properties and reactivity of gas-phase molecules and aerosol particles. The complexity of the interconnected chemical and photochemical reactions which determine the composition of the atmosphere challenges the capacity of laboratory studies to provide the spectroscopic, photochemical, and kinetic data required for computer models. Here, we consider whether predictions from computational chemistry using modern electronic structure theory and nonadiabatic dynamics simulations are becoming sufficiently accurate to supplement quantitative laboratory data for wavelength-dependent absorption cross-sections, photochemical quantum yields, and reaction rate coefficients. Drawing on presentations and discussions from the CECAM workshop on Theoretical and Experimental Advances in Atmospheric Photochemistry held in March 2024, we describe key concepts in the theory of photochemistry, survey the state-of-the-art in computational photochemistry methods, and compare their capabilities with modern experimental laboratory techniques. From such considerations, we offer a perspective on the scope of computational (photo)chemistry methods based on rigorous electronic structure theory to become a fourth core domain of research in atmospheric chemistry.
Collapse
|
5
|
Vandaele E, Mališ M, Luber S. A Local Diabatisation Method for Two-State Adiabatic Conical Intersections. J Chem Theory Comput 2024; 20:856-872. [PMID: 38174710 DOI: 10.1021/acs.jctc.3c01008] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/05/2024]
Abstract
A methodology to locally characterize conical intersections (CIs) between two adiabatic electronic states for which no nonadiabatic coupling (NAC) vectors are available is presented. Based on the Hessian and gradient at the CI, the branching space coordinates are identified. The potential energy surface around the CI in the branching space is expressed in the diabatic representation, from which the NAC vectors can be calculated in a wave-function-free, energy-based approach. To demonstrate the universality of the developed methodology, the minimum-energy CI (MECI) between the first (S1) and second (S2) singlet excited states of formamide is investigated at the state-averaged complete active space self-consistent field (SA-CASSCF) and extended multistate complete active space second-order perturbation theory (XMS-CASPT2) levels of theory. In addition, the asymmetrical MECI between the ground state (S0) and S1 of cyclopropanone is evaluated using SA-CASSCF, as well as (ME)CIs between the S1 and S2 states of benzene using SA-CASSCF and time-dependent density functional theory (TDDFT). Finally, a CI between the S1 and S2 excited states of thiophene was analyzed using TDDFT.
Collapse
Affiliation(s)
- Eva Vandaele
- Department of Chemistry, University of Zürich, Winterthurerstrasse 190, 8057 Zürich, Switzerland
| | - Momir Mališ
- Department of Chemistry, University of Zürich, Winterthurerstrasse 190, 8057 Zürich, Switzerland
| | - Sandra Luber
- Department of Chemistry, University of Zürich, Winterthurerstrasse 190, 8057 Zürich, Switzerland
| |
Collapse
|
6
|
Janoš J, Slavíček P. What Controls the Quality of Photodynamical Simulations? Electronic Structure Versus Nonadiabatic Algorithm. J Chem Theory Comput 2023; 19:8273-8284. [PMID: 37939301 PMCID: PMC10688183 DOI: 10.1021/acs.jctc.3c00908] [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: 08/18/2023] [Revised: 10/18/2023] [Accepted: 10/19/2023] [Indexed: 11/10/2023]
Abstract
The field of nonadiabatic dynamics has matured over the last decade with a range of algorithms and electronic structure methods available at the moment. While the community currently focuses more on developing and benchmarking new nonadiabatic dynamics algorithms, the underlying electronic structure controls the outcome of nonadiabatic simulations. Yet, the electronic-structure sensitivity analysis is typically neglected. In this work, we present a sensitivity analysis of the nonadiabatic dynamics of cyclopropanone to electronic structure methods and nonadiabatic dynamics algorithms. In particular, we compare wave function-based CASSCF, FOMO-CASCI, MS- and XMS-CASPT2, density-functional REKS, and semiempirical MRCI-OM3 electronic structure methods with the Landau-Zener surface hopping, fewest switches surface hopping, and ab initio multiple spawning with informed stochastic selection algorithms. The results clearly demonstrate that the electronic structure choice significantly influences the accuracy of nonadiabatic dynamics for cyclopropanone even when the potential energy surfaces exhibit qualitative and quantitative similarities. Thus, selecting the electronic structure solely on the basis of the mapping of potential energy surfaces can be misleading. Conversely, we observe no discernible differences in the performance of the nonadiabatic dynamics algorithms across the various methods. Based on the above results, we discuss the present-day practice in computational photodynamics.
Collapse
Affiliation(s)
- Jiří Janoš
- Department of Physical Chemistry, University of Chemistry and Technology, Technická 5, 16628 Prague 6, Czech Republic
| | - Petr Slavíček
- Department of Physical Chemistry, University of Chemistry and Technology, Technická 5, 16628 Prague 6, Czech Republic
| |
Collapse
|
7
|
Wang Z, Jing R, Li Y, Song D, Wan Y, Fukui N, Shinokubo H, Kuang Z, Xia A. Intrinsic Photostability in Dithiolonaphthalimide Achieved by Disulfide Bond-Induced Excited-State Quenching. J Phys Chem Lett 2023; 14:8485-8492. [PMID: 37721763 DOI: 10.1021/acs.jpclett.3c02260] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 09/19/2023]
Abstract
Disulfide bridges common in proteins show excellent photostability achieved by ultrafast internal conversion and maintain the stability of the tertiary structure. When disulfide bonds exist in aromatic compounds, the rigid chemical structure may affect the cleavage and reforming dynamics of disulfide bonds. In this work, a model compound with a disulfide five-membered-ring structure, 4,5-dithiolo-N-(2,6-dimethylphenyl)-1,8-naphthalimide (DTDPNI), is selected to elaborate the effect of disulfide modification on the excited-state deactivation mechanism. Quantum chemical calculations show that the S-S stretching leads to a dramatic decrease in the energy gap between the S1 and S0 states, similar to the situation in 1,2-dithiane. Due to the efficient nonradiative process, the excited-state lifetime of DTDPNI resolved by ultrafast spectroscopy is determined to be ∼20 ps. It is found that the excellent photostability is achieved by ultrafast excited-state quenching induced by the S-S stretching, rather than the cleavage of the disulfide bond; even the disulfide bridge is in a very rigid aromatic molecular system.
Collapse
Affiliation(s)
- Zeming Wang
- State Key Laboratory of Information Photonic and Optical Communications, and School of Science, Beijing University of Posts and Telecommunications (BUPT), Beijing 100876, P. R. China
| | - Rui Jing
- State Key Laboratory of Information Photonic and Optical Communications, and School of Science, Beijing University of Posts and Telecommunications (BUPT), Beijing 100876, P. R. China
| | - Yang Li
- State Key Laboratory of Information Photonic and Optical Communications, and School of Science, Beijing University of Posts and Telecommunications (BUPT), Beijing 100876, P. R. China
| | - Di Song
- State Key Laboratory of Information Photonic and Optical Communications, and School of Science, Beijing University of Posts and Telecommunications (BUPT), Beijing 100876, P. R. China
| | - Yan Wan
- College of Chemistry, Beijing Normal University, Beijing 100875, P. R. China
| | - Norihito Fukui
- Department of Molecular and Macromolecular Chemistry, Graduate School of Engineering, Nagoya University, Furo-chi, Chikusa-ku, Nagoya 464-8603, Japan
| | - Hiroshi Shinokubo
- Department of Molecular and Macromolecular Chemistry, Graduate School of Engineering, Nagoya University, Furo-chi, Chikusa-ku, Nagoya 464-8603, Japan
| | - Zhuoran Kuang
- State Key Laboratory of Information Photonic and Optical Communications, and School of Science, Beijing University of Posts and Telecommunications (BUPT), Beijing 100876, P. R. China
| | - Andong Xia
- State Key Laboratory of Information Photonic and Optical Communications, and School of Science, Beijing University of Posts and Telecommunications (BUPT), Beijing 100876, P. R. China
| |
Collapse
|
8
|
Mandal A, Taylor MA, Weight BM, Koessler ER, Li X, Huo P. Theoretical Advances in Polariton Chemistry and Molecular Cavity Quantum Electrodynamics. Chem Rev 2023; 123:9786-9879. [PMID: 37552606 PMCID: PMC10450711 DOI: 10.1021/acs.chemrev.2c00855] [Citation(s) in RCA: 31] [Impact Index Per Article: 31.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/07/2022] [Indexed: 08/10/2023]
Abstract
When molecules are coupled to an optical cavity, new light-matter hybrid states, so-called polaritons, are formed due to quantum light-matter interactions. With the experimental demonstrations of modifying chemical reactivities by forming polaritons under strong light-matter interactions, theorists have been encouraged to develop new methods to simulate these systems and discover new strategies to tune and control reactions. This review summarizes some of these exciting theoretical advances in polariton chemistry, in methods ranging from the fundamental framework to computational techniques and applications spanning from photochemistry to vibrational strong coupling. Even though the theory of quantum light-matter interactions goes back to the midtwentieth century, the gaps in the knowledge of molecular quantum electrodynamics (QED) have only recently been filled. We review recent advances made in resolving gauge ambiguities, the correct form of different QED Hamiltonians under different gauges, and their connections to various quantum optics models. Then, we review recently developed ab initio QED approaches which can accurately describe polariton states in a realistic molecule-cavity hybrid system. We then discuss applications using these method advancements. We review advancements in polariton photochemistry where the cavity is made resonant to electronic transitions to control molecular nonadiabatic excited state dynamics and enable new photochemical reactivities. When the cavity resonance is tuned to the molecular vibrations instead, ground-state chemical reaction modifications have been demonstrated experimentally, though its mechanistic principle remains unclear. We present some recent theoretical progress in resolving this mystery. Finally, we review the recent advances in understanding the collective coupling regime between light and matter, where many molecules can collectively couple to a single cavity mode or many cavity modes. We also lay out the current challenges in theory to explain the observed experimental results. We hope that this review will serve as a useful document for anyone who wants to become familiar with the context of polariton chemistry and molecular cavity QED and thus significantly benefit the entire community.
Collapse
Affiliation(s)
- Arkajit Mandal
- Department
of Chemistry, University of Rochester, 120 Trustee Road, Rochester, New York 14627, United States
- Department
of Chemistry, Columbia University, New York, New York 10027, United States
| | - Michael A.D. Taylor
- The
Institute of Optics, Hajim School of Engineering, University of Rochester, Rochester, New York 14627, United States
| | - Braden M. Weight
- Department
of Physics and Astronomy, University of
Rochester, Rochester, New York 14627, United
States
| | - Eric R. Koessler
- Department
of Chemistry, University of Rochester, 120 Trustee Road, Rochester, New York 14627, United States
| | - Xinyang Li
- Department
of Chemistry, University of Rochester, 120 Trustee Road, Rochester, New York 14627, United States
- Theoretical
Division, Los Alamos National Laboratory, Los Alamos, New Mexico 87545, United States
| | - Pengfei Huo
- Department
of Chemistry, University of Rochester, 120 Trustee Road, Rochester, New York 14627, United States
- The
Institute of Optics, Hajim School of Engineering, University of Rochester, Rochester, New York 14627, United States
| |
Collapse
|
9
|
Lassmann Y, Curchod BFE. Probing the sensitivity of ab initio multiple spawning to its parameters. Theor Chem Acc 2023; 142:66. [PMID: 37520272 PMCID: PMC10382418 DOI: 10.1007/s00214-023-03004-w] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/23/2023] [Accepted: 07/13/2023] [Indexed: 08/01/2023]
Abstract
Full multiple spawning (FMS) offers a strategy to simulate the nonadiabatic dynamics of molecular systems by describing their nuclear wavefunctions by a linear combination of coupled trajectory basis functions (TBFs). Applying a series of controlled approximations to the full multiple spawning (FMS) equations leads to the ab initio multiple spawning (AIMS), which is compatible with an on-the-fly propagation of the TBFs and an accurate description of nonadiabatic processes. The AIMS strategy and its numerical implementations, however, rely on a series of user-defined parameters. Herein, we investigate the influence of these parameters on the electronic-state population of two molecular systems- trans-azomethane and a two-dimensional model of the butatriene cation. This work highlights the stability of AIMS with respect to most of its parameters, underlines the specific parameters that require particular attention from the user of the method, and offers prescriptions for an informed selection of their value. Supplementary Information The online version contains supplementary material available at 10.1007/s00214-023-03004-w.
Collapse
Affiliation(s)
- Yorick Lassmann
- Centre for Computational Chemistry, School of Chemistry, Cantock’s Close, University of Bristol, Bristol, BS8 1TS UK
| | - Basile F. E. Curchod
- Centre for Computational Chemistry, School of Chemistry, Cantock’s Close, University of Bristol, Bristol, BS8 1TS UK
| |
Collapse
|
10
|
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
|
11
|
Mukherjee S, Hazra S, Ghosh S, Mukherjee S, Adhikari S. Trajectory Surface Hopping vs. Quantum Scattering Calculations on D+ + H2 and H + H2+ Reactions using Ab Initio Surfaces and Couplings. Chem Phys 2022. [DOI: 10.1016/j.chemphys.2022.111588] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/16/2022]
|
12
|
Liang R, Bakhtiiari A. Effects of Enzyme-Ligand Interactions on the Photoisomerization of a Light-Regulated Chemotherapeutic Drug. J Phys Chem B 2022; 126:2382-2393. [PMID: 35297246 DOI: 10.1021/acs.jpcb.1c10819] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Molecular photoswitches permit using light to control protein activity with high spatiotemporal resolutions, thereby alleviating the side effects of conventional chemotherapy. However, due to the challenges in probing ultrafast photoisomerization reactions in biological environments, it remains elusive how the protein influences the photochemistry of the photoswitches, which hampers the rational design of light-regulated therapeutics. To overcome this challenge, we employed first-principles nonadiabatic dynamics simulations to characterize the photodynamics of the phototrexate (PTX), a recently developed photoswitchable anticancer chemotherapeutic that reversibly inhibits its target enzyme dihydrofolate reductase (DHFR). Our simulations show that the protein environment impedes the trans to cis photoisomerization of the PTX. The confinement in the ligand-binding cavity slows down the isomerization kinetics and quantum yield of the photoswitch by reshaping its conical intersection, increasing its excited-state free-energy barrier and quenching its local density fluctuations. Also, the protein environment results in a suboptimal binding mode of the photoproduct that needs to undergo large structural rearrangement to effectively inhibit the enzyme. Therefore, we predict that the PTX's trans → cis photoisomerization in solution precedes its binding with the protein, despite the favorable binding energy of the trans isomer. Our findings highlight the importance of the protein environment on the photochemical reactions of the molecular photoswitches. As such, our work represents an important step toward the rational design of light-regulated drugs in photopharmacology.
Collapse
Affiliation(s)
- Ruibin Liang
- Department of Chemistry and Biochemistry, Texas Tech University, Lubbock, Texas 79409, United States
| | - Amirhossein Bakhtiiari
- Department of Chemistry and Biochemistry, Texas Tech University, Lubbock, Texas 79409, United States
| |
Collapse
|
13
|
Vandaele E, Mališ M, Luber S. The photodissociation of solvated cyclopropanone and its hydrate explored via non-adiabatic molecular dynamics using ΔSCF. Phys Chem Chem Phys 2022; 24:5669-5679. [PMID: 35179527 PMCID: PMC8890323 DOI: 10.1039/d1cp05187c] [Citation(s) in RCA: 8] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
The decay of cyclopropanone is a typical example of a photodecomposition process. Ethylene and carbon monoxide are formed following the excitation to the first singlet excited state through a symmetrical or asymmetrical pathway. The results obtained with non-adiabatic molecular dynamics (NAMD) using the delta self-consistent field (ΔSCF) method correspond well to previous experimental and multireference theoretical studies carried out in the gas phase. Moreover, this efficient methodology allows NAMD simulations of cyclopropanone in aqueous solution to be performed, which reveal analogue deactivation mechanisms, but a shorter lifetime and reduced photodissociation as compared to the gas-phase. The excited state dynamics of cyclopropanone hydrate, an enzyme inhibitor, in an aqueous environment are reported as well. Cyclopropanone hydrate strongly interacts with the surrounding solvent via the formation of hydrogen bonds. Excitation to the first singlet excited state shows an asymmetric pathway with cyclopropanone hydrate and propionic acid as the main photoproducts. The lifetime and photodissociation of cyclopropanone are reduced in aqueous solution, while the excitation of solvated cyclopropanone hydrate yields a range of photoproducts.![]()
Collapse
Affiliation(s)
- Eva Vandaele
- Department of Chemistry, University of Zürich, Winterthurerstrasse 190, 8057 Zürich, Switzerland.
| | - Momir Mališ
- Department of Chemistry, University of Zürich, Winterthurerstrasse 190, 8057 Zürich, Switzerland.
| | - Sandra Luber
- Department of Chemistry, University of Zürich, Winterthurerstrasse 190, 8057 Zürich, Switzerland.
| |
Collapse
|
14
|
Abstract
In this paper, we discuss coupled-trajectory schemes for molecular-dynamics simulations of excited-state processes. New coupled-trajectory strategies to capture decoherence effects, revival of coherence and nonadiabatic interferences in long-time dynamics are proposed, and compared to independent-trajectory schemes. The working framework is provided by the exact factorization of the electron-nuclear wave function, and it exploits ideas emanating from various surface-hopping schemes. The new coupled-trajectory algorithms are tested on a one-dimensional two-state system using different model parameters which allow one to induce different dynamics. The benchmark is provided by the numerically exact solution of the time-dependent Schrödinger equation.
Collapse
Affiliation(s)
- Carlotta Pieroni
- CNRS, Institut de Chimie Physique UMR8000, Université Paris-Saclay, 91405 Orsay, France.,Dipartimento di Chimica e Chimica Industriale, Università di Pisa, via G. Moruzzi 13, 56124 Pisa, Italy
| | - Federica Agostini
- CNRS, Institut de Chimie Physique UMR8000, Université Paris-Saclay, 91405 Orsay, France
| |
Collapse
|
15
|
Guo H, Worth G, Domcke W. Quantum dynamics with ab initio potentials. J Chem Phys 2021; 155:080401. [PMID: 34470339 DOI: 10.1063/5.0066234] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Hua Guo
- Department of Chemistry and Chemical Biology, University of New Mexico, Albuquerque, New Mexico 87131, USA
| | - Graham Worth
- Department of Chemistry, University College London, London WC1H 0AJ, United Kingdom
| | - Wolfgang Domcke
- Department of Chemistry, Technical University of Munich, D-85747 Garching, Germany
| |
Collapse
|
16
|
Zobel JP, González L. The Quest to Simulate Excited-State Dynamics of Transition Metal Complexes. JACS AU 2021; 1:1116-1140. [PMID: 34467353 PMCID: PMC8397362 DOI: 10.1021/jacsau.1c00252] [Citation(s) in RCA: 24] [Impact Index Per Article: 8.0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 06/07/2021] [Indexed: 05/15/2023]
Abstract
This Perspective describes current computational efforts in the field of simulating photodynamics of transition metal complexes. We present the typical workflows and feature the strengths and limitations of the different contemporary approaches. From electronic structure methods suitable to describe transition metal complexes to approaches able to simulate their nuclear dynamics under the effect of light, we give particular attention to build a bridge between theory and experiment by critically discussing the different models commonly adopted in the interpretation of spectroscopic experiments and the simulation of particular observables. Thereby, we review all the studies of excited-state dynamics on transition metal complexes, both in gas phase and in solution from reduced to full dimensionality.
Collapse
Affiliation(s)
- J. Patrick Zobel
- Institute
of Theoretical Chemistry, Faculty of Chemistry, University of Vienna, Währingerstr. 19, 1090 Vienna Austria
| | - Leticia González
- Institute
of Theoretical Chemistry, Faculty of Chemistry, University of Vienna, Währingerstr. 19, 1090 Vienna Austria
- Vienna
Research Platform on Accelerating Photoreaction Discovery, University of Vienna, Währingerstr. 19, 1090 Vienna Austria
| |
Collapse
|
17
|
Vindel-Zandbergen P, Ibele LM, Ha JK, Min SK, Curchod BFE, Maitra NT. Study of the Decoherence Correction Derived from the Exact Factorization Approach for Nonadiabatic Dynamics. J Chem Theory Comput 2021; 17:3852-3862. [PMID: 34138553 PMCID: PMC8280698 DOI: 10.1021/acs.jctc.1c00346] [Citation(s) in RCA: 22] [Impact Index Per Article: 7.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/03/2023]
Abstract
![]()
We present a detailed
study of the decoherence correction to surface
hopping that was recently derived from the exact factorization approach.
Ab initio multiple spawning calculations that use the same initial
conditions and the same electronic structure method are used as a
reference for three molecules: ethylene, the methaniminium cation,
and fulvene, for which nonadiabatic dynamics follows a photoexcitation.
A comparison with the Granucci–Persico energy-based decoherence
correction and the augmented fewest-switches surface-hopping scheme
shows that the three decoherence-corrected methods operate on individual
trajectories in a qualitatively different way, but the results averaged
over trajectories are similar for these systems.
Collapse
Affiliation(s)
| | - Lea M Ibele
- Department of Chemistry, Durham University, South Road, Durham DH1 3LE, U.K
| | - Jong-Kwon Ha
- Department of Chemistry, Ulsan National Institute of Science and Technology (UNIST), 50 UNIST-gil, Ulsan 44919, Republic of Korea
| | - Seung Kyu Min
- Department of Chemistry, Ulsan National Institute of Science and Technology (UNIST), 50 UNIST-gil, Ulsan 44919, Republic of Korea
| | - Basile F E Curchod
- Department of Chemistry, Durham University, South Road, Durham DH1 3LE, U.K
| | - Neepa T Maitra
- Department of Physics, Rutgers University, Newark, New Jersey 07102, United States
| |
Collapse
|
18
|
Lassmann Y, Curchod BFE. AIMSWISS-Ab initio multiple spawning with informed stochastic selections. J Chem Phys 2021; 154:211106. [PMID: 34240975 DOI: 10.1063/5.0052118] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/31/2023] Open
Abstract
Ab initio multiple spawning (AIMS) offers a reliable strategy to describe the excited-state dynamics and nonadiabatic processes of molecular systems. AIMS represents nuclear wavefunctions as linear combinations of traveling, coupled Gaussians called trajectory basis functions (TBFs) and uses a spawning algorithm to increase as needed the size of this basis set during nonadiabatic transitions. While the success of AIMS resides in this spawning algorithm, the dramatic increase in TBFs generated by multiple crossings between electronic states can rapidly lead to intractable dynamics. In this Communication, we introduce a new flavor of AIMS, coined ab initio multiple spawning with informed stochastic selections (AIMSWISS), which proposes a parameter-free strategy to beat the growing number of TBFs in an AIMS dynamics while preserving its accurate description of nonadiabatic transitions. The performance of AIMSWISS is validated against the photodynamics of ethylene, cyclopropanone, and fulvene. This technique, built upon the recently developed stochastic-selection AIMS, is intended to serve as a computationally affordable starting point for multiple spawning simulations.
Collapse
Affiliation(s)
- Yorick Lassmann
- 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
|
19
|
Tang D, Shen L, Fang WH. Evaluation of mixed quantum-classical molecular dynamics on cis-azobenzene photoisomerization. Phys Chem Chem Phys 2021; 23:13951-13964. [PMID: 34142685 DOI: 10.1039/d1cp01374b] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
The quantitative prediction of nonadiabatic transitions between different electronic states is important to understand ultrafast processes in photochemistry. A variety of mixed quantum-classical molecular dynamics methods such as surface hopping and Ehrenfest mean-field have been developed. However, how to choose an appropriate one from a wide diversity of dynamics algorithms to study a realistic photochemical process is still unclear. In this work, we implemented 30 combinations of different mixed quantum-classical dynamics methods, including 24 surface hopping models with 8 decoherence corrections and 3 momentum rescaling strategies as well as 6 mean-field models. Then we performed numerical investigations by simulating the photoisomerization of cis-azobenzene combined with on-the-fly electronic structure calculations. Predictions of the S1 lifetime and the quantum yield of the photoproduct using different models are distinct. Surface hopping is more robust than mean-field in our test system. Moreover, the choice of momentum rescaling methods in surface hopping brings more significant changes than decoherence corrections, while a large discrepancy between simulation results with different mean-field algorithms has been observed.
Collapse
Affiliation(s)
- Diandong Tang
- Key Laboratory of Theoretical and Computational Photochemistry of Ministry of Education, College of Chemistry, Beijing Normal University, Beijing 100875, China.
| | - Lin Shen
- Key Laboratory of Theoretical and Computational Photochemistry of Ministry of Education, College of Chemistry, Beijing Normal University, Beijing 100875, China.
| | - Wei-Hai Fang
- Key Laboratory of Theoretical and Computational Photochemistry of Ministry of Education, College of Chemistry, Beijing Normal University, Beijing 100875, China.
| |
Collapse
|