Wall clock time and date at job start Tue Mar 31 2020 23:47:31 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.08998 * 109.47983 * 2 1 4 4 C 1.53088 * 109.47644 * 239.97588 * 2 1 3 5 5 C 1.53178 * 109.15977 * 182.40760 * 4 2 1 6 6 H 1.08999 * 109.54579 * 176.85056 * 5 4 2 7 7 C 1.50697 * 109.54511 * 297.02551 * 5 4 2 8 8 O 1.21286 * 120.00061 * 9.22123 * 7 5 4 9 9 N 1.34777 * 120.00077 * 189.21797 * 7 5 4 10 10 C 1.46924 * 120.63131 * 5.46559 * 9 7 5 11 11 C 1.52763 * 109.00823 * 233.58302 * 10 9 7 12 12 C 1.53043 * 109.38345 * 305.58817 * 11 10 9 13 13 N 1.46505 * 109.44952 * 181.25382 * 12 11 10 14 14 C 1.36938 * 119.99986 * 299.93772 * 13 12 11 15 15 H 1.07997 * 119.99644 * 319.94783 * 14 13 12 16 16 C 1.38806 * 120.00063 * 139.94403 * 14 13 12 17 17 N 1.31614 * 120.00176 * 343.18844 * 16 14 13 18 18 Si 1.86796 * 120.00029 * 163.18885 * 16 14 13 19 19 H 1.48501 * 109.47378 * 60.00059 * 18 16 14 20 20 H 1.48502 * 109.47442 * 179.97438 * 18 16 14 21 21 H 1.48500 * 109.47369 * 299.99990 * 18 16 14 22 22 C 1.46493 * 119.99887 * 119.93612 * 13 12 11 23 23 C 1.52997 * 109.47207 * 84.54903 * 22 13 12 24 24 O 1.42898 * 109.47479 * 64.99734 * 23 22 13 25 25 C 1.53043 * 109.53850 * 61.24068 * 12 11 10 26 26 C 1.46931 * 120.63099 * 185.74293 * 9 7 5 27 27 C 1.53177 * 109.03858 * 56.98424 * 5 4 2 28 28 C 1.53088 * 109.16096 * 303.02023 * 27 5 4 29 29 O 1.42903 * 109.42759 * 57.60096 * 28 27 5 30 30 H 1.09004 * 109.47213 * 180.02562 * 1 2 3 31 31 H 1.09002 * 109.47054 * 300.00065 * 1 2 3 32 32 H 1.09002 * 109.47529 * 59.99769 * 1 2 3 33 33 H 1.08998 * 109.52725 * 302.31250 * 4 2 1 34 34 H 1.09001 * 109.52357 * 62.50351 * 4 2 1 35 35 H 1.09002 * 109.58668 * 353.41339 * 10 9 7 36 36 H 1.08999 * 109.59001 * 113.84295 * 10 9 7 37 37 H 1.09001 * 109.56474 * 185.55069 * 11 10 9 38 38 H 1.09005 * 109.36332 * 65.49977 * 11 10 9 39 39 H 1.09005 * 109.45581 * 301.21636 * 12 11 10 40 40 H 0.97001 * 120.00353 * 179.97438 * 17 16 14 41 41 H 1.09005 * 109.47186 * 204.55017 * 22 13 12 42 42 H 1.09001 * 109.47688 * 324.54290 * 22 13 12 43 43 H 1.08999 * 109.47166 * 184.99928 * 23 22 13 44 44 H 1.08999 * 109.47161 * 304.99894 * 23 22 13 45 45 H 0.96694 * 114.00664 * 180.02562 * 24 23 22 46 46 H 1.08998 * 109.49348 * 178.68776 * 25 12 11 47 47 H 1.09000 * 109.52415 * 58.60048 * 25 12 11 48 48 H 1.09002 * 109.58692 * 246.15896 * 26 9 7 49 49 H 1.08999 * 109.58314 * 6.44162 * 26 9 7 50 50 H 1.09001 * 109.51980 * 62.92057 * 27 5 4 51 51 H 1.09000 * 109.52513 * 183.11155 * 27 5 4 52 52 H 1.09008 * 109.48240 * 177.59060 * 28 27 5 53 53 H 1.09003 * 109.48210 * 297.61465 * 28 27 5 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.8935 1.0276 0.0000 4 6 2.0404 -0.7222 -1.2496 5 6 3.5710 -0.7708 -1.2122 6 1 3.9391 -1.3324 -2.0708 7 6 4.1211 0.6316 -1.2539 8 8 3.3737 1.5799 -1.1396 9 7 5.4437 0.8320 -1.4182 10 6 6.3457 -0.2984 -1.6773 11 6 7.4795 -0.2750 -0.6537 12 6 8.1625 1.0942 -0.6835 13 7 9.2692 1.1083 0.2764 14 6 9.0251 0.8901 1.6061 15 1 8.1516 1.3212 2.0725 16 6 9.9002 0.1149 2.3543 17 7 10.7631 -0.6734 1.7492 18 14 9.8531 0.1903 4.2202 19 1 8.5093 -0.2226 4.6986 20 1 10.8798 -0.7257 4.7788 21 1 10.1324 1.5788 4.6667 22 6 10.6370 1.3558 -0.1862 23 6 10.8789 2.8636 -0.2793 24 8 10.8120 3.4384 1.0273 25 6 7.1501 2.1795 -0.3103 26 6 6.0119 2.1846 -1.3357 27 6 4.0176 -1.4601 0.0807 28 6 3.4307 -0.7095 1.2789 29 8 2.0066 -0.6741 1.1665 30 1 -0.3634 -1.0277 0.0005 31 1 -0.3633 0.5139 0.8900 32 1 -0.3634 0.5139 -0.8900 33 1 1.6442 -1.7374 -1.2711 34 1 1.7147 -0.1841 -2.1399 35 1 5.7939 -1.2342 -1.5878 36 1 6.7616 -0.2116 -2.6810 37 1 8.2067 -1.0498 -0.8967 38 1 7.0722 -0.4568 0.3410 39 1 8.5480 1.2859 -1.6850 40 1 11.3749 -1.2149 2.2721 41 1 11.3438 0.9171 0.5182 42 1 10.7757 0.9042 -1.1686 43 1 11.8644 3.0477 -0.7071 44 1 10.1169 3.3154 -0.9143 45 1 10.9560 4.3944 1.0444 46 1 7.6429 3.1517 -0.3091 47 1 6.7458 1.9752 0.6811 48 1 6.3999 2.4781 -2.3111 49 1 5.2391 2.8870 -1.0234 50 1 3.6612 -2.4902 0.0860 51 1 5.1058 -1.4499 0.1418 52 1 3.7100 -1.2206 2.2004 53 1 3.8194 0.3087 1.2959 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 ZINC001098512185.mol2 53 1SCF WAS SPECIFIED, SO GEOMETRY OPTIMIZATION WAS NOT USED SCF FIELD WAS ACHIEVED Wall clock time and date at job start Tue Mar 31 2020 23:47:31 Heat of formation + Delta-G solvation = -112.278242 kcal Electronic energy + Delta-G solvation = -32784.369023 eV Core-core repulsion = 28656.538014 eV Total energy + Delta-G solvation = -4127.831009 eV No. of doubly occupied orbitals = 66 Molecular weight (most abundant/longest-lived isotopes) = 340.220 amu Computer time = 1.25 seconds Orbital eigenvalues (eV) -40.84674 -38.69564 -37.38668 -36.17092 -34.68682 -33.81940 -32.74362 -31.89327 -30.62299 -30.35393 -27.95027 -27.67262 -25.35310 -24.46484 -23.32042 -22.72481 -22.33148 -21.34020 -20.36926 -19.77454 -18.45014 -16.98393 -16.59144 -16.28241 -16.21707 -16.02673 -15.51651 -15.19869 -14.89133 -14.50304 -14.44648 -13.96909 -13.81780 -13.71729 -13.42391 -13.28308 -13.09924 -12.67129 -12.54520 -12.43213 -12.35031 -12.11939 -12.01901 -11.83457 -11.71618 -11.51764 -11.20247 -11.05669 -10.91708 -10.78784 -10.70874 -10.68593 -10.53880 -10.41643 -10.05780 -9.91159 -9.87834 -9.62657 -9.50874 -9.11093 -9.03974 -8.40211 -8.37531 -6.48887 -5.72745 -3.31792 2.70749 3.17325 3.42347 3.47206 3.53017 3.87463 4.46674 4.53446 4.62205 4.72898 4.91649 4.99166 5.13224 5.17128 5.26218 5.29713 5.38688 5.41872 5.50486 5.51483 5.54828 5.56900 5.61295 5.69855 5.72823 5.76545 5.80644 5.85682 5.94174 5.99528 6.00554 6.00718 6.05735 6.12965 6.16812 6.24245 6.30929 6.37684 6.45905 6.61201 6.67040 6.79061 6.82492 7.07904 7.15752 7.36691 7.59851 7.65644 7.88865 8.30792 8.45028 8.64165 9.24891 9.86874 10.43894 11.29701 Molecular weight = 340.22amu Principal moments of inertia in cm(-1) A = 0.014043 B = 0.003640 C = 0.003236 Principal moments of inertia in units of 10**(-40)*gram-cm**2 A = 1993.381948 B = 7690.867678 C = 8651.550408 Net atomic charges, atomic populations, and dipole contributions calculated with CM2 Atom NO. Type Charge No. of electrons 1 C -0.141 4.141 2 C 0.069 3.931 3 H 0.093 0.907 4 C -0.125 4.125 5 C -0.111 4.111 6 H 0.096 0.904 7 C 0.519 3.481 8 O -0.548 6.548 9 N -0.601 5.601 10 C 0.101 3.899 11 C -0.135 4.135 12 C 0.179 3.821 13 N -0.571 5.571 14 C -0.234 4.234 15 H 0.053 0.947 16 C 0.018 3.982 17 N -0.896 5.896 18 Si 0.895 3.105 19 H -0.284 1.284 20 H -0.293 1.293 21 H -0.282 1.282 22 C 0.147 3.853 23 C 0.058 3.942 24 O -0.544 6.544 25 C -0.138 4.138 26 C 0.114 3.886 27 C -0.151 4.151 28 C 0.062 3.938 29 O -0.384 6.384 30 H 0.057 0.943 31 H 0.060 0.940 32 H 0.067 0.933 33 H 0.073 0.927 34 H 0.075 0.925 35 H 0.079 0.921 36 H 0.063 0.937 37 H 0.073 0.927 38 H 0.072 0.928 39 H 0.045 0.955 40 H 0.252 0.748 41 H 0.102 0.898 42 H 0.035 0.965 43 H 0.019 0.981 44 H 0.021 0.979 45 H 0.349 0.651 46 H 0.065 0.935 47 H 0.067 0.933 48 H 0.059 0.941 49 H 0.083 0.917 50 H 0.081 0.919 51 H 0.089 0.911 52 H 0.096 0.904 53 H 0.051 0.949 For the charges calculated by CM2: Dipole moment (debyes) X Y Z Total from point charges -11.103 0.282 -9.294 14.482 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.198 4.198 2 C 0.012 3.988 3 H 0.110 0.890 4 C -0.163 4.163 5 C -0.133 4.133 6 H 0.114 0.886 7 C 0.307 3.693 8 O -0.427 6.427 9 N -0.333 5.333 10 C -0.022 4.022 11 C -0.175 4.175 12 C 0.071 3.929 13 N -0.290 5.290 14 C -0.364 4.364 15 H 0.071 0.929 16 C -0.181 4.181 17 N -0.540 5.540 18 Si 0.725 3.275 19 H -0.211 1.211 20 H -0.219 1.219 21 H -0.208 1.208 22 C 0.020 3.980 23 C -0.022 4.022 24 O -0.345 6.345 25 C -0.178 4.178 26 C -0.009 4.009 27 C -0.188 4.188 28 C -0.014 4.014 29 O -0.303 6.303 30 H 0.076 0.924 31 H 0.079 0.921 32 H 0.086 0.914 33 H 0.092 0.908 34 H 0.093 0.907 35 H 0.098 0.902 36 H 0.081 0.919 37 H 0.092 0.908 38 H 0.091 0.909 39 H 0.063 0.937 40 H 0.062 0.938 41 H 0.119 0.881 42 H 0.053 0.947 43 H 0.037 0.963 44 H 0.039 0.961 45 H 0.192 0.808 46 H 0.084 0.916 47 H 0.086 0.914 48 H 0.077 0.923 49 H 0.101 0.899 50 H 0.100 0.900 51 H 0.108 0.892 52 H 0.115 0.885 53 H 0.069 0.931 Dipole moment (debyes) X Y Z Total from point charges -11.575 -0.071 -8.867 14.581 hybrid contribution 0.653 0.657 0.501 1.053 sum -10.922 0.586 -8.365 13.770 Atomic orbital electron populations 1.21734 0.93764 1.02172 1.02160 1.22480 0.94490 0.97646 0.84175 0.88973 1.21992 0.93023 1.01547 0.99765 1.21392 0.97261 0.94213 1.00393 0.88575 1.20545 0.81883 0.91129 0.75706 1.90572 1.57146 1.44064 1.50939 1.48154 1.06426 1.08825 1.69903 1.21942 0.92178 0.88418 0.99657 1.22249 0.95949 1.00586 0.98703 1.20296 0.88758 0.91523 0.92336 1.43652 1.00626 1.81279 1.03404 1.18683 1.09567 1.23853 0.84301 0.92946 1.25111 0.99713 0.95477 0.97777 1.69911 1.23663 1.26077 1.34320 0.86509 0.77075 0.77513 0.86391 1.21065 1.21890 1.20845 1.20980 0.89733 0.88876 0.98403 1.23114 0.95997 0.99297 0.83813 1.86849 1.95223 1.23500 1.28956 1.22353 0.95826 1.00771 0.98869 1.21672 0.96694 0.83123 0.99412 1.22293 1.02625 1.00232 0.93690 1.22512 0.81002 0.99877 0.98054 1.87867 1.20828 1.78222 1.43421 0.92409 0.92122 0.91367 0.90835 0.90681 0.90246 0.91877 0.90814 0.90903 0.93674 0.93771 0.88054 0.94677 0.96326 0.96058 0.80788 0.91594 0.91392 0.92283 0.89897 0.90017 0.89216 0.88542 0.93082 Number of geometries 1 Number of calculations of the screened coulomb radii 1 The total number of SCF iterations 15. 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.14 -1.53 9.56 37.16 0.36 -1.17 16 2 C 0.07 0.92 2.98 -26.68 -0.08 0.84 16 3 H 0.09 1.45 5.49 -51.93 -0.29 1.17 16 4 C -0.13 -1.45 4.59 -26.59 -0.12 -1.57 16 5 C -0.11 -1.33 2.16 -91.59 -0.20 -1.53 16 6 H 0.10 0.96 7.19 -51.93 -0.37 0.59 16 7 C 0.52 7.63 5.74 -10.99 -0.06 7.56 16 8 O -0.55 -9.44 13.03 5.55 0.07 -9.36 16 9 N -0.60 -8.63 2.98 -173.03 -0.52 -9.15 16 10 C 0.10 1.33 5.49 -3.93 -0.02 1.31 16 11 C -0.13 -2.31 5.16 -26.83 -0.14 -2.44 16 12 C 0.18 3.31 2.60 -67.88 -0.18 3.13 16 13 N -0.57 -12.92 2.31 -165.61 -0.38 -13.31 16 14 C -0.23 -6.03 8.00 -15.06 -0.12 -6.15 16 15 H 0.05 1.41 6.90 -52.49 -0.36 1.04 16 16 C 0.02 0.50 5.00 -17.43 -0.09 0.41 16 17 N -0.90 -24.82 11.45 55.49 0.64 -24.19 16 18 Si 0.90 21.82 29.59 -169.99 -5.03 16.79 16 19 H -0.28 -6.84 7.11 56.52 0.40 -6.44 16 20 H -0.29 -7.18 7.11 56.52 0.40 -6.78 16 21 H -0.28 -7.02 7.11 56.52 0.40 -6.62 16 22 C 0.15 3.24 5.55 -4.05 -0.02 3.22 16 23 C 0.06 1.10 7.39 37.16 0.27 1.37 16 24 O -0.54 -11.13 13.39 -35.23 -0.47 -11.60 16 25 C -0.14 -2.40 5.30 -26.61 -0.14 -2.54 16 26 C 0.11 1.70 6.42 -3.71 -0.02 1.68 16 27 C -0.15 -1.83 5.25 -26.59 -0.14 -1.97 16 28 C 0.06 0.87 6.62 37.20 0.25 1.12 16 29 O -0.38 -5.64 10.41 -35.23 -0.37 -6.01 16 30 H 0.06 0.58 8.14 -51.93 -0.42 0.15 16 31 H 0.06 0.65 8.14 -51.93 -0.42 0.23 16 32 H 0.07 0.64 8.14 -51.93 -0.42 0.22 16 33 H 0.07 0.75 8.14 -51.93 -0.42 0.33 16 34 H 0.07 0.84 8.13 -51.93 -0.42 0.42 16 35 H 0.08 0.90 4.66 -51.93 -0.24 0.66 16 36 H 0.06 0.75 8.14 -51.93 -0.42 0.33 16 37 H 0.07 1.36 8.14 -51.93 -0.42 0.94 16 38 H 0.07 1.42 7.11 -51.93 -0.37 1.05 16 39 H 0.05 0.79 7.87 -51.93 -0.41 0.38 16 40 H 0.25 6.89 8.31 -40.82 -0.34 6.55 16 41 H 0.10 2.64 5.61 -51.93 -0.29 2.35 16 42 H 0.03 0.72 8.04 -51.93 -0.42 0.30 16 43 H 0.02 0.31 8.14 -51.93 -0.42 -0.11 16 44 H 0.02 0.37 8.14 -51.93 -0.42 -0.05 16 45 H 0.35 5.20 9.12 45.56 0.42 5.62 16 46 H 0.07 1.17 8.14 -51.93 -0.42 0.74 16 47 H 0.07 1.35 7.08 -51.93 -0.37 0.98 16 48 H 0.06 0.80 8.14 -51.93 -0.42 0.38 16 49 H 0.08 1.31 7.04 -51.93 -0.37 0.94 16 50 H 0.08 0.88 8.14 -51.93 -0.42 0.46 16 51 H 0.09 1.14 5.47 -51.93 -0.28 0.86 16 52 H 0.10 1.26 8.14 -51.92 -0.42 0.83 16 53 H 0.05 0.81 7.38 -51.93 -0.38 0.43 16 LS Contribution 395.36 15.07 5.96 5.96 Total: -1.00 -30.72 395.36 -8.92 -39.64 By element: Atomic # 1 Polarization: 16.33 SS G_CDS: -8.36 Total: 7.96 kcal Atomic # 6 Polarization: 3.73 SS G_CDS: -0.46 Total: 3.27 kcal Atomic # 7 Polarization: -46.38 SS G_CDS: -0.26 Total: -46.64 kcal Atomic # 8 Polarization: -26.21 SS G_CDS: -0.77 Total: -26.98 kcal Atomic # 14 Polarization: 21.82 SS G_CDS: -5.03 Total: 16.79 kcal Total LS contribution 5.96 Total: 5.96 kcal Total: -30.72 -8.92 -39.64 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. ZINC001098512185.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 -72.637 kcal (2) G-P(sol) polarization free energy of solvation -30.720 kcal (3) G-ENP(sol) elect.-nuc.-pol. free energy of system -103.357 kcal (4) G-CDS(sol) cavity-dispersion-solvent structure free energy -8.921 kcal (5) G-P-CDS(sol) = G-P(sol) + G-CDS(sol) = (2) + (4) -39.641 kcal (6) G-S(sol) free energy of system = (1) + (5) -112.278 kcal 1SCF WAS SPECIFIED, SO GEOMETRY OPTIMIZATION WAS NOT USED Total computer time = 1.26 seconds