Wall clock time and date at job start Sat Apr 11 2020 06:10:09 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.53004 * 1 3 3 C 1.52994 * 109.18840 * 2 1 4 4 C 1.47092 * 109.32802 * 119.57260 * 2 1 3 5 5 C 1.52572 * 114.32543 * 183.06355 * 4 2 1 6 6 H 1.09001 * 109.48846 * 34.97109 * 5 4 2 7 7 N 1.46496 * 109.49004 * 155.01722 * 5 4 2 8 8 C 1.34784 * 119.99719 * 205.27321 * 7 5 4 9 9 O 1.21284 * 119.99744 * 359.72581 * 8 7 5 10 10 C 1.50692 * 119.99878 * 179.72346 * 8 7 5 11 11 O 1.42902 * 109.47450 * 179.97438 * 10 8 7 12 12 C 1.42891 * 114.00312 * 180.02562 * 11 10 8 13 13 H 1.09004 * 110.46325 * 27.39210 * 12 11 10 14 14 C 1.54514 * 110.51298 * 150.00315 * 12 11 10 15 15 C 1.54480 * 104.06318 * 118.64351 * 14 12 11 16 16 O 1.44424 * 104.80239 * 23.98054 * 15 14 12 17 17 C 1.44425 * 105.27539 * 319.48495 * 16 15 14 18 18 C 1.50880 * 109.35559 * 274.98463 * 5 4 2 19 19 C 1.53974 * 110.10992 * 91.87502 * 18 5 4 20 20 N 1.47600 * 112.12821 * 319.45908 * 19 18 5 21 21 C 1.36931 * 118.64588 * 137.85882 * 20 19 18 22 22 H 1.08002 * 120.00241 * 205.57213 * 21 20 19 23 23 C 1.38818 * 120.00199 * 25.56483 * 21 20 19 24 24 N 1.31620 * 119.99429 * 17.90362 * 23 21 20 25 25 Si 1.86795 * 120.00259 * 197.90611 * 23 21 20 26 26 H 1.48505 * 109.47054 * 90.00171 * 25 23 21 27 27 H 1.48496 * 109.47150 * 209.99947 * 25 23 21 28 28 H 1.48502 * 109.47129 * 330.00178 * 25 23 21 29 29 C 1.44581 * 122.71174 * 317.83557 * 20 19 18 30 30 H 1.09003 * 109.46421 * 300.43470 * 1 2 3 31 31 H 1.08998 * 109.46983 * 60.42861 * 1 2 3 32 32 H 1.08990 * 109.47282 * 180.43661 * 1 2 3 33 33 H 1.08997 * 109.47311 * 179.56791 * 3 2 1 34 34 H 1.09002 * 109.47407 * 299.57024 * 3 2 1 35 35 H 1.09006 * 109.47072 * 59.56876 * 3 2 1 36 36 H 1.08993 * 108.51945 * 61.80366 * 4 2 1 37 37 H 1.09003 * 108.50699 * 304.31059 * 4 2 1 38 38 H 0.96999 * 120.00559 * 24.99506 * 7 5 4 39 39 H 1.09007 * 109.47522 * 300.00402 * 10 8 7 40 40 H 1.09008 * 109.47431 * 59.99936 * 10 8 7 41 41 H 1.08995 * 110.50847 * 237.29084 * 14 12 11 42 42 H 1.08999 * 110.50668 * 359.86022 * 14 12 11 43 43 H 1.09006 * 110.37083 * 142.81301 * 15 14 12 44 44 H 1.09002 * 110.37090 * 265.14946 * 15 14 12 45 45 H 1.08998 * 110.36961 * 281.68119 * 17 16 15 46 46 H 1.09003 * 110.36702 * 159.34219 * 17 16 15 47 47 H 1.09003 * 109.36440 * 212.08007 * 18 5 4 48 48 H 1.09000 * 109.36826 * 331.68866 * 18 5 4 49 49 H 1.09004 * 108.96302 * 198.76833 * 19 18 5 50 50 H 1.08999 * 108.94520 * 80.13815 * 19 18 5 51 51 H 0.97005 * 119.99544 * 179.97438 * 24 23 21 52 52 H 1.08999 * 109.45060 * 330.71236 * 29 20 19 53 53 H 1.08995 * 109.42982 * 210.77191 * 29 20 19 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.0329 1.4449 0.0000 4 6 2.0169 -0.6850 1.2072 5 6 3.5349 -0.6863 1.3600 6 1 3.9399 0.2498 0.9755 7 7 3.8835 -0.8246 2.7762 8 6 5.0644 -0.3580 3.2282 9 8 5.8406 0.1718 2.4614 10 6 5.4204 -0.4946 4.6861 11 8 6.7125 0.0707 4.9160 12 6 7.1480 -0.0079 6.2747 13 1 6.6877 -0.8608 6.7736 14 6 8.6884 -0.1064 6.3447 15 6 9.1099 1.1562 7.1286 16 8 7.9461 1.5141 7.9053 17 6 6.8212 1.3025 7.0246 18 6 4.1147 -1.8428 0.5836 19 6 4.4799 -1.3928 -0.8429 20 7 3.4624 -0.4923 -1.4194 21 6 3.8804 0.6049 -2.1241 22 1 3.2423 1.4729 -2.2005 23 6 5.1244 0.6016 -2.7401 24 7 5.7486 -0.5368 -2.9562 25 14 5.9124 2.2127 -3.2621 26 1 6.7411 2.7419 -2.1492 27 1 6.7706 1.9814 -4.4517 28 1 4.8519 3.1957 -3.6004 29 6 2.0425 -0.7211 -1.2712 30 1 -0.3632 0.5206 0.8861 31 1 -0.3633 0.5072 -0.8938 32 1 -0.3633 -1.0275 0.0078 33 1 3.1228 1.4504 0.0077 34 1 1.6734 1.9548 -0.8938 35 1 1.6607 1.9592 0.8861 36 1 1.6707 -1.7182 1.1824 37 1 1.5816 -0.1988 2.0803 38 1 3.2627 -1.2482 3.3894 39 1 4.6809 0.0299 5.2914 40 1 5.4318 -1.5496 4.9601 41 1 8.9895 -1.0067 6.8802 42 1 9.1175 -0.0975 5.3428 43 1 9.9504 0.9320 7.7856 44 1 9.3692 1.9610 6.4407 45 1 6.7295 2.1298 6.3208 46 1 5.9033 1.1906 7.6018 47 1 5.0101 -2.2047 1.0890 48 1 3.3812 -2.6474 0.5319 49 1 4.5746 -2.2731 -1.4787 50 1 5.4377 -0.8732 -0.8157 51 1 6.6177 -0.5392 -3.3870 52 1 1.8547 -1.7911 -1.1820 53 1 1.5194 -0.3323 -2.1448 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 ZINC001489342107.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:10:09 Heat of formation + Delta-G solvation = -100.303025 kcal Electronic energy + Delta-G solvation = -31732.283097 eV Core-core repulsion = 27604.971373 eV Total energy + Delta-G solvation = -4127.311724 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) -40.16612 -39.70093 -38.49273 -37.02563 -34.68589 -34.34964 -31.86803 -31.44434 -30.55216 -29.15197 -27.92063 -27.03091 -26.66974 -25.43064 -23.74882 -23.24629 -22.11305 -21.48479 -20.41418 -19.52988 -18.17854 -17.55447 -17.21058 -16.62513 -16.35510 -16.26613 -15.79412 -15.33890 -15.15460 -14.83221 -14.48332 -13.98962 -13.96620 -13.76472 -13.56804 -13.40364 -13.20266 -12.82943 -12.68275 -12.47481 -12.26403 -12.00719 -11.88329 -11.61330 -11.54507 -11.32114 -11.23879 -11.06563 -11.01932 -10.95898 -10.70984 -10.55168 -10.24454 -10.22746 -10.10049 -9.85532 -9.82240 -9.67493 -9.56090 -9.49403 -9.38341 -8.70628 -8.26811 -6.58994 -5.69085 -3.10592 2.78482 2.84322 3.47255 3.53588 3.55866 3.62631 3.91531 4.21586 4.44386 4.58071 4.59884 4.77186 4.90336 4.92380 4.98018 5.03466 5.10835 5.23943 5.33247 5.40728 5.42473 5.46652 5.53542 5.57195 5.58496 5.64899 5.72200 5.96562 6.00989 6.14098 6.18720 6.25054 6.27064 6.33367 6.36954 6.40641 6.50882 6.61738 6.66227 6.78927 6.81397 6.84805 6.89007 6.98149 7.12104 7.35177 7.49563 7.75637 7.79486 7.94450 8.50258 8.56186 9.32809 10.00426 10.44347 11.38623 Molecular weight = 340.22amu Principal moments of inertia in cm(-1) A = 0.011415 B = 0.003084 C = 0.002681 Principal moments of inertia in units of 10**(-40)*gram-cm**2 A = 2452.286841 B = 9077.615219 C =10443.145475 Net atomic charges, atomic populations, and dipole contributions calculated with CM2 Atom NO. Type Charge No. of electrons 1 C -0.133 4.133 2 C -0.093 4.093 3 C -0.131 4.131 4 C -0.096 4.096 5 C 0.151 3.849 6 H 0.110 0.890 7 N -0.727 5.727 8 C 0.505 3.495 9 O -0.504 6.504 10 C 0.042 3.958 11 O -0.346 6.346 12 C 0.054 3.946 13 H 0.079 0.921 14 C -0.156 4.156 15 C 0.045 3.955 16 O -0.374 6.374 17 C 0.023 3.977 18 C -0.147 4.147 19 C 0.182 3.818 20 N -0.587 5.587 21 C -0.258 4.258 22 H 0.070 0.930 23 C 0.007 3.993 24 N -0.901 5.901 25 Si 0.894 3.106 26 H -0.286 1.286 27 H -0.293 1.293 28 H -0.280 1.280 29 C 0.169 3.831 30 H 0.043 0.957 31 H 0.053 0.947 32 H 0.047 0.953 33 H 0.082 0.918 34 H 0.047 0.953 35 H 0.026 0.974 36 H 0.053 0.947 37 H 0.048 0.952 38 H 0.390 0.610 39 H 0.077 0.923 40 H 0.078 0.922 41 H 0.084 0.916 42 H 0.096 0.904 43 H 0.094 0.906 44 H 0.065 0.935 45 H 0.074 0.926 46 H 0.094 0.906 47 H 0.054 0.946 48 H 0.045 0.955 49 H 0.025 0.975 50 H 0.078 0.922 51 H 0.252 0.748 52 H 0.031 0.969 53 H 0.046 0.954 For the charges calculated by CM2: Dipole moment (debyes) X Y Z Total from point charges -3.591 -2.755 20.075 20.579 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.191 4.191 2 C -0.097 4.097 3 C -0.188 4.188 4 C -0.134 4.134 5 C 0.048 3.952 6 H 0.128 0.872 7 N -0.380 5.380 8 C 0.289 3.711 9 O -0.377 6.377 10 C -0.038 4.038 11 O -0.264 6.264 12 C -0.005 4.005 13 H 0.097 0.903 14 C -0.194 4.194 15 C -0.031 4.031 16 O -0.292 6.292 17 C -0.055 4.055 18 C -0.187 4.187 19 C 0.056 3.944 20 N -0.307 5.307 21 C -0.388 4.388 22 H 0.088 0.912 23 C -0.192 4.192 24 N -0.545 5.545 25 Si 0.724 3.276 26 H -0.213 1.213 27 H -0.219 1.219 28 H -0.205 1.205 29 C 0.040 3.960 30 H 0.062 0.938 31 H 0.072 0.928 32 H 0.066 0.934 33 H 0.101 0.899 34 H 0.066 0.934 35 H 0.045 0.955 36 H 0.071 0.929 37 H 0.067 0.933 38 H 0.223 0.777 39 H 0.095 0.905 40 H 0.096 0.904 41 H 0.103 0.897 42 H 0.114 0.886 43 H 0.112 0.888 44 H 0.083 0.917 45 H 0.092 0.908 46 H 0.112 0.888 47 H 0.073 0.927 48 H 0.064 0.936 49 H 0.043 0.957 50 H 0.096 0.904 51 H 0.062 0.938 52 H 0.049 0.951 53 H 0.065 0.935 Dipole moment (debyes) X Y Z Total from point charges -2.713 -2.253 19.471 19.788 hybrid contribution -0.052 0.280 0.109 0.304 sum -2.765 -1.973 19.580 19.872 Atomic orbital electron populations 1.21803 0.95672 1.00679 1.00943 1.20464 0.93882 0.97346 0.97973 1.21930 1.03564 0.94290 0.99027 1.20593 0.94140 1.00362 0.98277 1.20650 0.94764 1.02104 0.77680 0.87195 1.46212 1.18721 1.61512 1.11542 1.20085 0.82376 0.79456 0.89192 1.90825 1.48581 1.41109 1.57220 1.22348 0.86850 0.99147 0.95434 1.88133 1.32336 1.80663 1.25245 1.23565 0.95714 0.97891 0.83338 0.90341 1.23032 0.94048 0.99340 1.02976 1.23433 0.89426 0.95641 0.94630 1.89129 1.15545 1.88401 1.36165 1.23971 0.90609 0.97356 0.93519 1.22448 0.99614 0.96772 0.99873 1.20946 0.93216 0.93591 0.86625 1.43645 0.98967 1.28423 1.59694 1.19259 1.01414 0.92587 1.25503 0.91227 1.24980 0.99695 0.96262 0.98246 1.70076 1.15742 1.23644 1.45053 0.86415 0.79501 0.83922 0.77718 1.21324 1.21948 1.20520 1.20715 0.86135 0.95681 0.93480 0.93787 0.92806 0.93430 0.89946 0.93402 0.95510 0.92853 0.93301 0.77687 0.90542 0.90390 0.89711 0.88564 0.88767 0.91659 0.90769 0.88782 0.92698 0.93643 0.95685 0.90410 0.93797 0.95086 0.93519 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.13 -1.80 8.37 37.16 0.31 -1.49 16 2 C -0.09 -1.53 0.72 -156.36 -0.11 -1.64 16 3 C -0.13 -2.38 6.77 37.15 0.25 -2.13 16 4 C -0.10 -1.43 3.40 -29.79 -0.10 -1.53 16 5 C 0.15 2.60 1.98 -69.22 -0.14 2.47 16 6 H 0.11 2.41 4.73 -51.93 -0.25 2.16 16 7 N -0.73 -9.99 4.93 -55.06 -0.27 -10.26 16 8 C 0.51 7.92 7.74 -10.99 -0.08 7.84 16 9 O -0.50 -10.98 16.16 -14.35 -0.23 -11.21 16 10 C 0.04 0.46 6.25 36.00 0.23 0.68 16 11 O -0.35 -4.54 9.51 -55.14 -0.52 -5.07 16 12 C 0.05 0.51 3.61 -25.16 -0.09 0.42 16 13 H 0.08 0.58 7.82 -51.93 -0.41 0.17 16 14 C -0.16 -1.39 6.90 -25.16 -0.17 -1.56 16 15 C 0.05 0.40 7.63 37.93 0.29 0.69 16 16 O -0.37 -4.02 11.65 -35.23 -0.41 -4.43 16 17 C 0.02 0.22 7.38 37.93 0.28 0.50 16 18 C -0.15 -2.81 5.39 -27.29 -0.15 -2.95 16 19 C 0.18 4.37 4.68 -2.93 -0.01 4.36 16 20 N -0.59 -14.04 1.86 -172.07 -0.32 -14.36 16 21 C -0.26 -6.89 7.64 -15.06 -0.12 -7.00 16 22 H 0.07 1.83 6.44 -52.49 -0.34 1.49 16 23 C 0.01 0.20 4.81 -17.43 -0.08 0.11 16 24 N -0.90 -26.51 11.10 55.49 0.62 -25.89 16 25 Si 0.89 22.47 29.59 -169.99 -5.03 17.44 16 26 H -0.29 -7.38 7.11 56.52 0.40 -6.98 16 27 H -0.29 -7.42 7.11 56.52 0.40 -7.02 16 28 H -0.28 -6.82 7.11 56.52 0.40 -6.42 16 29 C 0.17 3.36 4.70 -4.08 -0.02 3.34 16 30 H 0.04 0.54 8.14 -51.93 -0.42 0.12 16 31 H 0.05 0.73 8.14 -51.93 -0.42 0.31 16 32 H 0.05 0.59 8.14 -51.93 -0.42 0.16 16 33 H 0.08 1.78 4.48 -51.93 -0.23 1.55 16 34 H 0.05 0.88 7.50 -51.93 -0.39 0.49 16 35 H 0.03 0.43 8.14 -51.93 -0.42 0.00 16 36 H 0.05 0.74 7.36 -51.93 -0.38 0.36 16 37 H 0.05 0.64 8.01 -51.93 -0.42 0.23 16 38 H 0.39 3.77 8.60 -40.82 -0.35 3.42 16 39 H 0.08 0.61 8.14 -51.93 -0.42 0.18 16 40 H 0.08 0.63 8.11 -51.92 -0.42 0.21 16 41 H 0.08 0.60 8.14 -51.93 -0.42 0.17 16 42 H 0.10 0.97 8.00 -51.93 -0.42 0.55 16 43 H 0.09 0.64 8.14 -51.93 -0.42 0.21 16 44 H 0.07 0.58 8.14 -51.93 -0.42 0.16 16 45 H 0.07 0.79 8.08 -51.93 -0.42 0.37 16 46 H 0.09 0.70 8.14 -51.93 -0.42 0.28 16 47 H 0.05 1.06 7.86 -51.93 -0.41 0.66 16 48 H 0.04 0.79 7.32 -51.93 -0.38 0.41 16 49 H 0.03 0.65 8.05 -51.93 -0.42 0.23 16 50 H 0.08 2.18 6.33 -51.93 -0.33 1.86 16 51 H 0.25 7.27 8.31 -40.82 -0.34 6.93 16 52 H 0.03 0.59 8.07 -51.93 -0.42 0.17 16 53 H 0.05 0.96 7.95 -51.93 -0.41 0.55 16 LS Contribution 400.35 15.07 6.03 6.03 Total: -1.00 -33.47 400.35 -9.18 -42.66 By element: Atomic # 1 Polarization: 12.32 SS G_CDS: -9.32 Total: 2.99 kcal Atomic # 6 Polarization: 1.82 SS G_CDS: 0.28 Total: 2.10 kcal Atomic # 7 Polarization: -50.54 SS G_CDS: 0.02 Total: -50.51 kcal Atomic # 8 Polarization: -19.54 SS G_CDS: -1.17 Total: -20.71 kcal Atomic # 14 Polarization: 22.47 SS G_CDS: -5.03 Total: 17.44 kcal Total LS contribution 6.03 Total: 6.03 kcal Total: -33.47 -9.18 -42.66 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. ZINC001489342107.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 -57.646 kcal (2) G-P(sol) polarization free energy of solvation -33.474 kcal (3) G-ENP(sol) elect.-nuc.-pol. free energy of system -91.120 kcal (4) G-CDS(sol) cavity-dispersion-solvent structure free energy -9.183 kcal (5) G-P-CDS(sol) = G-P(sol) + G-CDS(sol) = (2) + (4) -42.657 kcal (6) G-S(sol) free energy of system = (1) + (5) -100.303 kcal 1SCF WAS SPECIFIED, SO GEOMETRY OPTIMIZATION WAS NOT USED Total computer time = 1.08 seconds