Dynamic Electron Correlation Effects on the Ground State Potential Energy Surface of a Retinal Chromophore Model.

The ground state potential energy surface of the retinal chromophore of visual pigments (e.g., bovine rhodopsin) features a low-lying conical intersection surrounded by regions with variable charge-transfer and diradical electronic structures. This implies that dynamic electron correlation may have a large effect on the shape of the force fields driving its reactivity. To investigate this effect, we focus on mapping the potential energy for three paths located along the ground state CASSCF potential energy surface of the penta-2,4-dieniminium cation taken as a minimal model of the retinal chromophore. The first path spans the bond length alternation coordinate and intercepts a conical intersection point. The other two are minimum energy paths along two distinct but kinetically competitive thermal isomerization coordinates. We show that the effect of introducing the missing dynamic electron correlation variationally (with MRCISD) and perturbatively (with the CASPT2, NEVPT2, and XMCQDPT2 methods) leads, invariably, to a stabilization of the regions with charge transfer character and to a significant reshaping of the reference CASSCF potential energy surface and suggesting a change in the dominating isomerization mechanism. The possible impact of such a correction on the photoisomerization of the retinal chromophore is discussed.

[1]  Riccardo Basosi,et al.  Properties of the emitting state of the green fluorescent protein resolved at the CASPT2//CASSCF/CHARMM level. , 2005, Journal of the American Chemical Society.

[2]  Celestino Angeli,et al.  A quasidegenerate formulation of the second order n-electron valence state perturbation theory approach. , 2004, The Journal of chemical physics.

[3]  A. Bottoni,et al.  Product formation in rhodopsin by fast hydrogen motions. , 2011, Physical chemistry chemical physics : PCCP.

[4]  N. Ferré,et al.  Characterization of the conical intersection of the visual pigment rhodopsin at the CASPT2//CASSCF/AMBER level of theory , 2006 .

[5]  Luis Serrano-Andrés,et al.  The multi-state CASPT2 method , 1998 .

[6]  M. Barbatti,et al.  Ionization potentials of adenine along the internal conversion pathways. , 2011, Physical chemistry chemical physics : PCCP.

[7]  FRANCESCO AQUILANTE,et al.  MOLCAS 7: The Next Generation , 2010, J. Comput. Chem..

[8]  Isaiah Shavitt,et al.  New implementation of the graphical unitary group approach for multireference direct configuration interaction calculations , 2009 .

[9]  B. Roos,et al.  A modified definition of the zeroth-order Hamiltonian in multiconfigurational perturbation theory (CASPT2) , 2004 .

[10]  Toru Shiozaki,et al.  Communication: extended multi-state complete active space second-order perturbation theory: energy and nuclear gradients. , 2011, The Journal of chemical physics.

[11]  W. Jiang,et al.  Multireference composite approaches for the accurate study of ground and excited electronic states: C2, N2, and O2. , 2011, The Journal of chemical physics.

[12]  H. L. Thanh,et al.  Effects of solvents on a model of rhodopsin consisting of a conjugated imine and substituted acetic acids , 1984 .

[13]  Celestino Angeli,et al.  Introduction of n-electron valence states for multireference perturbation theory , 2001 .

[14]  Klaus Ruedenberg,et al.  Potential energy surfaces near intersections , 1991 .

[15]  Roland Lindh,et al.  The ultrafast photoisomerizations of rhodopsin and bathorhodopsin are modulated by bond length alternation and HOOP driven electronic effects. , 2011, Journal of the American Chemical Society.

[16]  F. Bernardi,et al.  Reaction path analysis of the "tunable" photoisomerization selectivity of free and locked retinal chromophores. , 2002, Journal of the American Chemical Society.

[17]  J. Hynes,et al.  Conical intersections in solution: non-equilibrium versus equilibrium solvation , 2006 .

[18]  Ranbir Singh,et al.  J. Mol. Struct. (Theochem) , 1996 .

[19]  Mikhail N. Ryazantsev,et al.  Computational Photobiology and Beyond , 2010 .

[20]  Thomas Müller,et al.  High-level multireference methods in the quantum-chemistry program system COLUMBUS: Analytic MR-CISD and MR-AQCC gradients and MR-AQCC-LRT for excited states, GUGA spin–orbit CI and parallel CI density , 2001 .

[21]  Toshifumi Mori,et al.  Conical intersections of free energy surfaces in solution: effect of electron correlation on a protonated Schiff base in methanol solution. , 2010, The Journal of chemical physics.

[22]  Russell M. Pitzer,et al.  A progress report on the status of the COLUMBUS MRCI program system , 1988 .

[23]  Massimo Olivucci,et al.  CASPT2//CASSCF and TDDFT//CASSCF mapping of the excited state isomerization path of a minimal model of the retinal chromophore , 2004 .

[24]  A. Granovsky,et al.  Extended multi-configuration quasi-degenerate perturbation theory: the new approach to multi-state multi-reference perturbation theory. , 2011, The Journal of chemical physics.

[25]  Omar Valsson,et al.  Photoisomerization of Model Retinal Chromophores: Insight from Quantum Monte Carlo and Multiconfigurational Perturbation Theory , 2010 .

[26]  Daniel Escudero,et al.  Progress and challenges in the calculation of electronic excited states. , 2012, Chemphyschem : a European journal of chemical physics and physical chemistry.

[27]  Massimo Olivucci,et al.  Ground and excited state CASPT2 geometry optimizations of small organic molecules , 2003, J. Comput. Chem..

[28]  Massimo Olivucci,et al.  Relationship between photoisomerization path and intersection space in a retinal chromophore model. , 2003, Journal of the American Chemical Society.

[29]  M. Olivucci,et al.  Origin of the absorption maxima of the photoactive yellow protein resolved via ab initio multiconfigurational methods. , 2008, The journal of physical chemistry. B.

[30]  Marco Garavelli,et al.  The C 5 H 6 NH 2 + Protonated Shiff Base: An ab Initio Minimal Model for Retinal Photoisomerization , 1997 .

[31]  Massimo Olivucci,et al.  The Molecular Mechanism of Thermal Noise in Rod Photoreceptors , 2012, Science.

[32]  C. Angeli,et al.  The lowest singlet states of octatetraene revisited. , 2011, The Journal of chemical physics.

[33]  Ernest R. Davidson,et al.  The Importance of Including Dynamic Electron Correlation in ab initio Calculations , 1996 .

[34]  M. Olivucci,et al.  Photoisomerization acceleration in retinal protonated Schiff-base models. , 2003, Photochemical & photobiological sciences : Official journal of the European Photochemistry Association and the European Society for Photobiology.

[35]  Manuel Calderón Sánchez,et al.  Solvent Effects on the Radiative and Nonradiative Decay of a Model of the Rhodopsin Chromophore. , 2011, Journal of chemical theory and computation.

[36]  Björn O. Roos,et al.  Second-order perturbation theory with a complete active space self-consistent field reference function , 1992 .

[37]  L. Cederbaum,et al.  Environmental effects on a conical intersection: a model study. , 2004, Faraday discussions.

[38]  R. Cimiraglia,et al.  New perspectives in multireference perturbation theory: the n-electron valence state approach , 2007 .

[39]  Massimo Olivucci,et al.  Probing the Photochemical Funnel of a Retinal Chromophore Model via Zero-Point Energy Sampling Semiclassical Dynamics. , 2004 .

[40]  C. Filippi,et al.  Excitation energies of retinal chromophores: critical role of the structural model. , 2012, Physical chemistry chemical physics : PCCP.

[41]  A. Zaitsevskii,et al.  Multiconfigurational second-order perturbative methods: Overview and comparison of basic properties , 1995 .

[42]  N. Ferré,et al.  The color of rhodopsins at the ab initio multiconfigurational perturbation theory resolution , 2006, Proceedings of the National Academy of Sciences.

[43]  M. Barbatti,et al.  Nonadiabatic excited-state dynamics of polar pi-systems and related model compounds of biological relevance. , 2008, Physical chemistry chemical physics : PCCP.

[44]  P. Knowles,et al.  An efficient internally contracted multiconfiguration–reference configuration interaction method , 1988 .

[45]  M. Garavelli,et al.  Modelling vibrational coherence in the primary rhodopsin photoproduct. , 2012, The Journal of chemical physics.

[46]  F. Bernardi,et al.  The short-chain acroleiniminium and pentadieniminium cations: towards a model for retinal photoisomerization. A CASSCF/PT2 study , 1999 .

[47]  Mihály Kállay,et al.  The barrier height of the F+H2 reaction revisited: coupled-cluster and multireference configuration-interaction benchmark calculations. , 2008, The Journal of chemical physics.

[48]  M. Képénekian,et al.  What zeroth-order Hamiltonian for CASPT2 adiabatic energetics of Fe(II)N(6) architectures? , 2009, The Journal of chemical physics.

[49]  N. Ferré,et al.  Modeling the fluorescence of protein-embedded tryptophans with ab initio multiconfigurational quantum chemistry: the limiting cases of parvalbumin and monellin. , 2009, The journal of physical chemistry. B.

[50]  H. Bernhard Schlegel,et al.  An improved algorithm for reaction path following , 1989 .

[51]  Massimo Olivucci,et al.  Structure, initial excited-state relaxation, and energy storage of rhodopsin resolved at the multiconfigurational perturbation theory level , 2004, Proceedings of the National Academy of Sciences of the United States of America.

[52]  N. Ferré,et al.  Tracking the excited-state time evolution of the visual pigment with multiconfigurational quantum chemistry , 2007, Proceedings of the National Academy of Sciences.

[53]  宁北芳,et al.  疟原虫var基因转换速率变化导致抗原变异[英]／Paul H, Robert P, Christodoulou Z, et al//Proc Natl Acad Sci U S A , 2005 .