Wall clock time and date at job start Wed Apr 1 2020 00:29:59 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 O 1.42896 * 1 3 3 C 1.42907 * 114.00009 * 2 1 4 4 C 1.53003 * 109.47082 * 180.02562 * 3 2 1 5 5 C 1.50705 * 110.74940 * 297.98004 * 4 3 2 6 6 O 1.21281 * 119.99617 * 242.01876 * 5 4 3 7 7 N 1.34774 * 119.99934 * 62.01941 * 5 4 3 8 8 C 1.46498 * 120.00123 * 179.97438 * 7 5 4 9 9 C 1.52991 * 109.47441 * 180.02562 * 8 7 5 10 10 C 1.53007 * 109.47232 * 179.97438 * 9 8 7 11 11 H 1.09007 * 110.88288 * 54.30204 * 10 9 8 12 12 C 1.54909 * 111.00563 * 290.41941 * 10 9 8 13 13 N 1.47465 * 104.80749 * 204.94896 * 12 10 9 14 14 C 1.36943 * 125.62438 * 204.04442 * 13 12 10 15 15 H 1.07993 * 119.99912 * 0.02562 * 14 13 12 16 16 C 1.38814 * 119.99761 * 180.02562 * 14 13 12 17 17 N 1.31617 * 120.00063 * 0.02562 * 16 14 13 18 18 Si 1.86798 * 119.99946 * 180.02562 * 16 14 13 19 19 H 1.48502 * 109.47229 * 59.99872 * 18 16 14 20 20 H 1.48500 * 109.47163 * 179.97438 * 18 16 14 21 21 H 1.48504 * 109.46845 * 300.00131 * 18 16 14 22 22 C 1.47012 * 108.72808 * 24.07173 * 13 12 10 23 23 C 1.54332 * 107.27183 * 358.97198 * 22 13 12 24 24 H 1.09001 * 110.71981 * 96.17329 * 23 22 13 25 25 C 1.52997 * 110.72196 * 219.44641 * 23 22 13 26 26 C 1.55566 * 111.01264 * 174.40783 * 4 3 2 27 27 C 1.54944 * 100.82489 * 207.60464 * 26 4 3 28 28 O 1.43645 * 105.22260 * 36.43849 * 27 26 4 29 29 C 1.43235 * 109.97529 * 335.52671 * 28 27 26 30 30 H 1.09005 * 109.47199 * 59.99864 * 1 2 3 31 31 H 1.09000 * 109.47340 * 179.97438 * 1 2 3 32 32 H 1.09005 * 109.47455 * 300.00164 * 1 2 3 33 33 H 1.08993 * 109.46866 * 59.99799 * 3 2 1 34 34 H 1.09005 * 109.46780 * 299.99587 * 3 2 1 35 35 H 0.97001 * 120.00393 * 359.97438 * 7 5 4 36 36 H 1.09010 * 109.47124 * 300.00436 * 8 7 5 37 37 H 1.09006 * 109.47111 * 59.99703 * 8 7 5 38 38 H 1.09003 * 109.47351 * 299.99856 * 9 8 7 39 39 H 1.08994 * 109.47366 * 60.00124 * 9 8 7 40 40 H 1.09006 * 110.36261 * 323.75429 * 12 10 9 41 41 H 1.08992 * 110.36592 * 86.07320 * 12 10 9 42 42 H 0.97004 * 119.99878 * 179.97438 * 17 16 14 43 43 H 1.09005 * 109.88678 * 239.56096 * 22 13 12 44 44 H 1.08996 * 110.01854 * 118.46167 * 22 13 12 45 45 H 1.08998 * 109.47184 * 173.51915 * 25 23 22 46 46 H 1.08999 * 109.47151 * 293.52222 * 25 23 22 47 47 H 1.09001 * 109.47310 * 53.52045 * 25 23 22 48 48 H 1.08997 * 111.22347 * 89.69902 * 26 4 3 49 49 H 1.09003 * 111.13190 * 325.46499 * 26 4 3 50 50 H 1.09002 * 110.30953 * 277.47426 * 27 26 4 51 51 H 1.09000 * 110.30834 * 155.40499 * 27 26 4 52 52 H 1.08995 * 109.83514 * 121.22466 * 29 28 27 53 53 H 1.09003 * 109.82666 * 242.20160 * 29 28 27 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.0102 1.3055 0.0000 4 6 3.5355 1.1847 -0.0006 5 6 4.0160 0.4832 1.2436 6 8 4.7400 1.0627 2.0251 7 7 3.6412 -0.7882 1.4878 8 6 4.1087 -1.4703 2.6970 9 6 3.5431 -2.8914 2.7320 10 6 4.0319 -3.6040 3.9946 11 1 3.7932 -3.0233 4.8857 12 6 5.5508 -3.8978 3.9136 13 7 5.7595 -5.0520 4.8074 14 6 6.9385 -5.3831 5.4204 15 1 7.8173 -4.7750 5.2651 16 6 7.0039 -6.5006 6.2414 17 7 5.9330 -7.2420 6.4302 18 14 8.6118 -6.9519 7.0782 19 1 9.0430 -5.8352 7.9571 20 1 8.4181 -8.1772 7.8945 21 1 9.6525 -7.2017 6.0486 22 6 4.4932 -5.7885 4.9309 23 6 3.4413 -5.0357 4.0891 24 1 3.3493 -5.4811 3.0986 25 6 2.0886 -5.0171 4.8037 26 6 4.2056 2.5815 -0.1422 27 6 5.5423 2.1952 -0.8240 28 8 5.2124 1.1217 -1.7197 29 6 4.0282 0.4542 -1.2683 30 1 -0.3634 0.5139 -0.8900 31 1 -0.3634 -1.0276 -0.0005 32 1 -0.3634 0.5139 0.8900 33 1 1.6886 1.8464 0.8899 34 1 1.6886 1.8464 -0.8901 35 1 3.0618 -1.2516 0.8629 36 1 5.1980 -1.5132 2.6918 37 1 3.7710 -0.9226 3.5768 38 1 2.4539 -2.8489 2.7370 39 1 3.8812 -3.4389 1.8523 40 1 5.8338 -4.1530 2.8924 41 1 6.1225 -3.0382 4.2630 42 1 5.9788 -8.0231 7.0037 43 1 4.1811 -5.8170 5.9749 44 1 4.6135 -6.8035 4.5523 45 1 1.3938 -4.3856 4.2500 46 1 1.6929 -6.0313 4.8588 47 1 2.2153 -4.6213 5.8114 48 1 4.3790 3.0365 0.8330 49 1 3.6144 3.2381 -0.7805 50 1 6.2634 1.8583 -0.0792 51 1 5.9420 3.0429 -1.3804 52 1 3.2615 0.4977 -2.0418 53 1 4.2572 -0.5847 -1.0309 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 ZINC001103276446.mol2 53 1SCF WAS SPECIFIED, SO GEOMETRY OPTIMIZATION WAS NOT USED SCF FIELD WAS ACHIEVED Wall clock time and date at job start Wed Apr 1 2020 00:29:59 Heat of formation + Delta-G solvation = -113.361729 kcal Electronic energy + Delta-G solvation = -31079.373683 eV Core-core repulsion = 26951.495691 eV Total energy + Delta-G solvation = -4127.877992 eV No. of doubly occupied orbitals = 66 Molecular weight (most abundant/longest-lived isotopes) = 340.220 amu Computer time = 1.13 seconds Orbital eigenvalues (eV) -41.06934 -38.66732 -38.11437 -37.99388 -35.22669 -33.51920 -33.12419 -31.16617 -30.61701 -29.91383 -27.51441 -26.86840 -25.68062 -25.64249 -24.32804 -22.91497 -22.68645 -21.78202 -20.11962 -19.54089 -19.14578 -17.94197 -16.94668 -16.74366 -16.73011 -16.14627 -15.91526 -15.53087 -15.41398 -15.18820 -14.80535 -14.60705 -14.21439 -13.87181 -13.65709 -13.51705 -13.43780 -13.01324 -12.87028 -12.52943 -12.29208 -12.16258 -12.08644 -11.83538 -11.80790 -11.67475 -11.60319 -11.35373 -11.28685 -11.17467 -10.99786 -10.73738 -10.70212 -10.29096 -10.05605 -9.96525 -9.84664 -9.81374 -9.74554 -9.58838 -9.42895 -9.01802 -8.53654 -6.75337 -5.64769 -2.97586 2.42012 2.97456 3.16047 3.45007 3.48866 3.50223 3.55630 3.86754 4.30480 4.43614 4.54394 4.56047 4.73728 4.81036 4.88970 4.91242 4.92747 4.98787 5.06821 5.16490 5.27024 5.32516 5.36908 5.37767 5.44702 5.51299 5.51772 5.56138 5.65445 5.69291 5.73832 5.80022 5.98794 6.04605 6.14730 6.26780 6.30662 6.40744 6.46369 6.52741 6.56809 6.62808 6.68216 6.71406 6.85034 7.07008 7.48563 7.63716 7.71059 7.90290 8.14778 8.33299 9.29371 9.89787 10.56561 11.46849 Molecular weight = 340.22amu Principal moments of inertia in cm(-1) A = 0.015236 B = 0.002450 C = 0.002195 Principal moments of inertia in units of 10**(-40)*gram-cm**2 A = 1837.303711 B =11424.528239 C =12753.677097 Net atomic charges, atomic populations, and dipole contributions calculated with CM2 Atom NO. Type Charge No. of electrons 1 C 0.034 3.966 2 O -0.391 6.391 3 C 0.082 3.918 4 C -0.117 4.117 5 C 0.528 3.472 6 O -0.539 6.539 7 N -0.714 5.714 8 C 0.118 3.882 9 C -0.114 4.114 10 C -0.095 4.095 11 H 0.075 0.925 12 C 0.139 3.861 13 N -0.607 5.607 14 C -0.226 4.226 15 H 0.076 0.924 16 C -0.010 4.010 17 N -0.923 5.923 18 Si 0.911 3.089 19 H -0.285 1.285 20 H -0.291 1.291 21 H -0.285 1.285 22 C 0.184 3.816 23 C -0.114 4.114 24 H 0.064 0.936 25 C -0.132 4.132 26 C -0.131 4.131 27 C 0.052 3.948 28 O -0.382 6.382 29 C 0.072 3.928 30 H 0.043 0.957 31 H 0.091 0.909 32 H 0.043 0.957 33 H 0.061 0.939 34 H 0.066 0.934 35 H 0.407 0.593 36 H 0.072 0.928 37 H 0.066 0.934 38 H 0.065 0.935 39 H 0.068 0.932 40 H 0.017 0.983 41 H 0.037 0.963 42 H 0.252 0.748 43 H 0.047 0.953 44 H 0.039 0.961 45 H 0.042 0.958 46 H 0.051 0.949 47 H 0.053 0.947 48 H 0.103 0.897 49 H 0.086 0.914 50 H 0.070 0.930 51 H 0.090 0.910 52 H 0.082 0.918 53 H 0.075 0.925 For the charges calculated by CM2: Dipole moment (debyes) X Y Z Total from point charges -13.929 18.529 -16.437 28.417 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.062 4.062 2 O -0.310 6.310 3 C 0.005 3.995 4 C -0.120 4.120 5 C 0.312 3.688 6 O -0.416 6.416 7 N -0.365 5.365 8 C -0.005 4.005 9 C -0.153 4.153 10 C -0.115 4.115 11 H 0.094 0.906 12 C 0.013 3.987 13 N -0.333 5.333 14 C -0.355 4.355 15 H 0.094 0.906 16 C -0.206 4.206 17 N -0.571 5.571 18 Si 0.742 3.258 19 H -0.212 1.212 20 H -0.216 1.216 21 H -0.212 1.212 22 C 0.059 3.941 23 C -0.134 4.134 24 H 0.082 0.918 25 C -0.190 4.190 26 C -0.169 4.169 27 C -0.025 4.025 28 O -0.301 6.301 29 C -0.005 4.005 30 H 0.062 0.938 31 H 0.109 0.891 32 H 0.062 0.938 33 H 0.079 0.921 34 H 0.084 0.916 35 H 0.243 0.757 36 H 0.090 0.910 37 H 0.084 0.916 38 H 0.084 0.916 39 H 0.086 0.914 40 H 0.035 0.965 41 H 0.055 0.945 42 H 0.062 0.938 43 H 0.065 0.935 44 H 0.057 0.943 45 H 0.061 0.939 46 H 0.070 0.930 47 H 0.072 0.928 48 H 0.121 0.879 49 H 0.105 0.895 50 H 0.088 0.912 51 H 0.108 0.892 52 H 0.100 0.900 53 H 0.094 0.906 Dipole moment (debyes) X Y Z Total from point charges -13.094 18.511 -15.879 27.681 hybrid contribution 0.148 -0.299 0.433 0.546 sum -12.945 18.211 -15.447 27.163 Atomic orbital electron populations 1.22777 0.79937 1.03366 1.00162 1.87759 1.20580 1.27225 1.95472 1.22368 0.89959 0.85573 1.01627 1.22504 0.95799 0.98011 0.95710 1.19870 0.80424 0.82095 0.86401 1.90707 1.36868 1.65779 1.48244 1.46072 1.51290 1.13475 1.25709 1.21504 0.99981 0.92573 0.86485 1.21722 1.00398 0.93754 0.99399 1.22326 0.97268 0.95742 0.96127 0.90629 1.21623 0.92423 0.90945 0.93700 1.44472 1.07846 1.27305 1.53666 1.19035 0.87321 1.07822 1.21289 0.90571 1.24788 0.97499 0.97883 1.00397 1.69960 1.26923 1.18285 1.41944 0.86210 0.82859 0.77549 0.79225 1.21211 1.21647 1.21217 1.21196 0.82621 0.93289 0.96991 1.22480 0.95815 0.95454 0.99617 0.91784 1.21669 0.96289 1.00874 1.00153 1.23226 0.96181 0.93977 1.03530 1.23419 0.96434 0.89877 0.92731 1.88412 1.44194 1.47446 1.50065 1.23115 0.87696 0.96417 0.93280 0.93836 0.89051 0.93812 0.92103 0.91597 0.75664 0.91000 0.91580 0.91631 0.91372 0.96523 0.94499 0.93819 0.93525 0.94286 0.93867 0.93034 0.92801 0.87894 0.89504 0.91216 0.89227 0.89998 0.90641 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.03 0.17 11.01 101.05 1.11 1.28 16 2 O -0.39 -2.83 8.31 -35.23 -0.29 -3.12 16 3 C 0.08 0.55 4.74 37.16 0.18 0.73 16 4 C -0.12 -1.05 0.63 -153.57 -0.10 -1.14 16 5 C 0.53 6.57 6.31 -10.98 -0.07 6.50 16 6 O -0.54 -8.60 15.62 5.56 0.09 -8.51 16 7 N -0.71 -8.30 4.96 -60.29 -0.30 -8.60 16 8 C 0.12 1.57 5.55 -4.05 -0.02 1.55 16 9 C -0.11 -1.50 4.84 -26.73 -0.13 -1.63 16 10 C -0.09 -1.62 2.48 -88.47 -0.22 -1.84 16 11 H 0.08 1.33 8.05 -51.93 -0.42 0.91 16 12 C 0.14 2.94 6.04 -2.51 -0.02 2.92 16 13 N -0.61 -15.91 3.40 -170.04 -0.58 -16.48 16 14 C -0.23 -6.50 10.28 -15.06 -0.15 -6.65 16 15 H 0.08 2.12 7.83 -52.49 -0.41 1.71 16 16 C -0.01 -0.28 5.49 -17.43 -0.10 -0.38 16 17 N -0.92 -27.64 10.31 55.49 0.57 -27.07 16 18 Si 0.91 22.66 29.55 -169.99 -5.02 17.63 16 19 H -0.29 -7.00 7.11 56.52 0.40 -6.60 16 20 H -0.29 -7.25 7.11 56.52 0.40 -6.85 16 21 H -0.29 -7.00 7.11 56.52 0.40 -6.60 16 22 C 0.18 4.51 4.97 -3.06 -0.02 4.49 16 23 C -0.11 -2.10 3.02 -88.78 -0.27 -2.36 16 24 H 0.06 1.12 8.04 -51.93 -0.42 0.71 16 25 C -0.13 -2.07 9.16 37.16 0.34 -1.73 16 26 C -0.13 -1.05 5.87 -24.34 -0.14 -1.20 16 27 C 0.05 0.48 7.57 38.18 0.29 0.77 16 28 O -0.38 -4.27 11.36 -35.23 -0.40 -4.67 16 29 C 0.07 0.67 6.71 37.88 0.25 0.92 16 30 H 0.04 0.16 8.14 -51.93 -0.42 -0.26 16 31 H 0.09 0.39 8.14 -51.93 -0.42 -0.04 16 32 H 0.04 0.20 8.14 -51.93 -0.42 -0.22 16 33 H 0.06 0.41 8.14 -51.93 -0.42 -0.01 16 34 H 0.07 0.32 8.13 -51.93 -0.42 -0.10 16 35 H 0.41 4.05 6.56 -40.82 -0.27 3.78 16 36 H 0.07 1.09 7.59 -51.92 -0.39 0.70 16 37 H 0.07 0.93 8.14 -51.93 -0.42 0.51 16 38 H 0.07 0.76 8.14 -51.93 -0.42 0.34 16 39 H 0.07 0.88 8.03 -51.93 -0.42 0.47 16 40 H 0.02 0.36 7.96 -51.93 -0.41 -0.06 16 41 H 0.04 0.78 8.10 -51.93 -0.42 0.36 16 42 H 0.25 7.32 8.31 -40.82 -0.34 6.99 16 43 H 0.05 1.26 7.01 -51.93 -0.36 0.90 16 44 H 0.04 1.07 7.50 -51.93 -0.39 0.68 16 45 H 0.04 0.56 8.14 -51.93 -0.42 0.14 16 46 H 0.05 0.78 8.14 -51.93 -0.42 0.36 16 47 H 0.05 0.88 8.06 -51.93 -0.42 0.46 16 48 H 0.10 0.92 7.28 -51.93 -0.38 0.54 16 49 H 0.09 0.48 8.04 -51.93 -0.42 0.06 16 50 H 0.07 0.78 8.13 -51.93 -0.42 0.35 16 51 H 0.09 0.65 8.14 -51.93 -0.42 0.22 16 52 H 0.08 0.62 7.81 -51.93 -0.41 0.21 16 53 H 0.08 0.74 7.44 -51.93 -0.39 0.36 16 LS Contribution 412.64 15.07 6.22 6.22 Total: -1.00 -33.86 412.64 -8.48 -42.34 By element: Atomic # 1 Polarization: 9.72 SS G_CDS: -9.70 Total: 0.01 kcal Atomic # 6 Polarization: 1.31 SS G_CDS: 0.94 Total: 2.25 kcal Atomic # 7 Polarization: -51.85 SS G_CDS: -0.30 Total: -52.15 kcal Atomic # 8 Polarization: -15.70 SS G_CDS: -0.61 Total: -16.30 kcal Atomic # 14 Polarization: 22.66 SS G_CDS: -5.02 Total: 17.63 kcal Total LS contribution 6.22 Total: 6.22 kcal Total: -33.86 -8.48 -42.34 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. ZINC001103276446.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 -71.024 kcal (2) G-P(sol) polarization free energy of solvation -33.863 kcal (3) G-ENP(sol) elect.-nuc.-pol. free energy of system -104.886 kcal (4) G-CDS(sol) cavity-dispersion-solvent structure free energy -8.475 kcal (5) G-P-CDS(sol) = G-P(sol) + G-CDS(sol) = (2) + (4) -42.338 kcal (6) G-S(sol) free energy of system = (1) + (5) -113.362 kcal 1SCF WAS SPECIFIED, SO GEOMETRY OPTIMIZATION WAS NOT USED Total computer time = 1.14 seconds