[molpro-user] RHF collapsed due to scfocc patch?
H. -J. Werner
werner at theochem.uni-stuttgart.de
Mon Jan 3 17:00:47 GMT 2005
Dear Prof. Hirano,
I have reproduced this problem and fixed it already.
Until there is a patch, you can use "closed" rather
than "open", i.e.:
rhf; occ,8,3,3,1;!open,8.1,3.2;
closed,7,2,3,1
I am sorry for the bug and the trouble.
Best regards
Joachim Werner
On Di, 04 Jan 2005, Tsuneo Hirano wrote:
>Dear Molpro users,
>
>After I applied the scfocc patch to the MOLPRO 2002.6, MOLPRO becomes to
>fail for an rhf open-shell job.
>
>Even if I specified the electron configuration by using 'open' command as
>--------------------------
> rhf; occ,8,3,3,1; open,8.1,3.2;
> sym,1,1,1,1,1,1,1,1,2,1;
> wf,28,2,2;
>--------------------------
>
>the output still complains as follows.
>-------------------------------------------------------------------
> .......
> Singly occupied orbitals: 8.1 3.2
>
> Extra symmetries of mos in symmetry 1: 1 1 1 1 1 1 1 2
> 1
>
> Orbital guess generated from atomic densities. Occupancy: 8 3 3
> 1
> Molecular orbital dump at record 2100.2
>
> ITERATION DDIFF GRAD ENERGY 2-EL.EN.
>DIPOLE MOMENTS DIIS
> 1 0.000D+00 0.000D+00 -1381.46180982 1049.816576
>0.000000 0.000000 -2.117146 0
> sorry, failure in assigning trial configuration
>
> assumed occupation
> 1- 4 8 3 3 1
>
> assumed closed shell orbitals
> 1- 4 7 2 3 1
> You should specify open shell configuration
> more precisely using CLOSED or OPEN commands
>
> ERROR EXIT
> CURRENT STACK: MAIN.
>---------------------------------------------------------------------
>My machine is operated by SunOS-5.8, and this failure did not happen
>before. I would like to know whether this is a machine dependent failure
>or a bug in the patch.
>
>Thank you in advance.
> Tsuneo Hirano
>
>P.S. My input file will follow.
>
>***,CoH as Co+ 3d8s0 (3F) + H-(1s2)(1S), Roos(Co) and aVQZ(H)
>memory,42,m;
>gprint,orbital;
>rcoh = 1.52 ang;
>
>basis={
>spdfg,Co,Bauschlicher-ANO;c;
>spdf,H,avqz;c;
>}
>
>geomtyp=xyz
>geometry={
>2
>CoH z-axis is the molecular axis
>Co, 0.0, 0.0, 0.0
>H, 0.0, 0.0, rcoh
>}
>
>!=== RHF initial MOs for Triplet-Phi (3d8s0, B1)
>text,rhf, to obtain ^3Phi (3d8s0, B1) initial orbitals
>rhf; occ,8,3,3,1;open,8.1,3.2;
>sym,1,1,1,1,1,1,1,1,2,1;
>wf,28,2,2;
>
>---
>
>
>
>
>**********************************
>Tsuneo Hirano
>
>Emeritus Professor of Chemistry, Ochanomizu University
>Visiting Professor at Hiroshima University
>
>e-mail: hirano at nccsk.net
>**********************************
>
--
Prof. Hans-Joachim Werner
Institute for Theoretical Chemistry
University of Stuttgart
Pfaffenwaldring 55
D-70569 Stuttgart, Germany
Tel.: (0049) 711 / 685 4400
Fax.: (0049) 711 / 685 4442
e-mail: werner at theochem.uni-stuttgart.de
More information about the Molpro-user
mailing list