Dear Molpro Users,<div><br></div><div>I'm having some convergence issues when trying to apply a DIP or QUAD perturbation to my structures. The error arises when MOLPRO gets to performing MP2/CCSD etc calculations:</div>
<div><br></div><div><i>1PROGRAM * MP2 (Closed-shell)     Authors: C. Hampel, H.-J. Werner, 1991, M. Deegan, P.J. Knowles, 1992</i></div><div><i><br></i></div><div><i><br></i></div><div><i> Convergence thresholds:  THRVAR = 1.00D-08  THRDEN = 1.00D-06</i></div>
<div><i> ?NO</i></div><div><i><br></i></div><div><i> GLOBAL ERROR fehler on processor   0</i></div><div><br></div><div>I tried to solve the problem, but couldn't manage to overcome it or find anything in the literature. The input file I'm using is:</div>
<div><br></div><div><i>GEOMETRY={</i></div><div><i>6</i></div><div><i>Cl:H:NH3</i></div><div><i>Cl,0.0,0.0,0.913142</i></div><div><i>H,0.0, 0.595015,0.456571</i></div><div><i>N,0.0,0.0,0.0</i></div><div><i>H,0.0,0.947552,-0.380044</i></div>
<div><i>H,-0.820604,-0.473776,-0.380044</i></div><div><i>H,0.820604,-0.473776,-0.380044</i></div><div><i>}</i></div><div><i>BASIS,Cl=AVDZ,N=AVDZ,H=VDZ</i></div><div><i>HF</i></div><div><i>i=1</i></div><div><i>DIP+,0.0,0.0,0.001</i></div>
<div><i>{MP2</i></div><div><i>CORE,0}</i></div><div><i>           </i></div><div>Changing the energy thresholds (THRDEN/THRVAR and ONEINT) seemed to do nothing. Also, when I set the DIP+ perturbation to 0 but left it in the code (ie DIP+,0.0,0.0,0.0) the same problem occurs. This problem doesn't occur when I remove the DIP+ line from the code (ie there is no problem and everything works). The same deal when QUAD+ is used. This suggests there is something happening with the DIP/QUAD perturbations that it doesn't like. This structure is quite high in energy (compared to its global minimum). Any thoughts?</div>
<div><br></div><div>Kind Regards,</div><div>Michael Morris</div>