Wall clock time and date at job start Tue Mar 31 2020 21:38:19 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.52996 * 1 3 3 H 1.09000 * 109.51217 * 2 1 4 4 C 1.53288 * 109.55341 * 239.87794 * 2 1 3 5 5 N 1.47154 * 108.37808 * 187.70830 * 4 2 1 6 6 C 1.34776 * 120.98798 * 230.61638 * 5 4 2 7 7 O 1.21287 * 119.99780 * 179.97438 * 6 5 4 8 8 C 1.50700 * 120.00056 * 359.97438 * 6 5 4 9 9 C 1.53004 * 109.46743 * 179.97438 * 8 6 5 10 10 C 1.52999 * 109.46845 * 180.02562 * 9 8 6 11 11 C 1.52997 * 117.49881 * 243.75023 * 10 9 8 12 12 C 1.52999 * 117.49899 * 175.08992 * 10 9 8 13 13 C 1.47158 * 118.02035 * 50.59389 * 5 4 2 14 14 C 1.53258 * 108.40619 * 309.51347 * 13 5 4 15 15 C 1.52997 * 109.35339 * 172.09736 * 14 13 5 16 16 C 1.53193 * 109.34492 * 181.13996 * 15 14 13 17 17 N 1.46934 * 108.77292 * 305.41334 * 16 15 14 18 18 C 1.36935 * 120.64434 * 233.66395 * 17 16 15 19 19 H 1.07995 * 119.99817 * 151.84246 * 18 17 16 20 20 C 1.38813 * 119.99897 * 331.83657 * 18 17 16 21 21 N 1.31616 * 120.00152 * 343.07000 * 20 18 17 22 22 Si 1.86796 * 119.99721 * 163.06102 * 20 18 17 23 23 H 1.48503 * 109.47213 * 90.00760 * 22 20 18 24 24 H 1.48499 * 109.47048 * 210.00609 * 22 20 18 25 25 H 1.48498 * 109.47533 * 330.00336 * 22 20 18 26 26 C 1.46923 * 118.70989 * 53.68292 * 17 16 15 27 27 C 1.52984 * 109.38939 * 292.18005 * 14 13 5 28 28 O 1.42891 * 109.32351 * 52.15841 * 14 13 5 29 29 H 1.09003 * 109.47218 * 299.91107 * 1 2 3 30 30 H 1.09007 * 109.46678 * 59.90699 * 1 2 3 31 31 H 1.08996 * 109.47155 * 179.90340 * 1 2 3 32 32 H 1.09001 * 109.68208 * 307.42388 * 4 2 1 33 33 H 1.09000 * 109.68633 * 67.94591 * 4 2 1 34 34 H 1.08990 * 109.47135 * 299.99606 * 8 6 5 35 35 H 1.09000 * 109.47163 * 60.00009 * 8 6 5 36 36 H 1.08994 * 109.47633 * 300.00004 * 9 8 6 37 37 H 1.09000 * 109.47087 * 60.00203 * 9 8 6 38 38 H 1.08998 * 115.55153 * 29.41420 * 10 9 8 39 39 H 1.08999 * 117.50338 * 359.97438 * 11 10 9 40 40 H 1.09001 * 117.50227 * 145.02263 * 11 10 9 41 41 H 1.09000 * 117.49614 * 214.97709 * 12 10 9 42 42 H 1.09001 * 117.50042 * 359.97438 * 12 10 9 43 43 H 1.08997 * 109.67574 * 69.24413 * 13 5 4 44 44 H 1.09001 * 109.67244 * 189.74689 * 13 5 4 45 45 H 1.08996 * 109.49220 * 301.09683 * 15 14 13 46 46 H 1.09000 * 109.49217 * 61.18302 * 15 14 13 47 47 H 1.08997 * 109.58740 * 185.62464 * 16 15 14 48 48 H 1.09008 * 109.70470 * 65.27043 * 16 15 14 49 49 H 0.97001 * 120.00317 * 179.97438 * 21 20 18 50 50 H 1.09002 * 109.58682 * 186.54399 * 26 17 16 51 51 H 1.09001 * 109.70574 * 66.19241 * 26 17 16 52 52 H 1.09000 * 109.49389 * 298.85914 * 27 14 13 53 53 H 1.08993 * 109.49399 * 58.91408 * 27 14 13 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.8940 1.0274 0.0000 4 6 2.0430 -0.7249 -1.2494 5 7 3.5023 -0.8739 -1.1330 6 6 4.3204 -0.4992 -2.1364 7 8 5.5207 -0.6297 -2.0211 8 6 3.7437 0.0817 -3.4017 9 6 4.8783 0.4208 -4.3705 10 6 4.2929 1.0112 -5.6548 11 6 4.5007 0.2266 -6.9518 12 6 5.2762 1.5238 -6.7090 13 6 4.0419 -1.4436 0.1119 14 6 3.4304 -0.6863 1.2957 15 6 3.8205 -1.3841 2.6001 16 6 3.1857 -0.6435 3.7814 17 7 3.5726 0.7725 3.7147 18 6 4.1552 1.3882 4.7902 19 1 4.8282 2.2189 4.6375 20 6 3.8807 0.9446 6.0766 21 7 2.8344 0.1777 6.2990 22 14 4.9963 1.4361 7.4919 23 1 4.5063 2.7032 8.0915 24 1 4.9899 0.3662 8.5217 25 1 6.3796 1.6294 6.9878 26 6 3.3219 1.5149 2.4719 27 6 3.9548 0.7509 1.3046 28 8 2.0074 -0.6756 1.1660 29 1 -0.3634 0.5125 0.8908 30 1 -0.3633 0.5153 -0.8892 31 1 -0.3633 -1.0276 -0.0017 32 1 1.5781 -1.7083 -1.3202 33 1 1.8028 -0.1405 -2.1376 34 1 3.0760 -0.6453 -3.8639 35 1 3.1861 0.9877 -3.1643 36 1 5.5458 1.1480 -3.9082 37 1 5.4363 -0.4849 -4.6079 38 1 3.3532 1.5511 -5.5384 39 1 5.0548 -0.7100 -6.8902 40 1 3.6979 0.2507 -7.6887 41 1 4.9832 2.4008 -7.2863 42 1 6.3402 1.4408 -6.4873 43 1 3.7791 -2.4995 0.1765 44 1 5.1263 -1.3339 0.1259 45 1 3.4637 -2.4139 2.5825 46 1 4.9053 -1.3762 2.7065 47 1 3.5406 -1.0756 4.7169 48 1 2.1003 -0.7286 3.7281 49 1 2.6428 -0.1326 7.1978 50 1 3.7669 2.5074 2.5430 51 1 2.2481 1.6069 2.3090 52 1 5.0386 0.7419 1.4209 53 1 3.6938 1.2391 0.3658 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 ZINC001098577850.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 21:38:19 Heat of formation + Delta-G solvation = -36.923025 kcal Electronic energy + Delta-G solvation = -31108.803157 eV Core-core repulsion = 27174.695848 eV Total energy + Delta-G solvation = -3934.107310 eV No. of doubly occupied orbitals = 65 Molecular weight (most abundant/longest-lived isotopes) = 336.220 amu Computer time = 1.43 seconds Orbital eigenvalues (eV) -40.73025 -39.00904 -37.71248 -36.38136 -34.58262 -33.61208 -32.40805 -31.57432 -30.13398 -29.63628 -27.48615 -26.43708 -25.08091 -23.64129 -22.68423 -22.59936 -21.67076 -20.84540 -20.58753 -19.08935 -18.90697 -17.22463 -16.93562 -16.28795 -16.06003 -15.66776 -15.47969 -14.95919 -14.77659 -14.46090 -14.36221 -14.08764 -13.60445 -13.41667 -13.30023 -13.23419 -12.92601 -12.68297 -12.51208 -12.25665 -12.00365 -11.88147 -11.80513 -11.62871 -11.58962 -11.43777 -11.39807 -11.28936 -10.99286 -10.78022 -10.63330 -10.59291 -10.35490 -10.23258 -10.19144 -10.06373 -9.82714 -9.73346 -9.48856 -9.29666 -8.41406 -8.36948 -6.59932 -5.73948 -3.15984 2.57563 2.85336 3.39939 3.47005 3.47977 3.49834 3.88580 4.20741 4.47223 4.53481 4.59936 4.62868 4.67793 4.84285 4.89565 4.99660 5.07269 5.14253 5.17426 5.23972 5.31863 5.32718 5.40627 5.42648 5.48101 5.54792 5.59132 5.60521 5.65164 5.68277 5.78112 5.80186 5.86067 5.89629 5.98896 6.06308 6.09441 6.17349 6.32089 6.34282 6.37050 6.56239 6.59542 6.70298 6.77900 7.03551 7.32670 7.44289 7.63099 7.77336 7.84127 8.11463 8.50142 9.26813 9.90320 10.48358 11.34365 Molecular weight = 336.22amu Principal moments of inertia in cm(-1) A = 0.018911 B = 0.002717 C = 0.002560 Principal moments of inertia in units of 10**(-40)*gram-cm**2 A = 1480.245157 B =10303.589132 C =10935.403838 Net atomic charges, atomic populations, and dipole contributions calculated with CM2 Atom NO. Type Charge No. of electrons 1 C -0.142 4.142 2 C 0.069 3.931 3 H 0.066 0.934 4 C 0.079 3.921 5 N -0.612 5.612 6 C 0.513 3.487 7 O -0.542 6.542 8 C -0.136 4.136 9 C -0.074 4.074 10 C -0.154 4.154 11 C -0.182 4.182 12 C -0.181 4.181 13 C 0.088 3.912 14 C 0.120 3.880 15 C -0.126 4.126 16 C 0.179 3.821 17 N -0.590 5.590 18 C -0.251 4.251 19 H 0.063 0.937 20 C 0.002 3.998 21 N -0.901 5.901 22 Si 0.900 3.100 23 H -0.288 1.288 24 H -0.290 1.290 25 H -0.280 1.280 26 C 0.153 3.847 27 C -0.152 4.152 28 O -0.353 6.353 29 H 0.070 0.930 30 H 0.067 0.933 31 H 0.063 0.937 32 H 0.075 0.925 33 H 0.089 0.911 34 H 0.092 0.908 35 H 0.093 0.907 36 H 0.074 0.926 37 H 0.071 0.929 38 H 0.105 0.895 39 H 0.094 0.906 40 H 0.092 0.908 41 H 0.093 0.907 42 H 0.095 0.905 43 H 0.070 0.930 44 H 0.097 0.903 45 H 0.058 0.942 46 H 0.070 0.930 47 H 0.086 0.914 48 H 0.028 0.972 49 H 0.253 0.747 50 H 0.049 0.951 51 H 0.020 0.980 52 H 0.067 0.933 53 H 0.049 0.951 For the charges calculated by CM2: Dipole moment (debyes) X Y Z Total from point charges -3.070 -1.747 -23.480 23.745 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.199 4.199 2 C 0.010 3.990 3 H 0.084 0.916 4 C -0.045 4.045 5 N -0.345 5.345 6 C 0.302 3.698 7 O -0.419 6.419 8 C -0.176 4.176 9 C -0.111 4.111 10 C -0.172 4.172 11 C -0.219 4.219 12 C -0.218 4.218 13 C -0.035 4.035 14 C 0.080 3.920 15 C -0.166 4.166 16 C 0.053 3.947 17 N -0.313 5.313 18 C -0.381 4.381 19 H 0.081 0.919 20 C -0.196 4.196 21 N -0.546 5.546 22 Si 0.730 3.270 23 H -0.215 1.215 24 H -0.216 1.216 25 H -0.205 1.205 26 C 0.027 3.973 27 C -0.191 4.191 28 O -0.272 6.272 29 H 0.090 0.910 30 H 0.086 0.914 31 H 0.082 0.918 32 H 0.094 0.906 33 H 0.107 0.893 34 H 0.111 0.889 35 H 0.111 0.889 36 H 0.093 0.907 37 H 0.089 0.911 38 H 0.123 0.877 39 H 0.112 0.888 40 H 0.110 0.890 41 H 0.111 0.889 42 H 0.113 0.887 43 H 0.088 0.912 44 H 0.116 0.884 45 H 0.076 0.924 46 H 0.088 0.912 47 H 0.104 0.896 48 H 0.046 0.954 49 H 0.063 0.937 50 H 0.068 0.932 51 H 0.038 0.962 52 H 0.085 0.915 53 H 0.068 0.932 Dipole moment (debyes) X Y Z Total from point charges -2.592 -1.968 -23.113 23.341 hybrid contribution 0.544 0.316 0.500 0.804 sum -2.048 -1.652 -22.614 22.766 Atomic orbital electron populations 1.21699 0.92600 1.02938 1.02695 1.22384 0.95226 0.96121 0.85279 0.91622 1.22538 0.80188 1.01903 0.99850 1.48184 1.06794 1.62913 1.16631 1.20709 0.86695 0.77688 0.84701 1.90697 1.14883 1.55169 1.81143 1.21672 1.00064 1.02249 0.93569 1.20926 0.96479 1.01411 0.92324 1.22356 1.00688 1.03763 0.90389 1.23064 1.03301 0.95976 0.99588 1.22975 0.98072 0.97801 1.02948 1.22479 1.00817 0.97276 0.82892 1.21739 0.80454 0.95421 0.94389 1.22200 1.01093 0.98628 0.94664 1.20673 0.96227 0.83188 0.94631 1.44022 1.71078 1.06025 1.10154 1.19157 1.19310 1.15080 0.84547 0.91890 1.25035 0.96691 1.00551 0.97364 1.70038 1.25414 1.38769 1.20336 0.86335 0.80712 0.77642 0.82310 1.21468 1.21574 1.20541 1.20752 0.96400 0.93860 0.86315 1.22514 1.01415 0.95763 0.99419 1.88072 1.20974 1.76297 1.41817 0.91045 0.91405 0.91806 0.90632 0.89271 0.88921 0.88882 0.90700 0.91061 0.87721 0.88774 0.88963 0.88861 0.88697 0.91205 0.88428 0.92382 0.91170 0.89603 0.95417 0.93694 0.93222 0.96201 0.91455 0.93207 Number of geometries 1 Number of calculations of the screened coulomb radii 1 The total number of SCF iterations 18. 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.48 9.57 37.16 0.36 -1.12 16 2 C 0.07 0.88 3.01 -26.58 -0.08 0.80 16 3 H 0.07 0.89 6.29 -51.93 -0.33 0.57 16 4 C 0.08 0.81 5.83 -3.54 -0.02 0.79 16 5 N -0.61 -7.97 2.52 -172.30 -0.43 -8.41 16 6 C 0.51 6.80 7.71 -10.99 -0.08 6.71 16 7 O -0.54 -9.22 15.45 5.55 0.09 -9.14 16 8 C -0.14 -1.22 4.03 -27.88 -0.11 -1.33 16 9 C -0.07 -0.69 4.78 -26.73 -0.13 -0.81 16 10 C -0.15 -1.15 5.73 -90.62 -0.52 -1.67 16 11 C -0.18 -1.25 10.18 -26.69 -0.27 -1.53 16 12 C -0.18 -1.30 10.25 -26.69 -0.27 -1.57 16 13 C 0.09 1.33 5.29 -3.55 -0.02 1.31 16 14 C 0.12 2.07 0.68 -90.50 -0.06 2.00 16 15 C -0.13 -2.44 5.10 -26.63 -0.14 -2.58 16 16 C 0.18 4.32 5.07 -3.71 -0.02 4.30 16 17 N -0.59 -15.00 2.99 -172.31 -0.52 -15.52 16 18 C -0.25 -6.99 9.58 -15.06 -0.14 -7.13 16 19 H 0.06 1.73 7.88 -52.49 -0.41 1.32 16 20 C 0.00 0.06 4.92 -17.43 -0.09 -0.03 16 21 N -0.90 -26.36 11.20 55.49 0.62 -25.74 16 22 Si 0.90 21.91 29.59 -169.99 -5.03 16.88 16 23 H -0.29 -6.99 7.11 56.52 0.40 -6.59 16 24 H -0.29 -7.07 7.11 56.52 0.40 -6.67 16 25 H -0.28 -6.72 7.11 56.52 0.40 -6.31 16 26 C 0.15 3.31 6.28 -3.71 -0.02 3.29 16 27 C -0.15 -2.77 4.33 -26.64 -0.12 -2.89 16 28 O -0.35 -5.98 8.65 -35.23 -0.30 -6.28 16 29 H 0.07 0.86 8.14 -51.93 -0.42 0.44 16 30 H 0.07 0.51 8.14 -51.93 -0.42 0.09 16 31 H 0.06 0.64 8.14 -51.93 -0.42 0.22 16 32 H 0.08 0.74 8.14 -51.93 -0.42 0.31 16 33 H 0.09 0.66 6.06 -51.93 -0.31 0.34 16 34 H 0.09 0.66 8.01 -51.93 -0.42 0.25 16 35 H 0.09 0.70 7.66 -51.93 -0.40 0.30 16 36 H 0.07 0.80 8.10 -51.93 -0.42 0.38 16 37 H 0.07 0.74 7.83 -51.93 -0.41 0.33 16 38 H 0.10 0.68 8.08 -51.93 -0.42 0.26 16 39 H 0.09 0.66 8.03 -51.93 -0.42 0.24 16 40 H 0.09 0.54 8.14 -51.93 -0.42 0.12 16 41 H 0.09 0.58 8.14 -51.93 -0.42 0.16 16 42 H 0.09 0.70 8.14 -51.93 -0.42 0.28 16 43 H 0.07 1.01 8.14 -51.93 -0.42 0.58 16 44 H 0.10 1.62 7.07 -51.93 -0.37 1.25 16 45 H 0.06 1.02 8.14 -51.93 -0.42 0.60 16 46 H 0.07 1.36 8.14 -51.93 -0.42 0.94 16 47 H 0.09 2.31 6.06 -51.93 -0.31 2.00 16 48 H 0.03 0.72 7.83 -51.92 -0.41 0.32 16 49 H 0.25 7.15 8.31 -40.82 -0.34 6.81 16 50 H 0.05 1.06 7.93 -51.93 -0.41 0.65 16 51 H 0.02 0.44 7.83 -51.93 -0.41 0.04 16 52 H 0.07 1.31 8.14 -51.93 -0.42 0.89 16 53 H 0.05 0.83 6.15 -51.93 -0.32 0.51 16 LS Contribution 402.78 15.07 6.07 6.07 Total: -1.00 -32.19 402.78 -10.79 -42.98 By element: Atomic # 1 Polarization: 10.14 SS G_CDS: -9.54 Total: 0.59 kcal Atomic # 6 Polarization: 0.30 SS G_CDS: -1.74 Total: -1.44 kcal Atomic # 7 Polarization: -49.33 SS G_CDS: -0.33 Total: -49.66 kcal Atomic # 8 Polarization: -15.20 SS G_CDS: -0.22 Total: -15.42 kcal Atomic # 14 Polarization: 21.91 SS G_CDS: -5.03 Total: 16.88 kcal Total LS contribution 6.07 Total: 6.07 kcal Total: -32.19 -10.79 -42.98 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. ZINC001098577850.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 6.060 kcal (2) G-P(sol) polarization free energy of solvation -32.193 kcal (3) G-ENP(sol) elect.-nuc.-pol. free energy of system -26.132 kcal (4) G-CDS(sol) cavity-dispersion-solvent structure free energy -10.791 kcal (5) G-P-CDS(sol) = G-P(sol) + G-CDS(sol) = (2) + (4) -42.983 kcal (6) G-S(sol) free energy of system = (1) + (5) -36.923 kcal 1SCF WAS SPECIFIED, SO GEOMETRY OPTIMIZATION WAS NOT USED Total computer time = 1.43 seconds