Wall clock time and date at job start Sat Apr 11 2020 06:15:22 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 N 1.46501 * 1 3 3 C 1.34771 * 119.99978 * 2 1 4 4 O 1.21586 * 120.00004 * 359.97438 * 3 2 1 5 5 N 1.34783 * 120.00207 * 179.97438 * 3 2 1 6 6 C 1.46497 * 120.00115 * 179.97438 * 5 3 2 7 7 C 1.53000 * 109.47049 * 180.02562 * 6 5 3 8 8 C 1.50697 * 109.47384 * 180.02562 * 7 6 5 9 9 O 1.21289 * 119.99965 * 0.02562 * 8 7 6 10 10 N 1.34770 * 120.00219 * 180.02562 * 8 7 6 11 11 C 1.46924 * 120.63195 * 354.36536 * 10 8 7 12 12 C 1.53624 * 108.54424 * 126.37067 * 11 10 8 13 13 C 1.53039 * 109.24416 * 54.85598 * 12 11 10 14 14 C 1.53038 * 109.53584 * 298.51509 * 13 12 11 15 15 H 1.08999 * 109.49864 * 181.31947 * 14 13 12 16 16 C 1.53004 * 109.52244 * 301.39182 * 14 13 12 17 17 C 1.46927 * 120.63036 * 174.64468 * 10 8 7 18 18 H 1.09005 * 109.58786 * 113.79867 * 17 10 8 19 19 C 1.53004 * 109.58707 * 353.37383 * 17 10 8 20 20 C 1.53007 * 109.46945 * 181.87105 * 19 17 10 21 21 N 1.46504 * 109.46781 * 185.10640 * 20 19 17 22 22 C 1.36935 * 119.99714 * 179.97438 * 21 20 19 23 23 H 1.08001 * 120.00281 * 0.02562 * 22 21 20 24 24 C 1.38815 * 119.99814 * 179.97438 * 22 21 20 25 25 N 1.31623 * 120.00034 * 0.02562 * 24 22 21 26 26 Si 1.86797 * 120.00153 * 179.97438 * 24 22 21 27 27 H 1.48501 * 109.46699 * 90.00535 * 26 24 22 28 28 H 1.48496 * 109.47680 * 210.00848 * 26 24 22 29 29 H 1.48506 * 109.47185 * 330.01041 * 26 24 22 30 30 H 1.08997 * 109.47099 * 300.00390 * 1 2 3 31 31 H 1.09002 * 109.46917 * 180.02562 * 1 2 3 32 32 H 1.09000 * 109.47071 * 60.00723 * 1 2 3 33 33 H 0.97003 * 119.99734 * 179.97438 * 2 1 3 34 34 H 0.97002 * 119.99163 * 0.02562 * 5 3 2 35 35 H 1.08992 * 109.47440 * 300.00540 * 6 5 3 36 36 H 1.09000 * 109.46925 * 60.00559 * 6 5 3 37 37 H 1.09001 * 109.46839 * 300.00408 * 7 6 5 38 38 H 1.08999 * 109.47339 * 60.00265 * 7 6 5 39 39 H 1.09004 * 109.58724 * 246.19959 * 11 10 8 40 40 H 1.08996 * 109.58347 * 6.62536 * 11 10 8 41 41 H 1.08999 * 109.62974 * 294.86683 * 12 11 10 42 42 H 1.09001 * 109.43088 * 174.72904 * 12 11 10 43 43 H 1.08992 * 109.45929 * 178.49429 * 13 12 11 44 44 H 1.08999 * 109.44795 * 58.52622 * 13 12 11 45 45 H 1.09002 * 109.46883 * 179.85590 * 16 14 13 46 46 H 1.08995 * 109.47739 * 299.86058 * 16 14 13 47 47 H 1.08999 * 109.46917 * 59.85863 * 16 14 13 48 48 H 1.09003 * 109.46950 * 301.86714 * 19 17 10 49 49 H 1.08997 * 109.47280 * 61.87330 * 19 17 10 50 50 H 1.08997 * 109.46802 * 305.10689 * 20 19 17 51 51 H 1.08997 * 109.47163 * 65.11092 * 20 19 17 52 52 H 0.97001 * 119.99753 * 359.97438 * 21 20 19 53 53 H 0.96997 * 119.99691 * 179.97438 * 25 24 22 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 7 1.4650 0.0000 0.0000 3 6 2.1389 1.1672 0.0000 4 8 1.5309 2.2201 -0.0005 5 7 3.4867 1.1672 0.0005 6 6 4.2191 2.4359 0.0000 7 6 5.7234 2.1564 0.0001 8 6 6.4769 3.4614 0.0001 9 8 5.8705 4.5118 0.0005 10 7 7.8246 3.4614 -0.0003 11 6 8.5733 2.2033 0.1233 12 6 9.5650 2.1079 -1.0460 13 6 10.4426 3.3616 -1.0666 14 6 9.5649 4.5971 -1.2792 15 1 10.1929 5.4873 -1.3172 16 6 8.7963 4.4581 -2.5949 17 6 8.5732 4.7202 -0.1182 18 1 9.1162 4.9105 0.8075 19 6 7.6023 5.8705 -0.3926 20 6 8.3886 7.1722 -0.5612 21 7 7.4546 8.2920 -0.7023 22 6 7.9318 9.5651 -0.8660 23 1 8.9966 9.7436 -0.8915 24 6 7.0468 10.6261 -1.0002 25 7 5.7490 10.4085 -0.9696 26 14 7.6976 12.3628 -1.2227 27 1 7.8681 12.6457 -2.6706 28 1 6.7347 13.3310 -0.6392 29 1 9.0077 12.4923 -0.5355 30 1 -0.3633 0.5139 0.8899 31 1 -0.3633 -1.0277 0.0005 32 1 -0.3633 0.5137 -0.8901 33 1 1.9500 -0.8401 0.0004 34 1 3.9716 0.3271 0.0005 35 1 3.9561 3.0069 -0.8904 36 1 3.9555 3.0080 0.8895 37 1 5.9867 1.5851 0.8903 38 1 5.9868 1.5844 -0.8897 39 1 9.1164 2.1925 1.0684 40 1 7.8817 1.3616 0.0874 41 1 9.0174 2.0317 -1.9853 42 1 10.1930 1.2263 -0.9179 43 1 11.1654 3.2855 -1.8787 44 1 10.9708 3.4504 -0.1173 45 1 8.1696 5.3373 -2.7442 46 1 9.5022 4.3705 -3.4207 47 1 8.1694 3.5672 -2.5569 48 1 7.0423 5.6636 -1.3046 49 1 6.9107 5.9696 0.4440 50 1 9.0183 7.3322 0.3140 51 1 9.0139 7.1062 -1.4516 52 1 6.4982 8.1316 -0.6798 53 1 5.1307 11.1500 -1.0631 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 ZINC001707149062.mol2 53 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:15:22 Heat of formation + Delta-G solvation = -43.605689 kcal Electronic energy + Delta-G solvation = -30849.659641 eV Core-core repulsion = 26757.099249 eV Total energy + Delta-G solvation = -4092.560392 eV No. of doubly occupied orbitals = 66 Molecular weight (most abundant/longest-lived isotopes) = 340.226 amu Computer time = 1.51 seconds Orbital eigenvalues (eV) -40.86569 -40.06220 -37.95899 -35.66043 -35.12729 -34.68446 -33.66191 -32.78478 -31.06555 -30.74127 -28.35137 -27.00962 -26.66102 -25.17296 -23.78345 -23.52457 -22.11848 -21.22939 -20.74917 -20.01745 -18.90214 -18.27659 -17.26842 -16.72916 -16.39786 -16.12918 -15.92519 -15.42506 -15.22002 -14.81000 -14.63759 -14.43551 -14.34127 -14.18277 -13.83669 -13.53239 -13.31652 -13.17525 -13.00250 -12.77585 -12.75455 -12.53108 -12.42241 -12.07838 -11.83138 -11.69713 -11.53154 -11.28446 -11.07219 -11.00569 -10.81741 -10.73090 -10.51314 -10.43187 -10.37636 -10.10974 -9.90394 -9.54454 -9.54271 -9.08595 -8.77665 -8.68299 -8.48460 -6.73669 -5.60059 -2.88289 2.38160 2.65683 2.98652 3.03404 3.52652 3.56639 3.59383 3.96930 4.42845 4.45645 4.60409 4.67097 4.72610 4.78405 4.82267 4.87565 4.93733 5.06078 5.12225 5.21102 5.25293 5.33156 5.50889 5.53138 5.59328 5.65589 5.69859 5.72138 5.73950 5.77131 5.83603 5.87652 6.04880 6.07544 6.13418 6.14077 6.22522 6.28694 6.35707 6.46663 6.50373 6.60824 6.64396 6.72901 7.13200 7.26060 7.45823 7.73973 7.82617 8.27611 8.40079 8.58779 9.69073 10.31695 10.66342 11.62752 Molecular weight = 340.23amu Principal moments of inertia in cm(-1) A = 0.008632 B = 0.002824 C = 0.002189 Principal moments of inertia in units of 10**(-40)*gram-cm**2 A = 3242.953468 B = 9911.434713 C =12788.548164 Net atomic charges, atomic populations, and dipole contributions calculated with CM2 Atom NO. Type Charge No. of electrons 1 C 0.107 3.893 2 N -0.754 5.754 3 C 0.699 3.301 4 O -0.580 6.580 5 N -0.747 5.747 6 C 0.151 3.849 7 C -0.153 4.153 8 C 0.513 3.487 9 O -0.507 6.507 10 N -0.608 5.608 11 C 0.108 3.892 12 C -0.137 4.137 13 C -0.115 4.115 14 C -0.105 4.105 15 H 0.088 0.912 16 C -0.146 4.146 17 C 0.158 3.842 18 H 0.070 0.930 19 C -0.119 4.119 20 C 0.149 3.851 21 N -0.731 5.731 22 C -0.236 4.236 23 H 0.068 0.932 24 C -0.017 4.017 25 N -0.931 5.931 26 Si 0.905 3.095 27 H -0.291 1.291 28 H -0.292 1.292 29 H -0.283 1.283 30 H 0.058 0.942 31 H 0.074 0.926 32 H 0.058 0.942 33 H 0.396 0.604 34 H 0.393 0.607 35 H 0.078 0.922 36 H 0.078 0.922 37 H 0.095 0.905 38 H 0.096 0.904 39 H 0.068 0.932 40 H 0.081 0.919 41 H 0.071 0.929 42 H 0.070 0.930 43 H 0.064 0.936 44 H 0.062 0.938 45 H 0.068 0.932 46 H 0.049 0.951 47 H 0.051 0.949 48 H 0.072 0.928 49 H 0.092 0.908 50 H 0.010 0.990 51 H 0.015 0.985 52 H 0.386 0.614 53 H 0.254 0.746 For the charges calculated by CM2: Dipole moment (debyes) X Y Z Total from point charges 2.770 -32.236 1.832 32.406 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.038 4.038 2 N -0.410 5.410 3 C 0.399 3.601 4 O -0.457 6.457 5 N -0.403 5.403 6 C 0.027 3.973 7 C -0.193 4.193 8 C 0.299 3.701 9 O -0.381 6.381 10 N -0.343 5.343 11 C -0.016 4.016 12 C -0.175 4.175 13 C -0.153 4.153 14 C -0.124 4.124 15 H 0.106 0.894 16 C -0.203 4.203 17 C 0.055 3.945 18 H 0.088 0.912 19 C -0.156 4.156 20 C 0.020 3.980 21 N -0.375 5.375 22 C -0.368 4.368 23 H 0.086 0.914 24 C -0.213 4.213 25 N -0.580 5.580 26 Si 0.737 3.263 27 H -0.219 1.219 28 H -0.218 1.218 29 H -0.209 1.209 30 H 0.077 0.923 31 H 0.092 0.908 32 H 0.077 0.923 33 H 0.228 0.772 34 H 0.226 0.774 35 H 0.097 0.903 36 H 0.096 0.904 37 H 0.113 0.887 38 H 0.114 0.886 39 H 0.087 0.913 40 H 0.099 0.901 41 H 0.090 0.910 42 H 0.089 0.911 43 H 0.082 0.918 44 H 0.081 0.919 45 H 0.087 0.913 46 H 0.068 0.932 47 H 0.070 0.930 48 H 0.091 0.909 49 H 0.110 0.890 50 H 0.028 0.972 51 H 0.033 0.967 52 H 0.219 0.781 53 H 0.064 0.936 Dipole moment (debyes) X Y Z Total from point charges 1.943 -30.498 1.826 30.615 hybrid contribution 0.686 -0.063 -0.150 0.705 sum 2.629 -30.561 1.675 30.720 Atomic orbital electron populations 1.21392 0.78113 1.01822 1.02487 1.44914 1.07795 1.07019 1.81225 1.15958 0.81742 0.84323 0.78095 1.91000 1.67688 1.30216 1.56788 1.45082 1.04771 1.09689 1.80799 1.20924 0.87851 0.84023 1.04511 1.21872 0.95565 0.97533 1.04375 1.20240 0.82135 0.90744 0.76970 1.90917 1.67868 1.31812 1.47544 1.48114 1.05157 1.07774 1.73233 1.21810 0.94925 0.86096 0.98732 1.22039 0.98041 0.98430 0.99024 1.21562 0.97551 0.94279 1.01892 1.21567 0.96750 0.98501 0.95578 0.89419 1.21878 1.00114 1.02461 0.95872 1.21333 0.92720 0.83054 0.97401 0.91242 1.22300 0.96911 0.91885 1.04551 1.20791 0.92792 0.87506 0.96874 1.42359 1.05168 1.02178 1.87760 1.19227 0.91884 0.84300 1.41357 0.91394 1.24755 0.95321 0.97859 1.03347 1.69687 1.00411 1.30491 1.57401 0.86211 0.77507 0.85023 0.77591 1.21858 1.21835 1.20888 0.92322 0.90764 0.92296 0.77176 0.77393 0.90339 0.90388 0.88672 0.88590 0.91339 0.90056 0.91023 0.91147 0.91762 0.91906 0.91336 0.93227 0.93034 0.90933 0.88954 0.97177 0.96737 0.78056 0.93606 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.11 0.70 10.70 59.85 0.64 1.34 16 2 N -0.75 -5.46 5.72 -72.07 -0.41 -5.87 16 3 C 0.70 8.08 8.57 -86.92 -0.74 7.34 16 4 O -0.58 -10.00 16.34 5.29 0.09 -9.92 16 5 N -0.75 -6.38 5.56 -65.23 -0.36 -6.74 16 6 C 0.15 1.69 4.96 -4.04 -0.02 1.67 16 7 C -0.15 -1.23 4.17 -27.88 -0.12 -1.34 16 8 C 0.51 6.50 7.53 -10.99 -0.08 6.42 16 9 O -0.51 -9.07 11.53 5.55 0.06 -9.00 16 10 N -0.61 -6.72 2.41 -164.96 -0.40 -7.12 16 11 C 0.11 0.74 6.38 -3.49 -0.02 0.72 16 12 C -0.14 -0.92 5.75 -26.39 -0.15 -1.07 16 13 C -0.12 -0.99 5.43 -26.69 -0.14 -1.13 16 14 C -0.10 -1.24 2.25 -90.50 -0.20 -1.45 16 15 H 0.09 1.14 6.86 -51.93 -0.36 0.79 16 16 C -0.15 -1.84 8.19 37.16 0.30 -1.54 16 17 C 0.16 2.15 2.74 -67.60 -0.18 1.96 16 18 H 0.07 0.94 8.14 -51.93 -0.42 0.52 16 19 C -0.12 -2.21 2.76 -26.73 -0.07 -2.29 16 20 C 0.15 3.25 5.69 -4.04 -0.02 3.22 16 21 N -0.73 -20.36 5.60 -59.90 -0.34 -20.69 16 22 C -0.24 -7.08 9.99 -15.06 -0.15 -7.23 16 23 H 0.07 1.96 7.83 -52.49 -0.41 1.55 16 24 C -0.02 -0.52 5.50 -17.43 -0.10 -0.62 16 25 N -0.93 -30.20 13.06 55.49 0.72 -29.47 16 26 Si 0.91 23.61 29.55 -169.99 -5.02 18.58 16 27 H -0.29 -7.51 7.11 56.52 0.40 -7.11 16 28 H -0.29 -7.64 7.11 56.52 0.40 -7.24 16 29 H -0.28 -7.14 7.11 56.52 0.40 -6.74 16 30 H 0.06 0.46 8.14 -51.93 -0.42 0.03 16 31 H 0.07 0.23 8.14 -51.93 -0.42 -0.19 16 32 H 0.06 0.46 8.14 -51.93 -0.42 0.04 16 33 H 0.40 1.38 8.84 -40.82 -0.36 1.02 16 34 H 0.39 1.72 8.49 -40.82 -0.35 1.37 16 35 H 0.08 1.10 8.10 -51.93 -0.42 0.68 16 36 H 0.08 1.07 8.10 -51.93 -0.42 0.65 16 37 H 0.09 0.52 7.71 -51.93 -0.40 0.12 16 38 H 0.10 0.55 7.92 -51.93 -0.41 0.14 16 39 H 0.07 0.43 8.14 -51.93 -0.42 0.01 16 40 H 0.08 0.36 5.93 -51.93 -0.31 0.05 16 41 H 0.07 0.53 6.33 -51.93 -0.33 0.20 16 42 H 0.07 0.34 8.14 -51.93 -0.42 -0.08 16 43 H 0.06 0.52 8.14 -51.93 -0.42 0.10 16 44 H 0.06 0.51 8.14 -51.93 -0.42 0.09 16 45 H 0.07 1.08 6.15 -51.93 -0.32 0.76 16 46 H 0.05 0.56 8.14 -51.93 -0.42 0.14 16 47 H 0.05 0.60 5.97 -51.93 -0.31 0.29 16 48 H 0.07 1.46 5.52 -51.93 -0.29 1.18 16 49 H 0.09 1.98 6.24 -51.93 -0.32 1.65 16 50 H 0.01 0.23 8.14 -51.93 -0.42 -0.20 16 51 H 0.01 0.32 6.63 -51.93 -0.34 -0.03 16 52 H 0.39 11.57 7.67 -40.82 -0.31 11.26 16 53 H 0.25 7.93 8.31 -40.82 -0.34 7.59 16 LS Contribution 405.73 15.07 6.11 6.11 Total: -1.00 -39.85 405.73 -9.64 -49.48 By element: Atomic # 1 Polarization: 17.66 SS G_CDS: -9.02 Total: 8.63 kcal Atomic # 6 Polarization: 7.07 SS G_CDS: -1.07 Total: 6.01 kcal Atomic # 7 Polarization: -69.12 SS G_CDS: -0.78 Total: -69.90 kcal Atomic # 8 Polarization: -19.07 SS G_CDS: 0.15 Total: -18.92 kcal Atomic # 14 Polarization: 23.61 SS G_CDS: -5.02 Total: 18.58 kcal Total LS contribution 6.11 Total: 6.11 kcal Total: -39.85 -9.64 -49.48 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. ZINC001707149062.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 5.878 kcal (2) G-P(sol) polarization free energy of solvation -39.849 kcal (3) G-ENP(sol) elect.-nuc.-pol. free energy of system -33.970 kcal (4) G-CDS(sol) cavity-dispersion-solvent structure free energy -9.635 kcal (5) G-P-CDS(sol) = G-P(sol) + G-CDS(sol) = (2) + (4) -49.484 kcal (6) G-S(sol) free energy of system = (1) + (5) -43.606 kcal 1SCF WAS SPECIFIED, SO GEOMETRY OPTIMIZATION WAS NOT USED Total computer time = 1.51 seconds