Wall clock time and date at job start Tue Mar 31 2020 21:05: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.52995 * 1 3 3 C 1.52999 * 109.47175 * 2 1 4 4 C 1.53004 * 109.47125 * 180.02562 * 3 2 1 5 5 H 1.09001 * 110.05405 * 306.61648 * 4 3 2 6 6 C 1.54258 * 110.03542 * 68.03336 * 4 3 2 7 7 C 1.54889 * 104.19753 * 142.89651 * 6 4 3 8 8 C 1.54888 * 102.74849 * 322.05835 * 7 6 4 9 9 C 1.53878 * 110.03300 * 185.20202 * 4 3 2 10 10 H 1.08998 * 110.07419 * 121.40338 * 9 4 3 11 11 C 1.50707 * 110.03194 * 359.97438 * 9 4 3 12 12 O 1.21280 * 119.99635 * 272.44250 * 11 9 4 13 13 N 1.34775 * 119.99855 * 92.43976 * 11 9 4 14 14 C 1.47123 * 121.06481 * 0.02562 * 13 11 9 15 15 C 1.53222 * 108.73364 * 232.28813 * 14 13 11 16 16 O 1.43131 * 109.34950 * 305.80363 * 15 14 13 17 17 C 1.42356 * 113.98816 * 59.77754 * 16 15 14 18 18 C 1.52921 * 114.29614 * 74.95129 * 17 16 15 19 19 N 1.47174 * 87.84764 * 88.27309 * 18 17 16 20 20 C 1.36934 * 136.00082 * 206.64799 * 19 18 17 21 21 H 1.08003 * 119.99978 * 359.97438 * 20 19 18 22 22 C 1.38811 * 119.99690 * 180.27528 * 20 19 18 23 23 N 1.31613 * 120.00370 * 359.72163 * 22 20 19 24 24 Si 1.86802 * 119.99650 * 179.72569 * 22 20 19 25 25 H 1.48498 * 109.47559 * 89.99509 * 24 22 20 26 26 H 1.48501 * 109.47170 * 209.99628 * 24 22 20 27 27 H 1.48500 * 109.47324 * 329.99541 * 24 22 20 28 28 C 1.47175 * 87.99673 * 26.90062 * 19 18 17 29 29 C 1.47363 * 121.05929 * 180.02562 * 13 11 9 30 30 H 1.08999 * 109.47271 * 180.02562 * 1 2 3 31 31 H 1.09001 * 109.47065 * 300.00294 * 1 2 3 32 32 H 1.09000 * 109.47214 * 59.99715 * 1 2 3 33 33 H 1.08996 * 109.47351 * 240.00052 * 2 1 3 34 34 H 1.09000 * 109.46921 * 119.99561 * 2 1 3 35 35 H 1.09001 * 109.47364 * 59.99889 * 3 2 1 36 36 H 1.09001 * 109.47356 * 299.99693 * 3 2 1 37 37 H 1.08995 * 110.48572 * 24.21665 * 6 4 3 38 38 H 1.09002 * 110.48537 * 261.57347 * 6 4 3 39 39 H 1.08998 * 110.75775 * 203.74255 * 7 6 4 40 40 H 1.09002 * 110.91471 * 80.44364 * 7 6 4 41 41 H 1.09001 * 110.48616 * 279.26623 * 8 7 6 42 42 H 1.09001 * 110.60198 * 156.67851 * 8 7 6 43 43 H 1.08999 * 109.59681 * 112.51267 * 14 13 11 44 44 H 1.09001 * 109.71372 * 352.13476 * 14 13 11 45 45 H 1.09008 * 109.49170 * 185.84577 * 15 14 13 46 46 H 1.08995 * 109.49725 * 65.76775 * 15 14 13 47 47 H 1.08998 * 113.49346 * 333.52856 * 18 17 16 48 48 H 1.09000 * 113.55028 * 203.02695 * 18 17 16 49 49 H 0.96998 * 119.99925 * 180.02562 * 23 22 20 50 50 H 1.08998 * 113.48920 * 87.84623 * 28 19 18 51 51 H 1.09002 * 113.48696 * 218.35059 * 28 19 18 52 52 H 1.09004 * 109.74391 * 12.11978 * 29 13 11 53 53 H 1.09000 * 109.74748 * 251.39901 * 29 13 11 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 6 2.0400 1.4425 0.0000 4 6 3.5700 1.4425 -0.0006 5 1 3.9441 0.8316 0.8210 6 6 4.0979 0.9006 -1.3450 7 6 5.3531 1.7606 -1.6347 8 6 4.9271 3.1648 -1.1389 9 6 4.0972 2.8822 0.1304 10 1 4.7269 2.9709 1.0157 11 6 2.9439 3.8483 0.2184 12 8 2.4510 4.2980 -0.7944 13 7 2.4602 4.2124 1.4225 14 6 3.0406 3.6816 2.6658 15 6 1.9140 3.0981 3.5249 16 8 0.9009 4.0909 3.7166 17 6 0.3148 4.5606 2.5073 18 6 -0.6634 3.5796 1.8597 19 7 -1.7251 4.1813 2.6823 20 6 -2.9237 3.7682 3.1997 21 1 -3.2732 2.7628 3.0167 22 6 -3.6926 4.6437 3.9540 23 7 -3.2671 5.8691 4.1767 24 14 -5.3314 4.0826 4.6534 25 1 -5.1201 3.4983 6.0021 26 1 -6.2499 5.2450 4.7559 27 1 -5.9271 3.0585 3.7581 28 6 -0.9363 5.4237 2.6995 29 6 1.3330 5.1548 1.5357 30 1 -0.3634 -1.0276 0.0005 31 1 -0.3633 0.5139 0.8900 32 1 -0.3634 0.5139 -0.8900 33 1 1.8933 -0.5138 -0.8899 34 1 1.8933 -0.5138 0.8900 35 1 1.6766 1.9564 0.8900 36 1 1.6767 1.9563 -0.8900 37 1 3.3535 1.0357 -2.1295 38 1 4.3692 -0.1511 -1.2522 39 1 5.5715 1.7778 -2.7024 40 1 6.2126 1.3966 -1.0719 41 1 4.3165 3.6657 -1.8901 42 1 5.8022 3.7677 -0.8961 43 1 3.5364 4.4854 3.2101 44 1 3.7628 2.9000 2.4301 45 1 2.3141 2.7953 4.4927 46 1 1.4840 2.2322 3.0217 47 1 -0.4607 2.5374 2.1063 48 1 -0.7982 3.7411 0.7902 49 1 -3.8042 6.4807 4.7041 50 1 -1.1439 6.0818 1.8558 51 1 -0.9648 5.9431 3.6574 52 1 0.8710 5.2955 0.5584 53 1 1.6899 6.1119 1.9161 RHF calculation, no. of doubly occupied orbitals= 65 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 ZINC001053320342.mol2 53 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 21:05:40 Heat of formation + Delta-G solvation = -27.462570 kcal Electronic energy + Delta-G solvation = -32214.979413 eV Core-core repulsion = 28281.282340 eV Total energy + Delta-G solvation = -3933.697073 eV No. of doubly occupied orbitals = 65 Molecular weight (most abundant/longest-lived isotopes) = 336.220 amu Computer time = 0.97 seconds Orbital eigenvalues (eV) -40.97314 -39.42516 -37.47908 -36.32773 -34.57382 -34.37234 -32.27268 -30.48281 -29.82262 -29.35076 -27.75363 -26.06955 -25.11880 -24.88077 -23.31967 -22.33841 -21.76098 -20.87720 -20.48440 -18.94633 -18.38569 -17.24526 -16.80675 -16.19684 -16.00285 -15.61902 -15.39141 -14.94105 -14.70252 -14.60752 -14.41011 -14.19796 -13.81569 -13.55647 -13.44890 -13.37342 -13.01583 -12.69901 -12.38356 -12.32631 -12.16239 -11.98160 -11.71440 -11.41274 -11.36196 -11.23231 -11.10846 -10.97114 -10.88816 -10.82881 -10.80627 -10.78953 -10.72734 -10.50007 -10.38625 -10.13877 -9.97920 -9.63484 -9.48247 -9.08763 -8.81835 -8.26539 -6.86325 -5.73890 -3.05769 2.69009 2.96530 3.41764 3.45906 3.48103 4.05345 4.22573 4.36460 4.40251 4.46500 4.71643 4.87170 4.94678 4.97780 5.00380 5.10068 5.14644 5.20088 5.22127 5.30338 5.31632 5.34681 5.40625 5.44095 5.46034 5.59442 5.60032 5.69025 5.72837 5.76296 5.80820 5.82268 5.84159 5.88323 5.97138 6.02423 6.04583 6.07092 6.19266 6.23176 6.29748 6.33409 6.37518 6.45982 6.52620 6.59990 6.95297 7.18392 7.45675 7.51625 7.83605 8.00038 8.24868 9.05322 9.72600 10.63661 11.43957 Molecular weight = 336.22amu Principal moments of inertia in cm(-1) A = 0.017263 B = 0.003540 C = 0.003217 Principal moments of inertia in units of 10**(-40)*gram-cm**2 A = 1621.558250 B = 7907.941264 C = 8700.903067 Net atomic charges, atomic populations, and dipole contributions calculated with CM2 Atom NO. Type Charge No. of electrons 1 C -0.155 4.155 2 C -0.117 4.117 3 C -0.129 4.129 4 C -0.075 4.075 5 H 0.078 0.922 6 C -0.117 4.117 7 C -0.123 4.123 8 C -0.113 4.113 9 C -0.113 4.113 10 H 0.090 0.910 11 C 0.529 3.471 12 O -0.540 6.540 13 N -0.612 5.612 14 C 0.071 3.929 15 C 0.045 3.955 16 O -0.331 6.331 17 C 0.046 3.954 18 C 0.104 3.896 19 N -0.599 5.599 20 C -0.215 4.215 21 H 0.080 0.920 22 C -0.006 4.006 23 N -0.919 5.919 24 Si 0.913 3.087 25 H -0.288 1.288 26 H -0.290 1.290 27 H -0.282 1.282 28 C 0.174 3.826 29 C 0.102 3.898 30 H 0.051 0.949 31 H 0.058 0.942 32 H 0.053 0.947 33 H 0.059 0.941 34 H 0.057 0.943 35 H 0.070 0.930 36 H 0.072 0.928 37 H 0.070 0.930 38 H 0.066 0.934 39 H 0.067 0.933 40 H 0.065 0.935 41 H 0.086 0.914 42 H 0.065 0.935 43 H 0.076 0.924 44 H 0.091 0.909 45 H 0.100 0.900 46 H 0.054 0.946 47 H 0.040 0.960 48 H 0.039 0.961 49 H 0.254 0.746 50 H 0.042 0.958 51 H 0.083 0.917 52 H 0.101 0.899 53 H 0.072 0.928 For the charges calculated by CM2: Dipole moment (debyes) X Y Z Total from point charges 19.824 -9.072 -7.600 23.088 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.212 4.212 2 C -0.155 4.155 3 C -0.167 4.167 4 C -0.094 4.094 5 H 0.096 0.904 6 C -0.154 4.154 7 C -0.160 4.160 8 C -0.151 4.151 9 C -0.135 4.135 10 H 0.108 0.892 11 C 0.318 3.682 12 O -0.417 6.417 13 N -0.345 5.345 14 C -0.052 4.052 15 C -0.033 4.033 16 O -0.249 6.249 17 C 0.004 3.996 18 C -0.022 4.022 19 N -0.324 5.324 20 C -0.342 4.342 21 H 0.098 0.902 22 C -0.203 4.203 23 N -0.566 5.566 24 Si 0.743 3.257 25 H -0.215 1.215 26 H -0.216 1.216 27 H -0.207 1.207 28 C 0.048 3.952 29 C -0.021 4.021 30 H 0.070 0.930 31 H 0.077 0.923 32 H 0.072 0.928 33 H 0.077 0.923 34 H 0.076 0.924 35 H 0.088 0.912 36 H 0.091 0.909 37 H 0.089 0.911 38 H 0.084 0.916 39 H 0.086 0.914 40 H 0.084 0.916 41 H 0.104 0.896 42 H 0.084 0.916 43 H 0.094 0.906 44 H 0.109 0.891 45 H 0.118 0.882 46 H 0.072 0.928 47 H 0.058 0.942 48 H 0.056 0.944 49 H 0.064 0.936 50 H 0.060 0.940 51 H 0.101 0.899 52 H 0.119 0.881 53 H 0.090 0.910 Dipole moment (debyes) X Y Z Total from point charges 19.147 -8.725 -7.761 22.427 hybrid contribution -0.808 -0.327 0.292 0.919 sum 18.339 -9.052 -7.469 21.772 Atomic orbital electron populations 1.21783 0.95573 1.01358 1.02490 1.21477 0.96138 0.96418 1.01447 1.22003 0.93735 0.97469 1.03469 1.21707 0.96048 0.93061 0.98598 0.90406 1.22440 0.97289 0.98916 0.96756 1.22472 0.97598 0.95018 1.00957 1.22367 0.98880 0.97514 0.96320 1.22108 0.95614 0.94961 1.00792 0.89163 1.20377 0.84180 0.80983 0.82671 1.90668 1.60527 1.55983 1.34517 1.48262 1.34568 1.46131 1.05590 1.22470 0.97011 0.99771 0.85960 1.22892 0.88576 0.91855 0.99997 1.88029 1.52380 1.56532 1.27920 1.23384 0.91654 0.99098 0.85428 1.22855 0.88383 0.96343 0.94649 1.45376 1.20207 1.09767 1.57081 1.18869 0.93049 0.96774 1.25477 0.90197 1.24765 0.98909 0.95872 1.00747 1.69804 1.38301 1.05161 1.43309 0.86111 0.83226 0.77322 0.79032 1.21508 1.21625 1.20725 1.22062 0.81755 0.90116 1.01317 1.22161 0.86585 0.91840 1.01495 0.93010 0.92295 0.92756 0.92252 0.92376 0.91182 0.90935 0.91111 0.91568 0.91416 0.91587 0.89560 0.91613 0.90572 0.89133 0.88163 0.92809 0.94225 0.94361 0.93618 0.94012 0.89856 0.88104 0.90955 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.87 9.91 37.16 0.37 -1.50 16 2 C -0.12 -1.23 5.79 -26.73 -0.15 -1.39 16 3 C -0.13 -1.51 3.40 -26.73 -0.09 -1.61 16 4 C -0.08 -0.70 2.92 -89.50 -0.26 -0.97 16 5 H 0.08 0.62 8.14 -51.93 -0.42 0.20 16 6 C -0.12 -0.98 5.54 -25.08 -0.14 -1.12 16 7 C -0.12 -0.95 6.91 -24.74 -0.17 -1.12 16 8 C -0.11 -1.08 6.31 -25.07 -0.16 -1.24 16 9 C -0.11 -1.15 2.47 -90.65 -0.22 -1.37 16 10 H 0.09 0.74 6.23 -51.93 -0.32 0.42 16 11 C 0.53 7.28 4.59 -10.98 -0.05 7.23 16 12 O -0.54 -9.31 14.57 5.56 0.08 -9.22 16 13 N -0.61 -8.40 2.44 -172.49 -0.42 -8.82 16 14 C 0.07 0.79 6.17 -3.59 -0.02 0.77 16 15 C 0.04 0.61 6.90 37.27 0.26 0.86 16 16 O -0.33 -6.36 10.55 -35.23 -0.37 -6.73 16 17 C 0.05 0.90 1.29 -90.68 -0.12 0.78 16 18 C 0.10 2.23 7.90 -8.21 -0.06 2.17 16 19 N -0.60 -15.78 3.78 -184.58 -0.70 -16.48 16 20 C -0.21 -6.15 10.74 -15.06 -0.16 -6.31 16 21 H 0.08 2.22 7.83 -52.48 -0.41 1.81 16 22 C -0.01 -0.18 5.50 -17.43 -0.10 -0.28 16 23 N -0.92 -27.59 11.51 55.49 0.64 -26.95 16 24 Si 0.91 22.48 29.55 -169.99 -5.02 17.45 16 25 H -0.29 -7.05 7.11 56.52 0.40 -6.65 16 26 H -0.29 -7.17 7.11 56.52 0.40 -6.77 16 27 H -0.28 -6.74 7.11 56.52 0.40 -6.34 16 28 C 0.17 4.17 6.83 -8.06 -0.06 4.12 16 29 C 0.10 1.67 5.83 -3.69 -0.02 1.65 16 30 H 0.05 0.56 8.14 -51.93 -0.42 0.14 16 31 H 0.06 0.82 8.00 -51.93 -0.42 0.40 16 32 H 0.05 0.68 8.14 -51.93 -0.42 0.26 16 33 H 0.06 0.56 7.52 -51.93 -0.39 0.17 16 34 H 0.06 0.56 8.14 -51.93 -0.42 0.14 16 35 H 0.07 0.94 5.65 -51.93 -0.29 0.65 16 36 H 0.07 1.02 7.17 -51.93 -0.37 0.65 16 37 H 0.07 0.68 7.75 -51.93 -0.40 0.28 16 38 H 0.07 0.48 8.14 -51.93 -0.42 0.05 16 39 H 0.07 0.49 8.14 -51.93 -0.42 0.07 16 40 H 0.07 0.43 8.14 -51.93 -0.42 0.01 16 41 H 0.09 1.03 6.87 -51.93 -0.36 0.67 16 42 H 0.07 0.55 8.14 -51.93 -0.42 0.13 16 43 H 0.08 0.79 8.14 -51.93 -0.42 0.36 16 44 H 0.09 0.74 5.61 -51.93 -0.29 0.45 16 45 H 0.10 1.12 8.14 -51.92 -0.42 0.69 16 46 H 0.05 0.75 6.39 -51.93 -0.33 0.42 16 47 H 0.04 0.83 7.22 -51.93 -0.37 0.45 16 48 H 0.04 0.83 7.93 -51.93 -0.41 0.41 16 49 H 0.25 7.37 8.31 -40.82 -0.34 7.03 16 50 H 0.04 0.99 8.09 -51.93 -0.42 0.57 16 51 H 0.08 2.24 7.14 -51.93 -0.37 1.87 16 52 H 0.10 1.77 6.90 -51.93 -0.36 1.42 16 53 H 0.07 1.12 8.14 -51.93 -0.42 0.70 16 LS Contribution 396.88 15.07 5.98 5.98 Total: -1.00 -33.16 396.88 -10.28 -43.44 By element: Atomic # 1 Polarization: 9.97 SS G_CDS: -9.31 Total: 0.66 kcal Atomic # 6 Polarization: 1.84 SS G_CDS: -1.16 Total: 0.68 kcal Atomic # 7 Polarization: -51.77 SS G_CDS: -0.48 Total: -52.25 kcal Atomic # 8 Polarization: -15.66 SS G_CDS: -0.29 Total: -15.96 kcal Atomic # 14 Polarization: 22.48 SS G_CDS: -5.02 Total: 17.45 kcal Total LS contribution 5.98 Total: 5.98 kcal Total: -33.16 -10.28 -43.44 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. ZINC001053320342.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 15.978 kcal (2) G-P(sol) polarization free energy of solvation -33.156 kcal (3) G-ENP(sol) elect.-nuc.-pol. free energy of system -17.178 kcal (4) G-CDS(sol) cavity-dispersion-solvent structure free energy -10.284 kcal (5) G-P-CDS(sol) = G-P(sol) + G-CDS(sol) = (2) + (4) -43.440 kcal (6) G-S(sol) free energy of system = (1) + (5) -27.463 kcal 1SCF WAS SPECIFIED, SO GEOMETRY OPTIMIZATION WAS NOT USED Total computer time = 0.97 seconds