Wall clock time and date at job start Sat Apr 11 2020 06:10:10 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= -1 * 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 H 1.09005 * 109.47035 * 2 1 4 4 C 1.52999 * 109.47195 * 119.99294 * 2 1 3 5 5 C 1.50699 * 109.47378 * 295.00711 * 4 2 1 6 6 O 1.21277 * 119.99869 * 359.97438 * 5 4 2 7 7 N 1.34772 * 119.99926 * 179.72429 * 5 4 2 8 8 C 1.46503 * 120.00014 * 180.27532 * 7 5 4 9 9 C 1.53006 * 109.47157 * 176.60858 * 8 7 5 10 10 N 1.46500 * 109.54324 * 54.52714 * 9 8 7 11 11 C 1.36939 * 119.99642 * 62.05117 * 10 9 8 12 12 H 1.07999 * 120.00426 * 354.01982 * 11 10 9 13 13 C 1.38817 * 119.99943 * 174.01613 * 11 10 9 14 14 N 1.31615 * 119.99753 * 0.02562 * 13 11 10 15 15 Si 1.86791 * 120.00080 * 179.97438 * 13 11 10 16 16 H 1.48502 * 109.47536 * 0.02562 * 15 13 11 17 17 H 1.48501 * 109.47102 * 120.00179 * 15 13 11 18 18 H 1.48499 * 109.47292 * 240.00198 * 15 13 11 19 19 C 1.53174 * 109.53979 * 294.26256 * 9 8 7 20 20 C 1.53098 * 109.16349 * 183.15471 * 19 9 8 21 21 O 1.42900 * 109.41241 * 57.60707 * 20 19 9 22 22 C 1.42904 * 114.09898 * 298.84269 * 21 20 19 23 23 C 1.53093 * 109.41259 * 61.16112 * 22 21 20 24 24 C 1.52997 * 109.47301 * 239.99776 * 2 1 3 25 25 H 1.09003 * 109.82021 * 185.56649 * 24 2 1 26 26 C 1.54385 * 109.83196 * 64.58597 * 24 2 1 27 27 C 1.55299 * 102.25402 * 219.12635 * 26 24 2 28 28 C 1.54954 * 100.93990 * 34.13672 * 27 26 24 29 29 O 1.43240 * 109.83082 * 306.45028 * 24 2 1 30 30 H 1.08999 * 109.46656 * 179.97438 * 1 2 3 31 31 H 1.08991 * 109.47155 * 300.00327 * 1 2 3 32 32 H 1.09003 * 109.46858 * 60.00780 * 1 2 3 33 33 H 1.08997 * 109.47183 * 174.99823 * 4 2 1 34 34 H 1.09007 * 109.46762 * 55.00664 * 4 2 1 35 35 H 0.97001 * 120.00324 * 0.26884 * 7 5 4 36 36 H 1.09007 * 109.46825 * 296.60407 * 8 7 5 37 37 H 1.08992 * 109.47215 * 56.60436 * 8 7 5 38 38 H 0.96994 * 119.99945 * 242.04518 * 10 9 8 39 39 H 0.97001 * 120.00495 * 179.97438 * 14 13 11 40 40 H 1.08997 * 109.52816 * 303.05984 * 19 9 8 41 41 H 1.09001 * 109.52495 * 63.24916 * 19 9 8 42 42 H 1.08998 * 109.47648 * 297.62774 * 20 19 9 43 43 H 1.08998 * 109.49265 * 177.60033 * 20 19 9 44 44 H 1.09002 * 109.48174 * 301.17892 * 22 21 20 45 45 H 1.08994 * 109.48861 * 181.14828 * 22 21 20 46 46 H 1.08999 * 109.51973 * 182.48120 * 23 22 21 47 47 H 1.09005 * 109.57413 * 62.32892 * 23 22 21 48 48 H 1.08995 * 110.84742 * 337.32419 * 26 24 2 49 49 H 1.09008 * 110.84329 * 100.93758 * 26 24 2 50 50 H 1.08998 * 111.12803 * 276.21485 * 27 26 24 51 51 H 1.08994 * 111.12539 * 152.05054 * 27 26 24 52 52 H 1.09007 * 110.30717 * 204.63986 * 28 27 26 53 53 H 1.08993 * 110.39704 * 82.51016 * 28 27 26 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 1 1.8934 1.0277 0.0000 4 6 2.0400 -0.7211 1.2493 5 6 1.6413 0.0572 2.4766 6 8 1.0116 1.0875 2.3632 7 7 1.9800 -0.3943 3.7004 8 6 1.5871 0.3592 4.8938 9 6 2.0184 -0.4073 6.1459 10 7 3.4528 -0.6963 6.0738 11 6 4.3569 0.3313 6.0312 12 1 4.0115 1.3516 5.9533 13 6 5.7174 0.0616 6.0884 14 7 6.1383 -1.1818 6.1829 15 14 6.9507 1.4633 6.0296 16 1 6.2268 2.7555 5.9227 17 1 7.7682 1.4556 7.2693 18 1 7.8376 1.2933 4.8507 19 6 1.2335 -1.7194 6.2390 20 6 1.6138 -2.4417 7.5342 21 8 1.3472 -1.5889 8.6494 22 6 2.0890 -0.3674 8.6404 23 6 1.7277 0.4381 7.3897 24 6 2.0401 -0.7213 -1.2492 25 1 3.1255 -0.8115 -1.2066 26 6 1.6271 0.0639 -2.5127 27 6 1.2542 -1.0745 -3.5010 28 6 0.6557 -2.1230 -2.5296 29 8 1.4473 -2.0224 -1.3351 30 1 -0.3632 -1.0277 -0.0005 31 1 -0.3633 0.5138 0.8899 32 1 -0.3633 0.5137 -0.8901 33 1 3.1263 -0.7993 1.2057 34 1 1.6055 -1.7197 1.2952 35 1 2.4839 -1.2182 3.7911 36 1 2.0700 1.3364 4.8818 37 1 0.5050 0.4887 4.9018 38 1 3.7621 -1.6154 6.0563 39 1 7.0890 -1.3704 6.2225 40 1 0.1648 -1.5053 6.2417 41 1 1.4772 -2.3511 5.3847 42 1 2.6746 -2.6915 7.5130 43 1 1.0270 -3.3555 7.6269 44 1 3.1560 -0.5902 8.6337 45 1 1.8444 0.2131 9.5299 46 1 2.3248 1.3495 7.3584 47 1 0.6691 0.6968 7.4151 48 1 0.7652 0.6995 -2.3099 49 1 2.4622 0.6512 -2.8946 50 1 2.1395 -1.4645 -4.0032 51 1 0.5105 -0.7393 -4.2238 52 1 0.7381 -3.1232 -2.9549 53 1 -0.3867 -1.8892 -2.3134 RHF calculation, no. of doubly occupied orbitals= 66 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=-1 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 ZINC001484190882.mol2 53 1SCF WAS SPECIFIED, SO GEOMETRY OPTIMIZATION WAS NOT USED SCF FIELD WAS ACHIEVED Wall clock time and date at job start Sat Apr 11 2020 06:10:10 Heat of formation + Delta-G solvation = -132.969773 kcal Electronic energy + Delta-G solvation = -31582.134038 eV Core-core repulsion = 27453.405778 eV Total energy + Delta-G solvation = -4128.728261 eV No. of doubly occupied orbitals = 66 Molecular weight (most abundant/longest-lived isotopes) = 340.220 amu Computer time = 1.09 seconds Orbital eigenvalues (eV) -40.15941 -39.51811 -38.67599 -35.87070 -34.39404 -33.38612 -32.29054 -32.15138 -30.62074 -29.18839 -28.95038 -27.95351 -26.20752 -24.09684 -23.54081 -22.97138 -22.13024 -21.70245 -20.12603 -19.48303 -18.57594 -17.28387 -17.12773 -16.70343 -15.96612 -15.66506 -15.58089 -15.34943 -15.26200 -14.62507 -14.46218 -14.27970 -14.11654 -14.00905 -13.79978 -13.58145 -13.13292 -12.75711 -12.64594 -12.39320 -12.13239 -12.04315 -11.95948 -11.78124 -11.66907 -11.54676 -11.46428 -11.21451 -10.97715 -10.94360 -10.89617 -10.77097 -10.70270 -10.61769 -10.33137 -10.01594 -9.75502 -9.69181 -9.59238 -9.24800 -9.19914 -8.77565 -8.38688 -7.10607 -5.88075 -3.18449 3.14361 3.36454 3.44707 3.44997 3.65102 3.78307 4.28775 4.40193 4.50159 4.58872 4.67247 4.92477 4.96611 5.00904 5.13434 5.18658 5.27310 5.34789 5.39676 5.45587 5.49267 5.52743 5.54898 5.56997 5.58331 5.66770 5.69195 5.73085 5.80263 5.83802 5.89356 6.00616 6.03310 6.11020 6.12396 6.14930 6.25810 6.33562 6.34556 6.43435 6.52365 6.56683 6.60618 6.74925 6.81201 6.93444 7.09319 7.64591 7.73994 8.11235 8.29228 8.76592 9.41930 9.92658 10.29499 11.33669 Molecular weight = 340.22amu Principal moments of inertia in cm(-1) A = 0.011474 B = 0.003053 C = 0.002655 Principal moments of inertia in units of 10**(-40)*gram-cm**2 A = 2439.703533 B = 9169.369538 C =10544.855465 Net atomic charges, atomic populations, and dipole contributions calculated with CM2 Atom NO. Type Charge No. of electrons 1 C -0.149 4.149 2 C -0.081 4.081 3 H 0.085 0.915 4 C -0.133 4.133 5 C 0.510 3.490 6 O -0.558 6.558 7 N -0.702 5.702 8 C 0.119 3.881 9 C 0.189 3.811 10 N -0.713 5.713 11 C -0.235 4.235 12 H 0.070 0.930 13 C -0.006 4.006 14 N -0.929 5.929 15 Si 0.901 3.099 16 H -0.281 1.281 17 H -0.292 1.292 18 H -0.292 1.292 19 C -0.150 4.150 20 C 0.066 3.934 21 O -0.394 6.394 22 C 0.070 3.930 23 C -0.152 4.152 24 C 0.083 3.917 25 H 0.081 0.919 26 C -0.154 4.154 27 C -0.157 4.157 28 C 0.050 3.950 29 O -0.373 6.373 30 H 0.052 0.948 31 H 0.076 0.924 32 H 0.046 0.954 33 H 0.091 0.909 34 H 0.094 0.906 35 H 0.403 0.597 36 H 0.063 0.937 37 H 0.052 0.948 38 H 0.384 0.616 39 H 0.255 0.745 40 H 0.068 0.932 41 H 0.069 0.931 42 H 0.060 0.940 43 H 0.081 0.919 44 H 0.069 0.931 45 H 0.077 0.923 46 H 0.076 0.924 47 H 0.054 0.946 48 H 0.083 0.917 49 H 0.080 0.920 50 H 0.076 0.924 51 H 0.079 0.921 52 H 0.081 0.919 53 H 0.058 0.942 For the charges calculated by CM2: Dipole moment (debyes) X Y Z Total from point charges -10.969 -2.329 -12.962 17.139 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.206 4.206 2 C -0.100 4.100 3 H 0.103 0.897 4 C -0.173 4.173 5 C 0.296 3.704 6 O -0.438 6.438 7 N -0.352 5.352 8 C -0.008 4.008 9 C 0.097 3.903 10 N -0.355 5.355 11 C -0.366 4.366 12 H 0.088 0.912 13 C -0.203 4.203 14 N -0.577 5.577 15 Si 0.732 3.268 16 H -0.206 1.206 17 H -0.219 1.219 18 H -0.220 1.220 19 C -0.189 4.189 20 C -0.011 4.011 21 O -0.313 6.313 22 C -0.006 4.006 23 C -0.193 4.193 24 C 0.025 3.975 25 H 0.098 0.902 26 C -0.191 4.191 27 C -0.194 4.194 28 C -0.026 4.026 29 O -0.292 6.292 30 H 0.071 0.929 31 H 0.095 0.905 32 H 0.065 0.935 33 H 0.109 0.891 34 H 0.112 0.888 35 H 0.238 0.762 36 H 0.081 0.919 37 H 0.070 0.930 38 H 0.219 0.781 39 H 0.065 0.935 40 H 0.087 0.913 41 H 0.087 0.913 42 H 0.078 0.922 43 H 0.099 0.901 44 H 0.088 0.912 45 H 0.096 0.904 46 H 0.095 0.905 47 H 0.073 0.927 48 H 0.101 0.899 49 H 0.099 0.901 50 H 0.095 0.905 51 H 0.098 0.902 52 H 0.099 0.901 53 H 0.076 0.924 Dipole moment (debyes) X Y Z Total from point charges -11.402 -1.658 -12.707 17.153 hybrid contribution 1.450 0.411 -0.429 1.567 sum -9.952 -1.247 -13.136 16.527 Atomic orbital electron populations 1.21910 0.94533 1.01657 1.02483 1.21212 0.97237 0.99245 0.92335 0.89704 1.21644 1.02816 1.00393 0.92471 1.20578 0.79328 0.85648 0.84848 1.90602 1.41268 1.25014 1.86868 1.45801 1.56297 1.27816 1.05321 1.21767 0.97774 0.96026 0.85205 1.19828 0.86375 0.92363 0.91715 1.42485 0.98240 1.06712 1.88050 1.19106 0.86669 0.88363 1.42419 0.91183 1.24692 0.96923 0.96072 1.02564 1.69650 1.14045 1.16885 1.57114 0.86214 0.81741 0.81659 0.77158 1.20583 1.21897 1.21959 1.22342 1.00200 0.97696 0.98613 1.22463 1.00883 0.92296 0.85470 1.87792 1.67092 1.31280 1.45104 1.22448 0.96905 0.84500 0.96794 1.22454 1.00187 1.00255 0.96378 1.22855 0.97403 0.82796 0.94406 0.90157 1.23237 1.02745 0.98398 0.94738 1.23250 1.01628 0.96215 0.98335 1.23346 0.95051 0.99331 0.84902 1.88432 1.66106 1.32119 1.42573 0.92874 0.90520 0.93513 0.89109 0.88758 0.76198 0.91906 0.92958 0.78088 0.93527 0.91302 0.91257 0.92201 0.90103 0.91243 0.90434 0.90503 0.92702 0.89858 0.90145 0.90516 0.90189 0.90091 0.92363 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.15 -1.85 8.47 37.16 0.31 -1.54 16 2 C -0.08 -1.01 1.86 -90.62 -0.17 -1.18 16 3 H 0.08 1.14 7.95 -51.93 -0.41 0.73 16 4 C -0.13 -1.81 5.04 -27.88 -0.14 -1.95 16 5 C 0.51 8.59 7.55 -10.99 -0.08 8.50 16 6 O -0.56 -10.87 13.97 5.56 0.08 -10.79 16 7 N -0.70 -11.94 5.13 -60.29 -0.31 -12.24 16 8 C 0.12 2.13 4.36 -4.04 -0.02 2.12 16 9 C 0.19 3.35 0.41 -131.64 -0.05 3.30 16 10 N -0.71 -15.20 3.60 -47.06 -0.17 -15.37 16 11 C -0.23 -5.79 8.67 -15.06 -0.13 -5.92 16 12 H 0.07 1.76 6.78 -52.49 -0.36 1.41 16 13 C -0.01 -0.16 5.50 -17.43 -0.10 -0.26 16 14 N -0.93 -25.81 13.07 55.49 0.73 -25.09 16 15 Si 0.90 21.51 29.55 -169.99 -5.02 16.49 16 16 H -0.28 -6.61 7.11 56.52 0.40 -6.21 16 17 H -0.29 -7.08 7.11 56.52 0.40 -6.68 16 18 H -0.29 -7.12 7.11 56.52 0.40 -6.72 16 19 C -0.15 -2.32 4.80 -26.59 -0.13 -2.45 16 20 C 0.07 1.05 6.67 37.21 0.25 1.30 16 21 O -0.39 -6.93 10.77 -35.23 -0.38 -7.31 16 22 C 0.07 1.23 6.68 37.21 0.25 1.48 16 23 C -0.15 -2.66 4.78 -26.59 -0.13 -2.79 16 24 C 0.08 0.89 3.21 -26.01 -0.08 0.80 16 25 H 0.08 0.85 8.14 -51.93 -0.42 0.43 16 26 C -0.15 -1.31 6.36 -24.79 -0.16 -1.46 16 27 C -0.16 -1.14 7.16 -24.48 -0.18 -1.31 16 28 C 0.05 0.44 7.71 38.19 0.29 0.73 16 29 O -0.37 -4.32 9.82 -35.23 -0.35 -4.66 16 30 H 0.05 0.65 7.54 -51.93 -0.39 0.26 16 31 H 0.08 1.14 5.95 -51.93 -0.31 0.83 16 32 H 0.05 0.51 6.35 -51.93 -0.33 0.18 16 33 H 0.09 1.20 8.14 -51.93 -0.42 0.78 16 34 H 0.09 1.21 8.01 -51.92 -0.42 0.80 16 35 H 0.40 6.68 7.43 -40.82 -0.30 6.38 16 36 H 0.06 1.29 7.47 -51.93 -0.39 0.90 16 37 H 0.05 0.91 8.14 -51.93 -0.42 0.49 16 38 H 0.38 8.81 6.82 -40.82 -0.28 8.53 16 39 H 0.25 7.00 8.31 -40.82 -0.34 6.66 16 40 H 0.07 1.00 8.14 -51.93 -0.42 0.58 16 41 H 0.07 1.01 7.40 -51.93 -0.38 0.62 16 42 H 0.06 1.06 7.15 -51.93 -0.37 0.69 16 43 H 0.08 1.10 8.14 -51.93 -0.42 0.67 16 44 H 0.07 1.42 7.75 -51.93 -0.40 1.01 16 45 H 0.08 1.21 8.14 -51.93 -0.42 0.79 16 46 H 0.08 1.43 7.20 -51.93 -0.37 1.06 16 47 H 0.05 0.87 8.14 -51.93 -0.42 0.45 16 48 H 0.08 0.71 6.07 -51.93 -0.32 0.40 16 49 H 0.08 0.61 8.14 -51.92 -0.42 0.18 16 50 H 0.08 0.53 8.14 -51.93 -0.42 0.10 16 51 H 0.08 0.47 8.14 -51.93 -0.42 0.04 16 52 H 0.08 0.60 8.14 -51.93 -0.42 0.18 16 53 H 0.06 0.48 8.14 -51.93 -0.42 0.06 16 LS Contribution 402.35 15.07 6.06 6.06 Total: -1.00 -29.08 402.35 -8.86 -37.94 By element: Atomic # 1 Polarization: 24.84 SS G_CDS: -9.24 Total: 15.60 kcal Atomic # 6 Polarization: -0.37 SS G_CDS: -0.25 Total: -0.63 kcal Atomic # 7 Polarization: -52.95 SS G_CDS: 0.25 Total: -52.70 kcal Atomic # 8 Polarization: -22.11 SS G_CDS: -0.65 Total: -22.76 kcal Atomic # 14 Polarization: 21.51 SS G_CDS: -5.02 Total: 16.49 kcal Total LS contribution 6.06 Total: 6.06 kcal Total: -29.08 -8.86 -37.94 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. ZINC001484190882.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 -95.031 kcal (2) G-P(sol) polarization free energy of solvation -29.082 kcal (3) G-ENP(sol) elect.-nuc.-pol. free energy of system -124.113 kcal (4) G-CDS(sol) cavity-dispersion-solvent structure free energy -8.857 kcal (5) G-P-CDS(sol) = G-P(sol) + G-CDS(sol) = (2) + (4) -37.938 kcal (6) G-S(sol) free energy of system = (1) + (5) -132.970 kcal 1SCF WAS SPECIFIED, SO GEOMETRY OPTIMIZATION WAS NOT USED Total computer time = 1.09 seconds