Wall clock time and date at job start Sun Mar 7 2021 00:05:52 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= 0 * AM1 - THE AM1 HAMILTONIAN TO BE USED * SM5.42R - SM5.42R CALCULATIONS WILL BE PERFORMED * - CHARGE MODEL 2 WILL BE USED (CM2A) * - THE SOLVENT IS WATER ******************************************************************************* 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 C 1.52997 * 109.47113 * 2 1 4 4 C 1.53001 * 109.47421 * 294.25543 * 3 2 1 5 5 C 1.53003 * 109.47317 * 174.25620 * 3 2 1 6 6 C 1.50699 * 109.47151 * 179.23869 * 5 3 2 7 7 O 1.21287 * 119.99854 * 236.33201 * 6 5 3 8 8 N 1.34776 * 120.00570 * 56.59879 * 6 5 3 9 9 C 1.46499 * 120.00335 * 184.68876 * 8 6 5 10 10 H 1.08995 * 112.84724 * 336.19552 * 9 8 6 11 11 C 1.53777 * 113.61270 * 107.50763 * 9 8 6 12 12 C 1.53778 * 87.08051 * 89.11543 * 11 9 8 13 13 H 1.09003 * 113.61526 * 270.87956 * 12 11 9 14 14 C 1.50705 * 113.61420 * 139.97879 * 12 11 9 15 15 C 1.38462 * 119.98746 * 142.50109 * 14 12 11 16 16 F 1.35095 * 120.00888 * 359.96005 * 15 14 12 17 17 C 1.38431 * 119.97498 * 179.97438 * 15 14 12 18 18 C 1.38404 * 119.95144 * 0.02562 * 17 15 14 19 19 F 1.35101 * 120.01717 * 179.97438 * 18 17 15 20 20 C 1.38495 * 119.97162 * 359.97438 * 18 17 15 21 21 C 1.38190 * 120.02850 * 359.75116 * 20 18 17 22 22 C 1.53780 * 87.08356 * 25.43364 * 12 11 9 23 23 N 1.46499 * 109.47335 * 59.23546 * 5 3 2 24 24 C 1.36341 * 125.63422 * 66.08881 * 23 5 3 25 25 C 1.34854 * 108.21566 * 179.83341 * 24 23 5 26 26 C 1.40893 * 107.41729 * 0.42559 * 25 24 23 27 27 C 1.34855 * 107.41468 * 359.73655 * 26 25 24 28 28 C 1.50697 * 125.88588 * 179.97438 * 27 26 25 29 29 O 1.42906 * 109.46825 * 150.02532 * 28 27 26 30 30 Si 1.86298 * 109.47294 * 30.02569 * 28 27 26 31 31 H 1.09003 * 109.46998 * 299.99616 * 1 2 3 32 32 H 1.08995 * 109.47217 * 59.99910 * 1 2 3 33 33 H 1.08997 * 109.46970 * 180.02562 * 1 2 3 34 34 H 1.09003 * 109.46998 * 240.00384 * 2 1 3 35 35 H 1.08995 * 109.47217 * 120.00090 * 2 1 3 36 36 H 1.08997 * 109.47222 * 54.24883 * 3 2 1 37 37 H 1.08999 * 109.47804 * 187.81315 * 4 3 2 38 38 H 1.09004 * 109.46721 * 307.81143 * 4 3 2 39 39 H 1.08997 * 109.47309 * 67.80462 * 4 3 2 40 40 H 1.09002 * 109.46519 * 299.23158 * 5 3 2 41 41 H 0.97000 * 119.99470 * 4.69751 * 8 6 5 42 42 H 1.08999 * 113.61084 * 203.66443 * 11 9 8 43 43 H 1.08993 * 113.61496 * 334.56749 * 11 9 8 44 44 H 1.07995 * 120.02850 * 180.02562 * 17 15 14 45 45 H 1.07998 * 119.98584 * 179.97438 * 20 18 17 46 46 H 1.07996 * 119.98016 * 180.25041 * 21 20 18 47 47 H 1.09008 * 113.61378 * 220.02199 * 22 12 11 48 48 H 1.08995 * 113.69450 * 89.12552 * 22 12 11 49 49 H 1.08000 * 125.89221 * 0.02562 * 24 23 5 50 50 H 1.07999 * 126.29303 * 180.16361 * 25 24 23 51 51 H 1.07991 * 126.29262 * 179.73451 * 26 25 24 52 52 H 1.08999 * 109.47298 * 270.02384 * 28 27 26 53 53 H 0.96694 * 113.99899 * 299.99466 * 29 28 27 54 54 H 1.48504 * 109.99856 * 179.97438 * 30 28 27 55 55 H 1.48502 * 109.99964 * 301.31433 * 30 28 27 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.4425 0.0000 4 6 1.6513 2.1209 -1.3151 5 6 3.5632 1.4449 0.1444 6 6 4.0673 2.8650 0.1256 7 8 4.8986 3.2026 -0.6905 8 7 3.5984 3.7585 1.0190 9 6 4.0101 5.1620 0.9366 10 1 4.9622 5.2943 0.4227 11 6 2.9008 6.1009 0.4338 12 6 2.4752 6.2566 1.9033 13 1 1.7492 5.5117 2.2289 14 6 2.0776 7.6631 2.2707 15 6 1.0265 7.8788 3.1458 16 9 0.3569 6.8288 3.6695 17 6 0.6618 9.1708 3.4837 18 6 1.3481 10.2455 2.9457 19 9 0.9930 11.5070 3.2741 20 6 2.3991 10.0273 2.0706 21 6 2.7660 8.7371 1.7385 22 6 3.9254 5.9087 2.2783 23 7 4.1603 0.7000 -0.9668 24 6 4.1248 1.0673 -2.2794 25 6 4.7833 0.1340 -2.9962 26 6 5.2336 -0.8503 -2.0944 27 6 4.8382 -0.4780 -0.8600 28 6 5.1000 -1.2341 0.4169 29 8 5.1896 -0.3131 1.5060 30 14 6.7014 -2.1739 0.2651 31 1 -0.3633 0.5138 0.8900 32 1 -0.3633 0.5138 -0.8899 33 1 -0.3633 -1.0276 0.0005 34 1 1.8933 -0.5138 -0.8900 35 1 1.8934 -0.5138 0.8899 36 1 1.5951 1.9852 0.8340 37 1 2.1305 3.0981 -1.3755 38 1 1.9776 1.5040 -2.1525 39 1 0.5691 2.2438 -1.3546 40 1 3.8394 0.9738 1.0877 41 1 2.9883 3.4722 1.7167 42 1 3.2770 7.0218 -0.0118 43 1 2.1554 5.6032 -0.1863 44 1 -0.1577 9.3400 4.1664 45 1 2.9346 10.8658 1.6504 46 1 3.5883 8.5678 1.0592 47 1 4.0070 5.2557 3.1473 48 1 4.5827 6.7761 2.3375 49 1 3.6535 1.9544 -2.6761 50 1 4.9338 0.1380 -4.0657 51 1 5.7927 -1.7395 -2.3455 52 1 4.2841 -1.9336 0.5989 53 1 5.8960 0.3402 1.4113 54 1 6.9713 -2.9234 1.5184 55 1 6.6325 -3.1125 -0.8837 There are 75 doubly filled levels RHF calculation, no. of doubly occupied orbitals= 75 No. of singly occupied orbitals= 1 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) F: (AM1): M.J.S. DEWAR AND E. G. ZOEBISCH, THEOCHEM, 180, 1 (1988). 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=0 AM1 1SCF TLIMIT=15 GEO-OK SM5.42R & SOLVNT=WATER s_22_13428180_12360566.mol2 55 1SCF WAS SPECIFIED, SO GEOMETRY OPTIMIZATION WAS NOT USED SCF FIELD WAS ACHIEVED Wall clock time and date at job start Sun Mar 7 2021 00:05:52 Heat of formation + Delta-G solvation = -40.817761 kcal Electronic energy + Delta-G solvation = -40028.048461 eV Core-core repulsion = 34862.923280 eV Total energy + Delta-G solvation = -5165.125181 eV No. of doubly occupied orbitals = 75 and no. of open levels = 1 Molecular weight (most abundant/longest-lived isotopes) = 405.190 amu Computer time = 1.90 seconds Orbital eigenvalues (eV) -50.07382 -50.04620 -42.85920 -41.44651 -40.21073 -38.73709 -37.82432 -37.65603 -34.58264 -33.32639 -31.89940 -31.48342 -31.22157 -30.34553 -28.25018 -26.70976 -26.27227 -25.88929 -25.17116 -24.23748 -22.96584 -22.17763 -21.55765 -21.48970 -20.37317 -20.19613 -19.35080 -18.71076 -18.22587 -18.10626 -17.38940 -17.05362 -16.99525 -16.83559 -16.68896 -16.53017 -16.37301 -16.32183 -16.07428 -15.53954 -15.45603 -15.34117 -15.28882 -15.10383 -14.68245 -14.64967 -14.39067 -14.09169 -13.98301 -13.95196 -13.76656 -13.72702 -13.55536 -13.42890 -13.28466 -13.08892 -12.98548 -12.82089 -12.56226 -12.49670 -12.10446 -12.05056 -12.01907 -11.97058 -11.79966 -11.70619 -11.59863 -11.30232 -11.03984 -10.58750 -10.14794 -10.06936 -9.59361 -9.51173 -8.85997 -4.92132 -0.12587 -0.03449 1.04417 1.38037 1.43989 1.63838 1.79398 1.94182 2.14373 2.28081 2.29761 2.82315 3.23382 3.30696 3.36543 3.45355 3.49988 3.58844 3.63693 3.66666 3.69410 3.81823 3.87805 3.94433 3.95394 4.08581 4.17219 4.18619 4.25373 4.38460 4.44075 4.50723 4.53972 4.61718 4.62173 4.68308 4.70988 4.74103 4.75736 4.83346 4.85588 4.89528 4.91567 4.92785 4.96031 5.01930 5.06644 5.09038 5.11872 5.12816 5.26601 5.29774 5.33226 5.38358 5.40087 5.57798 5.72693 6.03739 6.27435 6.30933 6.49801 6.52139 7.09616 Net atomic charges, atomic populations, and dipole contributions calculated with CM2 Atom NO. Type Charge No. of electrons 1 C -0.151 4.151 2 C -0.111 4.111 3 C -0.088 4.088 4 C -0.160 4.160 5 C 0.101 3.899 6 C 0.516 3.484 7 O -0.555 6.555 8 N -0.696 5.696 9 C 0.119 3.881 10 H 0.095 0.905 11 C -0.138 4.138 12 C -0.071 4.071 13 H 0.104 0.896 14 C -0.111 4.111 15 C 0.116 3.884 16 F -0.137 7.137 17 C -0.199 4.199 18 C 0.106 3.894 19 F -0.143 7.143 20 C -0.155 4.155 21 C -0.056 4.056 22 C -0.124 4.124 23 N -0.406 5.406 24 C 0.016 3.984 25 C -0.224 4.224 26 C -0.208 4.208 27 C -0.017 4.017 28 C -0.040 4.040 29 O -0.547 6.547 30 Si 0.785 3.215 31 H 0.065 0.935 32 H 0.055 0.945 33 H 0.051 0.949 34 H 0.048 0.952 35 H 0.067 0.933 36 H 0.112 0.888 37 H 0.049 0.951 38 H 0.058 0.942 39 H 0.072 0.928 40 H 0.164 0.836 41 H 0.414 0.586 42 H 0.107 0.893 43 H 0.070 0.930 44 H 0.149 0.851 45 H 0.161 0.839 46 H 0.171 0.829 47 H 0.086 0.914 48 H 0.110 0.890 49 H 0.158 0.842 50 H 0.137 0.863 51 H 0.138 0.862 52 H 0.117 0.883 53 H 0.392 0.608 54 H -0.287 1.287 55 H -0.286 1.286 For the charges calculated by CM2: Dipole moment (debyes) X Y Z Total from point charges 0.319 1.909 2.384 3.071 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.208 4.208 2 C -0.149 4.149 3 C -0.107 4.107 4 C -0.217 4.217 5 C -0.004 4.004 6 C 0.301 3.699 7 O -0.432 6.432 8 N -0.348 5.348 9 C 0.014 3.986 10 H 0.113 0.887 11 C -0.176 4.176 12 C -0.090 4.090 13 H 0.122 0.878 14 C -0.112 4.112 15 C 0.097 3.903 16 F -0.115 7.115 17 C -0.219 4.219 18 C 0.087 3.913 19 F -0.122 7.122 20 C -0.174 4.174 21 C -0.073 4.073 22 C -0.162 4.162 23 N -0.082 5.082 24 C -0.115 4.115 25 C -0.248 4.248 26 C -0.233 4.233 27 C -0.128 4.128 28 C -0.151 4.151 29 O -0.355 6.355 30 Si 0.689 3.311 31 H 0.084 0.916 32 H 0.074 0.926 33 H 0.071 0.929 34 H 0.067 0.933 35 H 0.086 0.914 36 H 0.130 0.870 37 H 0.069 0.931 38 H 0.077 0.923 39 H 0.090 0.910 40 H 0.181 0.819 41 H 0.251 0.749 42 H 0.126 0.874 43 H 0.089 0.911 44 H 0.167 0.833 45 H 0.179 0.821 46 H 0.189 0.811 47 H 0.104 0.896 48 H 0.129 0.871 49 H 0.176 0.824 50 H 0.155 0.845 51 H 0.156 0.844 52 H 0.134 0.866 53 H 0.242 0.758 54 H -0.214 1.214 55 H -0.213 1.213 Dipole moment (debyes) X Y Z Total from point charges 0.600 1.230 2.344 2.714 hybrid contribution -1.475 0.105 0.143 1.486 sum -0.875 1.335 2.487 2.955 Atomic orbital electron populations 1.21780 0.94755 1.01918 1.02366 1.21465 0.96603 0.94942 1.01886 1.21514 0.90469 0.98273 1.00424 1.22173 1.02560 0.99761 0.97201 1.22045 0.95033 0.90139 0.93222 1.20663 0.80980 0.87395 0.80817 1.90791 1.33099 1.80062 1.39289 1.45502 1.49119 1.07805 1.32376 1.21836 0.98641 0.79717 0.98433 0.88718 1.23584 0.98330 1.00613 0.95089 1.21553 0.98033 0.91884 0.97502 0.87782 1.19989 0.95540 0.97717 0.97924 1.17306 0.93404 0.83010 0.96601 1.91501 1.79274 1.55750 1.85016 1.21169 1.04294 0.91021 1.05416 1.17581 0.96816 0.78557 0.98358 1.91524 1.91404 1.37866 1.91368 1.21520 0.97147 0.99981 0.98766 1.21067 0.98957 0.89514 0.97777 1.23418 0.94874 1.01036 0.96823 1.43471 1.41577 1.17475 1.05670 1.21828 1.05442 1.01233 0.82984 1.21005 1.06569 0.97467 0.99724 1.20770 1.08222 1.02154 0.92157 1.20835 1.04023 0.90960 0.97011 1.24821 1.11946 0.95717 0.82581 1.86308 1.58670 1.47315 1.43157 0.93016 0.81301 0.80300 0.76457 0.91593 0.92607 0.92945 0.93344 0.91450 0.87033 0.93150 0.92295 0.90952 0.81931 0.74860 0.87427 0.91148 0.83274 0.82126 0.81120 0.89570 0.87125 0.82443 0.84506 0.84405 0.86608 0.75808 1.21399 1.21288 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.15 0.50 9.69 71.98 0.70 1.20 16 2 C -0.11 0.28 5.14 30.59 0.16 0.44 16 3 C -0.09 0.42 2.14 -10.79 -0.02 0.40 16 4 C -0.16 0.39 7.32 71.98 0.53 0.92 16 5 C 0.10 -0.27 2.44 44.25 0.11 -0.16 16 6 C 0.52 -0.53 5.68 87.66 0.50 -0.04 16 7 O -0.56 -4.10 14.85 -3.05 -0.05 -4.14 16 8 N -0.70 4.76 5.10 -440.04 -2.24 2.51 16 9 C 0.12 -0.79 4.36 45.52 0.20 -0.59 16 10 H 0.09 -0.26 7.42 -2.39 -0.02 -0.28 16 11 C -0.14 1.04 7.25 31.12 0.23 1.27 16 12 C -0.07 0.55 3.78 -11.01 -0.04 0.50 16 13 H 0.10 -0.77 7.45 -2.39 -0.02 -0.79 16 14 C -0.11 0.49 5.30 -19.81 -0.10 0.39 16 15 C 0.12 -0.04 7.24 22.37 0.16 0.12 16 16 F -0.14 -0.20 16.36 44.97 0.74 0.54 16 17 C -0.20 -0.10 10.01 22.36 0.22 0.12 16 18 C 0.11 0.06 7.30 22.37 0.16 0.22 16 19 F -0.14 -0.56 17.27 44.97 0.78 0.22 16 20 C -0.16 0.63 10.02 22.32 0.22 0.85 16 21 C -0.06 0.37 8.60 22.26 0.19 0.56 16 22 C -0.12 1.16 7.29 31.12 0.23 1.39 16 23 N -0.41 -0.89 1.75 -491.51 -0.86 -1.75 16 24 C 0.02 0.09 7.50 83.21 0.62 0.71 16 25 C -0.22 -1.54 10.86 22.14 0.24 -1.30 16 26 C -0.21 -1.16 9.43 22.14 0.21 -0.95 16 27 C -0.02 -0.05 6.16 41.08 0.25 0.21 16 28 C -0.04 0.05 2.99 71.23 0.21 0.26 16 29 O -0.55 1.96 9.79 -148.98 -1.46 0.50 16 30 Si 0.78 -0.99 31.31 68.60 2.15 1.16 16 31 H 0.07 -0.34 8.14 -2.39 -0.02 -0.36 16 32 H 0.05 -0.15 6.77 -2.39 -0.02 -0.17 16 33 H 0.05 -0.11 8.14 -2.39 -0.02 -0.13 16 34 H 0.05 0.10 7.76 -2.39 -0.02 0.09 16 35 H 0.07 -0.18 8.14 -2.39 -0.02 -0.20 16 36 H 0.11 -1.03 7.50 -2.39 -0.02 -1.04 16 37 H 0.05 -0.06 6.91 -2.39 -0.02 -0.07 16 38 H 0.06 0.05 5.40 -2.39 -0.01 0.04 16 39 H 0.07 -0.26 6.98 -2.39 -0.02 -0.28 16 40 H 0.16 -0.76 5.02 -2.39 -0.01 -0.77 16 41 H 0.41 -4.70 8.07 -92.71 -0.75 -5.45 16 42 H 0.11 -1.01 6.32 -2.39 -0.02 -1.03 16 43 H 0.07 -0.39 8.14 -2.39 -0.02 -0.41 16 44 H 0.15 -0.11 8.06 -2.91 -0.02 -0.13 16 45 H 0.16 -1.24 8.06 -2.91 -0.02 -1.26 16 46 H 0.17 -2.08 6.67 -2.91 -0.02 -2.10 16 47 H 0.09 -0.85 8.14 -2.38 -0.02 -0.87 16 48 H 0.11 -1.19 7.72 -2.39 -0.02 -1.21 16 49 H 0.16 0.82 6.48 -2.91 -0.02 0.80 16 50 H 0.14 0.71 8.06 -2.91 -0.02 0.69 16 51 H 0.14 0.60 6.67 -2.92 -0.02 0.58 16 52 H 0.12 -0.13 7.96 -2.39 -0.02 -0.15 16 53 H 0.39 -3.24 8.70 -74.06 -0.64 -3.89 16 54 H -0.29 -1.20 7.11 99.48 0.71 -0.49 16 55 H -0.29 -1.23 6.78 99.48 0.67 -0.56 16 Total: 0.00 -17.45 435.46 3.59 -13.86 By element: Atomic # 1 Polarization: -18.99 SS G_CDS: -0.43 Total: -19.43 kcal Atomic # 6 Polarization: 1.56 SS G_CDS: 4.97 Total: 6.53 kcal Atomic # 7 Polarization: 3.87 SS G_CDS: -3.10 Total: 0.76 kcal Atomic # 8 Polarization: -2.13 SS G_CDS: -1.50 Total: -3.64 kcal Atomic # 9 Polarization: -0.76 SS G_CDS: 1.51 Total: 0.76 kcal Atomic # 14 Polarization: -0.99 SS G_CDS: 2.15 Total: 1.16 kcal Total: -17.45 3.59 -13.86 kcal The number of atoms in the molecule is 55 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. s_22_13428180_12360566.mol2 55 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 -26.959 kcal (2) G-P(sol) polarization free energy of solvation -17.451 kcal (3) G-ENP(sol) elect.-nuc.-pol. free energy of system -44.410 kcal (4) G-CDS(sol) cavity-dispersion-solvent structure free energy 3.592 kcal (5) G-P-CDS(sol) = G-P(sol) + G-CDS(sol) = (2) + (4) -13.859 kcal (6) G-S(sol) free energy of system = (1) + (5) -40.818 kcal 1SCF WAS SPECIFIED, SO GEOMETRY OPTIMIZATION WAS NOT USED Total computer time = 1.90 seconds