Wall clock time and date at job start Wed Apr 1 2020 00:33:29 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.52997 * 1 3 3 H 1.09001 * 110.91470 * 2 1 4 4 C 1.54910 * 111.00454 * 236.10274 * 2 1 3 5 5 N 1.47426 * 104.83373 * 204.91811 * 4 2 1 6 6 C 1.36936 * 125.64602 * 204.16653 * 5 4 2 7 7 H 1.07998 * 119.99933 * 179.97438 * 6 5 4 8 8 C 1.38810 * 120.00201 * 359.95099 * 6 5 4 9 9 N 1.31616 * 120.00069 * 0.02562 * 8 6 5 10 10 Si 1.86807 * 119.99974 * 180.02562 * 8 6 5 11 11 H 1.48502 * 109.47126 * 0.02562 * 10 8 6 12 12 H 1.48499 * 109.47030 * 120.00067 * 10 8 6 13 13 H 1.48501 * 109.47106 * 240.00199 * 10 8 6 14 14 C 1.47022 * 108.69967 * 24.12539 * 5 4 2 15 15 C 1.54314 * 107.27471 * 358.94486 * 14 5 4 16 16 H 1.09000 * 110.72176 * 96.34615 * 15 14 5 17 17 C 1.53000 * 110.72156 * 219.44242 * 15 14 5 18 18 C 1.52997 * 109.47273 * 293.54759 * 17 15 14 19 19 N 1.46499 * 109.47555 * 179.97438 * 18 17 15 20 20 C 1.34777 * 120.00199 * 180.02562 * 19 18 17 21 21 O 1.21283 * 119.99750 * 359.97438 * 20 19 18 22 22 C 1.50702 * 119.99928 * 180.27479 * 20 19 18 23 23 O 1.42899 * 109.46876 * 179.72354 * 22 20 19 24 24 C 1.42901 * 113.99607 * 179.97438 * 23 22 20 25 25 C 1.53006 * 109.46807 * 179.97438 * 24 23 22 26 26 H 1.08996 * 110.43192 * 303.90616 * 25 24 23 27 27 C 1.54812 * 110.48642 * 181.22087 * 25 24 23 28 28 C 1.54993 * 102.09134 * 220.44163 * 27 25 24 29 29 O 1.44383 * 103.48863 * 37.02437 * 28 27 25 30 30 C 1.43743 * 106.97614 * 319.93942 * 29 28 27 31 31 H 1.08999 * 109.47372 * 296.09698 * 1 2 3 32 32 H 1.09001 * 109.47435 * 56.09579 * 1 2 3 33 33 H 1.08997 * 109.47297 * 176.09118 * 1 2 3 34 34 H 1.08999 * 110.36573 * 86.07594 * 4 2 1 35 35 H 1.09005 * 110.45929 * 323.80289 * 4 2 1 36 36 H 0.97003 * 120.00344 * 179.97438 * 9 8 6 37 37 H 1.08999 * 109.87818 * 239.51944 * 14 5 4 38 38 H 1.09004 * 110.00736 * 118.43464 * 14 5 4 39 39 H 1.08998 * 109.47307 * 173.54235 * 17 15 14 40 40 H 1.09005 * 109.46653 * 53.54169 * 17 15 14 41 41 H 1.08998 * 109.47382 * 59.99521 * 18 17 15 42 42 H 1.09002 * 109.47216 * 299.99958 * 18 17 15 43 43 H 0.97000 * 119.99484 * 359.97438 * 19 18 17 44 44 H 1.09004 * 109.46965 * 299.72211 * 22 20 19 45 45 H 1.08994 * 109.46990 * 59.72065 * 22 20 19 46 46 H 1.09003 * 109.47075 * 299.99807 * 24 23 22 47 47 H 1.09001 * 109.47248 * 60.00362 * 24 23 22 48 48 H 1.09006 * 110.90915 * 338.64333 * 27 25 24 49 49 H 1.09003 * 110.90596 * 102.23995 * 27 25 24 50 50 H 1.08995 * 110.62115 * 278.52072 * 28 27 25 51 51 H 1.09003 * 110.61562 * 155.51971 * 28 27 25 52 52 H 1.08998 * 109.87597 * 267.09157 * 30 29 28 53 53 H 1.09000 * 110.00265 * 146.01280 * 30 29 28 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.9191 1.0182 0.0000 4 6 2.0852 -0.8065 -1.2004 5 7 3.4271 -1.2431 -0.7734 6 6 4.4644 -1.5652 -1.6073 7 1 5.4130 -1.8749 -1.1944 8 6 4.2964 -1.4944 -2.9834 9 7 3.1403 -1.1170 -3.4867 10 14 5.7113 -1.9346 -4.1210 11 1 6.8938 -2.3249 -3.3120 12 1 5.3116 -3.0691 -4.9918 13 1 6.0529 -0.7603 -4.9634 14 6 3.4673 -1.2742 0.6959 15 6 2.0861 -0.8081 1.2021 16 1 1.4478 -1.6634 1.4238 17 6 2.2349 0.0874 2.4338 18 6 2.7848 -0.7365 3.5998 19 7 2.9279 0.1210 4.7790 20 6 3.3989 -0.3926 5.9327 21 8 3.7038 -1.5648 5.9947 22 6 3.5402 0.4881 7.1474 23 8 4.0495 -0.2825 8.2377 24 6 4.2225 0.4623 9.4449 25 6 4.7687 -0.4583 10.5382 26 1 4.1114 -1.3168 10.6763 27 6 4.9364 0.3185 11.8668 28 6 6.2826 -0.2354 12.3990 29 8 7.0736 -0.4046 11.2030 30 6 6.1960 -0.9224 10.1892 31 1 -0.3634 0.4520 0.9229 32 1 -0.3634 0.5732 -0.8529 33 1 -0.3634 -1.0252 -0.0701 34 1 2.1559 -0.1718 -2.0837 35 1 1.4524 -1.6703 -1.4045 36 1 3.0230 -1.0671 -4.4483 37 1 4.2426 -0.5999 1.0598 38 1 3.6661 -2.2885 1.0421 39 1 1.2618 0.4967 2.7051 40 1 2.9224 0.9029 2.2091 41 1 3.7579 -1.1458 3.3285 42 1 2.0975 -1.5521 3.8245 43 1 2.6843 1.0586 4.7292 44 1 4.2283 1.3041 6.9265 45 1 2.5659 0.8968 7.4151 46 1 4.9250 1.2775 9.2714 47 1 3.2625 0.8708 9.7606 48 1 5.0019 1.3912 11.6846 49 1 4.1227 0.0914 12.5556 50 1 6.1332 -1.1924 12.8989 51 1 6.7532 0.4811 13.0722 52 1 6.2402 -2.0114 10.1818 53 1 6.4839 -0.5333 9.2126 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 ZINC001103487734.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:33:29 Heat of formation + Delta-G solvation = -116.053576 kcal Electronic energy + Delta-G solvation = -29157.629283 eV Core-core repulsion = 25029.634564 eV Total energy + Delta-G solvation = -4127.994719 eV No. of doubly occupied orbitals = 66 Molecular weight (most abundant/longest-lived isotopes) = 340.220 amu Computer time = 1.20 seconds Orbital eigenvalues (eV) -40.15293 -39.71711 -38.25265 -37.66877 -35.23093 -33.36384 -33.24802 -31.36425 -30.67801 -30.08898 -28.51795 -26.60128 -25.80259 -25.43526 -23.98733 -22.94021 -22.03984 -21.15305 -20.70860 -19.67199 -19.22018 -17.58213 -17.24645 -17.11381 -16.79031 -16.37878 -16.29034 -15.68695 -15.43526 -15.33987 -14.89862 -14.50288 -14.33414 -13.84607 -13.70420 -13.49412 -13.28494 -12.98666 -12.85905 -12.72750 -12.32239 -12.19334 -12.05405 -11.99735 -11.73617 -11.65654 -11.50672 -11.47625 -11.17480 -11.11913 -11.00870 -10.83438 -10.73936 -10.16555 -10.11482 -10.01980 -9.99331 -9.85784 -9.78247 -9.69924 -9.30627 -9.14459 -8.62309 -6.78718 -5.62174 -3.03394 2.47354 2.53330 3.41526 3.47389 3.49761 3.51428 3.56346 4.03059 4.34849 4.36856 4.51532 4.52430 4.66807 4.72142 4.75049 4.80452 4.84508 4.92107 4.98301 5.10248 5.13667 5.23728 5.28470 5.32888 5.38058 5.47478 5.53829 5.55061 5.57170 5.66597 5.67437 5.77985 5.93721 5.97607 6.14460 6.24342 6.26404 6.32051 6.35102 6.44399 6.56396 6.61916 6.63380 6.71646 6.84916 7.03840 7.34720 7.57966 7.67808 7.89757 8.08934 8.39997 9.26071 9.87020 10.52378 11.41569 Molecular weight = 340.22amu Principal moments of inertia in cm(-1) A = 0.017758 B = 0.001624 C = 0.001509 Principal moments of inertia in units of 10**(-40)*gram-cm**2 A = 1576.356434 B =17234.632351 C =18550.345550 Net atomic charges, atomic populations, and dipole contributions calculated with CM2 Atom NO. Type Charge No. of electrons 1 C -0.133 4.133 2 C -0.115 4.115 3 H 0.068 0.932 4 C 0.158 3.842 5 N -0.606 5.606 6 C -0.220 4.220 7 H 0.077 0.923 8 C -0.013 4.013 9 N -0.923 5.923 10 Si 0.913 3.087 11 H -0.280 1.280 12 H -0.291 1.291 13 H -0.291 1.291 14 C 0.139 3.861 15 C -0.093 4.093 16 H 0.073 0.927 17 C -0.115 4.115 18 C 0.119 3.881 19 N -0.728 5.728 20 C 0.509 3.491 21 O -0.518 6.518 22 C 0.045 3.955 23 O -0.362 6.362 24 C 0.069 3.931 25 C -0.121 4.121 26 H 0.094 0.906 27 C -0.148 4.148 28 C 0.046 3.954 29 O -0.380 6.380 30 C 0.053 3.947 31 H 0.042 0.958 32 H 0.052 0.948 33 H 0.051 0.949 34 H 0.104 0.896 35 H 0.012 0.988 36 H 0.250 0.750 37 H 0.022 0.978 38 H 0.020 0.980 39 H 0.065 0.935 40 H 0.069 0.931 41 H 0.074 0.926 42 H 0.068 0.932 43 H 0.399 0.601 44 H 0.080 0.920 45 H 0.079 0.921 46 H 0.056 0.944 47 H 0.055 0.945 48 H 0.080 0.920 49 H 0.084 0.916 50 H 0.062 0.938 51 H 0.094 0.906 52 H 0.063 0.937 53 H 0.091 0.909 For the charges calculated by CM2: Dipole moment (debyes) X Y Z Total from point charges -3.373 8.627 29.344 30.772 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.191 4.191 2 C -0.134 4.134 3 H 0.086 0.914 4 C 0.033 3.967 5 N -0.331 5.331 6 C -0.348 4.348 7 H 0.095 0.905 8 C -0.208 4.208 9 N -0.572 5.572 10 Si 0.744 3.256 11 H -0.205 1.205 12 H -0.218 1.218 13 H -0.218 1.218 14 C 0.013 3.987 15 C -0.112 4.112 16 H 0.091 0.909 17 C -0.154 4.154 18 C -0.004 4.004 19 N -0.381 5.381 20 C 0.294 3.706 21 O -0.393 6.393 22 C -0.035 4.035 23 O -0.280 6.280 24 C -0.008 4.008 25 C -0.140 4.140 26 H 0.113 0.887 27 C -0.186 4.186 28 C -0.031 4.031 29 O -0.299 6.299 30 C -0.023 4.023 31 H 0.061 0.939 32 H 0.071 0.929 33 H 0.070 0.930 34 H 0.122 0.878 35 H 0.029 0.971 36 H 0.061 0.939 37 H 0.040 0.960 38 H 0.038 0.962 39 H 0.084 0.916 40 H 0.088 0.912 41 H 0.092 0.908 42 H 0.086 0.914 43 H 0.233 0.767 44 H 0.098 0.902 45 H 0.097 0.903 46 H 0.074 0.926 47 H 0.073 0.927 48 H 0.099 0.901 49 H 0.103 0.897 50 H 0.080 0.920 51 H 0.112 0.888 52 H 0.081 0.919 53 H 0.109 0.891 Dipole moment (debyes) X Y Z Total from point charges -2.801 7.560 29.125 30.221 hybrid contribution 0.183 0.405 -0.617 0.761 sum -2.618 7.965 28.508 29.715 Atomic orbital electron populations 1.21634 0.95778 1.01345 1.00329 1.22632 0.94359 0.99323 0.97106 0.91388 1.21970 0.84112 0.95851 0.94718 1.44506 1.10494 1.78937 0.99164 1.19147 0.91064 1.37356 0.87251 0.90488 1.24779 0.97273 1.02161 0.96609 1.69813 1.22350 1.54144 1.10863 0.86012 0.81206 0.77582 0.80818 1.20536 1.21790 1.21792 1.21368 0.93696 0.97265 0.86364 1.22144 0.98166 0.97637 0.93286 0.90886 1.21781 1.00932 0.98684 0.93956 1.21525 1.02198 0.93290 0.83382 1.45960 1.66972 1.14033 1.11173 1.20191 0.77844 0.88604 0.83968 1.90758 1.47776 1.16749 1.83968 1.22205 1.00162 0.93930 0.87184 1.88010 1.85626 1.35597 1.18773 1.22597 0.99781 0.93971 0.84419 1.22544 0.98218 0.99397 0.93821 0.88744 1.23156 0.98753 1.00925 0.95774 1.23545 0.90444 1.00359 0.88731 1.88897 1.34085 1.85904 1.21058 1.23218 0.87695 0.98897 0.92537 0.93866 0.92903 0.93009 0.87799 0.97051 0.93869 0.96007 0.96249 0.91615 0.91194 0.90817 0.91413 0.76686 0.90172 0.90272 0.92556 0.92688 0.90130 0.89714 0.92012 0.88828 0.91871 0.89086 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.13 -2.00 9.31 37.16 0.35 -1.65 16 2 C -0.11 -2.07 3.09 -88.47 -0.27 -2.35 16 3 H 0.07 1.23 8.05 -51.93 -0.42 0.81 16 4 C 0.16 3.82 5.36 -2.53 -0.01 3.81 16 5 N -0.61 -15.68 3.39 -170.08 -0.58 -16.25 16 6 C -0.22 -6.23 10.27 -15.06 -0.15 -6.39 16 7 H 0.08 2.13 7.83 -52.49 -0.41 1.72 16 8 C -0.01 -0.37 5.49 -17.43 -0.10 -0.47 16 9 N -0.92 -27.09 10.09 55.49 0.56 -26.53 16 10 Si 0.91 22.39 29.55 -169.99 -5.02 17.37 16 11 H -0.28 -6.70 7.11 56.52 0.40 -6.30 16 12 H -0.29 -7.18 7.11 56.52 0.40 -6.78 16 13 H -0.29 -7.13 7.11 56.52 0.40 -6.73 16 14 C 0.14 2.91 5.99 -3.07 -0.02 2.89 16 15 C -0.09 -1.54 2.45 -88.79 -0.22 -1.75 16 16 H 0.07 1.20 8.04 -51.93 -0.42 0.78 16 17 C -0.12 -1.43 4.72 -26.73 -0.13 -1.55 16 18 C 0.12 1.50 5.52 -4.04 -0.02 1.47 16 19 N -0.73 -7.35 5.56 -60.29 -0.33 -7.68 16 20 C 0.51 6.14 7.91 -10.98 -0.09 6.05 16 21 O -0.52 -8.93 16.50 -14.35 -0.24 -9.16 16 22 C 0.05 0.37 6.42 36.01 0.23 0.60 16 23 O -0.36 -3.73 9.62 -55.14 -0.53 -4.26 16 24 C 0.07 0.47 5.03 37.16 0.19 0.66 16 25 C -0.12 -0.83 3.34 -89.09 -0.30 -1.13 16 26 H 0.09 0.66 8.14 -51.93 -0.42 0.24 16 27 C -0.15 -0.71 6.71 -24.72 -0.17 -0.88 16 28 C 0.05 0.27 7.72 38.21 0.29 0.57 16 29 O -0.38 -3.67 11.52 -35.23 -0.41 -4.08 16 30 C 0.05 0.47 7.23 37.73 0.27 0.74 16 31 H 0.04 0.54 8.14 -51.93 -0.42 0.11 16 32 H 0.05 0.78 8.14 -51.93 -0.42 0.36 16 33 H 0.05 0.79 8.07 -51.93 -0.42 0.37 16 34 H 0.10 2.84 5.87 -51.93 -0.30 2.53 16 35 H 0.01 0.30 7.99 -51.93 -0.41 -0.12 16 36 H 0.25 7.12 8.31 -40.82 -0.34 6.78 16 37 H 0.02 0.45 7.76 -51.93 -0.40 0.05 16 38 H 0.02 0.42 8.14 -51.93 -0.42 0.00 16 39 H 0.07 0.69 8.14 -51.93 -0.42 0.27 16 40 H 0.07 0.87 8.06 -51.93 -0.42 0.45 16 41 H 0.07 1.09 7.49 -51.93 -0.39 0.70 16 42 H 0.07 0.90 8.14 -51.93 -0.42 0.48 16 43 H 0.40 2.65 8.47 -40.82 -0.35 2.30 16 44 H 0.08 0.46 8.14 -51.93 -0.42 0.04 16 45 H 0.08 0.42 8.14 -51.93 -0.42 0.00 16 46 H 0.06 0.34 8.14 -51.93 -0.42 -0.09 16 47 H 0.05 0.28 8.14 -51.93 -0.42 -0.15 16 48 H 0.08 0.34 7.87 -51.93 -0.41 -0.07 16 49 H 0.08 0.28 8.14 -51.93 -0.42 -0.15 16 50 H 0.06 0.32 8.14 -51.93 -0.42 -0.10 16 51 H 0.09 0.42 8.14 -51.93 -0.42 0.00 16 52 H 0.06 0.56 8.14 -51.93 -0.42 0.14 16 53 H 0.09 0.88 7.60 -51.93 -0.39 0.48 16 LS Contribution 419.49 15.07 6.32 6.32 Total: -1.00 -35.35 419.49 -10.17 -45.52 By element: Atomic # 1 Polarization: 7.93 SS G_CDS: -9.80 Total: -1.87 kcal Atomic # 6 Polarization: 0.76 SS G_CDS: -0.14 Total: 0.63 kcal Atomic # 7 Polarization: -50.11 SS G_CDS: -0.35 Total: -50.47 kcal Atomic # 8 Polarization: -16.32 SS G_CDS: -1.17 Total: -17.50 kcal Atomic # 14 Polarization: 22.39 SS G_CDS: -5.02 Total: 17.37 kcal Total LS contribution 6.32 Total: 6.32 kcal Total: -35.35 -10.17 -45.52 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. ZINC001103487734.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 -70.538 kcal (2) G-P(sol) polarization free energy of solvation -35.350 kcal (3) G-ENP(sol) elect.-nuc.-pol. free energy of system -105.888 kcal (4) G-CDS(sol) cavity-dispersion-solvent structure free energy -10.166 kcal (5) G-P-CDS(sol) = G-P(sol) + G-CDS(sol) = (2) + (4) -45.516 kcal (6) G-S(sol) free energy of system = (1) + (5) -116.054 kcal 1SCF WAS SPECIFIED, SO GEOMETRY OPTIMIZATION WAS NOT USED Total computer time = 1.20 seconds