Wall clock time and date at job start Tue Mar 31 2020 06:53:13 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.50704 * 1 3 3 N 1.29946 * 126.28738 * 2 1 4 4 O 1.21302 * 110.25120 * 179.97438 * 3 2 1 5 5 C 1.33643 * 109.81254 * 0.02562 * 4 3 2 6 6 C 1.50706 * 126.55398 * 179.97438 * 5 4 3 7 7 H 1.08994 * 109.46953 * 195.02426 * 6 5 4 8 8 C 1.52996 * 109.46904 * 75.02632 * 6 5 4 9 9 N 1.46901 * 109.46898 * 315.02210 * 6 5 4 10 10 C 1.46950 * 110.99804 * 179.97438 * 9 6 5 11 11 C 1.53053 * 109.36589 * 183.19091 * 10 9 6 12 12 N 1.46951 * 109.36240 * 301.60134 * 11 10 9 13 13 C 1.46901 * 111.00160 * 183.18681 * 12 11 10 14 14 C 1.52993 * 109.47476 * 170.00059 * 13 12 11 15 15 C 1.50700 * 109.47277 * 180.02562 * 14 13 12 16 16 O 1.21287 * 120.00075 * 359.97438 * 15 14 13 17 17 N 1.34783 * 120.00195 * 179.72358 * 15 14 13 18 18 C 1.37100 * 120.00185 * 180.27153 * 17 15 14 19 19 H 1.08001 * 119.99500 * 17.72019 * 18 17 15 20 20 C 1.38943 * 119.99647 * 197.72103 * 18 17 15 21 21 N 1.31408 * 120.00658 * 21.38805 * 20 18 17 22 22 Si 1.86807 * 119.99682 * 201.39062 * 20 18 17 23 23 H 1.48496 * 109.47238 * 89.99773 * 22 20 18 24 24 H 1.48499 * 109.46998 * 210.00254 * 22 20 18 25 25 H 1.48501 * 109.46883 * 329.99624 * 22 20 18 26 26 C 1.46502 * 119.99501 * 0.27093 * 17 15 14 27 27 C 1.47205 * 109.47195 * 84.23457 * 26 17 15 28 28 N 9.22447 * 113.92768 * 277.75702 * 2 1 3 29 29 C 1.46964 * 110.88374 * 59.31840 * 12 11 10 30 30 C 1.46951 * 111.00160 * 303.80064 * 9 6 5 31 31 N 1.31131 * 106.89448 * 0.22510 * 5 4 3 32 32 H 1.09003 * 109.46361 * 89.97786 * 1 2 3 33 33 H 1.08996 * 109.47031 * 209.97457 * 1 2 3 34 34 H 1.08992 * 109.46874 * 329.97944 * 1 2 3 35 35 H 1.09000 * 109.46892 * 180.02562 * 8 6 5 36 36 H 1.09003 * 109.46977 * 299.99916 * 8 6 5 37 37 H 1.09000 * 109.47725 * 59.99806 * 8 6 5 38 38 H 1.08995 * 109.48851 * 63.22317 * 10 9 6 39 39 H 1.08999 * 109.49056 * 303.15711 * 10 9 6 40 40 H 1.08999 * 109.48846 * 181.63194 * 11 10 9 41 41 H 1.09005 * 109.48544 * 61.56336 * 11 10 9 42 42 H 1.09001 * 109.47334 * 49.99538 * 13 12 11 43 43 H 1.09001 * 109.47069 * 290.00314 * 13 12 11 44 44 H 1.09001 * 109.47702 * 59.99743 * 14 13 12 45 45 H 1.09010 * 109.47058 * 299.99966 * 14 13 12 46 46 H 0.96999 * 120.00569 * 179.97438 * 21 20 18 47 47 H 1.09001 * 109.46693 * 204.22592 * 26 17 15 48 48 H 1.08998 * 109.47521 * 324.23176 * 26 17 15 49 49 H 1.09001 * 109.48344 * 60.64654 * 29 12 11 50 50 H 1.08993 * 109.48919 * 180.70949 * 29 12 11 51 51 H 1.08997 * 109.48951 * 296.77619 * 30 9 6 52 52 H 1.09004 * 109.49195 * 56.82315 * 30 9 6 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 7 2.2761 1.0474 0.0000 4 8 3.4419 0.7124 0.0005 5 6 3.5299 -0.6212 0.0003 6 6 4.7970 -1.4371 0.0003 7 1 4.5652 -2.4682 -0.2663 8 6 5.4275 -1.4000 1.3938 9 7 5.7391 -0.8783 -0.9786 10 6 6.9895 -1.6499 -1.0059 11 6 7.9716 -0.9906 -1.9772 12 7 7.3657 -0.9273 -3.3145 13 6 8.3078 -0.3685 -4.2934 14 6 7.7431 -0.5467 -5.7041 15 6 8.7099 0.0260 -6.7083 16 8 9.7482 0.5278 -6.3326 17 7 8.4231 -0.0238 -8.0243 18 6 9.3052 0.4922 -8.9382 19 1 10.3314 0.6739 -8.6550 20 6 8.8788 0.7812 -10.2286 21 7 7.6011 0.9691 -10.4720 22 14 10.1218 0.9028 -11.6178 23 1 10.2988 -0.4330 -12.2419 24 1 9.6320 1.8647 -12.6375 25 1 11.4228 1.3734 -11.0780 26 6 7.1687 -0.6295 -8.4780 27 6 6.0885 0.3681 -8.4080 28 7 5.2483 1.1380 -8.3545 29 6 6.1153 -0.1557 -3.2872 30 6 5.1331 -0.8145 -2.3158 31 7 2.3031 -1.0841 -0.0053 32 1 -0.3632 0.0004 -1.0277 33 1 -0.3633 -0.8902 0.5134 34 1 -0.3633 0.8897 0.5141 35 1 6.3442 -1.9897 1.3936 36 1 5.6593 -0.3688 1.6603 37 1 4.7285 -1.8146 2.1201 38 1 6.7813 -2.6681 -1.3345 39 1 7.4253 -1.6716 -0.0071 40 1 8.8890 -1.5775 -2.0210 41 1 8.2007 0.0182 -1.6337 42 1 9.2630 -0.8880 -4.2161 43 1 8.4542 0.6927 -4.0920 44 1 6.7880 -0.0271 -5.7815 45 1 7.5968 -1.6079 -5.9055 46 1 7.3035 1.1712 -11.3727 47 1 7.2815 -0.9711 -9.5070 48 1 6.9223 -1.4772 -7.8387 49 1 6.3235 0.8626 -2.9586 50 1 5.6796 -0.1338 -4.2860 51 1 4.9036 -1.8230 -2.6595 52 1 4.2156 -0.2276 -2.2717 RHF calculation, no. of doubly occupied orbitals= 71 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 ZINC000069370378.mol2 52 1SCF WAS SPECIFIED, SO GEOMETRY OPTIMIZATION WAS NOT USED SCF FIELD WAS ACHIEVED Wall clock time and date at job start Tue Mar 31 2020 06:53:13 Heat of formation + Delta-G solvation = 93.860235 kcal Electronic energy + Delta-G solvation = -34753.383842 eV Core-core repulsion = 30178.565655 eV Total energy + Delta-G solvation = -4574.818187 eV No. of doubly occupied orbitals = 71 Molecular weight (most abundant/longest-lived isotopes) = 376.201 amu Computer time = 1.56 seconds Orbital eigenvalues (eV) -45.86258 -39.86607 -38.86045 -36.86853 -36.27384 -35.71975 -34.31452 -33.51975 -32.02221 -31.68996 -31.24821 -29.78062 -29.35568 -27.46445 -26.64577 -25.18379 -24.43083 -23.65238 -23.37844 -22.06558 -20.93582 -19.91582 -18.95572 -18.44258 -17.91205 -17.67470 -16.91352 -16.41969 -16.23698 -16.08424 -15.36728 -15.09783 -14.93948 -14.72881 -14.38694 -14.37260 -14.14874 -14.01136 -13.92877 -13.67659 -13.34548 -12.97140 -12.95458 -12.81035 -12.68951 -12.62015 -12.45165 -12.18097 -12.00586 -11.85879 -11.68382 -11.64071 -11.53636 -11.42391 -11.32584 -11.26573 -11.17457 -11.10304 -10.87032 -10.58465 -10.45301 -10.39766 -10.27697 -10.15245 -9.27273 -8.85077 -8.60653 -8.30291 -8.05646 -6.51282 -4.00665 0.84740 1.89288 2.48059 3.00996 3.09521 3.21656 3.25645 3.52214 3.53944 3.73990 3.79870 3.85668 3.95002 4.12640 4.35134 4.48785 4.58800 4.74146 4.79799 4.87292 4.92954 5.01578 5.03272 5.06807 5.10676 5.14004 5.16713 5.19695 5.22937 5.28837 5.36147 5.45753 5.48575 5.54686 5.55454 5.65777 5.67534 5.68710 5.78911 5.83155 5.88423 5.98975 6.12717 6.28229 6.32886 6.47051 6.50121 6.56262 7.06810 7.09492 7.46463 7.89686 7.90367 8.27312 8.72583 8.99613 9.19231 9.73262 10.54880 Molecular weight = 376.20amu Principal moments of inertia in cm(-1) A = 0.014118 B = 0.001919 C = 0.001735 Principal moments of inertia in units of 10**(-40)*gram-cm**2 A = 1982.776761 B =14589.621740 C =16133.473653 Net atomic charges, atomic populations, and dipole contributions calculated with CM2 Atom NO. Type Charge No. of electrons 1 C -0.062 4.062 2 C 0.164 3.836 3 N -0.386 5.386 4 O 0.095 5.905 5 C 0.195 3.805 6 C 0.175 3.825 7 H 0.090 0.910 8 C -0.148 4.148 9 N -0.534 5.534 10 C 0.058 3.942 11 C 0.058 3.942 12 N -0.541 5.541 13 C 0.078 3.922 14 C -0.132 4.132 15 C 0.457 3.543 16 O -0.576 6.576 17 N -0.444 5.444 18 C -0.313 4.313 19 H 0.109 0.891 20 C 0.031 3.969 21 N -0.878 5.878 22 Si 0.927 3.073 23 H -0.276 1.276 24 H -0.286 1.286 25 H -0.269 1.269 26 C 0.222 3.778 27 C 0.245 3.755 28 N -0.474 5.474 29 C 0.055 3.945 30 C 0.058 3.942 31 N -0.483 5.483 32 H 0.089 0.911 33 H 0.091 0.909 34 H 0.086 0.914 35 H 0.081 0.919 36 H 0.075 0.925 37 H 0.069 0.931 38 H 0.050 0.950 39 H 0.088 0.912 40 H 0.085 0.915 41 H 0.045 0.955 42 H 0.083 0.917 43 H 0.042 0.958 44 H 0.076 0.924 45 H 0.089 0.911 46 H 0.273 0.727 47 H 0.127 0.873 48 H 0.085 0.915 49 H 0.046 0.954 50 H 0.085 0.915 51 H 0.051 0.949 52 H 0.068 0.932 For the charges calculated by CM2: Dipole moment (debyes) X Y Z Total from point charges -8.141 -10.463 21.260 25.055 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.121 4.121 2 C -0.109 4.109 3 N -0.093 5.093 4 O 0.036 5.964 5 C -0.027 4.027 6 C 0.065 3.935 7 H 0.108 0.892 8 C -0.206 4.206 9 N -0.261 5.261 10 C -0.069 4.069 11 C -0.071 4.071 12 N -0.267 5.267 13 C -0.050 4.050 14 C -0.172 4.172 15 C 0.246 3.754 16 O -0.459 6.459 17 N -0.164 5.164 18 C -0.440 4.440 19 H 0.126 0.874 20 C -0.173 4.173 21 N -0.517 5.517 22 Si 0.754 3.246 23 H -0.201 1.201 24 H -0.212 1.212 25 H -0.194 1.194 26 C 0.095 3.905 27 C -0.073 4.073 28 N -0.155 5.155 29 C -0.074 4.074 30 C -0.069 4.069 31 N -0.218 5.218 32 H 0.108 0.892 33 H 0.109 0.891 34 H 0.105 0.895 35 H 0.100 0.900 36 H 0.094 0.906 37 H 0.088 0.912 38 H 0.068 0.932 39 H 0.106 0.894 40 H 0.103 0.897 41 H 0.063 0.937 42 H 0.102 0.898 43 H 0.061 0.939 44 H 0.094 0.906 45 H 0.108 0.892 46 H 0.084 0.916 47 H 0.145 0.855 48 H 0.103 0.897 49 H 0.064 0.936 50 H 0.104 0.896 51 H 0.069 0.931 52 H 0.086 0.914 Dipole moment (debyes) X Y Z Total from point charges -9.610 -8.504 21.335 24.897 hybrid contribution 0.685 -0.012 -1.120 1.313 sum -8.925 -8.517 20.215 23.682 Atomic orbital electron populations 1.20139 0.82888 1.04586 1.04439 1.25913 1.00626 0.83469 1.00847 1.74692 0.89584 1.23990 1.21038 1.85367 1.21831 1.24427 1.64742 1.28542 0.96239 0.82400 0.95470 1.20100 0.83837 0.93773 0.95817 0.89230 1.21959 1.01508 1.04341 0.92840 1.62231 1.04581 1.52668 1.06583 1.22283 0.89038 0.95518 1.00042 1.22125 0.96644 1.00318 0.87983 1.62183 1.05232 1.52801 1.06440 1.21741 0.94794 0.98600 0.89820 1.21574 0.99160 1.02441 0.94073 1.20105 0.88632 0.83944 0.82754 1.90540 1.27542 1.49067 1.78760 1.43657 1.16157 1.51870 1.04752 1.19725 0.93112 1.41305 0.89880 0.87383 1.25278 0.96704 0.96014 0.99269 1.69855 1.14420 1.45039 1.22342 0.85899 0.80526 0.78095 0.80089 1.20139 1.21161 1.19356 1.18740 0.79154 0.88882 1.03771 1.29550 0.92905 0.94232 0.90593 1.81107 1.13375 1.09289 1.11696 1.22170 0.90475 0.95565 0.99150 1.22202 0.96889 1.00465 0.87384 1.71888 0.99522 1.29626 1.20759 0.89234 0.89061 0.89494 0.90010 0.90586 0.91246 0.93181 0.89397 0.89688 0.93657 0.89832 0.93926 0.90551 0.89231 0.91640 0.85546 0.89718 0.93616 0.89633 0.93102 0.91384 Number of geometries 1 Number of calculations of the screened coulomb radii 1 The total number of SCF iterations 17. 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.43 10.51 36.01 0.38 -0.05 16 2 C 0.16 1.79 8.22 -156.37 -1.28 0.51 16 3 N -0.39 -4.86 12.16 30.90 0.38 -4.48 16 4 O 0.10 1.16 11.12 18.10 0.20 1.36 16 5 C 0.19 2.02 5.65 -13.05 -0.07 1.95 16 6 C 0.18 1.38 2.62 -68.86 -0.18 1.20 16 7 H 0.09 0.65 8.14 -51.93 -0.42 0.23 16 8 C -0.15 -0.83 8.95 37.16 0.33 -0.50 16 9 N -0.53 -4.86 4.57 -227.89 -1.04 -5.91 16 10 C 0.06 0.46 5.32 -3.77 -0.02 0.44 16 11 C 0.06 0.57 5.61 -3.77 -0.02 0.55 16 12 N -0.54 -6.65 4.71 -234.47 -1.10 -7.76 16 13 C 0.08 1.17 4.25 -3.83 -0.02 1.16 16 14 C -0.13 -2.17 4.24 -27.89 -0.12 -2.29 16 15 C 0.46 10.01 7.57 -10.98 -0.08 9.93 16 16 O -0.58 -14.60 15.63 5.55 0.09 -14.52 16 17 N -0.44 -10.14 2.60 -117.93 -0.31 -10.45 16 18 C -0.31 -7.96 9.32 -14.93 -0.14 -8.10 16 19 H 0.11 2.91 7.39 -52.49 -0.39 2.52 16 20 C 0.03 0.76 4.87 -17.47 -0.09 0.67 16 21 N -0.88 -22.03 9.48 55.50 0.53 -21.50 16 22 Si 0.93 19.53 29.60 -169.99 -5.03 14.50 16 23 H -0.28 -5.61 7.11 56.52 0.40 -5.21 16 24 H -0.29 -6.12 7.11 56.52 0.40 -5.72 16 25 H -0.27 -5.75 7.11 56.52 0.40 -5.35 16 26 C 0.22 4.31 4.91 -6.83 -0.03 4.28 16 27 C 0.25 5.35 12.98 -20.60 -0.27 5.08 16 28 N -0.47 -11.20 18.57 42.35 0.79 -10.42 16 29 C 0.05 0.68 5.41 -3.77 -0.02 0.66 16 30 C 0.06 0.63 5.34 -3.77 -0.02 0.61 16 31 N -0.48 -5.41 11.55 -9.86 -0.11 -5.53 16 32 H 0.09 0.54 8.14 -51.93 -0.42 0.11 16 33 H 0.09 0.49 8.14 -51.93 -0.42 0.07 16 34 H 0.09 0.50 8.14 -51.93 -0.42 0.07 16 35 H 0.08 0.32 6.33 -51.93 -0.33 0.00 16 36 H 0.08 0.48 8.14 -51.93 -0.42 0.06 16 37 H 0.07 0.35 8.14 -51.93 -0.42 -0.07 16 38 H 0.05 0.39 8.14 -51.93 -0.42 -0.04 16 39 H 0.09 0.50 6.14 -51.93 -0.32 0.18 16 40 H 0.08 0.76 8.07 -51.93 -0.42 0.34 16 41 H 0.05 0.49 8.14 -51.93 -0.42 0.07 16 42 H 0.08 1.27 7.98 -51.93 -0.41 0.86 16 43 H 0.04 0.72 8.08 -51.93 -0.42 0.30 16 44 H 0.08 1.27 5.88 -51.93 -0.31 0.97 16 45 H 0.09 1.28 7.41 -51.92 -0.38 0.89 16 46 H 0.27 6.66 8.31 -40.82 -0.34 6.32 16 47 H 0.13 2.52 6.46 -51.93 -0.34 2.19 16 48 H 0.09 1.25 6.67 -51.93 -0.35 0.90 16 49 H 0.05 0.62 8.14 -51.93 -0.42 0.20 16 50 H 0.09 1.13 6.34 -51.93 -0.33 0.80 16 51 H 0.05 0.54 8.14 -51.93 -0.42 0.12 16 52 H 0.07 0.80 6.72 -51.93 -0.35 0.45 16 LS Contribution 420.27 15.07 6.33 6.33 Total: -1.00 -32.36 420.27 -8.64 -41.01 By element: Atomic # 1 Polarization: 8.97 SS G_CDS: -7.70 Total: 1.27 kcal Atomic # 6 Polarization: 17.74 SS G_CDS: -1.65 Total: 16.09 kcal Atomic # 7 Polarization: -65.16 SS G_CDS: -0.88 Total: -66.04 kcal Atomic # 8 Polarization: -13.44 SS G_CDS: 0.29 Total: -13.15 kcal Atomic # 14 Polarization: 19.53 SS G_CDS: -5.03 Total: 14.50 kcal Total LS contribution 6.33 Total: 6.33 kcal Total: -32.36 -8.64 -41.01 kcal The number of atoms in the molecule is 52 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. ZINC000069370378.mol2 52 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 134.867 kcal (2) G-P(sol) polarization free energy of solvation -32.364 kcal (3) G-ENP(sol) elect.-nuc.-pol. free energy of system 102.503 kcal (4) G-CDS(sol) cavity-dispersion-solvent structure free energy -8.643 kcal (5) G-P-CDS(sol) = G-P(sol) + G-CDS(sol) = (2) + (4) -41.007 kcal (6) G-S(sol) free energy of system = (1) + (5) 93.860 kcal 1SCF WAS SPECIFIED, SO GEOMETRY OPTIMIZATION WAS NOT USED Total computer time = 1.56 seconds