Wall clock time and date at job start Sat Mar 6 2021 13:13:28 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= 0 * AM1 - THE AM1 HAMILTONIAN TO BE USED * SM5.42R - SM5.42R CALCULATIONS WILL BE PERFORMED * - CHARGE MODEL 2 WILL BE USED (CM2A) * - THE SOLVENT IS WATER ******************************************************************************* 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.52995 * 1 3 3 C 1.52993 * 109.47585 * 2 1 4 4 C 1.53003 * 109.47198 * 66.18412 * 3 2 1 5 5 C 1.53002 * 109.47296 * 186.18176 * 3 2 1 6 6 C 1.50703 * 109.46993 * 60.69721 * 5 3 2 7 7 O 1.21278 * 120.00148 * 237.49950 * 6 5 3 8 8 N 1.34783 * 119.99850 * 57.49411 * 6 5 3 9 9 C 1.46494 * 119.99936 * 185.04605 * 8 6 5 10 10 H 1.09002 * 112.85271 * 336.20430 * 9 8 6 11 11 C 1.53778 * 113.61354 * 107.50606 * 9 8 6 12 12 C 1.53781 * 87.08072 * 89.11858 * 11 9 8 13 13 H 1.09004 * 113.61100 * 270.87874 * 12 11 9 14 14 C 1.50700 * 113.61194 * 139.98499 * 12 11 9 15 15 C 1.38335 * 120.00487 * 142.50367 * 14 12 11 16 16 Cl 1.73596 * 120.01541 * 359.97438 * 15 14 12 17 17 C 1.38373 * 119.97218 * 179.97438 * 15 14 12 18 18 C 1.38199 * 119.98627 * 0.02562 * 17 15 14 19 19 C 1.38253 * 120.01346 * 0.02562 * 18 17 15 20 20 C 1.38230 * 120.02913 * 359.70273 * 19 18 17 21 21 C 1.53780 * 87.08203 * 25.43324 * 12 11 9 22 22 N 1.46494 * 109.47243 * 300.69660 * 5 3 2 23 23 C 1.36344 * 125.63293 * 64.52251 * 22 5 3 24 24 C 1.34855 * 108.21545 * 179.82819 * 23 22 5 25 25 C 1.40896 * 107.41745 * 0.42771 * 24 23 22 26 26 C 1.34848 * 107.41297 * 359.73045 * 25 24 23 27 27 C 1.50702 * 125.88423 * 180.02562 * 26 25 24 28 28 O 1.42908 * 109.46827 * 209.97815 * 27 26 25 29 29 Si 1.86300 * 109.47444 * 329.97367 * 27 26 25 30 30 H 1.08998 * 109.47572 * 180.02562 * 1 2 3 31 31 H 1.09007 * 109.46923 * 300.00263 * 1 2 3 32 32 H 1.08998 * 109.47385 * 59.99405 * 1 2 3 33 33 H 1.09003 * 109.46910 * 120.00041 * 2 1 3 34 34 H 1.09007 * 109.46923 * 239.99737 * 2 1 3 35 35 H 1.08999 * 109.47548 * 306.18541 * 3 2 1 36 36 H 1.08998 * 109.47449 * 186.78606 * 4 3 2 37 37 H 1.08996 * 109.46856 * 306.79298 * 4 3 2 38 38 H 1.09003 * 109.47034 * 66.78599 * 4 3 2 39 39 H 1.08997 * 109.47107 * 180.69594 * 5 3 2 40 40 H 0.96996 * 119.99772 * 5.05297 * 8 6 5 41 41 H 1.09006 * 113.61378 * 203.66845 * 11 9 8 42 42 H 1.08998 * 113.61808 * 334.56517 * 11 9 8 43 43 H 1.07999 * 120.00847 * 179.97438 * 17 15 14 44 44 H 1.07997 * 119.99017 * 179.97438 * 18 17 15 45 45 H 1.08007 * 119.98741 * 179.97438 * 19 18 17 46 46 H 1.07997 * 119.99540 * 180.29768 * 20 19 18 47 47 H 1.09001 * 113.61332 * 220.02044 * 21 12 11 48 48 H 1.08998 * 113.68697 * 89.12544 * 21 12 11 49 49 H 1.07999 * 125.89298 * 0.02562 * 23 22 5 50 50 H 1.08007 * 126.29155 * 180.15842 * 24 23 22 51 51 H 1.07999 * 126.29326 * 179.74140 * 25 24 23 52 52 H 1.08996 * 109.47470 * 89.97777 * 27 26 25 53 53 H 0.96705 * 113.99825 * 59.99560 * 28 27 26 54 54 H 1.48492 * 110.00117 * 180.02562 * 29 27 26 55 55 H 1.48500 * 109.99861 * 301.32675 * 29 27 26 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.0400 1.4424 0.0000 4 6 1.6609 2.1174 1.3197 5 6 3.5622 1.4451 -0.1553 6 6 3.9348 0.7955 -1.4631 7 8 4.6734 -0.1663 -1.4737 8 7 3.4470 1.2835 -2.6210 9 6 3.7186 0.5877 -3.8812 10 1 4.6162 -0.0293 -3.8392 11 6 2.4940 -0.1416 -4.4585 12 6 2.1585 1.1677 -5.1921 13 1 1.5298 1.8439 -4.6127 14 6 1.6570 0.9605 -6.5980 15 6 0.6543 1.7699 -7.1012 16 17 -0.0270 3.0242 -6.1132 17 6 0.1948 1.5794 -8.3924 18 6 0.7372 0.5803 -9.1782 19 6 1.7384 -0.2290 -8.6744 20 6 2.2019 -0.0355 -7.3866 21 6 3.6574 1.5027 -5.1157 22 7 4.1672 0.6987 0.9505 23 6 4.0085 -0.6341 1.1901 24 6 4.7203 -0.9523 2.2903 25 6 5.3333 0.2292 2.7524 26 6 4.9774 1.2231 1.9135 27 6 5.4025 2.6649 2.0216 28 8 4.3801 3.5047 1.4815 29 14 5.6803 3.0980 3.8121 30 1 -0.3634 -1.0276 0.0005 31 1 -0.3633 0.5139 0.8900 32 1 -0.3634 0.5139 -0.8899 33 1 1.8933 -0.5139 0.8900 34 1 1.8933 -0.5139 -0.8900 35 1 1.5892 1.9873 -0.8294 36 1 1.9164 3.1760 1.2729 37 1 0.5896 2.0093 1.4884 38 1 2.2068 1.6478 2.1380 39 1 3.9268 2.4722 -0.1429 40 1 2.9163 2.0954 -2.6198 41 1 2.7504 -0.9634 -5.1272 42 1 1.7573 -0.4214 -3.7055 43 1 -0.5883 2.2107 -8.7857 44 1 0.3782 0.4313 -10.1858 45 1 2.1618 -1.0099 -9.2888 46 1 2.9874 -0.6650 -6.9952 47 1 3.8598 2.5504 -4.8932 48 1 4.2328 1.1311 -5.9636 49 1 3.4174 -1.3177 0.5987 50 1 4.8018 -1.9320 2.7378 51 1 5.9710 0.3217 3.6191 52 1 6.3272 2.8139 1.4641 53 1 3.5272 3.4241 1.9302 54 1 6.0993 4.5175 3.9322 55 1 4.4357 2.8687 4.5891 There are 72 doubly filled levels RHF calculation, no. of doubly occupied orbitals= 72 No. of singly occupied orbitals= 1 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). Cl: (AM1): M.J.S. DEWAR AND E. G. ZOEBISCH, THEOCHEM, 180, 1 (1988). 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=0 AM1 1SCF TLIMIT=15 GEO-OK SM5.42R & SOLVNT=WATER s_22_10565092_12360566.mol2 55 1SCF WAS SPECIFIED, SO GEOMETRY OPTIMIZATION WAS NOT USED SCF FIELD WAS ACHIEVED Wall clock time and date at job start Sat Mar 6 2021 13:13:28 Heat of formation + Delta-G solvation = 42.525911 kcal Electronic energy + Delta-G solvation = -37583.951984 eV Core-core repulsion = 33001.549015 eV Total energy + Delta-G solvation = -4582.402970 eV No. of doubly occupied orbitals = 72 and no. of open levels = 1 Molecular weight (most abundant/longest-lived isotopes) = 403.170 amu Computer time = 1.51 seconds Orbital eigenvalues (eV) -42.92080 -41.58615 -40.41174 -38.94294 -38.05717 -37.65277 -37.43589 -34.69475 -33.36943 -31.83437 -31.69244 -31.09007 -30.54644 -28.43171 -26.78526 -26.28022 -26.10001 -24.85912 -23.80369 -22.87223 -22.41840 -21.61115 -21.25134 -20.49727 -20.17544 -18.95162 -18.34632 -18.00115 -17.39824 -17.02774 -16.78627 -16.52374 -16.41285 -16.35955 -15.99366 -15.66161 -15.36549 -15.28702 -14.97890 -14.89337 -14.78970 -14.60459 -14.41537 -14.27244 -14.14108 -14.03672 -13.65500 -13.50064 -13.38269 -13.37304 -13.17115 -12.99582 -12.96903 -12.72320 -12.60234 -12.50376 -12.36954 -12.18960 -12.01443 -11.90338 -11.83010 -11.76800 -11.60682 -11.55234 -11.31271 -11.10738 -10.78192 -10.21409 -9.86141 -9.60922 -9.51305 -8.85712 -4.86970 0.15600 0.26122 1.06592 1.14305 1.28170 1.34357 1.67862 1.79434 1.98003 2.12240 2.22183 2.60321 3.10572 3.30789 3.36026 3.47028 3.55590 3.60856 3.63395 3.76065 3.82338 3.90458 3.91496 4.07374 4.11984 4.13606 4.23565 4.34965 4.36113 4.37532 4.39244 4.49963 4.53789 4.57249 4.58853 4.68169 4.71703 4.72211 4.73514 4.77957 4.85444 4.86679 4.93201 4.94353 4.97454 4.99308 5.06128 5.07207 5.26135 5.27289 5.34372 5.36098 5.36884 5.42042 5.46851 5.53844 5.84052 6.02331 6.23388 6.25211 6.35429 6.46223 7.00352 Net atomic charges, atomic populations, and dipole contributions calculated with CM2 Atom NO. Type Charge No. of electrons 1 C -0.148 4.148 2 C -0.131 4.131 3 C -0.087 4.087 4 C -0.155 4.155 5 C 0.107 3.893 6 C 0.511 3.489 7 O -0.559 6.559 8 N -0.695 5.695 9 C 0.119 3.881 10 H 0.089 0.911 11 C -0.135 4.135 12 C -0.073 4.073 13 H 0.105 0.895 14 C -0.070 4.070 15 C -0.050 4.050 16 Cl -0.074 7.074 17 C -0.121 4.121 18 C -0.117 4.117 19 C -0.116 4.116 20 C -0.112 4.112 21 C -0.124 4.124 22 N -0.413 5.413 23 C 0.008 3.992 24 C -0.222 4.222 25 C -0.194 4.194 26 C -0.030 4.030 27 C -0.059 4.059 28 O -0.545 6.545 29 Si 0.793 3.207 30 H 0.052 0.948 31 H 0.063 0.937 32 H 0.060 0.940 33 H 0.069 0.931 34 H 0.049 0.951 35 H 0.110 0.890 36 H 0.074 0.926 37 H 0.078 0.922 38 H 0.048 0.952 39 H 0.153 0.847 40 H 0.414 0.586 41 H 0.102 0.898 42 H 0.069 0.931 43 H 0.139 0.861 44 H 0.149 0.851 45 H 0.148 0.852 46 H 0.146 0.854 47 H 0.083 0.917 48 H 0.105 0.895 49 H 0.156 0.844 50 H 0.142 0.858 51 H 0.164 0.836 52 H 0.123 0.877 53 H 0.381 0.619 54 H -0.291 1.291 55 H -0.287 1.287 For the charges calculated by CM2: Dipole moment (debyes) X Y Z Total from point charges -1.979 1.106 -2.209 3.165 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.205 4.205 2 C -0.168 4.168 3 C -0.106 4.106 4 C -0.212 4.212 5 C 0.001 3.999 6 C 0.297 3.703 7 O -0.436 6.436 8 N -0.348 5.348 9 C 0.013 3.987 10 H 0.107 0.893 11 C -0.174 4.174 12 C -0.092 4.092 13 H 0.123 0.877 14 C -0.070 4.070 15 C -0.076 4.076 16 Cl -0.045 7.045 17 C -0.140 4.140 18 C -0.135 4.135 19 C -0.134 4.134 20 C -0.130 4.130 21 C -0.162 4.162 22 N -0.089 5.089 23 C -0.122 4.122 24 C -0.246 4.246 25 C -0.218 4.218 26 C -0.140 4.140 27 C -0.167 4.167 28 O -0.355 6.355 29 Si 0.695 3.305 30 H 0.071 0.929 31 H 0.082 0.918 32 H 0.079 0.921 33 H 0.088 0.912 34 H 0.068 0.932 35 H 0.128 0.872 36 H 0.093 0.907 37 H 0.097 0.903 38 H 0.067 0.933 39 H 0.170 0.830 40 H 0.251 0.749 41 H 0.120 0.880 42 H 0.088 0.912 43 H 0.157 0.843 44 H 0.166 0.834 45 H 0.166 0.834 46 H 0.164 0.836 47 H 0.102 0.898 48 H 0.123 0.877 49 H 0.173 0.827 50 H 0.160 0.840 51 H 0.182 0.818 52 H 0.141 0.859 53 H 0.232 0.768 54 H -0.218 1.218 55 H -0.213 1.213 Dipole moment (debyes) X Y Z Total from point charges -1.441 1.393 -2.143 2.934 hybrid contribution -1.475 0.356 -0.628 1.642 sum -2.916 1.749 -2.771 4.386 Atomic orbital electron populations 1.21772 0.94279 1.01999 1.02461 1.21898 0.97303 0.95740 1.01872 1.21502 0.89785 0.98032 1.01261 1.21992 1.03077 1.00534 0.95605 1.21874 0.96301 0.98559 0.83166 1.20855 0.79462 0.83756 0.86244 1.90772 1.36714 1.28742 1.87422 1.45489 1.54530 1.29479 1.05252 1.21823 0.98790 0.95249 0.82808 0.89332 1.23529 0.97507 0.96809 0.99541 1.21662 0.96824 0.98928 0.91740 0.87659 1.19691 0.95241 0.95446 0.96670 1.21033 0.98920 0.94256 0.93396 1.98360 1.84221 1.51744 1.70216 1.21303 1.00910 0.99450 0.92343 1.21442 0.96651 0.94615 1.00810 1.21593 0.96461 0.99552 0.95800 1.21248 1.00341 0.97338 0.94061 1.23371 0.94856 1.00767 0.97190 1.43488 1.37810 1.07386 1.20242 1.21872 1.05323 0.84747 1.00270 1.21085 1.04235 0.99542 0.99725 1.20889 1.06245 0.92012 1.02696 1.20965 1.02327 0.96638 0.94112 1.24978 0.93393 0.84223 1.14148 1.86504 1.23727 1.65331 1.59964 0.93210 0.85610 0.79872 0.71773 0.92912 0.91839 0.92055 0.91244 0.93202 0.87171 0.90709 0.90278 0.93326 0.82999 0.74905 0.87986 0.91194 0.84343 0.83359 0.83436 0.83601 0.89842 0.87671 0.82682 0.84017 0.81830 0.85912 0.76759 1.21828 1.21332 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.15 0.43 9.70 71.98 0.70 1.13 16 2 C -0.13 0.15 3.90 30.59 0.12 0.27 16 3 C -0.09 0.34 2.14 -10.79 -0.02 0.32 16 4 C -0.16 0.71 7.82 71.98 0.56 1.27 16 5 C 0.11 0.03 2.24 44.25 0.10 0.13 16 6 C 0.51 1.82 5.77 87.66 0.51 2.33 16 7 O -0.56 -6.96 14.90 -3.02 -0.05 -7.01 16 8 N -0.70 1.60 5.09 -440.04 -2.24 -0.64 16 9 C 0.12 -0.21 4.36 45.52 0.20 -0.01 16 10 H 0.09 0.25 7.42 -2.39 -0.02 0.23 16 11 C -0.14 0.41 7.25 31.12 0.23 0.64 16 12 C -0.07 0.26 3.64 -11.01 -0.04 0.22 16 13 H 0.11 -0.38 6.58 -2.39 -0.02 -0.39 16 14 C -0.07 0.07 5.29 -19.84 -0.10 -0.04 16 15 C -0.05 -0.05 6.31 22.33 0.14 0.09 16 16 Cl -0.07 -0.07 27.16 -2.72 -0.07 -0.15 16 17 C -0.12 0.01 9.83 22.29 0.22 0.23 16 18 C -0.12 0.27 10.05 22.26 0.22 0.50 16 19 C -0.12 0.36 10.05 22.27 0.22 0.59 16 20 C -0.11 0.30 8.60 22.27 0.19 0.49 16 21 C -0.12 0.55 7.29 31.12 0.23 0.78 16 22 N -0.41 -1.79 1.63 -491.52 -0.80 -2.59 16 23 C 0.01 0.06 7.49 83.21 0.62 0.68 16 24 C -0.22 -1.47 10.86 22.14 0.24 -1.23 16 25 C -0.19 -0.80 9.43 22.14 0.21 -0.59 16 26 C -0.03 -0.11 5.91 41.08 0.24 0.13 16 27 C -0.06 -0.02 2.94 71.24 0.21 0.19 16 28 O -0.54 0.01 9.86 -148.98 -1.47 -1.46 16 29 Si 0.79 -1.09 31.14 68.60 2.14 1.05 16 30 H 0.05 -0.08 8.14 -2.39 -0.02 -0.10 16 31 H 0.06 -0.22 6.47 -2.38 -0.02 -0.24 16 32 H 0.06 -0.25 8.14 -2.39 -0.02 -0.27 16 33 H 0.07 0.12 5.18 -2.39 -0.01 0.11 16 34 H 0.05 0.08 7.00 -2.38 -0.02 0.06 16 35 H 0.11 -0.87 7.52 -2.39 -0.02 -0.88 16 36 H 0.07 -0.54 6.47 -2.39 -0.02 -0.56 16 37 H 0.08 -0.43 6.32 -2.39 -0.02 -0.44 16 38 H 0.05 -0.03 7.10 -2.39 -0.02 -0.05 16 39 H 0.15 -0.12 5.37 -2.39 -0.01 -0.13 16 40 H 0.41 -2.98 8.11 -92.71 -0.75 -3.73 16 41 H 0.10 -0.42 6.32 -2.38 -0.02 -0.43 16 42 H 0.07 -0.10 8.14 -2.39 -0.02 -0.12 16 43 H 0.14 -0.30 8.06 -2.91 -0.02 -0.32 16 44 H 0.15 -0.94 8.06 -2.91 -0.02 -0.96 16 45 H 0.15 -1.06 8.06 -2.91 -0.02 -1.09 16 46 H 0.15 -0.89 6.67 -2.91 -0.02 -0.91 16 47 H 0.08 -0.42 8.14 -2.39 -0.02 -0.44 16 48 H 0.10 -0.55 7.72 -2.39 -0.02 -0.56 16 49 H 0.16 1.07 6.32 -2.91 -0.02 1.05 16 50 H 0.14 0.61 8.06 -2.91 -0.02 0.58 16 51 H 0.16 -0.22 6.89 -2.91 -0.02 -0.24 16 52 H 0.12 -0.06 7.96 -2.39 -0.02 -0.08 16 53 H 0.38 -1.26 5.34 -74.05 -0.40 -1.66 16 54 H -0.29 -1.30 7.11 99.48 0.71 -0.60 16 55 H -0.29 -1.05 7.11 99.48 0.71 -0.35 16 Total: 0.00 -17.50 430.46 2.33 -15.18 By element: Atomic # 1 Polarization: -12.33 SS G_CDS: -0.17 Total: -12.50 kcal Atomic # 6 Polarization: 3.13 SS G_CDS: 4.99 Total: 8.12 kcal Atomic # 7 Polarization: -0.19 SS G_CDS: -3.04 Total: -3.23 kcal Atomic # 8 Polarization: -6.95 SS G_CDS: -1.51 Total: -8.47 kcal Atomic # 14 Polarization: -1.09 SS G_CDS: 2.14 Total: 1.05 kcal Atomic # 17 Polarization: -0.07 SS G_CDS: -0.07 Total: -0.15 kcal Total: -17.50 2.33 -15.18 kcal The number of atoms in the molecule is 55 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. s_22_10565092_12360566.mol2 55 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 57.703 kcal (2) G-P(sol) polarization free energy of solvation -17.505 kcal (3) G-ENP(sol) elect.-nuc.-pol. free energy of system 40.198 kcal (4) G-CDS(sol) cavity-dispersion-solvent structure free energy 2.328 kcal (5) G-P-CDS(sol) = G-P(sol) + G-CDS(sol) = (2) + (4) -15.177 kcal (6) G-S(sol) free energy of system = (1) + (5) 42.526 kcal 1SCF WAS SPECIFIED, SO GEOMETRY OPTIMIZATION WAS NOT USED Total computer time = 1.51 seconds