The density, including the cores?

MULTIBINIT, aTDEP, ANADDB, optics, cut3d, mrg* postprocessors

Moderators: MMNSchmitt, gonze

Locked
Dr Eigen
Posts: 2
Joined: Wed Jun 08, 2011 6:00 pm

The density, including the cores?

Post by Dr Eigen » Wed Jun 08, 2011 6:11 pm

Hi folks

I would like to get the total electron density of my system. The density in the DEN file integrates to the number of valence electrons.

If I'm using pre-rolled pseudopotentials (not PAW) downloaded from the abinit website, is there an easy way to get a density including the core electrons?

Ciao
Terry

nnikola
Posts: 8
Joined: Wed Oct 05, 2011 11:26 am

Re: The density, including the cores?

Post by nnikola » Wed Oct 05, 2011 11:44 am

Well, according to manual, the way to achieve this is to use prtden keyword (PAW only):

'prtden=3 causes generation of a file _PAWDEN that contains the bulk full charge density (valence+core)'.

The catch is that you need to provide external files in the working directory, with atomic core and valence charge (two-column format) for every distinct atom type in your structure .
I would be grateful if someone can answer me from just how to generate these files?

ilukacevic
Posts: 271
Joined: Sat Jan 16, 2010 12:05 pm
Location: Dept. of Physics, University J. J. Strossmayer, Osijek, Croatia
Contact:

Re: The density, including the cores?

Post by ilukacevic » Wed Oct 05, 2011 12:08 pm

Hi!

I'm not completely sure, but I think that you can find these at
http://www.abinit.org/downloads/all_core_electron

Cheers!

Igor L.

nnikola
Posts: 8
Joined: Wed Oct 05, 2011 11:26 am

Re: The density, including the cores?

Post by nnikola » Fri Oct 07, 2011 11:30 am

Igor,

These are the all-electron files for Hirschfeld charge analysis.
Manual clearly states there must be two (core + valence) files per atom for prtden .eq. 3 to work.
Abinit stops with error if these files are not present in working directory.
Also, can someone explain why the separate valence-charge file is needed?

Regards,

Nikola

Locked