[Wien] ghost bands
Dr. B.R.Sahu
sahu at matter3.ph.utexas.edu
Fri Aug 1 16:49:31 CEST 2003
Dear W2k users
I have a supercell(of one K + Ti) self consistent calculation with
distance of Ti from first K, d(K-Ti)=8.01 au. In this run, I had set
d(K-Ti)=5.0 au(no core charge leakage) to begin with (ie RMT=2.5 each)
and the calculation converged without any ghost bands(u-dot smaller than
u).I get the expected values.
I wanted to increase the RMT of K and Ti to 4.0 a.u each so that they
almost touch each other ie d(TM-A1)=8.0 au. I started the calculation all
over from scratch. This time the calculation is faster than the previous
run(due to small number of plane waves in the interstitials) but I get the
QTL-B value very large ie
QTL-B VALUE .EQ. 32.30382 !!!!!!
I understand that such a large value for QTL-B is indeed ghost band
signature in the scf calculation.
The parameters used for both runs are SAME:
RMTKMAX= 8.0
number of k-points in IBZ = 120( a dense mesh is not needed for such a
large cell but I believe it is not one of the reason for appearance of
ghost band).
unit cell size = about 20 au.
RMT(TM)=2.5 and RMT(A1)=2.5 in the first case(NO GHOST BANDS).
RMT(TM)=4.0 and RMT(A1)=4.0 in the second case(GHOST BANDS).
Is RMTKMAX=8.0 is large enough to cause ghost bands?
Here is the case.in1 file:
_______________________________________________________________________________
WFFIL (WFPRI, SUPWF)
8.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)
0 -2.29 0.010 CONT 1
0 0.30 0.000 CONT 1
1 -1.07 0.010 CONT 1
1 0.30 0.000 CONT 1
0.30 4 0 (GLOBAL E-PARAMETER WITH n OTHER CHOICES, global
APW/LAPW)
0 -2.29 0.010 CONT 1
0 0.30 0.000 CONT 1
1 -1.07 0.010 CONT 1
1 0.30 0.000 CONT 1
0.30 4 0 (GLOBAL E-PARAMETER WITH n OTHER CHOICES, global
APW/LAPW)
0 -2.29 0.010 CONT 1
0 0.30 0.000 CONT 1
1 -1.07 0.010 CONT 1
1 0.30 0.000 CONT 1
0.30 4 0 (GLOBAL E-PARAMETER WITH n OTHER CHOICES, global
APW/LAPW)
0 -2.29 0.010 CONT 1
0 0.30 0.000 CONT 1
1 -1.07 0.010 CONT 1
1 0.30 0.000 CONT 1
0.30 5 0 (GLOBAL E-PARAMETER WITH n OTHER CHOICES, global
APW/LAPW)
0 0.30 0.000 CONT 1
0 -4.35 0.005 STOP 1
1 -2.58 0.010 CONT 1
1 0.30 0.000 CONT 1
2 0.30 0.010 CONT 1
K-VECTORS FROM UNIT:4 -9.0 2.5 emin/emax window
_______________________________________________________
where I am going wrong?pl. suggest.
best regards
sahu
More information about the Wien
mailing list