1
|
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
|
2
|
Cuzzocrea A, Scemama A, Briels WJ, Moroni S, Filippi C. Variational Principles in Quantum Monte Carlo: The Troubled Story of Variance Minimization. J Chem Theory Comput 2020; 16:4203-4212. [PMID: 32419451 PMCID: PMC7365558 DOI: 10.1021/acs.jctc.0c00147] [Citation(s) in RCA: 14] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/26/2023]
Abstract
![]()
We
investigate the use of different variational principles in quantum
Monte Carlo, namely, energy and variance minimization, prompted by
the interest in the robust and accurate estimation of electronic excited
states. For two prototypical, challenging molecules, we readily reach
the accuracy of the best available reference excitation energies using
energy minimization in a state-specific or state-average fashion for
states of different or equal symmetry, respectively. On the other
hand, in variance minimization, where the use of suitable functionals
is expected to target specific states regardless of the symmetry,
we encounter severe problems for a variety of wave functions: as the
variance converges, the energy drifts away from that of the selected
state. This unexpected behavior is sometimes observed even when the
target is the ground state and generally prevents the robust estimation
of total and excitation energies. We analyze this problem using a
very simple wave function and infer that the optimization finds little
or no barrier to escape from a local minimum or local plateau, eventually
converging to a lower-variance state instead of the target state.
For the increasingly complex systems becoming in reach of quantum
Monte Carlo simulations, variance minimization with current functionals
appears to be an impractical route.
Collapse
Affiliation(s)
- Alice Cuzzocrea
- MESA+ Institute for Nanotechnology, University of Twente, P.O. Box 217, 7500 AE Enschede, The Netherlands
| | - Anthony Scemama
- Laboratoire de Chimie et Physique Quantiques, Université de Toulouse, CNRS, UPS, 118 route de Narbonne, 31062 Toulouse Cedex 09, France
| | - Wim J Briels
- MESA+ Institute for Nanotechnology, University of Twente, P.O. Box 217, 7500 AE Enschede, The Netherlands
| | - Saverio Moroni
- CNR-IOM DEMOCRITOS, Istituto Officina dei Materiali, Via Bonomea 265, I-34136 Trieste, Italy.,SISSA Scuola Internazionale Superiore di Studi Avanzati, Via Bonomea 265, I-34136 Trieste, Italy
| | - Claudia Filippi
- MESA+ Institute for Nanotechnology, University of Twente, P.O. Box 217, 7500 AE Enschede, The Netherlands
| |
Collapse
|
3
|
A DFT/TDDFT investigation on structure–photophysical properties relationship of phenothiazine derivatives with substitutions on C-3/N-10 sites. COMPUT THEOR CHEM 2019. [DOI: 10.1016/j.comptc.2019.112512] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/18/2022]
|
4
|
Dash M, Feldt J, Moroni S, Scemama A, Filippi C. Excited States with Selected Configuration Interaction-Quantum Monte Carlo: Chemically Accurate Excitation Energies and Geometries. J Chem Theory Comput 2019; 15:4896-4906. [PMID: 31348645 PMCID: PMC6740157 DOI: 10.1021/acs.jctc.9b00476] [Citation(s) in RCA: 24] [Impact Index Per Article: 4.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
We employ quantum Monte Carlo to obtain chemically accurate vertical and adiabatic excitation energies, and equilibrium excited-state structures for the small, yet challenging, formaldehyde and thioformaldehyde molecules. A key ingredient is a robust protocol to obtain balanced ground- and excited-state Jastrow-Slater wave functions at a given geometry, and to maintain such a balanced description as we relax the structure in the excited state. We use determinantal components generated via a selected configuration interaction scheme which targets the same second-order perturbation energy correction for all states of interest at different geometries, and fully optimize all variational parameters in the resultant Jastrow-Slater wave functions. Importantly, the excitation energies as well as the structural parameters in the ground and excited states are converged with very compact wave functions comprising few thousand determinants in a minimally augmented double-ζ basis set. These results are obtained already at the variational Monte Carlo level, the more accurate diffusion Monte Carlo method yielding only a small improvement in the adiabatic excitation energies. We find that matching Jastrow-Slater wave functions with similar variances can yield excitation energies compatible with our best estimates; however, the variance-matching procedure requires somewhat larger determinantal expansions to achieve the same accuracy, and it is less straightforward to adapt during structural optimization in the excited state.
Collapse
Affiliation(s)
- Monika Dash
- MESA+ Institute for Nanotechnology , University of Twente , P.O. Box 217, 7500 AE Enschede , The Netherlands
| | - Jonas Feldt
- MESA+ Institute for Nanotechnology , University of Twente , P.O. Box 217, 7500 AE Enschede , The Netherlands
| | - Saverio Moroni
- CNR-IOM DEMOCRITOS , Istituto Officina dei Materiali and SISSA Scuola Internazionale Superiore di Studi Avanzati , Via Bonomea 265 , I-34136 Trieste , Italy
| | - Anthony Scemama
- Laboratoire de Chimie et Physique Quantiques , Université de Toulouse, CNRS, UPS , 31062 Toulouse , France
| | - Claudia Filippi
- MESA+ Institute for Nanotechnology , University of Twente , P.O. Box 217, 7500 AE Enschede , The Netherlands
| |
Collapse
|
5
|
Lee S, Shostak S, Filatov M, Choi CH. Conical Intersections in Organic Molecules: Benchmarking Mixed-Reference Spin-Flip Time-Dependent DFT (MRSF-TD-DFT) vs Spin-Flip TD-DFT. J Phys Chem A 2019; 123:6455-6462. [PMID: 31283235 DOI: 10.1021/acs.jpca.9b06142] [Citation(s) in RCA: 22] [Impact Index Per Article: 4.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
The mixed-reference spin-flip time-dependent density functional theory (MRSF-TD-DFT) method eliminates the erroneous spin contamination of the SF-TD-DFT methodology, while retaining the conceptual and practical simplicity of the latter. The availability of the analytic gradient of the energy of the MRSF-TD-DFT response states enables automatic geometry optimization of the targeted states. Here, we apply the new method to optimize the geometry of several S1/S0 conical intersections occurring in typical organic molecules. We demonstrate that MRSF-TD-DFT is capable of producing the correct double-cone topology of the intersections and describing the geometry of the lowest-energy conical intersections and their relative energies with accuracy matching that of the best multireference wavefunction ab initio methods. In this regard, MRSF-TD-DFT differs from many popular single-reference methods, such as, e.g., the linear response TD-DFT method, which fail to produce the correct topology of the intersections. As the new methodology completely eliminates the ambiguity with the identification of the response states as proper singlets or triplets, which is plaguing the SF-TD-DFT calculations, it can be used for automatic geometry optimization and molecular dynamic simulations not requiring constant human intervention.
Collapse
Affiliation(s)
- Seunghoon Lee
- Department of Chemistry , Seoul National University , Seoul 151-747 , South Korea
| | - Svetlana Shostak
- Department of Chemistry , Kyungpook National University , Daegu 702-701 , South Korea
| | - Michael Filatov
- Department of Chemistry , Kyungpook National University , Daegu 702-701 , South Korea
| | - Cheol Ho Choi
- Department of Chemistry , Kyungpook National University , Daegu 702-701 , South Korea
| |
Collapse
|
6
|
Borin VA, Wiebeler C, Schapiro I. A QM/MM study of the initial excited state dynamics of green-absorbing proteorhodopsin. Faraday Discuss 2019; 207:137-152. [PMID: 29393940 DOI: 10.1039/c7fd00198c] [Citation(s) in RCA: 11] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/11/2023]
Abstract
The primary photochemical reaction of the green-absorbing proteorhodopsin is studied by means of a hybrid quantum mechanics/molecular mechanics (QM/MM) approach. The simulations are based on a homology model derived from the blue-absorbing proteorhodopsin crystal structure. The geometry of retinal and the surrounding sidechains in the protein binding pocket were optimized using the QM/MM method. Starting from this geometry the isomerization was studied with a relaxed scan along the C13[double bond, length as m-dash]C14 dihedral. It revealed an "aborted bicycle pedal" mechanism of isomerization that was originally proposed by Warshel for bovine rhodopsin and bacteriorhodopsin. However, the isomerization involved the concerted rotation about C13[double bond, length as m-dash]C14 and C15[double bond, length as m-dash]N, with the latter being highly twisted but not isomerized. Further, the simulation showed an increased steric interaction between the hydrogen at the C14 of the isomerizing bond and the hydroxyl group at the neighbouring tyrosine 200. In addition, we have simulated a nonadiabatic trajectory which showed the timing of the isomerization. In the first 20 fs upon excitation the order of the conjugated double and single bonds is inverted, consecutively the C13[double bond, length as m-dash]C14 rotation is activated for 200 fs until the S1-S0 transition is detected. However, the isomerization is reverted due to the specific interaction with the tyrosine as observed along the relaxed scan calculation. Our simulations indicate that the retinal - tyrosine 200 interaction plays an important role in the outcome of the photoisomerization.
Collapse
Affiliation(s)
- Veniamin A Borin
- Fritz Haber Center for Molecular Dynamics, Institute of Chemistry, The Hebrew University of Jerusalem, Jerusalem 91904, Israel.
| | | | | |
Collapse
|
7
|
Pradhan E, Sato K, Akimov AV. Non-adiabatic molecular dynamics with ΔSCF excited states. JOURNAL OF PHYSICS. CONDENSED MATTER : AN INSTITUTE OF PHYSICS JOURNAL 2018; 30:484002. [PMID: 30407924 DOI: 10.1088/1361-648x/aae864] [Citation(s) in RCA: 22] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/08/2023]
Abstract
Accurate modelling of nonadiabatic transitions and electron-phonon interactions in extended systems is essential for understanding the charge and energy transfer in photovoltaic and photocatalytic materials. The extensive computational costs of the advanced excited state methods have stimulated the development of many approximations to study the nonadiabatic molecular dynamics (NA-MD) in solid-state and molecular materials. In this work, we present a novel ▵SCF-NA-MD methodology that aims to account for electron-hole interactions and electron-phonon back-reaction critical in modelling photoinduced nuclear dynamics. The excited states dynamics is described using the delta self-consistent field (▵SCF) technique within the density functional formalism and the trajectory surface hopping. The technique is implemented in the open-source Libra-X package freely available on the Internet (https://github.com/Quantum-Dynamics-Hub/Libra-X). This work illustrates the general utility of the developed ▵SCF-NA-MD methodology by characterizing the excited state energies and lifetimes, reorganization energies, photoisomerization quantum yields, and by providing the mechanistic details of reactive processes in a number of organic molecules.
Collapse
Affiliation(s)
- Ekadashi Pradhan
- Department of Chemistry, University at Buffalo, The State University of New York, Buffalo, NY 14260-3000, United States of America
| | | | | |
Collapse
|
8
|
Barborini M, Sorella S, Rontani M, Corni S. Angle-resolved photoemission spectroscopy from first-principles quantum Monte Carlo. J Chem Phys 2018; 149:154102. [DOI: 10.1063/1.5038864] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Affiliation(s)
| | - Sandro Sorella
- Scuola Internazionale Superiore di Studi Avanzati (SISSA) and CNR-IOM Democritos National Simulation Center, Via Bonomea 265, 34136 Trieste, Italy
| | | | - Stefano Corni
- CNR-NANO, Via Campi 213/a, 41125 Modena, Italy
- Dipartimento di Scienze Chimiche—Università degli Studi di Padova, Via Marzolo 1, 35131 Padova, Italy
| |
Collapse
|
9
|
Grabarek D, Andruniów T. Initial excited-state relaxation of locked retinal protonated schiff base chromophore. An insight from coupled cluster and multireference perturbation theory calculations. J Comput Chem 2018; 39:1720-1727. [PMID: 29727036 DOI: 10.1002/jcc.25346] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/16/2018] [Revised: 04/06/2018] [Accepted: 04/12/2018] [Indexed: 11/07/2022]
Abstract
The initial S1 excited-state relaxation of retinal protonated Schiff base (RPSB) analog with central C11C12 double bond locked by eight-membered ring (locked-11.8) was investigated by means of multireference perturbation theory methods (XMCQDPT2, XMS-CASPT2, MS-CASPT2) as well as single-reference coupled-cluster CC2 method. The analysis of XMCQDPT2-based geometries reveals rather weak coupling between in-plane and out-of-plane structural evolution and minor energetical relaxation of three locked-11.8 conformers. Therefore, a strong coupling between bonds length inversion and backbone out-of-plane deformation resulting in a very steep S1 energy profile predicted by CASSCF/CASPT2 calculations is in clear contradiction with the reference XMCQDPT2 results. Even though CC2 method predicts good quality ground-state structures, the excited-state structures display more advanced torsional deformation leading to ca. 0.2 eV exaggerated energy relaxation and significantly red shifted (0.4-0.7 eV) emission maxima. According to our findings, the initial photoisomerization process in locked-11.8, and possibly in other RPSB analogs, studied fully (both geometries and energies) by multireference perturbation theory may be somewhat slower than predicted by CASSCF/CASPT2 or CC2 methods. © 2018 Wiley Periodicals, Inc.
Collapse
Affiliation(s)
- 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
|
10
|
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
|
11
|
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
|
12
|
Mussard B, Coccia E, Assaraf R, Otten M, Umrigar CJ, Toulouse J. Time-Dependent Linear-Response Variational Monte Carlo. ADVANCES IN QUANTUM CHEMISTRY 2018. [DOI: 10.1016/bs.aiq.2017.05.005] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/05/2022]
|
13
|
Mato J, Gordon MS. A general spin-complete spin-flip configuration interaction method. Phys Chem Chem Phys 2018; 20:2615-2626. [DOI: 10.1039/c7cp06837a] [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
A new spin-complete spin-flip configuration interaction (SF-CI) method was developed using the ORMAS-CI algorithm, along with a perturbative correction for capturing dynamic and non-dynamic correlation.
Collapse
Affiliation(s)
- Joani Mato
- Department of Chemistry
- Iowa State University
- Ames
- USA
| | | |
Collapse
|
14
|
Isomerization of the RPSB chromophore in the gas phase along the torsional pathways using QTAIM. Chem Phys Lett 2017. [DOI: 10.1016/j.cplett.2017.07.066] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
|
15
|
Tenti L, Giner E, Malrieu JP, Angeli C. Strongly localized approaches for delocalized systems. I. Ground state of linear polyenes. COMPUT THEOR CHEM 2017. [DOI: 10.1016/j.comptc.2017.01.021] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
|
16
|
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
|
17
|
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
|
18
|
Introduction to the Variational Monte Carlo Method in Quantum Chemistry and Physics. VARIATIONAL METHODS IN MOLECULAR MODELING 2017. [DOI: 10.1007/978-981-10-2502-0_10] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 02/10/2023]
|
19
|
Barborini M, Guidoni L. Geometries of low spin states of multi-centre transition metal complexes through extended broken symmetry variational Monte Carlo. J Chem Phys 2016; 145:124107. [DOI: 10.1063/1.4963015] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Matteo Barborini
- Dipartimento di Ingegneria, Scienze dell’Informazione e Matematica, Università degli studi dell’Aquila, Via Vetoio 2, 67100 Coppito, L’Aquila, Italy
- Dipartimento di Scienze Fisiche e Chimiche, Università degli studi dell’Aquila, Via Vetoio 2, 67100 Coppito, L’Aquila, Italy
| | - Leonardo Guidoni
- Dipartimento di Scienze Fisiche e Chimiche, Università degli studi dell’Aquila, Via Vetoio 2, 67100 Coppito, L’Aquila, Italy
| |
Collapse
|
20
|
Grabarek D, Walczak E, Andruniów T. Assessing the Accuracy of Various Ab Initio Methods for Geometries and Excitation Energies of Retinal Chromophore Minimal Model by Comparison with CASPT3 Results. J Chem Theory Comput 2016; 12:2346-56. [DOI: 10.1021/acs.jctc.6b00108] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Dawid Grabarek
- Advanced Materials Engineering
and Modelling Group, Wroclaw University of Technology, Wyb. Wyspianskiego
27, 50-370 Wroclaw, Poland
| | - Elżbieta Walczak
- Advanced Materials Engineering
and Modelling Group, Wroclaw University of Technology, Wyb. Wyspianskiego
27, 50-370 Wroclaw, Poland
| | - Tadeusz Andruniów
- Advanced Materials Engineering
and Modelling Group, Wroclaw University of Technology, Wyb. Wyspianskiego
27, 50-370 Wroclaw, Poland
| |
Collapse
|
21
|
Pastore M, De Angelis F, Angeli C. Optical absorption spectrum of the N3 solar cell sensitizer by second-order multireference perturbation theory. Theor Chem Acc 2016. [DOI: 10.1007/s00214-016-1868-y] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/22/2022]
|
22
|
Guareschi R, Zulfikri H, Daday C, Floris FM, Amovilli C, Mennucci B, Filippi C. Introducing QMC/MMpol: Quantum Monte Carlo in Polarizable Force Fields for Excited States. J Chem Theory Comput 2016; 12:1674-83. [DOI: 10.1021/acs.jctc.6b00044] [Citation(s) in RCA: 24] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/23/2022]
Affiliation(s)
- Riccardo Guareschi
- MESA+
Institute for Nanotechnology, University of Twente, P.O. Box 217, 7500 AE Enschede, The Netherlands
| | - Habiburrahman Zulfikri
- MESA+
Institute for Nanotechnology, University of Twente, P.O. Box 217, 7500 AE Enschede, The Netherlands
| | - Csaba Daday
- MESA+
Institute for Nanotechnology, University of Twente, P.O. Box 217, 7500 AE Enschede, The Netherlands
| | - Franca Maria Floris
- Dipartimento
di Chimica e Chimica Industriale, Università di Pisa, Via Giuseppe
Moruzzi 13, 56124 Pisa, Italy
| | - Claudio Amovilli
- Dipartimento
di Chimica e Chimica Industriale, Università di Pisa, Via Giuseppe
Moruzzi 13, 56124 Pisa, Italy
| | - Benedetta Mennucci
- Dipartimento
di Chimica e Chimica Industriale, Università di Pisa, Via Giuseppe
Moruzzi 13, 56124 Pisa, Italy
| | - Claudia Filippi
- MESA+
Institute for Nanotechnology, University of Twente, P.O. Box 217, 7500 AE Enschede, The Netherlands
| |
Collapse
|
23
|
Kuisma MJ, Lundin AM, Moth-Poulsen K, Hyldgaard P, Erhart P. Comparative Ab-Initio Study of Substituted Norbornadiene-Quadricyclane Compounds for Solar Thermal Storage. THE JOURNAL OF PHYSICAL CHEMISTRY. C, NANOMATERIALS AND INTERFACES 2016; 120:3635-3645. [PMID: 26966476 PMCID: PMC4780837 DOI: 10.1021/acs.jpcc.5b11489] [Citation(s) in RCA: 25] [Impact Index Per Article: 3.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 11/24/2015] [Revised: 01/25/2016] [Indexed: 05/05/2023]
Abstract
Molecular photoswitches that are capable of storing solar energy, so-called molecular solar thermal storage systems, are interesting candidates for future renewable energy applications. In this context, substituted norbornadiene-quadricyclane systems have received renewed interest due to recent advances in their synthesis. The optical, thermodynamic, and kinetic properties of these systems can vary dramatically depending on the chosen substituents. The molecular design of optimal compounds therefore requires a detailed understanding of the effect of individual substituents as well as their interplay. Here, we model absorption spectra, potential energy storage, and thermal barriers for back-conversion of several substituted systems using both single-reference (density functional theory using PBE, B3LYP, CAM-B3LYP, M06, M06-2x, and M06-L functionals as well as MP2 calculations) and multireference methods (complete active space techniques). Already the diaryl substituted compound displays a strong red-shift compared to the unsubstituted system, which is shown to result from the extension of the conjugated π-system upon substitution. Using specific donor/acceptor groups gives rise to a further albeit relatively smaller red-shift. The calculated storage energy is found to be rather insensitive to the specific substituents, although solvent effects are likely to be important and require further study. The barrier for thermal back-conversion exhibits strong multireference character and as a result is noticeably correlated with the red-shift. Two possible reaction paths for the thermal back-conversion of diaryl substituted quadricyclane are identified and it is shown that among the compounds considered the path via the acceptor side is systematically favored. Finally, the present study establishes the basis for high-throughput screening of norbornadiene-quadricyclane compounds as it provides guidelines for the level of accuracy that can be expected for key properties from several different techniques.
Collapse
Affiliation(s)
- Mikael J. Kuisma
- Department
of Physics, Chalmers University of Technology, SE-412 96 Gothenburg, Sweden
- E-mail:
| | - Angelica M. Lundin
- Department
of Chemistry and Chemical Engineering, Chalmers
University of Technology, SE-412 96 Gothenburg, Sweden
| | - Kasper Moth-Poulsen
- Department
of Chemistry and Chemical Engineering, Chalmers
University of Technology, SE-412 96 Gothenburg, Sweden
| | - Per Hyldgaard
- Department
of Microtechnology and Nanoscience, Chalmers
University of Technology, SE-412 96 Gothenburg, Sweden
| | - Paul Erhart
- Department
of Physics, Chalmers University of Technology, SE-412 96 Gothenburg, Sweden
- E-mail:
| |
Collapse
|
24
|
Nile blue and Nile red optical properties predicted by TD-DFT and CASPT2 methods: static and dynamic solvent effects. Theor Chem Acc 2016. [DOI: 10.1007/s00214-016-1814-z] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/22/2022]
|
25
|
Cheng C, Kamiya M, Uchida Y, Hayashi S. Molecular Mechanism of Wide Photoabsorption Spectral Shifts of Color Variants of Human Cellular Retinol Binding Protein II. J Am Chem Soc 2015; 137:13362-70. [PMID: 26418582 DOI: 10.1021/jacs.5b08316] [Citation(s) in RCA: 25] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Color variants of human cellular retinol binding protein II (hCRBPII) created by protein engineering were recently shown to exhibit anomalously wide photoabsorption spectral shifts over ∼200 nm across the visible region. The remarkable phenomenon provides a unique opportunity to gain insight into the molecular basis of the color tuning of retinal binding proteins for understanding of color vision as well as for engineering of novel color variants of retinal binding photoreceptor proteins employed in optogenetics. Here, we report a theoretical investigation of the molecular mechanism underlying the anomalously wide spectral shifts of the color variants of hCRBPII. Computational modeling of the color variants with hybrid molecular simulations of free energy geometry optimization succeeded in reproducing the experimentally observed wide spectral shifts, and revealed that protein flexibility, through which the active site structure of the protein and bound water molecules is altered by remote mutations, plays a significant role in inducing the large spectral shifts.
Collapse
Affiliation(s)
- Cheng Cheng
- Department of Chemistry, Graduate School of Science, Kyoto University , Kyoto 606-8502, Japan
| | - Motoshi Kamiya
- Department of Chemistry, Graduate School of Science, Kyoto University , Kyoto 606-8502, Japan
| | - Yoshihiro Uchida
- Department of Chemistry, Graduate School of Science, Kyoto University , Kyoto 606-8502, Japan
| | - Shigehiko Hayashi
- Department of Chemistry, Graduate School of Science, Kyoto University , Kyoto 606-8502, Japan
| |
Collapse
|
26
|
Barborini M, Guidoni L. Ground State Geometries of Polyacetylene Chains from Many-Particle Quantum Mechanics. J Chem Theory Comput 2015; 11:4109-18. [PMID: 26405437 PMCID: PMC4570661 DOI: 10.1021/acs.jctc.5b00427] [Citation(s) in RCA: 17] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/08/2015] [Indexed: 11/28/2022]
Abstract
Due to the crucial role played by electron correlation, the accurate determination of ground state geometries of π-conjugated molecules is still a challenge for many quantum chemistry methods. Because of the high parallelism of the algorithms and their explicit treatment of electron correlation effects, Quantum Monte Carlo calculations can offer an accurate and reliable description of the electronic states and of the geometries of such systems, competing with traditional quantum chemistry approaches. Here, we report the structural properties of polyacetylene chains H-(C₂H₂)(N)-H up to N = 12 acetylene units, by means of Variational Monte Carlo (VMC) calculations based on the multi-determinant Jastrow Antisymmetrized Geminal Power (JAGP) wave function. This compact ansatz can provide for such systems an accurate description of the dynamical electronic correlation as recently detailed for the 1,3-butadiene molecule [J. Chem. Theory Comput. 2015 11 (2), 508-517]. The calculated Bond Length Alternation (BLA), namely the difference between the single and double carbon bonds, extrapolates, for N → ∞, to a value of 0.0910(7) Å, compatible with the experimental data. An accurate analysis was able to distinguish between the influence of the multi-determinantal AGP expansion and of the Jastrow factor on the geometrical properties of the fragments. Our size-extensive and self-interaction-free results provide new and accurate ab initio references for the structures of the ground state of polyenes.
Collapse
Affiliation(s)
- Matteo Barborini
- Dipartimento
di Ingegneria e Scienze dell’Informazione e Matematica, Università degli studi dell’Aquila, Via Vetoio (Coppito), 67100 L’Aquila, Italy
- Dipartimento
di Scienze Fisiche e Chimiche, Università
degli studi dell’Aquila, Via Vetoio (Coppito), 67100 L’Aquila, Italy
| | - Leonardo Guidoni
- Dipartimento
di Scienze Fisiche e Chimiche, Università
degli studi dell’Aquila, Via Vetoio (Coppito), 67100 L’Aquila, Italy
| |
Collapse
|
27
|
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
|
28
|
Valsson O, Filippi C, Casida ME. Regarding the use and misuse of retinal protonated Schiff base photochemistry as a test case for time-dependent density-functional theory. J Chem Phys 2015; 142:144104. [DOI: 10.1063/1.4916354] [Citation(s) in RCA: 14] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
29
|
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
|
30
|
Guareschi R, Floris FM, Amovilli C, Filippi C. Solvent Effects on Excited-State Structures: A Quantum Monte Carlo and Density Functional Study. J Chem Theory Comput 2014; 10:5528-37. [DOI: 10.1021/ct500723s] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/27/2023]
Affiliation(s)
- Riccardo Guareschi
- MESA+
Institute for Nanotechnology, University of Twente, P.O. Box 217, 7500 AE Enschede, The Netherlands
| | - Franca Maria Floris
- Dipartimento
di Chimica e Chimica Industriale, Università di Pisa, Via Giuseppe
Moruzzi 3, 56124 Pisa, Italy
| | - Claudio Amovilli
- Dipartimento
di Chimica e Chimica Industriale, Università di Pisa, Via Giuseppe
Moruzzi 3, 56124 Pisa, Italy
| | - Claudia Filippi
- MESA+
Institute for Nanotechnology, University of Twente, P.O. Box 217, 7500 AE Enschede, The Netherlands
| |
Collapse
|
31
|
Di Donato M, Segado Centellas M, Lapini A, Lima M, Avila F, Santoro F, Cappelli C, Righini R. Combination of transient 2D-IR experiments and ab initio computations sheds light on the formation of the charge-transfer state in photoexcited carbonyl carotenoids. J Phys Chem B 2014; 118:9613-30. [PMID: 25050938 DOI: 10.1021/jp505473j] [Citation(s) in RCA: 14] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
The excited state dynamics of carbonyl carotenoids is very complex because of the coupling of single- and doubly excited states and the possible involvement of intramolecular charge-transfer (ICT) states. In this contribution we employ ultrafast infrared spectroscopy and theoretical computations to investigate the relaxation dynamics of trans-8'-apo-β-carotenal occurring on the picosecond time scale, after excitation in the S2 state. In a (slightly) polar solvent like chloroform, one-dimensional (T1D-IR) and two-dimensional (T2D-IR) transient infrared spectroscopy reveal spectral components with characteristic frequencies and lifetimes that are not observed in nonpolar solvents (cyclohexane). Combining experimental evidence with an analysis of CASPT2//CASSCF ground and excited state minima and energy profiles, complemented with TDDFT calculations in gas phase and in solvent, we propose a photochemical decay mechanism for this system where only the bright single-excited 1Bu(+) and the dark double-excited 2Ag(-) states are involved. Specifically, the initially populated 1Bu(+) relaxes toward 2Ag(-) in 200 fs. In a nonpolar solvent 2Ag(-) decays to the ground state (GS) in 25 ps. In polar solvents, distortions along twisting modes of the chain promote a repopulation of the 1Bu(+) state which then quickly relaxes to the GS (18 ps in chloroform). The 1Bu(+) state has a high electric dipole and is the main contributor to the charge-transfer state involved in the dynamics in polar solvents. The 2Ag(-) → 1Bu(+) population transfer is evidenced by a cross peak on the T2D-IR map revealing that the motions along the same stretching of the conjugated chain on the 2Ag(-) and 1Bu(+) states are coupled.
Collapse
Affiliation(s)
- Mariangela Di Donato
- LENS (European Laboratory for Nonlinear Spectroscopy) via N. Carrara 1, 50019 Sesto Fiorentino (FI), Italy
| | | | | | | | | | | | | | | |
Collapse
|
32
|
Varsano D, Coccia E, Pulci O, Conte AM, Guidoni L. Ground state structures and electronic excitations of biological chromophores at Quantum Monte Carlo/Many Body Green’s Function Theory level. COMPUT THEOR CHEM 2014. [DOI: 10.1016/j.comptc.2014.03.011] [Citation(s) in RCA: 24] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/18/2022]
|
33
|
Varsano D, Barborini M, Guidoni L. Kohn-Sham orbitals and potentials from quantum Monte Carlo molecular densities. J Chem Phys 2014; 140:054102. [DOI: 10.1063/1.4863213] [Citation(s) in RCA: 8] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
34
|
Guareschi R, Filippi C. Ground- and Excited-State Geometry Optimization of Small Organic Molecules with Quantum Monte Carlo. J Chem Theory Comput 2013; 9:5513-25. [DOI: 10.1021/ct400876y] [Citation(s) in RCA: 47] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/26/2023]
Affiliation(s)
- Riccardo Guareschi
- MESA+ Institute for Nanotechnology, University of Twente,
P.O. Box 217, 7500 AE Enschede, The Netherlands
| | - Claudia Filippi
- MESA+ Institute for Nanotechnology, University of Twente,
P.O. Box 217, 7500 AE Enschede, The Netherlands
| |
Collapse
|
35
|
Zhou P, Liu J, Han K, He G. The photoisomerization of 11-cis-retinal protonated schiff base in gas phase: Insight from spin-flip density functional theory. J Comput Chem 2013; 35:109-20. [DOI: 10.1002/jcc.23463] [Citation(s) in RCA: 33] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/21/2013] [Revised: 09/12/2013] [Accepted: 09/24/2013] [Indexed: 11/10/2022]
Affiliation(s)
- Panwang Zhou
- State Key Laboratory of Molecular Reaction Dynamics; Dalian Institute of Chemical Physics; Chinese Academy of Sciences; Dalian 116023 Liaoning People's Republic of China
| | - Jianyong Liu
- State Key Laboratory of Molecular Reaction Dynamics; Dalian Institute of Chemical Physics; Chinese Academy of Sciences; Dalian 116023 Liaoning People's Republic of China
| | - Keli Han
- State Key Laboratory of Molecular Reaction Dynamics; Dalian Institute of Chemical Physics; Chinese Academy of Sciences; Dalian 116023 Liaoning People's Republic of China
| | - Guozhong He
- State Key Laboratory of Molecular Reaction Dynamics; Dalian Institute of Chemical Physics; Chinese Academy of Sciences; Dalian 116023 Liaoning People's Republic of China
| |
Collapse
|
36
|
Walczak E, Szefczyk B, Andruniów T. Geometries and Vertical Excitation Energies in Retinal Analogues Resolved at the CASPT2 Level of Theory: Critical Assessment of the Performance of CASSCF, CC2, and DFT Methods. J Chem Theory Comput 2013; 9:4915-27. [DOI: 10.1021/ct400423u] [Citation(s) in RCA: 25] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Elżbieta Walczak
- Wroclaw University of Technology, Institute of Physical & Theoretical Chemistry, Wyb. Wyspianskiego 27, 50-370 Wroclaw, Poland
| | - Borys Szefczyk
- Wroclaw University of Technology, Institute of Physical & Theoretical Chemistry, Wyb. Wyspianskiego 27, 50-370 Wroclaw, Poland
| | - Tadeusz Andruniów
- Wroclaw University of Technology, Institute of Physical & Theoretical Chemistry, Wyb. Wyspianskiego 27, 50-370 Wroclaw, Poland
| |
Collapse
|
37
|
Wang W, Geiger JH, Borhan B. The photochemical determinants of color vision: revealing how opsins tune their chromophore's absorption wavelength. Bioessays 2013; 36:65-74. [PMID: 24323922 DOI: 10.1002/bies.201300094] [Citation(s) in RCA: 28] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/08/2023]
Abstract
The evolution of a variety of important chromophore-dependent biological processes, including microbial light sensing and mammalian color vision, relies on protein modifications that alter the spectral characteristics of a bound chromophore. Three different color opsins share the same chromophore, but have three distinct absorptions that together cover the entire visible spectrum, giving rise to trichromatic vision. The influence of opsins on the absorbance of the chromophore has been studied through methods such as model compounds, opsin mutagenesis, and computational modeling. The recent development of rhodopsin mimic that uses small soluble proteins to recapitulate the binding and wavelength tuning of the native opsins provides a new platform for studying protein-regulated spectral tuning. The ability to achieve far-red shifted absorption in the rhodopsin mimic system was attributed to a combination of the lack of a counteranion proximal to the iminium, and a uniformly neutral electrostatic environment surrounding the chromophore.
Collapse
Affiliation(s)
- Wenjing Wang
- Department of Chemistry, Massachusetts Institute of Technology, Cambridge, MA, USA
| | | | | |
Collapse
|
38
|
Malrieu JP, Caballol R, Calzado CJ, de Graaf C, Guihéry N. Magnetic interactions in molecules and highly correlated materials: physical content, analytical derivation, and rigorous extraction of magnetic Hamiltonians. Chem Rev 2013; 114:429-92. [PMID: 24102410 DOI: 10.1021/cr300500z] [Citation(s) in RCA: 282] [Impact Index Per Article: 25.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/07/2023]
Affiliation(s)
- Jean Paul Malrieu
- Laboratoire de Chimie et Physique Quantiques, Université de Toulouse 3 , 118 route de Narbonne, 31062 Toulouse, France
| | | | | | | | | |
Collapse
|
39
|
Huntress MM, Gozem S, Malley KR, Jailaubekov AE, Vasileiou C, Vengris M, Geiger JH, Borhan B, Schapiro I, Larsen DS, Olivucci M. Toward an Understanding of the Retinal Chromophore in Rhodopsin Mimics. J Phys Chem B 2013; 117:10053-70. [DOI: 10.1021/jp305935t] [Citation(s) in RCA: 34] [Impact Index Per Article: 3.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Mark M. Huntress
- Department of Chemistry, Bowling Green State University, Bowling Green, Ohio
43402, United States
| | - Samer Gozem
- Department of Chemistry, Bowling Green State University, Bowling Green, Ohio
43402, United States
| | - Konstantin R. Malley
- Department
of Chemistry, University of California Davis, One Shields Avenure,
Davis, California 95616, United States
| | - Askat E. Jailaubekov
- Department
of Chemistry, University of California Davis, One Shields Avenure,
Davis, California 95616, United States
| | - Chrysoula Vasileiou
- Department of Chemistry, Michigan State University, Lansing, Michigan 48824,
United States
| | - Mikas Vengris
- Department
of Chemistry, University of California Davis, One Shields Avenure,
Davis, California 95616, United States
- Faculty of
Physics, Vilnius University, Sauletekio
10 LT10223 Vilnius,
Lithuania
| | - James H. Geiger
- Department of Chemistry, Michigan State University, Lansing, Michigan 48824,
United States
| | - Babak Borhan
- Department of Chemistry, Michigan State University, Lansing, Michigan 48824,
United States
| | - Igor Schapiro
- Department of Chemistry, Bowling Green State University, Bowling Green, Ohio
43402, United States
| | - Delmar S. Larsen
- Department
of Chemistry, University of California Davis, One Shields Avenure,
Davis, California 95616, United States
| | - Massimo Olivucci
- Department of Chemistry, Bowling Green State University, Bowling Green, Ohio
43402, United States
| |
Collapse
|
40
|
Huix-Rotllant M, Filatov M, Gozem S, Schapiro I, Olivucci M, Ferré N. Assessment of Density Functional Theory for Describing the Correlation Effects on the Ground and Excited State Potential Energy Surfaces of a Retinal Chromophore Model. J Chem Theory Comput 2013; 9:3917-32. [PMID: 26592387 DOI: 10.1021/ct4003465] [Citation(s) in RCA: 78] [Impact Index Per Article: 7.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/14/2023]
Abstract
In the quest for a cost-effective level of theory able to describe a large portion of the ground and excited potential energy surfaces of large chromophores, promising approaches are rooted in various approximations to the exact density functional theory (DFT). In the present work, we investigate how generalized Kohn-Sham DFT (GKS-DFT), time-dependent DFT (TDDFT), and spin-restricted ensemble-DFT (REKS) methods perform along three important paths characterizing a model retinal chromophore (the penta-2,4-dieniminium cation) in a region of near-degeneracy (close to a conical intersection) with respect to reference high-level multiconfigurational wave function methods. If GKS-DFT correctly describes the closed-shell charge transfer state, only TDDFT and REKS approaches give access to the open-shell diradical, one which sometimes corresponds to the electronic ground state. It is demonstrated that the main drawback of the usual DFT-based methods lies in the absence of interactions between the charge transfer and the diradicaloid configurations. Hence, we test a new computational scheme based on the State-averaged REKS (SA-REKS) approach, which explicitly includes these interactions into account. The State-Interaction SA-REKS (SI-SA-REKS) method significantly improves on the REKS and the SA-REKS results for the target system. The similarities and differences between DFT and wave function-based approaches are analyzed according to (1) the active space dimensions of the wave function-based methods and (2) the relative electronegativities of the allyl and protonated Schiff base moieties.
Collapse
Affiliation(s)
- Miquel Huix-Rotllant
- Aix-Marseille Université, CNRS, Institut de Chimie Radicalaire , 13397 Marseille Cedex 20, France
| | - Michael Filatov
- Institut für Physicalische und Theoretische Chemie, Universität Bonn , Beringstr. 4, 53115 Bonn, Germany
| | - Samer Gozem
- Department of Chemistry, Bowling Green State University , Bowling Green, Ohio 43402, United States
| | - Igor Schapiro
- Department of Chemistry, Bowling Green State University , Bowling Green, Ohio 43402, United States.,Max Planck Institute for Chemical Energy Conversion , Stiftstr. 34 - 36, Mülheim an der Ruhr, Germany
| | - Massimo Olivucci
- Department of Chemistry, Bowling Green State University , Bowling Green, Ohio 43402, United States.,Dipartimento di Chimica, Università di Siena , Siena, Italy
| | - Nicolas Ferré
- Aix-Marseille Université, CNRS, Institut de Chimie Radicalaire , 13397 Marseille Cedex 20, France
| |
Collapse
|
41
|
Schapiro I, Sivalingam K, Neese F. Assessment of n-Electron Valence State Perturbation Theory for Vertical Excitation Energies. J Chem Theory Comput 2013; 9:3567-80. [DOI: 10.1021/ct400136y] [Citation(s) in RCA: 102] [Impact Index Per Article: 9.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/08/2023]
Affiliation(s)
- Igor Schapiro
- Max Planck Institute for Chemical Energy Conversion, Stiftstr. 34-36, 45470 Mülheim an
der Ruhr, Germany
| | - Kantharuban Sivalingam
- Max Planck Institute for Chemical Energy Conversion, Stiftstr. 34-36, 45470 Mülheim an
der Ruhr, Germany
| | - Frank Neese
- Max Planck Institute for Chemical Energy Conversion, Stiftstr. 34-36, 45470 Mülheim an
der Ruhr, Germany
| |
Collapse
|
42
|
Valsson O, Campomanes P, Tavernelli I, Rothlisberger U, Filippi C. Rhodopsin Absorption from First Principles: Bypassing Common Pitfalls. J Chem Theory Comput 2013; 9:2441-54. [PMID: 26583734 DOI: 10.1021/ct3010408] [Citation(s) in RCA: 71] [Impact Index Per Article: 6.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Bovine rhodopsin is the most extensively studied retinal protein and is considered the prototype of this important class of photosensitive biosystems involved in the process of vision. Many theoretical investigations have attempted to elucidate the role of the protein matrix in modulating the absorption of retinal chromophore in rhodopsin, but, while generally agreeing in predicting the correct location of the absorption maximum, they often reached contradicting conclusions on how the environment tunes the spectrum. To address this controversial issue, we combine here a thorough structural and dynamical characterization of rhodopsin with a careful validation of its excited-state properties via the use of a wide range of state-of-the-art quantum chemical approaches including various flavors of time-dependent density functional theory (TDDFT), different multireference perturbative schemes (CASPT2 and NEVPT2), and quantum Monte Carlo (QMC) methods. Through extensive quantum mechanical/molecular mechanical (QM/MM) molecular dynamics simulations, we obtain a comprehensive structural description of the chromophore-protein system and sample a wide range of thermally accessible configurations. We show that, in order to obtain reliable excitation properties, it is crucial to employ a sufficient number of representative configurations of the system. In fact, the common use of a single, ad hoc structure can easily lead to an incorrect model and an agreement with experimental absorption spectra due to cancelation of errors. Finally, we show that, to properly account for polarization effects on the chromophore and to quench the large blue-shift induced by the counterion on the excitation energies, it is necessary to adopt an enhanced description of the protein environment as given by a large quantum region including as many as 250 atoms.
Collapse
Affiliation(s)
- Omar Valsson
- MESA+ Institute for Nanotechnology, University of Twente , P.O. Box 217, 7500 AE Enschede, The Netherlands
| | - Pablo Campomanes
- Laboratory of Computational Chemistry and Biochemistry, Ecole Polytechnique Fédérale de Lausanne (EPFL) , CH-1015 Lausanne, Switzerland
| | - Ivano Tavernelli
- Laboratory of Computational Chemistry and Biochemistry, Ecole Polytechnique Fédérale de Lausanne (EPFL) , CH-1015 Lausanne, Switzerland
| | - Ursula Rothlisberger
- Laboratory of Computational Chemistry and Biochemistry, Ecole Polytechnique Fédérale de Lausanne (EPFL) , CH-1015 Lausanne, Switzerland
| | - Claudia Filippi
- MESA+ Institute for Nanotechnology, University of Twente , P.O. Box 217, 7500 AE Enschede, The Netherlands
| |
Collapse
|
43
|
Muñoz-Losa A, Fdez. Galván I, Aguilar MA, Martín ME. Simultaneous Solvent and Counterion Effects on the Absorption Properties of a Model of the Rhodopsin Chromophore. J Chem Theory Comput 2013; 9:1548-56. [DOI: 10.1021/ct301090v] [Citation(s) in RCA: 14] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Aurora Muñoz-Losa
- Química
Física, Universidad de Extremadura, Avda. de Elvas s/n, Badajoz,
06071, Spain
| | - Ignacio Fdez. Galván
- Química
Física, Universidad de Extremadura, Avda. de Elvas s/n, Badajoz,
06071, Spain
| | - Manuel A. Aguilar
- Química
Física, Universidad de Extremadura, Avda. de Elvas s/n, Badajoz,
06071, Spain
| | - M. Elena Martín
- Química
Física, Universidad de Extremadura, Avda. de Elvas s/n, Badajoz,
06071, Spain
| |
Collapse
|
44
|
Coccia E, Varsano D, Guidoni L. Protein Field Effect on the Dark State of 11- cis Retinal in Rhodopsin by Quantum Monte Carlo/Molecular Mechanics. J Chem Theory Comput 2013; 9:8-12. [PMID: 24611033 PMCID: PMC3943175 DOI: 10.1021/ct3007502] [Citation(s) in RCA: 30] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
The accurate determination of the geometrical details of the dark state of 11-cis retinal in rhodopsin represents a fundamental step for the rationalization of the protein role in the optical spectral tuning in the vision mechanism. We have calculated geometries of the full retinal protonated Schiff base chromophore in the gas phase and in the protein environment using the correlated variational Monte Carlo method. The bond length alternation of the conjugated carbon chain of the chromophore in the gas phase shows a significant reduction when moving from the β-ionone ring to the nitrogen, whereas, as expected, the protein environment reduces the electronic conjugation. The proposed dark state structure is fully compatible with solid-state NMR data reported by Carravetta et al. [J. Am. Chem. Soc. 2004, 126, 3948-3953]. TDDFT/B3LYP calculations on such geometries show a blue opsin shift of 0.28 and 0.24 eV induced by the protein for S1 and S2 states, consistently with literature spectroscopic data. The effect of the geometrical distortion alone is a red shift of 0.21 and 0.16 eV with respect to the optimized gas phase chromophore. Our results open new perspectives for the study of the properties of chromophores in their biological environment using correlated methods.
Collapse
Affiliation(s)
- Emanuele Coccia
- Dipartimento di Scienze Fisiche e Chimiche, Università degli Studi dell’Aquila, via Vetoio, 67100, L’Aquila, Italy
| | - Daniele Varsano
- Dipartimento di Fisica, “Sapienza” - Università di Roma, piazzale Aldo Moro 5, 00185, Rome, Italy
| | - Leonardo Guidoni
- Dipartimento di Scienze Fisiche e Chimiche, Università degli Studi dell’Aquila, via Vetoio, 67100, L’Aquila, Italy
| |
Collapse
|
45
|
Barborini M, Guidoni L. Reaction pathways by quantum Monte Carlo: insight on the torsion barrier of 1,3-butadiene, and the conrotatory ring opening of cyclobutene. J Chem Phys 2012; 137:224309. [PMID: 23249005 PMCID: PMC3925821 DOI: 10.1063/1.4769791] [Citation(s) in RCA: 32] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Quantum Monte Carlo (QMC) methods are used to investigate the intramolecular reaction pathways of 1,3-butadiene. The ground state geometries of the three conformers s-trans, s-cis, and gauche, as well as the cyclobutene structure are fully optimised at the variational Monte Carlo (VMC) level, obtaining an excellent agreement with the experimental results and other quantum chemistry high level calculations. Transition state geometries are also estimated at the VMC level for the s-trans to gauche torsion barrier of 1,3-butadiene and for the conrotatory ring opening of cyclobutene to the gauche-1,3-butadiene conformer. The energies of the conformers and the reaction barriers are calculated at both variational and diffusional Monte Carlo levels providing a precise picture of the potential energy surface of 1,3-butadiene and supporting one of the two model profiles recently obtained by Raman spectroscopy [Boopalachandran et al., J. Phys. Chem. A 115, 8920 (2011)]. Considering the good scaling of QMC techniques with the system's size, our results also demonstrate how variational Monte Carlo calculations can be applied in the future to properly investigate the reaction pathways of large and correlated molecular systems.
Collapse
Affiliation(s)
- Matteo Barborini
- Dipartimento di Matematica Pura ed Applicata, Università degli studi dell'Aquila, via Vetoio (Coppito 2), 67100 L'Aquila, Italy
| | | |
Collapse
|
46
|
Gozem S, Huntress M, Schapiro I, Lindh R, Granovsky AA, Angeli C, Olivucci M. Dynamic Electron Correlation Effects on the Ground State Potential Energy Surface of a Retinal Chromophore Model. J Chem Theory Comput 2012; 8:4069-80. [DOI: 10.1021/ct3003139] [Citation(s) in RCA: 129] [Impact Index Per Article: 10.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Samer Gozem
- Department of Chemistry, Bowling
Green State University, Bowling Green, Ohio 43403, United States
| | - Mark Huntress
- Department of Chemistry, Bowling
Green State University, Bowling Green, Ohio 43403, United States
| | - Igor Schapiro
- Department of Chemistry, Bowling
Green State University, Bowling Green, Ohio 43403, United States
| | - Roland Lindh
- Department of Chemistry −
Ångström, the Theoretical Chemistry Programme, POB 518,
SE-751 20 Uppsala, Sweden
| | | | - Celestino Angeli
- Dipartimento di
Chimica, Università
di Ferrara, via Borsari 46, I-44121 Ferrara, Italy
| | - Massimo Olivucci
- Department of Chemistry, Bowling
Green State University, Bowling Green, Ohio 43403, United States
- Dipartimento di Chimica, Università
di Siena, via De Gasperi 2, I-53100 Siena, Italy
| |
Collapse
|