<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:st1="urn:schemas:contacts" xmlns:st2="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] Problems in phonon calculations (3)</title>
<o:SmartTagType namespaceuri="urn:schemas:contacts" name="middlename"/>
<o:SmartTagType namespaceuri="urn:schemas-microsoft-com:office:smarttags"
 name="City"/>
<o:SmartTagType namespaceuri="urn:schemas:contacts" name="title"/>
<o:SmartTagType namespaceuri="urn:schemas-microsoft-com:office:smarttags"
 name="PersonName"/>
<o:SmartTagType namespaceuri="urn:schemas-microsoft-com:office:smarttags"
 name="country-region"/>
<o:SmartTagType namespaceuri="urn:schemas:contacts" name="GivenName"/>
<o:SmartTagType namespaceuri="urn:schemas:contacts" name="Sn"/>
<o:SmartTagType namespaceuri="urn:schemas-microsoft-com:office:smarttags"
 name="place" downloadurl="http://www.5iantlavalamp.com/"/>
<!--[if !mso]>
<style>
st1\:*{behavior:url(#default#ieooui) }st2\:*{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'>Dear <st2:PersonName w:st="on"><st1:title w:st="on">Prof.</st1:title> <st1:Sn
 w:st="on">Blaha</st1:Sn></st2:PersonName>, <o:p></o:p></span></font></p>

<p class=MsoNormal><font size=2 face="Times New Roman"><span style='font-size:
11.0pt'>Thank you very much for the detailed reply <o:p></o:p></span></font></p>

<p class=MsoNormal><font size=2 face="Times New Roman"><span style='font-size:
11.0pt'>I had taken a 48 atom supercell with Pnma as space group in the phonon
calculations<o:p></o:p></span></font></p>

<p class=MsoNormal><font size=2 face="Times New Roman"><span style='font-size:
11.0pt'>I had used a about 1400 k&#8211;points in the IBZ with rkmax 7.00 <o:p></o:p></span></font></p>

<p class=MsoNormal><font size=2 face="Times New Roman"><span style='font-size:
11.0pt'>However, I had used rkmax of 8.00 when I optimized the lattice parameters.
<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'>After it took more than 40 hours to complete lapw1, I used a 12 atom
supercell and completed case_1. I need some help here. I copy here the
case_1.scf<o:p></o:p></span></font></p>

<p class=MsoNormal><font size=2 face="Times New Roman"><span style='font-size:
11.0pt'>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; TOTAL FORCE IN mRy/a.u. = |F|&nbsp;&nbsp;&nbsp;&nbsp;
Fx&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
Fy&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; Fz<o:p></o:p></span></font></p>

<p class=MsoNormal><font size=2 face="Times New Roman"><span style='font-size:
11.0pt'>&nbsp;with/without FOR in case.in2<o:p></o:p></span></font></p>

<p class=MsoNormal><font size=2 face="Times New Roman"><span style='font-size:
11.0pt'>:FOR001:&nbsp;&nbsp; 1.ATOM&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
7.217&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; -7.180&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
0.732&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0.000 pa<o:p></o:p></span></font></p>

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

<p class=MsoNormal><font size=2 face="Times New Roman"><span style='font-size:
11.0pt'>:FOR002:&nbsp;&nbsp; 2.ATOM&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
11.249&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; -7.074&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
-8.746&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0.000 pa<o:p></o:p></span></font></p>

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

<p class=MsoNormal><font size=2 face="Times New Roman"><span style='font-size:
11.0pt'>:FOR003:&nbsp;&nbsp; 3.ATOM&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
7.417&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 7.416&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
0.109&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0.000 pa<o:p></o:p></span></font></p>

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

<p class=MsoNormal><font size=2 face="Times New Roman"><span style='font-size:
11.0pt'>:FOR004:&nbsp;&nbsp; 4.ATOM&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
6.958&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 6.957&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
0.120&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0.000 pa<o:p></o:p></span></font></p>

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

<p class=MsoNormal><font size=2 face="Times New Roman"><span style='font-size:
11.0pt'>:FOR005:&nbsp;&nbsp; 5.ATOM&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
14.967&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 14.967&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
0.080&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0.000 pa<o:p></o:p></span></font></p>

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

<p class=MsoNormal><font size=2 face="Times New Roman"><span style='font-size:
11.0pt'>:FOR006:&nbsp;&nbsp; 6.ATOM&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
14.973&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 14.959&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
0.654&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0.000 pa<o:p></o:p></span></font></p>

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

<p class=MsoNormal><font size=2 face="Times New Roman"><span style='font-size:
11.0pt'>:FOR007:&nbsp;&nbsp; 7.ATOM&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
14.113&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; -14.028&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
1.543&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0.000 pa<o:p></o:p></span></font></p>

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

<p class=MsoNormal><font size=2 face="Times New Roman"><span style='font-size:
11.0pt'>:FOR008:&nbsp;&nbsp; 8.ATOM&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
15.810&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; -15.749&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
1.386&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0.000 pa<o:p></o:p></span></font></p>

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

<p class=MsoNormal><font size=2 face="Times New Roman"><span style='font-size:
11.0pt'>:FOR009:&nbsp;&nbsp; 9.ATOM&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
83.030&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 82.549&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
8.927&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0.000 pa<o:p></o:p></span></font></p>

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

<p class=MsoNormal><font size=2 face="Times New Roman"><span style='font-size:
11.0pt'>:FOR010:&nbsp; 10.ATOM&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 95.213&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
94.803&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 8.826&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
0.000 pa<o:p></o:p></span></font></p>

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

<p class=MsoNormal><font size=2 face="Times New Roman"><span style='font-size:
11.0pt'>:FOR011:&nbsp; 11.ATOM&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 88.802&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
-88.782&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 1.906&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
0.000 pa<o:p></o:p></span></font></p>

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

<p class=MsoNormal><font size=2 face="Times New Roman"><span style='font-size:
11.0pt'>:FOR012:&nbsp; 12.ATOM&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 88.672&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
-88.672&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; -0.236&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
0.000 pa<o:p></o:p></span></font></p>

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

<p class=MsoNormal><font size=2 face="Times New Roman"><span style='font-size:
11.0pt'>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; TOTAL FORCE WITH RESPECT TO THE
GLOBAL COORDINATE SYSTEM:<o:p></o:p></span></font></p>

<p class=MsoNormal><font size=2 face="Times New Roman"><span style='font-size:
11.0pt'>:FGL001:&nbsp;&nbsp; 1.ATOM&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
0.732000000&nbsp;&nbsp;&nbsp;&nbsp; 0.000000000&nbsp;&nbsp;&nbsp; -7.180000000
pa<o:p></o:p></span></font></p>

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

<p class=MsoNormal><font size=2 face="Times New Roman"><span style='font-size:
11.0pt'>:FGL002:&nbsp;&nbsp; 2.ATOM&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
-8.746000000&nbsp;&nbsp;&nbsp;&nbsp; 0.000000000&nbsp;&nbsp;&nbsp; -7.074000000
pa<o:p></o:p></span></font></p>

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

<p class=MsoNormal><font size=2 face="Times New Roman"><span style='font-size:
11.0pt'>:FGL003:&nbsp;&nbsp; 3.ATOM&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
0.109000000&nbsp;&nbsp;&nbsp;&nbsp; 0.000000000&nbsp;&nbsp;&nbsp;&nbsp; 7.416000000
pa<o:p></o:p></span></font></p>

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

<p class=MsoNormal><font size=2 face="Times New Roman"><span style='font-size:
11.0pt'>:FGL004:&nbsp;&nbsp; 4.ATOM&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
0.120000000&nbsp;&nbsp;&nbsp;&nbsp; 0.000000000&nbsp;&nbsp;&nbsp;&nbsp; 6.957000000
pa<o:p></o:p></span></font></p>

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

<p class=MsoNormal><font size=2 face="Times New Roman"><span style='font-size:
11.0pt'>:FGL005:&nbsp;&nbsp; 5.ATOM&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
0.080000000&nbsp;&nbsp;&nbsp;&nbsp; 0.000000000&nbsp;&nbsp;&nbsp; 14.967000000
pa<o:p></o:p></span></font></p>

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

<p class=MsoNormal><font size=2 face="Times New Roman"><span style='font-size:
11.0pt'>:FGL006:&nbsp;&nbsp; 6.ATOM&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
0.654000000&nbsp;&nbsp;&nbsp;&nbsp; 0.000000000&nbsp;&nbsp;&nbsp; 14.959000000
pa<o:p></o:p></span></font></p>

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

<p class=MsoNormal><font size=2 face="Times New Roman"><span style='font-size:
11.0pt'>:FGL007:&nbsp;&nbsp; 7.ATOM&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
1.543000000&nbsp;&nbsp;&nbsp;&nbsp; 0.000000000&nbsp;&nbsp; -14.028000000 pa<o:p></o:p></span></font></p>

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

<p class=MsoNormal><font size=2 face="Times New Roman"><span style='font-size:
11.0pt'>:FGL008:&nbsp;&nbsp; 8.ATOM&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
1.386000000&nbsp;&nbsp;&nbsp;&nbsp; 0.000000000&nbsp;&nbsp; -15.749000000 pa<o:p></o:p></span></font></p>

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

<p class=MsoNormal><font size=2 face="Times New Roman"><span style='font-size:
11.0pt'>:FGL009:&nbsp;&nbsp; 9.ATOM&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
8.927000000&nbsp;&nbsp;&nbsp;&nbsp; 0.000000000&nbsp;&nbsp;&nbsp; 82.549000000
pa<o:p></o:p></span></font></p>

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

<p class=MsoNormal><font size=2 face="Times New Roman"><span style='font-size:
11.0pt'>:FGL010:&nbsp; 10.ATOM&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
8.826000000&nbsp;&nbsp;&nbsp;&nbsp; 0.000000000&nbsp;&nbsp;&nbsp; 94.803000000
pa<o:p></o:p></span></font></p>

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

<p class=MsoNormal><font size=2 face="Times New Roman"><span style='font-size:
11.0pt'>:FGL011:&nbsp; 11.ATOM&nbsp; &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;1.906000000&nbsp;&nbsp;&nbsp;&nbsp;
0.000000000&nbsp;&nbsp; -88.782000000 pa<o:p></o:p></span></font></p>

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

<p class=MsoNormal><font size=2 face="Times New Roman"><span style='font-size:
11.0pt'>:FGL012:&nbsp; 12.ATOM&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
-0.236000000&nbsp;&nbsp;&nbsp;&nbsp; 0.000000000&nbsp;&nbsp; -88.672000000 pa<o:p></o:p></span></font></p>

<p class=MsoNormal><font size=2 face="Times New Roman"><span style='font-size:
11.0pt'>rtial forces<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'>May I kindly request you to suggest on the forces. Is it ok or do I
have to change &#8211;fc 0.1 to &#8211;fc 0.01? <o:p></o:p></span></font></p>

<p class=MsoNormal><font size=2 face="Times New Roman"><span style='font-size:
11.0pt'>Thank you <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'>Suddhasattwa<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><st2:PersonName w:st="on"><st1:GivenName w:st="on">Peter</st1:GivenName>
 <st1:Sn w:st="on">Blaha</st1:Sn></st2:PersonName><br>
<b><span style='font-weight:bold'>Sent:</span></b> Saturday, January 08, 2011
7:18 PM<br>
<b><span style='font-weight:bold'>To:</span></b> <st2:PersonName w:st="on"><st1:GivenName
 w:st="on">A</st1:GivenName> <st1:middlename w:st="on">Mailing list</st1:middlename>
 <st1:middlename w:st="on">for</st1:middlename> <st1:middlename w:st="on">WIEN2k</st1:middlename>
 <st1:Sn w:st="on">users</st1:Sn></st2:PersonName><br>
<b><span style='font-weight:bold'>Subject:</span></b> Re: [<st1:Sn w:st="on">Wien</st1:Sn>]
Problems in phonon calculations (3)</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'>You need
to understand what you are doing:<br>
<br>
I guess you submit your job in a directory&nbsp;&nbsp; AB2_phonon.<br>
<br>
The job creates a .machines file in THIS directory (look at it with some<br>
editor and verify that it looks as expected (with 24 lines)).<br>
<br>
Now understand the run_phonon job file.<br>
<br>
Inside the loop over all 18 cases it has a:<br>
<br>
cd case_$i&nbsp;&nbsp; and inside this directory it makes a run_lapw ...<br>
<br>
So change into case_1 and examine the .machines file.<br>
<br>
It will of course NOT show the 24 lines, but is some dummy-default .machines
file with just 2 lines.<br>
<br>
so before the&nbsp;&nbsp;&nbsp;&nbsp; cd case_$i&nbsp;&nbsp;&nbsp;&nbsp; insert
a line:<br>
<br>
cp .machines case_$i<br>
<br>
into the run_phonon job.<br>
<br>
<br>
PS: If your job is still running, you can also do this manually on the
commandline (for all case_1 ... case_18)<br>
and the proper .machines file will be used for all new iterations.<br>
<br>
PPS: If I look into your cpu-time (more than 40 h !! for lapw1) I doubt that
you want to wait<br>
until this job has finished (even on 24 cpus). How large is your supercell -
this determines how you can<br>
reduce the k-mesh!<br>
<br>
PPS: You should learn on small examples how to make meaningful calculations.
You must examine for<br>
YOUR system (nobody can help you here) and your properties, what a good RKMAX
or k-mesh is.<br>
In other words: caclculate phonons with a &quot;low RKMAX and few
k-points&quot;, and increase the parameters<br>
until the frequencies did not change significantly anymore.<br>
<br>
<br>
<br>
Am 08.01.2011 07:42, schrieb <st2:PersonName w:st="on"><st1:GivenName w:st="on">Ghosh</st1:GivenName>
 <st1:Sn w:st="on">SUDDHASATTWA</st1:Sn></st2:PersonName>:<br>
&gt; Dear Wien2k users, and <st2:PersonName w:st="on"><st1:title w:st="on">Prof.</st1:title>
 <st1:Sn w:st="on">Blaha</st1:Sn></st2:PersonName>,<br>
&gt;<br>
&gt; Can you please help me out in the following querie. As suggested earlier,
I had initialized my phonon calculations with 48 atoms per lattice.<br>
&gt;<br>
&gt; I created and initialized case_1<br>
&gt;<br>
&gt; My phonon job file is<br>
&gt;<br>
&gt; #!/bin/csh -f<br>
&gt;<br>
&gt; #<br>
&gt;<br>
&gt; set file=AB2_phonon<br>
&gt;<br>
&gt; #<br>
&gt;<br>
&gt; foreach i ( \<br>
&gt;<br>
&gt; 1 \<br>
&gt; 17 \<br>
&gt;<br>
&gt; 18 \<br>
&gt;<br>
&gt; )<br>
&gt;<br>
&gt; cd case_$i<br>
&gt;<br>
&gt; echo running case_$i<br>
&gt;<br>
&gt; #<br>
&gt;<br>
&gt; # select other options if necessary<br>
&gt;<br>
&gt; run_lapw -I -i 200 -cc 0.0001 -in1ef -p -fc 0.1<br>
&gt;<br>
&gt; #<br>
&gt;<br>
&gt; # select other save-name if necessary<br>
&gt;<br>
&gt; save_lapw case_${i}_gga_rkm7.00_3000k<br>
&gt;<br>
&gt; cd ..<br>
&gt;<br>
&gt; end<br>
&gt;<br>
&gt; I submitted the job by qsub &#8211;pe kpoint 24 kpoint.sh<br>
&gt;<br>
&gt; My machine file is<br>
&gt;<br>
&gt; 1:ibnx70<br>
&gt;<br>
&gt; 1:ibnx70<br>
&gt;<br>
&gt; 1:ibnx70<br>
&gt;<br>
&gt; 1:ibnx81<br>
&gt;<br>
&gt; 1:ibnx81<br>
&gt;<br>
&gt; 1:ibnx81<br>
&gt;<br>
&gt; 1:ibnx81<br>
&gt;<br>
&gt; 1:ibnx81<br>
&gt;<br>
&gt; 1:ibnx79<br>
&gt;<br>
&gt; 1:ibnx79<br>
&gt;<br>
&gt; 1:ibnx79<br>
&gt;<br>
&gt; 1:ibnx65<br>
&gt;<br>
&gt; 1:ibnx65<br>
&gt;<br>
&gt; 1:ibnx60<br>
&gt;<br>
&gt; 1:ibnx60<br>
&gt;<br>
&gt; 1:ibnx60<br>
&gt;<br>
&gt; 1:ibnx60<br>
&gt;<br>
&gt; 1:ibnx60<br>
&gt;<br>
&gt; 1:ibnx60<br>
&gt;<br>
&gt; 1:ibnx60<br>
&gt;<br>
&gt; 1:ibnx60<br>
&gt;<br>
&gt; 1:ibnx96<br>
&gt;<br>
&gt; 1:ibnx96<br>
&gt;<br>
&gt; 1:ibnx96<br>
&gt;<br>
&gt; granularity:1<br>
&gt;<br>
&gt; extrafine:1<br>
&gt;<br>
&gt; For the last three days, the queuing shows the job is running<br>
&gt;<br>
&gt; Surprisingly, there is no &#8220;running&#8221; mode in w2web.<br>
&gt;<br>
&gt; I checked the case_1 directory<br>
&gt;<br>
&gt; case_1.dayfile shows (surprisingly)<br>
&gt;<br>
&gt; Calculating case_1 in
/group5/cg/sghosh/WIEN2k/lapw/phonon/UZr2_phonon/case_1<br>
&gt;<br>
&gt; on nx70.igcar.gov.in with PID 10356<br>
&gt;<br>
&gt; start (Thu Jan 6 19:11:06 IST 2011) with lapw0 (-in1ef/99 to go)<br>
&gt;<br>
&gt; cycle 1 (Thu Jan 6 19:11:06 IST 2011) (-in1ef/99 to go)<br>
&gt;<br>
&gt;&gt;&nbsp; lapw0 -p (19:11:06) starting parallel lapw0 at Thu Jan 6
19:11:06 IST 20<br>
&gt;<br>
&gt; 11<br>
&gt;<br>
&gt; -------- .machine0 : processors<br>
&gt;<br>
&gt; running lapw0 in single mode<br>
&gt;<br>
&gt; 177.410u 5.691s 3:11.36 95.6% 0+0k 0+0io 17pf+0w<br>
&gt;<br>
&gt;&gt;&nbsp; lapw1 -c -p (19:14:17) starting parallel lapw1 at Thu Jan 6
19:14:1<br>
&gt;<br>
&gt; 8 IST 2011<br>
&gt;<br>
&gt; -&gt; starting parallel LAPW1 jobs at Thu Jan 6 19:14:19 IST 2011<br>
&gt;<br>
&gt; running LAPW1 in parallel mode (using .machines)<br>
&gt;<br>
&gt; 2 number_of_parallel_jobs<br>
&gt;<br>
&gt; localhost(672) 90833.210u 359.642s 1+01:28:03.65 99.46% 0+0k 0+0io 0pf<br>
&gt;<br>
&gt; +0w<br>
&gt;<br>
&gt; 105493.144u 437.929s 1+13:49:58.03 77.78% 0+0k 0+0io 0pf+0w<br>
&gt;<br>
&gt; 120959.238u 433.725s 1+09:50:20.14 99.65% 0+0k 0+0io 0pf+0w<br>
&gt;<br>
&gt; localhost(672) 90579.555u 334.655s 1+01:22:25.95 99.53% 0+0k 0+0io 0pf<br>
&gt;<br>
&gt; +0w<br>
&gt;<br>
&gt; 105170.999u 500.369s 1+13:44:24.32 77.78% 0+0k 0+0io 0pf+0w<br>
&gt;<br>
&gt; 120363.259u 426.543s 1+09:40:41.44 99.63% 0+0k 0+0io 0pf+0w<br>
&gt;<br>
&gt; Summary of lapw1para:<br>
&gt;<br>
&gt; localhost k=1344 user=181413 wallclock=170<br>
&gt;<br>
&gt; 73.491u 124.896s 33:50:25.33 0.1% 0+0k 0+0io 0pf+0w<br>
&gt;<br>
&gt;&gt;&nbsp; lapw2 -c -p (05:04:43) running LAPW2 in parallel mode<br>
&gt;<br>
&gt; [1] 4250<br>
&gt;<br>
&gt; [2] 4292<br>
&gt;<br>
&gt; That is 2 number of parallel jobs. I thought it should be 24<br>
&gt;<br>
&gt; Can you please why is this happening?<br>
&gt;<br>
&gt; Thank you<br>
&gt;<br>
&gt; Suddhasattwa<br>
&gt;<br>
&gt; SUDDHASATTWAGHOSH<br>
&gt;<br>
&gt; Scientific Officer (D)<br>
&gt;<br>
&gt; Pyrochemical Process Studies Section<br>
&gt;<br>
&gt; Fuel Chemistry Division<br>
&gt;<br>
&gt; Chemistry Group<br>
&gt;<br>
&gt; IndiraGandhi Centre for Atomic Research<br>
&gt;<br>
&gt; Kalpakkam<br>
&gt;<br>
&gt; Tamilnadu<br>
&gt;<br>
&gt; 603102<br>
&gt;<br>
&gt; <st2:country-region w:st="on"><st2:place w:st="on">India</st2:place></st2:country-region><br>
&gt;<br>
&gt; Phone: 91-44-27480500 (Ext:24283)<br>
&gt;<br>
&gt; Fax: 91-44-27480065<br>
&gt;<br>
&gt;<br>
&gt;<br>
&gt; _______________________________________________<br>
&gt; <st1:Sn w:st="on">Wien</st1: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>
<br>
--<br>
-----------------------------------------<br>
<st2:PersonName w:st="on"><st1:GivenName w:st="on">Peter</st1:GivenName> <st1:Sn
 w:st="on">Blaha</st1:Sn></st2:PersonName><br>
Inst. Materials Chemistry, TU Vienna<br>
Getreidemarkt 9, A-1060 <st2:place w:st="on"><st2:City w:st="on">Vienna</st2:City>,
 <st2:country-region w:st="on">Austria</st2:country-region></st2:place><br>
Tel: +43-1-5880115671<br>
Fax: +43-1-5880115698<br>
email: pblaha@theochem.tuwien.ac.at<br>
-----------------------------------------<br>
_______________________________________________<br>
<st1:Sn w:st="on">Wien</st1: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/3366 - Release Date: 01/07/11<o:p></o:p></span></font></p>

</div>

</body>

</html>