Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userScreening in 2D: GW calculations for surfaces and thin films using the repeated-slab approach
397
0
0.0
(
0
)
Added by
Christoph Freysoldt
Publication date
2008
fields
Physics
and research's language is
English
Ask ChatGPT about the research
No Arabic abstract
In the context of photoelectron spectroscopy, the $GW$ approach has developed into the method of choice for computing excitation spectra of weakly correlated bulk systems and their surfaces. To employ the established computational schemes that have been developed for three-dimensional crystals, two-dimensional systems are typically treated in the repeated-slab approach. In this work we critically examine this approach and identify three important aspects for which the treatment of long-range screening in two dimensions differs from the bulk: (1) anisotropy of the macroscopic screening (2) $mathbf k$-point sampling parallel to the surface (3) periodic repetition and slab-slab interaction. For prototypical semiconductor (silicon) and ionic (NaCl) thin films we quantify the individual contributions of points (1) to (3) and develop robust and efficient correction schemes derived from the classic theory of dielectric screening.
rate research
Read More
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. Quasiparticle $GW$ calculations 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 dense $k$-grid to properly converge the calculated quasiparticle energies. In this work, we present a combined non-uniform sub-sampling and analytical integration method that can drastically improve the efficiency of the BZ integration in 2D $GW$ calculations. 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 $Sigma(mathbf{q})$ in the small $mathbf{q}$ limit. This method, when combined with another accelerated $GW$ method that we developed recently, can drastically speed-up (by over three orders of magnitude) $GW$ calculations for 2D materials. Our method allows fully converged $GW$ calculations 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 the capability and performance of our new method, we have carried out fully converged $GW$ calculations for monolayer C$_2$N, a recently discovered 2D material with a large unit cell, and investigate its quasiparticle band structure in detail.
The search for new materials, based on computational screening, relies on methods that accurately predict, in an automatic manner, total energy, atomic-scale geometries, and other fundamental characteristics of materials. Many technologically important material properties directly stem from the electronic structure of a material, but the usual workhorse for total energies, namely density-functional theory, is plagued by fundamental shortcomings and errors from approximate exchange-correlation functionals in its prediction of the electronic structure. At variance, the $GW$ method is currently the state-of-the-art {em ab initio} approach for accurate electronic structure. It is mostly used to perturbatively correct density-functional theory results, but is however computationally demanding and also requires expert knowledge to give accurate results. Accordingly, it is not presently used in high-throughput screening: fully automatized algorithms for setting up the calculations and determining convergence are lacking. In this work we develop such a method and, as a first application, use it to validate the accuracy of $G_0W_0$ using the PBE starting point, and the Godby-Needs plasmon pole model ($G_0W_0^textrm{GN}$@PBE), on a set of about 80 solids. The results of the automatic convergence study utilized provides valuable insights. Indeed, we find correlations between computational parameters that can be used to further improve the automatization of $GW$ calculations. Moreover, we find that $G_0W_0^textrm{GN}$@PBE shows a correlation between the PBE and the $G_0W_0^textrm{GN}$@PBE gaps that is much stronger than that between $GW$ and experimental gaps. However, the $G_0W_0^textrm{GN}$@PBE gaps still describe the experimental gaps more accurately than a linear model based on the PBE gaps.
By means of quasiparticle-energy calculations in the G0W0 approach, we show for the prototypical insulator/semiconductor system NaCl/Ge(001) that polarization effects at the interfaces noticeably affect the excitation spectrum of molecules adsorbed on the surface of the NaCl films. The magnitude of the effect can be controlled by varying the thickness of the film, offering new opportunities for tuning electronic excitations in e.g. molecular electronics or quantum transport. Polarization effects are visible even for the excitation spectrum of the NaCl films themselves, which has important implications for the interpretation of surface science experiments for the characterization of insulator surfaces.
In view of their immensely intriguing properties, two dimensional materials are being intensely researched in search of novel phenomena and diverse application interests, however, studies on the realization of nanocomposites in the application-worthy thin-film platform are rare. Here we have grown MoS2-hBN composite thin films on different substrates by the pulsed laser deposition technique and made comparative studies with the pristine MoS2 and hBN films. The Raman, XPS and HRTEM confirm the concomitant presence of both the 1T (conducting) and 2H (semiconducting) polymorphs of MoS2 in the composite film. Interestingly, a peculiar reentrant semiconductor-metal-insulator transition is seen in the composite film which is absent in the MoS2 film, and it correlates well with the signatures of phonon softening seen in temperature-dependent Raman spectroscopy. Furthermore, electrostatic force microscopy reveals the presence of three distinct regions (metallic, semiconducting, and insulating) in the composite film with differing contact potentials and enhanced propensity for charge transfer with respect to pristine MoS2. A triboelectric nanogenerator device containing biphasic composite film as an electron acceptor exhibits more than twofold (sixfold) enhancement in peak-to-peak output voltage as compared to the pristine MoS2 (hBN) film. These observations bring out the potential of nanocomposite thin films for unfolding emergent phenomena and technological applications.
Segregation at surfaces of metal-covalent binary liquids is often non-classical and in extreme cases such as AuSi, the surface crystallizes above the melting point. In this study, we employ atomic-scale computational frameworks to study the surface crystallization of AuSi films and droplets as a function of composition, temperature and size. For temperatures in the range $T_s^ast=765-780$K above the melting point $(T_s^astapprox1.3,T_m)$, both thin film and droplet surfaces undergo a first order transition, from a 2D Au$_2$Si crystalline phase to a laterally disordered yet stratified layer. The thin film surfaces exhibit an effective surface tension that increases with temperature and decreases with Si concentration. On the other hand, for droplets in the size range $10-30$ nm, the bulk Laplace pressure alters the surface segregation as it occurs with respect to a strained bulk. Above $T_s^ast$ the size effect on the surface tension is small, while for $T<T_s^ast$ the surface layer is strained and composed of 2D crystallites separated by extended grain boundary scars that lead to large fluctuations in its energetics. As a specific application, all-atom simulations of AuSi droplets on Si(111) substrate subject to Si surface flux show that the supersaturation dependent surface tension destabilizes the contact line via formation of a precursor wetting film on the solid-vapor interface, and has ramifications for size selection during VLS-based routes for nanowire growth. Our study sheds light on the interplay between stability and energetics of surfaces in these unique class of binary alloys and offers pathways for exploiting their surface structure for varied applications such as catalytic nanocrystal growth, dealloying, and polymer crystallization.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
