1
|
Li S, Xie BB, Yin BW, Liu L, Shen L, Fang WH. Construction of Highly Accurate Machine Learning Potential Energy Surfaces for Excited-State Dynamics Simulations Based on Low-Level Data Sets. J Phys Chem A 2024; 128:5516-5524. [PMID: 38954640 DOI: 10.1021/acs.jpca.4c02028] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 07/04/2024]
Abstract
Machine learning is capable of effectively predicting the potential energies of molecules in the presence of high-quality data sets. Its application in the construction of ground- and excited-state potential energy surfaces is attractive to accelerate nonadiabatic molecular dynamics simulations of photochemical reactions. Because of the huge computational cost of excited-state electronic structure calculations, the construction of a high-quality data set becomes a bottleneck. In the present work, we first built two data sets. One was obtained from surface hopping dynamics simulations at the semiempirical OM2/MRCI level. Another was extracted from the dynamics trajectories at the CASSCF level, which was reported previously. The ground- and excited-state potential energy surfaces of ethylene-bridged azobenzene at the CASSCF computational level were constructed based on the former low-level data set. Although non-neural network machine learning methods can achieve good or modest performance during the training process, only neural network models provide reliable predictions on the latter external test data set. The BPNN and SchNet combined with the Δ-ML scheme and the force term in the loss functions are recommended for dynamics simulations. Then, we performed excited-state dynamics simulations of the photoisomerization of ethylene-bridged azobenzene on machine learning potential energy surfaces. Compared with the lifetimes of the first excited state (S1) estimated at different computational levels, our results on the E isomer are in good agreement with the high-level estimation. However, the overestimation of the Z isomer is unimproved. It suggests that smaller errors during the training process do not necessarily translate to more accurate predictions on high-level potential energies or better performance on nonadiabatic dynamics simulations, at least in the present case.
Collapse
Affiliation(s)
- Shuai Li
- Key Laboratory of Theoretical and Computational Photochemistry of Ministry of Education, College of Chemistry, Beijing Normal University, Beijing 100875, P. R. China
| | - Bin-Bin Xie
- Hangzhou Institute of Advanced Studies, Zhejiang Normal University, Hangzhou 311231, Zhejiang, P. R. China
| | - Bo-Wen Yin
- Hangzhou Institute of Advanced Studies, Zhejiang Normal University, Hangzhou 311231, Zhejiang, P. R. China
| | - Lihong Liu
- Key Laboratory of Theoretical and Computational Photochemistry of Ministry of Education, College of Chemistry, Beijing Normal University, Beijing 100875, P. R. China
| | - Lin Shen
- Key Laboratory of Theoretical and Computational Photochemistry of Ministry of Education, College of Chemistry, Beijing Normal University, Beijing 100875, P. R. China
- Yantai-Jingshi Institute of Material Genome Engineering, Yantai 265505, Shandong, P. R. China
| | - Wei-Hai Fang
- Key Laboratory of Theoretical and Computational Photochemistry of Ministry of Education, College of Chemistry, Beijing Normal University, Beijing 100875, P. R. China
- Shandong Laboratory of Yantai Advanced Materials and Green Manufacturing, Yantai 264006, Shandong, P. R. China
| |
Collapse
|
2
|
Polonius S, Lehrner D, González L, Mai S. Resolving Photoinduced Femtosecond Three-Dimensional Solute-Solvent Dynamics through Surface Hopping Simulations. J Chem Theory Comput 2024; 20:4738-4750. [PMID: 38768386 PMCID: PMC11171268 DOI: 10.1021/acs.jctc.4c00169] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/09/2024] [Revised: 04/22/2024] [Accepted: 05/07/2024] [Indexed: 05/22/2024]
Abstract
Photoinduced dynamics in solution is governed by mutual solute-solvent interactions, which give rise to phenomena like solvatochromism, the Stokes shift, dual fluorescence, or charge transfer. Understanding these phenomena requires simulating the solute's photoinduced dynamics and simultaneously resolving the three-dimensional solvent distribution dynamics. If using trajectory surface hopping (TSH) to this aim, thousands of trajectories are required to adequately sample the time-dependent three-dimensional solvent distribution functions, and thus resolve the solvent dynamics with sub-Ångstrom and femtosecond accuracy and sufficiently low noise levels. Unfortunately, simulating thousands of trajectories with TSH in the framework of hybrid quantum mechanical/molecular mechanical (QM/MM) can be prohibitively expensive when employing ab initio electronic structure methods. To tackle this challenge, we recently introduced a computationally efficient approach that combines efficient linear vibronic coupling models with molecular mechanics (LVC/MM) via electrostatic embedding [Polonius et al., JCTC 2023, 19, 7171-7186]. This method provides solvent-embedded, nonadiabatically coupled potential energy surfaces while scaling similarly to MM force fields. Here, we employ TSH with LVC/MM to unravel the photoinduced dynamics of two small thiocarbonyl compounds solvated in water. We describe how to estimate the number of trajectories required to produce nearly noise-free three-dimensional solvent distribution functions and present an analysis based on approximately 10,000 trajectories propagated for 3 ps. In the electronic ground state, both molecules exhibit in-plane hydrogen bonds to the sulfur atom. Shortly after excitation, these bonds are broken and reform perpendicular to the molecular plane on timescales that differ by an order of magnitude due to steric effects. We also show that the solvent relaxation dynamics is coupled to the electronic dynamics, including intersystem crossing. These findings are relevant to advance the understanding of the coupled solute-solvent dynamics of solvated photoexcited molecules, e.g., biologically relevant thio-nucleobases.
Collapse
Affiliation(s)
- Severin Polonius
- Institute
of Theoretical Chemistry, Faculty of Chemistry, University of Vienna, Währinger Str. 17, 1090 Vienna, Austria
- Vienna
Doctoral School in Chemistry (DoSChem), University of Vienna, Währinger Str. 42, 1090 Vienna, Austria
| | - David Lehrner
- Institute
of Theoretical Chemistry, Faculty of Chemistry, University of Vienna, Währinger Str. 17, 1090 Vienna, Austria
| | - Leticia González
- Institute
of Theoretical Chemistry, Faculty of Chemistry, University of Vienna, Währinger Str. 17, 1090 Vienna, Austria
- Vienna
Research Platform on Accelerating Photoreaction Discovery, University of Vienna, Währinger Straße 17, 1090 Vienna, Austria
| | - Sebastian Mai
- Institute
of Theoretical Chemistry, Faculty of Chemistry, University of Vienna, Währinger Str. 17, 1090 Vienna, Austria
| |
Collapse
|
3
|
Papineau TV, Jacquemin D, Vacher M. Which Electronic Structure Method to Choose in Trajectory Surface Hopping Dynamics Simulations? Azomethane as a Case Study. J Phys Chem Lett 2024; 15:636-643. [PMID: 38205955 DOI: 10.1021/acs.jpclett.3c03014] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/12/2024]
Abstract
Nonadiabatic dynamics simulations have become a standard approach to explore photochemical reactions. Such simulations require underlying potential energy surfaces and couplings between them, calculated at a chosen level of theory, yet this aspect is rarely assessed. Here, in combination with the popular trajectory surface hopping dynamics method, we use a high-accuracy XMS-CASPT2 electronic structure level as a benchmark for assessing the performances of various post-Hartree-Fock methods (namely, CIS, ADC(2), CC2, and CASSCF) and exchange-correlation functionals (PBE, PBE0, and CAM-B3LYP) in a TD-DFT/TDA context, using the isomerization around a double bond as test case. Different relaxation pathways are identified, and the ability of the different methods to reproduce their relative importance and time scale is discussed. The results show that multireference electronic structure methods should be preferred, when studying nonadiabatic decay between excited and ground states. If not affordable, TD-DFT with TDA and hybrid functionals and ADC(2) are efficient alternatives but overestimate the nonradiative decay yield and thus may miss deexcitation pathways.
Collapse
Affiliation(s)
| | - Denis Jacquemin
- Nantes Université, CNRS, CEISAM UMR 6230, Nantes F-44000, France
- Institut Universitaire de France (IUF), Paris F-75005, France
| | - Morgane Vacher
- Nantes Université, CNRS, CEISAM UMR 6230, Nantes F-44000, France
| |
Collapse
|
4
|
Li Manni G, Fdez. Galván I, Alavi A, Aleotti F, Aquilante F, Autschbach J, Avagliano D, Baiardi A, Bao JJ, Battaglia S, Birnoschi L, Blanco-González A, Bokarev SI, Broer R, Cacciari R, Calio PB, Carlson RK, Carvalho Couto R, Cerdán L, Chibotaru LF, Chilton NF, Church JR, Conti I, Coriani S, Cuéllar-Zuquin J, Daoud RE, Dattani N, Decleva P, de Graaf C, Delcey M, De Vico L, Dobrautz W, Dong SS, Feng R, Ferré N, Filatov(Gulak) M, Gagliardi L, Garavelli M, González L, Guan Y, Guo M, Hennefarth MR, Hermes MR, Hoyer CE, Huix-Rotllant M, Jaiswal VK, Kaiser A, Kaliakin DS, Khamesian M, King DS, Kochetov V, Krośnicki M, Kumaar AA, Larsson ED, Lehtola S, Lepetit MB, Lischka H, López Ríos P, Lundberg M, Ma D, Mai S, Marquetand P, Merritt ICD, Montorsi F, Mörchen M, Nenov A, Nguyen VHA, Nishimoto Y, Oakley MS, Olivucci M, Oppel M, Padula D, Pandharkar R, Phung QM, Plasser F, Raggi G, Rebolini E, Reiher M, Rivalta I, Roca-Sanjuán D, Romig T, Safari AA, Sánchez-Mansilla A, Sand AM, Schapiro I, Scott TR, Segarra-Martí J, Segatta F, Sergentu DC, Sharma P, Shepard R, Shu Y, Staab JK, Straatsma TP, Sørensen LK, Tenorio BNC, Truhlar DG, Ungur L, Vacher M, Veryazov V, Voß TA, Weser O, Wu D, Yang X, Yarkony D, Zhou C, Zobel JP, Lindh R. The OpenMolcas Web: A Community-Driven Approach to Advancing Computational Chemistry. J Chem Theory Comput 2023; 19:6933-6991. [PMID: 37216210 PMCID: PMC10601490 DOI: 10.1021/acs.jctc.3c00182] [Citation(s) in RCA: 56] [Impact Index Per Article: 56.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/13/2023] [Indexed: 05/24/2023]
Abstract
The developments of the open-source OpenMolcas chemistry software environment since spring 2020 are described, with a focus on novel functionalities accessible in the stable branch of the package or via interfaces with other packages. These developments span a wide range of topics in computational chemistry and are presented in thematic sections: electronic structure theory, electronic spectroscopy simulations, analytic gradients and molecular structure optimizations, ab initio molecular dynamics, and other new features. This report offers an overview of the chemical phenomena and processes OpenMolcas can address, while showing that OpenMolcas is an attractive platform for state-of-the-art atomistic computer simulations.
Collapse
Affiliation(s)
- Giovanni Li Manni
- Electronic
Structure Theory Department, Max Planck
Institute for Solid State Research, Heisenbergstraße 1, 70569 Stuttgart, Germany
| | - Ignacio Fdez. Galván
- Department
of Chemistry − BMC, Uppsala University, P.O. Box 576, SE-75123 Uppsala, Sweden
| | - Ali Alavi
- Electronic
Structure Theory Department, Max Planck
Institute for Solid State Research, Heisenbergstraße 1, 70569 Stuttgart, Germany
- Yusuf Hamied
Department of Chemistry, University of Cambridge, Lensfield Road, Cambridge CB2 1EW, United Kingdom
| | - Flavia Aleotti
- Department
of Industrial Chemistry “Toso Montanari”, University of Bologna, 40136 Bologna, Italy
| | - Francesco Aquilante
- Theory and
Simulation of Materials (THEOS) and National Centre for Computational
Design and Discovery of Novel Materials (MARVEL), École Polytechnique Fédérale de Lausanne (EPFL), CH-1015 Lausanne, Switzerland
| | - Jochen Autschbach
- Department
of Chemistry, University at Buffalo, State
University of New York, Buffalo, New York 14260-3000, United States
| | - Davide Avagliano
- Department
of Industrial Chemistry “Toso Montanari”, University of Bologna, 40136 Bologna, Italy
| | - Alberto Baiardi
- ETH Zurich, Laboratory for Physical Chemistry, Vladimir-Prelog-Weg 2, 8093 Zurich, Switzerland
| | - Jie J. Bao
- Department
of Chemistry, Chemical Theory Center, and Minnesota Supercomputing
Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, United
States
| | - Stefano Battaglia
- Department
of Chemistry − BMC, Uppsala University, P.O. Box 576, SE-75123 Uppsala, Sweden
| | - Letitia Birnoschi
- The Department
of Chemistry, The University of Manchester, M13 9PL, Manchester, U.K.
| | - Alejandro Blanco-González
- Chemistry
Department, Bowling Green State University, Overmann Hall, Bowling Green, Ohio 43403, United States
| | - Sergey I. Bokarev
- Institut
für Physik, Universität Rostock, Albert-Einstein-Str. 23-24, 18059 Rostock, Germany
- Chemistry
Department, School of Natural Sciences, Technical University of Munich, Lichtenbergstr. 4, 85748 Garching, Germany
| | - Ria Broer
- Theoretical
Chemistry, Zernike Institute for Advanced Materials, University of Groningen, Nijenborgh 4, 9747AG Groningen, The Netherlands
| | - Roberto Cacciari
- Dipartimento
di Biotecnologie, Chimica e Farmacia, Università
di Siena, Via A. Moro 2, 53100 Siena, Italy
| | - Paul B. Calio
- Department
of Chemistry, Pritzker School of Molecular Engineering, James Franck
Institute, Chicago Center for Theoretical Chemistry, The University of Chicago, Chicago, Illinois 60637, United States
| | - Rebecca K. Carlson
- Department
of Chemistry, Chemical Theory Center, and Minnesota Supercomputing
Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, United
States
| | - Rafael Carvalho Couto
- Division
of Theoretical Chemistry and Biology, School of Engineering Sciences
in Chemistry, Biotechnology and Health, KTH Royal Institute of Technology, SE-106 91 Stockholm, Sweden
| | - Luis Cerdán
- Instituto
de Ciencia Molecular, Universitat de València, Catedrático José Beltrán
Martínez n. 2, 46980 Paterna, Spain
- Instituto
de Óptica (IO−CSIC), Consejo
Superior de Investigaciones Científicas, 28006, Madrid, Spain
| | - Liviu F. Chibotaru
- Department
of Chemistry, KU Leuven, Celestijnenlaan 200F, 3001 Leuven, Belgium
| | - Nicholas F. Chilton
- The Department
of Chemistry, The University of Manchester, M13 9PL, Manchester, U.K.
| | | | - Irene Conti
- Department
of Industrial Chemistry “Toso Montanari”, University of Bologna, 40136 Bologna, Italy
| | - Sonia Coriani
- Department
of Chemistry, Technical University of Denmark, Kemitorvet Bldg 207, 2800 Kongens Lyngby, Denmark
| | - Juliana Cuéllar-Zuquin
- Instituto
de Ciencia Molecular, Universitat de València, Catedrático José Beltrán
Martínez n. 2, 46980 Paterna, Spain
| | - Razan E. Daoud
- Dipartimento
di Biotecnologie, Chimica e Farmacia, Università
di Siena, Via A. Moro 2, 53100 Siena, Italy
| | - Nike Dattani
- HPQC Labs, Waterloo, N2T 2K9 Ontario Canada
- HPQC College, Waterloo, N2T 2K9 Ontario Canada
| | - Piero Decleva
- Istituto
Officina dei Materiali IOM-CNR and Dipartimento di Scienze Chimiche
e Farmaceutiche, Università degli
Studi di Trieste, I-34121 Trieste, Italy
| | - Coen de Graaf
- Department
of Physical and Inorganic Chemistry, Universitat
Rovira i Virgili, Tarragona 43007, Spain
- ICREA, Pg. Lluís
Companys 23, 08010 Barcelona, Spain
| | - Mickaël
G. Delcey
- Division
of Theoretical Chemistry and Biology, School of Engineering Sciences
in Chemistry, Biotechnology and Health, KTH Royal Institute of Technology, SE-106 91 Stockholm, Sweden
| | - Luca De Vico
- Dipartimento
di Biotecnologie, Chimica e Farmacia, Università
di Siena, Via A. Moro 2, 53100 Siena, Italy
| | - Werner Dobrautz
- Chalmers
University of Technology, Department of Chemistry
and Chemical Engineering, 41296 Gothenburg, Sweden
| | - Sijia S. Dong
- Department
of Chemistry, Chemical Theory Center, and Minnesota Supercomputing
Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, United
States
- Department
of Chemistry and Chemical Biology, Department of Physics, and Department
of Chemical Engineering, Northeastern University, Boston, Massachusetts 02115, United States
| | - Rulin Feng
- Department
of Chemistry, University at Buffalo, State
University of New York, Buffalo, New York 14260-3000, United States
- Department
of Chemistry, Fudan University, Shanghai 200433, China
| | - Nicolas Ferré
- Institut
de Chimie Radicalaire (UMR-7273), Aix-Marseille
Univ, CNRS, ICR 13013 Marseille, France
| | | | - Laura Gagliardi
- Department
of Chemistry, Chemical Theory Center, and Minnesota Supercomputing
Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, United
States
- Department
of Chemistry, Pritzker School of Molecular Engineering, James Franck
Institute, Chicago Center for Theoretical Chemistry, The University of Chicago, Chicago, Illinois 60637, United States
| | - Marco Garavelli
- Department
of Industrial Chemistry “Toso Montanari”, University of Bologna, 40136 Bologna, Italy
| | - Leticia González
- Institute
of Theoretical Chemistry, Faculty of Chemistry, University of Vienna, Währinger Straße 17, A-1090 Vienna, Austria
| | - Yafu Guan
- State Key
Laboratory of Molecular Reaction Dynamics and Center for Theoretical
Computational Chemistry, Dalian Institute
of Chemical Physics, Chinese Academy of Sciences, Dalian 116023, People’s Republic of China
| | - Meiyuan Guo
- SSRL, SLAC National Accelerator Laboratory, Menlo Park, California 94025, United States
| | - Matthew R. Hennefarth
- Department
of Chemistry, Pritzker School of Molecular Engineering, James Franck
Institute, Chicago Center for Theoretical Chemistry, The University of Chicago, Chicago, Illinois 60637, United States
| | - Matthew R. Hermes
- Department
of Chemistry, Chemical Theory Center, and Minnesota Supercomputing
Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, United
States
- Department
of Chemistry, Pritzker School of Molecular Engineering, James Franck
Institute, Chicago Center for Theoretical Chemistry, The University of Chicago, Chicago, Illinois 60637, United States
| | - Chad E. Hoyer
- Department
of Chemistry, Chemical Theory Center, and Minnesota Supercomputing
Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, United
States
- Department
of Chemistry, University of Washington, Seattle, Washington 98195, United States
| | - Miquel Huix-Rotllant
- Institut
de Chimie Radicalaire (UMR-7273), Aix-Marseille
Univ, CNRS, ICR 13013 Marseille, France
| | - Vishal Kumar Jaiswal
- Department
of Industrial Chemistry “Toso Montanari”, University of Bologna, 40136 Bologna, Italy
| | - Andy Kaiser
- Institut
für Physik, Universität Rostock, Albert-Einstein-Str. 23-24, 18059 Rostock, Germany
| | - Danil S. Kaliakin
- Chemistry
Department, Bowling Green State University, Overmann Hall, Bowling Green, Ohio 43403, United States
| | - Marjan Khamesian
- Department
of Chemistry − BMC, Uppsala University, P.O. Box 576, SE-75123 Uppsala, Sweden
| | - Daniel S. King
- Department
of Chemistry, Pritzker School of Molecular Engineering, James Franck
Institute, Chicago Center for Theoretical Chemistry, The University of Chicago, Chicago, Illinois 60637, United States
| | - Vladislav Kochetov
- Institut
für Physik, Universität Rostock, Albert-Einstein-Str. 23-24, 18059 Rostock, Germany
| | - Marek Krośnicki
- Institute
of Theoretical Physics and Astrophysics, Faculty of Mathematics, Physics
and Informatics, University of Gdańsk, ul Wita Stwosza 57, 80-952, Gdańsk, Poland
| | | | - Ernst D. Larsson
- Division
of Theoretical Chemistry, Chemical Centre, Lund University, P.O. Box 124, SE-22100, Lund, Sweden
| | - Susi Lehtola
- Molecular
Sciences Software Institute, Blacksburg, Virginia 24061, United States
- Department
of Chemistry, University of Helsinki, P.O. Box 55, FI-00014 University of Helsinki, Finland
| | - Marie-Bernadette Lepetit
- Condensed
Matter Theory Group, Institut Néel, CNRS UPR 2940, 38042 Grenoble, France
- Theory
Group, Institut Laue Langevin, 38042 Grenoble, France
| | - Hans Lischka
- Department
of Chemistry and Biochemistry, Texas Tech
University, Lubbock, Texas 79409-1061, United States
| | - Pablo López Ríos
- Electronic
Structure Theory Department, Max Planck
Institute for Solid State Research, Heisenbergstraße 1, 70569 Stuttgart, Germany
| | - Marcus Lundberg
- Department
of Chemistry − Ångström Laboratory, Uppsala University, SE-75120 Uppsala, Sweden
| | - Dongxia Ma
- Electronic
Structure Theory Department, Max Planck
Institute for Solid State Research, Heisenbergstraße 1, 70569 Stuttgart, Germany
- Department
of Chemistry, Chemical Theory Center, and Minnesota Supercomputing
Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, United
States
| | - Sebastian Mai
- Institute
of Theoretical Chemistry, Faculty of Chemistry, University of Vienna, Währinger Straße 17, A-1090 Vienna, Austria
| | - Philipp Marquetand
- Institute
of Theoretical Chemistry, Faculty of Chemistry, University of Vienna, Währinger Straße 17, A-1090 Vienna, Austria
| | | | - Francesco Montorsi
- Department
of Industrial Chemistry “Toso Montanari”, University of Bologna, 40136 Bologna, Italy
| | - Maximilian Mörchen
- ETH Zurich, Laboratory for Physical Chemistry, Vladimir-Prelog-Weg 2, 8093 Zurich, Switzerland
| | - Artur Nenov
- Department
of Industrial Chemistry “Toso Montanari”, University of Bologna, 40136 Bologna, Italy
| | - Vu Ha Anh Nguyen
- Department
of Chemistry, National University of Singapore, 3 Science Drive 3, 117543 Singapore
| | - Yoshio Nishimoto
- Graduate
School of Science, Kyoto University, Kyoto 606-8502, Japan
| | - Meagan S. Oakley
- Department
of Chemistry, Chemical Theory Center, and Minnesota Supercomputing
Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, United
States
| | - Massimo Olivucci
- Chemistry
Department, Bowling Green State University, Overmann Hall, Bowling Green, Ohio 43403, United States
- Dipartimento
di Biotecnologie, Chimica e Farmacia, Università
di Siena, Via A. Moro 2, 53100 Siena, Italy
| | - Markus Oppel
- Institute
of Theoretical Chemistry, Faculty of Chemistry, University of Vienna, Währinger Straße 17, A-1090 Vienna, Austria
| | - Daniele Padula
- Dipartimento
di Biotecnologie, Chimica e Farmacia, Università
di Siena, Via A. Moro 2, 53100 Siena, Italy
| | - Riddhish Pandharkar
- Department
of Chemistry, Chemical Theory Center, and Minnesota Supercomputing
Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, United
States
- Department
of Chemistry, Pritzker School of Molecular Engineering, James Franck
Institute, Chicago Center for Theoretical Chemistry, The University of Chicago, Chicago, Illinois 60637, United States
| | - Quan Manh Phung
- Department
of Chemistry, Graduate School of Science, Nagoya University, Furo-cho, Chikusa-ku, Nagoya, Aichi 464-8602, Japan
- Institute
of Transformative Bio-Molecules (WPI-ITbM), Nagoya University, Furo-cho, Chikusa-ku, Nagoya, Aichi 464-8601, Japan
| | - Felix Plasser
- Department
of Chemistry, Loughborough University, Loughborough, LE11 3TU, U.K.
| | - Gerardo Raggi
- Department
of Chemistry − BMC, Uppsala University, P.O. Box 576, SE-75123 Uppsala, Sweden
- Quantum
Materials and Software LTD, 128 City Road, London, EC1V 2NX, United Kingdom
| | - Elisa Rebolini
- Scientific
Computing Group, Institut Laue Langevin, 38042 Grenoble, France
| | - Markus Reiher
- ETH Zurich, Laboratory for Physical Chemistry, Vladimir-Prelog-Weg 2, 8093 Zurich, Switzerland
| | - Ivan Rivalta
- Department
of Industrial Chemistry “Toso Montanari”, University of Bologna, 40136 Bologna, Italy
| | - Daniel Roca-Sanjuán
- Instituto
de Ciencia Molecular, Universitat de València, Catedrático José Beltrán
Martínez n. 2, 46980 Paterna, Spain
| | - Thies Romig
- Institut
für Physik, Universität Rostock, Albert-Einstein-Str. 23-24, 18059 Rostock, Germany
| | - Arta Anushirwan Safari
- Electronic
Structure Theory Department, Max Planck
Institute for Solid State Research, Heisenbergstraße 1, 70569 Stuttgart, Germany
| | - Aitor Sánchez-Mansilla
- Department
of Physical and Inorganic Chemistry, Universitat
Rovira i Virgili, Tarragona 43007, Spain
| | - Andrew M. Sand
- Department
of Chemistry, Chemical Theory Center, and Minnesota Supercomputing
Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, United
States
- Department
of Chemistry and Biochemistry, Butler University, Indianapolis, Indiana 46208, United States
| | - Igor Schapiro
- Institute
of Chemistry, The Hebrew University of Jerusalem, Jerusalem 91904, Israel
| | - Thais R. Scott
- Department
of Chemistry, Chemical Theory Center, and Minnesota Supercomputing
Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, United
States
- Department
of Chemistry, Pritzker School of Molecular Engineering, James Franck
Institute, Chicago Center for Theoretical Chemistry, The University of Chicago, Chicago, Illinois 60637, United States
- Department
of Chemistry, University of California, Irvine, California 92697, United States
| | - Javier Segarra-Martí
- Instituto
de Ciencia Molecular, Universitat de València, Catedrático José Beltrán
Martínez n. 2, 46980 Paterna, Spain
| | - Francesco Segatta
- Department
of Industrial Chemistry “Toso Montanari”, University of Bologna, 40136 Bologna, Italy
| | - Dumitru-Claudiu Sergentu
- Department
of Chemistry, University at Buffalo, State
University of New York, Buffalo, New York 14260-3000, United States
- Laboratory
RA-03, RECENT AIR, A. I. Cuza University of Iaşi, RA-03 Laboratory (RECENT AIR), Iaşi 700506, Romania
| | - Prachi Sharma
- Department
of Chemistry, Chemical Theory Center, and Minnesota Supercomputing
Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, United
States
| | - Ron Shepard
- Chemical
Sciences and Engineering Division, Argonne
National Laboratory, Lemont, Illinois 60439, USA
| | - Yinan Shu
- Department
of Chemistry, Chemical Theory Center, and Minnesota Supercomputing
Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, United
States
| | - Jakob K. Staab
- The Department
of Chemistry, The University of Manchester, M13 9PL, Manchester, U.K.
| | - Tjerk P. Straatsma
- National
Center for Computational Sciences, Oak Ridge
National Laboratory, Oak Ridge, Tennessee 37831-6373, United States
- Department
of Chemistry and Biochemistry, University
of Alabama, Tuscaloosa, Alabama 35487-0336, United States
| | | | - Bruno Nunes Cabral Tenorio
- Department
of Chemistry, Technical University of Denmark, Kemitorvet Bldg 207, 2800 Kongens Lyngby, Denmark
| | - Donald G. Truhlar
- Department
of Chemistry, Chemical Theory Center, and Minnesota Supercomputing
Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, United
States
| | - Liviu Ungur
- Department
of Chemistry, National University of Singapore, 3 Science Drive 3, 117543 Singapore
| | - Morgane Vacher
- Nantes
Université, CNRS, CEISAM, UMR 6230, F-44000 Nantes, France
| | - Valera Veryazov
- Division
of Theoretical Chemistry, Chemical Centre, Lund University, P.O. Box 124, SE-22100, Lund, Sweden
| | - Torben Arne Voß
- Institut
für Physik, Universität Rostock, Albert-Einstein-Str. 23-24, 18059 Rostock, Germany
| | - Oskar Weser
- Electronic
Structure Theory Department, Max Planck
Institute for Solid State Research, Heisenbergstraße 1, 70569 Stuttgart, Germany
| | - Dihua Wu
- Department
of Chemistry, Chemical Theory Center, and Minnesota Supercomputing
Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, United
States
| | - Xuchun Yang
- Chemistry
Department, Bowling Green State University, Overmann Hall, Bowling Green, Ohio 43403, United States
| | - David Yarkony
- Department
of Chemistry, Johns Hopkins University, Baltimore, Maryland 21218, United States
| | - Chen Zhou
- Department
of Chemistry, Chemical Theory Center, and Minnesota Supercomputing
Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, United
States
| | - J. Patrick Zobel
- Institute
of Theoretical Chemistry, Faculty of Chemistry, University of Vienna, Währinger Straße 17, A-1090 Vienna, Austria
| | - Roland Lindh
- Department
of Chemistry − BMC, Uppsala University, P.O. Box 576, SE-75123 Uppsala, Sweden
- Uppsala
Center for Computational Chemistry (UC3), Uppsala University, PO Box 576, SE-751 23 Uppsala. Sweden
| |
Collapse
|
5
|
Polonius S, Zhuravel O, Bachmair B, Mai S. LVC/MM: A Hybrid Linear Vibronic Coupling/Molecular Mechanics Model with Distributed Multipole-Based Electrostatic Embedding for Highly Efficient Surface Hopping Dynamics in Solution. J Chem Theory Comput 2023; 19:7171-7186. [PMID: 37788824 PMCID: PMC10601485 DOI: 10.1021/acs.jctc.3c00805] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/25/2023] [Indexed: 10/05/2023]
Abstract
We present a theoretical framework for a hybrid linear vibronic coupling model electrostatically embedded into a molecular mechanics environment, termed the linear vibronic coupling/molecular mechanics (LVC/MM) method, for the surface hopping including arbitrary coupling (SHARC) molecular dynamics package. Electrostatic embedding is realized through the computation of interactions between environment point charges and distributed multipole expansions (DMEs, up to quadrupoles) that represent each electronic state and transition densities in the diabatic basis. The DME parameters are obtained through a restrained electrostatic potential (RESP) fit, which we extended to yield higher-order multipoles. We also implemented in SHARC a scheme for achieving roto-translational invariance of LVC models as well as a general quantum mechanics/molecular mechanics (QM/MM) interface, an OpenMM interface, and restraining potentials for simulating liquid droplets. Using thioformaldehyde in water as a test case, we demonstrate that LVC/MM can accurately reproduce the solvation structure and energetics of rigid solutes, with errors on the order of 1-2 kcal/mol compared to a BP86/MM reference. The implementation in SHARC is shown to be very efficient, enabling the simulation of trajectories on the nanosecond time scale in a matter of days.
Collapse
Affiliation(s)
- Severin Polonius
- Institute
of Theoretical Chemistry, Faculty of Chemistry, University of Vienna, Währinger Str. 17, 1090 Vienna, Austria
- Vienna
Doctoral School in Chemistry (DoSChem), University of Vienna, Währinger Str. 42, 1090 Vienna, Austria
| | - Oleksandra Zhuravel
- Institute
of Theoretical Chemistry, Faculty of Chemistry, University of Vienna, Währinger Str. 17, 1090 Vienna, Austria
| | - Brigitta Bachmair
- Vienna
Doctoral School in Chemistry (DoSChem), University of Vienna, Währinger Str. 42, 1090 Vienna, Austria
- Research
Platform on Accelerating Photoreaction Discovery (ViRAPID), University of Vienna, Währinger Str. 17, 1090 Vienna, Austria
| | - Sebastian Mai
- Institute
of Theoretical Chemistry, Faculty of Chemistry, University of Vienna, Währinger Str. 17, 1090 Vienna, Austria
| |
Collapse
|
6
|
Guan Y, Xie C, Guo H, Yarkony DR. Toward a Unified Analytical Description of Internal Conversion and Intersystem Crossing in the Photodissociation of Thioformaldehyde. I. Diabatic Singlet States. J Chem Theory Comput 2023; 19:6414-6424. [PMID: 37698839 DOI: 10.1021/acs.jctc.3c00628] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 09/13/2023]
Abstract
The photodissociation of thioformaldehyde is an archetypal system for the study of competition between internal conversion and intersystem crossing, which involves its two singlet states (S0 and S1) and two triplet states (T1 and T2). In order to perform accurate dynamic simulations, either quantum or quasi-classical, it is essential to construct an analytical representation for all necessary electronic structure data. In this work, a diabatic potential energy matrix (DPEM), Hd, for the two singlet states (S0 and S1) is reported. The analytical form of DPEM is symmetrized and constructed to reproduce adiabatic energies, energy gradients, and derivative couplings obtained from high-level multireference configuration interaction wave functions. The Hd is fully saturated in the molecular configuration space with a trajectory-guided point sampling approach. This Hd can provide the accurate description of the photodissociation of thioformaldehyde on its singlet states and is also a necessary part for incorporating the spin-orbit couplings into a unified diabatic framework. Preliminary quasi-classical trajectory simulations show that a roaming mechanism also exists in the molecular dissociation channel of thioformaldehyde.
Collapse
Affiliation(s)
- Yafu Guan
- State Key Laboratory of Molecular Reaction Dynamics and Center for Theoretical Computational Chemistry, Dalian Institute of Chemical Physics, Chinese Academy of Sciences, Dalian 116023, People's Republic of China
| | - Changjian Xie
- Institute of Modern Physics, Northwest University, Xi'an, Shaanxi 710069, People's Republic of China
| | - Hua Guo
- Department of Chemistry and Chemical Biology, University of New Mexico, Albuquerque, New Mexico 87131, United States
| | - David R Yarkony
- Department of Chemistry, Johns Hopkins University, Baltimore, Maryland 21218, United States
| |
Collapse
|
7
|
Hutton L, Curchod BFE. Photodynamics of Gas‐Phase Pyruvic Acid Following Light Absorption in the Actinic Region. CHEMPHOTOCHEM 2022. [DOI: 10.1002/cptc.202200151] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/09/2022]
Affiliation(s)
- Lewis Hutton
- Centre for Computational Chemistry School of Chemistry University of Bristol Bristol BS8 1TS UK
| | - Basile F. E. Curchod
- Centre for Computational Chemistry School of Chemistry University of Bristol Bristol BS8 1TS UK
| |
Collapse
|
8
|
Calio PB, Truhlar DG, Gagliardi L. Nonadiabatic Molecular Dynamics by Multiconfiguration Pair-Density Functional Theory. J Chem Theory Comput 2022; 18:614-622. [PMID: 35030306 DOI: 10.1021/acs.jctc.1c01048] [Citation(s) in RCA: 8] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/19/2022]
Abstract
We present the first implementation of multiconfiguration pair-density functional theory (MC-PDFT) ab initio molecular dynamics. MC-PDFT is a multireference electronic structure method that in many cases has a similar accuracy (or even better accuracy) the complete active space second-order perturbation theory (CASPT2) at a significantly lower computational cost. In this study, we introduced MC-PDFT analytical gradients into the SHARC molecular dynamics program for ab initio, nonadiabatic molecular dynamics simulations. We verify our implementation by examining the intersystem crossing dynamics of thioformaldehyde, and we observe excellent agreement with recent CASPT2 and experimental findings. Moreover, with MC-PDFT, we could perform dynamics simulations with the 12 electron in 10 orbitals active space that was computationally too expensive for direct dynamics with CASPT2.
Collapse
Affiliation(s)
- Paul B Calio
- Department of Chemistry, Chicago Center for Theoretical Chemistry, James Franck Institute, University of Chicago, Chicago, Illinois 60637-5418, United States
| | - Donald G Truhlar
- Department of Chemistry, Chemical Theory Center, and Minnesota Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431 United States
| | - Laura Gagliardi
- Department of Chemistry, Center for Theoretical Chemistry, Pritzker School of Molecular Engineering, Chicago James Franck Institute, University of Chicago, Chicago, Illinois 60637-5418, United States.,Argonne National Laboratory, Lemont, Illinois 60439 United States
| |
Collapse
|
9
|
Xie C, Guan Y, Yarkony DR, Guo H. Vibrational energy levels of the S0 and S1 states of formaldehyde using an accurate ab initio based global diabatic potential energy matrix. Mol Phys 2021. [DOI: 10.1080/00268976.2021.1918775] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/21/2022]
Affiliation(s)
- Changjian Xie
- Institute of Modern Physics, Northwest University, Xi’an, Shaanxi, People’s Republic of China
| | - Yafu Guan
- Department of Chemistry, Johns Hopkins University, Baltimore, MD, USA
| | - David R. Yarkony
- Department of Chemistry, Johns Hopkins University, Baltimore, MD, USA
| | - Hua Guo
- Department of Chemistry and Chemical Biology, University of New Mexico, Albuquerque, NM, USA
| |
Collapse
|
10
|
Kiss FL, Corbet BP, Simeth NA, Feringa BL, Crespi S. Predicting the substituent effects in the optical and electrochemical properties of N,N'-substituted isoindigos. Photochem Photobiol Sci 2021; 20:927-938. [PMID: 34227039 PMCID: PMC8550769 DOI: 10.1007/s43630-021-00071-5] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/29/2021] [Accepted: 06/24/2021] [Indexed: 11/30/2022]
Abstract
Isoindigo, the structural isomer of the well-known dye indigo, has seen a major revival recently because of the increasing interest of its use as a potential drug core structure and for the development of organic photovoltaic materials. Highly beneficial for diverse applications are its facile synthesis, straightforward functionalisation and the broad absorption band in the visible range. Moreover, its intrinsic electron deficiency renders isoindigo a promising acceptor structure in bulk heterojunction architectures. Here we present new insights into the substituent effects of N-functionalised isoindigos, developing a reliable and fast in silico screening approach of a library of compounds. Using experimental UV-Vis and electrochemical data increased the accuracy of the TD-DFT method employed. This procedure allowed us to accurately predict the optical and electrochemical properties of N-functionalised isoindigos and the elucidation of the relationship between substituent effects and electronic properties.
Collapse
Affiliation(s)
- Ferdinand L Kiss
- Faculty for Science and Engineering, Stratingh Institute for Chemistry, University of Groningen, Nijenborgh 4, 9747 AG, Groningen, The Netherlands
- Department Chemie, Ludwig-Maximilians-Universität München, 81377, Munich, Germany
| | - Brian P Corbet
- Faculty for Science and Engineering, Stratingh Institute for Chemistry, University of Groningen, Nijenborgh 4, 9747 AG, Groningen, The Netherlands
| | - Nadja A Simeth
- Faculty for Science and Engineering, Stratingh Institute for Chemistry, University of Groningen, Nijenborgh 4, 9747 AG, Groningen, The Netherlands
| | - Ben L Feringa
- Faculty for Science and Engineering, Stratingh Institute for Chemistry, University of Groningen, Nijenborgh 4, 9747 AG, Groningen, The Netherlands.
| | - Stefano Crespi
- Faculty for Science and Engineering, Stratingh Institute for Chemistry, University of Groningen, Nijenborgh 4, 9747 AG, Groningen, The Netherlands.
| |
Collapse
|
11
|
Marsili E, Prlj A, Curchod BFE. Caveat when using ADC(2) for studying the photochemistry of carbonyl-containing molecules. Phys Chem Chem Phys 2021; 23:12945-12949. [PMID: 34085679 PMCID: PMC8207513 DOI: 10.1039/d1cp02185k] [Citation(s) in RCA: 13] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/19/2021] [Accepted: 05/28/2021] [Indexed: 12/21/2022]
Abstract
Several electronic-structure methods are available to study the photochemistry and photophysics of organic molecules. Among them, ADC(2) stands as a sweet spot between computational efficiency and accuracy. As a result, ADC(2) has recently seen its number of applications booming, in particular to unravel the deactivation pathways and photodynamics of organic molecules. Despite this growing success, we demonstrate here that care has to be taken when studying the nonradiative pathways of carbonyl-containing molecules, as ADC(2) appears to suffer from a systematic flaw.
Collapse
Affiliation(s)
| | - Antonio Prlj
- Department of Chemistry, Durham University, Durham DH1 3LE, UK.
| | | |
Collapse
|
12
|
Abstract
In this article, we review nonadiabatic molecular dynamics (NAMD) methods for modeling spin-crossover transitions. First, we discuss different representations of electronic states employed in the grid-based and direct NAMD simulations. The nature of interstate couplings in different representations is highlighted, with the main focus on nonadiabatic and spin-orbit couplings. Second, we describe three NAMD methods that have been used to simulate spin-crossover dynamics, including trajectory surface hopping, ab initio multiple spawning, and multiconfiguration time-dependent Hartree. Some aspects of employing different electronic structure methods to obtain information about potential energy surfaces and interstate couplings for NAMD simulations are also discussed. Third, representative applications of NAMD to spin crossovers in molecular systems of different sizes and complexities are highlighted. Finally, we pose several fundamental questions related to spin-dependent processes. These questions should be possible to address with future methodological developments in NAMD.
Collapse
Affiliation(s)
- Saikat Mukherjee
- Institut de Chimie Radicalaire, CNRS 7273, Aix-Marseille University, 13013 Marseille, France;
| | - Dmitry A Fedorov
- Oak Ridge Associated Universities, Oak Ridge, Tennessee 37830, USA;
| | - Sergey A Varganov
- Department of Chemistry, University of Nevada, Reno, Nevada 89557-0216, USA;
| |
Collapse
|
13
|
Zhang L, Shu Y, Sun S, Truhlar DG. Direct coherent switching with decay of mixing for intersystem crossing dynamics of thioformaldehyde: The effect of decoherence. J Chem Phys 2021; 154:094310. [PMID: 33685154 DOI: 10.1063/5.0037878] [Citation(s) in RCA: 15] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/14/2023] Open
Abstract
We evaluate the effect of electronic decoherence on intersystem crossing in the photodynamics of thioformaldehyde. First, we show that the state-averaged complete-active-space self-consistent field electronic structure calculations with a properly chosen active space of 12 active electrons in 10 active orbitals can predict the potential energy surfaces and the singlet-triplet spin-orbit couplings quite well for CH2S, and we use this method for direct dynamics by coherent switching with decay of mixing (CSDM). We obtain similar dynamical results with CSDM or by adding energy-based decoherence to trajectory surface hopping, with the population of triplet states tending to a small steady-state value over 500 fs. Without decoherence, the state populations calculated by the conventional trajectory surface hopping method or the semiclassical Ehrenfest method gradually increase. This difference shows that decoherence changes the nature of the results not just quantitatively but qualitatively.
Collapse
Affiliation(s)
- Linyao Zhang
- School of Energy Science and Engineering, Harbin Institute of Technology, Harbin 150001, People's Republic of China
| | - Yinan Shu
- Department of Chemistry and Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, USA
| | - 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, USA
| |
Collapse
|
14
|
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
|
15
|
Guan Y, Xie C, Guo H, Yarkony DR. Neural Network Based Quasi-diabatic Representation for S0 and S1 States of Formaldehyde. J Phys Chem A 2020; 124:10132-10142. [DOI: 10.1021/acs.jpca.0c08948] [Citation(s) in RCA: 12] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Yafu Guan
- Department of Chemistry, Johns Hopkins University, Baltimore, Maryland 21218, United States
| | - Changjian Xie
- Institute of Modern Physics, Northwest University, Xi’an, Shaanxi 710069, People’s Republic of China
| | - Hua Guo
- Department of Chemistry and Chemical Biology, University of New Mexico, Albuquerque, New Mexico 87131, United States
| | - David R. Yarkony
- Department of Chemistry, Johns Hopkins University, Baltimore, Maryland 21218, United States
| |
Collapse
|
16
|
Fumanal M, Corminboeuf C. Direct, Mediated, and Delayed Intramolecular Singlet Fission Mechanism in Donor-Acceptor Copolymers. J Phys Chem Lett 2020; 11:9788-9794. [PMID: 33147966 DOI: 10.1021/acs.jpclett.0c03076] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/11/2023]
Abstract
Donor-acceptor (D-A) extended copolymers have shown great potential to be exploited for intramolecular singlet fission (iSF) because of their modular tunability and intrinsic ability to incorporate low-lying charge-transfer (CT) and a triplet-pair (1TT) states. While the SF mechanism has been widely debated in homo- and heterodimers, little is known about the singlet splitting process in A-D-A copolymer trimers. Unlike traditional two-site SF, the process of iSF in D-A copolymers involves three molecular units consisting of two A's and one D following an A-D-A polymeric chain. This scenario is, therefore, different from that of the homodimer analogues in terms of which states (if any) may drive the SF process. In this work, we identify how singlet splitting occurs in prototypical iSF D-A copolymer poly(benzodithiophene-alt-thiophene-1,1-dioxide) by means of wave packet propagations on the basis of the linear vibronic coupling model Hamiltonian. Our results reveal that three different mechanisms drive the S1 → 1TT population transfer via antisymmetric and symmetric vibrational motion, including two favorable mechanisms of direct and mediated interactions, as well as a parasitic decay pathway that potentially delays the process. Remarkably, we uncover the interplay between an upper state of marked multiexcitonic character and a low-lying CT state in balancing the splitting efficiency, which anticipates their major role in defining future guidelines for the molecular design of D-A copolymers for iSF.
Collapse
Affiliation(s)
- Maria Fumanal
- Laboratory for Computational Molecular Design, Institute of Chemical Sciences and Engineering, École Polytechnique Fédérale de Lausanne (EPFL), CH-1015 Lausanne, Switzerland
| | - Clémence Corminboeuf
- Laboratory for Computational Molecular Design, Institute of Chemical Sciences and Engineering, École Polytechnique Fédérale de Lausanne (EPFL), CH-1015 Lausanne, Switzerland
| |
Collapse
|
17
|
Mai S, González L. Molecular Photochemistry: Recent Developments in Theory. Angew Chem Int Ed Engl 2020; 59:16832-16846. [PMID: 32052547 PMCID: PMC7540682 DOI: 10.1002/anie.201916381] [Citation(s) in RCA: 77] [Impact Index Per Article: 19.3] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/20/2019] [Revised: 02/12/2020] [Indexed: 12/16/2022]
Abstract
Photochemistry is a fascinating branch of chemistry that is concerned with molecules and light. However, the importance of simulating light-induced processes is reflected also in fields as diverse as biology, material science, and medicine. This Minireview highlights recent progress achieved in theoretical chemistry to calculate electronically excited states of molecules and simulate their photoinduced dynamics, with the aim of reaching experimental accuracy. We focus on emergent methods and give selected examples that illustrate the progress in recent years towards predicting complex electronic structures with strong correlation, calculations on large molecules, describing multichromophoric systems, and simulating non-adiabatic molecular dynamics over long time scales, for molecules in the gas phase or in complex biological environments.
Collapse
Affiliation(s)
- Sebastian Mai
- Photonics InstituteVienna University of TechnologyGusshausstrasse 27–291040ViennaAustria
| | - Leticia González
- Institute of Theoretical ChemistryFaculty of ChemistryUniversity of ViennaWähringer Strasse 171090ViennaAustria
| |
Collapse
|
18
|
Khalili K, Inhester L, Arnold C, Gertsen AS, Andreasen JW, Santra R. Simulation of time-resolved x-ray absorption spectroscopy of ultrafast dynamics in particle-hole-excited 4-(2-thienyl)-2,1,3-benzothiadiazole. STRUCTURAL DYNAMICS (MELVILLE, N.Y.) 2020; 7:044101. [PMID: 32665964 PMCID: PMC7340508 DOI: 10.1063/4.0000016] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 03/27/2020] [Accepted: 06/15/2020] [Indexed: 06/01/2023]
Abstract
To date, alternating co-polymers based on electron-rich and electron-poor units are the most attractive materials to control functionality of organic semiconductor layers in which ultrafast excited-state processes play a key role. We present a computational study of the photoinduced excited-state dynamics of the 4-(2-thienyl)-2,1,3-benzothiadiazole (BT-1T) molecule, which is a common building block in the backbone of π-conjugated polymers used for organic electronics. In contrast to homo-polymer materials, such as oligothiophene, BT-1T has two non-identical units, namely, thiophene and benzothiadiazole, making it attractive for intramolecular charge transfer studies. To gain a thorough understanding of the coupling of excited-state dynamics with nuclear motion, we consider a scenario based on femtosecond time-resolved x-ray absorption spectroscopy using an x-ray free-electron laser in combination with a synchronized ultraviolet femtosecond laser. Using Tully's fewest switches surface hopping approach in combination with excited-state calculations at the level of configuration interaction singles, we calculate the gas-phase x-ray absorption spectrum at the carbon and nitrogen K edges as a function of time after excitation to the lowest electronically excited state. The results of our time-resolved calculations exhibit the charge transfer driven by non-Born-Oppenheimer physics from the benzothiadiazole to thiophene units during relaxation to the ground state. Furthermore, our ab initio molecular dynamics simulations indicate that the excited-state relaxation processes involve bond elongation in the benzothiadiazole unit as well as thiophene ring puckering at a time scale of 100 fs. We show that these dynamical trends can be identified from the time-dependent x-ray absorption spectrum.
Collapse
Affiliation(s)
- Khadijeh Khalili
- Department of Energy Conversion and Storage, Technical University of Denmark, Fysikvej 310, 2800 Kgs. Lyngby, Denmark
| | | | | | - Anders S. Gertsen
- Department of Energy Conversion and Storage, Technical University of Denmark, Fysikvej 310, 2800 Kgs. Lyngby, Denmark
| | - Jens Wenzel Andreasen
- Department of Energy Conversion and Storage, Technical University of Denmark, Fysikvej 310, 2800 Kgs. Lyngby, Denmark
| | | |
Collapse
|
19
|
Mai S, González L. Molekulare Photochemie: Moderne Entwicklungen in der theoretischen Chemie. Angew Chem Int Ed Engl 2020. [DOI: 10.1002/ange.201916381] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/17/2022]
Affiliation(s)
- Sebastian Mai
- Institut für Photonik Technische Universität Wien Gußhausstraße 27–29 1040 Wien Österreich
| | - Leticia González
- Institut für theoretische Chemie Fakultät für Chemie Universität Wien Währinger Straße 17 1090 Wien Österreich
| |
Collapse
|
20
|
Westermayr J, Gastegger M, Marquetand P. Combining SchNet and SHARC: The SchNarc Machine Learning Approach for Excited-State Dynamics. J Phys Chem Lett 2020; 11:3828-3834. [PMID: 32311258 PMCID: PMC7246974 DOI: 10.1021/acs.jpclett.0c00527] [Citation(s) in RCA: 86] [Impact Index Per Article: 21.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/17/2020] [Accepted: 04/20/2020] [Indexed: 05/26/2023]
Abstract
In recent years, deep learning has become a part of our everyday life and is revolutionizing quantum chemistry as well. In this work, we show how deep learning can be used to advance the research field of photochemistry by learning all important properties-multiple energies, forces, and different couplings-for photodynamics simulations. We simplify such simulations substantially by (i) a phase-free training skipping costly preprocessing of raw quantum chemistry data; (ii) rotationally covariant nonadiabatic couplings, which can either be trained or (iii) alternatively be approximated from only ML potentials, their gradients, and Hessians; and (iv) incorporating spin-orbit couplings. As the deep-learning method, we employ SchNet with its automatically determined representation of molecular structures and extend it for multiple electronic states. In combination with the molecular dynamics program SHARC, our approach termed SchNarc is tested on two polyatomic molecules and paves the way toward efficient photodynamics simulations of complex systems.
Collapse
Affiliation(s)
- Julia Westermayr
- Institute
of Theoretical Chemistry, Faculty of Chemistry, University of Vienna, Währinger Str. 17, 1090 Vienna, Austria
| | - Michael Gastegger
- Machine
Learning Group, Technical University of
Berlin, 10587 Berlin, Germany
| | - Philipp Marquetand
- Institute
of Theoretical Chemistry, Faculty of Chemistry, University of Vienna, Währinger Str. 17, 1090 Vienna, Austria
- Vienna
Research Platform on Accelerating Photoreaction Discovery, University of Vienna, Währinger Str. 17, 1090 Vienna, Austria
- Data
Science @ Uni Vienna, University of Vienna, Währinger Str. 29, 1090 Vienna, Austria
| |
Collapse
|
21
|
Park JW, Al-Saadon R, MacLeod MK, Shiozaki T, Vlaisavljevich B. Multireference Electron Correlation Methods: Journeys along Potential Energy Surfaces. Chem Rev 2020; 120:5878-5909. [PMID: 32239929 DOI: 10.1021/acs.chemrev.9b00496] [Citation(s) in RCA: 54] [Impact Index Per Article: 13.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/16/2022]
Abstract
Multireference electron correlation methods describe static and dynamical electron correlation in a balanced way and, therefore, can yield accurate and predictive results even when single-reference methods or multiconfigurational self-consistent field theory fails. One of their most prominent applications in quantum chemistry is the exploration of potential energy surfaces. This includes the optimization of molecular geometries, such as equilibrium geometries and conical intersections and on-the-fly photodynamics simulations, both of which depend heavily on the ability of the method to properly explore the potential energy surface. Because such applications require nuclear gradients and derivative couplings, the availability of analytical nuclear gradients greatly enhances the scope of quantum chemical methods. This review focuses on the developments and advances made in the past two decades. A detailed account of the analytical nuclear gradient and derivative coupling theories is presented. Emphasis is given to the software infrastructure that allows one to make use of these methods. Notable applications of multireference electron correlation methods to chemistry, including geometry optimizations and on-the-fly dynamics, are summarized at the end followed by a discussion of future prospects.
Collapse
Affiliation(s)
- Jae Woo Park
- Department of Chemistry, Chungbuk National University, Chungdae-ro 1, Cheongju 28644, Korea
| | - Rachael Al-Saadon
- Department of Chemistry, Northwestern University, 2145 Sheridan Road, Evanston, Illinois 60208, United States
| | - Matthew K MacLeod
- Workday, 4900 Pearl Circle East, Suite 100, Boulder, Colorado 80301, United States
| | - Toru Shiozaki
- Department of Chemistry, Northwestern University, 2145 Sheridan Road, Evanston, Illinois 60208, United States.,Quantum Simulation Technologies, Inc., 625 Massachusetts Avenue, Cambridge, Massachusetts 02139, United States
| | - Bess Vlaisavljevich
- Department of Chemistry, University of South Dakota, 414 East Clark Street, Vermillion, South Dakota 57069, United States
| |
Collapse
|
22
|
Guan Y, Yarkony DR. Accurate Neural Network Representation of the Ab Initio Determined Spin-Orbit Interaction in the Diabatic Representation Including the Effects of Conical Intersections. J Phys Chem Lett 2020; 11:1848-1858. [PMID: 32062966 DOI: 10.1021/acs.jpclett.0c00074] [Citation(s) in RCA: 21] [Impact Index Per Article: 5.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/26/2023]
Abstract
A method for fitting ab initio determined spin-orbit coupling interactions, in the Breit-Pauli approximation, based on quasidiabatic representations using neural network fits is reported. The algorithm generalizes our recently reported neural network approach for representing the dipole interaction. The S0, S1, and T1 states of formaldehyde are used as an example. First, the two singlet states S0 and S1 are diabatized with a modified Boys Localization diabatization method. Second, the spin-orbit coupling between singlet and triplet states is transformed to the diabatic representation. This removes the discontinuities in the adiabatic representation. The diabatized spin-orbit couplings are then fit with smooth neural network functions. The analytic representation of spin-orbit coupling interactions in a diabatic basis by neural networks will make accurate full-dimensional quantum dynamical treatment of both internal conversion and intersystem crossing possible, which will help us to gain better understanding of both processes.
Collapse
Affiliation(s)
- Yafu Guan
- Department of Chemistry, Johns Hopkins University, Baltimore, Maryland 21218, United States
| | - David R Yarkony
- Department of Chemistry, Johns Hopkins University, Baltimore, Maryland 21218, United States
| |
Collapse
|
23
|
Smith B, Akimov AV. Modeling nonadiabatic dynamics in condensed matter materials: some recent advances and applications. JOURNAL OF PHYSICS. CONDENSED MATTER : AN INSTITUTE OF PHYSICS JOURNAL 2020; 32:073001. [PMID: 31661681 DOI: 10.1088/1361-648x/ab5246] [Citation(s) in RCA: 33] [Impact Index Per Article: 8.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/10/2023]
Abstract
This review focuses on recent developments in the field of nonadiabatic molecular dynamics (NA-MD), with particular attention given to condensed-matter systems. NA-MD simulations for small molecular systems can be performed using high-level electronic structure (ES) calculations, methods accounting for the quantization of nuclear motion, and using fewer approximations in the dynamical methodology itself. Modeling condensed-matter systems imposes many limitations on various aspects of NA-MD computations, requiring approximations at various levels of theory-from the ES, to the ways in which the coupling of electrons and nuclei are accounted for. Nonetheless, the approximate treatment of NA-MD in condensed-phase materials has gained a spin lately in many applied studies. A number of advancements of the methodology and computational tools have been undertaken, including general-purpose methods, as well as those tailored to nanoscale and condensed matter systems. This review summarizes such methodological and software developments, puts them into the broader context of existing approaches, and highlights some of the challenges that remain to be solved.
Collapse
Affiliation(s)
- Brendan Smith
- Department of Chemistry, University at Buffalo, The State University of New York, Buffalo, New York 14260-3000, United States of America
| | | |
Collapse
|
24
|
Nelson TR, White AJ, Bjorgaard JA, Sifain AE, Zhang Y, Nebgen B, Fernandez-Alberti S, Mozyrsky D, Roitberg AE, Tretiak S. Non-adiabatic Excited-State Molecular Dynamics: Theory and Applications for Modeling Photophysics in Extended Molecular Materials. Chem Rev 2020; 120:2215-2287. [PMID: 32040312 DOI: 10.1021/acs.chemrev.9b00447] [Citation(s) in RCA: 198] [Impact Index Per Article: 49.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/22/2022]
Abstract
Optically active molecular materials, such as organic conjugated polymers and biological systems, are characterized by strong coupling between electronic and vibrational degrees of freedom. Typically, simulations must go beyond the Born-Oppenheimer approximation to account for non-adiabatic coupling between excited states. Indeed, non-adiabatic dynamics is commonly associated with exciton dynamics and photophysics involving charge and energy transfer, as well as exciton dissociation and charge recombination. Understanding the photoinduced dynamics in such materials is vital to providing an accurate description of exciton formation, evolution, and decay. This interdisciplinary field has matured significantly over the past decades. Formulation of new theoretical frameworks, development of more efficient and accurate computational algorithms, and evolution of high-performance computer hardware has extended these simulations to very large molecular systems with hundreds of atoms, including numerous studies of organic semiconductors and biomolecules. In this Review, we will describe recent theoretical advances including treatment of electronic decoherence in surface-hopping methods, the role of solvent effects, trivial unavoided crossings, analysis of data based on transition densities, and efficient computational implementations of these numerical methods. We also emphasize newly developed semiclassical approaches, based on the Gaussian approximation, which retain phase and width information to account for significant decoherence and interference effects while maintaining the high efficiency of surface-hopping approaches. The above developments have been employed to successfully describe photophysics in a variety of molecular materials.
Collapse
Affiliation(s)
- Tammie R Nelson
- Theoretical Division , Los Alamos National Laboratory , Los Alamos , New Mexico 87545 , United States
| | - Alexander J White
- Theoretical Division , Los Alamos National Laboratory , Los Alamos , New Mexico 87545 , United States
| | - Josiah A Bjorgaard
- Theoretical Division , Los Alamos National Laboratory , Los Alamos , New Mexico 87545 , United States
| | - Andrew E Sifain
- Theoretical Division , Los Alamos National Laboratory , Los Alamos , New Mexico 87545 , United States.,U.S. Army Research Laboratory , Aberdeen Proving Ground , Maryland 21005 , United States
| | - Yu Zhang
- Theoretical Division , Los Alamos National Laboratory , Los Alamos , New Mexico 87545 , United States
| | - Benjamin Nebgen
- Theoretical Division , Los Alamos National Laboratory , Los Alamos , New Mexico 87545 , United States
| | | | - Dmitry Mozyrsky
- Theoretical Division , Los Alamos National Laboratory , Los Alamos , New Mexico 87545 , United States
| | - Adrian E Roitberg
- Department of Chemistry , University of Florida , Gainesville , Florida 32611 , United States
| | - Sergei Tretiak
- Theoretical Division , Los Alamos National Laboratory , Los Alamos , New Mexico 87545 , United States
| |
Collapse
|
25
|
Bellonzi N, Alguire E, Fatehi S, Shao Y, Subotnik JE. TD-DFT spin-adiabats with analytic nonadiabatic derivative couplings. J Chem Phys 2020; 152:044112. [PMID: 32007078 PMCID: PMC7043850 DOI: 10.1063/1.5126440] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/10/2019] [Accepted: 12/22/2019] [Indexed: 11/14/2022] Open
Abstract
We present an algorithm for efficient calculation of analytic nonadiabatic derivative couplings between spin-adiabatic, time-dependent density functional theory states within the Tamm-Dancoff approximation. Our derivation is based on the direct differentiation of the Kohn-Sham pseudowavefunction using the framework of Ou et al. Our implementation is limited to the case of a system with an even number of electrons in a closed shell ground state, and we validate our algorithm against finite difference at an S1/T2 crossing of benzaldehyde. Through the introduction of a magnetic field spin-coupling operator, we break time-reversal symmetry to generate complex valued nonadiabatic derivative couplings. Although the nonadiabatic derivative couplings are complex valued, we find that a phase rotation can generate an almost entirely real-valued derivative coupling vector for the case of benzaldehyde.
Collapse
Affiliation(s)
- Nicole Bellonzi
- Department of Chemistry, University of Pennsylvania, Philadelphia, Pennsylvania 19104, USA
| | | | - Shervin Fatehi
- Department of Chemistry, The University of Texas Rio Grande Valley, Edinburg, Texas 78539, USA
| | - Yihan Shao
- Department of Chemistry and Biochemistry, University of Oklahoma, Norman, Oklahoma 73019, USA
| | - Joseph E. Subotnik
- Department of Chemistry, University of Pennsylvania, Philadelphia, Pennsylvania 19104, USA
| |
Collapse
|
26
|
Plasser F, Mai S, Fumanal M, Gindensperger E, Daniel C, González L. Strong Influence of Decoherence Corrections and Momentum Rescaling in Surface Hopping Dynamics of Transition Metal Complexes. J Chem Theory Comput 2019; 15:5031-5045. [DOI: 10.1021/acs.jctc.9b00525] [Citation(s) in RCA: 33] [Impact Index Per Article: 6.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/16/2023]
Affiliation(s)
- Felix Plasser
- Institute for Theoretical Chemistry, Faculty of Chemistry, University of Vienna, Währingerstr. 17, 1090 Vienna, Austria
- Department of Chemistry, Loughborough University, Loughborough, LE11 3TU, U.K
| | - Sebastian Mai
- Institute for Theoretical Chemistry, Faculty of Chemistry, University of Vienna, Währingerstr. 17, 1090 Vienna, Austria
| | - Maria Fumanal
- Laboratoire de Chimie Quantique, Institut de Chimie Strasbourg, UMR7177 CNRS/Université de Strasbourg 4 Rue Blaise Pascal BP296/R8, F-67008 Strasbourg, France
| | - Etienne Gindensperger
- Laboratoire de Chimie Quantique, Institut de Chimie Strasbourg, UMR7177 CNRS/Université de Strasbourg 4 Rue Blaise Pascal BP296/R8, F-67008 Strasbourg, France
| | - Chantal Daniel
- Laboratoire de Chimie Quantique, Institut de Chimie Strasbourg, UMR7177 CNRS/Université de Strasbourg 4 Rue Blaise Pascal BP296/R8, F-67008 Strasbourg, France
| | - Leticia González
- Institute for Theoretical Chemistry, Faculty of Chemistry, University of Vienna, Währingerstr. 17, 1090 Vienna, Austria
| |
Collapse
|