<div dir="ltr"><span style="font-family:arial,sans-serif;font-size:12.727272033691406px">Dear Pang Rui,</span><div><font face="arial, sans-serif"><br></font></div><div><font face="arial, sans-serif">starting_magnetization on atomic type i in a spin polarized calculation (such as nspin=2). Values range of starting_magnetization at -1 (all spin down for the valence electrons of atom) to 1 (all spin-up).<br>
</font></div><div>Following this link:</div><div><br></div><div><a href="http://www.quantum-espresso.org/wp-content/uploads/Doc/INPUT_PW.html#idp113952">http://www.quantum-espresso.org/wp-content/uploads/Doc/INPUT_PW.html#idp113952</a><br>
</div></div><div class="gmail_extra"><br><br><div class="gmail_quote">On Fri, Jul 4, 2014 at 3:36 PM, Gabriele Sclauzero <span dir="ltr"><<a href="mailto:gabriele.sclauzero@mat.ethz.ch" target="_blank">gabriele.sclauzero@mat.ethz.ch</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Dear Pang Rui,<br>
<br>
     It's difficult to judge just from the input. Has the scf converged?<br>
What are the final values of theta for the two atoms? Can you report the<br>
relevant part of the output?<br>
<br>
GS<br>
<div><div class="h5"><br>
On 06/28/2014 11:40 AM, Pang Rui wrote:<br>
> Dear all<br>
> Could anyone gave me some suggestion to do the constraint magnetic moment<br>
> calculation?<br>
> I followed the constraint magnetic moment calculation of the following<br>
> handson with PWSCF.<br>
> <a href="http://www.vasp.at/vasp-workshop/slides/handsonIV.pdf" target="_blank">http://www.vasp.at/vasp-workshop/slides/handsonIV.pdf</a><br>
> VASP uses the same method of PWSCF on constraint MM. However, I found the<br>
> constraint energy is difficult to go zero. In the handson, they use<br>
> lambda=50eV, got a constraint energy of 0.22591E-03eV. I used lambda=150Ry,<br>
> but got a constraint energy=0.35Ry. Even PWSCF is different with VASP on<br>
> some detail, I think it is not reasonable for two codes showing so large<br>
> difference. So could anyone point the unreasonable part of my input file?<br>
> Or give me some suggestion? This puzzled me for quite a long time.<br>
> The following is the input file.<br>
> &system<br>
>      ibrav=  0, nat=  2, ntyp= 2,<br>
>      ecutwfc = 80.0,<br>
>      occupations='smearing',<br>
>      degauss=0.001,<br>
>      smearing='gauss'<br>
>      starting_magnetization(1)=2.5<br>
>      starting_magnetization(2)=2.5<br>
>      angle1(1)=0.0<br>
>      angle1(2)=45.0<br>
>      angle2(1)=0.0<br>
>      angle2(2)=90.0<br>
>      constrained_magnetization="atomic direction"<br>
>      noncolin=.ture.<br>
>      lambda=150.0<br>
>      nosym=.true.<br>
> /<br>
>   &electrons<br>
>      mixing_beta = 0.1<br>
>      electron_maxstep=200<br>
>      startingwfc='file'<br>
>      startingpot='file'<br>
>   /<br>
> ATOMIC_SPECIES<br>
>   Fe1 56 Fe.pbe-spn-kjpaw_psl.0.2.1.UPF<br>
>   Fe2 56 Fe.pbe-spn-kjpaw_psl.0.2.1.UPF<br>
> ATOMIC_POSITIONS angstrom<br>
> Fe1 3.0 0.0 0.0<br>
> Fe2 5.0 0.0 0.0<br>
> K_POINTS automatic<br>
> 1 1 1 0 0 0<br>
> CELL_PARAMETERS angstrom<br>
>      8.0    0.0000000000000000    0.0000000000000000<br>
>      0.0 8.0    0.0000000000000000<br>
>      0.0 0.0 8.0<br>
><br>
> PS1: I found the PWSCF used starting_magnetization as magnetization/charge<br>
> for generating initial wave function,but used it as magnetization in<br>
> constraint MM. So I set 2.5 for them.<br>
> PS2: The startingwfc is of a smaller lambda calculation.<br>
<br>
<br>
--<br>
</div></div>Dr. Gabriele Sclauzero<br>
ETH Zurich<br>
Materials Theory<br>
HIT G 43.2<br>
Wolfgang-Pauli-Str. 27<br>
8093 Zürich, Switzerland<br>
<br>
Phone +41 44 633 94 10<br>
Fax +41 44 633 14 59<br>
<a href="mailto:gabriele.sclauzero@mat.ethz.ch">gabriele.sclauzero@mat.ethz.ch</a><br>
<a href="http://www.theory.mat.ethz.ch/people/postdocs/gsclauze" target="_blank">http://www.theory.mat.ethz.ch/people/postdocs/gsclauze</a><br>
<div class="HOEnZb"><div class="h5"><br>
_______________________________________________<br>
Pw_forum mailing list<br>
<a href="mailto:Pw_forum@pwscf.org">Pw_forum@pwscf.org</a><br>
<a href="http://pwscf.org/mailman/listinfo/pw_forum" target="_blank">http://pwscf.org/mailman/listinfo/pw_forum</a><br>
</div></div></blockquote></div><br><br clear="all"><div><br></div>-- <br><div dir="ltr">Nguyen Van Chuong<div>PhD Student,<br><div>Department of Physics,</div><div>Don State Technical University. </div><div>Rostov on Don, Russia</div>
<div>Group Graphene</div><div>Email: <a href="mailto:chuongnguyen11@gmail.com" target="_blank">chuongnguyen11@gmail.com</a> or <a href="mailto:cpuphys86@gmail.com" target="_blank">c</a><a href="mailto:hnguyen@dstu.edu.ru" target="_blank">hnguyen@d</a><a href="http://stu.edu.ru" target="_blank">stu.edu.ru</a></div>
<div>Website: <a href="http://cpuphys86.webs.com/" target="_blank">http://cpuphys86.webs.com/</a></div><div>Phone mobile: +7 905 45 94 888</div><div><br></div></div></div>
</div>