From 94de9393e00638a0921377c807293117056bfe8d Mon Sep 17 00:00:00 2001 From: eginer Date: Tue, 2 Apr 2019 12:28:53 +0200 Subject: [PATCH] added Li2.out BeH.out LiH.out Na2.out for mixed core in vqz --- G09/Mixed_core/Molecules/vqz/BeH.out | 282 +++++++++++------------ G09/Mixed_core/Molecules/vqz/Li2.out | 295 ++++++++++++------------ G09/Mixed_core/Molecules/vqz/LiH.out | 270 +++++++++++----------- G09/Mixed_core/Molecules/vqz/Na2.out | 321 ++++++++++++--------------- G09/Mixed_core/Molecules/vqz/copy.sh | 7 +- Manuscript/G2-srDFT.pdf | Bin 476877 -> 477055 bytes Manuscript/G2-srDFT.tex | 15 +- 7 files changed, 575 insertions(+), 615 deletions(-) diff --git a/G09/Mixed_core/Molecules/vqz/BeH.out b/G09/Mixed_core/Molecules/vqz/BeH.out index 2ea09f3..4e92371 100644 --- a/G09/Mixed_core/Molecules/vqz/BeH.out +++ b/G09/Mixed_core/Molecules/vqz/BeH.out @@ -2,8 +2,8 @@ Input=BeH.inp Output=BeH.out Initial command: - /share/apps/gaussian/g09d01/nehalem/g09/l1.exe "/mnt/beegfs/tmpdir/42176/Gau-41101.inp" -scrdir="/mnt/beegfs/tmpdir/42176/" - Entering Link 1 = /share/apps/gaussian/g09d01/nehalem/g09/l1.exe PID= 41102. + /share/apps/gaussian/g09d01/nehalem/g09/l1.exe "/mnt/beegfs/tmpdir/42179/Gau-51886.inp" -scrdir="/mnt/beegfs/tmpdir/42179/" + Entering Link 1 = /share/apps/gaussian/g09d01/nehalem/g09/l1.exe PID= 51887. Copyright (c) 1988,1990,1992,1993,1995,1998,2003,2009,2013, Gaussian, Inc. All Rights Reserved. @@ -76,21 +76,21 @@ ****************************************** Gaussian 09: ES64L-G09RevD.01 24-Apr-2013 - 29-Mar-2019 + 30-Mar-2019 ****************************************** - ------------------------------------- - #p ROCCSD(T) cc-pVQZ pop=full gfprint - ------------------------------------- + ------------------------------------------------------------- + #p ROCCSD(T,FreezeInnerNobleGasCore) cc-pVQZ pop=full gfprint + ------------------------------------------------------------- 1/38=1/1; 2/12=2,17=6,18=5,40=1/2; 3/5=16,6=2,11=2,16=1,24=100,25=1,30=1,116=101/1,2,3; 4//1; 5/5=2,38=5/2; - 8/5=-1,6=4,9=120000,10=1/1,4; + 8/5=-1,6=4,9=120000,10=3/1,4; 9/5=7,14=2/13; 6/7=3/1; 99/5=1,9=1/99; - Leave Link 1 at Fri Mar 29 11:44:14 2019, MaxMem= 0 cpu: 0.0 + Leave Link 1 at Sat Mar 30 23:17:43 2019, MaxMem= 0 cpu: 0.0 (Enter /share/apps/gaussian/g09d01/nehalem/g09/l101.exe) -- G2 @@ -116,7 +116,7 @@ NQMom= 5.2880000 0.0000000 NMagM= -1.1779000 2.7928460 AtZNuc= 4.0000000 1.0000000 - Leave Link 101 at Fri Mar 29 11:44:14 2019, MaxMem= 33554432 cpu: 0.0 + Leave Link 101 at Sat Mar 30 23:17:44 2019, MaxMem= 33554432 cpu: 0.1 (Enter /share/apps/gaussian/g09d01/nehalem/g09/l202.exe) Input orientation: --------------------------------------------------------------------- @@ -141,7 +141,7 @@ 2 1 0 0.000000 0.000000 -1.075046 --------------------------------------------------------------------- Rotational constants (GHZ): 0.0000000 308.7422835 308.7422835 - Leave Link 202 at Fri Mar 29 11:44:14 2019, MaxMem= 33554432 cpu: 0.0 + Leave Link 202 at Sat Mar 30 23:17:44 2019, MaxMem= 33554432 cpu: 0.0 (Enter /share/apps/gaussian/g09d01/nehalem/g09/l301.exe) Standard basis: CC-pVQZ (5D, 7F) Ernie: Thresh= 0.10000D-02 Tol= 0.10000D-05 Strict=F. @@ -231,17 +231,17 @@ Integral buffers will be 131072 words long. Raffenetti 2 integral format. Two-electron integral symmetry is turned on. - Leave Link 301 at Fri Mar 29 11:44:14 2019, MaxMem= 33554432 cpu: 0.0 + Leave Link 301 at Sat Mar 30 23:17:44 2019, MaxMem= 33554432 cpu: 0.0 (Enter /share/apps/gaussian/g09d01/nehalem/g09/l302.exe) NPDir=0 NMtPBC= 1 NCelOv= 1 NCel= 1 NClECP= 1 NCelD= 1 NCelK= 1 NCelE2= 1 NClLst= 1 CellRange= 0.0. One-electron integrals computed using PRISM. NBasis= 85 RedAO= T EigKep= 1.99D-03 NBF= 35 10 20 20 NBsUse= 85 1.00D-06 EigRej= -1.00D+00 NBFU= 35 10 20 20 - Leave Link 302 at Fri Mar 29 11:44:14 2019, MaxMem= 33554432 cpu: 0.1 + Leave Link 302 at Sat Mar 30 23:17:44 2019, MaxMem= 33554432 cpu: 0.1 (Enter /share/apps/gaussian/g09d01/nehalem/g09/l303.exe) DipDrv: MaxL=1. - Leave Link 303 at Fri Mar 29 11:44:14 2019, MaxMem= 33554432 cpu: 0.0 + Leave Link 303 at Sat Mar 30 23:17:44 2019, MaxMem= 33554432 cpu: 0.0 (Enter /share/apps/gaussian/g09d01/nehalem/g09/l401.exe) ExpMin= 4.26D-02 ExpMax= 1.46D+04 ExpMxC= 1.41D+02 IAcc=2 IRadAn= 4 AccDes= 0.00D+00 Harris functional with IExCor= 205 and IRadAn= 4 diagonalized for initial guess. @@ -267,7 +267,7 @@ (SG) (PI) (PI) (SG) (SG) (PI) (PI) (DLTA) (DLTA) (PI) (PI) (SG) (SG) The electronic state of the initial guess is 2-SG. - Leave Link 401 at Fri Mar 29 11:44:15 2019, MaxMem= 33554432 cpu: 0.1 + Leave Link 401 at Sat Mar 30 23:17:44 2019, MaxMem= 33554432 cpu: 0.2 (Enter /share/apps/gaussian/g09d01/nehalem/g09/l502.exe) Restricted open shell SCF: Using DIIS extrapolation, IDIIS= 1040. @@ -444,7 +444,7 @@ KE= 1.514824875785D+01 PE=-3.793381586314D+01 EE= 6.057350472991D+00 Annihilation of the first spin contaminant: S**2 before annihilation 0.7500, after 0.7500 - Leave Link 502 at Fri Mar 29 11:44:17 2019, MaxMem= 33554432 cpu: 1.6 + Leave Link 502 at Sat Mar 30 23:17:46 2019, MaxMem= 33554432 cpu: 1.4 (Enter /share/apps/gaussian/g09d01/nehalem/g09/l801.exe) Windowed orbitals will be sorted by symmetry type. GenMOA: NOpAll= 4 NOp2=4 NOpUse= 4 JSym2X=1 @@ -456,59 +456,51 @@ NMat0= 1 NMatS0= 1 NMatT0= 0 NMatD0= 1 NMtDS0= 0 NMtDT0= 0 Petite list used in FoFCou. = 0.0000 = 0.0000 = 0.5000 = 0.7500 S= 0.5000 - ExpMin= 4.26D-02 ExpMax= 1.46D+04 ExpMxC= 1.41D+02 IAcc=3 IRadAn= 5 AccDes= 0.00D+00 - HarFok: IExCor= 205 AccDes= 0.00D+00 IRadAn= 5 IDoV=-2 UseB2=F ITyADJ=14 - ICtDFT= 12500011 ScaDFX= 1.000000 1.000000 1.000000 1.000000 - Largest valence mixing into a core orbital is 4.51D-04 - Largest core mixing into a valence orbital is 7.51D-05 - Largest valence mixing into a core orbital is 4.77D-04 - Largest core mixing into a valence orbital is 6.53D-05 - Range of M.O.s used for correlation: 2 85 - NBasis= 85 NAE= 3 NBE= 2 NFC= 1 NFV= 0 - NROrb= 84 NOA= 2 NOB= 1 NVA= 82 NVB= 83 - Singles contribution to E2= -0.1449771289D-03 - Leave Link 801 at Fri Mar 29 11:44:18 2019, MaxMem= 33554432 cpu: 1.0 + Range of M.O.s used for correlation: 1 85 + NBasis= 85 NAE= 3 NBE= 2 NFC= 0 NFV= 0 + NROrb= 85 NOA= 3 NOB= 2 NVA= 82 NVB= 83 + Singles contribution to E2= -0.1556399356D-03 + Leave Link 801 at Sat Mar 30 23:17:47 2019, MaxMem= 33554432 cpu: 0.9 (Enter /share/apps/gaussian/g09d01/nehalem/g09/l804.exe) Open-shell transformation, MDV= 33554432 ITran=4 ISComp=2. Semi-Direct transformation. - ModeAB= 2 MOrb= 2 LenV= 33119195 - LASXX= 149226 LTotXX= 149226 LenRXX= 149226 - LTotAB= 156890 MaxLAS= 934920 LenRXY= 934920 - NonZer= 991368 LenScr= 2162688 LnRSAI= 0 - LnScr1= 0 LExtra= 0 Total= 3246834 + ModeAB= 2 MOrb= 3 LenV= 33116444 + LASXX= 231189 LTotXX= 231189 LenRXX= 231189 + LTotAB= 243210 MaxLAS= 1419075 LenRXY= 1419075 + NonZer= 1504755 LenScr= 2883584 LnRSAI= 0 + LnScr1= 0 LExtra= 0 Total= 4533848 MaxDsk= -1 SrtSym= F ITran= 4 DoSDTr: NPSUse= 1 - JobTyp=1 Pass 1: I= 1 to 2. + JobTyp=1 Pass 1: I= 1 to 3. (rs|ai) integrals will be sorted in core. Complete sort for first half transformation. First half transformation complete. Complete sort for second half transformation. Second half transformation complete. - ModeAB= 2 MOrb= 1 LenV= 33119195 - LASXX= 75138 LTotXX= 75138 LenRXX= 467460 - LTotAB= 74902 MaxLAS= 467460 LenRXY= 74902 - NonZer= 495684 LenScr= 1441792 LnRSAI= 0 - LnScr1= 0 LExtra= 0 Total= 1984154 + ModeAB= 2 MOrb= 2 LenV= 33116444 + LASXX= 155274 LTotXX= 155274 LenRXX= 946050 + LTotAB= 150807 MaxLAS= 946050 LenRXY= 150807 + NonZer= 1003170 LenScr= 2162688 LnRSAI= 0 + LnScr1= 0 LExtra= 0 Total= 3259545 MaxDsk= -1 SrtSym= F ITran= 4 DoSDTr: NPSUse= 1 - JobTyp=2 Pass 1: I= 1 to 1. + JobTyp=2 Pass 1: I= 1 to 2. (rs|ai) integrals will be sorted in core. Complete sort for first half transformation. First half transformation complete. Complete sort for second half transformation. Second half transformation complete. Spin components of T(2) and E(2): - alpha-alpha T2 = 0.9270092877D-03 E2= -0.1600379917D-02 - alpha-beta T2 = 0.1443890931D-01 E2= -0.3341011378D-01 - beta-beta T2 = 0.0000000000D+00 E2= 0.0000000000D+00 - ANorm= 0.1007737814D+01 - E2 = -0.3515547082D-01 EUMP2 = -0.15188214126894D+02 + alpha-alpha T2 = 0.1033220160D-02 E2= -0.2505054873D-02 + alpha-beta T2 = 0.1563870010D-01 E2= -0.5293459499D-01 + beta-beta T2 = 0.4383376928D-04 E2= -0.3798079455D-03 + ANorm= 0.1008408049D+01 + E2 = -0.5597509774D-01 EUMP2 = -0.15209033753818D+02 (S**2,0)= 0.75000D+00 (S**2,1)= 0.75000D+00 - E(PUHF)= -0.15153058656D+02 E(PMP2)= -0.15188214127D+02 - Leave Link 804 at Fri Mar 29 11:44:22 2019, MaxMem= 33554432 cpu: 3.8 + E(PUHF)= -0.15153058656D+02 E(PMP2)= -0.15209033754D+02 + Leave Link 804 at Sat Mar 30 23:17:51 2019, MaxMem= 33554432 cpu: 4.1 (Enter /share/apps/gaussian/g09d01/nehalem/g09/l913.exe) CIDS: MDV= 33554432. - Frozen-core window: NFC= 1 NFV= 0. IFCWin=0 IBDFC=1 NFBD= 0 0 NFCmp= 0 0 NFFFC= 0 0 Using original routines for 1st iteration, S=T. Using DD4UQ or CC4UQ for 2nd and later iterations. @@ -523,107 +515,115 @@ Iterations= 50 Convergence= 0.100D-06 Iteration Nr. 1 ********************** - DD1Dir will call FoFMem 1 times, MxPair= 6 - NAB= 2 NAA= 1 NBB= 0. - DD1Dir will call FoFMem 1 times, MxPair= 6 - NAB= 2 NAA= 1 NBB= 0. - MP4(R+Q)= 0.68639559D-02 + DD1Dir will call FoFMem 1 times, MxPair= 16 + NAB= 6 NAA= 3 NBB= 1. + DD1Dir will call FoFMem 1 times, MxPair= 16 + NAB= 6 NAA= 3 NBB= 1. + MP4(R+Q)= 0.88987974D-02 Maximum subspace dimension= 5 - Norm of the A-vectors is 4.2947926D-03 conv= 1.00D-05. - RLE energy= -0.0346245806 - E3= -0.63224568D-02 EROMP3= -0.15194536584D+02 - E4(SDQ)= -0.18578919D-02 ROMP4(SDQ)= -0.15196394476D+02 + Norm of the A-vectors is 4.8629723D-03 conv= 1.00D-05. + RLE energy= -0.0553690138 + E3= -0.82838103D-02 EROMP3= -0.15217317564D+02 + E4(SDQ)= -0.20179625D-02 ROMP4(SDQ)= -0.15219335527D+02 VARIATIONAL ENERGIES WITH THE FIRST-ORDER WAVEFUNCTION: - DE(Corr)= -0.34616265E-01 E(Corr)= -15.187674921 - NORM(A)= 0.10074396D+01 + DE(Corr)= -0.55362253E-01 E(Corr)= -15.208420909 + NORM(A)= 0.10080847D+01 Iteration Nr. 2 ********************** - DD1Dir will call FoFMem 1 times, MxPair= 6 - NAB= 2 NAA= 1 NBB= 0. - Norm of the A-vectors is 5.7205604D-02 conv= 1.00D-05. - RLE energy= -0.0350072695 - DE(Corr)= -0.40828656E-01 E(CORR)= -15.193887312 Delta=-6.21D-03 - NORM(A)= 0.10076665D+01 + DD1Dir will call FoFMem 1 times, MxPair= 16 + NAB= 6 NAA= 3 NBB= 1. + Norm of the A-vectors is 6.5031618D-02 conv= 1.00D-05. + RLE energy= -0.0561633433 + DE(Corr)= -0.63525452E-01 E(CORR)= -15.216584108 Delta=-8.16D-03 + NORM(A)= 0.10084509D+01 Iteration Nr. 3 ********************** - DD1Dir will call FoFMem 1 times, MxPair= 6 - NAB= 2 NAA= 1 NBB= 0. - Norm of the A-vectors is 5.5124680D-02 conv= 1.00D-05. - RLE energy= -0.0373513727 - DE(Corr)= -0.40960682E-01 E(CORR)= -15.194019339 Delta=-1.32D-04 - NORM(A)= 0.10092446D+01 + DD1Dir will call FoFMem 1 times, MxPair= 16 + NAB= 6 NAA= 3 NBB= 1. + Norm of the A-vectors is 6.1048855D-02 conv= 1.00D-05. + RLE energy= -0.0587906226 + DE(Corr)= -0.63753835E-01 E(CORR)= -15.216812491 Delta=-2.28D-04 + NORM(A)= 0.10099211D+01 Iteration Nr. 4 ********************** - DD1Dir will call FoFMem 1 times, MxPair= 6 - NAB= 2 NAA= 1 NBB= 0. - Norm of the A-vectors is 4.2102983D-02 conv= 1.00D-05. - RLE energy= -0.0442574873 - DE(Corr)= -0.41792032E-01 E(CORR)= -15.194850688 Delta=-8.31D-04 - NORM(A)= 0.10161980D+01 + DD1Dir will call FoFMem 1 times, MxPair= 16 + NAB= 6 NAA= 3 NBB= 1. + Norm of the A-vectors is 4.7250585D-02 conv= 1.00D-05. + RLE energy= -0.0192579678 + DE(Corr)= -0.64566077E-01 E(CORR)= -15.217624733 Delta=-8.12D-04 + NORM(A)= 0.10150651D+01 Iteration Nr. 5 ********************** - DD1Dir will call FoFMem 1 times, MxPair= 6 - NAB= 2 NAA= 1 NBB= 0. - Norm of the A-vectors is 2.8345381D-03 conv= 1.00D-05. - RLE energy= -0.0452151015 - DE(Corr)= -0.44381988E-01 E(CORR)= -15.197440644 Delta=-2.59D-03 - NORM(A)= 0.10177326D+01 + DD1Dir will call FoFMem 1 times, MxPair= 16 + NAB= 6 NAA= 3 NBB= 1. + Norm of the A-vectors is 2.5253002D-01 conv= 1.00D-05. + RLE energy= -0.0633623006 + DE(Corr)= -0.51173133E-01 E(CORR)= -15.204231789 Delta= 1.34D-02 + NORM(A)= 0.10139960D+01 Iteration Nr. 6 ********************** - DD1Dir will call FoFMem 1 times, MxPair= 6 - NAB= 2 NAA= 1 NBB= 0. - Norm of the A-vectors is 4.3015053D-03 conv= 1.00D-05. - RLE energy= -0.0444954946 - DE(Corr)= -0.44759183E-01 E(CORR)= -15.197817839 Delta=-3.77D-04 - NORM(A)= 0.10168545D+01 + DD1Dir will call FoFMem 1 times, MxPair= 16 + NAB= 6 NAA= 3 NBB= 1. + Norm of the A-vectors is 2.1365947D-02 conv= 1.00D-05. + RLE energy= -0.0675217238 + DE(Corr)= -0.66165651E-01 E(CORR)= -15.219224307 Delta=-1.50D-02 + NORM(A)= 0.10176796D+01 Iteration Nr. 7 ********************** - DD1Dir will call FoFMem 1 times, MxPair= 6 - NAB= 2 NAA= 1 NBB= 0. - Norm of the A-vectors is 2.7188685D-04 conv= 1.00D-05. - RLE energy= -0.0445094685 - DE(Corr)= -0.44501145E-01 E(CORR)= -15.197559801 Delta= 2.58D-04 - NORM(A)= 0.10168842D+01 + DD1Dir will call FoFMem 1 times, MxPair= 16 + NAB= 6 NAA= 3 NBB= 1. + Norm of the A-vectors is 1.1848706D-03 conv= 1.00D-05. + RLE energy= -0.0673402142 + DE(Corr)= -0.67447047E-01 E(CORR)= -15.220505703 Delta=-1.28D-03 + NORM(A)= 0.10173206D+01 Iteration Nr. 8 ********************** - DD1Dir will call FoFMem 1 times, MxPair= 6 - NAB= 2 NAA= 1 NBB= 0. - Norm of the A-vectors is 1.1345976D-04 conv= 1.00D-05. - RLE energy= -0.0445066949 - DE(Corr)= -0.44507049E-01 E(CORR)= -15.197565705 Delta=-5.90D-06 - NORM(A)= 0.10168792D+01 + DD1Dir will call FoFMem 1 times, MxPair= 16 + NAB= 6 NAA= 3 NBB= 1. + Norm of the A-vectors is 4.1237335D-04 conv= 1.00D-05. + RLE energy= -0.0674029875 + DE(Corr)= -0.67369733E-01 E(CORR)= -15.220428389 Delta= 7.73D-05 + NORM(A)= 0.10174465D+01 Iteration Nr. 9 ********************** - DD1Dir will call FoFMem 1 times, MxPair= 6 - NAB= 2 NAA= 1 NBB= 0. - Norm of the A-vectors is 3.5813154D-05 conv= 1.00D-05. - RLE energy= -0.0445062365 - DE(Corr)= -0.44506240E-01 E(CORR)= -15.197564896 Delta= 8.10D-07 - NORM(A)= 0.10168789D+01 + DD1Dir will call FoFMem 1 times, MxPair= 16 + NAB= 6 NAA= 3 NBB= 1. + Norm of the A-vectors is 8.8143018D-05 conv= 1.00D-05. + RLE energy= -0.0673955792 + DE(Corr)= -0.67397746E-01 E(CORR)= -15.220456402 Delta=-2.80D-05 + NORM(A)= 0.10174314D+01 Iteration Nr. 10 ********************** - DD1Dir will call FoFMem 1 times, MxPair= 6 - NAB= 2 NAA= 1 NBB= 0. - Norm of the A-vectors is 1.2327548D-05 conv= 1.00D-05. - RLE energy= -0.0445060648 - DE(Corr)= -0.44506099E-01 E(CORR)= -15.197564755 Delta= 1.41D-07 - NORM(A)= 0.10168787D+01 + DD1Dir will call FoFMem 1 times, MxPair= 16 + NAB= 6 NAA= 3 NBB= 1. + Norm of the A-vectors is 2.9502179D-05 conv= 1.00D-05. + RLE energy= -0.0673936367 + DE(Corr)= -0.67394539E-01 E(CORR)= -15.220453195 Delta= 3.21D-06 + NORM(A)= 0.10174269D+01 Iteration Nr. 11 ********************** - DD1Dir will call FoFMem 1 times, MxPair= 6 - NAB= 2 NAA= 1 NBB= 0. - Norm of the A-vectors is 3.4979351D-06 conv= 1.00D-05. - RLE energy= -0.0445061360 - DE(Corr)= -0.44506067E-01 E(CORR)= -15.197564723 Delta= 3.14D-08 - NORM(A)= 0.10168791D+01 - CI/CC converged in 11 iterations to DelEn= 3.14D-08 Conv= 1.00D-07 ErrA1= 3.50D-06 Conv= 1.00D-05 - Largest amplitude= 3.26D-02 - Time for triples= 5.21 seconds. - T4(CCSD)= -0.57412290D-03 - T5(CCSD)= -0.17865262D-05 - CCSD(T)= -0.15198140633D+02 + DD1Dir will call FoFMem 1 times, MxPair= 16 + NAB= 6 NAA= 3 NBB= 1. + Norm of the A-vectors is 6.7788554D-06 conv= 1.00D-05. + RLE energy= -0.0673936553 + DE(Corr)= -0.67393643E-01 E(CORR)= -15.220452299 Delta= 8.96D-07 + NORM(A)= 0.10174269D+01 + Iteration Nr. 12 + ********************** + DD1Dir will call FoFMem 1 times, MxPair= 16 + NAB= 6 NAA= 3 NBB= 1. + Norm of the A-vectors is 2.4159716D-06 conv= 1.00D-05. + RLE energy= -0.0673936473 + DE(Corr)= -0.67393638E-01 E(CORR)= -15.220452294 Delta= 5.49D-09 + NORM(A)= 0.10174270D+01 + CI/CC converged in 12 iterations to DelEn= 5.49D-09 Conv= 1.00D-07 ErrA1= 2.42D-06 Conv= 1.00D-05 + Largest amplitude= 3.20D-02 + Time for triples= 11.24 seconds. + T4(CCSD)= -0.72605184D-03 + T5(CCSD)= 0.42959536D-05 + CCSD(T)= -0.15221174050D+02 Discarding MO integrals. - Leave Link 913 at Fri Mar 29 11:45:02 2019, MaxMem= 33554432 cpu: 8.3 + Leave Link 913 at Sat Mar 30 23:19:07 2019, MaxMem= 33554432 cpu: 15.7 (Enter /share/apps/gaussian/g09d01/nehalem/g09/l601.exe) Copying SCF densities to generalized density rwf, IOpCl= 0 IROHF=1. @@ -4761,19 +4761,27 @@ --------------------------------------------------------------------------------- No NMR shielding tensors so no spin-rotation constants. - Leave Link 601 at Fri Mar 29 11:45:02 2019, MaxMem= 33554432 cpu: 0.3 + Leave Link 601 at Sat Mar 30 23:19:08 2019, MaxMem= 33554432 cpu: 0.2 (Enter /share/apps/gaussian/g09d01/nehalem/g09/l9999.exe) - 1\1\GINC-COMPUTE-7-0\SP\ROCCSD(T)-FC\CC-pVQZ\Be1H1(2)\LOOS\29-Mar-2019 - \0\\#p ROCCSD(T) cc-pVQZ pop=full gfprint\\G2\\0,2\Be\H,1,1.34380733\\ - Version=ES64L-G09RevD.01\State=2-SG\HF=-15.1530587\MP2=-15.1882141\MP3 - =-15.1945366\PUHF=-15.1530587\PMP2-0=-15.1882141\MP4SDQ=-15.1963945\CC - SD=-15.1975647\CCSD(T)=-15.1981406\RMSD=2.914e-09\PG=C*V [C*(H1Be1)]\\ - @ + 1\1\GINC-COMPUTE-7-0\SP\ROCCSD(T)-FC1\CC-pVQZ\Be1H1(2)\LOOS\30-Mar-201 + 9\0\\#p ROCCSD(T,FreezeInnerNobleGasCore) cc-pVQZ pop=full gfprint\\G2 + \\0,2\Be\H,1,1.34380733\\Version=ES64L-G09RevD.01\State=2-SG\HF=-15.15 + 30587\MP2=-15.2090338\MP3=-15.2173176\PUHF=-15.1530587\PMP2-0=-15.2090 + 338\MP4SDQ=-15.2193355\CCSD=-15.2204523\CCSD(T)=-15.221174\RMSD=2.914e + -09\PG=C*V [C*(H1Be1)]\\@ - THE RED LIGHT IS ALWAYS LONGER THAN THE GREEN LIGHT. + THESE ARE THE TIMES THAT TRY MENS SOULS. + THE SUMMER SOLDIER AND THE SUNSHINE PATRIOT WILL IN THIS CRISIS, + SHRINK FROM THE SERVICE OF HIS COUNTRY. BUT HE THAT STANDS NOW, + DESERVES THE LOVE AND THANKS OF MAN AND WOMAN. + TYRANNY, LIKE HELL, IS NOT EASILY CONQUERED, + YET WE HAVE THIS CONSOLATION WITH US, + THAT THE HARDER THE CONFLICT, THE MORE GLORIOUS THE TRIUMPH. + WHAT WE OBTAIN TOO CHEAP, WE ESTEEM TOO LIGHTLY, + 'TIS DEARNESS ONLY THAT GIVES EVERYTHING ITS VALUE. - -- PETER'S THEORY OF RELATIVITY - Job cpu time: 0 days 0 hours 0 minutes 15.9 seconds. - File lengths (MBytes): RWF= 100 Int= 0 D2E= 0 Chk= 2 Scr= 1 - Normal termination of Gaussian 09 at Fri Mar 29 11:45:02 2019. + -- TOM PAINE + Job cpu time: 0 days 0 hours 0 minutes 23.1 seconds. + File lengths (MBytes): RWF= 112 Int= 0 D2E= 0 Chk= 2 Scr= 1 + Normal termination of Gaussian 09 at Sat Mar 30 23:19:08 2019. diff --git a/G09/Mixed_core/Molecules/vqz/Li2.out b/G09/Mixed_core/Molecules/vqz/Li2.out index 20b440c..fa65611 100644 --- a/G09/Mixed_core/Molecules/vqz/Li2.out +++ b/G09/Mixed_core/Molecules/vqz/Li2.out @@ -2,8 +2,8 @@ Input=Li2.inp Output=Li2.out Initial command: - /share/apps/gaussian/g09d01/nehalem/g09/l1.exe "/mnt/beegfs/tmpdir/42176/Gau-48548.inp" -scrdir="/mnt/beegfs/tmpdir/42176/" - Entering Link 1 = /share/apps/gaussian/g09d01/nehalem/g09/l1.exe PID= 48549. + /share/apps/gaussian/g09d01/nehalem/g09/l1.exe "/mnt/beegfs/tmpdir/42179/Gau-61549.inp" -scrdir="/mnt/beegfs/tmpdir/42179/" + Entering Link 1 = /share/apps/gaussian/g09d01/nehalem/g09/l1.exe PID= 61550. Copyright (c) 1988,1990,1992,1993,1995,1998,2003,2009,2013, Gaussian, Inc. All Rights Reserved. @@ -76,21 +76,21 @@ ****************************************** Gaussian 09: ES64L-G09RevD.01 24-Apr-2013 - 30-Mar-2019 + 1-Apr-2019 ****************************************** - ------------------------------------- - #p ROCCSD(T) cc-pVQZ pop=full gfprint - ------------------------------------- + ------------------------------------------------------------- + #p ROCCSD(T,FreezeInnerNobleGasCore) cc-pVQZ pop=full gfprint + ------------------------------------------------------------- 1/38=1/1; 2/12=2,17=6,18=5,40=1/2; 3/5=16,6=2,11=2,16=1,24=100,25=1,30=1,116=101/1,2,3; 4//1; 5/5=2,38=5/2; - 8/5=-1,6=4,9=120000,10=1/1,4; + 8/5=-1,6=4,9=120000,10=3/1,4; 9/5=7,14=2/13; 6/7=3/1; 99/5=1,9=1/99; - Leave Link 1 at Sat Mar 30 12:16:13 2019, MaxMem= 0 cpu: 0.0 + Leave Link 1 at Mon Apr 1 08:37:51 2019, MaxMem= 0 cpu: 0.0 (Enter /share/apps/gaussian/g09d01/nehalem/g09/l101.exe) -- G2 @@ -116,7 +116,7 @@ NQMom= -4.0100000 -4.0100000 NMagM= 3.2564240 3.2564240 AtZNuc= 3.0000000 3.0000000 - Leave Link 101 at Sat Mar 30 12:16:13 2019, MaxMem= 33554432 cpu: 0.0 + Leave Link 101 at Mon Apr 1 08:37:51 2019, MaxMem= 33554432 cpu: 0.1 (Enter /share/apps/gaussian/g09d01/nehalem/g09/l202.exe) Input orientation: --------------------------------------------------------------------- @@ -141,7 +141,7 @@ 2 3 0 0.000000 0.000000 -1.360640 --------------------------------------------------------------------- Rotational constants (GHZ): 0.0000000 19.4540897 19.4540897 - Leave Link 202 at Sat Mar 30 12:16:13 2019, MaxMem= 33554432 cpu: 0.0 + Leave Link 202 at Mon Apr 1 08:37:51 2019, MaxMem= 33554432 cpu: 0.0 (Enter /share/apps/gaussian/g09d01/nehalem/g09/l301.exe) Standard basis: CC-pVQZ (5D, 7F) Ernie: Thresh= 0.10000D-02 Tol= 0.10000D-05 Strict=F. @@ -261,7 +261,7 @@ Integral buffers will be 131072 words long. Raffenetti 2 integral format. Two-electron integral symmetry is turned on. - Leave Link 301 at Sat Mar 30 12:16:13 2019, MaxMem= 33554432 cpu: 0.0 + Leave Link 301 at Mon Apr 1 08:37:51 2019, MaxMem= 33554432 cpu: 0.0 (Enter /share/apps/gaussian/g09d01/nehalem/g09/l302.exe) NPDir=0 NMtPBC= 1 NCelOv= 1 NCel= 1 NClECP= 1 NCelD= 1 NCelK= 1 NCelE2= 1 NClLst= 1 CellRange= 0.0. @@ -269,10 +269,10 @@ One-electron integral symmetry used in STVInt NBasis= 110 RedAO= T EigKep= 4.29D-04 NBF= 22 7 13 13 7 22 13 13 NBsUse= 110 1.00D-06 EigRej= -1.00D+00 NBFU= 22 7 13 13 7 22 13 13 - Leave Link 302 at Sat Mar 30 12:16:13 2019, MaxMem= 33554432 cpu: 0.1 + Leave Link 302 at Mon Apr 1 08:37:51 2019, MaxMem= 33554432 cpu: 0.2 (Enter /share/apps/gaussian/g09d01/nehalem/g09/l303.exe) DipDrv: MaxL=1. - Leave Link 303 at Sat Mar 30 12:16:13 2019, MaxMem= 33554432 cpu: 0.0 + Leave Link 303 at Mon Apr 1 08:37:51 2019, MaxMem= 33554432 cpu: 0.0 (Enter /share/apps/gaussian/g09d01/nehalem/g09/l401.exe) ExpMin= 1.80D-02 ExpMax= 6.60D+03 ExpMxC= 6.43D+01 IAcc=3 IRadAn= 5 AccDes= 0.00D+00 Harris functional with IExCor= 205 and IRadAn= 5 diagonalized for initial guess. @@ -302,7 +302,7 @@ (PIG) (PIG) (PIG) (PIG) (SGU) (PIU) (PIU) (DLTU) (DLTU) (SGU) (SGG) (PIG) (PIG) (SGG) (SGU) The electronic state of the initial guess is 1-SGG. - Leave Link 401 at Sat Mar 30 12:16:13 2019, MaxMem= 33554432 cpu: 0.2 + Leave Link 401 at Mon Apr 1 08:37:52 2019, MaxMem= 33554432 cpu: 0.3 (Enter /share/apps/gaussian/g09d01/nehalem/g09/l502.exe) Restricted open shell SCF: Using DIIS extrapolation, IDIIS= 1040. @@ -438,7 +438,7 @@ KE= 1.488198932834D+01 PE=-3.794812619041D+01 EE= 6.444277515144D+00 Annihilation of the first spin contaminant: S**2 before annihilation 0.0000, after 0.0000 - Leave Link 502 at Sat Mar 30 12:16:15 2019, MaxMem= 33554432 cpu: 1.9 + Leave Link 502 at Mon Apr 1 08:37:54 2019, MaxMem= 33554432 cpu: 2.2 (Enter /share/apps/gaussian/g09d01/nehalem/g09/l801.exe) Windowed orbitals will be sorted by symmetry type. GenMOA: NOpAll= 8 NOp2=8 NOpUse= 8 JSym2X=1 @@ -450,66 +450,57 @@ NMat0= 1 NMatS0= 1 NMatT0= 0 NMatD0= 1 NMtDS0= 0 NMtDT0= 0 Petite list used in FoFCou. = 0.0000 = 0.0000 = 0.0000 = 0.0000 S= 0.0000 - ExpMin= 1.80D-02 ExpMax= 6.60D+03 ExpMxC= 6.43D+01 IAcc=3 IRadAn= 5 AccDes= 0.00D+00 - HarFok: IExCor= 205 AccDes= 0.00D+00 IRadAn= 5 IDoV=-2 UseB2=F ITyADJ=14 - ICtDFT= 12500011 ScaDFX= 1.000000 1.000000 1.000000 1.000000 - DSYEVD-2 returned Info= 1623 IAlg= 4 N= 55 NDim= 55 NE2= 609099 trying DSYEV. - Largest valence mixing into a core orbital is 6.69D-04 - Largest core mixing into a valence orbital is 2.64D-04 - Largest valence mixing into a core orbital is 6.69D-04 - Largest core mixing into a valence orbital is 2.64D-04 - Range of M.O.s used for correlation: 3 110 - NBasis= 110 NAE= 3 NBE= 3 NFC= 2 NFV= 0 - NROrb= 108 NOA= 1 NOB= 1 NVA= 107 NVB= 107 + Range of M.O.s used for correlation: 1 110 + NBasis= 110 NAE= 3 NBE= 3 NFC= 0 NFV= 0 + NROrb= 110 NOA= 3 NOB= 3 NVA= 107 NVB= 107 **** Warning!!: The largest alpha MO coefficient is 0.14662816D+02 **** Warning!!: The largest beta MO coefficient is 0.14662816D+02 - Singles contribution to E2= -0.2281768523D-16 - Leave Link 801 at Sat Mar 30 12:16:17 2019, MaxMem= 33554432 cpu: 1.3 + Singles contribution to E2= -0.2373569626D-16 + Leave Link 801 at Mon Apr 1 08:37:56 2019, MaxMem= 33554432 cpu: 1.3 (Enter /share/apps/gaussian/g09d01/nehalem/g09/l804.exe) Open-shell transformation, MDV= 33554432 ITran=4 ISComp=2. Semi-Direct transformation. - ModeAB= 4 MOrb= 1 LenV= 32985718 - LASXX= 79572 LTotXX= 79572 LenRXX= 162186 - LTotAB= 82614 MaxLAS= 591192 LenRXY= 0 - NonZer= 241758 LenScr= 720896 LnRSAI= 591192 - LnScr1= 1441792 LExtra= 0 Total= 2916066 + ModeAB= 4 MOrb= 3 LenV= 32978604 + LASXX= 250906 LTotXX= 250906 LenRXX= 511596 + LTotAB= 260690 MaxLAS= 1806420 LenRXY= 0 + NonZer= 762502 LenScr= 1572864 LnRSAI= 1806420 + LnScr1= 3276800 LExtra= 0 Total= 7167680 MaxDsk= -1 SrtSym= T ITran= 4 DoSDTr: NPSUse= 1 - JobTyp=1 Pass 1: I= 1 to 1. + JobTyp=1 Pass 1: I= 1 to 3. (rs|ai) integrals will be sorted in core. Complete sort for first half transformation. First half transformation complete. Complete sort for second half transformation. Second half transformation complete. - ModeAB= 4 MOrb= 1 LenV= 32985718 - LASXX= 79572 LTotXX= 79572 LenRXX= 160530 - LTotAB= 80958 MaxLAS= 591192 LenRXY= 0 - NonZer= 240102 LenScr= 720896 LnRSAI= 591192 - LnScr1= 1441792 LExtra= 0 Total= 2914410 + ModeAB= 4 MOrb= 3 LenV= 32978604 + LASXX= 250906 LTotXX= 250906 LenRXX= 496309 + LTotAB= 245403 MaxLAS= 1806420 LenRXY= 0 + NonZer= 747215 LenScr= 1572864 LnRSAI= 1806420 + LnScr1= 3276800 LExtra= 0 Total= 7152393 MaxDsk= -1 SrtSym= T ITran= 4 DoSDTr: NPSUse= 1 - JobTyp=2 Pass 1: I= 1 to 1. + JobTyp=2 Pass 1: I= 1 to 3. (rs|ai) integrals will be sorted in core. Complete sort for first half transformation. First half transformation complete. Complete sort for second half transformation. Second half transformation complete. Spin components of T(2) and E(2): - alpha-alpha T2 = 0.0000000000D+00 E2= 0.0000000000D+00 - alpha-beta T2 = 0.3401366332D-01 E2= -0.2208696457D-01 - beta-beta T2 = 0.0000000000D+00 E2= 0.0000000000D+00 - ANorm= 0.1016864624D+01 - E2 = -0.2208696457D-01 EUMP2 = -0.14893815819656D+02 + alpha-alpha T2 = 0.1164740706D-03 E2= -0.4855718283D-03 + alpha-beta T2 = 0.3759439006D-01 E2= -0.5342057546D-01 + beta-beta T2 = 0.1164740706D-03 E2= -0.4855718283D-03 + ANorm= 0.1018738111D+01 + E2 = -0.5439171912D-01 EUMP2 = -0.14926120574205D+02 (S**2,0)= 0.00000D+00 (S**2,1)= 0.00000D+00 - E(PUHF)= -0.14871728855D+02 E(PMP2)= -0.14893815820D+02 - Leave Link 804 at Sat Mar 30 12:16:21 2019, MaxMem= 33554432 cpu: 4.4 + E(PUHF)= -0.14871728855D+02 E(PMP2)= -0.14926120574D+02 + Leave Link 804 at Mon Apr 1 08:38:02 2019, MaxMem= 33554432 cpu: 6.0 (Enter /share/apps/gaussian/g09d01/nehalem/g09/l913.exe) CIDS: MDV= 33554432. - Frozen-core window: NFC= 2 NFV= 0. IFCWin=0 IBDFC=1 NFBD= 0 0 NFCmp= 0 0 NFFFC= 0 0 Using original routines for 1st iteration, S=T. Using DD4UQ or CC4UQ for 2nd and later iterations. @@ -524,115 +515,131 @@ Iterations= 50 Convergence= 0.100D-06 Iteration Nr. 1 ********************** - DD1Dir will call FoFMem 1 times, MxPair= 2 - NAB= 1 NAA= 0 NBB= 0. - DD1Dir will call FoFMem 1 times, MxPair= 2 - NAB= 1 NAA= 0 NBB= 0. - MP4(R+Q)= 0.55146737D-02 + DD1Dir will call FoFMem 1 times, MxPair= 24 + NAB= 9 NAA= 3 NBB= 3. + DD1Dir will call FoFMem 1 times, MxPair= 24 + NAB= 9 NAA= 3 NBB= 3. + MP4(R+Q)= 0.98325051D-02 Maximum subspace dimension= 5 - Norm of the A-vectors is 7.6873760D-03 conv= 1.00D-05. - RLE energy= -0.0213604185 - E3= -0.47634151D-02 EROMP3= -0.14898579235D+02 - E4(SDQ)= -0.26495728D-02 ROMP4(SDQ)= -0.14901228808D+02 + Norm of the A-vectors is 9.7801446D-03 conv= 1.00D-05. + RLE energy= -0.0534649533 + E3= -0.88896746D-02 EROMP3= -0.14935010249D+02 + E4(SDQ)= -0.31238862D-02 ROMP4(SDQ)= -0.14938134135D+02 VARIATIONAL ENERGIES WITH THE FIRST-ORDER WAVEFUNCTION: - DE(Corr)= -0.21335706E-01 E(Corr)= -14.893064561 - NORM(A)= 0.10156028D+01 + DE(Corr)= -0.53448889E-01 E(Corr)= -14.925177744 + NORM(A)= 0.10172529D+01 Iteration Nr. 2 ********************** - DD1Dir will call FoFMem 1 times, MxPair= 2 - NAB= 1 NAA= 0 NBB= 0. - Norm of the A-vectors is 7.1969522D-02 conv= 1.00D-05. - RLE energy= -0.0216425617 - DE(Corr)= -0.25924668E-01 E(CORR)= -14.897653523 Delta=-4.59D-03 - NORM(A)= 0.10161903D+01 + DD1Dir will call FoFMem 1 times, MxPair= 24 + NAB= 9 NAA= 3 NBB= 3. + Norm of the A-vectors is 9.5650662D-02 conv= 1.00D-05. + RLE energy= -0.0550659774 + DE(Corr)= -0.62133678E-01 E(CORR)= -14.933862533 Delta=-8.68D-03 + NORM(A)= 0.10190667D+01 Iteration Nr. 3 ********************** - DD1Dir will call FoFMem 1 times, MxPair= 2 - NAB= 1 NAA= 0 NBB= 0. - Norm of the A-vectors is 6.9982664D-02 conv= 1.00D-05. - RLE energy= -0.0248325931 - DE(Corr)= -0.26107640E-01 E(CORR)= -14.897836495 Delta=-1.83D-04 - NORM(A)= 0.10237008D+01 + DD1Dir will call FoFMem 1 times, MxPair= 24 + NAB= 9 NAA= 3 NBB= 3. + Norm of the A-vectors is 8.5991169D-02 conv= 1.00D-05. + RLE energy= -0.0567895820 + DE(Corr)= -0.62782025E-01 E(CORR)= -14.934510880 Delta=-6.48D-04 + NORM(A)= 0.10218889D+01 Iteration Nr. 4 ********************** - DD1Dir will call FoFMem 1 times, MxPair= 2 - NAB= 1 NAA= 0 NBB= 0. - Norm of the A-vectors is 4.9584809D-02 conv= 1.00D-05. - RLE energy= 0.0150330347 - DE(Corr)= -0.27852083E-01 E(CORR)= -14.899580938 Delta=-1.74D-03 - NORM(A)= 0.10965681D+01 + DD1Dir will call FoFMem 1 times, MxPair= 24 + NAB= 9 NAA= 3 NBB= 3. + Norm of the A-vectors is 7.5770898D-02 conv= 1.00D-05. + RLE energy= -0.4455317343 + DE(Corr)= -0.63484404E-01 E(CORR)= -14.935213260 Delta=-7.02D-04 + NORM(A)= 0.47380196D+01 Iteration Nr. 5 ********************** - DD1Dir will call FoFMem 1 times, MxPair= 2 - NAB= 1 NAA= 0 NBB= 0. - Norm of the A-vectors is 2.6786384D-01 conv= 1.00D-05. - RLE energy= -0.0254620713 - DE(Corr)= -0.19467389E-02 E(CORR)= -14.873675594 Delta= 2.59D-02 - NORM(A)= 0.10258493D+01 + DD1Dir will call FoFMem 1 times, MxPair= 24 + NAB= 9 NAA= 3 NBB= 3. + Norm of the A-vectors is 4.9259173D+00 conv= 1.00D-05. + RLE energy= -0.0061400624 + DE(Corr)= -0.69960765E-01 E(CORR)= -14.941689621 Delta=-6.48D-03 + NORM(A)= 0.10522188D+01 Iteration Nr. 6 ********************** - DD1Dir will call FoFMem 1 times, MxPair= 2 - NAB= 1 NAA= 0 NBB= 0. - Norm of the A-vectors is 4.5042192D-02 conv= 1.00D-05. - RLE energy= -0.0330316563 - DE(Corr)= -0.28258836E-01 E(CORR)= -14.899987691 Delta=-2.63D-02 - NORM(A)= 0.10552789D+01 + DD1Dir will call FoFMem 1 times, MxPair= 24 + NAB= 9 NAA= 3 NBB= 3. + Norm of the A-vectors is 3.5697592D-01 conv= 1.00D-05. + RLE energy= -0.0662615660 + DE(Corr)= -0.40506157E-01 E(CORR)= -14.912235012 Delta= 2.95D-02 + NORM(A)= 0.10404569D+01 Iteration Nr. 7 ********************** - DD1Dir will call FoFMem 1 times, MxPair= 2 - NAB= 1 NAA= 0 NBB= 0. - Norm of the A-vectors is 8.2107374D-03 conv= 1.00D-05. - RLE energy= -0.0318283206 - DE(Corr)= -0.32493469E-01 E(CORR)= -14.904222324 Delta=-4.23D-03 - NORM(A)= 0.10495588D+01 + DD1Dir will call FoFMem 1 times, MxPair= 24 + NAB= 9 NAA= 3 NBB= 3. + Norm of the A-vectors is 1.8172244D-02 conv= 1.00D-05. + RLE energy= -0.0665816060 + DE(Corr)= -0.67482825E-01 E(CORR)= -14.939211680 Delta=-2.70D-02 + NORM(A)= 0.10417437D+01 Iteration Nr. 8 ********************** - DD1Dir will call FoFMem 1 times, MxPair= 2 - NAB= 1 NAA= 0 NBB= 0. - Norm of the A-vectors is 9.1553446D-04 conv= 1.00D-05. - RLE energy= -0.0319613803 - DE(Corr)= -0.31881330E-01 E(CORR)= -14.903610185 Delta= 6.12D-04 - NORM(A)= 0.10501870D+01 + DD1Dir will call FoFMem 1 times, MxPair= 24 + NAB= 9 NAA= 3 NBB= 3. + Norm of the A-vectors is 1.5639926D-02 conv= 1.00D-05. + RLE energy= -0.0673831567 + DE(Corr)= -0.67643863E-01 E(CORR)= -14.939372718 Delta=-1.61D-04 + NORM(A)= 0.10449693D+01 Iteration Nr. 9 ********************** - DD1Dir will call FoFMem 1 times, MxPair= 2 - NAB= 1 NAA= 0 NBB= 0. - Norm of the A-vectors is 1.9823672D-04 conv= 1.00D-05. - RLE energy= -0.0319722299 - DE(Corr)= -0.31954686E-01 E(CORR)= -14.903683541 Delta=-7.34D-05 - NORM(A)= 0.10502354D+01 + DD1Dir will call FoFMem 1 times, MxPair= 24 + NAB= 9 NAA= 3 NBB= 3. + Norm of the A-vectors is 9.7460080D-03 conv= 1.00D-05. + RLE energy= -0.0684042917 + DE(Corr)= -0.68041817E-01 E(CORR)= -14.939770673 Delta=-3.98D-04 + NORM(A)= 0.10493494D+01 Iteration Nr. 10 ********************** - DD1Dir will call FoFMem 1 times, MxPair= 2 - NAB= 1 NAA= 0 NBB= 0. - Norm of the A-vectors is 1.9640378D-04 conv= 1.00D-05. - RLE energy= -0.0319490165 - DE(Corr)= -0.31961161E-01 E(CORR)= -14.903690016 Delta=-6.47D-06 - NORM(A)= 0.10501224D+01 + DD1Dir will call FoFMem 1 times, MxPair= 24 + NAB= 9 NAA= 3 NBB= 3. + Norm of the A-vectors is 2.1700507D-03 conv= 1.00D-05. + RLE energy= -0.0686821972 + DE(Corr)= -0.68549142E-01 E(CORR)= -14.940277997 Delta=-5.07D-04 + NORM(A)= 0.10505889D+01 Iteration Nr. 11 ********************** - DD1Dir will call FoFMem 1 times, MxPair= 2 - NAB= 1 NAA= 0 NBB= 0. - Norm of the A-vectors is 1.9804701D-05 conv= 1.00D-05. - RLE energy= -0.0319489107 - DE(Corr)= -0.31948999E-01 E(CORR)= -14.903677854 Delta= 1.22D-05 - NORM(A)= 0.10501218D+01 + DD1Dir will call FoFMem 1 times, MxPair= 24 + NAB= 9 NAA= 3 NBB= 3. + Norm of the A-vectors is 8.9501224D-05 conv= 1.00D-05. + RLE energy= -0.0686845807 + DE(Corr)= -0.68684253E-01 E(CORR)= -14.940413108 Delta=-1.35D-04 + NORM(A)= 0.10505985D+01 Iteration Nr. 12 ********************** - DD1Dir will call FoFMem 1 times, MxPair= 2 - NAB= 1 NAA= 0 NBB= 0. - Norm of the A-vectors is 6.9131188D-06 conv= 1.00D-05. - RLE energy= -0.0319490734 - DE(Corr)= -0.31949017E-01 E(CORR)= -14.903677872 Delta=-1.79D-08 - NORM(A)= 0.10501224D+01 - CI/CC converged in 12 iterations to DelEn=-1.79D-08 Conv= 1.00D-07 ErrA1= 6.91D-06 Conv= 1.00D-05 - Largest amplitude= 8.20D-02 - Time for triples= 0.00 seconds. - T4(CCSD)= 0.00000000D+00 - T5(CCSD)= 0.00000000D+00 - CCSD(T)= -0.14903677872D+02 + DD1Dir will call FoFMem 1 times, MxPair= 24 + NAB= 9 NAA= 3 NBB= 3. + Norm of the A-vectors is 2.8109103D-05 conv= 1.00D-05. + RLE energy= -0.0686857858 + DE(Corr)= -0.68685315E-01 E(CORR)= -14.940414170 Delta=-1.06D-06 + NORM(A)= 0.10506027D+01 + Iteration Nr. 13 + ********************** + DD1Dir will call FoFMem 1 times, MxPair= 24 + NAB= 9 NAA= 3 NBB= 3. + Norm of the A-vectors is 9.7913017D-06 conv= 1.00D-05. + RLE energy= -0.0686857050 + DE(Corr)= -0.68685668E-01 E(CORR)= -14.940414523 Delta=-3.53D-07 + NORM(A)= 0.10506024D+01 + Iteration Nr. 14 + ********************** + DD1Dir will call FoFMem 1 times, MxPair= 24 + NAB= 9 NAA= 3 NBB= 3. + Norm of the A-vectors is 2.9603690D-06 conv= 1.00D-05. + RLE energy= -0.0686857010 + DE(Corr)= -0.68685698E-01 E(CORR)= -14.940414553 Delta=-3.01D-08 + NORM(A)= 0.10506024D+01 + CI/CC converged in 14 iterations to DelEn=-3.01D-08 Conv= 1.00D-07 ErrA1= 2.96D-06 Conv= 1.00D-05 + Largest amplitude= 7.97D-02 + Time for triples= 62.48 seconds. + T4(CCSD)= -0.49202587D-03 + T5(CCSD)= 0.35203850D-05 + CCSD(T)= -0.14940903058D+02 Discarding MO integrals. - Leave Link 913 at Sat Mar 30 12:16:31 2019, MaxMem= 33554432 cpu: 3.8 + Leave Link 913 at Mon Apr 1 08:43:36 2019, MaxMem= 33554432 cpu: 71.5 (Enter /share/apps/gaussian/g09d01/nehalem/g09/l601.exe) Copying SCF densities to generalized density rwf, IOpCl= 0 IROHF=1. @@ -7339,18 +7346,18 @@ --------------------------------------------------------------------------------- No NMR shielding tensors so no spin-rotation constants. - Leave Link 601 at Sat Mar 30 12:16:31 2019, MaxMem= 33554432 cpu: 0.2 + Leave Link 601 at Mon Apr 1 08:43:37 2019, MaxMem= 33554432 cpu: 0.3 (Enter /share/apps/gaussian/g09d01/nehalem/g09/l9999.exe) - 1\1\GINC-COMPUTE-7-0\SP\ROCCSD(T)-FC\CC-pVQZ\Li2\LOOS\30-Mar-2019\0\\# - p ROCCSD(T) cc-pVQZ pop=full gfprint\\G2\\0,1\Li\Li,1,2.72127987\\Vers - ion=ES64L-G09RevD.01\State=1-SGG\HF=-14.8717289\MP2=-14.8938158\MP3=-1 - 4.8985792\PUHF=-14.8717289\PMP2-0=-14.8938158\MP4SDQ=-14.9012288\CCSD= - -14.9036779\CCSD(T)=-14.9036779\RMSD=1.558e-09\PG=D*H [C*(Li1.Li1)]\\@ + 1\1\GINC-COMPUTE-7-0\SP\ROCCSD(T)-FC1\CC-pVQZ\Li2\LOOS\01-Apr-2019\0\\ + #p ROCCSD(T,FreezeInnerNobleGasCore) cc-pVQZ pop=full gfprint\\G2\\0,1 + \Li\Li,1,2.72127987\\Version=ES64L-G09RevD.01\State=1-SGG\HF=-14.87172 + 89\MP2=-14.9261206\MP3=-14.9350102\PUHF=-14.8717289\PMP2-0=-14.9261206 + \MP4SDQ=-14.9381341\CCSD=-14.9404146\CCSD(T)=-14.9409031\RMSD=1.558e-0 + 9\PG=D*H [C*(Li1.Li1)]\\@ - THE MORE ACCURATE THE CALCULATIONS BECOME, THE MORE THE CONCEPTS - TEND TO VANISH INTO THIN AIR. - -- R.S. MULLIKEN, J.C.P. 43,S2(1965) - Job cpu time: 0 days 0 hours 0 minutes 12.2 seconds. - File lengths (MBytes): RWF= 66 Int= 0 D2E= 0 Chk= 2 Scr= 1 - Normal termination of Gaussian 09 at Sat Mar 30 12:16:31 2019. + NATURE REVEALS EVERY SECRET ONCE. + - RALPH WALDO EMERSON + Job cpu time: 0 days 0 hours 1 minutes 22.1 seconds. + File lengths (MBytes): RWF= 98 Int= 0 D2E= 0 Chk= 2 Scr= 1 + Normal termination of Gaussian 09 at Mon Apr 1 08:43:37 2019. diff --git a/G09/Mixed_core/Molecules/vqz/LiH.out b/G09/Mixed_core/Molecules/vqz/LiH.out index 6125559..5310661 100644 --- a/G09/Mixed_core/Molecules/vqz/LiH.out +++ b/G09/Mixed_core/Molecules/vqz/LiH.out @@ -2,8 +2,8 @@ Input=LiH.inp Output=LiH.out Initial command: - /share/apps/gaussian/g09d01/nehalem/g09/l1.exe "/mnt/beegfs/tmpdir/42176/Gau-48636.inp" -scrdir="/mnt/beegfs/tmpdir/42176/" - Entering Link 1 = /share/apps/gaussian/g09d01/nehalem/g09/l1.exe PID= 48637. + /share/apps/gaussian/g09d01/nehalem/g09/l1.exe "/mnt/beegfs/tmpdir/42179/Gau-61633.inp" -scrdir="/mnt/beegfs/tmpdir/42179/" + Entering Link 1 = /share/apps/gaussian/g09d01/nehalem/g09/l1.exe PID= 61634. Copyright (c) 1988,1990,1992,1993,1995,1998,2003,2009,2013, Gaussian, Inc. All Rights Reserved. @@ -76,21 +76,21 @@ ****************************************** Gaussian 09: ES64L-G09RevD.01 24-Apr-2013 - 30-Mar-2019 + 1-Apr-2019 ****************************************** - ------------------------------------- - #p ROCCSD(T) cc-pVQZ pop=full gfprint - ------------------------------------- + ------------------------------------------------------------- + #p ROCCSD(T,FreezeInnerNobleGasCore) cc-pVQZ pop=full gfprint + ------------------------------------------------------------- 1/38=1/1; 2/12=2,17=6,18=5,40=1/2; 3/5=16,6=2,11=2,16=1,24=100,25=1,30=1,116=101/1,2,3; 4//1; 5/5=2,38=5/2; - 8/5=-1,6=4,9=120000,10=1/1,4; + 8/5=-1,6=4,9=120000,10=3/1,4; 9/5=7,14=2/13; 6/7=3/1; 99/5=1,9=1/99; - Leave Link 1 at Sat Mar 30 12:34:31 2019, MaxMem= 0 cpu: 0.0 + Leave Link 1 at Mon Apr 1 08:54:46 2019, MaxMem= 0 cpu: 0.0 (Enter /share/apps/gaussian/g09d01/nehalem/g09/l101.exe) -- G2 @@ -116,7 +116,7 @@ NQMom= -4.0100000 0.0000000 NMagM= 3.2564240 2.7928460 AtZNuc= 3.0000000 1.0000000 - Leave Link 101 at Sat Mar 30 12:34:31 2019, MaxMem= 33554432 cpu: 0.1 + Leave Link 101 at Mon Apr 1 08:54:46 2019, MaxMem= 33554432 cpu: 0.1 (Enter /share/apps/gaussian/g09d01/nehalem/g09/l202.exe) Input orientation: --------------------------------------------------------------------- @@ -141,7 +141,7 @@ 2 1 0 0.000000 0.000000 -1.210897 --------------------------------------------------------------------- Rotational constants (GHZ): 0.0000000 220.0046195 220.0046195 - Leave Link 202 at Sat Mar 30 12:34:31 2019, MaxMem= 33554432 cpu: 0.0 + Leave Link 202 at Mon Apr 1 08:54:46 2019, MaxMem= 33554432 cpu: 0.0 (Enter /share/apps/gaussian/g09d01/nehalem/g09/l301.exe) Standard basis: CC-pVQZ (5D, 7F) Ernie: Thresh= 0.10000D-02 Tol= 0.10000D-05 Strict=F. @@ -231,17 +231,17 @@ Integral buffers will be 131072 words long. Raffenetti 2 integral format. Two-electron integral symmetry is turned on. - Leave Link 301 at Sat Mar 30 12:34:31 2019, MaxMem= 33554432 cpu: 0.0 + Leave Link 301 at Mon Apr 1 08:54:46 2019, MaxMem= 33554432 cpu: 0.1 (Enter /share/apps/gaussian/g09d01/nehalem/g09/l302.exe) NPDir=0 NMtPBC= 1 NCelOv= 1 NCel= 1 NClECP= 1 NCelD= 1 NCelK= 1 NCelE2= 1 NClLst= 1 CellRange= 0.0. One-electron integrals computed using PRISM. NBasis= 85 RedAO= T EigKep= 2.09D-03 NBF= 35 10 20 20 NBsUse= 85 1.00D-06 EigRej= -1.00D+00 NBFU= 35 10 20 20 - Leave Link 302 at Sat Mar 30 12:34:31 2019, MaxMem= 33554432 cpu: 0.1 + Leave Link 302 at Mon Apr 1 08:54:46 2019, MaxMem= 33554432 cpu: 0.1 (Enter /share/apps/gaussian/g09d01/nehalem/g09/l303.exe) DipDrv: MaxL=1. - Leave Link 303 at Sat Mar 30 12:34:31 2019, MaxMem= 33554432 cpu: 0.0 + Leave Link 303 at Mon Apr 1 08:54:46 2019, MaxMem= 33554432 cpu: 0.0 (Enter /share/apps/gaussian/g09d01/nehalem/g09/l401.exe) ExpMin= 1.80D-02 ExpMax= 6.60D+03 ExpMxC= 8.26D+01 IAcc=3 IRadAn= 5 AccDes= 0.00D+00 Harris functional with IExCor= 205 and IRadAn= 5 diagonalized for initial guess. @@ -267,7 +267,7 @@ (DLTA) (DLTA) (PI) (PI) (SG) (SG) (PI) (PI) (SG) (DLTA) (DLTA) (SG) (PI) (PI) The electronic state of the initial guess is 1-SG. - Leave Link 401 at Sat Mar 30 12:34:32 2019, MaxMem= 33554432 cpu: 0.2 + Leave Link 401 at Mon Apr 1 08:54:47 2019, MaxMem= 33554432 cpu: 0.2 (Enter /share/apps/gaussian/g09d01/nehalem/g09/l502.exe) Restricted open shell SCF: Using DIIS extrapolation, IDIIS= 1040. @@ -403,7 +403,7 @@ KE= 7.982129617549D+00 PE=-2.042752621232D+01 EE= 3.474937601709D+00 Annihilation of the first spin contaminant: S**2 before annihilation 0.0000, after 0.0000 - Leave Link 502 at Sat Mar 30 12:34:33 2019, MaxMem= 33554432 cpu: 1.4 + Leave Link 502 at Mon Apr 1 08:54:48 2019, MaxMem= 33554432 cpu: 1.5 (Enter /share/apps/gaussian/g09d01/nehalem/g09/l801.exe) Windowed orbitals will be sorted by symmetry type. GenMOA: NOpAll= 4 NOp2=4 NOpUse= 4 JSym2X=1 @@ -415,59 +415,51 @@ NMat0= 1 NMatS0= 1 NMatT0= 0 NMatD0= 1 NMtDS0= 0 NMtDT0= 0 Petite list used in FoFCou. = 0.0000 = 0.0000 = 0.0000 = 0.0000 S= 0.0000 - ExpMin= 1.80D-02 ExpMax= 6.60D+03 ExpMxC= 8.26D+01 IAcc=3 IRadAn= 5 AccDes= 0.00D+00 - HarFok: IExCor= 205 AccDes= 0.00D+00 IRadAn= 5 IDoV=-2 UseB2=F ITyADJ=14 - ICtDFT= 12500011 ScaDFX= 1.000000 1.000000 1.000000 1.000000 - Largest valence mixing into a core orbital is 6.45D-04 - Largest core mixing into a valence orbital is 9.55D-05 - Largest valence mixing into a core orbital is 6.45D-04 - Largest core mixing into a valence orbital is 9.55D-05 - Range of M.O.s used for correlation: 2 85 - NBasis= 85 NAE= 2 NBE= 2 NFC= 1 NFV= 0 - NROrb= 84 NOA= 1 NOB= 1 NVA= 83 NVB= 83 - Singles contribution to E2= -0.8521802631D-16 - Leave Link 801 at Sat Mar 30 12:34:34 2019, MaxMem= 33554432 cpu: 1.0 + Range of M.O.s used for correlation: 1 85 + NBasis= 85 NAE= 2 NBE= 2 NFC= 0 NFV= 0 + NROrb= 85 NOA= 2 NOB= 2 NVA= 83 NVB= 83 + Singles contribution to E2= -0.9145810081D-16 + Leave Link 801 at Mon Apr 1 08:54:49 2019, MaxMem= 33554432 cpu: 0.9 (Enter /share/apps/gaussian/g09d01/nehalem/g09/l804.exe) Open-shell transformation, MDV= 33554432 ITran=4 ISComp=2. Semi-Direct transformation. - ModeAB= 2 MOrb= 1 LenV= 33118976 - LASXX= 75138 LTotXX= 75138 LenRXX= 75138 - LTotAB= 78980 MaxLAS= 467460 LenRXY= 467460 - NonZer= 495684 LenScr= 1441792 LnRSAI= 0 - LnScr1= 0 LExtra= 0 Total= 1984390 + ModeAB= 2 MOrb= 2 LenV= 33116235 + LASXX= 155274 LTotXX= 155274 LenRXX= 155274 + LTotAB= 163245 MaxLAS= 946050 LenRXY= 946050 + NonZer= 1003170 LenScr= 2162688 LnRSAI= 0 + LnScr1= 0 LExtra= 0 Total= 3264012 MaxDsk= -1 SrtSym= F ITran= 4 DoSDTr: NPSUse= 1 - JobTyp=1 Pass 1: I= 1 to 1. + JobTyp=1 Pass 1: I= 1 to 2. (rs|ai) integrals will be sorted in core. Complete sort for first half transformation. First half transformation complete. Complete sort for second half transformation. Second half transformation complete. - ModeAB= 2 MOrb= 1 LenV= 33118976 - LASXX= 75138 LTotXX= 75138 LenRXX= 75138 - LTotAB= 76924 MaxLAS= 467460 LenRXY= 467460 - NonZer= 495684 LenScr= 1441792 LnRSAI= 0 - LnScr1= 0 LExtra= 0 Total= 1984390 + ModeAB= 2 MOrb= 2 LenV= 33116235 + LASXX= 155274 LTotXX= 155274 LenRXX= 946050 + LTotAB= 154884 MaxLAS= 946050 LenRXY= 154884 + NonZer= 1003170 LenScr= 2162688 LnRSAI= 0 + LnScr1= 0 LExtra= 0 Total= 3263622 MaxDsk= -1 SrtSym= F ITran= 4 DoSDTr: NPSUse= 1 - JobTyp=2 Pass 1: I= 1 to 1. + JobTyp=2 Pass 1: I= 1 to 2. (rs|ai) integrals will be sorted in core. Complete sort for first half transformation. First half transformation complete. Complete sort for second half transformation. Second half transformation complete. Spin components of T(2) and E(2): - alpha-alpha T2 = 0.0000000000D+00 E2= 0.0000000000D+00 - alpha-beta T2 = 0.1534863668D-01 E2= -0.2891514999D-01 - beta-beta T2 = 0.0000000000D+00 E2= 0.0000000000D+00 - ANorm= 0.1007645095D+01 - E2 = -0.2891514999D-01 EUMP2 = -0.80160960986076D+01 + alpha-alpha T2 = 0.6000134307D-04 E2= -0.2866376339D-03 + alpha-beta T2 = 0.1711797398D-01 E2= -0.4446848811D-01 + beta-beta T2 = 0.6000134307D-04 E2= -0.2866376339D-03 + ANorm= 0.1008582162D+01 + E2 = -0.4504176337D-01 EUMP2 = -0.80322227119962D+01 (S**2,0)= 0.00000D+00 (S**2,1)= 0.00000D+00 - E(PUHF)= -0.79871809486D+01 E(PMP2)= -0.80160960986D+01 - Leave Link 804 at Sat Mar 30 12:34:38 2019, MaxMem= 33554432 cpu: 3.5 + E(PUHF)= -0.79871809486D+01 E(PMP2)= -0.80322227120D+01 + Leave Link 804 at Mon Apr 1 08:54:54 2019, MaxMem= 33554432 cpu: 4.5 (Enter /share/apps/gaussian/g09d01/nehalem/g09/l913.exe) CIDS: MDV= 33554432. - Frozen-core window: NFC= 1 NFV= 0. IFCWin=0 IBDFC=1 NFBD= 0 0 NFCmp= 0 0 NFFFC= 0 0 Using original routines for 1st iteration, S=T. Using DD4UQ or CC4UQ for 2nd and later iterations. @@ -482,107 +474,107 @@ Iterations= 50 Convergence= 0.100D-06 Iteration Nr. 1 ********************** - DD1Dir will call FoFMem 1 times, MxPair= 2 - NAB= 1 NAA= 0 NBB= 0. - DD1Dir will call FoFMem 1 times, MxPair= 2 - NAB= 1 NAA= 0 NBB= 0. - MP4(R+Q)= 0.57708721D-02 + DD1Dir will call FoFMem 1 times, MxPair= 10 + NAB= 4 NAA= 1 NBB= 1. + DD1Dir will call FoFMem 1 times, MxPair= 10 + NAB= 4 NAA= 1 NBB= 1. + MP4(R+Q)= 0.79652696D-02 Maximum subspace dimension= 5 - Norm of the A-vectors is 2.7697089D-03 conv= 1.00D-05. - RLE energy= -0.0284780507 - E3= -0.53270640D-02 EROMP3= -0.80214231626D+01 - E4(SDQ)= -0.17157620D-02 ROMP4(SDQ)= -0.80231389246D+01 + Norm of the A-vectors is 3.3594639D-03 conv= 1.00D-05. + RLE energy= -0.0445364536 + E3= -0.74542266D-02 EROMP3= -0.80396769386D+01 + E4(SDQ)= -0.19961790D-02 ROMP4(SDQ)= -0.80416731176D+01 VARIATIONAL ENERGIES WITH THE FIRST-ORDER WAVEFUNCTION: - DE(Corr)= -0.28471342E-01 E(Corr)= -8.0156522905 - NORM(A)= 0.10073435D+01 + DE(Corr)= -0.44530720E-01 E(Corr)= -8.0317116690 + NORM(A)= 0.10082549D+01 Iteration Nr. 2 ********************** - DD1Dir will call FoFMem 1 times, MxPair= 2 - NAB= 1 NAA= 0 NBB= 0. - Norm of the A-vectors is 3.6547574D-02 conv= 1.00D-05. - RLE energy= -0.0288667090 - DE(Corr)= -0.33701614E-01 E(CORR)= -8.0208825625 Delta=-5.23D-03 - NORM(A)= 0.10076153D+01 + DD1Dir will call FoFMem 1 times, MxPair= 10 + NAB= 4 NAA= 1 NBB= 1. + Norm of the A-vectors is 4.7321597D-02 conv= 1.00D-05. + RLE energy= -0.0453785124 + DE(Corr)= -0.51877799E-01 E(CORR)= -8.0390587475 Delta=-7.35D-03 + NORM(A)= 0.10086989D+01 Iteration Nr. 3 ********************** - DD1Dir will call FoFMem 1 times, MxPair= 2 - NAB= 1 NAA= 0 NBB= 0. - Norm of the A-vectors is 3.4830862D-02 conv= 1.00D-05. - RLE energy= -0.0314430112 - DE(Corr)= -0.33846406E-01 E(CORR)= -8.0210273549 Delta=-1.45D-04 - NORM(A)= 0.10096111D+01 + DD1Dir will call FoFMem 1 times, MxPair= 10 + NAB= 4 NAA= 1 NBB= 1. + Norm of the A-vectors is 4.3660847D-02 conv= 1.00D-05. + RLE energy= -0.0475731641 + DE(Corr)= -0.52137636E-01 E(CORR)= -8.0393185849 Delta=-2.60D-04 + NORM(A)= 0.10100370D+01 Iteration Nr. 4 ********************** - DD1Dir will call FoFMem 1 times, MxPair= 2 - NAB= 1 NAA= 0 NBB= 0. - Norm of the A-vectors is 2.3793039D-02 conv= 1.00D-05. - RLE energy= -0.0346443770 - DE(Corr)= -0.34774565E-01 E(CORR)= -8.0219555140 Delta=-9.28D-04 - NORM(A)= 0.10126805D+01 + DD1Dir will call FoFMem 1 times, MxPair= 10 + NAB= 4 NAA= 1 NBB= 1. + Norm of the A-vectors is 3.4256936D-02 conv= 1.00D-05. + RLE energy= -0.0552266528 + DE(Corr)= -0.52822446E-01 E(CORR)= -8.0400033944 Delta=-6.85D-04 + NORM(A)= 0.10159919D+01 Iteration Nr. 5 ********************** - DD1Dir will call FoFMem 1 times, MxPair= 2 - NAB= 1 NAA= 0 NBB= 0. - Norm of the A-vectors is 9.5539492D-03 conv= 1.00D-05. - RLE energy= -0.0335422402 - DE(Corr)= -0.35966451E-01 E(CORR)= -8.0231473999 Delta=-1.19D-03 - NORM(A)= 0.10116514D+01 + DD1Dir will call FoFMem 1 times, MxPair= 10 + NAB= 4 NAA= 1 NBB= 1. + Norm of the A-vectors is 1.9798167D-03 conv= 1.00D-05. + RLE energy= -0.0544401998 + DE(Corr)= -0.55255763E-01 E(CORR)= -8.0424367120 Delta=-2.43D-03 + NORM(A)= 0.10154192D+01 Iteration Nr. 6 ********************** - DD1Dir will call FoFMem 1 times, MxPair= 2 - NAB= 1 NAA= 0 NBB= 0. - Norm of the A-vectors is 1.4277513D-02 conv= 1.00D-05. - RLE energy= -0.0367642958 - DE(Corr)= -0.35570829E-01 E(CORR)= -8.0227517775 Delta= 3.96D-04 - NORM(A)= 0.10152043D+01 + DD1Dir will call FoFMem 1 times, MxPair= 10 + NAB= 4 NAA= 1 NBB= 1. + Norm of the A-vectors is 3.9700148D-03 conv= 1.00D-05. + RLE energy= -0.0552803576 + DE(Corr)= -0.55016601E-01 E(CORR)= -8.0421975500 Delta= 2.39D-04 + NORM(A)= 0.10162474D+01 Iteration Nr. 7 ********************** - DD1Dir will call FoFMem 1 times, MxPair= 2 - NAB= 1 NAA= 0 NBB= 0. - Norm of the A-vectors is 3.1086415D-04 conv= 1.00D-05. - RLE energy= -0.0367331129 - DE(Corr)= -0.36746054E-01 E(CORR)= -8.0239270025 Delta=-1.18D-03 - NORM(A)= 0.10151633D+01 + DD1Dir will call FoFMem 1 times, MxPair= 10 + NAB= 4 NAA= 1 NBB= 1. + Norm of the A-vectors is 3.2159425D-04 conv= 1.00D-05. + RLE energy= -0.0552876915 + DE(Corr)= -0.55283786E-01 E(CORR)= -8.0424647344 Delta=-2.67D-04 + NORM(A)= 0.10162700D+01 Iteration Nr. 8 ********************** - DD1Dir will call FoFMem 1 times, MxPair= 2 - NAB= 1 NAA= 0 NBB= 0. - Norm of the A-vectors is 9.6565027D-05 conv= 1.00D-05. - RLE energy= -0.0367348331 - DE(Corr)= -0.36735207E-01 E(CORR)= -8.0239161556 Delta= 1.08D-05 - NORM(A)= 0.10151654D+01 + DD1Dir will call FoFMem 1 times, MxPair= 10 + NAB= 4 NAA= 1 NBB= 1. + Norm of the A-vectors is 1.2019544D-04 conv= 1.00D-05. + RLE energy= -0.0552860497 + DE(Corr)= -0.55287481E-01 E(CORR)= -8.0424684300 Delta=-3.70D-06 + NORM(A)= 0.10162679D+01 Iteration Nr. 9 ********************** - DD1Dir will call FoFMem 1 times, MxPair= 2 - NAB= 1 NAA= 0 NBB= 0. - Norm of the A-vectors is 4.4921349D-05 conv= 1.00D-05. - RLE energy= -0.0367361590 - DE(Corr)= -0.36735522E-01 E(CORR)= -8.0239164702 Delta=-3.15D-07 - NORM(A)= 0.10151686D+01 + DD1Dir will call FoFMem 1 times, MxPair= 10 + NAB= 4 NAA= 1 NBB= 1. + Norm of the A-vectors is 5.0725875D-05 conv= 1.00D-05. + RLE energy= -0.0552864750 + DE(Corr)= -0.55286192E-01 E(CORR)= -8.0424671408 Delta= 1.29D-06 + NORM(A)= 0.10162678D+01 Iteration Nr. 10 ********************** - DD1Dir will call FoFMem 1 times, MxPair= 2 - NAB= 1 NAA= 0 NBB= 0. - Norm of the A-vectors is 9.3225217D-06 conv= 1.00D-05. - RLE energy= -0.0367361736 - DE(Corr)= -0.36736120E-01 E(CORR)= -8.0239170691 Delta=-5.99D-07 - NORM(A)= 0.10151691D+01 + DD1Dir will call FoFMem 1 times, MxPair= 10 + NAB= 4 NAA= 1 NBB= 1. + Norm of the A-vectors is 1.1832212D-05 conv= 1.00D-05. + RLE energy= -0.0552863826 + DE(Corr)= -0.55286403E-01 E(CORR)= -8.0424673520 Delta=-2.11D-07 + NORM(A)= 0.10162680D+01 Iteration Nr. 11 ********************** - DD1Dir will call FoFMem 1 times, MxPair= 2 - NAB= 1 NAA= 0 NBB= 0. - Norm of the A-vectors is 3.1481720D-06 conv= 1.00D-05. - RLE energy= -0.0367361959 - DE(Corr)= -0.36736196E-01 E(CORR)= -8.0239171451 Delta=-7.60D-08 - NORM(A)= 0.10151690D+01 - CI/CC converged in 11 iterations to DelEn=-7.60D-08 Conv= 1.00D-07 ErrA1= 3.15D-06 Conv= 1.00D-05 + DD1Dir will call FoFMem 1 times, MxPair= 10 + NAB= 4 NAA= 1 NBB= 1. + Norm of the A-vectors is 4.0992894D-06 conv= 1.00D-05. + RLE energy= -0.0552863772 + DE(Corr)= -0.55286381E-01 E(CORR)= -8.0424673293 Delta= 2.27D-08 + NORM(A)= 0.10162680D+01 + CI/CC converged in 11 iterations to DelEn= 2.27D-08 Conv= 1.00D-07 ErrA1= 4.10D-06 Conv= 1.00D-05 Largest amplitude= 3.64D-02 - Time for triples= 0.00 seconds. - T4(CCSD)= 0.00000000D+00 - T5(CCSD)= 0.00000000D+00 - CCSD(T)= -0.80239171451D+01 + Time for triples= 11.69 seconds. + T4(CCSD)= -0.11638087D-03 + T5(CCSD)= 0.94387118D-06 + CCSD(T)= -0.80425827663D+01 Discarding MO integrals. - Leave Link 913 at Sat Mar 30 12:34:43 2019, MaxMem= 33554432 cpu: 2.4 + Leave Link 913 at Mon Apr 1 08:56:06 2019, MaxMem= 33554432 cpu: 15.5 (Enter /share/apps/gaussian/g09d01/nehalem/g09/l601.exe) Copying SCF densities to generalized density rwf, IOpCl= 0 IROHF=1. @@ -4720,25 +4712,17 @@ --------------------------------------------------------------------------------- No NMR shielding tensors so no spin-rotation constants. - Leave Link 601 at Sat Mar 30 12:34:43 2019, MaxMem= 33554432 cpu: 0.1 + Leave Link 601 at Mon Apr 1 08:56:07 2019, MaxMem= 33554432 cpu: 0.2 (Enter /share/apps/gaussian/g09d01/nehalem/g09/l9999.exe) - 1\1\GINC-COMPUTE-7-0\SP\ROCCSD(T)-FC\CC-pVQZ\H1Li1\LOOS\30-Mar-2019\0\ - \#p ROCCSD(T) cc-pVQZ pop=full gfprint\\G2\\0,1\Li\H,1,1.61452972\\Ver - sion=ES64L-G09RevD.01\State=1-SG\HF=-7.9871809\MP2=-8.0160961\MP3=-8.0 - 214232\PUHF=-7.9871809\PMP2-0=-8.0160961\MP4SDQ=-8.0231389\CCSD=-8.023 - 9171\CCSD(T)=-8.0239171\RMSD=3.547e-09\PG=C*V [C*(H1Li1)]\\@ + 1\1\GINC-COMPUTE-7-0\SP\ROCCSD(T)-FC1\CC-pVQZ\H1Li1\LOOS\01-Apr-2019\0 + \\#p ROCCSD(T,FreezeInnerNobleGasCore) cc-pVQZ pop=full gfprint\\G2\\0 + ,1\Li\H,1,1.61452972\\Version=ES64L-G09RevD.01\State=1-SG\HF=-7.987180 + 9\MP2=-8.0322227\MP3=-8.0396769\PUHF=-7.9871809\PMP2-0=-8.0322227\MP4S + DQ=-8.0416731\CCSD=-8.0424673\CCSD(T)=-8.0425828\RMSD=3.547e-09\PG=C*V + [C*(H1Li1)]\\@ - UPON JULIA'S CLOTHES - - WHENAS IN SILKS MY JULIA GOES, - THEN, THEN, METHINKS, HOW SWEETLY FLOWS - THAT LIQUEFACTION OF HER CLOTHES. - - NEXT, WHEN I CAST MINE EYES, AND SEE - THAT BRAVE VIBRATION, EACH WAY FREE, - O, HOW THAT GLITTERING TAKETH ME! - -- ROBERT HERRICK, 1648 - Job cpu time: 0 days 0 hours 0 minutes 9.1 seconds. - File lengths (MBytes): RWF= 63 Int= 0 D2E= 0 Chk= 2 Scr= 1 - Normal termination of Gaussian 09 at Sat Mar 30 12:34:43 2019. + ASKING DUMB QUESTIONS IS EASIER THAN CORECTING DUMB MISTAKES. + Job cpu time: 0 days 0 hours 0 minutes 23.5 seconds. + File lengths (MBytes): RWF= 111 Int= 0 D2E= 0 Chk= 2 Scr= 1 + Normal termination of Gaussian 09 at Mon Apr 1 08:56:07 2019. diff --git a/G09/Mixed_core/Molecules/vqz/Na2.out b/G09/Mixed_core/Molecules/vqz/Na2.out index 9a029e7..7ccea4a 100644 --- a/G09/Mixed_core/Molecules/vqz/Na2.out +++ b/G09/Mixed_core/Molecules/vqz/Na2.out @@ -2,8 +2,8 @@ Input=Na2.inp Output=Na2.out Initial command: - /share/apps/gaussian/g09d01/nehalem/g09/l1.exe "/mnt/beegfs/tmpdir/42176/Gau-49388.inp" -scrdir="/mnt/beegfs/tmpdir/42176/" - Entering Link 1 = /share/apps/gaussian/g09d01/nehalem/g09/l1.exe PID= 49389. + /share/apps/gaussian/g09d01/nehalem/g09/l1.exe "/mnt/beegfs/tmpdir/42179/Gau-62618.inp" -scrdir="/mnt/beegfs/tmpdir/42179/" + Entering Link 1 = /share/apps/gaussian/g09d01/nehalem/g09/l1.exe PID= 62619. Copyright (c) 1988,1990,1992,1993,1995,1998,2003,2009,2013, Gaussian, Inc. All Rights Reserved. @@ -76,21 +76,21 @@ ****************************************** Gaussian 09: ES64L-G09RevD.01 24-Apr-2013 - 30-Mar-2019 + 1-Apr-2019 ****************************************** - ------------------------------------- - #p ROCCSD(T) cc-pVQZ pop=full gfprint - ------------------------------------- + ------------------------------------------------------------- + #p ROCCSD(T,FreezeInnerNobleGasCore) cc-pVQZ pop=full gfprint + ------------------------------------------------------------- 1/38=1/1; 2/12=2,17=6,18=5,40=1/2; 3/5=16,6=2,11=2,16=1,24=100,25=1,30=1,116=101/1,2,3; 4//1; 5/5=2,38=5/2; - 8/5=-1,6=4,9=120000,10=1/1,4; + 8/5=-1,6=4,9=120000,10=3/1,4; 9/5=7,14=2/13; 6/7=3/1; 99/5=1,9=1/99; - Leave Link 1 at Sat Mar 30 15:07:52 2019, MaxMem= 0 cpu: 0.0 + Leave Link 1 at Mon Apr 1 12:13:25 2019, MaxMem= 0 cpu: 0.0 (Enter /share/apps/gaussian/g09d01/nehalem/g09/l101.exe) -- G2 @@ -116,7 +116,7 @@ NQMom= 10.4000000 10.4000000 NMagM= 2.2175200 2.2175200 AtZNuc= 11.0000000 11.0000000 - Leave Link 101 at Sat Mar 30 15:07:52 2019, MaxMem= 33554432 cpu: 0.0 + Leave Link 101 at Mon Apr 1 12:13:25 2019, MaxMem= 33554432 cpu: 0.1 (Enter /share/apps/gaussian/g09d01/nehalem/g09/l202.exe) Input orientation: --------------------------------------------------------------------- @@ -141,7 +141,7 @@ 2 11 0 0.000000 0.000000 -1.507479 --------------------------------------------------------------------- Rotational constants (GHZ): 0.0000000 4.8367086 4.8367086 - Leave Link 202 at Sat Mar 30 15:07:52 2019, MaxMem= 33554432 cpu: 0.0 + Leave Link 202 at Mon Apr 1 12:13:25 2019, MaxMem= 33554432 cpu: 0.0 (Enter /share/apps/gaussian/g09d01/nehalem/g09/l301.exe) Standard basis: CC-pVQZ (5D, 7F) Ernie: Thresh= 0.10000D-02 Tol= 0.10000D-05 Strict=F. @@ -325,7 +325,7 @@ Integral buffers will be 131072 words long. Raffenetti 2 integral format. Two-electron integral symmetry is turned on. - Leave Link 301 at Sat Mar 30 15:07:52 2019, MaxMem= 33554432 cpu: 0.0 + Leave Link 301 at Mon Apr 1 12:13:25 2019, MaxMem= 33554432 cpu: 0.0 (Enter /share/apps/gaussian/g09d01/nehalem/g09/l302.exe) NPDir=0 NMtPBC= 1 NCelOv= 1 NCel= 1 NClECP= 1 NCelD= 1 NCelK= 1 NCelE2= 1 NClLst= 1 CellRange= 0.0. @@ -333,10 +333,10 @@ One-electron integral symmetry used in STVInt NBasis= 118 RedAO= T EigKep= 1.92D-04 NBF= 24 7 14 14 7 24 14 14 NBsUse= 118 1.00D-06 EigRej= -1.00D+00 NBFU= 24 7 14 14 7 24 14 14 - Leave Link 302 at Sat Mar 30 15:07:52 2019, MaxMem= 33554432 cpu: 0.1 + Leave Link 302 at Mon Apr 1 12:13:26 2019, MaxMem= 33554432 cpu: 0.2 (Enter /share/apps/gaussian/g09d01/nehalem/g09/l303.exe) DipDrv: MaxL=1. - Leave Link 303 at Sat Mar 30 15:07:52 2019, MaxMem= 33554432 cpu: 0.0 + Leave Link 303 at Mon Apr 1 12:13:26 2019, MaxMem= 33554432 cpu: 0.0 (Enter /share/apps/gaussian/g09d01/nehalem/g09/l401.exe) ExpMin= 1.34D-02 ExpMax= 1.22D+06 ExpMxC= 5.43D+02 IAcc=3 IRadAn= 5 AccDes= 0.00D+00 Harris functional with IExCor= 205 and IRadAn= 5 diagonalized for initial guess. @@ -367,7 +367,7 @@ (DLTU) (SGG) (SGU) (DLTU) (DLTU) (PIU) (PIU) (PIG) (PIG) (SGG) (SGU) (PIG) (PIG) (SGG) (SGU) The electronic state of the initial guess is 1-SGG. - Leave Link 401 at Sat Mar 30 15:07:53 2019, MaxMem= 33554432 cpu: 0.2 + Leave Link 401 at Mon Apr 1 12:13:26 2019, MaxMem= 33554432 cpu: 0.3 (Enter /share/apps/gaussian/g09d01/nehalem/g09/l502.exe) Restricted open shell SCF: Using DIIS extrapolation, IDIIS= 1040. @@ -503,7 +503,7 @@ KE= 3.237479996651D+02 PE=-8.222570061205D+02 EE= 1.535554747141D+02 Annihilation of the first spin contaminant: S**2 before annihilation 0.0000, after 0.0000 - Leave Link 502 at Sat Mar 30 15:07:56 2019, MaxMem= 33554432 cpu: 3.0 + Leave Link 502 at Mon Apr 1 12:13:30 2019, MaxMem= 33554432 cpu: 3.2 (Enter /share/apps/gaussian/g09d01/nehalem/g09/l801.exe) Windowed orbitals will be sorted by symmetry type. GenMOA: NOpAll= 8 NOp2=8 NOpUse= 8 JSym2X=1 @@ -518,62 +518,62 @@ ExpMin= 1.34D-02 ExpMax= 1.22D+06 ExpMxC= 5.43D+02 IAcc=3 IRadAn= 5 AccDes= 0.00D+00 HarFok: IExCor= 205 AccDes= 0.00D+00 IRadAn= 5 IDoV=-2 UseB2=F ITyADJ=14 ICtDFT= 12500011 ScaDFX= 1.000000 1.000000 1.000000 1.000000 - Largest valence mixing into a core orbital is 2.39D-04 - Largest core mixing into a valence orbital is 7.88D-05 - Largest valence mixing into a core orbital is 2.39D-04 - Largest core mixing into a valence orbital is 7.88D-05 - Range of M.O.s used for correlation: 11 118 - NBasis= 118 NAE= 11 NBE= 11 NFC= 10 NFV= 0 - NROrb= 108 NOA= 1 NOB= 1 NVA= 107 NVB= 107 + Largest valence mixing into a core orbital is 2.84D-05 + Largest core mixing into a valence orbital is 1.38D-05 + Largest valence mixing into a core orbital is 2.84D-05 + Largest core mixing into a valence orbital is 1.38D-05 + Range of M.O.s used for correlation: 3 118 + NBasis= 118 NAE= 11 NBE= 11 NFC= 2 NFV= 0 + NROrb= 116 NOA= 9 NOB= 9 NVA= 107 NVB= 107 **** Warning!!: The largest alpha MO coefficient is 0.23725999D+02 **** Warning!!: The largest beta MO coefficient is 0.23725999D+02 - Singles contribution to E2= -0.1363659726D-15 - Leave Link 801 at Sat Mar 30 15:07:58 2019, MaxMem= 33554432 cpu: 2.1 + Singles contribution to E2= -0.1541949552D-15 + Leave Link 801 at Mon Apr 1 12:13:32 2019, MaxMem= 33554432 cpu: 2.2 (Enter /share/apps/gaussian/g09d01/nehalem/g09/l804.exe) Open-shell transformation, MDV= 33554432 ITran=4 ISComp=2. Semi-Direct transformation. - ModeAB= 4 MOrb= 1 LenV= 32899201 - LASXX= 79572 LTotXX= 79572 LenRXX= 162186 - LTotAB= 82614 MaxLAS= 654372 LenRXY= 0 - NonZer= 241758 LenScr= 720896 LnRSAI= 654372 - LnScr1= 1441792 LExtra= 0 Total= 2979246 + ModeAB= 4 MOrb= 9 LenV= 32869529 + LASXX= 844618 LTotXX= 844618 LenRXX= 1721496 + LTotAB= 876878 MaxLAS= 6325596 LenRXY= 0 + NonZer= 2566114 LenScr= 4325376 LnRSAI= 6325596 + LnScr1= 10092544 LExtra= 0 Total= 22465012 MaxDsk= -1 SrtSym= T ITran= 4 DoSDTr: NPSUse= 1 - JobTyp=1 Pass 1: I= 1 to 1. + JobTyp=1 Pass 1: I= 1 to 9. (rs|ai) integrals will be sorted in core. Complete sort for first half transformation. First half transformation complete. Complete sort for second half transformation. Second half transformation complete. - ModeAB= 4 MOrb= 1 LenV= 32899201 - LASXX= 79572 LTotXX= 79572 LenRXX= 160530 - LTotAB= 80958 MaxLAS= 654372 LenRXY= 0 - NonZer= 240102 LenScr= 720896 LnRSAI= 654372 - LnScr1= 1441792 LExtra= 0 Total= 2977590 + ModeAB= 4 MOrb= 9 LenV= 32869529 + LASXX= 844618 LTotXX= 844618 LenRXX= 1588696 + LTotAB= 744078 MaxLAS= 6325596 LenRXY= 0 + NonZer= 2433314 LenScr= 4325376 LnRSAI= 6325596 + LnScr1= 10092544 LExtra= 0 Total= 22332212 MaxDsk= -1 SrtSym= T ITran= 4 DoSDTr: NPSUse= 1 - JobTyp=2 Pass 1: I= 1 to 1. + JobTyp=2 Pass 1: I= 1 to 9. (rs|ai) integrals will be sorted in core. Complete sort for first half transformation. First half transformation complete. Complete sort for second half transformation. Second half transformation complete. Spin components of T(2) and E(2): - alpha-alpha T2 = 0.0000000000D+00 E2= 0.0000000000D+00 - alpha-beta T2 = 0.3401125020D-01 E2= -0.1985060912D-01 - beta-beta T2 = 0.0000000000D+00 E2= 0.0000000000D+00 - ANorm= 0.1016863437D+01 - E2 = -0.1985060912D-01 EUMP2 = -0.32373579504585D+03 + alpha-alpha T2 = 0.1028504344D-02 E2= -0.4527483160D-02 + alpha-beta T2 = 0.3880206724D-01 E2= -0.4557051157D-01 + beta-beta T2 = 0.1028504344D-02 E2= -0.4527483160D-02 + ANorm= 0.1020225012D+01 + E2 = -0.5462547789D-01 EUMP2 = -0.32377056991462D+03 (S**2,0)= 0.00000D+00 (S**2,1)= 0.00000D+00 - E(PUHF)= -0.32371594444D+03 E(PMP2)= -0.32373579505D+03 - Leave Link 804 at Sat Mar 30 15:08:06 2019, MaxMem= 33554432 cpu: 8.2 + E(PUHF)= -0.32371594444D+03 E(PMP2)= -0.32377056991D+03 + Leave Link 804 at Mon Apr 1 12:13:46 2019, MaxMem= 33554432 cpu: 13.2 (Enter /share/apps/gaussian/g09d01/nehalem/g09/l913.exe) CIDS: MDV= 33554432. - Frozen-core window: NFC= 10 NFV= 0. + Frozen-core window: NFC= 2 NFV= 0. IFCWin=0 IBDFC=1 NFBD= 0 0 NFCmp= 0 0 NFFFC= 0 0 Using original routines for 1st iteration, S=T. Using DD4UQ or CC4UQ for 2nd and later iterations. @@ -588,155 +588,115 @@ Iterations= 50 Convergence= 0.100D-06 Iteration Nr. 1 ********************** - DD1Dir will call FoFMem 1 times, MxPair= 2 - NAB= 1 NAA= 0 NBB= 0. - DD1Dir will call FoFMem 1 times, MxPair= 2 - NAB= 1 NAA= 0 NBB= 0. - MP4(R+Q)= 0.43884995D-02 + DD1Dir will call FoFMem 1 times, MxPair= 234 + NAB= 81 NAA= 36 NBB= 36. + DD1Dir will call FoFMem 1 times, MxPair= 234 + NAB= 81 NAA= 36 NBB= 36. + MP4(R+Q)= 0.54052112D-02 Maximum subspace dimension= 5 - Norm of the A-vectors is 7.5880349D-03 conv= 1.00D-05. - RLE energy= -0.0191976723 - E3= -0.37133555D-02 EROMP3= -0.32373950840D+03 - E4(SDQ)= -0.21234154D-02 ROMP4(SDQ)= -0.32374163182D+03 + Norm of the A-vectors is 1.2772787D-02 conv= 1.00D-05. + RLE energy= -0.0537245945 + E3= -0.44892214D-02 EROMP3= -0.32377505914D+03 + E4(SDQ)= -0.29522642D-02 ROMP4(SDQ)= -0.32377801140D+03 VARIATIONAL ENERGIES WITH THE FIRST-ORDER WAVEFUNCTION: - DE(Corr)= -0.19175465E-01 E(Corr)= -323.73511990 - NORM(A)= 0.10156179D+01 + DE(Corr)= -0.53709488E-01 E(Corr)= -323.76965392 + NORM(A)= 0.10185829D+01 Iteration Nr. 2 ********************** - DD1Dir will call FoFMem 1 times, MxPair= 2 - NAB= 1 NAA= 0 NBB= 0. - Norm of the A-vectors is 6.7903864D-02 conv= 1.00D-05. - RLE energy= -0.0194071523 - DE(Corr)= -0.22754797E-01 E(CORR)= -323.73869923 Delta=-3.58D-03 - NORM(A)= 0.10161175D+01 + DD1Dir will call FoFMem 1 times, MxPair= 234 + NAB= 81 NAA= 36 NBB= 36. + Norm of the A-vectors is 1.5757331D-01 conv= 1.00D-05. + RLE energy= -0.0544053149 + DE(Corr)= -0.57996151E-01 E(CORR)= -323.77394059 Delta=-4.29D-03 + NORM(A)= 0.10199081D+01 Iteration Nr. 3 ********************** - DD1Dir will call FoFMem 1 times, MxPair= 2 - NAB= 1 NAA= 0 NBB= 0. - Norm of the A-vectors is 6.6017986D-02 conv= 1.00D-05. - RLE energy= -0.0216715304 - DE(Corr)= -0.22899343E-01 E(CORR)= -323.73884378 Delta=-1.45D-04 - NORM(A)= 0.10219533D+01 + DD1Dir will call FoFMem 1 times, MxPair= 234 + NAB= 81 NAA= 36 NBB= 36. + Norm of the A-vectors is 1.3847674D-01 conv= 1.00D-05. + RLE energy= -0.0567541944 + DE(Corr)= -0.58523914E-01 E(CORR)= -323.77446835 Delta=-5.28D-04 + NORM(A)= 0.10244303D+01 Iteration Nr. 4 ********************** - DD1Dir will call FoFMem 1 times, MxPair= 2 - NAB= 1 NAA= 0 NBB= 0. - Norm of the A-vectors is 4.8512481D-02 conv= 1.00D-05. - RLE energy= 0.5274153290 - DE(Corr)= -0.24130071E-01 E(CORR)= -323.74007451 Delta=-1.23D-03 - NORM(A)= 0.88172211D+01 + DD1Dir will call FoFMem 1 times, MxPair= 234 + NAB= 81 NAA= 36 NBB= 36. + Norm of the A-vectors is 1.0806968D-01 conv= 1.00D-05. + RLE energy= -0.0651329656 + DE(Corr)= -0.59639289E-01 E(CORR)= -323.77558373 Delta=-1.12D-03 + NORM(A)= 0.10523749D+01 Iteration Nr. 5 ********************** - DD1Dir will call FoFMem 1 times, MxPair= 2 - NAB= 1 NAA= 0 NBB= 0. - Norm of the A-vectors is 7.4726707D+00 conv= 1.00D-05. - RLE energy= -0.0137310363 - DE(Corr)= 1.2104798 E(CORR)= -322.50546459 Delta= 1.23D+00 - NORM(A)= 0.10080247D+01 + DD1Dir will call FoFMem 1 times, MxPair= 234 + NAB= 81 NAA= 36 NBB= 36. + Norm of the A-vectors is 2.9368612D-02 conv= 1.00D-05. + RLE energy= -0.0634626252 + DE(Corr)= -0.64336969E-01 E(CORR)= -323.78028141 Delta=-4.70D-03 + NORM(A)= 0.10460210D+01 Iteration Nr. 6 ********************** - DD1Dir will call FoFMem 1 times, MxPair= 2 - NAB= 1 NAA= 0 NBB= 0. - Norm of the A-vectors is 1.1238829D-01 conv= 1.00D-05. - RLE energy= -0.0235876785 - DE(Corr)= -0.19239089E-01 E(CORR)= -323.73518353 Delta=-1.23D+00 - NORM(A)= 0.10284678D+01 + DD1Dir will call FoFMem 1 times, MxPair= 234 + NAB= 81 NAA= 36 NBB= 36. + Norm of the A-vectors is 2.4030426D-03 conv= 1.00D-05. + RLE energy= -0.0634241779 + DE(Corr)= -0.63452584E-01 E(CORR)= -323.77939702 Delta= 8.84D-04 + NORM(A)= 0.10459648D+01 Iteration Nr. 7 ********************** - DD1Dir will call FoFMem 1 times, MxPair= 2 - NAB= 1 NAA= 0 NBB= 0. - Norm of the A-vectors is 3.5616956D-02 conv= 1.00D-05. - RLE energy= -0.0225957877 - DE(Corr)= -0.25468394E-01 E(CORR)= -323.74141283 Delta=-6.23D-03 - NORM(A)= 0.10253695D+01 + DD1Dir will call FoFMem 1 times, MxPair= 234 + NAB= 81 NAA= 36 NBB= 36. + Norm of the A-vectors is 7.2581575D-04 conv= 1.00D-05. + RLE energy= -0.0634445103 + DE(Corr)= -0.63429775E-01 E(CORR)= -323.77937421 Delta= 2.28D-05 + NORM(A)= 0.10460949D+01 Iteration Nr. 8 ********************** - DD1Dir will call FoFMem 1 times, MxPair= 2 - NAB= 1 NAA= 0 NBB= 0. - Norm of the A-vectors is 4.4157577D-02 conv= 1.00D-05. - RLE energy= -0.0183305183 - DE(Corr)= -0.24969117E-01 E(CORR)= -323.74091355 Delta= 4.99D-04 - NORM(A)= 0.10169406D+01 + DD1Dir will call FoFMem 1 times, MxPair= 234 + NAB= 81 NAA= 36 NBB= 36. + Norm of the A-vectors is 3.0343948D-04 conv= 1.00D-05. + RLE energy= -0.0634349134 + DE(Corr)= -0.63441626E-01 E(CORR)= -323.77938606 Delta=-1.19D-05 + NORM(A)= 0.10460542D+01 Iteration Nr. 9 ********************** - DD1Dir will call FoFMem 1 times, MxPair= 2 - NAB= 1 NAA= 0 NBB= 0. - Norm of the A-vectors is 8.1339530D-02 conv= 1.00D-05. - RLE energy= -0.0395291831 - DE(Corr)= -0.22701643E-01 E(CORR)= -323.73864608 Delta= 2.27D-03 - NORM(A)= 0.11366219D+01 + DD1Dir will call FoFMem 1 times, MxPair= 234 + NAB= 81 NAA= 36 NBB= 36. + Norm of the A-vectors is 1.0969098D-04 conv= 1.00D-05. + RLE energy= -0.0634370384 + DE(Corr)= -0.63436482E-01 E(CORR)= -323.77938092 Delta= 5.14D-06 + NORM(A)= 0.10460603D+01 Iteration Nr. 10 ********************** - DD1Dir will call FoFMem 1 times, MxPair= 2 - NAB= 1 NAA= 0 NBB= 0. - Norm of the A-vectors is 1.2365186D-01 conv= 1.00D-05. - RLE energy= -0.0269817872 - DE(Corr)= -0.33166721E-01 E(CORR)= -323.74911116 Delta=-1.05D-02 - NORM(A)= 0.10423981D+01 + DD1Dir will call FoFMem 1 times, MxPair= 234 + NAB= 81 NAA= 36 NBB= 36. + Norm of the A-vectors is 3.6902429D-05 conv= 1.00D-05. + RLE energy= -0.0634368383 + DE(Corr)= -0.63436982E-01 E(CORR)= -323.77938142 Delta=-5.01D-07 + NORM(A)= 0.10460594D+01 Iteration Nr. 11 ********************** - DD1Dir will call FoFMem 1 times, MxPair= 2 - NAB= 1 NAA= 0 NBB= 0. - Norm of the A-vectors is 4.8735812D-03 conv= 1.00D-05. - RLE energy= -0.0272582911 - DE(Corr)= -0.27187932E-01 E(CORR)= -323.74313237 Delta= 5.98D-03 - NORM(A)= 0.10438715D+01 + DD1Dir will call FoFMem 1 times, MxPair= 234 + NAB= 81 NAA= 36 NBB= 36. + Norm of the A-vectors is 1.4952661D-05 conv= 1.00D-05. + RLE energy= -0.0634369208 + DE(Corr)= -0.63436854E-01 E(CORR)= -323.77938129 Delta= 1.28D-07 + NORM(A)= 0.10460599D+01 Iteration Nr. 12 ********************** - DD1Dir will call FoFMem 1 times, MxPair= 2 - NAB= 1 NAA= 0 NBB= 0. - Norm of the A-vectors is 1.8116780D-03 conv= 1.00D-05. - RLE energy= -0.0274179164 - DE(Corr)= -0.27336889E-01 E(CORR)= -323.74328133 Delta=-1.49D-04 - NORM(A)= 0.10447323D+01 - Iteration Nr. 13 - ********************** - DD1Dir will call FoFMem 1 times, MxPair= 2 - NAB= 1 NAA= 0 NBB= 0. - Norm of the A-vectors is 7.1205872D-05 conv= 1.00D-05. - RLE energy= -0.0274344242 - DE(Corr)= -0.27420114E-01 E(CORR)= -323.74336455 Delta=-8.32D-05 - NORM(A)= 0.10448189D+01 - Iteration Nr. 14 - ********************** - DD1Dir will call FoFMem 1 times, MxPair= 2 - NAB= 1 NAA= 0 NBB= 0. - Norm of the A-vectors is 1.4943057D-04 conv= 1.00D-05. - RLE energy= -0.0274211660 - DE(Corr)= -0.27428216E-01 E(CORR)= -323.74337265 Delta=-8.10D-06 - NORM(A)= 0.10447456D+01 - Iteration Nr. 15 - ********************** - DD1Dir will call FoFMem 1 times, MxPair= 2 - NAB= 1 NAA= 0 NBB= 0. - Norm of the A-vectors is 1.2040241D-05 conv= 1.00D-05. - RLE energy= -0.0274215001 - DE(Corr)= -0.27421233E-01 E(CORR)= -323.74336567 Delta= 6.98D-06 - NORM(A)= 0.10447470D+01 - Iteration Nr. 16 - ********************** - DD1Dir will call FoFMem 1 times, MxPair= 2 - NAB= 1 NAA= 0 NBB= 0. - Norm of the A-vectors is 3.4337654D-06 conv= 1.00D-05. - RLE energy= -0.0274214303 - DE(Corr)= -0.27421439E-01 E(CORR)= -323.74336588 Delta=-2.06D-07 - NORM(A)= 0.10447466D+01 - Iteration Nr. 17 - ********************** - DD1Dir will call FoFMem 1 times, MxPair= 2 - NAB= 1 NAA= 0 NBB= 0. - Norm of the A-vectors is 1.2995941D-06 conv= 1.00D-05. - RLE energy= -0.0274214089 - DE(Corr)= -0.27421417E-01 E(CORR)= -323.74336585 Delta= 2.20D-08 - NORM(A)= 0.10447465D+01 - CI/CC converged in 17 iterations to DelEn= 2.20D-08 Conv= 1.00D-07 ErrA1= 1.30D-06 Conv= 1.00D-05 - Largest amplitude= 9.11D-02 - Time for triples= 0.00 seconds. - T4(CCSD)= 0.00000000D+00 - T5(CCSD)= 0.00000000D+00 - CCSD(T)= -0.32374336585D+03 + DD1Dir will call FoFMem 1 times, MxPair= 234 + NAB= 81 NAA= 36 NBB= 36. + Norm of the A-vectors is 5.1542774D-06 conv= 1.00D-05. + RLE energy= -0.0634369469 + DE(Corr)= -0.63436919E-01 E(CORR)= -323.77938136 Delta=-6.44D-08 + NORM(A)= 0.10460601D+01 + CI/CC converged in 12 iterations to DelEn=-6.44D-08 Conv= 1.00D-07 ErrA1= 5.15D-06 Conv= 1.00D-05 + Largest amplitude= 8.86D-02 + Time for triples= 681.31 seconds. + T4(CCSD)= -0.14211246D-02 + T5(CCSD)= 0.72680615D-04 + CCSD(T)= -0.32378072980D+03 Discarding MO integrals. - Leave Link 913 at Sat Mar 30 15:08:18 2019, MaxMem= 33554432 cpu: 5.2 + Leave Link 913 at Mon Apr 1 13:07:30 2019, MaxMem= 33554432 cpu: 715.7 (Enter /share/apps/gaussian/g09d01/nehalem/g09/l601.exe) Copying SCF densities to generalized density rwf, IOpCl= 0 IROHF=1. @@ -8447,19 +8407,18 @@ --------------------------------------------------------------------------------- No NMR shielding tensors so no spin-rotation constants. - Leave Link 601 at Sat Mar 30 15:08:18 2019, MaxMem= 33554432 cpu: 0.2 + Leave Link 601 at Mon Apr 1 13:07:31 2019, MaxMem= 33554432 cpu: 0.3 (Enter /share/apps/gaussian/g09d01/nehalem/g09/l9999.exe) - 1\1\GINC-COMPUTE-7-0\SP\ROCCSD(T)-FC\CC-pVQZ\Na2\LOOS\30-Mar-2019\0\\# - p ROCCSD(T) cc-pVQZ pop=full gfprint\\G2\\0,1\Na\Na,1,3.0149584\\Versi - on=ES64L-G09RevD.01\State=1-SGG\HF=-323.7159444\MP2=-323.735795\MP3=-3 - 23.7395084\PUHF=-323.7159444\PMP2-0=-323.735795\MP4SDQ=-323.7416318\CC - SD=-323.7433659\CCSD(T)=-323.7433659\RMSD=8.165e-09\PG=D*H [C*(Na1.Na1 - )]\\@ + 1\1\GINC-COMPUTE-7-0\SP\ROCCSD(T)-FC1\CC-pVQZ\Na2\LOOS\01-Apr-2019\0\\ + #p ROCCSD(T,FreezeInnerNobleGasCore) cc-pVQZ pop=full gfprint\\G2\\0,1 + \Na\Na,1,3.0149584\\Version=ES64L-G09RevD.01\State=1-SGG\HF=-323.71594 + 44\MP2=-323.7705699\MP3=-323.7750591\PUHF=-323.7159444\PMP2-0=-323.770 + 5699\MP4SDQ=-323.7780114\CCSD=-323.7793814\CCSD(T)=-323.7807298\RMSD=8 + .165e-09\PG=D*H [C*(Na1.Na1)]\\@ - There ain't no surer way to find out whether you - like people or hate them than to travel with them. - -- Mark Twain - Job cpu time: 0 days 0 hours 0 minutes 19.4 seconds. - File lengths (MBytes): RWF= 70 Int= 0 D2E= 0 Chk= 2 Scr= 1 - Normal termination of Gaussian 09 at Sat Mar 30 15:08:18 2019. + The lyf so short, the craft so long to lerne. + -- Chaucer + Job cpu time: 0 days 0 hours 12 minutes 15.6 seconds. + File lengths (MBytes): RWF= 227 Int= 0 D2E= 0 Chk= 2 Scr= 1 + Normal termination of Gaussian 09 at Mon Apr 1 13:07:31 2019. diff --git a/G09/Mixed_core/Molecules/vqz/copy.sh b/G09/Mixed_core/Molecules/vqz/copy.sh index 027ded3..7fd2669 100755 --- a/G09/Mixed_core/Molecules/vqz/copy.sh +++ b/G09/Mixed_core/Molecules/vqz/copy.sh @@ -1,8 +1,9 @@ -for i in NaCl LiF +for i in Li2 BeH LiH Na2 do - for j in vdz vtz vqz + #for j in vdz vtz vqz + for j in vqz do - cp Mixed_core/${j}/${i}.out ../${j}/ + cp ../../../Small_core/Molecules/${j}/${i}.out ../${j}/ done done diff --git a/Manuscript/G2-srDFT.pdf b/Manuscript/G2-srDFT.pdf index be4398cdc4d1026ccb293b2e39c7476ef17060f4..b3d02c3e2f51576a7fe86d1bae50b77c44ffe44b 100644 GIT binary patch delta 19779 zcmV)BK*PVyjUE4v9k47i0W*;i69YFiHj|NLDStg{bK^#G-}hJW?!IsiPz z>ylh5u{V`VY=597u7ooQ>PTvJn_r*qc>puu0H7plSC0LOAUMFxbocbTn`Cx9BeQQ_ zeg5XvCtp(`X1215so9&iGbwc?6q{)&E2H4n)$IE@fBoZ|e||;6Q!`RoW8jj^=uC>4 zW`FoNpS{_@&>t8TtJe!A*&IIf{BL*5L%nVA1u)t_Ly4)3?LSB)3De-1-E|oJEll|R zoa%fDZaj5wO;L|GbPN{yU72cBU(1tQ(;#&FG9*iUwqVKCoKK>1>gViwNjb z5CEW^@zWBhIIT210nn@%vwq+(h|Eg?5$ubK6F`)Vb0AtUW9R?Wy*GJ++KMwiqkjVK z1gPFRI{)g;tG~T!kBi|Ifj4UcOHz7vxq0>dk7RZQKO*%DVrKVl=mu8GB!{c@?7LV0 z5sr6GrtmH*6}*d71U5A|r-p5CIa5|v`gl%myE!TBTiE`jDg3Fnm8N#4jH-lSeo_{3 zk=9&^ze!Yz@q}Dyi5Ou~#Ku8Y=>P%eRf48u&q&7YBCCPJDCtaT*dhzO?z7jF&4J%h zE$8>I4e~sy)%+q5tykpz!FEW;oYXuTkid0GkwiwpOB$+SQ`>nr%Y26FlO7AY?@BX8m z5Wnn}&Gl-rU)?Tu%gCO|9KUT=u=McYzQlnP2PgsoF;N*Ia}ETk=mf(M_*u1rM$6-uy5DmlWgOl2*wD{d-&Jm&-1l@UMy zz(|;JWO5R_5;?odsK`1E1p$3TL>XIK02hMoFAR)S6MT)H#MqQYlBE2_ur^X@le0FO zyp7u#fLsBL5gMQzz(x)5dE9yc@RzTdg?BF~A!Pp9?X71Pdj9u%Dq-FK^(*rTrqe zmVri*>3~Wq1)L(zJPZ6&6Y;API*3$_j2(`&l&rGC_$U52HH|?Io6izp06hQhaI@X* z{g<2NyU2j$+{5t}xqt3EMd-bIvAR(JM+CT(#X(1Zgzo!&X;#pzSWEuRz>(CZ zRcL=4KxGR%sz6gll1TyhQvloDK3FKNblf}ic`00@yuZZzbbo%iUhc!O*uM3j-+lY# z*KfRl1|5nexB>%8&_HWQUVmO}u8Ihz5+nix-i(ZEoK(yMNaabK8|~2v&MksW%+(H-@pj8q9C3m{ z`y-&wl2(z(W@AWntu)!xO;H$1FjN3*8kAu8x$NEn3XQd8 zSa+}CKn11%>+Ut&nM}u^&b26@a&GXw!TP!x(((YY?hb|hy(jMiRS#aXNXGrOlkhhc ze>X*38(4me!G9>8P+m!(app9n@ku#YU<@a1hYqMwXa&O(MLh8<&+tl+0TYp+KF(DS zF7lXTKOD-nXAG18PvGdvfGyZo%Ob?{HXw7fGLq*WhHG@gwYxvk>Nf*yqp#IJg#pIBKW#d7bT^PqF0&C zDCpx>!!1g2&^MNRxkYQp(hB0&^Wt7eM zi1^~x<24swBMuq^*Z`aeLV&_d8w`G;hIB?W>q#_AJB;Tc!Gk7%0;Ma1my0;7AX@NZ zy=|@+E?|M*jSrS67c9|p5Ieo^{-HDi(-TYFXfU>eN65SFVCZ36j3-dse-j62lgdhY zF5UITVCKjHqk#NCsn%huMInT$oYcS#b#yM!LR5i(G*EiRl>&w@$D#_r$(>(&IH(Zc z1saO&g+LR3v&IK7(88i`|9{Zq%;2M}mnwnk1*?P%hySW1$U$#VvJiAiFy*3LQYXMT z4(u}{X340*snqmTI!f|le_%yd6qs%@>b^u7y5b<35@qPDM0l)uu}t&ol!GIwo@h;f zEx{&d>#d7RHpcdo`?Q)@C=iW!AXB#g4h%(e=R~P+uU%d2xBTBvxJ;fL%8>$ z#F0g{4b)kgh1e48WE;iwhZZ@kth$-f^81_AS`U!JY4P4ET=zsge=|Hq^~KOsk3n-^ zqsI_!GHk?YEc?f*&z{f$mZ}`0RXM;+K)*UVXJ}^}-Pebk5*bJedk}7{jUwk)$`i~K zPf)XD*eFksACD(wX$CpfSg_gN`l|a$E7@c3JNhSF2`UAKv>rlcqz9ejh_+R#Zf> z3G1#C*v+h!O){{d6r<&32-fM zHaW~wI_4QSY0k9C6EJEAl(0T;6DfBfMHvXnmU5{i3y_(jXTKdyZ90G6xl#EPdH4ZU z+ZRnbC!7yEa6bBp`XW!)bJ;+>d^M)=rbwzgVp>m<669}O)h?GfWZu1X`C^-18s7c@ za%qV_g%HZOf8Sv$Qb^I;IXB~a+2xU?(oNb*qFXr45Ta@Rm@Clm2~spp2?Lgwhj?pw z$m%`l+)en_yEj8>F%1q9*-76XJXr3=5rjb4-z!8x7be;HStN0)7Vd4sVSJKrpSUhVxQ z>;nf4g>>qg(;h_Xn3#3bT)BMhvt#X*GOehLrq&o0pizS_R@(<)C%JRuyC04vwLa zs(4fgfI!7KZFB9?YQZ4<;>y2}zj*!;S7hoJG_+C)^ zq*4ih1_)A6n+bbBPN|;EB1{76y(~fv3N-Dqe+U;P`~yuThV5U9DAO+Eh*v^9I6I>+ z>J&dm>fdBut7-FE(Qrq9mL-*mqGi4O35fu>|GJB?i@?o2Y8U00TT*wR70uuJm&vt& zDQXDIpQ)i(9x1ucu;o2b!vm9zqDpl1O*IfitaI7a*EfN&U9G&q)m(XT60=-e zfB6!vvr@4^TkF8!*yU3^Q7Pw-OCdc5+ zbDB8UXBFp1she58jQ6ZyD3{3-%7_$^`LX~V!1K{t;){}qmM$He27?0CuW}8Xf0EjX z3p1`swD@I~RRV-dhTnRBq~O9d+&H_&pLs+rm-D;*cuE3uwQ^LD%Kh?76Yt_A@opx? zyEt=Z{{WwAD|NJ6RfaJa*Ks&QHZ#H|P`Y8kTK4q03+Nxoh!g3i(Zo)d*yp+&D`+!| z2yOu`bpF^x*8@4)^`PZhKnAGzf21$=QkKUGEiU4zP18uQU4oS??qE>$;0y>>oXT^= zNjyiK$a8d`=k)V9gjKE%NJaK^9_#B7^|Mc4*Y-4aeHynGCv$66;?}teB~y*)2g2z< zh93iKB!TQBkZ{Cj=}0lwdHF2n#_3-M5;M44(2iJ&!?SSs50kwu-~g#?e@>>-xK`>f zFc{ZLeURR57d%&BHC_YX>F5&cZdV1EVEGgelnDc*{@=Rwwuq zG%zDA9<6T-xPtVFy5rgC_?h|mknX5Q)C`^$4pyatZH9TSY6dU!{dZh2AXaVvtE(RW ztV|+p$?4;3#|Kpgj)GQWf7=z&m1PeSvQWV}Dr6v0Arb`xzxlfocf<+oi0k=cUH=?@ z&HNPjeLT7=xvezzV{C3Brd7U{qTB&X(?r$We&RTx>ILsKzg@q2f{N`uC>ocOM(-+Jv#!x0eJ!G!4JQu$s60(0dOqV2oHiV;?< z>+&a-DEtIQ5MEf{sw@7o6e?P)p0&E4MsH?NA=6$WhC5Y6v^~#EDs# z5HEeZAWcs%Uf}{b1Mz)^9D#49P8!Q{lArsevC)%8Cz?`#4C4xmL4`KQsa9f$*rvlp z7iS{;!GJ=GSc|E|JV$Jd(zn4A7fD}aL3F&mzuE2%ZS@T{e?-~tj)IaLj70(fRUFVt zXbcg!Y=od1YT9)~(QF=4=ll=1%iW>A#K;wXb1roqUFvWH(YD=f{N3YX9~g1H-2025 ztKr({5nWT`OCZw-KWK|f{}KG~yI7E#Dv&C%aRfGVyK@dAxU9p#Z&z=tt{YK7OC-YS z9oZuOIN;%Jf0r(?(gx@ud&+t6HFiJp>vmO6mqGuH=AY?fhUki#;3M`ry3jF$v03bI zR_hoAc)vfaHoi}dZrjYkrL*i@Y~Uu`+27X9!Uv+(O8{JJ`N3~4uC&_+NE*ayL`dI}i3X5|AL4^T9pXDF|e?fsF~^gB#ms9iH}!go|c- zh(w12-gZByKF4-N&dVG4ti!=wQ^JC?(f1x=1MRsucR5SGPi-{x{tYuy?)?k?{qwKd z)`F{Df1Krz&e3@U#w|KTz3iEE*tUoons>hi#6(4!tIjf0J(4(T%oxsN3NitpQlNF` z^NhNUxfygDr#lc98PUXtLHGXl_U%|)0ws`w4QH?1Ok~Ok0^$FkZu8wX`W ze`Sxr(G-k&o{w|76u?zjnBy?EpZFCP?Xxhf?{j7`GvI9ypcg)i{XVL9nad>tLYa!o z!2J4c?iiEvE-8(9NecIx$UopZFbCETVTHxG3*6Yfwm855b+KUTqg;HOpxFdr?V!K)8e|K8>k!9>*Ksbfk(Qk=4n~a zRV9JUvz>G~G9@Yny<6~@LYDJqJQ7?-wdHUh;f#r3_=C?h<0zKyTp}ObbG#lek|W}T z5{}DxyVEgfPr(wLJC8#&9GcOCA!ua*7}-PR0wYryanzWiA3*7RI%Ftr-cGG)elyZU#Q?nx0(7qo^WK;M);Ugl(H7CXwor%_szMc^$&xPCPpE zvk>-5Qa~QxgAjnJ)k|#wHc?Lce|(7wDX9qT*#LG%Mk-!;x9%J$e-1=QgBHUIEec5* z8KqY+3Om5q1h<852Pt)nj=IaOTqX zK1cxf!*|ns&Q!?f>85E;>86RTP$hGZN;OuM=v`deHJ2Z1qzx*XQ}(zaf6meJl;+*| zE(I&&n6S#yhT*xx5X?DD_{klHq};p+-O8E6kPdXtA%gddqGA;hCL;1-2q;)BY!hBp zNzTV0Fxm8@eOC$rb_@ZJnUc%2w>E;Z zjG-+F&vM1veDG^$i=w}Hz>4EQW5Psl-og3SIOR@}#!qjl!G(9sq9?$T_OD8^5uM3` z-@2LoLw zLL;ah_6!Gq*nYVvwqr`|m%+{1Y;%8j`sN-?q?mj<#NM#v#EMKYvZcN%Q7^*MAQk_H z29I+E}wwq;hH3|VsT-&bGAYjf-F7mt(5I_i;_aTsG zZi24m18C`Af3buD#r6Bw()d9%>N{e7ck4fWvmUs(3SOrHea7|ikhBY&#+ttn7*0$2 zc&EboC|m51;Kf;}2t14-P=E%DHfc8PZbDhWs1NNwbTK*&=693c?p};!u%}A~_S}-8 z$EPJX>0!7Snod3&ub$Ez@&B3x=ph+udp6LqeusYze``Nm=Wxb5Z z$-Ru8Kp6Hk!qAh1@w}>Gk=<}xsZ5d+PRg2`XD`$MW1uIhF!&b^g6GdX8FBV}h&$Vf zJ`Hi^f1HR*aRVv<`nobW`?MGlI~dC2us-~V?VWsTkMJHRBI5(X3lF^Te@9mmHdJX? zugzK+Bq@}GWnwXf(q)zvgLE8y11AOPc;?ut!oN%5fALwGSF1#?Sr+IIxut$3?sFnR zj208bdTcB~ZIu!Pf3T-c*TQr#*4cWd2g-0jmEf&Q0hS26 z|D-q{t$3FvJrboz;LL}+oTuoxI9A5JMFY}Q=4^wO+c`8!c<(di9ZF651}U1Tdz+Sar5=kB3%|!DeER~W{9nFtO7Pk416(_ ze?Eh**6gvV#4-*$aoi@G%Z16E#^86=Wu6XAj*&DZ&*J1q|Lja#Z%A1qq3(|rCX#nr zGZHmnVFtI5bVN>9xBBuie@{8Sh`Ike|E7V|yujStjzrGfG!W#{o%Ho*d8aXtNwB;kom^U|J-&qgsE@4#L{}- zd1J5+;r?~_H<2UVz?#Z3J&-$^MhIG#-!?0Bk*tt#Es)X!~E= z7yiGtk11AtT7pAf912Hfs5rFOrk1e0PwlEiDk3ElJl26uU|mQi%$!Un^fAWGNtE&d zV;;x<;V_OKs8$`an~ZeF%@dlDyRF zukSPD>3pZyq>vu^;jmn4cIMgH*?DFtMUosDMV_LhL`o@Al9cwO3@H~#SyCP|JSF@h zX{$*q(!ne=o{}XE4ftV=)e>*WIz-p{u@pO38mbaT zNF}`2e-KL$*N{sdV~h)IN*nk{V?A3{627wRUWpZU*bMyX*aB3!!36u6%>#u_nnmDi z$s2>K49qOx0XW&Rz6du5K0`qDd|*&+JbMcmVUpd0Rl)KbtO|)8hgH$AEgVlo5Fi*# z(r%&n^b-3Lr)?$7D4e&I4r=&-B!>g`jV)k>f5R9De|=cj&^)vhz6bb!6P0JLI#ll=unEl(5yG$!LtM+ z#xY&e~BD|t3EJxryIic)|vhSPghO%>d zf6%! zy!^U8J?gux`?k5;3@<07%_x2Sv25BYt#&6z=NH%XKJS*@tFG*+XxerZqGe4v)tkQi zSXR}9F85g%>nWh4F3NS^O|0?bn)C90e_1aFYt22ZwPjwt>ehKt;&I)UU0L_}xtlm?zF5A3sTI$j zD>-UEU9a2b({eSa>UJA0n$@~0S9})j1Kq>n{P%n+CX_mKaaSPZCn0R({D<>X%J*UX zPbjw6;ivN7^!D|e(bemw8sg>6*^lSXRI|M;e;?28RZ~|FM+CXOKkv$JQ%znPUdweCBvRUo*?;K*jv7D%?>z2VJ|<=rxW-F%`Oe-`nI$yDDi zprFM{P>vbpPV2gXq8lcCwh7K3hA<4biv#O}@FKg~yzhtk56k-3MfSRB@5*+tBX6_6 zWS805jT$Bl3`cSU-V@B&K}mQLtY?xi3`rYdt%ckvJ$2$eT+<`e5WZW2ISLf5l#rTh zH3!qERD$rRv9$NXnlB?MfB&@^j#}smAYZZ!z9qG@3%C54%y}RV!i&eg?`sUA7({$f0+b&-fVzcgE^;K7j`}XB$czbd5 z%Rm0fr6VdR5#auMQ&qQL5i?PXq4u*u1QZHWvM58ncW|ek6|+*Hf2f)jQ%(doYt$zA z%uH)T1tqNWT$Krp;%!PcN}*A38H2rqIpgpliqvTzBB&&4PA@$a`e(@+2fw2Zo0JV! z7?h>TPD?N-~<=1b(r* z2|Aa@1je0H^t%ZZe}u=0$C|jCpy^4w37Q4@I6>*0c1F?t?1a!+JX+DXRw?Wzfl}N} z5~Wcw_LIGBo8k&p%ndT>`2}Uy<){8>l6#V3mC{dhELZSVLQfS z=3VqMLUQL~&2SmeO$ej8pS76JCOdb$w2OXH3|+<+h~wMof1U}%9uJ3sL2Wo8Igt<> z#!C?qmfm^H%^Z>=dIW4EBwAX_9Edi{0ue4P@q4TE?BN{IDHeH!DanK&;0qf)QikC! z7%6F<6ZJllbB;RoEZRAc!_7b>Kv+c))2y%;fj^!&Bu|{;q_Af~8+ut`b(~X-?2gRh zfW~M$Ee!W)e{R^BDV>{9tZ>Tc+)cDY1vDuPmVsgr<_rLD1%E7v2%MRH8HidP;9VtB z3^T%y^o7XOX-;zhd3yfvfDrP}Q*5y&S@k?&V+;zB8Ceiw0GVh6kD3F_Ctx`hCE?V} zNgxIW#awVW|5ai zSoqV{IBw@g9+&Bm$2|fZrA%`CwT%-Kf8yk7JI>~?aKsJ+COs0lLOE%NPnw+3_ye*R zCfWZd>vCfE{6P=be(1a6$a;>wje)*gXG>!{dD6 zx#wwi4xBU>33D~E({{m@MbY_#WY7-bWF(5D?Z(?c+ zJUj|7Ol59obZ8(mIX5?xkz^=;-CRj?8#fld>sPQQmu8I`gB^`bk%G(aq{H+Kp(Fe7}o zosm*RDU9V9Vc{oyMTH7uocaz*x+Nl6w+RXRCQ zqbg_TC;vSm@EDH>?202?!AWY+1i?m=`N=P@@#qp>emR1fbTztvvqR@Fv4SNHhi0Rf zC*OJnO*YoeR)O?lf1D&sn?e{E=xO_VPGn%Pz@5xoYJ zjhGrmpZq1XiyhQ1d0*|y%#&<+zNIBTam}o4QreATTfRhQRU%Rmm<}tkJ!D3Q)-5Se z%4Ibw6+@izl;wSYp_V)r>i6Is`ye{_4wmr<-l;v7F5r%q}_Vbl4I3`jR*jLCD|EQigR2vB16G)M?*<= ztBId@g|VH@44{D5DkKBK%bkEW-iX2paLda*C%|D_WG{2VTsh|K=vD6oNKOR7nSxR$ zClCoDyfqT)8x%eX^AY zupvr+xlKY`2y~WUHv(DDpT5X$KqCpFW#EDz_&O(GV5o7YX`kjeP*xzT!UesG(c_sQ z4>AF?Q?qopJ@hEB?)&cKa0yt~PCfRXtK!wpxkgqDydDfI4Xz9dITR5PjR?Tm$czn2 zhCYlLI^lT@>Y#Cl(u{L7SkNpetq&o9hlD?WuAg(fGy{m)L^ROn3@G^j#Ougiz>}yQ z;&2=$b-|C%a`1y$4*s#S9DB}vJ5PE@$vmq8Yfq@i1xcw3B3wTN09B8-DSpJJ#6yQE z4@YoGSr4xyu_jNBiScMC(~;hl6oUUsxiFJbZ~D#G4t+! zEblLkpR@iFE}XxF3+FErk7%?4J2MnQpB?|0>O81l3IUh~J8f*Aaku{&E9d8xC>DmM zJ990IX0diug44g;Os4K>?FPO{H)URAHw7=;rr;a1sXitF-_MjrjCwqh*{?$q-3~1{ z(Z$_OGr93VrEYfsOM6q#W^SaxD=GnhHWtTBS3uaj?IRH;!2_v)L^hKf$A=5U&|(U3 ziB3Bn)s}5!%dJhjsb_;cZrz)v|J7u^yCI1@uxyCQSwyf+IV zw`lIyEG^91id}59n}ed;%!<|SY}5hzU4fa43ldwin63=&Qh^CI=w3q9B3iV z8P>~sDGH8Wy>2keMy>sE%+ zt;9wpf_KpiT<)Y>1q&vSLp(K`*)Yw*d-JLnjKP>*^pWOUJIcUR2<)T4C_p(X3Oyn( z3(8{BK9t$~^aCb;W*;{vsE^ADh6>JC&&E_2pPzRADlfKR2ZX6}@PNktKN5s#Im(oc z@lifk?1wf09$>0CBKpwTD+pF`HHsRuV$n?+4KR{~B|Ti9Dwl@9;x@hkC!m5BME%urNi3 zO!r4_zh&=#P?*}Wwd|wV+Q+M6$-gRSc(;<*dw3iEfZj$NT^|FT#@p18wZ+xp0vS`b zwa0$*;>C+lBKz*_<@aBmzYK{;%8_J#L!XB+=UpIDx6`gjV*ftf@&Qy)vC29%g1Av^ zf@Ekar8(jl*e2t!yHVbH{yR4*cq#&x>76J2dSFt2Y0cg{PWn7IsRmPe1fj;^sCdqL zaoG#}n)z)#S*Iu7y`0sP-;#|LQ}+9Wj1y(^8uMYDKgR%zGGbttUS%5FLt=D)wG z-!vwFOg4|)m?|1$c3k$0cp>Zt*bNp;InvrGoGxRc+od^;iC>?|;^M7)?dFDycD=du z&+NJ|rhAu`1>r#`3%YT;b;XV(MgS-I60GB|60vGU=p}s0G%`(upxs>hvLrpk9BzW3 z-f~8eD;ld@W-(#%hxc4ef_anq#lGIVm;(NP+;c4@tr^a(^?dJI6^JYNjQg#o!5Zjv zv76<5SuQq{*|hf6_RGeK!^GV2V!oU?(b~+6IG6!8r5}~$P4cQQ7sO<-eKGvmX5lKz zlz~!C%G*6T#ts zm;$<6f%=-YdG_?hH!ojg6IcR6l&(Wfe0IHDHS3#sYD$(D>pB!KoEd5AiVQK!h(Pm!(~w!U)>GUmEnPj5~C1=X^6-mKd7sT*4d<o>(&9#GeCl2f@U^3oV`h_+$>cBLtlah} z>=Q2M&7xkc+f(=G${xhr-Poofmox9O?X2EumWwEH& zo0aYUfv3yWMYFDFRvs!{Jnx8qi2A0bbN+VtmHV5W<8;?L!x=+l!zn{dcn9G3LzK?R z2>}4)Wz3971>ZdJLS)8yf3Ra}0uw4?e+ZPUe!wcqLspR|R?#wKz^+8+Ak41^F-s|l z#MExA_8m@}^N=EvtZAsNFiWGm zYTttCm^;%UikH6SjAlNEDOXK7&s5}jTYHs=pq0k}80Zu3oOMf*ry=c4y>Ozx1OIDia)j{}B$02ubffLbVO zpkZ}T3{2O^;U8TcAqH5dLVDPK9M!qQ9y)KD{fUXF}qsi^PMmhV8EF z-~5&7s+YPj(eKPZ>Ts>B0`C{4c{T8aZfumpaPK0_m@x*Tho`>^_OHrKnxa$ArfpR0qt_WW9lPz5D?etGSmk+ z^|=lAdj%ICcK>`!_12Ql$x97(JG4%m5hZ(&Bw;8a7+siI5r9|Ol59obZ8(nIW#qskz^^CDJusKf9)I4xPY@b zdna++H=P8HK(Q5Db7jjd$(^Ks{XM6Oq)Ba&avOtrm;@w?WG&x!zP)HE>nd3(>nk*r zQB{nUao|)d z)nFvQ)@X3TpfG}q1;++RjZz+{fRWMS&V$z&FVaZ-DoHn7c~}KC5qWWcetGKm<=E*6J@MGwlWC>04Zgq zCwtU@ln%JM#jvcm)+w@NsLy0tTPh*6nL1!2L`8|)=!ra4O+zhZ)*K5}W&;E|YA2NX zz|s^n5U`ewEM*!4REL6eSgy)ufnk+qXGgVRXK$cL>>OZOe~O(WCBV)xxPza6J`fM8 z4?qCMcT$Pp-o3|^Vc>zDR(g;yv}xBUQ|E-T!~`^ zx`BKibWEuw*WADaQa#-$>l^4}yz6F(czpv^i+8gn%DovADCzo! zBCea3q~{w@e=T0kre|+Pi|RuqUOhOj#OdhkqG}mB{dP0Li{ZuSK->p&qwB?-65{X! zaWa}uZ)aztd3GK7(TmafczAF6wfaCtJWvlw|BhvbGmOm&!hRiSoxFn<&MrhX2+)}o z5p??nI|wvSsRsIcUCls3W4sQ7ko@+(F$^@2b}9}Ge++bMU;1Nf2V`h#KMzu&QAWvu zg&J<{t7IC`Tfn3T7G8I2XQQCA5O|PwK*zGZv$g{V_PsX)v#ho2W|lBLuc6^G?bun$ zfrsDPsk4)Apu=~_uu2W!7)(e5slPXI@-Rk(OydWP1Upfotr-OEcQoOiZFViNYtKuD z+BN%*e+70;yYB%uSo4n>sy8c!>eY`Ls@J<2YPX+)avBRuH82$fk*%G*8UP8S{2-;U zQ&$N=2DPIAyE5xJ*Pdwm(rjJmk4hlORSCo&mO$E30;wwj2fH*|B#@S2bw&^119$4G z(KYG6RRX(${=rzi4A6b~gK7RD6%uz+kOz)&1a*1^#S zf{jkbNVp@shpd=kCIpFqOeMffT`E@(^LLiH4de;8?3w>aq%CZCOyn(dnUp8WSAAEc!ZTtw_KOz+uwqC8LOD+KWE8B5#U`3*nx95Y!#qjRjFi&(P17InX~M ze>Y>=>EK^11mZ}vEpzF_N+T2gJ1}($XJ&BC!IV6&U0<2%(|7agF)yNEzVhBMUPAt+^>G0pU1>CLz(qI`|q=E1yJTy%uk#R^^?*ADSLevBu$*$Az8(9h*49_ z%X1>}Y=w!MBKLL;AHZ3|!SdclVgvkfe=VGNtV3j6mW~0J7(XyoeTI(<%{l0SUUSA= z15=;7&w(YMJD7FlLKwv6oB~JXb6^;(&ravJ<=mOGd(L@!=YYc6fo(1!!Ozs^?kFVM zgS6%x>2lVbQIC*_`U0HM6>so6gGSt~CmIIod#`xSXjoA-^qsip1;w)TbQ}-I=4D%dA~( zDo?W9B{ehX{1j;-JpxiX_$%!6nm7>$P-Zo`}%$TERU^cFM619s^@R@ zOt;HBSqtTk(cA(!n#;1vW}n&xf9(W-lkJ%1piu_yj$}u{8iH)t)-?t0l^B{?aeYh> zeXj1I#s>W(=dtBVm#wK%CeH`vSo6zjSO4j_yHop5i}Ig+NCf@6W<5%ujeE8ueG8sD zO9DVFs%Mk^Jpz-?+dy~*9p@Ov5m(k)UcY#@mI7|v$PQ+C2ZnK*j2SSTe^TuG#JXd{ zOCHTp2@GUnxRa$ramqd{Ew(&7V{8jd&GXJPNS}-}i+7wkWE66qP}V4CO2~VeBewiM z^uHV4;&*9-!0>4bh-20eK(}7?W%)o9Mk$hnsv0K6!~0w)Hu>%>rSvJ4nD6I zQba8C&gT@vqaA`8*gD+^Dsq-&xlWFHE}wmwhb068MH+B-PD*ytf09!q#tWZv=8)q` zw1LBN7Qrce-kTs1j7=$)#d1#KOBA!;%_F*?*_`LBo?EBbiO&5vT*lu;+Y;j$%yNr) zABu()nQx;viwPDNh8re`^r!^_92qw#4UHrhYB{e%ZlCg)V7Q9SSMDS^1RuoYJT7E0 z>x?uYvH38l%dD08e`*+1@Qx)9x67H#T2~4RON&V*k@YMXV67}L7|jerx@2VYxlD5K z6%NE%tYzWbjf=VBaXKdfc|PtRTM))Vqgz~91qLi5Xx?IHYJfyh<%GR_Tu@-yDn^Z* z3ZXBHB12Hf1Lw9pOs=j1=(3nAVrW!Z8+mO?JG%b`J1K4Q%CSmW{$`DrJ)T&R2GOf2px9vk_| z$7Orsmo}<>7Kx_y0&3cw)LG^7HI^|Hg=Zb^xT{yve+K-vt<>`&<@}Ip zS#n110k{&JIgVzVH#mIuElTsNleBTG`MF|#Srss!WTnG(Rc(G$J=Q=n%Y|khl&sB9 z=A$(~e_bc^`^s9&`iOn=qpWFyG`$e=y2O__2JpcoUMJ?sGWR2H{P;+>N%|eXlUh~w zQnH(r`7xIJ8DK{fGo{*d znz&GYj$7H7m-KUYB8OhgMa=0O%YN7Y?|c2wf4zBisu#@Dx3r+7UXUd$_7`rXUBwuefBoI}#@`N=ybgnMx`KjZE# zJja2!KfHN78ed#42*82ha(s6PcXc41EQXWu+0pgIWK_um@n$i)`X`j)K>SvZu^u<~ ze=mnKZh!ulI1=~7Lvbt~iKpVZcqQJ5Q!y07v)jc;oQbp1`FJuJinHm})li&^k+>9> z-)=5P*J3Ot;!0eLshD1m#7xY^LM*;a#jW`E--Vy!F=$Z9clVpWhtliErzh_o^-#LE z9i~uVLcY?&N%e3tsfSbCfD`hD4p2HRe*<-2JgC8XES>`&SFvjkp=j#?$j2a$m(a@mC+cxDR?yPTxI!_pFcJ%T4gXND6%UoBuxW z`Tq#MM;W>o0I$(DLRP@~{Ba4eu8u83b~zi3YT*79UkF|oIxXcs&^>2@+1EwUm`rq>;H%nOOpAGgCvjW&*tV^Iz~`?wIF$A3Z;)7x2J693aG?=`Tw6A3>n zO-7#=r8Ub7)N;C++|I=x;*Zf*O;z!|F|8^7h?#gX?$!N zeSYme{62f}_Vnm!pWi3z`DGb!Sn9InAmdZ$bE|vY>OQxJO1H}FH!Aaee<*Qd^shhO z_LR2?H%p$ADoItO7LL1qbDhS8nrAp04EM(E5R@7nj~cCNX=sMT4z)+kqKY(XBsH2) zMQcV-YCl?}widPb-nDAQDz$}5QAN~xj2eIM`+u+Z>-*(g=iJwQUFXv|AAZ+;Vq(Uw zSM8YK*mZrKKRhN@2&oisHhOZIi(u-G4AqI_@?6XiD6a#m&f7Z5$_J_6ydp5x6wswodC2h>|kV~0>jVSK}boN+I4AjkV`RWO*Lvoksi(w2WAAeoaCs_Pn~XPy&K{0 z;jz%sT3+-yVC%j40kmKF{qW$K+$yz>6h=5)JfKuQRcZd~6yXm=q3|O)#}i?3!l(cx ztJ^QbY`~EL@a!2Mgm4F|N^EI&q^n&(q&DajwXrlb-{iY>^s{b>Kr7>}PVc>Ui*etx zE_~woip(#A7M-vhGa96Nxgz|VQ5N&HUB{rEpIBvVz;jO!{%-doRuW4xmHQr9Xq&Ox zd~I`ewe-wycQ@C@h|==c$A9F6jWUEODTafxow)8YsWri04h`5F_Cxw=;F!Z`%OV*; z5E~M7Xh1$n-B{Wq#H6RDE1?(H9efb8pRX1S7uEC()DKJz3?-~72BsTlQOsDT@mD7! zhm5%6OqQdcKCQLLxSdsMAR>fDMPWUsCk;5dS-UMdDz9eG-we?Q%pCi9Anw+$=+#Ut zdxPrJDE+}LerLOFGyDFK-aQtGlaSW7FY`6eiUP(@w23i7zJ;l8b+4M*dS()?h^l{o zF}+gohrIq!QfkIGep}g;mfre)k}{8r-$KV)<*NMR-8TY+QzbQ##SDM{vYIzQ2HXPp z9xp|b&+f7DP(}X(@^;T#iH|uJC9b?q3XDykromq@qev`9~`iEYAEp@3Y!oxpWIO7jxK`z1V zG8Hnvs!y@~)da&O3p(~cyo4Gjd#5B=t$va$vb|l<7AnPiiFz7>o2IwA zwXP`Rq;My3EJA$gG(^fZSCY5NE~-M^=$GTeV|t|3H2#s5b@8;d>+?~|+4a@~_x#qT zhN&<0c3YH^Q^w>5D-m#M>&GlfZz~4){+1O(z~|a$MYg2XM<%+XJ3i;K!q}+?;h6-8 zBw;zybwxHpqonR#l?u)p*daDuH7!)7c3JGt%i)@wL90(V&3FDooP{N~4j(&GhtSK% zbNy2LCl@*^x8!*Wc)r(xAC~Nuil+$MBU*pv6D%s?zcw{J7TI?LDu(yk=6Glb}YfOLub5AL8Cfb%mlk zXG4wV#F3O}qj4j`c;=FG;e_Ek?Gmxf*JGUn`fEb;E$n3*4C_aW(kIneefF|d#-he; zX_o_Osn1$?p?v!=J*q|Sh}86SY1*eKlzJeZe%yy$A6eyVRnlmM9<8CjH5Z|QyQtd? zvK;^FmEiF0WAS+5i>vioY37kNZv_U5B0c9b1^Obib}cl;3>*%Ez-~hzwGs1TuNp5-6v8SuBKsSx->kToX@06@0`%}P27Hx zD8sVK=V+e_6}?v&wZBe{*Og+MnF;eH%3oLCA#y;b;-sA93XKnsMFF5Gy^!FHmFBc zsWJ{HMFIsKv41YAdy25T=KStySM0r<^j@)p^!jB*c@LqzQTotQg~Y)Z)I1N+Mc-R`ODXTp1H`5_fI z8|_R>b}Bd;4{j7yy6$K$_&c-o@TPq+DT|RCdeips%Ryn4l!zg772WtmE%MP+HmyVX zcv+iZizO_kjVhIIh&iSW>VzyK-54ZXcLbD(@s&xjygCVTZ4m`dTE|QeG|hRZ9Tb%m z1{vdZf#AszI`6g^1yOhw>?mM1rp0=_2h$KRyVXoi6quSI1%0rg+X_<#e~k4i*aQo)^FL%T;E~dAqsPs_^db^E$Ayc`f`EP%Vvl?b~=^;Y2d?E#W3P zWpmRLW&l}vK`C{^Wzv+}Y<{0Qustu6w~9kA7uvO1^f5<=jy#$Z5WcYFTz~z?tL`G~ zeZcw_uz^)2q84&e%@zHEX>)s`YZH06rJk zc}DL65Uh*Gx}J;D|CuyItd(H08gMuSAq#~d)Meygsv2t2s#+RavQUJon!FlBRaxYJ zTh7=0C-{eumnB`i2aq5L#_QruO}I4+(QnV37>L}gP=lWW-Pz@%fHLARfK^xD=%?6C z7S?diz)oPD?-RDDD>usNLpg*%PeQJUs2Uo?;1fc1{6iA-wuW;>hG008l;4?yvNw(D+JuF-dYHaTZ{$nq7c5|nF@$i4upW4iR& zz>QOKyOTA1DVt#OmESGRJ<5X-#k3BO4rYVsx`x9Uv2EAn*bM2>@drfs#`u}yvLa`4 zCi7rmW#Dt(;M-WkCK&XE-BG7m_$qI3?t3u0g=6w>wK~mO9|W2&C5peBM*Y2LMy=}z z=Dzt1xaamj3B>fUmZT2O>g4pSt4=YoD|fS~14Z#i5P+Jm->b=eNODv6$AeiGE#rCg z*y2KT)acD;(a6}+Co%2U{7eqo$N3V=AlKhKDsShR<;Jthh6t#!Uw_linka?OdsI_L zTIA!ag_g{h&8(?F+urz2#lU}^Zm4HNjck;)3w)n$;Iczs79U5zX+%GWSVs2 z5UPPp?{Gon1zi#$lswRnswU6%UDI1Ob}uCj6nM+qRG+=$ickk)@D49L+%dMee|7&# doxbp}b@zDXj==)qvNA9T3@9e1X{1FD{6Fkw)A;}Z delta 19647 zcmYJ3Q*_Zu0NeALXase(XN$9|X!Ww$>k#BnpKPGbC=4+cuWH7x}6O%J#dZ z9!>*9i>l;bCr)(<>%Y%eFRsdX-+ascU%L@w0$Tn#zRW(|D{gPS@@sz_cY6B5fFtoR z_8|h7`r}$}D`-;S6%{F??%@GzVbINiwc1-M(aIIkGKjN0Fr(XmdropTpV|r>0ZgbC z^lEpB&!>WeKVK&n1b4zFq0R6HxMtfh1~yXI-D(^p9}h{(YD(DE=K2yM1;HR7SFr`@ zBEvSNd$56)!`B5DzJhkuj;BG=JJv zKxDz)x=eny3l#?E%NzljVeDmzIIx#J3EYqx76Q!Gv4-{0c@(qt#E1JVWk?5i_?giYur<{-^`ly#w+ zC7=nBnYWMsdHPOrTleiBeNkT$JM{Yh71cFIeI6mRptD|ZRX@bL080mYl zlZx;OZK++40pAWt@5|Zk_cGHiBMXG1)S*nvHP$X`#9G7T#w4D36Y;cwKmlPvvn0rx zGSw!t(#FNqEiWViLQAAJv|u`uu6kcl z4`dYsPX%AA02NBmLk9W2=WV*Z#-t2o5bo-Fm=a_D`E>&zz~*^(yUD~fgcyzw#B1x? zbo2zi)g~p3A8oSS_qr%B6@@0bqYkmcV$RDSO#QRR+l|-zi9T3n>~3#1I$U9!qhTXd zku42WAIh={0$7s;c!)tKRV6hHCZ#89YT^70-OMVQB@P}pLEwPQ{uSQ7Mwb3L^8Dgs zGwyhF?iL5sbqPzc1_^sZ;u2j2@d>Qe=HAWY!~U8h2w#Z zVT77MTJrG{kWbpct8W~JR4qOt!uM{BI@KMcJ%2!F_eNd9p3bGw+w4v6*ZoC?p-2eQ zFq(c4D5W)!v*ZpxS7{AGo@9MH59;ue>+^gZgsJy^U)kurn!7DN(YZC}xxViU7`;@j zFC17YTDJzobHDEv3iw3n*2lK#)$K>(LR}>U3d#SH zB+H2-oo`t%*XEWqq3Krvt+7!FKMHLn&uc{8dCEar+s;YUkUMGgg{I+HNjKT<2?xi_ z1gq&f`GK z#PkE`@>L9FfiA}!C1rdp4iWJcs@9N$Xp#bV__MciiJ<67K%p5;WG2_K24bY6 z4R$EGKx$`qi#{y-Y75kTAN)k#U%2CL$b+1ndmG!x&F|YVo2Pb9Hv@r}VK|9nN|1ro zm;wGW@QG5`WmbYI`RhKvTLHwvf;q?{vd`jCjS(tpf~2pVs?2h7xRQU^^i<|XE0yfA zS2`_rE1I^JNhz$=7-*`T9J!(Y*t$(_6Kn_mBQCt;RAS*xKLt^@UMJ~RO&Tcx=K>(`y+g}YqEw{9 z{R|wO)z{_B@R3z?Ep@&BteANpz+j?)V$Aig*Ra!C<;%X^0jr0$sau4%B)lQ)MQziK z$7twBG(>pZopmnjGPUNG$jbuP&2UUaR$2RF%|(Goh`NyownG1gn9_Ojp$Fu zYqs+oP9M~+1PU02<%aPszpV_mM3YdArOoJ_0G0O~gbL8z;?@<_3R42!(Eeq5dIFN+ zw(G7R?I@V?H>_#-%eHW;l<{xVD15{b&s!LuHbMl!tC%~8Xn^lNh=ZCgEpusC07oV+ zk1<9fG0ifceG5&fN?HnDdtv+xr_o@Gic&qas^N7!k!2FNoiwOOoXqJqhIy_miRV!(k8ZG1Sx;Xe z?uQRFJjHCo(vJd<+Y{5*!cv&hB+m?*p>-Xse836`yc8UO0av%u@|g z<&HE0YThql35_c!NPne)PeD9Z7rTXfj?Cx0yX%pG%bFk$#I>j%eumsbxW!XfNV3@{ zQ;#8!qBnAc* zDpZKUT>^YJlbMCW!+lCO=u3THp3Y6MQgYhb?EE+ zHDfdAas|K4s{^1ex}xmeq|9Sdk9WGp3}Gjv^#RUPx(Tr1L}8OL5HXg0E#2#fECq|K z%$0tT}d;7|9o}S*Xm!#TvrM8=u-=PigNag~ofkDJ|&xiuR41i!7S}$(OOxj{haHoDV z6_wX&mxu@vKOt#o&12@t3flat(PzI=C<3Pq|GS0YJtQW&9rD6Hg-Us6F2^St8h~fJ zWs=|s4j3TNk*k*JFUD>P9G0!#LuW5J{O4e3if4xzsh1YA0;#7N>%h+wE>^X#N(uso zdC)u^CtnR2avLK^Q@%PTN^#0$_h~$tvq2a<0#7^`Dh+ zgR8c0c_&H99|(Fo=$ry9>oipo(*HETMTht@4E){jaW*R|fBl&t7Dx?@4XZpKq&h~DS#I5<(W1@d+mv)J3yMdZ z)c0*MY_ME5sL5&j-hcY-3m~R5^jGp$1LLQuiOx2O21~q2hBBzB?#QNzZ7EYLdD(T8 zSU|~<*rix~Q&00E3J8BYq*7$Mf#5alnt*&I&|*@Ea)bTM1LaQ)w4AYugkV?By@-i^ z5Iwg{XFd?5$C7}6f_n@7miJRdn_nn3TQpVyzH8P60o2v_T};Y2qpBS_8G>MG5k21f;`~Cm~KAtnZs^F3i-t! z0_6pxq&2P{B6N-AlR8ET6$Sp#o(?As7d~6RrjA#thHBj~<;w8pkwX`g5QL#tWWd&? zNlf~z&Lmh@EZ!K7N6!8tQdl~-0SOy0BCz*vZ>W9s`aspz0l5mLF11TW#w>Md?<75& zOx|^UwtlDtGKv-P45`&n4w)l(l8@}QAzzpZvpTw2W!?p1X2H8i8&kK`*2Yr51b+R) zae+L`>O2@eeT{m1&AhIlUP8#nLeeS?b_n;Wvjt^-`0~-G(ZJ|=vlh&e+SrAgJRWtPi@K-8ymzV8ZIzcS#Jgf7ZQ;yJ|YZ@ji~j{%Ci!|hT1__ z%DOUe+T45Fnm~`TkJor0=`)yJ0#OrWSpSV^Y3L0N10<`Je-D$5n~Y3Z9GusM9(&QssZ_wzabQkWv%Z246?~q#Mdhk>i%I#^qESe_Q4rY>ow&#|^|1k*zTafuiroMw6sp2@ZSdye>Qf40XgEi zN2xRO9S~)Td_z5#>+Wc^RFn_TZGo8;vTn;+mQ~VdAv3XiB4s4#cvJQ&T@OA~Un5(w z?{SNSFfqY}FP5_pwi&F42&(@MWFZ?h80Dfa+B1qZbrQWOfKLz%FR!+4IzJjtf-U2H zcTULVgoMiyU{r}uw*-$#rj08|+i76uMT8=0&Np%1kB^=_Qw(J*3*JK(T8%B%@ARtlbDm}O?UtE|Op)VZvbt89 z1PWrxv)FrYQ8I9)K#R{!fH5>iOi)Udv6;?g7ZHn{b@Q1sxb7&$0Zyvgwz!(gEb5q8 z*-QSm8eVTcUHk?BPFZ!(aO@{SJ_FfV`q;1>vA{M5p21teaop3NoTB|5c``%uEI@2+ zV3DXCq&6s-9`cQ`GR-R_m01o%HBf189Y8QlHo2l?YdXz2;Y9T=Y-)EKT^=VQZ*C-FWX(9M)*byTgIrYX zMd7Mzpi>m3r<=J+)p>m1AWNlA(%cC1W|I7u1O(VbNC{s!=0M=JZIZHK!l2 zd%Zsa*~(c6ZyaA7)42iJYY3WbFBKZ<9#K2vr3Sul4O$yO0Fc)L837zc_cL$9>!1sx zCYC-iWrf1+m`S|;gyjhOE70#0ab^=6)m}Lnsg%6OEpqp(u|h{Wv@+6<04!P)uXUWW z3}TWNCfDI(zL!1>1YfMrtPDCxjIfR*Qjq86RhZnb6A62Ucr_*FVCwC>9R*l(cKFv2 z3h^N?dA@iBkdKefe{{gx^5%sl^xP~Nt%JJ4k@4{a>_Ax*g+ zIiuh2v_18Gcerdq?tR6@X8OVoe_Q@+A7Cau3LFeFAItwZjaiLQ5H}38DcN0G@mQYf z+ds_&jH!^rV$a)4jB+{r_(K6Y#9yHICm|Bx?})JRp?Ox3la4gLg36|2odK%H8K@4M zZS&?alA)ETCrrWW4fO2x0NHqvJznLu0<4flAd;~=UhOSTI6t|?p>GD;{dmcg%8_ZQ z!u(~$d^7dj(0Ve^I*OW4@9W#uTA-#LV59<=OH`ApLA76>I*Gg%BWYQ!(fHUvMPpY6 z@ITMe!yy7a*U$+*T{e#-8*tTayJ8(@EOS1#yPU53C88_K%`Oq)m1@S^$@I>v&@{;y z>O6lY*o~-@>wql2&Q$dQAE1GpSJTp76d{4%Nw&^z4GLDeO_yX4`D(U5NFP*Uzf$;8(BBJGe3@fX}8aOPBp}C zxp78-wAM$gY{?TegM-G%%Hxlb&GFZ~!pfG@ZWTpsGR*}1l7a3RxMq3Gb@K}wt@V3j zP9t_KN**T}Srb00zl`9$zTT{K1>iWDAOc@IF3;L7H?@Ux)|dv*(rp>zGX((E|6nE6 zP;F_tD#YpO%agKI68NculCeMHd@?<;d--4}aj8lZVdDFWiqG@5MYz`^CouvYicRdj z_v2hg;aD+epmoC{>)MRMKqcF`)L@DdU;I_j86*|rEa0}IvG-Ga`|r8t#XaQ}nPZ;) zSReU?^}}X+h3d+*iLn6Hy$Jy8qw27(Hg%UI8=pX+=(HdNbY>!3UznEZKydj6poz29 z7avg%zgZ7=xgSTVh=DN*t;7jH+g5;1u?I;MU}e6rLe#r1~v{M0vZT8rqp8xK~OGyna^2W^}DLx~Q8v=C^Bx z@0AYJ;Fif3ISQ*k#^##xT(uz-nBg*QwIx3T(BcnuDsD^h3 zyAWeDvE<4Un9${mQQIkoc3+$I=erkZN=eOYqIEe_G@zig;Mj`&iA?cGJi03IYvVn; zd#KKs{U(NR->=CaHH1P>SyA8sW8=)|Cx)T{^yun2ZH{3Etp4Fikex?3380OYB;?=0 zTib?e@a6pMn{y{)kgTF;iTk+tAXH0cbP~hUD?}b}K6RS$p9#~V7*A3?q?MB^NhPC{ zwNai(bPYm%Xv|T(v`(eXKd!X7zpA_CnR&{cHn#F@$ZhZ(_eYqJD z$q^Kmy_F?0t@{E(P%R}rlsay`ZXXy+qs4DjOQcKj$n$j*5=Q|&c(!259O&}kBW@l{ zG<>A;A`?SrWeF4#9EAnP!ZqJ02UsUPST4Bl)`cJhei>ID5vPhkbR;{1%QzAMiy4IO zG*)=;3JsS0et%9F81x0fj!bB@GD5t0u_gDO;(gMWc)X@M0q>;(B$g7XEdm$_6a-_< zB^k_3EF2L3wy!_%)m5v}E2xF5(f6*86grZfU>SxOn^PwFA=_wo9r|>{9b9iV*zNYW zi>uf=ro}(*RYZ%E$13`mfb}K7U-WZM|CjCk<(F)w56hNM<(WTc1_rGK=@!xez2}0n zOSg_+SU!^Fk_|Qw`+HAqXJ<(}7Tayp+m_zEL0shT?{4eO>Ka!>(kIxfkE}Q8+c|ZY zmDU@c5c7>?x#x=i{F$n4gVTLWhTY-jDq>u>z0j>by^p)QjaDXu^xTjEjt*JcPMscx zsowf+7TPQRy&rDKP?AlR*~`+lD;7B~it_$ll`-ASH&QOdm_Gddg7*Fz1N658+l3gM zX)iAbU!@T9M2Yi6FH?@y=F)Bg`>ULvlz-%OPA(Xz#!8)V82$*S&3Wau73aO0Xa;EuqiC#gN91YKY3P@cmRmIwbf{&Xq>#fltRSELSX6;-A31M2uZG|s7jC>sF&G*v14>Df82_EVNpy_ z5r1`rpHh{WmxHbB&~`J9Q=)`PAd>fGu6gUEIQq*KQXQ5LKsmI_Rh^T);w=t_M<#zI zNhbVFSz|qt9GQ>nq>0NhuB5C4^}c?FBDb#joPnJ2>Dw}qq51E#&^G%zVX}2nJ8Z3{ zi`#!WZNJz5i)JT<=lw(Yd%TvMj|SEF`?tDg+rPl9is8ABbrP`O8?+1?DmR2vl~yKN za2%RVEJqT~fF=cp;OQw2PIr|7zFC8;IqRLWdR4f%@JXCgwByXdaBH=4he%-ZDR0js z$*EB&Cz7YMlTZ|*77j-opR9=$1!ub{S}ix9OAwboUN!2c0D(*I!ylb22Meef#>(`}d zCWEIwN1Z0zhawQ;B}c-B0>QuW#X()aFs4J4O7Ji6vh#C6SlM+}tw!zLC;kZOyNyye z-$(~bDg5BP&2FvkF4G6QDK=4{Lx)`-O2L1?I^T-%xoLcckon3qe<;l#W9aV6of%_@ zNErzn5We9o=1NA&6`Mgu#?bSymcuI>Z;h&fPs+ORt1bvfh>>ETlrq36E#G4XHO@3i zyPYBn1;XB>SbBm*0Fs6k3={}fuji(s>)&oaf{^z z_l6n7UYr5CWZH>*8lZ;2zH9I%W`52g%?tGY9XYif3&1srVIQx`2#RC1mawd+i3A;% zAopI%?Yjh>tfYZ0w-gA%!GUob-23E)KM;N6XRuI?*%GDYuLx7V*hty;Qdvc!+c7Z! z`2TWAP`R7>L!bY4h3S+UChQPr$LAO({{0MN-^?<8r_H4#^J$aPb;*;{HGoPDU;ikx z3+GTmC}iPEb#P;4m}&k~u$+$GdBXe~tqm>7^g!o|pi#xuRtW>Wz{>q$Aod$$w2lQq ziV1OKCeUS*aYHk~yh=P53OOYhV4 zJE6G@tB}Y3Ui@@?vx1YF7vr|Z)B2p5-kk+{vtZo^!DRiz9cG3-C(iefr@2 z!E^u5in8NKn`g5UypTeTs}?KK@RyTNy~xlN_N`O9|G~)Lh@gcfB+n`f`y}%yEO?rFIPrV3KT71TVK}=pBwEzXM(YVVR*|vgCjy? zS_DtJSDu^$68*Y<0FG)Z2P`sUMB~Bd`_%|5L^w`zqXn<5CS7(hv^5Za+lZJ0*C1i+ z)QDPnM)~_L7LJvzR8=8v=8u%xxf%`Ef?9f7uyjUo0UT?9&-WYbolzWB@@?b3JBLXl zM<_p_Xm$D+>daW;)sdgfjM7YmUPiWr_|}Tb0~lvjFqho79Ga>jBhB~V-dD<0T9&R9 zDacQRC`VRvsh&(}6`UCugV{{ZC!jwO^{Hrs9?WRo0atD!IRW z2$j+8BBfM~aJaHOkt!bvIw{CJFQ<%-zBmC;a9j`+Oo7hbbC9F8n$lzg-OgLQZ>KSnV7acBPJaV8)LY8h#7Mxo2e%{^Xeo>ySGc=avhU3Z^@SeKl%wG!JS?C#IkChUkq&udvT!-X4v%#59 z**<3E$HdD77k~g3kLHr-9XE4+39jhO*0C$9uQ&?Yjkn@WnAB3L=DAoTA9RV9ky~o80Pt^NYCA!OUPF zc{GA2kqziAsgy_5lRF2>j2I1Yr>9uof#f4Qh5+{M0jpi6Fc%^M;(55l*i%bhZhfZv!KQ8dEC=Rhb60mj$ep#=BGs z8pC$JPRves+vlDF^TSmwv+fNYHosn4dp&MOJ6~&cmvSsy+?kyflrI3#*>zZL59>pB z$8X7C-|Q$pETfC2ty@dCq-n$WuX4%z=Z?I?suwEe(S!_?o*tC5XK zn>}MfE0AAv12tjn1wR0nKAxbvL?SIH2;Ku3ypRV#@onCB9_k=j$WA%~nQC5sLzr_Q ziUcj9vi9a+cyea1l@eD*@YQAD)Q4fm{gE5cq8m(0SvS&%uUEaDTo$1IcEj# zF$>vg>Io$Qv^XC^cpY*Qky1rX1TfIEv~KLak^Qa)GD*P##`30dTi%yz zh8yPjSzQ1u_6)#iAnYL{N_t3LG#3{xXgVICFe^at7X!>0+!N1%xO4L-zUATJ7nc0t z8K+@&k**^_#_l8YO&R-3vuEg^1g!_jN@rvMCKdXo?)j%`1VjPHR#{+a?| zYS5)F-iVHKEcIC3S$+`3ZM#-n{ITu)=JHZqU9ZTnFCVN{#4g%+ z^*5Js(-Lstft%jy81P`{*hI?-#e`>~26>`UUf!IcQ|w$qo*+cXl^kuHfdYusdK=~Q@6I9G7Xh=Je!aQq7GVWN zWreA#bG0>10&n{*#xFRIFy=?LpO3dNYfoEE0-Ck@c9~9!M@seOA}-Pk%Sc{6YQG4;N_soJw2BqcFqiDMu2OCLD+iDN33*3hPQ z-UW(6!XgMsW&4!*p5(V6A=dQ3);&%t{%vagv0On@Elp}n8O5zlk!T{w3*s0E2#`eR zZdkjjS5g=tR45AwXjL%@h@iDv^=LTAz6=b>!jeK`?YN(IN8Ui9@`zCRI5y3mJ78iy z{oG%oJ4Es?!$CcG=t=Go(gBJ`kM;<96O2e^)zGT*au)AC`+HZwYX@G{Jzmy zhvjMzgr^=_aZ1QhX^tQR9pNd98IT{JMJptD^{Hv-(2pd8B1sE|_hCl?vuK)v(m4AW zIGgrSh!_ZsLl+;p}P)LO0Mk-`P0MmLdyRGZ*xccfKFw(PL#v%PF9?WHcE zqgBQ(;w29qHIE(aAZn$=u`Qcgc0@K%AFDs%@H!&Xdo39^p$-4TU@DLq3Q%}w#CluC zok;-Ob@W+hDv;LCBlxvFSVlpX@2SP55gqTtP08J3`e`#mQb<{Zv8)m-xsdd4S$vbt z1jDK1NsysAR514j70~=!<5ti}keC0^OC)iE;K-{A`h|Z#@t?`gO(b}LG!cVB1L@*d zP>$W39!8m`J>w5a>+@Ze4&Zzi4sWOK@Cj*0N&wiFB)vMitB~SgJGQZZ?i>^=zSJ+s z;q(3JLzCq1_PoEzJYy6pL4jVnA0S!)mHohkQn{(CP0;`l-kbo9FFH>1ux-ZhuD7Wa zo1~d!9qs2T9K5dRe)?2)RJql_U|7^cen091PGtAA5j^TTs(_b00$`-(_fEa1e*3J$ ze`XIm2yU0~@ATxZE*hmoW8yg3_K2zFc2{hcDIQ1y^o&;8@b&WGAL_)qJCOS8dm6ml zv~MH~qoaY&g_86|pI*`IRwbmu6&A&f#-Zeu^ z*o0g;G$lgl^B-fp0i>5uYlv71{3{!{mrN&(dX(2AOS7|@E_Wk6o@RrRT^(JIvzOf> zi#_s_!M7x-kc8Wg1`TFW>V2}k5EqBl7JkiSk7z{&yfdxF&J#yc7Tjm$4fWuihax$^XO;cjk%oBZ04fO9yZgMcR7dMK*AOUo7= z#dNR5x0>hW=gE<(L5v%T2A%0V+wEyp!oYTl*Q?AlzaeFLKwI5W$Zj0W4C(B7s6#g(vbVC23EVn0h6L=L3I&$*XbhZ8 zTD)M533+1l;3q*8)!aew!sGUFrrb$k5d;23?f8E40#560ipQp~?f$%gie)8h?aT%6 zy~N>)n{7~YLF2#WSE!p`lv`^fZgzVjvUzViTDGkL`eW(xIWIij&aE$;+ZB|3Ni!K+ zX5qaQ<}+zuc3T{rf?_}|uhVlZ^P5xW_Px3ymyaz?1pNwy5A31~v0QXhCK_YQ;avbK4N2I<#5wXaIZ}?+MwGGEbA1 z8|Ti=OU+FZyS{^P?+9r*z9G?a^4*ZQr8BP>06%6<9YKk{QA9E)d<&+dkGu92J_{i{ zz>3fOZlEog)_FM&PB0vA@0@+qZN zDDEtwbnBO*dfvZ#czuG$op(8N(a%AhGZr_a=>w$VJ|MV6&w8ENhjixiKe7hoo;9?3T@PdSkz!QF z2|_Wvmu7KrK9F%g04#XRaPL&zUkZFMw0e>5l76gurW1ZkAya5bmpmqMCjUV8f=%bY z+Rk?-F9}Vdf2~faf1ufGC?Wr!lmGulo{O!Cz(N)ZjElF)$}R;0l(Xr>Nk0RWyG;ig zl(8S2m6x|IatxI936wQsR|EZxc;{T>n9X1^M7=9^w1ptVS)6}Jb;F0I-m=vmCLize!CMNEMAZ)d+YH; zx?VEB?9aQe5~gg*5*lqAYi2}{Xrs)^z5mH(Ogcv^4SeO~z^Pca4ox+jjYVPHs3s=g zZPgv#b`5PTcUVhU8faVlFGGVpIS-OrW`YRZEosd)9q%~l&RbE;8EjFPIu(dnPE*uU zdO$ZrQ3_6o&%^>|GB?{0c$AA8WC%4+v*}>b9vR-L6$O*e;w|6vc64cS-|8P1I=$2&gxkEkGpMFl!(T4MS=GxCRr(y?_Cy+_MUYs2RHO z3J$3j3~3051L?F5BgxPZ50)|yVht%#G*|lKG@}rGC?=MYN@%2G0T~sq8aDeVY)^q7lR_1mGDX zG>Q9pPK~7LVAotMDp7xg`$DU*~nXaYx&W_(?I~^p#u@AjE(V^{2{40 z`0@Kn*4K@>(k?|dgW+r$L*>0|0nDI8{Jsl)GoVmIy{cUy1ppxwoo=2)nSJ5D0Zn<%|E2zwZ*3p1eJK7y`aKP)3c@fMuei3$$eF|3n-`SpL!6B2{EN z#yfLzr`NA=8w+oF3@#OR%qQ*JTgG76u#c!mJ|UKBb60|7;ixbgPFqMkuTY)s>&!TO>jGYP=K^&Q4%i5J!pc!+BxloL~{ATXw^xX2= zCxjIEUJ)?LcECGAwccbTC|JmOhxAJxnL-*v5_#1`*hwlCKR)R_KpS#)H>&n6cPzgr(_~Ue&Q2x!ly}Le*BFahyz#@vsA^>|w~ITA z)T{rM)Z;O~YGEvq&ygp`UcQ4TW4+xhxcF7J|0q~j+4c(ve$AUSB4D`^Q{3Yy0k0~h z3)i-d*mOmPJkS?WUDSSOV|^Soac&{b^pW3#kyM7+3&U93Qt)I0-N7ik7XHfwc1RQU zl{j5wO%<$=9u&T7Nsj_0y#80MlmfojS=Ah8)n84puHbJ`c0DnUMJQbSo(ZG~Fvl~I zW1A#FcJ^ITd6v;kF#z##5;qHOLlSK>tw6sNb=;K5gP^%ZznnzkZKNNOrL9>nU*=y( z(Xpm^y*<9IIX_FqtC8V~=e zwVQSA!^YFHso=Y#Az=!OOpx}DUF1^ZSfAiWIjcv`naR2KU*%cx!)&41_kypE5?9i|p@G?@lWUq(kmB}Zf>YO(Tvgk|$ak%KDvS3vbXb;t`o`H~`X;EQyi*B#tC$EiqTjTm<_q57mc8Bbghx%D$o^OEDKlnQXUN#{I@7J+AueCS0OI$t0zmAA@ zSWnBHFR$&%)`Qf}dJzx2XHSxz#SNS;BUaxH+-4gpiR=EcgIE0xn$B8iqm|G$T627_ z{0%UuTZYh%Kgx9HI2MPlN|XoO($Cz?bnLWv^1j5HznJhVWGh`1JTthjyFRXZK zZV6Av+xU|x@@1<>yUz&P5}eg*QtOyQ7ic}OlBDxqs73G1JsZX=nmi5M+4~ly^-W~A*2slzI1vnzMffA}^p$4`JB$PT_$*Sa~Vp3&AGUyoheI1C%c~$h+*!Y797%Dhf z68YqUXQry9;(WO?{<=D68Ul9*Fo3p%4)B5=$Ug{mE3r2!y-es`jb4i6Ef*M9t(G!^MhCyZW zM5*1F84cVh-b5G{j^`HDRCNJj18{Pp; zD)iDgb@Nbt-<1=}mthttGEhQ|{m?1i0SrqxlL14F5E)xkXKkNgw&pwXvW|@k#8PKvrU&1zHmOVt z<3m=R&r`}}Pwg&FZ&~S<5)TK*C?6eclbuUvn-0LJGyN-FNqrjO$y)m<`S%X~Br&Ic z=5k<6ddhRxH{jnz1mL+R;KqK*HRWEkg(pYZF}=xWr}@*uFN_;G-iWw||LyGU@O{Ue zD^LHQVc!0eyI%v1!nGIk)dtZ7+95%5z-Yo6DcM#`@i4BJ(;hF<(AVAcfvu}c2L%Xp z@Z{Y#>DF?PcO!u^wUzWL_F@$c@=%V2$CEP6aL{Wlkr{#GHNYtTWb6Br*@fadDq>$p z08M$=?M|*6wcM0p#HwvlE;p*_+>b3Uet7-z+?HWRl9Flos4Zvm{FlOzhJKVCD|*v~ z+|;lb!|?s2!dA4?Pu?4Eejrqqu;yc16L<0B``y`im$YE&wL`?XRc7L!N9Sb0T=&a8 z46fbK?n3a?Xn;8IItH3kCt(=9?^{hfIAj0eP&v}nlmAA3AMop_2a3y!Y|pCH=E)~Z@1a?Qu{`JOO6-&L zyFHJw8t@sMtfv7m9T+~W)=?OZk(awa_+=CQcgSejZw0f8?np?i=jZ-yJf*l-rL7aa zR_s8jq$7ncxvh7zMbI}9|Elo9tL0kmC1KOHAywyc#RjcAaj%kMuU>l}Lw(Sn#*^=Xd!BGtoEIu*d4>DfcpC!U@%M@}R ziQooNgV-C}n9Bj~-riRZ;)1X0U&(jADKRrC5Fo^gbs5XUegc98431?7PF^B+7QPQQ z&H%UhpoXi=YazrNrlEr6RSOUQ-ps^%q4xh%sFg=Uy|7{XL?Ks{BuYeceTEf2FPi${&dKJ;5WMA;7uYrl3C^P?aY4v6jLIBk=k{dA7O*GQYPC)%I@NCGKu-QvvnR{lSXs z-u4NmVVNeFUhegt!g;|NlIw`iFxCAQFngmMw@<|;%F=XdgqIH(}lut9!wAC1EW(ZeP$eR)3af4TDZG)QTVs~}6oorhrz5=Iw>G!?a^hnm%TH2BU`v(Ikei0`aOp8ZGGUZ@nVY7MV!nRN%mJjv7lDIs<<36!QrL>o}}wqtvct;jC%d7 z6!($Cf#;oI;_aF?%M6!-%`wxVoDy;q%}h~qIrQWQ<7sG6R7rWs5tUk-q#>SY6JVq} z(*OB#?cwfMOP-d{;5~~J2H1CT>AUzcCjBIqpr8Y6y4T*>wePfGz_b8ppQtPv^&efN z$}9?-C0*ZYd!hbY5?AY%uYaE0+<52SrXMa`5WlZ^Ya=z4=MMxo% zu+-;z#6?kh(7VHK=1RVmw#vi5c53U#x`<~BotiAGm~Xht(z%q&YQ(wtY!gV5?Ud<) zHwROh=vEB=2fQQz>x|V)7eLV;CdmQFyQj38qo%Sejo}(v&5h1@N}kg@+bhy*CL=y@ zSG~Sx@~J|hPROce?uCAbI>$BW0orx zBn-1s8QLD#r$z{Wml5a>r8%Ce>Y+FOOuQ~jFcoe5zOJu5yFvPA9l>VpKq%G??|;o)h3;!F zZNF@9s)T5F5AQcSE_X>DsINy~-SivQ#+Z5_N<7CxR7MkDNt@#LE`AYpaVA?F=d=&y z+4T?+rva+l)dbIQqQALvP2e~1;cA?fhwfNl{!-G4yJppW`kSmG%*-85OiIslrESa# zq<>3K<|Q_e`eeZ;0=KieQV0u!I%m*Qh?!&j=y4=9WfJGnrzk^de;_a#Jysq$_Ti3T zbzp{qpA4f+zN#n2dgs&vg_25Qsa%aW*`*EuGEEfj%B3!oifNM?hGoAkJw8=)_Zrm% z_kK&SRfUjm7n#gYSE6hjmflPb`V9?jW{^1LN-VGW#p`86;D=F|moB5ZGE}1fWwP#T z*LadNEaTq!^&}b zreWS{;V$=mC)|Ezq=-hpBO%~{{KKMLP>{tF2K`L2ua3vz`JS-2wBQEqZKzJeZfghR zh82S)pPDjjV@|)+tmM;MD(hHji%+6-gZaz|YwY)@On1ia8t9mh&4shvPx-+rYPmRI zYLO4>DZm*z=p)RVFBmi4ylV(VzEMMVGvs}x=Pf=|2(eb^&XMLhRLPi-m}E;DQhFDk>qAOQ*NzDHwwrgKW zg5)-drJ9@o8whrpE&CR4^d)P$yWNKTuA_Y%x-*JcX8YY?^EgM~ z$*qxKDFDutxAFc-nz!t3+#(&NR5gSWswi1aK9w!T`&2>BLZ-woe!*NI5$9`jYQiUJ zWQ`?kYIqrAp=dnP{-l!FSVE7?K7CUoYhgPsDxj5&RJd&M{J6_mBrsZ?X~aG#^V8S; z^8TP{82DFodK|RKenZ)?Pqki@x6-H^HF@+g5MU@N+>q2W2&I9N|JW{j#GsGk19h_} z;lD96h~${AmOuSf!9QPW>@tFXA*4-m*3+jx;;Zg%>XPUpD7lRUmvSpSYMoe zu=TOTj_?j%4LzT}=~iTNPP*yv#t#Wae*xte?)^Huf|T^r@V8H zTtk397byhT``rIH!kigHe9Q^YU|BtG0A(L}MP*D_8dtxae^iiZ=CM#;XR5yPYR<{H zPe_!v_rnuX{=~kA;sEP=VZ2WeS`;)pb$A2L`g*xrY`dTThN^UN?TDy?r&*J075Bbl!fV89TVYcbv`~GR^O%& zrgS%LHeY^L#y0l(j#hKg?vdx{Iq-!G_gvlcEgnUr_p=|e@iM!vx(uy_uTPE-xxFT! z437N_RCF?I4Du5>&lNEUJoM54>%GztQZ`&yDYh~>V|sAjjTU*-UIP}rE8 zVEbVTX$Zbg{M zyeZ^QlV6GbdS!sx)U`CGK?gQDEQY=k^qusjHvJQn*|qocm(e_V3p*BRB)``c6c95t zEf!j5SJ&GETEdf<0$Z+*sxHSC%`m#hT7RbSrb-~Kk5+z*k8W#5PnRcKNZdHb*USz~ z|1Rd9xV@PG#M|DIe3NokKlP2w0CD~3nNahk(MykbUY$#gjqWV1baIUIWcNA+`(3b^ z^G)