Wall clock time and date at job start Sun Mar 7 2021 05:27:12 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= 0 * 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 N 1.46499 * 1 3 3 C 1.46501 * 120.00138 * 2 1 4 4 C 1.53001 * 109.47440 * 270.00373 * 3 2 1 5 5 H 1.09003 * 112.84663 * 316.19466 * 4 3 2 6 6 C 1.53778 * 113.61342 * 87.49804 * 4 3 2 7 7 C 1.53781 * 87.07900 * 139.98549 * 6 4 3 8 8 H 1.08999 * 113.62044 * 89.11449 * 7 6 4 9 9 O 1.42906 * 113.61538 * 220.02289 * 7 6 4 10 10 C 1.42906 * 113.99380 * 247.49826 * 9 7 6 11 11 C 1.52991 * 109.47361 * 59.99817 * 10 9 7 12 12 C 1.53004 * 109.46444 * 179.97438 * 10 9 7 13 13 C 1.52999 * 109.47109 * 299.99080 * 10 9 7 14 14 C 1.53775 * 113.61656 * 184.99979 * 4 3 2 15 15 C 1.34780 * 120.00328 * 180.02562 * 2 1 3 16 16 O 1.21278 * 120.00440 * 179.97438 * 15 2 1 17 17 C 1.50700 * 119.99560 * 359.97438 * 15 2 1 18 18 O 1.42904 * 109.46825 * 180.02562 * 17 15 2 19 19 C 1.35919 * 116.99538 * 180.02562 * 18 17 15 20 20 C 1.38798 * 120.05602 * 179.97438 * 19 18 17 21 21 Br 1.89097 * 120.03518 * 359.97438 * 20 19 18 22 22 C 1.38232 * 119.92934 * 179.97438 * 20 19 18 23 23 C 1.38248 * 120.04418 * 359.97438 * 22 20 19 24 24 C 1.50693 * 119.94191 * 179.97438 * 23 22 20 25 25 O 1.42907 * 109.47661 * 300.00763 * 24 23 22 26 26 Si 1.86298 * 109.47618 * 180.02562 * 24 23 22 27 27 C 1.38285 * 120.12000 * 0.26272 * 23 22 20 28 28 C 1.38147 * 120.06930 * 359.48249 * 27 23 22 29 29 H 1.09003 * 109.46857 * 269.99698 * 1 2 3 30 30 H 1.08995 * 109.47144 * 29.99496 * 1 2 3 31 31 H 1.08996 * 109.46826 * 150.00120 * 1 2 3 32 32 H 1.09003 * 109.46925 * 30.00068 * 3 2 1 33 33 H 1.08996 * 109.47323 * 150.00008 * 3 2 1 34 34 H 1.08997 * 113.61628 * 254.52928 * 6 4 3 35 35 H 1.08999 * 113.61392 * 25.44004 * 6 4 3 36 36 H 1.08997 * 109.47758 * 59.99920 * 11 10 9 37 37 H 1.09000 * 109.46861 * 179.97438 * 11 10 9 38 38 H 1.09003 * 109.47001 * 299.99403 * 11 10 9 39 39 H 1.08998 * 109.46951 * 60.00043 * 12 10 9 40 40 H 1.09000 * 109.46823 * 179.97438 * 12 10 9 41 41 H 1.08991 * 109.47521 * 299.99461 * 12 10 9 42 42 H 1.09003 * 109.47075 * 60.00706 * 13 10 9 43 43 H 1.09000 * 109.46907 * 180.02562 * 13 10 9 44 44 H 1.09002 * 109.47233 * 300.00818 * 13 10 9 45 45 H 1.09001 * 113.61382 * 334.56678 * 14 4 3 46 46 H 1.09000 * 113.61789 * 105.38549 * 14 4 3 47 47 H 1.08999 * 109.47327 * 60.00151 * 17 15 2 48 48 H 1.08996 * 109.47015 * 299.99727 * 17 15 2 49 49 H 1.07996 * 119.97973 * 180.02562 * 22 20 19 50 50 H 1.09000 * 109.46936 * 60.00386 * 24 23 22 51 51 H 0.96697 * 114.00527 * 179.97438 * 25 24 23 52 52 H 1.48506 * 109.99834 * 298.68714 * 26 24 23 53 53 H 1.48502 * 110.00011 * 59.99978 * 26 24 23 54 54 H 1.07993 * 119.96404 * 179.72002 * 27 23 22 55 55 H 1.08000 * 120.02340 * 180.26314 * 28 27 23 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.4650 0.0000 0.0000 3 6 2.1975 1.2687 0.0000 4 6 2.4525 1.7105 -1.4425 5 1 2.7834 0.8931 -2.0832 6 6 1.3071 2.5420 -2.0436 7 6 2.3902 3.4193 -2.6934 8 1 2.7033 3.0697 -3.6772 9 8 2.0861 4.8155 -2.6700 10 6 1.8374 5.3777 -3.9600 11 6 0.6460 4.6669 -4.6051 12 6 1.5238 6.8682 -3.8147 13 6 3.0755 5.2001 -4.8411 14 6 3.3232 2.9730 -1.5553 15 6 2.1390 -1.1672 -0.0005 16 8 3.3517 -1.1672 -0.0010 17 6 1.3854 -2.4723 -0.0005 18 8 2.3141 -3.5584 -0.0005 19 6 1.7945 -4.8144 -0.0011 20 6 2.6389 -5.9160 -0.0016 21 35 4.5140 -5.6713 -0.0024 22 6 2.1077 -7.1922 -0.0016 23 6 0.7368 -7.3713 -0.0006 24 6 0.1602 -8.7635 -0.0012 25 8 0.5993 -9.4606 -1.1688 26 14 -1.7002 -8.6653 -0.0002 27 6 -0.1062 -6.2751 -0.0049 28 6 0.4195 -4.9976 0.0004 29 1 -0.3633 -0.0001 1.0277 30 1 -0.3633 0.8900 -0.5137 31 1 -0.3633 -0.8900 -0.5138 32 1 1.6084 2.0284 0.5139 33 1 3.1499 1.1384 0.5138 34 1 0.7001 1.9930 -2.7635 35 1 0.7100 3.0664 -1.2976 36 1 -0.2361 4.7935 -3.9774 37 1 0.4561 5.0961 -5.5888 38 1 0.8690 3.6050 -4.7085 39 1 2.3724 7.3748 -3.3550 40 1 1.3345 7.2972 -4.7987 41 1 0.6417 6.9945 -3.1872 42 1 3.2985 4.1381 -4.9446 43 1 2.8859 5.6293 -5.8249 44 1 3.9243 5.7065 -4.3815 45 1 3.2773 3.6151 -0.6758 46 1 4.3442 2.7717 -1.8796 47 1 0.7584 -2.5298 -0.8903 48 1 0.7592 -2.5302 0.8897 49 1 2.7639 -8.0499 -0.0017 50 1 0.4967 -9.2965 0.8881 51 1 0.2660 -10.3659 -1.2352 52 1 -2.1715 -7.9703 1.2247 53 1 -2.1712 -7.9410 -1.2080 54 1 -1.1766 -6.4179 -0.0086 55 1 -0.2397 -4.1420 0.0013 There are 72 doubly filled levels RHF calculation, no. of doubly occupied orbitals= 72 No. of singly occupied orbitals= 1 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). Br: (AM1): M.J.S. DEWAR AND E. G. ZOEBISCH, THEOCHEM, 180, 1 (1988). 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=0 AM1 1SCF TLIMIT=15 GEO-OK SM5.42R & SOLVNT=WATER s_22_15093190_6201210.mol2 55 1SCF WAS SPECIFIED, SO GEOMETRY OPTIMIZATION WAS NOT USED SCF FIELD WAS ACHIEVED Wall clock time and date at job start Sun Mar 7 2021 05:27:12 Heat of formation + Delta-G solvation = -78.168758 kcal Electronic energy + Delta-G solvation = -34986.699514 eV Core-core repulsion = 30231.784569 eV Total energy + Delta-G solvation = -4754.914945 eV No. of doubly occupied orbitals = 72 and no. of open levels = 1 Molecular weight (most abundant/longest-lived isotopes) = 442.123 amu Computer time = 1.72 seconds Orbital eigenvalues (eV) -41.87661 -41.21823 -40.53218 -39.20085 -38.41555 -37.27897 -36.46529 -33.91435 -33.76425 -32.16824 -31.09015 -30.40905 -28.96969 -28.41200 -28.32248 -27.22252 -26.58208 -26.27761 -24.35060 -23.55686 -22.33265 -21.03323 -20.63861 -20.34566 -19.78049 -18.44384 -18.13462 -17.84256 -17.74066 -17.56861 -17.02373 -16.74801 -16.71924 -16.03648 -15.95914 -15.88964 -15.71443 -15.51458 -14.88709 -14.69763 -14.62952 -14.59229 -14.25286 -14.11807 -13.95790 -13.85779 -13.61746 -13.57527 -13.50708 -13.37992 -13.27745 -13.22178 -12.89576 -12.81780 -12.69030 -12.64044 -12.35669 -12.23639 -12.21690 -12.06652 -12.03259 -11.91550 -11.83680 -11.39459 -11.29754 -11.25336 -10.98094 -10.38277 -10.32714 -10.09705 -9.78892 -9.43059 -4.89862 -0.23587 -0.01957 0.51549 1.27473 1.32715 1.35695 1.68174 1.96785 2.22485 2.27711 2.70985 3.19753 3.32462 3.47538 3.51909 3.53927 3.59283 3.75844 3.77128 3.79512 3.88580 3.89105 3.91019 3.94334 3.96082 4.07875 4.12793 4.19070 4.21629 4.27192 4.42486 4.52290 4.53017 4.58111 4.59442 4.64387 4.69939 4.73292 4.75948 4.78232 4.78867 4.79643 4.83955 4.87473 4.89763 4.91589 4.99358 5.07798 5.09553 5.10951 5.28737 5.29894 5.31972 5.33938 5.46857 5.90756 6.18230 6.39576 6.91721 6.95455 Net atomic charges, atomic populations, and dipole contributions calculated with CM2 Atom NO. Type Charge No. of electrons 1 C 0.061 3.939 2 N -0.604 5.604 3 C 0.130 3.870 4 C -0.117 4.117 5 H 0.076 0.924 6 C -0.134 4.134 7 C 0.059 3.941 8 H 0.089 0.911 9 O -0.378 6.378 10 C 0.098 3.902 11 C -0.171 4.171 12 C -0.138 4.138 13 C -0.178 4.178 14 C -0.132 4.132 15 C 0.508 3.492 16 O -0.565 6.565 17 C 0.034 3.966 18 O -0.301 6.301 19 C 0.148 3.852 20 C -0.140 4.140 21 Br -0.044 7.044 22 C -0.030 4.030 23 C -0.119 4.119 24 C -0.081 4.081 25 O -0.566 6.566 26 Si 0.776 3.224 27 C -0.076 4.076 28 C -0.190 4.190 29 H 0.078 0.922 30 H 0.091 0.909 31 H 0.091 0.909 32 H 0.105 0.895 33 H 0.064 0.936 34 H 0.075 0.925 35 H 0.097 0.903 36 H 0.062 0.938 37 H 0.074 0.926 38 H 0.058 0.942 39 H 0.056 0.944 40 H 0.083 0.917 41 H 0.061 0.939 42 H 0.058 0.942 43 H 0.077 0.923 44 H 0.057 0.943 45 H 0.092 0.908 46 H 0.065 0.935 47 H 0.125 0.875 48 H 0.127 0.873 49 H 0.148 0.852 50 H 0.094 0.906 51 H 0.404 0.596 52 H -0.279 1.279 53 H -0.290 1.290 54 H 0.140 0.860 55 H 0.175 0.825 For the charges calculated by CM2: Dipole moment (debyes) X Y Z Total from point charges -10.960 -1.597 -0.377 11.082 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.081 4.081 2 N -0.338 5.338 3 C 0.007 3.993 4 C -0.136 4.136 5 H 0.094 0.906 6 C -0.171 4.171 7 C 0.001 3.999 8 H 0.107 0.893 9 O -0.298 6.298 10 C 0.059 3.941 11 C -0.229 4.229 12 C -0.195 4.195 13 C -0.235 4.235 14 C -0.170 4.170 15 C 0.296 3.704 16 O -0.443 6.443 17 C -0.043 4.043 18 O -0.214 6.214 19 C 0.100 3.900 20 C -0.220 4.220 21 Br 0.042 6.958 22 C -0.049 4.049 23 C -0.120 4.120 24 C -0.190 4.190 25 O -0.376 6.376 26 Si 0.680 3.320 27 C -0.093 4.093 28 C -0.209 4.209 29 H 0.096 0.904 30 H 0.109 0.891 31 H 0.109 0.891 32 H 0.124 0.876 33 H 0.083 0.917 34 H 0.093 0.907 35 H 0.115 0.885 36 H 0.081 0.919 37 H 0.093 0.907 38 H 0.077 0.923 39 H 0.075 0.925 40 H 0.102 0.898 41 H 0.080 0.920 42 H 0.077 0.923 43 H 0.096 0.904 44 H 0.076 0.924 45 H 0.111 0.889 46 H 0.084 0.916 47 H 0.143 0.857 48 H 0.144 0.856 49 H 0.165 0.835 50 H 0.112 0.888 51 H 0.256 0.744 52 H -0.206 1.206 53 H -0.217 1.217 54 H 0.157 0.843 55 H 0.192 0.808 Dipole moment (debyes) X Y Z Total from point charges -9.820 -0.235 -0.303 9.828 hybrid contribution 0.541 0.194 -0.196 0.607 sum -9.279 -0.042 -0.499 9.292 Atomic orbital electron populations 1.22262 0.76449 1.05418 1.03930 1.47815 1.08705 1.04620 1.72622 1.21501 0.97525 0.84336 0.95900 1.22548 0.98035 0.96085 0.96952 0.90583 1.23528 0.96400 0.96625 1.00559 1.23645 0.96535 0.81311 0.98394 0.89294 1.88014 1.89812 1.22887 1.29076 1.21674 0.94386 0.94679 0.83362 1.22416 0.98090 1.00692 1.01656 1.21783 1.02151 0.92500 1.03115 1.22461 0.98044 1.02251 1.00787 1.23233 0.98185 0.96075 0.99493 1.20443 0.87592 0.85893 0.76482 1.90773 1.15849 1.87601 1.50124 1.22900 0.91636 0.83080 1.06729 1.86254 1.35421 1.11719 1.87962 1.18747 0.92071 0.83875 0.95309 1.21328 0.91015 0.94460 1.15228 1.96758 1.03451 1.97850 1.97775 1.20943 0.95436 0.94466 0.94098 1.19672 0.92963 0.95309 1.04090 1.24889 1.12048 0.90425 0.91641 1.86447 1.76677 1.32379 1.42126 0.94025 0.72929 0.88221 0.76812 1.21047 0.99613 0.92139 0.96542 1.21436 0.92421 0.99629 1.07401 0.90388 0.89062 0.89087 0.87639 0.91732 0.90655 0.88502 0.91853 0.90742 0.92334 0.92488 0.89829 0.92008 0.92335 0.90386 0.92435 0.88899 0.91631 0.85716 0.85591 0.83453 0.88796 0.74446 1.20612 1.21747 0.84299 0.80782 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.06 -0.64 10.24 127.77 1.31 0.67 16 2 N -0.60 1.00 2.87 -846.77 -2.43 -1.44 16 3 C 0.13 -0.03 5.24 86.38 0.45 0.43 16 4 C -0.12 -0.02 3.96 -10.27 -0.04 -0.06 16 5 H 0.08 0.32 8.14 -2.39 -0.02 0.30 16 6 C -0.13 0.35 7.04 31.12 0.22 0.57 16 7 C 0.06 0.05 2.82 31.12 0.09 0.14 16 8 H 0.09 0.05 5.55 -2.39 -0.01 0.04 16 9 O -0.38 -1.41 10.00 -148.98 -1.49 -2.90 16 10 C 0.10 0.06 1.11 -10.80 -0.01 0.04 16 11 C -0.17 0.32 8.38 71.98 0.60 0.92 16 12 C -0.14 0.07 8.84 71.98 0.64 0.71 16 13 C -0.18 0.03 8.38 71.98 0.60 0.64 16 14 C -0.13 -0.21 7.47 31.12 0.23 0.02 16 15 C 0.51 2.56 7.75 87.66 0.68 3.24 16 16 O -0.57 -9.24 16.08 16.02 0.26 -8.98 16 17 C 0.03 -0.03 3.77 71.24 0.27 0.24 16 18 O -0.30 -2.45 9.01 -73.60 -0.66 -3.12 16 19 C 0.15 0.94 6.68 22.52 0.15 1.09 16 20 C -0.14 -1.23 6.24 22.40 0.14 -1.09 16 21 Br -0.04 -0.49 33.16 -20.37 -0.68 -1.16 16 22 C -0.03 -0.18 9.35 22.27 0.21 0.03 16 23 C -0.12 -0.40 5.42 -19.85 -0.11 -0.51 16 24 C -0.08 0.19 2.81 71.23 0.20 0.39 16 25 O -0.57 0.78 12.74 -148.98 -1.90 -1.12 16 26 Si 0.78 -4.06 29.78 68.60 2.04 -2.02 16 27 C -0.08 -0.12 8.66 22.26 0.19 0.07 16 28 C -0.19 -0.08 9.04 22.36 0.20 0.13 16 29 H 0.08 -0.98 8.14 -2.39 -0.02 -1.00 16 30 H 0.09 -1.09 7.64 -2.39 -0.02 -1.11 16 31 H 0.09 -1.30 6.44 -2.39 -0.02 -1.32 16 32 H 0.11 -0.49 7.77 -2.39 -0.02 -0.50 16 33 H 0.06 0.34 7.42 -2.39 -0.02 0.32 16 34 H 0.07 -0.30 8.14 -2.39 -0.02 -0.32 16 35 H 0.10 -0.34 7.90 -2.39 -0.02 -0.36 16 36 H 0.06 -0.09 8.14 -2.39 -0.02 -0.11 16 37 H 0.07 -0.26 8.14 -2.39 -0.02 -0.28 16 38 H 0.06 -0.11 7.02 -2.39 -0.02 -0.13 16 39 H 0.06 0.05 8.14 -2.39 -0.02 0.03 16 40 H 0.08 -0.28 8.14 -2.39 -0.02 -0.30 16 41 H 0.06 -0.01 8.14 -2.39 -0.02 -0.03 16 42 H 0.06 -0.01 6.25 -2.39 -0.01 -0.02 16 43 H 0.08 -0.18 8.14 -2.39 -0.02 -0.20 16 44 H 0.06 0.09 8.14 -2.39 -0.02 0.07 16 45 H 0.09 0.15 8.10 -2.39 -0.02 0.13 16 46 H 0.06 0.20 8.14 -2.39 -0.02 0.19 16 47 H 0.13 -0.93 6.92 -2.39 -0.02 -0.94 16 48 H 0.13 -0.95 7.63 -2.39 -0.02 -0.97 16 49 H 0.15 0.66 8.06 -2.91 -0.02 0.64 16 50 H 0.09 -0.56 7.96 -2.39 -0.02 -0.58 16 51 H 0.40 -4.80 8.74 -74.06 -0.65 -5.45 16 52 H -0.28 0.29 7.11 99.48 0.71 1.00 16 53 H -0.29 -0.42 7.11 99.48 0.71 0.28 16 54 H 0.14 -0.20 4.87 -2.91 -0.01 -0.21 16 55 H 0.17 -1.35 6.30 -2.91 -0.02 -1.37 16 Total: 0.00 -26.72 455.20 1.46 -25.26 By element: Atomic # 1 Polarization: -12.48 SS G_CDS: 0.29 Total: -12.19 kcal Atomic # 6 Polarization: 1.64 SS G_CDS: 6.03 Total: 7.66 kcal Atomic # 7 Polarization: 1.00 SS G_CDS: -2.43 Total: -1.44 kcal Atomic # 8 Polarization: -12.33 SS G_CDS: -3.79 Total: -16.12 kcal Atomic # 14 Polarization: -4.06 SS G_CDS: 2.04 Total: -2.02 kcal Atomic # 35 Polarization: -0.49 SS G_CDS: -0.68 Total: -1.16 kcal Total: -26.72 1.46 -25.26 kcal The number of atoms in the molecule is 55 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. s_22_15093190_6201210.mol2 55 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 -52.904 kcal (2) G-P(sol) polarization free energy of solvation -26.721 kcal (3) G-ENP(sol) elect.-nuc.-pol. free energy of system -79.625 kcal (4) G-CDS(sol) cavity-dispersion-solvent structure free energy 1.457 kcal (5) G-P-CDS(sol) = G-P(sol) + G-CDS(sol) = (2) + (4) -25.265 kcal (6) G-S(sol) free energy of system = (1) + (5) -78.169 kcal 1SCF WAS SPECIFIED, SO GEOMETRY OPTIMIZATION WAS NOT USED Total computer time = 1.72 seconds