<div dir="ltr">Dear Susi,<div><br></div><div>Thank you.</div><div>Yes, it is Fe2+, not Fe. However, since I put 4 negative point charges, the nuclear repulsive energy should be negative.</div><div>Prof. Fred Manby has answered my question. In Molpro, the interaction between point charges are ignored, as it is normally handled in an MM code, while other QM codes include this interaction.</div>

<div><br></div><div>The input I gave is just a toy model, as CASSCF will be done to check the wavefunction.</div><div><br></div><div>Best regards,</div><div>Quan </div><div><br></div></div><div class="gmail_extra"><br><br>

<div class="gmail_quote">On Tue, Apr 8, 2014 at 8:36 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 class="">On Mon, 7 Apr 2014 21:20:22 +0200<br>
Quan Phung <<a href="mailto:manhquank9@gmail.com">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>
</div>First, you seem to be calculating Fe^{2+}, not Fe. Second, the nuclear<br>
repulsion energy should be positive. Anyway...<br>
<br>
Transition metals are notoriously difficult for single-determinant<br>
calculations. The wave function space is laden with saddle point<br>
solutions, and it is extremely difficult to converge even to a local<br>
minimum. I remember trying to toy around with point charges to lift<br>
this degeneracy, but didn't get anywhere.<br>
<br>
You really need to run a stability analysis before you can trust the<br>
numbers you are getting.<br>
<span class="HOEnZb"><font color="#888888">--<br>
---------------------------------------------------------------<br>
Mr. Susi Lehtola, PhD             Research Associate<br>
<a href="mailto:susi.lehtola@alumni.helsinki.fi">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">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>
</font></span><div class="HOEnZb"><div class="h5">_______________________________________________<br>
Molpro-user mailing list<br>
<a href="mailto:Molpro-user@molpro.net">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>