-
Polaronic Enhancement of Second-Harmonic Generation in Lithium Niobate
Authors:
Agnieszka L. Kozub,
Arno Schindlmayr,
Uwe Gerstmann,
Wolf Gero Schmidt
Abstract:
Density-functional theory within a Berry-phase formulation of the dynamical polarization is used to determine the second-order susceptibility $χ^{(2)}$ of lithium niobate (LiNbO$_3$). Defect trapped polarons and bipolarons are found to strongly enhance the nonlinear susceptibility of the material, in particular if localized at Nb$_\mathrm{V}$-V$_{\mathrm{Li}}$ defect pairs. This is essentially a c…
▽ More
Density-functional theory within a Berry-phase formulation of the dynamical polarization is used to determine the second-order susceptibility $χ^{(2)}$ of lithium niobate (LiNbO$_3$). Defect trapped polarons and bipolarons are found to strongly enhance the nonlinear susceptibility of the material, in particular if localized at Nb$_\mathrm{V}$-V$_{\mathrm{Li}}$ defect pairs. This is essentially a consequence of the polaronic excitation resulting in relaxation-induced gap states. The occupation of these levels leads to strongly enhanced $χ^{(2)}$ coefficients and allows for the spatial and transient modification of the second-harmonic generation of macroscopic samples.
△ Less
Submitted 2 June, 2021;
originally announced June 2021.
-
Analytic evaluation of the electronic self-energy in the GW approximation for two electrons on a sphere
Authors:
Arno Schindlmayr
Abstract:
The GW approximation for the electronic self-energy is an important tool for the quantitative prediction of excited states in solids, but its mathematical exploration is hampered by the fact that it must, in general, be evaluated numerically even for very simple systems. In this paper I describe a nontrivial model consisting of two electrons on the surface of a sphere, interacting with the normal…
▽ More
The GW approximation for the electronic self-energy is an important tool for the quantitative prediction of excited states in solids, but its mathematical exploration is hampered by the fact that it must, in general, be evaluated numerically even for very simple systems. In this paper I describe a nontrivial model consisting of two electrons on the surface of a sphere, interacting with the normal long-range Coulomb potential, and show that the GW self-energy, in the absence of self-consistency, can in fact be derived completely analytically in this case. The resulting expression is subsequently used to analyze the convergence of the energy gap between the highest occupied and the lowest unoccupied quasiparticle orbital with respect to the total number of states included in the spectral summations. The asymptotic formula for the truncation error obtained in this way, whose dominant contribution is proportional to the cutoff energy to the power -3/2, may be adapted to extrapolate energy gaps in other systems.
△ Less
Submitted 26 February, 2013;
originally announced February 2013.
-
First-principles calculation of electronic excitations in solids with SPEX
Authors:
Arno Schindlmayr,
Christoph Friedrich,
Ersoy Sasioglu,
Stefan Blügel
Abstract:
We describe the software package SPEX, which allows first-principles calculations of quasiparticle and collective electronic excitations in solids using techniques from many-body perturbation theory. The implementation is based on the full-potential linearized augmented-plane-wave (FLAPW) method, which treats core and valence electrons on an equal footing and can be applied to a wide range of mate…
▽ More
We describe the software package SPEX, which allows first-principles calculations of quasiparticle and collective electronic excitations in solids using techniques from many-body perturbation theory. The implementation is based on the full-potential linearized augmented-plane-wave (FLAPW) method, which treats core and valence electrons on an equal footing and can be applied to a wide range of materials, including transition metals and rare earths. After a discussion of essential features that contribute to the high numerical efficiency of the code, we present illustrative results for quasiparticle band structures calculated within the GW approximation for the electronic self-energy, electron-energy-loss spectra with inter- and intraband transitions as well as local-field effects, and spin-wave spectra of itinerant ferromagnets. In all cases the inclusion of many-body correlation terms leads to very good quantitative agreement with experimental spectroscopies.
△ Less
Submitted 7 October, 2011;
originally announced October 2011.
-
Optical conductivity of metals from first principles
Authors:
Arno Schindlmayr
Abstract:
A computational method to obtain optical conductivities from first principles is presented. It exploits a relation between the conductivity and the complex dielectric function, which is constructed from the full electronic band structure within the random-phase approximation. In contrast to the Drude model, no empirical parameters are used. As interband transitions as well as local-field effects a…
▽ More
A computational method to obtain optical conductivities from first principles is presented. It exploits a relation between the conductivity and the complex dielectric function, which is constructed from the full electronic band structure within the random-phase approximation. In contrast to the Drude model, no empirical parameters are used. As interband transitions as well as local-field effects are properly included, the calculated spectra are valid over a wide frequency range. As an illustration I present quantitative results for selected simple metals, noble metals, and ferromagnetic transition metals. The implementation is based on the full-potential linearized augmented-plane-wave method.
△ Less
Submitted 13 September, 2011;
originally announced September 2011.
-
Efficient implementation of the GW approximation within the all-electron FLAPW method
Authors:
Christoph Friedrich,
Stefan Blügel,
Arno Schindlmayr
Abstract:
We present an implementation of the GW approximation for the electronic self-energy within the full-potential linearized augmented-plane-wave (FLAPW) method. The algorithm uses an all-electron mixed product basis for the representation of response matrices and related quantities. This basis is derived from the FLAPW basis and is exact for wave-function products. The singularity of the bare and scr…
▽ More
We present an implementation of the GW approximation for the electronic self-energy within the full-potential linearized augmented-plane-wave (FLAPW) method. The algorithm uses an all-electron mixed product basis for the representation of response matrices and related quantities. This basis is derived from the FLAPW basis and is exact for wave-function products. The singularity of the bare and screened interaction potentials gives rise to a numerically important self-energy contribution, which we treat analytically to achieve good convergence with respect to the k-point sampling. As numerical realizations of the GW approximation typically suffer from the high computational expense required for the evaluation of the nonlocal and frequency-dependent self-energy, we demonstrate how the algorithm can be made very efficient by exploiting spatial and time-reversal symmetry as well as by applying an optimization of the mixed product basis that retains only the numerically important contributions of the electron-electron interaction. Furthermore, we demonstrate that one can employ an extrapolar approximation for high-lying states to reduce the number of empty states that must be taken into account explicitly in the construction of the polarization function and the self-energy. We show convergence tests, CPU timings, and results for prototype semiconductors and insulators as well as ferromagnetic nickel.
△ Less
Submitted 12 November, 2010; v1 submitted 1 March, 2010;
originally announced March 2010.
-
Wannier-function approach to spin excitations in solids
Authors:
Ersoy Sasioglu,
Arno Schindlmayr,
Christoph Friedrich,
Frank Freimuth,
Stefan Blügel
Abstract:
We present a computational scheme to study spin excitations in magnetic materials from first principles. The central quantity is the transverse spin susceptibility, from which the complete excitation spectrum, including single-particle spin-flip Stoner excitations and collective spin-wave modes, can be obtained. The susceptibility is derived from many-body perturbation theory and includes dynami…
▽ More
We present a computational scheme to study spin excitations in magnetic materials from first principles. The central quantity is the transverse spin susceptibility, from which the complete excitation spectrum, including single-particle spin-flip Stoner excitations and collective spin-wave modes, can be obtained. The susceptibility is derived from many-body perturbation theory and includes dynamic correlation through a summation over ladder diagrams that describe the coupling of electrons and holes with opposite spins. In contrast to earlier studies, we do not use a model potential with adjustable parameters for the electron-hole interaction but employ the random-phase approximation. To reduce the numerical cost for the calculation of the four-point scattering matrix we perform a projection onto maximally localized Wannier functions, which allows us to truncate the matrix efficiently by exploiting the short spatial range of electronic correlation in the partially filled d or f orbitals. Our implementation is based on the FLAPW method. Starting from a ground-state calculation within the LSDA, we first analyze the matrix elements of the screened Coulomb potential in the Wannier basis for the 3d transition-metal series. In particular, we discuss the differences between a constrained nonmagnetic and a proper spin-polarized treatment for the ferromagnets Fe, Co, and Ni. The spectrum of single-particle and collective spin excitations in fcc Ni is then studied in detail. The calculated spin-wave dispersion is in good overall agreement with experimental data and contains both an acoustic and an optical branch for intermediate wave vectors along the [100] direction. In addition, we find evidence for a similar double-peak structure in the spectral function along the [111] direction.
△ Less
Submitted 25 February, 2010;
originally announced February 2010.
-
Efficient calculation of the Coulomb matrix and its expansion around k=0 within the FLAPW method
Authors:
Christoph Friedrich,
Arno Schindlmayr,
Stefan Bluegel
Abstract:
We derive formulas for the Coulomb matrix within the full-potential linearized augmented-plane-wave (FLAPW) method. The Coulomb matrix is a central ingredient in implementations of many-body perturbation theory, such as the Hartree-Fock and GW approximations for the electronic self-energy or the random-phase approximation for the dielectric function. It is represented in the mixed product basis,…
▽ More
We derive formulas for the Coulomb matrix within the full-potential linearized augmented-plane-wave (FLAPW) method. The Coulomb matrix is a central ingredient in implementations of many-body perturbation theory, such as the Hartree-Fock and GW approximations for the electronic self-energy or the random-phase approximation for the dielectric function. It is represented in the mixed product basis, which combines numerical muffin-tin functions and interstitial plane waves that are here expanded with the Rayleigh formula. The resulting algorithm is very efficient in terms of both computational cost and accuracy and is superior to an implementation with the Fourier transform of the step function. In order to allow an analytic treatment of the divergence at k=0 in reciprocal space, we expand the Coulomb matrix analytically around this point without resorting to a projection onto plane waves. We then apply a basis transformation that diagonalizes the Coulomb matrix and confines the divergence to a single eigenvalue. At the same time, response matrices like the dielectric function separate into head, wings, and body with the same mathematical properties as in a plane-wave basis. As an illustration we apply the formulas to electron-energy-loss spectra for nickel at different k vectors including k=0. The convergence of the spectra towards the result at k=0 is clearly seen. Our treatment also allows to include transitions from core states that give rise to a shallow peak at high energies and lead to good agreement with experiment.
△ Less
Submitted 14 November, 2008;
originally announced November 2008.
-
Screening in 2D: GW calculations for surfaces and thin films using the repeated-slab approach
Authors:
Christoph Freysoldt,
Philipp Eggert,
Patrick Rinke,
Arno Schindlmayr,
Matthias Scheffler
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 criticall…
▽ More
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.
△ Less
Submitted 11 June, 2008; v1 submitted 11 January, 2008;
originally announced January 2008.
-
Quasiparticle Corrections to the Electronic Properties of Anion Vacancies at GaAs(110) and InP(110)
Authors:
Magnus Hedström,
Arno Schindlmayr,
Günther Schwarz,
Matthias Scheffler
Abstract:
We propose a new method for calculating optical defect levels and thermodynamic charge-transition levels of point defects in semiconductors, which includes quasiparticle corrections to the Kohn-Sham eigenvalues of density-functional theory. Its applicability is demonstrated for anion vacancies at the (110) surfaces of III-V semiconductors. We find the (+/0) charge-transition level to be 0.49 eV…
▽ More
We propose a new method for calculating optical defect levels and thermodynamic charge-transition levels of point defects in semiconductors, which includes quasiparticle corrections to the Kohn-Sham eigenvalues of density-functional theory. Its applicability is demonstrated for anion vacancies at the (110) surfaces of III-V semiconductors. We find the (+/0) charge-transition level to be 0.49 eV above the surface valence-band maximum for GaAs(110) and 0.82 eV for InP(110). The results show a clear improvement over the local-density approximation and agree closely with an experimental analysis.
△ Less
Submitted 27 November, 2006;
originally announced November 2006.
-
Dielectric anisotropy in the GW space-time method
Authors:
Christoph Freysoldt,
Philipp Eggert,
Patrick Rinke,
Arno Schindlmayr,
R. W. Godby,
Matthias Scheffler
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 C…
▽ More
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.
△ Less
Submitted 9 August, 2006;
originally announced August 2006.
-
Elimination of the linearization error in GW calculations based on the linearized augmented-plane-wave method
Authors:
Christoph Friedrich,
Arno Schindlmayr,
Stefan Blügel,
Takao Kotani
Abstract:
This paper investigates the influence of the basis set on the GW self-energy correction in the full-potential linearized augmented-plane-wave (LAPW) approach and similar linearized all-electron methods. A systematic improvement is achieved by including local orbitals that are defined as second and higher energy derivatives of solutions to the radial scalar-relativistic Dirac equation and thus co…
▽ More
This paper investigates the influence of the basis set on the GW self-energy correction in the full-potential linearized augmented-plane-wave (LAPW) approach and similar linearized all-electron methods. A systematic improvement is achieved by including local orbitals that are defined as second and higher energy derivatives of solutions to the radial scalar-relativistic Dirac equation and thus constitute a natural extension of the LAPW basis set. Within this approach linearization errors can be eliminated, and the basis set becomes complete. While the exchange contribution to the self-energy is little affected by the increased basis-set flexibility, the correlation contribution benefits from the better description of the unoccupied states, as do the quasiparticle energies. The resulting band gaps remain relatively unaffected, however; for Si we find an increase of 0.03 eV.
△ Less
Submitted 23 June, 2006;
originally announced June 2006.
-
Quasiparticle Calculations for Point Defects on Semiconductor Surfaces
Authors:
Magnus Hedstrom,
Arno Schindlmayr,
Matthias Scheffler
Abstract:
We discuss the implementation of quasiparticle calculations for point defects on semiconductor surfaces and, as a specific example, present an ab initio study of the electronic structure of the As vacancy in the +1 charge state on the GaAs(110) surface. The structural properties are calculated with the plane-wave pseudopotential method, and the quasiparticle energies are obtained from Hedin's GW…
▽ More
We discuss the implementation of quasiparticle calculations for point defects on semiconductor surfaces and, as a specific example, present an ab initio study of the electronic structure of the As vacancy in the +1 charge state on the GaAs(110) surface. The structural properties are calculated with the plane-wave pseudopotential method, and the quasiparticle energies are obtained from Hedin's GW approximation. Our calculations show that the 1a" vacancy state in the band gap is shifted from 0.06 to 0.65 eV above the valence-band maximum after the self-energy correction to the Kohn-Sham eigenvalues. The GW result is in close agreement with a recent surface photovoltage imaging measurement.
△ Less
Submitted 30 September, 2002;
originally announced September 2002.
-
Self-consistency and vertex corrections beyond the GW approximation
Authors:
Arno Schindlmayr
Abstract:
The good performance of the GW approximation for band-structure calculations in solids was long taken as a sign that the sum of self-energy diagrams is converged and that all omitted terms are small. However, with modern computational resources it has now become possible to evaluate self-consistency and vertex corrections explicitly, and the numerical results show that they are, in general, not…
▽ More
The good performance of the GW approximation for band-structure calculations in solids was long taken as a sign that the sum of self-energy diagrams is converged and that all omitted terms are small. However, with modern computational resources it has now become possible to evaluate self-consistency and vertex corrections explicitly, and the numerical results show that they are, in general, not individually negligible. In this review the available data is examined, and the implications for practical calculations and the theoretical foundation of the GW approximation are discussed.
△ Less
Submitted 26 June, 2002;
originally announced June 2002.
-
Diagrammatic self-energy approximations and the total particle number
Authors:
Arno Schindlmayr,
P. Garcia-Gonzalez,
R. W. Godby
Abstract:
There is increasing interest in many-body perturbation theory as a practical tool for the calculation of ground-state properties. As a consequence, unambiguous sum rules such as the conservation of particle number under the influence of the Coulomb interaction have acquired an importance that did not exist for calculations of excited-state properties. In this paper we obtain a rigorous, simple r…
▽ More
There is increasing interest in many-body perturbation theory as a practical tool for the calculation of ground-state properties. As a consequence, unambiguous sum rules such as the conservation of particle number under the influence of the Coulomb interaction have acquired an importance that did not exist for calculations of excited-state properties. In this paper we obtain a rigorous, simple relation whose fulfilment guarantees particle-number conservation in a given diagrammatic self-energy approximation. Hedin's G_0W_0 approximation does not satisfy this relation and hence violates the particle-number sum rule. Very precise calculations for the homogeneous electron gas and a model inhomogeneous electron system allow the extent of the nonconservation to be estimated.
△ Less
Submitted 20 October, 2001;
originally announced October 2001.
-
Exchange-correlation kernels for excited states in solids
Authors:
Krzysztof Tatarczyk,
Arno Schindlmayr,
Matthias Scheffler
Abstract:
The performance of several common approximations for the exchange-correlation kernel within time-dependent density-functional theory is tested for elementary excitations in the homogeneous electron gas. Although the adiabatic local-density approximation gives a reasonably good account of the plasmon dispersion, systematic errors are pointed out and traced to the neglect of the wavevector depende…
▽ More
The performance of several common approximations for the exchange-correlation kernel within time-dependent density-functional theory is tested for elementary excitations in the homogeneous electron gas. Although the adiabatic local-density approximation gives a reasonably good account of the plasmon dispersion, systematic errors are pointed out and traced to the neglect of the wavevector dependence. Kernels optimized for atoms are found to perform poorly in extended systems due to an incorrect behavior in the long-wavelength limit, leading to quantitative deviations that significantly exceed the experimental error bars for the plasmon dispersion in the alkali metals.
△ Less
Submitted 16 March, 2001;
originally announced March 2001.
-
Decay properties of the one-particle Green function in real space and imaginary time
Authors:
Arno Schindlmayr
Abstract:
The decay properties of the one-particle Green function in real space and imaginary time are systematically studied for solids. I present an analytic solution for the homogeneous electron gas at finite and at zero temperature as well as asymptotic formulas for real metals and insulators that allow an analytic treatment in electronic-structure calculations based on a space-time representation. Th…
▽ More
The decay properties of the one-particle Green function in real space and imaginary time are systematically studied for solids. I present an analytic solution for the homogeneous electron gas at finite and at zero temperature as well as asymptotic formulas for real metals and insulators that allow an analytic treatment in electronic-structure calculations based on a space-time representation. The generic dependence of the decay constants on known system parameters is used to compare the scaling of reciprocal-space algorithms for the GW approximation and the space-time method.
△ Less
Submitted 28 August, 2000;
originally announced August 2000.
-
Universality of the Hohenberg-Kohn functional
Authors:
Arno Schindlmayr
Abstract:
We comment on a recent article by H. L. Neal [Am. J. Phys. 66, 512 (1998)], in which an analytic expression for the Hohenberg-Kohn functional was derived for one-dimensional two-particle systems with the harmonic interaction. We argue that the functional given by Neal is not exact but really constitutes an approximation in the same spirit as the local-density approximation.
We comment on a recent article by H. L. Neal [Am. J. Phys. 66, 512 (1998)], in which an analytic expression for the Hohenberg-Kohn functional was derived for one-dimensional two-particle systems with the harmonic interaction. We argue that the functional given by Neal is not exact but really constitutes an approximation in the same spirit as the local-density approximation.
△ Less
Submitted 12 March, 1999;
originally announced March 1999.
-
Spectra and total energies from self-consistent many-body perturbation theory
Authors:
Arno Schindlmayr,
Thomas J. Pollehn,
R. W. Godby
Abstract:
With the aim of identifying universal trends, we compare fully self-consistent electronic spectra and total energies obtained from the GW approximation with those from an extended GWGamma scheme that includes a nontrivial vertex function and the fundamentally distinct Bethe-Goldstone approach based on the T-matrix. The self-consistent Green's function G, as derived from Dyson's equation, is used…
▽ More
With the aim of identifying universal trends, we compare fully self-consistent electronic spectra and total energies obtained from the GW approximation with those from an extended GWGamma scheme that includes a nontrivial vertex function and the fundamentally distinct Bethe-Goldstone approach based on the T-matrix. The self-consistent Green's function G, as derived from Dyson's equation, is used not only in the self-energy but also to construct the screened interaction W for a model system. For all approximations we observe a similar deterioration of the spectrum, which is not removed by vertex corrections. In particular, satellite peaks are systematically broadened and move closer to the chemical potential. The corresponding total energies are universally raised, independent of the system parameters. Our results therefore suggest that any improvement in total energy due to self-consistency, such as for the electron gas in the GW approximation, may be fortuitous.
△ Less
Submitted 9 June, 1998;
originally announced June 1998.
-
Assessment of the GW approximation using Hubbard chains
Authors:
Thomas J. Pollehn,
Arno Schindlmayr,
R. W. Godby
Abstract:
We investigate the performance of the GW approximation by comparison to exact results for small model systems. The role of the chemical potentials in Dyson's equation as well as the consequences of numerical resonance broadening are examined, and we show how a proper treatment can improve computational implementations of many-body perturbation theory in general. GW and exchange-only calculations…
▽ More
We investigate the performance of the GW approximation by comparison to exact results for small model systems. The role of the chemical potentials in Dyson's equation as well as the consequences of numerical resonance broadening are examined, and we show how a proper treatment can improve computational implementations of many-body perturbation theory in general. GW and exchange-only calculations are performed over a wide range of fractional band fillings and correlation strengths. We thus identify the physical situations where these schemes are applicable.
△ Less
Submitted 13 November, 1997;
originally announced November 1997.
-
Systematic vertex corrections through iterative solution of Hedin's equations beyond the GW approximation
Authors:
Arno Schindlmayr,
R. W. Godby
Abstract:
We present a general procedure for obtaining progressively more accurate functional expressions for the electron self-energy by iterative solution of Hedin's coupled equations. The iterative process starting from Hartree theory, which gives rise to the GW approximation, is continued further, and an explicit formula for the vertex function from the second full cycle is given. Calculated excitatio…
▽ More
We present a general procedure for obtaining progressively more accurate functional expressions for the electron self-energy by iterative solution of Hedin's coupled equations. The iterative process starting from Hartree theory, which gives rise to the GW approximation, is continued further, and an explicit formula for the vertex function from the second full cycle is given. Calculated excitation energies for a Hubbard Hamiltonian demonstrate the convergence of the iterative process and provide further strong justification for the GW approximation.
△ Less
Submitted 28 October, 1997;
originally announced October 1997.
-
Excitons with anisotropic effective mass
Authors:
Arno Schindlmayr
Abstract:
We present a simple analytic scheme for calculating the binding energy of excitons in semiconductors that takes full account of the existing anisotropy in the effective mass, as a complement to the qualitative treatment in most textbooks. Results obtained for excitons in gallium nitride form the basis for a discussion of the accuracy of this approach.
We present a simple analytic scheme for calculating the binding energy of excitons in semiconductors that takes full account of the existing anisotropy in the effective mass, as a complement to the qualitative treatment in most textbooks. Results obtained for excitons in gallium nitride form the basis for a discussion of the accuracy of this approach.
△ Less
Submitted 25 September, 1997;
originally announced September 1997.
-
Violation of particle number conservation in the GW approximation
Authors:
Arno Schindlmayr
Abstract:
We present a nontrivial model system of interacting electrons that can be solved analytically in the GW approximation. We obtain the particle number from the GW Green's function strictly analytically, and prove that there is a genuine violation of particle number conservation if the self-energy is calculated non-self-consistently from a zeroth order Green's function, as done in virtually all pra…
▽ More
We present a nontrivial model system of interacting electrons that can be solved analytically in the GW approximation. We obtain the particle number from the GW Green's function strictly analytically, and prove that there is a genuine violation of particle number conservation if the self-energy is calculated non-self-consistently from a zeroth order Green's function, as done in virtually all practical implementations. We also show that a simple shift of the self-energy that partially restores self-consistency reduces the numerical deviation significantly.
△ Less
Submitted 25 September, 1997;
originally announced September 1997.
-
Density-functional theory and the v-representability problem for model strongly correlated electron systems
Authors:
Arno Schindlmayr,
R. W. Godby
Abstract:
Inspired by earlier work on the band-gap problem in insulators, we reexamine the treatment of strongly correlated Hubbard-type models within density-functional theory. In contrast to previous studies, the density is fully parametrized by occupation numbers and overlap of orbitals centered at neighboring atomic sites, as is the local potential by the hopping matrix. This corresponds to a good for…
▽ More
Inspired by earlier work on the band-gap problem in insulators, we reexamine the treatment of strongly correlated Hubbard-type models within density-functional theory. In contrast to previous studies, the density is fully parametrized by occupation numbers and overlap of orbitals centered at neighboring atomic sites, as is the local potential by the hopping matrix. This corresponds to a good formal agreement between density-functional theory in real space and second quantization. It is shown that density-functional theory is formally applicable to such systems and the theoretical framework is provided. The question of noninteracting v representability is studied numerically for finite one-dimensional clusters, for which exact results are available, and qualitatively for infinite systems. This leads to the conclusion that the electron density corresponding to interacting systems of the type studied here is in fact not noninteracting v representable because the Kohn-Sham electrons are unable to reproduce the correlation-induced localization correctly.
△ Less
Submitted 24 September, 1997;
originally announced September 1997.