<div dir="ltr">Hi,<div><br></div><div>individual total energy are almost always meaningless, especially when comparing different levels of theory. But, relative energies that correspond to physical quantities should give the same answer, if the level of theory is adequate.</div><div><br></div><div>Your input looks OK to me, except EPS_SCF 1.0E-4 is maybe too large, check that the energy is well converged.</div><div><br></div><div>Is your density of vacancies the same in each case? Did you relax the structures? Is LDA a good enough density functional. Is a DZVP basis set good enough?</div><div><br></div><div>Matt<br><br>On Thursday, June 1, 2017 at 4:31:36 PM UTC+1, Sebastian Hütter wrote:<blockquote class="gmail_quote" style="margin: 0;margin-left: 0.8ex;border-left: 1px #ccc solid;padding-left: 1ex;"><div dir="ltr">Hi all,<br><br>this may be a simple question, but somehow I feel like I am missing something very basic here.<br><br>I'm simply comparing bulk energies of metals (here: Al) computed over XYZ-periodic unit cells obtained via DFT in cp2k and several EAM potentials in LAMMPS. Different published potentials in LAMMPS all agree on an energy of somewhere around -14 eV/unit cell. However, when I run the same computation in cp2k (input file attached), I find energies (printed as the "FORCE_EVAL ( QS ) energy") of -8.3 Ha/unit cell, or -226 eV/unit cell. As this is more than one order of magnitude off, something can't be right here.<br>To see if the relative energies are usable, I then tried to calculate vacancy formation energy by comparing the supercell energies with and without vacancies (cf. <a href="https://icme.hpc.msstate.edu/mediawiki/index.php/LAMMPS_Vacancy_Formation_Energy" target="_blank" rel="nofollow" onmousedown="this.href='https://www.google.com/url?q\x3dhttps%3A%2F%2Ficme.hpc.msstate.edu%2Fmediawiki%2Findex.php%2FLAMMPS_Vacancy_Formation_Energy\x26sa\x3dD\x26sntz\x3d1\x26usg\x3dAFQjCNEmY2G2i453HWfJFU0yuzgA8rBmNA';return true;" onclick="this.href='https://www.google.com/url?q\x3dhttps%3A%2F%2Ficme.hpc.msstate.edu%2Fmediawiki%2Findex.php%2FLAMMPS_Vacancy_Formation_Energy\x26sa\x3dD\x26sntz\x3d1\x26usg\x3dAFQjCNEmY2G2i453HWfJFU0yuzgA8rBmNA';return true;">https://icme.hpc.msstate.edu/<wbr>mediawiki/index.php/LAMMPS_<wbr>Vacancy_Formation_Energy</a>) , but that is wrong as well, while the LAMMPS solution agrees reasonably well with experimental data.<br><br>It appears as if every energy value calculated is too large by a factor of 16? Why would that be the case?<br><br><br>Thank your for your help,<br><br>Sebastian<br></div></blockquote></div></div>