[Wien] Deformation of Diamond

Hui Wang hwang at imr.ac.cn
Sun May 9 09:42:57 CEST 2010


Dear wien2k master:
    I am running wien2k09.2 on Suse linux OS, with ifort 11 and mkl 10, inter mpi. 
    Here, there are 3 problem confusing me:
    (1)I wanna calculate the bulk modulus of Diamond. Diamond have cubic sysmetry, and there are 8 atoms in unit cell, 2 atoms in primitive cell. My question is: I wanna calculate the unit cell with 8 atoms, I don't wanna calculate the primitive cell with 2 atoms, because "x sgroup" will automatically change to primitive cell with 2 atoms if I set up a unit cell with 8 atoms, so how can I calculate the unit cell with 8 atoms ?
 
    (2)I wanna calculate the elastic constant such as C11 of Diamond, so I wanna deform the unit cell with 8 atoms in x axes, for instance: lattice constant is a=b=c=3.56 Ang, I wanna calculate the unit cell with 8 atoms when a=3.54, a=3.55, a=3.56, a=3.57, a=3.58 (b=c=3.56 all the time). When I calculate a=3.54, the sgroup will find a new spacegroup "141 (I 41/a m d) [origin choice 2]", and "x sgroup" find:
=====================================   
Bravais lattice: Tetragonal body centred 
  a             b            c
 4.73936977    4.73936977   6.72248100
     alpha         beta         gamma
 90.00000000    90.00000000   90.00000000
=====================================
I think this result is too ridiculus, far away from my requirment. If I don't accept the "x sgroup" result, there are errors when "x symetry" ends. So what shoud I do when I wanna deform Diamond along x axes ?
 
    (3)I used to use vasp to calculate bulk modulus of Diamond, when I use Murnaghan fit formula or Mosrse fit formula to get the bulk modulus, the energy unit is eV and lattice constant is Angstron. Here my questions is: wien2k's energy unit is Ry and lattice constant unit is Bohr, so do I need to convert Ry to eV and Angstron to Bohr when I use Murnaghan fit and Morse fit ? 
 
    I am looking forward to your reply, any help will be appreicated:)
 
Yours sincerely:
Hui 

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://zeus.theochem.tuwien.ac.at/pipermail/wien/attachments/20100509/daecb0d2/attachment.htm>


More information about the Wien mailing list