The ABINIT Users Mailing List ( CLOSED )

[Z_Lin <applelinux@gmail.com>]

> Ok - excellent question. Yes: the states will be different due to the
> differing high T density you get, and the perturbed potential will be too.
> This is definitely the best approximation you can get from the standard
> code. There could be additional terms coming from the perturbation of the
> entropic term, but this is probably less of a worry (small and probably
> fine). Setting the smearing in elphon to the same value is a bit off, as
> it uses a gaussian smearing (only supposed to be doing kpoint integration,
> not real temperature). You could modify the routine (probably elphon.F90)
> to use a Fermi Dirac smearing function, or adjust the smearing width to
> get a rough equivalent to the thermal function.
I think you mean the occupation number from previous run at some high
electron temperature? The Fermi function will have to be modified for
a change in the chemical potential (fermi level) which is a function
of electron temperature. Is it taken into account in the code?

Zhibin