[Wien] supercell calc / impurities in V
wiener at arcscluster.caltech.edu
wiener at arcscluster.caltech.edu
Wed Sep 1 02:38:03 CEST 2004
Dear Wien users,
I am trying to calculate a 16 atom 2x2x2 bcc supercell of V
with a central transition metal impurity. When introducing a
Pd or Pt impurity atom in the supercell, I expect an expansion of
the lattice constant. However, my calculations give a
contraction of about -3% by volume as a result.
Another surprising result is a very atomic-like 4d (resp 5d)
band for the Pd impurity (resp Pt) in the DOS, which shows up
as a very sharp peak in the middle of the DOS.
The calculations seem well converged otherwise
(DIS in the range of 1e-3 to 1e-4) and I have no ghost
band warnings. I use 1e3 k-points and the
following case.in1 file (constructed by following the example
of S. Cottenier).
If someone has insights into what is happening and how to
go about it, help would be appreciated.
Thanks!
Olivier.
======================================================
Here is the case.in1 file. The first atom is Pd in this case,
following ones are respectively V1nn, V2nn, V3nn and V5nn.
I am using a core level cutoff energy of -8.0 Ry.
WFFIL (WFPRI, SUPWF)
7.00 10 4 (R-MT*K-MAX; MAX L IN WF, V-NMT
0.45 6 0 (GLOBAL E-PARAMETER WITH n OTHER CHOICES, global
APW/LAPW)
0 0.40 0.000 CONT 0
0 -6.55 0.001 STOP 0
1 0.50 0.000 CONT 0
1 -4.30 0.001 STOP 0
2 0.37 0.000 CONT 1
2 1.50 0.000 CONT 1
0.45 5 0 (GLOBAL E-PARAMETER WITH n OTHER CHOICES, global
APW/LAPW)
0 0.45 0.000 CONT 0
0 -4.89 0.005 STOP 0
1 0.45 0.000 CONT 0
1 -2.94 0.005 STOP 0
2 0.50 0.000 CONT 1
0.45 5 0 (GLOBAL E-PARAMETER WITH n OTHER CHOICES, global
APW/LAPW)
0 0.45 0.000 CONT 0
0 -4.89 0.005 STOP 0
1 0.45 0.000 CONT 0
1 -2.94 0.005 STOP 0
2 0.50 0.000 CONT 1
0.45 5 0 (GLOBAL E-PARAMETER WITH n OTHER CHOICES, global
APW/LAPW)
0 0.45 0.000 CONT 0
0 -4.89 0.005 STOP 0
1 0.45 0.000 CONT 0
1 -2.94 0.005 STOP 0
2 0.50 0.000 CONT 1
0.45 5 0 (GLOBAL E-PARAMETER WITH n OTHER CHOICES, global
APW/LAPW)
0 0.45 0.000 CONT 0
0 -4.89 0.005 STOP 0
1 0.45 0.000 CONT 0
1 -2.94 0.005 STOP 0
2 0.50 0.000 CONT 1
K-VECTORS FROM UNIT:4 -8.0 2.5 emin/emax window
More information about the Wien
mailing list