Hello Peter and David,<br><br>I had a similar problem, and exploring the code, I found a way to, somehow, go over it.<br><br>It's possible to restart an IRC with the keyword idir=0 instead of idir=+1/-1, within the optg command. It's not reasonable to determine a direction if you don't start at a saddle point, which is the case. Setting idir=0 makes it simply go downhill.<br>
<br>Well, it will not give exactly the same path as if the job hadn't died, since the algorithm that determines a step in the IRC uses some information (curvature and gradient, if I'm not mistaken) from a previous point. But from the second point on (of your restart), the algorithm is the same used from the third point on of the irc starting from a TS geometry.<br>
<br>In other words, despite the paths are not exactly the same (and it's not possible to really guarantee which one is the best), comparing the path obtained from a restart with the corresponding that would be got in the original job, the result is pretty much the same. In principle, the same minimum will be reached, and the path will be valid.<br>
<br>It must be noticed, however, that it seems that Molpro doesn't read hessians from a restart or external file. It implies that approximate hessian will be calculated for the first point in the restart, and henceforth the update scheme chosen will be applicable.<br>
<br>I recommend to check if the geometry read by the restart job (printed in the .out) matches the geometry of the last step of the killed job. Although I think it is, it's possible that the geometry read is the last one used in the gradient procedure, which is not the one desired.<br>
<br>Hope that helps.<br><br>Regards,<br>Gabriel Freitas<br><a href="mailto:gabrielnfreitas@gmail.com">gabrielnfreitas@gmail.com</a><br>Ph.D Student<br>UFRJ - Rio de Janeiro - Brazil<br><br><div class="gmail_quote">2010/12/29 Peter Ruprecht <span dir="ltr"><<a href="mailto:ruprech@jilau1.colorado.edu">ruprech@jilau1.colorado.edu</a>></span><br>
<blockquote class="gmail_quote" style="margin: 0pt 0pt 0pt 0.8ex; border-left: 1px solid rgb(204, 204, 204); padding-left: 1ex;">Dear Molpro community,<br>
<br>
(Disclaimer: I'm not much of a Molpro user, just the system administrator, so please take that into consideration if replying ;)<br>
<br>
One of my users had a long-running calculation that was killed after 100 days upon hitting our cluster's walltime limit. Since the job was very close to finishing at that point, we wanted to be able to restart it where it left off. So, the user resubmitted the identical input file to the same cluster node where the job had been running. On this node, the .wfu file and various other temp files from the original job were still waiting in the local $TMPDIR. Our understanding is that Molpro is smart enough to pick up again where the original calculation was killed.<br>
<br>
However, it looks like the calculation has started over from the beginning. Any ideas what we might have done wrong? The input file and output file from the attempted restart are attached. (Don't be confused by the fact that the input file has "restart" in it; so did the original job, which had been restarted in a different way from still another previous job.)<br>
<br>
If my explanation is not clear enough or you need more info, please let me know. Any help at all toward not losing 100 days of runtime will be gratefully appreciated!<br>
<br>
Thanks,<br><font color="#888888">
Peter Ruprecht<br>
JILA / University of Colorado<br>
</font><br><br>
Primary working directories : /data/nesbitt/djn<br>
Secondary working directories : /data/nesbitt/djn<br>
Wavefunction directory : /users/nesbitt/djn/wfu/<br>
Main file repository : /data/nesbitt/djn/<br>
<br>
cpu : P4 2003.000 MHz<br>
FC : /opt/intel/fce/10.1.008/bin/ifort<br>
FCVERSION : 10.1<br>
BLASLIB :<br>
id : colorado<br>
<br>
Nodes nprocs<br>
node14 8<br>
ga_uses_ma=false, calling ma_init with nominal heap.<br>
GA-space will be limited to 8.0 MW (determined by -G option)<br>
<br>
MPP tuning parameters: Latency= 0 Microseconds, Broadcast speed= 0 MB/sec<br>
default implementation of scratch files=sf<br>
<br>
<br>
memory,128,m<br>
<br>
!file,1,<a href="http://filename.int" target="_blank">filename.int</a> !allocate permanent integral file<br>
<br>
file,2,mal_vtzf12_irc_full_restart.wfu !allocate permanent wave-function (dump) file<br>
<br>
basis=vtz-f12<br>
<br>
geometry<br>
c1<br>
h2,1,rcht<br>
c3,1,rcclt,2,thcclt<br>
c4,1,rccrt,2,thccrt,3,180.<br>
o5,3,rcolt,1,tocclt,2,180.<br>
o6,4,rcort,1,toccrt,2,180.<br>
h7,3,rchlt,5,thcolt,1,180.<br>
h8,4,rchrt,6,thcort,1,180.<br>
x9,1,5.0,2,90.,3,-90.<br>
x10,1,rxct,9,90.,2,180.<br>
x11,10,5.0,1,90.,9,0.<br>
h12,10,rhxtunt,11,90.,1,dhxct<br>
end<br>
<br>
rxct = 1.5 ang<br>
<br>
readvar,mal_vtzf12_irc_full.act<br>
<br>
! start coordinates at last irc point<br>
<br>
i = 10<br>
<br>
rcclt = rccl(i)<br>
rchlt = rchl(i)<br>
rcolt = rcol(i)<br>
thcolt = thcol(i)<br>
thcclt = thccl(i)<br>
tocclt = toccl(i)<br>
<br>
rccrt = rccr(i)<br>
rchrt = rchr(i)<br>
rcort = rcor(i)<br>
thcort = thcor(i)<br>
thccrt = thccr(i)<br>
toccrt = toccr(i)<br>
<br>
rcht = rch(i)<br>
rhxtunt = rhxtun(i)<br>
dhxct = dhxc(i)<br>
<br>
rhf<br>
{ccsd(t)-f12b,ri_basis=cc-pvtz-f12/optri}<br>
<br>
{optg,stepmax=0.05,maxit=200,method=qsdpath,dir=1,numhess=5,hesscentral,saveact=mal_vtzf12_irc_full_restart;inactive,rxct} ! find IRC in positive direction<br>
<br>
table,irc,rccl,rchl,rcol,thcol,thccl,toccl,rccr,rchr,rcor,thcor,thccr,toccr,rch,rxc,rhxtun,dhxc,e_opt !tabulate results<br>
<br>
readvar,mal_vtzf12_irc_full_restart.act<br>
<br>
table,irc,rccl,rchl,rcol,thcol,thccl,toccl,rccr,rchr,rcor,thcor,thccr,toccr,rch,rxc,rhxtun,dhxc,e_opt !tabulate results<br>
<br>
<br>
<br>
Fortran compiler diagnostic<br>
***************************<br>
<br>
integer size 8<br>
double precision size 8<br>
real size 8<br>
logical size 8<br>
<br>
diagnostic completed successfully<br>
<br>
Variables initialized (631), CPU time= 0.00 sec<br>
Commands initialized (426), CPU time= 0.01 sec, 469 directives.<br>
Default parameters read. Elapsed time= 0.12 sec<br>
Checking input...<br>
Passed<br>
1<br>
<br>
<br>
*** PROGRAM SYSTEM MOLPRO ***<br>
Copyright, University College Cardiff Consultants Limited, 2008<br>
<br>
Version 2009.1 linked 7 Oct 2009 04:36:48<br>
<br>
<br>
**********************************************************************************************************************************<br>
LABEL *<br>
Linux-2.6.18-92.1.1.el5/node14(x86_64) 64 bit mpp version DATE: 29-Dec-10 TIME: 05:54:20<br>
**********************************************************************************************************************************<br>
<br>
Patch level: 20<br>
**********************************************************************************************************************************<br>
<br>
Variable memory set to 128000000 words, buffer space 230000 words<br>
<br>
<br>
<br>
Permanent file 2 mal_vtzf12_irc_full_restart.wfu assigned. Implementation=df Size= 3.74 MB<br>
<br>
<br>
PROGRAM * RESTART<br>
<br>
<br>
Reading variables from file 2<br>
<br>
_NUMVAR = 700.00000000<br>
_DKROLL_DONE = 0.00000000<br>
_FOCKDONE = 1.00000000<br>
_INTDONE = 1.00000000<br>
_SCFDONE = 1.00000000<br>
_DIRECT = 0.00000000<br>
_BASIS = VTZ-F12<br>
_NELEC = 38.00000000<br>
_OPTCONV = 0.01256362<br>
_QSDSTEP = 0.01256362<br>
_QSDIRC = 0.52406178<br>
_PROGRAM = RHF-SCF<br>
_CPUSTEP = 188.26000000<br>
_SYSSTEP = 12.94000000<br>
_WALLSTEP = 309.28000000<br>
_BASINP = 610.00000000<br>
_CPUTOT = 0.71492238D+07 SEC<br>
_DMX = 1.21705632<br>
_DMZ = -0.22020963<br>
_EMP2_SING = -0.61445288<br>
_EMP2_TRIP = -0.37099670<br>
_EMP2_SINGLES = 0.00000000<br>
_EMP2_SCS = -266.72370242<br>
_ECSING = -0.65159085<br>
_ECTRIP = -0.33008902<br>
_EF12_SING = -0.07462324<br>
_EF12_TRIP = -0.00896734<br>
_EF12_RHFRELAX = -0.00232285<br>
_EF12_SINGLES = -0.00232285<br>
_EF12_SCS = -0.09512757<br>
_EF12_STRONG = -0.08359058<br>
_EF12_CLOSE = 0.00000000<br>
_EF12_WEAK = 0.00000000<br>
_EF12_DIST = 0.00000000<br>
_EF12S_SING = -0.07468898<br>
_EF12S_TRIP = -0.00865961<br>
_EF12S_SCS = -0.09501497<br>
_EF12S_STRONG = -0.08334859<br>
_EF12S_CLOSE = 0.00000000<br>
_EF12S_WEAK = 0.00000000<br>
_EF12S_DIST = 0.00000000<br>
_EF12D_SING = -0.07918959<br>
_EF12D_TRIP = -0.00928509<br>
_EF12D_SCS = -0.10080490<br>
_EF12D_STRONG = -0.08847468<br>
_EF12D_CLOSE = 0.00000000<br>
_EF12D_WEAK = 0.00000000<br>
_EF12D_DIST = 0.00000000<br>
_ENERG_CC = 0.00000000<br>
_ENERG_CV = 0.00000000<br>
_ENERG_VV = -1.04968965<br>
_ENERGY = -265.75102472<br>
_ENUC = 164.78978903<br>
_GRADEXFAC = 9999.00000000<br>
_GRADMETHOD = 1.00000000<br>
_ITERATIONS = 11.00000000<br>
_LASTNELEC = 38.00000000<br>
_LASTSPIN = 0.00000000<br>
_LASTSYM = 1.00000000<br>
_OPTCONVMX = 0.00866340<br>
_OPTCONVRMS = 0.00083757<br>
_OPTGRAD = 0.00139213<br>
_OPTGRADMX = 0.00090381<br>
_OPTGRADRMS = 0.00009281<br>
_ORBITAL = 2100.10000000<br>
_EPDIAG = -0.14617726<br>
_SEIG_MIN = 0.00004183<br>
_SEWPROP = 1.00000000<br>
_STATUS = 1.00000000<br>
_SYSTOT = 0.23005290D+06 SEC<br>
_T1DIAG = 0.01726582<br>
_D1DIAG = 0.06630695<br>
_WALLTOT = 0.86391642D+07 SEC<br>
!DFBASIS_MP2 = AUG-CC-PVTZ-MP2F<br>
!DFBASIS_F12K = CC-PVTZ-JKFIT<br>
!RIBASIS_MP2 = CC-PVTZ-F12_OPT<br>
!RIBASIS_EXCH =<br>
_DATE = 15-Sep-10<br>
_LASTORB = RHF<br>
_PGROUP = Cs<br>
_TIME = 22:53:22<br>
_SVDDEL = 0.00000000<br>
RXCT = 1.50000000<br>
E_OPT(1:10) = -266.84674708 -266.84692062 -266.84717443 -266.84752392 -266.84795564 -266.84843843<br>
-266.84893857 -266.84942034 -266.84986487 -266.85026250<br>
IRC(1:10) = 0.00320006 0.02210900 0.03497764 0.04809850 0.06173091 0.07584387<br>
0.09016285 0.10443742 0.11867072 0.13282827<br>
RCH(1:10) = 1.07755169 1.07755294 1.07756399 1.07757970 1.07760243 1.07763558<br>
1.07768201 1.07774421 1.07782572 1.07792782<br>
RCCL(1:10) = 1.39710342 1.39939918 1.40105897 1.40280543 1.40469641 1.40675103<br>
1.40896484 1.41133402 1.41384800 1.41647919<br>
THCCL(1:10) = 121.81106172 121.25806403 120.85029858 120.42013914 119.96653450 119.48568208<br>
119.00791089 118.54013563 118.09917602 117.69659884<br>
RCCR(1:10) = 1.39636244 1.39407100 1.39249157 1.39087129 1.38913246 1.38725457<br>
1.38523430 1.38304495 1.38068707 1.37817260<br>
THCCR(1:10) = 121.98474342 122.52897075 122.84891143 123.14971423 123.43549837 123.69211078<br>
123.90087779 124.06700155 124.18246533 124.25894610<br>
RCOL(1:10) = 1.27667505 1.27432879 1.27263822 1.27086285 1.26892336 1.26679256<br>
1.26445941 1.26193086 1.25923877 1.25645120<br>
TOCCL(1:10) = 121.71795942 121.97438788 122.16475258 122.36465740 122.57350299 122.79083841<br>
123.00300733 123.20460550 123.38844280 123.54974039<br>
RCOR(1:10) = 1.27743854 1.27978209 1.28144052 1.28316757 1.28504195 1.28707239<br>
1.28928867 1.29168692 1.29425249 1.29694643<br>
TOCCR(1:10) = 121.63571532 121.38253679 121.23730102 121.10410036 120.98185382 120.87589228<br>
120.80060798 120.75262465 120.73614794 120.74439212<br>
RCHL(1:10) = 1.09059617 1.09089483 1.09111309 1.09134533 1.09160484 1.09189829<br>
1.09223468 1.09261917 1.09305686 1.09354699<br>
THCOL(1:10) = 117.00864432 117.05415013 117.08726895 117.12246978 117.16182094 117.20650632<br>
117.25813064 117.31811039 117.38797917 117.46884981<br>
RCHR(1:10) = 1.09049923 1.09020079 1.08998387 1.08975405 1.08949879 1.08921346<br>
1.08889026 1.08852716 1.08812303 1.08768286<br>
THCOR(1:10) = 116.99372312 116.94814255 116.91430019 116.87804148 116.83704376 116.79028436<br>
116.73542215 116.67091528 116.59476305 116.50554598<br>
RHXTUN(1:10) = 0.75289829 0.75290545 0.75308513 0.75333831 0.75369158 0.75333818<br>
0.75451309 0.75527174 0.75634104 0.75768516<br>
DHXC(1:10) = 179.91559673 179.48732667 179.19296076 178.89475927 178.58622554 178.27267758<br>
177.95737205 177.64535148 177.33958394 177.03685172<br>
I = 10.00000000<br>
RCCLT = 1.44141069<br>
RCHLT = 1.09951248<br>
RCOLT = 1.23562213<br>
THCOLT = 119.09765026<br>
THCCLT = 120.67309570<br>
TOCCLT = 123.26567742<br>
RCCRT = 1.35930678<br>
RCHRT = 1.08372450<br>
RCORT = 1.32262356<br>
THCORT = 113.68779853<br>
THCCRT = 118.39512575<br>
TOCCRT = 123.83006367<br>
RCHT = 1.07935371<br>
RHXTUNT = 0.89985896<br>
DHXCT = 155.75087965<br>
<br>
Geometry written to block 41 of record 700<br>
<br>
<br>
**********************************************************************************************************************************<br>
DATASETS * FILE NREC LENGTH (MB) RECORD NAMES<br>
2 10 3.93 500 610 700 1000 2100 7360 7350 5300 5600 5500<br>
VAR BASINP GEOM BASIS RHF F12ABS EF12 HESS GINFO HINFO<br>
<br>
PROGRAMS * TOTAL RESTART<br>
CPU TIMES * 0.20 0.08<br>
REAL TIME * 0.49 SEC<br>
DISK USED * 0.26 MB<br>
**********************************************************************************************************************************<br>
SETTING BASIS = VTZ-F12<br>
SETTING RXCT = 1.50000000 ANG<br>
<br>
Reading variables from file mal_vtzf12_irc_full.act<br>
SETTING E_OPT(1) = -266.84674708 AU<br>
SETTING IRC(1) = 0.00320006 AU<br>
SETTING RCH(1) = 1.07755169 ANG<br>
SETTING RCCL(1) = 1.39710342 ANG<br>
SETTING THCCL(1) = 121.81106172 DEG<br>
SETTING RCCR(1) = 1.39636244 ANG<br>
SETTING THCCR(1) = 121.98474342 DEG<br>
SETTING RCOL(1) = 1.27667505 ANG<br>
SETTING TOCCL(1) = 121.71795942 DEG<br>
SETTING RCOR(1) = 1.27743854 ANG<br>
SETTING TOCCR(1) = 121.63571532 DEG<br>
SETTING RCHL(1) = 1.09059617 ANG<br>
SETTING THCOL(1) = 117.00864432 DEG<br>
SETTING RCHR(1) = 1.09049923 ANG<br>
SETTING THCOR(1) = 116.99372312 DEG<br>
SETTING RHXTUN(1) = 0.75289829 ANG<br>
SETTING DHXC(1) = 179.91559673 DEG<br>
SETTING E_OPT(2) = -266.84692062 AU<br>
SETTING IRC(2) = 0.02210900 AU<br>
SETTING RCH(2) = 1.07755294 ANG<br>
SETTING RCCL(2) = 1.39939918 ANG<br>
SETTING THCCL(2) = 121.25806403 DEG<br>
SETTING RCCR(2) = 1.39407100 ANG<br>
SETTING THCCR(2) = 122.52897075 DEG<br>
SETTING RCOL(2) = 1.27432879 ANG<br>
SETTING TOCCL(2) = 121.97438788 DEG<br>
SETTING RCOR(2) = 1.27978209 ANG<br>
SETTING TOCCR(2) = 121.38253679 DEG<br>
SETTING RCHL(2) = 1.09089483 ANG<br>
SETTING THCOL(2) = 117.05415013 DEG<br>
SETTING RCHR(2) = 1.09020079 ANG<br>
SETTING THCOR(2) = 116.94814255 DEG<br>
SETTING RHXTUN(2) = 0.75290545 ANG<br>
SETTING DHXC(2) = 179.48732667 DEG<br>
SETTING E_OPT(3) = -266.84717443 AU<br>
SETTING IRC(3) = 0.03497764 AU<br>
SETTING RCH(3) = 1.07756399 ANG<br>
SETTING RCCL(3) = 1.40105897 ANG<br>
SETTING THCCL(3) = 120.85029858 DEG<br>
SETTING RCCR(3) = 1.39249157 ANG<br>
SETTING THCCR(3) = 122.84891143 DEG<br>
SETTING RCOL(3) = 1.27263822 ANG<br>
SETTING TOCCL(3) = 122.16475258 DEG<br>
SETTING RCOR(3) = 1.28144052 ANG<br>
SETTING TOCCR(3) = 121.23730102 DEG<br>
SETTING RCHL(3) = 1.09111309 ANG<br>
SETTING THCOL(3) = 117.08726895 DEG<br>
SETTING RCHR(3) = 1.08998387 ANG<br>
SETTING THCOR(3) = 116.91430019 DEG<br>
SETTING RHXTUN(3) = 0.75308513 ANG<br>
SETTING DHXC(3) = 179.19296076 DEG<br>
SETTING E_OPT(4) = -266.84752392 AU<br>
SETTING IRC(4) = 0.04809850 AU<br>
SETTING RCH(4) = 1.07757970 ANG<br>
SETTING RCCL(4) = 1.40280543 ANG<br>
SETTING THCCL(4) = 120.42013914 DEG<br>
SETTING RCCR(4) = 1.39087129 ANG<br>
SETTING THCCR(4) = 123.14971423 DEG<br>
SETTING RCOL(4) = 1.27086285 ANG<br>
SETTING TOCCL(4) = 122.36465740 DEG<br>
SETTING RCOR(4) = 1.28316757 ANG<br>
SETTING TOCCR(4) = 121.10410036 DEG<br>
SETTING RCHL(4) = 1.09134533 ANG<br>
SETTING THCOL(4) = 117.12246978 DEG<br>
SETTING RCHR(4) = 1.08975405 ANG<br>
SETTING THCOR(4) = 116.87804148 DEG<br>
SETTING RHXTUN(4) = 0.75333831 ANG<br>
SETTING DHXC(4) = 178.89475927 DEG<br>
SETTING E_OPT(5) = -266.84795564 AU<br>
SETTING IRC(5) = 0.06173091 AU<br>
SETTING RCH(5) = 1.07760243 ANG<br>
SETTING RCCL(5) = 1.40469641 ANG<br>
SETTING THCCL(5) = 119.96653450 DEG<br>
SETTING RCCR(5) = 1.38913246 ANG<br>
SETTING THCCR(5) = 123.43549837 DEG<br>
SETTING RCOL(5) = 1.26892336 ANG<br>
SETTING TOCCL(5) = 122.57350299 DEG<br>
SETTING RCOR(5) = 1.28504195 ANG<br>
SETTING TOCCR(5) = 120.98185382 DEG<br>
SETTING RCHL(5) = 1.09160484 ANG<br>
SETTING THCOL(5) = 117.16182094 DEG<br>
SETTING RCHR(5) = 1.08949879 ANG<br>
SETTING THCOR(5) = 116.83704376 DEG<br>
SETTING RHXTUN(5) = 0.75369158 ANG<br>
SETTING DHXC(5) = 178.58622554 DEG<br>
SETTING E_OPT(6) = -266.84843843 AU<br>
SETTING IRC(6) = 0.07584387 AU<br>
SETTING RCH(6) = 1.07763558 ANG<br>
SETTING RCCL(6) = 1.40675103 ANG<br>
SETTING THCCL(6) = 119.48568208 DEG<br>
SETTING RCCR(6) = 1.38725457 ANG<br>
SETTING THCCR(6) = 123.69211078 DEG<br>
SETTING RCOL(6) = 1.26679256 ANG<br>
SETTING TOCCL(6) = 122.79083841 DEG<br>
SETTING RCOR(6) = 1.28707239 ANG<br>
SETTING TOCCR(6) = 120.87589228 DEG<br>
SETTING RCHL(6) = 1.09189829 ANG<br>
SETTING THCOL(6) = 117.20650632 DEG<br>
SETTING RCHR(6) = 1.08921346 ANG<br>
SETTING THCOR(6) = 116.79028436 DEG<br>
SETTING RHXTUN(6) = 0.75333818 ANG<br>
SETTING DHXC(6) = 178.27267758 DEG<br>
SETTING E_OPT(7) = -266.84893857 AU<br>
SETTING IRC(7) = 0.09016285 AU<br>
SETTING RCH(7) = 1.07768201 ANG<br>
SETTING RCCL(7) = 1.40896484 ANG<br>
SETTING THCCL(7) = 119.00791089 DEG<br>
SETTING RCCR(7) = 1.38523430 ANG<br>
SETTING THCCR(7) = 123.90087779 DEG<br>
SETTING RCOL(7) = 1.26445941 ANG<br>
SETTING TOCCL(7) = 123.00300733 DEG<br>
SETTING RCOR(7) = 1.28928867 ANG<br>
SETTING TOCCR(7) = 120.80060798 DEG<br>
SETTING RCHL(7) = 1.09223468 ANG<br>
SETTING THCOL(7) = 117.25813064 DEG<br>
SETTING RCHR(7) = 1.08889026 ANG<br>
SETTING THCOR(7) = 116.73542215 DEG<br>
SETTING RHXTUN(7) = 0.75451309 ANG<br>
SETTING DHXC(7) = 177.95737205 DEG<br>
SETTING E_OPT(8) = -266.84942034 AU<br>
SETTING IRC(8) = 0.10443742 AU<br>
SETTING RCH(8) = 1.07774421 ANG<br>
SETTING RCCL(8) = 1.41133402 ANG<br>
SETTING THCCL(8) = 118.54013563 DEG<br>
SETTING RCCR(8) = 1.38304495 ANG<br>
SETTING THCCR(8) = 124.06700155 DEG<br>
SETTING RCOL(8) = 1.26193086 ANG<br>
SETTING TOCCL(8) = 123.20460550 DEG<br>
SETTING RCOR(8) = 1.29168692 ANG<br>
SETTING TOCCR(8) = 120.75262465 DEG<br>
SETTING RCHL(8) = 1.09261917 ANG<br>
SETTING THCOL(8) = 117.31811039 DEG<br>
SETTING RCHR(8) = 1.08852716 ANG<br>
SETTING THCOR(8) = 116.67091528 DEG<br>
SETTING RHXTUN(8) = 0.75527174 ANG<br>
SETTING DHXC(8) = 177.64535148 DEG<br>
SETTING E_OPT(9) = -266.84986487 AU<br>
SETTING IRC(9) = 0.11867072 AU<br>
SETTING RCH(9) = 1.07782572 ANG<br>
SETTING RCCL(9) = 1.41384800 ANG<br>
SETTING THCCL(9) = 118.09917602 DEG<br>
SETTING RCCR(9) = 1.38068707 ANG<br>
SETTING THCCR(9) = 124.18246533 DEG<br>
SETTING RCOL(9) = 1.25923877 ANG<br>
SETTING TOCCL(9) = 123.38844280 DEG<br>
SETTING RCOR(9) = 1.29425249 ANG<br>
SETTING TOCCR(9) = 120.73614794 DEG<br>
SETTING RCHL(9) = 1.09305686 ANG<br>
SETTING THCOL(9) = 117.38797917 DEG<br>
SETTING RCHR(9) = 1.08812303 ANG<br>
SETTING THCOR(9) = 116.59476305 DEG<br>
SETTING RHXTUN(9) = 0.75634104 ANG<br>
SETTING DHXC(9) = 177.33958394 DEG<br>
SETTING E_OPT(10) = -266.85026250 AU<br>
SETTING IRC(10) = 0.13282827 AU<br>
SETTING RCH(10) = 1.07792782 ANG<br>
SETTING RCCL(10) = 1.41647919 ANG<br>
SETTING THCCL(10) = 117.69659884 DEG<br>
SETTING RCCR(10) = 1.37817260 ANG<br>
SETTING THCCR(10) = 124.25894610 DEG<br>
SETTING RCOL(10) = 1.25645120 ANG<br>
SETTING TOCCL(10) = 123.54974039 DEG<br>
SETTING RCOR(10) = 1.29694643 ANG<br>
SETTING TOCCR(10) = 120.74439212 DEG<br>
SETTING RCHL(10) = 1.09354699 ANG<br>
SETTING THCOL(10) = 117.46884981 DEG<br>
SETTING RCHR(10) = 1.08768286 ANG<br>
SETTING THCOR(10) = 116.50554598 DEG<br>
SETTING RHXTUN(10) = 0.75768516 ANG<br>
SETTING DHXC(10) = 177.03685172 DEG<br>
SETTING I = 10.00000000<br>
SETTING RCCLT = 1.41647919 ANG<br>
SETTING RCHLT = 1.09354699 ANG<br>
SETTING RCOLT = 1.25645120 ANG<br>
SETTING THCOLT = 117.46884981 DEG<br>
SETTING THCCLT = 117.69659884 DEG<br>
SETTING TOCCLT = 123.54974039 DEG<br>
SETTING RCCRT = 1.37817260 ANG<br>
SETTING RCHRT = 1.08768286 ANG<br>
SETTING RCORT = 1.29694643 ANG<br>
SETTING THCORT = 116.50554598 DEG<br>
SETTING THCCRT = 124.25894610 DEG<br>
SETTING TOCCRT = 120.74439212 DEG<br>
SETTING RCHT = 1.07792782 ANG<br>
SETTING RHXTUNT = 0.75768516 ANG<br>
SETTING DHXCT = 177.03685172 DEG<br>
<br>
<br>
Recomputing integrals since basis changed<br>
<br>
<br>
Using spherical harmonics<br>
<br>
Library entry C S cc-pVTZ-F12 selected for orbital group 1<br>
Library entry C P cc-pVTZ-F12 selected for orbital group 1<br>
Library entry C D cc-pVTZ-F12 selected for orbital group 1<br>
Library entry C F cc-pVTZ-F12 selected for orbital group 1<br>
Library entry H S cc-pVTZ-F12 selected for orbital group 2<br>
Library entry H P cc-pVTZ-F12 selected for orbital group 2<br>
Library entry H D cc-pVTZ-F12 selected for orbital group 2<br>
Library entry O S cc-pVTZ-F12 selected for orbital group 4<br>
Library entry O P cc-pVTZ-F12 selected for orbital group 4<br>
Library entry O D cc-pVTZ-F12 selected for orbital group 4<br>
Library entry O F cc-pVTZ-F12 selected for orbital group 4<br>
<br>
1PROGRAM * SEWARD (Integral evaluation for generally contracted gaussian basis sets) Author: Roland Lindh, 1990<br>
<br>
Geometry written to block 41 of record 700<br>
<br>
Orientation using atomic masses<br>
Molecule type: Asymmetric top<br>
Symmetry elements: Y<br>
Rotational constants: 5.5020530 3.5322560 9.8663261 GHz<br>
<br>
Point group Cs<br>
<br>
<br>
<br>
ATOMIC COORDINATES<br>
<br>
NR ATOM CHARGE X Y Z<br>
<br>
1 C1 6.00 2.168357216 0.000000000 -0.075826690<br>
2 H2 1.00 4.203364336 0.000000000 0.013993353<br>
3 C3 6.00 0.820931849 0.000000000 2.237066926<br>
4 C4 6.00 0.798609662 0.000000000 -2.290895167<br>
5 O5 8.00 -1.549417705 0.000000000 2.374815160<br>
6 O6 8.00 -1.651900036 0.000000000 -2.248664521<br>
7 H7 1.00 1.878907728 0.000000000 4.012208372<br>
8 H8 1.00 1.684082566 0.000000000 -4.145807647<br>
9 H12 1.00 -2.088723576 0.000000000 -0.337811341<br>
<br>
Bond lengths in Bohr (Angstrom)<br>
<br>
1-2 2.036988369 1-3 2.676757740 1-4 2.604368776 3-5 2.374348665 3-7 2.066504323<br>
(1.077927820) (1.416479190) (1.378172600) (1.256451200) (1.093546990)<br>
<br>
4-6 2.450873560 4-8 2.055422723 6-9 1.960146597<br>
(1.296946430) (1.087682860) (1.037264905)<br>
<br>
Bond angles<br>
<br>
1-3-5 123.54974039 1-3-7 118.98140980 1-4-6 120.74439212 1-4-8 122.75006190<br>
<br>
2-1-3 117.69659884 2-1-4 124.25894610 3-1-4 118.04445506 4-6-9 103.86393275<br>
<br>
5-3-7 117.46884981 6-4-8 116.50554598<br>
<br>
NUCLEAR CHARGE: 38<br>
NUMBER OF PRIMITIVE AOS: 444<br>
NUMBER OF SYMMETRY AOS: 395<br>
NUMBER OF CONTRACTIONS: 337 ( 227A' + 110A" )<br>
NUMBER OF CORE ORBITALS: 5 ( 5A' + 0A" )<br>
NUMBER OF VALENCE ORBITALS: 24 ( 19A' + 5A" )<br>
<br>
<br>
NUCLEAR REPULSION ENERGY 167.30092648<br>
<br>
<br>
Eigenvalues of metric<br>
<br>
1 0.419E-04 0.758E-04 0.840E-04 0.909E-04 0.102E-03 0.110E-03 0.131E-03 0.135E-03<br>
2 0.461E-03 0.506E-03 0.544E-03 0.686E-03 0.688E-03 0.175E-02 0.234E-02 0.552E-02<br>
<br>
<br>
Contracted 2-electron integrals neglected if value below 1.0D-12<br>
AO integral compression algorithm 1 Integral accuracy 1.0D-12<br>
<br>
3828.089 MB (compressed) written to integral file ( 52.1%)<br>
<br>
Node minimum: 434.110 MB, node maximum: 524.026 MB<br>
<br>
<br>
NUMBER OF SORTED TWO-ELECTRON INTEGRALS: 103414430. BUFFER LENGTH: 32768<br>
NUMBER OF SEGMENTS: 4 SEGMENT LENGTH: 31997334 RECORD LENGTH: 524288<br>
<br>
Memory used in sort: 32.55 MW<br>
<br>
SORT1 READ 919382790. AND WROTE 98240967. INTEGRALS IN 283 RECORDS. CPU TIME: 73.45 SEC, REAL TIME: 84.49 SEC<br>
SORT2 READ 786705259. AND WROTE 827442516. INTEGRALS IN 13752 RECORDS. CPU TIME: 21.04 SEC, REAL TIME: 56.83 SEC<br>
<br>
Node minimum: 103410433. Node maximum: 103496063. integrals<br>
<br>
OPERATOR DM FOR CENTER 0 COORDINATES: 0.000000 0.000000 0.000000<br>
<br>
<br>
**********************************************************************************************************************************<br>
DATASETS * FILE NREC LENGTH (MB) RECORD NAMES<br>
1 18 4.71 500 700 610 900 950 970 1000 129 960 1100<br>
VAR GEOM BASINP SYMINP ZMAT AOBASIS BASIS P2S ABASIS S<br>
1400 1410 1200 1210 1080 1600 1650 1700<br>
T V H0 H01 AOSYM SMH MOLCAS OPER<br>
<br>
2 10 3.93 500 610 700 1000 2100 7360 7350 5300 5600 5500<br>
VAR BASINP GEOM BASIS RHF F12ABS EF12 HESS GINFO HINFO<br>
<br>
PROGRAMS * TOTAL INT RESTART<br>
CPU TIMES * 113.85 113.63 0.08<br>
REAL TIME * 171.95 SEC<br>
DISK USED * 13.27 GB<br>
GA USED * 0.00 MB (max) 0.00 MB (current)<br>
**********************************************************************************************************************************<br>
<br>
1PROGRAM * RHF-SCF (CLOSED SHELL) Authors: W. Meyer, H.-J. Werner<br>
<br>
<br>
NUMBER OF ELECTRONS: 19+ 19-<br>
CONVERGENCE THRESHOLDS: 1.00E-06 (Density) 1.00E-07 (Energy)<br>
MAX. NUMBER OF ITERATIONS: 60<br>
INTERPOLATION TYPE: DIIS<br>
INTERPOLATION STEPS: 2 (START) 1 (STEP)<br>
LEVEL SHIFTS: 0.00 (CLOSED) 0.00 (OPEN)<br>
<br>
<br>
<br>
Molecular orbitals read from record 2100.2 Type=RHF/CANONICAL (state 1.1)<br>
<br>
Molecular orbital dump at record 2100.2<br>
<br>
ITERATION DDIFF GRAD ENERGY 2-EL.EN. DIPOLE MOMENTS DIIS<br>
1 0.000D+00 0.000D+00 -265.73256008 522.551542 1.280075 0.000000 -0.183275 0<br>
2 0.000D+00 0.117D-02 -265.74362170 522.497473 1.272549 0.000000 -0.205194 1<br>
3 0.247D-02 0.288D-03 -265.74499249 522.536734 1.285853 0.000000 -0.218236 2<br>
4 0.502D-03 0.107D-03 -265.74523988 522.501201 1.275387 0.000000 -0.249712 3<br>
5 0.189D-03 0.665D-04 -265.74533370 522.518128 1.280294 0.000000 -0.259703 4<br>
6 0.121D-03 0.270D-04 -265.74536848 522.511558 1.279259 0.000000 -0.275102 5<br>
7 0.998D-04 0.100D-04 -265.74537240 522.513998 1.279604 0.000000 -0.281295 6<br>
8 0.419D-04 0.254D-05 -265.74537254 522.513821 1.279858 0.000000 -0.281927 7<br>
9 0.755D-05 0.102D-05 -265.74537255 522.513443 1.279765 0.000000 -0.282188 8<br>
10 0.260D-05 0.323D-06 -265.74537255 522.513770 1.279852 0.000000 -0.282218 9<br>
11 0.638D-06 0.216D-06 -265.74537255 522.513674 1.279830 0.000000 -0.282211 0<br>
<br>
Final occupancy: 16 3<br>
<br>
!RHF STATE 1.1 Energy -265.745372550858<br>
Nuclear energy 167.30092648<br>
One-electron energy -694.30313587<br>
Two-electron energy 261.25683684<br>
Virial quotient -1.00108852<br>
!RHF STATE 1.1 Dipole moment 1.27982953 0.00000000 -0.28221143<br>
Dipole moment /Debye 3.25278913 0.00000000 -0.71726293<br>
<br>
Orbital energies:<br>
<br>
1.1 2.1 3.1 4.1 5.1 6.1 7.1 8.1 9.1 10.1<br>
-20.597629 -20.558284 -11.353773 -11.344793 -11.229189 -1.445689 -1.380910 -1.068325 -0.903281 -0.834999<br>
<br>
11.1 12.1 13.1 14.1 15.1 16.1 17.1 18.1<br>
-0.723005 -0.671719 -0.648016 -0.626524 -0.559653 -0.438054 0.041669 0.047640<br>
<br>
1.2 2.2 3.2 4.2 5.2<br>
-0.607984 -0.538466 -0.364438 0.059273 0.072260<br>
<br>
<br>
**********************************************************************************************************************************<br>
DATASETS * FILE NREC LENGTH (MB) RECORD NAMES<br>
1 18 4.76 500 700 610 900 950 970 1000 129 960 1100<br>
VAR GEOM BASINP SYMINP ZMAT AOBASIS BASIS P2S ABASIS S<br>
1400 1410 1200 1210 1080 1600 1650 1700<br>
T V H0 H01 AOSYM SMH MOLCAS OPER<br>
<br>
2 10 3.93 500 610 700 1000 2100 7360 7350 5300 5600 5500<br>
VAR BASINP GEOM BASIS RHF F12ABS EF12 HESS GINFO HINFO<br>
<br>
PROGRAMS * TOTAL RHF INT RESTART<br>
CPU TIMES * 140.12 26.26 113.63 0.08<br>
REAL TIME * 202.70 SEC<br>
DISK USED * 13.39 GB<br>
SF USED * 0.26 MB<br>
GA USED * 0.00 MB (max) 0.00 MB (current)<br>
**********************************************************************************************************************************<br>
<br>
<br>
1PROGRAM * CCSD (Closed-shell coupled cluster) Authors: C. Hampel, H.-J. Werner, 1991, M. Deegan, P.J. Knowles, 1992<br>
<br>
CCSD-F12 implementation by H.-J. Werner, 2007<br>
<br>
Density fitting integral evaluation by F. R. Manby, 2003,2007<br>
<br>
Basis set CC-PVTZ/JKFIT generated. Number of basis functions: 515<br>
Basis set CC-PVTZ-F12/OPTRI generated. Number of basis functions: 539<br>
Basis set AUG-CC-PVTZ/MP2FIT generated. Number of basis functions: 714<br>
<br>
Convergence thresholds: THRVAR = 1.00D-12 THRDEN = 1.00D-08<br>
<br>
CCSD(T) terms to be evaluated (factor= 1.000)<br>
<br>
<br>
Number of core orbitals: 5 ( 5 0 )<br>
Number of closed-shell orbitals: 14 ( 11 3 )<br>
Number of external orbitals: 318 ( 211 107 )<br>
<br>
Molecular orbitals read from record 2100.2 Type=RHF/CANONICAL (state 1.1)<br>
<br>
MP2-F12 correlation treatment (H.-J. Werner, 2006)<br>
--------------------------------------------------<br>
<br>
Using MP2-F12 with ansatz 3C(FIX)<br>
<br>
Using projected zeroth-order Hamiltonian (+Z)<br>
<br>
FOCKRIB=T FOCKRIC=T FOCKRIP=T CABSP=T CABSA=T CABSK=T CABSF=T GBC=F EBC=F DMAT=T NOFIK=T NOPAO=1 SOLVE=0 USEPAO=0<br>
EXCH_A= T EXCH_B= F EXCH_C= F EXCH_P= F<br>
<br>
Geminal basis: OPTFULL GEM_TYPE=SLATER BETA=1.0 NGEM=6<br>
Geminal optimization for beta= 1.0000<br>
Weight function: m=0, omega= 1.4646<br>
<br>
Augmented Hessian optimization of geminal fit. Trust ratio= 0.40000<br>
Convergence reached after 2 iterations. Final gradient= 8.00D-16, Step= 4.26D-06, Delta= 1.28D-09<br>
<br>
Alpha: 0.19532 0.81920 2.85917 9.50073 35.69989 197.79328<br>
Coeff: 0.27070 0.30552 0.18297 0.10986 0.06810 0.04224<br>
<br>
<br>
Symmetry turned off. Reference energy: -265.745372550870<br>
<br>
All pairs explicitly correlated. Number of r12-pairs: 105<br>
<br>
AO(A)-basis ORBITAL loaded. Number of functions: 337<br>
RI(R)-basis CC-PVTZ-F12/OPTRI loaded. Number of functions: 539<br>
DF-basis CC-PVTZ/JKFIT loaded. Number of functions: 515<br>
<br>
Screening thresholds: THRAO= 1.00D-10 THRMO= 1.00D-10 THRPROD= 1.00D-10<br>
THRSW= 1.00D-07 THROV= 1.00D-12 THRF12= 1.00D-08<br>
<br>
CPU time for Fock operators 28.23 sec<br>
<br>
Using frozen CA basis dimensions: nDel= 0 nDelAux= -1<br>
<br>
Construction of ABS:<br>
Smallest eigenvalue of S 7.11E-05 (threshold= 1.00E-08)<br>
Ratio eigmin/eigmax 3.26E-07 (threshold= 1.00E-09)<br>
Smallest eigenvalue of S kept 7.11E-05 (threshold= 7.11E-05, 0 functions deleted, 539 kept)<br>
<br>
Construction of CABS:<br>
Smallest eigenvalue of S 4.49E-07 (threshold= 1.00E-08)<br>
Ratio eigmin/eigmax 4.49E-07 (threshold= 1.00E-09)<br>
Smallest eigenvalue of S kept 4.49E-07 (threshold= 4.49E-07, 0 functions deleted, 539 kept)<br>
<br>
CPU time for CABS singles 21.68 sec<br>
<br>
CABS-singles contribution of -0.00230469 patched into reference energy.<br>
New reference energy -265.74767724<br>
<br>
AO(A)-basis ORBITAL loaded. Number of functions: 337<br>
RI(R)-basis CC-PVTZ-F12/OPTRI loaded. Number of functions: 539<br>
DF-basis AUG-CC-PVTZ/MP2FIT loaded. Number of functions: 714<br>
<br>
Screening thresholds: THRAO= 1.00D-10 THRMO= 1.00D-10 THRPROD= 1.00D-10<br>
THRSW= 1.00D-07 THROV= 1.00D-12 THRF12= 1.00D-08<br>
<br>
CPU time for 3-index integral evaluation 26.44 sec<br>
CPU time for first half transformation 1.42 sec ( 2683.2 MFLOP/sec)<br>
CPU time for second half transformation 0.35 sec ( 312.0 MFLOP/sec)<br>
CPU time for sorting 3.36 sec<br>
CPU time for fitting 0.81 sec ( 2348.8 MFLOP/sec)<br>
CPU time for tilde quantities 2.13 sec ( 1197.7 MFLOP/sec)<br>
CPU time for assembly 12.99 sec ( 2699.7 MFLOP/sec)<br>
CPU time for tranop_f12 12.94 sec ( 3555.0 MFLOP/sec)<br>
CPU time for f12 integrals (total) 85.41 sec<br>
F12-matrices built in 8 passes.<br>
F12-matrices built in 8 passes.<br>
F12-matrices built in 8 passes.<br>
F12-matrices built in 8 passes.<br>
FC-matrix built in 8 passes.<br>
CPU time for f12 matrices (total) 12.41 sec<br>
<br>
Diagonal F12 approximation with fixed coefficients: TSING= 0.500, TTRIP= 0.250 (scaled by -1/beta)<br>
<br>
DF-MP2-F12 energy corrections:<br>
------------------------------<br>
Approx. Singlet Triplet Total<br>
DF-MP2-F12/3*C(DX,FIX) -0.079179498845 -0.009329084332 -0.088508583177<br>
DF-MP2-F12/3*C(FIX) -0.074686043380 -0.008703962849 -0.083390006229<br>
DF-MP2-F12/3C(FIX) -0.074619258083 -0.009015376255 -0.083634634338<br>
<br>
DF-MP2-F12 correlation energies:<br>
--------------------------------<br>
Approx. Singlet Triplet Ecorr Total Energy<br>
DF-MP2 -0.615520653077 -0.376272709269 -0.991793362346 -266.739470602214<br>
DF-MP2-F12/3*C(DX,FIX) -0.694700151923 -0.385601793600 -1.080301945523 -266.827979185390<br>
DF-MP2-F12/3*C(FIX) -0.690206696457 -0.384976672118 -1.075183368575 -266.822860608443<br>
DF-MP2-F12/3C(FIX) -0.690139911160 -0.385288085524 -1.075427996684 -266.823105236552<br>
<br>
SCS-DF-MP2 energies (F_SING= 1.20000 F_TRIP= 0.62222 F_PARALLEL= 0.33333):<br>
----------------------------------------------------------------------------<br>
SCS-DF-MP2 -0.972750025015 -266.720427264883<br>
SCS-DF-MP2-F12/3*C(DX,FIX) -1.073570187214 -266.821247427082<br>
SCS-DF-MP2-F12/3*C(FIX) -1.067789076177 -266.815466316045<br>
SCS-DF-MP2-F12/3C(FIX) -1.067902702163 -266.815579942030<br>
<br>
<br>
Symmetry restored.<br>
Symmetry transformation completed.<br>
<br>
Number of N-1 electron functions: 14<br>
Number of N-2 electron functions: 105<br>
Number of singly external CSFs: 2642<br>
Number of doubly external CSFs: 5130358<br>
Total number of CSFs: 5133001<br>
<br>
Pair and operator lists are different<br>
<br>
Length of J-op integral file: 1.29 GB<br>
Length of K-op integral file: 1.45 GB<br>
Length of 3-ext integral record: 0.00 MB<br>
<br>
For full I/O caching in triples, increase memory by155552668 words to 283.8 Mword<br>
NPASS IN TRSOR3 2<br>
<br>
Integral transformation finished. Total CPU: 298.80 sec, npass= 1 Memory used: 82.01 MW<br>
<br>
Reference energy: -265.74767724<br>
<br>
Adding F12 terms to K(Cij), methodcc=6, factor=1.0<br>
<br>
ITER. SQ.NORM CORR.ENERGY TOTAL ENERGY ENERGY CHANGE DEN1 VAR(S) VAR(P) DIIS TIME<br>
1 1.26432057 -0.98651566 -266.73419290 0.00538108 0.03537158 0.84D-02 0.57D-02 1 1 516.83<br>
2 1.28388838 -0.98526539 -266.73294263 0.00125027 -0.02303067 0.20D-03 0.80D-03 2 2 618.24<br>
3 1.29273234 -0.98585371 -266.73353095 -0.00058831 -0.00335794 0.30D-03 0.54D-04 3 3 719.76<br>
4 1.29809943 -0.98579533 -266.73347257 0.00005838 -0.00207146 0.42D-04 0.14D-04 4 4 821.84<br>
5 1.30111023 -0.98580438 -266.73348162 -0.00000905 -0.00063593 0.11D-04 0.18D-05 5 5 923.60<br>
6 1.30233570 -0.98575943 -266.73343667 0.00004495 -0.00008957 0.86D-06 0.40D-06 6 6 1025.71<br>
7 1.30266405 -0.98574494 -266.73342218 0.00001449 -0.00002586 0.18D-06 0.56D-07 6 1 1128.45<br>
8 1.30270775 -0.98574589 -266.73342313 -0.00000094 0.00000315 0.38D-07 0.12D-07 6 3 1231.21<br>
9 1.30275207 -0.98574262 -266.73341986 0.00000326 -0.00001171 0.91D-08 0.20D-08 6 2 1333.53<br>
10 1.30275499 -0.98574293 -266.73342017 -0.00000031 0.00000035 0.29D-08 0.62D-09 6 4 1435.98<br>
11 1.30275653 -0.98574311 -266.73342035 -0.00000018 0.00000019 0.82D-09 0.19D-09 6 5 1538.35<br>
12 1.30275340 -0.98574339 -266.73342063 -0.00000028 0.00000095 0.15D-09 0.53D-10 6 1 1640.82<br>
13 1.30275306 -0.98574348 -266.73342072 -0.00000009 0.00000037 0.16D-10 0.78D-11 6 6 1742.90<br>
14 1.30275288 -0.98574345 -266.73342069 0.00000003 -0.00000010 0.34D-11 0.16D-11 6 2 1845.79<br>
15 1.30275303 -0.98574345 -266.73342069 0.00000000 -0.00000003 0.39D-12 0.29D-12 6 3 1948.55<br>
16 1.30275305 -0.98574344 -266.73342068 0.00000001 -0.00000003 0.73D-13 0.39D-13 6 1 2050.74<br>
<br>
Norm of t1 vector: 0.13806819 S-energy: 0.00000048 T1 diagnostic: 0.01845014<br>
D1 diagnostic: 0.07307945<br>
<br>
Singles amplitudes (print threshold = 0.500E-01):<br>
<br>
I SYM. A A T(IA)<br>
<br>
13 2 1 -0.05421522<br>
<br>
Adding energy correction from Wij-K(Fij), methodcc=6 factor=1.0<br>
<br>
F12b singlet correction 0.013617243080<br>
F12b triplet correction 0.001048026430<br>
F12b total correction 0.014665269510<br>
<br>
CPU time for triples: 2032.89 sec for i = 1 - 13 / 14 (~ 81.0 % done)<br>
CPU time for triples: 510.27 sec for i = 14 - 14 / 14 ( 100.0 % done)<br>
<br>
Total CPU time for triples: 2543.66 sec<br>
<br>
<br>
RESULTS<br>
=======<br>
<br>
Reference energy -265.745372550870<br>
F12 singles correction -0.002304688998<br>
<br>
F12 singles corrections added to reference energy<br>
<br>
New reference energy -265.747677239868<br>
<br>
F12b singlet correction -0.060401462859<br>
F12b triplet correction -0.007630120151<br>
F12b total correction -0.068031583010<br>
<br>
F12b corrections for ansatz F12/3C(FIX) added to CCSD energy<br>
<br>
CCSD-F12b singlet pair energy -0.711797917864<br>
CCSD-F12b triplet pair energy -0.341977585683<br>
<br>
CCSD-F12b correlation energy -1.053775024658<br>
Triples (T) contribution -0.048810207392<br>
Total correlation energy -1.102585232050<br>
<br>
CCSD-F12b total energy -266.801452264526<br>
CCSD[T]-F12b energy -266.852943141286<br>
CCSD-T-F12b energy -266.849416345160<br>
!CCSD(T)-F12b total energy -266.850262471918<br>
<br>
Program statistics:<br>
<br>
Available memory in ccsd: 127998387<br>
Min. memory needed in ccsd: 12893024<br>
Max. memory used in ccsd: 17563958<br>
Max. memory used in cckext: 18713946 (17 integral passes)<br>
<br>
<br>
<br>
**********************************************************************************************************************************<br>
DATASETS * FILE NREC LENGTH (MB) RECORD NAMES<br>
1 20 1002.68 500 700 610 900 950 970 1000 129 960 1100<br>
VAR GEOM BASINP SYMINP ZMAT AOBASIS BASIS P2S ABASIS S<br>
1400 1410 1200 1210 1080 1600 1650 1700 960(1) 1380<br>
T V H0 H01 AOSYM SMH MOLCAS OPER ABASIS JKOP<br>
<br>
2 10 3.93 500 610 700 1000 2100 7360 7350 5300 5600 5500<br>
VAR BASINP GEOM BASIS RHF F12ABS EF12 HESS GINFO HINFO<br>
<br>
PROGRAMS * TOTAL CCSD(T) RHF INT RESTART<br>
CPU TIMES * 4819.11 4678.98 26.26 113.63 0.08<br>
REAL TIME * 5323.02 SEC<br>
DISK USED * 14.92 GB<br>
SF USED * 4.48 GB<br>
GA USED * 2.40 MB (max) 0.00 MB (current)<br>
**********************************************************************************************************************************<br>
<br>
<br>
PROGRAM * OPT (Geometry optimization) Authors: F. Eckert and H.-J. Werner<br>
<br>
<br>
*** Long output written to logfile /users/nesbitt/djn/molpro/malonaldehyde/newgeom_vnz12/vtzf12_irc_full_restart/mal_vtzf12_irc_full_restart.log ***<br>
<br>
<br>
Geometry optimization using default procedure for command CCSD(T)-F12B<br>
<br>
Inactive variables: RXCT<br>
<br>
Active variables: RCHT RCCLT THCCLT RCCRT THCCRT RCOLT<br>
TOCCLT RCORT TOCCRT RCHLT THCOLT RCHRT<br>
THCORT RHXTUNT DHXCT<br>
<br>
Geometry written to block 1 of record 700<br>
<br>
Optimized active variables along the reaction path will be saved in file mal_vtzf12_irc_full_restart.act<br>
<br>
Numerically approximating hessian using central energy differences<br>
<br>
Task list generated. Total number of displacements: 240<br>
<br>
<br>
<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>
<br></blockquote></div><br>