Wall clock time and date at job start Tue Mar 31 2020 21:12:31 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 N 2 2 C 1.31617 * 1 3 3 Si 1.86806 * 119.99945 * 2 1 4 4 H 1.48499 * 109.46965 * 240.00187 * 3 2 1 5 5 H 1.48500 * 109.46860 * 359.97438 * 3 2 1 6 6 H 1.48497 * 109.46952 * 119.99977 * 3 2 1 7 7 C 1.38815 * 120.00210 * 180.02562 * 2 1 3 8 8 H 1.07999 * 120.00039 * 179.97438 * 7 2 1 9 9 N 1.36931 * 120.00053 * 359.97438 * 7 2 1 10 10 C 1.47432 * 125.64856 * 36.29236 * 9 7 2 11 11 C 1.54904 * 104.82914 * 204.09544 * 10 9 7 12 12 C 1.55153 * 101.58402 * 322.98943 * 11 10 9 13 13 C 1.47021 * 125.65163 * 216.25542 * 9 7 2 14 14 H 1.09002 * 109.87948 * 59.39936 * 13 9 7 15 15 C 1.53005 * 109.88030 * 298.39439 * 13 9 7 16 16 H 1.08992 * 109.69943 * 194.50300 * 15 13 9 17 17 C 1.53187 * 109.58472 * 314.85658 * 15 13 9 18 18 C 1.53038 * 109.31502 * 174.42155 * 17 15 13 19 19 C 1.53038 * 109.53889 * 298.64208 * 18 17 15 20 20 C 1.53180 * 109.31735 * 61.35673 * 19 18 17 21 21 N 1.46923 * 109.58469 * 74.14611 * 15 13 9 22 22 C 1.34775 * 120.63795 * 6.88273 * 21 15 13 23 23 O 1.21286 * 119.99503 * 354.04137 * 22 21 15 24 24 C 1.50695 * 120.00650 * 174.04863 * 22 21 15 25 25 C 1.53001 * 115.54790 * 304.44393 * 24 22 21 26 26 C 1.53001 * 109.46935 * 293.95658 * 25 24 22 27 27 F 1.39907 * 109.46827 * 180.61006 * 26 25 24 28 28 C 1.52999 * 117.50131 * 158.77142 * 24 22 21 29 29 C 1.52996 * 117.49988 * 90.10749 * 24 22 21 30 30 H 0.96995 * 120.00074 * 0.02562 * 1 2 3 31 31 H 1.09001 * 110.36703 * 322.93557 * 10 9 7 32 32 H 1.08997 * 110.36330 * 85.25623 * 10 9 7 33 33 H 1.09003 * 111.00671 * 81.07337 * 11 10 9 34 34 H 1.09009 * 110.99922 * 204.91654 * 11 10 9 35 35 H 1.08997 * 110.72343 * 277.13955 * 12 11 10 36 36 H 1.09005 * 110.71925 * 153.87235 * 12 11 10 37 37 H 1.09008 * 109.49815 * 54.47167 * 17 15 13 38 38 H 1.09007 * 109.52909 * 294.39373 * 17 15 13 39 39 H 1.09004 * 109.46386 * 178.61354 * 18 17 15 40 40 H 1.09004 * 109.44545 * 58.64688 * 18 17 15 41 41 H 1.09008 * 109.49651 * 301.40510 * 19 18 17 42 42 H 1.09002 * 109.52690 * 181.33131 * 19 18 17 43 43 H 1.09006 * 109.59096 * 185.57850 * 20 19 18 44 44 H 1.08999 * 109.70249 * 65.22686 * 20 19 18 45 45 H 1.09000 * 109.47225 * 53.95539 * 25 24 22 46 46 H 1.08997 * 109.47558 * 173.95697 * 25 24 22 47 47 H 1.08993 * 109.47190 * 300.60877 * 26 25 24 48 48 H 1.09001 * 109.47405 * 60.61320 * 26 25 24 49 49 H 1.08997 * 117.49566 * 0.02562 * 28 24 22 50 50 H 1.09002 * 117.49446 * 145.01919 * 28 24 22 51 51 H 1.08999 * 117.50323 * 214.97843 * 29 24 22 52 52 H 1.09000 * 117.49930 * 359.97438 * 29 24 22 Note: An asterisk (*) indicates this parameter was or will be optimized. Cartesian coordinates (angstroms) NO. Atom X Y Z 1 7 0.0000 0.0000 0.0000 2 6 1.3162 0.0000 0.0000 3 14 2.2502 1.6178 0.0000 4 1 3.1039 1.6965 -1.2125 5 1 1.2851 2.7465 -0.0006 6 1 3.1039 1.6964 1.2125 7 6 2.0103 -1.2021 -0.0005 8 1 3.0903 -1.2021 -0.0010 9 7 1.3257 -2.3880 -0.0005 10 6 0.0595 -2.6493 -0.7092 11 6 -0.5932 -3.8219 0.0644 12 6 0.6465 -4.6761 0.4393 13 6 1.7318 -3.6120 0.7055 14 1 1.8071 -3.4153 1.7750 15 6 3.0801 -4.0984 0.1701 16 1 3.2578 -5.1225 0.4981 17 6 3.0748 -4.0432 -1.3608 18 6 4.4684 -4.3970 -1.8849 19 6 5.4841 -3.3757 -1.3681 20 6 5.5167 -3.4234 0.1626 21 7 4.1521 -3.2277 0.6711 22 6 3.8901 -2.2713 1.5838 23 8 2.7456 -2.0506 1.9192 24 6 5.0189 -1.4761 2.1875 25 6 6.1222 -2.2990 2.8555 26 6 5.5542 -3.0108 4.0850 27 9 6.5665 -3.7522 4.7038 28 6 4.6753 -0.1201 2.8073 29 6 5.4318 -0.1850 1.4781 30 1 -0.4850 0.8400 -0.0004 31 1 -0.5823 -1.7689 -0.6761 32 1 0.2555 -2.9383 -1.7418 33 1 -1.1056 -3.4637 0.9573 34 1 -1.2742 -4.3818 -0.5767 35 1 0.9349 -5.3221 -0.3899 36 1 0.4532 -5.2629 1.3374 37 1 2.8075 -3.0382 -1.6876 38 1 2.3480 -4.7571 -1.7487 39 1 4.4610 -4.3825 -2.9748 40 1 4.7450 -5.3923 -1.5369 41 1 5.1943 -2.3770 -1.6949 42 1 6.4726 -3.6143 -1.7607 43 1 6.1641 -2.6314 0.5393 44 1 5.8941 -4.3914 0.4922 45 1 6.5013 -3.0387 2.1503 46 1 6.9339 -1.6387 3.1607 47 1 5.1653 -2.2717 4.7853 48 1 4.7498 -3.6794 3.7785 49 1 3.6327 0.1966 2.7788 50 1 5.2357 0.1858 3.6908 51 1 6.4894 0.0785 1.4877 52 1 4.8873 0.0891 0.5745 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) F: (AM1): M.J.S. DEWAR AND E. G. ZOEBISCH, THEOCHEM, 180, 1 (1988). 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 ZINC001491280984.mol2 52 1SCF WAS SPECIFIED, SO GEOMETRY OPTIMIZATION WAS NOT USED SCF FIELD WAS ACHIEVED Wall clock time and date at job start Tue Mar 31 2020 21:12:31 Heat of formation + Delta-G solvation = 10.611672 kcal Electronic energy + Delta-G solvation = -33946.360819 eV Core-core repulsion = 29863.542985 eV Total energy + Delta-G solvation = -4082.817834 eV No. of doubly occupied orbitals = 65 Molecular weight (most abundant/longest-lived isotopes) = 338.211 amu Computer time = 0.98 seconds Orbital eigenvalues (eV) -47.55394 -40.52524 -38.00093 -36.50036 -35.90355 -32.44318 -32.30472 -31.17433 -30.72447 -29.03810 -27.20291 -25.02144 -24.30159 -24.06981 -22.60508 -22.18007 -21.78872 -20.11816 -19.73946 -19.33612 -17.91732 -16.85735 -16.32497 -16.17794 -15.63904 -15.44282 -15.16239 -14.71056 -14.45574 -14.29886 -14.09339 -13.86131 -13.79493 -13.49525 -13.11579 -12.99649 -12.53932 -12.36412 -12.09653 -11.86646 -11.73532 -11.61075 -11.44223 -11.27329 -11.08175 -10.93239 -10.82549 -10.72778 -10.68833 -10.60743 -10.43231 -10.20183 -10.03385 -9.97449 -9.92946 -9.76217 -9.69818 -9.59834 -9.36054 -8.30166 -8.02127 -7.62008 -6.17160 -5.45971 -3.01858 3.30172 3.42729 3.46438 3.48808 3.57035 4.04738 4.57637 4.60009 4.82155 4.89411 4.93848 5.08092 5.15460 5.19812 5.22845 5.31723 5.41409 5.51428 5.54408 5.62947 5.69079 5.73537 5.74472 5.82401 5.86353 5.90697 5.92601 6.00442 6.05563 6.08238 6.13102 6.15716 6.23330 6.29615 6.34013 6.43600 6.46056 6.55981 6.58536 6.59533 6.75446 6.87205 6.91817 7.14091 7.29188 7.49138 7.64313 7.87330 8.06670 8.24193 8.70087 9.20048 9.38434 10.00276 10.81890 11.50256 Molecular weight = 338.21amu Principal moments of inertia in cm(-1) A = 0.009748 B = 0.006719 C = 0.005343 Principal moments of inertia in units of 10**(-40)*gram-cm**2 A = 2871.709610 B = 4166.068135 C = 5239.155997 Net atomic charges, atomic populations, and dipole contributions calculated with CM2 Atom NO. Type Charge No. of electrons 1 N -0.904 5.904 2 C -0.001 4.001 3 Si 0.879 3.121 4 H -0.290 1.290 5 H -0.295 1.295 6 H -0.272 1.272 7 C -0.168 4.168 8 H 0.061 0.939 9 N -0.565 5.565 10 C 0.121 3.879 11 C -0.130 4.130 12 C -0.134 4.134 13 C 0.140 3.860 14 H 0.088 0.912 15 C 0.144 3.856 16 H 0.049 0.951 17 C -0.124 4.124 18 C -0.117 4.117 19 C -0.133 4.133 20 C 0.110 3.890 21 N -0.603 5.603 22 C 0.514 3.486 23 O -0.438 6.438 24 C -0.154 4.154 25 C -0.114 4.114 26 C 0.064 3.936 27 F -0.224 7.224 28 C -0.144 4.144 29 C -0.159 4.159 30 H 0.249 0.751 31 H 0.106 0.894 32 H -0.006 1.006 33 H 0.059 0.941 34 H 0.047 0.953 35 H 0.053 0.947 36 H 0.051 0.949 37 H 0.092 0.908 38 H 0.057 0.943 39 H 0.054 0.946 40 H 0.051 0.949 41 H 0.068 0.932 42 H 0.062 0.938 43 H 0.083 0.917 44 H 0.054 0.946 45 H 0.086 0.914 46 H 0.076 0.924 47 H 0.076 0.924 48 H 0.076 0.924 49 H 0.114 0.886 50 H 0.093 0.907 51 H 0.095 0.905 52 H 0.100 0.900 For the charges calculated by CM2: Dipole moment (debyes) X Y Z Total from point charges 11.580 -7.519 1.715 13.913 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 N -0.550 5.550 2 C -0.198 4.198 3 Si 0.710 3.290 4 H -0.217 1.217 5 H -0.221 1.221 6 H -0.199 1.199 7 C -0.297 4.297 8 H 0.079 0.921 9 N -0.288 5.288 10 C -0.004 4.004 11 C -0.168 4.168 12 C -0.173 4.173 13 C 0.033 3.967 14 H 0.105 0.895 15 C 0.039 3.961 16 H 0.067 0.933 17 C -0.163 4.163 18 C -0.154 4.154 19 C -0.171 4.171 20 C -0.014 4.014 21 N -0.339 5.339 22 C 0.301 3.699 23 O -0.309 6.309 24 C -0.157 4.157 25 C -0.151 4.151 26 C 0.009 3.991 27 F -0.205 7.205 28 C -0.182 4.182 29 C -0.196 4.196 30 H 0.060 0.940 31 H 0.123 0.877 32 H 0.012 0.988 33 H 0.078 0.922 34 H 0.066 0.934 35 H 0.072 0.928 36 H 0.070 0.930 37 H 0.111 0.889 38 H 0.076 0.924 39 H 0.073 0.927 40 H 0.070 0.930 41 H 0.087 0.913 42 H 0.080 0.920 43 H 0.101 0.899 44 H 0.072 0.928 45 H 0.104 0.896 46 H 0.095 0.905 47 H 0.095 0.905 48 H 0.094 0.906 49 H 0.132 0.868 50 H 0.112 0.888 51 H 0.113 0.887 52 H 0.118 0.882 Dipole moment (debyes) X Y Z Total from point charges 10.878 -7.157 1.948 13.166 hybrid contribution 0.672 0.850 -0.261 1.115 sum 11.550 -6.307 1.687 13.267 Atomic orbital electron populations 1.69571 1.05213 1.24500 1.55760 1.25271 0.95241 0.97445 1.01878 0.86755 0.79276 0.84796 0.78179 1.21696 1.22106 1.19851 1.18649 0.89637 0.81758 1.39679 0.92120 1.44105 1.24371 1.07054 1.53238 1.21991 0.90168 0.97314 0.90911 1.22723 0.96103 0.97783 1.00192 1.22701 0.97615 0.97558 0.99421 1.21889 0.88620 0.87006 0.99137 0.89501 1.21895 0.89375 0.91161 0.93655 0.93305 1.21905 0.94735 1.04081 0.95557 1.21515 0.96279 0.97929 0.99715 1.21885 0.98826 1.00733 0.95674 1.21699 0.84768 1.00847 0.94052 1.47645 1.07908 1.35437 1.42902 1.18921 0.88395 0.80938 0.81683 1.90878 1.19208 1.62055 1.58731 1.21400 0.98578 0.90123 1.05643 1.21847 0.96768 0.99602 0.96917 1.21950 0.88017 0.94770 0.94396 1.92738 1.63649 1.79410 1.84700 1.23051 1.02989 0.96250 0.95871 1.23059 1.00559 0.98190 0.97807 0.94031 0.87655 0.98774 0.92231 0.93421 0.92830 0.92999 0.88934 0.92416 0.92691 0.93003 0.91301 0.91977 0.89899 0.92768 0.89609 0.90498 0.90533 0.90614 0.86819 0.88828 0.88660 0.88151 Number of geometries 1 Number of calculations of the screened coulomb radii 1 The total number of SCF iterations 15. 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 N -0.90 -26.52 10.78 55.49 0.60 -25.92 16 2 C 0.00 -0.04 5.44 -17.43 -0.09 -0.13 16 3 Si 0.88 21.56 28.63 -169.99 -4.87 16.69 16 4 H -0.29 -7.05 7.11 56.52 0.40 -6.64 16 5 H -0.29 -7.50 7.11 56.52 0.40 -7.10 16 6 H -0.27 -6.32 6.44 56.52 0.36 -5.95 16 7 C -0.17 -4.37 5.91 -15.06 -0.09 -4.46 16 8 H 0.06 1.55 4.60 -52.49 -0.24 1.30 16 9 N -0.56 -13.67 1.80 -162.57 -0.29 -13.96 16 10 C 0.12 2.91 6.83 -2.53 -0.02 2.89 16 11 C -0.13 -2.62 7.08 -24.59 -0.17 -2.80 16 12 C -0.13 -2.47 6.33 -24.90 -0.16 -2.63 16 13 C 0.14 2.92 2.31 -66.96 -0.15 2.77 16 14 H 0.09 2.06 5.35 -55.03 -0.29 1.76 16 15 C 0.14 2.54 2.53 -67.60 -0.17 2.36 16 16 H 0.05 0.80 8.11 -51.93 -0.42 0.38 16 17 C -0.12 -2.12 4.94 -26.62 -0.13 -2.25 16 18 C -0.12 -1.59 5.88 -26.69 -0.16 -1.75 16 19 C -0.13 -1.74 6.08 -26.62 -0.16 -1.90 16 20 C 0.11 1.46 4.91 -3.71 -0.02 1.44 16 21 N -0.60 -10.08 2.58 -165.23 -0.43 -10.50 16 22 C 0.51 9.52 4.00 -10.99 -0.04 9.47 16 23 O -0.44 -10.05 9.14 5.55 0.05 -10.00 16 24 C -0.15 -2.33 0.85 -155.66 -0.13 -2.46 16 25 C -0.11 -1.37 5.05 -26.73 -0.14 -1.50 16 26 C 0.06 0.80 7.42 37.16 0.28 1.07 16 27 F -0.22 -3.02 17.36 2.25 0.04 -2.98 16 28 C -0.14 -2.29 9.73 -26.69 -0.26 -2.55 16 29 C -0.16 -2.40 9.67 -26.69 -0.26 -2.66 16 30 H 0.25 7.21 8.31 -40.82 -0.34 6.87 16 31 H 0.11 2.91 5.38 -51.93 -0.28 2.63 16 32 H -0.01 -0.13 8.14 -51.93 -0.42 -0.55 16 33 H 0.06 1.28 8.14 -51.93 -0.42 0.86 16 34 H 0.05 0.85 8.14 -51.92 -0.42 0.43 16 35 H 0.05 0.90 6.67 -51.93 -0.35 0.56 16 36 H 0.05 0.87 8.14 -51.93 -0.42 0.45 16 37 H 0.09 1.95 6.49 -51.92 -0.34 1.61 16 38 H 0.06 0.96 7.01 -51.93 -0.36 0.60 16 39 H 0.05 0.71 8.14 -51.93 -0.42 0.29 16 40 H 0.05 0.64 8.14 -51.93 -0.42 0.22 16 41 H 0.07 1.05 8.14 -51.92 -0.42 0.62 16 42 H 0.06 0.66 8.14 -51.93 -0.42 0.24 16 43 H 0.08 1.02 4.84 -51.93 -0.25 0.77 16 44 H 0.05 0.64 8.12 -51.93 -0.42 0.22 16 45 H 0.09 0.97 5.27 -51.93 -0.27 0.70 16 46 H 0.08 0.81 8.09 -51.93 -0.42 0.39 16 47 H 0.08 0.92 8.14 -51.93 -0.42 0.49 16 48 H 0.08 1.02 7.76 -51.93 -0.40 0.62 16 49 H 0.11 2.24 7.20 -51.93 -0.37 1.86 16 50 H 0.09 1.20 8.14 -51.93 -0.42 0.78 16 51 H 0.10 1.14 8.14 -51.93 -0.42 0.72 16 52 H 0.10 1.81 6.75 -51.93 -0.35 1.46 16 LS Contribution 375.44 15.07 5.66 5.66 Total: -1.00 -29.79 375.44 -9.72 -39.51 By element: Atomic # 1 Polarization: 15.18 SS G_CDS: -8.60 Total: 6.58 kcal Atomic # 6 Polarization: -3.19 SS G_CDS: -1.88 Total: -5.07 kcal Atomic # 7 Polarization: -50.26 SS G_CDS: -0.12 Total: -50.38 kcal Atomic # 8 Polarization: -10.05 SS G_CDS: 0.05 Total: -10.00 kcal Atomic # 9 Polarization: -3.02 SS G_CDS: 0.04 Total: -2.98 kcal Atomic # 14 Polarization: 21.56 SS G_CDS: -4.87 Total: 16.69 kcal Total LS contribution 5.66 Total: 5.66 kcal Total: -29.79 -9.72 -39.51 kcal The number of atoms in the molecule is 52 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. ZINC001491280984.mol2 52 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 50.122 kcal (2) G-P(sol) polarization free energy of solvation -29.789 kcal (3) G-ENP(sol) elect.-nuc.-pol. free energy of system 20.333 kcal (4) G-CDS(sol) cavity-dispersion-solvent structure free energy -9.721 kcal (5) G-P-CDS(sol) = G-P(sol) + G-CDS(sol) = (2) + (4) -39.510 kcal (6) G-S(sol) free energy of system = (1) + (5) 10.612 kcal 1SCF WAS SPECIFIED, SO GEOMETRY OPTIMIZATION WAS NOT USED Total computer time = 0.98 seconds