Wall clock time and date at job start Sat Apr 11 2020 06:11:28 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.52998 * 1 3 3 O 1.42900 * 109.47212 * 2 1 4 4 C 1.42899 * 114.00245 * 180.02562 * 3 2 1 5 5 C 1.53780 * 113.61494 * 209.99718 * 4 3 2 6 6 C 1.53778 * 87.08457 * 220.01885 * 5 4 3 7 7 C 1.53003 * 113.61616 * 139.97916 * 6 5 4 8 8 C 1.50701 * 109.46871 * 176.53074 * 7 6 5 9 9 O 1.21277 * 119.99936 * 0.02562 * 8 7 6 10 10 N 1.34775 * 119.99403 * 179.97438 * 8 7 6 11 11 C 1.46932 * 120.62573 * 179.97438 * 10 8 7 12 12 C 1.53189 * 108.77592 * 233.63105 * 11 10 8 13 13 C 1.53043 * 109.31439 * 305.37036 * 12 11 10 14 14 O 1.42898 * 109.46024 * 301.35607 * 13 12 11 15 15 C 1.53006 * 109.45814 * 181.38435 * 13 12 11 16 16 N 1.46500 * 109.46998 * 184.99823 * 15 13 12 17 17 C 1.36931 * 120.00396 * 180.02562 * 16 15 13 18 18 H 1.08003 * 119.99774 * 359.97438 * 17 16 15 19 19 C 1.38806 * 120.00545 * 180.02562 * 17 16 15 20 20 N 1.31624 * 120.00512 * 179.97438 * 19 17 16 21 21 Si 1.86799 * 119.99852 * 359.97438 * 19 17 16 22 22 H 1.48499 * 109.47136 * 90.00264 * 21 19 17 23 23 H 1.48499 * 109.47548 * 210.00043 * 21 19 17 24 24 H 1.48496 * 109.47428 * 330.00419 * 21 19 17 25 25 C 1.53037 * 109.54135 * 61.36523 * 13 12 11 26 26 C 1.46924 * 120.63345 * 0.02585 * 10 8 7 27 27 C 1.53788 * 87.16965 * 25.38925 * 6 5 4 28 28 H 1.09001 * 109.46812 * 180.02562 * 1 2 3 29 29 H 1.09000 * 109.47025 * 300.00027 * 1 2 3 30 30 H 1.08997 * 109.47078 * 60.00208 * 1 2 3 31 31 H 1.09000 * 109.47025 * 239.99973 * 2 1 3 32 32 H 1.08997 * 109.47078 * 119.99792 * 2 1 3 33 33 H 1.08992 * 112.90685 * 341.25126 * 4 3 2 34 34 H 1.08997 * 113.61124 * 105.46926 * 5 4 3 35 35 H 1.09006 * 113.69328 * 334.58095 * 5 4 3 36 36 H 1.08997 * 113.61450 * 270.80467 * 6 5 4 37 37 H 1.09003 * 109.46594 * 56.53471 * 7 6 5 38 38 H 1.08996 * 109.47053 * 296.53150 * 7 6 5 39 39 H 1.08997 * 109.58232 * 113.84551 * 11 10 8 40 40 H 1.08992 * 109.58821 * 353.42548 * 11 10 8 41 41 H 1.08995 * 109.50151 * 185.41691 * 12 11 10 42 42 H 1.09003 * 109.49700 * 65.31889 * 12 11 10 43 43 H 0.96698 * 114.00272 * 180.02562 * 14 13 12 44 44 H 1.08993 * 109.46955 * 305.00048 * 15 13 12 45 45 H 1.09000 * 109.47050 * 64.99877 * 15 13 12 46 46 H 0.97010 * 119.99909 * 359.97438 * 16 15 13 47 47 H 0.96997 * 120.00065 * 179.97438 * 20 19 17 48 48 H 1.08998 * 109.49633 * 58.58307 * 25 13 12 49 49 H 1.08991 * 109.50061 * 178.68396 * 25 13 12 50 50 H 1.08998 * 109.58360 * 246.14876 * 26 10 8 51 51 H 1.08999 * 109.58623 * 6.57990 * 26 10 8 52 52 H 1.09000 * 113.60890 * 89.16316 * 27 6 5 53 53 H 1.08995 * 113.61378 * 220.05704 * 27 6 5 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 8 2.0063 1.3473 0.0000 4 6 3.4309 1.4601 -0.0006 5 6 3.9489 2.7252 0.7038 6 6 4.9896 2.8315 -0.4233 7 6 5.2307 4.2647 -0.9017 8 6 6.3221 4.2705 -1.9409 9 8 6.8542 3.2306 -2.2670 10 7 6.7059 5.4311 -2.5084 11 6 7.7688 5.4528 -3.5226 12 6 8.8401 6.4599 -3.0926 13 6 8.1843 7.8195 -2.8400 14 8 7.5273 8.2596 -4.0302 15 6 9.2564 8.8367 -2.4437 16 7 8.6173 10.1056 -2.0864 17 6 9.3828 11.1719 -1.6964 18 1 10.4577 11.0773 -1.6517 19 6 8.7772 12.3740 -1.3574 20 7 9.5130 13.3988 -0.9820 21 14 6.9180 12.5374 -1.4340 22 1 6.3334 12.1589 -0.1224 23 1 6.5571 13.9429 -1.7496 24 1 6.3852 11.6404 -2.4907 25 6 7.1602 7.6918 -1.7100 26 6 6.0702 6.6986 -2.1240 27 6 4.0038 2.0456 -1.3040 28 1 -0.3633 -1.0277 0.0005 29 1 -0.3633 0.5138 0.8900 30 1 -0.3633 0.5138 -0.8900 31 1 1.8933 -0.5138 -0.8900 32 1 1.8933 -0.5138 0.8900 33 1 3.9290 0.5459 0.3218 34 1 4.3849 2.5316 1.6839 35 1 3.2312 3.5456 0.7144 36 1 5.9144 2.2933 -0.2158 37 1 5.5305 4.8842 -0.0564 38 1 4.3133 4.6619 -1.3360 39 1 7.3507 5.7520 -4.4837 40 1 8.2128 4.4611 -3.6079 41 1 9.5860 6.5567 -3.8814 42 1 9.3202 6.1121 -2.1779 43 1 7.0888 9.1171 -3.9433 44 1 9.8172 8.4589 -1.5888 45 1 9.9348 8.9949 -3.2821 46 1 7.6518 10.1907 -2.1269 47 1 9.0898 14.2389 -0.7455 48 1 7.6561 7.3323 -0.8084 49 1 6.7104 8.6652 -1.5149 50 1 5.5134 7.0997 -2.9709 51 1 5.3935 6.5288 -1.2865 52 1 3.3005 2.6792 -1.8444 53 1 4.4842 1.3044 -1.9426 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 ZINC001570419507.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:11:28 Heat of formation + Delta-G solvation = -102.908622 kcal Electronic energy + Delta-G solvation = -30430.726855 eV Core-core repulsion = 26303.302143 eV Total energy + Delta-G solvation = -4127.424711 eV No. of doubly occupied orbitals = 66 Molecular weight (most abundant/longest-lived isotopes) = 340.220 amu Computer time = 0.72 seconds Orbital eigenvalues (eV) -40.65531 -39.69864 -38.52270 -36.76388 -35.59572 -34.22589 -32.35633 -31.73888 -31.34404 -30.65548 -29.76676 -26.99474 -25.99801 -24.24813 -23.49669 -23.31701 -22.85870 -21.54000 -20.39503 -19.78673 -19.01235 -18.06852 -17.24691 -17.21984 -16.19112 -15.99862 -15.86288 -15.76450 -15.43451 -14.96573 -14.59688 -14.52104 -14.29757 -14.04163 -13.65879 -13.39192 -13.27327 -13.11943 -13.06777 -12.72191 -12.67333 -12.61821 -12.45366 -12.29181 -12.05726 -11.75467 -11.67043 -11.57084 -11.44244 -11.22686 -11.19955 -11.08998 -10.96301 -10.93959 -10.60372 -10.34817 -10.26713 -9.81234 -9.71539 -9.66730 -9.56450 -9.01603 -8.63676 -7.35801 -5.85413 -3.22719 2.56886 2.90689 3.23959 3.28174 3.35591 3.36822 4.13049 4.27660 4.36906 4.43202 4.48474 4.49319 4.60560 4.66764 4.68726 4.76539 4.84329 4.88873 4.96943 4.97756 5.04390 5.16680 5.24138 5.27420 5.33105 5.40733 5.41966 5.44382 5.48956 5.49808 5.50860 5.52026 5.58367 5.70090 5.71040 5.75928 5.84166 5.87481 5.89702 5.99171 6.11504 6.13782 6.22335 6.47284 6.52132 6.57519 6.94874 7.35234 7.83256 7.91408 8.11069 8.32103 9.09093 9.80936 9.98771 11.32807 Molecular weight = 340.22amu Principal moments of inertia in cm(-1) A = 0.016937 B = 0.002127 C = 0.001971 Principal moments of inertia in units of 10**(-40)*gram-cm**2 A = 1652.757060 B =13162.929520 C =14199.966614 Net atomic charges, atomic populations, and dipole contributions calculated with CM2 Atom NO. Type Charge No. of electrons 1 C -0.152 4.152 2 C 0.055 3.945 3 O -0.376 6.376 4 C 0.061 3.939 5 C -0.133 4.133 6 C -0.092 4.092 7 C -0.129 4.129 8 C 0.515 3.485 9 O -0.547 6.547 10 N -0.609 5.609 11 C 0.118 3.882 12 C -0.116 4.116 13 C 0.126 3.874 14 O -0.559 6.559 15 C 0.133 3.867 16 N -0.753 5.753 17 C -0.247 4.247 18 H 0.091 0.909 19 C 0.006 3.994 20 N -0.931 5.931 21 Si 0.885 3.115 22 H -0.284 1.284 23 H -0.290 1.290 24 H -0.279 1.279 25 C -0.160 4.160 26 C 0.110 3.890 27 C -0.146 4.146 28 H 0.069 0.931 29 H 0.062 0.938 30 H 0.062 0.938 31 H 0.047 0.953 32 H 0.047 0.953 33 H 0.086 0.914 34 H 0.073 0.927 35 H 0.092 0.908 36 H 0.102 0.898 37 H 0.097 0.903 38 H 0.094 0.906 39 H 0.074 0.926 40 H 0.089 0.911 41 H 0.083 0.917 42 H 0.079 0.921 43 H 0.378 0.622 44 H 0.035 0.965 45 H 0.038 0.962 46 H 0.357 0.643 47 H 0.260 0.740 48 H 0.079 0.921 49 H 0.083 0.917 50 H 0.075 0.925 51 H 0.083 0.917 52 H 0.087 0.913 53 H 0.071 0.929 For the charges calculated by CM2: Dipole moment (debyes) X Y Z Total from point charges -14.798 -25.460 0.941 29.464 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.209 4.209 2 C -0.021 4.021 3 O -0.295 6.295 4 C 0.004 3.996 5 C -0.170 4.170 6 C -0.111 4.111 7 C -0.169 4.169 8 C 0.304 3.696 9 O -0.425 6.425 10 N -0.342 5.342 11 C -0.005 4.005 12 C -0.154 4.154 13 C 0.085 3.915 14 O -0.365 6.365 15 C 0.005 3.995 16 N -0.400 5.400 17 C -0.376 4.376 18 H 0.108 0.892 19 C -0.190 4.190 20 N -0.579 5.579 21 Si 0.715 3.285 22 H -0.210 1.210 23 H -0.216 1.216 24 H -0.205 1.205 25 C -0.198 4.198 26 C -0.013 4.013 27 C -0.183 4.183 28 H 0.088 0.912 29 H 0.081 0.919 30 H 0.082 0.918 31 H 0.066 0.934 32 H 0.065 0.935 33 H 0.104 0.896 34 H 0.092 0.908 35 H 0.110 0.890 36 H 0.120 0.880 37 H 0.115 0.885 38 H 0.112 0.888 39 H 0.092 0.908 40 H 0.107 0.893 41 H 0.102 0.898 42 H 0.097 0.903 43 H 0.226 0.774 44 H 0.053 0.947 45 H 0.056 0.944 46 H 0.191 0.809 47 H 0.071 0.929 48 H 0.097 0.903 49 H 0.101 0.899 50 H 0.093 0.907 51 H 0.101 0.899 52 H 0.105 0.895 53 H 0.090 0.910 Dipole moment (debyes) X Y Z Total from point charges -14.205 -26.275 0.597 29.875 hybrid contribution -1.273 0.893 0.272 1.579 sum -15.478 -25.381 0.869 29.741 Atomic orbital electron populations 1.21751 0.93851 1.02191 1.03140 1.22498 0.95567 0.84138 0.99934 1.87952 1.20808 1.26139 1.94618 1.23705 0.80834 0.97902 0.97205 1.23312 0.97717 0.97977 0.98003 1.22068 0.97930 0.94705 0.96349 1.21679 0.98303 0.98623 0.98321 1.20650 0.83785 0.83017 0.82140 1.90660 1.56786 1.31486 1.63522 1.48160 1.38020 1.06917 1.41101 1.21510 0.87882 1.00173 0.90891 1.21838 0.97855 0.92801 1.02934 1.21381 0.91850 0.93983 0.84321 1.86414 1.73111 1.41228 1.35773 1.21308 0.94027 0.85922 0.98199 1.42687 1.04723 1.09692 1.82909 1.19032 0.92361 0.88926 1.37321 0.89166 1.24690 0.97224 0.95862 1.01272 1.69695 1.34044 1.02667 1.51470 0.86564 0.87088 0.77195 0.77623 1.21015 1.21579 1.20545 1.22556 0.97633 0.99967 0.99670 1.21736 0.95288 0.83553 1.00677 1.23554 0.97774 1.00135 0.96882 0.91205 0.91877 0.91842 0.93433 0.93493 0.89630 0.90828 0.88963 0.88017 0.88519 0.88751 0.90761 0.89306 0.89812 0.90282 0.77359 0.94727 0.94430 0.80939 0.92948 0.90267 0.89893 0.90701 0.89886 0.89468 0.91027 Number of geometries 1 Number of calculations of the screened coulomb radii 1 The total number of SCF iterations 14. 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 -0.88 10.28 37.16 0.38 -0.50 16 2 C 0.06 0.39 6.50 37.16 0.24 0.63 16 3 O -0.38 -3.53 10.36 -35.23 -0.36 -3.89 16 4 C 0.06 0.51 4.31 -25.83 -0.11 0.39 16 5 C -0.13 -1.02 7.47 -25.92 -0.19 -1.21 16 6 C -0.09 -0.81 3.11 -89.80 -0.28 -1.09 16 7 C -0.13 -1.15 3.87 -27.88 -0.11 -1.26 16 8 C 0.52 6.28 7.70 -10.99 -0.08 6.20 16 9 O -0.55 -8.32 15.20 5.56 0.08 -8.23 16 10 N -0.61 -7.21 2.97 -172.79 -0.51 -7.73 16 11 C 0.12 1.46 6.19 -3.71 -0.02 1.44 16 12 C -0.12 -1.47 5.11 -26.61 -0.14 -1.61 16 13 C 0.13 1.81 0.71 -90.58 -0.06 1.74 16 14 O -0.56 -8.27 11.70 -35.23 -0.41 -8.68 16 15 C 0.13 2.43 5.43 -4.04 -0.02 2.40 16 16 N -0.75 -16.53 4.87 -59.91 -0.29 -16.83 16 17 C -0.25 -6.73 10.47 -15.06 -0.16 -6.89 16 18 H 0.09 2.68 8.06 -52.48 -0.42 2.25 16 19 C 0.01 0.16 5.50 -17.43 -0.10 0.06 16 20 N -0.93 -28.57 13.97 55.49 0.78 -27.79 16 21 Si 0.88 20.91 27.74 -169.99 -4.71 16.20 16 22 H -0.28 -6.63 7.11 56.52 0.40 -6.23 16 23 H -0.29 -7.04 7.11 56.52 0.40 -6.64 16 24 H -0.28 -5.96 6.52 56.52 0.37 -5.59 16 25 C -0.16 -2.07 4.87 -26.63 -0.13 -2.20 16 26 C 0.11 1.18 6.10 -3.73 -0.02 1.16 16 27 C -0.15 -1.32 7.48 -25.83 -0.19 -1.51 16 28 H 0.07 0.29 8.14 -51.93 -0.42 -0.14 16 29 H 0.06 0.38 8.14 -51.93 -0.42 -0.04 16 30 H 0.06 0.39 8.14 -51.93 -0.42 -0.03 16 31 H 0.05 0.31 8.14 -51.93 -0.42 -0.12 16 32 H 0.05 0.28 8.14 -51.93 -0.42 -0.14 16 33 H 0.09 0.61 7.87 -51.93 -0.41 0.20 16 34 H 0.07 0.50 8.14 -51.93 -0.42 0.08 16 35 H 0.09 0.71 8.10 -51.93 -0.42 0.29 16 36 H 0.10 1.00 7.73 -51.93 -0.40 0.60 16 37 H 0.10 0.74 7.69 -51.93 -0.40 0.34 16 38 H 0.09 0.73 7.67 -51.93 -0.40 0.33 16 39 H 0.07 0.95 7.81 -51.93 -0.41 0.54 16 40 H 0.09 1.14 7.00 -51.93 -0.36 0.78 16 41 H 0.08 1.03 8.14 -51.93 -0.42 0.61 16 42 H 0.08 0.98 8.14 -51.93 -0.42 0.56 16 43 H 0.38 5.39 8.19 45.56 0.37 5.77 16 44 H 0.03 0.64 8.14 -51.93 -0.42 0.22 16 45 H 0.04 0.73 8.14 -51.93 -0.42 0.31 16 46 H 0.36 7.38 3.41 -40.82 -0.14 7.24 16 47 H 0.26 7.60 8.31 -40.82 -0.34 7.26 16 48 H 0.08 1.04 8.14 -51.93 -0.42 0.61 16 49 H 0.08 1.16 6.53 -51.93 -0.34 0.82 16 50 H 0.07 0.80 7.81 -51.93 -0.41 0.40 16 51 H 0.08 0.70 5.93 -51.93 -0.31 0.39 16 52 H 0.09 0.82 7.92 -51.93 -0.41 0.40 16 53 H 0.07 0.69 8.14 -51.93 -0.42 0.27 16 LS Contribution 410.40 15.07 6.18 6.18 Total: -1.00 -32.70 410.40 -8.94 -41.65 By element: Atomic # 1 Polarization: 20.06 SS G_CDS: -8.69 Total: 11.37 kcal Atomic # 6 Polarization: -1.24 SS G_CDS: -1.00 Total: -2.24 kcal Atomic # 7 Polarization: -52.32 SS G_CDS: -0.03 Total: -52.35 kcal Atomic # 8 Polarization: -20.12 SS G_CDS: -0.69 Total: -20.81 kcal Atomic # 14 Polarization: 20.91 SS G_CDS: -4.71 Total: 16.20 kcal Total LS contribution 6.18 Total: 6.18 kcal Total: -32.70 -8.94 -41.65 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. ZINC001570419507.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 -61.261 kcal (2) G-P(sol) polarization free energy of solvation -32.704 kcal (3) G-ENP(sol) elect.-nuc.-pol. free energy of system -93.965 kcal (4) G-CDS(sol) cavity-dispersion-solvent structure free energy -8.944 kcal (5) G-P-CDS(sol) = G-P(sol) + G-CDS(sol) = (2) + (4) -41.648 kcal (6) G-S(sol) free energy of system = (1) + (5) -102.909 kcal 1SCF WAS SPECIFIED, SO GEOMETRY OPTIMIZATION WAS NOT USED Total computer time = 0.72 seconds