Wall clock time and date at job start Sat Apr 11 2020 06:12:23 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.45197 * 1 3 3 C 1.34230 * 117.00112 * 2 1 4 4 O 1.20836 * 119.99673 * 359.97438 * 3 2 1 5 5 C 1.50693 * 120.00395 * 179.97438 * 3 2 1 6 6 N 1.46903 * 108.85474 * 316.41314 * 5 3 2 7 7 C 1.46899 * 110.99844 * 56.28376 * 6 5 3 8 8 C 1.53004 * 109.46768 * 170.97723 * 7 6 5 9 9 C 1.50701 * 109.46831 * 180.02562 * 8 7 6 10 10 O 1.21286 * 120.00302 * 359.97438 * 9 8 7 11 11 N 1.34774 * 120.00241 * 180.02562 * 9 8 7 12 12 C 1.46506 * 119.99602 * 359.97438 * 11 9 8 13 13 C 1.46494 * 120.00277 * 180.02562 * 11 9 8 14 14 C 1.50696 * 109.47338 * 89.99777 * 13 11 9 15 15 H 1.07995 * 120.00251 * 359.97438 * 14 13 11 16 16 C 1.37886 * 119.99995 * 180.02562 * 14 13 11 17 17 N 1.31510 * 119.99692 * 359.97438 * 16 14 13 18 18 Si 1.86796 * 120.00059 * 179.97438 * 16 14 13 19 19 H 1.48500 * 109.47137 * 89.99463 * 18 16 14 20 20 H 1.48493 * 109.47361 * 209.99727 * 18 16 14 21 21 H 1.48504 * 109.47035 * 329.99957 * 18 16 14 22 22 C 1.50480 * 108.79861 * 74.92418 * 5 3 2 23 23 C 1.54240 * 108.52759 * 214.87342 * 22 5 3 24 24 C 1.45544 * 112.40124 * 280.01904 * 23 22 5 25 25 C 1.54237 * 112.41058 * 61.93413 * 24 23 22 26 26 C 1.50472 * 108.52068 * 280.00572 * 25 24 23 27 27 C 1.54062 * 112.49514 * 94.15025 * 26 25 24 28 28 H 1.09005 * 109.47275 * 299.99779 * 1 2 3 29 29 H 1.08994 * 109.47145 * 60.00585 * 1 2 3 30 30 H 1.08997 * 109.47045 * 180.02562 * 1 2 3 31 31 H 1.00891 * 110.99653 * 292.33124 * 6 5 3 32 32 H 1.08996 * 109.47213 * 290.97631 * 7 6 5 33 33 H 1.09007 * 109.47216 * 50.97581 * 7 6 5 34 34 H 1.08996 * 109.47179 * 299.99930 * 8 7 6 35 35 H 1.08997 * 109.46891 * 60.00574 * 8 7 6 36 36 H 1.08996 * 109.46900 * 90.00161 * 12 11 9 37 37 H 1.09005 * 109.46916 * 210.00537 * 12 11 9 38 38 H 1.08999 * 109.46888 * 330.00204 * 12 11 9 39 39 H 1.09000 * 109.47042 * 209.99791 * 13 11 9 40 40 H 1.09001 * 109.47306 * 329.99421 * 13 11 9 41 41 H 0.97001 * 119.99865 * 180.02562 * 17 16 14 42 42 H 1.09003 * 109.62425 * 334.58648 * 22 5 3 43 43 H 1.09006 * 109.62605 * 95.16102 * 22 5 3 44 44 H 1.08994 * 108.89526 * 159.26402 * 23 22 5 45 45 H 1.09002 * 108.89434 * 40.76833 * 23 22 5 46 46 H 1.08999 * 108.89327 * 301.17912 * 24 23 22 47 47 H 1.09003 * 108.81649 * 182.64109 * 24 23 22 48 48 H 1.09007 * 109.62902 * 160.29326 * 25 24 23 49 49 H 1.08998 * 109.62815 * 39.71988 * 25 24 23 50 50 H 1.09002 * 108.87390 * 333.37794 * 26 25 24 51 51 H 1.08993 * 108.87698 * 214.92677 * 26 25 24 52 52 H 1.09001 * 108.81346 * 197.42316 * 27 26 25 53 53 H 1.09002 * 108.72384 * 79.08883 * 27 26 25 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.4520 0.0000 0.0000 3 6 2.0614 1.1960 0.0000 4 8 1.4032 2.2094 -0.0005 5 6 3.5662 1.2749 0.0006 6 7 4.0927 0.2944 0.9595 7 6 3.6127 -1.0582 0.6463 8 6 4.3518 -2.0759 1.5175 9 6 3.8599 -3.4636 1.1959 10 8 3.0022 -3.6229 0.3532 11 7 4.3722 -4.5276 1.8453 12 6 5.4078 -4.3352 2.8636 13 6 3.8944 -5.8766 1.5322 14 6 4.7353 -6.4638 0.4281 15 1 5.5444 -5.8894 0.0019 16 6 4.4710 -7.7340 -0.0387 17 7 3.4853 -8.4332 0.4800 18 14 5.5128 -8.4615 -1.4078 19 1 6.6674 -9.1867 -0.8192 20 1 4.6910 -9.4052 -2.2072 21 1 6.0086 -7.3710 -2.2854 22 6 4.0704 0.9300 -1.3746 23 6 5.4057 0.1715 -1.2311 24 6 6.5290 1.0659 -0.9936 25 6 6.3810 1.8595 0.3206 26 6 5.4502 3.0166 0.0781 27 6 3.9912 2.6868 0.4473 28 1 -0.3634 0.5138 0.8900 29 1 -0.3633 0.5137 -0.8900 30 1 -0.3633 -1.0276 0.0005 31 1 3.8567 0.5517 1.9060 32 1 2.5427 -1.1191 0.8449 33 1 3.8001 -1.2759 -0.4052 34 1 5.4218 -2.0151 1.3190 35 1 4.1643 -1.8583 2.5689 36 1 4.9383 -4.1969 3.8375 37 1 6.0552 -5.2116 2.8934 38 1 6.0002 -3.4541 2.6173 39 1 3.9708 -6.5043 2.4201 40 1 2.8542 -5.8280 1.2101 41 1 3.2996 -9.3270 0.1519 42 1 3.3428 0.2976 -1.8834 43 1 4.2284 1.8444 -1.9467 44 1 5.5910 -0.3930 -2.1450 45 1 5.3275 -0.5261 -0.3973 46 1 6.6025 1.7684 -1.8238 47 1 7.4447 0.4762 -0.9487 48 1 7.3558 2.2332 0.6345 49 1 5.9674 1.2127 1.0943 50 1 5.4940 3.2894 -0.9764 51 1 5.7814 3.8666 0.6744 52 1 3.3352 3.4176 -0.0257 53 1 3.8784 2.7621 1.5288 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 ZINC000179384606.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:12:23 Heat of formation + Delta-G solvation = -104.044097 kcal Electronic energy + Delta-G solvation = -32265.143692 eV Core-core repulsion = 28137.669743 eV Total energy + Delta-G solvation = -4127.473949 eV No. of doubly occupied orbitals = 66 Molecular weight (most abundant/longest-lived isotopes) = 340.220 amu Computer time = 0.70 seconds Orbital eigenvalues (eV) -41.48622 -39.26856 -38.77434 -36.72641 -35.29218 -34.69701 -32.87934 -32.59446 -31.84202 -29.29510 -27.98641 -27.55578 -25.80490 -25.18562 -24.94481 -23.30114 -21.97041 -21.43287 -20.32189 -19.35016 -18.28736 -17.95108 -17.29169 -16.83158 -16.50987 -16.21942 -15.85393 -15.27341 -14.98794 -14.85753 -14.54071 -14.42511 -14.10238 -13.73284 -13.52474 -13.41188 -13.30343 -13.04519 -12.78138 -12.32282 -12.29257 -12.27446 -12.20257 -12.10801 -11.81884 -11.61619 -11.44213 -11.26348 -11.08592 -10.96097 -10.90228 -10.88303 -10.79094 -10.56193 -10.48641 -10.38993 -10.14655 -10.14394 -10.03923 -9.60891 -8.85066 -8.55686 -8.42703 -7.60596 -6.03379 -4.07929 1.99554 2.54797 3.30510 3.37115 3.40692 3.50118 3.69203 4.22911 4.28353 4.45433 4.47677 4.58390 4.88021 4.97323 5.01323 5.05015 5.10762 5.16074 5.19127 5.23570 5.28798 5.32370 5.35190 5.45268 5.50647 5.59423 5.64652 5.69196 5.74245 5.86745 5.91276 5.94588 5.97434 5.99862 6.10432 6.14450 6.20112 6.23247 6.25211 6.33410 6.35693 6.52673 6.62938 6.80159 7.05932 7.10348 7.46082 7.57087 7.58005 7.93426 8.25116 8.46004 8.92182 9.14041 9.57269 11.04673 Molecular weight = 340.22amu Principal moments of inertia in cm(-1) A = 0.013288 B = 0.003088 C = 0.002854 Principal moments of inertia in units of 10**(-40)*gram-cm**2 A = 2106.642638 B = 9063.756153 C = 9808.429470 Net atomic charges, atomic populations, and dipole contributions calculated with CM2 Atom NO. Type Charge No. of electrons 1 C 0.024 3.976 2 O -0.334 6.334 3 C 0.424 3.576 4 O -0.513 6.513 5 C 0.143 3.857 6 N -0.658 5.658 7 C 0.044 3.956 8 C -0.135 4.135 9 C 0.504 3.496 10 O -0.541 6.541 11 N -0.604 5.604 12 C 0.084 3.916 13 C 0.245 3.755 14 C -0.525 4.525 15 H 0.057 0.943 16 C 0.060 3.940 17 N -0.896 5.896 18 Si 0.897 3.103 19 H -0.280 1.280 20 H -0.285 1.285 21 H -0.274 1.274 22 C -0.107 4.107 23 C -0.126 4.126 24 C -0.106 4.106 25 C -0.128 4.128 26 C -0.108 4.108 27 C -0.100 4.100 28 H 0.059 0.941 29 H 0.060 0.940 30 H 0.104 0.896 31 H 0.341 0.659 32 H 0.080 0.920 33 H 0.094 0.906 34 H 0.093 0.907 35 H 0.082 0.918 36 H 0.042 0.958 37 H 0.072 0.928 38 H 0.058 0.942 39 H 0.030 0.970 40 H 0.043 0.957 41 H 0.263 0.737 42 H 0.069 0.931 43 H 0.074 0.926 44 H 0.053 0.947 45 H 0.105 0.895 46 H 0.046 0.954 47 H 0.055 0.945 48 H 0.051 0.949 49 H 0.101 0.899 50 H 0.058 0.942 51 H 0.059 0.941 52 H 0.080 0.920 53 H 0.066 0.934 For the charges calculated by CM2: Dipole moment (debyes) X Y Z Total from point charges 3.015 22.293 2.842 22.675 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.070 4.070 2 O -0.250 6.250 3 C 0.264 3.736 4 O -0.399 6.399 5 C 0.049 3.951 6 N -0.290 5.290 7 C -0.083 4.083 8 C -0.175 4.175 9 C 0.290 3.710 10 O -0.419 6.419 11 N -0.336 5.336 12 C -0.060 4.060 13 C 0.126 3.874 14 C -0.564 4.564 15 H 0.075 0.925 16 C -0.142 4.142 17 N -0.535 5.535 18 Si 0.725 3.275 19 H -0.207 1.207 20 H -0.211 1.211 21 H -0.199 1.199 22 C -0.145 4.145 23 C -0.163 4.163 24 C -0.144 4.144 25 C -0.166 4.166 26 C -0.145 4.145 27 C -0.138 4.138 28 H 0.077 0.923 29 H 0.079 0.921 30 H 0.123 0.877 31 H 0.159 0.841 32 H 0.097 0.903 33 H 0.112 0.888 34 H 0.111 0.889 35 H 0.100 0.900 36 H 0.061 0.939 37 H 0.091 0.909 38 H 0.077 0.923 39 H 0.048 0.952 40 H 0.061 0.939 41 H 0.072 0.928 42 H 0.088 0.912 43 H 0.093 0.907 44 H 0.072 0.928 45 H 0.123 0.877 46 H 0.065 0.935 47 H 0.074 0.926 48 H 0.070 0.930 49 H 0.119 0.881 50 H 0.077 0.923 51 H 0.078 0.922 52 H 0.098 0.902 53 H 0.085 0.915 Dipole moment (debyes) X Y Z Total from point charges 2.809 22.207 2.058 22.479 hybrid contribution -0.601 -0.404 -0.415 0.835 sum 2.208 21.803 1.643 21.976 Atomic orbital electron populations 1.23336 0.77553 1.04116 1.01962 1.87025 1.22917 1.32006 1.83066 1.25829 0.93681 0.81854 0.72208 1.90618 1.67490 1.36607 1.45222 1.20888 0.82783 0.94753 0.96646 1.60989 1.60210 1.02194 1.05563 1.22444 0.99041 0.85162 1.01641 1.21607 1.02023 0.92984 1.00930 1.20264 0.82756 0.86020 0.81939 1.90675 1.32328 1.86212 1.32696 1.48155 1.38676 1.08860 1.37940 1.21629 0.91456 1.02149 0.90801 1.19260 0.96418 0.75741 0.95981 1.21444 1.15858 0.99795 1.19261 0.92460 1.24940 0.95379 0.96953 0.96886 1.69933 1.29445 1.11512 1.42589 0.86506 0.80015 0.78933 0.82030 1.20663 1.21120 1.19865 1.21447 0.98886 1.01396 0.92812 1.20627 0.92618 0.99543 1.03536 1.20360 0.98336 0.98470 0.97194 1.21547 0.98232 0.97400 0.99397 1.21243 0.93371 0.98549 1.01372 1.21947 0.97669 0.92600 1.01550 0.92253 0.92135 0.87737 0.84055 0.90276 0.88793 0.88887 0.89963 0.93923 0.90927 0.92288 0.95176 0.93913 0.92759 0.91236 0.90710 0.92783 0.87725 0.93497 0.92601 0.93038 0.88121 0.92332 0.92214 0.90174 0.91491 Number of geometries 1 Number of calculations of the screened coulomb radii 1 The total number of SCF iterations 12. 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.27 10.56 101.05 1.07 1.33 16 2 O -0.33 -4.41 6.59 -39.26 -0.26 -4.67 16 3 C 0.42 5.34 6.43 36.00 0.23 5.58 16 4 O -0.51 -7.43 15.30 -18.76 -0.29 -7.72 16 5 C 0.14 1.46 0.53 -133.46 -0.07 1.39 16 6 N -0.66 -6.37 2.37 -101.86 -0.24 -6.61 16 7 C 0.04 0.57 2.67 -3.82 -0.01 0.56 16 8 C -0.13 -1.81 4.50 -27.88 -0.13 -1.94 16 9 C 0.50 10.00 7.67 -10.99 -0.08 9.91 16 10 O -0.54 -13.11 15.56 5.55 0.09 -13.02 16 11 N -0.60 -12.48 2.86 -182.95 -0.52 -13.00 16 12 C 0.08 1.28 10.24 59.85 0.61 1.89 16 13 C 0.25 6.37 5.16 -5.20 -0.03 6.35 16 14 C -0.52 -14.81 9.39 -34.44 -0.32 -15.14 16 15 H 0.06 1.61 7.81 -52.49 -0.41 1.20 16 16 C 0.06 1.72 5.46 -17.82 -0.10 1.62 16 17 N -0.90 -26.59 13.29 55.43 0.74 -25.85 16 18 Si 0.90 21.80 29.54 -169.99 -5.02 16.78 16 19 H -0.28 -6.61 7.11 56.52 0.40 -6.20 16 20 H -0.29 -6.89 7.11 56.52 0.40 -6.49 16 21 H -0.27 -6.63 7.11 56.52 0.40 -6.22 16 22 C -0.11 -1.13 4.26 -27.34 -0.12 -1.25 16 23 C -0.13 -1.33 4.46 -29.60 -0.13 -1.46 16 24 C -0.11 -0.97 5.26 -29.60 -0.16 -1.13 16 25 C -0.13 -1.06 5.61 -27.35 -0.15 -1.21 16 26 C -0.11 -0.84 5.31 -27.44 -0.15 -0.99 16 27 C -0.10 -0.89 4.68 -25.62 -0.12 -1.01 16 28 H 0.06 0.61 8.13 -51.93 -0.42 0.19 16 29 H 0.06 0.64 8.13 -51.93 -0.42 0.22 16 30 H 0.10 1.06 8.14 -51.93 -0.42 0.64 16 31 H 0.34 2.87 8.53 -39.30 -0.34 2.54 16 32 H 0.08 1.16 5.08 -52.19 -0.26 0.90 16 33 H 0.09 1.39 5.45 -51.92 -0.28 1.11 16 34 H 0.09 1.14 7.14 -51.93 -0.37 0.77 16 35 H 0.08 0.92 8.11 -51.93 -0.42 0.50 16 36 H 0.04 0.56 8.14 -51.93 -0.42 0.14 16 37 H 0.07 1.16 8.07 -51.93 -0.42 0.74 16 38 H 0.06 0.72 6.44 -51.93 -0.33 0.39 16 39 H 0.03 0.79 8.07 -51.93 -0.42 0.37 16 40 H 0.04 1.25 7.42 -51.93 -0.39 0.87 16 41 H 0.26 7.36 8.31 -40.82 -0.34 7.02 16 42 H 0.07 0.82 7.49 -51.93 -0.39 0.43 16 43 H 0.07 0.72 6.98 -51.93 -0.36 0.36 16 44 H 0.05 0.60 8.14 -51.93 -0.42 0.18 16 45 H 0.10 1.21 3.88 -51.93 -0.20 1.01 16 46 H 0.05 0.41 7.47 -51.93 -0.39 0.02 16 47 H 0.06 0.52 8.14 -51.93 -0.42 0.09 16 48 H 0.05 0.40 8.14 -51.93 -0.42 -0.03 16 49 H 0.10 0.87 5.57 -51.93 -0.29 0.58 16 50 H 0.06 0.46 6.63 -51.93 -0.34 0.12 16 51 H 0.06 0.41 8.14 -51.93 -0.42 -0.01 16 52 H 0.08 0.80 7.12 -51.93 -0.37 0.43 16 53 H 0.07 0.55 7.85 -51.93 -0.41 0.14 16 LS Contribution 397.56 15.07 5.99 5.99 Total: -1.00 -33.51 397.56 -8.08 -41.58 By element: Atomic # 1 Polarization: 10.91 SS G_CDS: -8.91 Total: 2.00 kcal Atomic # 6 Polarization: 4.18 SS G_CDS: 0.35 Total: 4.53 kcal Atomic # 7 Polarization: -45.44 SS G_CDS: -0.03 Total: -45.46 kcal Atomic # 8 Polarization: -24.96 SS G_CDS: -0.46 Total: -25.42 kcal Atomic # 14 Polarization: 21.80 SS G_CDS: -5.02 Total: 16.78 kcal Total LS contribution 5.99 Total: 5.99 kcal Total: -33.51 -8.08 -41.58 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. ZINC000179384606.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.460 kcal (2) G-P(sol) polarization free energy of solvation -33.506 kcal (3) G-ENP(sol) elect.-nuc.-pol. free energy of system -95.966 kcal (4) G-CDS(sol) cavity-dispersion-solvent structure free energy -8.078 kcal (5) G-P-CDS(sol) = G-P(sol) + G-CDS(sol) = (2) + (4) -41.584 kcal (6) G-S(sol) free energy of system = (1) + (5) -104.044 kcal 1SCF WAS SPECIFIED, SO GEOMETRY OPTIMIZATION WAS NOT USED Total computer time = 0.70 seconds