SPIN-Orbit calculation dumps only one spin channel in EIGENVAL, DOSCAR files
Posted: Mon Oct 16, 2006 7:45 pm
Hello,
I am new with the code so I might be seeing this as a bug that is not really present. I'm using VASP 4.6.28, compiled with pgi fortran and running an a xeon cluster with cent-Os (free fedora)
I have made a calculation on bulk gold both with ISPIN=1 and 2. ISPIN=2 gives me results for two spin channels, as expected.
Then I turn on the SPIN-ORBIT coupling. And the output dumps only one eigenchannel. My INCAR file follows:
===========================================================
Bulk gold test calculation
ISPIN=2
SYSTEM = AuBulk
NWRITE = 2
PREC=Accurate
IALGO=48
GGA=91
LSORBIT = .TRUE.
SAXIS = 0 0 1 ! direction of the magnetic field
NBANDS = 50
NUPDOWN = 1
===========================================================
Standard output is just normal:
[salva@tbird AuBulk]$ ./vasp
running on 1 nodes
distr: one band on 1 nodes, 1 groups
vasp.4.6.28 25Jul05 complex
POSCAR found : 1 types and 1 ions
LDA part: xc-table for Ceperly-Alder, standard interpolation
generate k-points for: 12 12 12
|| W W A A R R N N II N N GGGG !!! || || Full k-point grid generated || Inversion symmetry is not applied || | -----------------------------------------------------------------------------
POSCAR, INCAR and KPOINTS ok, starting setup
WARNING: wrap around errors must be expected
FFT: planning ... 1
reading WAVECAR
entering main loop
N E dE d eps ncg rms rms(c)
.
.
.
RMM: 18 -0.350209045025E+01 0.16626E-04 -0.10425E-06 32337 0.238E-03 1 F= -.35020905E+01 E0= -.35020912E+01 d E =0.224021E-05 mag= 0.0003 0.0003 0.0003
writing wavefunctions
FORTRAN STOP
===========================================================
These are the first entries of the EIGENVAL file:
1 1 1 1
0.1825866E+02 0.2955706E-09 0.2955706E-09 0.2955706E-09 0.5000000E-15
1.0000000000000000E-004
CAR
AuBulk
11 364 50
0.3742002E-15 0.3742002E-15 0.3742002E-15 0.5787037E-03
1 -4.2850
2 -4.2849
3 -0.1777
4 -0.1776
5 -0.1774
6 -0.1774
===========================================================
...showing only one spin component.
What do I need to do to override the "Startparameters" as dumped by OUTCAR?
Startparameter for this run:
NWRITE = 2 write-flag & timer
PREC = accura medium, high low
ISTART = 0 job : 0-new 1-cont 2-samecut
ICHARG = 2 charge: 1-file 2-atom 10-const
ISPIN = 1 spin polarized calculation? <=======This particular line!
LNONCOLLINEAR = T non collinear calculations
LSORBIT = T spin-orbit coupling
==========================================================
I haven't found a way around the problem. Your assistance will be greatly apprecciated.
Thank you in advance,
S.
I am new with the code so I might be seeing this as a bug that is not really present. I'm using VASP 4.6.28, compiled with pgi fortran and running an a xeon cluster with cent-Os (free fedora)
I have made a calculation on bulk gold both with ISPIN=1 and 2. ISPIN=2 gives me results for two spin channels, as expected.
Then I turn on the SPIN-ORBIT coupling. And the output dumps only one eigenchannel. My INCAR file follows:
===========================================================
Bulk gold test calculation
ISPIN=2
SYSTEM = AuBulk
NWRITE = 2
PREC=Accurate
IALGO=48
GGA=91
LSORBIT = .TRUE.
SAXIS = 0 0 1 ! direction of the magnetic field
NBANDS = 50
NUPDOWN = 1
===========================================================
Standard output is just normal:
[salva@tbird AuBulk]$ ./vasp
running on 1 nodes
distr: one band on 1 nodes, 1 groups
vasp.4.6.28 25Jul05 complex
POSCAR found : 1 types and 1 ions
LDA part: xc-table for Ceperly-Alder, standard interpolation
generate k-points for: 12 12 12
|| W W A A R R N N II N N GGGG !!! || || Full k-point grid generated || Inversion symmetry is not applied || | -----------------------------------------------------------------------------
POSCAR, INCAR and KPOINTS ok, starting setup
WARNING: wrap around errors must be expected
FFT: planning ... 1
reading WAVECAR
entering main loop
N E dE d eps ncg rms rms(c)
.
.
.
RMM: 18 -0.350209045025E+01 0.16626E-04 -0.10425E-06 32337 0.238E-03 1 F= -.35020905E+01 E0= -.35020912E+01 d E =0.224021E-05 mag= 0.0003 0.0003 0.0003
writing wavefunctions
FORTRAN STOP
===========================================================
These are the first entries of the EIGENVAL file:
1 1 1 1
0.1825866E+02 0.2955706E-09 0.2955706E-09 0.2955706E-09 0.5000000E-15
1.0000000000000000E-004
CAR
AuBulk
11 364 50
0.3742002E-15 0.3742002E-15 0.3742002E-15 0.5787037E-03
1 -4.2850
2 -4.2849
3 -0.1777
4 -0.1776
5 -0.1774
6 -0.1774
===========================================================
...showing only one spin component.
What do I need to do to override the "Startparameters" as dumped by OUTCAR?
Startparameter for this run:
NWRITE = 2 write-flag & timer
PREC = accura medium, high low
ISTART = 0 job : 0-new 1-cont 2-samecut
ICHARG = 2 charge: 1-file 2-atom 10-const
ISPIN = 1 spin polarized calculation? <=======This particular line!
LNONCOLLINEAR = T non collinear calculations
LSORBIT = T spin-orbit coupling
==========================================================
I haven't found a way around the problem. Your assistance will be greatly apprecciated.
Thank you in advance,
S.