Wall clock time and date at job start Tue Mar 31 2020 23:54:21 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.42899 * 1 3 3 C 1.42898 * 114.00525 * 2 1 4 4 C 1.53001 * 109.47293 * 179.97438 * 3 2 1 5 5 N 1.46500 * 109.46559 * 65.00743 * 4 3 2 6 6 C 1.36939 * 119.99794 * 264.42931 * 5 4 3 7 7 H 1.08006 * 120.00064 * 146.14050 * 6 5 4 8 8 C 1.38814 * 120.00531 * 326.14863 * 6 5 4 9 9 N 1.31622 * 119.99273 * 337.66455 * 8 6 5 10 10 Si 1.86790 * 120.00353 * 157.66549 * 8 6 5 11 11 H 1.48500 * 109.47323 * 60.00179 * 10 8 6 12 12 H 1.48502 * 109.47299 * 180.02562 * 10 8 6 13 13 H 1.48498 * 109.47477 * 299.99739 * 10 8 6 14 14 C 1.46505 * 120.00425 * 84.43667 * 5 4 3 15 15 C 1.53030 * 109.45871 * 239.99717 * 14 5 4 16 16 C 1.53194 * 109.31141 * 178.61297 * 15 14 5 17 17 N 1.46921 * 108.77623 * 54.63660 * 16 15 14 18 18 C 1.34774 * 120.63489 * 126.36749 * 17 16 15 19 19 O 1.21287 * 119.99804 * 0.02562 * 18 17 16 20 20 C 1.50700 * 120.00178 * 179.97438 * 18 17 16 21 21 H 1.08997 * 110.45704 * 57.40291 * 20 18 17 22 22 C 1.54472 * 110.51267 * 294.64254 * 20 18 17 23 23 O 1.44430 * 104.80520 * 217.37612 * 22 20 18 24 24 C 1.44422 * 105.27036 * 40.51328 * 23 22 20 25 25 H 1.09002 * 110.37519 * 78.32144 * 24 23 22 26 26 C 1.53002 * 110.46472 * 200.60943 * 24 23 22 27 27 C 1.54473 * 104.80644 * 319.48638 * 24 23 22 28 28 C 1.46932 * 118.73541 * 306.64559 * 17 16 15 29 29 C 1.52752 * 109.01025 * 53.30819 * 28 17 16 30 30 H 1.08993 * 109.47567 * 179.97438 * 1 2 3 31 31 H 1.09006 * 109.46577 * 299.99620 * 1 2 3 32 32 H 1.09003 * 109.47077 * 59.99244 * 1 2 3 33 33 H 1.09001 * 109.46550 * 60.00565 * 3 2 1 34 34 H 1.08999 * 109.47344 * 300.00827 * 3 2 1 35 35 H 1.08997 * 109.47748 * 305.00560 * 4 3 2 36 36 H 1.08999 * 109.46810 * 185.00243 * 4 3 2 37 37 H 0.97010 * 119.99656 * 179.97438 * 9 8 6 38 38 H 1.08993 * 109.45579 * 359.97383 * 14 5 4 39 39 H 1.08999 * 109.49867 * 58.66195 * 15 14 5 40 40 H 1.08998 * 109.53046 * 298.57903 * 15 14 5 41 41 H 1.08993 * 109.58652 * 294.84445 * 16 15 14 42 42 H 1.09004 * 109.70127 * 174.49580 * 16 15 14 43 43 H 1.08999 * 110.37126 * 336.20834 * 22 20 18 44 44 H 1.09006 * 110.37208 * 98.54609 * 22 20 18 45 45 H 1.08993 * 109.47390 * 175.43841 * 26 24 23 46 46 H 1.09005 * 109.47330 * 295.44069 * 26 24 23 47 47 H 1.09000 * 109.47327 * 55.43269 * 26 24 23 48 48 H 1.09002 * 110.51647 * 265.33613 * 27 24 23 49 49 H 1.09007 * 110.62607 * 142.68410 * 27 24 23 50 50 H 1.09005 * 109.58166 * 173.14433 * 28 17 16 51 51 H 1.09003 * 109.58562 * 293.56449 * 28 17 16 52 52 H 1.09001 * 109.56732 * 185.55139 * 29 28 17 53 53 H 1.09004 * 109.35958 * 65.49878 * 29 28 17 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.0103 1.3054 0.0000 4 6 3.5356 1.1845 0.0006 5 7 3.9766 0.5639 -1.2510 6 6 4.2075 -0.7849 -1.3023 7 1 3.9961 -1.3349 -2.2074 8 6 4.7130 -1.4454 -0.1909 9 7 5.3725 -0.7745 0.7296 10 14 4.4574 -3.2860 -0.0018 11 1 5.0771 -3.9926 -1.1516 12 1 5.0885 -3.7463 1.2611 13 1 3.0028 -3.5826 0.0357 14 6 4.1709 1.3863 -2.4478 15 6 3.2428 0.8943 -3.5606 16 6 3.4236 1.7780 -4.7988 17 7 4.8467 1.7954 -5.1636 18 6 5.2394 1.4765 -6.4128 19 8 4.4148 1.1718 -7.2485 20 6 6.7028 1.4997 -6.7717 21 1 7.2663 0.8354 -6.1165 22 6 7.2566 2.9396 -6.6936 23 8 8.1890 3.0365 -7.7923 24 6 7.5555 2.3454 -8.8909 25 1 6.7953 2.9780 -9.3493 26 6 8.5965 1.9267 -9.9310 27 6 6.9040 1.0998 -8.2506 28 6 5.8342 2.1612 -4.1390 29 6 5.6259 1.2786 -2.9097 30 1 -0.3634 -1.0276 -0.0005 31 1 -0.3633 0.5138 0.8901 32 1 -0.3633 0.5140 -0.8899 33 1 1.6888 1.8462 0.8901 34 1 1.6887 1.8464 -0.8899 35 1 3.8509 0.5683 0.8426 36 1 3.9786 2.1764 0.0903 37 1 5.7255 -1.2360 1.5065 38 1 3.9400 2.4258 -2.2154 39 1 2.2081 0.9493 -3.2223 40 1 3.4892 -0.1375 -3.8110 41 1 3.0910 2.7917 -4.5759 42 1 2.8391 1.3754 -5.6261 43 1 6.4506 3.6632 -6.8150 44 1 7.7701 3.0978 -5.7451 45 1 8.1143 1.3375 -10.7110 46 1 9.0466 2.8157 -10.3731 47 1 9.3706 1.3287 -9.4501 48 1 7.5688 0.2393 -8.3265 49 1 5.9451 0.8830 -8.7215 50 1 6.8402 2.0090 -4.5302 51 1 5.7023 3.2069 -3.8611 52 1 6.2874 1.6071 -2.1081 53 1 5.8505 0.2432 -3.1661 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 ZINC001269644559.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:54:21 Heat of formation + Delta-G solvation = -101.740795 kcal Electronic energy + Delta-G solvation = -31862.273596 eV Core-core repulsion = 27734.899525 eV Total energy + Delta-G solvation = -4127.374071 eV No. of doubly occupied orbitals = 66 Molecular weight (most abundant/longest-lived isotopes) = 340.220 amu Computer time = 1.58 seconds Orbital eigenvalues (eV) -40.70564 -38.95012 -37.62407 -36.85791 -35.01256 -32.71879 -32.35131 -31.68675 -30.59445 -30.39950 -28.16997 -26.24002 -25.69144 -25.11985 -23.73735 -22.79707 -22.48092 -21.43237 -20.56853 -19.82362 -18.77500 -17.66329 -16.44481 -16.22814 -15.86728 -15.77937 -15.59701 -15.17221 -15.05999 -14.84302 -14.66390 -14.04715 -13.90747 -13.79299 -13.53841 -13.50744 -13.21677 -13.04067 -12.63792 -12.41975 -12.30478 -12.21688 -11.83338 -11.79118 -11.65640 -11.47788 -11.36405 -11.04800 -10.94422 -10.91934 -10.79391 -10.74857 -10.63477 -10.48155 -10.08252 -9.98045 -9.90959 -9.82846 -9.64586 -9.44437 -8.86054 -8.44760 -8.37117 -6.57497 -5.73022 -3.23190 2.65567 3.07637 3.44416 3.49580 3.56183 3.69104 4.27977 4.52365 4.54235 4.56005 4.65976 4.80489 4.94386 5.03955 5.09621 5.17845 5.19394 5.37975 5.40529 5.47665 5.49012 5.49120 5.55049 5.69207 5.70823 5.72184 5.76580 5.76952 5.82803 5.89700 5.94399 5.95321 5.99718 6.04709 6.10993 6.12014 6.24588 6.36927 6.43411 6.59208 6.61674 6.65824 6.72626 6.80692 7.06793 7.13657 7.54235 7.67024 7.88408 8.04658 8.23997 8.57300 9.23845 9.89130 10.41759 11.29924 Molecular weight = 340.22amu Principal moments of inertia in cm(-1) A = 0.015291 B = 0.002954 C = 0.002757 Principal moments of inertia in units of 10**(-40)*gram-cm**2 A = 1830.699255 B = 9477.266983 C =10152.266373 Net atomic charges, atomic populations, and dipole contributions calculated with CM2 Atom NO. Type Charge No. of electrons 1 C 0.030 3.970 2 O -0.356 6.356 3 C 0.045 3.955 4 C 0.142 3.858 5 N -0.571 5.571 6 C -0.234 4.234 7 H 0.057 0.943 8 C 0.015 3.985 9 N -0.894 5.894 10 Si 0.894 3.106 11 H -0.286 1.286 12 H -0.293 1.293 13 H -0.282 1.282 14 C 0.178 3.822 15 C -0.139 4.139 16 C 0.114 3.886 17 N -0.600 5.600 18 C 0.523 3.477 19 O -0.548 6.548 20 C -0.146 4.146 21 H 0.111 0.889 22 C 0.060 3.940 23 O -0.370 6.370 24 C 0.069 3.931 25 H 0.069 0.931 26 C -0.146 4.146 27 C -0.141 4.141 28 C 0.102 3.898 29 C -0.136 4.136 30 H 0.074 0.926 31 H 0.026 0.974 32 H 0.028 0.972 33 H 0.024 0.976 34 H 0.026 0.974 35 H 0.110 0.890 36 H 0.037 0.963 37 H 0.252 0.748 38 H 0.047 0.953 39 H 0.066 0.934 40 H 0.068 0.932 41 H 0.060 0.940 42 H 0.084 0.916 43 H 0.063 0.937 44 H 0.102 0.898 45 H 0.071 0.929 46 H 0.061 0.939 47 H 0.064 0.936 48 H 0.080 0.920 49 H 0.100 0.900 50 H 0.078 0.922 51 H 0.065 0.935 52 H 0.074 0.926 53 H 0.072 0.928 For the charges calculated by CM2: Dipole moment (debyes) X Y Z Total from point charges 1.082 8.896 -11.861 14.866 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.067 4.067 2 O -0.273 6.273 3 C -0.034 4.034 4 C 0.016 3.984 5 N -0.290 5.290 6 C -0.364 4.364 7 H 0.074 0.926 8 C -0.184 4.184 9 N -0.538 5.538 10 Si 0.724 3.276 11 H -0.212 1.212 12 H -0.220 1.220 13 H -0.209 1.209 14 C 0.070 3.930 15 C -0.179 4.179 16 C -0.009 4.009 17 N -0.332 5.332 18 C 0.310 3.690 19 O -0.426 6.426 20 C -0.167 4.167 21 H 0.129 0.871 22 C -0.018 4.018 23 O -0.289 6.289 24 C 0.012 3.988 25 H 0.087 0.913 26 C -0.203 4.203 27 C -0.179 4.179 28 C -0.021 4.021 29 C -0.176 4.176 30 H 0.092 0.908 31 H 0.045 0.955 32 H 0.047 0.953 33 H 0.042 0.958 34 H 0.044 0.956 35 H 0.127 0.873 36 H 0.055 0.945 37 H 0.062 0.938 38 H 0.065 0.935 39 H 0.085 0.915 40 H 0.087 0.913 41 H 0.079 0.921 42 H 0.102 0.898 43 H 0.081 0.919 44 H 0.120 0.880 45 H 0.090 0.910 46 H 0.081 0.919 47 H 0.083 0.917 48 H 0.099 0.901 49 H 0.119 0.881 50 H 0.096 0.904 51 H 0.084 0.916 52 H 0.093 0.907 53 H 0.090 0.910 Dipole moment (debyes) X Y Z Total from point charges 1.060 8.612 -12.117 14.904 hybrid contribution -0.704 -0.879 0.438 1.208 sum 0.356 7.733 -11.679 14.012 Atomic orbital electron populations 1.22628 0.84233 1.01470 0.98387 1.88208 1.16555 1.27752 1.94818 1.23473 0.99109 0.84052 0.96744 1.21081 0.87136 0.96825 0.93407 1.43750 1.75658 1.04369 1.05206 1.18747 1.35590 0.82726 0.99317 0.92556 1.24912 0.96984 0.97420 0.99066 1.69949 1.36075 1.35347 1.12423 0.86511 0.78206 0.86114 0.76737 1.21241 1.21953 1.20866 1.20294 0.91487 0.93394 0.87832 1.22355 1.00931 1.00090 0.94504 1.21692 0.80592 0.99995 0.98628 1.48123 1.08339 1.66893 1.09884 1.20232 0.90391 0.76469 0.81887 1.90637 1.54320 1.49054 1.48627 1.22674 0.94470 0.99915 0.99650 0.87073 1.23429 0.92970 0.93718 0.91639 1.89137 1.41518 1.76815 1.21437 1.23027 0.94748 0.92282 0.88754 0.91290 1.21797 0.98837 1.01860 0.97783 1.23020 1.03354 0.98505 0.93051 1.21964 0.91429 0.99342 0.89394 1.22254 1.00149 1.00228 0.94990 0.90759 0.95532 0.95330 0.95824 0.95602 0.87281 0.94471 0.93771 0.93493 0.91505 0.91319 0.92129 0.89779 0.91893 0.87970 0.91044 0.91948 0.91656 0.90129 0.88117 0.90405 0.91631 0.90708 0.90952 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 C 0.03 0.49 11.01 101.05 1.11 1.61 16 2 O -0.36 -7.62 10.44 -35.23 -0.37 -7.99 16 3 C 0.05 0.88 6.28 37.16 0.23 1.11 16 4 C 0.14 3.13 5.57 -4.04 -0.02 3.10 16 5 N -0.57 -12.89 2.35 -165.62 -0.39 -13.28 16 6 C -0.23 -6.03 8.06 -15.06 -0.12 -6.15 16 7 H 0.06 1.50 6.87 -52.48 -0.36 1.14 16 8 C 0.02 0.40 4.91 -17.43 -0.09 0.32 16 9 N -0.89 -24.72 11.34 55.49 0.63 -24.09 16 10 Si 0.89 21.90 29.61 -169.99 -5.03 16.86 16 11 H -0.29 -6.95 7.11 56.52 0.40 -6.55 16 12 H -0.29 -7.23 7.11 56.52 0.40 -6.83 16 13 H -0.28 -7.04 7.11 56.52 0.40 -6.64 16 14 C 0.18 3.27 2.60 -67.89 -0.18 3.10 16 15 C -0.14 -2.41 5.30 -26.62 -0.14 -2.55 16 16 C 0.11 1.69 6.44 -3.71 -0.02 1.66 16 17 N -0.60 -8.57 2.98 -173.03 -0.52 -9.08 16 18 C 0.52 7.70 6.95 -10.99 -0.08 7.62 16 19 O -0.55 -9.68 14.69 5.55 0.08 -9.60 16 20 C -0.15 -1.72 2.54 -90.21 -0.23 -1.95 16 21 H 0.11 1.25 7.57 -51.93 -0.39 0.86 16 22 C 0.06 0.64 6.12 37.93 0.23 0.87 16 23 O -0.37 -4.40 11.35 -35.23 -0.40 -4.80 16 24 C 0.07 0.72 3.92 -25.96 -0.10 0.62 16 25 H 0.07 0.71 8.14 -51.93 -0.42 0.28 16 26 C -0.15 -1.20 9.52 37.16 0.35 -0.85 16 27 C -0.14 -1.58 6.21 -25.16 -0.16 -1.73 16 28 C 0.10 1.31 5.71 -3.93 -0.02 1.29 16 29 C -0.14 -2.32 5.29 -26.84 -0.14 -2.46 16 30 H 0.07 1.13 8.14 -51.93 -0.42 0.71 16 31 H 0.03 0.38 8.14 -51.93 -0.42 -0.04 16 32 H 0.03 0.42 8.14 -51.93 -0.42 0.00 16 33 H 0.02 0.41 8.14 -51.93 -0.42 -0.01 16 34 H 0.03 0.47 8.14 -51.93 -0.42 0.05 16 35 H 0.11 2.81 5.34 -51.93 -0.28 2.53 16 36 H 0.04 0.75 8.04 -51.93 -0.42 0.34 16 37 H 0.25 6.88 8.31 -40.82 -0.34 6.54 16 38 H 0.05 0.81 7.87 -51.93 -0.41 0.40 16 39 H 0.07 1.19 8.14 -51.93 -0.42 0.76 16 40 H 0.07 1.37 7.03 -51.93 -0.37 1.01 16 41 H 0.06 0.80 8.14 -51.93 -0.42 0.38 16 42 H 0.08 1.32 7.00 -51.93 -0.36 0.96 16 43 H 0.06 0.68 7.86 -51.93 -0.41 0.27 16 44 H 0.10 0.93 6.88 -51.93 -0.36 0.57 16 45 H 0.07 0.49 8.14 -51.93 -0.42 0.07 16 46 H 0.06 0.48 8.14 -51.93 -0.42 0.06 16 47 H 0.06 0.54 8.14 -51.93 -0.42 0.12 16 48 H 0.08 0.81 7.87 -51.93 -0.41 0.40 16 49 H 0.10 1.27 6.50 -51.93 -0.34 0.93 16 50 H 0.08 0.86 4.52 -51.93 -0.23 0.62 16 51 H 0.07 0.75 8.14 -51.93 -0.42 0.32 16 52 H 0.07 1.36 8.14 -51.93 -0.42 0.93 16 53 H 0.07 1.43 7.71 -51.93 -0.40 1.03 16 LS Contribution 405.74 15.07 6.11 6.11 Total: -1.00 -30.44 405.74 -8.61 -39.05 By element: Atomic # 1 Polarization: 10.58 SS G_CDS: -9.36 Total: 1.21 kcal Atomic # 6 Polarization: 4.96 SS G_CDS: 0.63 Total: 5.59 kcal Atomic # 7 Polarization: -46.17 SS G_CDS: -0.27 Total: -46.45 kcal Atomic # 8 Polarization: -21.70 SS G_CDS: -0.69 Total: -22.39 kcal Atomic # 14 Polarization: 21.90 SS G_CDS: -5.03 Total: 16.86 kcal Total LS contribution 6.11 Total: 6.11 kcal Total: -30.44 -8.61 -39.05 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. ZINC001269644559.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 -62.688 kcal (2) G-P(sol) polarization free energy of solvation -30.443 kcal (3) G-ENP(sol) elect.-nuc.-pol. free energy of system -93.131 kcal (4) G-CDS(sol) cavity-dispersion-solvent structure free energy -8.610 kcal (5) G-P-CDS(sol) = G-P(sol) + G-CDS(sol) = (2) + (4) -39.053 kcal (6) G-S(sol) free energy of system = (1) + (5) -101.741 kcal 1SCF WAS SPECIFIED, SO GEOMETRY OPTIMIZATION WAS NOT USED Total computer time = 1.58 seconds