Wall clock time and date at job start Sun Mar 7 2021 13:24:30 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 O 2 2 C 1.42905 * 1 3 3 Si 1.86298 * 109.47434 * 2 1 4 4 C 1.50704 * 109.47012 * 120.00478 * 2 1 3 5 5 C 1.38234 * 119.93404 * 245.01623 * 4 2 1 6 6 C 1.38200 * 120.07312 * 179.76642 * 5 4 2 7 7 C 1.38670 * 119.95524 * 0.46516 * 6 5 4 8 8 O 1.35925 * 120.05981 * 179.78893 * 7 6 5 9 9 C 1.42903 * 116.99722 * 359.72423 * 8 7 6 10 10 C 1.53000 * 109.46726 * 179.97438 * 9 8 7 11 11 N 1.46499 * 109.46840 * 179.97438 * 10 9 8 12 12 C 1.34818 * 126.59197 * 125.00117 * 11 10 9 13 13 C 1.36599 * 106.83805 * 180.02562 * 12 11 10 14 14 C 1.47520 * 125.76967 * 180.02562 * 13 12 11 15 15 O 1.21598 * 120.00232 * 359.97438 * 14 13 12 16 16 N 1.34772 * 120.00287 * 180.02562 * 14 13 12 17 17 C 1.46502 * 120.00615 * 179.97438 * 16 14 13 18 18 C 1.52992 * 109.47596 * 179.97438 * 17 16 14 19 19 H 1.09001 * 112.85089 * 316.19651 * 18 17 16 20 20 C 1.53786 * 113.61130 * 87.50143 * 18 17 16 21 21 C 1.53783 * 87.08172 * 139.98469 * 20 18 17 22 22 H 1.08995 * 113.61599 * 220.02177 * 21 20 18 23 23 N 1.46502 * 113.69216 * 89.12310 * 21 20 18 24 24 C 1.34429 * 126.82286 * 347.47237 * 23 21 20 25 25 N 1.30825 * 108.07373 * 179.97438 * 24 23 21 26 26 C 1.33238 * 109.83457 * 0.02562 * 25 24 23 27 27 N 1.30721 * 108.94699 * 359.75104 * 26 25 24 28 28 C 1.53776 * 113.61534 * 184.99853 * 18 17 16 29 29 N 1.34398 * 108.46219 * 0.02562 * 13 12 11 30 30 N 1.28767 * 109.51642 * 359.70839 * 29 13 12 31 31 C 1.38861 * 119.87901 * 359.51320 * 7 6 5 32 32 Cl 1.73603 * 120.03650 * 180.27326 * 31 7 6 33 33 C 1.38207 * 119.92328 * 0.25422 * 31 7 6 34 34 H 0.96695 * 113.99715 * 59.99412 * 1 2 3 35 35 H 1.09008 * 109.46606 * 240.00042 * 2 1 3 36 36 H 1.48500 * 110.00058 * 298.68270 * 3 2 1 37 37 H 1.48499 * 109.99738 * 60.00153 * 3 2 1 38 38 H 1.07991 * 119.96548 * 0.02562 * 5 4 2 39 39 H 1.07994 * 120.02056 * 180.25738 * 6 5 4 40 40 H 1.09001 * 109.46938 * 59.99959 * 9 8 7 41 41 H 1.08992 * 109.47388 * 299.99951 * 9 8 7 42 42 H 1.09002 * 109.47398 * 300.00244 * 10 9 8 43 43 H 1.09001 * 109.47194 * 60.00425 * 10 9 8 44 44 H 1.07997 * 126.57877 * 0.02610 * 12 11 10 45 45 H 0.97001 * 119.99603 * 359.97438 * 16 14 13 46 46 H 1.08997 * 109.46812 * 300.00223 * 17 16 14 47 47 H 1.09008 * 109.46847 * 59.99310 * 17 16 14 48 48 H 1.08999 * 113.61225 * 254.53198 * 20 18 17 49 49 H 1.08993 * 113.61414 * 25.43223 * 20 18 17 50 50 H 1.07995 * 125.95788 * 359.92191 * 24 23 21 51 51 H 1.07997 * 125.53374 * 179.97438 * 26 25 24 52 52 H 1.08993 * 113.61791 * 334.56575 * 28 18 17 53 53 H 1.09001 * 113.61670 * 105.39287 * 28 18 17 54 54 H 1.07999 * 119.97655 * 179.97438 * 33 31 7 Note: An asterisk (*) indicates this parameter was or will be optimized. Cartesian coordinates (angstroms) NO. Atom X Y Z 1 8 0.0000 0.0000 0.0000 2 6 1.4290 0.0000 0.0000 3 14 2.0501 1.7564 0.0000 4 6 1.9314 -0.7105 1.2304 5 6 2.6383 -0.0111 2.1906 6 6 3.0961 -0.6588 3.3223 7 6 2.8543 -2.0138 3.4902 8 8 3.3073 -2.6539 4.6005 9 6 4.0309 -1.8700 5.5513 10 6 4.4560 -2.7556 6.7243 11 7 5.1983 -1.9522 7.6988 12 6 4.9198 -1.8214 9.0114 13 6 5.8634 -0.9837 9.5345 14 6 5.9508 -0.5561 10.9437 15 8 5.1349 -0.9556 11.7520 16 7 6.9342 0.2781 11.3351 17 6 7.0207 0.7032 12.7345 18 6 8.2175 1.6400 12.9099 19 1 9.1105 1.2871 12.3942 20 6 7.8888 3.1171 12.6359 21 6 8.7833 3.4580 13.8395 22 1 8.3875 4.2580 14.4651 23 7 10.1922 3.6539 13.4891 24 6 10.7833 3.3418 12.3228 25 7 12.0464 3.6721 12.4068 26 6 12.2982 4.1951 13.6061 27 7 11.1925 4.2098 14.3033 28 6 8.4708 2.0489 14.3706 29 7 6.6988 -0.6148 8.5485 30 7 6.3377 -1.1750 7.4467 31 6 2.1435 -2.7145 2.5247 32 17 1.8323 -4.4095 2.7342 33 6 1.6838 -2.0608 1.3972 34 1 -0.3933 0.4418 -0.7650 35 1 1.7923 -0.5139 -0.8901 36 1 1.5880 2.4585 -1.2242 37 1 1.5617 2.4678 1.2085 38 1 2.8295 1.0435 2.0582 39 1 3.6447 -0.1106 4.0739 40 1 3.3948 -1.0639 5.9169 41 1 4.9157 -1.4472 5.0755 42 1 3.5709 -3.1781 7.2000 43 1 5.0920 -3.5619 6.3590 44 1 4.1067 -2.2892 9.5464 45 1 7.5849 0.5971 10.6903 46 1 6.1057 1.2261 13.0125 47 1 7.1476 -0.1714 13.3727 48 1 8.2781 3.4820 11.6854 49 1 6.8387 3.3662 12.7879 50 1 10.3000 2.8954 11.4664 51 1 13.2590 4.5525 13.9458 52 1 7.5798 2.0060 14.9968 53 1 9.3291 1.5503 14.8210 54 1 1.1313 -2.6050 0.6455 There are 81 doubly filled levels RHF calculation, no. of doubly occupied orbitals= 81 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) Si: (AM1): M.J.S.DEWAR, C. JIE, ORGANOMETALLICS, 6, 1486-1490 (1987). Cl: (AM1): M.J.S. DEWAR AND E. G. ZOEBISCH, THEOCHEM, 180, 1 (1988). 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_16946164_7143556.mol2 54 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 13:24:30 Heat of formation + Delta-G solvation = 151.686862 kcal Electronic energy + Delta-G solvation = -40589.448545 eV Core-core repulsion = 34981.288981 eV Total energy + Delta-G solvation = -5608.159564 eV No. of doubly occupied orbitals = 81 and no. of open levels = 1 Molecular weight (most abundant/longest-lived isotopes) = 460.146 amu Computer time = 2.16 seconds Orbital eigenvalues (eV) -43.33174 -43.19860 -41.49898 -40.78686 -39.83542 -39.39086 -37.79237 -37.23113 -36.17895 -35.71456 -33.30877 -32.57404 -32.20014 -32.04063 -31.99171 -30.79784 -28.99802 -26.91610 -26.72494 -25.74306 -25.50612 -24.20248 -23.87316 -23.04587 -22.81696 -22.36394 -22.30705 -21.50214 -20.61189 -19.78863 -18.99759 -18.58946 -17.86897 -17.85966 -17.55641 -17.35861 -17.23202 -17.05430 -16.99366 -16.67855 -16.32197 -16.08742 -15.95120 -15.77804 -15.67779 -15.56502 -15.39305 -15.09011 -15.05846 -14.96758 -14.79202 -14.52574 -14.06681 -13.96270 -13.88586 -13.78959 -13.46550 -13.42393 -13.18995 -13.14325 -12.97551 -12.77574 -12.57690 -12.50869 -12.21056 -12.10217 -12.04153 -11.97255 -11.72818 -11.69724 -11.65765 -11.46490 -11.23266 -11.03684 -10.84771 -10.56077 -10.52679 -10.49693 -10.05766 -10.02661 -9.34499 -4.93669 -0.03933 0.04336 0.07206 0.56674 0.62405 0.73888 1.05492 1.30022 1.31191 1.59051 1.93502 1.96155 1.97339 2.08166 2.20473 2.34476 3.08356 3.16445 3.22738 3.31711 3.33009 3.45099 3.49059 3.55833 3.62849 3.67341 3.71552 3.76350 3.85641 3.90565 3.90788 4.00033 4.05615 4.09338 4.13454 4.20474 4.24893 4.38838 4.41544 4.46781 4.52493 4.57186 4.62364 4.68685 4.70646 4.72937 4.79648 4.82091 4.93232 4.95194 4.96654 5.16171 5.22568 5.26947 5.32713 5.38369 5.43966 5.44315 5.91294 6.00978 6.37164 6.57408 6.82627 7.08596 7.10607 Net atomic charges, atomic populations, and dipole contributions calculated with CM2 Atom NO. Type Charge No. of electrons 1 O -0.558 6.558 2 C -0.082 4.082 3 Si 0.763 3.237 4 C -0.105 4.105 5 C -0.072 4.072 6 C -0.188 4.188 7 C 0.138 3.862 8 O -0.305 6.305 9 C 0.045 3.955 10 C 0.124 3.876 11 N -0.351 5.351 12 C 0.116 3.884 13 C -0.123 4.123 14 C 0.605 3.395 15 O -0.569 6.569 16 N -0.705 5.705 17 C 0.133 3.867 18 C -0.119 4.119 19 H 0.104 0.896 20 C -0.123 4.123 21 C 0.116 3.884 22 H 0.117 0.883 23 N -0.362 5.362 24 C 0.256 3.744 25 N -0.579 5.579 26 C 0.190 3.810 27 N -0.351 5.351 28 C -0.129 4.129 29 N -0.213 5.213 30 N -0.081 5.081 31 C -0.081 4.081 32 Cl -0.060 7.060 33 C -0.054 4.054 34 H 0.396 0.604 35 H 0.098 0.902 36 H -0.287 1.287 37 H -0.294 1.294 38 H 0.153 0.847 39 H 0.162 0.838 40 H 0.082 0.918 41 H 0.086 0.914 42 H 0.147 0.853 43 H 0.123 0.877 44 H 0.235 0.765 45 H 0.416 0.584 46 H 0.081 0.919 47 H 0.066 0.934 48 H 0.093 0.907 49 H 0.129 0.871 50 H 0.250 0.750 51 H 0.236 0.764 52 H 0.124 0.876 53 H 0.064 0.936 54 H 0.145 0.855 For the charges calculated by CM2: Dipole moment (debyes) X Y Z Total from point charges -2.114 0.614 -3.602 4.222 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 O -0.369 6.369 2 C -0.191 4.191 3 Si 0.668 3.332 4 C -0.106 4.106 5 C -0.090 4.090 6 C -0.207 4.207 7 C 0.091 3.909 8 O -0.217 6.217 9 C -0.031 4.031 10 C 0.003 3.997 11 N -0.136 5.136 12 C -0.041 4.041 13 C -0.262 4.262 14 C 0.390 3.610 15 O -0.447 6.447 16 N -0.356 5.356 17 C 0.010 3.990 18 C -0.137 4.137 19 H 0.122 0.878 20 C -0.161 4.161 21 C 0.013 3.987 22 H 0.135 0.865 23 N -0.152 5.152 24 C -0.040 4.040 25 N -0.306 5.306 26 C -0.113 4.113 27 N -0.179 5.179 28 C -0.167 4.167 29 N -0.087 5.087 30 N -0.062 5.062 31 C -0.109 4.109 32 Cl -0.032 7.032 33 C -0.072 4.072 34 H 0.247 0.753 35 H 0.116 0.884 36 H -0.214 1.214 37 H -0.222 1.222 38 H 0.170 0.830 39 H 0.179 0.821 40 H 0.100 0.900 41 H 0.104 0.896 42 H 0.165 0.835 43 H 0.141 0.859 44 H 0.251 0.749 45 H 0.253 0.747 46 H 0.099 0.901 47 H 0.084 0.916 48 H 0.112 0.888 49 H 0.147 0.853 50 H 0.266 0.734 51 H 0.252 0.748 52 H 0.142 0.858 53 H 0.083 0.917 54 H 0.163 0.837 Dipole moment (debyes) X Y Z Total from point charges -1.757 1.164 -2.988 3.656 hybrid contribution -1.582 -0.529 -0.163 1.676 sum -3.339 0.636 -3.150 4.634 Atomic orbital electron populations 1.86548 1.18459 1.80716 1.51152 1.24966 0.84140 1.11775 0.98173 0.94184 0.90249 0.72054 0.76691 1.19326 1.00944 0.93189 0.97145 1.21149 0.95274 1.00515 0.92036 1.21390 1.05832 0.92365 1.01075 1.18657 0.94469 0.91157 0.86601 1.86000 1.68928 1.37914 1.28901 1.22999 0.97161 0.94090 0.88829 1.21898 0.98375 0.91743 0.87670 1.50330 1.21197 1.34130 1.07960 1.23910 0.99997 0.96108 0.84127 1.21487 1.00204 1.08832 0.95710 1.16832 0.81123 0.78383 0.84647 1.90833 1.42357 1.57888 1.53633 1.45716 1.30462 1.47420 1.12022 1.21417 0.98431 0.98073 0.81035 1.22880 0.98353 0.94818 0.97667 0.87809 1.23556 1.01191 0.94222 0.97157 1.22469 0.75049 1.00176 1.01011 0.86530 1.50397 1.09380 1.44631 1.10810 1.26625 0.89039 0.96185 0.92103 1.70592 1.15896 1.21970 1.22188 1.26540 0.99015 0.97446 0.88280 1.78026 0.98162 1.22027 1.19680 1.23364 1.03466 0.93336 0.96577 1.74608 1.21222 1.16021 0.96831 1.81999 0.99270 1.08215 1.16682 1.20335 1.06491 0.86688 0.97353 1.98381 1.94641 1.12875 1.97264 1.20881 0.96766 0.94203 0.95387 0.75336 0.88391 1.21400 1.22171 0.83017 0.82072 0.90016 0.89615 0.83471 0.85897 0.74853 0.74697 0.90125 0.91600 0.88811 0.85334 0.73450 0.74765 0.85756 0.91749 0.83736 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 O -0.56 -0.04 12.81 -148.98 -1.91 -1.95 16 2 C -0.08 0.25 2.79 71.24 0.20 0.45 16 3 Si 0.76 -3.47 29.85 68.60 2.05 -1.42 16 4 C -0.11 0.03 5.39 -19.85 -0.11 -0.07 16 5 C -0.07 0.17 8.66 22.26 0.19 0.37 16 6 C -0.19 0.53 9.04 22.36 0.20 0.74 16 7 C 0.14 0.02 6.68 22.52 0.15 0.17 16 8 O -0.30 0.68 9.60 -92.71 -0.89 -0.21 16 9 C 0.05 -0.34 4.95 71.98 0.36 0.02 16 10 C 0.12 -1.25 6.91 86.38 0.60 -0.65 16 11 N -0.35 1.78 3.69 -347.09 -1.28 0.50 16 12 C 0.12 -0.55 11.27 83.85 0.95 0.39 16 13 C -0.12 0.00 6.99 40.82 0.29 0.28 16 14 C 0.61 0.93 7.82 86.69 0.68 1.61 16 15 O -0.57 -4.21 16.86 -4.02 -0.07 -4.28 16 16 N -0.71 2.74 5.55 -466.29 -2.59 0.15 16 17 C 0.13 -0.84 5.51 86.38 0.48 -0.36 16 18 C -0.12 1.14 3.35 -10.27 -0.03 1.11 16 19 H 0.10 -1.01 7.03 -2.39 -0.02 -1.03 16 20 C -0.12 1.38 6.89 31.12 0.21 1.59 16 21 C 0.12 -0.66 5.21 45.52 0.24 -0.43 16 22 H 0.12 -0.66 8.14 -2.39 -0.02 -0.67 16 23 N -0.36 -0.31 3.81 -328.92 -1.25 -1.56 16 24 C 0.26 -0.02 10.83 180.36 1.95 1.93 16 25 N -0.58 -3.59 11.19 -204.01 -2.28 -5.87 16 26 C 0.19 1.05 13.03 180.76 2.36 3.40 16 27 N -0.35 -2.47 12.70 -61.24 -0.78 -3.24 16 28 C -0.13 0.96 7.44 31.12 0.23 1.19 16 29 N -0.21 -0.38 12.04 -50.00 -0.60 -0.98 16 30 N -0.08 -0.08 13.65 37.02 0.51 0.43 16 31 C -0.08 -0.21 6.33 22.41 0.14 -0.07 16 32 Cl -0.06 -0.21 28.51 -2.72 -0.08 -0.29 16 33 C -0.05 -0.08 9.46 22.27 0.21 0.13 16 34 H 0.40 -3.05 8.74 -74.06 -0.65 -3.70 16 35 H 0.10 -0.68 7.96 -2.38 -0.02 -0.70 16 36 H -0.29 -0.04 7.11 99.48 0.71 0.67 16 37 H -0.29 -0.61 7.05 99.48 0.70 0.09 16 38 H 0.15 -0.95 4.86 -2.92 -0.01 -0.96 16 39 H 0.16 -1.21 6.30 -2.91 -0.02 -1.23 16 40 H 0.08 -0.71 7.65 -2.39 -0.02 -0.73 16 41 H 0.09 -0.61 7.67 -2.39 -0.02 -0.63 16 42 H 0.15 -2.22 8.08 -2.39 -0.02 -2.23 16 43 H 0.12 -1.34 8.14 -2.39 -0.02 -1.36 16 44 H 0.24 -2.57 8.06 -2.91 -0.02 -2.60 16 45 H 0.42 -2.77 8.16 -92.71 -0.76 -3.53 16 46 H 0.08 -0.51 7.84 -2.39 -0.02 -0.53 16 47 H 0.07 -0.24 8.14 -2.38 -0.02 -0.25 16 48 H 0.09 -1.16 7.04 -2.39 -0.02 -1.18 16 49 H 0.13 -1.72 7.90 -2.39 -0.02 -1.74 16 50 H 0.25 -2.25 6.76 -2.91 -0.02 -2.27 16 51 H 0.24 -0.18 8.06 -2.91 -0.02 -0.20 16 52 H 0.12 -1.28 8.10 -2.39 -0.02 -1.30 16 53 H 0.06 -0.20 8.14 -2.39 -0.02 -0.22 16 54 H 0.14 0.06 8.06 -2.91 -0.02 0.03 16 Total: 0.00 -32.95 473.79 -0.25 -33.21 By element: Atomic # 1 Polarization: -25.91 SS G_CDS: -0.36 Total: -26.28 kcal Atomic # 6 Polarization: 2.51 SS G_CDS: 9.28 Total: 11.80 kcal Atomic # 7 Polarization: -2.29 SS G_CDS: -8.28 Total: -10.57 kcal Atomic # 8 Polarization: -3.58 SS G_CDS: -2.87 Total: -6.44 kcal Atomic # 14 Polarization: -3.47 SS G_CDS: 2.05 Total: -1.42 kcal Atomic # 17 Polarization: -0.21 SS G_CDS: -0.08 Total: -0.29 kcal Total: -32.95 -0.25 -33.21 kcal The number of atoms in the molecule is 54 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_16946164_7143556.mol2 54 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 184.892 kcal (2) G-P(sol) polarization free energy of solvation -32.951 kcal (3) G-ENP(sol) elect.-nuc.-pol. free energy of system 151.941 kcal (4) G-CDS(sol) cavity-dispersion-solvent structure free energy -0.254 kcal (5) G-P-CDS(sol) = G-P(sol) + G-CDS(sol) = (2) + (4) -33.205 kcal (6) G-S(sol) free energy of system = (1) + (5) 151.687 kcal 1SCF WAS SPECIFIED, SO GEOMETRY OPTIMIZATION WAS NOT USED Total computer time = 2.16 seconds