Wall clock time and date at job start Tue Mar 31 2020 12:20:19 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.52997 * 1 3 3 C 1.52999 * 109.47317 * 2 1 4 4 C 1.53006 * 109.47127 * 240.00489 * 2 1 3 5 5 H 1.08996 * 109.46945 * 304.99450 * 4 2 1 6 6 N 1.46901 * 109.47113 * 184.99881 * 4 2 1 7 7 C 1.46900 * 110.99814 * 84.99689 * 6 4 2 8 8 C 1.53000 * 109.47049 * 179.97438 * 7 6 4 9 9 H 1.09001 * 109.47441 * 60.00162 * 8 7 6 10 10 O 1.42901 * 109.47075 * 300.00347 * 8 7 6 11 11 C 1.50696 * 109.47499 * 180.02562 * 8 7 6 12 12 H 1.07998 * 120.00064 * 0.02562 * 11 8 7 13 13 C 1.37880 * 120.00541 * 180.02562 * 11 8 7 14 14 N 1.31512 * 119.99703 * 359.97438 * 13 11 8 15 15 Si 1.86797 * 120.00333 * 180.02562 * 13 11 8 16 16 H 1.48499 * 109.47150 * 59.99676 * 15 13 11 17 17 H 1.48499 * 109.47105 * 179.97438 * 15 13 11 18 18 H 1.48500 * 109.47152 * 299.99804 * 15 13 11 19 19 C 1.50695 * 109.47196 * 64.99596 * 4 2 1 20 20 C 1.38243 * 119.99850 * 79.73089 * 19 4 2 21 21 C 1.38238 * 119.99755 * 179.78734 * 20 19 4 22 22 C 1.38230 * 119.99838 * 0.46249 * 21 20 19 23 23 C 1.50704 * 119.99608 * 179.79125 * 22 21 20 24 24 C 1.53002 * 109.47221 * 239.97507 * 23 22 21 25 25 C 1.52996 * 109.47349 * 359.97438 * 23 22 21 26 26 C 1.53003 * 109.47002 * 119.97793 * 23 22 21 27 27 C 1.38230 * 120.00408 * 359.76538 * 22 21 20 28 28 C 1.38229 * 120.00234 * 260.00222 * 19 4 2 29 29 H 1.09002 * 109.47153 * 299.99641 * 1 2 3 30 30 H 1.09001 * 109.46832 * 179.97438 * 1 2 3 31 31 H 1.09004 * 109.47261 * 59.99813 * 1 2 3 32 32 H 1.08996 * 109.47853 * 120.00331 * 2 1 3 33 33 H 1.09002 * 109.46716 * 60.00091 * 3 2 1 34 34 H 1.09000 * 109.47486 * 180.02562 * 3 2 1 35 35 H 1.09001 * 109.47362 * 300.00732 * 3 2 1 36 36 H 1.00894 * 111.00080 * 208.95076 * 6 4 2 37 37 H 1.08999 * 109.47444 * 299.99825 * 7 6 4 38 38 H 1.08999 * 109.47319 * 60.00032 * 7 6 4 39 39 H 0.96703 * 114.00207 * 299.99647 * 10 8 7 40 40 H 0.97004 * 120.00229 * 180.02562 * 14 13 11 41 41 H 1.07997 * 119.99964 * 359.97438 * 20 19 4 42 42 H 1.08002 * 119.99574 * 180.25017 * 21 20 19 43 43 H 1.09007 * 109.46763 * 59.99727 * 24 23 22 44 44 H 1.08993 * 109.47256 * 179.97438 * 24 23 22 45 45 H 1.09002 * 109.46933 * 300.00370 * 24 23 22 46 46 H 1.09002 * 109.46878 * 59.99975 * 25 23 22 47 47 H 1.08996 * 109.47239 * 179.97438 * 25 23 22 48 48 H 1.08998 * 109.47141 * 300.00574 * 25 23 22 49 49 H 1.08993 * 109.47255 * 59.99724 * 26 23 22 50 50 H 1.09006 * 109.46717 * 180.02562 * 26 23 22 51 51 H 1.09002 * 109.47018 * 299.99898 * 26 23 22 52 52 H 1.08004 * 120.00146 * 180.02562 * 27 22 21 53 53 H 1.08002 * 120.00177 * 359.97438 * 28 19 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.5300 0.0000 0.0000 3 6 2.0400 1.4425 0.0000 4 6 2.0400 -0.7212 -1.2494 5 1 1.6055 -1.7197 -1.2953 6 7 3.5040 -0.8266 -1.1905 7 6 3.9163 -2.0062 -0.4181 8 6 5.4439 -2.0823 -0.3798 9 1 5.8409 -1.1829 0.0910 10 8 5.9493 -2.1850 -1.7125 11 6 5.8669 -3.2927 0.4120 12 1 5.1208 -3.9429 0.8443 13 6 7.2064 -3.5701 0.5846 14 7 8.1150 -2.7786 0.0578 15 14 7.7307 -5.0708 1.5655 16 1 7.1906 -6.2957 0.9228 17 1 9.2134 -5.1446 1.6028 18 1 7.2049 -4.9655 2.9503 19 6 1.6413 0.0571 -2.4766 20 6 2.3695 1.1690 -2.8568 21 6 2.0064 1.8803 -3.9851 22 6 0.9093 1.4856 -4.7277 23 6 0.5107 2.2639 -5.9550 24 6 -0.9110 2.8012 -5.7788 25 6 1.4769 3.4334 -6.1535 26 6 0.5612 1.3472 -7.1790 27 6 0.1783 0.3766 -4.3449 28 6 0.5440 -0.3373 -3.2190 29 1 -0.3633 0.5138 0.8900 30 1 -0.3633 -1.0277 -0.0005 31 1 -0.3634 0.5139 -0.8900 32 1 1.8934 -0.5138 0.8899 33 1 1.6767 1.9563 0.8900 34 1 3.1300 1.4425 -0.0005 35 1 1.6766 1.9564 -0.8899 36 1 3.9028 -0.8437 -2.1172 37 1 3.5304 -1.9280 0.5983 38 1 3.5198 -2.9056 -0.8891 39 1 5.6368 -2.9640 -2.1927 40 1 9.0574 -2.9739 0.1789 41 1 3.2240 1.4800 -2.2741 42 1 2.5778 2.7466 -4.2844 43 1 -0.9470 3.4543 -4.9067 44 1 -1.1991 3.3645 -6.6662 45 1 -1.5994 1.9680 -5.6373 46 1 2.4897 3.0506 -6.2791 47 1 1.1888 3.9961 -7.0414 48 1 1.4409 4.0864 -5.2815 49 1 -0.1271 0.5140 -7.0376 50 1 0.2733 1.9103 -8.0668 51 1 1.5741 0.9644 -7.3046 52 1 -0.6786 0.0679 -4.9253 53 1 -0.0274 -1.2036 -2.9197 RHF calculation, no. of doubly occupied orbitals= 62 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 ZINC000929859878.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 12:20:19 Heat of formation + Delta-G solvation = -54.826235 kcal Electronic energy + Delta-G solvation = -27791.264335 eV Core-core repulsion = 24240.488787 eV Total energy + Delta-G solvation = -3550.775548 eV No. of doubly occupied orbitals = 62 Molecular weight (most abundant/longest-lived isotopes) = 319.225 amu Computer time = 0.73 seconds Orbital eigenvalues (eV) -40.21978 -38.13083 -36.15996 -35.16666 -32.79990 -32.06861 -31.11484 -30.10977 -27.91896 -27.64485 -27.45559 -26.83161 -25.36158 -23.50202 -22.40918 -21.98644 -20.46729 -19.37582 -19.07425 -16.89645 -16.59823 -15.54632 -15.45827 -14.96843 -14.76283 -14.55242 -14.43007 -14.30724 -13.96808 -13.74825 -13.66788 -13.50352 -13.21559 -12.86091 -12.76508 -12.59226 -12.21374 -12.08020 -12.02581 -11.88870 -11.81058 -11.64280 -11.43022 -11.19954 -11.10100 -11.05879 -10.87668 -10.82653 -10.79240 -10.61842 -10.59581 -10.45364 -10.15936 -10.00237 -8.99665 -8.74456 -8.55355 -8.40595 -8.13423 -7.74954 -5.98576 -4.16825 1.43600 1.53740 3.29711 3.35108 3.36734 3.84017 4.37675 4.44767 4.57620 4.68260 4.81118 4.87526 5.06480 5.12967 5.16274 5.23609 5.28472 5.31150 5.36456 5.39688 5.52842 5.54982 5.57789 5.62183 5.66002 5.69164 5.79228 5.82372 5.85172 5.86934 5.99045 6.07177 6.09582 6.09880 6.10775 6.14327 6.18894 6.21024 6.30639 6.37410 6.45930 6.53573 6.63100 6.69224 6.70976 6.79563 6.92750 7.12952 7.59552 7.78075 7.89210 8.04014 8.26040 8.86569 9.40186 9.60250 11.05601 Molecular weight = 319.23amu Principal moments of inertia in cm(-1) A = 0.018159 B = 0.002823 C = 0.002622 Principal moments of inertia in units of 10**(-40)*gram-cm**2 A = 1541.553851 B = 9916.027052 C =10677.586443 Net atomic charges, atomic populations, and dipole contributions calculated with CM2 Atom NO. Type Charge No. of electrons 1 C -0.147 4.147 2 C -0.079 4.079 3 C -0.142 4.142 4 C 0.129 3.871 5 H 0.042 0.958 6 N -0.677 5.677 7 C 0.037 3.963 8 C 0.233 3.767 9 H 0.063 0.937 10 O -0.565 6.565 11 C -0.573 4.573 12 H 0.053 0.947 13 C 0.080 3.920 14 N -0.858 5.858 15 Si 0.889 3.111 16 H -0.277 1.277 17 H -0.286 1.286 18 H -0.276 1.276 19 C -0.079 4.079 20 C -0.092 4.092 21 C -0.126 4.126 22 C -0.077 4.077 23 C -0.029 4.029 24 C -0.137 4.137 25 C -0.144 4.144 26 C -0.138 4.138 27 C -0.121 4.121 28 C -0.113 4.113 29 H 0.051 0.949 30 H 0.048 0.952 31 H 0.055 0.945 32 H 0.072 0.928 33 H 0.040 0.960 34 H 0.070 0.930 35 H 0.054 0.946 36 H 0.345 0.655 37 H 0.064 0.936 38 H 0.015 0.985 39 H 0.340 0.660 40 H 0.261 0.739 41 H 0.134 0.866 42 H 0.123 0.877 43 H 0.055 0.945 44 H 0.055 0.945 45 H 0.055 0.945 46 H 0.060 0.940 47 H 0.055 0.945 48 H 0.058 0.942 49 H 0.055 0.945 50 H 0.055 0.945 51 H 0.057 0.943 52 H 0.114 0.886 53 H 0.115 0.885 For the charges calculated by CM2: Dipole moment (debyes) X Y Z Total from point charges -18.405 10.621 -10.568 23.732 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.205 4.205 2 C -0.098 4.098 3 C -0.199 4.199 4 C 0.019 3.981 5 H 0.060 0.940 6 N -0.309 5.309 7 C -0.092 4.092 8 C 0.176 3.824 9 H 0.081 0.919 10 O -0.367 6.367 11 C -0.612 4.612 12 H 0.072 0.928 13 C -0.127 4.127 14 N -0.494 5.494 15 Si 0.716 3.284 16 H -0.202 1.202 17 H -0.212 1.212 18 H -0.202 1.202 19 C -0.079 4.079 20 C -0.110 4.110 21 C -0.144 4.144 22 C -0.077 4.077 23 C -0.029 4.029 24 C -0.194 4.194 25 C -0.201 4.201 26 C -0.195 4.195 27 C -0.139 4.139 28 C -0.132 4.132 29 H 0.070 0.930 30 H 0.067 0.933 31 H 0.074 0.926 32 H 0.091 0.909 33 H 0.059 0.941 34 H 0.089 0.911 35 H 0.073 0.927 36 H 0.165 0.835 37 H 0.083 0.917 38 H 0.034 0.966 39 H 0.183 0.817 40 H 0.070 0.930 41 H 0.152 0.848 42 H 0.141 0.859 43 H 0.075 0.925 44 H 0.074 0.926 45 H 0.074 0.926 46 H 0.079 0.921 47 H 0.074 0.926 48 H 0.077 0.923 49 H 0.075 0.925 50 H 0.074 0.926 51 H 0.076 0.924 52 H 0.132 0.868 53 H 0.133 0.867 Dipole moment (debyes) X Y Z Total from point charges -17.828 11.660 -9.881 23.483 hybrid contribution -0.229 -2.256 -0.042 2.268 sum -18.057 9.404 -9.923 22.649 Atomic orbital electron populations 1.21875 0.95960 1.01152 1.01516 1.21006 0.94171 0.97200 0.97433 1.21892 1.02534 0.95079 1.00429 1.20503 0.90064 0.95930 0.91650 0.93964 1.60554 1.05561 1.43337 1.21428 1.22887 1.00895 0.88423 0.97028 1.19396 0.84821 0.97153 0.81060 0.91903 1.86590 1.75033 1.46715 1.28341 1.21597 0.92105 1.11524 1.35941 0.92842 1.25027 0.96772 0.96215 0.94636 1.70072 1.00642 1.36565 1.42075 0.86742 0.77459 0.83445 0.80727 1.20243 1.21202 1.20160 1.19944 0.95543 0.96347 0.96097 1.21206 0.99586 0.94309 0.95924 1.21187 0.97155 1.00564 0.95469 1.19722 0.95840 0.95915 0.96192 1.19476 0.95773 0.94750 0.92884 1.21878 0.95023 1.00916 1.01591 1.21958 0.99224 0.97350 1.01573 1.21886 1.02234 0.99027 0.96369 1.20814 1.00028 0.95829 0.97197 1.20976 0.97107 0.99825 0.95290 0.92978 0.93254 0.92635 0.90927 0.94075 0.91076 0.92734 0.83495 0.91737 0.96637 0.81662 0.92965 0.84813 0.85925 0.92549 0.92590 0.92557 0.92111 0.92618 0.92255 0.92545 0.92587 0.92388 0.86750 0.86682 Number of geometries 1 Number of calculations of the screened coulomb radii 1 The total number of SCF iterations 12. 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.72 8.91 37.16 0.33 -1.39 16 2 C -0.08 -1.09 2.82 -90.62 -0.26 -1.34 16 3 C -0.14 -1.92 7.84 37.16 0.29 -1.63 16 4 C 0.13 1.89 2.06 -68.86 -0.14 1.74 16 5 H 0.04 0.65 7.66 -51.93 -0.40 0.25 16 6 N -0.68 -11.82 5.16 -107.95 -0.56 -12.38 16 7 C 0.04 0.75 5.13 -3.82 -0.02 0.73 16 8 C 0.23 5.88 3.10 -27.89 -0.09 5.79 16 9 H 0.06 1.78 8.14 -51.93 -0.42 1.36 16 10 O -0.56 -14.34 12.07 -35.23 -0.43 -14.76 16 11 C -0.57 -15.75 8.76 -34.44 -0.30 -16.05 16 12 H 0.05 1.43 7.74 -52.49 -0.41 1.03 16 13 C 0.08 2.24 5.40 -17.82 -0.10 2.15 16 14 N -0.86 -25.50 12.84 55.43 0.71 -24.79 16 15 Si 0.89 21.27 29.54 -169.99 -5.02 16.25 16 16 H -0.28 -6.53 7.11 56.52 0.40 -6.13 16 17 H -0.29 -6.81 7.11 56.52 0.40 -6.41 16 18 H -0.28 -6.51 7.11 56.52 0.40 -6.11 16 19 C -0.08 -1.05 3.62 -104.63 -0.38 -1.43 16 20 C -0.09 -1.22 6.66 -39.58 -0.26 -1.48 16 21 C -0.13 -1.48 8.98 -39.58 -0.36 -1.84 16 22 C -0.08 -0.83 4.76 -104.63 -0.50 -1.33 16 23 C -0.03 -0.25 0.85 -155.66 -0.13 -0.38 16 24 C -0.14 -1.06 8.29 37.16 0.31 -0.75 16 25 C -0.14 -1.16 7.13 37.16 0.27 -0.89 16 26 C -0.14 -1.13 8.29 37.16 0.31 -0.82 16 27 C -0.12 -1.33 8.71 -39.58 -0.34 -1.68 16 28 C -0.11 -1.37 8.76 -39.58 -0.35 -1.72 16 29 H 0.05 0.56 8.14 -51.93 -0.42 0.14 16 30 H 0.05 0.55 8.14 -51.93 -0.42 0.13 16 31 H 0.05 0.64 6.76 -51.93 -0.35 0.29 16 32 H 0.07 1.10 7.35 -51.93 -0.38 0.72 16 33 H 0.04 0.51 8.14 -51.93 -0.42 0.09 16 34 H 0.07 1.08 7.57 -51.93 -0.39 0.69 16 35 H 0.05 0.71 5.64 -51.93 -0.29 0.42 16 36 H 0.35 6.30 7.39 -39.29 -0.29 6.01 16 37 H 0.06 1.32 7.36 -51.93 -0.38 0.94 16 38 H 0.02 0.31 8.05 -51.93 -0.42 -0.10 16 39 H 0.34 7.74 9.04 45.56 0.41 8.16 16 40 H 0.26 7.36 8.31 -40.82 -0.34 7.02 16 41 H 0.13 1.89 6.84 -52.49 -0.36 1.53 16 42 H 0.12 1.28 5.90 -52.49 -0.31 0.97 16 43 H 0.06 0.46 8.14 -51.92 -0.42 0.04 16 44 H 0.06 0.37 8.14 -51.93 -0.42 -0.05 16 45 H 0.06 0.43 7.39 -51.93 -0.38 0.04 16 46 H 0.06 0.51 7.51 -51.93 -0.39 0.12 16 47 H 0.05 0.39 8.14 -51.93 -0.42 -0.04 16 48 H 0.06 0.47 7.51 -51.93 -0.39 0.08 16 49 H 0.06 0.45 7.39 -51.93 -0.38 0.07 16 50 H 0.06 0.39 8.14 -51.93 -0.42 -0.03 16 51 H 0.06 0.52 8.14 -51.93 -0.42 0.10 16 52 H 0.11 1.07 7.10 -52.48 -0.37 0.70 16 53 H 0.12 1.27 8.04 -52.49 -0.42 0.85 16 LS Contribution 404.81 15.07 6.10 6.10 Total: -1.00 -29.30 404.81 -9.76 -39.07 By element: Atomic # 1 Polarization: 21.71 SS G_CDS: -8.85 Total: 12.85 kcal Atomic # 6 Polarization: -20.62 SS G_CDS: -1.72 Total: -22.34 kcal Atomic # 7 Polarization: -37.33 SS G_CDS: 0.16 Total: -37.17 kcal Atomic # 8 Polarization: -14.34 SS G_CDS: -0.43 Total: -14.76 kcal Atomic # 14 Polarization: 21.27 SS G_CDS: -5.02 Total: 16.25 kcal Total LS contribution 6.10 Total: 6.10 kcal Total: -29.30 -9.76 -39.07 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. ZINC000929859878.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.761 kcal (2) G-P(sol) polarization free energy of solvation -29.303 kcal (3) G-ENP(sol) elect.-nuc.-pol. free energy of system -45.064 kcal (4) G-CDS(sol) cavity-dispersion-solvent structure free energy -9.762 kcal (5) G-P-CDS(sol) = G-P(sol) + G-CDS(sol) = (2) + (4) -39.065 kcal (6) G-S(sol) free energy of system = (1) + (5) -54.826 kcal 1SCF WAS SPECIFIED, SO GEOMETRY OPTIMIZATION WAS NOT USED Total computer time = 0.73 seconds