Wall clock time and date at job start Tue Mar 31 2020 05:50:45 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 * 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 N 2 2 C 1.31613 * 1 3 3 Si 1.86805 * 120.00066 * 2 1 4 4 H 1.48498 * 109.46896 * 274.97848 * 3 2 1 5 5 H 1.48499 * 109.46870 * 34.97769 * 3 2 1 6 6 H 1.48500 * 109.47003 * 154.98028 * 3 2 1 7 7 C 1.38815 * 120.00011 * 180.02562 * 2 1 3 8 8 H 1.08001 * 119.99974 * 162.58214 * 7 2 1 9 9 N 1.36934 * 120.00057 * 342.57912 * 7 2 1 10 10 C 1.47223 * 121.39476 * 332.52776 * 9 7 2 11 11 C 1.53106 * 111.02630 * 115.45301 * 10 9 7 12 12 N 1.46985 * 109.55638 * 43.83950 * 11 10 9 13 13 C 1.36586 * 129.16596 * 165.02327 * 12 11 10 14 14 C 1.39271 * 131.77244 * 356.57507 * 13 12 11 15 15 C 1.38337 * 117.21152 * 178.93510 * 14 13 12 16 16 C 1.38532 * 122.14763 * 0.49639 * 15 14 13 17 17 C 1.50705 * 118.94380 * 179.97438 * 16 15 14 18 18 C 1.53003 * 109.47130 * 299.72480 * 17 16 15 19 19 C 1.53002 * 109.46886 * 180.02562 * 18 17 16 20 20 C 1.53001 * 109.46946 * 300.00275 * 19 18 17 21 21 C 1.53000 * 109.47205 * 60.00054 * 20 19 18 22 22 C 1.52997 * 109.47315 * 299.99556 * 21 20 19 23 23 C 1.37175 * 122.10084 * 359.65428 * 16 15 14 24 24 C 1.39909 * 117.21921 * 0.25995 * 23 16 15 25 25 C 1.43678 * 133.10841 * 179.39878 * 24 23 16 26 26 C 1.34869 * 106.38125 * 180.81886 * 25 24 23 27 27 C 1.47693 * 121.39409 * 152.54357 * 9 7 2 28 28 H 1.09002 * 110.74300 * 106.77106 * 27 9 7 29 29 C 1.51852 * 111.35197 * 343.13596 * 27 9 7 30 30 C 1.51105 * 115.76242 * 205.06756 * 29 27 9 31 31 C 1.51271 * 135.13818 * 2.63178 * 25 24 23 32 32 H 0.97001 * 120.00100 * 359.97438 * 1 2 3 33 33 H 1.09000 * 109.16553 * 355.05182 * 10 9 7 34 34 H 1.09006 * 109.16571 * 235.85559 * 10 9 7 35 35 H 1.08998 * 109.45820 * 283.81396 * 11 10 9 36 36 H 1.08996 * 109.45889 * 163.86446 * 11 10 9 37 37 H 1.08002 * 121.39476 * 359.19856 * 14 13 12 38 38 H 1.07994 * 118.92742 * 180.24088 * 15 14 13 39 39 H 1.09004 * 109.46509 * 179.72923 * 17 16 15 40 40 H 1.08993 * 109.46585 * 60.00673 * 18 17 16 41 41 H 1.08994 * 109.47108 * 300.00013 * 18 17 16 42 42 H 1.09004 * 109.47001 * 60.00291 * 19 18 17 43 43 H 1.08994 * 109.46695 * 180.02562 * 19 18 17 44 44 H 1.09000 * 109.46699 * 180.02562 * 20 19 18 45 45 H 1.08994 * 109.47368 * 299.99785 * 20 19 18 46 46 H 1.09001 * 109.46897 * 60.00438 * 21 20 19 47 47 H 1.09003 * 109.46965 * 179.97438 * 21 20 19 48 48 H 1.08991 * 109.47220 * 179.97438 * 22 21 20 49 49 H 1.08992 * 109.47352 * 300.00581 * 22 21 20 50 50 H 1.07996 * 121.39069 * 180.23297 * 23 16 15 51 51 H 1.09001 * 108.22384 * 326.69887 * 29 27 9 52 52 H 1.08999 * 108.22654 * 83.43275 * 29 27 9 53 53 H 1.09005 * 107.98151 * 121.73913 * 30 29 27 54 54 H 1.08995 * 107.98270 * 238.07695 * 30 29 27 55 55 H 1.09004 * 109.77494 * 17.19519 * 31 25 24 56 56 H 1.09000 * 108.78454 * 257.37343 * 31 25 24 Note: An asterisk (*) indicates this parameter was or will be optimized. Cartesian coordinates (angstroms) NO. Atom X Y Z 1 7 0.0000 0.0000 0.0000 2 6 1.3161 0.0000 0.0000 3 14 2.2502 1.6178 0.0000 4 1 2.3924 2.1071 -1.3948 5 1 1.5042 2.6200 0.8026 6 1 3.5964 1.4121 0.5921 7 6 2.0102 -1.2022 -0.0005 8 1 3.0531 -1.2235 -0.2805 9 7 1.3727 -2.3610 0.3540 10 6 0.2224 -2.3533 1.2729 11 6 -1.0495 -2.8278 0.5649 12 7 -0.7447 -4.0091 -0.2549 13 6 -1.6008 -4.9092 -0.8228 14 6 -2.9863 -5.0373 -0.7626 15 6 -3.5584 -6.0678 -1.4869 16 6 -2.7996 -6.9510 -2.2375 17 6 -3.4982 -8.0483 -2.9986 18 6 -4.2618 -8.9420 -2.0192 19 6 -4.9705 -10.0563 -2.7919 20 6 -5.9515 -9.4397 -3.7910 21 6 -5.1878 -8.5461 -4.7704 22 6 -4.4792 -7.4317 -3.9978 23 6 -1.4335 -6.8548 -2.3152 24 6 -0.8227 -5.8221 -1.5954 25 6 0.5460 -5.4195 -1.4259 26 6 0.5297 -4.3554 -0.5974 27 6 1.8104 -3.6710 -0.1692 28 1 2.3165 -4.2467 0.6058 29 6 2.7028 -3.5097 -1.3872 30 6 2.7341 -4.6833 -2.3385 31 6 1.8768 -5.8877 -1.9719 32 1 -0.4850 0.8400 0.0004 33 1 0.0684 -1.3404 1.6449 34 1 0.4318 -3.0145 2.1139 35 1 -1.4283 -2.0303 -0.0742 36 1 -1.8037 -3.0864 1.3080 37 1 -3.5893 -4.3579 -0.1785 38 1 -4.6311 -6.1903 -1.4631 39 1 -2.7603 -8.6447 -3.5352 40 1 -4.9996 -8.3455 -1.4827 41 1 -3.5628 -9.3815 -1.3077 42 1 -4.2323 -10.6527 -3.3282 43 1 -5.5148 -10.6927 -2.0942 44 1 -6.4567 -10.2336 -4.3412 45 1 -6.6896 -8.8434 -3.2548 46 1 -4.4498 -9.1425 -5.3067 47 1 -5.8867 -8.1072 -5.4825 48 1 -3.9349 -6.7954 -4.6953 49 1 -5.2173 -6.8354 -3.4614 50 1 -0.8535 -7.5475 -2.9069 51 1 2.3694 -2.6315 -1.9402 52 1 3.7198 -3.3232 -1.0424 53 1 2.4161 -4.3299 -3.3194 54 1 3.7672 -5.0198 -2.4239 55 1 1.7164 -6.5034 -2.8569 56 1 2.3831 -6.4755 -1.2062 RHF calculation, no. of doubly occupied orbitals= 69 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) 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=WATER ZINC000019201937.mol2 56 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 05:50:45 Heat of formation + Delta-G solvation = 46.472611 kcal Electronic energy + Delta-G solvation = -33849.981428 eV Core-core repulsion = 29915.724106 eV Total energy + Delta-G solvation = -3934.257322 eV No. of doubly occupied orbitals = 69 Molecular weight (most abundant/longest-lived isotopes) = 364.221 amu Computer time = 1.80 seconds Orbital eigenvalues (eV) -42.60892 -40.17062 -38.58654 -37.09889 -36.76548 -33.85687 -32.78424 -32.40095 -31.53168 -31.46505 -29.89086 -29.22115 -27.71284 -25.82043 -25.36731 -23.83926 -23.36660 -23.32963 -21.98940 -21.53826 -21.04540 -20.13912 -19.15625 -18.53901 -17.55587 -17.31311 -16.80942 -16.60243 -16.11496 -16.00911 -15.89219 -15.39526 -15.06077 -14.98176 -14.90327 -14.79842 -14.47022 -14.29071 -14.05757 -13.84835 -13.76867 -13.63275 -13.50377 -13.27274 -13.08384 -12.97542 -12.81851 -12.71821 -12.48621 -12.45890 -12.14430 -12.09929 -12.01942 -11.93304 -11.68650 -11.55314 -11.52345 -11.40633 -11.20558 -11.07954 -10.84633 -10.83162 -10.65198 -10.32094 -8.84171 -8.69785 -8.20799 -7.95400 -5.39375 0.27460 0.81366 1.44311 1.47304 1.56049 1.86606 2.05994 2.47675 3.01527 3.22995 3.35308 3.71634 3.77904 3.88243 3.91649 3.99694 4.02998 4.06787 4.16456 4.17361 4.31439 4.35899 4.39372 4.40914 4.46912 4.49249 4.57961 4.63346 4.65380 4.70148 4.71735 4.72308 4.76069 4.83549 4.89275 4.97755 4.99081 5.06314 5.11458 5.16074 5.19457 5.24859 5.28279 5.31165 5.32727 5.34329 5.38282 5.42128 5.43535 5.46877 5.55512 5.56019 5.63340 5.65098 5.71379 5.92845 6.03463 6.16257 6.30485 6.45803 6.60908 7.02538 7.68049 8.25242 9.12470 Molecular weight = 364.22amu Principal moments of inertia in cm(-1) A = 0.016644 B = 0.002496 C = 0.002334 Principal moments of inertia in units of 10**(-40)*gram-cm**2 A = 1681.881878 B =11215.087734 C =11991.268319 Net atomic charges, atomic populations, and dipole contributions calculated with CM2 Atom NO. Type Charge No. of electrons 1 N -0.927 5.927 2 C -0.045 4.045 3 Si 0.966 3.034 4 H -0.292 1.292 5 H -0.293 1.293 6 H -0.269 1.269 7 C -0.251 4.251 8 H 0.074 0.926 9 N -0.605 5.605 10 C 0.162 3.838 11 C 0.111 3.889 12 N -0.462 5.462 13 C 0.060 3.940 14 C -0.142 4.142 15 C -0.099 4.099 16 C -0.109 4.109 17 C -0.033 4.033 18 C -0.108 4.108 19 C -0.116 4.116 20 C -0.118 4.118 21 C -0.116 4.116 22 C -0.110 4.110 23 C -0.076 4.076 24 C -0.104 4.104 25 C -0.207 4.207 26 C 0.047 3.953 27 C 0.225 3.775 28 H 0.061 0.939 29 C -0.100 4.100 30 C -0.103 4.103 31 C -0.025 4.025 32 H 0.250 0.750 33 H 0.086 0.914 34 H 0.022 0.978 35 H 0.066 0.934 36 H 0.078 0.922 37 H 0.123 0.877 38 H 0.134 0.866 39 H 0.078 0.922 40 H 0.060 0.940 41 H 0.056 0.944 42 H 0.065 0.935 43 H 0.067 0.933 44 H 0.069 0.931 45 H 0.059 0.941 46 H 0.064 0.936 47 H 0.065 0.935 48 H 0.051 0.949 49 H 0.059 0.941 50 H 0.141 0.859 51 H 0.041 0.959 52 H 0.069 0.931 53 H 0.048 0.952 54 H 0.089 0.911 55 H 0.085 0.915 56 H 0.077 0.923 For the charges calculated by CM2: Dipole moment (debyes) X Y Z Total from point charges -7.016 -18.967 -5.838 21.049 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 N -0.572 5.572 2 C -0.241 4.241 3 Si 0.795 3.205 4 H -0.219 1.219 5 H -0.219 1.219 6 H -0.194 1.194 7 C -0.381 4.381 8 H 0.092 0.908 9 N -0.328 5.328 10 C 0.035 3.965 11 C -0.012 4.012 12 N -0.141 5.141 13 C -0.047 4.047 14 C -0.162 4.162 15 C -0.118 4.118 16 C -0.110 4.110 17 C -0.052 4.052 18 C -0.145 4.145 19 C -0.154 4.154 20 C -0.156 4.156 21 C -0.154 4.154 22 C -0.147 4.147 23 C -0.094 4.094 24 C -0.109 4.109 25 C -0.213 4.213 26 C -0.064 4.064 27 C 0.119 3.881 28 H 0.078 0.922 29 C -0.138 4.138 30 C -0.141 4.141 31 C -0.062 4.062 32 H 0.059 0.941 33 H 0.104 0.896 34 H 0.040 0.960 35 H 0.084 0.916 36 H 0.096 0.904 37 H 0.141 0.859 38 H 0.152 0.848 39 H 0.097 0.903 40 H 0.078 0.922 41 H 0.075 0.925 42 H 0.083 0.917 43 H 0.086 0.914 44 H 0.088 0.912 45 H 0.078 0.922 46 H 0.083 0.917 47 H 0.084 0.916 48 H 0.070 0.930 49 H 0.078 0.922 50 H 0.159 0.841 51 H 0.060 0.940 52 H 0.088 0.912 53 H 0.067 0.933 54 H 0.108 0.892 55 H 0.103 0.897 56 H 0.095 0.905 Dipole moment (debyes) X Y Z Total from point charges -7.174 -18.477 -5.596 20.596 hybrid contribution 0.044 1.935 0.342 1.965 sum -7.130 -16.542 -5.254 18.764 Atomic orbital electron populations 1.70074 1.06367 1.25791 1.54988 1.25710 0.96076 1.01167 1.01150 0.85280 0.78361 0.79823 0.77051 1.21878 1.21922 1.19399 1.19330 0.93167 0.84880 1.40691 0.90798 1.44676 1.34750 1.01947 1.51458 1.20846 0.85148 1.00506 0.90006 1.22318 1.00361 0.84424 0.94069 1.43723 1.05255 1.22336 1.42831 1.19101 0.92620 0.92886 1.00117 1.20516 0.90137 1.01473 1.04117 1.20702 0.99955 0.94207 0.96985 1.19482 0.92849 0.98443 1.00243 1.20073 0.96697 0.93338 0.95055 1.21567 0.98852 0.96783 0.97343 1.21531 0.98496 0.97216 0.98124 1.21538 0.97790 0.99058 0.97191 1.21527 0.99345 0.96780 0.97734 1.21577 0.97802 0.97849 0.97506 1.20619 0.93985 0.96343 0.98455 1.18900 0.91907 0.98452 1.01609 1.20252 0.94466 1.01184 1.05432 1.21767 0.87539 0.97193 0.99932 1.19662 0.88759 0.85221 0.94424 0.92172 1.20882 0.98836 0.98868 0.95217 1.20885 1.01067 0.94181 0.97922 1.19774 0.88996 0.96750 1.00688 0.94063 0.89608 0.96022 0.91605 0.90423 0.85863 0.84779 0.90341 0.92169 0.92539 0.91668 0.91413 0.91184 0.92239 0.91680 0.91642 0.93005 0.92228 0.84091 0.94005 0.91216 0.93314 0.89212 0.89680 0.90510 Number of geometries 1 Number of calculations of the screened coulomb radii 1 The total number of SCF iterations 18. 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 N -0.93 -51.63 10.20 21.65 0.22 -51.41 16 2 C -0.05 -2.41 4.94 84.86 0.42 -1.99 16 3 Si 0.97 43.59 29.59 68.60 2.03 45.62 16 4 H -0.29 -13.26 7.11 99.48 0.71 -12.55 16 5 H -0.29 -13.28 7.11 99.48 0.71 -12.57 16 6 H -0.27 -11.64 7.11 99.48 0.71 -10.93 16 7 C -0.25 -13.17 8.22 84.03 0.69 -12.48 16 8 H 0.07 3.77 7.27 -2.91 -0.02 3.75 16 9 N -0.61 -29.83 2.99 -678.16 -2.03 -31.86 16 10 C 0.16 7.63 5.39 86.28 0.47 8.10 16 11 C 0.11 4.97 7.15 86.32 0.62 5.58 16 12 N -0.46 -20.41 3.19 -497.98 -1.59 -21.99 16 13 C 0.06 2.42 6.87 39.76 0.27 2.70 16 14 C -0.14 -4.87 10.23 22.54 0.23 -4.64 16 15 C -0.10 -2.81 8.55 22.36 0.19 -2.62 16 16 C -0.11 -2.90 4.77 -20.04 -0.10 -3.00 16 17 C -0.03 -0.66 2.32 -11.54 -0.03 -0.68 16 18 C -0.11 -1.88 5.01 30.60 0.15 -1.72 16 19 C -0.12 -1.58 5.92 30.60 0.18 -1.40 16 20 C -0.12 -1.52 5.92 30.59 0.18 -1.34 16 21 C -0.12 -1.67 5.92 30.59 0.18 -1.49 16 22 C -0.11 -2.03 5.00 30.59 0.15 -1.87 16 23 C -0.08 -2.27 9.79 22.42 0.22 -2.05 16 24 C -0.10 -3.87 6.27 -20.39 -0.13 -4.00 16 25 C -0.21 -7.79 6.38 -19.08 -0.12 -7.92 16 26 C 0.05 2.03 6.10 41.29 0.25 2.28 16 27 C 0.23 9.59 3.61 43.86 0.16 9.75 16 28 H 0.06 2.45 8.14 -2.39 -0.02 2.44 16 29 C -0.10 -3.76 4.42 29.60 0.13 -3.63 16 30 C -0.10 -3.04 5.80 29.75 0.17 -2.87 16 31 C -0.02 -0.71 6.17 29.80 0.18 -0.53 16 32 H 0.25 13.65 8.31 -92.71 -0.77 12.88 16 33 H 0.09 4.43 5.97 -2.39 -0.01 4.41 16 34 H 0.02 0.97 8.14 -2.38 -0.02 0.95 16 35 H 0.07 3.37 7.18 -2.39 -0.02 3.35 16 36 H 0.08 3.03 8.14 -2.39 -0.02 3.02 16 37 H 0.12 3.95 8.06 -2.91 -0.02 3.92 16 38 H 0.13 3.10 6.84 -2.91 -0.02 3.08 16 39 H 0.08 1.46 7.83 -2.39 -0.02 1.44 16 40 H 0.06 1.08 7.27 -2.39 -0.02 1.06 16 41 H 0.06 1.05 8.14 -2.39 -0.02 1.03 16 42 H 0.06 0.84 8.14 -2.39 -0.02 0.82 16 43 H 0.07 0.82 8.14 -2.39 -0.02 0.80 16 44 H 0.07 0.75 8.14 -2.39 -0.02 0.73 16 45 H 0.06 0.79 8.14 -2.39 -0.02 0.77 16 46 H 0.06 0.89 8.14 -2.39 -0.02 0.87 16 47 H 0.06 0.87 8.14 -2.39 -0.02 0.85 16 48 H 0.05 1.06 8.14 -2.39 -0.02 1.05 16 49 H 0.06 1.13 7.26 -2.39 -0.02 1.12 16 50 H 0.14 3.49 7.93 -2.91 -0.02 3.47 16 51 H 0.04 1.88 7.07 -2.39 -0.02 1.87 16 52 H 0.07 2.41 8.08 -2.39 -0.02 2.39 16 53 H 0.05 1.49 8.14 -2.38 -0.02 1.47 16 54 H 0.09 2.06 8.14 -2.39 -0.02 2.04 16 55 H 0.08 2.08 8.14 -2.38 -0.02 2.06 16 56 H 0.08 1.99 8.14 -2.39 -0.02 1.97 16 Total: -1.00 -61.88 413.26 3.97 -57.91 By element: Atomic # 1 Polarization: 26.70 SS G_CDS: 0.85 Total: 27.55 kcal Atomic # 6 Polarization: -30.30 SS G_CDS: 4.48 Total: -25.82 kcal Atomic # 7 Polarization: -101.87 SS G_CDS: -3.39 Total: -105.26 kcal Atomic # 14 Polarization: 43.59 SS G_CDS: 2.03 Total: 45.62 kcal Total: -61.88 3.97 -57.91 kcal The number of atoms in the molecule is 56 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. ZINC000019201937.mol2 56 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 104.380 kcal (2) G-P(sol) polarization free energy of solvation -61.877 kcal (3) G-ENP(sol) elect.-nuc.-pol. free energy of system 42.503 kcal (4) G-CDS(sol) cavity-dispersion-solvent structure free energy 3.970 kcal (5) G-P-CDS(sol) = G-P(sol) + G-CDS(sol) = (2) + (4) -57.907 kcal (6) G-S(sol) free energy of system = (1) + (5) 46.473 kcal 1SCF WAS SPECIFIED, SO GEOMETRY OPTIMIZATION WAS NOT USED Total computer time = 1.81 seconds