1
|
Yoshida Y, Takemori N, Mizukami W. Ab initio extended Hubbard model of short polyenes for efficient quantum computing. J Chem Phys 2024; 161:084303. [PMID: 39193941 DOI: 10.1063/5.0213525] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/11/2024] [Accepted: 08/05/2024] [Indexed: 08/29/2024] Open
Abstract
We propose introducing an extended Hubbard Hamiltonian derived via the ab initio downfolding method, which was originally formulated for periodic materials, toward efficient quantum computing of molecular electronic structure calculations. By utilizing this method, the first-principles Hamiltonian of chemical systems can be coarse-grained by eliminating the electronic degrees of freedom in higher energy space and reducing the number of terms of electron repulsion integral from O(N4) to O(N2). Our approach is validated numerically on the vertical excitation energies and excitation characters of ethylene, butadiene, and hexatriene. The dynamical electron correlation is incorporated within the framework of the constrained random phase approximation in advance of quantum computations, and the constructed models capture the trend of experimental and high-level quantum chemical calculation results. As expected, the L1-norm of the fermion-to-qubit mapped model Hamiltonians is significantly lower than that of conventional ab initio Hamiltonians, suggesting improved scalability of quantum computing. Those numerical outcomes and the results of the simulation of excited-state sampling demonstrate that the ab initio extended Hubbard Hamiltonian may hold significant potential for quantum chemical calculations using quantum computers.
Collapse
Affiliation(s)
- Yuichiro Yoshida
- Center for Quantum Information and Quantum Biology, Osaka University, 1-2 Machikaneyama, Toyonaka, Osaka 560-0043, Japan
| | - Nayuta Takemori
- Center for Quantum Information and Quantum Biology, Osaka University, 1-2 Machikaneyama, Toyonaka, Osaka 560-0043, Japan
- Center for Emergent Matter Science, RIKEN, Wako, Saitama 351-0198, Japan
| | - Wataru Mizukami
- Center for Quantum Information and Quantum Biology, Osaka University, 1-2 Machikaneyama, Toyonaka, Osaka 560-0043, Japan
- Graduate School of Engineering Science, Osaka University, 1-3 Machikaneyama, Toyonaka, Osaka 560-8531, Japan
| |
Collapse
|
2
|
Glebov IO, Poddubnyy VV, Khokhlov D. Perturbation theory in the complete degenerate active space (CDAS-PT2). J Chem Phys 2024; 161:024114. [PMID: 38995081 DOI: 10.1063/5.0211210] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/29/2024] [Accepted: 06/17/2024] [Indexed: 07/13/2024] Open
Abstract
Methods based on the multireference perturbation theory (MRPT) with the one-electron zeroth-order Hamiltonian are widely used for the description of excited states, for example, due to their relatively low computational cost. However, current methods have a common drawback-use of a model space with low size. In this article, we propose the MRPT method with the model space extended to the complete active space. The one-electron zeroth-order Hamiltonian suitable for this extension is formulated. The proposed method was applied to common models, such as LiF, ethylene, and trans-butadiene. It was shown to have accuracy superior to XMCQDPT2 in most cases, especially in the case of the small active space.
Collapse
Affiliation(s)
- Ilya O Glebov
- Chemistry Department, Moscow State University, Leninskie Gory 1-3, Moscow 119991, Russia
| | - Vladimir V Poddubnyy
- Chemistry Department, Moscow State University, Leninskie Gory 1-3, Moscow 119991, Russia
| | | |
Collapse
|
3
|
Kossoski F, Boggio-Pasqua M, Loos PF, Jacquemin D. Reference Energies for Double Excitations: Improvement and Extension. J Chem Theory Comput 2024; 20:5655-5678. [PMID: 38885174 DOI: 10.1021/acs.jctc.4c00410] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/20/2024]
Abstract
In the realm of photochemistry, the significance of double excitations (also known as doubly excited states), where two electrons are concurrently elevated to higher energy levels, lies in their involvement in key electronic transitions essential in light-induced chemical reactions as well as their challenging nature from the computational theoretical chemistry point of view. Based on state-of-the-art electronic structure methods (such as high-order coupled-cluster, selected configuration interaction, and multiconfigurational methods), we improve and expand our prior set of accurate reference excitation energies for electronic states exhibiting a substantial amount of double excitations [Loos et al. J. Chem. Theory Comput. 2019, 15, 1939]. This extended collection encompasses 47 electronic transitions across 26 molecular systems that we separate into two distinct subsets: (i) 28 "genuine" doubly excited states where the transitions almost exclusively involve doubly excited configurations and (ii) 19 "partial" doubly excited states which exhibit a more balanced character between singly and doubly excited configurations. For each subset, we assess the performance of high-order coupled-cluster (CC3, CCSDT, CC4, and CCSDTQ) and multiconfigurational methods (CASPT2, CASPT3, PC-NEVPT2, and SC-NEVPT2). Using as a probe the percentage of single excitations involved in a given transition (%T1) computed at the CC3 level, we also propose a simple correction that reduces the errors of CC3 by a factor of 3, for both sets of excitations. We hope that this more complete and diverse compilation of double excitations will help future developments of electronic excited-state methodologies.
Collapse
Affiliation(s)
- Fábris Kossoski
- Laboratoire de Chimie et Physique Quantiques (UMR 5626), Université de Toulouse, CNRS, UPS, 31062 Toulouse, France
| | - Martial Boggio-Pasqua
- Laboratoire de Chimie et Physique Quantiques (UMR 5626), Université de Toulouse, CNRS, UPS, 31062 Toulouse, France
| | - Pierre-François Loos
- Laboratoire de Chimie et Physique Quantiques (UMR 5626), Université de Toulouse, CNRS, UPS, 31062 Toulouse, France
| | - Denis Jacquemin
- Nantes Université, CNRS, CEISAM UMR 6230, F-44000 Nantes, France
- Institut Universitaire de France (IUF), F-75005 Paris, France
| |
Collapse
|
4
|
Marie A, Loos PF. Reference Energies for Valence Ionizations and Satellite Transitions. J Chem Theory Comput 2024; 20:4751-4777. [PMID: 38776293 PMCID: PMC11171335 DOI: 10.1021/acs.jctc.4c00216] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/21/2024] [Revised: 04/10/2024] [Accepted: 04/11/2024] [Indexed: 05/24/2024]
Abstract
Upon ionization of an atom or a molecule, another electron (or more) can be simultaneously excited. These concurrently generated states are called "satellites" (or shakeup transitions) as they appear in ionization spectra as higher-energy peaks with weaker intensity and larger width than the main peaks associated with single-particle ionizations. Satellites, which correspond to electronically excited states of the cationic species, are notoriously challenging to model using conventional single-reference methods due to their high excitation degree compared to the neutral reference state. This work reports 42 satellite transition energies and 58 valence ionization potentials (IPs) of full configuration interaction quality computed in small molecular systems. Following the protocol developed for the quest database [Véril, M.; Scemama, A.; Caffarel, M.; Lipparini, F.; Boggio-Pasqua, M.; Jacquemin, D.; and Loos, P.-F. Wiley Interdiscip. Rev.: Comput. Mol. Sci. 2021, 11, e1517], these reference energies are computed using the configuration interaction using a perturbative selection made iteratively (CIPSI) method. In addition, the accuracy of the well-known coupled-cluster (CC) hierarchy (CC2, CCSD, CC3, CCSDT, CC4, and CCSDTQ) is gauged against these new accurate references. The performances of various approximations based on many-body Green's functions (GW, GF2, and T-matrix) for IPs are also analyzed. Their limitations in correctly modeling satellite transitions are discussed.
Collapse
Affiliation(s)
- Antoine Marie
- Laboratoire de Chimie et Physique
Quantiques (UMR 5626), Université de Toulouse, CNRS, UPS, Toulouse 31062, France
| | - Pierre-François Loos
- Laboratoire de Chimie et Physique
Quantiques (UMR 5626), Université de Toulouse, CNRS, UPS, Toulouse 31062, France
| |
Collapse
|
5
|
Folkestad S, Paul AC, Paul née Matveeva R, Reinholdt P, Coriani S, Odelius M, Koch H. Quantum Mechanical Versus Polarizable Embedding Schemes: A Study of the Xray Absorption Spectra of Aqueous Ammonia and Ammonium. J Chem Theory Comput 2024; 20:4161-4169. [PMID: 38713524 PMCID: PMC11137810 DOI: 10.1021/acs.jctc.4c00105] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/26/2024] [Revised: 04/16/2024] [Accepted: 04/17/2024] [Indexed: 05/09/2024]
Abstract
The X-ray absorption spectra of aqueous ammonia and ammonium are computed using a combination of coupled cluster singles and doubles (CCSD) with different quantum mechanical and molecular mechanical embedding schemes. Specifically, we compare frozen Hartree-Fock (HF) density embedding, polarizable embedding (PE), and polarizable density embedding (PDE). Integrating CCSD with frozen HF density embedding is possible within the CC-in-HF framework, which circumvents the conventional system-size limitations of standard coupled cluster methods. We reveal similarities between PDE and frozen HF density descriptions, while PE spectra differ significantly. By including approximate triple excitations, we also investigate the effect of improving the electronic structure theory. The spectra computed using this approach show an improved intensity ratio compared to CCSD-in-HF. Charge transfer analysis of the excitations shows the local character of the pre-edge and main-edge, while the post-edge is formed by excitations delocalized over the first solvation shell and beyond.
Collapse
Affiliation(s)
- Sarai
Dery Folkestad
- Department
of Chemistry, Norwegian University of Science
and Technology, NTNU, 7491 Trondheim, Norway
| | - Alexander C. Paul
- Department
of Chemistry, Norwegian University of Science
and Technology, NTNU, 7491 Trondheim, Norway
| | - Regina Paul née Matveeva
- Department
of Chemistry, Norwegian University of Science
and Technology, NTNU, 7491 Trondheim, Norway
| | - Peter Reinholdt
- Department
of Physics, Chemistry and Pharmacy, University
of Southern Denmark, SDU, Campusvej 55, 5230 Odense, Denmark
| | - Sonia Coriani
- Department
of Chemistry, Technical University of Denmark,
DTU, Kemitorvet Bldg 207, 2800 Kongens Lyngby, Denmark
| | - Michael Odelius
- Department
of Physics, Stockholm University, 10691 Stockholm, Sweden
| | - Henrik Koch
- Department
of Chemistry, Norwegian University of Science
and Technology, NTNU, 7491 Trondheim, Norway
| |
Collapse
|
6
|
Folkestad SD, Paul AC, Paul Née Matveeva R, Coriani S, Odelius M, Iannuzzi M, Koch H. Understanding X-ray absorption in liquid water using triple excitations in multilevel coupled cluster theory. Nat Commun 2024; 15:3551. [PMID: 38670938 PMCID: PMC11053016 DOI: 10.1038/s41467-024-47690-x] [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/04/2023] [Accepted: 04/03/2024] [Indexed: 04/28/2024] Open
Abstract
X-ray absorption (XA) spectroscopy is an essential experimental tool to investigate the local structure of liquid water. Interpretation of the experiment poses a significant challenge and requires a quantitative theoretical description. High-quality theoretical XA spectra require reliable molecular dynamics simulations and accurate electronic structure calculations. Here, we present the first successful application of coupled cluster theory to model the XA spectrum of liquid water. We overcome the computational limitations on system size by employing a multilevel coupled cluster framework for large molecular systems. Excellent agreement with the experimental spectrum is achieved by including triple excitations in the wave function and using molecular structures from state-of-the-art path-integral molecular dynamics. We demonstrate that an accurate description of the electronic structure within the first solvation shell is sufficient to successfully model the XA spectrum of liquid water within the multilevel framework. Furthermore, we present a rigorous charge transfer analysis of the XA spectrum, which is reliable due to the accuracy and robustness of the electronic structure methodology. This analysis aligns with previous studies regarding the character of the prominent features of the XA spectrum of liquid water.
Collapse
Affiliation(s)
- Sarai Dery Folkestad
- Department of Chemistry, Norwegian University of Science and Technology, NTNU, 7491, Trondheim, Norway
| | - Alexander C Paul
- Department of Chemistry, Norwegian University of Science and Technology, NTNU, 7491, Trondheim, Norway
| | - Regina Paul Née Matveeva
- Department of Chemistry, Norwegian University of Science and Technology, NTNU, 7491, Trondheim, Norway
| | - Sonia Coriani
- Department of Chemistry, Technical University of Denmark, DTU, 2800, Kongens Lyngby, Denmark
| | - Michael Odelius
- Department of Physics, Stockholm University, 10691, Stockholm, Sweden
| | - Marcella Iannuzzi
- Department of Chemistry, University of Zurich, 8057, Zürich, Switzerland
| | - Henrik Koch
- Department of Chemistry, Norwegian University of Science and Technology, NTNU, 7491, Trondheim, Norway.
| |
Collapse
|
7
|
Janoš J, Figueira Nunes JP, Hollas D, Slavíček P, Curchod BFE. Predicting the photodynamics of cyclobutanone triggered by a laser pulse at 200 nm and its MeV-UED signals-A trajectory surface hopping and XMS-CASPT2 perspective. J Chem Phys 2024; 160:144305. [PMID: 38591685 DOI: 10.1063/5.0203105] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/09/2024] [Accepted: 03/14/2024] [Indexed: 04/10/2024] Open
Abstract
This work is part of a prediction challenge that invited theoretical/computational chemists to predict the photochemistry of cyclobutanone in the gas phase, excited at 200 nm by a laser pulse, and the expected signal that will be recorded during a time-resolved megaelectronvolt ultrafast electron diffraction (MeV-UED). We present here our theoretical predictions based on a combination of trajectory surface hopping with XMS-CASPT2 (for the nonadiabatic molecular dynamics) and Born-Oppenheimer molecular dynamics with MP2 (for the athermal ground-state dynamics following internal conversion), coined (NA+BO)MD. The initial conditions were sampled from Born-Oppenheimer molecular dynamics coupled to a quantum thermostat. Our simulations indicate that the main photoproducts after 2 ps of dynamics are CO + cyclopropane (50%), CO + propene (10%), and ethene and ketene (34%). The photoexcited cyclobutanone in its second excited electronic state S2 can follow two pathways for its nonradiative decay: (i) a ring-opening in S2 and a subsequent rapid decay to the ground electronic state, where the photoproducts are formed, or (ii) a transfer through a closed-ring conical intersection to S1, where cyclobutanone ring opens and then funnels to the ground state. Lifetimes for the photoproduct and electronic populations were determined. We calculated a stationary MeV-UED signal [difference pair distribution function-ΔPDF(r)] for each (interpolated) pathway as well as a time-resolved signal [ΔPDF(r,t) and ΔI/I(s,t)] for the full swarm of (NA+BO)MD trajectories. Furthermore, our analysis provides time-independent basis functions that can be used to fit the time-dependent experimental UED signals [both ΔPDF(r,t) and ΔI/I(s,t)] and potentially recover the population of photoproducts. We also offer a detailed analysis of the limitations of our model and their potential impact on the predicted experimental signals.
Collapse
Affiliation(s)
- Jiří Janoš
- Department of Physical Chemistry, University of Chemistry and Technology, Technická 5, Prague 6 166 28, Czech Republic
- Centre for Computational Chemistry, School of Chemistry, University of Bristol, Bristol BS8 1TS, United Kingdom
| | | | - Daniel Hollas
- Centre for Computational Chemistry, School of Chemistry, University of Bristol, Bristol BS8 1TS, United Kingdom
| | - Petr Slavíček
- Department of Physical Chemistry, University of Chemistry and Technology, Technická 5, Prague 6 166 28, Czech Republic
| | - Basile F E Curchod
- Centre for Computational Chemistry, School of Chemistry, University of Bristol, Bristol BS8 1TS, United Kingdom
| |
Collapse
|
8
|
Kitsaras MP, Grazioli L, Stopkowicz S. The approximate coupled-cluster methods CC2 and CC3 in a finite magnetic field. J Chem Phys 2024; 160:094112. [PMID: 38441261 DOI: 10.1063/5.0189350] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/28/2023] [Accepted: 01/18/2024] [Indexed: 08/16/2024] Open
Abstract
In this paper, we report on the implementation of CC2 and CC3 in the context of molecules in finite magnetic fields. The methods are applied to the investigation of atoms and molecules through spectroscopic predictions and geometry optimizations for the study of the atmosphere of highly magnetized White Dwarf stars. We show that ground-state finite-field (ff) CC2 is a reasonable alternative to CCSD for energies and, in particular, for geometrical properties. For excited states, ff-CC2 is shown to perform well for states with predominant single-excitation character. Yet, for cases in which the excited state wavefunction has double-excitation character with respect to the reference, ff-CC2 can easily lead to completely unphysical results. Ff-CC3, however, is shown to reproduce the CCSDT behavior very well and enables the treatment of larger systems at a high accuracy.
Collapse
Affiliation(s)
- Marios-Petros Kitsaras
- Fachrichtung Chemie, Universität des Saarlandes, Campus B2.2, D-66123 Saarbrücken, Germany
- Department Chemie, Johannes Gutenberg-Universität Mainz, Duesbergweg 10-14, D-55128 Mainz, Germany
| | - Laura Grazioli
- Department Chemie, Johannes Gutenberg-Universität Mainz, Duesbergweg 10-14, D-55128 Mainz, Germany
| | - Stella Stopkowicz
- Fachrichtung Chemie, Universität des Saarlandes, Campus B2.2, D-66123 Saarbrücken, Germany
- Department Chemie, Johannes Gutenberg-Universität Mainz, Duesbergweg 10-14, D-55128 Mainz, Germany
- Hylleraas Centre for Quantum Molecular Sciences, Department of Chemistry, University of Oslo, P.O. Box 1033 Blindern, N-0315 Oslo, Norway
| |
Collapse
|
9
|
Hait D, Martínez TJ. Predicting the X-ray Absorption Spectrum of Ozone with Single Configuration State Functions. J Chem Theory Comput 2024; 20:873-881. [PMID: 38175153 DOI: 10.1021/acs.jctc.3c01035] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/05/2024]
Abstract
X-ray absorption spectra (XAS) of biradicaloid species are often thought to represent a challenge to theoretical methods. This has led to the testing of recently developed multireference techniques on the XAS of ozone, but reproduction of the experimental spectral profile has proven difficult. We utilize a minimal model consisting of a single configuration state function (CSF) per excited state to model core-level excitations of ozone, with the orbitals of each CSF optimized using the restricted open-shell Kohn-Sham (ROKS) method. This protocol leads to semiquantitative agreement with experimental XAS. In fact, we find that low-lying core-hole excited states in biradicaloids can be approximated with individual CSFs, despite the presence of multireference character in the ground state. We also report that the 1s → π* and 1s → σ* transitions have quite distinct widths for O3. This reveals the importance of sampling over a representative range of geometries from the vibrational ground state for properly assessing the accuracy of electronic structure methods against experiments instead of the popular procedure of uniformly broadening stick spectra at the equilibrium geometry.
Collapse
Affiliation(s)
- Diptarka Hait
- Department of Chemistry and The PULSE Institute, Stanford University, Stanford, California 94305, United States
| | - Todd J Martínez
- Department of Chemistry and The PULSE Institute, Stanford University, Stanford, California 94305, United States
- SLAC National Accelerator Laboratory, 2575 Sand Hill Road, Menlo Park, California 94025, United States
| |
Collapse
|
10
|
Dey S, Folkestad SD, Paul AC, Koch H, Krylov AI. Core-ionization spectrum of liquid water. Phys Chem Chem Phys 2024; 26:1845-1859. [PMID: 38174659 DOI: 10.1039/d3cp02499g] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/05/2024]
Abstract
We present state-of-the-art calculations of the core-ionization spectrum of water. Despite significant progress in procedures developed to mitigate various experimental complications and uncertainties, the experimental determination of ionization energies of solvated species involves several non-trivial steps such as assessing the effect of the surface potential, electrolytes, and finite escape depths of photoelectrons. This provides a motivation to obtain robust theoretical values of the intrinsic bulk ionization energy and the corresponding solvent-induced shift. Here we develop theoretical protocols based on coupled-cluster theory and electrostatic embedding. Our value of the intrinsic solvent-induced shift of the 1sO ionization energy of water is -1.79 eV. The computed absolute position and the width of the 1sO peak in photoelectron spectrum of water are 538.47 eV and 1.44 eV, respectively, agreeing well with the best experimental values.
Collapse
Affiliation(s)
- Sourav Dey
- Department of Chemistry, University of Southern California, Los Angeles, CA 90089, USA.
| | - Sarai Dery Folkestad
- Department of Chemistry, University of Southern California, Los Angeles, CA 90089, USA.
| | - Alexander C Paul
- Department of Chemistry, Norwegian University of Science and Technology, 7491 Trondheim, Norway
| | - Henrik Koch
- Department of Chemistry, Norwegian University of Science and Technology, 7491 Trondheim, Norway
| | - Anna I Krylov
- Department of Chemistry, University of Southern California, Los Angeles, CA 90089, USA.
| |
Collapse
|
11
|
Folkestad SD, Paul AC, Ponzi A, Grazioli C, Coreno M, de Simone M, Koch H, Coriani S. Electronic Characterization of Glycolaldehyde: Experimental and Theoretical Insights from the Core- and Valence-Level Spectroscopy. J Phys Chem A 2023; 127:10621-10631. [PMID: 38084657 DOI: 10.1021/acs.jpca.3c06703] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/22/2023]
Abstract
The core-level electron excitation and ionization spectra of glycolaldehyde have been investigated by photoabsorption and photoemission spectroscopy at both carbon and oxygen K-edges; the valence ionization spectra were also recorded by photoelectron spectroscopy in the UV-vis region. The spectra are interpreted by means of ab initio calculations based on the equation-of-motion coupled cluster singles and doubles (EOM-CCSD) and coupled cluster singles, doubles, and perturbative are in good agreement with the experimental results, and many of the observed features are assigned. The photoabsorption spectra are not only dominated by transitions from core-level orbitals to unoccupied π and σ orbitals but also show structures due to Rydberg transitions.
Collapse
Affiliation(s)
- Sarai Dery Folkestad
- Department of Chemistry, NTNU─Norwegian University of Science and Technology, Trondheim 7491, Norway
| | - Alexander C Paul
- Department of Chemistry, NTNU─Norwegian University of Science and Technology, Trondheim 7491, Norway
| | - Aurora Ponzi
- Department of Physical Chemistry, R. Bošković Institute, Bijenička Cesta 54, Zagreb 10000, Croatia
| | - Cesare Grazioli
- IOM-CNR, Istituto Officina dei Materiali del Consiglio Nazionale delle Ricerche, Basovizza Area Science Park, Trieste I-34149, Italy
| | - Marcello Coreno
- ISM-CNR, Istituto di Struttura della Materia del Consiglio Nazionale delle Ricerche, Basovizza Area Science Park, Trieste I-34149, Italy
| | - Monica de Simone
- IOM-CNR, Istituto Officina dei Materiali del Consiglio Nazionale delle Ricerche, Basovizza Area Science Park, Trieste I-34149, Italy
| | - Henrik Koch
- Department of Chemistry, NTNU─Norwegian University of Science and Technology, Trondheim 7491, Norway
| | - Sonia Coriani
- DTU Chemistry, Technical University of Denmark, Kemitorvet Bldg 207, Kgs. Lyngby DK-2800, Denmark
| |
Collapse
|
12
|
Jacquemin D, Kossoski F, Gam F, Boggio-Pasqua M, Loos PF. Reference Vertical Excitation Energies for Transition Metal Compounds. J Chem Theory Comput 2023. [PMID: 37965941 DOI: 10.1021/acs.jctc.3c01080] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/16/2023]
Abstract
To enrich and enhance the diversity of the quest database of highly accurate excitation energies [Véril, M.; et al. Wiley Interdiscip. Rev.: Comput. Mol. Sci. 2021, 11, e1517], we report vertical transition energies in transition metal compounds. Eleven diatomic molecules with a singlet or doublet ground state containing a fourth-row transition metal (CuCl, CuF, CuH, ScF, ScH, ScO, ScS, TiN, ZnH, ZnO, and ZnS) are considered, and the corresponding excitation energies are computed using high-level coupled-cluster (CC) methods, namely, CC3, CCSDT, CC4, and CCSDTQ, as well as multiconfigurational methods such as CASPT2 and NEVPT2. In many cases, to provide more comprehensive benchmark data, we also provide full configuration interaction estimates computed with the configuration interaction using a perturbative selection made iteratively (CIPSI) method. Based on these calculations, theoretical best estimates of the transition energies are established in both the aug-cc-pVDZ and aug-cc-pVTZ basis sets. This allows us to accurately assess the performance of the CC and multiconfigurational methods for this specific set of challenging transitions. Furthermore, comparisons with experimental data and previous theoretical results are also reported.
Collapse
Affiliation(s)
- Denis Jacquemin
- Nantes Université, CNRS, CEISAM UMR 6230, F-44000 Nantes, France
- Institut Universitaire de France (IUF), F-75005 Paris, France
| | - Fábris Kossoski
- Laboratoire de Chimie et Physique Quantiques (UMR 5626), Université de Toulouse, CNRS, UPS, F-31062 Toulouse, France
| | - Franck Gam
- Nantes Université, CNRS, CEISAM UMR 6230, F-44000 Nantes, France
| | - Martial Boggio-Pasqua
- Laboratoire de Chimie et Physique Quantiques (UMR 5626), Université de Toulouse, CNRS, UPS, F-31062 Toulouse, France
| | - Pierre-François Loos
- Laboratoire de Chimie et Physique Quantiques (UMR 5626), Université de Toulouse, CNRS, UPS, F-31062 Toulouse, France
| |
Collapse
|
13
|
Dreuw A, Papapostolou A, Dempwolff AL. Algebraic Diagrammatic Construction Schemes Employing the Intermediate State Formalism: Theory, Capabilities, and Interpretation. J Phys Chem A 2023; 127:6635-6646. [PMID: 37498297 DOI: 10.1021/acs.jpca.3c02761] [Citation(s) in RCA: 4] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 07/28/2023]
Abstract
Algebraic diagrammatic construction (ADC) schemes represent a family of ab initio methods for the calculation of excited electronic states and electron-detached and -attached states. All ADC methods have been demonstrated to possess great potential for molecular applications, e.g., for the calculation of absorption or photoelectron spectra or electron attachment processes. ADC originates from Green's function or propagator theory; however, most recent ADC developments heavily rely on the intermediate state representation or effective Liouvillian formalisms, which comprise new ADC methods and computational schemes for high-order properties. The different approaches for the calculation of excitation energies, ionization potentials, and electron affinities are intimately related, and they provide a coherent description of these quantities at equivalent levels of theory and with comparable errors. Most quantum chemical program packages contain ADC methods; however, the most complete ADC suite of methods can be found in the recent release of Q-Chem.
Collapse
Affiliation(s)
- Andreas Dreuw
- Interdisciplinary Center for Scientific Computing, Ruprecht-Karls University, 69120 Heidelberg, Germany
| | - Antonia Papapostolou
- Interdisciplinary Center for Scientific Computing, Ruprecht-Karls University, 69120 Heidelberg, Germany
| | - Adrian L Dempwolff
- Interdisciplinary Center for Scientific Computing, Ruprecht-Karls University, 69120 Heidelberg, Germany
| |
Collapse
|
14
|
Huang M, Evangelista FA. A study of core-excited states of organic molecules computed with the generalized active space driven similarity renormalization group. J Chem Phys 2023; 158:124112. [PMID: 37003756 DOI: 10.1063/5.0137096] [Citation(s) in RCA: 3] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 04/03/2023] Open
Abstract
This work examines the accuracy and precision of x-ray absorption spectra computed with a multireference approach that combines generalized active space (GAS) references with the driven similarity renormalization group (DSRG). We employ the x-ray absorption benchmark of organic molecule (XABOOM) set, consisting of 116 transitions from mostly organic molecules [Fransson et al., J. Chem. Theory Comput. 17, 1618 (2021)]. Several approximations to a full-valence active space are examined and benchmarked. Absolute excitation energies and intensities computed with the GAS-DSRG truncated to second-order in perturbation theory are found to systematically underestimate experimental and reference theoretical values. Third-order perturbative corrections significantly improve the accuracy of GAS-DSRG absolute excitation energies, bringing the mean absolute deviation from experimental values down to 0.32 eV. The ozone molecule and glyoxylic acid are particularly challenging for second-order perturbation theory and are examined in detail to assess the importance of active space truncation and intruder states.
Collapse
Affiliation(s)
- Meng Huang
- Department of Chemistry and Cherry Emerson Center for Scientific Computation, Emory University, Atlanta, Georgia 30322, USA
| | - Francesco A Evangelista
- Department of Chemistry and Cherry Emerson Center for Scientific Computation, Emory University, Atlanta, Georgia 30322, USA
| |
Collapse
|
15
|
Schnack-Petersen AK, Moitra T, Folkestad SD, Coriani S. New Implementation of an Equation-of-Motion Coupled-Cluster Damped-Response Framework with Illustrative Applications to Resonant Inelastic X-ray Scattering. J Phys Chem A 2023; 127:1775-1793. [PMID: 36763003 DOI: 10.1021/acs.jpca.2c08181] [Citation(s) in RCA: 3] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/11/2023]
Abstract
We present an implementation of a damped response framework for calculating resonant inelastic X-ray scattering (RIXS) at the equation-of-motion coupled-cluster singles and doubles (CCSD) and second-order approximate coupled-cluster singles and doubles (CC2) levels of theory in the open-source program eT. This framework lays the foundation for future extension to higher excitation methods (notably, the coupled-cluster singles and doubles with perturbative triples, CC3) and to multilevel approaches. Our implementation adopts a fully relaxed ground state and different variants of the core-valence separation projection technique to address convergence issues. Illustrative results are compared with those obtained within the frozen-core core-valence separated approach, available in Q-Chem, as well as with experiment. The performance of the CC2 method is evaluated in comparison with that of CCSD. It is found that, while the CC2 method is noticeably inferior to CCSD for X-ray absorption spectra, the quality of the CC2 RIXS spectra is often comparable to that of the CCSD level of theory, when the same valence excited states are probed. Finally, we present preliminary RIXS results for a solvated molecule in aqueous solution.
Collapse
Affiliation(s)
| | - Torsha Moitra
- DTU Chemistry, Technical University of Denmark, DK-2800 Kongens Lyngby, Denmark.,Hylleraas Centre for Quantum Molecular Sciences, Department of Chemistry, UiTThe Arctic University of Norway, 9037 Tromsø, Norway
| | - Sarai Dery Folkestad
- Department of Chemistry, Norwegian University of Science and Technology, NO-7491 Trondheim, Norway
| | - Sonia Coriani
- DTU Chemistry, Technical University of Denmark, DK-2800 Kongens Lyngby, Denmark.,Department of Chemistry, Norwegian University of Science and Technology, NO-7491 Trondheim, Norway
| |
Collapse
|
16
|
Carlini L, Montorsi F, Wu Y, Bolognesi P, Borrego-Varillas R, Casavola AR, Castrovilli MC, Chiarinelli J, Mocci D, Vismarra F, Lucchini M, Nisoli M, Mukamel S, Garavelli M, Richter R, Nenov A, Avaldi L. Electron and ion spectroscopy of azobenzene in the valence and core shells. J Chem Phys 2023; 158:054201. [PMID: 36754795 DOI: 10.1063/5.0133824] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/19/2023] Open
Abstract
Azobenzene is a prototype and a building block of a class of molecules of extreme technological interest as molecular photo-switches. We present a joint experimental and theoretical study of its response to irradiation with light across the UV to x-ray spectrum. The study of valence and inner shell photo-ionization and excitation processes combined with measurement of valence photoelectron-photoion coincidence and mass spectra across the core thresholds provides a detailed insight into the site- and state-selected photo-induced processes. Photo-ionization and excitation measurements are interpreted via the multi-configurational restricted active space self-consistent field method corrected by second order perturbation theory. Using static modeling, we demonstrate that the carbon and nitrogen K edges of azobenzene are suitable candidates for exploring its photoinduced dynamics thanks to the transient signals appearing in background-free regions of the NEXAFS and XPS.
Collapse
Affiliation(s)
- L Carlini
- CNR-Istituto di Struttura Della Materia, CNR-ISM, Area Della Ricerca di Roma 1, Monterotondo, Italy
| | - F Montorsi
- Dipartimento di Chimica Industriale, Università Degli Studi di Bologna, Bologna, Italy
| | - Y Wu
- Dipartimento di Fisica, Politecnico di Milano, Piazza Leonardo da Vinci 32, Milano, Italy
| | - P Bolognesi
- CNR-Istituto di Struttura Della Materia, CNR-ISM, Area Della Ricerca di Roma 1, Monterotondo, Italy
| | - R Borrego-Varillas
- CNR-Istituto di Fotonica e Nanotecnologie, CNR-IFN, Piazza Leonardo da Vinci 32, Milano, Italy
| | - A R Casavola
- CNR-Istituto di Struttura Della Materia, CNR-ISM, Area Della Ricerca di Roma 1, Monterotondo, Italy
| | - M C Castrovilli
- CNR-Istituto di Struttura Della Materia, CNR-ISM, Area Della Ricerca di Roma 1, Monterotondo, Italy
| | - J Chiarinelli
- CNR-Istituto di Struttura Della Materia, CNR-ISM, Area Della Ricerca di Roma 1, Monterotondo, Italy
| | - D Mocci
- Dipartimento di Fisica, Politecnico di Milano, Piazza Leonardo da Vinci 32, Milano, Italy
| | - F Vismarra
- Dipartimento di Fisica, Politecnico di Milano, Piazza Leonardo da Vinci 32, Milano, Italy
| | - M Lucchini
- Dipartimento di Fisica, Politecnico di Milano, Piazza Leonardo da Vinci 32, Milano, Italy
| | - M Nisoli
- Dipartimento di Fisica, Politecnico di Milano, Piazza Leonardo da Vinci 32, Milano, Italy
| | - S Mukamel
- Department of Chemistry and Department of Physics and Astronomy, University of California, Irvine, California 92697, USA
| | - M Garavelli
- Dipartimento di Chimica Industriale, Università Degli Studi di Bologna, Bologna, Italy
| | - R Richter
- Elettra Sincrotrone Trieste, Area Science Park, Basovizza, Italy
| | - A Nenov
- Dipartimento di Chimica Industriale, Università Degli Studi di Bologna, Bologna, Italy
| | - L Avaldi
- CNR-Istituto di Struttura Della Materia, CNR-ISM, Area Della Ricerca di Roma 1, Monterotondo, Italy
| |
Collapse
|
17
|
Damour Y, Quintero-Monsebaiz R, Caffarel M, Jacquemin D, Kossoski F, Scemama A, Loos PF. Ground- and Excited-State Dipole Moments and Oscillator Strengths of Full Configuration Interaction Quality. J Chem Theory Comput 2023; 19:221-234. [PMID: 36548519 DOI: 10.1021/acs.jctc.2c01111] [Citation(s) in RCA: 6] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/24/2022]
Abstract
We report ground- and excited-state dipole moments and oscillator strengths (computed in different "gauges" or representations) of full configuration interaction (FCI) quality using the selected configuration interaction method known as Configuration Interaction using a Perturbative Selection made Iteratively (CIPSI). Thanks to a set encompassing 35 ground- and excited-state properties computed in 11 small molecules, the present near-FCI estimates allow us to assess the accuracy of high-order coupled-cluster (CC) calculations including up to quadruple excitations. In particular, we show that incrementing the excitation degree of the CC expansion (from CC with singles and doubles (CCSD) to CC with singles, doubles, and triples (CCSDT) or from CCSDT to CC with singles, doubles, triples, and quadruples (CCSDTQ)) reduces the average error with respect to the near-FCI reference values by approximately 1 order of magnitude.
Collapse
Affiliation(s)
- Yann Damour
- Laboratoire de Chimie et Physique Quantiques (UMR 5626), Université de Toulouse, CNRS, UPS, 31062 Toulouse, France
| | - Raúl Quintero-Monsebaiz
- Laboratoire de Chimie et Physique Quantiques (UMR 5626), Université de Toulouse, CNRS, UPS, 31062 Toulouse, France
| | - Michel Caffarel
- Laboratoire de Chimie et Physique Quantiques (UMR 5626), Université de Toulouse, CNRS, UPS, 31062 Toulouse, France
| | - Denis Jacquemin
- Nantes Université, CNRS, CEISAM UMR 6230, F-44000 Nantes, France.,Institut Universitaire de France (IUF), F-75005 Paris, France
| | - Fábris Kossoski
- Laboratoire de Chimie et Physique Quantiques (UMR 5626), Université de Toulouse, CNRS, UPS, 31062 Toulouse, France
| | - Anthony Scemama
- Laboratoire de Chimie et Physique Quantiques (UMR 5626), Université de Toulouse, CNRS, UPS, 31062 Toulouse, France
| | - Pierre-François Loos
- Laboratoire de Chimie et Physique Quantiques (UMR 5626), Université de Toulouse, CNRS, UPS, 31062 Toulouse, France
| |
Collapse
|
18
|
Glebov IO, Poddubnyy VV, Khokhlov DV. Perturbative Expansion of Nonorthogonal Product Approach for Charge Transfer States. J Phys Chem A 2022; 126:5800-5813. [PMID: 35998639 DOI: 10.1021/acs.jpca.2c02612] [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
Modeling of the excited states of multichromophoric systems is crucial for the understanding of photosynthesis functioning. The excitonic Hamiltonian method is widely used for such calculations. Excited states of the combined system are constructed from the wave functions of individual chromophores while their interactions are described by excitonic couplings. In the current study, we enhance a previously proposed nonorthogonal product approach to incorporate dynamic correlation effects accounted for by the multireference perturbation theory. We discuss the problems of constructing the excitonic Hamiltonian including charge transfer states for the molecular systems where the overlap contribution to the excitonic couplings is non-negligible. The benchmark calculations were performed for a model system. It was shown that the overlap component of the excitonic coupling is of great importance. The enhanced method provides an accurate description of the excited state energies and other properties.
Collapse
Affiliation(s)
- Ilya O Glebov
- Chemistry Department, Moscow State University, Leninskie gory 1-3, Moscow 119991, Russia
| | - Vladimir V Poddubnyy
- Chemistry Department, Moscow State University, Leninskie gory 1-3, Moscow 119991, Russia
| | - Daniil V Khokhlov
- Chemistry Department, Moscow State University, Leninskie gory 1-3, Moscow 119991, Russia
| |
Collapse
|
19
|
Paul AC, Folkestad SD, Myhre RH, Koch H. Oscillator Strengths in the Framework of Equation of Motion Multilevel CC3. J Chem Theory Comput 2022; 18:5246-5258. [PMID: 35921447 PMCID: PMC9476665 DOI: 10.1021/acs.jctc.2c00164] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
![]()
We present an efficient implementation of the equation
of motion
oscillator strengths for the closed-shell multilevel coupled cluster
singles and doubles with perturbative triples method (MLCC3) in the
electronic structure program eT. The orbital space is split into an active part treated with
CC3 and an inactive part computed at the coupled cluster singles and
doubles (CCSD) level of theory. Asymptotically, the CC3 contribution
scales as floating-point operations, where nV is the total number of virtual orbitals while nv and no are the
number of active virtual and occupied orbitals, respectively. The
CC3 contribution, thus, only scales linearly with the full system
size and can become negligible compared to the cost of CCSD. We demonstrate
the capabilities of our implementation by calculating the ultraviolet–visible
spectrum of azobenzene and a core excited state of betaine 30 with
more than 1000 molecular orbitals.
Collapse
Affiliation(s)
- Alexander C Paul
- Department of Chemistry, Norwegian University of Science and Technology, NTNU, 7491 Trondheim, Norway
| | - Sarai Dery Folkestad
- 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
|
20
|
Grabarz AM, Ośmiałowski B. Benchmarking Density Functional Approximations for Excited-State Properties of Fluorescent Dyes. Molecules 2021; 26:7434. [PMID: 34946515 PMCID: PMC8703901 DOI: 10.3390/molecules26247434] [Citation(s) in RCA: 10] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/25/2021] [Revised: 11/28/2021] [Accepted: 12/04/2021] [Indexed: 12/05/2022] Open
Abstract
This study presents an extensive analysis of the predictive power of time-dependent density functional theory in determining the excited-state properties of two groups of important fluorescent dyes, difluoroboranes and hydroxyphenylimidazo[1,2-a]pyridine derivatives. To ensure statistically meaningful results, the data set is comprised of 85 molecules manifesting diverse photophysical properties. The vertical excitation energies and dipole moments (in the electronic ground and excited states) of the aforementioned dyes were determined using the RI-CC2 method (reference) and with 18 density functional approximations (DFA). The set encompasses DFAs with varying amounts of exact exchange energy (EEX): from 0% (e.g., SVWN, BLYP), through a medium (e.g., TPSSh, B3LYP), up to a major contribution of EEX (e.g., BMK, MN15). It also includes range-separated hybrids (CAM-B3LYP, LC-BLYP). Similar error profiles of vertical energy were obtained for both dye groups, although the errors related to hydroxyphenylimidazopiridines are significantly larger. Overall, functionals including 40-55% of EEX (SOGGA11-X, BMK, M06-2X) ensure satisfactory agreement with the reference vertical excitation energies obtained using the RI-CC2 method; however, MN15 significantly outperforms them, providing a mean absolute error of merely 0.04 eV together with a very high correlation coefficient (R2 = 0.98). Within the investigated set of functionals, there is no single functional that would equally accurately determine ground- and excited-state dipole moments of difluoroboranes and hydroxyphenylimidazopiridine derivatives. Depending on the chosen set of dyes, the most accurate μGS predictions were delivered by MN15 incorporating a major EEX contribution (difluoroboranes) and by PBE0 containing a minor EEX fraction (hydroxyphenylimidazopiridines). Reverse trends are observed for μES, i.e., for difluoroboranes the best results were obtained with functionals including a minor fraction of EEX, specifically PBE0, while in the case of hydroxyphenylimidazopiridines, much more accurate predictions were provided by functionals incorporating a major EEX contribution (BMK, MN15).
Collapse
Affiliation(s)
- Anna M. Grabarz
- Faculty of Chemistry, Wrocław University of Science and Technology, Wyb. Wyspiańskiego 27, PL-50370 Wrocław, Poland
| | - Borys Ośmiałowski
- Faculty of Chemistry, Nicolaus Copernicus University, Gagarina 7, PL-87100 Toruń, Poland;
| |
Collapse
|
21
|
Hohenstein EG, Yu JK, Bannwarth C, List NH, Paul AC, Folkestad SD, Koch H, Martínez TJ. Predictions of Pre-edge Features in Time-Resolved Near-Edge X-ray Absorption Fine Structure Spectroscopy from Hole-Hole Tamm-Dancoff-Approximated Density Functional Theory. J Chem Theory Comput 2021; 17:7120-7133. [PMID: 34623139 DOI: 10.1021/acs.jctc.1c00478] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/16/2022]
Abstract
Time-resolved near-edge X-ray absorption fine structure (TR-NEXAFS) spectroscopy is a powerful technique for studying photochemical reaction dynamics with femtosecond time resolution. In order to avoid ambiguity in TR-NEXAFS spectra from nonadiabatic dynamics simulations, core- and valence-excited states must be evaluated on equal footing and those valence states must also define the potential energy surfaces used in the nonadiabatic dynamics simulation. In this work, we demonstrate that hole-hole Tamm-Dancoff-approximated density functional theory (hh-TDA) is capable of directly simulating TR-NEXAFS spectroscopies. We apply hh-TDA to the excited-state dynamics of acrolein. We identify two pre-edge features in the oxygen K-edge TR-NEXAFS spectrum associated with the S2 (ππ*) and S1 (nπ*) excited states. We show that these features can be used to follow the internal conversion dynamics between the lowest three electronic states of acrolein. Due to the low, O(N2) apparent computational complexity of hh-TDA and our GPU-accelerated implementation, this method is promising for the simulation of pre-edge features in TR-NEXAFS spectra of large molecules and molecules in the condensed phase.
Collapse
Affiliation(s)
- Edward G Hohenstein
- Department of Chemistry and The PULSE Institute, Stanford University, Stanford, California 94305, United States.,SLAC National Accelerator Laboratory, 2575 Sand Hill Road, Menlo Park, California 94025, United States
| | - Jimmy K Yu
- Department of Chemistry and The PULSE Institute, Stanford University, Stanford, California 94305, United States.,SLAC National Accelerator Laboratory, 2575 Sand Hill Road, Menlo Park, California 94025, United States.,Biophysics Program, Stanford University, Stanford, California 94305, United States
| | - Christoph Bannwarth
- Department of Chemistry and The PULSE Institute, Stanford University, Stanford, California 94305, United States.,SLAC National Accelerator Laboratory, 2575 Sand Hill Road, Menlo Park, California 94025, United States
| | - Nanna Holmgaard List
- Department of Chemistry and The PULSE Institute, Stanford University, Stanford, California 94305, United States.,SLAC National Accelerator Laboratory, 2575 Sand Hill Road, Menlo Park, California 94025, United States
| | - Alexander C Paul
- Department of Chemistry, Norwegian University of Science and Technology, NTNU, 7491 Trondheim, Norway
| | - Sarai D Folkestad
- 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
| | - Todd J Martínez
- Department of Chemistry and The PULSE Institute, Stanford University, Stanford, California 94305, United States.,SLAC National Accelerator Laboratory, 2575 Sand Hill Road, Menlo Park, California 94025, United States
| |
Collapse
|
22
|
Halbert L, Vidal ML, Shee A, Coriani S, Severo Pereira Gomes A. Relativistic EOM-CCSD for Core-Excited and Core-Ionized State Energies Based on the Four-Component Dirac-Coulomb(-Gaunt) Hamiltonian. J Chem Theory Comput 2021; 17:3583-3598. [PMID: 33944570 DOI: 10.1021/acs.jctc.0c01203] [Citation(s) in RCA: 20] [Impact Index Per Article: 6.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/06/2023]
Abstract
We report an implementation of the core-valence separation approach to the four-component relativistic Hamiltonian-based equation-of-motion coupled-cluster with singles and doubles theory (CVS-EOM-CCSD) for the calculation of relativistic core-ionization potentials and core-excitation energies. With this implementation, which is capable of exploiting double group symmetry, we investigate the effects of the different CVS-EOM-CCSD variants and the use of different Hamiltonians based on the exact two-component (X2C) framework on the energies of different core-ionized and -excited states in halogen- (CH3I, HX, and X-, X = Cl-At) and xenon-containing (Xe, XeF2) species. Our results show that the X2C molecular mean-field approach [Sikkema, J.; J. Chem. Phys. 2009, 131, 124116], based on four-component Dirac-Coulomb mean-field calculations (2DCM), is capable of providing core excitations and ionization energies that are nearly indistinguishable from the reference four-component energies for up to and including fifth-row elements. We observe that two-electron integrals over the small-component basis sets lead to non-negligible contributions to core binding energies for the K and L edges for atoms such as iodine or astatine and that the approach based on Dirac-Coulomb-Gaunt mean-field calculations (2DCGM) are significantly more accurate than X2C calculations for which screened two-electron spin-orbit interactions are included via atomic mean-field integrals.
Collapse
Affiliation(s)
- Loïc Halbert
- CNRS, UMR 8523-PhLAM-Physique des Lasers, Atomes et Molécules, Université de Lille, F-59000 Lille, France
| | - Marta L Vidal
- DTU Chemistry-Department of Chemistry, Technical University of Denmark, DK-2800 Kongens Lyngby, Denmark
| | - Avijit Shee
- Department of Chemistry, University of Michigan, Ann Arbor, Michigan 48109, United States
| | - Sonia Coriani
- DTU Chemistry-Department of Chemistry, Technical University of Denmark, DK-2800 Kongens Lyngby, Denmark
| | - André Severo Pereira Gomes
- CNRS, UMR 8523-PhLAM-Physique des Lasers, Atomes et Molécules, Université de Lille, F-59000 Lille, France
| |
Collapse
|
23
|
Ambroise MA, Dreuw A, Jensen F. Probing Basis Set Requirements for Calculating Core Ionization and Core Excitation Spectra Using Correlated Wave Function Methods. J Chem Theory Comput 2021; 17:2832-2842. [PMID: 33900755 DOI: 10.1021/acs.jctc.1c00042] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/10/2023]
Abstract
We investigate the basis set requirements for the accurate calculation of core excitations and core ionizations using correlated wave functions of coupled cluster type and linear response methods for describing the excitation. When a core excitation is described as an energy difference calculated using density functional theory, the basis set can be tailored to provide a balanced description of the reference- and excited-hole states. When the core excitation process is described by coupled cluster linear response methods, however, the basis set requirements are somewhat different. A systematic study of the sensitivity of the result to the basis set parameters suggests that a relatively large set of s- and p-type basis functions in combination with a careful selection of valence and core polarization functions is required. Based on these results, we propose a hierarchical sequence of basis sets, denoted ccX-nZ (n = D, T, Q, 5) for the atoms B-Ne, which are suitable for the calculation of core excitations by the correlated wave function linear response and equation-of-motion methods. The ccX-nZ series provides lower basis set errors for a given cardinal number or number of basis functions than other existing basis sets. For large systems, the ccX-nZ basis sets can be combined with the standard basis sets by placing the ccX-nZ only on the atoms where core excitations are of interest, but the accuracy of such mixed basis sets appears to be system-dependent.
Collapse
Affiliation(s)
- Maximilien A Ambroise
- Interdisciplinary Center for Scientific Computing, University of Heidelberg, 69117 Heidelberg, Germany
| | - Andreas Dreuw
- Interdisciplinary Center for Scientific Computing, University of Heidelberg, 69117 Heidelberg, Germany
| | - Frank Jensen
- Department of Chemistry, Aarhus University, DK-8000 Aarhus, Denmark
| |
Collapse
|
24
|
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
|
25
|
Gulania S, Kjønstad EF, Stanton JF, Koch H, Krylov AI. Equation-of-motion coupled-cluster method with double electron-attaching operators: Theory, implementation, and benchmarks. J Chem Phys 2021; 154:114115. [PMID: 33752380 DOI: 10.1063/5.0041822] [Citation(s) in RCA: 10] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
We report a production-level implementation of the equation-of-motion (EOM) coupled-cluster (CC) method with double electron-attaching (DEA) EOM operators of 2p and 3p1h types, EOM-DEA-CCSD. This ansatz, suitable for treating electronic structure patterns that can be described as two-electrons-in-many orbitals, represents a useful addition to the EOM-CC family of methods. We analyze the performance of EOM-DEA-CCSD for energy differences and molecular properties. By considering reduced quantities, such as state and transition one-particle density matrices, we compare EOM-DEA-CCSD wave functions with wave functions computed by other EOM-CCSD methods. The benchmarks illustrate that EOM-DEA-CCSD is capable of treating diradicals, bond-breaking, and some types of conical intersections.
Collapse
Affiliation(s)
- Sahil Gulania
- Department of Chemistry, University of Southern California, Los Angeles, California 90089, USA
| | - Eirik F Kjønstad
- Department of Chemistry, Norwegian University of Science and Technology, 7491 Trondheim, Norway
| | - John F Stanton
- Quantum Theory Project, Departments of Chemistry and Physics, University of Florida, Gainesville, Florida 32611, USA
| | - Henrik Koch
- Scuola Normale Superiore, Piazza dei Cavaleri 7, 56126 Pisa, Italy
| | - Anna I Krylov
- Department of Chemistry, University of Southern California, Los Angeles, California 90089, USA
| |
Collapse
|
26
|
Bilalbegović G, Maksimović A, Valencic LA, Lehtola S. Sulfur Molecules in Space by X-rays: A Computational Study. ACS EARTH & SPACE CHEMISTRY 2021; 5:436-448. [PMID: 33842801 PMCID: PMC8028330 DOI: 10.1021/acsearthspacechem.0c00238] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 08/29/2020] [Revised: 02/12/2021] [Accepted: 02/15/2021] [Indexed: 06/12/2023]
Abstract
X-ray astronomy lacks high resolution spectra of interstellar dust analogues and molecules, severely hampering interstellar medium studies based on upcoming X-ray missions. Various theoretical approaches may be used to address this problem, but they must first be shown to reproduce reliable spectra compared to the experiment. In this work, we calculate the sulfur K edge X-ray absorption spectra of H2S, SO2, and OCS, whose spectra are already known from X-ray experiments and predict the X-ray spectrum of CS, which as far as we are aware has not been measured, thereby hampering its detection by X-ray telescopes. We chose these four molecules as the astrochemistry of sulfur is an unsolved problem and as the four molecules are already known to exist in space. We consider three types of methods for modeling the X-ray spectra: more accurate calculations with the algebraic-diagrammatic construction (ADC) and the CC2, CCSD, and CC3 coupled cluster (CC) approaches as well as more affordable ones with transition potential density functional theory (TP-DFT). A comparison of our computational results to previously reported experimental spectra shows that the core-valence separation (CVS) approaches CVS-ADC(2)-x and CVS-CC3 generally yield a good qualitative level of agreement with the experiment, suggesting that they can be used for interpreting measured spectra, while the TP-DFT method is not reliable for these molecules. However, quantitative agreement with the experiment is still outside the reach of the computational methods studied in this work.
Collapse
Affiliation(s)
- Goranka Bilalbegović
- Department
of Physics, Faculty of Science, University
of Zagreb, Bijenička
cesta 32, 10000 Zagreb, Croatia
| | - Aleksandar Maksimović
- Center
of Excellence for Advanced Materials and Sensing Devices, Rudjer Bošković Institute, Bijenička cesta 54, 10000 Zagreb, Croatia
| | - Lynne A. Valencic
- NASA
Goddard Space Flight Center, Greenbelt, 20771 Maryland, United States
- Department
of Physics & Astronomy, The Johns Hopkins
University, 366 Bloomberg
Center, 3400 N. Charles Street, Baltimore, 21218 Maryland, United States
| | - Susi Lehtola
- Department
of Chemistry, University of Helsinki, P.O. Box 55, A. I. Virtasen aukio
1, FI-00014 Helsinki, Finland
| |
Collapse
|
27
|
Fransson T, Brumboiu IE, Vidal ML, Norman P, Coriani S, Dreuw A. XABOOM: An X-ray Absorption Benchmark of Organic Molecules Based on Carbon, Nitrogen, and Oxygen 1s → π* Transitions. J Chem Theory Comput 2021; 17:1618-1637. [PMID: 33544612 PMCID: PMC8023667 DOI: 10.1021/acs.jctc.0c01082] [Citation(s) in RCA: 31] [Impact Index Per Article: 10.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/14/2020] [Indexed: 01/05/2023]
Abstract
The performance of several standard and popular approaches for calculating X-ray absorption spectra at the carbon, nitrogen, and oxygen K-edges of 40 primarily organic molecules up to the size of guanine has been evaluated, focusing on the low-energy and intense 1s → π* transitions. Using results obtained with CVS-ADC(2)-x and fc-CVS-EOM-CCSD as benchmark references, we investigate the performance of CC2, ADC(2), ADC(3/2), and commonly adopted density functional theory (DFT)-based approaches. Here, focus is on precision rather than on accuracy of transition energies and intensities-in other words, we target relative energies and intensities and the spread thereof, rather than absolute values. The use of exchange-correlation functionals tailored for time-dependent DFT calculations of core excitations leads to error spreads similar to those seen for more standard functionals, despite yielding superior absolute energies. Long-range corrected functionals are shown to perform particularly well compared to our reference data, showing error spreads in energy and intensity of 0.2-0.3 eV and ∼10%, respectively, as compared to 0.3-0.6 eV and ∼20% for a typical pure hybrid. In comparing intensities, state mixing can complicate matters, and techniques to avoid this issue are discussed. Furthermore, the influence of basis sets in high-level ab initio calculations is investigated, showing that reasonably accurate results are obtained with the use of 6-311++G**. We name this benchmark suite as XABOOM (X-ray absorption benchmark of organic molecules) and provide molecular structures and ground-state self-consistent field energies and spectroscopic data. We believe that it provides a good assessment of electronic structure theory methods for calculating X-ray absorption spectra and will become useful for future developments in this field.
Collapse
Affiliation(s)
- Thomas Fransson
- Interdisciplinary
Center for Scientific Computing, Ruprecht-Karls
University, Im Neuenheimer
Feld 205, 69120 Heidelberg, Germany
- Fysikum, Stockholm University, Albanova, 10691 Stockholm, Sweden
| | - Iulia E. Brumboiu
- Department
of Theoretical Chemistry and Biology, KTH
Royal Institute of Technology, 10691 Stockholm, Sweden
- Department
of Chemistry, Korea Advanced Institute of
Science and Technology, 34141 Daejeon, Korea
| | - Marta L. Vidal
- DTU
Chemistry, Technical University of Denmark, Kemitorvet Bldg 207, DK-2800 Kongens Lyngby, Denmark
| | - Patrick Norman
- Department
of Theoretical Chemistry and Biology, KTH
Royal Institute of Technology, 10691 Stockholm, Sweden
| | - Sonia Coriani
- DTU
Chemistry, Technical University of Denmark, Kemitorvet Bldg 207, DK-2800 Kongens Lyngby, Denmark
- Department
of Chemistry, NTNU-Norwegian University
of Science and Technology, N-7991 Trondheim, Norway
| | - Andreas Dreuw
- Interdisciplinary
Center for Scientific Computing, Ruprecht-Karls
University, Im Neuenheimer
Feld 205, 69120 Heidelberg, Germany
| |
Collapse
|
28
|
Tsuru S, Vidal ML, Pápai M, Krylov AI, Møller KB, Coriani S. An assessment of different electronic structure approaches for modeling time-resolved x-ray absorption spectroscopy. STRUCTURAL DYNAMICS (MELVILLE, N.Y.) 2021; 8:024101. [PMID: 33786337 PMCID: PMC7986275 DOI: 10.1063/4.0000070] [Citation(s) in RCA: 16] [Impact Index Per Article: 5.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 12/14/2020] [Accepted: 02/11/2021] [Indexed: 05/06/2023]
Abstract
We assess the performance of different protocols for simulating excited-state x-ray absorption spectra. We consider three different protocols based on equation-of-motion coupled-cluster singles and doubles, two of them combined with the maximum overlap method. The three protocols differ in the choice of a reference configuration used to compute target states. Maximum-overlap-method time-dependent density functional theory is also considered. The performance of the different approaches is illustrated using uracil, thymine, and acetylacetone as benchmark systems. The results provide guidance for selecting an electronic structure method for modeling time-resolved x-ray absorption spectroscopy.
Collapse
Affiliation(s)
- Shota Tsuru
- DTU Chemistry, Technical University of Denmark, Kemitorvet Building 207, DK-2800 Kgs. Lyngby, Denmark
| | - Marta L. Vidal
- DTU Chemistry, Technical University of Denmark, Kemitorvet Building 207, DK-2800 Kgs. Lyngby, Denmark
| | - Mátyás Pápai
- DTU Chemistry, Technical University of Denmark, Kemitorvet Building 207, DK-2800 Kgs. Lyngby, Denmark
| | - Anna I. Krylov
- Department of Chemistry, University of Southern California, Los Angeles, California 90089, USA
| | - Klaus B. Møller
- DTU Chemistry, Technical University of Denmark, Kemitorvet Building 207, DK-2800 Kgs. Lyngby, Denmark
| | - Sonia Coriani
- DTU Chemistry, Technical University of Denmark, Kemitorvet Building 207, DK-2800 Kgs. Lyngby, Denmark
| |
Collapse
|
29
|
Fedotov DA, Paul AC, Posocco P, Santoro F, Garavelli M, Koch H, Coriani S, Improta R. Excited-State Absorption of Uracil in the Gas Phase: Mapping the Main Decay Paths by Different Electronic Structure Methods. J Chem Theory Comput 2021; 17:1638-1652. [PMID: 33529532 DOI: 10.1021/acs.jctc.0c01150] [Citation(s) in RCA: 9] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/21/2022]
Abstract
We present a computational study of the one-photon and excited-state absorption (ESA) from the two lowest energy excited states of uracil in the gas phase: an nπ* dark state (1n) and the lowest energy bright ππ* state (1π). The predictions of six different linear response electronic structure methods, namely, TD-CAM-B3LYP, EOM-CCSD, EOM-CC3, ADC(2), ADC(2)-x, and ADC(3) are critically compared. In general, the spectral shapes predicted by TD-CAM-B3LYP, EOM-CCSD, EOM-CC3, and ADC(3) are fairly similar, though the quality of TD-CAM-B3LYP slightly deteriorates in the high-energy region. By computing the spectra at some key structures on different potential energy surfaces (PES), that is, the Franck-Condon point, the 1n minimum, and structures representative of different regions of the 1π PES, we obtain important insights into the shift of the ESA spectra, following the motion of the wavepacket on the excited-state PES. Though 1π has larger ESA than 1n, some spectral regions are dominated by these latter signals. Aside from its methodological interest, we thus obtain interesting indications to interpret transient absorption spectra to disentangle the photoactivated dynamics of nucleobases.
Collapse
Affiliation(s)
- Daniil A Fedotov
- DTU Chemistry, Technical University of Denmark, Kemitorvet Bldg 207, DK-2800 Kongens Lyngby, Denmark
| | - Alexander C Paul
- Department of Chemistry, NTNU-Norwegian University of Science and Technology, N-7491 Trondheim, Norway
| | - Paolo Posocco
- Dipartimento di Scienze Chimiche e Farmaceutiche, Università degli Studi di Trieste, via L. Giorgieri 1, I-34127 Trieste, Italy
| | - Fabrizio Santoro
- Istituto di Chimica dei Composti Organometallici (ICCOM-CNR), Area della Ricerca del CNR, Via Moruzzi 1, I-56124 Pisa, Italy
| | - Marco Garavelli
- Department of Industrial Chemistry "Toso Montanari", Università degli Studi di Bologna, I-40126 Bologna, Italy
| | - Henrik Koch
- Department of Chemistry, NTNU-Norwegian University of Science and Technology, N-7491 Trondheim, Norway.,Scuola Normale Superiore, Piazza dei Cavalieri 7, I-56100 Pisa, Italy
| | - Sonia Coriani
- DTU Chemistry, Technical University of Denmark, Kemitorvet Bldg 207, DK-2800 Kongens Lyngby, Denmark.,Department of Chemistry, NTNU-Norwegian University of Science and Technology, N-7491 Trondheim, Norway
| | - Roberto Improta
- Istituto di Biostrutture e Bioimmagini-CNR, Via Mezzocannone 6, I-80134 Napoli, Italy
| |
Collapse
|
30
|
Folkestad SD, Koch H. Equation-of-Motion MLCCSD and CCSD-in-HF Oscillator Strengths and Their Application to Core Excitations. J Chem Theory Comput 2020; 16:6869-6879. [PMID: 32955866 PMCID: PMC8011930 DOI: 10.1021/acs.jctc.0c00707] [Citation(s) in RCA: 15] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/31/2023]
Abstract
We present an implementation of equation-of-motion oscillator strengths for the multilevel CCSD (MLCCSD) model where CCS is used as the lower level method (CCS/CCSD). In this model, the double excitations of the cluster operator are restricted to an active orbital space, whereas the single excitations are unrestricted. Calculated nitrogen K-edge spectra of adenosine, adenosine triphosphate (ATP), and an ATP-water system are used to demonstrate the performance of the model. Projected atomic orbitals (PAOs) are used to partition the virtual space into active and inactive orbital sets. Cholesky decomposition of the Hartree-Fock density is used to partition the occupied orbitals. This Cholesky-PAO partitioning is cheap, scaling as O(N3), and is suitable for the calculation of core excitations, which are localized in character. By restricting the single excitations of the cluster operator to the active space, as well as the double excitations, the CCSD-in-HF model is obtained. A comparison of the two models-MLCCSD and CCSD-in-HF-is presented for the core excitation spectra of the adenosine and ATP systems.
Collapse
Affiliation(s)
- Sarai Dery Folkestad
- Department of Chemistry, Norwegian University of Science and Technology, Trondheim N-7491, Norway
| | - Henrik Koch
- Department of Chemistry, Norwegian University of Science and Technology, Trondheim N-7491, Norway.,Scuola Normale Superiore, Piazza dei Cavaleri 7, Pisa 56126, Italy
| |
Collapse
|