<div dir="ltr"><div class="gmail_default" style="font-family:verdana,sans-serif;color:#000000">If it is in the first cycle, this is a nasty problem that arises because the initial density is not always very good. I think you are going to have to adjust (by hand) the RMTs until you find something which avoids this problem. Later you can use expand_rmt_lapw (which works for changes as well) to adjust. My hunch is to use a larger RMT for Sn, e.g. </div><div class="gmail_default" style="font-family:verdana,sans-serif;color:#000000"><br></div><div class="gmail_default" style="font-family:verdana,sans-serif;color:#000000">setrmt case -a Na:2.0,Fe:1.9,Sn:2.0,O:1.7<br></div><div class="gmail_default" style="font-family:verdana,sans-serif;color:#000000">cp case.struct_setrmt case.struct</div><div class="gmail_default" style="font-family:verdana,sans-serif;color:#000000"><br></div><div class="gmail_default" style="font-family:verdana,sans-serif;color:#000000">(replace "case" by whatever"). Then run "x lstart ; x dstart ; x dstart -up ; x dstart -dn" and try again. So long as you can get past the first interaction you will be up and running.</div><div class="gmail_default" style="font-family:verdana,sans-serif;color:#000000"><br></div><div class="gmail_default" style="font-family:verdana,sans-serif;color:#000000">N.B., Peter may have some other ideas.</div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Tue, Feb 4, 2020 at 12:30 PM Biswajit Mandal <<a href="mailto:biswajit.mandal095@gmail.com" target="_blank">biswajit.mandal095@gmail.com</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr"><div dir="ltr"><div>Dear sir,<br></div><div>Thank you very much for your suggestions.</div><div>I have again simulated with -7.0Ry separation energy and encountered the same error in first cycle. For this calculation, the parent material is NaFeSnO4 which is an anti-ferromagnetic material confirmed experimentally. The anti-ferromagnetic nature arises due to antiparallel spin orientation of Fe atoms, distributed in two different 4c sites. Each of these 4c sites are equally shared by Fe and Sn atom. In this case, we have broken symmetry using vasp and run spin-polarized calculation with GGA+U scheme. For Fe atoms we have taken Ueff=6eV. The case.inst for this simulation is as follows;</div><div>Na<br>Ne 1  <br>3,-1,0.5  N<br>3,-1,0.5  N<br>Na<br>Ne 1  <br>3,-1,0.5  N<br>3,-1,0.5  N<br>Fe<br>Ar 3  <br>3, 2,2.0  N<br>3, 2,2.0  N<br>3,-3,0.0  N<br>3,-3,2.5  N<br>4,-1,0.5  N<br>4,-1,1.0  N<br>Fe<br>Ar 3  <br>3, 2,2.0  N<br>3, 2,2.0  N<br>3,-3,2.5  N<br>3,-3,0.0  N<br>4,-1,1.0  N<br>4,-1,0.5  N<br>Sn<br>Kr 5  <br>4, 2,2.0  N<br>4, 2,2.0  N<br>4,-3,3.0  N<br>4,-3,3.0  N<br>5,-1,1.0  N<br>5,-1,1.0  N<br>5, 1,1.0  N<br>5, 1,1.0  N<br>5,-2,0.0  N<br>5,-2,0.0  N<br>Sn<br>Kr 5  <br>4, 2,2.0  N<br>4, 2,2.0  N<br>4,-3,3.0  N<br>4,-3,3.0  N<br>5,-1,1.0  N<br>5,-1,1.0  N<br>5, 1,1.0  N<br>5, 1,1.0  N<br>5,-2,0.0  N<br>5,-2,0.0  N<br>Fe<br>Ar 3  <br>3, 2,2.0  N<br>3, 2,2.0  N<br>3,-3,2.5  N<br>3,-3,0.0  N<br>4,-1,1.0  N<br>4,-1,0.5  N<br>Fe<br>Ar 3  <br>3, 2,2.0  N<br>3, 2,2.0  N<br>3,-3,0.0  N<br>3,-3,2.5  N<br>4,-1,0.5  N<br>4,-1,1.0  N<br>Sn<br>Kr 5  <br>4, 2,2.0  N<br>4, 2,2.0  N<br>4,-3,3.0  N<br>4,-3,3.0  N<br>5,-1,1.0  N<br>5,-1,1.0  N<br>5, 1,1.0  N<br>5, 1,1.0  N<br>5,-2,0.0  N<br>5,-2,0.0  N<br>Sn<br>Kr 5  <br>4, 2,2.0  N<br>4, 2,2.0  N<br>4,-3,3.0  N<br>4,-3,3.0  N<br>5,-1,1.0  N<br>5,-1,1.0  N<br>5, 1,1.0  N<br>5, 1,1.0  N<br>5,-2,0.0  N<br>5,-2,0.0  N<br>O<br>He 3  <br>2,-1,1.0  N<br>2,-1,1.0  N<br>2, 1,1.0  N<br>2, 1,1.0  N<br>2,-2,1.0  N<br>2,-2,1.0  N<br>O<br>He 3  <br>2,-1,1.0  N<br>2,-1,1.0  N<br>2, 1,1.0  N<br>2, 1,1.0  N<br>2,-2,1.0  N<br>2,-2,1.0  N<br>O<br>He 3  <br>2,-1,1.0  N<br>2,-1,1.0  N<br>2, 1,1.0  N<br>2, 1,1.0  N<br>2,-2,1.0  N<br>2,-2,1.0  N<br>O<br>He 3  <br>2,-1,1.0  N<br>2,-1,1.0  N<br>2, 1,1.0  N<br>2, 1,1.0  N<br>2,-2,1.0  N<br>2,-2,1.0  N<br>O<br>He 3  <br>2,-1,1.0  N<br>2,-1,1.0  N<br>2, 1,1.0  N<br>2, 1,1.0  N<br>2,-2,1.0  N<br>2,-2,1.0  N<br>O<br>He 3  <br>2,-1,1.0  N<br>2,-1,1.0  N<br>2, 1,1.0  N<br>2, 1,1.0  N<br>2,-2,1.0  N<br>2,-2,1.0  N<br>O<br>He 3  <br>2,-1,1.0  N<br>2,-1,1.0  N<br>2, 1,1.0  N<br>2, 1,1.0  N<br>2,-2,1.0  N<br>2,-2,1.0  N<br>O<br>He 3  <br>2,-1,1.0  N<br>2,-1,1.0  N<br>2, 1,1.0  N<br>2, 1,1.0  N<br>2,-2,1.0  N<br>2,-2,1.0  N<br>O<br>He 3  <br>2,-1,1.0  N<br>2,-1,1.0  N<br>2, 1,1.0  N<br>2, 1,1.0  N<br>2,-2,1.0  N<br>2,-2,1.0  N<br>O<br>He 3  <br>2,-1,1.0  N<br>2,-1,1.0  N<br>2, 1,1.0  N<br>2, 1,1.0  N<br>2,-2,1.0  N<br>2,-2,1.0  N<br>O<br>He 3  <br>2,-1,1.0  N<br>2,-1,1.0  N<br>2, 1,1.0  N<br>2, 1,1.0  N<br>2,-2,1.0  N<br>2,-2,1.0  N<br>O<br>He 3  <br>2,-1,1.0  N<br>2,-1,1.0  N<br>2, 1,1.0  N<br>2, 1,1.0  N<br>2,-2,1.0  N<br>2,-2,1.0  N<br>O<br>He 3  <br>2,-1,1.0  N<br>2,-1,1.0  N<br>2, 1,1.0  N<br>2, 1,1.0  N<br>2,-2,1.0  N<br>2,-2,1.0  N<br>O<br>He 3  <br>2,-1,1.0  N<br>2,-1,1.0  N<br>2, 1,1.0  N<br>2, 1,1.0  N<br>2,-2,1.0  N<br>2,-2,1.0  N<br>O<br>He 3  <br>2,-1,1.0  N<br>2,-1,1.0  N<br>2, 1,1.0  N<br>2, 1,1.0  N<br>2,-2,1.0  N<br>2,-2,1.0  N<br>O<br>He 3  <br>2,-1,1.0  N<br>2,-1,1.0  N<br>2, 1,1.0  N<br>2, 1,1.0  N<br>2,-2,1.0  N<br>2,-2,1.0  N<br>****<br>****         END of input (instgen_lapw)<br></div><div>Is this approach wrong? If wrong, then kindly help me to break the symmetry so that I could get some meaningful result. <br></div><div>Thanks in advance</div><div>With regards<br></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Tue, Feb 4, 2020 at 9:33 PM Laurence Marks <<a href="mailto:laurence.marks@gmail.com" target="_blank">laurence.marks@gmail.com</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr"><div dir="ltr"><div style="font-family:verdana,sans-serif;color:rgb(0,0,0)">Several other points:</div><div style="font-family:verdana,sans-serif;color:rgb(0,0,0)">1) It appears that you did not to follow the standard initiation. Your atom positions have errors, for instance atom 2 is</div><div style="font-family:verdana,sans-serif;color:rgb(0,0,0)"><span style="color:rgb(34,34,34)">ATOM  -2: X=0.77850183 Y=0.75000002 Z=0.15072562 </span> and should be<br></div><div style="font-family:verdana,sans-serif;color:rgb(0,0,0)"><span style="color:rgb(34,34,34)">ATOM  -2: X=0.77850183 Y=0.75000000 Z=0.15072562 </span><br></div><div style="font-family:verdana,sans-serif;color:rgb(0,0,0)"><span style="color:rgb(34,34,34)">This can matter. With this corrected Wien will (properly) switch axes to put the mirror down z.</span></div><div style="font-family:verdana,sans-serif;color:rgb(0,0,0)"><span style="color:rgb(34,34,34)"><br></span></div><div style="font-family:verdana,sans-serif;color:rgb(0,0,0)"><span style="color:rgb(34,34,34)">2) It looks like you used a non-standard energy in lstart, to avoid core leaking. These are warnings, and you should not adjust too much particularly when you have a complex structure with (perhaps) inaccurate atomic positions. I suggest just using -7, your use of -10 (?) was not needed.</span></div><div style="font-family:verdana,sans-serif;color:rgb(0,0,0)"><span style="color:rgb(34,34,34)"><br></span></div><div style="font-family:verdana,sans-serif;color:rgb(0,0,0)"><span style="color:rgb(34,34,34)">3) Beware that Vasp can reduce symmetry inappropriately. I would start of with the published Pnma symmetry and occupy one of the Fe sites and one of the Sn sites.</span></div><div style="font-family:verdana,sans-serif;color:rgb(0,0,0)"><span style="color:rgb(34,34,34)"><br></span></div><div style="font-family:verdana,sans-serif;color:rgb(0,0,0)"><span style="color:rgb(34,34,34)">4) You will need to look into the magnetism. I don't know if this is ferromagnetic or anti-ferromagnetic. If the later you will need to break the symmetry appropriately (carefully).</span></div><div style="font-family:verdana,sans-serif;color:rgb(0,0,0)"><span style="color:rgb(34,34,34)"><br></span></div><div style="font-family:verdana,sans-serif;color:rgb(0,0,0)"><span style="color:rgb(34,34,34)">5) I doubt that you will get meaningful results unless you use either -eece or +U for the Fe sites.</span></div><div style="font-family:verdana,sans-serif;color:rgb(0,0,0)"><br></div></div><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Tue, Feb 4, 2020 at 9:25 AM Peter Blaha <<a href="mailto:pblaha@theochem.tuwien.ac.at" target="_blank">pblaha@theochem.tuwien.ac.at</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">Does the problem happen in the first scf cycle or later on ???<br>
<br>
If it is later, make sure :NEC01 is ok in the scf fle.<br>
<br>
Check case.scf1 to see where the actual E-parameters for atom 5 l=2 are.<br>
If the E-parameters for APW and for the LO are too close ---> ghostbands<br>
<br>
<br>
Are you sure you need the Sn-s semicore states ??<br>
<br>
<br><br></blockquote></div>-- <br><div dir="ltr"><div dir="ltr">Professor Laurence Marks<br>Department of Materials Science and Engineering<br>Northwestern University<br><a href="http://www.numis.northwestern.edu/" target="_blank">www.numis.northwestern.edu</a><div>Corrosion in 4D: <a href="http://www.numis.northwestern.edu/MURI" target="_blank">www.numis.northwestern.edu/MURI</a><br>Co-Editor, Acta Cryst A<br>"Research is to see what everybody else has seen, and to think what nobody else has thought"<br>Albert Szent-Gyorgi</div></div></div></div>
_______________________________________________<br>
Wien mailing list<br>
<a href="mailto:Wien@zeus.theochem.tuwien.ac.at" target="_blank">Wien@zeus.theochem.tuwien.ac.at</a><br>
<a href="https://urldefense.proofpoint.com/v2/url?u=http-3A__zeus.theochem.tuwien.ac.at_mailman_listinfo_wien&d=DwMFaQ&c=yHlS04HhBraes5BQ9ueu5zKhE7rtNXt_d012z2PA6ws&r=U_T4PL6jwANfAy4rnxTj8IUxm818jnvqKFdqWLwmqg0&m=yJHoJ2zB-ZmPMjx4PujXltOpTcpFcWYIO1WaiwqIZlo&s=bJepWz-EvJEJ0BJXaZJcxQ1dj9Rsy4bs2rm9-JixDFA&e=" rel="noreferrer" target="_blank">http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien</a><br>
SEARCH the MAILING-LIST at:  <a href="https://urldefense.proofpoint.com/v2/url?u=http-3A__www.mail-2Darchive.com_wien-40zeus.theochem.tuwien.ac.at_index.html&d=DwMFaQ&c=yHlS04HhBraes5BQ9ueu5zKhE7rtNXt_d012z2PA6ws&r=U_T4PL6jwANfAy4rnxTj8IUxm818jnvqKFdqWLwmqg0&m=yJHoJ2zB-ZmPMjx4PujXltOpTcpFcWYIO1WaiwqIZlo&s=qnJ_6hNMuS6fCr8hxP3uKtEHysBsUs8Ekpo9LAWW1N8&e=" rel="noreferrer" target="_blank">http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html</a><br>
</blockquote></div><br clear="all"><br>-- <br><div dir="ltr"><div dir="ltr"><div><div><div><div><div><div>Yours Sincerely<br>------------------------<br></div>Biswajit Mandal<br></div>PhD Research Scholar<br></div>Department of Physics<br></div>IIT Patna<br></div>Bihta, Patna-801103<br></div>INDIA<br></div></div>
</div>
_______________________________________________<br>
Wien mailing list<br>
<a href="mailto:Wien@zeus.theochem.tuwien.ac.at" target="_blank">Wien@zeus.theochem.tuwien.ac.at</a><br>
<a href="https://urldefense.proofpoint.com/v2/url?u=http-3A__zeus.theochem.tuwien.ac.at_mailman_listinfo_wien&d=DwICAg&c=yHlS04HhBraes5BQ9ueu5zKhE7rtNXt_d012z2PA6ws&r=U_T4PL6jwANfAy4rnxTj8IUxm818jnvqKFdqWLwmqg0&m=yJHoJ2zB-ZmPMjx4PujXltOpTcpFcWYIO1WaiwqIZlo&s=bJepWz-EvJEJ0BJXaZJcxQ1dj9Rsy4bs2rm9-JixDFA&e=" rel="noreferrer" target="_blank">https://urldefense.proofpoint.com/v2/url?u=http-3A__zeus.theochem.tuwien.ac.at_mailman_listinfo_wien&d=DwICAg&c=yHlS04HhBraes5BQ9ueu5zKhE7rtNXt_d012z2PA6ws&r=U_T4PL6jwANfAy4rnxTj8IUxm818jnvqKFdqWLwmqg0&m=yJHoJ2zB-ZmPMjx4PujXltOpTcpFcWYIO1WaiwqIZlo&s=bJepWz-EvJEJ0BJXaZJcxQ1dj9Rsy4bs2rm9-JixDFA&e=</a> <br>
SEARCH the MAILING-LIST at:  <a href="https://urldefense.proofpoint.com/v2/url?u=http-3A__www.mail-2Darchive.com_wien-40zeus.theochem.tuwien.ac.at_index.html&d=DwICAg&c=yHlS04HhBraes5BQ9ueu5zKhE7rtNXt_d012z2PA6ws&r=U_T4PL6jwANfAy4rnxTj8IUxm818jnvqKFdqWLwmqg0&m=yJHoJ2zB-ZmPMjx4PujXltOpTcpFcWYIO1WaiwqIZlo&s=qnJ_6hNMuS6fCr8hxP3uKtEHysBsUs8Ekpo9LAWW1N8&e=" rel="noreferrer" target="_blank">https://urldefense.proofpoint.com/v2/url?u=http-3A__www.mail-2Darchive.com_wien-40zeus.theochem.tuwien.ac.at_index.html&d=DwICAg&c=yHlS04HhBraes5BQ9ueu5zKhE7rtNXt_d012z2PA6ws&r=U_T4PL6jwANfAy4rnxTj8IUxm818jnvqKFdqWLwmqg0&m=yJHoJ2zB-ZmPMjx4PujXltOpTcpFcWYIO1WaiwqIZlo&s=qnJ_6hNMuS6fCr8hxP3uKtEHysBsUs8Ekpo9LAWW1N8&e=</a> <br>
</blockquote></div><br clear="all"><div><br></div>-- <br><div dir="ltr"><div dir="ltr">Professor Laurence Marks<br>Department of Materials Science and Engineering<br>Northwestern University<br><a href="http://www.numis.northwestern.edu/" target="_blank">www.numis.northwestern.edu</a><div>Corrosion in 4D: <a href="http://www.numis.northwestern.edu/MURI" target="_blank">www.numis.northwestern.edu/MURI</a><br>Co-Editor, Acta Cryst A<br>"Research is to see what everybody else has seen, and to think what nobody else has thought"<br>Albert Szent-Gyorgi</div></div></div>