Wall clock time and date at job start Sat Apr 11 2020 06:12:03 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 C 1.50703 * 1 3 3 O 1.21284 * 119.99680 * 2 1 4 4 N 1.34775 * 119.99854 * 180.02562 * 2 1 3 5 5 C 1.47700 * 120.47658 * 359.97438 * 4 2 1 6 6 C 1.51058 * 110.42732 * 141.27666 * 5 4 2 7 7 O 1.42953 * 108.84616 * 89.24929 * 6 5 4 8 8 C 1.37218 * 118.04735 * 277.44570 * 7 6 5 9 9 C 1.54043 * 111.61885 * 63.02834 * 8 7 6 10 10 C 1.53255 * 108.57291 * 164.84660 * 9 8 7 11 11 N 1.46980 * 109.25370 * 184.45927 * 10 9 8 12 12 C 1.46893 * 111.00364 * 182.50485 * 11 10 9 13 13 H 1.09006 * 109.47121 * 55.56793 * 12 11 10 14 14 C 1.52995 * 109.47564 * 175.57280 * 12 11 10 15 15 C 1.53001 * 109.47221 * 295.57664 * 12 11 10 16 16 C 1.52998 * 109.47215 * 176.32664 * 15 12 11 17 17 C 1.50705 * 109.47151 * 185.20205 * 16 15 12 18 18 H 1.08002 * 119.99783 * 0.02562 * 17 16 15 19 19 C 1.37877 * 120.00191 * 179.97438 * 17 16 15 20 20 N 1.31504 * 120.00069 * 180.02562 * 19 17 16 21 21 Si 1.86807 * 119.99868 * 0.02562 * 19 17 16 22 22 H 1.48496 * 109.46982 * 89.99760 * 21 19 17 23 23 H 1.48496 * 109.46891 * 210.00047 * 21 19 17 24 24 H 1.48500 * 109.46830 * 329.99965 * 21 19 17 25 25 C 1.47012 * 110.77382 * 58.71516 * 11 10 9 26 26 C 1.53103 * 109.27078 * 301.56024 * 25 11 10 27 27 O 1.43026 * 109.27150 * 56.76806 * 26 25 11 28 28 C 1.48258 * 120.46865 * 180.02562 * 4 2 1 29 29 H 1.08998 * 109.46600 * 240.00118 * 1 2 3 30 30 H 1.09000 * 109.47132 * 359.97438 * 1 2 3 31 31 H 1.08998 * 109.47391 * 120.00087 * 1 2 3 32 32 H 1.09007 * 109.30877 * 20.98424 * 5 4 2 33 33 H 1.08999 * 109.31318 * 261.57633 * 5 4 2 34 34 H 1.09001 * 109.57756 * 329.44230 * 6 5 4 35 35 H 1.09000 * 109.58113 * 209.05768 * 6 5 4 36 36 H 1.09004 * 109.05000 * 302.48052 * 8 7 6 37 37 H 1.08999 * 108.99550 * 183.54120 * 8 7 6 38 38 H 1.08996 * 109.50741 * 64.52408 * 10 9 8 39 39 H 1.09000 * 109.51055 * 304.39146 * 10 9 8 40 40 H 1.09001 * 109.47495 * 180.02562 * 14 12 11 41 41 H 1.09003 * 109.46847 * 59.99871 * 14 12 11 42 42 H 1.09006 * 109.46923 * 300.00237 * 14 12 11 43 43 H 1.08997 * 109.47012 * 296.32242 * 15 12 11 44 44 H 1.08998 * 109.47205 * 56.32174 * 15 12 11 45 45 H 1.09001 * 109.47328 * 305.20178 * 16 15 12 46 46 H 1.09001 * 109.47254 * 65.20599 * 16 15 12 47 47 H 0.97008 * 120.00353 * 180.02562 * 20 19 17 48 48 H 1.09000 * 109.50655 * 61.49916 * 25 11 10 49 49 H 1.09000 * 109.50632 * 181.62608 * 25 11 10 50 50 H 1.08997 * 109.50519 * 176.70535 * 26 25 11 51 51 H 1.08998 * 109.50503 * 296.79088 * 26 25 11 52 52 H 1.09002 * 108.99903 * 19.92160 * 28 4 2 53 53 H 1.08990 * 109.00328 * 261.14980 * 28 4 2 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.5070 0.0000 0.0000 3 8 2.1134 1.0504 0.0000 4 7 2.1809 -1.1672 -0.0005 5 6 1.4530 -2.4524 -0.0005 6 6 2.1538 -3.4552 -0.8866 7 8 3.1172 -4.1670 -0.1063 8 6 4.3415 -3.5819 0.0978 9 6 4.2170 -2.2666 0.8900 10 6 5.6060 -1.8736 1.4049 11 7 5.4883 -0.6785 2.2524 12 6 6.8076 -0.2324 2.7196 13 1 7.4599 -0.0671 1.8620 14 6 6.6566 1.0721 3.5045 15 6 7.4176 -1.3046 3.6248 16 6 8.8287 -0.8813 4.0374 17 6 9.4771 -1.9924 4.8225 18 1 8.9368 -2.9089 5.0081 19 6 10.7604 -1.8392 5.3028 20 7 11.3263 -2.8089 5.9874 21 14 11.6952 -0.2542 4.9810 22 1 11.4340 0.7068 6.0825 23 1 13.1500 -0.5432 4.9106 24 1 11.2444 0.3336 3.6939 25 6 4.5856 -0.9295 3.3853 26 6 3.2090 -1.3374 2.8536 27 8 3.3498 -2.4822 2.0079 28 6 3.6634 -1.1794 0.0000 29 1 -0.3632 -0.5138 0.8900 30 1 -0.3633 1.0277 0.0005 31 1 -0.3634 -0.5138 -0.8899 32 1 0.4394 -2.2926 -0.3685 33 1 1.4069 -2.8406 1.0170 34 1 2.6568 -2.9341 -1.7012 35 1 1.4245 -4.1555 -1.2937 36 1 4.8016 -3.3733 -0.8681 37 1 4.9767 -4.2725 0.6525 38 1 6.2597 -1.6576 0.5599 39 1 6.0235 -2.6940 1.9885 40 1 7.6356 1.4033 3.8507 41 1 6.0043 0.9066 4.3620 42 1 6.2221 1.8359 2.8596 43 1 6.7993 -1.4236 4.5145 44 1 7.4650 -2.2510 3.0860 45 1 9.4207 -0.6728 3.1462 46 1 8.7732 0.0153 4.6547 47 1 12.2291 -2.7011 6.3257 48 1 4.9900 -1.7321 4.0020 49 1 4.4908 -0.0229 3.9829 50 1 2.5535 -1.5810 3.6897 51 1 2.7804 -0.5134 2.2832 52 1 4.0261 -0.2136 0.3520 53 1 4.0174 -1.3401 -1.0182 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 ZINC001172495386.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:12:03 Heat of formation + Delta-G solvation = -103.420216 kcal Electronic energy + Delta-G solvation = -32787.855428 eV Core-core repulsion = 28660.408532 eV Total energy + Delta-G solvation = -4127.446896 eV No. of doubly occupied orbitals = 66 Molecular weight (most abundant/longest-lived isotopes) = 340.220 amu Computer time = 0.82 seconds Orbital eigenvalues (eV) -41.30491 -39.39085 -38.72825 -36.84653 -35.35911 -33.59403 -33.17527 -31.77542 -30.48094 -29.93424 -29.36217 -27.78824 -26.55183 -25.36161 -24.77839 -23.23985 -22.73615 -21.14208 -20.31301 -19.40096 -18.71371 -17.59088 -17.35610 -17.13970 -16.75613 -16.14044 -16.07503 -15.67881 -15.27042 -14.97500 -14.64504 -14.49067 -14.31567 -14.28268 -13.93509 -13.71811 -13.32142 -13.22612 -12.93410 -12.83367 -12.61120 -12.52579 -12.44549 -12.11254 -11.96440 -11.84766 -11.72955 -11.58844 -11.39155 -11.13627 -11.05018 -10.79595 -10.55340 -10.48829 -10.18119 -9.95866 -9.93575 -9.79987 -9.69691 -9.51188 -9.21923 -8.75236 -8.55071 -8.00828 -6.05941 -4.03955 2.43055 2.63237 3.28998 3.29998 3.36637 3.39980 3.80363 3.89985 4.27887 4.43451 4.45298 4.50177 4.61434 4.68195 4.82256 4.86932 4.92625 5.02892 5.04899 5.15803 5.17831 5.23727 5.27844 5.31959 5.37150 5.45831 5.49973 5.53385 5.58408 5.59208 5.69784 5.77285 5.78036 5.86076 5.89038 6.02077 6.11515 6.14791 6.20996 6.39626 6.48135 6.57315 6.82260 7.02018 7.12489 7.15516 7.22441 7.26453 7.28250 7.41269 7.78145 7.96584 8.52395 8.98743 9.56672 11.07638 Molecular weight = 340.22amu Principal moments of inertia in cm(-1) A = 0.019016 B = 0.002919 C = 0.002694 Principal moments of inertia in units of 10**(-40)*gram-cm**2 A = 1472.094499 B = 9590.637263 C =10390.254109 Net atomic charges, atomic populations, and dipole contributions calculated with CM2 Atom NO. Type Charge No. of electrons 1 C -0.172 4.172 2 C 0.510 3.490 3 O -0.543 6.543 4 N -0.607 5.607 5 C 0.073 3.927 6 C 0.053 3.947 7 O -0.365 6.365 8 C 0.087 3.913 9 C 0.074 3.926 10 C 0.030 3.970 11 N -0.528 5.528 12 C 0.077 3.923 13 H 0.080 0.920 14 C -0.142 4.142 15 C -0.131 4.131 16 C -0.009 4.009 17 C -0.520 4.520 18 H 0.069 0.931 19 C 0.066 3.934 20 N -0.904 5.904 21 Si 0.879 3.121 22 H -0.277 1.277 23 H -0.288 1.288 24 H -0.265 1.265 25 C 0.013 3.987 26 C 0.061 3.939 27 O -0.354 6.354 28 C 0.089 3.911 29 H 0.087 0.913 30 H 0.090 0.910 31 H 0.086 0.914 32 H 0.084 0.916 33 H 0.095 0.905 34 H 0.058 0.942 35 H 0.106 0.894 36 H 0.055 0.945 37 H 0.106 0.894 38 H 0.085 0.915 39 H 0.061 0.939 40 H 0.066 0.934 41 H 0.052 0.948 42 H 0.054 0.946 43 H 0.054 0.946 44 H 0.054 0.946 45 H 0.014 0.986 46 H 0.014 0.986 47 H 0.263 0.737 48 H 0.052 0.948 49 H 0.086 0.914 50 H 0.092 0.908 51 H 0.057 0.943 52 H 0.113 0.887 53 H 0.064 0.936 For the charges calculated by CM2: Dipole moment (debyes) X Y Z Total from point charges -23.803 1.331 -15.199 28.273 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.232 4.232 2 C 0.299 3.701 3 O -0.420 6.420 4 N -0.343 5.343 5 C -0.049 4.049 6 C -0.024 4.024 7 O -0.282 6.282 8 C 0.008 3.992 9 C 0.032 3.968 10 C -0.099 4.099 11 N -0.252 5.252 12 C -0.031 4.031 13 H 0.098 0.902 14 C -0.199 4.199 15 C -0.170 4.170 16 C -0.046 4.046 17 C -0.558 4.558 18 H 0.087 0.913 19 C -0.135 4.135 20 N -0.543 5.543 21 Si 0.707 3.293 22 H -0.203 1.203 23 H -0.214 1.214 24 H -0.190 1.190 25 C -0.116 4.116 26 C -0.016 4.016 27 O -0.273 6.273 28 C -0.033 4.033 29 H 0.106 0.894 30 H 0.109 0.891 31 H 0.104 0.896 32 H 0.102 0.898 33 H 0.113 0.887 34 H 0.076 0.924 35 H 0.124 0.876 36 H 0.073 0.927 37 H 0.124 0.876 38 H 0.103 0.897 39 H 0.079 0.921 40 H 0.085 0.915 41 H 0.071 0.929 42 H 0.073 0.927 43 H 0.073 0.927 44 H 0.073 0.927 45 H 0.032 0.968 46 H 0.032 0.968 47 H 0.072 0.928 48 H 0.070 0.930 49 H 0.104 0.896 50 H 0.110 0.890 51 H 0.075 0.925 52 H 0.131 0.869 53 H 0.082 0.918 Dipole moment (debyes) X Y Z Total from point charges -23.674 1.885 -15.296 28.249 hybrid contribution 1.256 0.589 0.396 1.443 sum -22.418 2.474 -14.900 27.031 Atomic orbital electron populations 1.21941 0.91941 1.05298 1.03983 1.20741 0.90503 0.82955 0.75935 1.90750 1.70646 1.30896 1.49736 1.48612 1.06269 1.05960 1.73446 1.22311 0.95624 0.83798 1.03162 1.22618 0.92036 0.94898 0.92810 1.86533 1.23973 1.45688 1.72037 1.21968 0.86074 0.93390 0.97737 1.22567 0.90553 0.95931 0.87726 1.22640 0.95324 0.93603 0.98320 1.62506 1.12140 1.31081 1.19453 1.21889 0.87928 0.95481 0.97784 0.90213 1.21862 1.02155 0.96417 0.99512 1.22085 0.98939 0.97524 0.98427 1.19536 0.91894 0.97005 0.96213 1.21158 0.99421 1.01913 1.33357 0.91297 1.25040 0.97499 0.96370 0.94609 1.69916 1.14797 1.32589 1.37038 0.86780 0.78995 0.84187 0.79353 1.20299 1.21408 1.18959 1.22812 0.93395 1.01659 0.93696 1.22680 0.96212 0.89006 0.93708 1.88101 1.66130 1.35923 1.37099 1.21943 0.79406 1.02718 0.99238 0.89439 0.89115 0.89570 0.89761 0.88676 0.92388 0.87559 0.92748 0.87616 0.89688 0.92104 0.91509 0.92869 0.92731 0.92733 0.92707 0.96776 0.96817 0.92755 0.92972 0.89577 0.89022 0.92502 0.86857 0.91825 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.17 -1.09 8.77 36.01 0.32 -0.78 16 2 C 0.51 5.55 7.08 -10.98 -0.08 5.47 16 3 O -0.54 -8.42 16.58 5.55 0.09 -8.33 16 4 N -0.61 -5.92 2.13 -173.60 -0.37 -6.29 16 5 C 0.07 0.53 5.34 -4.36 -0.02 0.50 16 6 C 0.05 0.37 6.41 36.18 0.23 0.60 16 7 O -0.37 -3.62 9.74 -40.26 -0.39 -4.02 16 8 C 0.09 0.86 5.47 37.70 0.21 1.06 16 9 C 0.07 0.83 0.63 -90.94 -0.06 0.77 16 10 C 0.03 0.39 3.88 -3.65 -0.01 0.37 16 11 N -0.53 -7.47 3.88 -226.26 -0.88 -8.34 16 12 C 0.08 1.23 2.45 -67.72 -0.17 1.06 16 13 H 0.08 1.29 8.13 -51.93 -0.42 0.87 16 14 C -0.14 -2.17 8.71 37.15 0.32 -1.85 16 15 C -0.13 -2.54 3.91 -26.73 -0.10 -2.64 16 16 C -0.01 -0.20 3.63 -27.88 -0.10 -0.30 16 17 C -0.52 -14.90 9.66 -34.44 -0.33 -15.23 16 18 H 0.07 2.14 7.99 -52.49 -0.42 1.72 16 19 C 0.07 1.93 5.47 -17.82 -0.10 1.83 16 20 N -0.90 -28.12 13.96 55.43 0.77 -27.34 16 21 Si 0.88 21.83 27.47 -169.99 -4.67 17.16 16 22 H -0.28 -6.76 7.11 56.52 0.40 -6.36 16 23 H -0.29 -7.10 7.11 56.52 0.40 -6.69 16 24 H -0.26 -6.29 6.54 56.52 0.37 -5.92 16 25 C 0.01 0.17 4.86 -3.71 -0.02 0.15 16 26 C 0.06 0.74 6.57 37.21 0.24 0.98 16 27 O -0.35 -4.39 7.49 -51.73 -0.39 -4.78 16 28 C 0.09 1.00 3.87 -4.07 -0.02 0.98 16 29 H 0.09 0.43 8.12 -51.93 -0.42 0.01 16 30 H 0.09 0.60 8.12 -51.93 -0.42 0.17 16 31 H 0.09 0.36 7.48 -51.93 -0.39 -0.03 16 32 H 0.08 0.37 6.22 -51.92 -0.32 0.04 16 33 H 0.10 0.84 6.52 -51.93 -0.34 0.50 16 34 H 0.06 0.40 7.18 -51.93 -0.37 0.02 16 35 H 0.11 0.49 8.14 -51.93 -0.42 0.06 16 36 H 0.05 0.45 7.73 -51.93 -0.40 0.05 16 37 H 0.11 1.08 7.94 -51.93 -0.41 0.66 16 38 H 0.08 1.02 8.13 -51.93 -0.42 0.59 16 39 H 0.06 0.88 6.37 -51.93 -0.33 0.55 16 40 H 0.07 1.08 6.64 -51.93 -0.35 0.73 16 41 H 0.05 0.78 6.37 -51.93 -0.33 0.45 16 42 H 0.05 0.78 8.14 -51.93 -0.42 0.36 16 43 H 0.05 1.08 6.65 -51.93 -0.35 0.74 16 44 H 0.05 1.09 6.39 -51.93 -0.33 0.76 16 45 H 0.01 0.33 7.20 -51.93 -0.37 -0.04 16 46 H 0.01 0.32 6.22 -51.93 -0.32 0.00 16 47 H 0.26 7.64 8.31 -40.82 -0.34 7.30 16 48 H 0.05 0.80 6.50 -51.93 -0.34 0.46 16 49 H 0.09 1.15 6.21 -51.93 -0.32 0.83 16 50 H 0.09 0.97 8.14 -51.93 -0.42 0.55 16 51 H 0.06 0.74 5.99 -51.93 -0.31 0.43 16 52 H 0.11 1.60 6.15 -51.93 -0.32 1.28 16 53 H 0.06 0.65 7.13 -51.93 -0.37 0.28 16 LS Contribution 382.86 15.07 5.77 5.77 Total: -1.00 -34.24 382.86 -8.57 -42.80 By element: Atomic # 1 Polarization: 9.19 SS G_CDS: -8.82 Total: 0.37 kcal Atomic # 6 Polarization: -7.33 SS G_CDS: 0.31 Total: -7.02 kcal Atomic # 7 Polarization: -41.50 SS G_CDS: -0.47 Total: -41.97 kcal Atomic # 8 Polarization: -16.43 SS G_CDS: -0.69 Total: -17.12 kcal Atomic # 14 Polarization: 21.83 SS G_CDS: -4.67 Total: 17.16 kcal Total LS contribution 5.77 Total: 5.77 kcal Total: -34.24 -8.57 -42.80 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. ZINC001172495386.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 -60.616 kcal (2) G-P(sol) polarization free energy of solvation -34.236 kcal (3) G-ENP(sol) elect.-nuc.-pol. free energy of system -94.852 kcal (4) G-CDS(sol) cavity-dispersion-solvent structure free energy -8.568 kcal (5) G-P-CDS(sol) = G-P(sol) + G-CDS(sol) = (2) + (4) -42.804 kcal (6) G-S(sol) free energy of system = (1) + (5) -103.420 kcal 1SCF WAS SPECIFIED, SO GEOMETRY OPTIMIZATION WAS NOT USED Total computer time = 0.82 seconds