<p dir="ltr">Peter may be busy.....</p>
<p dir="ltr">The RKM value is critically important, and you cannot compare different values. I guess you are running on a single node without mpi, and the RKM value has been automatically reduced to avoid going beyond the memory limits (set via NMATMAX). Unless you have small minimum RMT's such as for H (~0.5) an RKM value around 3 is way too small. You need to repair this, either using more cores & mpi or something else.</p>
<p dir="ltr">You sound like you have some DFT experience, but I suspect not much with WIEN2K. Have you tried simpler systems first?</p>
<p dir="ltr">---<br>
Professor Laurence Marks<br>
Department of Materials Science and Engineering<br>
Northwestern University<br>
<a href="http://www.numis.northwestern.edu">http://www.numis.northwestern.edu</a><br>
Corrosion in 4D <a href="http://MURI4D.numis.northwestern.edu">http://MURI4D.numis.northwestern.edu</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</p>
<div class="gmail_quote">On Dec 1, 2015 01:40, "Bruno Landeros" <<a href="mailto:brunolanderos@hotmail.com">brunolanderos@hotmail.com</a>> wrote:<br type="attribution"><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<div dir="ltr">
<div dir="ltr" style="font-family:'Calibri','Segoe UI','Meiryo','Microsoft YaHei UI','Microsoft JhengHei UI','Malgun Gothic','sans-serif';font-size:12pt">
<div>Yes, there was a WARNING in the total energies.</div>
<div><br>
</div>
<div>RKM for the NOSYM is 3.18 and 3.57 for the symmetric.</div>
<div>Is this the reason? Should I increase it?<br>
</div>
<div>
<div><br>
</div>
<div>Enviado desde Correo de Windows</div>
<div><br>
</div>
</div>
<div style="padding-top:5px;border-top-color:rgb(229,229,229);border-top-width:1px;border-top-style:solid">
<div><font face=" 'Calibri', 'Segoe UI', 'Meiryo', 'Microsoft YaHei UI', 'Microsoft JhengHei UI', 'Malgun Gothic', 'sans-serif'" style="line-height:15pt;letter-spacing:0.02em;font-family:"Calibri","Segoe UI","Meiryo","Microsoft YaHei UI","Microsoft JhengHei UI","Malgun Gothic","sans-serif";font-size:12pt"><b>De:</b> <a href="mailto:pblaha@theochem.tuwien.ac.at" target="_blank">pblaha@theochem.tuwien.ac.at</a><br>
<b>Enviado el:</b> martes, 1 de diciembre de 2015 01:23 a. m.<br>
<b>Para:</b> <a href="mailto:wien@zeus.theochem.tuwien.ac.at" target="_blank">wien@zeus.theochem.tuwien.ac.at</a></font></div>
</div>
<div><br>
</div>
<div dir="">
<div>One more thought:<br>
<br>
Since these are large cells and one has inversion and the other not, <br>
could it be that due to size-limitations of NMATMAX the complex case (no <br>
symmetry) used a much lower RKMAX ? You also should see a WARNING in <br>
your total energies.<br>
<br>
grep :RKM case.scf for the two cases.<br>
<br>
On 12/01/2015 08:03 AM, Bruno Landeros wrote:<br>
> The first time I generated the Nosymm structure file was by using the 76<br>
> posiciones generated by a first symmetric scf calculation (no coordinate<br>
> optimization), so the same coordinates where used in both cases.<br>
><br>
> To be sure, I run again the Nosymm structure file but accepted the new<br>
> structfile suggested by the SGROUP program, which detected the corrected<br>
> symmetry and this gave the energy I just mentioned here, so in principle<br>
> they are equivalent.<br>
><br>
> Would you like me to send you the case.struct files with the<br>
> specifications I gave with init_lapw?<br>
><br>
> Enviado desde Correo de Windows<br>
><br>
> *De:* <a href="mailto:pblaha@theochem.tuwien.ac.at" target="_blank">pblaha@theochem.tuwien.ac.at</a> <<a href="mailto:pblaha@theochem.tuwien.ac.at" target="_blank">mailto:pblaha@theochem.tuwien.ac.at</a>><br>
> *Enviado el:* martes, 1 de diciembre de 2015 12:53 a. m.<br>
> *Para:* <a href="mailto:wien@zeus.theochem.tuwien.ac.at" target="_blank">wien@zeus.theochem.tuwien.ac.at</a><br>
> <<a href="mailto:wien@zeus.theochem.tuwien.ac.at" target="_blank">mailto:wien@zeus.theochem.tuwien.ac.at</a>><br>
><br>
> No. This is a huge E-difference which must come from something very<br>
> severe. Did you compare the distances in the 2 case.outputnn file or<br>
> compare the structures in xcrysden.<br>
><br>
> On 12/01/2015 07:37 AM, Bruno Landeros wrote:<br>
> > Dear Peter:<br>
> ><br>
> > Yes, my mistake. First energy (no symmetry) is<br>
> > -5752.55335845 Ry.<br>
> > while second energy is<br>
> > -5752.17735109 Ry.<br>
> ><br>
> > Since cell parameters are note small (12.994513, 25.002400, 17.381701<br>
> Bohr)<br>
> > and I was testing convergence I asked for just 1 K point in both cases.<br>
> > May be this the origin of the discrepancy?<br>
> ><br>
> > Greetings,<br>
> ><br>
> > Bruno<br>
> ><br>
> ><br>
> > Enviado desde Correo de Windows<br>
> ><br>
> > *De:* <a href="mailto:pblaha@theochem.tuwien.ac.at" target="_blank">pblaha@theochem.tuwien.ac.at</a> <<a href="mailto:pblaha@theochem.tuwien.ac.at" target="_blank">mailto:pblaha@theochem.tuwien.ac.at</a>><br>
> > *Enviado el:* martes, 1 de diciembre de 2015 12:16 a. m.<br>
> > *Para:* <a href="mailto:wien@zeus.theochem.tuwien.ac.at" target="_blank">wien@zeus.theochem.tuwien.ac.at</a><br>
> > <<a href="mailto:wien@zeus.theochem.tuwien.ac.at" target="_blank">mailto:wien@zeus.theochem.tuwien.ac.at</a>><br>
> ><br>
> > The energies you posted are identical !<br>
> ><br>
> > Anyway: I hope you did not just copy the case.klist file from the low to<br>
> > the high-symmetry case ??<br>
> ><br>
> > Otherwise: send me your 2 struct files together with the description of<br>
> > the chosen calculational parameters (everything which is non-default) to<br>
> > my private email.<br>
> ><br>
> > Am 01.12.2015 um 03:00 schrieb Bruno Landeros:<br>
> > > I'm studying a molecular crystal in which two different molecules form<br>
> > > an 1:1 adduct. There are 2 molecules of each kind (4 in total) in the<br>
> > > unit cell. There are 19 atoms in the asymmetric unit, and since it<br>
> > > belongs to the P21/c spacegroup, a total of 76 atoms are contained in<br>
> > > the unit cell. I was trying to find a possible mistake when<br>
> calculating<br>
> > > interaction energies and I found the next problem:<br>
> > ><br>
> > > If I calculate the SCF by manually inserting the 76 positions and<br>
> giving<br>
> > > it a P1 symmetry the energy converged to<br>
> > > -5752.17735109 Ry.<br>
> > ><br>
> > > If I calculate the SCF by inserting the 19 atoms of the asymmetric<br>
> unit<br>
> > > with the correct symmetry and let the program to generate the rest of<br>
> > > the positions (which where exactly the same as when inserted<br>
> manually),<br>
> > > the same functional, same RMT, same K points and same convergence<br>
> > > criteria the energy converged to<br>
> > > -5752.17735109 Ry.<br>
> > ><br>
> > > Transforming the energy difference between the 2 calculations<br>
> gives 195<br>
> > > kcal/mol. The same happened when doing the same for other systems.<br>
> > ><br>
> > > My question is, ¿shouldn't I get the same energy in both calculations<br>
> > > since it's the same system, same coordinates and the input was<br>
> prepared<br>
> > > equally except for the symmetry?<br>
> > > Since I'm trying to calculate interaction energies using different<br>
> > > subsystems, some with and some without symmetry, this is a relevant<br>
> > > question to my research.<br>
> > ><br>
> > ><br>
> > > Greetings,<br>
> > ><br>
> > > Bruno L<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>
> > > SEARCH the MAILING-LIST at:<br>
> ><br>
> <a href="http://www.mail-archive.com/wien%40zeus.theochem.tuwien.ac.at/index.html" target="_blank">
http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html</a><br>
> <<a href="http://www.mail-archive.com/wien%40zeus.theochem.tuwien.ac.at/index.html" target="_blank">http://www.mail-archive.com/wien%40zeus.theochem.tuwien.ac.at/index.html</a>><br>
> ><br>
> <<a href="http://www.mail-archive.com/wien%40zeus.theochem.tuwien.ac.at/index.html" target="_blank">http://www.mail-archive.com/wien%40zeus.theochem.tuwien.ac.at/index.html</a>><br>
> > ><br>
> ><br>
> > --<br>
> ><br>
> --------------------------------------------------------------------------<br>
> > Peter BLAHA, Inst.f. Materials Chemistry, TU Vienna, A-1060 Vienna<br>
> > Phone: +43-1-58801-165300 FAX: +43-1-58801-165982<br>
> > Email: <a href="mailto:blaha@theochem.tuwien.ac.at" target="_blank">blaha@theochem.tuwien.ac.at</a> WIEN2k: <a href="http://www.wien2k.at" target="_blank">
http://www.wien2k.at</a><br>
> > WWW: <a href="http://www.imc.tuwien.ac.at/staff/tc_group_e.php" target="_blank">
http://www.imc.tuwien.ac.at/staff/tc_group_e.php</a><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>
> > SEARCH the MAILING-LIST at:<br>
> ><br>
> <a href="http://www.mail-archive.com/wien%40zeus.theochem.tuwien.ac.at/index.html" target="_blank">
http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html</a><br>
> <<a href="http://www.mail-archive.com/wien%40zeus.theochem.tuwien.ac.at/index.html" target="_blank">http://www.mail-archive.com/wien%40zeus.theochem.tuwien.ac.at/index.html</a>><br>
> ><br>
> <<a href="http://www.mail-archive.com/wien%40zeus.theochem.tuwien.ac.at/index.html" target="_blank">http://www.mail-archive.com/wien%40zeus.theochem.tuwien.ac.at/index.html</a>><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>
> > SEARCH the MAILING-LIST at:<br>
> <a href="http://www.mail-archive.com/wien%40zeus.theochem.tuwien.ac.at/index.html" target="_blank">
http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html</a><br>
> <<a href="http://www.mail-archive.com/wien%40zeus.theochem.tuwien.ac.at/index.html" target="_blank">http://www.mail-archive.com/wien%40zeus.theochem.tuwien.ac.at/index.html</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-165300 FAX: +43-1-58801-165982<br>
> Email: <a href="mailto:blaha@theochem.tuwien.ac.at" target="_blank">blaha@theochem.tuwien.ac.at</a> WIEN2k: <a href="http://www.wien2k.at" target="_blank">
http://www.wien2k.at</a><br>
> WWW: <a href="http://www.imc.tuwien.ac.at/staff/tc_group_e.php" target="_blank">
http://www.imc.tuwien.ac.at/staff/tc_group_e.php</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>
> SEARCH the MAILING-LIST at:<br>
> <a href="http://www.mail-archive.com/wien%40zeus.theochem.tuwien.ac.at/index.html" target="_blank">
http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html</a><br>
> <<a href="http://www.mail-archive.com/wien%40zeus.theochem.tuwien.ac.at/index.html" target="_blank">http://www.mail-archive.com/wien%40zeus.theochem.tuwien.ac.at/index.html</a>><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>
> SEARCH the MAILING-LIST at: <a href="http://www.mail-archive.com/wien%40zeus.theochem.tuwien.ac.at/index.html" target="_blank">
http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html</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-165300 FAX: +43-1-58801-165982<br>
Email: <a href="mailto:blaha@theochem.tuwien.ac.at" target="_blank">blaha@theochem.tuwien.ac.at</a> WIEN2k: <a href="http://www.wien2k.at" target="_blank">
http://www.wien2k.at</a><br>
WWW: <a href="http://www.imc.tuwien.ac.at/staff/tc_group_e.php" target="_blank">
http://www.imc.tuwien.ac.at/staff/tc_group_e.php</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>
SEARCH the MAILING-LIST at: <a href="http://www.mail-archive.com/wien%40zeus.theochem.tuwien.ac.at/index.html" target="_blank">
http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html</a><br>
</div>
</div>
</div>
</div>
</blockquote></div>