Wall clock time and date at job start Wed Apr 1 2020 04:22:44 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.17401 * 1 3 3 C 1.47199 * 179.97438 * 2 1 4 4 N 1.46905 * 109.46838 * 271.63633 * 3 2 1 5 5 C 1.46854 * 110.99690 * 170.00361 * 4 3 2 6 6 C 1.53031 * 109.52164 * 185.84038 * 5 4 3 7 7 C 1.53118 * 109.27792 * 300.81960 * 6 5 4 8 8 C 1.50708 * 109.52284 * 177.56260 * 7 6 5 9 9 O 1.21277 * 119.99886 * 336.94336 * 8 7 6 10 10 N 1.34780 * 119.99868 * 156.94867 * 8 7 6 11 11 C 1.47845 * 122.22849 * 355.59022 * 10 8 7 12 12 C 1.52711 * 111.06682 * 232.46667 * 11 10 8 13 13 C 1.53872 * 111.52287 * 311.47420 * 12 11 10 14 14 C 1.51082 * 107.33704 * 56.92444 * 13 12 11 15 15 H 1.08996 * 108.73504 * 50.83409 * 14 13 12 16 16 C 1.46103 * 122.73400 * 175.18376 * 10 8 7 17 17 H 1.09001 * 111.69330 * 240.23164 * 16 10 8 18 18 C 1.53437 * 115.85131 * 9.72408 * 16 10 8 19 19 C 1.55524 * 101.76109 * 152.09026 * 18 16 10 20 20 N 1.45024 * 119.19122 * 176.09931 * 14 13 12 21 21 C 1.36935 * 126.56281 * 330.59435 * 20 14 13 22 22 H 1.08003 * 120.00084 * 359.71807 * 21 20 14 23 23 C 1.38809 * 119.99920 * 180.02562 * 21 20 14 24 24 N 1.31616 * 120.00407 * 359.72211 * 23 21 20 25 25 Si 1.86808 * 119.99971 * 179.72271 * 23 21 20 26 26 H 1.48497 * 109.46844 * 89.99565 * 25 23 21 27 27 H 1.48498 * 109.47001 * 210.00029 * 25 23 21 28 28 H 1.48500 * 109.46823 * 329.99660 * 25 23 21 29 29 C 1.53116 * 109.21504 * 57.62289 * 7 6 5 30 30 C 1.46847 * 110.99868 * 294.21036 * 4 3 2 31 31 H 4.40885 * 5.52223 * 179.97438 * 4 1 2 32 32 H 1.09008 * 109.46977 * 151.64019 * 3 2 1 33 33 H 1.08996 * 109.47504 * 31.63860 * 3 2 1 34 34 H 1.08998 * 109.46530 * 65.82118 * 5 4 3 35 35 H 1.08999 * 109.45868 * 305.85295 * 5 4 3 36 36 H 1.08997 * 109.50548 * 180.87888 * 6 5 4 37 37 H 1.08998 * 109.53610 * 60.77904 * 6 5 4 38 38 H 1.08999 * 109.52364 * 297.68650 * 7 6 5 39 39 H 1.09000 * 109.15964 * 112.05663 * 11 10 8 40 40 H 1.08998 * 109.12010 * 352.84544 * 11 10 8 41 41 H 1.09004 * 109.07559 * 190.95156 * 12 11 10 42 42 H 1.09000 * 109.01619 * 71.96342 * 12 11 10 43 43 H 1.08999 * 109.88719 * 297.46470 * 13 12 11 44 44 H 1.09000 * 109.94787 * 176.41916 * 13 12 11 45 45 H 1.08999 * 110.98873 * 270.38928 * 18 16 10 46 46 H 1.09005 * 110.98545 * 34.19888 * 18 16 10 47 47 H 1.08999 * 110.08472 * 218.32199 * 19 18 16 48 48 H 1.09002 * 109.86765 * 97.06240 * 19 18 16 49 49 H 0.96995 * 120.00529 * 180.02562 * 24 23 21 50 50 H 1.09007 * 109.50244 * 62.32024 * 29 7 6 51 51 H 1.09001 * 109.50380 * 182.43696 * 29 7 6 52 52 H 1.09004 * 109.46566 * 294.17646 * 30 4 3 53 53 H 1.08992 * 109.46500 * 54.14558 * 30 4 3 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 7 3.1356 -0.0387 1.3845 5 6 4.5841 -0.2772 1.4238 6 6 5.0336 -0.4567 2.8755 7 6 4.6831 0.8027 3.6727 8 6 5.0692 0.6087 5.1165 9 8 5.1739 -0.5108 5.5709 10 7 5.2981 1.6793 5.9026 11 6 5.2781 3.0614 5.3779 12 6 6.5573 3.8082 5.7496 13 6 6.8888 3.6510 7.2439 14 6 6.9932 2.1709 7.5286 15 1 7.6658 1.7206 6.7987 16 6 5.5748 1.5648 7.3326 17 1 4.8310 2.0916 7.9304 18 6 5.8208 0.1431 7.8546 19 6 6.7474 0.4094 9.0750 20 7 7.3840 1.7315 8.8543 21 6 8.1909 2.4050 9.7321 22 1 8.6077 3.3625 9.4567 23 6 8.4782 1.8541 10.9733 24 7 7.9708 0.6869 11.3088 25 14 9.5852 2.7693 12.1678 26 1 8.7600 3.6545 13.0283 27 1 10.3128 1.7919 13.0165 28 1 10.5628 3.5890 11.4077 29 6 3.1764 1.0574 3.5757 30 6 2.7804 1.1915 2.1034 31 1 -1.0500 0.0002 0.0000 32 1 3.0089 0.9052 -0.4882 33 1 3.0098 -0.8741 -0.5390 34 1 5.1045 0.5741 0.9851 35 1 4.8193 -1.1779 0.8568 36 1 6.1112 -0.6182 2.9056 37 1 4.5245 -1.3165 3.3111 38 1 5.2238 1.6557 3.2625 39 1 4.4211 3.5901 5.7952 40 1 5.1849 3.0299 4.2924 41 1 6.4306 4.8667 5.5224 42 1 7.3835 3.4136 5.1581 43 1 6.0946 4.0896 7.8480 44 1 7.8372 4.1394 7.4678 45 1 4.8871 -0.3257 8.1653 46 1 6.3311 -0.4644 7.1071 47 1 7.5115 -0.3651 9.1409 48 1 6.1570 0.4290 9.9911 49 1 8.1713 0.3022 12.1763 50 1 2.6367 0.2223 4.0222 51 1 2.9279 1.9770 4.1055 52 1 3.3092 2.0361 1.6616 53 1 1.7057 1.3576 2.0301 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 ZINC001049917973.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 04:22:44 Heat of formation + Delta-G solvation = 62.784099 kcal Electronic energy + Delta-G solvation = -31665.442390 eV Core-core repulsion = 27732.462706 eV Total energy + Delta-G solvation = -3932.979684 eV No. of doubly occupied orbitals = 66 Molecular weight (most abundant/longest-lived isotopes) = 345.216 amu Computer time = 1.21 seconds Orbital eigenvalues (eV) -40.64219 -39.20465 -37.62589 -35.44318 -34.70851 -32.60494 -32.56952 -31.90932 -31.07506 -29.02129 -28.23699 -27.53496 -24.94310 -24.43999 -23.45806 -23.35754 -21.47514 -21.06879 -19.97468 -19.56013 -18.71715 -18.45756 -16.85554 -16.47231 -16.09792 -15.59864 -15.56183 -15.17635 -14.96051 -14.72639 -14.35768 -14.12856 -13.97320 -13.78005 -13.64312 -13.36788 -13.16076 -13.05542 -12.71656 -12.42637 -12.30177 -12.16468 -11.89788 -11.59047 -11.55277 -11.35776 -11.18299 -10.92393 -10.84372 -10.79006 -10.66415 -10.60363 -10.46233 -10.42415 -10.34489 -10.06873 -9.94080 -9.81527 -9.58672 -8.86023 -8.71254 -8.53284 -8.18127 -6.79732 -5.67170 -3.01575 2.37668 2.54109 2.86605 3.18281 3.41439 3.46146 3.48482 3.68821 4.06338 4.37901 4.51381 4.57282 4.60692 4.75982 4.89711 4.90376 5.01970 5.05195 5.11554 5.16369 5.23074 5.27602 5.34627 5.39057 5.45424 5.48469 5.55329 5.57684 5.63062 5.65715 5.74808 5.79463 5.82819 5.89735 5.99938 6.02725 6.10967 6.23361 6.28657 6.31652 6.42293 6.59221 6.64349 6.71544 6.75450 6.81891 6.93905 7.24277 7.31564 7.61427 7.69653 7.77780 7.88057 8.28322 8.50101 9.22197 9.86145 10.54825 11.44689 Molecular weight = 345.22amu Principal moments of inertia in cm(-1) A = 0.027503 B = 0.002336 C = 0.002259 Principal moments of inertia in units of 10**(-40)*gram-cm**2 A = 1017.833982 B =11983.029204 C =12392.424952 Net atomic charges, atomic populations, and dipole contributions calculated with CM2 Atom NO. Type Charge No. of electrons 1 C -0.198 4.198 2 C -0.180 4.180 3 C 0.162 3.838 4 N -0.536 5.536 5 C 0.058 3.942 6 C -0.099 4.099 7 C -0.109 4.109 8 C 0.522 3.478 9 O -0.529 6.529 10 N -0.601 5.601 11 C 0.109 3.891 12 C -0.129 4.129 13 C -0.116 4.116 14 C 0.159 3.841 15 H 0.017 0.983 16 C 0.116 3.884 17 H 0.080 0.920 18 C -0.134 4.134 19 C 0.172 3.828 20 N -0.593 5.593 21 C -0.225 4.225 22 H 0.078 0.922 23 C -0.008 4.008 24 N -0.919 5.919 25 Si 0.910 3.090 26 H -0.289 1.289 27 H -0.291 1.291 28 H -0.282 1.282 29 C -0.103 4.103 30 C 0.059 3.941 31 H 0.217 0.783 32 H 0.057 0.943 33 H 0.098 0.902 34 H 0.036 0.964 35 H 0.081 0.919 36 H 0.073 0.927 37 H 0.086 0.914 38 H 0.087 0.913 39 H 0.065 0.935 40 H 0.078 0.922 41 H 0.071 0.929 42 H 0.062 0.938 43 H 0.063 0.937 44 H 0.066 0.934 45 H 0.058 0.942 46 H 0.087 0.913 47 H 0.038 0.962 48 H 0.056 0.944 49 H 0.252 0.748 50 H 0.079 0.921 51 H 0.073 0.927 52 H 0.036 0.964 53 H 0.078 0.922 For the charges calculated by CM2: Dipole moment (debyes) X Y Z Total from point charges -11.534 3.525 -23.565 26.472 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.216 4.216 2 C -0.180 4.180 3 C 0.036 3.964 4 N -0.261 5.261 5 C -0.069 4.069 6 C -0.137 4.137 7 C -0.130 4.130 8 C 0.310 3.690 9 O -0.405 6.405 10 N -0.334 5.334 11 C -0.014 4.014 12 C -0.168 4.168 13 C -0.154 4.154 14 C 0.049 3.951 15 H 0.035 0.965 16 C 0.012 3.988 17 H 0.097 0.903 18 C -0.172 4.172 19 C 0.048 3.952 20 N -0.318 5.318 21 C -0.354 4.354 22 H 0.096 0.904 23 C -0.205 4.205 24 N -0.567 5.567 25 Si 0.741 3.259 26 H -0.217 1.217 27 H -0.217 1.217 28 H -0.208 1.208 29 C -0.141 4.141 30 C -0.068 4.068 31 H 0.234 0.766 32 H 0.075 0.925 33 H 0.116 0.884 34 H 0.055 0.945 35 H 0.099 0.901 36 H 0.092 0.908 37 H 0.105 0.895 38 H 0.105 0.895 39 H 0.083 0.917 40 H 0.096 0.904 41 H 0.090 0.910 42 H 0.081 0.919 43 H 0.082 0.918 44 H 0.084 0.916 45 H 0.077 0.923 46 H 0.105 0.895 47 H 0.056 0.944 48 H 0.074 0.926 49 H 0.063 0.937 50 H 0.098 0.902 51 H 0.091 0.909 52 H 0.054 0.946 53 H 0.097 0.903 Dipole moment (debyes) X Y Z Total from point charges -11.561 2.466 -22.854 25.730 hybrid contribution 0.779 1.323 0.326 1.569 sum -10.782 3.789 -22.528 25.261 Atomic orbital electron populations 1.24911 0.95920 1.01300 0.99494 1.22825 0.94187 1.00193 1.00839 1.18987 0.85234 1.01985 0.90225 1.61883 1.14481 1.36986 1.12742 1.22133 0.86115 1.00883 0.97762 1.21435 1.01544 0.97715 0.92985 1.21427 0.97840 1.00343 0.93422 1.20372 0.76882 0.84249 0.87520 1.90685 1.51264 1.23723 1.74871 1.48110 1.69927 1.05848 1.09469 1.21735 0.98522 0.82316 0.98819 1.21901 0.97256 0.99688 0.97962 1.21437 1.01927 0.98266 0.93814 1.20757 0.92326 0.92616 0.89361 0.96534 1.22491 0.97867 0.99681 0.78736 0.90256 1.23037 1.00440 0.94206 0.99533 1.21753 0.93425 0.86596 0.93454 1.44086 1.56149 1.18932 1.12597 1.18976 1.20987 1.01140 0.94257 0.90421 1.24809 1.00709 0.97143 0.97800 1.69947 1.47937 1.17294 1.21485 0.86151 0.80098 0.78540 0.81154 1.21651 1.21714 1.20773 1.21507 0.93713 1.03009 0.95892 1.22098 1.01068 0.89798 0.93870 0.76645 0.92477 0.88382 0.94543 0.90068 0.90832 0.89512 0.89529 0.91674 0.90420 0.91000 0.91902 0.91816 0.91551 0.92336 0.89477 0.94376 0.92590 0.93741 0.90220 0.90888 0.94561 0.90323 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.20 -1.54 19.25 29.55 0.57 -0.97 16 2 C -0.18 -1.57 12.34 -37.12 -0.46 -2.03 16 3 C 0.16 1.19 5.46 -6.61 -0.04 1.15 16 4 N -0.54 -4.70 4.64 -238.17 -1.10 -5.81 16 5 C 0.06 0.47 5.46 -3.83 -0.02 0.45 16 6 C -0.10 -1.04 4.93 -26.65 -0.13 -1.18 16 7 C -0.11 -1.18 1.85 -91.65 -0.17 -1.35 16 8 C 0.52 7.98 6.86 -10.98 -0.08 7.90 16 9 O -0.53 -10.09 12.79 5.56 0.07 -10.01 16 10 N -0.60 -9.15 2.91 -163.73 -0.48 -9.62 16 11 C 0.11 1.24 6.21 -3.45 -0.02 1.22 16 12 C -0.13 -1.59 5.92 -26.43 -0.16 -1.75 16 13 C -0.12 -1.96 5.45 -27.24 -0.15 -2.11 16 14 C 0.16 3.32 2.59 -68.36 -0.18 3.14 16 15 H 0.02 0.39 8.04 -51.93 -0.42 -0.03 16 16 C 0.12 2.21 3.34 -66.60 -0.22 1.99 16 17 H 0.08 1.53 8.11 -51.93 -0.42 1.11 16 18 C -0.13 -2.88 5.13 -25.16 -0.13 -3.01 16 19 C 0.17 4.41 5.48 -1.66 -0.01 4.40 16 20 N -0.59 -15.32 3.35 -163.25 -0.55 -15.86 16 21 C -0.23 -6.25 9.95 -15.06 -0.15 -6.40 16 22 H 0.08 2.06 7.32 -52.48 -0.38 1.68 16 23 C -0.01 -0.23 5.49 -17.43 -0.10 -0.33 16 24 N -0.92 -27.47 10.46 55.49 0.58 -26.89 16 25 Si 0.91 22.24 29.55 -169.99 -5.02 17.21 16 26 H -0.29 -7.02 7.11 56.52 0.40 -6.62 16 27 H -0.29 -7.21 7.11 56.52 0.40 -6.81 16 28 H -0.28 -6.70 7.11 56.52 0.40 -6.30 16 29 C -0.10 -1.00 5.36 -26.65 -0.14 -1.15 16 30 C 0.06 0.50 5.30 -3.83 -0.02 0.48 16 31 H 0.22 0.98 7.80 -54.09 -0.42 0.56 16 32 H 0.06 0.35 8.12 -51.92 -0.42 -0.07 16 33 H 0.10 0.60 7.99 -51.93 -0.41 0.18 16 34 H 0.04 0.26 8.14 -51.93 -0.42 -0.17 16 35 H 0.08 0.57 8.02 -51.93 -0.42 0.15 16 36 H 0.07 0.79 8.14 -51.93 -0.42 0.37 16 37 H 0.09 1.11 7.80 -51.93 -0.40 0.70 16 38 H 0.09 0.80 6.27 -51.93 -0.33 0.47 16 39 H 0.06 0.67 8.14 -51.93 -0.42 0.25 16 40 H 0.08 0.68 5.61 -51.93 -0.29 0.38 16 41 H 0.07 0.70 8.14 -51.93 -0.42 0.28 16 42 H 0.06 0.79 8.14 -51.93 -0.42 0.37 16 43 H 0.06 1.08 8.09 -51.93 -0.42 0.66 16 44 H 0.07 1.18 7.70 -51.93 -0.40 0.78 16 45 H 0.06 1.22 8.14 -51.93 -0.42 0.80 16 46 H 0.09 1.97 6.00 -51.93 -0.31 1.66 16 47 H 0.04 1.07 7.83 -51.93 -0.41 0.67 16 48 H 0.06 1.56 7.16 -51.93 -0.37 1.19 16 49 H 0.25 7.33 8.31 -40.82 -0.34 6.99 16 50 H 0.08 0.93 8.14 -51.92 -0.42 0.51 16 51 H 0.07 0.64 8.07 -51.93 -0.42 0.22 16 52 H 0.04 0.27 8.14 -51.93 -0.42 -0.15 16 53 H 0.08 0.68 7.11 -51.93 -0.37 0.31 16 LS Contribution 401.87 15.07 6.06 6.06 Total: -1.00 -33.12 401.87 -11.17 -44.29 By element: Atomic # 1 Polarization: 9.30 SS G_CDS: -9.13 Total: 0.17 kcal Atomic # 6 Polarization: 2.06 SS G_CDS: -1.60 Total: 0.47 kcal Atomic # 7 Polarization: -56.63 SS G_CDS: -1.55 Total: -58.18 kcal Atomic # 8 Polarization: -10.09 SS G_CDS: 0.07 Total: -10.01 kcal Atomic # 14 Polarization: 22.24 SS G_CDS: -5.02 Total: 17.21 kcal Total LS contribution 6.06 Total: 6.06 kcal Total: -33.12 -11.17 -44.29 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. ZINC001049917973.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 107.074 kcal (2) G-P(sol) polarization free energy of solvation -33.117 kcal (3) G-ENP(sol) elect.-nuc.-pol. free energy of system 73.957 kcal (4) G-CDS(sol) cavity-dispersion-solvent structure free energy -11.173 kcal (5) G-P-CDS(sol) = G-P(sol) + G-CDS(sol) = (2) + (4) -44.290 kcal (6) G-S(sol) free energy of system = (1) + (5) 62.784 kcal 1SCF WAS SPECIFIED, SO GEOMETRY OPTIMIZATION WAS NOT USED Total computer time = 1.21 seconds