<div dir="ltr">Typos in your energies?</div><div class="gmail_extra"><br><div class="gmail_quote">On Mon, Nov 30, 2015 at 8:00 PM, Bruno Landeros <span dir="ltr"><<a href="mailto:brunolanderos@hotmail.com" target="_blank">brunolanderos@hotmail.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 dir="ltr">I'm studying a molecular crystal in which two different molecules form an 1:1 adduct. There are 2 molecules of each kind (4 in total) in the unit cell. There are 19 atoms in the asymmetric unit, and since it belongs to the P21/c spacegroup, a
 total of 76 atoms are contained in the unit cell. I was trying to find a possible mistake when calculating interaction energies and I found the next problem:
<div><br>
</div>
<div>If I calculate the SCF by manually inserting the 76 positions and giving it a P1 symmetry the energy converged to</div>
<div>
<div>-5752.17735109 Ry.</div>
</div>
<div><br>
</div>
<div>If I calculate the SCF by inserting the 19 atoms of the asymmetric unit with the correct symmetry and let the program to generate the rest of the positions (which where exactly the same as when inserted manually), the same functional,  same RMT, same K
 points and same convergence criteria the energy converged to</div>
<div>
<div>-5752.17735109 Ry. </div>
</div>
<div><br>
</div>
<div>Transforming the energy difference between the 2 calculations gives 195 kcal/mol. The same happened when doing the same for other systems. </div>
<div><br>
</div>
<div>My question is, ¿shouldn't I get the same energy in both calculations since it's the same system, same coordinates and the input was prepared equally except for the symmetry? </div>
<div>Since I'm trying to calculate interaction energies using different subsystems, some with and some without symmetry, this is a relevant question to my research. </div>
<div><br>
</div>
<div><br>
</div>
<div>Greetings, </div>
<div><br>
</div>
<div>Bruno L</div>
</div>
</div>

</blockquote></div><br><br clear="all"><div><br></div>-- <br><div class="gmail_signature"><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://MURI4D.numis.northwestern.edu" target="_blank">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</div></div></div>
</div>