Wall clock time and date at job start Fri Apr 10 2020 22:20:37 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 * SM5.42R - SM5.42R CALCULATIONS WILL BE PERFORMED * - CHARGE MODEL 2 WILL BE USED (CM2A) * - THE SOLVENT IS WATER ******************************************************************************* 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 C 1.53001 * 1 3 3 C 1.53007 * 109.47171 * 2 1 4 4 C 1.52994 * 109.47394 * 180.02562 * 3 2 1 5 5 C 1.52997 * 109.47247 * 180.02562 * 4 3 2 6 6 C 1.53006 * 109.47082 * 184.91480 * 5 4 3 7 7 C 1.53000 * 109.46986 * 175.37964 * 6 5 4 8 8 C 1.52996 * 109.46946 * 185.74924 * 7 6 5 9 9 C 1.53004 * 109.47140 * 184.69995 * 8 7 6 10 10 C 1.52996 * 109.47140 * 180.02562 * 9 8 7 11 11 C 1.53004 * 109.47140 * 180.02562 * 10 9 8 12 12 N 1.46895 * 109.46811 * 295.38168 * 6 5 4 13 13 C 1.46904 * 111.00121 * 100.25144 * 12 6 5 14 14 H 1.08999 * 109.41280 * 323.51110 * 13 12 6 15 15 C 1.52958 * 109.35531 * 83.33101 * 13 12 6 16 16 C 1.53228 * 109.40281 * 183.67087 * 15 13 12 17 17 C 1.51266 * 108.53720 * 305.37259 * 16 15 13 18 18 C 1.37884 * 120.85053 * 231.07337 * 17 16 15 19 19 N 1.31515 * 119.99704 * 180.28075 * 18 17 16 20 20 Si 1.86797 * 119.99936 * 0.27839 * 18 17 16 21 21 H 1.48506 * 109.46940 * 60.00304 * 20 18 17 22 22 H 1.48495 * 109.46806 * 179.97438 * 20 18 17 23 23 H 1.48495 * 109.47324 * 299.99993 * 20 18 17 24 24 C 1.51274 * 118.30086 * 51.34797 * 17 16 15 25 25 C 1.52955 * 109.41184 * 203.65123 * 13 12 6 26 26 H 1.08995 * 109.46949 * 180.02562 * 1 2 3 27 27 H 1.09007 * 109.46762 * 300.00093 * 1 2 3 28 28 H 1.08998 * 109.47104 * 60.00159 * 1 2 3 29 29 H 1.08997 * 109.47163 * 240.00177 * 2 1 3 30 30 H 1.09001 * 109.47243 * 119.99088 * 2 1 3 31 31 H 1.08998 * 109.47097 * 59.99511 * 3 2 1 32 32 H 1.09007 * 109.47003 * 300.00191 * 3 2 1 33 33 H 1.08999 * 109.47149 * 59.99243 * 4 3 2 34 34 H 1.08998 * 109.46846 * 299.99819 * 4 3 2 35 35 H 1.08997 * 109.47388 * 64.91641 * 5 4 3 36 36 H 1.09001 * 109.46923 * 304.91368 * 5 4 3 37 37 H 1.09001 * 109.46879 * 55.38425 * 6 5 4 38 38 H 1.09000 * 109.47265 * 305.75001 * 7 6 5 39 39 H 1.08994 * 109.47038 * 65.74993 * 7 6 5 40 40 H 1.08998 * 109.47205 * 304.70179 * 8 7 6 41 41 H 1.09001 * 109.47079 * 64.69975 * 8 7 6 42 42 H 1.08998 * 109.47222 * 300.00174 * 9 8 7 43 43 H 1.09001 * 109.47150 * 60.00024 * 9 8 7 44 44 H 1.09007 * 109.47351 * 300.00286 * 10 9 8 45 45 H 1.09001 * 109.47079 * 59.99980 * 10 9 8 46 46 H 1.09001 * 109.46601 * 180.02562 * 11 10 9 47 47 H 1.08998 * 109.47222 * 300.00174 * 11 10 9 48 48 H 1.09000 * 109.46475 * 60.00345 * 11 10 9 49 49 H 1.00900 * 111.00653 * 336.29448 * 12 6 5 50 50 H 1.09001 * 109.48552 * 303.64393 * 15 13 12 51 51 H 1.08994 * 109.48420 * 63.68038 * 15 13 12 52 52 H 1.09003 * 109.63136 * 185.65807 * 16 15 13 53 53 H 1.09008 * 109.62753 * 65.08517 * 16 15 13 54 54 H 0.97004 * 119.99570 * 180.02562 * 19 18 17 55 55 H 1.09004 * 109.62547 * 68.36645 * 24 17 16 56 56 H 1.08999 * 109.62987 * 188.93890 * 24 17 16 57 57 H 1.08998 * 109.48510 * 56.38602 * 25 13 12 58 58 H 1.09003 * 109.49251 * 296.34259 * 25 13 12 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 6 1.5300 0.0000 0.0000 3 6 2.0400 1.4426 0.0000 4 6 3.5700 1.4426 -0.0006 5 6 4.0800 2.8851 0.0000 6 6 5.6050 2.8870 -0.1242 7 6 6.1055 4.3282 -0.2395 8 6 7.6118 4.3250 -0.5071 9 6 8.1331 5.7635 -0.5059 10 6 9.6393 5.7604 -0.7742 11 6 10.1607 7.1988 -0.7724 12 7 6.1938 2.2577 1.0653 13 6 6.5799 0.8675 0.7889 14 1 5.8490 0.4149 0.1188 15 6 7.9585 0.8464 0.1265 16 6 8.3225 -0.5948 -0.2455 17 6 8.2138 -1.4537 0.9949 18 6 9.2502 -2.2792 1.3764 19 7 9.1356 -3.0229 2.4550 20 14 10.8171 -2.3410 0.3614 21 1 11.4320 -0.9899 0.3174 22 1 11.7697 -3.2980 0.9793 23 1 10.4955 -2.7837 -1.0191 24 6 6.9364 -1.3886 1.8026 25 6 6.6226 0.0814 2.1003 26 1 -0.3633 -1.0276 0.0005 27 1 -0.3633 0.5139 0.8901 28 1 -0.3633 0.5138 -0.8900 29 1 1.8933 -0.5138 -0.8900 30 1 1.8934 -0.5137 0.8901 31 1 1.6766 1.9564 0.8899 32 1 1.6767 1.9564 -0.8900 33 1 3.9330 0.9288 -0.8908 34 1 3.9337 0.9287 0.8891 35 1 3.7907 3.3723 0.9311 36 1 3.6461 3.4242 -0.8422 37 1 5.8977 2.3285 -1.0134 38 1 5.5915 4.8267 -1.0613 39 1 5.9029 4.8586 0.6909 40 1 8.1180 3.7545 0.2716 41 1 7.8075 3.8684 -1.4774 42 1 7.6269 6.3340 -1.2846 43 1 7.9374 6.2201 0.4643 44 1 10.1459 5.1897 0.0042 45 1 9.8346 5.3039 -1.7446 46 1 11.2338 7.1965 -0.9640 47 1 9.6543 7.7698 -1.5505 48 1 9.9655 7.6547 0.1983 49 1 5.5629 2.3083 1.8511 50 1 7.9385 1.4602 -0.7741 51 1 8.7010 1.2417 0.8196 52 1 9.3435 -0.6282 -0.6255 53 1 7.6335 -0.9617 -1.0064 54 1 9.8649 -3.6034 2.7236 55 1 6.1191 -1.8290 1.2313 56 1 7.0668 -1.9319 2.7385 57 1 5.6560 0.1531 2.5988 58 1 7.3968 0.4945 2.7470 RHF calculation, no. of doubly occupied orbitals= 62 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=WATER ZINC001257550451.mol2 58 1SCF WAS SPECIFIED, SO GEOMETRY OPTIMIZATION WAS NOT USED SCF FIELD WAS ACHIEVED Wall clock time and date at job start Fri Apr 10 2020 22:20:37 Heat of formation + Delta-G solvation = -89.435493 kcal Electronic energy + Delta-G solvation = -27311.190761 eV Core-core repulsion = 23995.662199 eV Total energy + Delta-G solvation = -3315.528562 eV No. of doubly occupied orbitals = 62 Molecular weight (most abundant/longest-lived isotopes) = 309.273 amu Computer time = 1.16 seconds Orbital eigenvalues (eV) -40.00051 -37.79378 -36.89164 -35.83061 -34.22683 -32.64654 -32.10439 -31.34397 -29.65242 -27.62470 -26.90275 -24.21664 -23.45771 -23.17019 -22.28395 -21.94992 -21.58684 -20.96499 -20.67043 -17.97174 -17.09803 -16.58489 -16.41955 -16.27019 -15.79380 -15.36452 -15.04259 -14.74930 -14.66470 -14.50337 -14.43411 -14.25834 -14.12993 -13.93973 -13.67662 -13.49309 -13.32605 -12.91991 -12.75970 -12.69985 -12.67645 -12.55898 -12.39420 -12.19868 -12.01702 -11.87733 -11.79812 -11.75269 -11.71013 -11.68162 -11.58081 -11.47971 -11.37582 -11.36499 -11.21375 -10.95675 -10.71038 -10.53838 -10.43752 -9.04392 -8.16743 -6.09298 1.37551 1.38504 1.46476 2.40244 3.08312 3.42847 3.65321 3.99745 4.03747 4.11971 4.13141 4.19992 4.23120 4.29830 4.32543 4.41473 4.42102 4.47950 4.63292 4.64417 4.73292 4.75489 4.79555 4.82068 4.86909 4.90356 4.92891 4.96044 4.98429 5.05789 5.07170 5.09314 5.11444 5.12971 5.17641 5.20119 5.22200 5.26752 5.28316 5.30860 5.32135 5.36585 5.39062 5.42536 5.46533 5.49369 5.53279 5.62525 5.64723 5.65326 5.75348 5.78832 5.92661 6.37964 6.42355 6.73903 6.96658 7.48033 8.94611 Molecular weight = 309.27amu Principal moments of inertia in cm(-1) A = 0.006088 B = 0.005505 C = 0.003052 Principal moments of inertia in units of 10**(-40)*gram-cm**2 A = 4597.854094 B = 5084.647690 C = 9171.006496 Net atomic charges, atomic populations, and dipole contributions calculated with CM2 Atom NO. Type Charge No. of electrons 1 C -0.148 4.148 2 C -0.131 4.131 3 C -0.109 4.109 4 C -0.122 4.122 5 C -0.131 4.131 6 C 0.086 3.914 7 C -0.093 4.093 8 C -0.113 4.113 9 C -0.111 4.111 10 C -0.129 4.129 11 C -0.149 4.149 12 N -0.678 5.678 13 C 0.117 3.883 14 H 0.033 0.967 15 C -0.081 4.081 16 C -0.015 4.015 17 C -0.494 4.494 18 C 0.021 3.979 19 N -0.933 5.933 20 Si 0.940 3.060 21 H -0.267 1.267 22 H -0.293 1.293 23 H -0.268 1.268 24 C 0.011 3.989 25 C -0.095 4.095 26 H 0.053 0.947 27 H 0.057 0.943 28 H 0.060 0.940 29 H 0.053 0.947 30 H 0.049 0.951 31 H 0.065 0.935 32 H 0.068 0.932 33 H 0.046 0.954 34 H 0.046 0.954 35 H 0.071 0.929 36 H 0.078 0.922 37 H 0.055 0.945 38 H 0.070 0.930 39 H 0.063 0.937 40 H 0.055 0.945 41 H 0.044 0.956 42 H 0.066 0.934 43 H 0.059 0.941 44 H 0.048 0.952 45 H 0.056 0.944 46 H 0.054 0.946 47 H 0.060 0.940 48 H 0.055 0.945 49 H 0.339 0.661 50 H 0.049 0.951 51 H 0.038 0.962 52 H 0.022 0.978 53 H 0.017 0.983 54 H 0.262 0.738 55 H -0.005 1.005 56 H 0.022 0.978 57 H 0.045 0.955 58 H 0.038 0.962 For the charges calculated by CM2: Dipole moment (debyes) X Y Z Total from point charges -10.878 17.474 -7.081 21.767 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.206 4.206 2 C -0.169 4.169 3 C -0.147 4.147 4 C -0.159 4.159 5 C -0.170 4.170 6 C -0.024 4.024 7 C -0.131 4.131 8 C -0.150 4.150 9 C -0.148 4.148 10 C -0.167 4.167 11 C -0.206 4.206 12 N -0.310 5.310 13 C 0.008 3.992 14 H 0.051 0.949 15 C -0.120 4.120 16 C -0.053 4.053 17 C -0.513 4.513 18 C -0.177 4.177 19 N -0.573 5.573 20 Si 0.767 3.233 21 H -0.193 1.193 22 H -0.220 1.220 23 H -0.193 1.193 24 C -0.027 4.027 25 C -0.134 4.134 26 H 0.072 0.928 27 H 0.076 0.924 28 H 0.079 0.921 29 H 0.071 0.929 30 H 0.067 0.933 31 H 0.084 0.916 32 H 0.087 0.913 33 H 0.065 0.935 34 H 0.065 0.935 35 H 0.090 0.910 36 H 0.096 0.904 37 H 0.073 0.927 38 H 0.089 0.911 39 H 0.081 0.919 40 H 0.073 0.927 41 H 0.063 0.937 42 H 0.084 0.916 43 H 0.078 0.922 44 H 0.067 0.933 45 H 0.074 0.926 46 H 0.073 0.927 47 H 0.079 0.921 48 H 0.074 0.926 49 H 0.158 0.842 50 H 0.067 0.933 51 H 0.057 0.943 52 H 0.040 0.960 53 H 0.035 0.965 54 H 0.071 0.929 55 H 0.014 0.986 56 H 0.041 0.959 57 H 0.064 0.936 58 H 0.057 0.943 Dipole moment (debyes) X Y Z Total from point charges -10.160 17.606 -7.089 21.528 hybrid contribution 0.330 -1.357 -0.084 1.399 sum -9.830 16.250 -7.172 20.301 Atomic orbital electron populations 1.21747 0.94813 1.01825 1.02174 1.21568 0.96750 0.97234 1.01306 1.21546 0.94340 0.96434 1.02359 1.21653 0.96220 0.97120 1.00936 1.22168 0.95431 0.96335 1.03076 1.21693 0.92194 0.95826 0.92680 1.21481 0.95310 0.94842 1.01516 1.21491 0.96316 0.96485 1.00747 1.21557 0.95041 0.96609 1.01629 1.21556 0.96365 0.97558 1.01210 1.21756 1.00851 0.96003 1.02002 1.60585 1.46758 1.16629 1.06984 1.21152 0.95541 0.87031 0.95496 0.94944 1.21401 0.96425 0.95857 0.98275 1.19687 0.96028 0.94241 0.95383 1.20644 1.00526 1.21609 1.08568 1.25578 1.01072 0.95132 0.95914 1.69967 1.33741 1.27782 1.25840 0.85843 0.80127 0.77937 0.79386 1.19254 1.21985 1.19267 1.19291 0.94172 0.92871 0.96407 1.21572 0.98478 0.96918 0.96406 0.92751 0.92353 0.92103 0.92862 0.93253 0.91584 0.91291 0.93500 0.93516 0.91025 0.90379 0.92650 0.91123 0.91858 0.92664 0.93700 0.91570 0.92237 0.93342 0.92564 0.92728 0.92134 0.92627 0.84241 0.93276 0.94311 0.96001 0.96507 0.92912 0.98621 0.95905 0.93625 0.94329 Number of geometries 1 Number of calculations of the screened coulomb radii 1 The total number of SCF iterations 15. 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.15 -2.46 9.91 71.98 0.71 -1.75 16 2 C -0.13 -2.68 5.93 30.60 0.18 -2.50 16 3 C -0.11 -2.08 5.21 30.59 0.16 -1.92 16 4 C -0.12 -2.77 3.90 30.59 0.12 -2.65 16 5 C -0.13 -2.60 4.53 30.60 0.14 -2.46 16 6 C 0.09 2.06 1.95 44.92 0.09 2.14 16 7 C -0.09 -1.95 4.89 30.59 0.15 -1.80 16 8 C -0.11 -2.66 5.14 30.59 0.16 -2.51 16 9 C -0.11 -2.20 5.21 30.59 0.16 -2.04 16 10 C -0.13 -2.64 5.93 30.59 0.18 -2.46 16 11 C -0.15 -2.49 9.91 71.98 0.71 -1.77 16 12 N -0.68 -19.66 5.29 -469.61 -2.49 -22.15 16 13 C 0.12 4.26 1.45 44.89 0.07 4.33 16 14 H 0.03 1.22 6.04 -2.39 -0.01 1.21 16 15 C -0.08 -3.27 5.41 30.66 0.17 -3.11 16 16 C -0.02 -0.72 5.59 30.10 0.17 -0.55 16 17 C -0.49 -26.65 5.51 -16.17 -0.09 -26.74 16 18 C 0.02 1.18 5.46 84.66 0.46 1.64 16 19 N -0.93 -56.21 12.73 21.45 0.27 -55.93 16 20 Si 0.94 44.94 26.80 68.60 1.84 46.77 16 21 H -0.27 -12.37 7.11 99.48 0.71 -11.67 16 22 H -0.29 -13.91 7.11 99.48 0.71 -13.21 16 23 H -0.27 -12.31 7.11 99.48 0.71 -11.60 16 24 C 0.01 0.54 5.97 30.11 0.18 0.72 16 25 C -0.10 -4.01 5.25 30.66 0.16 -3.85 16 26 H 0.05 0.88 8.14 -2.39 -0.02 0.86 16 27 H 0.06 0.90 8.14 -2.38 -0.02 0.88 16 28 H 0.06 0.90 8.14 -2.39 -0.02 0.88 16 29 H 0.05 1.16 8.14 -2.39 -0.02 1.14 16 30 H 0.05 1.13 8.14 -2.39 -0.02 1.11 16 31 H 0.07 1.17 8.14 -2.39 -0.02 1.15 16 32 H 0.07 1.17 8.14 -2.38 -0.02 1.15 16 33 H 0.05 1.20 7.90 -2.39 -0.02 1.18 16 34 H 0.05 1.25 6.83 -2.39 -0.02 1.23 16 35 H 0.07 1.32 8.04 -2.39 -0.02 1.30 16 36 H 0.08 1.33 8.14 -2.39 -0.02 1.31 16 37 H 0.06 1.52 7.05 -2.39 -0.02 1.51 16 38 H 0.07 1.24 8.14 -2.39 -0.02 1.22 16 39 H 0.06 1.30 8.14 -2.39 -0.02 1.28 16 40 H 0.05 1.61 7.65 -2.39 -0.02 1.59 16 41 H 0.04 1.11 8.14 -2.39 -0.02 1.09 16 42 H 0.07 1.15 8.14 -2.39 -0.02 1.13 16 43 H 0.06 1.17 8.14 -2.39 -0.02 1.15 16 44 H 0.05 1.12 8.14 -2.38 -0.02 1.10 16 45 H 0.06 1.16 8.14 -2.39 -0.02 1.14 16 46 H 0.05 0.87 8.14 -2.39 -0.02 0.85 16 47 H 0.06 0.89 8.14 -2.39 -0.02 0.87 16 48 H 0.05 0.89 8.14 -2.39 -0.02 0.87 16 49 H 0.34 9.26 7.60 -89.69 -0.68 8.58 16 50 H 0.05 1.78 7.50 -2.39 -0.02 1.76 16 51 H 0.04 1.68 8.14 -2.39 -0.02 1.66 16 52 H 0.02 1.00 5.17 -2.39 -0.01 0.99 16 53 H 0.02 0.76 8.14 -2.38 -0.02 0.74 16 54 H 0.26 15.06 8.31 -92.70 -0.77 14.29 16 55 H 0.00 -0.24 8.14 -2.38 -0.02 -0.26 16 56 H 0.02 1.29 7.07 -2.39 -0.02 1.27 16 57 H 0.05 1.70 8.05 -2.39 -0.02 1.68 16 58 H 0.04 1.74 8.14 -2.39 -0.02 1.72 16 Total: -1.00 -58.99 429.67 3.57 -55.41 By element: Atomic # 1 Polarization: 23.08 SS G_CDS: 0.07 Total: 23.15 kcal Atomic # 6 Polarization: -51.13 SS G_CDS: 3.87 Total: -47.26 kcal Atomic # 7 Polarization: -75.87 SS G_CDS: -2.21 Total: -78.08 kcal Atomic # 14 Polarization: 44.94 SS G_CDS: 1.84 Total: 46.77 kcal Total: -58.99 3.57 -55.41 kcal The number of atoms in the molecule is 58 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. ZINC001257550451.mol2 58 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 -34.021 kcal (2) G-P(sol) polarization free energy of solvation -58.987 kcal (3) G-ENP(sol) elect.-nuc.-pol. free energy of system -93.008 kcal (4) G-CDS(sol) cavity-dispersion-solvent structure free energy 3.573 kcal (5) G-P-CDS(sol) = G-P(sol) + G-CDS(sol) = (2) + (4) -55.415 kcal (6) G-S(sol) free energy of system = (1) + (5) -89.435 kcal 1SCF WAS SPECIFIED, SO GEOMETRY OPTIMIZATION WAS NOT USED Total computer time = 1.16 seconds