1
|
Vandaele E, Mališ M, Luber S. A Local Diabatisation Method for Two-State Adiabatic Conical Intersections. J Chem Theory Comput 2024; 20:856-872. [PMID: 38174710 DOI: 10.1021/acs.jctc.3c01008] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/05/2024]
Abstract
A methodology to locally characterize conical intersections (CIs) between two adiabatic electronic states for which no nonadiabatic coupling (NAC) vectors are available is presented. Based on the Hessian and gradient at the CI, the branching space coordinates are identified. The potential energy surface around the CI in the branching space is expressed in the diabatic representation, from which the NAC vectors can be calculated in a wave-function-free, energy-based approach. To demonstrate the universality of the developed methodology, the minimum-energy CI (MECI) between the first (S1) and second (S2) singlet excited states of formamide is investigated at the state-averaged complete active space self-consistent field (SA-CASSCF) and extended multistate complete active space second-order perturbation theory (XMS-CASPT2) levels of theory. In addition, the asymmetrical MECI between the ground state (S0) and S1 of cyclopropanone is evaluated using SA-CASSCF, as well as (ME)CIs between the S1 and S2 states of benzene using SA-CASSCF and time-dependent density functional theory (TDDFT). Finally, a CI between the S1 and S2 excited states of thiophene was analyzed using TDDFT.
Collapse
Affiliation(s)
- Eva Vandaele
- Department of Chemistry, University of Zürich, Winterthurerstrasse 190, 8057 Zürich, Switzerland
| | - Momir Mališ
- Department of Chemistry, University of Zürich, Winterthurerstrasse 190, 8057 Zürich, Switzerland
| | - Sandra Luber
- Department of Chemistry, University of Zürich, Winterthurerstrasse 190, 8057 Zürich, Switzerland
| |
Collapse
|
2
|
Temen S, Akimov AV. A Simple Solution to Trivial Crossings: A Stochastic State Tracking Approach. J Phys Chem Lett 2021; 12:850-860. [PMID: 33427475 DOI: 10.1021/acs.jpclett.0c03428] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/12/2023]
Abstract
We present a new state tracking algorithm based on a stochastic state reassignment that reflects the quantum mechanical interpretation of the state time-overlaps. We assess the new method with a range of model Hamiltonians and demonstrate that it yields the results generally consistent with the deterministic min-cost algorithm. However, the stochastic state tracking algorithm reduces magnitudes of the state population fluctuations as the quantum system evolves toward its equilibrium. The new algorithm facilitates the thermalization of quantum state populations and suppresses the population revivals and oscillations near the equilibrium in many-state systems. The new stochastic algorithm has a favorable computational scaling, is easy to implement due to its conceptual transparency, and treats various types of state identity changes (trivial or avoided crossings and any intermediate cases) on equal footing.
Collapse
Affiliation(s)
- Story Temen
- Department of Chemistry, University at Buffalo, The State University of New York, Buffalo, New York 14260, United States
| | - Alexey V Akimov
- Department of Chemistry, University at Buffalo, The State University of New York, Buffalo, New York 14260, United States
| |
Collapse
|
3
|
Rivera M, Dommett M, Crespo-Otero R. ONIOM(QM:QM′) Electrostatic Embedding Schemes for Photochemistry in Molecular Crystals. J Chem Theory Comput 2019; 15:2504-2516. [DOI: 10.1021/acs.jctc.8b01180] [Citation(s) in RCA: 17] [Impact Index Per Article: 3.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/18/2022]
Affiliation(s)
- Miguel Rivera
- School of Biological and Chemical Sciences, Queen Mary University of London, Mile End Road, London E1 4NS, U.K
| | - Michael Dommett
- School of Biological and Chemical Sciences, Queen Mary University of London, Mile End Road, London E1 4NS, U.K
| | - Rachel Crespo-Otero
- School of Biological and Chemical Sciences, Queen Mary University of London, Mile End Road, London E1 4NS, U.K
| |
Collapse
|
4
|
Three-state conical intersection optimization methods: development and implementation at QM/MM level. Theor Chem Acc 2016. [DOI: 10.1007/s00214-016-2029-z] [Citation(s) in RCA: 15] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/20/2022]
|
5
|
Improta R, Santoro F, Blancafort L. Quantum Mechanical Studies on the Photophysics and the Photochemistry of Nucleic Acids and Nucleobases. Chem Rev 2016; 116:3540-93. [PMID: 26928320 DOI: 10.1021/acs.chemrev.5b00444] [Citation(s) in RCA: 338] [Impact Index Per Article: 42.3] [Reference Citation Analysis] [Abstract] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/24/2022]
Abstract
The photophysics and photochemistry of DNA is of great importance due to the potential damage of the genetic code by UV light. Quantum mechanical studies have played a key role in interpretating the results of modern time-resolved pump-probe spectroscopy, and in elucidating the main photoactivated reactive paths. This review provides a concise, complete picture of the computational studies carried out, approximately, in the past decade. We start with an overview of the photophysics of the nucleobases in the gas phase and in solution. We discuss the proposed mechanisms for ultrafast decay to the ground state, that involve conical intersections, consider the role of triplet states, and analyze how the solvent modulates the photophysics. Then we move to larger systems, from dinucleotides to single- and double-stranded oligonucleotides. We focus on the possible role of charge transfer and delocalized or excitonic states in the photophysics of these systems and discuss the main photochemical paths. We finish with an outlook on the current challenges in the field and future directions of research.
Collapse
Affiliation(s)
- Roberto Improta
- Istituto di Biostrutture Biommagini (IBB-CNR), CNR-Consiglio Nazionale delle Ricerche , Via Mezzocannone 16, I-80134, Napoli, Italy
| | - Fabrizio Santoro
- Area della Ricerca di Pisa, Istituto di Chimica dei Composti Organo Metallici (ICCOM-CNR), CNR-Consiglio Nazionale delle Ricerche , Via G. Moruzzi 1, I-56124 Pisa, Italy
| | - Lluís Blancafort
- Institut de Química Computacional i Catàlisi and Departament de Química, Universitat de Girona, Campus de Montilivi , 17071 Girona, Spain
| |
Collapse
|
6
|
Xia SH, Liu XY, Fang Q, Cui G. Photodissociation dynamics of CH3C(O)SH in argon matrix: A QM/MM nonadiabatic dynamics simulation. J Chem Phys 2015; 143:194303. [DOI: 10.1063/1.4935598] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/13/2022] Open
|
7
|
Cui G, Thiel W. Generalized trajectory surface-hopping method for internal conversion and intersystem crossing. J Chem Phys 2015; 141:124101. [PMID: 25273406 DOI: 10.1063/1.4894849] [Citation(s) in RCA: 117] [Impact Index Per Article: 13.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/30/2023] Open
Abstract
Trajectory-based fewest-switches surface-hopping (FSSH) dynamics simulations have become a popular and reliable theoretical tool to simulate nonadiabatic photophysical and photochemical processes. Most available FSSH methods model internal conversion. We present a generalized trajectory surface-hopping (GTSH) method for simulating both internal conversion and intersystem crossing processes on an equal footing. We consider hops between adiabatic eigenstates of the non-relativistic electronic Hamiltonian (pure spin states), which is appropriate for sufficiently small spin-orbit coupling. This choice allows us to make maximum use of existing electronic structure programs and to minimize the changes to available implementations of the traditional FSSH method. The GTSH method is formulated within the quantum mechanics (QM)/molecular mechanics framework, but can of course also be applied at the pure QM level. The algorithm implemented in the GTSH code is specified step by step. As an initial GTSH application, we report simulations of the nonadiabatic processes in the lowest four electronic states (S0, S1, T1, and T2) of acrolein both in vacuo and in acetonitrile solution, in which the acrolein molecule is treated at the ab initio complete-active-space self-consistent-field level. These dynamics simulations provide detailed mechanistic insight by identifying and characterizing two nonadiabatic routes to the lowest triplet state, namely, direct S1 → T1 hopping as major pathway and sequential S1 → T2 → T1 hopping as minor pathway, with the T2 state acting as a relay state. They illustrate the potential of the GTSH approach to explore photoinduced processes in complex systems, in which intersystem crossing plays an important role.
Collapse
Affiliation(s)
- Ganglong Cui
- Max-Planck-Institut für Kohlenforschung, Kaiser-Wilhelm-Platz 1, 45470 Mülheim an der Ruhr, Germany
| | - Walter Thiel
- Max-Planck-Institut für Kohlenforschung, Kaiser-Wilhelm-Platz 1, 45470 Mülheim an der Ruhr, Germany
| |
Collapse
|
8
|
Minezawa N. Optimizing minimum free-energy crossing points in solution: Linear-response free energy/spin-flip density functional theory approach. J Chem Phys 2014; 141:164118. [DOI: 10.1063/1.4899049] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Noriyuki Minezawa
- Fukui Institute for Fundamental Chemistry, Kyoto University, Sakyo-ku, Kyoto 606-8103, Japan
| |
Collapse
|
9
|
Ruiz-Barragan S, Robb MA, Blancafort L. Conical Intersection Optimization Based on a Double Newton–Raphson Algorithm Using Composed Steps. J Chem Theory Comput 2013; 9:1433-42. [DOI: 10.1021/ct301059t] [Citation(s) in RCA: 18] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Sergi Ruiz-Barragan
- Institut de Química Computacional
and Department de Química, University of Girona, 17071 Girona,
Spain
| | - Michael A. Robb
- Department of Chemistry,
Imperial
College, London SW7 2AZ, United Kingdom
| | - Lluís Blancafort
- Institut de Química Computacional
and Department de Química, University of Girona, 17071 Girona,
Spain
| |
Collapse
|
10
|
Mori T, Martínez TJ. Exploring the Conical Intersection Seam: The Seam Space Nudged Elastic Band Method. J Chem Theory Comput 2013; 9:1155-63. [DOI: 10.1021/ct300892t] [Citation(s) in RCA: 39] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Toshifumi Mori
- PULSE Institute and
Department
of Chemistry, 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
- PULSE Institute and
Department
of Chemistry, Stanford University, Stanford, California 94305, United
States
- SLAC National Accelerator
Laboratory,
2575 Sand Hill Road, Menlo Park, California 94025, United States
| |
Collapse
|
11
|
Sparta M, Shirvanyants D, Ding F, Dokholyan NV, Alexandrova AN. Hybrid dynamics simulation engine for metalloproteins. Biophys J 2013; 103:767-76. [PMID: 22947938 DOI: 10.1016/j.bpj.2012.06.024] [Citation(s) in RCA: 24] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/12/2012] [Revised: 06/08/2012] [Accepted: 06/18/2012] [Indexed: 11/15/2022] Open
Abstract
Quality computational description of metalloproteins is a great challenge due to the vast span of time- and lengthscales characteristic of their existence. We present an efficient new method that allows for robust characterization of metalloproteins. It combines quantum mechanical (QM) description of the metal-containing active site, and extensive dynamics of the protein captured by discrete molecular dynamics (DMD) (QM/DMD). DMD samples the entire protein, including the backbone, and most of the active site, except for the immediate coordination region of the metal. QM operates on the part of the protein of electronic and chemical significance, which may include tens to hundreds of atoms. The breathing quantum-classical boundary provides a continuous mutual feedback between the two machineries. We test QM/DMD using the Fe-containing electron transporter protein, rubredoxin, and its three mutants as a model. QM/DMD can provide a reliable balanced description of metalloproteins' structure, dynamics, and electronic structure in a reasonable amount of time. As an illustration of QM/DMD capabilities, we then predict the structure of the Ca(2+) form of the enzyme catechol O-methyl transferase, which, unlike the native Mg(2+) form, is catalytically inactive. The Mg(2+) site is ochtahedral, but the Ca(2+) is 7-coordinate and features the misalignment of the reacting parts of the system. The change is facilitated by the backbone adjustment. QM/DMD is ideal and fast for providing this level of structural insight.
Collapse
Affiliation(s)
- Manuel Sparta
- Department of Chemistry and Biochemistry, University of California, Los Angeles, California, USA
| | | | | | | | | |
Collapse
|
12
|
Minezawa N, Gordon MS. Optimizing conical intersections of solvated molecules: The combined spin-flip density functional theory/effective fragment potential method. J Chem Phys 2012; 137:034116. [DOI: 10.1063/1.4734314] [Citation(s) in RCA: 28] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|