<div>thanks for your reply ,I changed the lambda=5.0 and used <br> occupations='tetrahedra' but the magnetization still did not converge to the initial value ,but the total energy converged.</div>
<div>I do know what is the problem. <br></div>
<div class="gmail_quote">2009/3/11 Gabriele Sclauzero <span dir="ltr"><<a href="mailto:sclauzer@sissa.it">sclauzer@sissa.it</a>></span><br>
<blockquote class="gmail_quote" style="PADDING-LEFT: 1ex; MARGIN: 0px 0px 0px 0.8ex; BORDER-LEFT: #ccc 1px solid">
<div class="im"><br><br>yumin qian wrote:<br>> Dear PWSCF users<br>> I am doing calculate on the antiferromagnetic ground state of<br>> LiNiO2, I want to plot the energy as a function of fixed spin moment of<br>
> Ni atom.<br>> the following is my in put file for pwscf<br>><br>> &system<br>> ibrav = 0, nat=16 , ntyp= 4,<br>> ecutwfc = 35.0, ecutrho = 200.0,<br>> occupations='smearing', smearing='gaussian', degauss=0.02<br>
<br></div>First of all a word of caution: this value of degauss may not be low enough to converge<br>the magnetization (even if it is more than enough to converge the total energy). Also<br>ecutrho may not be high enough to converge total energy, since you're using US PPs.<br>
<div>
<div></div>
<div class="h5"><br>> nspin=2,constrained_magnetization='atomic',<br>> starting_magnetization(1)= 0.0,<br>> starting_magnetization(2)= 1.8,<br>> starting_magnetization(3)=-1.8,<br>> starting_magnetization(4)= 0.0,<br>
> lda_plus_u=.true.<br>> Hubbard_U(1)=0.0<br>> Hubbard_U(2)=3.0<br>> Hubbard_U(3)=3.0<br>> Hubbard_U(4)=0.0<br>> /<br>> &electrons<br>> diagonalization='david',<br>> conv_thr = 1.0e-6,<br>
> mixing_beta = 0.5,<br>> /<br>> ATOMIC_SPECIES<br>> Li *6.914* Li.pbe-n-mt_bw.UPF<br>> Ni1 *58.6934* Ni.pbe-n-van.UPF<br>> Ni2 *58.6934* Ni.pbe-n-van.UPF<br>> O 15.9994 O.pbe-van_bm.UPF<br>
><br>> and the output file gives<br>><br>> ==============================================================================<br>> atom number 5 relative position : -0.1652 0.2839 0.5083<br>> charge : 2.373353<br>
> magnetization : 1.544423<br>> magnetization/charge: 0.650735<br>> constrained moment : 1.800000<br>><br>> ==============================================================================<br>
> the output file gives a magnetic moment of 1.5444 ,how should I do to<br>> make the calculation converge to the starting_magnetization 1.8 ?<br><br></div></div>You can increase the lambda parameter which premultiplies the penalty energy added in<br>
v_of_rho. In the limit your constraint gets satisfied, the energy should not depend on<br>which lambda you chose, but you may want to verify this.<br>The default should be 1. I think is not enough to force the magnetization on the atom.<br>
<br>However, I think that the magnetization of a single atom inside a bulk material or in<br>general when it forms covalent bonds is not a well defined quantity, so I would not care<br>too much of the precise number I get for the magnetization, since it may well change if I<br>
change which portion of charge is attributed to which atom (remember that you're using a<br>plane waves code, and PWs do not know anything about atomic orbitals a priori...)<br><br>HTH<br><br>GS<br>
<div class="im"><br>><br>> I change the starting_magnetization(2) and starting_magnetization(3)<br>> from 1.8 to 0.8 , then I will get the E(m) curve. how should I do ,<br>> does anyone give<br>> a suggestion ?<br>
> --<br>> Sincerely Y. M. Qian<br>> Lab.of Condensed Matter Theory and Materials Computation<br>> Institute of Physics<br>> Chinese Academy of Sciences<br>> Tel: + 8610 8264 9147<br></div>> <a href="mailto:E-Mail%3Ayuminqian@gmail.com">E-Mail:yuminqian@gmail.com</a> <mailto:<a href="mailto:E-Mail%253Ayuminqian@gmail.com">E-Mail%3Ayuminqian@gmail.com</a>><br>
<div class="im">> P.O.Box 603 Beijing 100190<br>> China<br>><br>><br></div>> ------------------------------------------------------------------------<br>><br>> _______________________________________________<br>
> Pw_forum mailing list<br>> <a href="mailto:Pw_forum@pwscf.org">Pw_forum@pwscf.org</a><br>> <a href="http://www.democritos.it/mailman/listinfo/pw_forum" target="_blank">http://www.democritos.it/mailman/listinfo/pw_forum</a><br>
<br>--<br><br><br>o ------------------------------------------------ o<br>| Gabriele Sclauzero, PhD Student |<br>| c/o: SISSA & CNR-INFM Democritos, |<br>| via Beirut 2-4, 34014 Trieste (Italy) |<br>
| email: <a href="mailto:sclauzer@sissa.it">sclauzer@sissa.it</a> |<br>| phone: +39 040 3787 511 |<br>| skype: gurlonotturno |<br>o ------------------------------------------------ o<br>
_______________________________________________<br>Pw_forum mailing list<br><a href="mailto:Pw_forum@pwscf.org">Pw_forum@pwscf.org</a><br><a href="http://www.democritos.it/mailman/listinfo/pw_forum" target="_blank">http://www.democritos.it/mailman/listinfo/pw_forum</a><br>
</blockquote></div><br><br clear="all"><br>-- <br>Sincerely Y. M. Qian <br>Lab.of Condensed Matter Theory and Materials Computation <br>Institute of Physics <br>Chinese Academy of Sciences<br>Tel: + 8610 8264 9147<br><a href="mailto:E-Mail%3Ayuminqian@gmail.com">E-Mail:yuminqian@gmail.com</a> <br>
P.O.Box 603 Beijing 100190<br>China <br><br>