Wall clock time and date at job start Fri Mar 5 2021 09:44:45 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 * 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.17394 * 1 3 3 C 1.43197 * 179.97438 * 2 1 4 4 C 1.39715 * 119.99905 * 278.37980 * 3 2 1 5 5 C 1.37596 * 120.00510 * 180.02562 * 4 3 2 6 6 C 1.39747 * 119.99691 * 359.97130 * 5 4 3 7 7 C 1.47693 * 120.00455 * 180.02562 * 6 5 4 8 8 O 1.21566 * 119.99841 * 0.02562 * 7 6 5 9 9 N 1.34775 * 119.99904 * 180.02562 * 7 6 5 10 10 C 1.46497 * 119.99724 * 180.02562 * 9 7 6 11 11 C 1.53001 * 109.46985 * 179.97438 * 10 9 7 12 12 C 1.53514 * 114.17088 * 184.99390 * 11 10 9 13 13 N 1.47479 * 87.93405 * 220.14043 * 12 11 10 14 14 C 1.46894 * 110.99953 * 138.70663 * 13 12 11 15 15 C 1.50703 * 109.47274 * 170.00304 * 14 13 12 16 16 C 1.37946 * 119.93063 * 89.99685 * 15 14 13 17 17 C 1.39556 * 119.85442 * 180.02562 * 16 15 14 18 18 C 1.43205 * 120.14159 * 180.02562 * 17 16 15 19 19 C 9.88275 * 109.10414 * 144.07239 * 7 1 2 20 20 C 1.39526 * 119.71906 * 0.27677 * 17 16 15 21 21 C 1.37992 * 119.84926 * 359.42184 * 20 17 16 22 22 C 1.50705 * 119.92612 * 180.31689 * 21 20 17 23 23 O 1.42901 * 109.46865 * 120.26171 * 22 21 20 24 24 Si 1.86304 * 109.46889 * 0.26306 * 22 21 20 25 25 H 1.48498 * 109.47201 * 59.99984 * 24 22 21 26 26 H 1.48502 * 109.46772 * 179.97438 * 24 22 21 27 27 H 1.48500 * 109.47111 * 299.99227 * 24 22 21 28 28 C 1.38321 * 120.14671 * 0.58069 * 21 20 17 29 29 C 1.47479 * 87.83553 * 27.01207 * 13 12 11 30 30 C 1.39745 * 119.99368 * 0.28098 * 6 5 4 31 31 C 1.37603 * 119.99192 * 359.48875 * 30 6 5 32 32 H 4.51954 * 14.06437 * 359.97438 * 4 2 1 33 33 H 1.08001 * 119.99836 * 359.97438 * 4 3 2 34 34 H 1.07999 * 120.00414 * 179.97438 * 5 4 3 35 35 H 0.96999 * 120.00108 * 359.97438 * 9 7 6 36 36 H 1.08996 * 109.47283 * 59.99774 * 10 9 7 37 37 H 1.09001 * 109.47291 * 299.99345 * 10 9 7 38 38 H 1.09000 * 113.25795 * 317.94974 * 11 10 9 39 39 H 1.09007 * 113.47206 * 105.37598 * 12 11 10 40 40 H 1.09001 * 113.47873 * 334.91151 * 12 11 10 41 41 H 1.09002 * 109.47576 * 50.00004 * 14 13 12 42 42 H 1.09004 * 109.47351 * 289.99958 * 14 13 12 43 43 H 1.07995 * 120.07468 * 0.02562 * 16 15 14 44 44 H 9.56436 * 106.99094 * 237.79909 * 9 1 2 45 45 H 1.07997 * 120.07215 * 179.69882 * 20 17 16 46 46 H 1.08998 * 109.47044 * 240.26153 * 22 21 20 47 47 H 0.96706 * 114.00080 * 299.99750 * 23 22 21 48 48 H 1.07992 * 119.86036 * 179.74035 * 28 21 20 49 49 H 1.09001 * 113.47651 * 87.75729 * 29 13 12 50 50 H 1.09004 * 113.47517 * 218.22377 * 29 13 12 51 51 H 1.07997 * 120.00066 * 179.70362 * 30 6 5 52 52 H 1.07999 * 119.99681 * 180.23483 * 31 30 6 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.1739 0.0000 0.0000 3 6 2.6059 0.0006 0.0000 4 6 3.3045 -0.1754 1.1971 5 6 4.6805 -0.1743 1.1973 6 6 5.3792 0.0024 0.0000 7 6 6.8562 0.0030 -0.0001 8 8 7.4641 -0.1505 1.0414 9 7 7.5299 0.1739 -1.1547 10 6 8.9949 0.1740 -1.1548 11 6 9.5048 0.3845 -2.5819 12 6 11.0297 0.2826 -2.7260 13 7 10.9641 1.4219 -3.6601 14 6 12.0669 2.3625 -3.4215 15 6 12.1144 3.3744 -4.5372 16 6 11.3967 4.5479 -4.4331 17 6 11.4414 5.4836 -5.4675 18 6 10.6992 6.7041 -5.3657 19 6 10.0873 7.5663 -5.4795 20 6 12.2082 5.2227 -6.6036 21 6 12.9264 4.0478 -6.6926 22 6 13.7608 3.7675 -7.9158 23 8 15.1277 3.6039 -7.5327 24 14 13.6159 5.1931 -9.1064 25 1 14.0987 6.4356 -8.4520 26 1 14.4386 4.9171 -10.3115 27 1 12.1954 5.3630 -9.5046 28 6 12.8779 3.1244 -5.6639 29 6 9.7403 1.8537 -2.9595 30 6 4.6804 0.1837 -1.1965 31 6 3.3043 0.1778 -1.1970 32 1 -1.0501 0.0005 0.0002 33 1 2.7646 -0.3115 2.1224 34 1 5.2205 -0.3099 2.1227 35 1 7.0448 0.2960 -1.9858 36 1 9.3583 -0.7815 -0.7768 37 1 9.3582 0.9793 -0.5164 38 1 8.9471 -0.1837 -3.3264 39 1 11.3606 -0.6444 -3.1944 40 1 11.5688 0.5062 -1.8054 41 1 13.0094 1.8160 -3.3863 42 1 11.9085 2.8751 -2.4727 43 1 10.8016 4.7420 -3.5532 44 1 9.5448 8.6085 -5.1888 45 1 12.2430 5.9407 -7.4096 46 1 13.4064 2.8554 -8.3961 47 1 15.5099 4.3787 -7.0981 48 1 13.4385 2.2045 -5.7402 49 1 9.9352 2.4968 -2.1012 50 1 8.9783 2.2583 -3.6257 51 1 5.2202 0.3250 -2.1211 52 1 2.7642 0.3141 -2.1222 RHF calculation, no. of doubly occupied orbitals= 71 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=0 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 m_275030_14113930_14126694_16591380.mol2 52 1SCF WAS SPECIFIED, SO GEOMETRY OPTIMIZATION WAS NOT USED SCF FIELD WAS ACHIEVED Wall clock time and date at job start Fri Mar 5 2021 09:44:45 Heat of formation + Delta-G solvation = 165.319609 kcal Electronic energy + Delta-G solvation = -32864.055861 eV Core-core repulsion = 28458.603834 eV Total energy + Delta-G solvation = -4405.452027 eV No. of doubly occupied orbitals = 71 Molecular weight (most abundant/longest-lived isotopes) = 388.170 amu Computer time = 1.24 seconds Orbital eigenvalues (eV) -41.78876 -41.15190 -40.56167 -39.43319 -37.84126 -36.33319 -35.81504 -34.89163 -33.93297 -32.27772 -31.93031 -31.63527 -30.53965 -27.83723 -27.20488 -26.75933 -25.51049 -24.38861 -23.81558 -23.44036 -21.92535 -21.65764 -20.86511 -20.05472 -19.31174 -18.88721 -17.84498 -17.52082 -17.37718 -17.01156 -16.81572 -16.74705 -16.46654 -16.35631 -16.18387 -16.00131 -15.80335 -15.64151 -15.47296 -14.92833 -14.78587 -14.73898 -14.53774 -14.44122 -14.21475 -14.11443 -13.85287 -13.69865 -13.32741 -13.14136 -13.02812 -12.77145 -12.70352 -12.49044 -12.37446 -12.27814 -12.19921 -12.16295 -12.01214 -11.84339 -11.66909 -11.63985 -11.51480 -11.00218 -10.47071 -10.24259 -10.16510 -10.09837 -9.77957 -9.63847 -9.57251 -0.56455 0.02369 0.19574 0.25068 1.12389 1.41240 1.47281 1.49292 1.69170 1.89145 1.99920 2.07544 2.27152 2.27606 2.82804 2.93714 3.02200 3.13614 3.16828 3.21461 3.40674 3.49381 3.61980 3.62677 3.68164 3.69444 3.81636 3.90563 3.99440 4.02086 4.13626 4.16058 4.18180 4.21894 4.32846 4.36766 4.38766 4.44819 4.47256 4.53052 4.54449 4.57112 4.62987 4.70679 4.81606 4.86467 4.92149 4.99577 5.01303 5.10813 5.15252 5.23601 5.30076 5.32527 5.36050 5.36921 5.52830 5.69155 5.81590 5.89658 6.37862 6.46832 6.49779 6.51296 7.13467 Net atomic charges, atomic populations, and dipole contributions calculated with CM2 Atom NO. Type Charge No. of electrons 1 C -0.154 4.154 2 C -0.165 4.165 3 C 0.040 3.960 4 C -0.104 4.104 5 C -0.056 4.056 6 C -0.130 4.130 7 C 0.559 3.441 8 O -0.533 6.533 9 N -0.725 5.725 10 C 0.130 3.870 11 C -0.130 4.130 12 C 0.008 3.992 13 N -0.463 5.463 14 C 0.057 3.943 15 C -0.073 4.073 16 C -0.081 4.081 17 C -0.004 4.004 18 C -0.139 4.139 19 C -0.174 4.174 20 C -0.103 4.103 21 C -0.119 4.119 22 C -0.077 4.077 23 O -0.537 6.537 24 Si 0.944 3.056 25 H -0.251 1.251 26 H -0.251 1.251 27 H -0.250 1.250 28 C -0.061 4.061 29 C 0.001 3.999 30 C -0.084 4.084 31 C -0.104 4.104 32 H 0.227 0.773 33 H 0.136 0.864 34 H 0.143 0.857 35 H 0.402 0.598 36 H 0.076 0.924 37 H 0.072 0.928 38 H 0.111 0.889 39 H 0.093 0.907 40 H 0.076 0.924 41 H 0.088 0.912 42 H 0.059 0.941 43 H 0.124 0.876 44 H 0.244 0.756 45 H 0.122 0.878 46 H 0.121 0.879 47 H 0.377 0.623 48 H 0.131 0.869 49 H 0.074 0.926 50 H 0.087 0.913 51 H 0.131 0.869 52 H 0.135 0.865 For the charges calculated by CM2: Dipole moment (debyes) X Y Z Total from point charges -2.300 1.378 -1.343 2.999 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.170 4.170 2 C -0.165 4.165 3 C 0.040 3.960 4 C -0.122 4.122 5 C -0.075 4.075 6 C -0.133 4.133 7 C 0.346 3.654 8 O -0.410 6.410 9 N -0.378 5.378 10 C 0.006 3.994 11 C -0.149 4.149 12 C -0.121 4.121 13 N -0.185 5.185 14 C -0.071 4.071 15 C -0.073 4.073 16 C -0.099 4.099 17 C -0.004 4.004 18 C -0.139 4.139 19 C -0.190 4.190 20 C -0.121 4.121 21 C -0.119 4.119 22 C -0.188 4.188 23 O -0.342 6.342 24 Si 0.767 3.233 25 H -0.175 1.175 26 H -0.175 1.175 27 H -0.174 1.174 28 C -0.079 4.079 29 C -0.128 4.128 30 C -0.102 4.102 31 C -0.122 4.122 32 H 0.244 0.756 33 H 0.154 0.846 34 H 0.161 0.839 35 H 0.237 0.763 36 H 0.095 0.905 37 H 0.090 0.910 38 H 0.129 0.871 39 H 0.112 0.888 40 H 0.094 0.906 41 H 0.106 0.894 42 H 0.077 0.923 43 H 0.142 0.858 44 H 0.260 0.740 45 H 0.140 0.860 46 H 0.138 0.862 47 H 0.224 0.776 48 H 0.149 0.851 49 H 0.092 0.908 50 H 0.106 0.894 51 H 0.149 0.851 52 H 0.153 0.847 Dipole moment (debyes) X Y Z Total from point charges -1.923 1.285 -2.159 3.164 hybrid contribution -0.374 0.589 0.977 1.200 sum -2.296 1.875 -1.182 3.191 Atomic orbital electron populations 1.25043 0.95886 0.97697 0.98413 1.19962 0.92013 1.02410 1.02149 1.16195 0.86862 0.99077 0.93909 1.21245 0.93773 0.99576 0.97570 1.21223 0.93453 0.94834 0.97975 1.20025 0.92794 1.05743 0.94771 1.18232 0.87144 0.76410 0.83598 1.90811 1.69832 1.48861 1.31484 1.45880 1.07414 1.74383 1.10084 1.21370 0.79255 1.02943 0.95810 1.23008 0.97315 0.97879 0.96706 1.23819 0.94963 0.97118 0.96207 1.68984 1.02648 1.05460 1.41405 1.21663 0.96745 0.93487 0.95194 1.19559 0.97349 0.94428 0.96008 1.21058 0.98126 0.92784 0.97951 1.16160 0.97807 0.91545 0.94915 1.16998 0.98940 0.96555 1.01363 1.25338 0.97331 0.99238 0.97073 1.21076 0.96934 0.96820 0.97286 1.20187 0.99468 0.94699 0.97566 1.24938 0.83131 1.11654 0.99091 1.86520 1.23076 1.43369 1.81242 0.85759 0.81489 0.77010 0.79076 1.17458 1.17533 1.17378 1.20943 0.96810 0.98501 0.91692 1.24082 0.94842 0.96016 0.97848 1.21192 0.94136 0.97816 0.97101 1.21155 0.93758 0.99755 0.97524 0.75597 0.84622 0.83910 0.76261 0.90549 0.90999 0.87054 0.88830 0.90596 0.89353 0.92322 0.85847 0.73990 0.85998 0.86152 0.77626 0.85080 0.90802 0.89434 0.85075 0.84740 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 0.03 19.25 29.55 0.57 0.60 16 2 C -0.17 -0.23 13.23 -38.84 -0.51 -0.75 16 3 C 0.04 0.04 5.73 -106.87 -0.61 -0.57 16 4 C -0.10 -0.07 9.73 -39.26 -0.38 -0.46 16 5 C -0.06 -0.07 9.57 -39.25 -0.38 -0.44 16 6 C -0.13 -0.12 5.88 -104.91 -0.62 -0.73 16 7 C 0.56 1.11 7.77 -12.40 -0.10 1.01 16 8 O -0.53 -3.06 16.33 5.31 0.09 -2.97 16 9 N -0.72 0.84 5.50 -61.30 -0.34 0.50 16 10 C 0.13 -0.20 5.59 -4.04 -0.02 -0.22 16 11 C -0.13 0.37 4.08 -90.09 -0.37 0.01 16 12 C 0.01 -0.02 7.26 -7.55 -0.05 -0.07 16 13 N -0.46 0.46 6.32 -246.08 -1.55 -1.10 16 14 C 0.06 -0.06 4.94 -4.98 -0.02 -0.09 16 15 C -0.07 -0.03 5.46 -104.68 -0.57 -0.61 16 16 C -0.08 -0.09 9.40 -39.20 -0.37 -0.46 16 17 C 0.00 -0.01 5.72 -107.00 -0.61 -0.62 16 18 C -0.14 -0.30 12.65 -40.38 -0.51 -0.81 16 19 C -0.17 -0.06 18.54 28.01 0.52 0.46 16 20 C -0.10 -0.18 8.37 -39.19 -0.33 -0.50 16 21 C -0.12 -0.10 5.39 -104.68 -0.56 -0.66 16 22 C -0.08 0.02 2.82 36.01 0.10 0.12 16 23 O -0.54 -0.11 12.77 -35.23 -0.45 -0.56 16 24 Si 0.94 -0.05 25.07 -169.99 -4.26 -4.31 16 25 H -0.25 -0.07 7.11 56.52 0.40 0.33 16 26 H -0.25 -0.29 7.11 56.52 0.40 0.12 16 27 H -0.25 -0.19 7.11 56.52 0.40 0.22 16 28 C -0.06 -0.05 9.25 -39.50 -0.37 -0.41 16 29 C 0.00 0.00 7.29 -7.55 -0.06 -0.06 16 30 C -0.08 0.09 9.52 -39.25 -0.37 -0.28 16 31 C -0.10 0.04 9.73 -39.26 -0.38 -0.35 16 32 H 0.23 -0.86 7.80 -54.09 -0.42 -1.29 16 33 H 0.14 -0.10 8.06 -52.49 -0.42 -0.52 16 34 H 0.14 0.16 7.64 -52.49 -0.40 -0.24 16 35 H 0.40 -1.46 6.49 -40.82 -0.26 -1.72 16 36 H 0.08 -0.05 8.14 -51.93 -0.42 -0.48 16 37 H 0.07 -0.05 7.82 -51.93 -0.41 -0.46 16 38 H 0.11 -0.35 8.14 -51.93 -0.42 -0.77 16 39 H 0.09 -0.23 8.14 -51.93 -0.42 -0.65 16 40 H 0.08 -0.23 7.87 -51.93 -0.41 -0.64 16 41 H 0.09 -0.14 8.02 -51.93 -0.42 -0.56 16 42 H 0.06 -0.10 7.24 -51.93 -0.38 -0.48 16 43 H 0.12 0.00 8.06 -52.49 -0.42 -0.43 16 44 H 0.24 -1.16 9.09 -44.17 -0.40 -1.56 16 45 H 0.12 0.13 4.87 -52.49 -0.26 -0.13 16 46 H 0.12 -0.12 7.96 -51.93 -0.41 -0.53 16 47 H 0.38 -1.08 8.70 45.56 0.40 -0.69 16 48 H 0.13 0.02 8.06 -52.49 -0.42 -0.40 16 49 H 0.07 -0.17 7.16 -51.93 -0.37 -0.54 16 50 H 0.09 -0.15 8.14 -51.93 -0.42 -0.57 16 51 H 0.13 -0.54 6.38 -52.49 -0.34 -0.87 16 52 H 0.13 -0.27 8.06 -52.49 -0.42 -0.69 16 LS Contribution 446.34 15.07 6.73 6.73 Total: 0.00 -9.11 446.34 -12.05 -21.16 By element: Atomic # 1 Polarization: -7.30 SS G_CDS: -6.26 Total: -13.56 kcal Atomic # 6 Polarization: 0.11 SS G_CDS: -6.01 Total: -5.89 kcal Atomic # 7 Polarization: 1.30 SS G_CDS: -1.89 Total: -0.59 kcal Atomic # 8 Polarization: -3.17 SS G_CDS: -0.36 Total: -3.53 kcal Atomic # 14 Polarization: -0.05 SS G_CDS: -4.26 Total: -4.31 kcal Total LS contribution 6.73 Total: 6.73 kcal Total: -9.11 -12.05 -21.16 kcal The number of atoms in the molecule is 52 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. m_275030_14113930_14126694_16591380.mol2 52 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 186.481 kcal (2) G-P(sol) polarization free energy of solvation -9.108 kcal (3) G-ENP(sol) elect.-nuc.-pol. free energy of system 177.373 kcal (4) G-CDS(sol) cavity-dispersion-solvent structure free energy -12.053 kcal (5) G-P-CDS(sol) = G-P(sol) + G-CDS(sol) = (2) + (4) -21.161 kcal (6) G-S(sol) free energy of system = (1) + (5) 165.320 kcal 1SCF WAS SPECIFIED, SO GEOMETRY OPTIMIZATION WAS NOT USED Total computer time = 1.24 seconds