Wall clock time and date at job start Sat Apr 11 2020 06:10:30 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.52994 * 1 3 3 C 1.52995 * 109.47489 * 2 1 4 4 C 1.53000 * 109.46833 * 240.00348 * 2 1 3 5 5 C 1.50699 * 109.47173 * 60.00124 * 4 2 1 6 6 H 1.07993 * 120.00080 * 359.97438 * 5 4 2 7 7 C 1.37884 * 119.99943 * 180.02562 * 5 4 2 8 8 N 1.31509 * 119.99957 * 359.97438 * 7 5 4 9 9 Si 1.86800 * 119.99777 * 179.97438 * 7 5 4 10 10 H 1.48493 * 109.47158 * 60.00085 * 9 7 5 11 11 H 1.48500 * 109.47137 * 180.02562 * 9 7 5 12 12 H 1.48501 * 109.47195 * 299.99630 * 9 7 5 13 13 C 1.52994 * 109.47387 * 120.00537 * 2 1 3 14 14 N 1.46501 * 109.47184 * 299.99745 * 13 2 1 15 15 C 1.34777 * 119.99898 * 180.02562 * 14 13 2 16 16 O 1.21581 * 120.00353 * 355.01471 * 15 14 13 17 17 N 1.34780 * 119.99586 * 175.01694 * 15 14 13 18 18 C 1.46924 * 120.63194 * 6.28077 * 17 15 14 19 19 C 1.53619 * 108.54736 * 126.36537 * 18 17 15 20 20 C 1.53047 * 109.24142 * 54.85276 * 19 18 17 21 21 C 1.53038 * 109.53906 * 298.51141 * 20 19 18 22 22 C 1.46921 * 120.62976 * 186.55278 * 17 15 14 23 23 H 1.08997 * 109.58915 * 113.80008 * 22 17 15 24 24 C 1.53001 * 109.58871 * 353.51915 * 22 17 15 25 25 O 1.43755 * 109.63703 * 186.70781 * 24 22 17 26 26 C 1.44581 * 104.91971 * 215.70914 * 25 24 22 27 27 C 1.54916 * 101.63907 * 36.26773 * 26 25 24 28 28 O 1.43371 * 109.63089 * 67.49052 * 24 22 17 29 29 H 1.09006 * 109.46868 * 299.99793 * 1 2 3 30 30 H 1.09000 * 109.47276 * 59.99381 * 1 2 3 31 31 H 1.08998 * 109.47674 * 179.97438 * 1 2 3 32 32 H 1.08999 * 109.47225 * 180.02562 * 3 2 1 33 33 H 1.09000 * 109.47367 * 300.00564 * 3 2 1 34 34 H 1.09006 * 109.46897 * 60.00224 * 3 2 1 35 35 H 1.08999 * 109.47520 * 179.97438 * 4 2 1 36 36 H 1.09004 * 109.47611 * 299.99956 * 4 2 1 37 37 H 0.97004 * 119.99759 * 179.97438 * 8 7 5 38 38 H 1.08999 * 109.47412 * 179.97438 * 13 2 1 39 39 H 1.09004 * 109.47386 * 59.99834 * 13 2 1 40 40 H 0.97000 * 119.99695 * 359.97438 * 14 13 2 41 41 H 1.09000 * 109.58194 * 246.20188 * 18 17 15 42 42 H 1.08998 * 109.58547 * 6.62138 * 18 17 15 43 43 H 1.09004 * 109.63613 * 294.86432 * 19 18 17 44 44 H 1.09003 * 109.45946 * 174.73893 * 19 18 17 45 45 H 1.08998 * 109.46110 * 178.48971 * 20 19 18 46 46 H 1.08999 * 109.44731 * 58.52689 * 20 19 18 47 47 H 1.08996 * 109.49993 * 301.41770 * 21 20 19 48 48 H 1.08994 * 109.52218 * 181.32871 * 21 20 19 49 49 H 1.09005 * 109.62706 * 307.10073 * 24 22 17 50 50 H 1.08993 * 111.08914 * 278.15599 * 26 25 24 51 51 H 1.08996 * 110.97452 * 154.32905 * 26 25 24 52 52 H 1.08989 * 110.69631 * 205.95905 * 27 26 25 53 53 H 1.09002 * 110.69173 * 82.75142 * 27 26 25 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.5299 0.0000 0.0000 3 6 2.0400 1.4424 0.0000 4 6 2.0399 -0.7212 -1.2493 5 6 1.5375 -0.0107 -2.4797 6 1 0.9164 0.8671 -2.3806 7 6 1.8702 -0.4820 -3.7321 8 7 2.6262 -1.5513 -3.8528 9 14 1.2469 0.3983 -5.2572 10 1 1.7419 1.7983 -5.2576 11 1 1.7409 -0.3023 -6.4698 12 1 -0.2381 0.3982 -5.2566 13 6 2.0400 -0.7213 1.2491 14 7 1.5517 -0.0308 2.4454 15 6 1.8779 -0.4911 3.6694 16 8 2.5004 -1.5293 3.7826 17 7 1.5137 0.1995 4.7681 18 6 0.6456 1.3796 4.6569 19 6 -0.5573 1.1914 5.5935 20 6 -0.0537 0.9134 7.0118 21 6 0.7687 -0.3773 7.0191 22 6 1.9821 -0.2100 6.0990 23 1 2.6452 0.5544 6.5042 24 6 2.7350 -1.5376 5.9918 25 8 3.0571 -2.0179 7.3079 26 6 4.3506 -2.6486 7.1678 27 6 5.0392 -1.7231 6.1338 28 8 3.9544 -1.3424 5.2634 29 1 -0.3633 0.5138 0.8901 30 1 -0.3634 0.5139 -0.8899 31 1 -0.3634 -1.0276 -0.0005 32 1 3.1300 1.4424 -0.0005 33 1 1.6767 1.9563 -0.8899 34 1 1.6767 1.9563 0.8901 35 1 3.1299 -0.7209 -1.2496 36 1 1.6765 -1.7489 -1.2494 37 1 2.8598 -1.8831 -4.7339 38 1 3.1300 -0.7217 1.2489 39 1 1.6767 -1.7490 1.2491 40 1 0.9944 0.7581 2.3564 41 1 1.2020 2.2710 4.9467 42 1 0.2977 1.4818 3.6290 43 1 -1.1606 0.3508 5.2505 44 1 -1.1629 2.0977 5.5936 45 1 -0.9043 0.8057 7.6847 46 1 0.5699 1.7431 7.3446 47 1 0.1530 -1.2030 6.6625 48 1 1.1075 -0.5876 8.0334 49 1 2.1139 -2.2702 5.4765 50 1 4.2527 -3.6624 6.7798 51 1 4.8906 -2.6457 8.1145 52 1 5.8047 -2.2665 5.5801 53 1 5.4651 -0.8476 6.6240 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 ZINC000807593713.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:10:30 Heat of formation + Delta-G solvation = -109.497400 kcal Electronic energy + Delta-G solvation = -32151.702436 eV Core-core repulsion = 28023.992013 eV Total energy + Delta-G solvation = -4127.710422 eV No. of doubly occupied orbitals = 66 Molecular weight (most abundant/longest-lived isotopes) = 340.220 amu Computer time = 0.71 seconds Orbital eigenvalues (eV) -41.02058 -39.66173 -37.88168 -36.81611 -34.83427 -33.92371 -32.80016 -32.34338 -31.49814 -29.40883 -27.14948 -26.49594 -26.47891 -25.41341 -24.42937 -23.20497 -22.90535 -21.91878 -21.01814 -19.70173 -18.58409 -17.59123 -17.14956 -16.47984 -16.33668 -16.17049 -15.61808 -15.41027 -14.86516 -14.58923 -14.51868 -14.42639 -14.12448 -13.87003 -13.61820 -13.44350 -13.21146 -13.04923 -12.93932 -12.83618 -12.59499 -12.38040 -12.04700 -12.03151 -11.66451 -11.61749 -11.22573 -11.08980 -11.04244 -10.97753 -10.91041 -10.84466 -10.60149 -10.36648 -10.32495 -10.23121 -9.84480 -9.74786 -9.22753 -9.05885 -8.90137 -8.67856 -8.37073 -8.33391 -5.97631 -3.98840 2.62827 3.33175 3.38581 3.39705 3.40018 3.56525 3.92978 4.19138 4.46722 4.53750 4.54550 4.58863 4.67784 4.80464 4.86368 5.01726 5.04070 5.15117 5.21181 5.33299 5.36589 5.37772 5.49676 5.55157 5.56443 5.63723 5.64254 5.69505 5.74498 5.79315 5.79692 5.84459 5.95181 6.13807 6.21392 6.25457 6.32532 6.46862 6.50011 6.65936 6.66470 6.88258 7.03311 7.10680 7.11077 7.19391 7.25383 7.67275 7.80370 8.01215 8.18645 8.43338 8.88530 9.05836 9.63968 11.13822 Molecular weight = 340.22amu Principal moments of inertia in cm(-1) A = 0.016911 B = 0.002735 C = 0.002437 Principal moments of inertia in units of 10**(-40)*gram-cm**2 A = 1655.282086 B =10234.542452 C =11487.630012 Net atomic charges, atomic populations, and dipole contributions calculated with CM2 Atom NO. Type Charge No. of electrons 1 C -0.125 4.125 2 C -0.073 4.073 3 C -0.125 4.125 4 C 0.037 3.963 5 C -0.511 4.511 6 H 0.058 0.942 7 C 0.056 3.944 8 N -0.895 5.895 9 Si 0.889 3.111 10 H -0.278 1.278 11 H -0.288 1.288 12 H -0.278 1.278 13 C 0.148 3.852 14 N -0.737 5.737 15 C 0.700 3.300 16 O -0.559 6.559 17 N -0.629 5.629 18 C 0.114 3.886 19 C -0.136 4.136 20 C -0.117 4.117 21 C -0.124 4.124 22 C 0.155 3.845 23 H 0.073 0.927 24 C 0.246 3.754 25 O -0.393 6.393 26 C 0.023 3.977 27 C 0.006 3.994 28 O -0.336 6.336 29 H 0.026 0.974 30 H 0.066 0.934 31 H 0.047 0.953 32 H 0.047 0.953 33 H 0.066 0.934 34 H 0.026 0.974 35 H 0.018 0.982 36 H 0.018 0.982 37 H 0.260 0.740 38 H 0.061 0.939 39 H 0.061 0.939 40 H 0.391 0.609 41 H 0.068 0.932 42 H 0.080 0.920 43 H 0.069 0.931 44 H 0.075 0.925 45 H 0.067 0.933 46 H 0.063 0.937 47 H 0.068 0.932 48 H 0.081 0.919 49 H 0.107 0.893 50 H 0.074 0.926 51 H 0.101 0.899 52 H 0.094 0.906 53 H 0.065 0.935 For the charges calculated by CM2: Dipole moment (debyes) X Y Z Total from point charges -2.421 4.696 25.031 25.582 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.182 4.182 2 C -0.073 4.073 3 C -0.182 4.182 4 C 0.000 4.000 5 C -0.550 4.550 6 H 0.077 0.923 7 C -0.145 4.145 8 N -0.534 5.534 9 Si 0.717 3.283 10 H -0.203 1.203 11 H -0.214 1.214 12 H -0.204 1.204 13 C 0.025 3.975 14 N -0.395 5.395 15 C 0.402 3.598 16 O -0.437 6.437 17 N -0.369 5.369 18 C -0.009 4.009 19 C -0.175 4.175 20 C -0.155 4.155 21 C -0.163 4.163 22 C 0.050 3.950 23 H 0.091 0.909 24 C 0.152 3.848 25 O -0.313 6.313 26 C -0.054 4.054 27 C -0.071 4.071 28 O -0.256 6.256 29 H 0.045 0.955 30 H 0.085 0.915 31 H 0.066 0.934 32 H 0.066 0.934 33 H 0.085 0.915 34 H 0.045 0.955 35 H 0.037 0.963 36 H 0.037 0.963 37 H 0.070 0.930 38 H 0.079 0.921 39 H 0.079 0.921 40 H 0.228 0.772 41 H 0.086 0.914 42 H 0.098 0.902 43 H 0.088 0.912 44 H 0.093 0.907 45 H 0.085 0.915 46 H 0.082 0.918 47 H 0.087 0.913 48 H 0.100 0.900 49 H 0.124 0.876 50 H 0.092 0.908 51 H 0.119 0.881 52 H 0.112 0.888 53 H 0.083 0.917 Dipole moment (debyes) X Y Z Total from point charges -1.970 3.898 24.906 25.286 hybrid contribution -0.407 0.822 -0.963 1.330 sum -2.377 4.720 23.943 24.519 Atomic orbital electron populations 1.21741 0.95385 1.01106 0.99991 1.20909 0.95383 0.95849 0.95166 1.21740 1.00553 0.95942 0.99988 1.19072 0.96355 0.95198 0.89393 1.21152 1.31187 1.12057 0.90555 0.92347 1.25010 0.95338 0.95083 0.99041 1.69956 1.36928 1.23006 1.23505 0.86718 0.79128 0.79766 0.82719 1.20347 1.21415 1.20354 1.21453 0.98571 0.95616 0.81875 1.45565 1.56146 1.33228 1.04530 1.16005 0.79780 0.83575 0.80453 1.90938 1.45384 1.23394 1.83938 1.47192 1.53610 1.31074 1.05058 1.21691 0.91058 0.87943 1.00220 1.22019 0.97412 1.02556 0.95478 1.21495 0.99726 0.97597 0.96647 1.21844 0.96171 0.98370 0.99883 1.21329 0.95157 0.96243 0.82258 0.90868 1.24255 0.88801 0.94652 0.77109 1.89105 1.35951 1.75184 1.31085 1.24239 0.83429 0.97313 1.00412 1.24188 0.92558 0.97846 0.92518 1.89397 1.13000 1.87619 1.35609 0.95488 0.91532 0.93376 0.93374 0.91547 0.95465 0.96341 0.96345 0.93003 0.92095 0.92098 0.77176 0.91422 0.90183 0.91234 0.90662 0.91476 0.91777 0.91284 0.90043 0.87629 0.90826 0.88078 0.88818 0.91727 Number of geometries 1 Number of calculations of the screened coulomb radii 1 The total number of SCF iterations 12. 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.12 -2.22 8.05 37.15 0.30 -1.92 16 2 C -0.07 -1.40 0.54 -154.52 -0.08 -1.48 16 3 C -0.12 -2.23 8.05 37.16 0.30 -1.93 16 4 C 0.04 0.92 4.04 -27.88 -0.11 0.81 16 5 C -0.51 -13.63 7.65 -34.44 -0.26 -13.89 16 6 H 0.06 1.49 5.68 -52.49 -0.30 1.19 16 7 C 0.06 1.57 5.46 -17.82 -0.10 1.48 16 8 N -0.89 -26.56 13.29 55.43 0.74 -25.83 16 9 Si 0.89 21.27 29.54 -169.99 -5.02 16.25 16 10 H -0.28 -6.53 7.11 56.52 0.40 -6.13 16 11 H -0.29 -6.92 7.11 56.52 0.40 -6.52 16 12 H -0.28 -6.53 7.11 56.52 0.40 -6.13 16 13 C 0.15 2.72 4.57 -4.04 -0.02 2.70 16 14 N -0.74 -11.37 4.87 -65.23 -0.32 -11.69 16 15 C 0.70 11.25 7.99 -86.91 -0.69 10.55 16 16 O -0.56 -10.68 10.51 -20.05 -0.21 -10.89 16 17 N -0.63 -7.90 2.97 -170.65 -0.51 -8.41 16 18 C 0.11 1.06 6.40 -3.49 -0.02 1.03 16 19 C -0.14 -1.03 6.10 -26.39 -0.16 -1.20 16 20 C -0.12 -0.87 5.88 -26.69 -0.16 -1.02 16 21 C -0.12 -1.17 5.07 -26.61 -0.13 -1.31 16 22 C 0.15 1.84 2.96 -67.60 -0.20 1.64 16 23 H 0.07 0.86 8.14 -51.93 -0.42 0.44 16 24 C 0.25 3.43 2.68 37.16 0.10 3.53 16 25 O -0.39 -4.80 10.56 -59.05 -0.62 -5.43 16 26 C 0.02 0.22 7.95 38.17 0.30 0.52 16 27 C 0.01 0.07 8.07 38.17 0.31 0.38 16 28 O -0.34 -5.29 8.31 -68.04 -0.57 -5.85 16 29 H 0.03 0.39 6.86 -51.93 -0.36 0.03 16 30 H 0.07 1.26 6.62 -51.93 -0.34 0.92 16 31 H 0.05 0.88 8.14 -51.93 -0.42 0.45 16 32 H 0.05 0.89 8.14 -51.93 -0.42 0.46 16 33 H 0.07 1.26 6.62 -51.93 -0.34 0.92 16 34 H 0.03 0.40 6.86 -51.93 -0.36 0.04 16 35 H 0.02 0.49 8.14 -51.93 -0.42 0.07 16 36 H 0.02 0.49 8.14 -51.93 -0.42 0.07 16 37 H 0.26 7.36 8.31 -40.82 -0.34 7.03 16 38 H 0.06 1.21 8.14 -51.93 -0.42 0.79 16 39 H 0.06 1.20 8.14 -51.93 -0.42 0.77 16 40 H 0.39 5.20 3.67 -40.82 -0.15 5.05 16 41 H 0.07 0.56 8.14 -51.93 -0.42 0.13 16 42 H 0.08 0.72 5.67 -51.93 -0.29 0.42 16 43 H 0.07 0.59 8.14 -51.93 -0.42 0.17 16 44 H 0.07 0.44 8.14 -51.93 -0.42 0.02 16 45 H 0.07 0.43 8.14 -51.93 -0.42 0.01 16 46 H 0.06 0.44 8.14 -51.93 -0.42 0.02 16 47 H 0.07 0.72 8.14 -51.93 -0.42 0.29 16 48 H 0.08 0.73 7.75 -51.93 -0.40 0.33 16 49 H 0.11 1.73 6.46 -51.93 -0.34 1.40 16 50 H 0.07 0.70 8.14 -51.93 -0.42 0.27 16 51 H 0.10 0.70 8.14 -51.93 -0.42 0.27 16 52 H 0.09 0.94 8.14 -51.94 -0.42 0.52 16 53 H 0.06 0.64 8.14 -51.93 -0.42 0.21 16 LS Contribution 395.72 15.07 5.96 5.96 Total: -1.00 -32.06 395.72 -10.38 -42.44 By element: Atomic # 1 Polarization: 12.74 SS G_CDS: -9.20 Total: 3.54 kcal Atomic # 6 Polarization: 0.54 SS G_CDS: -0.64 Total: -0.10 kcal Atomic # 7 Polarization: -45.84 SS G_CDS: -0.09 Total: -45.92 kcal Atomic # 8 Polarization: -20.77 SS G_CDS: -1.40 Total: -22.17 kcal Atomic # 14 Polarization: 21.27 SS G_CDS: -5.02 Total: 16.25 kcal Total LS contribution 5.96 Total: 5.96 kcal Total: -32.06 -10.38 -42.44 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. ZINC000807593713.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 -67.054 kcal (2) G-P(sol) polarization free energy of solvation -32.060 kcal (3) G-ENP(sol) elect.-nuc.-pol. free energy of system -99.115 kcal (4) G-CDS(sol) cavity-dispersion-solvent structure free energy -10.383 kcal (5) G-P-CDS(sol) = G-P(sol) + G-CDS(sol) = (2) + (4) -42.443 kcal (6) G-S(sol) free energy of system = (1) + (5) -109.497 kcal 1SCF WAS SPECIFIED, SO GEOMETRY OPTIMIZATION WAS NOT USED Total computer time = 0.71 seconds