Wall clock time and date at job start Sat Apr 11 2020 06:09:25 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.31002 * 1 3 3 C 1.50698 * 120.00135 * 2 1 4 4 C 1.52998 * 109.47018 * 119.99967 * 3 2 1 5 5 C 1.53001 * 109.47472 * 239.99962 * 3 2 1 6 6 C 1.53001 * 109.46909 * 300.00308 * 5 3 2 7 7 O 1.42903 * 109.46710 * 179.97438 * 6 5 3 8 8 C 1.42899 * 113.99733 * 179.97438 * 7 6 5 9 9 C 1.50694 * 109.47122 * 0.02562 * 3 2 1 10 10 O 1.21284 * 120.00460 * 232.80790 * 9 3 2 11 11 N 1.34777 * 120.00048 * 52.81071 * 9 3 2 12 12 C 1.47502 * 118.75855 * 3.69737 * 11 9 3 13 13 C 1.53947 * 110.82361 * 138.25183 * 12 11 9 14 14 C 1.51541 * 109.72528 * 320.31349 * 13 12 11 15 15 O 1.42146 * 108.76770 * 91.29277 * 14 13 12 16 16 C 1.38014 * 119.43738 * 272.78148 * 15 14 13 17 17 H 1.09005 * 109.46663 * 303.42523 * 16 15 14 18 18 C 1.53001 * 109.47227 * 183.43453 * 16 15 14 19 19 N 1.46493 * 109.47386 * 64.98332 * 18 16 15 20 20 C 1.36938 * 120.00026 * 180.02562 * 19 18 16 21 21 H 1.07996 * 119.99815 * 359.97438 * 20 19 18 22 22 C 1.38811 * 119.99692 * 180.02562 * 20 19 18 23 23 N 1.31622 * 119.99490 * 180.02562 * 22 20 19 24 24 Si 1.86800 * 120.00457 * 359.97438 * 22 20 19 25 25 H 1.48499 * 109.47025 * 90.00107 * 24 22 20 26 26 H 1.48503 * 109.47111 * 210.00523 * 24 22 20 27 27 H 1.48506 * 109.46852 * 330.00451 * 24 22 20 28 28 C 1.45129 * 118.75837 * 183.96808 * 11 9 3 29 29 H 1.07995 * 119.99826 * 0.02562 * 1 2 3 30 30 H 1.07992 * 119.99937 * 179.71896 * 1 2 3 31 31 H 1.07997 * 119.99766 * 180.02562 * 2 1 3 32 32 H 1.09002 * 109.47068 * 37.26912 * 4 3 2 33 33 H 1.09004 * 109.47314 * 157.26710 * 4 3 2 34 34 H 1.09003 * 109.47460 * 277.26867 * 4 3 2 35 35 H 1.09006 * 109.47197 * 60.00116 * 5 3 2 36 36 H 1.08998 * 109.47448 * 179.97438 * 5 3 2 37 37 H 1.08996 * 109.46853 * 299.99366 * 6 5 3 38 38 H 1.09000 * 109.47593 * 60.00017 * 6 5 3 39 39 H 1.09008 * 109.47175 * 60.00476 * 8 7 6 40 40 H 1.08999 * 109.47237 * 179.97438 * 8 7 6 41 41 H 1.08997 * 109.47363 * 299.99817 * 8 7 6 42 42 H 1.08999 * 109.17123 * 258.77415 * 12 11 9 43 43 H 1.09004 * 109.16245 * 17.97493 * 12 11 9 44 44 H 1.08994 * 109.46341 * 200.17864 * 13 12 11 45 45 H 1.09002 * 109.27924 * 80.32099 * 13 12 11 46 46 H 1.08993 * 109.58642 * 331.50154 * 14 13 12 47 47 H 1.08999 * 109.70196 * 211.14325 * 14 13 12 48 48 H 1.08996 * 109.47002 * 184.98774 * 18 16 15 49 49 H 1.09007 * 109.46698 * 304.98475 * 18 16 15 50 50 H 0.97005 * 119.99780 * 0.02562 * 19 18 16 51 51 H 0.97002 * 119.99560 * 179.97438 * 23 22 20 52 52 H 1.09003 * 109.46038 * 29.30918 * 28 11 9 53 53 H 1.08999 * 109.46394 * 149.27323 * 28 11 9 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.3100 0.0000 0.0000 3 6 2.0635 1.3051 0.0000 4 6 2.9431 1.3861 1.2492 5 6 2.9432 1.3862 -1.2492 6 6 3.9374 0.2232 -1.2492 7 8 4.7586 0.2986 -2.4163 8 6 5.7275 -0.7472 -2.5138 9 6 1.0843 2.4505 0.0006 10 8 1.1725 3.3242 -0.8359 11 7 0.1105 2.5041 0.9309 12 6 0.0782 1.4725 1.9846 13 6 -0.2471 2.0981 3.3531 14 6 0.4712 3.4251 3.4926 15 8 -0.3833 4.4593 3.0225 16 6 -0.3801 4.7846 1.6813 17 1 0.6318 5.0488 1.3739 18 6 -1.3123 5.9742 1.4427 19 7 -0.7828 7.1549 2.1295 20 6 -1.4536 8.3466 2.0582 21 1 -2.3727 8.4161 1.4954 22 6 -0.9516 9.4655 2.7086 23 7 -1.5966 10.6108 2.6405 24 14 0.6385 9.3455 3.6815 25 1 1.7873 9.6539 2.7924 26 1 0.6065 10.3170 4.8041 27 1 0.7858 7.9699 4.2214 28 6 -0.8702 3.5711 0.8544 29 1 -0.5399 0.9353 -0.0004 30 1 -0.5400 -0.9352 -0.0046 31 1 1.8500 -0.9353 -0.0004 32 1 3.3512 0.3999 1.4708 33 1 3.7599 2.0862 1.0738 34 1 2.3450 1.7298 2.0932 35 1 2.3166 1.3285 -2.1393 36 1 3.4879 2.3303 -1.2494 37 1 4.5640 0.2810 -0.3593 38 1 3.3924 -0.7208 -1.2493 39 1 5.2199 -1.7115 -2.5395 40 1 6.3098 -0.6192 -3.4263 41 1 6.3915 -0.7090 -1.6503 42 1 -0.6812 0.7309 1.7365 43 1 1.0511 0.9835 2.0367 44 1 0.0770 1.4247 4.1466 45 1 -1.3226 2.2578 3.4293 46 1 1.3860 3.4094 2.9003 47 1 0.7168 3.6018 4.5397 48 1 -1.3787 6.1750 0.3735 49 1 -2.3040 5.7421 1.8312 50 1 0.0429 7.0925 2.6347 51 1 -1.2460 11.3926 3.0952 52 1 -1.0018 3.8696 -0.1857 53 1 -1.8213 3.2196 1.2542 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 ZINC001201080021.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:09:25 Heat of formation + Delta-G solvation = -49.349080 kcal Electronic energy + Delta-G solvation = -32090.019171 eV Core-core repulsion = 27964.916976 eV Total energy + Delta-G solvation = -4125.102195 eV No. of doubly occupied orbitals = 66 Molecular weight (most abundant/longest-lived isotopes) = 340.220 amu Computer time = 1.07 seconds Orbital eigenvalues (eV) -41.68795 -39.11289 -38.19921 -36.81854 -35.64691 -34.80235 -32.89171 -32.10327 -31.13473 -29.84286 -28.97574 -27.00128 -26.55158 -25.84505 -24.65272 -22.75462 -22.24645 -21.60986 -21.06030 -20.74938 -18.52229 -17.75873 -17.27918 -17.07963 -16.65304 -16.41474 -16.33244 -15.88066 -15.45433 -15.28543 -14.88998 -14.43855 -14.32824 -14.21841 -13.96322 -13.58558 -13.45577 -13.36890 -13.06550 -12.78998 -12.65946 -12.47801 -12.39685 -12.21251 -12.14138 -12.02002 -11.80325 -11.69246 -11.31094 -11.28709 -11.11460 -10.94919 -10.86771 -10.80323 -10.39769 -10.06266 -9.99254 -9.74035 -9.67477 -9.51386 -9.36349 -8.82599 -8.64871 -7.13610 -5.68553 -3.02826 1.79285 2.40895 2.68201 3.40322 3.45881 3.50528 3.52080 3.82834 3.98655 4.26393 4.49041 4.54170 4.64430 4.68457 4.70023 4.77000 4.79910 4.84741 4.89215 4.98240 5.02614 5.07517 5.09009 5.16162 5.20110 5.27713 5.31676 5.34026 5.36753 5.41284 5.49937 5.53627 5.61701 5.62726 5.72398 5.78106 5.82082 5.90262 5.94569 5.99566 6.02314 6.14478 6.25778 6.40802 6.59229 6.70382 6.79627 7.08792 7.86491 7.97488 8.07358 8.50499 9.30114 10.02753 10.23531 11.49940 Molecular weight = 340.22amu Principal moments of inertia in cm(-1) A = 0.015000 B = 0.002865 C = 0.002663 Principal moments of inertia in units of 10**(-40)*gram-cm**2 A = 1866.229916 B = 9771.155326 C =10510.208033 Net atomic charges, atomic populations, and dipole contributions calculated with CM2 Atom NO. Type Charge No. of electrons 1 C -0.193 4.193 2 C -0.166 4.166 3 C -0.045 4.045 4 C -0.156 4.156 5 C -0.091 4.091 6 C 0.062 3.938 7 O -0.390 6.390 8 C 0.032 3.968 9 C 0.519 3.481 10 O -0.521 6.521 11 N -0.607 5.607 12 C 0.082 3.918 13 C -0.170 4.170 14 C 0.071 3.929 15 O -0.371 6.371 16 C 0.085 3.915 17 H 0.052 0.948 18 C 0.171 3.829 19 N -0.750 5.750 20 C -0.240 4.240 21 H 0.090 0.910 22 C 0.003 3.997 23 N -0.936 5.936 24 Si 0.883 3.117 25 H -0.285 1.285 26 H -0.292 1.292 27 H -0.275 1.275 28 C 0.080 3.920 29 H 0.119 0.881 30 H 0.104 0.896 31 H 0.120 0.880 32 H 0.072 0.928 33 H 0.070 0.930 34 H 0.074 0.926 35 H 0.078 0.922 36 H 0.083 0.917 37 H 0.051 0.949 38 H 0.051 0.949 39 H 0.040 0.960 40 H 0.086 0.914 41 H 0.040 0.960 42 H 0.081 0.919 43 H 0.102 0.898 44 H 0.084 0.916 45 H 0.083 0.917 46 H 0.047 0.953 47 H 0.088 0.912 48 H 0.044 0.956 49 H 0.035 0.965 50 H 0.363 0.637 51 H 0.258 0.742 52 H 0.095 0.905 53 H 0.086 0.914 For the charges calculated by CM2: Dipole moment (debyes) X Y Z Total from point charges 10.271 -31.769 -2.968 33.519 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.229 4.229 2 C -0.184 4.184 3 C -0.048 4.048 4 C -0.213 4.213 5 C -0.129 4.129 6 C -0.014 4.014 7 O -0.309 6.309 8 C -0.065 4.065 9 C 0.306 3.694 10 O -0.395 6.395 11 N -0.341 5.341 12 C -0.040 4.040 13 C -0.208 4.208 14 C -0.006 4.006 15 O -0.289 6.289 16 C 0.025 3.975 17 H 0.070 0.930 18 C 0.043 3.957 19 N -0.395 5.395 20 C -0.370 4.370 21 H 0.108 0.892 22 C -0.192 4.192 23 N -0.584 5.584 24 Si 0.714 3.286 25 H -0.212 1.212 26 H -0.218 1.218 27 H -0.201 1.201 28 C -0.044 4.044 29 H 0.137 0.863 30 H 0.122 0.878 31 H 0.138 0.862 32 H 0.091 0.909 33 H 0.089 0.911 34 H 0.093 0.907 35 H 0.097 0.903 36 H 0.102 0.898 37 H 0.070 0.930 38 H 0.069 0.931 39 H 0.058 0.942 40 H 0.104 0.896 41 H 0.059 0.941 42 H 0.099 0.901 43 H 0.119 0.881 44 H 0.103 0.897 45 H 0.101 0.899 46 H 0.065 0.935 47 H 0.106 0.894 48 H 0.061 0.939 49 H 0.053 0.947 50 H 0.196 0.804 51 H 0.068 0.932 52 H 0.113 0.887 53 H 0.104 0.896 Dipole moment (debyes) X Y Z Total from point charges 9.987 -30.880 -3.491 32.642 hybrid contribution 1.320 -0.547 0.826 1.651 sum 11.307 -31.428 -2.665 33.506 Atomic orbital electron populations 1.24646 0.96564 1.03746 0.97957 1.22663 0.95808 0.96535 1.03378 1.19814 0.94995 0.92253 0.97732 1.22235 0.99811 1.03363 0.95854 1.21288 0.96266 0.98891 0.96423 1.22614 0.92502 0.98487 0.87822 1.87901 1.49291 1.50041 1.43702 1.22739 0.91512 0.90215 1.02012 1.20428 0.81288 0.86671 0.80991 1.90760 1.71515 1.39292 1.37947 1.48040 1.29070 1.28748 1.28212 1.23082 1.03990 0.90272 0.86620 1.22633 1.03366 0.98305 0.96527 1.22157 0.92406 0.86397 0.99590 1.86426 1.63331 1.54322 1.24810 1.22324 1.00008 0.93801 0.81410 0.93019 1.20589 0.94727 0.84178 0.96216 1.42570 1.28291 1.06877 1.61800 1.19059 1.05742 0.85265 1.26886 0.89185 1.24645 0.98551 0.95515 1.00525 1.69709 1.40214 1.00255 1.48249 0.86540 0.84758 0.77137 0.80209 1.21192 1.21806 1.20054 1.22435 0.91443 0.87526 1.02982 0.86305 0.87785 0.86178 0.90880 0.91092 0.90716 0.90287 0.89810 0.93036 0.93131 0.94164 0.89584 0.94091 0.90093 0.88061 0.89734 0.89872 0.93518 0.89391 0.93872 0.94671 0.80375 0.93158 0.88691 0.89555 Number of geometries 1 Number of calculations of the screened coulomb radii 1 The total number of SCF iterations 14. 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.19 -1.36 10.47 29.01 0.30 -1.06 16 2 C -0.17 -1.14 6.79 -36.03 -0.24 -1.39 16 3 C -0.04 -0.40 0.51 -156.82 -0.08 -0.48 16 4 C -0.16 -1.12 6.97 37.16 0.26 -0.86 16 5 C -0.09 -0.92 3.84 -26.73 -0.10 -1.02 16 6 C 0.06 0.50 5.41 37.16 0.20 0.70 16 7 O -0.39 -3.89 10.70 -35.23 -0.38 -4.27 16 8 C 0.03 0.20 11.01 101.05 1.11 1.31 16 9 C 0.52 6.54 5.49 -10.99 -0.06 6.48 16 10 O -0.52 -8.83 15.45 5.55 0.09 -8.74 16 11 N -0.61 -6.78 0.93 -172.24 -0.16 -6.94 16 12 C 0.08 0.65 3.52 -3.00 -0.01 0.64 16 13 C -0.17 -1.46 5.79 -26.98 -0.16 -1.61 16 14 C 0.07 0.85 6.37 36.42 0.23 1.09 16 15 O -0.37 -6.16 9.38 -26.53 -0.25 -6.41 16 16 C 0.08 1.44 2.68 -25.78 -0.07 1.37 16 17 H 0.05 0.98 8.14 -51.93 -0.42 0.56 16 18 C 0.17 3.52 5.92 -4.05 -0.02 3.50 16 19 N -0.75 -19.14 5.34 -59.91 -0.32 -19.46 16 20 C -0.24 -7.12 10.47 -15.06 -0.16 -7.28 16 21 H 0.09 2.82 8.06 -52.49 -0.42 2.40 16 22 C 0.00 0.09 5.50 -17.43 -0.10 0.00 16 23 N -0.94 -30.38 13.97 55.49 0.78 -29.60 16 24 Si 0.88 23.22 27.74 -169.99 -4.71 18.51 16 25 H -0.29 -7.45 7.11 56.52 0.40 -7.05 16 26 H -0.29 -7.67 7.11 56.52 0.40 -7.27 16 27 H -0.27 -6.90 6.52 56.52 0.37 -6.53 16 28 C 0.08 1.07 4.83 -3.83 -0.02 1.05 16 29 H 0.12 1.07 4.38 -52.49 -0.23 0.84 16 30 H 0.10 0.56 8.06 -52.49 -0.42 0.13 16 31 H 0.12 0.60 7.00 -52.49 -0.37 0.23 16 32 H 0.07 0.35 7.54 -51.93 -0.39 -0.04 16 33 H 0.07 0.56 7.93 -51.93 -0.41 0.15 16 34 H 0.07 0.55 4.33 -51.93 -0.22 0.32 16 35 H 0.08 0.92 8.14 -51.93 -0.42 0.50 16 36 H 0.08 1.00 8.00 -51.93 -0.42 0.58 16 37 H 0.05 0.35 6.99 -51.93 -0.36 -0.02 16 38 H 0.05 0.33 6.62 -51.93 -0.34 -0.01 16 39 H 0.04 0.21 8.14 -51.92 -0.42 -0.22 16 40 H 0.09 0.46 8.14 -51.93 -0.42 0.04 16 41 H 0.04 0.22 8.14 -51.93 -0.42 -0.21 16 42 H 0.08 0.54 5.64 -51.93 -0.29 0.25 16 43 H 0.10 0.67 3.19 -51.93 -0.17 0.50 16 44 H 0.08 0.52 8.14 -51.93 -0.42 0.09 16 45 H 0.08 0.80 8.01 -51.93 -0.42 0.38 16 46 H 0.05 0.56 6.32 -51.93 -0.33 0.23 16 47 H 0.09 0.97 8.14 -51.93 -0.42 0.55 16 48 H 0.04 0.89 8.14 -51.93 -0.42 0.47 16 49 H 0.04 0.73 8.14 -51.92 -0.42 0.30 16 50 H 0.36 9.25 5.52 -40.82 -0.23 9.03 16 51 H 0.26 7.97 8.31 -40.82 -0.34 7.63 16 52 H 0.09 1.37 7.33 -51.93 -0.38 0.99 16 53 H 0.09 1.02 8.11 -51.93 -0.42 0.60 16 LS Contribution 394.42 15.07 5.94 5.94 Total: -1.00 -36.40 394.42 -6.72 -43.12 By element: Atomic # 1 Polarization: 14.23 SS G_CDS: -8.80 Total: 5.44 kcal Atomic # 6 Polarization: 1.33 SS G_CDS: 1.09 Total: 2.42 kcal Atomic # 7 Polarization: -56.30 SS G_CDS: 0.29 Total: -56.01 kcal Atomic # 8 Polarization: -18.88 SS G_CDS: -0.54 Total: -19.42 kcal Atomic # 14 Polarization: 23.22 SS G_CDS: -4.71 Total: 18.51 kcal Total LS contribution 5.94 Total: 5.94 kcal Total: -36.40 -6.72 -43.12 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. ZINC001201080021.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 -6.226 kcal (2) G-P(sol) polarization free energy of solvation -36.401 kcal (3) G-ENP(sol) elect.-nuc.-pol. free energy of system -42.626 kcal (4) G-CDS(sol) cavity-dispersion-solvent structure free energy -6.723 kcal (5) G-P-CDS(sol) = G-P(sol) + G-CDS(sol) = (2) + (4) -43.124 kcal (6) G-S(sol) free energy of system = (1) + (5) -49.349 kcal 1SCF WAS SPECIFIED, SO GEOMETRY OPTIMIZATION WAS NOT USED Total computer time = 1.08 seconds