1
|
Manna A, Jangid B, Pant R, Dutta AK. Efficient State-Specific Natural Orbital Based Equation of Motion Coupled Cluster Method for Core-Ionization Energies: Theory, Implementation, and Benchmark. J Chem Theory Comput 2024. [PMID: 39073757 DOI: 10.1021/acs.jctc.4c00546] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 07/30/2024]
Abstract
We have implemented a reduced-cost partial triples correction scheme to the equation of motion coupled cluster method for core-ionization energy based on state-specific natural orbitals. The second-order Algebraic Diagrammatic Construction (ADC) method is used to generate the state-specific natural orbital, which provides quicker convergence of the core-IP value with respect to the size of the virtual space than that observed in standard MP2-based natural orbitals. The error due to truncation of the virtual orbital can be reduced by using a perturbative correction. The accuracy of the method can be controlled by a single threshold, and there is a black box to use. The inclusion of the partial triples correction in the natural orbital based EOM-CCSD method greatly improves the agreement of the results with the experiment. The efficiency of the present implementation is demonstrated by calculating the core-ionization energy of a molecule containing 60 atoms and more than 2000 basis functions.
Collapse
Affiliation(s)
- Amrita Manna
- Department of Chemistry, Indian Institute of Technology Bombay, Powai, Mumbai 400076, India
| | - Bhavnesh Jangid
- Department of Chemistry, Indian Institute of Technology Bombay, Powai, Mumbai 400076, India
| | - Rakesh Pant
- Department of Chemistry, Indian Institute of Technology Bombay, Powai, Mumbai 400076, India
| | - Achintya Kumar Dutta
- Department of Chemistry, Indian Institute of Technology Bombay, Powai, Mumbai 400076, India
| |
Collapse
|
2
|
Majumder R, Sokolov AY. Consistent Second-Order Treatment of Spin-Orbit Coupling and Dynamic Correlation in Quasidegenerate N-Electron Valence Perturbation Theory. J Chem Theory Comput 2024; 20:4676-4688. [PMID: 38795071 DOI: 10.1021/acs.jctc.4c00458] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/27/2024]
Abstract
We present a formulation and implementation of second-order quasidegenerate N-electron valence perturbation theory (QDNEVPT2) that provides a balanced and accurate description of spin-orbit coupling and dynamic correlation effects in multiconfigurational electronic states. In our approach, the energies and wave functions of electronic states are computed by treating electron repulsion and spin-orbit coupling operators as equal perturbations to the nonrelativistic complete active-space wave functions, and their contributions are incorporated fully up to the second order. The spin-orbit effects are described using the Breit-Pauli (BP) or exact two-component Douglas-Kroll-Hess (DKH) Hamiltonians within spin-orbit mean-field approximation. The resulting second-order methods (BP2- and DKH2-QDNEVPT2) are capable of treating spin-orbit coupling effects in nearly degenerate electronic states by diagonalizing an effective Hamiltonian expanded in a compact non-relativistic basis. For a variety of atoms and small molecules across the entire periodic table, we demonstrate that DKH2-QDNEVPT2 is competitive in accuracy with variational two-component relativistic theories. BP2-QDNEVPT2 shows high accuracy for the second- and third-period elements, but its performance deteriorates for heavier atoms and molecules. We also consider the first-order spin-orbit QDNEVPT2 approximations (BP1- and DKH1-QDNEVPT2), among which DKH1-QDNEVPT2 is reliable but less accurate than DKH2-QDNEVPT2. Both DKH1- and DKH2-QDNEVPT2 hold promise as efficient and accurate electronic structure methods for treating electron correlation and spin-orbit coupling in a variety of applications.
Collapse
Affiliation(s)
- Rajat Majumder
- Department of Chemistry and Biochemistry, The Ohio State University, Columbus, Ohio 43210, United States
| | - Alexander Yu Sokolov
- Department of Chemistry and Biochemistry, The Ohio State University, Columbus, Ohio 43210, United States
| |
Collapse
|
3
|
Jangid B, Hermes MR, Gagliardi L. Core Binding Energy Calculations: A Scalable Approach with the Quantum Embedding-Based Equation-of-Motion Coupled-Cluster Method. J Phys Chem Lett 2024; 15:5954-5963. [PMID: 38810243 DOI: 10.1021/acs.jpclett.4c00957] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/31/2024]
Abstract
We investigated the use of density matrix embedding theory to facilitate the computation of core ionization energies (IPs) of large molecules at the equation-of-motion coupled-cluster singles doubles with perturbative triples (EOM-CCSD*) level in combination with the core-valence separation (CVS) approximation. The unembedded IP-CVS-EOM-CCSD* method with a triple-ζ basis set produced ionization energies within 1 eV of experiment with a standard deviation of ∼0.2 eV for the core65 data set. The embedded variant contributed very little systematic error relative to the unembedded method, with a mean unsigned error of 0.07 eV and a standard deviation of ∼0.1 eV, in exchange for accelerating the calculations by many orders of magnitude. By employing embedded EOM-CC methods, we computed the core ionization energies of the uracil hexamer, doped fullerene, and chlorophyll molecule, utilizing up to ∼4000 basis functions within 1 eV from experimental values. Such calculations are not currently possible with the unembedded EOM-CC method.
Collapse
Affiliation(s)
- Bhavnesh Jangid
- Department of Chemistry, The University of Chicago, Chicago, Illinois 60637, United States
| | - Matthew R Hermes
- Department of Chemistry, The University of Chicago, Chicago, Illinois 60637, United States
| | - Laura Gagliardi
- Department of Chemistry, The University of Chicago, Chicago, Illinois 60637, United States
- Pritzker School of Molecular Engineering, The University of Chicago, Chicago, Illinois 60637, United States
| |
Collapse
|
4
|
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
|
5
|
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
|
6
|
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
|
7
|
Datar A, Gudivada S, Matthews DA. Ab Initio Investigation of Intramolecular Charge Transfer States in DMABN by Calculation of Excited State X-ray Absorption Spectra. J Phys Chem A 2023. [PMID: 37209154 DOI: 10.1021/acs.jpca.3c01409] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/22/2023]
Abstract
Dual fluorescence in 4-(dimethylamino)benzonitrile (DMABN) and its derivatives in polar solvents has been studied extensively for the past several decades. An intramolecular charge transfer (ICT) minimum on the excited state potential energy surface, in addition to the localized low-energy (LE) minimum, has been proposed as a mechanism for this dual fluorescence, with large geometric relaxation and molecular orbital reorganization a key feature of the ICT pathway. Herein, we have used both equation-of-motion coupled-cluster with single and double excitations (EOM-CCSD) and time-dependent density functional (TDDFT) methods to investigate the landscape of excited state potential energy surfaces across a number of geometric conformations proposed as ICT structures. In order to correlate these geometries and valence excited states in terms of potential experimental observables, we have calculated the nitrogen K-edge ground and excited state absorption spectra for each of the predicted "signpost" structures and identified several key spectral features that could be used to interpret a future time-resolved X-ray absorption experiment.
Collapse
Affiliation(s)
- Avdhoot Datar
- Department of Chemistry, Southern Methodist University, Dallas, Texas 75275, United States
| | - Saisrinivas Gudivada
- Department of Chemistry, Southern Methodist University, Dallas, Texas 75275, United States
| | - Devin A Matthews
- Department of Chemistry, Southern Methodist University, Dallas, Texas 75275, United States
| |
Collapse
|
8
|
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
|
9
|
Majumder R, Sokolov AY. Simulating Spin-Orbit Coupling with Quasidegenerate N-Electron Valence Perturbation Theory. J Phys Chem A 2023; 127:546-559. [PMID: 36599072 DOI: 10.1021/acs.jpca.2c07952] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/06/2023]
Abstract
We present the first implementation of spin-orbit coupling effects in fully internally contracted second-order quasidegenerate N-electron valence perturbation theory (SO-QDNEVPT2). The SO-QDNEVPT2 approach enables the computations of ground- and excited-state energies and oscillator strengths combining the description of static electron correlation with an efficient treatment of dynamic correlation and spin-orbit coupling. In addition to SO-QDNEVPT2 with the full description of one- and two-body spin-orbit interactions at the level of two-component Breit-Pauli Hamiltonian, our implementation also features a simplified approach that takes advantage of spin-orbit mean-field approximation (SOMF-QDNEVPT2). The accuracy of these methods is tested for the group 14 and 16 hydrides, 3d and 4d transition metal ions, and two actinide dioxides (neptunyl and plutonyl dications). The zero-field splittings of group 14 and 16 molecules computed using SO-QDNEVPT2 and SOMF-QDNEVPT2 are in good agreement with the available experimental data. For the 3d transition metal ions, the SO-QDNEVPT2 method is significantly more accurate than SOMF-QDNEVPT2, while no substantial difference in the performance of two methods is observed for the 4d ions. Finally, we demonstrate that for the actinide dioxides the results of SO-QDNEVPT2 and SOMF-QDNEVPT2 are in good agreement with the data from previous theoretical studies of these systems. Overall, our results demonstrate that SO-QDNEVPT2 and SOMF-QDNEVPT2 are promising multireference methods for treating spin-orbit coupling with a relatively low computational cost.
Collapse
Affiliation(s)
- Rajat Majumder
- Department of Chemistry and Biochemistry, The Ohio State University, Columbus, Ohio43210, United States
| | - Alexander Yu Sokolov
- Department of Chemistry and Biochemistry, The Ohio State University, Columbus, Ohio43210, United States
| |
Collapse
|
10
|
Simons M, Matthews DA. Transition-potential coupled cluster II: optimisation of the core orbital occupation number. Mol Phys 2022. [DOI: 10.1080/00268976.2022.2088421] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/18/2022]
|
11
|
Hait D, Oosterbaan KJ, Carter-Fenk K, Head-Gordon M. Computing x-ray absorption spectra from linear-response particles atop optimized holes. J Chem Phys 2022; 156:201104. [PMID: 35649868 DOI: 10.1063/5.0092987] [Citation(s) in RCA: 6] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/24/2023] Open
Abstract
State specific orbital optimized density functional theory (OO-DFT) methods, such as restricted open-shell Kohn-Sham (ROKS), can attain semiquantitative accuracy for predicting x-ray absorption spectra of closed-shell molecules. OO-DFT methods, however, require that each state be individually optimized. In this Communication, we present an approach to generate an approximate core-excited state density for use with the ROKS energy ansatz, which is capable of giving reasonable accuracy without requiring state-specific optimization. This is achieved by fully optimizing the core-hole through the core-ionized state, followed by the use of electron-addition configuration interaction singles to obtain the particle level. This hybrid approach can be viewed as a DFT generalization of the static-exchange (STEX) method and can attain ∼0.6 eV rms error for the K-edges of C-F through the use of local functionals, such as PBE and OLYP. This ROKS(STEX) approach can also be used to identify important transitions for full OO ROKS treatment and can thus help reduce the computational cost of obtaining OO-DFT quality spectra. ROKS(STEX), therefore, appears to be a useful technique for the efficient prediction of x-ray absorption spectra.
Collapse
Affiliation(s)
- Diptarka Hait
- Kenneth S. Pitzer Center for Theoretical Chemistry, Department of Chemistry, University of California, Berkeley, California 94720, USA
| | - Katherine J Oosterbaan
- Kenneth S. Pitzer Center for Theoretical Chemistry, Department of Chemistry, University of California, Berkeley, California 94720, USA
| | - Kevin Carter-Fenk
- Kenneth S. Pitzer Center for Theoretical Chemistry, Department of Chemistry, University of California, Berkeley, California 94720, USA
| | - Martin Head-Gordon
- Kenneth S. Pitzer Center for Theoretical Chemistry, Department of Chemistry, University of California, Berkeley, California 94720, USA
| |
Collapse
|
12
|
Simons M, Matthews DA. Accurate Core-Excited States via Inclusion of Core Triple Excitations in Similarity-Transformed Equation-of-Motion Theory. J Chem Theory Comput 2022; 18:3759-3765. [PMID: 35536592 DOI: 10.1021/acs.jctc.2c00268] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
The phenomenon of orbital relaxation upon excitation of core electrons is a major problem in the linear-response treatment of core-hole spectroscopies. Rather than addressing relaxation through direct dynamical correlation of the excited state via equation-of-motion coupled cluster theory (EOMEE-CC), we extend the alternative similarity-transformed equation-of-motion coupled cluster theory (STEOMEE-CC) by including the core-valence separation (CVS) and correlation of triple excitations only within the calculation of core ionization energies. This new method, CVS-STEOMEE-CCSD+cT, significantly improves on CVS-EOMEE-CCSD and unmodified CVS-STEOMEE-CCSD when compared to full CVS-EOM-CCSDT for K-edge core-excitation energies of a set of small molecules. The improvement in both absolute and relative (shifted) peak positions is nearly as good as that for transition-potential coupled cluster (TP-CC), which includes an explicit treatment of orbital relaxation, and CVS-EOMEE-CCSD*, which includes a perturbative treatment of triple excitations.
Collapse
Affiliation(s)
- Megan Simons
- Department of Chemistry, Southern Methodist University, Dallas, Texas 75275, United States
| | - Devin A Matthews
- Department of Chemistry, Southern Methodist University, Dallas, Texas 75275, United States
| |
Collapse
|
13
|
Cunha LA, Hait D, Kang R, Mao Y, Head-Gordon M. Relativistic Orbital-Optimized Density Functional Theory for Accurate Core-Level Spectroscopy. J Phys Chem Lett 2022; 13:3438-3449. [PMID: 35412838 DOI: 10.1021/acs.jpclett.2c00578] [Citation(s) in RCA: 20] [Impact Index Per Article: 10.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/14/2023]
Abstract
Core-level spectra of 1s electrons of elements heavier than Ne show significant relativistic effects. We combine advances in orbital-optimized density functional theory (OO-DFT) with the spin-free exact two-component (X2C) model for scalar relativistic effects to study K-edge spectra of third period elements. OO-DFT/X2C is found to be quite accurate at predicting energies, yielding a ∼0.5 eV root-mean-square error versus experiment with the modern SCAN (and related) functionals. This marks a significant improvement over the >50 eV deviations that are typical for the popular time-dependent DFT (TDDFT) approach. Consequently, experimental spectra are quite well reproduced by OO-DFT/X2C, sans empirical shifts for alignment. OO-DFT/X2C combines high accuracy with ground state DFT cost and is thus a promising route for computing core-level spectra of third period elements. We also explored K and L edges of 3d transition metals to identify limitations of the OO-DFT/X2C approach in modeling the spectra of heavier atoms.
Collapse
Affiliation(s)
- Leonardo A Cunha
- Kenneth S. Pitzer Center for Theoretical Chemistry, Department of Chemistry, University of California, Berkeley, California 94720, United States
- Chemical Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, United States
| | - Diptarka Hait
- Kenneth S. Pitzer Center for Theoretical Chemistry, Department of Chemistry, University of California, Berkeley, California 94720, United States
- Chemical Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, United States
| | - Richard Kang
- Kenneth S. Pitzer Center for Theoretical Chemistry, Department of Chemistry, University of California, Berkeley, California 94720, United States
- Chemical Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, United States
| | - Yuezhi Mao
- Department of Chemistry, Stanford University, Stanford, California 94305, United States
| | - Martin Head-Gordon
- Kenneth S. Pitzer Center for Theoretical Chemistry, Department of Chemistry, University of California, Berkeley, California 94720, United States
- Chemical Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, United States
| |
Collapse
|
14
|
Ranga S, Dutta AK. A Core-Valence Separated Similarity Transformed EOM-CCSD Method for Core-Excitation Spectra. J Chem Theory Comput 2021; 17:7428-7446. [PMID: 34814683 DOI: 10.1021/acs.jctc.1c00402] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
We present the theory and implementation of a core-valence separated similarity transformed EOM-CCSD (STEOM-CCSD) method for K-edge core excitation spectra. The method can select an appropriate active space using CIS natural orbitals and near "black box" to use. The second similarity transformed Hamiltonian is diagonalized in the space of single excitation. Therefore, the final diagonalization step is free from the convergence problem arising due to the coupling of the core-excited states with the continuum of doubly excited states. Convergence trouble can appear for the preceding core-ionized state calculation in STEOM-CCSD. A core-valence separation (CVS) scheme compatible with the natural orbital based active space selection (CVS-STEOM-CCSD-NO) is implemented to overcome the problem. The CVS-STEOM-CCSD-NO has a similar accuracy to that of the standard CVS-EOM-CCSD method but comes with a lower computational cost. The modification required in the CVS scheme to make use of the CIS natural orbital is highlighted. The suitability of the CVS-STEOM-CCSD-NO method for chemical application is demonstrated by simulating the K-edge spectra of glycine and thymine.
Collapse
Affiliation(s)
- Santosh Ranga
- Department of Chemistry, Indian Institute of Technology Bombay, Powai, Mumbai 400076, India
| | - Achintya Kumar Dutta
- Department of Chemistry, Indian Institute of Technology Bombay, Powai, Mumbai 400076, India
| |
Collapse
|
15
|
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
|
16
|
Moitra T, Coriani S, Cabral Tenorio BN. Inner-shell photoabsorption and photoionisation cross-sections of valence excited states from asymmetric-Lanczos equation-of-motion coupled cluster singles and doubles theory. Mol Phys 2021. [DOI: 10.1080/00268976.2021.1980235] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/13/2022]
Affiliation(s)
- Torsha Moitra
- DTU Chemistry–Department of Chemistry, Technical University of Denmark, Kongens Lyngby, Denmark
| | - Sonia Coriani
- DTU Chemistry–Department of Chemistry, Technical University of Denmark, Kongens Lyngby, Denmark
| | | |
Collapse
|
17
|
Reinholdt P, Vidal ML, Kongsted J, Iannuzzi M, Coriani S, Odelius M. Nitrogen K-Edge X-ray Absorption Spectra of Ammonium and Ammonia in Water Solution: Assessing the Performance of Polarizable Embedding Coupled Cluster Methods. J Phys Chem Lett 2021; 12:8865-8871. [PMID: 34498464 PMCID: PMC8450933 DOI: 10.1021/acs.jpclett.1c02031] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/10/2023]
Abstract
The recent development of liquid jet and liquid leaf sample delivery systems allows for accurate measurements of soft X-ray absorption spectra in transmission mode of solutes in a liquid environment. As this type of measurement becomes increasingly accessible, there is a strong need for reliable theoretical methods for assisting in the interpretation of the experimental data. Coupled cluster methods have been extensively developed over the past decade to simulate X-ray absorption in the gas phase. Their performance for solvated species, on the contrary, remains largely unexplored. Here, we investigate the current state of the art of coupled cluster modeling of nitrogen K-edge X-ray absorption of aqueous ammonia and ammonium based on quantum mechanics/molecular mechanics, where both the level of coupled cluster calculations and polarizable embedding are scrutinized. The results are compared to existing experimental data as well as simulations based on transition potential density functional theory.
Collapse
Affiliation(s)
- Peter Reinholdt
- Institut
for Fysik, Kemi og Farmaci, Syddansk Universitet, DK-5230 Odense, Denmark
| | - Marta L. Vidal
- DTU
Chemistry, Technical University of Denmark, DK-2800 Kongens
Lyngby, Denmark
| | - Jacob Kongsted
- Institut
for Fysik, Kemi og Farmaci, Syddansk Universitet, DK-5230 Odense, Denmark
| | - Marcella Iannuzzi
- Physical
Chemistry Institute, University of Zürich, 8057 Zürich, Switzerland
| | - Sonia Coriani
- DTU
Chemistry, Technical University of Denmark, DK-2800 Kongens
Lyngby, Denmark
| | - Michael Odelius
- Department
of Physics, AlbaNova University Center, Stockholm University, SE-106 91 Stockholm, Sweden
| |
Collapse
|
18
|
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: 29] [Impact Index Per Article: 9.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
|
19
|
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
|
20
|
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: 34] [Impact Index Per Article: 11.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
|
21
|
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
|
22
|
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
|
23
|
Hait D, Haugen EA, Yang Z, Oosterbaan KJ, Leone SR, Head-Gordon M. Accurate prediction of core-level spectra of radicals at density functional theory cost via square gradient minimization and recoupling of mixed configurations. J Chem Phys 2020; 153:134108. [DOI: 10.1063/5.0018833] [Citation(s) in RCA: 17] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Diptarka Hait
- Department of Chemistry, University of California, Berkeley, California 94720, USA
- Chemical Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, USA
| | - Eric A. Haugen
- Department of Chemistry, University of California, Berkeley, California 94720, USA
- Chemical Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, USA
| | - Zheyue Yang
- Department of Chemistry, University of California, Berkeley, California 94720, USA
- Chemical Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, USA
| | - Katherine J. Oosterbaan
- Department of Chemistry, University of California, Berkeley, California 94720, USA
- Chemical Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, USA
| | - Stephen R. Leone
- Department of Chemistry, University of California, Berkeley, California 94720, USA
- Chemical Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, USA
- Department of Physics, University of California, Berkeley, California 94720, USA
| | - Martin Head-Gordon
- Department of Chemistry, University of California, Berkeley, California 94720, USA
- Chemical Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, USA
| |
Collapse
|
24
|
Vidal ML, Pokhilko P, Krylov AI, Coriani S. Equation-of-Motion Coupled-Cluster Theory to Model L-Edge X-ray Absorption and Photoelectron Spectra. J Phys Chem Lett 2020; 11:8314-8321. [PMID: 32897075 DOI: 10.1021/acs.jpclett.0c02027] [Citation(s) in RCA: 48] [Impact Index Per Article: 12.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/03/2023]
Abstract
We present an extension of the equation-of-motion coupled-cluster singles and doubles (EOM-CCSD) theory for computing X-ray L-edge spectra, both in the absorption (XAS) and in the photoelectron (XPS) regimes. The approach is based on the perturbative evaluation of spin-orbit couplings using the Breit-Pauli Hamiltonian and nonrelativistic wave functions described by the fc-CVS-EOM-CCSD ansatz (EOM-CCSD within the frozen-core core-valence separated (fc-CVS) scheme). The formalism is based on spinless one-particle density matrices. The approach is illustrated by modeling XAS and XPS of several model systems ranging from Ar to small molecules containing sulfur and silicon.
Collapse
Affiliation(s)
- Marta L Vidal
- DTU Chemistry - Department of Chemistry, Technical University of Denmark, Kemitorvet Bldg 207, DK-2800 Kongens Lyngby, Denmark
| | - Pavel Pokhilko
- Department of Chemistry, University of Southern California, Los Angeles, California 90089-0482, United States
| | - Anna I Krylov
- Department of Chemistry, University of Southern California, Los Angeles, California 90089-0482, United States
| | - Sonia Coriani
- DTU Chemistry - Department of Chemistry, Technical University of Denmark, Kemitorvet Bldg 207, DK-2800 Kongens Lyngby, Denmark
| |
Collapse
|
25
|
Sarangi R, Vidal ML, Coriani S, Krylov AI. On the basis set selection for calculations of core-level states: different strategies to balance cost and accuracy. Mol Phys 2020. [DOI: 10.1080/00268976.2020.1769872] [Citation(s) in RCA: 24] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/24/2022]
Affiliation(s)
- Ronit Sarangi
- Department of Chemistry, University of Southern California, Los Angeles, CA, USA
| | - Marta L. Vidal
- DTU Chemistry – Department of Chemistry, Technical University of Denmark, Lyngby, Denmark
| | - Sonia Coriani
- DTU Chemistry – Department of Chemistry, Technical University of Denmark, Lyngby, Denmark
| | - Anna I. Krylov
- Department of Chemistry, University of Southern California, Los Angeles, CA, USA
| |
Collapse
|
26
|
Matthews DA. EOM-CC methods with approximate triple excitations applied to core excitation and ionisation energies. Mol Phys 2020. [DOI: 10.1080/00268976.2020.1771448] [Citation(s) in RCA: 12] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/24/2022]
|
27
|
Balasubramani SG, Chen GP, Coriani S, Diedenhofen M, Frank MS, Franzke YJ, Furche F, Grotjahn R, Harding ME, Hättig C, Hellweg A, Helmich-Paris B, Holzer C, Huniar U, Kaupp M, Marefat Khah A, Karbalaei Khani S, Müller T, Mack F, Nguyen BD, Parker SM, Perlt E, Rappoport D, Reiter K, Roy S, Rückert M, Schmitz G, Sierka M, Tapavicza E, Tew DP, van Wüllen C, Voora VK, Weigend F, Wodyński A, Yu JM. TURBOMOLE: Modular program suite for ab initio quantum-chemical and condensed-matter simulations. J Chem Phys 2020; 152:184107. [PMID: 32414256 PMCID: PMC7228783 DOI: 10.1063/5.0004635] [Citation(s) in RCA: 563] [Impact Index Per Article: 140.8] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/12/2020] [Accepted: 04/07/2020] [Indexed: 01/30/2023] Open
Abstract
TURBOMOLE is a collaborative, multi-national software development project aiming to provide highly efficient and stable computational tools for quantum chemical simulations of molecules, clusters, periodic systems, and solutions. The TURBOMOLE software suite is optimized for widely available, inexpensive, and resource-efficient hardware such as multi-core workstations and small computer clusters. TURBOMOLE specializes in electronic structure methods with outstanding accuracy-cost ratio, such as density functional theory including local hybrids and the random phase approximation (RPA), GW-Bethe-Salpeter methods, second-order Møller-Plesset theory, and explicitly correlated coupled-cluster methods. TURBOMOLE is based on Gaussian basis sets and has been pivotal for the development of many fast and low-scaling algorithms in the past three decades, such as integral-direct methods, fast multipole methods, the resolution-of-the-identity approximation, imaginary frequency integration, Laplace transform, and pair natural orbital methods. This review focuses on recent additions to TURBOMOLE's functionality, including excited-state methods, RPA and Green's function methods, relativistic approaches, high-order molecular properties, solvation effects, and periodic systems. A variety of illustrative applications along with accuracy and timing data are discussed. Moreover, available interfaces to users as well as other software are summarized. TURBOMOLE's current licensing, distribution, and support model are discussed, and an overview of TURBOMOLE's development workflow is provided. Challenges such as communication and outreach, software infrastructure, and funding are highlighted.
Collapse
Affiliation(s)
- Sree Ganesh Balasubramani
- Department of Chemistry, University of California, Irvine, 1102 Natural Sciences II, Irvine, California 92697-2025, USA
| | - Guo P Chen
- Department of Chemistry, University of California, Irvine, 1102 Natural Sciences II, Irvine, California 92697-2025, USA
| | - Sonia Coriani
- DTU Chemistry, Technical University of Denmark, Kemitorvet Build. 207, DK-2800 Kongens Lyngby, Denmark
| | - Michael Diedenhofen
- Dassault Systèmes Deutschland GmbH, Imbacher Weg 46, 51379 Leverkusen, Germany
| | - Marius S Frank
- Lehrstuhl für Theoretische Chemie, Ruhr-Universität Bochum, 44801 Bochum, Germany
| | - Yannick J Franzke
- Institute of Physical Chemistry, Karlsruhe Institute of Technology (KIT), KIT Campus South, P.O. Box 6980, 76049 Karlsruhe, Germany
| | - Filipp Furche
- Department of Chemistry, University of California, Irvine, 1102 Natural Sciences II, Irvine, California 92697-2025, USA
| | - Robin Grotjahn
- Institut für Chemie, Theoretische Chemie/Quantenchemie, Technische Universität Berlin, Sekr. C7, Straße des 17. Juni 135, 10623 Berlin, Germany
| | | | - Christof Hättig
- Lehrstuhl für Theoretische Chemie, Ruhr-Universität Bochum, 44801 Bochum, Germany
| | - Arnim Hellweg
- Dassault Systèmes Deutschland GmbH, Imbacher Weg 46, 51379 Leverkusen, Germany
| | - Benjamin Helmich-Paris
- Max-Planck-Institut für Kohlenforschung, Kaiser-Wilhelm-Platz 1, 45470 Mülheim an der Ruhr, Germany
| | - Christof Holzer
- Institute of Physical Chemistry, Karlsruhe Institute of Technology (KIT), KIT Campus South, P.O. Box 6980, 76049 Karlsruhe, Germany
| | - Uwe Huniar
- Dassault Systèmes Deutschland GmbH, Imbacher Weg 46, 51379 Leverkusen, Germany
| | - Martin Kaupp
- Institut für Chemie, Theoretische Chemie/Quantenchemie, Technische Universität Berlin, Sekr. C7, Straße des 17. Juni 135, 10623 Berlin, Germany
| | - Alireza Marefat Khah
- Lehrstuhl für Theoretische Chemie, Ruhr-Universität Bochum, 44801 Bochum, Germany
| | | | - Thomas Müller
- Forschungszentrum Jülich, Jülich Supercomputer Centre, Wilhelm-Jonen Straße, 52425 Jülich, Germany
| | - Fabian Mack
- Institute of Physical Chemistry, Karlsruhe Institute of Technology (KIT), KIT Campus South, P.O. Box 6980, 76049 Karlsruhe, Germany
| | - Brian D Nguyen
- Department of Chemistry, University of California, Irvine, 1102 Natural Sciences II, Irvine, California 92697-2025, USA
| | - Shane M Parker
- Department of Chemistry, Case Western Reserve University, 10900 Euclid Avenue, Cleveland, Ohio 44106, USA
| | - Eva Perlt
- Department of Chemistry, University of California, Irvine, 1102 Natural Sciences II, Irvine, California 92697-2025, USA
| | - Dmitrij Rappoport
- Department of Chemistry, University of North Carolina at Chapel Hill, Chapel Hill, North Carolina 27599, USA
| | - Kevin Reiter
- Institute of Nanotechnology, Karlsruhe Institute of Technology (KIT), KIT Campus North, P.O. Box 3640, 76021 Karlsruhe, Germany
| | - Saswata Roy
- Department of Chemistry, University of California, Irvine, 1102 Natural Sciences II, Irvine, California 92697-2025, USA
| | - Matthias Rückert
- Lehrstuhl für Theoretische Chemie, Ruhr-Universität Bochum, 44801 Bochum, Germany
| | - Gunnar Schmitz
- Department of Chemistry, Aarhus Universitet, Langelandsgade 140, DK-8000 Aarhus, Denmark
| | - Marek Sierka
- TURBOMOLE GmbH, Litzenhardtstraße 19, 76135 Karlsruhe, Germany
| | - Enrico Tapavicza
- Department of Chemistry and Biochemistry, California State University, Long Beach, 1250 Bellflower Boulevard, Long Beach, California 90840, USA
| | - David P Tew
- Max Planck Institute for Solid State Research, Heisenbergstaße 1, 70569 Stuttgart, Germany
| | - Christoph van Wüllen
- Fachbereich Chemie and Forschungszentrum OPTIMAS, Technische Universität Kaiserslautern, Erwin-Schrödinger-Staße 52, 67663 Kaiserslautern, Germany
| | - Vamsee K Voora
- Department of Chemical Sciences, Tata Institute of Fundamental Research, Homi Bhabha Road, Colaba, Mumbai 400005, India
| | - Florian Weigend
- Institute of Nanotechnology, Karlsruhe Institute of Technology (KIT), KIT Campus North, P.O. Box 3640, 76021 Karlsruhe, Germany
| | - Artur Wodyński
- Institut für Chemie, Theoretische Chemie/Quantenchemie, Technische Universität Berlin, Sekr. C7, Straße des 17. Juni 135, 10623 Berlin, Germany
| | - Jason M Yu
- Department of Chemistry, University of California, Irvine, 1102 Natural Sciences II, Irvine, California 92697-2025, USA
| |
Collapse
|
28
|
Zheng X, Liu J, Doumy G, Young L, Cheng L. Hetero-site Double Core Ionization Energies with Sub-electronvolt Accuracy from Delta-Coupled-Cluster Calculations. J Phys Chem A 2020; 124:4413-4426. [DOI: 10.1021/acs.jpca.0c00901] [Citation(s) in RCA: 9] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Xuechen Zheng
- Department of Chemistry, The Johns Hopkins University, Baltimore, Maryland 21218, United States
| | - Junzi Liu
- Department of Chemistry, The Johns Hopkins University, Baltimore, Maryland 21218, United States
| | - Gilles Doumy
- Chemical Sciences and Engineering Division, Argonne National Laboratory, Lemont, Illinois 60439, United States
| | - Linda Young
- Chemical Sciences and Engineering Division, Argonne National Laboratory, Lemont, Illinois 60439, United States
| | - Lan Cheng
- Department of Chemistry, The Johns Hopkins University, Baltimore, Maryland 21218, United States
| |
Collapse
|
29
|
Myhre RH, Coriani S, Koch H. X-ray and UV Spectra of Glycine within Coupled Cluster Linear Response Theory. J Phys Chem A 2019; 123:9701-9711. [PMID: 31549830 DOI: 10.1021/acs.jpca.9b06590] [Citation(s) in RCA: 12] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
The coupled cluster models CCSD and CC3 are used to investigate the (core) excited states and ionization energies of glycine in the gas phase. Excited states and ionization energies in the UV spectral range are calculated using a standard coupled cluster linear response, while core-level excited states and ionization potentials are calculated using the core-valence separation approximation. The temperature dependence from different conformers is also assessed.
Collapse
Affiliation(s)
- Rolf H Myhre
- Department of Chemistry , Norwegian University of Science and Technology, NTNU , 7491 Trondheim , Norway.,Hylleraas Centre for Quantum Molecular Sciences, Department of Chemistry , University of Oslo , 0315 Oslo , Norway
| | - Sonia Coriani
- DTU Chemistry , Technical University of Denmark , DK-2800 Kongens Lyngby , Denmark.,Aarhus Institute of Advanced Studies , Aarhus University , DK-8000 Aarhus C , Denmark
| | - Henrik Koch
- Department of Chemistry , Norwegian University of Science and Technology, NTNU , 7491 Trondheim , Norway.,Scuola Normale Superiore , Piazza dei Cavalieri 7 , 56126 Pisa , Italy
| |
Collapse
|
30
|
Faber R, Kjønstad EF, Koch H, Coriani S. Spin adapted implementation of EOM-CCSD for triplet excited states: Probing intersystem crossings of acetylacetone at the carbon and oxygen K-edges. J Chem Phys 2019; 151:144107. [PMID: 31615219 DOI: 10.1063/1.5112164] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Abstract
We present an equation of motion coupled cluster singles and doubles approach for computing transient absorption spectra from a triplet excited state. The implementation determines the left and right excitation vectors by explicitly spin-adapting the triplet excitation space. As an illustrative application, we compute transient state X-ray absorption spectra at the carbon and oxygen K-edges for the acetylacetone molecule.
Collapse
Affiliation(s)
- Rasmus Faber
- DTU Chemistry, Technical University of Denmark, Kemitorvet Build. 207, DK-2800 Kongens Lyngby, Denmark
| | - Eirik F Kjønstad
- Department of Chemistry, Norwegian University of Science and Technology, N-7491 Trondheim, Norway
| | - Henrik Koch
- Department of Chemistry, Norwegian University of Science and Technology, N-7491 Trondheim, Norway
| | - Sonia Coriani
- DTU Chemistry, Technical University of Denmark, Kemitorvet Build. 207, DK-2800 Kongens Lyngby, Denmark
| |
Collapse
|
31
|
Faber R, Coriani S. Resonant Inelastic X-ray Scattering and Nonesonant X-ray Emission Spectra from Coupled-Cluster (Damped) Response Theory. J Chem Theory Comput 2018; 15:520-528. [DOI: 10.1021/acs.jctc.8b01020] [Citation(s) in RCA: 32] [Impact Index Per Article: 5.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Rasmus Faber
- Department of Chemistry, Technical University of Denmark, Kemitorvet Building 207, 2800 Kongens Lyngby, Denmark
| | - Sonia Coriani
- Department of Chemistry, Technical University of Denmark, Kemitorvet Building 207, 2800 Kongens Lyngby, Denmark
| |
Collapse
|