Precise response functions in all-electron methods: Application to the optimized-effective-potential approach
arXiv:1205.6335 · doi:10.1103/PhysRevB.85.245124
Abstract
The optimized-effective-potential (OEP) method is a special technique to construct local Kohn-Sham potentials from general orbital-dependent energy functionals. In a recent publication [M. Betzinger, C. Friedrich, S. Blügel, A. Görling, Phys. Rev. B 83, 045105 (2011)] we showed that uneconomically large basis sets were required to obtain a smooth local potential without spurious oscillations within the full-potential linearized augmented-plane-wave method (FLAPW). This could be attributed to the slow convergence behavior of the density response function. In this paper, we derive an incomplete-basis-set correction for the response, which consists of two terms: (1) a correction that is formally similar to the Pulay correction in atomic-force calculations and (2) a numerically more important basis response term originating from the potential dependence of the basis functions. The basis response term is constructed from the solutions of radial Sternheimer equations in the muffin-tin spheres. With these corrections the local potential converges at much smaller basis sets, at much fewer states, and its construction becomes numerically very stable. We analyze the improvements for rock-salt ScN and report results for BN, AlN, and GaN, as well as the perovskites CaTiO3, SrTiO3, and BaTiO3. The incomplete-basis-set correction can be applied to other electronic-structure methods with potential-dependent basis sets and opens the perspective to investigate a broad spectrum of problems in theoretical solid-state physics that involve response functions.
17 pages, 7 figures, 3 tables
References in corpus (11)
- The SIESTA method for ab initio order-N materials simulation
- Efficient index handling of multidimensional periodic boundary conditions
- Linear-response theory and lattice dynamics: a muffin-tin orbital approach
- Efficient implementation of the GW approximation within the all-electron FLAPW method
- Band convergence and linearization error correction of all-electron GW calculations: The extreme case of zinc oxide
- Exact-exchange based quasiparticle energy calculations for the band gap, effective masses and deformation potentials of ScN
- Linear scaling Krylov subspace method for large scale {\it ab initio} electronic structure calculations of metals
- Elimination of the linearization error in GW calculations based on the linearized augmented-plane-wave method
- Hybrid functionals within the all-electron FLAPW method: implementation and applications of PBE0
- Local exact exchange potentials within the all-electron FLAPW method and a comparison with pseudopotential results
- Efficient calculation of the Coulomb matrix and its expansion around k=0 within the FLAPW method
Cited by in corpus (14)
- The GW compendium: A practical guide to theoretical photoemission spectroscopy
- A self-interaction-free local hybrid functional: Accurate binding energies vis-à-vis accurate ionization potentials from Kohn-Sham eigenvalues
- Semilocal exchange-correlation potentials for solid-state calculations: Current status and future directions
- All-electron periodic implementation with numerical atomic orbital basis functions: algorithm and benchmarks
- Kohn-Sham band gaps and potentials of solids from the optimised effective potential method within the random phase approximation
- Comparison between exact and semilocal exchange potentials: An all-electron study for solids
- One-body reduced density-matrix functional theory in finite basis sets at elevated temperatures
- Electronic Band Structure of Cuprous and Silver Halides: a Numerically Accurate All-Electron Study
- Electronic structure of TiSe from a quasi-self-consistent GW approach
- Approximations to the exact exchange potential: KLI versus semilocal
- Ground state properties of 3d metals from self-consistent GW approach
- Phonons from Density-Functional Perturbation Theory using the All-Electron Full-Potential Linearized Augmented Plane-Wave Method FLEUR
- Local Exchange-Correlation Potentials by Density Inversion in Solids
- Self-consistent random phase approximation and optimized hybrid functionals for solids