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.53002 * 1 3 3 H 1.09007 * 109.46274 * 2 1 4 4 C 1.53039 * 109.45721 * 240.02905 * 2 1 3 5 5 C 1.53196 * 109.33748 * 178.59861 * 4 2 1 6 6 N 1.46995 * 108.73708 * 54.62190 * 5 4 2 7 7 C 1.34779 * 120.60954 * 126.39083 * 6 5 4 8 8 O 1.21511 * 119.99426 * 180.02562 * 7 6 5 9 9 O 1.34630 * 120.00346 * 0.02562 * 7 6 5 10 10 C 1.45200 * 116.99796 * 180.02562 * 9 7 6 11 11 C 1.53001 * 109.47436 * 59.99684 * 10 9 7 12 12 C 1.53003 * 109.46970 * 179.97438 * 10 9 7 13 13 C 1.53004 * 109.47330 * 299.99480 * 10 9 7 14 14 C 1.46925 * 118.75597 * 306.40484 * 6 5 4 15 15 H 1.08994 * 109.70376 * 173.46782 * 14 6 5 16 16 C 1.50701 * 109.58854 * 293.82776 * 14 6 5 17 17 O 1.21280 * 120.00310 * 19.28503 * 16 14 6 18 18 N 1.34775 * 119.99723 * 199.29061 * 16 14 6 19 19 C 1.46498 * 120.00364 * 359.97438 * 18 16 14 20 20 C 1.46504 * 119.99841 * 179.97438 * 18 16 14 21 21 C 1.50704 * 109.47160 * 89.99651 * 20 18 16 22 22 H 1.08000 * 119.99677 * 0.02562 * 21 20 18 23 23 C 1.37875 * 119.99804 * 179.97438 * 21 20 18 24 24 N 1.31509 * 120.00190 * 0.02562 * 23 21 20 25 25 Si 1.86801 * 119.99791 * 180.02562 * 23 21 20 26 26 H 1.48492 * 109.47274 * 59.99845 * 25 23 21 27 27 H 1.48497 * 109.47181 * 180.02562 * 25 23 21 28 28 H 1.48501 * 109.46736 * 300.00005 * 25 23 21 29 29 C 1.53033 * 109.46004 * 119.95865 * 2 1 3 30 30 H 1.08999 * 109.47143 * 179.97438 * 1 2 3 31 31 H 1.09008 * 109.46794 * 299.97314 * 1 2 3 32 32 H 1.08999 * 109.47146 * 59.96636 * 1 2 3 33 33 H 1.09002 * 109.49756 * 298.56584 * 4 2 1 34 34 H 1.09000 * 109.49883 * 58.63327 * 4 2 1 35 35 H 1.09001 * 109.58261 * 174.43824 * 5 4 2 36 36 H 1.08996 * 109.58319 * 294.86146 * 5 4 2 37 37 H 1.09000 * 109.46815 * 59.99950 * 11 10 9 38 38 H 1.09000 * 109.47277 * 179.97438 * 11 10 9 39 39 H 1.08994 * 109.47173 * 299.99950 * 11 10 9 40 40 H 1.09003 * 109.47421 * 59.99754 * 12 10 9 41 41 H 1.09000 * 109.47110 * 179.97438 * 12 10 9 42 42 H 1.09000 * 109.47220 * 300.00013 * 12 10 9 43 43 H 1.09003 * 109.46800 * 59.99710 * 13 10 9 44 44 H 1.08995 * 109.47741 * 180.02562 * 13 10 9 45 45 H 1.09000 * 109.46803 * 300.00665 * 13 10 9 46 46 H 1.08996 * 109.47273 * 89.99554 * 19 18 16 47 47 H 1.08996 * 109.47263 * 210.00483 * 19 18 16 48 48 H 1.09005 * 109.46479 * 330.00076 * 19 18 16 49 49 H 1.09000 * 109.47276 * 209.99731 * 20 18 16 50 50 H 1.08998 * 109.47167 * 329.99853 * 20 18 16 51 51 H 0.96999 * 120.00601 * 180.02562 * 24 23 21 52 52 H 1.08997 * 109.50670 * 301.33526 * 29 2 1 53 53 H 1.08998 * 109.50149 * 61.40361 * 29 2 1 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.0278 0.0000 4 6 2.0398 -0.7209 -1.2500 5 6 3.5714 -0.6887 -1.2650 6 7 4.0678 -1.2489 0.0001 7 6 4.9614 -2.2578 0.0000 8 8 5.3640 -2.7127 1.0524 9 8 5.4080 -2.7616 -1.1659 10 6 6.3689 -3.8475 -1.0897 11 6 5.7457 -5.0287 -0.3433 12 6 6.7565 -4.2855 -2.5036 13 6 7.6167 -3.3729 -0.3423 14 6 3.5714 -0.6884 1.2643 15 1 3.9417 -1.2822 2.0999 16 6 4.0455 0.7350 1.4066 17 8 4.4168 1.3510 0.4301 18 7 4.0570 1.3233 2.6191 19 6 3.6081 0.5794 3.7986 20 6 4.5184 2.7069 2.7576 21 6 3.3504 3.6448 2.5926 22 1 2.3631 3.2498 2.4039 23 6 3.5421 5.0072 2.6836 24 7 4.7442 5.4881 2.9139 25 14 2.0942 6.1697 2.4798 26 1 1.0765 5.8803 3.5217 27 1 2.5619 7.5721 2.6201 28 1 1.4934 5.9810 1.1350 29 6 2.0398 -0.7206 1.2501 30 1 -0.3633 -1.0277 -0.0005 31 1 -0.3633 0.5135 0.8903 32 1 -0.3633 0.5143 -0.8897 33 1 1.6977 -1.7557 -1.2382 34 1 1.6567 -0.2212 -2.1398 35 1 3.9411 -1.2846 -2.0995 36 1 3.9142 0.3408 -1.3683 37 1 4.8568 -5.3668 -0.8758 38 1 6.4669 -5.8440 -0.2865 39 1 5.4692 -4.7170 0.6639 40 1 7.2003 -3.4441 -3.0358 41 1 7.4782 -5.1004 -2.4465 42 1 5.8674 -4.6237 -3.0357 43 1 7.3402 -3.0610 0.6649 44 1 8.3383 -4.1878 -0.2855 45 1 8.0606 -2.5314 -0.8741 46 1 4.4559 0.0596 4.2447 47 1 3.1833 1.2720 4.5251 48 1 2.8509 -0.1471 3.5034 49 1 4.9593 2.8458 3.7447 50 1 5.2655 2.9190 1.9928 51 1 4.8791 6.4465 2.9783 52 1 1.6573 -0.2205 2.1399 53 1 1.6974 -1.7553 1.2388 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 ZINC000495491885.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 = -120.196290 kcal Electronic energy + Delta-G solvation = -32696.435989 eV Core-core repulsion = 28568.261628 eV Total energy + Delta-G solvation = -4128.174361 eV No. of doubly occupied orbitals = 66 Molecular weight (most abundant/longest-lived isotopes) = 340.220 amu Computer time = 0.62 seconds Orbital eigenvalues (eV) -41.03721 -39.74609 -37.59626 -36.60475 -34.89870 -34.36514 -33.43916 -32.71309 -31.23780 -28.49280 -27.61392 -27.53291 -27.45589 -25.93113 -24.39442 -23.32262 -23.01094 -21.57335 -20.42196 -19.42378 -18.30521 -17.72664 -16.62790 -16.52819 -16.27843 -15.80663 -15.53855 -15.15325 -14.86100 -14.57255 -14.39684 -14.09689 -13.94948 -13.81518 -13.55819 -13.41699 -13.01282 -12.81854 -12.60145 -12.32592 -12.21811 -12.18511 -12.01941 -11.99737 -11.93529 -11.85042 -11.64142 -11.53071 -11.37628 -11.23713 -11.17065 -10.90602 -10.80050 -10.57527 -10.44220 -10.39134 -10.36446 -10.14782 -9.89552 -9.65426 -9.35553 -8.44135 -8.21859 -7.65764 -6.04728 -4.09827 2.12507 2.91652 3.30292 3.32002 3.35496 3.39435 3.81298 4.03207 4.43117 4.61281 4.65487 4.77769 4.79403 4.81021 4.86593 5.05970 5.15598 5.29462 5.29802 5.32643 5.41769 5.52238 5.53594 5.55518 5.56550 5.64574 5.66547 5.69727 5.75465 5.76506 5.82908 5.85172 5.90394 5.95821 5.97223 6.11859 6.18704 6.21146 6.23435 6.27409 6.29751 6.48500 6.52992 6.58616 7.04061 7.43369 7.48727 7.52748 7.91025 8.22169 8.35616 8.44063 8.87975 9.05047 9.55058 11.03339 Molecular weight = 340.22amu Principal moments of inertia in cm(-1) A = 0.014803 B = 0.003411 C = 0.003119 Principal moments of inertia in units of 10**(-40)*gram-cm**2 A = 1890.991721 B = 8206.449980 C = 8976.170266 Net atomic charges, atomic populations, and dipole contributions calculated with CM2 Atom NO. Type Charge No. of electrons 1 C -0.142 4.142 2 C -0.100 4.100 3 H 0.087 0.913 4 C -0.123 4.123 5 C 0.104 3.896 6 N -0.570 5.570 7 C 0.657 3.343 8 O -0.575 6.575 9 O -0.360 6.360 10 C 0.135 3.865 11 C -0.182 4.182 12 C -0.140 4.140 13 C -0.184 4.184 14 C 0.137 3.863 15 H 0.118 0.882 16 C 0.501 3.499 17 O -0.532 6.532 18 N -0.597 5.597 19 C 0.076 3.924 20 C 0.247 3.753 21 C -0.525 4.525 22 H 0.057 0.943 23 C 0.062 3.938 24 N -0.895 5.895 25 Si 0.897 3.103 26 H -0.277 1.277 27 H -0.287 1.287 28 H -0.276 1.276 29 C -0.116 4.116 30 H 0.051 0.949 31 H 0.056 0.944 32 H 0.055 0.945 33 H 0.061 0.939 34 H 0.075 0.925 35 H 0.085 0.915 36 H 0.098 0.902 37 H 0.057 0.943 38 H 0.068 0.932 39 H 0.088 0.912 40 H 0.067 0.933 41 H 0.078 0.922 42 H 0.065 0.935 43 H 0.088 0.912 44 H 0.068 0.932 45 H 0.060 0.940 46 H 0.047 0.953 47 H 0.077 0.923 48 H 0.059 0.941 49 H 0.032 0.968 50 H 0.043 0.957 51 H 0.264 0.736 52 H 0.086 0.914 53 H 0.074 0.926 For the charges calculated by CM2: Dipole moment (debyes) X Y Z Total from point charges -2.705 -20.251 -7.097 21.629 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.200 4.200 2 C -0.118 4.118 3 H 0.105 0.895 4 C -0.161 4.161 5 C -0.018 4.018 6 N -0.309 5.309 7 C 0.412 3.588 8 O -0.465 6.465 9 O -0.275 6.275 10 C 0.098 3.902 11 C -0.239 4.239 12 C -0.197 4.197 13 C -0.241 4.241 14 C 0.031 3.969 15 H 0.136 0.864 16 C 0.286 3.714 17 O -0.409 6.409 18 N -0.329 5.329 19 C -0.068 4.068 20 C 0.128 3.872 21 C -0.564 4.564 22 H 0.075 0.925 23 C -0.141 4.141 24 N -0.533 5.533 25 Si 0.725 3.275 26 H -0.202 1.202 27 H -0.212 1.212 28 H -0.202 1.202 29 C -0.154 4.154 30 H 0.070 0.930 31 H 0.076 0.924 32 H 0.074 0.926 33 H 0.080 0.920 34 H 0.094 0.906 35 H 0.103 0.897 36 H 0.116 0.884 37 H 0.077 0.923 38 H 0.086 0.914 39 H 0.107 0.893 40 H 0.086 0.914 41 H 0.097 0.903 42 H 0.084 0.916 43 H 0.107 0.893 44 H 0.087 0.913 45 H 0.079 0.921 46 H 0.066 0.934 47 H 0.095 0.905 48 H 0.078 0.922 49 H 0.050 0.950 50 H 0.061 0.939 51 H 0.073 0.927 52 H 0.105 0.895 53 H 0.092 0.908 Dipole moment (debyes) X Y Z Total from point charges -2.207 -19.927 -7.297 21.336 hybrid contribution -0.703 0.726 0.065 1.013 sum -2.910 -19.202 -7.232 20.724 Atomic orbital electron populations 1.21782 0.94884 1.01266 1.02030 1.21177 0.95730 1.01014 0.93910 0.89483 1.21929 0.95192 1.00573 0.98431 1.21981 0.93893 1.01719 0.84187 1.46536 1.42950 1.33504 1.07862 1.17841 0.78813 0.79797 0.82313 1.90859 1.61060 1.62878 1.31716 1.86295 1.58282 1.50563 1.32334 1.22040 0.87045 0.84251 0.96871 1.22505 1.01106 0.97149 1.03130 1.21846 1.02054 1.01846 0.93928 1.22514 0.96428 1.01890 1.03273 1.21176 0.95143 0.93976 0.86634 0.86421 1.20442 0.79384 0.87648 0.83937 1.90610 1.48010 1.65138 1.37117 1.48135 1.63177 1.15890 1.05727 1.21788 0.98844 0.96731 0.89466 1.19293 0.95732 0.74361 0.97852 1.21367 0.94603 0.91037 1.49405 0.92468 1.24931 0.95905 0.98117 0.95106 1.69921 1.20916 1.13414 1.49086 0.86574 0.82472 0.79737 0.78729 1.20245 1.21213 1.20179 1.21963 0.92042 1.01720 0.99724 0.93039 0.92442 0.92602 0.92019 0.90623 0.89699 0.88395 0.92347 0.91351 0.89309 0.91434 0.90302 0.91588 0.89276 0.91346 0.92096 0.93439 0.90455 0.92236 0.94965 0.93909 0.92657 0.89545 0.90766 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.14 -1.62 9.19 37.16 0.34 -1.28 16 2 C -0.10 -1.37 2.52 -90.58 -0.23 -1.60 16 3 H 0.09 1.52 7.10 -51.92 -0.37 1.16 16 4 C -0.12 -1.51 5.36 -26.61 -0.14 -1.65 16 5 C 0.10 1.52 5.95 -3.67 -0.02 1.50 16 6 N -0.57 -9.15 2.54 -172.23 -0.44 -9.59 16 7 C 0.66 10.68 7.74 54.06 0.42 11.10 16 8 O -0.57 -10.05 11.30 -19.28 -0.22 -10.26 16 9 O -0.36 -5.13 8.53 -40.30 -0.34 -5.48 16 10 C 0.13 1.55 1.13 -90.62 -0.10 1.45 16 11 C -0.18 -1.92 8.37 37.16 0.31 -1.61 16 12 C -0.14 -1.18 8.85 37.16 0.33 -0.85 16 13 C -0.18 -2.19 8.37 37.16 0.31 -1.88 16 14 C 0.14 2.24 2.15 -68.75 -0.15 2.09 16 15 H 0.12 1.89 6.10 -51.93 -0.32 1.57 16 16 C 0.50 10.42 5.98 -10.99 -0.07 10.36 16 17 O -0.53 -12.94 13.54 5.56 0.08 -12.86 16 18 N -0.60 -12.70 2.86 -182.95 -0.52 -13.22 16 19 C 0.08 1.22 9.86 59.85 0.59 1.81 16 20 C 0.25 6.45 5.16 -5.19 -0.03 6.42 16 21 C -0.53 -14.63 9.39 -34.44 -0.32 -14.95 16 22 H 0.06 1.57 7.81 -52.49 -0.41 1.16 16 23 C 0.06 1.74 5.46 -17.82 -0.10 1.64 16 24 N -0.89 -26.51 13.29 55.43 0.74 -25.77 16 25 Si 0.90 21.59 29.54 -169.99 -5.02 16.57 16 26 H -0.28 -6.45 7.11 56.52 0.40 -6.05 16 27 H -0.29 -6.82 7.11 56.52 0.40 -6.42 16 28 H -0.28 -6.64 7.11 56.52 0.40 -6.24 16 29 C -0.12 -1.56 4.49 -26.61 -0.12 -1.68 16 30 H 0.05 0.50 8.14 -51.93 -0.42 0.08 16 31 H 0.06 0.67 8.14 -51.92 -0.42 0.24 16 32 H 0.05 0.60 8.14 -51.93 -0.42 0.18 16 33 H 0.06 0.69 8.14 -51.93 -0.42 0.26 16 34 H 0.08 0.83 8.14 -51.93 -0.42 0.41 16 35 H 0.08 1.10 6.93 -51.93 -0.36 0.74 16 36 H 0.10 1.71 6.16 -51.93 -0.32 1.39 16 37 H 0.06 0.57 8.14 -51.93 -0.42 0.15 16 38 H 0.07 0.58 8.14 -51.93 -0.42 0.15 16 39 H 0.09 1.15 5.88 -51.93 -0.31 0.85 16 40 H 0.07 0.58 8.14 -51.93 -0.42 0.15 16 41 H 0.08 0.51 8.14 -51.93 -0.42 0.08 16 42 H 0.07 0.52 8.14 -51.93 -0.42 0.10 16 43 H 0.09 1.30 5.88 -51.93 -0.31 1.00 16 44 H 0.07 0.64 8.14 -51.93 -0.42 0.22 16 45 H 0.06 0.74 8.14 -51.93 -0.42 0.31 16 46 H 0.05 0.70 8.14 -51.93 -0.42 0.27 16 47 H 0.08 1.25 8.07 -51.93 -0.42 0.83 16 48 H 0.06 0.79 4.84 -51.93 -0.25 0.54 16 49 H 0.03 0.86 8.07 -51.93 -0.42 0.44 16 50 H 0.04 1.27 7.42 -51.93 -0.39 0.88 16 51 H 0.26 7.41 8.31 -40.82 -0.34 7.07 16 52 H 0.09 1.14 6.44 -51.93 -0.33 0.80 16 53 H 0.07 0.87 8.14 -51.93 -0.42 0.44 16 LS Contribution 405.89 15.07 6.12 6.12 Total: -1.00 -33.01 405.89 -7.84 -40.85 By element: Atomic # 1 Polarization: 12.03 SS G_CDS: -9.25 Total: 2.78 kcal Atomic # 6 Polarization: 9.84 SS G_CDS: 1.03 Total: 10.87 kcal Atomic # 7 Polarization: -48.36 SS G_CDS: -0.22 Total: -48.58 kcal Atomic # 8 Polarization: -28.12 SS G_CDS: -0.49 Total: -28.60 kcal Atomic # 14 Polarization: 21.59 SS G_CDS: -5.02 Total: 16.57 kcal Total LS contribution 6.12 Total: 6.12 kcal Total: -33.01 -7.84 -40.85 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. ZINC000495491885.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 -79.347 kcal (2) G-P(sol) polarization free energy of solvation -33.011 kcal (3) G-ENP(sol) elect.-nuc.-pol. free energy of system -112.358 kcal (4) G-CDS(sol) cavity-dispersion-solvent structure free energy -7.839 kcal (5) G-P-CDS(sol) = G-P(sol) + G-CDS(sol) = (2) + (4) -40.849 kcal (6) G-S(sol) free energy of system = (1) + (5) -120.196 kcal 1SCF WAS SPECIFIED, SO GEOMETRY OPTIMIZATION WAS NOT USED Total computer time = 0.62 seconds