<html><head><meta http-equiv="Content-Type" content="text/html charset=utf-8"></head><body style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space;" class="">Thank you for both your emails.<div class=""><br class=""></div><div class="">Answering to the first one, yes, I forgot to write about the magnetic moments but they are zero in both calculations.</div><div class=""><br class=""></div><div class="">Answering to the second one, the forces were a little bit different (less than 1.0 mRy in difference) but I wasn’t expecting such a big difference in the EFG due to it, so I didn’t think that it was important to send them.</div><div class="">My mistake, I was wrong.</div><div class="">They were all in fact below 2 mRy, so the minimization stopped, but it stopped in different positions for each method.</div><div class=""><br class=""></div><div class="">I did a another minimization with both methods until the forces were all below 0.10 mRy (tried first with 1.0 mRy but it was not enough) and in this situation the results are much more similar:</div><div class=""><br class=""></div><div class="">:POS001 = <span style="font-family: Menlo; font-size: 11px;" class="">0.66667 0.33333 -0.00012</span> (no sp)</div><div class="">:POS001 = <span style="font-family: Menlo; font-size: 11px;" class="">0.66667 0.33333 -0.00022 </span>(sp)</div><div class=""><br class=""></div><div class=""><div class="">:POS002 = <span style="font-family: Menlo; font-size: 11px;" class="">0.66667 0.33333 0.37638 </span>(no sp)</div><div class="">:POS002 = <span style="font-family: Menlo; font-size: 11px;" class="">0.66667 0.33333 0.37627 </span>(sp)</div></div><div class=""><br class=""></div><div class="">:EFG001 = <span style="font-family: Menlo; font-size: 11px;" class="">0.41661 </span>(no sp)</div><div class="">:EFG001 = <span style="font-family: Menlo; font-size: 11px;" class="">0.41404 </span>(sp)</div><div class=""><br class=""></div><div class=""><div class="">:EFG002 = <span style="font-family: Menlo; font-size: 11px;" class="">0.10461 </span>(no sp)</div><div class="">:EFG002 = <span style="font-family: Menlo; font-size: 11px;" class="">0.10472 </span>(sp)</div></div><div class=""><br class=""></div><div class="">I was doing these tests to see if I should use spin-polarization or not in calculations of GaN, specially for the calculation of the EFGs.</div><div class="">That’s why I was worried about the different EFGs that I was obtaining.</div><div class="">However, since they give the “same" results after all, am I right to assume that for this material I don’t need to consider spin-polarization?</div><div class=""><br class=""></div><div class="">Thank you both very much for your help,</div><div class="">Marcelo</div><div class=""><br class=""><div><blockquote type="cite" class=""><div class="">On 08/04/2016, at 13:58, Peter Blaha <<a href="mailto:pblaha@theochem.tuwien.ac.at" class="">pblaha@theochem.tuwien.ac.at</a>> wrote:</div><br class="Apple-interchange-newline"><div class=""><span style="font-family: Helvetica; font-size: 12px; font-style: normal; font-variant-caps: normal; font-weight: normal; letter-spacing: normal; orphans: auto; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; widows: auto; word-spacing: 0px; -webkit-text-stroke-width: 0px; float: none; display: inline !important;" class="">What you did not show was that the minimization will give different forces too.</span><br style="font-family: Helvetica; font-size: 12px; font-style: normal; font-variant-caps: normal; font-weight: normal; letter-spacing: normal; orphans: auto; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; widows: auto; word-spacing: 0px; -webkit-text-stroke-width: 0px;" class=""><br style="font-family: Helvetica; font-size: 12px; font-style: normal; font-variant-caps: normal; font-weight: normal; letter-spacing: normal; orphans: auto; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; widows: auto; word-spacing: 0px; -webkit-text-stroke-width: 0px;" class=""><span style="font-family: Helvetica; font-size: 12px; font-style: normal; font-variant-caps: normal; font-weight: normal; letter-spacing: normal; orphans: auto; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; widows: auto; word-spacing: 0px; -webkit-text-stroke-width: 0px; float: none; display: inline !important;" class="">The minimization will stop once the forces are for a couple of times below 2.0 mRy (case.inM) and the movement is small enough.</span><br style="font-family: Helvetica; font-size: 12px; font-style: normal; font-variant-caps: normal; font-weight: normal; letter-spacing: normal; orphans: auto; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; widows: auto; word-spacing: 0px; -webkit-text-stroke-width: 0px;" class=""><br style="font-family: Helvetica; font-size: 12px; font-style: normal; font-variant-caps: normal; font-weight: normal; letter-spacing: normal; orphans: auto; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; widows: auto; word-spacing: 0px; -webkit-text-stroke-width: 0px;" class=""><span style="font-family: Helvetica; font-size: 12px; font-style: normal; font-variant-caps: normal; font-weight: normal; letter-spacing: normal; orphans: auto; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; widows: auto; word-spacing: 0px; -webkit-text-stroke-width: 0px; float: none; display: inline !important;" class="">Even very small differences between the sp and non-sp calculation will lead to a different minimization-path and thus mixer will stop by chance at different positions leading to different EFGs (and E-tot).</span><br style="font-family: Helvetica; font-size: 12px; font-style: normal; font-variant-caps: normal; font-weight: normal; letter-spacing: normal; orphans: auto; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; widows: auto; word-spacing: 0px; -webkit-text-stroke-width: 0px;" class=""><br style="font-family: Helvetica; font-size: 12px; font-style: normal; font-variant-caps: normal; font-weight: normal; letter-spacing: normal; orphans: auto; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; widows: auto; word-spacing: 0px; -webkit-text-stroke-width: 0px;" class=""><span style="font-family: Helvetica; font-size: 12px; font-style: normal; font-variant-caps: normal; font-weight: normal; letter-spacing: normal; orphans: auto; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; widows: auto; word-spacing: 0px; -webkit-text-stroke-width: 0px; float: none; display: inline !important;" class="">In other words: If you would really want to get out the same results in both runs you need to put the cut-off in case.inM very small, so that mixer will optimize the positions until the forces are really small in all cases.</span><br style="font-family: Helvetica; font-size: 12px; font-style: normal; font-variant-caps: normal; font-weight: normal; letter-spacing: normal; orphans: auto; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; widows: auto; word-spacing: 0px; -webkit-text-stroke-width: 0px;" class=""><br style="font-family: Helvetica; font-size: 12px; font-style: normal; font-variant-caps: normal; font-weight: normal; letter-spacing: normal; orphans: auto; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; widows: auto; word-spacing: 0px; -webkit-text-stroke-width: 0px;" class=""><br style="font-family: Helvetica; font-size: 12px; font-style: normal; font-variant-caps: normal; font-weight: normal; letter-spacing: normal; orphans: auto; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; widows: auto; word-spacing: 0px; -webkit-text-stroke-width: 0px;" class=""><span style="font-family: Helvetica; font-size: 12px; font-style: normal; font-variant-caps: normal; font-weight: normal; letter-spacing: normal; orphans: auto; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; widows: auto; word-spacing: 0px; -webkit-text-stroke-width: 0px; float: none; display: inline !important;" class="">On 04/08/2016 12:41 PM, Marcelo Barbosa wrote:</span><br style="font-family: Helvetica; font-size: 12px; font-style: normal; font-variant-caps: normal; font-weight: normal; letter-spacing: normal; orphans: auto; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; widows: auto; word-spacing: 0px; -webkit-text-stroke-width: 0px;" class=""><blockquote type="cite" style="font-family: Helvetica; font-size: 12px; font-style: normal; font-variant-caps: normal; font-weight: normal; letter-spacing: normal; orphans: auto; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; widows: auto; word-spacing: 0px; -webkit-text-stroke-width: 0px;" class="">Dear Sirs,<br class=""><br class="">I did some calculations with GaN (struct file in attachment) to compare the difference in energy and electric field gradient between the use of spin-polarization or not.<br class=""><br class="">The calculations were done using the commands<br class=""><br class="">init_lapw -b -ecut -8 -numk 480 -rkmax 7 (-sp)<br class="">run(sp)_lapw -p -ec 0.00001 -cc 0.0001 -fc 0.5<br class=""><br class="">and the results were the following<br class=""><br class="">:ENE = -7995.60614872 (no spin-polarization)<br class="">:ENE = -7995.60614871 (spin-polarization)<br class=""><br class="">:EFG001 = 0.57937 (no sp)<br class="">:EFG001 = 0.57931 (sp)<br class=""><br class="">so they basically give the same.<br class=""><br class="">Then, I relaxed the structure using<br class=""><br class="">run(sp)_lapw -p -ec 0.00001 -cc 0.0001 -fc 0.5 -min<br class=""><br class="">and the results were<br class=""><br class="">:POS001 = 0.66667 0.33333 0.00002 (no sp)<br class="">:POS001 = 0.66667 0.33333 0.00000 (sp)<br class=""><br class="">:POS002 = 0.66667 0.33333 0.37662 (no sp)<br class="">:POS002 = 0.66667 0.33333 0.37701 (sp)<br class=""><br class="">:ENE = -7995.60615476 (no sp)<br class="">:ENE = -7995.60615035 (sp)<br class=""><br class="">:EFG001 = 0.45038 (no sp)<br class="">:EFG001 = 0.58801 (sp)<br class=""><br class="">so they gave different final positions and very different :EFG, as a consequence.<br class=""><br class="">However, if I use the final structure from the minimization without spin-polarization and do a calculation with spin-polarization (both from the beginning, using the previous commands), I get equal values<br class=""><br class="">:ENE = -7995.60615467 (no sp)<br class="">:ENE = -7995.60615466 (sp)<br class=""><br class="">:EFG001 = 0.45036 (no sp)<br class="">:EFG001 = 0.45034 (sp)<br class=""><br class="">and if I use the final structure from the minimization with spin-polarization and do a calculation without spin-polarization, I also get equal values<br class=""><br class="">:ENE = -7995.60614842 (no sp)<br class="">:ENE = -7995.60614842 (sp)<br class=""><br class="">:EFG001 = 0.58272 (no sp)<br class="">:EFG001 = 0.58271 (sp)<br class=""><br class=""><br class="">This means that the only thing that I am getting different between them is the atomic positions during the minimization, since I always get the same :ENE and :EFG if both calculations are done with the same struct file.<br class=""><br class="">Is there any reason for this to happen?<br class="">Which values for :ENE, :EFG and :POS should I consider the best, the spin-polarized ones or the ones without spin-polarization?<br class=""><br class="">Best regards,<br class="">Marcelo<br class=""><br class=""><br class=""><br class=""><br class=""><br class=""><br class=""><br class="">_______________________________________________<br class="">Wien mailing list<br class=""><a href="mailto:Wien@zeus.theochem.tuwien.ac.at" class="">Wien@zeus.theochem.tuwien.ac.at</a><br class=""><a href="http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien" class="">http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien</a><br class="">SEARCH the MAILING-LIST at: <a href="http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html" class="">http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html</a><br class=""><br class=""></blockquote><br style="font-family: Helvetica; font-size: 12px; font-style: normal; font-variant-caps: normal; font-weight: normal; letter-spacing: normal; orphans: auto; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; widows: auto; word-spacing: 0px; -webkit-text-stroke-width: 0px;" class=""><span style="font-family: Helvetica; font-size: 12px; font-style: normal; font-variant-caps: normal; font-weight: normal; letter-spacing: normal; orphans: auto; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; widows: auto; word-spacing: 0px; -webkit-text-stroke-width: 0px; float: none; display: inline !important;" class="">--<span class="Apple-converted-space"> </span></span><br style="font-family: Helvetica; font-size: 12px; font-style: normal; font-variant-caps: normal; font-weight: normal; letter-spacing: normal; orphans: auto; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; widows: auto; word-spacing: 0px; -webkit-text-stroke-width: 0px;" class=""><br style="font-family: Helvetica; font-size: 12px; font-style: normal; font-variant-caps: normal; font-weight: normal; letter-spacing: normal; orphans: auto; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; widows: auto; word-spacing: 0px; -webkit-text-stroke-width: 0px;" class=""><span style="font-family: Helvetica; font-size: 12px; font-style: normal; font-variant-caps: normal; font-weight: normal; letter-spacing: normal; orphans: auto; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; widows: auto; word-spacing: 0px; -webkit-text-stroke-width: 0px; float: none; display: inline !important;" class=""> P.Blaha</span><br style="font-family: Helvetica; font-size: 12px; font-style: normal; font-variant-caps: normal; font-weight: normal; letter-spacing: normal; orphans: auto; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; widows: auto; word-spacing: 0px; -webkit-text-stroke-width: 0px;" class=""><span style="font-family: Helvetica; font-size: 12px; font-style: normal; font-variant-caps: normal; font-weight: normal; letter-spacing: normal; orphans: auto; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; widows: auto; word-spacing: 0px; -webkit-text-stroke-width: 0px; float: none; display: inline !important;" class="">--------------------------------------------------------------------------</span><br style="font-family: Helvetica; font-size: 12px; font-style: normal; font-variant-caps: normal; font-weight: normal; letter-spacing: normal; orphans: auto; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; widows: auto; word-spacing: 0px; -webkit-text-stroke-width: 0px;" class=""><span style="font-family: Helvetica; font-size: 12px; font-style: normal; font-variant-caps: normal; font-weight: normal; letter-spacing: normal; orphans: auto; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; widows: auto; word-spacing: 0px; -webkit-text-stroke-width: 0px; float: none; display: inline !important;" class="">Peter BLAHA, Inst.f. Materials Chemistry, TU Vienna, A-1060 Vienna</span><br style="font-family: Helvetica; font-size: 12px; font-style: normal; font-variant-caps: normal; font-weight: normal; letter-spacing: normal; orphans: auto; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; widows: auto; word-spacing: 0px; -webkit-text-stroke-width: 0px;" class=""><span style="font-family: Helvetica; font-size: 12px; font-style: normal; font-variant-caps: normal; font-weight: normal; letter-spacing: normal; orphans: auto; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; widows: auto; word-spacing: 0px; -webkit-text-stroke-width: 0px; float: none; display: inline !important;" class="">Phone: +43-1-58801-165300 FAX: +43-1-58801-165982</span><br style="font-family: Helvetica; font-size: 12px; font-style: normal; font-variant-caps: normal; font-weight: normal; letter-spacing: normal; orphans: auto; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; widows: auto; word-spacing: 0px; -webkit-text-stroke-width: 0px;" class=""><span style="font-family: Helvetica; font-size: 12px; font-style: normal; font-variant-caps: normal; font-weight: normal; letter-spacing: normal; orphans: auto; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; widows: auto; word-spacing: 0px; -webkit-text-stroke-width: 0px; float: none; display: inline !important;" class="">Email:<span class="Apple-converted-space"> </span></span><a href="mailto:blaha@theochem.tuwien.ac.at" style="font-family: Helvetica; font-size: 12px; font-style: normal; font-variant-caps: normal; font-weight: normal; letter-spacing: normal; orphans: auto; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; widows: auto; word-spacing: 0px; -webkit-text-stroke-width: 0px;" class="">blaha@theochem.tuwien.ac.at</a><span style="font-family: Helvetica; font-size: 12px; font-style: normal; font-variant-caps: normal; font-weight: normal; letter-spacing: normal; orphans: auto; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; widows: auto; word-spacing: 0px; -webkit-text-stroke-width: 0px; float: none; display: inline !important;" class=""><span class="Apple-converted-space"> </span> WIEN2k:<span class="Apple-converted-space"> </span></span><a href="http://www.wien2k.at/" style="font-family: Helvetica; font-size: 12px; font-style: normal; font-variant-caps: normal; font-weight: normal; letter-spacing: normal; orphans: auto; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; widows: auto; word-spacing: 0px; -webkit-text-stroke-width: 0px;" class="">http://www.wien2k.at</a><br style="font-family: Helvetica; font-size: 12px; font-style: normal; font-variant-caps: normal; font-weight: normal; letter-spacing: normal; orphans: auto; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; widows: auto; word-spacing: 0px; -webkit-text-stroke-width: 0px;" class=""><span style="font-family: Helvetica; font-size: 12px; font-style: normal; font-variant-caps: normal; font-weight: normal; letter-spacing: normal; orphans: auto; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; widows: auto; word-spacing: 0px; -webkit-text-stroke-width: 0px; float: none; display: inline !important;" class="">WWW: </span><a href="http://www.imc.tuwien.ac.at/staff/tc_group_e.php" style="font-family: Helvetica; font-size: 12px; font-style: normal; font-variant-caps: normal; font-weight: normal; letter-spacing: normal; orphans: auto; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; widows: auto; word-spacing: 0px; -webkit-text-stroke-width: 0px;" class="">http://www.imc.tuwien.ac.at/staff/tc_group_e.php</a><br style="font-family: Helvetica; font-size: 12px; font-style: normal; font-variant-caps: normal; font-weight: normal; letter-spacing: normal; orphans: auto; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; widows: auto; word-spacing: 0px; -webkit-text-stroke-width: 0px;" class=""><span style="font-family: Helvetica; font-size: 12px; font-style: normal; font-variant-caps: normal; font-weight: normal; letter-spacing: normal; orphans: auto; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; widows: auto; word-spacing: 0px; -webkit-text-stroke-width: 0px; float: none; display: inline !important;" class="">--------------------------------------------------------------------------</span><br style="font-family: Helvetica; font-size: 12px; font-style: normal; font-variant-caps: normal; font-weight: normal; letter-spacing: normal; orphans: auto; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; widows: auto; word-spacing: 0px; -webkit-text-stroke-width: 0px;" class=""><span style="font-family: Helvetica; font-size: 12px; font-style: normal; font-variant-caps: normal; font-weight: normal; letter-spacing: normal; orphans: auto; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; widows: auto; word-spacing: 0px; -webkit-text-stroke-width: 0px; float: none; display: inline !important;" class="">_______________________________________________</span><br style="font-family: Helvetica; font-size: 12px; font-style: normal; font-variant-caps: normal; font-weight: normal; letter-spacing: normal; orphans: auto; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; widows: auto; word-spacing: 0px; -webkit-text-stroke-width: 0px;" class=""><span style="font-family: Helvetica; font-size: 12px; font-style: normal; font-variant-caps: normal; font-weight: normal; letter-spacing: normal; orphans: auto; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; widows: auto; word-spacing: 0px; -webkit-text-stroke-width: 0px; float: none; display: inline !important;" class="">Wien mailing list</span><br style="font-family: Helvetica; font-size: 12px; font-style: normal; font-variant-caps: normal; font-weight: normal; letter-spacing: normal; orphans: auto; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; widows: auto; word-spacing: 0px; -webkit-text-stroke-width: 0px;" class=""><a href="mailto:Wien@zeus.theochem.tuwien.ac.at" style="font-family: Helvetica; font-size: 12px; font-style: normal; font-variant-caps: normal; font-weight: normal; letter-spacing: normal; orphans: auto; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; widows: auto; word-spacing: 0px; -webkit-text-stroke-width: 0px;" class="">Wien@zeus.theochem.tuwien.ac.at</a><br style="font-family: Helvetica; font-size: 12px; font-style: normal; font-variant-caps: normal; font-weight: normal; letter-spacing: normal; orphans: auto; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; widows: auto; word-spacing: 0px; -webkit-text-stroke-width: 0px;" class=""><a href="http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien" style="font-family: Helvetica; font-size: 12px; font-style: normal; font-variant-caps: normal; font-weight: normal; letter-spacing: normal; orphans: auto; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; widows: auto; word-spacing: 0px; -webkit-text-stroke-width: 0px;" class="">http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien</a><br style="font-family: Helvetica; font-size: 12px; font-style: normal; font-variant-caps: normal; font-weight: normal; letter-spacing: normal; orphans: auto; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; widows: auto; word-spacing: 0px; -webkit-text-stroke-width: 0px;" class=""><span style="font-family: Helvetica; font-size: 12px; font-style: normal; font-variant-caps: normal; font-weight: normal; letter-spacing: normal; orphans: auto; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; widows: auto; word-spacing: 0px; -webkit-text-stroke-width: 0px; float: none; display: inline !important;" class="">SEARCH the MAILING-LIST at: </span><a href="http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html" style="font-family: Helvetica; font-size: 12px; font-style: normal; font-variant-caps: normal; font-weight: normal; letter-spacing: normal; orphans: auto; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; widows: auto; word-spacing: 0px; -webkit-text-stroke-width: 0px;" class="">http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html</a></div></blockquote></div><br class=""></div></body></html>