Wall clock time and date at job start Tue Mar 31 2020 05:48:57 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 * SM5.42R - SM5.42R CALCULATIONS WILL BE PERFORMED * - CHARGE MODEL 2 WILL BE USED (CM2A) * - THE SOLVENT IS WATER ******************************************************************************* 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.52997 * 1 3 3 H 1.08999 * 109.46987 * 2 1 4 4 C 1.50692 * 109.47267 * 120.00602 * 2 1 3 5 5 H 1.08003 * 120.00316 * 114.36808 * 4 2 1 6 6 C 1.37875 * 120.00133 * 294.36950 * 4 2 1 7 7 N 1.31516 * 120.00510 * 156.90156 * 6 4 2 8 8 Si 1.86810 * 119.99756 * 337.17312 * 6 4 2 9 9 H 1.48501 * 109.47061 * 179.72434 * 8 6 4 10 10 H 1.48495 * 109.46894 * 299.72787 * 8 6 4 11 11 H 1.48496 * 109.46696 * 59.72479 * 8 6 4 12 12 N 1.46503 * 109.46749 * 240.00646 * 2 1 3 13 13 C 1.46505 * 120.00016 * 299.99254 * 12 2 1 14 14 C 1.34773 * 120.00310 * 119.99959 * 12 2 1 15 15 O 1.21276 * 119.99959 * 359.97438 * 14 12 2 16 16 C 1.50704 * 120.00250 * 179.97438 * 14 12 2 17 17 C 1.52995 * 109.46974 * 179.97438 * 16 14 12 18 18 C 1.53008 * 109.46819 * 180.02562 * 17 16 14 19 19 N 1.46903 * 109.46813 * 180.02562 * 18 17 16 20 20 C 1.46955 * 111.00022 * 293.80281 * 19 18 17 21 21 C 1.53059 * 109.36296 * 176.81386 * 20 19 18 22 22 N 1.46952 * 109.36407 * 58.39793 * 21 20 19 23 23 C 1.39675 * 111.00036 * 176.80882 * 22 21 20 24 24 C 1.38916 * 120.07307 * 123.52450 * 23 22 21 25 25 C 1.38098 * 119.92404 * 179.77003 * 24 23 22 26 26 C 1.38289 * 120.07325 * 0.51234 * 25 24 23 27 27 C 1.38282 * 120.15299 * 359.74042 * 26 25 24 28 28 C 1.38109 * 120.07050 * 0.02562 * 27 26 25 29 29 C 1.46956 * 110.88401 * 300.68375 * 22 21 20 30 30 C 1.46952 * 110.99674 * 170.00075 * 19 18 17 31 31 H 1.09005 * 109.47375 * 53.29697 * 1 2 3 32 32 H 1.09000 * 109.47483 * 173.29600 * 1 2 3 33 33 H 1.09002 * 109.47176 * 293.29860 * 1 2 3 34 34 H 0.96989 * 119.99527 * 355.20009 * 7 6 4 35 35 H 1.09007 * 109.46817 * 90.00493 * 13 12 2 36 36 H 1.08992 * 109.47593 * 210.00355 * 13 12 2 37 37 H 1.08998 * 109.47142 * 330.00868 * 13 12 2 38 38 H 1.08998 * 109.47200 * 300.00174 * 16 14 12 39 39 H 1.09003 * 109.46884 * 60.00320 * 16 14 12 40 40 H 1.08998 * 109.47315 * 300.00122 * 17 16 14 41 41 H 1.09003 * 109.47038 * 60.00433 * 17 16 14 42 42 H 1.08998 * 109.47275 * 300.00144 * 18 17 16 43 43 H 1.08995 * 109.47046 * 60.00710 * 18 17 16 44 44 H 1.09003 * 109.48818 * 296.77386 * 20 19 18 45 45 H 1.08997 * 109.49467 * 56.84835 * 20 19 18 46 46 H 1.08997 * 109.48922 * 298.43368 * 21 20 19 47 47 H 1.08996 * 109.48757 * 178.36700 * 21 20 19 48 48 H 1.07996 * 120.03624 * 0.02562 * 24 23 22 49 49 H 1.07998 * 119.96315 * 180.25019 * 25 24 23 50 50 H 1.08008 * 119.92323 * 179.76718 * 26 25 24 51 51 H 1.08000 * 119.96620 * 180.02562 * 27 26 25 52 52 H 1.07999 * 120.03366 * 179.97438 * 28 27 26 53 53 H 1.09003 * 109.48235 * 299.34829 * 29 22 21 54 54 H 1.08999 * 109.48903 * 179.29022 * 29 22 21 55 55 H 1.08996 * 109.49135 * 303.15553 * 30 19 18 56 56 H 1.09005 * 109.49153 * 63.22672 * 30 19 18 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 1 1.8933 1.0277 0.0000 4 6 2.0323 -0.7105 1.2303 5 1 2.5762 -1.6386 1.1338 6 6 1.7976 -0.1758 2.4793 7 7 1.8082 -0.9532 3.5401 8 14 1.4571 1.6517 2.6640 9 1 1.2325 1.9710 4.0968 10 1 2.6206 2.4252 2.1609 11 1 0.2477 2.0114 1.8811 12 7 2.0182 -0.6905 -1.1963 13 6 1.6644 -2.0931 -1.4282 14 6 2.7930 -0.0355 -2.0834 15 8 3.0857 1.1257 -1.8917 16 6 3.2957 -0.7460 -3.3138 17 6 4.1451 0.2160 -4.1469 18 6 4.6549 -0.5052 -5.3963 19 7 5.4709 0.4182 -6.1959 20 6 4.6515 1.5115 -6.7370 21 6 5.5499 2.4918 -7.4951 22 7 6.2491 1.7783 -8.5728 23 6 7.0250 2.6563 -9.3332 24 6 8.3871 2.4291 -9.4844 25 6 9.1504 3.2965 -10.2408 26 6 8.5622 4.3961 -10.8386 27 6 7.2075 4.6282 -10.6858 28 6 6.4375 3.7620 -9.9347 29 6 7.0685 0.6850 -8.0317 30 6 6.1702 -0.2952 -7.2737 31 1 -0.3634 0.6142 -0.8239 32 1 -0.3634 -1.0206 -0.1200 33 1 -0.3634 0.4065 0.9439 34 1 1.9116 -1.9121 3.4375 35 1 0.7374 -2.1443 -1.9995 36 1 2.4627 -2.5818 -1.9866 37 1 1.5295 -2.5968 -0.4710 38 1 3.9014 -1.6022 -3.0169 39 1 2.4481 -1.0888 -3.9072 40 1 3.5390 1.0719 -4.4437 41 1 4.9927 0.5592 -3.5537 42 1 5.2609 -1.3614 -5.0998 43 1 3.8071 -0.8480 -5.9893 44 1 3.9037 1.1033 -7.4170 45 1 4.1540 2.0324 -5.9190 46 1 6.2805 2.9205 -6.8091 47 1 4.9403 3.2882 -7.9219 48 1 8.8478 1.5729 -9.0144 49 1 10.2086 3.1183 -10.3621 50 1 9.1622 5.0750 -11.4265 51 1 6.7512 5.4876 -11.1543 52 1 5.3796 3.9440 -9.8154 53 1 7.8163 1.0932 -7.3519 54 1 7.5661 0.1642 -8.8498 55 1 6.7798 -1.0916 -6.8468 56 1 5.4396 -0.7240 -7.9598 RHF calculation, no. of doubly occupied orbitals= 69 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=WATER ZINC000071909759.mol2 56 1SCF WAS SPECIFIED, SO GEOMETRY OPTIMIZATION WAS NOT USED SCF FIELD WAS ACHIEVED Wall clock time and date at job start Tue Mar 31 2020 05:48:57 Heat of formation + Delta-G solvation = -13.988462 kcal Electronic energy + Delta-G solvation = -31696.574904 eV Core-core repulsion = 27604.728493 eV Total energy + Delta-G solvation = -4091.846411 eV No. of doubly occupied orbitals = 69 Molecular weight (most abundant/longest-lived isotopes) = 359.231 amu Computer time = 1.46 seconds Orbital eigenvalues (eV) -41.38246 -40.82887 -39.33911 -37.37029 -36.78072 -34.91749 -34.13785 -32.62693 -32.09642 -31.41392 -30.70815 -28.26413 -27.62128 -26.83039 -26.02184 -25.18536 -23.35191 -22.73588 -22.46311 -21.66146 -20.78603 -19.91722 -19.52500 -18.70141 -17.92774 -17.11413 -17.02744 -16.76553 -16.52770 -16.20625 -16.12205 -15.67012 -15.64861 -15.40959 -15.11067 -14.97852 -14.75245 -14.57878 -14.44077 -14.36260 -14.06570 -13.72136 -13.51115 -13.46735 -13.27698 -13.12772 -13.04319 -12.92128 -12.68005 -12.55832 -12.44457 -12.38415 -12.22323 -12.13378 -11.95540 -11.86446 -11.74844 -11.70908 -11.42403 -11.24208 -10.69466 -10.63358 -10.46742 -9.85948 -9.48287 -9.33339 -8.82205 -8.72507 -6.38235 0.44736 0.48495 1.41297 1.48433 1.68160 1.80036 2.32169 2.41471 2.60842 2.86971 2.96196 2.97146 3.67845 3.74293 3.88292 3.94769 4.03736 4.07688 4.11688 4.21750 4.26159 4.29594 4.31022 4.37084 4.45128 4.47961 4.55425 4.57075 4.58973 4.60742 4.63582 4.67817 4.70758 4.73751 4.75742 4.77482 4.81297 4.84664 4.92348 4.98214 5.00176 5.12202 5.13914 5.20282 5.30215 5.35012 5.40334 5.51020 5.65884 5.73130 5.81611 5.83250 5.94641 6.03823 6.05921 6.09930 6.34247 6.56167 6.77148 7.28190 7.43098 8.57838 Molecular weight = 359.23amu Principal moments of inertia in cm(-1) A = 0.019726 B = 0.001573 C = 0.001520 Principal moments of inertia in units of 10**(-40)*gram-cm**2 A = 1419.139363 B =17797.176188 C =18417.497658 Net atomic charges, atomic populations, and dipole contributions calculated with CM2 Atom NO. Type Charge No. of electrons 1 C -0.145 4.145 2 C 0.232 3.768 3 H 0.027 0.973 4 C -0.612 4.612 5 H 0.034 0.966 6 C -0.001 4.001 7 N -0.923 5.923 8 Si 1.065 2.935 9 H -0.280 1.280 10 H -0.284 1.284 11 H -0.247 1.247 12 N -0.601 5.601 13 C 0.082 3.918 14 C 0.491 3.509 15 O -0.603 6.603 16 C -0.122 4.122 17 C -0.091 4.091 18 C 0.067 3.933 19 N -0.555 5.555 20 C 0.056 3.944 21 C 0.066 3.934 22 N -0.496 5.496 23 C 0.096 3.904 24 C -0.183 4.183 25 C -0.102 4.102 26 C -0.155 4.155 27 C -0.102 4.102 28 C -0.155 4.155 29 C 0.055 3.945 30 C 0.059 3.941 31 H 0.045 0.955 32 H 0.058 0.942 33 H 0.049 0.951 34 H 0.259 0.741 35 H 0.072 0.928 36 H 0.092 0.908 37 H 0.054 0.946 38 H 0.111 0.889 39 H 0.119 0.881 40 H 0.055 0.945 41 H 0.049 0.951 42 H 0.090 0.910 43 H 0.050 0.950 44 H 0.046 0.954 45 H 0.081 0.919 46 H 0.038 0.962 47 H 0.092 0.908 48 H 0.134 0.866 49 H 0.140 0.860 50 H 0.140 0.860 51 H 0.141 0.859 52 H 0.127 0.873 53 H 0.047 0.953 54 H 0.098 0.902 55 H 0.098 0.902 56 H 0.044 0.956 For the charges calculated by CM2: Dipole moment (debyes) X Y Z Total from point charges 9.552 2.526 -32.548 34.014 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.204 4.204 2 C 0.129 3.871 3 H 0.045 0.955 4 C -0.650 4.650 5 H 0.052 0.948 6 C -0.203 4.203 7 N -0.557 5.557 8 Si 0.888 3.112 9 H -0.205 1.205 10 H -0.210 1.210 11 H -0.171 1.171 12 N -0.331 5.331 13 C -0.060 4.060 14 C 0.281 3.719 15 O -0.486 6.486 16 C -0.162 4.162 17 C -0.130 4.130 18 C -0.060 4.060 19 N -0.280 5.280 20 C -0.072 4.072 21 C -0.060 4.060 22 N -0.212 5.212 23 C -0.002 4.002 24 C -0.203 4.203 25 C -0.120 4.120 26 C -0.174 4.174 27 C -0.121 4.121 28 C -0.175 4.175 29 C -0.072 4.072 30 C -0.069 4.069 31 H 0.064 0.936 32 H 0.077 0.923 33 H 0.068 0.932 34 H 0.068 0.932 35 H 0.091 0.909 36 H 0.111 0.889 37 H 0.073 0.927 38 H 0.129 0.871 39 H 0.137 0.863 40 H 0.074 0.926 41 H 0.068 0.932 42 H 0.108 0.892 43 H 0.068 0.932 44 H 0.064 0.936 45 H 0.099 0.901 46 H 0.056 0.944 47 H 0.111 0.889 48 H 0.152 0.848 49 H 0.158 0.842 50 H 0.157 0.843 51 H 0.159 0.841 52 H 0.145 0.855 53 H 0.066 0.934 54 H 0.116 0.884 55 H 0.116 0.884 56 H 0.062 0.938 Dipole moment (debyes) X Y Z Total from point charges 9.350 4.204 -31.250 32.888 hybrid contribution 0.298 -1.482 -0.408 1.565 sum 9.648 2.722 -31.658 33.207 Atomic orbital electron populations 1.22245 0.96585 1.01259 1.00280 1.18881 0.91878 0.93277 0.83099 0.95484 1.22250 1.40207 1.09963 0.92621 0.94764 1.26711 0.93073 1.08541 0.91954 1.70961 1.46254 1.05888 1.32610 0.83765 0.78602 0.72190 0.76633 1.20542 1.21027 1.17114 1.48080 1.49581 1.13343 1.22128 1.21817 1.01663 0.80523 1.02005 1.21570 0.77794 0.86874 0.85614 1.90610 1.59867 1.19738 1.78379 1.21624 1.02101 0.99631 0.92842 1.21203 0.98675 0.97626 0.95489 1.22034 0.96264 0.94678 0.93066 1.62095 1.30493 1.12150 1.23236 1.22195 0.94454 0.92932 0.97643 1.22251 0.96941 0.95984 0.90792 1.59314 1.28699 1.10778 1.22378 1.19169 0.91332 0.92664 0.97001 1.21526 0.92527 1.02607 1.03679 1.21321 1.00608 0.94176 0.95929 1.21408 0.95418 0.98734 1.01848 1.21318 0.94294 0.99091 0.97372 1.20993 1.00491 0.95406 1.00641 1.22376 0.94724 0.91911 0.98183 1.22209 0.97117 0.95817 0.91744 0.93573 0.92304 0.93160 0.93219 0.90914 0.88936 0.92719 0.87138 0.86332 0.92639 0.93192 0.89194 0.93220 0.93616 0.90114 0.94367 0.88941 0.84843 0.84249 0.84250 0.84107 0.85520 0.93427 0.88424 0.88401 0.93752 Number of geometries 1 Number of calculations of the screened coulomb radii 1 The total number of SCF iterations 15. 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.15 -6.16 9.01 71.98 0.65 -5.51 16 2 C 0.23 11.29 2.23 44.24 0.10 11.39 16 3 H 0.03 1.47 5.48 -2.39 -0.01 1.46 16 4 C -0.61 -34.38 8.30 25.60 0.21 -34.16 16 5 H 0.03 2.06 7.38 -2.91 -0.02 2.04 16 6 C 0.00 -0.08 5.20 84.65 0.44 0.36 16 7 N -0.92 -56.51 13.70 21.45 0.29 -56.22 16 8 Si 1.07 51.95 26.44 68.60 1.81 53.77 16 9 H -0.28 -13.39 7.11 99.48 0.71 -12.68 16 10 H -0.28 -13.97 7.11 99.48 0.71 -13.27 16 11 H -0.25 -10.80 6.12 99.48 0.61 -10.19 16 12 N -0.60 -24.14 2.75 -826.81 -2.28 -26.41 16 13 C 0.08 2.48 9.10 127.77 1.16 3.64 16 14 C 0.49 18.87 7.67 87.66 0.67 19.54 16 15 O -0.60 -28.85 15.31 -3.02 -0.05 -28.89 16 16 C -0.12 -2.96 4.25 29.85 0.13 -2.83 16 17 C -0.09 -2.22 4.44 30.60 0.14 -2.09 16 18 C 0.07 1.16 4.96 86.31 0.43 1.59 16 19 N -0.55 -10.36 4.71 -895.70 -4.22 -14.58 16 20 C 0.06 1.05 5.41 86.31 0.47 1.51 16 21 C 0.07 1.20 5.38 86.31 0.46 1.67 16 22 N -0.50 -8.76 4.75 -694.79 -3.30 -12.06 16 23 C 0.10 1.71 5.64 38.21 0.22 1.92 16 24 C -0.18 -2.91 8.82 22.40 0.20 -2.71 16 25 C -0.10 -1.39 10.04 22.25 0.22 -1.17 16 26 C -0.16 -2.01 10.04 22.29 0.22 -1.79 16 27 C -0.10 -1.41 10.04 22.25 0.22 -1.19 16 28 C -0.16 -2.56 9.19 22.40 0.21 -2.35 16 29 C 0.06 0.84 4.29 86.31 0.37 1.21 16 30 C 0.06 0.87 5.61 86.31 0.48 1.36 16 31 H 0.05 1.74 8.14 -2.38 -0.02 1.72 16 32 H 0.06 2.30 7.46 -2.39 -0.02 2.28 16 33 H 0.05 2.22 5.44 -2.39 -0.01 2.20 16 34 H 0.26 15.84 8.73 -92.72 -0.81 15.03 16 35 H 0.07 1.77 8.14 -2.38 -0.02 1.75 16 36 H 0.09 2.09 6.44 -2.39 -0.02 2.07 16 37 H 0.05 1.94 6.92 -2.39 -0.02 1.92 16 38 H 0.11 2.25 7.40 -2.39 -0.02 2.23 16 39 H 0.12 2.14 8.11 -2.39 -0.02 2.12 16 40 H 0.05 1.49 6.45 -2.39 -0.02 1.48 16 41 H 0.05 1.50 8.10 -2.39 -0.02 1.48 16 42 H 0.09 1.34 8.03 -2.39 -0.02 1.32 16 43 H 0.05 0.75 8.12 -2.39 -0.02 0.73 16 44 H 0.05 0.83 8.14 -2.39 -0.02 0.81 16 45 H 0.08 1.67 6.34 -2.39 -0.02 1.66 16 46 H 0.04 0.81 8.14 -2.39 -0.02 0.79 16 47 H 0.09 1.50 6.82 -2.39 -0.02 1.48 16 48 H 0.13 1.80 5.81 -2.91 -0.02 1.78 16 49 H 0.14 1.27 8.06 -2.91 -0.02 1.25 16 50 H 0.14 1.16 8.06 -2.91 -0.02 1.14 16 51 H 0.14 1.31 8.06 -2.91 -0.02 1.29 16 52 H 0.13 1.89 7.18 -2.91 -0.02 1.86 16 53 H 0.05 0.80 7.57 -2.39 -0.02 0.79 16 54 H 0.10 1.22 7.18 -2.39 -0.02 1.20 16 55 H 0.10 1.17 8.07 -2.39 -0.02 1.15 16 56 H 0.04 0.66 8.14 -2.38 -0.02 0.64 16 Total: -1.00 -74.46 425.52 -0.02 -74.48 By element: Atomic # 1 Polarization: 18.81 SS G_CDS: 0.71 Total: 19.52 kcal Atomic # 6 Polarization: -16.61 SS G_CDS: 7.00 Total: -9.60 kcal Atomic # 7 Polarization: -99.77 SS G_CDS: -9.50 Total: -109.27 kcal Atomic # 8 Polarization: -28.85 SS G_CDS: -0.05 Total: -28.89 kcal Atomic # 14 Polarization: 51.95 SS G_CDS: 1.81 Total: 53.77 kcal Total: -74.46 -0.02 -74.48 kcal The number of atoms in the molecule is 56 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. ZINC000071909759.mol2 56 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 60.490 kcal (2) G-P(sol) polarization free energy of solvation -74.459 kcal (3) G-ENP(sol) elect.-nuc.-pol. free energy of system -13.969 kcal (4) G-CDS(sol) cavity-dispersion-solvent structure free energy -0.020 kcal (5) G-P-CDS(sol) = G-P(sol) + G-CDS(sol) = (2) + (4) -74.479 kcal (6) G-S(sol) free energy of system = (1) + (5) -13.988 kcal 1SCF WAS SPECIFIED, SO GEOMETRY OPTIMIZATION WAS NOT USED Total computer time = 1.46 seconds