Wall clock time and date at job start Sat Apr 11 2020 06:12:40 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.53001 * 1 3 3 O 1.42899 * 109.47240 * 2 1 4 4 C 1.42898 * 113.99997 * 184.39619 * 3 2 1 5 5 H 1.09004 * 109.47526 * 8.65231 * 4 3 2 6 6 C 1.50700 * 109.47195 * 128.65211 * 4 3 2 7 7 O 1.21286 * 119.99648 * 3.48789 * 6 4 3 8 8 N 1.34777 * 120.00266 * 183.49629 * 6 4 3 9 9 C 1.46499 * 119.99989 * 174.63760 * 8 6 4 10 10 C 1.53003 * 109.46915 * 179.97438 * 9 8 6 11 11 H 1.09004 * 110.88277 * 54.30223 * 10 9 8 12 12 C 1.55155 * 111.00418 * 178.17928 * 10 9 8 13 13 C 1.54317 * 102.94451 * 153.61540 * 12 10 9 14 14 N 1.47029 * 107.27138 * 337.80980 * 13 12 10 15 15 C 1.36930 * 125.64825 * 178.98738 * 14 13 12 16 16 H 1.08001 * 119.99959 * 179.97438 * 15 14 13 17 17 C 1.38808 * 120.00200 * 359.96992 * 15 14 13 18 18 N 1.31626 * 120.00175 * 0.02562 * 17 15 14 19 19 Si 1.86805 * 120.00261 * 179.97438 * 17 15 14 20 20 H 1.48499 * 109.46752 * 60.00312 * 19 17 15 21 21 H 1.48494 * 109.47461 * 179.97438 * 19 17 15 22 22 H 1.48502 * 109.47239 * 300.00512 * 19 17 15 23 23 C 1.47426 * 108.70042 * 358.93588 * 14 13 12 24 24 C 1.52998 * 109.47385 * 248.64927 * 4 3 2 25 25 C 1.53187 * 109.54177 * 61.27823 * 24 4 3 26 26 C 1.53101 * 109.15755 * 183.24290 * 25 24 4 27 27 O 1.42906 * 109.48814 * 57.47290 * 26 25 24 28 28 C 1.42909 * 114.09642 * 298.88797 * 27 26 25 29 29 C 1.53096 * 109.41488 * 61.16098 * 28 27 26 30 30 H 1.09002 * 109.46637 * 299.99552 * 1 2 3 31 31 H 1.09004 * 109.46978 * 59.99478 * 1 2 3 32 32 H 1.08995 * 109.47256 * 179.97438 * 1 2 3 33 33 H 1.09002 * 109.46637 * 240.00447 * 2 1 3 34 34 H 1.08997 * 109.47353 * 120.00582 * 2 1 3 35 35 H 0.97002 * 119.99554 * 354.63255 * 8 6 4 36 36 H 1.08999 * 109.47259 * 59.99793 * 9 8 6 37 37 H 1.09001 * 109.47388 * 299.99490 * 9 8 6 38 38 H 1.09001 * 110.72270 * 271.97938 * 12 10 9 39 39 H 1.09005 * 110.72165 * 35.07100 * 12 10 9 40 40 H 1.09001 * 109.88125 * 97.22389 * 13 12 10 41 41 H 1.08998 * 109.88536 * 218.39197 * 13 12 10 42 42 H 0.96992 * 119.99870 * 179.97438 * 18 17 15 43 43 H 1.08992 * 110.36771 * 265.29211 * 23 14 13 44 44 H 1.08999 * 110.36121 * 142.96434 * 23 14 13 45 45 H 1.08993 * 109.54824 * 181.45376 * 24 4 3 46 46 H 1.09003 * 109.51781 * 303.05850 * 25 24 4 47 47 H 1.08998 * 109.51885 * 63.24467 * 25 24 4 48 48 H 1.09004 * 109.40840 * 297.49385 * 26 25 24 49 49 H 1.09000 * 109.48441 * 177.50075 * 26 25 24 50 50 H 1.08996 * 109.48457 * 301.17708 * 28 27 26 51 51 H 1.09001 * 109.48962 * 181.14632 * 28 27 26 52 52 H 1.08993 * 109.52335 * 182.48400 * 29 28 27 53 53 H 1.08994 * 109.57551 * 62.32106 * 29 28 27 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 8 2.0064 1.3473 0.0000 4 6 3.4273 1.4613 -0.1001 5 1 3.8565 0.4819 -0.3113 6 6 3.7775 2.4141 -1.2139 7 8 2.8979 2.9917 -1.8168 8 7 5.0682 2.6241 -1.5401 9 6 5.4051 3.4576 -2.6969 10 6 6.9260 3.5340 -2.8448 11 1 7.3625 2.5363 -2.8910 12 6 7.3224 4.3726 -4.0886 13 6 8.7219 4.9048 -3.7150 14 7 8.8264 4.8666 -2.2489 15 6 9.9223 5.2384 -1.5170 16 1 9.8993 5.1734 -0.4392 17 6 11.0630 5.6987 -2.1603 18 7 11.0910 5.7784 -3.4739 19 14 12.5579 6.2066 -1.1618 20 1 13.0378 5.0487 -0.3654 21 1 13.6358 6.6569 -2.0785 22 1 12.1886 7.3171 -0.2477 23 6 7.5553 4.3630 -1.6974 24 6 3.9907 1.9910 1.2201 25 6 3.6631 1.0079 2.3482 26 6 4.1592 1.5844 3.6769 27 8 3.5561 2.8610 3.8976 28 6 3.8679 3.8341 2.8987 29 6 3.3580 3.3492 1.5390 30 1 -0.3633 0.5138 0.8901 31 1 -0.3633 0.5139 -0.8900 32 1 -0.3633 -1.0276 -0.0005 33 1 1.8933 -0.5138 -0.8901 34 1 1.8934 -0.5139 0.8899 35 1 5.7741 2.2235 -1.0089 36 1 5.0037 4.4605 -2.5514 37 1 4.9737 3.0207 -3.5975 38 1 7.3736 3.7437 -4.9774 39 1 6.6219 5.1941 -4.2391 40 1 9.4900 4.2708 -4.1580 41 1 8.8352 5.9297 -4.0682 42 1 11.8879 6.1004 -3.9233 43 1 6.9033 5.1936 -1.4277 44 1 7.7427 3.7295 -0.8304 45 1 5.0716 2.1030 1.1365 46 1 2.5847 0.8570 2.3981 47 1 4.1566 0.0553 2.1559 48 1 5.2430 1.6957 3.6413 49 1 3.8895 0.9096 4.4894 50 1 4.9478 3.9749 2.8528 51 1 3.3884 4.7801 3.1503 52 1 3.6337 4.0692 0.7686 53 1 2.2733 3.2471 1.5710 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 ZINC001714632118.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:12:40 Heat of formation + Delta-G solvation = -117.905039 kcal Electronic energy + Delta-G solvation = -31271.559959 eV Core-core repulsion = 27143.484954 eV Total energy + Delta-G solvation = -4128.075005 eV No. of doubly occupied orbitals = 66 Molecular weight (most abundant/longest-lived isotopes) = 340.220 amu Computer time = 0.87 seconds Orbital eigenvalues (eV) -40.63787 -39.11429 -38.06628 -37.41088 -34.76907 -34.31650 -32.49935 -31.53265 -31.24893 -29.99182 -28.76299 -27.64949 -25.50504 -24.24069 -23.86201 -23.22717 -22.87030 -20.83895 -20.28103 -19.96682 -19.68298 -17.49890 -17.16495 -16.78940 -16.70583 -16.13456 -15.91701 -15.60488 -15.15519 -14.99988 -14.79863 -14.62179 -14.41132 -14.06404 -13.98645 -13.61652 -13.29613 -13.01289 -12.72554 -12.53764 -12.44039 -12.37523 -12.07714 -11.99019 -11.78506 -11.64600 -11.48681 -11.42304 -11.28123 -10.99307 -10.93330 -10.80600 -10.54141 -10.43498 -10.25442 -10.17476 -10.02068 -9.90742 -9.80791 -9.68436 -9.30510 -8.98833 -8.62594 -6.86401 -5.73934 -3.07731 2.56025 2.82959 3.37979 3.42392 3.43773 3.63038 3.67817 4.25819 4.31641 4.46646 4.50564 4.67960 4.70455 4.79917 4.80302 4.98778 5.01527 5.06456 5.09487 5.15591 5.19550 5.25641 5.32442 5.38695 5.40326 5.43157 5.47192 5.51833 5.62440 5.63195 5.65495 5.78530 5.84829 5.86163 5.90940 6.01615 6.02765 6.16212 6.31362 6.40152 6.45357 6.50286 6.64183 6.82130 6.87699 6.98905 7.50283 7.52586 7.60472 7.79714 8.21165 8.23366 9.19235 9.79691 10.45696 11.37739 Molecular weight = 340.22amu Principal moments of inertia in cm(-1) A = 0.013011 B = 0.002732 C = 0.002357 Principal moments of inertia in units of 10**(-40)*gram-cm**2 A = 2151.541424 B =10246.591091 C =11877.937142 Net atomic charges, atomic populations, and dipole contributions calculated with CM2 Atom NO. Type Charge No. of electrons 1 C -0.151 4.151 2 C 0.053 3.947 3 O -0.349 6.349 4 C 0.080 3.920 5 H 0.088 0.912 6 C 0.505 3.495 7 O -0.522 6.522 8 N -0.724 5.724 9 C 0.134 3.866 10 C -0.118 4.118 11 H 0.079 0.921 12 C -0.135 4.135 13 C 0.178 3.822 14 N -0.606 5.606 15 C -0.230 4.230 16 H 0.076 0.924 17 C -0.008 4.008 18 N -0.922 5.922 19 Si 0.912 3.088 20 H -0.285 1.285 21 H -0.290 1.290 22 H -0.285 1.285 23 C 0.143 3.857 24 C -0.098 4.098 25 C -0.151 4.151 26 C 0.062 3.938 27 O -0.383 6.383 28 C 0.061 3.939 29 C -0.146 4.146 30 H 0.063 0.937 31 H 0.065 0.935 32 H 0.070 0.930 33 H 0.048 0.952 34 H 0.050 0.950 35 H 0.400 0.600 36 H 0.069 0.931 37 H 0.064 0.936 38 H 0.061 0.939 39 H 0.059 0.941 40 H 0.046 0.954 41 H 0.043 0.957 42 H 0.253 0.747 43 H 0.019 0.981 44 H 0.034 0.966 45 H 0.086 0.914 46 H 0.079 0.921 47 H 0.072 0.928 48 H 0.051 0.949 49 H 0.093 0.907 50 H 0.053 0.947 51 H 0.095 0.905 52 H 0.074 0.926 53 H 0.083 0.917 For the charges calculated by CM2: Dipole moment (debyes) X Y Z Total from point charges -17.346 -13.665 9.161 23.907 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.208 4.208 2 C -0.024 4.024 3 O -0.266 6.266 4 C 0.018 3.982 5 H 0.106 0.894 6 C 0.289 3.711 7 O -0.398 6.398 8 N -0.378 5.378 9 C 0.010 3.990 10 C -0.138 4.138 11 H 0.097 0.903 12 C -0.174 4.174 13 C 0.053 3.947 14 N -0.332 5.332 15 C -0.358 4.358 16 H 0.094 0.906 17 C -0.204 4.204 18 N -0.569 5.569 19 Si 0.742 3.258 20 H -0.212 1.212 21 H -0.216 1.216 22 H -0.212 1.212 23 C 0.017 3.983 24 C -0.117 4.117 25 C -0.188 4.188 26 C -0.015 4.015 27 O -0.302 6.302 28 C -0.016 4.016 29 C -0.184 4.184 30 H 0.083 0.917 31 H 0.084 0.916 32 H 0.089 0.911 33 H 0.066 0.934 34 H 0.069 0.931 35 H 0.234 0.766 36 H 0.087 0.913 37 H 0.082 0.918 38 H 0.080 0.920 39 H 0.078 0.922 40 H 0.064 0.936 41 H 0.061 0.939 42 H 0.063 0.937 43 H 0.037 0.963 44 H 0.052 0.948 45 H 0.104 0.896 46 H 0.098 0.902 47 H 0.091 0.909 48 H 0.070 0.930 49 H 0.111 0.889 50 H 0.071 0.929 51 H 0.113 0.887 52 H 0.093 0.907 53 H 0.102 0.898 Dipole moment (debyes) X Y Z Total from point charges -18.021 -12.959 8.816 23.883 hybrid contribution 1.623 -0.232 0.322 1.671 sum -16.398 -13.191 9.137 22.943 Atomic orbital electron populations 1.21722 0.93461 1.02138 1.03448 1.22674 0.96052 0.83703 0.99960 1.88042 1.18491 1.26134 1.93981 1.21971 0.84048 0.98584 0.93582 0.89422 1.20619 0.82895 0.83088 0.84457 1.90716 1.47565 1.48297 1.53173 1.46066 1.06977 1.52659 1.32062 1.21402 0.95350 0.94950 0.87267 1.22618 0.92206 1.00152 0.98776 0.90272 1.22859 1.00228 0.98465 0.95831 1.21145 0.90058 0.99492 0.84023 1.44470 1.10406 1.76309 1.02056 1.18991 0.90127 1.36660 0.90054 0.90569 1.24808 0.98215 1.01582 0.95827 1.69953 1.28443 1.52223 1.06320 0.86203 0.82403 0.78403 0.78760 1.21190 1.21577 1.21151 1.21554 0.88553 0.94399 0.93788 1.21585 1.00248 0.96577 0.93271 1.22265 1.03131 0.98371 0.95055 1.22539 0.96915 0.85277 0.96793 1.87932 1.73190 1.20971 1.48089 1.22617 0.99751 0.91788 0.87454 1.22143 1.02553 0.95873 0.97789 0.91743 0.91613 0.91143 0.93382 0.93133 0.76557 0.91261 0.91760 0.92007 0.92185 0.93559 0.93889 0.93673 0.96323 0.94771 0.89583 0.90193 0.90905 0.93022 0.88910 0.92913 0.88712 0.90713 0.89800 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 -1.12 10.28 37.16 0.38 -0.74 16 2 C 0.05 0.43 6.05 37.16 0.22 0.66 16 3 O -0.35 -3.97 8.17 -55.12 -0.45 -4.42 16 4 C 0.08 0.81 2.39 -27.88 -0.07 0.74 16 5 H 0.09 0.68 7.59 -51.93 -0.39 0.28 16 6 C 0.50 6.68 6.07 -10.99 -0.07 6.61 16 7 O -0.52 -9.15 16.51 -14.33 -0.24 -9.39 16 8 N -0.72 -8.74 5.39 -60.29 -0.32 -9.06 16 9 C 0.13 1.87 5.13 -4.04 -0.02 1.85 16 10 C -0.12 -1.88 3.07 -88.47 -0.27 -2.15 16 11 H 0.08 1.23 8.14 -51.93 -0.42 0.81 16 12 C -0.13 -2.41 6.42 -24.90 -0.16 -2.57 16 13 C 0.18 4.22 5.37 -3.06 -0.02 4.20 16 14 N -0.61 -15.09 3.39 -170.08 -0.58 -15.67 16 15 C -0.23 -6.31 10.28 -15.06 -0.15 -6.46 16 16 H 0.08 2.00 7.83 -52.49 -0.41 1.59 16 17 C -0.01 -0.22 5.49 -17.43 -0.10 -0.32 16 18 N -0.92 -26.87 10.31 55.49 0.57 -26.29 16 19 Si 0.91 21.85 29.55 -169.99 -5.02 16.83 16 20 H -0.29 -6.70 7.11 56.52 0.40 -6.30 16 21 H -0.29 -7.02 7.11 56.52 0.40 -6.61 16 22 H -0.28 -6.75 7.11 56.52 0.40 -6.35 16 23 C 0.14 2.78 6.56 -2.53 -0.02 2.76 16 24 C -0.10 -0.88 2.08 -90.43 -0.19 -1.07 16 25 C -0.15 -1.14 4.89 -26.58 -0.13 -1.27 16 26 C 0.06 0.50 6.84 37.21 0.25 0.76 16 27 O -0.38 -4.19 10.77 -35.23 -0.38 -4.57 16 28 C 0.06 0.61 6.84 37.21 0.25 0.87 16 29 C -0.15 -1.57 5.02 -26.59 -0.13 -1.70 16 30 H 0.06 0.49 8.14 -51.93 -0.42 0.07 16 31 H 0.06 0.56 8.14 -51.93 -0.42 0.14 16 32 H 0.07 0.38 8.14 -51.93 -0.42 -0.05 16 33 H 0.05 0.37 8.10 -51.93 -0.42 -0.05 16 34 H 0.05 0.33 7.06 -51.93 -0.37 -0.04 16 35 H 0.40 4.03 7.91 -40.82 -0.32 3.70 16 36 H 0.07 1.07 8.04 -51.93 -0.42 0.65 16 37 H 0.06 0.89 8.14 -51.93 -0.42 0.47 16 38 H 0.06 0.99 8.14 -51.93 -0.42 0.57 16 39 H 0.06 1.02 8.04 -51.93 -0.42 0.61 16 40 H 0.05 1.21 7.38 -51.93 -0.38 0.83 16 41 H 0.04 1.14 7.50 -51.93 -0.39 0.75 16 42 H 0.25 7.18 8.31 -40.82 -0.34 6.84 16 43 H 0.02 0.38 7.96 -51.93 -0.41 -0.03 16 44 H 0.03 0.62 8.14 -51.93 -0.42 0.20 16 45 H 0.09 0.71 7.62 -51.93 -0.40 0.31 16 46 H 0.08 0.67 6.74 -51.93 -0.35 0.32 16 47 H 0.07 0.41 8.14 -51.93 -0.42 -0.01 16 48 H 0.05 0.39 8.14 -51.93 -0.42 -0.03 16 49 H 0.09 0.62 8.14 -51.93 -0.42 0.20 16 50 H 0.05 0.51 8.14 -51.93 -0.42 0.09 16 51 H 0.09 0.89 8.14 -51.93 -0.42 0.47 16 52 H 0.07 0.91 7.39 -51.93 -0.38 0.53 16 53 H 0.08 1.02 7.55 -51.93 -0.39 0.63 16 LS Contribution 410.93 15.07 6.19 6.19 Total: -1.00 -33.53 410.93 -10.10 -43.63 By element: Atomic # 1 Polarization: 10.25 SS G_CDS: -9.67 Total: 0.58 kcal Atomic # 6 Polarization: 2.38 SS G_CDS: -0.20 Total: 2.18 kcal Atomic # 7 Polarization: -50.70 SS G_CDS: -0.33 Total: -51.03 kcal Atomic # 8 Polarization: -17.31 SS G_CDS: -1.07 Total: -18.38 kcal Atomic # 14 Polarization: 21.85 SS G_CDS: -5.02 Total: 16.83 kcal Total LS contribution 6.19 Total: 6.19 kcal Total: -33.53 -10.10 -43.63 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. ZINC001714632118.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 -74.273 kcal (2) G-P(sol) polarization free energy of solvation -33.534 kcal (3) G-ENP(sol) elect.-nuc.-pol. free energy of system -107.806 kcal (4) G-CDS(sol) cavity-dispersion-solvent structure free energy -10.099 kcal (5) G-P-CDS(sol) = G-P(sol) + G-CDS(sol) = (2) + (4) -43.632 kcal (6) G-S(sol) free energy of system = (1) + (5) -117.905 kcal 1SCF WAS SPECIFIED, SO GEOMETRY OPTIMIZATION WAS NOT USED Total computer time = 0.87 seconds