Wall clock time and date at job start Tue Mar 31 2020 09:00:33 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 N 2 2 C 1.31521 * 1 3 3 Si 1.86797 * 120.00291 * 2 1 4 4 H 1.48504 * 109.46787 * 270.00064 * 3 2 1 5 5 H 1.48503 * 109.47208 * 29.99764 * 3 2 1 6 6 H 1.48497 * 109.47218 * 149.99627 * 3 2 1 7 7 C 1.37874 * 119.99341 * 179.97438 * 2 1 3 8 8 H 1.08008 * 119.99702 * 180.02562 * 7 2 1 9 9 C 1.50695 * 120.00671 * 0.02562 * 7 2 1 10 10 C 1.53123 * 109.57273 * 299.71759 * 9 7 2 11 11 C 1.53035 * 109.27774 * 182.41108 * 10 9 7 12 12 N 1.46846 * 109.52744 * 59.17321 * 11 10 9 13 13 C 1.46904 * 110.99942 * 174.16069 * 12 11 10 14 14 H 1.09000 * 109.47065 * 55.00354 * 13 12 11 15 15 C 1.52994 * 109.47136 * 295.00288 * 13 12 11 16 16 C 1.53006 * 109.47037 * 180.02562 * 15 13 12 17 17 C 1.52998 * 109.47093 * 59.99574 * 16 15 13 18 18 C 1.52997 * 109.47491 * 300.00494 * 17 16 15 19 19 C 1.52998 * 109.47491 * 179.97438 * 18 17 16 20 20 C 1.52994 * 109.47527 * 179.97438 * 19 18 17 21 21 C 1.53004 * 109.47189 * 300.00061 * 20 19 18 22 22 C 1.53008 * 109.46956 * 59.99608 * 21 20 19 23 23 C 1.52991 * 109.47128 * 300.00208 * 22 21 20 24 24 C 1.52997 * 109.47335 * 59.99515 * 18 17 16 25 25 C 1.46844 * 111.19783 * 298.26158 * 12 11 10 26 26 C 1.53032 * 109.52527 * 61.73901 * 25 12 11 27 27 H 0.97007 * 119.99481 * 359.97438 * 1 2 3 28 28 H 1.08998 * 109.52068 * 59.86745 * 9 7 2 29 29 H 1.08995 * 109.51026 * 62.46381 * 10 9 7 30 30 H 1.09011 * 109.49940 * 302.34346 * 10 9 7 31 31 H 1.09008 * 109.45841 * 299.15478 * 11 10 9 32 32 H 1.09003 * 109.46290 * 179.19461 * 11 10 9 33 33 H 1.09000 * 109.47198 * 299.99985 * 15 13 12 34 34 H 1.09000 * 109.47573 * 60.00090 * 15 13 12 35 35 H 1.09000 * 109.47027 * 299.99982 * 16 15 13 36 36 H 1.08997 * 109.46848 * 179.97438 * 16 15 13 37 37 H 1.09000 * 109.46967 * 179.97438 * 17 16 15 38 38 H 1.09000 * 109.46851 * 60.00300 * 17 16 15 39 39 H 1.09000 * 109.46851 * 299.99806 * 19 18 17 40 40 H 1.09000 * 109.46967 * 59.99467 * 19 18 17 41 41 H 1.09000 * 109.47573 * 179.97438 * 20 19 18 42 42 H 1.09004 * 109.46926 * 59.99466 * 20 19 18 43 43 H 1.09000 * 109.46982 * 299.99966 * 21 20 19 44 44 H 1.08998 * 109.47220 * 179.97438 * 21 20 19 45 45 H 1.09000 * 109.46871 * 180.02562 * 22 21 20 46 46 H 1.09000 * 109.46983 * 60.00358 * 22 21 20 47 47 H 1.09000 * 109.47310 * 299.99949 * 23 22 21 48 48 H 1.09000 * 109.47209 * 180.02562 * 23 22 21 49 49 H 1.09001 * 109.47275 * 180.02562 * 24 18 17 50 50 H 1.09000 * 109.46913 * 60.00252 * 24 18 17 51 51 H 1.09000 * 109.46271 * 181.75698 * 25 12 11 52 52 H 1.08998 * 109.46026 * 301.72258 * 25 12 11 53 53 H 1.09000 * 109.50025 * 60.76400 * 26 25 12 54 54 H 1.08996 * 109.50621 * 180.87820 * 26 25 12 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.3152 0.0000 0.0000 3 14 2.2493 1.6177 0.0000 4 1 2.4967 2.0463 -1.4001 5 1 1.4467 2.6527 0.7000 6 1 3.5468 1.4401 0.7001 7 6 2.0044 -1.1941 0.0005 8 1 3.0845 -1.1942 0.0010 9 6 1.2510 -2.4992 0.0005 10 6 0.3746 -2.5855 1.2532 11 6 -0.3411 -3.9380 1.2769 12 7 0.6526 -5.0192 1.2773 13 6 0.0061 -6.3284 1.4392 14 1 -0.7535 -6.4563 0.6681 15 6 -0.6499 -6.4060 2.8192 16 6 -1.3228 -7.7697 2.9881 17 6 -0.2733 -8.8756 2.8592 18 6 0.3828 -8.7979 1.4793 19 6 1.4327 -9.9034 1.3507 20 6 2.0883 -9.8261 -0.0295 21 6 1.0219 -10.0073 -1.1116 22 6 -0.0280 -8.9017 -0.9830 23 6 -0.6836 -8.9791 0.3971 24 6 1.0555 -7.4342 1.3104 25 6 1.4749 -4.9851 0.0611 26 6 2.2448 -3.6639 0.0012 27 1 -0.4850 0.8401 0.0004 28 1 0.6226 -2.5569 -0.8882 29 1 -0.3633 -1.7835 1.2372 30 1 0.9990 -2.4887 2.1415 31 1 -0.9759 -4.0284 0.3954 32 1 -0.9548 -4.0085 2.1749 33 1 0.1097 -6.2780 3.5903 34 1 -1.3973 -5.6179 2.9111 35 1 -2.0823 -7.8981 2.2169 36 1 -1.7897 -7.8252 3.9714 37 1 -0.7526 -9.8471 2.9800 38 1 0.4863 -8.7472 3.6304 39 1 2.1922 -9.7746 2.1218 40 1 0.9537 -10.8751 1.4712 41 1 2.8366 -10.6134 -0.1209 42 1 2.5671 -8.8544 -0.1504 43 1 0.5430 -10.9789 -0.9907 44 1 1.4889 -9.9518 -2.0949 45 1 -0.7880 -9.0310 -1.7536 46 1 0.4508 -7.9300 -1.1039 47 1 -1.1625 -9.9507 0.5180 48 1 -1.4319 -8.1918 0.4886 49 1 1.5226 -7.3787 0.3271 50 1 1.8151 -7.3059 2.0816 51 1 2.1801 -5.8161 0.0761 52 1 0.8316 -5.0707 -0.8146 53 1 2.8986 -3.5840 0.8697 54 1 2.8430 -3.6323 -0.9094 RHF calculation, no. of doubly occupied orbitals= 60 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=GENORG IOFR=1.4345 ALPHA=0.00 BETA=0.00 GAMMA=38.93 & DIELEC=2.06 FACARB=0.00 FEHALO=0.00 DEV ZINC001173344580.mol2 54 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 09:00:33 Heat of formation + Delta-G solvation = -40.961305 kcal Electronic energy + Delta-G solvation = -27507.799017 eV Core-core repulsion = 24248.995411 eV Total energy + Delta-G solvation = -3258.803606 eV No. of doubly occupied orbitals = 60 Molecular weight (most abundant/longest-lived isotopes) = 305.241 amu Computer time = 0.68 seconds Orbital eigenvalues (eV) -40.69911 -38.03867 -35.79304 -33.65334 -31.61764 -31.10641 -30.72750 -29.77928 -28.61887 -25.94837 -25.35729 -23.16230 -22.77542 -22.04552 -21.32812 -20.62873 -19.74655 -18.57942 -17.99758 -16.54169 -16.03209 -15.28312 -14.84737 -14.51587 -14.45918 -14.21852 -14.02954 -13.96432 -13.39724 -13.17601 -12.92319 -12.88388 -12.52383 -12.43988 -12.27455 -12.21786 -11.94625 -11.66158 -11.40573 -11.19217 -11.14331 -10.97118 -10.91851 -10.81842 -10.66814 -10.63340 -10.40026 -10.34771 -10.12406 -10.04946 -9.88072 -9.70191 -9.65926 -9.32486 -8.82960 -8.62375 -8.24291 -7.60376 -5.92761 -3.97828 3.33836 3.40496 3.44285 4.32864 4.48911 4.50073 4.59507 4.68506 4.98126 5.02233 5.04995 5.22196 5.23469 5.31561 5.37255 5.42696 5.46188 5.49717 5.53544 5.54388 5.62483 5.70299 5.71085 5.81882 5.89480 5.93962 5.94645 5.98574 6.01359 6.06982 6.09546 6.11810 6.18170 6.22552 6.26113 6.30240 6.34676 6.36181 6.41103 6.42659 6.54246 6.63776 6.70946 6.73430 6.93497 7.15984 7.25280 7.52403 7.71181 7.82495 7.92017 8.11656 8.22463 8.53965 9.05954 9.66272 11.13327 Molecular weight = 305.24amu Principal moments of inertia in cm(-1) A = 0.023678 B = 0.003491 C = 0.003253 Principal moments of inertia in units of 10**(-40)*gram-cm**2 A = 1182.268408 B = 8017.921607 C = 8604.193639 Net atomic charges, atomic populations, and dipole contributions calculated with CM2 Atom NO. Type Charge No. of electrons 1 N -0.892 5.892 2 C 0.060 3.940 3 Si 0.889 3.111 4 H -0.282 1.282 5 H -0.287 1.287 6 H -0.275 1.275 7 C -0.509 4.509 8 H 0.053 0.947 9 C 0.048 3.952 10 C -0.103 4.103 11 C 0.057 3.943 12 N -0.531 5.531 13 C 0.085 3.915 14 H 0.035 0.965 15 C -0.111 4.111 16 C -0.119 4.119 17 C -0.107 4.107 18 C -0.053 4.053 19 C -0.106 4.106 20 C -0.123 4.123 21 C -0.116 4.116 22 C -0.123 4.123 23 C -0.110 4.110 24 C -0.105 4.105 25 C 0.058 3.942 26 C -0.106 4.106 27 H 0.258 0.742 28 H 0.019 0.981 29 H 0.090 0.910 30 H 0.042 0.958 31 H 0.019 0.981 32 H 0.057 0.943 33 H 0.068 0.932 34 H 0.061 0.939 35 H 0.061 0.939 36 H 0.059 0.941 37 H 0.056 0.944 38 H 0.061 0.939 39 H 0.060 0.940 40 H 0.058 0.942 41 H 0.057 0.943 42 H 0.066 0.934 43 H 0.057 0.943 44 H 0.057 0.943 45 H 0.057 0.943 46 H 0.066 0.934 47 H 0.056 0.944 48 H 0.067 0.933 49 H 0.071 0.929 50 H 0.067 0.933 51 H 0.061 0.939 52 H 0.019 0.981 53 H 0.055 0.945 54 H 0.047 0.953 For the charges calculated by CM2: Dipole moment (debyes) X Y Z Total from point charges -1.558 -21.129 2.840 21.376 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.532 5.532 2 C -0.140 4.140 3 Si 0.717 3.283 4 H -0.208 1.208 5 H -0.213 1.213 6 H -0.200 1.200 7 C -0.547 4.547 8 H 0.072 0.928 9 C 0.030 3.970 10 C -0.141 4.141 11 C -0.070 4.070 12 N -0.256 5.256 13 C -0.023 4.023 14 H 0.053 0.947 15 C -0.149 4.149 16 C -0.157 4.157 17 C -0.144 4.144 18 C -0.054 4.054 19 C -0.143 4.143 20 C -0.160 4.160 21 C -0.154 4.154 22 C -0.161 4.161 23 C -0.148 4.148 24 C -0.142 4.142 25 C -0.069 4.069 26 C -0.144 4.144 27 H 0.068 0.932 28 H 0.037 0.963 29 H 0.109 0.891 30 H 0.061 0.939 31 H 0.037 0.963 32 H 0.075 0.925 33 H 0.086 0.914 34 H 0.080 0.920 35 H 0.079 0.921 36 H 0.077 0.923 37 H 0.075 0.925 38 H 0.080 0.920 39 H 0.079 0.921 40 H 0.077 0.923 41 H 0.076 0.924 42 H 0.085 0.915 43 H 0.076 0.924 44 H 0.076 0.924 45 H 0.075 0.925 46 H 0.085 0.915 47 H 0.075 0.925 48 H 0.086 0.914 49 H 0.089 0.911 50 H 0.085 0.915 51 H 0.079 0.921 52 H 0.038 0.962 53 H 0.074 0.926 54 H 0.066 0.934 Dipole moment (debyes) X Y Z Total from point charges -1.396 -20.914 3.237 21.209 hybrid contribution -0.155 0.847 -0.721 1.123 sum -1.550 -20.067 2.516 20.284 Atomic orbital electron populations 1.69890 1.05850 1.26191 1.51241 1.24965 0.94824 0.98945 0.95280 0.86649 0.79946 0.83746 0.77972 1.20847 1.21322 1.19967 1.21226 0.92858 0.90744 1.49909 0.92827 1.18670 0.92699 0.91705 0.93944 1.21640 1.00101 0.93072 0.99290 1.22288 0.93192 0.93064 0.98500 1.62012 1.23641 1.01854 1.38053 1.21598 0.94883 0.89586 0.96281 0.94742 1.21581 0.99741 0.96186 0.97350 1.21564 0.97946 0.94956 1.01205 1.21578 0.98198 0.98179 0.96460 1.20510 0.95040 0.95113 0.94701 1.21553 0.97767 0.98049 0.96950 1.21579 0.97638 1.01975 0.94798 1.21506 0.96719 0.98704 0.98467 1.21585 0.98115 0.99410 0.96976 1.21668 0.97637 1.01412 0.94038 1.21809 0.97834 0.92172 1.02403 1.22166 0.95899 0.99121 0.89748 1.21672 0.98690 0.93961 1.00072 0.93181 0.96318 0.89146 0.93935 0.96320 0.92509 0.91361 0.91996 0.92062 0.92261 0.92509 0.91977 0.92107 0.92339 0.92446 0.91494 0.92429 0.92385 0.92454 0.91514 0.92504 0.91439 0.91075 0.91474 0.92077 0.96233 0.92596 0.93447 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 N -0.89 -25.72 11.37 55.43 0.63 -25.09 16 2 C 0.06 1.69 5.31 -17.82 -0.09 1.60 16 3 Si 0.89 21.55 29.54 -169.99 -5.02 16.53 16 4 H -0.28 -6.79 7.11 56.52 0.40 -6.39 16 5 H -0.29 -6.92 7.11 56.52 0.40 -6.52 16 6 H -0.27 -6.59 7.11 56.52 0.40 -6.19 16 7 C -0.51 -14.20 8.56 -34.44 -0.29 -14.49 16 8 H 0.05 1.48 7.75 -52.48 -0.41 1.08 16 9 C 0.05 1.22 2.16 -91.65 -0.20 1.02 16 10 C -0.10 -2.47 5.12 -26.65 -0.14 -2.60 16 11 C 0.06 1.13 5.29 -3.84 -0.02 1.11 16 12 N -0.53 -9.47 4.52 -226.70 -1.03 -10.49 16 13 C 0.08 1.22 1.68 -67.71 -0.11 1.11 16 14 H 0.03 0.53 6.42 -51.93 -0.33 0.19 16 15 C -0.11 -1.45 5.11 -26.73 -0.14 -1.59 16 16 C -0.12 -1.25 5.65 -26.73 -0.15 -1.41 16 17 C -0.11 -1.04 4.73 -26.73 -0.13 -1.17 16 18 C -0.05 -0.54 0.51 -154.51 -0.08 -0.62 16 19 C -0.11 -0.97 4.73 -26.74 -0.13 -1.10 16 20 C -0.12 -1.13 5.65 -26.73 -0.15 -1.28 16 21 C -0.12 -1.06 5.92 -26.72 -0.16 -1.22 16 22 C -0.12 -1.26 5.65 -26.73 -0.15 -1.41 16 23 C -0.11 -1.12 4.04 -26.73 -0.11 -1.22 16 24 C -0.10 -1.30 2.78 -26.73 -0.07 -1.38 16 25 C 0.06 1.06 5.00 -3.84 -0.02 1.04 16 26 C -0.11 -2.27 5.32 -26.66 -0.14 -2.42 16 27 H 0.26 7.12 8.31 -40.82 -0.34 6.78 16 28 H 0.02 0.52 8.14 -51.93 -0.42 0.10 16 29 H 0.09 2.41 6.17 -51.93 -0.32 2.09 16 30 H 0.04 1.05 8.14 -51.92 -0.42 0.63 16 31 H 0.02 0.39 8.14 -51.92 -0.42 -0.04 16 32 H 0.06 1.05 6.18 -51.93 -0.32 0.73 16 33 H 0.07 0.94 8.14 -51.93 -0.42 0.52 16 34 H 0.06 0.86 6.40 -51.93 -0.33 0.53 16 35 H 0.06 0.64 6.43 -51.93 -0.33 0.31 16 36 H 0.06 0.57 8.14 -51.93 -0.42 0.15 16 37 H 0.06 0.50 8.14 -51.93 -0.42 0.08 16 38 H 0.06 0.62 8.14 -51.93 -0.42 0.20 16 39 H 0.06 0.56 8.14 -51.93 -0.42 0.14 16 40 H 0.06 0.50 8.14 -51.93 -0.42 0.08 16 41 H 0.06 0.48 8.14 -51.93 -0.42 0.05 16 42 H 0.07 0.68 6.43 -51.93 -0.33 0.35 16 43 H 0.06 0.49 8.14 -51.93 -0.42 0.07 16 44 H 0.06 0.51 8.14 -51.93 -0.42 0.09 16 45 H 0.06 0.56 8.14 -51.93 -0.42 0.14 16 46 H 0.07 0.78 6.43 -51.93 -0.33 0.45 16 47 H 0.06 0.53 8.14 -51.93 -0.42 0.11 16 48 H 0.07 0.75 4.81 -51.93 -0.25 0.50 16 49 H 0.07 0.92 2.85 -51.93 -0.15 0.77 16 50 H 0.07 0.90 8.14 -51.93 -0.42 0.47 16 51 H 0.06 0.99 5.83 -51.93 -0.30 0.69 16 52 H 0.02 0.37 8.14 -51.93 -0.42 -0.06 16 53 H 0.06 1.23 8.14 -51.93 -0.42 0.80 16 54 H 0.05 0.98 8.14 -51.93 -0.42 0.56 16 LS Contribution 370.55 15.07 5.58 5.58 Total: -1.00 -27.74 370.55 -12.27 -40.02 By element: Atomic # 1 Polarization: 9.62 SS G_CDS: -10.16 Total: -0.55 kcal Atomic # 6 Polarization: -23.73 SS G_CDS: -2.28 Total: -26.01 kcal Atomic # 7 Polarization: -35.18 SS G_CDS: -0.39 Total: -35.58 kcal Atomic # 14 Polarization: 21.55 SS G_CDS: -5.02 Total: 16.53 kcal Total LS contribution 5.58 Total: 5.58 kcal Total: -27.74 -12.27 -40.02 kcal The number of atoms in the molecule is 54 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. ZINC001173344580.mol2 54 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 -0.942 kcal (2) G-P(sol) polarization free energy of solvation -27.745 kcal (3) G-ENP(sol) elect.-nuc.-pol. free energy of system -28.687 kcal (4) G-CDS(sol) cavity-dispersion-solvent structure free energy -12.275 kcal (5) G-P-CDS(sol) = G-P(sol) + G-CDS(sol) = (2) + (4) -40.020 kcal (6) G-S(sol) free energy of system = (1) + (5) -40.961 kcal 1SCF WAS SPECIFIED, SO GEOMETRY OPTIMIZATION WAS NOT USED Total computer time = 0.68 seconds