Wall clock time and date at job start Sat Apr 11 2020 06:09:33 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 N 2 2 C 1.31627 * 1 3 3 Si 1.86798 * 120.00083 * 2 1 4 4 H 1.48505 * 109.47267 * 269.99971 * 3 2 1 5 5 H 1.48504 * 109.47218 * 29.99841 * 3 2 1 6 6 H 1.48501 * 109.47291 * 149.99951 * 3 2 1 7 7 C 1.38803 * 119.99971 * 179.97438 * 2 1 3 8 8 H 1.08001 * 119.99846 * 180.02562 * 7 2 1 9 9 S 1.76197 * 120.00004 * 0.02562 * 7 2 1 10 10 C 1.80993 * 100.00056 * 180.02562 * 9 7 2 11 11 C 1.50703 * 109.47297 * 180.02562 * 10 9 7 12 12 O 1.21277 * 119.99568 * 359.97438 * 11 10 9 13 13 N 1.34778 * 120.00339 * 180.02562 * 11 10 9 14 14 C 1.46929 * 120.62869 * 359.97438 * 13 11 10 15 15 C 1.53191 * 108.77608 * 233.61236 * 14 13 11 16 16 C 1.53039 * 109.31134 * 305.36437 * 15 14 13 17 17 C 1.53040 * 109.53686 * 61.36605 * 16 15 14 18 18 H 1.09000 * 109.50167 * 58.62028 * 17 16 15 19 19 N 1.46905 * 109.49673 * 178.68309 * 17 16 15 20 20 C 1.46901 * 110.99870 * 204.93632 * 19 17 16 21 21 C 1.53004 * 109.46916 * 180.02562 * 20 19 17 22 22 C 1.52996 * 109.46700 * 179.97438 * 21 20 19 23 23 C 1.53000 * 109.46897 * 179.97438 * 22 21 20 24 24 C 1.53780 * 113.61835 * 272.50098 * 23 22 21 25 25 C 1.53785 * 87.08217 * 220.02024 * 24 23 22 26 26 C 1.53774 * 87.08081 * 25.42807 * 25 24 23 27 27 C 1.46928 * 120.60491 * 180.02562 * 13 11 10 28 28 H 0.97000 * 119.99928 * 0.02562 * 1 2 3 29 29 H 1.08998 * 109.47230 * 60.00044 * 10 9 7 30 30 H 1.08999 * 109.47071 * 299.99995 * 10 9 7 31 31 H 1.09004 * 109.58558 * 113.82409 * 14 13 11 32 32 H 1.09007 * 109.58285 * 353.39912 * 14 13 11 33 33 H 1.08998 * 109.49522 * 185.41254 * 15 14 13 34 34 H 1.09001 * 109.52233 * 65.33393 * 15 14 13 35 35 H 1.09002 * 109.46152 * 301.34280 * 16 15 14 36 36 H 1.08999 * 109.45932 * 181.38328 * 16 15 14 37 37 H 1.00901 * 110.99869 * 80.98343 * 19 17 16 38 38 H 1.09004 * 109.47031 * 299.99620 * 20 19 17 39 39 H 1.08995 * 109.47111 * 59.99900 * 20 19 17 40 40 H 1.08995 * 109.46615 * 299.99641 * 21 20 19 41 41 H 1.08996 * 109.46939 * 60.00272 * 21 20 19 42 42 H 1.08993 * 109.47740 * 299.99810 * 22 21 20 43 43 H 1.09005 * 109.47573 * 59.99976 * 22 21 20 44 44 H 1.09004 * 112.84864 * 43.80415 * 23 22 21 45 45 H 1.09008 * 113.61371 * 334.56366 * 24 23 22 46 46 H 1.09000 * 113.61599 * 105.46606 * 24 23 22 47 47 H 1.09003 * 113.61002 * 270.88119 * 25 24 23 48 48 H 1.08996 * 113.69418 * 139.99462 * 25 24 23 49 49 H 1.08996 * 113.61609 * 220.02031 * 26 25 24 50 50 H 1.09006 * 113.69219 * 89.13165 * 26 25 24 51 51 H 1.08997 * 109.61742 * 6.61625 * 27 13 11 52 52 H 1.09008 * 109.55955 * 246.18817 * 27 13 11 Note: An asterisk (*) indicates this parameter was or will be optimized. Cartesian coordinates (angstroms) NO. Atom X Y Z 1 7 0.0000 0.0000 0.0000 2 6 1.3163 0.0000 0.0000 3 14 2.2503 1.6177 0.0000 4 1 2.4978 2.0464 -1.4001 5 1 1.4477 2.6527 0.7000 6 1 3.5478 1.4402 0.7000 7 6 2.0103 -1.2021 0.0005 8 1 3.0903 -1.2021 0.0010 9 16 1.1293 -2.7280 0.0005 10 6 2.5158 -3.8914 0.0021 11 6 1.9873 -5.3027 0.0017 12 8 0.7910 -5.5022 0.0009 13 7 2.8440 -6.3432 0.0033 14 6 4.2959 -6.1174 0.0048 15 6 4.9166 -6.8925 -1.1617 16 6 4.5094 -8.3645 -1.0650 17 6 2.9866 -8.4817 -1.1631 18 1 2.6487 -8.0525 -2.1063 19 7 2.5989 -9.8974 -1.1018 20 6 1.3034 -10.1173 -1.7585 21 6 0.9339 -11.5996 -1.6740 22 6 -0.4157 -11.8284 -2.3574 23 6 -0.7850 -13.3108 -2.2734 24 6 -0.2230 -14.1523 -3.4314 25 6 -1.5682 -14.8975 -3.4387 26 6 -2.1846 -13.6293 -2.8252 27 6 2.3433 -7.7245 0.0023 28 1 -0.4850 0.8400 -0.0004 29 1 3.1237 -3.7310 0.8925 30 1 3.1251 -3.7317 -0.8874 31 1 4.7196 -6.4706 0.9449 32 1 4.5002 -5.0532 -0.1136 33 1 6.0026 -6.8116 -1.1152 34 1 4.5610 -6.4777 -2.1049 35 1 4.8443 -8.7720 -0.1111 36 1 4.9698 -8.9229 -1.8800 37 1 2.5814 -10.2265 -0.1481 38 1 0.5376 -9.5233 -1.2595 39 1 1.3709 -9.8183 -2.8045 40 1 1.6995 -12.1938 -2.1726 41 1 0.8662 -11.8980 -0.6278 42 1 -1.1811 -11.2342 -1.8584 43 1 -0.3486 -11.5298 -3.4036 44 1 -0.5959 -13.7390 -1.2890 45 1 -0.0373 -13.5776 -4.3388 46 1 0.6219 -14.7790 -3.1458 47 1 -1.5927 -15.7607 -2.7735 48 1 -1.9384 -15.1233 -4.4387 49 1 -2.9233 -13.8338 -2.0503 50 1 -2.5359 -12.9107 -3.5658 51 1 1.2600 -7.7209 -0.1180 52 1 2.6071 -8.2089 0.9425 RHF calculation, no. of doubly occupied orbitals= 63 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). S: (AM1): M.J.S. DEWAR, Y.-C. YUAN, INORGANIC CHEM., 29, 589 (1990) 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 ZINC001207430648.mol2 52 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:09:33 Heat of formation + Delta-G solvation = -39.445962 kcal Electronic energy + Delta-G solvation = -26772.476754 eV Core-core repulsion = 23090.556861 eV Total energy + Delta-G solvation = -3681.919893 eV No. of doubly occupied orbitals = 63 Molecular weight (most abundant/longest-lived isotopes) = 340.193 amu Computer time = 0.69 seconds Orbital eigenvalues (eV) -40.07339 -39.32821 -37.82507 -35.97868 -33.48853 -32.44129 -32.14972 -30.33834 -28.73074 -26.56286 -25.91427 -25.58117 -24.97926 -23.36239 -22.30499 -21.59694 -20.66627 -19.69316 -19.07450 -17.69480 -16.76634 -16.42644 -15.98261 -15.57774 -15.40784 -15.21626 -14.82375 -14.80720 -14.57310 -14.42930 -14.00093 -13.76943 -13.42468 -13.21183 -13.12108 -12.88226 -12.78543 -12.72191 -12.49123 -12.09656 -11.95493 -11.85958 -11.51544 -11.49086 -11.10291 -11.09277 -11.07030 -11.03033 -10.94377 -10.89119 -10.78284 -10.75360 -10.42840 -10.24784 -9.81906 -9.41585 -9.07669 -8.92286 -8.32749 -8.23564 -6.51379 -6.22537 -3.93313 2.89743 3.02968 3.21575 3.24833 3.29443 3.55992 3.84727 4.19273 4.26173 4.31971 4.37905 4.52135 4.59978 4.64321 4.79062 4.82890 4.85817 4.90076 5.02444 5.11366 5.12638 5.13684 5.27264 5.30818 5.37408 5.40496 5.41826 5.44825 5.49313 5.53446 5.56362 5.58898 5.62033 5.63202 5.66776 5.70794 5.75314 5.84959 5.87399 5.90932 5.98585 6.02086 6.07744 6.10480 6.13896 6.19781 6.52251 6.83250 7.11696 7.54030 7.58712 8.34160 8.53850 9.22695 10.84893 Molecular weight = 340.19amu Principal moments of inertia in cm(-1) A = 0.016338 B = 0.001831 C = 0.001681 Principal moments of inertia in units of 10**(-40)*gram-cm**2 A = 1713.391017 B =15286.744095 C =16649.879263 Net atomic charges, atomic populations, and dipole contributions calculated with CM2 Atom NO. Type Charge No. of electrons 1 N -0.867 5.867 2 C 0.074 3.926 3 Si 0.919 3.081 4 H -0.277 1.277 5 H -0.282 1.282 6 H -0.272 1.272 7 C -0.535 4.535 8 H 0.073 0.927 9 S -0.062 6.062 10 C -0.150 4.150 11 C 0.534 3.466 12 O -0.519 6.519 13 N -0.617 5.617 14 C 0.107 3.893 15 C -0.136 4.136 16 C -0.104 4.104 17 C 0.073 3.927 18 H 0.082 0.918 19 N -0.680 5.680 20 C 0.054 3.946 21 C -0.113 4.113 22 C -0.110 4.110 23 C -0.109 4.109 24 C -0.137 4.137 25 C -0.145 4.145 26 C -0.140 4.140 27 C 0.087 3.913 28 H 0.268 0.732 29 H 0.089 0.911 30 H 0.089 0.911 31 H 0.067 0.933 32 H 0.090 0.910 33 H 0.077 0.923 34 H 0.069 0.931 35 H 0.063 0.937 36 H 0.071 0.929 37 H 0.339 0.661 38 H 0.034 0.966 39 H 0.072 0.928 40 H 0.070 0.930 41 H 0.059 0.941 42 H 0.064 0.936 43 H 0.064 0.936 44 H 0.088 0.912 45 H 0.080 0.920 46 H 0.064 0.936 47 H 0.077 0.923 48 H 0.065 0.935 49 H 0.064 0.936 50 H 0.080 0.920 51 H 0.089 0.911 52 H 0.061 0.939 For the charges calculated by CM2: Dipole moment (debyes) X Y Z Total from point charges 6.509 -27.231 -4.091 28.296 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 N -0.503 5.503 2 C -0.136 4.136 3 Si 0.745 3.255 4 H -0.203 1.203 5 H -0.207 1.207 6 H -0.197 1.197 7 C -0.684 4.684 8 H 0.091 0.909 9 S 0.165 5.835 10 C -0.300 4.300 11 C 0.322 3.678 12 O -0.393 6.393 13 N -0.353 5.353 14 C -0.017 4.017 15 C -0.175 4.175 16 C -0.142 4.142 17 C -0.037 4.037 18 H 0.100 0.900 19 N -0.310 5.310 20 C -0.075 4.075 21 C -0.151 4.151 22 C -0.148 4.148 23 C -0.128 4.128 24 C -0.174 4.174 25 C -0.182 4.182 26 C -0.177 4.177 27 C -0.038 4.038 28 H 0.078 0.922 29 H 0.107 0.893 30 H 0.108 0.892 31 H 0.085 0.915 32 H 0.109 0.891 33 H 0.096 0.904 34 H 0.087 0.913 35 H 0.081 0.919 36 H 0.090 0.910 37 H 0.157 0.843 38 H 0.052 0.948 39 H 0.090 0.910 40 H 0.088 0.912 41 H 0.078 0.922 42 H 0.082 0.918 43 H 0.082 0.918 44 H 0.106 0.894 45 H 0.099 0.901 46 H 0.083 0.917 47 H 0.096 0.904 48 H 0.083 0.917 49 H 0.083 0.917 50 H 0.099 0.901 51 H 0.108 0.892 52 H 0.079 0.921 Dipole moment (debyes) X Y Z Total from point charges 4.917 -27.592 -4.648 28.410 hybrid contribution 1.407 1.678 0.700 2.299 sum 6.324 -25.915 -3.948 26.966 Atomic orbital electron populations 1.70018 1.05962 1.27455 1.46830 1.24992 0.95025 0.99025 0.94520 0.86167 0.79414 0.81557 0.78356 1.20262 1.20742 1.19683 1.22859 0.94890 0.91616 1.58994 0.90932 1.85403 1.05029 0.96367 1.96701 1.23550 1.06463 0.95312 1.04719 1.20374 0.86204 0.86270 0.74904 1.90683 1.15249 1.85241 1.48127 1.48039 1.06496 1.05331 1.75472 1.21759 0.79837 1.01640 0.98444 1.21942 1.00742 0.95563 0.99211 1.21471 0.94414 0.95805 1.02548 1.21964 0.95909 0.87202 0.98666 0.90022 1.60677 1.30656 1.17076 1.22572 1.22070 0.89495 0.96295 0.99676 1.21529 0.97117 0.95689 1.00805 1.21542 0.96880 0.95357 1.01025 1.22549 0.96247 0.95364 0.98612 1.22980 0.96935 0.98710 0.98820 1.22814 0.97172 0.97170 1.01080 1.22959 0.96561 0.97801 1.00425 1.22255 1.00966 0.82335 0.98196 0.92181 0.89251 0.89230 0.91518 0.89135 0.90408 0.91257 0.91869 0.91041 0.84266 0.94756 0.90997 0.91168 0.92239 0.91777 0.91758 0.89404 0.90125 0.91710 0.90408 0.91654 0.91679 0.90091 0.89243 0.92051 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 N -0.87 -26.19 13.22 55.49 0.73 -25.46 16 2 C 0.07 2.10 5.49 -17.43 -0.10 2.01 16 3 Si 0.92 21.38 29.55 -169.99 -5.02 16.35 16 4 H -0.28 -6.31 7.11 56.52 0.40 -5.91 16 5 H -0.28 -6.46 7.11 56.52 0.40 -6.06 16 6 H -0.27 -6.14 7.11 56.52 0.40 -5.73 16 7 C -0.54 -15.26 9.50 30.94 0.29 -14.97 16 8 H 0.07 1.95 7.80 -52.49 -0.41 1.54 16 9 S -0.06 -1.69 20.57 -107.50 -2.21 -3.90 16 10 C -0.15 -3.07 4.50 36.01 0.16 -2.91 16 11 C 0.53 10.33 7.72 -10.98 -0.08 10.25 16 12 O -0.52 -12.03 15.27 5.56 0.08 -11.94 16 13 N -0.62 -8.88 2.97 -172.78 -0.51 -9.39 16 14 C 0.11 1.17 6.34 -3.71 -0.02 1.14 16 15 C -0.14 -1.23 6.08 -26.61 -0.16 -1.39 16 16 C -0.10 -0.92 5.49 -26.69 -0.15 -1.06 16 17 C 0.07 0.75 2.61 -67.61 -0.18 0.58 16 18 H 0.08 0.98 7.80 -51.93 -0.41 0.58 16 19 N -0.68 -6.06 6.25 -106.61 -0.67 -6.73 16 20 C 0.05 0.48 5.48 -3.82 -0.02 0.46 16 21 C -0.11 -0.84 5.47 -26.73 -0.15 -0.99 16 22 C -0.11 -0.80 5.05 -26.73 -0.14 -0.93 16 23 C -0.11 -0.73 3.65 -89.81 -0.33 -1.06 16 24 C -0.14 -0.87 7.63 -25.92 -0.20 -1.06 16 25 C -0.14 -0.87 8.08 -25.92 -0.21 -1.08 16 26 C -0.14 -0.92 7.64 -25.92 -0.20 -1.11 16 27 C 0.09 1.08 5.69 -3.73 -0.02 1.06 16 28 H 0.27 7.52 8.31 -40.82 -0.34 7.19 16 29 H 0.09 1.70 7.94 -51.92 -0.41 1.28 16 30 H 0.09 1.70 7.69 -51.92 -0.40 1.30 16 31 H 0.07 0.63 8.14 -51.93 -0.42 0.21 16 32 H 0.09 1.02 5.93 -51.93 -0.31 0.71 16 33 H 0.08 0.56 8.14 -51.93 -0.42 0.14 16 34 H 0.07 0.70 8.14 -51.93 -0.42 0.28 16 35 H 0.06 0.52 8.14 -51.93 -0.42 0.10 16 36 H 0.07 0.62 8.14 -51.93 -0.42 0.20 16 37 H 0.34 2.70 8.41 -39.29 -0.33 2.37 16 38 H 0.03 0.35 7.56 -51.93 -0.39 -0.04 16 39 H 0.07 0.69 8.06 -51.93 -0.42 0.27 16 40 H 0.07 0.54 8.14 -51.93 -0.42 0.12 16 41 H 0.06 0.43 8.14 -51.93 -0.42 0.00 16 42 H 0.06 0.50 8.14 -51.93 -0.42 0.08 16 43 H 0.06 0.48 7.84 -51.93 -0.41 0.07 16 44 H 0.09 0.58 8.13 -51.93 -0.42 0.16 16 45 H 0.08 0.52 7.90 -51.92 -0.41 0.11 16 46 H 0.06 0.39 8.14 -51.93 -0.42 -0.03 16 47 H 0.08 0.43 8.14 -51.93 -0.42 0.01 16 48 H 0.06 0.37 8.14 -51.93 -0.42 -0.05 16 49 H 0.06 0.41 8.14 -51.93 -0.42 -0.01 16 50 H 0.08 0.54 8.10 -51.93 -0.42 0.12 16 51 H 0.09 1.32 6.48 -51.93 -0.34 0.98 16 52 H 0.06 0.67 8.04 -51.92 -0.42 0.26 16 LS Contribution 419.31 15.07 6.32 6.32 Total: -1.00 -33.13 419.31 -12.47 -45.59 By element: Atomic # 1 Polarization: 9.92 SS G_CDS: -9.70 Total: 0.22 kcal Atomic # 6 Polarization: -9.58 SS G_CDS: -1.49 Total: -11.07 kcal Atomic # 7 Polarization: -41.13 SS G_CDS: -0.45 Total: -41.58 kcal Atomic # 8 Polarization: -12.03 SS G_CDS: 0.08 Total: -11.94 kcal Atomic # 14 Polarization: 21.38 SS G_CDS: -5.02 Total: 16.35 kcal Atomic # 16 Polarization: -1.69 SS G_CDS: -2.21 Total: -3.90 kcal Total LS contribution 6.32 Total: 6.32 kcal Total: -33.13 -12.47 -45.59 kcal The number of atoms in the molecule is 52 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. ZINC001207430648.mol2 52 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 6.147 kcal (2) G-P(sol) polarization free energy of solvation -33.127 kcal (3) G-ENP(sol) elect.-nuc.-pol. free energy of system -26.980 kcal (4) G-CDS(sol) cavity-dispersion-solvent structure free energy -12.466 kcal (5) G-P-CDS(sol) = G-P(sol) + G-CDS(sol) = (2) + (4) -45.593 kcal (6) G-S(sol) free energy of system = (1) + (5) -39.446 kcal 1SCF WAS SPECIFIED, SO GEOMETRY OPTIMIZATION WAS NOT USED Total computer time = 0.69 seconds