Wall clock time and date at job start Sat Apr 11 2020 06:09:19 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.31515 * 1 3 3 Si 1.86796 * 120.00243 * 2 1 4 4 H 1.48504 * 109.46897 * 240.00077 * 3 2 1 5 5 H 1.48497 * 109.47213 * 359.97438 * 3 2 1 6 6 H 1.48501 * 109.47113 * 119.99919 * 3 2 1 7 7 C 1.37881 * 119.99742 * 179.97438 * 2 1 3 8 8 H 1.08001 * 119.99697 * 179.97438 * 7 2 1 9 9 C 1.50703 * 119.99984 * 0.02562 * 7 2 1 10 10 O 1.42896 * 109.46765 * 180.02562 * 9 7 2 11 11 C 1.42909 * 113.99925 * 180.02562 * 10 9 7 12 12 C 1.52992 * 109.47400 * 180.02562 * 11 10 9 13 13 N 1.46901 * 109.47088 * 64.99470 * 12 11 10 14 14 C 1.46904 * 111.00134 * 175.04818 * 13 12 11 15 15 C 1.50697 * 109.47281 * 305.08350 * 14 13 12 16 16 O 1.21283 * 120.00383 * 122.31056 * 15 14 13 17 17 N 1.34779 * 119.99941 * 302.03386 * 15 14 13 18 18 C 1.47154 * 120.97855 * 354.16302 * 17 15 14 19 19 C 1.53712 * 108.16177 * 129.43297 * 18 17 15 20 20 O 1.42963 * 109.15095 * 52.53591 * 19 18 17 21 21 C 1.42966 * 114.16368 * 297.29217 * 20 19 18 22 22 C 1.47147 * 120.98654 * 174.44315 * 17 15 14 23 23 C 1.53002 * 109.46752 * 65.07919 * 14 13 12 24 24 C 1.53004 * 109.47083 * 180.02562 * 23 14 13 25 25 C 1.53004 * 109.46763 * 299.99781 * 24 23 14 26 26 C 1.53002 * 109.47142 * 60.00255 * 25 24 23 27 27 C 1.52991 * 109.47513 * 185.07704 * 14 13 12 28 28 H 0.96999 * 119.99878 * 0.02562 * 1 2 3 29 29 H 1.09003 * 109.46825 * 60.00076 * 9 7 2 30 30 H 1.08999 * 109.47178 * 300.00102 * 9 7 2 31 31 H 1.08995 * 109.46819 * 59.99375 * 11 10 9 32 32 H 1.09000 * 109.46729 * 299.99852 * 11 10 9 33 33 H 1.09002 * 109.47108 * 305.00122 * 12 11 10 34 34 H 1.08997 * 109.47216 * 184.99925 * 12 11 10 35 35 H 1.00899 * 110.99743 * 51.10398 * 13 12 11 36 36 H 1.09004 * 109.68053 * 249.21099 * 18 17 15 37 37 H 1.08994 * 109.68555 * 9.76397 * 18 17 15 38 38 H 1.09004 * 109.65669 * 292.57294 * 19 18 17 39 39 H 1.09007 * 109.44596 * 172.37179 * 19 18 17 40 40 H 1.08990 * 109.51380 * 302.65235 * 21 20 19 41 41 H 1.09000 * 109.55007 * 182.52437 * 21 20 19 42 42 H 1.08998 * 109.67942 * 350.23751 * 22 17 15 43 43 H 1.09006 * 109.68032 * 110.75185 * 22 17 15 44 44 H 1.09007 * 109.47220 * 299.99743 * 23 14 13 45 45 H 1.09003 * 109.47204 * 60.00188 * 23 14 13 46 46 H 1.08999 * 109.47059 * 180.02562 * 24 23 14 47 47 H 1.09003 * 109.47298 * 59.99815 * 24 23 14 48 48 H 1.09003 * 109.47242 * 300.00069 * 25 24 23 49 49 H 1.09003 * 109.46846 * 179.97438 * 25 24 23 50 50 H 1.08995 * 109.46929 * 180.02562 * 26 25 24 51 51 H 1.09003 * 109.47204 * 59.99791 * 26 25 24 52 52 H 1.09003 * 109.47382 * 300.00813 * 27 14 13 53 53 H 1.09008 * 109.46770 * 60.00247 * 27 14 13 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.3151 0.0000 0.0000 3 14 2.2492 1.6177 0.0000 4 1 3.1029 1.6962 -1.2126 5 1 1.2843 2.7464 -0.0006 6 1 3.1029 1.6963 1.2125 7 6 2.0045 -1.1941 0.0005 8 1 3.0845 -1.1941 0.0001 9 6 1.2509 -2.4992 0.0005 10 8 2.1795 -3.5853 0.0017 11 6 1.5649 -4.8755 0.0013 12 6 2.6478 -5.9563 0.0033 13 7 3.4060 -5.8874 1.2596 14 6 4.3783 -6.9850 1.3486 15 6 3.6675 -8.3016 1.1689 16 8 4.0028 -9.0556 0.2801 17 7 2.6626 -8.6438 1.9994 18 6 2.3256 -7.8101 3.1643 19 6 0.8130 -7.5381 3.1372 20 8 0.1150 -8.7805 3.0231 21 6 0.3870 -9.4958 1.8155 22 6 1.8742 -9.8651 1.7717 23 6 5.4336 -6.8216 0.2529 24 6 6.4459 -7.9651 0.3451 25 6 7.1251 -7.9356 1.7158 26 6 6.0698 -8.0990 2.8115 27 6 5.0574 -6.9550 2.7192 28 1 -0.4850 0.8400 -0.0004 29 1 0.6246 -2.5571 -0.8897 30 1 0.6239 -2.5567 0.8903 31 1 0.9457 -4.9826 -0.8893 32 1 0.9438 -4.9824 0.8906 33 1 3.3226 -5.7960 -0.8376 34 1 2.1819 -6.9376 -0.0866 35 1 2.7867 -5.8876 2.0561 36 1 2.5904 -8.3364 4.0814 37 1 2.8702 -6.8675 3.1111 38 1 0.5696 -6.9042 2.2845 39 1 0.5176 -7.0380 4.0597 40 1 0.1422 -8.8691 0.9581 41 1 -0.2147 -10.4042 1.7862 42 1 2.1198 -10.2840 0.7959 43 1 2.0949 -10.5951 2.5506 44 1 4.9498 -6.8431 -0.7237 45 1 5.9478 -5.8693 0.3828 46 1 7.1980 -7.8484 -0.4352 47 1 5.9314 -8.9173 0.2149 48 1 7.6396 -6.9835 1.8459 49 1 7.8460 -8.7505 1.7817 50 1 6.5537 -8.0775 3.7880 51 1 5.5554 -9.0512 2.6814 52 1 5.5716 -6.0027 2.8492 53 1 4.3056 -7.0714 3.4998 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 ZINC000934055371.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:19 Heat of formation + Delta-G solvation = -107.804932 kcal Electronic energy + Delta-G solvation = -32661.804801 eV Core-core repulsion = 28534.167770 eV Total energy + Delta-G solvation = -4127.637031 eV No. of doubly occupied orbitals = 66 Molecular weight (most abundant/longest-lived isotopes) = 340.220 amu Computer time = 0.77 seconds Orbital eigenvalues (eV) -41.63408 -39.15847 -37.46503 -36.89369 -35.54738 -34.64137 -32.71035 -31.96091 -30.81959 -30.37931 -28.26625 -27.65120 -26.37164 -24.73330 -23.74581 -23.07749 -22.31242 -21.78298 -21.18937 -20.63312 -19.21433 -17.13624 -16.67283 -16.63112 -16.47531 -16.24994 -15.80074 -15.63961 -15.48334 -15.16249 -14.96604 -14.71445 -14.44781 -14.24502 -14.09592 -13.84336 -13.43828 -13.19750 -12.95216 -12.76482 -12.55361 -12.43821 -12.12872 -12.11704 -11.96658 -11.73557 -11.67604 -11.33635 -11.16185 -11.06766 -10.94032 -10.84013 -10.79697 -10.57545 -10.54196 -10.35763 -10.07038 -10.00078 -9.82096 -9.10818 -8.87477 -8.82053 -8.35298 -8.16384 -6.09963 -4.11184 2.35492 2.54116 3.33353 3.35875 3.40001 3.44489 3.85098 4.17110 4.32114 4.34904 4.43719 4.55557 4.60589 4.68358 4.78403 4.79317 4.87430 4.88894 5.02064 5.09256 5.15715 5.24165 5.26641 5.30464 5.32968 5.37180 5.42445 5.43667 5.46376 5.52260 5.58856 5.67070 5.79973 5.89286 5.93789 5.97994 6.02889 6.04459 6.07705 6.14599 6.18312 6.29810 6.48902 6.71408 6.83924 7.06694 7.08939 7.24842 7.42045 7.81231 7.91443 8.16359 8.58894 9.02477 9.53633 11.01630 Molecular weight = 340.22amu Principal moments of inertia in cm(-1) A = 0.011467 B = 0.003300 C = 0.002797 Principal moments of inertia in units of 10**(-40)*gram-cm**2 A = 2441.213564 B = 8483.217774 C =10006.859845 Net atomic charges, atomic populations, and dipole contributions calculated with CM2 Atom NO. Type Charge No. of electrons 1 N -0.900 5.900 2 C 0.056 3.944 3 Si 0.908 3.092 4 H -0.275 1.275 5 H -0.287 1.287 6 H -0.274 1.274 7 C -0.512 4.512 8 H 0.067 0.933 9 C 0.194 3.806 10 O -0.378 6.378 11 C 0.066 3.934 12 C 0.060 3.940 13 N -0.654 5.654 14 C 0.105 3.895 15 C 0.500 3.500 16 O -0.530 6.530 17 N -0.617 5.617 18 C 0.051 3.949 19 C 0.042 3.958 20 O -0.376 6.376 21 C 0.047 3.953 22 C 0.076 3.924 23 C -0.087 4.087 24 C -0.134 4.134 25 C -0.119 4.119 26 C -0.117 4.117 27 C -0.124 4.124 28 H 0.264 0.736 29 H 0.002 0.998 30 H 0.003 0.997 31 H 0.043 0.957 32 H 0.032 0.968 33 H 0.087 0.913 34 H 0.033 0.967 35 H 0.339 0.661 36 H 0.083 0.917 37 H 0.113 0.887 38 H 0.061 0.939 39 H 0.112 0.888 40 H 0.057 0.943 41 H 0.106 0.894 42 H 0.099 0.901 43 H 0.074 0.926 44 H 0.068 0.932 45 H 0.080 0.920 46 H 0.060 0.940 47 H 0.089 0.911 48 H 0.063 0.937 49 H 0.058 0.942 50 H 0.067 0.933 51 H 0.056 0.944 52 H 0.090 0.910 53 H 0.074 0.926 For the charges calculated by CM2: Dipole moment (debyes) X Y Z Total from point charges 6.379 -21.376 7.630 23.577 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.540 5.540 2 C -0.145 4.145 3 Si 0.735 3.265 4 H -0.200 1.200 5 H -0.212 1.212 6 H -0.200 1.200 7 C -0.550 4.550 8 H 0.085 0.915 9 C 0.120 3.880 10 O -0.298 6.298 11 C -0.012 4.012 12 C -0.069 4.069 13 N -0.290 5.290 14 C 0.011 3.989 15 C 0.287 3.713 16 O -0.405 6.405 17 N -0.353 5.353 18 C -0.072 4.072 19 C -0.036 4.036 20 O -0.294 6.294 21 C -0.031 4.031 22 C -0.046 4.046 23 C -0.125 4.125 24 C -0.171 4.171 25 C -0.157 4.157 26 C -0.155 4.155 27 C -0.161 4.161 28 H 0.073 0.927 29 H 0.020 0.980 30 H 0.021 0.979 31 H 0.061 0.939 32 H 0.050 0.950 33 H 0.106 0.894 34 H 0.051 0.949 35 H 0.163 0.837 36 H 0.101 0.899 37 H 0.131 0.869 38 H 0.079 0.921 39 H 0.130 0.870 40 H 0.075 0.925 41 H 0.124 0.876 42 H 0.118 0.882 43 H 0.092 0.908 44 H 0.086 0.914 45 H 0.099 0.901 46 H 0.079 0.921 47 H 0.107 0.893 48 H 0.082 0.918 49 H 0.077 0.923 50 H 0.085 0.915 51 H 0.075 0.925 52 H 0.109 0.891 53 H 0.092 0.908 Dipole moment (debyes) X Y Z Total from point charges 6.918 -21.222 7.033 23.403 hybrid contribution -0.579 0.211 0.039 0.618 sum 6.339 -21.011 7.071 23.058 Atomic orbital electron populations 1.69953 1.06114 1.27218 1.50686 1.24941 0.94951 0.99400 0.95186 0.86400 0.79306 0.81925 0.78873 1.20032 1.21230 1.19981 1.21015 0.94077 0.89273 1.50640 0.91507 1.20193 0.88831 0.83854 0.95142 1.87971 1.34606 1.12164 1.95009 1.22212 0.93594 0.86447 0.98936 1.22111 0.95893 0.95750 0.93132 1.60754 1.20590 1.42776 1.04891 1.21360 0.93713 0.85114 0.98756 1.22021 0.79244 0.89355 0.80666 1.90590 1.63402 1.53092 1.33369 1.48249 1.34272 1.25179 1.27608 1.23170 0.97537 0.97231 0.89288 1.23177 0.90264 0.86418 1.03724 1.87977 1.67607 1.33865 1.39989 1.22997 0.94821 0.97500 0.87765 1.22271 0.91200 0.87280 1.03892 1.21347 0.94887 0.99386 0.96876 1.21845 0.98065 1.00791 0.96435 1.21552 0.97431 1.01712 0.94963 1.21527 0.96900 0.97886 0.99233 1.22008 0.99031 1.01053 0.94049 0.92655 0.97987 0.97892 0.93925 0.95033 0.89438 0.94938 0.83679 0.89897 0.86935 0.92101 0.87021 0.92503 0.87590 0.88241 0.90794 0.91368 0.90126 0.92095 0.89319 0.91797 0.92278 0.91461 0.92479 0.89132 0.90806 Number of geometries 1 Number of calculations of the screened coulomb radii 1 The total number of SCF iterations 13. 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.81 13.29 55.43 0.74 -26.08 16 2 C 0.06 1.62 5.46 -17.82 -0.10 1.52 16 3 Si 0.91 21.80 29.54 -169.99 -5.02 16.78 16 4 H -0.27 -6.50 7.11 56.52 0.40 -6.10 16 5 H -0.29 -6.73 7.11 56.52 0.40 -6.33 16 6 H -0.27 -6.45 7.11 56.52 0.40 -6.05 16 7 C -0.51 -15.00 9.43 -34.44 -0.32 -15.33 16 8 H 0.07 1.95 7.81 -52.49 -0.41 1.54 16 9 C 0.19 5.27 5.69 36.01 0.20 5.47 16 10 O -0.38 -8.95 10.27 -35.23 -0.36 -9.31 16 11 C 0.07 1.20 6.17 37.15 0.23 1.43 16 12 C 0.06 0.91 5.42 -3.83 -0.02 0.89 16 13 N -0.65 -9.02 6.04 -101.15 -0.61 -9.63 16 14 C 0.11 1.27 0.42 -132.76 -0.06 1.21 16 15 C 0.50 6.23 3.37 -10.99 -0.04 6.19 16 16 O -0.53 -7.93 13.08 5.56 0.07 -7.85 16 17 N -0.62 -6.40 3.01 -172.02 -0.52 -6.92 16 18 C 0.05 0.44 4.54 -3.31 -0.02 0.43 16 19 C 0.04 0.38 6.95 37.53 0.26 0.64 16 20 O -0.38 -3.82 10.74 -35.23 -0.38 -4.20 16 21 C 0.05 0.41 6.99 37.31 0.26 0.67 16 22 C 0.08 0.71 6.61 -3.54 -0.02 0.68 16 23 C -0.09 -1.12 4.28 -26.73 -0.11 -1.23 16 24 C -0.13 -1.54 5.63 -26.73 -0.15 -1.69 16 25 C -0.12 -1.15 5.92 -26.73 -0.16 -1.31 16 26 C -0.12 -1.09 5.70 -26.71 -0.15 -1.24 16 27 C -0.12 -1.23 4.63 -26.71 -0.12 -1.35 16 28 H 0.26 7.38 8.31 -40.82 -0.34 7.04 16 29 H 0.00 0.07 8.14 -51.93 -0.42 -0.36 16 30 H 0.00 0.09 8.14 -51.93 -0.42 -0.33 16 31 H 0.04 0.75 8.14 -51.93 -0.42 0.33 16 32 H 0.03 0.53 8.12 -51.93 -0.42 0.11 16 33 H 0.09 1.43 6.69 -51.93 -0.35 1.08 16 34 H 0.03 0.47 6.76 -51.93 -0.35 0.12 16 35 H 0.34 4.24 5.29 -39.29 -0.21 4.04 16 36 H 0.08 0.61 8.14 -51.93 -0.42 0.19 16 37 H 0.11 0.99 2.49 -51.93 -0.13 0.86 16 38 H 0.06 0.62 8.14 -51.93 -0.42 0.19 16 39 H 0.11 0.81 8.14 -51.93 -0.42 0.39 16 40 H 0.06 0.53 8.14 -51.93 -0.42 0.11 16 41 H 0.11 0.69 8.14 -51.93 -0.42 0.27 16 42 H 0.10 1.00 7.13 -51.93 -0.37 0.63 16 43 H 0.07 0.62 8.14 -51.93 -0.42 0.19 16 44 H 0.07 0.97 6.61 -51.92 -0.34 0.63 16 45 H 0.08 1.06 8.14 -51.93 -0.42 0.64 16 46 H 0.06 0.65 8.14 -51.93 -0.42 0.23 16 47 H 0.09 1.18 5.43 -51.93 -0.28 0.90 16 48 H 0.06 0.60 8.14 -51.93 -0.42 0.17 16 49 H 0.06 0.52 8.14 -51.93 -0.42 0.10 16 50 H 0.07 0.53 8.14 -51.93 -0.42 0.10 16 51 H 0.06 0.59 7.22 -51.93 -0.37 0.22 16 52 H 0.09 0.93 8.14 -51.93 -0.42 0.51 16 53 H 0.07 0.63 4.72 -51.92 -0.25 0.38 16 LS Contribution 393.22 15.07 5.93 5.93 Total: -1.00 -33.08 393.22 -9.43 -42.52 By element: Atomic # 1 Polarization: 10.76 SS G_CDS: -8.96 Total: 1.80 kcal Atomic # 6 Polarization: -2.72 SS G_CDS: -0.32 Total: -3.04 kcal Atomic # 7 Polarization: -42.24 SS G_CDS: -0.39 Total: -42.63 kcal Atomic # 8 Polarization: -20.69 SS G_CDS: -0.67 Total: -21.36 kcal Atomic # 14 Polarization: 21.80 SS G_CDS: -5.02 Total: 16.78 kcal Total LS contribution 5.93 Total: 5.93 kcal Total: -33.08 -9.43 -42.52 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. ZINC000934055371.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 -65.287 kcal (2) G-P(sol) polarization free energy of solvation -33.085 kcal (3) G-ENP(sol) elect.-nuc.-pol. free energy of system -98.372 kcal (4) G-CDS(sol) cavity-dispersion-solvent structure free energy -9.433 kcal (5) G-P-CDS(sol) = G-P(sol) + G-CDS(sol) = (2) + (4) -42.518 kcal (6) G-S(sol) free energy of system = (1) + (5) -107.805 kcal 1SCF WAS SPECIFIED, SO GEOMETRY OPTIMIZATION WAS NOT USED Total computer time = 0.77 seconds