[Wien] PHONON calculation

Ghosh SUDDHASATTWA ssghosh at igcar.gov.in
Tue Jan 11 08:21:50 CET 2011


Do not accept the new structure file. 

It is also mentioned in the UG that

"...sgroup might change the unit cell and then the collection of forces into
the original supercell is not possible"

 

 

 

  _____  

From: wien-bounces at zeus.theochem.tuwien.ac.at
[mailto:wien-bounces at zeus.theochem.tuwien.ac.at] On Behalf Of
puday at iitk.ac.in
Sent: Tuesday, January 11, 2011 12:19 PM
To: Wien at zeus.theochem.tuwien.ac.at
Subject: [Wien] PHONON calculation

 

Dear Prof Blaha and WIEN user,
I am doing PHONON calculation of "AlNi" as given example in the manual.I
am making the supercell 2x2x2(16 atoms) of this cubic symmetry using
phonon and I get a file *.d45 and this file I am using the input structure
file in WIEN2K.
Now I am running init_phonon_lapw
I am putting nn=2
and this warning is coming.....

    ATOM 16  Ni         ATOM  8  Al
 RMT( 16)=2.28000 AND RMT(  8)=2.02000
 SUMS TO 4.30000  LT.  NN-DIST= 4.70099
 WARNING: Mult not equal. PLEASE CHECK outputnn-file
 WARNING: Mult not equal. PLEASE CHECK outputnn-file
 WARNING: ityp not equal. PLEASE CHECK outputnn-file
 WARNING: ityp not equal. PLEASE CHECK outputnn-file
 WARNING: Mult not equal. PLEASE CHECK outputnn-file
 WARNING: ityp not equal. PLEASE CHECK outputnn-file
 WARNING: Mult not equal. PLEASE CHECK outputnn-file
 WARNING: ityp not equal. PLEASE CHECK outputnn-file
 WARNING: Mult not equal. PLEASE CHECK outputnn-file
 WARNING: ityp not equal. PLEASE CHECK outputnn-file
 WARNING: Mult not equal. PLEASE CHECK outputnn-file
 WARNING: ityp not equal. PLEASE CHECK outputnn-file
 WARNING: Mult not equal. PLEASE CHECK outputnn-file
 WARNING: ityp not equal. PLEASE CHECK outputnn-file
 WARNING: Mult not equal. PLEASE CHECK outputnn-file
 WARNING: ityp not equal. PLEASE CHECK outputnn-file
 WARNING: Mult not equal. PLEASE CHECK outputnn-file
 WARNING: ityp not equal. PLEASE CHECK outputnn-file
 WARNING: Mult not equal. PLEASE CHECK outputnn-file
 WARNING: ityp not equal. PLEASE CHECK outputnn-file
 WARNING: Mult not equal. PLEASE CHECK outputnn-file
 WARNING: ityp not equal. PLEASE CHECK outputnn-file
 WARNING: Mult not equal. PLEASE CHECK outputnn-file
 WARNING: ityp not equal. PLEASE CHECK outputnn-file
FORTRAN STOP NN created a new CASE.STRUCT_NN FILE
0.021u 0.001s 0:01.22 1.6%      0+0k 0+0io 0pf+0w
-----> check in  case_1.outputnn  for overlapping spheres,
       coordination and nearest neighbor distances
-----> DO YOU WANT TO USE THE NEW case_1.struct_nn file (y/n)


If I accept this new structure file and new space group 99 (P 4 m m) then
the program is running but while calculating Hellmann-Feynman forces and
dispersion curve I am not getting same as given in the example. Forces and
Dispersion curve going into negative!!


One thing is written in the WIEN2K user manual I am putting that line here.

Please use mainly nn to reduce equivalent atoms. sgroup might change the
unitcell and than the collection of forces into the original supercell is
not possible (or quite difficult).

I don't understand these lines.
Any help will be highly appreciated.


Regards,
Uday
Dept.of Physics
Indian Institute of Technology, Kanpur
_______________________________________________
Wien mailing list
Wien at zeus.theochem.tuwien.ac.at
http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien 

  _____  

No virus found in this message.
Checked by AVG - www.avg.com
Version: 10.0.1170 / Virus Database: 1435/3372 - Release Date: 01/10/11

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


More information about the Wien mailing list