<div dir="ltr">Ah, I see. Very helpful!<div><br></div><div>Thanks a bunch!</div><div><br></div><div>~Alan Chien†</div></div><div class="gmail_extra"><br><br><div class="gmail_quote">On Wed, Jan 22, 2014 at 4:14 AM, Grant Hill <span dir="ltr"><<a href="mailto:Grant.Hill@glasgow.ac.uk" target="_blank">Grant.Hill@glasgow.ac.uk</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Dear Alan,<br>
<br>
Iíve no experience with the Qchem program, but Iím guessing that it does not use the frozen core approximation by default. If I change your MP2 directive to {mp2;core,0} I get:<br>
<br>
!MP2 total energy † † † † † † † † † †-78.429489771679<br>
<br>
Best regards,<br>
<br>
Grant<br>
<div class="im HOEnZb"><br>
<br>
<br>
On 21 Jan 2014, at 20:35, Alan Chien <<a href="mailto:alandc@umich.edu">alandc@umich.edu</a>> wrote:<br>
<br>
</div><div class="HOEnZb"><div class="h5">> Hello all. I'm a new user of Molpro and was running a few little tests on my newly installed 2012.1.9 version of Molpro :D. However, I'm getting some weird results.<br>

><br>
> So, I'm comparing an MP2 single point energy calculation of ethylene between Molpro and Qchem.<br>
><br>
> I get the same HF energy, but differing MP2 energies at the same geometry.<br>
><br>
> Molpro: HF= -<a href="tel:78.06321916" value="+17806321916">78.06321916</a>, MP2= -78.39895488<br>
> Qchem: HF= -78.0632189754, MP2= -78.4294907701<br>
><br>
> I looked over my inputs (provided below) and they seem fineÖany suggestions?<br>
><br>
> I'll provide outputs if people want them.<br>
><br>
><br>
> QCHEM IMPUT<br>
> -----------------------------------------------------------<br>
> $molecule<br>
> 0 1<br>
> C † † † -2.0293377388 † † †1.5831581629 † † †0.0000000000<br>
> C † † † -0.6979276316 † † †1.6261615210 † † -0.0000000000<br>
> H † † † -2.5974389181 † † †2.2373976708 † † †0.6547741027<br>
> H † † † -2.5540361699 † † †0.8936222030 † † -0.6547741027<br>
> H † † † -0.1732292005 † † †2.3156974809 † † †0.6547741027<br>
> H † † † -0.1298264523 † † †0.9719220131 † † -0.6547741027<br>
> $end<br>
><br>
> $rem<br>
> Exchange HF<br>
> Correlation MP2<br>
><br>
> basis cc-pvtz<br>
><br>
> symmetry false<br>
> sym_ignore true<br>
> unrestricted false<br>
><br>
> jobtype sp<br>
> $end<br>
> -----------------------------------------------------------<br>
><br>
><br>
><br>
> MOLPRO INPUT<br>
> -----------------------------------------------------------<br>
> ***, Ethylene Geom Opt<br>
> symmetry=nosym<br>
> geometry={<br>
> C † † † -2.0293377388 † † †1.5831581629 † † †0.0000000000<br>
> C † † † -0.6979276316 † † †1.6261615210 † † -0.0000000000<br>
> H † † † -2.5974389181 † † †2.2373976708 † † †0.6547741027<br>
> H † † † -2.5540361699 † † †0.8936222030 † † -0.6547741027<br>
> H † † † -0.1732292005 † † †2.3156974809 † † †0.6547741027<br>
> H † † † -0.1298264523 † † †0.9719220131 † † -0.6547741027<br>
> }<br>
> basis=VTZ<br>
> hf<br>
> mp2<br>
> -----------------------------------------------------------<br>
</div></div><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>
<br>
</div></div></blockquote></div><br></div>