[Wien] Phonon difficulties

puday at iitk.ac.in puday at iitk.ac.in
Sat Feb 12 04:30:34 CET 2011


Dear Prof. Blaha and wien2k users,

I am doing phonon calculation for the given example in the manual "AlNi" I
have used supercell 2x2x2 and displacement of atom 0.03 Angstrom by
default. It is taking the displacement of the atoms 8 and 9 in the
"AlNi.d45" as

      0.005186   0.000000   0.000000  Al       8
      0.005186   0.000000   0.000000  Ni       9

And see the forces bellow for the atom 8th Al  and 9th Ni.
It is giving huge forces in all direction particularly Fx. And force
constant plot gives negative values.Here is my AlNi.dat file please see
it.


----------------------------------------------------------------------------
! Phonon-file by phonon2wien2k from AlNi.d45
! :RKM : MATRIX SIZE 1365LOs: 152 RKM= 6.97 WEIGHT= 1.00 PGR:
! :KPT : NUMBER OF K-POINTS: 112
   5.784714   0.000000   0.000000
   0.000000   5.784714   0.000000
   0.000000   0.000000   5.784714
         0.2500000  0.2500000  0.2500000      0.0051860   0.0000000  
0.0000000
  1 A   0.7500000 0.7500000 0.7500000       0.2400   -0.0070    0.0020
  2 A   0.2500000 0.7500000 0.7500000       0.1280   -0.0560   -0.0420
  3 A   0.7500000 0.2500000 0.7500000      -0.2400    0.0090   -0.0380
  4 A   0.2500000 0.2500000 0.7500000       0.4160    0.0570   -0.0070
  5 A   0.7500000 0.7500000 0.2500000      -0.2440   -0.0490    0.0000
  6 A   0.2500000 0.7500000 0.2500000       0.4230   -0.0170    0.0480
  7 A   0.7500000 0.2500000 0.2500000       1.7240    0.0420    0.0300
  8 A   0.2500000 0.2500000 0.2500000      -8.2650    0.0130    0.0050
  9 N   0.0000000 0.0000000 0.0000000       0.7010    0.3900    0.3880
 10 N   0.5000000 0.0000000 0.0000000       0.7680   -0.5280   -0.5150
 11 N   0.0000000 0.5000000 0.0000000       0.7000   -0.3810    0.3810
 12 N   0.5000000 0.5000000 0.0000000       0.7610    0.5210   -0.5240
 13 N   0.0000000 0.0000000 0.5000000       0.7000    0.3790   -0.3840
 14 N   0.5000000 0.0000000 0.5000000       0.7590   -0.5280    0.5120
 15 N   0.0000000 0.5000000 0.5000000       0.7110   -0.3790   -0.3790
 16 N   0.5000000 0.5000000 0.5000000       0.7290    0.5260    0.5230
         0.0000000  0.0000000  0.0000000      0.0051860   0.0000000  
0.0000000
  1 A   0.7500000 0.7500000 0.7500000       0.7300    0.4470    0.4620
  2 A   0.2500000 0.7500000 0.7500000       0.7300   -0.6140   -0.5990
  3 A   0.7500000 0.2500000 0.7500000       0.7390   -0.4440    0.4630
  4 A   0.2500000 0.2500000 0.7500000       0.7140    0.6120   -0.5980
  5 A   0.7500000 0.7500000 0.2500000       0.7390    0.4520   -0.4590
  6 A   0.2500000 0.7500000 0.2500000       0.7160   -0.6110    0.5960
  7 A   0.7500000 0.2500000 0.2500000       0.7520   -0.4480   -0.4620
  8 A   0.2500000 0.2500000 0.2500000       0.6970    0.6090    0.5990
  9 N   0.0000000 0.0000000 0.0000000      -5.6470   -0.0200   -0.0170
 10 N   0.5000000 0.0000000 0.0000000      -0.2730    0.0150    0.0150
 11 N   0.0000000 0.5000000 0.0000000      -0.0710    0.0210    0.0390
 12 N   0.5000000 0.5000000 0.0000000       0.4970   -0.0240   -0.0400
 13 N   0.0000000 0.0000000 0.5000000      -0.0650    0.0400    0.0200
 14 N   0.5000000 0.0000000 0.5000000       0.4970   -0.0410   -0.0260
 15 N   0.0000000 0.5000000 0.5000000      -0.7700   -0.0340   -0.0330
 16 N   0.5000000 0.5000000 0.5000000      -0.0010    0.0420    0.0400
----------------------------------------------------------------------------


What should I do to resolve it?Please help me in this regard.Any help will
be very helpfull for me.

reds,
Uday
Research scholar
Dept.of Physics
IIT Kanpur
India


More information about the Wien mailing list