Density Functional Theory and the Local Density Approximation

In summary, the conversation revolved around calculating the ground state energy of Helium using density functional theory and the local density approximation. The individual seeking help had set up universal functionals and was struggling with finding a reasonably accurate exchange-correlation term. They were also working on solving a single-electron Schrodingers equation with an external Coulombic potential and a Hartree potential to get a new density. The expert suggested using a concrete expression for the LDA correlation energy given in previous papers, and using Slater/Dirac exchange as an xc functional for a simpler implementation. The individual had some difficulties with convergence and obtaining the correct energy values, but the expert offered potential solutions and reassured them that the density functional derivatives may be complex but it
  • #1
JMedley
8
0
Im trying to calculate the ground state energy of Helium using a density functional theory approach combined with the local density approximation. So far I have set up universal functionals and I mainly need help with the actual algorithm the evaluation of the Hartree energy functional.
 
Physics news on Phys.org
  • #2
You need to be more specific. I understand neither what you are trying to achieve, how you are trying to achieve it, nor what exactly is the problem.
 
  • #3
As to what I'm trying to achieve:
"I'm trying to calculate the ground state energy of Helium"
As to how I plan on doing it:
"density functional theory"
And admittedly I wasn't particularly clear on the problems but in particular finding a reasonably accurate exchange-correlation term. I'd read papers where people refer to Monte-Carlo methods for finding a good correlation energy but haven't been able to find any literature on it
 
  • #4
JMedley said:
As to what I'm trying to achieve:
"I'm trying to calculate the ground state energy of Helium"
As to how I plan on doing it:
"density functional theory"
Great. That totally clarifies it. So let me give you a solution to your problem:
(1) Get an account on a computer with molpro installed.
(2) Make a textfile he.com:
memory,50,m;
geometry={He};
basis=vtz
ks,lda
(3) execute "molpro he.com".
...tada, you've just calculated the LDA ground state energy of He. Sure you don't want to go into somewhat more detail about what you are trying to achive and how?

And admittedly I wasn't particularly clear on the problems but in particular finding a reasonably accurate exchange-correlation term.

Your topic says you want to use LDA. LDA implies a specific exchange-correlation term, you don't need to find it. Just get it from the VWN parameterization; or use the one of Perdew et al.

I'd read papers where people refer to Monte-Carlo methods for finding a good correlation energy but haven't been able to find any literature on it
Correlation energy of what? Monte Carlo methods are used to calculate the properties of the homogeneous and non-homegenous electron gas. This data is used in parameterizations of ab initio density functionals (like PBE, for example). If you intend to apply an existing functional, like LDA, you do not need that because it's already parameterized for you.
 
  • #5
Ok, I guess I underestimated this one! Thus far I have an exchange-correlation term of:
[tex] E_{xc}= -\frac{3}{4}\sqrt[3]{\frac{3n(r)}{\pi}}+\frac{0.44}{r_s+7.8} [/tex]
Where [tex] n(r) [/tex] is the density of electrons and [ tex ] r_s [/tex] is the Wigner radius. Is this an acceptable expression? I have then tried to use the functional derivative of this to get the exchange-correlation potential as:
[tex] \frac{\delta E_{xc}[n(r)]}{\delta n(r)} = \frac{d}{dn}(n(r)E_{xc}) [/tex]
This is then fed, along with an external Coulombic potential and a hartree potential, back into solving a single-electron Schrodingers equation to get a new density and so on and so on...
Is that the gist of it? Having calculated this differential I know it turns out very messy so I doubt I'm on track! Thanks for your help
 
  • #6
JMedley said:
Ok, I guess I underestimated this one! Thus far I have an exchange-correlation term of:
[tex] E_{xc}= -\frac{3}{4}\sqrt[3]{\frac{3n(r)}{\pi}}+\frac{0.44}{r_s+7.8} [/tex]
Where [tex] n(r) [/tex] is the density of electrons and [ tex ] r_s [/tex] is the Wigner radius. Is this an acceptable expression?
Your exchange correlation term has the right format, but it does look somewhat too simple to me. I don't think it has the correct asymptotics in the RPA limit (shouldn't there be a ln() in there?).

Concrete expressions for the LDA correlation energy were given in
http://dx.doi.org/10.1103/PhysRevB.45.13244 (Perdew-Wang parameterization)
http://dx.doi.org/10.1139/p80-159 (VWN parameterization)

Your expression may be a viable parameterization for your concrete task. In order to check that you could compare it with the expressions given in these^ papers. There was also a nice explanation (including giving all the leading terms) in Schwabl's textbook "Advanced Quantum Mechanics".

I have then tried to use the functional derivative of this to get the exchange-correlation potential as:
[tex] \frac{\delta E_{xc}[n(r)]}{\delta n(r)} = \frac{d}{dn}(n(r)E_{xc}) [/tex]
This is then fed, along with an external Coulombic potential and a hartree potential, back into solving a single-electron Schrodingers equation to get a new density and so on and so on...
Is that the gist of it? Having calculated this differential I know it turns out very messy so I doubt I'm on track! Thanks for your help

Your track seems right. Note that density functional derivatives tend to be very ugly when done correctly. So if you get long and cumbersome expressions I'd rather interpret that as a good sign :). Anyway: In practice density functional derivatives are usually implemented in an automatic fashion. That is, only the energy expression is explicitly programmed, and then code for evaluating it and its derivatives is generated from a computer algebra system or sythesized via automatic differentiation.

If you just want to see if you got the concept right, you could concentrate on just using Slater/Dirac exchange as xc functional. This expression is very simple and allows for a non-messy implementation, including derivatives. And if you get it right, you can just assume that it would work for other functionals too, in principle.
 
  • #7
Cool, So I've got it so that it converges on an energy (having used a simple mixing scheme). Unfortunately its the wrong energy! I've solved Possions equation in order to get the Hartree energy of the Helium atom and used the exchange-correlation functional of:

[tex] V_{xv}[n(r)] = -\left(\frac{3}{2\pi}\right)^{\frac{2}{3}}\frac{1}{r_s}[1+0.0545r_sln(1+11.4/r_s)] [/tex]

Where [tex]r_s[/tex] is the Wigner radius. So far as I know I should be getting a value of about 5.9Ry but its converging on about 12.1Ry. The ODE I'm integrating to get the Hartree energy is:

[tex] \frac{d^2S}{dr^2}=-4\pi rn(r) [/tex]

Where [tex] S=rV_H [/tex] and [tex]V_H[/tex] is the Hartree potential. The Coulombic potential is still:

[tex]V_{ext}=-\frac{4}{r}[/tex]

Do these look OK? I can't really figure out what the problem is. Cheers. Jack
 
Last edited:
  • #8
The He/LDA energy is supposed to be somewhat like -2.834 Hartree (that's the atomic unit; one Hartree is two Ry), with about -0.975 H coming from the functional.

The most common error people make then writing their first SCF program is to sum up the energy wrong. The energy is *not* the sum of the orbital eigenvalues, but rather it is <(1/2)(h + f)> where h is the core Hamiltonian (kinetic energy + nuclear potential) plus the fock matrix f (<f> is the sum of the eigenvalues, <> being the expectation value)

About the external potential: shouldn't this be -2/r for He?
 
  • #9
Well because we've been working in Rydberg units I thought the External potential was:
[tex]-\frac{2Z}{r}[/tex]
In that case I think we may be summing the wrong energy; at the moment we're inputting an initial density of double the Hydrogen 1s density. We then use this to get a hartree potential and an exchange-corellation potential. These are then input to Schrodingers equation with the Coulombic potential to get a new density and so on and so on...Presently we are just looking for convergence of the eigenenergy the comes from Schrodingers equation each time. Is this wrong?
 

What is Density Functional Theory (DFT)?

Density Functional Theory (DFT) is a method used in quantum mechanics to calculate the electronic structure of atoms, molecules, and solids. It is based on the concept of the electron density, which is a fundamental property of a system that contains all the information about the electronic structure.

What is the Local Density Approximation (LDA)?

The Local Density Approximation (LDA) is a commonly used approximation in Density Functional Theory (DFT) calculations. It assumes that the exchange-correlation energy, which describes the interactions between electrons, is only dependent on the local density of electrons at a particular point in space. This simplifies the calculations, making them more computationally efficient.

What are the limitations of DFT and LDA?

One limitation of DFT is that it is based on approximations and does not always accurately capture the behavior of strongly correlated systems. Additionally, the LDA approximation does not account for spatial variations in the exchange-correlation energy, leading to inaccuracies in some systems.

How is DFT and LDA used in materials science?

DFT and LDA are commonly used in materials science to study the electronic and structural properties of materials. They can be used to predict the properties of new materials and to understand the behavior of existing materials under different conditions. They are also used to guide experimental research and to design new materials with specific properties.

Are there any alternatives to DFT and LDA?

There are several alternative methods to DFT and LDA, such as Hartree-Fock theory and hybrid functionals. These methods may provide more accurate results for certain systems, but they are also more computationally demanding. Researchers often choose the appropriate method based on the specific system and the level of accuracy required.

Similar threads

Replies
1
Views
658
  • Quantum Physics
Replies
1
Views
737
  • Quantum Physics
Replies
21
Views
2K
Replies
1
Views
647
Replies
3
Views
902
  • Atomic and Condensed Matter
Replies
1
Views
515
Replies
27
Views
2K
  • Atomic and Condensed Matter
Replies
0
Views
958
  • Quantum Physics
Replies
24
Views
1K
  • Quantum Physics
Replies
11
Views
1K
Back
Top