Dear Gabriele and Andrea,<div><br></div><div>Thank you very much. I fixed the starting_magnetization and the k-point sampling. Indeed, it appears that the problem was the insufficient number of k-points. With <span class="Apple-style-span" style="font-family: arial, sans-serif; font-size: 13px; border-collapse: collapse; ">K-POINTS AUTOMATIC (4 4 4 1 1 1) I do get magnetic moment parallel, disregard of their starting orientations.</span></div>
<div><span class="Apple-style-span" style="font-family: arial, sans-serif; font-size: 13px; border-collapse: collapse; "><br></span></div><div><font class="Apple-style-span" face="arial, sans-serif"><span class="Apple-style-span" style="border-collapse: collapse;">But now I end up getting another problem. I am trying to impose a constraint on the magnetization angle (i.e. cosine of angle1, theta) by setting:</span></font></div>
<div><font class="Apple-style-span" face="arial, sans-serif"><span class="Apple-style-span" style="border-collapse: collapse;">constrained_magnetization='atomic direction'</span></font></div><div><font class="Apple-style-span" face="arial, sans-serif"><span class="Apple-style-span" style="border-collapse: collapse;">This doesn't seem to work, unless I start with angle1 (theta) being equal for both iron atoms. When I set them (theta's) different, they do not stay fixed anymore. I attach my input file along with the portion of the final output showing the magnetization angles. I thought that probably I have to constrain the directions separately for each atom, but I did not figure out how to do it. Simply setting </span></font></div>
<div><font class="Apple-style-span" face="arial, sans-serif"><span class="Apple-style-span" style="border-collapse: collapse;">constrained_magnetization(1)='atomic direction'</span></font></div><div><font class="Apple-style-span" face="arial, sans-serif"><span class="Apple-style-span" style="border-collapse: collapse;">constrained_magnetization(2)='atomic direction'</span></font></div>
<div><font class="Apple-style-span" face="arial, sans-serif"><span class="Apple-style-span" style="border-collapse: collapse;">results in error, and I couldn't find in the manual and/or examples any help on setting these constraints separately.</span></font></div>
<div><font class="Apple-style-span" face="arial, sans-serif"><span class="Apple-style-span" style="border-collapse: collapse;">Any help to resolve this issue is mostly appreciated (i.e. I want to constrain different magnetization directions on different atoms). I couldn't find the solution on PW forum archives either.</span></font></div>
<div><font class="Apple-style-span" face="arial, sans-serif"><span class="Apple-style-span" style="border-collapse: collapse; "><br></span></font></div><div><font class="Apple-style-span" face="arial, sans-serif"><span class="Apple-style-span" style="border-collapse: collapse;"><span class="Apple-style-span" style="border-collapse: separate; font-family: arial; ">Gabriele, as for your question, I am using version 4.2. And I include non-collinear scheme for this somewhat simple system mostly to test how non-collinearity works with constrained_magnetization, i.e. with fixed theta. I need this for a bit more complicated project I am currently working on.</span></span></font></div>
<div><br></div><div>Thank you,</div><div><br></div><div>Pavel </div><div><br></div><div><span class="Apple-style-span" style="font-family: arial, sans-serif; font-size: 13px; border-collapse: collapse; "><div>Pavel Lukashev, Ph.D.</div>
<div>Research Associate</div><div>UNL, Department of Physics & Astronomy</div><div>208 Theodore Jorgensen Hall</div><div>855 N. 16th Street</div><div>Lincoln, NE 68588-0299</div><div>Office: (402) 472-2400</div></span></div>