<div dir="ltr"><div>Werner,<br>Thank you for this suggestion. I should have mentioned that this was one of the permutations I tried. The problem is that after the first iteration optg defaults to the basis set defined in hessproc. This is most easily seen in the molpro.log file which shows the # of external orbitals for avtz as 33  21  21  10, and avdz is 17   9   9   4. This problem is what led me to try using proc and hessproc.<br>
</div>Thanks,<br>Benj<br></div><div class="gmail_extra"><br><br><div class="gmail_quote">On Wed, May 22, 2013 at 8:34 AM, Werner Gyorffy <span dir="ltr"><<a href="mailto:gyorffy@theochem.uni-stuttgart.de" target="_blank">gyorffy@theochem.uni-stuttgart.de</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Dear Benj,<br>
<br>
This is a slightly modified version of your input which does the job you wanted:<div class="im"><br>
<br>
***,h2o<br>
memory,150,M<br>
gthresh,oneint=1.d-14,twoint=<u></u>1.d-14,zero=1.d-14<br>
<br>
oc2    =     1.400000<br>
<br>
angstrom<br>
symmetry,auto<br>
geom={<br>
 c;<br>
 o,   1, oc2;<br>
}<br>
basis=avtz;<br>
<br></div>
hf;<br>
ccsd;<br>
{optg,energy=1.d-5,gradient=1.<u></u>d-5,numhess=5,hessproc=<u></u>b3lyphess}<br>
<br>
proc b3lyphess<br>
basis=avdz; ! for example<br>
hf;<br>
ks,b3lyp;<br>
endproc<br>
<br>
---<br>
<br>
Regards,<br>
<br>
Werner Gyorffy.<div><div class="h5"><br>
<br>
On 05/21/2013 11:38 PM, Benj FitzPatrick wrote:<br>
</div></div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div><div class="h5">
Hello all,<br>
I am trying to perform geometry optimizations using a hessian procedure<br>
with a smaller basis set, which then requires I use another procedure<br>
for the energy and gradient calculations (otherwise the hessian basis<br>
set becomes the default). The first point uses the analytic CCSD<br>
gradients, but after the hessian calculation numerical CCSD gradients<br>
are used. Have I specified one of my procedures incorrectly or is<br>
something more nefarious going on?<br>
Thanks,<br>
Benj FitzPatrick<br>
<br>
***,h2o<br>
memory,150,M<br>
gthresh,oneint=1.d-14,twoint=<u></u>1.d-14,zero=1.d-14<br>
<br>
oc2    =     1.400000<br>
<br>
angstrom<br>
symmetry,auto<br>
geom={<br>
  c;<br>
  o,   1, oc2;<br>
}<br>
basis=avtz;<br>
{optg,root=1,energy=1.d-5,<u></u>gradient=1.d-5,procedure=<u></u>ccsdavtz,numhess=5,hessproc=<u></u>b3lyphess}<br>
<br>
proc ccsdavtz<br>
hf;<br>
ccsd;<br>
endproc<br>
<br>
proc b3lyphess<br>
hf;<br>
ks,b3lyp;<br>
endproc<br>
<br>
<br>
<br></div></div>
______________________________<u></u>_________________<br>
Molpro-user mailing list<br>
<a href="mailto:Molpro-user@molpro.net" target="_blank">Molpro-user@molpro.net</a><br>
<a href="http://www.molpro.net/mailman/listinfo/molpro-user" target="_blank">http://www.molpro.net/mailman/<u></u>listinfo/molpro-user</a><br>
<br>
</blockquote>
______________________________<u></u>_________________<br>
Molpro-user mailing list<br>
<a href="mailto:Molpro-user@molpro.net" target="_blank">Molpro-user@molpro.net</a><br>
<a href="http://www.molpro.net/mailman/listinfo/molpro-user" target="_blank">http://www.molpro.net/mailman/<u></u>listinfo/molpro-user</a><br>
</blockquote></div><br></div>