Wall clock time and date at job start Tue Mar 9 2021 12:54: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 O 2 2 C 1.42899 * 1 3 3 Si 1.86298 * 109.47406 * 2 1 4 4 C 1.50702 * 109.46954 * 240.00547 * 2 1 3 5 5 C 1.35083 * 126.08095 * 105.17886 * 4 2 1 6 6 N 1.35035 * 107.98169 * 179.87227 * 5 4 2 7 7 C 1.46508 * 125.93181 * 180.11408 * 6 5 4 8 8 C 1.50694 * 109.47105 * 125.30878 * 7 6 5 9 9 O 1.21282 * 120.00066 * 359.97438 * 8 7 6 10 10 N 1.34775 * 119.99704 * 180.02562 * 8 7 6 11 11 C 1.46506 * 119.99635 * 179.97438 * 10 8 7 12 12 H 1.08997 * 109.47243 * 324.99635 * 11 10 8 13 13 C 1.52994 * 109.46949 * 85.00039 * 11 10 8 14 14 C 1.53006 * 109.47144 * 180.02562 * 13 11 10 15 15 C 1.53003 * 109.46892 * 299.99891 * 14 13 11 16 16 H 1.09005 * 109.47005 * 300.00111 * 15 14 13 17 17 C 1.52994 * 109.47139 * 179.97438 * 15 14 13 18 18 O 1.42904 * 109.47012 * 65.00405 * 17 15 14 19 19 C 1.53005 * 109.47089 * 59.99736 * 15 14 13 20 20 C 1.52994 * 109.47204 * 300.00051 * 19 15 14 21 21 N 1.40075 * 108.13896 * 0.40005 * 6 5 4 22 22 C 1.31616 * 108.10544 * 359.75018 * 21 6 5 23 23 C 1.48418 * 126.03383 * 179.97438 * 22 21 6 24 24 C 1.39332 * 120.11593 * 137.23427 * 23 22 21 25 25 C 1.38028 * 119.87702 * 179.74297 * 24 23 22 26 26 C 1.38321 * 120.11657 * 0.54040 * 25 24 23 27 27 C 1.38320 * 120.24367 * 359.73485 * 26 25 24 28 28 C 1.38016 * 120.12042 * 359.97438 * 27 26 25 29 29 H 0.96703 * 114.00123 * 59.99340 * 1 2 3 30 30 H 1.08991 * 109.47450 * 120.00478 * 2 1 3 31 31 H 1.48499 * 110.00274 * 298.68769 * 3 2 1 32 32 H 1.48505 * 109.99890 * 60.00278 * 3 2 1 33 33 H 1.08002 * 126.00694 * 359.97438 * 5 4 2 34 34 H 1.08996 * 109.46881 * 245.30898 * 7 6 5 35 35 H 1.09008 * 109.46750 * 5.30483 * 7 6 5 36 36 H 0.96993 * 120.00207 * 359.97438 * 10 8 7 37 37 H 1.08998 * 109.47198 * 299.99328 * 13 11 10 38 38 H 1.08997 * 109.47478 * 59.99740 * 13 11 10 39 39 H 1.09001 * 109.47499 * 179.97438 * 14 13 11 40 40 H 1.08997 * 109.47241 * 59.99478 * 14 13 11 41 41 H 1.09006 * 109.47014 * 184.99581 * 17 15 14 42 42 H 1.08997 * 109.47520 * 305.00058 * 17 15 14 43 43 H 0.96702 * 113.99895 * 180.02562 * 18 17 15 44 44 H 1.09009 * 109.46833 * 180.02562 * 19 15 14 45 45 H 1.08997 * 109.47201 * 60.00317 * 19 15 14 46 46 H 1.09005 * 109.47405 * 299.99940 * 20 19 15 47 47 H 1.09001 * 109.47233 * 179.97438 * 20 19 15 48 48 H 1.08000 * 120.06231 * 0.02562 * 24 23 22 49 49 H 1.08000 * 119.94016 * 180.29859 * 25 24 23 50 50 H 1.08004 * 119.87351 * 179.72107 * 26 25 24 51 51 H 1.07997 * 119.93733 * 180.02562 * 27 26 25 52 52 H 1.08005 * 120.06319 * 180.02562 * 28 27 26 Note: An asterisk (*) indicates this parameter was or will be optimized. Cartesian coordinates (angstroms) NO. Atom X Y Z 1 8 0.0000 0.0000 0.0000 2 6 1.4290 0.0000 0.0000 3 14 2.0501 1.7564 0.0000 4 6 1.9313 -0.7103 -1.2306 5 6 2.4660 -0.1251 -2.3244 6 7 2.7928 -1.0992 -3.2006 7 6 3.4035 -0.9015 -4.5175 8 6 4.6642 -1.7214 -4.6148 9 8 5.0180 -2.4033 -3.6764 10 7 5.3990 -1.6960 -5.7444 11 6 6.6249 -2.4926 -5.8387 12 1 7.1166 -2.5193 -4.8663 13 6 6.2750 -3.9174 -6.2726 14 6 7.5553 -4.7496 -6.3706 15 6 8.4956 -4.1193 -7.3999 16 1 8.0043 -4.0926 -8.3726 17 6 9.7760 -4.9511 -7.4974 18 8 9.4643 -6.2466 -8.0136 19 6 8.8454 -2.6944 -6.9660 20 6 7.5653 -1.8624 -6.8680 21 7 2.4406 -2.3311 -2.6345 22 6 1.9210 -2.1020 -1.4471 23 6 1.4132 -3.1288 -0.5034 24 6 1.7059 -3.0351 0.8556 25 6 1.2267 -3.9898 1.7296 26 6 0.4666 -5.0444 1.2569 27 6 0.1778 -5.1454 -0.0920 28 6 0.6469 -4.1930 -0.9739 29 1 -0.3933 0.4418 -0.7650 30 1 1.7924 -0.5139 0.8898 31 1 1.5879 2.4586 -1.2241 32 1 1.5616 2.4678 1.2086 33 1 2.6086 0.9353 -2.4715 34 1 2.7049 -1.2159 -5.2928 35 1 3.6446 0.1531 -4.6516 36 1 5.1162 -1.1503 -6.4947 37 1 5.6055 -4.3665 -5.5391 38 1 5.7833 -3.8907 -7.2450 39 1 7.3060 -5.7648 -6.6793 40 1 8.0466 -4.7762 -5.3980 41 1 10.4818 -4.4548 -8.1636 42 1 10.2203 -5.0517 -6.5072 43 1 10.2308 -6.8296 -8.1022 44 1 9.5151 -2.2453 -7.6996 45 1 9.3368 -2.7211 -5.9934 46 1 7.0739 -1.8357 -7.8407 47 1 7.8144 -0.8474 -6.5585 48 1 2.3028 -2.2143 1.2251 49 1 1.4489 -3.9155 2.7839 50 1 0.0974 -5.7917 1.9437 51 1 -0.4165 -5.9706 -0.4556 52 1 0.4205 -4.2726 -2.0270 There are 69 doubly filled levels RHF calculation, no. of doubly occupied orbitals= 69 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) N: (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) Si: (AM1): M.J.S.DEWAR, C. JIE, ORGANOMETALLICS, 6, 1486-1490 (1987). 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_22_7207542_12622748.mol2 52 1SCF WAS SPECIFIED, SO GEOMETRY OPTIMIZATION WAS NOT USED SCF FIELD WAS ACHIEVED Wall clock time and date at job start Tue Mar 9 2021 12:54:38 Heat of formation + Delta-G solvation = -14.722728 kcal Electronic energy + Delta-G solvation = -33622.879843 eV Core-core repulsion = 29169.223694 eV Total energy + Delta-G solvation = -4453.656148 eV No. of doubly occupied orbitals = 69 and no. of open levels = 1 Molecular weight (most abundant/longest-lived isotopes) = 372.188 amu Computer time = 1.07 seconds Orbital eigenvalues (eV) -42.42469 -41.08754 -39.73646 -39.44299 -37.75274 -37.45282 -35.24996 -34.72634 -32.34031 -31.96387 -31.85194 -31.44285 -31.08033 -28.67911 -26.38418 -25.73583 -24.04313 -23.83949 -23.26654 -22.38421 -21.66896 -21.50259 -21.05191 -20.71777 -18.81117 -18.47299 -18.12979 -17.60682 -17.49392 -16.97854 -16.64950 -16.29289 -16.12483 -15.99990 -15.76225 -15.53914 -15.48142 -15.43924 -15.12549 -14.92380 -14.61229 -14.33965 -14.15486 -14.10025 -14.00834 -13.77765 -13.71569 -13.46186 -13.26985 -13.15406 -13.06408 -12.47799 -12.30872 -12.24932 -12.21279 -11.99355 -11.74806 -11.67674 -11.54877 -11.50144 -11.33791 -11.08412 -10.89620 -10.54728 -10.25568 -10.11877 -9.66272 -9.65806 -9.16252 -4.93294 0.31800 0.68270 0.84091 1.22873 1.36845 1.48791 1.57168 1.84020 1.99797 2.03081 3.08167 3.14406 3.30455 3.46913 3.53898 3.56272 3.65099 3.70486 3.75738 3.85224 3.88536 3.92278 4.13490 4.17980 4.20310 4.25069 4.28668 4.36451 4.36980 4.39795 4.46829 4.49886 4.52246 4.56425 4.65601 4.71252 4.72893 4.85731 4.88354 4.94085 4.95089 5.01464 5.07614 5.14574 5.24031 5.29014 5.32169 5.38173 5.44197 5.46377 5.74549 5.77668 5.79744 5.86228 6.21321 6.33910 6.57000 6.65651 6.72659 7.10276 Net atomic charges, atomic populations, and dipole contributions calculated with CM2 Atom NO. Type Charge No. of electrons 1 O -0.526 6.526 2 C -0.035 4.035 3 Si 0.680 3.320 4 C -0.247 4.247 5 C 0.060 3.940 6 N -0.358 5.358 7 C 0.115 3.885 8 C 0.499 3.501 9 O -0.501 6.501 10 N -0.727 5.727 11 C 0.147 3.853 12 H 0.089 0.911 13 C -0.129 4.129 14 C -0.108 4.108 15 C -0.089 4.089 16 H 0.079 0.921 17 C 0.078 3.922 18 O -0.573 6.573 19 C -0.110 4.110 20 C -0.124 4.124 21 N -0.247 5.247 22 C 0.132 3.868 23 C -0.007 4.007 24 C -0.092 4.092 25 C -0.120 4.120 26 C -0.113 4.113 27 C -0.121 4.121 28 C -0.092 4.092 29 H 0.371 0.629 30 H 0.115 0.885 31 H -0.260 1.260 32 H -0.255 1.255 33 H 0.178 0.822 34 H 0.117 0.883 35 H 0.128 0.872 36 H 0.400 0.600 37 H 0.069 0.931 38 H 0.064 0.936 39 H 0.074 0.926 40 H 0.063 0.937 41 H 0.046 0.954 42 H 0.046 0.954 43 H 0.381 0.619 44 H 0.071 0.929 45 H 0.069 0.931 46 H 0.070 0.930 47 H 0.072 0.928 48 H 0.125 0.875 49 H 0.126 0.874 50 H 0.123 0.877 51 H 0.124 0.876 52 H 0.122 0.878 For the charges calculated by CM2: Dipole moment (debyes) X Y Z Total from point charges 4.019 4.362 -4.489 7.439 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 O -0.331 6.331 2 C -0.147 4.147 3 Si 0.581 3.419 4 C -0.259 4.259 5 C -0.084 4.084 6 N -0.139 5.139 7 C -0.010 4.010 8 C 0.282 3.718 9 O -0.374 6.374 10 N -0.380 5.380 11 C 0.043 3.957 12 H 0.107 0.893 13 C -0.168 4.168 14 C -0.146 4.146 15 C -0.108 4.108 16 H 0.097 0.903 17 C 0.001 3.999 18 O -0.379 6.379 19 C -0.148 4.148 20 C -0.162 4.162 21 N -0.067 5.067 22 C -0.034 4.034 23 C -0.009 4.009 24 C -0.111 4.111 25 C -0.138 4.138 26 C -0.132 4.132 27 C -0.139 4.139 28 C -0.112 4.112 29 H 0.218 0.782 30 H 0.133 0.867 31 H -0.184 1.184 32 H -0.179 1.179 33 H 0.195 0.805 34 H 0.134 0.866 35 H 0.146 0.854 36 H 0.234 0.766 37 H 0.088 0.912 38 H 0.083 0.917 39 H 0.093 0.907 40 H 0.081 0.919 41 H 0.064 0.936 42 H 0.065 0.935 43 H 0.228 0.772 44 H 0.090 0.910 45 H 0.087 0.913 46 H 0.089 0.911 47 H 0.091 0.909 48 H 0.143 0.857 49 H 0.144 0.856 50 H 0.141 0.859 51 H 0.142 0.858 52 H 0.140 0.860 Dipole moment (debyes) X Y Z Total from point charges 3.959 3.824 -4.269 6.966 hybrid contribution -0.757 0.570 -0.580 1.111 sum 3.203 4.394 -4.850 7.286 Atomic orbital electron populations 1.86632 1.17887 1.79787 1.48825 1.24348 0.83534 1.09704 0.97080 0.95273 0.91502 0.74545 0.80545 1.20389 1.13208 0.93233 0.99028 1.22276 1.00304 0.95414 0.90429 1.47840 1.47681 1.04372 1.14030 1.21520 0.94216 1.02297 0.82988 1.20494 0.85043 0.83035 0.83189 1.90839 1.69427 1.39977 1.37127 1.46001 1.26205 1.49121 1.16688 1.21038 0.84011 0.90296 1.00325 0.89279 1.21924 0.98318 0.94697 1.01850 1.21479 0.94822 0.99504 0.98760 1.21200 0.95497 0.93776 1.00325 0.90284 1.22226 0.96503 0.83816 0.97394 1.86474 1.38359 1.25422 1.87636 1.21460 0.96973 0.94098 1.02286 1.21773 0.95835 0.99426 0.99132 1.77816 1.17859 1.01596 1.09428 1.20595 0.97889 0.92356 0.92587 1.17701 0.98582 0.93230 0.91421 1.21118 0.99234 0.98888 0.91880 1.21077 0.98286 0.94666 0.99726 1.21252 0.98129 0.98124 0.95680 1.21089 1.00337 0.99008 0.93485 1.21161 0.96862 0.93316 0.99824 0.78151 0.86747 1.18440 1.17938 0.80538 0.86560 0.85432 0.76553 0.91228 0.91724 0.90743 0.91861 0.93618 0.93549 0.77241 0.91008 0.91269 0.91093 0.90923 0.85696 0.85606 0.85882 0.85775 0.85993 Number of geometries 1 Number of calculations of the screened coulomb radii 1 The total number of SCF iterations 13. 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 O -0.53 -1.28 12.88 -35.23 -0.45 -1.74 16 2 C -0.04 -0.04 2.50 36.01 0.09 0.05 16 3 Si 0.68 -0.58 30.80 -169.99 -5.24 -5.82 16 4 C -0.25 -0.54 5.35 -104.76 -0.56 -1.10 16 5 C 0.06 0.02 9.96 -17.26 -0.17 -0.15 16 6 N -0.36 -0.66 3.40 -62.58 -0.21 -0.87 16 7 C 0.12 -0.25 6.99 -5.19 -0.04 -0.28 16 8 C 0.50 0.33 7.82 -10.99 -0.09 0.25 16 9 O -0.50 -2.99 15.37 5.56 0.09 -2.90 16 10 N -0.73 1.56 5.08 -53.81 -0.27 1.29 16 11 C 0.15 -0.13 2.45 -67.93 -0.17 -0.29 16 12 H 0.09 0.13 7.58 -51.93 -0.39 -0.27 16 13 C -0.13 -0.05 5.49 -26.73 -0.15 -0.20 16 14 C -0.11 -0.05 5.10 -26.73 -0.14 -0.19 16 15 C -0.09 0.07 2.43 -90.62 -0.22 -0.15 16 16 H 0.08 -0.01 8.14 -51.93 -0.42 -0.43 16 17 C 0.08 -0.06 6.43 37.15 0.24 0.17 16 18 O -0.57 -0.37 13.20 -35.23 -0.46 -0.83 16 19 C -0.11 0.26 5.20 -26.73 -0.14 0.12 16 20 C -0.12 0.36 5.39 -26.73 -0.14 0.21 16 21 N -0.25 -1.34 11.67 33.43 0.39 -0.95 16 22 C 0.13 0.63 6.67 -82.78 -0.55 0.07 16 23 C -0.01 -0.03 5.74 -104.89 -0.60 -0.64 16 24 C -0.09 -0.31 8.89 -39.25 -0.35 -0.66 16 25 C -0.12 -0.29 10.03 -39.63 -0.40 -0.69 16 26 C -0.11 -0.27 10.04 -39.52 -0.40 -0.66 16 27 C -0.12 -0.35 10.03 -39.63 -0.40 -0.75 16 28 C -0.09 -0.39 9.78 -39.26 -0.38 -0.78 16 29 H 0.37 -0.20 8.70 45.56 0.40 0.20 16 30 H 0.11 0.10 5.57 -51.93 -0.29 -0.19 16 31 H -0.26 0.00 7.00 56.52 0.40 0.40 16 32 H -0.25 -0.63 7.11 56.52 0.40 -0.22 16 33 H 0.18 -0.55 6.33 -52.49 -0.33 -0.89 16 34 H 0.12 -0.46 8.14 -51.93 -0.42 -0.89 16 35 H 0.13 -0.81 8.07 -51.92 -0.42 -1.23 16 36 H 0.40 -2.18 8.60 -40.82 -0.35 -2.53 16 37 H 0.07 0.20 8.14 -51.93 -0.42 -0.22 16 38 H 0.06 -0.01 8.14 -51.93 -0.42 -0.43 16 39 H 0.07 0.14 7.82 -51.93 -0.41 -0.26 16 40 H 0.06 0.06 8.14 -51.93 -0.42 -0.36 16 41 H 0.05 -0.08 8.14 -51.93 -0.42 -0.51 16 42 H 0.05 -0.06 8.14 -51.93 -0.42 -0.48 16 43 H 0.38 -1.29 9.12 45.56 0.42 -0.87 16 44 H 0.07 -0.22 8.14 -51.92 -0.42 -0.64 16 45 H 0.07 -0.14 8.14 -51.93 -0.42 -0.57 16 46 H 0.07 -0.25 8.14 -51.93 -0.42 -0.67 16 47 H 0.07 -0.25 8.14 -51.93 -0.42 -0.67 16 48 H 0.13 0.22 6.21 -52.49 -0.33 -0.10 16 49 H 0.13 0.03 8.06 -52.49 -0.42 -0.39 16 50 H 0.12 0.03 8.06 -52.48 -0.42 -0.40 16 51 H 0.12 0.11 8.06 -52.49 -0.42 -0.31 16 52 H 0.12 0.47 8.06 -52.48 -0.42 0.05 16 LS Contribution 422.59 15.07 6.37 6.37 Total: 0.00 -12.38 422.59 -11.61 -23.98 By element: Atomic # 1 Polarization: -5.62 SS G_CDS: -7.25 Total: -12.87 kcal Atomic # 6 Polarization: -1.10 SS G_CDS: -4.56 Total: -5.66 kcal Atomic # 7 Polarization: -0.44 SS G_CDS: -0.10 Total: -0.53 kcal Atomic # 8 Polarization: -4.64 SS G_CDS: -0.83 Total: -5.47 kcal Atomic # 14 Polarization: -0.58 SS G_CDS: -5.24 Total: -5.82 kcal Total LS contribution 6.37 Total: 6.37 kcal Total: -12.38 -11.61 -23.98 kcal The number of atoms in the molecule is 52 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_22_7207542_12622748.mol2 52 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 9.262 kcal (2) G-P(sol) polarization free energy of solvation -12.379 kcal (3) G-ENP(sol) elect.-nuc.-pol. free energy of system -3.117 kcal (4) G-CDS(sol) cavity-dispersion-solvent structure free energy -11.606 kcal (5) G-P-CDS(sol) = G-P(sol) + G-CDS(sol) = (2) + (4) -23.985 kcal (6) G-S(sol) free energy of system = (1) + (5) -14.723 kcal 1SCF WAS SPECIFIED, SO GEOMETRY OPTIMIZATION WAS NOT USED Total computer time = 1.07 seconds