Wall clock time and date at job start Sat Apr 11 2020 06:11:27 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.52995 * 1 3 3 H 1.09004 * 110.63963 * 2 1 4 4 C 1.54334 * 110.85086 * 236.83407 * 2 1 3 5 5 N 1.47024 * 107.26646 * 219.67074 * 4 2 1 6 6 C 1.36932 * 125.64836 * 178.85202 * 5 4 2 7 7 H 1.08001 * 120.00172 * 347.84778 * 6 5 4 8 8 C 1.38810 * 119.99990 * 167.85563 * 6 5 4 9 9 N 1.31624 * 119.99867 * 326.24876 * 8 6 5 10 10 Si 1.86800 * 120.00129 * 145.98140 * 8 6 5 11 11 H 1.48506 * 109.47282 * 359.97438 * 10 8 6 12 12 H 1.48494 * 109.47535 * 119.99869 * 10 8 6 13 13 H 1.48503 * 109.47098 * 239.99494 * 10 8 6 14 14 C 1.47422 * 108.70201 * 358.58505 * 5 4 2 15 15 C 1.54913 * 104.83247 * 24.35686 * 14 5 4 16 16 H 1.08995 * 111.00610 * 81.06475 * 15 14 5 17 17 N 1.46497 * 111.00373 * 204.87955 * 15 14 5 18 18 C 1.34777 * 120.00030 * 273.54696 * 17 15 14 19 19 O 1.21279 * 120.00139 * 359.97438 * 18 17 15 20 20 C 1.50698 * 120.00191 * 180.02562 * 18 17 15 21 21 C 1.53001 * 109.47494 * 179.97438 * 20 18 17 22 22 C 1.52992 * 109.47633 * 179.97438 * 21 20 18 23 23 O 1.42901 * 109.47510 * 180.02562 * 22 21 20 24 24 C 1.42894 * 114.00165 * 179.97438 * 23 22 21 25 25 C 1.53175 * 109.54647 * 269.55727 * 24 23 22 26 26 C 1.53093 * 109.16577 * 183.15289 * 25 24 23 27 27 O 1.42906 * 109.41249 * 57.60617 * 26 25 24 28 28 C 1.42894 * 114.09882 * 298.84374 * 27 26 25 29 29 C 1.53087 * 109.41801 * 61.16042 * 28 27 26 30 30 H 1.09000 * 109.47522 * 296.85818 * 1 2 3 31 31 H 1.09001 * 109.47262 * 56.85915 * 1 2 3 32 32 H 1.09001 * 109.47218 * 176.85755 * 1 2 3 33 33 H 1.09007 * 109.99211 * 100.19376 * 4 2 1 34 34 H 1.09002 * 109.77986 * 339.02692 * 4 2 1 35 35 H 0.96994 * 119.99575 * 353.97716 * 9 8 6 36 36 H 1.08998 * 110.36625 * 143.19544 * 14 5 4 37 37 H 1.08997 * 110.36749 * 265.40909 * 14 5 4 38 38 H 0.97006 * 120.00267 * 93.54539 * 17 15 14 39 39 H 1.08998 * 109.47467 * 299.99613 * 20 18 17 40 40 H 1.08997 * 109.47611 * 59.99448 * 20 18 17 41 41 H 1.09008 * 109.46766 * 300.00361 * 21 20 18 42 42 H 1.08997 * 109.46928 * 59.99292 * 21 20 18 43 43 H 1.09006 * 109.47068 * 300.00465 * 22 21 20 44 44 H 1.09001 * 109.47041 * 59.99484 * 22 21 20 45 45 H 1.09004 * 109.54704 * 29.82637 * 24 23 22 46 46 H 1.09006 * 109.52169 * 303.05790 * 25 24 23 47 47 H 1.08995 * 109.52525 * 63.24176 * 25 24 23 48 48 H 1.08994 * 109.48245 * 297.62238 * 26 25 24 49 49 H 1.09000 * 109.48797 * 177.59225 * 26 25 24 50 50 H 1.09007 * 109.47952 * 301.17314 * 28 27 26 51 51 H 1.09007 * 109.48964 * 181.15153 * 28 27 26 52 52 H 1.09004 * 109.52767 * 182.48856 * 29 28 27 53 53 H 1.08999 * 109.52349 * 62.29683 * 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 6 1.5300 0.0000 0.0000 3 1 1.9142 1.0201 0.0000 4 6 2.0793 -0.7890 -1.2073 5 7 3.2445 -1.5519 -0.7364 6 6 4.0043 -2.3994 -1.4977 7 1 3.6582 -2.7053 -2.4739 8 6 5.2194 -2.8653 -1.0147 9 7 5.9612 -2.0954 -0.2469 10 14 5.7985 -4.5875 -1.4484 11 1 4.7870 -5.2468 -2.3131 12 1 7.0924 -4.5071 -2.1726 13 1 5.9761 -5.3781 -0.2039 14 6 3.4633 -1.2485 0.6896 15 6 2.0789 -0.7934 1.2151 16 1 1.4452 -1.6511 1.4403 17 7 2.2214 0.0730 2.3878 18 6 2.3156 -0.4687 3.6184 19 8 2.2819 -1.6733 3.7549 20 6 2.4615 0.4225 4.8248 21 6 2.5460 -0.4378 6.0872 22 6 2.6949 0.4669 7.3120 23 8 2.7744 -0.3366 8.4910 24 6 2.9163 0.4149 9.6981 25 6 4.4035 0.6197 10.0024 26 6 4.5466 1.3566 11.3366 27 8 3.8975 0.6067 12.3655 28 6 2.4955 0.4241 12.1582 29 6 2.2699 -0.3503 10.8571 30 1 -0.3634 0.4643 0.9168 31 1 -0.3634 0.5618 -0.8605 32 1 -0.3634 -1.0261 -0.0563 33 1 2.3814 -0.1003 -1.9964 34 1 1.3156 -1.4714 -1.5804 35 1 5.7025 -1.1740 -0.0896 36 1 3.7961 -2.1404 1.2205 37 1 4.1945 -0.4477 0.7992 38 1 2.2480 1.0365 2.2787 39 1 3.3695 1.0179 4.7291 40 1 1.5986 1.0848 4.8940 41 1 1.6377 -1.0329 6.1828 42 1 3.4088 -1.1003 6.0183 43 1 3.6033 1.0617 7.2161 44 1 1.8323 1.1298 7.3810 45 1 2.4293 1.3834 9.5841 46 1 4.8990 -0.3491 10.0663 47 1 4.8596 1.2114 9.2087 48 1 4.0863 2.3415 11.2592 49 1 5.6035 1.4674 11.5793 50 1 2.0082 1.3969 12.0916 51 1 2.0745 -0.1368 12.9928 52 1 1.2001 -0.4521 10.6745 53 1 2.7223 -1.3386 10.9383 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 ZINC001210193454.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:11:27 Heat of formation + Delta-G solvation = -126.180814 kcal Electronic energy + Delta-G solvation = -29682.952156 eV Core-core repulsion = 25554.518287 eV Total energy + Delta-G solvation = -4128.433869 eV No. of doubly occupied orbitals = 66 Molecular weight (most abundant/longest-lived isotopes) = 340.220 amu Computer time = 0.89 seconds Orbital eigenvalues (eV) -40.03954 -39.71829 -38.03448 -37.64017 -34.19490 -34.02994 -32.82531 -32.31263 -31.58759 -30.38687 -28.46027 -26.78067 -25.51885 -24.63484 -24.06857 -23.58046 -22.05346 -20.83683 -20.47946 -20.19313 -19.52464 -17.56147 -17.41621 -16.93098 -16.59612 -16.34209 -16.08387 -15.94270 -15.31958 -15.15158 -14.59316 -14.44239 -14.14959 -14.04928 -13.70547 -13.54700 -13.24763 -13.14904 -13.02020 -12.84245 -12.59695 -12.51608 -12.21953 -11.86885 -11.78748 -11.77504 -11.58599 -11.48777 -11.15389 -10.97802 -10.85071 -10.78007 -10.67318 -10.62321 -10.41257 -10.12785 -10.07606 -9.98031 -9.79010 -9.66003 -9.58958 -8.65984 -8.29829 -6.94673 -6.12731 -3.06112 2.86138 3.06499 3.40720 3.53587 3.56736 3.61647 3.73766 4.28761 4.33465 4.50389 4.60079 4.62773 4.69632 4.73034 4.77336 4.83189 4.91987 4.94943 5.02766 5.15916 5.17820 5.19309 5.20956 5.22840 5.27037 5.29267 5.34179 5.38554 5.48403 5.56248 5.65161 5.69523 5.78329 5.79641 5.86785 5.97214 6.00098 6.05648 6.17394 6.22446 6.29728 6.58286 6.62950 6.71602 6.80962 6.94368 7.01011 7.10707 7.74426 7.82034 8.46428 8.52045 9.02476 9.84409 10.38736 11.21728 Molecular weight = 340.22amu Principal moments of inertia in cm(-1) A = 0.014374 B = 0.001887 C = 0.001730 Principal moments of inertia in units of 10**(-40)*gram-cm**2 A = 1947.513069 B =14832.753378 C =16182.637871 Net atomic charges, atomic populations, and dipole contributions calculated with CM2 Atom NO. Type Charge No. of electrons 1 C -0.137 4.137 2 C -0.110 4.110 3 H 0.075 0.925 4 C 0.142 3.858 5 N -0.601 5.601 6 C -0.345 4.345 7 H 0.068 0.932 8 C -0.021 4.021 9 N -0.907 5.907 10 Si 1.061 2.939 11 H -0.291 1.291 12 H -0.295 1.295 13 H -0.289 1.289 14 C 0.141 3.859 15 C 0.130 3.870 16 H 0.093 0.907 17 N -0.718 5.718 18 C 0.510 3.490 19 O -0.533 6.533 20 C -0.136 4.136 21 C -0.104 4.104 22 C 0.063 3.937 23 O -0.376 6.376 24 C 0.090 3.910 25 C -0.181 4.181 26 C 0.064 3.936 27 O -0.382 6.382 28 C 0.064 3.936 29 C -0.155 4.155 30 H 0.051 0.949 31 H 0.054 0.946 32 H 0.057 0.943 33 H 0.024 0.976 34 H 0.030 0.970 35 H 0.267 0.733 36 H 0.056 0.944 37 H 0.024 0.976 38 H 0.396 0.604 39 H 0.094 0.906 40 H 0.091 0.909 41 H 0.078 0.922 42 H 0.081 0.919 43 H 0.049 0.951 44 H 0.048 0.952 45 H 0.060 0.940 46 H 0.083 0.917 47 H 0.076 0.924 48 H 0.050 0.950 49 H 0.098 0.902 50 H 0.051 0.949 51 H 0.100 0.900 52 H 0.080 0.920 53 H 0.082 0.918 For the charges calculated by CM2: Dipole moment (debyes) X Y Z Total from point charges -10.133 13.207 20.953 26.761 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.195 4.195 2 C -0.129 4.129 3 H 0.093 0.907 4 C 0.016 3.984 5 N -0.328 5.328 6 C -0.472 4.472 7 H 0.086 0.914 8 C -0.224 4.224 9 N -0.547 5.547 10 Si 0.890 3.110 11 H -0.218 1.218 12 H -0.223 1.223 13 H -0.215 1.215 14 C 0.014 3.986 15 C 0.025 3.975 16 H 0.111 0.889 17 N -0.369 5.369 18 C 0.296 3.704 19 O -0.409 6.409 20 C -0.177 4.177 21 C -0.142 4.142 22 C -0.013 4.013 23 O -0.295 6.295 24 C 0.032 3.968 25 C -0.219 4.219 26 C -0.013 4.013 27 O -0.301 6.301 28 C -0.013 4.013 29 C -0.193 4.193 30 H 0.070 0.930 31 H 0.073 0.927 32 H 0.076 0.924 33 H 0.042 0.958 34 H 0.048 0.952 35 H 0.078 0.922 36 H 0.075 0.925 37 H 0.041 0.959 38 H 0.229 0.771 39 H 0.112 0.888 40 H 0.110 0.890 41 H 0.097 0.903 42 H 0.100 0.900 43 H 0.067 0.933 44 H 0.066 0.934 45 H 0.078 0.922 46 H 0.101 0.899 47 H 0.095 0.905 48 H 0.069 0.931 49 H 0.116 0.884 50 H 0.069 0.931 51 H 0.118 0.882 52 H 0.099 0.901 53 H 0.101 0.899 Dipole moment (debyes) X Y Z Total from point charges -8.835 10.845 21.388 25.556 hybrid contribution -1.198 2.150 -0.837 2.600 sum -10.033 12.995 20.551 26.304 Atomic orbital electron populations 1.21734 0.94646 1.02042 1.01033 1.22419 0.95512 0.99164 0.95851 0.90694 1.21357 0.90558 0.92960 0.93555 1.44473 1.31969 1.51289 1.05107 1.19702 1.02080 1.22956 1.02486 0.91379 1.25678 0.91945 1.05143 0.99597 1.70736 1.35361 1.07812 1.40745 0.83635 0.75060 0.76340 0.75955 1.21754 1.22278 1.21541 1.21769 0.92139 0.98902 0.85752 1.22367 0.99264 0.92864 0.83029 0.88884 1.45929 1.74730 1.09989 1.06249 1.20496 0.76336 0.89219 0.84339 1.90746 1.48968 1.15162 1.86041 1.21529 1.04200 0.98843 0.93080 1.21212 1.05069 0.97820 0.90083 1.22570 0.99807 0.93664 0.85288 1.87965 1.94121 1.33758 1.13685 1.22064 0.94746 0.94685 0.85352 1.22752 0.97746 1.03034 0.98322 1.22566 0.98484 0.93503 0.86761 1.87953 1.21285 1.74914 1.45951 1.22584 0.82271 0.99225 0.97200 1.22311 1.02057 1.01119 0.93778 0.93017 0.92661 0.92376 0.95790 0.95212 0.92168 0.92540 0.95869 0.77070 0.88781 0.89042 0.90285 0.90021 0.93286 0.93392 0.92189 0.89850 0.90549 0.93129 0.88416 0.93068 0.88223 0.90106 0.89886 Number of geometries 1 Number of calculations of the screened coulomb radii 1 The total number of SCF iterations 18. 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.14 -1.78 9.16 37.16 0.34 -1.44 16 2 C -0.11 -1.78 3.16 -88.78 -0.28 -2.06 16 3 H 0.07 1.10 8.00 -51.93 -0.42 0.68 16 4 C 0.14 2.95 6.85 -3.06 -0.02 2.92 16 5 N -0.60 -15.63 3.33 -170.06 -0.57 -16.20 16 6 C -0.35 -10.05 9.97 -15.06 -0.15 -10.20 16 7 H 0.07 1.97 8.02 -52.49 -0.42 1.55 16 8 C -0.02 -0.63 4.96 -17.43 -0.09 -0.72 16 9 N -0.91 -28.19 12.86 55.49 0.71 -27.47 16 10 Si 1.06 26.43 30.03 -169.99 -5.10 21.32 16 11 H -0.29 -6.89 7.11 56.52 0.40 -6.49 16 12 H -0.30 -7.21 7.11 56.52 0.40 -6.81 16 13 H -0.29 -7.23 7.11 56.52 0.40 -6.83 16 14 C 0.14 3.42 5.03 -2.53 -0.01 3.41 16 15 C 0.13 2.41 2.77 -65.78 -0.18 2.23 16 16 H 0.09 1.90 7.63 -51.93 -0.40 1.51 16 17 N -0.72 -10.50 5.23 -51.73 -0.27 -10.77 16 18 C 0.51 7.87 7.83 -10.99 -0.09 7.78 16 19 O -0.53 -10.89 15.73 5.56 0.09 -10.80 16 20 C -0.14 -1.39 5.70 -27.88 -0.16 -1.55 16 21 C -0.10 -1.19 4.87 -26.73 -0.13 -1.32 16 22 C 0.06 0.58 5.75 37.15 0.21 0.79 16 23 O -0.38 -4.10 10.00 -35.23 -0.35 -4.45 16 24 C 0.09 0.74 2.70 -26.55 -0.07 0.67 16 25 C -0.18 -1.40 5.76 -26.59 -0.15 -1.55 16 26 C 0.06 0.43 6.84 37.21 0.25 0.69 16 27 O -0.38 -3.30 10.77 -35.23 -0.38 -3.68 16 28 C 0.06 0.41 6.84 37.20 0.25 0.67 16 29 C -0.15 -1.16 5.79 -26.59 -0.15 -1.31 16 30 H 0.05 0.59 8.14 -51.93 -0.42 0.17 16 31 H 0.05 0.63 8.14 -51.93 -0.42 0.21 16 32 H 0.06 0.80 8.07 -51.93 -0.42 0.38 16 33 H 0.02 0.50 8.14 -51.92 -0.42 0.08 16 34 H 0.03 0.61 7.77 -51.93 -0.40 0.21 16 35 H 0.27 8.32 6.43 -40.82 -0.26 8.06 16 36 H 0.06 1.60 7.81 -51.93 -0.41 1.20 16 37 H 0.02 0.59 7.02 -51.93 -0.36 0.23 16 38 H 0.40 4.37 8.52 -40.82 -0.35 4.02 16 39 H 0.09 0.81 8.14 -51.93 -0.42 0.39 16 40 H 0.09 0.67 8.14 -51.93 -0.42 0.25 16 41 H 0.08 1.00 8.10 -51.92 -0.42 0.58 16 42 H 0.08 1.13 8.10 -51.93 -0.42 0.71 16 43 H 0.05 0.39 7.99 -51.93 -0.42 -0.02 16 44 H 0.05 0.36 8.09 -51.93 -0.42 -0.06 16 45 H 0.06 0.41 7.79 -51.93 -0.40 0.01 16 46 H 0.08 0.77 8.14 -51.93 -0.42 0.34 16 47 H 0.08 0.53 7.99 -51.93 -0.42 0.12 16 48 H 0.05 0.28 8.14 -51.93 -0.42 -0.15 16 49 H 0.10 0.56 8.14 -51.93 -0.42 0.13 16 50 H 0.05 0.26 8.14 -51.92 -0.42 -0.16 16 51 H 0.10 0.51 8.14 -51.92 -0.42 0.08 16 52 H 0.08 0.54 8.14 -51.93 -0.42 0.12 16 53 H 0.08 0.73 8.14 -51.93 -0.42 0.30 16 LS Contribution 418.27 15.07 6.30 6.30 Total: -1.00 -36.14 418.27 -9.79 -45.94 By element: Atomic # 1 Polarization: 10.61 SS G_CDS: -9.80 Total: 0.81 kcal Atomic # 6 Polarization: -0.57 SS G_CDS: -0.42 Total: -1.00 kcal Atomic # 7 Polarization: -54.32 SS G_CDS: -0.12 Total: -54.45 kcal Atomic # 8 Polarization: -18.29 SS G_CDS: -0.64 Total: -18.93 kcal Atomic # 14 Polarization: 26.43 SS G_CDS: -5.10 Total: 21.32 kcal Total LS contribution 6.30 Total: 6.30 kcal Total: -36.14 -9.79 -45.94 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. ZINC001210193454.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 -80.245 kcal (2) G-P(sol) polarization free energy of solvation -36.144 kcal (3) G-ENP(sol) elect.-nuc.-pol. free energy of system -116.389 kcal (4) G-CDS(sol) cavity-dispersion-solvent structure free energy -9.792 kcal (5) G-P-CDS(sol) = G-P(sol) + G-CDS(sol) = (2) + (4) -45.935 kcal (6) G-S(sol) free energy of system = (1) + (5) -126.181 kcal 1SCF WAS SPECIFIED, SO GEOMETRY OPTIMIZATION WAS NOT USED Total computer time = 0.89 seconds