Wall clock time and date at job start Thu Mar 4 2021 15:22:59 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 N 1.46495 * 1 3 3 N 1.28466 * 124.88795 * 2 1 4 4 C 1.31808 * 110.07398 * 180.16171 * 3 2 1 5 5 C 1.47712 * 126.51399 * 179.81974 * 4 3 2 6 6 O 1.21573 * 119.99767 * 179.68436 * 5 4 3 7 7 N 1.34777 * 120.00297 * 359.68663 * 5 4 3 8 8 C 1.46498 * 119.99828 * 180.27851 * 7 5 4 9 9 C 1.52996 * 109.47059 * 179.72339 * 8 7 5 10 10 F 1.39902 * 112.85217 * 50.32848 * 9 8 7 11 11 C 1.53784 * 113.61654 * 181.52898 * 9 8 7 12 12 C 1.53786 * 87.08016 * 139.97812 * 11 9 8 13 13 H 1.08998 * 113.61069 * 219.93201 * 12 11 9 14 14 C 1.52999 * 113.61153 * 89.11785 * 12 11 9 15 15 N 1.46500 * 109.46640 * 174.99529 * 14 12 11 16 16 C 1.34789 * 119.99592 * 180.02562 * 15 14 12 17 17 O 1.21505 * 120.00001 * 359.97438 * 16 15 14 18 18 O 1.34634 * 119.99663 * 180.02562 * 16 15 14 19 19 C 1.45199 * 116.99801 * 179.97438 * 18 16 15 20 20 C 1.52997 * 109.47204 * 60.00213 * 19 18 16 21 21 C 1.52998 * 109.46985 * 180.02562 * 19 18 16 22 22 C 1.53004 * 109.47014 * 299.99876 * 19 18 16 23 23 C 1.53789 * 113.61485 * 279.13484 * 9 8 7 24 24 C 1.41500 * 106.96327 * 359.55459 * 4 3 2 25 25 C 1.35035 * 124.88356 * 179.65707 * 2 1 3 26 26 C 1.50708 * 126.38429 * 0.02562 * 25 2 1 27 27 O 1.42896 * 109.47066 * 329.99425 * 26 25 2 28 28 Si 1.86304 * 109.46792 * 209.99495 * 26 25 2 29 29 H 1.09003 * 109.47653 * 270.33831 * 1 2 3 30 30 H 1.09006 * 109.47425 * 30.34447 * 1 2 3 31 31 H 1.09000 * 109.47198 * 150.33976 * 1 2 3 32 32 H 0.97001 * 119.99727 * 359.97438 * 7 5 4 33 33 H 1.09004 * 109.46963 * 299.72135 * 8 7 5 34 34 H 1.08994 * 109.47567 * 59.72287 * 8 7 5 35 35 H 1.08995 * 113.61200 * 254.52803 * 11 9 8 36 36 H 1.08998 * 113.61142 * 25.42744 * 11 9 8 37 37 H 1.08995 * 109.47383 * 294.99920 * 14 12 11 38 38 H 1.09006 * 109.47272 * 54.99601 * 14 12 11 39 39 H 0.96997 * 120.00569 * 359.97438 * 15 14 12 40 40 H 1.08999 * 109.47296 * 59.99563 * 20 19 18 41 41 H 1.09002 * 109.46996 * 180.02562 * 20 19 18 42 42 H 1.09007 * 109.47039 * 299.99839 * 20 19 18 43 43 H 1.08998 * 109.47422 * 60.00186 * 21 19 18 44 44 H 1.09002 * 109.47216 * 180.02562 * 21 19 18 45 45 H 1.08994 * 109.47122 * 299.99718 * 21 19 18 46 46 H 1.08999 * 109.46927 * 59.99942 * 22 19 18 47 47 H 1.09002 * 109.46889 * 180.02562 * 22 19 18 48 48 H 1.08998 * 109.47480 * 299.99867 * 22 19 18 49 49 H 1.09001 * 113.61215 * 334.56632 * 23 9 8 50 50 H 1.08995 * 113.61486 * 105.46404 * 23 9 8 51 51 H 1.08003 * 127.24464 * 180.26893 * 24 4 3 52 52 H 1.08993 * 109.46984 * 90.00084 * 26 25 2 53 53 H 0.96700 * 114.00205 * 300.00232 * 27 26 25 54 54 H 1.48494 * 109.99817 * 179.97438 * 28 26 25 55 55 H 1.48491 * 110.00345 * 301.31963 * 28 26 25 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 7 1.4649 0.0000 0.0000 3 7 2.1997 1.0538 0.0000 4 6 3.4740 0.7168 -0.0035 5 6 4.6272 1.6398 -0.0084 6 8 5.7597 1.1976 -0.0055 7 7 4.4241 2.9721 -0.0161 8 6 5.5679 3.8874 -0.0271 9 6 5.0648 5.3322 -0.0411 10 9 4.1103 5.5848 0.9501 11 6 6.1923 6.3767 -0.0910 12 6 5.3094 7.1677 -1.0707 13 1 5.8655 7.6214 -1.8910 14 6 4.3584 8.1494 -0.3830 15 7 3.6508 8.9351 -1.3969 16 6 2.7554 9.8700 -1.0212 17 8 2.5350 10.0607 0.1584 18 8 2.1047 10.5917 -1.9530 19 6 1.1544 11.5805 -1.4759 20 6 1.8790 12.6018 -0.5969 21 6 0.5187 12.2947 -2.6703 22 6 0.0636 10.8862 -0.6578 23 6 4.6488 5.8262 -1.4368 24 6 3.5271 -0.6972 0.0046 25 6 2.2372 -1.1077 0.0066 26 6 1.7532 -2.5349 0.0146 27 8 0.4951 -2.6075 0.6883 28 14 2.9942 -3.6064 0.8994 29 1 -0.3634 0.0061 1.0276 30 1 -0.3634 0.8869 -0.5192 31 1 -0.3633 -0.8930 -0.5085 32 1 3.5205 3.3249 -0.0187 33 1 6.1739 3.7223 0.8638 34 1 6.1713 3.7042 -0.9161 35 1 6.3571 6.8870 0.8579 36 1 7.1139 6.0034 -0.5376 37 1 3.6363 7.5958 0.2171 38 1 4.9297 8.8178 0.2613 39 1 3.8270 8.7832 -2.3385 40 1 2.6561 13.0963 -1.1797 41 1 1.1655 13.3437 -0.2384 42 1 2.3316 12.0926 0.2540 43 1 0.0023 11.5669 -3.2962 44 1 -0.1944 13.0372 -2.3122 45 1 1.2954 12.7892 -3.2534 46 1 0.5162 10.3770 0.1930 47 1 -0.6501 11.6284 -0.3001 48 1 -0.4527 10.1586 -1.2840 49 1 5.1477 5.3021 -2.2520 50 1 3.5694 5.8926 -1.5727 51 1 4.4108 -1.3181 0.0085 52 1 1.6370 -2.8848 -1.0110 53 1 0.5239 -2.3107 1.6082 54 1 2.5283 -5.0163 0.9147 55 1 4.3151 -3.5148 0.2272 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 m_22_16498492_13259642.mol2 55 1SCF WAS SPECIFIED, SO GEOMETRY OPTIMIZATION WAS NOT USED SCF FIELD WAS ACHIEVED Wall clock time and date at job start Thu Mar 4 2021 15:22:59 Heat of formation + Delta-G solvation = -89.359539 kcal Electronic energy + Delta-G solvation = -38008.069395 eV Core-core repulsion = 32771.217857 eV Total energy + Delta-G solvation = -5236.851539 eV No. of doubly occupied orbitals = 75 and no. of open levels = 1 Molecular weight (most abundant/longest-lived isotopes) = 399.205 amu Computer time = 1.82 seconds Orbital eigenvalues (eV) -48.89828 -43.88466 -41.84615 -40.65720 -39.82197 -38.25979 -38.03150 -36.15385 -35.50129 -34.71629 -34.41341 -32.69382 -30.46002 -28.62113 -28.39490 -28.32209 -27.40815 -27.00605 -26.56563 -24.32247 -23.30184 -22.79782 -22.23840 -21.31275 -21.01082 -20.18000 -19.31257 -19.04653 -18.00446 -17.67711 -17.63958 -17.47417 -17.35350 -17.09802 -16.94680 -16.65325 -16.26537 -16.08201 -15.82178 -15.75764 -15.63532 -15.55596 -15.42143 -15.08157 -14.93150 -14.69298 -14.52415 -14.17521 -14.14246 -14.01684 -13.70415 -13.53257 -13.46411 -13.37166 -13.23929 -13.04674 -12.94104 -12.89030 -12.62898 -12.55944 -12.43508 -12.38114 -12.29226 -12.06518 -11.87010 -11.60755 -11.42360 -11.32601 -11.22618 -11.09436 -10.94309 -10.41163 -10.24897 -10.09092 -10.00407 -4.97214 0.38745 0.77824 1.12181 1.18079 1.52232 1.69885 1.72055 1.93936 2.01576 2.23033 3.02376 3.11403 3.15183 3.21938 3.26782 3.29812 3.35545 3.42227 3.58868 3.72080 3.77096 3.92022 3.94752 3.97537 4.28060 4.31243 4.34025 4.38620 4.42330 4.45440 4.51286 4.54483 4.57886 4.58756 4.60866 4.68187 4.69669 4.71725 4.73922 4.76163 4.82417 4.85436 4.87859 4.94719 5.00030 5.08455 5.15149 5.19672 5.24393 5.24600 5.30848 5.33549 5.56063 5.64388 6.24284 6.35878 6.54705 6.76257 7.00892 7.30687 Net atomic charges, atomic populations, and dipole contributions calculated with CM2 Atom NO. Type Charge No. of electrons 1 C 0.076 3.924 2 N -0.258 5.258 3 N -0.263 5.263 4 C -0.039 4.039 5 C 0.606 3.394 6 O -0.565 6.565 7 N -0.708 5.708 8 C 0.111 3.889 9 C 0.071 3.929 10 F -0.199 7.199 11 C -0.129 4.129 12 C -0.115 4.115 13 H 0.122 0.878 14 C 0.148 3.852 15 N -0.697 5.697 16 C 0.641 3.359 17 O -0.596 6.596 18 O -0.368 6.368 19 C 0.127 3.873 20 C -0.180 4.180 21 C -0.137 4.137 22 C -0.180 4.180 23 C -0.143 4.143 24 C -0.196 4.196 25 C -0.006 4.006 26 C -0.035 4.035 27 O -0.543 6.543 28 Si 0.787 3.213 29 H 0.109 0.891 30 H 0.098 0.902 31 H 0.117 0.883 32 H 0.413 0.587 33 H 0.076 0.924 34 H 0.104 0.896 35 H 0.087 0.913 36 H 0.135 0.865 37 H 0.063 0.937 38 H 0.068 0.932 39 H 0.418 0.582 40 H 0.063 0.937 41 H 0.087 0.913 42 H 0.062 0.938 43 H 0.069 0.931 44 H 0.098 0.902 45 H 0.070 0.930 46 H 0.061 0.939 47 H 0.087 0.913 48 H 0.062 0.938 49 H 0.126 0.874 50 H 0.077 0.923 51 H 0.162 0.838 52 H 0.129 0.871 53 H 0.398 0.602 54 H -0.280 1.280 55 H -0.290 1.290 For the charges calculated by CM2: Dipole moment (debyes) X Y Z Total from point charges -2.830 0.199 -2.430 3.735 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.062 4.062 2 N -0.047 5.047 3 N -0.101 5.101 4 C -0.197 4.197 5 C 0.391 3.609 6 O -0.442 6.442 7 N -0.359 5.359 8 C -0.013 4.013 9 C 0.052 3.948 10 F -0.180 7.180 11 C -0.168 4.168 12 C -0.133 4.133 13 H 0.140 0.860 14 C 0.026 3.974 15 N -0.357 5.357 16 C 0.395 3.605 17 O -0.487 6.487 18 O -0.283 6.283 19 C 0.090 3.910 20 C -0.237 4.237 21 C -0.194 4.194 22 C -0.237 4.237 23 C -0.181 4.181 24 C -0.224 4.224 25 C -0.131 4.131 26 C -0.146 4.146 27 O -0.351 6.351 28 Si 0.690 3.310 29 H 0.127 0.873 30 H 0.116 0.884 31 H 0.135 0.865 32 H 0.249 0.751 33 H 0.094 0.906 34 H 0.122 0.878 35 H 0.106 0.894 36 H 0.154 0.846 37 H 0.082 0.918 38 H 0.086 0.914 39 H 0.255 0.745 40 H 0.082 0.918 41 H 0.106 0.894 42 H 0.081 0.919 43 H 0.088 0.912 44 H 0.117 0.883 45 H 0.089 0.911 46 H 0.080 0.920 47 H 0.106 0.894 48 H 0.081 0.919 49 H 0.144 0.856 50 H 0.096 0.904 51 H 0.180 0.820 52 H 0.146 0.854 53 H 0.249 0.751 54 H -0.207 1.207 55 H -0.216 1.216 Dipole moment (debyes) X Y Z Total from point charges -2.724 0.107 -2.453 3.667 hybrid contribution -0.119 -0.230 -0.540 0.599 sum -2.843 -0.123 -2.993 4.130 Atomic orbital electron populations 1.22671 0.70376 1.06565 1.06577 1.45722 1.09078 1.01663 1.48281 1.73498 1.02620 1.14805 1.19209 1.21741 0.90985 0.94238 1.12761 1.16992 0.86426 0.81867 0.75609 1.90793 1.23754 1.77723 1.51964 1.45653 1.11277 1.04868 1.74095 1.21978 0.88632 0.85985 1.04670 1.22400 0.89992 0.94269 0.88140 1.92848 1.66166 1.95024 1.63933 1.23592 0.97478 0.94317 1.01396 1.22495 0.96632 0.94106 1.00090 0.85959 1.21269 0.94253 0.92396 0.89525 1.44642 1.41765 1.37307 1.12015 1.18360 0.79279 0.79573 0.83324 1.90915 1.68537 1.70957 1.18300 1.86305 1.53921 1.49580 1.38497 1.22362 0.86629 0.85861 0.96107 1.22473 1.01721 0.99039 1.00491 1.21887 1.01187 1.00478 0.95852 1.22476 0.98432 1.01984 1.00855 1.23781 1.02817 0.96091 0.95398 1.21619 0.97350 0.93613 1.09790 1.22165 0.90612 0.92292 1.08071 1.24821 0.94683 0.91923 1.03198 1.86476 1.28836 1.88810 1.30953 0.93156 0.73561 0.75657 0.88618 0.87288 0.88362 0.86462 0.75056 0.90589 0.87776 0.89427 0.84642 0.91839 0.91368 0.74479 0.91818 0.89429 0.91938 0.91157 0.88295 0.91091 0.91986 0.89402 0.91939 0.85605 0.90384 0.82032 0.85411 0.75143 1.20701 1.21623 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.08 -0.55 10.12 127.77 1.29 0.75 16 2 N -0.26 -0.10 3.70 -355.92 -1.32 -1.41 16 3 N -0.26 -1.34 11.67 -48.18 -0.56 -1.90 16 4 C -0.04 -0.31 7.13 42.46 0.30 0.00 16 5 C 0.61 5.58 7.82 86.75 0.68 6.25 16 6 O -0.56 -8.53 16.93 -3.95 -0.07 -8.60 16 7 N -0.71 -2.51 5.47 -466.11 -2.55 -5.06 16 8 C 0.11 -0.06 5.35 86.38 0.46 0.40 16 9 C 0.07 -0.22 1.63 -10.26 -0.02 -0.23 16 10 F -0.20 -0.70 15.08 44.97 0.68 -0.02 16 11 C -0.13 0.93 7.13 31.12 0.22 1.15 16 12 C -0.11 0.82 3.95 -10.21 -0.04 0.78 16 13 H 0.12 -1.36 8.14 -2.39 -0.02 -1.37 16 14 C 0.15 -0.06 5.28 86.38 0.46 0.40 16 15 N -0.70 -0.42 5.56 -590.84 -3.29 -3.71 16 16 C 0.64 4.88 8.06 129.79 1.05 5.92 16 17 O -0.60 -8.16 11.88 19.83 0.24 -7.92 16 18 O -0.37 -2.19 9.94 -55.43 -0.55 -2.74 16 19 C 0.13 0.44 1.13 -10.79 -0.01 0.43 16 20 C -0.18 -0.68 8.37 71.98 0.60 -0.08 16 21 C -0.14 0.26 8.85 71.98 0.64 0.90 16 22 C -0.18 -0.73 8.37 71.98 0.60 -0.12 16 23 C -0.14 0.94 7.08 31.18 0.22 1.16 16 24 C -0.20 -1.04 9.35 22.41 0.21 -0.83 16 25 C -0.01 0.00 6.18 41.40 0.26 0.25 16 26 C -0.03 0.15 2.88 71.24 0.20 0.35 16 27 O -0.54 3.81 10.40 -148.98 -1.55 2.26 16 28 Si 0.79 -2.83 31.37 68.60 2.15 -0.68 16 29 H 0.11 -1.20 8.14 -2.39 -0.02 -1.22 16 30 H 0.10 -0.76 8.14 -2.38 -0.02 -0.78 16 31 H 0.12 -1.00 6.77 -2.39 -0.02 -1.02 16 32 H 0.41 0.98 7.96 -92.71 -0.74 0.24 16 33 H 0.08 0.10 8.14 -2.38 -0.02 0.08 16 34 H 0.10 -0.30 7.98 -2.39 -0.02 -0.32 16 35 H 0.09 -0.49 8.10 -2.39 -0.02 -0.51 16 36 H 0.14 -1.56 8.10 -2.39 -0.02 -1.58 16 37 H 0.06 0.27 6.31 -2.39 -0.02 0.26 16 38 H 0.07 0.10 8.14 -2.38 -0.02 0.08 16 39 H 0.42 -1.96 8.58 -92.71 -0.80 -2.76 16 40 H 0.06 0.18 8.14 -2.39 -0.02 0.16 16 41 H 0.09 0.02 8.14 -2.39 -0.02 0.00 16 42 H 0.06 0.59 5.88 -2.38 -0.01 0.58 16 43 H 0.07 -0.15 8.14 -2.39 -0.02 -0.17 16 44 H 0.10 -0.51 8.14 -2.39 -0.02 -0.53 16 45 H 0.07 -0.16 8.14 -2.39 -0.02 -0.18 16 46 H 0.06 0.61 5.88 -2.39 -0.01 0.59 16 47 H 0.09 0.03 8.14 -2.39 -0.02 0.01 16 48 H 0.06 0.20 8.14 -2.39 -0.02 0.18 16 49 H 0.13 -1.25 7.99 -2.39 -0.02 -1.27 16 50 H 0.08 -0.33 8.10 -2.39 -0.02 -0.35 16 51 H 0.16 0.87 6.84 -2.91 -0.02 0.85 16 52 H 0.13 -0.58 7.96 -2.39 -0.02 -0.60 16 53 H 0.40 -5.54 8.70 -74.06 -0.64 -6.18 16 54 H -0.28 -0.53 7.11 99.48 0.71 0.18 16 55 H -0.29 -0.90 6.83 99.48 0.68 -0.22 16 Total: 0.00 -27.22 447.46 -0.91 -28.13 By element: Atomic # 1 Polarization: -14.62 SS G_CDS: -1.22 Total: -15.84 kcal Atomic # 6 Polarization: 10.37 SS G_CDS: 7.12 Total: 17.49 kcal Atomic # 7 Polarization: -4.36 SS G_CDS: -7.71 Total: -12.08 kcal Atomic # 8 Polarization: -15.06 SS G_CDS: -1.93 Total: -16.99 kcal Atomic # 9 Polarization: -0.70 SS G_CDS: 0.68 Total: -0.02 kcal Atomic # 14 Polarization: -2.83 SS G_CDS: 2.15 Total: -0.68 kcal Total: -27.22 -0.91 -28.13 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. m_22_16498492_13259642.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 -61.234 kcal (2) G-P(sol) polarization free energy of solvation -27.217 kcal (3) G-ENP(sol) elect.-nuc.-pol. free energy of system -88.451 kcal (4) G-CDS(sol) cavity-dispersion-solvent structure free energy -0.909 kcal (5) G-P-CDS(sol) = G-P(sol) + G-CDS(sol) = (2) + (4) -28.126 kcal (6) G-S(sol) free energy of system = (1) + (5) -89.360 kcal 1SCF WAS SPECIFIED, SO GEOMETRY OPTIMIZATION WAS NOT USED Total computer time = 1.82 seconds