Wall clock time and date at job start Sat Apr 11 2020 06:09:43 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.53001 * 1 3 3 C 1.53005 * 109.47114 * 2 1 4 4 H 1.09004 * 109.87070 * 54.13579 * 3 2 1 5 5 C 1.54325 * 109.88537 * 293.12016 * 3 2 1 6 6 H 1.08994 * 110.71985 * 100.02191 * 5 3 2 7 7 C 1.50703 * 110.86692 * 336.82069 * 5 3 2 8 8 O 1.21292 * 120.00379 * 301.14537 * 7 5 3 9 9 N 1.34776 * 119.99804 * 121.15481 * 7 5 3 10 10 C 1.37101 * 120.00187 * 184.57405 * 9 7 5 11 11 H 1.08006 * 119.99654 * 13.82483 * 10 9 7 12 12 C 1.38945 * 120.00137 * 193.82314 * 10 9 7 13 13 N 1.31411 * 120.00246 * 12.17814 * 12 10 9 14 14 Si 1.86799 * 119.99988 * 192.17350 * 12 10 9 15 15 H 1.48503 * 109.46948 * 90.00201 * 14 12 10 16 16 H 1.48498 * 109.47167 * 210.00250 * 14 12 10 17 17 H 1.48499 * 109.47266 * 330.00354 * 14 12 10 18 18 C 1.46499 * 119.99888 * 4.56461 * 9 7 5 19 19 C 1.53002 * 109.46973 * 83.10888 * 18 9 7 20 20 C 1.55154 * 102.94085 * 218.38934 * 5 3 2 21 21 C 1.54912 * 101.62040 * 35.48069 * 20 5 3 22 22 N 1.47025 * 109.88430 * 175.30782 * 3 2 1 23 23 C 1.34780 * 125.64853 * 298.39444 * 22 3 2 24 24 O 1.21513 * 119.99661 * 4.96138 * 23 22 3 25 25 O 1.34630 * 119.99669 * 184.95507 * 23 22 3 26 26 C 1.45207 * 116.99620 * 180.02562 * 25 23 22 27 27 C 1.53001 * 109.47328 * 59.99782 * 26 25 23 28 28 C 1.53002 * 109.47162 * 180.02562 * 26 25 23 29 29 C 1.53001 * 109.47102 * 299.99810 * 26 25 23 30 30 H 1.09003 * 109.46501 * 305.44477 * 1 2 3 31 31 H 1.08997 * 109.47173 * 65.44487 * 1 2 3 32 32 H 1.08993 * 109.46670 * 185.44628 * 1 2 3 33 33 H 1.09005 * 109.47320 * 240.00673 * 2 1 3 34 34 H 1.09002 * 109.47709 * 120.00372 * 2 1 3 35 35 H 0.96992 * 120.00213 * 179.97438 * 13 12 10 36 36 H 1.09000 * 109.46970 * 203.10218 * 18 9 7 37 37 H 1.08999 * 109.47643 * 323.10510 * 18 9 7 38 38 H 1.08999 * 109.47268 * 183.76445 * 19 18 9 39 39 H 1.08998 * 109.46660 * 303.76223 * 19 18 9 40 40 H 1.09000 * 109.47310 * 63.76323 * 19 18 9 41 41 H 1.09002 * 111.00509 * 153.57299 * 20 5 3 42 42 H 1.09001 * 111.00349 * 277.38742 * 20 5 3 43 43 H 1.08998 * 110.36504 * 81.84805 * 21 20 5 44 44 H 1.08998 * 110.36646 * 204.17364 * 21 20 5 45 45 H 1.09002 * 109.46571 * 60.00182 * 27 26 25 46 46 H 1.08996 * 109.47273 * 179.97438 * 27 26 25 47 47 H 1.08997 * 109.47308 * 300.00266 * 27 26 25 48 48 H 1.08994 * 109.47191 * 59.99699 * 28 26 25 49 49 H 1.08996 * 109.47289 * 180.02562 * 28 26 25 50 50 H 1.09005 * 109.46761 * 299.99560 * 28 26 25 51 51 H 1.09006 * 109.46794 * 59.99889 * 29 26 25 52 52 H 1.09002 * 109.47427 * 180.02562 * 29 26 25 53 53 H 1.08994 * 109.47322 * 300.00345 * 29 26 25 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 6 2.0400 1.4425 0.0000 4 1 1.5973 1.9921 0.8308 5 6 1.6777 2.1274 -1.3347 6 1 0.7947 2.7557 -1.2185 7 6 1.4635 1.1056 -2.4215 8 8 2.3564 0.3406 -2.7193 9 7 0.2798 1.0420 -3.0628 10 6 0.0521 0.0535 -3.9851 11 1 0.7160 -0.7966 -4.0407 12 6 -1.0329 0.1452 -4.8482 13 7 -1.6977 1.2742 -4.9504 14 14 -1.5466 -1.3362 -5.8635 15 1 -2.5372 -2.1387 -5.1020 16 1 -2.1562 -0.8762 -7.1371 17 1 -0.3547 -2.1715 -6.1585 18 6 -0.7635 2.0291 -2.7742 19 6 -0.5076 3.2959 -3.5931 20 6 2.9312 2.9918 -1.6331 21 6 4.0685 2.0890 -1.0935 22 7 3.5059 1.4547 0.1131 23 6 4.2179 0.9668 1.1483 24 8 3.6536 0.3829 2.0522 25 8 5.5543 1.1264 1.1806 26 6 6.2565 0.5748 2.3256 27 6 6.0212 -0.9357 2.3890 28 6 7.7547 0.8527 2.1875 29 6 5.7335 1.2267 3.6072 30 1 -0.3632 0.5960 0.8373 31 1 -0.3633 0.4271 -0.9347 32 1 -0.3632 -1.0230 0.0975 33 1 1.8934 -0.5137 -0.8901 34 1 1.8935 -0.5139 0.8899 35 1 -2.4548 1.3383 -5.5532 36 1 -1.7372 1.6163 -3.0383 37 1 -0.7486 2.2738 -1.7121 38 1 -1.2446 4.0544 -3.3290 39 1 0.4930 3.6713 -3.3784 40 1 -0.5888 3.0647 -4.6552 41 1 3.0424 3.1620 -2.7040 42 1 2.8921 3.9368 -1.0913 43 1 4.3340 1.3317 -1.8311 44 1 4.9401 2.6904 -0.8353 45 1 6.3938 -1.4000 1.4760 46 1 6.5487 -1.3498 3.2483 47 1 4.9538 -1.1333 2.4877 48 1 7.9220 1.9288 2.1425 49 1 8.2819 0.4382 3.0467 50 1 8.1277 0.3885 1.2744 51 1 4.6661 1.0291 3.7058 52 1 6.2609 0.8129 4.4667 53 1 5.9011 2.3027 3.5620 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 ZINC000496866060.mol2 53 1SCF WAS SPECIFIED, SO GEOMETRY OPTIMIZATION WAS NOT USED SCF FIELD WAS ACHIEVED Wall clock time and date at job start Sat Apr 11 2020 06:09:43 Heat of formation + Delta-G solvation = -95.550595 kcal Electronic energy + Delta-G solvation = -32989.391482 eV Core-core repulsion = 28862.285838 eV Total energy + Delta-G solvation = -4127.105643 eV No. of doubly occupied orbitals = 66 Molecular weight (most abundant/longest-lived isotopes) = 340.220 amu Computer time = 0.89 seconds Orbital eigenvalues (eV) -40.81890 -39.73486 -37.59936 -36.16838 -34.96198 -34.09833 -32.87474 -31.94596 -31.03303 -29.42242 -27.68720 -27.46841 -27.40633 -25.46911 -24.02736 -22.34476 -22.21520 -21.78093 -20.71281 -19.11602 -18.20035 -17.59304 -16.96451 -16.68541 -16.22712 -15.96511 -15.12423 -15.00719 -14.73991 -14.39675 -14.18686 -13.84424 -13.71263 -13.27792 -13.02886 -12.79032 -12.70894 -12.56493 -12.48098 -12.41914 -12.08763 -11.93661 -11.89642 -11.80181 -11.75927 -11.46414 -11.29719 -11.20203 -11.17423 -11.13140 -11.01789 -10.84746 -10.80913 -10.36770 -10.31249 -10.31047 -10.16647 -10.05868 -9.86287 -9.67614 -8.90896 -8.37693 -8.26159 -7.70157 -6.28760 -3.85480 2.29215 3.09193 3.13701 3.24211 3.29051 3.76544 4.07070 4.09020 4.67508 4.70572 4.72532 4.85213 4.85563 4.91721 5.18630 5.31273 5.36951 5.41723 5.46012 5.56135 5.59965 5.61486 5.67902 5.70077 5.80919 5.81566 5.87722 5.89795 5.92168 5.95775 6.12433 6.18008 6.23495 6.25994 6.30712 6.34461 6.44901 6.50852 6.55827 6.59275 6.62927 6.64915 6.76604 6.90274 7.02818 7.32745 7.49431 7.72506 7.80296 8.08198 8.50729 8.58830 9.16966 9.42151 9.96985 10.66537 Molecular weight = 340.22amu Principal moments of inertia in cm(-1) A = 0.020468 B = 0.003141 C = 0.003040 Principal moments of inertia in units of 10**(-40)*gram-cm**2 A = 1367.686994 B = 8911.043165 C = 9206.862714 Net atomic charges, atomic populations, and dipole contributions calculated with CM2 Atom NO. Type Charge No. of electrons 1 C -0.154 4.154 2 C -0.151 4.151 3 C 0.181 3.819 4 H 0.074 0.926 5 C -0.119 4.119 6 H 0.103 0.897 7 C 0.462 3.538 8 O -0.584 6.584 9 N -0.438 5.438 10 C -0.293 4.293 11 H 0.101 0.899 12 C 0.025 3.975 13 N -0.878 5.878 14 Si 0.925 3.075 15 H -0.279 1.279 16 H -0.285 1.285 17 H -0.271 1.271 18 C 0.115 3.885 19 C -0.148 4.148 20 C -0.113 4.113 21 C 0.108 3.892 22 N -0.590 5.590 23 C 0.668 3.332 24 O -0.572 6.572 25 O -0.363 6.363 26 C 0.135 3.865 27 C -0.183 4.183 28 C -0.139 4.139 29 C -0.182 4.182 30 H 0.044 0.956 31 H 0.074 0.926 32 H 0.052 0.948 33 H 0.093 0.907 34 H 0.070 0.930 35 H 0.265 0.735 36 H 0.113 0.887 37 H 0.030 0.970 38 H 0.035 0.965 39 H 0.022 0.978 40 H 0.091 0.909 41 H 0.087 0.913 42 H 0.075 0.925 43 H 0.087 0.913 44 H 0.071 0.929 45 H 0.059 0.941 46 H 0.066 0.934 47 H 0.089 0.911 48 H 0.065 0.935 49 H 0.077 0.923 50 H 0.066 0.934 51 H 0.089 0.911 52 H 0.067 0.933 53 H 0.057 0.943 For the charges calculated by CM2: Dipole moment (debyes) X Y Z Total from point charges 9.410 4.457 11.570 15.565 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.212 4.212 2 C -0.189 4.189 3 C 0.079 3.921 4 H 0.092 0.908 5 C -0.141 4.141 6 H 0.121 0.879 7 C 0.251 3.749 8 O -0.470 6.470 9 N -0.155 5.155 10 C -0.421 4.421 11 H 0.119 0.881 12 C -0.176 4.176 13 N -0.519 5.519 14 Si 0.753 3.247 15 H -0.205 1.205 16 H -0.211 1.211 17 H -0.196 1.196 18 C -0.008 4.008 19 C -0.206 4.206 20 C -0.152 4.152 21 C -0.014 4.014 22 N -0.330 5.330 23 C 0.424 3.576 24 O -0.463 6.463 25 O -0.278 6.278 26 C 0.099 3.901 27 C -0.240 4.240 28 C -0.197 4.197 29 C -0.239 4.239 30 H 0.063 0.937 31 H 0.092 0.908 32 H 0.071 0.929 33 H 0.111 0.889 34 H 0.088 0.912 35 H 0.076 0.924 36 H 0.131 0.869 37 H 0.048 0.952 38 H 0.054 0.946 39 H 0.041 0.959 40 H 0.110 0.890 41 H 0.106 0.894 42 H 0.094 0.906 43 H 0.105 0.895 44 H 0.090 0.910 45 H 0.078 0.922 46 H 0.085 0.915 47 H 0.108 0.892 48 H 0.084 0.916 49 H 0.096 0.904 50 H 0.085 0.915 51 H 0.107 0.893 52 H 0.085 0.915 53 H 0.076 0.924 Dipole moment (debyes) X Y Z Total from point charges 9.601 3.910 11.667 15.607 hybrid contribution -0.421 -0.498 -0.481 0.810 sum 9.180 3.412 11.186 14.867 Atomic orbital electron populations 1.21862 0.95217 1.01708 1.02428 1.22147 0.96321 0.95865 1.04572 1.21397 0.77033 0.97530 0.96181 0.90784 1.22503 0.99628 0.97652 0.94291 0.87930 1.19376 0.82947 0.86955 0.85650 1.90492 1.42244 1.45248 1.68988 1.43810 1.10771 1.25899 1.35067 1.19754 1.10910 1.00107 1.11350 0.88130 1.25230 0.97327 0.97654 0.97386 1.69997 1.21190 1.23335 1.37399 0.85924 0.78915 0.81127 0.78759 1.20482 1.21063 1.19568 1.21621 0.95474 0.89580 0.94150 1.22177 0.97433 0.95950 1.04996 1.22620 0.93297 0.96968 1.02296 1.22391 0.96094 0.96312 0.86625 1.47197 1.06916 1.60916 1.17950 1.17419 0.80433 0.78180 0.81609 1.90849 1.71068 1.46184 1.38158 1.86308 1.19231 1.74946 1.47325 1.21998 0.91158 0.93575 0.83417 1.22508 1.05345 0.93970 1.02210 1.21841 0.91957 1.02859 1.03000 1.22510 1.04667 1.00544 0.96182 0.93690 0.90750 0.92918 0.88877 0.91151 0.92387 0.86913 0.95205 0.94579 0.95903 0.89022 0.89431 0.90629 0.89496 0.91038 0.92178 0.91472 0.89240 0.91617 0.90384 0.91490 0.89274 0.91454 0.92389 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.15 -2.74 9.56 37.16 0.36 -2.38 16 2 C -0.15 -2.89 5.04 -26.73 -0.13 -3.02 16 3 C 0.18 3.16 2.80 -66.95 -0.19 2.97 16 4 H 0.07 1.20 8.14 -51.93 -0.42 0.78 16 5 C -0.12 -2.09 2.21 -89.94 -0.20 -2.29 16 6 H 0.10 1.68 5.38 -51.93 -0.28 1.40 16 7 C 0.46 10.24 4.08 -10.98 -0.04 10.20 16 8 O -0.58 -14.81 13.30 5.55 0.07 -14.73 16 9 N -0.44 -10.22 1.86 -114.81 -0.21 -10.44 16 10 C -0.29 -7.63 9.44 -14.93 -0.14 -7.77 16 11 H 0.10 2.80 7.31 -52.48 -0.38 2.42 16 12 C 0.03 0.63 5.16 -17.47 -0.09 0.53 16 13 N -0.88 -21.34 8.07 55.50 0.45 -20.89 16 14 Si 0.93 20.03 29.57 -169.99 -5.03 15.01 16 15 H -0.28 -5.99 7.11 56.52 0.40 -5.59 16 16 H -0.28 -6.12 7.11 56.52 0.40 -5.72 16 17 H -0.27 -5.99 7.11 56.52 0.40 -5.59 16 18 C 0.12 2.39 4.05 -4.04 -0.02 2.37 16 19 C -0.15 -2.90 9.55 37.16 0.35 -2.54 16 20 C -0.11 -1.72 6.53 -24.58 -0.16 -1.88 16 21 C 0.11 1.72 6.47 -2.50 -0.02 1.71 16 22 N -0.59 -10.18 3.07 -168.71 -0.52 -10.70 16 23 C 0.67 11.58 7.74 54.06 0.42 12.00 16 24 O -0.57 -10.58 10.67 -19.28 -0.21 -10.78 16 25 O -0.36 -5.56 9.39 -40.30 -0.38 -5.94 16 26 C 0.14 1.66 1.13 -90.62 -0.10 1.56 16 27 C -0.18 -2.26 8.37 37.16 0.31 -1.95 16 28 C -0.14 -1.27 8.85 37.16 0.33 -0.94 16 29 C -0.18 -2.08 8.37 37.16 0.31 -1.76 16 30 H 0.04 0.69 8.14 -51.93 -0.42 0.27 16 31 H 0.07 1.44 4.00 -51.93 -0.21 1.24 16 32 H 0.05 0.91 8.14 -51.93 -0.42 0.48 16 33 H 0.09 2.10 5.50 -51.93 -0.29 1.82 16 34 H 0.07 1.38 6.56 -51.93 -0.34 1.04 16 35 H 0.27 6.21 8.31 -40.82 -0.34 5.87 16 36 H 0.11 2.59 5.74 -51.93 -0.30 2.29 16 37 H 0.03 0.55 5.41 -51.93 -0.28 0.26 16 38 H 0.04 0.62 8.14 -51.93 -0.42 0.19 16 39 H 0.02 0.41 7.72 -51.93 -0.40 0.01 16 40 H 0.09 2.08 6.26 -51.93 -0.33 1.76 16 41 H 0.09 1.37 8.03 -51.93 -0.42 0.96 16 42 H 0.08 0.95 8.14 -51.93 -0.42 0.52 16 43 H 0.09 1.60 7.06 -51.93 -0.37 1.23 16 44 H 0.07 0.98 7.99 -51.93 -0.41 0.56 16 45 H 0.06 0.74 8.14 -51.93 -0.42 0.32 16 46 H 0.07 0.65 8.14 -51.93 -0.42 0.23 16 47 H 0.09 1.34 5.88 -51.93 -0.31 1.04 16 48 H 0.06 0.57 8.14 -51.93 -0.42 0.14 16 49 H 0.08 0.54 8.14 -51.93 -0.42 0.12 16 50 H 0.07 0.61 8.14 -51.93 -0.42 0.19 16 51 H 0.09 1.24 5.88 -51.93 -0.31 0.94 16 52 H 0.07 0.61 8.14 -51.93 -0.42 0.19 16 53 H 0.06 0.62 8.14 -51.93 -0.42 0.20 16 LS Contribution 391.34 15.07 5.90 5.90 Total: -1.00 -28.44 391.34 -7.75 -36.20 By element: Atomic # 1 Polarization: 18.40 SS G_CDS: -8.82 Total: 9.58 kcal Atomic # 6 Polarization: 5.81 SS G_CDS: 0.99 Total: 6.80 kcal Atomic # 7 Polarization: -41.74 SS G_CDS: -0.28 Total: -42.03 kcal Atomic # 8 Polarization: -30.95 SS G_CDS: -0.51 Total: -31.46 kcal Atomic # 14 Polarization: 20.03 SS G_CDS: -5.03 Total: 15.01 kcal Total LS contribution 5.90 Total: 5.90 kcal Total: -28.44 -7.75 -36.20 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. ZINC000496866060.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 -59.355 kcal (2) G-P(sol) polarization free energy of solvation -28.442 kcal (3) G-ENP(sol) elect.-nuc.-pol. free energy of system -87.797 kcal (4) G-CDS(sol) cavity-dispersion-solvent structure free energy -7.753 kcal (5) G-P-CDS(sol) = G-P(sol) + G-CDS(sol) = (2) + (4) -36.196 kcal (6) G-S(sol) free energy of system = (1) + (5) -95.551 kcal 1SCF WAS SPECIFIED, SO GEOMETRY OPTIMIZATION WAS NOT USED Total computer time = 0.90 seconds