Wall clock time and date at job start Tue Mar 9 2021 13:39:29 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 * SM5.42R - SM5.42R CALCULATIONS WILL BE PERFORMED * - CHARGE MODEL 2 WILL BE USED (CM2A) * - THE SOLVENT IS WATER ******************************************************************************* 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 O 2 2 C 1.42904 * 1 3 3 Si 1.86300 * 109.47298 * 2 1 4 4 C 1.50704 * 109.47024 * 120.00582 * 2 1 3 5 5 C 1.38260 * 119.94688 * 64.99703 * 4 2 1 6 6 C 1.38145 * 120.05447 * 179.97438 * 5 4 2 7 7 C 1.38726 * 119.94375 * 0.02562 * 6 5 4 8 8 O 1.35904 * 120.05571 * 180.02562 * 7 6 5 9 9 C 1.42896 * 117.00327 * 179.97438 * 8 7 6 10 10 C 1.50706 * 109.47239 * 179.97438 * 9 8 7 11 11 C 1.36414 * 125.79641 * 90.00303 * 10 9 8 12 12 C 1.47375 * 127.19495 * 359.96580 * 11 10 9 13 13 O 1.21620 * 120.00114 * 359.97438 * 12 11 10 14 14 N 1.34775 * 120.00179 * 179.97438 * 12 11 10 15 15 C 1.46503 * 120.00094 * 179.97438 * 14 12 11 16 16 C 1.50696 * 112.85356 * 309.67541 * 15 14 12 17 17 C 1.38435 * 120.01433 * 170.28715 * 16 15 14 18 18 F 1.35093 * 120.02480 * 0.04455 * 17 16 15 19 19 C 1.38490 * 119.94736 * 180.02562 * 17 16 15 20 20 C 1.38175 * 119.97331 * 359.97438 * 19 17 16 21 21 C 1.38263 * 120.02891 * 0.02562 * 20 19 17 22 22 C 1.38221 * 120.01038 * 350.32447 * 16 15 14 23 23 C 1.53782 * 113.61372 * 178.47373 * 15 14 12 24 24 C 1.53781 * 87.08058 * 270.80584 * 23 15 14 25 25 H 1.09007 * 113.61032 * 219.93797 * 24 23 15 26 26 O 1.42901 * 113.61900 * 89.11751 * 24 23 15 27 27 C 1.53779 * 113.61052 * 80.86668 * 15 14 12 28 28 C 1.47021 * 105.61405 * 180.22452 * 11 10 9 29 29 C 1.34179 * 105.95074 * 359.63026 * 28 11 10 30 30 O 1.33653 * 125.80282 * 270.00311 * 10 9 8 31 31 C 1.38723 * 119.89260 * 0.27105 * 7 6 5 32 32 C 1.38155 * 119.94270 * 359.47008 * 31 7 6 33 33 H 0.96698 * 114.00045 * 59.99790 * 1 2 3 34 34 H 1.09009 * 109.46900 * 240.00000 * 2 1 3 35 35 H 1.48502 * 109.99610 * 60.00176 * 3 2 1 36 36 H 1.48495 * 109.99791 * 181.31881 * 3 2 1 37 37 H 1.08001 * 119.97146 * 359.97438 * 5 4 2 38 38 H 1.08002 * 120.02659 * 180.02562 * 6 5 4 39 39 H 1.09003 * 109.46670 * 60.00034 * 9 8 7 40 40 H 1.08997 * 109.47218 * 300.00025 * 9 8 7 41 41 H 0.96999 * 119.99894 * 0.02562 * 14 12 11 42 42 H 1.08001 * 120.00780 * 180.24926 * 19 17 16 43 43 H 1.07998 * 119.98101 * 179.97438 * 20 19 17 44 44 H 1.08000 * 119.97593 * 179.97438 * 21 20 19 45 45 H 1.08002 * 119.98827 * 359.96689 * 22 16 15 46 46 H 1.09005 * 113.61250 * 25.35234 * 23 15 14 47 47 H 1.08995 * 113.61507 * 156.16917 * 23 15 14 48 48 H 0.96695 * 114.00296 * 179.97438 * 26 24 23 49 49 H 1.09000 * 113.61624 * 203.74549 * 27 15 14 50 50 H 1.09006 * 113.61596 * 334.65493 * 27 15 14 51 51 H 1.07995 * 127.02619 * 179.73938 * 28 11 10 52 52 H 1.08008 * 125.43712 * 180.12309 * 29 28 11 53 53 H 1.07997 * 120.03567 * 179.72444 * 31 7 6 54 54 H 1.07995 * 119.97214 * 180.22880 * 32 31 7 Note: An asterisk (*) indicates this parameter was or will be optimized. Cartesian coordinates (angstroms) NO. Atom X Y Z 1 8 0.0000 0.0000 0.0000 2 6 1.4290 0.0000 0.0000 3 14 2.0501 1.7564 0.0000 4 6 1.9314 -0.7106 1.2304 5 6 1.6840 -2.0606 1.3971 6 6 2.1439 -2.7143 2.5239 7 6 2.8540 -2.0148 3.4888 8 8 3.3079 -2.6557 4.5980 9 6 4.0319 -1.8727 5.5491 10 6 4.4517 -2.7461 6.7033 11 6 5.6192 -3.4463 6.7907 12 6 6.7027 -3.4999 5.7932 13 8 6.6057 -2.8728 4.7556 14 7 7.8011 -4.2438 6.0306 15 6 8.8780 -4.2975 5.0388 16 6 8.3809 -4.6135 3.6517 17 6 9.2355 -4.5038 2.5681 18 9 10.5157 -4.1143 2.7536 19 6 8.7770 -4.7940 1.2940 20 6 7.4676 -5.1931 1.1057 21 6 6.6150 -5.3024 2.1886 22 6 7.0713 -5.0122 3.4608 23 6 10.0751 -5.1586 5.4752 24 6 10.6569 -3.8863 6.1136 25 1 11.7338 -3.7836 5.9792 26 8 10.2450 -3.6850 7.4670 27 6 9.8259 -3.0878 5.0927 28 6 5.5678 -4.1414 8.0852 29 6 4.3995 -3.7921 8.6452 30 8 3.7435 -2.9663 7.8152 31 6 3.0955 -0.6590 3.3222 32 6 2.6381 -0.0109 2.1911 33 1 -0.3933 0.4417 -0.7650 34 1 1.7924 -0.5139 -0.8901 35 1 1.5616 2.4678 1.2086 36 1 3.5346 1.7702 -0.0321 37 1 1.1308 -2.6046 0.6458 38 1 1.9508 -3.7690 2.6535 39 1 3.3958 -1.0672 5.9160 40 1 4.9163 -1.4492 5.0731 41 1 7.8784 -4.7443 6.8580 42 1 9.4439 -4.7126 0.4484 43 1 7.1099 -5.4191 0.1120 44 1 5.5917 -5.6140 2.0405 45 1 6.4044 -5.0978 4.3060 46 1 9.8163 -5.9312 6.1992 47 1 10.6654 -5.5399 4.6421 48 1 10.6080 -2.8865 7.8741 49 1 10.3489 -2.9008 4.1548 50 1 9.3663 -2.1911 5.5084 51 1 6.3179 -4.7975 8.5013 52 1 4.0421 -4.1257 9.6082 53 1 3.6442 -0.1113 4.0741 54 1 2.8291 1.0438 2.0591 There are 80 doubly filled levels RHF calculation, no. of doubly occupied orbitals= 80 No. of singly occupied orbitals= 1 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) F: (AM1): M.J.S. DEWAR AND E. G. ZOEBISCH, THEOCHEM, 180, 1 (1988). 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=WATER s_22_17639084_4549634.mol2 54 1SCF WAS SPECIFIED, SO GEOMETRY OPTIMIZATION WAS NOT USED SCF FIELD WAS ACHIEVED Wall clock time and date at job start Tue Mar 9 2021 13:39:29 Heat of formation + Delta-G solvation = -62.643804 kcal Electronic energy + Delta-G solvation = -44960.214167 eV Core-core repulsion = 39326.684143 eV Total energy + Delta-G solvation = -5633.530024 eV No. of doubly occupied orbitals = 80 and no. of open levels = 1 Molecular weight (most abundant/longest-lived isotopes) = 440.156 amu Computer time = 1.93 seconds Orbital eigenvalues (eV) -50.16606 -43.80487 -42.05332 -41.41629 -39.66416 -39.44396 -39.21040 -37.26086 -36.30805 -34.59326 -33.62583 -32.93874 -32.15125 -32.03647 -31.56421 -31.37952 -28.11255 -27.08582 -26.38181 -26.05146 -24.71867 -23.98160 -23.69633 -22.93866 -22.39121 -22.14554 -21.05516 -20.62087 -19.65321 -19.30784 -18.80592 -18.36872 -18.22279 -17.74800 -17.66625 -17.33952 -17.09936 -17.05268 -16.90631 -16.80936 -16.73754 -16.64218 -16.09464 -15.96389 -15.76511 -15.57109 -15.41058 -15.21308 -14.83474 -14.71055 -14.56606 -14.40440 -14.32327 -14.03332 -13.99511 -13.88645 -13.82767 -13.68172 -13.58649 -13.49008 -13.16257 -13.02924 -12.56464 -12.55626 -12.40645 -12.29227 -12.07205 -11.97663 -11.91009 -11.67817 -11.22587 -11.02368 -10.87672 -10.30543 -10.27796 -10.21316 -9.92731 -9.86272 -9.63523 -9.37123 -4.98346 -0.11996 0.08683 0.08799 0.14355 0.20825 0.97324 1.28880 1.41706 1.61171 1.86763 1.94861 2.02762 2.35688 2.47417 2.52505 2.98247 3.17755 3.24380 3.30095 3.39186 3.40437 3.47569 3.51957 3.58666 3.65330 3.68643 3.70648 3.80003 3.83386 3.91289 3.99325 4.00075 4.11583 4.22569 4.28387 4.28672 4.36325 4.38497 4.39555 4.54756 4.55843 4.56712 4.65659 4.73226 4.77643 4.78264 4.81092 4.83237 4.86859 4.97862 5.00072 5.04063 5.06794 5.20175 5.26803 5.32443 5.40102 5.48583 5.53074 5.81494 5.85327 6.10467 6.37430 6.76454 6.92744 7.05374 Net atomic charges, atomic populations, and dipole contributions calculated with CM2 Atom NO. Type Charge No. of electrons 1 O -0.563 6.563 2 C -0.070 4.070 3 Si 0.770 3.230 4 C -0.128 4.128 5 C -0.051 4.051 6 C -0.153 4.153 7 C 0.128 3.872 8 O -0.314 6.314 9 C 0.137 3.863 10 C 0.041 3.959 11 C -0.219 4.219 12 C 0.584 3.416 13 O -0.552 6.552 14 N -0.678 5.678 15 C 0.195 3.805 16 C -0.114 4.114 17 C 0.088 3.912 18 F -0.135 7.135 19 C -0.159 4.159 20 C -0.091 4.091 21 C -0.137 4.137 22 C -0.069 4.069 23 C -0.121 4.121 24 C 0.095 3.905 25 H 0.092 0.908 26 O -0.559 6.559 27 C -0.186 4.186 28 C -0.197 4.197 29 C -0.009 4.009 30 O -0.177 6.177 31 C -0.209 4.209 32 C -0.100 4.100 33 H 0.403 0.597 34 H 0.096 0.904 35 H -0.297 1.297 36 H -0.293 1.293 37 H 0.154 0.846 38 H 0.146 0.854 39 H 0.093 0.907 40 H 0.104 0.896 41 H 0.417 0.583 42 H 0.154 0.846 43 H 0.156 0.844 44 H 0.136 0.864 45 H 0.112 0.888 46 H 0.095 0.905 47 H 0.116 0.884 48 H 0.398 0.602 49 H 0.110 0.890 50 H 0.061 0.939 51 H 0.188 0.812 52 H 0.253 0.747 53 H 0.139 0.861 54 H 0.119 0.881 For the charges calculated by CM2: Dipole moment (debyes) X Y Z Total from point charges 1.009 -1.420 -0.057 1.743 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 O -0.373 6.373 2 C -0.180 4.180 3 Si 0.677 3.323 4 C -0.128 4.128 5 C -0.069 4.069 6 C -0.172 4.172 7 C 0.083 3.917 8 O -0.228 6.228 9 C 0.062 3.938 10 C -0.012 4.012 11 C -0.224 4.224 12 C 0.373 3.627 13 O -0.432 6.432 14 N -0.330 5.330 15 C 0.108 3.892 16 C -0.116 4.116 17 C 0.068 3.932 18 F -0.114 7.114 19 C -0.178 4.178 20 C -0.109 4.109 21 C -0.155 4.155 22 C -0.087 4.087 23 C -0.158 4.158 24 C 0.035 3.965 25 H 0.110 0.890 26 O -0.367 6.367 27 C -0.224 4.224 28 C -0.216 4.216 29 C -0.074 4.074 30 O -0.072 6.072 31 C -0.227 4.227 32 C -0.117 4.117 33 H 0.255 0.745 34 H 0.114 0.886 35 H -0.225 1.225 36 H -0.221 1.221 37 H 0.172 0.828 38 H 0.164 0.836 39 H 0.111 0.889 40 H 0.122 0.878 41 H 0.255 0.745 42 H 0.171 0.829 43 H 0.174 0.826 44 H 0.154 0.846 45 H 0.130 0.870 46 H 0.114 0.886 47 H 0.134 0.866 48 H 0.249 0.751 49 H 0.128 0.872 50 H 0.079 0.921 51 H 0.205 0.795 52 H 0.269 0.731 53 H 0.157 0.843 54 H 0.137 0.863 Dipole moment (debyes) X Y Z Total from point charges 0.950 -1.318 0.218 1.639 hybrid contribution 0.796 0.244 1.061 1.348 sum 1.746 -1.074 1.279 2.416 Atomic orbital electron populations 1.86446 1.17098 1.82344 1.51405 1.24689 0.85903 1.10454 0.96990 0.94586 0.77179 0.75901 0.84611 1.19641 1.01253 0.95011 0.96890 1.21057 0.96030 0.93340 0.96506 1.20726 1.01329 1.01467 0.93628 1.19105 0.94516 0.91077 0.86985 1.86049 1.68865 1.38482 1.29452 1.21072 0.97898 0.90759 0.84046 1.21911 0.90771 1.00396 0.88107 1.21344 0.96984 1.05303 0.98746 1.17327 0.81472 0.78901 0.85022 1.90842 1.80189 1.44980 1.27196 1.45219 1.16894 1.49167 1.21719 1.19898 0.86250 0.98487 0.84598 1.19774 0.94660 1.02847 0.94307 1.17471 0.79026 1.02447 0.94225 1.91355 1.35862 1.88411 1.95748 1.21390 0.98055 1.02656 0.95693 1.21312 0.92271 0.97093 1.00254 1.21420 1.00934 1.00667 0.92497 1.20844 0.93887 0.96184 0.97829 1.23438 0.94193 0.96953 1.01261 1.22603 0.97539 0.94800 0.81547 0.89015 1.86341 1.75686 1.49757 1.24909 1.24454 0.98259 0.97129 1.02592 1.23225 1.00702 1.04466 0.93166 1.24912 0.92609 0.92732 0.97188 1.83795 1.44519 1.53836 1.25094 1.21286 1.06528 0.94086 1.00812 1.20832 0.98198 0.98876 0.93826 0.74544 0.88587 1.22528 1.22104 0.82794 0.83634 0.88854 0.87804 0.74491 0.82854 0.82590 0.84598 0.87002 0.88619 0.86617 0.75146 0.87167 0.92067 0.79516 0.73070 0.84303 0.86314 Number of geometries 1 Number of calculations of the screened coulomb radii 1 The total number of SCF iterations 15. 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 O -0.56 -0.52 12.81 -148.98 -1.91 -2.43 16 2 C -0.07 -0.09 2.79 71.24 0.20 0.11 16 3 Si 0.77 0.24 29.85 68.60 2.05 2.29 16 4 C -0.13 -0.97 5.39 -19.85 -0.11 -1.08 16 5 C -0.05 -0.33 9.68 22.25 0.22 -0.11 16 6 C -0.15 -1.23 9.99 22.36 0.22 -1.01 16 7 C 0.13 1.36 6.69 22.50 0.15 1.51 16 8 O -0.31 -3.52 10.12 -92.46 -0.94 -4.46 16 9 C 0.14 1.07 4.74 71.24 0.34 1.41 16 10 C 0.04 0.30 7.00 25.26 0.18 0.47 16 11 C -0.22 -1.26 6.31 -18.99 -0.12 -1.38 16 12 C 0.58 4.74 5.75 86.65 0.50 5.24 16 13 O -0.55 -7.81 12.23 -4.10 -0.05 -7.86 16 14 N -0.68 -2.36 4.36 -427.19 -1.86 -4.22 16 15 C 0.20 0.84 1.03 3.38 0.00 0.84 16 16 C -0.11 -0.88 4.57 -19.82 -0.09 -0.97 16 17 C 0.09 0.61 6.71 22.38 0.15 0.76 16 18 F -0.14 -0.78 12.93 44.97 0.58 -0.20 16 19 C -0.16 -0.70 10.03 22.31 0.22 -0.48 16 20 C -0.09 -0.32 10.04 22.26 0.22 -0.09 16 21 C -0.14 -0.85 10.04 22.27 0.22 -0.63 16 22 C -0.07 -0.61 7.48 22.26 0.17 -0.44 16 23 C -0.12 0.01 6.88 31.12 0.21 0.22 16 24 C 0.09 -0.10 5.28 31.18 0.16 0.07 16 25 H 0.09 -0.35 8.14 -2.38 -0.02 -0.37 16 26 O -0.56 1.08 12.96 -148.98 -1.93 -0.85 16 27 C -0.19 -0.75 6.92 31.18 0.22 -0.53 16 28 C -0.20 0.24 10.60 23.65 0.25 0.49 16 29 C -0.01 0.02 12.03 66.89 0.80 0.83 16 30 O -0.18 -0.79 10.52 4.51 0.05 -0.75 16 31 C -0.21 -2.04 9.04 22.37 0.20 -1.84 16 32 C -0.10 -0.96 8.66 22.26 0.19 -0.77 16 33 H 0.40 -3.99 8.74 -74.06 -0.65 -4.64 16 34 H 0.10 -0.31 7.96 -2.38 -0.02 -0.32 16 35 H -0.30 -2.31 7.05 99.48 0.70 -1.61 16 36 H -0.29 -1.91 7.11 99.48 0.71 -1.20 16 37 H 0.15 0.28 8.06 -2.91 -0.02 0.25 16 38 H 0.15 0.71 8.06 -2.91 -0.02 0.69 16 39 H 0.09 0.39 7.65 -2.39 -0.02 0.37 16 40 H 0.10 0.99 6.02 -2.39 -0.01 0.98 16 41 H 0.42 -0.21 7.65 -92.71 -0.71 -0.92 16 42 H 0.15 0.00 8.06 -2.91 -0.02 -0.02 16 43 H 0.16 -0.33 8.06 -2.91 -0.02 -0.36 16 44 H 0.14 0.45 8.06 -2.91 -0.02 0.43 16 45 H 0.11 1.14 5.65 -2.91 -0.02 1.12 16 46 H 0.10 -0.16 8.11 -2.38 -0.02 -0.18 16 47 H 0.12 -0.07 6.94 -2.39 -0.02 -0.08 16 48 H 0.40 -3.32 9.12 -74.06 -0.68 -3.99 16 49 H 0.11 0.57 5.22 -2.39 -0.01 0.55 16 50 H 0.06 0.45 8.14 -2.38 -0.02 0.43 16 51 H 0.19 -1.35 7.74 -2.91 -0.02 -1.37 16 52 H 0.25 -2.97 8.06 -2.91 -0.02 -3.00 16 53 H 0.14 0.91 6.30 -2.91 -0.02 0.89 16 54 H 0.12 1.02 4.86 -2.91 -0.01 1.01 16 Total: 0.00 -26.75 444.17 -0.47 -27.21 By element: Atomic # 1 Polarization: -10.36 SS G_CDS: -0.97 Total: -11.33 kcal Atomic # 6 Polarization: -1.91 SS G_CDS: 4.52 Total: 2.61 kcal Atomic # 7 Polarization: -2.36 SS G_CDS: -1.86 Total: -4.22 kcal Atomic # 8 Polarization: -11.57 SS G_CDS: -4.78 Total: -16.35 kcal Atomic # 9 Polarization: -0.78 SS G_CDS: 0.58 Total: -0.20 kcal Atomic # 14 Polarization: 0.24 SS G_CDS: 2.05 Total: 2.29 kcal Total: -26.75 -0.47 -27.21 kcal The number of atoms in the molecule is 54 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. s_22_17639084_4549634.mol2 54 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 -35.431 kcal (2) G-P(sol) polarization free energy of solvation -26.746 kcal (3) G-ENP(sol) elect.-nuc.-pol. free energy of system -62.177 kcal (4) G-CDS(sol) cavity-dispersion-solvent structure free energy -0.467 kcal (5) G-P-CDS(sol) = G-P(sol) + G-CDS(sol) = (2) + (4) -27.213 kcal (6) G-S(sol) free energy of system = (1) + (5) -62.644 kcal 1SCF WAS SPECIFIED, SO GEOMETRY OPTIMIZATION WAS NOT USED Total computer time = 1.93 seconds