Wall clock time and date at job start Sat Apr 11 2020 06:12:32 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 N 1.46507 * 1 3 3 C 1.34771 * 119.99555 * 2 1 4 4 O 1.21592 * 120.00433 * 359.97438 * 3 2 1 5 5 N 1.34783 * 120.00109 * 179.97438 * 3 2 1 6 6 C 1.46494 * 120.00223 * 180.02562 * 5 3 2 7 7 C 1.50700 * 109.46998 * 179.97438 * 6 5 3 8 8 O 1.21282 * 120.00058 * 0.02562 * 7 6 5 9 9 N 1.34777 * 119.99997 * 179.97438 * 7 6 5 10 10 C 1.46925 * 120.62980 * 180.02562 * 9 7 6 11 11 C 1.53187 * 108.77688 * 126.41364 * 10 9 7 12 12 C 1.53046 * 109.31043 * 54.63722 * 11 10 9 13 13 C 1.53001 * 109.46008 * 178.63085 * 12 11 10 14 14 C 1.53000 * 109.46981 * 300.02190 * 13 12 11 15 15 C 1.52997 * 109.47440 * 60.02370 * 13 12 11 16 16 C 1.53000 * 109.47072 * 180.02562 * 13 12 11 17 17 N 1.46499 * 109.47170 * 179.97438 * 16 13 12 18 18 C 1.36932 * 120.00350 * 179.97438 * 17 16 13 19 19 H 1.08005 * 119.99868 * 359.97438 * 18 17 16 20 20 C 1.38810 * 120.00619 * 180.02562 * 18 17 16 21 21 N 1.31622 * 119.99828 * 359.97438 * 20 18 17 22 22 Si 1.86800 * 120.00512 * 180.02562 * 20 18 17 23 23 H 1.48501 * 109.47357 * 89.99384 * 22 20 18 24 24 H 1.48494 * 109.47219 * 209.99807 * 22 20 18 25 25 H 1.48503 * 109.46989 * 329.99705 * 22 20 18 26 26 C 1.53043 * 109.53475 * 298.63386 * 12 11 10 27 27 C 1.46921 * 120.63225 * 0.02846 * 9 7 6 28 28 H 1.09001 * 109.46583 * 60.00726 * 1 2 3 29 29 H 1.09002 * 109.46681 * 179.97438 * 1 2 3 30 30 H 1.08994 * 109.47203 * 300.00488 * 1 2 3 31 31 H 0.96995 * 120.00195 * 179.97438 * 2 1 3 32 32 H 0.97002 * 119.99057 * 359.97438 * 5 3 2 33 33 H 1.08997 * 109.47249 * 59.99892 * 6 5 3 34 34 H 1.08999 * 109.47356 * 299.99870 * 6 5 3 35 35 H 1.09003 * 109.59222 * 6.62634 * 10 9 7 36 36 H 1.09006 * 109.70276 * 246.27263 * 10 9 7 37 37 H 1.08996 * 109.50293 * 294.68191 * 11 10 9 38 38 H 1.09001 * 109.50064 * 174.58451 * 11 10 9 39 39 H 1.09002 * 109.46068 * 58.65418 * 12 11 10 40 40 H 1.09005 * 109.47243 * 60.00407 * 14 13 12 41 41 H 1.09000 * 109.47231 * 180.02562 * 14 13 12 42 42 H 1.09002 * 109.47199 * 300.00391 * 14 13 12 43 43 H 1.08995 * 109.47302 * 60.00023 * 15 13 12 44 44 H 1.09000 * 109.47089 * 180.02562 * 15 13 12 45 45 H 1.09005 * 109.47083 * 300.00337 * 15 13 12 46 46 H 1.09008 * 109.46946 * 299.99481 * 16 13 12 47 47 H 1.08992 * 109.47337 * 59.99460 * 16 13 12 48 48 H 0.97007 * 120.00299 * 359.97438 * 17 16 13 49 49 H 0.96993 * 120.00331 * 180.02562 * 21 20 18 50 50 H 1.09004 * 109.49679 * 181.31796 * 26 12 11 51 51 H 1.08994 * 109.49744 * 301.41312 * 26 12 11 52 52 H 1.09002 * 109.58756 * 353.36576 * 27 9 7 53 53 H 1.09002 * 109.59282 * 113.79441 * 27 9 7 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.4651 0.0000 0.0000 3 6 2.1388 1.1672 0.0000 4 8 1.5309 2.2202 -0.0005 5 7 3.4867 1.1673 0.0005 6 6 4.2191 2.4360 0.0011 7 6 5.7007 2.1607 0.0022 8 8 6.1050 1.0173 0.0021 9 7 6.5765 3.1852 0.0024 10 6 8.0239 2.9326 0.0029 11 6 8.6580 3.6959 -1.1641 12 6 8.2781 5.1752 -1.0674 13 6 8.9376 5.9468 -2.2122 14 6 10.4577 5.8026 -2.1152 15 6 8.4563 5.3835 -3.5508 16 6 8.5584 7.4259 -2.1154 17 7 9.1905 8.1648 -3.2112 18 6 8.9929 9.5145 -3.3305 19 1 8.3707 10.0343 -2.6169 20 6 9.5913 10.2146 -4.3691 21 7 10.3491 9.5811 -5.2391 22 14 9.3223 12.0560 -4.5314 23 1 10.3725 12.7787 -3.7699 24 1 9.3946 12.4425 -5.9633 25 1 7.9857 12.4106 -3.9899 26 6 6.7576 5.3205 -1.1645 27 6 6.1020 4.5757 0.0021 28 1 -0.3632 0.5137 -0.8901 29 1 -0.3633 -1.0277 -0.0005 30 1 -0.3633 0.5139 0.8899 31 1 1.9501 -0.8400 0.0004 32 1 3.9716 0.3272 0.0013 33 1 3.9554 3.0074 0.8911 34 1 3.9561 3.0078 -0.8889 35 1 8.2087 1.8649 -0.1153 36 1 8.4565 3.2764 0.9425 37 1 8.2937 3.2875 -2.1067 38 1 9.7424 3.5951 -1.1184 39 1 8.6211 5.5767 -0.1138 40 1 10.7279 4.7488 -2.1842 41 1 10.9277 6.3526 -2.9305 42 1 10.8006 6.2039 -1.1615 43 1 7.3734 5.4862 -3.6200 44 1 8.9266 5.9329 -4.3664 45 1 8.7264 4.3297 -3.6198 46 1 7.4754 7.5286 -2.1846 47 1 8.9013 7.8272 -1.1618 48 1 9.7490 7.6979 -3.8524 49 1 10.7675 10.0702 -5.9646 50 1 6.4894 6.3760 -1.1191 51 1 6.4111 4.8969 -2.1071 52 1 5.0186 4.5921 -0.1169 53 1 6.3758 5.0555 0.9417 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 ZINC001101546977.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:32 Heat of formation + Delta-G solvation = -59.042311 kcal Electronic energy + Delta-G solvation = -30543.692380 eV Core-core repulsion = 26450.462606 eV Total energy + Delta-G solvation = -4093.229774 eV No. of doubly occupied orbitals = 66 Molecular weight (most abundant/longest-lived isotopes) = 340.226 amu Computer time = 1.87 seconds Orbital eigenvalues (eV) -40.87608 -40.05091 -38.84595 -36.34202 -35.08609 -34.66728 -33.81135 -31.81751 -31.33438 -30.03309 -28.67297 -26.90740 -26.29812 -25.94309 -24.84980 -23.21425 -22.51648 -21.38322 -21.03262 -19.64140 -19.01585 -18.05708 -17.28836 -17.07344 -16.81001 -16.37226 -15.85892 -15.51457 -14.82007 -14.69851 -14.65639 -14.58596 -14.26066 -14.12045 -13.91250 -13.72766 -13.63855 -13.11507 -13.09401 -12.80692 -12.63850 -12.61596 -12.58064 -12.31226 -11.83386 -11.67875 -11.51950 -11.19825 -10.96764 -10.88234 -10.79782 -10.69272 -10.63059 -10.29199 -10.18147 -10.14170 -9.98176 -9.71898 -9.67339 -9.53052 -9.26829 -8.89548 -8.82896 -6.86911 -5.74660 -3.01899 2.15789 2.22256 2.79709 2.86260 3.44299 3.49626 3.51123 3.91621 4.05404 4.57963 4.60884 4.75578 4.78042 4.80091 4.95047 4.97434 5.01444 5.13611 5.16165 5.20648 5.32531 5.41460 5.45653 5.49312 5.56553 5.60968 5.63923 5.75962 5.87337 5.90589 5.99169 6.00910 6.03618 6.07271 6.11488 6.18281 6.31609 6.34954 6.38003 6.39318 6.50661 6.59863 6.72837 6.77003 6.96250 6.98324 7.05358 7.61096 7.71773 7.75626 8.22740 8.45200 9.56213 10.16147 10.50239 11.48421 Molecular weight = 340.23amu Principal moments of inertia in cm(-1) A = 0.016403 B = 0.002081 C = 0.001916 Principal moments of inertia in units of 10**(-40)*gram-cm**2 A = 1706.593582 B =13453.858022 C =14607.834879 Net atomic charges, atomic populations, and dipole contributions calculated with CM2 Atom NO. Type Charge No. of electrons 1 C 0.107 3.893 2 N -0.750 5.750 3 C 0.705 3.295 4 O -0.590 6.590 5 N -0.735 5.735 6 C 0.120 3.880 7 C 0.512 3.488 8 O -0.546 6.546 9 N -0.601 5.601 10 C 0.112 3.888 11 C -0.132 4.132 12 C -0.072 4.072 13 C -0.065 4.065 14 C -0.138 4.138 15 C -0.138 4.138 16 C 0.162 3.838 17 N -0.729 5.729 18 C -0.236 4.236 19 H 0.069 0.931 20 C -0.015 4.015 21 N -0.934 5.934 22 Si 0.905 3.095 23 H -0.291 1.291 24 H -0.292 1.292 25 H -0.283 1.283 26 C -0.136 4.136 27 C 0.101 3.899 28 H 0.059 0.941 29 H 0.075 0.925 30 H 0.058 0.942 31 H 0.399 0.601 32 H 0.407 0.593 33 H 0.098 0.902 34 H 0.100 0.900 35 H 0.090 0.910 36 H 0.068 0.932 37 H 0.071 0.929 38 H 0.080 0.920 39 H 0.071 0.929 40 H 0.046 0.954 41 H 0.066 0.934 42 H 0.049 0.951 43 H 0.051 0.949 44 H 0.065 0.935 45 H 0.044 0.956 46 H 0.020 0.980 47 H 0.015 0.985 48 H 0.384 0.616 49 H 0.254 0.746 50 H 0.087 0.913 51 H 0.072 0.928 52 H 0.086 0.914 53 H 0.075 0.925 For the charges calculated by CM2: Dipole moment (debyes) X Y Z Total from point charges -12.480 -21.612 13.009 28.144 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.038 4.038 2 N -0.406 5.406 3 C 0.405 3.595 4 O -0.468 6.468 5 N -0.394 5.394 6 C -0.008 4.008 7 C 0.299 3.701 8 O -0.425 6.425 9 N -0.334 5.334 10 C -0.010 4.010 11 C -0.170 4.170 12 C -0.091 4.091 13 C -0.065 4.065 14 C -0.195 4.195 15 C -0.196 4.196 16 C 0.035 3.965 17 N -0.374 5.374 18 C -0.367 4.367 19 H 0.087 0.913 20 C -0.211 4.211 21 N -0.583 5.583 22 Si 0.737 3.263 23 H -0.218 1.218 24 H -0.219 1.219 25 H -0.208 1.208 26 C -0.174 4.174 27 C -0.021 4.021 28 H 0.078 0.922 29 H 0.094 0.906 30 H 0.077 0.923 31 H 0.232 0.768 32 H 0.244 0.756 33 H 0.116 0.884 34 H 0.118 0.882 35 H 0.108 0.892 36 H 0.087 0.913 37 H 0.089 0.911 38 H 0.098 0.902 39 H 0.089 0.911 40 H 0.065 0.935 41 H 0.085 0.915 42 H 0.068 0.932 43 H 0.070 0.930 44 H 0.084 0.916 45 H 0.063 0.937 46 H 0.038 0.962 47 H 0.033 0.967 48 H 0.219 0.781 49 H 0.064 0.936 50 H 0.106 0.894 51 H 0.091 0.909 52 H 0.104 0.896 53 H 0.093 0.907 Dipole moment (debyes) X Y Z Total from point charges -12.900 -21.054 13.031 27.919 hybrid contribution 0.561 0.316 -0.384 0.750 sum -12.339 -20.738 12.647 27.244 Atomic orbital electron populations 1.21414 0.77872 1.02051 1.02418 1.44873 1.07853 1.07108 1.80721 1.15840 0.81344 0.84454 0.77898 1.90954 1.67878 1.30215 1.57738 1.44743 1.04688 1.09618 1.80356 1.20935 0.88229 0.86760 1.04865 1.20452 0.88388 0.85084 0.76156 1.90669 1.79743 1.21069 1.50971 1.48068 1.06850 1.06372 1.72122 1.21692 0.78516 1.01799 0.99008 1.22099 1.01053 0.95482 0.98358 1.21101 0.94393 0.94416 0.99187 1.20684 0.95412 0.96018 0.94376 1.21851 0.94960 1.00279 1.02397 1.21893 1.00948 0.99413 0.97310 1.20734 0.95485 0.88483 0.91839 1.42430 1.58760 1.05612 1.30639 1.19229 1.24344 0.84547 1.08611 0.91344 1.24716 1.00433 0.97919 0.98053 1.69651 1.36619 1.34773 1.17298 0.86210 0.77623 0.85677 0.76810 1.21834 1.21854 1.20833 1.22131 0.95842 1.01420 0.98016 1.21956 0.99533 0.80532 1.00099 0.92243 0.90594 0.92288 0.76789 0.75565 0.88423 0.88234 0.89159 0.91337 0.91082 0.90174 0.91077 0.93541 0.91522 0.93191 0.92978 0.91638 0.93672 0.96223 0.96712 0.78057 0.93602 0.89424 0.90902 0.89606 0.90664 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.11 0.60 10.70 59.85 0.64 1.24 16 2 N -0.75 -5.14 5.72 -72.06 -0.41 -5.55 16 3 C 0.71 7.34 8.57 -86.92 -0.74 6.59 16 4 O -0.59 -8.43 16.34 5.28 0.09 -8.35 16 5 N -0.74 -6.66 5.53 -66.42 -0.37 -7.03 16 6 C 0.12 1.08 4.50 -5.20 -0.02 1.06 16 7 C 0.51 5.26 7.77 -10.98 -0.09 5.18 16 8 O -0.55 -6.98 15.08 5.56 0.08 -6.89 16 9 N -0.60 -5.66 2.97 -172.80 -0.51 -6.17 16 10 C 0.11 1.06 6.43 -3.71 -0.02 1.04 16 11 C -0.13 -1.34 4.65 -26.61 -0.12 -1.46 16 12 C -0.07 -0.82 1.51 -90.57 -0.14 -0.95 16 13 C -0.06 -0.96 0.33 -154.51 -0.05 -1.01 16 14 C -0.14 -2.08 7.95 37.16 0.30 -1.78 16 15 C -0.14 -2.12 7.60 37.16 0.28 -1.83 16 16 C 0.16 3.09 4.80 -4.04 -0.02 3.07 16 17 N -0.73 -17.88 5.00 -59.91 -0.30 -18.18 16 18 C -0.24 -6.52 9.99 -15.06 -0.15 -6.67 16 19 H 0.07 1.84 7.83 -52.48 -0.41 1.43 16 20 C -0.02 -0.44 5.50 -17.43 -0.10 -0.54 16 21 N -0.93 -28.06 13.06 55.49 0.72 -27.34 16 22 Si 0.91 22.40 29.55 -169.99 -5.02 17.38 16 23 H -0.29 -7.18 7.11 56.52 0.40 -6.77 16 24 H -0.29 -7.30 7.11 56.52 0.40 -6.90 16 25 H -0.28 -6.81 7.11 56.52 0.40 -6.41 16 26 C -0.14 -1.42 4.65 -26.63 -0.12 -1.54 16 27 C 0.10 0.84 6.34 -3.74 -0.02 0.81 16 28 H 0.06 0.39 8.14 -51.93 -0.42 -0.03 16 29 H 0.08 0.19 8.14 -51.93 -0.42 -0.24 16 30 H 0.06 0.38 8.14 -51.93 -0.42 -0.05 16 31 H 0.40 1.60 8.84 -40.82 -0.36 1.24 16 32 H 0.41 3.30 7.45 -40.82 -0.30 2.99 16 33 H 0.10 0.81 7.94 -51.93 -0.41 0.40 16 34 H 0.10 0.89 7.69 -51.93 -0.40 0.49 16 35 H 0.09 0.94 7.01 -51.93 -0.36 0.57 16 36 H 0.07 0.58 8.14 -51.93 -0.42 0.16 16 37 H 0.07 0.79 6.51 -51.93 -0.34 0.45 16 38 H 0.08 0.78 6.36 -51.93 -0.33 0.45 16 39 H 0.07 0.80 8.14 -51.93 -0.42 0.38 16 40 H 0.05 0.59 6.40 -51.93 -0.33 0.26 16 41 H 0.07 1.21 6.86 -51.93 -0.36 0.85 16 42 H 0.05 0.74 8.14 -51.93 -0.42 0.31 16 43 H 0.05 0.78 6.46 -51.93 -0.34 0.44 16 44 H 0.06 1.20 6.86 -51.93 -0.36 0.84 16 45 H 0.04 0.60 6.46 -51.93 -0.34 0.27 16 46 H 0.02 0.38 6.40 -51.92 -0.33 0.05 16 47 H 0.02 0.29 8.14 -51.93 -0.42 -0.13 16 48 H 0.38 9.90 5.36 -40.82 -0.22 9.68 16 49 H 0.25 7.41 8.31 -40.82 -0.34 7.07 16 50 H 0.09 0.98 6.36 -51.93 -0.33 0.65 16 51 H 0.07 0.81 6.51 -51.93 -0.34 0.48 16 52 H 0.09 0.59 5.93 -51.93 -0.31 0.29 16 53 H 0.08 0.56 8.14 -51.93 -0.42 0.13 16 LS Contribution 402.54 15.07 6.07 6.07 Total: -1.00 -34.78 402.54 -8.72 -43.50 By element: Atomic # 1 Polarization: 18.05 SS G_CDS: -8.68 Total: 9.37 kcal Atomic # 6 Polarization: 3.58 SS G_CDS: -0.38 Total: 3.20 kcal Atomic # 7 Polarization: -63.40 SS G_CDS: -0.87 Total: -64.27 kcal Atomic # 8 Polarization: -15.41 SS G_CDS: 0.17 Total: -15.24 kcal Atomic # 14 Polarization: 22.40 SS G_CDS: -5.02 Total: 17.38 kcal Total LS contribution 6.07 Total: 6.07 kcal Total: -34.78 -8.72 -43.50 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. ZINC001101546977.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 -15.543 kcal (2) G-P(sol) polarization free energy of solvation -34.780 kcal (3) G-ENP(sol) elect.-nuc.-pol. free energy of system -50.323 kcal (4) G-CDS(sol) cavity-dispersion-solvent structure free energy -8.719 kcal (5) G-P-CDS(sol) = G-P(sol) + G-CDS(sol) = (2) + (4) -43.499 kcal (6) G-S(sol) free energy of system = (1) + (5) -59.042 kcal 1SCF WAS SPECIFIED, SO GEOMETRY OPTIMIZATION WAS NOT USED Total computer time = 1.87 seconds