<div dir="ltr">Hello,<div><br></div><div>As I mentioned above, Molpro does not include the interaction between the point charges, which is 1.0145591 Ha.</div><div>The interaction between Fe (nucleus charge +26) and the point charge (-1) is <span style="font-size:13px;font-family:arial,sans-serif"><font color="#000000">-27.56070553 Ha, which is computed correctly by Molpro.</font></span></div>


<div><span style="font-size:13px;font-family:arial,sans-serif"><font color="#000000">So total = -26.54614640 Ha.</font></span></div><div><span style="font-size:13px;font-family:arial,sans-serif"><font color="#000000"><br>


</font></span></div><div><span style="font-size:13px;font-family:arial,sans-serif"><font color="#000000">Best regards,</font></span></div><div><span style="font-size:13px;font-family:arial,sans-serif"><font color="#000000">Quan</font></span></div>


<div><span style="font-size:13px;font-family:arial,sans-serif"><font color="#000000"><br></font></span></div><div class="gmail_extra"><br><br><div class="gmail_quote">On Tue, Apr 8, 2014 at 10:16 AM, Susi Lehtola <span dir="ltr"><<a href="mailto:susi.lehtola@alumni.helsinki.fi" target="_blank">susi.lehtola@alumni.helsinki.fi</a>></span> wrote:<br>


<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div>On Tue, 8 Apr 2014 09:36:09 +0300<br>
Susi Lehtola <<a href="mailto:susi.lehtola@alumni.helsinki.fi" target="_blank">susi.lehtola@alumni.helsinki.fi</a>> wrote:<br>
> On Mon, 7 Apr 2014 21:20:22 +0200<br>
> Quan Phung <<a href="mailto:manhquank9@gmail.com" target="_blank">manhquank9@gmail.com</a>> wrote:<br>
> > Dear Molpro users,<br>
> ><br>
> > I'm trying to calculate Fe atom with 4 surrounding point charges.<br>
> > The nuclear repulsion energy by Molpro is  -27.56070553 Ha.<br>
> > However, both Orca and Molcas give -26.54614640 Ha.<br>
> > Did I do something wrong?<br>
><br>
> First, you seem to be calculating Fe^{2+}, not Fe. Second, the nuclear<br>
> repulsion energy should be positive. Anyway...<br>
<br>
</div>The numbers do seem really weird. First, the nuclear repulsion energy<br>
doesn't depend on the electronic structure. Second, electronic<br>
contributions for iron would be much larger than a few dozen Hartree.<br>
But the number is really big to be purely Coulombic.<br>
<br>
Computing the numbers by hand I get<br>
<br>
- repulsion of the point charges: 1.0145 Ha<br>
- attraction of point charges and nucleus: -1.060 Ha<br>
total: -0.0455 Ha<br>
<br>
If you compute the difference between the Molpro and Orca numbers, you<br>
get the difference of 1.0145591 Ha, which you see is the repulsion of<br>
the point charges.<br>
<br>
I'm guessing Molpro doesn't include the static repulsion of the point<br>
charges in the energy it computes.<br>
<div><div>--<br>
---------------------------------------------------------------<br>
Mr. Susi Lehtola, PhD             Research Associate<br>
<a href="mailto:susi.lehtola@alumni.helsinki.fi" target="_blank">susi.lehtola@alumni.helsinki.fi</a>   Department of Applied Physics<br>
<a href="http://www.helsinki.fi/~jzlehtol" target="_blank">http://www.helsinki.fi/~jzlehtol</a>  Aalto University<br>
                                  Finland<br>
---------------------------------------------------------------<br>
Susi Lehtola, FT                  Tutkijatohtori<br>
<a href="mailto:susi.lehtola@alumni.helsinki.fi" target="_blank">susi.lehtola@alumni.helsinki.fi</a>   Fysiikan laitos<br>
<a href="http://www.helsinki.fi/~jzlehtol" target="_blank">http://www.helsinki.fi/~jzlehtol</a>  Aalto-yliopisto<br>
---------------------------------------------------------------<br>
_______________________________________________<br>
Molpro-user mailing list<br>
<a href="mailto:Molpro-user@molpro.net" target="_blank">Molpro-user@molpro.net</a><br>
<a href="http://www.molpro.net/mailman/listinfo/molpro-user" target="_blank">http://www.molpro.net/mailman/listinfo/molpro-user</a><br>
</div></div></blockquote></div><br></div></div>