Wall clock time and date at job start Tue Mar 31 2020 12:42:34 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.53004 * 1 3 3 C 1.53001 * 109.46750 * 2 1 4 4 C 1.52995 * 109.47123 * 119.99862 * 2 1 3 5 5 C 1.53006 * 109.46994 * 83.18836 * 4 2 1 6 6 C 1.53004 * 109.47055 * 302.75350 * 5 4 2 7 7 C 1.52997 * 109.47006 * 182.75371 * 5 4 2 8 8 C 1.53000 * 109.47127 * 62.75013 * 5 4 2 9 9 N 1.46894 * 109.47096 * 179.97438 * 8 5 4 10 10 C 1.46903 * 111.00009 * 180.02562 * 9 8 5 11 11 H 1.09001 * 109.47088 * 325.00185 * 10 9 8 12 12 C 1.52997 * 109.47124 * 84.99960 * 10 9 8 13 13 C 1.50695 * 109.47292 * 205.00415 * 10 9 8 14 14 C 1.38310 * 119.92839 * 319.72726 * 13 10 9 15 15 C 1.38096 * 120.06718 * 179.72802 * 14 13 10 16 16 C 1.38888 * 119.93163 * 0.54154 * 15 14 13 17 17 O 1.35750 * 120.06900 * 179.74653 * 16 15 14 18 18 C 1.34830 * 116.99831 * 174.42582 * 17 16 15 19 19 H 1.07999 * 120.00050 * 5.05323 * 18 17 16 20 20 C 1.38651 * 119.99846 * 185.05780 * 18 17 16 21 21 N 1.31559 * 120.00427 * 359.97438 * 20 18 17 22 22 Si 1.86798 * 119.99820 * 179.97438 * 20 18 17 23 23 H 1.48500 * 109.47340 * 59.99861 * 22 20 18 24 24 H 1.48504 * 109.47330 * 179.97438 * 22 20 18 25 25 H 1.48502 * 109.47296 * 299.99640 * 22 20 18 26 26 C 1.38877 * 119.85955 * 359.73039 * 16 15 14 27 27 C 1.38103 * 119.93474 * 359.97438 * 26 16 15 28 28 H 1.08991 * 109.47652 * 305.63991 * 1 2 3 29 29 H 1.09008 * 109.46893 * 65.64283 * 1 2 3 30 30 H 1.09002 * 109.46967 * 185.63687 * 1 2 3 31 31 H 1.09003 * 109.47332 * 240.00352 * 2 1 3 32 32 H 1.08997 * 109.47013 * 179.97438 * 3 2 1 33 33 H 1.08994 * 109.47177 * 299.99990 * 3 2 1 34 34 H 1.09003 * 109.46915 * 60.00097 * 3 2 1 35 35 H 1.08998 * 109.47134 * 323.19237 * 4 2 1 36 36 H 1.08992 * 109.47558 * 203.18845 * 4 2 1 37 37 H 1.09001 * 109.46623 * 189.64571 * 6 5 4 38 38 H 1.08996 * 109.47008 * 309.64411 * 6 5 4 39 39 H 1.09000 * 109.47253 * 69.64445 * 6 5 4 40 40 H 1.08998 * 109.46947 * 295.00449 * 7 5 4 41 41 H 1.09000 * 109.47366 * 55.00951 * 7 5 4 42 42 H 1.09005 * 109.46948 * 175.01111 * 7 5 4 43 43 H 1.09000 * 109.46839 * 59.99882 * 8 5 4 44 44 H 1.08998 * 109.47162 * 300.00288 * 8 5 4 45 45 H 1.00904 * 110.99946 * 56.05077 * 9 8 5 46 46 H 1.09001 * 109.46954 * 60.00143 * 12 10 9 47 47 H 1.08996 * 109.47076 * 180.02562 * 12 10 9 48 48 H 1.08998 * 109.47402 * 300.00047 * 12 10 9 49 49 H 1.07997 * 119.96309 * 359.97438 * 14 13 10 50 50 H 1.07997 * 120.03658 * 180.26879 * 15 14 13 51 51 H 0.96998 * 120.00193 * 180.02562 * 21 20 18 52 52 H 1.07998 * 120.03664 * 179.97438 * 26 16 15 53 53 H 1.07996 * 119.96653 * 180.02562 * 27 26 16 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 6 2.0399 1.4425 0.0000 4 6 2.0400 -0.7212 1.2492 5 6 2.0487 -2.2306 0.9989 6 6 0.6415 -2.6918 0.6141 7 6 2.4932 -2.9559 2.2705 8 6 3.0199 -2.5534 -0.1384 9 7 3.0287 -4.0025 -0.3784 10 6 3.9528 -4.3505 -1.4660 11 1 3.9425 -3.5614 -2.2179 12 6 5.3680 -4.5012 -0.9043 13 6 3.5212 -5.6498 -2.0956 14 6 3.0641 -6.6841 -1.2993 15 6 2.6625 -7.8741 -1.8736 16 6 2.7287 -8.0356 -3.2515 17 8 2.3395 -9.2059 -3.8186 18 6 2.3195 -9.2564 -5.1658 19 1 2.6824 -8.4189 -5.7430 20 6 1.8330 -10.3836 -5.8101 21 7 1.3912 -11.4041 -5.1070 22 14 1.8046 -10.4533 -7.6766 23 1 3.1876 -10.3107 -8.1984 24 1 1.2426 -11.7562 -8.1148 25 1 0.9615 -9.3487 -8.2005 26 6 3.1935 -6.9971 -4.0479 27 6 3.5878 -5.8073 -3.4680 28 1 -0.3634 0.5987 0.8351 29 1 -0.3633 0.4239 -0.9363 30 1 -0.3633 -1.0227 0.1009 31 1 1.8934 -0.5138 -0.8900 32 1 3.1299 1.4426 0.0005 33 1 1.6766 1.9563 -0.8899 34 1 1.6766 1.9563 0.8900 35 1 1.3854 -0.4964 2.0912 36 1 3.0517 -0.3847 1.4755 37 1 0.6067 -3.7812 0.6030 38 1 -0.0766 -2.3127 1.3411 39 1 0.3921 -2.3102 -0.3760 40 1 3.5274 -2.6948 2.4950 41 1 1.8543 -2.6570 3.1015 42 1 2.4147 -4.0328 2.1209 43 1 4.0224 -2.2247 0.1357 44 1 2.7032 -2.0366 -1.0444 45 1 2.0985 -4.3428 -0.5711 46 1 5.3783 -5.2903 -0.1524 47 1 6.0537 -4.7589 -1.7114 48 1 5.6802 -3.5614 -0.4490 49 1 3.0174 -6.5596 -0.2275 50 1 2.3016 -8.6798 -1.2517 51 1 1.0512 -12.1928 -5.5578 52 1 3.2465 -7.1197 -5.1196 53 1 3.9490 -4.9993 -4.0869 RHF calculation, no. of doubly occupied orbitals= 62 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 ZINC000391128680.mol2 53 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 12:42:34 Heat of formation + Delta-G solvation = -26.117343 kcal Electronic energy + Delta-G solvation = -27198.041663 eV Core-core repulsion = 23648.511026 eV Total energy + Delta-G solvation = -3549.530637 eV No. of doubly occupied orbitals = 62 Molecular weight (most abundant/longest-lived isotopes) = 319.225 amu Computer time = 0.90 seconds Orbital eigenvalues (eV) -39.56726 -38.96163 -37.24029 -35.40001 -33.73215 -32.23833 -30.31672 -29.95819 -28.51475 -27.90190 -26.90311 -26.41253 -25.44589 -24.94808 -22.26911 -21.73929 -20.57832 -19.90756 -18.78981 -17.92326 -17.01297 -15.97228 -15.82924 -15.25038 -14.85202 -14.54620 -14.44775 -14.17613 -14.06506 -13.71243 -13.63829 -13.48238 -13.28522 -12.77395 -12.74912 -12.52073 -12.46627 -12.34198 -12.15206 -12.03524 -11.94159 -11.84327 -11.74283 -11.65173 -11.34473 -11.17820 -10.96647 -10.91523 -10.77329 -10.66062 -10.57768 -10.45291 -10.41688 -10.21305 -10.08840 -9.33346 -8.82100 -8.25051 -8.15583 -7.16548 -6.44806 -4.03796 2.23777 2.32003 3.12756 3.13660 3.16548 3.96666 4.24344 4.47674 4.57701 4.60602 4.74277 4.82450 5.00062 5.00715 5.09421 5.15968 5.33179 5.37156 5.37289 5.38560 5.49754 5.53098 5.55512 5.63486 5.66264 5.68746 5.71337 5.72300 5.75956 5.76725 5.86388 5.89382 5.90842 5.97562 6.01811 6.06236 6.06967 6.14981 6.24361 6.25561 6.30694 6.38344 6.42602 6.50404 6.74608 6.89177 7.05445 7.10842 7.28425 7.44352 7.46755 7.60482 7.84694 8.00183 9.00740 9.71361 10.73156 Molecular weight = 319.23amu Principal moments of inertia in cm(-1) A = 0.019275 B = 0.002377 C = 0.002313 Principal moments of inertia in units of 10**(-40)*gram-cm**2 A = 1452.296360 B =11776.621994 C =12104.703271 Net atomic charges, atomic populations, and dipole contributions calculated with CM2 Atom NO. Type Charge No. of electrons 1 C -0.170 4.170 2 C -0.102 4.102 3 C -0.147 4.147 4 C -0.103 4.103 5 C -0.061 4.061 6 C -0.152 4.152 7 C -0.130 4.130 8 C 0.068 3.932 9 N -0.681 5.681 10 C 0.143 3.857 11 H 0.038 0.962 12 C -0.131 4.131 13 C -0.191 4.191 14 C -0.055 4.055 15 C -0.195 4.195 16 C 0.155 3.845 17 O -0.178 6.178 18 C -0.392 4.392 19 H 0.091 0.909 20 C 0.039 3.961 21 N -0.863 5.863 22 Si 0.940 3.060 23 H -0.270 1.270 24 H -0.278 1.278 25 H -0.270 1.270 26 C -0.235 4.235 27 C -0.066 4.066 28 H 0.049 0.951 29 H 0.048 0.952 30 H 0.079 0.921 31 H 0.074 0.926 32 H 0.052 0.948 33 H 0.051 0.949 34 H 0.051 0.949 35 H 0.062 0.938 36 H 0.060 0.940 37 H 0.059 0.941 38 H 0.056 0.944 39 H 0.062 0.938 40 H 0.049 0.951 41 H 0.043 0.957 42 H 0.068 0.932 43 H 0.067 0.933 44 H 0.029 0.971 45 H 0.344 0.656 46 H 0.061 0.939 47 H 0.052 0.948 48 H 0.044 0.956 49 H 0.119 0.881 50 H 0.116 0.884 51 H 0.270 0.730 52 H 0.125 0.875 53 H 0.103 0.897 For the charges calculated by CM2: Dipole moment (debyes) X Y Z Total from point charges 1.477 21.286 10.859 23.942 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.227 4.227 2 C -0.120 4.120 3 C -0.204 4.204 4 C -0.141 4.141 5 C -0.061 4.061 6 C -0.208 4.208 7 C -0.187 4.187 8 C -0.062 4.062 9 N -0.314 5.314 10 C 0.033 3.967 11 H 0.056 0.944 12 C -0.189 4.189 13 C -0.192 4.192 14 C -0.074 4.074 15 C -0.214 4.214 16 C 0.101 3.899 17 O -0.081 6.081 18 C -0.468 4.468 19 H 0.108 0.892 20 C -0.165 4.165 21 N -0.498 5.498 22 Si 0.765 3.235 23 H -0.195 1.195 24 H -0.203 1.203 25 H -0.195 1.195 26 C -0.254 4.254 27 C -0.085 4.085 28 H 0.068 0.932 29 H 0.067 0.933 30 H 0.098 0.902 31 H 0.092 0.908 32 H 0.071 0.929 33 H 0.070 0.930 34 H 0.070 0.930 35 H 0.081 0.919 36 H 0.079 0.921 37 H 0.078 0.922 38 H 0.075 0.925 39 H 0.080 0.920 40 H 0.068 0.932 41 H 0.062 0.938 42 H 0.087 0.913 43 H 0.086 0.914 44 H 0.047 0.953 45 H 0.164 0.836 46 H 0.080 0.920 47 H 0.071 0.929 48 H 0.063 0.937 49 H 0.137 0.863 50 H 0.134 0.866 51 H 0.080 0.920 52 H 0.143 0.857 53 H 0.121 0.879 Dipole moment (debyes) X Y Z Total from point charges 1.762 20.693 11.397 23.689 hybrid contribution -0.445 0.372 -1.856 1.945 sum 1.316 21.065 9.541 23.162 Atomic orbital electron populations 1.22405 0.95452 1.03653 1.01191 1.21229 0.95975 0.95833 0.98977 1.21864 1.00742 0.96369 1.01424 1.21725 1.01033 0.93930 0.97419 1.20553 0.95131 0.94525 0.95900 1.22218 0.95897 1.01473 1.01261 1.21725 1.00223 1.01401 0.95353 1.22117 0.99020 0.88701 0.96323 1.60616 1.13523 1.12551 1.44697 1.20124 0.92737 0.93518 0.90280 0.94407 1.21756 0.95156 1.01147 1.00803 1.19283 1.10479 0.96241 0.93245 1.20487 0.96471 0.91659 0.98735 1.20256 1.09790 0.99027 0.92373 1.19425 0.91876 0.86218 0.92383 1.83980 1.73821 1.27504 1.22826 1.21265 1.43450 1.03748 0.78380 0.89188 1.25280 0.96042 0.94697 1.00486 1.70028 1.38456 1.05196 1.36101 0.85794 0.78960 0.78061 0.80687 1.19525 1.20266 1.19531 1.21110 1.09120 0.94671 1.00500 1.20583 0.97786 0.96358 0.93734 0.93182 0.93326 0.90203 0.90803 0.92879 0.93007 0.92965 0.91898 0.92089 0.92189 0.92455 0.91959 0.93167 0.93779 0.91327 0.91439 0.95270 0.83558 0.92036 0.92871 0.93651 0.86315 0.86630 0.91975 0.85701 0.87885 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.17 -1.29 8.62 37.16 0.32 -0.97 16 2 C -0.10 -0.79 2.45 -90.62 -0.22 -1.01 16 3 C -0.15 -1.04 9.08 37.16 0.34 -0.71 16 4 C -0.10 -0.80 4.01 -26.73 -0.11 -0.91 16 5 C -0.06 -0.52 0.62 -154.50 -0.10 -0.62 16 6 C -0.15 -1.33 6.09 37.16 0.23 -1.11 16 7 C -0.13 -1.17 8.23 37.16 0.31 -0.87 16 8 C 0.07 0.67 4.20 -3.83 -0.02 0.66 16 9 N -0.68 -8.01 5.01 -107.96 -0.54 -8.55 16 10 C 0.14 1.91 2.81 -68.87 -0.19 1.71 16 11 H 0.04 0.51 7.66 -51.93 -0.40 0.11 16 12 C -0.13 -1.65 9.39 37.16 0.35 -1.30 16 13 C -0.19 -3.18 4.71 -104.57 -0.49 -3.67 16 14 C -0.05 -0.97 9.00 -39.61 -0.36 -1.32 16 15 C -0.20 -4.11 9.99 -39.39 -0.39 -4.50 16 16 C 0.15 3.73 6.65 -39.10 -0.26 3.47 16 17 O -0.18 -5.09 9.76 0.05 0.00 -5.09 16 18 C -0.39 -10.98 9.49 30.88 0.29 -10.69 16 19 H 0.09 2.33 5.20 -52.49 -0.27 2.06 16 20 C 0.04 1.08 5.49 -17.52 -0.10 0.98 16 21 N -0.86 -25.36 13.67 55.48 0.76 -24.60 16 22 Si 0.94 20.54 29.55 -169.99 -5.02 15.52 16 23 H -0.27 -5.71 7.11 56.52 0.40 -5.31 16 24 H -0.28 -5.93 7.11 56.52 0.40 -5.52 16 25 H -0.27 -5.72 7.11 56.52 0.40 -5.32 16 26 C -0.23 -5.08 9.17 -39.39 -0.36 -5.44 16 27 C -0.07 -1.20 9.66 -39.61 -0.38 -1.58 16 28 H 0.05 0.36 8.14 -51.93 -0.42 -0.06 16 29 H 0.05 0.37 8.14 -51.92 -0.42 -0.05 16 30 H 0.08 0.62 3.74 -51.93 -0.19 0.42 16 31 H 0.07 0.64 5.76 -51.93 -0.30 0.34 16 32 H 0.05 0.37 8.12 -51.93 -0.42 -0.06 16 33 H 0.05 0.36 8.14 -51.93 -0.42 -0.06 16 34 H 0.05 0.35 8.14 -51.93 -0.42 -0.07 16 35 H 0.06 0.48 7.94 -51.93 -0.41 0.07 16 36 H 0.06 0.47 7.92 -51.93 -0.41 0.06 16 37 H 0.06 0.58 7.04 -51.93 -0.37 0.21 16 38 H 0.06 0.46 7.39 -51.93 -0.38 0.08 16 39 H 0.06 0.56 6.00 -51.93 -0.31 0.24 16 40 H 0.05 0.46 8.14 -51.93 -0.42 0.03 16 41 H 0.04 0.37 8.14 -51.93 -0.42 -0.06 16 42 H 0.07 0.70 7.74 -51.93 -0.40 0.30 16 43 H 0.07 0.69 7.43 -51.93 -0.39 0.30 16 44 H 0.03 0.30 6.04 -51.93 -0.31 -0.01 16 45 H 0.34 4.30 6.27 -39.29 -0.25 4.05 16 46 H 0.06 0.82 8.14 -51.93 -0.42 0.40 16 47 H 0.05 0.67 8.14 -51.93 -0.42 0.25 16 48 H 0.04 0.51 7.45 -51.93 -0.39 0.12 16 49 H 0.12 1.85 7.72 -52.49 -0.41 1.45 16 50 H 0.12 2.37 8.06 -52.49 -0.42 1.95 16 51 H 0.27 7.30 8.31 -40.82 -0.34 6.96 16 52 H 0.13 2.72 5.64 -52.49 -0.30 2.43 16 53 H 0.10 1.62 8.03 -52.49 -0.42 1.20 16 LS Contribution 403.57 15.07 6.08 6.08 Total: -1.00 -28.86 403.57 -9.14 -38.00 By element: Atomic # 1 Polarization: 15.78 SS G_CDS: -9.27 Total: 6.51 kcal Atomic # 6 Polarization: -26.72 SS G_CDS: -1.15 Total: -27.86 kcal Atomic # 7 Polarization: -33.37 SS G_CDS: 0.22 Total: -33.15 kcal Atomic # 8 Polarization: -5.09 SS G_CDS: 0.00 Total: -5.09 kcal Atomic # 14 Polarization: 20.54 SS G_CDS: -5.02 Total: 15.52 kcal Total LS contribution 6.08 Total: 6.08 kcal Total: -28.86 -9.14 -38.00 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. ZINC000391128680.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 11.882 kcal (2) G-P(sol) polarization free energy of solvation -28.862 kcal (3) G-ENP(sol) elect.-nuc.-pol. free energy of system -16.981 kcal (4) G-CDS(sol) cavity-dispersion-solvent structure free energy -9.137 kcal (5) G-P-CDS(sol) = G-P(sol) + G-CDS(sol) = (2) + (4) -37.999 kcal (6) G-S(sol) free energy of system = (1) + (5) -26.117 kcal 1SCF WAS SPECIFIED, SO GEOMETRY OPTIMIZATION WAS NOT USED Total computer time = 0.90 seconds