Wall clock time and date at job start Wed Apr 1 2020 00:25:14 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.52999 * 1 3 3 H 1.08996 * 109.47116 * 2 1 4 4 C 1.50700 * 109.47379 * 239.99986 * 2 1 3 5 5 O 1.21291 * 119.99707 * 270.54447 * 4 2 1 6 6 N 1.34772 * 119.99918 * 90.54318 * 4 2 1 7 7 C 1.46494 * 120.00280 * 354.76875 * 6 4 2 8 8 C 1.46506 * 119.99709 * 174.77198 * 6 4 2 9 9 C 1.50699 * 109.47151 * 89.99642 * 8 6 4 10 10 H 1.07998 * 119.99884 * 0.02562 * 9 8 6 11 11 C 1.37876 * 120.00038 * 179.97438 * 9 8 6 12 12 N 1.31508 * 120.00217 * 180.02562 * 11 9 8 13 13 Si 1.86801 * 119.99776 * 0.02562 * 11 9 8 14 14 H 1.48502 * 109.47180 * 89.99781 * 13 11 9 15 15 H 1.48505 * 109.47105 * 209.99907 * 13 11 9 16 16 H 1.48502 * 109.47421 * 329.99874 * 13 11 9 17 17 N 1.46908 * 109.47039 * 120.00232 * 2 1 3 18 18 C 1.47009 * 110.99862 * 178.29694 * 17 2 1 19 19 C 1.53098 * 109.38248 * 182.82734 * 18 17 2 20 20 H 1.09000 * 109.51698 * 62.54445 * 19 18 17 21 21 C 1.53001 * 109.52246 * 182.71665 * 19 18 17 22 22 O 1.43151 * 109.19464 * 302.63072 * 19 18 17 23 23 C 1.42295 * 113.81277 * 57.52051 * 22 19 18 24 24 C 1.53995 * 110.32963 * 65.46343 * 23 22 19 25 25 C 1.54610 * 103.87927 * 118.89351 * 24 23 22 26 26 O 1.44474 * 104.81564 * 336.10213 * 25 24 23 27 27 C 1.44548 * 105.36948 * 40.16653 * 26 25 24 28 28 H 1.08994 * 110.48951 * 78.85703 * 27 26 25 29 29 C 1.53001 * 110.39091 * 201.19709 * 27 26 25 30 30 C 1.47140 * 110.99587 * 301.82758 * 17 2 1 31 31 H 1.09006 * 109.47067 * 185.70467 * 1 2 3 32 32 H 1.08993 * 109.47301 * 305.70584 * 1 2 3 33 33 H 1.08998 * 109.47272 * 65.70835 * 1 2 3 34 34 H 1.08998 * 109.46906 * 94.49633 * 7 6 4 35 35 H 1.08994 * 109.47215 * 214.49514 * 7 6 4 36 36 H 1.08997 * 109.47024 * 334.49872 * 7 6 4 37 37 H 1.08996 * 109.46789 * 210.00016 * 8 6 4 38 38 H 1.09001 * 109.47070 * 330.00064 * 8 6 4 39 39 H 0.97002 * 120.00076 * 180.02562 * 12 11 9 40 40 H 1.09003 * 109.48818 * 62.85709 * 18 17 2 41 41 H 1.08994 * 109.50354 * 302.81151 * 18 17 2 42 42 H 1.08995 * 109.47527 * 179.71332 * 21 19 18 43 43 H 1.08996 * 109.46829 * 299.71741 * 21 19 18 44 44 H 1.09008 * 109.46955 * 59.71197 * 21 19 18 45 45 H 1.09000 * 110.54997 * 237.49393 * 24 23 22 46 46 H 1.09000 * 110.54939 * 0.29067 * 24 23 22 47 47 H 1.08995 * 110.36876 * 94.93515 * 25 24 23 48 48 H 1.09003 * 110.36512 * 217.27089 * 25 24 23 49 49 H 1.09005 * 109.47113 * 175.27828 * 29 27 26 50 50 H 1.08999 * 109.47312 * 295.28184 * 29 27 26 51 51 H 1.08998 * 109.47157 * 55.27799 * 29 27 26 52 52 H 1.09003 * 109.49569 * 58.78858 * 30 17 2 53 53 H 1.08996 * 109.49694 * 298.69065 * 30 17 2 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.8933 1.0276 0.0000 4 6 2.0324 -0.7104 -1.2304 5 8 2.2251 -1.9076 -1.2032 6 7 2.2675 -0.0155 -2.3610 7 6 1.9271 1.4077 -2.4302 8 6 2.8634 -0.6834 -3.5208 9 6 1.7706 -1.2457 -4.3930 10 1 0.7357 -1.1179 -4.1117 11 6 2.0920 -1.9223 -5.5505 12 7 1.1384 -2.4126 -6.3119 13 14 3.8819 -2.1426 -6.0373 14 1 4.4030 -3.4033 -5.4504 15 1 3.9885 -2.2043 -7.5172 16 1 4.6774 -0.9953 -5.5312 17 7 2.0197 -0.6926 1.1995 18 6 3.4887 -0.7476 1.2132 19 6 3.9561 -1.5363 2.4393 20 1 3.5921 -2.5616 2.3732 21 6 5.4852 -1.5380 2.4937 22 8 3.4363 -0.9196 3.6219 23 6 2.0171 -0.8173 3.6402 24 6 1.3804 -2.2186 3.6886 25 6 0.5863 -2.2196 5.0152 26 8 1.2335 -1.2205 5.8339 27 6 1.5450 -0.1320 4.9353 28 1 0.6573 0.4708 4.7441 29 6 2.6626 0.7362 5.5169 30 6 1.5124 -0.0511 2.4227 31 1 -0.3633 -1.0226 0.1022 32 1 -0.3633 0.5997 0.8344 33 1 -0.3634 0.4227 -0.9367 34 1 2.8050 2.0037 -2.1810 35 1 1.5956 1.6534 -3.4390 36 1 1.1275 1.6249 -1.7220 37 1 3.4505 0.0356 -4.0921 38 1 3.5100 -1.4926 -3.1815 39 1 1.3645 -2.8889 -7.1261 40 1 3.8901 0.2648 1.2575 41 1 3.8429 -1.2398 0.3075 42 1 5.8170 -2.0954 3.3696 43 1 5.8802 -2.0078 1.5929 44 1 5.8476 -0.5119 2.5571 45 1 0.7100 -2.3637 2.8414 46 1 2.1514 -2.9889 3.7036 47 1 -0.4535 -1.9471 4.8351 48 1 0.6461 -3.1979 5.4923 49 1 2.9430 1.5008 4.7924 50 1 2.3128 1.2137 6.4321 51 1 3.5281 0.1126 5.7404 52 1 0.4224 -0.0603 2.4131 53 1 1.8676 0.9784 2.4677 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 ZINC000369303960.mol2 53 1SCF WAS SPECIFIED, SO GEOMETRY OPTIMIZATION WAS NOT USED SCF FIELD WAS ACHIEVED Wall clock time and date at job start Wed Apr 1 2020 00:25:14 Heat of formation + Delta-G solvation = -102.803501 kcal Electronic energy + Delta-G solvation = -32614.598855 eV Core-core repulsion = 28487.178702 eV Total energy + Delta-G solvation = -4127.420153 eV No. of doubly occupied orbitals = 66 Molecular weight (most abundant/longest-lived isotopes) = 340.220 amu Computer time = 1.03 seconds Orbital eigenvalues (eV) -41.10059 -39.51046 -38.06927 -35.98566 -34.90815 -33.41092 -32.63240 -31.73815 -31.52658 -30.05970 -28.48560 -27.66708 -26.06450 -25.35450 -24.63852 -23.14909 -22.68439 -21.63528 -21.10765 -19.41826 -18.41455 -17.08487 -16.97993 -16.77565 -16.66524 -15.96221 -15.84305 -15.56278 -15.20759 -14.94941 -14.48805 -14.44305 -14.15055 -13.78335 -13.71057 -13.18958 -13.03060 -12.84596 -12.65435 -12.56620 -12.45561 -12.42126 -12.31804 -12.27434 -12.08279 -11.79895 -11.68622 -11.59810 -11.49646 -11.37799 -11.30767 -11.10329 -11.00107 -10.90750 -10.56448 -10.41775 -10.17870 -9.92210 -9.67779 -9.56112 -8.96140 -8.61487 -8.36704 -7.74552 -6.10937 -4.09678 3.18343 3.26361 3.31928 3.33800 3.35952 3.43425 3.68853 3.94946 4.36502 4.36953 4.45095 4.61399 4.74618 4.82052 4.89699 4.98844 5.02884 5.08406 5.11203 5.14143 5.20358 5.22056 5.24125 5.25294 5.34672 5.44594 5.48544 5.51314 5.60239 5.64325 5.66235 5.71214 5.73919 5.76112 5.78077 5.87473 5.92017 5.95745 6.06606 6.08664 6.11621 6.25745 6.60499 6.62002 6.79886 6.84584 6.92001 7.14454 7.29062 7.57934 8.19343 8.44265 8.90460 8.95879 9.49144 11.02530 Molecular weight = 340.22amu Principal moments of inertia in cm(-1) A = 0.018794 B = 0.002985 C = 0.002871 Principal moments of inertia in units of 10**(-40)*gram-cm**2 A = 1489.494720 B = 9378.430028 C = 9750.498953 Net atomic charges, atomic populations, and dipole contributions calculated with CM2 Atom NO. Type Charge No. of electrons 1 C -0.141 4.141 2 C 0.066 3.934 3 H 0.057 0.943 4 C 0.511 3.489 5 O -0.517 6.517 6 N -0.600 5.600 7 C 0.081 3.919 8 C 0.218 3.782 9 C -0.527 4.527 10 H 0.074 0.926 11 C 0.072 3.928 12 N -0.897 5.897 13 Si 0.876 3.124 14 H -0.273 1.273 15 H -0.285 1.285 16 H -0.267 1.267 17 N -0.509 5.509 18 C 0.022 3.978 19 C 0.090 3.910 20 H 0.072 0.928 21 C -0.142 4.142 22 O -0.360 6.360 23 C 0.098 3.902 24 C -0.171 4.171 25 C 0.045 3.955 26 O -0.369 6.369 27 C 0.087 3.913 28 H 0.079 0.921 29 C -0.135 4.135 30 C 0.034 3.966 31 H 0.074 0.926 32 H 0.061 0.939 33 H 0.070 0.930 34 H 0.044 0.956 35 H 0.076 0.924 36 H 0.061 0.939 37 H 0.029 0.971 38 H 0.034 0.966 39 H 0.265 0.735 40 H 0.038 0.962 41 H 0.096 0.904 42 H 0.065 0.935 43 H 0.074 0.926 44 H 0.061 0.939 45 H 0.095 0.905 46 H 0.083 0.917 47 H 0.060 0.940 48 H 0.092 0.908 49 H 0.068 0.932 50 H 0.058 0.942 51 H 0.073 0.927 52 H 0.091 0.909 53 H 0.043 0.957 For the charges calculated by CM2: Dipole moment (debyes) X Y Z Total from point charges 2.130 8.494 21.815 23.507 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.199 4.199 2 C -0.046 4.046 3 H 0.075 0.925 4 C 0.296 3.704 5 O -0.392 6.392 6 N -0.333 5.333 7 C -0.063 4.063 8 C 0.098 3.902 9 C -0.566 4.566 10 H 0.092 0.908 11 C -0.130 4.130 12 N -0.536 5.536 13 Si 0.703 3.297 14 H -0.199 1.199 15 H -0.210 1.210 16 H -0.191 1.191 17 N -0.234 5.234 18 C -0.106 4.106 19 C 0.033 3.967 20 H 0.090 0.910 21 C -0.200 4.200 22 O -0.279 6.279 23 C 0.058 3.942 24 C -0.209 4.209 25 C -0.031 4.031 26 O -0.288 6.288 27 C 0.029 3.971 28 H 0.097 0.903 29 C -0.193 4.193 30 C -0.094 4.094 31 H 0.093 0.907 32 H 0.080 0.920 33 H 0.089 0.911 34 H 0.062 0.938 35 H 0.095 0.905 36 H 0.079 0.921 37 H 0.047 0.953 38 H 0.052 0.948 39 H 0.074 0.926 40 H 0.056 0.944 41 H 0.115 0.885 42 H 0.084 0.916 43 H 0.093 0.907 44 H 0.080 0.920 45 H 0.113 0.887 46 H 0.102 0.898 47 H 0.079 0.921 48 H 0.110 0.890 49 H 0.087 0.913 50 H 0.077 0.923 51 H 0.092 0.908 52 H 0.110 0.890 53 H 0.061 0.939 Dipole moment (debyes) X Y Z Total from point charges 2.201 7.319 21.983 23.274 hybrid contribution 0.822 0.617 -1.077 1.489 sum 3.023 7.936 20.906 22.565 Atomic orbital electron populations 1.21792 0.92054 1.03767 1.02270 1.21701 0.97432 0.95698 0.89736 0.92533 1.20041 0.78879 0.87671 0.83782 1.90687 1.50691 1.13980 1.83811 1.47868 1.61848 1.09670 1.13885 1.21706 0.99246 0.82025 1.03347 1.19639 0.94579 0.92961 0.83059 1.21306 0.94191 1.35439 1.05640 0.90765 1.24860 0.97237 0.94826 0.96101 1.69892 1.32175 1.37104 1.14380 0.86859 0.85675 0.79361 0.77822 1.19898 1.21042 1.19138 1.61993 1.08974 1.52229 1.00243 1.22877 0.87494 0.99615 1.00593 1.22185 0.94396 0.96653 0.83496 0.91040 1.21726 0.92638 1.02815 1.02798 1.87912 1.18672 1.81101 1.40248 1.22324 0.82473 0.94921 0.94531 1.23447 1.01201 0.97483 0.98778 1.23428 0.95880 0.93011 0.90831 1.89183 1.78926 1.28695 1.31978 1.22864 0.97704 0.89102 0.87443 0.90298 1.21569 0.98106 0.98892 1.00684 1.22691 1.01171 0.96041 0.89531 0.90656 0.92042 0.91129 0.93791 0.90484 0.92054 0.95274 0.94804 0.92552 0.94390 0.88550 0.91597 0.90700 0.92001 0.88656 0.89849 0.92136 0.89012 0.91346 0.92320 0.90794 0.89030 0.93884 Number of geometries 1 Number of calculations of the screened coulomb radii 1 The total number of SCF iterations 13. 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.80 8.44 37.16 0.31 -1.49 16 2 C 0.07 0.95 1.61 -68.86 -0.11 0.84 16 3 H 0.06 0.72 6.77 -51.93 -0.35 0.37 16 4 C 0.51 10.08 5.68 -10.99 -0.06 10.02 16 5 O -0.52 -12.44 15.28 5.55 0.08 -12.35 16 6 N -0.60 -12.03 2.86 -182.96 -0.52 -12.55 16 7 C 0.08 1.20 9.87 59.84 0.59 1.80 16 8 C 0.22 5.25 4.21 -5.19 -0.02 5.23 16 9 C -0.53 -15.02 9.94 -34.44 -0.34 -15.36 16 10 H 0.07 2.33 8.06 -52.49 -0.42 1.90 16 11 C 0.07 2.09 5.47 -17.82 -0.10 1.99 16 12 N -0.90 -27.56 13.96 55.43 0.77 -26.78 16 13 Si 0.88 21.39 27.47 -169.99 -4.67 16.72 16 14 H -0.27 -6.72 7.11 56.52 0.40 -6.32 16 15 H -0.28 -6.83 7.11 56.52 0.40 -6.43 16 16 H -0.27 -6.10 6.54 56.52 0.37 -5.73 16 17 N -0.51 -7.01 4.09 -228.02 -0.93 -7.95 16 18 C 0.02 0.31 4.57 -3.71 -0.02 0.29 16 19 C 0.09 1.12 3.13 -26.68 -0.08 1.04 16 20 H 0.07 0.96 6.80 -51.93 -0.35 0.61 16 21 C -0.14 -1.47 9.57 37.16 0.36 -1.12 16 22 O -0.36 -4.41 8.30 -35.23 -0.29 -4.70 16 23 C 0.10 1.08 0.90 -89.89 -0.08 1.00 16 24 C -0.17 -1.87 5.77 -25.37 -0.15 -2.02 16 25 C 0.05 0.44 7.65 38.00 0.29 0.73 16 26 O -0.37 -4.28 11.35 -35.23 -0.40 -4.68 16 27 C 0.09 0.87 3.46 -26.24 -0.09 0.78 16 28 H 0.08 0.67 7.84 -51.93 -0.41 0.26 16 29 C -0.14 -1.26 8.96 37.16 0.33 -0.93 16 30 C 0.03 0.37 3.60 -3.97 -0.01 0.36 16 31 H 0.07 1.08 8.14 -51.93 -0.42 0.66 16 32 H 0.06 0.59 6.68 -51.93 -0.35 0.25 16 33 H 0.07 0.92 7.22 -51.93 -0.37 0.55 16 34 H 0.04 0.55 8.14 -51.93 -0.42 0.13 16 35 H 0.08 1.20 8.10 -51.93 -0.42 0.78 16 36 H 0.06 0.76 5.39 -51.93 -0.28 0.48 16 37 H 0.03 0.67 7.13 -51.93 -0.37 0.30 16 38 H 0.03 0.87 6.93 -51.93 -0.36 0.51 16 39 H 0.26 7.59 8.31 -40.82 -0.34 7.25 16 40 H 0.04 0.50 8.14 -51.93 -0.42 0.08 16 41 H 0.10 1.59 6.29 -51.93 -0.33 1.27 16 42 H 0.06 0.63 8.14 -51.93 -0.42 0.21 16 43 H 0.07 0.75 8.14 -51.93 -0.42 0.33 16 44 H 0.06 0.61 8.14 -51.92 -0.42 0.19 16 45 H 0.09 1.09 7.65 -51.93 -0.40 0.70 16 46 H 0.08 0.98 6.86 -51.93 -0.36 0.62 16 47 H 0.06 0.50 8.14 -51.93 -0.42 0.07 16 48 H 0.09 0.78 8.14 -51.93 -0.42 0.36 16 49 H 0.07 0.58 8.14 -51.93 -0.42 0.16 16 50 H 0.06 0.48 8.14 -51.93 -0.42 0.06 16 51 H 0.07 0.79 7.48 -51.93 -0.39 0.40 16 52 H 0.09 0.89 6.34 -51.93 -0.33 0.56 16 53 H 0.04 0.45 8.14 -51.93 -0.42 0.02 16 LS Contribution 400.27 15.07 6.03 6.03 Total: -1.00 -34.11 400.27 -8.41 -42.52 By element: Atomic # 1 Polarization: 9.89 SS G_CDS: -9.30 Total: 0.59 kcal Atomic # 6 Polarization: 2.34 SS G_CDS: 0.82 Total: 3.15 kcal Atomic # 7 Polarization: -46.59 SS G_CDS: -0.68 Total: -47.28 kcal Atomic # 8 Polarization: -21.13 SS G_CDS: -0.61 Total: -21.74 kcal Atomic # 14 Polarization: 21.39 SS G_CDS: -4.67 Total: 16.72 kcal Total LS contribution 6.03 Total: 6.03 kcal Total: -34.11 -8.41 -42.52 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. ZINC000369303960.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.287 kcal (2) G-P(sol) polarization free energy of solvation -34.106 kcal (3) G-ENP(sol) elect.-nuc.-pol. free energy of system -94.392 kcal (4) G-CDS(sol) cavity-dispersion-solvent structure free energy -8.411 kcal (5) G-P-CDS(sol) = G-P(sol) + G-CDS(sol) = (2) + (4) -42.517 kcal (6) G-S(sol) free energy of system = (1) + (5) -102.804 kcal 1SCF WAS SPECIFIED, SO GEOMETRY OPTIMIZATION WAS NOT USED Total computer time = 1.03 seconds