Wall clock time and date at job start Sat Apr 11 2020 06:10:12 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.53000 * 1 3 3 O 1.42901 * 109.47348 * 2 1 4 4 C 1.42895 * 114.00431 * 179.97438 * 3 2 1 5 5 H 1.09007 * 115.54476 * 4.41272 * 4 3 2 6 6 C 1.52994 * 117.50268 * 218.74555 * 4 3 2 7 7 C 1.52992 * 60.00141 * 252.48389 * 6 4 3 8 8 H 1.08999 * 117.50077 * 107.48480 * 7 6 4 9 9 C 1.50701 * 117.49668 * 252.51032 * 7 6 4 10 10 O 1.21276 * 120.00504 * 68.62367 * 9 7 6 11 11 N 1.34777 * 119.99880 * 248.61612 * 9 7 6 12 12 C 1.46932 * 120.62955 * 359.72644 * 11 9 7 13 13 C 1.52766 * 109.00466 * 233.57777 * 12 11 9 14 14 C 1.53039 * 109.38598 * 305.58698 * 13 12 11 15 15 N 1.46495 * 109.45989 * 181.25735 * 14 13 12 16 16 C 1.36936 * 120.00512 * 304.73059 * 15 14 13 17 17 H 1.08005 * 119.99506 * 321.82561 * 16 15 14 18 18 C 1.38815 * 120.00492 * 141.82897 * 16 15 14 19 19 N 1.31609 * 120.00055 * 157.03669 * 18 16 15 20 20 Si 1.86798 * 119.99636 * 337.03562 * 18 16 15 21 21 H 1.48505 * 109.47180 * 95.45192 * 20 18 16 22 22 H 1.48496 * 109.46914 * 215.44762 * 20 18 16 23 23 H 1.48498 * 109.47387 * 335.45189 * 20 18 16 24 24 C 1.46504 * 119.99600 * 124.72456 * 15 14 13 25 25 C 1.53005 * 109.46955 * 90.00407 * 24 15 14 26 26 O 1.42901 * 109.46916 * 64.99833 * 25 24 15 27 27 C 1.42904 * 113.99387 * 180.02562 * 26 25 24 28 28 C 1.53047 * 109.53721 * 61.24191 * 14 13 12 29 29 C 1.46925 * 120.63147 * 179.97438 * 11 9 7 30 30 H 1.09000 * 109.47098 * 180.02562 * 1 2 3 31 31 H 1.08996 * 109.47479 * 300.00893 * 1 2 3 32 32 H 1.09008 * 109.46804 * 60.00240 * 1 2 3 33 33 H 1.08990 * 109.47339 * 239.99511 * 2 1 3 34 34 H 1.09003 * 109.46639 * 119.99312 * 2 1 3 35 35 H 1.08999 * 117.49639 * 359.97271 * 6 4 3 36 36 H 1.08996 * 117.50539 * 144.99429 * 6 4 3 37 37 H 1.09001 * 109.58706 * 353.41039 * 12 11 9 38 38 H 1.08994 * 109.58866 * 113.84169 * 12 11 9 39 39 H 1.09000 * 109.56255 * 185.55278 * 13 12 11 40 40 H 1.09003 * 109.38279 * 65.51671 * 13 12 11 41 41 H 1.09001 * 109.44787 * 301.22819 * 14 13 12 42 42 H 0.97001 * 119.99938 * 180.02562 * 19 18 16 43 43 H 1.09001 * 109.46918 * 210.00176 * 24 15 14 44 44 H 1.08993 * 109.47305 * 330.00019 * 24 15 14 45 45 H 1.09000 * 109.46916 * 184.99516 * 25 24 15 46 46 H 1.08998 * 109.46672 * 304.99515 * 25 24 15 47 47 H 1.08998 * 109.46658 * 180.02562 * 27 26 25 48 48 H 1.08999 * 109.46778 * 299.99413 * 27 26 25 49 49 H 1.08992 * 109.47495 * 60.00239 * 27 26 25 50 50 H 1.08994 * 109.49557 * 178.68313 * 28 14 13 51 51 H 1.09002 * 109.52286 * 58.60134 * 28 14 13 52 52 H 1.08997 * 109.58806 * 246.15777 * 29 11 9 53 53 H 1.08993 * 109.58663 * 6.58713 * 29 11 9 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 8 2.0064 1.3473 0.0000 4 6 3.4309 1.4601 0.0006 5 1 3.9770 0.5197 -0.0747 6 6 4.0512 2.5712 0.8500 7 6 4.0430 2.6886 -0.6754 8 1 3.3591 3.4110 -1.1209 9 6 5.3548 2.5059 -1.3944 10 8 5.8848 1.4155 -1.4214 11 7 5.9389 3.5545 -2.0074 12 6 5.3060 4.8803 -1.9886 13 6 5.1745 5.3892 -3.4230 14 6 6.5459 5.3566 -4.1013 15 7 6.4272 5.8725 -5.4673 16 6 5.5284 5.3121 -6.3352 17 1 5.3678 4.2441 -6.3312 18 6 4.8237 6.1167 -7.2201 19 7 4.3244 5.6013 -8.3234 20 14 4.5841 7.9308 -6.8449 21 1 5.6361 8.7217 -7.5326 22 1 3.2469 8.3631 -7.3247 23 1 4.6788 8.1499 -5.3792 24 6 7.2701 6.9881 -5.9048 25 6 8.5709 6.4426 -6.4974 26 8 8.2797 5.7009 -7.6836 27 6 9.4325 5.1446 -8.3189 28 6 7.0620 3.9164 -4.1434 29 6 7.2180 3.3922 -2.7119 30 1 -0.3633 -1.0277 0.0005 31 1 -0.3634 0.5139 0.8898 32 1 -0.3633 0.5138 -0.8901 33 1 1.8933 -0.5139 -0.8898 34 1 1.8933 -0.5138 0.8901 35 1 3.3727 3.2164 1.4081 36 1 5.0054 2.3619 1.3336 37 1 4.3170 4.8067 -1.5363 38 1 5.9225 5.5714 -1.4138 39 1 4.7971 6.4118 -3.4143 40 1 4.4820 4.7505 -3.9714 41 1 7.2440 5.9763 -3.5385 42 1 3.8322 6.1636 -8.9419 43 1 6.7413 7.5680 -6.6612 44 1 7.4991 7.6267 -5.0517 45 1 9.2349 7.2719 -6.7413 46 1 9.0557 5.7904 -5.7709 47 1 9.1273 4.5978 -9.2111 48 1 10.1147 5.9468 -8.6005 49 1 9.9348 4.4654 -7.6302 50 1 8.0281 3.8908 -4.6474 51 1 6.3524 3.2908 -4.6849 52 1 7.9933 3.9597 -2.1972 53 1 7.4902 2.3371 -2.7367 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 ZINC001716717439.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:12 Heat of formation + Delta-G solvation = -63.172566 kcal Electronic energy + Delta-G solvation = -31266.307049 eV Core-core repulsion = 27140.605422 eV Total energy + Delta-G solvation = -4125.701627 eV No. of doubly occupied orbitals = 66 Molecular weight (most abundant/longest-lived isotopes) = 340.220 amu Computer time = 0.98 seconds Orbital eigenvalues (eV) -40.29418 -38.61644 -37.72213 -36.90721 -35.23671 -35.04471 -32.01893 -31.18545 -30.97208 -30.71298 -28.37592 -25.70293 -25.56643 -24.10335 -23.10254 -22.95361 -22.56151 -22.19475 -21.02302 -19.78527 -18.67441 -17.87240 -16.54913 -16.28500 -16.05531 -15.81817 -15.73978 -15.37078 -15.20131 -14.74569 -14.60644 -14.28951 -14.10063 -13.94022 -13.48495 -13.37522 -13.00035 -12.93089 -12.74612 -12.63383 -12.45756 -12.25477 -12.14503 -11.90848 -11.82027 -11.59092 -11.41581 -11.38034 -11.14548 -10.78823 -10.70134 -10.67469 -10.65538 -10.53542 -10.27544 -9.92263 -9.69749 -9.67893 -9.38174 -9.31789 -8.86266 -8.47557 -8.32412 -6.75575 -5.85615 -3.26566 2.53282 2.93199 3.36652 3.42863 3.50457 3.53116 3.75312 4.52945 4.57270 4.60860 4.62849 4.78267 4.84477 4.85903 4.89647 5.05571 5.14554 5.19052 5.31251 5.37378 5.45962 5.47536 5.51770 5.60839 5.65968 5.68095 5.68727 5.72804 5.77087 5.81651 5.87863 5.89650 5.92956 5.95522 6.05753 6.13143 6.14729 6.16500 6.18785 6.29499 6.46410 6.56561 6.74669 6.81167 7.00989 7.01763 7.05082 7.56503 7.98272 8.18155 8.28965 8.51736 9.07223 9.84303 10.20428 11.30561 Molecular weight = 340.22amu Principal moments of inertia in cm(-1) A = 0.013295 B = 0.002907 C = 0.002660 Principal moments of inertia in units of 10**(-40)*gram-cm**2 A = 2105.554118 B = 9628.127352 C =10522.144230 Net atomic charges, atomic populations, and dipole contributions calculated with CM2 Atom NO. Type Charge No. of electrons 1 C -0.152 4.152 2 C 0.058 3.942 3 O -0.356 6.356 4 C 0.052 3.948 5 H 0.124 0.876 6 C -0.155 4.155 7 C -0.192 4.192 8 H 0.133 0.867 9 C 0.547 3.453 10 O -0.553 6.553 11 N -0.598 5.598 12 C 0.107 3.893 13 C -0.141 4.141 14 C 0.180 3.820 15 N -0.593 5.593 16 C -0.247 4.247 17 H 0.074 0.926 18 C 0.045 3.955 19 N -0.896 5.896 20 Si 0.869 3.131 21 H -0.288 1.288 22 H -0.304 1.304 23 H -0.271 1.271 24 C 0.134 3.866 25 C 0.047 3.953 26 O -0.346 6.346 27 C 0.032 3.968 28 C -0.139 4.139 29 C 0.114 3.886 30 H 0.072 0.928 31 H 0.062 0.938 32 H 0.065 0.935 33 H 0.054 0.946 34 H 0.051 0.949 35 H 0.109 0.891 36 H 0.105 0.895 37 H 0.081 0.919 38 H 0.064 0.936 39 H 0.071 0.929 40 H 0.079 0.921 41 H 0.049 0.951 42 H 0.259 0.741 43 H 0.062 0.938 44 H 0.052 0.948 45 H 0.026 0.974 46 H 0.036 0.964 47 H 0.078 0.922 48 H 0.025 0.975 49 H 0.029 0.971 50 H 0.067 0.933 51 H 0.072 0.928 52 H 0.061 0.939 53 H 0.085 0.915 For the charges calculated by CM2: Dipole moment (debyes) X Y Z Total from point charges 1.947 -4.025 16.945 17.525 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.209 4.209 2 C -0.019 4.019 3 O -0.275 6.275 4 C -0.006 4.006 5 H 0.141 0.859 6 C -0.193 4.193 7 C -0.213 4.213 8 H 0.151 0.849 9 C 0.335 3.665 10 O -0.432 6.432 11 N -0.330 5.330 12 C -0.016 4.016 13 C -0.181 4.181 14 C 0.073 3.927 15 N -0.314 5.314 16 C -0.375 4.375 17 H 0.092 0.908 18 C -0.155 4.155 19 N -0.537 5.537 20 Si 0.700 3.300 21 H -0.217 1.217 22 H -0.233 1.233 23 H -0.196 1.196 24 C 0.008 3.992 25 C -0.032 4.032 26 O -0.263 6.263 27 C -0.066 4.066 28 C -0.179 4.179 29 C -0.009 4.009 30 H 0.091 0.909 31 H 0.081 0.919 32 H 0.084 0.916 33 H 0.072 0.928 34 H 0.069 0.931 35 H 0.128 0.872 36 H 0.124 0.876 37 H 0.100 0.900 38 H 0.083 0.917 39 H 0.089 0.911 40 H 0.098 0.902 41 H 0.068 0.932 42 H 0.069 0.931 43 H 0.080 0.920 44 H 0.071 0.929 45 H 0.044 0.956 46 H 0.054 0.946 47 H 0.096 0.904 48 H 0.043 0.957 49 H 0.047 0.953 50 H 0.085 0.915 51 H 0.091 0.909 52 H 0.079 0.921 53 H 0.104 0.896 Dipole moment (debyes) X Y Z Total from point charges 1.796 -4.426 16.683 17.354 hybrid contribution 0.328 1.113 -0.138 1.168 sum 2.123 -3.313 16.546 17.007 Atomic orbital electron populations 1.21759 0.93522 1.02389 1.03256 1.22568 0.95550 0.83226 1.00544 1.87814 1.18554 1.27130 1.93963 1.22578 0.81970 0.94360 1.01725 0.85898 1.23411 1.02862 1.02156 0.90881 1.22715 0.97805 1.04384 0.96405 0.84945 1.19482 0.84913 0.83699 0.78378 1.90649 1.65398 1.25254 1.61933 1.48033 1.22956 1.09514 1.52462 1.21833 0.97930 0.82889 0.98921 1.22250 1.00769 1.01888 0.93154 1.20316 0.93786 0.92590 0.85987 1.44098 1.41481 1.38246 1.07567 1.18455 1.14635 0.92190 1.12239 0.90809 1.25338 0.97088 0.96894 0.96209 1.69887 1.38610 1.36070 1.09181 0.87055 0.75991 0.88209 0.78743 1.21676 1.23275 1.19573 1.21343 0.90261 0.90672 0.96884 1.23278 0.98920 0.93628 0.87422 1.88313 1.28647 1.69981 1.39394 1.22636 0.89622 0.96737 0.97567 1.22380 1.02225 0.99380 0.93950 1.21701 0.83877 1.02337 0.92962 0.90934 0.91884 0.91607 0.92799 0.93075 0.87226 0.87647 0.90036 0.91740 0.91075 0.90201 0.93238 0.93074 0.91968 0.92936 0.95602 0.94587 0.90364 0.95659 0.95270 0.91452 0.90899 0.92093 0.89642 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.15 -1.20 10.28 37.16 0.38 -0.82 16 2 C 0.06 0.55 6.53 37.16 0.24 0.79 16 3 O -0.36 -4.32 10.60 -35.23 -0.37 -4.70 16 4 C 0.05 0.62 5.95 -26.74 -0.16 0.47 16 5 H 0.12 1.49 7.45 -51.93 -0.39 1.10 16 6 C -0.16 -1.57 10.18 -26.74 -0.27 -1.84 16 7 C -0.19 -2.41 5.08 -91.78 -0.47 -2.88 16 8 H 0.13 1.62 6.55 -51.93 -0.34 1.28 16 9 C 0.55 8.40 7.48 -10.98 -0.08 8.32 16 10 O -0.55 -10.05 15.67 5.56 0.09 -9.96 16 11 N -0.60 -8.78 2.98 -173.00 -0.52 -9.29 16 12 C 0.11 1.41 6.32 -3.93 -0.02 1.39 16 13 C -0.14 -2.38 5.21 -26.83 -0.14 -2.52 16 14 C 0.18 3.23 2.60 -67.89 -0.18 3.06 16 15 N -0.59 -13.07 1.98 -165.61 -0.33 -13.40 16 16 C -0.25 -6.58 8.05 -15.06 -0.12 -6.70 16 17 H 0.07 2.20 7.24 -52.48 -0.38 1.82 16 18 C 0.04 1.27 5.29 -17.43 -0.09 1.18 16 19 N -0.90 -27.40 13.97 55.49 0.78 -26.63 16 20 Si 0.87 21.77 25.12 -169.99 -4.27 17.50 16 21 H -0.29 -7.18 6.67 56.52 0.38 -6.80 16 22 H -0.30 -7.86 7.11 56.52 0.40 -7.45 16 23 H -0.27 -6.42 7.11 56.52 0.40 -6.02 16 24 C 0.13 2.78 4.93 -4.04 -0.02 2.76 16 25 C 0.05 0.89 6.28 37.16 0.23 1.12 16 26 O -0.35 -7.62 10.18 -35.23 -0.36 -7.97 16 27 C 0.03 0.53 11.01 101.05 1.11 1.65 16 28 C -0.14 -2.41 5.38 -26.58 -0.14 -2.56 16 29 C 0.11 1.69 6.44 -3.69 -0.02 1.67 16 30 H 0.07 0.44 8.14 -51.93 -0.42 0.01 16 31 H 0.06 0.48 8.14 -51.93 -0.42 0.06 16 32 H 0.06 0.57 8.14 -51.92 -0.42 0.14 16 33 H 0.05 0.51 8.14 -51.93 -0.42 0.09 16 34 H 0.05 0.41 8.14 -51.93 -0.42 -0.01 16 35 H 0.11 0.96 8.14 -51.93 -0.42 0.53 16 36 H 0.11 1.00 8.14 -51.93 -0.42 0.58 16 37 H 0.08 0.94 5.51 -51.93 -0.29 0.66 16 38 H 0.06 0.73 8.14 -51.93 -0.42 0.31 16 39 H 0.07 1.23 8.14 -51.93 -0.42 0.81 16 40 H 0.08 1.58 7.61 -51.93 -0.40 1.19 16 41 H 0.05 0.79 7.77 -51.93 -0.40 0.39 16 42 H 0.26 7.60 8.31 -40.82 -0.34 7.26 16 43 H 0.06 1.46 3.40 -51.93 -0.18 1.29 16 44 H 0.05 0.92 7.88 -51.93 -0.41 0.51 16 45 H 0.03 0.43 8.14 -51.93 -0.42 0.01 16 46 H 0.04 0.63 8.14 -51.93 -0.42 0.21 16 47 H 0.08 1.27 8.14 -51.93 -0.42 0.85 16 48 H 0.02 0.37 8.14 -51.93 -0.42 -0.06 16 49 H 0.03 0.44 8.14 -51.93 -0.42 0.02 16 50 H 0.07 1.17 8.14 -51.93 -0.42 0.75 16 51 H 0.07 1.50 7.24 -51.93 -0.38 1.13 16 52 H 0.06 0.79 8.14 -51.93 -0.42 0.36 16 53 H 0.09 1.35 7.00 -51.93 -0.36 0.99 16 LS Contribution 414.62 15.07 6.25 6.25 Total: -1.00 -33.21 414.62 -7.93 -41.14 By element: Atomic # 1 Polarization: 11.43 SS G_CDS: -9.44 Total: 1.99 kcal Atomic # 6 Polarization: 4.83 SS G_CDS: 0.25 Total: 5.08 kcal Atomic # 7 Polarization: -49.25 SS G_CDS: -0.07 Total: -49.31 kcal Atomic # 8 Polarization: -21.99 SS G_CDS: -0.64 Total: -22.64 kcal Atomic # 14 Polarization: 21.77 SS G_CDS: -4.27 Total: 17.50 kcal Total LS contribution 6.25 Total: 6.25 kcal Total: -33.21 -7.93 -41.14 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. ZINC001716717439.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 -22.032 kcal (2) G-P(sol) polarization free energy of solvation -33.214 kcal (3) G-ENP(sol) elect.-nuc.-pol. free energy of system -55.246 kcal (4) G-CDS(sol) cavity-dispersion-solvent structure free energy -7.927 kcal (5) G-P-CDS(sol) = G-P(sol) + G-CDS(sol) = (2) + (4) -41.140 kcal (6) G-S(sol) free energy of system = (1) + (5) -63.173 kcal 1SCF WAS SPECIFIED, SO GEOMETRY OPTIMIZATION WAS NOT USED Total computer time = 0.98 seconds