<div dir="ltr">Dear Molpro users,<div><br></div><div>I'm trying to calculate Fe atom with 4 surrounding point charges.</div><div>The nuclear repulsion energy by Molpro is  -27.56070553 Ha.</div><div>However, both Orca and Molcas give -26.54614640 Ha.</div>


<div>Did I do something wrong? <br></div><div><br></div><div>Best regards,</div><div>Quan</div><div><br></div><div>Tthe Molpro input and the point charge files:<br></div><div>---------------------------------------------------------------------</div>


<div><div>Point charge in A</div><div>4</div><div>1.41198194, 1.41198194, 0.00000000, -1.0, 0</div><div>-1.41198194, 1.41198194, 0.00000000, -1.0, 0</div><div>1.41198194,-1.41198194, 0.00000000, -1.0, 0</div><div>-1.41198194,-1.41198194, 0.00000000, -1.0, 0</div>


<div><br></div></div><div>-------------------------------------------------------------------------</div><div><div>***,fe</div><div>memory,800,m</div><div>file,1,<a href="http://fe.int" target="_blank">fe.int</a>,new</div>

<div>
file,2,fe.wfu,new</div><div><br></div><div>SYMMETRY,NOSYM</div><div>orient,noorient</div><div>GPRINT,ORBITALS</div><div>geomtyp=xyz</div><div>geometry={</div><div>1</div><div>Fe</div><div>Fe     0.000000      0.000000    0.000000</div>


<div>}</div><div>LATTICE,INFILE=pc</div><div> </div><div>SET,DKROLL=1</div><div><br></div><div> basis=cc-pVDZ</div><div><br></div><div>{rhf</div><div>wf,24,1,4</div><div>}</div><div><br></div></div></div>