Wall clock time and date at job start Tue Mar 31 2020 19:38:20 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 C 1.52990 * 1 3 3 H 1.09008 * 109.46963 * 2 1 4 4 C 1.52999 * 109.47248 * 119.99675 * 2 1 3 5 5 C 1.50695 * 109.47443 * 294.99969 * 4 2 1 6 6 O 1.21285 * 119.99901 * 0.02562 * 5 4 2 7 7 N 1.34776 * 120.00347 * 179.97438 * 5 4 2 8 8 C 1.46950 * 120.44220 * 359.72953 * 7 5 4 9 9 C 1.54061 * 108.29289 * 125.71318 * 8 7 5 10 10 H 1.09003 * 110.52407 * 172.28382 * 9 8 7 11 11 C 1.53376 * 109.34037 * 294.24982 * 9 8 7 12 12 C 1.53008 * 109.25647 * 60.54321 * 11 9 8 13 13 C 1.54193 * 110.61527 * 59.71471 * 12 11 9 14 14 C 1.53364 * 109.08196 * 301.01562 * 13 12 11 15 15 H 1.09003 * 110.59191 * 176.63074 * 14 13 12 16 16 C 1.46943 * 120.44468 * 180.02562 * 7 5 4 17 17 N 1.47565 * 107.70451 * 51.34799 * 9 8 7 18 18 C 1.36929 * 121.61335 * 120.34044 * 17 9 8 19 19 H 1.07999 * 119.99807 * 332.08128 * 18 17 9 20 20 C 1.38816 * 120.00064 * 152.07595 * 18 17 9 21 21 N 1.31612 * 119.99901 * 341.79228 * 20 18 17 22 22 Si 1.86802 * 119.99870 * 161.79026 * 20 18 17 23 23 H 1.48504 * 109.47026 * 90.00044 * 22 20 18 24 24 H 1.48502 * 109.46980 * 209.99849 * 22 20 18 25 25 H 1.48501 * 109.47200 * 330.00040 * 22 20 18 26 26 C 1.52996 * 109.47459 * 240.00084 * 2 1 3 27 27 H 1.09005 * 110.54704 * 301.42433 * 26 2 1 28 28 C 1.55229 * 110.78026 * 64.28655 * 26 2 1 29 29 C 1.54814 * 101.98745 * 155.50538 * 28 26 2 30 30 C 1.54096 * 104.30454 * 339.26090 * 29 28 26 31 31 O 1.43747 * 107.29413 * 357.90699 * 30 29 28 32 32 H 1.08999 * 109.47397 * 299.99919 * 1 2 3 33 33 H 1.08998 * 109.47414 * 60.00232 * 1 2 3 34 34 H 1.09000 * 109.47483 * 180.02562 * 1 2 3 35 35 H 1.09003 * 109.46983 * 55.00370 * 4 2 1 36 36 H 1.09007 * 109.46916 * 175.00092 * 4 2 1 37 37 H 1.09002 * 109.62210 * 245.48015 * 8 7 5 38 38 H 1.09002 * 109.62284 * 6.03719 * 8 7 5 39 39 H 1.09003 * 109.50631 * 180.47882 * 11 9 8 40 40 H 1.09001 * 109.50749 * 300.56275 * 11 9 8 41 41 H 1.08998 * 109.41725 * 180.68100 * 12 11 9 42 42 H 1.09005 * 109.41766 * 300.54471 * 12 11 9 43 43 H 1.09002 * 110.80589 * 180.39226 * 13 12 11 44 44 H 1.09002 * 108.37387 * 60.17582 * 13 12 11 45 45 H 1.08996 * 109.62657 * 353.96064 * 16 7 5 46 46 H 1.08996 * 109.62895 * 114.54260 * 16 7 5 47 47 H 0.96999 * 119.99676 * 180.02562 * 21 20 18 48 48 H 1.08994 * 110.88890 * 37.37211 * 28 26 2 49 49 H 1.09005 * 111.01993 * 273.69803 * 28 26 2 50 50 H 1.08998 * 110.48183 * 220.51985 * 29 28 26 51 51 H 1.08992 * 110.48774 * 98.00611 * 29 28 26 52 52 H 1.09000 * 109.87661 * 117.32963 * 30 29 28 53 53 H 1.08998 * 109.88287 * 238.48048 * 30 29 28 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 6 1.5299 0.0000 0.0000 3 1 1.8932 1.0277 0.0000 4 6 2.0399 -0.7212 1.2493 5 6 1.6413 0.0571 2.4766 6 8 1.0109 1.0870 2.3633 7 7 1.9845 -0.3918 3.7002 8 6 2.7421 -1.6420 3.8491 9 6 4.0025 -1.3471 4.6846 10 1 4.5215 -2.2728 4.9332 11 6 4.9242 -0.4062 3.8986 12 6 4.1872 0.9053 3.6196 13 6 3.7721 1.5862 4.9393 14 6 2.8653 0.6354 5.7304 15 1 2.6049 1.0689 6.6961 16 6 1.6004 0.3591 4.9035 17 7 3.5743 -0.6462 5.9106 18 6 3.8224 -1.1605 7.1551 19 1 3.9115 -2.2283 7.2902 20 6 3.9595 -0.3093 8.2431 21 7 4.2099 0.9685 8.0516 22 14 3.7884 -0.9776 9.9791 23 1 5.1176 -1.4191 10.4728 24 1 3.2582 0.0864 10.8691 25 1 2.8535 -2.1314 9.9755 26 6 2.0400 -0.7212 -1.2492 27 1 1.6658 -1.7446 -1.2797 28 6 1.6299 0.0467 -2.5344 29 6 2.7079 -0.4176 -3.5440 30 6 3.8651 -0.9220 -2.6604 31 8 3.4830 -0.7060 -1.2915 32 1 -0.3634 0.5138 0.8900 33 1 -0.3634 0.5138 -0.8900 34 1 -0.3634 -1.0276 0.0005 35 1 1.6053 -1.7198 1.2952 36 1 3.1263 -0.7993 1.2056 37 1 2.1272 -2.3843 4.3582 38 1 3.0312 -2.0142 2.8662 39 1 5.8209 -0.2025 4.4839 40 1 5.2040 -0.8751 2.9552 41 1 4.8357 1.5714 3.0506 42 1 3.2851 0.6986 3.0436 43 1 3.2406 2.5176 4.7440 44 1 4.6729 1.7946 5.5166 45 1 1.1386 1.3031 4.6141 46 1 0.8981 -0.2295 5.4937 47 1 4.3060 1.5631 8.8119 48 1 1.6826 1.1239 -2.3771 49 1 0.6341 -0.2486 -2.8654 50 1 3.0362 0.4178 -4.1624 51 1 2.3240 -1.2249 -4.1676 52 1 4.0292 -1.9852 -2.8360 53 1 4.7735 -0.3628 -2.8843 RHF calculation, no. of doubly occupied orbitals= 65 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 ZINC001039984158.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 19:38:20 Heat of formation + Delta-G solvation = -45.357582 kcal Electronic energy + Delta-G solvation = -32046.562769 eV Core-core repulsion = 28112.089710 eV Total energy + Delta-G solvation = -3934.473059 eV No. of doubly occupied orbitals = 65 Molecular weight (most abundant/longest-lived isotopes) = 336.220 amu Computer time = 1.20 seconds Orbital eigenvalues (eV) -40.63044 -39.63079 -37.78326 -34.78252 -33.93859 -32.62572 -32.28957 -31.76930 -30.44022 -28.81075 -27.70811 -26.40697 -25.43686 -23.80965 -23.25857 -21.79667 -21.40126 -21.13246 -20.31460 -19.12264 -18.15498 -17.16818 -16.72029 -16.32011 -15.72502 -15.45253 -15.23156 -14.97466 -14.62550 -14.09064 -13.99329 -13.63488 -13.53572 -13.38507 -13.23869 -12.83796 -12.82162 -12.41067 -12.35613 -12.09999 -11.99593 -11.81708 -11.70062 -11.59783 -11.38862 -11.21497 -11.15290 -10.91264 -10.86900 -10.76871 -10.69151 -10.61539 -10.30865 -10.09610 -9.98675 -9.85855 -9.78048 -9.49119 -9.39801 -9.07622 -8.48545 -7.80154 -6.58892 -5.81141 -3.24386 3.06150 3.35663 3.43593 3.44690 3.52761 3.66996 4.33062 4.44879 4.47438 4.66567 4.82159 4.83546 4.99073 5.12475 5.17505 5.23115 5.33401 5.35180 5.47222 5.48613 5.56180 5.59031 5.66158 5.70790 5.71434 5.74410 5.80939 5.82204 5.89062 5.94172 5.95777 6.04579 6.12863 6.23358 6.31467 6.35984 6.40558 6.52663 6.54753 6.57412 6.68567 6.78727 6.80617 6.96237 7.00175 7.14289 7.20798 7.55426 7.65364 7.76554 8.00431 8.38687 8.65852 9.14793 9.77258 10.39477 11.25421 Molecular weight = 336.22amu Principal moments of inertia in cm(-1) A = 0.024613 B = 0.003041 C = 0.002933 Principal moments of inertia in units of 10**(-40)*gram-cm**2 A = 1137.344417 B = 9204.481330 C = 9545.541845 Net atomic charges, atomic populations, and dipole contributions calculated with CM2 Atom NO. Type Charge No. of electrons 1 C -0.151 4.151 2 C -0.065 4.065 3 H 0.085 0.915 4 C -0.136 4.136 5 C 0.515 3.485 6 O -0.554 6.554 7 N -0.590 5.590 8 C 0.094 3.906 9 C 0.151 3.849 10 H 0.076 0.924 11 C -0.121 4.121 12 C -0.136 4.136 13 C -0.112 4.112 14 C 0.179 3.821 15 H 0.113 0.887 16 C 0.092 3.908 17 N -0.579 5.579 18 C -0.251 4.251 19 H 0.066 0.934 20 C 0.004 3.996 21 N -0.895 5.895 22 Si 0.901 3.099 23 H -0.288 1.288 24 H -0.290 1.290 25 H -0.280 1.280 26 C 0.080 3.920 27 H 0.069 0.931 28 C -0.148 4.148 29 C -0.158 4.158 30 C 0.047 3.953 31 O -0.370 6.370 32 H 0.083 0.917 33 H 0.052 0.948 34 H 0.047 0.953 35 H 0.086 0.914 36 H 0.100 0.900 37 H 0.064 0.936 38 H 0.073 0.927 39 H 0.058 0.942 40 H 0.041 0.959 41 H 0.030 0.970 42 H 0.076 0.924 43 H 0.038 0.962 44 H 0.069 0.931 45 H 0.079 0.921 46 H 0.057 0.943 47 H 0.254 0.746 48 H 0.079 0.921 49 H 0.079 0.921 50 H 0.075 0.925 51 H 0.072 0.928 52 H 0.056 0.944 53 H 0.082 0.918 For the charges calculated by CM2: Dipole moment (debyes) X Y Z Total from point charges -4.374 -4.501 -20.053 21.012 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.208 4.208 2 C -0.084 4.084 3 H 0.104 0.896 4 C -0.176 4.176 5 C 0.304 3.696 6 O -0.434 6.434 7 N -0.323 5.323 8 C -0.030 4.030 9 C 0.043 3.957 10 H 0.094 0.906 11 C -0.160 4.160 12 C -0.174 4.174 13 C -0.151 4.151 14 C 0.071 3.929 15 H 0.131 0.869 16 C -0.031 4.031 17 N -0.300 5.300 18 C -0.380 4.380 19 H 0.084 0.916 20 C -0.195 4.195 21 N -0.538 5.538 22 Si 0.731 3.269 23 H -0.215 1.215 24 H -0.216 1.216 25 H -0.206 1.206 26 C 0.023 3.977 27 H 0.086 0.914 28 C -0.186 4.186 29 C -0.195 4.195 30 C -0.031 4.031 31 O -0.290 6.290 32 H 0.102 0.898 33 H 0.071 0.929 34 H 0.067 0.933 35 H 0.105 0.895 36 H 0.118 0.882 37 H 0.083 0.917 38 H 0.091 0.909 39 H 0.077 0.923 40 H 0.060 0.940 41 H 0.049 0.951 42 H 0.094 0.906 43 H 0.057 0.943 44 H 0.088 0.912 45 H 0.097 0.903 46 H 0.075 0.925 47 H 0.064 0.936 48 H 0.098 0.902 49 H 0.098 0.902 50 H 0.094 0.906 51 H 0.091 0.909 52 H 0.074 0.926 53 H 0.100 0.900 Dipole moment (debyes) X Y Z Total from point charges -4.535 -3.830 -19.740 20.614 hybrid contribution 0.120 -0.407 0.723 0.839 sum -4.415 -4.237 -19.017 19.977 Atomic orbital electron populations 1.21928 0.93856 1.01519 1.03481 1.20879 0.96227 1.00172 0.91111 0.89637 1.21745 1.03590 0.99668 0.92553 1.20306 0.78952 0.85703 0.84599 1.90616 1.41463 1.25057 1.86300 1.48169 1.53671 1.24477 1.05939 1.22483 0.94920 0.85357 1.00282 1.21288 0.91195 0.95387 0.87851 0.90607 1.21917 0.99182 0.96744 0.98202 1.22210 1.00344 0.96121 0.98758 1.21966 1.01023 0.96759 0.95402 1.21107 0.87811 0.86442 0.97495 0.86946 1.22492 0.99523 0.95903 0.85196 1.44402 1.63494 1.21829 1.00229 1.19106 1.44034 0.91209 0.83683 0.91589 1.25020 1.00072 0.96044 0.98322 1.69990 1.50580 1.03651 1.29620 0.86333 0.77842 0.77997 0.84769 1.21509 1.21614 1.20565 1.23289 0.81095 0.98399 0.94937 0.91352 1.23195 0.99037 1.01594 0.94756 1.22943 0.96249 1.01624 0.98717 1.23114 0.97854 1.00077 0.82010 1.88923 1.20849 1.93682 1.25505 0.89783 0.92896 0.93344 0.89538 0.88210 0.91726 0.90880 0.92302 0.94021 0.95097 0.90612 0.94293 0.91217 0.90262 0.92456 0.93599 0.90241 0.90211 0.90620 0.90880 0.92559 0.90018 Number of geometries 1 Number of calculations of the screened coulomb radii 1 The total number of SCF iterations 19. 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.15 -1.71 8.76 37.15 0.33 -1.38 16 2 C -0.06 -0.79 1.82 -90.63 -0.17 -0.96 16 3 H 0.09 1.22 7.96 -51.92 -0.41 0.81 16 4 C -0.14 -1.81 3.43 -27.89 -0.10 -1.90 16 5 C 0.51 8.71 4.77 -10.99 -0.05 8.66 16 6 O -0.55 -10.78 13.36 5.55 0.07 -10.71 16 7 N -0.59 -10.51 1.63 -171.93 -0.28 -10.79 16 8 C 0.09 1.55 5.86 -3.24 -0.02 1.53 16 9 C 0.15 2.94 3.58 -66.60 -0.24 2.70 16 10 H 0.08 1.44 8.01 -51.93 -0.42 1.02 16 11 C -0.12 -2.31 5.64 -26.53 -0.15 -2.46 16 12 C -0.14 -2.71 5.58 -26.10 -0.15 -2.86 16 13 C -0.11 -2.51 5.65 -25.92 -0.15 -2.65 16 14 C 0.18 4.08 2.62 -66.83 -0.18 3.91 16 15 H 0.11 2.93 5.89 -51.93 -0.31 2.62 16 16 C 0.09 1.86 5.89 -3.47 -0.02 1.84 16 17 N -0.58 -13.56 2.53 -155.13 -0.39 -13.95 16 18 C -0.25 -6.51 9.71 -15.06 -0.15 -6.66 16 19 H 0.07 1.68 7.89 -52.49 -0.41 1.27 16 20 C 0.00 0.12 4.98 -17.43 -0.09 0.03 16 21 N -0.89 -24.75 10.97 55.49 0.61 -24.14 16 22 Si 0.90 21.05 29.59 -169.99 -5.03 16.02 16 23 H -0.29 -6.79 7.11 56.52 0.40 -6.39 16 24 H -0.29 -6.88 7.11 56.52 0.40 -6.48 16 25 H -0.28 -6.46 7.11 56.52 0.40 -6.06 16 26 C 0.08 0.82 2.83 -25.55 -0.07 0.75 16 27 H 0.07 0.64 8.14 -51.93 -0.42 0.22 16 28 C -0.15 -1.24 5.69 -24.59 -0.14 -1.38 16 29 C -0.16 -1.13 7.00 -25.20 -0.18 -1.31 16 30 C 0.05 0.41 7.70 37.73 0.29 0.70 16 31 O -0.37 -4.43 10.60 -35.23 -0.37 -4.81 16 32 H 0.08 1.19 5.95 -51.93 -0.31 0.88 16 33 H 0.05 0.51 7.37 -51.93 -0.38 0.13 16 34 H 0.05 0.48 8.14 -51.93 -0.42 0.06 16 35 H 0.09 1.00 7.91 -51.93 -0.41 0.58 16 36 H 0.10 1.38 7.20 -51.93 -0.37 1.01 16 37 H 0.06 1.05 8.14 -51.93 -0.42 0.63 16 38 H 0.07 0.99 5.90 -51.93 -0.31 0.69 16 39 H 0.06 1.25 8.14 -51.93 -0.42 0.82 16 40 H 0.04 0.69 8.14 -51.93 -0.42 0.26 16 41 H 0.03 0.60 8.14 -51.93 -0.42 0.17 16 42 H 0.08 1.46 2.63 -51.93 -0.14 1.32 16 43 H 0.04 0.84 8.14 -51.93 -0.42 0.42 16 44 H 0.07 1.76 7.96 -51.93 -0.41 1.34 16 45 H 0.08 1.62 6.98 -51.93 -0.36 1.26 16 46 H 0.06 1.15 8.14 -51.93 -0.42 0.73 16 47 H 0.25 6.88 8.31 -40.82 -0.34 6.55 16 48 H 0.08 0.71 8.06 -51.93 -0.42 0.29 16 49 H 0.08 0.56 8.07 -51.93 -0.42 0.14 16 50 H 0.08 0.48 8.14 -51.93 -0.42 0.06 16 51 H 0.07 0.43 8.14 -51.93 -0.42 0.01 16 52 H 0.06 0.43 8.14 -51.93 -0.42 0.01 16 53 H 0.08 0.65 8.14 -51.93 -0.42 0.23 16 LS Contribution 385.31 15.07 5.81 5.81 Total: -1.00 -29.32 385.31 -10.09 -39.41 By element: Atomic # 1 Polarization: 13.90 SS G_CDS: -9.29 Total: 4.61 kcal Atomic # 6 Polarization: -0.23 SS G_CDS: -1.21 Total: -1.44 kcal Atomic # 7 Polarization: -48.82 SS G_CDS: -0.07 Total: -48.89 kcal Atomic # 8 Polarization: -15.22 SS G_CDS: -0.30 Total: -15.52 kcal Atomic # 14 Polarization: 21.05 SS G_CDS: -5.03 Total: 16.02 kcal Total LS contribution 5.81 Total: 5.81 kcal Total: -29.32 -10.09 -39.41 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. ZINC001039984158.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 -5.950 kcal (2) G-P(sol) polarization free energy of solvation -29.316 kcal (3) G-ENP(sol) elect.-nuc.-pol. free energy of system -35.266 kcal (4) G-CDS(sol) cavity-dispersion-solvent structure free energy -10.092 kcal (5) G-P-CDS(sol) = G-P(sol) + G-CDS(sol) = (2) + (4) -39.408 kcal (6) G-S(sol) free energy of system = (1) + (5) -45.358 kcal 1SCF WAS SPECIFIED, SO GEOMETRY OPTIMIZATION WAS NOT USED Total computer time = 1.20 seconds