[molpro-user] Restrict keyword for pair excitation

[Panida Surawatanawong]

Dear Molpro users,

Recently, I have some problems with calculation MCSCF by adding keyword
RESTRICT. I tried to do PEMCSCF for BeH2. From the manual, I think if I put
keyword: restrict,-1,-1,orb1,orb2,.... The program should omit singly
excited configurations in that orbital. After I calculated it, I got the
same energy as the one without restrict keyword. The result energies should
be different, isn't it? Did I do anything wrong? My input is following.

***,BeH2PE
file,2,beh2pe.wfu
geometry={
Be;
H1,Be,r1;
H2,Be,r1,H1,theta;
}
r1=2.52955373
theta=180.00
basis=vdz

hf;
config;
wf,6,1,0;
orbprint,nvirt=3

multi;
occ,3,0,0,0,2,0,0,0;
core,1,0,0,0,0,0,0,0;
wf,6,1,0;
restrict,-1,-1,2.1,1.5,3.1,2.5;
print,ref

Best Regards,
Panida S.
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://www.molpro.net/pipermail/molpro-user/attachments/20050927/1363e0b5/attachment.html>