1
|
Kato D, Suzuki H, Abe R, Kageyama H. Band engineering of layered oxyhalide photocatalysts for visible-light water splitting. Chem Sci 2024; 15:11719-11736. [PMID: 39092126 PMCID: PMC11290441 DOI: 10.1039/d4sc02093f] [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: 03/29/2024] [Accepted: 06/25/2024] [Indexed: 08/04/2024] Open
Abstract
The band structure offers fundamental information on electronic properties of solid state materials, and hence it is crucial for solid state chemists to understand and predict the relationship between the band structure and electronic structure to design chemical and physical properties. Here, we review layered oxyhalide photocatalysts for water splitting with a particular emphasis on band structure control. The unique feature of these materials including Sillén and Sillén-Aurivillius oxyhalides lies in their band structure including a remarkably high oxygen band, allowing them to exhibit both visible light responsiveness and photocatalytic stability unlike conventional mixed anion compounds, which show good light absorption, but frequently encounter stability issues. For band structure control, simple strategies effective in mixed-anion compounds, such as anion substitution forming high energy p orbitals in accordance with its electronegativity, is not effective for oxyhalides with high oxygen bands. We overview key concepts for band structure control of oxyhalide photocatalysts such as lone-pair interactions and electrostatic interactions. The control of the band structure of inorganic solid materials is a crucial challenge across a wide range of materials chemistry fields, and the insights obtained by the development of oxyhalide photocatalysts are expected to provide knowledge for diverse materials chemistry.
Collapse
Affiliation(s)
- Daichi Kato
- Department of Energy and Hydrocarbon Chemistry, Graduate School of Engineering, Kyoto University Nishikyo-ku Kyoto 615-8510 Japan
| | - Hajime Suzuki
- Department of Energy and Hydrocarbon Chemistry, Graduate School of Engineering, Kyoto University Nishikyo-ku Kyoto 615-8510 Japan
| | - Ryu Abe
- Department of Energy and Hydrocarbon Chemistry, Graduate School of Engineering, Kyoto University Nishikyo-ku Kyoto 615-8510 Japan
| | - Hiroshi Kageyama
- Department of Energy and Hydrocarbon Chemistry, Graduate School of Engineering, Kyoto University Nishikyo-ku Kyoto 615-8510 Japan
| |
Collapse
|
2
|
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
|
3
|
Vacondio S, Varsano D, Ruini A, Ferretti A. Going Beyond the GW Approximation Using the Time-Dependent Hartree-Fock Vertex. J Chem Theory Comput 2024; 20:4718-4737. [PMID: 38772396 DOI: 10.1021/acs.jctc.4c00100] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/23/2024]
Abstract
The time-dependent Hartree-Fock (TDHF) vertex of many-body perturbation theory (MBPT) makes it possible to extend TDHF theory to charged excitations. Here we assess its performance by applying it to spherical atoms in their neutral electronic configuration. On a theoretical level, we recast the TDHF vertex as a reducible vertex, highlighting the emergence of a self-energy expansion purely in orders of the bare Coulomb interaction; then, on a numerical level, we present results for polarizabilities, ionization energies (IEs), and photoemission satellites. We confirm the superiority of THDF over simpler methods such as the random phase approximation for the prediction of atomic polarizabilities. We then find that the TDHF vertex reliably provides better IEs than GW and low-order self-energies do in the light-atom, few-electron regime; its performance degrades in heavier, many-electron atoms instead, where an expansion in orders of an unscreened Coulomb interaction becomes less justified. New relevant features are introduced in the satellite spectrum by the TDHF vertex, but the experimental spectra are not fully reproduced due to a missing account of nonlinear effects connected to hole relaxation. We also explore various truncations of the self-energy given by the TDHF vertex, but do not find them to be more convenient than low-order approximations such as GW and second Born (2B), suggesting that vertex corrections should be carried out consistently both in the self-energy and in the polarizability.
Collapse
Affiliation(s)
- Simone Vacondio
- Dipartimento di Scienze Fisiche, Informatiche e Matematiche, Università degli Studi di Modena e Reggio Emilia, Via G. Campi 213/a, 41125 Modena, Italy
- Centro S3, CNR-Istituto Nanoscienze, Via G. Campi 213/a, 41125 Modena, Italy
| | - Daniele Varsano
- Centro S3, CNR-Istituto Nanoscienze, Via G. Campi 213/a, 41125 Modena, Italy
| | - Alice Ruini
- Dipartimento di Scienze Fisiche, Informatiche e Matematiche, Università degli Studi di Modena e Reggio Emilia, Via G. Campi 213/a, 41125 Modena, Italy
- Centro S3, CNR-Istituto Nanoscienze, Via G. Campi 213/a, 41125 Modena, Italy
| | - Andrea Ferretti
- Centro S3, CNR-Istituto Nanoscienze, Via G. Campi 213/a, 41125 Modena, Italy
| |
Collapse
|
4
|
Krumland J, Cocchi C. Ab Initio Modeling of Mixed-Dimensional Heterostructures: A Path Forward. J Phys Chem Lett 2024; 15:5350-5358. [PMID: 38728611 PMCID: PMC11129309 DOI: 10.1021/acs.jpclett.4c00803] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/15/2024] [Revised: 04/26/2024] [Accepted: 05/08/2024] [Indexed: 05/12/2024]
Abstract
Understanding the electronic structure of mixed-dimensional heterostructures is essential for maximizing their application potential. However, accurately modeling such interfaces is challenging due to the complex interplay between the subsystems. We employ a computational framework integrating first-principles methods, including GW, density functional theory (DFT), and the polarizable continuum model, to elucidate the electronic structure of mixed-dimensional heterojunctions formed by free-base phthalocyanines and monolayer molybdenum disulfide. We assess the impact of dielectric screening across various scenarios, from isolated molecules to organic films on a substrate-supported monolayer. Our findings show that while polarization effects cause significant renormalization of molecular energy levels, band energies and alignments in the most relevant setup can be accurately predicted through DFT simulations of the individual subsystems. Additionally, we analyze orbital hybridization, revealing potential pathways for interfacial charge transfer. This study offers new insights into hybrid inorganic/organic interfaces and provides a practical computational protocol suitable for scaled-up studies.
Collapse
Affiliation(s)
- Jannis Krumland
- Institute
of Physics, Carl von Ossietzky Universität
Oldenburg, 26129 Oldenburg, Germany
- Physics
Department and IRIS Adlershof, Humboldt-Universität
zu Berlin, 12489 Berlin, Germany
| | - Caterina Cocchi
- Institute
of Physics, Carl von Ossietzky Universität
Oldenburg, 26129 Oldenburg, Germany
- Physics
Department and IRIS Adlershof, Humboldt-Universität
zu Berlin, 12489 Berlin, Germany
| |
Collapse
|
5
|
Wang X, Gao S, Luo Y, Liu X, Tom R, Zhao K, Chang V, Marom N. Computational Discovery of Intermolecular Singlet Fission Materials Using Many-Body Perturbation Theory. THE JOURNAL OF PHYSICAL CHEMISTRY. C, NANOMATERIALS AND INTERFACES 2024; 128:7841-7864. [PMID: 38774154 PMCID: PMC11103713 DOI: 10.1021/acs.jpcc.4c01340] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 02/29/2024] [Revised: 04/16/2024] [Accepted: 04/17/2024] [Indexed: 05/24/2024]
Abstract
Intermolecular singlet fission (SF) is the conversion of a photogenerated singlet exciton into two triplet excitons residing on different molecules. SF has the potential to enhance the conversion efficiency of solar cells by harvesting two charge carriers from one high-energy photon, whose surplus energy would otherwise be lost to heat. The development of commercial SF-augmented modules is hindered by the limited selection of molecular crystals that exhibit intermolecular SF in the solid state. Computational exploration may accelerate the discovery of new SF materials. The GW approximation and Bethe-Salpeter equation (GW+BSE) within the framework of many-body perturbation theory is the current state-of-the-art method for calculating the excited-state properties of molecular crystals with periodic boundary conditions. In this Review, we discuss the usage of GW+BSE to assess candidate SF materials as well as its combination with low-cost physical or machine learned models in materials discovery workflows. We demonstrate three successful strategies for the discovery of new SF materials: (i) functionalization of known materials to tune their properties, (ii) finding potential polymorphs with improved crystal packing, and (iii) exploring new classes of materials. In addition, three new candidate SF materials are proposed here, which have not been published previously.
Collapse
Affiliation(s)
- Xiaopeng Wang
- School
of Foundational Education, University of
Health and Rehabilitation Sciences, Qingdao 266113, China
- Qingdao
Institute for Theoretical and Computational Sciences, Institute of
Frontier and Interdisciplinary Science, Shandong University, Qingdao, Shandong 266237, P. R. China
| | - Siyu Gao
- Department
of Materials Science and Engineering, Carnegie
Mellon University, Pittsburgh, Pennsylvania 15213, United States
| | - Yiqun Luo
- Department
of Physics, Carnegie Mellon University, Pittsburgh, Pennsylvania 15213, United States
| | - Xingyu Liu
- Department
of Materials Science and Engineering, Carnegie
Mellon University, Pittsburgh, Pennsylvania 15213, United States
| | - Rithwik Tom
- Department
of Physics, Carnegie Mellon University, Pittsburgh, Pennsylvania 15213, United States
| | - Kaiji Zhao
- Department
of Materials Science and Engineering, Carnegie
Mellon University, Pittsburgh, Pennsylvania 15213, United States
| | - Vincent Chang
- Department
of Materials Science and Engineering, Carnegie
Mellon University, Pittsburgh, Pennsylvania 15213, United States
| | - Noa Marom
- Department
of Materials Science and Engineering, Carnegie
Mellon University, Pittsburgh, Pennsylvania 15213, United States
- Department
of Physics, Carnegie Mellon University, Pittsburgh, Pennsylvania 15213, United States
- Department
of Chemistry, Carnegie Mellon University, Pittsburgh, Pennsylvania 15213, United States
| |
Collapse
|
6
|
Yeh CN, Morales MA. Low-Scaling Algorithms for GW and Constrained Random Phase Approximation Using Symmetry-Adapted Interpolative Separable Density Fitting. J Chem Theory Comput 2024; 20:3184-3198. [PMID: 38597496 DOI: 10.1021/acs.jctc.4c00085] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 04/11/2024]
Abstract
We present low-scaling algorithms for GW and constrained random phase approximation based on a symmetry-adapted interpolative separable density fitting (ISDF) procedure that incorporates the space-group symmetries of crystalline systems. The resulting formulations scale cubically, with respect to system size, and linearly with the number of k-points, regardless of the choice of single-particle basis and whether a quasiparticle approximation is employed. We validate these methods through comparisons with published literature and demonstrate their efficiency in treating large-scale systems through the construction of downfolded many-body Hamiltonians for carbon dimer defects embedded in hexagonal boron nitride supercells. Our work highlights the efficiency and general applicability of ISDF in the context of large-scale many-body calculations with k-point sampling beyond density functional theory.
Collapse
Affiliation(s)
- Chia-Nan Yeh
- Center for Computational Quantum Physics, Flatiron Institute, New York, New York 10010, United States
| | - Miguel A Morales
- Center for Computational Quantum Physics, Flatiron Institute, New York, New York 10010, United States
| |
Collapse
|
7
|
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
|
8
|
Wen M, Abraham V, Harsha G, Shee A, Whaley KB, Zgid D. Comparing Self-Consistent GW and Vertex-Corrected G0W0 ( G0W0Γ) Accuracy for Molecular Ionization Potentials. J Chem Theory Comput 2024; 20:3109-3120. [PMID: 38573104 DOI: 10.1021/acs.jctc.3c01279] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 04/05/2024]
Abstract
We test the performance of self-consistent GW and several representative implementations of vertex-corrected G0W0 (G0W0Γ). These approaches are tested on benchmark data sets covering full valence spectra (first ionization potentials and some inner valence shell excitations). For small molecules, when comparing against state-of-the-art wave function techniques, our results show that full self-consistency in the GW scheme either systematically outperforms vertex-corrected G0W0 or gives results of at least comparative quality. Moreover, G0W0Γ results in additional computational cost when compared to G0W0 or self-consistent GW. The dependency of G0W0Γ on the starting mean-field solution is frequently more dominant than the magnitude of the vertex correction itself. Consequently, for molecular systems, self-consistent GW performed on the imaginary axis (and then followed by modern analytical continuation techniques) offers a more reliable approach to make predictions of ionization potentials.
Collapse
Affiliation(s)
- Ming Wen
- Department of Chemistry, University of Michigan, Ann Arbor, Michigan 48109, United States
| | - 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
| | - Avijit Shee
- Department of Chemistry, University of Michigan, Ann Arbor, Michigan 48109, United States
| | - K Birgitta Whaley
- Department of Chemistry, University of California, Berkeley, California 94720-1460, United States
| | - Dominika Zgid
- Department of Chemistry, University of Michigan, Ann Arbor, Michigan 48109, United States
- Department of Physics, University of Michigan, Ann Arbor, Michigan 48109, United States
| |
Collapse
|
9
|
Hu R, Ho DQ, To DQ, Bryant GW, Janotti A. Fermi-Level Pinning in ErAs Nanoparticles Embedded in III-V Semiconductors. NANO LETTERS 2024; 24:4376-4382. [PMID: 38591335 DOI: 10.1021/acs.nanolett.3c04995] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 04/10/2024]
Abstract
Embedding rare-earth monopnictide nanoparticles into III-V semiconductors enables unique optical, electrical, and thermal properties for THz photoconductive switches, tunnel junctions, and thermoelectric devices. Despite the high structural quality and control over growth, particle size (<3 nm), and density, the underlying electronic structure of these nanocomposite materials has only been hypothesized. Structural and electronic properties of ErAs nanoparticles with different shapes and sizes (cubic to spherical, 1.14, 1.71, and 2.28 nm) in AlAs, GaAs, InAs, and their alloys are investigated using first-principles calculations, revealing that spherical nanoparticles have lower formation energies. For the lowest-energy nanoparticles, the Fermi level is pinned near midgap in GaAs and AlAs but resonant in the conduction band in InAs. The Fermi level is shifted down as the particle size increases and is pinned on an absolute energy scale considering the band alignment at AlAs/GaAs/InAs interfaces, offering insights into the rational design of these nanomaterials.
Collapse
Affiliation(s)
- Ruiqi Hu
- Department of Materials Science and Engineering, University of Delaware, Newark, Delaware 19716, United States
| | - Dai Q Ho
- Department of Materials Science and Engineering, University of Delaware, Newark, Delaware 19716, United States
- Faculty of Natural Sciences, Quy Nhon University, Quy Nhon 551130, Vietnam
| | - D Quang To
- Department of Materials Science and Engineering, University of Delaware, Newark, Delaware 19716, United States
| | - Garnett W Bryant
- Nanoscale Device Characterization Division, Joint Quantum Institute, National Institute of Standards and Technology, Gaithersburg, Maryland 20899-8423, United States
- University of Maryland, College Park, Maryland 20742, United States
| | - Anderson Janotti
- Department of Materials Science and Engineering, University of Delaware, Newark, Delaware 19716, United States
| |
Collapse
|
10
|
Kiyohara S, Hinuma Y, Oba F. Band Alignment of Oxides by Learnable Structural-Descriptor-Aided Neural Network and Transfer Learning. J Am Chem Soc 2024; 146:9697-9708. [PMID: 38546127 PMCID: PMC11009958 DOI: 10.1021/jacs.3c13574] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/02/2023] [Revised: 02/25/2024] [Accepted: 02/27/2024] [Indexed: 04/11/2024]
Abstract
The band alignment of semiconductors, insulators, and dielectrics is relevant to diverse material properties and device structures utilizing their surfaces and interfaces. In particular, the ionization potential and electron affinity are fundamental quantities that describe surface-dependent band-edge positions with respect to the vacuum level. Their accurate and systematic determination, however, demands elaborate experiments or simulations for well-characterized surfaces. Here, we report machine learning for the band alignment of nonmetallic oxides using a high-throughput first-principles calculation data set containing about 3000 oxide surfaces. Our neural network accurately predicts the band positions for relaxed surfaces of binary oxides simply by using the information on bulk structures and surface termination planes. Moreover, we extend the model to naturally include multiple-cation effects and transfer it to ternary oxides. The present approach enables the band alignment of a vast number of solid surfaces, thereby opening the way to a systematic understanding and materials screening.
Collapse
Affiliation(s)
- Shin Kiyohara
- Laboratory
for Materials and Structures, Institute of Innovative Research, Tokyo Institute of Technology, R3-7, 4259 Nagatsuta, Midori-ku, Yokohama 226-8501, Japan
- Institute
for Materials Research, Tohoku University, 2-2-1 Katahira,
Aoba-ku, Sendai 980-8577, Japan
| | - Yoyo Hinuma
- Department
of Energy and Environment, National Institute
of Advanced Industrial Science and Technology (AIST), 1-8-31 Midorigaoka, Ikeda, Osaka 563-8577, Japan
| | - Fumiyasu Oba
- Laboratory
for Materials and Structures, Institute of Innovative Research, Tokyo Institute of Technology, R3-7, 4259 Nagatsuta, Midori-ku, Yokohama 226-8501, Japan
- MDX
Research Center for Element Strategy, International Research Frontiers
Initiative, Tokyo Institute of Technology, SE-6, 4259 Nagatsuta, Midori-ku, Yokohama 226-8501, Japan
| |
Collapse
|
11
|
El-Sahili A, Sottile F, Reining L. Total Energy beyond GW: Exact Results and Guidelines for Approximations. J Chem Theory Comput 2024; 20:1972-1987. [PMID: 38324673 DOI: 10.1021/acs.jctc.3c01200] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/09/2024]
Abstract
The total energy and electron addition and removal spectra can, in principle, be obtained exactly from the one-body Green's function (GF). In practice, the GF is obtained from an approximate self-energy. In the framework of many-body perturbation theory, we derive different expressions that are based on an approximate self-energy, but that yield nevertheless, in principle, the exact exchange-correlation contribution to the total energy for any interaction strength. Response functions play a crucial role, which explains why, for example, ingredients from time-dependent density functional theory can be used to build these approximate self-energies. We show that the key requirement for obtaining exact results is the consistent combination of ingredients. Also when further approximations are made, as it is necessary in practice, this consistency remains the key to obtain good results. All findings are illustrated using the exactly solvable symmetric Hubbard dimer.
Collapse
Affiliation(s)
- Abdallah El-Sahili
- LSI, CNRS, CEA/DRF/IRAMIS, École Polytechnique, Institut Polytechnique de Paris, Palaiseau F-91120, France
- European Theoretical Spectroscopy Facility (ETSF), https://www.etsf.eu/
| | - Francesco Sottile
- LSI, CNRS, CEA/DRF/IRAMIS, École Polytechnique, Institut Polytechnique de Paris, Palaiseau F-91120, France
- European Theoretical Spectroscopy Facility (ETSF), https://www.etsf.eu/
| | - Lucia Reining
- LSI, CNRS, CEA/DRF/IRAMIS, École Polytechnique, Institut Polytechnique de Paris, Palaiseau F-91120, France
- European Theoretical Spectroscopy Facility (ETSF), https://www.etsf.eu/
| |
Collapse
|
12
|
Tal A, Bischoff T, Pasquarello A. Absolute energy levels of liquid water from many-body perturbation theory with effective vertex corrections. Proc Natl Acad Sci U S A 2024; 121:e2311472121. [PMID: 38427604 DOI: 10.1073/pnas.2311472121] [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/11/2023] [Accepted: 01/31/2024] [Indexed: 03/03/2024] Open
Abstract
We demonstrate the importance of addressing the Γ vertex and thus going beyond the GW approximation for achieving the energy levels of liquid water in many-body perturbation theory. In particular, we consider an effective vertex function in both the polarizability and the self-energy, which does not produce any computational overhead compared with the GW approximation. We yield the band gap, the ionization potential, and the electron affinity in good agreement with experiment and with a hybrid functional description. The achieved electronic structure and dielectric screening further lead to a good description of the optical absorption spectrum, as obtained through the solution of the Bethe-Salpeter equation. In particular, the experimental peak position of the exciton is accurately reproduced.
Collapse
Affiliation(s)
- Alexey Tal
- Chaire de Simulation à l'Echelle Atomique, Ecole Polytechnique Fédérale de Lausanne, Lausanne CH-1015, Switzerland
| | - Thomas Bischoff
- Chaire de Simulation à l'Echelle Atomique, Ecole Polytechnique Fédérale de Lausanne, Lausanne CH-1015, Switzerland
| | - Alfredo Pasquarello
- Chaire de Simulation à l'Echelle Atomique, Ecole Polytechnique Fédérale de Lausanne, Lausanne CH-1015, Switzerland
| |
Collapse
|
13
|
Orlando R, Romaniello P, Loos PF. The three channels of many-body perturbation theory: GW, particle-particle, and electron-hole T-matrix self-energies. J Chem Phys 2023; 159:184113. [PMID: 37962450 DOI: 10.1063/5.0176898] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/18/2023] [Accepted: 10/23/2023] [Indexed: 11/15/2023] Open
Abstract
We derive the explicit expression of the three self-energies that one encounters in many-body perturbation theory: the well-known GW self-energy, as well as the particle-particle and electron-hole T-matrix self-energies. Each of these can be easily computed via the eigenvalues and eigenvectors of a different random-phase approximation linear eigenvalue problem that completely defines their corresponding response function. For illustrative and comparative purposes, we report the principal ionization potentials of a set of small molecules computed at each level of theory. The performance of these schemes on strongly correlated systems (B2 and C2) is also discussed.
Collapse
Affiliation(s)
- Roberto Orlando
- Laboratoire de Chimie et Physique Quantiques, Université de Toulouse, CNRS, UPS, Toulouse, France
- Laboratoire de Physique Théorique, Université de Toulouse, CNRS, UPS, Toulouse, France
- European Theoretical Spectroscopy Facility (ETSF)
| | - Pina Romaniello
- Laboratoire de Physique Théorique, Université de Toulouse, CNRS, UPS, Toulouse, France
- European Theoretical Spectroscopy Facility (ETSF)
| | - Pierre-François Loos
- Laboratoire de Chimie et Physique Quantiques, Université de Toulouse, CNRS, UPS, Toulouse, France
| |
Collapse
|
14
|
Förster A, van Lenthe E, Spadetto E, Visscher L. Two-Component GW Calculations: Cubic Scaling Implementation and Comparison of Vertex-Corrected and Partially Self-Consistent GW Variants. J Chem Theory Comput 2023; 19:5958-5976. [PMID: 37594901 PMCID: PMC10501001 DOI: 10.1021/acs.jctc.3c00512] [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/17/2023] [Indexed: 08/20/2023]
Abstract
We report an all-electron, atomic orbital (AO)-based, two-component (2C) implementation of the GW approximation (GWA) for closed-shell molecules. Our algorithm is based on the space-time formulation of the GWA and uses analytical continuation (AC) of the self-energy, and pair-atomic density fitting (PADF) to switch between AO and auxiliary basis. By calculating the dynamical contribution to the GW self-energy at a quasi-one-component level, our 2C-GW algorithm is only about a factor of 2-3 slower than in the scalar relativistic case. Additionally, we present a 2C implementation of the simplest vertex correction to the self-energy, the statically screened G3W2 correction. Comparison of first ionization potentials (IPs) of a set of 67 molecules with heavy elements (a subset of the SOC81 set) calculated with our implementation against results from the WEST code reveals mean absolute deviations (MAD) of around 70 meV for G0W0@PBE and G0W0@PBE0. We check the accuracy of our AC treatment by comparison to full-frequency GW calculations, which shows that in the absence of multisolution cases, the errors due to AC are only minor. This implies that the main sources of the observed deviations between both implementations are the different single-particle bases and the pseudopotential approximation in the WEST code. Finally, we assess the performance of some (partially self-consistent) variants of the GWA for the calculation of first IPs by comparison to vertical experimental reference values. G0W0@PBE0 (25% exact exchange) and G0W0@BHLYP (50% exact exchange) perform best with mean absolute deviations (MAD) of about 200 meV. Explicit treatment of spin-orbit effects at the 2C level is crucial for systematic agreement with experiment. On the other hand, eigenvalue-only self-consistent GW (evGW) and quasi-particle self-consistent GW (qsGW) significantly overestimate the IPs. Perturbative G3W2 corrections increase the IPs and therefore improve the agreement with experiment in cases where G0W0 alone underestimates the IPs. With a MAD of only 140 meV, 2C-G0W0@PBE0 + G3W2 is in best agreement with the experimental reference values.
Collapse
Affiliation(s)
- Arno Förster
- Theoretical
Chemistry, Vrije Universiteit, De Boelelaan 1083, 1081 HV Amsterdam, The Netherlands
| | - Erik van Lenthe
- Software
for Chemistry and Materials NV, 1081 HV Amsterdam, The Netherlands
| | - Edoardo Spadetto
- Software
for Chemistry and Materials NV, 1081 HV Amsterdam, The Netherlands
| | - Lucas Visscher
- Theoretical
Chemistry, Vrije Universiteit, De Boelelaan 1083, 1081 HV Amsterdam, The Netherlands
| |
Collapse
|
15
|
Yasumura S, Kamachi T, Toyao T, Shimizu KI, Hinuma Y. Prediction of Stable Surfaces of Metal Oxides through the Unsaturated Coordination Index. ACS OMEGA 2023; 8:29779-29788. [PMID: 37599947 PMCID: PMC10433516 DOI: 10.1021/acsomega.3c04253] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 06/15/2023] [Accepted: 07/10/2023] [Indexed: 08/22/2023]
Abstract
This study proposes the unsaturated coordination index, σ, as a potential descriptor of the stability of metal-oxide surfaces cleaved from bulk. The value of σ, the number of missing bonds per unit area, can be obtained very quickly using only crystallographic data, namely, the bulk geometry. The surface energies of various binary oxides, with and without atom relaxation, were calculated. Their correlations with σ had good coefficients of determination (R2) values, particularly in high-symmetry crystals. The proposed descriptor is very useful for an initial evaluation of stable metal-oxide surfaces without conducting any surface model calculations.
Collapse
Affiliation(s)
- Shunsaku Yasumura
- Institute
of Industrial Science, The University of
Tokyo, Komaba 4-6-1, Meguro, Tokyo 153-8505, Japan
| | - Takashi Kamachi
- Department
of Life, Environment and Applied Chemistry, Fukuoka Institute of Technology, 3-30-1 Wajiro-Higashi, Higashi-ku, Fukuoka 811-0295, Japan
| | - Takashi Toyao
- Institute
for Catalysis, Hokkaido University, N-21, W-10, Kita, Sapporo 001-0021, Hokkaido, Japan
| | - Ken-ichi Shimizu
- Institute
for Catalysis, Hokkaido University, N-21, W-10, Kita, Sapporo 001-0021, Hokkaido, Japan
| | - Yoyo Hinuma
- Department
of Energy and Environment, National Institute
of Advanced Industrial Science and Technology (AIST), 1-8-31, Midorigaoka, Ikeda 563-8577, Osaka, Japan
| |
Collapse
|
16
|
Marie A, Loos PF. A Similarity Renormalization Group Approach to Green's Function Methods. J Chem Theory Comput 2023. [PMID: 37311565 DOI: 10.1021/acs.jctc.3c00281] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/15/2023]
Abstract
The family of Green's function methods based on the GW approximation has gained popularity in the electronic structure theory thanks to its accuracy in weakly correlated systems combined with its cost-effectiveness. Despite this, self-consistent versions still pose challenges in terms of convergence. A recent study [Monino and Loos J. Chem. Phys. 2022, 156, 231101.] has linked these convergence issues to the intruder-state problem. In this work, a perturbative analysis of the similarity renormalization group (SRG) approach is performed on Green's function methods. The SRG formalism enables us to derive, from first-principles, the expression of a naturally static and Hermitian form of the self-energy that can be employed in quasiparticle self-consistent GW (qsGW) calculations. The resulting SRG-based regularized self-energy significantly accelerates the convergence of qsGW calculations, slightly improves the overall accuracy, and is straightforward to implement in existing code.
Collapse
Affiliation(s)
- Antoine Marie
- Laboratoire de Chimie et Physique Quantiques (UMR 5626), Université de Toulouse, CNRS, UPS, 31400 Toulouse, France
| | - Pierre-François Loos
- Laboratoire de Chimie et Physique Quantiques (UMR 5626), Université de Toulouse, CNRS, UPS, 31400 Toulouse, France
| |
Collapse
|
17
|
Reshetnyak I, Lorin A, Pasquarello A. Many-body screening effects in liquid water. Nat Commun 2023; 14:2705. [PMID: 37169764 PMCID: PMC10175292 DOI: 10.1038/s41467-023-38420-w] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/30/2020] [Accepted: 04/27/2023] [Indexed: 05/13/2023] Open
Abstract
The screening arising from many-body excitations is a crucial quantity for describing absorption and inelastic X-ray scattering (IXS) of materials. Similarly, the electron screening plays a critical role in state-of-the-art approaches for determining the fundamental band gap. However, ab initio studies of the screening in liquid water have remained limited. Here, we use a combined analysis based on the Bethe-Salpeter equation and time-dependent density functional theory. We first show that absorption spectra at near-edge energies are insufficient to assess the accuracy by which the screening is described. Next, when the energy range under scrutiny is extended, we instead find that the IXS spectra are highly sensitive and allow for the selection of the optimal theoretical scheme. This leads to good agreement with experiment over a large range of transferred energies and momenta, and enables establishing the elusive fundamental band gap of liquid water at 9.3 eV.
Collapse
Affiliation(s)
- Igor Reshetnyak
- Chaire de Simulation à l'Echelle Atomique (CSEA), Ecole Polytechnique Fédérale de Lausanne (EPFL), CH-1015, Lausanne, Switzerland
| | - Arnaud Lorin
- Chaire de Simulation à l'Echelle Atomique (CSEA), Ecole Polytechnique Fédérale de Lausanne (EPFL), CH-1015, Lausanne, Switzerland
| | - Alfredo Pasquarello
- Chaire de Simulation à l'Echelle Atomique (CSEA), Ecole Polytechnique Fédérale de Lausanne (EPFL), CH-1015, Lausanne, Switzerland.
| |
Collapse
|
18
|
Weng G, Mallarapu R, Vlček V. Embedding vertex corrections in GW self-energy: Theory, implementation, and outlook. J Chem Phys 2023; 158:144105. [PMID: 37061461 DOI: 10.1063/5.0139117] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 04/17/2023] Open
Abstract
The vertex function (Γ) within the Green's function formalism encapsulates information about all higher-order electron-electron interaction beyond those mediated by density fluctuations. Herein, we present an efficient approach that embeds vertex corrections in the one-shot GW correlation self-energy for isolated and periodic systems. The vertex-corrected self-energy is constructed through the proposed separation-propagation-recombination procedure: the electronic Hilbert space is separated into an active space and its orthogonal complement denoted as the "rest;" the active component is propagated by a space-specific effective Hamiltonian different from the rest. The vertex corrections are introduced by a rescaled time-dependent nonlocal exchange interaction. The direct Γ correction to the self-energy is further updated by adjusting the rescaling factor in a self-consistent post-processing cycle. Our embedding method is tested mainly on donor-acceptor charge-transfer systems. The embedded vertex effects consistently and significantly correct the quasiparticle energies of the gap-edge states. The fundamental gap is generally improved by 1-3 eV upon the one-shot GW approximation. Furthermore, we provide an outlook for applications of (embedded) vertex corrections in calculations of extended solids.
Collapse
Affiliation(s)
- Guorong Weng
- Department of Chemistry and Biochemistry, University of California, Santa Barbara, California 93106-9510, USA
| | - Rushil Mallarapu
- Department of Chemistry and Biochemistry, University of California, Santa Barbara, California 93106-9510, USA
| | - Vojtěch Vlček
- Department of Chemistry and Biochemistry, University of California, Santa Barbara, California 93106-9510, USA
| |
Collapse
|
19
|
Monzel L, Holzer C, Klopper W. Natural virtual orbitals for the GW method in the random-phase approximation and beyond. J Chem Phys 2023; 158:144102. [PMID: 37061489 DOI: 10.1063/5.0144469] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 04/17/2023] Open
Abstract
The increasingly popular GW method is becoming a convenient tool to determine vertical ionization energies in molecular systems. However, depending on the formalism used and the range of orbitals investigated, it may be hampered by a steep computational scaling. To alleviate this issue, correlated natural virtual orbitals (NVOs) based on second-order Møller-Plesset (MP2) and direct MP2 correlation energies are implemented, and the resulting correlated NVOs are tested on GW quasiparticle energies. Test cases include the popular GW variants G0W0 and evGW0 as well as more elaborate vertex corrections. We find that for increasingly larger molecular systems and basis sets, NVOs considerably improve efficiency. Furthermore, we test the performance of the truncated (frozen) NVO ansatz on the GW100 test set. For the latter, it is demonstrated that, using a carefully chosen truncation threshold, NVOs lead to a negligible loss in accuracy while providing speedups of one order of magnitude. Furthermore, we compare the resulting quasiparticle energies to very accurate vertical ionization energies obtained from coupled-cluster theory with singles, doubles, and noniterative triples [CCSD(T)], confirming that the loss in accuracy introduced by truncating the NVOs is negligible compared to the methodical errors in the GW approximation. It is also demonstrated that the choice of basis set impacts results far more than using a suitably truncated NVO space. Therefore, at the same computational expense, more accurate results can be obtained using NVOs. Finally, we provide improved reference CCSD(T) values for the GW100 test set, which have been obtained using the def2-QZVPP basis set.
Collapse
Affiliation(s)
- Laurenz Monzel
- Institute of Physical Chemistry, Karlsruhe Institute of Technology (KIT), Fritz-Haber-Weg 2, 76131 Karlsruhe, Germany
| | - Christof Holzer
- Institute of Theoretical Solid State Physics, Karlsruhe Institute of Technology (KIT), Wolfgang-Gaede-Straße 1, 76131 Karlsruhe, Germany
| | - Wim Klopper
- Institute of Physical Chemistry, Karlsruhe Institute of Technology (KIT), Fritz-Haber-Weg 2, 76131 Karlsruhe, Germany
| |
Collapse
|
20
|
Chatratin I, Dou B, Wei SH, Janotti A. Doping Limits of Phosphorus, Arsenic, and Antimony in CdTe. J Phys Chem Lett 2023; 14:273-278. [PMID: 36595563 DOI: 10.1021/acs.jpclett.2c03233] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/17/2023]
Abstract
Low p-type doping is a limiting factor to increase CdTe thin-film solar-cell efficiency toward the theoretical Shockley-Queisser limit of 33%. Previous calculations predict relatively high ionization energies for group-V acceptors (P, As, and Sb), and they are plagued by self-compensation, forming AX centers, severely limiting hole concentration. However, recent experiments on CdTe single crystals indicate a much more favorable scenario, where P, As, and Sb behave as shallow acceptors. Using hybrid functional calculations, we solve this puzzle by showing that the ionization energies significantly decrease with the supercell size. When including the effects of spin-orbit coupling and extrapolating the results to the dilute limit, we find these impurities behave as hydrogenic-like shallow acceptors, and AX centers are unstable and do not limit p-type doping. We address the differences between our results and previous theoretical predictions and show that our ionization energies predict hole concentrations that agree with recent temperature-dependent Hall measurements.
Collapse
Affiliation(s)
- Intuon Chatratin
- Department of Materials Science & Engineering, University of Delaware, Newark, Delaware19716, United States
| | - Baoying Dou
- Beijing Computational Science Research Center, Beijing100193, China
| | - Su-Huai Wei
- Beijing Computational Science Research Center, Beijing100193, China
| | - Anderson Janotti
- Department of Materials Science & Engineering, University of Delaware, Newark, Delaware19716, United States
| |
Collapse
|
21
|
Li W, Cai XF, Valdes N, Wang T, Shafarman W, Wei SH, Janotti A. In 2Se 3, In 2Te 3, and In 2(Se,Te) 3 Alloys as Photovoltaic Materials. J Phys Chem Lett 2022; 13:12026-12031. [PMID: 36541824 DOI: 10.1021/acs.jpclett.2c02975] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/17/2023]
Abstract
In its lowest-energy three-dimensional (3D) hexagonal crystal structure (γ phase), In2Se3 has a direct band gap of ∼1.8 eV and displays high absorption coefficient, making it a promising semiconductor material for optoelectronics. Incorporation of Te allows for tuning the band gap, adding flexibility to device design and extending the application range. Here we report results of hybrid density functional theory calculations to assess the electronic and optical properties of γ-In2Se3, γ-In2Te3, and γ-In2(Se1-xTex)3 alloys, and initial experiments on the growth and characterization of γ-In2Se3 thin films. The predicted band gap of 1.84 eV for γ-In2Se3 is in good agreement with the absorption onset derived from transmission and reflection spectra of thin films. We show that incorporation of Te gives γ-In2(Se1-xTex)3 alloys with a band gap ranging from 1.84 eV down to 1.23 eV, thus covering the optimal band gap range for single-junction solar cells. In addition, the γ-In2Se3/γ-In2(Se1-xTex)3 bilayer could be employed in tandem solar-cell architectures absorbing at Eg ≈ 1.8 eV and at Eg ≤ 1.4 eV, toward overcoming the ∼33% efficiency set by the Shockley-Queisser limit for single junction solar cells. We also discuss band gap bowing and mixing enthalpies, aiming at adding γ-In2Se3, γ-In2Te3, and γ-In2(Se1-xTex)3 alloys to the available toolbox of materials for solar cells and other optoelectronic applications.
Collapse
Affiliation(s)
- Wei Li
- Department of Materials Science & Engineering, University of Delaware, Newark, Delaware19716, United States
- Computer, Computational, and Statistical Sciences Division, Los Alamos National Laboratory, Los Alamos, New Mexico87545, United States
| | - Xue-Fen Cai
- Department of Materials Science & Engineering, University of Delaware, Newark, Delaware19716, United States
| | - Nicholas Valdes
- Department of Materials Science & Engineering and Institute of Energy Conversion, University of Delaware, Newark, Delaware19716, United States
| | - Tianshi Wang
- Department of Materials Science & Engineering and Institute of Energy Conversion, University of Delaware, Newark, Delaware19716, United States
| | - William Shafarman
- Department of Materials Science & Engineering and Institute of Energy Conversion, University of Delaware, Newark, Delaware19716, United States
| | - Su-Huai Wei
- Beijing Computational Science Research Center, Beijing100193, China
| | - Anderson Janotti
- Department of Materials Science & Engineering, University of Delaware, Newark, Delaware19716, United States
| |
Collapse
|
22
|
Förster A. Assessment of the Second-Order Statically Screened Exchange Correction to the Random Phase Approximation for Correlation Energies. J Chem Theory Comput 2022; 18:5948-5965. [PMID: 36150190 PMCID: PMC9558381 DOI: 10.1021/acs.jctc.2c00366] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
![]()
With increasing interelectronic distance, the screening
of the
electron–electron interaction by the presence of other electrons
becomes the dominant source of electron correlation. This effect is
described by the random phase approximation (RPA) which is therefore
a promising method for the calculation of weak interactions. The success
of the RPA relies on the cancellation of errors, which can be traced
back to the violation of the crossing symmetry of the 4-point vertex,
leading to strongly overestimated total correlation energies. By the
addition of second-order screened exchange (SOSEX) to the correlation
energy, this issue is substantially reduced. In the adiabatic connection
(AC) SOSEX formalism, one of the two electron–electron interaction
lines in the second-order exchange term is dynamically screened (SOSEX(W, vc)). A
related SOSEX expression in which both electron–electron interaction
lines are statically screened (SOSEX(W(0), W(0))) is obtained from the G3W2 contribution to the electronic self-energy. In contrast to SOSEX(W, vc), the
evaluation of this correlation energy expression does not require
an expensive numerical frequency integration and is therefore advantageous
from a computational perspective. We compare the accuracy of the statically
screened variant to RPA and RPA+SOSEX(W, vc) for a wide range of chemical
reactions. While both methods fail for barrier heights, SOSEX(W(0), W(0)) agrees very well with SOSEX(W, vc) for
charged excitations and noncovalent interactions where they lead to
major improvements over RPA.
Collapse
Affiliation(s)
- Arno Förster
- Theoretical Chemistry, Vrije Universiteit, De Boelelaan 1083, NL-1081 HV, Amsterdam, The Netherlands
| |
Collapse
|
23
|
Ghosh A, Jana S, Rauch T, Tran F, Marques MAL, Botti S, Constantin L, Niranjan MK, Samal P. Efficient and improved prediction of the band offsets at semiconductorheterojunctions from meta-GGA density functionals: a benchmark study. J Chem Phys 2022; 157:124108. [DOI: 10.1063/5.0111693] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Abstract
Accurate theoretical prediction of the band offsets at interfaces of semiconductor heterostructures can of-ten be quite challenging. Although density functional theory has been reasonably successful to carry outsuch calculations and efficient and accurate semilocal functionals are desirable to reduce the computational cost. In general, the semilocal functionals based on the generalized gradient approximation (GGA) significantly underestimate the bulk band gaps. This, in turn, results in inaccurate estimates of the band offsets at the heterointerfaces. In this paper, we investigate the performance of several advanced meta-GGA functionals in the computational prediction of band offsets at semiconductor heterojunctions. In particular, we investigate the performance of r 2 SCAN (revised strongly-constrained and appropriately-normed functional), rMGGAC (revised semilocal functional based on cuspless hydrogen model and Pauli kinetic energy density functional), mTASK (modified Aschebrock and Kümmel meta-GGA functional), and LMBJ (local modified Becke-Johnson) exchange-correlation functionals. Our results strongly suggest that these meta-GGA functionals for supercell calculations perform quite well, especially, when compared to computationally more demanding GW calculations. We also present band offsets calculated using ionization potentials and electron affinities, as well as band alignment via the branch point energies. Overall, our study shows that the aforementioned meta-GGA functionals can be used within the DFT framework to estimate the band offsets in semiconductor heterostructures with predictive accuracy.
Collapse
Affiliation(s)
| | - Subrata Jana
- Department of Chemistry and Biochemistry, The Ohio State University, United States of America
| | - Tomas Rauch
- Friedrich Schiller Universität Jena Institut für Festkörpertheorie und -optik, Germany
| | - Fabien Tran
- Institute of Materials Chemistry, Vienna University of Technology, Austria
| | | | - Silvana Botti
- Institut für Festkörpertheorie und -optik, Friedrich Schiller Universität Jena Institut für Festkörpertheorie und -optik, Germany
| | - Lucian Constantin
- Department of Physics, Istituto di Nanoscienze, Consiglio Nazionale delle Ricerche CNR-NANO, 41125 Modena, Italy, Italy
| | | | - Prasanjit Samal
- School of Physical Sciences, National Institute of Science Education and Research, India
| |
Collapse
|
24
|
Vacondio S, Varsano D, Ruini A, Ferretti A. Numerically Precise Benchmark of Many-Body Self-Energies on Spherical Atoms. J Chem Theory Comput 2022; 18:3703-3717. [PMID: 35561415 PMCID: PMC9202310 DOI: 10.1021/acs.jctc.2c00048] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
![]()
We investigate the
performance of beyond-GW approaches in many-body
perturbation theory by addressing atoms described within the spherical
approximation via a dedicated numerical treatment based on B-splines
and spherical harmonics. We consider the GW, second Born (2B), and
GW + second order screened exchange (GW+SOSEX) self-energies and use
them to obtain ionization potentials from the quasi-particle equation
(QPE) solved perturbatively on top of independent-particle calculations.
We also solve the linearized Sham–Schlüter equation
(LSSE) and compare the resulting xc potentials against exact data.
We find that the LSSE provides consistent starting points for the
QPE but does not present any practical advantage in the present context.
Still, the features of the xc potentials obtained with it shed light
on possible strategies for the inclusion of beyond-GW diagrams in
the many-body self-energy. Our findings show that solving the QPE
with the GW+SOSEX self-energy on top of a PBE or PBE0 solution is
a viable scheme to go beyond GW in finite systems, even in the atomic
limit. However, GW shows a comparable performance if one agrees to
use a hybrid starting point. We also obtain promising results with
the 2B self-energy on top of Hartree–Fock, suggesting that
the full time-dependent Hartree–Fock vertex may be another
viable beyond-GW scheme for finite systems.
Collapse
Affiliation(s)
- S Vacondio
- Dipartimento di Scienze Fisiche, Informatiche e Matematiche, Università di Modena e Reggio Emilia, Via G. Campi 213/a, Modena 41121, Italy.,Centro S3, CNR-Istituto Nanoscienze, 41125 Modena, Italy
| | - D Varsano
- Centro S3, CNR-Istituto Nanoscienze, 41125 Modena, Italy
| | - A Ruini
- Dipartimento di Scienze Fisiche, Informatiche e Matematiche, Università di Modena e Reggio Emilia, Via G. Campi 213/a, Modena 41121, Italy.,Centro S3, CNR-Istituto Nanoscienze, 41125 Modena, Italy
| | - A Ferretti
- Centro S3, CNR-Istituto Nanoscienze, 41125 Modena, Italy
| |
Collapse
|
25
|
Santra G, Semidalas E, Mehta N, Karton A, Martin JML. S66x8 noncovalent interactions revisited: new benchmark and performance of composite localized coupled-cluster methods. Phys Chem Chem Phys 2022; 24:25555-25570. [DOI: 10.1039/d2cp03938a] [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/14/2022]
Abstract
The S66x8 noncovalent interactions benchmark has been re-evaluated at the “sterling silver” level. Against this, a selection of computationally more economical alternatives has been assayed, ranging from localized CC to double hybrids and SAPT(DFT).
Collapse
Affiliation(s)
- Golokesh Santra
- Department of Molecular Chemistry and Materials Science, Weizmann Institute of Science, 7610001 Reḥovot, Israel
| | - Emmanouil Semidalas
- Department of Molecular Chemistry and Materials Science, Weizmann Institute of Science, 7610001 Reḥovot, Israel
| | - Nisha Mehta
- Department of Molecular Chemistry and Materials Science, Weizmann Institute of Science, 7610001 Reḥovot, Israel
| | - Amir Karton
- School of Molecular Sciences, The University of Western Australia, Perth, WA 6009, Australia
- School of Science and Technology, University of New England, Armidale, NSW 2351, Australia
| | - Jan M. L. Martin
- Department of Molecular Chemistry and Materials Science, Weizmann Institute of Science, 7610001 Reḥovot, Israel
| |
Collapse
|
26
|
Rauch T, Marques MAL, Botti S. Electronic Structure of Molecules, Surfaces, and Molecules on Surfaces with the Local Modified Becke-Johnson Exchange-Correlation Potential. J Chem Theory Comput 2021; 17:4746-4755. [PMID: 34242509 DOI: 10.1021/acs.jctc.1c00255] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
The knowledge of electronic properties of matter is the key to the understanding of its properties and to propose useful applications. To model hybrid organic/inorganic systems with the plane-wave approach, large supercells with many atoms are usually necessary to minimize artificial interactions between periodic images. For such systems, accurate approximations to the exchange-correlation functional of density functional theory, such as hybrid functionals, become computationally expensive, and cheaper approaches need to be considered. Here, we apply the local modified Becke-Johnson exchange-correlation potential to free molecules and surfaces and study its accuracy for calculated ionization potentials. This quantity being important to understand the band alignment of composite heterogeneous systems, we demonstrate the application of the potential to the electronic structure calculation of an exemplary composite semiconductor/molecule system, namely, a F6-TCNNQ molecule adsorbed on a hydrogenated Si(111) surface.
Collapse
Affiliation(s)
- Tomáš Rauch
- Institut für Festkörpertheorie und -optik, Friedrich-Schiller-Universität Jena, Max-Wien-Platz 1, 07743 Jena, Germany
| | - Miguel A L Marques
- Institut für Physik, Martin-Luther-Universität Halle-Wittenberg, 06120 Halle/Saale, Germany.,European Theoretical Spectroscopy Facility, https://www.etsf.eu
| | - Silvana Botti
- Institut für Festkörpertheorie und -optik, Friedrich-Schiller-Universität Jena, Max-Wien-Platz 1, 07743 Jena, Germany.,European Theoretical Spectroscopy Facility, https://www.etsf.eu
| |
Collapse
|
27
|
Sun HY, Li SX, Jiang H. Pros and cons of the time-dependent hybrid density functional approach for calculating the optical spectra of solids: a case study of CeO 2. Phys Chem Chem Phys 2021; 23:16296-16306. [PMID: 34312647 DOI: 10.1039/d1cp02049h] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/20/2022]
Abstract
The prediction of optical spectra of complex solids remains a great challenge for first-principles calculations due to the huge computational cost of the state-of-the-art many-body perturbation theory based GW-Bethe Salpeter equation (BSE) approach. An alternative method is the time-dependent density-functional theory (TDDFT) based on hybrid exchange-correlation functionals, which involves the essential ingredients of electron-hole interactions in its formalism in contrast to its local/semi-local functional counterparts. In this work, we investigate the optical absorption spectra of ceria (CeO2), a prototypical lanthanide oxide with a 4f0 configuration, utilizing TDDFT based on four well-established hybrid functionals for ground state DFT calculations. All four functionals reproduce well the excitonic features of the experimental optical spectra, in spite of the significant differences in their band structures arising from different hybridization parameters (i.e. the fraction of the Hartree-Fock exchange and the screening parameter). It is demonstrated that the apparently weak dependence of the resulting optical spectra on the employed functionals is quite universal and applies to simple semiconductors such as Si and GaAs and insulator LiF as well. This study highlights the feasibility of TDDFT based on existing hybrids to describe optical spectra of solids, and also, points out the difficulty of obtaining accurate exciton binding energies using these hybrid functionals due to the strong functional dependence of quasi-particle band structures.
Collapse
Affiliation(s)
- Huai-Yang Sun
- Beijing National Laboratory for Molecular Sciences, Institute of Theoretical and Computational Chemistry, College of Chemistry and Molecular Engineering, Peking University, Beijing 100871, China.
| | | | | |
Collapse
|
28
|
Wang Y, Rinke P, Ren X. Assessing the G0W0Γ 0(1) Approach: Beyond G0W0 with Hedin's Full Second-Order Self-Energy Contribution. J Chem Theory Comput 2021; 17:5140-5154. [PMID: 34319724 DOI: 10.1021/acs.jctc.1c00488] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
We present and benchmark a self-energy approach for quasiparticle energy calculations that goes beyond Hedin's GW approximation by adding the full second-order self-energy (FSOS-W) contribution. The FSOS-W diagram involves two screened Coulomb interaction (W) lines, and adding the FSOS-W to the GW self-energy can be interpreted as first-order vertex correction to GW (GWΓ(1)). Our FSOS-W implementation is based on the resolution-of-identity technique and exhibits better than O(N5) scaling with system size for small- to medium-sized molecules. We then present one-shot GWΓ(1) (G0W0Γ0(1)) benchmarks for the GW100 test set and a set of 24 acceptor molecules. For semilocal or hybrid density functional theory starting points, G0W0Γ0(1) systematically outperforms G0W0 for the first vertical ionization potentials and electron affinities of both test sets. Finally, we demonstrate that a static FSOS-W self-energy significantly underestimates the quasiparticle energies.
Collapse
Affiliation(s)
- Yanyong Wang
- CAS Key Laboratory of Quantum Information, University of Science and Technology of China, Hefei, Anhui 230026, China
| | - Patrick Rinke
- Department of Applied Physics, School of Science, Aalto University, 00076 Aalto, Finland
| | - Xinguo Ren
- Beijing National Laboratory for Condensed Matter Physics, Institute of Physics, Chinese Academy of Sciences, Beijing 100190, China
| |
Collapse
|
29
|
Riemelmoser S, Kaltak M, Kresse G. Optimized effective potentials from the random-phase approximation: Accuracy of the quasiparticle approximation. J Chem Phys 2021; 154:154103. [PMID: 33887939 DOI: 10.1063/5.0045400] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/13/2023] Open
Abstract
The optimized effective potential (OEP) method presents an unambiguous way to construct the Kohn-Sham potential corresponding to a given diagrammatic approximation for the exchange-correlation functional. The OEP from the random-phase approximation (RPA) has played an important role ever since the conception of the OEP formalism. However, the solution of the OEP equation is computationally fairly expensive and has to be done in a self-consistent way. So far, large scale solid state applications have, therefore, been performed only using the quasiparticle approximation (QPA), neglecting certain dynamical screening effects. We obtain the exact RPA-OEP for 15 semiconductors and insulators by direct solution of the linearized Sham-Schlüter equation. We investigate the accuracy of the QPA on Kohn-Sham bandgaps and dielectric constants, and comment on the issue of self-consistency.
Collapse
Affiliation(s)
- Stefan Riemelmoser
- Faculty of Physics and Center for Computational Materials Science, University of Vienna, Kolingasse 14-16, A-1090 Vienna, Austria
| | - Merzuk Kaltak
- VASP Software GmbH, Sensengasse 8/17, A-1090 Vienna, Austria
| | - Georg Kresse
- Faculty of Physics and Center for Computational Materials Science, University of Vienna, Kolingasse 14-16, A-1090 Vienna, Austria
| |
Collapse
|
30
|
Mejuto-Zaera C, Weng G, Romanova M, Cotton SJ, Whaley KB, Tubman NM, Vlček V. Are multi-quasiparticle interactions important in molecular ionization? J Chem Phys 2021; 154:121101. [DOI: 10.1063/5.0044060] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/06/2023] Open
Affiliation(s)
| | - Guorong Weng
- University of California, Santa Barbara, California 93106, USA
| | - Mariya Romanova
- University of California, Santa Barbara, California 93106, USA
| | - Stephen J. Cotton
- Quantum Artificial Intelligence Laboratory (QuAIL), Exploration Technology Directorate, NASA Ames Research Center, Moffett Field, California 94035, USA
- KBR, 601 Jefferson St., Houston, Texas 77002, USA
| | | | - Norm M. Tubman
- Quantum Artificial Intelligence Laboratory (QuAIL), Exploration Technology Directorate, NASA Ames Research Center, Moffett Field, California 94035, USA
| | - Vojtěch Vlček
- University of California, Santa Barbara, California 93106, USA
| |
Collapse
|
31
|
Ma H, Sheng N, Govoni M, Galli G. Quantum Embedding Theory for Strongly Correlated States in Materials. J Chem Theory Comput 2021; 17:2116-2125. [DOI: 10.1021/acs.jctc.0c01258] [Citation(s) in RCA: 20] [Impact Index Per Article: 6.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/19/2022]
Affiliation(s)
- He Ma
- Department of Chemistry, University of Chicago, Chicago, Illinois 60637, United States
| | - Nan Sheng
- Department of Chemistry, University of Chicago, Chicago, Illinois 60637, United States
| | - Marco Govoni
- Pritzker School of Molecular Engineering, University of Chicago, Chicago, Illinois 60637, United States
- Materials Science Division and Center for Molecular Engineering, Argonne National Laboratory, Lemont, Illinois 60439, United States
| | - Giulia Galli
- Department of Chemistry, University of Chicago, Chicago, Illinois 60637, United States
- Pritzker School of Molecular Engineering, University of Chicago, Chicago, Illinois 60637, United States
- Materials Science Division and Center for Molecular Engineering, Argonne National Laboratory, Lemont, Illinois 60439, United States
| |
Collapse
|
32
|
Kim K, Siegel DJ. Predicting Wettability and the Electrochemical Window of Lithium-Metal/Solid Electrolyte Interfaces. ACS APPLIED MATERIALS & INTERFACES 2019; 11:39940-39950. [PMID: 31576739 DOI: 10.1021/acsami.9b13311] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/10/2023]
Abstract
The development of solid electrolytes (SEs) is expected to enhance the safety of lithium-ion batteries. Additionally, a viable SE could allow the use of a Li-metal negative electrode, which would increase energy density. Recently, several antiperovskites have been reported to exhibit high ionic conductivities, prompting investigations of their use as an SE. In addition to having a suitable conductivity, phenomena at the interface between an SE and an electrode are also of great importance in determining the viability of an SE. For example, interfacial interactions can change the positions of the band edges of the SE, altering its stability against undesirable oxidation or reduction. Furthermore, the wettability of the SE by the metallic anode is desired to enable low interfacial resistance and uniform metal plating and stripping during cycling. The present study probes several properties of the SE/electrode interface at the atomic scale. Adopting the antiperovskite SE Li3OCl (LOC)/Li-metal anode interface as a model system, the interfacial energy, work of adhesion, wettability, band edge shifts, and the electrochemical window are predicted computationally. The oxygen-terminated interface was determined to be the most thermodynamically stable. Moreover, the large calculated work of adhesion for this system implies that Li will wet LOC, suggesting the possibility for low interfacial resistance. Nevertheless, these strong interfacial interactions come at a cost to electrochemical stability: strong interfacial bonding lowers the energy of the conduction band minimum (CBM) significantly and narrows the local band gap by 30% in the vicinity of the interface. Despite this interface-induced reduction in electrochemical window, the CBM in LOC remains more negative than the Li/Li+ redox potential, implying stability against reduction by the anode. In sum, this study illustrates a comprehensive computational approach to assessing electrode/electrolyte interfacial properties in solid-state batteries.
Collapse
|
33
|
Golze D, Dvorak M, Rinke P. The GW Compendium: A Practical Guide to Theoretical Photoemission Spectroscopy. Front Chem 2019; 7:377. [PMID: 31355177 PMCID: PMC6633269 DOI: 10.3389/fchem.2019.00377] [Citation(s) in RCA: 153] [Impact Index Per Article: 30.6] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/14/2018] [Accepted: 05/08/2019] [Indexed: 12/22/2022] Open
Abstract
The GW approximation in electronic structure theory has become a widespread tool for predicting electronic excitations in chemical compounds and materials. In the realm of theoretical spectroscopy, the GW method provides access to charged excitations as measured in direct or inverse photoemission spectroscopy. The number of GW calculations in the past two decades has exploded with increased computing power and modern codes. The success of GW can be attributed to many factors: favorable scaling with respect to system size, a formal interpretation for charged excitation energies, the importance of dynamical screening in real systems, and its practical combination with other theories. In this review, we provide an overview of these formal and practical considerations. We expand, in detail, on the choices presented to the scientist performing GW calculations for the first time. We also give an introduction to the many-body theory behind GW, a review of modern applications like molecules and surfaces, and a perspective on methods which go beyond conventional GW calculations. This review addresses chemists, physicists and material scientists with an interest in theoretical spectroscopy. It is intended for newcomers to GW calculations but can also serve as an alternative perspective for experts and an up-to-date source of computational techniques.
Collapse
Affiliation(s)
- Dorothea Golze
- Department of Applied Physics, Aalto University, School of Science, Espoo, Finland
| | | | | |
Collapse
|
34
|
Lewis AM, Berkelbach TC. Ab Initio Lifetime and Concomitant Double-Excitation Character of Plasmons at Metallic Densities. PHYSICAL REVIEW LETTERS 2019; 122:226402. [PMID: 31283277 DOI: 10.1103/physrevlett.122.226402] [Citation(s) in RCA: 12] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 10/15/2018] [Indexed: 06/09/2023]
Abstract
The accurate calculation of excited state properties of interacting electrons in the condensed phase is an immense challenge in computational physics. Here, we use state-of-the-art equation-of-motion coupled-cluster theory with single and double excitations (EOM-CCSD) to calculate the dynamic structure factor, which can be experimentally measured by inelastic x-ray and electron scattering. Our calculations are performed on the uniform electron gas at densities corresponding to Wigner-Seitz radii of r_{s}=5, 4, and 3 corresponding to the valence electron densities of common metals. We compare our results to those obtained using the random-phase approximation (RPA), which is known to provide a reasonable description of the collective plasmon excitation and which resums only a small subset of the polarizability diagrams included in EOM-CCSD. We find that EOM-CCSD, instead of providing a perturbative improvement on the RPA plasmon, predicts a many-state plasmon resonance, where each contributing state has a double-excitation character of 80% or more. This finding amounts to an ab initio treatment of the plasmon linewidth, which is in good quantitative agreement with previous diagrammatic calculations, and highlights the strongly correlated nature of lifetime effects in condensed-phase electronic structure theory.
Collapse
Affiliation(s)
- Alan M Lewis
- Department of Chemistry and James Franck Institute, University of Chicago, Chicago, Illinois 60637, USA
| | - Timothy C Berkelbach
- Department of Chemistry, Columbia University, New York, New York 10027 USA
- Center for Computational Quantum Physics, Flatiron Institute, New York, New York 10010, USA
| |
Collapse
|
35
|
Li LF, Li YF, Liu ZP. CO2 Photoreduction via Quantum Tunneling: Thin TiO2-Coated GaP with Coherent Interface To Achieve Electron Tunneling. ACS Catal 2019. [DOI: 10.1021/acscatal.9b01645] [Citation(s) in RCA: 14] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Li-Fen Li
- Collaborative Innovation Center of Chemistry for Energy Material, Key Laboratory of Computational Physical Science (Ministry of Education), Shanghai Key Laboratory of Molecular Catalysis and Innovative Materials, Department of Chemistry, Fudan University, Shanghai 200433, China
| | - Ye-Fei Li
- Collaborative Innovation Center of Chemistry for Energy Material, Key Laboratory of Computational Physical Science (Ministry of Education), Shanghai Key Laboratory of Molecular Catalysis and Innovative Materials, Department of Chemistry, Fudan University, Shanghai 200433, China
| | - Zhi-Pan Liu
- Collaborative Innovation Center of Chemistry for Energy Material, Key Laboratory of Computational Physical Science (Ministry of Education), Shanghai Key Laboratory of Molecular Catalysis and Innovative Materials, Department of Chemistry, Fudan University, Shanghai 200433, China
| |
Collapse
|
36
|
Lewis AM, Berkelbach TC. Vertex Corrections to the Polarizability Do Not Improve the GW Approximation for the Ionization Potential of Molecules. J Chem Theory Comput 2019; 15:2925-2932. [DOI: 10.1021/acs.jctc.8b00995] [Citation(s) in RCA: 12] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Alan M. Lewis
- Department of Chemistry and James Franck Institute, University of Chicago, Chicago, Illinois 60637, United States
| | - Timothy C. Berkelbach
- Department of Chemistry, Columbia University, New York, New York 10027, United States
- Center for Computational Quantum Physics, Flatiron Institute, New York, New York 10010, United States
| |
Collapse
|
37
|
Ma H, Govoni M, Gygi F, Galli G. A Finite-Field Approach for GW Calculations beyond the Random Phase Approximation. J Chem Theory Comput 2018; 15:154-164. [DOI: 10.1021/acs.jctc.8b00864] [Citation(s) in RCA: 19] [Impact Index Per Article: 3.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- He Ma
- Institute for Molecular Engineering, University of Chicago, Chicago, Illinois 60637, United States
- Department of Chemistry, University of Chicago, Chicago, Illinois 60637, United States
| | - Marco Govoni
- Institute for Molecular Engineering, University of Chicago, Chicago, Illinois 60637, United States
- Materials Science Division, Argonne National Laboratory, Lemont, Illinois 60439, United States
| | - Francois Gygi
- Department of Computer Science, University of California Davis, Davis, California 95616, United States
| | - Giulia Galli
- Institute for Molecular Engineering, University of Chicago, Chicago, Illinois 60637, United States
- Department of Chemistry, University of Chicago, Chicago, Illinois 60637, United States
- Materials Science Division, Argonne National Laboratory, Lemont, Illinois 60439, United States
| |
Collapse
|
38
|
Vlček V, Baer R, Rabani E, Neuhauser D. Simple eigenvalue-self-consistent Δ ¯ G W 0 . J Chem Phys 2018; 149:174107. [PMID: 30409020 DOI: 10.1063/1.5042785] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
We show that a rigid scissors-like GW self-consistency approach, labeled here Δ ¯ G W 0 , can be trivially implemented at zero additional cost for large scale one-shot G 0 W 0 calculations. The method significantly improves one-shot G 0 W 0 and for large systems is very accurate. Δ ¯ G W 0 is similar in spirit to evGW 0 where the self-consistency is only applied on the eigenvalues entering Green's function, while both W and the eigenvectors of Green's function are held fixed. Δ ¯ G W 0 further assumes that the shift of the eigenvalues is rigid scissors-like so that all occupied states are shifted by the same amount and analogously for all the unoccupied states. We show that this results in a trivial modification of the time-dependent G 0 W 0 self-energy, enabling an a posteriori self-consistency cycle. The method is applicable for our recent stochastic-GW approach, thereby enabling self-consistent calculations for giant systems with thousands of electrons. The accuracy of Δ ¯ G W 0 increases with the system size. For molecules, it is up to 0.4-0.5 eV away from coupled-cluster single double triple (CCSD(T)), but for tetracene and hexacene, it matches the ionization energies from both CCSD(T) and evGW 0 to better than 0.05 eV. For solids, as exemplified here by periodic supercells of semiconductors and insulators with 6192 valence electrons, the method matches evGW 0 quite well and both methods are in good agreement with the experiment.
Collapse
Affiliation(s)
- Vojtěch Vlček
- Department of Chemistry and Biochemistry, University of California, Los Angeles, California 90095, USA
| | - Roi Baer
- Fritz Haber Center for Molecular Dynamics, Institute of Chemistry, The Hebrew University of Jerusalem, Jerusalem 91904, Israel
| | - Eran Rabani
- Department of Chemistry, University of California and Materials Science Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, USA
| | - Daniel Neuhauser
- Department of Chemistry and Biochemistry, University of California, Los Angeles, California 90095, USA
| |
Collapse
|
39
|
Ghosh S, Verma P, Cramer CJ, Gagliardi L, Truhlar DG. Combining Wave Function Methods with Density Functional Theory for Excited States. Chem Rev 2018; 118:7249-7292. [PMID: 30044618 DOI: 10.1021/acs.chemrev.8b00193] [Citation(s) in RCA: 128] [Impact Index Per Article: 21.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/29/2022]
Abstract
We review state-of-the-art electronic structure methods based both on wave function theory (WFT) and density functional theory (DFT). Strengths and limitations of both the wave function and density functional based approaches are discussed, and modern attempts to combine these two methods are presented. The challenges in modeling excited-state chemistry using both single-reference and multireference methods are described. Topics covered include background, combining density functional theory with single-configuration wave function theory, generalized Kohn-Sham (KS) theory, global hybrids, range-separated hybrids, local hybrids, using KS orbitals in many-body theory (including calculations of the self-energy and the GW approximation), Bethe-Salpeter equation, algorithms to accelerate GW calculations, combining DFT with multiconfigurational WFT, orbital-dependent correlation functionals based on multiconfigurational WFT, building multiconfigurational wave functions from KS configurations, adding correlation functionals to multiconfiguration self-consistent-field (MCSCF) energies, combining DFT with configuration-interaction singles by means of time-dependent DFT, using range separation to combine DFT with MCSCF, embedding multiconfigurational WFT in DFT, and multiconfiguration pair-density functional theory.
Collapse
Affiliation(s)
- Soumen Ghosh
- Department of Chemistry, Chemical Theory Center, and Minnesota Supercomputing Institute , University of Minnesota , 207 Pleasant Street SE , Minneapolis , Minnesota 55455-0431 , United States
| | - Pragya Verma
- Department of Chemistry, Chemical Theory Center, and Minnesota Supercomputing Institute , University of Minnesota , 207 Pleasant Street SE , Minneapolis , Minnesota 55455-0431 , United States
| | - Christopher J Cramer
- Department of Chemistry, Chemical Theory Center, and Minnesota Supercomputing Institute , University of Minnesota , 207 Pleasant Street SE , Minneapolis , Minnesota 55455-0431 , United States
| | - Laura Gagliardi
- Department of Chemistry, Chemical Theory Center, and Minnesota Supercomputing Institute , University of Minnesota , 207 Pleasant Street SE , Minneapolis , Minnesota 55455-0431 , United States
| | - Donald G Truhlar
- Department of Chemistry, Chemical Theory Center, and Minnesota Supercomputing Institute , University of Minnesota , 207 Pleasant Street SE , Minneapolis , Minnesota 55455-0431 , United States
| |
Collapse
|
40
|
Lange MF, Berkelbach TC. On the Relation between Equation-of-Motion Coupled-Cluster Theory and the GW Approximation. J Chem Theory Comput 2018; 14:4224-4236. [DOI: 10.1021/acs.jctc.8b00455] [Citation(s) in RCA: 45] [Impact Index Per Article: 7.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Malte F. Lange
- Department of Chemistry and James Franck Institute, University of Chicago, Chicago, Illinois 60637, United States
| | - Timothy C. Berkelbach
- Department of Chemistry and James Franck Institute, University of Chicago, Chicago, Illinois 60637, United States
| |
Collapse
|
41
|
Wang X, Liu X, Cook C, Schatschneider B, Marom N. On the possibility of singlet fission in crystalline quaterrylene. J Chem Phys 2018; 148:184101. [DOI: 10.1063/1.5027553] [Citation(s) in RCA: 14] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/09/2023] Open
Affiliation(s)
- Xiaopeng Wang
- Department of Materials Science and Engineering, Carnegie Mellon University, Pittsburgh, Pennsylvania 15213, USA
| | - Xingyu Liu
- Department of Materials Science and Engineering, Carnegie Mellon University, Pittsburgh, Pennsylvania 15213, USA
| | - Cameron Cook
- Department of Chemistry and Biochemistry, California State Polytechnic University at Pomona, Pomona, California 91768, USA
| | - Bohdan Schatschneider
- Department of Chemistry and Biochemistry, California State Polytechnic University at Pomona, Pomona, California 91768, USA
| | - Noa Marom
- Department of Materials Science and Engineering, Carnegie Mellon University, Pittsburgh, Pennsylvania 15213, USA
- Department of Chemistry, Carnegie Mellon University, Pittsburgh, Pennsylvania 15213, USA
- Department of Physics, Carnegie Mellon University, Pittsburgh, Pennsylvania 15213, USA
| |
Collapse
|
42
|
Gerosa M, Bottani CE, Di Valentin C, Onida G, Pacchioni G. Accuracy of dielectric-dependent hybrid functionals in the prediction of optoelectronic properties of metal oxide semiconductors: a comprehensive comparison with many-body GW and experiments. JOURNAL OF PHYSICS. CONDENSED MATTER : AN INSTITUTE OF PHYSICS JOURNAL 2018; 30:044003. [PMID: 29087359 DOI: 10.1088/1361-648x/aa9725] [Citation(s) in RCA: 20] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/07/2023]
Abstract
Understanding the electronic structure of metal oxide semiconductors is crucial to their numerous technological applications, such as photoelectrochemical water splitting and solar cells. The needed experimental and theoretical knowledge goes beyond that of pristine bulk crystals, and must include the effects of surfaces and interfaces, as well as those due to the presence of intrinsic defects (e.g. oxygen vacancies), or dopants for band engineering. In this review, we present an account of the recent efforts in predicting and understanding the optoelectronic properties of oxides using ab initio theoretical methods. In particular, we discuss the performance of recently developed dielectric-dependent hybrid functionals, providing a comparison against the results of many-body GW calculations, including G 0 W 0 as well as more refined approaches, such as quasiparticle self-consistent GW. We summarize results in the recent literature for the band gap, the band level alignment at surfaces, and optical transition energies in defective oxides, including wide gap oxide semiconductors and transition metal oxides. Correlated transition metal oxides are also discussed. For each method, we describe successes and drawbacks, emphasizing the challenges faced by the development of improved theoretical approaches. The theoretical section is preceded by a critical overview of the main experimental techniques needed to characterize the optoelectronic properties of semiconductors, including absorption and reflection spectroscopy, photoemission, and scanning tunneling spectroscopy (STS).
Collapse
Affiliation(s)
- M Gerosa
- Institute for Molecular Engineering, University of Chicago, Chicago, IL 60637, United States of America
| | | | | | | | | |
Collapse
|
43
|
Maggio E, Kresse G. GW Vertex Corrected Calculations for Molecular Systems. J Chem Theory Comput 2017; 13:4765-4778. [DOI: 10.1021/acs.jctc.7b00586] [Citation(s) in RCA: 36] [Impact Index Per Article: 5.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Emanuele Maggio
- Faculty of Physics
and Center
for Computational Materials Science, University of Vienna, Sensengasse
8/12, A-1090 Vienna, Austria
| | - Georg Kresse
- Faculty of Physics
and Center
for Computational Materials Science, University of Vienna, Sensengasse
8/12, A-1090 Vienna, Austria
| |
Collapse
|
44
|
Hung L, Bruneval F, Baishya K, Öğüt S. Benchmarking the GW Approximation and Bethe–Salpeter Equation for Groups IB and IIB Atoms and Monoxides. J Chem Theory Comput 2017; 13:2135-2146. [DOI: 10.1021/acs.jctc.7b00123] [Citation(s) in RCA: 26] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/19/2023]
Affiliation(s)
- Linda Hung
- NIST
Center for Neutron Research, National Institute of Standard and Technology, Gaithersburg, Maryland 20899, United States
| | - Fabien Bruneval
- CEA,
DEN, Service de Recherches de Métallurgie Physique, Université Paris-Saclay, F-91191 Gif-sur-Yvette, France
| | - Kopinjol Baishya
- Department
of Physics, University of Illinois at Chicago, Chicago, Illinois 60607, United States
| | - Serdar Öğüt
- Department
of Physics, University of Illinois at Chicago, Chicago, Illinois 60607, United States
| |
Collapse
|
45
|
Marom N. Accurate description of the electronic structure of organic semiconductors by GW methods. JOURNAL OF PHYSICS. CONDENSED MATTER : AN INSTITUTE OF PHYSICS JOURNAL 2017; 29:103003. [PMID: 28145283 DOI: 10.1088/1361-648x/29/10/103003] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/06/2023]
Abstract
Electronic properties associated with charged excitations, such as the ionization potential (IP), the electron affinity (EA), and the energy level alignment at interfaces, are critical parameters for the performance of organic electronic devices. To computationally design organic semiconductors and functional interfaces with tailored properties for target applications it is necessary to accurately predict these properties from first principles. Many-body perturbation theory is often used for this purpose within the GW approximation, where G is the one particle Green's function and W is the dynamically screened Coulomb interaction. Here, the formalism of GW methods at different levels of self-consistency is briefly introduced and some recent applications to organic semiconductors and interfaces are reviewed.
Collapse
Affiliation(s)
- Noa Marom
- Department of Materials Science and Engineering, Department of Chemistry, and Department of Physics, Carnegie Mellon University, Pittsburgh, PA 15213, United States of America
| |
Collapse
|
46
|
McClain J, Sun Q, Chan GKL, Berkelbach TC. Gaussian-Based Coupled-Cluster Theory for the Ground-State and Band Structure of Solids. J Chem Theory Comput 2017; 13:1209-1218. [DOI: 10.1021/acs.jctc.7b00049] [Citation(s) in RCA: 130] [Impact Index Per Article: 18.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- James McClain
- Department
of Chemistry, Princeton University, Princeton, New Jersey 08544, United States
| | - Qiming Sun
- Division
of Chemistry and Chemical Engineering, California Institute of Technology, Pasadena, California 91125, United States
| | - Garnet Kin-Lic Chan
- Division
of Chemistry and Chemical Engineering, California Institute of Technology, Pasadena, California 91125, United States
| | - Timothy C. Berkelbach
- Department
of Chemistry and James Franck Institute, University of Chicago, Chicago, Illinois 60637, United States
| |
Collapse
|
47
|
Viñes F, Illas F. Electronic structure of stoichiometric and reduced ZnO from periodic relativistic all electron hybrid density functional calculations using numeric atom-centered orbitals. J Comput Chem 2017; 38:523-529. [DOI: 10.1002/jcc.24705] [Citation(s) in RCA: 14] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/20/2016] [Revised: 11/16/2016] [Accepted: 11/17/2016] [Indexed: 01/28/2023]
Affiliation(s)
- Francesc Viñes
- Departament de Ciència dels Materials i Química Física & Institut de Química Teòrica i Computacional (IQTCUB); Universitat de Barcelona; c/Martí i Franquès 1 Barcelona 08028 Spain
| | - Francesc Illas
- Departament de Ciència dels Materials i Química Física & Institut de Química Teòrica i Computacional (IQTCUB); Universitat de Barcelona; c/Martí i Franquès 1 Barcelona 08028 Spain
| |
Collapse
|
48
|
Chen W, Ambrosio F, Miceli G, Pasquarello A. Ab initio Electronic Structure of Liquid Water. PHYSICAL REVIEW LETTERS 2016; 117:186401. [PMID: 27835004 DOI: 10.1103/physrevlett.117.186401] [Citation(s) in RCA: 45] [Impact Index Per Article: 5.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 12/15/2015] [Indexed: 05/26/2023]
Abstract
Self-consistent GW calculations with efficient vertex corrections are employed to determine the electronic structure of liquid water. Nuclear quantum effects are taken into account through ab initio path-integral molecular dynamics simulations. We reveal a sizable band-gap renormalization of up to 0.7 eV due to hydrogen-bond quantum fluctuations. Our calculations lead to a band gap of 8.9 eV, in accord with the experimental estimate. We further resolve the ambiguities in the band-edge positions of liquid water. The valence-band maximum and the conduction-band minimum are found at -9.4 and -0.5 eV with respect to the vacuum level, respectively.
Collapse
Affiliation(s)
- Wei Chen
- Chaire de Simulation à l'Echelle Atomique (CSEA), Ecole Polytechnique Fédérale de Lausanne (EPFL), CH-1015 Lausanne, Switzerland
| | - Francesco Ambrosio
- Chaire de Simulation à l'Echelle Atomique (CSEA), Ecole Polytechnique Fédérale de Lausanne (EPFL), CH-1015 Lausanne, Switzerland
| | - Giacomo Miceli
- Chaire de Simulation à l'Echelle Atomique (CSEA), Ecole Polytechnique Fédérale de Lausanne (EPFL), CH-1015 Lausanne, Switzerland
| | - Alfredo Pasquarello
- Chaire de Simulation à l'Echelle Atomique (CSEA), Ecole Polytechnique Fédérale de Lausanne (EPFL), CH-1015 Lausanne, Switzerland
| |
Collapse
|
49
|
Lüdtke T, Weber D, Schmidt A, Müller A, Reimann C, Becker N, Bredow T, Dronskowski R, Ressler T, Lerch M. Synthesis and characterization of metastable transition metal oxides and oxide nitrides. Z KRIST-CRYST MATER 2016. [DOI: 10.1515/zkri-2016-1961] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022]
Abstract
Abstract
New routes to vanadium sesquioxide and tantalum oxide nitride (γ- and δ-phase) are presented. Phase pure V2O3 with bixbyite-type structure, a metastable polymorph, was obtained from vanadium fluoride hydrates at ~750 K. It crystallizes in the cubic crystal system in space group
I
a
3
¯
$Ia\bar 3$
with lattice parameter a=939.30(5) pm. The catalytical properties of the corresponding oxide nitride phases and their oxidation and reduction solid-state kinetics were investigated. The preparation of γ-TaON as a phase pure sample can be realized by ammonolysis of X-ray amorphous tantalum oxide precursors at 1073 K. This metastable tantalum oxide nitride crystallizes in the monoclinic VO2(B)-type structure in space group C2/m. The same precursors can be used to synthesize the δ-modification with an anatase-type structure at 1023 K. It crystallizes in the tetragonal crystal system in space group I41/amd. A maximum yield of 82 m % could be obtained. The fundamental band gaps of the synthesized and of other metastable TaON polymorphs were calculated from first principles using the GW method. The present results are compared to experimental data and to previous calculations at hybrid DFT level.
Collapse
Affiliation(s)
- Tobias Lüdtke
- Institut für Chemie, Technische Universität Berlin, Straße des 17. Juni 135, D-10623 Berlin, Germany
| | - Dominik Weber
- Institut für Chemie, Technische Universität Berlin, Straße des 17. Juni 135, D-10623 Berlin, Germany
| | - Alexander Schmidt
- Institut für Chemie, Technische Universität Berlin, Straße des 17. Juni 135, D-10623 Berlin, Germany
| | - Alexander Müller
- Institut für Chemie, Technische Universität Berlin, Straße des 17. Juni 135, D-10623 Berlin, Germany
| | - Christoph Reimann
- Mulliken Center for Theoretical Chemistry, Institut für Physikalische und Theoretische Chemie, University of Bonn, Beringstraße 4, D-53115 Bonn, Germany
| | - Nils Becker
- Chair of Solid-State and Quantum Chemistry, Institute of Inorganic Chemistry, RWTH Aachen University, Landoltweg 1, D-52056 Aachen, Germany
| | - Thomas Bredow
- Mulliken Center for Theoretical Chemistry, Institut für Physikalische und Theoretische Chemie, University of Bonn, Beringstraße 4, D-53115 Bonn, Germany
| | - Richard Dronskowski
- Chair of Solid-State and Quantum Chemistry, Institute of Inorganic Chemistry, RWTH Aachen University, Landoltweg 1, D-52056 Aachen, Germany
| | - Thorsten Ressler
- Institut für Chemie, Technische Universität Berlin, Straße des 17. Juni 135, D-10623 Berlin, Germany
| | - Martin Lerch
- Institut für Chemie, Technische Universität Berlin, Straße des 17. Juni 135, D-10623 Berlin, Germany
| |
Collapse
|
50
|
Scherpelz P, Govoni M, Hamada I, Galli G. Implementation and Validation of Fully Relativistic GW Calculations: Spin–Orbit Coupling in Molecules, Nanocrystals, and Solids. J Chem Theory Comput 2016; 12:3523-44. [DOI: 10.1021/acs.jctc.6b00114] [Citation(s) in RCA: 124] [Impact Index Per Article: 15.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Peter Scherpelz
- Institute
for Molecular Engineering, The University of Chicago, Chicago, Illinois 60637, United States
| | - Marco Govoni
- Institute
for Molecular Engineering, The University of Chicago, Chicago, Illinois 60637, United States
- Materials
Science Division, Argonne National Laboratory, Argonne, Illinois 60439, United States
| | - Ikutaro Hamada
- International
Center for Materials Nanoarchitectonics, Global Research Center for
Environment and Energy based on Nanomaterials Science, and Center
for Materials Research by Information Integration, National Institute for Materials Science, 1-1 Namiki, Tsukuba 305-0044, Japan
| | - Giulia Galli
- Institute
for Molecular Engineering, The University of Chicago, Chicago, Illinois 60637, United States
- Materials
Science Division, Argonne National Laboratory, Argonne, Illinois 60439, United States
| |
Collapse
|