[molpro-user] Molpro-user Digest, Vol 101, Issue 8
dattani.nike at gmail.com
Wed Dec 21 23:06:46 CET 2016
Dear Jacky, and Kirk,
Both of your methods worked!
In summary, using MOLPRO 2012:
basis = cc-pVDZ
Gives: RHF STATE 1.1 Energy = -12.635662909597
If I wrap up the RHF part in brackets:
I get: RHF STATE 1.1 Energy = -14.821433808867
Even if I wrap up everything but the "occ" part:
I get the correct energy and occupations.
It appears that the SCF program simply does not remember the "occ" cards
defined before the RHF command, even though it DOES remember charge, spin,
sym, and closed.
Is this enough to be considered a bug?
Or was "occ" never designed to work the same way as "closed" ?
Does my first input work in MOLPRO 2015 ?
Another puzzle is that even with Kirk's input format:
basis = aug-cc-pCV5Z
The occupations are correct, but the RHF energy is few cm-1 higher than
and the CCSD(T) energy is 0.1 Hartree higher than CFOUR at aug-cc-pCV5Z
RHF__CFOUR = -14.824,356
RHF_MOLPRO = -14.824,341
CCSD(T)__CFOUR = -14.91920357
CCSD(T)_MOLPRO = -14.82979310
-------------- next part --------------
An HTML attachment was scrubbed...
More information about the Molpro-user