Wall clock time and date at job start Sat Apr 11 2020 06:10:46 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.52998 * 1 3 3 O 1.42895 * 109.46859 * 2 1 4 4 C 1.42908 * 114.00017 * 175.72484 * 3 2 1 5 5 H 1.08992 * 109.46861 * 30.60139 * 4 3 2 6 6 C 1.52997 * 109.47218 * 150.60331 * 4 3 2 7 7 C 1.50704 * 109.47231 * 295.00109 * 6 4 3 8 8 O 1.21286 * 119.99772 * 359.97438 * 7 6 4 9 9 N 1.34770 * 120.00103 * 180.02562 * 7 6 4 10 10 C 1.46500 * 119.99997 * 180.02562 * 9 7 6 11 11 H 1.09000 * 110.77595 * 31.45903 * 10 9 7 12 12 C 1.55053 * 110.87362 * 267.70128 * 10 9 7 13 13 O 1.44567 * 103.64534 * 205.24949 * 12 10 9 14 14 C 1.42969 * 106.77763 * 39.22755 * 13 12 10 15 15 C 1.53212 * 115.36135 * 203.79814 * 14 13 12 16 16 N 1.47727 * 85.89820 * 220.05736 * 15 14 13 17 17 C 1.36935 * 134.43312 * 204.09841 * 16 15 14 18 18 H 1.08000 * 119.99710 * 180.02562 * 17 16 15 19 19 C 1.38813 * 119.99970 * 0.02562 * 17 16 15 20 20 N 1.31613 * 120.00050 * 0.02562 * 19 17 16 21 21 Si 1.86800 * 119.99698 * 179.97438 * 19 17 16 22 22 H 1.48500 * 109.47146 * 60.00438 * 21 19 17 23 23 H 1.48505 * 109.47051 * 180.02562 * 21 19 17 24 24 H 1.48493 * 109.47527 * 300.00125 * 21 19 17 25 25 C 1.47721 * 91.12914 * 24.11533 * 16 15 14 26 26 C 1.53443 * 108.16801 * 334.19415 * 14 13 12 27 27 C 1.53005 * 109.46958 * 270.60135 * 4 3 2 28 28 C 1.53001 * 109.46611 * 174.61783 * 27 4 3 29 29 C 1.52995 * 109.47202 * 294.61821 * 27 4 3 30 30 H 1.09000 * 109.47589 * 300.00141 * 1 2 3 31 31 H 1.08996 * 109.47201 * 60.00028 * 1 2 3 32 32 H 1.09004 * 109.46889 * 180.02562 * 1 2 3 33 33 H 1.09001 * 109.47066 * 239.99802 * 2 1 3 34 34 H 1.09000 * 109.46850 * 119.99315 * 2 1 3 35 35 H 1.08999 * 109.47529 * 175.00176 * 6 4 3 36 36 H 1.08999 * 109.47154 * 54.99666 * 6 4 3 37 37 H 0.96995 * 120.00232 * 359.97438 * 9 7 6 38 38 H 1.09005 * 110.59277 * 86.70609 * 12 10 9 39 39 H 1.09006 * 110.59037 * 323.79077 * 12 10 9 40 40 H 1.09002 * 113.80500 * 334.29393 * 15 14 13 41 41 H 1.09000 * 113.80957 * 105.81106 * 15 14 13 42 42 H 0.97005 * 120.00136 * 180.02562 * 20 19 17 43 43 H 1.09005 * 113.80160 * 221.63456 * 25 16 15 44 44 H 1.09002 * 113.90958 * 90.14979 * 25 16 15 45 45 H 1.08999 * 110.55118 * 120.25958 * 26 14 13 46 46 H 1.09005 * 110.54627 * 242.93493 * 26 14 13 47 47 H 1.08999 * 109.46877 * 54.62115 * 27 4 3 48 48 H 1.08999 * 109.47229 * 60.00548 * 28 27 4 49 49 H 1.09003 * 109.46439 * 179.97438 * 28 27 4 50 50 H 1.08995 * 109.47049 * 299.99923 * 28 27 4 51 51 H 1.09004 * 109.47008 * 173.97964 * 29 27 4 52 52 H 1.08997 * 109.47429 * 293.98306 * 29 27 4 53 53 H 1.09009 * 109.47177 * 53.98094 * 29 27 4 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.0062 1.3473 0.0000 4 6 3.4274 1.4614 0.0973 5 1 3.8944 0.5960 -0.3728 6 6 3.8864 2.7368 -0.6121 7 6 3.5956 2.6240 -2.0865 8 8 3.0730 1.6219 -2.5266 9 7 3.9170 3.6363 -2.9162 10 6 3.6337 3.5269 -4.3493 11 1 2.7430 2.9229 -4.5224 12 6 4.8545 2.9495 -5.1113 13 8 4.6886 3.4562 -6.4550 14 6 4.2145 4.7996 -6.3339 15 6 3.4448 5.3356 -7.5454 16 7 4.1632 6.6041 -7.3069 17 6 3.9166 7.9169 -7.6087 18 1 4.6218 8.6802 -7.3145 19 6 2.7604 8.2664 -8.2926 20 7 1.9009 7.3362 -8.6506 21 14 2.4237 10.0573 -8.7036 22 1 3.5187 10.5834 -9.5577 23 1 1.1333 10.1642 -9.4307 24 1 2.3513 10.8480 -7.4488 25 6 5.2543 5.8897 -6.6131 26 6 3.4832 4.9296 -4.9912 27 6 3.8344 1.5228 1.5710 28 6 5.3605 1.5097 1.6795 29 6 3.2896 2.8080 2.1973 30 1 -0.3634 0.5138 0.8899 31 1 -0.3633 0.5138 -0.8899 32 1 -0.3633 -1.0277 0.0005 33 1 1.8933 -0.5139 -0.8900 34 1 1.8933 -0.5137 0.8901 35 1 4.9578 2.8705 -0.4624 36 1 3.3524 3.5935 -0.2010 37 1 4.3352 4.4376 -2.5643 38 1 4.8289 1.8597 -5.1066 39 1 5.7852 3.3138 -4.6763 40 1 3.7095 4.8501 -8.4847 41 1 2.3683 5.4033 -7.3880 42 1 1.0930 7.5805 -9.1289 43 1 5.6119 6.3967 -5.7168 44 1 6.0616 5.5695 -7.2718 45 1 3.9579 5.6881 -4.3688 46 1 2.4317 5.1693 -5.1499 47 1 3.4253 0.6604 2.0971 48 1 5.7696 2.3722 1.1533 49 1 5.6504 1.5540 2.7293 50 1 5.7486 0.5942 1.2332 51 1 3.4866 2.8015 3.2694 52 1 3.7790 3.6695 1.7431 53 1 2.2148 2.8685 2.0258 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 ZINC001755061285.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:46 Heat of formation + Delta-G solvation = -82.536106 kcal Electronic energy + Delta-G solvation = -30786.287040 eV Core-core repulsion = 26659.745748 eV Total energy + Delta-G solvation = -4126.541293 eV No. of doubly occupied orbitals = 66 Molecular weight (most abundant/longest-lived isotopes) = 340.220 amu Computer time = 0.82 seconds Orbital eigenvalues (eV) -40.40551 -39.61301 -38.36892 -36.21088 -35.60265 -34.43833 -32.07538 -31.64587 -30.79323 -29.36517 -28.00120 -27.52528 -25.70383 -24.95646 -24.77921 -22.96819 -22.13312 -21.60831 -20.64957 -19.64233 -18.77486 -17.63410 -17.05916 -16.77184 -16.39325 -15.78471 -15.42887 -15.31710 -15.06159 -14.93255 -14.52691 -14.24359 -14.09069 -13.98303 -13.53106 -13.46088 -13.37727 -12.98755 -12.93409 -12.63843 -12.43055 -12.41448 -12.14572 -11.99581 -11.81905 -11.56483 -11.51452 -11.36681 -11.18896 -11.12934 -11.01937 -10.93034 -10.75660 -10.68598 -10.56441 -10.40771 -10.08171 -9.81792 -9.60026 -9.55446 -9.06781 -8.87203 -8.74043 -6.92373 -5.79128 -3.12092 2.62852 3.06302 3.38464 3.42409 3.44275 3.50353 4.26260 4.32619 4.41504 4.52168 4.67137 4.74076 4.86235 4.89105 4.95522 4.97974 5.08029 5.13444 5.13951 5.18934 5.23905 5.26286 5.28135 5.31789 5.38119 5.45266 5.47783 5.55347 5.56713 5.59742 5.66241 5.73373 5.82085 5.90245 5.93544 5.96394 5.97466 6.02252 6.13487 6.30286 6.39464 6.46178 6.60726 6.67753 7.01810 7.12149 7.22001 7.46826 7.49904 7.62330 7.94144 8.21729 9.01617 9.67694 10.55870 11.38247 Molecular weight = 340.22amu Principal moments of inertia in cm(-1) A = 0.016042 B = 0.002196 C = 0.002181 Principal moments of inertia in units of 10**(-40)*gram-cm**2 A = 1745.025609 B =12747.456380 C =12836.470562 Net atomic charges, atomic populations, and dipole contributions calculated with CM2 Atom NO. Type Charge No. of electrons 1 C -0.150 4.150 2 C 0.063 3.937 3 O -0.370 6.370 4 C 0.113 3.887 5 H 0.067 0.933 6 C -0.144 4.144 7 C 0.515 3.485 8 O -0.532 6.532 9 N -0.720 5.720 10 C 0.110 3.890 11 H 0.105 0.895 12 C 0.041 3.959 13 O -0.361 6.361 14 C 0.051 3.949 15 C 0.176 3.824 16 N -0.604 5.604 17 C -0.223 4.223 18 H 0.081 0.919 19 C -0.007 4.007 20 N -0.920 5.920 21 Si 0.914 3.086 22 H -0.285 1.285 23 H -0.289 1.289 24 H -0.284 1.284 25 C 0.109 3.891 26 C -0.117 4.117 27 C -0.119 4.119 28 C -0.150 4.150 29 C -0.146 4.146 30 H 0.058 0.942 31 H 0.064 0.936 32 H 0.067 0.933 33 H 0.052 0.948 34 H 0.040 0.960 35 H 0.106 0.894 36 H 0.104 0.896 37 H 0.401 0.599 38 H 0.088 0.912 39 H 0.057 0.943 40 H 0.036 0.964 41 H 0.080 0.920 42 H 0.256 0.744 43 H 0.043 0.957 44 H 0.041 0.959 45 H 0.079 0.921 46 H 0.090 0.910 47 H 0.073 0.927 48 H 0.059 0.941 49 H 0.058 0.942 50 H 0.055 0.945 51 H 0.051 0.949 52 H 0.056 0.944 53 H 0.064 0.936 For the charges calculated by CM2: Dipole moment (debyes) X Y Z Total from point charges 4.452 -13.440 23.030 27.034 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.207 4.207 2 C -0.013 4.013 3 O -0.289 6.289 4 C 0.054 3.946 5 H 0.085 0.915 6 C -0.183 4.183 7 C 0.302 3.698 8 O -0.408 6.408 9 N -0.373 5.373 10 C 0.005 3.995 11 H 0.123 0.877 12 C -0.038 4.038 13 O -0.278 6.278 14 C 0.010 3.990 15 C 0.051 3.949 16 N -0.333 5.333 17 C -0.351 4.351 18 H 0.098 0.902 19 C -0.203 4.203 20 N -0.566 5.566 21 Si 0.744 3.256 22 H -0.211 1.211 23 H -0.215 1.215 24 H -0.211 1.211 25 C -0.018 4.018 26 C -0.155 4.155 27 C -0.138 4.138 28 C -0.207 4.207 29 C -0.203 4.203 30 H 0.078 0.922 31 H 0.083 0.917 32 H 0.086 0.914 33 H 0.071 0.929 34 H 0.058 0.942 35 H 0.124 0.876 36 H 0.123 0.877 37 H 0.235 0.765 38 H 0.106 0.894 39 H 0.075 0.925 40 H 0.054 0.946 41 H 0.098 0.902 42 H 0.066 0.934 43 H 0.062 0.938 44 H 0.059 0.941 45 H 0.097 0.903 46 H 0.108 0.892 47 H 0.091 0.909 48 H 0.078 0.922 49 H 0.077 0.923 50 H 0.074 0.926 51 H 0.070 0.930 52 H 0.075 0.925 53 H 0.083 0.917 Dipole moment (debyes) X Y Z Total from point charges 3.782 -13.866 22.344 26.568 hybrid contribution 0.270 1.300 -0.011 1.328 sum 4.052 -12.566 22.333 25.944 Atomic orbital electron populations 1.21720 0.94049 1.01954 1.02990 1.22402 0.95288 0.83971 0.99684 1.87898 1.20241 1.26343 1.94450 1.21732 0.82567 0.96953 0.93330 0.91490 1.21939 1.03439 1.00777 0.92176 1.20568 0.77835 0.83051 0.88360 1.90699 1.42003 1.30180 1.77903 1.45980 1.63506 1.20108 1.07714 1.22688 0.99394 0.98093 0.79368 0.87712 1.23814 0.96658 0.99448 0.83876 1.88853 1.86432 1.28064 1.24420 1.24220 0.96473 0.80999 0.97295 1.22328 0.97640 0.83207 0.91735 1.45686 1.19655 1.01657 1.66289 1.18913 1.03936 0.80878 1.31327 0.90171 1.24801 0.97201 0.97827 1.00511 1.69801 1.11882 1.31208 1.43745 0.86181 0.77308 0.84047 0.78091 1.21147 1.21508 1.21068 1.22983 0.92518 0.90698 0.95579 1.22670 1.02947 0.94743 0.95181 1.21635 0.96220 0.99638 0.96297 1.21896 0.94663 1.02087 1.02092 1.21821 1.01927 0.96615 0.99931 0.92244 0.91689 0.91416 0.92948 0.94166 0.87595 0.87746 0.76499 0.89367 0.92500 0.94633 0.90228 0.93438 0.93847 0.94129 0.90278 0.89169 0.90869 0.92179 0.92284 0.92555 0.93004 0.92476 0.91742 Number of geometries 1 Number of calculations of the screened coulomb radii 1 The total number of SCF iterations 14. 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.32 10.28 37.16 0.38 -0.94 16 2 C 0.06 0.62 6.44 37.16 0.24 0.86 16 3 O -0.37 -4.19 8.68 -49.63 -0.43 -4.62 16 4 C 0.11 1.03 1.42 -26.73 -0.04 0.99 16 5 H 0.07 0.66 7.64 -51.93 -0.40 0.27 16 6 C -0.14 -1.20 4.62 -27.88 -0.13 -1.33 16 7 C 0.52 6.20 7.71 -10.99 -0.08 6.11 16 8 O -0.53 -8.40 15.06 -8.85 -0.13 -8.53 16 9 N -0.72 -8.15 5.19 -51.74 -0.27 -8.42 16 10 C 0.11 1.58 3.51 -65.79 -0.23 1.35 16 11 H 0.10 1.70 7.54 -51.93 -0.39 1.31 16 12 C 0.04 0.60 7.42 38.24 0.28 0.88 16 13 O -0.36 -6.77 11.29 -35.23 -0.40 -7.17 16 14 C 0.05 1.00 1.56 -90.29 -0.14 0.86 16 15 C 0.18 4.30 7.00 -3.27 -0.02 4.28 16 16 N -0.60 -15.56 3.68 -181.58 -0.67 -16.23 16 17 C -0.22 -6.22 10.70 -15.06 -0.16 -6.38 16 18 H 0.08 2.14 7.83 -52.49 -0.41 1.73 16 19 C -0.01 -0.18 5.50 -17.43 -0.10 -0.28 16 20 N -0.92 -27.21 11.39 55.49 0.63 -26.58 16 21 Si 0.91 21.99 29.55 -169.99 -5.02 16.97 16 22 H -0.28 -6.77 7.11 56.52 0.40 -6.37 16 23 H -0.29 -7.02 7.11 56.52 0.40 -6.62 16 24 H -0.28 -6.68 7.11 56.52 0.40 -6.28 16 25 C 0.11 2.28 8.43 -3.26 -0.03 2.25 16 26 C -0.12 -1.94 6.20 -25.45 -0.16 -2.10 16 27 C -0.12 -0.82 2.89 -90.62 -0.26 -1.08 16 28 C -0.15 -0.83 8.96 37.16 0.33 -0.50 16 29 C -0.15 -0.98 8.45 37.15 0.31 -0.66 16 30 H 0.06 0.49 8.14 -51.93 -0.42 0.07 16 31 H 0.06 0.66 8.14 -51.93 -0.42 0.24 16 32 H 0.07 0.47 8.14 -51.93 -0.42 0.05 16 33 H 0.05 0.57 8.14 -51.93 -0.42 0.15 16 34 H 0.04 0.32 7.67 -51.93 -0.40 -0.07 16 35 H 0.11 0.62 6.36 -51.93 -0.33 0.29 16 36 H 0.10 0.80 6.79 -51.93 -0.35 0.45 16 37 H 0.40 3.48 8.41 -40.82 -0.34 3.14 16 38 H 0.09 1.25 8.14 -51.93 -0.42 0.83 16 39 H 0.06 0.75 8.14 -51.93 -0.42 0.33 16 40 H 0.04 0.92 8.14 -51.93 -0.42 0.49 16 41 H 0.08 2.09 6.89 -51.93 -0.36 1.73 16 42 H 0.26 7.34 8.31 -40.82 -0.34 7.00 16 43 H 0.04 0.83 7.85 -51.93 -0.41 0.42 16 44 H 0.04 0.86 8.14 -51.93 -0.42 0.44 16 45 H 0.08 1.20 7.37 -51.93 -0.38 0.82 16 46 H 0.09 1.64 7.75 -51.93 -0.40 1.23 16 47 H 0.07 0.52 7.68 -51.93 -0.40 0.12 16 48 H 0.06 0.30 6.47 -51.93 -0.34 -0.04 16 49 H 0.06 0.29 8.14 -51.93 -0.42 -0.13 16 50 H 0.06 0.31 8.14 -51.93 -0.42 -0.11 16 51 H 0.05 0.31 8.14 -51.93 -0.42 -0.12 16 52 H 0.06 0.34 6.83 -51.93 -0.35 -0.01 16 53 H 0.06 0.52 7.68 -51.92 -0.40 0.12 16 LS Contribution 415.89 15.07 6.27 6.27 Total: -1.00 -33.28 415.89 -9.27 -42.55 By element: Atomic # 1 Polarization: 10.91 SS G_CDS: -9.45 Total: 1.47 kcal Atomic # 6 Polarization: 4.10 SS G_CDS: 0.20 Total: 4.30 kcal Atomic # 7 Polarization: -50.93 SS G_CDS: -0.30 Total: -51.23 kcal Atomic # 8 Polarization: -19.36 SS G_CDS: -0.96 Total: -20.32 kcal Atomic # 14 Polarization: 21.99 SS G_CDS: -5.02 Total: 16.97 kcal Total LS contribution 6.27 Total: 6.27 kcal Total: -33.28 -9.27 -42.55 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. ZINC001755061285.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 -39.987 kcal (2) G-P(sol) polarization free energy of solvation -33.281 kcal (3) G-ENP(sol) elect.-nuc.-pol. free energy of system -73.267 kcal (4) G-CDS(sol) cavity-dispersion-solvent structure free energy -9.269 kcal (5) G-P-CDS(sol) = G-P(sol) + G-CDS(sol) = (2) + (4) -42.549 kcal (6) G-S(sol) free energy of system = (1) + (5) -82.536 kcal 1SCF WAS SPECIFIED, SO GEOMETRY OPTIMIZATION WAS NOT USED Total computer time = 0.82 seconds