Question about the CHGCAR file

Queries about input and output files, running specific calculations, etc.


Moderators: Global Moderator, Moderator

Locked
Message
Author
Shenlan.Yang@monash.edu
Newbie
Newbie
Posts: 4
Joined: Tue Aug 18, 2020 5:58 am

Question about the CHGCAR file

#1 Post by Shenlan.Yang@monash.edu » Wed Dec 02, 2020 3:52 am

Dear all,

I have a question about the CHGCAR file and I hope you could give me some help. Though I noticed that there are already some discussions on CHGCAR unit, I'd like to double check for that. Based on the VASP guide, the charge values inside the CHGCAR file indicate the charge density ρ(r) multiplied by the grid volume defined by NGX*NGY*NGZ, however, I'm not quite sure about the physical meaning of this charge density ρ(r), does it represents the charge for one grid point or it shows the charge density defined by the real supercell volume? Thanks for the help!

Regards,
Elina

martin.schlipf
Global Moderator
Global Moderator
Posts: 542
Joined: Fri Nov 08, 2019 7:18 am

Re: Question about the CHGCAR file

#2 Post by martin.schlipf » Wed Dec 02, 2020 8:21 am

The CHGCAR file is setup in such a way that the typical integral
N = int d^3r rho(r),
where N is the number of electrons, is replaced by a discrete sum
N = sum_xyz rho(x,y,z) / (N_x N_y N_z),
where N_x N_y N_z, are the dimensions of the grid along the 3 lattice vectors. The rho(x,y,z) are the one you find in the CHGCAR file.

Martin Schlipf
VASP developer


Shenlan.Yang@monash.edu
Newbie
Newbie
Posts: 4
Joined: Tue Aug 18, 2020 5:58 am

Re: Question about the CHGCAR file

#3 Post by Shenlan.Yang@monash.edu » Thu Dec 03, 2020 3:08 am

Thank you for the help. And there is one thing I'm not quite sure, based on what you said, rho(x,y,z) should be the charge at (x,y,z) multiplied by (N_x N_y N_z), am I right?

martin.schlipf
Global Moderator
Global Moderator
Posts: 542
Joined: Fri Nov 08, 2019 7:18 am

Re: Question about the CHGCAR file

#4 Post by martin.schlipf » Thu Dec 03, 2020 8:11 am

Yes.

Martin Schlipf
VASP developer


Locked