Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userDielectric anisotropy in the GW space-time method
73
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
Excited-state calculations, notably for quasiparticle band structures, are nowadays routinely performed within the GW approximation for the electronic self-energy. Nevertheless, certain numerical approximations and simplifications are still employed in practice to make the computations feasible. An important aspect for periodic systems is the proper treatment of the singularity of the screened Coulomb interaction in reciprocal space, which results from the slow 1/r decay in real space. This must be done without introducing artificial interactions between the quasiparticles and their periodic images in repeated cells, which occur when integrals of the screened Coulomb interaction are discretised in reciprocal space. An adequate treatment of both aspects is crucial for a numerically stable computation of the self-energy. In this article we build on existing schemes for isotropic screening and present an extension for anisotropic systems. We also show how the contributions to the dielectric function arising from the non-local part of the pseudopotentials can be computed efficiently. These improvements are crucial for obtaining a fast convergence with respect to the number of points used for the Brillouin zone integration and prove to be essential to make GW calculations for strongly anisotropic systems, such as slabs or multilayers, efficient.
rate research
Read More
We implement the GW space-time method at finite temperatures, in which the Greens function G and the screened Coulomb interaction W are represented in the real space on a suitable mesh and in imaginary time in terms of Chebyshev polynomials, paying particular attention to controlling systematic errors of the representation. Having validated the technique by the canonical application to silicon and germanium, we apply it to calculation of band gaps in hexagonal solid hydrogen with the bare Greens function obtained from density functional approximation and the interaction screened within the random phase approximation (RPA). The results, obtained from the asymptotic decay of the full Greens function without resorting to analytic continuation, suggest that the solid hydrogen above 250 GPa can not adopt the hexagonal-closed-pack (hcp) structure. The demonstrated ability of the method to store the full G and W functions in memory with sufficient accuracy is crucial for its subsequent extensions to include higher orders of the diagrammatic series by means of diagrammatic Monte Carlo algorithms.
Molecule-metal interfaces have a broad range of applications in nanoscale materials science. Accurate characterization of their electronic structures from first-principles is key in understanding material and device properties. The GW approach within many-body perturbation theory is state-of-the-art and can in principle yield accurate quasiparticle energy levels and interfacial level alignments that are in quantitative agreement with experiments. However, the interfaces are large heterogeneous systems that are currently challenging for first-principles GW calculations. In this work, we develop a GW-based dielectric embedding approach for molecule-metal interfaces, significantly reducing the computational cost of direct GW without sacrificing accuracy. To be specific, we perform explicit GW calculations only in the simulation cell of the molecular adsorbate, in which the dielectric effect of the metallic substrate is embedded. This is made possible via a real-space truncation of the substrate polarizability and the use of the interface plasma frequency in the adsorbate GW calculation. Here, we focus on the level alignment at weakly coupled molecule-metal interfaces, i.e., the energy difference between a molecular frontier orbital resonance and the substrate Fermi level. We demonstrate our method and assess a few GW-based approximations using two well-studied systems, benzene adsorbed on the Al (111) and on the graphite (0001) surfaces.
We apply a recently developed quasiparticle self-consistent $GW$ method (QSGW) to Gd, Er, EuN, GdN, ErAs, YbN and GdAs. We show that QSGW combines advantages separately found in conventional $GW$ and LDA+$U$ theory, in a simple and fully emph{ab initio} way. qsgw reproduces the experimental occupied $4f$ levels well, though unoccupied levels are systematically overestimated. Properties of the Fermi surface responsible for electronic properties are in good agreement with available experimental data. GdN is predicted to be very near a critical point of a first-order metal-insulator transition.
We study the spectral function of the homogeneous electron gas using many-body perturbation theory and the cumulant expansion. We compute the angle-resolved spectral function based on the GW approximation and the `GW plus cumulant approach. In agreement with previous studies, the GW spectral function exhibits a spurious plasmaron peak at energies 1.5$omega_{rm pl}$ below the quasiparticle peak, $omega_{rm pl}$ being the plasma energy. The GW plus cumulant approach, on the other hand, reduces significantly the intensity of the plasmon-induced spectral features and renormalizes their energy relative to the quasiparticle energy to $omega_{rm pl}$. Consistently with previous work on semiconductors, our results show that the HEG is characterized by the emergence of plasmonic polaron bands, that is, broadened replica of the quasiparticle bands, red-shifted by the plasmon energy.
We present a high sensitivity method allowing the measurement of the non linear dielectric susceptibility of an insulating material at finite frequency. It has been developped for the study of dynamic heterogeneities in supercooled liquids using dielectric spectroscopy at frequencies 0.05 Hz < f < 30000 Hz . It relies on the measurement of the third harmonics component of the current flowing out of a capacitor. We first show that standard laboratory electronics (amplifiers and voltage sources) nonlinearities lead to limits on the third harmonics measurements that preclude reaching the level needed by our physical goal, a ratio of the third harmonics to the fundamental signal about 7 orders of magnitude lower than 1. We show that reaching such a sensitivity needs a method able to get rid of the nonlinear contributions both of the measuring device (lock-in amplifier) and of the excitation voltage source. A bridge using two sources fulfills only the first of these two requirements, but allows to measure the nonlinearities of the sources. Our final method is based on a bridge with two plane capacitors characterized by different dielectric layer thicknesses. It gets rid of the source and amplifier nonlinearities because in spite of a strong frequency dependence of the capacitors impedance, it is equilibrated at any frequency. We present the first measurements of the physical nonlinear response using our method. Two extensions of the method are suggested.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
