Wall clock time and date at job start Tue Mar 31 2020 21:18:21 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.30999 * 1 3 3 C 1.50701 * 119.99474 * 2 1 4 4 C 1.53007 * 109.46998 * 128.00596 * 3 2 1 5 5 C 1.50697 * 109.54223 * 293.80657 * 4 3 2 6 6 O 1.21278 * 120.00508 * 250.95108 * 5 4 3 7 7 N 1.34784 * 119.99553 * 70.95406 * 5 4 3 8 8 C 1.46502 * 119.99635 * 174.87848 * 7 5 4 9 9 H 1.08998 * 109.52559 * 329.43739 * 8 7 5 10 10 C 1.53046 * 109.49267 * 89.51753 * 8 7 5 11 11 C 1.53041 * 109.52844 * 58.57820 * 10 8 7 12 12 C 1.53194 * 109.31171 * 61.36893 * 11 10 8 13 13 N 1.46920 * 108.77548 * 305.36533 * 12 11 10 14 14 C 1.36941 * 120.62885 * 233.62975 * 13 12 11 15 15 H 1.08003 * 120.00159 * 143.21722 * 14 13 12 16 16 C 1.38817 * 119.99969 * 323.21525 * 14 13 12 17 17 N 1.31614 * 120.00088 * 359.79064 * 16 14 13 18 18 Si 1.86802 * 119.99860 * 179.78357 * 16 14 13 19 19 H 1.48506 * 109.46664 * 60.00357 * 18 16 14 20 20 H 1.48492 * 109.46972 * 180.02562 * 18 16 14 21 21 H 1.48500 * 109.47246 * 300.00725 * 18 16 14 22 22 C 1.46919 * 118.74222 * 53.60978 * 13 12 11 23 23 H 1.09001 * 109.58957 * 66.17986 * 22 13 12 24 24 C 1.53001 * 109.59028 * 186.45943 * 22 13 12 25 25 C 1.53174 * 109.54168 * 53.98878 * 4 3 2 26 26 C 1.53288 * 109.06085 * 176.85386 * 25 4 3 27 27 O 1.42903 * 109.45780 * 302.44448 * 26 25 4 28 28 C 1.42902 * 114.10168 * 61.15727 * 27 26 25 29 29 C 1.53094 * 109.41376 * 298.84283 * 28 27 26 30 30 H 1.08000 * 120.00440 * 0.02562 * 1 2 3 31 31 H 1.08003 * 119.99245 * 180.02562 * 1 2 3 32 32 H 1.08000 * 120.00440 * 179.97438 * 2 1 3 33 33 H 1.08997 * 109.47500 * 8.00411 * 3 2 1 34 34 H 1.09000 * 109.47280 * 248.00037 * 3 2 1 35 35 H 0.97000 * 119.99930 * 354.88571 * 7 5 4 36 36 H 1.08999 * 109.45869 * 298.55704 * 10 8 7 37 37 H 1.09006 * 109.45249 * 178.59100 * 10 8 7 38 38 H 1.08998 * 109.49630 * 181.31847 * 11 10 8 39 39 H 1.08994 * 109.49561 * 301.41927 * 11 10 8 40 40 H 1.09001 * 109.58563 * 65.15072 * 12 11 10 41 41 H 1.09005 * 109.58359 * 185.57519 * 12 11 10 42 42 H 0.96998 * 120.00661 * 179.97438 * 17 16 14 43 43 H 1.09003 * 109.46611 * 60.00013 * 24 22 13 44 44 H 1.08996 * 109.47291 * 179.97438 * 24 22 13 45 45 H 1.09001 * 109.47114 * 300.00375 * 24 22 13 46 46 H 1.09003 * 109.52243 * 296.83940 * 25 4 3 47 47 H 1.08995 * 109.52350 * 57.03313 * 25 4 3 48 48 H 1.08997 * 109.45401 * 182.45000 * 26 25 4 49 49 H 1.08999 * 109.46751 * 62.44664 * 26 25 4 50 50 H 1.08997 * 109.47979 * 178.85693 * 28 27 26 51 51 H 1.09001 * 109.48390 * 58.82697 * 28 27 26 52 52 H 1.08999 * 109.52024 * 297.69631 * 29 28 27 53 53 H 1.08996 * 109.51747 * 177.51040 * 29 28 27 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.3100 0.0000 0.0000 3 6 2.0634 1.3052 0.0000 4 6 3.0876 1.3028 1.1367 5 6 2.3701 1.2869 2.4618 6 8 2.3759 2.2740 3.1665 7 7 1.7218 0.1757 2.8637 8 6 1.1201 0.1212 4.1983 9 1 1.6982 0.7403 4.8842 10 6 -0.3185 0.6392 4.1314 11 6 -1.1224 -0.2093 3.1435 12 6 -1.1535 -1.6612 3.6311 13 7 0.2268 -2.1205 3.8372 14 6 0.6776 -3.2664 3.2380 15 1 1.6927 -3.3216 2.8733 16 6 -0.1701 -4.3569 3.0994 17 7 -1.4058 -4.2910 3.5476 18 14 0.4431 -5.9178 2.2766 19 1 1.5995 -6.4566 3.0368 20 1 -0.6471 -6.9258 2.2526 21 1 0.8643 -5.6117 0.8859 22 6 1.1143 -1.3269 4.6981 23 1 0.7506 -1.3578 5.7252 24 6 2.5322 -1.8987 4.6388 25 6 3.9783 0.0620 1.0212 26 6 5.0492 0.1107 2.1169 27 8 5.8052 1.3167 1.9902 28 6 5.0309 2.5089 2.1358 29 6 3.9610 2.5576 1.0419 30 1 -0.5401 0.9353 -0.0004 31 1 -0.5399 -0.9354 0.0004 32 1 1.8501 -0.9353 0.0004 33 1 1.3637 2.1286 0.1431 34 1 2.5784 1.4274 -0.9528 35 1 1.6539 -0.5879 2.2694 36 1 -0.3155 1.6776 3.7998 37 1 -0.7738 0.5750 5.1198 38 1 -2.1405 0.1746 3.0787 39 1 -0.6532 -0.1655 2.1607 40 1 -1.7022 -1.7193 4.5712 41 1 -1.6392 -2.2876 2.8828 42 1 -1.9983 -5.0529 3.4505 43 1 2.8959 -1.8623 3.6119 44 1 3.1890 -1.3084 5.2776 45 1 2.5219 -2.9325 4.9841 46 1 3.3708 -0.8350 1.1423 47 1 4.4588 0.0483 0.0430 48 1 5.7134 -0.7473 2.0137 49 1 4.5694 0.0835 3.0952 50 1 5.6821 3.3784 2.0464 51 1 4.5508 2.5134 3.1144 52 1 4.4413 2.5943 0.0641 53 1 3.3413 3.4439 1.1778 RHF calculation, no. of doubly occupied orbitals= 65 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 ZINC001486186398.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 21:18:21 Heat of formation + Delta-G solvation = -54.140234 kcal Electronic energy + Delta-G solvation = -33391.069103 eV Core-core repulsion = 29456.215199 eV Total energy + Delta-G solvation = -3934.853904 eV No. of doubly occupied orbitals = 65 Molecular weight (most abundant/longest-lived isotopes) = 336.220 amu Computer time = 1.13 seconds Orbital eigenvalues (eV) -40.94381 -38.66264 -37.29780 -35.91660 -33.95446 -32.74589 -32.34024 -31.67401 -30.92377 -29.34394 -27.66646 -26.32224 -25.04390 -23.83303 -23.67805 -22.39433 -21.51647 -20.74208 -19.77772 -19.12685 -18.59728 -16.74576 -16.47321 -16.12564 -15.93436 -15.45131 -15.04666 -15.01196 -14.65927 -14.47774 -14.12088 -13.86929 -13.75548 -13.59173 -13.08805 -13.05275 -12.88983 -12.68457 -12.43428 -12.24006 -12.14666 -11.81653 -11.51578 -11.42010 -11.29489 -11.22581 -11.14957 -10.95700 -10.75917 -10.66257 -10.60425 -10.42932 -10.32600 -10.15253 -10.02175 -9.70646 -9.53322 -9.35674 -9.28354 -8.96552 -8.50398 -8.42392 -6.58343 -5.85973 -3.40542 2.41603 2.98773 3.29498 3.35859 3.37415 3.62701 3.92852 4.39091 4.47319 4.70355 4.76135 4.98022 5.09514 5.21152 5.24712 5.27730 5.31738 5.34892 5.49472 5.56917 5.61016 5.61612 5.65885 5.69380 5.72996 5.77162 5.83342 5.90847 5.98211 6.00889 6.04492 6.11428 6.16120 6.19491 6.24067 6.34333 6.42513 6.54414 6.69203 6.73592 6.80550 6.81008 6.91743 6.95710 6.96723 7.14683 7.23177 7.41220 7.69892 7.73126 8.07795 8.43798 8.66689 9.16877 9.75516 10.37175 11.20613 Molecular weight = 336.22amu Principal moments of inertia in cm(-1) A = 0.013133 B = 0.004847 C = 0.004594 Principal moments of inertia in units of 10**(-40)*gram-cm**2 A = 2131.475726 B = 5775.259309 C = 6093.746112 Net atomic charges, atomic populations, and dipole contributions calculated with CM2 Atom NO. Type Charge No. of electrons 1 C -0.163 4.163 2 C -0.180 4.180 3 C -0.074 4.074 4 C -0.075 4.075 5 C 0.520 3.480 6 O -0.556 6.556 7 N -0.693 5.693 8 C 0.131 3.869 9 H 0.074 0.926 10 C -0.119 4.119 11 C -0.151 4.151 12 C 0.173 3.827 13 N -0.584 5.584 14 C -0.256 4.256 15 H 0.064 0.936 16 C 0.015 3.985 17 N -0.899 5.899 18 Si 0.897 3.103 19 H -0.282 1.282 20 H -0.291 1.291 21 H -0.283 1.283 22 C 0.159 3.841 23 H 0.037 0.963 24 C -0.134 4.134 25 C -0.146 4.146 26 C 0.058 3.942 27 O -0.386 6.386 28 C 0.054 3.946 29 C -0.128 4.128 30 H 0.091 0.909 31 H 0.106 0.894 32 H 0.115 0.885 33 H 0.070 0.930 34 H 0.080 0.920 35 H 0.409 0.591 36 H 0.047 0.953 37 H 0.059 0.941 38 H 0.043 0.957 39 H 0.066 0.934 40 H 0.014 0.986 41 H 0.089 0.911 42 H 0.255 0.745 43 H 0.051 0.949 44 H 0.037 0.963 45 H 0.064 0.936 46 H 0.084 0.916 47 H 0.081 0.919 48 H 0.095 0.905 49 H 0.053 0.947 50 H 0.089 0.911 51 H 0.068 0.932 52 H 0.074 0.926 53 H 0.081 0.919 For the charges calculated by CM2: Dipole moment (debyes) X Y Z Total from point charges 7.459 10.142 -5.797 13.860 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.199 4.199 2 C -0.198 4.198 3 C -0.112 4.112 4 C -0.078 4.078 5 C 0.305 3.695 6 O -0.436 6.436 7 N -0.344 5.344 8 C 0.025 3.975 9 H 0.092 0.908 10 C -0.157 4.157 11 C -0.190 4.190 12 C 0.048 3.952 13 N -0.308 5.308 14 C -0.386 4.386 15 H 0.081 0.919 16 C -0.185 4.185 17 N -0.542 5.542 18 Si 0.727 3.273 19 H -0.209 1.209 20 H -0.217 1.217 21 H -0.210 1.210 22 C 0.052 3.948 23 H 0.055 0.945 24 C -0.191 4.191 25 C -0.184 4.184 26 C -0.019 4.019 27 O -0.304 6.304 28 C -0.023 4.023 29 C -0.166 4.166 30 H 0.109 0.891 31 H 0.124 0.876 32 H 0.132 0.868 33 H 0.088 0.912 34 H 0.098 0.902 35 H 0.247 0.753 36 H 0.066 0.934 37 H 0.078 0.922 38 H 0.061 0.939 39 H 0.085 0.915 40 H 0.032 0.968 41 H 0.107 0.893 42 H 0.065 0.935 43 H 0.070 0.930 44 H 0.056 0.944 45 H 0.083 0.917 46 H 0.102 0.898 47 H 0.099 0.901 48 H 0.114 0.886 49 H 0.071 0.929 50 H 0.108 0.892 51 H 0.086 0.914 52 H 0.092 0.908 53 H 0.099 0.901 Dipole moment (debyes) X Y Z Total from point charges 7.514 10.652 -5.304 14.074 hybrid contribution -0.337 -1.295 -0.572 1.455 sum 7.177 9.358 -5.876 13.176 Atomic orbital electron populations 1.23886 0.95925 1.02088 0.98028 1.22893 0.96228 0.98008 1.02624 1.20711 0.96891 0.96422 0.97170 1.20876 0.95595 0.97324 0.94004 1.20292 0.78783 0.83219 0.87182 1.90613 1.62850 1.33788 1.56336 1.46102 1.53892 1.16916 1.17526 1.21547 0.93203 0.99081 0.83681 0.90833 1.21740 0.95627 0.98448 0.99930 1.22356 0.97798 0.99370 0.99475 1.20829 0.85914 0.89445 0.99039 1.44026 1.08667 1.25267 1.52835 1.18962 0.97924 0.89578 1.32149 0.91889 1.25324 0.96614 0.97299 0.99230 1.69913 1.05434 1.30265 1.48574 0.86480 0.77965 0.83666 0.79231 1.20856 1.21668 1.20971 1.20745 0.91714 0.88642 0.93745 0.94510 1.21771 0.95469 1.01178 1.00726 1.22384 0.97593 0.97412 1.01005 1.22639 0.94060 0.87593 0.97590 1.87881 1.35333 1.13866 1.93364 1.22701 0.93321 0.87331 0.98900 1.22010 0.97340 0.97175 1.00040 0.89070 0.87599 0.86752 0.91167 0.90195 0.75251 0.93414 0.92245 0.93870 0.91539 0.96753 0.89278 0.93489 0.93013 0.94404 0.91720 0.89801 0.90058 0.88650 0.92900 0.89249 0.91407 0.90783 0.90083 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.16 -2.16 11.67 29.01 0.34 -1.82 16 2 C -0.18 -2.26 6.37 -36.03 -0.23 -2.49 16 3 C -0.07 -0.86 4.72 -27.88 -0.13 -0.99 16 4 C -0.07 -0.94 0.50 -155.48 -0.08 -1.02 16 5 C 0.52 8.09 4.87 -10.98 -0.05 8.04 16 6 O -0.56 -10.14 13.97 5.56 0.08 -10.06 16 7 N -0.69 -11.10 2.79 -53.68 -0.15 -11.25 16 8 C 0.13 2.30 2.36 -67.81 -0.16 2.14 16 9 H 0.07 1.32 7.52 -51.93 -0.39 0.93 16 10 C -0.12 -2.09 5.61 -26.69 -0.15 -2.24 16 11 C -0.15 -2.77 5.93 -26.61 -0.16 -2.93 16 12 C 0.17 3.82 5.22 -3.71 -0.02 3.80 16 13 N -0.58 -12.86 2.62 -164.78 -0.43 -13.29 16 14 C -0.26 -6.15 8.37 -15.06 -0.13 -6.27 16 15 H 0.06 1.47 6.44 -52.48 -0.34 1.13 16 16 C 0.02 0.38 5.24 -17.43 -0.09 0.29 16 17 N -0.90 -24.23 10.91 55.49 0.61 -23.62 16 18 Si 0.90 20.18 29.55 -169.99 -5.02 15.16 16 19 H -0.28 -6.36 7.11 56.52 0.40 -5.96 16 20 H -0.29 -6.74 7.11 56.52 0.40 -6.34 16 21 H -0.28 -6.26 7.11 56.52 0.40 -5.86 16 22 C 0.16 3.06 3.21 -67.60 -0.22 2.85 16 23 H 0.04 0.76 8.14 -51.93 -0.42 0.33 16 24 C -0.13 -2.36 7.81 37.16 0.29 -2.07 16 25 C -0.15 -1.73 4.30 -26.49 -0.11 -1.84 16 26 C 0.06 0.74 6.62 37.31 0.25 0.98 16 27 O -0.39 -5.28 10.77 -35.23 -0.38 -5.66 16 28 C 0.05 0.69 6.56 37.21 0.24 0.93 16 29 C -0.13 -1.51 4.56 -26.59 -0.12 -1.63 16 30 H 0.09 1.15 7.85 -52.49 -0.41 0.74 16 31 H 0.11 1.53 8.06 -52.48 -0.42 1.10 16 32 H 0.11 1.46 6.62 -52.49 -0.35 1.11 16 33 H 0.07 0.84 7.80 -51.93 -0.41 0.44 16 34 H 0.08 0.78 8.14 -51.93 -0.42 0.36 16 35 H 0.41 6.75 4.96 -40.82 -0.20 6.55 16 36 H 0.05 0.84 8.14 -51.93 -0.42 0.42 16 37 H 0.06 1.03 8.14 -51.93 -0.42 0.61 16 38 H 0.04 0.75 8.14 -51.93 -0.42 0.33 16 39 H 0.07 1.18 4.55 -51.93 -0.24 0.94 16 40 H 0.01 0.34 8.14 -51.93 -0.42 -0.08 16 41 H 0.09 2.22 6.19 -51.93 -0.32 1.89 16 42 H 0.25 6.78 8.31 -40.82 -0.34 6.44 16 43 H 0.05 0.89 6.49 -51.93 -0.34 0.55 16 44 H 0.04 0.59 8.14 -51.93 -0.42 0.17 16 45 H 0.06 1.23 7.68 -51.93 -0.40 0.83 16 46 H 0.08 1.03 5.45 -51.93 -0.28 0.75 16 47 H 0.08 0.86 8.14 -51.93 -0.42 0.43 16 48 H 0.10 1.06 8.14 -51.93 -0.42 0.64 16 49 H 0.05 0.75 7.38 -51.93 -0.38 0.37 16 50 H 0.09 0.97 8.14 -51.93 -0.42 0.54 16 51 H 0.07 1.01 6.22 -51.93 -0.32 0.69 16 52 H 0.07 0.76 8.14 -51.93 -0.42 0.34 16 53 H 0.08 0.98 7.88 -51.93 -0.41 0.57 16 LS Contribution 384.83 15.07 5.80 5.80 Total: -1.00 -27.20 384.83 -9.03 -36.22 By element: Atomic # 1 Polarization: 19.96 SS G_CDS: -9.00 Total: 10.96 kcal Atomic # 6 Polarization: -3.74 SS G_CDS: -0.53 Total: -4.27 kcal Atomic # 7 Polarization: -48.18 SS G_CDS: 0.02 Total: -48.16 kcal Atomic # 8 Polarization: -15.42 SS G_CDS: -0.30 Total: -15.72 kcal Atomic # 14 Polarization: 20.18 SS G_CDS: -5.02 Total: 15.16 kcal Total LS contribution 5.80 Total: 5.80 kcal Total: -27.20 -9.03 -36.22 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. ZINC001486186398.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 -17.917 kcal (2) G-P(sol) polarization free energy of solvation -27.196 kcal (3) G-ENP(sol) elect.-nuc.-pol. free energy of system -45.113 kcal (4) G-CDS(sol) cavity-dispersion-solvent structure free energy -9.027 kcal (5) G-P-CDS(sol) = G-P(sol) + G-CDS(sol) = (2) + (4) -36.223 kcal (6) G-S(sol) free energy of system = (1) + (5) -54.140 kcal 1SCF WAS SPECIFIED, SO GEOMETRY OPTIMIZATION WAS NOT USED Total computer time = 1.13 seconds