Wall clock time and date at job start Sat Apr 11 2020 06:11:34 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.42898 * 1 3 3 C 1.42904 * 113.99691 * 2 1 4 4 C 1.53005 * 109.46705 * 180.02562 * 3 2 1 5 5 N 1.46498 * 109.46830 * 65.00506 * 4 3 2 6 6 C 1.36934 * 119.99978 * 269.99155 * 5 4 3 7 7 H 1.08000 * 120.00189 * 145.91106 * 6 5 4 8 8 C 1.38813 * 120.00024 * 325.91672 * 6 5 4 9 9 N 1.31622 * 120.00102 * 157.02523 * 8 6 5 10 10 Si 1.86801 * 120.00143 * 337.02391 * 8 6 5 11 11 H 1.48500 * 109.47175 * 95.41430 * 10 8 6 12 12 H 1.48494 * 109.47443 * 215.41774 * 10 8 6 13 13 H 1.48504 * 109.47068 * 335.41929 * 10 8 6 14 14 C 1.46499 * 120.00300 * 90.00054 * 5 4 3 15 15 C 1.53038 * 109.45951 * 237.96627 * 14 5 4 16 16 C 1.53191 * 109.31237 * 178.59955 * 15 14 5 17 17 N 1.46921 * 108.76826 * 54.65278 * 16 15 14 18 18 C 1.34775 * 120.63038 * 126.40213 * 17 16 15 19 19 O 1.21278 * 120.00254 * 5.02035 * 18 17 16 20 20 C 1.50703 * 119.99912 * 185.01494 * 18 17 16 21 21 H 1.08996 * 109.52281 * 46.71995 * 20 18 17 22 22 C 1.53183 * 109.52016 * 166.91005 * 20 18 17 23 23 C 1.53179 * 109.03356 * 176.88783 * 22 20 18 24 24 C 1.53091 * 109.16458 * 303.01330 * 23 22 20 25 25 O 1.42903 * 109.41070 * 57.61110 * 24 23 22 26 26 C 1.42906 * 114.10112 * 298.84012 * 25 24 23 27 27 C 1.46924 * 118.74026 * 306.38831 * 17 16 15 28 28 C 1.53040 * 109.46544 * 117.90898 * 14 5 4 29 29 H 1.08996 * 109.47159 * 180.02562 * 1 2 3 30 30 H 1.09000 * 109.47516 * 59.99366 * 1 2 3 31 31 H 1.09001 * 109.46874 * 299.99245 * 1 2 3 32 32 H 1.08999 * 109.47117 * 300.00562 * 3 2 1 33 33 H 1.08997 * 109.46956 * 60.00713 * 3 2 1 34 34 H 1.08998 * 109.46622 * 184.99862 * 4 3 2 35 35 H 1.08995 * 109.47310 * 304.99732 * 4 3 2 36 36 H 0.97004 * 120.00036 * 179.97438 * 9 8 6 37 37 H 1.08999 * 109.45749 * 357.93631 * 14 5 4 38 38 H 1.09000 * 109.49491 * 298.56511 * 15 14 5 39 39 H 1.09000 * 109.49781 * 58.66137 * 15 14 5 40 40 H 1.09001 * 109.60041 * 294.86533 * 16 15 14 41 41 H 1.09003 * 109.58445 * 174.43614 * 16 15 14 42 42 H 1.09001 * 109.54154 * 296.75113 * 22 20 18 43 43 H 1.08995 * 109.55151 * 57.01592 * 22 20 18 44 44 H 1.08999 * 109.51934 * 183.10805 * 23 22 20 45 45 H 1.09001 * 109.57829 * 62.95712 * 23 22 20 46 46 H 1.08993 * 109.48257 * 297.63049 * 24 23 22 47 47 H 1.08995 * 109.49303 * 177.59474 * 24 23 22 48 48 H 1.08993 * 109.47929 * 301.18210 * 26 25 24 49 49 H 1.08997 * 109.48587 * 181.14098 * 26 25 24 50 50 H 1.09003 * 109.58367 * 293.81410 * 27 17 16 51 51 H 1.09001 * 109.58653 * 173.39431 * 27 17 16 52 52 H 1.08997 * 109.49625 * 301.34164 * 28 14 5 53 53 H 1.09006 * 109.49453 * 61.43215 * 28 14 5 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.0102 1.3055 0.0000 4 6 3.5354 1.1847 -0.0006 5 7 3.9755 0.5646 -1.2528 6 6 4.0967 -0.7968 -1.3370 7 1 3.8751 -1.3023 -2.2653 8 6 4.5044 -1.5272 -0.2292 9 7 4.2027 -2.8050 -0.1360 10 14 5.4812 -0.6967 1.1293 11 1 4.5580 -0.2629 2.2086 12 1 6.4734 -1.6535 1.6816 13 1 6.1886 0.4868 0.5778 14 6 4.2862 1.4008 -2.4148 15 6 3.4027 0.9849 -3.5932 16 6 3.7055 1.8843 -4.7957 17 7 5.1471 1.8312 -5.0745 18 6 5.5973 1.5193 -6.3059 19 8 4.8126 1.3685 -7.2183 20 6 7.0758 1.3581 -6.5492 21 1 7.5082 0.7363 -5.7654 22 6 7.3010 0.6962 -7.9122 23 6 8.8063 0.6086 -8.1823 24 6 9.4123 2.0126 -8.1111 25 8 9.1401 2.5829 -6.8293 26 6 7.7486 2.7332 -6.5411 27 6 6.0886 2.1250 -3.9855 28 6 5.7588 1.2232 -2.7919 29 1 -0.3633 -1.0276 0.0005 30 1 -0.3634 0.5139 -0.8899 31 1 -0.3633 0.5137 0.8901 32 1 1.6884 1.8465 -0.8899 33 1 1.6886 1.8463 0.8900 34 1 3.9784 2.1765 0.0889 35 1 3.8514 0.5683 0.8409 36 1 4.4879 -3.3155 0.6380 37 1 4.0975 2.4462 -2.1707 38 1 3.6089 -0.0532 -3.8538 39 1 2.3537 1.0894 -3.3160 40 1 3.4150 2.9098 -4.5677 41 1 3.1519 1.5303 -5.6653 42 1 6.8727 -0.3061 -7.9061 43 1 6.8251 1.2927 -8.6905 44 1 8.9745 0.1892 -9.1742 45 1 9.2753 -0.0291 -7.4330 46 1 8.9727 2.6381 -8.8879 47 1 10.4902 1.9515 -8.2602 48 1 7.2872 3.3681 -7.2974 49 1 7.6280 3.1909 -5.5592 50 1 5.9924 3.1702 -3.6914 51 1 7.1073 1.9298 -4.3205 52 1 6.3869 1.4972 -1.9444 53 1 5.9426 0.1828 -3.0603 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 ZINC001323746434.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:11:34 Heat of formation + Delta-G solvation = -96.438232 kcal Electronic energy + Delta-G solvation = -32108.160162 eV Core-core repulsion = 27981.016027 eV Total energy + Delta-G solvation = -4127.144134 eV No. of doubly occupied orbitals = 66 Molecular weight (most abundant/longest-lived isotopes) = 340.220 amu Computer time = 1.07 seconds Orbital eigenvalues (eV) -40.64211 -38.93869 -37.79768 -36.92435 -35.14438 -33.79930 -32.40236 -31.63486 -30.98016 -30.41220 -28.36277 -25.71164 -25.54878 -25.20963 -24.21520 -22.99013 -22.84298 -21.48099 -20.64457 -19.73077 -18.82369 -17.23356 -16.61204 -16.41968 -16.35815 -15.90469 -15.64281 -15.39785 -15.22737 -14.78636 -14.47113 -14.34263 -14.23700 -14.12416 -13.88465 -13.35531 -13.29869 -13.03449 -12.85808 -12.72086 -12.50169 -12.25762 -11.93109 -11.87217 -11.74578 -11.59188 -11.39270 -11.26840 -11.16439 -10.80632 -10.72945 -10.67625 -10.60455 -10.52707 -10.30133 -9.94873 -9.87880 -9.63480 -9.53077 -9.34021 -8.86972 -8.50863 -8.34737 -6.83952 -5.85316 -3.23251 2.65776 3.07329 3.35969 3.49577 3.52511 3.61783 4.38800 4.52959 4.55136 4.62124 4.65013 4.69587 4.81458 4.84998 4.97214 5.08419 5.15766 5.31378 5.38341 5.40973 5.45607 5.46284 5.52401 5.62182 5.63792 5.67571 5.71749 5.74484 5.77904 5.83783 5.89153 5.93285 5.95141 6.02935 6.10208 6.10998 6.12749 6.14068 6.15646 6.27565 6.43054 6.52714 6.69435 6.76885 6.88164 6.97924 7.00187 7.52608 7.96351 8.14878 8.24814 8.47868 9.05007 9.84148 10.17531 11.31022 Molecular weight = 340.22amu Principal moments of inertia in cm(-1) A = 0.015248 B = 0.003186 C = 0.002929 Principal moments of inertia in units of 10**(-40)*gram-cm**2 A = 1835.808210 B = 8786.649734 C = 9558.217033 Net atomic charges, atomic populations, and dipole contributions calculated with CM2 Atom NO. Type Charge No. of electrons 1 C 0.032 3.968 2 O -0.345 6.345 3 C 0.047 3.953 4 C 0.134 3.866 5 N -0.594 5.594 6 C -0.247 4.247 7 H 0.078 0.922 8 C 0.043 3.957 9 N -0.898 5.898 10 Si 0.872 3.128 11 H -0.288 1.288 12 H -0.306 1.306 13 H -0.272 1.272 14 C 0.181 3.819 15 C -0.137 4.137 16 C 0.114 3.886 17 N -0.603 5.603 18 C 0.522 3.478 19 O -0.547 6.547 20 C -0.145 4.145 21 H 0.106 0.894 22 C -0.102 4.102 23 C -0.157 4.157 24 C 0.057 3.943 25 O -0.377 6.377 26 C 0.071 3.929 27 C 0.101 3.899 28 C -0.143 4.143 29 H 0.078 0.922 30 H 0.029 0.971 31 H 0.025 0.975 32 H 0.036 0.964 33 H 0.026 0.974 34 H 0.053 0.947 35 H 0.062 0.938 36 H 0.259 0.741 37 H 0.051 0.949 38 H 0.074 0.926 39 H 0.068 0.932 40 H 0.060 0.940 41 H 0.085 0.915 42 H 0.074 0.926 43 H 0.074 0.926 44 H 0.077 0.923 45 H 0.074 0.926 46 H 0.051 0.949 47 H 0.095 0.905 48 H 0.056 0.944 49 H 0.102 0.898 50 H 0.066 0.934 51 H 0.083 0.917 52 H 0.070 0.930 53 H 0.078 0.922 For the charges calculated by CM2: Dipole moment (debyes) X Y Z Total from point charges 4.897 12.892 -10.491 17.327 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.066 4.066 2 O -0.263 6.263 3 C -0.033 4.033 4 C 0.008 3.992 5 N -0.314 5.314 6 C -0.375 4.375 7 H 0.095 0.905 8 C -0.157 4.157 9 N -0.539 5.539 10 Si 0.703 3.297 11 H -0.217 1.217 12 H -0.235 1.235 13 H -0.198 1.198 14 C 0.074 3.926 15 C -0.177 4.177 16 C -0.009 4.009 17 N -0.336 5.336 18 C 0.310 3.690 19 O -0.424 6.424 20 C -0.166 4.166 21 H 0.124 0.876 22 C -0.140 4.140 23 C -0.194 4.194 24 C -0.019 4.019 25 O -0.295 6.295 26 C -0.007 4.007 27 C -0.022 4.022 28 C -0.183 4.183 29 H 0.097 0.903 30 H 0.047 0.953 31 H 0.043 0.957 32 H 0.054 0.946 33 H 0.044 0.956 34 H 0.072 0.928 35 H 0.079 0.921 36 H 0.069 0.931 37 H 0.069 0.931 38 H 0.092 0.908 39 H 0.087 0.913 40 H 0.079 0.921 41 H 0.103 0.897 42 H 0.092 0.908 43 H 0.092 0.908 44 H 0.095 0.905 45 H 0.093 0.907 46 H 0.069 0.931 47 H 0.113 0.887 48 H 0.074 0.926 49 H 0.120 0.880 50 H 0.085 0.915 51 H 0.101 0.899 52 H 0.088 0.912 53 H 0.097 0.903 Dipole moment (debyes) X Y Z Total from point charges 4.755 12.526 -10.743 17.173 hybrid contribution -0.034 0.383 0.948 1.022 sum 4.721 12.908 -9.795 16.878 Atomic orbital electron populations 1.22637 0.83867 1.01859 0.98198 1.88319 1.16834 1.26352 1.94768 1.23278 0.97393 0.85152 0.97460 1.21348 0.90414 0.96452 0.91004 1.44116 1.75928 1.03754 1.07644 1.18499 1.39698 0.79658 0.99641 0.90496 1.25231 0.98322 0.94743 0.97380 1.69882 1.46027 1.08193 1.29841 0.87047 0.80681 0.79038 0.82901 1.21713 1.23461 1.19763 1.20294 0.91533 0.93120 0.87670 1.22317 1.00566 1.00462 0.94380 1.21689 0.80006 0.99746 0.99439 1.48156 1.07754 1.68169 1.09495 1.20332 0.90293 0.76189 0.82205 1.90661 1.57838 1.51319 1.42631 1.22002 0.93383 1.00164 1.01078 0.87556 1.21409 0.95503 1.01879 0.95238 1.22226 0.96132 0.98128 1.02948 1.22636 1.00870 0.93138 0.85287 1.87969 1.23258 1.82602 1.35660 1.22615 0.81843 0.93935 1.02280 1.22025 0.92744 0.99542 0.87910 1.22389 0.99341 1.00748 0.95802 0.90348 0.95265 0.95669 0.94571 0.95569 0.92840 0.92061 0.93095 0.93112 0.90778 0.91292 0.92135 0.89655 0.90773 0.90777 0.90487 0.90741 0.93098 0.88736 0.92607 0.87979 0.91546 0.89881 0.91185 0.90302 Number of geometries 1 Number of calculations of the screened coulomb radii 1 The total number of SCF iterations 17. 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.53 11.01 101.05 1.11 1.65 16 2 O -0.35 -7.58 10.18 -35.23 -0.36 -7.94 16 3 C 0.05 0.88 6.28 37.16 0.23 1.11 16 4 C 0.13 2.75 4.91 -4.04 -0.02 2.73 16 5 N -0.59 -12.92 2.06 -165.62 -0.34 -13.27 16 6 C -0.25 -6.51 8.29 -15.06 -0.12 -6.63 16 7 H 0.08 2.25 6.92 -52.49 -0.36 1.89 16 8 C 0.04 1.21 5.32 -17.43 -0.09 1.12 16 9 N -0.90 -27.32 13.97 55.49 0.78 -26.55 16 10 Si 0.87 21.80 24.93 -169.99 -4.24 17.56 16 11 H -0.29 -7.17 6.66 56.52 0.38 -6.79 16 12 H -0.31 -7.90 7.11 56.52 0.40 -7.50 16 13 H -0.27 -6.45 7.11 56.52 0.40 -6.05 16 14 C 0.18 3.17 2.59 -67.89 -0.18 3.00 16 15 C -0.14 -2.34 5.26 -26.61 -0.14 -2.48 16 16 C 0.11 1.65 6.42 -3.71 -0.02 1.62 16 17 N -0.60 -8.42 2.97 -172.78 -0.51 -8.94 16 18 C 0.52 7.62 6.95 -10.98 -0.08 7.54 16 19 O -0.55 -9.67 15.37 5.56 0.09 -9.58 16 20 C -0.14 -1.70 1.86 -91.63 -0.17 -1.87 16 21 H 0.11 1.26 7.20 -51.93 -0.37 0.88 16 22 C -0.10 -1.13 4.75 -26.54 -0.13 -1.26 16 23 C -0.16 -1.37 6.14 -26.59 -0.16 -1.53 16 24 C 0.06 0.49 6.84 37.20 0.25 0.74 16 25 O -0.38 -4.09 10.77 -35.23 -0.38 -4.46 16 26 C 0.07 0.73 5.77 37.21 0.21 0.94 16 27 C 0.10 1.25 5.35 -3.72 -0.02 1.23 16 28 C -0.14 -2.32 5.34 -26.61 -0.14 -2.46 16 29 H 0.08 1.27 8.14 -51.93 -0.42 0.85 16 30 H 0.03 0.44 8.14 -51.93 -0.42 0.02 16 31 H 0.02 0.36 8.14 -51.93 -0.42 -0.06 16 32 H 0.04 0.63 8.14 -51.93 -0.42 0.20 16 33 H 0.03 0.43 8.14 -51.93 -0.42 0.01 16 34 H 0.05 0.93 7.94 -51.93 -0.41 0.52 16 35 H 0.06 1.44 2.85 -51.93 -0.15 1.29 16 36 H 0.26 7.57 8.31 -40.82 -0.34 7.23 16 37 H 0.05 0.79 7.81 -51.93 -0.41 0.39 16 38 H 0.07 1.50 6.97 -51.93 -0.36 1.13 16 39 H 0.07 1.20 8.14 -51.93 -0.42 0.77 16 40 H 0.06 0.77 8.14 -51.93 -0.42 0.34 16 41 H 0.09 1.33 7.04 -51.93 -0.37 0.96 16 42 H 0.07 0.89 8.14 -51.93 -0.42 0.46 16 43 H 0.07 0.86 7.86 -51.93 -0.41 0.45 16 44 H 0.08 0.55 8.14 -51.93 -0.42 0.13 16 45 H 0.07 0.67 8.14 -51.93 -0.42 0.25 16 46 H 0.05 0.40 8.14 -51.93 -0.42 -0.02 16 47 H 0.09 0.65 8.14 -51.93 -0.42 0.23 16 48 H 0.06 0.57 8.14 -51.93 -0.42 0.14 16 49 H 0.10 0.96 6.01 -51.93 -0.31 0.64 16 50 H 0.07 0.70 8.14 -51.93 -0.42 0.28 16 51 H 0.08 0.91 4.64 -51.93 -0.24 0.67 16 52 H 0.07 1.16 8.14 -51.93 -0.42 0.73 16 53 H 0.08 1.50 7.75 -51.93 -0.40 1.09 16 LS Contribution 397.62 15.07 5.99 5.99 Total: -1.00 -32.83 397.62 -7.73 -40.56 By element: Atomic # 1 Polarization: 10.45 SS G_CDS: -9.30 Total: 1.16 kcal Atomic # 6 Polarization: 4.92 SS G_CDS: 0.54 Total: 5.46 kcal Atomic # 7 Polarization: -48.67 SS G_CDS: -0.08 Total: -48.75 kcal Atomic # 8 Polarization: -21.33 SS G_CDS: -0.65 Total: -21.99 kcal Atomic # 14 Polarization: 21.80 SS G_CDS: -4.24 Total: 17.56 kcal Total LS contribution 5.99 Total: 5.99 kcal Total: -32.83 -7.73 -40.56 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. ZINC001323746434.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 -55.874 kcal (2) G-P(sol) polarization free energy of solvation -32.831 kcal (3) G-ENP(sol) elect.-nuc.-pol. free energy of system -88.705 kcal (4) G-CDS(sol) cavity-dispersion-solvent structure free energy -7.733 kcal (5) G-P-CDS(sol) = G-P(sol) + G-CDS(sol) = (2) + (4) -40.564 kcal (6) G-S(sol) free energy of system = (1) + (5) -96.438 kcal 1SCF WAS SPECIFIED, SO GEOMETRY OPTIMIZATION WAS NOT USED Total computer time = 1.07 seconds