jasp
jasp copied to clipboard
Spin Polarized charge densities
Hi John,
What is the role of spin=0 in get_charge_density? https://github.com/jkitchin/jasp/blob/master/jasp/volumetric_data.py#L45
If I call calc.get_charge_density(spin=1) for a spin polarized calculation, I get
Traceback (most recent call last):
File "
I think the charge density is the spin up + spin down density (atleast according the the CHGCAR documentation).
"For spinpolarized calculations, two sets of data can be found in the CHGCAR file. The first set contains the total charge density (spin up plus spin down), the second one the magnetization density (spin up minus spin down). For non collinear calculations the CHGCAR file contains the total charge density and the magnetisation density in the x, y and z direction in this order."
There is a self.chgdiff in VaspChargeDensity, that stores the difference between spin-up and spin-down charge densities. https://trac.fysik.dtu.dk/projects/ase/browser/trunk/ase/calculators/vasp.py#L1475
once upon a time in Dacapo you could get the spin up and down densities this way. I guess it doesn't work that way in VASP.