<div dir="ltr">Thank you so much, that makes a lot more sense (I'm very new to the program, especially variables).  I tried running this, but I ran into problems with the CCSD(T): it crashed, stating "?Error: RHF not converged.  This error exit can be avoided using the IGNORE_ERROR option on the ORBITAL directive," which I've run into a few times before, and I finally got it fixed just now!  Thank you so much for all your help, I really appreciate it!  I'll let you know if I run into any problems with this, but it appears to have begun going through each different step.</div>
<div class="gmail_extra"><br><br><div class="gmail_quote">On Mon, Jul 15, 2013 at 1:15 PM, Peterson, Kirk <span dir="ltr"><<a href="mailto:kipeters@wsu.edu" target="_blank">kipeters@wsu.edu</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">




<div style="word-wrap:break-word">
<div>Brandon,</div>
<div><br>
</div>
Try the following (gets rid of R1(i,j,k,l), saves you an integral evaluation per bond length, and simplifies your tables):
<div><br>
</div>
<div>
<div dir="ltr">****,I3<br>
<div>memory,50,M</div>
<div><br>
R1=[1.3 1.4 1.5 1.6 1.7 1.8 1.9 2.0 2.1 2.2 2.3 2.4 2.5 2.6 2.7 2.8 2.9 3.0 3.1 3.2 3.3 3.4 3.5] angstrom</div>
<div>R2=3.2446 angstrom</div>
<div>Theta=180.0<br>
<br>
geometry={<br>
</div>
<div>I1;<br>
</div>
<div>I2,I1,R1(i);<br>
I3,I2,R2,I1,Theta<br>
}</div>
<div><br>
DO i=1,#R1</div>
<div><br>
basis=vtz<br>
hf<br>
ehf(i)=energy<br>
ccsd(t)<br>
eccsdpt(i)=energy</div>
<div><br>
</div>
<div>casscf<br>
</div>
<div>ecas(i)=energy<br>
mrci<br>
emrci(i)=energy<br>
<br>
ENDDO<br>
table,ehf,eccsdpt,ecas,emrci</div>
<div>head,scf,ccsd(t),casscf,mrci<div class="im"><br>
title,Results for I3 3_2446 $basis,basis=$basis<br>
<br>
</div></div>
<div>---</div>
</div><div><div class="h5">
<div><br>
</div>
<br>
<div>
<div>On Jul 15, 2013, at 7:08 PM, Brandon Buncher <<a href="mailto:brandon.buncher.12@gmail.com" target="_blank">brandon.buncher.12@gmail.com</a>> wrote:</div>
<br>
<blockquote type="cite">
<div dir="ltr">
<div>I changed the input file to the following and still received the same error as before.  Also, as confirmation, would this create two tables: one with the HF and CCSDT results, and then one with the CASSCF and MRCI results?  Or am I setting up the tables
 incorrectly.  Thank you<br>
<br>
</div>
****,I3<br>
<div>memory,50,M<br>
R1=[1.3 1.4 1.5 1.6 1.7 1.8 1.9 2.0 2.1 2.2 2.3 2.4 2.5 2.6 2.7 2.8 2.9 3.0 3.1 3.2 3.3 3.4 3.5] ang, R2=3.2446 ang,Theta=180.0<br>
i=1<br>
geometry={<br>
</div>
<div>I1;<br>
</div>
<div>I2,I1,R1(i,j,k,l);<br>
I3,I2,R2,I1,Theta<br>
}<br>
do i=1,#R1<br>
basis=vtz<br>
hf<br>
e(i)=energy<br>
method(i)=program<br>
ccsd(t)<br>
e(i)=energy<br>
method(i)=program<br>
ENDDO<br>
table,method,e<br>
title,Results for I3 3_2446 $basis,basis=$basis<br>
<br>
do i=1,#R1<br>
hf<br>
casscf<br>
e(i)=energy<br>
method(i)=program<br>
mrci<br>
e(i)=energy<br>
method(i)=program<br>
ENDDO<br>
table,method,e<br>
title,Results for I3 3_2446 $basis,basis=$basis<br>
---<br>
<br>
</div>
</div>
<div class="gmail_extra"><br>
<br>
<div class="gmail_quote">On Mon, Jul 15, 2013 at 1:02 PM, Brandon Buncher <span dir="ltr">
<<a href="mailto:brandon.buncher.12@gmail.com" target="_blank">brandon.buncher.12@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">I used the SEWARD_MEM variable because my runs were terminating in errors, and after doing some research, I found out that I should use it, and it worked (URL below).  I will try a run with and without it to see if it will change anything.  With
 regards to the R1(I,J,K,L), I wasn't sure if having multiple variables named e(i) or method(i) would be a problem for the tables.  I will try changing what you have suggested and let you know what happens.  Thank you!
<div><br>
</div>
<div><a href="http://www.molpro.net/pipermail/molpro-user/2012-August/004997.html" target="_blank">http://www.molpro.net/pipermail/molpro-user/2012-August/004997.html</a><br>
</div>
<div><br>
</div>
</div>
<div class="gmail_extra">
<div>
<div><br>
<br>
<div class="gmail_quote">On Mon, Jul 15, 2013 at 12:36 PM, Peterson, Kirk <span dir="ltr">
<<a href="mailto:kipeters@wsu.edu" target="_blank">kipeters@wsu.edu</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
Brandon,<br>
<br>
two things I see right off:  The amount of memory you give on the main memory card, 20mw, is pretty small.  I would go with at least 50mw to start.  Also, why do you need the seward_mem variable?  I personally have never used this.  The main problem though
 is that Molpro's arrays are never larger than one-dimensional.  Thus in your z-matrix,  R1(i,j,k,l) makes no sense.  Why do you need different do loop counters when the arrays are all the same size, i.e., #R1 ?  I would also suggest to combine loops 1 and
 2 and combine loops 3 and 4 since in the first case you have to do the HF before the CCSD(T) anyway and in loop 4 presumably you should be doing a casscf before the mrci?<br>
<br>
best,<br>
<br>
-Kirk<br>
<div><br>
<br>
On Jul 15, 2013, at 5:35 PM, Brandon Buncher <<a href="mailto:brandon.buncher.12@gmail.com" target="_blank">brandon.buncher.12@gmail.com</a>> wrote:<br>
<br>
> Hi!  I've been trying to perform a variety of energy calculations on several different bond lengths of a molecule, beginning with HF and moving through CCSD(T), CASSCF, and MRCI.  However, whenever I try to do this, the calculation becomes stuck on "recomputing
 integrals since basis changed" right after setting all variables.  Is there a way to get around this?  I have tried changing the number of bond length options and removing calculation types, but they all become stuck here.  The input and output files are below.
  Thanks so much!<br>
><br>
> INPUT:<br>
> ****,I3<br>
> memory,20,M<br>
> set,seward_mem=3000000<br>
> R1=[1.3 1.4 1.5 1.6 1.7 1.8 1.9 2.0 2.1 2.2 2.3 2.4 2.5 2.6 2.7 2.8 2.9 3.0 3.1 3.2 3.3 3.4 3.5] ang, R2=3.2446 ang,Theta=180.0<br>
> i=1<br>
> j=1<br>
> k=1<br>
> l=1<br>
> geometry={<br>
> I1;<br>
> I2,I1,R1(i,j,k,l);<br>
> I3,I2,R2,I1,Theta<br>
> }<br>
> basis=vtz<br>
> do i=1,#R1<br>
> hf<br>
> e(i)=energy<br>
> method(i)=program<br>
> ENDDO<br>
> table,method,e<br>
> title,Results for I3 3_2446 $basis,basis=$basis<br>
><br>
> do j=1,#R1<br>
> hf<br>
> ccsd(t)<br>
> e(j)=energy<br>
> method(j)=program<br>
> ENDDO<br>
> table,method,e<br>
> title,Results for I3 3_2446 $basis,basis=$basis<br>
><br>
> do k=1,#R1<br>
> hf<br>
> casscf<br>
> e(k)=energy<br>
> method(k)=program<br>
> ENDDO<br>
> table,method,e<br>
> title,Results for I3 3_2446 $basis,basis=$basis<br>
><br>
> do l=1,#R1<br>
> hf<br>
> mrci<br>
> e(l)=energy<br>
> method(l)=program<br>
> ENDDO<br>
> table,method,e<br>
> title,Results for I3 3_2446 $basis,basis=$basis<br>
> ---<br>
><br>
><br>
><br>
> OUTPUT:<br>
>  Variables initialized (766), CPU time= 0.00 sec<br>
>  Commands  initialized (545), CPU time= 0.01 sec, 513 directives.<br>
>  Default parameters read. Elapsed time= 0.08 sec<br>
><br>
>  Checking input...<br>
>  Passed<br>
> 1<br>
>                                          ***  PROGRAM SYSTEM MOLPRO  ***<br>
>                          Copyright, University College Cardiff Consultants Limited, 2008<br>
><br>
>                                     Version 2012.1 linked 5 Sep 2012 15:26:37<br>
><br>
><br>
>  **********************************************************************************************************************************<br>
>  LABEL *   I3<br>
>  Linux-2.6.18-308.16.1.el5/molecular07.cluster(x86_64) 64 bit mpp version                DATE: 10-Jul-13          TIME: 16:41:43<br>
>  **********************************************************************************************************************************<br>
><br>
>  SHA1:             2c68d29c09da70e1723824271fadde4bcd5f07a0<br>
>  **********************************************************************************************************************************<br>
>  SETTING SEWARD_MEM     =     0.30000000D+07<br>
>  SETTING R1(1)          =         1.30000000  ANG<br>
>  SETTING R1(2)          =         1.40000000  ANG<br>
>  SETTING R1(3)          =         1.50000000  ANG<br>
>  SETTING R1(4)          =         1.60000000  ANG<br>
>  SETTING R1(5)          =         1.70000000  ANG<br>
>  SETTING R1(6)          =         1.80000000  ANG<br>
>  SETTING R1(7)          =         1.90000000  ANG<br>
>  SETTING R1(8)          =         2.00000000  ANG<br>
>  SETTING R1(9)          =         2.10000000  ANG<br>
>  SETTING R1(10)         =         2.20000000  ANG<br>
>  SETTING R1(11)         =         2.30000000  ANG<br>
>  SETTING R1(12)         =         2.40000000  ANG<br>
>  SETTING R1(13)         =         2.50000000  ANG<br>
>  SETTING R1(14)         =         2.60000000  ANG<br>
>  SETTING R1(15)         =         2.70000000  ANG<br>
>  SETTING R1(16)         =         2.80000000  ANG<br>
>  SETTING R1(17)         =         2.90000000  ANG<br>
>  SETTING R1(18)         =         3.00000000  ANG<br>
>  SETTING R1(19)         =         3.10000000  ANG<br>
>  SETTING R1(20)         =         3.20000000  ANG<br>
>  SETTING R1(21)         =         3.30000000  ANG<br>
>  SETTING R1(22)         =         3.40000000  ANG<br>
>  SETTING R1(23)         =         3.50000000  ANG<br>
>  SETTING R2             =         3.24460000  ANG<br>
>  SETTING THETA          =       180.00000000<br>
>  SETTING I              =         1.00000000<br>
>  SETTING J              =         1.00000000<br>
>  SETTING K              =         1.00000000<br>
>  SETTING L              =         1.00000000<br>
><br>
>  Variable memory set to   20000000 words,  buffer space   230000 words<br>
><br>
>  SETTING BASIS          =    VTZ<br>
><br>
>  DO I                   =         1.00000000<br>
><br>
><br>
>  Recomputing integrals since basis changed<br>
><br>
><br>
> --<br>
> Brandon Buncher<br>
</div>
> _______________________________________________<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/listinfo/molpro-user</a><br>
<br>
</blockquote>
</div>
<br>
<br clear="all">
<div><br>
</div>
</div>
</div>
<span><font color="#888888">-- <br>
Brandon Buncher </font></span></div>
</blockquote>
</div>
<br>
<br clear="all">
<br>
-- <br>
Brandon Buncher </div>
</blockquote>
</div>
<br>
</div></div></div>
</div>

</blockquote></div><br><br clear="all"><div><br></div>-- <br>Brandon Buncher
</div>