Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userEfficient dielectric matrix calculations using the Lanczos algorithm for fast many-body $G_0W_0$ implementations
124
0
0.0
(
0
)
Added by
Jonathan Laflamme Janssen
Publication date
2014
fields
Physics
and research's language is
English
Ask ChatGPT about the research
No Arabic abstract
We present a $G_0W_0$ implementation that assesses the two major bottlenecks of traditional plane-waves implementations, the summations over conduction states and the inversion of the dielectric matrix, without introducing new approximations in the formalism. The first bottleneck is circumvented by converting the summations into Sternheimer equations. Then, the novel avenue of expressing the dielectric matrix in a Lanczos basis is developed, which reduces the matrix size by orders of magnitude while being computationally efficient. We also develop a model dielectric operator that allows us to further reduce the size of the dielectric matrix without accuracy loss. Furthermore, we develop a scheme that reduces the numerical cost of the contour deformation technique to the level of the lightest plasmon pole model. Finally, the use of the simplified quasi-minimal residual scheme in replacement of the conjugate gradients algorithm allows a direct evaluation of the $G_0W_0$ corrections at the desired real frequencies, without need for analytical continuation. The performance of the resulting $G_0W_0$ implementation is demonstrated by comparison with a traditional plane-waves implementation, which reveals a 500-fold speedup for the silane molecule. Finally, the accuracy of our $G_0W_0$ implementation is demonstrated by comparison with other $G_0W_0$ calculations and experimental results.
rate research
Read More
Ab initio many-body perturbation theory within the $GW$ approximation is a Greens function formalism widely used in the calculation of quasiparticle excitation energies of solids. In what has become an increasingly standard approach, Kohn-Sham eigenenergies, generated from a DFT calculation with a strategically-chosen exchange correlation functional ``starting point, are used to construct $G$ and $W$, and then perturbatively corrected by the resultant $GW$ self-energy. In practice, there are several ways to construct the $GW$ self-energy, and these can lead to variations in predicted quasiparticle energies. For example, for ZnO and TiO$_2$, reported $GW$ fundamental gaps can vary by more than 1 eV. In this work, we address the convergence and key approximations in contemporary $G_0W_0$ calculations, including frequency-integration schemes and the treatment of the Coulomb divergence in the exact-exchange term. We study several systems,and compare three different $GW$ codes: BerkeleyGW, Abinit and Yambo. We demonstrate, for the first time, that the same quasiparticle energies for systems in the condensed phase can be obtained with different codes, and we provide a comprehensive assessment of implementations of the $GW$ approximation.
It is often computationally advantageous to model space as a discrete set of points forming a lattice grid. This technique is particularly useful for computationally difficult problems such as quantum many-body systems. For reasons of simplicity and familiarity, nearly all quantum many-body calculations have been performed on simple cubic lattices. Since the removal of lattice artifacts is often an important concern, it would be useful to perform calculations using more than one lattice geometry. In this work we show how to perform quantum many-body calculations using auxiliary-field Monte Carlo simulations on a three-dimensional body-centered cubic (BCC) lattice. As a benchmark test we compute the ground state energy of 33 spin-up and 33 spin-down fermions in the unitary limit, which is an idealized limit where the interaction range is zero and scattering length is infinite. As a fraction of the free Fermi gas energy $E_{rm FG}$, we find that the ground state energy is $E_0/E_{rm FG}= 0.369(2), 0.371(2),$ using two different definitions of the finite-system energy ratio. This is in excellent agreement with recent results obtained on a cubic lattice cite{He:2019ipt}. We find that the computational effort and performance on a BCC lattice is approximately the same as that for a cubic lattice with the same number of lattice points. We discuss how the lattice simulations with different geometries can be used to constrain the size lattice artifacts in simulations of continuum quantum many-body systems.
The formation energies and electronic structure of europium doped zinc oxide has been determined using DFT and many-body $GW$ methods. In the absence of intrisic defects we find that the europium-$f$ states are located in the ZnO band gap with europium possessing a formal charge of 2+. On the other hand, the presence of intrinsic defects in ZnO allows intraband $f-f$ transitions otherwise forbidden in atomic europium. This result coorroborates with recently observed photoluminescence in the visible red region [1].
The interaction of electrons with crystal lattice vibrations (phonons) and collective charge-density fluctuations (plasmons) influences profoundly the spectral properties of solids revealed by photoemission spectroscopy experiments. Photoemission satellites, for instance, are a prototypical example of quantum emergent behavior that may result from the strong coupling of electronic states to plasmons and phonons. The existence of these spectral features has been verified over energy scales spanning several orders of magnitude (from 50 meV to 15-20 eV) and for a broad class of compounds such as simple metals, semiconductors, and highly-doped oxides. During the past few years the cumulant expansion approach, alongside with the GW approximation and the theory of electron-phonon and electron-plasmon coupling in solids, has evolved into a predictive and quantitatively accurate approach for the description of the spectral signatures of electron-boson coupling entirely from first principles, and it has thus become the state-of-the-art theoretical tool for the description of these phenomena. In this chapter we introduce the fundamental concepts needed to interpret plasmon and phonon satellites in photoelectron spectra, and we review recent progress on first-principles calculations of these features using the cumulant expansion method.
Matrix Product States form the basis of powerful simulation methods for ground state problems in one dimension. Their power stems from the fact that they faithfully approximate states with a low amount of entanglement, the area law. In this work, we establish the mixed state analogue of this result: We show that one-dimensional mixed states with a low amount of entanglement, quantified by the entanglement of purification, can be efficiently approximated by Matrix Product Density Operators (MPDOs). In combination with results establishing area laws for thermal states, this helps to put the use of MPDOs in the simulation of thermal states on a formal footing.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
