1
|
Seki S, Yoshida K, Sugisaki M, Yamano N, Fujii R. Characterization of the Ultraviolet-B Absorption Band of Carotenoids Using Solvent-dependent Shifts in Steady-State and Transient Absorption Spectra. J Phys Chem B 2024; 128:5623-5629. [PMID: 38833602 DOI: 10.1021/acs.jpcb.4c02212] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/06/2024]
Abstract
The versatile functions of carotenoids in biological systems are associated with the extended π-electron conjugation system. Strong visible absorption resulting from the optically allowed S2 (1Bu+) state and the low-lying optically forbidden S1 (2Ag-) state examined. Carotenoids also exhibit an absorption band in the ultraviolet-B region; however, the origin of this band (hereafter referred to as Suv state) is not well characterized. The Suv state is a candidate for the destination level of the well-known S1 → Sn transient absorption; however, an obvious energy mismatch has been observed. In this study, we examined the steady-state and picosecond transient absorption spectra of lycopene in various solvents. The Suv absorption of carotenoids with diverse conjugation lengths was also examined. The dependence of the energies on solvent polarizability and conjugation length revealed that both Suv and Sn are the "second" Bu+ state. The absorption spectrum for lycopene at 200 K revealed an additional vibrational band, which may be the vibrational origin of the S0 → Suv band. Considering the slow vibrational relaxation of the 2Ag- state, the S1 → Sn transition may represent the 2Ag- (v = 1) → 2Bu+ (v = 0) transition, and the energetic contradiction can be resolved.
Collapse
Affiliation(s)
- Soichiro Seki
- Graduate School of Science, Osaka Metropolitan University, 3-3-138 Sugimoto, Sumiyoshi-ku, Osaka 558-8585, Japan
| | - Kazuhiro Yoshida
- Graduate School of Science, Osaka Metropolitan University, 3-3-138 Sugimoto, Sumiyoshi-ku, Osaka 558-8585, Japan
| | - Mitsuru Sugisaki
- Graduate School of Science, Osaka Metropolitan University, 3-3-138 Sugimoto, Sumiyoshi-ku, Osaka 558-8585, Japan
| | - Nami Yamano
- Key Laboratory of Advanced Light Conversion Materials and Biophotonics, Department of Chemistry, Renmin University of China, Beijing 100872, P. R. China
| | - Ritsuko Fujii
- Graduate School of Science, Osaka Metropolitan University, 3-3-138 Sugimoto, Sumiyoshi-ku, Osaka 558-8585, Japan
- Research Center for Artificial Photosynthesis, Osaka Metropolitan University, 3-3-138 Sugimoto, Sumiyoshi-ku, Osaka 558-8585, Japan
| |
Collapse
|
2
|
Arcidiacono A, Cignoni E, Mazzeo P, Cupellini L, Mennucci B. Predicting Solvatochromism of Chromophores in Proteins through QM/MM and Machine Learning. J Phys Chem A 2024; 128:3646-3658. [PMID: 38683801 PMCID: PMC11089512 DOI: 10.1021/acs.jpca.4c00249] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/12/2024] [Revised: 03/03/2024] [Accepted: 04/01/2024] [Indexed: 05/02/2024]
Abstract
Solvatochromism occurs in both homogeneous solvents and more complex biological environments, such as proteins. While in both cases the solvatochromic effects report on the surroundings of the chromophore, their interpretation in proteins becomes more complicated not only because of structural effects induced by the protein pocket but also because the protein environment is highly anisotropic. This is particularly evident for highly conjugated and flexible molecules such as carotenoids, whose excitation energy is strongly dependent on both the geometry and the electrostatics of the environment. Here, we introduce a machine learning (ML) strategy trained on quantum mechanics/molecular mechanics calculations of geometrical and electrochromic contributions to carotenoids' excitation energies. We employ this strategy to compare solvatochromism in protein and solvent environments. Despite the important specifities of the protein, ML models trained on solvents can faithfully predict excitation energies in the protein environment, demonstrating the robustness of the chosen descriptors.
Collapse
Affiliation(s)
- Amanda Arcidiacono
- Department of Chemistry and Industrial
Chemistry, University of Pisa, Via G. Moruzzi 13, 56124 Pisa, Italy
| | - Edoardo Cignoni
- Department of Chemistry and Industrial
Chemistry, University of Pisa, Via G. Moruzzi 13, 56124 Pisa, Italy
| | - Patrizia Mazzeo
- Department of Chemistry and Industrial
Chemistry, University of Pisa, Via G. Moruzzi 13, 56124 Pisa, Italy
| | - Lorenzo Cupellini
- Department of Chemistry and Industrial
Chemistry, University of Pisa, Via G. Moruzzi 13, 56124 Pisa, Italy
| | - Benedetta Mennucci
- Department of Chemistry and Industrial
Chemistry, University of Pisa, Via G. Moruzzi 13, 56124 Pisa, Italy
| |
Collapse
|
3
|
Barford W. Singlet Fission in Lycopene H-Aggregates. J Phys Chem Lett 2023; 14:9842-9847. [PMID: 37890074 PMCID: PMC10641873 DOI: 10.1021/acs.jpclett.3c02435] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/30/2023] [Revised: 10/09/2023] [Accepted: 10/16/2023] [Indexed: 10/29/2023]
Abstract
A theory of singlet fission (SF) in carotenoid dimers is applied to explain the SF in lycopene H-aggregates observed after high-energy photoexcitation. The explanation proposed here is that a high energy, delocalized bright 1Bu+ state first relaxes and localizes onto a single lycopene monomer. The high-energy intramonomer state then undergoes internal conversion to the 11Bu- state. Once populated, the 11Bu- state allows exothermic bimolecular singlet fission, while its internal conversion to the 21Ag- state is symmetry forbidden. The simulation of SF predicts that the intramonomer triplet-pair state undergoes almost complete population transfer to the intermonomer singlet-pair state within 100 ps. Simultaneously, ZFS interactions begin to partially populate the intermonomer quintet triplet-pair state up to ca. 2 ns, after which hyperfine interactions thermally equilibrate the triplet-pair states, thus forming free single triplets within 50 ns.
Collapse
Affiliation(s)
- William Barford
- Department of Chemistry, Physical and
Theoretical Chemistry Laboratory,University
of Oxford, Oxford, OX1 3QZ, United
Kingdom
| |
Collapse
|
4
|
Menger MFSJ, Köppel H. On the Fluorescence Properties and Nonradiative Transitions in Medium-Sized All-Trans Polyenes. J Phys Chem A 2023; 127:8501-8507. [PMID: 37815131 DOI: 10.1021/acs.jpca.3c03117] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/11/2023]
Abstract
The nonadiabatic photodynamics of all-trans linear polyenes with N = 4-8 conjugated double bonds is studied from an electronic structure perspective. Excitation energies and stationary points for the 1Bu and 2Ag singlet states have been computed by using the state-average complete active space (SA-CASSCF) method and its second-order perturbation theory variant (MS-CASPT2). The dependence of the two low-lying excited states on the "chain length" N has been elucidated. In addition, the 1Bu-2Ag crossing seam has been mapped out in a suitable two-dimensional coordinate space and its minimum within the subspace has been determined. This minimum is found to increase substantially and monotonously in energy with increasing N. This increase is discussed and interpreted in relation to the fluorescence properties of these systems. In particular, it allows to understand the crossover from S1(2Ag) fluorescence for smaller N to S2(1Bu) (or dual) fluorescence for larger N.
Collapse
Affiliation(s)
- Maximilian F S J Menger
- Theoretische Chemie, Physikalisch-Chemisches Institut, University Heidelberg, INF 229, 69120 Heidelberg, Germany
| | - Horst Köppel
- Theoretische Chemie, Physikalisch-Chemisches Institut, University Heidelberg, INF 229, 69120 Heidelberg, Germany
| |
Collapse
|
5
|
Yin BW, Wang JL, Xue PJ, Zhang TS, Xie BB, Shen L, Fang WH. Understanding the Excited-State Relaxation Mechanisms of Xanthophyll Lutein by Multi-configurational Electronic Structure Calculations. J Chem Inf Model 2023; 63:4679-4690. [PMID: 37489739 DOI: 10.1021/acs.jcim.3c00640] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 07/26/2023]
Abstract
The contradictory behaviors in light harvesting and non-photochemical quenching make xanthophyll lutein the most attractive functional molecule in photosynthesis. Despite several theoretical simulations on the spectral properties and excited-state dynamics, the atomic-level photophysical mechanisms need to be further studied and established, especially for an accurate description of geometric and electronic structures of conical intersections for the lowest several electronic states of lutein. In the present work, semiempirical OM2/MRCI and multi-configurational restricted active space self-consistent field methods were performed to optimize the minima and conical intersections in and between the 1Ag-, 2Ag-, 1Bu+, and 1Bu- states. Meanwhile, the relative energies were refined by MS-CASPT2(10,8)/6-31G*, which can reproduce correct electronic state properties as those in the spectroscopic experiments. Based on the above calculation results, we proposed a possible excited-state relaxation mechanism for lutein from its initially populated 1Bu+ state. Once excited to the optically bright 1Bu+ state, the system will propagate along the key reaction coordinate, i.e., the stretching vibration of the conjugated carbon chain. During this period of time, the 1Bu- state will participate in and forms a resonance state between the 1Bu- and 1Bu+ states. Later, the system will rapidly hop to the 2Ag- state via the 1Bu+/2Ag- conical intersection. Finally, the lutein molecule will survive in the 2Ag- state for a relatively long time before it internally converts to the ground state directly or via a twisted S1/S0 conical intersection. Notably, though the photophysical picture may be very different in solvents and proteins, the current theoretical study proposed a promising calculation protocol and also provided many valuable mechanistic insights for lutein and similar carotenoids.
Collapse
Affiliation(s)
- Bo-Wen Yin
- Hangzhou Institute of Advanced Studies, Zhejiang Normal University, 1108 Gengwen Road, Hangzhou 311231 Zhejiang, P. R. China
| | - Jie-Lei Wang
- Hangzhou Institute of Advanced Studies, Zhejiang Normal University, 1108 Gengwen Road, Hangzhou 311231 Zhejiang, P. R. China
| | - Pu-Jie Xue
- Key Laboratory of Theoretical and Computational Photochemistry of Ministry of Education, College of Chemistry, Beijing Normal University, Beijing 100875, P. R. China
| | - Teng-Shuo Zhang
- College of Chemical Engineering, Zhejiang University of Technology, Hangzhou 310014, P. R. China
| | - Bin-Bin Xie
- Hangzhou Institute of Advanced Studies, Zhejiang Normal University, 1108 Gengwen Road, Hangzhou 311231 Zhejiang, P. R. China
| | - Lin Shen
- Key Laboratory of Theoretical and Computational Photochemistry of Ministry of Education, College of Chemistry, Beijing Normal University, Beijing 100875, P. R. China
| | - Wei-Hai Fang
- Key Laboratory of Theoretical and Computational Photochemistry of Ministry of Education, College of Chemistry, Beijing Normal University, Beijing 100875, P. R. China
| |
Collapse
|
6
|
Manawadu D, Valentine DJ, Barford W. Dynamical Simulations of Carotenoid Photoexcited States Using Density Matrix Renormalization Group Techniques. J Phys Chem A 2023; 127:3714-3727. [PMID: 37054397 PMCID: PMC10150368 DOI: 10.1021/acs.jpca.3c00988] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/03/2023]
Abstract
We present a dynamical simulation scheme to model the highly correlated excited state dynamics of linear polyenes. We apply it to investigate the internal conversion processes of carotenoids following their photoexcitation. We use the extended Hubbard-Peierls model, H^UVP, to describe the π-electronic system coupled to nuclear degrees of freedom. This is supplemented by a Hamiltonian, H^ϵ, that explicitly breaks both the particle-hole and two-fold rotation symmetries of idealized carotenoid structures. The electronic degrees of freedom are treated quantum mechanically by solving the time-dependent Schrödinger equation using the adaptive time-dependent DMRG (tDMRG) method, while nuclear dynamics are treated via the Ehrenfest equations of motion. By defining adiabatic excited states as the eigenstates of the full Hamiltonian, H^=H^UVP+H^ϵ, and diabatic excited states as eigenstates of H^UVP, we present a computational framework to monitor the internal conversion process from the initial photoexcited 11Bu+ state to the singlet triplet-pair states of carotenoids. We further incorporate Lanczos-DMRG to the tDMRG-Ehrenfest method to calculate transient absorption spectra from the evolving photoexcited state. We describe in detail the accuracy and convergence criteria for DMRG, and show that this method accurately describes the dynamical processes of carotenoid excited states. We also discuss the effect of the symmetry-breaking term, H^ϵ, on the internal conversion process, and show that its effect on the extent of internal conversion can be described by a Landau-Zener-type transition. This methodological paper is a companion to our more explanatory discussion of carotenoid excited state dynamics in Manawadu, D.; Georges, T. N.; Barford, W. Photoexcited State Dynamics and Singlet Fission in Carotenoids. J. Phys. Chem. A 2023, 127, 1342.
Collapse
Affiliation(s)
- Dilhan Manawadu
- Department of Chemistry, Physical and Theoretical Chemistry Laboratory, University of Oxford, Oxford OX1 3QZ, United Kingdom
- Linacre College, University of Oxford, Oxford OX1 3JA, United Kingdom
| | - Darren J Valentine
- Department of Chemistry, Physical and Theoretical Chemistry Laboratory, University of Oxford, Oxford OX1 3QZ, United Kingdom
- Balliol College, University of Oxford, Oxford OX1 3BJ, United Kingdom
| | - William Barford
- Department of Chemistry, Physical and Theoretical Chemistry Laboratory, University of Oxford, Oxford OX1 3QZ, United Kingdom
| |
Collapse
|
7
|
Dombrowski DR, Schulz T, Kleinschmidt M, Marian CM. R2022: A DFT/MRCI Ansatz with Improved Performance for Double Excitations. J Phys Chem A 2023; 127:2011-2025. [PMID: 36799533 DOI: 10.1021/acs.jpca.2c07951] [Citation(s) in RCA: 5] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/18/2023]
Abstract
A reformulation of the combined density functional theory and multireference configuration interaction method (DFT/MRCI) is presented. Expressions for ab initio matrix elements are used to derive correction terms for a new effective Hamiltonian. On the example of diatomic carbon, the correction terms are derived, focusing on the doubly excited 1Δg state, which was problematic in previous formulations of the method, as were double excitations in general. The derivation shows that a splitting of the parameters for intra- and interorbital interactions is necessary for a concise description of the underlying physics. Results for 1La and 1Lb states in polyacenes and 1Au and 1Ag states in mini-β-carotenoids suggest that the presented formulation is superior to former effective Hamiltonians. Furthermore, statistical analysis reveals that all the benefits of the previous DFT/MRCI Hamiltonians are retained. Consequently, the here presented formulation should be considered as the new standard for DFT/MRCI calculations.
Collapse
Affiliation(s)
- Dennis R Dombrowski
- Institute of Theoretical and Computational Chemistry, Heinrich-Heine-University Düsseldorf, 40225 Düsseldorf, Germany
| | - Timo Schulz
- Institute of Theoretical and Computational Chemistry, Heinrich-Heine-University Düsseldorf, 40225 Düsseldorf, Germany
| | - Martin Kleinschmidt
- Institute of Theoretical and Computational Chemistry, Heinrich-Heine-University Düsseldorf, 40225 Düsseldorf, Germany
| | - Christel M Marian
- Institute of Theoretical and Computational Chemistry, Heinrich-Heine-University Düsseldorf, 40225 Düsseldorf, Germany
| |
Collapse
|
8
|
Manawadu D, Georges TN, Barford W. Photoexcited State Dynamics and Singlet Fission in Carotenoids. J Phys Chem A 2023; 127:1342-1352. [PMID: 36701532 PMCID: PMC9923744 DOI: 10.1021/acs.jpca.2c07781] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/27/2023]
Abstract
We describe our simulations of the excited state dynamics of the carotenoid neurosporene, following its photoexcitation into the "bright" (nominally 11Bu+) state. To account for the experimental and theoretical uncertainty in the relative energetic ordering of the nominal 11Bu+ and 21Ag- states at the Franck-Condon point, we consider two parameter sets. In both cases, there is ultrafast internal conversion from the "bright" state to a "dark" singlet triplet-pair state, i.e., to one member of the "2Ag" family of states. For one parameter set, internal conversion from the 11Bu+ to 21Ag- states occurs via the dark, intermediate 11Bu- state. In this case, there is a cross over of the 11Bu+ and 11Bu- diabatic energies within 5 fs and an associated avoided crossing of the S2 and S3 adiabatic energies. After the adiabatic evolution of the S2 state from predominately 11Bu+ character to predominately 11Bu- character, there is a slower nonadiabatic transition from S2 to S1, accompanied by an increase in the population of the 21Ag- state. For the other parameter set, the 21Ag- energy lies higher than the 11Bu+ energy at the Franck-Condon point. In this case, there is cross over of the 21Ag- and 11Bu+ energies and an avoided crossing of the S1 and S2 energies, as the S1 state evolves adiabatically from being of 11Bu+ character to 21Ag- character. We make a direct connection from our predictions to experimental observables by calculating the time-resolved excited state absorption. For the case of direct 11Bu+ to 21Ag- internal conversion, we show that the dominant transition at ca. 2 eV, being close to but lower in energy than the T1 to T1* transition, can be attributed to the 21Ag- component of S1. Moreover, we show that it is the charge-transfer exciton component of the 21Ag- state that is responsible for this transition (to a higher-lying exciton state), and not its triplet-pair component. These simulations are performed using the adaptive tDMRG method on the extended Hubbard model of π-conjugated electrons. The Ehrenfest equations of motion are used to simulate the coupled nuclei dynamics. We next discuss the microscopic mechanism of "bright" to "dark" state internal conversion and emphasize that this occurs via the exciton components of both states. Finally, we describe a mechanism relying on torsional relaxation whereby the strongly bound intrachain triplet-pairs of the "dark" state may undergo interchain exothermic dissociation.
Collapse
Affiliation(s)
- Dilhan Manawadu
- Department
of Chemistry, Physical and Theoretical Chemistry Laboratory, University of Oxford, Oxford OX1 3QZ, United Kingdom,Linacre
College, University of Oxford, Oxford OX1 3JA, United Kingdom,E-mail:
| | - Timothy N. Georges
- Department
of Chemistry, Physical and Theoretical Chemistry Laboratory, University of Oxford, Oxford OX1 3QZ, United Kingdom,Brasenose
College, University of Oxford, Oxford OX1 4AJ, United Kingdom
| | - William Barford
- Department
of Chemistry, Physical and Theoretical Chemistry Laboratory, University of Oxford, Oxford OX1 3QZ, United Kingdom,E-mail:
| |
Collapse
|
9
|
Nikoobakht B, Hakim R, Menger MF, Köppel H. On the UV spectroscopy and photodynamics of octatetraene. Mol Phys 2022. [DOI: 10.1080/00268976.2022.2132186] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/17/2022]
Affiliation(s)
- Behnam Nikoobakht
- Theoretische Chemie, Physikalisch-Chemisches Institut, Universität Heidelberg, Heidelberg, Germany
| | - Raymond Hakim
- Theoretische Chemie, Physikalisch-Chemisches Institut, Universität Heidelberg, Heidelberg, Germany
| | - Maximilian F.S.J. Menger
- Theoretische Chemie, Physikalisch-Chemisches Institut, Universität Heidelberg, Heidelberg, Germany
| | - Horst Köppel
- Theoretische Chemie, Physikalisch-Chemisches Institut, Universität Heidelberg, Heidelberg, Germany
| |
Collapse
|
10
|
Khokhlov D, Belov A. Low-Lying Excited States of Natural Carotenoids Viewed by Ab Initio Methods. J Phys Chem A 2022; 126:4376-4391. [PMID: 35767689 DOI: 10.1021/acs.jpca.2c02485] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Low-lying excited states of carotenoids (the optically dark 2Ag- and bright 1Bu+) are deeply involved in energy transfer processes in photosynthetic antennas, such as light harvesting and non-photochemical quenching. Though any ab initio modeling of these phenomena has to rely on precise energies of the carotenoid electronic states, the accurate evaluation of these states remains a challenging problem due to their different natures. The paper aims to study the accuracy of the excitation energies of the low-lying excited states of certain open- and closed-chain carotenoids obtained by a state-of-the-art multireference approach for electronic structure calculation. Here, density matrix renormalization group SCF (DMRGSCF) and a perturbative approach based on driven similarity renormalization group second-order multireference perturbation theory (DSRG-MRPT2) were used to treat the static and dynamic correlation, respectively. Nuclear geometries of the electronic states were optimized with DFT-based approaches. It is demonstrated that spin-flip TD-DFT can replace multiconfigurational methods for the geometry optimization of the 2Ag- state but not for the calculation of the excitation energy. Adiabatic excitation energies to the 1Bu+ state were shown to be within a margin of 1000 cm-1 with an appropriate flow parameter value. Adiabatic excitation energies to the 2Ag- state for the open-chain carotenoids lie within a range of experimental values (taking into account the broad range of experimental estimates); for the closed-chain ones, the error does not exceed 2000 cm-1. Ab initio stationary (1Ag- → 1Bu+) and transient (2Ag- → 1Bu+) absorption spectra were modeled for violaxanthin and lycopene, and these spectra showed good agreement with the experimental ones both in terms of the vibronic structure and the transition energies.
Collapse
Affiliation(s)
- Daniil Khokhlov
- Department of Chemistry, Lomonosov Moscow State University, Moscow 119991, Russia
| | - Aleksandr Belov
- Department of Chemistry, Lomonosov Moscow State University, Moscow 119991, Russia
| |
Collapse
|
11
|
Manawadu D, Valentine DJ, Marcus M, Barford W. Singlet Triplet-Pair Production and Possible Singlet-Fission in Carotenoids. J Phys Chem Lett 2022; 13:1344-1349. [PMID: 35108016 PMCID: PMC9084603 DOI: 10.1021/acs.jpclett.1c03812] [Citation(s) in RCA: 9] [Impact Index Per Article: 4.5] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/19/2021] [Accepted: 02/01/2022] [Indexed: 05/19/2023]
Abstract
Internal conversion from the photoexcited state to a correlated singlet triplet-pair state is believed to be the precursor of singlet fission in carotenoids. We present numerical simulations of this process using a π-electron model that fully accounts for electron-electron interactions and electron-nuclear coupling. The time-evolution of the electrons is determined rigorously using the time-dependent density matrix renormalization group method, while the nuclei are evolved via the Ehrenfest equations of motion. We apply this to zeaxanthin, a carotenoid chain with 18 fully conjugated carbon atoms. We show that the internal conversion of the primary photoexcited state, S2, to the singlet triplet-pair state occurs adiabatically via an avoided crossing within ∼50 fs with a yield of ∼60%. We further discuss whether this singlet triplet-pair state will undergo exothermic versus endothermic intra- or interchain singlet fission.
Collapse
Affiliation(s)
- Dilhan Manawadu
- Department
of Chemistry, Physical and Theoretical Chemistry Laboratory, University of Oxford, Oxford OX1 3QZ, United Kingdom
- Linacre College, University of Oxford, Oxford OX1 3JA, United Kingdom
| | - Darren J. Valentine
- Department
of Chemistry, Physical and Theoretical Chemistry Laboratory, University of Oxford, Oxford OX1 3QZ, United Kingdom
- Balliol College, University of Oxford, Oxford OX1 3BJ, United Kingdom
| | - Max Marcus
- Department
of Chemistry, Physical and Theoretical Chemistry Laboratory, University of Oxford, Oxford OX1 3QZ, United Kingdom
| | - William Barford
- Department
of Chemistry, Physical and Theoretical Chemistry Laboratory, University of Oxford, Oxford OX1 3QZ, United Kingdom
| |
Collapse
|
12
|
Cheng Y, Xie Z, Ma H. Post-Density Matrix Renormalization Group Methods for Describing Dynamic Electron Correlation with Large Active Spaces. J Phys Chem Lett 2022; 13:904-915. [PMID: 35049302 DOI: 10.1021/acs.jpclett.1c04078] [Citation(s) in RCA: 15] [Impact Index Per Article: 7.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/14/2023]
Abstract
The ab initio density matrix renormalization group (DMRG) method has been well-established and has become one of the most accurate numerical methods for the precise electronic structure solution of large active spaces. In the past few years, to capture the missing dynamic correlation, various post-DMRG approaches have been proposed through the combination of DMRG and multireference quantum chemical methods or density functional theory. With this in mind, this work provides a brief overview of ab initio DMRG principles and the new developments within post-DMRG methods. For clarity, post-DMRG methods are classified into two main categories depending on whether high-order n-electron reduced density matrices are used, and their merits and disadvantages are properly discussed. Finally, we conclude by discussing unsolved bottlenecks and giving development perspectives of post-DMRG approaches, which are expected to yield quantitative descriptions of complex electronic structures in large strongly correlated molecules and materials.
Collapse
Affiliation(s)
- Yifan Cheng
- School of Chemistry and Chemical Engineering, Jiangsu Key Laboratory of Vehicle Emissions Control, Nanjing University, Nanjing 210023, China
| | - Zhaoxuan Xie
- School of Chemistry and Chemical Engineering, Jiangsu Key Laboratory of Vehicle Emissions Control, Nanjing University, Nanjing 210023, China
| | - Haibo Ma
- School of Chemistry and Chemical Engineering, Jiangsu Key Laboratory of Vehicle Emissions Control, Nanjing University, Nanjing 210023, China
- Chemistry and Biomedicine Innovation Center (ChemBIC), Nanjing University, Nanjing 210023, China
| |
Collapse
|
13
|
Bondanza M, Jacquemin D, Mennucci B. Excited States of Xanthophylls Revisited: Toward the Simulation of Biologically Relevant Systems. J Phys Chem Lett 2021; 12:6604-6612. [PMID: 34251826 PMCID: PMC8311646 DOI: 10.1021/acs.jpclett.1c01929] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/16/2021] [Accepted: 07/08/2021] [Indexed: 06/13/2023]
Abstract
Xanthophylls are a class of oxygen-containing carotenoids, which play a fundamental role in light-harvesting pigment-protein complexes and in many photoresponsive proteins. The complexity of the manifold of the electronic states and the large sensitivity to the environment still prevent a clear and coherent interpretation of their photophysics and photochemistry. In this Letter, we compare cutting-edge ab initio methods (CC3 and DMRG/NEVPT2) with time-dependent DFT and semiempirical CI (SECI) on model keto-carotenoids and show that SECI represents the right compromise between accuracy and computational cost to be applied to real xanthophylls in their biological environment. As an example, we investigate canthaxanthin in the orange carotenoid protein and show that the conical intersections between excited states and excited-ground states are mostly determined by the effective bond length alternation coordinate, which is significantly tuned by the protein through geometrical constraints and electrostatic effects.
Collapse
Affiliation(s)
- Mattia Bondanza
- Dipartimento
di Chimica e Chimica Industriale, University
of Pisa, via G. Moruzzi 13, 56124 Pisa, Italy
| | - Denis Jacquemin
- Université
de Nantes, CNRS, CEISAM UMR 6230, F-44000 Nantes, France
| | - Benedetta Mennucci
- Dipartimento
di Chimica e Chimica Industriale, University
of Pisa, via G. Moruzzi 13, 56124 Pisa, Italy
| |
Collapse
|
14
|
Khokhlov D, Belov A. Toward an Accurate Ab Initio Description of Low-Lying Singlet Excited States of Polyenes. J Chem Theory Comput 2021; 17:4301-4315. [PMID: 34125516 DOI: 10.1021/acs.jctc.0c01293] [Citation(s) in RCA: 12] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/24/2023]
Abstract
The low-lying excited states of carotenoids play a crucial role in many important biophysical processes such as photosynthesis. Most of these excited states are strongly correlated, which makes them both challenging for a qualitative ab initio description and an engaging model system for trying out emerging multireference methods. Among these methods, driven similarity renormalization group (DSRG) and its perturbative version (DSRG-MRPT2) are especially attractive in terms of both accuracy and moderate numerical complexity. In this paper, we applied density matrix renormalization group (DMRG) followed by DSRG-MRPT2 for the calculation of vertical and adiabatic excitation energies into the 2Ag-, 1Bu-, and 1Bu+ electronic states of polyenes containing from 8 to 13 conjugating double bonds acting as a model for natural carotenoids. It was shown that the DSRG flow parameter should be adjusted to ensure both the energy convergence with respect to it and the agreement with the experimental data. With the increased flow parameter, the proposed combination of methods provides a reasonable agreement with the experiment. The deviations of the adiabatic excitation energies are less than 1000 cm-1 for the 2Ag- and less than 3000 cm-1 for the excited states of the Bu symmetry, which in terms of accuracy significantly outperforms the N-electron valence state perturbation theory. At the same time, DSRG-MRPT2 is shown to be robust with respect to variation of quality of the DMRG reference wave function such as the orbital optimization or the number of electronic states in the averaging.
Collapse
Affiliation(s)
- Daniil Khokhlov
- Department of Chemistry, Lomonosov Moscow State University, Moscow, 119991, Russia
| | - Aleksandr Belov
- Department of Chemistry, Lomonosov Moscow State University, Moscow, 119991, Russia
| |
Collapse
|