<div><br></div><div><div>Dear <span style="font-size: 12.8px; line-height: 19.2px;">Professor Laurence Marks,</span></div><div><span style="font-size: 12.8px; line-height: 19.2px;">Thanks for your kind reply.</span></div><div><span style="font-size: 12.8px; line-height: 19.2px;">Mingcui</span></div><div><span style="font-size: 12.8px; line-height: 19.2px;"><br></span></div><div><br></div><div style="font-size: 12px;font-family: Arial Narrow;padding:2px 0 2px 0;">------------------ 原始邮件 ------------------</div><div style="font-size: 12px;background:#efefef;padding:8px;"><div><b>发件人:</b> "Laurence Marks";<L-marks@northwestern.edu>;</div><div><b>发送时间:</b> 2016年4月22日(星期五) 上午6:57</div><div><b>收件人:</b> "A Mailing list for WIEN2k users"<wien@zeus.theochem.tuwien.ac.at>; <wbr></div><div></div><div><b>主题:</b> Re: [Wien] lsda</div></div><div><br></div><div dir="ltr">In general, to do this calculation you need to find another code.<div><br></div><div>If you are an expert in using Wien2k there is a method that might work, although I doubt that it will converge. You would need to do something like write a fortran (or C, or C++) routine that will average the density within the RMT for a given atom, for instance using the files in SRC_clmaddsub as a guide. Then edit the runsp_lapw script to insert your program.</div><div><br></div><div>Since this is not a very physically realistic thing to do I doubt that you will find many people prepared to help.</div></div><div class="gmail_extra"><br><div class="gmail_quote">On Thu, Apr 21, 2016 at 5:44 PM, ding <span dir="ltr"><<a href="mailto:dingmingcui@qq.com" target="_blank">dingmingcui@qq.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">



<div>
<div><span style="font-family:Tahoma;font-size:13.3333px;line-height:22.6667px">Dear WIEN2K users,</span></div>
<div><span style="font-family:Tahoma;font-size:13.3333px;line-height:22.6667px">I want to  </span><font face="Tahoma"><span style="font-size:13.3333px;line-height:22.6667px">force everything but  only one atom d states to be non-spinpolarized. For
 example, there are two Fe atoms in a unit cell, and I want to keep only one Fe atom to be spinpolarized, while the other to be non-spinpolarized. Would anyone tell me how to carry out this calculation? </span></font></div>
<div><font face="Tahoma"><span style="font-size:13.3333px;line-height:22.6667px"> </span></font></div>
<div><font face="Tahoma"><span style="font-size:13.3333px;line-height:22.6667px">Thanks and best regards</span></font></div>
<div><font face="Tahoma"><span style="font-size:13.3333px;line-height:22.6667px">Mingcui</span></font></div>
<div><font face="Tahoma"><span style="font-size:13.3333px;line-height:22.6667px"><br>
</span></font></div>
</div>

</blockquote></div><br><br clear="all"><div><br></div>-- <br><div class="gmail_signature"><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><span style="font-size:12.8px">Professor Laurence Marks</span><br></div><div dir="ltr"><span style="font-size:12.8px">"Research is to see what everybody else has seen, and to think what nobody else has thought", </span><span style="font-size:12.8px">Albert Szent-Gyorgi</span><br><a href="http://www.numis.northwestern.edu" target="_blank">www.numis.northwestern.edu</a> ; <span style="font-size:12.8px">Corrosion in 4D: </span><a href="http://MURI4D.numis.northwestern.edu" style="font-size:12.8px" target="_blank">MURI4D.numis.northwestern.edu</a><div><span style="font-size:12.8px">Partner of the CFW 100% program for gender equity, </span><a href="http://www.cfw.org/100-percent" style="font-size:12.8px" target="_blank">www.cfw.org/100-percent</a></div><div>Co-Editor, Acta Cryst A</div></div></div></div></div></div></div></div></div></div></div>
</div></div>