<html xmlns:v="urn:schemas-microsoft-com:vml" xmlns:o="urn:schemas-microsoft-com:office:office" xmlns:w="urn:schemas-microsoft-com:office:word" xmlns:m="http://schemas.microsoft.com/office/2004/12/omml" xmlns="http://www.w3.org/TR/REC-html40">
<head>
<meta http-equiv="Content-Type" content="text/html; charset=utf-8">
<meta name="Generator" content="Microsoft Word 14 (filtered medium)">
<!--[if !mso]><style>v\:* {behavior:url(#default#VML);}
o\:* {behavior:url(#default#VML);}
w\:* {behavior:url(#default#VML);}
.shape {behavior:url(#default#VML);}
</style><![endif]--><style><!--
/* Font Definitions */
@font-face
{font-family:Calibri;
panose-1:2 15 5 2 2 2 4 3 2 4;}
@font-face
{font-family:Tahoma;
panose-1:2 11 6 4 3 5 4 4 2 4;}
@font-face
{font-family:Consolas;
panose-1:2 11 6 9 2 2 4 3 2 4;}
/* Style Definitions */
p.MsoNormal, li.MsoNormal, div.MsoNormal
{margin:0in;
margin-bottom:.0001pt;
font-size:12.0pt;
font-family:"Times New Roman","serif";}
a:link, span.MsoHyperlink
{mso-style-priority:99;
color:blue;
text-decoration:underline;}
a:visited, span.MsoHyperlinkFollowed
{mso-style-priority:99;
color:purple;
text-decoration:underline;}
p
{mso-style-priority:99;
mso-margin-top-alt:auto;
margin-right:0in;
mso-margin-bottom-alt:auto;
margin-left:0in;
font-size:12.0pt;
font-family:"Times New Roman","serif";}
pre
{mso-style-priority:99;
mso-style-link:"HTML Preformatted Char";
margin:0in;
margin-bottom:.0001pt;
font-size:10.0pt;
font-family:"Courier New";}
span.HTMLPreformattedChar
{mso-style-name:"HTML Preformatted Char";
mso-style-priority:99;
mso-style-link:"HTML Preformatted";
font-family:Consolas;}
span.EmailStyle20
{mso-style-type:personal;
font-family:"Calibri","sans-serif";
color:#1F497D;}
span.EmailStyle21
{mso-style-type:personal-reply;
font-family:"Calibri","sans-serif";
color:#1F497D;}
.MsoChpDefault
{mso-style-type:export-only;
font-size:10.0pt;}
@page WordSection1
{size:8.5in 11.0in;
margin:1.0in 1.0in 1.0in 1.0in;}
div.WordSection1
{page:WordSection1;}
--></style><!--[if gte mso 9]><xml>
<o:shapedefaults v:ext="edit" spidmax="1026" />
</xml><![endif]--><!--[if gte mso 9]><xml>
<o:shapelayout v:ext="edit">
<o:idmap v:ext="edit" data="1" />
</o:shapelayout></xml><![endif]-->
</head>
<body lang="EN-US" link="blue" vlink="purple">
<div class="WordSection1">
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D">Dear Juanjo,<o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D"><o:p> </o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D">Thank you for the advice. We did try a DFT+U calculation only on the neodymium atoms, but it was also with the default number of bands. We saw similar results
in terms of total magnetization as the other calculations we ran without the +U correction.
<o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D"><o:p> </o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D">Thanks,<o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D">Michael<o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D"><o:p> </o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D">Graduate Student in Nuclear Engineering<o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D">Department of Mechanical and Aerospace Engineering<o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D">The Ohio State University<o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D"><o:p> </o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D"><o:p> </o:p></span></p>
<div>
<div style="border:none;border-top:solid #B5C4DF 1.0pt;padding:3.0pt 0in 0in 0in">
<p class="MsoNormal"><b><span style="font-size:10.0pt;font-family:"Tahoma","sans-serif"">From:</span></b><span style="font-size:10.0pt;font-family:"Tahoma","sans-serif""> pw_forum-bounces@pwscf.org [mailto:pw_forum-bounces@pwscf.org]
<b>On Behalf Of </b>Juanjo Meléndez<br>
<b>Sent:</b> Monday, June 16, 2014 1:37 PM<br>
<b>To:</b> PWSCF Forum<br>
<b>Subject:</b> Re: [Pw_forum] Nd2Fe14B Magnetization<o:p></o:p></span></p>
</div>
</div>
<p class="MsoNormal"><o:p> </o:p></p>
<div>
<div>
<div>
<p class="MsoNormal"><span style="font-family:"Calibri","sans-serif";color:black">Dear Michael<o:p></o:p></span></p>
</div>
<div>
<p class="MsoNormal"><span style="font-family:"Calibri","sans-serif";color:black"> <o:p></o:p></span></p>
</div>
<div>
<p class="MsoNormal"><span style="font-family:"Calibri","sans-serif";color:black">I think that the problem is the presence of Nd in your lattice. As far as I know, f-electrons are not only tricky to handle within a DFT formalism, but they constitute the border
of its validity. I would not expect a great degree of accuracy in your calculations for any PP you may use.
<o:p></o:p></span></p>
</div>
<div>
<p class="MsoNormal"><span style="font-family:"Calibri","sans-serif";color:black"> <o:p></o:p></span></p>
</div>
<div>
<p class="MsoNormal"><span style="font-family:"Calibri","sans-serif";color:black">Anyway, you could try with DFT+U with the proper correction for correlated f-electrons (and maybe d-electrons from Fe as well). Sometimes it is hard to converge, but it could
help. Anyway, I would not expect any miracle. <o:p></o:p></span></p>
</div>
<div>
<p class="MsoNormal"><span style="font-family:"Calibri","sans-serif";color:black"> <o:p></o:p></span></p>
</div>
<div>
<p class="MsoNormal"><span style="font-family:"Calibri","sans-serif";color:black">BTW. Do not waste your time trying to change ecutrho. You are using a norm-conserving pseudopotential, so that ecutrho is meaningless.
<o:p></o:p></span></p>
</div>
<div>
<p class="MsoNormal"><span style="font-family:"Calibri","sans-serif";color:black"> <o:p></o:p></span></p>
</div>
<div>
<p class="MsoNormal"><span style="font-family:"Calibri","sans-serif";color:black">Good luck<o:p></o:p></span></p>
</div>
<div>
<p class="MsoNormal"><span style="font-family:"Calibri","sans-serif";color:black"> <o:p></o:p></span></p>
</div>
<div>
<p class="MsoNormal"><span style="font-family:"Calibri","sans-serif";color:black">Juanjo<o:p></o:p></span></p>
</div>
<div>
<p class="MsoNormal"><span style="font-family:"Calibri","sans-serif";color:black"> <o:p></o:p></span></p>
</div>
<div>
<p class="MsoNormal"><span style="font-family:"Calibri","sans-serif";color:black">Juan J. Meléndez
<br>
Associate Professor<br>
Department of Physics · University of Extremadura<br>
Avda. de Elvas, s/n 06006 Badajoz (Spain)<br>
Phone: +34 924 28 96 55<br>
Fax: +34 924 28 96 51<br>
Email: melendez@unex.es<br>
Web: http://materiales.unex.es/miembros/personal/jj-melendez/Index.html<o:p></o:p></span></p>
</div>
<div>
<div>
<p class="MsoNormal"><span style="font-size:10.0pt;font-family:"Tahoma","sans-serif";color:black"> <o:p></o:p></span></p>
</div>
<div>
<div>
<p class="MsoNormal" style="background:whitesmoke"><b><span style="font-size:10.0pt;font-family:"Tahoma","sans-serif";color:black">From:</span></b><span style="font-size:10.0pt;font-family:"Tahoma","sans-serif";color:black">
<a href="mailto:kurth.17@buckeyemail.osu.edu" title="kurth.17@buckeyemail.osu.edu">
Kurth, Michael R.</a> <o:p></o:p></span></p>
</div>
<div>
<p class="MsoNormal" style="background:whitesmoke"><b><span style="font-size:10.0pt;font-family:"Tahoma","sans-serif";color:black">Sent:</span></b><span style="font-size:10.0pt;font-family:"Tahoma","sans-serif";color:black"> Monday, June 16, 2014 5:52 PM<o:p></o:p></span></p>
</div>
<div>
<p class="MsoNormal" style="background:whitesmoke"><b><span style="font-size:10.0pt;font-family:"Tahoma","sans-serif";color:black">To:</span></b><span style="font-size:10.0pt;font-family:"Tahoma","sans-serif";color:black">
<a href="mailto:pw_forum@pwscf.org" title="pw_forum@pwscf.org">PWSCF Forum</a> <o:p>
</o:p></span></p>
</div>
<div>
<p class="MsoNormal" style="background:whitesmoke"><b><span style="font-size:10.0pt;font-family:"Tahoma","sans-serif";color:black">Subject:</span></b><span style="font-size:10.0pt;font-family:"Tahoma","sans-serif";color:black"> Re: [Pw_forum] Nd2Fe14B Magnetization<o:p></o:p></span></p>
</div>
</div>
<div>
<p class="MsoNormal"><span style="font-family:"Calibri","sans-serif";color:black"> <o:p></o:p></span></p>
</div>
</div>
<div>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D">Dear Dr. Rebaza,<o:p></o:p></span></p>
<p class="MsoNormal"><span style="color:black"> </span><span style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D"><o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D">Thank you for the reply. We did increase ecutrho to approximately 6 times the ecutwfc value and did not see a notable difference in the calculation. As far
as the value of ecutwfc=260, we arrived at this value after running calculations at ecutwfc=300 and not observing a 1 mRy/atom difference in the total energy. We also have tested PBE and PW pseudopotentials but have not observed a difference in the total magnetization
value. <o:p></o:p></span></p>
<p class="MsoNormal"><span style="color:black"> </span><span style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D"><o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D">Sincerely,<o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D">Michael Kurth<o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D"><o:p> </o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D">Graduate Student in Nuclear Engineering<o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D">Department of Mechanical and Aerospace Engineering<o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D">The Ohio State University<o:p></o:p></span></p>
<p class="MsoNormal"><b><span style="font-size:10.0pt;font-family:"Tahoma","sans-serif";color:black">From:</span></b><span style="font-size:10.0pt;font-family:"Tahoma","sans-serif";color:black"> pw_forum-bounces@pwscf.org [mailto:pw_forum-bounces@pwscf.org]
<b>On Behalf Of </b>Arles V. Gil Rebaza<br>
<b>Sent:</b> Monday, June 16, 2014 11:28 AM<br>
<b>To:</b> PWSCF Forum<br>
<b>Subject:</b> Re: [Pw_forum] Nd2Fe14B Magnetization<o:p></o:p></span></p>
<p class="MsoNormal"><span style="color:black"> <o:p></o:p></span></p>
<div>
<p class="MsoNormal"><span style="color:black">Dear Michael, are you sure that "</span><span style="font-size:10.0pt;font-family:"Courier New";color:black">ecutwfc = 260.0" is a enough value for your system, may be you should increse, and you're using the default
value of "ecutrho", try to converge with this tag.!!!</span><span style="color:black"><o:p></o:p></span></p>
<div>
<p class="MsoNormal"><span style="color:black"> <o:p></o:p></span></p>
</div>
<div>
<p class="MsoNormal"><span style="font-size:10.0pt;font-family:"Courier New";color:black">Best</span><span style="color:black"><o:p></o:p></span></p>
</div>
<div>
<p class="MsoNormal"><span style="color:black"> <o:p></o:p></span></p>
</div>
<div>
<p class="MsoNormal"><span style="font-size:10.0pt;font-family:"Courier New";color:black">PhD. Arles V. GIl Rebaza</span><span style="color:black"><o:p></o:p></span></p>
</div>
<div>
<p class="MsoNormal"><span style="font-size:10.0pt;font-family:"Courier New";color:black">Instituto de Física La Plata</span><span style="color:black"><o:p></o:p></span></p>
</div>
<div>
<p class="MsoNormal"><span style="font-size:10.0pt;font-family:"Courier New";color:black">La Plata - Argentina</span><span style="color:black"><o:p></o:p></span></p>
</div>
</div>
<div>
<p class="MsoNormal" style="margin-bottom:12.0pt"><span style="color:black"> <o:p></o:p></span></p>
<div>
<p class="MsoNormal"><span style="color:black">2014-06-16 12:16 GMT-03:00 Kurth, Michael R. <<a href="mailto:kurth.17@buckeyemail.osu.edu" target="_blank">kurth.17@buckeyemail.osu.edu</a>>:<o:p></o:p></span></p>
<div>
<div>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto;background:white">
<span style="color:#282828"> </span><span style="color:black"><o:p></o:p></span></p>
<p style="background:white"><span style="font-family:"Calibri","sans-serif";color:black">Hello,</span><span style="color:black"><o:p></o:p></span></p>
<p style="background:white"><span style="font-family:"Calibri","sans-serif";color:black"> </span><span style="color:black"><o:p></o:p></span></p>
<p style="background:white"><span style="font-family:"Calibri","sans-serif";color:black">I am trying to calculate the magnetization of NdFeB which is a ferromagnetic material and a metal. This calculation has been done before using different methods than plane
waves and pseudopotentials. Moreover, the magnetization per formula unit has been experimentally measured. The previous calculations and experimental values agree and report a magnetization per formula unit of around 37 Bohr Magneton. The previous calculations
used the Generalized Gradient Approximation for the exchange correlation functional.</span><span style="color:black"><o:p></o:p></span></p>
<p style="background:white"><span style="font-family:"Calibri","sans-serif";color:black"> </span><span style="color:black"><o:p></o:p></span></p>
<p style="background:white"><span style="font-family:"Calibri","sans-serif";color:black">We first tried the calculation using the default number of bands and the input file is attached below:</span><span style="color:black"><o:p></o:p></span></p>
<p style="background:white"><span style="font-family:"Calibri","sans-serif";color:black"> </span><span style="color:black"><o:p></o:p></span></p>
<pre style="background:white;WORD-WRAP: break-word;white-space:pre-wrap"><span style="color:black">&control<o:p></o:p></span></pre>
<pre style="background:white"><span style="color:black"> calculation='scf'<o:p></o:p></span></pre>
<pre style="background:white"><span style="color:black"> restart_mode='from_scratch',<o:p></o:p></span></pre>
<pre style="background:white"><span style="color:black"> wf_collect=.TRUE.,<o:p></o:p></span></pre>
<pre style="background:white"><span style="color:black"> prefix='NdFeB',<o:p></o:p></span></pre>
<pre style="background:white"><span style="color:black"> pseudo_dir = '/usr/local/espresso/espresso-5.0.3/pseudo/',<o:p></o:p></span></pre>
<pre style="background:white"><span style="color:black"> outdir = '/fs/lustre/osu7834/NdFeB/'<o:p></o:p></span></pre>
<pre style="background:white"><span style="color:black">/<o:p></o:p></span></pre>
<pre style="background:white"><span style="color:black">&system<o:p></o:p></span></pre>
<pre style="background:white"><span style="color:black"> ibrav=0, nat=68, ntyp=3,<o:p></o:p></span></pre>
<pre style="background:white"><span style="color:black"> nspin=2, starting_magnetization(1)=0.6, starting_magnetization(2)=0.6, starting_magnetization(3)=0.6,<o:p></o:p></span></pre>
<pre style="background:white"><span style="color:black"> ecutwfc = 260.0<o:p></o:p></span></pre>
<pre style="background:white"><span style="color:black"> occupations='smearing', smearing='mv', degauss=0.002<o:p></o:p></span></pre>
<pre style="background:white"><span style="color:black">/<o:p></o:p></span></pre>
<pre style="background:white"><span style="color:black">&electrons<o:p></o:p></span></pre>
<pre style="background:white"><span style="color:black"> conv_thr=1.0e-6<o:p></o:p></span></pre>
<pre style="background:white"><span style="color:black"> mixing_beta=0.02<o:p></o:p></span></pre>
<pre style="background:white"><span style="color:black"> mixing_mode = 'local-TF'<o:p></o:p></span></pre>
<pre style="background:white"><span style="color:black">/<o:p></o:p></span></pre>
<pre style="background:white"><span style="color:black">CELL_PARAMETERS {angstrom}<o:p></o:p></span></pre>
<pre style="background:white"><span style="color:black">8.8 0.0 0.0<o:p></o:p></span></pre>
<pre style="background:white"><span style="color:black">0.0 12.19 0.0<o:p></o:p></span></pre>
<pre style="background:white"><span style="color:black">0.0 0.0 8.8<o:p></o:p></span></pre>
<pre style="background:white"><span style="color:black">ATOMIC_SPECIES<o:p></o:p></span></pre>
<pre style="background:white"><span style="color:black">Nd 144.24 Nd.pbe-mt_fhi.UPF<o:p></o:p></span></pre>
<pre style="background:white"><span style="color:black">Fe 55.845 Fe.pbe-mt_fhi.UPF <o:p></o:p></span></pre>
<pre style="background:white"><span style="color:black">B 10.811 B.pbe-mt_fhi.UPF <o:p></o:p></span></pre>
<pre style="background:white"><span style="color:black">ATOMIC_POSITIONS {angstrom}<o:p></o:p></span></pre>
<pre style="background:white"><span style="color:black">B 5.54400 0.00000 3.25600<o:p></o:p></span></pre>
<pre style="background:white"><span style="color:black">B 3.25600 0.00000 5.54400<o:p></o:p></span></pre>
<pre style="background:white"><span style="color:black">B 7.65600 6.09500 7.65600<o:p></o:p></span></pre>
<pre style="background:white"><span style="color:black">B 1.14400 6.09500 1.14400<o:p></o:p></span></pre>
<pre style="background:white"><span style="color:black">Fe 4.40000 0.00000 0.00000<o:p></o:p></span></pre>
<pre style="background:white"><span style="color:black">Fe 4.40000 6.09500 0.00000<o:p></o:p></span></pre>
<pre style="background:white"><span style="color:black">Fe 0.00000 6.09500 4.40000<o:p></o:p></span></pre>
<pre style="background:white"><span style="color:black">Fe 0.00000 0.00000 4.40000<o:p></o:p></span></pre>
<pre style="background:white"><span style="color:black">Fe 4.40000 1.37747 4.40000<o:p></o:p></span></pre>
<pre style="background:white"><span style="color:black">Fe 0.00000 7.47247 0.00000<o:p></o:p></span></pre>
<pre style="background:white"><span style="color:black">Fe 0.00000 4.71753 0.00000<o:p></o:p></span></pre>
<pre style="background:white"><span style="color:black">Fe 4.40000 10.81253 4.40000<o:p></o:p></span></pre>
<pre style="background:white"><span style="color:black">Fe 0.85360 2.49895 0.85360<o:p></o:p></span></pre>
<pre style="background:white"><span style="color:black">Fe 7.94640 2.49895 7.94640<o:p></o:p></span></pre>
<pre style="background:white"><span style="color:black">Fe 5.25360 8.59395 3.54640<o:p></o:p></span></pre>
<pre style="background:white"><span style="color:black">Fe 3.54640 8.59395 5.25360<o:p></o:p></span></pre>
<pre style="background:white"><span style="color:black">Fe 5.25360 3.59605 3.54640<o:p></o:p></span></pre>
<pre style="background:white"><span style="color:black">Fe 3.54640 3.59605 5.25360<o:p></o:p></span></pre>
<pre style="background:white"><span style="color:black">Fe 0.85360 9.69105 0.85360<o:p></o:p></span></pre>
<pre style="background:white"><span style="color:black">Fe 7.94640 9.69105 7.94640<o:p></o:p></span></pre>
<pre style="background:white"><span style="color:black">Fe 2.79840 3.01093 2.79840<o:p></o:p></span></pre>
<pre style="background:white"><span style="color:black">Fe 6.00160 3.01093 6.00160<o:p></o:p></span></pre>
<pre style="background:white"><span style="color:black">Fe 7.19840 9.10593 1.60160<o:p></o:p></span></pre>
<pre style="background:white"><span style="color:black">Fe 1.60160 9.10593 7.19840<o:p></o:p></span></pre>
<pre style="background:white"><span style="color:black">Fe 7.19840 3.08407 1.60160<o:p></o:p></span></pre>
<pre style="background:white"><span style="color:black">Fe 1.60160 3.08407 7.19840<o:p></o:p></span></pre>
<pre style="background:white"><span style="color:black">Fe 2.79840 9.17907 2.79840<o:p></o:p></span></pre>
<pre style="background:white"><span style="color:black">Fe 6.00160 9.17907 6.00160<o:p></o:p></span></pre>
<pre style="background:white"><span style="color:black">Fe 4.99840 1.56032 1.97120<o:p></o:p></span></pre>
<pre style="background:white"><span style="color:black">Fe 3.80160 1.56032 6.82880<o:p></o:p></span></pre>
<pre style="background:white"><span style="color:black">Fe 6.37120 7.65532 8.20160<o:p></o:p></span></pre>
<pre style="background:white"><span style="color:black">Fe 2.42880 7.65532 0.59840<o:p></o:p></span></pre>
<pre style="background:white"><span style="color:black">Fe 0.59840 4.53468 2.42880<o:p></o:p></span></pre>
<pre style="background:white"><span style="color:black">Fe 8.20160 4.53468 6.37120<o:p></o:p></span></pre>
<pre style="background:white"><span style="color:black">Fe 1.97120 10.62968 4.99840<o:p></o:p></span></pre>
<pre style="background:white"><span style="color:black">Fe 6.82880 1.56032 3.80160<o:p></o:p></span></pre>
<pre style="background:white"><span style="color:black">Fe 3.80160 10.62968 6.82880<o:p></o:p></span></pre>
<pre style="background:white"><span style="color:black">Fe 4.99840 10.62968 1.97120<o:p></o:p></span></pre>
<pre style="background:white"><span style="color:black">Fe 2.42880 4.53468 0.59840<o:p></o:p></span></pre>
<pre style="background:white"><span style="color:black">Fe 6.37120 4.53468 8.20160<o:p></o:p></span></pre>
<pre style="background:white"><span style="color:black">Fe 8.20160 7.65532 6.37120<o:p></o:p></span></pre>
<pre style="background:white"><span style="color:black">Fe 0.59840 7.65532 2.42880<o:p></o:p></span></pre>
<pre style="background:white"><span style="color:black">Fe 6.82880 10.62968 3.80160<o:p></o:p></span></pre>
<pre style="background:white"><span style="color:black">Fe 1.97120 1.56032 4.99840<o:p></o:p></span></pre>
<pre style="background:white"><span style="color:black">Fe 3.15920 2.14544 0.34320<o:p></o:p></span></pre>
<pre style="background:white"><span style="color:black">Fe 5.64080 2.14544 8.45680<o:p></o:p></span></pre>
<pre style="background:white"><span style="color:black">Fe 4.74320 8.24044 1.24080<o:p></o:p></span></pre>
<pre style="background:white"><span style="color:black">Fe 4.05680 8.24044 7.55920<o:p></o:p></span></pre>
<pre style="background:white"><span style="color:black">Fe 7.55920 3.94956 4.05680<o:p></o:p></span></pre>
<pre style="background:white"><span style="color:black">Fe 1.24080 3.94956 4.74320<o:p></o:p></span></pre>
<pre style="background:white"><span style="color:black">Fe 0.34320 10.04456 3.15920<o:p></o:p></span></pre>
<pre style="background:white"><span style="color:black">Fe 8.45680 2.14544 5.64080<o:p></o:p></span></pre>
<pre style="background:white"><span style="color:black">Fe 5.64080 10.04456 8.45680<o:p></o:p></span></pre>
<pre style="background:white"><span style="color:black">Fe 3.15920 10.04456 0.34320<o:p></o:p></span></pre>
<pre style="background:white"><span style="color:black">Fe 4.05680 3.94956 7.55920<o:p></o:p></span></pre>
<pre style="background:white"><span style="color:black">Fe 4.74320 3.94956 1.24080<o:p></o:p></span></pre>
<pre style="background:white"><span style="color:black">Fe 1.24080 8.24044 4.74320<o:p></o:p></span></pre>
<pre style="background:white"><span style="color:black">Fe 7.55920 8.24044 4.05680<o:p></o:p></span></pre>
<pre style="background:white"><span style="color:black">Fe 8.45680 10.04456 5.64080<o:p></o:p></span></pre>
<pre style="background:white"><span style="color:black">Fe 0.34320 2.14544 3.15920<o:p></o:p></span></pre>
<pre style="background:white"><span style="color:black">Nd 2.34080 0.00000 2.34080<o:p></o:p></span></pre>
<pre style="background:white"><span style="color:black">Nd 6.45920 0.00000 6.45920<o:p></o:p></span></pre>
<pre style="background:white"><span style="color:black">Nd 6.74080 6.09500 2.05920<o:p></o:p></span></pre>
<pre style="background:white"><span style="color:black">Nd 2.05920 6.09500 6.74080<o:p></o:p></span></pre>
<pre style="background:white"><span style="color:black">Nd 7.56800 0.00000 1.23200<o:p></o:p></span></pre>
<pre style="background:white"><span style="color:black">Nd 1.23200 0.00000 7.56800<o:p></o:p></span></pre>
<pre style="background:white"><span style="color:black">Nd 5.63200 6.09500 5.63200<o:p></o:p></span></pre>
<pre style="background:white"><span style="color:black">Nd 3.16800 6.09500 3.16800<o:p></o:p></span></pre>
<pre style="background:white"><span style="color:black"> <o:p></o:p></span></pre>
<pre style="background:white"><span style="color:black">K_POINTS {automatic}<o:p></o:p></span></pre>
<pre style="background:white"><span style="color:black">4 4 4 0 0 0<o:p></o:p></span></pre>
<pre style="background:white"><span style="color:black"> <o:p></o:p></span></pre>
<pre style="background:white"><span style="color:black"> <o:p></o:p></span></pre>
<p style="background:white"><span style="font-family:"Calibri","sans-serif";color:black">This input file was prepared after checking convergence with respect to the energy cutoff and the number of points and did converge after 2 restart runs. Unfortunately,
we obtained a magnetization of 25 Bohr Magneton/f.u. (which is too small).</span><span style="color:black"><o:p></o:p></span></p>
<p style="background:white"><span style="font-family:"Calibri","sans-serif";color:black"> </span><span style="color:black"><o:p></o:p></span></p>
<p style="background:white"><span style="font-family:"Calibri","sans-serif";color:black">Furthermore, and after checking for convergence with respect to the number of bands (by adding more bands) we found that we end up with a magnetization of 50 Bohr Magneton/f.u.
(which is too big).</span><span style="color:black"><o:p></o:p></span></p>
<p style="background:white"><span style="font-family:"Calibri","sans-serif";color:black">I also think that the degauss value used for smearing is small enough and can't be the problem.</span><span style="color:black"><o:p></o:p></span></p>
<p style="background:white"><span style="font-family:"Calibri","sans-serif";color:black">Are there any recommendations for what other parameters I could examine to help me get the right answer?</span><span style="color:black"><o:p></o:p></span></p>
<p style="background:white"><span style="font-family:"Calibri","sans-serif";color:black"> </span><span style="color:black"><o:p></o:p></span></p>
<p style="background:white"><span style="font-family:"Calibri","sans-serif";color:black">Your help is greatly appreciated!</span><span style="color:black"><o:p></o:p></span></p>
<p style="background:white"><span style="font-family:"Calibri","sans-serif";color:black"> </span><span style="color:black"><o:p></o:p></span></p>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"><span style="color:black"> <o:p></o:p></span></p>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"><span style="color:black">Sincerely,<o:p></o:p></span></p>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"><span style="color:black">Michael Kurth<o:p></o:p></span></p>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"><span style="color:black"> <o:p></o:p></span></p>
<p style="background:white"><span style="font-family:"Calibri","sans-serif";color:black">Graduate Student in Nuclear Engineering</span><span style="color:black"><o:p></o:p></span></p>
<p style="background:white"><span style="font-family:"Calibri","sans-serif";color:black">Department of Mechanical and Aerospace Engineering</span><span style="color:black"><o:p></o:p></span></p>
<p style="background:white"><span style="font-family:"Calibri","sans-serif";color:black">The Ohio State University
</span><span style="color:black"><o:p></o:p></span></p>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"><span style="color:black"> <o:p></o:p></span></p>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"><span style="color:black"> <o:p></o:p></span></p>
</div>
</div>
<p class="MsoNormal"><span style="color:black"><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><o:p></o:p></span></p>
</div>
<p class="MsoNormal"><span style="color:black"><br>
<br clear="all">
<o:p></o:p></span></p>
<div>
<p class="MsoNormal"><span style="color:black"> <o:p></o:p></span></p>
</div>
<p class="MsoNormal"><span style="color:black">-- <br>
###---------> Arles V. <---------### <o:p></o:p></span></p>
</div>
<div class="MsoNormal" align="center" style="text-align:center"><span style="color:black">
<hr size="1" width="100%" noshade="" style="color:#A0A0A0" align="center">
</span></div>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"><span style="color:black">No virus found in this message.<br>
Checked by AVG - <a href="http://www.avg.com">www.avg.com</a><br>
Version: 2014.0.4592 / Virus Database: 3964/7686 - Release Date: 06/16/14<o:p></o:p></span></p>
<div class="MsoNormal" align="center" style="text-align:center"><span style="font-family:"Calibri","sans-serif";color:black">
<hr size="3" width="100%" align="center">
</span></div>
<p class="MsoNormal"><span style="font-family:"Calibri","sans-serif";color:black">_______________________________________________<br>
Pw_forum mailing list<br>
Pw_forum@pwscf.org<br>
http://pwscf.org/mailman/listinfo/pw_forum<o:p></o:p></span></p>
</div>
</div>
</div>
<p class="MsoNormal" style="margin-bottom:12.0pt"><o:p> </o:p></p>
<div class="MsoNormal" align="center" style="text-align:center">
<hr size="1" width="99%" noshade="" style="color:#909090" align="center">
</div>
<table class="MsoNormalTable" border="0" cellspacing="0" cellpadding="0" style="border-collapse:collapse">
<tbody>
<tr>
<td style="padding:0in 11.25pt 0in 6.0pt">
<p class="MsoNormal"><a href="http://www.avast.com/"><span style="border:solid windowtext 1.0pt;padding:0in;text-decoration:none"><img border="0" width="100" height="100" id="_x0000_i1028" src="cid:~WRD000.jpg" alt="Image removed by sender."></span></a><o:p></o:p></p>
</td>
<td style="padding:.75pt .75pt .75pt .75pt">
<p><span style="font-family:"Calibri","sans-serif";color:#3D4D5A">Este mensaje no contiene virus ni malware porque la protección de
<a href="http://www.avast.com/">avast! Antivirus</a> está activa. <o:p></o:p></span></p>
</td>
</tr>
</tbody>
</table>
<p class="MsoNormal"><o:p> </o:p></p>
<div class="MsoNormal" align="center" style="text-align:center">
<hr size="1" width="100%" noshade="" style="color:#A0A0A0" align="center">
</div>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto">No virus found in this message.<br>
Checked by AVG - <a href="http://www.avg.com">www.avg.com</a><br>
Version: 2014.0.4592 / Virus Database: 3964/7686 - Release Date: 06/16/14<o:p></o:p></p>
</div>
</body>
</html>