1
|
Morgunov A, Tran HK, Meitei OR, Chien YC, Van Voorhis T. MP2-Based Composite Extrapolation Schemes Can Predict Core-Ionization Energies for First-Row Elements with Coupled-Cluster Level Accuracy. J Phys Chem A 2024; 128:6989-6998. [PMID: 39121455 DOI: 10.1021/acs.jpca.4c01606] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 08/11/2024]
Abstract
X-ray photoelectron spectroscopy (XPS) measures core-electron binding energies (CEBEs) to reveal element-specific insights into the chemical environment and bonding. Accurate theoretical CEBE prediction aids XPS interpretation but requires proper modeling of orbital relaxation and electron correlation upon core-ionization. This work systematically investigates basis set selection for extrapolation to the complete basis set limit of CEBEs from ΔMP2 and ΔCC energies across 94 K-edges in diverse organic molecules. We demonstrate that an alternative composite scheme using ΔMP2 in a large basis corrected by ΔCC-ΔMP2 difference in a small basis can quantitatively recover optimally extrapolated ΔCC CEBEs within 0.02 eV. Unlike ΔCC, MP2 calculations do not suffer from convergence issues and are computationally cheaper, and thus, the composite ΔMP2/ΔCC scheme balances accuracy and cost, overcoming limitations of solely using either method. We conclude by providing a comprehensive analysis of the choice of small and large basis sets for the composite schemes and provide practical recommendations for highly accurate (within 0.10-0.15 eV MAE) ab initio prediction of XPS data.
Collapse
Affiliation(s)
- Anton Morgunov
- Department of Chemistry, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, United States
| | - Henry K Tran
- Department of Chemistry, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, United States
| | - Oinam Romesh Meitei
- Department of Chemistry, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, United States
| | - Yu-Che Chien
- Department of Chemistry, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, United States
| | - Troy Van Voorhis
- Department of Chemistry, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, United States
| |
Collapse
|
2
|
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
|
3
|
Mandal A, Berquist EJ, Herbert JM. A new parameterization of the DFT/CIS method with applications to core-level spectroscopy. J Chem Phys 2024; 161:044114. [PMID: 39051834 DOI: 10.1063/5.0220535] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/26/2024] [Accepted: 07/08/2024] [Indexed: 07/27/2024] Open
Abstract
Time-dependent density functional theory (TD-DFT) within a restricted excitation space is an efficient means to compute core-level excitation energies using only a small subset of the occupied orbitals. However, core-to-valence excitation energies are significantly underestimated when standard exchange-correlation functionals are used, which is partly traceable to systemic issues with TD-DFT's description of Rydberg and charge-transfer excited states. To mitigate this, we have implemented an empirically modified combination of configuration interaction with single substitutions (CIS) based on Kohn-Sham orbitals, which is known as "DFT/CIS." This semi-empirical approach is well-suited for simulating x-ray near-edge spectra, as it contains sufficient exact exchange to model charge-transfer excitations yet retains DFT's low-cost description of dynamical electron correlation. Empirical corrections to the matrix elements enable semi-quantitative simulation of near-edge x-ray spectra without the need for significant a posteriori shifts; this should be useful in complex molecules and materials with multiple overlapping x-ray edges. Parameter optimization for use with a specific range-separated hybrid functional makes this a black-box method intended for both core and valence spectroscopy. Results herein demonstrate that realistic K-edge absorption and emission spectra can be obtained for second- and third-row elements and 3d transition metals, with promising results for L-edge spectra as well. DFT/CIS calculations require absolute shifts that are considerably smaller than what is typical in TD-DFT.
Collapse
Affiliation(s)
- Aniket Mandal
- Department of Chemistry and Biochemistry, The Ohio State University, Columbus, Ohio 43210, USA
| | | | - John M Herbert
- Department of Chemistry and Biochemistry, The Ohio State University, Columbus, Ohio 43210, USA
| |
Collapse
|
4
|
Liu L, Xu Q, dos Anjos Cunha L, Xin H, Head-Gordon M, Qian J. Real-Space Pseudopotential Method for the Calculation of Third-Row Elements X-ray Photoelectron Spectroscopic Signatures. J Chem Theory Comput 2024; 20:6134-6143. [PMID: 38970155 PMCID: PMC11270745 DOI: 10.1021/acs.jctc.4c00535] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/21/2024] [Revised: 06/18/2024] [Accepted: 06/20/2024] [Indexed: 07/08/2024]
Abstract
X-ray photoelectron spectroscopy (XPS) is a powerful characterization technique that unveils subtle chemical environment differences via core-electron binding energy (CEBE) analysis. We extend the development of real-space pseudopotential methods to calculating 1s, 2s, and 2p3/2 CEBEs of third-row elements (S, P, and Si) within the framework of Kohn-Sham density-functional theory (KS-DFT). The new approach systematically prevents variational collapse and simplifies core-excited orbital selection within dense energy level distributions. However, careful error cancellation analysis is required to achieve accuracy comparable to all-electron methods and experiments. Combined with real-space KS-DFT implementation, this development enables large-scale simulations with both Dirichlet boundary conditions and periodic boundary conditions.
Collapse
Affiliation(s)
- Liping Liu
- Department
of Chemical Engineering, Virginia Polytechnic
Institute and State University, Blacksburg, Virginia 24060, United States
| | - Qiang Xu
- Chemical
Sciences Division, Lawrence Berkeley National
Laboratory, Berkeley, California 94720, United States
| | - Leonardo dos Anjos Cunha
- Chemical
Sciences Division, Lawrence Berkeley National
Laboratory, Berkeley, California 94720, United States
- Department
of Chemistry, University of California,
Berkeley, Berkeley, California 94720, United States
| | - Hongliang Xin
- Department
of Chemical Engineering, Virginia Polytechnic
Institute and State University, Blacksburg, Virginia 24060, United States
| | - Martin Head-Gordon
- Chemical
Sciences Division, Lawrence Berkeley National
Laboratory, Berkeley, California 94720, United States
- Department
of Chemistry, University of California,
Berkeley, Berkeley, California 94720, United States
| | - Jin Qian
- Chemical
Sciences Division, Lawrence Berkeley National
Laboratory, Berkeley, California 94720, United States
| |
Collapse
|
5
|
de Moura CEV, Sokolov AY. Efficient Spin-Adapted Implementation of Multireference Algebraic Diagrammatic Construction Theory. I. Core-Ionized States and X-ray Photoelectron Spectra. J Phys Chem A 2024; 128:5816-5831. [PMID: 38962857 DOI: 10.1021/acs.jpca.4c03161] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 07/05/2024]
Abstract
We present an efficient implementation of multireference algebraic diagrammatic construction theory (MR-ADC) for simulating core-ionized states and X-ray photoelectron spectra (XPS). Taking advantage of spin adaptation, automatic code generation, and density fitting, our implementation can perform calculations for molecules with more than 1500 molecular orbitals, incorporating static and dynamic correlation in the ground and excited electronic states. We demonstrate the capabilities of MR-ADC methods by simulating the XPS spectra of substituted ferrocene complexes and azobenzene isomers. For the ground electronic states of these molecules, the XPS spectra computed using the extended second-order MR-ADC method (MR-ADC(2)-X) are in a very good agreement with available experimental results. We further show that MR-ADC can be used as a tool for interpreting or predicting the results of time-resolved XPS measurements by simulating the core ionization spectra of azobenzene along its photoisomerization, including the XPS signatures of excited states and the minimum energy conical intersection. This work is the first in a series of publications reporting the efficient implementations of MR-ADC methods.
Collapse
Affiliation(s)
- Carlos E V de Moura
- 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
|
6
|
Nykänen A, Thiessen L, Borrelli EM, Krishna V, Knecht S, Pavošević F. Toward Accurate Calculation of Excitation Energies on Quantum Computers with ΔADAPT-VQE: A Case Study of BODIPY Derivatives. J Phys Chem Lett 2024; 15:7111-7117. [PMID: 38954795 DOI: 10.1021/acs.jpclett.4c01301] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 07/04/2024]
Abstract
Quantum chemistry simulations offer a cost-effective way to computationally design BODIPY photosensitizers. However, accurate predictions of excitation energies pose a challenge for time-dependent density functional theory and equation-of-motion coupled-cluster singles and doubles methods. By contrast, reliable predictions can be achieved by multireference quantum chemistry methods; unfortunately, their computational cost increases exponentially with the number of electrons. Alternatively, quantum computing holds potential for an exact simulation of the photophysical properties in a computationally more efficient way. Herein, we introduce the state-specific ΔUCCSD-VQE (unitary coupled-cluster singles and doubles-variational quantum eigensolver) and ΔADAPT-VQE methods in which the electronically excited state is calculated via a non-Aufbau configuration. We show for six BODIPY derivatives that the proposed methods predict accurate excitation energies that are in good agreement with those from experiments. Due to its performance and simplicity, we believe that ΔADAPT will become a useful approach for the simulation of BODIPY photosensitizers on near-term quantum devices.
Collapse
Affiliation(s)
- Anton Nykänen
- Algorithmiq Ltd., Kanavakatu 3C, FI-00160 Helsinki, Finland
| | | | | | - Vijay Krishna
- Department of Biomedical Engineering, Lerner Research Institute, Cleveland Clinic, Cleveland, Ohio 44195, United States
- Department of Biomedical Engineering, Cleveland Clinic Lerner College of Medicine, Case Western Reserve University, Cleveland, Ohio 44106, United States
| | - Stefan Knecht
- Algorithmiq Ltd., Kanavakatu 3C, FI-00160 Helsinki, Finland
| | | |
Collapse
|
7
|
Park W, Lashkaripour A, Komarov K, Lee S, Huix-Rotllant M, Choi CH. Toward Consistent Predictions of Core/Valence Ionization Potentials and Valence Excitation Energies by MRSF-TDDFT. J Chem Theory Comput 2024; 20:5679-5694. [PMID: 38902891 DOI: 10.1021/acs.jctc.4c00640] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/22/2024]
Abstract
Optimizing exchange-correlation functionals for both core/valence ionization potentials (cIPs/vIPs) and valence excitation energies (VEEs) at the same time in the framework of MRSF-TDDFT is self-contradictory. To overcome the challenge, within the previous "adaptive exact exchange" or double-tuning strategy on Coulomb-attenuating XC functionals (CAM), a new XC functional specifically for cIPs and vIPs was first developed by enhancing exact exchange to both short- and long-range regions. The resulting DTCAM-XI functional achieved remarkably high accuracy in its predictions with errors of less than half eV. An additional concept of "valence attenuation", where the amount of exact exchange for the frontier orbital regions is selectively suppressed, was introduced to consistently predict both VEEs and IPs at the same time. The second functional, DTCAM-XIV, exhibits consistent overall prediction accuracy at ∼0.64 eV. By preferentially optimizing VEEs within the same "valence attenuation" concept, a third functional, DTCAM-VAEE, was obtained, which exhibits improved performance as compared to that of the previous DTCAM-VEE and DTCAM-AEE in the prediction of VEEs, making it an attractive alternative to BH&HLYP. As the combination of "adaptive exchange" and "valence attenuation" is operative, it would be exciting to explore its potential with a more tunable framework in the future.
Collapse
Affiliation(s)
- Woojin Park
- Department of Chemistry, Kyungpook National University, Daegu 41566, South Korea
| | - Alireza Lashkaripour
- Department of Chemistry, Kyungpook National University, Daegu 41566, South Korea
| | - Konstantin Komarov
- Center for Quantum Dynamics, Pohang University of Science and Technology, Pohang 37673, South Korea
- Department of Chemistry, University of Zürich, Zürich 8057, Switzerland
| | - Seunghoon Lee
- Department of Chemistry, Seoul National University, Seoul 151-747, South Korea
| | | | - Cheol Ho Choi
- Department of Chemistry, Kyungpook National University, Daegu 41566, South Korea
| |
Collapse
|
8
|
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
|
9
|
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
|
10
|
Damour Y, Scemama A, Jacquemin D, Kossoski F, Loos PF. State-Specific Coupled-Cluster Methods for Excited States. J Chem Theory Comput 2024; 20:4129-4145. [PMID: 38749498 PMCID: PMC11137840 DOI: 10.1021/acs.jctc.4c00034] [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/10/2024] [Revised: 03/06/2024] [Accepted: 03/06/2024] [Indexed: 05/29/2024]
Abstract
We reexamine ΔCCSD, a state-specific coupled-cluster (CC) with single and double excitations (CCSD) approach that targets excited states through the utilization of non-Aufbau determinants. This methodology is particularly efficient when dealing with doubly excited states, a domain in which the standard equation-of-motion CCSD (EOM-CCSD) formalism falls short. Our goal here to evaluate the effectiveness of ΔCCSD when applied to other types of excited states, comparing its consistency and accuracy with EOM-CCSD. To this end, we report a benchmark on excitation energies computed with the ΔCCSD and EOM-CCSD methods for a set of molecular excited-state energies that encompasses not only doubly excited states but also doublet-doublet transitions and (singlet and triplet) singly excited states of closed-shell systems. In the latter case, we rely on a minimalist version of multireference CC known as the two-determinant CCSD method to compute the excited states. Our data set, consisting of 276 excited states stemming from the quest database [Véril et al., WIREs Comput. Mol. Sci. 2021, 11, e1517], provides a significant base to draw general conclusions concerning the accuracy of ΔCCSD. Except for the doubly excited states, we found that ΔCCSD underperforms EOM-CCSD. For doublet-doublet transitions, the difference between the mean absolute errors (MAEs) of the two methodologies (of 0.10 and 0.07 eV) is less pronounced than that obtained for singly excited states of closed-shell systems (MAEs of 0.15 and 0.08 eV). This discrepancy is largely attributed to a greater number of excited states in the latter set exhibiting multiconfigurational characters, which are more challenging for ΔCCSD. We also found typically small improvements by employing state-specific optimized orbitals.
Collapse
Affiliation(s)
- Yann Damour
- Laboratoire
de Chimie et Physique Quantiques (UMR 5626), Université de Toulouse, CNRS, UPS, 31000 Toulouse, France
| | - Anthony Scemama
- Laboratoire
de Chimie et Physique Quantiques (UMR 5626), Université de Toulouse, CNRS, UPS, 31000 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, 31000 Toulouse, France
| | - Pierre-François Loos
- Laboratoire
de Chimie et Physique Quantiques (UMR 5626), Université de Toulouse, CNRS, UPS, 31000 Toulouse, France
| |
Collapse
|
11
|
Schraivogel T, Kats D. Two determinant distinguishable cluster. J Chem Phys 2024; 160:124109. [PMID: 38526108 DOI: 10.1063/5.0199274] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/22/2024] [Accepted: 03/06/2024] [Indexed: 03/26/2024] Open
Abstract
A two reference determinant version of the distinguishable cluster with singles and doubles (DCSD) has been developed. We have implemented the two determinant distinguishable cluster (2D-DCSD) and the corresponding traditional 2D-CCSD method in a new open-source package written in Julia called ElemCo.jl. The methods were benchmarked on singlet and triplet excited states of valence and Rydberg character, as well as for singlet-triplet gaps of diradicals. It is demonstrated that the distinguishable cluster approximation improves the accuracy of 2D-CCSD.
Collapse
Affiliation(s)
- Thomas Schraivogel
- Max Planck Institute for Solid State Research, Heisenbergstraße 1, 70569 Stuttgart, Germany
| | - Daniel Kats
- Max Planck Institute for Solid State Research, Heisenbergstraße 1, 70569 Stuttgart, Germany
| |
Collapse
|
12
|
Haasler M, Maier TM, Kaupp M. Toward a correct treatment of core properties with local hybrid functionals. J Comput Chem 2023; 44:2461-2477. [PMID: 37635647 DOI: 10.1002/jcc.27211] [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: 07/04/2023] [Revised: 08/01/2023] [Accepted: 08/03/2023] [Indexed: 08/29/2023]
Abstract
In local hybrid functionals (LHs), a local mixing function (LMF) determines the position-dependent exact-exchange admixture. We report new LHs that focus on an improvement of the LMF in the core region while retaining or partly improving upon the high accuracy in the valence region exhibited by the LH20t functional. The suggested new pt-LMFs are based on a Padé form and modify the previously used ratio between von Weizsäcker and Kohn-Sham local kinetic energies by different powers of the density to enable flexibly improved approximations to the correct high-density and iso-orbital limits relevant for the innermost core region. Using TDDFT calculations for a set of K-shell core excitations of second- and third-period systems including accurate state-of-the-art relativistic orbital corrections, the core part of the LMF is optimized, while the valence part is optimized as previously reported for test sets of atomization energies and reaction barriers (Haasler et al., J Chem Theory Comput 2020, 16, 5645). The LHs are completed by a calibration function that minimizes spurious nondynamical correlation effects caused by the gauge ambiguities of exchange-energy densities, as well as by B95c meta-GGA correlation. The resulting LH23pt functional relates to the previous LH20t functional but specifically improves upon the core region.
Collapse
Affiliation(s)
- Matthias Haasler
- Technische Universität Berlin, Institute of Chemistry Theoretical Chemistry/Quantum Chemistry, Berlin, Germany
| | - Toni M Maier
- Technische Universität Braunschweig, Institute of Physical and Theoretical Chemistry, Braunschweig, Germany
| | - Martin Kaupp
- Technische Universität Berlin, Institute of Chemistry Theoretical Chemistry/Quantum Chemistry, Berlin, Germany
| |
Collapse
|
13
|
Ferté A, Giner E, Taïeb R, Carniato S. Unraveling the variational breakdown of core valence separation calculations: Diagnostic and cure to the over relaxation error of double core hole states. J Chem Phys 2023; 159:144104. [PMID: 37811825 DOI: 10.1063/5.0159493] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/24/2023] [Accepted: 09/21/2023] [Indexed: 10/10/2023] Open
Abstract
The core valence separation (CVS) approximation is the most employed strategy to prevent the variational collapse of standard wave function optimization when attempting to compute electronic states bearing one or more electronic vacancies in core orbitals. Here, we explore the spurious consequences of this approximation on the properties of the computed core hole states. We especially focus on the less studied case of double core hole (DCH) states, whose spectroscopic interest has recently been rapidly growing. We show that the CVS error leads to a systematic underestimation of DCH energies, a property in stark contrast with the case of single core hole states. We highlight that the CVS error can then be interpreted as an over relaxation effect and design a new correction strategy adapted to these specificities.
Collapse
Affiliation(s)
- Anthony Ferté
- Laboratoire de Chimie Physique-Matière et Rayonnement, Sorbonne Université and CNRS, F-75005 Paris, France
- Laboratoire de Chimie et Interdisciplinarité, Synthèse, Analyse, Modélisation, Nantes Université and CNRS, F-44000 Nantes, France
| | - Emmanuel Giner
- Laboratoire de Chimie Théorique, Sorbonne Université and CNRS, F-75005 Paris, France
| | - Richard Taïeb
- Laboratoire de Chimie Physique-Matière et Rayonnement, Sorbonne Université and CNRS, F-75005 Paris, France
| | - Stéphane Carniato
- Laboratoire de Chimie Physique-Matière et Rayonnement, Sorbonne Université and CNRS, F-75005 Paris, France
| |
Collapse
|
14
|
Misael WA, Severo Pereira Gomes A. Core Excitations of Uranyl in Cs 2UO 2Cl 4 from Relativistic Embedded Damped Response Time-Dependent Density Functional Theory Calculations. Inorg Chem 2023; 62:11589-11601. [PMID: 37432868 DOI: 10.1021/acs.inorgchem.3c01302] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 07/13/2023]
Abstract
X-ray spectroscopies, by their high selectivity and sensitivity to the chemical environment around the atoms probed, provide significant insights into the electronic structures of molecules and materials. Interpreting experimental results requires reliable theoretical models, accounting for environmental, relativistic, electron correlation, and orbital relaxation effects in a balanced manner. In this work, we present a protocol for the simulation of core excited spectra with damped response time-dependent density functional theory based on the Dirac-Coulomb Hamiltonian (4c-DR-TD-DFT), in which environmental effects are accounted for through the frozen density embedding (FDE) method. We showcase this approach for the uranium M4- and L3-edges and oxygen K-edge of the uranyl tetrachloride (UO2Cl42-) unit as found in a host Cs2UO2Cl4 crystal. We have found that the 4c-DR-TD-DFT simulations yield excitation spectra that very closely match the experiment for the uranium M4-edge and the oxygen K-edge, with good agreement for the broad experimental spectra for the L3-edge. By decomposing the complex polarizability in terms of its components, we have been able to correlate our results with angle-resolved spectra. We have observed that for all edges, but in particular the uranium M4-edge, an embedded model in which the chloride ligands are replaced by an embedding potential reproduces rather well the spectral profile obtained for UO2Cl42-. Our results underscore the importance of the equatorial ligands to simulating core spectra at both uranium and oxygen edges.
Collapse
Affiliation(s)
- Wilken Aldair Misael
- Univ. Lille, CNRS, UMR 8523-PhLAM-Physique des Lasers Atomes et Molécules, F-59000 Lille, France
| | | |
Collapse
|
15
|
Hou SJ, Yang YF, Cui ZH, Cederbaum LS. Can anions possess bound doubly-excited electronic states? Chem Sci 2023; 14:7230-7236. [PMID: 37416703 PMCID: PMC10321500 DOI: 10.1039/d3sc00370a] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/21/2023] [Accepted: 05/27/2023] [Indexed: 07/08/2023] Open
Abstract
Anions play an important role in many fields of chemistry. Many molecules possess stable anions, but these anions often do not have stable electronic excited states and the anion loses its excess electron once excited. All the known stable valence excited states of anions are singly-excited states, i.e., valence doubly-excited states have not been reported. As excited states are relevant for numerous applications, and constitute basic properties, we searched for valence doubly-excited states which are stable, i.e., exhibit energies below that of the ground state of the respective neutral molecule. We concentrated on two promising prototype candidates, the anions of the smallest endocircular carbon ring Li@C12 and of the smallest endohedral fullerene Li@C20. By employing accurate state-of-the-art many-electron quantum chemistry methods, we investigated the low-lying excited states of these anions and found that they possess several low-lying stable singly-excited states and, in particular, a stable doubly-excited state each. It is noteworthy that the found doubly-excited state of Li@C12- possesses a cumulenic carbon ring in sharp contrast to the ground and singly-excited states. The findings shed light on how to design anions with stable valence singly- and doubly-excited states. Possible applications are mentioned.
Collapse
Affiliation(s)
- Shi-Jie Hou
- Institute of Atomic and Molecular Physics, Jilin University Changchun 130023 China
| | - Yi-Fan Yang
- Quantum Theory Project, Departments of Physics and Chemistry, University of Florida Gainesville Florida 32611 USA
| | - Zhong-Hua Cui
- Institute of Atomic and Molecular Physics, Jilin University Changchun 130023 China
| | - Lorenz S Cederbaum
- Theoretical Chemistry, Institute of Physical Chemistry, Universität Heidelberg Im Neuenheimer Feld 229 D-69120 Heidelberg Germany
| |
Collapse
|
16
|
Mester D, Kállay M. Vertical Ionization Potentials and Electron Affinities at the Double-Hybrid Density Functional Level. J Chem Theory Comput 2023. [PMID: 37326360 DOI: 10.1021/acs.jctc.3c00363] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/17/2023]
Abstract
The double-hybrid (DH) time-dependent density functional theory is extended to vertical ionization potentials (VIPs) and electron affinities (VEAs). Utilizing the density fitting approximation, efficient implementations are presented for the genuine DH ansatz relying on the perturbative second-order correction, while an iterative analogue is also elaborated using our second-order algebraic-diagrammatic construction [ADC(2)]-based DH approach. The favorable computational requirements of the present schemes are discussed in detail. The performance of the recently proposed spin-component-scaled and spin-opposite-scaled (SOS) range-separated (RS) and long-range corrected (LC) DH functionals is comprehensively assessed, while popular hybrid and global DH approaches are also discussed. For the benchmark calculations, up-to-date test sets are selected with high-level coupled-cluster references. Our results show that the ADC(2)-based SOS-RS-PBE-P86 approach is the most accurate and robust functional. This method consistently outperforms the excellent SOS-ADC(2) approach for VIPs, although the results are somewhat less satisfactory for VEAs. Among the genuine DH functionals, the SOS-ωPBEPP86 approach is also recommended for describing ionization processes, but its performance is even less reliable for electron-attached states. In addition, surprisingly good results are attained by the LC hybrid ωB97X-D functional, where the corresponding occupied (unoccupied) orbital energies are retrieved as VIPs (VEAs) within the present formalism.
Collapse
Affiliation(s)
- Dávid Mester
- Department of Physical Chemistry and Materials Science, Faculty of Chemical Technology and Biotechnology, Budapest University of Technology and Economics, Műegyetem rkp. 3., H-1111 Budapest, Hungary
- ELKH-BME Quantum Chemistry Research Group, Műegyetem rkp. 3., H-1111 Budapest, Hungary
- MTA-BME Lendület Quantum Chemistry Research Group, Műegyetem rkp. 3., H-1111 Budapest, Hungary
| | - Mihály Kállay
- Department of Physical Chemistry and Materials Science, Faculty of Chemical Technology and Biotechnology, Budapest University of Technology and Economics, Műegyetem rkp. 3., H-1111 Budapest, Hungary
- ELKH-BME Quantum Chemistry Research Group, Műegyetem rkp. 3., H-1111 Budapest, Hungary
- MTA-BME Lendület Quantum Chemistry Research Group, Műegyetem rkp. 3., H-1111 Budapest, Hungary
| |
Collapse
|
17
|
Banerjee S, Sokolov AY. Algebraic Diagrammatic Construction Theory for Simulating Charged Excited States and Photoelectron Spectra. J Chem Theory Comput 2023. [PMID: 37191264 DOI: 10.1021/acs.jctc.3c00251] [Citation(s) in RCA: 3] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/17/2023]
Abstract
Charged excitations are electronic transitions that involve a change in the total charge of a molecule or material. Understanding the properties and reactivity of charged species requires insights from theoretical calculations that can accurately describe orbital relaxation and electron correlation effects in open-shell electronic states. In this Review, we describe the current state of algebraic diagrammatic construction (ADC) theory for simulating charged excitations and its recent developments. We start with a short overview of ADC formalism for the one-particle Green's function, including its single- and multireference formulations and extension to periodic systems. Next, we focus on the capabilities of ADC methods and discuss recent findings about their accuracy for calculating a wide range of excited-state properties. We conclude our Review by outlining possible directions for future developments of this theoretical approach.
Collapse
Affiliation(s)
- Samragni Banerjee
- 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
|
18
|
Kahk JM, Lischner J. Combining the Δ-Self-Consistent-Field and GW Methods for Predicting Core Electron Binding Energies in Periodic Solids. J Chem Theory Comput 2023. [PMID: 37163299 DOI: 10.1021/acs.jctc.3c00121] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/11/2023]
Abstract
For the computational prediction of core electron binding energies in solids, two distinct kinds of modeling strategies have been pursued: the Δ-Self-Consistent-Field method based on density functional theory (DFT), and the GW method. In this study, we examine the formal relationship between these two approaches and establish a link between them. The link arises from the equivalence, in DFT, between the total energy difference result for the first ionization energy, and the eigenvalue of the highest occupied state, in the limit of infinite supercell size. This link allows us to introduce a new formalism, which highlights how in DFT─even if the total energy difference method is used to calculate core electron binding energies─the accuracy of the results still implicitly depends on the accuracy of the eigenvalue at the valence band maximum in insulators, or at the Fermi level in metals. We examine whether incorporating a quasiparticle correction for this eigenvalue from GW theory improves the accuracy of the calculated core electron binding energies, and find that the inclusion of vertex corrections is required for achieving quantitative agreement with experiment.
Collapse
Affiliation(s)
- J Matthias Kahk
- Institute of Physics, University of Tartu, W. Ostwaldi 1, 50411 Tartu, Estonia
| | - Johannes Lischner
- Department of Physics, Department of Materials, and the Thomas Young Centre for Theory and Simulation of Materials, Imperial College London, London SW7 2AZ, United Kingdom
| |
Collapse
|
19
|
Jana S, Herbert JM. Slater transition methods for core-level electron binding energies. J Chem Phys 2023; 158:094111. [PMID: 36889976 DOI: 10.1063/5.0134459] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/18/2023] Open
Abstract
Methods for computing core-level ionization energies using self-consistent field (SCF) calculations are evaluated and benchmarked. These include a "full core hole" (or "ΔSCF") approach that fully accounts for orbital relaxation upon ionization, but also methods based on Slater's transition concept in which the binding energy is estimated from an orbital energy level that is obtained from a fractional-occupancy SCF calculation. A generalization that uses two different fractional-occupancy SCF calculations is also considered. The best of the Slater-type methods afford mean errors of 0.3-0.4 eV with respect to experiment for a dataset of K-shell ionization energies, a level of accuracy that is competitive with more expensive many-body techniques. An empirical shifting procedure with one adjustable parameter reduces the average error below 0.2 eV. This shifted Slater transition method is a simple and practical way to compute core-level binding energies using only initial-state Kohn-Sham eigenvalues. It requires no more computational effort than ΔSCF and may be especially useful for simulating transient x-ray experiments where core-level spectroscopy is used to probe an excited electronic state, for which the ΔSCF approach requires a tedious state-by-state calculation of the spectrum. As an example, we use Slater-type methods to model x-ray emission spectroscopy.
Collapse
Affiliation(s)
- Subrata Jana
- Department of Chemistry and Biochemistry, The Ohio State University, Columbus, Ohio 43210, USA
| | - John M Herbert
- Department of Chemistry and Biochemistry, The Ohio State University, Columbus, Ohio 43210, USA
| |
Collapse
|
20
|
Datar A, Wright C, Matthews DA. Theoretical Investigation of the X-ray Stark Effect in Small Molecules. J Phys Chem A 2023; 127:1576-1587. [PMID: 36787229 DOI: 10.1021/acs.jpca.2c08311] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/15/2023]
Abstract
We have studied the Stark effect in the soft x-ray region for various small molecules by calculating the field-dependent x-ray absorption spectra. This effect is explained in terms of the response of molecular orbitals (core and valence), the molecular dipole moment, and the molecular geometry to the applied electric field. A number of consistent trends are observed linking the computed shifts in absorption energies and intensities with specific features of the molecular electronic structure. We find that both the virtual molecular orbitals (valence and/or Rydberg) as well as the core orbitals contribute to observed trends in a complementary fashion. This initial study highlights the potential impact of x-ray Stark spectroscopy as a tool to study electronic structure and environmental perturbations at a submolecular scale.
Collapse
Affiliation(s)
- Avdhoot Datar
- Department of Chemistry, Southern Methodist University, Dallas, Texas 75275, United States
| | - Catherine Wright
- 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
|
21
|
Mukatayev I, Moevus F, Sklénard B, Olevano V, Li J. XPS Core-Level Chemical Shift by Ab Initio Many-Body Theory. J Phys Chem A 2023; 127:1642-1648. [PMID: 36787463 DOI: 10.1021/acs.jpca.3c00173] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/16/2023]
Abstract
X-ray photoemission spectroscopy (XPS) provides direct information on atomic composition and stoichiometry by measuring core-electron binding energies. Moreover, from the shift of the binding energy, the so-called chemical shift, the precise chemical type of bonds can be inferred, which brings additional information on the local structure. In this work, we present a theoretical study of the chemical shift first by comparing different theories, from Hartree-Fock and density functional theory to many-body perturbation theory approaches like the GW approximation and its static version (COHSEX). The accuracy of each theory is assessed with respect to a carbon 1s chemical shift experimental benchmark measured on a set of gas-phase molecules. More importantly, by decomposing the chemical shift into different contributions according to terms in the total Hamiltonian, classical electrostatics is identified as the major contributor to the chemical shift, one order of magnitude larger than the correlation.
Collapse
Affiliation(s)
| | - Florient Moevus
- Université Grenoble Alpes, CEA, Leti, F-38000, Grenoble, France
| | - Benoît Sklénard
- Université Grenoble Alpes, CEA, Leti, F-38000, Grenoble, France.,European Theoretical Spectroscopy Facility (ETSF), bâtiment B5a Université de Liège Allée du 6 août, numéro 17 Sart-Tilman, F-38000 Grenoble, France
| | - Valerio Olevano
- European Theoretical Spectroscopy Facility (ETSF), bâtiment B5a Université de Liège Allée du 6 août, numéro 17 Sart-Tilman, F-38000 Grenoble, France.,Université Grenoble Alpes, F-38000 Grenoble, France.,CNRS, Institut Néel, F-38042 Grenoble, France
| | - Jing Li
- Université Grenoble Alpes, CEA, Leti, F-38000, Grenoble, France.,European Theoretical Spectroscopy Facility (ETSF), bâtiment B5a Université de Liège Allée du 6 août, numéro 17 Sart-Tilman, F-38000 Grenoble, France
| |
Collapse
|
22
|
Li J, Jin Y, Rinke P, Yang W, Golze D. Benchmark of GW Methods for Core-Level Binding Energies. J Chem Theory Comput 2022; 18:7570-7585. [PMID: 36322136 PMCID: PMC9753590 DOI: 10.1021/acs.jctc.2c00617] [Citation(s) in RCA: 19] [Impact Index Per Article: 9.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/14/2022]
Abstract
The GW approximation has recently gained increasing attention as a viable method for the computation of deep core-level binding energies as measured by X-ray photoelectron spectroscopy. We present a comprehensive benchmark study of different GW methodologies (starting point optimized, partial and full eigenvalue-self-consistent, Hedin shift, and renormalized singles) for molecular inner-shell excitations. We demonstrate that all methods yield a unique solution and apply them to the CORE65 benchmark set and ethyl trifluoroacetate. Three GW schemes clearly outperform the other methods for absolute core-level energies with a mean absolute error of 0.3 eV with respect to experiment. These are partial eigenvalue self-consistency, in which the eigenvalues are only updated in the Green's function, single-shot GW calculations based on an optimized hybrid functional starting point, and a Hedin shift in the Green's function. While all methods reproduce the experimental relative binding energies well, the eigenvalue self-consistent schemes and the Hedin shift yield with mean absolute errors <0.2 eV the best results.
Collapse
Affiliation(s)
- Jiachen Li
- Department
of Chemistry, Duke University, Durham, North Carolina27708, United States
| | - Ye Jin
- Department
of Chemistry, Duke University, Durham, North Carolina27708, United States
| | - Patrick Rinke
- Department
of Applied Physics, Aalto University, Otakaari 1, FI-02150Espoo, Finland
| | - Weitao Yang
- Department
of Chemistry, Duke University, Durham, North Carolina27708, United States
| | - Dorothea Golze
- Department
of Applied Physics, Aalto University, Otakaari 1, FI-02150Espoo, Finland,Faculty
of Chemistry and Food Chemistry, Technische
Universität Dresden, 01062Dresden, Germany,
| |
Collapse
|
23
|
Galleni L, Sajjadian FS, Conard T, Escudero D, Pourtois G, van Setten MJ. Modeling X-ray Photoelectron Spectroscopy of Macromolecules Using GW. J Phys Chem Lett 2022; 13:8666-8672. [PMID: 36084286 DOI: 10.1021/acs.jpclett.2c01935] [Citation(s) in RCA: 7] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/15/2023]
Abstract
We propose a simple additive approach to simulate X-ray photoelectron spectra (XPS) of macromolecules based on the GW method. Single-shot GW (G0W0) is a promising technique to compute accurate core-electron binding energies (BEs). However, its application to large molecules is still unfeasible. To circumvent the computational cost of G0W0, we break the macromolecule into tractable building blocks, such as isolated monomers, and sum up the theoretical spectra of each component, weighted by their molar ratio. In this work, we provide a first proof of concept by applying the method to four test polymers and one copolymer and show that it leads to an excellent agreement with experiments. The method could be used to retrieve the composition of unknown materials and study chemical reactions, by comparing the simulated spectra with experimental ones.
Collapse
Affiliation(s)
- Laura Galleni
- Department of Chemistry, KU Leuven, Celestijnenlaan 200F, 3001 Leuven, Belgium
- Imec, Kapeldreef 75, 3001 Leuven, Belgium
| | - Faegheh S Sajjadian
- Department of Chemistry, KU Leuven, Celestijnenlaan 200F, 3001 Leuven, Belgium
- Imec, Kapeldreef 75, 3001 Leuven, Belgium
| | | | - Daniel Escudero
- Department of Chemistry, KU Leuven, Celestijnenlaan 200F, 3001 Leuven, Belgium
| | | | | |
Collapse
|
24
|
Andersen JH, Nanda KD, Krylov AI, Coriani S. Cherry-Picking Resolvents: Recovering the Valence Contribution in X-ray Two-Photon Absorption within the Core-Valence-Separated Equation-of-Motion Coupled-Cluster Response Theory. J Chem Theory Comput 2022; 18:6189-6202. [PMID: 36084326 DOI: 10.1021/acs.jctc.2c00541] [Citation(s) in RCA: 6] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
Calculations of first-order response wave functions in the X-ray regime often diverge within correlated frameworks such as equation-of-motion coupled-cluster singles and doubles (EOM-CCSD), a consequence of the coupling with the valence ionization continuum. Here, we extend our strategy of introducing a hierarchy of approximations to the EOM-EE-CCSD resolvent (or, inversely, the model Hamiltonian) involved in the response equations for the calculation of X-ray two-photon absorption (X2PA) cross sections. We exploit the frozen-core core-valence separation (fc-CVS) scheme to first decouple the core and valence Fock spaces, followed by a separate approximate treatment of the valence resolvent. We demonstrate the robust convergence of X-ray response calculations within this framework and compare X2PA spectra of small benchmark molecules with the previously reported density functional theory results.
Collapse
Affiliation(s)
- Josefine H Andersen
- DTU Chemistry, Technical University of Denmark, Kemitorvet Bldg 207, DK-2800 Kongens Lyngby, Denmark
| | - Kaushik D Nanda
- Department of Chemistry, University of Southern California, Los Angeles, California 90089, United States
| | - Anna I Krylov
- Department of Chemistry, University of Southern California, Los Angeles, California 90089, United States
| | - Sonia Coriani
- DTU Chemistry, Technical University of Denmark, Kemitorvet Bldg 207, DK-2800 Kongens Lyngby, Denmark
| |
Collapse
|
25
|
Arias-Martinez JE, Cunha LA, Oosterbaan KJ, Lee J, Head-Gordon M. Accurate core excitation and ionization energies from a state-specific coupled-cluster singles and doubles approach. Phys Chem Chem Phys 2022; 24:20728-20741. [PMID: 36004629 DOI: 10.1039/d2cp01998a] [Citation(s) in RCA: 5] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
We investigate the use of orbital-optimized references in conjunction with single-reference coupled-cluster theory with single and double substitutions (CCSD) for the study of core excitations and ionizations of 18 small organic molecules, without the use of response theory or equation-of-motion (EOM) formalisms. Three schemes are employed to successfully address the convergence difficulties associated with the coupled-cluster equations, and the spin contamination resulting from the use of a spin symmetry-broken reference, in the case of excitations. In order to gauge the inherent potential of the methods studied, an effort is made to provide reasonable basis set limit estimates for the transition energies. Overall, we find that the two best-performing schemes studied here for ΔCCSD are capable of predicting excitation and ionization energies with errors comparable to experimental accuracies. The proposed ΔCCSD schemes reduces statistical errors against experimental excitation energies by more than a factor of two when compared to the frozen-core core-valence separated (FC-CVS) EOM-CCSD approach - a successful variant of EOM-CCSD tailored towards core excitations.
Collapse
Affiliation(s)
- Juan E Arias-Martinez
- Kenneth S. Pitzer Center for Theoretical Chemistry, Department of Chemistry, University of California, Berkeley, California 94720, USA. .,Chemical Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, USA
| | - Leonardo A Cunha
- Kenneth S. Pitzer Center for Theoretical Chemistry, Department of Chemistry, University of California, Berkeley, California 94720, USA. .,Chemical Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, USA
| | - Katherine J Oosterbaan
- Kenneth S. Pitzer Center for Theoretical Chemistry, Department of Chemistry, University of California, Berkeley, California 94720, USA.
| | - Joonho Lee
- Department of Chemistry, Columbia University, New York 10027, USA
| | - Martin Head-Gordon
- Kenneth S. Pitzer Center for Theoretical Chemistry, Department of Chemistry, University of California, Berkeley, California 94720, USA. .,Chemical Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, USA
| |
Collapse
|
26
|
Matz F, Jagau TC. Channel-specific core-valence projectors for determining partial Auger decay widths. Mol Phys 2022. [DOI: 10.1080/00268976.2022.2105270] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/15/2022]
Affiliation(s)
- Florian Matz
- Division of Quantum Chemistry and Physical Chemistry, KU Leuven, Leuven, Belgium
| | | |
Collapse
|
27
|
Opoku RA, Toubin C, Gomes ASP. Simulating core electron binding energies of halogenated species adsorbed on ice surfaces and in solution via relativistic quantum embedding calculations. Phys Chem Chem Phys 2022; 24:14390-14407. [PMID: 35647703 DOI: 10.1039/d1cp05836c] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
In this work, we investigate the effects of the environment on the X-ray photoelectron spectra of hydrogen chloride and chloride ions adsorbed on ice surfaces, as well as of chloride ions in water droplets. In our approach, we combine a density functional theory (DFT) description of the ice surface with that of halogen species using the recently developed relativistic core-valence separation equation of motion coupled cluster (CVS-EOM-IP-CCSD) via the frozen density embedding formalism (FDE), to determine the K and L1,2,3 edges of chlorine. Our calculations, which incorporate temperature effects through snapshots from classical molecular dynamics simulations, are shown to reproduce experimental trends in the change of the core binding energies of Cl- upon moving from a liquid (water droplets) to an interfacial (ice quasi-liquid layer) environment. Our simulations yield water valence band binding energies in good agreement with experiment, which vary little between the droplets and the ice surface. For halide core binding energies there is an overall trend for overestimating experimental values, though good agreement between theory and experiment is found for Cl- in water droplets and on ice. For HCl on the other hand there are significant discrepancies between experimental and calculated core binding energies when we consider structural models that maintain the H-Cl bond more or less intact. An analysis of models that allow for pre-dissociated and dissociated structures suggests that experimentally observed chemical shifts in binding energies between Cl- and HCl would require that H+ (in the form of H3O+) and Cl- are separated by roughly 4-6 Å.
Collapse
Affiliation(s)
- Richard A Opoku
- Univ. Lille, CNRS, UMR 8523 - PhLAM - Physique des Lasers Atomes et Molécules, F-59000 Lille, France.
| | - Céline Toubin
- Univ. Lille, CNRS, UMR 8523 - PhLAM - Physique des Lasers Atomes et Molécules, F-59000 Lille, France.
| | | |
Collapse
|
28
|
Zheng X, Zhang C, Jin Z, Southworth SH, Cheng L. Benchmark relativistic delta-coupled-cluster calculations of K-edge core-ionization energies of third-row elements. Phys Chem Chem Phys 2022; 24:13587-13596. [PMID: 35616685 DOI: 10.1039/d2cp00993e] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/22/2022]
Abstract
A benchmark computational study of K-edge core-ionization energies of third-row elements using relativistic delta-coupled-cluster (ΔCC) methods and a revised core-valence separation (CVS) scheme is reported. High-level relativistic (HLR) corrections beyond the spin-free exact two-component theory in its one-electron variant (SFX2C-1e), including the contributions from two-electron picture-change effects, spin-orbit coupling, the Breit term, and quantum electrodynamics effects, have been taken into account and demonstrated to play an important role. Relativistic ΔCC calculations are shown to provide accurate results for core-ionization energies of third-row elements. The SFX2C-1e-CVS-ΔCC results augmented with HLR corrections show a maximum deviation of less than 0.5 eV with respect to experimental values.
Collapse
Affiliation(s)
- Xuechen Zheng
- Department of Chemistry, The Johns Hopkins University, Baltimore, MD 21218, USA.
| | - Chaoqun Zhang
- Department of Chemistry, The Johns Hopkins University, Baltimore, MD 21218, USA.
| | - Zheqi Jin
- Department of Chemistry, University College London, London, WC1E 6BT, UK
| | - Stephen H Southworth
- Chemical Sciences and Engineering Division, Argonne National Laboratory, Lemont, Illinois 60439, USA
| | - Lan Cheng
- Department of Chemistry, The Johns Hopkins University, Baltimore, MD 21218, USA.
| |
Collapse
|
29
|
Dreuw A, Fransson T. Using core-hole reference states for calculating X-ray photoelectron and emission spectra. Phys Chem Chem Phys 2022; 24:11259-11267. [PMID: 35481477 DOI: 10.1039/d2cp00584k] [Citation(s) in RCA: 8] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
For the calculation of core-ionization energies (IEs), X-ray photoelectron spectra (XPS), and X-ray emission spectra (XES), a commonly applied approach is to use non-Aufbau reference states with a core-hole as either final (IE and XPS) or initial (XES) state. However, such reference states can introduce numerical instabilities in post-HF methods, relating to the denominator of the energy corrections involved. This may become arbitrarily close to zero if a negative virtual orbital is present, e.g. a core-hole, leading to near-singularities. The resulting instabilities lead to severe convergence issues of the calculation schemes and, in addition, can strongly affect both energies and intensities, with oscillator strengths seen to reach values up to 4 × 107. For the K-edge we propose freezing the highest-energy virtual orbitals which contribute to any denominator below a threshold of 0.1 Hartree. Stable and reliable spectra are then produced, with minimal influence due to freezing energetically high-lying virtual orbitals (typically removing <5% of the total number of MOs). The developed protocol is here tested for Møller-Plesset perturbation theory and for the algebraic diagrammatic construction scheme for the polarization propagator, and it is also relevant for coupled cluster theory and other related methods.
Collapse
Affiliation(s)
- Andreas Dreuw
- Interdisciplinary Center for Scientific Computing, Heidelberg University, Im Neuenheimer Feld 205, Heidelberg 69120, Germany
| | - Thomas Fransson
- Department of Theoretical Chemistry and Biology, KTH Royal Institute of Technology, Stockholm 10691, Sweden.
| |
Collapse
|
30
|
Jagau TC. Theory of electronic resonances: fundamental aspects and recent advances. Chem Commun (Camb) 2022; 58:5205-5224. [PMID: 35395664 DOI: 10.1039/d1cc07090h] [Citation(s) in RCA: 17] [Impact Index Per Article: 8.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/21/2023]
Abstract
Electronic resonances are states that are unstable towards loss of electrons. They play critical roles in high-energy environments across chemistry, physics, and biology but are also relevant to processes under ambient conditions that involve unbound electrons. This feature article focuses on complex-variable techniques such as complex scaling and complex absorbing potentials that afford a treatment of electronic resonances in terms of discrete square-integrable eigenstates of non-Hermitian Hamiltonians with complex energy. Fundamental aspects of these techniques as well as their integration into molecular electronic-structure theory are discussed and an overview of some recent developments is given: analytic gradient theory for electronic resonances, the application of rank-reduction techniques and quantum embedding to them, as well as approaches for evaluating partial decay widths.
Collapse
Affiliation(s)
- Thomas-C Jagau
- Department of Chemistry, KU Leuven, Celestijnenlaan 200F, B-3001 Leuven, Belgium.
| |
Collapse
|
31
|
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
|
32
|
Matz F, Jagau TC. Molecular Auger Decay Rates from Complex-Variable Coupled-Cluster Theory. J Chem Phys 2022; 156:114117. [DOI: 10.1063/5.0075646] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
The emission of an Auger electron is the predominant relaxation mechanism of core-vacant states in molecules composed of light nuclei. In this non-radiative decay process, one valence electron fills the core vacancy while a second valence electron is emitted into the ionization continuum. Because of this coupling to the continuum, core-vacant states represent electronic resonances that can be tackled with standard quantum-chemical methods only if they are approximated as bound states, meaning that Auger decay is neglected. Here, we present an approach to compute Auger decay rates of core-vacant states from coupled-cluster and equation-of-motion coupled-cluster wave functions combined with complex scaling of the Hamiltonian or, alternatively, complex-scaled basis functions. Through energy decomposition analysis, we illustrate how complex-scaled methods are capable of describing the coupling to the ionization continuum without the need to model the wave function of the Auger electron explicitly. In addition, we introduce in this work several approaches for the determination of partial decay widths and Auger branching ratios from complex-scaled coupled-cluster wave functions. We demonstrate the capabilities of our new approach by computations on core-ionized states of neon, water, dinitrogen, and benzene. Coupled-cluster and equation-of-motion coupled-cluster theory in the singles and doubles approximation both deliver excellent results for total decay widths, whereas we find partial widths more straightforward to evaluate with the former method. We also observe that the requirements towards the basis set are less arduous for Auger decay than for other types of resonances so that extensions to larger molecules are readily possible.
Collapse
Affiliation(s)
- Florian Matz
- Katholieke Universiteit Leuven Departement Chemie, Belgium
| | | |
Collapse
|
33
|
de Moura CEV, Sokolov AY. Simulating X-ray photoelectron spectra with strong electron correlation using multireference algebraic diagrammatic construction theory. Phys Chem Chem Phys 2022; 24:4769-4784. [DOI: 10.1039/d1cp05476g] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/01/2023]
Abstract
A new theoretical approach for the simulations of X-ray photoelectron spectra of strongly correlated molecular systems that combines multireference algebraic diagrammatic construction theory (MR-ADC) with a core–valence separation (CVS) technique.
Collapse
Affiliation(s)
- Carlos E. V. de Moura
- Department of Chemistry and Biochemistry, The Ohio State University, Columbus, Ohio, 43210, USA
| | - Alexander Yu. Sokolov
- Department of Chemistry and Biochemistry, The Ohio State University, Columbus, Ohio, 43210, USA
| |
Collapse
|
34
|
Kahk JM, Lischner J. Predicting Core Electron Binding Energies in Elements of the First Transition Series Using the Δ-Self-Consistent-Field Method. Faraday Discuss 2022; 236:364-373. [DOI: 10.1039/d1fd00103e] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
The Δ-Self-Consistent-Field (ΔSCF) method has been established as an accurate and computationally efficient approach for calculating absolute core electron binding energies for light elements up to chlorine, but relatively little...
Collapse
|
35
|
Mejia-Rodriguez D, Kunitsa A, Aprà E, Govind N. Scalable Molecular GW Calculations: Valence and Core Spectra. J Chem Theory Comput 2021; 17:7504-7517. [PMID: 34855381 DOI: 10.1021/acs.jctc.1c00738] [Citation(s) in RCA: 14] [Impact Index Per Article: 4.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
We present a scalable implementation of the GW approximation using Gaussian atomic orbitals to study the valence and core ionization spectroscopies of molecules. The implementation of the standard spectral decomposition approach to the screened-Coulomb interaction, as well as a contour-deformation method, is described. We have implemented both of these approaches using the robust variational fitting approximation to the four-center electron repulsion integrals. We have utilized the MINRES solver with the contour-deformation approach to reduce the computational scaling by 1 order of magnitude. A complex heuristic in the quasiparticle equation solver further allows a speed-up of the computation of core and semicore ionization energies. Benchmark tests using the GW100 and CORE65 data sets and the carbon 1s binding energy of the well-studied ethyl trifluoroacetate, or ESCA molecule, were performed to validate the accuracy of our implementation. We also demonstrate and discuss the parallel performance and computational scaling of our implementation using a range of water clusters of increasing size.
Collapse
Affiliation(s)
- Daniel Mejia-Rodriguez
- Environmental Molecular Sciences Laboratory, Pacific Northwest National Laboratory, Richland, Washington 99352, United States
| | - Alexander Kunitsa
- Department of Chemistry, University of Illinois at Urbana-Champaign, 600 South Mathews Avenue, Urbana, Illinois 61801, United States
| | - Edoardo Aprà
- Environmental Molecular Sciences Laboratory, Pacific Northwest National Laboratory, Richland, Washington 99352, United States
| | - Niranjan Govind
- Physical and Computational Sciences Directorate, Pacific Northwest National Laboratory, Richland, Washington 99352, United States
| |
Collapse
|
36
|
Kahk JM, Michelitsch GS, Maurer RJ, Reuter K, Lischner J. Core Electron Binding Energies in Solids from Periodic All-Electron Δ-Self-Consistent-Field Calculations. J Phys Chem Lett 2021; 12:9353-9359. [PMID: 34549969 DOI: 10.1021/acs.jpclett.1c02380] [Citation(s) in RCA: 9] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/13/2023]
Abstract
Theoretical calculations of core electron binding energies are required for the interpretation of experimental X-ray photoelectron spectra, but achieving accurate results for solids has proven difficult. In this work, we demonstrate that accurate absolute core electron binding energies in both metallic and insulating solids can be obtained from periodic all-electron Δ-self-consistent-field (ΔSCF) calculations. In particular, we show that core electron binding energies referenced to the valence band maximum can be obtained as total energy differences between two (N - 1)-electron systems: one with a core hole and one with an electron removed from the highest occupied valence state. To achieve convergence with respect to the supercell size, the analogy between localized core holes and charged defects is exploited. Excellent agreement between calculated and experimental core electron binding energies is found for both metals and insulators, with a mean absolute error of 0.24 eV for the systems considered.
Collapse
Affiliation(s)
- J Matthias Kahk
- Department of Materials, Imperial College London, South Kensington, London SW7 2AZ, United Kingdom
- Institute of Physics, University of Tartu, W. Ostwaldi 1, 50411 Tartu, Estonia
| | - Georg S Michelitsch
- Chair for Theoretical Chemistry, Technische Universität München, Lichtenbergstr. 4, D-85747 Garching, Germany
| | - Reinhard J Maurer
- Chair for Theoretical Chemistry, Technische Universität München, Lichtenbergstr. 4, D-85747 Garching, Germany
- Department of Chemistry, University of Warwick, Gibbet Hill Rd., Coventry CV4 7AL, United Kingdom
| | - Karsten Reuter
- Chair for Theoretical Chemistry, Technische Universität München, Lichtenbergstr. 4, D-85747 Garching, Germany
| | - Johannes Lischner
- Department of Physics and Department of Materials, and the Thomas Young Centre for Theory and Simulation of Materials, Imperial College London, London SW7 2AZ, United Kingdom
| |
Collapse
|
37
|
Epifanovsky E, Gilbert ATB, Feng X, Lee J, Mao Y, Mardirossian N, Pokhilko P, White AF, Coons MP, Dempwolff AL, Gan Z, Hait D, Horn PR, Jacobson LD, Kaliman I, Kussmann J, Lange AW, Lao KU, Levine DS, Liu J, McKenzie SC, Morrison AF, Nanda KD, Plasser F, Rehn DR, Vidal ML, You ZQ, Zhu Y, Alam B, Albrecht BJ, Aldossary A, Alguire E, Andersen JH, Athavale V, Barton D, Begam K, Behn A, Bellonzi N, Bernard YA, Berquist EJ, Burton HGA, Carreras A, Carter-Fenk K, Chakraborty R, Chien AD, Closser KD, Cofer-Shabica V, Dasgupta S, de Wergifosse M, Deng J, Diedenhofen M, Do H, Ehlert S, Fang PT, Fatehi S, Feng Q, Friedhoff T, Gayvert J, Ge Q, Gidofalvi G, Goldey M, Gomes J, González-Espinoza CE, Gulania S, Gunina AO, Hanson-Heine MWD, Harbach PHP, Hauser A, Herbst MF, Hernández Vera M, Hodecker M, Holden ZC, Houck S, Huang X, Hui K, Huynh BC, Ivanov M, Jász Á, Ji H, Jiang H, Kaduk B, Kähler S, Khistyaev K, Kim J, Kis G, Klunzinger P, Koczor-Benda Z, Koh JH, Kosenkov D, Koulias L, Kowalczyk T, Krauter CM, Kue K, Kunitsa A, Kus T, Ladjánszki I, Landau A, Lawler KV, Lefrancois D, Lehtola S, Li RR, Li YP, Liang J, Liebenthal M, Lin HH, Lin YS, Liu F, Liu KY, Loipersberger M, Luenser A, Manjanath A, Manohar P, Mansoor E, Manzer SF, Mao SP, Marenich AV, Markovich T, Mason S, Maurer SA, McLaughlin PF, Menger MFSJ, Mewes JM, Mewes SA, Morgante P, Mullinax JW, Oosterbaan KJ, Paran G, Paul AC, Paul SK, Pavošević F, Pei Z, Prager S, Proynov EI, Rák Á, Ramos-Cordoba E, Rana B, Rask AE, Rettig A, Richard RM, Rob F, Rossomme E, Scheele T, Scheurer M, Schneider M, Sergueev N, Sharada SM, Skomorowski W, Small DW, Stein CJ, Su YC, Sundstrom EJ, Tao Z, Thirman J, Tornai GJ, Tsuchimochi T, Tubman NM, Veccham SP, Vydrov O, Wenzel J, Witte J, Yamada A, Yao K, Yeganeh S, Yost SR, Zech A, Zhang IY, Zhang X, Zhang Y, Zuev D, Aspuru-Guzik A, Bell AT, Besley NA, Bravaya KB, Brooks BR, Casanova D, Chai JD, Coriani S, Cramer CJ, Cserey G, DePrince AE, DiStasio RA, Dreuw A, Dunietz BD, Furlani TR, Goddard WA, Hammes-Schiffer S, Head-Gordon T, Hehre WJ, Hsu CP, Jagau TC, Jung Y, Klamt A, Kong J, Lambrecht DS, Liang W, Mayhall NJ, McCurdy CW, Neaton JB, Ochsenfeld C, Parkhill JA, Peverati R, Rassolov VA, Shao Y, Slipchenko LV, Stauch T, Steele RP, Subotnik JE, Thom AJW, Tkatchenko A, Truhlar DG, Van Voorhis T, Wesolowski TA, Whaley KB, Woodcock HL, Zimmerman PM, Faraji S, Gill PMW, Head-Gordon M, Herbert JM, Krylov AI. Software for the frontiers of quantum chemistry: An overview of developments in the Q-Chem 5 package. J Chem Phys 2021; 155:084801. [PMID: 34470363 PMCID: PMC9984241 DOI: 10.1063/5.0055522] [Citation(s) in RCA: 480] [Impact Index Per Article: 160.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/04/2023] Open
Abstract
This article summarizes technical advances contained in the fifth major release of the Q-Chem quantum chemistry program package, covering developments since 2015. A comprehensive library of exchange-correlation functionals, along with a suite of correlated many-body methods, continues to be a hallmark of the Q-Chem software. The many-body methods include novel variants of both coupled-cluster and configuration-interaction approaches along with methods based on the algebraic diagrammatic construction and variational reduced density-matrix methods. Methods highlighted in Q-Chem 5 include a suite of tools for modeling core-level spectroscopy, methods for describing metastable resonances, methods for computing vibronic spectra, the nuclear-electronic orbital method, and several different energy decomposition analysis techniques. High-performance capabilities including multithreaded parallelism and support for calculations on graphics processing units are described. Q-Chem boasts a community of well over 100 active academic developers, and the continuing evolution of the software is supported by an "open teamware" model and an increasingly modular design.
Collapse
Affiliation(s)
- Evgeny Epifanovsky
- Q-Chem, Inc., 6601 Owens Drive, Suite 105, Pleasanton, California 94588, USA
| | | | | | - Joonho Lee
- Department of Chemistry, University of California, Berkeley, California 94720, USA
| | - Yuezhi Mao
- Department of Chemistry, University of California, Berkeley, California 94720, USA
| | | | - Pavel Pokhilko
- Department of Chemistry, University of Southern California, Los Angeles, California 90089, USA
| | - Alec F. White
- Department of Chemistry, University of California, Berkeley, California 94720, USA
| | - Marc P. Coons
- Department of Chemistry and Biochemistry, The Ohio State University, Columbus, Ohio 43210, USA
| | - Adrian L. Dempwolff
- Interdisciplinary Center for Scientific Computing, Ruprecht-Karls University, Im Neuenheimer Feld 205, 69120 Heidelberg, Germany
| | - Zhengting Gan
- Q-Chem, Inc., 6601 Owens Drive, Suite 105, Pleasanton, California 94588, USA
| | - Diptarka Hait
- Department of Chemistry, University of California, Berkeley, California 94720, USA
| | - Paul R. Horn
- Department of Chemistry, University of California, Berkeley, California 94720, USA
| | - Leif D. Jacobson
- Department of Chemistry and Biochemistry, The Ohio State University, Columbus, Ohio 43210, USA
| | | | - Jörg Kussmann
- Department of Chemistry, Ludwig Maximilian University, Butenandtstr. 7, D-81377 München, Germany
| | - Adrian W. Lange
- Department of Chemistry and Biochemistry, The Ohio State University, Columbus, Ohio 43210, USA
| | - Ka Un Lao
- Department of Chemistry and Biochemistry, The Ohio State University, Columbus, Ohio 43210, USA
| | - Daniel S. Levine
- Department of Chemistry, University of California, Berkeley, California 94720, USA
| | | | - Simon C. McKenzie
- Research School of Chemistry, Australian National University, Canberra, Australia
| | | | - Kaushik D. Nanda
- Department of Chemistry, University of Southern California, Los Angeles, California 90089, USA
| | | | - Dirk R. Rehn
- Interdisciplinary Center for Scientific Computing, Ruprecht-Karls University, Im Neuenheimer Feld 205, 69120 Heidelberg, Germany
| | - Marta L. Vidal
- Department of Chemistry, Technical University of Denmark, Kemitorvet Bldg. 207, DK-2800 Kgs Lyngby, Denmark
| | | | - Ying Zhu
- Department of Chemistry and Biochemistry, The Ohio State University, Columbus, Ohio 43210, USA
| | - Bushra Alam
- Department of Chemistry and Biochemistry, The Ohio State University, Columbus, Ohio 43210, USA
| | - Benjamin J. Albrecht
- Department of Chemistry, University of Pittsburgh, Pittsburgh, Pennsylvania 15260, USA
| | | | - Ethan Alguire
- Department of Chemistry, University of Pennsylvania, Philadelphia, Pennsylvania 19104, USA
| | - Josefine H. Andersen
- Department of Chemistry, Technical University of Denmark, Kemitorvet Bldg. 207, DK-2800 Kgs Lyngby, Denmark
| | - Vishikh Athavale
- Department of Chemistry, University of Pennsylvania, Philadelphia, Pennsylvania 19104, USA
| | - Dennis Barton
- Department of Physics and Materials Science, University of Luxembourg, L-1511 Luxembourg, Luxembourg
| | - Khadiza Begam
- Department of Physics, Kent State University, Kent, Ohio 44242, USA
| | - Andrew Behn
- Department of Chemistry, University of California, Berkeley, California 94720, USA
| | - Nicole Bellonzi
- Department of Chemistry, University of Pennsylvania, Philadelphia, Pennsylvania 19104, USA
| | - Yves A. Bernard
- Department of Chemistry, University of Southern California, Los Angeles, California 90089, USA
| | | | - Hugh G. A. Burton
- Department of Chemistry, University of Cambridge, Cambridge, United Kingdom
| | - Abel Carreras
- Donostia International Physics Center, 20080 Donostia, Euskadi, Spain
| | - Kevin Carter-Fenk
- Department of Chemistry and Biochemistry, The Ohio State University, Columbus, Ohio 43210, USA
| | | | - Alan D. Chien
- Department of Chemistry, University of Michigan, Ann Arbor, Michigan 48109, USA
| | | | - Vale Cofer-Shabica
- Department of Chemistry, University of Pennsylvania, Philadelphia, Pennsylvania 19104, USA
| | - Saswata Dasgupta
- Department of Chemistry and Biochemistry, The Ohio State University, Columbus, Ohio 43210, USA
| | - Marc de Wergifosse
- Department of Chemistry, University of Southern California, Los Angeles, California 90089, USA
| | - Jia Deng
- Research School of Chemistry, Australian National University, Canberra, Australia
| | | | - Hainam Do
- School of Chemistry, University of Nottingham, Nottingham, United Kingdom
| | - Sebastian Ehlert
- Mulliken Center for Theoretical Chemistry, Institut für Physikalische und Theoretische Chemie, Beringstr. 4, 53115 Bonn, Germany
| | - Po-Tung Fang
- Department of Physics, National Taiwan University, Taipei 10617, Taiwan
| | | | - Qingguo Feng
- Department of Chemistry and Biochemistry, Kent State University, Kent, Ohio 44240, USA
| | - Triet Friedhoff
- Department of Chemistry and Biochemistry, University of Notre Dame, Notre Dame, Indiana 46556, USA
| | - James Gayvert
- Department of Chemistry, Boston University, Boston, Massachusetts 02215, USA
| | - Qinghui Ge
- Department of Chemistry, University of California, Berkeley, California 94720, USA
| | - Gergely Gidofalvi
- Department of Chemistry and Biochemistry, Gonzaga University, Spokane, Washington 99258, USA
| | - Matthew Goldey
- Department of Chemistry, University of California, Berkeley, California 94720, USA
| | - Joe Gomes
- Department of Chemistry, University of California, Berkeley, California 94720, USA
| | | | - Sahil Gulania
- Department of Chemistry, University of Southern California, Los Angeles, California 90089, USA
| | - Anastasia O. Gunina
- Department of Chemistry, University of Southern California, Los Angeles, California 90089, USA
| | | | - Phillip H. P. Harbach
- Interdisciplinary Center for Scientific Computing, Ruprecht-Karls University, Im Neuenheimer Feld 205, 69120 Heidelberg, Germany
| | - Andreas Hauser
- Institute of Experimental Physics, Graz University of Technology, Graz, Austria
| | | | - Mario Hernández Vera
- Department of Chemistry, Ludwig Maximilian University, Butenandtstr. 7, D-81377 München, Germany
| | - Manuel Hodecker
- Interdisciplinary Center for Scientific Computing, Ruprecht-Karls University, Im Neuenheimer Feld 205, 69120 Heidelberg, Germany
| | - Zachary C. Holden
- Department of Chemistry and Biochemistry, The Ohio State University, Columbus, Ohio 43210, USA
| | - Shannon Houck
- Department of Chemistry, Virginia Tech, Blacksburg, Virginia 24061, USA
| | - Xunkun Huang
- Department of Chemistry, Xiamen University, Xiamen 361005, China
| | - Kerwin Hui
- Department of Physics, National Taiwan University, Taipei 10617, Taiwan
| | - Bang C. Huynh
- Department of Chemistry, University of Cambridge, Cambridge, United Kingdom
| | - Maxim Ivanov
- Department of Chemistry, University of Southern California, Los Angeles, California 90089, USA
| | - Ádám Jász
- Stream Novation Ltd., Práter utca 50/a, H-1083 Budapest, Hungary
| | - Hyunjun Ji
- Graduate School of Energy, Environment, Water and Sustainability (EEWS), Korea Advanced Institute of Science and Technology (KAIST), Daejeon 34141, Republic of Korea
| | - Hanjie Jiang
- Department of Chemistry, University of Michigan, Ann Arbor, Michigan 48109, USA
| | - Benjamin Kaduk
- Department of Chemistry, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, USA
| | - Sven Kähler
- Department of Chemistry, University of Southern California, Los Angeles, California 90089, USA
| | - Kirill Khistyaev
- Department of Chemistry, University of Southern California, Los Angeles, California 90089, USA
| | - Jaehoon Kim
- Graduate School of Energy, Environment, Water and Sustainability (EEWS), Korea Advanced Institute of Science and Technology (KAIST), Daejeon 34141, Republic of Korea
| | - Gergely Kis
- Stream Novation Ltd., Práter utca 50/a, H-1083 Budapest, Hungary
| | | | - Zsuzsanna Koczor-Benda
- Department of Chemistry, Ludwig Maximilian University, Butenandtstr. 7, D-81377 München, Germany
| | - Joong Hoon Koh
- Department of Chemistry and Biochemistry, University of Notre Dame, Notre Dame, Indiana 46556, USA
| | - Dimitri Kosenkov
- Department of Chemistry, Purdue University, West Lafayette, Indiana 47907, USA
| | - Laura Koulias
- Department of Chemistry and Biochemistry, Florida State University, Tallahassee, Florida 32306, USA
| | | | - Caroline M. Krauter
- Interdisciplinary Center for Scientific Computing, Ruprecht-Karls University, Im Neuenheimer Feld 205, 69120 Heidelberg, Germany
| | - Karl Kue
- Institute of Chemistry, Academia Sinica, 128, Academia Road Section 2, Nangang District, Taipei 11529, Taiwan
| | - Alexander Kunitsa
- Department of Chemistry, Boston University, Boston, Massachusetts 02215, USA
| | - Thomas Kus
- Department of Chemistry, University of Southern California, Los Angeles, California 90089, USA
| | | | - Arie Landau
- Department of Chemistry, University of Southern California, Los Angeles, California 90089, USA
| | - Keith V. Lawler
- Department of Chemistry, University of California, Berkeley, California 94720, USA
| | - Daniel Lefrancois
- Interdisciplinary Center for Scientific Computing, Ruprecht-Karls University, Im Neuenheimer Feld 205, 69120 Heidelberg, Germany
| | | | - Run R. Li
- Department of Chemistry and Biochemistry, Florida State University, Tallahassee, Florida 32306, USA
| | - Yi-Pei Li
- Department of Chemistry, University of California, Berkeley, California 94720, USA
| | - Jiashu Liang
- Department of Chemistry, University of California, Berkeley, California 94720, USA
| | - Marcus Liebenthal
- Department of Chemistry and Biochemistry, Florida State University, Tallahassee, Florida 32306, USA
| | - Hung-Hsuan Lin
- Institute of Chemistry, Academia Sinica, 128, Academia Road Section 2, Nangang District, Taipei 11529, Taiwan
| | - You-Sheng Lin
- Department of Physics, National Taiwan University, Taipei 10617, Taiwan
| | - Fenglai Liu
- Q-Chem, Inc., 6601 Owens Drive, Suite 105, Pleasanton, California 94588, USA
| | | | | | - Arne Luenser
- Department of Chemistry, Ludwig Maximilian University, Butenandtstr. 7, D-81377 München, Germany
| | - Aaditya Manjanath
- Institute of Chemistry, Academia Sinica, 128, Academia Road Section 2, Nangang District, Taipei 11529, Taiwan
| | - Prashant Manohar
- Department of Chemistry, University of Southern California, Los Angeles, California 90089, USA
| | - Erum Mansoor
- Department of Chemistry, University of California, Berkeley, California 94720, USA
| | - Sam F. Manzer
- Department of Chemistry, University of California, Berkeley, California 94720, USA
| | - Shan-Ping Mao
- Department of Physics, National Taiwan University, Taipei 10617, Taiwan
| | | | - Thomas Markovich
- Department of Chemistry and Chemical Biology, Harvard University, Cambridge, Massachusetts 02138, USA
| | - Stephen Mason
- School of Chemistry, University of Nottingham, Nottingham, United Kingdom
| | - Simon A. Maurer
- Department of Chemistry, Ludwig Maximilian University, Butenandtstr. 7, D-81377 München, Germany
| | - Peter F. McLaughlin
- Q-Chem, Inc., 6601 Owens Drive, Suite 105, Pleasanton, California 94588, USA
| | | | - Jan-Michael Mewes
- Interdisciplinary Center for Scientific Computing, Ruprecht-Karls University, Im Neuenheimer Feld 205, 69120 Heidelberg, Germany
| | - Stefanie A. Mewes
- Interdisciplinary Center for Scientific Computing, Ruprecht-Karls University, Im Neuenheimer Feld 205, 69120 Heidelberg, Germany
| | - Pierpaolo Morgante
- Department of Chemistry, Florida Institute of Technology, Melbourne, Florida 32901, USA
| | - J. Wayne Mullinax
- Department of Chemistry, Florida Institute of Technology, Melbourne, Florida 32901, USA
| | | | | | - Alexander C. Paul
- Interdisciplinary Center for Scientific Computing, Ruprecht-Karls University, Im Neuenheimer Feld 205, 69120 Heidelberg, Germany
| | - Suranjan K. Paul
- Department of Chemistry and Biochemistry, The Ohio State University, Columbus, Ohio 43210, USA
| | - Fabijan Pavošević
- Department of Chemistry, Yale University, New Haven, Connecticut 06520, USA
| | - Zheng Pei
- School of Electrical and Computer Engineering, University of Oklahoma, Norman, Oklahoma 73019, USA
| | - Stefan Prager
- Interdisciplinary Center for Scientific Computing, Ruprecht-Karls University, Im Neuenheimer Feld 205, 69120 Heidelberg, Germany
| | - Emil I. Proynov
- Q-Chem, Inc., 6601 Owens Drive, Suite 105, Pleasanton, California 94588, USA
| | - Ádám Rák
- Stream Novation Ltd., Práter utca 50/a, H-1083 Budapest, Hungary
| | - Eloy Ramos-Cordoba
- Department of Chemistry, University of California, Berkeley, California 94720, USA
| | - Bhaskar Rana
- Department of Chemistry and Biochemistry, The Ohio State University, Columbus, Ohio 43210, USA
| | - Alan E. Rask
- Department of Chemistry, University of Michigan, Ann Arbor, Michigan 48109, USA
| | - Adam Rettig
- Department of Chemistry, University of California, Berkeley, California 94720, USA
| | - Ryan M. Richard
- Department of Chemistry and Biochemistry, The Ohio State University, Columbus, Ohio 43210, USA
| | - Fazle Rob
- Q-Chem, Inc., 6601 Owens Drive, Suite 105, Pleasanton, California 94588, USA
| | - Elliot Rossomme
- Department of Chemistry, University of California, Berkeley, California 94720, USA
| | - Tarek Scheele
- Institute for Physical and Theoretical Chemistry, University of Bremen, Bremen, Germany
| | - Maximilian Scheurer
- Interdisciplinary Center for Scientific Computing, Ruprecht-Karls University, Im Neuenheimer Feld 205, 69120 Heidelberg, Germany
| | - Matthias Schneider
- Interdisciplinary Center for Scientific Computing, Ruprecht-Karls University, Im Neuenheimer Feld 205, 69120 Heidelberg, Germany
| | - Nickolai Sergueev
- Department of Chemistry and Biochemistry, Kent State University, Kent, Ohio 44240, USA
| | - Shaama M. Sharada
- Department of Chemistry, University of California, Berkeley, California 94720, USA
| | - Wojciech Skomorowski
- Department of Chemistry, University of Southern California, Los Angeles, California 90089, USA
| | - David W. Small
- Department of Chemistry, University of California, Berkeley, California 94720, USA
| | - Christopher J. Stein
- Department of Chemistry, University of California, Berkeley, California 94720, USA
| | - Yu-Chuan Su
- Department of Physics, National Taiwan University, Taipei 10617, Taiwan
| | - Eric J. Sundstrom
- Department of Chemistry, University of California, Berkeley, California 94720, USA
| | - Zhen Tao
- Department of Chemistry, Yale University, New Haven, Connecticut 06520, USA
| | - Jonathan Thirman
- Department of Chemistry, University of California, Berkeley, California 94720, USA
| | - Gábor J. Tornai
- Stream Novation Ltd., Práter utca 50/a, H-1083 Budapest, Hungary
| | - Takashi Tsuchimochi
- Department of Chemistry, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, USA
| | - Norm M. Tubman
- Department of Chemistry, University of California, Berkeley, California 94720, USA
| | | | - Oleg Vydrov
- Department of Chemistry, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, USA
| | - Jan Wenzel
- Interdisciplinary Center for Scientific Computing, Ruprecht-Karls University, Im Neuenheimer Feld 205, 69120 Heidelberg, Germany
| | - Jon Witte
- Department of Chemistry, University of California, Berkeley, California 94720, USA
| | - Atsushi Yamada
- Department of Chemistry and Biochemistry, Kent State University, Kent, Ohio 44240, USA
| | - Kun Yao
- Department of Chemistry and Biochemistry, University of Notre Dame, Notre Dame, Indiana 46556, USA
| | - Sina Yeganeh
- Department of Chemistry, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, USA
| | - Shane R. Yost
- Department of Chemistry, University of California, Berkeley, California 94720, USA
| | - Alexander Zech
- Department of Physical Chemistry, University of Geneva, 30, Quai Ernest-Ansermet, CH-1211 Geneva 4, Switzerland
| | - Igor Ying Zhang
- Department of Chemistry, Fudan University, Shanghai 200433, China
| | - Xing Zhang
- Department of Chemistry and Biochemistry, The Ohio State University, Columbus, Ohio 43210, USA
| | - Yu Zhang
- Q-Chem, Inc., 6601 Owens Drive, Suite 105, Pleasanton, California 94588, USA
| | - Dmitry Zuev
- Department of Chemistry, University of Southern California, Los Angeles, California 90089, USA
| | - Alán Aspuru-Guzik
- Department of Chemistry and Chemical Biology, Harvard University, Cambridge, Massachusetts 02138, USA
| | - Alexis T. Bell
- Department of Chemical Engineering, University of California, Berkeley, California 94720, USA
| | - Nicholas A. Besley
- School of Chemistry, University of Nottingham, Nottingham, United Kingdom
| | - Ksenia B. Bravaya
- Department of Chemistry, Boston University, Boston, Massachusetts 02215, USA
| | - Bernard R. Brooks
- Laboratory of Computational Biophysics, National Institute of Health, Bethesda, Maryland 20892, USA
| | - David Casanova
- Donostia International Physics Center, 20080 Donostia, Euskadi, Spain
| | | | - Sonia Coriani
- Department of Chemistry, Technical University of Denmark, Kemitorvet Bldg. 207, DK-2800 Kgs Lyngby, Denmark
| | | | | | - A. Eugene DePrince
- Department of Chemistry and Biochemistry, Florida State University, Tallahassee, Florida 32306, USA
| | - Robert A. DiStasio
- Department of Chemistry and Chemical Biology, Cornell University, Ithaca, New York 14853, USA
| | - Andreas Dreuw
- Interdisciplinary Center for Scientific Computing, Ruprecht-Karls University, Im Neuenheimer Feld 205, 69120 Heidelberg, Germany
| | - Barry D. Dunietz
- Department of Chemistry and Biochemistry, Kent State University, Kent, Ohio 44240, USA
| | - Thomas R. Furlani
- Department of Chemistry, University at Buffalo, State University of New York, Buffalo, New York 14260, USA
| | - William A. Goddard
- Materials and Process Simulation Center, California Institute of Technology, Pasadena, California 91125, USA
| | | | - Teresa Head-Gordon
- Department of Chemistry, University of California, Berkeley, California 94720, USA
| | | | | | | | - Yousung Jung
- Graduate School of Energy, Environment, Water and Sustainability (EEWS), Korea Advanced Institute of Science and Technology (KAIST), Daejeon 34141, Republic of Korea
| | - Andreas Klamt
- COSMOlogic GmbH & Co. KG, Imbacher Weg 46, D-51379 Leverkusen, Germany
| | - Jing Kong
- Q-Chem, Inc., 6601 Owens Drive, Suite 105, Pleasanton, California 94588, USA
| | - Daniel S. Lambrecht
- Department of Chemistry, University of Pittsburgh, Pittsburgh, Pennsylvania 15260, USA
| | | | | | - C. William McCurdy
- Department of Chemistry, University of California, Davis, California 95616, USA
| | - Jeffrey B. Neaton
- Department of Physics, University of California, Berkeley, California 94720, USA
| | - Christian Ochsenfeld
- Department of Chemistry, Ludwig Maximilian University, Butenandtstr. 7, D-81377 München, Germany
| | - John A. Parkhill
- Department of Chemistry and Biochemistry, University of Notre Dame, Notre Dame, Indiana 46556, USA
| | - Roberto Peverati
- Department of Chemistry, Florida Institute of Technology, Melbourne, Florida 32901, USA
| | - Vitaly A. Rassolov
- Department of Chemistry and Biochemistry, University of South Carolina, Columbia, South Carolina 29208, USA
| | | | | | | | - Ryan P. Steele
- Department of Chemistry and Henry Eyring Center for Theoretical Chemistry, University of Utah, Salt Lake City, Utah 84112, USA
| | - Joseph E. Subotnik
- Department of Chemistry, University of Pennsylvania, Philadelphia, Pennsylvania 19104, USA
| | - Alex J. W. Thom
- Department of Chemistry, University of Cambridge, Cambridge, United Kingdom
| | - Alexandre Tkatchenko
- Department of Physics and Materials Science, University of Luxembourg, L-1511 Luxembourg, Luxembourg
| | - Donald G. Truhlar
- Department of Chemistry, University of Minnesota, Minneapolis, Minnesota 55455, USA
| | - Troy Van Voorhis
- Department of Chemistry, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, USA
| | - Tomasz A. Wesolowski
- Department of Physical Chemistry, University of Geneva, 30, Quai Ernest-Ansermet, CH-1211 Geneva 4, Switzerland
| | - K. Birgitta Whaley
- Department of Chemistry, University of California, Berkeley, California 94720, USA
| | - H. Lee Woodcock
- Department of Chemistry, University of South Florida, Tampa, Florida 33620, USA
| | - Paul M. Zimmerman
- Department of Chemistry, University of Michigan, Ann Arbor, Michigan 48109, USA
| | - Shirin Faraji
- Zernike Institute for Advanced Materials, University of Groningen, 9774AG Groningen, The Netherlands
| | | | - Martin Head-Gordon
- Department of Chemistry, University of California, Berkeley, California 94720, USA
| | - John M. Herbert
- Department of Chemistry and Biochemistry, The Ohio State University, Columbus, Ohio 43210, USA
| | - Anna I. Krylov
- Department of Chemistry, University of Southern California, Los Angeles, California 90089, USA,Author to whom correspondence should be addressed:
| |
Collapse
|
38
|
Schraivogel T, Kats D. Accuracy of the distinguishable cluster approximation for triple excitations for open-shell molecules and excited states. J Chem Phys 2021; 155:064101. [PMID: 34391360 DOI: 10.1063/5.0059181] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
The distinguishable cluster approximation for triple excitations has been applied to calculate thermochemical properties and excited states involving closed-shell and open-shell species, such as small molecules, 3d transition metal atoms, ozone, and an iron-porphyrin model. Excitation energies have been computed using the ΔCC approach by directly optimizing the excited states. A fixed-reference technique has been introduced to target selected spin-states for open-shell molecular systems. The distinguishable cluster approximation consistently improves coupled cluster with singles doubles and triples results for absolute and relative energies. For excited states dominated by a single configuration state function, the fixed-reference approach combined with high-level coupled-cluster methods has a comparable accuracy to the corresponding equation-of-motion coupled-cluster methods with a negligible amount of spin contamination.
Collapse
Affiliation(s)
- Thomas Schraivogel
- Max Planck Institute for Solid State Research, Heisenbergstraße 1, 70569 Stuttgart, Germany
| | - Daniel Kats
- Max Planck Institute for Solid State Research, Heisenbergstraße 1, 70569 Stuttgart, Germany
| |
Collapse
|
39
|
Zhao R, Grofe A, Wang Z, Bao P, Chen X, Liu W, Gao J. Dynamic-then-Static Approach for Core Excitations of Open-Shell Molecules. J Phys Chem Lett 2021; 12:7409-7417. [PMID: 34328742 DOI: 10.1021/acs.jpclett.1c02039] [Citation(s) in RCA: 19] [Impact Index Per Article: 6.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/13/2023]
Abstract
Delta self-consistent-field methods are widely used in studies of electronically excited states. However, the nonaufbau determinants are generally spin-contaminated. Here, we describe a general approach for spin-coupling interactions of open-shell molecules, making use of multistate density functional theory (MSDFT). In particular, the effective exchange integrals that determine spin coupling are obtained by enforcing the multiplet degeneracy of the S+1 state in the MS = S manifold. Consequently, they are consistent with the energy of the high-spin state that is adequately treated by Kohn-Sham density functional theory (DFT) and, thereby, free of double counting of correlation. The method was applied to core excitations of open-shell molecules and compared with those by spin-adapted time-dependent DFT. An excellent agreement with experiment was found employing the BLYP functional and aug-cc-pCVQZ basis set. Overall, MSDFT provides an effective combination of the strengths of DFT and wave function theory to achieve efficiency and accuracy.
Collapse
Affiliation(s)
- Ruoqi Zhao
- Institute of Systems and Physical Biology, Shenzhen Bay Laboratory, Shenzhen 518055, Guangdong, China
- Institute of Theoretical Chemistry, Jilin University Changchun 130023, Jilin, China
| | - Adam Grofe
- Institute of Theoretical Chemistry, Jilin University Changchun 130023, Jilin, China
| | - Zikuan Wang
- Qingdao Institute for Theoretical and Computational Sciences, Institute of Frontier and Interdisciplinary Science, Shandong University, Qingdao 266237, Shandong, China
| | - Peng Bao
- Beijing National Laboratory for Molecular Sciences, State Key Laboratory for Structural Chemistry of Unstable and Stable Species, Institute of Chemistry, Chinese Academy of Sciences, Beijing 100190, China
| | - Xin Chen
- Institute of Systems and Physical Biology, Shenzhen Bay Laboratory, Shenzhen 518055, Guangdong, China
| | - Wenjian Liu
- Qingdao Institute for Theoretical and Computational Sciences, Institute of Frontier and Interdisciplinary Science, Shandong University, Qingdao 266237, Shandong, China
| | - Jiali Gao
- Institute of Systems and Physical Biology, Shenzhen Bay Laboratory, Shenzhen 518055, Guangdong, China
- Department of Chemistry and Supercomputing Institute, University of Minnesota, Minneapolis 55455, Minnesota, United States
| |
Collapse
|
40
|
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: 27] [Impact Index Per Article: 9.0] [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
|
41
|
Besley NA. Density Functional Theory Calculations of Core-Electron Binding Energies at the K-Edge of Heavier Elements. J Chem Theory Comput 2021; 17:3644-3651. [PMID: 34042434 DOI: 10.1021/acs.jctc.1c00171] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/23/2022]
Abstract
The capability to determine core-electron binding energies (CEBEs) is vital in the analysis of X-ray photoelectron spectroscopy, and the continued development of light sources has made inner shell spectroscopy of heavier elements increasingly accessible. Density functional theory is widely used to determine CEBEs of lighter elements (boron-fluorine). It is shown that good performance of exchange-correlation functionals for these elements does not necessarily translate to the calculation of CEBEs for the heavier elements from the next row of the periodic table, and in general, larger errors are observed. Two strategies are explored that improve the accuracy of the calculated CEBEs. The first is to apply element and functional dependent energy corrections, and the second is a reparametrization of a short-range corrected functional. This functional is able to reproduce experimental phosphorus and sulfur K-edge CEBEs with an average error of 0.15 eV demonstrating the importance of reducing the self-interaction error associated with the core electrons and represents progress toward a density functional theory calculation that performs equally well for ionization at the K-edge of all elements.
Collapse
Affiliation(s)
- Nicholas A Besley
- School of Chemistry, University of Nottingham, University Park, Nottingham NG7 2RD, U.K
| |
Collapse
|
42
|
Klein BP, Hall SJ, Maurer RJ. The nuts and bolts of core-hole constrained ab initiosimulation for K-shell x-ray photoemission and absorption spectra. JOURNAL OF PHYSICS. CONDENSED MATTER : AN INSTITUTE OF PHYSICS JOURNAL 2021; 33. [PMID: 33682682 DOI: 10.1088/1361-648x/abdf00] [Citation(s) in RCA: 16] [Impact Index Per Article: 5.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 10/20/2020] [Accepted: 01/22/2021] [Indexed: 05/13/2023]
Abstract
X-ray photoemission (XPS) and near edge x-ray absorption fine structure (NEXAFS) spectroscopy play an important role in investigating the structure and electronic structure of materials and surfaces.Ab initiosimulations provide crucial support for the interpretation of complex spectra containing overlapping signatures. Approximate core-hole simulation methods based on density functional theory (DFT) such as the delta-self-consistent-field (ΔSCF) method or the transition potential (TP) method are widely used to predictK-shell XPS and NEXAFS signatures of organic molecules, inorganic materials and metal-organic interfaces at reliable accuracy and affordable computational cost. We present the numerical and technical details of our variants of the ΔSCF and TP method (coined ΔIP-TP) to simulate XPS and NEXAFS transitions. Using exemplary molecules in gas-phase, in bulk crystals, and at metal-organic interfaces, we systematically assess how practical simulation choices affect the stability and accuracy of simulations. These include the choice of exchange-correlation functional, basis set, the method of core-hole localization, and the use of periodic boundary conditions (PBC). We particularly focus on the choice of aperiodic or periodic description of systems and how spurious charge effects in periodic calculations affect the simulation outcomes. For the benefit of practitioners in the field, we discuss sensible default choices, limitations of the methods, and future prospects.
Collapse
Affiliation(s)
- Benedikt P Klein
- Department of Chemistry, University of Warwick, Gibbet Hill Rd, Coventry, CV4 7AL, United Kingdom
- Diamond Light Source, Harwell Science and Innovation Campus, Didcot, OX11 0DE, United Kingdom
| | - Samuel J Hall
- Department of Chemistry, University of Warwick, Gibbet Hill Rd, Coventry, CV4 7AL, United Kingdom
- MAS CDT, Senate House, University of Warwick, Gibbet Hill Rd, Coventry, CV4 7AL, United Kingdom
| | - Reinhard J Maurer
- Department of Chemistry, University of Warwick, Gibbet Hill Rd, Coventry, CV4 7AL, United Kingdom
| |
Collapse
|
43
|
Keller L, Blum V, Rinke P, Golze D. Relativistic correction scheme for core-level binding energies from GW. J Chem Phys 2020; 153:114110. [DOI: 10.1063/5.0018231] [Citation(s) in RCA: 11] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Levi Keller
- Department of Applied Physics, Aalto University, Otakaari 1, FI-02150 Espoo, Finland
| | - Volker Blum
- Department of Mechanical Engineering and Materials Science, Duke University, Durham, North Carolina 27708,
USA
| | - Patrick Rinke
- Department of Applied Physics, Aalto University, Otakaari 1, FI-02150 Espoo, Finland
| | - Dorothea Golze
- Department of Applied Physics, Aalto University, Otakaari 1, FI-02150 Espoo, Finland
| |
Collapse
|
44
|
Matthews DA, Cheng L, Harding ME, Lipparini F, Stopkowicz S, Jagau TC, Szalay PG, Gauss J, Stanton JF. Coupled-cluster techniques for computational chemistry: The CFOUR program package. J Chem Phys 2020; 152:214108. [DOI: 10.1063/5.0004837] [Citation(s) in RCA: 214] [Impact Index Per Article: 53.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/30/2022] Open
Affiliation(s)
- Devin A. Matthews
- Department of Chemistry, Southern Methodist University, Dallas, Texas 75275, USA
| | - Lan Cheng
- Department of Chemistry, Johns Hopkins University, Baltimore, Maryland 21218, USA
| | - Michael E. Harding
- Institut für Physikalische Chemie, Karlsruher Institut für Technologie (KIT), Kaiserstr. 12, D-76131 Karlsruhe, Germany
| | - Filippo Lipparini
- Dipartimento di Chimica e Chimica Industriale, Università di Pisa, Via G. Moruzzi 13, I-56124 Pisa, Italy
| | - Stella Stopkowicz
- Department Chemie, Johannes Gutenberg-Universität Mainz, Duesbergweg 10-14, D-55128 Mainz, Germany
| | - Thomas-C. Jagau
- Department of Chemistry, University of Munich (LMU), Butenandtstr. 5-13, D-81377 Munich, Germany
| | - Péter G. Szalay
- ELTE Eötvös Loránd University, Institute of Chemistry, Pázmány Péter sétány 1/A, H-1117 Budapest, Hungary
| | - Jürgen Gauss
- Department Chemie, Johannes Gutenberg-Universität Mainz, Duesbergweg 10-14, D-55128 Mainz, Germany
| | - John F. Stanton
- Quantum Theory Project, Departments of Chemistry and Physics, University of Florida, Gainesville, Florida 32611, USA
| |
Collapse
|
45
|
Maier TM, Ikabata Y, Nakai H. Relativistic local hybrid functionals and their impact on 1s core orbital energies. J Chem Phys 2020; 152:214103. [DOI: 10.1063/5.0010400] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Toni M. Maier
- Department of Chemistry and Biochemistry, School of Advanced Science and Engineering, Waseda University, 3-4-1 Okubo, Shinjuku-ku, Tokyo 169-8555, Japan
- Technische Universität Berlin, Institut für Chemie, Theoretische Chemie/Quantenchemie, Sekr. C7, Straße des 17. Juni 135, D-10623 Berlin, Germany
| | - Yasuhiro Ikabata
- Waseda Research Institute for Science and Engineering, Waseda University, 3-4-1 Okubo, Shinjuku-ku, Tokyo 169-8555, Japan
| | - Hiromi Nakai
- Department of Chemistry and Biochemistry, School of Advanced Science and Engineering, Waseda University, 3-4-1 Okubo, Shinjuku-ku, Tokyo 169-8555, Japan
- Waseda Research Institute for Science and Engineering, Waseda University, 3-4-1 Okubo, Shinjuku-ku, Tokyo 169-8555, Japan
- Elements Strategy Initiative for Catalysts and Batteries (ESICB), Kyoto University, Katsura, Kyoto 615-8520, Japan
| |
Collapse
|
46
|
Zhao L, Neuscamman E. Excited state mean-field theory without automatic differentiation. J Chem Phys 2020; 152:204112. [DOI: 10.1063/5.0003438] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/16/2023] Open
Affiliation(s)
- Luning Zhao
- Department of Chemistry, University of Washington, Seattle, Washington 98195, USA
| | - Eric Neuscamman
- Department of Chemistry, University of California, Berkeley, California 94720, USA
- Chemical Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, USA
| |
Collapse
|
47
|
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]
|
48
|
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
|
49
|
Golze D, Keller L, Rinke P. Accurate Absolute and Relative Core-Level Binding Energies from GW. J Phys Chem Lett 2020; 11:1840-1847. [PMID: 32043890 PMCID: PMC7735733 DOI: 10.1021/acs.jpclett.9b03423] [Citation(s) in RCA: 50] [Impact Index Per Article: 12.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/20/2019] [Accepted: 02/11/2020] [Indexed: 05/13/2023]
Abstract
We present an accurate approach to compute X-ray photoelectron spectra based on the GW Green's function method that overcomes the shortcomings of common density functional theory approaches. GW has become a popular tool to compute valence excitations for a wide range of materials. However, core-level spectroscopy is thus far almost uncharted in GW. We show that single-shot perturbation calculations in the G0W0 approximation, which are routinely used for valence states, cannot be applied for core levels and suffer from an extreme, erroneous transfer of spectral weight to the satellite spectrum. The correct behavior can be restored by partial self-consistent GW schemes or by using hybrid functionals with almost 50% of exact exchange as a starting point for G0W0. We also include relativistic corrections and present a benchmark study for 65 molecular 1s excitations. Our absolute and relative GW core-level binding energies agree within 0.3 and 0.2 eV with experiment, respectively.
Collapse
Affiliation(s)
- Dorothea Golze
- Department of Applied Physics, Aalto University, Otakaari 1, FI-02150 Espoo, Finland
| | - Levi Keller
- Department of Applied Physics, Aalto University, Otakaari 1, FI-02150 Espoo, Finland
| | - Patrick Rinke
- Department of Applied Physics, Aalto University, Otakaari 1, FI-02150 Espoo, Finland
| |
Collapse
|
50
|
Hait D, Head-Gordon M. Excited State Orbital Optimization via Minimizing the Square of the Gradient: General Approach and Application to Singly and Doubly Excited States via Density Functional Theory. J Chem Theory Comput 2020; 16:1699-1710. [PMID: 32017554 DOI: 10.1021/acs.jctc.9b01127] [Citation(s) in RCA: 99] [Impact Index Per Article: 24.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
We present a general approach to converge excited state solutions to any quantum chemistry orbital optimization process, without the risk of variational collapse. The resulting square gradient minimization (SGM) approach only requires analytic energy/Lagrangian orbital gradients and merely costs 3 times as much as ground state orbital optimization (per iteration), when implemented via a finite difference approach. SGM is applied to both single determinant ΔSCF and spin-purified restricted open-shell Kohn-Sham (ROKS) approaches to study the accuracy of orbital optimized DFT excited states. It is found that SGM can converge challenging states where the maximum overlap method (MOM) or analogues either collapse to the ground state or fail to converge. We also report that ΔSCF/ROKS predict highly accurate excitation energies for doubly excited states (which are inaccessible via TDDFT). Singly excited states obtained via ROKS are also found to be quite accurate, especially for Rydberg states that frustrate (semi)local TDDFT. Our results suggest that orbital optimized excited state DFT methods can be used to push past the limitations of TDDFT to doubly excited, charge-transfer, or Rydberg states, making them a useful tool for the practical quantum chemist's toolbox for studying excited states in large systems.
Collapse
Affiliation(s)
- 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
| | - 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
|