Wall clock time and date at job start Sat Apr 11 2020 06:13:12 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.52996 * 1 3 3 H 1.08998 * 109.47526 * 2 1 4 4 C 1.53006 * 109.47034 * 119.99760 * 2 1 3 5 5 N 1.46498 * 109.47138 * 175.94892 * 4 2 1 6 6 C 1.46495 * 120.00558 * 264.48580 * 5 4 2 7 7 C 1.34782 * 119.99685 * 84.48601 * 5 4 2 8 8 O 1.21281 * 119.99755 * 359.97438 * 7 5 4 9 9 C 1.50699 * 119.99994 * 180.27811 * 7 5 4 10 10 H 1.08999 * 109.46221 * 39.72270 * 9 7 5 11 11 C 1.53043 * 109.45950 * 279.74342 * 9 7 5 12 12 C 1.53111 * 109.27749 * 180.79707 * 11 9 7 13 13 C 1.53123 * 109.15538 * 57.66887 * 12 11 9 14 14 C 1.53032 * 109.27584 * 302.34192 * 13 12 11 15 15 N 1.46845 * 109.45866 * 159.70315 * 9 7 5 16 16 C 1.46896 * 111.00224 * 305.84951 * 15 9 7 17 17 N 1.46903 * 109.46995 * 239.99617 * 2 1 3 18 18 C 1.46898 * 110.99855 * 80.60148 * 17 2 1 19 19 C 1.50703 * 109.47068 * 180.02562 * 18 17 2 20 20 O 1.21294 * 120.00003 * 0.02562 * 19 18 17 21 21 N 1.34778 * 119.99865 * 180.02562 * 19 18 17 22 22 C 1.37100 * 119.99928 * 180.02562 * 21 19 18 23 23 H 1.07999 * 119.99998 * 359.97438 * 22 21 19 24 24 C 1.38943 * 119.99869 * 179.97438 * 22 21 19 25 25 N 1.31405 * 120.00336 * 359.97438 * 24 22 21 26 26 Si 1.86803 * 119.99940 * 179.97438 * 24 22 21 27 27 H 1.48502 * 109.47319 * 59.99787 * 26 24 22 28 28 H 1.48504 * 109.47263 * 180.02562 * 26 24 22 29 29 H 1.48505 * 109.47182 * 300.00077 * 26 24 22 30 30 H 1.09000 * 109.47265 * 300.00369 * 1 2 3 31 31 H 1.08998 * 109.47286 * 59.99953 * 1 2 3 32 32 H 1.09000 * 109.47439 * 180.02562 * 1 2 3 33 33 H 1.09005 * 109.46816 * 55.95937 * 4 2 1 34 34 H 1.08997 * 109.47135 * 295.95569 * 4 2 1 35 35 H 1.08999 * 109.47092 * 264.59170 * 6 5 4 36 36 H 1.09007 * 109.47126 * 24.59201 * 6 5 4 37 37 H 1.09000 * 109.47527 * 144.58649 * 6 5 4 38 38 H 1.08998 * 109.50131 * 300.74970 * 11 9 7 39 39 H 1.08996 * 109.50099 * 60.86298 * 11 9 7 40 40 H 1.09004 * 109.53141 * 297.76305 * 12 11 9 41 41 H 1.09001 * 109.52541 * 177.57037 * 12 11 9 42 42 H 1.08996 * 109.50303 * 182.39745 * 13 12 11 43 43 H 1.09002 * 109.49882 * 62.28204 * 13 12 11 44 44 H 1.08999 * 109.46762 * 299.16066 * 14 13 12 45 45 H 1.09009 * 109.45253 * 179.18890 * 14 13 12 46 46 H 1.08995 * 109.46968 * 60.00366 * 16 15 9 47 47 H 1.08997 * 109.47298 * 180.02562 * 16 15 9 48 48 H 1.08997 * 109.47169 * 300.00479 * 16 15 9 49 49 H 1.00896 * 111.00102 * 204.55905 * 17 2 1 50 50 H 1.09007 * 109.47243 * 299.99771 * 18 17 2 51 51 H 1.08999 * 109.46974 * 59.99835 * 18 17 2 52 52 H 0.97001 * 119.99686 * 359.97438 * 21 19 18 53 53 H 0.96999 * 120.00394 * 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 6 1.5300 0.0000 0.0000 3 1 1.8934 1.0276 0.0000 4 6 2.0400 -0.7212 1.2493 5 7 3.5017 -0.8063 1.2016 6 6 4.1486 -2.0130 0.6807 7 6 4.2513 0.2257 1.6369 8 8 3.7157 1.2249 2.0677 9 6 5.7549 0.1408 1.5820 10 1 6.0576 -0.3446 0.6541 11 6 6.2634 -0.6744 2.7733 12 6 7.7923 -0.7413 2.7244 13 6 8.3579 0.6817 2.7316 14 6 7.7887 1.4576 1.5417 15 7 6.3236 1.4936 1.6347 16 6 5.7503 2.3470 0.5855 17 7 2.0196 -0.6926 -1.1994 18 6 1.9839 0.1936 -2.3704 19 6 2.4916 -0.5501 -3.5789 20 8 2.8477 -1.7049 -3.4749 21 7 2.5500 0.0680 -4.7751 22 6 3.0114 -0.6087 -5.8746 23 1 3.3285 -1.6369 -5.7820 24 6 3.0712 0.0284 -7.1079 25 7 2.6849 1.2794 -7.2207 26 14 3.7005 -0.8934 -8.6058 27 1 2.8407 -2.0782 -8.8553 28 1 3.6650 0.0027 -9.7895 29 1 5.0967 -1.3373 -8.3633 30 1 -0.3634 0.5139 0.8899 31 1 -0.3634 0.5138 -0.8900 32 1 -0.3634 -1.0276 0.0005 33 1 1.6186 -1.7259 1.2863 34 1 1.7370 -0.1672 2.1378 35 1 4.4033 -2.6749 1.5084 36 1 3.4676 -2.5261 0.0015 37 1 5.0563 -1.7361 0.1443 38 1 5.9493 -0.1970 3.7015 39 1 5.8534 -1.6832 2.7255 40 1 8.1057 -1.2532 1.8146 41 1 8.1612 -1.2843 3.5945 42 1 9.4444 0.6417 2.6540 43 1 8.0768 1.1802 3.6593 44 1 8.0817 0.9662 0.6140 45 1 8.1784 2.4756 1.5523 46 1 4.6712 2.4159 0.7226 47 1 6.1890 3.3428 0.6470 48 1 5.9655 1.9152 -0.3919 49 1 2.9479 -1.0593 -1.0516 50 1 0.9590 0.5190 -2.5487 51 1 2.6146 1.0635 -2.1871 52 1 2.2649 0.9915 -4.8583 53 1 2.7270 1.7242 -8.0816 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 ZINC001277906964.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:13:12 Heat of formation + Delta-G solvation = -60.285720 kcal Electronic energy + Delta-G solvation = -31601.390425 eV Core-core repulsion = 27508.106733 eV Total energy + Delta-G solvation = -4093.283693 eV No. of doubly occupied orbitals = 66 Molecular weight (most abundant/longest-lived isotopes) = 340.226 amu Computer time = 0.96 seconds Orbital eigenvalues (eV) -40.78677 -38.73322 -38.12355 -36.75301 -34.94629 -33.90676 -33.64382 -31.77724 -31.24311 -30.71336 -27.53787 -27.19084 -26.60522 -25.84529 -24.14659 -23.54026 -23.01034 -20.84576 -20.08310 -19.30729 -19.10773 -17.91121 -17.34222 -16.60824 -16.45655 -15.73621 -15.54146 -15.30198 -15.09070 -14.87323 -14.57289 -14.45372 -14.17382 -13.97451 -13.80100 -13.46314 -13.29370 -13.18749 -12.93055 -12.73991 -12.68985 -12.46875 -12.18965 -11.89806 -11.88366 -11.83434 -11.67461 -11.61677 -11.52803 -11.38941 -11.15319 -11.08082 -10.99854 -10.80508 -10.67842 -10.46227 -10.23230 -9.63969 -9.51837 -8.78830 -8.48165 -8.44841 -8.07356 -7.97520 -6.43091 -3.84278 2.47436 2.75488 3.18107 3.24094 3.28139 3.69864 3.89208 4.08516 4.22700 4.34432 4.47026 4.62047 4.69192 4.83723 4.89656 4.93570 5.02749 5.07745 5.18935 5.22024 5.24757 5.29460 5.32133 5.35184 5.39800 5.41735 5.53089 5.57321 5.59056 5.65002 5.67661 5.70689 5.78702 5.80484 5.82924 5.86270 5.99543 6.01327 6.06090 6.07572 6.22131 6.28158 6.45497 6.73699 6.86109 7.18381 7.28225 7.53032 7.68277 8.04294 8.07983 8.59483 9.15767 9.52666 10.03863 10.74550 Molecular weight = 340.23amu Principal moments of inertia in cm(-1) A = 0.010947 B = 0.002967 C = 0.002597 Principal moments of inertia in units of 10**(-40)*gram-cm**2 A = 2557.231058 B = 9434.765185 C =10779.920076 Net atomic charges, atomic populations, and dipole contributions calculated with CM2 Atom NO. Type Charge No. of electrons 1 C -0.137 4.137 2 C 0.065 3.935 3 H 0.045 0.955 4 C 0.118 3.882 5 N -0.617 5.617 6 C 0.066 3.934 7 C 0.518 3.482 8 O -0.519 6.519 9 C 0.067 3.933 10 H 0.065 0.935 11 C -0.106 4.106 12 C -0.120 4.120 13 C -0.118 4.118 14 C 0.054 3.946 15 N -0.520 5.520 16 C 0.013 3.987 17 N -0.663 5.663 18 C 0.050 3.950 19 C 0.441 3.559 20 O -0.585 6.585 21 N -0.563 5.563 22 C -0.301 4.301 23 H 0.102 0.898 24 C 0.023 3.977 25 N -0.901 5.901 26 Si 0.931 3.069 27 H -0.273 1.273 28 H -0.284 1.284 29 H -0.274 1.274 30 H 0.055 0.945 31 H 0.058 0.942 32 H 0.063 0.937 33 H 0.088 0.912 34 H 0.086 0.914 35 H 0.060 0.940 36 H 0.092 0.908 37 H 0.071 0.929 38 H 0.077 0.923 39 H 0.073 0.927 40 H 0.063 0.937 41 H 0.069 0.931 42 H 0.068 0.932 43 H 0.072 0.928 44 H 0.035 0.965 45 H 0.079 0.921 46 H 0.087 0.913 47 H 0.070 0.930 48 H 0.023 0.977 49 H 0.346 0.654 50 H 0.086 0.914 51 H 0.039 0.961 52 H 0.392 0.608 53 H 0.272 0.728 For the charges calculated by CM2: Dipole moment (debyes) X Y Z Total from point charges 6.409 -1.276 20.430 21.450 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.196 4.196 2 C -0.046 4.046 3 H 0.063 0.937 4 C -0.004 4.004 5 N -0.351 5.351 6 C -0.077 4.077 7 C 0.305 3.695 8 O -0.394 6.394 9 C -0.044 4.044 10 H 0.082 0.918 11 C -0.144 4.144 12 C -0.157 4.157 13 C -0.156 4.156 14 C -0.073 4.073 15 N -0.245 5.245 16 C -0.135 4.135 17 N -0.298 5.298 18 C -0.083 4.083 19 C 0.227 3.773 20 O -0.470 6.470 21 N -0.199 5.199 22 C -0.430 4.430 23 H 0.120 0.880 24 C -0.179 4.179 25 N -0.542 5.542 26 Si 0.758 3.242 27 H -0.198 1.198 28 H -0.210 1.210 29 H -0.199 1.199 30 H 0.074 0.926 31 H 0.077 0.923 32 H 0.082 0.918 33 H 0.106 0.894 34 H 0.104 0.896 35 H 0.079 0.921 36 H 0.110 0.890 37 H 0.089 0.911 38 H 0.095 0.905 39 H 0.091 0.909 40 H 0.081 0.919 41 H 0.088 0.912 42 H 0.087 0.913 43 H 0.091 0.909 44 H 0.053 0.947 45 H 0.097 0.903 46 H 0.106 0.894 47 H 0.088 0.912 48 H 0.041 0.959 49 H 0.168 0.832 50 H 0.104 0.896 51 H 0.057 0.943 52 H 0.225 0.775 53 H 0.082 0.918 Dipole moment (debyes) X Y Z Total from point charges 5.640 -1.474 20.864 21.663 hybrid contribution 1.067 0.439 -1.709 2.062 sum 6.707 -1.036 19.155 20.321 Atomic orbital electron populations 1.21729 0.93892 1.02503 1.01457 1.21922 0.96291 0.96460 0.89957 0.93713 1.21707 0.78513 1.02092 0.98121 1.47947 1.06858 1.17081 1.63230 1.22083 0.99393 0.86835 0.99425 1.20457 0.89071 0.82687 0.77238 1.90665 1.73815 1.30007 1.44918 1.21668 0.91558 0.91404 0.99791 0.91784 1.21532 0.96171 0.99098 0.97584 1.21581 0.94708 0.97173 1.02280 1.21612 0.99539 0.95404 0.99057 1.22156 0.86926 1.00896 0.97364 1.61675 1.05111 1.09159 1.48563 1.22594 1.02639 0.96972 0.91262 1.60287 1.31845 1.36022 1.01681 1.21709 1.01519 0.93997 0.91055 1.19803 0.85095 0.88079 0.84323 1.90496 1.53697 1.17599 1.85203 1.41863 1.56842 1.14377 1.06836 1.19646 1.39385 0.98211 0.85783 0.87988 1.25262 0.97558 0.95872 0.99196 1.69602 1.47567 1.10728 1.26347 0.85924 0.79383 0.78927 0.80011 1.19848 1.20952 1.19941 0.92618 0.92285 0.91761 0.89355 0.89581 0.92147 0.88964 0.91077 0.90461 0.90874 0.91855 0.91234 0.91317 0.90917 0.94704 0.90273 0.89433 0.91157 0.95856 0.83236 0.89554 0.94298 0.77469 0.91787 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.14 -1.64 9.36 37.16 0.35 -1.29 16 2 C 0.06 0.88 2.90 -67.71 -0.20 0.68 16 3 H 0.04 0.68 7.87 -51.93 -0.41 0.27 16 4 C 0.12 1.37 5.02 -4.04 -0.02 1.35 16 5 N -0.62 -7.14 2.46 -181.57 -0.45 -7.59 16 6 C 0.07 0.64 8.30 59.84 0.50 1.13 16 7 C 0.52 6.81 6.07 -10.98 -0.07 6.74 16 8 O -0.52 -8.39 14.19 5.56 0.08 -8.31 16 9 C 0.07 0.76 1.54 -68.87 -0.11 0.65 16 10 H 0.06 0.72 6.35 -51.93 -0.33 0.39 16 11 C -0.11 -0.94 5.28 -26.65 -0.14 -1.08 16 12 C -0.12 -0.89 6.00 -26.61 -0.16 -1.05 16 13 C -0.12 -0.92 6.09 -26.65 -0.16 -1.08 16 14 C 0.05 0.52 5.46 -3.83 -0.02 0.50 16 15 N -0.52 -6.21 4.51 -236.33 -1.07 -7.28 16 16 C 0.01 0.17 9.10 60.06 0.55 0.71 16 17 N -0.66 -10.89 5.92 -107.99 -0.64 -11.53 16 18 C 0.05 0.92 5.97 -4.98 -0.03 0.89 16 19 C 0.44 10.50 7.82 -10.98 -0.09 10.41 16 20 O -0.59 -15.65 15.78 5.55 0.09 -15.56 16 21 N -0.56 -14.39 5.29 -10.96 -0.06 -14.45 16 22 C -0.30 -8.39 9.68 -14.93 -0.14 -8.54 16 23 H 0.10 2.97 7.16 -52.49 -0.38 2.59 16 24 C 0.02 0.61 5.50 -17.47 -0.10 0.52 16 25 N -0.90 -24.63 13.05 55.49 0.72 -23.90 16 26 Si 0.93 21.01 29.55 -169.99 -5.02 15.99 16 27 H -0.27 -6.13 7.11 56.52 0.40 -5.72 16 28 H -0.28 -6.27 7.11 56.52 0.40 -5.86 16 29 H -0.27 -6.12 7.11 56.52 0.40 -5.72 16 30 H 0.05 0.57 8.14 -51.93 -0.42 0.14 16 31 H 0.06 0.75 7.19 -51.93 -0.37 0.37 16 32 H 0.06 0.75 8.14 -51.93 -0.42 0.33 16 33 H 0.09 0.89 8.03 -51.93 -0.42 0.48 16 34 H 0.09 1.00 7.59 -51.93 -0.39 0.60 16 35 H 0.06 0.44 7.38 -51.93 -0.38 0.06 16 36 H 0.09 1.01 6.97 -51.92 -0.36 0.64 16 37 H 0.07 0.69 6.17 -51.93 -0.32 0.37 16 38 H 0.08 0.78 8.14 -51.93 -0.42 0.35 16 39 H 0.07 0.55 6.78 -51.93 -0.35 0.20 16 40 H 0.06 0.46 8.14 -51.93 -0.42 0.04 16 41 H 0.07 0.43 8.14 -51.93 -0.42 0.01 16 42 H 0.07 0.44 8.14 -51.93 -0.42 0.01 16 43 H 0.07 0.60 8.14 -51.93 -0.42 0.18 16 44 H 0.03 0.34 8.14 -51.93 -0.42 -0.08 16 45 H 0.08 0.70 8.09 -51.92 -0.42 0.28 16 46 H 0.09 1.31 5.62 -51.93 -0.29 1.02 16 47 H 0.07 0.80 8.09 -51.93 -0.42 0.38 16 48 H 0.02 0.30 8.14 -51.93 -0.42 -0.12 16 49 H 0.35 5.77 4.93 -39.29 -0.19 5.58 16 50 H 0.09 1.50 7.16 -51.93 -0.37 1.12 16 51 H 0.04 0.69 8.06 -51.93 -0.42 0.27 16 52 H 0.39 9.68 7.90 -40.82 -0.32 9.35 16 53 H 0.27 7.05 8.31 -40.82 -0.34 6.71 16 LS Contribution 409.10 15.07 6.17 6.17 Total: -1.00 -32.56 409.10 -9.11 -41.67 By element: Atomic # 1 Polarization: 23.35 SS G_CDS: -9.09 Total: 14.25 kcal Atomic # 6 Polarization: 10.38 SS G_CDS: 0.16 Total: 10.54 kcal Atomic # 7 Polarization: -63.25 SS G_CDS: -1.48 Total: -64.74 kcal Atomic # 8 Polarization: -24.04 SS G_CDS: 0.17 Total: -23.88 kcal Atomic # 14 Polarization: 21.01 SS G_CDS: -5.02 Total: 15.99 kcal Total LS contribution 6.17 Total: 6.17 kcal Total: -32.56 -9.11 -41.67 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. ZINC001277906964.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 -18.615 kcal (2) G-P(sol) polarization free energy of solvation -32.561 kcal (3) G-ENP(sol) elect.-nuc.-pol. free energy of system -51.176 kcal (4) G-CDS(sol) cavity-dispersion-solvent structure free energy -9.109 kcal (5) G-P-CDS(sol) = G-P(sol) + G-CDS(sol) = (2) + (4) -41.671 kcal (6) G-S(sol) free energy of system = (1) + (5) -60.286 kcal 1SCF WAS SPECIFIED, SO GEOMETRY OPTIMIZATION WAS NOT USED Total computer time = 0.97 seconds