Wall clock time and date at job start Tue Mar 31 2020 23:43:07 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 H 1.08997 * 109.47120 * 2 1 4 4 C 1.50698 * 109.47107 * 240.00029 * 2 1 3 5 5 H 1.08005 * 120.00134 * 335.15284 * 4 2 1 6 6 C 1.37879 * 120.00282 * 155.15314 * 4 2 1 7 7 N 1.31513 * 120.00041 * 5.34423 * 6 4 2 8 8 Si 1.86802 * 120.00077 * 185.33800 * 6 4 2 9 9 H 1.48503 * 109.47024 * 180.02562 * 8 6 4 10 10 H 1.48502 * 109.47210 * 299.99904 * 8 6 4 11 11 H 1.48498 * 109.47037 * 60.00424 * 8 6 4 12 12 N 1.46497 * 109.47212 * 120.00218 * 2 1 3 13 13 C 1.34780 * 119.99771 * 274.99543 * 12 2 1 14 14 O 1.21280 * 120.00140 * 0.02562 * 13 12 2 15 15 C 1.50703 * 120.00143 * 179.97438 * 13 12 2 16 16 H 1.08996 * 109.46036 * 299.97768 * 15 13 12 17 17 C 1.53038 * 109.45767 * 59.93474 * 15 13 12 18 18 C 1.53186 * 109.31407 * 181.38478 * 17 15 13 19 19 N 1.46928 * 108.77799 * 305.36951 * 18 17 15 20 20 C 1.34775 * 120.63117 * 233.63068 * 19 18 17 21 21 O 1.21509 * 120.00243 * 185.62592 * 20 19 18 22 22 O 1.34643 * 119.99638 * 5.62478 * 20 19 18 23 23 C 1.45207 * 116.99653 * 179.97438 * 22 20 19 24 24 C 1.53005 * 109.46821 * 60.00319 * 23 22 20 25 25 C 1.52996 * 109.46815 * 179.97438 * 23 22 20 26 26 C 1.53001 * 109.47178 * 300.00090 * 23 22 20 27 27 C 1.46926 * 118.73896 * 53.60494 * 19 18 17 28 28 H 1.08995 * 109.58670 * 66.18056 * 27 19 18 29 29 C 1.53002 * 109.70058 * 186.53927 * 27 19 18 30 30 C 1.53038 * 109.46257 * 180.02562 * 15 13 12 31 31 H 1.09002 * 109.47399 * 295.05133 * 1 2 3 32 32 H 1.09006 * 109.47088 * 55.04990 * 1 2 3 33 33 H 1.09004 * 109.47391 * 175.05122 * 1 2 3 34 34 H 0.97003 * 119.99606 * 355.03863 * 7 6 4 35 35 H 0.96998 * 120.00763 * 94.99825 * 12 2 1 36 36 H 1.09000 * 109.49917 * 301.34097 * 17 15 13 37 37 H 1.09003 * 109.50138 * 61.43335 * 17 15 13 38 38 H 1.08998 * 109.59192 * 185.57360 * 18 17 15 39 39 H 1.09010 * 109.70367 * 65.22270 * 18 17 15 40 40 H 1.08993 * 109.47469 * 59.99663 * 24 23 22 41 41 H 1.08999 * 109.46965 * 180.02562 * 24 23 22 42 42 H 1.09006 * 109.46900 * 299.99532 * 24 23 22 43 43 H 1.09002 * 109.46703 * 60.00778 * 25 23 22 44 44 H 1.09008 * 109.46705 * 179.97438 * 25 23 22 45 45 H 1.08993 * 109.47251 * 300.00200 * 25 23 22 46 46 H 1.09000 * 109.47005 * 60.00108 * 26 23 22 47 47 H 1.09000 * 109.46785 * 179.97438 * 26 23 22 48 48 H 1.08997 * 109.47058 * 300.00330 * 26 23 22 49 49 H 1.09009 * 109.46840 * 172.26756 * 29 27 19 50 50 H 1.08995 * 109.47431 * 292.26251 * 29 27 19 51 51 H 1.09001 * 109.47317 * 52.27142 * 29 27 19 52 52 H 1.08999 * 109.50135 * 298.57341 * 30 15 13 53 53 H 1.09009 * 109.49258 * 58.63057 * 30 15 13 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 1 1.8933 1.0276 0.0000 4 6 2.0323 -0.7104 -1.2304 5 1 1.4121 -1.4468 -1.7199 6 6 3.2836 -0.4203 -1.7315 7 7 4.0802 0.3971 -1.0782 8 14 3.8476 -1.1885 -3.3381 9 1 5.2205 -0.7200 -3.6557 10 1 3.8480 -2.6677 -3.2063 11 1 2.9229 -0.7890 -4.4292 12 7 2.0183 -0.6906 1.1961 13 6 2.1471 -0.0182 2.3571 14 8 1.8583 1.1584 2.4117 15 6 2.6489 -0.7288 3.5877 16 1 1.9756 -1.5502 3.8327 17 6 4.0510 -1.2813 3.3213 18 6 4.5438 -2.0324 4.5621 19 7 4.4591 -1.1353 5.7227 20 6 5.5365 -0.9183 6.5029 21 8 5.4292 -0.2570 7.5166 22 8 6.7316 -1.4344 6.1590 23 6 7.8480 -1.1597 7.0459 24 6 7.5436 -1.7190 8.4371 25 6 9.1118 -1.8233 6.4954 26 6 8.0646 0.3520 7.1401 27 6 3.1778 -0.4789 6.0162 28 1 2.4407 -1.2285 6.3039 29 6 3.3616 0.5241 7.1569 30 6 2.6999 0.2547 4.7592 31 1 -0.3634 0.4351 0.9310 32 1 -0.3633 0.5887 -0.8424 33 1 -0.3634 -1.0238 -0.0887 34 1 3.7618 0.8519 -0.2827 35 1 2.2490 -1.6318 1.1526 36 1 4.0178 -1.9638 2.4721 37 1 4.7306 -0.4581 3.1006 38 1 5.5782 -2.3426 4.4147 39 1 3.9200 -2.9100 4.7324 40 1 7.3892 -2.7959 8.3701 41 1 8.3814 -1.5124 9.1031 42 1 6.6430 -1.2466 8.8296 43 1 9.3289 -1.4251 5.5041 44 1 9.9500 -1.6167 7.1610 45 1 8.9572 -2.9001 6.4286 46 1 7.1641 0.8243 7.5325 47 1 8.9024 0.5581 7.8062 48 1 8.2814 0.7504 6.1489 49 1 2.4478 1.1055 7.2801 50 1 3.5784 -0.0120 8.0808 51 1 4.1892 1.1935 6.9223 52 1 3.3918 1.0637 4.5252 53 1 1.7052 0.6655 4.9328 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 ZINC000492660545.mol2 53 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 23:43:07 Heat of formation + Delta-G solvation = -130.028362 kcal Electronic energy + Delta-G solvation = -31905.052856 eV Core-core repulsion = 27776.452145 eV Total energy + Delta-G solvation = -4128.600712 eV No. of doubly occupied orbitals = 66 Molecular weight (most abundant/longest-lived isotopes) = 340.220 amu Computer time = 0.99 seconds Orbital eigenvalues (eV) -41.14005 -39.53856 -38.21295 -36.46000 -34.95386 -34.35766 -32.95601 -32.35156 -31.95153 -29.38001 -27.62690 -27.54499 -27.53222 -25.16539 -24.66906 -23.17718 -22.36912 -21.17059 -20.44154 -19.38609 -19.12076 -17.64929 -17.37643 -16.70509 -16.32905 -16.04966 -15.52525 -14.94374 -14.77259 -14.65015 -14.25906 -13.91847 -13.88924 -13.83627 -13.63040 -13.39983 -13.32761 -12.71808 -12.68392 -12.46182 -12.42916 -12.31262 -12.05131 -11.89301 -11.85657 -11.77108 -11.62781 -11.53341 -11.22020 -11.14854 -11.06662 -11.04325 -10.84791 -10.62768 -10.54263 -10.41909 -10.24567 -10.04705 -9.38584 -9.10898 -8.75778 -8.69653 -8.53800 -7.81097 -6.56078 -4.36325 2.08327 2.81410 3.32546 3.35538 3.47908 3.54473 3.71033 4.15727 4.37933 4.55094 4.59700 4.71421 4.74158 4.82368 4.97491 5.01917 5.19721 5.22768 5.27812 5.35407 5.36682 5.40715 5.45468 5.50132 5.56854 5.59787 5.65184 5.68966 5.77953 5.81012 5.85397 5.93379 6.04178 6.10600 6.14185 6.17146 6.26158 6.32841 6.34799 6.42491 6.45896 6.51223 6.59282 6.63270 7.03118 7.23726 7.38549 7.67186 7.89680 8.24579 8.33258 8.59970 8.79804 9.23173 9.43963 10.81527 Molecular weight = 340.22amu Principal moments of inertia in cm(-1) A = 0.013619 B = 0.002994 C = 0.002593 Principal moments of inertia in units of 10**(-40)*gram-cm**2 A = 2055.397641 B = 9350.300810 C =10794.928077 Net atomic charges, atomic populations, and dipole contributions calculated with CM2 Atom NO. Type Charge No. of electrons 1 C -0.140 4.140 2 C 0.237 3.763 3 H 0.050 0.950 4 C -0.609 4.609 5 H 0.055 0.945 6 C 0.031 3.969 7 N -0.888 5.888 8 Si 1.069 2.931 9 H -0.282 1.282 10 H -0.286 1.286 11 H -0.286 1.286 12 N -0.702 5.702 13 C 0.501 3.499 14 O -0.565 6.565 15 C -0.110 4.110 16 H 0.089 0.911 17 C -0.125 4.125 18 C 0.115 3.885 19 N -0.591 5.591 20 C 0.654 3.346 21 O -0.565 6.565 22 O -0.365 6.365 23 C 0.136 3.864 24 C -0.181 4.181 25 C -0.141 4.141 26 C -0.183 4.183 27 C 0.159 3.841 28 H 0.075 0.925 29 C -0.137 4.137 30 C -0.108 4.108 31 H 0.031 0.969 32 H 0.050 0.950 33 H 0.039 0.961 34 H 0.277 0.723 35 H 0.386 0.614 36 H 0.085 0.915 37 H 0.074 0.926 38 H 0.094 0.906 39 H 0.071 0.929 40 H 0.057 0.943 41 H 0.067 0.933 42 H 0.088 0.912 43 H 0.067 0.933 44 H 0.078 0.922 45 H 0.065 0.935 46 H 0.088 0.912 47 H 0.067 0.933 48 H 0.060 0.940 49 H 0.048 0.952 50 H 0.076 0.924 51 H 0.068 0.932 52 H 0.074 0.926 53 H 0.083 0.917 For the charges calculated by CM2: Dipole moment (debyes) X Y Z Total from point charges 5.856 -7.273 17.117 19.498 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.198 4.198 2 C 0.137 3.863 3 H 0.068 0.932 4 C -0.648 4.648 5 H 0.074 0.926 6 C -0.178 4.178 7 N -0.522 5.522 8 Si 0.897 3.103 9 H -0.207 1.207 10 H -0.213 1.213 11 H -0.212 1.212 12 N -0.351 5.351 13 C 0.286 3.714 14 O -0.445 6.445 15 C -0.132 4.132 16 H 0.108 0.892 17 C -0.163 4.163 18 C -0.008 4.008 19 N -0.332 5.332 20 C 0.409 3.591 21 O -0.455 6.455 22 O -0.280 6.280 23 C 0.100 3.900 24 C -0.239 4.239 25 C -0.198 4.198 26 C -0.241 4.241 27 C 0.057 3.943 28 H 0.093 0.907 29 C -0.195 4.195 30 C -0.147 4.147 31 H 0.050 0.950 32 H 0.069 0.931 33 H 0.058 0.942 34 H 0.089 0.911 35 H 0.219 0.781 36 H 0.104 0.896 37 H 0.093 0.907 38 H 0.112 0.888 39 H 0.090 0.910 40 H 0.076 0.924 41 H 0.086 0.914 42 H 0.106 0.894 43 H 0.086 0.914 44 H 0.097 0.903 45 H 0.084 0.916 46 H 0.107 0.893 47 H 0.086 0.914 48 H 0.079 0.921 49 H 0.067 0.933 50 H 0.095 0.905 51 H 0.087 0.913 52 H 0.093 0.907 53 H 0.102 0.898 Dipole moment (debyes) X Y Z Total from point charges 6.977 -6.335 16.422 18.934 hybrid contribution -1.168 -0.423 0.881 1.523 sum 5.809 -6.758 17.303 19.463 Atomic orbital electron populations 1.21996 0.98632 0.99851 0.99337 1.18852 0.89844 0.94898 0.82707 0.93162 1.22110 1.04901 1.26964 1.10783 0.92645 1.26004 0.93199 0.96017 1.02629 1.70749 1.36770 1.25698 1.18965 0.83547 0.76387 0.75682 0.74734 1.20703 1.21281 1.21185 1.47627 1.64113 1.13547 1.09769 1.20521 0.78381 0.87559 0.84978 1.90612 1.52228 1.16246 1.85438 1.21378 0.97667 0.99142 0.95034 0.89239 1.21895 0.95711 1.01209 0.97501 1.21763 1.03066 0.91713 0.84249 1.47135 1.08543 1.48401 1.29078 1.18062 0.80627 0.78241 0.82137 1.90924 1.84623 1.43548 1.26411 1.86320 1.16885 1.71890 1.52897 1.22017 0.83294 0.95934 0.88801 1.22501 1.04324 1.01356 0.95676 1.21860 0.94462 1.01297 1.02227 1.22514 1.04204 0.94597 1.02740 1.21339 0.82188 0.94758 0.95973 0.90673 1.22199 1.00858 0.98645 0.97752 1.21817 1.02987 0.97749 0.92128 0.95011 0.93137 0.94174 0.91122 0.78112 0.89623 0.90691 0.88805 0.91033 0.92362 0.91393 0.89356 0.91361 0.90319 0.91551 0.89334 0.91384 0.92114 0.93316 0.90546 0.91332 0.90696 0.89819 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.14 -3.03 9.46 37.15 0.35 -2.68 16 2 C 0.24 5.81 2.13 -69.09 -0.15 5.66 16 3 H 0.05 1.43 6.37 -51.93 -0.33 1.10 16 4 C -0.61 -15.82 8.50 -34.44 -0.29 -16.11 16 5 H 0.06 1.42 7.81 -52.48 -0.41 1.01 16 6 C 0.03 0.81 5.45 -17.82 -0.10 0.72 16 7 N -0.89 -25.53 13.63 55.43 0.76 -24.78 16 8 Si 1.07 23.46 29.91 -169.99 -5.08 18.38 16 9 H -0.28 -6.29 7.11 56.52 0.40 -5.89 16 10 H -0.29 -5.94 7.11 56.52 0.40 -5.54 16 11 H -0.29 -6.06 7.11 56.52 0.40 -5.66 16 12 N -0.70 -14.88 4.95 -54.94 -0.27 -15.15 16 13 C 0.50 10.18 7.12 -10.98 -0.08 10.10 16 14 O -0.56 -13.66 15.64 5.56 0.09 -13.57 16 15 C -0.11 -1.66 2.49 -91.73 -0.23 -1.89 16 16 H 0.09 1.14 8.14 -51.93 -0.42 0.72 16 17 C -0.12 -1.75 5.27 -26.62 -0.14 -1.89 16 18 C 0.11 1.33 6.46 -3.71 -0.02 1.31 16 19 N -0.59 -7.59 2.97 -170.97 -0.51 -8.10 16 20 C 0.65 8.95 7.53 54.05 0.41 9.36 16 21 O -0.56 -8.13 7.46 -19.28 -0.14 -8.27 16 22 O -0.36 -4.66 8.53 -40.33 -0.34 -5.01 16 23 C 0.14 1.43 1.13 -90.62 -0.10 1.32 16 24 C -0.18 -1.74 8.37 37.16 0.31 -1.43 16 25 C -0.14 -1.15 8.85 37.16 0.33 -0.82 16 26 C -0.18 -1.94 8.37 37.16 0.31 -1.63 16 27 C 0.16 1.97 3.48 -67.60 -0.24 1.74 16 28 H 0.08 0.82 8.14 -51.93 -0.42 0.40 16 29 C -0.14 -1.72 7.62 37.16 0.28 -1.44 16 30 C -0.11 -1.57 3.94 -26.61 -0.10 -1.68 16 31 H 0.03 0.67 8.14 -51.93 -0.42 0.25 16 32 H 0.05 1.07 8.14 -51.93 -0.42 0.64 16 33 H 0.04 0.81 8.14 -51.93 -0.42 0.38 16 34 H 0.28 8.23 6.77 -40.82 -0.28 7.96 16 35 H 0.39 7.45 7.99 -40.82 -0.33 7.12 16 36 H 0.09 1.19 7.52 -51.93 -0.39 0.80 16 37 H 0.07 1.24 8.14 -51.93 -0.42 0.82 16 38 H 0.09 1.06 6.92 -51.93 -0.36 0.70 16 39 H 0.07 0.69 8.14 -51.92 -0.42 0.27 16 40 H 0.06 0.53 8.14 -51.93 -0.42 0.11 16 41 H 0.07 0.52 8.14 -51.93 -0.42 0.10 16 42 H 0.09 1.03 5.88 -51.93 -0.31 0.72 16 43 H 0.07 0.58 8.14 -51.93 -0.42 0.16 16 44 H 0.08 0.49 8.14 -51.92 -0.42 0.07 16 45 H 0.07 0.51 8.14 -51.93 -0.42 0.09 16 46 H 0.09 1.12 5.88 -51.93 -0.31 0.81 16 47 H 0.07 0.58 8.14 -51.93 -0.42 0.15 16 48 H 0.06 0.67 8.14 -51.93 -0.42 0.24 16 49 H 0.05 0.55 8.14 -51.92 -0.42 0.13 16 50 H 0.08 0.99 6.70 -51.93 -0.35 0.64 16 51 H 0.07 0.98 6.81 -51.93 -0.35 0.63 16 52 H 0.07 1.27 8.13 -51.93 -0.42 0.84 16 53 H 0.08 1.22 8.07 -51.92 -0.42 0.80 16 LS Contribution 407.58 15.07 6.14 6.14 Total: -1.00 -30.92 407.58 -8.21 -39.13 By element: Atomic # 1 Polarization: 19.96 SS G_CDS: -9.38 Total: 10.58 kcal Atomic # 6 Polarization: 0.09 SS G_CDS: 0.54 Total: 0.64 kcal Atomic # 7 Polarization: -48.00 SS G_CDS: -0.02 Total: -48.03 kcal Atomic # 8 Polarization: -26.44 SS G_CDS: -0.40 Total: -26.84 kcal Atomic # 14 Polarization: 23.46 SS G_CDS: -5.08 Total: 18.38 kcal Total LS contribution 6.14 Total: 6.14 kcal Total: -30.92 -8.21 -39.13 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. ZINC000492660545.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 -90.897 kcal (2) G-P(sol) polarization free energy of solvation -30.923 kcal (3) G-ENP(sol) elect.-nuc.-pol. free energy of system -121.820 kcal (4) G-CDS(sol) cavity-dispersion-solvent structure free energy -8.208 kcal (5) G-P-CDS(sol) = G-P(sol) + G-CDS(sol) = (2) + (4) -39.131 kcal (6) G-S(sol) free energy of system = (1) + (5) -130.028 kcal 1SCF WAS SPECIFIED, SO GEOMETRY OPTIMIZATION WAS NOT USED Total computer time = 0.99 seconds