Wall clock time and date at job start Sat Apr 11 2020 06:13:25 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.42894 * 1 3 3 C 1.42898 * 113.99952 * 2 1 4 4 C 1.53006 * 109.47414 * 180.02562 * 3 2 1 5 5 C 1.50690 * 109.47411 * 180.02562 * 4 3 2 6 6 O 1.21286 * 119.99911 * 0.02562 * 5 4 3 7 7 N 1.34783 * 120.00261 * 179.97438 * 5 4 3 8 8 C 1.46923 * 120.62923 * 179.97438 * 7 5 4 9 9 C 1.53196 * 108.77351 * 233.63020 * 8 7 5 10 10 C 1.53038 * 109.31174 * 305.36629 * 9 8 7 11 11 C 1.53001 * 109.45995 * 181.38407 * 10 9 8 12 12 H 1.09002 * 109.48301 * 300.00196 * 11 10 9 13 13 C 1.53093 * 109.48321 * 179.97438 * 11 10 9 14 14 C 1.53178 * 109.16373 * 182.40028 * 13 11 10 15 15 H 1.08997 * 109.54056 * 297.02111 * 14 13 11 16 16 N 1.46498 * 109.54151 * 176.84819 * 14 13 11 17 17 C 1.36943 * 119.99727 * 155.17693 * 16 14 13 18 18 H 1.07998 * 120.00392 * 359.97438 * 17 16 14 19 19 C 1.38811 * 119.99586 * 179.97438 * 17 16 14 20 20 N 1.31618 * 120.00460 * 0.02562 * 19 17 16 21 21 Si 1.86798 * 119.99694 * 179.97438 * 19 17 16 22 22 H 1.48506 * 109.47297 * 179.97438 * 21 19 17 23 23 H 1.48503 * 109.47261 * 60.00129 * 21 19 17 24 24 H 1.48500 * 109.47454 * 300.00122 * 21 19 17 25 25 C 1.53178 * 109.11069 * 56.93509 * 14 13 11 26 26 C 1.53099 * 109.15641 * 303.15102 * 25 14 13 27 27 O 1.42902 * 109.48802 * 60.02490 * 11 10 9 28 28 C 1.53045 * 109.52498 * 61.37705 * 10 9 8 29 29 C 1.46920 * 120.62544 * 0.02562 * 7 5 4 30 30 H 1.09000 * 109.47115 * 179.97438 * 1 2 3 31 31 H 1.08997 * 109.47113 * 299.99310 * 1 2 3 32 32 H 1.09002 * 109.47411 * 59.99720 * 1 2 3 33 33 H 1.08999 * 109.47440 * 60.00377 * 3 2 1 34 34 H 1.08998 * 109.47446 * 300.00242 * 3 2 1 35 35 H 1.08999 * 109.46745 * 59.99663 * 4 3 2 36 36 H 1.08998 * 109.46776 * 300.00364 * 4 3 2 37 37 H 1.09005 * 109.58272 * 113.84747 * 8 7 5 38 38 H 1.08993 * 109.70997 * 353.49240 * 8 7 5 39 39 H 1.09002 * 109.49774 * 185.41761 * 9 8 7 40 40 H 1.08999 * 109.49646 * 65.32085 * 9 8 7 41 41 H 1.09001 * 109.46412 * 301.35833 * 10 9 8 42 42 H 1.08996 * 109.52059 * 62.49057 * 13 11 10 43 43 H 1.09002 * 109.57345 * 302.33753 * 13 11 10 44 44 H 0.96998 * 120.00095 * 335.17152 * 16 14 13 45 45 H 0.97001 * 119.99952 * 359.97438 * 20 19 17 46 46 H 1.08998 * 109.52100 * 62.96721 * 25 14 13 47 47 H 1.09000 * 109.52137 * 183.16122 * 25 14 13 48 48 H 1.08999 * 109.40898 * 297.49349 * 26 25 14 49 49 H 1.09004 * 109.49285 * 177.50477 * 26 25 14 50 50 H 1.08994 * 109.50141 * 58.57699 * 28 10 9 51 51 H 1.09003 * 109.49385 * 178.67610 * 28 10 9 52 52 H 1.08998 * 109.58717 * 246.15327 * 29 7 5 53 53 H 1.09001 * 109.59136 * 6.43759 * 29 7 5 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.4289 0.0000 0.0000 3 6 2.0101 1.3054 0.0000 4 6 3.5354 1.1847 -0.0006 5 6 4.1483 2.5613 0.0000 6 8 3.4354 3.5425 0.0013 7 7 5.4888 2.7023 0.0001 8 6 6.1010 4.0379 0.0002 9 6 7.0884 4.1277 -1.1677 10 6 8.0919 2.9762 -1.0717 11 6 9.1008 3.0825 -2.2170 12 1 8.5754 3.0280 -3.1705 13 6 10.1053 1.9312 -2.1211 14 6 11.1466 2.0783 -3.2347 15 1 10.6568 1.9971 -4.2050 16 7 12.1551 1.0236 -3.1051 17 6 12.8679 0.6193 -4.2023 18 1 12.6863 1.0777 -5.1632 19 6 13.8238 -0.3796 -4.0793 20 7 14.0448 -0.9386 -2.9084 21 14 14.7956 -0.9316 -5.5761 22 1 15.7464 -2.0021 -5.1818 23 1 15.5507 0.2214 -6.1289 24 1 13.8628 -1.4527 -6.6074 25 6 11.8200 3.4489 -3.1138 26 6 10.7490 4.5415 -3.1693 27 8 9.7956 4.3279 -2.1263 28 6 7.3495 1.6414 -1.1684 29 6 6.3653 1.5232 -0.0010 30 1 -0.3633 -1.0277 -0.0005 31 1 -0.3633 0.5137 0.8900 32 1 -0.3634 0.5139 -0.8900 33 1 1.6886 1.8464 0.8900 34 1 1.6885 1.8464 -0.8899 35 1 3.8567 0.6439 -0.8908 36 1 3.8574 0.6436 0.8891 37 1 6.6306 4.1970 0.9395 38 1 5.3274 4.7969 -0.1162 39 1 7.6200 5.0782 -1.1224 40 1 6.5447 4.0589 -2.1099 41 1 8.6176 3.0311 -0.1184 42 1 9.5826 0.9814 -2.2343 43 1 10.6023 1.9600 -1.1514 44 1 12.3180 0.6118 -2.2422 45 1 13.5399 -0.6522 -2.1313 46 1 12.3526 3.5105 -2.1648 47 1 12.5225 3.5834 -3.9363 48 1 10.2451 4.5061 -4.1352 49 1 11.2175 5.5170 -3.0392 50 1 6.8034 1.5962 -2.1106 51 1 8.0670 0.8220 -1.1236 52 1 6.9165 1.4769 0.9382 53 1 5.7645 0.6212 -0.1169 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 ZINC001186041285.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:13:25 Heat of formation + Delta-G solvation = -130.654673 kcal Electronic energy + Delta-G solvation = -30483.040763 eV Core-core repulsion = 26354.412892 eV Total energy + Delta-G solvation = -4128.627870 eV No. of doubly occupied orbitals = 66 Molecular weight (most abundant/longest-lived isotopes) = 340.220 amu Computer time = 0.97 seconds Orbital eigenvalues (eV) -40.67011 -39.46202 -38.22427 -37.41128 -35.28462 -33.97310 -32.65036 -32.32558 -31.05306 -30.02475 -29.00562 -27.13388 -25.94584 -24.40257 -23.62119 -23.56050 -22.65651 -21.91313 -20.52692 -19.86434 -19.11595 -17.81054 -17.35069 -16.72774 -16.53878 -16.08123 -15.99322 -15.76293 -15.49559 -15.16096 -14.76576 -14.48446 -14.17760 -14.09987 -14.02448 -13.68729 -13.59191 -13.56007 -13.42937 -12.94769 -12.79962 -12.60280 -12.28213 -12.15810 -11.84945 -11.59324 -11.54181 -11.44763 -11.17225 -10.97911 -10.89048 -10.78428 -10.67493 -10.55080 -10.25122 -10.08603 -10.07264 -10.05730 -9.91510 -9.42532 -9.08349 -8.92902 -8.91205 -7.08389 -6.14150 -3.14324 2.29666 2.54712 3.41757 3.48580 3.51804 3.69435 3.81584 4.25635 4.44959 4.49566 4.61315 4.64070 4.67530 4.79088 4.83500 4.96308 4.97160 5.01913 5.05839 5.16066 5.19065 5.25473 5.26552 5.28410 5.36818 5.39246 5.39557 5.42360 5.47492 5.52599 5.60846 5.63822 5.69878 5.70066 5.80869 5.85611 5.91856 5.96511 6.14870 6.18168 6.27464 6.29199 6.31578 6.50290 6.71777 6.71944 6.94621 7.06285 7.44806 7.76208 7.82965 8.71547 9.18254 10.04665 10.19887 11.34100 Molecular weight = 340.22amu Principal moments of inertia in cm(-1) A = 0.014386 B = 0.002137 C = 0.001924 Principal moments of inertia in units of 10**(-40)*gram-cm**2 A = 1945.929556 B =13098.296381 C =14553.187290 Net atomic charges, atomic populations, and dipole contributions calculated with CM2 Atom NO. Type Charge No. of electrons 1 C 0.034 3.966 2 O -0.395 6.395 3 C 0.076 3.924 4 C -0.141 4.141 5 C 0.521 3.479 6 O -0.547 6.547 7 N -0.609 5.609 8 C 0.111 3.889 9 C -0.124 4.124 10 C -0.077 4.077 11 C 0.102 3.898 12 H 0.054 0.946 13 C -0.154 4.154 14 C 0.183 3.817 15 H 0.046 0.954 16 N -0.734 5.734 17 C -0.332 4.332 18 H 0.063 0.937 19 C -0.045 4.045 20 N -0.943 5.943 21 Si 1.079 2.921 22 H -0.285 1.285 23 H -0.295 1.295 24 H -0.295 1.295 25 C -0.171 4.171 26 C 0.065 3.935 27 O -0.382 6.382 28 C -0.133 4.133 29 C 0.106 3.894 30 H 0.087 0.913 31 H 0.042 0.958 32 H 0.042 0.958 33 H 0.057 0.943 34 H 0.058 0.942 35 H 0.104 0.896 36 H 0.103 0.897 37 H 0.068 0.932 38 H 0.092 0.908 39 H 0.085 0.915 40 H 0.064 0.936 41 H 0.080 0.920 42 H 0.072 0.928 43 H 0.066 0.934 44 H 0.367 0.633 45 H 0.255 0.745 46 H 0.062 0.938 47 H 0.069 0.931 48 H 0.043 0.957 49 H 0.085 0.915 50 H 0.070 0.930 51 H 0.087 0.913 52 H 0.074 0.926 53 H 0.089 0.911 For the charges calculated by CM2: Dipole moment (debyes) X Y Z Total from point charges -22.624 7.051 7.003 24.710 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.063 4.063 2 O -0.314 6.314 3 C -0.001 4.001 4 C -0.181 4.181 5 C 0.309 3.691 6 O -0.425 6.425 7 N -0.343 5.343 8 C -0.011 4.011 9 C -0.163 4.163 10 C -0.095 4.095 11 C 0.044 3.956 12 H 0.072 0.928 13 C -0.192 4.192 14 C 0.073 3.927 15 H 0.064 0.936 16 N -0.376 5.376 17 C -0.463 4.463 18 H 0.082 0.918 19 C -0.243 4.243 20 N -0.588 5.588 21 Si 0.910 3.090 22 H -0.210 1.210 23 H -0.223 1.223 24 H -0.223 1.223 25 C -0.212 4.212 26 C -0.012 4.012 27 O -0.302 6.302 28 C -0.171 4.171 29 C -0.017 4.017 30 H 0.105 0.895 31 H 0.060 0.940 32 H 0.061 0.939 33 H 0.076 0.924 34 H 0.076 0.924 35 H 0.122 0.878 36 H 0.121 0.879 37 H 0.087 0.913 38 H 0.110 0.890 39 H 0.104 0.896 40 H 0.083 0.917 41 H 0.099 0.901 42 H 0.091 0.909 43 H 0.085 0.915 44 H 0.200 0.800 45 H 0.066 0.934 46 H 0.081 0.919 47 H 0.088 0.912 48 H 0.061 0.939 49 H 0.103 0.897 50 H 0.089 0.911 51 H 0.106 0.894 52 H 0.093 0.907 53 H 0.107 0.893 Dipole moment (debyes) X Y Z Total from point charges -21.751 6.551 6.388 23.598 hybrid contribution -0.985 0.971 0.455 1.456 sum -22.736 7.523 6.843 24.907 Atomic orbital electron populations 1.22738 0.80522 1.02825 1.00204 1.87807 1.20129 1.27942 1.95566 1.22161 0.91887 0.84264 1.01738 1.21517 0.95450 0.95361 1.05744 1.20381 0.81928 0.91052 0.75730 1.90700 1.60505 1.40890 1.50419 1.48146 1.05276 1.07905 1.72945 1.21660 0.96622 0.84031 0.98785 1.21903 0.96343 0.99781 0.98227 1.21103 0.94745 0.94495 0.99202 1.22007 0.91789 0.84329 0.97486 0.92828 1.22548 0.97679 0.98267 1.00745 1.20142 0.89022 0.87899 0.95671 0.93622 1.42420 1.43733 1.42328 1.09098 1.20028 1.17651 1.19306 0.89274 0.91849 1.25771 0.98806 0.96220 1.03473 1.70886 1.43907 1.44572 0.99450 0.83254 0.75994 0.76085 0.73678 1.21003 1.22289 1.22258 1.22908 0.98597 0.96988 1.02664 1.22573 0.89819 0.98903 0.89883 1.87813 1.54688 1.37003 1.50675 1.22019 0.98333 0.98097 0.98682 1.21834 0.92732 0.87690 0.99410 0.89458 0.93953 0.93919 0.92424 0.92393 0.87804 0.87917 0.91345 0.89006 0.89602 0.91731 0.90140 0.90914 0.91485 0.80011 0.93404 0.91930 0.91208 0.93856 0.89652 0.91111 0.89417 0.90742 0.89271 Number of geometries 1 Number of calculations of the screened coulomb radii 1 The total number of SCF iterations 15. 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.21 11.01 101.05 1.11 1.33 16 2 O -0.40 -3.64 10.70 -35.23 -0.38 -4.02 16 3 C 0.08 0.68 5.23 37.16 0.19 0.87 16 4 C -0.14 -1.12 4.42 -27.88 -0.12 -1.24 16 5 C 0.52 5.42 7.70 -10.99 -0.08 5.34 16 6 O -0.55 -7.51 15.40 5.55 0.09 -7.43 16 7 N -0.61 -5.66 2.97 -172.80 -0.51 -6.18 16 8 C 0.11 1.06 6.43 -3.71 -0.02 1.04 16 9 C -0.12 -1.24 5.23 -26.61 -0.14 -1.38 16 10 C -0.08 -0.81 1.95 -90.58 -0.18 -0.98 16 11 C 0.10 1.35 2.22 -26.68 -0.06 1.29 16 12 H 0.05 0.71 8.14 -51.93 -0.42 0.29 16 13 C -0.15 -2.41 4.47 -26.59 -0.12 -2.53 16 14 C 0.18 3.57 2.86 -67.75 -0.19 3.37 16 15 H 0.05 0.93 8.08 -51.93 -0.42 0.51 16 16 N -0.73 -18.02 5.13 -53.25 -0.27 -18.30 16 17 C -0.33 -9.19 9.93 -15.06 -0.15 -9.34 16 18 H 0.06 1.73 7.83 -52.49 -0.41 1.32 16 19 C -0.05 -1.30 5.50 -17.43 -0.10 -1.39 16 20 N -0.94 -28.86 13.70 55.49 0.76 -28.10 16 21 Si 1.08 25.68 29.92 -169.99 -5.09 20.59 16 22 H -0.28 -6.89 7.11 56.52 0.40 -6.49 16 23 H -0.30 -6.73 7.11 56.52 0.40 -6.33 16 24 H -0.30 -6.70 7.11 56.52 0.40 -6.30 16 25 C -0.17 -3.05 5.73 -26.59 -0.15 -3.20 16 26 C 0.07 0.94 6.84 37.21 0.25 1.20 16 27 O -0.38 -5.54 9.71 -35.23 -0.34 -5.88 16 28 C -0.13 -1.20 5.04 -26.63 -0.13 -1.33 16 29 C 0.11 0.80 6.34 -3.73 -0.02 0.77 16 30 H 0.09 0.43 8.14 -51.93 -0.42 0.00 16 31 H 0.04 0.23 8.14 -51.93 -0.42 -0.19 16 32 H 0.04 0.24 8.14 -51.93 -0.42 -0.18 16 33 H 0.06 0.52 8.10 -51.93 -0.42 0.10 16 34 H 0.06 0.53 8.10 -51.93 -0.42 0.11 16 35 H 0.10 0.73 7.69 -51.93 -0.40 0.33 16 36 H 0.10 0.67 7.94 -51.93 -0.41 0.26 16 37 H 0.07 0.60 8.14 -51.93 -0.42 0.18 16 38 H 0.09 0.98 7.01 -51.93 -0.36 0.61 16 39 H 0.09 0.89 7.64 -51.93 -0.40 0.49 16 40 H 0.06 0.68 8.14 -51.93 -0.42 0.26 16 41 H 0.08 0.88 8.14 -51.93 -0.42 0.45 16 42 H 0.07 1.11 6.54 -51.93 -0.34 0.77 16 43 H 0.07 1.09 8.14 -51.93 -0.42 0.67 16 44 H 0.37 9.31 6.63 -40.82 -0.27 9.04 16 45 H 0.25 7.66 6.77 -40.82 -0.28 7.39 16 46 H 0.06 1.18 8.14 -51.93 -0.42 0.76 16 47 H 0.07 1.33 8.14 -51.93 -0.42 0.91 16 48 H 0.04 0.59 8.14 -51.93 -0.42 0.17 16 49 H 0.09 1.09 8.14 -51.93 -0.42 0.66 16 50 H 0.07 0.66 8.14 -51.93 -0.42 0.23 16 51 H 0.09 0.81 6.42 -51.93 -0.33 0.48 16 52 H 0.07 0.51 8.14 -51.93 -0.42 0.09 16 53 H 0.09 0.52 5.93 -51.93 -0.31 0.21 16 LS Contribution 408.43 15.07 6.15 6.15 Total: -1.00 -33.54 408.43 -8.99 -42.53 By element: Atomic # 1 Polarization: 16.30 SS G_CDS: -9.49 Total: 6.81 kcal Atomic # 6 Polarization: -6.28 SS G_CDS: 0.09 Total: -6.19 kcal Atomic # 7 Polarization: -52.54 SS G_CDS: -0.03 Total: -52.57 kcal Atomic # 8 Polarization: -16.70 SS G_CDS: -0.63 Total: -17.33 kcal Atomic # 14 Polarization: 25.68 SS G_CDS: -5.09 Total: 20.59 kcal Total LS contribution 6.15 Total: 6.15 kcal Total: -33.54 -8.99 -42.53 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. ZINC001186041285.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 -88.120 kcal (2) G-P(sol) polarization free energy of solvation -33.542 kcal (3) G-ENP(sol) elect.-nuc.-pol. free energy of system -121.663 kcal (4) G-CDS(sol) cavity-dispersion-solvent structure free energy -8.992 kcal (5) G-P-CDS(sol) = G-P(sol) + G-CDS(sol) = (2) + (4) -42.534 kcal (6) G-S(sol) free energy of system = (1) + (5) -130.655 kcal 1SCF WAS SPECIFIED, SO GEOMETRY OPTIMIZATION WAS NOT USED Total computer time = 0.97 seconds