Core charge

Bader charge density analysis

Moderator: moderators

danilo
Posts: 2
Joined: Thu Oct 07, 2010 12:55 pm

Re: Core charge

Post by danilo »

[/quote]You should add AECCAR0 and AECCAR2.[/quote]

There is a mistake in my post. I add AECCAR0 and AECCAR2.
But the result of Bader analysis on CHGCAR_sum is 1621.4448 core electrons
graeme
Site Admin
Posts: 2253
Joined: Tue Apr 26, 2005 4:25 am
Contact:

Re: Core charge

Post by graeme »

You won't be able to accurately calculate the total charge of the core electrons without an extremely fine grid. But I don't think this is what you want to do anyway. You probably want to do the Bader partitioning on the charge density sum (that you calculated) and the integration on the valance electrons:

bader CHGCAR -ref CHGCAR_sum

A little more information about this can be found here: http://theory.cm.utexas.edu/vtsttools/bader/vasp.php
cenwanglai
Posts: 46
Joined: Sun Jun 13, 2010 5:46 am

Re: Core charge

Post by cenwanglai »

"It is essential to increase the FFT grid (NGXF, NGYF, and NGZF) until the core charge converges to the expected value. "

From where I can see that the core charge converges to the expected value?

The standard output on the screen? or by

bader AECCAR0

??
graeme
Site Admin
Posts: 2253
Joined: Tue Apr 26, 2005 4:25 am
Contact:

Re: Core charge

Post by graeme »

This thread started many years ago and this statement about requiring convergence of the core charge is no longer necessary. Instead follow the instructions at:
http://theory.cm.utexas.edu/bader/vasp.php
yxcai
Posts: 3
Joined: Thu Apr 07, 2011 5:50 am

Re: Core charge

Post by yxcai »

Why is the total electrons not correct? I have set LAECHG=TRUE. and chosen PAW psudopotential. I also find that the total electrons chang with the cutoffenergy increaseing.

H C N Si
32 6 4 192

32+ 6*6+4*7+192*14=2784

but the total electron is 1674

why does this error happen?
graeme
Site Admin
Posts: 2253
Joined: Tue Apr 26, 2005 4:25 am
Contact:

Re: Core charge

Post by graeme »

Are you integrating the valence charge? I'm going to cut-and-paste from: http://theory.cm.utexas.edu/bader/vasp.php
Note that the partitioning is done on the total charge (CHGCAR_sum) but the integration is done on the normal CHGCAR file from vasp, containing only the valence electrons.

-------

By adding the LAECHG=.TRUE. to the INCAR file, the core charge is written to AECCAR0 and the valance charge to AECCAR2. These two charge density files can be summed using the chgsum.pl script:

chgsum.pl AECCAR0 AECCAR2
The total charge will be written to CHGCAR_sum.
The bader analysis can then be done on this total charge density file:

bader CHGCAR -ref CHGCAR_sum

-------
lgxyz
Posts: 51
Joined: Tue Sep 07, 2010 1:52 pm

Re: Core charge

Post by lgxyz »

On the webpage: http://theory.cm.utexas.edu/bader/vasp.php, the last sentence says "One finally note is that you need a fine fft grid to accurately reproduce the correct total core charge. It is essential to do a few calculations, increasing NG(X,Y,Z)F until the total charge is correct. " How to increasing the NG(X,Y,Z)F ? Pls give some suggestions. Thanks!
graeme
Site Admin
Posts: 2253
Joined: Tue Apr 26, 2005 4:25 am
Contact:

Re: Core charge

Post by graeme »

set:
NGXF=xxx
NGYF=yyy
NGZF=zzz
in the INCAR file. Choose values which are (for example) twice as large as those reported in the OUTCAR from a standard calculation.
lgxyz
Posts: 51
Joined: Tue Sep 07, 2010 1:52 pm

Re: Core charge

Post by lgxyz »

I have set NGX(Y,Z)F to be twice as the values in the standard calculations. After the calculation of "bader CHGCAR -ref CHGCAR_sum", I found the "NUMBER OF ELECTRONS" in ACF.dat is the number of valence electrons of the system. Is that correct? Where is the number of the total electron, including the core and valence electron. Thanks!
graeme
Site Admin
Posts: 2253
Joined: Tue Apr 26, 2005 4:25 am
Contact:

Re: Core charge

Post by graeme »

read this: http://theory.cm.utexas.edu/bader/vasp.php

You are integrating the valence charge density, so yes, the total number of electrons is the total number of valence electrons.

If you ran 'bader CHGCAR_sum' you would be analyzing the total charge, but the values would have large errors because of the cusps in the total charge density.
protonic
Posts: 3
Joined: Fri Apr 20, 2012 7:43 am

Re: Core charge

Post by protonic »

To graeme: you posted above (Mar 31, 2011) "..this statement about requiring convergence of the core charge is no longer necessary". Yet on Sep 29, 2011 you posted "Choose values [of NG(X,Y,Z)F] which are (for example) twice as large as those reported in the OUTCAR from a standard calculation." Are you saying that the tests changing NG(X,Y,Z)F are no longer necessary, or are they necessary?
graeme
Site Admin
Posts: 2253
Joined: Tue Apr 26, 2005 4:25 am
Contact:

Re: Core charge

Post by graeme »

The best procedure is to include the core charges in the total density and then integrate over the valence charge. This does not require sufficient resolution to resolve the cusps at the core; you should, however, increase the grid density to see convergence in the integrated valence charges.
emma
Posts: 1
Joined: Mon Nov 12, 2012 12:44 pm

Re: Core charge

Post by emma »

This forum is really awesome!
zwmiracle
Posts: 3
Joined: Sat Jan 05, 2013 7:11 am

Re: Core charge

Post by zwmiracle »

[quote="graeme"]It is now possible to get core charges from the VASP code. This is important for Bader calculations in which most or all of the valance electrons participate with bonding.

Getting the core charge density from a PAW calculation requires updating to the latest version of the vasp code (4.6.31 08Feb07). With this code, the core and valance charge density will be written if the LAECHG = .TRUE. tag is set in the INCAR file. It is essential to increase the FFT grid (NGXF, NGYF, and NGZF) until the core charge converges to the expected value.

The core charge is written to the AECCAR0 and the valance to AECCAR2. These can be added using our chgsum.pl script. The bader analysis on this total charge density will not have any pseudopotential errors.[/quote]

HELLO~
What is the difference between:
1. "bader CHGCAR " directly, after self-consistence in VASP without LAECHG setting
2. set LAECHG=ture, then use AECCAR...blablabla
?

the first one also gives me the ACF.dat, which is sometimes the same with the second one.
What do the two ways differ? The easier way and the more complex way...

THANK YOU !!! ^_^
zwmiracle
Posts: 3
Joined: Sat Jan 05, 2013 7:11 am

Re: Core charge

Post by zwmiracle »

[quote="zwmiracle"][quote="graeme"]It is now possible to get core charges from the VASP code. This is important for Bader calculations in which most or all of the valance electrons participate with bonding.

Getting the core charge density from a PAW calculation requires updating to the latest version of the vasp code (4.6.31 08Feb07). With this code, the core and valance charge density will be written if the LAECHG = .TRUE. tag is set in the INCAR file. It is essential to increase the FFT grid (NGXF, NGYF, and NGZF) until the core charge converges to the expected value.

The core charge is written to the AECCAR0 and the valance to AECCAR2. These can be added using our chgsum.pl script. The bader analysis on this total charge density will not have any pseudopotential errors.[/quote]

HELLO~
What is the difference between:
1. "bader CHGCAR " directly, after self-consistence in VASP without LAECHG setting
2. set LAECHG=ture, then use AECCAR...blablabla
?

the first one also gives me the ACF.dat, which is sometimes the same with the second one.
What do the two ways differ? The easier way and the more complex way...

THANK YOU !!! ^_^[/quote]


OH .the former only contains valent electrons, and the latter contains core electrons too?
so if the test shows that my system is identital to these two ways, then i can use the easier one without worry?
sorry it is really not easy to get along with all your discussions....
Post Reply