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

Systematic Excited State Studies of Reversibly Switchable Fluorescent Proteins

The reversibly switchable fluorescent proteins Dronpa, rsFastLime, rsKame, Padron, and bsDronpa feature the same chromophore but display a 40 nm variation in absorption maxima and an only 18 nm variation in emission maxima. In the present contribution, we employ QM/MM models to investigate the mechanism of such remarkably different spectral variations, which are caused by just a few amino acid replacements. We show that the models, which are based on CASPT2//CASSCF level of QM theory, reproduce the observed trends in absorption maxima, with only a 3.5 kcal/mol blue-shift, and in emission maxima, with an even smaller 1.5 kcal/mol blue-shift with respect to the observed quantities. In order to explain the variations across the series, we look at the chromophore's electronic structure change during absorption and emission. Such analysis indicates that a change in charge-transfer character, which is more pronounced during absorption, triggers a cascade of hydrogen-bond-network rearrangements, suggesting preparation to an isomerization event. We also show how the contribution of Arg 89 and Arg 64 residues to the chromophore conformational changes correlate with the spectral variations in absorption and emission. Furthermore, we describe how the conical intersection stability is related to the protein's photophysical properties. While for the Dronpa, rsFastLime, and rsKame triad, the stability correlates with the photoswitching speed, this does not happen for bsDronpa and Padron, suggesting a less obvious photoisomerization mechanism.

Theoretical Modeling of Singlet Fission

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.

Impacts of hydroxylation on the photophysics of chalcones: insights into the relation between the chemical composition and the electronic structure

A combined theoretical/experimental study of the photoreactivity of two flavylium-derived chalcones, 2,4,4'-trihydroxychalcone and 2,4'-dihydroxychalcone, at the multiconfigurational wavefunction level of theory (CASSCF//CASPT2) in vacuo and in an implicit solvent (water, treated as a polarisable continuum) and by means of linear absorption spectroscopy is presented. The photosensitivity of flavium salts is expressed in the ability of their chalcone form to undergo a cis-trans isomerisation which has found application in logical networks. Despite a considerable amount of experimental data documenting the dependence of the isomerisation on solvent, pH and temperature, the knowledge of how chalcones process energy under various conditions at the molecular level is still scarce. On the example of 2,4,4'-trihydroxychalcone we unravel the complex excited state deactivation mechanism in vacuo involving ultrafast decay through conical intersections, formation of twisted intramolecular charge transfer species, intramolecular proton transfer and inter system crossings. Furthermore, we rationalise the observed discrepancies in the linear absorption spectra of 2,4,4'-trihydroxychalcone and 2,4'-dihydroxychalcone, thereby establishing a link between the functionalisation pattern and the observed spectral properties.

Accurate Excited-State Geometries: A CASPT2 and Coupled-Cluster Reference Database for Small Molecules

We present an investigation of the excited-state structural parameters determined for a large set of small compounds with the dual goals of defining reference values for further works and assessing the quality of the geometries obtained with relatively cheap computational approaches. In the first stage, we compare the excited-state geometries obtained with ADC(2), CC2, CCSD, CCSDR(3), CC3, and CASPT2 and large atomic basis sets. It is found that CASPT2 and CC3 results are generally in very good agreement with one another (typical differences of ca. 3 × 10^–3 Å) when all electrons are correlated and when the aug-cc-pVTZ atomic basis set is employed with both methods. In a second stage, a statistical analysis reveals that, on the one hand, the excited-state (ES) bond lengths are much more sensitive to the selected level of theory than their ground-state (GS) counterparts and, on the other hand, that CCSDR(3) is probably the most cost-effective method delivering accurate structures. Indeed, CCSD tends to provide too compact multiple bond lengths on an almost systematic basis, whereas both CC2 and ADC(2) tend to exaggerate these bond distances, with more erratic error patterns, especially for the latter method. The deviations are particularly marked for the polarized CO and CN bonds, as well as for the puckering angle in formaldehyde homologues. In the last part of this contribution, we provide a series of CCSDR(3) GS and ES geometries of medium-sized molecules to be used as references in further investigations.

Excited-state minima and emission energies of retinal chromophore analogues: Performance of CASSCF and CC2 methods as compared with CASPT2

This study provides gas-phase S₁ 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 S₁ 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 S₁ minima and their spectroscopic and electronic properties. However, this picture deteriorates for twisted minima. As expected, the CC2 description of the S₂ 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.

An experimental comparison of human and bovine rhodopsin provides insight into the molecular basis of retinal disease

Rhodopsin is the visual pigment that mediates dim-light vision in vertebrates and is a model system for the study of retinal disease. The majority of rhodopsin experiments are performed using bovine rhodopsin; however, recent evidence suggests that significant functional differences exist among mammalian rhodopsins. In this study, we identify differences in both thermal decay and light-activated retinal release rates between bovine and human rhodopsin and perform mutagenesis studies to highlight two clusters of substitutions that contribute to these differences. We also demonstrate that the retinitis pigmentosa-associated mutation G51A behaves differently in human rhodopsin compared to bovine rhodopsin and determine that the thermal decay rate of an ancestrally reconstructed mammalian rhodopsin displays an intermediate phenotype compared to the two extant pigments.

pH-dependent absorption spectra of rhodopsin mutant E113Q: On the role of counterions and protein

The absorption spectra of bovine rhodopsin mutant E113Q in solutions were investigated at the molecular level by using a hybrid quantum mechanics/molecular mechanics (QM/MM) method. The calculations suggest the mechanism of the absorption variations of E113Q at different pH values. The results indicate that the polarizations of the counterions in the vicinity of Schiff base under protonation and unprotonation states of the mutant E113Q would be a crucial factor to change the energy gap of the retinal to tune the absorption spectra. Glu-181 residue, which is close to the chromophore, cannot serve as the counterion of the protonated Schiff base of E113Q in dark state. Moreover, the results of the absorption maximum in mutant E113Q with the various anions (Cl^-, Br^-, I^- and NO₃^-) manifested that the mutant E113Q could have the potential for use as a template of anion biosensors at visible wavelength.

Theoretical description of protein field effects on electronic excitations of biological chromophores

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.

Assessment of the Potential Energy Hypersurfaces in Thymine within Multiconfigurational Theory: CASSCF vs. CASPT2

The present study provides new insights into the topography of the potential energy hypersurfaces (PEHs) of the thymine nucleobase in order to rationalize its main ultrafast photochemical decay paths by employing two methodologies based on the complete active space self-consistent field (CASSCF) and the complete active space second-order perturbation theory (CASPT2) methods: (i) CASSCF optimized structures and energies corrected with the CASPT2 method at the CASSCF geometries and (ii) CASPT2 optimized geometries and energies. A direct comparison between these strategies is drawn, yielding qualitatively similar results within a static framework. A number of analyses are performed to assess the accuracy of these different computational strategies under study based on a variety of numerical thresholds and optimization methods. Several basis sets and active spaces have also been calibrated to understand to what extent they can influence the resulting geometries and subsequent interpretation of the photochemical decay channels. The study shows small discrepancies between CASSCF and CASPT2 PEHs, displaying a shallow planar or twisted ¹(ππ*) minimum, respectively, and thus featuring a qualitatively similar scenario for supporting the ultrafast bi-exponential deactivation registered in thymine upon UV-light exposure. A deeper knowledge of the PEHs at different levels of theory provides useful insight into its correct characterization and subsequent interpretation of the experimental observations. The discrepancies displayed by the different methods studied here are then discussed and framed within their potential consequences in on-the-fly non-adiabatic molecular dynamics simulations, where qualitatively diverse outcomes are expected.

Role of the Molecular Environment in Flavoprotein Color and Redox Tuning: QM Cluster versus QM/MM Modeling

We investigate the origin of the excitation energy shifts induced by the apoprotein in the active site of the bacterial photoreceptor BLUF (Blue Light sensor Using Flavin adenine dinucleotide). In order to compute the vertical excitation energies of three low-lying electronic states, including two π-π* states of flavin (S1 and S2) and a π-π* tyrosine-flavin electron-transfer state (ET), with respect to the energy of the closed-shell ground state (S0), we prepared alternative quantum mechanical (QM) cluster and quantum mechanics/molecular mechanics (QM/MM) models. We found that the excitation energies computed with both types of models correlate with the magnitude of the charge transfer character of the excitation. Accordingly, we conclude that the small charge transfer character of the light absorbing S0-S1 transition and the substantial charge transfer character of the nonabsorbing but redox active S0-ET transition explain the small color changes but substantial redox tuning in BLUF and also in other flavoproteins. Further analysis showed that redox tuning is governed by the electrostatic interaction in the QM/MM model and transfer of charge between the active site and its environment in the QM cluster. Moreover, the wave function polarization of the QM subsystem by the MM subsystem influences the magnitude of the charge transfer, resulting in the QM/MM and QM excitation energies that are not entirely consistent.

Quantum Monte Carlo Treatment of the Charge Transfer and Diradical Electronic Character in a Retinal Chromophore Minimal Model

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.

Intersystem Crossing Pathways in the Noncanonical Nucleobase 2-Thiouracil: A Time-Dependent Picture

The deactivation mechanism after ultraviolet irradiation of 2-thiouracil has been investigated using nonadiabatic dynamics simulations at the MS-CASPT2 level of theory. It is found that after excitation the S2 quickly relaxes to S1, and from there intersystem crossing takes place to both T2 and T1 with a time constant of 400 fs and a triplet yield above 80%, in very good agreement with recent femtosecond experiments in solution. Both indirect S1 → T2 → T1 and direct S1 → T1 pathways contribute to intersystem crossing, with the former being predominant. The results contribute to the understanding of how some noncanonical nucleobases respond to harmful ultraviolet light, which could be relevant for prospective photochemotherapeutic applications.

Beyond Time-Dependent Density Functional Theory Using Only Single Excitations: Methods for Computational Studies of Excited States in Complex Systems

Single-excitation methods, namely, configuration interaction singles and time-dependent density functional theory (TDDFT), along with semiempirical versions thereof, represent the most computationally affordable electronic structure methods for describing electronically excited states, scaling as [Formula: see text] absent further approximations. This relatively low cost, combined with a treatment of electron correlation, has made TDDFT the most widely used excited-state quantum chemistry method over the past 20+ years. Nevertheless, certain inherent problems (beyond just the accuracy of this or that exchange-correlation functional) limit the utility of traditional TDDFT. For one, it affords potential energy surfaces whose topology is incorrect in the vicinity of any conical intersection (CI) that involves the ground state. Since CIs are the conduits for transitions between electronic states, the TDDFT description of photochemistry (internal conversion and intersystem crossing) is therefore suspect. Second, the [Formula: see text] cost can become prohibitive in large systems, especially those that involve multiple electronically coupled chromophores, for example, the antennae structures of light-harvesting complexes or the conjugated polymers used in organic photovoltaics. In such cases, the smallest realistic mimics might already be quite large from the standpoint of ab initio quantum chemistry. This Account describes several new computational methods that address these problems. Topology around a CI can be rigorously corrected using a "spin-flip" version of TDDFT, which involves an α → β spin-flipping transition in addition to occupied → virtual excitation of one electron. Within this formalism, singlet states are generated via excitation from a high-spin triplet reference state, doublets from a quartet, etc. This provides a more balanced treatment of electron correlation between ground and excited states. Spin contamination is problematic away from the Franck-Condon region, but we describe a "spin-complete" version of the theory in which proper spin eigenstates are obtained by construction. For systems of coupled chromophores, we have developed an ab initio version of the Frenkel-Davydov exciton model in which collective excitations of the system are expanded in a basis of excited states computed for individual chromophores. The monomer calculations are trivially parallelizable, as is computation of the coupling matrix elements needed to construct the exciton Hamiltonian, and systems containing hundreds of chromophores can be tackled on commodity hardware. This enables calculations on organic semiconductors, where even small model systems exhibit a semicontinuum of excited states that renders traditional TDDFT computationally challenging. Despite including only single excitations on each monomer, the exciton model can describe entangled spins on two or more monomers, an effect that is responsible for excitation energy transfer between chromophores, for example, in singlet fission. Excitonic approximations can also be applied to the TDDFT equations themselves, and a particularly promising application is to describe the effects of environment on an excitation that is localized on a single chromophore. This "local excitation approximation" to TDDFT allows an essentially arbitrary number of solvent molecules to be included in the calculation in a highly parallelizable way such that the time-to-solution increases only very slowly as additional solvent molecules are added. It is therefore possible to converge the calculation with respect to describing an ever-larger portion of the environment at a quantum-mechanical level.

Multiconfigurational Second-Order Perturbation Theory with Frozen Natural Orbitals Extended to the Treatment of Photochemical Problems

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.

Efficient and Flexible Computation of Many-Electron Wave Function Overlaps

A new algorithm for the computation of the overlap between many-electron wave functions is described. This algorithm allows for the extensive use of recurring intermediates and thus provides high computational efficiency. Because of the general formalism employed, overlaps can be computed for varying wave function types, molecular orbitals, basis sets, and molecular geometries. This paves the way for efficiently computing nonadiabatic interaction terms for dynamics simulations. In addition, other application areas can be envisaged, such as the comparison of wave functions constructed at different levels of theory. Aside from explaining the algorithm and evaluating the performance, a detailed analysis of the numerical stability of wave function overlaps is carried out, and strategies for overcoming potential severe pitfalls due to displaced atoms and truncated wave functions are presented.

Rhodopsins carrying modified chromophores--the 'making of', structural modelling and their light-induced reactivity

A series of vitamin-A aldehydes (retinals) with modified alkyl group substituents (9-demethyl-, 9-ethyl-, 9-isopropyl-, 10-methyl, 10-methyl-13-demethyl-, and 13-demethyl retinal) was synthesized and their 11-cis isomers were used as chromophores to reconstitute the visual pigment rhodopsin. Structural changes were selectively introduced around the photoisomerizing C11=C12 bond. The effect of these structural changes on rhodopsin formation and bleaching was determined. Global fit of assembly kinetics yielded lifetimes and spectral features of the assembly intermediates. Rhodopsin formation proceeds stepwise with prolonged lifetimes especially for 9-demethyl retinal (longest lifetime τ3 = 7500 s, cf., 3500 s for retinal), and for 10-methyl retinal (τ3 = 7850 s). These slowed-down processes are interpreted as either a loss of fixation (9dm) or an increased steric hindrance (10me) during the conformational adjustment within the protein. Combined quantum mechanics and molecular mechanics (QM/MM) simulations provided structural insight into the retinal analogues-assembled, full-length rhodopsins. Extinction coefficients, quantum yields and kinetics of the bleaching process (μs-to-ms time range) were determined. Global fit analysis yielded lifetimes and spectral features of bleaching intermediates, revealing remarkably altered kinetics: whereas the slowest process of wild-type rhodopsin and of bleached and 11-cis retinal assembled rhodopsin takes place with lifetimes of 7 and 3.8 s, respectively, this process for 10-methyl-13-demethyl retinal was nearly 10 h (34670 s), coming to completion only after ca. 50 h. The structural changes in retinal derivatives clearly identify the precise interactions between chromophore and protein during the light-induced changes that yield the outstanding efficiency of rhodopsin.

Probing the Photodynamics of Rhodopsins with Reduced Retinal Chromophores

While the light-induced population dynamics of different photoresponsive proteins has been investigated spectroscopically, systematic computational studies have not yet been possible due to the phenomenally high cost of suitable high level quantum chemical methods and the need of propagating hundreds, if not thousands, of nonadiabatic trajectories. Here we explore the possibility of studying the photodynamics of rhodopsins by constructing and investigating quantum mechanics/molecular mechanics (QM/MM) models featuring reduced retinal chromophores. In order to do so we use the sensory rhodopsin found in the cyanobacterium Anabaena PCC7120 (ASR) as a benchmark system. We find that the basic mechanistic features associated with the excited state dynamics of ASR QM/MM models are reproduced using models incorporating a minimal (i.e., three double-bond) chromophore. Furthermore, we show that ensembles of nonadiabatic ASR trajectories computed using the same abridged models replicate, at both the CASPT2 and CASSCF levels of theory, the trends in spectroscopy and lifetimes estimated using unabridged models and observed experimentally at room temperature. We conclude that a further expansion of these studies may lead to low-cost QM/MM rhodopsin models that may be used as effective tools in high-throughput in silico mutant screening.

Molecular bases for the selection of the chromophore of animal rhodopsins

The functions of microbial and animal rhodopsins are triggered by the isomerization of their all-trans and 11-cis retinal chromophores, respectively. To lay the molecular basis driving the evolutionary transition from the all-trans to the 11-cis chromophore, multiconfigurational quantum chemistry is used to compare the isomerization mechanisms of the sensory rhodopsin from the cyanobacterium Anabaena PCC 7120 (ASR) and of the bovine rhodopsin (Rh). It is found that, despite their evolutionary distance, these eubacterial and vertebrate rhodopsins start to isomerize via distinct implementations of the same bicycle-pedal mechanism originally proposed by Warshel [Warshel A (1976) Nature 260:678-683]. However, by following the electronic structure changes of ASR (featuring the all-trans chromophore) during the isomerization, we find that ASR enters a region of degeneracy between the first and second excited states not found in Rh (featuring the 11-cis chromophore). We show that such degeneracy is modulated by the preorganized structure of the chromophore and by the position of the reactive double bond. It is argued that the optimization of the electronic properties of the chromophore, which affects the photoisomerization efficiency and the thermal isomerization barrier, provided a key factor for the emergence of the striking amino acid sequence divergence observed between the microbial and animal rhodopsins.

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

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.

A Static Picture of the Relaxation and Intersystem Crossing Mechanisms of Photoexcited 2-Thiouracil

Accurate excited-state quantum chemical calculations on 2-thiouracil, employing large active spaces and up to quadruple-ζ quality basis sets in multistate complete active space perturbation theory calculations, are reported. The results suggest that the main relaxation path for 2-thiouracil after photoexcitation should be S2 → S1 → T2 → T1, and that this relaxation occurs on a subpicosecond time scale. There are two deactivation pathways from the initially excited bright S2 state to S1, one of which is nearly barrierless and should promote ultrafast internal conversion. After relaxation to the S1 minimum, small singlet-triplet energy gaps and spin-orbit couplings of about 130 cm(-1) are expected to facilitate intersystem crossing to T2, from where very fast internal conversion to T1 occurs. An important finding is that 2-thiouracil shows strong pyramidalization at the carbon atom of the thiocarbonyl group in several excited states.

Excited States and Photochemistry of Chromophores in the Photoactive Proteins Explored by the Combined Quantum Mechanical and Molecular Mechanical Calculations

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.

A Light-Induced Reaction with Oxygen Leads to Chromophore Decomposition and Irreversible Photobleaching in GFP-Type Proteins

Photobleaching and photostability of proteins of the green fluorescent protein (GFP) family are crucially important for practical applications of these widely used biomarkers. On the basis of simulations, we propose a mechanism for irreversible bleaching in GFP-type proteins under intense light illumination. The key feature of the mechanism is a photoinduced reaction of the chromophore with molecular oxygen (O2) inside the protein barrel leading to the chromophore's decomposition. Using quantum mechanics/molecular mechanics (QM/MM) modeling we show that a model system comprising the protein-bound Chro(-) and O2 can be excited to an electronic state of the intermolecular charge-transfer (CT) character (Chro(•)···O2(-•)). Once in the CT state, the system undergoes a series of chemical reactions with low activation barriers resulting in the cleavage of the bridging bond between the phenolic and imidazolinone rings and disintegration of the chromophore.

100 fs photo-isomerization with vibrational coherences but low quantum yield in Anabaena Sensory Rhodopsin

Counter-intuitive photochemistry: in Anabaena Sensory Rhodopsin, the retinal 13-cis isomer isomerizes much faster than all-trans ASR, but with a 3-times lower quantum yield.

Description of Conical Intersections with Density Functional Methods

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.

Spectral properties and isomerisation path of retinal in C1C2 channelrhodopsin

Computed torsion profiles along the reactive coordinate in S₁reveal a two-path deactivation mechanism for retinal in C1C2 channelrhodopsin.

Quantum dynamics of the photostability of pyrazine

We propose a new mechanism for the radiationless decay of photoexcited pyrazine to its ground electronic state involving a conical intersection between the dark A_u(nπ) state and the ground state.

Impacts of retinal polyene (de)methylation on the photoisomerization mechanism and photon energy storage of rhodopsin

Similar to native rhodopsin, a two-mode space-saving isomerization mechanism drives the photoreaction in (de)methylated rhodopsin analogues.