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 O 1.42898 * 1 3 3 C 1.42895 * 114.00067 * 2 1 4 4 C 1.53002 * 109.54505 * 60.47316 * 3 2 1 5 5 H 1.08993 * 109.47301 * 169.39671 * 4 3 2 6 6 C 1.53002 * 109.46970 * 49.39833 * 4 3 2 7 7 N 1.46501 * 109.46862 * 289.40179 * 4 3 2 8 8 C 1.34776 * 119.99534 * 199.94706 * 7 4 3 9 9 O 1.21588 * 120.00788 * 4.60892 * 8 7 4 10 10 N 1.34776 * 119.99850 * 184.61414 * 8 7 4 11 11 C 1.46925 * 120.62872 * 355.08840 * 10 8 7 12 12 C 1.53193 * 108.77674 * 126.39720 * 11 10 8 13 13 C 1.53034 * 109.31589 * 54.63201 * 12 11 10 14 14 C 1.53042 * 109.54012 * 298.63800 * 13 12 11 15 15 H 1.08999 * 109.49751 * 301.38673 * 14 13 12 16 16 C 1.50701 * 109.50014 * 181.31484 * 14 13 12 17 17 H 1.07995 * 119.99831 * 0.06114 * 16 14 13 18 18 C 1.37874 * 119.99917 * 180.34475 * 16 14 13 19 19 N 1.31516 * 120.00254 * 359.71806 * 18 16 14 20 20 Si 1.86795 * 119.99952 * 179.72508 * 18 16 14 21 21 H 1.48502 * 109.47308 * 59.99658 * 20 18 16 22 22 H 1.48506 * 109.47365 * 179.97438 * 20 18 16 23 23 H 1.48500 * 109.47310 * 299.99712 * 20 18 16 24 24 C 1.46937 * 120.60708 * 175.07062 * 10 8 7 25 25 C 1.53179 * 109.54312 * 180.64836 * 3 2 1 26 26 C 1.53106 * 109.15525 * 63.06848 * 25 3 2 27 27 O 1.42894 * 109.49149 * 57.47137 * 26 25 3 28 28 C 1.42904 * 114.09783 * 298.89196 * 27 26 25 29 29 C 1.53090 * 109.41681 * 61.15647 * 28 27 26 30 30 H 1.09000 * 109.47108 * 170.23714 * 1 2 3 31 31 H 1.08994 * 109.47529 * 290.23611 * 1 2 3 32 32 H 1.08999 * 109.46792 * 50.24094 * 1 2 3 33 33 H 1.09002 * 109.46895 * 180.02562 * 6 4 3 34 34 H 1.09005 * 109.47121 * 300.00359 * 6 4 3 35 35 H 1.08997 * 109.47114 * 60.00315 * 6 4 3 36 36 H 0.96996 * 120.00299 * 19.94336 * 7 4 3 37 37 H 1.09004 * 109.58968 * 6.60403 * 11 10 8 38 38 H 1.09008 * 109.70714 * 246.25101 * 11 10 8 39 39 H 1.09005 * 109.49145 * 294.68135 * 12 11 10 40 40 H 1.09000 * 109.49754 * 174.59109 * 12 11 10 41 41 H 1.09002 * 109.46294 * 178.61766 * 13 12 11 42 42 H 1.09003 * 109.44670 * 58.65242 * 13 12 11 43 43 H 0.96999 * 120.00520 * 180.02562 * 19 18 16 44 44 H 1.09004 * 109.55718 * 113.80267 * 24 10 8 45 45 H 1.08995 * 109.61163 * 353.38010 * 24 10 8 46 46 H 1.09010 * 109.51884 * 182.88142 * 25 3 2 47 47 H 1.09000 * 109.52629 * 303.06962 * 25 3 2 48 48 H 1.08998 * 109.40081 * 297.48280 * 26 25 3 49 49 H 1.09001 * 109.48100 * 177.50658 * 26 25 3 50 50 H 1.08998 * 109.47980 * 301.17474 * 28 27 26 51 51 H 1.08995 * 109.48179 * 181.13868 * 28 27 26 52 52 H 1.08999 * 109.52612 * 182.48932 * 29 28 27 53 53 H 1.09001 * 109.52590 * 62.29696 * 29 28 27 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.4290 0.0000 0.0000 3 6 2.0102 1.3054 0.0000 4 6 1.5692 2.0620 1.2546 5 1 2.1541 2.9766 1.3517 6 6 1.7916 1.1814 2.4858 7 7 0.1479 2.4004 1.1465 8 6 -0.3694 3.3869 1.9051 9 8 0.3582 4.0517 2.6171 10 7 -1.6935 3.6365 1.8746 11 6 -2.5699 2.9154 0.9415 12 6 -3.7174 2.2805 1.7335 13 6 -4.4257 3.3615 2.5530 14 6 -3.4386 3.9807 3.5452 15 1 -3.0433 3.2028 4.1985 16 6 -4.1452 5.0225 4.3738 17 1 -5.1937 5.2225 4.2102 18 6 -3.4514 5.7242 5.3368 19 7 -2.1745 5.4806 5.5360 20 14 -4.3253 7.0209 6.3586 21 1 -5.4417 6.3909 7.1083 22 1 -3.3655 7.6259 7.3168 23 1 -4.8637 8.0765 5.4636 24 6 -2.2876 4.6348 2.7744 25 6 3.5373 1.1864 -0.0163 26 6 3.9739 0.5002 -1.3135 27 8 3.4905 1.2452 -2.4330 28 6 2.0665 1.3475 -2.4949 29 6 1.5534 2.0697 -1.2464 30 1 -0.3633 -1.0128 -0.1743 31 1 -0.3634 0.3554 0.9642 32 1 -0.3633 0.6573 -0.7900 33 1 1.4771 1.7202 3.3796 34 1 1.2067 0.2666 2.3887 35 1 2.8491 0.9296 2.5663 36 1 -0.4207 1.9144 0.5289 37 1 -2.0011 2.1362 0.4340 38 1 -2.9746 3.6104 0.2057 39 1 -3.3189 1.5185 2.4035 40 1 -4.4267 1.8238 1.0433 41 1 -5.2580 2.9167 3.0985 42 1 -4.8021 4.1356 1.8842 43 1 -1.6865 5.9740 6.2137 44 1 -2.6685 5.4721 2.1897 45 1 -1.5329 4.9898 3.4762 46 1 3.9807 2.1809 0.0353 47 1 3.8658 0.5953 0.8385 48 1 3.5645 -0.5095 -1.3446 49 1 5.0622 0.4525 -1.3511 50 1 1.6321 0.3489 -2.5415 51 1 1.7808 1.9102 -3.3835 52 1 0.4644 2.1097 -1.2703 53 1 1.9555 3.0825 -1.2203 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 ZINC000595643421.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 = -104.236844 kcal Electronic energy + Delta-G solvation = -32841.129178 eV Core-core repulsion = 28713.646871 eV Total energy + Delta-G solvation = -4127.482307 eV No. of doubly occupied orbitals = 66 Molecular weight (most abundant/longest-lived isotopes) = 340.220 amu Computer time = 0.69 seconds Orbital eigenvalues (eV) -41.12719 -39.47495 -37.78032 -36.78823 -34.62443 -33.55291 -32.67294 -31.73664 -31.54639 -30.01961 -29.39129 -27.20174 -25.56586 -25.08844 -24.25047 -23.66970 -21.93834 -21.67029 -20.54510 -19.90980 -18.79020 -17.71432 -16.69513 -16.47109 -16.24349 -16.07585 -15.87751 -15.66980 -15.17514 -14.73451 -14.62331 -14.57751 -14.07185 -13.80585 -13.56063 -13.40825 -13.06772 -12.95705 -12.84624 -12.67298 -12.61625 -12.55341 -12.12061 -11.98667 -11.79821 -11.71935 -11.49702 -11.24748 -11.00080 -10.91561 -10.77585 -10.58966 -10.53050 -10.43972 -10.23987 -9.99149 -9.90235 -9.57062 -9.38728 -9.11938 -8.96908 -8.77669 -8.42175 -7.67327 -5.98751 -4.07905 3.12944 3.32622 3.38315 3.39407 3.41764 3.64389 3.82785 4.25617 4.35831 4.41154 4.44770 4.48681 4.71854 4.86966 4.98817 5.10512 5.13746 5.18661 5.24055 5.30825 5.31086 5.37145 5.40055 5.43554 5.45080 5.46624 5.55259 5.58355 5.71093 5.75737 5.82723 5.88057 5.92535 6.06884 6.12157 6.27288 6.39516 6.55267 6.62743 6.76653 6.79705 6.90051 6.94838 6.96336 7.03391 7.32144 7.46538 7.72977 7.91565 8.05483 8.40173 8.54719 8.95268 9.23197 9.57637 11.06245 Molecular weight = 340.22amu Principal moments of inertia in cm(-1) A = 0.021222 B = 0.002941 C = 0.002773 Principal moments of inertia in units of 10**(-40)*gram-cm**2 A = 1319.039575 B = 9518.124286 C =10094.118268 Net atomic charges, atomic populations, and dipole contributions calculated with CM2 Atom NO. Type Charge No. of electrons 1 C 0.005 3.995 2 O -0.368 6.368 3 C 0.111 3.889 4 C 0.142 3.858 5 H 0.092 0.908 6 C -0.159 4.159 7 N -0.760 5.760 8 C 0.699 3.301 9 O -0.581 6.581 10 N -0.622 5.622 11 C 0.116 3.884 12 C -0.134 4.134 13 C -0.109 4.109 14 C 0.040 3.960 15 H 0.028 0.972 16 C -0.523 4.523 17 H 0.057 0.943 18 C 0.063 3.937 19 N -0.879 5.879 20 Si 0.892 3.108 21 H -0.278 1.278 22 H -0.287 1.287 23 H -0.277 1.277 24 C 0.130 3.870 25 C -0.138 4.138 26 C 0.065 3.935 27 O -0.383 6.383 28 C 0.067 3.933 29 C -0.181 4.181 30 H 0.089 0.911 31 H 0.061 0.939 32 H 0.045 0.955 33 H 0.067 0.933 34 H 0.059 0.941 35 H 0.061 0.939 36 H 0.391 0.609 37 H 0.053 0.947 38 H 0.061 0.939 39 H 0.059 0.941 40 H 0.058 0.942 41 H 0.051 0.949 42 H 0.049 0.951 43 H 0.264 0.736 44 H 0.038 0.962 45 H 0.123 0.877 46 H 0.089 0.911 47 H 0.088 0.912 48 H 0.059 0.941 49 H 0.094 0.906 50 H 0.055 0.945 51 H 0.093 0.907 52 H 0.077 0.923 53 H 0.086 0.914 For the charges calculated by CM2: Dipole moment (debyes) X Y Z Total from point charges 7.899 -12.450 -15.755 21.578 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.092 4.092 2 O -0.286 6.286 3 C 0.070 3.930 4 C 0.036 3.964 5 H 0.110 0.890 6 C -0.218 4.218 7 N -0.425 5.425 8 C 0.401 3.599 9 O -0.458 6.458 10 N -0.359 5.359 11 C -0.009 4.009 12 C -0.172 4.172 13 C -0.147 4.147 14 C 0.021 3.979 15 H 0.046 0.954 16 C -0.561 4.561 17 H 0.076 0.924 18 C -0.140 4.140 19 N -0.516 5.516 20 Si 0.720 3.280 21 H -0.203 1.203 22 H -0.213 1.213 23 H -0.203 1.203 24 C 0.009 3.991 25 C -0.176 4.176 26 C -0.012 4.012 27 O -0.302 6.302 28 C -0.010 4.010 29 C -0.218 4.218 30 H 0.108 0.892 31 H 0.080 0.920 32 H 0.064 0.936 33 H 0.086 0.914 34 H 0.078 0.922 35 H 0.080 0.920 36 H 0.232 0.768 37 H 0.071 0.929 38 H 0.079 0.921 39 H 0.078 0.922 40 H 0.077 0.923 41 H 0.070 0.930 42 H 0.068 0.932 43 H 0.074 0.926 44 H 0.057 0.943 45 H 0.141 0.859 46 H 0.107 0.893 47 H 0.107 0.893 48 H 0.077 0.923 49 H 0.112 0.888 50 H 0.073 0.927 51 H 0.111 0.889 52 H 0.095 0.905 53 H 0.105 0.895 Dipole moment (debyes) X Y Z Total from point charges 8.810 -11.957 -15.224 21.269 hybrid contribution -1.314 0.596 0.542 1.541 sum 7.497 -11.362 -14.682 20.021 Atomic orbital electron populations 1.23222 0.81893 1.03047 1.01007 1.87940 1.19243 1.26661 1.94795 1.21810 0.94321 0.83434 0.93388 1.21272 0.82221 0.98602 0.94277 0.88978 1.22140 1.01596 0.99833 0.98196 1.46115 1.10254 1.39603 1.46528 1.16003 0.82685 0.80479 0.80730 1.90914 1.58532 1.51633 1.44687 1.47796 1.04322 1.43198 1.40589 1.21483 0.92276 0.95003 0.92116 1.21955 0.96963 0.99218 0.99060 1.21610 0.98828 0.96269 0.98010 1.18977 0.91204 0.95645 0.92062 0.95405 1.21029 0.96061 1.20848 1.18181 0.92447 1.25051 0.94898 0.97093 0.96952 1.69878 1.07562 1.41012 1.33158 0.86656 0.79336 0.81869 0.80116 1.20323 1.21259 1.20261 1.22196 0.98797 0.87571 0.90535 1.22198 0.92084 1.03357 0.99968 1.22530 0.99354 0.94627 0.84707 1.87916 1.21054 1.74903 1.46334 1.22487 0.81391 0.99839 0.97290 1.22764 1.02150 1.01300 0.95634 0.89228 0.92046 0.93640 0.91392 0.92191 0.91972 0.76798 0.92917 0.92055 0.92180 0.92312 0.93010 0.93174 0.92599 0.94328 0.85888 0.89286 0.89339 0.92278 0.88761 0.92700 0.88932 0.90481 0.89495 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.00 0.04 8.81 101.05 0.89 0.93 16 2 O -0.37 -3.81 8.39 -35.23 -0.30 -4.11 16 3 C 0.11 1.16 0.37 -90.43 -0.03 1.12 16 4 C 0.14 1.91 1.76 -67.93 -0.12 1.79 16 5 H 0.09 1.40 7.64 -51.93 -0.40 1.00 16 6 C -0.16 -2.16 8.88 37.16 0.33 -1.83 16 7 N -0.76 -11.29 3.44 -58.68 -0.20 -11.49 16 8 C 0.70 14.13 8.18 -86.92 -0.71 13.41 16 9 O -0.58 -14.37 15.40 5.29 0.08 -14.29 16 10 N -0.62 -12.42 2.97 -178.48 -0.53 -12.95 16 11 C 0.12 1.75 6.43 -3.71 -0.02 1.72 16 12 C -0.13 -2.08 6.08 -26.61 -0.16 -2.24 16 13 C -0.11 -2.13 5.11 -26.69 -0.14 -2.27 16 14 C 0.04 0.99 2.25 -91.65 -0.21 0.78 16 15 H 0.03 0.77 8.14 -51.93 -0.42 0.34 16 16 C -0.52 -14.10 8.57 -34.44 -0.30 -14.39 16 17 H 0.06 1.51 7.74 -52.49 -0.41 1.10 16 18 C 0.06 1.75 5.30 -17.82 -0.09 1.65 16 19 N -0.88 -25.66 11.39 55.43 0.63 -25.03 16 20 Si 0.89 21.27 29.54 -169.99 -5.02 16.25 16 21 H -0.28 -6.53 7.11 56.52 0.40 -6.12 16 22 H -0.29 -6.85 7.11 56.52 0.40 -6.45 16 23 H -0.28 -6.50 7.11 56.52 0.40 -6.10 16 24 C 0.13 3.14 5.38 -3.70 -0.02 3.12 16 25 C -0.14 -1.26 4.91 -26.58 -0.13 -1.39 16 26 C 0.06 0.56 6.69 37.21 0.25 0.81 16 27 O -0.38 -4.02 10.77 -35.23 -0.38 -4.40 16 28 C 0.07 0.60 6.69 37.20 0.25 0.85 16 29 C -0.18 -1.78 4.31 -26.59 -0.11 -1.90 16 30 H 0.09 0.66 8.14 -51.93 -0.42 0.24 16 31 H 0.06 0.66 5.11 -51.93 -0.27 0.39 16 32 H 0.05 0.36 5.38 -51.93 -0.28 0.08 16 33 H 0.07 1.15 8.03 -51.93 -0.42 0.73 16 34 H 0.06 0.78 6.46 -51.93 -0.34 0.45 16 35 H 0.06 0.71 6.86 -51.93 -0.36 0.35 16 36 H 0.39 4.84 1.48 -40.82 -0.06 4.77 16 37 H 0.05 0.63 5.56 -51.93 -0.29 0.34 16 38 H 0.06 0.87 8.14 -51.92 -0.42 0.45 16 39 H 0.06 0.97 8.14 -51.93 -0.42 0.55 16 40 H 0.06 0.73 8.14 -51.93 -0.42 0.31 16 41 H 0.05 0.99 8.14 -51.93 -0.42 0.57 16 42 H 0.05 0.95 8.14 -51.93 -0.42 0.53 16 43 H 0.26 7.46 8.31 -40.82 -0.34 7.12 16 44 H 0.04 0.94 8.14 -51.93 -0.42 0.51 16 45 H 0.12 3.52 5.14 -51.93 -0.27 3.25 16 46 H 0.09 0.80 8.14 -51.92 -0.42 0.38 16 47 H 0.09 0.74 6.93 -51.93 -0.36 0.38 16 48 H 0.06 0.52 7.82 -51.93 -0.41 0.11 16 49 H 0.09 0.66 8.14 -51.93 -0.42 0.24 16 50 H 0.05 0.48 7.82 -51.93 -0.41 0.07 16 51 H 0.09 0.72 8.14 -51.93 -0.42 0.29 16 52 H 0.08 0.73 5.22 -51.93 -0.27 0.46 16 53 H 0.09 0.94 8.14 -51.93 -0.42 0.52 16 LS Contribution 386.19 15.07 5.82 5.82 Total: -1.00 -32.19 386.19 -8.95 -41.14 By element: Atomic # 1 Polarization: 15.60 SS G_CDS: -8.72 Total: 6.88 kcal Atomic # 6 Polarization: 2.51 SS G_CDS: -0.33 Total: 2.18 kcal Atomic # 7 Polarization: -49.36 SS G_CDS: -0.10 Total: -49.46 kcal Atomic # 8 Polarization: -22.21 SS G_CDS: -0.59 Total: -22.80 kcal Atomic # 14 Polarization: 21.27 SS G_CDS: -5.02 Total: 16.25 kcal Total LS contribution 5.82 Total: 5.82 kcal Total: -32.19 -8.95 -41.14 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. ZINC000595643421.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 -63.095 kcal (2) G-P(sol) polarization free energy of solvation -32.192 kcal (3) G-ENP(sol) elect.-nuc.-pol. free energy of system -95.287 kcal (4) G-CDS(sol) cavity-dispersion-solvent structure free energy -8.950 kcal (5) G-P-CDS(sol) = G-P(sol) + G-CDS(sol) = (2) + (4) -41.142 kcal (6) G-S(sol) free energy of system = (1) + (5) -104.237 kcal 1SCF WAS SPECIFIED, SO GEOMETRY OPTIMIZATION WAS NOT USED Total computer time = 0.69 seconds