Wall clock time and date at job start Tue Mar 31 2020 20:48:39 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.17405 * 1 3 3 C 1.47195 * 179.97438 * 2 1 4 4 O 1.42907 * 109.47201 * 156.04315 * 3 2 1 5 5 C 1.42895 * 113.99813 * 179.97438 * 4 3 2 6 6 C 1.53003 * 109.46772 * 179.97438 * 5 4 3 7 7 C 1.50696 * 109.47202 * 180.02562 * 6 5 4 8 8 O 1.21278 * 119.99961 * 359.97438 * 7 6 5 9 9 N 1.34779 * 120.00221 * 180.02562 * 7 6 5 10 10 C 1.46930 * 120.62586 * 5.44909 * 9 7 6 11 11 C 1.53200 * 108.77257 * 126.37117 * 10 9 7 12 12 C 1.53039 * 109.31253 * 54.63210 * 11 10 9 13 13 C 1.53036 * 109.54038 * 298.63320 * 12 11 10 14 14 H 1.08998 * 109.49930 * 181.31565 * 13 12 11 15 15 C 1.52994 * 109.49991 * 301.41716 * 13 12 11 16 16 C 1.46923 * 120.63330 * 185.42623 * 9 7 6 17 17 H 1.09004 * 109.58530 * 113.70009 * 16 9 7 18 18 C 1.52994 * 109.58468 * 353.41846 * 16 9 7 19 19 C 1.53006 * 109.47034 * 180.02562 * 18 16 9 20 20 C 1.53002 * 109.47285 * 175.07095 * 19 18 16 21 21 N 1.46495 * 109.47389 * 180.02562 * 20 19 18 22 22 C 1.36936 * 120.00029 * 180.02562 * 21 20 19 23 23 H 1.07995 * 119.99847 * 359.97438 * 22 21 20 24 24 C 1.38816 * 119.99445 * 179.97438 * 22 21 20 25 25 N 1.31610 * 120.00451 * 359.97438 * 24 22 21 26 26 Si 1.86806 * 119.99370 * 179.97438 * 24 22 21 27 27 H 1.48503 * 109.47242 * 269.99697 * 26 24 22 28 28 H 1.48498 * 109.47491 * 30.00011 * 26 24 22 29 29 H 1.48497 * 109.47021 * 150.00424 * 26 24 22 30 30 H 4.38455 * 5.44472 * 180.02562 * 4 1 2 31 31 H 1.09001 * 109.47339 * 36.04465 * 3 2 1 32 32 H 1.08997 * 109.47422 * 276.04355 * 3 2 1 33 33 H 1.09004 * 109.47212 * 60.00343 * 5 4 3 34 34 H 1.08993 * 109.47460 * 299.99548 * 5 4 3 35 35 H 1.09001 * 109.46978 * 59.99781 * 6 5 4 36 36 H 1.09002 * 109.47306 * 300.00457 * 6 5 4 37 37 H 1.08995 * 109.58408 * 246.15304 * 10 9 7 38 38 H 1.08996 * 109.59200 * 6.58356 * 10 9 7 39 39 H 1.09001 * 109.49448 * 174.58336 * 11 10 9 40 40 H 1.08999 * 109.49843 * 294.67886 * 11 10 9 41 41 H 1.08997 * 109.45893 * 58.65806 * 12 11 10 42 42 H 1.09003 * 109.45674 * 178.60979 * 12 11 10 43 43 H 1.09003 * 109.47755 * 179.83555 * 15 13 12 44 44 H 1.09001 * 109.47039 * 299.83937 * 15 13 12 45 45 H 1.09009 * 109.47284 * 59.83516 * 15 13 12 46 46 H 1.08997 * 109.47527 * 60.05149 * 18 16 9 47 47 H 1.08997 * 109.47304 * 300.04529 * 18 16 9 48 48 H 1.08999 * 109.46987 * 55.07295 * 19 18 16 49 49 H 1.09000 * 109.47417 * 295.07348 * 19 18 16 50 50 H 1.08999 * 109.46944 * 60.00175 * 20 19 18 51 51 H 1.09000 * 109.47100 * 300.00448 * 20 19 18 52 52 H 0.97001 * 120.00115 * 359.97438 * 21 20 19 53 53 H 0.96999 * 120.00476 * 0.02562 * 25 24 22 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.1740 0.0000 0.0000 3 6 2.6460 0.0007 0.0000 4 8 3.1218 1.2321 -0.5471 5 6 4.5463 1.3361 -0.5923 6 6 4.9402 2.6869 -1.1932 7 6 6.4424 2.7963 -1.2416 8 8 7.1302 1.8840 -0.8350 9 7 7.0215 3.9085 -1.7357 10 6 6.2066 5.0691 -2.1202 11 6 6.7451 6.3098 -1.4006 12 6 8.2377 6.4620 -1.7023 13 6 8.9926 5.2365 -1.1825 14 1 10.0584 5.3562 -1.3769 15 6 8.7602 5.0942 0.3230 16 6 8.4796 3.9837 -1.8996 17 1 8.7257 4.0405 -2.9600 18 6 9.1291 2.7406 -1.2883 19 6 10.6474 2.8218 -1.4599 20 6 11.3056 1.6418 -0.7421 21 7 12.7593 1.7199 -0.9058 22 6 13.5637 0.7636 -0.3458 23 1 13.1265 -0.0478 0.2169 24 6 14.9411 0.8374 -0.5015 25 7 15.4739 1.8261 -1.1877 26 14 16.0386 -0.4667 0.2631 27 1 16.4208 -0.0492 1.6360 28 1 15.3057 -1.7568 0.3240 29 1 17.2617 -0.6307 -0.5629 30 1 -1.0502 -0.0003 0.0003 31 1 3.0097 -0.8301 -0.6047 32 1 3.0094 -0.1074 1.0219 33 1 4.9498 0.5330 -1.2091 34 1 4.9487 1.2560 0.4174 35 1 4.5366 3.4902 -0.5769 36 1 4.5377 2.7666 -2.2031 37 1 6.2672 5.2170 -3.1984 38 1 5.1697 4.8998 -1.8302 39 1 6.2112 7.1934 -1.7502 40 1 6.6007 6.1984 -0.3259 41 1 8.3834 6.5486 -2.7790 42 1 8.6178 7.3579 -1.2115 43 1 9.2957 4.2195 0.6923 44 1 9.1252 5.9858 0.8328 45 1 7.6940 4.9759 0.5165 46 1 8.7541 1.8497 -1.7920 47 1 8.8854 2.6879 -0.2273 48 1 11.0118 3.7560 -1.0327 49 1 10.8956 2.7862 -2.5206 50 1 10.9412 0.7076 -1.1693 51 1 11.0574 1.6775 0.3187 52 1 13.1520 2.4486 -1.4116 53 1 14.9041 2.5035 -1.5844 RHF calculation, no. of doubly occupied orbitals= 65 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 ZINC001489063663.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 20:48:39 Heat of formation + Delta-G solvation = -12.787638 kcal Electronic energy + Delta-G solvation = -29613.686547 eV Core-core repulsion = 25680.625826 eV Total energy + Delta-G solvation = -3933.060721 eV No. of doubly occupied orbitals = 65 Molecular weight (most abundant/longest-lived isotopes) = 336.220 amu Computer time = 1.21 seconds Orbital eigenvalues (eV) -41.05404 -38.69658 -37.85051 -35.81999 -34.31206 -33.77512 -33.36824 -31.78874 -31.36890 -29.39010 -27.63559 -27.03331 -26.03674 -24.56858 -23.80233 -23.03230 -21.67368 -20.84263 -20.04273 -18.91753 -18.49542 -17.49876 -17.08494 -16.83012 -16.37287 -16.03172 -15.81176 -15.44947 -15.25776 -14.65148 -14.55698 -14.42361 -13.95316 -13.75311 -13.67052 -13.43847 -13.26568 -13.19820 -12.94557 -12.73032 -12.19040 -12.15596 -11.96581 -11.77793 -11.69154 -11.45219 -11.25733 -11.04858 -10.94217 -10.86320 -10.85270 -10.69063 -10.61296 -10.40588 -10.33638 -10.06957 -9.95140 -9.93540 -9.12287 -9.00346 -8.76918 -8.66256 -6.88198 -6.01735 -2.99369 2.03885 2.38774 2.44612 2.71626 3.48817 3.56887 3.62487 3.76383 4.06796 4.32154 4.39821 4.56646 4.60161 4.67995 4.74798 4.75427 4.81161 4.86269 4.96799 4.98272 5.00557 5.16888 5.18461 5.35063 5.38921 5.42740 5.46054 5.49620 5.51496 5.56158 5.58950 5.61315 5.68515 5.73112 5.81108 5.92583 5.95426 6.00193 6.16443 6.26401 6.27858 6.36719 6.46374 6.49744 6.69189 6.80163 7.00709 7.02203 7.14091 7.24128 7.91949 8.03438 8.85402 9.34912 10.21231 10.37458 11.47520 Molecular weight = 336.22amu Principal moments of inertia in cm(-1) A = 0.010241 B = 0.002375 C = 0.001977 Principal moments of inertia in units of 10**(-40)*gram-cm**2 A = 2733.359883 B =11786.897882 C =14157.739485 Net atomic charges, atomic populations, and dipole contributions calculated with CM2 Atom NO. Type Charge No. of electrons 1 C -0.177 4.177 2 C -0.191 4.191 3 C 0.173 3.827 4 O -0.388 6.388 5 C 0.077 3.923 6 C -0.140 4.140 7 C 0.516 3.484 8 O -0.522 6.522 9 N -0.606 5.606 10 C 0.106 3.894 11 C -0.138 4.138 12 C -0.116 4.116 13 C -0.106 4.106 14 H 0.087 0.913 15 C -0.146 4.146 16 C 0.153 3.847 17 H 0.071 0.929 18 C -0.102 4.102 19 C -0.135 4.135 20 C 0.163 3.837 21 N -0.746 5.746 22 C -0.327 4.327 23 H 0.064 0.936 24 C -0.050 4.050 25 N -0.948 5.948 26 Si 1.083 2.917 27 H -0.296 1.296 28 H -0.295 1.295 29 H -0.288 1.288 30 H 0.222 0.778 31 H 0.077 0.923 32 H 0.077 0.923 33 H 0.063 0.937 34 H 0.063 0.937 35 H 0.099 0.901 36 H 0.100 0.900 37 H 0.069 0.931 38 H 0.085 0.915 39 H 0.073 0.927 40 H 0.071 0.929 41 H 0.063 0.937 42 H 0.065 0.935 43 H 0.065 0.935 44 H 0.051 0.949 45 H 0.052 0.948 46 H 0.082 0.918 47 H 0.080 0.920 48 H 0.059 0.941 49 H 0.054 0.946 50 H 0.018 0.982 51 H 0.017 0.983 52 H 0.365 0.635 53 H 0.254 0.746 For the charges calculated by CM2: Dipole moment (debyes) X Y Z Total from point charges -30.356 7.960 -1.202 31.406 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.194 4.194 2 C -0.191 4.191 3 C 0.098 3.902 4 O -0.308 6.308 5 C 0.000 4.000 6 C -0.180 4.180 7 C 0.303 3.697 8 O -0.398 6.398 9 N -0.340 5.340 10 C -0.017 4.017 11 C -0.176 4.176 12 C -0.153 4.153 13 C -0.125 4.125 14 H 0.105 0.895 15 C -0.203 4.203 16 C 0.051 3.949 17 H 0.089 0.911 18 C -0.140 4.140 19 C -0.173 4.173 20 C 0.036 3.964 21 N -0.392 5.392 22 C -0.458 4.458 23 H 0.082 0.918 24 C -0.247 4.247 25 N -0.594 5.594 26 Si 0.914 3.086 27 H -0.224 1.224 28 H -0.222 1.222 29 H -0.214 1.214 30 H 0.239 0.761 31 H 0.095 0.905 32 H 0.095 0.905 33 H 0.082 0.918 34 H 0.081 0.919 35 H 0.117 0.883 36 H 0.118 0.882 37 H 0.088 0.912 38 H 0.103 0.897 39 H 0.092 0.908 40 H 0.090 0.910 41 H 0.082 0.918 42 H 0.084 0.916 43 H 0.084 0.916 44 H 0.070 0.930 45 H 0.071 0.929 46 H 0.100 0.900 47 H 0.098 0.902 48 H 0.077 0.923 49 H 0.073 0.927 50 H 0.035 0.965 51 H 0.035 0.965 52 H 0.198 0.802 53 H 0.065 0.935 Dipole moment (debyes) X Y Z Total from point charges -28.625 6.913 -0.723 29.457 hybrid contribution -1.608 0.723 -0.333 1.794 sum -30.233 7.635 -1.056 31.200 Atomic orbital electron populations 1.24987 0.95898 0.98149 1.00410 1.22582 0.93916 1.01841 1.00810 1.19060 0.84570 0.87483 0.99109 1.87844 1.21838 1.36957 1.84131 1.22184 0.79638 0.96862 1.01275 1.21544 0.93344 0.99024 1.04117 1.20298 0.89673 0.82164 0.77516 1.90843 1.63578 1.39253 1.46165 1.48100 1.07518 1.14431 1.63907 1.21775 0.94917 0.85456 0.99587 1.21966 0.96815 0.97290 1.01500 1.21563 0.94623 0.98068 1.01080 1.21576 1.00625 0.94064 0.96259 0.89455 1.21876 1.01576 1.02064 0.94809 1.21417 0.77844 0.96464 0.99193 0.91144 1.22088 0.93210 0.95904 1.02793 1.21951 0.95140 0.99657 1.00565 1.20502 0.84586 0.94839 0.96517 1.42415 1.00753 1.30899 1.65107 1.20107 0.90977 1.04896 1.29806 0.91820 1.25770 0.94648 1.03068 1.01178 1.70835 1.37108 1.14097 1.37381 0.83186 0.75594 0.74201 0.75571 1.22391 1.22197 1.21424 0.76145 0.90517 0.90528 0.91833 0.91851 0.88264 0.88203 0.91226 0.89705 0.90844 0.90997 0.91821 0.91585 0.91643 0.93013 0.92865 0.89985 0.90178 0.92252 0.92682 0.96457 0.96523 0.80209 0.93469 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.18 -1.30 19.25 29.55 0.57 -0.73 16 2 C -0.19 -1.71 13.73 -37.12 -0.51 -2.22 16 3 C 0.17 1.40 6.41 34.37 0.22 1.62 16 4 O -0.39 -3.74 10.30 -35.23 -0.36 -4.11 16 5 C 0.08 0.70 5.23 37.16 0.19 0.89 16 6 C -0.14 -1.13 4.48 -27.88 -0.12 -1.25 16 7 C 0.52 5.66 7.61 -10.99 -0.08 5.58 16 8 O -0.52 -7.64 11.39 5.56 0.06 -7.58 16 9 N -0.61 -5.92 2.41 -165.22 -0.40 -6.32 16 10 C 0.11 0.69 6.33 -3.70 -0.02 0.67 16 11 C -0.14 -0.89 5.74 -26.61 -0.15 -1.05 16 12 C -0.12 -0.93 5.43 -26.69 -0.14 -1.08 16 13 C -0.11 -1.13 2.31 -90.50 -0.21 -1.34 16 14 H 0.09 1.03 6.66 -51.93 -0.35 0.68 16 15 C -0.15 -1.65 8.20 37.15 0.30 -1.35 16 16 C 0.15 1.78 2.72 -67.60 -0.18 1.60 16 17 H 0.07 0.81 8.14 -51.93 -0.42 0.39 16 18 C -0.10 -1.54 2.65 -26.73 -0.07 -1.61 16 19 C -0.13 -2.36 5.13 -26.73 -0.14 -2.50 16 20 C 0.16 3.68 6.26 -4.04 -0.03 3.65 16 21 N -0.75 -20.35 5.60 -59.91 -0.34 -20.69 16 22 C -0.33 -9.86 9.99 -15.06 -0.15 -10.01 16 23 H 0.06 1.90 7.83 -52.49 -0.41 1.49 16 24 C -0.05 -1.52 5.50 -17.43 -0.10 -1.61 16 25 N -0.95 -30.49 13.70 55.49 0.76 -29.73 16 26 Si 1.08 27.15 29.92 -169.99 -5.09 22.06 16 27 H -0.30 -7.18 7.11 56.52 0.40 -6.78 16 28 H -0.30 -7.00 7.11 56.52 0.40 -6.60 16 29 H -0.29 -7.22 7.11 56.52 0.40 -6.82 16 30 H 0.22 0.76 7.80 -54.08 -0.42 0.34 16 31 H 0.08 0.52 8.14 -51.93 -0.42 0.09 16 32 H 0.08 0.51 8.14 -51.93 -0.42 0.09 16 33 H 0.06 0.60 8.10 -51.93 -0.42 0.18 16 34 H 0.06 0.59 8.10 -51.93 -0.42 0.17 16 35 H 0.10 0.66 7.46 -51.93 -0.39 0.27 16 36 H 0.10 0.68 8.12 -51.93 -0.42 0.25 16 37 H 0.07 0.41 8.14 -51.93 -0.42 -0.02 16 38 H 0.08 0.41 5.96 -51.93 -0.31 0.10 16 39 H 0.07 0.36 8.14 -51.93 -0.42 -0.06 16 40 H 0.07 0.52 6.33 -51.93 -0.33 0.19 16 41 H 0.06 0.49 8.14 -51.93 -0.42 0.06 16 42 H 0.07 0.51 8.14 -51.93 -0.42 0.09 16 43 H 0.06 0.89 6.34 -51.93 -0.33 0.56 16 44 H 0.05 0.54 8.14 -51.93 -0.42 0.12 16 45 H 0.05 0.57 5.97 -51.92 -0.31 0.26 16 46 H 0.08 1.38 6.45 -51.93 -0.34 1.04 16 47 H 0.08 1.31 5.21 -51.93 -0.27 1.04 16 48 H 0.06 1.02 6.33 -51.93 -0.33 0.69 16 49 H 0.05 0.97 8.14 -51.93 -0.42 0.54 16 50 H 0.02 0.41 8.14 -51.93 -0.42 -0.01 16 51 H 0.02 0.39 8.14 -51.93 -0.42 -0.03 16 52 H 0.37 10.02 6.50 -40.82 -0.27 9.76 16 53 H 0.25 8.05 6.77 -40.82 -0.28 7.77 16 LS Contribution 411.10 15.07 6.20 6.20 Total: -1.00 -36.21 411.10 -8.82 -45.02 By element: Atomic # 1 Polarization: 14.91 SS G_CDS: -9.03 Total: 5.88 kcal Atomic # 6 Polarization: -10.11 SS G_CDS: -0.62 Total: -10.73 kcal Atomic # 7 Polarization: -56.77 SS G_CDS: 0.03 Total: -56.74 kcal Atomic # 8 Polarization: -11.38 SS G_CDS: -0.30 Total: -11.68 kcal Atomic # 14 Polarization: 27.15 SS G_CDS: -5.09 Total: 22.06 kcal Total LS contribution 6.20 Total: 6.20 kcal Total: -36.21 -8.82 -45.02 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. ZINC001489063663.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 32.235 kcal (2) G-P(sol) polarization free energy of solvation -36.207 kcal (3) G-ENP(sol) elect.-nuc.-pol. free energy of system -3.971 kcal (4) G-CDS(sol) cavity-dispersion-solvent structure free energy -8.816 kcal (5) G-P-CDS(sol) = G-P(sol) + G-CDS(sol) = (2) + (4) -45.023 kcal (6) G-S(sol) free energy of system = (1) + (5) -12.788 kcal 1SCF WAS SPECIFIED, SO GEOMETRY OPTIMIZATION WAS NOT USED Total computer time = 1.21 seconds