Wall clock time and date at job start Sat Apr 11 2020 04:53:48 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.53004 * 109.47105 * 2 1 4 4 N 1.46498 * 109.47177 * 180.02562 * 3 2 1 5 5 C 1.34265 * 125.75636 * 125.27413 * 4 3 2 6 6 C 1.38042 * 107.41640 * 180.18208 * 5 4 3 7 7 C 1.47430 * 126.41127 * 179.86848 * 6 5 4 8 8 O 1.21611 * 120.00205 * 7.88988 * 7 6 5 9 9 N 1.34776 * 119.99868 * 187.88623 * 7 6 5 10 10 C 1.46493 * 120.00106 * 6.44439 * 9 7 6 11 11 C 1.46506 * 119.99601 * 186.44658 * 9 7 6 12 12 C 1.53000 * 109.46787 * 269.99630 * 11 9 7 13 13 C 1.53002 * 109.47069 * 180.02562 * 12 11 9 14 14 H 1.09000 * 109.47290 * 55.00093 * 13 12 11 15 15 N 1.46509 * 109.47003 * 295.00180 * 13 12 11 16 16 C 1.36937 * 119.99964 * 154.99747 * 15 13 12 17 17 H 1.08004 * 119.99485 * 0.02562 * 16 15 13 18 18 C 1.38806 * 120.00226 * 179.97438 * 16 15 13 19 19 N 1.31628 * 119.99505 * 359.73076 * 18 16 15 20 20 Si 1.86798 * 120.00363 * 180.02562 * 18 16 15 21 21 H 1.48501 * 109.46875 * 269.72846 * 20 18 16 22 22 H 1.48498 * 109.46759 * 29.72443 * 20 18 16 23 23 H 1.48499 * 109.47068 * 149.72642 * 20 18 16 24 24 C 1.53001 * 109.47252 * 175.00399 * 13 12 11 25 25 C 1.53002 * 109.46865 * 299.99959 * 24 13 12 26 26 C 1.52996 * 109.47355 * 59.99937 * 24 13 12 27 27 C 1.41287 * 107.17443 * 359.60875 * 6 5 4 28 28 N 1.30450 * 108.02060 * 0.23767 * 27 6 5 29 29 H 1.09003 * 109.47005 * 299.99353 * 1 2 3 30 30 H 1.08996 * 109.47370 * 59.99747 * 1 2 3 31 31 H 1.09004 * 109.46962 * 180.02562 * 1 2 3 32 32 H 1.09004 * 109.47351 * 240.00237 * 2 1 3 33 33 H 1.08996 * 109.47370 * 120.00253 * 2 1 3 34 34 H 1.09003 * 109.46661 * 60.00437 * 3 2 1 35 35 H 1.08996 * 109.47213 * 300.00348 * 3 2 1 36 36 H 1.07996 * 126.28719 * 0.03296 * 5 4 3 37 37 H 1.08996 * 109.47086 * 275.86493 * 10 9 7 38 38 H 1.09003 * 109.47266 * 35.86618 * 10 9 7 39 39 H 1.09007 * 109.47268 * 155.86174 * 10 9 7 40 40 H 1.08996 * 109.46719 * 150.00381 * 11 9 7 41 41 H 1.08992 * 109.47439 * 29.99931 * 11 9 7 42 42 H 1.09001 * 109.46942 * 60.00364 * 12 11 9 43 43 H 1.08998 * 109.47417 * 300.00780 * 12 11 9 44 44 H 0.96996 * 119.99982 * 335.00166 * 15 13 12 45 45 H 0.96993 * 119.99975 * 0.27094 * 19 18 16 46 46 H 1.08996 * 109.46950 * 180.02562 * 24 13 12 47 47 H 1.08996 * 109.47205 * 60.00317 * 25 24 13 48 48 H 1.08997 * 109.46970 * 180.02562 * 25 24 13 49 49 H 1.09000 * 109.46950 * 299.99861 * 25 24 13 50 50 H 1.08999 * 109.46986 * 59.99771 * 26 24 13 51 51 H 1.08997 * 109.47682 * 180.02562 * 26 24 13 52 52 H 1.09002 * 109.46941 * 300.00336 * 26 24 13 53 53 H 1.07991 * 125.99075 * 180.23798 * 27 6 5 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.0400 1.4425 0.0000 4 7 3.5049 1.4426 -0.0006 5 6 4.2891 2.0719 -0.8904 6 6 5.5996 1.8176 -0.5388 7 6 6.8033 2.3056 -1.2363 8 8 6.6979 3.1184 -2.1347 9 7 8.0205 1.8510 -0.8782 10 6 8.1390 0.7707 0.1040 11 6 9.2250 2.4374 -1.4713 12 6 9.5939 1.6686 -2.7416 13 6 10.8522 2.2804 -3.3608 14 1 11.6496 2.3007 -2.6180 15 7 10.5633 3.6465 -3.8045 16 6 11.5769 4.5614 -3.9085 17 1 12.5894 4.2757 -3.6638 18 6 11.3034 5.8555 -4.3294 19 7 10.0705 6.2017 -4.6337 20 14 12.6862 7.1033 -4.4720 21 1 12.8358 7.8257 -3.1832 22 1 13.9543 6.4018 -4.7963 23 1 12.3677 8.0740 -5.5496 24 6 11.2927 1.4374 -4.5591 25 6 11.5943 0.0108 -4.0957 26 6 10.1734 1.4090 -5.6018 27 6 5.5714 1.0072 0.6182 28 7 4.3230 0.7907 0.9289 29 1 -0.3633 0.5137 0.8901 30 1 -0.3634 0.5138 -0.8899 31 1 -0.3633 -1.0277 0.0005 32 1 1.8934 -0.5138 -0.8900 33 1 1.8933 -0.5138 0.8899 34 1 1.6767 1.9563 0.8901 35 1 1.6766 1.9564 -0.8899 36 1 3.9558 2.6657 -1.7286 37 1 8.0995 1.1884 1.1100 38 1 7.3178 0.0664 -0.0295 39 1 9.0880 0.2533 -0.0376 40 1 10.0472 2.3770 -0.7584 41 1 9.0360 3.4815 -1.7205 42 1 8.7716 1.7289 -3.4546 43 1 9.7828 0.6244 -2.4925 44 1 9.6541 3.9029 -4.0246 45 1 9.3526 5.5537 -4.5600 46 1 12.1893 1.8731 -4.9999 47 1 10.6979 -0.4252 -3.6547 48 1 11.9085 -0.5897 -4.9493 49 1 12.3918 0.0310 -3.3529 50 1 9.9585 2.4253 -5.9319 51 1 10.4870 0.8080 -6.4554 52 1 9.2770 0.9731 -5.1607 53 1 6.4323 0.6325 1.1517 RHF calculation, no. of doubly occupied orbitals= 65 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 ZINC001485437576.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 04:53:48 Heat of formation + Delta-G solvation = 52.314858 kcal Electronic energy + Delta-G solvation = -29751.613697 eV Core-core repulsion = 25853.668783 eV Total energy + Delta-G solvation = -3897.944913 eV No. of doubly occupied orbitals = 65 Molecular weight (most abundant/longest-lived isotopes) = 336.227 amu Computer time = 0.80 seconds Orbital eigenvalues (eV) -41.60306 -39.28021 -37.40300 -35.81740 -35.33729 -33.35592 -32.23232 -31.23007 -30.26814 -29.22456 -27.14297 -26.75640 -26.55759 -24.89676 -23.09873 -22.66928 -22.44016 -21.41568 -20.87698 -19.32880 -18.40046 -17.47577 -17.16425 -16.36823 -15.97830 -15.84555 -15.47395 -15.01996 -14.78644 -14.43352 -14.21037 -13.96032 -13.64631 -13.62669 -13.50769 -13.26845 -13.15629 -12.89710 -12.79718 -12.74189 -12.43827 -11.94149 -11.85901 -11.82327 -11.61284 -11.53575 -11.44475 -11.17048 -11.05048 -10.79975 -10.69835 -10.34258 -10.24257 -10.07467 -9.84464 -9.78990 -9.76210 -9.50253 -9.01036 -8.82633 -8.76623 -8.35289 -6.84947 -5.93988 -2.93144 1.43620 2.08470 2.10060 3.22123 3.57783 3.58621 3.64001 3.77787 4.03339 4.33456 4.52471 4.58763 4.71633 4.78590 4.81436 4.90127 4.97517 5.06157 5.13442 5.26913 5.30166 5.38637 5.46187 5.47941 5.48756 5.52172 5.55527 5.60120 5.69772 5.72907 5.80893 5.83304 5.95651 6.00536 6.12544 6.20276 6.22837 6.34293 6.35834 6.44546 6.49154 6.61590 6.69517 6.80532 6.86835 6.94984 7.11952 7.14479 7.35350 7.63921 7.95517 8.49702 8.92168 9.40423 10.26696 10.42880 11.54556 Molecular weight = 336.23amu Principal moments of inertia in cm(-1) A = 0.014043 B = 0.002485 C = 0.002309 Principal moments of inertia in units of 10**(-40)*gram-cm**2 A = 1993.462545 B =11266.657239 C =12120.959988 Net atomic charges, atomic populations, and dipole contributions calculated with CM2 Atom NO. Type Charge No. of electrons 1 C -0.151 4.151 2 C -0.136 4.136 3 C 0.129 3.871 4 N -0.363 5.363 5 C 0.125 3.875 6 C -0.286 4.286 7 C 0.587 3.413 8 O -0.533 6.533 9 N -0.613 5.613 10 C 0.068 3.932 11 C 0.123 3.877 12 C -0.138 4.138 13 C 0.191 3.809 14 H 0.045 0.955 15 N -0.730 5.730 16 C -0.327 4.327 17 H 0.062 0.938 18 C -0.046 4.046 19 N -0.945 5.945 20 Si 1.080 2.920 21 H -0.297 1.297 22 H -0.297 1.297 23 H -0.290 1.290 24 C -0.114 4.114 25 C -0.144 4.144 26 C -0.135 4.135 27 C 0.061 3.939 28 N -0.300 5.300 29 H 0.059 0.941 30 H 0.058 0.942 31 H 0.062 0.938 32 H 0.068 0.932 33 H 0.078 0.922 34 H 0.088 0.912 35 H 0.099 0.901 36 H 0.187 0.813 37 H 0.057 0.943 38 H 0.076 0.924 39 H 0.073 0.927 40 H 0.070 0.930 41 H 0.094 0.906 42 H 0.066 0.934 43 H 0.059 0.941 44 H 0.370 0.630 45 H 0.258 0.742 46 H 0.067 0.933 47 H 0.050 0.950 48 H 0.042 0.958 49 H 0.047 0.953 50 H 0.063 0.937 51 H 0.035 0.965 52 H 0.043 0.957 53 H 0.204 0.796 For the charges calculated by CM2: Dipole moment (debyes) X Y Z Total from point charges -14.641 -16.623 11.794 25.096 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.208 4.208 2 C -0.174 4.174 3 C 0.008 3.992 4 N -0.142 5.142 5 C -0.026 4.026 6 C -0.301 4.301 7 C 0.375 3.625 8 O -0.410 6.410 9 N -0.346 5.346 10 C -0.075 4.075 11 C 0.002 3.998 12 C -0.177 4.177 13 C 0.080 3.920 14 H 0.063 0.937 15 N -0.373 5.373 16 C -0.459 4.459 17 H 0.080 0.920 18 C -0.244 4.244 19 N -0.590 5.590 20 Si 0.911 3.089 21 H -0.225 1.225 22 H -0.224 1.224 23 H -0.216 1.216 24 C -0.135 4.135 25 C -0.202 4.202 26 C -0.192 4.192 27 C -0.126 4.126 28 N -0.120 5.120 29 H 0.078 0.922 30 H 0.077 0.923 31 H 0.081 0.919 32 H 0.087 0.913 33 H 0.097 0.903 34 H 0.106 0.894 35 H 0.117 0.883 36 H 0.204 0.796 37 H 0.075 0.925 38 H 0.094 0.906 39 H 0.092 0.908 40 H 0.089 0.911 41 H 0.113 0.887 42 H 0.085 0.915 43 H 0.078 0.922 44 H 0.204 0.796 45 H 0.070 0.930 46 H 0.085 0.915 47 H 0.069 0.931 48 H 0.061 0.939 49 H 0.066 0.934 50 H 0.082 0.918 51 H 0.054 0.946 52 H 0.062 0.938 53 H 0.220 0.780 Dipole moment (debyes) X Y Z Total from point charges -15.061 -15.203 11.796 24.436 hybrid contribution 0.059 -1.222 -0.216 1.242 sum -15.002 -16.425 11.580 25.079 Atomic orbital electron populations 1.21714 0.94281 1.02616 1.02238 1.21799 0.96507 0.95186 1.03911 1.21934 0.73577 0.99352 1.04313 1.48591 1.06840 1.37957 1.20805 1.23012 0.90171 0.94970 0.94444 1.20335 0.93580 1.12919 1.03265 1.17294 0.82063 0.80895 0.82241 1.90853 1.86052 1.34007 1.30122 1.47669 1.06289 1.37238 1.43372 1.22240 1.03167 0.90203 0.91888 1.21556 0.84470 0.99482 0.94334 1.22267 0.99673 0.98608 0.97147 1.20094 0.93671 0.85366 0.92828 0.93701 1.42478 1.07697 1.06237 1.80874 1.20019 0.89606 0.93395 1.42844 0.92017 1.25655 1.01918 0.96901 0.99926 1.70794 0.99606 1.37270 1.51331 0.83239 0.73672 0.75636 0.76305 1.22543 1.22364 1.21596 1.21754 0.98951 0.94400 0.98428 1.21898 1.00621 0.97420 1.00274 1.21698 0.97914 1.01960 0.97666 1.24986 0.95529 0.96509 0.95527 1.77125 1.01862 1.19645 1.13399 0.92239 0.92289 0.91863 0.91283 0.90282 0.89350 0.88261 0.79582 0.92497 0.90610 0.90784 0.91117 0.88724 0.91534 0.92234 0.79643 0.93003 0.91498 0.93108 0.93880 0.93402 0.91816 0.94614 0.93810 0.77982 Number of geometries 1 Number of calculations of the screened coulomb radii 1 The total number of SCF iterations 14. 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.63 9.91 37.16 0.37 -0.26 16 2 C -0.14 -0.81 6.11 -26.73 -0.16 -0.97 16 3 C 0.13 0.83 6.49 -4.04 -0.03 0.81 16 4 N -0.36 -3.70 3.65 -60.81 -0.22 -3.92 16 5 C 0.13 1.46 11.26 -16.56 -0.19 1.27 16 6 C -0.29 -3.82 5.92 -105.07 -0.62 -4.44 16 7 C 0.59 9.74 7.65 -12.52 -0.10 9.64 16 8 O -0.53 -11.27 16.56 5.27 0.09 -11.19 16 9 N -0.61 -8.89 2.87 -182.85 -0.53 -9.42 16 10 C 0.07 0.65 8.64 59.84 0.52 1.17 16 11 C 0.12 2.10 5.22 -4.04 -0.02 2.08 16 12 C -0.14 -2.37 4.36 -26.73 -0.12 -2.49 16 13 C 0.19 3.82 2.11 -67.92 -0.14 3.68 16 14 H 0.05 0.92 8.08 -51.93 -0.42 0.50 16 15 N -0.73 -18.35 3.81 -53.42 -0.20 -18.55 16 16 C -0.33 -9.24 9.93 -15.06 -0.15 -9.39 16 17 H 0.06 1.72 7.83 -52.48 -0.41 1.31 16 18 C -0.05 -1.37 5.50 -17.43 -0.10 -1.46 16 19 N -0.94 -30.20 13.70 55.49 0.76 -29.44 16 20 Si 1.08 26.76 29.92 -169.99 -5.09 21.67 16 21 H -0.30 -7.18 7.11 56.52 0.40 -6.78 16 22 H -0.30 -6.94 7.11 56.52 0.40 -6.54 16 23 H -0.29 -7.29 7.11 56.52 0.40 -6.89 16 24 C -0.11 -2.10 2.91 -90.62 -0.26 -2.36 16 25 C -0.14 -2.11 8.85 37.16 0.33 -1.78 16 26 C -0.14 -2.50 8.68 37.16 0.32 -2.17 16 27 C 0.06 0.68 9.86 -16.98 -0.17 0.51 16 28 N -0.30 -3.46 12.33 30.86 0.38 -3.08 16 29 H 0.06 0.22 8.14 -51.93 -0.42 -0.20 16 30 H 0.06 0.21 8.14 -51.93 -0.42 -0.21 16 31 H 0.06 0.25 8.14 -51.93 -0.42 -0.17 16 32 H 0.07 0.47 8.14 -51.93 -0.42 0.04 16 33 H 0.08 0.57 8.14 -51.93 -0.42 0.15 16 34 H 0.09 0.50 8.14 -51.93 -0.42 0.08 16 35 H 0.10 0.44 8.06 -51.93 -0.42 0.02 16 36 H 0.19 1.87 8.06 -52.49 -0.42 1.44 16 37 H 0.06 0.45 7.26 -51.93 -0.38 0.07 16 38 H 0.08 0.63 5.44 -51.93 -0.28 0.35 16 39 H 0.07 0.61 8.03 -51.92 -0.42 0.20 16 40 H 0.07 1.09 8.07 -51.93 -0.42 0.67 16 41 H 0.09 2.04 7.13 -51.93 -0.37 1.67 16 42 H 0.07 1.29 6.71 -51.93 -0.35 0.94 16 43 H 0.06 0.88 6.17 -51.93 -0.32 0.56 16 44 H 0.37 10.05 6.63 -40.82 -0.27 9.78 16 45 H 0.26 8.30 6.77 -40.82 -0.28 8.02 16 46 H 0.07 1.38 8.14 -51.93 -0.42 0.96 16 47 H 0.05 0.69 6.32 -51.93 -0.33 0.36 16 48 H 0.04 0.58 8.14 -51.93 -0.42 0.16 16 49 H 0.05 0.66 8.14 -51.93 -0.42 0.23 16 50 H 0.06 1.38 7.60 -51.93 -0.39 0.98 16 51 H 0.03 0.59 8.14 -51.93 -0.42 0.17 16 52 H 0.04 0.76 6.54 -51.93 -0.34 0.42 16 53 H 0.20 1.65 4.53 -52.49 -0.24 1.42 16 LS Contribution 418.21 15.07 6.30 6.30 Total: -1.00 -35.99 418.21 -8.10 -44.09 By element: Atomic # 1 Polarization: 18.78 SS G_CDS: -9.08 Total: 9.70 kcal Atomic # 6 Polarization: -5.66 SS G_CDS: -0.51 Total: -6.17 kcal Atomic # 7 Polarization: -64.60 SS G_CDS: 0.19 Total: -64.41 kcal Atomic # 8 Polarization: -11.27 SS G_CDS: 0.09 Total: -11.19 kcal Atomic # 14 Polarization: 26.76 SS G_CDS: -5.09 Total: 21.67 kcal Total LS contribution 6.30 Total: 6.30 kcal Total: -35.99 -8.10 -44.09 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. ZINC001485437576.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 96.407 kcal (2) G-P(sol) polarization free energy of solvation -35.994 kcal (3) G-ENP(sol) elect.-nuc.-pol. free energy of system 60.413 kcal (4) G-CDS(sol) cavity-dispersion-solvent structure free energy -8.098 kcal (5) G-P-CDS(sol) = G-P(sol) + G-CDS(sol) = (2) + (4) -44.092 kcal (6) G-S(sol) free energy of system = (1) + (5) 52.315 kcal 1SCF WAS SPECIFIED, SO GEOMETRY OPTIMIZATION WAS NOT USED Total computer time = 0.81 seconds