1
|
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
|
2
|
Jin Y, Yu VWZ, Govoni M, Xu AC, Galli G. Excited State Properties of Point Defects in Semiconductors and Insulators Investigated with Time-Dependent Density Functional Theory. J Chem Theory Comput 2023. [PMID: 38039161 DOI: 10.1021/acs.jctc.3c00986] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/03/2023]
Abstract
We present a formulation of spin-conserving and spin-flip hybrid time-dependent density functional theory (TDDFT), including the calculation of analytical forces, which allows for efficient calculations of excited state properties of solid-state systems with hundreds to thousands of atoms. We discuss an implementation on both GPU- and CPU-based architectures along with several acceleration techniques. We then apply our formulation to the study of several point defects in semiconductors and insulators, specifically the negatively charged nitrogen-vacancy and neutral silicon-vacancy centers in diamond, the neutral divacancy center in 4H silicon carbide, and the neutral oxygen-vacancy center in magnesium oxide. Our results highlight the importance of taking into account structural relaxations in excited states in order to interpret and predict optical absorption and emission mechanisms in spin defects.
Collapse
Affiliation(s)
- Yu Jin
- Department of Chemistry, University of Chicago, Chicago, Illinois 60637, United States
| | - Victor Wen-Zhe Yu
- Materials Science Division, Argonne National Laboratory, Lemont, Illinois 60439, United States
| | - Marco Govoni
- Materials Science Division, Argonne National Laboratory, Lemont, Illinois 60439, United States
- Department of Physics, Computer Science, and Mathematics, University of Modena and Reggio Emilia, Modena 41125, Italy
- Pritzker School of Molecular Engineering, University of Chicago, Chicago, Illinois 60637, United States
| | - Andrew C Xu
- Pritzker School of Molecular Engineering, University of Chicago, Chicago, Illinois 60637, United States
| | - Giulia Galli
- Department of Chemistry, University of Chicago, Chicago, Illinois 60637, United States
- Materials Science Division, Argonne National Laboratory, Lemont, Illinois 60439, United States
- Pritzker School of Molecular Engineering, University of Chicago, Chicago, Illinois 60637, United States
| |
Collapse
|
3
|
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
|
4
|
Knysh I, Villalobos-Castro JDJ, Duchemin I, Blase X, Jacquemin D. Exploring Bethe-Salpeter Excited-State Dipoles: The Challenging Case of Increasingly Long Push-Pull Oligomers. J Phys Chem Lett 2023; 14:3727-3734. [PMID: 37042642 DOI: 10.1021/acs.jpclett.3c00699] [Citation(s) in RCA: 3] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/17/2023]
Abstract
The change of molecular dipole moment induced by photon absorption is key to interpret the measured optical spectra. Except for compact molecules, time-dependent density functional theory (TD-DFT) remains the only theory allowing to quickly predict excited-state dipoles (μES), albeit with a strong dependency on the selected exchange-correlation functional. This Letter presents the first assessment of the performances of the many-body Green's function Bethe-Salpeter equation (BSE) formalism for the evaluation of the μES. We explore increasingly long push-pull oligomers as they present an excited-state nature evolving with system size. This work shows that BSE's μES do present the same evolution with oligomeric length as their CC2 and CCSD counterparts, with a dependency on the starting exchange-correlation functional that is strongly decreased as compared to TD-DFT. This Letter demonstrates that BSE is a valuable alternative to TD-DFT for properties related to the excited-state density and not only for transition energies and oscillator strengths.
Collapse
Affiliation(s)
- Iryna Knysh
- Nantes Université, CNRS, CEISAM UMR 6230, F-44000 Nantes, France
| | | | - Ivan Duchemin
- Université Grenoble Alpes, CEA, IRIG-MEM-L Sim, 38054 Grenoble, France
| | - Xavier Blase
- Université Grenoble Alpes, CNRS, Institut NEEL, F-38042 Grenoble, France
| | - Denis Jacquemin
- Nantes Université, CNRS, CEISAM UMR 6230, F-44000 Nantes, France
- Institut Universitaire de France, 75005 Paris, France
| |
Collapse
|
5
|
Yang H, Govoni M, Kundu A, Galli G. Computational Protocol to Evaluate Electron-Phonon Interactions Within Density Matrix Perturbation Theory. J Chem Theory Comput 2022; 18:6031-6042. [PMID: 36126338 PMCID: PMC9558376 DOI: 10.1021/acs.jctc.2c00579] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
![]()
We present a computational protocol, based on density
matrix perturbation
theory, to obtain non-adiabatic, frequency-dependent electron–phonon
self-energies for molecules and solids. Our approach enables the evaluation
of electron–phonon interaction using hybrid functionals, for
spin-polarized systems, and the computational overhead to include
dynamical and non-adiabatic terms in the evaluation of electron–phonon
self-energies is negligible. We discuss results for molecules, as
well as pristine and defective solids.
Collapse
Affiliation(s)
- Han Yang
- Department of Chemistry, University of Chicago, Chicago, Illinois 60637, United States.,Pritzker School of Molecular Engineering, The University of Chicago, Chicago, Illinois 60637, United States
| | - Marco Govoni
- Pritzker School of Molecular Engineering, The University of Chicago, Chicago, Illinois 60637, United States.,Materials Science Division and Center for Molecular Engineering, Argonne National Laboratory, Lemont, Illinois 60439, United States
| | - Arpan Kundu
- Pritzker School of Molecular Engineering, The University of Chicago, Chicago, Illinois 60637, United States
| | - Giulia Galli
- Department of Chemistry, University of Chicago, Chicago, Illinois 60637, United States.,Pritzker School of Molecular Engineering, The 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
|
6
|
Yu VWZ, Govoni M. GPU Acceleration of Large-Scale Full-Frequency GW Calculations. J Chem Theory Comput 2022; 18:4690-4707. [PMID: 35913080 DOI: 10.1021/acs.jctc.2c00241] [Citation(s) in RCA: 10] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Many-body perturbation theory is a powerful method to simulate electronic excitations in molecules and materials starting from the output of density functional theory calculations. By implementing the theory efficiently so as to run at scale on the latest leadership high-performance computing systems it is possible to extend the scope of GW calculations. We present a GPU acceleration study of the full-frequency GW method as implemented in the WEST code. Excellent performance is achieved through the use of (i) optimized GPU libraries, e.g., cuFFT and cuBLAS, (ii) a hierarchical parallelization strategy that minimizes CPU-CPU, CPU-GPU, and GPU-GPU data transfer operations, (iii) nonblocking MPI communications that overlap with GPU computations, and (iv) mixed precision in selected portions of the code. A series of performance benchmarks has been carried out on leadership high-performance computing systems, showing a substantial speedup of the GPU-accelerated version of WEST with respect to its CPU version. Good strong and weak scaling is demonstrated using up to 25 920 GPUs. Finally, we showcase the capability of the GPU version of WEST for large-scale, full-frequency GW calculations of realistic systems, e.g., a nanostructure, an interface, and a defect, comprising up to 10 368 valence electrons.
Collapse
Affiliation(s)
- Victor Wen-Zhe Yu
- Materials Science Division, Argonne National Laboratory, Lemont, Illinois 60439, United States
| | - Marco Govoni
- Materials Science Division, Argonne National Laboratory, Lemont, Illinois 60439, United States.,Pritzker School of Molecular Engineering, The University of Chicago, Chicago, Illinois 60637, United States
| |
Collapse
|
7
|
Bradbury NC, Nguyen M, Caram J, Neuhauser D. Bethe Salpeter Equation Spectra for Very Large Systems. J Chem Phys 2022; 157:031104. [DOI: 10.1063/5.0100213] [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
We present a highly efficient method for the extraction of optical properties of very large molecules via the Bethe-Salpeter equation. The crutch of this approach is the calculation of the action of the effective Coulombic interaction, $W$, through a stochastic TD Hartree propagation, which uses only 10 stochastic orbitals rather than propagating the full sea of occupied states. This leads to a scaling that is at most cubic in system size, with trivial parallelization of the calculation. We apply this new method to calculate the spectra and electronic density of the dominant excitons of a carbon-nanohoop bound fullerene system with 520 electrons, using less than 4000 core hours.
Collapse
Affiliation(s)
- Nadine Claire Bradbury
- Chemistry and Biochemistry, University of California Los Angeles Department of Chemistry and Biochemistry, United States of America
| | - Minh Nguyen
- Chemistry and Biochemistry, University of California Los Angeles Department of Chemistry and Biochemistry, United States of America
| | - Justin Caram
- UCLA, UCLA Division of Physical Sciences, United States of America
| | - Daniel Neuhauser
- Department of Chemistry and Biochemistry, UCLA, United States of America
| |
Collapse
|
8
|
Umari P. A Fully Linear Response G0W0 Method That Scales Linearly up to Tens of Thousands of Cores. J Phys Chem A 2022; 126:3384-3391. [PMID: 35584057 PMCID: PMC9169063 DOI: 10.1021/acs.jpca.2c01328] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
![]()
We present a G0W0 approach that is based
on the evaluation of the linear response of the actions of the G0 and W0 operators.
In this way we avoid sums over empty one-particle orbitals and do
not have to explicitly develop the screened Coulomb interaction W0 on a dedicated basis. For a given orbital,
the self-energy is found by summing terms relative to a set of points
in the real-space simulation cell. This permits us to easily control
the ratio of the accuracy to the computational cost. A trivial parallelization
strategy allows strong linear scaling up to tens of thousands of computing
cores.
Collapse
Affiliation(s)
- Paolo Umari
- Dipartimento di Fisica e Astronomia, Università di Padova, Padova I-35131, Italy.,CNR-IOM Democritos, Istituto Officina dei Materiali, Consiglio Nazionale delle Ricerchei, Trieste I-34149, Italy
| |
Collapse
|
9
|
Ping Y, Smart TJ. Computational design of quantum defects in two-dimensional materials. NATURE COMPUTATIONAL SCIENCE 2021; 1:646-654. [PMID: 38217204 DOI: 10.1038/s43588-021-00140-w] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 03/07/2021] [Accepted: 09/15/2021] [Indexed: 01/15/2024]
Abstract
Missing atoms or atom substitutions (point defects) in crystal lattices in two-dimensional (2D) materials are potential hosts for emerging quantum technologies, such as single-photon emitters and spin quantum bits (qubits). First-principles-guided design of quantum defects in 2D materials is paving the way for rational spin qubit discovery. Here we discuss the frontier of first-principles theory development and the challenges in predicting the critical physical properties of point defects in 2D materials for quantum information technology, in particular for optoelectronic and spin-optotronic properties. Strong many-body interactions at reduced dimensionality require advanced electronic structure methods beyond mean-field theory. The great challenges for developing theoretical methods that are appropriate for strongly correlated defect states, as well as general approaches for predicting spin relaxation and the decoherence time of spin defects, are yet to be addressed.
Collapse
Affiliation(s)
- Yuan Ping
- Department of Chemistry and Biochemistry, University of California, Santa Cruz, CA, USA.
| | - Tyler J Smart
- Department of Physics, University of California, Santa Cruz, CA, USA
- Lawrence Livermore National Laboratory, Livermore, CA, USA
| |
Collapse
|
10
|
Bhang J, Ma H, Yim D, Galli G, Seo H. First-Principles Predictions of Out-of-Plane Group IV and V Dimers as High-Symmetry, High-Spin Defects in Hexagonal Boron Nitride. ACS APPLIED MATERIALS & INTERFACES 2021; 13:45768-45777. [PMID: 34541839 DOI: 10.1021/acsami.1c16988] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/13/2023]
Abstract
Hexagonal boron nitride (h-BN) has been recently found to host a variety of quantum point defects, which are promising candidates as single-photon sources for solid-state quantum nanophotonic applications. Most recently, optically addressable spin qubits in h-BN have been the focus of intensive research due to their unique potential in quantum computation, communication, and sensing. However, the number of high-symmetry, high-spin defects that are desirable for developing spin qubits in h-BN is highly limited. Here, we combine density functional theory (DFT) and quantum embedding theories to show that out-of-plane XNYi dimer defects (X, Y = C, N, P, and Si) form a new class of stable C3v spin-triplet defects in h-BN. We find that the dimer defects have a robust 3A2 ground state and 3E excited state, both of which are isolated from the h-BN bulk states. We show that 1E and 1A shelving states exist and they are positioned between the 3E and 3A2 states for all the dimer defects considered in this study. To support future experimental identification of the XNYi dimer defects, we provide extensive characterization of the defects in terms of their spin and optical properties. We predict that the zero-phonon line of the spin-triplet XNYi defects lies in the visible range (800 nm to 500 nm). We compute the zero-field splitting of the dimers' spin to range from 1.79 GHz (SiNPi0) to 29.5 GHz (CNNi0). Our results broaden the scope of high-spin defect candidates that would be useful for the development of spin-based solid-state quantum technologies in two-dimensional hexagonal boron nitride.
Collapse
Affiliation(s)
- Jooyong Bhang
- Department of Energy Systems Research and Department of Physics, Ajou University, Suwon, Gyeonggi 16499, Korea
| | - He Ma
- Pritzker School of Molecular Engineering, University of Chicago, Chicago, Illinois 60637, United States
- Department of Chemistry, University of Chicago, Chicago, Illinois 60637, United States
| | - Donggyu Yim
- Department of Energy Systems Research and Department of Physics, Ajou University, Suwon, Gyeonggi 16499, Korea
| | - Giulia Galli
- Pritzker School of Molecular Engineering, University of Chicago, Chicago, Illinois 60637, United States
- Department of Chemistry, University of Chicago, Chicago, Illinois 60637, United States
- Materials Science Division and Center for Molecular Engineering, Argonne National Laboratory, Lemont, Illinois 60439, United States
| | - Hosung Seo
- Department of Energy Systems Research and Department of Physics, Ajou University, Suwon, Gyeonggi 16499, Korea
| |
Collapse
|
11
|
Elliott JD, Mosconi E, De Angelis F, Ambrosetti A, Umari P. Real Space-Real Time Evolution of Excitonic States Based on the Bethe-Salpeter Equation Method. J Phys Chem Lett 2021; 12:7261-7269. [PMID: 34314589 DOI: 10.1021/acs.jpclett.1c01742] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/13/2023]
Abstract
We introduce a method for constructing localized excitations and simulating the real time dynamics of excitons at the Many-Body Perturbation Theory Bethe-Salpeter Equation level. We track, on the femto-seconds scale, electron injection from a photoexcited dye into a semiconducting slab. From the time-dependent many-body wave function we compute the spatial evolution of the electron and of the hole; full electron injection is attained within 5 fs. Time-resolved analysis of the electron density and electron-hole interaction energy hints at a two-step charge transfer mechanism through an intermediary partially injected state. We adopt the Von-Neumann entropy for analyzing how the electron and hole entangle. We find that the excitation of the dye-semiconductor model may be represented by a four-level system and register a decrease in entanglement upon electron injection. At full injection, the electron and the hole exhibit only a small degree of entanglement indicative of pure electron and hole states.
Collapse
Affiliation(s)
- Joshua D Elliott
- School of Chemical Engineering and Analytical Science, University of Manchester, Manchester, United Kingdom
| | - Edoardo Mosconi
- Istituto CNR di Scienze e Technologie Molecolare, via Elce di Sotto, 8, I-06123, Perugia, Italy
| | - Filippo De Angelis
- Istituto CNR di Scienze e Technologie Molecolare, via Elce di Sotto, 8, I-06123, Perugia, Italy
- CompuNet, Istituto Italiano di Technologie, via Morego 30, 16163, Genova, Italy
| | - Alberto Ambrosetti
- Dipartimento di Fisica e Astronomia, University of Padova, Padova, Italy
| | - Paolo Umari
- Dipartimento di Fisica e Astronomia, University of Padova, Padova, Italy
- CNR-IOM DEMOCRITOS, Consiglio Nazionale delle Ricerche-Istituto Officina dei Materiali, c/o SISSA, Via Bonomea 265, 34136, Trieste, Italy
| |
Collapse
|
12
|
Loos PF, Comin M, Blase X, Jacquemin D. Reference Energies for Intramolecular Charge-Transfer Excitations. J Chem Theory Comput 2021; 17:3666-3686. [DOI: 10.1021/acs.jctc.1c00226] [Citation(s) in RCA: 20] [Impact Index Per Article: 6.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/23/2022]
Affiliation(s)
- Pierre-François Loos
- Laboratoire de Chimie et Physique Quantiques, Université de Toulouse, CNRS, UPS, F-31400 Toulouse, France
| | | | - Xavier Blase
- Univ. Grenoble Alpes, CNRS, Inst NEEL, F-38042 Grenoble, France
| | - Denis Jacquemin
- Université de Nantes, CNRS, CEISAM UMR 6230, F-44000 Nantes, France
| |
Collapse
|
13
|
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
|
14
|
Dong SS, Govoni M, Galli G. Machine learning dielectric screening for the simulation of excited state properties of molecules and materials. Chem Sci 2021; 12:4970-4980. [PMID: 34163744 PMCID: PMC8179553 DOI: 10.1039/d1sc00503k] [Citation(s) in RCA: 7] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/27/2021] [Accepted: 02/12/2021] [Indexed: 11/21/2022] Open
Abstract
Accurate and efficient calculations of absorption spectra of molecules and materials are essential for the understanding and rational design of broad classes of systems. Solving the Bethe-Salpeter equation (BSE) for electron-hole pairs usually yields accurate predictions of absorption spectra, but it is computationally expensive, especially if thermal averages of spectra computed for multiple configurations are required. We present a method based on machine learning to evaluate a key quantity entering the definition of absorption spectra: the dielectric screening. We show that our approach yields a model for the screening that is transferable between multiple configurations sampled during first principles molecular dynamics simulations; hence it leads to a substantial improvement in the efficiency of calculations of finite temperature spectra. We obtained computational gains of one to two orders of magnitude for systems with 50 to 500 atoms, including liquids, solids, nanostructures, and solid/liquid interfaces. Importantly, the models of dielectric screening derived here may be used not only in the solution of the BSE but also in developing functionals for time-dependent density functional theory (TDDFT) calculations of homogeneous and heterogeneous systems. Overall, our work provides a strategy to combine machine learning with electronic structure calculations to accelerate first principles simulations of excited-state properties.
Collapse
Affiliation(s)
- Sijia S Dong
- Materials Science Division and Center for Molecular Engineering, Argonne National Laboratory Lemont IL 60439 USA
- Pritzker School of Molecular Engineering, The University of Chicago Chicago IL 60637 USA
| | - Marco Govoni
- Materials Science Division and Center for Molecular Engineering, Argonne National Laboratory Lemont IL 60439 USA
- Pritzker School of Molecular Engineering, The University of Chicago Chicago IL 60637 USA
| | - Giulia Galli
- Materials Science Division and Center for Molecular Engineering, Argonne National Laboratory Lemont IL 60439 USA
- Pritzker School of Molecular Engineering, The University of Chicago Chicago IL 60637 USA
| |
Collapse
|
15
|
Head-Marsden K, Flick J, Ciccarino CJ, Narang P. Quantum Information and Algorithms for Correlated Quantum Matter. Chem Rev 2020; 121:3061-3120. [PMID: 33326218 DOI: 10.1021/acs.chemrev.0c00620] [Citation(s) in RCA: 37] [Impact Index Per Article: 9.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/13/2022]
Abstract
Discoveries in quantum materials, which are characterized by the strongly quantum-mechanical nature of electrons and atoms, have revealed exotic properties that arise from correlations. It is the promise of quantum materials for quantum information science superimposed with the potential of new computational quantum algorithms to discover new quantum materials that inspires this Review. We anticipate that quantum materials to be discovered and developed in the next years will transform the areas of quantum information processing including communication, storage, and computing. Simultaneously, efforts toward developing new quantum algorithmic approaches for quantum simulation and advanced calculation methods for many-body quantum systems enable major advances toward functional quantum materials and their deployment. The advent of quantum computing brings new possibilities for eliminating the exponential complexity that has stymied simulation of correlated quantum systems on high-performance classical computers. Here, we review new algorithms and computational approaches to predict and understand the behavior of correlated quantum matter. The strongly interdisciplinary nature of the topics covered necessitates a common language to integrate ideas from these fields. We aim to provide this common language while weaving together fields across electronic structure theory, quantum electrodynamics, algorithm design, and open quantum systems. Our Review is timely in presenting the state-of-the-art in the field toward algorithms with nonexponential complexity for correlated quantum matter with applications in grand-challenge problems. Looking to the future, at the intersection of quantum information science and algorithms for correlated quantum matter, we envision seminal advances in predicting many-body quantum states and describing excitonic quantum matter and large-scale entangled states, a better understanding of high-temperature superconductivity, and quantifying open quantum system dynamics.
Collapse
Affiliation(s)
- Kade Head-Marsden
- John A. Paulson School of Engineering and Applied Sciences, Harvard University, Cambridge, Massachusetts 02138, United States
| | - Johannes Flick
- Center for Computational Quantum Physics, Flatiron Institute, New York, New York 10010, United States
| | - Christopher J Ciccarino
- John A. Paulson School of Engineering and Applied Sciences, Harvard University, Cambridge, Massachusetts 02138, United States.,Department of Chemistry and Chemical Biology, Harvard University, Cambridge, Massachusetts 02138, United States
| | - Prineha Narang
- John A. Paulson School of Engineering and Applied Sciences, Harvard University, Cambridge, Massachusetts 02138, United States
| |
Collapse
|
16
|
DiRisio RJ, Jones CM, Ma H, Rousseau BJG. Viewpoints on the 2020 Virtual Conference on Theoretical Chemistry. J Phys Chem A 2020; 124:8875-8883. [PMID: 33054223 DOI: 10.1021/acs.jpca.0c08955] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/03/2023]
Affiliation(s)
- Ryan J DiRisio
- Department of Chemistry, University of Washington, Seattle, Washington 98195, United States
| | - Chey M Jones
- Department of Chemistry, Stanford University, Stanford, California 94305, United States
| | - He Ma
- Institute for Molecular engineering, University of Chicago, 5640 S. Ellis Avenue, Chicago, Illinois 60637, United States
| | - Benjamin J G Rousseau
- Department of Chemistry, Yale University, 225 Prospect Street, New Haven, Connecticut 06520, United States
| |
Collapse
|
17
|
Blase X, Duchemin I, Jacquemin D, Loos PF. The Bethe-Salpeter Equation Formalism: From Physics to Chemistry. J Phys Chem Lett 2020; 11:7371-7382. [PMID: 32787315 DOI: 10.1021/acs.jpclett.0c01875] [Citation(s) in RCA: 70] [Impact Index Per Article: 17.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/11/2023]
Abstract
The Bethe-Salpeter equation (BSE) formalism is steadily asserting itself as a new efficient and accurate tool in the ensemble of computational methods available to chemists in order to predict optical excitations in molecular systems. In particular, the combination of the so-called GW approximation, giving access to reliable ionization energies and electron affinities, and the BSE formalism, able to model UV/vis spectra, has shown to provide accurate singlet excitation energies with a typical error of 0.1-0.3 eV. With a similar computational cost as time-dependent density-functional theory (TD-DFT), BSE is able to provide an accuracy on par with the most accurate global and range-separated hybrid functionals without the unsettling choice of the exchange-correlation functional, resolving further known issues (e.g., charge-transfer excitations). In this Perspective, we provide a historical overview of BSE, with a particular focus on its condensed-matter roots. We also propose a critical review of its strengths and weaknesses in different chemical situations.
Collapse
Affiliation(s)
- Xavier Blase
- Université Grenoble Alpes, CNRS, Institut NEEL, F-38042 Grenoble, France
| | - Ivan Duchemin
- Université Grenoble Alpes, CEA, IRIG-MEM-L Sim, 38054 Grenoble, France
| | - Denis Jacquemin
- Université de Nantes, CNRS, CEISAM UMR 6230, F-44000 Nantes, France
| | - Pierre-François Loos
- Laboratoire de Chimie et Physique Quantiques, Université de Toulouse, CNRS, UPS, Toulouse, France
| |
Collapse
|
18
|
Ma H, Wang W, Kim S, Cheng MH, Govoni M, Galli G. PyCDFT: A Python package for constrained density functional theory. J Comput Chem 2020; 41:1859-1867. [PMID: 32497321 DOI: 10.1002/jcc.26354] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/16/2020] [Revised: 05/14/2020] [Accepted: 05/15/2020] [Indexed: 12/12/2022]
Abstract
We present PyCDFT, a Python package to compute diabatic states using constrained density functional theory (CDFT). PyCDFT provides an object-oriented, customizable implementation of CDFT, and allows for both single-point self-consistent-field calculations and geometry optimizations. PyCDFT is designed to interface with existing density functional theory (DFT) codes to perform CDFT calculations where constraint potentials are added to the Kohn-Sham Hamiltonian. Here, we demonstrate the use of PyCDFT by performing calculations with a massively parallel first-principles molecular dynamics code, Qbox, and we benchmark its accuracy by computing the electronic coupling between diabatic states for a set of organic molecules. We show that PyCDFT yields results in agreement with existing implementations and is a robust and flexible package for performing CDFT calculations. The program is available at https://dx.doi.org/10.5281/zenodo.3821097.
Collapse
Affiliation(s)
- He Ma
- Department of Chemistry, University of Chicago, Chicago, Illinois, USA.,Materials Science Division and Center for Molecular Engineering, Argonne National Laboratory, Lemont, Illinois, USA
| | - Wennie Wang
- Pritzker School of Molecular Engineering, University of Chicago, Chicago, Illinois, USA
| | - Siyoung Kim
- Pritzker School of Molecular Engineering, University of Chicago, Chicago, Illinois, USA
| | - Man-Hin Cheng
- Pritzker School of Molecular Engineering, University of Chicago, Chicago, Illinois, USA
| | - Marco Govoni
- Materials Science Division and Center for Molecular Engineering, Argonne National Laboratory, Lemont, Illinois, USA.,Pritzker School of Molecular Engineering, University of Chicago, Chicago, Illinois, USA
| | - Giulia Galli
- Department of Chemistry, University of Chicago, Chicago, Illinois, USA.,Materials Science Division and Center for Molecular Engineering, Argonne National Laboratory, Lemont, Illinois, USA.,Pritzker School of Molecular Engineering, University of Chicago, Chicago, Illinois, USA
| |
Collapse
|
19
|
Ma H, Govoni M, Gygi F, Galli G. Correction to A Finite-Field Approach for GW Calculations beyond the Random Phase Approximation. J Chem Theory Comput 2020; 16:2877-2879. [PMID: 32181645 DOI: 10.1021/acs.jctc.0c00221] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
|
20
|
Ma H, Sheng N, Govoni M, Galli G. First-principles studies of strongly correlated states in defect spin qubits in diamond. Phys Chem Chem Phys 2020; 22:25522-25527. [PMID: 33084673 DOI: 10.1039/d0cp04585c] [Citation(s) in RCA: 12] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/20/2022]
Abstract
Using a recently developed quantum embedding theory, we present first-principles calculations of strongly correlated states of spin defects in diamond. Using this theory, effective Hamiltonians are constructed, which can be solved by classical and quantum computers; the latter promise a much more favorable scaling as a function of system size than the former. In particular, we report a study on the neutral group-IV vacancy complexes in diamond, and we discuss their strongly correlated spin-singlet and spin-triplet excited states. Our results provide valuable predictions for experiments aimed at optical manipulation of these defects for quantum information technology applications.
Collapse
Affiliation(s)
- He Ma
- Department of Chemistry, University of Chicago, Chicago, IL 60637, USA
| | | | | | | |
Collapse
|