Wall clock time and date at job start Sat Mar 6 2021 16:13:43 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= 0 * 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 O 2 2 C 1.42899 * 1 3 3 Si 1.86298 * 109.47392 * 2 1 4 4 C 1.50703 * 109.46969 * 120.00532 * 2 1 3 5 5 C 1.32683 * 124.68388 * 64.99362 * 4 2 1 6 6 C 1.44520 * 113.09732 * 179.97438 * 5 4 2 7 7 C 1.47010 * 123.49134 * 180.02562 * 6 5 4 8 8 O 1.21657 * 119.99613 * 179.97438 * 7 6 5 9 9 N 1.34774 * 120.00268 * 359.97438 * 7 6 5 10 10 C 1.39191 * 119.99999 * 175.56257 * 9 7 6 11 11 C 1.40292 * 120.16840 * 215.52507 * 10 9 7 12 12 C 1.47497 * 120.19914 * 0.02562 * 11 10 9 13 13 O 1.21589 * 120.00147 * 5.44089 * 12 11 10 14 14 N 1.34777 * 119.99885 * 185.43438 * 12 11 10 15 15 C 1.46505 * 119.99957 * 179.97438 * 14 12 11 16 16 H 1.08995 * 109.47211 * 325.00667 * 15 14 12 17 17 C 1.52995 * 109.46716 * 85.00670 * 15 14 12 18 18 C 1.53005 * 109.47360 * 180.02562 * 17 15 14 19 19 C 1.53007 * 109.46566 * 300.00130 * 18 17 15 20 20 H 1.09002 * 109.46652 * 299.99134 * 19 18 17 21 21 O 1.42897 * 109.46799 * 179.97438 * 19 18 17 22 22 C 1.52995 * 109.47388 * 59.99887 * 19 18 17 23 23 C 1.52999 * 109.46934 * 205.00153 * 15 14 12 24 24 C 1.39691 * 119.59832 * 180.02562 * 11 10 9 25 25 C 1.37791 * 119.92552 * 359.97438 * 24 11 10 26 26 C 1.38574 * 120.33289 * 0.02776 * 25 24 11 27 27 C 1.38011 * 120.41217 * 359.67377 * 26 25 24 28 28 C 1.36251 * 113.00938 * 359.97438 * 6 5 4 29 29 Cl 1.73601 * 124.60547 * 180.02562 * 28 6 5 30 30 S 1.70450 * 110.78533 * 359.75512 * 28 6 5 31 31 H 0.96703 * 113.99848 * 300.00077 * 1 2 3 32 32 H 1.09001 * 109.47304 * 239.99988 * 2 1 3 33 33 H 1.48493 * 110.00094 * 299.99839 * 3 2 1 34 34 H 1.48505 * 109.99677 * 61.32192 * 3 2 1 35 35 H 1.07999 * 123.44973 * 359.94492 * 5 4 2 36 36 H 0.97004 * 119.99761 * 355.57098 * 9 7 6 37 37 H 0.97001 * 119.99981 * 0.02562 * 14 12 11 38 38 H 1.09004 * 109.47544 * 299.99575 * 17 15 14 39 39 H 1.09002 * 109.47839 * 59.99702 * 17 15 14 40 40 H 1.08998 * 109.47420 * 179.97438 * 18 17 15 41 41 H 1.09001 * 109.46595 * 59.99524 * 18 17 15 42 42 H 0.96703 * 114.00369 * 59.99460 * 21 19 18 43 43 H 1.08997 * 109.46996 * 180.02562 * 22 19 18 44 44 H 1.08995 * 109.47283 * 59.99635 * 22 19 18 45 45 H 1.08995 * 109.47591 * 300.00281 * 23 15 14 46 46 H 1.08998 * 109.47083 * 60.00562 * 23 15 14 47 47 H 1.07999 * 120.03822 * 180.02562 * 24 11 10 48 48 H 1.08001 * 119.83787 * 180.02562 * 25 24 11 49 49 H 1.08005 * 119.79310 * 179.97438 * 26 25 24 50 50 H 1.07992 * 119.97147 * 180.30405 * 27 26 25 Note: An asterisk (*) indicates this parameter was or will be optimized. Cartesian coordinates (angstroms) NO. Atom X Y Z 1 8 0.0000 0.0000 0.0000 2 6 1.4290 0.0000 0.0000 3 14 2.0501 1.7564 0.0000 4 6 1.9313 -0.7105 1.2304 5 6 1.7481 -1.9997 1.4855 6 6 2.3362 -2.4191 2.7373 7 6 2.2700 -3.7978 3.2431 8 8 2.7885 -4.0811 4.3065 9 7 1.6349 -4.7480 2.5286 10 6 1.6511 -6.0710 2.9608 11 6 1.7037 -7.1170 2.0274 12 6 1.7428 -6.8211 0.5830 13 8 1.8293 -5.6703 0.2002 14 7 1.6827 -7.8264 -0.3127 15 6 1.7209 -7.5326 -1.7474 16 1 1.1986 -6.5956 -1.9404 17 6 3.1759 -7.4088 -2.2040 18 6 3.2159 -7.1012 -3.7023 19 6 2.5338 -8.2330 -4.4735 20 1 3.0562 -9.1701 -4.2809 21 8 2.5707 -7.9455 -5.8728 22 6 1.0788 -8.3568 -4.0170 23 6 1.0391 -8.6648 -2.5182 24 6 1.7190 -8.4404 2.4743 25 6 1.6828 -8.7102 3.8251 26 6 1.6308 -7.6747 4.7445 27 6 1.6214 -6.3620 4.3185 28 6 2.9593 -1.4043 3.3994 29 17 3.7649 -1.5725 4.9280 30 16 2.8466 0.0512 2.5195 31 1 -0.3933 0.4417 0.7651 32 1 1.7924 -0.5138 -0.8900 33 1 1.5616 2.4678 -1.2084 34 1 1.5880 2.4585 1.2243 35 1 1.2138 -2.6671 0.8256 36 1 1.1666 -4.5098 1.7132 37 1 1.6134 -8.7445 -0.0073 38 1 3.6617 -6.6022 -1.6549 39 1 3.6983 -8.3458 -2.0110 40 1 4.2525 -7.0126 -4.0275 41 1 2.6936 -6.1640 -3.8948 42 1 3.4637 -7.8531 -6.2320 43 1 0.5931 -9.1634 -4.5661 44 1 0.5566 -7.4197 -4.2096 45 1 1.5614 -9.6017 -2.3253 46 1 0.0025 -8.7530 -2.1930 47 1 1.7598 -9.2511 1.7620 48 1 1.6951 -9.7339 4.1691 49 1 1.6024 -7.8962 5.8012 50 1 1.5857 -5.5605 5.0413 There are 74 doubly filled levels RHF calculation, no. of doubly occupied orbitals= 74 No. of singly occupied orbitals= 1 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) Cl: (AM1): M.J.S. DEWAR AND E. G. ZOEBISCH, THEOCHEM, 180, 1 (1988). 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=0 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 s_22_11560734_13431040.mol2 50 1SCF WAS SPECIFIED, SO GEOMETRY OPTIMIZATION WAS NOT USED SCF FIELD WAS ACHIEVED Wall clock time and date at job start Sat Mar 6 2021 16:13:43 Heat of formation + Delta-G solvation = -88.667735 kcal Electronic energy + Delta-G solvation = -38267.128023 eV Core-core repulsion = 33185.446889 eV Total energy + Delta-G solvation = -5081.681134 eV No. of doubly occupied orbitals = 74 and no. of open levels = 1 Molecular weight (most abundant/longest-lived isotopes) = 437.094 amu Computer time = 1.45 seconds Orbital eigenvalues (eV) -42.05958 -41.00721 -39.97968 -39.19411 -38.75718 -37.92237 -37.06629 -37.03703 -35.79296 -34.54208 -32.29474 -31.97443 -31.82728 -30.74188 -29.31225 -26.96984 -26.17996 -25.00703 -23.97860 -23.76039 -22.82866 -22.14437 -21.64621 -21.10968 -20.58335 -19.40241 -19.16171 -18.48550 -18.25033 -18.01234 -17.50662 -17.38047 -16.74569 -16.53219 -16.36568 -16.17542 -15.95511 -15.79843 -15.74992 -15.56388 -15.40513 -15.18727 -14.81563 -14.70340 -14.53913 -14.40402 -14.28482 -14.15180 -13.95855 -13.60409 -13.43462 -13.34545 -13.15572 -13.02918 -12.96301 -12.65291 -12.56111 -12.50527 -12.45298 -12.27923 -11.96459 -11.84682 -11.84182 -11.49492 -11.38409 -10.99481 -10.88774 -10.80882 -10.63159 -10.23919 -9.95461 -9.89818 -9.74629 -9.14979 -5.17657 -0.64079 -0.21842 -0.19023 0.27314 0.45025 0.76497 1.07222 1.08657 1.42454 1.79962 1.82026 1.98664 2.06907 2.24150 2.84693 3.12471 3.13294 3.38729 3.46284 3.60327 3.61880 3.70863 3.78441 3.89878 3.92206 4.01863 4.02782 4.09414 4.17999 4.24628 4.28648 4.31058 4.34105 4.36976 4.43189 4.50043 4.55834 4.73608 4.77739 4.81433 4.86554 4.92692 4.93537 4.99387 5.02544 5.11082 5.16348 5.21148 5.23780 5.38971 5.40524 5.48967 5.91409 6.20917 6.30856 6.45783 6.64511 7.04476 7.14888 Net atomic charges, atomic populations, and dipole contributions calculated with CM2 Atom NO. Type Charge No. of electrons 1 O -0.520 6.520 2 C -0.032 4.032 3 Si 0.646 3.354 4 C -0.235 4.235 5 C -0.063 4.063 6 C -0.159 4.159 7 C 0.579 3.421 8 O -0.504 6.504 9 N -0.652 5.652 10 C 0.213 3.787 11 C -0.174 4.174 12 C 0.568 3.432 13 O -0.541 6.541 14 N -0.721 5.721 15 C 0.144 3.856 16 H 0.093 0.907 17 C -0.127 4.127 18 C -0.149 4.149 19 C 0.103 3.897 20 H 0.056 0.944 21 O -0.568 6.568 22 C -0.109 4.109 23 C -0.124 4.124 24 C -0.067 4.067 25 C -0.143 4.143 26 C -0.070 4.070 27 C -0.132 4.132 28 C -0.112 4.112 29 Cl 0.027 6.973 30 S 0.241 5.759 31 H 0.378 0.622 32 H 0.132 0.868 33 H -0.229 1.229 34 H -0.235 1.235 35 H 0.161 0.839 36 H 0.422 0.578 37 H 0.401 0.599 38 H 0.073 0.927 39 H 0.066 0.934 40 H 0.069 0.931 41 H 0.074 0.926 42 H 0.380 0.620 43 H 0.077 0.923 44 H 0.074 0.926 45 H 0.072 0.928 46 H 0.078 0.922 47 H 0.129 0.871 48 H 0.131 0.869 49 H 0.133 0.867 50 H 0.148 0.852 For the charges calculated by CM2: Dipole moment (debyes) X Y Z Total from point charges 0.517 -2.376 -1.832 3.045 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 O -0.327 6.327 2 C -0.144 4.144 3 Si 0.544 3.456 4 C -0.352 4.352 5 C -0.092 4.092 6 C -0.177 4.177 7 C 0.366 3.634 8 O -0.379 6.379 9 N -0.301 5.301 10 C 0.115 3.885 11 C -0.180 4.180 12 C 0.357 3.643 13 O -0.420 6.420 14 N -0.371 5.371 15 C 0.040 3.960 16 H 0.111 0.889 17 C -0.166 4.166 18 C -0.187 4.187 19 C 0.044 3.956 20 H 0.074 0.926 21 O -0.374 6.374 22 C -0.147 4.147 23 C -0.163 4.163 24 C -0.086 4.086 25 C -0.162 4.162 26 C -0.089 4.089 27 C -0.152 4.152 28 C -0.259 4.259 29 Cl 0.056 6.944 30 S 0.506 5.494 31 H 0.227 0.773 32 H 0.150 0.850 33 H -0.152 1.152 34 H -0.159 1.159 35 H 0.178 0.822 36 H 0.265 0.735 37 H 0.236 0.764 38 H 0.091 0.909 39 H 0.085 0.915 40 H 0.088 0.912 41 H 0.093 0.907 42 H 0.227 0.773 43 H 0.095 0.905 44 H 0.093 0.907 45 H 0.090 0.910 46 H 0.097 0.903 47 H 0.146 0.854 48 H 0.149 0.851 49 H 0.151 0.849 50 H 0.166 0.834 Dipole moment (debyes) X Y Z Total from point charges 0.768 0.300 -1.497 1.709 hybrid contribution -0.901 -2.184 0.206 2.371 sum -0.133 -1.884 -1.291 2.287 Atomic orbital electron populations 1.86654 1.16929 1.79872 1.49250 1.24420 0.84000 1.10112 0.95891 0.95218 0.91800 0.74995 0.83575 1.24738 1.09005 0.98600 1.02831 1.21082 0.98573 0.94988 0.94582 1.20041 1.06592 0.94017 0.97085 1.17503 0.77995 0.83945 0.83953 1.90860 1.44860 1.81355 1.20847 1.43446 1.52646 1.05061 1.28922 1.17613 0.95500 0.82948 0.92435 1.19686 1.12298 0.93747 0.92268 1.18100 0.75613 0.85068 0.85495 1.90644 1.50368 1.20320 1.80695 1.45944 1.75053 1.09807 1.06337 1.21071 0.96479 0.99814 0.78646 0.88864 1.21865 0.96429 1.02674 0.95604 1.22234 1.01432 1.00059 0.95017 1.21850 0.96623 0.96305 0.80837 0.92623 1.86519 1.33640 1.94361 1.22845 1.21407 0.95481 1.03552 0.94309 1.21786 1.01667 0.98123 0.94677 1.21086 0.96868 0.94252 0.96432 1.21109 1.02516 0.99871 0.92673 1.21329 0.95202 0.91956 1.00380 1.21388 1.01594 0.98050 0.94191 1.25542 1.03596 1.00591 0.96185 1.98440 1.75691 1.97303 1.22939 1.83355 1.50076 1.06093 1.09913 0.77348 0.85035 1.15220 1.15917 0.82201 0.73484 0.76358 0.90853 0.91548 0.91243 0.90686 0.77280 0.90471 0.90688 0.90967 0.90327 0.85353 0.85066 0.84935 0.83427 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 O -0.52 -0.01 12.88 -35.23 -0.45 -0.46 16 2 C -0.03 0.04 2.94 36.01 0.11 0.14 16 3 Si 0.65 -1.13 30.92 -169.99 -5.26 -6.39 16 4 C -0.23 0.06 5.65 -36.20 -0.20 -0.14 16 5 C -0.06 -0.06 9.75 -39.00 -0.38 -0.44 16 6 C -0.16 -0.43 5.98 -104.58 -0.63 -1.06 16 7 C 0.58 2.44 7.62 -12.71 -0.10 2.34 16 8 O -0.50 -3.12 12.35 5.23 0.06 -3.06 16 9 N -0.65 -1.96 4.51 -10.77 -0.05 -2.01 16 10 C 0.21 0.62 6.30 -83.51 -0.53 0.09 16 11 C -0.17 -0.24 5.89 -104.81 -0.62 -0.85 16 12 C 0.57 0.83 7.71 -12.49 -0.10 0.74 16 13 O -0.54 -1.98 13.65 5.29 0.07 -1.91 16 14 N -0.72 0.95 5.03 -54.86 -0.28 0.67 16 15 C 0.14 -0.17 2.45 -67.93 -0.17 -0.34 16 16 H 0.09 0.06 7.58 -51.93 -0.39 -0.34 16 17 C -0.13 0.12 5.49 -26.73 -0.15 -0.03 16 18 C -0.15 0.18 5.55 -26.72 -0.15 0.03 16 19 C 0.10 -0.10 3.47 -26.73 -0.09 -0.20 16 20 H 0.06 -0.09 8.14 -51.93 -0.42 -0.52 16 21 O -0.57 -0.23 13.32 -35.23 -0.47 -0.69 16 22 C -0.11 0.13 5.57 -26.70 -0.15 -0.02 16 23 C -0.12 0.33 5.39 -26.70 -0.14 0.18 16 24 C -0.07 0.06 9.52 -39.20 -0.37 -0.31 16 25 C -0.14 0.12 10.02 -39.62 -0.40 -0.27 16 26 C -0.07 -0.03 10.03 -39.54 -0.40 -0.42 16 27 C -0.13 -0.33 8.68 -39.43 -0.34 -0.68 16 28 C -0.11 -0.27 6.88 30.02 0.21 -0.07 16 29 Cl 0.03 0.09 27.20 -51.86 -1.41 -1.32 16 30 S 0.24 -0.05 20.18 -107.50 -2.17 -2.22 16 31 H 0.38 -1.08 8.70 45.56 0.40 -0.69 16 32 H 0.13 -0.26 7.96 -51.93 -0.41 -0.67 16 33 H -0.23 -0.22 7.11 56.52 0.40 0.18 16 34 H -0.24 0.08 7.11 56.52 0.40 0.48 16 35 H 0.16 0.05 7.10 -52.49 -0.37 -0.32 16 36 H 0.42 0.95 5.15 -40.82 -0.21 0.74 16 37 H 0.40 -1.47 6.66 -40.82 -0.27 -1.74 16 38 H 0.07 0.08 8.14 -51.93 -0.42 -0.34 16 39 H 0.07 -0.10 8.14 -51.93 -0.42 -0.52 16 40 H 0.07 -0.13 8.14 -51.93 -0.42 -0.55 16 41 H 0.07 -0.01 8.14 -51.93 -0.42 -0.43 16 42 H 0.38 -1.27 9.08 45.56 0.41 -0.86 16 43 H 0.08 -0.08 8.14 -51.93 -0.42 -0.50 16 44 H 0.07 -0.02 8.14 -51.93 -0.42 -0.45 16 45 H 0.07 -0.26 8.14 -51.93 -0.42 -0.68 16 46 H 0.08 -0.23 8.14 -51.93 -0.42 -0.65 16 47 H 0.13 -0.51 6.44 -52.49 -0.34 -0.85 16 48 H 0.13 -0.40 8.06 -52.49 -0.42 -0.83 16 49 H 0.13 -0.20 8.06 -52.48 -0.42 -0.63 16 50 H 0.15 0.52 6.23 -52.49 -0.33 0.19 16 LS Contribution 433.44 15.07 6.53 6.53 Total: 0.00 -8.74 433.44 -13.37 -22.11 By element: Atomic # 1 Polarization: -4.59 SS G_CDS: -5.36 Total: -9.95 kcal Atomic # 6 Polarization: 3.30 SS G_CDS: -4.59 Total: -1.30 kcal Atomic # 7 Polarization: -1.02 SS G_CDS: -0.32 Total: -1.34 kcal Atomic # 8 Polarization: -5.33 SS G_CDS: -0.79 Total: -6.12 kcal Atomic # 14 Polarization: -1.13 SS G_CDS: -5.26 Total: -6.39 kcal Atomic # 16 Polarization: -0.05 SS G_CDS: -2.17 Total: -2.22 kcal Atomic # 17 Polarization: 0.09 SS G_CDS: -1.41 Total: -1.32 kcal Total LS contribution 6.53 Total: 6.53 kcal Total: -8.74 -13.37 -22.11 kcal The number of atoms in the molecule is 50 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. s_22_11560734_13431040.mol2 50 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 -66.561 kcal (2) G-P(sol) polarization free energy of solvation -8.736 kcal (3) G-ENP(sol) elect.-nuc.-pol. free energy of system -75.297 kcal (4) G-CDS(sol) cavity-dispersion-solvent structure free energy -13.370 kcal (5) G-P-CDS(sol) = G-P(sol) + G-CDS(sol) = (2) + (4) -22.106 kcal (6) G-S(sol) free energy of system = (1) + (5) -88.668 kcal 1SCF WAS SPECIFIED, SO GEOMETRY OPTIMIZATION WAS NOT USED Total computer time = 1.45 seconds