[Wien] Problem in SO calculation

Stefaan Cottenier Stefaan.Cottenier at UGent.be
Thu May 30 08:56:13 CEST 2013



> -in1ef was introduced in WIEN2k_09, but since WIEN2k_10   -in1ef
> does not appear as switch, because this is the default behavior
> anyway.

As I never noticed that, I explicitly tested this. And no, -in1ef does 
not seem to be the default behaviour in my installation:

cat $WIENROOT/VERSION
WIEN2k_12.1 (Release 22/7/2012)

OK -- this is the most recent version.

Now run a small case without any switches:

cp $WIENROOT/example_struct_files/fe4n.struct ./case.struct
init_lapw -b
run_lapw -i 3

This is the top of case.in1 after these 3 iterations -- all 0.30 values 
are still there:

WFFIL  EF=.877157622425   (WFFIL, WFPRI, ENFIL, SUPWF)
   7.00       10    4 (R-MT*K-MAX; MAX L IN WF, V-NMT
   0.30    4  0      (GLOBAL E-PARAMETER WITH n OTHER CHOICES, global
APW/LAPW)
  1    0.30      0.000 CONT 1
  1   -3.97      0.001 STOP 1
  2    0.30      0.005 CONT 1
  0    0.30      0.000 CONT 1

Repeat that with the -in1ef switch:

cp $WIENROOT/example_struct_files/fe4n.struct ./case.struct
init_lapw -b
run_lapw -i 3 -in1ef

Now the 0.2 Ry subtraction did happen:

WFFIL  EF=.877063513400   (WFFIL, WFPRI, ENFIL, SUPWF)
   7.00       10    4 (R-MT*K-MAX; MAX L IN WF, V-NMT
  .68859   4  0      (GLOBAL E-PARAMETER WITH n OTHER CHOICES, global
APW/LAPW)
  1   .68859     0.000 CONT 1
  1   -3.97      0.001 STOP 1
  2   .68859     0.005 CONT 1
  0   .68859     0.000 CONT 1

> Therefore it has also been removed from the UG (but unfortunately not
> from a faq-page).

It looks like it has indeed been removed from the UG, yet not from the 
official distribution of the code...?

Stefaan


More information about the Wien mailing list