1
|
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
|
2
|
Gałyńska M, Boguslawski K. Benchmarking Ionization Potentials from pCCD Tailored Coupled Cluster Models. J Chem Theory Comput 2024; 20:4182-4195. [PMID: 38752491 PMCID: PMC11137826 DOI: 10.1021/acs.jctc.4c00172] [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/09/2024] [Revised: 05/04/2024] [Accepted: 05/07/2024] [Indexed: 05/29/2024]
Abstract
The ionization potential (IP) is an important parameter providing essential insights into the reactivity of chemical systems. IPs are also crucial for designing, optimizing, and understanding the functionality of modern technological devices. We recently showed that limiting the CC ansatz to the seniority-zero sector proves insufficient in predicting reliable and accurate ionization potentials within an IP equation-of-motion coupled-cluster formalism. Specifically, the absence of dynamical correlation in the seniority-zero pair coupled cluster doubles (pCCD) model led to unacceptably significant errors of approximately 1.5 eV. In this work, we aim to explore the impact of dynamical correlation and the choice of the molecular orbital basis (canonical vs localized) in CC-type methods targeting 230 ionized states in 70 molecules, comprising small organic molecules, medium-sized organic acceptors, and nucleobases. We focus on pCCD-based approaches as well as the conventional IP-EOM-CCD and IP-EOM-CCSD. Their performance is compared to the CCSD(T) or CCSDT equivalent and experimental reference data. Our statistical analysis reveals that all investigated frozen-pair coupled cluster methods exhibit similar performance, with differences in errors typically within chemical accuracy (1 kcal/mol or 0.05 eV). Notably, the effect of the molecular orbital basis, such as canonical Hartree-Fock or natural pCCD-optimized orbitals, on the IPs is marginal if dynamical correlation is accounted for. Our study suggests that triple excitations are crucial in achieving chemical accuracy in IPs when modeling electron detachment processes with pCCD-based methods.
Collapse
Affiliation(s)
- Marta Gałyńska
- Institute of Physics, Faculty of Physics,
Astronomy, and Informatics, Nicolaus Copernicus
University in Toruń, Grudziadzka 5, 87-100 Toruń, Poland
| | - Katharina Boguslawski
- Institute of Physics, Faculty of Physics,
Astronomy, and Informatics, Nicolaus Copernicus
University in Toruń, Grudziadzka 5, 87-100 Toruń, Poland
| |
Collapse
|
3
|
Mukhopadhyay T, Jangid B, Dutta AK. State-specific frozen natural orbital for reduced-cost algebraic diagrammatic construction calculations: The application to ionization problem. J Chem Phys 2023; 159:084113. [PMID: 37638624 DOI: 10.1063/5.0160024] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/29/2023] [Accepted: 08/03/2023] [Indexed: 08/29/2023] Open
Abstract
We have developed a reduced-cost algebraic diagrammatic construction (ADC) method based on state-specific frozen natural orbital and natural auxiliary functions. The newly developed method has been benchmarked on the GW100 test set for the ionization problem. The use of state-specific natural orbitals drastically reduces the size of the virtual space with a systematically controllable accuracy and offers a significant speedup over the standard ionization potential (IP)-ADC(3) method. The accuracy of the method can be controlled by two thresholds and nearly a black box to use. The inclusion of the perturbative correction significantly improves the accuracy of the calculated IP values, and the efficiency of the method has been demonstrated by calculating the IP of a molecule with 60 atoms and more than 2216 basis functions.
Collapse
Affiliation(s)
- Tamoghna Mukhopadhyay
- 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
| | - Achintya Kumar Dutta
- Department of Chemistry, Indian Institute of Technology Bombay, Powai, Mumbai 400076, India
| |
Collapse
|
4
|
Dreuw A, Papapostolou A, Dempwolff AL. Algebraic Diagrammatic Construction Schemes Employing the Intermediate State Formalism: Theory, Capabilities, and Interpretation. J Phys Chem A 2023; 127:6635-6646. [PMID: 37498297 DOI: 10.1021/acs.jpca.3c02761] [Citation(s) in RCA: 11] [Impact Index Per Article: 11.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 07/28/2023]
Abstract
Algebraic diagrammatic construction (ADC) schemes represent a family of ab initio methods for the calculation of excited electronic states and electron-detached and -attached states. All ADC methods have been demonstrated to possess great potential for molecular applications, e.g., for the calculation of absorption or photoelectron spectra or electron attachment processes. ADC originates from Green's function or propagator theory; however, most recent ADC developments heavily rely on the intermediate state representation or effective Liouvillian formalisms, which comprise new ADC methods and computational schemes for high-order properties. The different approaches for the calculation of excitation energies, ionization potentials, and electron affinities are intimately related, and they provide a coherent description of these quantities at equivalent levels of theory and with comparable errors. Most quantum chemical program packages contain ADC methods; however, the most complete ADC suite of methods can be found in the recent release of Q-Chem.
Collapse
Affiliation(s)
- Andreas Dreuw
- Interdisciplinary Center for Scientific Computing, Ruprecht-Karls University, 69120 Heidelberg, Germany
| | - Antonia Papapostolou
- Interdisciplinary Center for Scientific Computing, Ruprecht-Karls University, 69120 Heidelberg, Germany
| | - Adrian L Dempwolff
- Interdisciplinary Center for Scientific Computing, Ruprecht-Karls University, 69120 Heidelberg, Germany
| |
Collapse
|
5
|
Mamache S, Gałyńska M, Boguslawski K. Benchmarking ionization potentials using the simple pCCD model. Phys Chem Chem Phys 2023. [PMID: 37378457 DOI: 10.1039/d3cp01963b] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/29/2023]
Abstract
The electron-detachment energy is measured by the ionization potential (IP). As a result, it is a fundamental, observable and important molecular electronic signature in photoelectron spectroscopy. A precise theoretical prediction of electron-detachment energies or ionization potentials is essential for organic optoelectronic systems like transistors, solar cells, or light-emitting diodes. In this work, we benchmark the performance of the recently presented IP variant of the equation-of-motion pair coupled cluster doubles (IP-EOM-pCCD) model to determine IPs. Specifically, the predicted ionization energies are compared to experimental results and higher-order coupled cluster theories based on statistically assessing 201 electron-detached states of 41 organic molecules for three different molecular orbital basis sets and two sets of particle-hole operators. While IP-EOM-pCCD features a reasonable spread and skewness of ionization energies, its mean error and standard deviation differ by up to 1.5 eV from reference data. Our study, thus, highlights the importance of dynamical correlation to reliably predict IPs from a pCCD reference function in small organic molecules.
Collapse
Affiliation(s)
- Saddem Mamache
- Institute of Physics, Faculty of Physics, Astronomy and Informatics, Nicolaus Copernicus University in Toruń, Grudziadzka 5, 87-100 Toruń, Poland.
| | - Marta Gałyńska
- Institute of Physics, Faculty of Physics, Astronomy and Informatics, Nicolaus Copernicus University in Toruń, Grudziadzka 5, 87-100 Toruń, Poland.
| | - Katharina Boguslawski
- Institute of Physics, Faculty of Physics, Astronomy and Informatics, Nicolaus Copernicus University in Toruń, Grudziadzka 5, 87-100 Toruń, Poland.
| |
Collapse
|
6
|
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
|
7
|
Ceriotti M, Jensen L, Manolopoulos DE, Martinez T, Reichman DR, Sciortino F, Sherrill CD, Shi Q, Vega C, Wang LS, Weiss EA, Zhu X, Stein J, Lian T. 2021 JCP Emerging Investigator Special Collection. J Chem Phys 2023; 158:060401. [PMID: 36792492 DOI: 10.1063/5.0143234] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/15/2023] Open
|
8
|
Banerjee S, Sokolov AY. Non-Dyson Algebraic Diagrammatic Construction Theory for Charged Excitations in Solids. J Chem Theory Comput 2022; 18:5337-5348. [PMID: 35976918 DOI: 10.1021/acs.jctc.2c00565] [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/29/2022]
Abstract
We present the first implementation and applications of non-Dyson algebraic diagrammatic construction theory for charged excitations in three-dimensional periodic solids (EA/IP-ADC). The EA/IP-ADC approach has a computational cost similar to the ground-state Møller-Plesset perturbation theory, enabling efficient calculations of a variety of crystalline excited-state properties (e.g., band structure, band gap, density of states) sampled in the Brillouin zone. We use EA/IP-ADC to compute the quasiparticle band structures and band gaps of several materials (from large-gap atomic and ionic solids to small-gap semiconductors) and analyze the errors of EA/IP-ADC approximations up to the third order in perturbation theory. Our work also reports the first-ever calculations of ground-state properties (equation-of-state and lattice constants) of three-dimensional crystalline systems using a periodic implementation of third-order Møller-Plesset perturbation theory (MP3).
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
|
9
|
Opoku E, Pawłowski F, Ortiz JV. Electron Propagator Self-Energies versus Improved GW100 Vertical Ionization Energies. J Chem Theory Comput 2022; 18:4927-4944. [PMID: 35822816 DOI: 10.1021/acs.jctc.2c00502] [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/28/2022]
Abstract
Ab initio electron propagator (EP) methods that are free of adjustable parameters in their self-energy formulae and in the generation of their orbital bases have been applied to the calculation of the lowest vertical ionization energies (VIEs) of the GW100 set. An improved set of standard results accompanied by irreducible representation assignments has been produced indirectly with coupled-cluster singles and doubles plus perturbative triples, i.e., CCSD(T), total energy differences at initial-state geometries reoptimized (in 28 cases) with the largest applicable point groups. The best compromises of accuracy and efficiency belong to a new generation of EP self-energies, several members of which may be derived from an intermediately normalized, Hermitized super-operator metric. The following diagonal self-energy methods are optimal: opposite-spin non-Dyson second order (os-nD-D2), approximately renormalized partial third order (P3+), approximately renormalized quasiparticle third order (Q3+), and non-Dyson approximately renormalized linear third order version B (nD-L3+B). Their mean absolute errors (MAEs) in electron volts and arithmetic scaling factors expressed in terms of occupied (O) and virtual (V) orbital dimensions are, respectively, (0.18, OV2), (0.14, O2V3), (0.15, O2V3), and (0.11, OV4). The 0.06 eV MAE for the non-diagonal, sixth-power (O2V4) Brueckner doubles, triple-field operator (BD-T1) EP method is exceeded by the 0.1 eV MAE with respect to experiments in seventh-power, ΔCCSD(T) calculations and indicates that BD-T1 may serve as a direct, spin-symmetry-conserving alternative in the generation of standard results for VIEs of larger, closed-shell molecules.
Collapse
Affiliation(s)
- Ernest Opoku
- Department of Chemistry and Biochemistry, Auburn University, Auburn, Alabama 36849-5312, United States
| | - Filip Pawłowski
- Department of Chemistry and Biochemistry, Auburn University, Auburn, Alabama 36849-5312, United States
| | - J V Ortiz
- Department of Chemistry and Biochemistry, Auburn University, Auburn, Alabama 36849-5312, United States
| |
Collapse
|
10
|
Liu J, Matthews DA, Cheng L. Quadratic Unitary Coupled-Cluster Singles and Doubles Scheme: Efficient Implementation, Benchmark Study, and Formulation of an Extended Version. J Chem Theory Comput 2022; 18:2281-2291. [PMID: 35312299 DOI: 10.1021/acs.jctc.1c01210] [Citation(s) in RCA: 4] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
An efficient implementation of the quadratic unitary coupled-cluster singles and doubles (qUCCSD) scheme for calculations of electronic ground and excited states using an unrestricted molecular spin-orbital formulation and an efficient tensor contraction library is reported. The accuracy of the qUCCSD scheme and the efficiency of the present implementation are demonstrated using extensive benchmark calculations of excitation energies and an application to S0 → S1 vertical excitation energies for cis- and trans-4a,4b-dihydrotriphenylene. The qUCCSD scheme has been shown to provide improved excitation energies compared with the UCC3 scheme formulated based on perturbation theory. A UCC truncation scheme that can provide excitation energies correct through the fourth order is also presented to further improve the accuracy of the qUCCSD scheme.
Collapse
Affiliation(s)
- Junzi Liu
- School of Chemistry and Biological Engineering, University of Science and Technology Beijing, Beijing, 100083, China
| | - Devin A Matthews
- Department of Chemistry, Southern Methodist University, Dallas, Texas 75275, United States
| | - Lan Cheng
- Department of Chemistry, The Johns Hopkins University, Baltimore, Maryland 21218, United States
| |
Collapse
|
11
|
Hodecker M, Dempwolff AL, Schirmer J, Dreuw A. Theoretical analysis and comparison of unitary coupled-cluster and algebraic-diagrammatic construction methods for ionization. J Chem Phys 2022; 156:074104. [DOI: 10.1063/5.0070967] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/01/2023] Open
Affiliation(s)
- Manuel Hodecker
- Interdisciplinary Center for Scientific Computing, Heidelberg University, Im Neuenheimer Feld 205, 69120 Heidelberg, Germany
| | - Adrian L. Dempwolff
- Interdisciplinary Center for Scientific Computing, Heidelberg University, Im Neuenheimer Feld 205, 69120 Heidelberg, Germany
| | - Jochen Schirmer
- Theoretical Chemistry, Heidelberg University, Im Neuenheimer Feld 229, 69120 Heidelberg, Germany
| | - Andreas Dreuw
- Interdisciplinary Center for Scientific Computing, Heidelberg University, Im Neuenheimer Feld 205, 69120 Heidelberg, Germany
| |
Collapse
|
12
|
Dempwolff AL, Hodecker M, Dreuw A. Vertical ionization potential benchmark for unitary coupled-cluster and algebraic-diagrammatic construction methods. J Chem Phys 2022; 156:054114. [DOI: 10.1063/5.0079047] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Adrian L. Dempwolff
- Interdisciplinary Center for Scientific Computing, Heidelberg University, Im Neuenheimer Feld 205, 69120 Heidelberg, Germany
| | - Manuel Hodecker
- Interdisciplinary Center for Scientific Computing, Heidelberg University, Im Neuenheimer Feld 205, 69120 Heidelberg, Germany
- Department of Theoretical Chemistry and Biology, KTH Royal Institute of Technology, Malvinas Väg 10, 114 28 Stockholm, Sweden
| | - Andreas Dreuw
- Interdisciplinary Center for Scientific Computing, Heidelberg University, Im Neuenheimer Feld 205, 69120 Heidelberg, Germany
| |
Collapse
|
13
|
Totani R, Ljubić I, Ciavardini A, Grazioli C, Galdenzi F, de Simone M, Coreno M. Frontier orbital stability of nitroxyl organic radicals probed by means of inner shell resonantly enhanced valence band photoelectron spectroscopy. Phys Chem Chem Phys 2022; 24:1993-2003. [PMID: 35018901 DOI: 10.1039/d1cp05264k] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
We have investigated the frontier orbitals of persistent organic radicals known as nitroxyls by resonant photoelectron spectroscopy (ResPES) under inner shell excitation. By means of this site-specific technique, we were able to disentangle the different atomic contributions to the outer valence molecular orbitals and examine several core-hole relaxation pathways involving the singly occupied molecular orbital (SOMO) localized on the nitroxyl group. To interpret the ResPES intensity trends, especially the strong enhancement of the SOMO ionized state at the N K-edge, we computed the Dyson spin orbitals (DSOs) pertaining to the transitions between the core-excited initial states and several of the singly ionized valence final states. We found that the computed vertical valence ionization potentials and norms of the DSOs are reasonably reliable when based on the long-range corrected CAM-B3LYP density functional. Thanks to their unpaired electrons, nitroxyls have recently found application in technological fields implying a spin control, such as spintronics and quantum computing. The present findings on the electronic structure of nitroxyl persistent radicals furnish important hints for their implementation in technological devices and, more in general, for the synthesis of new and stable organic radicals with tailored properties.
Collapse
Affiliation(s)
- R Totani
- ISM-CNR, LD2 Unit, Basovizza Area Science Park, 34149 Trieste, Italy.
| | - I Ljubić
- Division of Physical Chemistry, Ruđer Bošković Institute, Bijenička cesta 54, HR-10000 Zagreb, Croatia.
| | - A Ciavardini
- University of Nova Gorica, Vipavska 13, SI-5000 Nova Gorica, Slovenia
| | - C Grazioli
- IOM-CNR, Laboratorio TASC, Basovizza SS-14, km 163.5, 34149 Trieste, Italy
| | - F Galdenzi
- University of Nova Gorica, Vipavska 13, SI-5000 Nova Gorica, Slovenia
| | - M de Simone
- IOM-CNR, Laboratorio TASC, Basovizza SS-14, km 163.5, 34149 Trieste, Italy
| | - M Coreno
- ISM-CNR, LD2 Unit, Basovizza Area Science Park, 34149 Trieste, Italy.
| |
Collapse
|
14
|
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
|
15
|
Recent progress in electron-propagator, extended-Koopmans-theorem and self-consistent-field approaches to the interpretation and prediction of electron binding energies. ADVANCES IN QUANTUM CHEMISTRY 2022. [DOI: 10.1016/bs.aiq.2022.03.004] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/18/2022]
|
16
|
Opoku E, Pawłowski F, Ortiz JV. A new generation of diagonal self-energies for the calculation of electron removal energies. J Chem Phys 2021; 155:204107. [PMID: 34852487 DOI: 10.1063/5.0070849] [Citation(s) in RCA: 16] [Impact Index Per Article: 5.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Abstract
A new generation of diagonal self-energy approximations in ab initio electron propagator theory for the calculation of electron removal energies of molecules and molecular ions has been derived from an intermediately normalized, Hermitized super-operator metric. These methods and widely used antecedents such as the outer valence Green's function and the approximately renormalized partial third order method are tested with respect to a dataset of vertical ionization energies generated with a valence, triple-ζ, correlation-consistent basis set and a converged series of many-body calculations whose accuracy approaches that of full configuration interaction. Several modifications of the diagonal second-order self-energy, a version of G0W0 theory based on Tamm-Dancoff excitations and several non-diagonal self-energies are also included in the tests. All new methods employ canonical Hartree-Fock orbitals. No adjustable or empirical parameters appear. A hierarchy of methods with optimal accuracy for a given level of computational efficiency is established. Several widely used diagonal self-energy methods are rendered obsolete by the new hierarchy whose members, in order of increasing accuracy, are (1) the opposite-spin non-Dyson diagonal second-order or os-nD-D2, (2) the approximately renormalized third-order quasiparticle or Q3+, (3) the renormalized third-order quasiparticle or RQ3, (4) the approximately renormalized linear third-order or L3+, and (5) the renormalized linear third-order or RL3 self-energies.
Collapse
Affiliation(s)
- Ernest Opoku
- Department of Chemistry and Biochemistry, Auburn University, Auburn, Alabama 36849-5312, USA
| | - Filip Pawłowski
- Department of Chemistry and Biochemistry, Auburn University, Auburn, Alabama 36849-5312, USA
| | - J V Ortiz
- Department of Chemistry and Biochemistry, Auburn University, Auburn, Alabama 36849-5312, USA
| |
Collapse
|
17
|
Liu J, Cheng L. Unitary coupled-cluster based self-consistent polarization propagator theory: A quadratic unitary coupled-cluster singles and doubles scheme. J Chem Phys 2021; 155:174102. [PMID: 34742195 DOI: 10.1063/5.0062090] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
The development of a quadratic unitary coupled-cluster singles and doubles (qUCCSD) based self-consistent polarization propagator method is reported. We present a simple strategy for truncating the commutator expansion of the unitary version of coupled-cluster transformed Hamiltonian H̄. The qUCCSD method for the electronic ground state includes up to double commutators for the amplitude equations and up to cubic commutators for the energy expression. The qUCCSD excited-state eigenvalue equations include up to double commutators for the singles-singles block of H̄, single commutators for the singles-doubles and doubles-singles blocks, and the bare Hamiltonian for the doubles-doubles block. Benchmark qUCCSD calculations of the ground-state properties and excitation energies for representative molecules demonstrate significant improvement of the accuracy and robustness over the previous UCC3 scheme derived using Møller-Plesset perturbation theory.
Collapse
Affiliation(s)
- Junzi Liu
- Department of Chemistry, The Johns Hopkins University, Baltimore, Maryland 21218, USA
| | - Lan Cheng
- Department of Chemistry, The Johns Hopkins University, Baltimore, Maryland 21218, USA
| |
Collapse
|
18
|
Mazin IM, Sokolov AY. Multireference Algebraic Diagrammatic Construction Theory for Excited States: Extended Second-Order Implementation and Benchmark. J Chem Theory Comput 2021; 17:6152-6165. [PMID: 34553937 DOI: 10.1021/acs.jctc.1c00684] [Citation(s) in RCA: 7] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
We present an implementation and benchmark of new approximations in multireference algebraic diagrammatic construction theory for simulations of neutral electronic excitations and UV/vis spectra of strongly correlated molecular systems (MR-ADC). Following our work on the first-order MR-ADC approximation [J. Chem. Phys. 2018, 149, 204113], we report the strict and extended second-order MR-ADC methods (MR-ADC(2) and MR-ADC(2)-X) that combine the description of static and dynamic electron correlation in the ground and excited electronic states without relying on state-averaged reference wave functions. We present an extensive benchmark of the new MR-ADC methods for excited states in several small molecules, including the carbon dimer, ethylene, and butadiene. Our results demonstrate that, for weakly correlated electronic states, the MR-ADC(2) and MR-ADC(2)-X methods outperform the third-order single-reference ADC approximation and are competitive with the results from equation-of-motion coupled cluster theory. For states with multireference character, the performance of the MR-ADC methods is similar to that of an N-electron valence perturbation theory. In contrast to conventional multireference perturbation theories, the MR-ADC methods have many attractive features, such as a straightforward and efficient calculation of excited-state properties and a direct access to excitations outside of the frontier (active) orbitals.
Collapse
Affiliation(s)
- Ilia M Mazin
- 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
|
19
|
Backhouse OJ, Santana-Bonilla A, Booth GH. Scalable and Predictive Spectra of Correlated Molecules with Moment Truncated Iterated Perturbation Theory. J Phys Chem Lett 2021; 12:7650-7658. [PMID: 34351782 DOI: 10.1021/acs.jpclett.1c02383] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/13/2023]
Abstract
A reliable and efficient computation of the entire single-particle spectrum of correlated molecules is an outstanding challenge in the field of quantum chemistry, with standard density functional theory approaches often giving an inadequate description of excitation energies and gaps. In this work, we expand upon a recently introduced approach that relies on a fully self-consistent many-body perturbation theory coupled to a nonperturbative truncation of the effective dynamics at each step. We show that this yields a low-scaling and accurate method across a diverse benchmark test set that is capable of treating moderate levels of strong correlation effects, and we detail an efficient implementation for applications involving up to ∼1000 orbitals on parallel resources. We then use this method to characterize the spectral properties of the antimalarial drug molecule artemisinin, resolving discrepancies in previous works concerning the active sites of the lowest-energy fundamental excitations of the system.
Collapse
Affiliation(s)
- Oliver J Backhouse
- Department of Physics, King's College London, Strand, London WC2R 2LS, U.K
| | | | - George H Booth
- Department of Physics, King's College London, Strand, London WC2R 2LS, U.K
| |
Collapse
|