<html xmlns:v="urn:schemas-microsoft-com:vml" xmlns:o="urn:schemas-microsoft-com:office:office" xmlns:w="urn:schemas-microsoft-com:office:word" xmlns:m="http://schemas.microsoft.com/office/2004/12/omml" xmlns="http://www.w3.org/TR/REC-html40">
<head>
<meta http-equiv=Content-Type content="text/html; charset=us-ascii">
<meta name=Generator content="Microsoft Word 12 (filtered medium)">
<style>
<!--
/* Font Definitions */
@font-face
{font-family:Calibri;
panose-1:2 15 5 2 2 2 4 3 2 4;}
@font-face
{font-family:Consolas;
panose-1:2 11 6 9 2 2 4 3 2 4;}
/* Style Definitions */
p.MsoNormal, li.MsoNormal, div.MsoNormal
{margin:0in;
margin-bottom:.0001pt;
font-size:11.0pt;
font-family:"Calibri","sans-serif";}
a:link, span.MsoHyperlink
{mso-style-priority:99;
color:blue;
text-decoration:underline;}
a:visited, span.MsoHyperlinkFollowed
{mso-style-priority:99;
color:purple;
text-decoration:underline;}
p.MsoPlainText, li.MsoPlainText, div.MsoPlainText
{mso-style-priority:99;
mso-style-link:"Plain Text Char";
margin:0in;
margin-bottom:.0001pt;
font-size:10.5pt;
font-family:Consolas;}
span.EmailStyle17
{mso-style-type:personal-compose;
font-family:"Calibri","sans-serif";
color:windowtext;}
span.PlainTextChar
{mso-style-name:"Plain Text Char";
mso-style-priority:99;
mso-style-link:"Plain Text";
font-family:Consolas;}
.MsoChpDefault
{mso-style-type:export-only;}
@page Section1
{size:8.5in 11.0in;
margin:1.0in 1.0in 1.0in 1.0in;}
div.Section1
{page:Section1;}
-->
</style>
<!--[if gte mso 9]><xml>
<o:shapedefaults v:ext="edit" spidmax="1026" />
</xml><![endif]--><!--[if gte mso 9]><xml>
<o:shapelayout v:ext="edit">
<o:idmap v:ext="edit" data="1" />
</o:shapelayout></xml><![endif]-->
</head>
<body lang=EN-US link=blue vlink=purple>
<div class=Section1>
<p class=MsoNormal>Dear MOLPRO Users,<o:p></o:p></p>
<p class=MsoPlainText>in I have been having a problem with SP energy
calculations with<o:p></o:p></p>
<p class=MsoNormal>CCSD(T) for specific values of coordinates of the O2m(superoxide
radical) —H2 open-shell complex. Here is the input file that
I use for calculating the energy of the complex at various RCM values.
The Z-matrix is defined in the Jacobi coordinates, where RX3=the distance
between centers of masses of the H2 and O2m monomers. The other
coordinates are those obtained from full geometry optimization of the complex.
In order to generate one-d radial cut through the PES I need to vary RX3
from 1.9 -12 Ang. When RX3 greater than 2.3 Ang I get errors.
When I use Gaussian, I can scan the Jacobi radial coordinate from 1.9
Ang. It would be appreciated if anyone can help. Here is my
input file and the errors that I get<o:p></o:p></p>
<p class=MsoNormal><o:p> </o:p></p>
<p class=MsoNormal>***, O2m_H2 SP<o:p></o:p></p>
<p class=MsoNormal>Memory,90,m<o:p></o:p></p>
<p class=MsoNormal>print, basis, orbital<o:p></o:p></p>
<p class=MsoNormal>RX1=0.39036354 Ang<o:p></o:p></p>
<p class=MsoNormal>RX2=0.68220689 Ang<o:p></o:p></p>
<p class=MsoNormal>RX3= 2.30 Ang<o:p></o:p></p>
<p class=MsoNormal>AN1=105.11873358 deg<o:p></o:p></p>
<p class=MsoNormal>AN2=15.44582222 deg<o:p></o:p></p>
<p class=MsoNormal><o:p> </o:p></p>
<p class=MsoNormal>geometry<o:p></o:p></p>
<p class=MsoNormal>Ang<o:p></o:p></p>
<p class=MsoNormal>X1<o:p></o:p></p>
<p class=MsoNormal>X2 X1 1.0<o:p></o:p></p>
<p class=MsoNormal>H3 X2 RX1 X1 90.0<o:p></o:p></p>
<p class=MsoNormal>H4 X2 RX1 X1 90.0 H3 180.0<o:p></o:p></p>
<p class=MsoNormal>X5 X2 RX3 H3 AN2 X1 180.0<o:p></o:p></p>
<p class=MsoNormal>O6 X5 RX2 X2 AN1 H3 180.0<o:p></o:p></p>
<p class=MsoNormal>X7 X5 1.0 O6 90.0 X1 180.0<o:p></o:p></p>
<p class=MsoNormal>O8 X5 RX2 X7 90.0 O6 180.0<o:p></o:p></p>
<p class=MsoNormal>end<o:p></o:p></p>
<p class=MsoNormal><o:p> </o:p></p>
<p class=MsoNormal>basis=aug-cc-pvdz<o:p></o:p></p>
<p class=MsoNormal>gthresh,energy=1.d-9<o:p></o:p></p>
<p class=MsoNormal><o:p> </o:p></p>
<p class=MsoNormal>{hf,maxit=200<o:p></o:p></p>
<p class=MsoNormal>wf,19,2,1,-1<o:p></o:p></p>
<p class=MsoNormal>occ,8,2<o:p></o:p></p>
<p class=MsoNormal>open,2.2}<o:p></o:p></p>
<p class=MsoNormal>uccsd(t)<o:p></o:p></p>
<p class=MsoNormal>E_dim= energy<o:p></o:p></p>
<p class=MsoNormal>---<o:p></o:p></p>
<p class=MsoNormal>ERRORS for specific RX3 values<o:p></o:p></p>
<p class=MsoNormal>RX3=2.22<o:p></o:p></p>
<p class=MsoNormal>gthresh,energy=1.d-6<o:p></o:p></p>
<p class=MsoNormal>calculation is completed without error<o:p></o:p></p>
<p class=MsoNormal><o:p> </o:p></p>
<p class=MsoNormal>RX3=2.22<o:p></o:p></p>
<p class=MsoNormal>Gthresh,energy=1.d-9<o:p></o:p></p>
<p class=MsoNormal>?Error: RHF not converged. This error exit can be avoided
using the IGNORE_ERROR option on the ORBITAL directive<o:p></o:p></p>
<p class=MsoNormal> ERROR EXIT<o:p></o:p></p>
<p class=MsoNormal> CURRENT STACK: MAIN<o:p></o:p></p>
<p class=MsoNormal><o:p> </o:p></p>
<p class=MsoNormal>RX3=2.2<o:p></o:p></p>
<p class=MsoNormal>gthresh,energy=1.d-6<o:p></o:p></p>
<p class=MsoNormal>?Error: RHF not converged. This error exit can be avoided
using the IGNORE_ERROR option on the ORBITAL directive<o:p></o:p></p>
<p class=MsoNormal> ERROR EXIT<o:p></o:p></p>
<p class=MsoNormal> CURRENT STACK: MAIN<o:p></o:p></p>
<p class=MsoNormal><o:p> </o:p></p>
<p class=MsoNormal>RX3=2.3<o:p></o:p></p>
<p class=MsoNormal>Gthresh,energy=1.d-9<o:p></o:p></p>
<p class=MsoNormal>Calculation is completed without error<o:p></o:p></p>
<p class=MsoNormal><o:p> </o:p></p>
</div>
</body>
</html>