1
|
Alvertis AM, Williams-Young DB, Bruneval F, Neaton JB. Influence of Electronic Correlations on Electron-Phonon Interactions of Molecular Systems with the GW and Coupled Cluster Methods. J Chem Theory Comput 2024; 20:6175-6183. [PMID: 38954597 DOI: 10.1021/acs.jctc.4c00327] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 07/04/2024]
Abstract
Electron-phonon interactions are of great importance to a variety of physical phenomena, and their accurate description is an important goal for first-principles calculations. Isolated examples of materials and molecular systems have emerged where electron-phonon coupling is enhanced over density functional theory (DFT) when using the Green's-function-based ab initio GW method, which provides a more accurate description of electronic correlations. It is, however, unclear how general this enhancement is and how employing high-end quantum chemistry methods, which further improve the description of electronic correlations, might further alter electron-phonon interactions over GW or DFT. Here, we address these questions by computing the renormalization of the highest occupied molecular orbital energies of Thiel's set of organic molecules by harmonic vibrations using DFT, GW, and equation-of-motion coupled-cluster calculations. We find that, depending on the amount of exact exchange included in the DFT starting point, GW can increase the magnitude of the electron-phonon coupling across Thiel's set of molecules by an average factor of 1.1-1.8 compared to the underlying DFT, while equation-of-motion coupled-cluster leads to an increase of 1.4-2. The electron-phonon coupling predicted with the ab initio GW method is generally in much closer agreement to coupled cluster values compared to DFT, establishing GW as a promising route for accurately computing electron-phonon phenomena in molecules and beyond at a much lower computational cost than higher-end quantum chemistry techniques.
Collapse
Affiliation(s)
- Antonios M Alvertis
- KBR, Inc., NASA Ames Research Center, Moffett Field, California 94035, United States
- Materials Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, United States
| | - David B Williams-Young
- Applied Mathematics and Computational Research Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, United States
| | - Fabien Bruneval
- Université Paris-Saclay, CEA, Service de Corrosion et de Comportement des Matériaux, SRMP, 91191 Gif-sur-Yvette, France
| | - Jeffrey B Neaton
- Materials Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, United States
- Department of Physics, University of California Berkeley, Berkeley, California 94720, United States
- Kavli Energy NanoScience Institute at Berkeley, Berkeley, California 94720, United States
| |
Collapse
|
2
|
Kundu A, Galli G. Quantum Vibronic Effects on the Excitation Energies of the Nitrogen-Vacancy Center in Diamond. J Phys Chem Lett 2024; 15:802-810. [PMID: 38232151 DOI: 10.1021/acs.jpclett.3c03269] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/19/2024]
Abstract
We investigated the impact of quantum vibronic coupling on the electronic properties of solid-state spin defects using stochastic methods and first-principles molecular dynamics with a quantum thermostat. Focusing on the negatively charged nitrogen-vacancy center in diamond as an exemplary case, we found a significant dynamic Jahn-Teller splitting of the doubly degenerate single-particle levels within the diamond's band gap, even at 0 K, with a magnitude exceeding 180 meV. This pronounced splitting leads to substantial renormalizations of these levels and, subsequently, of the vertical excitation energies of the doubly degenerate singlet and triplet excited states. Our findings underscore the pressing need to incorporate quantum vibronic effects into first-principles calculations, particularly when comparing computed vertical excitation energies with experimental data. Our study also reveals the efficiency of stochastic thermal line sampling for studying phonon renormalizations of solid-state spin defects.
Collapse
Affiliation(s)
- Arpan Kundu
- Pritzker School of Molecular Engineering, The University of Chicago, Chicago, Illinois 60637, United States
| | - Giulia Galli
- Pritzker School of Molecular Engineering, The University of Chicago, Chicago, Illinois 60637, United States
- Department of Chemistry, 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
|
3
|
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
|
4
|
Kundu A, Galli G. Quantum Vibronic Effects on the Electronic Properties of Molecular Crystals. J Chem Theory Comput 2023. [PMID: 37378491 DOI: 10.1021/acs.jctc.3c00424] [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/29/2023]
Abstract
We present a study of molecular crystals, focused on the effect of nuclear quantum motion and anharmonicity on their electronic properties. We consider a system composed of relatively rigid molecules, a diamondoid crystal, and one composed of floppier molecules, NAI-DMAC, a thermally activated delayed fluorescence compound. We compute fundamental electronic gaps at the density functional theory (DFT) level of theory, with the Perdew-Burke-Erzenhof (PBE) and strongly constrained and approximately normed (SCAN) functionals, by coupling first-principles molecular dynamics with a nuclear quantum thermostat. We find a sizable zero-point renormalization (ZPR) of the band gaps, which is much larger in the case of diamondoids (0.6 eV) than for NAI-DMAC (0.22 eV). We show that the frozen phonon (FP) approximation, which neglects intermolecular anharmonic effects, leads to a large error (∼50%) in the calculation of the band gap ZPR. Instead, when using a stochastic method, we obtain results in good agreement with those of our quantum simulations for the diamondoid crystal. However, the agreement is worse for NAI-DMAC where intramolecular anharmonicities contribute to the ZPR. Our results highlight the importance of accurately including nuclear and anharmonic quantum effects to predict the electronic properties of molecular crystals.
Collapse
Affiliation(s)
- Arpan Kundu
- Pritzker School of Molecular Engineering, The University of Chicago, Chicago, Illinois 60637, United States
| | - Giulia Galli
- Pritzker School of Molecular Engineering, The 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
| |
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
|