Wall clock time and date at job start Mon Mar 8 2021 14:58:14 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= 0 * AM1 - THE AM1 HAMILTONIAN TO BE USED * SM5.42R - SM5.42R CALCULATIONS WILL BE PERFORMED * - CHARGE MODEL 2 WILL BE USED (CM2A) * - THE SOLVENT IS WATER ******************************************************************************* 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.53005 * 1 3 3 O 1.45200 * 109.47163 * 2 1 4 4 C 1.34229 * 117.00003 * 179.97438 * 3 2 1 5 5 O 1.20820 * 119.99829 * 0.02562 * 4 3 2 6 6 C 1.50698 * 119.99999 * 180.02562 * 4 3 2 7 7 H 1.09005 * 112.84737 * 311.19880 * 6 4 3 8 8 C 1.53784 * 113.61692 * 82.50275 * 6 4 3 9 9 C 1.53784 * 87.07950 * 139.98073 * 8 6 4 10 10 F 1.39897 * 113.68989 * 220.00729 * 9 8 6 11 11 C 1.53003 * 113.61054 * 89.11734 * 9 8 6 12 12 S 1.81401 * 109.46957 * 82.50518 * 11 9 8 13 13 O 1.42092 * 110.54713 * 68.36801 * 12 11 9 14 14 O 1.42104 * 110.53680 * 291.62050 * 12 11 9 15 15 N 1.65598 * 104.44908 * 179.97438 * 12 11 9 16 16 C 1.47083 * 120.88301 * 89.99597 * 15 12 11 17 17 C 1.53237 * 108.52218 * 231.21097 * 16 15 12 18 18 N 1.46905 * 109.33607 * 306.16557 * 17 16 15 19 19 C 1.39614 * 110.99730 * 187.25406 * 18 17 16 20 20 C 1.34101 * 125.23742 * 180.02562 * 19 18 17 21 21 C 1.38580 * 114.72268 * 179.71843 * 20 19 18 22 22 C 1.33228 * 114.98079 * 0.49379 * 21 20 19 23 23 C 1.50694 * 125.12978 * 179.79447 * 22 21 20 24 24 O 1.42901 * 109.47118 * 64.71977 * 23 22 21 25 25 Si 1.86298 * 109.47775 * 184.72222 * 23 22 21 26 26 S 1.75855 * 109.74796 * 359.53363 * 22 21 20 27 27 C 1.46903 * 111.26793 * 63.10872 * 18 17 16 28 28 C 1.47074 * 120.88211 * 270.01993 * 15 12 11 29 29 C 1.53779 * 113.61452 * 180.02562 * 6 4 3 30 30 H 1.08998 * 109.47224 * 179.97438 * 1 2 3 31 31 H 1.08997 * 109.46814 * 299.99777 * 1 2 3 32 32 H 1.08991 * 109.47211 * 59.99465 * 1 2 3 33 33 H 1.08997 * 109.46814 * 240.00223 * 2 1 3 34 34 H 1.08998 * 109.47476 * 120.00251 * 2 1 3 35 35 H 1.09002 * 113.60778 * 25.43203 * 8 6 4 36 36 H 1.08998 * 113.61324 * 254.53364 * 8 6 4 37 37 H 1.08995 * 109.47203 * 322.50168 * 11 9 8 38 38 H 1.09005 * 109.46511 * 202.49921 * 11 9 8 39 39 H 1.08996 * 109.62998 * 111.49561 * 16 15 12 40 40 H 1.09005 * 109.78671 * 351.01455 * 16 15 12 41 41 H 1.09006 * 109.49277 * 186.21382 * 17 16 15 42 42 H 1.08989 * 109.51912 * 66.14191 * 17 16 15 43 43 H 1.08005 * 122.63965 * 0.02562 * 20 19 18 44 44 H 1.07997 * 122.51420 * 180.28730 * 21 20 19 45 45 H 1.09004 * 109.47407 * 304.72334 * 23 22 21 46 46 H 0.96708 * 114.00246 * 180.02562 * 24 23 22 47 47 H 1.48506 * 109.99834 * 300.00486 * 25 23 22 48 48 H 1.48502 * 110.00011 * 61.31751 * 25 23 22 49 49 H 1.08999 * 109.49223 * 56.85452 * 27 18 17 50 50 H 1.09000 * 109.49292 * 176.93384 * 27 18 17 51 51 H 1.09000 * 109.62900 * 248.50495 * 28 15 12 52 52 H 1.09000 * 109.62897 * 9.08115 * 28 15 12 53 53 H 1.08995 * 113.61845 * 105.47101 * 29 6 4 54 54 H 1.08999 * 113.61071 * 334.56167 * 29 6 4 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.5301 0.0000 0.0000 3 8 2.0141 1.3690 0.0000 4 6 3.3448 1.5448 0.0005 5 8 4.0807 0.5866 0.0005 6 6 3.9208 2.9374 0.0012 7 1 3.4714 3.5816 -0.7547 8 6 3.9859 3.5765 1.3984 9 6 5.3562 4.1422 0.9895 10 9 6.3100 4.0988 2.0121 11 6 5.2763 5.5103 0.3093 12 16 5.1599 6.8075 1.5720 13 8 6.3776 6.9080 2.2973 14 8 3.9101 6.7378 2.2446 15 7 5.0816 8.2185 0.7088 16 6 3.7869 8.7753 0.2880 17 6 3.8487 9.0696 -1.2146 18 7 5.0121 9.9214 -1.4959 19 6 5.0115 10.3456 -2.8260 20 6 5.9407 11.1328 -3.3875 21 6 5.7288 11.4144 -4.7277 22 6 4.6310 10.8696 -5.2504 23 6 4.1794 11.0296 -6.6791 24 8 5.1346 10.4223 -7.5515 25 14 2.5229 10.2070 -6.9029 26 16 3.7923 9.9294 -4.0235 27 6 6.2677 9.2356 -1.1626 28 6 6.3049 8.9480 0.3421 29 6 5.4583 2.9680 0.0011 30 1 -0.3633 -1.0276 -0.0005 31 1 -0.3633 0.5138 0.8900 32 1 -0.3633 0.5139 -0.8899 33 1 1.8933 -0.5138 -0.8900 34 1 1.8934 -0.5139 0.8899 35 1 4.0377 2.8496 2.2089 36 1 3.2287 4.3430 1.5635 37 1 4.3952 5.5467 -0.3313 38 1 6.1706 5.6691 -0.2934 39 1 3.5902 9.6979 0.8341 40 1 2.9932 8.0553 0.4878 41 1 2.9383 9.5841 -1.5222 42 1 3.9415 8.1337 -1.7654 43 1 6.7901 11.5133 -2.8396 44 1 6.4012 12.0320 -5.3046 45 1 4.0934 12.0900 -6.9165 46 1 4.9081 10.4877 -8.4894 47 1 2.6268 8.7606 -6.5824 48 1 1.5134 10.8506 -6.0243 49 1 6.3270 8.2971 -1.7139 50 1 7.1116 9.8704 -1.4326 51 1 6.3503 9.8872 0.8933 52 1 7.1783 8.3397 0.5771 53 1 5.8895 3.2435 -0.9613 54 1 5.9128 2.0748 0.4297 There are 78 doubly filled levels RHF calculation, no. of doubly occupied orbitals= 78 No. of singly occupied orbitals= 1 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) F: (AM1): M.J.S. DEWAR AND E. G. ZOEBISCH, THEOCHEM, 180, 1 (1988). Si: (AM1): M.J.S.DEWAR, C. JIE, ORGANOMETALLICS, 6, 1486-1490 (1987). S: (AM1): M.J.S. DEWAR, Y.-C. YUAN, INORGANIC CHEM., 29, 589 (1990) 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=0 AM1 1SCF TLIMIT=15 GEO-OK SM5.42R & SOLVNT=WATER s_40_80748_15494190.mol2 54 1SCF WAS SPECIFIED, SO GEOMETRY OPTIMIZATION WAS NOT USED SCF FIELD WAS ACHIEVED Wall clock time and date at job start Mon Mar 8 2021 14:58:14 Heat of formation + Delta-G solvation = -146.224532 kcal Electronic energy + Delta-G solvation = -41455.942431 eV Core-core repulsion = 35950.646883 eV Total energy + Delta-G solvation = -5505.295548 eV No. of doubly occupied orbitals = 78 and no. of open levels = 1 Molecular weight (most abundant/longest-lived isotopes) = 449.127 amu Computer time = 1.81 seconds Orbital eigenvalues (eV) -49.13535 -41.98105 -40.77928 -39.64042 -38.93340 -38.45311 -37.67299 -37.46625 -37.20636 -35.62848 -32.80313 -32.49838 -31.49668 -31.00869 -28.18411 -27.37591 -27.10209 -25.61795 -25.15431 -24.06327 -23.58164 -22.52154 -22.04111 -20.57510 -20.33974 -20.04680 -19.42201 -19.10852 -18.47457 -18.26738 -18.15037 -17.84643 -17.55772 -17.35521 -17.13979 -16.86613 -16.53940 -16.31657 -16.07905 -15.92032 -15.84762 -15.74589 -15.50094 -15.23223 -15.18915 -14.92638 -14.54313 -14.50323 -14.22400 -14.19670 -13.83781 -13.73648 -13.73004 -13.70021 -13.54950 -13.39113 -13.23002 -13.14267 -13.07496 -12.97477 -12.80339 -12.67822 -12.55265 -12.54160 -12.52851 -12.49712 -12.29736 -12.13520 -11.99797 -11.86574 -11.63358 -11.58240 -11.50173 -11.00423 -10.33746 -10.18725 -9.60230 -8.87082 -5.19879 -1.09588 -0.14447 0.34525 0.70082 0.92119 1.23223 1.33416 1.48563 1.54957 1.71011 1.90851 2.43827 2.48488 2.59906 2.78745 2.86602 2.92988 3.09617 3.25588 3.37625 3.49898 3.65821 3.74655 3.84011 3.91053 3.92048 3.93436 3.94488 3.98423 4.01340 4.09785 4.19393 4.26573 4.28084 4.34058 4.36066 4.42052 4.43978 4.49903 4.51697 4.53319 4.61791 4.66565 4.76821 4.78061 4.81338 4.85285 4.89972 4.94950 4.95909 5.08230 5.15513 5.27296 5.34853 5.38658 5.63278 5.76357 6.53021 6.87816 Net atomic charges, atomic populations, and dipole contributions calculated with CM2 Atom NO. Type Charge No. of electrons 1 C -0.149 4.149 2 C 0.053 3.947 3 O -0.357 6.357 4 C 0.472 3.528 5 O -0.535 6.535 6 C -0.099 4.099 7 H 0.172 0.828 8 C -0.154 4.154 9 C 0.089 3.911 10 F -0.204 7.204 11 C -0.632 4.632 12 S 2.610 3.390 13 O -0.975 6.975 14 O -0.969 6.969 15 N -0.998 5.998 16 C 0.126 3.874 17 C 0.024 3.976 18 N -0.453 5.453 19 C 0.029 3.971 20 C -0.183 4.183 21 C -0.079 4.079 22 C -0.211 4.211 23 C -0.036 4.036 24 O -0.562 6.562 25 Si 0.749 3.251 26 S 0.034 5.966 27 C 0.040 3.960 28 C 0.121 3.879 29 C -0.120 4.120 30 H 0.096 0.904 31 H 0.068 0.932 32 H 0.071 0.929 33 H 0.070 0.930 34 H 0.066 0.934 35 H 0.103 0.897 36 H 0.122 0.878 37 H 0.193 0.807 38 H 0.179 0.821 39 H 0.085 0.915 40 H 0.118 0.882 41 H 0.097 0.903 42 H 0.071 0.929 43 H 0.158 0.842 44 H 0.154 0.846 45 H 0.101 0.899 46 H 0.404 0.596 47 H -0.280 1.280 48 H -0.275 1.275 49 H 0.065 0.935 50 H 0.109 0.891 51 H 0.086 0.914 52 H 0.111 0.889 53 H 0.124 0.876 54 H 0.099 0.901 For the charges calculated by CM2: Dipole moment (debyes) X Y Z Total from point charges -3.943 2.743 -9.407 10.562 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.206 4.206 2 C -0.021 4.021 3 O -0.273 6.273 4 C 0.316 3.684 5 O -0.425 6.425 6 C -0.117 4.117 7 H 0.190 0.810 8 C -0.191 4.191 9 C 0.069 3.931 10 F -0.184 7.184 11 C -0.762 4.762 12 S 2.706 3.294 13 O -0.965 6.965 14 O -0.959 6.959 15 N -0.837 5.837 16 C 0.000 4.000 17 C -0.105 4.105 18 N -0.170 5.170 19 C -0.185 4.185 20 C -0.215 4.215 21 C -0.110 4.110 22 C -0.328 4.328 23 C -0.148 4.148 24 O -0.372 6.372 25 Si 0.653 3.347 26 S 0.293 5.707 27 C -0.087 4.087 28 C -0.004 4.004 29 C -0.159 4.159 30 H 0.114 0.886 31 H 0.087 0.913 32 H 0.090 0.910 33 H 0.089 0.911 34 H 0.085 0.915 35 H 0.121 0.879 36 H 0.140 0.860 37 H 0.210 0.790 38 H 0.197 0.803 39 H 0.103 0.897 40 H 0.136 0.864 41 H 0.116 0.884 42 H 0.089 0.911 43 H 0.176 0.824 44 H 0.171 0.829 45 H 0.118 0.882 46 H 0.256 0.744 47 H -0.206 1.206 48 H -0.201 1.201 49 H 0.084 0.916 50 H 0.127 0.873 51 H 0.104 0.896 52 H 0.129 0.871 53 H 0.142 0.858 54 H 0.118 0.882 Dipole moment (debyes) X Y Z Total from point charges -5.444 1.319 -7.067 9.018 hybrid contribution 1.856 -0.056 -0.426 1.905 sum -3.587 1.263 -7.493 8.403 Atomic orbital electron populations 1.21831 0.90607 1.04313 1.03829 1.23290 0.96993 0.79180 1.02592 1.86516 1.23906 1.33632 1.83271 1.23745 0.80749 0.92599 0.71273 1.90691 1.61194 1.44091 1.46478 1.22462 0.97005 0.90909 1.01363 0.81038 1.23902 0.96256 1.02219 0.96771 1.21920 0.88363 0.95275 0.87564 1.92960 1.66699 1.97425 1.61325 1.31032 1.17080 1.08859 1.19193 1.04936 0.73505 0.75036 0.75962 1.93596 1.47978 1.85584 1.69353 1.93587 1.44490 1.86373 1.71421 1.57776 1.13080 1.38557 1.74257 1.21643 0.83993 1.01654 0.92676 1.22891 0.93139 0.96673 0.97747 1.58890 1.05191 1.39609 1.13290 1.23771 0.99985 1.03388 0.91366 1.20165 1.02860 1.03204 0.95279 1.20119 0.96368 1.00026 0.94469 1.23509 1.01946 1.08666 0.98644 1.24014 1.02258 1.02850 0.85723 1.86457 1.48644 1.80422 1.21679 0.94091 0.74307 0.77409 0.88869 1.84413 1.31205 1.52394 1.02707 1.22892 0.90679 0.98942 0.96206 1.21721 0.85971 0.98732 0.94017 1.23609 0.93864 0.97035 1.01366 0.88566 0.91297 0.91027 0.91110 0.91511 0.87872 0.85982 0.78976 0.80321 0.89679 0.86353 0.88446 0.91127 0.82425 0.82876 0.88174 0.74449 1.20608 1.20088 0.91638 0.87257 0.89608 0.87081 0.85761 0.88234 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 0.61 10.29 71.98 0.74 1.36 16 2 C 0.05 0.02 6.19 71.98 0.45 0.46 16 3 O -0.36 -0.79 10.22 -51.41 -0.53 -1.32 16 4 C 0.47 1.39 7.48 71.24 0.53 1.93 16 5 O -0.54 -5.16 15.65 21.98 0.34 -4.82 16 6 C -0.10 0.43 3.91 -11.01 -0.04 0.38 16 7 H 0.17 -2.09 7.66 -2.38 -0.02 -2.11 16 8 C -0.15 0.33 7.09 31.12 0.22 0.55 16 9 C 0.09 -0.15 1.79 -10.25 -0.02 -0.16 16 10 F -0.20 -1.33 16.31 44.97 0.73 -0.60 16 11 C -0.63 5.09 5.76 71.98 0.41 5.50 16 12 S 2.61 1.16 5.14 -56.49 -0.29 0.87 16 13 O -0.97 -8.74 16.23 -127.65 -2.07 -10.82 16 14 O -0.97 -5.29 16.22 -127.65 -2.07 -7.36 16 15 N -1.00 3.71 3.26 -512.17 -1.67 2.04 16 16 C 0.13 -0.66 6.46 86.35 0.56 -0.11 16 17 C 0.02 -0.12 4.70 86.38 0.41 0.28 16 18 N -0.45 1.04 4.71 -698.00 -3.29 -2.25 16 19 C 0.03 0.02 6.03 82.51 0.50 0.52 16 20 C -0.18 -0.12 9.48 21.41 0.20 0.08 16 21 C -0.08 -0.14 10.04 21.22 0.21 0.07 16 22 C -0.21 -0.60 5.72 24.55 0.14 -0.46 16 23 C -0.04 0.07 2.95 71.23 0.21 0.28 16 24 O -0.56 1.12 12.86 -148.98 -1.92 -0.80 16 25 Si 0.75 -2.19 30.97 68.60 2.12 -0.06 16 26 S 0.03 0.09 18.70 -56.49 -1.06 -0.97 16 27 C 0.04 -0.24 5.47 86.38 0.47 0.24 16 28 C 0.12 -0.59 6.46 86.35 0.56 -0.03 16 29 C -0.12 0.40 6.89 31.13 0.21 0.61 16 30 H 0.10 -0.65 8.14 -2.39 -0.02 -0.67 16 31 H 0.07 -0.26 8.14 -2.39 -0.02 -0.28 16 32 H 0.07 -0.32 8.14 -2.39 -0.02 -0.34 16 33 H 0.07 0.03 8.13 -2.39 -0.02 0.01 16 34 H 0.07 0.08 8.13 -2.39 -0.02 0.06 16 35 H 0.10 0.01 8.09 -2.39 -0.02 -0.01 16 36 H 0.12 -0.56 7.51 -2.39 -0.02 -0.58 16 37 H 0.19 -3.12 7.34 -2.39 -0.02 -3.14 16 38 H 0.18 -2.36 8.14 -2.39 -0.02 -2.38 16 39 H 0.08 -0.33 8.14 -2.39 -0.02 -0.35 16 40 H 0.12 -0.99 7.20 -2.38 -0.02 -1.00 16 41 H 0.10 -0.45 7.43 -2.38 -0.02 -0.46 16 42 H 0.07 -0.55 8.10 -2.39 -0.02 -0.57 16 43 H 0.16 -0.41 7.58 -2.91 -0.02 -0.43 16 44 H 0.15 -0.12 8.06 -2.91 -0.02 -0.14 16 45 H 0.10 -0.59 7.96 -2.38 -0.02 -0.61 16 46 H 0.40 -5.05 8.74 -74.05 -0.65 -5.69 16 47 H -0.28 -1.07 7.11 99.48 0.71 -0.36 16 48 H -0.27 -0.69 7.11 99.48 0.71 0.01 16 49 H 0.07 -0.52 8.14 -2.39 -0.02 -0.54 16 50 H 0.11 -0.80 7.14 -2.39 -0.02 -0.81 16 51 H 0.09 -0.33 8.14 -2.39 -0.02 -0.34 16 52 H 0.11 -0.69 7.19 -2.39 -0.02 -0.71 16 53 H 0.12 -1.07 8.10 -2.39 -0.02 -1.09 16 54 H 0.10 0.14 7.43 -2.39 -0.02 0.12 16 Total: 0.00 -33.41 460.01 -3.59 -37.00 By element: Atomic # 1 Polarization: -22.76 SS G_CDS: 0.33 Total: -22.43 kcal Atomic # 6 Polarization: 5.74 SS G_CDS: 5.77 Total: 11.50 kcal Atomic # 7 Polarization: 4.76 SS G_CDS: -4.96 Total: -0.20 kcal Atomic # 8 Polarization: -18.87 SS G_CDS: -6.24 Total: -25.11 kcal Atomic # 9 Polarization: -1.33 SS G_CDS: 0.73 Total: -0.60 kcal Atomic # 14 Polarization: -2.19 SS G_CDS: 2.12 Total: -0.06 kcal Atomic # 16 Polarization: 1.25 SS G_CDS: -1.35 Total: -0.10 kcal Total: -33.41 -3.59 -37.00 kcal The number of atoms in the molecule is 54 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. s_40_80748_15494190.mol2 54 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 -109.223 kcal (2) G-P(sol) polarization free energy of solvation -33.410 kcal (3) G-ENP(sol) elect.-nuc.-pol. free energy of system -142.633 kcal (4) G-CDS(sol) cavity-dispersion-solvent structure free energy -3.592 kcal (5) G-P-CDS(sol) = G-P(sol) + G-CDS(sol) = (2) + (4) -37.002 kcal (6) G-S(sol) free energy of system = (1) + (5) -146.225 kcal 1SCF WAS SPECIFIED, SO GEOMETRY OPTIMIZATION WAS NOT USED Total computer time = 1.81 seconds