Wall clock time and date at job start Tue Mar 31 2020 23:55:15 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 O 1.42896 * 1 3 3 C 1.42891 * 114.00139 * 2 1 4 4 C 1.53004 * 109.47174 * 180.02562 * 3 2 1 5 5 C 1.53000 * 109.45961 * 299.97387 * 4 3 2 6 6 N 1.46491 * 109.47263 * 60.02592 * 5 4 3 7 7 C 1.34779 * 120.00362 * 179.97438 * 6 5 4 8 8 O 1.21281 * 119.99793 * 0.02562 * 7 6 5 9 9 C 1.50701 * 120.00190 * 179.97438 * 7 6 5 10 10 H 1.09001 * 110.98530 * 55.80625 * 9 7 6 11 11 C 1.54943 * 111.12625 * 291.55870 * 9 7 6 12 12 O 1.43642 * 105.13575 * 205.63789 * 11 9 7 13 13 C 1.43239 * 109.97252 * 24.42737 * 12 11 9 14 14 H 1.09000 * 109.82980 * 238.77711 * 13 12 11 15 15 C 1.53004 * 109.83038 * 117.80269 * 13 12 11 16 16 C 1.54385 * 107.58682 * 358.28883 * 13 12 11 17 17 C 1.53050 * 109.45804 * 180.02562 * 4 3 2 18 18 C 1.53193 * 109.30912 * 178.61474 * 17 4 3 19 19 N 1.46927 * 108.77611 * 54.63249 * 18 17 4 20 20 C 1.36934 * 120.62780 * 126.36533 * 19 18 17 21 21 H 1.08000 * 120.00322 * 332.25176 * 20 19 18 22 22 C 1.38814 * 119.99759 * 152.25025 * 20 19 18 23 23 N 1.31620 * 119.99729 * 143.92630 * 22 20 19 24 24 Si 1.86802 * 120.00099 * 323.92471 * 22 20 19 25 25 H 1.48505 * 109.46754 * 270.00125 * 24 22 20 26 26 H 1.48500 * 109.47020 * 150.00301 * 24 22 20 27 27 H 1.48496 * 109.46992 * 29.99793 * 24 22 20 28 28 C 1.46922 * 118.73907 * 306.64224 * 19 18 17 29 29 C 1.52760 * 109.00969 * 53.31149 * 28 19 18 30 30 H 1.08994 * 109.47501 * 59.99922 * 1 2 3 31 31 H 1.09003 * 109.47190 * 179.97438 * 1 2 3 32 32 H 1.09003 * 109.47221 * 299.99620 * 1 2 3 33 33 H 1.09003 * 109.47229 * 60.00132 * 3 2 1 34 34 H 1.08997 * 109.47376 * 300.00213 * 3 2 1 35 35 H 1.09000 * 109.47042 * 180.02562 * 5 4 3 36 36 H 1.09006 * 109.47140 * 300.02495 * 5 4 3 37 37 H 0.97000 * 119.99969 * 359.97438 * 6 5 4 38 38 H 1.09002 * 110.30984 * 324.55168 * 11 9 7 39 39 H 1.08995 * 110.30524 * 86.72245 * 11 9 7 40 40 H 1.09004 * 109.46950 * 301.86852 * 15 13 12 41 41 H 1.09004 * 109.46895 * 61.86152 * 15 13 12 42 42 H 1.08995 * 109.47454 * 181.86768 * 15 13 12 43 43 H 1.08998 * 110.84233 * 220.44827 * 16 13 12 44 44 H 1.08999 * 110.84269 * 96.83348 * 16 13 12 45 45 H 1.09002 * 109.49708 * 58.66801 * 17 4 3 46 46 H 1.08998 * 109.49389 * 298.56415 * 17 4 3 47 47 H 1.09003 * 109.58966 * 294.84963 * 18 17 4 48 48 H 1.08999 * 109.70198 * 174.49211 * 18 17 4 49 49 H 0.96998 * 120.00804 * 354.92389 * 23 22 20 50 50 H 1.09007 * 109.58776 * 173.13990 * 28 19 18 51 51 H 1.08995 * 109.59076 * 293.56539 * 28 19 18 52 52 H 1.09000 * 109.56762 * 185.55184 * 29 28 19 53 53 H 1.09003 * 109.36309 * 65.49758 * 29 28 19 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 8 1.4290 0.0000 0.0000 3 6 2.0102 1.3054 0.0000 4 6 3.5354 1.1845 -0.0006 5 6 3.9871 0.4256 1.2487 6 7 3.3917 -0.9129 1.2490 7 6 3.6514 -1.7659 2.2596 8 8 4.3777 -1.4252 3.1692 9 6 3.0383 -3.1425 2.2602 10 1 3.2835 -3.6794 1.3438 11 6 1.5053 -3.0711 2.4734 12 8 1.1478 -4.3361 3.0524 13 6 2.2805 -4.9073 3.7177 14 1 2.5110 -5.8809 3.2853 15 6 1.9839 -5.0577 5.2111 16 6 3.4733 -3.9484 3.5145 17 6 4.1577 2.5828 0.0000 18 6 5.6841 2.4580 0.0343 19 7 6.1142 1.6066 -1.0832 20 6 7.0593 2.0459 -1.9714 21 1 7.1750 3.1034 -2.1579 22 6 7.8687 1.1321 -2.6321 23 7 8.2272 1.3518 -3.8793 24 14 8.4492 -0.4095 -1.7511 25 1 9.7400 -0.1342 -1.0704 26 1 8.6335 -1.5038 -2.7379 27 1 7.4357 -0.8158 -0.7447 28 6 5.5066 0.2774 -1.2331 29 6 3.9862 0.4243 -1.2501 30 1 -0.3634 0.5138 -0.8899 31 1 -0.3634 -1.0277 -0.0005 32 1 -0.3634 0.5138 0.8900 33 1 1.6886 1.8463 0.8900 34 1 1.6886 1.8463 -0.8899 35 1 5.0738 0.3398 1.2486 36 1 3.6654 0.9666 2.1387 37 1 2.8104 -1.1852 0.5217 38 1 1.2557 -2.2586 3.1559 39 1 0.9966 -2.9345 1.5192 40 1 1.7382 -4.0829 5.6327 41 1 1.1412 -5.7354 5.3477 42 1 2.8609 -5.4616 5.7169 43 1 4.3897 -4.5053 3.3193 44 1 3.5934 -3.2909 4.3755 45 1 3.8177 3.1319 0.8781 46 1 3.8563 3.1155 -0.9020 47 1 5.9928 2.0067 0.9773 48 1 6.1358 3.4452 -0.0635 49 1 7.8672 2.1131 -4.3606 50 1 5.8391 -0.1731 -2.1684 51 1 5.8004 -0.3554 -0.3957 52 1 3.5241 -0.5628 -1.2617 53 1 3.6864 0.9773 -2.1403 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 ZINC001162335888.mol2 53 1SCF WAS SPECIFIED, SO GEOMETRY OPTIMIZATION WAS NOT USED SCF FIELD WAS ACHIEVED Wall clock time and date at job start Tue Mar 31 2020 23:55:15 Heat of formation + Delta-G solvation = -106.190191 kcal Electronic energy + Delta-G solvation = -32496.689378 eV Core-core repulsion = 28369.122367 eV Total energy + Delta-G solvation = -4127.567011 eV No. of doubly occupied orbitals = 66 Molecular weight (most abundant/longest-lived isotopes) = 340.220 amu Computer time = 1.46 seconds Orbital eigenvalues (eV) -40.58996 -39.67623 -38.28763 -37.29772 -34.68247 -33.22369 -32.53686 -31.86953 -31.19201 -30.41475 -27.81858 -26.81209 -26.18198 -25.68001 -23.41876 -23.12320 -22.69811 -21.28169 -20.87590 -19.99290 -18.65184 -17.63044 -17.26625 -16.69718 -16.32188 -15.99269 -15.93458 -15.59906 -15.25256 -14.93090 -14.40516 -14.35653 -14.20315 -13.98239 -13.93087 -13.53203 -13.38526 -12.98275 -12.94102 -12.84592 -12.81361 -12.42216 -12.16090 -12.05538 -11.89359 -11.66533 -11.49950 -11.39942 -11.31922 -11.14773 -11.03705 -10.70020 -10.51044 -10.36683 -10.10967 -9.86693 -9.80980 -9.68578 -9.64856 -9.60403 -9.21446 -8.75661 -8.26230 -6.97650 -6.35911 -3.24351 2.70838 3.16368 3.52117 3.53325 3.60233 3.73582 3.76350 4.15692 4.57431 4.65350 4.66489 4.78103 4.85609 4.88941 4.92236 4.96128 5.00505 5.12620 5.15990 5.17094 5.26720 5.36968 5.39234 5.41143 5.44203 5.46274 5.49095 5.51071 5.55621 5.63594 5.76669 5.83364 5.85238 5.86747 6.02059 6.03385 6.10773 6.18866 6.34966 6.41682 6.46410 6.59608 6.67153 6.73780 6.75330 6.80255 6.98153 7.03428 7.78401 8.06930 8.40743 8.46506 9.03457 9.65826 10.25161 11.04549 Molecular weight = 340.22amu Principal moments of inertia in cm(-1) A = 0.012516 B = 0.003654 C = 0.003123 Principal moments of inertia in units of 10**(-40)*gram-cm**2 A = 2236.560631 B = 7661.196111 C = 8962.879059 Net atomic charges, atomic populations, and dipole contributions calculated with CM2 Atom NO. Type Charge No. of electrons 1 C 0.033 3.967 2 O -0.393 6.393 3 C 0.073 3.927 4 C -0.059 4.059 5 C 0.130 3.870 6 N -0.706 5.706 7 C 0.519 3.481 8 O -0.546 6.546 9 C -0.148 4.148 10 H 0.111 0.889 11 C 0.070 3.930 12 O -0.378 6.378 13 C 0.074 3.926 14 H 0.082 0.918 15 C -0.162 4.162 16 C -0.132 4.132 17 C -0.119 4.119 18 C 0.142 3.858 19 N -0.598 5.598 20 C -0.361 4.361 21 H 0.064 0.936 22 C 0.005 3.995 23 N -0.899 5.899 24 Si 1.013 2.987 25 H -0.298 1.298 26 H -0.295 1.295 27 H -0.269 1.269 28 C 0.142 3.858 29 C -0.122 4.122 30 H 0.041 0.959 31 H 0.084 0.916 32 H 0.039 0.961 33 H 0.048 0.952 34 H 0.052 0.948 35 H 0.081 0.919 36 H 0.059 0.941 37 H 0.410 0.590 38 H 0.063 0.937 39 H 0.092 0.908 40 H 0.060 0.940 41 H 0.061 0.939 42 H 0.067 0.933 43 H 0.086 0.914 44 H 0.093 0.907 45 H 0.057 0.943 46 H 0.068 0.932 47 H 0.023 0.977 48 H 0.060 0.940 49 H 0.267 0.733 50 H 0.064 0.936 51 H 0.032 0.968 52 H 0.058 0.942 53 H 0.063 0.937 For the charges calculated by CM2: Dipole moment (debyes) X Y Z Total from point charges -18.154 -7.601 12.307 23.212 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.064 4.064 2 O -0.312 6.312 3 C -0.004 4.004 4 C -0.060 4.060 5 C 0.006 3.994 6 N -0.363 5.363 7 C 0.303 3.697 8 O -0.424 6.424 9 C -0.170 4.170 10 H 0.129 0.871 11 C -0.007 4.007 12 O -0.297 6.297 13 C 0.017 3.983 14 H 0.100 0.900 15 C -0.219 4.219 16 C -0.170 4.170 17 C -0.159 4.159 18 C 0.016 3.984 19 N -0.322 5.322 20 C -0.489 4.489 21 H 0.081 0.919 22 C -0.199 4.199 23 N -0.536 5.536 24 Si 0.842 3.158 25 H -0.226 1.226 26 H -0.221 1.221 27 H -0.195 1.195 28 C 0.017 3.983 29 C -0.161 4.161 30 H 0.060 0.940 31 H 0.102 0.898 32 H 0.057 0.943 33 H 0.066 0.934 34 H 0.071 0.929 35 H 0.099 0.901 36 H 0.077 0.923 37 H 0.251 0.749 38 H 0.081 0.919 39 H 0.110 0.890 40 H 0.079 0.921 41 H 0.080 0.920 42 H 0.086 0.914 43 H 0.105 0.895 44 H 0.112 0.888 45 H 0.075 0.925 46 H 0.086 0.914 47 H 0.041 0.959 48 H 0.079 0.921 49 H 0.076 0.924 50 H 0.082 0.918 51 H 0.049 0.951 52 H 0.077 0.923 53 H 0.081 0.919 Dipole moment (debyes) X Y Z Total from point charges -16.758 -8.907 12.729 22.852 hybrid contribution -1.172 1.658 -0.246 2.045 sum -17.930 -7.249 12.484 23.019 Atomic orbital electron populations 1.22719 0.81155 1.02559 0.99950 1.87530 1.19485 1.28779 1.95364 1.22735 0.93535 0.84008 1.00084 1.20800 0.91844 0.97271 0.96042 1.21473 0.99257 0.82023 0.96634 1.46276 1.50009 1.14686 1.25297 1.20056 0.81309 0.85983 0.82303 1.90690 1.37767 1.77656 1.36310 1.23014 0.97097 0.95184 1.01660 0.87113 1.23225 0.91951 0.86790 0.98736 1.88425 1.34741 1.34287 1.72246 1.22786 0.85266 0.96974 0.93270 0.90027 1.22136 1.02170 1.02624 0.94998 1.23091 0.98465 0.97614 0.97819 1.21971 0.99064 0.93505 1.01355 1.20978 0.90799 0.95728 0.90894 1.44304 1.44601 1.14601 1.28658 1.20022 1.15019 0.94100 1.19713 0.91866 1.26086 0.95672 1.02253 0.95873 1.70889 1.46839 1.28534 1.07336 0.84684 0.77438 0.77105 0.76567 1.22563 1.22119 1.19540 1.21314 0.89686 0.87773 0.99570 1.22002 0.99300 0.98982 0.95831 0.94045 0.89777 0.94265 0.93402 0.92938 0.90056 0.92290 0.74854 0.91897 0.88968 0.92126 0.91991 0.91369 0.89532 0.88849 0.92453 0.91360 0.95887 0.92142 0.92353 0.91759 0.95083 0.92318 0.91851 Number of geometries 1 Number of calculations of the screened coulomb radii 1 The total number of SCF iterations 16. 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.03 0.26 11.01 101.05 1.11 1.37 16 2 O -0.39 -4.18 6.73 -35.23 -0.24 -4.42 16 3 C 0.07 0.84 4.80 37.16 0.18 1.02 16 4 C -0.06 -0.84 0.51 -154.46 -0.08 -0.92 16 5 C 0.13 1.96 3.80 -4.05 -0.02 1.95 16 6 N -0.71 -9.64 4.63 -60.30 -0.28 -9.92 16 7 C 0.52 7.43 7.15 -10.98 -0.08 7.35 16 8 O -0.55 -9.81 15.71 5.56 0.09 -9.73 16 9 C -0.15 -1.63 3.49 -89.53 -0.31 -1.94 16 10 H 0.11 1.08 8.14 -51.93 -0.42 0.66 16 11 C 0.07 0.70 7.25 38.18 0.28 0.97 16 12 O -0.38 -4.15 11.02 -35.23 -0.39 -4.54 16 13 C 0.07 0.70 4.16 -26.00 -0.11 0.59 16 14 H 0.08 0.67 8.14 -51.93 -0.42 0.25 16 15 C -0.16 -1.33 9.44 37.16 0.35 -0.98 16 16 C -0.13 -1.38 5.70 -24.78 -0.14 -1.52 16 17 C -0.12 -1.81 4.91 -26.61 -0.13 -1.94 16 18 C 0.14 2.73 6.05 -3.71 -0.02 2.71 16 19 N -0.60 -13.85 2.72 -172.58 -0.47 -14.32 16 20 C -0.36 -9.91 9.70 -15.06 -0.15 -10.06 16 21 H 0.06 1.86 8.06 -52.49 -0.42 1.44 16 22 C 0.01 0.15 4.82 -17.43 -0.08 0.07 16 23 N -0.90 -27.34 13.70 55.49 0.76 -26.58 16 24 Si 1.01 25.50 27.01 -169.99 -4.59 20.90 16 25 H -0.30 -7.55 7.11 56.52 0.40 -7.15 16 26 H -0.29 -7.48 7.11 56.52 0.40 -7.08 16 27 H -0.27 -6.33 4.73 56.52 0.27 -6.06 16 28 C 0.14 3.01 3.88 -3.94 -0.02 2.99 16 29 C -0.12 -2.03 4.57 -26.83 -0.12 -2.16 16 30 H 0.04 0.30 8.14 -51.93 -0.42 -0.12 16 31 H 0.08 0.53 8.14 -51.93 -0.42 0.11 16 32 H 0.04 0.28 8.14 -51.93 -0.42 -0.14 16 33 H 0.05 0.51 8.14 -51.93 -0.42 0.09 16 34 H 0.05 0.59 8.14 -51.93 -0.42 0.17 16 35 H 0.08 1.49 4.99 -51.93 -0.26 1.23 16 36 H 0.06 0.88 8.14 -51.93 -0.42 0.46 16 37 H 0.41 5.05 4.71 -40.82 -0.19 4.86 16 38 H 0.06 0.65 8.02 -51.93 -0.42 0.24 16 39 H 0.09 0.77 8.14 -51.93 -0.42 0.35 16 40 H 0.06 0.53 8.14 -51.93 -0.42 0.11 16 41 H 0.06 0.48 8.14 -51.93 -0.42 0.05 16 42 H 0.07 0.47 8.14 -51.93 -0.42 0.05 16 43 H 0.09 0.84 8.14 -51.93 -0.42 0.41 16 44 H 0.09 1.11 7.06 -51.93 -0.37 0.74 16 45 H 0.06 0.73 8.14 -51.93 -0.42 0.31 16 46 H 0.07 1.08 8.14 -51.93 -0.42 0.66 16 47 H 0.02 0.44 6.69 -51.93 -0.35 0.10 16 48 H 0.06 1.18 7.93 -51.93 -0.41 0.77 16 49 H 0.27 8.02 8.73 -40.82 -0.36 7.66 16 50 H 0.06 1.50 5.83 -51.93 -0.30 1.20 16 51 H 0.03 0.68 5.16 -51.93 -0.27 0.41 16 52 H 0.06 0.89 6.50 -51.93 -0.34 0.55 16 53 H 0.06 1.12 8.14 -51.93 -0.42 0.70 16 LS Contribution 395.63 15.07 5.96 5.96 Total: -1.00 -32.28 395.63 -7.87 -40.14 By element: Atomic # 1 Polarization: 12.38 SS G_CDS: -9.38 Total: 3.00 kcal Atomic # 6 Polarization: -1.17 SS G_CDS: 0.66 Total: -0.50 kcal Atomic # 7 Polarization: -50.83 SS G_CDS: 0.01 Total: -50.82 kcal Atomic # 8 Polarization: -18.15 SS G_CDS: -0.54 Total: -18.69 kcal Atomic # 14 Polarization: 25.50 SS G_CDS: -4.59 Total: 20.90 kcal Total LS contribution 5.96 Total: 5.96 kcal Total: -32.28 -7.87 -40.14 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. ZINC001162335888.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 -66.048 kcal (2) G-P(sol) polarization free energy of solvation -32.275 kcal (3) G-ENP(sol) elect.-nuc.-pol. free energy of system -98.324 kcal (4) G-CDS(sol) cavity-dispersion-solvent structure free energy -7.867 kcal (5) G-P-CDS(sol) = G-P(sol) + G-CDS(sol) = (2) + (4) -40.142 kcal (6) G-S(sol) free energy of system = (1) + (5) -106.190 kcal 1SCF WAS SPECIFIED, SO GEOMETRY OPTIMIZATION WAS NOT USED Total computer time = 1.46 seconds