Wall clock time and date at job start Sat Apr 11 2020 06:09:31 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.53000 * 1 3 3 C 1.53006 * 109.47070 * 2 1 4 4 C 1.50703 * 109.46570 * 240.00644 * 2 1 3 5 5 H 1.08007 * 119.99213 * 239.99671 * 4 2 1 6 6 C 1.37875 * 120.00425 * 59.99862 * 4 2 1 7 7 N 1.31515 * 119.99686 * 0.02562 * 6 4 2 8 8 Si 1.86801 * 119.99919 * 179.97438 * 6 4 2 9 9 H 1.48493 * 109.47051 * 90.00279 * 8 6 4 10 10 H 1.48502 * 109.46814 * 209.99925 * 8 6 4 11 11 H 1.48498 * 109.47224 * 329.99639 * 8 6 4 12 12 C 1.52993 * 109.47114 * 120.00569 * 2 1 3 13 13 N 1.46502 * 109.46971 * 299.99447 * 12 2 1 14 14 C 1.34779 * 119.99811 * 179.97438 * 13 12 2 15 15 O 1.21591 * 119.99763 * 5.07839 * 14 13 12 16 16 N 1.34773 * 120.00080 * 185.06976 * 14 13 12 17 17 C 1.46928 * 120.63095 * 174.97128 * 16 14 13 18 18 C 1.53188 * 108.89312 * 233.67873 * 17 16 14 19 19 C 1.53047 * 109.31426 * 305.29667 * 18 17 16 20 20 O 1.42901 * 109.45924 * 181.37603 * 19 18 17 21 21 C 1.42897 * 114.00153 * 149.99825 * 20 19 18 22 22 H 1.09005 * 110.45462 * 330.70138 * 21 20 19 23 23 C 1.54507 * 110.51252 * 208.09105 * 21 20 19 24 24 C 1.54480 * 104.06353 * 241.35380 * 23 21 20 25 25 O 1.44426 * 104.80440 * 23.97837 * 24 23 21 26 26 C 1.44790 * 105.29727 * 319.45512 * 25 24 23 27 27 C 1.53034 * 109.53861 * 61.36442 * 19 18 17 28 28 C 1.47206 * 120.81292 * 355.11790 * 16 14 13 29 29 H 1.09004 * 109.47114 * 299.99844 * 1 2 3 30 30 H 1.09006 * 109.47514 * 59.99260 * 1 2 3 31 31 H 1.09000 * 109.47376 * 179.97438 * 1 2 3 32 32 H 1.09003 * 109.46506 * 180.02562 * 3 2 1 33 33 H 1.09000 * 109.47458 * 300.00299 * 3 2 1 34 34 H 1.08999 * 109.46897 * 60.00570 * 3 2 1 35 35 H 0.97003 * 119.99858 * 180.02562 * 7 6 4 36 36 H 1.09000 * 109.46670 * 179.97438 * 12 2 1 37 37 H 1.08996 * 109.47132 * 59.99869 * 12 2 1 38 38 H 0.96990 * 120.00214 * 0.02562 * 13 12 2 39 39 H 1.09001 * 109.58684 * 113.82376 * 17 16 14 40 40 H 1.08993 * 109.58878 * 353.53749 * 17 16 14 41 41 H 1.09003 * 109.49926 * 65.24247 * 18 17 16 42 42 H 1.09002 * 109.49814 * 185.31110 * 18 17 16 43 43 H 1.09002 * 109.45505 * 301.34337 * 19 18 17 44 44 H 1.08999 * 110.50972 * 359.97438 * 23 21 20 45 45 H 1.09002 * 110.50888 * 122.57593 * 23 21 20 46 46 H 1.08995 * 110.37611 * 142.81103 * 24 23 21 47 47 H 1.08998 * 110.37526 * 265.14802 * 24 23 21 48 48 H 1.09006 * 110.42177 * 281.77771 * 26 25 24 49 49 H 1.08995 * 110.48937 * 159.32195 * 26 25 24 50 50 H 1.09002 * 109.49962 * 58.58997 * 27 19 18 51 51 H 1.08998 * 109.50062 * 178.68687 * 27 19 18 52 52 H 1.08996 * 109.53754 * 6.52841 * 28 16 14 53 53 H 1.09001 * 109.67412 * 246.22715 * 28 16 14 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 6 2.0400 1.4426 0.0000 4 6 2.0322 -0.7103 -1.2306 5 1 2.6532 -1.5885 -1.1315 6 6 1.6991 -0.2392 -2.4828 7 7 0.9434 0.8304 -2.6033 8 14 2.3222 -1.1193 -4.0081 9 1 1.3445 -2.1612 -4.4125 10 1 2.4870 -0.1399 -5.1121 11 1 3.6306 -1.7563 -3.7122 12 6 2.0400 -0.7213 1.2491 13 7 1.5517 -0.0307 2.4453 14 6 1.8769 -0.4916 3.6694 15 8 2.6495 -1.4237 3.7826 16 7 1.3459 0.0807 4.7680 17 6 1.6033 -0.4784 6.1021 18 6 0.2638 -0.7293 6.8017 19 6 -0.5519 0.5657 6.8127 20 8 -1.7888 0.3434 7.4929 21 6 -2.3056 1.5009 8.1525 22 1 -1.4927 2.1579 8.4616 23 6 -3.1643 1.0920 9.3701 24 6 -4.5688 1.6451 9.0417 25 8 -4.5961 1.7514 7.6016 26 6 -3.2901 2.2546 7.2309 27 6 -0.8332 1.0056 5.3742 28 6 0.4939 1.2762 4.6584 29 1 -0.3633 0.5138 0.8900 30 1 -0.3634 0.5140 -0.8899 31 1 -0.3634 -1.0276 -0.0005 32 1 3.1300 1.4425 -0.0005 33 1 1.6767 1.9565 -0.8899 34 1 1.6767 1.9563 0.8900 35 1 0.7087 1.1617 -3.4843 36 1 3.1300 -0.7217 1.2488 37 1 1.6766 -1.7490 1.2490 38 1 0.9939 0.7577 2.3564 39 1 2.1941 0.2272 6.6863 40 1 2.1469 -1.4182 6.0058 41 1 -0.2875 -1.5019 6.2656 42 1 0.4438 -1.0552 7.8261 43 1 0.0106 1.3447 7.3273 44 1 -3.1950 0.0069 9.4687 45 1 -2.7759 1.5469 10.2814 46 1 -5.3386 0.9544 9.3858 47 1 -4.7066 2.6257 9.4972 48 1 -3.2313 3.3283 7.4101 49 1 -3.0778 2.0314 6.1854 50 1 -1.3741 0.2166 4.8516 51 1 -1.4339 1.9151 5.3827 52 1 0.3045 1.4936 3.6073 53 1 0.9930 2.1264 5.1232 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 ZINC001427711396.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:09:31 Heat of formation + Delta-G solvation = -118.110436 kcal Electronic energy + Delta-G solvation = -31054.804473 eV Core-core repulsion = 26926.720562 eV Total energy + Delta-G solvation = -4128.083911 eV No. of doubly occupied orbitals = 66 Molecular weight (most abundant/longest-lived isotopes) = 340.220 amu Computer time = 0.78 seconds Orbital eigenvalues (eV) -40.48329 -39.82773 -38.20426 -36.64677 -35.92168 -33.46178 -32.35759 -31.44752 -31.06159 -30.43803 -29.02311 -25.80410 -25.71696 -25.25132 -24.34168 -23.36495 -23.12842 -21.62038 -21.16411 -19.44177 -18.90135 -18.11396 -17.40397 -16.82051 -16.35207 -15.95769 -15.70832 -15.45146 -15.33067 -14.87412 -14.73765 -14.49956 -14.09267 -13.99682 -13.62516 -13.45849 -13.22416 -12.87751 -12.68991 -12.65757 -12.47545 -12.32105 -12.11722 -11.97850 -11.73048 -11.47640 -11.37441 -11.28428 -11.18656 -11.10023 -10.85685 -10.72083 -10.67670 -10.38734 -10.37378 -10.18750 -9.77072 -9.63268 -9.61779 -9.28905 -8.78078 -8.73380 -8.39376 -8.14059 -6.00954 -4.10641 2.88446 3.15820 3.29574 3.36463 3.40834 3.54509 3.61457 4.11761 4.22696 4.38344 4.41902 4.53135 4.54741 4.56242 4.65447 4.77529 4.86727 4.91466 4.98989 5.07588 5.17596 5.23870 5.25937 5.31850 5.32726 5.40052 5.49657 5.56219 5.60973 5.61933 5.79112 5.89263 6.03861 6.21012 6.33676 6.38118 6.51090 6.57346 6.80946 6.82634 6.92297 6.93470 7.04269 7.10388 7.24344 7.39713 7.56861 7.71050 8.02154 8.13328 8.37581 8.42875 8.88980 9.01904 9.55023 11.00558 Molecular weight = 340.22amu Principal moments of inertia in cm(-1) A = 0.017988 B = 0.002267 C = 0.002087 Principal moments of inertia in units of 10**(-40)*gram-cm**2 A = 1556.225278 B =12348.110555 C =13411.188364 Net atomic charges, atomic populations, and dipole contributions calculated with CM2 Atom NO. Type Charge No. of electrons 1 C -0.129 4.129 2 C 0.056 3.944 3 C -0.129 4.129 4 C -0.504 4.504 5 H 0.058 0.942 6 C 0.066 3.934 7 N -0.887 5.887 8 Si 0.894 3.106 9 H -0.281 1.281 10 H -0.287 1.287 11 H -0.274 1.274 12 C 0.149 3.851 13 N -0.733 5.733 14 C 0.703 3.297 15 O -0.597 6.597 16 N -0.637 5.637 17 C 0.135 3.865 18 C -0.131 4.131 19 C 0.091 3.909 20 O -0.363 6.363 21 C 0.056 3.944 22 H 0.077 0.923 23 C -0.155 4.155 24 C 0.044 3.956 25 O -0.369 6.369 26 C 0.019 3.981 27 C -0.160 4.160 28 C 0.115 3.885 29 H 0.004 0.996 30 H 0.085 0.915 31 H 0.027 0.973 32 H 0.027 0.973 33 H 0.085 0.915 34 H 0.004 0.996 35 H 0.259 0.741 36 H 0.053 0.947 37 H 0.051 0.949 38 H 0.393 0.607 39 H 0.062 0.938 40 H 0.099 0.901 41 H 0.077 0.923 42 H 0.080 0.920 43 H 0.063 0.937 44 H 0.090 0.910 45 H 0.085 0.915 46 H 0.094 0.906 47 H 0.060 0.940 48 H 0.066 0.934 49 H 0.102 0.898 50 H 0.081 0.919 51 H 0.076 0.924 52 H 0.084 0.916 53 H 0.068 0.932 For the charges calculated by CM2: Dipole moment (debyes) X Y Z Total from point charges -7.373 6.246 24.403 26.246 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.187 4.187 2 C 0.055 3.945 3 C -0.187 4.187 4 C -0.542 4.542 5 H 0.076 0.924 6 C -0.136 4.136 7 N -0.525 5.525 8 Si 0.722 3.278 9 H -0.207 1.207 10 H -0.213 1.213 11 H -0.199 1.199 12 C 0.026 3.974 13 N -0.389 5.389 14 C 0.404 3.596 15 O -0.476 6.476 16 N -0.376 5.376 17 C 0.011 3.989 18 C -0.169 4.169 19 C 0.034 3.966 20 O -0.282 6.282 21 C -0.003 4.003 22 H 0.095 0.905 23 C -0.193 4.193 24 C -0.033 4.033 25 O -0.287 6.287 26 C -0.058 4.058 27 C -0.198 4.198 28 C -0.009 4.009 29 H 0.023 0.977 30 H 0.104 0.896 31 H 0.046 0.954 32 H 0.046 0.954 33 H 0.104 0.896 34 H 0.023 0.977 35 H 0.070 0.930 36 H 0.071 0.929 37 H 0.070 0.930 38 H 0.230 0.770 39 H 0.080 0.920 40 H 0.117 0.883 41 H 0.096 0.904 42 H 0.098 0.902 43 H 0.081 0.919 44 H 0.109 0.891 45 H 0.103 0.897 46 H 0.112 0.888 47 H 0.078 0.922 48 H 0.085 0.915 49 H 0.120 0.880 50 H 0.099 0.901 51 H 0.094 0.906 52 H 0.103 0.897 53 H 0.086 0.914 Dipole moment (debyes) X Y Z Total from point charges -7.142 5.420 24.134 25.745 hybrid contribution 0.221 0.345 -0.749 0.853 sum -6.921 5.765 23.385 25.059 Atomic orbital electron populations 1.21820 0.97788 0.99186 0.99906 1.18366 0.92770 0.92774 0.90592 1.21818 0.98620 0.98363 0.99892 1.20981 1.31388 1.12052 0.89826 0.92398 1.24881 0.94904 0.94925 0.98850 1.69805 1.36791 1.22740 1.23193 0.86566 0.79335 0.79371 0.82571 1.20724 1.21300 1.19900 1.21599 0.98028 0.95348 0.82447 1.45678 1.55548 1.33069 1.04606 1.15959 0.80770 0.82195 0.80632 1.90872 1.39261 1.31773 1.85646 1.47348 1.53559 1.30836 1.05898 1.21327 0.95437 1.00583 0.81557 1.22037 0.95972 0.97218 1.01721 1.22216 0.84690 0.97835 0.91906 1.88040 1.38324 1.33052 1.68795 1.23592 0.96353 0.88028 0.92295 0.90508 1.23087 0.95162 1.03882 0.97123 1.23414 0.97117 1.00428 0.82312 1.89247 1.32516 1.86328 1.20650 1.24233 0.83155 0.97347 1.01086 1.22571 0.97898 1.02811 0.96535 1.21582 0.90093 0.87732 1.01479 0.97725 0.89615 0.95418 0.95404 0.89611 0.97734 0.93036 0.92853 0.93017 0.76996 0.92000 0.88261 0.90433 0.90165 0.91927 0.89099 0.89666 0.88808 0.92202 0.91543 0.87968 0.90074 0.90555 0.89748 0.91400 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.13 -2.93 8.08 37.16 0.30 -2.63 16 2 C 0.06 1.30 0.50 -155.66 -0.08 1.22 16 3 C -0.13 -2.92 8.08 37.16 0.30 -2.62 16 4 C -0.50 -13.25 8.01 -34.44 -0.28 -13.53 16 5 H 0.06 1.54 7.74 -52.48 -0.41 1.14 16 6 C 0.07 1.73 5.15 -17.82 -0.09 1.64 16 7 N -0.89 -23.64 9.35 55.43 0.52 -23.12 16 8 Si 0.89 20.82 29.54 -169.99 -5.02 15.80 16 9 H -0.28 -6.56 7.11 56.52 0.40 -6.15 16 10 H -0.29 -6.66 7.11 56.52 0.40 -6.26 16 11 H -0.27 -6.40 7.11 56.52 0.40 -5.99 16 12 C 0.15 3.24 4.51 -4.04 -0.02 3.22 16 13 N -0.73 -13.49 4.87 -65.23 -0.32 -13.81 16 14 C 0.70 12.82 8.30 -86.92 -0.72 12.10 16 15 O -0.60 -13.08 15.78 5.29 0.08 -12.99 16 16 N -0.64 -8.67 2.98 -178.28 -0.53 -9.20 16 17 C 0.13 1.52 6.41 -3.71 -0.02 1.50 16 18 C -0.13 -1.27 5.73 -26.61 -0.15 -1.42 16 19 C 0.09 0.80 2.57 -26.69 -0.07 0.74 16 20 O -0.36 -3.58 9.63 -35.23 -0.34 -3.91 16 21 C 0.06 0.42 3.57 -25.16 -0.09 0.33 16 22 H 0.08 0.43 7.42 -51.93 -0.39 0.04 16 23 C -0.15 -0.98 6.90 -25.16 -0.17 -1.16 16 24 C 0.04 0.30 7.63 37.93 0.29 0.59 16 25 O -0.37 -3.85 11.66 -35.23 -0.41 -4.26 16 26 C 0.02 0.15 7.15 37.94 0.27 0.42 16 27 C -0.16 -1.50 5.30 -26.61 -0.14 -1.64 16 28 C 0.12 1.21 6.44 -3.56 -0.02 1.18 16 29 H 0.00 0.08 6.86 -51.93 -0.36 -0.28 16 30 H 0.09 2.14 6.07 -51.93 -0.32 1.83 16 31 H 0.03 0.62 8.14 -51.93 -0.42 0.20 16 32 H 0.03 0.62 8.14 -51.93 -0.42 0.20 16 33 H 0.09 2.13 6.07 -51.93 -0.32 1.82 16 34 H 0.00 0.07 6.86 -51.93 -0.36 -0.28 16 35 H 0.26 6.69 8.31 -40.82 -0.34 6.35 16 36 H 0.05 1.23 8.14 -51.93 -0.42 0.80 16 37 H 0.05 1.18 8.14 -51.93 -0.42 0.76 16 38 H 0.39 6.53 3.64 -40.82 -0.15 6.38 16 39 H 0.06 0.58 8.14 -51.93 -0.42 0.16 16 40 H 0.10 1.31 7.03 -51.93 -0.37 0.94 16 41 H 0.08 0.90 8.14 -51.93 -0.42 0.48 16 42 H 0.08 0.66 8.14 -51.93 -0.42 0.23 16 43 H 0.06 0.45 7.39 -51.93 -0.38 0.07 16 44 H 0.09 0.66 8.00 -51.93 -0.42 0.24 16 45 H 0.08 0.36 8.14 -51.93 -0.42 -0.07 16 46 H 0.09 0.55 8.14 -51.93 -0.42 0.13 16 47 H 0.06 0.31 8.14 -51.93 -0.42 -0.11 16 48 H 0.07 0.41 8.14 -51.93 -0.42 -0.01 16 49 H 0.10 0.88 5.85 -51.93 -0.30 0.58 16 50 H 0.08 0.93 8.14 -51.93 -0.42 0.51 16 51 H 0.08 0.59 6.11 -51.93 -0.32 0.27 16 52 H 0.08 0.93 5.59 -51.93 -0.29 0.64 16 53 H 0.07 0.58 8.14 -51.93 -0.42 0.16 16 LS Contribution 398.27 15.07 6.00 6.00 Total: -1.00 -31.08 398.27 -9.70 -40.78 By element: Atomic # 1 Polarization: 13.76 SS G_CDS: -8.99 Total: 4.77 kcal Atomic # 6 Polarization: 0.64 SS G_CDS: -0.70 Total: -0.06 kcal Atomic # 7 Polarization: -45.80 SS G_CDS: -0.33 Total: -46.13 kcal Atomic # 8 Polarization: -20.50 SS G_CDS: -0.67 Total: -21.17 kcal Atomic # 14 Polarization: 20.82 SS G_CDS: -5.02 Total: 15.80 kcal Total LS contribution 6.00 Total: 6.00 kcal Total: -31.08 -9.70 -40.78 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. ZINC001427711396.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 -77.329 kcal (2) G-P(sol) polarization free energy of solvation -31.077 kcal (3) G-ENP(sol) elect.-nuc.-pol. free energy of system -108.406 kcal (4) G-CDS(sol) cavity-dispersion-solvent structure free energy -9.705 kcal (5) G-P-CDS(sol) = G-P(sol) + G-CDS(sol) = (2) + (4) -40.782 kcal (6) G-S(sol) free energy of system = (1) + (5) -118.110 kcal 1SCF WAS SPECIFIED, SO GEOMETRY OPTIMIZATION WAS NOT USED Total computer time = 0.78 seconds