Wall clock time and date at job start Sat Apr 11 2020 06:10: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= -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.52995 * 1 3 3 H 1.09001 * 109.46687 * 2 1 4 4 C 1.50699 * 109.47759 * 119.99888 * 2 1 3 5 5 H 1.07998 * 120.00236 * 119.99901 * 4 2 1 6 6 C 1.37881 * 119.99524 * 300.00561 * 4 2 1 7 7 N 1.31504 * 120.00621 * 359.97438 * 6 4 2 8 8 Si 1.86805 * 119.99439 * 179.97438 * 6 4 2 9 9 H 1.48500 * 109.47059 * 90.00241 * 8 6 4 10 10 H 1.48501 * 109.47135 * 209.99972 * 8 6 4 11 11 H 1.48494 * 109.47464 * 329.99713 * 8 6 4 12 12 N 1.46501 * 109.47146 * 240.00156 * 2 1 3 13 13 C 1.46497 * 120.00077 * 299.26328 * 12 2 1 14 14 C 1.34781 * 119.99952 * 119.54278 * 12 2 1 15 15 O 1.21276 * 120.00045 * 175.10469 * 14 12 2 16 16 C 1.50701 * 120.00103 * 355.10361 * 14 12 2 17 17 H 1.09008 * 109.47025 * 311.27881 * 16 14 12 18 18 C 1.52999 * 109.47314 * 71.27494 * 16 14 12 19 19 C 1.52995 * 109.47460 * 191.27888 * 16 14 12 20 20 C 1.53939 * 113.74751 * 158.73936 * 19 16 14 21 21 N 1.47811 * 86.04100 * 222.08611 * 20 19 16 22 22 C 1.34775 * 134.45894 * 204.64498 * 21 20 19 23 23 O 1.21508 * 120.00032 * 179.86012 * 22 21 20 24 24 O 1.34638 * 119.99687 * 359.86515 * 22 21 20 25 25 C 1.45207 * 116.99748 * 179.97438 * 24 22 21 26 26 C 1.52996 * 109.46958 * 59.99764 * 25 24 22 27 27 C 1.53002 * 109.46676 * 179.97438 * 25 24 22 28 28 C 1.52997 * 109.46845 * 299.99834 * 25 24 22 29 29 C 1.47578 * 91.04863 * 24.54553 * 21 20 19 30 30 H 1.08991 * 109.47456 * 60.00116 * 1 2 3 31 31 H 1.09001 * 109.46687 * 179.97438 * 1 2 3 32 32 H 1.09006 * 109.47115 * 299.99478 * 1 2 3 33 33 H 0.97003 * 120.00577 * 180.02562 * 7 6 4 34 34 H 1.08993 * 109.47203 * 264.99557 * 13 12 2 35 35 H 1.09005 * 109.46763 * 24.99575 * 13 12 2 36 36 H 1.09003 * 109.47079 * 144.98887 * 13 12 2 37 37 H 1.09008 * 109.47267 * 293.21445 * 18 16 14 38 38 H 1.09000 * 109.47053 * 53.21150 * 18 16 14 39 39 H 1.08993 * 109.47465 * 173.21401 * 18 16 14 40 40 H 1.09003 * 112.94448 * 290.33277 * 19 16 14 41 41 H 1.09009 * 113.76611 * 336.34758 * 20 19 16 42 42 H 1.08997 * 113.77233 * 107.73942 * 20 19 16 43 43 H 1.09003 * 109.46991 * 60.00109 * 26 25 24 44 44 H 1.08998 * 109.46913 * 179.97438 * 26 25 24 45 45 H 1.08994 * 109.47150 * 300.00220 * 26 25 24 46 46 H 1.08997 * 109.47224 * 60.00567 * 27 25 24 47 47 H 1.09008 * 109.46796 * 179.97438 * 27 25 24 48 48 H 1.08992 * 109.47486 * 299.99619 * 27 25 24 49 49 H 1.08991 * 109.47511 * 60.00189 * 28 25 24 50 50 H 1.09008 * 109.46939 * 180.02562 * 28 25 24 51 51 H 1.09001 * 109.47135 * 299.99399 * 28 25 24 52 52 H 1.09000 * 113.76950 * 221.15441 * 29 21 20 53 53 H 1.09000 * 113.86939 * 89.77168 * 29 21 20 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.8932 1.0277 0.0000 4 6 2.0324 -0.7104 1.2304 5 1 2.6534 -1.5884 1.1313 6 6 1.6993 -0.2393 2.4827 7 7 0.9435 0.8301 2.6034 8 14 2.3216 -1.1202 4.0079 9 1 3.6298 -0.5456 4.4123 10 1 1.3434 -0.9493 5.1120 11 1 2.4862 -2.5661 3.7120 12 7 2.0183 -0.6906 -1.1962 13 6 1.6778 -2.0977 -1.4200 14 6 2.7856 -0.0329 -2.0880 15 8 3.2631 -0.6266 -3.0315 16 6 3.0500 1.4404 -1.9128 17 1 2.1093 1.9588 -1.7270 18 6 3.9919 1.6519 -0.7258 19 6 3.6968 1.9970 -3.1826 20 6 3.5444 3.5208 -3.3392 21 7 3.2605 3.2430 -4.7630 22 6 3.3718 3.9523 -5.9035 23 8 3.0559 3.4507 -6.9642 24 8 3.8331 5.2167 -5.8677 25 6 3.9270 5.9203 -7.1344 26 6 2.5397 6.0068 -7.7737 27 6 4.4666 7.3319 -6.8948 28 6 4.8731 5.1638 -8.0690 29 6 2.8088 1.8770 -4.4344 30 1 -0.3634 0.5138 -0.8899 31 1 -0.3633 -1.0277 -0.0005 32 1 -0.3634 0.5138 0.8901 33 1 0.7095 1.1618 3.4844 34 1 2.4818 -2.7303 -1.0441 35 1 0.7526 -2.3356 -0.8949 36 1 1.5453 -2.2742 -2.4874 37 1 3.4847 1.3665 0.1959 38 1 4.8830 1.0377 -0.8552 39 1 4.2790 2.7020 -0.6724 40 1 4.7160 1.6404 -3.3318 41 1 2.6959 3.9357 -2.7949 42 1 4.4688 4.0742 -3.1741 43 1 1.8657 6.5457 -7.1078 44 1 2.6102 6.5354 -8.7244 45 1 2.1550 5.0014 -7.9442 46 1 5.4549 7.2706 -6.4392 47 1 4.5375 7.8599 -7.8458 48 1 3.7923 7.8705 -6.2292 49 1 4.4883 4.1585 -8.2395 50 1 4.9439 5.6923 -9.0198 51 1 5.8615 5.1023 -7.6135 52 1 3.1471 1.1203 -5.1422 53 1 1.7440 1.8060 -4.2125 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 ZINC000495488702.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:10:28 Heat of formation + Delta-G solvation = -82.095340 kcal Electronic energy + Delta-G solvation = -31319.531364 eV Core-core repulsion = 27193.009185 eV Total energy + Delta-G solvation = -4126.522180 eV No. of doubly occupied orbitals = 66 Molecular weight (most abundant/longest-lived isotopes) = 340.220 amu Computer time = 0.67 seconds Orbital eigenvalues (eV) -40.91904 -39.94494 -38.56837 -36.55439 -35.60164 -34.17999 -33.14817 -32.69523 -29.22542 -28.88203 -27.61113 -27.55275 -26.61420 -26.03002 -25.39271 -24.04959 -22.11686 -21.23160 -20.39027 -19.25529 -18.34734 -17.67012 -16.79453 -16.68511 -16.44740 -16.16785 -15.37591 -15.21632 -14.79325 -14.65524 -14.42991 -14.20914 -13.99545 -13.49485 -13.36555 -12.96119 -12.82752 -12.70744 -12.49477 -12.43729 -12.38276 -12.07612 -12.05920 -12.01623 -11.70818 -11.65528 -11.51810 -11.46240 -11.37068 -11.16454 -11.03944 -10.94953 -10.79177 -10.68254 -10.49843 -10.47393 -10.24219 -10.16852 -9.95500 -9.36592 -8.93713 -8.76231 -8.58152 -7.66968 -6.23474 -4.31689 1.97900 2.87351 3.18217 3.26256 3.27830 3.35470 3.68572 3.87036 4.26780 4.49866 4.53404 4.57146 4.71363 4.91123 4.95158 5.07888 5.19770 5.25184 5.25636 5.33083 5.39650 5.47865 5.48128 5.54806 5.56173 5.64002 5.66913 5.70815 5.77034 5.83814 5.87647 6.05293 6.11301 6.15888 6.19800 6.22978 6.26327 6.41076 6.43912 6.45252 6.53443 6.68263 6.74398 6.94334 7.15810 7.36585 7.61930 7.69275 7.81130 8.05956 8.34030 8.41358 8.71078 8.96803 9.37129 10.80891 Molecular weight = 340.22amu Principal moments of inertia in cm(-1) A = 0.021183 B = 0.002382 C = 0.002283 Principal moments of inertia in units of 10**(-40)*gram-cm**2 A = 1321.518615 B =11751.742164 C =12260.446075 Net atomic charges, atomic populations, and dipole contributions calculated with CM2 Atom NO. Type Charge No. of electrons 1 C -0.154 4.154 2 C 0.246 3.754 3 H 0.065 0.935 4 C -0.527 4.527 5 H 0.061 0.939 6 C 0.068 3.932 7 N -0.888 5.888 8 Si 0.896 3.104 9 H -0.279 1.279 10 H -0.286 1.286 11 H -0.274 1.274 12 N -0.590 5.590 13 C 0.099 3.901 14 C 0.513 3.487 15 O -0.572 6.572 16 C -0.096 4.096 17 H 0.097 0.903 18 C -0.151 4.151 19 C -0.095 4.095 20 C 0.111 3.889 21 N -0.633 5.633 22 C 0.672 3.328 23 O -0.572 6.572 24 O -0.363 6.363 25 C 0.135 3.865 26 C -0.183 4.183 27 C -0.140 4.140 28 C -0.183 4.183 29 C 0.145 3.855 30 H 0.012 0.988 31 H 0.033 0.967 32 H 0.095 0.905 33 H 0.265 0.735 34 H 0.044 0.956 35 H 0.070 0.930 36 H 0.054 0.946 37 H 0.097 0.903 38 H 0.055 0.945 39 H 0.058 0.942 40 H 0.108 0.892 41 H 0.091 0.909 42 H 0.085 0.915 43 H 0.059 0.941 44 H 0.068 0.932 45 H 0.087 0.913 46 H 0.066 0.934 47 H 0.079 0.921 48 H 0.066 0.934 49 H 0.088 0.912 50 H 0.068 0.932 51 H 0.059 0.941 52 H 0.090 0.910 53 H 0.085 0.915 For the charges calculated by CM2: Dipole moment (debyes) X Y Z Total from point charges 4.253 10.649 -12.543 16.994 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.212 4.212 2 C 0.145 3.855 3 H 0.082 0.918 4 C -0.566 4.566 5 H 0.079 0.921 6 C -0.135 4.135 7 N -0.526 5.526 8 Si 0.724 3.276 9 H -0.205 1.205 10 H -0.211 1.211 11 H -0.199 1.199 12 N -0.320 5.320 13 C -0.044 4.044 14 C 0.300 3.700 15 O -0.454 6.454 16 C -0.118 4.118 17 H 0.115 0.885 18 C -0.209 4.209 19 C -0.115 4.115 20 C -0.012 4.012 21 N -0.376 5.376 22 C 0.428 3.572 23 O -0.462 6.462 24 O -0.278 6.278 25 C 0.098 3.902 26 C -0.240 4.240 27 C -0.197 4.197 28 C -0.240 4.240 29 C 0.025 3.975 30 H 0.031 0.969 31 H 0.052 0.948 32 H 0.113 0.887 33 H 0.075 0.925 34 H 0.063 0.937 35 H 0.089 0.911 36 H 0.072 0.928 37 H 0.115 0.885 38 H 0.074 0.926 39 H 0.077 0.923 40 H 0.127 0.873 41 H 0.109 0.891 42 H 0.104 0.896 43 H 0.078 0.922 44 H 0.087 0.913 45 H 0.106 0.894 46 H 0.085 0.915 47 H 0.098 0.902 48 H 0.085 0.915 49 H 0.107 0.893 50 H 0.087 0.913 51 H 0.078 0.922 52 H 0.109 0.891 53 H 0.103 0.897 Dipole moment (debyes) X Y Z Total from point charges 4.320 10.480 -13.431 17.575 hybrid contribution 0.128 -0.415 0.980 1.072 sum 4.448 10.065 -12.451 16.617 Atomic orbital electron populations 1.22159 0.98264 0.99937 1.00852 1.19281 0.91392 0.93921 0.80875 0.91804 1.21321 1.31890 1.13346 0.90003 0.92053 1.24863 0.94723 0.94954 0.98955 1.69840 1.36099 1.22619 1.24027 0.86520 0.78651 0.80129 0.82316 1.20521 1.21131 1.19900 1.48250 1.49073 1.11469 1.23228 1.21335 1.00471 0.82561 1.00080 1.20370 0.78551 0.88992 0.82056 1.90601 1.52007 1.69710 1.33061 1.21214 1.00978 0.93972 0.95615 0.88503 1.21956 0.99156 1.01482 0.98274 1.22937 1.01125 0.93944 0.93465 1.23429 1.01801 0.96940 0.78989 1.48880 1.67398 1.16519 1.04759 1.17540 0.76686 0.79817 0.83155 1.90881 1.55683 1.71315 1.28352 1.86315 1.79273 1.25234 1.36973 1.22045 0.96408 0.91211 0.80560 1.22510 0.95777 1.04882 1.00846 1.21847 1.01937 0.93007 1.02883 1.22514 0.99132 1.03273 0.99101 1.23221 0.99152 0.80381 0.94761 0.96884 0.94826 0.88685 0.92534 0.93750 0.91131 0.92769 0.88470 0.92617 0.92270 0.87341 0.89127 0.89619 0.92219 0.91260 0.89389 0.91504 0.90248 0.91486 0.89350 0.91285 0.92215 0.89148 0.89710 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.24 8.97 37.16 0.33 -2.91 16 2 C 0.25 5.26 1.85 -69.08 -0.13 5.13 16 3 H 0.06 1.37 4.07 -51.93 -0.21 1.16 16 4 C -0.53 -12.62 8.08 -34.44 -0.28 -12.90 16 5 H 0.06 1.52 7.62 -52.49 -0.40 1.12 16 6 C 0.07 1.65 5.30 -17.82 -0.09 1.56 16 7 N -0.89 -22.13 11.32 55.43 0.63 -21.50 16 8 Si 0.90 19.46 29.54 -169.99 -5.02 14.44 16 9 H -0.28 -6.07 7.11 56.52 0.40 -5.67 16 10 H -0.29 -6.22 7.11 56.52 0.40 -5.82 16 11 H -0.27 -6.00 7.11 56.52 0.40 -5.60 16 12 N -0.59 -11.73 2.75 -175.27 -0.48 -12.21 16 13 C 0.10 1.83 9.51 59.85 0.57 2.40 16 14 C 0.51 9.64 6.66 -10.98 -0.07 9.57 16 15 O -0.57 -11.81 15.18 5.56 0.08 -11.73 16 16 C -0.10 -1.46 1.60 -91.77 -0.15 -1.61 16 17 H 0.10 1.43 7.26 -51.92 -0.38 1.06 16 18 C -0.15 -2.30 8.78 37.16 0.33 -1.97 16 19 C -0.10 -1.28 2.76 -89.70 -0.25 -1.53 16 20 C 0.11 1.29 8.38 -2.84 -0.02 1.27 16 21 N -0.63 -9.02 3.61 -183.97 -0.66 -9.69 16 22 C 0.67 10.35 8.22 54.05 0.44 10.79 16 23 O -0.57 -9.92 12.61 -19.28 -0.24 -10.16 16 24 O -0.36 -4.92 9.87 -40.32 -0.40 -5.32 16 25 C 0.13 1.46 1.13 -90.62 -0.10 1.36 16 26 C -0.18 -1.92 8.37 37.16 0.31 -1.61 16 27 C -0.14 -1.09 8.85 37.16 0.33 -0.76 16 28 C -0.18 -1.93 8.37 37.16 0.31 -1.62 16 29 C 0.15 2.10 7.34 -2.97 -0.02 2.08 16 30 H 0.01 0.24 8.14 -51.93 -0.42 -0.19 16 31 H 0.03 0.69 6.69 -51.93 -0.35 0.34 16 32 H 0.09 2.20 6.07 -51.93 -0.32 1.89 16 33 H 0.26 6.38 8.31 -40.82 -0.34 6.04 16 34 H 0.04 0.82 8.14 -51.93 -0.42 0.40 16 35 H 0.07 1.26 6.48 -51.93 -0.34 0.92 16 36 H 0.05 0.95 7.66 -51.93 -0.40 0.55 16 37 H 0.10 1.79 5.42 -51.92 -0.28 1.50 16 38 H 0.05 0.85 8.14 -51.93 -0.42 0.43 16 39 H 0.06 0.73 8.08 -51.93 -0.42 0.31 16 40 H 0.11 1.50 8.14 -51.93 -0.42 1.08 16 41 H 0.09 0.96 7.99 -51.92 -0.41 0.54 16 42 H 0.09 0.88 8.14 -51.93 -0.42 0.46 16 43 H 0.06 0.60 8.14 -51.93 -0.42 0.18 16 44 H 0.07 0.57 8.14 -51.93 -0.42 0.15 16 45 H 0.09 1.16 5.88 -51.93 -0.31 0.85 16 46 H 0.07 0.51 8.14 -51.93 -0.42 0.08 16 47 H 0.08 0.46 8.14 -51.92 -0.42 0.04 16 48 H 0.07 0.51 8.14 -51.93 -0.42 0.09 16 49 H 0.09 1.17 5.88 -51.93 -0.31 0.86 16 50 H 0.07 0.57 8.14 -51.92 -0.42 0.14 16 51 H 0.06 0.60 8.14 -51.93 -0.42 0.18 16 52 H 0.09 1.42 8.14 -51.93 -0.42 1.00 16 53 H 0.08 1.19 8.10 -51.93 -0.42 0.77 16 LS Contribution 411.75 15.07 6.21 6.21 Total: -1.00 -28.28 411.75 -7.55 -35.82 By element: Atomic # 1 Polarization: 14.06 SS G_CDS: -9.16 Total: 4.89 kcal Atomic # 6 Polarization: 7.74 SS G_CDS: 1.51 Total: 9.25 kcal Atomic # 7 Polarization: -42.88 SS G_CDS: -0.52 Total: -43.40 kcal Atomic # 8 Polarization: -26.65 SS G_CDS: -0.56 Total: -27.21 kcal Atomic # 14 Polarization: 19.46 SS G_CDS: -5.02 Total: 14.44 kcal Total LS contribution 6.21 Total: 6.21 kcal Total: -28.28 -7.55 -35.82 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. ZINC000495488702.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 -46.273 kcal (2) G-P(sol) polarization free energy of solvation -28.277 kcal (3) G-ENP(sol) elect.-nuc.-pol. free energy of system -74.549 kcal (4) G-CDS(sol) cavity-dispersion-solvent structure free energy -7.546 kcal (5) G-P-CDS(sol) = G-P(sol) + G-CDS(sol) = (2) + (4) -35.823 kcal (6) G-S(sol) free energy of system = (1) + (5) -82.095 kcal 1SCF WAS SPECIFIED, SO GEOMETRY OPTIMIZATION WAS NOT USED Total computer time = 0.67 seconds