1
|
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
|
2
|
Abraham V, Harsha G, Zgid D. Relativistic Fully Self-Consistent GW for Molecules: Total Energies and Ionization Potentials. J Chem Theory Comput 2024; 20:4579-4590. [PMID: 38778459 DOI: 10.1021/acs.jctc.4c00075] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/25/2024]
Abstract
The fully self-consistent GW (scGW) method with an iterative solution of the Dyson equation provides a consistent approach for describing the ground and excited states without any dependence on the mean-field reference. In this work, we present a relativistic version of scGW for molecules containing heavy elements using the exact two-component (X2C) Coulomb approximation. We benchmark SOC-81 data set containing closed shell heavy elements for the first ionization potential using the fully self-consistent GW as well as one-shot GW. The self-consistent GW provides superior results compared to G0W0 with PBE reference and comparable results to G0W0 with PBE0 while also removing the starting point dependence. The photoelectron spectra obtained at the X2C level demonstrate very good agreement with the experimental spectra. We also observe that scGW provides very good estimation of ionization potential for the inner d-shell orbitals. Additionally, using the well-conserved total energy, we investigate the equilibrium bond length and harmonic frequencies of a few halogen dimers using scGW. Overall, our findings demonstrate the applicability of the fully self-consistent GW method for accurate ionization potential, photoelectron spectra, and total energies in finite systems with heavy elements with a reasonable computational scaling.
Collapse
Affiliation(s)
- Vibin Abraham
- Department of Chemistry, University of Michigan, Ann Arbor, Michigan 48109, United States
| | - Gaurav Harsha
- Department of Chemistry, University of Michigan, Ann Arbor, Michigan 48109, United States
| | - Dominika Zgid
- Department of Chemistry, University of Michigan, Ann Arbor, Michigan 48109, United States
- Department of Physics and Astronomy, University of Michigan, Ann Arbor, Michigan 48109, United States
| |
Collapse
|
3
|
Bruneval F, Förster A. Fully Dynamic G3 W2 Self-Energy for Finite Systems: Formulas and Benchmark. J Chem Theory Comput 2024; 20:3218-3230. [PMID: 38603811 DOI: 10.1021/acs.jctc.4c00090] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 04/13/2024]
Abstract
Over the years, Hedin's GW self-energy has been proven to be a rather accurate and simple approximation to evaluate electronic quasiparticle energies in solids and in molecules. Attempts to improve over the simple GW approximation, the so-called vertex corrections, have been constantly proposed in the literature. Here, we derive, analyze, and benchmark the complete second-order term in the screened Coulomb interaction W for finite systems. This self-energy named G3W2 contains all the possible time orderings that combine 3 Green's functions G and 2 dynamic W. We present the analytic formula and its imaginary frequency counterpart, with the latter allowing us to treat larger molecules. The accuracy of the G3W2 self-energy is evaluated on well-established benchmarks (GW100, Acceptor 24, and Core 65) for valence and core quasiparticle energies. Its link with the simpler static approximation, named SOSEX for static screened second-order exchange, is analyzed, which leads us to propose a more consistent approximation named 2SOSEX. In the end, we find that neither the G3W2 self-energy nor any of the investigated approximations to it improve over one-shot G0W0 with a good starting point. Only quasi-particle self-consistent GW HOMO energies are slightly improved by addition of the G3W2 self-energy correction. We show that this is due to the self-consistent update of the screened Coulomb interaction, leading to an overall sign change of the vertex correction to the frontier quasiparticle energies.
Collapse
Affiliation(s)
- Fabien Bruneval
- Université Paris-Saclay, CEA, Service de recherche en Corrosion et Comportement des Matériaux, SRMP, 91191 Gif-sur-Yvette, France
| | - Arno Förster
- Theoretical Chemistry, Vrije Universiteit, De Boelelaan 1108, 1081 HZ Amsterdam, The Netherlands
| |
Collapse
|
4
|
Amblard D, Blase X, Duchemin I. Static versus dynamically polarizable environments within the many-body GW formalism. J Chem Phys 2024; 160:154104. [PMID: 38624115 DOI: 10.1063/5.0203637] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/15/2024] [Accepted: 03/28/2024] [Indexed: 04/17/2024] Open
Abstract
Continuum- or discrete-polarizable models for the study of optoelectronic processes in embedded subsystems rely mostly on the restriction of the surrounding electronic dielectric response to its low frequency limit. Such a description hinges on the assumption that the electrons in the surrounding medium react instantaneously to any excitation in the central subsystem, thus treating the environment in the adiabatic limit. Exploiting a recently developed embedded GW formalism with an environment described at the fully ab initio level, we assess the merits of the adiabatic limit with respect to an environment where the full dynamics of the dielectric response are considered. Furthermore, we show how to properly take the static limit of the environment's susceptibility by introducing the so-called Coulomb-hole and screened-exchange contributions to the reaction field. As a first application, we consider a C60 molecule at the surface of a C60 crystal, namely, a case where the dynamics of the embedded and embedding subsystems are similar. The common adiabatic assumption, when properly treated, generates errors below 10% on the polarization energy associated with frontier energy levels and associated energy gaps. Finally, we consider a water molecule inside a metallic nanotube, the worst case for the environment's adiabatic limit. The error on the gap polarization energy remains below 10%, even though the error on the frontier orbital polarization energies can reach a few tenths of an electronvolt.
Collapse
Affiliation(s)
- David Amblard
- University Grenoble Alpes, CNRS, Inst NEEL, F-38042 Grenoble, France
| | - Xavier Blase
- University Grenoble Alpes, CNRS, Inst NEEL, F-38042 Grenoble, France
| | - Ivan Duchemin
- University Grenoble Alpes, CEA, IRIG-MEM-L_Sim, 38054 Grenoble, France
| |
Collapse
|
5
|
Mejia-Rodriguez D, Kunitsa AA, Fulton JL, Aprà E, Govind N. G0W0 Ionization Potentials of First-Row Transition Metal Aqua Ions. J Phys Chem A 2023; 127:9684-9694. [PMID: 37938891 DOI: 10.1021/acs.jpca.3c04419] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/10/2023]
Abstract
We report computations of the vertical ionization potentials within the GW approximation of the near-complete series of first-row transition metal (V-Cu) aqua ions in their most common oxidation states, i.e., V3+, Cr3+, Cr2+, Mn2+, Fe3+, Fe2+, Co2+, Ni2+, and Cu2+. The d-orbital occupancy of these systems spans a broad range from d2 to d9. All of the structures were first optimized at the density functional theory level using a large cluster of explicit water molecules that are embedded in a continuum solvation model. Vertical ionization potentials were computed with the one-shot G0W0 approach on a range of transition metal ion clusters (6, 18, 40, and 60 explicit water molecules), wherein the convergence with respect to the basis set size was evaluated using the systems with 40 water molecules. We assess the results using three different density functional approximations as starting points for the vertical ionization potential calculations, namely, G0W0@PBE, G0W0@PBE0, and G0W0@r2SCAN. While the predicted ground-state structures are similar to all three exchange-correlation functionals, the vertical ionization potentials were in closer agreement with experiment when using the G0W0@PBE0 and G0W0@r2SCAN approaches, with the r2SCAN-based calculations being significantly less expensive. Computed bond distances and vertical ionization potentials for all structures are in good agreement with available experimental data.
Collapse
Affiliation(s)
- Daniel Mejia-Rodriguez
- Environmental Molecular Sciences Laboratory, Pacific Northwest National Laboratory, Richland, Washington 99352, United States
- Physical and Computational Sciences Directorate, Pacific Northwest National Laboratory, Richland, Washington 99352, United States
| | - Alexander A Kunitsa
- Zapata Computing, Inc., 100 Federal Street, Boston, Massachusetts 02110, United States
| | - John L Fulton
- Physical and Computational Sciences Directorate, Pacific Northwest National Laboratory, Richland, Washington 99352, 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
- Department of Chemistry, University of Washington, Seattle, Washington 98195, United States
| |
Collapse
|
6
|
Mejia-Rodriguez D, Aprà E, Autschbach J, Bauman NP, Bylaska EJ, Govind N, Hammond JR, Kowalski K, Kunitsa A, Panyala A, Peng B, Rehr JJ, Song H, Tretiak S, Valiev M, Vila FD. NWChem: Recent and Ongoing Developments. J Chem Theory Comput 2023; 19:7077-7096. [PMID: 37458314 DOI: 10.1021/acs.jctc.3c00421] [Citation(s) in RCA: 4] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/25/2023]
Abstract
This paper summarizes developments in the NWChem computational chemistry suite since the last major release (NWChem 7.0.0). Specifically, we focus on functionality, along with input blocks, that is accessible in the current stable release (NWChem 7.2.0) and in the "master" development branch, interfaces to quantum computing simulators, interfaces to external libraries, the NWChem github repository, and containerization of NWChem executable images. Some ongoing developments that will be available in the near future are also discussed.
Collapse
Affiliation(s)
- Daniel Mejia-Rodriguez
- Physical Sciences Division, Pacific Northwest National Laboratory, Richland, Washington 99352, United States
| | - Edoardo Aprà
- Environmental Molecular Sciences Laboratory, Pacific Northwest National Laboratory, Richland, Washington 99352, United States
| | - Jochen Autschbach
- Department of Chemistry, University at Buffalo, State University of New York, Buffalo, New York 14260-3000, United States
| | - Nicholas P Bauman
- Physical Sciences Division, Pacific Northwest National Laboratory, Richland, Washington 99352, United States
| | - Eric J Bylaska
- Physical Sciences Division, Pacific Northwest National Laboratory, Richland, Washington 99352, United States
| | - Niranjan Govind
- Physical Sciences Division, Pacific Northwest National Laboratory, Richland, Washington 99352, United States
| | - Jeff R Hammond
- Accelerated Computing, NVIDIA Helsinki Oy, Porkkalankatu 1, 00180 Helsinki, Finland
| | - Karol Kowalski
- Physical Sciences Division, Pacific Northwest National Laboratory, Richland, Washington 99352, United States
| | - Alexander Kunitsa
- Zapata Computing, Inc., 100 Federal Street, Boston, Massachusetts 02110, United States
| | - Ajay Panyala
- Advanced Computing, Mathematics, and Data Division, Pacific Northwest National Laboratory, Richland, Washington 99352, United States
| | - Bo Peng
- Physical Sciences Division, Pacific Northwest National Laboratory, Richland, Washington 99352, United States
| | - John J Rehr
- Department of Physics, University of Washington, Seattle, Washington 98195, United States
| | - Huajing Song
- Physics and Chemistry of Materials, Theoretical Division, Los Alamos National Laboratory, Los Alamos, New Mexico 87545, United States
| | - Sergei Tretiak
- Physics and Chemistry of Materials, Theoretical Division, Los Alamos National Laboratory, Los Alamos, New Mexico 87545, United States
- Center for Integrated Nanotechnologies, Los Alamos National Laboratory, Los Alamos, New Mexico 87545, United States
| | - Marat Valiev
- Physical Sciences Division, Pacific Northwest National Laboratory, Richland, Washington 99352, United States
| | - Fernando D Vila
- Department of Physics, University of Washington, Seattle, Washington 98195, United States
| |
Collapse
|
7
|
Panadés-Barrueta RL, Golze D. Accelerating Core-Level GW Calculations by Combining the Contour Deformation Approach with the Analytic Continuation of W. J Chem Theory Comput 2023; 19:5450-5464. [PMID: 37566917 PMCID: PMC10448726 DOI: 10.1021/acs.jctc.3c00555] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/25/2023] [Indexed: 08/13/2023]
Abstract
In recent years, the GW method has emerged as a reliable tool for computing core-level binding energies. The contour deformation (CD) technique has been established as an efficient, scalable, and numerically stable approach to compute the GW self-energy for deep core excitations. However, core-level GW calculations with CD face the challenge of higher scaling with respect to system size N compared to the conventional quartic scaling in valence-state algorithms. In this work, we present the CD-WAC method [CD with W analytic continuation (AC)], which reduces the scaling of CD applied to the inner shells from O(N5) to O(N4) by employing an AC of the screened Coulomb interaction W. Our proposed method retains the numerical accuracy of CD for the computationally challenging deep core case, yielding mean absolute errors <5 meV for well-established benchmark sets, such as CORE65, for single-shot GW calculations. More extensive testing for different GW flavors proves the reliability of the method. We have confirmed the theoretical scaling by performing scaling experiments on large acene chains and amorphous carbon clusters, achieving speedups of up to 10× for structures of only 116 atoms. This improvement in computational efficiency paves the way for more accurate and efficient core-level GW calculations on larger and more complex systems.
Collapse
Affiliation(s)
| | - Dorothea Golze
- Faculty of Chemistry and Food Chemistry, Technische Universität Dresden, Dresden 01062, Germany
| |
Collapse
|
8
|
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
|
9
|
Fupeng Z, Junxi L, Bomiao Q, Mengmeng L, Shaofeng P, Yanbin W, Qiong S. A comparable DFT study on reaction of CHCl •- with O 3 and S 2O. J Mol Model 2023; 29:85. [PMID: 36867312 DOI: 10.1007/s00894-023-05483-x] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/27/2022] [Accepted: 02/21/2023] [Indexed: 03/04/2023]
Abstract
CONTEXT In this discussion, we began building two model, S2O + CHCl•- and O3 + CHCl•-, using DFT-BHandHLYP method, to study their reactions mechanisms on singlet PES. For this purpose, we hope to explore the effects of the difference between sulfur and oxygen atoms on the CHCl•- anion. Experimentalists and computer scientists may utilize the collected data to generate a wide range of hypotheses for experimental phenomena and predictions, allowing them to realize their full potential. METHODS The ion-molecule reaction mechanism of CHCl•- with S2O and O3 was studied using the DFT-BHandHLYP level of theory with the aug-cc-pVDZ basis set. Our theoretical findings show that Path 6 is the favored reaction pathway for CHCl•- + O3 reaction as identified by the O-abstraction reaction pattern. Comparing to the direct H- and Cl-abstraction mechanisms, the reaction (CHCl•- + S2O) prefers the intramolecular SN2 reaction pattern. Moreover, the calculated results demonstrated that the CHCl•- + S2O reaction is thermodynamically more favorable than the CHCl•- + O3 reaction, which is kinetically more advantageous. As a result, if the required reaction condition in the atmospheric process is met, the O3 reaction will happen more effectively. In terms of kinetics and thermodynamics viewpoints, the CHCl•- anion was very effective in eliminating S2O and O3.
Collapse
Affiliation(s)
- Zhang Fupeng
- Key Laboratory of Environment-Friendly Composite Materials of the State Ethnic Affairs Commission, Key Laboratory for Utility of Environment-Friendly Composite Materials and Biomass in University of Gansu Province, College of Chemical Engineering, Gansu Provincial Biomass Function Composites Engineering Research Center, Northwest Minzu University, Gansu, 730030, Lanzhou, China
| | - Liang Junxi
- Key Laboratory of Environment-Friendly Composite Materials of the State Ethnic Affairs Commission, Key Laboratory for Utility of Environment-Friendly Composite Materials and Biomass in University of Gansu Province, College of Chemical Engineering, Gansu Provincial Biomass Function Composites Engineering Research Center, Northwest Minzu University, Gansu, 730030, Lanzhou, China.
| | - Qi Bomiao
- Key Laboratory of Environment-Friendly Composite Materials of the State Ethnic Affairs Commission, Key Laboratory for Utility of Environment-Friendly Composite Materials and Biomass in University of Gansu Province, College of Chemical Engineering, Gansu Provincial Biomass Function Composites Engineering Research Center, Northwest Minzu University, Gansu, 730030, Lanzhou, China
| | - Lu Mengmeng
- Key Laboratory of Environment-Friendly Composite Materials of the State Ethnic Affairs Commission, Key Laboratory for Utility of Environment-Friendly Composite Materials and Biomass in University of Gansu Province, College of Chemical Engineering, Gansu Provincial Biomass Function Composites Engineering Research Center, Northwest Minzu University, Gansu, 730030, Lanzhou, China
| | - Pang Shaofeng
- Key Laboratory of Environment-Friendly Composite Materials of the State Ethnic Affairs Commission, Key Laboratory for Utility of Environment-Friendly Composite Materials and Biomass in University of Gansu Province, College of Chemical Engineering, Gansu Provincial Biomass Function Composites Engineering Research Center, Northwest Minzu University, Gansu, 730030, Lanzhou, China
| | - Wang Yanbin
- Key Laboratory of Environment-Friendly Composite Materials of the State Ethnic Affairs Commission, Key Laboratory for Utility of Environment-Friendly Composite Materials and Biomass in University of Gansu Province, College of Chemical Engineering, Gansu Provincial Biomass Function Composites Engineering Research Center, Northwest Minzu University, Gansu, 730030, Lanzhou, China
| | - Su Qiong
- Key Laboratory of Environment-Friendly Composite Materials of the State Ethnic Affairs Commission, Key Laboratory for Utility of Environment-Friendly Composite Materials and Biomass in University of Gansu Province, College of Chemical Engineering, Gansu Provincial Biomass Function Composites Engineering Research Center, Northwest Minzu University, Gansu, 730030, Lanzhou, China.
| |
Collapse
|
10
|
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
|
11
|
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
|
12
|
Junxi L, Fupeng Z, Bomiao Q, Weimin J, Hongqiang L, Qiong S. Reaction of CHCl•- with HCHO and H2O: A theoretical study. COMPUT THEOR CHEM 2022. [DOI: 10.1016/j.comptc.2022.113932] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/06/2022]
|
13
|
Li J, Yang W. Renormalized Singles with Correlation in GW Green's Function Theory for Accurate Quasiparticle Energies. J Phys Chem Lett 2022; 13:9372-9380. [PMID: 36190273 DOI: 10.1021/acs.jpclett.2c02051] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/16/2023]
Abstract
We apply the renormalized singles with the correlation (RSc) Green function in the GW approximation for accurate quasiparticle (QP) energies and orbitals. The RSc Green function includes singles contributions from the associated density functional approximation (DFA) and considers correlation contributions perturbatively. GRScWRSc uses the RSc Green function as the new starting point and in the formulation of the screened interaction. GRScW0 fixes the screened interaction at the DFA level. For the calculations of ionization potentials, GRScWRSc and GRScW0 significantly reduce the starting point dependence and provide accurate results with errors around 0.2 eV. For the calculations of core-level binding energies, GRScWRSc slightly overestimates the results because of underscreening, but GRScW0 with GGA functionals provides the optimal accuracy with errors of 0.40 eV. We also show that GRScWRSc predicts accurate dipole moments. GRScWRSc and GRScW0, are computationally favorable compared with any self-consistent GW methods. The RSc approach is promising for making GW and other Green function methods efficient and robust.
Collapse
Affiliation(s)
- Jiachen Li
- Department of Chemistry, Duke University, Durham, North Carolina27708, United States
| | - Weitao Yang
- Department of Chemistry, Duke University, Durham, North Carolina27708, United States
| |
Collapse
|