Wall clock time and date at job start Wed Apr 1 2020 00:58:40 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 O 1.45201 * 1 3 3 C 1.34231 * 117.00561 * 2 1 4 4 O 1.20829 * 119.99529 * 0.02562 * 3 2 1 5 5 C 1.50699 * 120.00129 * 180.02562 * 3 2 1 6 6 C 1.52999 * 109.46732 * 180.02562 * 5 3 2 7 7 N 1.46908 * 109.46936 * 179.97438 * 6 5 3 8 8 C 1.46894 * 110.99932 * 293.95552 * 7 6 5 9 9 C 1.53235 * 109.34173 * 187.25618 * 8 7 6 10 10 N 1.47081 * 108.51816 * 306.16751 * 9 8 7 11 11 C 1.34777 * 120.88228 * 231.16246 * 10 9 8 12 12 O 1.21286 * 119.99791 * 180.02562 * 11 10 9 13 13 C 1.50693 * 119.99927 * 0.02562 * 11 10 9 14 14 C 1.53002 * 109.46950 * 180.02562 * 13 11 10 15 15 C 1.53004 * 109.47341 * 180.02562 * 14 13 11 16 16 C 1.52989 * 109.47531 * 179.97438 * 15 14 13 17 17 C 1.50703 * 109.47449 * 180.02562 * 16 15 14 18 18 H 1.08002 * 119.99770 * 359.97438 * 17 16 15 19 19 C 1.37877 * 120.00012 * 180.02562 * 17 16 15 20 20 N 1.31515 * 119.99719 * 179.97438 * 19 17 16 21 21 Si 1.86791 * 120.00251 * 359.97438 * 19 17 16 22 22 H 1.48500 * 109.47130 * 90.00178 * 21 19 17 23 23 H 1.48501 * 109.47254 * 210.00268 * 21 19 17 24 24 H 1.48504 * 109.47168 * 330.00114 * 21 19 17 25 25 C 1.47086 * 118.23796 * 51.18621 * 10 9 8 26 26 C 1.46894 * 110.99613 * 58.25684 * 7 6 5 27 27 H 1.08997 * 109.49367 * 292.70105 * 26 7 6 28 28 C 1.52998 * 109.49654 * 52.78750 * 26 7 6 29 29 H 1.09003 * 109.46842 * 59.99706 * 1 2 3 30 30 H 1.08993 * 109.47249 * 179.97438 * 1 2 3 31 31 H 1.09005 * 109.46729 * 300.00031 * 1 2 3 32 32 H 1.09001 * 109.47038 * 60.00544 * 5 3 2 33 33 H 1.08996 * 109.46855 * 300.00088 * 5 3 2 34 34 H 1.08997 * 109.47149 * 59.99949 * 6 5 3 35 35 H 1.08996 * 109.47443 * 299.99631 * 6 5 3 36 36 H 1.09004 * 109.49408 * 67.28971 * 8 7 6 37 37 H 1.09005 * 109.49423 * 307.21355 * 8 7 6 38 38 H 1.08998 * 109.63341 * 186.46272 * 9 8 7 39 39 H 1.08997 * 109.63256 * 65.87637 * 9 8 7 40 40 H 1.08995 * 109.47366 * 299.99609 * 13 11 10 41 41 H 1.09003 * 109.47058 * 60.00251 * 13 11 10 42 42 H 1.09001 * 109.47231 * 299.99464 * 14 13 11 43 43 H 1.09002 * 109.47510 * 59.99904 * 14 13 11 44 44 H 1.09001 * 109.46965 * 300.00323 * 15 14 13 45 45 H 1.09002 * 109.47025 * 59.99450 * 15 14 13 46 46 H 1.09001 * 109.47313 * 300.00175 * 16 15 14 47 47 H 1.09002 * 109.47362 * 59.99722 * 16 15 14 48 48 H 0.97002 * 119.99875 * 180.02562 * 20 19 17 49 49 H 1.08995 * 109.62590 * 68.52272 * 25 10 9 50 50 H 1.09001 * 109.62973 * 189.10591 * 25 10 9 51 51 H 1.08999 * 109.47174 * 179.86114 * 28 26 7 52 52 H 1.08996 * 109.47207 * 299.86091 * 28 26 7 53 53 H 1.09006 * 109.47059 * 59.86547 * 28 26 7 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 8 1.4520 0.0000 0.0000 3 6 2.0615 1.1959 0.0000 4 8 1.4035 2.2093 0.0005 5 6 3.5664 1.2747 -0.0006 6 6 4.0003 2.7419 0.0002 7 7 5.4673 2.8187 0.0002 8 6 6.0221 2.2895 1.2532 9 6 7.5314 2.5520 1.2908 10 7 8.1314 2.0043 0.0647 11 6 9.1950 1.1786 0.1242 12 8 9.6774 0.7366 -0.8970 13 6 9.7847 0.8046 1.4596 14 6 10.9732 -0.1352 1.2470 15 6 11.5716 -0.5156 2.6028 16 6 12.7596 -1.4558 2.3903 17 6 13.3486 -1.8309 3.7258 18 1 12.9188 -1.4376 4.6352 19 6 14.4368 -2.6754 3.7866 20 7 14.9511 -3.0023 4.9520 21 14 15.1806 -3.3552 2.2138 22 1 16.2374 -2.4325 1.7269 23 1 15.7741 -4.6895 2.4832 24 1 14.1216 -3.4839 1.1807 25 6 7.5404 2.3821 -1.2281 26 6 6.0307 2.1263 -1.1664 27 1 5.8465 1.0557 -1.0779 28 6 5.3708 2.6548 -2.4416 29 1 -0.3633 0.5139 -0.8900 30 1 -0.3633 -1.0276 -0.0005 31 1 -0.3633 0.5139 0.8900 32 1 3.9558 0.7808 -0.8909 33 1 3.9564 0.7804 0.8891 34 1 3.6109 3.2354 0.8906 35 1 3.6103 3.2366 -0.8893 36 1 5.8380 1.2166 1.3080 37 1 5.5457 2.7844 2.0996 38 1 7.9667 2.0622 2.1618 39 1 7.7153 3.6252 1.3407 40 1 10.1215 1.7051 1.9731 41 1 9.0281 0.3034 2.0633 42 1 10.6368 -1.0359 0.7333 43 1 11.7296 0.3663 0.6434 44 1 11.9078 0.3848 3.1170 45 1 10.8149 -1.0172 3.2061 46 1 12.4230 -2.3562 1.8763 47 1 13.5164 -0.9546 1.7869 48 1 15.7169 -3.5961 4.9948 49 1 7.7252 3.4389 -1.4203 50 1 7.9816 1.7799 -2.0224 51 1 5.7909 2.1425 -3.3072 52 1 4.2971 2.4727 -2.3974 53 1 5.5545 3.7257 -2.5288 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 ZINC000589844898.mol2 53 1SCF WAS SPECIFIED, SO GEOMETRY OPTIMIZATION WAS NOT USED SCF FIELD WAS ACHIEVED Wall clock time and date at job start Wed Apr 1 2020 00:58:40 Heat of formation + Delta-G solvation = -128.815221 kcal Electronic energy + Delta-G solvation = -29834.471680 eV Core-core repulsion = 25705.923574 eV Total energy + Delta-G solvation = -4128.548106 eV No. of doubly occupied orbitals = 66 Molecular weight (most abundant/longest-lived isotopes) = 340.220 amu Computer time = 0.98 seconds Orbital eigenvalues (eV) -41.17397 -40.50293 -38.89336 -37.05711 -35.44682 -35.05449 -33.30254 -32.50834 -31.68304 -30.22952 -28.74211 -27.83424 -27.17854 -25.50153 -23.89696 -23.31543 -22.06228 -21.59908 -20.07119 -19.96707 -18.87003 -18.30229 -17.29562 -17.20475 -17.05568 -16.55003 -16.45162 -15.98694 -15.87435 -15.28399 -14.90218 -14.78490 -14.46193 -14.19632 -13.90620 -13.82669 -13.43518 -13.38792 -13.04337 -12.73542 -12.64434 -12.56814 -12.49167 -12.39487 -12.23998 -12.01447 -11.95159 -11.76294 -11.53760 -11.39344 -11.34031 -11.17986 -11.10364 -11.01433 -10.71934 -10.53844 -9.92984 -9.80570 -9.64481 -9.41555 -9.12599 -9.11462 -8.69517 -8.49009 -6.00280 -3.97677 1.65846 2.16800 2.61268 2.74714 3.30823 3.35580 3.42561 3.48575 4.16617 4.21386 4.28214 4.37973 4.44142 4.50396 4.54794 4.59489 4.66956 4.68516 4.72910 4.81141 4.89494 4.91999 5.02724 5.15298 5.16297 5.22398 5.24928 5.28257 5.34644 5.38033 5.41430 5.49609 5.51977 5.61580 5.76887 5.79604 5.85162 6.02469 6.12213 6.23641 6.41333 6.44859 6.51148 6.68280 6.81945 6.99335 7.19378 7.24092 7.38555 7.51311 7.85129 8.13383 8.58582 9.04667 9.63170 11.13185 Molecular weight = 340.22amu Principal moments of inertia in cm(-1) A = 0.015680 B = 0.001822 C = 0.001731 Principal moments of inertia in units of 10**(-40)*gram-cm**2 A = 1785.229848 B =15363.863330 C =16169.814511 Net atomic charges, atomic populations, and dipole contributions calculated with CM2 Atom NO. Type Charge No. of electrons 1 C 0.035 3.965 2 O -0.361 6.361 3 C 0.466 3.534 4 O -0.506 6.506 5 C -0.158 4.158 6 C 0.066 3.934 7 N -0.525 5.525 8 C 0.033 3.967 9 C 0.111 3.889 10 N -0.623 5.623 11 C 0.516 3.484 12 O -0.537 6.537 13 C -0.138 4.138 14 C -0.106 4.106 15 C -0.102 4.102 16 C -0.016 4.016 17 C -0.519 4.519 18 H 0.071 0.929 19 C 0.064 3.936 20 N -0.906 5.906 21 Si 0.880 3.120 22 H -0.277 1.277 23 H -0.289 1.289 24 H -0.265 1.265 25 C 0.121 3.879 26 C 0.062 3.938 27 H 0.051 0.949 28 C -0.150 4.150 29 H 0.063 0.937 30 H 0.102 0.898 31 H 0.064 0.936 32 H 0.108 0.892 33 H 0.107 0.893 34 H 0.088 0.912 35 H 0.093 0.907 36 H 0.044 0.956 37 H 0.092 0.908 38 H 0.095 0.905 39 H 0.078 0.922 40 H 0.092 0.908 41 H 0.091 0.909 42 H 0.063 0.937 43 H 0.063 0.937 44 H 0.053 0.947 45 H 0.052 0.948 46 H 0.015 0.985 47 H 0.011 0.989 48 H 0.262 0.738 49 H 0.073 0.927 50 H 0.095 0.905 51 H 0.066 0.934 52 H 0.060 0.940 53 H 0.070 0.930 For the charges calculated by CM2: Dipole moment (debyes) X Y Z Total from point charges -30.125 13.111 -12.092 35.009 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.059 4.059 2 O -0.278 6.278 3 C 0.308 3.692 4 O -0.394 6.394 5 C -0.198 4.198 6 C -0.062 4.062 7 N -0.248 5.248 8 C -0.096 4.096 9 C -0.012 4.012 10 N -0.358 5.358 11 C 0.305 3.695 12 O -0.413 6.413 13 C -0.177 4.177 14 C -0.144 4.144 15 C -0.140 4.140 16 C -0.053 4.053 17 C -0.557 4.557 18 H 0.089 0.911 19 C -0.136 4.136 20 N -0.545 5.545 21 Si 0.707 3.293 22 H -0.203 1.203 23 H -0.215 1.215 24 H -0.189 1.189 25 C -0.002 4.002 26 C -0.047 4.047 27 H 0.069 0.931 28 C -0.208 4.208 29 H 0.082 0.918 30 H 0.121 0.879 31 H 0.082 0.918 32 H 0.126 0.874 33 H 0.125 0.875 34 H 0.106 0.894 35 H 0.111 0.889 36 H 0.062 0.938 37 H 0.110 0.890 38 H 0.114 0.886 39 H 0.097 0.903 40 H 0.111 0.889 41 H 0.109 0.891 42 H 0.081 0.919 43 H 0.082 0.918 44 H 0.072 0.928 45 H 0.071 0.929 46 H 0.034 0.966 47 H 0.030 0.970 48 H 0.072 0.928 49 H 0.092 0.908 50 H 0.114 0.886 51 H 0.085 0.915 52 H 0.079 0.921 53 H 0.089 0.911 Dipole moment (debyes) X Y Z Total from point charges -29.779 13.549 -12.591 35.055 hybrid contribution -0.095 -1.249 -0.528 1.359 sum -29.874 12.301 -13.118 34.869 Atomic orbital electron populations 1.23298 0.76018 1.03915 1.02638 1.86530 1.23740 1.33531 1.83956 1.23659 0.91919 0.81389 0.72185 1.90688 1.67249 1.36247 1.45184 1.22261 0.88266 1.02932 1.06292 1.22166 0.87104 0.91906 1.05049 1.62281 1.04683 1.55628 1.02221 1.22609 0.95539 0.98782 0.92628 1.21853 0.91819 1.00121 0.87417 1.48264 1.31780 1.48510 1.07231 1.21101 0.78942 0.77051 0.92391 1.90692 1.60282 1.53344 1.36955 1.21799 0.99988 1.02073 0.93883 1.21493 0.95323 0.98274 0.99327 1.21474 0.98553 0.99774 0.94190 1.19575 0.93503 0.94650 0.97568 1.21185 1.13714 1.28286 0.92544 0.91107 1.25043 0.96617 0.95995 0.95972 1.69922 1.24239 1.36355 1.23995 0.86774 0.79746 0.78378 0.84401 1.20317 1.21464 1.18916 1.21682 0.90684 1.01286 0.86517 1.22105 0.94991 0.96597 0.91029 0.93121 1.21907 1.00453 1.02577 0.95845 0.91803 0.87943 0.91754 0.87390 0.87463 0.89375 0.88900 0.93771 0.88970 0.88650 0.90349 0.88945 0.89102 0.91858 0.91778 0.92844 0.92879 0.96643 0.97049 0.92829 0.90843 0.88629 0.91520 0.92114 0.91142 Number of geometries 1 Number of calculations of the screened coulomb radii 1 The total number of SCF iterations 13. 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.03 0.18 10.56 101.05 1.07 1.24 16 2 O -0.36 -2.71 10.56 -39.26 -0.41 -3.13 16 3 C 0.47 3.60 7.96 36.00 0.29 3.88 16 4 O -0.51 -5.09 15.73 -18.75 -0.29 -5.39 16 5 C -0.16 -0.89 4.98 -27.88 -0.14 -1.03 16 6 C 0.07 0.39 4.08 -3.82 -0.02 0.37 16 7 N -0.52 -3.35 4.68 -226.35 -1.06 -4.41 16 8 C 0.03 0.19 5.16 -3.71 -0.02 0.18 16 9 C 0.11 0.78 6.50 -3.60 -0.02 0.76 16 10 N -0.62 -6.23 3.01 -172.46 -0.52 -6.75 16 11 C 0.52 6.87 7.71 -10.99 -0.08 6.78 16 12 O -0.54 -9.29 15.44 5.55 0.09 -9.21 16 13 C -0.14 -1.71 4.07 -27.89 -0.11 -1.82 16 14 C -0.11 -1.77 4.50 -26.73 -0.12 -1.89 16 15 C -0.10 -1.99 5.15 -26.73 -0.14 -2.13 16 16 C -0.02 -0.38 4.04 -27.89 -0.11 -0.49 16 17 C -0.52 -15.08 9.66 -34.44 -0.33 -15.41 16 18 H 0.07 2.18 7.99 -52.49 -0.42 1.76 16 19 C 0.06 1.92 5.47 -17.82 -0.10 1.82 16 20 N -0.91 -28.34 13.96 55.43 0.77 -27.57 16 21 Si 0.88 22.37 27.47 -169.99 -4.67 17.70 16 22 H -0.28 -6.93 7.11 56.52 0.40 -6.53 16 23 H -0.29 -7.23 7.11 56.52 0.40 -6.83 16 24 H -0.26 -6.53 6.54 56.52 0.37 -6.16 16 25 C 0.12 1.10 5.88 -3.59 -0.02 1.07 16 26 C 0.06 0.43 1.97 -67.59 -0.13 0.30 16 27 H 0.05 0.39 6.70 -51.93 -0.35 0.04 16 28 C -0.15 -0.84 8.61 37.16 0.32 -0.52 16 29 H 0.06 0.31 8.13 -51.93 -0.42 -0.11 16 30 H 0.10 0.27 8.14 -51.93 -0.42 -0.15 16 31 H 0.06 0.31 8.13 -51.93 -0.42 -0.11 16 32 H 0.11 0.57 6.17 -51.93 -0.32 0.25 16 33 H 0.11 0.53 6.71 -51.93 -0.35 0.19 16 34 H 0.09 0.54 8.00 -51.93 -0.42 0.13 16 35 H 0.09 0.59 6.28 -51.93 -0.33 0.26 16 36 H 0.04 0.27 6.83 -51.93 -0.35 -0.08 16 37 H 0.09 0.42 8.04 -51.93 -0.42 0.00 16 38 H 0.10 0.58 6.02 -51.93 -0.31 0.27 16 39 H 0.08 0.52 8.14 -51.93 -0.42 0.09 16 40 H 0.09 0.99 8.00 -51.93 -0.42 0.58 16 41 H 0.09 0.96 7.66 -51.93 -0.40 0.56 16 42 H 0.06 1.12 8.10 -51.93 -0.42 0.70 16 43 H 0.06 1.13 8.10 -51.93 -0.42 0.71 16 44 H 0.05 1.01 8.14 -51.93 -0.42 0.59 16 45 H 0.05 1.00 8.14 -51.93 -0.42 0.58 16 46 H 0.02 0.38 7.20 -51.93 -0.37 0.01 16 47 H 0.01 0.28 7.57 -51.93 -0.39 -0.12 16 48 H 0.26 7.67 8.31 -40.82 -0.34 7.33 16 49 H 0.07 0.63 8.14 -51.93 -0.42 0.21 16 50 H 0.10 1.02 7.06 -51.93 -0.37 0.65 16 51 H 0.07 0.35 8.14 -51.93 -0.42 -0.07 16 52 H 0.06 0.32 5.96 -51.93 -0.31 0.01 16 53 H 0.07 0.40 8.14 -51.93 -0.42 -0.02 16 LS Contribution 411.86 15.07 6.21 6.21 Total: -1.00 -35.81 411.86 -8.90 -44.71 By element: Atomic # 1 Polarization: 4.06 SS G_CDS: -9.33 Total: -5.27 kcal Atomic # 6 Polarization: -7.22 SS G_CDS: 0.32 Total: -6.90 kcal Atomic # 7 Polarization: -37.92 SS G_CDS: -0.80 Total: -38.73 kcal Atomic # 8 Polarization: -17.10 SS G_CDS: -0.62 Total: -17.72 kcal Atomic # 14 Polarization: 22.37 SS G_CDS: -4.67 Total: 17.70 kcal Total LS contribution 6.21 Total: 6.21 kcal Total: -35.81 -8.90 -44.71 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. ZINC000589844898.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 -84.105 kcal (2) G-P(sol) polarization free energy of solvation -35.812 kcal (3) G-ENP(sol) elect.-nuc.-pol. free energy of system -119.917 kcal (4) G-CDS(sol) cavity-dispersion-solvent structure free energy -8.899 kcal (5) G-P-CDS(sol) = G-P(sol) + G-CDS(sol) = (2) + (4) -44.710 kcal (6) G-S(sol) free energy of system = (1) + (5) -128.815 kcal 1SCF WAS SPECIFIED, SO GEOMETRY OPTIMIZATION WAS NOT USED Total computer time = 0.98 seconds