Wall clock time and date at job start Thu Mar 4 2021 15:30:07 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.46494 * 1 3 3 N 1.28466 * 124.89068 * 2 1 4 4 C 1.31805 * 110.07774 * 179.97438 * 3 2 1 5 5 C 1.47703 * 126.52039 * 179.97438 * 4 3 2 6 6 O 1.21570 * 120.00016 * 180.02562 * 5 4 3 7 7 N 1.34778 * 120.00400 * 0.02562 * 5 4 3 8 8 C 1.46493 * 120.00552 * 180.02562 * 7 5 4 9 9 C 1.53005 * 109.47192 * 179.97438 * 8 7 5 10 10 F 1.39897 * 112.84883 * 48.80191 * 9 8 7 11 11 C 1.53786 * 113.61119 * 179.97438 * 9 8 7 12 12 C 1.53780 * 87.07735 * 139.98301 * 11 9 8 13 13 H 1.09000 * 113.61818 * 89.19676 * 12 11 9 14 14 C 1.53002 * 113.60890 * 220.01716 * 12 11 9 15 15 N 1.46504 * 109.46793 * 184.99857 * 14 12 11 16 16 C 1.34779 * 119.99823 * 179.97438 * 15 14 12 17 17 O 1.21510 * 120.00026 * 0.02562 * 16 15 14 18 18 O 1.34638 * 119.99882 * 179.97438 * 16 15 14 19 19 C 1.45197 * 116.99871 * 180.02562 * 18 16 15 20 20 C 1.52998 * 109.47326 * 60.00172 * 19 18 16 21 21 C 1.53002 * 109.47286 * 179.97438 * 19 18 16 22 22 C 1.52998 * 109.47503 * 300.00083 * 19 18 16 23 23 C 1.53775 * 113.61228 * 277.49420 * 9 8 7 24 24 C 1.41503 * 106.96348 * 359.76134 * 4 3 2 25 25 C 1.35037 * 124.87810 * 179.97438 * 2 1 3 26 26 C 1.50699 * 126.38862 * 0.02562 * 25 2 1 27 27 O 1.42892 * 109.47249 * 330.00334 * 26 25 2 28 28 Si 1.86298 * 109.47032 * 209.99720 * 26 25 2 29 29 H 1.09001 * 109.47477 * 270.05107 * 1 2 3 30 30 H 1.09008 * 109.46530 * 30.04399 * 1 2 3 31 31 H 1.09002 * 109.47624 * 150.04098 * 1 2 3 32 32 H 0.97003 * 119.99465 * 0.02562 * 7 5 4 33 33 H 1.09005 * 109.46999 * 299.99988 * 8 7 5 34 34 H 1.08997 * 109.47332 * 59.99885 * 8 7 5 35 35 H 1.09004 * 113.61119 * 254.52644 * 11 9 8 36 36 H 1.08999 * 113.61458 * 25.34564 * 11 9 8 37 37 H 1.08995 * 109.47424 * 304.99415 * 14 12 11 38 38 H 1.09001 * 109.47547 * 65.00414 * 14 12 11 39 39 H 0.96999 * 120.00275 * 0.02562 * 15 14 12 40 40 H 1.09001 * 109.47027 * 59.99536 * 20 19 18 41 41 H 1.09000 * 109.47502 * 179.97438 * 20 19 18 42 42 H 1.08998 * 109.47213 * 299.99814 * 20 19 18 43 43 H 1.09004 * 109.46728 * 59.99648 * 21 19 18 44 44 H 1.08993 * 109.47723 * 179.97438 * 21 19 18 45 45 H 1.09006 * 109.47216 * 300.00631 * 21 19 18 46 46 H 1.09002 * 109.47244 * 59.99801 * 22 19 18 47 47 H 1.08996 * 109.47723 * 180.02562 * 22 19 18 48 48 H 1.09004 * 109.46771 * 300.00511 * 22 19 18 49 49 H 1.09003 * 113.61991 * 334.56851 * 23 9 8 50 50 H 1.08999 * 113.61521 * 105.47023 * 23 9 8 51 51 H 1.08001 * 127.25358 * 180.24670 * 24 4 3 52 52 H 1.09004 * 109.46453 * 90.00402 * 26 25 2 53 53 H 0.96703 * 114.00176 * 299.99748 * 27 26 25 54 54 H 1.48506 * 109.99834 * 58.68765 * 28 26 25 55 55 H 1.48502 * 110.00011 * 180.02562 * 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.1998 1.0537 0.0000 4 6 3.4740 0.7168 0.0006 5 6 4.6273 1.6396 0.0002 6 8 5.7597 1.1973 0.0003 7 7 4.4243 2.9720 -0.0002 8 6 5.5681 3.8873 -0.0011 9 6 5.0650 5.3323 -0.0021 10 9 4.0846 5.5668 0.9679 11 6 6.1932 6.3773 -0.0025 12 6 5.3367 7.1863 -0.9907 13 1 4.6470 7.8777 -0.5065 14 6 6.1476 7.8424 -2.1101 15 7 5.2655 8.6822 -2.9243 16 6 5.7645 9.3632 -3.9749 17 8 6.9459 9.2803 -4.2470 18 8 4.9538 10.1346 -4.7236 19 6 5.5596 10.8387 -5.8396 20 6 6.1737 9.8264 -6.8085 21 6 4.4908 11.6566 -6.5673 22 6 6.6525 11.7746 -5.3195 23 6 4.6886 5.8525 -1.3994 24 6 3.5271 -0.6973 0.0066 25 6 2.2371 -1.1078 0.0005 26 6 1.7531 -2.5350 0.0006 27 8 0.4950 -2.6112 0.6737 28 14 2.9940 -3.6114 0.8794 29 1 -0.3634 0.0009 1.0277 30 1 -0.3633 0.8897 -0.5146 31 1 -0.3634 -0.8903 -0.5132 32 1 3.5207 3.3249 0.0000 33 1 6.1732 3.7135 0.8887 34 1 6.1724 3.7128 -0.8913 35 1 6.3322 6.8701 0.9598 36 1 7.1270 6.0116 -0.4296 37 1 6.9354 8.4576 -1.6754 38 1 6.5937 7.0699 -2.7365 39 1 4.3226 8.7487 -2.7069 40 1 5.3951 9.1596 -7.1790 41 1 6.6281 10.3549 -7.6465 42 1 6.9351 9.2433 -6.2904 43 1 4.0528 12.3778 -5.8772 44 1 4.9455 12.1856 -7.4048 45 1 3.7124 10.9895 -6.9379 46 1 7.4139 11.1916 -4.8013 47 1 7.1069 12.3035 -6.1572 48 1 6.2149 12.4958 -4.6291 49 1 5.2105 5.3438 -2.2100 50 1 3.6135 5.9217 -1.5646 51 1 4.4107 -1.3183 0.0123 52 1 1.6369 -2.8792 -1.0271 53 1 0.5236 -2.3196 1.5952 54 1 3.1852 -3.1329 2.2722 55 1 2.5286 -5.0216 0.8858 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_16498590_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:30:07 Heat of formation + Delta-G solvation = -91.392277 kcal Electronic energy + Delta-G solvation = -37479.507284 eV Core-core repulsion = 32242.567600 eV Total energy + Delta-G solvation = -5236.939685 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.87 seconds Orbital eigenvalues (eV) -48.84357 -43.85841 -41.81896 -40.63339 -39.80112 -38.24863 -38.17009 -36.11011 -35.48404 -34.71499 -34.40702 -32.68566 -30.46145 -28.68317 -28.37630 -28.31086 -27.49989 -27.00095 -26.58163 -24.24459 -23.35433 -22.28143 -22.20732 -21.35283 -20.97026 -20.47634 -19.36287 -19.10457 -17.82237 -17.68636 -17.60210 -17.48237 -17.32624 -17.12604 -16.89888 -16.60150 -16.32952 -16.04447 -15.97180 -15.76686 -15.66192 -15.52920 -15.43531 -15.07725 -14.88665 -14.68909 -14.46058 -14.20283 -14.08422 -13.93000 -13.68283 -13.58295 -13.41424 -13.35466 -13.32015 -13.12305 -12.92365 -12.87581 -12.69120 -12.54041 -12.37035 -12.29556 -12.19455 -12.05738 -12.00543 -11.60309 -11.40219 -11.35423 -11.22918 -11.10533 -11.08154 -10.38962 -10.21992 -10.07432 -9.98626 -4.89441 0.41029 0.80983 1.16283 1.21590 1.52854 1.67081 1.74994 1.92636 2.02452 2.26844 3.02756 3.08286 3.12306 3.17637 3.28000 3.30845 3.34475 3.49538 3.57021 3.73598 3.78688 3.93645 3.97255 4.05473 4.23989 4.28303 4.34624 4.36009 4.43301 4.47836 4.50682 4.58742 4.60394 4.62035 4.64084 4.66874 4.68892 4.73051 4.74918 4.75569 4.82730 4.86674 4.88787 4.91486 5.00045 5.07266 5.18949 5.21559 5.26218 5.26687 5.33420 5.41421 5.58472 5.65567 6.25131 6.29164 6.37712 6.77840 7.03855 7.32780 Net atomic charges, atomic populations, and dipole contributions calculated with CM2 Atom NO. Type Charge No. of electrons 1 C 0.077 3.923 2 N -0.264 5.264 3 N -0.267 5.267 4 C -0.042 4.042 5 C 0.605 3.395 6 O -0.561 6.561 7 N -0.709 5.709 8 C 0.110 3.890 9 C 0.079 3.921 10 F -0.197 7.197 11 C -0.132 4.132 12 C -0.110 4.110 13 H 0.111 0.889 14 C 0.150 3.850 15 N -0.700 5.700 16 C 0.641 3.359 17 O -0.595 6.595 18 O -0.371 6.371 19 C 0.127 3.873 20 C -0.180 4.180 21 C -0.138 4.138 22 C -0.180 4.180 23 C -0.144 4.144 24 C -0.185 4.185 25 C -0.013 4.013 26 C -0.044 4.044 27 O -0.537 6.537 28 Si 0.783 3.217 29 H 0.100 0.900 30 H 0.101 0.899 31 H 0.118 0.882 32 H 0.414 0.586 33 H 0.077 0.923 34 H 0.101 0.899 35 H 0.093 0.907 36 H 0.123 0.877 37 H 0.071 0.929 38 H 0.074 0.926 39 H 0.416 0.584 40 H 0.061 0.939 41 H 0.086 0.914 42 H 0.064 0.936 43 H 0.069 0.931 44 H 0.098 0.902 45 H 0.069 0.931 46 H 0.064 0.936 47 H 0.086 0.914 48 H 0.061 0.939 49 H 0.113 0.887 50 H 0.088 0.912 51 H 0.183 0.817 52 H 0.147 0.853 53 H 0.388 0.612 54 H -0.296 1.296 55 H -0.283 1.283 For the charges calculated by CM2: Dipole moment (debyes) X Y Z Total from point charges -4.640 0.269 -1.536 4.895 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.061 4.061 2 N -0.053 5.053 3 N -0.105 5.105 4 C -0.200 4.200 5 C 0.391 3.609 6 O -0.438 6.438 7 N -0.361 5.361 8 C -0.014 4.014 9 C 0.060 3.940 10 F -0.178 7.178 11 C -0.170 4.170 12 C -0.129 4.129 13 H 0.129 0.871 14 C 0.028 3.972 15 N -0.360 5.360 16 C 0.395 3.605 17 O -0.486 6.486 18 O -0.285 6.285 19 C 0.091 3.909 20 C -0.237 4.237 21 C -0.195 4.195 22 C -0.237 4.237 23 C -0.182 4.182 24 C -0.213 4.213 25 C -0.138 4.138 26 C -0.154 4.154 27 O -0.345 6.345 28 Si 0.686 3.314 29 H 0.118 0.882 30 H 0.120 0.880 31 H 0.136 0.864 32 H 0.250 0.750 33 H 0.095 0.905 34 H 0.119 0.881 35 H 0.112 0.888 36 H 0.141 0.859 37 H 0.089 0.911 38 H 0.092 0.908 39 H 0.253 0.747 40 H 0.080 0.920 41 H 0.105 0.895 42 H 0.083 0.917 43 H 0.088 0.912 44 H 0.117 0.883 45 H 0.088 0.912 46 H 0.083 0.917 47 H 0.105 0.895 48 H 0.080 0.920 49 H 0.131 0.869 50 H 0.107 0.893 51 H 0.200 0.800 52 H 0.164 0.836 53 H 0.238 0.762 54 H -0.222 1.222 55 H -0.210 1.210 Dipole moment (debyes) X Y Z Total from point charges -4.181 0.169 -1.619 4.487 hybrid contribution -1.177 0.120 1.068 1.594 sum -5.358 0.289 -0.550 5.394 Atomic orbital electron populations 1.22616 0.70756 1.06755 1.06019 1.45697 1.08994 1.01907 1.48748 1.73476 1.02814 1.14818 1.19361 1.21795 0.90847 0.94554 1.12766 1.16989 0.86478 0.81816 0.75630 1.90796 1.23611 1.77655 1.51752 1.45648 1.11205 1.04947 1.74254 1.21965 0.88722 0.86068 1.04599 1.22130 0.89151 0.94091 0.88589 1.92861 1.64293 1.95489 1.65146 1.23563 0.97671 0.94484 1.01327 1.22706 0.98466 0.96110 0.95629 0.87069 1.21380 0.91934 0.93084 0.90789 1.44733 1.14239 1.48381 1.28662 1.18368 0.83411 0.78267 0.80481 1.90911 1.16649 1.68455 1.72559 1.86308 1.41933 1.58041 1.42231 1.22342 0.94327 0.91284 0.82955 1.22474 1.01852 0.99739 0.99652 1.21893 0.97719 0.99484 1.00405 1.22474 0.98934 0.99399 1.02926 1.23796 1.03532 0.95782 0.95081 1.21729 0.97489 0.93370 1.08714 1.22206 0.90993 0.92100 1.08536 1.24855 0.93403 0.91984 1.05195 1.86579 1.29305 1.88553 1.30056 0.93242 0.83910 0.74828 0.79458 0.88206 0.88039 0.86415 0.75014 0.90489 0.88096 0.88810 0.85858 0.91075 0.90764 0.74702 0.91952 0.89505 0.91730 0.91209 0.88339 0.91186 0.91746 0.89490 0.91966 0.86918 0.89350 0.79956 0.83577 0.76212 1.22228 1.20966 Number of geometries 1 Number of calculations of the screened coulomb radii 1 The total number of SCF iterations 14. 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.50 10.12 127.77 1.29 0.79 16 2 N -0.26 -0.14 3.70 -355.91 -1.32 -1.46 16 3 N -0.27 -1.27 11.67 -48.18 -0.56 -1.83 16 4 C -0.04 -0.28 7.13 42.46 0.30 0.02 16 5 C 0.61 4.96 7.82 86.74 0.68 5.64 16 6 O -0.56 -7.67 16.93 -3.94 -0.07 -7.74 16 7 N -0.71 -2.03 5.46 -466.11 -2.55 -4.57 16 8 C 0.11 -0.08 5.37 86.38 0.46 0.38 16 9 C 0.08 -0.30 1.67 -10.27 -0.02 -0.32 16 10 F -0.20 -0.26 16.80 44.97 0.76 0.49 16 11 C -0.13 0.93 7.10 31.12 0.22 1.15 16 12 C -0.11 0.79 3.73 -10.25 -0.04 0.75 16 13 H 0.11 -0.83 8.14 -2.39 -0.02 -0.85 16 14 C 0.15 -0.55 5.51 86.38 0.48 -0.07 16 15 N -0.70 0.65 5.56 -590.84 -3.29 -2.64 16 16 C 0.64 3.88 8.06 129.79 1.05 4.93 16 17 O -0.59 -6.82 11.88 19.81 0.24 -6.58 16 18 O -0.37 -1.93 9.94 -55.48 -0.55 -2.48 16 19 C 0.13 0.36 1.13 -10.79 -0.01 0.35 16 20 C -0.18 -0.56 8.37 71.98 0.60 0.04 16 21 C -0.14 0.29 8.85 71.98 0.64 0.93 16 22 C -0.18 -0.57 8.37 71.98 0.60 0.03 16 23 C -0.14 0.95 7.08 31.13 0.22 1.17 16 24 C -0.19 -0.58 9.35 22.41 0.21 -0.37 16 25 C -0.01 0.00 6.18 41.40 0.26 0.26 16 26 C -0.04 0.21 2.87 71.24 0.20 0.41 16 27 O -0.54 2.55 10.40 -148.98 -1.55 1.00 16 28 Si 0.78 -2.92 31.22 68.60 2.14 -0.78 16 29 H 0.10 -0.93 8.14 -2.39 -0.02 -0.94 16 30 H 0.10 -0.78 8.14 -2.38 -0.02 -0.80 16 31 H 0.12 -0.96 6.77 -2.39 -0.02 -0.97 16 32 H 0.41 0.57 7.91 -92.71 -0.73 -0.16 16 33 H 0.08 0.07 8.14 -2.38 -0.02 0.05 16 34 H 0.10 -0.23 7.94 -2.39 -0.02 -0.25 16 35 H 0.09 -0.63 8.14 -2.39 -0.02 -0.65 16 36 H 0.12 -1.17 8.06 -2.39 -0.02 -1.19 16 37 H 0.07 -0.14 8.14 -2.39 -0.02 -0.16 16 38 H 0.07 -0.20 7.85 -2.39 -0.02 -0.22 16 39 H 0.42 -2.10 8.58 -92.71 -0.80 -2.90 16 40 H 0.06 0.15 8.14 -2.39 -0.02 0.13 16 41 H 0.09 -0.02 8.14 -2.39 -0.02 -0.04 16 42 H 0.06 0.54 5.88 -2.39 -0.01 0.52 16 43 H 0.07 -0.16 8.14 -2.38 -0.02 -0.18 16 44 H 0.10 -0.53 8.14 -2.39 -0.02 -0.55 16 45 H 0.07 -0.16 8.14 -2.38 -0.02 -0.18 16 46 H 0.06 0.54 5.88 -2.39 -0.01 0.53 16 47 H 0.09 -0.02 8.14 -2.39 -0.02 -0.04 16 48 H 0.06 0.15 8.14 -2.38 -0.02 0.14 16 49 H 0.11 -0.84 7.72 -2.39 -0.02 -0.86 16 50 H 0.09 -0.57 8.14 -2.39 -0.02 -0.59 16 51 H 0.18 -0.07 7.01 -2.91 -0.02 -0.09 16 52 H 0.15 -1.13 7.96 -2.38 -0.02 -1.15 16 53 H 0.39 -3.48 8.70 -74.06 -0.64 -4.13 16 54 H -0.30 -0.88 7.11 99.48 0.71 -0.17 16 55 H -0.28 -0.57 7.11 99.48 0.71 0.14 16 Total: 0.00 -25.27 450.74 -0.79 -26.06 By element: Atomic # 1 Polarization: -14.39 SS G_CDS: -1.19 Total: -15.58 kcal Atomic # 6 Polarization: 8.96 SS G_CDS: 7.15 Total: 16.11 kcal Atomic # 7 Polarization: -2.79 SS G_CDS: -7.71 Total: -10.50 kcal Atomic # 8 Polarization: -13.87 SS G_CDS: -1.93 Total: -15.80 kcal Atomic # 9 Polarization: -0.26 SS G_CDS: 0.76 Total: 0.49 kcal Atomic # 14 Polarization: -2.92 SS G_CDS: 2.14 Total: -0.78 kcal Total: -25.27 -0.79 -26.06 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_16498590_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 -65.335 kcal (2) G-P(sol) polarization free energy of solvation -25.269 kcal (3) G-ENP(sol) elect.-nuc.-pol. free energy of system -90.605 kcal (4) G-CDS(sol) cavity-dispersion-solvent structure free energy -0.788 kcal (5) G-P-CDS(sol) = G-P(sol) + G-CDS(sol) = (2) + (4) -26.057 kcal (6) G-S(sol) free energy of system = (1) + (5) -91.392 kcal 1SCF WAS SPECIFIED, SO GEOMETRY OPTIMIZATION WAS NOT USED Total computer time = 1.88 seconds