skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.

Explore Research Products in the PAR
It may take a few hours for recently added research products to appear in PAR search results.


Title: All electron GW with linearized augmented plane waves for metals and semiconductors
GW approximation is one of the most popular parameter-free many-body methods that go beyond the limitations of the standard density functional theory (DFT) to determine the excitation spectra for moderately correlated materials and in particular the semiconductors. It is also the first step in developing the diagrammatic Monte Carlo method into an electronic structure tool, which would offer a numerically exact solution to the solid-state problem. While most electronic structure packages offer support for GW calculations for band-insulating materials, the level of support for metallic systems is somewhat limited. This limitation can be partly attributed to the relatively minor differences often observed between GW and DFT results in treating metallic systems, which is not expected to persist to higher orders in perturbation theory. Describing metals within the GW framework presents a challenge, as it requires accurate resolution of Fermi surface singularities, which, in turn, calls for a dense momentum mesh. Here we implement the GW algorithm within the all-electron Linear Augmented Plane Wave framework, where we pay special attention to the metallic systems, the convergence with respect to momentum mesh, and proper treatment of the deep laying core states, as needed for the future variational diagrammatic Monte Carlo implementation. Our improved algorithm for resolving Fermi surface singularities allows us a stable and accurate analytic continuation of imaginary axis data, which is carried out for GW excitation spectra throughout the Brillouin zone in both the metallic and insulating materials and is compared to numerically more stable contour deformation integration technique. We compute band structures for elemental metallic systems Li, Na, and Mg as well as for various narrow and wide bandgap insulators such as Si, BN, SiC, MgO, LiF, ZnS, and CdS and compare our results with previous GW calculations and available experiments data. Our results are in good agreement with the available literature. Thus our software allows users to compute full bandstructures for metals and insulators using all-electron potential without downfolding to Wannier orbital basis.  more » « less
Award ID(s):
2233892 2311557
PAR ID:
10511761
Author(s) / Creator(s):
;
Publisher / Repository:
Elsevier
Date Published:
Journal Name:
Computer Physics Communications
Volume:
295
Issue:
C
ISSN:
0010-4655
Page Range / eLocation ID:
108986
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. ; ; (, npj Computational Materials)
    Abstract Transition metald-electron oxides with an odd number of electrons per unit cell are expected to form metals with partially occupied energy bands, but exhibit in fact a range of behaviors, being either insulators, or metals, or having insulator-metal transitions. Traditional explanations involved predominantly electron-electron interactions in fixed structural symmetry. The present work focuses instead on the role of symmetry breaking local structural motifs. Viewing the previously observed V-V dimerization in VO2as a continuous knob, reveals in density functional calculations the splitting of an isolated flat band from the broad conduction band. This leads past a critical percent dimerization to the formation of the insulating phase while lowering the total energy. In VO2this transition is found to have a rather low energy barrier approaching the thermal energy at room temperature, suggesting energy-efficient switching in neuromorphic computing. Interestingly, sufficient V-V dimerization suppresses magnetism, leading to the nonmagnetic insulating state, whereas magnetism appears when dimerization is reduced, forming a metallic state. This study opens the way to design novel functional quantum materials with symmetry breaking-induced flat bands. 
    more » « less
  2. ; ; ; ; ; ; (, Monthly Notices of the Royal Astronomical Society)
    ABSTRACT The association of GRB170817A with a binary neutron star (BNS) merger has revealed that BNSs produce at least a fraction of short gamma-ray bursts (SGRBs). As gravitational wave (GW) detectors push their horizons, it is important to assess coupled electromagnetic (EM)/GW probabilities and maximize observational prospects. Here, we perform BNS population synthesis calculations with the code mobse, seeding the binaries in galaxies at three representative redshifts, $$z$$ = 0.01, 0.1, and 1 of the Illustris TNG50 simulation. The binaries are evolved and their locations numerically tracked in the host galactic potentials until merger. Adopting the microphysics parameters of GRB170817A, we numerically compute the broad-band light curves of jets from BNS mergers, with the afterglow brightness dependent on the local medium density at the merger site. We perform Monte Carlo simulations of the resulting EM population assuming either a random viewing angle with respect to the jet, or a jet aligned with the orbital angular momentum of the binary, which biases the viewing angle probability for GW-triggered events. We find a gamma-ray detection probability of $$\sim\!2{{\rm per\ cent}},10{{\rm per\ cent}},\mathrm{and}\ 40{{\rm per\ cent}}$$ for BNSs at $$z$$ = 1, 0.1, and 0.01, respectively, for the random case, rising to $$\sim\!75{{\rm per\ cent}}$$ for the $$z$$ = 0.01, GW-triggered aligned case. Afterglow detection probabilities of GW-triggered BNS mergers vary in the range of $$\sim \! 0.3 \!-\! 0.5{{\rm per\ cent}}$$, with higher values for aligned jets, and are comparable across the high- and low-energy bands, unlike gamma-ray-triggered events (cosmological SGRBs) which are significantly brighter at higher energies. We further quantify observational biases with respect to host galaxy masses. 
    more » « less
  3. ; (, Scientific Reports)
    Abstract We calculate the single-particle excitation spectrum and the Landau liquid parameters for the archetypal model of solids, the three-dimensional uniform electron gas, with the variational diagrammatic Monte Carlo method, which gives numerically controlled results without systematic error. In the metallic range of density, we establish benchmark values for the wave-function renormalization factor Z , the effective mass $$m^*/m$$ m ∗ / m , and the Landau parameters $$F_0^s$$ F 0 s and $$F_0^a$$ F 0 a with unprecedented accuracy, and we resolve the long-standing puzzle of non-monotonic dependence of mass on density. We also exclude the possibility that experimentally measured large reduction of bandwidth in Na metal can originate from the charge and spin fluctuations contained in the model of the uniform electron gas. 
    more » « less
  4. ; ; ; ; ; ; ; (, ChemRxiv)
    The excited-state properties of molecular crystals are important for applications in organic electronic devices. The GW approximation and Bethe-Salpeter equation (GW+BSE) is the state-of-the-art method for calculating the excited-state properties of crystalline solids with periodic boundary conditions. We present the PAH101 dataset of GW +BSE calculations for 101 molecular crystals of polycyclic aromatic hydrocarbons (PAHs) with up to ∼500 atoms in the unit cell. The data records include the GW quasiparticle band structure, the fundamental band gap, the static dielectric constant, the first singlet exciton energy (optical gap), the first triplet exciton energy, the dielectric function, and optical absorption spectra for light polarized along the three lattice vectors. In addition, the dataset includes the density functional theory (DFT) single-molecule and crystal features used in Liu et al. [npj Computational Materials, 8, 70 (2022)]. We envision the dataset being used to (i) identify correlations between DFT and GW +BSE quantities, (ii) discover materials with desired electronic/ optical properties in the dataset itself, and (iii) train machine-learned models to help in materials discovery efforts. We provide examples to illustrate these three use cases. 
    more » « less
  5. ; ; ; ; ; ; ; (, ChemRxiv)
    The excited-state properties of molecular crystals are important for applications in organic electronic devices. The GW approximation and Bethe-Salpeter equation (GW +BSE) is the state-of-the-art method for calculating the excited-state properties of crystalline solids with periodic boundary conditions. We present the PAH101 dataset of GW +BSE calculations for 101 molecular crystals of polycyclic aromatic hydrocarbons (PAHs) with up to ∼500 atoms in the unit cell. The data records include the GW quasiparticle band structure, the fundamental band gap, the static dielectric constant, the first singlet exciton energy (optical gap), the first triplet exciton energy, the dielectric function, and optical absorption spectra for light polarized along the three lattice vectors. In addition, the dataset includes the density functional theory (DFT) single-molecule and crystal features used in Liu et al. [npj Computational Materials, 8, 70 (2022)]. We envision the dataset being used to (i) identify correlations between DFT and GW +BSE quantities, (ii) discover materials with desired electronic/ optical properties in the dataset itself, and (iii) train machine-learned models to help in materials discovery efforts. We provide examples to illustrate these three use cases. 
    more » « less
  • Citation Formats
  • MLA
  • APA
  • Chicago
  • BibTeX
  • Save / Share this Record
  • Send to Email