Wall clock time and date at job start Sat Apr 11 2020 06:10:39 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.52998 * 1 3 3 N 1.46505 * 109.46994 * 2 1 4 4 C 1.46503 * 119.99696 * 90.00032 * 3 2 1 5 5 C 1.50703 * 109.46875 * 90.00020 * 4 3 2 6 6 H 1.07992 * 120.00379 * 0.02562 * 5 4 3 7 7 C 1.37880 * 119.99512 * 179.97438 * 5 4 3 8 8 N 1.31514 * 120.00366 * 359.97438 * 7 5 4 9 9 Si 1.86806 * 119.99595 * 180.02562 * 7 5 4 10 10 H 1.48503 * 109.47105 * 60.00173 * 9 7 5 11 11 H 1.48503 * 109.46640 * 179.97438 * 9 7 5 12 12 H 1.48495 * 109.46820 * 299.99789 * 9 7 5 13 13 C 1.34774 * 120.00064 * 270.00077 * 3 2 1 14 14 O 1.21284 * 119.99704 * 179.97438 * 13 3 2 15 15 C 1.50704 * 120.00055 * 0.02562 * 13 3 2 16 16 N 1.46896 * 109.46752 * 173.74750 * 15 13 3 17 17 C 1.46848 * 110.99455 * 169.40141 * 16 15 13 18 18 C 1.53042 * 109.52153 * 185.84205 * 17 16 15 19 19 C 1.53126 * 109.27514 * 300.81944 * 18 17 16 20 20 C 1.50703 * 109.57230 * 177.59207 * 19 18 17 21 21 O 1.20820 * 120.00108 * 240.28341 * 20 19 18 22 22 O 1.34230 * 119.99870 * 60.28604 * 20 19 18 23 23 C 1.45206 * 116.99976 * 179.97438 * 22 20 19 24 24 C 1.52998 * 109.47119 * 149.99649 * 23 22 20 25 25 C 1.52997 * 109.46988 * 270.00339 * 23 22 20 26 26 C 1.53118 * 109.16066 * 57.65495 * 19 18 17 27 27 C 1.46846 * 111.00425 * 293.60905 * 16 15 13 28 28 H 1.09002 * 109.47041 * 300.00009 * 1 2 3 29 29 H 1.09000 * 109.47470 * 59.99758 * 1 2 3 30 30 H 1.08998 * 109.47064 * 180.02562 * 1 2 3 31 31 H 1.08997 * 109.47375 * 240.00307 * 2 1 3 32 32 H 1.09004 * 109.47026 * 120.00010 * 2 1 3 33 33 H 1.09000 * 109.47085 * 209.99628 * 4 3 2 34 34 H 1.09001 * 109.47143 * 330.00067 * 4 3 2 35 35 H 0.96999 * 120.00026 * 180.02562 * 8 7 5 36 36 H 1.08995 * 109.46787 * 293.75180 * 15 13 3 37 37 H 1.09000 * 109.47013 * 53.74864 * 15 13 3 38 38 H 1.09000 * 109.46501 * 305.85172 * 17 16 15 39 39 H 1.09001 * 109.46489 * 65.82722 * 17 16 15 40 40 H 1.09001 * 109.50491 * 180.87976 * 18 17 16 41 41 H 1.08995 * 109.50563 * 60.75990 * 18 17 16 42 42 H 1.09004 * 109.51903 * 297.74941 * 19 18 17 43 43 H 1.09003 * 109.46805 * 29.99895 * 23 22 20 44 44 H 1.09003 * 109.47081 * 60.00415 * 24 23 22 45 45 H 1.08996 * 109.47127 * 179.97438 * 24 23 22 46 46 H 1.08994 * 109.47709 * 300.00349 * 24 23 22 47 47 H 1.08991 * 109.47596 * 59.99598 * 25 23 22 48 48 H 1.08996 * 109.47061 * 180.02562 * 25 23 22 49 49 H 1.09003 * 109.46986 * 299.99446 * 25 23 22 50 50 H 1.08999 * 109.49963 * 182.40006 * 26 19 18 51 51 H 1.09000 * 109.49803 * 62.28812 * 26 19 18 52 52 H 1.09004 * 109.46151 * 294.18036 * 27 16 15 53 53 H 1.09000 * 109.46097 * 54.15336 * 27 16 15 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.5300 0.0000 0.0000 3 7 2.0183 1.3813 0.0000 4 6 2.2624 2.0718 1.2688 5 6 3.6858 1.8350 1.7038 6 1 4.3428 1.2310 1.0958 7 6 4.1493 2.3905 2.8775 8 7 3.3496 3.1273 3.6172 9 14 5.9138 2.0975 3.4164 10 1 6.8409 2.6209 2.3811 11 1 6.1615 2.7984 4.7019 12 1 6.1416 0.6410 3.5942 13 6 2.2429 2.0166 -1.1672 14 8 2.6467 3.1603 -1.1671 15 6 1.9912 1.3064 -2.4723 16 7 2.1649 2.2501 -3.5846 17 6 2.1763 1.5473 -4.8739 18 6 2.5035 2.5380 -5.9936 19 6 1.4518 3.6509 -6.0070 20 6 1.8004 4.6621 -7.0687 21 8 2.0076 5.8124 -6.7626 22 8 1.8806 4.2848 -8.3543 23 6 2.2203 5.3119 -9.3229 24 6 2.9747 4.6792 -10.4939 25 6 0.9392 5.9710 -9.8378 26 6 1.4266 4.3350 -4.6373 27 6 1.1308 3.2923 -3.5571 28 1 -0.3633 0.5138 0.8900 29 1 -0.3634 0.5139 -0.8899 30 1 -0.3633 -1.0276 0.0005 31 1 1.8933 -0.5138 -0.8900 32 1 1.8933 -0.5139 0.8900 33 1 2.0955 3.1411 1.1385 34 1 1.5816 1.6877 2.0284 35 1 3.6758 3.5184 4.4427 36 1 2.6976 0.4835 -2.5802 37 1 0.9737 0.9156 -2.4829 38 1 2.9312 0.7613 -4.8525 39 1 1.1967 1.1049 -5.0552 40 1 2.4965 2.0192 -6.9522 41 1 3.4888 2.9705 -5.8204 42 1 0.4718 3.2234 -6.2188 43 1 2.8506 6.0642 -8.8486 44 1 2.3445 3.9268 -10.9682 45 1 3.2293 5.4501 -11.2211 46 1 3.8874 4.2096 -10.1270 47 1 0.4017 6.4217 -9.0036 48 1 1.1941 6.7417 -10.5652 49 1 0.3089 5.2186 -10.3121 50 1 0.6503 5.1001 -4.6270 51 1 2.3951 4.7958 -4.4429 52 1 0.1560 2.8417 -3.7445 53 1 1.1263 3.7737 -2.5792 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 ZINC000886723607.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:10:39 Heat of formation + Delta-G solvation = -128.574804 kcal Electronic energy + Delta-G solvation = -30648.412135 eV Core-core repulsion = 26519.874455 eV Total energy + Delta-G solvation = -4128.537681 eV No. of doubly occupied orbitals = 66 Molecular weight (most abundant/longest-lived isotopes) = 340.220 amu Computer time = 0.62 seconds Orbital eigenvalues (eV) -40.96003 -39.54900 -38.19469 -36.94805 -35.22179 -33.92076 -33.21992 -32.69248 -30.90561 -30.07357 -28.57091 -27.69307 -26.66273 -24.64176 -23.33180 -22.99909 -22.50264 -21.89501 -20.82291 -19.80734 -18.27879 -17.93314 -16.78978 -16.62295 -16.37472 -16.14750 -15.81345 -15.20779 -14.99985 -14.56638 -14.43878 -14.36878 -14.22024 -13.75731 -13.62379 -13.44980 -13.34113 -13.04135 -12.84095 -12.69818 -12.62152 -12.46748 -12.21056 -12.08375 -11.96921 -11.75461 -11.52102 -11.38709 -11.28946 -11.22225 -11.15805 -10.89976 -10.83320 -10.82569 -10.70310 -10.52213 -10.49280 -10.30557 -10.08827 -9.59608 -8.83837 -8.39412 -8.15159 -7.59704 -6.00432 -4.04808 2.05668 2.43154 3.34611 3.39219 3.41484 3.50389 3.63694 4.19493 4.33642 4.40856 4.48609 4.50963 4.56082 4.62461 4.76341 4.94942 5.01854 5.11980 5.16077 5.24773 5.26512 5.30689 5.33800 5.35490 5.36137 5.41644 5.45297 5.51239 5.54055 5.61240 5.63865 5.66906 5.73737 5.75629 5.79799 5.84160 5.98114 6.02489 6.16393 6.23846 6.26838 6.38151 6.47374 6.95327 7.15210 7.19198 7.48977 7.60605 7.61057 7.96513 8.28077 8.48648 8.95291 9.18628 9.60243 11.07529 Molecular weight = 340.22amu Principal moments of inertia in cm(-1) A = 0.017230 B = 0.002207 C = 0.002089 Principal moments of inertia in units of 10**(-40)*gram-cm**2 A = 1624.718151 B =12686.085316 C =13401.126817 Net atomic charges, atomic populations, and dipole contributions calculated with CM2 Atom NO. Type Charge No. of electrons 1 C -0.169 4.169 2 C 0.107 3.893 3 N -0.601 5.601 4 C 0.248 3.752 5 C -0.524 4.524 6 H 0.056 0.944 7 C 0.060 3.940 8 N -0.896 5.896 9 Si 0.897 3.103 10 H -0.276 1.276 11 H -0.287 1.287 12 H -0.277 1.277 13 C 0.505 3.495 14 O -0.521 6.521 15 C 0.038 3.962 16 N -0.521 5.521 17 C 0.059 3.941 18 C -0.105 4.105 19 C -0.080 4.080 20 C 0.467 3.533 21 O -0.503 6.503 22 O -0.354 6.354 23 C 0.077 3.923 24 C -0.149 4.149 25 C -0.169 4.169 26 C -0.097 4.097 27 C 0.056 3.944 28 H 0.061 0.939 29 H 0.048 0.952 30 H 0.059 0.941 31 H 0.062 0.938 32 H 0.082 0.918 33 H 0.042 0.958 34 H 0.030 0.970 35 H 0.263 0.737 36 H 0.098 0.902 37 H 0.055 0.945 38 H 0.081 0.919 39 H 0.035 0.965 40 H 0.079 0.921 41 H 0.079 0.921 42 H 0.104 0.896 43 H 0.085 0.915 44 H 0.067 0.933 45 H 0.077 0.923 46 H 0.065 0.935 47 H 0.061 0.939 48 H 0.074 0.926 49 H 0.065 0.935 50 H 0.074 0.926 51 H 0.085 0.915 52 H 0.028 0.972 53 H 0.100 0.900 For the charges calculated by CM2: Dipole moment (debyes) X Y Z Total from point charges -7.619 -4.266 -25.728 27.169 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.227 4.227 2 C -0.017 4.017 3 N -0.334 5.334 4 C 0.129 3.871 5 C -0.563 4.563 6 H 0.075 0.925 7 C -0.142 4.142 8 N -0.535 5.535 9 Si 0.725 3.275 10 H -0.202 1.202 11 H -0.212 1.212 12 H -0.203 1.203 13 C 0.290 3.710 14 O -0.397 6.397 15 C -0.093 4.093 16 N -0.245 5.245 17 C -0.069 4.069 18 C -0.143 4.143 19 C -0.099 4.099 20 C 0.309 3.691 21 O -0.390 6.390 22 O -0.271 6.271 23 C 0.023 3.977 24 C -0.206 4.206 25 C -0.227 4.227 26 C -0.135 4.135 27 C -0.071 4.071 28 H 0.080 0.920 29 H 0.067 0.933 30 H 0.078 0.922 31 H 0.081 0.919 32 H 0.100 0.900 33 H 0.060 0.940 34 H 0.048 0.952 35 H 0.073 0.927 36 H 0.116 0.884 37 H 0.073 0.927 38 H 0.100 0.900 39 H 0.053 0.947 40 H 0.098 0.902 41 H 0.098 0.902 42 H 0.122 0.878 43 H 0.103 0.897 44 H 0.086 0.914 45 H 0.095 0.905 46 H 0.084 0.916 47 H 0.080 0.920 48 H 0.093 0.907 49 H 0.084 0.916 50 H 0.092 0.908 51 H 0.104 0.896 52 H 0.046 0.954 53 H 0.118 0.882 Dipole moment (debyes) X Y Z Total from point charges -6.882 -2.934 -25.417 26.495 hybrid contribution 0.174 -0.499 0.270 0.593 sum -6.708 -3.433 -25.147 26.251 Atomic orbital electron populations 1.22225 0.96967 1.01307 1.02208 1.21384 0.93347 0.84425 1.02557 1.48278 1.64166 1.13203 1.07730 1.19229 0.96156 0.93587 0.78098 1.21403 0.96957 1.31830 1.06141 0.92548 1.24936 0.97794 0.94970 0.96452 1.69921 1.35712 1.33428 1.14414 0.86570 0.84941 0.78173 0.77797 1.20165 1.21239 1.20306 1.20142 0.78952 0.87200 0.84692 1.90668 1.44437 1.17784 1.86779 1.22065 1.01086 0.95309 0.90841 1.61520 1.48997 1.12829 1.01199 1.22079 1.00378 0.95798 0.88679 1.21573 0.98747 0.95204 0.98775 1.21231 1.01744 0.94715 0.92249 1.23846 0.73156 0.89199 0.82868 1.90664 1.45217 1.20885 1.82195 1.86575 1.80737 1.42421 1.17415 1.22726 0.97166 0.89711 0.88115 1.21882 1.00492 1.01828 0.96417 1.22203 0.97417 1.01065 1.01997 1.21490 1.04480 0.97339 0.90219 1.22199 0.91051 0.92740 1.01143 0.92018 0.93283 0.92249 0.91949 0.89970 0.93971 0.95162 0.92723 0.88415 0.92669 0.90041 0.94695 0.90199 0.90246 0.87847 0.89733 0.91413 0.90452 0.91587 0.91968 0.90741 0.91568 0.90784 0.89602 0.95435 0.88193 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 C -0.17 -2.18 10.18 37.16 0.38 -1.81 16 2 C 0.11 1.67 5.85 -4.04 -0.02 1.64 16 3 N -0.60 -12.46 2.46 -175.05 -0.43 -12.89 16 4 C 0.25 6.47 5.16 -5.19 -0.03 6.44 16 5 C -0.52 -14.93 9.39 -34.44 -0.32 -15.26 16 6 H 0.06 1.61 7.81 -52.49 -0.41 1.20 16 7 C 0.06 1.74 5.46 -17.82 -0.10 1.64 16 8 N -0.90 -26.82 13.29 55.43 0.74 -26.09 16 9 Si 0.90 22.12 29.54 -169.99 -5.02 17.10 16 10 H -0.28 -6.87 7.11 56.52 0.40 -6.47 16 11 H -0.29 -6.96 7.11 56.52 0.40 -6.56 16 12 H -0.28 -6.61 7.11 56.52 0.40 -6.21 16 13 C 0.50 10.32 7.20 -10.98 -0.08 10.24 16 14 O -0.52 -13.04 14.29 5.55 0.08 -12.96 16 15 C 0.04 0.56 3.92 -4.98 -0.02 0.54 16 16 N -0.52 -7.44 4.35 -236.12 -1.03 -8.47 16 17 C 0.06 0.61 5.46 -3.83 -0.02 0.59 16 18 C -0.11 -1.03 5.32 -26.64 -0.14 -1.18 16 19 C -0.08 -0.85 2.69 -91.64 -0.25 -1.10 16 20 C 0.47 5.57 7.23 36.01 0.26 5.83 16 21 O -0.50 -7.47 15.20 -18.74 -0.28 -7.76 16 22 O -0.35 -3.66 9.82 -39.26 -0.39 -4.05 16 23 C 0.08 0.67 3.18 -26.73 -0.09 0.58 16 24 C -0.15 -0.90 9.57 37.16 0.36 -0.54 16 25 C -0.17 -1.24 9.55 37.16 0.35 -0.88 16 26 C -0.10 -1.27 4.68 -26.65 -0.12 -1.40 16 27 C 0.06 0.79 5.12 -3.84 -0.02 0.77 16 28 H 0.06 0.90 8.14 -51.93 -0.42 0.47 16 29 H 0.05 0.60 7.12 -51.93 -0.37 0.23 16 30 H 0.06 0.61 8.14 -51.93 -0.42 0.19 16 31 H 0.06 0.79 6.61 -51.93 -0.34 0.45 16 32 H 0.08 1.34 8.07 -51.93 -0.42 0.92 16 33 H 0.04 1.25 7.42 -51.93 -0.39 0.86 16 34 H 0.03 0.80 8.07 -51.93 -0.42 0.38 16 35 H 0.26 7.47 8.31 -40.82 -0.34 7.13 16 36 H 0.10 1.28 7.55 -51.93 -0.39 0.89 16 37 H 0.06 0.69 6.83 -51.93 -0.35 0.33 16 38 H 0.08 0.76 8.01 -51.93 -0.42 0.35 16 39 H 0.03 0.31 8.14 -51.93 -0.42 -0.11 16 40 H 0.08 0.63 8.14 -51.93 -0.42 0.21 16 41 H 0.08 0.90 8.14 -51.93 -0.42 0.48 16 42 H 0.10 0.90 8.14 -51.93 -0.42 0.48 16 43 H 0.08 0.84 7.16 -51.93 -0.37 0.47 16 44 H 0.07 0.37 8.14 -51.93 -0.42 -0.05 16 45 H 0.08 0.34 8.14 -51.93 -0.42 -0.08 16 46 H 0.07 0.41 8.14 -51.93 -0.42 -0.01 16 47 H 0.06 0.54 8.14 -51.93 -0.42 0.12 16 48 H 0.07 0.41 8.14 -51.93 -0.42 -0.01 16 49 H 0.07 0.43 8.14 -51.93 -0.42 0.01 16 50 H 0.07 0.93 8.14 -51.93 -0.42 0.50 16 51 H 0.09 1.34 8.04 -51.93 -0.42 0.92 16 52 H 0.03 0.33 8.14 -51.93 -0.42 -0.09 16 53 H 0.10 1.75 6.05 -51.93 -0.31 1.44 16 LS Contribution 421.24 15.07 6.35 6.35 Total: -1.00 -34.70 421.24 -9.51 -44.21 By element: Atomic # 1 Polarization: 8.10 SS G_CDS: -9.67 Total: -1.56 kcal Atomic # 6 Polarization: 5.98 SS G_CDS: 0.14 Total: 6.12 kcal Atomic # 7 Polarization: -46.73 SS G_CDS: -0.72 Total: -47.45 kcal Atomic # 8 Polarization: -24.17 SS G_CDS: -0.59 Total: -24.76 kcal Atomic # 14 Polarization: 22.12 SS G_CDS: -5.02 Total: 17.10 kcal Total LS contribution 6.35 Total: 6.35 kcal Total: -34.70 -9.51 -44.21 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. ZINC000886723607.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 -84.366 kcal (2) G-P(sol) polarization free energy of solvation -34.699 kcal (3) G-ENP(sol) elect.-nuc.-pol. free energy of system -119.065 kcal (4) G-CDS(sol) cavity-dispersion-solvent structure free energy -9.509 kcal (5) G-P-CDS(sol) = G-P(sol) + G-CDS(sol) = (2) + (4) -44.209 kcal (6) G-S(sol) free energy of system = (1) + (5) -128.575 kcal 1SCF WAS SPECIFIED, SO GEOMETRY OPTIMIZATION WAS NOT USED Total computer time = 0.62 seconds