Wall clock time and date at job start Sat Apr 11 2020 06:12:54 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.42899 * 1 3 3 C 1.42897 * 113.99986 * 2 1 4 4 H 1.09002 * 110.01573 * 331.39034 * 3 2 1 5 5 C 1.53874 * 110.07329 * 92.80599 * 3 2 1 6 6 C 1.54467 * 106.68370 * 119.46766 * 5 3 2 7 7 C 1.54886 * 104.10130 * 23.34032 * 6 5 3 8 8 H 1.09004 * 110.75834 * 80.71231 * 7 6 5 9 9 C 1.50707 * 110.75582 * 203.89399 * 7 6 5 10 10 O 1.21281 * 119.99951 * 15.38334 * 9 7 6 11 11 N 1.34778 * 119.99760 * 195.37928 * 9 7 6 12 12 C 1.46931 * 120.62496 * 185.44627 * 11 9 7 13 13 C 1.53196 * 108.77470 * 233.62900 * 12 11 9 14 14 C 1.53035 * 109.31242 * 305.37020 * 13 12 11 15 15 N 1.46498 * 109.45913 * 181.38262 * 14 13 12 16 16 C 1.36938 * 120.00316 * 299.97690 * 15 14 13 17 17 H 1.07996 * 119.99786 * 41.88992 * 16 15 14 18 18 C 1.38814 * 119.99764 * 221.89838 * 16 15 14 19 19 N 1.31625 * 119.99518 * 201.89382 * 18 16 15 20 20 Si 1.86798 * 120.00316 * 21.89341 * 18 16 15 21 21 H 1.48496 * 109.47124 * 264.63685 * 20 18 16 22 22 H 1.48503 * 109.46732 * 24.63896 * 20 18 16 23 23 H 1.48503 * 109.47101 * 144.63518 * 20 18 16 24 24 C 1.46506 * 120.00175 * 119.97108 * 15 14 13 25 25 C 1.53004 * 109.47213 * 275.39388 * 24 15 14 26 26 O 1.42902 * 109.47280 * 64.99230 * 25 24 15 27 27 C 1.53042 * 109.53085 * 61.37382 * 14 13 12 28 28 C 1.46928 * 120.63317 * 5.47181 * 11 9 7 29 29 C 1.54269 * 110.03425 * 209.99746 * 3 2 1 30 30 H 1.08998 * 109.47202 * 300.00033 * 1 2 3 31 31 H 1.08997 * 109.46795 * 60.00607 * 1 2 3 32 32 H 1.09000 * 109.46702 * 180.02562 * 1 2 3 33 33 H 1.09003 * 110.03351 * 238.58908 * 5 3 2 34 34 H 1.09001 * 110.03429 * 0.02562 * 5 3 2 35 35 H 1.09002 * 110.40598 * 264.75646 * 6 5 3 36 36 H 1.09003 * 110.66337 * 142.05595 * 6 5 3 37 37 H 1.09003 * 109.57987 * 113.84470 * 12 11 9 38 38 H 1.08991 * 109.70579 * 353.48986 * 12 11 9 39 39 H 1.08994 * 109.49688 * 185.41940 * 13 12 11 40 40 H 1.09001 * 109.51989 * 65.33509 * 13 12 11 41 41 H 1.08996 * 109.46544 * 301.35662 * 14 13 12 42 42 H 0.96997 * 119.99581 * 179.97438 * 19 18 16 43 43 H 1.09002 * 109.46573 * 35.38679 * 24 15 14 44 44 H 1.08993 * 109.47431 * 155.38632 * 24 15 14 45 45 H 1.08994 * 109.47308 * 184.99950 * 25 24 15 46 46 H 1.09000 * 109.46631 * 304.99926 * 25 24 15 47 47 H 0.96693 * 113.99940 * 179.97438 * 26 25 24 48 48 H 1.09000 * 109.49553 * 58.57340 * 27 14 13 49 49 H 1.08994 * 109.50238 * 178.67366 * 27 14 13 50 50 H 1.08998 * 109.58460 * 246.15574 * 28 11 9 51 51 H 1.08996 * 109.58382 * 6.44074 * 28 11 9 52 52 H 1.08999 * 110.48309 * 98.37174 * 29 3 2 53 53 H 1.09001 * 110.48101 * 335.73506 * 29 3 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 8 1.4290 0.0000 0.0000 3 6 2.0102 1.3054 0.0000 4 1 1.3406 2.0120 -0.4904 5 6 2.2896 1.7591 1.4435 6 6 3.8125 1.9961 1.5470 7 6 4.2535 2.2599 0.0858 8 1 4.0461 3.2912 -0.2001 9 6 5.7150 1.9410 -0.0968 10 8 6.3014 1.2846 0.7376 11 7 6.3704 2.3851 -1.1876 12 6 7.7613 1.9841 -1.4400 13 6 8.6086 3.2442 -1.6427 14 6 7.9893 4.0976 -2.7518 15 7 8.8154 5.2874 -2.9710 16 6 9.0236 6.1743 -1.9487 17 1 9.2206 5.8115 -0.9508 18 6 8.9813 7.5398 -2.1951 19 7 9.5825 8.3754 -1.3748 20 14 8.0714 8.1913 -3.6907 21 1 9.0305 8.3863 -4.8075 22 1 7.0280 7.2165 -4.0986 23 1 7.4322 9.4890 -3.3547 24 6 9.4190 5.5283 -4.2840 25 6 10.7394 4.7620 -4.3868 26 8 10.4797 3.3575 -4.3419 27 6 6.5785 4.5253 -2.3408 28 6 5.7080 3.2798 -2.1465 29 6 3.3718 1.2777 -0.7246 30 1 -0.3633 0.5138 0.8900 31 1 -0.3633 0.5137 -0.8900 32 1 -0.3633 -1.0277 0.0005 33 1 1.7557 2.6859 1.6537 34 1 1.9824 0.9834 2.1449 35 1 4.3092 1.1095 1.9411 36 1 4.0253 2.8632 2.1722 37 1 7.8061 1.3654 -2.3363 38 1 8.1418 1.4210 -0.5879 39 1 9.6220 2.9595 -1.9257 40 1 8.6370 3.8170 -0.7158 41 1 7.9387 3.5160 -3.6722 42 1 9.5532 9.3295 -1.5471 43 1 8.7391 5.1868 -5.0646 44 1 9.6073 6.5948 -4.4072 45 1 11.2316 5.0107 -5.3269 46 1 11.3857 5.0388 -3.5539 47 1 11.2735 2.8087 -4.4027 48 1 6.6258 5.0857 -1.4071 49 1 6.1468 5.1525 -3.1207 50 1 5.5855 2.7668 -3.1003 51 1 4.7313 3.5707 -1.7600 52 1 3.2636 1.6218 -1.7532 53 1 3.7967 0.2742 -0.6998 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 ZINC001098413860.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:12:54 Heat of formation + Delta-G solvation = -113.074243 kcal Electronic energy + Delta-G solvation = -31762.549182 eV Core-core repulsion = 27634.683656 eV Total energy + Delta-G solvation = -4127.865526 eV No. of doubly occupied orbitals = 66 Molecular weight (most abundant/longest-lived isotopes) = 340.220 amu Computer time = 0.88 seconds Orbital eigenvalues (eV) -40.35975 -38.84577 -37.69216 -36.39762 -35.73444 -34.78819 -32.10583 -31.12841 -30.67633 -29.04743 -28.70724 -27.79967 -26.37872 -25.00037 -23.36128 -22.84545 -21.91125 -20.81123 -20.59417 -19.62862 -18.58086 -17.57466 -16.92072 -16.62569 -16.28694 -16.22373 -15.52690 -15.16226 -14.96340 -14.96047 -14.86327 -14.26887 -13.92674 -13.75670 -13.49216 -13.31264 -13.21995 -13.01629 -12.69746 -12.46967 -12.28629 -12.10322 -12.05275 -11.98049 -11.69009 -11.59679 -11.50175 -11.23670 -11.13972 -10.97627 -10.76346 -10.66535 -10.56617 -10.44461 -10.19044 -10.10824 -9.83893 -9.73423 -9.56851 -9.40015 -9.27561 -8.54955 -8.30335 -6.73276 -6.03225 -3.47756 2.78974 3.21881 3.25066 3.41718 3.46267 3.56519 4.27927 4.35882 4.47457 4.71520 4.78014 4.85254 4.94903 5.00125 5.03712 5.07031 5.19606 5.25792 5.26864 5.28701 5.34756 5.40620 5.50623 5.55617 5.56775 5.64275 5.75281 5.82336 5.83101 5.89128 5.91401 5.99935 6.01742 6.05569 6.09605 6.20166 6.22993 6.26742 6.36042 6.56390 6.58049 6.65921 6.74456 6.89724 6.97939 7.14599 7.16612 7.68497 8.05472 8.15302 8.41136 8.43402 8.99936 9.69048 10.14908 11.12411 Molecular weight = 340.22amu Principal moments of inertia in cm(-1) A = 0.015204 B = 0.002924 C = 0.002745 Principal moments of inertia in units of 10**(-40)*gram-cm**2 A = 1841.204233 B = 9573.715665 C =10196.398135 Net atomic charges, atomic populations, and dipole contributions calculated with CM2 Atom NO. Type Charge No. of electrons 1 C 0.030 3.970 2 O -0.383 6.383 3 C 0.085 3.915 4 H 0.067 0.933 5 C -0.147 4.147 6 C -0.104 4.104 7 C -0.116 4.116 8 H 0.102 0.898 9 C 0.524 3.476 10 O -0.547 6.547 11 N -0.603 5.603 12 C 0.114 3.886 13 C -0.133 4.133 14 C 0.191 3.809 15 N -0.595 5.595 16 C -0.262 4.262 17 H 0.076 0.924 18 C 0.037 3.963 19 N -0.899 5.899 20 Si 0.871 3.129 21 H -0.288 1.288 22 H -0.266 1.266 23 H -0.301 1.301 24 C 0.140 3.860 25 C 0.057 3.943 26 O -0.566 6.566 27 C -0.145 4.145 28 C 0.108 3.892 29 C -0.105 4.105 30 H 0.040 0.960 31 H 0.041 0.959 32 H 0.084 0.916 33 H 0.068 0.932 34 H 0.076 0.924 35 H 0.085 0.915 36 H 0.071 0.929 37 H 0.061 0.939 38 H 0.083 0.917 39 H 0.074 0.926 40 H 0.065 0.935 41 H 0.056 0.944 42 H 0.262 0.738 43 H 0.040 0.960 44 H 0.069 0.931 45 H 0.025 0.975 46 H 0.050 0.950 47 H 0.362 0.638 48 H 0.076 0.924 49 H 0.066 0.934 50 H 0.063 0.937 51 H 0.078 0.922 52 H 0.078 0.922 53 H 0.084 0.916 For the charges calculated by CM2: Dipole moment (debyes) X Y Z Total from point charges -13.537 -12.970 -2.227 18.880 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.066 4.066 2 O -0.301 6.301 3 C 0.028 3.972 4 H 0.085 0.915 5 C -0.184 4.184 6 C -0.142 4.142 7 C -0.138 4.138 8 H 0.120 0.880 9 C 0.312 3.688 10 O -0.425 6.425 11 N -0.336 5.336 12 C -0.008 4.008 13 C -0.172 4.172 14 C 0.084 3.916 15 N -0.315 5.315 16 C -0.391 4.391 17 H 0.094 0.906 18 C -0.163 4.163 19 N -0.540 5.540 20 Si 0.701 3.299 21 H -0.216 1.216 22 H -0.191 1.191 23 H -0.229 1.229 24 C 0.014 3.986 25 C -0.023 4.023 26 O -0.369 6.369 27 C -0.185 4.185 28 C -0.015 4.015 29 C -0.143 4.143 30 H 0.059 0.941 31 H 0.059 0.941 32 H 0.102 0.898 33 H 0.087 0.913 34 H 0.095 0.905 35 H 0.103 0.897 36 H 0.090 0.910 37 H 0.079 0.921 38 H 0.102 0.898 39 H 0.092 0.908 40 H 0.084 0.916 41 H 0.074 0.926 42 H 0.072 0.928 43 H 0.058 0.942 44 H 0.087 0.913 45 H 0.043 0.957 46 H 0.068 0.932 47 H 0.206 0.794 48 H 0.095 0.905 49 H 0.085 0.915 50 H 0.082 0.918 51 H 0.096 0.904 52 H 0.096 0.904 53 H 0.103 0.897 Dipole moment (debyes) X Y Z Total from point charges -13.698 -13.066 -1.691 19.006 hybrid contribution -0.010 0.919 -0.976 1.341 sum -13.709 -12.147 -2.668 18.509 Atomic orbital electron populations 1.22776 0.81354 1.02630 0.99890 1.87968 1.19404 1.27807 1.94970 1.22822 0.94562 0.84400 0.95460 0.91528 1.22730 0.97102 1.01558 0.97056 1.22244 0.94245 1.03273 0.94403 1.22270 0.92102 1.01156 0.98225 0.88017 1.20296 0.88539 0.79203 0.80787 1.90652 1.68960 1.40704 1.42212 1.48142 1.13984 1.47784 1.23642 1.21670 0.80360 0.97671 1.01128 1.22233 0.98950 0.94933 1.01131 1.20207 0.88101 0.87908 0.95408 1.44069 1.54473 1.24287 1.08721 1.18686 1.43883 0.85647 0.90873 0.90608 1.25610 0.97753 0.96238 0.96736 1.69883 1.42243 1.02949 1.38947 0.86978 0.81562 0.77977 0.83404 1.21597 1.19073 1.22881 1.21262 0.91709 0.97389 0.88271 1.22993 1.00306 0.80958 0.98002 1.86524 1.35550 1.17635 1.97238 1.22411 0.97168 0.95654 1.03284 1.21882 0.96248 0.91277 0.92141 1.22332 0.93159 1.00430 0.98370 0.94103 0.94079 0.89782 0.91311 0.90513 0.89664 0.91020 0.92113 0.89837 0.90782 0.91590 0.92558 0.92762 0.94162 0.91295 0.95655 0.93167 0.79381 0.90541 0.91499 0.91836 0.90422 0.90375 0.89724 Number of geometries 1 Number of calculations of the screened coulomb radii 1 The total number of SCF iterations 17. 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.03 0.22 10.85 101.05 1.10 1.31 16 2 O -0.38 -4.04 10.26 -35.23 -0.36 -4.41 16 3 C 0.09 0.81 3.38 -25.61 -0.09 0.73 16 4 H 0.07 0.53 7.80 -51.93 -0.41 0.12 16 5 C -0.15 -1.38 6.62 -25.51 -0.17 -1.55 16 6 C -0.10 -1.21 6.07 -24.97 -0.15 -1.36 16 7 C -0.12 -1.39 2.31 -89.78 -0.21 -1.60 16 8 H 0.10 1.09 6.45 -51.93 -0.33 0.76 16 9 C 0.52 8.19 6.94 -10.98 -0.08 8.11 16 10 O -0.55 -10.38 15.14 5.56 0.08 -10.30 16 11 N -0.60 -9.17 2.97 -172.76 -0.51 -9.68 16 12 C 0.11 1.82 6.44 -3.71 -0.02 1.80 16 13 C -0.13 -2.43 5.30 -26.61 -0.14 -2.57 16 14 C 0.19 3.61 1.96 -67.89 -0.13 3.47 16 15 N -0.59 -13.18 1.92 -165.62 -0.32 -13.50 16 16 C -0.26 -6.88 8.61 -15.06 -0.13 -7.01 16 17 H 0.08 2.22 7.10 -52.49 -0.37 1.85 16 18 C 0.04 1.02 5.24 -17.43 -0.09 0.93 16 19 N -0.90 -26.59 13.97 55.49 0.78 -25.81 16 20 Si 0.87 21.05 25.26 -169.99 -4.29 16.75 16 21 H -0.29 -6.80 6.73 56.52 0.38 -6.42 16 22 H -0.27 -6.20 7.11 56.52 0.40 -5.80 16 23 H -0.30 -7.49 7.11 56.52 0.40 -7.09 16 24 C 0.14 2.88 4.70 -4.04 -0.02 2.87 16 25 C 0.06 1.02 7.39 37.16 0.27 1.30 16 26 O -0.57 -9.10 12.20 -35.23 -0.43 -9.53 16 27 C -0.15 -2.54 5.31 -26.61 -0.14 -2.68 16 28 C 0.11 1.46 6.32 -3.71 -0.02 1.44 16 29 C -0.10 -1.10 6.00 -25.07 -0.15 -1.25 16 30 H 0.04 0.26 8.14 -51.93 -0.42 -0.17 16 31 H 0.04 0.25 8.10 -51.93 -0.42 -0.17 16 32 H 0.08 0.50 8.14 -51.93 -0.42 0.08 16 33 H 0.07 0.55 8.14 -51.93 -0.42 0.13 16 34 H 0.08 0.73 7.96 -51.93 -0.41 0.32 16 35 H 0.08 1.16 7.26 -51.93 -0.38 0.78 16 36 H 0.07 0.81 8.14 -51.93 -0.42 0.39 16 37 H 0.06 0.88 8.14 -51.93 -0.42 0.46 16 38 H 0.08 1.41 7.01 -51.93 -0.36 1.05 16 39 H 0.07 1.38 7.61 -51.93 -0.40 0.99 16 40 H 0.07 1.40 7.10 -51.93 -0.37 1.03 16 41 H 0.06 1.04 7.60 -51.93 -0.39 0.65 16 42 H 0.26 7.39 8.31 -40.82 -0.34 7.05 16 43 H 0.04 0.78 8.04 -51.93 -0.42 0.36 16 44 H 0.07 1.54 4.12 -51.93 -0.21 1.33 16 45 H 0.03 0.39 8.14 -51.93 -0.42 -0.03 16 46 H 0.05 0.98 8.14 -51.93 -0.42 0.56 16 47 H 0.36 4.14 9.12 45.56 0.42 4.56 16 48 H 0.08 1.54 7.71 -51.93 -0.40 1.13 16 49 H 0.07 1.20 8.14 -51.93 -0.42 0.77 16 50 H 0.06 0.77 8.14 -51.93 -0.42 0.35 16 51 H 0.08 0.88 5.52 -51.93 -0.29 0.59 16 52 H 0.08 0.71 8.14 -51.93 -0.42 0.28 16 53 H 0.08 1.03 7.98 -51.93 -0.41 0.62 16 LS Contribution 402.26 15.07 6.06 6.06 Total: -1.00 -32.24 402.26 -7.71 -39.95 By element: Atomic # 1 Polarization: 15.07 SS G_CDS: -8.55 Total: 6.52 kcal Atomic # 6 Polarization: 4.10 SS G_CDS: -0.17 Total: 3.92 kcal Atomic # 7 Polarization: -48.93 SS G_CDS: -0.06 Total: -48.99 kcal Atomic # 8 Polarization: -23.53 SS G_CDS: -0.71 Total: -24.23 kcal Atomic # 14 Polarization: 21.05 SS G_CDS: -4.29 Total: 16.75 kcal Total LS contribution 6.06 Total: 6.06 kcal Total: -32.24 -7.71 -39.95 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. ZINC001098413860.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 -73.120 kcal (2) G-P(sol) polarization free energy of solvation -32.241 kcal (3) G-ENP(sol) elect.-nuc.-pol. free energy of system -105.362 kcal (4) G-CDS(sol) cavity-dispersion-solvent structure free energy -7.713 kcal (5) G-P-CDS(sol) = G-P(sol) + G-CDS(sol) = (2) + (4) -39.954 kcal (6) G-S(sol) free energy of system = (1) + (5) -113.074 kcal 1SCF WAS SPECIFIED, SO GEOMETRY OPTIMIZATION WAS NOT USED Total computer time = 0.88 seconds