Wall clock time and date at job start Sat Apr 11 2020 06:12: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.53003 * 1 3 3 H 1.09003 * 109.46818 * 2 1 4 4 C 1.50703 * 109.47060 * 240.00579 * 2 1 3 5 5 H 1.07997 * 119.99516 * 275.00246 * 4 2 1 6 6 C 1.37879 * 120.00030 * 94.99692 * 4 2 1 7 7 N 1.31516 * 119.99935 * 184.91508 * 6 4 2 8 8 Si 1.86802 * 119.99850 * 4.92125 * 6 4 2 9 9 H 1.48493 * 109.47228 * 84.82826 * 8 6 4 10 10 H 1.48502 * 109.46526 * 204.83161 * 8 6 4 11 11 H 1.48503 * 109.46888 * 324.82472 * 8 6 4 12 12 N 1.46496 * 109.47233 * 120.00292 * 2 1 3 13 13 C 1.34773 * 120.00377 * 275.00096 * 12 2 1 14 14 O 1.21279 * 119.99694 * 0.02562 * 13 12 2 15 15 C 1.50701 * 120.00322 * 180.02562 * 13 12 2 16 16 C 1.52999 * 109.49265 * 299.93381 * 15 13 12 17 17 C 1.53044 * 109.52346 * 179.85378 * 15 13 12 18 18 C 1.53042 * 109.53664 * 178.66416 * 17 15 13 19 19 C 1.52758 * 109.37933 * 61.23950 * 18 17 15 20 20 N 1.46928 * 109.01299 * 305.58598 * 19 18 17 21 21 C 1.34776 * 120.63465 * 233.58654 * 20 19 18 22 22 O 1.21506 * 120.00387 * 179.72046 * 21 20 19 23 23 O 1.34643 * 119.99598 * 359.71861 * 21 20 19 24 24 C 1.45201 * 116.99793 * 179.97438 * 23 21 20 25 25 C 1.52999 * 109.46782 * 60.00242 * 24 23 21 26 26 C 1.52994 * 109.47084 * 179.97438 * 24 23 21 27 27 C 1.52997 * 109.47413 * 300.00451 * 24 23 21 28 28 C 1.46922 * 118.73670 * 53.30854 * 20 19 18 29 29 H 1.09000 * 109.47251 * 55.03656 * 1 2 3 30 30 H 1.09007 * 109.46797 * 175.03194 * 1 2 3 31 31 H 1.08997 * 109.47206 * 295.02713 * 1 2 3 32 32 H 0.96997 * 119.99776 * 180.02562 * 7 6 4 33 33 H 0.97002 * 119.99853 * 95.00417 * 12 2 1 34 34 H 1.09009 * 109.46897 * 180.02562 * 16 15 13 35 35 H 1.09003 * 109.47236 * 300.06133 * 16 15 13 36 36 H 1.08993 * 109.47229 * 60.06758 * 16 15 13 37 37 H 1.08991 * 109.45891 * 298.68704 * 17 15 13 38 38 H 1.09000 * 109.45464 * 58.64468 * 17 15 13 39 39 H 1.08998 * 109.49532 * 181.31750 * 18 17 15 40 40 H 1.09006 * 109.49690 * 301.38386 * 18 17 15 41 41 H 1.09006 * 109.51059 * 185.70339 * 19 18 17 42 42 H 1.09005 * 109.43404 * 65.42577 * 19 18 17 43 43 H 1.08996 * 109.47506 * 59.99904 * 25 24 23 44 44 H 1.09005 * 109.46962 * 179.97438 * 25 24 23 45 45 H 1.09006 * 109.47147 * 299.99689 * 25 24 23 46 46 H 1.09001 * 109.47441 * 60.00106 * 26 24 23 47 47 H 1.08996 * 109.47274 * 180.02562 * 26 24 23 48 48 H 1.09003 * 109.47577 * 300.00432 * 26 24 23 49 49 H 1.08994 * 109.47162 * 59.99371 * 27 24 23 50 50 H 1.08996 * 109.47159 * 180.02562 * 27 24 23 51 51 H 1.09007 * 109.47003 * 300.00030 * 27 24 23 52 52 H 1.09007 * 109.58528 * 66.43118 * 28 20 19 53 53 H 1.08997 * 109.58753 * 186.71090 * 28 20 19 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 1 1.8933 1.0277 0.0000 4 6 2.0324 -0.7103 -1.2306 5 1 2.1354 -1.7853 -1.2292 6 6 2.3602 0.0123 -2.3581 7 7 2.7059 -0.6129 -3.4623 8 14 2.3135 1.8793 -2.3173 9 1 3.5991 2.3950 -1.7820 10 1 2.1057 2.3994 -3.6926 11 1 1.1986 2.3276 -1.4448 12 7 2.0184 -0.6906 1.1961 13 6 2.1471 -0.0184 2.3571 14 8 1.8582 1.1582 2.4118 15 6 2.6500 -0.7288 3.5874 16 6 1.7060 -1.8822 3.9327 17 6 2.7053 0.2544 4.7590 18 6 3.1843 -0.4745 6.0165 19 6 4.5896 -1.0243 5.7795 20 7 4.5807 -1.8543 4.5671 21 6 5.0412 -3.1206 4.5984 22 8 5.0237 -3.7985 3.5901 23 8 5.5206 -3.6344 5.7468 24 6 5.9992 -5.0047 5.7068 25 6 7.1436 -5.1157 4.6974 26 6 6.5021 -5.4099 7.0937 27 6 4.8570 -5.9319 5.2868 28 6 4.0545 -1.2797 3.3215 29 1 -0.3634 0.5889 -0.8422 30 1 -0.3633 -1.0239 -0.0890 31 1 -0.3633 0.4347 0.9311 32 1 2.9369 -0.1045 -4.2555 33 1 2.2491 -1.6318 1.1524 34 1 2.0698 -2.3955 4.8229 35 1 0.7068 -1.4897 4.1218 36 1 1.6692 -2.5836 3.0992 37 1 3.3969 1.0628 4.5222 38 1 1.7116 0.6664 4.9343 39 1 3.2011 0.2208 6.8558 40 1 2.5064 -1.2983 6.2400 41 1 4.8929 -1.6291 6.6341 42 1 5.2862 -0.1959 5.6500 43 1 7.9574 -4.4552 4.9965 44 1 7.5033 -6.1443 4.6676 45 1 6.7858 -4.8268 3.7091 46 1 5.6872 -5.3307 7.8132 47 1 6.8618 -6.4384 7.0636 48 1 7.3160 -4.7493 7.3925 49 1 4.4992 -5.6429 4.2986 50 1 5.2162 -6.9605 5.2572 51 1 4.0416 -5.8527 6.0059 52 1 4.7064 -0.4725 2.9875 53 1 4.0071 -2.0523 2.5541 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 ZINC000179275321.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:02 Heat of formation + Delta-G solvation = -131.137849 kcal Electronic energy + Delta-G solvation = -32176.543311 eV Core-core repulsion = 28047.894489 eV Total energy + Delta-G solvation = -4128.648823 eV No. of doubly occupied orbitals = 66 Molecular weight (most abundant/longest-lived isotopes) = 340.220 amu Computer time = 0.78 seconds Orbital eigenvalues (eV) -41.16362 -39.98478 -38.12325 -36.63057 -34.92858 -34.44326 -33.27342 -32.48863 -31.81263 -29.45538 -27.64063 -27.55377 -27.25538 -25.72327 -25.30525 -23.94273 -22.32299 -21.45897 -19.97555 -19.45046 -18.72815 -18.10430 -17.18942 -16.78440 -16.31602 -15.90480 -15.59306 -15.32280 -14.97580 -14.72628 -14.59772 -14.26505 -13.91859 -13.72057 -13.51494 -13.46330 -13.30267 -12.83240 -12.68815 -12.50339 -12.46616 -12.39614 -12.19567 -12.09988 -12.08063 -11.85399 -11.61037 -11.51556 -11.44754 -11.22679 -11.11424 -10.97806 -10.94582 -10.80453 -10.57724 -10.51735 -10.50949 -10.06464 -9.86938 -9.18979 -9.10412 -8.79830 -8.68104 -8.08770 -6.21653 -4.25295 1.98661 2.76706 3.22789 3.27923 3.34595 3.41483 3.64546 3.93790 4.33445 4.50301 4.54348 4.56619 4.69195 4.71878 4.79816 4.88837 5.03138 5.14924 5.21970 5.24124 5.27783 5.31723 5.40244 5.46960 5.49552 5.54543 5.62087 5.66011 5.66561 5.69397 5.70841 5.76035 5.89372 5.92677 6.06228 6.12304 6.13765 6.20653 6.30434 6.33350 6.36625 6.39253 6.52919 6.68774 6.93809 7.19819 7.27349 7.38115 7.56063 8.14217 8.26715 8.57773 8.79393 8.95843 9.40123 10.90136 Molecular weight = 340.22amu Principal moments of inertia in cm(-1) A = 0.014449 B = 0.003156 C = 0.002819 Principal moments of inertia in units of 10**(-40)*gram-cm**2 A = 1937.371550 B = 8869.878901 C = 9929.467734 Net atomic charges, atomic populations, and dipole contributions calculated with CM2 Atom NO. Type Charge No. of electrons 1 C -0.147 4.147 2 C 0.231 3.769 3 H 0.050 0.950 4 C -0.530 4.530 5 H 0.079 0.921 6 C 0.081 3.919 7 N -0.888 5.888 8 Si 0.874 3.126 9 H -0.273 1.273 10 H -0.286 1.286 11 H -0.261 1.261 12 N -0.699 5.699 13 C 0.513 3.487 14 O -0.555 6.555 15 C -0.095 4.095 16 C -0.136 4.136 17 C -0.094 4.094 18 C -0.137 4.137 19 C 0.111 3.889 20 N -0.590 5.590 21 C 0.653 3.347 22 O -0.562 6.562 23 O -0.364 6.364 24 C 0.136 3.864 25 C -0.184 4.184 26 C -0.139 4.139 27 C -0.184 4.184 28 C 0.140 3.860 29 H 0.052 0.948 30 H 0.047 0.953 31 H 0.028 0.972 32 H 0.267 0.733 33 H 0.394 0.606 34 H 0.063 0.937 35 H 0.061 0.939 36 H 0.063 0.937 37 H 0.072 0.928 38 H 0.075 0.925 39 H 0.079 0.921 40 H 0.071 0.929 41 H 0.091 0.909 42 H 0.070 0.930 43 H 0.058 0.942 44 H 0.069 0.931 45 H 0.090 0.910 46 H 0.065 0.935 47 H 0.079 0.921 48 H 0.065 0.935 49 H 0.090 0.910 50 H 0.069 0.931 51 H 0.058 0.942 52 H 0.078 0.922 53 H 0.105 0.895 For the charges calculated by CM2: Dipole moment (debyes) X Y Z Total from point charges 5.812 -7.224 23.532 25.292 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.206 4.206 2 C 0.128 3.872 3 H 0.068 0.932 4 C -0.568 4.568 5 H 0.097 0.903 6 C -0.123 4.123 7 N -0.524 5.524 8 Si 0.701 3.299 9 H -0.199 1.199 10 H -0.212 1.212 11 H -0.185 1.185 12 N -0.347 5.347 13 C 0.297 3.703 14 O -0.434 6.434 15 C -0.099 4.099 16 C -0.193 4.193 17 C -0.132 4.132 18 C -0.175 4.175 19 C -0.012 4.012 20 N -0.330 5.330 21 C 0.407 3.593 22 O -0.451 6.451 23 O -0.279 6.279 24 C 0.099 3.901 25 C -0.241 4.241 26 C -0.197 4.197 27 C -0.241 4.241 28 C 0.018 3.982 29 H 0.071 0.929 30 H 0.066 0.934 31 H 0.047 0.953 32 H 0.077 0.923 33 H 0.227 0.773 34 H 0.082 0.918 35 H 0.080 0.920 36 H 0.082 0.918 37 H 0.091 0.909 38 H 0.093 0.907 39 H 0.097 0.903 40 H 0.089 0.911 41 H 0.109 0.891 42 H 0.089 0.911 43 H 0.077 0.923 44 H 0.088 0.912 45 H 0.109 0.891 46 H 0.084 0.916 47 H 0.098 0.902 48 H 0.084 0.916 49 H 0.109 0.891 50 H 0.088 0.912 51 H 0.077 0.923 52 H 0.096 0.904 53 H 0.124 0.876 Dipole moment (debyes) X Y Z Total from point charges 5.696 -6.655 23.038 24.647 hybrid contribution 0.150 0.463 -0.208 0.530 sum 5.846 -6.192 22.830 24.367 Atomic orbital electron populations 1.22214 0.98614 1.00570 0.99182 1.19136 0.90917 0.94287 0.82828 0.93244 1.21017 1.45179 0.95686 0.94913 0.90276 1.24942 0.93650 0.98562 0.95167 1.69886 1.44715 1.36152 1.01687 0.86900 0.79753 0.85947 0.77323 1.19872 1.21224 1.18507 1.46700 1.64339 1.13958 1.09735 1.20151 0.78626 0.87192 0.84326 1.90639 1.51779 1.15962 1.85003 1.21207 0.96034 0.97225 0.95468 1.21826 0.98567 0.96658 1.02270 1.21342 1.03040 0.96757 0.92069 1.21953 0.96382 1.00735 0.98475 1.21671 0.97714 0.95029 0.86775 1.46886 1.63969 1.10122 1.12033 1.17935 0.77135 0.82206 0.81980 1.90926 1.60580 1.59127 1.34461 1.86273 1.75005 1.31713 1.34924 1.22050 0.94130 0.77062 0.96860 1.22524 0.97427 1.02641 1.01552 1.21848 1.01826 1.01798 0.94195 1.22516 0.97893 0.99278 1.04395 1.21659 0.92399 0.99155 0.84980 0.92946 0.93414 0.95277 0.92260 0.77252 0.91795 0.91983 0.91826 0.90924 0.90658 0.90285 0.91052 0.89109 0.91131 0.92281 0.91223 0.89147 0.91563 0.90222 0.91570 0.89144 0.91232 0.92281 0.90362 0.87640 Number of geometries 1 Number of calculations of the screened coulomb radii 1 The total number of SCF iterations 12. 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 -3.23 9.47 37.16 0.35 -2.88 16 2 C 0.23 5.49 2.22 -69.08 -0.15 5.33 16 3 H 0.05 1.28 4.73 -51.93 -0.25 1.03 16 4 C -0.53 -14.30 9.05 -34.44 -0.31 -14.62 16 5 H 0.08 2.26 8.06 -52.49 -0.42 1.84 16 6 C 0.08 2.28 5.47 -17.82 -0.10 2.18 16 7 N -0.89 -26.31 13.96 55.43 0.77 -25.54 16 8 Si 0.87 21.65 26.99 -169.99 -4.59 17.06 16 9 H -0.27 -6.77 7.11 56.52 0.40 -6.37 16 10 H -0.29 -6.95 7.11 56.52 0.40 -6.55 16 11 H -0.26 -6.35 6.82 56.52 0.39 -5.97 16 12 N -0.70 -14.03 4.87 -54.93 -0.27 -14.30 16 13 C 0.51 9.59 6.45 -10.99 -0.07 9.52 16 14 O -0.55 -12.25 15.65 5.56 0.09 -12.16 16 15 C -0.10 -1.33 0.62 -155.54 -0.10 -1.42 16 16 C -0.14 -1.70 7.88 37.16 0.29 -1.41 16 17 C -0.09 -1.15 4.01 -26.69 -0.11 -1.25 16 18 C -0.14 -1.23 5.72 -26.83 -0.15 -1.39 16 19 C 0.11 1.03 6.42 -3.93 -0.03 1.01 16 20 N -0.59 -7.15 2.42 -178.82 -0.43 -7.58 16 21 C 0.65 8.72 7.74 54.05 0.42 9.14 16 22 O -0.56 -8.88 11.29 -19.28 -0.22 -9.10 16 23 O -0.36 -4.17 8.53 -40.33 -0.34 -4.51 16 24 C 0.14 1.31 1.13 -90.63 -0.10 1.21 16 25 C -0.18 -1.83 8.37 37.16 0.31 -1.52 16 26 C -0.14 -0.94 8.85 37.15 0.33 -0.61 16 27 C -0.18 -1.79 8.37 37.16 0.31 -1.48 16 28 C 0.14 1.93 4.97 -3.71 -0.02 1.92 16 29 H 0.05 1.16 8.14 -51.93 -0.42 0.74 16 30 H 0.05 0.99 8.14 -51.92 -0.42 0.57 16 31 H 0.03 0.60 8.14 -51.93 -0.42 0.18 16 32 H 0.27 7.45 8.31 -40.82 -0.34 7.11 16 33 H 0.39 7.26 7.56 -40.82 -0.31 6.95 16 34 H 0.06 0.73 5.97 -51.92 -0.31 0.42 16 35 H 0.06 0.77 8.14 -51.93 -0.42 0.34 16 36 H 0.06 0.83 7.57 -51.93 -0.39 0.44 16 37 H 0.07 0.94 8.13 -51.93 -0.42 0.52 16 38 H 0.07 0.97 8.07 -51.93 -0.42 0.55 16 39 H 0.08 0.56 8.14 -51.93 -0.42 0.14 16 40 H 0.07 0.65 6.33 -51.93 -0.33 0.32 16 41 H 0.09 0.76 6.93 -51.93 -0.36 0.40 16 42 H 0.07 0.61 8.14 -51.93 -0.42 0.19 16 43 H 0.06 0.56 8.14 -51.93 -0.42 0.13 16 44 H 0.07 0.56 8.14 -51.93 -0.42 0.14 16 45 H 0.09 1.16 5.88 -51.93 -0.31 0.85 16 46 H 0.07 0.42 8.14 -51.93 -0.42 0.00 16 47 H 0.08 0.41 8.14 -51.93 -0.42 -0.01 16 48 H 0.07 0.42 8.14 -51.93 -0.42 0.00 16 49 H 0.09 1.13 5.88 -51.93 -0.31 0.83 16 50 H 0.07 0.55 8.14 -51.93 -0.42 0.13 16 51 H 0.06 0.54 8.14 -51.93 -0.42 0.12 16 52 H 0.08 1.11 8.14 -51.92 -0.42 0.69 16 53 H 0.11 1.61 6.60 -51.93 -0.34 1.27 16 LS Contribution 405.47 15.07 6.11 6.11 Total: -1.00 -32.07 405.47 -7.23 -39.31 By element: Atomic # 1 Polarization: 16.23 SS G_CDS: -9.23 Total: 6.99 kcal Atomic # 6 Polarization: 2.84 SS G_CDS: 0.88 Total: 3.72 kcal Atomic # 7 Polarization: -47.49 SS G_CDS: 0.07 Total: -47.42 kcal Atomic # 8 Polarization: -25.30 SS G_CDS: -0.47 Total: -25.77 kcal Atomic # 14 Polarization: 21.65 SS G_CDS: -4.59 Total: 17.06 kcal Total LS contribution 6.11 Total: 6.11 kcal Total: -32.07 -7.23 -39.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. ZINC000179275321.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 -91.832 kcal (2) G-P(sol) polarization free energy of solvation -32.072 kcal (3) G-ENP(sol) elect.-nuc.-pol. free energy of system -123.904 kcal (4) G-CDS(sol) cavity-dispersion-solvent structure free energy -7.234 kcal (5) G-P-CDS(sol) = G-P(sol) + G-CDS(sol) = (2) + (4) -39.306 kcal (6) G-S(sol) free energy of system = (1) + (5) -131.138 kcal 1SCF WAS SPECIFIED, SO GEOMETRY OPTIMIZATION WAS NOT USED Total computer time = 0.78 seconds