1
|
Huang H, Peng J, Zhang Y, Gu FL, Lan Z, Xu C. The development of the QM/MM interface and its application for the on-the-fly QM/MM nonadiabatic dynamics in JADE package: Theory, implementation, and applications. J Chem Phys 2024; 160:234101. [PMID: 38884395 DOI: 10.1063/5.0215036] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/22/2024] [Accepted: 05/15/2024] [Indexed: 06/18/2024] Open
Abstract
Understanding the nonadiabatic dynamics of complex systems is a challenging task in computational photochemistry. Herein, we present an efficient and user-friendly quantum mechanics/molecular mechanics (QM/MM) interface to run on-the-fly nonadiabatic dynamics. Currently, this interface consists of an independent set of codes designed for general-purpose use. Herein, we demonstrate the ability and feasibility of the QM/MM interface by integrating it with our long-term developed JADE package. Tailored to handle nonadiabatic processes in various complex systems, especially condensed phases and protein environments, we delve into the theories, implementations, and applications of on-the-fly QM/MM nonadiabatic dynamics. The QM/MM approach is established within the framework of the additive QM/MM scheme, employing electrostatic embedding, link-atom inclusion, and charge-redistribution schemes to treat the QM/MM boundary. Trajectory surface-hopping dynamics are facilitated using the fewest switches algorithm, encompassing classical and quantum treatments for nuclear and electronic motions, respectively. Finally, we report simulations of nonadiabatic dynamics for two typical systems: azomethane in water and the retinal chromophore PSB3 in a protein environment. Our results not only illustrate the power of the QM/MM program but also reveal the important roles of environmental factors in nonadiabatic processes.
Collapse
Affiliation(s)
- Haiyi Huang
- MOE Key Laboratory of Environmental Theoretical Chemistry and Guangdong Provincial Key Laboratory of Chemical Pollution and Environmental Safety, SCNU Environmental Research Institute, School of Environment, South China Normal University, Guangzhou 510006, China
- Center for Advanced Materials Research, Beijing Normal University, Zhuhai 519087, China
- MOE Key Laboratory of Theoretical and Computational Photochemistry, College of Chemistry, Beijing Normal University, Beijing 100875, China
| | - Jiawei Peng
- MOE Key Laboratory of Environmental Theoretical Chemistry and Guangdong Provincial Key Laboratory of Chemical Pollution and Environmental Safety, SCNU Environmental Research Institute, School of Environment, South China Normal University, Guangzhou 510006, China
| | - Yulin Zhang
- MOE Key Laboratory of Environmental Theoretical Chemistry and Guangdong Provincial Key Laboratory of Chemical Pollution and Environmental Safety, SCNU Environmental Research Institute, School of Environment, South China Normal University, Guangzhou 510006, China
| | - Feng Long Gu
- MOE Key Laboratory of Environmental Theoretical Chemistry and Guangdong Provincial Key Laboratory of Chemical Pollution and Environmental Safety, SCNU Environmental Research Institute, School of Environment, South China Normal University, Guangzhou 510006, China
| | - Zhenggang Lan
- MOE Key Laboratory of Environmental Theoretical Chemistry and Guangdong Provincial Key Laboratory of Chemical Pollution and Environmental Safety, SCNU Environmental Research Institute, School of Environment, South China Normal University, Guangzhou 510006, China
| | - Chao Xu
- MOE Key Laboratory of Environmental Theoretical Chemistry and Guangdong Provincial Key Laboratory of Chemical Pollution and Environmental Safety, SCNU Environmental Research Institute, School of Environment, South China Normal University, Guangzhou 510006, China
| |
Collapse
|
2
|
Palombo R, Barneschi L, Pedraza-González L, Yang X, Olivucci M. Picosecond quantum-classical dynamics reveals that the coexistence of light-induced microbial and animal chromophore rotary motion modulates the isomerization quantum yield of heliorhodopsin. Phys Chem Chem Phys 2024; 26:10343-10356. [PMID: 38501246 DOI: 10.1039/d4cp00193a] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 03/20/2024]
Abstract
Rhodopsins are light-responsive proteins forming two vast and evolutionary distinct superfamilies whose functions are invariably triggered by the photoisomerization of a single retinal chromophore. In 2018 a third widespread superfamily of rhodopsins called heliorhodopsins was discovered using functional metagenomics. Heliorhodopsins, with their markedly different structural features with respect to the animal and microbial superfamilies, offer an opportunity to study how evolution has manipulated the chromophore photoisomerization to achieve adaptation. One question is related to the mechanism of such a reaction and how it differs from that of animal and microbial rhodopsins. To address this question, we use hundreds of quantum-classical trajectories to simulate the spectroscopically documented picosecond light-induced dynamics of a heliorhodopsin from the archaea thermoplasmatales archaeon (TaHeR). We show that, consistently with the observations, the trajectories reveal two excited state decay channels. However, inconsistently with previous hypotheses, only one channel is associated with the -C13C14- rotation of microbial rhodopsins while the second channel is characterized by the -C11C12- rotation typical of animal rhodopsins. The fact that such -C11C12- rotation is aborted upon decay and ground state relaxation, explains why illumination of TaHeR only produces the 13-cis isomer with a low quantum efficiency. We argue that the documented lack of regioselectivity in double-bond excited state twisting motion is the result of an "adaptation" that could be completely lost via specific residue substitutions modulating the steric hindrance experienced along the isomerization motion.
Collapse
Affiliation(s)
- Riccardo Palombo
- Dipartimento di Biotecnologie, Chimica e Farmacia, Università di Siena, via A. Moro 2, I-53100 Siena, Siena, Italy.
- Department of Chemistry, Bowling Green State University, Bowling Green, Ohio 43403, USA.
| | - Leonardo Barneschi
- Dipartimento di Biotecnologie, Chimica e Farmacia, Università di Siena, via A. Moro 2, I-53100 Siena, Siena, Italy.
| | - Laura Pedraza-González
- Dipartimento di Chimica e Chimica Industriale, Università di Pisa, Via Giuseppe Moruzzi, 13, I-56124 Pisa, Italy
| | - Xuchun Yang
- Department of Chemistry, Bowling Green State University, Bowling Green, Ohio 43403, USA.
| | - Massimo Olivucci
- Dipartimento di Biotecnologie, Chimica e Farmacia, Università di Siena, via A. Moro 2, I-53100 Siena, Siena, Italy.
- Department of Chemistry, Bowling Green State University, Bowling Green, Ohio 43403, USA.
| |
Collapse
|
3
|
Taylor JT, Tozer DJ, Curchod BFE. On the description of conical intersections between excited electronic states with LR-TDDFT and ADC(2). J Chem Phys 2023; 159:214115. [PMID: 38059547 DOI: 10.1063/5.0176140] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/11/2023] [Accepted: 11/14/2023] [Indexed: 12/08/2023] Open
Abstract
Conical intersections constitute the conceptual bedrock of our working understanding of ultrafast, nonadiabatic processes within photochemistry (and photophysics). Accurate calculation of potential energy surfaces within the vicinity of conical intersections, however, still poses a serious challenge to many popular electronic structure methods. Multiple works have reported on the deficiency of methods like linear-response time-dependent density functional theory within the adiabatic approximation (AA LR-TDDFT) or algebraic diagrammatic construction to second-order [ADC(2)]-approaches often used in excited-state molecular dynamics simulations-to describe conical intersections between the ground and excited electronic states. In the present study, we focus our attention on conical intersections between excited electronic states and probe the ability of AA LR-TDDFT and ADC(2) to describe their topology and topography, using protonated formaldimine and pyrazine as two exemplar molecules. We also take the opportunity to revisit the performance of these methods in describing conical intersections involving the ground electronic state in protonated formaldimine-highlighting in particular how the intersection ring exhibited by AA LR-TDDFT can be perceived either as a (near-to-linear) seam of intersection or two interpenetrating cones, depending on the magnitude of molecular distortions within the branching space.
Collapse
Affiliation(s)
- Jack T Taylor
- Department of Chemistry, Durham University, South Road, Durham DH1 3LE, United Kingdom
| | - David J Tozer
- Department of Chemistry, Durham University, South Road, Durham DH1 3LE, United Kingdom
| | - Basile F E Curchod
- Centre for Computational Chemistry, School of Chemistry, University of Bristol, Cantock's Close, Bristol BS8 1TS, United Kingdom
| |
Collapse
|
4
|
Barneschi L, Kaliakin D, Huix-Rotllant M, Ferré N, Filatov Gulak M, Olivucci M. Assessment of the Electron Correlation Treatment on the Quantum-Classical Dynamics of Retinal Protonated Schiff Base Models: XMS-CASPT2, RMS-CASPT2, and REKS Methods. J Chem Theory Comput 2023; 19:8189-8200. [PMID: 37937990 DOI: 10.1021/acs.jctc.3c00879] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/09/2023]
Abstract
We compare the performance of three different multiconfigurational wave function-based electronic structure methods and two implementations of the spin-restricted ensemble-referenced Kohn-Sham (REKS) method. The study is characterized by three features: (i) it uses a small set of quantum-classical trajectories rather than potential energy surface mapping, (ii) it focuses, exclusively, on the photoisomerization of retinal protonated Schiff base models, and (iii) it probes the effect of both methyl substitution and the increase in length of the conjugate π-system. For each tested method, the corresponding analytical gradients are used to drive the quantum-classical (Tully's FSSH method) trajectory propagation, including the recent multistate XMS-CASPT2 and RMS-CASPT2 gradients. It is shown that while CASSCF, XMS-CASPT2, and RMS-CASPT2 yield consistent photoisomerization dynamics descriptions, REKS produces, in some of these systems, qualitatively different behavior that is attributed to a flatter and topographically different excited state potential energy surface. The origin of this behavior can be traced back to the effect of the employed density functional approximation. The above studies are further expanded by benchmarking, at the CASSCF and REKS levels, the electronic structure methods using a QM/MM model of the visual pigment rhodopsin.
Collapse
Affiliation(s)
- Leonardo Barneschi
- Dipartimento di Biotecnologie, Chimica e Farmacia, Università di Siena, I-53100 Siena, Italy
| | - Danil Kaliakin
- Department of Chemistry, Bowling Green State University, Bowling Green, Ohio 43403, United States
| | - Miquel Huix-Rotllant
- Aix-Marseille Université, CNRS, Institut Chimie Radicalaire, 13013 Marseille, France
| | - Nicolas Ferré
- Aix-Marseille Université, CNRS, Institut Chimie Radicalaire, 13013 Marseille, France
| | - Michael Filatov Gulak
- Department of Chemistry, Southern Methodist University, Dallas, Texas 75275, United States
| | - Massimo Olivucci
- Dipartimento di Biotecnologie, Chimica e Farmacia, Università di Siena, I-53100 Siena, Italy
- Department of Chemistry, Bowling Green State University, Bowling Green, Ohio 43403, United States
| |
Collapse
|
5
|
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
|
6
|
Filiba O, Borin VA, Schapiro I. The involvement of triplet states in the isomerization of retinaloids. Phys Chem Chem Phys 2022; 24:26223-26231. [PMID: 36278932 DOI: 10.1039/d2cp03791b] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/28/2023]
Abstract
Rhodopsins form a family of photoreceptor proteins which utilize the retinal chromophore for light energy conversion. Upon light absorption the retinal chromophore undergoes a photoisomerization. This reaction involves a non-radiative relaxation through a conical intersection between the singlet excited state and the ground state. In this work we studied the possible involvement of triplet states in the photoisomerization of retinaloids using the extended multistate (XMS) version of CASPT2. To this end, truncated models of three retinaloids were considered: protonated Schiff base, deprotonated Schiff base and the aldehyde form. The optimized geometries of the reactant, the product and the conical intersection were connected by a linear interpolation of internal coordinates to describe the isomerization. The energetic position of the low-lying singlet and triplet states as well as their spin-orbit coupling matrix elements (SOCME) were calculated along the isomerization profile. The SOCME values peaked in vicinity of the conical intersection for all the retinaloids. Furthermore, the magnitude of SOCME is invariant to the number of double bonds in the model. The SOCME for the protonated Schiff base is negligible (1.5 cm-1) which renders the involvement of the triplet state as improbable. However, the largest SOCME value of 30 cm-1 was found for the aldehyde form, followed by 15 cm-1 for the deprotonated Schiff base.
Collapse
Affiliation(s)
- Ofer Filiba
- Fritz Haber Research Center for Molecular Dynamics and Institute of Chemistry, The Hebrew University of Jerusalem, 9190401, Jerusalem, Israel.
| | - Veniamin A Borin
- Fritz Haber Research Center for Molecular Dynamics and Institute of Chemistry, The Hebrew University of Jerusalem, 9190401, Jerusalem, Israel.
| | - Igor Schapiro
- Fritz Haber Research Center for Molecular Dynamics and Institute of Chemistry, The Hebrew University of Jerusalem, 9190401, Jerusalem, Israel.
| |
Collapse
|
7
|
Quantum-classical simulations of rhodopsin reveal excited-state population splitting and its effects on quantum efficiency. Nat Chem 2022; 14:441-449. [PMID: 35241801 PMCID: PMC8983576 DOI: 10.1038/s41557-022-00892-6] [Citation(s) in RCA: 13] [Impact Index Per Article: 6.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/28/2020] [Accepted: 01/12/2022] [Indexed: 01/12/2023]
Abstract
The activation of rhodopsin, the light-sensitive G-protein coupled receptor responsible for dim-light vision in vertebrates, is driven by an ultrafast excited-state double-bond isomerization with a quantum efficiency of almost 70%. The origin of such light sensitivity is not understood and a key question is whether in-phase nuclear motion controls the quantum efficiency value. Here, we use hundreds of quantum-classical trajectories to show that, 15 femtoseconds after light absorption, a degeneracy between the reactive excited state and a neighboring state causes the splitting of the rhodopsin population into subpopulations. These subpopulations propagate with different velocities and lead to distinct contributions to the quantum efficiency. We also show that such splitting is modulated by protein electrostatics, thus linking amino-acid sequence variations to quantum efficiency modulation. Finally, we discuss how such a linkage that in principle could be exploited to achieve higher quantum efficiencies, would simultaneously increase the receptor thermal noise leading to a trade-off that may have played a role in rhodopsin evolution.
Collapse
|
8
|
de Castro Araujo Valente D, Borges I, Cardozo TM. Nonradiative relaxation mechanisms of the elusive silole molecule. Phys Chem Chem Phys 2021; 23:26561-26574. [PMID: 34811562 DOI: 10.1039/d1cp03803f] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Silole derivatives have been extensively employed for developing organic optoelectronics, but few studies focused on the photophysical properties of the silole molecule. In this work, we investigate these properties by computing the absorption spectra and performing nonadiabatic molecular dynamics of silole employing the algebraic diagrammatic construction [ADC(2)] and extended multi-state XMS-CASPT2 ab initio electronic structure methods. For vertical excitations and excited state optimizations, the equation of motion coupled-cluster singles and doubles (EOM-CCSD) was also used. The nuclear ensemble and the fewest-switches surface hopping molecular dynamics methods coupled with the first two high-level electronic structure methods were applied to probe the relaxation mechanisms of silole. We could reproduce the experimental first absorption maximum value and found an ultrafast relaxation process occurring exclusively through ring-puckering distortions without breaking ring bonds or hydrogen elimination. Minimum energy conical intersection optimizations were carried out and potential energy curves, including triplet states, were calculated to further elucidate the relaxation process of silole.
Collapse
Affiliation(s)
| | - Itamar Borges
- Departamento de Química, Instituto Militar de Engenharia, Rio de Janeiro, RJ, Brazil
| | - Thiago Messias Cardozo
- Instituto de Química, Universidade Federal do Rio de Janeiro, Rio de Janeiro, RJ, Brazil.
| |
Collapse
|
9
|
Segarra‐Martí J, Bearpark MJ. Modelling Photoionisation in Isocytosine: Potential Formation of Longer-Lived Excited State Cations in its Keto Form. Chemphyschem 2021; 22:2172-2181. [PMID: 34370368 PMCID: PMC8597144 DOI: 10.1002/cphc.202100402] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/25/2021] [Revised: 08/09/2021] [Indexed: 11/24/2022]
Abstract
Studying the effects of UV and VUV radiation on non-canonical DNA/RNA nucleobases allows us to compare how they release excess energy following absorption with respect to their canonical counterparts. This has attracted much research attention in recent years because of its likely influence on the origin of our genetic lexicon in prebiotic times. Here we present a CASSCF and XMS-CASPT2 theoretical study of the photoionisation of non-canonical pyrimidine nucleobase isocytosine in both its keto and enol tautomeric forms. We analyse their lowest energy cationic excited states including 2 π + , 2 n O + and 2 n N + and compare these to the corresponding electronic states in cytosine. Investigating lower-energy decay pathways we find - unexpectedly - that keto-isocytosine+ presents a sizeable energy barrier potentially inhibiting decay to its cationic ground state, whereas enol-isocytosine+ features a barrierless and consequently ultrafast pathway analogous to the one previously found for the canonical (keto) form of cytosine+ . Dynamic electron correlation reduces the energy barrier in the keto form substantially (by ∼1 eV) but it is nevertheless still present. We additionally compute the UV/Vis absorption signals of the structures encountered along these decay channels to provide spectroscopic fingerprints to assist future experiments in monitoring these intricate photo-processes.
Collapse
Affiliation(s)
- Javier Segarra‐Martí
- Department of ChemistryMolecular Sciences Research HubImperial College LondonWhite City Campus, 82 Wood LaneLondonW12 0BZUK
- Present address: Instituto de Ciencia MolecularUniversitat de ValenciaP.O. Box 22085ES-46071ValenciaSpain
| | - Michael J. Bearpark
- Department of ChemistryMolecular Sciences Research HubImperial College LondonWhite City Campus, 82 Wood LaneLondonW12 0BZUK
| |
Collapse
|
10
|
Kochman MA, Palczewski K, Kubas A. Theoretical Study of the Photoisomerization Mechanism of All- Trans-Retinyl Acetate. J Phys Chem A 2021; 125:8358-8372. [PMID: 34546761 PMCID: PMC8488936 DOI: 10.1021/acs.jpca.1c05533] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/04/2022]
Abstract
![]()
The compound 9-cis-retinyl acetate (9-cis-RAc) is a precursor
to 9-cis-retinal,
which has potential application in the treatment of some hereditary
diseases of the retina. An attractive synthetic route to 9-cis-RAc is based on the photoisomerization reaction of the
readily available all-trans-RAc. In the present study,
we examine the mechanism of the photoisomerization reaction with the
use of state-of-the-art electronic structure calculations for two
polyenic model compounds: tEtEt-octatetraene and tEtEtEc-2,6-dimethyl-1,3,5,7,9-decapentaene. The occurrence
of photoisomerization is attributed to a chain-kinking mechanism,
whereby a series of S1/S0 conical intersections
associated with kinking deformations at different positions along
the polyenic chain mediate internal conversion to the S0 state, and subsequent isomerization around one of the double bonds.
Two other possible photoisomerization mechanisms are taken into account,
but they are rejected as incompatible with simulation results and/or
the available spectroscopic data.
Collapse
Affiliation(s)
- Michał Andrzej Kochman
- Institute of Physical Chemistry, Polish Academy of Sciences, Ul. Marcina Kasprzaka 44/52, 01-224 Warszawa, Poland
| | - Krzysztof Palczewski
- Department of Ophthalmology, Gavin Herbert Eye Institute, University of California, Irvine, California 92697, United States.,Department of Physiology and Biophysics, University of California, Irvine, California 92697, United States.,Department of Chemistry, University of California, Irvine, California 92697, United States
| | - Adam Kubas
- Institute of Physical Chemistry, Polish Academy of Sciences, Ul. Marcina Kasprzaka 44/52, 01-224 Warszawa, Poland
| |
Collapse
|
11
|
Ariyageadsakul P, Baeck KK. Wavepacket propagations for the early time dynamics of proton-coupled electron transfer in the charge-transfer state of NH 3Cl complex. J Chem Phys 2021; 154:154305. [PMID: 33887927 DOI: 10.1063/5.0046247] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
A charge-transfer (CT) excited state of NH3Cl, generated by photo-detachment of an electron from the anionic NH3Cl- precursor, can be represented as H2N+-H-Cl- and proceeds to two chemical reactions: one reaction generating NH2 and HCl resulting from a proton transfer (PT) and the other reaction producing NH3 and a Cl atom resulting from an electron transfer (ET); both are coupled to form a typical proton-coupled electron transfer (PCET) process. The early time dynamics of this CT were studied using time-dependent wavepacket propagation on three nonadiabatically coupled electronic states in a reduced three-dimensional space. The electronic states were treated using the XMS-CASPT2/aug-cc-pVTZ ab initio methodology. The population dynamics of the three coupled electronic states were analyzed in detail to reveal the initial stage of the PCET process up to ∼100 fs, while the branching ratio, χ = PT/(ET+PT), was determined after wavepacket propagations of up to 2000 fs. Another main result is the dependence of χ on the vibration levels of the initial precursor anion and the isotope substitution of the connecting H atom with deuterium and tritium. Our study reveals the detailed microscopic features of the PCET process embedded in the CT state of the NH3Cl complex and certain systematic dependences of the branching ratio χ on the above factors.
Collapse
Affiliation(s)
- Pinit Ariyageadsakul
- Department of Chemistry, Gangneung-Wonju National University, Gangneung, Gangwon-do 25457, Republic of Korea
| | - Kyoung Koo Baeck
- Department of Chemistry, Gangneung-Wonju National University, Gangneung, Gangwon-do 25457, Republic of Korea
| |
Collapse
|
12
|
Aleotti F, Aranda D, Yaghoubi Jouybari M, Garavelli M, Nenov A, Santoro F. Parameterization of a linear vibronic coupling model with multiconfigurational electronic structure methods to study the quantum dynamics of photoexcited pyrene. J Chem Phys 2021; 154:104106. [PMID: 33722019 DOI: 10.1063/5.0044693] [Citation(s) in RCA: 13] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/20/2022] Open
Abstract
With this work, we present a protocol for the parameterization of a Linear Vibronic Coupling (LVC) Hamiltonian for quantum dynamics using highly accurate multiconfigurational electronic structure methods such as RASPT2/RASSCF, combined with a maximum-overlap diabatization technique. Our approach is fully portable and can be applied to many medium-size rigid molecules whose excited state dynamics requires a quantum description. We present our model and discuss the details of the electronic structure calculations needed for the parameterization, analyzing critical situations that could arise in the case of strongly interacting excited states. The protocol was applied to the simulation of the excited state dynamics of the pyrene molecule, starting from either the first or the second bright state (S2 or S5). The LVC model was benchmarked against state-of-the-art quantum mechanical calculations with optimizations and energy scans and turned out to be very accurate. The dynamics simulations, performed including all active normal coordinates with the multilayer multiconfigurational time-dependent Hartree method, show good agreement with the available experimental data, endorsing prediction of the excited state mechanism, especially for S5, whose ultrafast deactivation mechanism was not yet clearly understood.
Collapse
Affiliation(s)
- Flavia Aleotti
- Dipartimento di Chimica Industriale "Toso Montanari," Università di Bologna, Viale del Risorgimento 4, 40136 Bologna, Italy
| | - Daniel Aranda
- Istituto di Chimica dei Composti Organometallici (ICCOM-CNR), Area della Ricerca del CNR, Via Moruzzi 1, I-56124 Pisa, Italy
| | - Martha Yaghoubi Jouybari
- Istituto di Chimica dei Composti Organometallici (ICCOM-CNR), Area della Ricerca del CNR, Via Moruzzi 1, I-56124 Pisa, Italy
| | - Marco Garavelli
- Dipartimento di Chimica Industriale "Toso Montanari," Università di Bologna, Viale del Risorgimento 4, 40136 Bologna, Italy
| | - Artur Nenov
- Dipartimento di Chimica Industriale "Toso Montanari," Università di Bologna, Viale del Risorgimento 4, 40136 Bologna, Italy
| | - Fabrizio Santoro
- Istituto di Chimica dei Composti Organometallici (ICCOM-CNR), Area della Ricerca del CNR, Via Moruzzi 1, I-56124 Pisa, Italy
| |
Collapse
|
13
|
Park JW, Al-Saadon R, MacLeod MK, Shiozaki T, Vlaisavljevich B. Multireference Electron Correlation Methods: Journeys along Potential Energy Surfaces. Chem Rev 2020; 120:5878-5909. [PMID: 32239929 DOI: 10.1021/acs.chemrev.9b00496] [Citation(s) in RCA: 54] [Impact Index Per Article: 13.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/16/2022]
Abstract
Multireference electron correlation methods describe static and dynamical electron correlation in a balanced way and, therefore, can yield accurate and predictive results even when single-reference methods or multiconfigurational self-consistent field theory fails. One of their most prominent applications in quantum chemistry is the exploration of potential energy surfaces. This includes the optimization of molecular geometries, such as equilibrium geometries and conical intersections and on-the-fly photodynamics simulations, both of which depend heavily on the ability of the method to properly explore the potential energy surface. Because such applications require nuclear gradients and derivative couplings, the availability of analytical nuclear gradients greatly enhances the scope of quantum chemical methods. This review focuses on the developments and advances made in the past two decades. A detailed account of the analytical nuclear gradient and derivative coupling theories is presented. Emphasis is given to the software infrastructure that allows one to make use of these methods. Notable applications of multireference electron correlation methods to chemistry, including geometry optimizations and on-the-fly dynamics, are summarized at the end followed by a discussion of future prospects.
Collapse
Affiliation(s)
- Jae Woo Park
- Department of Chemistry, Chungbuk National University, Chungdae-ro 1, Cheongju 28644, Korea
| | - Rachael Al-Saadon
- Department of Chemistry, Northwestern University, 2145 Sheridan Road, Evanston, Illinois 60208, United States
| | - Matthew K MacLeod
- Workday, 4900 Pearl Circle East, Suite 100, Boulder, Colorado 80301, United States
| | - Toru Shiozaki
- Department of Chemistry, Northwestern University, 2145 Sheridan Road, Evanston, Illinois 60208, United States.,Quantum Simulation Technologies, Inc., 625 Massachusetts Avenue, Cambridge, Massachusetts 02139, United States
| | - Bess Vlaisavljevich
- Department of Chemistry, University of South Dakota, 414 East Clark Street, Vermillion, South Dakota 57069, United States
| |
Collapse
|
14
|
Gromov EV, Domratcheva T. Four resonance structures elucidate double-bond isomerisation of a biological chromophore. Phys Chem Chem Phys 2020; 22:8535-8544. [DOI: 10.1039/d0cp00814a] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/13/2022]
Abstract
Four resonance structures determining the electronic structure of the chromophore’s ground and first excited states. Changing the relative energies of the structures by hydrogen-bonding interactions tunes all chromophore’s photochemical properties.
Collapse
Affiliation(s)
- Evgeniy V. Gromov
- Max-Planck Institute for Medical Research
- Jahnstraße 29
- 69120 Heidelberg
- Germany
| | - Tatiana Domratcheva
- Max-Planck Institute for Medical Research
- Jahnstraße 29
- 69120 Heidelberg
- Germany
| |
Collapse
|
15
|
Park JW. Analytical Gradient Theory for Quasidegenerate N-Electron Valence State Perturbation Theory (QD-NEVPT2). J Chem Theory Comput 2019; 16:326-339. [DOI: 10.1021/acs.jctc.9b00919] [Citation(s) in RCA: 15] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/17/2022]
Affiliation(s)
- Jae Woo Park
- Department of Chemistry, Chungbuk National University (CBNU), Cheongju 28644, Korea
| |
Collapse
|
16
|
Segarra‐Martí J, Tran T, Bearpark MJ. Computing the Ultrafast and Radiationless Electronic Excited State Decay of Cytosine and 5‐methyl‐cytosine Cations: Uncovering the Role of Dynamic Electron Correlation. CHEMPHOTOCHEM 2019. [DOI: 10.1002/cptc.201900105] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/14/2022]
Affiliation(s)
- Javier Segarra‐Martí
- Department of Chemistry, Molecular Sciences Research HubImperial College London White City Campus, 80 Wood Lane W12 0BZ London UK
| | - Thierry Tran
- Department of Chemistry, Molecular Sciences Research HubImperial College London White City Campus, 80 Wood Lane W12 0BZ London UK
| | - Michael J. Bearpark
- Department of Chemistry, Molecular Sciences Research HubImperial College London White City Campus, 80 Wood Lane W12 0BZ London UK
| |
Collapse
|
17
|
Park JW, Al-Saadon R, Strand NE, Shiozaki T. Imaginary Shift in CASPT2 Nuclear Gradient and Derivative Coupling Theory. J Chem Theory Comput 2019; 15:4088-4098. [DOI: 10.1021/acs.jctc.9b00368] [Citation(s) in RCA: 20] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/03/2023]
Affiliation(s)
- Jae Woo Park
- Department of Chemistry, Northwestern University, 2145 Sheridan Rd., Evanston, Illinois 60208, United States
- Department of Chemistry, Chungbuk National University, Chungdae-ro 1, Cheongju, Chungbuk 28644, Korea
| | - Rachael Al-Saadon
- Department of Chemistry, Northwestern University, 2145 Sheridan Rd., Evanston, Illinois 60208, United States
| | - Nils E. Strand
- Department of Chemistry, Northwestern University, 2145 Sheridan Rd., Evanston, Illinois 60208, United States
| | - Toru Shiozaki
- Department of Chemistry, Northwestern University, 2145 Sheridan Rd., Evanston, Illinois 60208, United States
| |
Collapse
|
18
|
Park JW. Single-State Single-Reference and Multistate Multireference Zeroth-Order Hamiltonians in MS-CASPT2 and Conical Intersections. J Chem Theory Comput 2019; 15:3960-3973. [DOI: 10.1021/acs.jctc.9b00067] [Citation(s) in RCA: 13] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Jae Woo Park
- Department of Chemistry, Chungbuk National University (CBNU), Cheongju 28644, Korea
| |
Collapse
|
19
|
Segarra-Martí J, Tran T, Bearpark MJ. Ultrafast and radiationless electronic excited state decay of uracil and thymine cations: computing the effects of dynamic electron correlation. Phys Chem Chem Phys 2019; 21:14322-14330. [PMID: 30698175 DOI: 10.1039/c8cp07189f] [Citation(s) in RCA: 13] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/02/2023]
Abstract
In this article we characterise the radiationless decay of the first few electronic excited states of the cations of DNA/RNA nucleobases uracil and thymine, including the effects of dynamic electron correlation on energies and geometries (optimised with XMS-CASPT2). In both systems, we find that one state of 2n and another two of 2π+ character can be populated following photoionisation, and their different minima and interstate crossings are located. We find strong similarities between uracil and thymine cations: with accessible conical intersections suggesting that depopulation of their electronic excited states takes place on ultrafast timescales in both systems, suggesting that they are photostable in agreement with previous theoretical (uracil+) evidence. We find that dynamic electron correlation separates the energy levels of the "3-state" conical intersection (D2/D1/D0)CI previously located with CASSCF for uracil+, which will therefore have a different geometry and higher energy. Simulating the electronic and vibrational absorptions allows us to characterise spectral fingerprints that could be used to monitor these cation photo-processes experimentally.
Collapse
Affiliation(s)
- Javier Segarra-Martí
- Department of Chemistry, Molecular Sciences Research Hub, Imperial College London, White City Campus, 80 Wood Lane, W12 0BZ, London, UK.
| | - Thierry Tran
- Department of Chemistry, Molecular Sciences Research Hub, Imperial College London, White City Campus, 80 Wood Lane, W12 0BZ, London, UK.
| | - Michael J Bearpark
- Department of Chemistry, Molecular Sciences Research Hub, Imperial College London, White City Campus, 80 Wood Lane, W12 0BZ, London, UK.
| |
Collapse
|
20
|
Zhou C, Gagliardi L, Truhlar DG. State-interaction pair density functional theory for locally avoided crossings of potential energy surfaces in methylamine. Phys Chem Chem Phys 2019; 21:13486-13493. [DOI: 10.1039/c9cp02240f] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
State-interaction pair-density functional theory agrees well with extended MS-CASPT2 in regions of strong state coupling near conical intersections.
Collapse
Affiliation(s)
- Chen Zhou
- Department of Chemistry
- Chemical Theory Center, and Minnesota Supercomputing Institute
- University of Minnesota
- Minneapolis
- USA
| | - Laura Gagliardi
- Department of Chemistry
- Chemical Theory Center, and Minnesota Supercomputing Institute
- University of Minnesota
- Minneapolis
- USA
| | - Donald G. Truhlar
- Department of Chemistry
- Chemical Theory Center, and Minnesota Supercomputing Institute
- University of Minnesota
- Minneapolis
- USA
| |
Collapse
|