Wall clock time and date at job start Sat Apr 11 2020 06:13:09 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 N 1.46898 * 1 3 3 C 1.46899 * 111.00143 * 2 1 4 4 C 1.46905 * 111.00058 * 123.95242 * 2 1 3 5 5 C 1.52997 * 109.47088 * 66.04514 * 4 2 1 6 6 C 1.50696 * 109.47364 * 180.02562 * 5 4 2 7 7 O 1.21283 * 120.00072 * 359.97438 * 6 5 4 8 8 N 1.34779 * 120.00186 * 179.97438 * 6 5 4 9 9 C 1.46925 * 120.63041 * 359.97438 * 8 6 5 10 10 C 1.53191 * 108.77423 * 233.62788 * 9 8 6 11 11 C 1.53037 * 109.31499 * 305.36620 * 10 9 8 12 12 C 1.53043 * 109.53919 * 61.36673 * 11 10 9 13 13 H 1.09000 * 109.49902 * 58.58363 * 12 11 10 14 14 N 1.46894 * 109.49819 * 178.68485 * 12 11 10 15 15 C 1.46902 * 111.00300 * 65.88485 * 14 12 11 16 16 C 1.46903 * 111.00196 * 189.83590 * 14 12 11 17 17 C 1.50694 * 109.47444 * 169.99850 * 16 14 12 18 18 O 1.21292 * 120.00041 * 0.02562 * 17 16 14 19 19 N 1.34777 * 120.00196 * 179.97438 * 17 16 14 20 20 C 1.37097 * 120.00161 * 179.97438 * 19 17 16 21 21 H 1.08004 * 120.00071 * 0.02562 * 20 19 17 22 22 C 1.38944 * 120.00058 * 179.97438 * 20 19 17 23 23 N 1.31415 * 120.00025 * 0.02562 * 22 20 19 24 24 Si 1.86804 * 120.00122 * 180.02562 * 22 20 19 25 25 H 1.48503 * 109.46702 * 60.00320 * 24 22 20 26 26 H 1.48500 * 109.47223 * 180.02562 * 24 22 20 27 27 H 1.48500 * 109.47159 * 300.00429 * 24 22 20 28 28 C 1.46926 * 120.63044 * 180.02562 * 8 6 5 29 29 H 1.09004 * 109.46918 * 296.04609 * 1 2 3 30 30 H 1.08994 * 109.47430 * 56.04774 * 1 2 3 31 31 H 1.09003 * 109.47402 * 176.04787 * 1 2 3 32 32 H 1.09000 * 109.47467 * 183.95415 * 3 2 1 33 33 H 1.09003 * 109.47055 * 303.95130 * 3 2 1 34 34 H 1.09003 * 109.47298 * 63.95196 * 3 2 1 35 35 H 1.08998 * 109.47053 * 306.04702 * 4 2 1 36 36 H 1.08997 * 109.47368 * 186.04751 * 4 2 1 37 37 H 1.08995 * 109.47395 * 59.99838 * 5 4 2 38 38 H 1.09006 * 109.47337 * 300.00211 * 5 4 2 39 39 H 1.08998 * 109.58340 * 353.41744 * 9 8 6 40 40 H 1.09003 * 109.58494 * 113.70309 * 9 8 6 41 41 H 1.08997 * 109.49512 * 65.31959 * 10 9 8 42 42 H 1.09005 * 109.49659 * 185.37937 * 10 9 8 43 43 H 1.08997 * 109.46383 * 181.38282 * 11 10 9 44 44 H 1.09001 * 109.46057 * 301.34715 * 11 10 9 45 45 H 1.08994 * 109.47236 * 59.99501 * 15 14 12 46 46 H 1.09000 * 109.47303 * 180.02562 * 15 14 12 47 47 H 1.09007 * 109.46732 * 299.99949 * 15 14 12 48 48 H 1.09003 * 109.46842 * 289.99818 * 16 14 12 49 49 H 1.09000 * 109.46972 * 49.99441 * 16 14 12 50 50 H 0.96995 * 119.99887 * 0.02562 * 19 17 16 51 51 H 0.97000 * 119.99829 * 180.02562 * 23 22 20 52 52 H 1.09006 * 109.58406 * 246.15786 * 28 8 6 53 53 H 1.09002 * 109.59019 * 6.57775 * 28 8 6 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.4690 0.0000 0.0000 3 6 1.9955 1.3714 0.0000 4 6 1.9955 -0.7660 1.1376 5 6 1.6315 -2.2426 0.9704 6 6 2.1709 -3.0283 2.1377 7 8 2.7865 -2.4630 3.0166 8 7 1.9695 -4.3592 2.2050 9 6 1.2253 -5.0581 1.1483 10 6 0.0880 -5.8572 1.7922 11 6 0.6638 -6.7785 2.8700 12 6 1.3379 -5.9370 3.9561 13 1 0.6101 -5.2500 4.3882 14 7 1.8626 -6.8204 5.0059 15 6 2.9483 -7.6646 4.4895 16 6 2.3074 -6.0458 6.1721 17 6 2.6065 -6.9832 7.3135 18 8 2.4637 -8.1791 7.1695 19 7 3.0337 -6.4930 8.4940 20 6 3.3053 -7.3458 9.5325 21 1 3.1773 -8.4105 9.4045 22 6 3.7461 -6.8404 10.7494 23 7 3.9023 -5.5448 10.9050 24 14 4.1169 -8.0024 12.1642 25 1 5.1874 -8.9488 11.7598 26 1 4.5671 -7.2194 13.3430 27 1 2.8904 -8.7627 12.5147 28 6 2.4888 -5.1387 3.3370 29 1 -0.3633 0.4513 0.9233 30 1 -0.3634 0.5739 -0.8524 31 1 -0.3634 -1.0252 -0.0708 32 1 3.0828 1.3432 -0.0709 33 1 1.5898 1.9163 -0.8525 34 1 1.7044 1.8724 0.9233 35 1 1.5611 -0.3872 2.0628 36 1 3.0797 -0.6611 1.1757 37 1 2.0657 -2.6214 0.0453 38 1 0.5471 -2.3475 0.9324 39 1 0.8112 -4.3302 0.4506 40 1 1.8945 -5.7345 0.6164 41 1 -0.6279 -5.1710 2.2446 42 1 -0.4126 -6.4555 1.0309 43 1 -0.1403 -7.3662 3.3129 44 1 1.3979 -7.4478 2.4214 45 1 3.7590 -7.0325 4.1271 46 1 3.3182 -8.3090 5.2870 47 1 2.5729 -8.2789 3.6710 48 1 3.2075 -5.4871 5.9154 49 1 1.5211 -5.3516 6.4686 50 1 3.1482 -5.5367 8.6091 51 1 4.2105 -5.1921 11.7544 52 1 3.2607 -5.8235 2.9856 53 1 2.9084 -4.4641 4.0834 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 ZINC001491947108.mol2 53 1SCF WAS SPECIFIED, SO GEOMETRY OPTIMIZATION WAS NOT USED SCF FIELD WAS ACHIEVED Wall clock time and date at job start Sat Apr 11 2020 06:13:09 Heat of formation + Delta-G solvation = -51.917587 kcal Electronic energy + Delta-G solvation = -30549.827462 eV Core-core repulsion = 26456.906639 eV Total energy + Delta-G solvation = -4092.920823 eV No. of doubly occupied orbitals = 66 Molecular weight (most abundant/longest-lived isotopes) = 340.226 amu Computer time = 0.67 seconds Orbital eigenvalues (eV) -40.64329 -38.98549 -37.96744 -37.41078 -35.73463 -33.84000 -32.65094 -31.86094 -31.20985 -29.87703 -27.98136 -27.83017 -26.36047 -25.43799 -24.62908 -24.06831 -22.20399 -21.62449 -19.96855 -19.10447 -19.00616 -17.79372 -17.59353 -16.45536 -16.31241 -15.95772 -15.49870 -15.25963 -15.10543 -14.98552 -14.78735 -14.67363 -14.36633 -14.04558 -13.85109 -13.65543 -13.52940 -13.09384 -13.02644 -12.94330 -12.74936 -12.70098 -12.20084 -12.08853 -12.07876 -11.84050 -11.58863 -11.46100 -11.34855 -11.32153 -11.25880 -11.08768 -10.95303 -10.81695 -10.51460 -10.41674 -10.18032 -9.95659 -9.43604 -9.03540 -8.78798 -8.37592 -7.93618 -7.75161 -6.37887 -3.78895 2.44384 2.69401 3.22300 3.27834 3.32074 3.56227 3.93447 3.96733 4.25075 4.34362 4.50272 4.60156 4.67114 4.70330 4.75745 4.78167 4.86881 4.99453 5.06473 5.12116 5.15265 5.15595 5.21839 5.25991 5.28943 5.37463 5.48976 5.49280 5.51123 5.53537 5.54675 5.63559 5.63978 5.71716 5.79900 5.82664 5.97050 6.00269 6.14720 6.16264 6.25214 6.32663 6.57559 6.57665 6.59072 7.17438 7.27700 7.39345 7.58262 7.98981 8.09609 8.65520 9.21391 9.58333 10.10130 10.79753 Molecular weight = 340.23amu Principal moments of inertia in cm(-1) A = 0.014898 B = 0.002202 C = 0.001972 Principal moments of inertia in units of 10**(-40)*gram-cm**2 A = 1879.032692 B =12714.905868 C =14193.628107 Net atomic charges, atomic populations, and dipole contributions calculated with CM2 Atom NO. Type Charge No. of electrons 1 C 0.018 3.982 2 N -0.550 5.550 3 C 0.025 3.975 4 C 0.068 3.932 5 C -0.135 4.135 6 C 0.514 3.486 7 O -0.540 6.540 8 N -0.608 5.608 9 C 0.105 3.895 10 C -0.133 4.133 11 C -0.111 4.111 12 C 0.073 3.927 13 H 0.083 0.917 14 N -0.505 5.505 15 C 0.014 3.986 16 C 0.042 3.958 17 C 0.445 3.555 18 O -0.578 6.578 19 N -0.564 5.564 20 C -0.299 4.299 21 H 0.103 0.897 22 C 0.021 3.979 23 N -0.902 5.902 24 Si 0.931 3.069 25 H -0.274 1.274 26 H -0.285 1.285 27 H -0.274 1.274 28 C 0.075 3.925 29 H 0.027 0.973 30 H 0.068 0.932 31 H 0.069 0.931 32 H 0.070 0.930 33 H 0.069 0.931 34 H 0.027 0.973 35 H 0.044 0.956 36 H 0.083 0.917 37 H 0.101 0.899 38 H 0.093 0.907 39 H 0.085 0.915 40 H 0.070 0.930 41 H 0.068 0.932 42 H 0.079 0.921 43 H 0.078 0.922 44 H 0.066 0.934 45 H 0.019 0.981 46 H 0.097 0.903 47 H 0.062 0.938 48 H 0.049 0.951 49 H 0.087 0.913 50 H 0.392 0.608 51 H 0.272 0.728 52 H 0.073 0.927 53 H 0.093 0.907 For the charges calculated by CM2: Dipole moment (debyes) X Y Z Total from point charges -7.229 9.356 -25.391 28.009 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.130 4.130 2 N -0.274 5.274 3 C -0.123 4.123 4 C -0.060 4.060 5 C -0.175 4.175 6 C 0.301 3.699 7 O -0.417 6.417 8 N -0.342 5.342 9 C -0.018 4.018 10 C -0.172 4.172 11 C -0.149 4.149 12 C -0.037 4.037 13 H 0.101 0.899 14 N -0.227 5.227 15 C -0.133 4.133 16 C -0.089 4.089 17 C 0.230 3.770 18 O -0.463 6.463 19 N -0.200 5.200 20 C -0.428 4.428 21 H 0.120 0.880 22 C -0.180 4.180 23 N -0.543 5.543 24 Si 0.758 3.242 25 H -0.199 1.199 26 H -0.210 1.210 27 H -0.199 1.199 28 C -0.049 4.049 29 H 0.046 0.954 30 H 0.087 0.913 31 H 0.088 0.912 32 H 0.088 0.912 33 H 0.087 0.913 34 H 0.046 0.954 35 H 0.063 0.937 36 H 0.101 0.899 37 H 0.119 0.881 38 H 0.111 0.889 39 H 0.103 0.897 40 H 0.089 0.911 41 H 0.087 0.913 42 H 0.097 0.903 43 H 0.096 0.904 44 H 0.085 0.915 45 H 0.037 0.963 46 H 0.115 0.885 47 H 0.081 0.919 48 H 0.067 0.933 49 H 0.105 0.895 50 H 0.225 0.775 51 H 0.082 0.918 52 H 0.091 0.909 53 H 0.111 0.889 Dipole moment (debyes) X Y Z Total from point charges -7.143 8.424 -25.097 27.420 hybrid contribution 0.599 1.009 1.299 1.750 sum -6.544 9.433 -23.798 26.423 Atomic orbital electron populations 1.22468 0.87493 1.02347 1.00687 1.61981 1.08024 1.14161 1.43208 1.22371 1.00959 0.88418 1.00600 1.21855 1.01199 0.91514 0.91433 1.21645 1.02583 0.95523 0.97749 1.20653 0.79772 0.81836 0.87613 1.90729 1.41293 1.69093 1.40570 1.48086 1.55390 1.07041 1.23719 1.21811 0.93846 0.96573 0.89549 1.21921 0.97520 0.98892 0.98882 1.21559 1.00832 0.97592 0.94896 1.22108 0.95937 0.93663 0.92025 0.89894 1.61872 1.35369 1.20839 1.04657 1.22645 0.91749 0.97290 1.01645 1.21911 0.99940 0.95567 0.91460 1.19695 0.84111 0.89131 0.84028 1.90476 1.57395 1.15834 1.82571 1.41866 1.58135 1.09938 1.10104 1.19640 1.41399 0.92083 0.89716 0.87965 1.25270 0.97417 0.96398 0.98927 1.69613 1.46750 1.11425 1.26489 0.85918 0.78977 0.79516 0.79788 1.19937 1.21011 1.19912 1.22500 0.97781 0.95344 0.89322 0.95429 0.91307 0.91186 0.91168 0.91270 0.95429 0.93736 0.89857 0.88083 0.88900 0.89654 0.91138 0.91348 0.90290 0.90372 0.91497 0.96308 0.88473 0.91919 0.93277 0.89457 0.77462 0.91831 0.90889 0.88882 Number of geometries 1 Number of calculations of the screened coulomb radii 1 The total number of SCF iterations 13. 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.02 0.11 9.32 60.06 0.56 0.67 16 2 N -0.55 -4.29 5.02 -251.26 -1.26 -5.55 16 3 C 0.02 0.17 9.57 60.06 0.57 0.74 16 4 C 0.07 0.62 4.25 -3.82 -0.02 0.60 16 5 C -0.13 -1.08 3.97 -27.89 -0.11 -1.19 16 6 C 0.51 6.00 7.70 -10.99 -0.08 5.91 16 7 O -0.54 -8.63 15.40 5.56 0.09 -8.55 16 8 N -0.61 -6.38 2.97 -172.80 -0.51 -6.90 16 9 C 0.11 0.73 6.34 -3.71 -0.02 0.71 16 10 C -0.13 -0.96 6.08 -26.61 -0.16 -1.12 16 11 C -0.11 -1.12 5.05 -26.69 -0.13 -1.26 16 12 C 0.07 0.94 2.34 -67.62 -0.16 0.78 16 13 H 0.08 1.13 7.94 -51.93 -0.41 0.71 16 14 N -0.50 -8.35 4.45 -236.21 -1.05 -9.40 16 15 C 0.01 0.23 8.43 60.07 0.51 0.73 16 16 C 0.04 0.79 5.02 -4.98 -0.03 0.76 16 17 C 0.44 10.68 7.48 -10.99 -0.08 10.60 16 18 O -0.58 -15.60 13.91 5.55 0.08 -15.52 16 19 N -0.56 -14.62 5.29 -10.97 -0.06 -14.68 16 20 C -0.30 -8.45 9.68 -14.94 -0.14 -8.59 16 21 H 0.10 3.01 7.16 -52.48 -0.38 2.63 16 22 C 0.02 0.59 5.50 -17.47 -0.10 0.49 16 23 N -0.90 -25.06 13.05 55.50 0.72 -24.33 16 24 Si 0.93 21.38 29.55 -169.99 -5.02 16.35 16 25 H -0.27 -6.23 7.11 56.52 0.40 -5.83 16 26 H -0.28 -6.39 7.11 56.52 0.40 -5.99 16 27 H -0.27 -6.24 7.11 56.52 0.40 -5.84 16 28 C 0.07 0.94 5.01 -3.74 -0.02 0.92 16 29 H 0.03 0.18 8.14 -51.93 -0.42 -0.24 16 30 H 0.07 0.36 8.12 -51.93 -0.42 -0.06 16 31 H 0.07 0.37 6.47 -51.93 -0.34 0.04 16 32 H 0.07 0.46 8.12 -51.93 -0.42 0.04 16 33 H 0.07 0.38 8.12 -51.93 -0.42 -0.04 16 34 H 0.03 0.19 8.14 -51.93 -0.42 -0.23 16 35 H 0.04 0.47 8.08 -51.93 -0.42 0.05 16 36 H 0.08 0.85 8.02 -51.93 -0.42 0.43 16 37 H 0.10 0.68 7.94 -51.93 -0.41 0.27 16 38 H 0.09 0.59 6.10 -51.93 -0.32 0.27 16 39 H 0.09 0.41 5.93 -51.93 -0.31 0.10 16 40 H 0.07 0.45 8.14 -51.93 -0.42 0.02 16 41 H 0.07 0.52 8.14 -51.93 -0.42 0.09 16 42 H 0.08 0.42 8.14 -51.93 -0.42 0.00 16 43 H 0.08 0.86 8.14 -51.93 -0.42 0.44 16 44 H 0.07 0.66 6.61 -51.93 -0.34 0.32 16 45 H 0.02 0.27 6.10 -51.93 -0.32 -0.04 16 46 H 0.10 1.91 5.82 -51.93 -0.30 1.61 16 47 H 0.06 0.86 6.50 -51.93 -0.34 0.52 16 48 H 0.05 0.87 7.20 -51.93 -0.37 0.50 16 49 H 0.09 1.52 7.93 -51.93 -0.41 1.11 16 50 H 0.39 9.85 7.90 -40.82 -0.32 9.53 16 51 H 0.27 7.17 8.31 -40.82 -0.34 6.83 16 52 H 0.07 0.91 6.04 -51.93 -0.31 0.60 16 53 H 0.09 1.42 6.05 -51.93 -0.31 1.10 16 LS Contribution 406.01 15.07 6.12 6.12 Total: -1.00 -33.48 406.01 -9.29 -42.77 By element: Atomic # 1 Polarization: 17.89 SS G_CDS: -8.97 Total: 8.92 kcal Atomic # 6 Polarization: 10.19 SS G_CDS: 0.58 Total: 10.77 kcal Atomic # 7 Polarization: -58.70 SS G_CDS: -2.16 Total: -60.86 kcal Atomic # 8 Polarization: -24.23 SS G_CDS: 0.16 Total: -24.07 kcal Atomic # 14 Polarization: 21.38 SS G_CDS: -5.02 Total: 16.35 kcal Total LS contribution 6.12 Total: 6.12 kcal Total: -33.48 -9.29 -42.77 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. ZINC001491947108.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 -9.151 kcal (2) G-P(sol) polarization free energy of solvation -33.479 kcal (3) G-ENP(sol) elect.-nuc.-pol. free energy of system -42.631 kcal (4) G-CDS(sol) cavity-dispersion-solvent structure free energy -9.287 kcal (5) G-P-CDS(sol) = G-P(sol) + G-CDS(sol) = (2) + (4) -42.766 kcal (6) G-S(sol) free energy of system = (1) + (5) -51.918 kcal 1SCF WAS SPECIFIED, SO GEOMETRY OPTIMIZATION WAS NOT USED Total computer time = 0.67 seconds