Amrit,<br><br> My guess is that the accu is smaller than the default (accu=10), so probaly the defaults kick in. Can you try 11 or 12 and see what it does? Some of the thresholds can be changed by "gthresh". see the manual for different options..<br>
<br><div class="gmail_quote">On Sat, May 12, 2012 at 2:37 PM, Amrit Jalan <span dir="ltr"><<a href="mailto:amrit.jalan@gmail.com" target="_blank">amrit.jalan@gmail.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<div dir="ltr"><p class="MsoNormal">Hi,</p>
<p class="MsoNormal">I am trying to change the convergence thresholds in a
UHF-SCF calculation but haven’t had much luck. I tried the “accu” keyword after
the UHF directive but the program reverted to the default value. My input is
given below:</p>
<p class="MsoNormal"> </p>
<p class="MsoNormal"><span style="font-size:10pt;font-family:'Lucida Console';background-color:rgb(204,204,204)">geomtyp=xyz</span></p>
<p class="MsoNormal"><span style="font-size:10pt;font-family:'Lucida Console';background-color:rgb(204,204,204)">geometry=ts2.xyz</span></p>
<p class="MsoNormal"><span style="font-size:10pt;font-family:'Lucida Console';background-color:rgb(204,204,204)"> </span></p>
<p class="MsoNormal"><span style="font-size:10pt;font-family:'Lucida Console';background-color:rgb(204,204,204)">basis
= aug-cc-pVDZ</span></p>
<p class="MsoNormal"><span style="font-size:10.0pt;font-family:"Lucida Console""><span style="background-color:rgb(204,204,204)">uhf,accu=5</span></span></p>
<p class="MsoNormal"> </p>
<p class="MsoNormal">The output has this:</p>
<p class="MsoNormal"><span style="font-size:10.0pt;font-family:"Lucida Console""> </span></p>
<p class="MsoNormal"><span style="font-size:10pt;font-family:'Lucida Console';background-color:rgb(204,204,204)">1PROGRAM
* UHF-SCF Authors: W. Meyer, H.-J.
Werner</span></p>
<p class="MsoNormal"><span style="font-size:10pt;font-family:'Lucida Console';background-color:rgb(204,204,204)"> </span></p>
<p class="MsoNormal"><span style="font-size:10pt;font-family:'Lucida Console';background-color:rgb(204,204,204)"> </span></p>
<p class="MsoNormal"><span style="font-size:10pt;font-family:'Lucida Console';background-color:rgb(204,204,204)"> NUMBER OF ELECTRONS: 24+
24- SPACE SYMMETRY=1 SPIN SYMMETRY=Singlet</span></p>
<p class="MsoNormal"><span style="font-size:10pt;font-family:'Lucida Console';background-color:rgb(204,204,204)"> </span><span style="font-size:10pt;font-family:'Lucida Console'"><span style="background-color:rgb(255,255,102)">CONVERGENCE THRESHOLDS: 1.00E-06 (Density) 1.00E-07 (Energy)</span><span style="background-color:rgb(204,204,204)"></span></span></p>
<p class="MsoNormal"><span style="font-size:10pt;font-family:'Lucida Console';background-color:rgb(204,204,204)"> MAX. NUMBER OF ITERATIONS: 60</span></p>
<p class="MsoNormal"><span style="font-size:10pt;font-family:'Lucida Console';background-color:rgb(204,204,204)"> INTERPOLATION TYPE: DIIS</span></p>
<p class="MsoNormal"><span style="font-size:10pt;font-family:'Lucida Console';background-color:rgb(204,204,204)"> INTERPOLATION STEPS: 2 (START) 1 (STEP)</span></p>
<p class="MsoNormal"><span style="font-size:10.0pt;font-family:"Lucida Console""><span style="background-color:rgb(204,204,204)"> LEVEL SHIFTS: 0.30 (ALPHA) 0.30 (BETA)</span></span></p><p class="MsoNormal">
<font face="'Lucida Console'"><br></font></p><p class="MsoNormal"><font face="'Lucida Console'"><br></font></p>Any info on where the above convergence thresholds are coming from and how they can be changed would be very helpful. Also, I have noticed that the defaults change depending on the calculation being run. If this is so, is there someplace where this information is available?<p class="MsoNormal">
<font face="'Lucida Console'"><br></font></p><div><br></div>Thanks,<br>Amrit</div>
<br>_______________________________________________<br>
Molpro-user mailing list<br>
<a href="mailto:Molpro-user@molpro.net">Molpro-user@molpro.net</a><br>
<a href="http://www.molpro.net/mailman/listinfo/molpro-user" target="_blank">http://www.molpro.net/mailman/listinfo/molpro-user</a><br></blockquote></div><br><br clear="all"><br>-- <br>Regards,<br>Neeraj.<br><br><br>