Contribution to the total Energy in VASP

Question on input files/tags, interpreting output, etc.

Please check whether the answer to your question is given in the VASP online manual or has been discussed in this forum previously!

Moderators: Global Moderator, Moderator

Post Reply
Message
Author
florent.boucher
Newbie
Newbie
Posts: 3
Joined: Wed Jan 12, 2005 3:31 pm
License Nr.: 123+paw

Contribution to the total Energy in VASP

#1 Post by florent.boucher » Mon Jul 11, 2005 4:06 pm

Dear Vasp users, I would like to have more details about the different crontributions to the total energy in the OUTCAR and espcially the "alpha Z" and "Ewald energy". Does any one knows exactly what they are ? Is the "Ewald energy" similar to a Madelung contribution ?
Regards
Florent
Last edited by florent.boucher on Mon Jul 11, 2005 4:06 pm, edited 1 time in total.

admin
Administrator
Administrator
Posts: 2922
Joined: Tue Aug 03, 2004 8:18 am
License Nr.: 458

Contribution to the total Energy in VASP

#2 Post by admin » Tue Jul 12, 2005 10:19 am

alpha Z and the Ewald energy define the electrostatic
interaction of the ions in a compensating electron gas.
The alpha Z component deals with the divergent parts (G=0).

The following parts are the Hartree and exchange correlation energy as
defined in the Kohn-Sham Hamiltonian. The entropy part stems from the
smearing (using the free energy as variational parameter, electronic entropy),
EBANDS from Kohn-Sham eigenvalues, and EATOM is the reference energy for the potential (which is defined in the POTCAR file).

Furthermore, you find a documentation of the respective terms of the Schroedinger
equation in the handouts of the VASP-workshop,
http://cms.mpi.univie.ac.at/vasp-workshop
especially in the session about dft (download especially
comput_mat.pdf, dft_introd.pdf, and accuracy.pdf (p.4 for EATOM))
Last edited by admin on Tue Jul 12, 2005 10:19 am, edited 1 time in total.

user
Newbie
Newbie
Posts: 30
Joined: Wed Mar 23, 2005 10:18 pm
License Nr.: 74

Contribution to the total Energy in VASP

#3 Post by user » Wed Mar 22, 2006 10:46 pm

Dear VASP Admin/Users,

1. I am wondering if the "PAW double counting" term in the OUTCAR has any physical meaning?

In PRB59, 1758 (1999), it is said that "in many band-structure codes, the total energy is evaluated as the sum of the Kohn-Sham eigenvalues minus double counting corrections". Is this the "PAW double counting" term that we see in the OUTCAR?

From the equation (48) in the above article, one can see that the "double counting" term has something to do with Hartree, exchange-correlation energies... Then,

2. How should we define the Madelung energy? i.e., what terms should we add up to get the Madelung energy?

I would be very glad if someone could tell me more about these matters.

Thanks!



<span class='smallblacktext'>[ Edited Wed Mar 22 2006, 11:50PM ]</span>
Last edited by user on Wed Mar 22, 2006 10:46 pm, edited 1 time in total.

User avatar
eariel
Newbie
Newbie
Posts: 16
Joined: Sat Sep 13, 2008 6:45 pm
License Nr.: 490

Contribution to the total Energy in VASP

#4 Post by eariel » Thu Jan 21, 2010 3:42 pm

[quote author=0).

The following parts are the Hartree and exchange correlation energy as
defined in the Kohn-Sham Hamiltonian. The entropy part stems from the
smearing (using the free energy as variational parameter, electronic entropy),
EBANDS from Kohn-Sham eigenvalues, and EATOM is the reference energy for the potential (which is defined in the POTCAR file).

Furthermore, you find a documentation of the respective terms of the Schroedinger
equation in the handouts of the VASP-workshop,
http://cms.mpi.univie.ac.at/vasp-workshop
especially in the session about dft (download especially
comput_mat.pdf, dft_introd.pdf, and accuracy.pdf (p.4 for EATOM))

[/quote]
Last edited by eariel on Thu Jan 21, 2010 3:42 pm, edited 1 time in total.

User avatar
eariel
Newbie
Newbie
Posts: 16
Joined: Sat Sep 13, 2008 6:45 pm
License Nr.: 490

Contribution to the total Energy in VASP

#5 Post by eariel » Thu Jan 21, 2010 3:46 pm

Hi,

[quote="admin"]alpha Z and the Ewald energy define the electrostatic
interaction of the ions in a compensating electron gas.
The alpha Z component deals with the divergent parts (G=0).

What means "The alpha Z component deals with the divergent parts (G=0)". Isn't the divergent term infinity?
Last edited by eariel on Thu Jan 21, 2010 3:46 pm, edited 1 time in total.

tlchan

Contribution to the total Energy in VASP

#6 Post by tlchan » Thu Jan 21, 2010 9:00 pm

The electron-ion, Hartree, and ion-ion interaction energies diverge individually. However, if you add the G=0 divergences up, it is a finite number. It can be shown that if the ionic potential is Coulomb, then the cancellation of the three G=0 terms is exact, i.e. the finite number is zero. Since VASP uses pseudopotentials, the summation of the three G=0 terms results in a non-zero finite number, which is called alpha.
Last edited by tlchan on Thu Jan 21, 2010 9:00 pm, edited 1 time in total.

ffdeli
Newbie
Newbie
Posts: 3
Joined: Thu Dec 16, 2010 10:29 pm

Contribution to the total Energy in VASP

#7 Post by ffdeli » Wed Jan 26, 2011 4:29 pm

[quote author=0).

[/quote]

Hi Admin,

Do you mean alpha Z and Ewald energy give the total electrostatic energy of the system? or I think of it in a wrong way? I want to know whether we can decompose the total energy from vasp into Kinetic, Coulomb, exchange and correlation energies, separately.

I was calculating adsorption of noble gas (Xe) on metal surface (111) at two different sites, of which the total energy is only slightly different:

The energies at site 1:

---------------------------------------------------
alpha Z PSCENC = 1057.06630392
Ewald energy TEWEN = 128376.65817811
-1/2 Hartree DENC = -145347.23431284
-exchange EXHF = 0.00000000
-V(xc)+E(xc) XCENC = 41.64037121
PAW double counting = 12658.25300008 -12074.80631272
entropy T*S EENTRO = -0.00031446
eigenvalues EBANDS = -575.59189787
atomic energy EATOM = 15740.57257888
---------------------------------------------------
free energy TOTEN = -123.44240568 eV


The energy at site 2:
---------------------------------------------------
alpha Z PSCENC = 1056.91493097
Ewald energy TEWEN = 128376.45871064
-1/2 Hartree DENC = -145344.66902862
-exchange EXHF = 0.00000000
-V(xc)+E(xc) XCENC = 41.63693035
PAW double counting = 12658.30538857 -12074.85364512
entropy T*S EENTRO = -0.00032545
eigenvalues EBANDS = -577.80127102
atomic energy EATOM = 15740.57257888
---------------------------------------------------
free energy TOTEN = -123.43573079 eV


The difference of the total energies at two different sites is only 0.0067 eV, which is reasonable. However, the difference of the (alpha Z+Ewald energy) is about 0.35 eV, which does not make sense to me. Since it is a weakly bound system, I expect the difference of the electrostatic energies at these two different sites to be very small, several meV or less. Therefore, I do not understand the (alpha Z+Ewald energy) energy shown above.

Any comments will be appreciated!


Chen
Last edited by ffdeli on Wed Jan 26, 2011 4:29 pm, edited 1 time in total.

tlchan

Contribution to the total Energy in VASP

#8 Post by tlchan » Thu Jan 27, 2011 3:10 pm

The electrostatic energy of a system should include the electrostatic interaction between ion-ion (Ewald), electron-electron (Hartree), and ion-electron. The information about ion-electron electrostatic interaction is embedded inside the eigenvalues.

Note that the individual number for the Ewald energy, and Hartree energy by themselves is not really meaningful. Only the total electrostatic energy of the system is a meaningful number. This is because in a plane wave code, the reference zero of the electrostatic potential is defined to be the average electrostatic potential of the unit cell, rather than defining the zero to be the "vacuum."
Last edited by tlchan on Thu Jan 27, 2011 3:10 pm, edited 1 time in total.

ffdeli
Newbie
Newbie
Posts: 3
Joined: Thu Dec 16, 2010 10:29 pm

Contribution to the total Energy in VASP

#9 Post by ffdeli » Fri Jan 28, 2011 4:06 pm

This is really helpful, thanks very much!
Last edited by ffdeli on Fri Jan 28, 2011 4:06 pm, edited 1 time in total.

juhL

Contribution to the total Energy in VASP

#10 Post by juhL » Mon May 30, 2011 4:21 pm

Dear all,

i am in principal trying to do the same as ffdeli, so i want to compare kinetic energy, electrostatic energy and exchange-correlation energy. however i am still not exactly sure, how to do this, and what each contribution to the free energy mean.

alphaZ + Ewald: electrostatic energy of ion-ion interaction
-1/2 Hartree: one part of the first double counting corrections, i assume to get E_H, one takes this value times (-1/2) ?
-exchange: I can't find out what this is about. Is it E_XC (according to the admin's statement)? then why is it included here? it should be part of the eigenvalues and there is no double counting correction for this contribution.
-V(xc)+E(xc): rest of the first double counting corrections
PAW double counting: double counting corrections 2&3
entropy T*S: is clear
eigenvalues: the one electron energies; here E_kin, E_electrostatic (electron electron and electron core) and E_XC are embedded
atomic energy: in principal E_kin of the cores

my questions now are:

1) is my understanding of the different contributions right?

2) if i am right so far, i would like to know where this (in my opinion additional) "-exchange" is coming from and why it is included into the calculation of the free energy.

3) is there any possibility to separate now the energy into E_kin, E_electrostatic and E_XC? Or am i right, that E_kin as well as the ion-electron electrostatic energy is embedded into the eigenvalues and one just can't do this separation into this three energy contributions?







<span class='smallblacktext'>[ Edited Mon May 30 2011, 09:22PM ]</span>
Last edited by juhL on Mon May 30, 2011 4:21 pm, edited 1 time in total.

tlchan

Contribution to the total Energy in VASP

#11 Post by tlchan » Tue Jun 07, 2011 6:42 pm

[quote="juhL"]
atomic energy: in principal E_kin of the cores
[/quote]
3) is there any possibility to separate now the energy into E_kin, E_electrostatic and E_XC? Or am i right, that E_kin as well as the ion-electron electrostatic energy is embedded into the eigenvalues and one just can't do this separation into this three energy contributions?
[/quote]
I think VASP already calculated E_XC. You can modify the source code slightly to print it out.

The electrostatic energy is the sum of ion-electron electrostatic interaction, the Hartree energy, the Ewald ion-ion, and the alpha term. As already noted by you, the terms are printed out already, except the ion-electron term. The ion-electron term is the expectation value of the pseudopotential. I don't know if VASP calculated it somewhere. If it doesn't, you've to do some coding to obtain the value.

For the kinetic energy, it is the expectation value of the Laplacian operator. I believe VASP should have calculated it somehow. You've to look at the source code to find out. Note the this kinetic energy is the single-particle kinetic energy. The many-particle contribution of the kinetic energy is within the correlation energy.
Last edited by tlchan on Tue Jun 07, 2011 6:42 pm, edited 1 time in total.

juhL

Contribution to the total Energy in VASP

#12 Post by juhL » Wed Jun 08, 2011 12:49 pm

hi,

okay, i was kind of hoping to get around of getting into the source code, but thanks a lot, that helps!
Last edited by juhL on Wed Jun 08, 2011 12:49 pm, edited 1 time in total.

Post Reply