# [molpro-user] Molpro-user Digest, Vol 101, Issue 8

Nike Dattani 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:

==========================
geometry={angstrom;Li1;Li2,Li1,3.065}
basis = cc-pVDZ
wf,charge=0,spin=2,sym=1
occ,   3,0,0,0,1,0,0,0
closed,1,0,0,0,1,0,0,0
core
rhf
==========================

Gives: RHF STATE 1.1 Energy =  -12.635662909597

If I wrap up the RHF part in brackets:

==========================
{rhf
wf,charge=0,spin=2,sym=1
occ,3,0,0,0,1,0,0,0
closed,1,0,0,0,1,0,0,0
core
}
==========================

I get: RHF STATE 1.1 Energy =  -14.821433808867

Even if I wrap up everything but the "occ" part:

==========================
wf,charge=0,spin=2,sym=1
closed,1,0,0,0,1,0,0,0
core

{rhf
occ,3,0,0,0,1,0,0,0
}
==========================

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:

angstrom
geometry={Li1;Li2,Li1,3.065}
basis = aug-cc-pCV5Z
{rhf;occ,3,0,0,0,1,0,0,0
closed,1,0,0,0,1,0,0,0
wf,6,1,2
}

The occupations are correct, but the RHF energy is few cm-1 higher than
CFOUR,
and the CCSD(T) energy is 0.1 Hartree higher than CFOUR at aug-cc-pCV5Z
level:

RHF__CFOUR = -14.824,356
RHF_MOLPRO = -14.824,341

CCSD(T)__CFOUR = -14.91920357
CCSD(T)_MOLPRO = -14.82979310

Respectfully yours,
Nike Dattani
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://www.molpro.net/pipermail/molpro-user/attachments/20161221/94ca9066/attachment.html>
```