Default is to use a total of 4 processors: 4 via shared-memory 1 via Linda Entering Link 1 = C:\G09W\l1.exe PID= 5160. Copyright (c) 1988,1990,1992,1993,1995,1998,2003,2009,2013, Gaussian, Inc. All Rights Reserved. This is part of the Gaussian(R) 09 program. It is based on the Gaussian(R) 03 system (copyright 2003, Gaussian, Inc.), the Gaussian(R) 98 system (copyright 1998, Gaussian, Inc.), the Gaussian(R) 94 system (copyright 1995, Gaussian, Inc.), the Gaussian 92(TM) system (copyright 1992, Gaussian, Inc.), the Gaussian 90(TM) system (copyright 1990, Gaussian, Inc.), the Gaussian 88(TM) system (copyright 1988, Gaussian, Inc.), the Gaussian 86(TM) system (copyright 1986, Carnegie Mellon University), and the Gaussian 82(TM) system (copyright 1983, Carnegie Mellon University). Gaussian is a federally registered trademark of Gaussian, Inc. This software contains proprietary and confidential information, including trade secrets, belonging to Gaussian, Inc. This software is provided under written license and may be used, copied, transmitted, or stored only in accord with that written license. The following legend is applicable only to US Government contracts under FAR: RESTRICTED RIGHTS LEGEND Use, reproduction and disclosure by the US Government is subject to restrictions as set forth in subparagraphs (a) and (c) of the Commercial Computer Software - Restricted Rights clause in FAR 52.227-19. Gaussian, Inc. 340 Quinnipiac St., Bldg. 40, Wallingford CT 06492 --------------------------------------------------------------- Warning -- This program may not be used in any manner that competes with the business of Gaussian, Inc. or will provide assistance to any competitor of Gaussian, Inc. The licensee of this program is prohibited from giving any competitor of Gaussian, Inc. access to this program. By using this program, the user acknowledges that Gaussian, Inc. is engaged in the business of creating and licensing software in the field of computational chemistry and represents and warrants to the licensee that it is not a competitor of Gaussian, Inc. and that it will not use this program in any manner prohibited above. --------------------------------------------------------------- Cite this work as: Gaussian 09, Revision D.01, M. J. Frisch, G. W. Trucks, H. B. Schlegel, G. E. Scuseria, M. A. Robb, J. R. Cheeseman, G. Scalmani, V. Barone, B. Mennucci, G. A. Petersson, H. Nakatsuji, M. Caricato, X. Li, H. P. Hratchian, A. F. Izmaylov, J. Bloino, G. Zheng, J. L. Sonnenberg, M. Hada, M. Ehara, K. Toyota, R. Fukuda, J. Hasegawa, M. Ishida, T. Nakajima, Y. Honda, O. Kitao, H. Nakai, T. Vreven, J. A. Montgomery, Jr., J. E. Peralta, F. Ogliaro, M. Bearpark, J. J. Heyd, E. Brothers, K. N. Kudin, V. N. Staroverov, T. Keith, R. Kobayashi, J. Normand, K. Raghavachari, A. Rendell, J. C. Burant, S. S. Iyengar, J. Tomasi, M. Cossi, N. Rega, J. M. Millam, M. Klene, J. E. Knox, J. B. Cross, V. Bakken, C. Adamo, J. Jaramillo, R. Gomperts, R. E. Stratmann, O. Yazyev, A. J. Austin, R. Cammi, C. Pomelli, J. W. Ochterski, R. L. Martin, K. Morokuma, V. G. Zakrzewski, G. A. Voth, P. Salvador, J. J. Dannenberg, S. Dapprich, A. D. Daniels, O. Farkas, J. B. Foresman, J. V. Ortiz, J. Cioslowski, and D. J. Fox, Gaussian, Inc., Wallingford CT, 2013. ****************************************** Gaussian 09: EM64W-G09RevD.01 13-Apr-2013 13-Oct-2013 ****************************************** %chk=\\ic.ac.uk\homes\ak7611\3rdyearlab\NEW_3RDYEARLAB\New NH3\AK_NH3BH3_OPT_631 g_dp.chk Default route: MaxDisk=10GB ---------------------------------------- # opt b3lyp/6-31g(d,p) geom=connectivity ---------------------------------------- 1/14=-1,18=20,19=15,26=3,38=1,57=2/1,3; 2/9=110,12=2,17=6,18=5,40=1/2; 3/5=1,6=6,7=101,11=2,16=1,25=1,30=1,71=1,74=-5/1,2,3; 4//1; 5/5=2,38=5/2; 6/7=2,8=2,9=2,10=2,28=1/1; 7//1,2,3,16; 1/14=-1,18=20,19=15,26=3/3(2); 2/9=110/2; 99//99; 2/9=110/2; 3/5=1,6=6,7=101,11=2,16=1,25=1,30=1,71=1,74=-5/1,2,3; 4/5=5,16=3,69=1/1; 5/5=2,38=5/2; 7//1,2,3,16; 1/14=-1,18=20,19=15,26=3/3(-5); 2/9=110/2; 6/7=2,8=2,9=2,10=2,19=2,28=1/1; 99/9=1/99; ------------------- NH3BH3 Optimisation ------------------- Symbolic Z-matrix: Charge = 0 Multiplicity = 1 H -2.95249 1.32362 0. H -2.95271 -0.24365 -0.90492 H -2.95277 -0.2437 0.90486 H -0.66163 -0.76599 0.00006 H -0.66211 0.80123 -0.90498 H -0.66209 0.80131 0.90492 N -2.55738 0.27869 0. B -1.05713 0.27869 0. GradGradGradGradGradGradGradGradGradGradGradGradGradGradGradGradGradGrad Berny optimization. Initialization pass. ---------------------------- ! Initial Parameters ! ! (Angstroms and Degrees) ! -------------------------- -------------------------- ! Name Definition Value Derivative Info. ! -------------------------------------------------------------------------------- ! R1 R(1,7) 1.1171 estimate D2E/DX2 ! ! R2 R(2,7) 1.1171 estimate D2E/DX2 ! ! R3 R(3,7) 1.1171 estimate D2E/DX2 ! ! R4 R(4,8) 1.117 estimate D2E/DX2 ! ! R5 R(5,8) 1.1172 estimate D2E/DX2 ! ! R6 R(6,8) 1.1172 estimate D2E/DX2 ! ! R7 R(7,8) 1.5003 estimate D2E/DX2 ! ! A1 A(1,7,2) 108.191 estimate D2E/DX2 ! ! A2 A(1,7,3) 108.1922 estimate D2E/DX2 ! ! A3 A(1,7,8) 110.7128 estimate D2E/DX2 ! ! A4 A(2,7,3) 108.193 estimate D2E/DX2 ! ! A5 A(2,7,8) 110.7248 estimate D2E/DX2 ! ! A6 A(3,7,8) 110.7282 estimate D2E/DX2 ! ! A7 A(4,8,5) 108.1971 estimate D2E/DX2 ! ! A8 A(4,8,6) 108.1957 estimate D2E/DX2 ! ! A9 A(4,8,7) 110.7359 estimate D2E/DX2 ! ! A10 A(5,8,6) 108.1987 estimate D2E/DX2 ! ! A11 A(5,8,7) 110.7071 estimate D2E/DX2 ! ! A12 A(6,8,7) 110.7081 estimate D2E/DX2 ! ! D1 D(1,7,8,4) 179.9967 estimate D2E/DX2 ! ! D2 D(1,7,8,5) -59.9975 estimate D2E/DX2 ! ! D3 D(1,7,8,6) 59.992 estimate D2E/DX2 ! ! D4 D(2,7,8,4) -60.0088 estimate D2E/DX2 ! ! D5 D(2,7,8,5) 59.997 estimate D2E/DX2 ! ! D6 D(2,7,8,6) 179.9866 estimate D2E/DX2 ! ! D7 D(3,7,8,4) 59.9984 estimate D2E/DX2 ! ! D8 D(3,7,8,5) -179.9958 estimate D2E/DX2 ! ! D9 D(3,7,8,6) -60.0063 estimate D2E/DX2 ! -------------------------------------------------------------------------------- Trust Radius=3.00D-01 FncErr=1.00D-07 GrdErr=1.00D-06 Number of steps in this run= 38 maximum allowed number of steps= 100. GradGradGradGradGradGradGradGradGradGradGradGradGradGradGradGradGradGrad Input orientation: --------------------------------------------------------------------- Center Atomic Atomic Coordinates (Angstroms) Number Number Type X Y Z --------------------------------------------------------------------- 1 1 0 -2.952491 1.323620 0.000000 2 1 0 -2.952713 -0.243653 -0.904922 3 1 0 -2.952772 -0.243698 0.904863 4 1 0 -0.661627 -0.765990 0.000060 5 1 0 -0.662105 0.801231 -0.904978 6 1 0 -0.662087 0.801313 0.904921 7 7 0 -2.557377 0.278689 0.000000 8 5 0 -1.057127 0.278689 0.000000 --------------------------------------------------------------------- Distance matrix (angstroms): 1 2 3 4 5 1 H 0.000000 2 H 1.809759 0.000000 3 H 1.809769 1.809785 0.000000 4 H 3.100730 2.518115 2.518095 0.000000 5 H 2.517488 2.517671 3.100735 1.809772 0.000000 6 H 2.517466 3.100718 2.517794 1.809755 1.809899 7 N 1.117137 1.117146 1.117140 2.164537 2.164276 8 B 2.164321 2.164479 2.164517 1.117038 1.117174 6 7 8 6 H 0.000000 7 N 2.164288 0.000000 8 B 1.117173 1.500250 0.000000 Stoichiometry BH6N Framework group C1[X(BH6N)] Deg. of freedom 18 Full point group C1 NOp 1 Largest Abelian subgroup C1 NOp 1 Largest concise Abelian subgroup C1 NOp 1 Standard orientation: --------------------------------------------------------------------- Center Atomic Atomic Coordinates (Angstroms) Number Number Type X Y Z --------------------------------------------------------------------- 1 1 0 -1.061861 0.077248 1.042048 2 1 0 -1.062108 0.863819 -0.587840 3 1 0 -1.062162 -0.941016 -0.454080 4 1 0 1.228973 -0.077301 -1.041871 5 1 0 1.228515 0.941131 0.454148 6 1 0 1.228538 -0.863809 0.588043 7 7 0 -0.666762 -0.000007 -0.000029 8 5 0 0.833488 -0.000005 -0.000050 --------------------------------------------------------------------- Rotational constants (GHZ): 76.5501021 20.0984297 20.0981220 Standard basis: 6-31G(d,p) (6D, 7F) There are 60 symmetry adapted cartesian basis functions of A symmetry. There are 60 symmetry adapted basis functions of A symmetry. 60 basis functions, 98 primitive gaussians, 60 cartesian basis functions 9 alpha electrons 9 beta electrons nuclear repulsion energy 41.7274192595 Hartrees. NAtoms= 8 NActive= 8 NUniq= 8 SFac= 1.00D+00 NAtFMM= 60 NAOKFM=F Big=F Integral buffers will be 131072 words long. Raffenetti 2 integral format. Two-electron integral symmetry is turned on. One-electron integrals computed using PRISM. NBasis= 60 RedAO= T EigKep= 5.69D-03 NBF= 60 NBsUse= 60 1.00D-06 EigRej= -1.00D+00 NBFU= 60 ExpMin= 1.27D-01 ExpMax= 4.17D+03 ExpMxC= 6.27D+02 IAcc=1 IRadAn= 1 AccDes= 0.00D+00 Harris functional with IExCor= 402 and IRadAn= 1 diagonalized for initial guess. HarFok: IExCor= 402 AccDes= 0.00D+00 IRadAn= 1 IDoV= 1 UseB2=F ITyADJ=14 ICtDFT= 3500011 ScaDFX= 1.000000 1.000000 1.000000 1.000000 FoFCou: FMM=F IPFlag= 0 FMFlag= 100000 FMFlg1= 0 NFxFlg= 0 DoJE=T BraDBF=F KetDBF=T FulRan=T wScrn= 0.000000 ICntrl= 500 IOpCl= 0 I1Cent= 200000004 NGrid= 0 NMat0= 1 NMatS0= 1 NMatT0= 0 NMatD0= 1 NMtDS0= 0 NMtDT0= 0 Petite list used in FoFCou. Keep R1 ints in memory in canonical form, NReq=2587035. Requested convergence on RMS density matrix=1.00D-08 within 128 cycles. Requested convergence on MAX density matrix=1.00D-06. Requested convergence on energy=1.00D-06. No special actions if energy rises. Integral accuracy reduced to 1.0D-05 until final iterations. Initial convergence to 1.0D-05 achieved. Increase integral accuracy. SCF Done: E(RB3LYP) = -83.1836748756 A.U. after 11 cycles NFock= 11 Conv=0.65D-08 -V/T= 2.0093 ********************************************************************** Population analysis using the SCF density. ********************************************************************** Orbital symmetries: Occupied (A) (A) (A) (A) (A) (A) (A) (A) (A) Virtual (A) (A) (A) (A) (A) (A) (A) (A) (A) (A) (A) (A) (A) (A) (A) (A) (A) (A) (A) (A) (A) (A) (A) (A) (A) (A) (A) (A) (A) (A) (A) (A) (A) (A) (A) (A) (A) (A) (A) (A) (A) (A) (A) (A) (A) (A) (A) (A) (A) (A) (A) The electronic state is 1-A. Alpha occ. eigenvalues -- -14.43760 -6.62437 -0.92390 -0.52662 -0.52661 Alpha occ. eigenvalues -- -0.51884 -0.36533 -0.25520 -0.25520 Alpha virt. eigenvalues -- -0.00031 0.06732 0.06732 0.23255 0.24781 Alpha virt. eigenvalues -- 0.24785 0.29802 0.45128 0.45129 0.49988 Alpha virt. eigenvalues -- 0.67092 0.69322 0.69326 0.73652 0.75664 Alpha virt. eigenvalues -- 0.75665 0.86743 0.97680 0.97680 1.13698 Alpha virt. eigenvalues -- 1.20114 1.20117 1.43830 1.58540 1.58548 Alpha virt. eigenvalues -- 1.78212 1.94176 1.94191 1.95622 2.01272 Alpha virt. eigenvalues -- 2.01277 2.12759 2.25395 2.25399 2.34313 Alpha virt. eigenvalues -- 2.45714 2.45729 2.57999 2.68565 2.73393 Alpha virt. eigenvalues -- 2.73397 2.87487 2.87496 2.94155 3.25565 Alpha virt. eigenvalues -- 3.25566 3.28269 3.48937 3.48948 3.63269 Alpha virt. eigenvalues -- 4.07179 Condensed to atoms (all electrons): 1 2 3 4 5 6 1 H 0.454176 -0.021390 -0.021390 0.005355 -0.003239 -0.003240 2 H -0.021390 0.454233 -0.021390 -0.003235 -0.003239 0.005361 3 H -0.021390 -0.021390 0.454233 -0.003236 0.005361 -0.003237 4 H 0.005355 -0.003235 -0.003236 0.747931 -0.017992 -0.017991 5 H -0.003239 -0.003239 0.005361 -0.017992 0.748122 -0.017985 6 H -0.003240 0.005361 -0.003237 -0.017991 -0.017985 0.748114 7 N 0.321700 0.321661 0.321659 -0.028894 -0.028929 -0.028925 8 B -0.033114 -0.033113 -0.033109 0.422671 0.422612 0.422612 7 8 1 H 0.321700 -0.033114 2 H 0.321661 -0.033113 3 H 0.321659 -0.033109 4 H -0.028894 0.422671 5 H -0.028929 0.422612 6 H -0.028925 0.422612 7 N 6.402957 0.250717 8 B 0.250717 3.638116 Mulliken charges: 1 1 H 0.301142 2 H 0.301112 3 H 0.301110 4 H -0.104609 5 H -0.104710 6 H -0.104709 7 N -0.531945 8 B -0.057390 Sum of Mulliken charges = 0.00000 Mulliken charges with hydrogens summed into heavy atoms: 1 7 N 0.371418 8 B -0.371418 Electronic spatial extent (au): = 109.4364 Charge= 0.0000 electrons Dipole moment (field-independent basis, Debye): X= -5.6586 Y= 0.0000 Z= -0.0002 Tot= 5.6586 Quadrupole moment (field-independent basis, Debye-Ang): XX= -16.1928 YY= -14.9505 ZZ= -14.9501 XY= 0.0000 XZ= -0.0003 YZ= 0.0001 Traceless Quadrupole moment (field-independent basis, Debye-Ang): XX= -0.8283 YY= 0.4139 ZZ= 0.4144 XY= 0.0000 XZ= -0.0003 YZ= 0.0001 Octapole moment (field-independent basis, Debye-Ang**2): XXX= -15.6014 YYY= -0.4284 ZZZ= 1.8975 XYY= -7.5298 XXY= -0.0002 XXZ= -0.0010 XZZ= -7.5291 YZZ= 0.4284 YYZ= -1.8984 XYZ= 0.0000 Hexadecapole moment (field-independent basis, Debye-Ang**3): XXXX= -95.5973 YYYY= -31.3567 ZZZZ= -31.3519 XXXY= 0.0003 XXXZ= 0.0004 YYYX= 0.2660 YYYZ= 0.0000 ZZZX= -1.1811 ZZZY= 0.0003 XXYY= -20.6864 XXZZ= -20.6863 YYZZ= -10.4517 XXYZ= 0.0001 YYXZ= 1.1793 ZZXY= -0.2662 N-N= 4.172741925952D+01 E-N=-2.756667790343D+02 KE= 8.241615083968D+01 Calling FoFJK, ICntrl= 2127 FMM=F ISym2X=0 I1Cent= 0 IOpClX= 0 NMat=1 NMatS=1 NMatT=0. ***** Axes restored to original set ***** ------------------------------------------------------------------- Center Atomic Forces (Hartrees/Bohr) Number Number X Y Z ------------------------------------------------------------------- 1 1 0.015263168 -0.058058301 -0.000007351 2 1 0.015277569 0.029018253 0.050276130 3 1 0.015274616 0.029018830 -0.050271435 4 1 0.011236977 -0.052360319 0.000004891 5 1 0.011215649 0.026141742 -0.045274008 6 1 0.011202653 0.026148108 0.045271391 7 7 -0.096508317 0.000001979 -0.000003393 8 5 0.017037686 0.000089709 0.000003775 ------------------------------------------------------------------- Cartesian Forces: Max 0.096508317 RMS 0.034760589 GradGradGradGradGradGradGradGradGradGradGradGradGradGradGradGradGradGrad Berny optimization. FormGI is forming the generalized inverse of G from B-inverse, IUseBI=4. Internal Forces: Max 0.059703792 RMS 0.028230653 Search for a local minimum. Step number 1 out of a maximum of 38 All quantities printed in internal units (Hartrees-Bohrs-Radians) Mixed Optimization -- RFO/linear search Second derivative matrix not updated -- first step. ITU= 0 Eigenvalues --- 0.00895 0.05926 0.05926 0.05928 0.05929 Eigenvalues --- 0.16000 0.16000 0.16000 0.16000 0.16000 Eigenvalues --- 0.16000 0.31852 0.31852 0.31855 0.31856 Eigenvalues --- 0.31856 0.31866 0.32351 RFO step: Lambda=-6.01003865D-02 EMin= 8.94965559D-03 Linear search not attempted -- first point. Maximum step size ( 0.300) exceeded in Quadratic search. -- Step size scaled by 0.740 Iteration 1 RMS(Cart)= 0.06059289 RMS(Int)= 0.00109892 Iteration 2 RMS(Cart)= 0.00154965 RMS(Int)= 0.00020257 Iteration 3 RMS(Cart)= 0.00000125 RMS(Int)= 0.00020257 Iteration 4 RMS(Cart)= 0.00000000 RMS(Int)= 0.00020257 Variable Old X -DE/DX Delta X Delta X Delta X New X (Linear) (Quad) (Total) R1 2.11108 -0.05970 0.00000 -0.11673 -0.11673 1.99435 R2 2.11110 -0.05970 0.00000 -0.11673 -0.11673 1.99437 R3 2.11109 -0.05969 0.00000 -0.11671 -0.11671 1.99438 R4 2.11090 0.05295 0.00000 0.10349 0.10349 2.21439 R5 2.11115 0.05287 0.00000 0.10338 0.10338 2.21453 R6 2.11115 0.05286 0.00000 0.10337 0.10337 2.21452 R7 2.83506 0.05069 0.00000 0.09784 0.09784 2.93290 A1 1.88829 -0.00660 0.00000 -0.02219 -0.02240 1.86589 A2 1.88831 -0.00660 0.00000 -0.02218 -0.02239 1.86592 A3 1.93230 0.00632 0.00000 0.02125 0.02102 1.95333 A4 1.88832 -0.00660 0.00000 -0.02220 -0.02241 1.86592 A5 1.93251 0.00631 0.00000 0.02123 0.02101 1.95352 A6 1.93257 0.00631 0.00000 0.02123 0.02101 1.95358 A7 1.88840 0.00845 0.00000 0.02842 0.02803 1.91643 A8 1.88837 0.00846 0.00000 0.02844 0.02805 1.91643 A9 1.93271 -0.00812 0.00000 -0.02731 -0.02769 1.90501 A10 1.88842 0.00843 0.00000 0.02838 0.02800 1.91642 A11 1.93220 -0.00806 0.00000 -0.02711 -0.02749 1.90471 A12 1.93222 -0.00807 0.00000 -0.02716 -0.02754 1.90469 D1 3.14154 0.00000 0.00000 0.00001 0.00001 3.14154 D2 -1.04715 0.00000 0.00000 -0.00002 -0.00002 -1.04717 D3 1.04706 0.00001 0.00000 0.00004 0.00004 1.04709 D4 -1.04735 0.00000 0.00000 0.00000 0.00001 -1.04735 D5 1.04715 0.00000 0.00000 -0.00002 -0.00002 1.04713 D6 3.14136 0.00001 0.00000 0.00003 0.00003 3.14139 D7 1.04717 0.00000 0.00000 -0.00001 -0.00001 1.04716 D8 -3.14152 -0.00001 0.00000 -0.00003 -0.00003 -3.14155 D9 -1.04731 0.00000 0.00000 0.00002 0.00002 -1.04729 Item Value Threshold Converged? Maximum Force 0.059704 0.000450 NO RMS Force 0.028231 0.000300 NO Maximum Displacement 0.124440 0.001800 NO RMS Displacement 0.060113 0.001200 NO Predicted change in Energy=-3.072840D-02 GradGradGradGradGradGradGradGradGradGradGradGradGradGradGradGradGradGrad Input orientation: --------------------------------------------------------------------- Center Atomic Atomic Coordinates (Angstroms) Number Number Type X Y Z --------------------------------------------------------------------- 1 1 0 -2.973542 1.257769 -0.000004 2 1 0 -2.973702 -0.210750 -0.847895 3 1 0 -2.973752 -0.210798 0.847853 4 1 0 -0.643141 -0.828220 0.000056 5 1 0 -0.643501 0.832341 -0.958819 6 1 0 -0.643524 0.832429 0.958756 7 7 0 -2.579580 0.278693 0.000002 8 5 0 -1.027558 0.278735 -0.000005 --------------------------------------------------------------------- Distance matrix (angstroms): 1 2 3 4 5 1 H 0.000000 2 H 1.695721 0.000000 3 H 1.695745 1.695748 0.000000 4 H 3.127639 2.555740 2.555722 0.000000 5 H 2.555271 2.555422 3.127662 1.917526 0.000000 6 H 2.555215 3.127610 2.555501 1.917515 1.917575 7 N 1.055365 1.055375 1.055378 2.230483 2.230307 8 B 2.178385 2.178529 2.178572 1.171805 1.171880 6 7 8 6 H 0.000000 7 N 2.230280 0.000000 8 B 1.171871 1.552022 0.000000 Stoichiometry BH6N Framework group C1[X(BH6N)] Deg. of freedom 18 Full point group C1 NOp 1 Largest Abelian subgroup C1 NOp 1 Largest concise Abelian subgroup C1 NOp 1 Standard orientation: --------------------------------------------------------------------- Center Atomic Atomic Coordinates (Angstroms) Number Number Type X Y Z --------------------------------------------------------------------- 1 1 0 -1.082066 0.187810 0.960883 2 1 0 -1.082269 0.738264 -0.643009 3 1 0 -1.082326 -0.926004 -0.317776 4 1 0 1.248269 -0.212410 -1.086430 5 1 0 1.247967 1.047189 0.359363 6 1 0 1.247935 -0.834759 0.727280 7 7 0 -0.688135 -0.000007 -0.000023 8 5 0 0.863887 -0.000008 -0.000030 --------------------------------------------------------------------- Rotational constants (GHZ): 76.5290475 19.1468055 19.1466672 Standard basis: 6-31G(d,p) (6D, 7F) There are 60 symmetry adapted cartesian basis functions of A symmetry. There are 60 symmetry adapted basis functions of A symmetry. 60 basis functions, 98 primitive gaussians, 60 cartesian basis functions 9 alpha electrons 9 beta electrons nuclear repulsion energy 41.3771112369 Hartrees. NAtoms= 8 NActive= 8 NUniq= 8 SFac= 1.00D+00 NAtFMM= 60 NAOKFM=F Big=F Integral buffers will be 131072 words long. Raffenetti 2 integral format. Two-electron integral symmetry is turned on. One-electron integrals computed using PRISM. NBasis= 60 RedAO= T EigKep= 7.04D-03 NBF= 60 NBsUse= 60 1.00D-06 EigRej= -1.00D+00 NBFU= 60 Initial guess from the checkpoint file: "\\ic.ac.uk\homes\ak7611\3rdyearlab\NEW_3RDYEARLAB\New NH3\AK_NH3BH3_OPT_631g_dp.chk" B after Tr= 0.000000 0.000000 0.000000 Rot= 0.998230 0.059473 -0.000005 0.000004 Ang= 6.82 deg. ExpMin= 1.27D-01 ExpMax= 4.17D+03 ExpMxC= 6.27D+02 IAcc=1 IRadAn= 1 AccDes= 0.00D+00 Harris functional with IExCor= 402 and IRadAn= 1 diagonalized for initial guess. HarFok: IExCor= 402 AccDes= 0.00D+00 IRadAn= 1 IDoV= 1 UseB2=F ITyADJ=14 ICtDFT= 3500011 ScaDFX= 1.000000 1.000000 1.000000 1.000000 FoFCou: FMM=F IPFlag= 0 FMFlag= 100000 FMFlg1= 0 NFxFlg= 0 DoJE=T BraDBF=F KetDBF=T FulRan=T wScrn= 0.000000 ICntrl= 500 IOpCl= 0 I1Cent= 200000004 NGrid= 0 NMat0= 1 NMatS0= 1 NMatT0= 0 NMatD0= 1 NMtDS0= 0 NMtDT0= 0 Petite list used in FoFCou. Keep R1 ints in memory in canonical form, NReq=2587035. Requested convergence on RMS density matrix=1.00D-08 within 128 cycles. Requested convergence on MAX density matrix=1.00D-06. Requested convergence on energy=1.00D-06. No special actions if energy rises. Integral accuracy reduced to 1.0D-05 until final iterations. Initial convergence to 1.0D-05 achieved. Increase integral accuracy. SCF Done: E(RB3LYP) = -83.2156707159 A.U. after 10 cycles NFock= 10 Conv=0.87D-08 -V/T= 2.0098 Calling FoFJK, ICntrl= 2127 FMM=F ISym2X=0 I1Cent= 0 IOpClX= 0 NMat=1 NMatS=1 NMatT=0. ***** Axes restored to original set ***** ------------------------------------------------------------------- Center Atomic Forces (Hartrees/Bohr) Number Number X Y Z ------------------------------------------------------------------- 1 1 0.008843517 -0.024578475 -0.000004952 2 1 0.008843262 0.012290600 0.021286574 3 1 0.008839306 0.012302194 -0.021285654 4 1 0.000934721 -0.020454201 0.000008434 5 1 0.000950094 0.010203658 -0.017686551 6 1 0.000951058 0.010213003 0.017684332 7 7 -0.054455130 -0.000019139 0.000001324 8 5 0.025093173 0.000042360 -0.000003506 ------------------------------------------------------------------- Cartesian Forces: Max 0.054455130 RMS 0.016952509 GradGradGradGradGradGradGradGradGradGradGradGradGradGradGradGradGradGrad Berny optimization. Using GEDIIS/GDIIS optimizer. FormGI is forming the generalized inverse of G from B-inverse, IUseBI=4. Internal Forces: Max 0.027929045 RMS 0.012297729 Search for a local minimum. Step number 2 out of a maximum of 38 All quantities printed in internal units (Hartrees-Bohrs-Radians) Mixed Optimization -- RFO/linear search Update second derivatives using D2CorX and points 1 2 DE= -3.20D-02 DEPred=-3.07D-02 R= 1.04D+00 TightC=F SS= 1.41D+00 RLast= 3.00D-01 DXNew= 5.0454D-01 8.9998D-01 Trust test= 1.04D+00 RLast= 3.00D-01 DXMaxT set to 5.05D-01 ITU= 1 0 Use linear search instead of GDIIS. Eigenvalues --- 0.00895 0.05714 0.05715 0.06207 0.06209 Eigenvalues --- 0.15282 0.16000 0.16000 0.16000 0.16000 Eigenvalues --- 0.16274 0.28460 0.31852 0.31855 0.31856 Eigenvalues --- 0.31861 0.31992 0.34155 RFO step: Lambda=-1.86581264D-03 EMin= 8.94965560D-03 Quartic linear search produced a step of 0.70703. Iteration 1 RMS(Cart)= 0.04352227 RMS(Int)= 0.00174695 Iteration 2 RMS(Cart)= 0.00177184 RMS(Int)= 0.00085013 Iteration 3 RMS(Cart)= 0.00000504 RMS(Int)= 0.00085012 Iteration 4 RMS(Cart)= 0.00000001 RMS(Int)= 0.00085012 Variable Old X -DE/DX Delta X Delta X Delta X New X (Linear) (Quad) (Total) R1 1.99435 -0.02611 -0.08253 -0.00710 -0.08963 1.90472 R2 1.99437 -0.02610 -0.08253 -0.00709 -0.08963 1.90475 R3 1.99438 -0.02611 -0.08252 -0.00712 -0.08964 1.90474 R4 2.21439 0.01963 0.07317 -0.01248 0.06070 2.27509 R5 2.21453 0.01960 0.07309 -0.01243 0.06067 2.27520 R6 2.21452 0.01961 0.07308 -0.01240 0.06068 2.27520 R7 2.93290 0.02793 0.06917 0.03623 0.10540 3.03830 A1 1.86589 -0.00097 -0.01584 0.02010 0.00412 1.87001 A2 1.86592 -0.00097 -0.01583 0.02005 0.00408 1.87000 A3 1.95333 0.00089 0.01486 -0.01857 -0.00386 1.94947 A4 1.86592 -0.00097 -0.01584 0.02017 0.00419 1.87010 A5 1.95352 0.00089 0.01485 -0.01853 -0.00383 1.94969 A6 1.95358 0.00090 0.01485 -0.01844 -0.00375 1.94983 A7 1.91643 0.00643 0.01982 0.03071 0.04857 1.96500 A8 1.91643 0.00643 0.01984 0.03065 0.04852 1.96495 A9 1.90501 -0.00654 -0.01958 -0.03126 -0.05254 1.85248 A10 1.91642 0.00642 0.01980 0.03066 0.04851 1.96494 A11 1.90471 -0.00649 -0.01944 -0.03101 -0.05214 1.85257 A12 1.90469 -0.00649 -0.01947 -0.03091 -0.05207 1.85262 D1 3.14154 0.00000 0.00001 -0.00001 0.00000 3.14154 D2 -1.04717 0.00000 -0.00001 -0.00006 -0.00007 -1.04724 D3 1.04709 0.00000 0.00003 0.00005 0.00007 1.04717 D4 -1.04735 0.00000 0.00000 -0.00006 -0.00006 -1.04741 D5 1.04713 -0.00001 -0.00001 -0.00012 -0.00012 1.04700 D6 3.14139 0.00000 0.00002 0.00000 0.00001 3.14141 D7 1.04716 0.00000 0.00000 0.00005 0.00005 1.04721 D8 -3.14155 0.00000 -0.00002 0.00000 -0.00002 -3.14157 D9 -1.04729 0.00001 0.00001 0.00011 0.00012 -1.04716 Item Value Threshold Converged? Maximum Force 0.027929 0.000450 NO RMS Force 0.012298 0.000300 NO Maximum Displacement 0.093331 0.001800 NO RMS Displacement 0.043657 0.001200 NO Predicted change in Energy=-8.221531D-03 GradGradGradGradGradGradGradGradGradGradGradGradGradGradGradGradGradGrad Input orientation: --------------------------------------------------------------------- Center Atomic Atomic Coordinates (Angstroms) Number Number Type X Y Z --------------------------------------------------------------------- 1 1 0 -2.969643 1.215136 0.000001 2 1 0 -2.969704 -0.189502 -0.811054 3 1 0 -2.969825 -0.189533 0.810991 4 1 0 -0.654334 -0.877609 0.000069 5 1 0 -0.654412 0.857072 -1.001592 6 1 0 -0.654345 0.857153 1.001522 7 7 0 -2.596918 0.278650 0.000008 8 5 0 -0.989119 0.278832 0.000000 --------------------------------------------------------------------- Distance matrix (angstroms): 1 2 3 4 5 1 H 0.000000 2 H 1.621979 0.000000 3 H 1.621974 1.622045 0.000000 4 H 3.120935 2.548009 2.548046 0.000000 5 H 2.547880 2.547980 3.121214 2.003109 0.000000 6 H 2.547902 3.121187 2.548196 2.003074 2.003114 7 N 1.007934 1.007948 1.007943 2.260656 2.260775 8 B 2.190695 2.190859 2.190952 1.203925 1.203984 6 7 8 6 H 0.000000 7 N 2.260815 0.000000 8 B 1.203983 1.607799 0.000000 Stoichiometry BH6N Framework group C1[X(BH6N)] Deg. of freedom 18 Full point group C1 NOp 1 Largest Abelian subgroup C1 NOp 1 Largest concise Abelian subgroup C1 NOp 1 Standard orientation: --------------------------------------------------------------------- Center Atomic Atomic Coordinates (Angstroms) Number Number Type X Y Z --------------------------------------------------------------------- 1 1 0 -1.080844 0.151399 0.924207 2 1 0 -1.081080 0.724711 -0.593071 3 1 0 -1.081202 -0.876004 -0.330881 4 1 0 1.234203 -0.187032 -1.141286 5 1 0 1.234343 1.081885 0.408649 6 1 0 1.234409 -0.894869 0.732554 7 7 0 -0.708235 -0.000009 -0.000006 8 5 0 0.899564 -0.000005 -0.000026 --------------------------------------------------------------------- Rotational constants (GHZ): 75.4829487 18.4353390 18.4351896 Standard basis: 6-31G(d,p) (6D, 7F) There are 60 symmetry adapted cartesian basis functions of A symmetry. There are 60 symmetry adapted basis functions of A symmetry. 60 basis functions, 98 primitive gaussians, 60 cartesian basis functions 9 alpha electrons 9 beta electrons nuclear repulsion energy 41.2023561874 Hartrees. NAtoms= 8 NActive= 8 NUniq= 8 SFac= 1.00D+00 NAtFMM= 60 NAOKFM=F Big=F Integral buffers will be 131072 words long. Raffenetti 2 integral format. Two-electron integral symmetry is turned on. One-electron integrals computed using PRISM. NBasis= 60 RedAO= T EigKep= 7.88D-03 NBF= 60 NBsUse= 60 1.00D-06 EigRej= -1.00D+00 NBFU= 60 Initial guess from the checkpoint file: "\\ic.ac.uk\homes\ak7611\3rdyearlab\NEW_3RDYEARLAB\New NH3\AK_NH3BH3_OPT_631g_dp.chk" B after Tr= 0.000000 0.000000 0.000000 Rot= 0.999883 -0.015321 -0.000033 0.000004 Ang= -1.76 deg. ExpMin= 1.27D-01 ExpMax= 4.17D+03 ExpMxC= 6.27D+02 IAcc=1 IRadAn= 1 AccDes= 0.00D+00 Harris functional with IExCor= 402 and IRadAn= 1 diagonalized for initial guess. HarFok: IExCor= 402 AccDes= 0.00D+00 IRadAn= 1 IDoV= 1 UseB2=F ITyADJ=14 ICtDFT= 3500011 ScaDFX= 1.000000 1.000000 1.000000 1.000000 FoFCou: FMM=F IPFlag= 0 FMFlag= 100000 FMFlg1= 0 NFxFlg= 0 DoJE=T BraDBF=F KetDBF=T FulRan=T wScrn= 0.000000 ICntrl= 500 IOpCl= 0 I1Cent= 200000004 NGrid= 0 NMat0= 1 NMatS0= 1 NMatT0= 0 NMatD0= 1 NMtDS0= 0 NMtDT0= 0 Petite list used in FoFCou. Keep R1 ints in memory in canonical form, NReq=2587035. Requested convergence on RMS density matrix=1.00D-08 within 128 cycles. Requested convergence on MAX density matrix=1.00D-06. Requested convergence on energy=1.00D-06. No special actions if energy rises. Integral accuracy reduced to 1.0D-05 until final iterations. Initial convergence to 1.0D-05 achieved. Increase integral accuracy. SCF Done: E(RB3LYP) = -83.2234695419 A.U. after 10 cycles NFock= 10 Conv=0.93D-08 -V/T= 2.0090 Calling FoFJK, ICntrl= 2127 FMM=F ISym2X=0 I1Cent= 0 IOpClX= 0 NMat=1 NMatS=1 NMatT=0. ***** Axes restored to original set ***** ------------------------------------------------------------------- Center Atomic Forces (Hartrees/Bohr) Number Number X Y Z ------------------------------------------------------------------- 1 1 -0.002660515 0.009930886 -0.000003699 2 1 -0.002655993 -0.004958909 -0.008586554 3 1 -0.002657897 -0.004959466 0.008586705 4 1 -0.000296149 -0.003496308 0.000002808 5 1 -0.000293087 0.001729596 -0.003012692 6 1 -0.000305214 0.001735611 0.003008776 7 7 -0.005438830 0.000012240 0.000005948 8 5 0.014307684 0.000006350 -0.000001291 ------------------------------------------------------------------- Cartesian Forces: Max 0.014307684 RMS 0.004947145 GradGradGradGradGradGradGradGradGradGradGradGradGradGradGradGradGradGrad Berny optimization. Using GEDIIS/GDIIS optimizer. FormGI is forming the generalized inverse of G from B-inverse, IUseBI=4. Internal Forces: Max 0.013413233 RMS 0.004412420 Search for a local minimum. Step number 3 out of a maximum of 38 All quantities printed in internal units (Hartrees-Bohrs-Radians) Mixed Optimization -- RFO/linear search Update second derivatives using D2CorX and points 1 2 3 DE= -7.80D-03 DEPred=-8.22D-03 R= 9.49D-01 TightC=F SS= 1.41D+00 RLast= 2.48D-01 DXNew= 8.4853D-01 7.4465D-01 Trust test= 9.49D-01 RLast= 2.48D-01 DXMaxT set to 7.45D-01 ITU= 1 1 0 Use linear search instead of GDIIS. Eigenvalues --- 0.00895 0.05752 0.05754 0.06735 0.06735 Eigenvalues --- 0.13995 0.16000 0.16000 0.16000 0.16000 Eigenvalues --- 0.16295 0.25170 0.31852 0.31855 0.31856 Eigenvalues --- 0.31861 0.32662 0.45652 RFO step: Lambda=-1.64315315D-03 EMin= 8.94965559D-03 Quartic linear search produced a step of -0.01474. Iteration 1 RMS(Cart)= 0.01458461 RMS(Int)= 0.00014443 Iteration 2 RMS(Cart)= 0.00018158 RMS(Int)= 0.00005347 Iteration 3 RMS(Cart)= 0.00000003 RMS(Int)= 0.00005347 Variable Old X -DE/DX Delta X Delta X Delta X New X (Linear) (Quad) (Total) R1 1.90472 0.01021 0.00132 0.01904 0.02036 1.92508 R2 1.90475 0.01020 0.00132 0.01900 0.02033 1.92507 R3 1.90474 0.01020 0.00132 0.01901 0.02033 1.92507 R4 2.27509 0.00327 -0.00089 0.01928 0.01839 2.29347 R5 2.27520 0.00326 -0.00089 0.01923 0.01834 2.29354 R6 2.27520 0.00325 -0.00089 0.01922 0.01832 2.29352 R7 3.03830 0.01341 -0.00155 0.05238 0.05083 3.08913 A1 1.87001 0.00114 -0.00006 0.00594 0.00587 1.87587 A2 1.87000 0.00114 -0.00006 0.00596 0.00588 1.87589 A3 1.94947 -0.00105 0.00006 -0.00542 -0.00538 1.94409 A4 1.87010 0.00114 -0.00006 0.00589 0.00581 1.87591 A5 1.94969 -0.00106 0.00006 -0.00553 -0.00549 1.94421 A6 1.94983 -0.00106 0.00006 -0.00553 -0.00549 1.94434 A7 1.96500 0.00141 -0.00072 0.01322 0.01238 1.97738 A8 1.96495 0.00141 -0.00072 0.01325 0.01240 1.97735 A9 1.85248 -0.00162 0.00077 -0.01524 -0.01457 1.83790 A10 1.96494 0.00142 -0.00072 0.01325 0.01241 1.97734 A11 1.85257 -0.00162 0.00077 -0.01520 -0.01454 1.83803 A12 1.85262 -0.00164 0.00077 -0.01532 -0.01466 1.83796 D1 3.14154 0.00000 0.00000 -0.00005 -0.00004 3.14150 D2 -1.04724 0.00000 0.00000 -0.00003 -0.00003 -1.04727 D3 1.04717 0.00000 0.00000 -0.00003 -0.00003 1.04714 D4 -1.04741 0.00000 0.00000 0.00000 0.00000 -1.04741 D5 1.04700 0.00000 0.00000 0.00001 0.00002 1.04702 D6 3.14141 0.00000 0.00000 0.00002 0.00002 3.14142 D7 1.04721 0.00000 0.00000 -0.00011 -0.00011 1.04710 D8 -3.14157 0.00000 0.00000 -0.00009 -0.00009 3.14153 D9 -1.04716 0.00000 0.00000 -0.00009 -0.00009 -1.04726 Item Value Threshold Converged? Maximum Force 0.013413 0.000450 NO RMS Force 0.004412 0.000300 NO Maximum Displacement 0.034739 0.001800 NO RMS Displacement 0.014524 0.001200 NO Predicted change in Energy=-8.264540D-04 GradGradGradGradGradGradGradGradGradGradGradGradGradGradGradGradGradGrad Input orientation: --------------------------------------------------------------------- Center Atomic Atomic Coordinates (Angstroms) Number Number Type X Y Z --------------------------------------------------------------------- 1 1 0 -2.977026 1.227198 0.000031 2 1 0 -2.977020 -0.195506 -0.821461 3 1 0 -2.977114 -0.195585 0.821402 4 1 0 -0.650293 -0.891754 0.000078 5 1 0 -0.650317 0.864106 -1.013836 6 1 0 -0.650362 0.864225 1.013733 7 7 0 -2.605432 0.278681 0.000016 8 5 0 -0.970736 0.278833 -0.000019 --------------------------------------------------------------------- Distance matrix (angstroms): 1 2 3 4 5 1 H 0.000000 2 H 1.642844 0.000000 3 H 1.642852 1.642864 0.000000 4 H 3.147005 2.563854 2.563849 0.000000 5 H 2.563852 2.563851 3.147225 2.027576 0.000000 6 H 2.563729 3.147110 2.563974 2.027549 2.027569 7 N 1.018708 1.018703 1.018701 2.278703 2.278835 8 B 2.219143 2.219221 2.219316 1.213655 1.213687 6 7 8 6 H 0.000000 7 N 2.278767 0.000000 8 B 1.213678 1.634696 0.000000 Stoichiometry BH6N Framework group C1[X(BH6N)] Deg. of freedom 18 Full point group C1 NOp 1 Largest Abelian subgroup C1 NOp 1 Largest concise Abelian subgroup C1 NOp 1 Standard orientation: --------------------------------------------------------------------- Center Atomic Atomic Coordinates (Angstroms) Number Number Type X Y Z --------------------------------------------------------------------- 1 1 0 -1.089494 0.445545 0.837401 2 1 0 -1.089614 0.502502 -0.804456 3 1 0 -1.089746 -0.947882 -0.032842 4 1 0 1.237024 -0.549965 -1.033405 5 1 0 1.237201 1.169894 0.040442 6 1 0 1.237109 -0.620021 0.992933 7 7 0 -0.717997 -0.000002 -0.000002 8 5 0 0.916700 -0.000011 -0.000012 --------------------------------------------------------------------- Rotational constants (GHZ): 73.6352786 17.9777405 17.9776813 Standard basis: 6-31G(d,p) (6D, 7F) There are 60 symmetry adapted cartesian basis functions of A symmetry. There are 60 symmetry adapted basis functions of A symmetry. 60 basis functions, 98 primitive gaussians, 60 cartesian basis functions 9 alpha electrons 9 beta electrons nuclear repulsion energy 40.7244743813 Hartrees. NAtoms= 8 NActive= 8 NUniq= 8 SFac= 1.00D+00 NAtFMM= 60 NAOKFM=F Big=F Integral buffers will be 131072 words long. Raffenetti 2 integral format. Two-electron integral symmetry is turned on. One-electron integrals computed using PRISM. NBasis= 60 RedAO= T EigKep= 8.24D-03 NBF= 60 NBsUse= 60 1.00D-06 EigRej= -1.00D+00 NBFU= 60 Initial guess from the checkpoint file: "\\ic.ac.uk\homes\ak7611\3rdyearlab\NEW_3RDYEARLAB\New NH3\AK_NH3BH3_OPT_631g_dp.chk" B after Tr= 0.000000 0.000000 0.000000 Rot= 0.986696 0.162574 0.000011 0.000006 Ang= 18.71 deg. ExpMin= 1.27D-01 ExpMax= 4.17D+03 ExpMxC= 6.27D+02 IAcc=1 IRadAn= 1 AccDes= 0.00D+00 Harris functional with IExCor= 402 and IRadAn= 1 diagonalized for initial guess. HarFok: IExCor= 402 AccDes= 0.00D+00 IRadAn= 1 IDoV= 1 UseB2=F ITyADJ=14 ICtDFT= 3500011 ScaDFX= 1.000000 1.000000 1.000000 1.000000 FoFCou: FMM=F IPFlag= 0 FMFlag= 100000 FMFlg1= 0 NFxFlg= 0 DoJE=T BraDBF=F KetDBF=T FulRan=T wScrn= 0.000000 ICntrl= 500 IOpCl= 0 I1Cent= 200000004 NGrid= 0 NMat0= 1 NMatS0= 1 NMatT0= 0 NMatD0= 1 NMtDS0= 0 NMtDT0= 0 Petite list used in FoFCou. Keep R1 ints in memory in canonical form, NReq=2587035. Requested convergence on RMS density matrix=1.00D-08 within 128 cycles. Requested convergence on MAX density matrix=1.00D-06. Requested convergence on energy=1.00D-06. No special actions if energy rises. Integral accuracy reduced to 1.0D-05 until final iterations. Initial convergence to 1.0D-05 achieved. Increase integral accuracy. SCF Done: E(RB3LYP) = -83.2244528670 A.U. after 10 cycles NFock= 10 Conv=0.62D-08 -V/T= 2.0102 Calling FoFJK, ICntrl= 2127 FMM=F ISym2X=0 I1Cent= 0 IOpClX= 0 NMat=1 NMatS=1 NMatT=0. ***** Axes restored to original set ***** ------------------------------------------------------------------- Center Atomic Forces (Hartrees/Bohr) Number Number X Y Z ------------------------------------------------------------------- 1 1 0.000241006 0.000468257 -0.000005704 2 1 0.000240153 -0.000239638 -0.000408410 3 1 0.000260925 -0.000239276 0.000409292 4 1 -0.000715810 0.001121112 0.000001908 5 1 -0.000732994 -0.000564467 0.000984183 6 1 -0.000722506 -0.000565366 -0.000974228 7 7 -0.007348687 0.000020408 0.000003788 8 5 0.008777914 -0.000001030 -0.000010830 ------------------------------------------------------------------- Cartesian Forces: Max 0.008777914 RMS 0.002391759 GradGradGradGradGradGradGradGradGradGradGradGradGradGradGradGradGradGrad Berny optimization. Using GEDIIS/GDIIS optimizer. FormGI is forming the generalized inverse of G from B-inverse, IUseBI=4. Internal Forces: Max 0.006606602 RMS 0.001352964 Search for a local minimum. Step number 4 out of a maximum of 38 All quantities printed in internal units (Hartrees-Bohrs-Radians) Mixed Optimization -- RFO/linear search Update second derivatives using D2CorX and points 1 2 3 4 DE= -9.83D-04 DEPred=-8.26D-04 R= 1.19D+00 TightC=F SS= 1.41D+00 RLast= 7.83D-02 DXNew= 1.2523D+00 2.3481D-01 Trust test= 1.19D+00 RLast= 7.83D-02 DXMaxT set to 7.45D-01 ITU= 1 1 1 0 Use linear search instead of GDIIS. Eigenvalues --- 0.00895 0.05807 0.05809 0.06878 0.06879 Eigenvalues --- 0.11243 0.16000 0.16000 0.16000 0.16000 Eigenvalues --- 0.16280 0.21350 0.31852 0.31855 0.31856 Eigenvalues --- 0.31861 0.37486 0.45650 RFO step: Lambda=-2.02097851D-04 EMin= 8.94966143D-03 Quartic linear search produced a step of 0.27502. Iteration 1 RMS(Cart)= 0.00675790 RMS(Int)= 0.00005857 Iteration 2 RMS(Cart)= 0.00004416 RMS(Int)= 0.00004446 Iteration 3 RMS(Cart)= 0.00000001 RMS(Int)= 0.00004446 Variable Old X -DE/DX Delta X Delta X Delta X New X (Linear) (Quad) (Total) R1 1.92508 0.00035 0.00560 -0.00541 0.00019 1.92527 R2 1.92507 0.00035 0.00559 -0.00538 0.00021 1.92528 R3 1.92507 0.00034 0.00559 -0.00542 0.00017 1.92524 R4 2.29347 -0.00127 0.00506 -0.00390 0.00115 2.29463 R5 2.29354 -0.00129 0.00504 -0.00398 0.00106 2.29460 R6 2.29352 -0.00128 0.00504 -0.00393 0.00111 2.29463 R7 3.08913 0.00661 0.01398 0.02575 0.03973 3.12886 A1 1.87587 0.00038 0.00161 0.00234 0.00394 1.87981 A2 1.87589 0.00039 0.00162 0.00244 0.00404 1.87993 A3 1.94409 -0.00035 -0.00148 -0.00212 -0.00361 1.94048 A4 1.87591 0.00039 0.00160 0.00236 0.00395 1.87986 A5 1.94421 -0.00036 -0.00151 -0.00218 -0.00370 1.94050 A6 1.94434 -0.00038 -0.00151 -0.00238 -0.00390 1.94044 A7 1.97738 0.00045 0.00341 0.00414 0.00744 1.98483 A8 1.97735 0.00045 0.00341 0.00404 0.00735 1.98470 A9 1.83790 -0.00053 -0.00401 -0.00484 -0.00893 1.82898 A10 1.97734 0.00045 0.00341 0.00410 0.00740 1.98474 A11 1.83803 -0.00055 -0.00400 -0.00499 -0.00907 1.82896 A12 1.83796 -0.00055 -0.00403 -0.00494 -0.00905 1.82890 D1 3.14150 0.00000 -0.00001 0.00012 0.00011 -3.14158 D2 -1.04727 0.00000 -0.00001 0.00019 0.00018 -1.04709 D3 1.04714 0.00000 -0.00001 0.00015 0.00014 1.04728 D4 -1.04741 0.00000 0.00000 0.00018 0.00018 -1.04723 D5 1.04702 0.00001 0.00000 0.00024 0.00024 1.04726 D6 3.14142 0.00000 0.00000 0.00020 0.00021 -3.14156 D7 1.04710 0.00000 -0.00003 0.00008 0.00005 1.04715 D8 3.14153 0.00000 -0.00003 0.00014 0.00012 -3.14154 D9 -1.04726 0.00000 -0.00003 0.00011 0.00008 -1.04717 Item Value Threshold Converged? Maximum Force 0.006607 0.000450 NO RMS Force 0.001353 0.000300 NO Maximum Displacement 0.024762 0.001800 NO RMS Displacement 0.006738 0.001200 NO Predicted change in Energy=-1.415974D-04 GradGradGradGradGradGradGradGradGradGradGradGradGradGradGradGradGradGrad Input orientation: --------------------------------------------------------------------- Center Atomic Atomic Coordinates (Angstroms) Number Number Type X Y Z --------------------------------------------------------------------- 1 1 0 -2.981513 1.228698 -0.000081 2 1 0 -2.981482 -0.196294 -0.822693 3 1 0 -2.981401 -0.196254 0.822706 4 1 0 -0.647553 -0.895197 0.000048 5 1 0 -0.647660 0.865861 -1.016729 6 1 0 -0.647704 0.865839 1.016727 7 7 0 -2.613354 0.278735 -0.000009 8 5 0 -0.957632 0.278809 -0.000025 --------------------------------------------------------------------- Distance matrix (angstroms): 1 2 3 4 5 1 H 0.000000 2 H 1.645385 0.000000 3 H 1.645439 1.645400 0.000000 4 H 3.155677 2.571496 2.571407 0.000000 5 H 2.571399 2.571487 3.155623 2.033510 0.000000 6 H 2.571426 3.155647 2.571347 2.033441 2.033457 7 N 1.018808 1.018813 1.018791 2.289649 2.289626 8 B 2.235706 2.235726 2.235668 1.214265 1.214249 6 7 8 6 H 0.000000 7 N 2.289590 0.000000 8 B 1.214267 1.655722 0.000000 Stoichiometry BH6N Framework group C1[X(BH6N)] Deg. of freedom 18 Full point group C1 NOp 1 Largest Abelian subgroup C1 NOp 1 Largest concise Abelian subgroup C1 NOp 1 Standard orientation: --------------------------------------------------------------------- Center Atomic Atomic Coordinates (Angstroms) Number Number Type X Y Z --------------------------------------------------------------------- 1 1 0 -1.094308 0.760736 0.568977 2 1 0 -1.094344 0.112403 -0.943291 3 1 0 -1.094252 -0.873145 0.374296 4 1 0 1.239560 -0.940134 -0.703186 5 1 0 1.239522 1.079083 -0.462519 6 1 0 1.239492 -0.138956 1.165771 7 7 0 -0.726190 -0.000006 0.000004 8 5 0 0.929532 0.000012 -0.000016 --------------------------------------------------------------------- Rotational constants (GHZ): 73.2868224 17.6690589 17.6688889 Standard basis: 6-31G(d,p) (6D, 7F) There are 60 symmetry adapted cartesian basis functions of A symmetry. There are 60 symmetry adapted basis functions of A symmetry. 60 basis functions, 98 primitive gaussians, 60 cartesian basis functions 9 alpha electrons 9 beta electrons nuclear repulsion energy 40.5181090373 Hartrees. NAtoms= 8 NActive= 8 NUniq= 8 SFac= 1.00D+00 NAtFMM= 60 NAOKFM=F Big=F Integral buffers will be 131072 words long. Raffenetti 2 integral format. Two-electron integral symmetry is turned on. One-electron integrals computed using PRISM. NBasis= 60 RedAO= T EigKep= 8.39D-03 NBF= 60 NBsUse= 60 1.00D-06 EigRej= -1.00D+00 NBFU= 60 Initial guess from the checkpoint file: "\\ic.ac.uk\homes\ak7611\3rdyearlab\NEW_3RDYEARLAB\New NH3\AK_NH3BH3_OPT_631g_dp.chk" B after Tr= 0.000000 0.000000 0.000000 Rot= 0.975943 0.218027 0.000007 -0.000023 Ang= 25.19 deg. ExpMin= 1.27D-01 ExpMax= 4.17D+03 ExpMxC= 6.27D+02 IAcc=1 IRadAn= 1 AccDes= 0.00D+00 Harris functional with IExCor= 402 and IRadAn= 1 diagonalized for initial guess. HarFok: IExCor= 402 AccDes= 0.00D+00 IRadAn= 1 IDoV= 1 UseB2=F ITyADJ=14 ICtDFT= 3500011 ScaDFX= 1.000000 1.000000 1.000000 1.000000 FoFCou: FMM=F IPFlag= 0 FMFlag= 100000 FMFlg1= 0 NFxFlg= 0 DoJE=T BraDBF=F KetDBF=T FulRan=T wScrn= 0.000000 ICntrl= 500 IOpCl= 0 I1Cent= 200000004 NGrid= 0 NMat0= 1 NMatS0= 1 NMatT0= 0 NMatD0= 1 NMtDS0= 0 NMtDT0= 0 Petite list used in FoFCou. Keep R1 ints in memory in canonical form, NReq=2587035. Requested convergence on RMS density matrix=1.00D-08 within 128 cycles. Requested convergence on MAX density matrix=1.00D-06. Requested convergence on energy=1.00D-06. No special actions if energy rises. Integral accuracy reduced to 1.0D-05 until final iterations. Initial convergence to 1.0D-05 achieved. Increase integral accuracy. SCF Done: E(RB3LYP) = -83.2246395572 A.U. after 9 cycles NFock= 9 Conv=0.41D-08 -V/T= 2.0105 Calling FoFJK, ICntrl= 2127 FMM=F ISym2X=0 I1Cent= 0 IOpClX= 0 NMat=1 NMatS=1 NMatT=0. ***** Axes restored to original set ***** ------------------------------------------------------------------- Center Atomic Forces (Hartrees/Bohr) Number Number X Y Z ------------------------------------------------------------------- 1 1 0.000193255 0.000005403 0.000002126 2 1 0.000209607 -0.000003592 -0.000003212 3 1 0.000189458 0.000001109 0.000014081 4 1 -0.000472350 0.001721918 -0.000002770 5 1 -0.000466510 -0.000862813 0.001485926 6 1 -0.000480611 -0.000863376 -0.001491445 7 7 -0.002812380 -0.000011983 -0.000015418 8 5 0.003639531 0.000013334 0.000010713 ------------------------------------------------------------------- Cartesian Forces: Max 0.003639531 RMS 0.001133445 GradGradGradGradGradGradGradGradGradGradGradGradGradGradGradGradGradGrad Berny optimization. Using GEDIIS/GDIIS optimizer. FormGI is forming the generalized inverse of G from B-inverse, IUseBI=4. Internal Forces: Max 0.002220061 RMS 0.000724204 Search for a local minimum. Step number 5 out of a maximum of 38 All quantities printed in internal units (Hartrees-Bohrs-Radians) Mixed Optimization -- En-DIIS/RFO-DIIS Update second derivatives using D2CorX and points 1 2 3 4 5 DE= -1.87D-04 DEPred=-1.42D-04 R= 1.32D+00 TightC=F SS= 1.41D+00 RLast= 4.56D-02 DXNew= 1.2523D+00 1.3683D-01 Trust test= 1.32D+00 RLast= 4.56D-02 DXMaxT set to 7.45D-01 ITU= 1 1 1 1 0 Use linear search instead of GDIIS. Eigenvalues --- 0.00895 0.05846 0.05846 0.06966 0.06966 Eigenvalues --- 0.09155 0.16000 0.16000 0.16000 0.16001 Eigenvalues --- 0.16235 0.19499 0.31852 0.31855 0.31856 Eigenvalues --- 0.31861 0.35062 0.47001 RFO step: Lambda=-3.51371899D-05 EMin= 8.94959436D-03 Quartic linear search produced a step of 0.48791. Iteration 1 RMS(Cart)= 0.00368194 RMS(Int)= 0.00001799 Iteration 2 RMS(Cart)= 0.00000531 RMS(Int)= 0.00001721 Iteration 3 RMS(Cart)= 0.00000000 RMS(Int)= 0.00001721 Variable Old X -DE/DX Delta X Delta X Delta X New X (Linear) (Quad) (Total) R1 1.92527 -0.00006 0.00009 -0.00004 0.00006 1.92533 R2 1.92528 -0.00007 0.00010 -0.00009 0.00001 1.92529 R3 1.92524 -0.00006 0.00008 0.00000 0.00008 1.92532 R4 2.29463 -0.00178 0.00056 -0.00548 -0.00492 2.28971 R5 2.29460 -0.00178 0.00052 -0.00544 -0.00492 2.28967 R6 2.29463 -0.00179 0.00054 -0.00551 -0.00497 2.28966 R7 3.12886 0.00222 0.01939 0.00229 0.02167 3.15053 A1 1.87981 0.00018 0.00192 0.00059 0.00250 1.88232 A2 1.87993 0.00017 0.00197 0.00042 0.00239 1.88232 A3 1.94048 -0.00017 -0.00176 -0.00052 -0.00229 1.93819 A4 1.87986 0.00018 0.00193 0.00059 0.00251 1.88237 A5 1.94050 -0.00018 -0.00181 -0.00058 -0.00240 1.93811 A6 1.94044 -0.00016 -0.00190 -0.00040 -0.00231 1.93813 A7 1.98483 0.00001 0.00363 -0.00080 0.00279 1.98762 A8 1.98470 0.00002 0.00358 -0.00067 0.00288 1.98758 A9 1.82898 -0.00003 -0.00436 0.00080 -0.00359 1.82538 A10 1.98474 0.00002 0.00361 -0.00064 0.00293 1.98767 A11 1.82896 -0.00002 -0.00443 0.00097 -0.00349 1.82547 A12 1.82890 -0.00003 -0.00442 0.00085 -0.00360 1.82530 D1 -3.14158 0.00000 0.00005 0.00019 0.00024 -3.14134 D2 -1.04709 0.00000 0.00009 0.00009 0.00018 -1.04691 D3 1.04728 0.00000 0.00007 0.00019 0.00026 1.04754 D4 -1.04723 0.00000 0.00009 0.00019 0.00027 -1.04696 D5 1.04726 0.00000 0.00012 0.00009 0.00021 1.04747 D6 -3.14156 0.00000 0.00010 0.00019 0.00029 -3.14126 D7 1.04715 0.00000 0.00002 0.00027 0.00030 1.04745 D8 -3.14154 0.00000 0.00006 0.00018 0.00023 -3.14131 D9 -1.04717 0.00000 0.00004 0.00028 0.00032 -1.04686 Item Value Threshold Converged? Maximum Force 0.002220 0.000450 NO RMS Force 0.000724 0.000300 NO Maximum Displacement 0.012704 0.001800 NO RMS Displacement 0.003682 0.001200 NO Predicted change in Energy=-3.560362D-05 GradGradGradGradGradGradGradGradGradGradGradGradGradGradGradGradGradGrad Input orientation: --------------------------------------------------------------------- Center Atomic Atomic Coordinates (Angstroms) Number Number Type X Y Z --------------------------------------------------------------------- 1 1 0 -2.984106 1.229567 -0.000190 2 1 0 -2.983957 -0.196863 -0.823376 3 1 0 -2.983951 -0.196547 0.823572 4 1 0 -0.645694 -0.893758 -0.000106 5 1 0 -0.645712 0.865216 -1.015462 6 1 0 -0.645870 0.865009 1.015548 7 7 0 -2.618099 0.278741 -0.000004 8 5 0 -0.950910 0.278832 -0.000038 --------------------------------------------------------------------- Distance matrix (angstroms): 1 2 3 4 5 1 H 0.000000 2 H 1.646919 0.000000 3 H 1.646932 1.646949 0.000000 4 H 3.158588 2.575055 2.575266 0.000000 5 H 2.575193 2.575325 3.158600 2.030994 0.000000 6 H 2.575262 3.158459 2.574952 2.030964 2.031011 7 N 1.018838 1.018820 1.018834 2.294589 2.294651 8 B 2.244501 2.244427 2.244455 1.211662 1.211643 6 7 8 6 H 0.000000 7 N 2.294503 0.000000 8 B 1.211638 1.667190 0.000000 Stoichiometry BH6N Framework group C1[X(BH6N)] Deg. of freedom 18 Full point group C1 NOp 1 Largest Abelian subgroup C1 NOp 1 Largest concise Abelian subgroup C1 NOp 1 Standard orientation: --------------------------------------------------------------------- Center Atomic Atomic Coordinates (Angstroms) Number Number Type X Y Z --------------------------------------------------------------------- 1 1 0 -1.096814 0.291634 0.905009 2 1 0 -1.096703 0.637990 -0.705077 3 1 0 -1.096741 -0.929592 -0.199977 4 1 0 1.241509 -0.359347 -1.116167 5 1 0 1.241592 1.146274 0.246933 6 1 0 1.241379 -0.787045 0.869252 7 7 0 -0.730848 0.000009 -0.000007 8 5 0 0.936342 0.000005 0.000016 --------------------------------------------------------------------- Rotational constants (GHZ): 73.3410241 17.5079689 17.5078660 Standard basis: 6-31G(d,p) (6D, 7F) There are 60 symmetry adapted cartesian basis functions of A symmetry. There are 60 symmetry adapted basis functions of A symmetry. 60 basis functions, 98 primitive gaussians, 60 cartesian basis functions 9 alpha electrons 9 beta electrons nuclear repulsion energy 40.4284235045 Hartrees. NAtoms= 8 NActive= 8 NUniq= 8 SFac= 1.00D+00 NAtFMM= 60 NAOKFM=F Big=F Integral buffers will be 131072 words long. Raffenetti 2 integral format. Two-electron integral symmetry is turned on. One-electron integrals computed using PRISM. NBasis= 60 RedAO= T EigKep= 8.43D-03 NBF= 60 NBsUse= 60 1.00D-06 EigRej= -1.00D+00 NBFU= 60 Initial guess from the checkpoint file: "\\ic.ac.uk\homes\ak7611\3rdyearlab\NEW_3RDYEARLAB\New NH3\AK_NH3BH3_OPT_631g_dp.chk" B after Tr= 0.000000 0.000000 0.000000 Rot= 0.952784 -0.303647 0.000010 0.000011 Ang= -35.35 deg. Keep R1 ints in memory in canonical form, NReq=2587035. Requested convergence on RMS density matrix=1.00D-08 within 128 cycles. Requested convergence on MAX density matrix=1.00D-06. Requested convergence on energy=1.00D-06. No special actions if energy rises. SCF Done: E(RB3LYP) = -83.2246851695 A.U. after 8 cycles NFock= 8 Conv=0.10D-08 -V/T= 2.0105 Calling FoFJK, ICntrl= 2127 FMM=F ISym2X=0 I1Cent= 0 IOpClX= 0 NMat=1 NMatS=1 NMatT=0. ***** Axes restored to original set ***** ------------------------------------------------------------------- Center Atomic Forces (Hartrees/Bohr) Number Number X Y Z ------------------------------------------------------------------- 1 1 0.000132188 -0.000259718 0.000009901 2 1 0.000126045 0.000118796 0.000214008 3 1 0.000132119 0.000122405 -0.000220329 4 1 -0.000161857 0.000791967 -0.000007812 5 1 -0.000172722 -0.000389827 0.000676628 6 1 -0.000160000 -0.000388364 -0.000681282 7 7 -0.000433302 0.000018972 -0.000008440 8 5 0.000537528 -0.000014232 0.000017327 ------------------------------------------------------------------- Cartesian Forces: Max 0.000791967 RMS 0.000332484 GradGradGradGradGradGradGradGradGradGradGradGradGradGradGradGradGradGrad Berny optimization. Using GEDIIS/GDIIS optimizer. FormGI is forming the generalized inverse of G from B-inverse, IUseBI=4. Internal Forces: Max 0.000806394 RMS 0.000279622 Search for a local minimum. Step number 6 out of a maximum of 38 All quantities printed in internal units (Hartrees-Bohrs-Radians) Mixed Optimization -- En-DIIS/RFO-DIIS Swapping is turned off. Update second derivatives using D2CorX and points 2 3 4 5 6 DE= -4.56D-05 DEPred=-3.56D-05 R= 1.28D+00 TightC=F SS= 1.41D+00 RLast= 2.53D-02 DXNew= 1.2523D+00 7.5937D-02 Trust test= 1.28D+00 RLast= 2.53D-02 DXMaxT set to 7.45D-01 ITU= 1 1 1 1 1 0 Eigenvalues --- 0.00895 0.05869 0.05869 0.07000 0.07001 Eigenvalues --- 0.08397 0.16000 0.16000 0.16000 0.16001 Eigenvalues --- 0.16111 0.18797 0.27892 0.31852 0.31856 Eigenvalues --- 0.31856 0.31862 0.47982 En-DIIS/RFO-DIIS IScMMF= 0 using points: 6 5 RFO step: Lambda=-4.05042776D-06. DidBck=F Rises=F RFO-DIIS coefs: 1.18752 -0.18752 Iteration 1 RMS(Cart)= 0.00102234 RMS(Int)= 0.00000085 Iteration 2 RMS(Cart)= 0.00000019 RMS(Int)= 0.00000084 Variable Old X -DE/DX Delta X Delta X Delta X New X (Linear) (Quad) (Total) R1 1.92533 -0.00029 0.00001 -0.00045 -0.00044 1.92488 R2 1.92529 -0.00027 0.00000 -0.00039 -0.00039 1.92490 R3 1.92532 -0.00028 0.00002 -0.00044 -0.00042 1.92490 R4 2.28971 -0.00081 -0.00092 -0.00207 -0.00299 2.28672 R5 2.28967 -0.00080 -0.00092 -0.00205 -0.00297 2.28670 R6 2.28966 -0.00080 -0.00093 -0.00204 -0.00297 2.28670 R7 3.15053 0.00004 0.00406 -0.00244 0.00163 3.15216 A1 1.88232 0.00003 0.00047 -0.00009 0.00038 1.88270 A2 1.88232 0.00003 0.00045 -0.00014 0.00030 1.88262 A3 1.93819 -0.00003 -0.00043 0.00014 -0.00029 1.93790 A4 1.88237 0.00003 0.00047 -0.00017 0.00030 1.88267 A5 1.93811 -0.00003 -0.00045 0.00016 -0.00029 1.93781 A6 1.93813 -0.00003 -0.00043 0.00009 -0.00034 1.93779 A7 1.98762 -0.00004 0.00052 -0.00073 -0.00020 1.98741 A8 1.98758 -0.00004 0.00054 -0.00066 -0.00012 1.98745 A9 1.82538 0.00006 -0.00067 0.00095 0.00027 1.82565 A10 1.98767 -0.00004 0.00055 -0.00073 -0.00018 1.98749 A11 1.82547 0.00004 -0.00065 0.00075 0.00009 1.82556 A12 1.82530 0.00006 -0.00068 0.00095 0.00028 1.82558 D1 -3.14134 0.00000 0.00005 -0.00010 -0.00006 -3.14139 D2 -1.04691 0.00000 0.00003 -0.00016 -0.00012 -1.04703 D3 1.04754 -0.00001 0.00005 -0.00021 -0.00016 1.04738 D4 -1.04696 0.00001 0.00005 -0.00002 0.00003 -1.04693 D5 1.04747 0.00000 0.00004 -0.00007 -0.00004 1.04743 D6 -3.14126 0.00000 0.00005 -0.00013 -0.00008 -3.14134 D7 1.04745 0.00000 0.00006 -0.00007 -0.00002 1.04743 D8 -3.14131 0.00000 0.00004 -0.00013 -0.00008 -3.14139 D9 -1.04686 0.00000 0.00006 -0.00018 -0.00012 -1.04698 Item Value Threshold Converged? Maximum Force 0.000806 0.000450 NO RMS Force 0.000280 0.000300 YES Maximum Displacement 0.002997 0.001800 NO RMS Displacement 0.001022 0.001200 YES Predicted change in Energy=-3.850284D-06 GradGradGradGradGradGradGradGradGradGradGradGradGradGradGradGradGradGrad Input orientation: --------------------------------------------------------------------- Center Atomic Atomic Coordinates (Angstroms) Number Number Type X Y Z --------------------------------------------------------------------- 1 1 0 -2.984236 1.229483 -0.000143 2 1 0 -2.984097 -0.196798 -0.823338 3 1 0 -2.984107 -0.196517 0.823443 4 1 0 -0.645457 -0.892172 -0.000084 5 1 0 -0.645604 0.864374 -1.014086 6 1 0 -0.645621 0.864257 1.014170 7 7 0 -2.618614 0.278758 -0.000024 8 5 0 -0.950564 0.278811 0.000007 --------------------------------------------------------------------- Distance matrix (angstroms): 1 2 3 4 5 1 H 0.000000 2 H 1.646793 0.000000 3 H 1.646746 1.646781 0.000000 4 H 3.157738 2.574982 2.575153 0.000000 5 H 2.574992 2.575078 3.157626 2.028215 0.000000 6 H 2.575138 3.157648 2.574904 2.028239 2.028256 7 N 1.018605 1.018615 1.018612 2.294434 2.294349 8 B 2.244905 2.244852 2.244834 1.210079 1.210070 6 7 8 6 H 0.000000 7 N 2.294363 0.000000 8 B 1.210067 1.668051 0.000000 Stoichiometry BH6N Framework group C1[X(BH6N)] Deg. of freedom 18 Full point group C1 NOp 1 Largest Abelian subgroup C1 NOp 1 Largest concise Abelian subgroup C1 NOp 1 Standard orientation: --------------------------------------------------------------------- Center Atomic Atomic Coordinates (Angstroms) Number Number Type X Y Z --------------------------------------------------------------------- 1 1 0 -1.096868 -0.439186 0.843210 2 1 0 -1.096772 0.949887 -0.041325 3 1 0 -1.096759 -0.510690 -0.801983 4 1 0 1.241864 0.541139 -1.038439 5 1 0 1.241742 0.628782 0.987882 6 1 0 1.241752 -1.169923 0.050608 7 7 0 -0.731268 0.000006 0.000000 8 5 0 0.936783 -0.000010 0.000010 --------------------------------------------------------------------- Rotational constants (GHZ): 73.4667591 17.4992917 17.4991831 Standard basis: 6-31G(d,p) (6D, 7F) There are 60 symmetry adapted cartesian basis functions of A symmetry. There are 60 symmetry adapted basis functions of A symmetry. 60 basis functions, 98 primitive gaussians, 60 cartesian basis functions 9 alpha electrons 9 beta electrons nuclear repulsion energy 40.4347193778 Hartrees. NAtoms= 8 NActive= 8 NUniq= 8 SFac= 1.00D+00 NAtFMM= 60 NAOKFM=F Big=F Integral buffers will be 131072 words long. Raffenetti 2 integral format. Two-electron integral symmetry is turned on. One-electron integrals computed using PRISM. NBasis= 60 RedAO= T EigKep= 8.41D-03 NBF= 60 NBsUse= 60 1.00D-06 EigRej= -1.00D+00 NBFU= 60 Initial guess from the checkpoint file: "\\ic.ac.uk\homes\ak7611\3rdyearlab\NEW_3RDYEARLAB\New NH3\AK_NH3BH3_OPT_631g_dp.chk" B after Tr= 0.000000 0.000000 0.000000 Rot= 0.922643 -0.385654 0.000007 -0.000013 Ang= -45.37 deg. Keep R1 ints in memory in canonical form, NReq=2587035. Requested convergence on RMS density matrix=1.00D-08 within 128 cycles. Requested convergence on MAX density matrix=1.00D-06. Requested convergence on energy=1.00D-06. No special actions if energy rises. SCF Done: E(RB3LYP) = -83.2246900711 A.U. after 7 cycles NFock= 7 Conv=0.14D-08 -V/T= 2.0104 Calling FoFJK, ICntrl= 2127 FMM=F ISym2X=0 I1Cent= 0 IOpClX= 0 NMat=1 NMatS=1 NMatT=0. ***** Axes restored to original set ***** ------------------------------------------------------------------- Center Atomic Forces (Hartrees/Bohr) Number Number X Y Z ------------------------------------------------------------------- 1 1 0.000061587 -0.000101858 0.000001757 2 1 0.000066618 0.000058029 0.000096703 3 1 0.000051550 0.000048553 -0.000094636 4 1 -0.000058718 0.000128568 -0.000001779 5 1 -0.000044289 -0.000059687 0.000106291 6 1 -0.000051240 -0.000058268 -0.000112525 7 7 -0.000085893 -0.000014983 -0.000008544 8 5 0.000060385 -0.000000353 0.000012732 ------------------------------------------------------------------- Cartesian Forces: Max 0.000128568 RMS 0.000068394 GradGradGradGradGradGradGradGradGradGradGradGradGradGradGradGradGradGrad Berny optimization. Using GEDIIS/GDIIS optimizer. FormGI is forming the generalized inverse of G from B-inverse, IUseBI=4. Internal Forces: Max 0.000138998 RMS 0.000063456 Search for a local minimum. Step number 7 out of a maximum of 38 All quantities printed in internal units (Hartrees-Bohrs-Radians) Mixed Optimization -- En-DIIS/RFO-DIIS Swapping is turned off. Update second derivatives using D2CorX and points 2 3 4 5 6 7 DE= -4.90D-06 DEPred=-3.85D-06 R= 1.27D+00 TightC=F SS= 1.41D+00 RLast= 5.54D-03 DXNew= 1.2523D+00 1.6620D-02 Trust test= 1.27D+00 RLast= 5.54D-03 DXMaxT set to 7.45D-01 ITU= 1 1 1 1 1 1 0 Eigenvalues --- 0.00895 0.05873 0.05873 0.06995 0.06999 Eigenvalues --- 0.07947 0.15958 0.16000 0.16000 0.16003 Eigenvalues --- 0.16115 0.20102 0.23374 0.31853 0.31856 Eigenvalues --- 0.31861 0.31868 0.45632 En-DIIS/RFO-DIIS IScMMF= 0 using points: 7 6 5 RFO step: Lambda=-2.10625407D-07. DidBck=F Rises=F RFO-DIIS coefs: 1.24197 -0.27187 0.02990 Iteration 1 RMS(Cart)= 0.00033849 RMS(Int)= 0.00000011 Iteration 2 RMS(Cart)= 0.00000003 RMS(Int)= 0.00000011 Variable Old X -DE/DX Delta X Delta X Delta X New X (Linear) (Quad) (Total) R1 1.92488 -0.00012 -0.00011 -0.00016 -0.00027 1.92461 R2 1.92490 -0.00013 -0.00009 -0.00023 -0.00032 1.92458 R3 1.92490 -0.00012 -0.00010 -0.00017 -0.00027 1.92463 R4 2.28672 -0.00014 -0.00058 -0.00003 -0.00061 2.28611 R5 2.28670 -0.00013 -0.00057 0.00001 -0.00057 2.28613 R6 2.28670 -0.00014 -0.00057 -0.00003 -0.00060 2.28609 R7 3.15216 -0.00009 -0.00025 -0.00008 -0.00034 3.15182 A1 1.88270 0.00002 0.00002 0.00015 0.00016 1.88286 A2 1.88262 0.00002 0.00000 0.00014 0.00014 1.88277 A3 1.93790 -0.00002 0.00000 -0.00018 -0.00018 1.93772 A4 1.88267 0.00001 0.00000 0.00013 0.00013 1.88280 A5 1.93781 -0.00002 0.00000 -0.00017 -0.00017 1.93765 A6 1.93779 -0.00001 -0.00001 -0.00005 -0.00006 1.93773 A7 1.98741 0.00002 -0.00013 0.00023 0.00010 1.98751 A8 1.98745 0.00003 -0.00012 0.00030 0.00019 1.98764 A9 1.82565 -0.00004 0.00017 -0.00043 -0.00025 1.82540 A10 1.98749 0.00001 -0.00013 0.00021 0.00008 1.98757 A11 1.82556 -0.00002 0.00013 -0.00024 -0.00011 1.82545 A12 1.82558 -0.00002 0.00017 -0.00027 -0.00010 1.82548 D1 -3.14139 0.00000 -0.00002 -0.00011 -0.00013 -3.14153 D2 -1.04703 0.00000 -0.00004 -0.00015 -0.00018 -1.04722 D3 1.04738 0.00000 -0.00005 -0.00014 -0.00019 1.04719 D4 -1.04693 0.00000 0.00000 -0.00016 -0.00016 -1.04709 D5 1.04743 0.00000 -0.00001 -0.00020 -0.00021 1.04722 D6 -3.14134 0.00000 -0.00003 -0.00019 -0.00021 -3.14155 D7 1.04743 0.00000 -0.00001 -0.00014 -0.00015 1.04728 D8 -3.14139 0.00000 -0.00003 -0.00017 -0.00020 3.14159 D9 -1.04698 0.00000 -0.00004 -0.00016 -0.00020 -1.04718 Item Value Threshold Converged? Maximum Force 0.000139 0.000450 YES RMS Force 0.000063 0.000300 YES Maximum Displacement 0.000771 0.001800 YES RMS Displacement 0.000338 0.001200 YES Predicted change in Energy=-2.028054D-07 Optimization completed. -- Stationary point found. ---------------------------- ! Optimized Parameters ! ! (Angstroms and Degrees) ! -------------------------- -------------------------- ! Name Definition Value Derivative Info. ! -------------------------------------------------------------------------------- ! R1 R(1,7) 1.0186 -DE/DX = -0.0001 ! ! R2 R(2,7) 1.0186 -DE/DX = -0.0001 ! ! R3 R(3,7) 1.0186 -DE/DX = -0.0001 ! ! R4 R(4,8) 1.2101 -DE/DX = -0.0001 ! ! R5 R(5,8) 1.2101 -DE/DX = -0.0001 ! ! R6 R(6,8) 1.2101 -DE/DX = -0.0001 ! ! R7 R(7,8) 1.6681 -DE/DX = -0.0001 ! ! A1 A(1,7,2) 107.8706 -DE/DX = 0.0 ! ! A2 A(1,7,3) 107.8664 -DE/DX = 0.0 ! ! A3 A(1,7,8) 111.0335 -DE/DX = 0.0 ! ! A4 A(2,7,3) 107.869 -DE/DX = 0.0 ! ! A5 A(2,7,8) 111.0286 -DE/DX = 0.0 ! ! A6 A(3,7,8) 111.0274 -DE/DX = 0.0 ! ! A7 A(4,8,5) 113.8703 -DE/DX = 0.0 ! ! A8 A(4,8,6) 113.8726 -DE/DX = 0.0 ! ! A9 A(4,8,7) 104.6022 -DE/DX = 0.0 ! ! A10 A(5,8,6) 113.8747 -DE/DX = 0.0 ! ! A11 A(5,8,7) 104.597 -DE/DX = 0.0 ! ! A12 A(6,8,7) 104.598 -DE/DX = 0.0 ! ! D1 D(1,7,8,4) -179.9885 -DE/DX = 0.0 ! ! D2 D(1,7,8,5) -59.9906 -DE/DX = 0.0 ! ! D3 D(1,7,8,6) 60.0105 -DE/DX = 0.0 ! ! D4 D(2,7,8,4) -59.9844 -DE/DX = 0.0 ! ! D5 D(2,7,8,5) 60.0136 -DE/DX = 0.0 ! ! D6 D(2,7,8,6) -179.9854 -DE/DX = 0.0 ! ! D7 D(3,7,8,4) 60.0135 -DE/DX = 0.0 ! ! D8 D(3,7,8,5) 180.0115 -DE/DX = 0.0 ! ! D9 D(3,7,8,6) -59.9875 -DE/DX = 0.0 ! -------------------------------------------------------------------------------- GradGradGradGradGradGradGradGradGradGradGradGradGradGradGradGradGradGrad Input orientation: --------------------------------------------------------------------- Center Atomic Atomic Coordinates (Angstroms) Number Number Type X Y Z --------------------------------------------------------------------- 1 1 0 -2.984236 1.229483 -0.000143 2 1 0 -2.984097 -0.196798 -0.823338 3 1 0 -2.984107 -0.196517 0.823443 4 1 0 -0.645457 -0.892172 -0.000084 5 1 0 -0.645604 0.864374 -1.014086 6 1 0 -0.645621 0.864257 1.014170 7 7 0 -2.618614 0.278758 -0.000024 8 5 0 -0.950564 0.278811 0.000007 --------------------------------------------------------------------- Distance matrix (angstroms): 1 2 3 4 5 1 H 0.000000 2 H 1.646793 0.000000 3 H 1.646746 1.646781 0.000000 4 H 3.157738 2.574982 2.575153 0.000000 5 H 2.574992 2.575078 3.157626 2.028215 0.000000 6 H 2.575138 3.157648 2.574904 2.028239 2.028256 7 N 1.018605 1.018615 1.018612 2.294434 2.294349 8 B 2.244905 2.244852 2.244834 1.210079 1.210070 6 7 8 6 H 0.000000 7 N 2.294363 0.000000 8 B 1.210067 1.668051 0.000000 Stoichiometry BH6N Framework group C1[X(BH6N)] Deg. of freedom 18 Full point group C1 NOp 1 Largest Abelian subgroup C1 NOp 1 Largest concise Abelian subgroup C1 NOp 1 Standard orientation: --------------------------------------------------------------------- Center Atomic Atomic Coordinates (Angstroms) Number Number Type X Y Z --------------------------------------------------------------------- 1 1 0 -1.096868 -0.439186 0.843210 2 1 0 -1.096772 0.949887 -0.041325 3 1 0 -1.096759 -0.510690 -0.801983 4 1 0 1.241864 0.541139 -1.038439 5 1 0 1.241742 0.628782 0.987882 6 1 0 1.241752 -1.169923 0.050608 7 7 0 -0.731268 0.000006 0.000000 8 5 0 0.936783 -0.000010 0.000010 --------------------------------------------------------------------- Rotational constants (GHZ): 73.4667591 17.4992917 17.4991831 ********************************************************************** Population analysis using the SCF density. ********************************************************************** Orbital symmetries: Occupied (A) (A) (A) (A) (A) (A) (A) (A) (A) Virtual (A) (A) (A) (A) (A) (A) (A) (A) (A) (A) (A) (A) (A) (A) (A) (A) (A) (A) (A) (A) (A) (A) (A) (A) (A) (A) (A) (A) (A) (A) (A) (A) (A) (A) (A) (A) (A) (A) (A) (A) (A) (A) (A) (A) (A) (A) (A) (A) (A) (A) (A) The electronic state is 1-A. Alpha occ. eigenvalues -- -14.41344 -6.67466 -0.94739 -0.54785 -0.54784 Alpha occ. eigenvalues -- -0.50377 -0.34682 -0.26699 -0.26698 Alpha virt. eigenvalues -- 0.02811 0.10580 0.10580 0.18568 0.22063 Alpha virt. eigenvalues -- 0.22063 0.24955 0.45501 0.45501 0.47855 Alpha virt. eigenvalues -- 0.65293 0.65295 0.66862 0.78871 0.80132 Alpha virt. eigenvalues -- 0.80133 0.88737 0.95652 0.95654 0.99940 Alpha virt. eigenvalues -- 1.18497 1.18499 1.44146 1.54899 1.54902 Alpha virt. eigenvalues -- 1.66067 1.76068 1.76069 2.00514 2.08658 Alpha virt. eigenvalues -- 2.18090 2.18091 2.27025 2.27028 2.29434 Alpha virt. eigenvalues -- 2.44305 2.44309 2.44800 2.69147 2.69149 Alpha virt. eigenvalues -- 2.72445 2.90638 2.90641 3.04012 3.16336 Alpha virt. eigenvalues -- 3.21872 3.21873 3.40163 3.40165 3.63709 Alpha virt. eigenvalues -- 4.11334 Condensed to atoms (all electrons): 1 2 3 4 5 6 1 H 0.418967 -0.021356 -0.021359 0.003399 -0.001440 -0.001438 2 H -0.021356 0.418964 -0.021355 -0.001440 -0.001438 0.003400 3 H -0.021359 -0.021355 0.418973 -0.001438 0.003400 -0.001440 4 H 0.003399 -0.001440 -0.001438 0.766720 -0.020041 -0.020040 5 H -0.001440 -0.001438 0.003400 -0.020041 0.766724 -0.020033 6 H -0.001438 0.003400 -0.001440 -0.020040 -0.020033 0.766733 7 N 0.338485 0.338483 0.338480 -0.027538 -0.027547 -0.027546 8 B -0.017533 -0.017536 -0.017538 0.417338 0.417339 0.417334 7 8 1 H 0.338485 -0.017533 2 H 0.338483 -0.017536 3 H 0.338480 -0.017538 4 H -0.027538 0.417338 5 H -0.027547 0.417339 6 H -0.027546 0.417334 7 N 6.475920 0.182855 8 B 0.182855 3.582085 Mulliken charges: 1 1 H 0.302276 2 H 0.302279 3 H 0.302277 4 H -0.116960 5 H -0.116964 6 H -0.116971 7 N -0.591593 8 B 0.035655 Sum of Mulliken charges = 0.00000 Mulliken charges with hydrogens summed into heavy atoms: 1 7 N 0.315240 8 B -0.315240 Electronic spatial extent (au): = 117.9549 Charge= 0.0000 electrons Dipole moment (field-independent basis, Debye): X= -5.5653 Y= 0.0000 Z= -0.0002 Tot= 5.5653 Quadrupole moment (field-independent basis, Debye-Ang): XX= -16.1085 YY= -15.5751 ZZ= -15.5753 XY= 0.0001 XZ= -0.0001 YZ= -0.0001 Traceless Quadrupole moment (field-independent basis, Debye-Ang): XX= -0.3555 YY= 0.1778 ZZ= 0.1776 XY= 0.0001 XZ= -0.0001 YZ= -0.0001 Octapole moment (field-independent basis, Debye-Ang**2): XXX= -18.3945 YYY= 1.5787 ZZZ= 0.2061 XYY= -8.1091 XXY= -0.0001 XXZ= 0.0002 XZZ= -8.1090 YZZ= -1.5784 YYZ= -0.2067 XYZ= 0.0001 Hexadecapole moment (field-independent basis, Debye-Ang**3): XXXX= -106.7251 YYYY= -34.2975 ZZZZ= -34.2972 XXXY= 0.0002 XXXZ= -0.0004 YYYX= -0.7776 YYYZ= -0.0001 ZZZX= -0.1023 ZZZY= 0.0000 XXYY= -23.5239 XXZZ= -23.5238 YYZZ= -11.4324 XXYZ= -0.0003 YYXZ= 0.1022 ZZXY= 0.7777 N-N= 4.043471937775D+01 E-N=-2.729558806190D+02 KE= 8.236626618078D+01 1|1| IMPERIAL COLLEGE-CHWS-271|FOpt|RB3LYP|6-31G(d,p)|B1H6N1|AK7611|13 -Oct-2013|0||# opt b3lyp/6-31g(d,p) geom=connectivity||NH3BH3 Optimisa tion||0,1|H,-2.9842356289,1.2294830826,-0.0001433227|H,-2.9840974538,- 0.1967978196,-0.8233379345|H,-2.9841065356,-0.1965172253,0.8234427552| H,-0.6454567181,-0.8921722002,-0.0000835123|H,-0.6456035598,0.86437447 37,-1.0140859666|H,-0.6456214954,0.864257318,1.0141695614|N,-2.6186143 988,0.2787583449,-0.0000244041|B,-0.9505635295,0.2788111859,0.00000682 37||Version=EM64W-G09RevD.01|State=1-A|HF=-83.2246901|RMSD=1.352e-009| RMSF=6.839e-005|Dipole=-2.1895743,-0.0001108,-0.0000143|Quadrupole=-0. 2642922,0.1321413,0.1321509,-0.000084,-0.0000234,0.0000875|PG=C01 [X(B 1H6N1)]||@ WOE UNTO THEM THAT CALL EVIL GOOD, AND GOOD EVIL Job cpu time: 0 days 0 hours 0 minutes 29.0 seconds. File lengths (MBytes): RWF= 5 Int= 0 D2E= 0 Chk= 1 Scr= 1 Normal termination of Gaussian 09 at Sun Oct 13 16:53:48 2013.