Wall clock time and date at job start Sat Apr 11 2020 06:11:01 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.53005 * 1 3 3 N 1.46500 * 109.46911 * 2 1 4 4 C 1.34779 * 120.00223 * 80.46475 * 3 2 1 5 5 O 1.21509 * 120.00010 * 179.97438 * 4 3 2 6 6 O 1.34634 * 119.99663 * 359.72949 * 4 3 2 7 7 C 1.45201 * 117.00044 * 180.27494 * 6 4 3 8 8 C 1.52997 * 109.47315 * 59.99691 * 7 6 4 9 9 C 1.53002 * 109.47014 * 179.97438 * 7 6 4 10 10 C 1.52998 * 109.46750 * 299.99458 * 7 6 4 11 11 C 1.46503 * 120.00066 * 260.46841 * 3 2 1 12 12 H 1.08998 * 109.47271 * 173.83560 * 11 3 2 13 13 C 1.53001 * 109.46953 * 293.83543 * 11 3 2 14 14 C 1.50703 * 109.47152 * 53.83878 * 11 3 2 15 15 O 1.21278 * 119.99731 * 315.75772 * 14 11 3 16 16 N 1.34769 * 119.99968 * 135.76600 * 14 11 3 17 17 C 1.46500 * 119.99784 * 185.26671 * 16 14 11 18 18 C 1.53009 * 109.47273 * 184.90793 * 17 16 14 19 19 C 1.50703 * 115.54508 * 65.35482 * 18 17 16 20 20 H 1.08001 * 119.99854 * 268.61723 * 19 18 17 21 21 C 1.37870 * 119.99716 * 88.62236 * 19 18 17 22 22 N 1.31513 * 120.00301 * 3.11746 * 21 19 18 23 23 Si 1.86803 * 119.99883 * 183.11836 * 21 19 18 24 24 H 1.48497 * 109.47263 * 0.02562 * 23 21 19 25 25 H 1.48499 * 109.47026 * 120.00035 * 23 21 19 26 26 H 1.48504 * 109.47090 * 239.99850 * 23 21 19 27 27 C 1.52998 * 117.49613 * 279.69076 * 18 17 16 28 28 C 1.52993 * 117.50183 * 211.02974 * 18 17 16 29 29 H 1.09005 * 109.46647 * 305.48948 * 1 2 3 30 30 H 1.08996 * 109.47446 * 65.49428 * 1 2 3 31 31 H 1.09003 * 109.46800 * 185.49871 * 1 2 3 32 32 H 1.09003 * 109.46983 * 240.00409 * 2 1 3 33 33 H 1.08998 * 109.47603 * 119.99795 * 2 1 3 34 34 H 1.08994 * 109.47132 * 59.99819 * 8 7 6 35 35 H 1.09001 * 109.47052 * 180.02562 * 8 7 6 36 36 H 1.09003 * 109.46646 * 300.00255 * 8 7 6 37 37 H 1.09007 * 109.47257 * 60.00044 * 9 7 6 38 38 H 1.09001 * 109.46937 * 179.97438 * 9 7 6 39 39 H 1.08993 * 109.47319 * 299.99488 * 9 7 6 40 40 H 1.08997 * 109.47306 * 60.00185 * 10 7 6 41 41 H 1.09001 * 109.47090 * 180.02562 * 10 7 6 42 42 H 1.09007 * 109.47534 * 299.99964 * 10 7 6 43 43 H 1.09000 * 109.46691 * 177.56281 * 13 11 3 44 44 H 1.09000 * 109.47263 * 297.56554 * 13 11 3 45 45 H 1.09002 * 109.46808 * 57.57047 * 13 11 3 46 46 H 0.97003 * 120.00342 * 5.26719 * 16 14 11 47 47 H 1.08998 * 109.47194 * 304.90967 * 17 16 14 48 48 H 1.09001 * 109.47402 * 64.90973 * 17 16 14 49 49 H 0.96999 * 120.00127 * 185.00703 * 22 21 19 50 50 H 1.08994 * 117.50268 * 359.97438 * 27 18 17 51 51 H 1.09003 * 117.49756 * 145.02393 * 27 18 17 52 52 H 1.09004 * 117.49871 * 214.98265 * 28 18 17 53 53 H 1.09011 * 117.49831 * 359.97136 * 28 18 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 7 2.0183 1.3812 0.0000 4 6 2.0607 2.0811 1.1511 5 8 2.4661 3.2265 1.1511 6 8 1.6491 1.5126 2.3000 7 6 1.7170 2.3273 3.5000 8 6 0.8450 3.5717 3.3216 9 6 1.2109 1.5176 4.6955 10 6 3.1662 2.7511 3.7472 11 6 2.4606 2.0018 -1.2512 12 1 2.6828 3.0545 -1.0760 13 6 3.7186 1.2921 -1.7557 14 6 1.3690 1.8840 -2.2836 15 8 0.7737 0.8362 -2.4198 16 7 1.0546 2.9430 -3.0557 17 6 -0.0811 2.8656 -3.9778 18 6 -0.1457 4.1442 -4.8158 19 6 1.0120 4.3647 -5.7551 20 1 1.8741 4.9264 -5.4270 21 6 0.9704 3.8496 -7.0333 22 7 -0.1043 3.2186 -7.4534 23 14 2.4407 4.0477 -8.1684 24 1 3.5155 4.7895 -7.4614 25 1 2.9434 2.7081 -8.5657 26 1 2.0323 4.8026 -9.3803 27 6 -0.7784 5.3795 -4.1721 28 6 -1.5261 4.6212 -5.2717 29 1 -0.3633 0.5967 0.8368 30 1 -0.3634 0.4262 -0.9350 31 1 -0.3633 -1.0230 0.0985 32 1 1.8934 -0.5138 -0.8900 33 1 1.8935 -0.5138 0.8900 34 1 -0.1873 3.2698 3.1454 35 1 0.8964 4.1835 4.2222 36 1 1.2052 4.1486 2.4697 37 1 1.8319 0.6308 4.8228 38 1 1.2624 2.1291 5.5963 39 1 0.1785 1.2160 4.5192 40 1 3.5264 3.3281 2.8954 41 1 3.2173 3.3624 4.6482 42 1 3.7875 1.8645 3.8743 43 1 4.0227 1.7254 -2.7085 44 1 4.5211 1.4143 -1.0283 45 1 3.5081 0.2311 -1.8898 46 1 1.5803 3.7564 -3.0008 47 1 0.0431 2.0057 -4.6359 48 1 -1.0045 2.7572 -3.4088 49 1 -0.1073 2.7934 -8.3252 50 1 -1.1491 5.2814 -3.1518 51 1 -0.3577 6.3519 -4.4284 52 1 -1.5969 5.0948 -6.2509 53 1 -2.3883 4.0238 -4.9749 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 ZINC000861498172.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:01 Heat of formation + Delta-G solvation = -94.806833 kcal Electronic energy + Delta-G solvation = -31875.189238 eV Core-core repulsion = 27748.115846 eV Total energy + Delta-G solvation = -4127.073391 eV No. of doubly occupied orbitals = 66 Molecular weight (most abundant/longest-lived isotopes) = 340.220 amu Computer time = 0.70 seconds Orbital eigenvalues (eV) -41.03643 -39.57618 -37.59454 -36.70294 -35.71092 -34.10743 -33.49321 -32.05978 -31.87055 -29.48384 -27.60570 -27.53497 -26.70634 -25.39880 -24.13941 -23.58420 -21.80984 -21.13448 -19.85834 -19.49086 -18.96892 -17.91687 -17.13301 -16.62334 -16.17512 -15.99795 -15.65764 -15.11470 -14.74730 -14.63808 -14.45058 -14.29422 -14.01370 -13.69379 -13.42710 -13.29844 -13.01132 -12.73685 -12.66556 -12.47201 -12.31916 -12.15595 -12.07856 -12.06195 -11.95753 -11.71877 -11.59465 -11.55427 -11.38442 -11.28435 -11.02832 -10.94088 -10.66216 -10.53124 -10.45723 -10.24960 -10.02909 -9.66407 -9.47516 -9.09690 -8.70445 -8.47724 -8.30463 -8.10003 -6.09995 -4.19135 1.95965 2.77393 3.06781 3.23886 3.33157 3.35221 3.69990 3.84116 4.34318 4.53724 4.57913 4.72310 4.93680 5.02116 5.13652 5.17327 5.18920 5.26032 5.34412 5.37506 5.38722 5.48605 5.49947 5.57460 5.59188 5.62922 5.68545 5.76901 5.78353 5.93493 6.00049 6.04458 6.11373 6.12544 6.18344 6.31371 6.39094 6.42355 6.47945 6.48642 6.61453 6.69967 6.81892 7.10562 7.21068 7.34351 7.37802 7.53640 7.78989 8.16133 8.30330 8.31518 8.73362 8.82889 9.48299 10.97726 Molecular weight = 340.22amu Principal moments of inertia in cm(-1) A = 0.016636 B = 0.002832 C = 0.002763 Principal moments of inertia in units of 10**(-40)*gram-cm**2 A = 1682.672637 B = 9885.830252 C =10132.348389 Net atomic charges, atomic populations, and dipole contributions calculated with CM2 Atom NO. Type Charge No. of electrons 1 C -0.174 4.174 2 C 0.121 3.879 3 N -0.578 5.578 4 C 0.655 3.345 5 O -0.570 6.570 6 O -0.357 6.357 7 C 0.134 3.866 8 C -0.184 4.184 9 C -0.139 4.139 10 C -0.183 4.183 11 C 0.136 3.864 12 H 0.118 0.882 13 C -0.154 4.154 14 C 0.496 3.504 15 O -0.540 6.540 16 N -0.695 5.695 17 C 0.164 3.836 18 C 0.022 3.978 19 C -0.497 4.497 20 H 0.060 0.940 21 C 0.073 3.927 22 N -0.875 5.875 23 Si 0.889 3.111 24 H -0.272 1.272 25 H -0.283 1.283 26 H -0.284 1.284 27 C -0.213 4.213 28 C -0.138 4.138 29 H 0.045 0.955 30 H 0.087 0.913 31 H 0.058 0.942 32 H 0.092 0.908 33 H 0.068 0.932 34 H 0.059 0.941 35 H 0.068 0.932 36 H 0.089 0.911 37 H 0.065 0.935 38 H 0.079 0.921 39 H 0.066 0.934 40 H 0.088 0.912 41 H 0.068 0.932 42 H 0.058 0.942 43 H 0.079 0.921 44 H 0.067 0.933 45 H 0.068 0.932 46 H 0.401 0.599 47 H 0.072 0.928 48 H 0.036 0.964 49 H 0.261 0.739 50 H 0.064 0.936 51 H 0.078 0.922 52 H 0.088 0.912 53 H 0.062 0.938 For the charges calculated by CM2: Dipole moment (debyes) X Y Z Total from point charges 4.544 -4.220 20.621 21.533 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.232 4.232 2 C -0.001 4.001 3 N -0.318 5.318 4 C 0.410 3.590 5 O -0.460 6.460 6 O -0.271 6.271 7 C 0.097 3.903 8 C -0.241 4.241 9 C -0.197 4.197 10 C -0.240 4.240 11 C 0.030 3.970 12 H 0.136 0.864 13 C -0.212 4.212 14 C 0.280 3.720 15 O -0.419 6.419 16 N -0.344 5.344 17 C 0.041 3.959 18 C 0.021 3.979 19 C -0.534 4.534 20 H 0.078 0.922 21 C -0.131 4.131 22 N -0.511 5.511 23 Si 0.716 3.284 24 H -0.197 1.197 25 H -0.209 1.209 26 H -0.210 1.210 27 C -0.250 4.250 28 C -0.176 4.176 29 H 0.064 0.936 30 H 0.106 0.894 31 H 0.077 0.923 32 H 0.110 0.890 33 H 0.086 0.914 34 H 0.078 0.922 35 H 0.087 0.913 36 H 0.107 0.893 37 H 0.084 0.916 38 H 0.097 0.903 39 H 0.085 0.915 40 H 0.107 0.893 41 H 0.087 0.913 42 H 0.077 0.923 43 H 0.098 0.902 44 H 0.086 0.914 45 H 0.087 0.913 46 H 0.236 0.764 47 H 0.090 0.910 48 H 0.054 0.946 49 H 0.071 0.929 50 H 0.082 0.918 51 H 0.097 0.903 52 H 0.107 0.893 53 H 0.081 0.919 Dipole moment (debyes) X Y Z Total from point charges 4.092 -4.642 20.407 21.325 hybrid contribution 0.880 0.442 -0.688 1.201 sum 4.972 -4.200 19.719 20.766 Atomic orbital electron populations 1.22357 0.95419 1.01798 1.03615 1.21488 0.95364 0.80598 1.02608 1.46641 1.66441 1.14685 1.04031 1.17962 0.76894 0.83005 0.81109 1.90880 1.52268 1.16948 1.85924 1.86339 1.79301 1.45800 1.15705 1.22057 0.96502 0.90443 0.81273 1.22514 0.99677 0.97524 1.04377 1.21845 1.02193 1.00450 0.95177 1.22512 0.94977 1.02414 1.04068 1.20968 0.93496 1.00481 0.82048 0.86411 1.22057 0.95374 1.01673 1.02145 1.20557 0.84342 0.84442 0.82653 1.90579 1.53851 1.28039 1.69388 1.46010 1.32883 1.15630 1.39847 1.21041 0.85068 1.00291 0.89496 1.19054 0.94222 0.88605 0.96009 1.20389 1.01723 1.35558 0.95769 0.92155 1.24978 0.95495 0.95114 0.97499 1.70009 1.26649 1.34577 1.19901 0.86597 0.82985 0.77909 0.80882 1.19670 1.20948 1.21043 1.23051 1.08068 0.97988 0.95894 1.22446 0.94944 1.02845 0.97414 0.93554 0.89416 0.92293 0.89018 0.91357 0.92185 0.91268 0.89262 0.91563 0.90271 0.91477 0.89277 0.91292 0.92299 0.90173 0.91351 0.91264 0.76438 0.90992 0.94597 0.92949 0.91793 0.90333 0.89325 0.91912 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.17 -2.39 9.78 37.16 0.36 -2.03 16 2 C 0.12 1.57 5.01 -4.04 -0.02 1.55 16 3 N -0.58 -8.10 2.39 -173.18 -0.41 -8.51 16 4 C 0.66 9.48 7.54 54.06 0.41 9.89 16 5 O -0.57 -9.05 11.22 -19.28 -0.22 -9.27 16 6 O -0.36 -4.57 8.60 -40.31 -0.35 -4.92 16 7 C 0.13 1.39 1.13 -90.62 -0.10 1.29 16 8 C -0.18 -1.95 8.37 37.16 0.31 -1.64 16 9 C -0.14 -1.06 8.85 37.16 0.33 -0.73 16 10 C -0.18 -1.79 8.37 37.16 0.31 -1.48 16 11 C 0.14 1.95 3.04 -69.08 -0.21 1.74 16 12 H 0.12 1.68 6.81 -51.93 -0.35 1.33 16 13 C -0.15 -1.86 9.28 37.16 0.34 -1.52 16 14 C 0.50 8.63 6.80 -10.99 -0.07 8.56 16 15 O -0.54 -10.49 11.83 5.56 0.07 -10.42 16 16 N -0.70 -12.90 5.46 -60.29 -0.33 -13.23 16 17 C 0.16 3.57 5.53 -4.04 -0.02 3.54 16 18 C 0.02 0.50 1.50 -155.66 -0.23 0.27 16 19 C -0.50 -12.61 8.74 -34.44 -0.30 -12.91 16 20 H 0.06 1.51 7.82 -52.49 -0.41 1.10 16 21 C 0.07 1.91 5.42 -17.82 -0.10 1.82 16 22 N -0.88 -23.96 12.42 55.43 0.69 -23.28 16 23 Si 0.89 20.37 29.54 -169.99 -5.02 15.35 16 24 H -0.27 -6.12 7.11 56.52 0.40 -5.72 16 25 H -0.28 -6.54 7.11 56.52 0.40 -6.14 16 26 H -0.28 -6.56 7.11 56.52 0.40 -6.16 16 27 C -0.21 -4.41 9.96 -26.69 -0.27 -4.68 16 28 C -0.14 -3.15 8.99 -26.70 -0.24 -3.39 16 29 H 0.05 0.64 7.77 -51.93 -0.40 0.23 16 30 H 0.09 1.48 5.30 -51.93 -0.28 1.21 16 31 H 0.06 0.71 8.14 -51.93 -0.42 0.28 16 32 H 0.09 1.21 5.98 -51.93 -0.31 0.90 16 33 H 0.07 0.77 7.70 -51.93 -0.40 0.37 16 34 H 0.06 0.64 8.14 -51.93 -0.42 0.21 16 35 H 0.07 0.58 8.14 -51.93 -0.42 0.16 16 36 H 0.09 1.18 5.88 -51.93 -0.31 0.87 16 37 H 0.07 0.47 8.14 -51.93 -0.42 0.05 16 38 H 0.08 0.45 8.14 -51.93 -0.42 0.03 16 39 H 0.07 0.51 8.14 -51.93 -0.42 0.09 16 40 H 0.09 1.09 5.88 -51.93 -0.31 0.78 16 41 H 0.07 0.54 8.14 -51.93 -0.42 0.11 16 42 H 0.06 0.54 8.14 -51.93 -0.42 0.12 16 43 H 0.08 0.95 8.14 -51.93 -0.42 0.53 16 44 H 0.07 0.73 8.14 -51.93 -0.42 0.31 16 45 H 0.07 0.84 6.94 -51.93 -0.36 0.48 16 46 H 0.40 6.95 8.38 -40.82 -0.34 6.61 16 47 H 0.07 1.79 8.10 -51.93 -0.42 1.37 16 48 H 0.04 0.77 8.09 -51.93 -0.42 0.35 16 49 H 0.26 6.87 8.32 -40.82 -0.34 6.53 16 50 H 0.06 1.22 8.14 -51.93 -0.42 0.80 16 51 H 0.08 1.62 8.14 -51.93 -0.42 1.19 16 52 H 0.09 2.19 8.03 -51.93 -0.42 1.77 16 53 H 0.06 1.36 8.14 -51.92 -0.42 0.93 16 LS Contribution 417.94 15.07 6.30 6.30 Total: -1.00 -28.86 417.94 -8.13 -36.99 By element: Atomic # 1 Polarization: 20.07 SS G_CDS: -9.36 Total: 10.71 kcal Atomic # 6 Polarization: -0.22 SS G_CDS: 0.50 Total: 0.28 kcal Atomic # 7 Polarization: -44.96 SS G_CDS: -0.05 Total: -45.02 kcal Atomic # 8 Polarization: -24.11 SS G_CDS: -0.50 Total: -24.61 kcal Atomic # 14 Polarization: 20.37 SS G_CDS: -5.02 Total: 15.35 kcal Total LS contribution 6.30 Total: 6.30 kcal Total: -28.86 -8.13 -36.99 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. ZINC000861498172.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.821 kcal (2) G-P(sol) polarization free energy of solvation -28.857 kcal (3) G-ENP(sol) elect.-nuc.-pol. free energy of system -86.678 kcal (4) G-CDS(sol) cavity-dispersion-solvent structure free energy -8.129 kcal (5) G-P-CDS(sol) = G-P(sol) + G-CDS(sol) = (2) + (4) -36.986 kcal (6) G-S(sol) free energy of system = (1) + (5) -94.807 kcal 1SCF WAS SPECIFIED, SO GEOMETRY OPTIMIZATION WAS NOT USED Total computer time = 0.70 seconds