<div dir="ltr"><div class="gmail_default" style="font-family:monospace,monospace;font-size:small">Dear Ajit,</div><div class="gmail_default" style="font-family:monospace,monospace;font-size:small">please bear in mind that I AM NOT AN EXPERT and I would instead really appreciate comments from more experienced users. I will try to give some advice but take it with a pinch of salt. I would also like to follow up on this by asking: </div><div class="gmail_default" style="font-family:monospace,monospace;font-size:small"><br></div><div class="gmail_default" style="font-family:monospace,monospace;font-size:small">can noncolinear calculations and DFT+U go together now?</div><div class="gmail_default" style="font-family:monospace,monospace;font-size:small"><br></div><div class="gmail_default" style="font-family:monospace,monospace;font-size:small">Referring to your points:</div><div class="gmail_default" style="font-family:monospace,monospace;font-size:small"><br></div><div class="gmail_default" style="font-family:monospace,monospace;font-size:small">1)the AFM initial guess is set up fine. The system might still converge to a FM situation, though. If instead you'd rather force it to stay AFM I'd go for 'tot_magnetization'</div><div class="gmail_default" style="font-family:monospace,monospace;font-size:small">2)yes. You may also need to specify lda_plus_u=.true. beforehand.</div><div class="gmail_default" style="font-family:monospace,monospace;font-size:small">3)I'd say so, and I also think that the 'angle' isn't a fixed value but rather will be optimised during the scf. Instead if you'd prefer the two (or total) magnetizations to be kept fixed along y I'd go for 'constrained_magnetization'</div><div class="gmail_default" style="font-family:monospace,monospace;font-size:small">4)I think that the default 'domag' (whatever it means) is set to false, hence there should be no need for your flag to be there. Also, 'nspin=4' may be redundant considering that you have set 'noncolin=.true.'</div><div class="gmail_default" style="font-family:monospace,monospace;font-size:small"><br></div><div class="gmail_default" style="font-family:monospace,monospace;font-size:small">hope some more hard-core expert will comment on this.</div><div class="gmail_default" style="font-family:monospace,monospace;font-size:small"><br></div><div class="gmail_default" style="font-family:monospace,monospace;font-size:small">good luck (:</div><div class="gmail_default" style="font-family:monospace,monospace;font-size:small">Roberta</div></div><div class="gmail_extra"><br clear="all"><div><div class="gmail_signature"><div dir="ltr"><div><div dir="ltr"><pre><span style="font-family:monospace,monospace"><font size="2">Roberta Pigliapochi
PhD Student | Grey Research Group<br>University of Cambridge</font></span></pre></div></div></div></div></div>
<br><div class="gmail_quote">On 8 September 2015 at 12:48, Ajit Kumar Jena <span dir="ltr"><<a href="mailto:jenapsi@gmail.com" target="_blank">jenapsi@gmail.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><div><div>Dear All,<br></div>            I wish to do a <b>spin-orbit calculation</b> of our system (LiFePO4) which involves iron (Fe) as a magnetic element. Her<i>e </i>I am providing my input variables containing the system name lists and atomic species.<br><br>My questions are the following :<br><br>1) As 
it is mentioned in the input, should I consider two different types of 
spins ( up and down ) to do a anti-ferromagnetic calculation ? <br>2) The way I have incorporated Hubbard 'U' is appropriate or not ?<br>3) If I am giving the inputs in this way, will it take the starting spin quantization axis along Y-direction or not ? <br>4) Is there any problem if i will not include " starting_spin_angle = .true. " ?<br><br>Any kind of help is highly appreciated. I am extremely sorry if my questions are very trivial.<br><br><b><u>Input</u></b><br>&system<br>    ibrav=  8,<br>    celldm(1) = 19.326236,<br>    celldm(2) = 0.587151,<br>    celldm(3) = 0.458766,<br>    nat=  28, <br>    ntyp= 7,<br>!    nspin = 4,<br>    starting_magnetization(2) = 0.7,<br>    starting_magnetization(3) = -0.7,<br>    ecutwfc =30.0,<br>    ecutrho = 250.0,<br>    occupations = 'smearing' ,<br>    degauss = 0.001 ,<br>    smearing = 'gaussian' ,<br>    noncolin = .true.<br>    lspinorb = .true.<br>    starting_spin_angle = .true.<br>    angle1(2) = 90.0<br>    angle2(2) = 90.0<br>    angle1(3) = 90.0<br>    angle2(3) = 90.0<br>    lda_plus_u_kind = 0,<br>    Hubbard_U(2)=3.0,<br>    Hubbard_U(3)=3.0,<br><br>ATOMIC_SPECIES<br> Li1  6.94      Li.pbe-s-van_ak.UPF<br> Fe1  55.845    Fe.rel-pbe-spn-rrkjus_psl.0.2.1.UPF<br> Fe2  55.845    Fe.rel-pbe-spn-rrkjus_psl.0.2.1.UPF<br> P1   30.9737   P.pbe-van_ak.UPF<br> O1   15.999    O.pbe-van_ak.UPF<br> O2   15.999    O.pbe-van_ak.UPF<br> O3   15.999    O.pbe-van_ak.UPF<br><br></div><div>Thanks & Regards,<br></div><div>Ajit,<br></div>IIT Madras, Chennai</div>
<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" rel="noreferrer" target="_blank">http://pwscf.org/mailman/listinfo/pw_forum</a><br></blockquote></div><br></div>