Wall clock time and date at job start Mon Mar 8 2021 12:09:38 AMSOL-version 7.1 by G. D. Hawkins, D. J. Giesen, G. C. Lynch, C. C. Chambers, I. Rossi, J. W. Storer, J. Li, J. D. Thompson, P. Winget, B. J. Lynch, D. Rinaldi, D. A. Liotard, C. J. Cramer, and D. G. Truhlar Copyright 2004 by Regents of the University of Minnesota. All rights reserved. Notice: recipients of this code are asked to comply with the user agreement in Section 1 of the documentation file. ******************************************************************************* * 1SCF - SCF CALCULATION WITHOUT GEOMETRY OPTIMIZATION * - USE EF ROUTINE FOR MINIMUM SEARCH (DEFAULT) * GEO-OK - OVERRIDE INTERATOMIC DISTANCE CHECK * TLIMIT= - A TIME OF 15. SECONDS REQUESTED * CHARGE - CHARGE ON SYSTEM= 0 * AM1 - THE AM1 HAMILTONIAN TO BE USED * DEV - DEVELOPER OPTIONS ARE ALLOWED * SM5.42R - SM5.42R CALCULATIONS WILL BE PERFORMED * - CHARGE MODEL 2 WILL BE USED (CM2A) * - THE SOLVENT IS USER-SPECIFIED * DIELEC - THE SOLVENT DIELECTRIC CONSTANT IS 2.06 * IOFR - THE SOLVENT INDEX OF REFRACTION IS 1.4345 * ALPHA - THE SOLVENT ALPHA IS 0.00 * BETA - THE SOLVENT BETA IS 0.00 * GAMMA - THE SOLVENT MACROSCOPIC SURFACE TENSION IS * 38.93 CAL MOL^-1 ANGSTROM^-2 * FACARB - THE FRACTION OF AROMATIC CARBONS IS 0.00 * FEHALO - THE FRACTION OF ELECTRONEGATIVE HALOGENS IS 0.00 ******************************************************************************* Atom NO. Chemical Bond length Bond angle Dihedral angle number (k) symbol (angstroms) (degrees) (degrees) (I) NA:I NB:NA:I NC:NB:NA:I NA NB NC 1 1 C 2 2 C 1.53002 * 1 3 3 O 1.45196 * 109.47293 * 2 1 4 4 C 1.34227 * 117.00335 * 180.02562 * 3 2 1 5 5 O 1.20830 * 119.99627 * 0.02562 * 4 3 2 6 6 C 1.50697 * 120.00352 * 179.97438 * 4 3 2 7 7 H 1.08999 * 112.84741 * 311.20063 * 6 4 3 8 8 C 1.53779 * 113.61006 * 82.50323 * 6 4 3 9 9 C 1.53788 * 87.08111 * 139.98209 * 8 6 4 10 10 F 1.39898 * 113.69375 * 220.00803 * 9 8 6 11 11 C 1.53004 * 113.61388 * 89.11233 * 9 8 6 12 12 S 1.81398 * 109.46828 * 82.50775 * 11 9 8 13 13 O 1.42098 * 110.54549 * 63.05084 * 12 11 9 14 14 O 1.42100 * 110.54085 * 286.29950 * 12 11 9 15 15 O 1.52204 * 104.44811 * 174.67735 * 12 11 9 16 16 C 1.35833 * 113.99794 * 79.01932 * 15 12 11 17 17 C 1.40951 * 120.21312 * 272.02227 * 16 15 12 18 18 C 1.40418 * 121.28821 * 0.02562 * 17 16 15 19 19 C 1.36131 * 119.66332 * 179.72333 * 18 17 16 20 20 C 1.39367 * 120.97800 * 0.02562 * 19 18 17 21 21 C 1.36114 * 120.97209 * 359.97438 * 20 19 18 22 22 C 1.40527 * 119.65093 * 0.03926 * 21 20 19 23 23 C 1.40613 * 121.26912 * 180.02562 * 22 21 20 24 24 Cl 1.73601 * 120.12953 * 359.97438 * 23 22 21 25 25 C 1.36269 * 119.74491 * 179.97438 * 23 22 21 26 26 C 1.36647 * 120.21774 * 92.26178 * 16 15 12 27 27 C 1.50703 * 119.53169 * 0.02562 * 26 16 15 28 28 O 1.42898 * 109.47040 * 283.66722 * 27 26 16 29 29 Si 1.86301 * 109.47178 * 163.66929 * 27 26 16 30 30 C 1.53780 * 113.61418 * 179.97438 * 6 4 3 31 31 H 1.09004 * 109.46638 * 299.99641 * 1 2 3 32 32 H 1.08998 * 109.47414 * 59.99801 * 1 2 3 33 33 H 1.08997 * 109.47237 * 180.02562 * 1 2 3 34 34 H 1.09004 * 109.46638 * 240.00359 * 2 1 3 35 35 H 1.09004 * 109.47078 * 120.00492 * 2 1 3 36 36 H 1.09001 * 113.61905 * 25.43841 * 8 6 4 37 37 H 1.09004 * 113.61930 * 254.52623 * 8 6 4 38 38 H 1.08994 * 109.46938 * 322.50849 * 11 9 8 39 39 H 1.08995 * 109.46736 * 202.50681 * 11 9 8 40 40 H 1.08001 * 120.17310 * 359.69747 * 18 17 16 41 41 H 1.07999 * 119.51399 * 179.97438 * 19 18 17 42 42 H 1.07998 * 119.51475 * 179.97438 * 20 19 18 43 43 H 1.07995 * 120.17147 * 180.02562 * 21 20 19 44 44 H 1.08000 * 119.48411 * 179.97438 * 25 23 22 45 45 H 1.08995 * 109.46938 * 43.66978 * 27 26 16 46 46 H 0.96701 * 114.00094 * 300.00215 * 28 27 26 47 47 H 1.48504 * 110.00092 * 58.68076 * 29 27 26 48 48 H 1.48499 * 109.99947 * 180.02562 * 29 27 26 49 49 H 1.09008 * 113.60719 * 105.47468 * 30 6 4 50 50 H 1.08992 * 113.69413 * 334.58618 * 30 6 4 Note: An asterisk (*) indicates this parameter was or will be optimized. Cartesian coordinates (angstroms) NO. Atom X Y Z 1 6 0.0000 0.0000 0.0000 2 6 1.5300 0.0000 0.0000 3 8 2.0140 1.3689 0.0000 4 6 3.3447 1.5448 -0.0005 5 8 4.0808 0.5865 -0.0015 6 6 3.9208 2.9373 -0.0012 7 1 3.4709 3.5811 -0.7570 8 6 3.9868 3.5772 1.3957 9 6 5.3569 4.1426 0.9856 10 9 6.3115 4.0997 2.0075 11 6 5.2766 5.5104 0.3047 12 16 5.1611 6.8083 1.5667 13 8 6.3324 6.8274 2.3710 14 8 3.8679 6.8174 2.1556 15 8 5.2246 8.1097 0.7801 16 6 4.0521 8.4546 0.1872 17 6 3.7438 7.9806 -1.1039 18 6 4.6297 7.1474 -1.8057 19 6 4.2959 6.7042 -3.0488 20 6 3.0858 7.0654 -3.6382 21 6 2.2028 7.8714 -2.9875 22 6 2.5134 8.3485 -1.7028 23 6 1.6260 9.1825 -0.9997 24 17 0.1128 9.6441 -1.7143 25 6 1.9585 9.6269 0.2449 26 6 3.1661 9.2673 0.8366 27 6 3.4992 9.7784 2.2147 28 8 3.3060 8.7333 3.1699 29 14 2.3864 11.2135 2.6306 30 6 5.4583 2.9680 -0.0011 31 1 -0.3633 0.5138 0.8901 32 1 -0.3634 0.5138 -0.8899 33 1 -0.3633 -1.0276 0.0005 34 1 1.8933 -0.5138 -0.8901 35 1 1.8934 -0.5139 0.8900 36 1 4.0393 2.8508 2.2067 37 1 3.2299 4.3439 1.5611 38 1 4.3950 5.5464 -0.3352 39 1 6.1703 5.6689 -0.2987 40 1 5.5714 6.8600 -1.3618 41 1 4.9793 6.0638 -3.5866 42 1 2.8451 6.7005 -4.6258 43 1 1.2685 8.1428 -3.4563 44 1 1.2740 10.2676 0.7810 45 1 4.5386 10.1055 2.2393 46 1 2.4015 8.3928 3.2028 47 1 0.9649 10.7883 2.5692 48 1 2.7054 11.7270 3.9870 49 1 5.8895 3.2427 -0.9638 50 1 5.9141 2.0756 0.4277 There are 78 doubly filled levels RHF calculation, no. of doubly occupied orbitals= 78 No. of singly occupied orbitals= 1 REFERENCES FOR PARAMETERS IN GAS-PHASE HAMILTONIAN: H: (AM1): M.J.S. DEWAR ET AL, J. AM. CHEM. SOC. 107 3902-3909 (1985) C: (AM1): M.J.S. DEWAR ET AL, J. AM. CHEM. SOC. 107 3902-3909 (1985) O: (AM1): M.J.S. DEWAR ET AL, J. AM. CHEM. SOC. 107 3902-3909 (1985) F: (AM1): M.J.S. DEWAR AND E. G. ZOEBISCH, THEOCHEM, 180, 1 (1988). Si: (AM1): M.J.S.DEWAR, C. JIE, ORGANOMETALLICS, 6, 1486-1490 (1987). S: (AM1): M.J.S. DEWAR, Y.-C. YUAN, INORGANIC CHEM., 29, 589 (1990) Cl: (AM1): M.J.S. DEWAR AND E. G. ZOEBISCH, THEOCHEM, 180, 1 (1988). REFERENCE FOR CHARGE MODEL 2: J. Li, J. Xing, C. J. Cramer, and D. G. Truhlar, J. Chem. Phys. 111 (1999) 885. REFERENCE FOR THE PARAMETERIZATION OF THE SM5.42R SOLVATION MODEL: J. Li, T. Zhu, G. D. Hawkins, P. Winget, D. A. Liotard, C. J. Cramer, and D. G. Truhlar, Theor. Chem. Acc. 103 (1999) 9-63 P. Winget, J. D. Thompson C. J. Cramer,and D. G. Truhlar, J. Phys. Chem. B. submitted. CHARGE=0 AM1 1SCF TLIMIT=15 GEO-OK SM5.42R & SOLVNT=GENORG IOFR=1.4345 ALPHA=0.00 BETA=0.00 GAMMA=38.93 & DIELEC=2.06 FACARB=0.00 FEHALO=0.00 DEV s_40_8272896_15494190.mol2 50 1SCF WAS SPECIFIED, SO GEOMETRY OPTIMIZATION WAS NOT USED SCF FIELD WAS ACHIEVED Wall clock time and date at job start Mon Mar 8 2021 12:09:38 Heat of formation + Delta-G solvation = -172.009976 kcal Electronic energy + Delta-G solvation = -44796.820345 eV Core-core repulsion = 39019.322442 eV Total energy + Delta-G solvation = -5777.497903 eV No. of doubly occupied orbitals = 78 and no. of open levels = 1 Molecular weight (most abundant/longest-lived isotopes) = 459.078 amu Computer time = 3.28 seconds Orbital eigenvalues (eV) -49.01717 -43.03672 -41.91069 -40.65690 -39.68034 -38.92321 -38.23765 -37.98632 -37.54145 -37.12945 -36.77749 -33.74192 -32.84651 -32.48150 -31.64504 -30.57002 -27.46595 -27.40172 -27.12940 -25.23560 -24.15259 -24.06227 -23.54694 -22.50776 -21.74021 -20.53923 -20.03402 -19.53191 -19.19214 -18.88355 -18.81346 -18.13600 -17.92148 -17.83223 -17.57793 -17.11729 -17.04359 -16.75231 -16.33895 -16.21588 -16.04832 -16.00770 -15.68288 -15.50235 -15.42438 -15.22315 -15.13471 -14.92373 -14.60962 -14.43248 -14.26693 -14.09292 -13.86617 -13.83754 -13.81514 -13.69066 -13.46337 -13.21727 -13.12405 -13.11496 -13.02239 -12.88518 -12.71857 -12.65851 -12.61183 -12.55779 -12.42430 -12.23574 -12.10421 -11.96751 -11.82549 -11.60290 -11.53293 -11.47027 -11.09009 -10.62168 -10.17933 -9.29908 -5.05515 -1.48560 -0.84123 -0.62355 0.35637 0.64547 0.90373 1.01458 1.11134 1.33235 1.46919 1.48004 1.88053 2.30411 2.40814 2.48954 2.68607 2.79574 2.84042 2.95193 3.01420 3.11335 3.18406 3.29526 3.31459 3.43959 3.48528 3.63461 3.69774 3.71254 3.80055 3.88693 3.91542 3.92573 3.93683 4.01884 4.07977 4.22364 4.26655 4.30628 4.34977 4.35850 4.40324 4.40513 4.52462 4.54036 4.61956 4.65656 4.73882 4.74296 4.87821 4.95496 4.98356 5.08095 5.13072 5.35809 5.40782 6.50572 6.61966 Net atomic charges, atomic populations, and dipole contributions calculated with CM2 Atom NO. Type Charge No. of electrons 1 C -0.154 4.154 2 C 0.060 3.940 3 O -0.354 6.354 4 C 0.475 3.525 5 O -0.499 6.499 6 C -0.107 4.107 7 H 0.134 0.866 8 C -0.148 4.148 9 C 0.104 3.896 10 F -0.178 7.178 11 C -0.684 4.684 12 S 2.593 3.407 13 O -0.900 6.900 14 O -0.912 6.912 15 O -0.653 6.653 16 C 0.099 3.901 17 C -0.050 4.050 18 C -0.139 4.139 19 C -0.086 4.086 20 C -0.104 4.104 21 C -0.091 4.091 22 C -0.032 4.032 23 C 0.004 3.996 24 Cl -0.036 7.036 25 C -0.106 4.106 26 C -0.058 4.058 27 C -0.099 4.099 28 O -0.498 6.498 29 Si 0.693 3.307 30 C -0.121 4.121 31 H 0.068 0.932 32 H 0.068 0.932 33 H 0.082 0.918 34 H 0.072 0.928 35 H 0.072 0.928 36 H 0.116 0.884 37 H 0.110 0.890 38 H 0.159 0.841 39 H 0.157 0.843 40 H 0.137 0.863 41 H 0.137 0.863 42 H 0.137 0.863 43 H 0.144 0.856 44 H 0.156 0.844 45 H 0.126 0.874 46 H 0.384 0.616 47 H -0.251 1.251 48 H -0.247 1.247 49 H 0.099 0.901 50 H 0.121 0.879 For the charges calculated by CM2: Dipole moment (debyes) X Y Z Total from point charges -5.021 1.489 -7.389 9.057 Note: The Mulliken population analysis charges presented below were not used in the solvation calculation but are provided for completeness. The chosen solvation model uses CM2 partial charges for calculating solvation energies. Net atomic charges, atomic populations, and dipole contributions using Mulliken population analysis Atom NO. Type Charge No. of electrons 1 C -0.211 4.211 2 C -0.013 4.013 3 O -0.271 6.271 4 C 0.317 3.683 5 O -0.386 6.386 6 C -0.126 4.126 7 H 0.152 0.848 8 C -0.185 4.185 9 C 0.085 3.915 10 F -0.158 7.158 11 C -0.807 4.807 12 S 2.685 3.315 13 O -0.899 6.899 14 O -0.910 6.910 15 O -0.611 6.611 16 C 0.057 3.943 17 C -0.052 4.052 18 C -0.157 4.157 19 C -0.104 4.104 20 C -0.121 4.121 21 C -0.109 4.109 22 C -0.033 4.033 23 C -0.024 4.024 24 Cl -0.007 7.007 25 C -0.124 4.124 26 C -0.062 4.062 27 C -0.209 4.209 28 O -0.306 6.306 29 Si 0.593 3.407 30 C -0.159 4.159 31 H 0.087 0.913 32 H 0.087 0.913 33 H 0.101 0.899 34 H 0.090 0.910 35 H 0.090 0.910 36 H 0.134 0.866 37 H 0.128 0.872 38 H 0.177 0.823 39 H 0.175 0.825 40 H 0.155 0.845 41 H 0.155 0.845 42 H 0.155 0.845 43 H 0.161 0.839 44 H 0.173 0.827 45 H 0.144 0.856 46 H 0.233 0.767 47 H -0.176 1.176 48 H -0.171 1.171 49 H 0.117 0.883 50 H 0.139 0.861 Dipole moment (debyes) X Y Z Total from point charges -4.753 2.008 -6.154 8.031 hybrid contribution -0.078 -1.702 -0.386 1.747 sum -4.831 0.305 -6.541 8.137 Atomic orbital electron populations 1.21828 0.92190 1.03355 1.03728 1.23131 0.95713 0.79580 1.02891 1.86507 1.23970 1.33008 1.83574 1.23131 0.80917 0.91757 0.72482 1.90708 1.60525 1.42652 1.44750 1.22261 0.97270 0.92186 1.00855 0.84831 1.23844 0.95949 1.02059 0.96667 1.21509 0.88909 0.93354 0.87760 1.93020 1.65596 1.97193 1.59973 1.31859 1.13753 1.14008 1.21033 1.05859 0.75795 0.74661 0.75202 1.93693 1.49856 1.83905 1.62410 1.93660 1.40849 1.83210 1.73295 1.85746 1.41976 1.59725 1.73690 1.17306 0.82740 1.02497 0.91794 1.17874 0.94979 0.99565 0.92774 1.21282 0.99516 0.98520 0.96389 1.21059 0.96770 0.98253 0.94322 1.21082 0.93522 0.96938 1.00607 1.21227 1.00523 0.96508 0.92615 1.17563 0.92406 0.98026 0.95340 1.20758 0.89196 0.98265 0.94204 1.98423 1.28583 1.90500 1.83194 1.20706 0.96377 1.00978 0.94340 1.20086 0.92629 0.95910 0.97602 1.25289 1.10802 0.99518 0.85245 1.86720 1.35454 1.51086 1.57328 0.95051 0.77600 0.81787 0.86269 1.23604 0.93487 0.98185 1.00635 0.91280 0.91335 0.89911 0.90970 0.90951 0.86580 0.87184 0.82323 0.82483 0.84509 0.84468 0.84489 0.83876 0.82692 0.85578 0.76655 1.17590 1.17148 0.88255 0.86097 Number of geometries 1 Number of calculations of the screened coulomb radii 1 The total number of SCF iterations 26. Note: The number of geometries may not correspond to the number of cycles due to rejected geometry changes. In the following table subtotal= G_P + SS G_CDS. Atom Chem. CM2 G_P Area Sigma k SS G_CDS Subtotal M number symbol chg. (kcal) (Ang**2) cal/(Ang**2) (kcal) (kcal) value 1 C -0.15 0.32 10.29 37.16 0.38 0.70 16 2 C 0.06 -0.01 6.19 37.16 0.23 0.22 16 3 O -0.35 -0.35 10.22 -39.25 -0.40 -0.75 16 4 C 0.48 0.61 7.48 36.00 0.27 0.88 16 5 O -0.50 -2.04 15.65 -18.75 -0.29 -2.34 16 6 C -0.11 0.16 3.91 -90.96 -0.36 -0.20 16 7 H 0.13 -0.55 7.66 -51.93 -0.40 -0.94 16 8 C -0.15 0.11 7.09 -25.92 -0.18 -0.08 16 9 C 0.10 0.00 1.79 -89.82 -0.16 -0.16 16 10 F -0.18 -0.61 16.09 2.25 0.04 -0.57 16 11 C -0.68 1.06 3.60 37.16 0.13 1.19 16 12 S 2.59 6.43 4.89 -107.50 -0.53 5.91 16 13 O -0.90 -6.04 17.25 -67.04 -1.16 -7.19 16 14 O -0.91 -2.61 12.30 -70.84 -0.87 -3.48 16 15 O -0.65 -1.16 8.64 -51.54 -0.45 -1.61 16 16 C 0.10 0.08 4.66 -39.13 -0.18 -0.10 16 17 C -0.05 0.03 5.44 -106.73 -0.58 -0.55 16 18 C -0.14 0.35 6.76 -39.52 -0.27 0.08 16 19 C -0.09 0.33 9.95 -39.92 -0.40 -0.07 16 20 C -0.10 0.38 9.95 -39.92 -0.40 -0.01 16 21 C -0.09 0.23 9.44 -39.48 -0.37 -0.14 16 22 C -0.03 0.03 5.74 -106.82 -0.61 -0.58 16 23 C 0.00 0.00 6.32 -39.39 -0.25 -0.25 16 24 Cl -0.04 0.06 27.50 -51.86 -1.43 -1.37 16 25 C -0.11 0.11 8.33 -39.93 -0.33 -0.22 16 26 C -0.06 -0.01 4.82 -104.84 -0.51 -0.52 16 27 C -0.10 -0.06 2.28 36.01 0.08 0.03 16 28 O -0.50 -1.20 9.42 -59.74 -0.56 -1.77 16 29 Si 0.69 -0.94 30.27 -169.99 -5.15 -6.09 16 30 C -0.12 0.15 6.88 -25.94 -0.18 -0.03 16 31 H 0.07 -0.13 8.14 -51.93 -0.42 -0.55 16 32 H 0.07 -0.15 8.14 -51.93 -0.42 -0.57 16 33 H 0.08 -0.27 8.14 -51.93 -0.42 -0.69 16 34 H 0.07 -0.03 8.13 -51.93 -0.42 -0.45 16 35 H 0.07 -0.01 8.14 -51.93 -0.42 -0.43 16 36 H 0.12 -0.07 8.09 -51.93 -0.42 -0.49 16 37 H 0.11 -0.19 7.69 -51.93 -0.40 -0.59 16 38 H 0.16 -0.63 5.02 -51.93 -0.26 -0.89 16 39 H 0.16 -0.66 6.73 -51.93 -0.35 -1.01 16 40 H 0.14 -0.55 4.69 -52.49 -0.25 -0.80 16 41 H 0.14 -0.84 8.06 -52.49 -0.42 -1.26 16 42 H 0.14 -0.78 8.06 -52.49 -0.42 -1.20 16 43 H 0.14 -0.52 6.98 -52.49 -0.37 -0.89 16 44 H 0.16 -0.40 4.78 -52.49 -0.25 -0.65 16 45 H 0.13 0.04 7.74 -51.93 -0.40 -0.36 16 46 H 0.38 -0.07 8.49 45.56 0.39 0.32 16 47 H -0.25 0.15 6.61 56.52 0.37 0.53 16 48 H -0.25 -0.45 7.11 56.52 0.40 -0.05 16 49 H 0.10 -0.35 8.10 -51.92 -0.42 -0.77 16 50 H 0.12 0.02 7.44 -51.93 -0.39 -0.37 16 LS Contribution 427.09 15.07 6.44 6.44 Total: 0.00 -11.01 427.09 -13.73 -24.74 By element: Atomic # 1 Polarization: -6.42 SS G_CDS: -5.70 Total: -12.11 kcal Atomic # 6 Polarization: 3.87 SS G_CDS: -3.68 Total: 0.19 kcal Atomic # 8 Polarization: -13.40 SS G_CDS: -3.73 Total: -17.13 kcal Atomic # 9 Polarization: -0.61 SS G_CDS: 0.04 Total: -0.57 kcal Atomic # 14 Polarization: -0.94 SS G_CDS: -5.15 Total: -6.09 kcal Atomic # 16 Polarization: 6.43 SS G_CDS: -0.53 Total: 5.91 kcal Atomic # 17 Polarization: 0.06 SS G_CDS: -1.43 Total: -1.37 kcal Total LS contribution 6.44 Total: 6.44 kcal Total: -11.01 -13.73 -24.74 kcal The number of atoms in the molecule is 50 The average number of expansion shells was 16.00 The maximum number of expansion shells was 16 The minimum number of expansion shells was 16 **** NOTA BENE **** This is the net solvation energy for this exact molecular structure (nuclear and electronic)! The standard-state solvation energy should be obtained as the difference between the heat of formation plus delta-G solvation for the relaxed solvated system and that for the relaxed gas-phase system. s_40_8272896_15494190.mol2 50 1SCF run This is a breakdown of the solvation energy calculated without geometric relaxation in solution: (1) E-EN(sol) electronic-nuclear energy of solute -147.270 kcal (2) G-P(sol) polarization free energy of solvation -11.008 kcal (3) G-ENP(sol) elect.-nuc.-pol. free energy of system -158.279 kcal (4) G-CDS(sol) cavity-dispersion-solvent structure free energy -13.731 kcal (5) G-P-CDS(sol) = G-P(sol) + G-CDS(sol) = (2) + (4) -24.740 kcal (6) G-S(sol) free energy of system = (1) + (5) -172.010 kcal 1SCF WAS SPECIFIED, SO GEOMETRY OPTIMIZATION WAS NOT USED Total computer time = 3.28 seconds