Wall clock time and date at job start Sat Apr 11 2020 06:09:34 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 C 1.52996 * 109.51264 * 2 1 4 4 N 1.46511 * 109.46688 * 59.90282 * 3 2 1 5 5 C 1.34773 * 119.99660 * 179.97438 * 4 3 2 6 6 O 1.21277 * 120.00622 * 0.02562 * 5 4 3 7 7 C 1.50700 * 119.99676 * 179.97438 * 5 4 3 8 8 C 1.53001 * 109.47123 * 180.02562 * 7 5 4 9 9 C 1.53002 * 109.47200 * 180.02562 * 8 7 5 10 10 H 1.09004 * 109.52175 * 55.12067 * 9 8 7 11 11 C 1.53177 * 109.52412 * 175.31088 * 9 8 7 12 12 C 1.53176 * 109.04005 * 176.88622 * 11 9 8 13 13 C 1.53092 * 109.16116 * 303.02083 * 12 11 9 14 14 O 1.42906 * 109.41751 * 57.60530 * 13 12 11 15 15 C 1.42906 * 114.09576 * 298.84462 * 14 13 12 16 16 C 1.53305 * 109.51270 * 239.90524 * 2 1 3 17 17 N 1.47153 * 108.38117 * 292.32487 * 16 2 1 18 18 C 1.36937 * 120.98344 * 129.44135 * 17 16 2 19 19 H 1.08007 * 119.99657 * 333.66082 * 18 17 16 20 20 C 1.38812 * 120.00149 * 153.65563 * 18 17 16 21 21 N 1.31608 * 120.00041 * 334.96609 * 20 18 17 22 22 Si 1.86804 * 119.99759 * 154.96275 * 20 18 17 23 23 H 1.48504 * 109.47099 * 0.02562 * 22 20 18 24 24 H 1.48498 * 109.46978 * 240.00020 * 22 20 18 25 25 H 1.48505 * 109.47009 * 120.00164 * 22 20 18 26 26 C 1.47154 * 118.03853 * 309.71667 * 17 16 2 27 27 C 1.52877 * 108.60853 * 50.22658 * 26 17 16 28 28 O 1.42956 * 109.55848 * 120.13321 * 2 1 3 29 29 H 1.09004 * 109.46437 * 180.09922 * 1 2 3 30 30 H 1.08999 * 109.47094 * 300.09138 * 1 2 3 31 31 H 1.08996 * 109.47387 * 60.09713 * 1 2 3 32 32 H 1.09001 * 109.47227 * 179.89582 * 3 2 1 33 33 H 1.08993 * 109.47632 * 299.90450 * 3 2 1 34 34 H 0.96995 * 119.99679 * 359.97438 * 4 3 2 35 35 H 1.09001 * 109.47176 * 299.99504 * 7 5 4 36 36 H 1.08992 * 109.47344 * 60.00247 * 7 5 4 37 37 H 1.09001 * 109.46785 * 299.99328 * 8 7 5 38 38 H 1.08993 * 109.47655 * 59.99939 * 8 7 5 39 39 H 1.09003 * 109.54375 * 296.75567 * 11 9 8 40 40 H 1.09000 * 109.54430 * 57.01495 * 11 9 8 41 41 H 1.08998 * 109.52555 * 183.11189 * 12 11 9 42 42 H 1.08995 * 109.57161 * 62.95050 * 12 11 9 43 43 H 1.09003 * 109.48224 * 297.62744 * 13 12 11 44 44 H 1.08993 * 109.48522 * 177.59023 * 13 12 11 45 45 H 1.08994 * 109.48116 * 301.17132 * 15 14 13 46 46 H 1.09000 * 109.47949 * 181.13825 * 15 14 13 47 47 H 1.08995 * 109.68083 * 172.60201 * 16 2 1 48 48 H 1.08992 * 109.71074 * 52.06571 * 16 2 1 49 49 H 0.97006 * 120.00245 * 354.23460 * 21 20 18 50 50 H 1.08997 * 109.68421 * 170.01151 * 26 17 16 51 51 H 1.09000 * 109.68138 * 290.56353 * 26 17 16 52 52 H 1.08999 * 109.58500 * 187.92992 * 27 26 17 53 53 H 1.09004 * 109.41073 * 67.83873 * 27 26 17 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.0410 1.4421 0.0000 4 7 1.5512 2.1337 1.1951 5 6 1.8778 3.4237 1.4085 6 8 2.5767 4.0123 0.6110 7 6 1.3735 4.1353 2.6375 8 6 1.8846 5.5774 2.6376 9 6 1.3731 6.2996 3.8857 10 1 1.6610 5.7387 4.7748 11 6 1.9781 7.7055 3.9476 12 6 1.3993 8.4458 5.1573 13 6 -0.1273 8.4746 5.0455 14 8 -0.6228 7.1368 4.9608 15 6 -0.1524 6.4114 3.8230 16 6 2.0421 -0.7246 -1.2502 17 7 1.6971 -2.1504 -1.1342 18 6 1.0619 -2.8072 -2.1541 19 1 0.4395 -2.2580 -2.8452 20 6 1.2180 -4.1788 -2.2999 21 7 2.2940 -4.7724 -1.8287 22 14 -0.0993 -5.1820 -3.1648 23 1 -1.1998 -4.2844 -3.5989 24 1 -0.6339 -6.2060 -2.2317 25 1 0.4891 -5.8530 -4.3518 26 6 2.0492 -2.8480 0.1127 27 6 1.5417 -2.0214 1.2943 28 8 2.0086 -0.6762 1.1650 29 1 -0.3632 -1.0277 0.0018 30 1 -0.3633 0.5152 0.8892 31 1 -0.3634 0.5123 -0.8908 32 1 3.1311 1.4413 0.0019 33 1 1.6796 1.9557 -0.8907 34 1 0.9919 1.6631 1.8327 35 1 0.2835 4.1366 2.6350 36 1 1.7342 3.6215 3.5285 37 1 2.9746 5.5765 2.6400 38 1 1.5235 6.0910 1.7467 39 1 3.0609 7.6324 4.0488 40 1 1.7318 8.2490 3.0355 41 1 1.7823 9.4660 5.1776 42 1 1.6877 7.9301 6.0732 43 1 -0.4159 9.0263 4.1508 44 1 -0.5471 8.9626 5.9250 45 1 -0.4413 6.9362 2.9124 46 1 -0.5902 5.4132 3.8219 47 1 3.1240 -0.6133 -1.3217 48 1 1.5703 -0.3038 -2.1381 49 1 3.0156 -4.2429 -1.4547 50 1 1.5818 -3.8326 0.1271 51 1 3.1319 -2.9542 0.1808 52 1 1.9142 -2.4481 2.2256 53 1 0.4517 -2.0293 1.3004 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 ZINC001108312180.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:34 Heat of formation + Delta-G solvation = -127.220168 kcal Electronic energy + Delta-G solvation = -30466.675702 eV Core-core repulsion = 26338.196762 eV Total energy + Delta-G solvation = -4128.478939 eV No. of doubly occupied orbitals = 66 Molecular weight (most abundant/longest-lived isotopes) = 340.220 amu Computer time = 1.85 seconds Orbital eigenvalues (eV) -40.19230 -39.67920 -38.91365 -35.57983 -35.12382 -34.47597 -33.11715 -32.20159 -30.86373 -29.99977 -28.32818 -26.81565 -25.69847 -25.43680 -24.53714 -23.37486 -22.07030 -21.18390 -20.45593 -20.26324 -18.76748 -17.67149 -17.24065 -16.75032 -16.38316 -16.20440 -15.76447 -15.50053 -15.11613 -14.94177 -14.72340 -14.43281 -14.25179 -14.03292 -13.91303 -13.63924 -13.61597 -13.21118 -13.03280 -12.84650 -12.74555 -12.52371 -12.12357 -12.08300 -11.70804 -11.67985 -11.53396 -11.38907 -11.10255 -11.09299 -10.91776 -10.71112 -10.63050 -10.62665 -10.59078 -10.16177 -9.99729 -9.82297 -9.79766 -9.47285 -9.20621 -8.98388 -8.43395 -6.97345 -6.14987 -3.32481 2.65089 3.18819 3.45679 3.47855 3.49834 3.61260 4.07371 4.28424 4.34784 4.39647 4.49300 4.52061 4.65673 4.80128 4.89651 4.92543 4.95914 5.07104 5.11370 5.19979 5.24753 5.27948 5.29630 5.33391 5.37077 5.44427 5.46621 5.48591 5.57337 5.71624 5.74649 5.78533 5.86410 5.87341 5.89939 5.93874 6.01777 6.04503 6.18955 6.26514 6.31356 6.48785 6.54018 6.61237 6.77699 6.81000 6.99393 7.54267 7.68809 7.87669 8.24420 8.46093 8.97929 9.73395 10.16363 11.13102 Molecular weight = 340.22amu Principal moments of inertia in cm(-1) A = 0.030404 B = 0.001818 C = 0.001798 Principal moments of inertia in units of 10**(-40)*gram-cm**2 A = 920.697938 B =15401.753225 C =15573.295240 Net atomic charges, atomic populations, and dipole contributions calculated with CM2 Atom NO. Type Charge No. of electrons 1 C -0.165 4.165 2 C 0.088 3.912 3 C 0.138 3.862 4 N -0.719 5.719 5 C 0.509 3.491 6 O -0.548 6.548 7 C -0.144 4.144 8 C -0.096 4.096 9 C -0.124 4.124 10 H 0.081 0.919 11 C -0.118 4.118 12 C -0.157 4.157 13 C 0.059 3.941 14 O -0.384 6.384 15 C 0.062 3.938 16 C 0.130 3.870 17 N -0.596 5.596 18 C -0.350 4.350 19 H 0.059 0.941 20 C -0.024 4.024 21 N -0.915 5.915 22 Si 1.069 2.931 23 H -0.292 1.292 24 H -0.290 1.290 25 H -0.293 1.293 26 C 0.118 3.882 27 C 0.042 3.958 28 O -0.382 6.382 29 H 0.077 0.923 30 H 0.047 0.953 31 H 0.061 0.939 32 H 0.072 0.928 33 H 0.079 0.921 34 H 0.403 0.597 35 H 0.092 0.908 36 H 0.094 0.906 37 H 0.074 0.926 38 H 0.071 0.929 39 H 0.069 0.931 40 H 0.065 0.935 41 H 0.072 0.928 42 H 0.071 0.929 43 H 0.049 0.951 44 H 0.092 0.908 45 H 0.051 0.949 46 H 0.095 0.905 47 H 0.026 0.974 48 H 0.065 0.935 49 H 0.270 0.730 50 H 0.086 0.914 51 H 0.025 0.975 52 H 0.083 0.917 53 H 0.049 0.951 For the charges calculated by CM2: Dipole moment (debyes) X Y Z Total from point charges -4.422 22.674 14.092 27.060 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.222 4.222 2 C 0.047 3.953 3 C 0.014 3.986 4 N -0.370 5.370 5 C 0.295 3.705 6 O -0.426 6.426 7 C -0.183 4.183 8 C -0.134 4.134 9 C -0.143 4.143 10 H 0.099 0.901 11 C -0.155 4.155 12 C -0.194 4.194 13 C -0.018 4.018 14 O -0.302 6.302 15 C -0.015 4.015 16 C 0.003 3.997 17 N -0.321 5.321 18 C -0.479 4.479 19 H 0.077 0.923 20 C -0.226 4.226 21 N -0.555 5.555 22 Si 0.899 3.101 23 H -0.218 1.218 24 H -0.216 1.216 25 H -0.220 1.220 26 C -0.008 4.008 27 C -0.036 4.036 28 O -0.301 6.301 29 H 0.096 0.904 30 H 0.066 0.934 31 H 0.080 0.920 32 H 0.090 0.910 33 H 0.097 0.903 34 H 0.238 0.762 35 H 0.110 0.890 36 H 0.113 0.887 37 H 0.093 0.907 38 H 0.090 0.910 39 H 0.088 0.912 40 H 0.084 0.916 41 H 0.091 0.909 42 H 0.090 0.910 43 H 0.068 0.932 44 H 0.110 0.890 45 H 0.070 0.930 46 H 0.113 0.887 47 H 0.044 0.956 48 H 0.083 0.917 49 H 0.081 0.919 50 H 0.104 0.896 51 H 0.043 0.957 52 H 0.101 0.899 53 H 0.068 0.932 Dipole moment (debyes) X Y Z Total from point charges -4.397 21.529 13.316 25.694 hybrid contribution 0.339 1.175 0.278 1.254 sum -4.058 22.704 13.594 26.772 Atomic orbital electron populations 1.22235 0.95389 1.03671 1.00944 1.22199 0.93716 0.90480 0.88938 1.21333 0.99794 0.89770 0.87697 1.45987 1.53362 1.11806 1.25861 1.20669 0.80393 0.81813 0.87630 1.90696 1.38166 1.66308 1.47409 1.21745 1.02536 0.97056 0.97010 1.21310 1.01807 0.93015 0.97248 1.21850 0.96201 0.96792 0.99410 0.90086 1.21630 1.00598 0.95146 0.98169 1.22198 0.96565 1.01283 0.99388 1.22584 0.93424 0.84041 1.01763 1.87945 1.70381 1.25473 1.46416 1.22751 0.93252 0.96958 0.88520 1.21420 0.97270 0.86639 0.94371 1.44067 1.69503 1.05049 1.13474 1.20005 1.24822 0.93518 1.09544 0.92281 1.25860 1.00644 0.95458 1.00626 1.70721 1.04764 1.39429 1.40603 0.83472 0.75059 0.75506 0.76057 1.21798 1.21626 1.22038 1.21517 0.97629 0.96704 0.84974 1.23114 0.98536 0.82514 0.99387 1.87828 1.72886 1.26358 1.43064 0.90400 0.93433 0.92041 0.90967 0.90314 0.76158 0.88959 0.88741 0.90693 0.90984 0.91220 0.91644 0.90902 0.91023 0.93248 0.89019 0.93034 0.88674 0.95575 0.91674 0.91880 0.89583 0.95709 0.89926 0.93237 Number of geometries 1 Number of calculations of the screened coulomb radii 1 The total number of SCF iterations 16. 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.16 -2.55 7.92 37.16 0.29 -2.25 16 2 C 0.09 1.47 0.76 -90.46 -0.07 1.41 16 3 C 0.14 2.06 4.82 -4.04 -0.02 2.04 16 4 N -0.72 -9.01 5.15 -60.29 -0.31 -9.32 16 5 C 0.51 6.40 7.83 -10.99 -0.09 6.31 16 6 O -0.55 -8.99 15.98 5.56 0.09 -8.90 16 7 C -0.14 -1.20 5.46 -27.88 -0.15 -1.35 16 8 C -0.10 -0.80 3.98 -26.73 -0.11 -0.91 16 9 C -0.12 -0.88 2.39 -90.48 -0.22 -1.09 16 10 H 0.08 0.59 8.14 -51.93 -0.42 0.17 16 11 C -0.12 -0.77 5.31 -26.55 -0.14 -0.91 16 12 C -0.16 -0.92 6.14 -26.59 -0.16 -1.09 16 13 C 0.06 0.38 6.84 37.21 0.25 0.64 16 14 O -0.38 -3.43 10.77 -35.23 -0.38 -3.81 16 15 C 0.06 0.46 5.82 37.20 0.22 0.67 16 16 C 0.13 2.67 5.41 -3.53 -0.02 2.65 16 17 N -0.60 -14.33 2.44 -172.08 -0.42 -14.75 16 18 C -0.35 -9.40 9.54 -15.06 -0.14 -9.54 16 19 H 0.06 1.59 7.94 -52.48 -0.42 1.18 16 20 C -0.02 -0.65 4.71 -17.44 -0.08 -0.73 16 21 N -0.92 -26.80 12.44 55.49 0.69 -26.11 16 22 Si 1.07 24.64 29.99 -169.99 -5.10 19.54 16 23 H -0.29 -6.38 7.11 56.52 0.40 -5.98 16 24 H -0.29 -6.60 7.11 56.52 0.40 -6.20 16 25 H -0.29 -6.78 7.11 56.52 0.40 -6.38 16 26 C 0.12 2.66 4.85 -3.75 -0.02 2.64 16 27 C 0.04 0.77 6.58 37.09 0.24 1.01 16 28 O -0.38 -6.78 9.98 -35.23 -0.35 -7.13 16 29 H 0.08 1.40 6.11 -51.93 -0.32 1.08 16 30 H 0.05 0.61 6.86 -51.93 -0.36 0.25 16 31 H 0.06 0.92 8.14 -51.93 -0.42 0.49 16 32 H 0.07 1.18 8.14 -51.93 -0.42 0.76 16 33 H 0.08 1.16 8.14 -51.93 -0.42 0.73 16 34 H 0.40 4.42 7.20 -40.82 -0.29 4.12 16 35 H 0.09 0.60 6.63 -51.93 -0.34 0.26 16 36 H 0.09 0.67 8.14 -51.93 -0.42 0.25 16 37 H 0.07 0.72 8.10 -51.93 -0.42 0.30 16 38 H 0.07 0.67 8.10 -51.93 -0.42 0.25 16 39 H 0.07 0.45 8.14 -51.93 -0.42 0.02 16 40 H 0.06 0.43 8.14 -51.93 -0.42 0.01 16 41 H 0.07 0.34 8.14 -51.93 -0.42 -0.08 16 42 H 0.07 0.45 8.14 -51.93 -0.42 0.02 16 43 H 0.05 0.29 8.14 -51.93 -0.42 -0.13 16 44 H 0.09 0.48 8.14 -51.93 -0.42 0.06 16 45 H 0.05 0.37 8.14 -51.93 -0.42 -0.06 16 46 H 0.10 0.62 6.69 -51.93 -0.35 0.28 16 47 H 0.03 0.56 8.14 -51.93 -0.42 0.14 16 48 H 0.06 1.31 7.95 -51.93 -0.41 0.90 16 49 H 0.27 8.01 6.61 -40.82 -0.27 7.74 16 50 H 0.09 2.15 6.44 -51.93 -0.33 1.81 16 51 H 0.02 0.59 7.76 -51.93 -0.40 0.19 16 52 H 0.08 1.31 8.14 -51.93 -0.42 0.89 16 53 H 0.05 0.88 6.25 -51.93 -0.32 0.55 16 LS Contribution 403.06 15.07 6.07 6.07 Total: -1.00 -32.01 403.06 -9.29 -41.29 By element: Atomic # 1 Polarization: 13.00 SS G_CDS: -9.38 Total: 3.63 kcal Atomic # 6 Polarization: -0.30 SS G_CDS: -0.21 Total: -0.51 kcal Atomic # 7 Polarization: -50.14 SS G_CDS: -0.04 Total: -50.18 kcal Atomic # 8 Polarization: -19.20 SS G_CDS: -0.64 Total: -19.84 kcal Atomic # 14 Polarization: 24.64 SS G_CDS: -5.10 Total: 19.54 kcal Total LS contribution 6.07 Total: 6.07 kcal Total: -32.01 -9.29 -41.29 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. ZINC001108312180.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 -85.926 kcal (2) G-P(sol) polarization free energy of solvation -32.006 kcal (3) G-ENP(sol) elect.-nuc.-pol. free energy of system -117.932 kcal (4) G-CDS(sol) cavity-dispersion-solvent structure free energy -9.288 kcal (5) G-P-CDS(sol) = G-P(sol) + G-CDS(sol) = (2) + (4) -41.294 kcal (6) G-S(sol) free energy of system = (1) + (5) -127.220 kcal 1SCF WAS SPECIFIED, SO GEOMETRY OPTIMIZATION WAS NOT USED Total computer time = 1.85 seconds