Hai Abinit users,
I would like to know how can I plot the charge density around the Fermi level. Using cut3d, I generated the .xsf file and I got only the total density. I wish to know whether I could get the energy projected density using the cut3d utility of Abinit.
Any comments on the above are greatly appreciated.
Thanks in advance
Seba
Electronic charge density around the Fermi level
Moderators: MMNSchmitt, gonze
-
Seba Darshan
- Posts: 63
- Joined: Mon Feb 17, 2014 5:19 am
Re: Electronic charge density around the Fermi level
Hi Seba,
You can get the density associated with each band by doing a multi-dataset run, using
then you will get <n> different _DEN files; each _DEN will be associated with a given band. Now you can process them with CUT3D and get the "charge density-per-band". I guess what you need is that associated with the upper occupied band.
Cheers,
Temok
You can get the density associated with each band by doing a multi-dataset run, using
Code: Select all
ndataset <n> # as many datasets as (occupied bands * nspinor)
charge # set such that you establish charge balance in your system
occopt 0
# Set occ such that you turn off all but one band.
# put as many columns as (occupied bands * nspinor)
# for example, for nspinor = 1,
occ1 2.0 0.0 0.00 ...
occ2 0.0 2.0 0.00 ...
occ3 0.0 0.0 2.00 ...
.
.
.
occ<n> 0.0 0.0 .... 2.0
then you will get <n> different _DEN files; each _DEN will be associated with a given band. Now you can process them with CUT3D and get the "charge density-per-band". I guess what you need is that associated with the upper occupied band.
Cheers,
Temok