<html xmlns:v="urn:schemas-microsoft-com:vml" xmlns:o="urn:schemas-microsoft-com:office:office" xmlns:w="urn:schemas-microsoft-com:office:word" xmlns:st2="urn:schemas:contacts" xmlns:st1="urn:schemas-microsoft-com:office:smarttags" xmlns="http://www.w3.org/TR/REC-html40">

<head>
<meta http-equiv=Content-Type content="text/html; charset=us-ascii">
<meta name=Generator content="Microsoft Word 11 (filtered medium)">
<!--[if !mso]>
<style>
v\:* {behavior:url(#default#VML);}
o\:* {behavior:url(#default#VML);}
w\:* {behavior:url(#default#VML);}
.shape {behavior:url(#default#VML);}
</style>
<![endif]-->
<title>Re: [Wien] PHONON calculation</title>
<o:SmartTagType namespaceuri="urn:schemas:contacts" name="title"/>
<o:SmartTagType namespaceuri="urn:schemas-microsoft-com:office:smarttags"
 name="City"/>
<o:SmartTagType namespaceuri="urn:schemas:contacts" name="middlename"/>
<o:SmartTagType namespaceuri="urn:schemas-microsoft-com:office:smarttags"
 name="place"/>
<o:SmartTagType namespaceuri="urn:schemas:contacts" name="Sn"/>
<o:SmartTagType namespaceuri="urn:schemas:contacts" name="GivenName"/>
<o:SmartTagType namespaceuri="urn:schemas-microsoft-com:office:smarttags"
 name="PersonName"/>
<!--[if !mso]>
<style>
st2\:*{behavior:url(#default#ieooui) }st1\:*{behavior:url(#default#ieooui) }
</style>
<![endif]-->
<style>
<!--
 /* Font Definitions */
 @font-face
        {font-family:Tahoma;
        panose-1:2 11 6 4 3 5 4 4 2 4;}
 /* Style Definitions */
 p.MsoNormal, li.MsoNormal, div.MsoNormal
        {margin:0in;
        margin-bottom:.0001pt;
        font-size:12.0pt;
        font-family:"Times New Roman";}
a:link, span.MsoHyperlink
        {color:blue;
        text-decoration:underline;}
a:visited, span.MsoHyperlinkFollowed
        {color:blue;
        text-decoration:underline;}
p
        {mso-margin-top-alt:auto;
        margin-right:0in;
        mso-margin-bottom-alt:auto;
        margin-left:0in;
        font-size:12.0pt;
        font-family:"Times New Roman";}
p.avgcert, li.avgcert, div.avgcert
        {mso-margin-top-alt:auto;
        margin-right:0in;
        mso-margin-bottom-alt:auto;
        margin-left:0in;
        font-size:12.0pt;
        font-family:"Times New Roman";}
span.EmailStyle19
        {mso-style-type:personal-reply;
        font-family:"Times New Roman";
        color:windowtext;
        font-weight:normal;
        font-style:normal;
        text-decoration:none none;}
@page Section1
        {size:8.5in 11.0in;
        margin:1.0in 1.25in 1.0in 1.25in;}
div.Section1
        {page:Section1;}
-->
</style>

</head>

<body lang=EN-US link=blue vlink=blue>

<div class=Section1>

<p class=MsoNormal><font size=2 face="Times New Roman"><span style='font-size:
11.0pt'>Please try to use 2x2x1 or 1x1x2 type of supercells and see if the same
thing happens. You should also go by what <st1:PersonName w:st="on"><st2:title
 w:st="on">Prof.</st2:title> <st2:Sn w:st="on">Blaha</st2:Sn></st1:PersonName>
has suggested. <o:p></o:p></span></font></p>

<p class=MsoNormal><font size=2 face="Times New Roman"><span style='font-size:
11.0pt'><o:p>&nbsp;</o:p></span></font></p>

<p class=MsoNormal><font size=2 face="Times New Roman"><span style='font-size:
11.0pt'>SG <o:p></o:p></span></font></p>

<p class=MsoNormal><font size=2 face="Times New Roman"><span style='font-size:
11.0pt'><o:p>&nbsp;</o:p></span></font></p>

<p class=MsoNormal><font size=2 face="Times New Roman"><span style='font-size:
11.0pt'><o:p>&nbsp;</o:p></span></font></p>

<div>

<div class=MsoNormal align=center style='text-align:center'><font size=3
face="Times New Roman"><span style='font-size:12.0pt'>

<hr size=2 width="100%" align=center tabindex=-1>

</span></font></div>

<p class=MsoNormal><b><font size=2 face=Tahoma><span style='font-size:10.0pt;
font-family:Tahoma;font-weight:bold'>From:</span></font></b><font size=2
face=Tahoma><span style='font-size:10.0pt;font-family:Tahoma'>
wien-bounces@zeus.theochem.tuwien.ac.at
[mailto:wien-bounces@zeus.theochem.tuwien.ac.at] <b><span style='font-weight:
bold'>On Behalf Of </span></b>puday@iitk.ac.in<br>
<b><span style='font-weight:bold'>Sent:</span></b> Friday, January 14, 2011
6:29 PM<br>
<b><span style='font-weight:bold'>To:</span></b> <st1:PersonName w:st="on"><st2:GivenName
 w:st="on">A</st2:GivenName> <st2:middlename w:st="on">Mailing list</st2:middlename>
 <st2:middlename w:st="on">for</st2:middlename> <st2:middlename w:st="on">WIEN2k</st2:middlename>
 <st2:Sn w:st="on">users</st2:Sn></st1:PersonName><br>
<b><span style='font-weight:bold'>Subject:</span></b> Re: [<st2:Sn w:st="on">Wien</st2:Sn>]
PHONON calculation</span></font><o:p></o:p></p>

</div>

<p class=MsoNormal><font size=3 face="Times New Roman"><span style='font-size:
12.0pt'><o:p>&nbsp;</o:p></span></font></p>

<p><font size=2 face="Times New Roman"><span style='font-size:10.0pt'>Dear <st2:GivenName
w:st="on">SUDDHASATTWA</st2:GivenName><br>
I am sending this mail to you only.<br>
I am doing PHONON calculation using WIEN2K.I am making super cell like<br>
2x2x2 and making *.D45 file using phonon.Now while doing init_phono_lapw I<br>
put nn=2. It ask for the new structure file created.I don't accept the new<br>
structure file.Next comes the space group option as given bellow<br>
<br>
-----&gt; DO YOU WANT TO USE THE NEW case_1.struct_nn file (y/n)<br>
n<br>
<br>
warning: !!! Struct file is not consistent with space group found.<br>
Number and name of space group: 99 (P 4 m m)<br>
warning: !!! Number of inequivalent atoms has changed.<br>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; !!! Old value=
16&nbsp;&nbsp;&nbsp;&nbsp; New value= 8<br>
warning: !!! Struct file is not consistent with space group found.<br>
&nbsp;-----&gt;&nbsp; check in&nbsp; case_1.outputsgroup&nbsp; for proper
symmetry, compare<br>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; with your struct file and
later with&nbsp; case_1.outputs<br>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; sgroup has also produced a new
struct file based on your old one.<br>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; If you see warnings above,
consider to use the newly generated<br>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; struct file, which you can
view (edit) now.<br>
&nbsp;-----&gt;&nbsp; continue with symmetry or edit case_1.struct_sgroup ?
(c/e)<br>
<br>
If I continue with this then I see the space group have changed for the<br>
super cell than the original structure<br>
<br>
Is it right way to do that?<br>
<br>
<br>
In fact I have done phonon calculation for other 3-4 examples but all of<br>
them giving negative force and negative dispersion.<br>
<br>
<br>
Have you done PHONON calculation successfully for any structure using WIEN2K?<br>
Please help me.<br>
<br>
Regards,<br>
<st1:PersonName w:st="on"><st2:GivenName w:st="on">Uday</st2:GivenName> <st2:middlename
 w:st="on">Bhanu</st2:middlename> <st2:Sn w:st="on">Paramanik</st2:Sn></st1:PersonName><br>
Dept.of Physics<br>
IIT <st1:City w:st="on"><st1:place w:st="on">Kanpur</st1:place></st1:City><br>
<br>
<br>
&gt; Do not accept the new structure file.<br>
&gt;<br>
&gt; It is also mentioned in the UG that<br>
&gt;<br>
&gt; &quot;...sgroup might change the unit cell and then the collection of forces<br>
&gt; into<br>
&gt; the original supercell is not possible&quot;<br>
&gt;<br>
&gt;<br>
&gt;<br>
&gt;<br>
&gt;<br>
&gt;<br>
&gt;<br>
&gt;&nbsp;&nbsp; _____<br>
&gt;<br>
&gt; From: wien-bounces@zeus.theochem.tuwien.ac.at<br>
&gt; [<a href="mailto:wien-bounces@zeus.theochem.tuwien.ac.at">mailto:wien-bounces@zeus.theochem.tuwien.ac.at</a>]
On Behalf Of<br>
&gt; puday@iitk.ac.in<br>
&gt; Sent: Tuesday, January 11, 2011 12:19 PM<br>
&gt; To: Wien@zeus.theochem.tuwien.ac.at<br>
&gt; Subject: [<st2:Sn w:st="on">Wien</st2:Sn>] PHONON calculation<br>
&gt;<br>
&gt;<br>
&gt;<br>
&gt; Dear <st1:PersonName w:st="on"><st2:title w:st="on">Prof</st2:title> <st2:Sn
 w:st="on">Blaha</st2:Sn></st1:PersonName> and <st2:Sn w:st="on">WIEN</st2:Sn>
user,<br>
&gt; I am doing PHONON calculation of &quot;AlNi&quot; as given example in the
manual.I<br>
&gt; am making the supercell 2x2x2(16 atoms) of this cubic symmetry using<br>
&gt; phonon and I get a file *.d45 and this file I am using the input structure<br>
&gt; file in WIEN2K.<br>
&gt; Now I am running init_phonon_lapw<br>
&gt; I am putting nn=2<br>
&gt; and this warning is coming.....<br>
&gt;<br>
&gt;&nbsp;&nbsp;&nbsp;&nbsp; ATOM 16&nbsp;
Ni&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; ATOM&nbsp; 8&nbsp; <st2:GivenName
w:st="on">Al</st2:GivenName><br>
&gt;&nbsp; RMT( 16)=2.28000 AND RMT(&nbsp; 8)=2.02000<br>
&gt;&nbsp; SUMS TO 4.30000&nbsp; LT.&nbsp; NN-DIST= 4.70099<br>
&gt;&nbsp; WARNING: Mult not equal. PLEASE CHECK outputnn-file<br>
&gt;&nbsp; WARNING: Mult not equal. PLEASE CHECK outputnn-file<br>
&gt;&nbsp; WARNING: ityp not equal. PLEASE CHECK outputnn-file<br>
&gt;&nbsp; WARNING: ityp not equal. PLEASE CHECK outputnn-file<br>
&gt;&nbsp; WARNING: Mult not equal. PLEASE CHECK outputnn-file<br>
&gt;&nbsp; WARNING: ityp not equal. PLEASE CHECK outputnn-file<br>
&gt;&nbsp; WARNING: Mult not equal. PLEASE CHECK outputnn-file<br>
&gt;&nbsp; WARNING: ityp not equal. PLEASE CHECK outputnn-file<br>
&gt;&nbsp; WARNING: Mult not equal. PLEASE CHECK outputnn-file<br>
&gt;&nbsp; WARNING: ityp not equal. PLEASE CHECK outputnn-file<br>
&gt;&nbsp; WARNING: Mult not equal. PLEASE CHECK outputnn-file<br>
&gt;&nbsp; WARNING: ityp not equal. PLEASE CHECK outputnn-file<br>
&gt;&nbsp; WARNING: Mult not equal. PLEASE CHECK outputnn-file<br>
&gt;&nbsp; WARNING: ityp not equal. PLEASE CHECK outputnn-file<br>
&gt;&nbsp; WARNING: Mult not equal. PLEASE CHECK outputnn-file<br>
&gt;&nbsp; WARNING: ityp not equal. PLEASE CHECK outputnn-file<br>
&gt;&nbsp; WARNING: Mult not equal. PLEASE CHECK outputnn-file<br>
&gt;&nbsp; WARNING: ityp not equal. PLEASE CHECK outputnn-file<br>
&gt;&nbsp; WARNING: Mult not equal. PLEASE CHECK outputnn-file<br>
&gt;&nbsp; WARNING: ityp not equal. PLEASE CHECK outputnn-file<br>
&gt;&nbsp; WARNING: Mult not equal. PLEASE CHECK outputnn-file<br>
&gt;&nbsp; WARNING: ityp not equal. PLEASE CHECK outputnn-file<br>
&gt;&nbsp; WARNING: Mult not equal. PLEASE CHECK outputnn-file<br>
&gt;&nbsp; WARNING: ityp not equal. PLEASE CHECK outputnn-file<br>
&gt; FORTRAN STOP NN created a new CASE.STRUCT_NN FILE<br>
&gt; 0.021u 0.001s 0:01.22 1.6%&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0+0k 0+0io 0pf+0w<br>
&gt; -----&gt; check in&nbsp; case_1.outputnn&nbsp; for overlapping spheres,<br>
&gt;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; coordination and nearest
neighbor distances<br>
&gt; -----&gt; DO YOU WANT TO USE THE NEW case_1.struct_nn file (y/n)<br>
&gt;<br>
&gt;<br>
&gt; If I accept this new structure file and new space group 99 (P 4 m m) then<br>
&gt; the program is running but while calculating Hellmann-Feynman forces and<br>
&gt; dispersion curve I am not getting same as given in the example. Forces and<br>
&gt; Dispersion curve going into negative!!<br>
&gt;<br>
&gt;<br>
&gt; One thing is written in the WIEN2K user manual I am putting that line<br>
&gt; here.<br>
&gt;<br>
&gt; Please use mainly nn to reduce equivalent atoms. sgroup might change the<br>
&gt; unitcell and than the collection of forces into the original supercell is<br>
&gt; not possible (or quite difficult).<br>
&gt;<br>
&gt; I don't understand these lines.<br>
&gt; Any help will be highly appreciated.<br>
&gt;<br>
&gt;<br>
&gt; Regards,<br>
&gt; Uday<br>
&gt; Dept.of Physics<br>
&gt; Indian Institute of Technology, <st1:City w:st="on"><st1:place w:st="on">Kanpur</st1:place></st1:City><br>
&gt; _______________________________________________<br>
&gt; <st2:Sn w:st="on">Wien</st2:Sn> mailing list<br>
&gt; Wien@zeus.theochem.tuwien.ac.at<br>
&gt; <a href="http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien">http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien</a><br>
&gt;<br>
&gt;&nbsp;&nbsp; _____<br>
&gt;<br>
&gt; No virus found in this message.<br>
&gt; Checked by AVG - www.avg.com<br>
&gt; Version: 10.0.1170 / Virus Database: 1435/3372 - Release Date: 01/10/11<br>
&gt;<br>
&gt; _______________________________________________<br>
&gt; <st2:Sn w:st="on">Wien</st2:Sn> mailing list<br>
&gt; Wien@zeus.theochem.tuwien.ac.at<br>
&gt; <a href="http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien">http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien</a><br>
&gt;<br>
<br>
<br>
<br>
_______________________________________________<br>
<st2:Sn w:st="on">Wien</st2:Sn> mailing list<br>
Wien@zeus.theochem.tuwien.ac.at<br>
<a href="http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien">http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien</a></span></font>
<o:p></o:p></p>

<div class=MsoNormal align=center style='text-align:center'><font size=3
face="Times New Roman"><span style='font-size:12.0pt'>

<hr size=1 width="100%" noshade color="#aca899" align=center>

</span></font></div>

<p class=avgcert color="#000000"><font size=3 face="Times New Roman"><span
style='font-size:12.0pt'>No virus found in this message.<br>
Checked by AVG - <a href="http://www.avg.com">www.avg.com</a><br>
Version: 10.0.1170 / Virus Database: 1435/3378 - Release Date: 01/13/11<o:p></o:p></span></font></p>

</div>

</body>

</html>