1
|
Graml M, Zollner K, Hernangómez-Pérez D, Faria Junior PE, Wilhelm J. Low-Scaling GW Algorithm Applied to Twisted Transition-Metal Dichalcogenide Heterobilayers. J Chem Theory Comput 2024; 20:2202-2208. [PMID: 38353944 PMCID: PMC10938508 DOI: 10.1021/acs.jctc.3c01230] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/06/2023] [Revised: 12/30/2023] [Accepted: 01/08/2024] [Indexed: 03/13/2024]
Abstract
The GW method is widely used for calculating the electronic band structure of materials. The high computational cost of GW algorithms prohibits their application to many systems of interest. We present a periodic, low-scaling, and highly efficient GW algorithm that benefits from the locality of the Gaussian basis and the polarizability. The algorithm enables G0W0 calculations on a MoSe2/WS2 bilayer with 984 atoms per unit cell, in 42 h using 1536 cores. This is 4 orders of magnitude faster than a plane-wave G0W0 algorithm, allowing for unprecedented computational studies of electronic excitations at the nanoscale.
Collapse
Affiliation(s)
- Maximilian Graml
- Institute
of Theoretical Physics, University of Regensburg, 93053 Regensburg, Germany
- Regensburg
Center for Ultrafast Nanoscopy (RUN), University
of Regensburg, 93053 Regensburg, Germany
| | - Klaus Zollner
- Institute
of Theoretical Physics, University of Regensburg, 93053 Regensburg, Germany
| | - Daniel Hernangómez-Pérez
- Department
of Molecular Chemistry and Materials Science, Weizmann Institute of Science, Rehovot 7610001, Israel
| | - Paulo E. Faria Junior
- Institute
of Theoretical Physics, University of Regensburg, 93053 Regensburg, Germany
| | - Jan Wilhelm
- Institute
of Theoretical Physics, University of Regensburg, 93053 Regensburg, Germany
- Regensburg
Center for Ultrafast Nanoscopy (RUN), University
of Regensburg, 93053 Regensburg, Germany
| |
Collapse
|
2
|
Alhassan A, Yu M. Insight into the stacking effect on shifted patterns of bilayer phosphorene: a comprehensive first-principles study. NANOTECHNOLOGY 2024; 35:155701. [PMID: 38198734 DOI: 10.1088/1361-6528/ad1d14] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 07/18/2023] [Accepted: 01/10/2024] [Indexed: 01/12/2024]
Abstract
It is crucial to deeply understand how the interlayer interaction acts on controlling the structural and electronic properties of shifted patterns of bilayer phosphorene. A comprehensive first-principles study on the bilayer phosphorene through relative translation along different directions has revealed that there is a direct correlation between the potential energy surface and the interlayer equilibrium distance. The shorter the interlayer distance, the lower the potential energy surface. The shifted patterns with the most stable state, the metastable state, and the transition state (with energy barrier of ∼1.3 meV/atom) were found associated with the AB, the Aδ, and the TS stacking configurations, respectively. The high energy barriers, on the other hand, are ∼9.3 meV/atom at the AA stacking configuration along the zigzag pathway, ∼5.3 meV/atom at the AB' stacking configuration along the armchair pathway, and ∼11.2 meV/atom at the AA' stacking configuration along the diagonal pathway, respectively. The character of electronic bandgap with respect to the shifting shows an anisotropic behavior (with the value of 0.69-1.22 eV). A transition from the indirect to the direct bandgap occurs under the shifting, implying a tunable bandgap by stacking engineering. Furthermore, the orbital hybridization at the interfacial region induces a redistribution of the net charge (∼0.002-0.011e) associated with the relative shifting between layers, leading to a strong polarization with stripe-like electron depletion near the lone pairs and accumulation in the middle of the interfacial region. It is expected that such interesting findings will provide a fundamental reference to deeply understand and analyze the complex local structural and electronic properties of twisted bilayer phosphorene and will make the shifted patterns of bilayer phosphorene promising for nanoelectronics as versatile shiftronics materials.
Collapse
Affiliation(s)
- Aswad Alhassan
- Department of Physics and Astronomy, University of Louisville, Louisville, KY 40292, United States of America
| | - Ming Yu
- Department of Physics and Astronomy, University of Louisville, Louisville, KY 40292, United States of America
| |
Collapse
|
3
|
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
|
4
|
Weng G, Pang A, Vlček V. Spatial Decay and Limits of Quantum Solute-Solvent Interactions. J Phys Chem Lett 2023; 14:2473-2480. [PMID: 36867592 DOI: 10.1021/acs.jpclett.3c00208] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/18/2023]
Abstract
Molecular excitations in the liquid-phase environment are renormalized by the surrounding solvent molecules. Herein, we employ the GW approximation to investigate the solvation effects on the ionization energy of phenol in various solvent environments. The electronic effects differ by up to 0.4 eV among the five investigated solvents. This difference depends on both the macroscopic solvent polarizability and the spatial decay of the solvation effects. The latter is probed by separating the electronic subspace and the GW correlation self-energy into fragments. The fragment correlation energy decays with increasing intermolecular distance and vanishes at ∼9 Å, and this pattern is independent of the type of solvent environment. The 9 Å cutoff defines an effective interacting volume within which the ionization energy shift per solvent molecule is proportional to the macroscopic solvent polarizability. Finally, we propose a simple model for computing the ionization energies of molecules in an arbitrary solvent environment.
Collapse
Affiliation(s)
- Guorong Weng
- Department of Chemistry and Biochemistry, University of California, Santa Barbara, California 93106-9510, United States
| | - Amanda Pang
- Department of Chemistry and Biochemistry, University of California, Santa Barbara, California 93106-9510, United States
| | - Vojtěch Vlček
- Department of Chemistry and Biochemistry, University of California, Santa Barbara, California 93106-9510, United States
- Materials Department, University of California, Santa Barbara, California 93106-5050, United States
| |
Collapse
|
5
|
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
|
6
|
Baer R, Neuhauser D, Rabani E. Stochastic Vector Techniques in Ground-State Electronic Structure. Annu Rev Phys Chem 2022; 73:255-272. [PMID: 35081326 DOI: 10.1146/annurev-physchem-090519-045916] [Citation(s) in RCA: 6] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/09/2022]
Abstract
We review a suite of stochastic vector computational approaches for studying the electronic structure of extended condensed matter systems. These techniques help reduce algorithmic complexity, facilitate efficient parallelization, simplify computational tasks, accelerate calculations, and diminish memory requirements. While their scope is vast, we limit our study to ground-state and finite temperature density functional theory (DFT) and second-order perturbation theory. More advanced topics, such as quasiparticle (charge) and optical (neutral) excitations and higher-order processes, are covered elsewhere. We start by explaining how to use stochastic vectors in computations, characterizing the associated statistical errors. Next, we show how to estimate the electron density in DFT and discuss highly effective techniques to reduce statistical errors. Finally, we review the use of stochastic vector techniques for calculating correlation energies within the second-order Møller-Plesset perturbation theory and its finite temperature variational form. Example calculation results are presented and used to demonstrate the efficacy of the methods. Expected final online publication date for the Annual Review of Physical Chemistry, Volume 73 is April 2022. Please see http://www.annualreviews.org/page/journal/pubdates for revised estimates.
Collapse
Affiliation(s)
- Roi Baer
- Fritz Haber Center of Molecular Dynamics and Institute of Chemistry, The Hebrew University of Jerusalem, Jerusalem, Israel;
| | - Daniel Neuhauser
- Department of Chemistry and Biochemistry, University of California, Los Angeles, California, USA;
| | - Eran Rabani
- Department of Chemistry, University of California, Berkeley, California, USA; .,Materials Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, California, USA.,The Raymond and Beverly Sackler Center of Computational Molecular and Materials Science, Tel Aviv University, Tel Aviv, Israel
| |
Collapse
|
7
|
Zhu Y, Long R. Density Functional Theory Half-Electron Self-Energy Correction for Fast and Accurate Nonadiabatic Molecular Dynamics. J Phys Chem Lett 2021; 12:10886-10892. [PMID: 34730966 DOI: 10.1021/acs.jpclett.1c03077] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/13/2023]
Abstract
The nonadiabatic (NA) process is crucial to photochemistry and photophysics and requires an atomistic understanding. However, conventional NA molecular dynamics (MD) for condensed-phase materials on the nanoscale are generally limited to the semilocal exchange-correlation functional, which suffers from the bandgap and thus NA coupling (NAC) problems. We consider TiO2 and a black phosphorus monolayer as two prototypical systems, perform NA-MD simulations of nonradiative electron-hole recombination, and demonstrate for the first time that density functional theory (DFT) half-electron self-energy correction can reproduce the bandgap, effective masses of carriers, luminescence line widths, NAC, and excited-state lifetimes of the two systems at the hybrid functional level while the computational cost remains at that of the Predew-Burke-Ernzerhof functional. Our study indicates that the DFT-1/2 method can greatly accelerate NA-MD simulations while maintaining the accuracy of the hybrid functional, providing an advantage for studying photoexcitation dynamics for large-scale condensed-phase materials.
Collapse
Affiliation(s)
- Yonghao Zhu
- College of Chemistry, Key Laboratory of Theoretical & Computational Photochemistry of Ministry of Education, Beijing Normal University, Beijing 100875, People's Republic of China
| | - Run Long
- College of Chemistry, Key Laboratory of Theoretical & Computational Photochemistry of Ministry of Education, Beijing Normal University, Beijing 100875, People's Republic of China
| |
Collapse
|
8
|
Twisted bilayer arsenene sheets as a chemical sensor for toluene and M-xylene vapours - A DFT investigation. J Mol Graph Model 2021; 109:108034. [PMID: 34555724 DOI: 10.1016/j.jmgm.2021.108034] [Citation(s) in RCA: 12] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/17/2021] [Revised: 09/02/2021] [Accepted: 09/09/2021] [Indexed: 12/14/2022]
Abstract
2D (two-dimensional) materials are emerging in today's world. Among the 2D materials, arsenene sheets are prominently used as chemical and biosensors. In the present work, the twisted bilayer arsenene sheets (TB-AsNS) are used to adsorb toluene and M-xylene vapours. Moreover, the band gap of pristine TB-AsNS is calculated to be 0.437 eV. Besides, the surface adsorption of toluene and M-xylene vapours modify the electronic properties of TB-AsNS noticed from the band structure, density of states, and electron density difference diagrams. The surface assimilation of target toluene and M-xylene on TB-AsNS falls in the physisorption regime facilitating the adsorption and desorption of molecules. Also, the charge transfer analysis infers that TB-AsNS acts as acceptor and target molecules play as donors. The findings support that TB-AsNS can be used as a sensing medium towards M-xylene and toluene.
Collapse
|
9
|
Weng G, Vlček V. Efficient treatment of molecular excitations in the liquid phase environment via stochastic many-body theory. J Chem Phys 2021; 155:054104. [PMID: 34364336 DOI: 10.1063/5.0058410] [Citation(s) in RCA: 10] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Accurate predictions of charge excitation energies of molecules in the disordered condensed phase are central to the chemical reactivity, stability, and optoelectronic properties of molecules and critically depend on the specific environment. Herein, we develop a stochastic GW method for calculating these charge excitation energies. The approach employs maximally localized electronic states to define the electronic subspace of a molecule and the rest of the system, both of which are randomly sampled. We test the method on three solute-solvent systems: phenol, thymine, and phenylalanine in water. The results are in excellent agreement with the previous high-level calculations and available experimental data. The stochastic calculations for supercells containing up to 1000 electrons representing the solvated systems are inexpensive and require ≤1000 central processing unit hrs. We find that the coupling with the environment accounts for ∼40% of the total correlation energy. The solvent-to-solute feedback mechanism incorporated in the molecular correlation term causes up to 0.6 eV destabilization of the quasiparticle energy. Simulated photo-emission spectra exhibit red shifts, state-degeneracy lifting, and lifetime shortening. Our method provides an efficient approach for an accurate study of excitations of large molecules in realistic condensed phase environments.
Collapse
Affiliation(s)
- Guorong Weng
- 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
|
10
|
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
|
11
|
Wilhelm J, Seewald P, Golze D. Low-Scaling GW with Benchmark Accuracy and Application to Phosphorene Nanosheets. J Chem Theory Comput 2021; 17:1662-1677. [DOI: 10.1021/acs.jctc.0c01282] [Citation(s) in RCA: 19] [Impact Index Per Article: 6.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/11/2022]
Affiliation(s)
- Jan Wilhelm
- Institute of Theoretical Physics, University of Regensburg, D-93053 Regensburg, Germany
| | - Patrick Seewald
- Department of Chemistry, University of Zurich, CH-8057 Zurich, Switzerland
| | - Dorothea Golze
- Department of Applied Physics, Aalto University, FI-00076 Aalto, Finland
| |
Collapse
|
12
|
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
|
13
|
Romanova M, Vlček V. Decomposition and embedding in the stochastic GW self-energy. J Chem Phys 2020; 153:134103. [DOI: 10.1063/5.0020430] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Mariya Romanova
- 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
|