<html><head></head><body style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space; ">Hello German,<div><br><div><div>Il giorno 01/giu/2011, alle ore 23.02, German Samolyuk ha scritto:</div><br class="Apple-interchange-newline"><blockquote type="cite"><div>I withdraw my question. It works with nosym=true.<br>German<br><br><br>---------- Forwarded message ----------<br>From: German Samolyuk <<a href="mailto:samolyuk@gmail.com">samolyuk@gmail.com</a>><br>Date: Wed, May 18, 2011 at 10:25 AM<br>Subject: pwscf: energy convergence of noncollinear with spin-orbit calculations<br>To: PWSCF Forum <<a href="mailto:pw_forum@pwscf.org">pw_forum@pwscf.org</a>><br><br><br>Dear qe users,<br><br>I'm trying to calculate total energy of artificial FCC PtFe compound<br>with magnetic moment (MM) along z- and x-directionsi (MAE).<br><br>In order to supress oscillation of moment direction I used for MM||z calculation<br>input instructions<br><br> constrained_magnetization = 'total direction' ,<br> fixed_magnetization(1) = 0.0,<br> fixed_magnetization(2) = 0.0,<br> fixed_magnetization(3) = 1.0,<br><br></div></blockquote><div><br></div><div>If it is v.>4.2, then a</div><div>Are you sure that this is what you want? </div><div>From the Doc/INPUT_PW:</div><div><br></div><div><div> 'total direction':</div><div> the angle theta of the total magnetization</div><div> with the z axis (theta = fixed_magnetization(3))</div><div> is constrained:</div><div> </div><div> LAMBDA * ( arccos(magnetization(3)/mag_tot) - theta )**2</div><div> </div><div> where mag_tot is the modulus of the total magnetization.</div><div><br></div></div>Then fixed_magnetization(3) = 1.0 does not correspond to the magnetization along the z axis (but you should reach the corresponding constrained solution anyway, even if different from what you actually want).</div><div><br></div><div><div>Can you please specify which version of the code you've been using? </div></div><div>Could you try to compare v.4.1.1 with v.4.3.1 (using the correct value for theta=fixed_magnetization(3))</div><div><div>The constraining potential implemented is different: in the most recent version it should work also for theta=90° (magnetization along x), while in the previous you could get into troubles, since it was defined as:</div><div><div> LAMBDA * ( magnetization(1) - magnetization(3)*tan(theta) )**2</div><div><br></div></div><div> </div><div>Also convergency should be better in the 4.3.1, though I'm not sure of that. You should try to start with local magnetizations on Pt and Fe close to the desired target direction. It would be nice if you perform this test.</div><div><br></div><div>Also, reporting only the total energy oscillations is not that useful. You should at least report the magnitude and angles of the total magnetization while the calculation is going.</div><div><br></div><div>Further advice: from your input, perhaps degauss=0.02 is not enough to get accurate magnetizations. I would use at least 0.01, better 0.005.</div><div><br></div><div><br></div><div>HTH</div><div><br></div><div><br></div><div>GS</div><div><br></div></div><div><br><blockquote type="cite"><div>However I can not reach total energy convergence and it looks like I<br>still have direction oscillation of MM.<br><br>What can be done to resolve this problem?<br><br>Thank you in advance.<br><br>The complete input file for MM||z is and energy accuracy output file<br>are attached.<br><br>German Samolyuk,<br>ORNL, USA<br><span><ptfe.pw.in></span><span><accu.out></span>_______________________________________________<br>Pw_forum mailing list<br><a href="mailto:Pw_forum@pwscf.org">Pw_forum@pwscf.org</a><br>http://www.democritos.it/mailman/listinfo/pw_forum<br></div></blockquote></div><br><div>
<div><span class="Apple-style-span" style="color: rgb(126, 126, 126); font-size: 16px; font-style: italic; "><br class="Apple-interchange-newline"></span><span class="Apple-style-span" style="color: rgb(126, 126, 126); font-size: 16px; font-style: italic; ">§ Gabriele Sclauzero, EPFL SB ITP CSEA</span></div><div><font class="Apple-style-span" color="#7E7E7E"><i> PH H2 462, Station 3, CH-1015 Lausanne</i></font></div>
</div>
<br></div></body></html>