<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>