<div dir="ltr">Dear molpro users,<div><br></div><div style>I am trying to perform a constrained optimization with slapaf but keep getting an error message about the atom specifications. I have tried many types of specification and geometry input but they all result in the same message:</div>

<div style><br></div><div style>*****************</div><div style><div>PROGRAM * OPT (Geometry optimization)     Authors: F. Eckert and H.-J. Werner</div><div><br></div><div> ? Error</div><div> ? Wrong atoms specification</div>

<div> ? The problem occurs in CONS_PARSESTR</div><div>*****************</div><div><br></div><div style>The input I use is as follows:</div><div style><br></div><div style>**********************</div><div style><div>geomtyp=xyz</div>

<div>geometry={</div><div>   14</div><div> HF-SCF000/AVDZ  ENERGY=-339.52736086</div><div> N1         -0.6008305860        1.3734844428       -0.0000004728</div><div> N2          0.6008227340        1.3734877869       -0.0000006381</div>

<div> C1         -0.6599690213       -0.7035806493       -1.2331042232</div><div> C2         -1.2241519460       -0.0204813392       -0.0000001266</div><div> C3          1.2241520754       -0.0204745467       -0.0000004633</div>

<div> C4          0.6599716914       -0.7035769499       -1.2331050567</div><div> C5          0.6599730600       -0.7035736683        1.2331056606</div><div> C6         -0.6599676800       -0.7035771336        1.2331064822</div>

<div> H1          2.2999833641        0.0941966335       -0.0000006327</div><div> H2         -2.2999839058        0.0941835441        0.0000000000</div><div> H3         -1.2996459936       -1.0842250081       -2.0128910281</div>

<div> H4          1.2996497841       -1.0842183815       -2.0128923852</div><div> H5         -1.2996437212       -1.0842204495        2.0128945600</div><div> H6          1.2996521153       -1.0842138455        2.0128927979</div>

<div>}</div><div>basis=AVDZ</div><div><br></div><div>hf</div><div>{optg;</div><div>method,SLNR;</div><div>primitive,d1,dissoc,group1=[3,4,5,6,7,8,9,10,11,12,13,14],group2=[1,2];</div><div>constraint,1.0,d1;</div><div>}</div>

<div><br></div><div>***************************</div><div><br></div><div style>I also tried writing N1,N2 etc. in the group but that didn't change anything. </div><div style><br></div><div style>Any ideas?</div><div style>

<br></div><div style>Best regards,</div><div style>Shachar</div></div></div></div>