Wall clock time and date at job start Tue Mar 31 2020 21:10:23 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 N 2 2 C 1.31620 * 1 3 3 Si 1.86795 * 120.00111 * 2 1 4 4 H 1.48501 * 109.46990 * 270.00007 * 3 2 1 5 5 H 1.48502 * 109.47320 * 29.99829 * 3 2 1 6 6 H 1.48500 * 109.47306 * 150.00394 * 3 2 1 7 7 C 1.38818 * 120.00015 * 180.02562 * 2 1 3 8 8 H 1.07999 * 120.00064 * 164.75871 * 7 2 1 9 9 N 1.36936 * 119.99748 * 344.75239 * 7 2 1 10 10 C 1.46929 * 120.63432 * 331.60856 * 9 7 2 11 11 C 1.53190 * 108.77498 * 233.58800 * 10 9 7 12 12 C 1.53042 * 109.31338 * 305.36132 * 11 10 9 13 13 C 1.53038 * 109.53754 * 61.36657 * 12 11 10 14 14 H 1.09005 * 109.49601 * 58.62101 * 13 12 11 15 15 C 1.53002 * 109.49852 * 178.68337 * 13 12 11 16 16 H 1.08997 * 109.52214 * 59.92128 * 15 13 12 17 17 C 1.53042 * 109.49528 * 180.02562 * 15 13 12 18 18 C 1.53042 * 109.53428 * 181.31660 * 17 15 13 19 19 C 1.53190 * 109.31056 * 298.63515 * 18 17 15 20 20 N 1.46920 * 108.78016 * 54.63910 * 19 18 17 21 21 C 1.34777 * 120.63519 * 126.40569 * 20 19 18 22 22 O 1.21289 * 119.99974 * 173.84996 * 21 20 19 23 23 C 1.50702 * 119.99884 * 353.85764 * 21 20 19 24 24 C 1.53000 * 115.54845 * 291.75318 * 23 21 20 25 25 F 1.39895 * 109.47147 * 294.32906 * 24 23 21 26 26 C 1.53002 * 117.49797 * 146.08409 * 23 21 20 27 27 C 1.53000 * 117.49437 * 77.42562 * 23 21 20 28 28 C 1.46925 * 118.73607 * 306.38785 * 20 19 18 29 29 C 1.46926 * 120.62936 * 151.58905 * 9 7 2 30 30 H 0.96996 * 120.00141 * 0.02562 * 1 2 3 31 31 H 1.08993 * 109.58520 * 113.79710 * 10 9 7 32 32 H 1.09006 * 109.70235 * 353.44365 * 10 9 7 33 33 H 1.09006 * 109.49423 * 185.41314 * 11 10 9 34 34 H 1.08995 * 109.52329 * 65.33556 * 11 10 9 35 35 H 1.09004 * 109.46141 * 301.34817 * 12 11 10 36 36 H 1.08999 * 109.44907 * 181.38217 * 12 11 10 37 37 H 1.08997 * 109.45919 * 61.30026 * 17 15 13 38 38 H 1.08992 * 109.44828 * 301.33064 * 17 15 13 39 39 H 1.09000 * 109.49720 * 58.58419 * 18 17 15 40 40 H 1.08998 * 109.49372 * 178.68460 * 18 17 15 41 41 H 1.09009 * 109.58268 * 174.42556 * 19 18 17 42 42 H 1.08996 * 109.58951 * 294.84365 * 19 18 17 43 43 H 1.08999 * 109.47012 * 54.32905 * 24 23 21 44 44 H 1.09004 * 109.46964 * 174.32568 * 24 23 21 45 45 H 1.09005 * 117.49641 * 359.97438 * 26 23 21 46 46 H 1.08997 * 117.49815 * 145.01611 * 26 23 21 47 47 H 1.09007 * 117.49264 * 214.97613 * 27 23 21 48 48 H 1.08995 * 117.50333 * 359.97438 * 27 23 21 49 49 H 1.09000 * 109.59146 * 293.81964 * 28 20 19 50 50 H 1.09000 * 109.70778 * 173.46260 * 28 20 19 51 51 H 1.09001 * 109.58962 * 6.61981 * 29 9 7 52 52 H 1.09008 * 109.58257 * 246.20311 * 29 9 7 Note: An asterisk (*) indicates this parameter was or will be optimized. Cartesian coordinates (angstroms) NO. Atom X Y Z 1 7 0.0000 0.0000 0.0000 2 6 1.3162 0.0000 0.0000 3 14 2.2502 1.6177 0.0000 4 1 2.4977 2.0463 -1.4001 5 1 1.4477 2.6527 0.7000 6 1 3.5478 1.4401 0.6999 7 6 2.0103 -1.2022 -0.0005 8 1 3.0618 -1.2185 -0.2464 9 7 1.3617 -2.3673 0.3108 10 6 0.1978 -2.3597 1.2076 11 6 0.4433 -3.3590 2.3425 12 6 0.7692 -4.7308 1.7474 13 6 2.0434 -4.6337 0.9053 14 1 2.8655 -4.2820 1.5287 15 6 2.3867 -6.0125 0.3377 16 1 1.5647 -6.3647 -0.2855 17 6 3.6613 -5.9154 -0.5039 18 6 3.9872 -7.2872 -1.0990 19 6 4.2320 -8.2866 0.0360 20 7 3.0678 -8.2789 0.9322 21 6 2.4291 -9.4256 1.2382 22 8 1.5411 -9.4233 2.0644 23 6 2.8184 -10.7107 0.5540 24 6 2.4866 -10.7909 -0.9374 25 9 1.0981 -10.7820 -1.1074 26 6 2.7238 -12.0019 1.3693 27 6 4.0937 -11.4075 1.0327 28 6 2.6133 -6.9974 1.4889 29 6 1.8162 -3.6489 -0.2458 30 1 -0.4850 0.8400 -0.0004 31 1 -0.6936 -2.6499 0.6516 32 1 0.0605 -1.3614 1.6230 33 1 -0.4512 -3.4339 2.9610 34 1 1.2800 -3.0189 2.9526 35 1 -0.0578 -5.0594 1.1179 36 1 0.9202 -5.4493 2.5530 37 1 4.4881 -5.5871 0.1259 38 1 3.5105 -5.1967 -1.3093 39 1 3.1503 -7.6273 -1.7090 40 1 4.8819 -7.2123 -1.7170 41 1 4.3668 -9.2854 -0.3794 42 1 5.1233 -7.9967 0.5923 43 1 2.9213 -9.9345 -1.4529 44 1 2.8976 -11.7114 -1.3520 45 1 2.3803 -11.9218 2.4007 46 1 2.4254 -12.9086 0.8431 47 1 4.6961 -11.9236 0.2850 48 1 4.6512 -10.9364 1.8421 49 1 3.3727 -6.6022 2.1637 50 1 1.6807 -7.1435 2.0338 51 1 2.7492 -3.5008 -0.7896 52 1 1.0566 -4.0441 -0.9204 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) F: (AM1): M.J.S. DEWAR AND E. G. ZOEBISCH, THEOCHEM, 180, 1 (1988). 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 ZINC001053038474.mol2 52 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:10:23 Heat of formation + Delta-G solvation = -53.365793 kcal Electronic energy + Delta-G solvation = -31173.717944 eV Core-core repulsion = 27088.125839 eV Total energy + Delta-G solvation = -4085.592105 eV No. of doubly occupied orbitals = 65 Molecular weight (most abundant/longest-lived isotopes) = 338.211 amu Computer time = 1.06 seconds Orbital eigenvalues (eV) -47.89202 -40.72036 -38.74593 -37.08011 -36.67793 -33.94121 -32.61462 -31.48122 -30.39820 -28.77387 -27.90309 -26.15492 -25.08196 -23.79511 -23.07814 -22.80324 -21.65841 -21.32367 -19.87013 -19.15187 -18.73827 -17.45152 -16.69475 -16.60428 -16.28411 -15.89931 -15.54963 -14.94933 -14.76929 -14.67212 -14.45119 -14.18774 -13.97645 -13.71251 -13.41744 -13.28650 -13.05627 -12.95036 -12.55213 -12.39402 -12.20707 -12.10255 -11.98305 -11.77140 -11.61382 -11.36770 -11.17138 -10.93916 -10.92528 -10.81731 -10.70132 -10.61473 -10.48001 -10.26550 -10.12559 -10.05753 -9.88133 -9.79808 -9.36413 -9.28727 -8.71141 -8.42602 -6.60365 -5.74496 -3.18128 2.33321 2.73418 3.11852 3.38293 3.45246 3.47268 3.86751 4.08143 4.21640 4.33257 4.50726 4.53953 4.62858 4.77111 4.80251 4.88214 4.98938 5.07512 5.13254 5.16120 5.23572 5.32603 5.40752 5.42377 5.45646 5.54240 5.61500 5.64760 5.65472 5.72486 5.79092 5.88167 5.95411 6.01904 6.13526 6.18264 6.31404 6.37200 6.46287 6.49106 6.56336 6.71962 6.78135 6.89920 6.97450 7.07674 7.24407 7.47779 7.74583 7.82549 8.00946 8.49844 9.25792 9.89555 10.46857 11.33556 Molecular weight = 338.21amu Principal moments of inertia in cm(-1) A = 0.022982 B = 0.002597 C = 0.002509 Principal moments of inertia in units of 10**(-40)*gram-cm**2 A = 1218.044860 B =10779.912573 C =11155.341632 Net atomic charges, atomic populations, and dipole contributions calculated with CM2 Atom NO. Type Charge No. of electrons 1 N -0.904 5.904 2 C 0.001 3.999 3 Si 0.902 3.098 4 H -0.287 1.287 5 H -0.290 1.290 6 H -0.280 1.280 7 C -0.250 4.250 8 H 0.064 0.936 9 N -0.592 5.592 10 C 0.170 3.830 11 C -0.142 4.142 12 C -0.106 4.106 13 C -0.094 4.094 14 H 0.067 0.933 15 C -0.092 4.092 16 H 0.071 0.929 17 C -0.117 4.117 18 C -0.137 4.137 19 C 0.105 3.895 20 N -0.605 5.605 21 C 0.552 3.448 22 O -0.530 6.530 23 C -0.202 4.202 24 C 0.105 3.895 25 F -0.206 7.206 26 C -0.134 4.134 27 C -0.150 4.150 28 C 0.116 3.884 29 C 0.156 3.844 30 H 0.253 0.747 31 H 0.015 0.985 32 H 0.087 0.913 33 H 0.047 0.953 34 H 0.052 0.948 35 H 0.053 0.947 36 H 0.045 0.955 37 H 0.064 0.936 38 H 0.076 0.924 39 H 0.066 0.934 40 H 0.075 0.925 41 H 0.082 0.918 42 H 0.071 0.929 43 H 0.084 0.916 44 H 0.085 0.915 45 H 0.108 0.892 46 H 0.106 0.894 47 H 0.109 0.891 48 H 0.104 0.896 49 H 0.066 0.934 50 H 0.094 0.906 51 H 0.053 0.947 52 H 0.015 0.985 For the charges calculated by CM2: Dipole moment (debyes) X Y Z Total from point charges 10.865 -23.628 -0.558 26.013 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 N -0.548 5.548 2 C -0.197 4.197 3 Si 0.731 3.269 4 H -0.214 1.214 5 H -0.216 1.216 6 H -0.205 1.205 7 C -0.380 4.380 8 H 0.081 0.919 9 N -0.314 5.314 10 C 0.044 3.956 11 C -0.181 4.181 12 C -0.144 4.144 13 C -0.115 4.115 14 H 0.086 0.914 15 C -0.112 4.112 16 H 0.090 0.910 17 C -0.154 4.154 18 C -0.175 4.175 19 C -0.018 4.018 20 N -0.337 5.337 21 C 0.339 3.661 22 O -0.406 6.406 23 C -0.206 4.206 24 C 0.050 3.950 25 F -0.187 7.187 26 C -0.171 4.171 27 C -0.187 4.187 28 C -0.005 4.005 29 C 0.029 3.971 30 H 0.063 0.937 31 H 0.034 0.966 32 H 0.105 0.895 33 H 0.065 0.935 34 H 0.071 0.929 35 H 0.071 0.929 36 H 0.064 0.936 37 H 0.083 0.917 38 H 0.094 0.906 39 H 0.084 0.916 40 H 0.094 0.906 41 H 0.101 0.899 42 H 0.089 0.911 43 H 0.102 0.898 44 H 0.103 0.897 45 H 0.126 0.874 46 H 0.125 0.875 47 H 0.127 0.873 48 H 0.123 0.877 49 H 0.085 0.915 50 H 0.112 0.888 51 H 0.071 0.929 52 H 0.033 0.967 Dipole moment (debyes) X Y Z Total from point charges 10.240 -23.445 -0.120 25.584 hybrid contribution 0.352 0.964 -0.115 1.033 sum 10.592 -22.481 -0.236 24.853 Atomic orbital electron populations 1.70020 1.05250 1.25382 1.54181 1.25059 0.95303 0.98335 1.01052 0.86312 0.79113 0.83773 0.77654 1.21449 1.21594 1.20541 1.19158 0.91921 0.84608 1.42282 0.91858 1.44035 1.34043 1.01715 1.51656 1.20787 0.88256 0.98056 0.88492 1.22179 0.99824 0.96335 0.99793 1.21497 0.97247 0.97353 0.98261 1.21651 0.97783 0.93138 0.98950 0.91433 1.21527 0.97871 0.95556 0.96237 0.91031 1.21622 0.97627 0.96130 1.00047 1.21979 1.02183 0.96095 0.97291 1.21842 0.88040 1.00419 0.91511 1.48069 1.30770 1.06988 1.47906 1.19274 0.81463 0.83535 0.81833 1.90709 1.30872 1.85628 1.33366 1.22065 1.07057 0.94549 0.96882 1.21493 0.77423 1.04057 0.92031 1.92829 1.31547 1.98130 1.96176 1.23033 0.98512 0.94102 1.01500 1.23125 0.90132 1.01945 1.03492 1.21824 1.01202 0.81510 0.95972 1.20913 0.97289 0.86640 0.92257 0.93680 0.96647 0.89528 0.93476 0.92909 0.92867 0.93613 0.91709 0.90587 0.91568 0.90617 0.89937 0.91085 0.89754 0.89718 0.87380 0.87522 0.87320 0.87724 0.91523 0.88772 0.92875 0.96749 Number of geometries 1 Number of calculations of the screened coulomb radii 1 The total number of SCF iterations 17. 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 N -0.90 -26.56 11.10 55.49 0.62 -25.94 16 2 C 0.00 0.02 4.96 -17.43 -0.09 -0.06 16 3 Si 0.90 21.77 29.58 -169.99 -5.03 16.75 16 4 H -0.29 -6.92 7.11 56.52 0.40 -6.52 16 5 H -0.29 -7.05 7.11 56.52 0.40 -6.65 16 6 H -0.28 -6.62 7.11 56.52 0.40 -6.21 16 7 C -0.25 -6.88 9.59 -15.06 -0.14 -7.02 16 8 H 0.06 1.70 7.87 -52.49 -0.41 1.29 16 9 N -0.59 -14.95 2.99 -172.31 -0.52 -15.46 16 10 C 0.17 4.17 5.29 -3.71 -0.02 4.15 16 11 C -0.14 -2.84 6.08 -26.61 -0.16 -3.00 16 12 C -0.11 -1.87 4.86 -26.69 -0.13 -2.00 16 13 C -0.09 -1.60 1.92 -90.50 -0.17 -1.77 16 14 H 0.07 1.19 8.14 -51.93 -0.42 0.77 16 15 C -0.09 -1.26 1.92 -90.50 -0.17 -1.43 16 16 H 0.07 1.08 8.14 -51.93 -0.42 0.65 16 17 C -0.12 -1.34 4.86 -26.69 -0.13 -1.47 16 18 C -0.14 -1.18 6.03 -26.61 -0.16 -1.34 16 19 C 0.10 0.84 5.67 -3.72 -0.02 0.82 16 20 N -0.60 -6.93 2.97 -172.77 -0.51 -7.44 16 21 C 0.55 7.15 6.91 -10.98 -0.08 7.07 16 22 O -0.53 -9.12 16.22 5.55 0.09 -9.03 16 23 C -0.20 -1.81 1.07 -155.66 -0.17 -1.97 16 24 C 0.10 0.82 6.98 37.16 0.26 1.08 16 25 F -0.21 -2.50 17.33 2.25 0.04 -2.46 16 26 C -0.13 -1.03 9.79 -26.69 -0.26 -1.29 16 27 C -0.15 -0.85 9.12 -26.69 -0.24 -1.09 16 28 C 0.12 1.53 5.17 -3.71 -0.02 1.51 16 29 C 0.16 3.21 5.28 -3.71 -0.02 3.19 16 30 H 0.25 7.19 8.31 -40.82 -0.34 6.85 16 31 H 0.02 0.40 8.14 -51.93 -0.42 -0.02 16 32 H 0.09 2.37 6.03 -51.93 -0.31 2.06 16 33 H 0.05 0.86 8.14 -51.93 -0.42 0.44 16 34 H 0.05 1.06 8.14 -51.93 -0.42 0.64 16 35 H 0.05 0.97 8.14 -51.93 -0.42 0.55 16 36 H 0.05 0.74 6.49 -51.93 -0.34 0.40 16 37 H 0.06 0.73 8.14 -51.93 -0.42 0.30 16 38 H 0.08 0.95 6.50 -51.93 -0.34 0.62 16 39 H 0.07 0.61 7.89 -51.93 -0.41 0.20 16 40 H 0.08 0.51 8.14 -51.93 -0.42 0.09 16 41 H 0.08 0.48 4.53 -51.92 -0.24 0.25 16 42 H 0.07 0.51 8.14 -51.93 -0.42 0.09 16 43 H 0.08 0.57 6.18 -51.93 -0.32 0.25 16 44 H 0.08 0.40 8.09 -51.93 -0.42 -0.02 16 45 H 0.11 0.99 8.05 -51.93 -0.42 0.57 16 46 H 0.11 0.64 8.14 -51.93 -0.42 0.21 16 47 H 0.11 0.34 8.14 -51.92 -0.42 -0.08 16 48 H 0.10 0.62 8.14 -51.93 -0.42 0.19 16 49 H 0.07 0.84 8.14 -51.93 -0.42 0.41 16 50 H 0.09 1.45 5.64 -51.93 -0.29 1.15 16 51 H 0.05 1.06 6.52 -51.93 -0.34 0.72 16 52 H 0.01 0.31 8.14 -51.92 -0.42 -0.11 16 LS Contribution 393.10 15.07 5.92 5.92 Total: -1.00 -33.22 393.10 -10.01 -43.22 By element: Atomic # 1 Polarization: 7.96 SS G_CDS: -8.89 Total: -0.93 kcal Atomic # 6 Polarization: -2.91 SS G_CDS: -1.73 Total: -4.64 kcal Atomic # 7 Polarization: -48.43 SS G_CDS: -0.41 Total: -48.84 kcal Atomic # 8 Polarization: -9.12 SS G_CDS: 0.09 Total: -9.03 kcal Atomic # 9 Polarization: -2.50 SS G_CDS: 0.04 Total: -2.46 kcal Atomic # 14 Polarization: 21.77 SS G_CDS: -5.03 Total: 16.75 kcal Total LS contribution 5.92 Total: 5.92 kcal Total: -33.22 -10.01 -43.22 kcal The number of atoms in the molecule is 52 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. ZINC001053038474.mol2 52 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 -10.142 kcal (2) G-P(sol) polarization free energy of solvation -33.219 kcal (3) G-ENP(sol) elect.-nuc.-pol. free energy of system -43.361 kcal (4) G-CDS(sol) cavity-dispersion-solvent structure free energy -10.005 kcal (5) G-P-CDS(sol) = G-P(sol) + G-CDS(sol) = (2) + (4) -43.224 kcal (6) G-S(sol) free energy of system = (1) + (5) -53.366 kcal 1SCF WAS SPECIFIED, SO GEOMETRY OPTIMIZATION WAS NOT USED Total computer time = 1.06 seconds