1
|
Rishi V, Cole-Filipiak NC, Ramasesha K, McCaslin LM. Excited state electronic structure of dimethyl disulfide involved in photodissociation at ∼200 nm. Phys Chem Chem Phys 2024. [PMID: 39240347 DOI: 10.1039/d4cp02505a] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 09/07/2024]
Abstract
Dimethyl disulfide (DMDS), one of the smallest organic molecules with an S-S bond, serves as a model system for understanding photofragmentation in polypeptides and proteins. Prior studies of DMDS photodissociation excited at ∼266 nm and ∼248 nm have elucidated the mechanisms of S-S and C-S bond cleavage, which involve the lowest excited electronic states S1 and S2. Far less is known about the dissociation mechanisms and electronic structure of relevant excited states of DMDS excited at ∼200 nm. Herein we present calculations of the electronic structure and properties of electronic states S1-S6 accessed when DMDS is excited at ∼200 nm. Our analysis includes a comparison of theoretical and experimental UV spectra, as well as theoretically predicted one-dimensional cuts through the singlet and triplet potential energy surfaces along the S-S and C-S bond dissociation coordinates. Finally, we present calculations of spin-orbit coupling constants at the Franck-Condon geometry to assess the likelihood of ultrafast intersystem crossing. We show that choosing an accurate yet computationally efficient electronic structure method for calculating the S0-S6 potential energy surfaces along relevant dissociation coordinates is challenging due to excited states with doubly excited character and/or mixed Rydberg-valence character. Our findings demonstrate that the extended multi-state complete active space second-order perturbation theory (XMS-CASPT2) balances this computational efficiency and accuracy, as it captures both the Rydberg character of states in the Franck-Condon region and multiconfigurational character toward the bond-dissociation limits. We compare the performance of XMS-CASPT2 to a new variant of equation of motion coupled cluster theory with single, double, and perturbative triple corrections, EOM-CCSD(T)(a)*, finding that EOM-CCSD(T)(a)* significantly improves the treatment of doubly excited states compared to EOM-CCSD, but struggles to quantitatively capture asymptotic energies along bond dissociation coordinates for these states.
Collapse
Affiliation(s)
- Varun Rishi
- Combustion Research Facility, Sandia National Laboratories, Livermore, California 94550, USA.
| | - Neil C Cole-Filipiak
- Combustion Research Facility, Sandia National Laboratories, Livermore, California 94550, USA.
| | - Krupa Ramasesha
- Combustion Research Facility, Sandia National Laboratories, Livermore, California 94550, USA.
| | - Laura M McCaslin
- Combustion Research Facility, Sandia National Laboratories, Livermore, California 94550, USA.
| |
Collapse
|
2
|
Shaalan Alag A, Szalay PG, Tajti A. Ab initio investigation of excited state charge transfer pathways in differently capped bithiophene cages. J Comput Chem 2024; 45:1078-1086. [PMID: 38241483 DOI: 10.1002/jcc.27307] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/25/2023] [Revised: 11/27/2023] [Accepted: 12/15/2023] [Indexed: 01/21/2024]
Abstract
The electronic excitations of conformationally constrained bithiophene cage systems as previously investigated by Lewis et al. (J. Am. Chem. Soc. 143, 18548 (2021)) are revisited, employing the correlated ab initio Scaled Opposite-Spin Algebraic Diagrammatic Construction Second Order electronic structure method. Quantitative descriptors are determined to assess the extent of charge transfer between the bithiophene moieties and the capping domains, represented by either phenyl or triazine groups. The investigation substantiates intrinsic differences in the photophysical behavior of these two structural variants and reveals the presence of lower-energy excited states characterized by noteworthy charge transfer contributions in the triazine cage system. The manifestation of this charge transfer character is discernible even at the Franck-Condon geometry, persisting throughout the relaxation of the excited state. By examining isolated monomer building blocks, we confirm the existence of analogous charge transfer contributions in their excitations. Employing this methodological approach facilitates the prospective identification of potential wall/cap chromophore pairs, wherein charge transfer pathways can be accessed within the energetically favorable regime.
Collapse
Affiliation(s)
- Ahmed Shaalan Alag
- Laboratory of Theoretical Chemistry, Institute of Chemistry, ELTE Eötvös Loránd University, Budapest, Hungary
| | - Péter G Szalay
- Laboratory of Theoretical Chemistry, Institute of Chemistry, ELTE Eötvös Loránd University, Budapest, Hungary
| | - Attila Tajti
- Laboratory of Theoretical Chemistry, Institute of Chemistry, ELTE Eötvös Loránd University, Budapest, Hungary
| |
Collapse
|
3
|
Watrous AG, Davis MC, Fortenberry RC. Performance of EOM-CCSD(T)(a)*-Based Quartic Force Fields in Computing Fundamental, Anharmonic Vibrational Frequencies of Molecular Electronically Excited States with Application to the Ã1A″ State of :CCH 2 (Vinylidene). J Phys Chem A 2024; 128:2150-2161. [PMID: 38466814 DOI: 10.1021/acs.jpca.3c08168] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 03/13/2024]
Abstract
Highly accurate anharmonic vibrational frequencies of electronically excited states are not as easily computed as their ground electronic state counterparts, but recently developed approximate triple excited state methods may be changing that. One emerging excited state method is equation of motion coupled cluster theory at the singles and doubles level with perturbative triples computed via the (a)* formalism, EOMEE-CCSD(T)(a)*. One of the most employed means for the ready computation of vibrational anharmonic frequencies for ground electronic states is second-order vibrational perturbation theory (VPT2), a theory based on quartic force fields (QFFs),fourth-order Taylor series expansions of the potential portion of the internuclear Watson Hamiltonian. The combination of these two is herein benchmarked for its performance for use as a means of computing rovibrational spectra of electronically excited states. Specifically, the EOMEE-CCSD(T)(a)* approach employing a complete basis set extrapolation along with core electron inclusion and relativity (the so-called "CcCR" approach) defining the QFF produces anharmonic fundamental vibrational frequencies within 2.83%, on the average, of reported gas-phase experimentally assigned values for the test set including the A ~ 1 A ″ states of HCF, HCCl, HSiF, HNO, and HPO. However, some states have exceptional accuracy in the fundamentals, most notably for ν2 of A ~ 1 A ″ HCCl in which the CcCR QFF value is within 1.8 cm-1 at 927.9 cm-1 (or 0.2%) of the experiment. Additionally, this approach produces rotational constants to, on the absolute average, within 0.41% of available experimental data, showcasing notable accuracy in the computation of rovibronic spectral data. Furthermore, utilizing a hybrid approach composed of harmonic CcCR force constants along with a set of simple EOMEE-CCSD(T)(a)*/aug-cc-pVQZ QFF cubic and quartic force constants is faster than using pure CcCR and better represents those modes that suffer from numerical instability in the anharmonic portion of the QFF, implying that this so-called "CcCR + QZ" QFF approach may be the best for future applications. Finally, complete, rovibrational spectral data are provided for A ~ 1 A 2 :CCH2, a molecule of potential astrochemical interest, in order to aid in its potential future experimental rovibronic characterization.
Collapse
Affiliation(s)
- Alexandria G Watrous
- Department of Chemistry & Biochemistry, University of Mississippi, University, Mississippi 38677-1848, United States
| | - Megan C Davis
- Department of Chemistry & Biochemistry, University of Mississippi, University, Mississippi 38677-1848, United States
| | - Ryan C Fortenberry
- Department of Chemistry & Biochemistry, University of Mississippi, University, Mississippi 38677-1848, United States
| |
Collapse
|
4
|
Garrett NR, Davis MC, Fortenberry RC. DFT + F12 QFFs for Cost-Effective Rovibrational Spectral Data Predictions of Ground and Excited Electronic States. J Chem Theory Comput 2024. [PMID: 38230913 DOI: 10.1021/acs.jctc.3c01179] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/18/2024]
Abstract
The quest for faster computation of anharmonic vibrational frequencies of both ground and excited electronic states has led to combining coupled cluster theory harmonic force constants with density functional theory cubic and quartic force constants for defining a quartic force field (QFF) utilized in conjunction with vibrational perturbation theory at second order (VPT2). This work shows that explicitly correlated coupled cluster theory at the singles, doubles, and perturbative triples levels [CCSD(T)-F12] provides accurate anharmonic vibrational frequencies and rotational constants when conjoined with any of B3LYP, CAM-B3LYP, BHandHLYP, PBE0, and ωB97XD for roughly one-quarter of the computational time of the CCSD(T)-F12 QFF alone for our test set. As the number of atoms in the molecule increases, however, the anharmonic terms become a greater portion of the QFF, and the cost comparison improves with HOCO+ and formic acid, requiring less than 15 and 10% of the time, respectively. In electronically excited states, PBE0 produces more consistently accurate results. Additionally, as the size of the molecule and, in turn, QFF increase, the cost savings for utilizing such a hybrid approach for both ground- and excited-state computations grows. As such, these methods are promising for predicting accurate rovibrational spectral properties for electronically excited states. In cases where well-behaved potentials for a small selection of targeted excited states are needed, such an approach should reduce the computational cost compared to that of methods requiring semiglobal potential surfaces or variational treatments of the rovibronic Hamiltonian. Such applications include spectral characterization of comets, exoplanets, or any situation in which gas phase molecules are being excited by UV-vis radiation.
Collapse
Affiliation(s)
- Noah R Garrett
- Department of Chemistry & Biochemistry, University of Mississippi, University, Mississippi 38677-1848, United States
| | - Megan C Davis
- Theoretical Division, T-1, Los Alamos National Laboratory, Los Alamos, New Mexico 87545, United States
- Center for Nonlinear Studies, Los Alamos National Laboratory, Los Alamos, New Mexico 87545, United States
| | - Ryan C Fortenberry
- Department of Chemistry & Biochemistry, University of Mississippi, University, Mississippi 38677-1848, United States
| |
Collapse
|
5
|
Parsons SW, Hucek DG, Mishra P, Plusquellic DF, Zwier TS, Drucker S. Jet-Cooled Phosphorescence Excitation Spectrum of the T 1(n,π*) ← S 0 Transition of 4 H-Pyran-4-one. J Phys Chem A 2023; 127:3636-3647. [PMID: 37067071 PMCID: PMC10150392 DOI: 10.1021/acs.jpca.3c01059] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/15/2023] [Revised: 03/31/2023] [Indexed: 04/18/2023]
Abstract
The 4H-pyran-4-one (4PN) molecule is a cyclic conjugated enone with spectroscopically accessible singlet and triplet (n,π*)excited states. Vibronic spectra of 4PN provide a stringent test of electronic-structure calculations, through comparison of predicted vs measured vibrational frequencies in the excited state. We report here the T1(n,π*) ← S0 phosphorescence excitation spectrum of 4PN, recorded under the cooling conditions of a supersonic free-jet expansion. The jet cooling has eliminated congestion appearing in previous room-temperature measurements of the T1 ← S0 band system and has enabled us to determine precise fundamental frequencies for seven vibrational modes of the molecule in its T1(n,π*) state. We have also analyzed the rotational contour of the 000 band, obtaining experimental values for spin-spin and spin-rotation constants of the T1(n,π*) state. We used the experimental results to test predictions from two commonly used computational methods, equation-of-motion excitation energies with dynamical correlation incorporated at the level of coupled cluster singles doubles (EOM-EE-CCSD) and time-dependent density functional theory (TDDFT). We find that each method predicts harmonic frequencies within a few percent of observed fundamentals, for in-plane vibrational modes. However, for out-of-plane modes, each method has specific liabilities that result in frequency errors on the order of 20-30%. The calculations have helped to identify a perturbation from the T2(π,π*) state that leads to unexpected features observed in the T1(n,π*) ← S0 origin band rotational contour.
Collapse
Affiliation(s)
- Sean W. Parsons
- Department
of Chemistry and Biochemistry, University
of Wisconsin-Eau Claire, 105 Garfield Avenue, Eau Claire, Wisconsin 54701, United States
| | - Devon G. Hucek
- Department
of Chemistry and Biochemistry, University
of Wisconsin-Eau Claire, 105 Garfield Avenue, Eau Claire, Wisconsin 54701, United States
| | - Piyush Mishra
- Department
of Chemistry, Purdue University, 560 Oval Drive, West Lafayette, Indiana 47907, United States
| | - David F. Plusquellic
- Applied
Physics Division, National Institute of
Standards and Technology, 325 Broadway Avenue, Boulder, Colorado 80305, United
States
| | - Timothy S. Zwier
- Department
of Chemistry, Purdue University, 560 Oval Drive, West Lafayette, Indiana 47907, United States
| | - Stephen Drucker
- Department
of Chemistry and Biochemistry, University
of Wisconsin-Eau Claire, 105 Garfield Avenue, Eau Claire, Wisconsin 54701, United States
| |
Collapse
|
6
|
Shaalan Alag A, Jelenfi DP, Tajti A, Szalay PG. Accurate Prediction of Vertical Ionization Potentials and Electron Affinities from Spin-Component Scaled CC2 and ADC(2) Models. J Chem Theory Comput 2022; 18:6794-6801. [PMID: 36269873 PMCID: PMC9890482 DOI: 10.1021/acs.jctc.2c00624] [Citation(s) in RCA: 5] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/04/2023]
Abstract
The CC2 and ADC(2) wave function models and their spin-component scaled modifications are adopted for predicting vertical ionization potentials (VIPs) and electron affinities (VEAs). The ionic solutions are obtained as electronic excitations in the continuum orbital formalism, making possible the use of existing, widespread quantum chemistry codes with minimal modifications, in full consistency with the treatment of charge transfer excitations. The performance of different variants is evaluated via benchmark calculations on various sets from previous works, containing small- and medium-sized systems, including the nucleobases. It is shown that with the spin-scaled approximate methods, in particular the scaled opposite-spin variant of the ADC(2) method, the accuracy of EOM-CCSD is achievable at a fraction of the computational cost, also outperforming many common electron propagator approaches.
Collapse
Affiliation(s)
- Ahmed Shaalan Alag
- György
Hevesy Doctoral School, Institute of Chemistry,
ELTE Eötvös Loránd University, H-1117Budapest, Hungary
| | - Dávid P. Jelenfi
- György
Hevesy Doctoral School, Institute of Chemistry,
ELTE Eötvös Loránd University, H-1117Budapest, Hungary
| | - Attila Tajti
- Laboratory
of Theoretical Chemistry, Institute of Chemistry,
ELTE Eötvös Loránd University, P.O. Box 32, H-1518Budapest 112, Hungary,E-mail:
| | - Péter G. Szalay
- Laboratory
of Theoretical Chemistry, Institute of Chemistry,
ELTE Eötvös Loránd University, P.O. Box 32, H-1518Budapest 112, Hungary
| |
Collapse
|
7
|
Barboza CA, Morawski O, Olas J, Gawrys P, Banasiewicz M, Suwinska K, Shova S, Kozankiewicz B, Sobolewski AL. Unravelling the ambiguity of the emission pattern of donor–acceptor salicylaldimines. J Mol Liq 2021. [DOI: 10.1016/j.molliq.2021.117532] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/20/2022]
|
8
|
Köse ME. How to Predict Excited State Geometry by Using Empirical Parameters Obtained from Franck-Condon Analysis of Optical Spectrum. Chemphyschem 2021; 22:2078-2092. [PMID: 34351030 DOI: 10.1002/cphc.202100437] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/07/2021] [Revised: 07/15/2021] [Indexed: 11/09/2022]
Abstract
Excited state geometries of molecules can be calculated with highly reliable wavefunction schemes. Most of such schemes, however, are applicable to small molecules and can hardly be viewed as error-free for excited state geometries. In this study, a theoretical approach is presented in which the excited state geometries of molecules can be predicted by using vibrationally resolved experimental absorption spectrum in combination with the theoretical modelling of vibrational pattern based on Franck-Condon approximation. Huang-Rhys factors have been empirically determined and used as input for revealing the structural changes occurring between the ground and the excited state geometries upon photoexcitation. Naphthalene molecule has been chosen as a test case to show the robustness of the proposed theoretical approach. Predicted 1B2u excited state geometry of the naphthalene has similar but slightly different bond length alternation pattern when compared with the geometries calculated with CIS, B3LYP, and CC2 methods. Excited state geometries of perylene and pyrene molecules are also determined with the presented theoretical approach. This powerful method can be applied to other molecules and specifically to relatively large molecules rather easily as long as vibrationally resolved experimental spectra are available to use.
Collapse
|
9
|
Jankowska J, Sobolewski AL. Modern Theoretical Approaches to Modeling the Excited-State Intramolecular Proton Transfer: An Overview. Molecules 2021; 26:molecules26175140. [PMID: 34500574 PMCID: PMC8434569 DOI: 10.3390/molecules26175140] [Citation(s) in RCA: 22] [Impact Index Per Article: 7.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/15/2021] [Revised: 08/20/2021] [Accepted: 08/23/2021] [Indexed: 02/02/2023] Open
Abstract
The excited-state intramolecular proton transfer (ESIPT) phenomenon is nowadays widely acknowledged to play a crucial role in many photobiological and photochemical processes. It is an extremely fast transformation, often taking place at sub-100 fs timescales. While its experimental characterization can be highly challenging, a rich manifold of theoretical approaches at different levels is nowadays available to support and guide experimental investigations. In this perspective, we summarize the state-of-the-art quantum-chemical methods, as well as molecular- and quantum-dynamics tools successfully applied in ESIPT process studies, focusing on a critical comparison of their specific properties.
Collapse
Affiliation(s)
- Joanna Jankowska
- Faculty of Chemistry, University of Warsaw, 02-093 Warsaw, Poland
- Correspondence:
| | | |
Collapse
|
10
|
Davis MC, Fortenberry RC. (T)+EOM Quartic Force Fields for Theoretical Vibrational Spectroscopy of Electronically Excited States. J Chem Theory Comput 2021; 17:4374-4382. [PMID: 34165980 DOI: 10.1021/acs.jctc.1c00307] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
(T)+EOM quartic force fields (QFFs) are proposed for ab initio rovibrational properties of electronically excited states of small molecules. The (T)+EOM method is a simple treatment of the potential surface of the excited state using a composite energy from the CCSD(T) energy for the ground-state configuration and the EOM-CCSD excitation energy for the target state. The method is benchmarked with two open-shell species, HOO and HNF, and two closed-shell species, HNO and HCF. A (T)+EOM QFF with a complete basis set extrapolation (C) and corrections for core correlation (cC) and scalar relativity (R), dubbed (T)+EOM/CcCR, achieves a mean absolute error (MAE) as low as 1.6 cm-1 for the à 2A' state of HOO versus an established benchmark QFF with CCSD(T)-F12b/cc-pVTZ-F12 (F12-TZ) for this variationally accessible electronically excited state. The MAE for anharmonic frequencies for (T)+EOM/CcCR versus F12-TZ for HNF is 7.5 cm-1. The closed-shell species are compared directly with the experiment, where a simpler (T)+EOM QFF using the aug-cc-pVTZ basis set compares more favorably than the more costly (T)+EOM/CcCR, suggesting a possible influence of decreasing accuracy with basis set size. Scans along internal coordinates are also provided which show reasonable modeling of the potential surface by (T)+EOM compared to benchmark QFFs computed for variationally accessible electronic states. The agreement between (T)+EOM/CcCR with F12-TZ and CcCR benchmarks is also shown to be quite accurate for rotational constants and geometries, with an MAE of 0.008 MHz for the rotational constants of (T)+EOM/CcCR versus CcCR for à 2A' HOO and agreement within 0.003 Šfor bond lengths.
Collapse
Affiliation(s)
- Megan C Davis
- Department of Chemistry & Biochemistry, University of Mississippi, University, Mississippi 38677-1848, United States
| | - Ryan C Fortenberry
- Department of Chemistry & Biochemistry, University of Mississippi, University, Mississippi 38677-1848, United States
| |
Collapse
|
11
|
Gozem S, Krylov AI. The
ezSpectra
suite: An easy‐to‐use toolkit for spectroscopy modeling. WILEY INTERDISCIPLINARY REVIEWS-COMPUTATIONAL MOLECULAR SCIENCE 2021. [DOI: 10.1002/wcms.1546] [Citation(s) in RCA: 26] [Impact Index Per Article: 8.7] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/10/2022]
Affiliation(s)
- Samer Gozem
- Department of Chemistry Georgia State University Atlanta Georgia USA
| | - Anna I. Krylov
- Department of Chemistry University of Southern California Los Angeles California USA
| |
Collapse
|
12
|
Tajti A, Kozma B, Szalay PG. Improved Description of Charge-Transfer Potential Energy Surfaces via Spin-Component-Scaled CC2 and ADC(2) Methods. J Chem Theory Comput 2021; 17:439-449. [PMID: 33326229 DOI: 10.1021/acs.jctc.0c01146] [Citation(s) in RCA: 11] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
The molecular level understanding of electronic transport properties depends on the reliable theoretical description of charge-transfer (CT)-type electronic states. In this paper, the performance of spin-component-scaled variants of the popular CC2 and ADC(2) methods is evaluated for CT states, following benchmark strategies of earlier studies that revealed a compromised accuracy of the unmodified models. In addition to statistics on the accuracy of vertical excitation energies at equilibrium and infinite separation of bimolecular complexes, potential energy surfaces of the ammonia-fluorine complex are also reported. The results show the capability of spin-component-scaled approaches to reduce the large errors of their regular counterparts to a significant extent, outperforming even the coupled-cluster single and double method in many cases. The cost-effective scaled-opposite-spin variants are found to provide a remarkably good agreement with the CCSDT-3 reference data, thereby being recommended methods of choice in the study of charge-transfer states.
Collapse
Affiliation(s)
- Attila Tajti
- Institute of Chemistry, Laboratory of Theoretical Chemistry, ELTE Eötvös Loránd University, P.O. Box 32, H-1518, Budapest 112, Hungary
| | - Balázs Kozma
- Institute of Chemistry, Laboratory of Theoretical Chemistry, ELTE Eötvös Loránd University, P.O. Box 32, H-1518, Budapest 112, Hungary
| | - Péter G Szalay
- Institute of Chemistry, Laboratory of Theoretical Chemistry, ELTE Eötvös Loránd University, P.O. Box 32, H-1518, Budapest 112, Hungary
| |
Collapse
|
13
|
Abstract
We present a new and efficient implementation of the closed shell coupled cluster singles and doubles with perturbative triples method (CC3) in the electronic structure program eT. Asymptotically, a ground state calculation has an iterative cost of 4nV4nO3 floating point operations (FLOP), where nV and nO are the number of virtual and occupied orbitals, respectively. The Jacobian and transpose Jacobian transformations, required to iteratively solve for excitation energies and transition moments, both require 8nV4nO3 FLOP. We have also implemented equation of motion (EOM) transition moments for CC3. The EOM transition densities require recalculation of triples amplitudes, as nV3nO3 tensors are not stored in memory. This results in a noniterative computational cost of 10nV4nO3 FLOP for the ground state density and 26nV4nO3 FLOP per state for the transition densities. The code is compared to the CC3 implementations in CFOUR, DALTON, and PSI4. We demonstrate the capabilities of our implementation by calculating valence and core excited states of l-proline.
Collapse
Affiliation(s)
- Alexander
C. Paul
- Department
of Chemistry, Norwegian University of Science
and Technology, NTNU, 7491 Trondheim, Norway
| | - Rolf H. Myhre
- Department
of Chemistry, Norwegian University of Science
and Technology, NTNU, 7491 Trondheim, Norway
| | - Henrik Koch
- Department
of Chemistry, Norwegian University of Science
and Technology, NTNU, 7491 Trondheim, Norway
- Scuola
Normale Superiore, Piazza dei Cavaleri 7, 56126 Pisa, Italy
| |
Collapse
|
14
|
Rishi V, Perera A, Bartlett RJ. A route to improving RPA excitation energies through its connection to equation-of-motion coupled cluster theory. J Chem Phys 2020; 153:234101. [DOI: 10.1063/5.0023862] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/13/2022] Open
Affiliation(s)
- Varun Rishi
- Chemistry and Chemical Engineering, California Institute of Technology, Pasadena, California 91125, USA
| | - Ajith Perera
- Quantum Theory Project, University of Florida, Gainesville, Florida 32611, USA
| | - Rodney J. Bartlett
- Quantum Theory Project, University of Florida, Gainesville, Florida 32611, USA
| |
Collapse
|
15
|
Chrayteh A, Blondel A, Loos PF, Jacquemin D. Mountaineering Strategy to Excited States: Highly Accurate Oscillator Strengths and Dipole Moments of Small Molecules. J Chem Theory Comput 2020; 17:416-438. [DOI: 10.1021/acs.jctc.0c01111] [Citation(s) in RCA: 12] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/13/2023]
Affiliation(s)
- Amara Chrayteh
- Université de Nantes, CNRS, CEISAM UMR 6230, F-44000 Nantes, France
| | - Aymeric Blondel
- Université de Nantes, CNRS, CEISAM UMR 6230, F-44000 Nantes, France
| | - Pierre-François Loos
- Laboratoire de Chimie et Physique Quantiques, Université de Toulouse, CNRS, UPS, 31062 Toulouse, France
| | - Denis Jacquemin
- Université de Nantes, CNRS, CEISAM UMR 6230, F-44000 Nantes, France
| |
Collapse
|
16
|
Kozma B, Tajti A, Demoulin B, Izsák R, Nooijen M, Szalay PG. A New Benchmark Set for Excitation Energy of Charge Transfer States: Systematic Investigation of Coupled Cluster Type Methods. J Chem Theory Comput 2020; 16:4213-4225. [PMID: 32502351 PMCID: PMC7467641 DOI: 10.1021/acs.jctc.0c00154] [Citation(s) in RCA: 42] [Impact Index Per Article: 10.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
![]()
The
numerous existing publications on benchmarking quantum chemistry
methods for excited states rarely include Charge Transfer (CT) states,
although many interesting phenomena in, e.g., biochemistry and material
physics involve the transfer of electrons between fragments of the
system. Therefore, it is timely to test the accuracy of quantum chemical
methods for CT states, as well. In this study we first propose a new
benchmark set consisting of dimers having low-energy CT states. On
this set, the vertical excitation energy has been calculated with
Coupled Cluster methods including triple excitations (CC3, CCSDT-3,
CCSD(T)(a)*), as well as with methods including full or approximate
doubles (CCSD, STEOM-CCSD, CC2, ADC(2), EOM-CCSD(2)). The results
show that the popular CC2 and ADC(2) methods are much less accurate
for CT states than for valence states. On the other hand, EOM-CCSD
seems to have similar systematic overestimation of the excitation
energies for both types of states. Among the triples methods the novel
EOM-CCSD(T)(a)* method including noniterative triple excitations is
found to stand out with its consistently good performance for all
types of states, delivering essentially EOM-CCSDT quality results.
Collapse
Affiliation(s)
- Balázs Kozma
- Institute of Chemistry, Laboratory of Theoretical Chemistry, ELTE Eötvös Loránd University, P.O. Box 32, H-1518, Budapest 112, Hungary
| | - Attila Tajti
- Institute of Chemistry, Laboratory of Theoretical Chemistry, ELTE Eötvös Loránd University, P.O. Box 32, H-1518, Budapest 112, Hungary
| | - Baptiste Demoulin
- Max-Planck-Institut für Kohlenforschung, 45470 Mülheim an der Ruhr, Germany
| | - Róbert Izsák
- Max-Planck-Institut für Kohlenforschung, 45470 Mülheim an der Ruhr, Germany
| | - Marcel Nooijen
- Department of Chemistry, University of Waterloo, Waterloo, Ontario N2L 3G1, Canada
| | - Péter G Szalay
- Institute of Chemistry, Laboratory of Theoretical Chemistry, ELTE Eötvös Loránd University, P.O. Box 32, H-1518, Budapest 112, Hungary
| |
Collapse
|
17
|
Matthews DA, Cheng L, Harding ME, Lipparini F, Stopkowicz S, Jagau TC, Szalay PG, Gauss J, Stanton JF. Coupled-cluster techniques for computational chemistry: The CFOUR program package. J Chem Phys 2020; 152:214108. [DOI: 10.1063/5.0004837] [Citation(s) in RCA: 214] [Impact Index Per Article: 53.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/30/2022] Open
Affiliation(s)
- Devin A. Matthews
- Department of Chemistry, Southern Methodist University, Dallas, Texas 75275, USA
| | - Lan Cheng
- Department of Chemistry, Johns Hopkins University, Baltimore, Maryland 21218, USA
| | - Michael E. Harding
- Institut für Physikalische Chemie, Karlsruher Institut für Technologie (KIT), Kaiserstr. 12, D-76131 Karlsruhe, Germany
| | - Filippo Lipparini
- Dipartimento di Chimica e Chimica Industriale, Università di Pisa, Via G. Moruzzi 13, I-56124 Pisa, Italy
| | - Stella Stopkowicz
- Department Chemie, Johannes Gutenberg-Universität Mainz, Duesbergweg 10-14, D-55128 Mainz, Germany
| | - Thomas-C. Jagau
- Department of Chemistry, University of Munich (LMU), Butenandtstr. 5-13, D-81377 Munich, Germany
| | - Péter G. Szalay
- ELTE Eötvös Loránd University, Institute of Chemistry, Pázmány Péter sétány 1/A, H-1117 Budapest, Hungary
| | - Jürgen Gauss
- Department Chemie, Johannes Gutenberg-Universität Mainz, Duesbergweg 10-14, D-55128 Mainz, Germany
| | - John F. Stanton
- Quantum Theory Project, Departments of Chemistry and Physics, University of Florida, Gainesville, Florida 32611, USA
| |
Collapse
|
18
|
Loos PF, Scemama A, Jacquemin D. The Quest for Highly Accurate Excitation Energies: A Computational Perspective. J Phys Chem Lett 2020; 11:2374-2383. [PMID: 32125872 DOI: 10.1021/acs.jpclett.0c00014] [Citation(s) in RCA: 87] [Impact Index Per Article: 21.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/11/2023]
Abstract
We provide an overview of the successive steps that made it possible to obtain increasingly accurate excitation energies with computational chemistry tools, eventually leading to chemically accurate vertical transition energies for small- and medium-size molecules. First, we describe the evolution of ab initio methods employed to define benchmark values, with the original Roos CASPT2 method, then the CC3 method as in the renowned Thiel set, and more recently the resurgence of selected configuration interaction methods. The latter method has been able to deliver consistently, for both single and double excitations, highly accurate excitation energies for small molecules, as well as medium-size molecules with compact basis sets. Second, we describe how these high-level methods and the creation of representative benchmark sets of excitation energies have allowed the fair and accurate assessment of the performance of computationally lighter methods. We conclude by discussing possible future theoretical and technological developments in the field.
Collapse
Affiliation(s)
- Pierre-François Loos
- Laboratoire de Chimie et Physique Quantiques, Université de Toulouse, CNRS, UPS, 31062 Toulouse, France
| | - Anthony Scemama
- Laboratoire de Chimie et Physique Quantiques, Université de Toulouse, CNRS, UPS, 31062 Toulouse, France
| | - Denis Jacquemin
- Université de Nantes, CNRS, CEISAM UMR 6230, F-44000 Nantes, France
| |
Collapse
|
19
|
Barboza CA, Gawrys P, Banasiewicz M, Suwinska K, Sobolewski AL. Photophysical transformations induced by chemical substitution to salicylaldimines. Phys Chem Chem Phys 2020; 22:6698-6705. [DOI: 10.1039/d0cp00110d] [Citation(s) in RCA: 18] [Impact Index Per Article: 4.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/02/2023]
Abstract
The role of electron-acceptor strength and microenvironment polarity on the photophysical properties of salicylaldimines.
Collapse
Affiliation(s)
| | - Pawel Gawrys
- Physics Institute
- Polish Academy of Sciences
- Warsaw
- Poland
| | | | - Kinga Suwinska
- Faculty of Mathematics and Natural Sciences
- Cardinal Stefan Wyszynski University in Warsaw
- Warsaw
- Poland
| | | |
Collapse
|
20
|
Tajti A, Tulipán L, Szalay PG. Accuracy of Spin-Component Scaled ADC(2) Excitation Energies and Potential Energy Surfaces. J Chem Theory Comput 2019; 16:468-474. [DOI: 10.1021/acs.jctc.9b01065] [Citation(s) in RCA: 16] [Impact Index Per Article: 3.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/26/2023]
Affiliation(s)
- Attila Tajti
- Laboratory of Theoretical Chemistry, Institute of Chemistry, ELTE Eötvös Loránd University, 112, P.O. Box 32, H-1518 Budapest, Hungary
| | - Levente Tulipán
- Laboratory of Theoretical Chemistry, Institute of Chemistry, ELTE Eötvös Loránd University, 112, P.O. Box 32, H-1518 Budapest, Hungary
| | - Péter G. Szalay
- Laboratory of Theoretical Chemistry, Institute of Chemistry, ELTE Eötvös Loránd University, 112, P.O. Box 32, H-1518 Budapest, Hungary
| |
Collapse
|
21
|
Cheng L. A study of non-iterative triples contributions in relativistic equation-of-motion coupled-cluster calculations using an exact two-component Hamiltonian with atomic mean-field spin-orbit integrals: Application to uranyl and other heavy-element compounds. J Chem Phys 2019; 151:104103. [DOI: 10.1063/1.5113796] [Citation(s) in RCA: 11] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Lan Cheng
- Department of Chemistry, The Johns Hopkins University, Baltimore, Maryland 21218, USA
| |
Collapse
|
22
|
Tajti A, Szalay PG. Accuracy of Spin-Component-Scaled CC2 Excitation Energies and Potential Energy Surfaces. J Chem Theory Comput 2019; 15:5523-5531. [DOI: 10.1021/acs.jctc.9b00676] [Citation(s) in RCA: 20] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Attila Tajti
- ELTE Eötvös Loránd University, Laboratory of Theoretical Chemistry, Pázmány Péter sétány 1/A, H-1117 Budapest, Hungary
| | - 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, H-1117 Budapest, Hungary
| |
Collapse
|
23
|
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]
|