Wall clock time and date at job start Sat Apr 11 2020 06:14:45 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 N 1.46906 * 1 3 3 C 1.46895 * 110.99949 * 2 1 4 4 C 1.46896 * 110.99584 * 123.95417 * 2 1 3 5 5 C 1.53008 * 109.47008 * 169.99867 * 4 2 1 6 6 C 1.50690 * 109.47320 * 179.97438 * 5 4 2 7 7 O 1.21287 * 120.00098 * 0.02562 * 6 5 4 8 8 N 1.34774 * 120.00056 * 180.02562 * 6 5 4 9 9 C 1.47018 * 125.65029 * 359.95931 * 8 6 5 10 10 C 1.54372 * 107.29103 * 181.13095 * 9 8 6 11 11 C 1.55157 * 102.92010 * 22.13544 * 10 9 8 12 12 H 1.09003 * 111.00188 * 82.64231 * 11 10 9 13 13 C 1.53004 * 111.00251 * 206.45481 * 11 10 9 14 14 N 1.46896 * 109.47166 * 290.41738 * 13 11 10 15 15 C 1.46897 * 111.00124 * 293.95737 * 14 13 11 16 16 C 1.46902 * 110.99925 * 169.99935 * 14 13 11 17 17 C 1.50696 * 109.46978 * 190.00261 * 16 14 13 18 18 O 1.21295 * 120.00238 * 359.97438 * 17 16 14 19 19 N 1.34777 * 120.00076 * 180.02562 * 17 16 14 20 20 C 1.37105 * 120.00244 * 180.02562 * 19 17 16 21 21 H 1.07998 * 119.99708 * 359.97438 * 20 19 17 22 22 C 1.38945 * 120.00017 * 179.97438 * 20 19 17 23 23 N 1.31413 * 119.99901 * 359.71836 * 22 20 19 24 24 Si 1.86800 * 120.00194 * 179.97438 * 22 20 19 25 25 H 1.48501 * 109.47364 * 179.72210 * 24 22 20 26 26 H 1.48508 * 109.46883 * 299.72333 * 24 22 20 27 27 H 1.48497 * 109.46870 * 59.72143 * 24 22 20 28 28 C 1.47431 * 125.64959 * 180.02562 * 8 6 5 29 29 H 1.08996 * 109.47558 * 296.04558 * 1 2 3 30 30 H 1.09002 * 109.46845 * 56.04992 * 1 2 3 31 31 H 1.08997 * 109.47021 * 176.04370 * 1 2 3 32 32 H 1.09001 * 109.47186 * 183.94821 * 3 2 1 33 33 H 1.09001 * 109.47529 * 303.95824 * 3 2 1 34 34 H 1.09007 * 109.47109 * 63.95467 * 3 2 1 35 35 H 1.09004 * 109.47240 * 50.00609 * 4 2 1 36 36 H 1.08996 * 109.47754 * 289.99773 * 4 2 1 37 37 H 1.09004 * 109.46513 * 60.00297 * 5 4 2 38 38 H 1.08996 * 109.46786 * 300.00499 * 5 4 2 39 39 H 1.09001 * 109.88040 * 300.51265 * 9 8 6 40 40 H 1.08994 * 109.88145 * 61.67407 * 9 8 6 41 41 H 1.08999 * 110.69934 * 263.79182 * 10 9 8 42 42 H 1.08997 * 110.75905 * 140.50575 * 10 9 8 43 43 H 1.09000 * 109.46595 * 170.42279 * 13 11 10 44 44 H 1.08995 * 109.47112 * 50.42008 * 13 11 10 45 45 H 1.08997 * 109.47316 * 296.04141 * 15 14 13 46 46 H 1.09003 * 109.47554 * 56.04624 * 15 14 13 47 47 H 1.09007 * 109.46965 * 176.04316 * 15 14 13 48 48 H 1.08999 * 109.46697 * 70.00700 * 16 14 13 49 49 H 1.08993 * 109.47253 * 310.00350 * 16 14 13 50 50 H 0.97004 * 119.99652 * 359.97438 * 19 17 16 51 51 H 0.97000 * 120.00371 * 0.28325 * 23 22 20 52 52 H 1.09003 * 110.36790 * 36.99410 * 28 8 6 53 53 H 1.08997 * 110.36608 * 274.67919 * 28 8 6 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 7 1.4691 0.0000 0.0000 3 6 1.9955 1.3714 0.0000 4 6 1.9954 -0.7660 1.1376 5 6 3.5045 -0.9554 0.9704 6 6 4.0446 -1.7407 2.1376 7 8 3.2962 -2.1127 3.0166 8 7 5.3593 -2.0296 2.2049 9 6 6.3734 -1.6431 1.2131 10 6 7.7251 -2.2139 1.6927 11 6 7.2986 -3.3696 2.6360 12 1 7.0727 -4.2731 2.0696 13 6 8.3731 -3.6414 3.6909 14 7 9.5532 -4.2329 3.0465 15 6 9.2470 -5.5650 2.5083 16 6 10.6883 -4.2894 3.9773 17 6 11.9374 -4.6675 3.2238 18 8 11.8870 -4.8704 2.0290 19 7 13.1105 -4.7811 3.8776 20 6 14.2470 -5.1245 3.1920 21 1 14.2022 -5.3044 2.1281 22 6 15.4565 -5.2411 3.8661 23 7 15.5120 -5.0168 5.1597 24 14 17.0049 -5.7097 2.9321 25 1 18.1536 -5.7603 3.8718 26 1 16.8199 -7.0428 2.3043 27 1 17.2720 -4.6988 1.8777 28 6 6.0185 -2.7868 3.2846 29 1 -0.3634 0.4512 0.9232 30 1 -0.3633 0.5739 -0.8525 31 1 -0.3633 -1.0252 -0.0709 32 1 3.0828 1.3432 -0.0708 33 1 1.5898 1.9164 -0.8524 34 1 1.7044 1.8723 0.9234 35 1 1.5089 -1.7407 1.1755 36 1 1.7975 -0.2249 2.0628 37 1 3.9909 0.0193 0.9324 38 1 3.7023 -1.4966 0.0452 39 1 6.1189 -2.0627 0.2399 40 1 6.4319 -0.5567 1.1469 41 1 8.2902 -1.4587 2.2390 42 1 8.3033 -2.5962 0.8515 43 1 7.9823 -4.3324 4.4378 44 1 8.6535 -2.7052 4.1734 45 1 8.9943 -6.2385 3.3272 46 1 8.4029 -5.4945 1.8222 47 1 10.1166 -5.9505 1.9759 48 1 10.4895 -5.0342 4.7480 49 1 10.8259 -3.3132 4.4422 50 1 13.1508 -4.6191 4.8332 51 1 14.7079 -4.7739 5.6448 52 1 5.3707 -3.5907 3.6343 53 1 6.2794 -2.1230 4.1088 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 ZINC001567736199.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:14:45 Heat of formation + Delta-G solvation = -54.875978 kcal Electronic energy + Delta-G solvation = -29433.828592 eV Core-core repulsion = 25340.779484 eV Total energy + Delta-G solvation = -4093.049109 eV No. of doubly occupied orbitals = 66 Molecular weight (most abundant/longest-lived isotopes) = 340.226 amu Computer time = 0.74 seconds Orbital eigenvalues (eV) -40.45468 -38.94514 -38.23750 -37.49805 -36.20492 -34.22610 -32.59062 -31.40430 -31.08836 -29.68352 -28.03953 -27.91794 -26.95908 -25.42112 -24.38536 -23.29513 -22.31291 -21.48715 -20.19966 -19.55705 -19.32838 -17.89327 -17.55049 -16.92401 -16.35779 -16.05646 -15.75886 -15.53153 -15.31203 -15.17234 -14.84925 -14.51905 -14.30655 -13.98938 -13.72276 -13.71562 -13.58180 -13.40562 -12.98080 -12.82432 -12.80212 -12.56492 -12.29004 -12.24365 -12.13130 -11.93881 -11.81951 -11.48824 -11.38343 -11.30397 -11.22989 -11.14458 -11.02550 -10.87313 -10.61819 -10.21324 -10.11400 -9.86144 -9.51616 -8.92568 -8.84639 -8.52586 -8.05252 -7.92641 -6.75635 -3.88898 2.35490 2.62414 3.30478 3.34372 3.51796 3.54553 3.81521 3.85031 4.23664 4.32813 4.39471 4.44284 4.58877 4.62556 4.72963 4.80160 4.91120 4.93099 4.96505 5.06844 5.07759 5.12828 5.16392 5.21674 5.24901 5.33369 5.43248 5.44227 5.46099 5.48182 5.49696 5.56660 5.57700 5.65326 5.67949 5.76320 5.85111 5.90059 5.95031 6.02679 6.06389 6.15334 6.50041 6.51691 6.52959 7.15767 7.22135 7.30046 7.80523 7.87320 8.07766 8.60044 8.80231 9.47415 9.88378 10.67773 Molecular weight = 340.23amu Principal moments of inertia in cm(-1) A = 0.033630 B = 0.001550 C = 0.001522 Principal moments of inertia in units of 10**(-40)*gram-cm**2 A = 832.398274 B =18060.515528 C =18392.216892 Net atomic charges, atomic populations, and dipole contributions calculated with CM2 Atom NO. Type Charge No. of electrons 1 C 0.025 3.975 2 N -0.550 5.550 3 C 0.018 3.982 4 C 0.070 3.930 5 C -0.133 4.133 6 C 0.529 3.471 7 O -0.546 6.546 8 N -0.618 5.618 9 C 0.096 3.904 10 C -0.114 4.114 11 C -0.096 4.096 12 H 0.088 0.912 13 C 0.065 3.935 14 N -0.515 5.515 15 C 0.019 3.981 16 C 0.046 3.954 17 C 0.447 3.553 18 O -0.571 6.571 19 N -0.577 5.577 20 C -0.382 4.382 21 H 0.099 0.901 22 C -0.009 4.009 23 N -0.901 5.901 24 Si 1.091 2.909 25 H -0.279 1.279 26 H -0.284 1.284 27 H -0.284 1.284 28 C 0.111 3.889 29 H 0.026 0.974 30 H 0.069 0.931 31 H 0.070 0.930 32 H 0.070 0.930 33 H 0.069 0.931 34 H 0.026 0.974 35 H 0.083 0.917 36 H 0.043 0.957 37 H 0.093 0.907 38 H 0.102 0.898 39 H 0.070 0.930 40 H 0.068 0.932 41 H 0.078 0.922 42 H 0.092 0.908 43 H 0.034 0.966 44 H 0.076 0.924 45 H 0.014 0.986 46 H 0.061 0.939 47 H 0.096 0.904 48 H 0.046 0.954 49 H 0.087 0.913 50 H 0.375 0.625 51 H 0.261 0.739 52 H 0.077 0.923 53 H 0.069 0.931 For the charges calculated by CM2: Dipole moment (debyes) X Y Z Total from point charges -25.188 10.314 -6.121 27.898 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.123 4.123 2 N -0.274 5.274 3 C -0.130 4.130 4 C -0.059 4.059 5 C -0.173 4.173 6 C 0.318 3.682 7 O -0.424 6.424 8 N -0.352 5.352 9 C -0.026 4.026 10 C -0.152 4.152 11 C -0.115 4.115 12 H 0.106 0.894 13 C -0.064 4.064 14 N -0.238 5.238 15 C -0.129 4.129 16 C -0.085 4.085 17 C 0.233 3.767 18 O -0.454 6.454 19 N -0.215 5.215 20 C -0.511 4.511 21 H 0.117 0.883 22 C -0.212 4.212 23 N -0.541 5.541 24 Si 0.917 3.083 25 H -0.204 1.204 26 H -0.210 1.210 27 H -0.210 1.210 28 C -0.012 4.012 29 H 0.045 0.955 30 H 0.088 0.912 31 H 0.088 0.912 32 H 0.089 0.911 33 H 0.088 0.912 34 H 0.045 0.955 35 H 0.101 0.899 36 H 0.061 0.939 37 H 0.111 0.889 38 H 0.120 0.880 39 H 0.088 0.912 40 H 0.087 0.913 41 H 0.097 0.903 42 H 0.110 0.890 43 H 0.053 0.947 44 H 0.095 0.905 45 H 0.033 0.967 46 H 0.080 0.920 47 H 0.114 0.886 48 H 0.064 0.936 49 H 0.105 0.895 50 H 0.208 0.792 51 H 0.072 0.928 52 H 0.095 0.905 53 H 0.087 0.913 Dipole moment (debyes) X Y Z Total from point charges -23.765 9.788 -7.140 26.675 hybrid contribution -2.097 0.356 1.871 2.833 sum -25.862 10.144 -5.270 28.275 Atomic orbital electron populations 1.22371 0.86846 1.02514 1.00578 1.61963 1.08256 1.13988 1.43181 1.22475 1.01067 0.88799 1.00657 1.21841 0.93461 0.99196 0.91366 1.21605 0.96789 1.01234 0.97655 1.20298 0.80303 0.79972 0.87613 1.90711 1.56564 1.54127 1.40995 1.48524 1.06825 1.55272 1.24600 1.22205 0.88414 0.99880 0.92150 1.22587 0.94675 0.97397 1.00526 1.22364 0.93204 0.98989 0.96965 0.89397 1.22039 0.90247 0.99683 0.94453 1.61539 1.03792 1.23599 1.34912 1.22596 1.03174 0.88205 0.98876 1.21865 0.91420 1.01513 0.93750 1.19651 0.84104 0.84633 0.88347 1.90483 1.85369 1.54759 1.14836 1.41948 1.07630 1.62050 1.09865 1.20522 0.86838 1.47515 0.96220 0.88332 1.26493 1.01647 0.96370 0.96661 1.70903 1.31311 1.44339 1.07514 0.83136 0.73987 0.76696 0.74484 1.20439 1.21049 1.21004 1.22224 0.93363 0.95581 0.90033 0.95513 0.91209 0.91166 0.91099 0.91226 0.95505 0.89884 0.93882 0.88914 0.87953 0.91179 0.91338 0.90293 0.88977 0.94746 0.90548 0.96741 0.92027 0.88580 0.93572 0.89502 0.79212 0.92788 0.90455 0.91292 Number of geometries 1 Number of calculations of the screened coulomb radii 1 The total number of SCF iterations 13. 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.15 9.57 60.07 0.57 0.73 16 2 N -0.55 -4.06 5.02 -251.25 -1.26 -5.32 16 3 C 0.02 0.11 9.32 60.06 0.56 0.67 16 4 C 0.07 0.58 4.25 -3.82 -0.02 0.57 16 5 C -0.13 -1.02 4.82 -27.88 -0.13 -1.16 16 6 C 0.53 5.72 7.85 -10.99 -0.09 5.64 16 7 O -0.55 -7.87 16.24 5.55 0.09 -7.78 16 8 N -0.62 -6.05 3.33 -170.51 -0.57 -6.62 16 9 C 0.10 0.75 6.49 -3.04 -0.02 0.73 16 10 C -0.11 -1.06 6.32 -24.87 -0.16 -1.22 16 11 C -0.10 -0.96 2.94 -88.48 -0.26 -1.22 16 12 H 0.09 0.93 6.34 -51.93 -0.33 0.60 16 13 C 0.06 0.74 4.40 -3.82 -0.02 0.72 16 14 N -0.51 -7.88 4.39 -244.38 -1.07 -8.95 16 15 C 0.02 0.30 9.03 60.06 0.54 0.84 16 16 C 0.05 0.78 5.45 -4.98 -0.03 0.75 16 17 C 0.45 10.25 7.51 -10.99 -0.08 10.17 16 18 O -0.57 -14.99 13.88 5.54 0.08 -14.91 16 19 N -0.58 -14.25 5.28 -10.96 -0.06 -14.31 16 20 C -0.38 -10.58 9.68 -14.93 -0.14 -10.73 16 21 H 0.10 2.89 7.16 -52.49 -0.38 2.51 16 22 C -0.01 -0.25 5.50 -17.47 -0.10 -0.35 16 23 N -0.90 -25.04 13.68 55.50 0.76 -24.28 16 24 Si 1.09 24.02 29.92 -169.99 -5.09 18.94 16 25 H -0.28 -6.10 7.11 56.52 0.40 -5.70 16 26 H -0.28 -6.03 7.11 56.52 0.40 -5.63 16 27 H -0.28 -6.03 7.11 56.52 0.40 -5.63 16 28 C 0.11 1.11 6.38 -2.54 -0.02 1.09 16 29 H 0.03 0.16 8.14 -51.93 -0.42 -0.26 16 30 H 0.07 0.35 8.12 -51.93 -0.42 -0.07 16 31 H 0.07 0.41 8.12 -51.93 -0.42 -0.01 16 32 H 0.07 0.37 6.47 -51.93 -0.34 0.04 16 33 H 0.07 0.35 8.12 -51.93 -0.42 -0.07 16 34 H 0.03 0.16 8.14 -51.92 -0.42 -0.26 16 35 H 0.08 0.77 8.02 -51.93 -0.42 0.35 16 36 H 0.04 0.40 8.08 -51.93 -0.42 -0.01 16 37 H 0.09 0.58 6.55 -51.93 -0.34 0.24 16 38 H 0.10 0.71 8.14 -51.93 -0.42 0.29 16 39 H 0.07 0.48 8.14 -51.93 -0.42 0.06 16 40 H 0.07 0.44 8.14 -51.93 -0.42 0.02 16 41 H 0.08 0.75 7.97 -51.93 -0.41 0.34 16 42 H 0.09 0.96 8.14 -51.93 -0.42 0.53 16 43 H 0.03 0.38 8.12 -51.93 -0.42 -0.04 16 44 H 0.08 0.82 7.83 -51.93 -0.41 0.41 16 45 H 0.01 0.20 8.14 -51.93 -0.42 -0.22 16 46 H 0.06 0.81 6.23 -51.93 -0.32 0.49 16 47 H 0.10 1.84 5.88 -51.93 -0.31 1.53 16 48 H 0.05 0.73 8.12 -51.93 -0.42 0.30 16 49 H 0.09 1.32 7.93 -51.93 -0.41 0.91 16 50 H 0.37 8.58 6.73 -40.82 -0.27 8.30 16 51 H 0.26 6.96 6.78 -40.82 -0.28 6.68 16 52 H 0.08 0.84 8.04 -51.93 -0.42 0.42 16 53 H 0.07 0.67 8.03 -51.93 -0.42 0.25 16 LS Contribution 418.19 15.07 6.30 6.30 Total: -1.00 -33.82 418.19 -9.52 -43.34 By element: Atomic # 1 Polarization: 15.69 SS G_CDS: -9.33 Total: 6.36 kcal Atomic # 6 Polarization: 6.61 SS G_CDS: 0.62 Total: 7.23 kcal Atomic # 7 Polarization: -57.28 SS G_CDS: -2.20 Total: -59.48 kcal Atomic # 8 Polarization: -22.86 SS G_CDS: 0.17 Total: -22.70 kcal Atomic # 14 Polarization: 24.02 SS G_CDS: -5.09 Total: 18.94 kcal Total LS contribution 6.30 Total: 6.30 kcal Total: -33.82 -9.52 -43.34 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. ZINC001567736199.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 -11.533 kcal (2) G-P(sol) polarization free energy of solvation -33.819 kcal (3) G-ENP(sol) elect.-nuc.-pol. free energy of system -45.352 kcal (4) G-CDS(sol) cavity-dispersion-solvent structure free energy -9.524 kcal (5) G-P-CDS(sol) = G-P(sol) + G-CDS(sol) = (2) + (4) -43.343 kcal (6) G-S(sol) free energy of system = (1) + (5) -54.876 kcal 1SCF WAS SPECIFIED, SO GEOMETRY OPTIMIZATION WAS NOT USED Total computer time = 0.74 seconds