Wall clock time and date at job start Sat Apr 11 2020 06:13:04 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.42899 * 1 3 3 C 1.42908 * 113.99702 * 2 1 4 4 H 1.08997 * 109.47296 * 29.99762 * 3 2 1 5 5 C 1.50700 * 109.46961 * 269.99842 * 3 2 1 6 6 O 1.21287 * 119.99667 * 354.99950 * 5 3 2 7 7 N 1.34766 * 120.00030 * 174.99949 * 5 3 2 8 8 C 1.46931 * 120.63392 * 0.02562 * 7 5 3 9 9 C 1.53191 * 108.77528 * 233.61025 * 8 7 5 10 10 C 1.53045 * 109.31249 * 305.36071 * 9 8 7 11 11 C 1.53039 * 109.53585 * 61.36393 * 10 9 8 12 12 H 1.09004 * 109.49569 * 58.62353 * 11 10 9 13 13 N 1.46499 * 109.49376 * 178.68266 * 11 10 9 14 14 C 1.36939 * 120.00328 * 59.93922 * 13 11 10 15 15 H 1.08003 * 119.99470 * 325.01956 * 14 13 11 16 16 C 1.38807 * 120.00382 * 145.01147 * 14 13 11 17 17 N 1.31621 * 120.00179 * 156.84059 * 16 14 13 18 18 Si 1.86803 * 119.99936 * 336.83313 * 16 14 13 19 19 H 1.48508 * 109.46692 * 95.14397 * 18 16 14 20 20 H 1.48504 * 109.47089 * 215.14211 * 18 16 14 21 21 H 1.48496 * 109.47448 * 335.14623 * 18 16 14 22 22 C 1.46504 * 120.00152 * 239.93279 * 13 11 10 23 23 C 1.53006 * 109.46606 * 84.34636 * 22 13 11 24 24 O 1.42898 * 109.46910 * 63.57741 * 23 22 13 25 25 C 1.46936 * 120.60825 * 180.02562 * 7 5 3 26 26 C 1.53005 * 109.46883 * 150.00150 * 3 2 1 27 27 C 1.53778 * 113.61736 * 39.13383 * 26 3 2 28 28 C 1.53958 * 87.01918 * 220.04833 * 27 26 3 29 29 C 1.53782 * 87.11145 * 25.36401 * 28 27 26 30 30 H 1.08997 * 109.46670 * 180.02562 * 1 2 3 31 31 H 1.09001 * 109.47156 * 299.99914 * 1 2 3 32 32 H 1.08999 * 109.47140 * 60.00529 * 1 2 3 33 33 H 1.08993 * 109.58767 * 113.81927 * 8 7 5 34 34 H 1.08998 * 109.58128 * 353.39350 * 8 7 5 35 35 H 1.08999 * 109.50150 * 185.40878 * 9 8 7 36 36 H 1.09006 * 109.52606 * 65.32732 * 9 8 7 37 37 H 1.08999 * 109.46083 * 301.34481 * 10 9 8 38 38 H 1.08999 * 109.44709 * 181.38060 * 10 9 8 39 39 H 0.97002 * 120.00049 * 180.02562 * 17 16 14 40 40 H 1.08996 * 109.47143 * 324.34696 * 22 13 11 41 41 H 1.09001 * 109.46945 * 204.34261 * 22 13 11 42 42 H 1.09000 * 109.46531 * 183.57439 * 23 22 13 43 43 H 1.09000 * 109.47105 * 303.57397 * 23 22 13 44 44 H 0.96700 * 113.99783 * 180.02562 * 24 23 22 45 45 H 1.09011 * 109.61051 * 6.62374 * 25 7 5 46 46 H 1.08997 * 109.56195 * 246.19995 * 25 7 5 47 47 H 1.09002 * 112.84262 * 170.33342 * 26 3 2 48 48 H 1.08998 * 113.61180 * 334.57113 * 27 26 3 49 49 H 1.09004 * 113.61530 * 105.47061 * 27 26 3 50 50 H 1.08999 * 113.53585 * 270.81635 * 28 27 26 51 51 H 1.09000 * 113.66572 * 139.93513 * 28 27 26 52 52 H 1.09006 * 113.61422 * 220.09183 * 29 28 27 53 53 H 1.08997 * 113.61672 * 89.18776 * 29 28 27 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.3056 0.0000 4 1 1.3449 1.9995 0.5138 5 6 2.2145 1.7644 -1.4208 6 8 1.9791 1.0105 -2.3413 7 7 2.6587 3.0124 -1.6686 8 6 2.9487 3.9370 -0.5641 9 6 2.1735 5.2381 -0.7938 10 6 2.5111 5.7903 -2.1807 11 6 2.0881 4.7808 -3.2503 12 1 1.0188 4.5873 -3.1654 13 7 2.3829 5.3236 -4.5787 14 6 3.6704 5.6450 -4.9169 15 1 4.3409 6.0448 -4.1705 16 6 4.1137 5.4568 -6.2187 17 7 5.1371 6.1493 -6.6720 18 14 3.2580 4.2206 -7.3273 19 1 2.2409 4.9184 -8.1543 20 1 4.2590 3.5750 -8.2142 21 1 2.5938 3.1858 -6.4948 22 6 1.3003 5.5225 -5.5455 23 6 0.6286 6.8729 -5.2882 24 8 1.5709 7.9238 -5.5111 25 6 2.8621 3.4739 -3.0487 26 6 3.3589 1.2634 0.7212 27 6 4.1773 -0.0041 0.4237 28 6 4.5679 -0.0044 1.9129 29 6 3.2650 0.7695 2.1746 30 1 -0.3632 -1.0277 0.0005 31 1 -0.3633 0.5138 0.8900 32 1 -0.3633 0.5137 -0.8900 33 1 4.0175 4.1485 -0.5357 34 1 2.6378 3.4881 0.3792 35 1 2.4535 5.9676 -0.0338 36 1 1.1033 5.0405 -0.7315 37 1 3.5847 5.9650 -2.2506 38 1 1.9795 6.7290 -2.3368 39 1 5.4467 6.0181 -7.5819 40 1 0.5655 4.7247 -5.4375 41 1 1.7083 5.5064 -6.5562 42 1 -0.2168 6.9931 -5.9657 43 1 0.2766 6.9131 -4.2574 44 1 1.2131 8.8103 -5.3658 45 1 2.4939 2.7200 -3.7446 46 1 3.9240 3.6472 -3.2228 47 1 3.9275 2.1866 0.6096 48 1 3.5674 -0.8591 0.1320 49 1 5.0178 0.1667 -0.2491 50 1 5.4683 0.5714 2.1268 51 1 4.5831 -0.9980 2.3610 52 1 3.3701 1.5641 2.9134 53 1 2.4053 0.1264 2.3626 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 ZINC001479657288.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:13:04 Heat of formation + Delta-G solvation = -70.594947 kcal Electronic energy + Delta-G solvation = -32540.791883 eV Core-core repulsion = 28414.768398 eV Total energy + Delta-G solvation = -4126.023485 eV No. of doubly occupied orbitals = 66 Molecular weight (most abundant/longest-lived isotopes) = 340.220 amu Computer time = 0.93 seconds Orbital eigenvalues (eV) -40.64810 -38.81444 -37.64849 -36.86024 -36.20792 -34.75112 -32.32857 -31.45461 -30.61759 -29.33907 -27.51418 -26.42583 -25.75171 -25.21297 -24.09595 -23.17263 -22.58956 -21.07543 -19.88286 -19.55702 -18.86861 -17.23076 -16.83599 -16.71810 -16.41080 -16.27959 -15.46899 -15.39051 -15.19520 -14.77777 -14.35699 -14.29544 -13.98903 -13.67734 -13.58200 -13.42895 -13.15291 -12.85462 -12.75017 -12.48772 -12.46518 -12.27798 -12.20062 -11.98063 -11.80227 -11.61569 -11.43320 -11.11188 -10.82215 -10.76656 -10.57683 -10.56520 -10.49916 -10.41542 -10.32753 -10.27803 -10.03077 -9.66499 -9.56319 -9.17979 -9.02818 -8.40619 -8.30444 -6.83205 -5.87803 -3.26237 2.77947 3.09512 3.37586 3.54675 3.59120 3.66152 4.28491 4.50562 4.55406 4.61086 4.68246 4.73076 4.77747 4.98080 5.03115 5.06063 5.11811 5.17062 5.23253 5.32129 5.36521 5.45729 5.49949 5.56489 5.60033 5.63931 5.70704 5.74602 5.75631 5.81861 5.87006 5.92749 5.96736 6.00882 6.08401 6.13032 6.13781 6.23047 6.34130 6.40845 6.47325 6.57042 6.70243 6.74592 6.88216 6.97852 6.98316 7.68473 8.09348 8.40016 8.42923 8.48466 9.02849 9.81773 10.16320 11.28040 Molecular weight = 340.22amu Principal moments of inertia in cm(-1) A = 0.015057 B = 0.003414 C = 0.003321 Principal moments of inertia in units of 10**(-40)*gram-cm**2 A = 1859.128422 B = 8199.348322 C = 8428.787046 Net atomic charges, atomic populations, and dipole contributions calculated with CM2 Atom NO. Type Charge No. of electrons 1 C 0.030 3.970 2 O -0.354 6.354 3 C 0.085 3.915 4 H 0.083 0.917 5 C 0.489 3.511 6 O -0.531 6.531 7 N -0.604 5.604 8 C 0.103 3.897 9 C -0.133 4.133 10 C -0.125 4.125 11 C 0.170 3.830 12 H 0.050 0.950 13 N -0.589 5.589 14 C -0.251 4.251 15 H 0.079 0.921 16 C 0.043 3.957 17 N -0.897 5.897 18 Si 0.867 3.133 19 H -0.288 1.288 20 H -0.303 1.303 21 H -0.269 1.269 22 C 0.131 3.869 23 C 0.063 3.937 24 O -0.535 6.535 25 C 0.101 3.899 26 C -0.117 4.117 27 C -0.127 4.127 28 C -0.141 4.141 29 C -0.128 4.128 30 H 0.088 0.912 31 H 0.040 0.960 32 H 0.045 0.955 33 H 0.072 0.928 34 H 0.079 0.921 35 H 0.070 0.930 36 H 0.059 0.941 37 H 0.071 0.929 38 H 0.060 0.940 39 H 0.260 0.740 40 H 0.049 0.951 41 H 0.068 0.932 42 H 0.023 0.977 43 H 0.030 0.970 44 H 0.353 0.647 45 H 0.083 0.917 46 H 0.081 0.919 47 H 0.110 0.890 48 H 0.094 0.906 49 H 0.065 0.935 50 H 0.078 0.922 51 H 0.067 0.933 52 H 0.066 0.934 53 H 0.086 0.914 For the charges calculated by CM2: Dipole moment (debyes) X Y Z Total from point charges -8.627 -5.793 19.045 21.696 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.067 4.067 2 O -0.271 6.271 3 C 0.023 3.977 4 H 0.101 0.899 5 C 0.276 3.724 6 O -0.407 6.407 7 N -0.339 5.339 8 C -0.021 4.021 9 C -0.171 4.171 10 C -0.164 4.164 11 C 0.062 3.938 12 H 0.069 0.931 13 N -0.310 5.310 14 C -0.379 4.379 15 H 0.097 0.903 16 C -0.158 4.158 17 N -0.538 5.538 18 Si 0.698 3.302 19 H -0.217 1.217 20 H -0.231 1.231 21 H -0.194 1.194 22 C 0.005 3.995 23 C -0.018 4.018 24 O -0.336 6.336 25 C -0.022 4.022 26 C -0.135 4.135 27 C -0.165 4.165 28 C -0.179 4.179 29 C -0.166 4.166 30 H 0.106 0.894 31 H 0.058 0.942 32 H 0.063 0.937 33 H 0.090 0.910 34 H 0.097 0.903 35 H 0.089 0.911 36 H 0.077 0.923 37 H 0.090 0.910 38 H 0.079 0.921 39 H 0.070 0.930 40 H 0.067 0.933 41 H 0.086 0.914 42 H 0.041 0.959 43 H 0.049 0.951 44 H 0.196 0.804 45 H 0.101 0.899 46 H 0.099 0.901 47 H 0.128 0.872 48 H 0.113 0.887 49 H 0.084 0.916 50 H 0.097 0.903 51 H 0.085 0.915 52 H 0.085 0.915 53 H 0.104 0.896 Dipole moment (debyes) X Y Z Total from point charges -8.115 -6.888 18.252 21.129 hybrid contribution -1.119 0.246 -0.749 1.369 sum -9.235 -6.642 17.503 20.875 Atomic orbital electron populations 1.22833 0.80929 1.02981 0.99924 1.87939 1.19668 1.25076 1.94419 1.21849 0.96683 0.86316 0.92819 0.89927 1.21200 0.77661 0.81743 0.91802 1.90697 1.50612 1.53302 1.46074 1.48304 1.65182 1.12154 1.08303 1.21900 1.00426 0.90510 0.89217 1.21952 0.99893 0.96259 0.99039 1.21941 1.02888 0.99104 0.92467 1.20611 0.96056 0.90691 0.86422 0.93145 1.44013 1.04221 1.72140 1.10602 1.18520 0.86936 1.36426 0.96035 0.90342 1.25278 0.95666 0.97666 0.97141 1.69878 1.28356 1.39290 1.16298 0.87121 0.78419 0.83736 0.80949 1.21673 1.23091 1.19390 1.21363 0.90663 0.92083 0.95397 1.22869 0.90735 0.91050 0.97147 1.86932 1.38832 1.15506 1.92356 1.22326 1.02436 0.97312 0.80087 1.22890 0.93838 1.00507 0.96268 1.22850 0.99347 0.98015 0.96261 1.22817 0.97432 1.00588 0.97027 1.22886 0.98791 0.99930 0.94970 0.89398 0.94188 0.93673 0.90980 0.90325 0.91129 0.92255 0.90988 0.92102 0.92991 0.93274 0.91431 0.95868 0.95145 0.80377 0.89863 0.90055 0.87191 0.88735 0.91613 0.90308 0.91456 0.91481 0.89557 Number of geometries 1 Number of calculations of the screened coulomb radii 1 The total number of SCF iterations 16. 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.31 10.81 101.05 1.09 1.40 16 2 O -0.35 -4.59 7.99 -55.09 -0.44 -5.03 16 3 C 0.08 1.03 1.70 -27.88 -0.05 0.98 16 4 H 0.08 0.82 7.08 -51.93 -0.37 0.45 16 5 C 0.49 7.97 7.05 -10.99 -0.08 7.89 16 6 O -0.53 -11.04 15.92 -14.31 -0.23 -11.26 16 7 N -0.60 -9.43 2.97 -172.74 -0.51 -9.94 16 8 C 0.10 1.20 4.90 -3.71 -0.02 1.18 16 9 C -0.13 -1.59 6.08 -26.61 -0.16 -1.75 16 10 C -0.12 -2.04 5.11 -26.69 -0.14 -2.18 16 11 C 0.17 3.24 2.78 -67.80 -0.19 3.05 16 12 H 0.05 0.92 7.87 -51.93 -0.41 0.51 16 13 N -0.59 -13.46 2.01 -165.51 -0.33 -13.79 16 14 C -0.25 -6.77 8.35 -15.06 -0.13 -6.89 16 15 H 0.08 2.28 7.00 -52.48 -0.37 1.91 16 16 C 0.04 1.23 5.27 -17.43 -0.09 1.14 16 17 N -0.90 -27.63 13.97 55.49 0.78 -26.86 16 18 Si 0.87 22.74 24.99 -169.99 -4.25 18.49 16 19 H -0.29 -7.42 6.64 56.52 0.38 -7.05 16 20 H -0.30 -8.12 7.11 56.52 0.40 -7.72 16 21 H -0.27 -6.98 7.11 56.52 0.40 -6.58 16 22 C 0.13 2.85 4.66 -4.04 -0.02 2.83 16 23 C 0.06 1.16 7.39 37.16 0.27 1.44 16 24 O -0.54 -11.09 13.34 -35.23 -0.47 -11.56 16 25 C 0.10 1.97 5.76 -3.69 -0.02 1.95 16 26 C -0.12 -1.20 2.88 -89.80 -0.26 -1.46 16 27 C -0.13 -1.44 7.48 -25.83 -0.19 -1.63 16 28 C -0.14 -1.22 8.09 -25.82 -0.21 -1.43 16 29 C -0.13 -1.06 7.54 -25.92 -0.20 -1.26 16 30 H 0.09 0.77 8.14 -51.93 -0.42 0.34 16 31 H 0.04 0.32 8.09 -51.93 -0.42 -0.10 16 32 H 0.04 0.53 8.14 -51.93 -0.42 0.11 16 33 H 0.07 0.83 8.13 -51.93 -0.42 0.41 16 34 H 0.08 0.70 4.83 -51.93 -0.25 0.45 16 35 H 0.07 0.70 8.14 -51.93 -0.42 0.28 16 36 H 0.06 0.68 8.14 -51.93 -0.42 0.26 16 37 H 0.07 1.33 7.03 -51.93 -0.37 0.96 16 38 H 0.06 1.03 8.14 -51.93 -0.42 0.61 16 39 H 0.26 7.69 8.31 -40.82 -0.34 7.35 16 40 H 0.05 0.96 8.04 -51.93 -0.42 0.55 16 41 H 0.07 1.67 3.27 -51.93 -0.17 1.50 16 42 H 0.02 0.38 8.14 -51.93 -0.42 -0.05 16 43 H 0.03 0.51 8.14 -51.93 -0.42 0.09 16 44 H 0.35 5.28 9.12 45.56 0.42 5.70 16 45 H 0.08 1.86 7.01 -51.92 -0.36 1.50 16 46 H 0.08 1.74 7.50 -51.93 -0.39 1.35 16 47 H 0.11 1.06 5.83 -51.93 -0.30 0.76 16 48 H 0.09 1.21 7.03 -51.93 -0.36 0.84 16 49 H 0.07 0.79 8.14 -51.93 -0.42 0.37 16 50 H 0.08 0.62 8.14 -51.93 -0.42 0.20 16 51 H 0.07 0.54 8.14 -51.93 -0.42 0.12 16 52 H 0.07 0.46 8.14 -51.93 -0.42 0.03 16 53 H 0.09 0.74 7.72 -51.93 -0.40 0.34 16 LS Contribution 401.32 15.07 6.05 6.05 Total: -1.00 -34.97 401.32 -8.19 -43.16 By element: Atomic # 1 Polarization: 13.89 SS G_CDS: -8.41 Total: 5.48 kcal Atomic # 6 Polarization: 5.65 SS G_CDS: -0.38 Total: 5.27 kcal Atomic # 7 Polarization: -50.52 SS G_CDS: -0.07 Total: -50.59 kcal Atomic # 8 Polarization: -26.72 SS G_CDS: -1.14 Total: -27.86 kcal Atomic # 14 Polarization: 22.74 SS G_CDS: -4.25 Total: 18.49 kcal Total LS contribution 6.05 Total: 6.05 kcal Total: -34.97 -8.19 -43.16 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. ZINC001479657288.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 -27.436 kcal (2) G-P(sol) polarization free energy of solvation -34.965 kcal (3) G-ENP(sol) elect.-nuc.-pol. free energy of system -62.402 kcal (4) G-CDS(sol) cavity-dispersion-solvent structure free energy -8.193 kcal (5) G-P-CDS(sol) = G-P(sol) + G-CDS(sol) = (2) + (4) -43.159 kcal (6) G-S(sol) free energy of system = (1) + (5) -70.595 kcal 1SCF WAS SPECIFIED, SO GEOMETRY OPTIMIZATION WAS NOT USED Total computer time = 0.93 seconds