Wall clock time and date at job start Tue Mar 31 2020 19:44:32 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 O 1.42902 * 1 3 3 C 1.42902 * 114.00256 * 2 1 4 4 C 1.53001 * 109.46933 * 59.99807 * 3 2 1 5 5 C 1.53000 * 109.47066 * 299.99696 * 3 2 1 6 6 C 1.52996 * 109.47084 * 179.97438 * 3 2 1 7 7 C 1.52993 * 109.47473 * 64.99859 * 6 3 2 8 8 N 1.46906 * 109.47378 * 179.97438 * 7 6 3 9 9 C 1.47652 * 111.00095 * 266.32303 * 8 7 6 10 10 C 1.52653 * 87.64987 * 222.00226 * 9 8 7 11 11 C 1.52115 * 129.73560 * 164.47229 * 10 9 8 12 12 C 1.52118 * 129.73085 * 246.28356 * 10 9 8 13 13 H 1.09000 * 117.52308 * 133.59083 * 12 10 9 14 14 C 1.50697 * 117.52191 * 348.44269 * 12 10 9 15 15 O 1.21283 * 120.00204 * 359.97438 * 14 12 10 16 16 N 1.34775 * 120.00174 * 179.97438 * 14 12 10 17 17 C 1.47070 * 136.05405 * 359.63065 * 16 14 12 18 18 C 1.53241 * 88.08007 * 207.20104 * 17 16 14 19 19 C 1.50706 * 114.32331 * 220.47088 * 18 17 16 20 20 H 1.07998 * 119.99687 * 93.76197 * 19 18 17 21 21 C 1.37873 * 119.99577 * 273.75602 * 19 18 17 22 22 N 1.31514 * 120.00636 * 0.02562 * 21 19 18 23 23 Si 1.86803 * 119.99868 * 180.02562 * 21 19 18 24 24 H 1.48508 * 109.47044 * 0.02562 * 23 21 19 25 25 H 1.48497 * 109.47205 * 120.00003 * 23 21 19 26 26 H 1.48499 * 109.47351 * 240.00435 * 23 21 19 27 27 C 1.47072 * 136.05972 * 179.97438 * 16 14 12 28 28 C 1.47654 * 110.99850 * 170.00188 * 8 7 6 29 29 H 1.09005 * 109.47018 * 179.97438 * 1 2 3 30 30 H 1.08994 * 109.47165 * 300.00294 * 1 2 3 31 31 H 1.09001 * 109.46848 * 59.99967 * 1 2 3 32 32 H 1.08998 * 109.47077 * 179.97438 * 4 3 2 33 33 H 1.09002 * 109.47083 * 300.00043 * 4 3 2 34 34 H 1.09000 * 109.46998 * 60.00028 * 4 3 2 35 35 H 1.08997 * 109.47108 * 300.00091 * 5 3 2 36 36 H 1.09002 * 109.47044 * 59.99899 * 5 3 2 37 37 H 1.08998 * 109.46865 * 179.97438 * 5 3 2 38 38 H 1.09004 * 109.47108 * 304.99971 * 6 3 2 39 39 H 1.09004 * 109.47368 * 185.00072 * 6 3 2 40 40 H 1.09004 * 109.46862 * 59.99962 * 7 6 3 41 41 H 1.09000 * 109.47250 * 299.99985 * 7 6 3 42 42 H 1.09006 * 113.52703 * 336.69250 * 9 8 7 43 43 H 1.09001 * 113.52266 * 107.30642 * 9 8 7 44 44 H 1.08994 * 117.51972 * 226.39860 * 11 10 9 45 45 H 1.08992 * 117.51839 * 11.55177 * 11 10 9 46 46 H 1.08993 * 113.46623 * 321.92833 * 17 16 14 47 47 H 1.08995 * 113.46992 * 92.37523 * 17 16 14 48 48 H 1.09001 * 114.01829 * 87.71387 * 18 17 16 49 49 H 0.96996 * 120.00300 * 179.97438 * 22 21 19 50 50 H 1.08996 * 113.46467 * 267.62725 * 27 16 14 51 51 H 1.08998 * 113.46309 * 38.00804 * 27 16 14 52 52 H 1.09005 * 113.52349 * 252.68932 * 28 8 7 53 53 H 1.08997 * 113.53114 * 23.23676 * 28 8 7 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 8 1.4290 0.0000 0.0000 3 6 2.0103 1.3054 0.0000 4 6 1.5588 2.0647 1.2492 5 6 1.5589 2.0647 -1.2493 6 6 3.5355 1.1846 0.0006 7 6 3.9963 0.5363 -1.3063 8 7 5.4608 0.4208 -1.3059 9 6 6.0836 1.5846 -1.9676 10 6 7.0457 0.6083 -2.6396 11 6 7.9649 0.8156 -3.8337 12 6 8.5485 0.4491 -2.4654 13 1 8.9364 -0.5615 -2.3375 14 6 9.2799 1.5276 -1.7087 15 8 8.6700 2.4776 -1.2654 16 7 10.6120 1.4367 -1.5246 17 6 11.6678 0.4787 -1.8855 18 6 12.4498 0.9565 -0.6573 19 6 13.9322 1.1050 -0.8844 20 1 14.3414 2.0627 -1.1703 21 6 14.7659 0.0182 -0.7270 22 7 14.2675 -1.1481 -0.3792 23 14 16.6035 0.2024 -1.0077 24 1 16.9081 1.6054 -1.3876 25 1 17.3318 -0.1414 0.2399 26 1 17.0316 -0.7102 -2.0983 27 6 11.6483 2.2486 -0.8690 28 6 5.9264 -0.4094 -2.4347 29 1 -0.3633 -1.0277 -0.0005 30 1 -0.3633 0.5138 0.8899 31 1 -0.3633 0.5138 -0.8900 32 1 2.0025 3.0603 1.2495 33 1 0.4722 2.1508 1.2492 34 1 1.8804 1.5238 2.1392 35 1 1.8806 1.5237 -2.1392 36 1 0.4723 2.1508 -1.2493 37 1 2.0020 3.0606 -1.2490 38 1 3.8508 0.5684 0.8427 39 1 3.9786 2.1764 0.0902 40 1 3.6810 1.1525 -2.1484 41 1 3.5533 -0.4556 -1.3959 42 1 5.4208 2.0933 -2.6677 43 1 6.5722 2.2692 -1.2744 44 1 7.9691 0.0460 -4.6056 45 1 8.1102 1.8367 -4.1861 46 1 11.3510 -0.5621 -1.8185 47 1 12.1596 0.7102 -2.8303 48 1 12.2035 0.4229 0.2607 49 1 14.8540 -1.9127 -0.2689 50 1 12.1347 2.9561 -1.5406 51 1 11.3145 2.7108 0.0600 52 1 6.2836 -1.3939 -2.1325 53 1 5.2197 -0.4530 -3.2634 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 ZINC001277131305.mol2 53 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 19:44:32 Heat of formation + Delta-G solvation = 20.122858 kcal Electronic energy + Delta-G solvation = -29178.827392 eV Core-core repulsion = 25247.193777 eV Total energy + Delta-G solvation = -3931.633614 eV No. of doubly occupied orbitals = 65 Molecular weight (most abundant/longest-lived isotopes) = 336.220 amu Computer time = 0.79 seconds Orbital eigenvalues (eV) -40.71965 -39.49616 -39.01864 -37.12940 -35.21239 -33.79452 -32.67501 -30.31616 -29.22336 -27.86200 -27.29130 -26.01398 -25.72769 -24.95945 -24.93919 -23.16582 -21.63784 -20.33482 -20.03812 -18.79360 -17.83622 -16.89278 -16.71383 -16.45971 -16.23280 -15.61388 -15.37726 -15.34230 -15.05104 -14.76691 -14.45503 -14.33502 -14.06147 -13.64387 -13.33739 -13.21380 -12.92371 -12.71550 -12.58866 -12.41071 -12.33705 -12.17598 -11.94905 -11.85264 -11.71124 -11.50543 -11.40666 -11.29434 -10.96025 -10.93325 -10.77000 -10.69024 -10.43773 -10.31509 -10.22916 -10.05662 -9.84393 -9.53157 -9.19783 -8.97320 -8.69490 -8.57567 -7.82304 -6.17638 -4.23255 2.97132 3.23186 3.32106 3.35301 3.45808 3.47782 3.73758 4.16427 4.30574 4.36168 4.49074 4.53255 4.62151 4.74593 4.84940 4.89976 4.98406 5.00699 5.11277 5.17581 5.19483 5.24042 5.27727 5.33987 5.37230 5.39618 5.43434 5.49188 5.57858 5.64502 5.70315 5.70973 5.74077 5.79303 5.87733 5.91874 5.99497 6.11238 6.12398 6.19052 6.26198 6.29525 6.43362 6.46695 6.60673 6.69327 6.86091 6.95418 7.09584 7.56904 7.80437 8.19155 8.56751 8.80941 8.82822 9.43635 10.92143 Molecular weight = 336.22amu Principal moments of inertia in cm(-1) A = 0.024951 B = 0.001855 C = 0.001835 Principal moments of inertia in units of 10**(-40)*gram-cm**2 A = 1121.924639 B =15094.301166 C =15254.592691 Net atomic charges, atomic populations, and dipole contributions calculated with CM2 Atom NO. Type Charge No. of electrons 1 C 0.025 3.975 2 O -0.377 6.377 3 C 0.111 3.889 4 C -0.178 4.178 5 C -0.181 4.181 6 C -0.095 4.095 7 C 0.022 3.978 8 N -0.467 5.467 9 C 0.036 3.964 10 C -0.137 4.137 11 C -0.135 4.135 12 C -0.137 4.137 13 H 0.123 0.877 14 C 0.563 3.437 15 O -0.561 6.561 16 N -0.640 5.640 17 C 0.100 3.900 18 C 0.033 3.967 19 C -0.516 4.516 20 H 0.064 0.936 21 C 0.068 3.932 22 N -0.885 5.885 23 Si 0.895 3.105 24 H -0.271 1.271 25 H -0.283 1.283 26 H -0.283 1.283 27 C 0.101 3.899 28 C 0.041 3.959 29 H 0.083 0.917 30 H 0.042 0.958 31 H 0.042 0.958 32 H 0.068 0.932 33 H 0.058 0.942 34 H 0.066 0.934 35 H 0.066 0.934 36 H 0.057 0.943 37 H 0.067 0.933 38 H 0.079 0.921 39 H 0.077 0.923 40 H 0.044 0.956 41 H 0.078 0.922 42 H 0.066 0.934 43 H 0.113 0.887 44 H 0.105 0.895 45 H 0.100 0.900 46 H 0.059 0.941 47 H 0.074 0.926 48 H 0.068 0.932 49 H 0.263 0.737 50 H 0.072 0.928 51 H 0.057 0.943 52 H 0.086 0.914 53 H 0.072 0.928 For the charges calculated by CM2: Dipole moment (debyes) X Y Z Total from point charges -24.809 2.293 -4.587 25.333 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.072 4.072 2 O -0.296 6.296 3 C 0.071 3.929 4 C -0.235 4.235 5 C -0.238 4.238 6 C -0.132 4.132 7 C -0.107 4.107 8 N -0.188 5.188 9 C -0.092 4.092 10 C -0.138 4.138 11 C -0.173 4.173 12 C -0.158 4.158 13 H 0.141 0.859 14 C 0.353 3.647 15 O -0.441 6.441 16 N -0.373 5.373 17 C -0.022 4.022 18 C 0.014 3.986 19 C -0.554 4.554 20 H 0.082 0.918 21 C -0.136 4.136 22 N -0.522 5.522 23 Si 0.723 3.277 24 H -0.196 1.196 25 H -0.209 1.209 26 H -0.209 1.209 27 C -0.021 4.021 28 C -0.088 4.088 29 H 0.101 0.899 30 H 0.061 0.939 31 H 0.060 0.940 32 H 0.087 0.913 33 H 0.077 0.923 34 H 0.085 0.915 35 H 0.085 0.915 36 H 0.076 0.924 37 H 0.086 0.914 38 H 0.098 0.902 39 H 0.096 0.904 40 H 0.063 0.937 41 H 0.096 0.904 42 H 0.084 0.916 43 H 0.131 0.869 44 H 0.123 0.877 45 H 0.118 0.882 46 H 0.078 0.922 47 H 0.092 0.908 48 H 0.086 0.914 49 H 0.073 0.927 50 H 0.091 0.909 51 H 0.075 0.925 52 H 0.105 0.895 53 H 0.090 0.910 Dipole moment (debyes) X Y Z Total from point charges -25.057 2.271 -3.628 25.420 hybrid contribution 0.444 0.264 -1.074 1.192 sum -24.614 2.535 -4.702 25.187 Atomic orbital electron populations 1.22768 0.81921 1.02477 1.00009 1.88001 1.19150 1.27439 1.94989 1.21544 0.94052 0.83171 0.94128 1.22405 1.01526 1.01152 0.98419 1.22496 1.01532 1.01170 0.98640 1.21425 0.90861 1.02017 0.98931 1.22725 0.89522 1.01307 0.97136 1.68967 1.07721 1.07058 1.35089 1.23574 0.94172 0.94531 0.96926 1.21793 0.90185 1.00004 1.01811 1.22785 1.01326 1.03146 0.90054 1.21881 0.93591 1.00211 1.00071 0.85880 1.19020 0.80477 0.85709 0.79496 1.90555 1.66449 1.34748 1.52392 1.49929 1.04317 1.27878 1.55130 1.23540 0.84566 0.92392 1.01740 1.20321 0.93914 0.90058 0.94306 1.20944 0.90521 0.96260 1.47697 0.91781 1.25034 0.99017 0.95260 0.94273 1.69961 1.35540 1.00420 1.46308 0.86483 0.84968 0.78348 0.77909 1.19564 1.20882 1.20917 1.23340 0.85467 0.97602 0.95722 1.23333 0.94336 0.97956 0.93181 0.89896 0.93917 0.93982 0.91321 0.92299 0.91537 0.91460 0.92377 0.91350 0.90227 0.90431 0.93739 0.90399 0.91619 0.86901 0.87678 0.88168 0.92245 0.90805 0.91393 0.92685 0.90944 0.92488 0.89538 0.90951 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.02 0.17 10.21 101.05 1.03 1.20 16 2 O -0.38 -3.76 9.88 -35.23 -0.35 -4.10 16 3 C 0.11 0.98 0.82 -90.62 -0.07 0.90 16 4 C -0.18 -1.36 8.25 37.16 0.31 -1.06 16 5 C -0.18 -1.36 7.76 37.16 0.29 -1.07 16 6 C -0.09 -0.99 4.83 -26.74 -0.13 -1.12 16 7 C 0.02 0.24 4.22 -3.82 -0.02 0.22 16 8 N -0.47 -6.23 6.32 -243.97 -1.54 -7.77 16 9 C 0.04 0.50 7.31 -6.79 -0.05 0.45 16 10 C -0.14 -1.87 3.13 -155.77 -0.49 -2.36 16 11 C -0.14 -1.71 10.24 -27.08 -0.28 -1.99 16 12 C -0.14 -2.17 6.02 -92.12 -0.55 -2.72 16 13 H 0.12 1.88 8.14 -51.93 -0.42 1.46 16 14 C 0.56 11.20 7.86 -10.99 -0.09 11.11 16 15 O -0.56 -12.33 15.02 5.56 0.08 -12.25 16 16 N -0.64 -13.25 3.69 -185.35 -0.68 -13.94 16 17 C 0.10 2.12 7.98 -8.62 -0.07 2.05 16 18 C 0.03 0.78 3.77 -91.38 -0.34 0.44 16 19 C -0.52 -13.14 8.82 -34.44 -0.30 -13.44 16 20 H 0.06 1.60 7.81 -52.49 -0.41 1.19 16 21 C 0.07 1.76 5.46 -17.82 -0.10 1.66 16 22 N -0.89 -24.13 13.26 55.43 0.74 -23.40 16 23 Si 0.90 20.04 29.54 -169.99 -5.02 15.01 16 24 H -0.27 -5.97 7.11 56.52 0.40 -5.57 16 25 H -0.28 -6.36 7.11 56.52 0.40 -5.95 16 26 H -0.28 -6.31 7.11 56.52 0.40 -5.91 16 27 C 0.10 2.13 8.19 -8.48 -0.07 2.06 16 28 C 0.04 0.48 7.71 -6.78 -0.05 0.43 16 29 H 0.08 0.49 8.14 -51.93 -0.42 0.07 16 30 H 0.04 0.26 6.42 -51.93 -0.33 -0.07 16 31 H 0.04 0.25 6.42 -51.93 -0.33 -0.08 16 32 H 0.07 0.49 8.14 -51.93 -0.42 0.06 16 33 H 0.06 0.39 6.54 -51.93 -0.34 0.05 16 34 H 0.07 0.55 8.14 -51.93 -0.42 0.12 16 35 H 0.07 0.53 6.35 -51.93 -0.33 0.20 16 36 H 0.06 0.38 6.54 -51.93 -0.34 0.04 16 37 H 0.07 0.48 8.14 -51.93 -0.42 0.06 16 38 H 0.08 0.95 8.14 -51.93 -0.42 0.53 16 39 H 0.08 0.82 8.14 -51.93 -0.42 0.40 16 40 H 0.04 0.44 5.51 -51.93 -0.29 0.16 16 41 H 0.08 0.88 7.76 -51.93 -0.40 0.47 16 42 H 0.07 0.79 7.42 -51.93 -0.39 0.40 16 43 H 0.11 1.93 5.97 -51.93 -0.31 1.62 16 44 H 0.10 1.08 8.14 -51.93 -0.42 0.65 16 45 H 0.10 1.30 8.14 -51.93 -0.42 0.87 16 46 H 0.06 1.25 8.14 -51.93 -0.42 0.82 16 47 H 0.07 1.52 8.10 -51.93 -0.42 1.10 16 48 H 0.07 1.76 8.11 -51.93 -0.42 1.34 16 49 H 0.26 6.82 8.31 -40.82 -0.34 6.49 16 50 H 0.07 1.47 8.10 -51.93 -0.42 1.05 16 51 H 0.06 1.19 8.14 -51.93 -0.42 0.77 16 52 H 0.09 1.02 8.14 -51.93 -0.42 0.60 16 53 H 0.07 0.71 7.91 -51.93 -0.41 0.30 16 LS Contribution 416.59 15.07 6.28 6.28 Total: -1.00 -29.34 416.59 -10.84 -40.18 By element: Atomic # 1 Polarization: 12.57 SS G_CDS: -9.35 Total: 3.22 kcal Atomic # 6 Polarization: -2.25 SS G_CDS: -0.98 Total: -3.23 kcal Atomic # 7 Polarization: -43.62 SS G_CDS: -1.49 Total: -45.11 kcal Atomic # 8 Polarization: -16.09 SS G_CDS: -0.26 Total: -16.35 kcal Atomic # 14 Polarization: 20.04 SS G_CDS: -5.02 Total: 15.01 kcal Total LS contribution 6.28 Total: 6.28 kcal Total: -29.34 -10.84 -40.18 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. ZINC001277131305.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 60.303 kcal (2) G-P(sol) polarization free energy of solvation -29.344 kcal (3) G-ENP(sol) elect.-nuc.-pol. free energy of system 30.959 kcal (4) G-CDS(sol) cavity-dispersion-solvent structure free energy -10.836 kcal (5) G-P-CDS(sol) = G-P(sol) + G-CDS(sol) = (2) + (4) -40.180 kcal (6) G-S(sol) free energy of system = (1) + (5) 20.123 kcal 1SCF WAS SPECIFIED, SO GEOMETRY OPTIMIZATION WAS NOT USED Total computer time = 0.79 seconds