Wall clock time and date at job start Sat Apr 11 2020 04:51:01 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.31622 * 1 3 3 Si 1.86803 * 119.99805 * 2 1 4 4 H 1.48503 * 109.47015 * 270.00279 * 3 2 1 5 5 H 1.48495 * 109.47164 * 30.00215 * 3 2 1 6 6 H 1.48503 * 109.46701 * 150.00279 * 3 2 1 7 7 C 1.38805 * 119.99650 * 180.02562 * 2 1 3 8 8 H 1.07995 * 120.00514 * 359.97438 * 7 2 1 9 9 N 1.36946 * 119.99766 * 179.97438 * 7 2 1 10 10 C 1.46500 * 120.00127 * 179.97438 * 9 7 2 11 11 C 1.53002 * 109.49223 * 60.07031 * 10 9 7 12 12 O 1.42898 * 109.47507 * 304.93428 * 11 10 9 13 13 C 1.53031 * 109.49852 * 300.00743 * 10 9 7 14 14 C 1.53038 * 109.53889 * 301.40902 * 13 10 9 15 15 C 1.53193 * 109.31495 * 298.63639 * 14 13 10 16 16 N 1.46930 * 108.77113 * 54.63457 * 15 14 13 17 17 C 1.34773 * 120.62999 * 126.37377 * 16 15 14 18 18 O 1.21280 * 119.99979 * 180.02562 * 17 16 15 19 19 C 1.50702 * 120.00021 * 359.97438 * 17 16 15 20 20 H 1.09001 * 109.47464 * 300.00888 * 19 17 16 21 21 C 1.52997 * 109.46960 * 60.00361 * 19 17 16 22 22 C 1.53006 * 109.47164 * 120.00088 * 21 19 17 23 23 C 1.52997 * 109.47127 * 59.99836 * 22 21 19 24 24 C 1.53000 * 109.47060 * 300.00449 * 23 22 21 25 25 C 1.52999 * 109.47043 * 180.02562 * 19 17 16 26 26 C 1.52998 * 109.47086 * 299.99851 * 25 19 17 27 27 C 1.52999 * 109.46623 * 299.99900 * 22 21 19 28 28 C 1.46923 * 118.73728 * 306.64097 * 16 15 14 29 29 H 0.97002 * 119.99722 * 359.97438 * 1 2 3 30 30 H 0.97002 * 119.99493 * 0.02562 * 9 7 2 31 31 H 1.08998 * 109.47100 * 64.93627 * 11 10 9 32 32 H 1.09000 * 109.46984 * 184.93356 * 11 10 9 33 33 H 0.96696 * 114.00708 * 180.02562 * 12 11 10 34 34 H 1.09006 * 109.46153 * 181.39115 * 13 10 9 35 35 H 1.09004 * 109.46375 * 61.43286 * 13 10 9 36 36 H 1.09000 * 109.49806 * 178.68187 * 14 13 10 37 37 H 1.08997 * 109.52608 * 58.60803 * 14 13 10 38 38 H 1.08993 * 109.58388 * 294.84994 * 15 14 13 39 39 H 1.08993 * 109.70346 * 174.48798 * 15 14 13 40 40 H 1.09001 * 109.47375 * 240.00122 * 21 19 17 41 41 H 1.09000 * 109.47140 * 0.02562 * 21 19 17 42 42 H 1.08996 * 109.46884 * 180.02562 * 22 21 19 43 43 H 1.09002 * 109.47153 * 180.02562 * 23 22 21 44 44 H 1.09001 * 109.46832 * 60.00321 * 23 22 21 45 45 H 1.09000 * 109.47517 * 120.00106 * 24 23 22 46 46 H 1.09004 * 109.46796 * 240.00014 * 24 23 22 47 47 H 1.08993 * 109.47655 * 59.99873 * 25 19 17 48 48 H 1.09002 * 109.47216 * 60.00424 * 26 25 19 49 49 H 1.09004 * 109.47344 * 180.02562 * 26 25 19 50 50 H 1.09002 * 109.46716 * 179.97438 * 27 22 21 51 51 H 1.09000 * 109.47486 * 299.99829 * 27 22 21 52 52 H 1.09004 * 109.58683 * 173.14937 * 28 16 15 53 53 H 1.09003 * 109.59183 * 293.57031 * 28 16 15 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.6178 0.0000 4 1 2.4976 2.0465 -1.4001 5 1 1.4476 2.6527 0.7001 6 1 3.5478 1.4401 0.7000 7 6 2.0102 -1.2021 -0.0005 8 1 1.4702 -2.1374 -0.0005 9 7 3.3796 -1.2022 -0.0011 10 6 4.1121 -2.4709 -0.0011 11 6 3.7437 -3.2731 -1.2507 12 8 2.3253 -3.4358 -1.3114 13 6 3.7432 -3.2737 1.2484 14 6 4.0804 -2.4582 2.4988 15 6 5.5863 -2.1788 2.5307 16 7 5.9748 -1.5359 1.2680 17 6 6.6450 -0.3666 1.2706 18 8 6.9596 0.1543 0.2215 19 6 7.0025 0.2941 2.5770 20 1 6.0920 0.5112 3.1355 21 6 7.8918 -0.6435 3.3960 22 6 9.2365 0.0351 3.6655 23 6 9.0066 1.3309 4.4457 24 6 8.1172 2.2684 3.6267 25 6 7.7537 1.5979 2.3003 26 6 9.0318 1.2936 1.5161 27 6 9.9216 0.3553 2.3354 28 6 5.6134 -2.1888 0.0024 29 1 -0.4850 0.8401 0.0004 30 1 3.8646 -0.3621 -0.0007 31 1 4.0867 -2.7412 -2.1381 32 1 4.2202 -4.2525 -1.2073 33 1 2.0193 -3.9366 -2.0799 34 1 4.3071 -4.2065 1.2604 35 1 2.6760 -3.4949 1.2360 36 1 3.7958 -3.0209 3.3879 37 1 3.5357 -1.5143 2.4760 38 1 6.1295 -3.1167 2.6454 39 1 5.8196 -1.5170 3.3646 40 1 7.4037 -0.8717 4.3436 41 1 8.0560 -1.5665 2.8399 42 1 9.8702 -0.6331 4.2486 43 1 9.9645 1.8145 4.6373 44 1 8.5185 1.1027 5.3932 45 1 7.2066 2.4856 4.1851 46 1 8.6527 3.1971 3.4295 47 1 7.1199 2.2657 1.7169 48 1 8.7732 0.8158 0.5712 49 1 9.5672 2.2225 1.3194 50 1 10.8793 0.8391 2.5273 51 1 10.0854 -0.5679 1.7795 52 1 5.8625 -1.5325 -0.8315 53 1 6.1590 -3.1276 -0.0933 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 ZINC001327608564.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 04:51:01 Heat of formation + Delta-G solvation = -75.000386 kcal Electronic energy + Delta-G solvation = -32682.413029 eV Core-core repulsion = 28746.654561 eV Total energy + Delta-G solvation = -3935.758468 eV No. of doubly occupied orbitals = 65 Molecular weight (most abundant/longest-lived isotopes) = 336.220 amu Computer time = 0.81 seconds Orbital eigenvalues (eV) -41.29311 -39.38034 -37.63185 -35.80134 -34.00302 -32.27952 -31.65487 -31.35903 -30.99358 -30.81922 -29.20444 -25.52639 -25.45937 -23.39426 -22.97753 -22.59822 -22.06533 -21.70608 -19.34830 -18.36561 -17.75456 -17.17501 -16.73931 -16.58613 -16.35750 -15.87729 -15.53742 -14.99620 -14.57660 -14.32391 -14.07654 -14.01159 -13.48620 -13.26450 -12.98014 -12.85381 -12.66218 -12.54868 -12.24917 -12.17820 -12.14627 -12.10975 -12.04990 -12.02091 -11.87110 -11.63198 -11.21217 -10.85633 -10.83874 -10.72994 -10.40555 -10.36211 -10.24457 -10.16628 -9.97712 -9.94196 -9.89570 -9.39630 -9.16926 -8.94332 -8.75038 -8.29799 -7.13272 -5.57090 -2.91106 2.86266 3.18023 3.47124 3.52760 3.56808 4.35312 4.43210 4.53098 4.72287 4.79741 4.91412 4.99410 5.10950 5.13262 5.16874 5.19782 5.20917 5.25103 5.31486 5.35784 5.37370 5.45158 5.47686 5.56511 5.60167 5.63533 5.69076 5.76585 5.83003 5.85726 5.88383 5.93806 5.97739 5.98584 6.04835 6.08313 6.10932 6.16218 6.23068 6.29037 6.31612 6.35922 6.42086 6.54412 6.66645 6.79686 6.92254 7.04463 7.65477 8.27979 8.40472 8.44517 8.65429 9.43227 10.11800 10.29092 11.61274 Molecular weight = 336.22amu Principal moments of inertia in cm(-1) A = 0.014134 B = 0.004511 C = 0.003939 Principal moments of inertia in units of 10**(-40)*gram-cm**2 A = 1980.579581 B = 6205.291754 C = 7106.338210 Net atomic charges, atomic populations, and dipole contributions calculated with CM2 Atom NO. Type Charge No. of electrons 1 N -0.939 5.939 2 C -0.003 4.003 3 Si 0.869 3.131 4 H -0.289 1.289 5 H -0.295 1.295 6 H -0.281 1.281 7 C -0.208 4.208 8 H 0.116 0.884 9 N -0.731 5.731 10 C 0.195 3.805 11 C 0.075 3.925 12 O -0.529 6.529 13 C -0.116 4.116 14 C -0.130 4.130 15 C 0.106 3.894 16 N -0.599 5.599 17 C 0.523 3.477 18 O -0.540 6.540 19 C -0.101 4.101 20 H 0.095 0.905 21 C -0.105 4.105 22 C -0.086 4.086 23 C -0.112 4.112 24 C -0.115 4.115 25 C -0.065 4.065 26 C -0.125 4.125 27 C -0.112 4.112 28 C 0.112 3.888 29 H 0.254 0.746 30 H 0.360 0.640 31 H 0.035 0.965 32 H 0.025 0.975 33 H 0.352 0.648 34 H 0.043 0.957 35 H 0.079 0.921 36 H 0.062 0.938 37 H 0.085 0.915 38 H 0.062 0.938 39 H 0.082 0.918 40 H 0.071 0.929 41 H 0.068 0.932 42 H 0.070 0.930 43 H 0.060 0.940 44 H 0.061 0.939 45 H 0.062 0.938 46 H 0.063 0.937 47 H 0.075 0.925 48 H 0.092 0.908 49 H 0.055 0.945 50 H 0.057 0.943 51 H 0.057 0.943 52 H 0.092 0.908 53 H 0.067 0.933 For the charges calculated by CM2: Dipole moment (debyes) X Y Z Total from point charges 19.595 -4.758 8.323 21.815 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.589 5.589 2 C -0.196 4.196 3 Si 0.701 3.299 4 H -0.216 1.216 5 H -0.222 1.222 6 H -0.207 1.207 7 C -0.337 4.337 8 H 0.133 0.867 9 N -0.373 5.373 10 C 0.102 3.898 11 C -0.006 4.006 12 O -0.331 6.331 13 C -0.156 4.156 14 C -0.168 4.168 15 C -0.017 4.017 16 N -0.331 5.331 17 C 0.311 3.689 18 O -0.418 6.418 19 C -0.123 4.123 20 H 0.114 0.886 21 C -0.142 4.142 22 C -0.104 4.104 23 C -0.150 4.150 24 C -0.152 4.152 25 C -0.084 4.084 26 C -0.162 4.162 27 C -0.150 4.150 28 C -0.011 4.011 29 H 0.064 0.936 30 H 0.193 0.807 31 H 0.053 0.947 32 H 0.043 0.957 33 H 0.197 0.803 34 H 0.062 0.938 35 H 0.098 0.902 36 H 0.081 0.919 37 H 0.103 0.897 38 H 0.081 0.919 39 H 0.100 0.900 40 H 0.090 0.910 41 H 0.087 0.913 42 H 0.088 0.912 43 H 0.079 0.921 44 H 0.080 0.920 45 H 0.080 0.920 46 H 0.082 0.918 47 H 0.094 0.906 48 H 0.110 0.890 49 H 0.074 0.926 50 H 0.076 0.924 51 H 0.076 0.924 52 H 0.110 0.890 53 H 0.085 0.915 Dipole moment (debyes) X Y Z Total from point charges 19.566 -4.175 8.332 21.672 hybrid contribution 0.394 0.952 -0.552 1.169 sum 19.960 -3.222 7.780 21.663 Atomic orbital electron populations 1.69684 1.04829 1.25012 1.59387 1.24617 0.95074 0.96361 1.03586 0.86778 0.79109 0.86317 0.77714 1.21631 1.22205 1.20713 1.19408 0.80114 0.95082 1.39139 0.86654 1.42805 1.04304 1.00842 1.89327 1.19828 0.89116 0.87977 0.92885 1.22606 0.83067 0.96386 0.98565 1.86996 1.20234 1.77569 1.48293 1.21840 1.02798 0.96279 0.94642 1.21822 0.94218 1.02714 0.98045 1.21927 0.96719 0.97401 0.85612 1.47896 1.55479 1.22492 1.07273 1.20568 0.77092 0.79997 0.91208 1.90596 1.54140 1.63223 1.33841 1.21305 0.99997 0.96451 0.94517 0.88644 1.21683 0.95187 0.97807 0.99557 1.21042 0.96124 0.96652 0.96581 1.21542 0.99697 0.95601 0.98126 1.21567 0.99627 0.98223 0.95804 1.20775 0.96817 0.94712 0.96065 1.21787 0.94756 0.99386 1.00307 1.21563 0.98223 0.99635 0.95580 1.21634 0.94942 0.97577 0.86911 0.93581 0.80681 0.94703 0.95667 0.80312 0.93810 0.90222 0.91919 0.89673 0.91948 0.89987 0.91022 0.91291 0.91186 0.92140 0.92043 0.91971 0.91815 0.90646 0.89006 0.92642 0.92404 0.92436 0.88984 0.91470 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 N -0.94 -30.05 13.97 55.49 0.78 -29.27 16 2 C 0.00 -0.08 5.50 -17.43 -0.10 -0.18 16 3 Si 0.87 23.32 27.74 -169.99 -4.71 18.61 16 4 H -0.29 -7.93 7.11 56.52 0.40 -7.53 16 5 H -0.30 -8.02 7.11 56.52 0.40 -7.62 16 6 H -0.28 -7.29 6.52 56.52 0.37 -6.92 16 7 C -0.21 -5.74 8.62 -15.06 -0.13 -5.87 16 8 H 0.12 3.44 5.52 -52.49 -0.29 3.15 16 9 N -0.73 -17.05 4.02 -47.33 -0.19 -17.24 16 10 C 0.19 3.57 0.68 -131.93 -0.09 3.48 16 11 C 0.07 1.29 6.09 37.16 0.23 1.52 16 12 O -0.53 -10.35 9.92 -35.23 -0.35 -10.70 16 13 C -0.12 -1.96 4.42 -26.70 -0.12 -2.08 16 14 C -0.13 -2.01 5.93 -26.61 -0.16 -2.17 16 15 C 0.11 1.35 5.18 -3.71 -0.02 1.33 16 16 N -0.60 -9.14 2.81 -173.04 -0.49 -9.63 16 17 C 0.52 8.49 6.55 -10.99 -0.07 8.42 16 18 O -0.54 -10.42 13.07 5.56 0.07 -10.35 16 19 C -0.10 -1.32 1.76 -91.77 -0.16 -1.48 16 20 H 0.10 1.26 7.70 -51.93 -0.40 0.86 16 21 C -0.10 -1.05 3.73 -26.73 -0.10 -1.15 16 22 C -0.09 -0.79 3.20 -90.62 -0.29 -1.08 16 23 C -0.11 -0.98 5.71 -26.73 -0.15 -1.13 16 24 C -0.11 -1.18 5.71 -26.73 -0.15 -1.33 16 25 C -0.06 -0.83 2.36 -90.62 -0.21 -1.05 16 26 C -0.13 -1.56 5.39 -26.69 -0.14 -1.70 16 27 C -0.11 -1.16 5.71 -26.69 -0.15 -1.31 16 28 C 0.11 1.80 5.10 -3.94 -0.02 1.78 16 29 H 0.25 7.86 8.31 -40.82 -0.34 7.52 16 30 H 0.36 8.90 5.50 -40.82 -0.22 8.67 16 31 H 0.03 0.61 8.14 -51.93 -0.42 0.18 16 32 H 0.03 0.36 8.14 -51.93 -0.42 -0.06 16 33 H 0.35 5.23 9.12 45.56 0.42 5.65 16 34 H 0.04 0.62 8.14 -51.93 -0.42 0.20 16 35 H 0.08 1.60 6.97 -51.93 -0.36 1.24 16 36 H 0.06 0.83 8.14 -51.93 -0.42 0.41 16 37 H 0.08 1.65 7.53 -51.93 -0.39 1.26 16 38 H 0.06 0.66 8.14 -51.93 -0.42 0.24 16 39 H 0.08 0.91 4.58 -51.93 -0.24 0.67 16 40 H 0.07 0.61 7.51 -51.93 -0.39 0.22 16 41 H 0.07 0.71 6.84 -51.93 -0.36 0.36 16 42 H 0.07 0.56 8.14 -51.93 -0.42 0.14 16 43 H 0.06 0.48 8.14 -51.93 -0.42 0.06 16 44 H 0.06 0.49 8.14 -51.93 -0.42 0.07 16 45 H 0.06 0.63 8.14 -51.93 -0.42 0.20 16 46 H 0.06 0.62 8.14 -51.93 -0.42 0.19 16 47 H 0.08 1.14 8.10 -51.93 -0.42 0.72 16 48 H 0.09 1.40 5.51 -51.93 -0.29 1.11 16 49 H 0.05 0.65 8.14 -51.93 -0.42 0.22 16 50 H 0.06 0.53 8.14 -51.93 -0.42 0.11 16 51 H 0.06 0.61 8.14 -51.93 -0.42 0.19 16 52 H 0.09 1.65 7.01 -51.93 -0.36 1.29 16 53 H 0.07 0.89 8.14 -51.93 -0.42 0.47 16 LS Contribution 378.10 15.07 5.70 5.70 Total: -1.00 -34.16 378.10 -9.43 -43.59 By element: Atomic # 1 Polarization: 21.67 SS G_CDS: -8.39 Total: 13.27 kcal Atomic # 6 Polarization: -2.14 SS G_CDS: -1.84 Total: -3.98 kcal Atomic # 7 Polarization: -56.24 SS G_CDS: 0.10 Total: -56.14 kcal Atomic # 8 Polarization: -20.77 SS G_CDS: -0.28 Total: -21.04 kcal Atomic # 14 Polarization: 23.32 SS G_CDS: -4.71 Total: 18.61 kcal Total LS contribution 5.70 Total: 5.70 kcal Total: -34.16 -9.43 -43.59 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. ZINC001327608564.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 -31.413 kcal (2) G-P(sol) polarization free energy of solvation -34.157 kcal (3) G-ENP(sol) elect.-nuc.-pol. free energy of system -65.570 kcal (4) G-CDS(sol) cavity-dispersion-solvent structure free energy -9.430 kcal (5) G-P-CDS(sol) = G-P(sol) + G-CDS(sol) = (2) + (4) -43.588 kcal (6) G-S(sol) free energy of system = (1) + (5) -75.000 kcal 1SCF WAS SPECIFIED, SO GEOMETRY OPTIMIZATION WAS NOT USED Total computer time = 0.82 seconds