[Wien] Modify the supercell-structure file!

Stefaan Cottenier Stefaan.Cottenier at fys.kuleuven.ac.be
Mon Oct 20 15:01:01 CEST 2003


(First of all: your mail was labeled as spam by the spam blocker we 
use here. Maybe some other mailing list readers might have deleted it 
automatically. Your very long added case.struct might have confused 
this spam blocker.)

>           I want to make a 72-atom supercell with 12 unit cells of
> TiO2 where one Ti atom has been replaced by a Cd atom as a
> exercise.Based on the U.G. , we should put the impurity into the
> cell. I do not know how to modify the created supercell-structure 
> file by replacing one atom with an impurity. 

1) increase the number of positions by 1
     
 P   LATTICE,NONEQUIV. ATOMS 25 36_P

2) Take any of the Ti-positions (MULT=n), and make one inequivalent 
position with MULT=n-1 and one with MULT=1 out of it. Change the atome 
of the one with MULT=1 from Ti to Cd. For instance, with this position:

> ATOM   1: X=0.00000000 Y=0.00000000 Z=0.00000000
>           MULT= 2          ISPLIT= 8
> ATOM   1: X=0.25000000 Y=0.25000000 Z=0.16666667
> Ti         NPT=  781  R0=0.00010000 RMT=    1.8000   Z: 22.0
> LOCAL ROT MATRIX:    1.0000000 0.0000000 0.0000000
>                      0.0000000 1.0000000 0.0000000
>                      0.0000000 0.0000000 1.0000000

You make

 ATOM   1: X=0.00000000 Y=0.00000000 Z=0.00000000
           MULT= 1          ISPLIT= 8
 Ti         NPT=  781  R0=0.00010000 RMT=    1.8000   Z: 22.0
 LOCAL ROT MATRIX:    1.0000000 0.0000000 0.0000000
                      0.0000000 1.0000000 0.0000000
                      0.0000000 0.0000000 1.0000000
 ATOM   2: X=0.25000000 Y=0.25000000 Z=0.16666667
           MULT= 1          ISPLIT= 8
 Cd         NPT=  781  R0=0.00010000 RMT=    1.8000   Z: 48.0
 LOCAL ROT MATRIX:    1.0000000 0.0000000 0.0000000
                      0.0000000 1.0000000 0.0000000
                      0.0000000 0.0000000 1.0000000

(mind the "ATOM=  x"-label, which should be a new, not yet present 
number)

3) Choose a suitable Rmt and R0for your impurity.

4) Go through init_lapw, and accept all changes that nn and sgroup 
propose.

5) You might check visually with xcrysden whether you have obtained 
the desired result.

Stefaan





More information about the Wien mailing list