Wall clock time and date at job start Sat Apr 11 2020 06:12:40 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.52992 * 1 3 3 C 1.53003 * 109.47066 * 2 1 4 4 C 1.52999 * 109.46984 * 240.00356 * 2 1 3 5 5 O 1.45199 * 109.47503 * 120.00344 * 2 1 3 6 6 C 1.34639 * 117.00000 * 59.99932 * 5 2 1 7 7 O 1.21511 * 119.99714 * 0.02562 * 6 5 2 8 8 N 1.34776 * 119.99974 * 180.02562 * 6 5 2 9 9 C 1.47155 * 120.98974 * 359.97438 * 8 6 5 10 10 C 1.53292 * 108.38071 * 230.60338 * 9 8 6 11 11 O 1.42969 * 109.23581 * 307.66979 * 10 9 8 12 12 C 1.42892 * 114.19470 * 62.49154 * 11 10 9 13 13 C 1.52960 * 109.57672 * 177.70809 * 12 11 10 14 14 C 1.52997 * 109.51340 * 299.56226 * 13 12 11 15 15 C 1.53000 * 109.47060 * 300.07539 * 14 13 12 16 16 N 1.46507 * 109.47099 * 180.02562 * 15 14 13 17 17 C 1.36943 * 119.99574 * 155.00276 * 16 15 14 18 18 H 1.07995 * 119.99945 * 359.97438 * 17 16 15 19 19 C 1.38815 * 119.99331 * 179.97438 * 17 16 15 20 20 N 1.31624 * 119.99339 * 179.97438 * 19 17 16 21 21 Si 1.86795 * 120.00380 * 359.97438 * 19 17 16 22 22 H 1.48496 * 109.47531 * 89.99912 * 21 19 17 23 23 H 1.48497 * 109.47172 * 210.00245 * 21 19 17 24 24 H 1.48505 * 109.47190 * 329.99945 * 21 19 17 25 25 C 1.53004 * 109.45117 * 60.05665 * 15 14 13 26 26 C 1.52946 * 109.58143 * 57.46052 * 12 11 10 27 27 C 1.47158 * 120.98550 * 180.02562 * 8 6 5 28 28 H 1.09002 * 109.46913 * 299.99785 * 1 2 3 29 29 H 1.08998 * 109.47722 * 59.99673 * 1 2 3 30 30 H 1.09002 * 109.47400 * 179.97438 * 1 2 3 31 31 H 1.09004 * 109.47227 * 180.02562 * 3 2 1 32 32 H 1.08989 * 109.47415 * 300.00215 * 3 2 1 33 33 H 1.09001 * 109.47220 * 60.01087 * 3 2 1 34 34 H 1.08993 * 109.47358 * 299.99843 * 4 2 1 35 35 H 1.09002 * 109.47175 * 60.00178 * 4 2 1 36 36 H 1.09000 * 109.47652 * 179.97438 * 4 2 1 37 37 H 1.09001 * 109.67603 * 110.88138 * 9 8 6 38 38 H 1.08993 * 109.71282 * 350.34569 * 9 8 6 39 39 H 1.08995 * 109.50854 * 187.74250 * 10 9 8 40 40 H 1.09003 * 109.55201 * 67.62415 * 10 9 8 41 41 H 1.09000 * 109.46397 * 179.54940 * 13 12 11 42 42 H 1.09000 * 109.46093 * 59.57623 * 13 12 11 43 43 H 1.09002 * 109.47153 * 180.02562 * 14 13 12 44 44 H 1.09001 * 109.46832 * 60.07410 * 14 13 12 45 45 H 1.08995 * 109.47704 * 300.06073 * 15 14 13 46 46 H 0.96996 * 120.00207 * 334.99250 * 16 15 14 47 47 H 0.97002 * 119.99331 * 179.97438 * 20 19 17 48 48 H 1.09000 * 109.46878 * 59.93585 * 25 15 14 49 49 H 1.08989 * 109.47330 * 179.97438 * 25 15 14 50 50 H 1.09010 * 109.46276 * 300.41250 * 26 12 11 51 51 H 1.09003 * 109.46732 * 180.44203 * 26 12 11 52 52 H 1.09001 * 109.68114 * 9.76268 * 27 8 6 53 53 H 1.09000 * 109.67889 * 249.20775 * 27 8 6 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.5299 0.0000 0.0000 3 6 2.0399 1.4425 0.0000 4 6 2.0399 -0.7212 -1.2493 5 8 2.0140 -0.6845 1.1855 6 6 1.6523 -1.9724 1.3382 7 8 0.9516 -2.5097 0.5034 8 7 2.0673 -2.6655 2.4170 9 6 2.9168 -2.0366 3.4409 10 6 4.1252 -2.9458 3.6915 11 8 3.6675 -4.2706 3.9734 12 6 2.9525 -4.8841 2.8991 13 6 2.5637 -6.3107 3.2908 14 6 3.8273 -7.1241 3.5781 15 6 4.7066 -7.1594 2.3266 16 7 5.9163 -7.9388 2.6016 17 6 6.5895 -8.5436 1.5737 18 1 6.2284 -8.4463 0.5606 19 6 7.7360 -9.2816 1.8343 20 7 8.3827 -9.8633 0.8464 21 14 8.3609 -9.4492 3.5866 22 1 7.7474 -10.6447 4.2188 23 1 9.8381 -9.5996 3.5731 24 1 7.9939 -8.2364 4.3611 25 6 5.0974 -5.7321 1.9379 26 6 3.8343 -4.9183 1.6499 27 6 1.6854 -4.0742 2.6053 28 1 -0.3633 0.5138 0.8900 29 1 -0.3634 0.5139 -0.8899 30 1 -0.3634 -1.0277 -0.0005 31 1 3.1300 1.4426 -0.0005 32 1 1.6766 1.9564 -0.8899 33 1 1.6767 1.9563 0.8901 34 1 1.6766 -1.7488 -1.2494 35 1 1.6766 -0.2073 -2.1393 36 1 3.1299 -0.7208 -1.2496 37 1 2.3502 -1.9160 4.3642 38 1 3.2568 -1.0627 3.0889 39 1 4.6943 -2.5688 4.5412 40 1 4.7604 -2.9608 2.8058 41 1 2.0115 -6.7741 2.4732 42 1 1.9378 -6.2846 4.1828 43 1 3.5495 -8.1406 3.8568 44 1 4.3787 -6.6615 4.3967 45 1 4.1549 -7.6209 1.5077 46 1 6.2405 -8.0264 3.5116 47 1 9.1840 -10.3787 1.0284 48 1 5.6481 -5.2703 2.7574 49 1 5.7242 -5.7569 1.0467 50 1 4.1126 -3.9011 1.3738 51 1 3.2843 -5.3795 0.8295 52 1 1.2113 -4.4504 1.6988 53 1 0.9940 -4.1580 3.4438 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 ZINC001254712930.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:12:40 Heat of formation + Delta-G solvation = -112.831662 kcal Electronic energy + Delta-G solvation = -32215.468387 eV Core-core repulsion = 28087.613380 eV Total energy + Delta-G solvation = -4127.855007 eV No. of doubly occupied orbitals = 66 Molecular weight (most abundant/longest-lived isotopes) = 340.220 amu Computer time = 1.11 seconds Orbital eigenvalues (eV) -41.29655 -40.15301 -37.55386 -36.69556 -35.05328 -34.06735 -33.82648 -31.72097 -30.58467 -30.01020 -28.34665 -27.71809 -27.64429 -25.90620 -23.79964 -23.05782 -22.55014 -21.64795 -20.25615 -19.54202 -19.27077 -18.22363 -17.54978 -16.94786 -16.54599 -15.97732 -15.68849 -15.56678 -15.06880 -14.82561 -14.73645 -14.47684 -14.06511 -13.97155 -13.74327 -13.52455 -13.04008 -12.88106 -12.74343 -12.57586 -12.46981 -12.39705 -12.18786 -12.07417 -12.03480 -11.90982 -11.85307 -11.70283 -11.57039 -11.23906 -11.01736 -10.89601 -10.75768 -10.63617 -10.43927 -10.36164 -10.17999 -10.02764 -10.00146 -9.54931 -9.53089 -8.85016 -8.70836 -7.20797 -5.72824 -3.07387 1.84362 2.63065 3.41100 3.44724 3.48418 3.50434 3.65160 4.41240 4.45992 4.47725 4.60445 4.63341 4.68220 4.75468 4.81308 4.95701 5.06954 5.13005 5.15799 5.24202 5.29320 5.35495 5.38147 5.44763 5.46020 5.50684 5.55299 5.57411 5.59910 5.62278 5.67964 5.69839 5.81616 5.82834 5.93564 6.00301 6.03973 6.08493 6.15572 6.20416 6.22514 6.32484 6.38359 6.58318 6.66997 6.95245 7.08649 7.10964 7.20131 8.04835 8.11246 8.46959 9.26250 9.99131 10.19820 11.45794 Molecular weight = 340.22amu Principal moments of inertia in cm(-1) A = 0.019104 B = 0.002686 C = 0.002574 Principal moments of inertia in units of 10**(-40)*gram-cm**2 A = 1465.305657 B =10420.562322 C =10875.688366 Net atomic charges, atomic populations, and dipole contributions calculated with CM2 Atom NO. Type Charge No. of electrons 1 C -0.183 4.183 2 C 0.136 3.864 3 C -0.140 4.140 4 C -0.184 4.184 5 O -0.364 6.364 6 C 0.655 3.345 7 O -0.566 6.566 8 N -0.595 5.595 9 C 0.079 3.921 10 C 0.043 3.957 11 O -0.359 6.359 12 C 0.109 3.891 13 C -0.096 4.096 14 C -0.122 4.122 15 C 0.184 3.816 16 N -0.749 5.749 17 C -0.239 4.239 18 H 0.089 0.911 19 C 0.002 3.998 20 N -0.936 5.936 21 Si 0.884 3.116 22 H -0.285 1.285 23 H -0.292 1.292 24 H -0.276 1.276 25 C -0.129 4.129 26 C -0.137 4.137 27 C 0.106 3.894 28 H 0.058 0.942 29 H 0.070 0.930 30 H 0.088 0.912 31 H 0.067 0.933 32 H 0.080 0.920 33 H 0.065 0.935 34 H 0.088 0.912 35 H 0.070 0.930 36 H 0.060 0.940 37 H 0.076 0.924 38 H 0.098 0.902 39 H 0.105 0.895 40 H 0.063 0.937 41 H 0.079 0.921 42 H 0.069 0.931 43 H 0.064 0.936 44 H 0.067 0.933 45 H 0.052 0.948 46 H 0.360 0.640 47 H 0.258 0.742 48 H 0.058 0.942 49 H 0.063 0.937 50 H 0.056 0.944 51 H 0.070 0.930 52 H 0.105 0.895 53 H 0.077 0.923 For the charges calculated by CM2: Dipole moment (debyes) X Y Z Total from point charges -16.784 23.198 3.622 28.861 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.240 4.240 2 C 0.099 3.901 3 C -0.198 4.198 4 C -0.242 4.242 5 O -0.279 6.279 6 C 0.411 3.589 7 O -0.455 6.455 8 N -0.336 5.336 9 C -0.044 4.044 10 C -0.035 4.035 11 O -0.278 6.278 12 C 0.069 3.931 13 C -0.134 4.134 14 C -0.161 4.161 15 C 0.075 3.925 16 N -0.391 5.391 17 C -0.369 4.369 18 H 0.107 0.893 19 C -0.194 4.194 20 N -0.585 5.585 21 Si 0.715 3.285 22 H -0.211 1.211 23 H -0.219 1.219 24 H -0.203 1.203 25 C -0.169 4.169 26 C -0.175 4.175 27 C -0.016 4.016 28 H 0.077 0.923 29 H 0.089 0.911 30 H 0.106 0.894 31 H 0.086 0.914 32 H 0.099 0.901 33 H 0.084 0.916 34 H 0.107 0.893 35 H 0.089 0.911 36 H 0.079 0.921 37 H 0.095 0.905 38 H 0.116 0.884 39 H 0.123 0.877 40 H 0.081 0.919 41 H 0.098 0.902 42 H 0.087 0.913 43 H 0.083 0.917 44 H 0.085 0.915 45 H 0.070 0.930 46 H 0.193 0.807 47 H 0.068 0.932 48 H 0.077 0.923 49 H 0.082 0.918 50 H 0.075 0.925 51 H 0.089 0.911 52 H 0.123 0.877 53 H 0.095 0.905 Dipole moment (debyes) X Y Z Total from point charges -16.898 22.909 3.245 28.652 hybrid contribution 0.599 0.168 0.889 1.086 sum -16.299 23.077 4.135 28.553 Atomic orbital electron populations 1.22508 0.94033 1.05069 1.02371 1.22077 0.94323 0.91498 0.82190 1.21858 1.02341 0.92501 1.03050 1.22522 1.01733 1.03572 0.96331 1.86275 1.71563 1.23609 1.46420 1.17903 0.78777 0.80554 0.81704 1.90895 1.41426 1.70331 1.42865 1.46906 1.53108 1.12139 1.21438 1.22447 0.92636 0.97642 0.91708 1.22969 0.94749 0.83132 1.02662 1.88016 1.73618 1.27673 1.38449 1.21874 0.90075 0.93031 0.88124 1.21436 0.98666 0.91430 1.01918 1.21749 0.95455 0.99302 0.99573 1.20174 0.87509 0.89982 0.94852 1.42796 1.28823 1.61053 1.06413 1.19090 1.04994 1.23427 0.89387 0.89311 1.24674 0.97874 1.00061 0.96753 1.69711 1.20347 1.38471 1.29926 0.86542 0.78097 0.77869 0.86013 1.21146 1.21899 1.20257 1.22133 0.95523 0.97500 1.01745 1.22068 0.98150 1.00545 0.96741 1.22428 0.97463 0.77980 1.03701 0.92301 0.91143 0.89367 0.91408 0.90124 0.91552 0.89278 0.91140 0.92061 0.90533 0.88413 0.87720 0.91912 0.90230 0.91277 0.91686 0.91451 0.92967 0.80691 0.93162 0.92340 0.91801 0.92521 0.91097 0.87714 0.90457 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 C -0.18 -1.60 8.37 37.15 0.31 -1.29 16 2 C 0.14 1.25 1.13 -90.62 -0.10 1.15 16 3 C -0.14 -0.91 8.85 37.16 0.33 -0.58 16 4 C -0.18 -1.77 8.37 37.16 0.31 -1.46 16 5 O -0.36 -4.04 8.58 -40.32 -0.35 -4.39 16 6 C 0.66 8.24 7.75 54.05 0.42 8.65 16 7 O -0.57 -8.17 11.35 -19.28 -0.22 -8.39 16 8 N -0.60 -6.73 2.54 -178.11 -0.45 -7.18 16 9 C 0.08 0.72 6.65 -3.53 -0.02 0.70 16 10 C 0.04 0.44 6.60 37.31 0.25 0.69 16 11 O -0.36 -4.61 8.93 -35.23 -0.31 -4.92 16 12 C 0.11 1.36 0.68 -90.55 -0.06 1.30 16 13 C -0.10 -1.19 4.94 -26.75 -0.13 -1.32 16 14 C -0.12 -1.91 5.25 -26.73 -0.14 -2.05 16 15 C 0.18 3.53 2.77 -67.92 -0.19 3.34 16 16 N -0.75 -17.85 4.92 -53.42 -0.26 -18.11 16 17 C -0.24 -6.78 10.41 -15.06 -0.16 -6.94 16 18 H 0.09 2.70 8.06 -52.49 -0.42 2.28 16 19 C 0.00 0.05 5.50 -17.43 -0.10 -0.04 16 20 N -0.94 -29.71 13.97 55.49 0.78 -28.94 16 21 Si 0.88 22.40 27.73 -169.99 -4.71 17.69 16 22 H -0.28 -7.16 7.11 56.52 0.40 -6.75 16 23 H -0.29 -7.52 7.11 56.52 0.40 -7.12 16 24 H -0.28 -6.62 6.52 56.52 0.37 -6.25 16 25 C -0.13 -2.28 5.36 -26.73 -0.14 -2.42 16 26 C -0.14 -1.96 4.19 -26.76 -0.11 -2.07 16 27 C 0.11 1.19 5.29 -3.57 -0.02 1.17 16 28 H 0.06 0.48 8.14 -51.93 -0.42 0.06 16 29 H 0.07 0.49 8.14 -51.93 -0.42 0.06 16 30 H 0.09 0.99 5.88 -51.93 -0.31 0.68 16 31 H 0.07 0.44 8.14 -51.93 -0.42 0.02 16 32 H 0.08 0.39 8.14 -51.94 -0.42 -0.03 16 33 H 0.07 0.40 8.14 -51.93 -0.42 -0.02 16 34 H 0.09 1.08 5.88 -51.93 -0.31 0.78 16 35 H 0.07 0.54 8.14 -51.93 -0.42 0.11 16 36 H 0.06 0.59 8.14 -51.93 -0.42 0.17 16 37 H 0.08 0.62 8.14 -51.93 -0.42 0.20 16 38 H 0.10 0.80 7.01 -51.93 -0.36 0.44 16 39 H 0.10 0.87 8.14 -51.93 -0.42 0.45 16 40 H 0.06 0.72 6.26 -51.93 -0.33 0.39 16 41 H 0.08 0.91 8.14 -51.93 -0.42 0.49 16 42 H 0.07 0.75 8.14 -51.93 -0.42 0.33 16 43 H 0.06 1.01 8.14 -51.93 -0.42 0.59 16 44 H 0.07 1.10 7.70 -51.93 -0.40 0.70 16 45 H 0.05 1.05 8.08 -51.93 -0.42 0.63 16 46 H 0.36 8.51 5.66 -40.82 -0.23 8.28 16 47 H 0.26 7.82 8.31 -40.82 -0.34 7.48 16 48 H 0.06 1.06 7.70 -51.93 -0.40 0.66 16 49 H 0.06 1.27 8.14 -51.94 -0.42 0.85 16 50 H 0.06 0.81 6.11 -51.92 -0.32 0.49 16 51 H 0.07 1.03 8.14 -51.93 -0.42 0.61 16 52 H 0.10 1.25 7.07 -51.93 -0.37 0.89 16 53 H 0.08 0.77 8.14 -51.93 -0.42 0.35 16 LS Contribution 396.75 15.07 5.98 5.98 Total: -1.00 -33.19 396.75 -8.48 -41.67 By element: Atomic # 1 Polarization: 17.16 SS G_CDS: -9.37 Total: 7.80 kcal Atomic # 6 Polarization: -1.64 SS G_CDS: 0.44 Total: -1.19 kcal Atomic # 7 Polarization: -54.29 SS G_CDS: 0.06 Total: -54.23 kcal Atomic # 8 Polarization: -16.82 SS G_CDS: -0.88 Total: -17.70 kcal Atomic # 14 Polarization: 22.40 SS G_CDS: -4.71 Total: 17.69 kcal Total LS contribution 5.98 Total: 5.98 kcal Total: -33.19 -8.48 -41.67 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. ZINC001254712930.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 -71.163 kcal (2) G-P(sol) polarization free energy of solvation -33.186 kcal (3) G-ENP(sol) elect.-nuc.-pol. free energy of system -104.349 kcal (4) G-CDS(sol) cavity-dispersion-solvent structure free energy -8.482 kcal (5) G-P-CDS(sol) = G-P(sol) + G-CDS(sol) = (2) + (4) -41.668 kcal (6) G-S(sol) free energy of system = (1) + (5) -112.832 kcal 1SCF WAS SPECIFIED, SO GEOMETRY OPTIMIZATION WAS NOT USED Total computer time = 1.11 seconds