1
|
Choudhury A, Santra S, Ghosh D. Understanding the Photoprocesses in Biological Systems: Need for Accurate Multireference Treatment. J Chem Theory Comput 2024; 20:4951-4964. [PMID: 38864715 DOI: 10.1021/acs.jctc.4c00027] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/13/2024]
Abstract
Light-matter interaction is crucial to life itself and revolves around many of the central processes in biology. The need for understanding these photochemical and photophysical processes cannot be overemphasized. Interaction of light with biological systems starts with the absorption of light and subsequent phenomena that occur in the excited states of the system. However, excited states are typically difficult to understand within the mean field approximation of quantum chemical methods. Therefore, suitable multireference methods and methodologies have been developed to understand these phenomena. In this Perspective, we will describe a few methods and methodologies suitable for these descriptions and discuss some persisting difficulties.
Collapse
Affiliation(s)
- Arpan Choudhury
- School of Chemical Sciences, Indian Association for the Cultivation of Science, Jadavpur, Kolkata 700032, India
| | - Supriyo Santra
- School of Chemical Sciences, Indian Association for the Cultivation of Science, Jadavpur, Kolkata 700032, India
| | - Debashree Ghosh
- School of Chemical Sciences, Indian Association for the Cultivation of Science, Jadavpur, Kolkata 700032, India
| |
Collapse
|
2
|
Ozaki Y, Morisawa Y, Tanabe I. ATR-far-ultraviolet spectroscopy: a challenge to new σ chemistry. Chem Soc Rev 2024; 53:1730-1768. [PMID: 38287893 DOI: 10.1039/d3cs00437f] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/31/2024]
Abstract
This review reports the recent progress on ATR-far ultraviolet (FUV) spectroscopy in the condensed phase. ATR-FUV spectroscopy for liquids and solids enables one to explore various topics in physical chemistry, analytical chemistry, nanoscience and technology, materials science, electrochemistry, and organic chemistry. In this review, we put particular emphasis on the three major topics: (1) studies on electronic transitions and structures of various molecules, which one cannot investigate via ordinary UV spectroscopy. The combined use of ATR-FUV spectroscopy and quantum chemical calculations allows for the investigation of various electronic transitions, including σ, n-Rydberg transitions. ATR-FUV spectroscopy may open a new avenue for σ-chemistry. (2) ATR-FUV spectroscopy enables one to measure the first electronic transition of water at approximately 160 nm without peak saturation. Using this band, one can study the electronic structure of water, aqueous solutions, and adsorbed water. (3) ATR-FUV spectroscopy has its own advantages of the ATR method as a surface analysis method. ATR-FUV spectroscopy is a powerful technique for exploring a variety of top surface phenomena (∼50 nm) in adsorbed water, polymers, graphene, organic materials, ionic liquids, and so on. This review briefly describes the principles, characteristics, and instrumentation of ATR-FUV spectroscopy. Next, a detailed description about quantum chemical calculation methods for FUV and UV regions is given. The recent application of ATR-FUV-UV spectroscopy studies on electronic transitions from σ orbitals in various saturated molecules is introduced first, followed by a discussion on the applications of ATR-FUV spectroscopy to studies on water, aqueous solutions, and adsorbed water. Applications of ATR-FUV spectroscopy in the analysis of other materials such as polymers, ionic liquids, inorganic semiconductors, graphene, and carbon nanocomposites are elucidated. In addition, ATR-FUV-UV-vis spectroscopy focusing on electrochemical interfaces is outlined. Finally, FUV-UV-surface plasmon resonance studies are discussed.
Collapse
Affiliation(s)
- Yukihiro Ozaki
- School of Biological and Environmental Sciences, Kwansei Gakuin University, Sanda, Hyogo 669-1330, Japan.
- Toyota Physical and Chemical Research Institute, Nagakute, Aichi 480-1192, Japan
| | - Yusuke Morisawa
- School of Science and Engineering, Kindai University, Higashi-Osaka, Osaka 577-8502, Japan.
| | - Ichiro Tanabe
- Department of Chemistry, School of Science, Rikkyo University, Toshima, Tokyo 171-8501, Japan.
| |
Collapse
|
3
|
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
|
4
|
Komarov K, Park W, Lee S, Huix-Rotllant M, Choi CH. Doubly Tuned Exchange-Correlation Functionals for Mixed-Reference Spin-Flip Time-Dependent Density Functional Theory. J Chem Theory Comput 2023; 19:7671-7684. [PMID: 37844129 DOI: 10.1021/acs.jctc.3c00884] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/18/2023]
Abstract
It is demonstrated that significant accuracy improvements in MRSF-TDDFT can be achieved by introducing two different exchange-correlation (XC) functionals for the reference Kohn-Sham DFT and the response part of the calculations, respectively. Accordingly, two new XC functionals of doubly tuned Coulomb attenuated method-vertical excitation energy (DTCAM-VEE) and DTCAM-AEE were developed on the basis of the "adaptive exact exchange (AEE)" concept in the framework of the Coulomb-attenuating XC functionals. The values by DTCAM-VEE are in excellent agreement with those of Thiel's set [mean absolute errors (MAEs) and the interquartile range (IQR) values of 0.218 and 0.327 eV, respectively]. On the other hand, DTCAM-AEE faithfully reproduced the qualitative aspects of conical intersections (CIs) of trans-butadiene and thymine and the nonadiabatic molecular dynamics (NAMD) simulations on thymine. The latter functional also remarkably exhibited the exact 1/R asymptotic behavior of the charge-transfer state of an ethylene-tetrafluoroethylene dimer and the accurate potential energy surfaces (PESs) along the two torsional angles of retinal protonated Schiff base model with six double bonds (rPSB6). Overall, DTCAM-AEE generally performs well, as its MAE (0.237) and IQR (0.41 eV) are much improved as compared to BH&HLYP. The current idea can also be applied to other XC functionals as well as other variants of linear response theories, opening a new way of developing XC functionals.
Collapse
Affiliation(s)
- Konstantin Komarov
- Center for Quantum Dynamics, Pohang University of Science and Technology, Pohang 37673, South Korea
| | - Woojin Park
- Department of Chemistry, Kyungpook National University, Daegu 41566, South Korea
| | - Seunghoon Lee
- Department of Chemistry, Seoul National University, Seoul, 151-747, South Korea
| | | | - Cheol Ho Choi
- Department of Chemistry, Kyungpook National University, Daegu 41566, South Korea
| |
Collapse
|
5
|
Fdez Galván I, Lindh R. Smooth Things Come in Threes: A Diabatic Surrogate Model for Conical Intersection Optimization. J Chem Theory Comput 2023. [PMID: 37192531 DOI: 10.1021/acs.jctc.3c00389] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/18/2023]
Abstract
The optimization of conical intersection structures is complicated by the nondifferentiability of the adiabatic potential energy surfaces. In this work, we build a pseudodiabatic surrogate model, based on Gaussian process regression, formed by three smooth and differentiable surfaces that can adequately reproduce the adiabatic surfaces. Using this model with the restricted variance optimization method results in a notable decrease of the overall computational effort required to obtain minimum energy crossing points.
Collapse
Affiliation(s)
- Ignacio Fdez Galván
- Department of Chemistry-BMC, Uppsala University, P.O. Box 576, SE-75123 Uppsala, Sweden
| | - Roland Lindh
- Department of Chemistry-BMC, Uppsala University, P.O. Box 576, SE-75123 Uppsala, Sweden
- Uppsala Center for Computational Chemistry (UC3), Uppsala University, P.O. Box 576, SE-75123 Uppsala, Sweden
| |
Collapse
|
6
|
Sirimatayanant S, Andruniów T. Benchmarking two-photon absorption strengths of rhodopsin chromophore models with CC3 and CCSD methodologies: An assessment of popular density functional approximations. J Chem Phys 2023; 158:094106. [PMID: 36889953 DOI: 10.1063/5.0135594] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/11/2023] Open
Abstract
This work presents the investigations of the impact of an increasing electron correlation in the hierarchy of coupled-cluster methods, i.e., CC2, CCSD, and CC3, on two-photon absorption (2PA) strengths for the lowest excited state of the minimal rhodopsin's chromophore model-cis-penta-2,4-dieniminium cation (PSB3). For a larger chromophore's model [4-cis-hepta-2,4,6-trieniminium cation (PSB4)], CC2 and CCSD calculations of 2PA strengths were performed. Additionally, 2PA strengths predicted by some popular density functional theory (DFT) functionals differing in HF exchange contribution were assessed against the reference CC3/CCSD data. For PSB3, the accuracy of 2PA strengths increases in the following order: CC2 < CCSD < CC3, with the CC2 deviation from both higher-level methods exceeding 10% at 6-31+G* basis sets and 2% at aug-cc-pVDZ basis set. However, for PSB4, this trend is reversed and CC2-based 2PA strength is larger than the corresponding CCSD value. Among the DFT functionals investigated, CAM-B3LYP and BHandHLYP provide 2PA strengths in best compliance with reference data, however, with the error approaching an order of magnitude.
Collapse
Affiliation(s)
- Saruti Sirimatayanant
- Institute of Advanced Materials, Department of Chemistry, Wroclaw University of Science and Technology, Wyb. Wyspianskiego 27, 50-370 Wroclaw, Poland
| | - Tadeusz Andruniów
- Institute of Advanced Materials, Department of Chemistry, Wroclaw University of Science and Technology, Wyb. Wyspianskiego 27, 50-370 Wroclaw, Poland
| |
Collapse
|
7
|
Kotaru S, Kähler S, Alessio M, Krylov AI. Magnetic exchange interactions in binuclear and tetranuclear iron(III) complexes described by spin-flip DFT and Heisenberg effective Hamiltonians. J Comput Chem 2023; 44:367-380. [PMID: 35699152 PMCID: PMC10084445 DOI: 10.1002/jcc.26941] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/20/2021] [Revised: 04/07/2022] [Accepted: 04/22/2022] [Indexed: 12/31/2022]
Abstract
Low-energy spectra of single-molecule magnets (SMMs) are often described by Heisenberg Hamiltonians. Within this formalism, exchange interactions between magnetic centers determine the ground-state multiplicity and energy separation between the ground and excited states. In this contribution, we extract exchange coupling constants (J) for a set of iron (III) binuclear and tetranuclear complexes from all-electron calculations using non-collinear spin-flip time-dependent density functional theory (NC-SF-TDDFT). For 12 binuclear complexes with J-values ranging from -6 to -132 cm-1 , our benchmark calculations using the short-range hybrid ωPBEh functional and 6-31G(d,p) basis set agree well with the experimentally derived values (mean absolute error of 4.7 cm-1 ). For the tetranuclear SMMs, the computed J constants are within 6 cm-1 from the experimentally derived values. We explore the range of applicability of the Heisenberg model by analyzing bonding patterns in these Fe(III) complexes using natural orbitals (NO), their occupations, and the number of effectively unpaired electrons. The results illustrate the efficiency of the spin-flip protocol for computing the exchange couplings and the utility of the NO analysis in assessing the validity of effective spin Hamiltonians.
Collapse
Affiliation(s)
- Saikiran Kotaru
- Department of ChemistryUniversity of Southern CaliforniaLos AngelesCalifornia
| | - Sven Kähler
- Department of ChemistryUniversity of Southern CaliforniaLos AngelesCalifornia
| | - Maristella Alessio
- Department of ChemistryUniversity of Southern CaliforniaLos AngelesCalifornia
| | - Anna I. Krylov
- Department of ChemistryUniversity of Southern CaliforniaLos AngelesCalifornia
| |
Collapse
|
8
|
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
|
9
|
Hutcheson A, Paul AC, Myhre RH, Koch H, Høyvik IM. Describing ground and excited state potential energy surfaces for molecular photoswitches using coupled cluster models. J Comput Chem 2021; 42:1419-1429. [PMID: 33973669 DOI: 10.1002/jcc.26553] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/16/2020] [Revised: 03/19/2021] [Accepted: 04/20/2021] [Indexed: 11/08/2022]
Abstract
In this article, we use two extensively studied systems, a retinal model system and azobenzene, to explore the use of coupled cluster models for describing ground and singlet excited state potential energy surfaces of photoswitchable systems. While not being suitable for describing nuclear dynamics of photoisomerization, coupled cluster models have useful attributes, such as the inclusion of dynamical correlation, their black box nature, and the systematic improvement offered by truncation level. Results for the studied systems show that when triple excitations (here through the CC3 model) are included, ground and excited state potential energy surfaces for isomerization paths may reliably be generated, also for states of doubly excited character. For ground state equilibrium cis- and trans-azobenzene, the molecular geometry and basis set is seen to significantly impact the vertical excitation energies for the two lowest excited states. Efficient implementations of coupled cluster models can therefore constitute valuable tools for investigating photoswitchable systems and can be used for preliminary black box studies to gather information before more complicated excited state dynamics approaches are pursued.
Collapse
Affiliation(s)
- Anders Hutcheson
- Department of Chemistry, The Norwegian University of Science and Technology, Trondheim, Norway
| | | | - Rolf H Myhre
- Department of Chemistry, The Norwegian University of Science and Technology, Trondheim, Norway
| | | | - Ida-Marie Høyvik
- Department of Chemistry, The Norwegian University of Science and Technology, Trondheim, Norway
| |
Collapse
|
10
|
Monino E, Loos PF. Spin-Conserved and Spin-Flip Optical Excitations from the Bethe-Salpeter Equation Formalism. J Chem Theory Comput 2021; 17:2852-2867. [PMID: 33724811 PMCID: PMC8154368 DOI: 10.1021/acs.jctc.1c00074] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/01/2023]
Abstract
![]()
Like adiabatic time-dependent
density-functional theory (TD-DFT),
the Bethe–Salpeter equation (BSE) formalism of many-body perturbation
theory, in its static approximation, is “blind” to double
(and higher) excitations, which are ubiquitous, for example, in conjugated
molecules like polyenes. Here, we apply the spin-flip ansatz (which considers the lowest triplet state as the reference configuration
instead of the singlet ground state) to the BSE formalism in order
to access, in particular, double excitations. The present scheme is
based on a spin-unrestricted version of the GW approximation
employed to compute the charged excitations and screened Coulomb potential
required for the BSE calculations. Dynamical corrections to the static
BSE optical excitations are taken into account via an unrestricted
generalization of our recently developed (renormalized) perturbative
treatment. The performance of the present spin-flip BSE formalism
is illustrated by computing excited-state energies of the beryllium
atom, the hydrogen molecule at various bond lengths, and cyclobutadiene
in its rectangular and square-planar geometries.
Collapse
Affiliation(s)
- Enzo Monino
- Laboratoire de Chimie et Physique Quantiques (UMR 5626), Université de Toulouse, CNRS, UPS, Route de Narbonne, 31062 Toulouse, France
| | - Pierre-François Loos
- Laboratoire de Chimie et Physique Quantiques (UMR 5626), Université de Toulouse, CNRS, UPS, Route de Narbonne, 31062 Toulouse, France
| |
Collapse
|
11
|
Gulania S, Kjønstad EF, Stanton JF, Koch H, Krylov AI. Equation-of-motion coupled-cluster method with double electron-attaching operators: Theory, implementation, and benchmarks. J Chem Phys 2021; 154:114115. [PMID: 33752380 DOI: 10.1063/5.0041822] [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/14/2022] Open
Abstract
We report a production-level implementation of the equation-of-motion (EOM) coupled-cluster (CC) method with double electron-attaching (DEA) EOM operators of 2p and 3p1h types, EOM-DEA-CCSD. This ansatz, suitable for treating electronic structure patterns that can be described as two-electrons-in-many orbitals, represents a useful addition to the EOM-CC family of methods. We analyze the performance of EOM-DEA-CCSD for energy differences and molecular properties. By considering reduced quantities, such as state and transition one-particle density matrices, we compare EOM-DEA-CCSD wave functions with wave functions computed by other EOM-CCSD methods. The benchmarks illustrate that EOM-DEA-CCSD is capable of treating diradicals, bond-breaking, and some types of conical intersections.
Collapse
Affiliation(s)
- Sahil Gulania
- Department of Chemistry, University of Southern California, Los Angeles, California 90089, USA
| | - Eirik F Kjønstad
- Department of Chemistry, Norwegian University of Science and Technology, 7491 Trondheim, Norway
| | - John F Stanton
- Quantum Theory Project, Departments of Chemistry and Physics, University of Florida, Gainesville, Florida 32611, USA
| | - Henrik Koch
- Scuola Normale Superiore, Piazza dei Cavaleri 7, 56126 Pisa, Italy
| | - Anna I Krylov
- Department of Chemistry, University of Southern California, Los Angeles, California 90089, USA
| |
Collapse
|
12
|
Kretz B, Egger DA. Accurate Molecular Geometries in Complex Excited-State Potential Energy Surfaces from Time-Dependent Density Functional Theory. J Chem Theory Comput 2021; 17:357-366. [PMID: 33284603 DOI: 10.1021/acs.jctc.0c00858] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/17/2022]
Abstract
The interplay of electronic excitations and structural changes in molecules impacts nonradiative decay and charge transfer in the excited state, thus influencing excited-state lifetimes and photocatalytic reaction rates in optoelectronic and energy devices. To capture such effects requires computational methods providing an accurate description of excited-state potential energy surfaces and geometries. We suggest time-dependent density functional theory using optimally tuned range-separated hybrid (OT-RSH) functionals as an accurate approach to obtain excited-state molecular geometries. We show that OT-RSH provides accurate molecular geometries in excited-state potential energy surfaces that are complex and involve an interplay of local and charge-transfer excitations, for which conventional semilocal and hybrid functionals fail. At the same time, the nonempirical OT-RSH approach maintains the high accuracy of parametrized functionals (e.g., B3LYP) for predicting excited-state geometries of small organic molecules showing valence excited states.
Collapse
Affiliation(s)
- Bernhard Kretz
- Department of Physics, Technical University of Munich, James-Franck-Str. 1, 85748 Garching, Germany
| | - David A Egger
- Department of Physics, Technical University of Munich, James-Franck-Str. 1, 85748 Garching, Germany
| |
Collapse
|
13
|
Qu Z, Ma Y. Variational Multistate Density Functional Theory for a Balanced Treatment of Static and Dynamic Correlations. J Chem Theory Comput 2020; 16:4912-4922. [PMID: 32672966 DOI: 10.1021/acs.jctc.0c00208] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
Here, an approach to variational multistate density functional theory (vMSDFT) is explored. In this approach, the Kohn-Sham orbitals as well as configuration coefficients were simultaneously optimized, thus yielding a full variational minimum. Furthermore, this work also proposes two important improvements on the MSDFT framework. First, a "point-to-point correction" is used to correct the static correlation present in the DFT framework. Therefore, double counting of static correlation in vMSDFT is mitigated. Second, a general form to construct the transition density functional in the vMSDFT framework is proposed, which allows for the properties of vMSDFT wave functions to be standardized to the complete active space self-consistent field properties. The utility of vMSDFT is illustrated on molecular systems of interest including bond breaking, diradicals, excited states, and conical intersections. The numerical results suggest that the accuracy of vMSDFT is in close agreement with the high-level multireference methods.
Collapse
Affiliation(s)
- Zexing Qu
- Institute of Theoretical Chemistry and Laboratory of Theoretical & Computational Chemistry, Jilin University, Changchun 130023, China
| | - Yingjin Ma
- Computer Network Information Center, Chinese Academy of Sciences, Beijing 100190, China.,Center of Scientific Computing Applications & Research, Chinese Academy of Sciences, Beijing 100190, China
| |
Collapse
|
14
|
Song C, Martínez TJ. Reduced scaling extended multi-state CASPT2 (XMS-CASPT2) using supporting subspaces and tensor hyper-contraction. J Chem Phys 2020; 152:234113. [PMID: 32571032 DOI: 10.1063/5.0007417] [Citation(s) in RCA: 12] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/31/2023] Open
Abstract
We present a reduced scaling formulation of the extended multi-state CASPT2 (XMS-CASPT2) method, which is based on our recently developed state-specific CASPT2 (SS-CASPT2) formulation using supporting subspaces and tensor hyper-contraction. By using these two techniques, the off-diagonal elements of the effective Hamiltonian can be computed with only O(N3) operations and O(N2) memory, where N is the number of basis functions. This limits the overall computational scaling to O(N4) operations and O(N2) memory. Thus, excited states can now be obtained at the same reduced (relative to previous algorithms) scaling we achieved for SS-CASPT2. In addition, we also investigate how the energy denominators can be factorized with the Laplace quadrature when some of the denominators are negative, which is critical for excited state calculations. An efficient implementation of the method has been developed using graphical processing units while also exploiting spatial sparsity in tensor operations. We benchmark the accuracy of the new method by comparison to non-THC formulated XMS-CASPT2 for the excited states of various molecules. In our tests, the THC approximation introduces negligible errors (≈0.01 eV) compared to the non-THC reference method. Scaling behavior and computational timings are presented to demonstrate performance. The new method is also interfaced with quantum mechanics/molecular mechanics (QM/MM). In an example study of green fluorescent protein, we show how the XMS-CASPT2 potential energy surfaces and excitation energies are affected by increasing the size of the QM region up to 278 QM atoms with more than 2300 basis functions.
Collapse
Affiliation(s)
- Chenchen Song
- 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
|
15
|
Abstract
This Perspective discusses salient features of the spin-flip approach to strong correlation and describes different methods that sprung from this idea. The spin-flip treatment exploits the different physics of low-spin and high-spin states and is based on the observation that correlation is small for same-spin electrons. By using a well-behaved high-spin state as a reference, one can access problematic low-spin states by deploying the same formal tools as in the excited-state treatments (i.e., linear response, propagator, or equation-of-motion theories). The Perspective reviews applications of this strategy within wave function and density functional theory frameworks as well as the extensions for molecular properties and spectroscopy. The utility of spin-flip methods is illustrated by examples. Limitations and proposed future directions are also discussed.
Collapse
Affiliation(s)
- David Casanova
- Donostia International Physics Center (DIPC), 20018 Donostia, Euskadi, Spain. and IKERBASQUE, Basque Foundation for Science, 48013 Bilbao, Euskadi, Spain
| | - Anna I Krylov
- Department of Chemistry, University of Southern California, Los Angeles, California 90089, USA
| |
Collapse
|
16
|
Mato J, Gordon MS. Analytic non-adiabatic couplings for the spin-flip ORMAS method. Phys Chem Chem Phys 2020; 22:1475-1484. [DOI: 10.1039/c9cp05849d] [Citation(s) in RCA: 9] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Analytic non-adiabatic coupling matrix elements (NACME) are derived and implemented for the spin-flip occupation restricted multiple active space configuration interaction (SF-ORMAS-CI) method.
Collapse
Affiliation(s)
- Joani Mato
- Department of Chemistry and Ames Laboratory
- Iowa State University
- Ames
- USA
| | - Mark S. Gordon
- Department of Chemistry and Ames Laboratory
- Iowa State University
- Ames
- USA
| |
Collapse
|
17
|
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
|
18
|
Loos PF, Boggio-Pasqua M, Scemama A, Caffarel M, Jacquemin D. Reference Energies for Double Excitations. J Chem Theory Comput 2019; 15:1939-1956. [DOI: 10.1021/acs.jctc.8b01205] [Citation(s) in RCA: 79] [Impact Index Per Article: 15.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/19/2022]
Affiliation(s)
- Pierre-François Loos
- Laboratoire de Chimie et Physique Quantiques (UMR 5626), Université de Toulouse, CNRS, UPS, 31062 Toulouse, France
| | - Martial Boggio-Pasqua
- Laboratoire de Chimie et Physique Quantiques (UMR 5626), Université de Toulouse, CNRS, UPS, 31062 Toulouse, France
| | - Anthony Scemama
- Laboratoire de Chimie et Physique Quantiques (UMR 5626), Université de Toulouse, CNRS, UPS, 31062 Toulouse, France
| | - Michel Caffarel
- Laboratoire de Chimie et Physique Quantiques (UMR 5626), Université de Toulouse, CNRS, UPS, 31062 Toulouse, France
| | - Denis Jacquemin
- Laboratoire CEISAM (UMR 6230), CNRS, Université de Nantes, 44399 Cedex 3 Nantes, France
| |
Collapse
|
19
|
Segalina A, Francés-Monerris A, Pastore M, Leininger T, Evangelisti S, Monari A. Conical intersection properties unraveled by the position spread tensor. Theor Chem Acc 2018. [DOI: 10.1007/s00214-018-2377-y] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/01/2022]
|
20
|
Sen S, Schapiro I. A comprehensive benchmark of the XMS-CASPT2 method for the photochemistry of a retinal chromophore model. Mol Phys 2018. [DOI: 10.1080/00268976.2018.1501112] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/03/2023]
Affiliation(s)
- Saumik Sen
- Fritz Haber Center for Molecular Dynamics, Institute of Chemistry, The Hebrew University of Jerusalem, Jerusalem, Israel
| | - Igor Schapiro
- Fritz Haber Center for Molecular Dynamics, Institute of Chemistry, The Hebrew University of Jerusalem, Jerusalem, Israel
| |
Collapse
|
21
|
Wolański Ł, Grabarek D, Andruniów T. Is the choice of a standard zeroth-order hamiltonian in CASPT2 ansatz optimal in calculations of excitation energies in protonated and unprotonated schiff bases of retinal? J Comput Chem 2018; 39:1470-1480. [PMID: 29635695 DOI: 10.1002/jcc.25217] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/29/2017] [Revised: 03/09/2018] [Accepted: 03/10/2018] [Indexed: 11/07/2022]
Abstract
To account for systematic error of CASPT2 method empirical modification of the zeroth-order Hamiltonian with Ionization Potential-Electron Affinity (IPEA) shift was introduced. The optimized IPEA value (0.25 a.u.), called standard IPEA (S-IPEA), was recommended but due to its unsatisfactory performance in multiple metallic and organic compounds it has been questioned lately as a general parameter working properly for all molecules under CASPT2 study. As we are interested in Schiff bases of retinal, an important question emerging from this conflict of choice, to use or not to use S-IPEA, is whether the introduction of the modified zeroth-order Hamiltonian into CASPT2 ansatz does really improve their energetics. To achieve this goal, we assessed an impact of the IPEA shift value, in a range of 0-0.35 a.u., on vertical excitation energies to low-lying singlet states of two protonated (RPSBs) and two unprotonated (RSBs) Schiff bases of retinal for which experimental data in gas phase are available. In addition, an effect of geometry, basis set, and active space on computed VEEs is also reported. We find, that for these systems, the choice of S-IPEA significantly overestimates both S0 →S1 and S0 →S2 energies and the best theoretical estimate, in reference to the experimental data, is provided with either unmodified zeroth-order Hamiltonian or small value of the IPEA shift in a range of 0.05-0.15 a.u., depending on active space and basis set size, equilibrium geometry, and character of the excited state. © 2018 Wiley Periodicals, Inc.
Collapse
Affiliation(s)
- Łukasz Wolański
- Advanced Materials Engineering and Modelling Group, Wroclaw University of Science and Technology, Wyb. Wyspianskiego 27, Wroclaw, 50-370, Poland
| | - Dawid Grabarek
- Advanced Materials Engineering and Modelling Group, Wroclaw University of Science and Technology, Wyb. Wyspianskiego 27, Wroclaw, 50-370, Poland
| | - Tadeusz Andruniów
- Advanced Materials Engineering and Modelling Group, Wroclaw University of Science and Technology, Wyb. Wyspianskiego 27, Wroclaw, 50-370, Poland
| |
Collapse
|
22
|
Filatov M, Min SK, Kim KS. Direct Nonadiabatic Dynamics by Mixed Quantum-Classical Formalism Connected with Ensemble Density Functional Theory Method: Application to trans-Penta-2,4-dieniminium Cation. J Chem Theory Comput 2018; 14:4499-4512. [DOI: 10.1021/acs.jctc.8b00217] [Citation(s) in RCA: 23] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/28/2022]
Affiliation(s)
- Michael Filatov
- Department of Chemistry, School of Natural Sciences, Ulsan National Institute of Science and Technology (UNIST), Ulsan 44919, Korea
| | - Seung Kyu Min
- Department of Chemistry, School of Natural Sciences, Ulsan National Institute of Science and Technology (UNIST), Ulsan 44919, Korea
| | - Kwang S. Kim
- Department of Chemistry, School of Natural Sciences, Ulsan National Institute of Science and Technology (UNIST), Ulsan 44919, Korea
| |
Collapse
|
23
|
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
|
24
|
Filatov M, Liu F, Martínez TJ. Analytical derivatives of the individual state energies in ensemble density functional theory method. I. General formalism. J Chem Phys 2018; 147:034113. [PMID: 28734302 DOI: 10.1063/1.4994542] [Citation(s) in RCA: 46] [Impact Index Per Article: 7.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Abstract
The state-averaged (SA) spin restricted ensemble referenced Kohn-Sham (REKS) method and its state interaction (SI) extension, SI-SA-REKS, enable one to describe correctly the shape of the ground and excited potential energy surfaces of molecules undergoing bond breaking/bond formation reactions including features such as conical intersections crucial for theoretical modeling of non-adiabatic reactions. Until recently, application of the SA-REKS and SI-SA-REKS methods to modeling the dynamics of such reactions was obstructed due to the lack of the analytical energy derivatives. In this work, the analytical derivatives of the individual SA-REKS and SI-SA-REKS energies are derived. The final analytic gradient expressions are formulated entirely in terms of traces of matrix products and are presented in the form convenient for implementation in the traditional quantum chemical codes employing basis set expansions of the molecular orbitals. The implementation and benchmarking of the derived formalism will be described in a subsequent article of this series.
Collapse
Affiliation(s)
- Michael Filatov
- Department of Chemistry, School of Natural Sciences, Ulsan National Institute of Science and Technology (UNIST), Ulsan 689-798, South Korea
| | - Fang Liu
- Department of Chemistry and the PULSE Institute, Stanford University, Stanford, California 94305, USA and SLAC National Accelerator Laboratory, Menlo Park, California 94025, USA
| | - Todd J Martínez
- Department of Chemistry and the PULSE Institute, Stanford University, Stanford, California 94305, USA and SLAC National Accelerator Laboratory, Menlo Park, California 94025, USA
| |
Collapse
|
25
|
Abstract
Singlet fission is a photophysical reaction in which a singlet excited electronic state splits into two spin-triplet states. Singlet fission was discovered more than 50 years ago, but the interest in this process has gained a lot of momentum in the past decade due to its potential as a way to boost solar cell efficiencies. This review presents and discusses the most recent advances with respect to the theoretical and computational studies on the singlet fission phenomenon. The work revisits important aspects regarding electronic states involved in the process, the evaluation of fission rates and interstate couplings, the study of the excited state dynamics in singlet fission, and the advances in the design and characterization of singlet fission compounds and materials such as molecular dimers, polymers, or extended structures. Finally, the review tries to pinpoint some aspects that need further improvement and proposes future lines of research for theoretical and computational chemists and physicists in order to further push the understanding and applicability of singlet fission.
Collapse
Affiliation(s)
- David Casanova
- Kimika Fakultatea , Euskal Herriko Unibertsitatea (UPV/EHU) and Donostia International Physics Center (DIPC) , P.K. 1072, 20080 Donostia , Euskadi, Spain.,IKERBASQUE, Basque, Foundation for Science , 48013 Bilbao , Euskadi, Spain
| |
Collapse
|
26
|
Faraji S, Matsika S, Krylov AI. Calculations of non-adiabatic couplings within equation-of-motion coupled-cluster framework: Theory, implementation, and validation against multi-reference methods. J Chem Phys 2018; 148:044103. [DOI: 10.1063/1.5009433] [Citation(s) in RCA: 31] [Impact Index Per Article: 5.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/31/2022] Open
Affiliation(s)
- Shirin Faraji
- Department of Chemistry, University of Southern California, Los Angeles, California 90089, USA
- Zernike Institute for Advanced Materials, Groningen, The Netherlands
| | - Spiridoula Matsika
- Department of Chemistry, Temple University, Philadelphia, Pennsylvania 19122, USA
| | - Anna I. Krylov
- Department of Chemistry, University of Southern California, Los Angeles, California 90089, USA
| |
Collapse
|
27
|
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
|
28
|
Farag MH, Jansen TLC, Knoester J. The origin of absorptive features in the two-dimensional electronic spectra of rhodopsin. Phys Chem Chem Phys 2018; 20:12746-12754. [DOI: 10.1039/c8cp00638e] [Citation(s) in RCA: 12] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
A three-state three-mode model Hamiltonian reveals the origin of the absorptive features in the two-dimensional electronic spectra of rhodopsin.
Collapse
Affiliation(s)
- Marwa H. Farag
- University of Groningen
- Zernike Institute for Advanced Materials
- 9747 AG Groningen
- The Netherlands
| | - Thomas L. C. Jansen
- University of Groningen
- Zernike Institute for Advanced Materials
- 9747 AG Groningen
- The Netherlands
| | - Jasper Knoester
- University of Groningen
- Zernike Institute for Advanced Materials
- 9747 AG Groningen
- The Netherlands
| |
Collapse
|
29
|
Hollas D, Šištík L, Hohenstein EG, Martínez TJ, Slavíček P. Nonadiabatic Ab Initio Molecular Dynamics with the Floating Occupation Molecular Orbital-Complete Active Space Configuration Interaction Method. J Chem Theory Comput 2017; 14:339-350. [DOI: 10.1021/acs.jctc.7b00958] [Citation(s) in RCA: 42] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Daniel Hollas
- Department
of Physical Chemistry, University of Chemistry and Technology, Prague, Technická 5, 16628 Prague 6, Czech Republic
| | - Lukáš Šištík
- Department
of Physical Chemistry, University of Chemistry and Technology, Prague, Technická 5, 16628 Prague 6, Czech Republic
| | - Edward G. Hohenstein
- Department
of Chemistry and Biochemistry, The City College of New York, New York, New York 10031, United States
- PhD
Program in Chemistry, The Graduate Center of the City University of New York, New York, New York 10016, United States
| | - Todd J. Martínez
- Department
of Chemistry and the PULSE Institute, Stanford University, Stanford, California 94305, United States
- SLAC National Accelerator Laboratory, Menlo Park, California 94025, United States
| | - Petr Slavíček
- Department
of Physical Chemistry, University of Chemistry and Technology, Prague, Technická 5, 16628 Prague 6, Czech Republic
- J. Heyrovský Institute of Physical Chemistry, Dolejškova 3, 18223 Prague 8, Czech Republic
| |
Collapse
|
30
|
Gozem S, Luk HL, Schapiro I, Olivucci M. Theory and Simulation of the Ultrafast Double-Bond Isomerization of Biological Chromophores. Chem Rev 2017; 117:13502-13565. [DOI: 10.1021/acs.chemrev.7b00177] [Citation(s) in RCA: 175] [Impact Index Per Article: 25.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/03/2023]
Affiliation(s)
- Samer Gozem
- Department
of Chemistry, Georgia State University, Atlanta, Georgia 30302, United States
| | - Hoi Ling Luk
- Chemistry
Department, Bowling Green State University, Overman Hall, Bowling Green, Ohio 43403, United States
| | - Igor Schapiro
- Fritz
Haber Center for Molecular Dynamics, Institute of Chemistry, The Hebrew University of Jerusalem, Jerusalem 91904, Israel
| | - Massimo Olivucci
- Chemistry
Department, Bowling Green State University, Overman Hall, Bowling Green, Ohio 43403, United States
- Dipartimento
di Biotecnologie, Chimica e Farmacia, Università di Siena, via A. Moro
2, 53100 Siena, Italy
| |
Collapse
|
31
|
Kjønstad EF, Myhre RH, Martínez TJ, Koch H. Crossing conditions in coupled cluster theory. J Chem Phys 2017; 147:164105. [DOI: 10.1063/1.4998724] [Citation(s) in RCA: 32] [Impact Index Per Article: 4.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Eirik F. Kjønstad
- Department of Chemistry, Norwegian University of Science and Technology, 7491 Trondheim, Norway
- Department of Chemistry and the PULSE Institute, Stanford University, Stanford, California 94305, USA
| | - Rolf H. Myhre
- Department of Chemistry, Norwegian University of Science and Technology, 7491 Trondheim, Norway
- 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
| | - Henrik Koch
- Department of Chemistry, Norwegian University of Science and Technology, 7491 Trondheim, Norway
- Department of Chemistry and the PULSE Institute, Stanford University, Stanford, California 94305, USA
| |
Collapse
|
32
|
Lefrancois D, Tuna D, Martínez TJ, Dreuw A. The Spin-Flip Variant of the Algebraic-Diagrammatic Construction Yields the Correct Topology of S1/S0 Conical Intersections. J Chem Theory Comput 2017; 13:4436-4441. [DOI: 10.1021/acs.jctc.7b00634] [Citation(s) in RCA: 22] [Impact Index Per Article: 3.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Daniel Lefrancois
- Interdisciplinary
Center for Scientific Computing, Ruprecht-Karls University, Im Neuenheimer
Feld 205, 69120 Heidelberg, Germany
| | - Deniz Tuna
- Max-Planck-Institut für Kohlenforschung, 45470 Mülheim an der Ruhr, Germany
| | - Todd J. Martínez
- Department
of Chemistry and The PULSE Institute, Stanford University, Stanford, California 94305, United States
- SLAC National Accelerator Laboratory, Menlo Park, California 94025, United States
| | - Andreas Dreuw
- Interdisciplinary
Center for Scientific Computing, Ruprecht-Karls University, Im Neuenheimer
Feld 205, 69120 Heidelberg, Germany
| |
Collapse
|
33
|
Szefczyk B, Grabarek D, Walczak E, Andruniów T. Excited-state minima and emission energies of retinal chromophore analogues: Performance of CASSCF and CC2 methods as compared with CASPT2. J Comput Chem 2017; 38:1799-1810. [PMID: 28512740 DOI: 10.1002/jcc.24821] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/06/2017] [Revised: 04/10/2017] [Accepted: 04/13/2017] [Indexed: 11/08/2022]
Abstract
This study provides gas-phase S1 excited-state geometries along with emission and adiabatic energies for methylated/demethylated and ring-locked analogues of protonated Schiff base retinal models comprising system of five conjugated double bonds (PSB5), using second order multiconfiguration perturbation theory (CASPT2). CASPT2 results serve as reference data to assess the performance of CC2 (second-order approximate coupled cluster singles and doubles) and a commonly used CASSCF/CASPT2 protocol, that is, complete active space self-consistent field (CASSCF) geometry optimization followed by CASPT2 energy calculation. We find that the CASSCF methodology fails to locate planar S1 minimum energy structures for four out of five investigated planar models in contrast to CC2 and CASPT2 methods. However, for those which were found: one planar and two twisted minima, there is an excellent agreement between CASSCF and CASPT2 results in terms of geometrical parameters, one-electron properties, as well as emission and adiabatic energies. CC2 performs well for in-plane S1 minima and their spectroscopic and electronic properties. However, this picture deteriorates for twisted minima. As expected, the CC2 description of the S2 electronic state, with strong multireference and significant double excitation character, is very poor, exhibiting errors in transition energies exceeding 1 eV. They may be substantially diminished by recalculating transition energies with CASPT2 method. Our work shows that CASSCF/CASPT2 and CC2 shortcomings may influence gas-phase retinal analogues' excited state description in a dramatic way. © 2017 Wiley Periodicals, Inc.
Collapse
Affiliation(s)
- Borys Szefczyk
- Advanced Materials Engineering and Modelling Group, Department of Chemistry, Wroclaw University of Science and Technology, Wyb. Wyspianskiego 27, Wroclaw, 50-370, Poland
| | - Dawid Grabarek
- Advanced Materials Engineering and Modelling Group, Department of Chemistry, Wroclaw University of Science and Technology, Wyb. Wyspianskiego 27, Wroclaw, 50-370, Poland
| | - Elżbieta Walczak
- Advanced Materials Engineering and Modelling Group, Department of Chemistry, Wroclaw University of Science and Technology, Wyb. Wyspianskiego 27, Wroclaw, 50-370, Poland
| | - Tadeusz Andruniów
- Advanced Materials Engineering and Modelling Group, Department of Chemistry, Wroclaw University of Science and Technology, Wyb. Wyspianskiego 27, Wroclaw, 50-370, Poland
| |
Collapse
|
34
|
Varsano D, Caprasecca S, Coccia E. Theoretical description of protein field effects on electronic excitations of biological chromophores. JOURNAL OF PHYSICS. CONDENSED MATTER : AN INSTITUTE OF PHYSICS JOURNAL 2017; 29:013002. [PMID: 27830666 DOI: 10.1088/0953-8984/29/1/013002] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/06/2023]
Abstract
Photoinitiated phenomena play a crucial role in many living organisms. Plants, algae, and bacteria absorb sunlight to perform photosynthesis, and convert water and carbon dioxide into molecular oxygen and carbohydrates, thus forming the basis for life on Earth. The vision of vertebrates is accomplished in the eye by a protein called rhodopsin, which upon photon absorption performs an ultrafast isomerisation of the retinal chromophore, triggering the signal cascade. Many other biological functions start with the photoexcitation of a protein-embedded pigment, followed by complex processes comprising, for example, electron or excitation energy transfer in photosynthetic complexes. The optical properties of chromophores in living systems are strongly dependent on the interaction with the surrounding environment (nearby protein residues, membrane, water), and the complexity of such interplay is, in most cases, at the origin of the functional diversity of the photoactive proteins. The specific interactions with the environment often lead to a significant shift of the chromophore excitation energies, compared with their absorption in solution or gas phase. The investigation of the optical response of chromophores is generally not straightforward, from both experimental and theoretical standpoints; this is due to the difficulty in understanding diverse behaviours and effects, occurring at different scales, with a single technique. In particular, the role played by ab initio calculations in assisting and guiding experiments, as well as in understanding the physics of photoactive proteins, is fundamental. At the same time, owing to the large size of the systems, more approximate strategies which take into account the environmental effects on the absorption spectra are also of paramount importance. Here we review the recent advances in the first-principle description of electronic and optical properties of biological chromophores embedded in a protein environment. We show their applications on paradigmatic systems, such as the light-harvesting complexes, rhodopsin and green fluorescent protein, emphasising the theoretical frameworks which are of common use in solid state physics, and emerging as promising tools for biomolecular systems.
Collapse
Affiliation(s)
- Daniele Varsano
- S3 Center, CNR Institute of Nanoscience, Via Campi 213/A, 41125 Modena, Italy
| | | | | |
Collapse
|
35
|
Farag MH, Jansen TLC, Knoester J. Probing the Interstate Coupling near a Conical Intersection by Optical Spectroscopy. J Phys Chem Lett 2016; 7:3328-3334. [PMID: 27509384 DOI: 10.1021/acs.jpclett.6b01463] [Citation(s) in RCA: 16] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/06/2023]
Abstract
Conical intersections are points where adiabatic potential energy surfaces cross. The interstate coupling between the potential energy surfaces plays a crucial role in many processes associated with conical intersections. Still no method exists to measure this coupling driving the chemical reactions between the potential energy surfaces involved. In this Letter, using a generic model for photoisomerization, we propose a novel experimental approach to estimate the coupling that mixes the electronic states near a conical intersection. The approach is based on analyzing the vibrational wavepacket of the reactant in the adiabatic ground and excited electronic states. The nuclear wavepacket dynamics are extracted from linear absorption and two-dimensional electronic spectroscopy. Comparing the frequencies of the coupling mode in the adiabatic ground and excited states from models with and without coupling between the potential energy surfaces suggests an experimental tool to determine the interstate coupling.
Collapse
Affiliation(s)
- Marwa H Farag
- Zernike Institute for Advanced Materials, University of Groningen , Nijenborgh 4, 9747 AG Groningen, The Netherlands
| | - Thomas L C Jansen
- Zernike Institute for Advanced Materials, University of Groningen , Nijenborgh 4, 9747 AG Groningen, The Netherlands
| | - Jasper Knoester
- Zernike Institute for Advanced Materials, University of Groningen , Nijenborgh 4, 9747 AG Groningen, The Netherlands
| |
Collapse
|
36
|
Zen A, Coccia E, Gozem S, Olivucci M, Guidoni L. Quantum Monte Carlo Treatment of the Charge Transfer and Diradical Electronic Character in a Retinal Chromophore Minimal Model. J Chem Theory Comput 2016; 11:992-1005. [PMID: 25821414 PMCID: PMC4357234 DOI: 10.1021/ct501122z] [Citation(s) in RCA: 36] [Impact Index Per Article: 4.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/12/2014] [Indexed: 01/22/2023]
Abstract
![]()
The
penta-2,4-dieniminium cation (PSB3) displays similar ground
state and first excited state potential energy features as those of
the retinal protonated Schiff base (RPSB) chromophore in rhodopsin.
Recently, PSB3 has been used to benchmark several electronic structure
methods, including highly correlated multireference wave function
approaches, highlighting the necessity to accurately describe the
electronic correlation in order to obtain reliable properties even
along the ground state (thermal) isomerization paths. In this work,
we apply two quantum Monte Carlo approaches, the variational Monte
Carlo and the lattice regularized diffusion Monte Carlo, to study
the energetics and electronic properties of PSB3 along representative
minimum energy paths and scans related to its thermal cis–trans isomerization. Quantum Monte Carlo
is used in combination with the Jastrow antisymmetrized geminal power
ansatz, which guarantees an accurate and balanced description of the
static electronic correlation thanks to the multiconfigurational nature
of the antisymmetrized geminal power term, and of the dynamical correlation,
due to the presence of the Jastrow factor explicitly depending on
electron–electron distances. Along the two ground state isomerization
minimum energy paths of PSB3, CASSCF calculations yield wave functions
having either charge transfer or diradical character in proximity
of the two transition state configurations. Here, we observe that
at the quantum Monte Carlo level of theory, only the transition state
with charge transfer character can be located. The conical intersection,
which becomes highly sloped, is observed only if the path connecting
the two original CASSCF transition states is extended beyond the diradical
one, namely by increasing the bond-length-alternation (BLA). These
findings are in good agreement with the results obtained by MRCISD+Q
calculations, and they demonstrate the importance of having an accurate
description of the static and dynamical correlation when studying
isomerization and transition states of conjugated systems.
Collapse
|
37
|
Segarra-Martí J, Garavelli M, Aquilante F. Multiconfigurational Second-Order Perturbation Theory with Frozen Natural Orbitals Extended to the Treatment of Photochemical Problems. J Chem Theory Comput 2016; 11:3772-84. [PMID: 26574459 DOI: 10.1021/acs.jctc.5b00479] [Citation(s) in RCA: 34] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
A new flavor of the frozen natural orbital complete active space second-order perturbation theory method (FNO-CASPT2, Aquilante et al., J. Chem. Phys. 131, 034113) is proposed herein. In this new implementation, the virtual space in Cholesky decomposition-based CASPT2 computations (CD-CASPT2) is truncated by excluding those orbitals that contribute the least toward preserving a predefined value of the trace of an approximate density matrix, as that represents a measure of the amount of dynamic correlation retained in the model. In this way, the amount of correlation included is practically constant at all nuclear arrangements, thus allowing for the computation of smooth electronic states surfaces and energy gradients-essential requirements for theoretical studies in photochemistry. The method has been benchmarked for a series of relevant biochromophores for which large speed-ups have been recorded while retaining the accuracy achieved in the corresponding CD-CASPT2 calculations. Both vertical excitation energies and gradient calculations have been carried out to establish general guidelines as to how much correlation needs to be retained in the calculation for the results to be consistent with the CD-CASPT2 findings. Our results feature errors within a tenth of an eV for the most difficult cases and have been validated to be used for gradient computations where an up to 3-fold speed-up is observed depending on the size of the system and the basis set employed.
Collapse
Affiliation(s)
- Javier Segarra-Martí
- Dipartimento di Chimica "G. Ciamician", Università di Bologna , Via Selmi 2, IT-40126 Bologna, Italy
| | - Marco Garavelli
- Dipartimento di Chimica "G. Ciamician", Università di Bologna , Via Selmi 2, IT-40126 Bologna, Italy.,Université de Lyon, CNRS , Institut de Chimie de Lyon, École Normale Supérieure de Lyon, 46 Allée d'Italie, F-69364 Lyon Cedex 07, France
| | - Francesco Aquilante
- Dipartimento di Chimica "G. Ciamician", Università di Bologna , Via Selmi 2, IT-40126 Bologna, Italy
| |
Collapse
|
38
|
Csehi A, Halász GJ, Cederbaum LS, Vibók Á. Tracking the photodissociation probability of D2+ induced by linearly chirped laser pulses. J Chem Phys 2016; 144:074309. [DOI: 10.1063/1.4941847] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
39
|
Saini P, Banerjee M, Chattopadhyay A. Computational Investigation of the Photochemical Reaction Path of Some Synthesized and Experimentally Analyzed Small-Chain Conjugated Nitrones. J Phys Chem A 2016; 120:396-406. [PMID: 26697926 DOI: 10.1021/acs.jpca.5b11069] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Praveen Saini
- Department
of Chemistry, Birla Institute of Technology and Science (BITS), Pilani−K.K.
Birla Goa Campus, Zuarinagar, Goa 403 726, India
| | - Mainak Banerjee
- Department
of Chemistry, Birla Institute of Technology and Science (BITS), Pilani−K.K.
Birla Goa Campus, Zuarinagar, Goa 403 726, India
| | - Anjan Chattopadhyay
- Department
of Chemistry, Birla Institute of Technology and Science (BITS), Pilani−K.K.
Birla Goa Campus, Zuarinagar, Goa 403 726, India
| |
Collapse
|
40
|
Casanova D, Krylov AI. Quantifying local exciton, charge resonance, and multiexciton character in correlated wave functions of multichromophoric systems. J Chem Phys 2016; 144:014102. [DOI: 10.1063/1.4939222] [Citation(s) in RCA: 27] [Impact Index Per Article: 3.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/12/2023] Open
Affiliation(s)
- David Casanova
- Kimika Fakultatea, Euskal Herriko Unibersitatea (UPV/EHU) and Donostia International Physics Center (DIPC), P.K. 1072, 20018 Donostia, Spain
- IKERBASQUE, Basque Foundation for Science, 48013 Bilbao, Spain
| | - Anna I. Krylov
- Department of Chemistry, University of Southern California, Los Angeles, California 90089-0482, USA
| |
Collapse
|
41
|
Tuna D, Lefrancois D, Wolański Ł, Gozem S, Schapiro I, Andruniów T, Dreuw A, Olivucci M. Assessment of Approximate Coupled-Cluster and Algebraic-Diagrammatic-Construction Methods for Ground- and Excited-State Reaction Paths and the Conical-Intersection Seam of a Retinal-Chromophore Model. J Chem Theory Comput 2015; 11:5758-81. [PMID: 26642989 DOI: 10.1021/acs.jctc.5b00022] [Citation(s) in RCA: 107] [Impact Index Per Article: 11.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/16/2022]
Abstract
As a minimal model of the chromophore of rhodopsin proteins, the penta-2,4-dieniminium cation (PSB3) poses a challenging test system for the assessment of electronic-structure methods for the exploration of ground- and excited-state potential-energy surfaces, the topography of conical intersections, and the dimensionality (topology) of the branching space. Herein, we report on the performance of the approximate linear-response coupled-cluster method of second order (CC2) and the algebraic-diagrammatic-construction scheme of the polarization propagator of second and third orders (ADC(2) and ADC(3)). For the ADC(2) method, we considered both the strict and extended variants (ADC(2)-s and ADC(2)-x). For both CC2 and ADC methods, we also tested the spin-component-scaled (SCS) and spin-opposite-scaled (SOS) variants. We have explored several ground- and excited-state reaction paths, a circular path centered around the S1/S0 surface crossing, and a 2D scan of the potential-energy surfaces along the branching space. We find that the CC2 and ADC methods yield a different dimensionality of the intersection space. While the ADC methods yield a linear intersection topology, we find a conical intersection topology for the CC2 method. We present computational evidence showing that the linear-response CC2 method yields a surface crossing between the reference state and the first response state featuring characteristics that are expected for a true conical intersection. Finally, we test the performance of these methods for the approximate geometry optimization of the S1/S0 minimum-energy conical intersection and compare the geometries with available data from multireference methods. The present study provides new insight into the performance of linear-response CC2 and polarization-propagator ADC methods for molecular electronic spectroscopy and applications in computational photochemistry.
Collapse
Affiliation(s)
- Deniz Tuna
- Max-Planck-Institut für Kohlenforschung , 45470 Mülheim an der Ruhr, Germany
| | - Daniel Lefrancois
- Interdisciplinary Center for Scientific Computing, University of Heidelberg , 69120 Heidelberg, Germany
| | - Łukasz Wolański
- Department of Chemistry, Wrocław University of Technology , 50370 Wrocław, Poland
| | - Samer Gozem
- Department of Chemistry, University of Southern California , Los Angeles, California 90089, United States
| | - Igor Schapiro
- Institut de Physique et Chimie des Matériaux de Strasbourg & Labex NIE, Université de Strasbourg, CNRS UMR 7504 , Strasbourg 67034, France
| | - Tadeusz Andruniów
- Department of Chemistry, Wrocław University of Technology , 50370 Wrocław, Poland
| | - Andreas Dreuw
- Interdisciplinary Center for Scientific Computing, University of Heidelberg , 69120 Heidelberg, Germany
| | - Massimo Olivucci
- Department of Chemistry, Bowling Green State University , Bowling Green, Ohio 43402, United States.,Dipartimento di Biotecnologie, Chimica e Farmacia, Universitá de Siena , 53100 Siena, Italy
| |
Collapse
|
42
|
Bernstein J, Landau A, Zemel E, Kolodney E. Tin-carbon clusters and the onset of microscopic level immiscibility: Experimental and computational study. J Chem Phys 2015; 143:114307. [PMID: 26395705 DOI: 10.1063/1.4930193] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/14/2022] Open
Abstract
We report the experimental observation and computational analysis of the binary tin-carbon gas phase species. These novel ionic compounds are generated by impact of C60(-) anions on a clean tin target at some kiloelectronvolts kinetic energies. Positive Sn(m)C(n)(+) (m = 1-12, 1 ≤ n ≤ 8) ions were detected mass spectrometrically following ejection from the surface. Impact induced shattering of the C60(-) ion followed by sub-surface penetration of the resulting atomic carbon flux forces efficient mixing between target and projectile atoms even though the two elements (Sn/C) are completely immiscible in the bulk. This approach of C60(-) ion beam induced synthesis can be considered as an effective way for producing novel metal-carbon species of the so-called non-carbide forming elements, thus exploring the possible onset of molecular level miscibility in these systems. Sn2C2(+) was found to be the most abundant carbide cluster ion. Its instantaneous formation kinetics and its measured kinetic energy distribution while exiting the surface demonstrate a single impact formation/emission event (on the sub-ps time scale). Optimal geometries were calculated for both neutral and positively charged species using Born-Oppenheimer molecular dynamics for identifying global minima, followed by density functional theory (DFT) structure optimization and energy calculations at the coupled cluster singles, doubles and perturbative triples [CCSD(T)] level. The calculated structures reflect two distinct binding tendencies. The carbon rich species exhibit polyynic/cummulenic nature (tin end capped carbon chains) while the more stoichiometrically balanced species have larger contributions of metal-metal bonding, sometimes resulting in distinct tin and carbon moieties attached to each other (segregated structures). The Sn2C(n) (n = 3-8) and Sn2C(n)(+) (n = 2-8) are polyynic/cummulenic while all neutral Sn(m)C(n) structures (m = 3-4) could be described as small tin clusters (dimer, trimer, and tetramer, correspondingly) attached to a nearly linear carbon chain. For example, the 1:1 (Sn:C) Sn3C3 and Sn4C4 clusters are composed of all-tin triangle and rhombus, correspondingly, with a short carbon chain (C3, C4) attached on top. The cationic Sn3C(n)(+) (n = 1-5) and Sn4C(n)(+) (n = 1-4) species exhibit various intermediate geometries. Structure calculations at the CCSD(T) level are essential since the segregation effect is not as easily evident based on the most stable structures calculated by DFT alone. Dependences of bond energies (per atom) reflect the evolution of the segregation effect. The mass spectral abundances could be reasonably rationalized in terms of calculated stabilities of the cluster ions with respect to various dissociation channels.
Collapse
Affiliation(s)
- J Bernstein
- Schulich Faculty of Chemistry, Technion, Haifa 3200008, Israel
| | - A Landau
- Schulich Faculty of Chemistry, Technion, Haifa 3200008, Israel
| | - E Zemel
- Schulich Faculty of Chemistry, Technion, Haifa 3200008, Israel
| | - E Kolodney
- Schulich Faculty of Chemistry, Technion, Haifa 3200008, Israel
| |
Collapse
|
43
|
Liu L, Cui G, Fang WH. Excited States and Photochemistry of Chromophores in the Photoactive Proteins Explored by the Combined Quantum Mechanical and Molecular Mechanical Calculations. ADVANCES IN PROTEIN CHEMISTRY AND STRUCTURAL BIOLOGY 2015; 100:255-84. [PMID: 26415847 DOI: 10.1016/bs.apcsb.2015.07.002] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/13/2022]
Abstract
A photoactive protein usually contains a unique chromophore that is responsible for the initial photoresponse and functions of the photoactive protein are determined by the interaction between the chromophore and its protein surroundings. The combined quantum mechanical and molecular mechanical (QM/MM) approach is demonstrated to be a very useful tool for exploring structures and functions of a photoactive protein with the chromophore and its protein surroundings treated by the QM and MM methods, respectively. In this review, we summarize the basic formulas of the QM/MM approach and emphasize its applications to excited states and photoreactions of chromophores in rhodopsin protein, photoactive yellow protein, and green fluorescent protein.
Collapse
Affiliation(s)
- Lihong Liu
- Key Laboratory of Theoretical and Computational Photochemistry, Ministry of Education, College of Chemistry, Beijing Normal University, Beijing, China.
| | - Ganglong Cui
- Key Laboratory of Theoretical and Computational Photochemistry, Ministry of Education, College of Chemistry, Beijing Normal University, Beijing, China
| | - Wei-Hai Fang
- Key Laboratory of Theoretical and Computational Photochemistry, Ministry of Education, College of Chemistry, Beijing Normal University, Beijing, China
| |
Collapse
|
44
|
Nikiforov A, Gamez JA, Thiel W, Huix-Rotllant M, Filatov M. Assessment of approximate computational methods for conical intersections and branching plane vectors in organic molecules. J Chem Phys 2015; 141:124122. [PMID: 25273427 DOI: 10.1063/1.4896372] [Citation(s) in RCA: 59] [Impact Index Per Article: 6.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/03/2023] Open
Abstract
Quantum-chemical computational methods are benchmarked for their ability to describe conical intersections in a series of organic molecules and models of biological chromophores. Reference results for the geometries, relative energies, and branching planes of conical intersections are obtained using ab initio multireference configuration interaction with single and double excitations (MRCISD). They are compared with the results from more approximate methods, namely, the state-interaction state-averaged restricted ensemble-referenced Kohn-Sham method, spin-flip time-dependent density functional theory, and a semiempirical MRCISD approach using an orthogonalization-corrected model. It is demonstrated that these approximate methods reproduce the ab initio reference data very well, with root-mean-square deviations in the optimized geometries of the order of 0.1 Å or less and with reasonable agreement in the computed relative energies. A detailed analysis of the branching plane vectors shows that all currently applied methods yield similar nuclear displacements for escaping the strong non-adiabatic coupling region near the conical intersections. Our comparisons support the use of the tested quantum-chemical methods for modeling the photochemistry of large organic and biological systems.
Collapse
Affiliation(s)
- Alexander Nikiforov
- Max-Planck-Institut für Kohlenforschung, Kaiser-Wilhelm-Platz 1, D-45470 Mülheim an der Ruhr, Germany
| | - Jose A Gamez
- Max-Planck-Institut für Kohlenforschung, Kaiser-Wilhelm-Platz 1, D-45470 Mülheim an der Ruhr, Germany
| | - Walter Thiel
- Max-Planck-Institut für Kohlenforschung, Kaiser-Wilhelm-Platz 1, D-45470 Mülheim an der Ruhr, Germany
| | - Miquel Huix-Rotllant
- Institute of Physical and Theoretical Chemistry, Goethe University Frankfurt, Max-von-Laue-Str. 7, D-60438 Frankfurt am Main, Germany
| | - Michael Filatov
- Mulliken Center for Theoretical Chemistry, Institut für Physikalische und Theoretische Chemie, Universität Bonn, Beringstr. 4, D-53115 Bonn, Germany
| |
Collapse
|
45
|
Bishop MM, Roscioli JD, Ghosh S, Mueller JJ, Shepherd NC, Beck WF. Vibrationally Coherent Preparation of the Transition State for Photoisomerization of the Cyanine Dye Cy5 in Water. J Phys Chem B 2015; 119:6905-15. [DOI: 10.1021/acs.jpcb.5b02391] [Citation(s) in RCA: 19] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Michael M. Bishop
- Department of Chemistry, Michigan State University, East Lansing, Michigan 48824-1322, United States
| | - Jerome D. Roscioli
- Department of Chemistry, Michigan State University, East Lansing, Michigan 48824-1322, United States
| | - Soumen Ghosh
- Department of Chemistry, Michigan State University, East Lansing, Michigan 48824-1322, United States
| | - Jenny Jo Mueller
- Department of Chemistry, Michigan State University, East Lansing, Michigan 48824-1322, United States
| | - Nolan C. Shepherd
- Department of Chemistry, Michigan State University, East Lansing, Michigan 48824-1322, United States
| | - Warren F. Beck
- Department of Chemistry, Michigan State University, East Lansing, Michigan 48824-1322, United States
| |
Collapse
|
46
|
Halász GJ, Vibók Á, Cederbaum LS. Direct Signature of Light-Induced Conical Intersections in Diatomics. J Phys Chem Lett 2015; 6:348-54. [PMID: 26261946 DOI: 10.1021/jz502468d] [Citation(s) in RCA: 57] [Impact Index Per Article: 6.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/20/2023]
Abstract
Nonadiabatic effects are ubiquitous in physics, chemistry, and biology. They are strongly amplified by conical intersections (CIs), which are degeneracies between electronic states of triatomic or larger molecules. A few years ago it was revealed that CIs in molecular systems can be formed by laser light, even in diatomics. Because of the prevailing strong nonadiabatic couplings, the existence of such laser-induced conical intersections (LICIs) may considerably change the dynamical behavior of molecular systems. By analyzing the photodissociation process of the D2+ molecule carefully, we found a robust effect in the angular distribution of the photofragments that serves as a direct signature of the LICI, providing undoubted evidence of its existence.
Collapse
Affiliation(s)
- G J Halász
- †Department of Information Technology, University of Debrecen, P.O. Box 12, H-4010 Debrecen, Hungary
| | - Á Vibók
- ‡Department of Theoretical Physics, University of Debrecen, P.O. Box 5, H-4010 Debrecen, Hungary
| | - L S Cederbaum
- §Theoretische Chemie, Physikalish-Chemisches Institut, Universität Heidelberg, Im Neuenheimer Feld 229, D-69120 Heidelberg, Germany
| |
Collapse
|
47
|
Schapiro I, Roca-Sanjuán D, Lindh R, Olivucci M. A surface hopping algorithm for nonadiabatic minimum energy path calculations. J Comput Chem 2015; 36:312-20. [DOI: 10.1002/jcc.23805] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/10/2014] [Revised: 10/21/2014] [Accepted: 11/16/2014] [Indexed: 11/08/2022]
Affiliation(s)
- Igor Schapiro
- Department of Chemistry; Bowling Green State University; Bowling Green Ohio 43403
| | - Daniel Roca-Sanjuán
- Instituto de Ciencia Molecular; Universitat de València; P. O. Box 22085 València 46071 Spain
| | - Roland Lindh
- Department of Chemistry-Ångström; Theoretical Chemistry Programme, Uppsala University; P. O. Box 518 Uppsala 75120 Sweden
- Uppsala Center for Computational Chemistry - UC 3; Uppsala University; P. O. Box 518 Uppsala 75120 Sweden
| | - Massimo Olivucci
- Department of Chemistry; Bowling Green State University; Bowling Green Ohio 43403
- Dipartimento di Biotechnologie, Chimica e Farmacia; Università di Siena; Siena 53100 Italy
| |
Collapse
|
48
|
Huix-Rotllant M, Nikiforov A, Thiel W, Filatov M. Description of Conical Intersections with Density Functional Methods. Top Curr Chem (Cham) 2015; 368:445-76. [PMID: 25896441 DOI: 10.1007/128_2015_631] [Citation(s) in RCA: 30] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/23/2022]
Abstract
Conical intersections are perhaps the most significant mechanistic features of chemical reactions occurring through excited states. By providing funnels for efficient non-adiabatic population transfer, conical intersections govern the branching ratio of products of such reactions, similar to what the transition states do for ground-state reactivity. In this regard, intersections between the ground and the lowest excited states play a special role, and the correct description of the potential energy surfaces in their vicinity is crucial for understanding the mechanism and dynamics of excited-state reactions. The methods of density functional theory, such as time-dependent density functional theory, are widely used to describe the excited states of large molecules. However, are these methods suitable for describing the conical intersections or do they lead to artifacts and, consequently, to erroneous description of reaction dynamics? Here we address the first part of this question and analyze the ability of several density functional approaches, including the linear-response time-dependent approach as well as the spin-flip and ensemble formalisms, to provide the correct description of conical intersections and the potential energy surfaces in their vicinity. It is demonstrated that the commonly used linear-response time-dependent theory does not yield a proper description of these features and that one should instead use alternative computational approaches.
Collapse
Affiliation(s)
- Miquel Huix-Rotllant
- Institute of Physical and Theoretical Chemistry, Goethe University Frankfurt, Max-von-Laue-Str. 7, 60438, Frankfurt am Main, Germany
| | | | | | | |
Collapse
|
49
|
Quantum Chemistry of Excited States in Polyhedral Boranes. CHALLENGES AND ADVANCES IN COMPUTATIONAL CHEMISTRY AND PHYSICS 2015. [DOI: 10.1007/978-3-319-22282-0_4] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 02/11/2023]
|
50
|
Kubli-Garfias C, Vázquez-Ramírez R, Cabrera-Vivas BM, Gómez-Reyes B, Ramírez JC. Atomic and molecular analysis highlights the biophysics of unprotonated and protonated retinal in UV and scotopic vision. Photochem Photobiol Sci 2015; 14:1660-72. [DOI: 10.1039/c5pp00091b] [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
Unprotonated (UR) and protonated (PR) retinal have marked atomic and molecular differences in cis and trans configurations. In conclusion, UR and PR uphold UV and light vision through their different biophysical properties.
Collapse
Affiliation(s)
- Carlos Kubli-Garfias
- Instituto de Investigaciones Biomédicas
- Universidad Nacional Autónoma de México
- México
- Mexico
| | | | | | | | - Juan Carlos Ramírez
- Facultad de Ciencias Químicas
- Benemérita Universidad Autónoma de Puebla
- Puebla 72530
- Mexico
| |
Collapse
|