<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:m="http://schemas.microsoft.com/office/2004/12/omml" xmlns="http://www.w3.org/TR/REC-html40">
<head>
<meta http-equiv="Content-Type" content="text/html; charset=utf-8">
<meta name="Generator" content="Microsoft Word 15 (filtered medium)">
<style><!--
/* Font Definitions */
@font-face
{font-family:Helvetica;
panose-1:2 11 6 4 2 2 2 2 2 4;}
@font-face
{font-family:"Cambria Math";
panose-1:2 4 5 3 5 4 6 3 2 4;}
@font-face
{font-family:Calibri;
panose-1:2 15 5 2 2 2 4 3 2 4;}
/* Style Definitions */
p.MsoNormal, li.MsoNormal, div.MsoNormal
{margin:0cm;
margin-bottom:.0001pt;
font-size:11.0pt;
font-family:"Calibri",sans-serif;}
a:link, span.MsoHyperlink
{mso-style-priority:99;
color:#0563C1;
text-decoration:underline;}
a:visited, span.MsoHyperlinkFollowed
{mso-style-priority:99;
color:#954F72;
text-decoration:underline;}
p.msonormal0, li.msonormal0, div.msonormal0
{mso-style-name:msonormal;
mso-margin-top-alt:auto;
margin-right:0cm;
mso-margin-bottom-alt:auto;
margin-left:0cm;
font-size:11.0pt;
font-family:"Calibri",sans-serif;}
span.E-mailStijl19
{mso-style-type:personal-reply;
font-family:"Calibri",sans-serif;
color:windowtext;}
.MsoChpDefault
{mso-style-type:export-only;
font-size:10.0pt;}
@page WordSection1
{size:612.0pt 792.0pt;
margin:70.85pt 70.85pt 70.85pt 70.85pt;}
div.WordSection1
{page:WordSection1;}
--></style><!--[if gte mso 9]><xml>
<o:shapedefaults v:ext="edit" spidmax="1026" />
</xml><![endif]--><!--[if gte mso 9]><xml>
<o:shapelayout v:ext="edit">
<o:idmap v:ext="edit" data="1" />
</o:shapelayout></xml><![endif]-->
</head>
<body lang="NL-BE" link="#0563C1" vlink="#954F72">
<div class="WordSection1">
<p class="MsoNormal"><span lang="EN-US" style="mso-fareast-language:EN-US">TOLF is needed only when you are creating your initial, entirely optimized cell that will be the start for your phonon calculations. It is not used during the force calculation in the
distorted cells.<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US" style="mso-fareast-language:EN-US"><o:p> </o:p></span></p>
<p class="MsoNormal"><span lang="EN-US" style="mso-fareast-language:EN-US">-fc 0.1 is used as a stopping criterium during one scf cycle: the scf cycle will stop when forces do not change any longer by more than 0.1.<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US" style="mso-fareast-language:EN-US"><o:p> </o:p></span></p>
<p class="MsoNormal"><span lang="EN-US" style="mso-fareast-language:EN-US">-fc must be smaller than TOLF, rule of thumb 3-10 times smaller. Imagine the opposite, say -fc 1.0 and TOLF=0.1 : the scf-cycle will stop as soon as forces vary within 1.0 only, but
the position optimization would stop only if forces are below 0.1. If that would happen, it would be by chance only.<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US" style="mso-fareast-language:EN-US"><o:p> </o:p></span></p>
<p class="MsoNormal"><span lang="EN-US" style="mso-fareast-language:EN-US">About the counting in your distorted cells: do you report the total number of atoms, or the inequivalent number of atoms? For more unambiguous counting, it’s easier to count the number
of electrons (NE in case.in2). Does it make more sense that way? If you have 2 times more electrons, the total number of k-points can be divided by 2 (e.g. from -numk 1000 to -numk 500), without having to care about divisions per axis.<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US" style="mso-fareast-language:EN-US"><o:p> </o:p></span></p>
<p class="MsoNormal"><span lang="EN-US" style="mso-fareast-language:EN-US">Stefaan<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US" style="mso-fareast-language:EN-US"><o:p> </o:p></span></p>
<p class="MsoNormal"><a name="_MailEndCompose"><span lang="EN-US" style="mso-fareast-language:EN-US"><o:p> </o:p></span></a></p>
<span style="mso-bookmark:_MailEndCompose"></span>
<div style="border:none;border-left:solid blue 1.5pt;padding:0cm 0cm 0cm 4.0pt">
<div>
<div style="border:none;border-top:solid #E1E1E1 1.0pt;padding:3.0pt 0cm 0cm 0cm">
<p class="MsoNormal"><b><span lang="NL">Van:</span></b><span lang="NL"> Wien [mailto:wien-bounces@zeus.theochem.tuwien.ac.at]
<b>Namens </b>chin Sabsu<br>
<b>Verzonden:</b> vrijdag 26 januari 2018 17:11<br>
<b>Aan:</b> wien@zeus.theochem.tuwien.ac.at<br>
<b>Onderwerp:</b> Re: [Wien] Optimization convergence for phonon calculations.<o:p></o:p></span></p>
</div>
</div>
<p class="MsoNormal"><o:p> </o:p></p>
<div>
<p class="MsoNormal" style="margin-bottom:12.0pt"><span style="font-size:12.0pt;font-family:"Helvetica",sans-serif">I would like to join this discussion and have two queries!<br>
<br>
I finished some scf with - fc 0.1. Previously I was not aware of TOLF but now I see the importance of TOLF factor. Now I see that TOLF is defined in case.inM.<br>
<br>
In UG, it is mentioned that "TOLF: force tolerance, geometry optimization will stop when all forces are below TOLF."<br>
<br>
In my case, TOLF value is 5 while - fc is 0.1 and total forces on the system is "4.598435 mRy/bohr. So according to Peter's reply, this scf calculation is not good for phonon calculations and the forces must be small as much as we can get (0.1mRy/bohr).<br>
<br>
<br>
My question is can we use the same value for TOLF and -fc? i.e. 0.1 .<o:p></o:p></span></p>
<div>
<p class="MsoNormal"><span style="font-size:12.0pt;font-family:"Helvetica",sans-serif">I have an additional question regarding a number of k-points:<o:p></o:p></span></p>
</div>
<div>
<p class="MsoNormal"><span style="font-size:12.0pt;font-family:"Helvetica",sans-serif"><o:p> </o:p></span></p>
</div>
<p class="MsoNormal" style="margin-bottom:12.0pt"><span style="font-size:12.0pt;font-family:"Helvetica",sans-serif">The pristine unit cell (225 SG has 3 non-equivalent atoms). Its P-cell is having 12 non-equivalent atoms.
<o:p></o:p></span></p>
<div>
<p class="MsoNormal"><span style="font-size:12.0pt;font-family:"Helvetica",sans-serif">If I create the distorted structure for Phonopy by 1x1x1 cell then the atoms become 16, 15, 28 and 28 with no change in lattice dimension.<o:p></o:p></span></p>
</div>
<div>
<p class="MsoNormal"><span style="font-size:12.0pt;font-family:"Helvetica",sans-serif"><o:p> </o:p></span></p>
</div>
<div>
<p class="MsoNormal"><span style="font-size:12.0pt;font-family:"Helvetica",sans-serif">When I create the distorted structure for by 2x2x2 cell then the atoms become 96, 72, 184 and 184 nonequivalent atoms in four *struct_00* files and lattice parameters are
doubled which is obviously should be!<o:p></o:p></span></p>
</div>
<div>
<p class="MsoNormal"><span style="font-size:12.0pt;font-family:"Helvetica",sans-serif"><o:p> </o:p></span></p>
</div>
<div>
<p class="MsoNormal"><span style="font-size:12.0pt;font-family:"Helvetica",sans-serif">I know the basic relation between cell dimension and mesh size, nxnxn supercell Kx/n, Ky/n, Kz/n mesh size, but here I stuck and not able to pick to suitable mesh size which
reduces my computational power.<o:p></o:p></span></p>
</div>
<div>
<p class="MsoNormal"><span style="font-size:12.0pt;font-family:"Helvetica",sans-serif"><o:p> </o:p></span></p>
</div>
<div>
<p class="MsoNormal"><span style="font-size:12.0pt;font-family:"Helvetica",sans-serif">In my bulk optimization, I used 12 12 12 mesh (for 3 non-equivalent atoms and total 9 atoms in 225 SG cell).<o:p></o:p></span></p>
</div>
<div>
<p class="MsoNormal"><span style="font-size:12.0pt;font-family:"Helvetica",sans-serif"><o:p> </o:p></span></p>
</div>
<div>
<p class="MsoNormal"><span style="font-size:12.0pt;font-family:"Helvetica",sans-serif">According to your expertise what should be the mesh size for the structure created by 1x1x1 supercell (four distorted structures with 16, 15, 28 and 28 atoms)?<o:p></o:p></span></p>
</div>
<div>
<p class="MsoNormal"><span style="font-size:12.0pt;font-family:"Helvetica",sans-serif"><o:p> </o:p></span></p>
</div>
<div>
<p class="MsoNormal"><span style="font-size:12.0pt;font-family:"Helvetica",sans-serif">What should be the mesh size for the cases (from 2x22x2 supercell) with atoms 96, 72, 184 and 184?<o:p></o:p></span></p>
</div>
<div>
<p class="MsoNormal"><span style="font-size:12.0pt;font-family:"Helvetica",sans-serif"><o:p> </o:p></span></p>
</div>
<div>
<p class="MsoNormal"><span style="font-size:12.0pt;font-family:"Helvetica",sans-serif">Why I asked this question on mesh size is: because with 6x6x6 mesh size and RKmax 6 one scf cycle is not completed even after two hrs on 8 cores (for a single case).<o:p></o:p></span></p>
</div>
<div>
<p class="MsoNormal"><span style="font-size:12.0pt;font-family:"Helvetica",sans-serif"><o:p> </o:p></span></p>
</div>
<div>
<p class="MsoNormal"><span style="font-size:12.0pt;font-family:"Helvetica",sans-serif">Hope, I explained my query nicely!<o:p></o:p></span></p>
</div>
<div>
<p class="MsoNormal"><span style="font-size:12.0pt;font-family:"Helvetica",sans-serif"><o:p> </o:p></span></p>
</div>
<p class="MsoNormal"><span style="font-size:12.0pt;font-family:"Helvetica",sans-serif"><br>
Regards<br>
<br>
Chin S.<br>
<br>
<br>
On Friday 26 January 2018, 8:55:18 PM IST, Peter Blaha <<a href="mailto:pblaha@theochem.tuwien.ac.at">pblaha@theochem.tuwien.ac.at</a>> wrote:
<br>
<br>
<br>
<br>
Yes, of course for phonon calculations you have first to optimize the <br>
structure. And for that, TOLF should be reduced, because for phonons you <br>
should start with a structure which has "zero" (typically 0.1 to 0.01) <br>
forces.<br>
tolf tells the code when to stop moving atoms.<br>
-fc xx tells the code to stop the scf when forces are converged to xx mRy.<br>
<br>
<br>
On 01/26/2018 02:30 PM, Luis Ogando wrote:<br>
> Dear Wien2k community,<br>
> <o:p></o:p></span></p>
</div>
</div>
</div>
</body>
</html>