Menu
Home
People
Places
Arts
History
Plants & Animals
Science
Life & Culture
Technology
Reference.org
Local-density approximation
open-in-new
<h2 id="applications">Applications</h2> <p>Local density approximations, as with GGAs are employed extensively by <a href="/facts/Solid-state_physics/Cmc6hvyR">solid state physicists</a> in ab-initio DFT studies to interpret electronic and magnetic interactions in semiconductor materials including semiconducting oxides and <a href="/facts/Spintronics/dAWl8rcI">spintronics</a>. The importance of these computational studies stems from the system complexities which bring about high sensitivity to synthesis parameters necessitating first-principles based analysis. The prediction of <a href="/facts/Fermi_level/WpyRA0V7">Fermi level</a> and band structure in doped semiconducting oxides is often carried out using LDA incorporated into simulation packages such as CASTEP and DMol3.<a class="footnote-ref" id="fnref:2" href="#fn:2"><sup>2</sup></a> However an underestimation in <a href="/facts/Band_gap/IjIeh95l">Band gap</a> values often associated with LDA and <a href="/facts/Density_functional_theory/wwz4Gqdg">GGA</a> approximations may lead to false predictions of impurity mediated conductivity and/or carrier mediated magnetism in such systems.<a class="footnote-ref" id="fnref:3" href="#fn:3"><sup>3</sup></a> Starting in 1998, the application of the <a href="/facts/Rayleigh_theorem_for_eigenvalues/9hL0ecpS">Rayleigh theorem for eigenvalues</a> has led to mostly accurate, calculated band gaps of materials, using LDA potentials.<a class="footnote-ref" id="fnref:4" href="#fn:4"><sup>4</sup></a><a class="footnote-ref" id="fnref:5" href="#fn:5"><sup>5</sup></a> A misunderstanding of the second theorem of DFT appears to explain most of the underestimation of band gap by LDA and GGA calculations, as explained in the description of <a href="/facts/Density_functional_theory/wwz4Gqdg">density functional theory</a>, in connection with the statements of the two theorems of DFT. </p> <h2 id="homogeneous-electron-gas">Homogeneous electron gas</h2> <p>Approximation for <i>є</i>xc depending only upon the density can be developed in numerous ways. The most successful approach is based on the homogeneous electron gas. This is constructed by placing <i>N</i> interacting electrons in to a volume, <i>V</i>, with a positive background charge keeping the system neutral. <i>N</i> and <i>V</i> are then taken to infinity in the manner that keeps the density (<i>ρ</i> = <i>N</i> / <i>V</i>) finite. This is a useful approximation, as the total energy consists of contributions only from the kinetic energy, electrostatic interaction energy and exchange-correlation energy, and that the wavefunction is expressible in terms of plane waves. In particular, for a constant density <i>ρ</i>, the exchange energy density is proportional to <i>ρ</i>⅓. </p> <h2 id="exchange-functional">Exchange functional</h2> <p>The exchange-energy density of a HEG is known analytically. The LDA for exchange employs this expression under the approximation that the exchange-energy in a system where the density is not homogeneous, is obtained by applying the HEG results pointwise, yielding the expression<a class="footnote-ref" id="fnref:6" href="#fn:6"><sup>6</sup></a><a class="footnote-ref" id="fnref:7" href="#fn:7"><sup>7</sup></a> </p> E x L D A [ ρ ] = − 3 e 2 16 π ε 0 ( 3 π ) 1 / 3 ∫ ρ ( r ) 4 / 3 d r = − 3 4 ( 3 π ) 1 / 3 ∫ ρ ( r ) 4 / 3 d r , {\displaystyle E_{\rm {x}}^{\mathrm {LDA} }[\rho ]=-{\frac {3e^{2}}{16\pi \varepsilon _{0}}}\left({\frac {3}{\pi }}\right)^{1/3}\int \rho (\mathbf {r} )^{4/3}\ \mathrm {d} \mathbf {r} =-{\frac {3}{4}}\left({\frac {3}{\pi }}\right)^{1/3}\int \rho (\mathbf {r} )^{4/3}\ \mathrm {d} \mathbf {r} \,,} <p>where the second formulation applies in <a href="/facts/Atomic_units/sHxrdE9Q">atomic units</a>. </p> <h2 id="correlation-functional">Correlation functional</h2> <p>Analytic expressions for the correlation energy of the HEG are available in the high- and low-density limits corresponding to infinitely-weak and infinitely-strong correlation. For a HEG with density <i>ρ</i>, the high-density limit of the correlation energy density is<a class="footnote-ref" id="fnref:8" href="#fn:8"><sup>8</sup></a> </p> ϵ c = A ln ( r s ) + B + r s ( C ln ( r s ) + D ) , {\displaystyle \epsilon _{\rm {c}}=A\ln(r_{\rm {s}})+B+r_{\rm {s}}(C\ln(r_{\rm {s}})+D)\ ,} <p>and the low limit </p> ϵ c = 1 2 ( g 0 r s + g 1 r s 3 / 2 + … ) , {\displaystyle \epsilon _{\rm {c}}={\frac {1}{2}}\left({\frac {g_{0}}{r_{\rm {s}}}}+{\frac {g_{1}}{r_{\rm {s}}^{3/2}}}+\dots \right)\ ,} <p>where the <a href="/facts/Wigner%25E2%2580%2593Seitz_cell/qU6azUs6">Wigner-Seitz parameter</a> r s {\displaystyle r_{\rm {s}}} is dimensionless.<a class="footnote-ref" id="fnref:9" href="#fn:9"><sup>9</sup></a> It is defined as the radius of a sphere which encompasses exactly one electron, divided by the Bohr radius <i>a</i>0. In terms of the density <i>ρ</i>, this means </p> 4 3 π r s 3 = 1 ρ a 0 3 . {\displaystyle {\frac {4}{3}}\pi r_{\rm {s}}^{3}={\frac {1}{\rho \,a_{0}^{3}}}\ .} <p>An analytical expression for the full range of densities has been proposed based on the many-body perturbation theory. The calculated correlation energies are in agreement with the results from <a href="/facts/Quantum_Monte_Carlo/qenhfw9u">quantum Monte Carlo</a> simulation to within 2 milli-Hartree. </p><p>Accurate <a href="/facts/Quantum_Monte_Carlo/qenhfw9u">quantum Monte Carlo</a> simulations for the energy of the HEG have been performed for several intermediate values of the density, in turn providing accurate values of the correlation energy density.<a class="footnote-ref" id="fnref:10" href="#fn:10"><sup>10</sup></a> </p> <h2 id="spin-polarization">Spin polarization</h2> <p>The extension of density functionals to <a href="/facts/Spin_polarization/S36kxcj2">spin-polarized</a> systems is straightforward for exchange, where the exact spin-scaling is known, but for correlation further approximations must be employed. A spin polarized system in DFT employs two spin-densities, <i>ρ</i>α and <i>ρ</i>β with <i>ρ</i> = <i>ρ</i>α + <i>ρ</i>β, and the form of the local-spin-density approximation (LSDA) is </p> E x c L S D A [ ρ α , ρ β ] = ∫ d r ρ ( r ) ϵ x c ( ρ α , ρ β ) . {\displaystyle E_{\rm {xc}}^{\mathrm {LSDA} }[\rho _{\alpha },\rho _{\beta }]=\int \mathrm {d} \mathbf {r} \ \rho (\mathbf {r} )\epsilon _{\rm {xc}}(\rho _{\alpha },\rho _{\beta })\ .} <p>For the exchange energy, the exact result (not just for local density approximations) is known in terms of the spin-unpolarized functional:<a class="footnote-ref" id="fnref:11" href="#fn:11"><sup>11</sup></a> </p> E x [ ρ α , ρ β ] = 1 2 ( E x [ 2 ρ α ] + E x [ 2 ρ β ] ) . {\displaystyle E_{\rm {x}}[\rho _{\alpha },\rho _{\beta }]={\frac {1}{2}}{\bigg (}E_{\rm {x}}[2\rho _{\alpha }]+E_{\rm {x}}[2\rho _{\beta }]{\bigg )}\ .} <p>The spin-dependence of the correlation energy density is approached by introducing the relative spin-polarization: </p> ζ ( r ) = ρ α ( r ) − ρ β ( r ) ρ α ( r ) + ρ β ( r ) . {\displaystyle \zeta (\mathbf {r} )={\frac {\rho _{\alpha }(\mathbf {r} )-\rho _{\beta }(\mathbf {r} )}{\rho _{\alpha }(\mathbf {r} )+\rho _{\beta }(\mathbf {r} )}}\ .} <p> ζ = 0 {\displaystyle \zeta =0\,} corresponds to the diamagnetic spin-unpolarized situation with equal α {\displaystyle \alpha \,} and β {\displaystyle \beta \,} spin densities whereas ζ = ± 1 {\displaystyle \zeta =\pm 1} corresponds to the ferromagnetic situation where one spin density vanishes. The spin correlation energy density for a given values of the total density and relative polarization, <i>є</i>c(<i>ρ</i>,<i>ζ</i>), is constructed so to interpolate the extreme values. Several forms have been developed in conjunction with LDA correlation functionals.<a class="footnote-ref" id="fnref:12" href="#fn:12"><sup>12</sup></a> </p> <h2 id="exchange-correlation-potential">Exchange-correlation potential</h2> <p>The exchange-correlation potential corresponding to the exchange-correlation energy for a local density approximation is given by<a class="footnote-ref" id="fnref:13" href="#fn:13"><sup>13</sup></a> </p> v x c L D A ( r ) = δ E L D A δ ρ ( r ) = ϵ x c ( ρ ( r ) ) + ρ ( r ) ∂ ϵ x c ( ρ ( r ) ) ∂ ρ ( r ) . {\displaystyle v_{\rm {xc}}^{\mathrm {LDA} }(\mathbf {r} )={\frac {\delta E^{\mathrm {LDA} }}{\delta \rho (\mathbf {r} )}}=\epsilon _{\rm {xc}}(\rho (\mathbf {r} ))+\rho (\mathbf {r} ){\frac {\partial \epsilon _{\rm {xc}}(\rho (\mathbf {r} ))}{\partial \rho (\mathbf {r} )}}\ .} <p>In finite systems, the LDA potential decays asymptotically with an exponential form. This result is in error; the true exchange-correlation potential decays much slower in a Coulombic manner. The artificially rapid decay manifests itself in the number of Kohn–Sham orbitals the potential can bind (that is, how many orbitals have energy less than zero). The LDA potential can not support a Rydberg series and those states it does bind are too high in energy. This results in the highest occupied molecular orbital (<a href="/facts/HOMO/XyofpH0p">HOMO</a>) energy being too high in energy, so that any predictions for the <a href="/facts/Ionization_potential/lezUajfA">ionization potential</a> based on <a href="/facts/Koopmans%2527_theorem/H1W3vFAQ">Koopmans' theorem</a> are poor. Further, the LDA provides a poor description of electron-rich species such as <a href="/facts/Anion/wrbwhF3z">anions</a> where it is often unable to bind an additional electron, erroneously predicating species to be unstable.<a class="footnote-ref" id="fnref:14" href="#fn:14"><sup>14</sup></a> In the case of spin polarization, the exchange-correlation potential acquires spin indices. However, if one only considers the exchange part of the exchange-correlation, one obtains a potential that is diagonal in spin indices (in atomic units):<a class="footnote-ref" id="fnref:15" href="#fn:15"><sup>15</sup></a> </p><p> v x c , α β L D A ( r ) = δ E L D A δ ρ α β ( r ) = 1 2 δ α β δ E L D A [ 2 ρ α ] δ ρ α = − δ α β ( 3 π ) 1 / 3 2 1 / 3 ρ α 1 / 3 {\displaystyle v_{\rm {xc,\alpha \beta }}^{\mathrm {LDA} }(\mathbf {r} )={\frac {\delta E^{\mathrm {LDA} }}{\delta \rho _{\alpha \beta }(\mathbf {r} )}}={\frac {1}{2}}\delta _{\alpha \beta }{\frac {\delta E^{\mathrm {LDA} }[2\rho _{\alpha }]}{\delta \rho _{\alpha }}}=-\delta _{\alpha \beta }{\Big (}{\frac {3}{\pi }}{\Big )}^{1/3}2^{1/3}\rho _{\alpha }^{1/3}} </p> <h2 id="references">References</h2> <ol> <li id="fn:1"><p>Bagayoko, Diola (December 2014). "Understanding density functional theory (DFT) and completing it in practice". AIP Advances. 4 (12): 127104. Bibcode:2014AIPA....4l7104B. doi:10.1063/1.4903408. ISSN 2158-3226. <a href="https://doi.org/10.1063%2F1.4903408" target="_blank">https://doi.org/10.1063%2F1.4903408</a> <a href="#fnref:1" class="footnote-back-ref">↩</a></p></li> <li id="fn:2"><p>Segall, M.D.; Lindan, P.J (2002). "First-principles simulation: ideas, illustrations and the CASTEP code". Journal of Physics: Condensed Matter. 14 (11): 2717. Bibcode:2002JPCM...14.2717S. doi:10.1088/0953-8984/14/11/301. S2CID 250828366. <a href="/wiki/Bibcode_(identifier)" target="_blank">/wiki/Bibcode_(identifier)</a> <a href="#fnref:2" class="footnote-back-ref">↩</a></p></li> <li id="fn:3"><p>Assadi, M.H.N; et al. (2013). "Theoretical study on copper's energetics and magnetism in TiO2 polymorphs". Journal of Applied Physics. 113 (23): 233913–233913–5. arXiv:1304.1854. Bibcode:2013JAP...113w3913A. doi:10.1063/1.4811539. S2CID 94599250. <a href="/wiki/ArXiv_(identifier)" target="_blank">/wiki/ArXiv_(identifier)</a> <a href="#fnref:3" class="footnote-back-ref">↩</a></p></li> <li id="fn:4"><p>Zhao, G. L.; Bagayoko, D.; Williams, T. D. (1999-07-15). "Local-density-approximation prediction of electronic properties of GaN, Si, C, and RuO2". Physical Review B. 60 (3): 1563–1572. Bibcode:1999PhRvB..60.1563Z. doi:10.1103/physrevb.60.1563. ISSN 0163-1829. <a href="/wiki/Bibcode_(identifier)" target="_blank">/wiki/Bibcode_(identifier)</a> <a href="#fnref:4" class="footnote-back-ref">↩</a></p></li> <li id="fn:5"><p>Bagayoko, Diola (December 2014). "Understanding density functional theory (DFT) and completing it in practice". AIP Advances. 4 (12): 127104. Bibcode:2014AIPA....4l7104B. doi:10.1063/1.4903408. ISSN 2158-3226. <a href="https://doi.org/10.1063%2F1.4903408" target="_blank">https://doi.org/10.1063%2F1.4903408</a> <a href="#fnref:5" class="footnote-back-ref">↩</a></p></li> <li id="fn:6"><p>Parr, Robert G; Yang, Weitao (1994). Density-Functional Theory of Atoms and Molecules. Oxford: Oxford University Press. ISBN 978-0-19-509276-9. <a href="978-0-19-509276-9" target="_blank">978-0-19-509276-9</a> <a href="#fnref:6" class="footnote-back-ref">↩</a></p></li> <li id="fn:7"><p>Dirac, P. A. M. (1930). "Note on exchange phenomena in the Thomas-Fermi atom". Mathematical Proceedings of the Cambridge Philosophical Society. 26 (3): 376–385. Bibcode:1930PCPS...26..376D. doi:10.1017/S0305004100016108. <a href="/wiki/Paul_Dirac" target="_blank">/wiki/Paul_Dirac</a> <a href="#fnref:7" class="footnote-back-ref">↩</a></p></li> <li id="fn:8"><p>Parr, Robert G; Yang, Weitao (1994). Density-Functional Theory of Atoms and Molecules. Oxford: Oxford University Press. ISBN 978-0-19-509276-9. <a href="978-0-19-509276-9" target="_blank">978-0-19-509276-9</a> <a href="#fnref:8" class="footnote-back-ref">↩</a></p></li> <li id="fn:9"><p>Murray Gell-Mann and Keith A. Brueckner (1957). "Correlation Energy of an Electron Gas at High Density" (PDF). Phys. Rev. 106 (2): 364–368. Bibcode:1957PhRv..106..364G. doi:10.1103/PhysRev.106.364. S2CID 120701027. <a href="https://authors.library.caltech.edu/3713/1/GELpr57b.pdf" target="_blank">https://authors.library.caltech.edu/3713/1/GELpr57b.pdf</a> <a href="#fnref:9" class="footnote-back-ref">↩</a></p></li> <li id="fn:10"><p>D. M. Ceperley and B. J. Alder (1980). "Ground State of the Electron Gas by a Stochastic Method". Phys. Rev. Lett. 45 (7): 566–569. Bibcode:1980PhRvL..45..566C. doi:10.1103/PhysRevLett.45.566. S2CID 55620379. <a href="https://digital.library.unt.edu/ark:/67531/metadc1059358/" target="_blank">https://digital.library.unt.edu/ark:/67531/metadc1059358/</a> <a href="#fnref:10" class="footnote-back-ref">↩</a></p></li> <li id="fn:11"><p>Oliver, G. L.; Perdew, J. P. (1979). "Spin-density gradient expansion for the kinetic energy". Phys. Rev. A. 20 (2): 397–403. Bibcode:1979PhRvA..20..397O. doi:10.1103/PhysRevA.20.397. <a href="/wiki/Bibcode_(identifier)" target="_blank">/wiki/Bibcode_(identifier)</a> <a href="#fnref:11" class="footnote-back-ref">↩</a></p></li> <li id="fn:12"><p>von Barth, U.; Hedin, L. (1972). "A local exchange-correlation potential for the spin polarized case". J. Phys. C: Solid State Phys. 5 (13): 1629–1642. Bibcode:1972JPhC....5.1629V. doi:10.1088/0022-3719/5/13/012. S2CID 120985795. <a href="/wiki/Bibcode_(identifier)" target="_blank">/wiki/Bibcode_(identifier)</a> <a href="#fnref:12" class="footnote-back-ref">↩</a></p></li> <li id="fn:13"><p>Parr, Robert G; Yang, Weitao (1994). Density-Functional Theory of Atoms and Molecules. Oxford: Oxford University Press. ISBN 978-0-19-509276-9. <a href="978-0-19-509276-9" target="_blank">978-0-19-509276-9</a> <a href="#fnref:13" class="footnote-back-ref">↩</a></p></li> <li id="fn:14"><p>Fiolhais, Carlos; Nogueira, Fernando; Marques Miguel (2003). A Primer in Density Functional Theory. Springer. p. 60. ISBN 978-3-540-03083-6. <a href="978-3-540-03083-6" target="_blank">978-3-540-03083-6</a> <a href="#fnref:14" class="footnote-back-ref">↩</a></p></li> <li id="fn:15"><p>Giustino, Feliciano (2014). Materials Modelling Using Density Functional Theory: Properties and Predictions. Oxford University Press. p. 229. <a href="#fnref:15" class="footnote-back-ref">↩</a></p></li> </ol>