1
|
Wang J, Durbeej B. Thermal Fluctuations in Conjugation and their Effect on Calculated Excitation Energies: A Case Study on the Astaxanthin Carotenoid. CHEMPHOTOCHEM 2021. [DOI: 10.1002/cptc.202100178] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/09/2022]
Affiliation(s)
- Jun Wang
- Jiangsu Key Laboratory for Chemistry of Low-Dimensional Materials Jiangsu Engineering Laboratory for Environment Functional Materials School of Chemistry and Chemical Engineering Huaiyin Normal University No. 111 West Changjiang Road 223300 Huaian Jiangsu Province China
| | - Bo Durbeej
- Division of Theoretical Chemistry IFM Linköping University 581 83 Linköping Sweden
| |
Collapse
|
2
|
Park W, Shen J, Lee S, Piecuch P, Filatov M, Choi CH. Internal Conversion between Bright (1 1Bu+) and Dark (2 1Ag-) States in s- trans-Butadiene and s- trans-Hexatriene. J Phys Chem Lett 2021; 12:9720-9729. [PMID: 34590847 DOI: 10.1021/acs.jpclett.1c02707] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/13/2023]
Abstract
Internal conversion (IC) between the two lowest singlet excited states, 11Bu+ and 21Ag-, of s-trans-butadiene and s-trans-hexatriene is investigated using a series of single- and multi- reference wave function and density functional theory (DFT) methodologies. Three independent types of the equation-of-motion coupled-cluster (EOMCC) theory capable of providing an accurate and balanced description of one- as well as two-electron transitions, abbreviated as δ-CR-EOMCC(2,3), DIP-EOMCC(4h2p){No}, and DEA-EOMCC(4p2h){Nu} or DEA-EOMCC(3p1h,4p2h){Nu}, consistently predict that the 11Bu+/21Ag- crossing in both molecules occurs along the bond length alternation coordinate. However, the analogous 11Bu+ and 21Ag- potentials obtained with some multireference approaches, such as CASSCF and MRCIS(D), as well as with the linear-response formulation of time-dependent DFT (TDDFT), do not cross. Hence, caution needs to be exercised when studying the low-lying singlet excited states of polyenes with conventional multiconfigurational methods and TDDFT. The multistate many-body perturbation theory methods, such as XMCQDPT2, do correctly reproduce the curve crossing. Among the simplest and least expensive computational methodologies, the DFT approaches that incorporate the contributions of doubly excited configurations, abbreviated as MRSF (mixed reference spin-flip) TDDFT and SSR(4,4), accurately reproduce our best EOMCC results. This is highly promising for nonadiabatic molecular dynamics simulations in larger systems.
Collapse
Affiliation(s)
- Woojin Park
- Department of Chemistry, Kyungpook National University, Daegu 41566, South Korea
| | - Jun Shen
- Department of Chemistry, Michigan State University, East Lansing, Michigan 48824, United States
| | - Seunghoon Lee
- Division of Chemistry and Chemical Engineering, California Institute of Technology, Pasadena, California 91125, United States
| | - Piotr Piecuch
- Department of Chemistry, Michigan State University, East Lansing, Michigan 48824, United States
- Department of Physics and Astronomy, Michigan State University, East Lansing, Michigan 48824, United States
| | - Michael Filatov
- Department of Chemistry, Kyungpook National University, Daegu 41566, South Korea
| | - Cheol Ho Choi
- Department of Chemistry, Kyungpook National University, Daegu 41566, South Korea
| |
Collapse
|
3
|
Hutcheson A, Paul AC, Myhre RH, Koch H, Høyvik IM. Describing ground and excited state potential energy surfaces for molecular photoswitches using coupled cluster models. J Comput Chem 2021; 42:1419-1429. [PMID: 33973669 DOI: 10.1002/jcc.26553] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/16/2020] [Revised: 03/19/2021] [Accepted: 04/20/2021] [Indexed: 11/08/2022]
Abstract
In this article, we use two extensively studied systems, a retinal model system and azobenzene, to explore the use of coupled cluster models for describing ground and singlet excited state potential energy surfaces of photoswitchable systems. While not being suitable for describing nuclear dynamics of photoisomerization, coupled cluster models have useful attributes, such as the inclusion of dynamical correlation, their black box nature, and the systematic improvement offered by truncation level. Results for the studied systems show that when triple excitations (here through the CC3 model) are included, ground and excited state potential energy surfaces for isomerization paths may reliably be generated, also for states of doubly excited character. For ground state equilibrium cis- and trans-azobenzene, the molecular geometry and basis set is seen to significantly impact the vertical excitation energies for the two lowest excited states. Efficient implementations of coupled cluster models can therefore constitute valuable tools for investigating photoswitchable systems and can be used for preliminary black box studies to gather information before more complicated excited state dynamics approaches are pursued.
Collapse
Affiliation(s)
- Anders Hutcheson
- Department of Chemistry, The Norwegian University of Science and Technology, Trondheim, Norway
| | | | - Rolf H Myhre
- Department of Chemistry, The Norwegian University of Science and Technology, Trondheim, Norway
| | | | - Ida-Marie Høyvik
- Department of Chemistry, The Norwegian University of Science and Technology, Trondheim, Norway
| |
Collapse
|
4
|
Park JW, Al-Saadon R, MacLeod MK, Shiozaki T, Vlaisavljevich B. Multireference Electron Correlation Methods: Journeys along Potential Energy Surfaces. Chem Rev 2020; 120:5878-5909. [PMID: 32239929 DOI: 10.1021/acs.chemrev.9b00496] [Citation(s) in RCA: 54] [Impact Index Per Article: 13.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/16/2022]
Abstract
Multireference electron correlation methods describe static and dynamical electron correlation in a balanced way and, therefore, can yield accurate and predictive results even when single-reference methods or multiconfigurational self-consistent field theory fails. One of their most prominent applications in quantum chemistry is the exploration of potential energy surfaces. This includes the optimization of molecular geometries, such as equilibrium geometries and conical intersections and on-the-fly photodynamics simulations, both of which depend heavily on the ability of the method to properly explore the potential energy surface. Because such applications require nuclear gradients and derivative couplings, the availability of analytical nuclear gradients greatly enhances the scope of quantum chemical methods. This review focuses on the developments and advances made in the past two decades. A detailed account of the analytical nuclear gradient and derivative coupling theories is presented. Emphasis is given to the software infrastructure that allows one to make use of these methods. Notable applications of multireference electron correlation methods to chemistry, including geometry optimizations and on-the-fly dynamics, are summarized at the end followed by a discussion of future prospects.
Collapse
Affiliation(s)
- Jae Woo Park
- Department of Chemistry, Chungbuk National University, Chungdae-ro 1, Cheongju 28644, Korea
| | - Rachael Al-Saadon
- Department of Chemistry, Northwestern University, 2145 Sheridan Road, Evanston, Illinois 60208, United States
| | - Matthew K MacLeod
- Workday, 4900 Pearl Circle East, Suite 100, Boulder, Colorado 80301, United States
| | - Toru Shiozaki
- Department of Chemistry, Northwestern University, 2145 Sheridan Road, Evanston, Illinois 60208, United States.,Quantum Simulation Technologies, Inc., 625 Massachusetts Avenue, Cambridge, Massachusetts 02139, United States
| | - Bess Vlaisavljevich
- Department of Chemistry, University of South Dakota, 414 East Clark Street, Vermillion, South Dakota 57069, United States
| |
Collapse
|
5
|
Theoretical studies of the ground and excited states of the NaLi molecule by the IOTC CASPT2 method. Struct Chem 2019. [DOI: 10.1007/s11224-019-01332-x] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/26/2022]
|
6
|
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
|
7
|
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
|
8
|
Lischka H, Nachtigallová D, Aquino AJA, Szalay PG, Plasser F, Machado FBC, Barbatti M. Multireference Approaches for Excited States of Molecules. Chem Rev 2018; 118:7293-7361. [DOI: 10.1021/acs.chemrev.8b00244] [Citation(s) in RCA: 197] [Impact Index Per Article: 32.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/18/2022]
Affiliation(s)
- Hans Lischka
- School of Pharmaceutical Sciences and Technology, Tianjin University, Tianjin 300072, P.R. China
- Department of Chemistry and Biochemistry, Texas Tech University, Lubbock, Texas 79409, United States
- Institute of Theoretical Chemistry, Faculty of Chemistry, University of Vienna, Währinger Straße 17, 1090 Vienna, Austria
| | - Dana Nachtigallová
- Institute of Organic Chemistry and Biochemistry v.v.i., The Czech Academy of Sciences, Flemingovo nám. 2, 16610 Prague 6, Czech Republic
- Regional Centre of Advanced Technologies and Materials, Palacký University, 78371 Olomouc, Czech Republic
| | - Adélia J. A. Aquino
- School of Pharmaceutical Sciences and Technology, Tianjin University, Tianjin 300072, P.R. China
- Department of Chemistry and Biochemistry, Texas Tech University, Lubbock, Texas 79409, United States
- Institute for Soil Research, University of Natural Resources and Life Sciences Vienna, Peter-Jordan-Strasse 82, A-1190 Vienna, Austria
| | - Péter G. Szalay
- ELTE Eötvös Loránd University, Laboratory of Theoretical Chemistry, Pázmány Péter sétány 1/A, 1117 Budapest, Hungary
| | - Felix Plasser
- Institute of Theoretical Chemistry, Faculty of Chemistry, University of Vienna, Währinger Straße 17, 1090 Vienna, Austria
- Department of Chemistry, Loughborough University, Leicestershire LE11 3TU, United Kingdom
| | - Francisco B. C. Machado
- Departamento de Química, Instituto Tecnológico de Aeronáutica, São José dos Campos 12228-900, São Paulo, Brazil
| | | |
Collapse
|
9
|
Cerezo J, Prampolini G, Cacelli I. Developing accurate intramolecular force fields for conjugated systems through explicit coupling terms. Theor Chem Acc 2018. [DOI: 10.1007/s00214-018-2254-8] [Citation(s) in RCA: 23] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/28/2022]
|
10
|
Gozem S, Luk HL, Schapiro I, Olivucci M. Theory and Simulation of the Ultrafast Double-Bond Isomerization of Biological Chromophores. Chem Rev 2017; 117:13502-13565. [DOI: 10.1021/acs.chemrev.7b00177] [Citation(s) in RCA: 175] [Impact Index Per Article: 25.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/03/2023]
Affiliation(s)
- Samer Gozem
- Department
of Chemistry, Georgia State University, Atlanta, Georgia 30302, United States
| | - Hoi Ling Luk
- Chemistry
Department, Bowling Green State University, Overman Hall, Bowling Green, Ohio 43403, United States
| | - Igor Schapiro
- Fritz
Haber Center for Molecular Dynamics, Institute of Chemistry, The Hebrew University of Jerusalem, Jerusalem 91904, Israel
| | - Massimo Olivucci
- Chemistry
Department, Bowling Green State University, Overman Hall, Bowling Green, Ohio 43403, United States
- Dipartimento
di Biotecnologie, Chimica e Farmacia, Università di Siena, via A. Moro
2, 53100 Siena, Italy
| |
Collapse
|
11
|
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
|
12
|
Momeni MR, Brown A. Why do TD-DFT excitation energies of BODIPY/Aza-BODIPY families largely deviate from experiment? Answers from electron correlated and multireference methods. J Chem Theory Comput 2016; 11:2619-32. [PMID: 26575559 DOI: 10.1021/ct500775r] [Citation(s) in RCA: 126] [Impact Index Per Article: 15.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/12/2022]
Abstract
The vertical excitation energies of 17 boron-dipyrromethene (BODIPY) core structures with a variety of substituents and ring sizes are benchmarked using time-dependent density functional theory (TD-DFT) with nine different functionals combined with the cc-pVTZ basis set. When compared to experimental measurements, all functionals provide mean absolute errors (mean AEs) greater than 0.3 eV, larger than the 0.1-0.3 eV differences typically expected from TD-DFT. Due to the high linear correlation of TD-DFT results with experiment, most functionals can be used to predict excitation energies if corrected empirically. Using the CAM-B3LYP functional, 0-0 transition energies are determined, and while the absolute difference is improved (mean AE = 0.478 eV compared to 0.579 eV), the correlation diminishes substantially (R(2) = 0.961 to 0.862). Two very recently introduced charge transfer (CT) indices, q(CT) and d(CT), and electron density difference (EDD) plots demonstrate that CT does not play a significant role for most of the BODIPYs examined and, thus, cannot be the source of error in TD-DFT. To assess TD-DFT methods, vertical excitation energies are determined utilizing TD-HF, configuration interaction CIS and CIS(D), equation of motion EOM-CCSD, SAC-CI, and Laplace-transform based local coupled-cluster singles and approximate doubles LCC2* methods. Moreover, multireference CASSCF and CASPT2 vertical excitation energies were also obtained for all species (except CASPT2 was not feasible for the four largest systems). The SAC-CI/cc-pVDZ, LCC2*/cc-pVDZ, and CASPT2/cc-pVDZ approaches are shown to have the smallest mean AEs of 0.154, 0.109, and 0.100 eV, respectively; the utility of the LCC2* approach is demonstrated for eight extended BODIPYs and aza-BODIPYs. We found that the problems with TD-DFT arise from difficulties in dealing with the differential electron correlation (as assessed by comparing CCS, CC2, LR-CCSD, CCSDR(T), and CCSDR(3) vertical excitation energies for five compounds) and from contributions of multireference character and double excitations (from analysis of the CASSCF wave functions).
Collapse
Affiliation(s)
- Mohammad R Momeni
- Department of Chemistry, University of Alberta , Edmonton, Alberta T6G 2G2, Canada
| | - Alex Brown
- Department of Chemistry, University of Alberta , Edmonton, Alberta T6G 2G2, Canada
| |
Collapse
|
13
|
Kupka T, Buczek A, Broda MA, Stachów M, Tarnowski P. DFT studies on the structural and vibrational properties of polyenes. J Mol Model 2016; 22:101. [PMID: 27048200 PMCID: PMC4821862 DOI: 10.1007/s00894-016-2969-1] [Citation(s) in RCA: 20] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/22/2015] [Accepted: 03/16/2016] [Indexed: 11/29/2022]
Abstract
Detailed density functional theory (DFT) calculations on the structure and harmonic frequencies of model all-trans and all-cis polyenes were undertaken. For the first time, we report on the convergence of selected B3LYP/6-311++G** and BLYP/6-311++G** calculated structural parameters resulting from a systematic increase in polyene size (chains containing 2 to 14 C = C units). The limiting values of the structural parameters for very long chains were estimated using simple three-parameter empirical formulae. BLYP/6-311++G** calculated ν(C = C) and ν(C-C) frequencies for all-trans and all-cis polyenes containing up to 14 carbon-carbon double bonds were used to estimate these values for very long chains. Correction of raw, unscaled vibrational data was performed by comparing theoretical and experimental wavenumbers for polyenes chains containing 3 to 12 conjugated C = C units with both ends substituted by tert-butyl groups. The corrected ν(C = C) and ν(C-C) wavenumbers for all-trans molecules were used to estimate the presence of 9 - 12 C = C units in all-trans polyene pigment in red coral.
Collapse
Affiliation(s)
- Teobald Kupka
- Faculty of Chemistry, University of Opole, 48 Oleska Street, 45-052, Opole, Poland.
| | - Aneta Buczek
- Faculty of Chemistry, University of Opole, 48 Oleska Street, 45-052, Opole, Poland.
| | - Małgorzata A Broda
- Faculty of Chemistry, University of Opole, 48 Oleska Street, 45-052, Opole, Poland
| | - Michał Stachów
- Faculty of Chemistry, University of Opole, 48 Oleska Street, 45-052, Opole, Poland
| | - Przemysław Tarnowski
- Faculty of Chemistry, University of Opole, 48 Oleska Street, 45-052, Opole, Poland
| |
Collapse
|
14
|
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
|
15
|
Casanova D, Krylov AI. Quantifying local exciton, charge resonance, and multiexciton character in correlated wave functions of multichromophoric systems. J Chem Phys 2016; 144:014102. [DOI: 10.1063/1.4939222] [Citation(s) in RCA: 27] [Impact Index Per Article: 3.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/12/2023] Open
Affiliation(s)
- David Casanova
- Kimika Fakultatea, Euskal Herriko Unibersitatea (UPV/EHU) and Donostia International Physics Center (DIPC), P.K. 1072, 20018 Donostia, Spain
- IKERBASQUE, Basque Foundation for Science, 48013 Bilbao, Spain
| | - Anna I. Krylov
- Department of Chemistry, University of Southern California, Los Angeles, California 90089-0482, USA
| |
Collapse
|
16
|
Tuna D, Lefrancois D, Wolański Ł, Gozem S, Schapiro I, Andruniów T, Dreuw A, Olivucci M. Assessment of Approximate Coupled-Cluster and Algebraic-Diagrammatic-Construction Methods for Ground- and Excited-State Reaction Paths and the Conical-Intersection Seam of a Retinal-Chromophore Model. J Chem Theory Comput 2015; 11:5758-81. [PMID: 26642989 DOI: 10.1021/acs.jctc.5b00022] [Citation(s) in RCA: 107] [Impact Index Per Article: 11.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/16/2022]
Abstract
As a minimal model of the chromophore of rhodopsin proteins, the penta-2,4-dieniminium cation (PSB3) poses a challenging test system for the assessment of electronic-structure methods for the exploration of ground- and excited-state potential-energy surfaces, the topography of conical intersections, and the dimensionality (topology) of the branching space. Herein, we report on the performance of the approximate linear-response coupled-cluster method of second order (CC2) and the algebraic-diagrammatic-construction scheme of the polarization propagator of second and third orders (ADC(2) and ADC(3)). For the ADC(2) method, we considered both the strict and extended variants (ADC(2)-s and ADC(2)-x). For both CC2 and ADC methods, we also tested the spin-component-scaled (SCS) and spin-opposite-scaled (SOS) variants. We have explored several ground- and excited-state reaction paths, a circular path centered around the S1/S0 surface crossing, and a 2D scan of the potential-energy surfaces along the branching space. We find that the CC2 and ADC methods yield a different dimensionality of the intersection space. While the ADC methods yield a linear intersection topology, we find a conical intersection topology for the CC2 method. We present computational evidence showing that the linear-response CC2 method yields a surface crossing between the reference state and the first response state featuring characteristics that are expected for a true conical intersection. Finally, we test the performance of these methods for the approximate geometry optimization of the S1/S0 minimum-energy conical intersection and compare the geometries with available data from multireference methods. The present study provides new insight into the performance of linear-response CC2 and polarization-propagator ADC methods for molecular electronic spectroscopy and applications in computational photochemistry.
Collapse
Affiliation(s)
- Deniz Tuna
- Max-Planck-Institut für Kohlenforschung , 45470 Mülheim an der Ruhr, Germany
| | - Daniel Lefrancois
- Interdisciplinary Center for Scientific Computing, University of Heidelberg , 69120 Heidelberg, Germany
| | - Łukasz Wolański
- Department of Chemistry, Wrocław University of Technology , 50370 Wrocław, Poland
| | - Samer Gozem
- Department of Chemistry, University of Southern California , Los Angeles, California 90089, United States
| | - Igor Schapiro
- Institut de Physique et Chimie des Matériaux de Strasbourg & Labex NIE, Université de Strasbourg, CNRS UMR 7504 , Strasbourg 67034, France
| | - Tadeusz Andruniów
- Department of Chemistry, Wrocław University of Technology , 50370 Wrocław, Poland
| | - Andreas Dreuw
- Interdisciplinary Center for Scientific Computing, University of Heidelberg , 69120 Heidelberg, Germany
| | - Massimo Olivucci
- Department of Chemistry, Bowling Green State University , Bowling Green, Ohio 43402, United States.,Dipartimento di Biotecnologie, Chimica e Farmacia, Universitá de Siena , 53100 Siena, Italy
| |
Collapse
|
17
|
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
|
18
|
Walczak E, Andruniów T. Impacts of retinal polyene (de)methylation on the photoisomerization mechanism and photon energy storage of rhodopsin. Phys Chem Chem Phys 2015; 17:17169-81. [DOI: 10.1039/c5cp01939g] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Similar to native rhodopsin, a two-mode space-saving isomerization mechanism drives the photoreaction in (de)methylated rhodopsin analogues.
Collapse
Affiliation(s)
- Elżbieta Walczak
- Department of Chemistry
- Wroclaw University of Technology
- 50-370 Wroclaw
- Poland
| | - Tadeusz Andruniów
- Department of Chemistry
- Wroclaw University of Technology
- 50-370 Wroclaw
- Poland
| |
Collapse
|
19
|
Champagne B, Liégeois V, Zutterman F. Pigment violet 19 - a test case to define a simple method to simulate the vibronic structure of absorption spectra of organic pigments and dyes in solution. Photochem Photobiol Sci 2014; 14:444-56. [PMID: 25501947 DOI: 10.1039/c4pp00317a] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
A typical quinacridone pigment, PV19, has been used to analyze the impact of several computational parameters on the UV/vis absorption band shape in solution, simulated using density functional theory and time-dependent density functional theory levels of approximation. These encompass, (i) the choice of exchange-correlation functional, (ii) the basis set, (iii) the method for non-equilibrium optimization of the excited state geometry, (iv) the approach for evaluating the vibronic band structure, (v) the peak broadening, and (vi) the scaling of the harmonic vibrational frequencies. Among these, the choice of exchange-correlation functional is certainly of the most importance because it can drastically modify the spectral shape. In the case of PV19, the M05-2X and to a lesser extent CAM-B3LYP XC functionals are the most efficient to reproduce the vibronic structure, confirming the important role of exact Hartree-Fock exchange. Still, these functionals are not the most reliable to predict the excitation energies and oscillator strengths, for which M05, a functional with less HF exchange, performs better. For evaluating the vibronic structure, the simple gradient method, where only one step of geometry optimization of the excited state is carried out and the gradients are used to evaluate the Huang-Rhys factors as well as to determine the excited state geometries produces a spectrum that is very similar to the ones obtained with the more involved Duschinsky and geometry methods, opening the way to a fast simulation of the UV/vis absorption spectra of pigments and dyes. Then, the effect of scaling the calculated vibrational frequencies to account for anharmonicity effects as well as for limitation of the method also impacts the shape of the vibronic spectrum and this effect depends on the method used to determine the Huang-Rhys factors. Indeed, scaling the vibrational frequencies by a factor which is typically smaller than 1.0 results in a relative decrease of the 0-1 peak intensity with respect to the 0-0 band when optimizing the geometry of the excited state whereas the effect is opposite and magnified if using the gradient method.
Collapse
Affiliation(s)
- Benoît Champagne
- Laboratoire de Chimie Théorique, Unité de Chimie-Physique Théorique et Structurale, Université de Namur, rue de Bruxelles, 61, B-5000 Namur, Belgium.
| | | | | |
Collapse
|
20
|
Tan YS, Urbančok D, Webster RD. Contrasting Voltammetric Behavior of Different Forms of Vitamin A in Aprotic Organic Solvents. J Phys Chem B 2014; 118:8591-600. [DOI: 10.1021/jp505456q] [Citation(s) in RCA: 12] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Ying Shan Tan
- Division of Chemistry and
Biological Chemistry, School of Physical and Mathematical Sciences, Nanyang Technological University, Singapore 637371
| | - Dejan Urbančok
- Division of Chemistry and
Biological Chemistry, School of Physical and Mathematical Sciences, Nanyang Technological University, Singapore 637371
| | - Richard D. Webster
- Division of Chemistry and
Biological Chemistry, School of Physical and Mathematical Sciences, Nanyang Technological University, Singapore 637371
| |
Collapse
|
21
|
Campomanes P, Neri M, Horta BAC, Röhrig UF, Vanni S, Tavernelli I, Rothlisberger U. Origin of the Spectral Shifts among the Early Intermediates of the Rhodopsin Photocycle. J Am Chem Soc 2014; 136:3842-51. [DOI: 10.1021/ja411303v] [Citation(s) in RCA: 36] [Impact Index Per Article: 3.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/04/2023]
Affiliation(s)
- Pablo Campomanes
- Laboratory
of Computational Chemistry and Biochemistry, Ecole Polytechnique Fédérale Lausanne, CH-1015 Lausanne, Switzerland
| | - Marilisa Neri
- Laboratory
of Computational Chemistry and Biochemistry, Ecole Polytechnique Fédérale Lausanne, CH-1015 Lausanne, Switzerland
| | - Bruno A. C. Horta
- Laboratory
of Computational Chemistry and Biochemistry, Ecole Polytechnique Fédérale Lausanne, CH-1015 Lausanne, Switzerland
| | - Ute F. Röhrig
- Molecular Modeling
Group, Swiss Institute of
Bioinformatics, CH-1015 Lausanne, Switzerland
| | - Stefano Vanni
- Laboratory
of Computational Chemistry and Biochemistry, Ecole Polytechnique Fédérale Lausanne, CH-1015 Lausanne, Switzerland
| | - Ivano Tavernelli
- Laboratory
of Computational Chemistry and Biochemistry, Ecole Polytechnique Fédérale Lausanne, CH-1015 Lausanne, Switzerland
| | - Ursula Rothlisberger
- Laboratory
of Computational Chemistry and Biochemistry, Ecole Polytechnique Fédérale Lausanne, CH-1015 Lausanne, Switzerland
| |
Collapse
|