Wall clock time and date at job start Sun Mar 7 2021 10:42:46 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 O 1.42903 * 1 3 3 C 1.35992 * 116.99565 * 2 1 4 4 C 1.38670 * 119.99783 * 0.02562 * 3 2 1 5 5 C 1.38154 * 119.99915 * 179.97438 * 4 3 2 6 6 C 1.38681 * 119.99691 * 0.02562 * 5 4 3 7 7 O 1.35986 * 119.99886 * 180.02562 * 6 5 4 8 8 C 1.42905 * 117.00016 * 0.02562 * 7 6 5 9 9 C 1.50700 * 109.47021 * 179.97438 * 8 7 6 10 10 O 1.21306 * 119.99782 * 0.02562 * 9 8 7 11 11 N 1.34778 * 119.99976 * 179.97438 * 9 8 7 12 12 C 1.40762 * 119.99649 * 180.02562 * 11 9 8 13 13 C 1.35758 * 126.07598 * 179.97438 * 12 11 9 14 14 N 1.35805 * 107.56427 * 179.97438 * 13 12 11 15 15 C 1.46504 * 126.05914 * 179.97438 * 14 13 12 16 16 Si 1.86296 * 109.46982 * 89.99826 * 15 14 13 17 17 C 1.86305 * 109.46706 * 60.00643 * 16 15 14 18 18 C 1.86299 * 109.47208 * 180.02562 * 16 15 14 19 19 C 1.86299 * 109.47658 * 300.00585 * 16 15 14 20 20 N 1.40078 * 107.89133 * 0.02562 * 14 13 12 21 21 C 1.30826 * 108.36573 * 359.71099 * 20 14 13 22 22 C 1.38712 * 120.00083 * 0.02562 * 6 5 4 23 23 C 1.50701 * 119.99957 * 179.97438 * 22 6 5 24 24 O 1.42895 * 109.47079 * 75.00898 * 23 22 6 25 25 Si 1.86301 * 109.47149 * 195.00393 * 23 22 6 26 26 H 1.48506 * 109.47302 * 59.99966 * 25 23 22 27 27 H 1.48503 * 109.47334 * 300.00099 * 25 23 22 28 28 H 1.48500 * 109.47376 * 180.02562 * 25 23 22 29 29 C 1.38110 * 119.99936 * 0.24058 * 22 6 5 30 30 H 1.08997 * 109.46779 * 180.02562 * 1 2 3 31 31 H 1.08998 * 109.46818 * 299.99390 * 1 2 3 32 32 H 1.08993 * 109.47435 * 59.99660 * 1 2 3 33 33 H 1.07996 * 120.00591 * 359.97438 * 4 3 2 34 34 H 1.08000 * 120.00545 * 179.97438 * 5 4 3 35 35 H 1.08999 * 109.47003 * 300.00171 * 8 7 6 36 36 H 1.09004 * 109.46705 * 60.00512 * 8 7 6 37 37 H 0.96998 * 119.99988 * 0.02562 * 11 9 8 38 38 H 1.07998 * 126.21820 * 359.97438 * 13 12 11 39 39 H 1.08992 * 109.47502 * 210.00599 * 15 14 13 40 40 H 1.09001 * 109.46709 * 330.00401 * 15 14 13 41 41 H 1.08999 * 109.47199 * 60.00334 * 17 16 15 42 42 H 1.09006 * 109.46870 * 179.97438 * 17 16 15 43 43 H 1.08997 * 109.47155 * 299.99840 * 17 16 15 44 44 H 1.09001 * 109.46724 * 60.00428 * 18 16 15 45 45 H 1.08997 * 109.47571 * 180.02562 * 18 16 15 46 46 H 1.09000 * 109.47560 * 300.00819 * 18 16 15 47 47 H 1.09002 * 109.46503 * 59.99879 * 19 16 15 48 48 H 1.08994 * 109.47040 * 179.97438 * 19 16 15 49 49 H 1.09001 * 109.47052 * 300.00049 * 19 16 15 50 50 H 1.08000 * 125.82726 * 180.26196 * 21 20 14 51 51 H 1.09001 * 109.47380 * 315.00295 * 23 22 6 52 52 H 0.96708 * 113.99763 * 299.99687 * 24 23 22 53 53 H 1.08008 * 119.99731 * 179.72155 * 29 22 6 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 8 1.4290 0.0000 0.0000 3 6 2.0463 1.2117 0.0000 4 6 1.2909 2.3746 0.0005 5 6 1.9180 3.6057 0.0011 6 6 3.3029 3.6783 0.0005 7 8 3.9202 4.8900 0.0016 8 6 3.0802 6.0461 0.0038 9 6 3.9343 7.2876 0.0053 10 8 5.1436 7.1918 0.0041 11 7 3.3546 8.5044 0.0070 12 6 4.1525 9.6640 0.0079 13 6 3.7017 10.9446 0.0090 14 7 4.7869 11.7611 0.0102 15 6 4.7639 13.2259 0.0111 16 14 4.7531 13.8456 1.7679 17 6 6.2835 13.2481 2.6465 18 6 4.7231 15.7084 1.7692 19 6 3.2416 13.2002 2.6451 20 7 5.9322 10.9546 0.0091 21 6 5.5544 9.7021 0.0136 22 6 4.0585 2.5150 -0.0011 23 6 5.5634 2.5940 -0.0022 24 8 6.0219 2.9658 1.2991 25 14 6.2721 0.9327 -0.4590 26 1 5.7957 0.5460 -1.8113 27 1 5.8305 -0.0829 0.5304 28 1 7.7551 1.0103 -0.4595 29 6 3.4315 1.2844 -0.0072 30 1 -0.3633 -1.0277 0.0005 31 1 -0.3633 0.5137 0.8900 32 1 -0.3634 0.5138 -0.8899 33 1 0.2125 2.3181 0.0005 34 1 1.3298 4.5114 0.0011 35 1 2.4502 6.0391 -0.8856 36 1 2.4517 6.0368 0.8944 37 1 2.3877 8.5809 0.0076 38 1 2.6679 11.2571 0.0091 39 1 5.6482 13.6036 -0.5020 40 1 3.8685 13.5755 -0.5029 41 1 7.1679 13.6258 2.1334 42 1 6.2774 13.6111 3.6743 43 1 6.3007 12.1583 2.6457 44 1 3.8275 16.0576 1.2553 45 1 4.7173 16.0710 2.7971 46 1 5.6072 16.0865 1.2559 47 1 3.2588 12.1103 2.6443 48 1 3.2356 13.5624 3.6730 49 1 2.3462 13.5498 2.1310 50 1 6.2102 8.8440 0.0181 51 1 5.8876 3.3394 -0.7285 52 1 5.6914 3.8225 1.6026 53 1 4.0199 0.3787 -0.0130 RHF calculation, no. of doubly occupied orbitals= 71 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_16608520_4974302.mol2 53 1SCF WAS SPECIFIED, SO GEOMETRY OPTIMIZATION WAS NOT USED SCF FIELD WAS ACHIEVED Wall clock time and date at job start Sun Mar 7 2021 10:42:46 Heat of formation + Delta-G solvation = -64.159965 kcal Electronic energy + Delta-G solvation = -33882.880099 eV Core-core repulsion = 29267.492656 eV Total energy + Delta-G solvation = -4615.387443 eV No. of doubly occupied orbitals = 71 Molecular weight (most abundant/longest-lived isotopes) = 393.172 amu Computer time = 1.43 seconds Orbital eigenvalues (eV) -42.33754 -41.75910 -40.66980 -39.76147 -38.01554 -36.82826 -35.90193 -33.72729 -32.78670 -31.40013 -30.56192 -29.75636 -29.24056 -28.59326 -28.45959 -28.08977 -26.82568 -26.48077 -24.00629 -22.82945 -22.38123 -22.09807 -21.03137 -19.53889 -18.85368 -18.64808 -18.09147 -17.84711 -17.58628 -17.36617 -16.74584 -16.69222 -16.44018 -16.17477 -15.72735 -15.57068 -15.48526 -15.23863 -15.07216 -14.73797 -14.39243 -14.33841 -14.22329 -14.16347 -14.04034 -13.95828 -13.82429 -13.65946 -13.57443 -13.48435 -13.43788 -13.27531 -13.18957 -13.06671 -12.70046 -12.29484 -12.19837 -12.06839 -11.68141 -11.61189 -11.51417 -11.39517 -11.09271 -11.06561 -10.95762 -10.89919 -10.83925 -10.05852 -10.05040 -8.81059 -8.64031 0.08665 0.36268 0.69667 1.12293 1.27871 1.36214 1.49752 1.51961 1.65936 1.66745 1.71913 1.92061 2.11746 2.13628 2.27053 2.38868 2.51671 2.57762 3.35149 3.55367 3.63264 3.67134 3.73485 3.81121 3.83786 3.97486 4.08583 4.12993 4.13765 4.17621 4.23629 4.34605 4.45678 4.48583 4.53722 4.54330 4.56029 4.58475 4.59612 4.61492 4.66838 4.71311 4.75405 4.78848 5.05023 5.08254 5.12865 5.19988 5.30316 5.37257 5.41487 5.56074 5.87627 5.98660 6.07950 6.27272 6.39430 6.77499 6.78250 7.03658 Net atomic charges, atomic populations, and dipole contributions calculated with CM2 Atom NO. Type Charge No. of electrons 1 C 0.026 3.974 2 O -0.315 6.315 3 C 0.104 3.896 4 C -0.165 4.165 5 C -0.165 4.165 6 C 0.103 3.897 7 O -0.288 6.288 8 C 0.046 3.954 9 C 0.516 3.484 10 O -0.500 6.500 11 N -0.639 5.639 12 C -0.006 4.006 13 C 0.017 3.983 14 N -0.343 5.343 15 C -0.169 4.169 16 Si 1.020 2.980 17 C -0.447 4.447 18 C -0.450 4.450 19 C -0.450 4.450 20 N -0.280 5.280 21 C 0.048 3.952 22 C -0.097 4.097 23 C -0.084 4.084 24 O -0.535 6.535 25 Si 1.013 2.987 26 H -0.248 1.248 27 H -0.251 1.251 28 H -0.255 1.255 29 C -0.116 4.116 30 H 0.102 0.898 31 H 0.056 0.944 32 H 0.056 0.944 33 H 0.135 0.865 34 H 0.134 0.866 35 H 0.096 0.904 36 H 0.097 0.903 37 H 0.411 0.589 38 H 0.177 0.823 39 H 0.105 0.895 40 H 0.099 0.901 41 H 0.071 0.929 42 H 0.069 0.931 43 H 0.071 0.929 44 H 0.071 0.929 45 H 0.075 0.925 46 H 0.072 0.928 47 H 0.069 0.931 48 H 0.074 0.926 49 H 0.070 0.930 50 H 0.205 0.795 51 H 0.082 0.918 52 H 0.371 0.629 53 H 0.142 0.858 For the charges calculated by CM2: Dipole moment (debyes) X Y Z Total from point charges -9.287 4.258 -0.369 10.223 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.069 4.069 2 O -0.229 6.229 3 C 0.059 3.941 4 C -0.184 4.184 5 C -0.184 4.184 6 C 0.058 3.942 7 O -0.201 6.201 8 C -0.033 4.033 9 C 0.302 3.698 10 O -0.373 6.373 11 N -0.286 5.286 12 C -0.110 4.110 13 C -0.131 4.131 14 N -0.126 5.126 15 C -0.335 4.335 16 Si 1.214 2.786 17 C -0.555 4.555 18 C -0.558 4.558 19 C -0.558 4.558 20 N -0.103 5.103 21 C -0.138 4.138 22 C -0.098 4.098 23 C -0.192 4.192 24 O -0.341 6.341 25 Si 0.833 3.167 26 H -0.171 1.171 27 H -0.175 1.175 28 H -0.179 1.179 29 C -0.135 4.135 30 H 0.121 0.879 31 H 0.075 0.925 32 H 0.075 0.925 33 H 0.153 0.847 34 H 0.152 0.848 35 H 0.114 0.886 36 H 0.114 0.886 37 H 0.246 0.754 38 H 0.194 0.806 39 H 0.124 0.876 40 H 0.117 0.883 41 H 0.090 0.910 42 H 0.088 0.912 43 H 0.090 0.910 44 H 0.090 0.910 45 H 0.094 0.906 46 H 0.092 0.908 47 H 0.089 0.911 48 H 0.093 0.907 49 H 0.089 0.911 50 H 0.221 0.779 51 H 0.100 0.900 52 H 0.219 0.781 53 H 0.159 0.841 Dipole moment (debyes) X Y Z Total from point charges -6.491 3.554 -0.701 7.433 hybrid contribution -2.644 0.872 0.015 2.784 sum -9.135 4.426 -0.686 10.174 Atomic orbital electron populations 1.23099 0.78056 1.04373 1.01343 1.86117 1.21570 1.26314 1.88902 1.18802 0.91897 0.84562 0.98845 1.21073 1.00305 0.91510 1.05505 1.21020 0.92845 0.98651 1.05842 1.18600 0.90744 0.84468 1.00420 1.86107 1.33712 1.10961 1.89284 1.22350 0.92912 0.83844 1.04237 1.19995 0.88866 0.83656 0.77315 1.90678 1.14066 1.85945 1.46567 1.42481 1.09562 1.03037 1.73512 1.17584 0.91342 0.85473 1.16605 1.21861 0.94076 0.89153 1.08040 1.48002 1.03252 1.05530 1.55793 1.26971 1.07135 0.77341 1.22057 0.79374 0.68274 0.66872 0.64105 1.26074 1.15053 1.06637 1.07724 1.26098 1.04805 1.19800 1.05121 1.26126 1.14618 1.06604 1.08417 1.77489 1.07514 1.01609 1.23698 1.24486 0.93519 0.95444 1.00380 1.19430 0.96725 0.89902 1.03740 1.25229 0.91315 1.13420 0.89274 1.86673 1.73715 1.42177 1.31522 0.84495 0.79744 0.72881 0.79616 1.17141 1.17490 1.17858 1.20231 0.92025 0.97404 1.03792 0.87921 0.92549 0.92543 0.84716 0.84826 0.88610 0.88560 0.75408 0.80638 0.87631 0.88280 0.91006 0.91205 0.90980 0.90976 0.90553 0.90841 0.91132 0.90714 0.91127 0.77859 0.89961 0.78059 0.84097 Number of geometries 1 Number of calculations of the screened coulomb radii 1 The total number of SCF iterations 14. 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.03 -0.05 9.81 101.05 0.99 0.94 16 2 O -0.32 -0.62 10.50 -19.71 -0.21 -0.82 16 3 C 0.10 0.21 6.69 -39.25 -0.26 -0.05 16 4 C -0.17 0.00 9.03 -39.45 -0.36 -0.35 16 5 C -0.16 0.09 9.03 -39.45 -0.36 -0.26 16 6 C 0.10 0.13 6.62 -39.24 -0.26 -0.13 16 7 O -0.29 -0.34 8.56 -39.61 -0.34 -0.68 16 8 C 0.05 -0.11 5.22 36.01 0.19 0.08 16 9 C 0.52 0.17 7.88 -10.98 -0.09 0.08 16 10 O -0.50 -2.28 13.66 -14.37 -0.20 -2.47 16 11 N -0.64 1.09 5.47 -10.04 -0.05 1.03 16 12 C -0.01 -0.01 6.90 -83.89 -0.58 -0.58 16 13 C 0.02 -0.01 11.44 -16.68 -0.19 -0.20 16 14 N -0.34 -0.23 3.01 -9.17 -0.03 -0.25 16 15 C -0.17 0.39 5.59 59.85 0.33 0.72 16 16 Si 1.02 -3.45 12.79 -169.99 -2.17 -5.62 16 17 C -0.45 0.18 8.65 101.05 0.87 1.05 16 18 C -0.45 1.42 8.65 101.05 0.87 2.29 16 19 C -0.45 1.24 8.65 101.05 0.87 2.12 16 20 N -0.28 -1.08 12.56 30.70 0.39 -0.70 16 21 C 0.05 0.18 10.84 -17.23 -0.19 -0.01 16 22 C -0.10 -0.24 5.28 -104.49 -0.55 -0.79 16 23 C -0.08 -0.22 2.33 36.01 0.08 -0.14 16 24 O -0.54 -1.96 12.65 -35.23 -0.45 -2.40 16 25 Si 1.01 2.16 25.73 -169.99 -4.37 -2.21 16 26 H -0.25 -0.40 7.11 56.52 0.40 0.01 16 27 H -0.25 -0.66 6.74 56.52 0.38 -0.28 16 28 H -0.25 -0.79 7.11 56.52 0.40 -0.39 16 29 C -0.12 -0.31 8.59 -39.45 -0.34 -0.65 16 30 H 0.10 -0.33 8.14 -51.93 -0.42 -0.76 16 31 H 0.06 -0.16 7.66 -51.93 -0.40 -0.56 16 32 H 0.06 -0.17 7.66 -51.93 -0.40 -0.56 16 33 H 0.14 -0.32 6.29 -52.49 -0.33 -0.65 16 34 H 0.13 -0.50 6.29 -52.49 -0.33 -0.83 16 35 H 0.10 -0.45 7.66 -51.93 -0.40 -0.85 16 36 H 0.10 -0.47 7.66 -51.93 -0.40 -0.87 16 37 H 0.41 -2.35 8.83 -40.82 -0.36 -2.71 16 38 H 0.18 -0.81 8.06 -52.49 -0.42 -1.23 16 39 H 0.11 -0.22 7.96 -51.93 -0.41 -0.63 16 40 H 0.10 -0.45 7.96 -51.93 -0.41 -0.86 16 41 H 0.07 -0.01 7.96 -51.93 -0.41 -0.42 16 42 H 0.07 -0.05 7.96 -51.93 -0.41 -0.47 16 43 H 0.07 0.11 7.96 -51.93 -0.41 -0.30 16 44 H 0.07 -0.26 7.96 -51.93 -0.41 -0.68 16 45 H 0.08 -0.21 7.96 -51.93 -0.41 -0.63 16 46 H 0.07 -0.22 7.96 -51.93 -0.41 -0.63 16 47 H 0.07 -0.13 7.96 -51.93 -0.41 -0.54 16 48 H 0.07 -0.18 7.96 -51.93 -0.41 -0.59 16 49 H 0.07 -0.26 7.96 -51.93 -0.41 -0.67 16 50 H 0.21 0.96 5.92 -52.49 -0.31 0.65 16 51 H 0.08 0.18 7.82 -51.93 -0.41 -0.23 16 52 H 0.37 0.72 8.60 45.56 0.39 1.12 16 53 H 0.14 0.30 4.88 -52.48 -0.26 0.05 16 LS Contribution 440.17 15.07 6.63 6.63 Total: 0.00 -10.75 440.17 -7.15 -17.90 By element: Atomic # 1 Polarization: -7.11 SS G_CDS: -7.40 Total: -14.52 kcal Atomic # 6 Polarization: 3.07 SS G_CDS: 1.05 Total: 4.12 kcal Atomic # 7 Polarization: -0.22 SS G_CDS: 0.30 Total: 0.08 kcal Atomic # 8 Polarization: -5.20 SS G_CDS: -1.19 Total: -6.39 kcal Atomic # 14 Polarization: -1.28 SS G_CDS: -6.55 Total: -7.83 kcal Total LS contribution 6.63 Total: 6.63 kcal Total: -10.75 -7.15 -17.90 kcal The number of atoms in the molecule is 53 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_16608520_4974302.mol2 53 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 -46.261 kcal (2) G-P(sol) polarization free energy of solvation -10.748 kcal (3) G-ENP(sol) elect.-nuc.-pol. free energy of system -57.009 kcal (4) G-CDS(sol) cavity-dispersion-solvent structure free energy -7.151 kcal (5) G-P-CDS(sol) = G-P(sol) + G-CDS(sol) = (2) + (4) -17.899 kcal (6) G-S(sol) free energy of system = (1) + (5) -64.160 kcal 1SCF WAS SPECIFIED, SO GEOMETRY OPTIMIZATION WAS NOT USED Total computer time = 1.43 seconds