Wall clock time and date at job start Tue Mar 9 2021 13:10:04 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.42903 * 1 3 3 Si 1.86297 * 109.47024 * 2 1 4 4 C 1.50699 * 109.47219 * 120.00381 * 2 1 3 5 5 C 1.38240 * 119.96941 * 245.02286 * 4 2 1 6 6 C 1.38289 * 120.11190 * 180.02562 * 5 4 2 7 7 C 1.38125 * 120.05452 * 359.67495 * 6 5 4 8 8 C 1.38755 * 119.94082 * 0.59884 * 7 6 5 9 9 O 1.35900 * 120.05506 * 179.70103 * 8 7 6 10 10 C 1.42902 * 117.00171 * 359.72458 * 9 8 7 11 11 C 1.50696 * 109.47173 * 180.02562 * 10 9 8 12 12 C 1.39555 * 120.07345 * 180.02562 * 11 10 9 13 13 C 1.47850 * 120.14235 * 359.74767 * 12 11 10 14 14 O 1.21550 * 120.00396 * 359.97438 * 13 12 11 15 15 N 1.34780 * 119.99976 * 180.02562 * 13 12 11 16 16 C 1.46498 * 120.00111 * 179.97438 * 15 13 12 17 17 H 1.09009 * 112.85041 * 23.79812 * 16 15 13 18 18 C 1.53778 * 113.61582 * 154.99145 * 16 15 13 19 19 C 1.53788 * 87.08089 * 139.97952 * 18 16 15 20 20 H 1.09005 * 113.61055 * 89.11576 * 19 18 16 21 21 C 1.52997 * 113.61898 * 220.02147 * 19 18 16 22 22 N 1.46496 * 109.47199 * 185.00167 * 21 19 18 23 23 C 1.34431 * 126.82246 * 124.71138 * 22 21 19 24 24 N 1.30838 * 108.07464 * 179.89246 * 23 22 21 25 25 C 1.33237 * 109.82741 * 0.37981 * 24 23 22 26 26 N 1.30716 * 108.94970 * 359.76388 * 25 24 23 27 27 C 1.53776 * 87.08043 * 334.57256 * 19 18 16 28 28 C 1.39593 * 119.71735 * 180.02562 * 12 11 10 29 29 C 1.37938 * 119.84737 * 359.97438 * 28 12 11 30 30 C 1.38377 * 120.14805 * 0.02562 * 29 28 12 31 31 C 1.37972 * 120.07587 * 0.02562 * 11 10 9 32 32 C 1.38160 * 119.97213 * 64.99942 * 4 2 1 33 33 H 0.96702 * 113.99733 * 59.99702 * 1 2 3 34 34 H 1.09000 * 109.46921 * 240.00382 * 2 1 3 35 35 H 1.48506 * 110.00028 * 60.00269 * 3 2 1 36 36 H 1.48501 * 110.00154 * 181.32033 * 3 2 1 37 37 H 1.07996 * 119.94878 * 359.97438 * 5 4 2 38 38 H 1.08004 * 119.97245 * 179.97438 * 6 5 4 39 39 H 1.08002 * 120.02949 * 180.32503 * 7 6 5 40 40 H 1.09005 * 109.46991 * 59.99897 * 10 9 8 41 41 H 1.09002 * 109.46767 * 300.00483 * 10 9 8 42 42 H 0.97000 * 120.00217 * 359.97438 * 15 13 12 43 43 H 1.08998 * 113.62086 * 254.52986 * 18 16 15 44 44 H 1.09004 * 113.61498 * 25.35251 * 18 16 15 45 45 H 1.09004 * 109.46859 * 304.99849 * 21 19 18 46 46 H 1.08999 * 109.47067 * 64.99979 * 21 19 18 47 47 H 1.07999 * 125.96113 * 0.02932 * 23 22 21 48 48 H 1.08005 * 125.52169 * 179.79939 * 25 24 23 49 49 H 1.09005 * 113.61165 * 270.88390 * 27 19 18 50 50 H 1.08998 * 113.61938 * 139.97936 * 27 19 18 51 51 H 1.07994 * 120.07972 * 179.97438 * 28 12 11 52 52 H 1.08006 * 119.92782 * 179.97438 * 29 28 12 53 53 H 1.08004 * 119.85611 * 179.97438 * 30 29 28 54 54 H 1.08004 * 119.92879 * 359.97438 * 31 11 10 55 55 H 1.07999 * 120.03138 * 0.04007 * 32 4 2 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.0500 1.7564 0.0000 4 6 1.9314 -0.7105 1.2304 5 6 2.6383 -0.0117 2.1911 6 6 3.1011 -0.6608 3.3210 7 6 2.8521 -2.0082 3.4956 8 6 2.1477 -2.7130 2.5301 9 8 1.9067 -4.0401 2.6966 10 6 2.4130 -4.6471 3.8871 11 6 2.0522 -6.1102 3.8975 12 6 2.4403 -6.9186 4.9668 13 6 3.2113 -6.3443 6.0900 14 8 3.5089 -5.1658 6.0833 15 7 3.5852 -7.1278 7.1210 16 6 4.3496 -6.5589 8.2338 17 1 4.2073 -5.4832 8.3389 18 6 4.1874 -7.3334 9.5523 19 6 5.6910 -7.1051 9.7801 20 1 5.9169 -6.1787 10.3082 21 6 6.4127 -8.3183 10.3703 22 7 7.8046 -7.9651 10.6601 23 6 8.4364 -8.0821 11.8409 24 7 9.6665 -7.6640 11.6864 25 6 9.8540 -7.2663 10.4286 26 7 8.7423 -7.4313 9.7611 27 6 5.8252 -6.9916 8.2524 28 6 2.1026 -8.2730 4.9685 29 6 1.3876 -8.8036 3.9150 30 6 1.0058 -7.9966 2.8578 31 6 1.3373 -6.6535 2.8499 32 6 1.6846 -2.0595 1.3975 33 1 -0.3933 0.4418 -0.7651 34 1 1.7923 -0.5138 -0.8900 35 1 1.5615 2.4679 1.2086 36 1 3.5346 1.7703 -0.0322 37 1 2.8303 1.0428 2.0586 38 1 3.6536 -0.1127 4.0700 39 1 3.2093 -2.5136 4.3807 40 1 1.9761 -4.1585 4.7581 41 1 3.4974 -4.5404 3.9164 42 1 3.3481 -8.0683 7.1263 43 1 3.5536 -6.8277 10.2808 44 1 3.9196 -8.3809 9.4133 45 1 5.9161 -8.6238 11.2912 46 1 6.3878 -9.1396 9.6542 47 1 8.0047 -8.4580 12.7567 48 1 10.7753 -6.8721 10.0258 49 1 6.0065 -7.9460 7.7579 50 1 6.5135 -6.2099 7.9314 51 1 2.4004 -8.9041 5.7928 52 1 1.1260 -9.8515 3.9148 53 1 0.4464 -8.4171 2.0351 54 1 1.0359 -6.0290 2.0220 55 1 1.1325 -2.6051 0.6466 There are 77 doubly filled levels RHF calculation, no. of doubly occupied orbitals= 77 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) 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_14644236_4557126.mol2 55 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:10:04 Heat of formation + Delta-G solvation = 83.421658 kcal Electronic energy + Delta-G solvation = -38708.400464 eV Core-core repulsion = 33680.669702 eV Total energy + Delta-G solvation = -5027.730761 eV No. of doubly occupied orbitals = 77 and no. of open levels = 1 Molecular weight (most abundant/longest-lived isotopes) = 421.183 amu Computer time = 1.92 seconds Orbital eigenvalues (eV) -42.98742 -41.90282 -41.31668 -40.33269 -39.23927 -38.87513 -37.35654 -35.62448 -34.27875 -33.38237 -32.75693 -32.18163 -32.11467 -31.18548 -29.09183 -27.96606 -26.76948 -26.70467 -25.42075 -24.26689 -23.78199 -23.15629 -22.86970 -22.30477 -21.89202 -21.57335 -20.75127 -19.00985 -18.93882 -18.42831 -18.34568 -17.67290 -17.60877 -17.40301 -17.27043 -16.80925 -16.54578 -16.31665 -16.14708 -16.07773 -15.92011 -15.65220 -15.50697 -15.38227 -15.06506 -14.72485 -14.68871 -14.59951 -14.53778 -14.31394 -14.12390 -13.70471 -13.66966 -13.52230 -13.47284 -13.27811 -13.18348 -12.91238 -12.80379 -12.48892 -12.35922 -12.25666 -12.24498 -12.17119 -11.94694 -11.89036 -11.43797 -11.00485 -10.94650 -10.80055 -10.57828 -10.42199 -10.29926 -10.07223 -10.00108 -9.85999 -9.46226 -4.96631 -0.30469 0.00329 0.16340 0.23111 0.69941 1.25402 1.45418 1.63395 1.70885 1.87058 1.96923 1.99610 2.01447 2.44068 2.69327 3.28296 3.32646 3.40162 3.46100 3.51554 3.53445 3.64480 3.65444 3.71445 3.75779 3.78211 3.82805 3.91770 4.00312 4.05400 4.09996 4.11115 4.15475 4.16818 4.24018 4.30614 4.33124 4.37677 4.39696 4.55126 4.55710 4.65288 4.67904 4.73416 4.77099 4.78102 4.87219 4.94620 5.01452 5.03176 5.11622 5.12767 5.18334 5.22162 5.29772 5.32980 5.40658 5.44761 5.48607 5.52120 5.79556 5.95170 6.17034 6.38840 6.88722 7.04465 7.13121 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.083 4.083 3 Si 0.786 3.214 4 C -0.050 4.050 5 C -0.178 4.178 6 C -0.076 4.076 7 C -0.211 4.211 8 C 0.123 3.877 9 O -0.321 6.321 10 C 0.117 3.883 11 C -0.056 4.056 12 C -0.148 4.148 13 C 0.553 3.447 14 O -0.576 6.576 15 N -0.697 5.697 16 C 0.114 3.886 17 H 0.090 0.910 18 C -0.112 4.112 19 C -0.111 4.111 20 H 0.083 0.917 21 C 0.135 3.865 22 N -0.379 5.379 23 C 0.256 3.744 24 N -0.586 5.586 25 C 0.183 3.817 26 N -0.345 5.345 27 C -0.135 4.135 28 C -0.075 4.075 29 C -0.122 4.122 30 C -0.089 4.089 31 C -0.123 4.123 32 C -0.134 4.134 33 H 0.405 0.595 34 H 0.097 0.903 35 H -0.296 1.296 36 H -0.289 1.289 37 H 0.124 0.876 38 H 0.148 0.852 39 H 0.136 0.864 40 H 0.075 0.925 41 H 0.075 0.925 42 H 0.417 0.583 43 H 0.090 0.910 44 H 0.123 0.877 45 H 0.145 0.855 46 H 0.101 0.899 47 H 0.260 0.740 48 H 0.228 0.772 49 H 0.093 0.907 50 H 0.042 0.958 51 H 0.164 0.836 52 H 0.164 0.836 53 H 0.153 0.847 54 H 0.133 0.867 55 H 0.141 0.859 For the charges calculated by CM2: Dipole moment (debyes) X Y Z Total from point charges -5.025 -7.319 0.364 8.886 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.192 4.192 3 Si 0.692 3.308 4 C -0.051 4.051 5 C -0.196 4.196 6 C -0.094 4.094 7 C -0.229 4.229 8 C 0.078 3.922 9 O -0.236 6.236 10 C 0.042 3.958 11 C -0.057 4.057 12 C -0.152 4.152 13 C 0.342 3.658 14 O -0.457 6.457 15 N -0.348 5.348 16 C 0.009 3.991 17 H 0.108 0.892 18 C -0.151 4.151 19 C -0.131 4.131 20 H 0.101 0.899 21 C 0.014 3.986 22 N -0.169 5.169 23 C -0.040 4.040 24 N -0.314 5.314 25 C -0.120 4.120 26 N -0.173 5.173 27 C -0.174 4.174 28 C -0.093 4.093 29 C -0.139 4.139 30 C -0.107 4.107 31 C -0.141 4.141 32 C -0.152 4.152 33 H 0.256 0.744 34 H 0.115 0.885 35 H -0.224 1.224 36 H -0.216 1.216 37 H 0.141 0.859 38 H 0.165 0.835 39 H 0.154 0.846 40 H 0.093 0.907 41 H 0.092 0.908 42 H 0.256 0.744 43 H 0.109 0.891 44 H 0.141 0.859 45 H 0.163 0.837 46 H 0.119 0.881 47 H 0.276 0.724 48 H 0.244 0.756 49 H 0.112 0.888 50 H 0.061 0.939 51 H 0.181 0.819 52 H 0.182 0.818 53 H 0.171 0.829 54 H 0.150 0.850 55 H 0.159 0.841 Dipole moment (debyes) X Y Z Total from point charges -4.287 -6.228 0.691 7.593 hybrid contribution 0.411 -0.383 0.513 0.761 sum -3.876 -6.611 1.204 7.758 Atomic orbital electron populations 1.86453 1.16942 1.82486 1.51440 1.24927 0.86128 1.11451 0.96689 0.94436 0.77180 0.74797 0.84417 1.19530 0.96662 0.93129 0.95743 1.21004 1.03198 0.99488 0.95907 1.21135 0.97409 0.93755 0.97100 1.21514 1.04479 0.96154 1.00788 1.19151 0.95240 0.82974 0.94841 1.85935 1.74584 1.19988 1.43052 1.21520 0.98757 0.90451 0.85024 1.19695 0.97826 0.93217 0.94922 1.19679 1.04448 0.95682 0.95355 1.18573 0.78642 0.85759 0.82778 1.90742 1.57612 1.17701 1.79628 1.45728 1.55556 1.13519 1.20037 1.22543 0.91145 0.98975 0.86449 0.89154 1.23308 0.94199 1.02768 0.94797 1.22491 0.96449 0.98480 0.95693 0.89881 1.22229 0.73969 0.97154 1.05225 1.50777 1.12064 1.46617 1.07477 1.26926 0.88402 0.95767 0.92874 1.70541 1.18451 1.22273 1.20136 1.26420 0.98413 0.98689 0.88484 1.78013 1.00415 1.21170 1.17719 1.23435 0.95740 1.01996 0.96260 1.21599 0.95899 0.92207 0.99630 1.21601 0.97427 1.01849 0.93044 1.21784 0.97037 0.92670 0.99174 1.21125 0.98256 0.96534 0.98158 1.20377 1.03333 0.94010 0.97519 0.74364 0.88486 1.22402 1.21614 0.85901 0.83473 0.84603 0.90722 0.90760 0.74402 0.89102 0.85927 0.83689 0.88056 0.72423 0.75587 0.88810 0.93872 0.81901 0.81810 0.82910 0.84955 0.84072 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 -1.37 12.81 -148.98 -1.91 -3.28 16 2 C -0.08 -0.21 2.79 71.24 0.20 -0.01 16 3 Si 0.79 0.11 29.85 68.60 2.05 2.16 16 4 C -0.05 -0.48 5.38 -19.88 -0.11 -0.59 16 5 C -0.18 -1.80 8.67 22.28 0.19 -1.60 16 6 C -0.08 -0.73 10.04 22.26 0.22 -0.50 16 7 C -0.21 -2.38 9.04 22.37 0.20 -2.18 16 8 C 0.12 1.51 6.69 22.51 0.15 1.66 16 9 O -0.32 -3.71 8.62 -92.43 -0.80 -4.51 16 10 C 0.12 1.19 3.38 71.23 0.24 1.44 16 11 C -0.06 -0.38 5.91 -19.61 -0.12 -0.49 16 12 C -0.15 -0.44 5.87 -20.09 -0.12 -0.56 16 13 C 0.55 2.69 7.75 86.79 0.67 3.36 16 14 O -0.58 -7.20 13.01 -3.87 -0.05 -7.25 16 15 N -0.70 2.13 5.31 -442.90 -2.35 -0.22 16 16 C 0.11 -0.29 4.24 45.52 0.19 -0.10 16 17 H 0.09 0.31 7.46 -2.38 -0.02 0.30 16 18 C -0.11 1.18 7.39 31.12 0.23 1.41 16 19 C -0.11 0.73 3.85 -10.27 -0.04 0.69 16 20 H 0.08 -0.21 8.14 -2.38 -0.02 -0.23 16 21 C 0.13 -1.28 6.39 86.38 0.55 -0.73 16 22 N -0.38 0.41 3.69 -352.07 -1.30 -0.89 16 23 C 0.26 -0.42 12.63 180.36 2.28 1.86 16 24 N -0.59 -3.50 11.19 -203.99 -2.28 -5.78 16 25 C 0.18 1.06 13.03 180.77 2.36 3.42 16 26 N -0.34 -2.10 12.57 -61.24 -0.77 -2.87 16 27 C -0.14 0.22 7.43 31.12 0.23 0.45 16 28 C -0.08 0.33 9.54 22.53 0.21 0.54 16 29 C -0.12 0.63 10.03 22.23 0.22 0.85 16 30 C -0.09 0.16 10.04 22.33 0.22 0.39 16 31 C -0.12 -0.48 9.37 22.23 0.21 -0.27 16 32 C -0.13 -1.44 9.62 22.37 0.22 -1.23 16 33 H 0.40 -3.61 8.74 -74.06 -0.65 -4.25 16 34 H 0.10 -0.18 7.96 -2.39 -0.02 -0.20 16 35 H -0.30 -2.25 7.05 99.48 0.70 -1.55 16 36 H -0.29 -1.63 7.11 99.48 0.71 -0.92 16 37 H 0.12 1.02 4.86 -2.91 -0.01 1.01 16 38 H 0.15 0.71 8.06 -2.91 -0.02 0.68 16 39 H 0.14 1.46 6.30 -2.91 -0.02 1.44 16 40 H 0.07 0.91 6.86 -2.38 -0.02 0.89 16 41 H 0.07 0.92 6.86 -2.39 -0.02 0.90 16 42 H 0.42 -4.30 6.76 -92.71 -0.63 -4.92 16 43 H 0.09 -0.98 8.14 -2.39 -0.02 -1.00 16 44 H 0.12 -1.97 8.08 -2.38 -0.02 -1.99 16 45 H 0.15 -2.49 8.09 -2.38 -0.02 -2.51 16 46 H 0.10 -1.10 7.84 -2.39 -0.02 -1.12 16 47 H 0.26 -3.17 8.06 -2.91 -0.02 -3.20 16 48 H 0.23 0.15 8.06 -2.91 -0.02 0.13 16 49 H 0.09 -0.29 7.90 -2.38 -0.02 -0.31 16 50 H 0.04 0.24 8.14 -2.39 -0.02 0.22 16 51 H 0.16 -2.00 6.41 -2.91 -0.02 -2.02 16 52 H 0.16 -1.86 8.06 -2.91 -0.02 -1.88 16 53 H 0.15 -1.00 8.06 -2.91 -0.02 -1.03 16 54 H 0.13 0.62 7.00 -2.91 -0.02 0.60 16 55 H 0.14 1.19 8.06 -2.91 -0.02 1.17 16 Total: 0.00 -35.36 454.18 0.73 -34.63 By element: Atomic # 1 Polarization: -19.52 SS G_CDS: -0.28 Total: -19.80 kcal Atomic # 6 Polarization: -0.62 SS G_CDS: 8.43 Total: 7.81 kcal Atomic # 7 Polarization: -3.06 SS G_CDS: -6.70 Total: -9.76 kcal Atomic # 8 Polarization: -12.28 SS G_CDS: -2.76 Total: -15.04 kcal Atomic # 14 Polarization: 0.11 SS G_CDS: 2.05 Total: 2.16 kcal Total: -35.36 0.73 -34.63 kcal The number of atoms in the molecule is 55 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_14644236_4557126.mol2 55 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 118.049 kcal (2) G-P(sol) polarization free energy of solvation -35.361 kcal (3) G-ENP(sol) elect.-nuc.-pol. free energy of system 82.688 kcal (4) G-CDS(sol) cavity-dispersion-solvent structure free energy 0.734 kcal (5) G-P-CDS(sol) = G-P(sol) + G-CDS(sol) = (2) + (4) -34.628 kcal (6) G-S(sol) free energy of system = (1) + (5) 83.422 kcal 1SCF WAS SPECIFIED, SO GEOMETRY OPTIMIZATION WAS NOT USED Total computer time = 1.92 seconds