Wall clock time and date at job start Sat Apr 11 2020 06:11:59 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 N 1.46498 * 1 3 3 C 1.34776 * 119.99934 * 2 1 4 4 O 1.21591 * 120.00146 * 0.02562 * 3 2 1 5 5 N 1.34779 * 119.99963 * 180.02562 * 3 2 1 6 6 C 1.46501 * 119.99753 * 180.02562 * 5 3 2 7 7 C 1.53001 * 109.47148 * 179.97438 * 6 5 3 8 8 C 1.50699 * 109.47598 * 179.97438 * 7 6 5 9 9 O 1.21286 * 119.99302 * 359.97438 * 8 7 6 10 10 N 1.34770 * 120.00187 * 179.97438 * 8 7 6 11 11 C 1.46505 * 120.00002 * 180.02562 * 10 8 7 12 12 C 1.53001 * 109.47161 * 180.02562 * 11 10 8 13 13 H 1.08999 * 109.69952 * 299.31236 * 12 11 10 14 14 C 1.52757 * 109.50726 * 179.19298 * 12 11 10 15 15 C 1.53042 * 109.38216 * 185.70762 * 14 12 11 16 16 C 1.53038 * 109.53901 * 61.23851 * 15 14 12 17 17 C 1.53001 * 109.50154 * 178.68295 * 16 15 14 18 18 C 1.53004 * 109.52444 * 58.60569 * 16 15 14 19 19 C 1.53187 * 109.31521 * 298.63658 * 16 15 14 20 20 N 1.46923 * 109.58514 * 59.66713 * 12 11 10 21 21 C 1.36945 * 120.63417 * 353.41495 * 20 12 11 22 22 H 1.08004 * 119.99143 * 30.96028 * 21 20 12 23 23 C 1.38804 * 120.00073 * 210.95649 * 21 20 12 24 24 N 1.31618 * 120.00327 * 21.01119 * 23 21 20 25 25 Si 1.86809 * 119.99939 * 201.00638 * 23 21 20 26 26 H 1.48503 * 109.46840 * 270.00287 * 25 23 21 27 27 H 1.48497 * 109.46819 * 30.00103 * 25 23 21 28 28 H 1.48497 * 109.47085 * 150.00458 * 25 23 21 29 29 H 1.09002 * 109.46814 * 179.97438 * 1 2 3 30 30 H 1.08994 * 109.47286 * 299.99683 * 1 2 3 31 31 H 1.09002 * 109.47188 * 59.99974 * 1 2 3 32 32 H 0.96997 * 120.00501 * 180.02562 * 2 1 3 33 33 H 0.96993 * 120.00050 * 359.97438 * 5 3 2 34 34 H 1.08993 * 109.47870 * 59.98886 * 6 5 3 35 35 H 1.09000 * 109.46719 * 299.98995 * 6 5 3 36 36 H 1.09009 * 109.47121 * 59.99225 * 7 6 5 37 37 H 1.09001 * 109.47193 * 299.99675 * 7 6 5 38 38 H 0.97002 * 120.00221 * 359.97438 * 10 8 7 39 39 H 1.08996 * 109.46932 * 60.00023 * 11 10 8 40 40 H 1.09000 * 109.46980 * 299.99990 * 11 10 8 41 41 H 1.08996 * 109.56522 * 65.67052 * 14 12 11 42 42 H 1.09006 * 109.36085 * 305.61625 * 14 12 11 43 43 H 1.08998 * 109.46104 * 301.21726 * 15 14 12 44 44 H 1.09000 * 109.45823 * 181.25908 * 15 14 12 45 45 H 1.08998 * 109.47757 * 179.83988 * 17 16 15 46 46 H 1.08999 * 109.47023 * 299.84404 * 17 16 15 47 47 H 1.08998 * 109.46757 * 59.83756 * 17 16 15 48 48 H 1.09002 * 109.46900 * 300.14715 * 18 16 15 49 49 H 1.08996 * 109.47092 * 60.14501 * 18 16 15 50 50 H 1.08998 * 109.47593 * 180.14905 * 18 16 15 51 51 H 1.08996 * 109.58681 * 294.84959 * 19 16 15 52 52 H 1.08999 * 109.70973 * 174.48860 * 19 16 15 53 53 H 0.97002 * 120.00122 * 4.98040 * 24 23 21 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 7 1.4650 0.0000 0.0000 3 6 2.1388 1.1672 0.0000 4 8 1.5309 2.2202 0.0005 5 7 3.4866 1.1672 -0.0005 6 6 4.2191 2.4360 0.0000 7 6 5.7233 2.1565 -0.0001 8 6 6.4769 3.4616 -0.0001 9 8 5.8703 4.5119 -0.0005 10 7 7.8246 3.4616 0.0003 11 6 8.5571 4.7303 0.0008 12 6 10.0613 4.4508 0.0005 13 1 10.3309 3.8893 0.8951 14 6 10.8259 5.7731 -0.0192 15 6 12.3234 5.4918 -0.1630 16 6 12.5827 4.7622 -1.4830 17 6 14.0835 4.5105 -1.6413 18 6 12.0831 5.6193 -2.6478 19 6 11.8377 3.4236 -1.4769 20 7 10.4177 3.6740 -1.1946 21 6 9.4452 3.1973 -2.0326 22 1 8.5208 3.7429 -2.1525 23 6 9.6475 2.0130 -2.7277 24 7 10.5414 1.1451 -2.3033 25 14 8.6508 1.6508 -4.2657 26 1 7.4081 0.9290 -3.8916 27 1 8.2986 2.9266 -4.9391 28 1 9.4535 0.8089 -5.1888 29 1 -0.3633 -1.0277 -0.0005 30 1 -0.3633 0.5138 0.8900 31 1 -0.3634 0.5138 -0.8900 32 1 1.9500 -0.8400 -0.0004 33 1 3.9716 0.3272 -0.0005 34 1 3.9560 3.0071 0.8903 35 1 3.9555 3.0080 -0.8896 36 1 5.9867 1.5852 -0.8903 37 1 5.9867 1.5846 0.8897 38 1 8.3096 2.6215 0.0010 39 1 8.2935 5.3016 0.8908 40 1 8.2940 5.3021 -0.8891 41 1 10.6448 6.3133 0.9100 42 1 10.4861 6.3744 -0.8625 43 1 12.6574 4.8698 0.6674 44 1 12.8722 6.4336 -0.1547 45 1 14.2675 3.9883 -2.5802 46 1 14.6132 5.4632 -1.6456 47 1 14.4392 3.9009 -0.8107 48 1 11.0141 5.7993 -2.5340 49 1 12.6134 6.5715 -2.6518 50 1 12.2657 5.0978 -3.5873 51 1 12.2551 2.7767 -0.7053 52 1 11.9375 2.9429 -2.4501 53 1 11.0071 1.2986 -1.4664 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 ZINC001099055373.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:11:59 Heat of formation + Delta-G solvation = -59.185777 kcal Electronic energy + Delta-G solvation = -31221.293854 eV Core-core repulsion = 27128.057859 eV Total energy + Delta-G solvation = -4093.235995 eV No. of doubly occupied orbitals = 66 Molecular weight (most abundant/longest-lived isotopes) = 340.226 amu Computer time = 1.02 seconds Orbital eigenvalues (eV) -40.01217 -39.42082 -38.57916 -35.51673 -35.00476 -34.53132 -33.50772 -32.31507 -30.79990 -29.83759 -27.67302 -26.74337 -26.58157 -25.51489 -24.03886 -22.59471 -22.14367 -21.02083 -20.28289 -19.76835 -18.64597 -17.95861 -17.27292 -16.66639 -16.02330 -15.89063 -15.52966 -15.28829 -14.87623 -14.42945 -14.42166 -14.18489 -13.95655 -13.67468 -13.62581 -13.43393 -13.24396 -12.99846 -12.74089 -12.45935 -12.40686 -12.13590 -11.93732 -11.72296 -11.58966 -11.44473 -11.29859 -11.22480 -10.82950 -10.80131 -10.72316 -10.32373 -10.30719 -10.16259 -10.01402 -9.92324 -9.65684 -9.39956 -9.39379 -9.12390 -8.89060 -8.55907 -8.24522 -6.91956 -6.10993 -3.36107 2.61262 3.15801 3.17905 3.38219 3.46822 3.60455 3.79545 4.26773 4.46188 4.88249 5.00198 5.01353 5.21977 5.22993 5.30435 5.34364 5.34961 5.36948 5.50306 5.59046 5.62887 5.66622 5.77364 5.89742 5.99500 6.04638 6.06440 6.09307 6.11900 6.19910 6.24567 6.28361 6.33889 6.35741 6.36842 6.41902 6.45233 6.53109 6.71648 6.79674 6.80183 6.84193 6.91090 7.11928 7.14177 7.19207 7.91578 7.99303 8.20166 8.50061 8.54643 8.82330 9.00151 9.74178 10.20386 11.13608 Molecular weight = 340.23amu Principal moments of inertia in cm(-1) A = 0.016154 B = 0.002727 C = 0.002456 Principal moments of inertia in units of 10**(-40)*gram-cm**2 A = 1732.907998 B =10266.020964 C =11398.435610 Net atomic charges, atomic populations, and dipole contributions calculated with CM2 Atom NO. Type Charge No. of electrons 1 C 0.108 3.892 2 N -0.754 5.754 3 C 0.701 3.299 4 O -0.589 6.589 5 N -0.746 5.746 6 C 0.156 3.844 7 C -0.154 4.154 8 C 0.512 3.488 9 O -0.562 6.562 10 N -0.690 5.690 11 C 0.138 3.862 12 C 0.163 3.837 13 H 0.028 0.972 14 C -0.131 4.131 15 C -0.095 4.095 16 C -0.089 4.089 17 C -0.135 4.135 18 C -0.125 4.125 19 C 0.166 3.834 20 N -0.592 5.592 21 C -0.341 4.341 22 H 0.062 0.938 23 C -0.018 4.018 24 N -0.918 5.918 25 Si 1.057 2.943 26 H -0.295 1.295 27 H -0.293 1.293 28 H -0.287 1.287 29 H 0.073 0.927 30 H 0.057 0.943 31 H 0.057 0.943 32 H 0.397 0.603 33 H 0.396 0.604 34 H 0.073 0.927 35 H 0.075 0.925 36 H 0.102 0.898 37 H 0.094 0.906 38 H 0.413 0.587 39 H 0.054 0.946 40 H 0.065 0.935 41 H 0.052 0.948 42 H 0.066 0.934 43 H 0.052 0.948 44 H 0.047 0.953 45 H 0.051 0.949 46 H 0.044 0.956 47 H 0.048 0.952 48 H 0.066 0.934 49 H 0.033 0.967 50 H 0.045 0.955 51 H 0.015 0.985 52 H 0.078 0.922 53 H 0.270 0.730 For the charges calculated by CM2: Dipole moment (debyes) X Y Z Total from point charges -5.110 -2.532 5.791 8.128 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.037 4.037 2 N -0.410 5.410 3 C 0.401 3.599 4 O -0.466 6.466 5 N -0.401 5.401 6 C 0.032 3.968 7 C -0.194 4.194 8 C 0.297 3.703 9 O -0.442 6.442 10 N -0.340 5.340 11 C 0.013 3.987 12 C 0.056 3.944 13 H 0.046 0.954 14 C -0.171 4.171 15 C -0.133 4.133 16 C -0.091 4.091 17 C -0.193 4.193 18 C -0.182 4.182 19 C 0.040 3.960 20 N -0.317 5.317 21 C -0.471 4.471 22 H 0.079 0.921 23 C -0.221 4.221 24 N -0.558 5.558 25 Si 0.888 3.112 26 H -0.223 1.223 27 H -0.219 1.219 28 H -0.213 1.213 29 H 0.092 0.908 30 H 0.075 0.925 31 H 0.076 0.924 32 H 0.230 0.770 33 H 0.229 0.771 34 H 0.091 0.909 35 H 0.093 0.907 36 H 0.121 0.879 37 H 0.112 0.888 38 H 0.251 0.749 39 H 0.072 0.928 40 H 0.084 0.916 41 H 0.070 0.930 42 H 0.085 0.915 43 H 0.070 0.930 44 H 0.066 0.934 45 H 0.071 0.929 46 H 0.063 0.937 47 H 0.067 0.933 48 H 0.085 0.915 49 H 0.052 0.948 50 H 0.065 0.935 51 H 0.033 0.967 52 H 0.096 0.904 53 H 0.082 0.918 Dipole moment (debyes) X Y Z Total from point charges -6.337 -1.881 4.528 8.012 hybrid contribution 1.049 -0.176 1.120 1.544 sum -5.288 -2.056 5.648 8.005 Atomic orbital electron populations 1.21365 0.78305 1.01763 1.02256 1.44921 1.07711 1.07176 1.81211 1.15946 0.81796 0.84235 0.77958 1.90947 1.67870 1.30299 1.57491 1.45146 1.04686 1.10097 1.80202 1.20883 0.88402 0.83764 1.03792 1.21771 0.94802 0.97671 1.05170 1.20110 0.82910 0.90347 0.76905 1.90681 1.68487 1.33452 1.51609 1.45955 1.05263 1.10787 1.72012 1.21230 0.91075 0.85750 1.00597 1.20619 0.91233 0.91886 0.90679 0.95404 1.21935 0.95897 0.98415 1.00828 1.21405 0.95737 0.98976 0.97190 1.21120 0.94471 0.98704 0.94790 1.21786 0.95876 1.00403 1.01274 1.21701 1.01824 0.97408 0.97272 1.21002 0.83347 0.91568 1.00087 1.44111 1.01455 1.62569 1.23552 1.19999 0.99327 1.04992 1.22751 0.92063 1.25705 1.01065 0.91507 1.03803 1.70768 1.27610 1.34332 1.23134 0.83624 0.75354 0.75475 0.76761 1.22258 1.21937 1.21300 0.90800 0.92470 0.92434 0.77009 0.77055 0.90850 0.90666 0.87946 0.88751 0.74900 0.92824 0.91639 0.92978 0.91528 0.92968 0.93385 0.92939 0.93673 0.93325 0.91469 0.94813 0.93538 0.96708 0.90392 0.91843 Number of geometries 1 Number of calculations of the screened coulomb radii 1 The total number of SCF iterations 17. 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.11 0.86 10.70 59.85 0.64 1.50 16 2 N -0.75 -7.04 5.72 -72.07 -0.41 -7.45 16 3 C 0.70 9.60 8.57 -86.92 -0.74 8.86 16 4 O -0.59 -10.85 16.34 5.29 0.09 -10.76 16 5 N -0.75 -8.84 5.56 -65.23 -0.36 -9.20 16 6 C 0.16 2.25 4.96 -4.04 -0.02 2.23 16 7 C -0.15 -2.09 5.81 -27.88 -0.16 -2.25 16 8 C 0.51 9.06 7.83 -10.99 -0.09 8.98 16 9 O -0.56 -11.77 15.98 5.55 0.09 -11.68 16 10 N -0.69 -12.36 3.81 -60.29 -0.23 -12.59 16 11 C 0.14 2.51 4.43 -4.04 -0.02 2.49 16 12 C 0.16 2.95 3.00 -67.82 -0.20 2.75 16 13 H 0.03 0.52 8.14 -51.93 -0.42 0.10 16 14 C -0.13 -2.10 5.00 -26.83 -0.13 -2.24 16 15 C -0.10 -1.48 4.71 -26.69 -0.13 -1.61 16 16 C -0.09 -1.52 0.77 -154.39 -0.12 -1.64 16 17 C -0.14 -2.03 8.48 37.16 0.32 -1.72 16 18 C -0.12 -2.16 7.99 37.16 0.30 -1.87 16 19 C 0.17 3.42 3.33 -3.71 -0.01 3.41 16 20 N -0.59 -12.45 2.26 -165.03 -0.37 -12.82 16 21 C -0.34 -8.00 6.86 -15.06 -0.10 -8.11 16 22 H 0.06 1.52 5.89 -52.48 -0.31 1.21 16 23 C -0.02 -0.45 4.77 -17.43 -0.08 -0.53 16 24 N -0.92 -24.27 12.50 55.49 0.69 -23.57 16 25 Si 1.06 23.55 29.97 -169.99 -5.09 18.46 16 26 H -0.29 -6.32 7.11 56.52 0.40 -5.92 16 27 H -0.29 -6.53 7.11 56.52 0.40 -6.13 16 28 H -0.29 -6.59 7.11 56.52 0.40 -6.19 16 29 H 0.07 0.33 8.14 -51.93 -0.42 -0.09 16 30 H 0.06 0.50 8.14 -51.93 -0.42 0.08 16 31 H 0.06 0.52 8.14 -51.93 -0.42 0.10 16 32 H 0.40 2.39 8.84 -40.82 -0.36 2.03 16 33 H 0.40 3.38 8.49 -40.82 -0.35 3.03 16 34 H 0.07 1.16 8.10 -51.93 -0.42 0.74 16 35 H 0.07 1.27 8.10 -51.93 -0.42 0.85 16 36 H 0.10 1.38 8.14 -51.92 -0.42 0.95 16 37 H 0.09 1.04 8.14 -51.93 -0.42 0.61 16 38 H 0.41 7.63 7.35 -40.82 -0.30 7.33 16 39 H 0.05 0.92 8.14 -51.93 -0.42 0.50 16 40 H 0.07 1.34 7.45 -51.93 -0.39 0.95 16 41 H 0.05 0.74 8.14 -51.93 -0.42 0.31 16 42 H 0.07 1.13 6.34 -51.93 -0.33 0.80 16 43 H 0.05 0.80 8.14 -51.93 -0.42 0.38 16 44 H 0.05 0.68 8.14 -51.93 -0.42 0.25 16 45 H 0.05 0.80 8.14 -51.93 -0.42 0.38 16 46 H 0.04 0.60 8.14 -51.93 -0.42 0.18 16 47 H 0.05 0.70 8.14 -51.93 -0.42 0.28 16 48 H 0.07 1.29 5.97 -51.93 -0.31 0.98 16 49 H 0.03 0.51 8.14 -51.93 -0.42 0.09 16 50 H 0.05 0.84 8.14 -51.93 -0.42 0.41 16 51 H 0.01 0.32 7.78 -51.93 -0.40 -0.08 16 52 H 0.08 1.87 6.43 -51.93 -0.33 1.53 16 53 H 0.27 7.16 6.45 -40.82 -0.26 6.89 16 LS Contribution 410.01 15.07 6.18 6.18 Total: -1.00 -31.32 410.01 -9.31 -40.62 By element: Atomic # 1 Polarization: 21.88 SS G_CDS: -9.32 Total: 12.55 kcal Atomic # 6 Polarization: 10.83 SS G_CDS: -0.56 Total: 10.27 kcal Atomic # 7 Polarization: -64.96 SS G_CDS: -0.68 Total: -65.64 kcal Atomic # 8 Polarization: -22.62 SS G_CDS: 0.18 Total: -22.44 kcal Atomic # 14 Polarization: 23.55 SS G_CDS: -5.09 Total: 18.46 kcal Total LS contribution 6.18 Total: 6.18 kcal Total: -31.32 -9.31 -40.62 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. ZINC001099055373.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 -18.561 kcal (2) G-P(sol) polarization free energy of solvation -31.319 kcal (3) G-ENP(sol) elect.-nuc.-pol. free energy of system -49.880 kcal (4) G-CDS(sol) cavity-dispersion-solvent structure free energy -9.306 kcal (5) G-P-CDS(sol) = G-P(sol) + G-CDS(sol) = (2) + (4) -40.625 kcal (6) G-S(sol) free energy of system = (1) + (5) -59.186 kcal 1SCF WAS SPECIFIED, SO GEOMETRY OPTIMIZATION WAS NOT USED Total computer time = 1.02 seconds