Wall clock time and date at job start Sat Apr 11 2020 06:13:02 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.52997 * 1 3 3 C 1.52999 * 109.46995 * 2 1 4 4 H 1.08993 * 109.47672 * 302.43141 * 3 2 1 5 5 C 1.50703 * 109.46956 * 182.42953 * 3 2 1 6 6 O 1.21277 * 120.00134 * 245.45669 * 5 3 2 7 7 N 1.34775 * 119.99773 * 65.46573 * 5 3 2 8 8 C 1.46504 * 120.00146 * 179.97438 * 7 5 3 9 9 H 1.08999 * 109.49747 * 324.91916 * 8 7 5 10 10 C 1.53124 * 109.50420 * 84.99236 * 8 7 5 11 11 C 1.53119 * 109.16057 * 177.63578 * 10 8 7 12 12 C 1.53040 * 109.27714 * 302.35126 * 11 10 8 13 13 N 1.46840 * 109.52594 * 59.17841 * 12 11 10 14 14 C 1.46903 * 111.00216 * 174.15895 * 13 12 11 15 15 C 1.50694 * 109.47319 * 191.26288 * 14 13 12 16 16 O 1.21294 * 120.00460 * 5.07147 * 15 14 13 17 17 N 1.34779 * 120.00014 * 185.07491 * 15 14 13 18 18 C 1.37097 * 120.00183 * 180.02562 * 17 15 14 19 19 H 1.07997 * 119.99790 * 359.97438 * 18 17 15 20 20 C 1.38954 * 120.00049 * 179.97438 * 18 17 15 21 21 N 1.31411 * 120.00050 * 359.97438 * 20 18 17 22 22 Si 1.86801 * 119.99717 * 179.97438 * 20 18 17 23 23 H 1.48505 * 109.46841 * 90.00152 * 22 20 18 24 24 H 1.48499 * 109.46837 * 209.99789 * 22 20 18 25 25 H 1.48493 * 109.47454 * 330.00165 * 22 20 18 26 26 C 1.46848 * 111.19968 * 298.25864 * 13 12 11 27 27 N 1.46902 * 109.47294 * 62.43491 * 3 2 1 28 28 C 1.46898 * 111.00086 * 59.30975 * 27 3 2 29 29 C 1.46900 * 111.00443 * 183.26716 * 27 3 2 30 30 H 1.09001 * 109.47605 * 64.88025 * 1 2 3 31 31 H 1.08999 * 109.47225 * 184.88427 * 1 2 3 32 32 H 1.09000 * 109.47062 * 304.88233 * 1 2 3 33 33 H 1.08999 * 109.46833 * 240.00532 * 2 1 3 34 34 H 1.08997 * 109.47305 * 119.99886 * 2 1 3 35 35 H 0.96999 * 120.00171 * 359.97438 * 7 5 3 36 36 H 1.09004 * 109.51767 * 57.73016 * 10 8 7 37 37 H 1.08993 * 109.57422 * 297.57790 * 10 8 7 38 38 H 1.09006 * 109.49830 * 62.28549 * 11 10 8 39 39 H 1.09000 * 109.50298 * 182.40561 * 11 10 8 40 40 H 1.09002 * 109.45784 * 179.19514 * 12 11 10 41 41 H 1.08997 * 109.45872 * 299.16070 * 12 11 10 42 42 H 1.09002 * 109.46963 * 311.26389 * 14 13 12 43 43 H 1.08996 * 109.46808 * 71.25777 * 14 13 12 44 44 H 0.97004 * 119.99697 * 359.97438 * 17 15 14 45 45 H 0.96994 * 120.00188 * 180.02562 * 21 20 18 46 46 H 1.09002 * 109.45843 * 301.72558 * 26 13 12 47 47 H 1.09000 * 109.45688 * 181.75410 * 26 13 12 48 48 H 1.08998 * 109.47086 * 66.80706 * 28 27 3 49 49 H 1.09001 * 109.47592 * 186.80650 * 28 27 3 50 50 H 1.09001 * 109.47370 * 306.80749 * 28 27 3 51 51 H 1.09005 * 109.46991 * 64.66759 * 29 27 3 52 52 H 1.08997 * 109.47173 * 184.66873 * 29 27 3 53 53 H 1.08997 * 109.46830 * 304.67433 * 29 27 3 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 1 1.6415 1.9688 -0.8673 5 6 3.5458 1.4429 -0.0602 6 8 4.1900 1.8658 0.8762 7 7 4.1765 0.9735 -1.1549 8 6 5.6403 0.9744 -1.2137 9 1 6.0247 1.8536 -0.6968 10 6 6.1798 -0.2908 -0.5407 11 6 7.7063 -0.3109 -0.6587 12 6 8.0998 -0.2460 -2.1362 13 7 7.5609 0.9820 -2.7342 14 6 8.0344 1.1471 -4.1150 15 6 7.6962 2.5332 -4.5998 16 8 7.2009 3.3378 -3.8391 17 7 7.9440 2.8790 -5.8787 18 6 7.6368 4.1403 -6.3196 19 1 7.1958 4.8566 -5.6423 20 6 7.8927 4.4969 -7.6380 21 7 8.4297 3.6255 -8.4620 22 14 7.4735 6.2153 -8.2389 23 1 8.6399 7.1106 -8.0310 24 1 7.1387 6.1660 -9.6848 25 1 6.3088 6.7351 -7.4785 26 6 6.0937 1.0019 -2.6751 27 7 1.5989 2.1178 1.2278 28 6 2.1047 1.4252 2.4204 29 6 2.0049 3.5296 1.2271 30 1 -0.3634 0.4362 -0.9304 31 1 -0.3633 -1.0239 0.0875 32 1 -0.3633 0.5877 0.8430 33 1 1.8932 -0.5138 -0.8900 34 1 1.8933 -0.5138 0.8900 35 1 3.6612 0.6350 -1.9038 36 1 5.7648 -1.1704 -1.0330 37 1 5.8956 -0.2942 0.5115 38 1 8.1232 0.5483 -0.1331 39 1 8.0932 -1.2302 -0.2190 40 1 9.1863 -0.2457 -2.2230 41 1 7.6941 -1.1126 -2.6584 42 1 9.1145 1.0040 -4.1482 43 1 7.5504 0.4099 -4.7555 44 1 8.3404 2.2357 -6.4869 45 1 8.6088 3.8745 -9.3822 46 1 5.6975 0.1296 -3.1948 47 1 5.7233 1.9088 -3.1532 48 1 3.1899 1.5192 2.4622 49 1 1.6668 1.8717 3.3131 50 1 1.8330 0.3708 2.3708 51 1 1.4939 4.0542 0.4196 52 1 1.7392 3.9845 2.1813 53 1 3.0827 3.5976 1.0796 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 ZINC001481538688.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:13:02 Heat of formation + Delta-G solvation = -59.052615 kcal Electronic energy + Delta-G solvation = -31133.938051 eV Core-core repulsion = 27040.707830 eV Total energy + Delta-G solvation = -4093.230221 eV No. of doubly occupied orbitals = 66 Molecular weight (most abundant/longest-lived isotopes) = 340.226 amu Computer time = 0.77 seconds Orbital eigenvalues (eV) -40.71246 -39.24650 -37.79804 -36.73295 -34.88912 -33.95823 -33.23861 -32.46056 -31.23267 -29.91293 -27.75705 -27.52512 -27.33764 -24.71834 -23.57810 -22.97421 -22.89189 -21.40427 -20.23379 -19.65782 -19.17207 -17.88136 -17.44020 -16.74811 -16.04953 -15.84576 -15.37321 -15.24494 -15.11825 -14.79661 -14.61345 -14.38975 -14.17844 -13.95701 -13.83211 -13.45398 -13.41349 -13.27299 -13.07163 -12.60759 -12.50409 -12.17320 -12.02178 -11.97670 -11.83974 -11.70832 -11.53541 -11.48293 -11.32122 -11.21712 -11.06913 -10.94437 -10.88672 -10.64316 -10.54712 -10.40923 -10.26648 -10.22781 -9.49802 -9.04627 -8.47375 -8.24391 -7.99897 -7.90928 -6.43543 -3.84289 2.50362 3.02869 3.16498 3.25316 3.30952 3.75482 3.88753 3.98493 4.23238 4.45743 4.63506 4.76004 4.83821 4.91282 4.95468 5.01071 5.04263 5.16062 5.21605 5.24554 5.29631 5.37257 5.41841 5.49871 5.55580 5.58598 5.59477 5.66670 5.69096 5.76526 5.77791 5.81277 5.82256 5.94305 5.95114 6.00095 6.06228 6.10757 6.12509 6.23212 6.28548 6.38775 6.48081 6.72602 6.82257 7.20253 7.33438 7.43657 7.53562 8.04638 8.15121 8.60089 9.16601 9.52671 10.04575 10.74225 Molecular weight = 340.23amu Principal moments of inertia in cm(-1) A = 0.011726 B = 0.002651 C = 0.002336 Principal moments of inertia in units of 10**(-40)*gram-cm**2 A = 2387.320693 B =10557.839442 C =11981.776322 Net atomic charges, atomic populations, and dipole contributions calculated with CM2 Atom NO. Type Charge No. of electrons 1 C -0.144 4.144 2 C -0.113 4.113 3 C 0.062 3.938 4 H 0.108 0.892 5 C 0.490 3.510 6 O -0.536 6.536 7 N -0.720 5.720 8 C 0.150 3.850 9 H 0.098 0.902 10 C -0.127 4.127 11 C -0.112 4.112 12 C 0.053 3.947 13 N -0.511 5.511 14 C 0.048 3.952 15 C 0.442 3.558 16 O -0.575 6.575 17 N -0.565 5.565 18 C -0.298 4.298 19 H 0.103 0.897 20 C 0.022 3.978 21 N -0.902 5.902 22 Si 0.929 3.071 23 H -0.278 1.278 24 H -0.284 1.284 25 H -0.270 1.270 26 C 0.052 3.948 27 N -0.513 5.513 28 C -0.002 4.002 29 C 0.009 3.991 30 H 0.054 0.946 31 H 0.051 0.949 32 H 0.070 0.930 33 H 0.069 0.931 34 H 0.070 0.930 35 H 0.398 0.602 36 H 0.057 0.943 37 H 0.065 0.935 38 H 0.074 0.926 39 H 0.068 0.932 40 H 0.076 0.924 41 H 0.030 0.970 42 H 0.091 0.909 43 H 0.047 0.953 44 H 0.392 0.608 45 H 0.271 0.729 46 H 0.026 0.974 47 H 0.102 0.898 48 H 0.067 0.933 49 H 0.061 0.939 50 H 0.063 0.937 51 H 0.063 0.937 52 H 0.066 0.934 53 H 0.053 0.947 For the charges calculated by CM2: Dipole moment (debyes) X Y Z Total from point charges -10.787 -12.339 14.657 21.987 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.201 4.201 2 C -0.151 4.151 3 C -0.051 4.051 4 H 0.125 0.875 5 C 0.275 3.725 6 O -0.414 6.414 7 N -0.370 5.370 8 C 0.045 3.955 9 H 0.116 0.884 10 C -0.165 4.165 11 C -0.150 4.150 12 C -0.075 4.075 13 N -0.236 5.236 14 C -0.083 4.083 15 C 0.227 3.773 16 O -0.459 6.459 17 N -0.202 5.202 18 C -0.428 4.428 19 H 0.121 0.879 20 C -0.179 4.179 21 N -0.544 5.544 22 Si 0.756 3.244 23 H -0.204 1.204 24 H -0.209 1.209 25 H -0.194 1.194 26 C -0.076 4.076 27 N -0.236 5.236 28 C -0.150 4.150 29 C -0.139 4.139 30 H 0.073 0.927 31 H 0.070 0.930 32 H 0.089 0.911 33 H 0.087 0.913 34 H 0.088 0.912 35 H 0.231 0.769 36 H 0.076 0.924 37 H 0.083 0.917 38 H 0.093 0.907 39 H 0.086 0.914 40 H 0.095 0.905 41 H 0.048 0.952 42 H 0.110 0.890 43 H 0.065 0.935 44 H 0.225 0.775 45 H 0.081 0.919 46 H 0.044 0.956 47 H 0.120 0.880 48 H 0.086 0.914 49 H 0.080 0.920 50 H 0.082 0.918 51 H 0.082 0.918 52 H 0.085 0.915 53 H 0.072 0.928 Dipole moment (debyes) X Y Z Total from point charges -10.962 -11.065 15.939 22.285 hybrid contribution 0.548 -0.393 -1.801 1.923 sum -10.414 -11.458 14.138 20.967 Atomic orbital electron populations 1.21669 0.94151 1.01420 1.02892 1.21573 0.97681 0.93211 1.02601 1.21993 0.89163 0.97787 0.96122 0.87459 1.21991 0.91413 0.76724 0.82392 1.90603 1.68147 1.46568 1.36037 1.46201 1.07666 1.63949 1.19221 1.20937 0.78619 0.99320 0.96642 0.88401 1.21947 0.96054 0.97525 1.01015 1.21515 0.95162 1.02291 0.96070 1.22166 1.00110 0.90340 0.94900 1.61332 1.08116 1.35801 1.18316 1.21855 1.01188 0.96893 0.88414 1.19755 0.85321 0.89041 0.83185 1.90519 1.48134 1.52312 1.54960 1.41881 1.55309 1.14032 1.08938 1.19690 1.38389 0.89617 0.95079 0.87945 1.25254 0.97112 0.99695 0.95860 1.69608 1.43509 1.35078 1.06178 0.85870 0.79020 0.81781 0.77705 1.20382 1.20940 1.19425 1.22325 0.88209 1.02204 0.94899 1.62492 1.57059 1.02877 1.01122 1.22761 1.01571 0.98639 0.91995 1.22556 1.00965 0.88068 1.02285 0.92732 0.93023 0.91101 0.91251 0.91188 0.76859 0.92397 0.91659 0.90712 0.91384 0.90512 0.95174 0.89046 0.93499 0.77526 0.91878 0.95581 0.87962 0.91428 0.91975 0.91842 0.91805 0.91503 0.92846 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.14 -1.01 9.84 37.16 0.37 -0.65 16 2 C -0.11 -0.87 5.10 -26.73 -0.14 -1.01 16 3 C 0.06 0.61 2.21 -68.86 -0.15 0.46 16 4 H 0.11 0.99 8.14 -51.93 -0.42 0.56 16 5 C 0.49 6.09 6.37 -10.98 -0.07 6.02 16 6 O -0.54 -8.37 11.35 5.56 0.06 -8.30 16 7 N -0.72 -8.40 5.03 -53.73 -0.27 -8.67 16 8 C 0.15 2.08 2.61 -67.85 -0.18 1.90 16 9 H 0.10 1.72 7.58 -51.93 -0.39 1.32 16 10 C -0.13 -1.50 5.58 -26.61 -0.15 -1.64 16 11 C -0.11 -1.28 6.09 -26.65 -0.16 -1.45 16 12 C 0.05 0.66 5.46 -3.84 -0.02 0.64 16 13 N -0.51 -8.29 4.34 -236.12 -1.02 -9.32 16 14 C 0.05 0.83 5.45 -4.98 -0.03 0.80 16 15 C 0.44 10.22 7.34 -10.99 -0.08 10.14 16 16 O -0.58 -15.32 14.27 5.55 0.08 -15.24 16 17 N -0.56 -14.08 5.29 -10.96 -0.06 -14.14 16 18 C -0.30 -8.22 9.68 -14.93 -0.14 -8.37 16 19 H 0.10 2.97 7.16 -52.49 -0.38 2.60 16 20 C 0.02 0.59 5.50 -17.47 -0.10 0.50 16 21 N -0.90 -24.33 13.05 55.50 0.72 -23.60 16 22 Si 0.93 20.94 29.55 -169.99 -5.02 15.92 16 23 H -0.28 -6.22 7.11 56.52 0.40 -5.82 16 24 H -0.28 -6.25 7.11 56.52 0.40 -5.85 16 25 H -0.27 -6.07 7.11 56.52 0.40 -5.67 16 26 C 0.05 0.80 4.59 -3.84 -0.02 0.79 16 27 N -0.51 -5.41 4.28 -244.59 -1.05 -6.45 16 28 C 0.00 -0.02 8.89 60.06 0.53 0.51 16 29 C 0.01 0.10 9.36 60.06 0.56 0.67 16 30 H 0.05 0.37 8.14 -51.93 -0.42 -0.06 16 31 H 0.05 0.32 8.14 -51.93 -0.42 -0.11 16 32 H 0.07 0.57 7.53 -51.93 -0.39 0.18 16 33 H 0.07 0.46 7.89 -51.93 -0.41 0.05 16 34 H 0.07 0.58 6.12 -51.93 -0.32 0.26 16 35 H 0.40 3.70 8.38 -40.82 -0.34 3.36 16 36 H 0.06 0.62 8.14 -51.93 -0.42 0.19 16 37 H 0.06 0.82 8.14 -51.93 -0.42 0.39 16 38 H 0.07 0.96 8.14 -51.93 -0.42 0.54 16 39 H 0.07 0.64 8.14 -51.93 -0.42 0.21 16 40 H 0.08 0.89 8.00 -51.93 -0.42 0.48 16 41 H 0.03 0.35 8.14 -51.93 -0.42 -0.08 16 42 H 0.09 1.44 7.96 -51.93 -0.41 1.02 16 43 H 0.05 0.74 8.11 -51.93 -0.42 0.32 16 44 H 0.39 9.33 7.90 -40.82 -0.32 9.01 16 45 H 0.27 6.93 8.31 -40.82 -0.34 6.59 16 46 H 0.03 0.36 8.14 -51.93 -0.42 -0.07 16 47 H 0.10 1.96 6.07 -51.93 -0.32 1.64 16 48 H 0.07 0.87 5.41 -51.93 -0.28 0.59 16 49 H 0.06 0.58 8.14 -51.93 -0.42 0.16 16 50 H 0.06 0.58 5.89 -51.93 -0.31 0.28 16 51 H 0.06 0.66 8.14 -51.93 -0.42 0.23 16 52 H 0.07 0.67 8.09 -51.93 -0.42 0.25 16 53 H 0.05 0.73 5.76 -51.93 -0.30 0.44 16 LS Contribution 408.29 15.07 6.15 6.15 Total: -1.00 -31.92 408.29 -9.39 -41.31 By element: Atomic # 1 Polarization: 22.24 SS G_CDS: -9.21 Total: 13.03 kcal Atomic # 6 Polarization: 9.08 SS G_CDS: 0.23 Total: 9.31 kcal Atomic # 7 Polarization: -60.50 SS G_CDS: -1.68 Total: -62.17 kcal Atomic # 8 Polarization: -23.68 SS G_CDS: 0.14 Total: -23.54 kcal Atomic # 14 Polarization: 20.94 SS G_CDS: -5.02 Total: 15.92 kcal Total LS contribution 6.15 Total: 6.15 kcal Total: -31.92 -9.39 -41.31 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. ZINC001481538688.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 -17.744 kcal (2) G-P(sol) polarization free energy of solvation -31.923 kcal (3) G-ENP(sol) elect.-nuc.-pol. free energy of system -49.667 kcal (4) G-CDS(sol) cavity-dispersion-solvent structure free energy -9.385 kcal (5) G-P-CDS(sol) = G-P(sol) + G-CDS(sol) = (2) + (4) -41.308 kcal (6) G-S(sol) free energy of system = (1) + (5) -59.053 kcal 1SCF WAS SPECIFIED, SO GEOMETRY OPTIMIZATION WAS NOT USED Total computer time = 0.77 seconds