The ABINIT Users Mailing List ( CLOSED )

[Shu Heng Pau <gpau@MIT.EDU>]

Hi all,

I am puzzled by the complex wavefunction (in real space) that I obtained while
performing 0 k-point calculation.  I have checked that the imaginary part is
not due to a shift of the form exp(i\alpha) and its' magnitude is of the same
order as the real component.  For a 0 k-point calculation, shouldn't the
wavefunction be real?

I attached my input file: it is basically a simulation of simple diamond
structure for Carbon.  The 4th band that I obtain is complex.

Thanks,
George

# Definition of the unit cell
acell 10.2 10.2 10.2
rprim  0.0  0.5  0.5
       0.5  0.0  0.5
       0.5  0.5  0.0
#Definition of the atom types
ntypat 1
znucl 6.0
#Definition of the atoms
natom 2
typat 1 1
xred 0.0 0.0 0.0
     1/4 1/4 1/4
#Definition of the planewave basis set
ecut 10
ecutsm 0.5
ngfft 50 50 50
nband 4
intxc 1
nkpt 1
kpt 0 0 0
nstep 100
toldfe 1.0d-7
occopt 1
iscf 5
ixc 6
istwfk 1
n1xccc  0
ionmov 2
optcell 2
diemac 12.0
dilatmx 2