1
|
Curchod BFE, Orr-Ewing AJ. Perspective on Theoretical and Experimental Advances in Atmospheric Photochemistry. J Phys Chem A 2024. [PMID: 39021090 DOI: 10.1021/acs.jpca.4c03481] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [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
|
2
|
Morreale D, Persico M. Topology of Conical Intersection Seams and the Geometric Phase. J Phys Chem A 2024; 128:1707-1714. [PMID: 38408203 DOI: 10.1021/acs.jpca.3c07042] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/28/2024]
Abstract
In this paper, we demonstrate two topological properties of crossing seams, that is, the sets of points in the N-dimensional space of nuclear coordinates where two electronic eigenstates are degenerate. We shall examine the typical case of states of the same spin with accidental degeneracies, whereby the crossing seam is of dimension N - 2. The first property we demonstrate is that a crossing seam has no boundary, therefore, it must either extend asymptotically to infinite values of one or more coordinates or wrap on itself. The second property is that two (or more) crossing seams can intersect each other but in such a way that neither of them ends at the intersection. When N = 3, the crossing seam is a line in a 3D space; this is so in triatomic molecules but also in reduced dimensionality treatments of larger polyatomics. The above-mentioned rules then mean that the crossing seam is a line of infinite length or a closed loop and can split into three branches but not in two. The example of the first two excited 1A' states of H2Cl+ illustrates these rules and shows their usefulness for computational search and characterization of crossing seams.
Collapse
Affiliation(s)
- Daniel Morreale
- Department of Chemistry and Industrial Chemistry, University of Pisa, v. G. Moruzzi 13, Pisa 56126, Italy
| | - Maurizio Persico
- Department of Chemistry and Industrial Chemistry, University of Pisa, v. G. Moruzzi 13, Pisa 56126, Italy
| |
Collapse
|
3
|
Rivero Santamaría A, Piris M. Time evolution of natural orbitals in ab initio molecular dynamics. J Chem Phys 2024; 160:071102. [PMID: 38364005 DOI: 10.1063/5.0188491] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/21/2023] [Accepted: 01/28/2024] [Indexed: 02/18/2024] Open
Abstract
This work combines for the first time ab initio molecular dynamics (AIMD) within the Born-Oppenheimer approximation with a global natural orbital functional (GNOF), an approximate functional of the one-particle reduced density matrix. The most prominent feature of GNOF-AIMD is its ability to display the real-time evolution of natural orbitals, providing detailed information on the time-dependent electronic structure of complex systems and processes, including reactive collisions. The quartet ground-state reaction N(4S) + H2(1Σ) → NH(3Σ) + H(2S) is taken as a validation test. Collision energy influences on integral cross sections for different initial rovibrational states of H2 and rotational-state distributions of the NH product are discussed, showing a good agreement with previous high-quality theoretical results.
Collapse
Affiliation(s)
| | - Mario Piris
- Donostia International Physics Center (DIPC), 20018 Donostia, Euskadi, Spain; Euskal Herriko Unibertsitatea (UPV/EHU), PK 1072, 20080 Donostia, Euskadi, Spain; and Basque Foundation for Science (IKERBASQUE), 48009 Bilbao, Euskadi, Spain
| |
Collapse
|
4
|
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
|
5
|
Janoš J, Vinklárek IS, Rakovský J, Mukhopadhyay DP, Curchod BFE, Fárník M, Slavíček P. On the Wavelength-Dependent Photochemistry of the Atmospheric Molecule CF 3COCl. ACS EARTH & SPACE CHEMISTRY 2023; 7:2275-2286. [PMID: 38026808 PMCID: PMC10658617 DOI: 10.1021/acsearthspacechem.3c00196] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 07/05/2023] [Revised: 09/29/2023] [Accepted: 09/29/2023] [Indexed: 12/01/2023]
Abstract
The wavelength control of photochemistry usually results from ultrafast dynamics following the excitation of different electronic states. Here, we investigate the CF3COCl molecule, exhibiting wavelength-dependent photochemistry both via (i) depositing increasing internal energy into a single state and (ii) populating different electronic states. We reveal the mechanism behind the photon-energy dependence by combining nonadiabatic ab initio molecular dynamics techniques with the velocity map imaging experiment. We describe a consecutive mechanism of photodissociation where an immediate release of Cl taking place in an excited electronic state is followed by a slower ground-state dissociation of the CO fragment. The CO release is subject to an activation barrier and is controlled by excess internal energy via the excitation wavelength. Therefore, a selective release of CO along with Cl can be achieved. The mechanism is fully supported by both the measured kinetic energy distributions and anisotropies of the angular distributions. Interestingly, the kinetic energy of the released Cl atom is sensitively modified by accounting for spin-orbit coupling. Given the atmospheric importance of CF3COCl, we discuss the consequences of our findings for atmospheric photochemistry.
Collapse
Affiliation(s)
- Jiří Janoš
- Department
of Physical Chemistry, University of Chemistry
and Technology, Technická 5, Prague 6 166 28, Czech Republic
| | - Ivo S. Vinklárek
- Department
of Dynamics of Molecules and Clusters, J. Heyrovský Institute
of Physical Chemistry, v.v.i., The Czech
Academy of Sciences, Dolejškova 2155/3, 182 23 Prague, Czech Republic
| | - Jozef Rakovský
- Department
of Dynamics of Molecules and Clusters, J. Heyrovský Institute
of Physical Chemistry, v.v.i., The Czech
Academy of Sciences, Dolejškova 2155/3, 182 23 Prague, Czech Republic
| | - Deb Pratim Mukhopadhyay
- Department
of Dynamics of Molecules and Clusters, J. Heyrovský Institute
of Physical Chemistry, v.v.i., The Czech
Academy of Sciences, Dolejškova 2155/3, 182 23 Prague, Czech Republic
| | - Basile F. E. Curchod
- Centre
for Computational Chemistry, School of Chemistry, University of Bristol, Bristol BS8 1TS, U.K.
| | - Michal Fárník
- Department
of Dynamics of Molecules and Clusters, J. Heyrovský Institute
of Physical Chemistry, v.v.i., The Czech
Academy of Sciences, Dolejškova 2155/3, 182 23 Prague, Czech Republic
| | - Petr Slavíček
- Department
of Physical Chemistry, University of Chemistry
and Technology, Technická 5, Prague 6 166 28, Czech Republic
| |
Collapse
|
6
|
Pracht P, Bannwarth C. Finding Excited-State Minimum Energy Crossing Points on a Budget: Non-Self-Consistent Tight-Binding Methods. J Phys Chem Lett 2023; 14:4440-4448. [PMID: 37144783 DOI: 10.1021/acs.jpclett.3c00494] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/06/2023]
Abstract
The automated exploration and identification of minimum energy conical intersections (MECIs) is a valuable computational strategy for the study of photochemical processes. Due to the immense computational effort involved in calculating non-adiabatic derivative coupling vectors, simplifications have been introduced focusing instead on minimum energy crossing points (MECPs), where promising attempts were made with semiempirical quantum mechanical methods. A simplified treatment for describing crossing points between almost arbitrary diabatic states based on a non-self-consistent extended tight-binding method, GFN0-xTB, is presented. Involving only a single diagonalization of the Hamiltonian, the method can provide energies and gradients for multiple electronic states, which can be used in a derivative coupling-vector-free scheme to calculate MECPs. By comparison with high-lying MECIs of benchmark systems, it is demonstrated that the identified geometries are good starting points for further MECI refinement with ab initio methods.
Collapse
Affiliation(s)
- Philipp Pracht
- Yusuf Hamied Department of Chemistry, University of Cambridge, Lensfield Road, Cambridge CB2 1EW, United Kingdom
| | - Christoph Bannwarth
- Institute for Physical Chemistry, RWTH Aachen University, Melatener Str. 20, 52074 Aachen, Germany
| |
Collapse
|
7
|
Cruzeiro VWD, Wang Y, Pieri E, Hohenstein EG, Martínez TJ. TeraChem protocol buffers (TCPB): Accelerating QM and QM/MM simulations with a client-server model. J Chem Phys 2023; 158:044801. [PMID: 36725506 DOI: 10.1063/5.0130886] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/06/2023] Open
Abstract
The routine use of electronic structures in many chemical simulation applications calls for efficient and easy ways to access electronic structure programs. We describe how the graphics processing unit (GPU) accelerated electronic structure program TeraChem can be set up as an electronic structure server, to be easily accessed by third-party client programs. We exploit Google's protocol buffer framework for data serialization and communication. The client interface, called TeraChem protocol buffers (TCPB), has been designed for ease of use and compatibility with multiple programming languages, such as C++, Fortran, and Python. To demonstrate the ease of coupling third-party programs with electronic structures using TCPB, we have incorporated the TCPB client into Amber for quantum mechanics/molecular mechanics (QM/MM) simulations. The TCPB interface saves time with GPU initialization and I/O operations, achieving a speedup of more than 2× compared to a prior file-based implementation for a QM region with ∼250 basis functions. We demonstrate the practical application of TCPB by computing the free energy profile of p-hydroxybenzylidene-2,3-dimethylimidazolinone (p-HBDI-)-a model chromophore in green fluorescent proteins-on the first excited singlet state using Hamiltonian replica exchange for enhanced sampling. All calculations in this work have been performed with the non-commercial freely-available version of TeraChem, which is sufficient for many QM region sizes in common use.
Collapse
Affiliation(s)
| | - Yuanheng Wang
- Department of Chemistry and The PULSE Institute, Stanford University, Stanford, California 94305, USA
| | - Elisa Pieri
- Department of Chemistry and The PULSE Institute, Stanford University, Stanford, California 94305, USA
| | - Edward G Hohenstein
- Department of Chemistry and The PULSE Institute, Stanford University, Stanford, California 94305, USA
| | - Todd J Martínez
- Department of Chemistry and The PULSE Institute, Stanford University, Stanford, California 94305, USA
| |
Collapse
|
8
|
Ovad T, Sapunar M, Sršeň Š, Slavíček P, Mašín Z, Jones NC, Hoffmann SV, Ranković M, Fedor J. Excitation and fragmentation of the dielectric gas C 4F 7N: Electrons vs photons. J Chem Phys 2023; 158:014303. [PMID: 36610949 DOI: 10.1063/5.0130216] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/13/2022] Open
Abstract
C4F7N is a promising candidate for the replacement of sulfur hexafluoride as an insulating medium, and it is important to understand the chemical changes initiated in the molecule by collision with free electrons, specifically the formation of neutral fragments. The first step of neutral fragmentation is electronic excitation, yet neither the absorption spectrum in the vacuum ultraviolet (VUV) region nor the electron energy loss spectrum have previously been reported. Here, we experimentally probed the excited states by VUV photoabsorption spectroscopy and electron energy loss spectroscopy (EELS). We found that the distribution of states populated upon electron impact with low-energy electrons is significantly different from that following photoabsorption. This difference was confirmed and interpreted with ab initio modeling of both VUV and EELS spectra. We propose here a new computational protocol for the simulation of EELS spectra combining the Born approximation with approximate forms of correlated wave functions, which allows us to calculate the (usually very expensive) scattering cross sections at a cost similar to the calculation of oscillator strengths. Finally, we perform semi-classical non-adiabatic dynamics simulations to investigate the possible neutral fragments of the molecule formed through electron-induced neutral dissociation. We show that the product distribution is highly non-statistical.
Collapse
Affiliation(s)
- Tomáš Ovad
- Department of Physical Chemistry, University of Chemistry and Technology, Technická 5, Prague 6, Czech Republic
| | - Marin Sapunar
- Department of Physical Chemistry, University of Chemistry and Technology, Technická 5, Prague 6, Czech Republic
| | - Štěpán Sršeň
- Department of Physical Chemistry, University of Chemistry and Technology, Technická 5, Prague 6, Czech Republic
| | - Petr Slavíček
- Department of Physical Chemistry, University of Chemistry and Technology, Technická 5, Prague 6, Czech Republic
| | - Zdeněk Mašín
- Faculty of Mathematics and Physics, Charles University, Institute of Theoretical Physics, V Holešovičkách 2, 18000 Prague, Czech Republic
| | - Nykola C Jones
- ISA, Department of Physics and Astronomy, Aarhus University, Ny Munkegade 120, 8000 Aarhus C, Denmark
| | - Søren Vrønning Hoffmann
- ISA, Department of Physics and Astronomy, Aarhus University, Ny Munkegade 120, 8000 Aarhus C, Denmark
| | - Miloš Ranković
- J. Heyrovský Institute of Physical Chemistry, Czech Academy of Sciences, Dolejškova 3, Prague 8, Czech Republic
| | - Juraj Fedor
- J. Heyrovský Institute of Physical Chemistry, Czech Academy of Sciences, Dolejškova 3, Prague 8, Czech Republic
| |
Collapse
|
9
|
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
|
10
|
Pracht P, Bannwarth C. Fast Screening of Minimum Energy Crossing Points with Semiempirical Tight-Binding Methods. J Chem Theory Comput 2022; 18:6370-6385. [PMID: 36121838 DOI: 10.1021/acs.jctc.2c00578] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
The investigation of photochemical processes is a highly active field in computational chemistry. One research direction is the automated exploration and identification of minimum energy conical intersection (MECI) geometries. However, due to the immense technical effort required to calculate nonadiabatic potential energy landscapes, the routine application of such computational protocols is severely limited. In this study, we will discuss the prospect of combining adiabatic potential energy surfaces from semiempirical quantum mechanical calculations with specialized confinement potential and metadynamics simulations to identify S0/T1 minimum energy crossing point (MECP) geometries. It is shown that MECPs calculated at the GFN2-xTB level can provide suitable approximations to high-level S0/S1ab initio conical intersection geometries at a fraction of the computational cost. Reference MECIs of benzene are studied to illustrate the basic concept. An example application of the presented protocol is demonstrated for a set of photoswitch molecules.
Collapse
Affiliation(s)
- Philipp Pracht
- Institute of Physical Chemistry, RWTH Aachen University, Melatener Str. 20, 52056Aachen, Germany
| | - Christoph Bannwarth
- Institute of Physical Chemistry, RWTH Aachen University, Melatener Str. 20, 52056Aachen, Germany
| |
Collapse
|
11
|
Wang Y, Seritan S, Lahana D, Ford JE, Valentini A, Hohenstein EG, Martínez TJ. InteraChem: Exploring Excited States in Virtual Reality with Ab Initio Interactive Molecular Dynamics. J Chem Theory Comput 2022; 18:3308-3317. [PMID: 35649124 DOI: 10.1021/acs.jctc.2c00005] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/14/2023]
Abstract
InteraChem is an ab initio interactive molecular dynamics (AI-IMD) visualizer that leverages recent advances in virtual reality hardware and software, as well as the graphical processing unit (GPU)-accelerated TeraChem electronic structure package, in order to render quantum chemistry in real time. We introduce the exploration of electronically excited states via AI-IMD using the floating occupation molecular orbital-complete active space configuration interaction method. The optimization tools in InteraChem enable identification of excited state minima as well as minimum energy conical intersections for further characterization of excited state chemistry in small- to medium-sized systems. We demonstrate that finite-temperature Hartree-Fock theory is an efficient method to perform ground state AI-IMD. InteraChem allows users to track electronic properties such as molecular orbitals and bond order in real time, resulting in an interactive visualization tool that aids in the interpretation of excited state chemistry data and makes quantum chemistry more accessible for both research and educational purposes.
Collapse
Affiliation(s)
- Yuanheng Wang
- Department of Chemistry and The PULSE Institute, Stanford University, Stanford, California 94305, United States.,SLAC National Accelerator Laboratory, 2575 Sand Hill Road, Menlo Park, California 94025, United States
| | - Stefan Seritan
- Department of Chemistry and The PULSE Institute, Stanford University, Stanford, California 94305, United States.,SLAC National Accelerator Laboratory, 2575 Sand Hill Road, Menlo Park, California 94025, United States
| | - Dean Lahana
- Department of Chemistry and The PULSE Institute, Stanford University, Stanford, California 94305, United States.,SLAC National Accelerator Laboratory, 2575 Sand Hill Road, Menlo Park, California 94025, United States
| | - Jason E Ford
- Department of Chemistry and The PULSE Institute, Stanford University, Stanford, California 94305, United States.,SLAC National Accelerator Laboratory, 2575 Sand Hill Road, Menlo Park, California 94025, United States
| | - Alessio Valentini
- Department of Chemistry and The PULSE Institute, Stanford University, Stanford, California 94305, United States.,SLAC National Accelerator Laboratory, 2575 Sand Hill Road, Menlo Park, California 94025, United States
| | - Edward G Hohenstein
- Department of Chemistry and The PULSE Institute, Stanford University, Stanford, California 94305, United States.,SLAC National Accelerator Laboratory, 2575 Sand Hill Road, Menlo Park, California 94025, United States
| | - Todd J Martínez
- Department of Chemistry and The PULSE Institute, Stanford University, Stanford, California 94305, United States.,SLAC National Accelerator Laboratory, 2575 Sand Hill Road, Menlo Park, California 94025, United States
| |
Collapse
|
12
|
Raucci U, Weir H, Bannwarth C, Sanchez DM, Martínez TJ. Chiral photochemistry of achiral molecules. Nat Commun 2022; 13:2091. [PMID: 35440559 PMCID: PMC9019051 DOI: 10.1038/s41467-022-29662-1] [Citation(s) in RCA: 7] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/04/2021] [Accepted: 03/25/2022] [Indexed: 11/23/2022] Open
Abstract
Chirality is a molecular property governed by the topography of the potential energy surface (PES). Thermally achiral molecules interconvert rapidly when the interconversion barrier between the two enantiomers is comparable to or lower than the thermal energy, in contrast to thermally stable chiral configurations. In principle, a change in the PES topography on the excited electronic state may diminish interconversion, leading to electronically prochiral molecules that can be converted from achiral to chiral by electronic excitation. Here we report that this is the case for two prototypical examples – cis-stilbene and cis-stiff stilbene. Both systems exhibit unidirectional photoisomerization for each enantiomer as a result of their electronic prochirality. We simulate an experiment to demonstrate this effect in cis-stilbene based on its interaction with circularly polarized light. Our results highlight the drastic change in chiral behavior upon electronic excitation, opening up the possibility for asymmetric photochemistry from an effectively nonchiral starting point. The authors report non-adiabatic first principles molecular dynamics to show how an achiral molecule can be converted to a chiral one upon photoexcitation. These results demonstrate the possibility of asymmetric photochemistry starting from achiral reactants.
Collapse
Affiliation(s)
- Umberto Raucci
- Department of Chemistry and The PULSE Institute, Stanford University, Stanford, CA, 94305, USA.,SLAC National Accelerator Laboratory, 2575 Sand Hill Road, Menlo Park, CA, 94025, USA.,Italian Institute of Technology, Genova, GE, Italy
| | - Hayley Weir
- Department of Chemistry and The PULSE Institute, Stanford University, Stanford, CA, 94305, USA.,SLAC National Accelerator Laboratory, 2575 Sand Hill Road, Menlo Park, CA, 94025, USA
| | - Christoph Bannwarth
- Department of Chemistry and The PULSE Institute, Stanford University, Stanford, CA, 94305, USA.,SLAC National Accelerator Laboratory, 2575 Sand Hill Road, Menlo Park, CA, 94025, USA.,Institute of Physical Chemistry, RWTH Aachen University, Aachen, Germany
| | - David M Sanchez
- Department of Chemistry and The PULSE Institute, Stanford University, Stanford, CA, 94305, USA.,SLAC National Accelerator Laboratory, 2575 Sand Hill Road, Menlo Park, CA, 94025, USA.,Design Physics Division, Lawrence Livermore National Laboratory, Livermore, CA, USA
| | - Todd J Martínez
- Department of Chemistry and The PULSE Institute, Stanford University, Stanford, CA, 94305, USA. .,SLAC National Accelerator Laboratory, 2575 Sand Hill Road, Menlo Park, CA, 94025, USA.
| |
Collapse
|
13
|
T. do Casal M, Toldo JM, Pinheiro Jr M, Barbatti M. Fewest switches surface hopping with Baeck-An couplings. OPEN RESEARCH EUROPE 2022; 1:49. [PMID: 37645211 PMCID: PMC10446015 DOI: 10.12688/openreseurope.13624.2] [Citation(s) in RCA: 6] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Accepted: 03/01/2022] [Indexed: 08/31/2023]
Abstract
In the Baeck-An (BA) approximation, first-order nonadiabatic coupling vectors are given in terms of adiabatic energy gaps and the second derivative of the gaps with respect to the coupling coordinate. In this paper, a time-dependent (TD) BA approximation is derived, where the couplings are computed from the energy gaps and their second time-derivatives. TD-BA couplings can be directly used in fewest switches surface hopping, enabling nonadiabatic dynamics with any electronic structure methods able to provide excitation energies and energy gradients. Test results of surface hopping with TD-BA couplings for ethylene and fulvene show that the TD-BA approximation delivers a qualitatively correct picture of the dynamics and a semiquantitative agreement with reference data computed with exact couplings. Nevertheless, TD-BA does not perform well in situations conjugating strong couplings and small velocities. Considered the uncertainties in the method, TD-BA couplings could be a competitive approach for inexpensive, exploratory dynamics with a small trajectories ensemble. We also assessed the potential use of TD-BA couplings for surface hopping dynamics with time-dependent density functional theory (TDDFT), but the results are not encouraging due to singlet instabilities near the crossing seam with the ground state.
Collapse
|
14
|
T. do Casal M, Toldo JM, Pinheiro Jr M, Barbatti M. Fewest switches surface hopping with Baeck-An couplings. OPEN RESEARCH EUROPE 2022; 1:49. [PMID: 37645211 PMCID: PMC10446015 DOI: 10.12688/openreseurope.13624.1] [Citation(s) in RCA: 11] [Impact Index Per Article: 5.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Accepted: 03/01/2022] [Indexed: 03/12/2024]
Abstract
In the Baeck-An (BA) approximation, first-order nonadiabatic coupling vectors are given in terms of adiabatic energy gaps and the second derivative of the gaps with respect to the coupling coordinate. In this paper, a time-dependent (TD) BA approximation is derived, where the couplings are computed from the energy gaps and their second time-derivatives. TD-BA couplings can be directly used in fewest switches surface hopping, enabling nonadiabatic dynamics with any electronic structure methods able to provide excitation energies and energy gradients. Test results of surface hopping with TD-BA couplings for ethylene and fulvene show that the TD-BA approximation delivers a qualitatively correct picture of the dynamics and a semiquantitative agreement with reference data computed with exact couplings. Nevertheless, TD-BA does not perform well in situations conjugating strong couplings and small velocities. Considered the uncertainties in the method, TD-BA couplings could be a competitive approach for inexpensive, exploratory dynamics with a small trajectories ensemble. We also assessed the potential use of TD-BA couplings for surface hopping dynamics with time-dependent density functional theory (TDDFT), but the results are not encouraging due to singlet instabilities near the crossing seam with the ground state.
Collapse
|
15
|
Poštulka J, Slavíček P, Pysanenko A, Poterya V, Fárník M. Bimolecular reactions on sticky and slippery clusters: Electron-induced reactions of hydrogen peroxide. J Chem Phys 2022; 156:054306. [DOI: 10.1063/5.0079283] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/15/2023] Open
Affiliation(s)
- Jan Poštulka
- Department of Physical Chemistry, University of Chemistry and Technology, Prague, Technická 5, 16628 Prague, Czech Republic
| | - Petr Slavíček
- Department of Physical Chemistry, University of Chemistry and Technology, Prague, Technická 5, 16628 Prague, Czech Republic
| | - Andriy Pysanenko
- J. Heyrovský Institute of Physical Chemistry, Czech Academy of Sciences, Dolejškova 2155/3, Prague 8, Czech Republic
| | - Viktoriya Poterya
- J. Heyrovský Institute of Physical Chemistry, Czech Academy of Sciences, Dolejškova 2155/3, Prague 8, Czech Republic
| | - Michal Fárník
- J. Heyrovský Institute of Physical Chemistry, Czech Academy of Sciences, Dolejškova 2155/3, Prague 8, Czech Republic
| |
Collapse
|
16
|
Westermayr J, Marquetand P. Machine Learning for Electronically Excited States of Molecules. Chem Rev 2021; 121:9873-9926. [PMID: 33211478 PMCID: PMC8391943 DOI: 10.1021/acs.chemrev.0c00749] [Citation(s) in RCA: 162] [Impact Index Per Article: 54.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/17/2020] [Indexed: 12/11/2022]
Abstract
Electronically excited states of molecules are at the heart of photochemistry, photophysics, as well as photobiology and also play a role in material science. Their theoretical description requires highly accurate quantum chemical calculations, which are computationally expensive. In this review, we focus on not only how machine learning is employed to speed up such excited-state simulations but also how this branch of artificial intelligence can be used to advance this exciting research field in all its aspects. Discussed applications of machine learning for excited states include excited-state dynamics simulations, static calculations of absorption spectra, as well as many others. In order to put these studies into context, we discuss the promises and pitfalls of the involved machine learning techniques. Since the latter are mostly based on quantum chemistry calculations, we also provide a short introduction into excited-state electronic structure methods and approaches for nonadiabatic dynamics simulations and describe tricks and problems when using them in machine learning for excited states of molecules.
Collapse
Affiliation(s)
- Julia Westermayr
- Institute
of Theoretical Chemistry, Faculty of Chemistry, University of Vienna, Währinger Strasse 17, 1090 Vienna, Austria
| | - Philipp Marquetand
- Institute
of Theoretical Chemistry, Faculty of Chemistry, University of Vienna, Währinger Strasse 17, 1090 Vienna, Austria
- Vienna
Research Platform on Accelerating Photoreaction Discovery, University of Vienna, Währinger Strasse 17, 1090 Vienna, Austria
- Data
Science @ Uni Vienna, University of Vienna, Währinger Strasse 29, 1090 Vienna, Austria
| |
Collapse
|
17
|
Abstract
Electronically excited states of molecules are at the heart of photochemistry, photophysics, as well as photobiology and also play a role in material science. Their theoretical description requires highly accurate quantum chemical calculations, which are computationally expensive. In this review, we focus on not only how machine learning is employed to speed up such excited-state simulations but also how this branch of artificial intelligence can be used to advance this exciting research field in all its aspects. Discussed applications of machine learning for excited states include excited-state dynamics simulations, static calculations of absorption spectra, as well as many others. In order to put these studies into context, we discuss the promises and pitfalls of the involved machine learning techniques. Since the latter are mostly based on quantum chemistry calculations, we also provide a short introduction into excited-state electronic structure methods and approaches for nonadiabatic dynamics simulations and describe tricks and problems when using them in machine learning for excited states of molecules.
Collapse
Affiliation(s)
- Julia Westermayr
- Institute of Theoretical Chemistry, Faculty of Chemistry, University of Vienna, Währinger Strasse 17, 1090 Vienna, Austria
| | - Philipp Marquetand
- Institute of Theoretical Chemistry, Faculty of Chemistry, University of Vienna, Währinger Strasse 17, 1090 Vienna, Austria
- Vienna Research Platform on Accelerating Photoreaction Discovery, University of Vienna, Währinger Strasse 17, 1090 Vienna, Austria
- Data Science @ Uni Vienna, University of Vienna, Währinger Strasse 29, 1090 Vienna, Austria
| |
Collapse
|
18
|
Matsika S. Electronic Structure Methods for the Description of Nonadiabatic Effects and Conical Intersections. Chem Rev 2021; 121:9407-9449. [PMID: 34156838 DOI: 10.1021/acs.chemrev.1c00074] [Citation(s) in RCA: 50] [Impact Index Per Article: 16.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/07/2023]
Abstract
Nonadiabatic effects are ubiquitous in photophysics and photochemistry, and therefore, many theoretical developments have been made to properly describe them. Conical intersections are central in nonadiabatic processes, as they promote efficient and ultrafast nonadiabatic transitions between electronic states. A proper theoretical description requires developments in electronic structure and specifically in methods that describe conical intersections between states and nonadiabatic coupling terms. This review focuses on the electronic structure aspects of nonadiabatic processes. We discuss the requirements of electronic structure methods to describe conical intersections and nonadiabatic couplings, how the most common excited state methods perform in describing these effects, and what the recent developments are in expanding the methodology and implementing nonadiabatic couplings.
Collapse
Affiliation(s)
- Spiridoula Matsika
- Department of Chemistry, Temple University, Philadelphia, Pennsylvania 19122, United States
| |
Collapse
|
19
|
Williams M, Forbes R, Weir H, Veyrinas K, MacDonell RJ, Boguslavskiy AE, Schuurman MS, Stolow A, Martinez TJ. Unmasking the cis-Stilbene Phantom State via Vacuum Ultraviolet Time-Resolved Photoelectron Spectroscopy and Ab Initio Multiple Spawning. J Phys Chem Lett 2021; 12:6363-6369. [PMID: 34231356 DOI: 10.1021/acs.jpclett.1c01227] [Citation(s) in RCA: 13] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/13/2023]
Abstract
We present the first vacuum ultraviolet time-resolved photoelectron spectroscopy (VUV-TRPES) study of photoisomerization dynamics in the paradigmatic molecule cis-stilbene. A key reaction intermediate in its dynamics, known as the phantom state, has often been invoked but never directly detected in the gas phase. We report direct spectral signatures of the phantom state in isolated cis-stilbene, observed and characterized through a combination of VUV-TRPES and ab initio multiple spawning (AIMS) nonadiabatic dynamics simulations of the channel-resolved observable. The high VUV probe photon energy tracks the complete excited-state dynamics via multiple photoionization channels, from initial excitation to its return to the "hot" ground state. The TRPES was compared with AIMS simulations of the dynamics from initial excitation, to the phantom-state intermediate (an S1 minimum), through to the ultimate electronic decay to the ground state. This combination revealed the unique spectral signatures and time-dependent dynamics of the phantom-state intermediate, permitting us to report here its direct observation.
Collapse
Affiliation(s)
- Monika Williams
- Department of Chemistry and The PULSE Institute, Stanford University, Stanford, California 94305, United States
- SLAC National Accelerator Laboratory, 2575 Sand Hill Road, Menlo Park, California 94025, United States
| | - Ruaridh Forbes
- SLAC National Accelerator Laboratory, 2575 Sand Hill Road, Menlo Park, California 94025, United States
- Department of Physics, University of Ottawa, 150 Louis Pasteur, Ottawa, ON K1N 6N5, Canada
| | - Hayley Weir
- Department of Chemistry and The PULSE Institute, Stanford University, Stanford, California 94305, United States
- SLAC National Accelerator Laboratory, 2575 Sand Hill Road, Menlo Park, California 94025, United States
| | - Kévin Veyrinas
- Department of Physics, University of Ottawa, 150 Louis Pasteur, Ottawa, ON K1N 6N5, Canada
| | - Ryan J MacDonell
- Department of Chemistry and Biomolecular Sciences, University of Ottawa, 150 Louis Pasteur, Ottawa, ON K1N 6N5, Canada
- National Research Council Canada, 100 Sussex Drive, Ottawa, ON K1A 0R6, Canada
| | - Andrey E Boguslavskiy
- Department of Physics, University of Ottawa, 150 Louis Pasteur, Ottawa, ON K1N 6N5, Canada
- Department of Chemistry and Biomolecular Sciences, University of Ottawa, 150 Louis Pasteur, Ottawa, ON K1N 6N5, Canada
- National Research Council Canada, 100 Sussex Drive, Ottawa, ON K1A 0R6, Canada
- Joint Centre for Extreme Photonics, NRC-University of Ottawa, Ottawa, ON K1A 0R6, Canada
| | - Michael S Schuurman
- Department of Chemistry and Biomolecular Sciences, University of Ottawa, 150 Louis Pasteur, Ottawa, ON K1N 6N5, Canada
- National Research Council Canada, 100 Sussex Drive, Ottawa, ON K1A 0R6, Canada
| | - Albert Stolow
- Department of Physics, University of Ottawa, 150 Louis Pasteur, Ottawa, ON K1N 6N5, Canada
- Department of Chemistry and Biomolecular Sciences, University of Ottawa, 150 Louis Pasteur, Ottawa, ON K1N 6N5, Canada
- National Research Council Canada, 100 Sussex Drive, Ottawa, ON K1A 0R6, Canada
- Joint Centre for Extreme Photonics, NRC-University of Ottawa, Ottawa, ON K1A 0R6, Canada
| | - Todd J Martinez
- Department of Chemistry and The PULSE Institute, Stanford University, Stanford, California 94305, United States
- SLAC National Accelerator Laboratory, 2575 Sand Hill Road, Menlo Park, California 94025, United States
| |
Collapse
|
20
|
Hu D, Xie Y, Peng J, Lan Z. On-the-Fly Symmetrical Quasi-Classical Dynamics with Meyer-Miller Mapping Hamiltonian for the Treatment of Nonadiabatic Dynamics at Conical Intersections. J Chem Theory Comput 2021; 17:3267-3279. [PMID: 34028268 DOI: 10.1021/acs.jctc.0c01249] [Citation(s) in RCA: 10] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
The on-the-fly version of the symmetrical quasi-classical dynamics method based on the Meyer-Miller mapping Hamiltonian (SQC/MM) is implemented to study the nonadiabatic dynamics at conical intersections of polyatomic systems. The current on-the-fly implementation of the SQC/MM method is based on the adiabatic representation and the dressed momentum. To include the zero-point energy (ZPE) correction of the electronic mapping variables, we employ both the γ-adjusted and γ-fixed approaches. Nonadiabatic dynamics of the methaniminium cation (CH2NH2+) and azomethane are simulated using the on-the-fly SQC/MM method. For CH2NH2+, both ZPE correction approaches give reasonable and consistent results. However, for azomethane, the γ-adjusted version of the SQC/MM dynamics behaves much better than the γ-fixed version. Further analysis indicates that it is always recommended to use the γ-adjusted SQC/MM dynamics in the on-the-fly simulation of photoinduced dynamics of polyatomic systems, particularly when the excited state is well separated from the ground state in the Franck-Condon region. This work indicates that the on-the-fly SQC/MM method is a powerful simulation protocol to deal with the nonadiabatic dynamics of realistic polyatomic systems.
Collapse
Affiliation(s)
- Deping Hu
- SCNU Environmental Research Institute, Guangdong Provincial Key Laboratory of Chemical Pollution and Environmental Safety & MOE Key Laboratory of Environmental Theoretical Chemistry, South China Normal University, Guangzhou 510006, China.,School of Environment, South China Normal University, Guangzhou 510006, China
| | - Yu Xie
- SCNU Environmental Research Institute, Guangdong Provincial Key Laboratory of Chemical Pollution and Environmental Safety & MOE Key Laboratory of Environmental Theoretical Chemistry, South China Normal University, Guangzhou 510006, China.,School of Environment, South China Normal University, Guangzhou 510006, China
| | - Jiawei Peng
- SCNU Environmental Research Institute, Guangdong Provincial Key Laboratory of Chemical Pollution and Environmental Safety & MOE Key Laboratory of Environmental Theoretical Chemistry, South China Normal University, Guangzhou 510006, China.,School of Environment, South China Normal University, Guangzhou 510006, China
| | - Zhenggang Lan
- SCNU Environmental Research Institute, Guangdong Provincial Key Laboratory of Chemical Pollution and Environmental Safety & MOE Key Laboratory of Environmental Theoretical Chemistry, South China Normal University, Guangzhou 510006, China.,School of Environment, South China Normal University, Guangzhou 510006, China
| |
Collapse
|
21
|
Barbatti M. Velocity Adjustment in Surface Hopping: Ethylene as a Case Study of the Maximum Error Caused by Direction Choice. J Chem Theory Comput 2021; 17:3010-3018. [PMID: 33844922 DOI: 10.1021/acs.jctc.1c00012] [Citation(s) in RCA: 21] [Impact Index Per Article: 7.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/22/2022]
Abstract
The most common surface hopping dynamics algorithms require velocity adjustment after hopping to ensure total-energy conservation. Based on the semiclassical analysis, this adjustment must be made parallel to the nonadiabatic coupling vector's direction. Nevertheless, this direction is not always known, and the common practice has been to adjust the velocity in either the linear momentum or velocity directions. This paper benchmarks surface hopping dynamics of photoexcited ethylene with velocity adjustment in several directions, including those of the nonadiabatic coupling vector, the momentum, and the energy gradient difference. It is shown that differences in time constants and structural evolution fall within the statistical uncertainty of the method considering up to 500 trajectories in each dynamics set, rendering the three approaches statistically equivalent. For larger ensembles beyond 1000 trajectories, significant differences between the results arise, limiting the validity of adjustment in alternative directions. Other possible adjustment directions (velocity, single-state gradients, angular momentum) are evaluated as well. Given the small size of ethylene, the results reported in this paper should be considered an upper limit for the error caused by the choice of the velocity-adjustment direction on surface hopping dynamics.
Collapse
|
22
|
Pieri E, Lahana D, Chang AM, Aldaz CR, Thompson KC, Martínez TJ. The non-adiabatic nanoreactor: towards the automated discovery of photochemistry. Chem Sci 2021; 12:7294-7307. [PMID: 34163820 PMCID: PMC8171323 DOI: 10.1039/d1sc00775k] [Citation(s) in RCA: 25] [Impact Index Per Article: 8.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022] Open
Abstract
The ab initio nanoreactor has previously been introduced to automate reaction discovery for ground state chemistry. In this work, we present the nonadiabatic nanoreactor, an analogous framework for excited state reaction discovery. We automate the study of nonadiabatic decay mechanisms of molecules by probing the intersection seam between adiabatic electronic states with hyper-real metadynamics, sampling the branching plane for relevant conical intersections, and performing seam-constrained path searches. We illustrate the effectiveness of the nonadiabatic nanoreactor by applying it to benzene, a molecule with rich photochemistry and a wide array of photochemical products. Our study confirms the existence of several types of S0/S1 and S1/S2 conical intersections which mediate access to a variety of ground state stationary points. We elucidate the connections between conical intersection energy/topography and the resulting photoproduct distribution, which changes smoothly along seam space segments. The exploration is performed with minimal user input, and the protocol requires no previous knowledge of the photochemical behavior of a target molecule. We demonstrate that the nonadiabatic nanoreactor is a valuable tool for the automated exploration of photochemical reactions and their mechanisms.
Collapse
Affiliation(s)
- Elisa Pieri
- Department of Chemistry, The PULSE Institute, Stanford University Stanford CA 94305 USA .,SLAC National Accelerator Laboratory 2575 Sand Hill Road Menlo Park CA 94025 USA
| | - Dean Lahana
- Department of Chemistry, The PULSE Institute, Stanford University Stanford CA 94305 USA .,SLAC National Accelerator Laboratory 2575 Sand Hill Road Menlo Park CA 94025 USA
| | - Alexander M Chang
- Department of Chemistry, The PULSE Institute, Stanford University Stanford CA 94305 USA .,SLAC National Accelerator Laboratory 2575 Sand Hill Road Menlo Park CA 94025 USA
| | - Cody R Aldaz
- Department of Chemistry, The PULSE Institute, Stanford University Stanford CA 94305 USA .,SLAC National Accelerator Laboratory 2575 Sand Hill Road Menlo Park CA 94025 USA
| | - Keiran C Thompson
- Department of Chemistry, The PULSE Institute, Stanford University Stanford CA 94305 USA .,SLAC National Accelerator Laboratory 2575 Sand Hill Road Menlo Park CA 94025 USA
| | - Todd J Martínez
- Department of Chemistry, The PULSE Institute, Stanford University Stanford CA 94305 USA .,SLAC National Accelerator Laboratory 2575 Sand Hill Road Menlo Park CA 94025 USA
| |
Collapse
|
23
|
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
|
24
|
Levine BG, Durden AS, Esch MP, Liang F, Shu Y. CAS without SCF-Why to use CASCI and where to get the orbitals. J Chem Phys 2021; 154:090902. [PMID: 33685182 DOI: 10.1063/5.0042147] [Citation(s) in RCA: 22] [Impact Index Per Article: 7.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/17/2022] Open
Abstract
The complete active space self-consistent field (CASSCF) method has seen broad adoption due to its ability to describe the electronic structure of both the ground and excited states of molecules over a broader swath of the potential energy surface than is possible with the simpler Hartree-Fock approximation. However, it also has a reputation for being unwieldy, computationally costly, and un-black-box. Here, we discuss a class of alternatives, complete active space configuration interaction (CASCI) methods, paying particular attention to their application to electronic excited states. The goal of this Perspective is fourfold. First, we argue that CASCI is not merely an approximation to CASSCF, in that it can be designed to have important qualitative advantages over CASSCF. Second, we present several insights drawn from our experience experimenting with different schemes for computing orbitals to be employed in CASCI. Third, we argue that CASCI is well suited for application to nanomaterials. Finally, we reason that, with the rise in new low-scaling approaches for describing multireference systems, there is a greater need than ever to develop new methods for defining orbitals that provide an efficient and accurate description of both static correlation and electronic excitations in a limited active space.
Collapse
Affiliation(s)
- Benjamin G Levine
- Institute for Advanced Computational Science and Department of Chemistry, Stony Brook University, Stony Brook, New York 11794, USA
| | - Andrew S Durden
- Institute for Advanced Computational Science and Department of Chemistry, Stony Brook University, Stony Brook, New York 11794, USA
| | - Michael P Esch
- Institute for Advanced Computational Science and Department of Chemistry, Stony Brook University, Stony Brook, New York 11794, USA
| | - Fangchun Liang
- Institute for Advanced Computational Science and Department of Chemistry, Stony Brook University, Stony Brook, New York 11794, USA
| | - Yinan Shu
- Department of Chemistry, University of Minnesota, Minneapolis, Minnesota 55455, USA
| |
Collapse
|
25
|
Westermayr J, Marquetand P. Deep learning for UV absorption spectra with SchNarc: First steps toward transferability in chemical compound space. J Chem Phys 2020; 153:154112. [DOI: 10.1063/5.0021915] [Citation(s) in RCA: 33] [Impact Index Per Article: 8.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/06/2023] Open
Affiliation(s)
- J. Westermayr
- Faculty of Chemistry, Institute of Theoretical Chemistry, University of Vienna, Währinger Str. 17, 1090 Vienna, Austria
| | - P. Marquetand
- Faculty of Chemistry, Institute of Theoretical 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
- Faculty of Chemistry, Data Science @ Uni Vienna, University of Vienna, Währinger Str. 29, 1090 Vienna, Austria
| |
Collapse
|
26
|
Farfan CA, Turner DB. A systematic model study quantifying how conical intersection topography modulates photochemical reactions. Phys Chem Chem Phys 2020; 22:20265-20283. [PMID: 32966428 DOI: 10.1039/d0cp03464a] [Citation(s) in RCA: 12] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/13/2022]
Abstract
Despite their important role in photochemistry and expected presence in most polyatomic molecules, conical intersections have been thoroughly characterized in a comparatively small number of systems. Conical intersections can confer molecular photoreactivity or photostability, often with remarkable efficacy, due to their unique structure: at a conical intersection, the adiabatic potential energy surfaces of two or more electronic states are degenerate, enabling ultrafast decay from an excited state without radiative emission, known as nonadiabatic transfer. Furthermore, the precise conical intersection topography determines fundamental properties of photochemical processes, including excited-state decay rate, efficacy, and molecular products that are formed. However, these relationships have yet to be defined comprehensively. In this article, we use an adaptable computational model to investigate a variety of conical intersection topographies, simulate resulting nonadiabatic dynamics, and calculate key photochemical observables. We varied the vibrational mode frequencies to modify conical intersection topography systematically in four primary classes of conical intersections and quantified the resulting rate, total yield, and product yield of nonadiabatic decay. The results reveal that higher vibrational mode frequencies reduce nonadiabatic transfer, but increase the transfer rate and resulting photoproduct formation. These trends can inform progress toward experimental control of photochemical reactions or tuning of molecules' photochemical properties based on conical intersections and their topography.
Collapse
Affiliation(s)
- Camille A Farfan
- Department of Chemistry, New York University, New York, NY 10003, USA
| | - Daniel B Turner
- Department of Chemistry, New York University, New York, NY 10003, USA
| |
Collapse
|
27
|
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
|
28
|
Westermayr J, Marquetand P. Machine learning and excited-state molecular dynamics. MACHINE LEARNING-SCIENCE AND TECHNOLOGY 2020. [DOI: 10.1088/2632-2153/ab9c3e] [Citation(s) in RCA: 27] [Impact Index Per Article: 6.8] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/21/2022]
|
29
|
Tran LN, Neuscamman E. Improving Excited-State Potential Energy Surfaces via Optimal Orbital Shapes. J Phys Chem A 2020; 124:8273-8279. [DOI: 10.1021/acs.jpca.0c07593] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/03/2023]
Affiliation(s)
- Lan Nguyen Tran
- Department of Chemistry, University of California, Berkeley, California 94720, United States
- Ho Chi Minh City Institute of Physics, VAST, Ho Chi Minh City 700000, Vietnam
| | - Eric Neuscamman
- Department of Chemistry, University of California, Berkeley, California 94720, United States
- Chemical Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, United States
| |
Collapse
|
30
|
Seritan S, Bannwarth C, Fales BS, Hohenstein EG, Isborn CM, Kokkila‐Schumacher SIL, Li X, Liu F, Luehr N, Snyder JW, Song C, Titov AV, Ufimtsev IS, Wang L, Martínez TJ. TeraChem
: A graphical processing unit
‐accelerated
electronic structure package for
large‐scale
ab initio molecular dynamics. WILEY INTERDISCIPLINARY REVIEWS-COMPUTATIONAL MOLECULAR SCIENCE 2020. [DOI: 10.1002/wcms.1494] [Citation(s) in RCA: 53] [Impact Index Per Article: 13.3] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/06/2022]
Affiliation(s)
- Stefan Seritan
- Department of Chemistry and the PULSE Institute Stanford University Stanford California USA
- SLAC National Accelerator Laboratory Menlo Park California USA
| | - Christoph Bannwarth
- Department of Chemistry and the PULSE Institute Stanford University Stanford California USA
- SLAC National Accelerator Laboratory Menlo Park California USA
| | - Bryan S. Fales
- Department of Chemistry and the PULSE Institute Stanford University Stanford California USA
- SLAC National Accelerator Laboratory Menlo Park California USA
| | - Edward G. Hohenstein
- Department of Chemistry and the PULSE Institute Stanford University Stanford California USA
- SLAC National Accelerator Laboratory Menlo Park California USA
| | - Christine M. Isborn
- Department of Chemistry University of California Merced Merced California USA
| | | | - Xin Li
- Division of Theoretical Chemistry and Biology, School of Engineering Sciences in Chemistry, Biotechnology and Health KTH Royal Institute of Technology Stockholm Sweden
| | - Fang Liu
- Department of Chemical Engineering Massachusetts Institute of Technology Cambridge Massachusetts USA
| | | | | | - Chenchen Song
- Department of Physics University of California Berkeley Berkeley California USA
- Molecular Foundry Lawrence Berkeley National Laboratory Berkeley California USA
| | | | - Ivan S. Ufimtsev
- Department of Structural Biology Stanford University School of Medicine Stanford California USA
| | - Lee‐Ping Wang
- Department of Chemistry University of California Davis Davis California USA
| | - Todd J. Martínez
- Department of Chemistry and the PULSE Institute Stanford University Stanford California USA
- SLAC National Accelerator Laboratory Menlo Park California USA
| |
Collapse
|
31
|
Suchan J, Janoš J, Slavíček P. Pragmatic Approach to Photodynamics: Mixed Landau–Zener Surface Hopping with Intersystem Crossing. J Chem Theory Comput 2020; 16:5809-5820. [DOI: 10.1021/acs.jctc.0c00512] [Citation(s) in RCA: 25] [Impact Index Per Article: 6.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/16/2022]
Affiliation(s)
- Jiří Suchan
- Department of Physical Chemistry, University of Chemistry and Technology, Prague, Technická 5, 16628 Prague, Czech Republic
| | - Jiří Janoš
- Department of Physical Chemistry, University of Chemistry and Technology, Prague, Technická 5, 16628 Prague, Czech Republic
| | - Petr Slavíček
- Department of Physical Chemistry, University of Chemistry and Technology, Prague, Technická 5, 16628 Prague, Czech Republic
| |
Collapse
|
32
|
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
|
33
|
Seritan S, Bannwarth C, Fales BS, Hohenstein EG, Kokkila-Schumacher SIL, Luehr N, Snyder JW, Song C, Titov AV, Ufimtsev IS, Martínez TJ. TeraChem: Accelerating electronic structure and ab initio molecular dynamics with graphical processing units. J Chem Phys 2020; 152:224110. [PMID: 32534542 PMCID: PMC7928072 DOI: 10.1063/5.0007615] [Citation(s) in RCA: 76] [Impact Index Per Article: 19.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/18/2020] [Accepted: 05/19/2020] [Indexed: 11/15/2022] Open
Abstract
Developed over the past decade, TeraChem is an electronic structure and ab initio molecular dynamics software package designed from the ground up to leverage graphics processing units (GPUs) to perform large-scale ground and excited state quantum chemistry calculations in the gas and the condensed phase. TeraChem's speed stems from the reformulation of conventional electronic structure theories in terms of a set of individually optimized high-performance electronic structure operations (e.g., Coulomb and exchange matrix builds, one- and two-particle density matrix builds) and rank-reduction techniques (e.g., tensor hypercontraction). Recent efforts have encapsulated these core operations and provided language-agnostic interfaces. This greatly increases the accessibility and flexibility of TeraChem as a platform to develop new electronic structure methods on GPUs and provides clear optimization targets for emerging parallel computing architectures.
Collapse
Affiliation(s)
| | | | | | | | | | | | | | | | | | - Ivan S. Ufimtsev
- Department of Structural Biology, Stanford University School of Medicine, Stanford, California 94305, USA
| | | |
Collapse
|
34
|
Peters LM, Kussmann J, Ochsenfeld C. Combining Graphics Processing Units, Simplified Time-Dependent Density Functional Theory, and Finite-Difference Couplings to Accelerate Nonadiabatic Molecular Dynamics. J Phys Chem Lett 2020; 11:3955-3961. [PMID: 32374606 PMCID: PMC7304892 DOI: 10.1021/acs.jpclett.0c00320] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/30/2020] [Accepted: 04/23/2020] [Indexed: 06/11/2023]
Abstract
Starting from our recently published implementation of nonadiabatic molecular dynamics (NAMD) on graphics processing units (GPUs), we explore further approaches to accelerate ab initio NAMD calculations at the time-dependent density functional theory (TDDFT) level of theory. We employ (1) the simplified TDDFT schemes of Grimme et al. and (2) the Hammes-Schiffer-Tully approach to obtain nonadiabatic couplings from finite-difference calculations. The resulting scheme delivers an accurate physical picture while virtually eliminating the two computationally most demanding steps of the algorithm. Combined with our GPU-based integral routines for SCF, TDDFT, and TDDFT derivative calculations, NAMD simulations of systems of a few hundreds of atoms at a reasonable time scale become accessible on a single compute node. To demonstrate this and to present a first, illustrative example, we perform TDDFT/MM-NAMD simulations of the rhodopsin protein.
Collapse
Affiliation(s)
- Laurens
D. M. Peters
- Chair
of Theoretical Chemistry, Department of Chemistry, University of Munich (LMU), Butenandtstr. 7, D-81377 München, Germany
| | - Jörg Kussmann
- Chair
of Theoretical Chemistry, Department of Chemistry, University of Munich (LMU), Butenandtstr. 7, D-81377 München, Germany
| | - Christian Ochsenfeld
- Chair
of Theoretical Chemistry, Department of Chemistry, University of Munich (LMU), Butenandtstr. 7, D-81377 München, Germany
- Max
Planck Institute for Solid State Research, Heisenbergstr. 1, D-70569 Stuttgart, Germany
| |
Collapse
|
35
|
Westermayr J, Faber FA, Christensen AS, von Lilienfeld OA, Marquetand P. Neural networks and kernel ridge regression for excited states dynamics of CH2NH$_2^+$: From single-state to multi-state representations and multi-property machine learning models. MACHINE LEARNING-SCIENCE AND TECHNOLOGY 2020. [DOI: 10.1088/2632-2153/ab88d0] [Citation(s) in RCA: 35] [Impact Index Per Article: 8.8] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/16/2022]
|
36
|
Pysanenko A, Gámez F, Fárník M, Chalabala J, Slavíček P. Photochemistry of Amylene Double Bond in Clusters on Free Argon Nanoparticles. J Phys Chem A 2020; 124:3038-3047. [PMID: 32240587 DOI: 10.1021/acs.jpca.0c00860] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
We have investigated reactivity of double bond in 2-methyl-2-butene (also trimethylethylene or amylene) in the excited and ionized states. In a combined experimental and theoretical study, we focused on both the intermolecular and intramolecular reactions. In a molecular beam experiment, we have sequentially picked up several amylene molecules on the surface of argon nanoparticles ArM, M̅ ≈ 90, acting as a cold support. Ionization with 70 eV electrons yields mass spectra strongly dominated by amylene cluster ions Am(Am)n+. Interestingly, upon multiphoton ionization with 193 nm (6.4 eV) photons, a new strong fragment series appears in the spectra, nominally corresponding to an addition of two carbon atoms, i.e., (Am)nC2+. This difference between electron and photoionization suggests a reaction in an excited state of amylene with a neighboring amylene molecule. We used techniques of nonadiabatic molecular dynamics to study the reactivity of amylene molecules both in the excited and in ionized states. Possible reaction pathways are proposed, substantiating the observed differences between the electron ionization and photoionization mass spectra.
Collapse
Affiliation(s)
- Andriy Pysanenko
- J. Heyrovský Institute of Physical Chemistry, The Czech Academy of Sciences, Dolejškova 3, 182 23 Prague, Czech Republic
| | - Francisco Gámez
- J. Heyrovský Institute of Physical Chemistry, The Czech Academy of Sciences, Dolejškova 3, 182 23 Prague, Czech Republic
| | - Michal Fárník
- J. Heyrovský Institute of Physical Chemistry, The Czech Academy of Sciences, Dolejškova 3, 182 23 Prague, Czech Republic
| | - Jan Chalabala
- University of Chemistry and Technology, 166 28 Prague, Czech Republic
| | - Petr Slavíček
- University of Chemistry and Technology, 166 28 Prague, Czech Republic
| |
Collapse
|
37
|
Peters LM, Kussmann J, Ochsenfeld C. Nonadiabatic Molecular Dynamics on Graphics Processing Units: Performance and Application to Rotary Molecular Motors. J Chem Theory Comput 2019; 15:6647-6659. [PMID: 31763834 PMCID: PMC6909237 DOI: 10.1021/acs.jctc.9b00859] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/29/2019] [Indexed: 11/29/2022]
Abstract
Nonadiabatic molecular dynamics (NAMD) simulations of molecular systems require the efficient evaluation of excited-state properties, such as energies, gradients, and nonadiabatic coupling vectors. Here, we investigate the use of graphics processing units (GPUs) in addition to central processing units (CPUs) to efficiently calculate these properties at the time-dependent density functional theory (TDDFT) level of theory. Our implementation in the FermiONs++ program package uses the J-engine and a preselective screening procedure for the calculation of Coulomb and exchange kernels, respectively. We observe good speed-ups for small and large molecular systems (comparable to those observed in ground-state calculations) and reduced (down to sublinear) scaling behavior with respect to the system size (depending on the spatial locality of the investigated excitation). As a first illustrative application, we present efficient NAMD simulations of a series of newly designed light-driven rotary molecular motors and compare their S1 lifetimes. Although all four rotors show different S1 excitation energies, their ability to rotate upon excitation is conserved, making the series an interesting starting point for rotary molecular motors with tunable excitation energies.
Collapse
Affiliation(s)
- Laurens
D. M. Peters
- Chair
of Theoretical Chemistry, Department of Chemistry, University of Munich (LMU), Butenandtstr. 7, D-81377 München, Germany
| | - Jörg Kussmann
- Chair
of Theoretical Chemistry, Department of Chemistry, University of Munich (LMU), Butenandtstr. 7, D-81377 München, Germany
| | - Christian Ochsenfeld
- Chair
of Theoretical Chemistry, Department of Chemistry, University of Munich (LMU), Butenandtstr. 7, D-81377 München, Germany
- Max
Planck Institute for Solid State Research, Heisenbergstr. 1, D-70569 Stuttgart, Germany
| |
Collapse
|
38
|
Zhou W, Mandal A, Huo P. Quasi-Diabatic Scheme for Nonadiabatic On-the-Fly Simulations. J Phys Chem Lett 2019; 10:7062-7070. [PMID: 31665889 DOI: 10.1021/acs.jpclett.9b02747] [Citation(s) in RCA: 21] [Impact Index Per Article: 4.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/15/2023]
Abstract
We use the quasi-diabatic (QD) propagation scheme to perform on-the-fly nonadiabatic simulations of the photodynamics of ethylene. The QD scheme enables a seamless interface between accurate diabatic-based quantum dynamics approaches and adiabatic electronic structure calculations, explicitly avoiding any efforts to construct global diabatic states or reformulate the diabatic dynamics approach to the adiabatic representation. Using the partial linearized path-integral approach and the symmetrical quasi-classical approach as the diabatic dynamics methods, the QD propagation scheme enables direct nonadiabatic simulation with complete active space self-consistent field on-the-fly electronic structure calculations. The population dynamics obtained from both approaches are in a close agreement with the quantum wavepacket-based method and outperform the widely used trajectory surface-hopping approach. Further analysis of the ethylene photodeactivation pathways demonstrates the correct predictions of competing processes of nonradiative relaxation mechanism through various conical intersections. This work provides the foundation of using accurate diabatic dynamics approaches and on-the-fly adiabatic electronic structure information to perform ab initio nonadiabatic simulation.
Collapse
Affiliation(s)
- Wanghuai Zhou
- Advanced Functional Material and Photoelectric Technology Research Institution, School of Science , Hubei University of Automotive Technology , Shiyan , Hubei 442002 , People's Republic of China
- Department of Chemistry , University of Rochester , 120 Trustee Road , Rochester , New York 14627 , United States
| | - Arkajit Mandal
- Department of Chemistry , University of Rochester , 120 Trustee Road , Rochester , New York 14627 , United States
| | - Pengfei Huo
- Department of Chemistry , University of Rochester , 120 Trustee Road , Rochester , New York 14627 , United States
| |
Collapse
|
39
|
Ibele LM, Nicolson A, Curchod BFE. Excited-state dynamics of molecules with classically driven trajectories and Gaussians. Mol Phys 2019. [DOI: 10.1080/00268976.2019.1665199] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/25/2022]
Affiliation(s)
- Lea M. Ibele
- Department of Chemistry, Durham University, Durham, UK
| | | | | |
Collapse
|
40
|
Farfan CA, Turner DB. Nonadiabatic Photochemistry Induced by Inaccessible Conical Intersections. J Phys Chem A 2019; 123:7768-7776. [DOI: 10.1021/acs.jpca.9b07739] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Camille A. Farfan
- Department of Chemistry, New York University, 100 Washington Square East, New York New York 10003, United States
| | - Daniel B. Turner
- Department of Chemistry, New York University, 100 Washington Square East, New York New York 10003, United States
| |
Collapse
|
41
|
Tran LN, Shea JAR, Neuscamman E. Tracking Excited States in Wave Function Optimization Using Density Matrices and Variational Principles. J Chem Theory Comput 2019; 15:4790-4803. [DOI: 10.1021/acs.jctc.9b00351] [Citation(s) in RCA: 17] [Impact Index Per Article: 3.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Lan Nguyen Tran
- Department of Chemistry, University of California, Berkeley, California 94720, United States
- Ho Chi Minh City Institute of Physics, VAST, Ho Chi Minh City 700000, Vietnam
| | - Jacqueline A. R. Shea
- Department of Chemistry, University of California, Berkeley, California 94720, United States
| | - Eric Neuscamman
- Department of Chemistry, University of California, Berkeley, California 94720, United States
- Chemical Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, United States
| |
Collapse
|
42
|
Thongyod W, Buranachai C, Pengpan T, Punwong C. Fluorescence quenching by photoinduced electron transfer between 7-methoxycoumarin and guanine base facilitated by hydrogen bonds: an in silico study. Phys Chem Chem Phys 2019; 21:16258-16269. [PMID: 31304496 DOI: 10.1039/c9cp02037c] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
In this study, the effects of hydrogen bond (H-bond) formation on fluorescence quenching of 7-methoxycoumarin (7MC) via photo-induced electron transfer from a guanine base (Gua) are investigated using a combined quantum mechanics/molecular mechanics simulation. The electronic structure is calculated by the floating occupation molecular orbital complete active space configuration interaction modification on a semiempirical method. Then the full multiple spawning method is employed for the dynamics simulations on multiple electronic states. The methods employed here are validated by simulating direct dynamics of 7MC (without Gua) and compared with available experimental results. Our computational results are in good agreement with the previously reported experimental results in terms of spectroscopic properties of 7MC. In the case of a H-bonded 7MC-Gua complex, the results from constrained dynamics simulations and single-point calculations suggest that the electron transfer occurs on the second excited state and it depends not only on the H-bond length but also on the intermolecular planarity between 7MC and Gua. Moreover, a proton coupled electron transfer can occur at ≈1 Å of H-bond length, where a proton from Gua is also transferred together with the electron to 7MC. The obtained simulations are expected to be greatly beneficial for designing effective fluorescently labeled nucleotide probes as well as providing information for precise fluorescence signal interpretation.
Collapse
Affiliation(s)
- Wutthinan Thongyod
- Department of Physics, Faculty of Science, Prince of Songkla University, Songkhla 90112, Thailand. and Center of Excellence for Trace Analysis and Biosensor, Prince of Songkla University, Songkhla 90112, Thailand
| | - Chittanon Buranachai
- Department of Physics, Faculty of Science, Prince of Songkla University, Songkhla 90112, Thailand. and Center of Excellence for Trace Analysis and Biosensor, Prince of Songkla University, Songkhla 90112, Thailand
| | - Teparksorn Pengpan
- Department of Physics, Faculty of Science, Prince of Songkla University, Songkhla 90112, Thailand.
| | - Chutintorn Punwong
- Department of Physics, Faculty of Science, Prince of Songkla University, Songkhla 90112, Thailand.
| |
Collapse
|
43
|
Accomasso D, Persico M, Granucci G. Diabatization by Localization in the Framework of Configuration Interaction Based on Floating Occupation Molecular Orbitals (FOMO−CI). CHEMPHOTOCHEM 2019. [DOI: 10.1002/cptc.201900056] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/22/2023]
Affiliation(s)
- Davide Accomasso
- Dipartimento di Chimica e Chimica IndustrialeUniversitá di Pisa v. G. Moruzzi 13 I-56124 Pisa Italy
| | - Maurizio Persico
- Dipartimento di Chimica e Chimica IndustrialeUniversitá di Pisa v. G. Moruzzi 13 I-56124 Pisa Italy
| | - Giovanni Granucci
- Dipartimento di Chimica e Chimica IndustrialeUniversitá di Pisa v. G. Moruzzi 13 I-56124 Pisa Italy
| |
Collapse
|
44
|
Parrish RM, Martínez TJ. Ab Initio Computation of Rotationally-Averaged Pump-Probe X-ray and Electron Diffraction Signals. J Chem Theory Comput 2019; 15:1523-1537. [PMID: 30702882 DOI: 10.1021/acs.jctc.8b01051] [Citation(s) in RCA: 21] [Impact Index Per Article: 4.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
We develop a new algorithm for the computation of the rotationally averaged elastic molecular diffraction signal for the cases of perpendicular or parallel pump-probe geometries. The algorithm first collocates the charge density from an arbitrary ab initio wave function onto a Becke quadrature grid [A. Becke, J. Chem. Phys. 1988 , 88 , 2457 ], providing a high-fidelity multiresolution representation of the charge density. A double sum is then performed over the Becke grid points, and the interaction between points computed using the scattering kernels of Williamson and Zewail [J. C. Williamson and A. H. Zewail, J. Phys. Chem. 1994 , 98 , 2766 ]. These kernels analytically average over the molecular orientations with the cos2 γ selection factor appropriate for one-photon dipole absorption in a perpendicular pump-probe geometry. We show that the method is converged with small grids containing <500 points/atom. We implement the algorithm on a GPU for increased efficiency and emonstrate the algorithm for molecules with up to a few dozen atoms. We explore the accuracy of the independent atom model (IAM) by comparison with our new and more accurate method. We also investigate the possibility of detecting signatures of electronic transitions in polyatomic pump-probe diffraction experiments.
Collapse
Affiliation(s)
- Robert M Parrish
- Department of Chemistry and The PULSE Institute , Stanford University , Stanford , California 94305 , United States.,SLAC National Accelerator Laboratory , 2575 Sand Hill Road , Menlo Park , California 94025 , United States
| | - Todd J Martínez
- Department of Chemistry and The PULSE Institute , Stanford University , Stanford , California 94305 , United States.,SLAC National Accelerator Laboratory , 2575 Sand Hill Road , Menlo Park , California 94025 , United States
| |
Collapse
|
45
|
Ranković M, Chalabala J, Zawadzki M, Kočišek J, Slavíček P, Fedor J. Dissociative ionization dynamics of dielectric gas C3F7CN. Phys Chem Chem Phys 2019; 21:16451-16458. [PMID: 31312828 DOI: 10.1039/c9cp02188d] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Fluoronitrile C3F7CN is a promising candidate for the replacement of SF6 dielectric gas in high-voltage insulation.
Collapse
Affiliation(s)
- M. Ranković
- J. Heyrovský Institute of Physical Chemistry
- Czech Academy of Sciences
- 18223 Prague
- Czech Republic
| | - J. Chalabala
- Department of Physical Chemistry
- University of Chemistry and Technology
- 16628 Prague
- Czech Republic
| | - M. Zawadzki
- J. Heyrovský Institute of Physical Chemistry
- Czech Academy of Sciences
- 18223 Prague
- Czech Republic
- Department of Atomic, Molecular, and Optical Physics
| | - J. Kočišek
- J. Heyrovský Institute of Physical Chemistry
- Czech Academy of Sciences
- 18223 Prague
- Czech Republic
| | - P. Slavíček
- J. Heyrovský Institute of Physical Chemistry
- Czech Academy of Sciences
- 18223 Prague
- Czech Republic
- Department of Physical Chemistry
| | - J. Fedor
- J. Heyrovský Institute of Physical Chemistry
- Czech Academy of Sciences
- 18223 Prague
- Czech Republic
| |
Collapse
|
46
|
Cozza C, Bonomi M, Pietropaolo A. A Versatile Computational Strategy To Characterize the Free-Energy Landscape of Excited States in Oligofluorenes. J Chem Theory Comput 2018; 14:5441-5445. [PMID: 30365313 DOI: 10.1021/acs.jctc.8b00949] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
The π-rich rings of conjugated polymers and molecular rotors shape their typical properties, allowing a variety of chemical and photoresponsive phenomena. Herein, we present a computational method in the framework of classical simulations to estimate the free-energy gap between ground and excited states of oligofluorenes, which is a class of molecular rotors widely used in optoelectronic devices, because of the inner torsional rotation triggered by light irradiation. We devised multiple sets of free-energy simulations in combination with free-energy perturbation theory to predict the free-energy gap between the ground state and the first excited state. The computed excitation energies show good agreement with experiments. The approach presented herein allows one to achieve comprehensive sampling of the conformational landscape and accurate estimates of the excited state free-energy landscapes at the same time.
Collapse
Affiliation(s)
- Concetta Cozza
- Dipartimento di Scienze della Salute , Università di Catanzaro , Viale Europa , 88100 Catanzaro , Italy
| | - Massimiliano Bonomi
- Department of Chemistry , University of Cambridge , Lensfield Road , Cambridge CB2 1EW , United Kingdom
| | - Adriana Pietropaolo
- Dipartimento di Scienze della Salute , Università di Catanzaro , Viale Europa , 88100 Catanzaro , Italy
| |
Collapse
|
47
|
|
48
|
Lischka H, Nachtigallová D, Aquino AJA, Szalay PG, Plasser F, Machado FBC, Barbatti M. Multireference Approaches for Excited States of Molecules. Chem Rev 2018; 118:7293-7361. [DOI: 10.1021/acs.chemrev.8b00244] [Citation(s) in RCA: 197] [Impact Index Per Article: 32.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/18/2022]
Affiliation(s)
- Hans Lischka
- School of Pharmaceutical Sciences and Technology, Tianjin University, Tianjin 300072, P.R. China
- Department of Chemistry and Biochemistry, Texas Tech University, Lubbock, Texas 79409, United States
- Institute of Theoretical Chemistry, Faculty of Chemistry, University of Vienna, Währinger Straße 17, 1090 Vienna, Austria
| | - Dana Nachtigallová
- Institute of Organic Chemistry and Biochemistry v.v.i., The Czech Academy of Sciences, Flemingovo nám. 2, 16610 Prague 6, Czech Republic
- Regional Centre of Advanced Technologies and Materials, Palacký University, 78371 Olomouc, Czech Republic
| | - Adélia J. A. Aquino
- School of Pharmaceutical Sciences and Technology, Tianjin University, Tianjin 300072, P.R. China
- Department of Chemistry and Biochemistry, Texas Tech University, Lubbock, Texas 79409, United States
- Institute for Soil Research, University of Natural Resources and Life Sciences Vienna, Peter-Jordan-Strasse 82, A-1190 Vienna, Austria
| | - Péter G. Szalay
- ELTE Eötvös Loránd University, Laboratory of Theoretical Chemistry, Pázmány Péter sétány 1/A, 1117 Budapest, Hungary
| | - Felix Plasser
- Institute of Theoretical Chemistry, Faculty of Chemistry, University of Vienna, Währinger Straße 17, 1090 Vienna, Austria
- Department of Chemistry, Loughborough University, Leicestershire LE11 3TU, United Kingdom
| | - Francisco B. C. Machado
- Departamento de Química, Instituto Tecnológico de Aeronáutica, São José dos Campos 12228-900, São Paulo, Brazil
| | | |
Collapse
|
49
|
Yue L, Liu Y, Zhu C. Performance of TDDFT with and without spin-flip in trajectory surface hopping dynamics: cis–trans azobenzene photoisomerization. Phys Chem Chem Phys 2018; 20:24123-24139. [DOI: 10.1039/c8cp03851a] [Citation(s) in RCA: 25] [Impact Index Per Article: 4.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Quantum yields simulated by LR-TDDFT and SF-TDDFT methods for azobenzene photoisomerizations.
Collapse
Affiliation(s)
- Ling Yue
- Department of Applied Chemistry and Institute of Molecular Science
- National Chiao-Tung University
- Hsinchu 30010
- Taiwan
- MOE Key Laboratory for Nonequilibrium Synthesis and Modulation of Condensed Matter
| | - Yajun Liu
- Key Laboratory of Theoretical and Computational Photochemistry
- Ministry of Education
- College of Chemistry
- Beijing Normal University
- Beijing 100875
| | - Chaoyuan Zhu
- Department of Applied Chemistry and Institute of Molecular Science
- National Chiao-Tung University
- Hsinchu 30010
- Taiwan
- Center for Emergent Functional Matter Science
| |
Collapse
|
50
|
Suchan J, Hollas D, Curchod BFE, Slavíček P. On the importance of initial conditions for excited-state dynamics. Faraday Discuss 2018; 212:307-330. [DOI: 10.1039/c8fd00088c] [Citation(s) in RCA: 27] [Impact Index Per Article: 4.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/18/2022]
Abstract
The vast majority of ab initio excited-state simulations are performed within semiclassical, trajectory-based approaches. Apart from the underlying electronic-structure theory, the reliability of the simulations is controlled by a selection of initial conditions for the classical trajectories. We discuss appropriate choices of initial conditions for simulations of different experimental arrangements: dynamics initiated by continuum-wave (CW) laser fields or triggered by ultrashort laser pulses.
Collapse
Affiliation(s)
- Jiří Suchan
- Department of Physical Chemistry
- University of Chemistry and Technology, Prague
- 16628 Prague
- Czech Republic
| | - Daniel Hollas
- Department of Physical Chemistry
- University of Chemistry and Technology, Prague
- 16628 Prague
- Czech Republic
| | | | - Petr Slavíček
- Department of Physical Chemistry
- University of Chemistry and Technology, Prague
- 16628 Prague
- Czech Republic
| |
Collapse
|