skip to main content

Title: Combined subsampling and analytical integration for efficient large-scale GW calculations for 2D systems

Accurate and efficient predictions of the quasiparticle properties of complex materials remain a major challenge due to the convergence issue and the unfavorable scaling of the computational cost with respect to the system size. QuasiparticleGWcalculations for two-dimensional (2D) materials are especially difficult. The unusual analytical behaviors of the dielectric screening and the electron self-energy of 2D materials make the conventional Brillouin zone (BZ) integration approach rather inefficient and require an extremely densek-grid to properly converge the calculated quasiparticle energies. In this work, we present a combined nonuniform subsampling and analytical integration method that can drastically improve the efficiency of the BZ integration in 2DGWcalculations. Our work is distinguished from previous work in that, instead of focusing on the intricate dielectric matrix or the screened Coulomb interaction matrix, we exploit the analytical behavior of various terms of the convolved self-energy Σ(q) in the smallqlimit. This method, when combined with another acceleratedGWmethod that we developed recently, can drastically speed up (by over three orders of magnitude)GWcalculations for 2D materials. Our method allows fully convergedGWcalculations for complex 2D systems at a fraction of computational cost, facilitating future high throughput screening of the quasiparticle properties of 2D semiconductors for various applications. To demonstrate more » the capability and performance of our new method, we have carried out fully convergedGWcalculations for monolayer C2N, a recently discovered 2D material with a large unit cell, and investigate its quasiparticle band structure in detail.

« less
; ; ; ; ; ;
Publication Date:
Journal Name:
npj Computational Materials
Nature Publishing Group
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract

    The discovery of high-dielectric materials is crucial to increasing the efficiency of electronic devices and batteries. Here, we report three previously unexplored materials with very high dielectric constants (69 < ϵ < 101) and large band gaps (2.9 < Eg(eV) < 5.5) obtained by screening materials databases using statistical optimization algorithms aided by artificial neural networks (ANN). Two of these new dielectrics are mixed-anion compounds (Eu5SiCl6O4and HoClO) and are shown to be thermodynamically stable against common semiconductors via phase diagram analysis. We also uncovered four other materials with relatively large dielectric constants (20 < ϵ < 40) and band gaps (2.3 < Eg(eV) < 2.7). While the ANN training-data are obtained from the Materials Project, the search-space consists of materials from the Open Quantum Materials Database (OQMD)—demonstrating a successful implementation of cross-database materials design. Overall, we report the dielectric properties of 17 materials calculated using ab initio calculations, that were selected in our design workflow. The dielectric materials with high-dielectric properties predicted in this work open up further experimental research opportunities.

  2. Abstract

    Complete theoretical understanding of the most complex superconductors requires a detailed knowledge of the symmetry of the superconducting energy-gap$${\mathrm{{\Delta}}}_{\mathbf{k}}^\alpha$$Δkα, for all momentakon the Fermi surface of every bandα. While there are a variety of techniques for determining$$|{\mathrm{{\Delta}}}_{\mathbf{k}}^\alpha |$$Δkα, no general method existed to measure the signed values of$${\mathrm{{\Delta}}}_{\mathbf{k}}^\alpha$$Δkα. Recently, however, a technique based on phase-resolved visualization of superconducting quasiparticle interference (QPI) patterns, centered on a single non-magnetic impurity atom, was introduced. In principle, energy-resolved and phase-resolved Fourier analysis of these images identifies wavevectors connecting allk-space regions where$${\mathrm{{\Delta}}}_{\mathbf{k}}^\alpha$$Δkαhas the same or opposite sign. But use of a single isolated impurity atom, from whose precise location the spatial phase of the scattering interference pattern must be measured, is technically difficult. Here we introduce a generalization of this approach for use with multiple impurity atoms, and demonstrate its validity by comparing the$${\mathrm{{\Delta}}}_{\mathbf{k}}^\alpha$$Δkαit generates to the$${\mathrm{{\Delta}}}_{\mathbf{k}}^\alpha$$Δkαdetermined from single-atom scattering in FeSe where s±energy-gap symmetry is established. Finally, to exemplify utility, we use the multi-atom technique on LiFeAs and find scattering interference between the hole-like and electron-like pockets as predicted for$${\mathrm{{\Delta}}}_{\mathbf{k}}^\alpha$$Δkαof opposite sign.

  3. Abstract

    For more than three decades, nearly free-electron elemental metals have been a topic of debate because the computed bandwidths are significantly wider in the local density approximation to density-functional theory (DFT) than indicated by angle-resolved photoemission (ARPES) experiments. Here, we systematically investigate this using first principles calculations for alkali and alkaline-earth metals using DFT and various beyond-DFT methods such as meta-GGA, G0W0, hybrid functionals (YS-PBE0, B3LYP), and LDA + eDMFT. We find that the static non-local exchange, as partly included in the hybrid functionals, significantly increase the bandwidths even compared to LDA, while the G0W0bands are only slightly narrower than in LDA. The agreement with the ARPES is best when the local approximation to the self-energy is used in the LDA + eDMFT method. We infer that even moderately correlated systems with partially occupiedsorbitals, which were assumed to approximate the uniform electron gas, are very well described in terms of short-range dynamical correlations that are only local to an atom.

  4. Discovering new materials with desired properties has been a dominant and crucial topic of interest in the field of materials science in the past few decades. In this work, novel carbon allotropes and ternary B–C–N structures were generated using the state-of-the-art RG 2 code. All structures were fully optimized using density functional theory with first-principles calculations. Several hundred carbon allotropes and ternary B–C–N structures were identified to be superhard materials. The thermodynamic stability of some randomly selected superhard materials was confirmed by evaluating the full phonon dispersions in the Brillouin zone. The new carbon allotropes and ternary B–C–N structures possess a wide range of mechanical properties generally and Vickers hardness specifically. Through 2D Pearson's correlation map, we first reproduced the well-accepted explanation and relationship of the Vickers hardness of the generated structures with other mechanical properties such as shear modulus, bulk modulus, Pugh's ratio, universal anisotropy, and Poisson's ratio. We then propose two fundamentally new descriptors from the electronic level, namely local potential and electron localization function averaged over a unit cell, both of which exhibit a strong correlation with Vickers hardness. More importantly, these descriptors are easy to access from first-principles calculations (at least two orders of magnitude fastermore »than the traditional calculation of elastic constants), and thus can serve as a fast and accurate approach for screening superhard materials. We also combined these new descriptors with known composition and structural descriptors in the machine learning training process. The new descriptors significantly enhance the performance of the trained machine learning model in predicting the Vickers hardness of unknown materials, which provides strong evidence for local potential and electron localization function to be considered in future high-throughput computation. This work unravels more fundamental but previously unexplored knowledge about superhard materials and the newly proposed electronic level descriptors are expected to accelerate the discovery of new superhard materials.« less
  5. Abstract

    Increasing demand for self-powered wearable sensors has spurred an urgent need to develop energy harvesting systems that can reliably and sufficiently power these devices. Within the last decade, reverse electrowetting-on-dielectric (REWOD)-based mechanical motion energy harvesting has been developed, where an electrolyte is modulated (repeatedly squeezed) between two dissimilar electrodes under an externally applied mechanical force to generate an AC current. In this work, we explored various combinations of electrolyte concentrations, dielectrics, and dielectric thicknesses to generate maximum output power employing REWOD energy harvester. With the objective of implementing a fully self-powered wearable sensor, a “zero applied-bias-voltage” approach was adopted. Three different concentrations of sodium chloride aqueous solutions (NaCl-0.1 M, NaCl-0.5 M, and NaCl-1.0 M) were used as electrolytes. Likewise, electrodes were fabricated with three different dielectric thicknesses (100 nm, 150 nm, and 200 nm) of Al2O3and SiO2with an additional layer of CYTOP for surface hydrophobicity. The REWOD energy harvester and its electrode–electrolyte layers were modeled using lumped components that include a resistor, a capacitor, and a current source representing the harvester. Without using any external bias voltage, AC current generation with a power density of 53.3 nW/cm2was demonstrated at an external excitation frequency of 3 Hz with an optimal external load. The experimental results were analytically verifiedmore »using the derived theoretical model. Superior performance of the harvester in terms of the figure-of-merit comparing previously reported works is demonstrated. The novelty of this work lies in the combination of an analytical modeling method and experimental validation that together can be used to increase the REWOD harvested power extensively without requiring any external bias voltage.

    « less