Wall clock time and date at job start Sat Mar 6 2021 03:31:11 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.42901 * 1 3 3 Si 1.86298 * 109.46997 * 2 1 4 4 C 1.50699 * 109.46830 * 240.00466 * 2 1 3 5 5 C 1.34459 * 126.48165 * 104.68373 * 4 2 1 6 6 N 1.36566 * 109.98310 * 179.87000 * 5 4 2 7 7 C 1.37826 * 109.86321 * 0.40966 * 6 5 4 8 8 C 1.38643 * 133.38914 * 179.76051 * 7 6 5 9 9 C 1.39319 * 119.67050 * 180.02562 * 8 7 6 10 10 C 1.47581 * 119.83011 * 179.97438 * 9 8 7 11 11 O 1.21589 * 119.99851 * 359.97438 * 10 9 8 12 12 N 1.34773 * 119.99791 * 180.02562 * 10 9 8 13 13 C 1.46506 * 119.99901 * 180.02562 * 12 10 9 14 14 H 1.09000 * 112.84692 * 336.18968 * 13 12 10 15 15 C 1.53777 * 113.61592 * 107.49937 * 13 12 10 16 16 C 1.53779 * 87.08097 * 89.11772 * 15 13 12 17 17 H 1.09002 * 113.60798 * 270.87923 * 16 15 13 18 18 C 1.50697 * 113.61522 * 139.98456 * 16 15 13 19 19 C 1.38003 * 119.94303 * 142.22978 * 18 16 15 20 20 C 1.39329 * 119.88156 * 179.73157 * 19 18 16 21 21 C 1.48461 * 120.11897 * 180.26964 * 20 19 18 22 22 N 1.31811 * 126.67493 * 179.97438 * 21 20 19 23 23 N 1.28876 * 107.81890 * 179.89487 * 22 21 20 24 24 N 1.28732 * 109.31235 * 0.37149 * 23 22 21 25 25 N 1.28984 * 108.94929 * 359.77024 * 24 23 22 26 26 C 1.39282 * 119.76317 * 0.50919 * 20 19 18 27 27 C 1.38047 * 119.88274 * 359.76120 * 26 20 19 28 28 C 1.38296 * 120.11734 * 359.97438 * 27 26 20 29 29 C 1.53780 * 87.08320 * 25.43170 * 16 15 13 30 30 C 1.40199 * 120.33619 * 359.97438 * 9 8 7 31 31 C 1.36226 * 120.37662 * 0.02562 * 30 9 8 32 32 C 1.39869 * 119.96149 * 359.96333 * 31 30 9 33 33 H 0.96708 * 114.00246 * 59.99398 * 1 2 3 34 34 H 1.08991 * 109.47274 * 120.00272 * 2 1 3 35 35 H 1.48506 * 109.99834 * 60.00805 * 3 2 1 36 36 H 1.48502 * 110.00011 * 181.32069 * 3 2 1 37 37 H 1.07999 * 125.00713 * 0.02562 * 5 4 2 38 38 H 0.97001 * 125.07229 * 180.14045 * 6 5 4 39 39 H 1.08003 * 120.16191 * 359.97438 * 8 7 6 40 40 H 0.96995 * 120.00559 * 359.97438 * 12 10 9 41 41 H 1.09006 * 113.60744 * 203.66420 * 15 13 12 42 42 H 1.09000 * 113.61219 * 334.56405 * 15 13 12 43 43 H 1.07995 * 120.06297 * 359.96882 * 19 18 16 44 44 H 0.96999 * 125.34066 * 180.14126 * 23 22 21 45 45 H 1.07998 * 120.05902 * 179.77814 * 26 20 19 46 46 H 1.07994 * 119.94594 * 180.02562 * 27 26 20 47 47 H 1.07993 * 119.88295 * 180.02562 * 28 27 26 48 48 H 1.09001 * 113.61282 * 220.01907 * 29 16 15 49 49 H 1.08999 * 113.61590 * 89.11938 * 29 16 15 50 50 H 1.08005 * 119.80992 * 180.02562 * 30 9 8 51 51 H 1.08002 * 120.01986 * 180.02562 * 31 30 9 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.0500 1.7565 0.0000 4 6 1.9313 -0.7103 -1.2305 5 6 2.4561 -0.1357 -2.3270 6 7 2.8025 -1.0944 -3.2359 7 6 2.4940 -2.3429 -2.7402 8 6 2.6302 -3.6277 -3.2433 9 6 2.2122 -4.7152 -2.4794 10 6 2.3550 -6.0843 -3.0117 11 8 2.8360 -6.2608 -4.1144 12 7 1.9513 -7.1388 -2.2759 13 6 2.0935 -8.4980 -2.8040 14 1 2.1348 -8.5268 -3.8928 15 6 3.2066 -9.3106 -2.1217 16 6 2.1458 -9.7602 -1.1032 17 1 2.0544 -9.0941 -0.2453 18 6 2.2628 -11.2103 -0.7101 19 6 2.0210 -11.5921 0.5939 20 6 2.1241 -12.9339 0.9547 21 6 1.8593 -13.3514 2.3546 22 7 1.9182 -14.5774 2.8351 23 7 1.6243 -14.5242 4.0887 24 7 1.3736 -13.3046 4.4159 25 7 1.5082 -12.5555 3.3745 26 6 2.4812 -13.8816 -0.0015 27 6 2.7264 -13.4881 -1.3018 28 6 2.6173 -12.1557 -1.6558 29 6 1.0979 -9.5015 -2.1985 30 6 1.6566 -4.5157 -1.2078 31 6 1.5192 -3.2568 -0.7059 32 6 1.9343 -2.1547 -1.4606 33 1 -0.3934 0.4418 -0.7650 34 1 1.7923 -0.5138 0.8899 35 1 1.5616 2.4678 1.2086 36 1 3.5346 1.7703 -0.0322 37 1 2.5869 0.9271 -2.4672 38 1 3.2021 -0.9221 -4.1028 39 1 3.0574 -3.7845 -4.2227 40 1 1.5673 -6.9981 -1.3964 41 1 3.5934 -10.1232 -2.7368 42 1 3.9943 -8.6942 -1.6886 43 1 1.7487 -10.8528 1.3324 44 1 1.5934 -15.2843 4.6905 45 1 2.5660 -14.9225 0.2738 46 1 3.0030 -14.2221 -2.0441 47 1 2.8095 -11.8521 -2.6742 48 1 0.1833 -9.0391 -1.8272 49 1 0.9082 -10.3662 -2.8344 50 1 1.3358 -5.3650 -0.6227 51 1 1.0908 -3.1097 0.2746 There are 75 doubly filled levels RHF calculation, no. of doubly occupied orbitals= 75 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). 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 m_280130_17857610_17904044.mol2 51 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 03:31:11 Heat of formation + Delta-G solvation = 176.798227 kcal Electronic energy + Delta-G solvation = -36936.703822 eV Core-core repulsion = 31999.937881 eV Total energy + Delta-G solvation = -4936.765941 eV No. of doubly occupied orbitals = 75 and no. of open levels = 1 Molecular weight (most abundant/longest-lived isotopes) = 417.159 amu Computer time = 1.69 seconds Orbital eigenvalues (eV) -45.39006 -42.47188 -41.58727 -40.39753 -38.97070 -37.82279 -37.40624 -35.38315 -34.83400 -34.17087 -33.23431 -32.11495 -31.44472 -30.69895 -30.38685 -26.77193 -26.68154 -26.39287 -25.61405 -25.42387 -24.48200 -23.60581 -23.07149 -22.14904 -21.51695 -20.81211 -20.39885 -19.80800 -19.42470 -18.42905 -18.07346 -17.94490 -17.60412 -17.08561 -16.93748 -16.57696 -16.36785 -16.20001 -15.89321 -15.73742 -15.55628 -15.38439 -15.18530 -15.02402 -14.81475 -14.69930 -14.59237 -14.47187 -14.33501 -14.15999 -13.99727 -13.56323 -13.37937 -13.30514 -13.10039 -12.87366 -12.59836 -12.50170 -12.32021 -12.29846 -12.22022 -12.12955 -12.01053 -11.90473 -11.79168 -11.60314 -11.26471 -10.89449 -10.73481 -10.34277 -10.05323 -9.79672 -9.57151 -9.33533 -8.79647 -5.09496 -0.26225 -0.24314 0.24934 0.36002 0.41732 1.29622 1.30256 1.38323 1.39178 1.62986 1.88357 1.98086 2.03020 2.22655 2.70388 2.75513 2.95364 3.21973 3.40966 3.55982 3.59184 3.65943 3.68524 3.76224 3.77950 3.86502 3.88309 3.95666 4.07665 4.11409 4.15171 4.23826 4.33271 4.40994 4.44106 4.46192 4.50432 4.56816 4.57395 4.64905 4.68517 4.74802 4.83014 4.84160 4.94393 4.99172 5.00408 5.04089 5.10873 5.25913 5.37079 5.43984 5.45670 5.47023 5.47325 5.55317 5.62654 5.65506 5.88678 6.00407 6.37651 6.51764 6.66202 7.19957 7.26476 Net atomic charges, atomic populations, and dipole contributions calculated with CM2 Atom NO. Type Charge No. of electrons 1 O -0.542 6.542 2 C -0.030 4.030 3 Si 0.684 3.316 4 C -0.211 4.211 5 C 0.077 3.923 6 N -0.584 5.584 7 C 0.081 3.919 8 C -0.056 4.056 9 C -0.122 4.122 10 C 0.567 3.433 11 O -0.536 6.536 12 N -0.709 5.709 13 C 0.120 3.880 14 H 0.107 0.893 15 C -0.144 4.144 16 C -0.074 4.074 17 H 0.100 0.900 18 C -0.072 4.072 19 C -0.072 4.072 20 C -0.017 4.017 21 C 0.217 3.783 22 N -0.293 5.293 23 N -0.226 5.226 24 N -0.046 5.046 25 N -0.268 5.268 26 C -0.083 4.083 27 C -0.118 4.118 28 C -0.095 4.095 29 C -0.136 4.136 30 C -0.113 4.113 31 C -0.073 4.073 32 C -0.063 4.063 33 H 0.377 0.623 34 H 0.111 0.889 35 H -0.254 1.254 36 H -0.256 1.256 37 H 0.179 0.821 38 H 0.422 0.578 39 H 0.136 0.864 40 H 0.402 0.598 41 H 0.096 0.904 42 H 0.076 0.924 43 H 0.135 0.865 44 H 0.477 0.523 45 H 0.132 0.868 46 H 0.132 0.868 47 H 0.130 0.870 48 H 0.078 0.922 49 H 0.100 0.900 50 H 0.126 0.874 51 H 0.130 0.870 For the charges calculated by CM2: Dipole moment (debyes) X Y Z Total from point charges -0.615 0.969 -1.615 1.981 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.347 6.347 2 C -0.141 4.141 3 Si 0.586 3.414 4 C -0.219 4.219 5 C -0.052 4.052 6 N -0.186 5.186 7 C -0.025 4.025 8 C -0.077 4.077 9 C -0.126 4.126 10 C 0.355 3.645 11 O -0.414 6.414 12 N -0.361 5.361 13 C 0.014 3.986 14 H 0.125 0.875 15 C -0.183 4.183 16 C -0.092 4.092 17 H 0.119 0.881 18 C -0.072 4.072 19 C -0.092 4.092 20 C -0.020 4.020 21 C -0.062 4.062 22 N -0.145 5.145 23 N -0.072 5.072 24 N -0.041 5.041 25 N -0.141 5.141 26 C -0.103 4.103 27 C -0.136 4.136 28 C -0.114 4.114 29 C -0.174 4.174 30 C -0.133 4.133 31 C -0.091 4.091 32 C -0.068 4.068 33 H 0.224 0.776 34 H 0.129 0.871 35 H -0.179 1.179 36 H -0.181 1.181 37 H 0.196 0.804 38 H 0.260 0.740 39 H 0.154 0.846 40 H 0.238 0.762 41 H 0.115 0.885 42 H 0.095 0.905 43 H 0.153 0.847 44 H 0.332 0.668 45 H 0.150 0.850 46 H 0.149 0.851 47 H 0.147 0.853 48 H 0.097 0.903 49 H 0.118 0.882 50 H 0.144 0.856 51 H 0.148 0.852 Dipole moment (debyes) X Y Z Total from point charges 0.191 1.264 -0.888 1.556 hybrid contribution 0.237 -1.175 0.331 1.243 sum 0.428 0.089 -0.557 0.708 Atomic orbital electron populations 1.86510 1.17420 1.81745 1.49015 1.24331 0.84756 1.08849 0.96186 0.95292 0.80752 0.77378 0.88002 1.20737 1.10900 0.92143 0.98114 1.21749 0.99764 0.94329 0.89327 1.41718 1.54620 1.04000 1.18282 1.18074 1.03229 0.86316 0.94892 1.20071 0.99430 0.88069 1.00131 1.19408 1.05098 0.92277 0.95839 1.18141 0.77295 0.84488 0.84544 1.90793 1.44274 1.86181 1.20153 1.45561 1.64230 1.05295 1.20978 1.21854 0.96397 0.80551 0.99762 0.87499 1.23587 0.97006 1.00020 0.97657 1.21615 0.97534 0.92390 0.97665 0.88146 1.19582 0.99611 0.95210 0.92831 1.21064 0.97519 0.94419 0.96164 1.17204 1.04705 0.94117 0.85955 1.23733 1.00936 0.85028 0.96534 1.74925 1.21207 1.23871 0.94479 1.47763 1.42775 1.10064 1.06611 1.78953 1.14607 0.89262 1.21319 1.76098 1.16952 1.25529 0.95550 1.21189 0.97192 0.99426 0.92537 1.21096 1.00335 0.96238 0.95925 1.21223 0.97785 0.92705 0.99678 1.23443 0.97301 0.99239 0.97376 1.20685 1.00557 0.97345 0.94699 1.20203 0.98716 0.90500 0.99691 1.19179 1.00976 0.93076 0.93597 0.77589 0.87127 1.17886 1.18092 0.80441 0.74000 0.84604 0.76194 0.88519 0.90538 0.84697 0.66845 0.84980 0.85052 0.85265 0.90283 0.88186 0.85640 0.85205 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.54 -0.08 12.88 -35.23 -0.45 -0.53 16 2 C -0.03 0.01 3.06 36.01 0.11 0.12 16 3 Si 0.68 -1.32 30.76 -169.99 -5.23 -6.55 16 4 C -0.21 0.00 5.57 -102.74 -0.57 -0.57 16 5 C 0.08 -0.16 10.43 -16.78 -0.17 -0.34 16 6 N -0.58 0.86 6.07 -10.68 -0.06 0.79 16 7 C 0.08 0.10 7.26 -84.00 -0.61 -0.51 16 8 C -0.06 -0.11 9.61 -39.03 -0.37 -0.49 16 9 C -0.12 -0.23 5.87 -104.95 -0.62 -0.85 16 10 C 0.57 1.56 7.76 -12.45 -0.10 1.46 16 11 O -0.54 -3.33 15.91 5.29 0.08 -3.25 16 12 N -0.71 0.16 5.21 -54.04 -0.28 -0.12 16 13 C 0.12 -0.08 4.32 -67.12 -0.29 -0.37 16 14 H 0.11 0.10 7.46 -51.93 -0.39 -0.29 16 15 C -0.14 0.18 7.25 -25.92 -0.19 -0.01 16 16 C -0.07 0.14 3.85 -90.96 -0.35 -0.21 16 17 H 0.10 -0.24 8.09 -51.93 -0.42 -0.66 16 18 C -0.07 0.04 5.28 -104.67 -0.55 -0.51 16 19 C -0.07 -0.06 9.36 -39.26 -0.37 -0.42 16 20 C -0.02 -0.03 5.88 -104.89 -0.62 -0.64 16 21 C 0.22 0.63 7.94 -156.30 -1.24 -0.61 16 22 N -0.29 -0.51 12.08 32.17 0.39 -0.13 16 23 N -0.23 0.06 7.90 60.35 0.48 0.53 16 24 N -0.05 -0.11 13.47 60.35 0.81 0.70 16 25 N -0.27 -1.03 12.23 32.16 0.39 -0.64 16 26 C -0.08 -0.04 9.70 -39.26 -0.38 -0.42 16 27 C -0.12 0.09 10.03 -39.63 -0.40 -0.31 16 28 C -0.10 0.10 8.60 -39.52 -0.34 -0.24 16 29 C -0.14 0.30 7.29 -25.92 -0.19 0.11 16 30 C -0.11 -0.04 9.46 -39.57 -0.37 -0.41 16 31 C -0.07 -0.08 10.03 -39.70 -0.40 -0.48 16 32 C -0.06 -0.11 6.35 -105.04 -0.67 -0.78 16 33 H 0.38 -1.42 8.70 45.56 0.40 -1.02 16 34 H 0.11 0.03 7.96 -51.93 -0.41 -0.38 16 35 H -0.25 -0.42 7.11 56.52 0.40 -0.02 16 36 H -0.26 -0.15 7.11 56.52 0.40 0.25 16 37 H 0.18 -1.04 6.35 -52.49 -0.33 -1.37 16 38 H 0.42 -1.78 8.96 -40.82 -0.37 -2.15 16 39 H 0.14 0.24 7.62 -52.48 -0.40 -0.16 16 40 H 0.40 -1.00 6.69 -40.82 -0.27 -1.27 16 41 H 0.10 -0.18 6.30 -51.93 -0.33 -0.51 16 42 H 0.08 -0.01 8.14 -51.93 -0.42 -0.44 16 43 H 0.14 0.06 7.88 -52.49 -0.41 -0.35 16 44 H 0.48 -2.33 8.96 -340.36 -3.05 -5.38 16 45 H 0.13 -0.04 7.88 -52.49 -0.41 -0.45 16 46 H 0.13 -0.35 8.06 -52.49 -0.42 -0.77 16 47 H 0.13 -0.35 6.65 -52.49 -0.35 -0.70 16 48 H 0.08 -0.19 8.14 -51.93 -0.42 -0.61 16 49 H 0.10 -0.26 7.72 -51.93 -0.40 -0.66 16 50 H 0.13 -0.30 6.36 -52.48 -0.33 -0.64 16 51 H 0.13 0.02 8.06 -52.49 -0.42 -0.41 16 LS Contribution 431.64 15.07 6.50 6.50 Total: 0.00 -12.72 431.64 -14.43 -27.15 By element: Atomic # 1 Polarization: -9.62 SS G_CDS: -8.37 Total: -17.99 kcal Atomic # 6 Polarization: 2.21 SS G_CDS: -8.69 Total: -6.47 kcal Atomic # 7 Polarization: -0.58 SS G_CDS: 1.73 Total: 1.14 kcal Atomic # 8 Polarization: -3.41 SS G_CDS: -0.37 Total: -3.77 kcal Atomic # 14 Polarization: -1.32 SS G_CDS: -5.23 Total: -6.55 kcal Total LS contribution 6.50 Total: 6.50 kcal Total: -12.72 -14.43 -27.15 kcal The number of atoms in the molecule is 51 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. m_280130_17857610_17904044.mol2 51 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 203.946 kcal (2) G-P(sol) polarization free energy of solvation -12.717 kcal (3) G-ENP(sol) elect.-nuc.-pol. free energy of system 191.229 kcal (4) G-CDS(sol) cavity-dispersion-solvent structure free energy -14.431 kcal (5) G-P-CDS(sol) = G-P(sol) + G-CDS(sol) = (2) + (4) -27.148 kcal (6) G-S(sol) free energy of system = (1) + (5) 176.798 kcal 1SCF WAS SPECIFIED, SO GEOMETRY OPTIMIZATION WAS NOT USED Total computer time = 1.69 seconds