Thank you very much for your suggestions!<br><br><div class="gmail_quote">On Fri, Oct 15, 2010 at 8:01 PM, Peter Blaha <span dir="ltr"><<a href="mailto:pblaha@theochem.tuwien.ac.at">pblaha@theochem.tuwien.ac.at</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin: 0pt 0pt 0pt 0.8ex; border-left: 1px solid rgb(204, 204, 204); padding-left: 1ex;">No, you are NOT using P1 symmetry (but you have 24 symops.)<br>
<br>
As I said, you have to find a "common" symmetry for all your cases. Your struct file definitely<br>
does not do this.<br>
<br>
If you are unable to manage this by yourself, eventually you may have to use a trick like<br>
"adding a 3rd atom at some arbitrary position", so that the WIEN2k initialization programs<br>
really find P1 symmetry with only ONE symop. (or the common one). Stop the initialization<br>
after kgen.<br>
After you have this, remove the 3rd atom from struct and in* files and run lstart and dstart<br>
to complete the initialization.<br>
The runsp<br>
save_lapw<br>
The do the spinorbit+lapw2 for the two directions.<br>
<br>
PS: Please note: As far as I remember, the MAE of Co is only a few mycroRy, i.e. one had to use<br>
enormous k-meshes to get anything meaningful (and LDA/GGA is wrong anyway ?). Check literature.<br>
<br>
<br>
Am 13.10.2010 10:27, schrieb Bin Shao:<br>
<blockquote class="gmail_quote" style="margin: 0pt 0pt 0pt 0.8ex; border-left: 1px solid rgb(204, 204, 204); padding-left: 1ex;"><div class="im">
Dear Prof. Peter Blaha,<br>
<br>
With your suggestion, I recalculate the hcp Co in P1 symmetry with experiment parameters, and the k-mesh is 40x40x24. Then I get the energy in (001), (100) and (010) direction.<br>
But I find that the energy in (100) is large than (001) and that in (010) is smaller than (001). In my opinion, they should be the same in (100) and in (010) direction.<br>
<br>
Please give me some comments, and the attachment is my struct file.<br>
<br>
Thank you in advance!<br>
<br>
Best,<br>
<br></div><div class="im">
On Tue, Oct 12, 2010 at 2:46 PM, Peter Blaha <<a href="mailto:pblaha@theochem.tuwien.ac.at" target="_blank">pblaha@theochem.tuwien.ac.at</a> <mailto:<a href="mailto:pblaha@theochem.tuwien.ac.at" target="_blank">pblaha@theochem.tuwien.ac.at</a>>> wrote:<br>
<br>
If you change the k-mesh, you need to recalculate also lapw1.<br>
<br>
However, the proper approach is to find a symmetry (usually the "lower" one), which accommodates<br>
both directions of the magnetization (usually you can run also the "higher" symmetry case with the<br>
low symmetry). Run already the non-SO calculation in this symmetry and then use only<br>
lapwso/lapw2 with the two directions in case.inso.<br>
<br>
<br>
<br>
Am 12.10.2010 05:24, schrieb Bin Shao:<br>
<br>
Dear all,<br>
<br>
I intend to calculate the MAE of hcp Co with force theorem. After the nosoc scf calculation, I add the soc non-selfconsistently with the directions of M || c and M || a,<br>
respectively. I use the initso_lapw to prepare the input files which creats new structs and new klists for spin-polarized case. Then I run the program by commands<br>
<br>
x lapwso -up -p<br>
x lapw2 -c -up -p<br>
x lapw2 -c -dn -p<br>
<br>
But here comes some errors in lapw2,<br>
<br>
"dnlapw2.error"<br>
----------------------------------------------------------------------------------------<br>
'FERMI' - number of k-points inconsistent when reading kgen<br>
'FERMI' - check IN1 and KGEN files!<br>
** testerror: Error in Parallel LAPW2<br>
----------------------------------------------------------------------------------------<br>
<br>
Before Co hcp, I calculated the case of Fe monolayer with the same approach, but no errors. How to deal with the problem or do I need to provide some more input files?<br>
<br>
Any suggestion will be appriciated, thank you in advance!<br>
<br>
Best regards,<br>
<br>
--<br>
Bin Shao, Ph.D. Candidate<br>
College of Information Technical Science, Nankai University<br>
94 Weijin Rd. Nankai Dist. Tianjin 300071, China<br></div>
Email: <a href="mailto:binshao1118@gmail.com" target="_blank">binshao1118@gmail.com</a> <mailto:<a href="mailto:binshao1118@gmail.com" target="_blank">binshao1118@gmail.com</a>> <mailto:<a href="mailto:binshao1118@gmail.com" target="_blank">binshao1118@gmail.com</a> <mailto:<a href="mailto:binshao1118@gmail.com" target="_blank">binshao1118@gmail.com</a>>><br>
<br>
<br>
<br>
_______________________________________________<br>
Wien mailing list<br>
<a href="mailto:Wien@zeus.theochem.tuwien.ac.at" target="_blank">Wien@zeus.theochem.tuwien.ac.at</a> <mailto:<a href="mailto:Wien@zeus.theochem.tuwien.ac.at" target="_blank">Wien@zeus.theochem.tuwien.ac.at</a>><div class="im">
<br>
<a href="http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien" target="_blank">http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien</a><br>
<br>
<br>
--<br>
<br>
P.Blaha<br>
--------------------------------------------------------------------------<br>
Peter BLAHA, Inst.f. Materials Chemistry, TU Vienna, A-1060 Vienna<br>
Phone: +43-1-58801-15671 FAX: +43-1-58801-15698<br></div>
Email: <a href="mailto:blaha@theochem.tuwien.ac.at" target="_blank">blaha@theochem.tuwien.ac.at</a> <mailto:<a href="mailto:blaha@theochem.tuwien.ac.at" target="_blank">blaha@theochem.tuwien.ac.at</a>> WWW: <a href="http://info.tuwien.ac.at/theochem/" target="_blank">http://info.tuwien.ac.at/theochem/</a><br>
--------------------------------------------------------------------------<br>
_______________________________________________<br>
Wien mailing list<br>
<a href="mailto:Wien@zeus.theochem.tuwien.ac.at" target="_blank">Wien@zeus.theochem.tuwien.ac.at</a> <mailto:<a href="mailto:Wien@zeus.theochem.tuwien.ac.at" target="_blank">Wien@zeus.theochem.tuwien.ac.at</a>><div class="im">
<br>
<a href="http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien" target="_blank">http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien</a><br>
<br>
<br>
<br>
<br>
--<br>
Bin Shao, Ph.D. Candidate<br>
College of Information Technical Science, Nankai University<br>
94 Weijin Rd. Nankai Dist. Tianjin 300071, China<br>
Email: <a href="mailto:binshao1118@gmail.com" target="_blank">binshao1118@gmail.com</a> <mailto:<a href="mailto:binshao1118@gmail.com" target="_blank">binshao1118@gmail.com</a>><br>
<br>
<br>
<br>
_______________________________________________<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="http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien" target="_blank">http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien</a><br>
</div></blockquote><div><div></div><div class="h5">
<br>
-- <br>
<br>
P.Blaha<br>
--------------------------------------------------------------------------<br>
Peter BLAHA, Inst.f. Materials Chemistry, TU Vienna, A-1060 Vienna<br>
Phone: +43-1-58801-15671 FAX: +43-1-58801-15698<br>
Email: <a href="mailto:blaha@theochem.tuwien.ac.at" target="_blank">blaha@theochem.tuwien.ac.at</a> WWW: <a href="http://info.tuwien.ac.at/theochem/" target="_blank">http://info.tuwien.ac.at/theochem/</a><br>
--------------------------------------------------------------------------<br>
_______________________________________________<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="http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien" target="_blank">http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien</a><br>
</div></div></blockquote></div><br><br clear="all"><br>-- <br>Bin Shao, Ph.D. Candidate<br>College of Information Technical Science, Nankai University<br>94 Weijin Rd. Nankai Dist. Tianjin 300071, China<br>Email: <a href="mailto:binshao1118@gmail.com" target="_blank">binshao1118@gmail.com</a><br>