Wall clock time and date at job start Sat Apr 11 2020 06:12:36 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.42894 * 1 3 3 C 1.42900 * 114.00153 * 2 1 4 4 C 1.53002 * 109.47349 * 179.97438 * 3 2 1 5 5 C 1.52992 * 109.47167 * 184.27470 * 4 3 2 6 6 C 1.50708 * 115.55268 * 304.97331 * 5 4 3 7 7 O 1.21278 * 120.00155 * 102.30461 * 6 5 4 8 8 N 1.34778 * 119.99969 * 282.30752 * 6 5 4 9 9 C 1.46919 * 120.63069 * 187.42971 * 8 6 5 10 10 C 1.53191 * 108.77574 * 233.63268 * 9 8 6 11 11 C 1.53041 * 109.31105 * 305.36347 * 10 9 8 12 12 N 1.46506 * 109.46291 * 181.38450 * 11 10 9 13 13 C 1.36935 * 119.99958 * 297.01031 * 12 11 10 14 14 H 1.08004 * 120.00207 * 324.77593 * 13 12 11 15 15 C 1.38814 * 120.00039 * 144.77595 * 13 12 11 16 16 N 1.31613 * 119.99856 * 156.20753 * 15 13 12 17 17 Si 1.86796 * 119.99977 * 336.20344 * 15 13 12 18 18 H 1.48510 * 109.47021 * 91.97615 * 17 15 13 19 19 H 1.48493 * 109.47580 * 211.97327 * 17 15 13 20 20 H 1.48497 * 109.47201 * 331.98099 * 17 15 13 21 21 C 1.46503 * 119.99737 * 117.00649 * 12 11 10 22 22 C 1.52994 * 109.47130 * 90.00591 * 21 12 11 23 23 O 1.42895 * 109.47423 * 64.99244 * 22 21 12 24 24 C 1.53046 * 109.53985 * 61.36632 * 11 10 9 25 25 C 1.46922 * 120.62756 * 7.45225 * 8 6 5 26 26 C 1.53005 * 117.50163 * 90.64140 * 5 4 3 27 27 C 1.52997 * 117.50110 * 159.30269 * 5 4 3 28 28 H 1.08997 * 109.47404 * 59.99856 * 1 2 3 29 29 H 1.09003 * 109.47505 * 180.02562 * 1 2 3 30 30 H 1.09007 * 109.47339 * 299.99991 * 1 2 3 31 31 H 1.09004 * 109.47246 * 59.99614 * 3 2 1 32 32 H 1.09003 * 109.47180 * 300.00200 * 3 2 1 33 33 H 1.09000 * 109.46737 * 64.28075 * 4 3 2 34 34 H 1.08994 * 109.46947 * 304.27755 * 4 3 2 35 35 H 1.08994 * 109.58996 * 113.84378 * 9 8 6 36 36 H 1.09005 * 109.70348 * 353.48517 * 9 8 6 37 37 H 1.09002 * 109.50060 * 185.41240 * 10 9 8 38 38 H 1.09000 * 109.50073 * 65.31028 * 10 9 8 39 39 H 1.08996 * 109.46345 * 301.34312 * 11 10 9 40 40 H 0.97010 * 119.99982 * 179.97438 * 16 15 13 41 41 H 1.09003 * 109.46875 * 210.00272 * 21 12 11 42 42 H 1.08992 * 109.47065 * 329.99997 * 21 12 11 43 43 H 1.09005 * 109.47237 * 184.99965 * 22 21 12 44 44 H 1.09006 * 109.47099 * 304.99549 * 22 21 12 45 45 H 0.96698 * 114.00362 * 179.97438 * 23 22 21 46 46 H 1.08998 * 109.49579 * 58.58007 * 24 11 10 47 47 H 1.08999 * 109.49708 * 178.68442 * 24 11 10 48 48 H 1.08999 * 109.58541 * 246.15494 * 25 8 6 49 49 H 1.08995 * 109.59047 * 6.58064 * 25 8 6 50 50 H 1.09000 * 117.50141 * 215.00535 * 26 5 4 51 51 H 1.08997 * 117.49027 * 0.02562 * 26 5 4 52 52 H 1.08996 * 117.49856 * 359.97438 * 27 5 4 53 53 H 1.09003 * 117.50223 * 144.99593 * 27 5 4 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.4289 0.0000 0.0000 3 6 2.0102 1.3054 0.0000 4 6 3.5354 1.1846 0.0006 5 6 4.1574 2.5783 0.1082 6 6 3.7166 3.5694 -0.9382 7 8 2.8994 4.4214 -0.6604 8 7 4.2304 3.5087 -2.1827 9 6 3.9371 4.5604 -3.1657 10 6 3.3834 3.9073 -4.4361 11 6 4.3674 2.8436 -4.9284 12 7 3.8568 2.2398 -6.1617 13 6 2.6775 1.5441 -6.1462 14 1 2.4093 0.9551 -5.2815 15 6 1.8266 1.5958 -7.2418 16 7 0.9679 0.6214 -7.4546 17 14 1.9004 3.0484 -8.4139 18 1 2.8386 2.7419 -9.5235 19 1 0.5470 3.3072 -8.9673 20 1 2.3745 4.2508 -7.6828 21 6 4.6079 2.3805 -7.4116 22 6 5.5896 1.2152 -7.5506 23 8 4.8620 -0.0075 -7.6823 24 6 4.5277 1.7605 -3.8590 25 6 5.0943 2.3894 -2.5821 26 6 4.4499 3.1152 1.5107 27 6 5.6014 2.6712 0.6051 28 1 -0.3634 0.5138 -0.8899 29 1 -0.3634 -1.0277 0.0005 30 1 -0.3634 0.5139 0.8900 31 1 1.6885 1.8464 0.8900 32 1 1.6886 1.8464 -0.8900 33 1 3.8627 0.7112 -0.9250 34 1 3.8516 0.5792 0.8500 35 1 4.8514 5.1041 -3.4032 36 1 3.1978 5.2495 -2.7572 37 1 3.2518 4.6660 -5.2075 38 1 2.4231 3.4407 -4.2164 39 1 5.3345 3.3064 -5.1243 40 1 0.3730 0.6577 -8.2200 41 1 3.9160 2.3751 -8.2539 42 1 5.1592 3.3207 -7.4004 43 1 6.2098 1.3660 -8.4342 44 1 6.2239 1.1670 -6.6654 45 1 5.4219 -0.7904 -7.7752 46 1 3.5566 1.3143 -3.6447 47 1 5.2106 0.9914 -4.2197 48 1 6.1030 2.7560 -2.7719 49 1 5.1177 1.6441 -1.7871 50 1 4.3235 4.1850 1.6767 51 1 4.1996 2.4752 2.3568 52 1 6.1082 1.7393 0.8559 53 1 6.2326 3.4489 0.1753 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 ZINC001098387896.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:12:36 Heat of formation + Delta-G solvation = -70.004432 kcal Electronic energy + Delta-G solvation = -32446.925687 eV Core-core repulsion = 28320.927809 eV Total energy + Delta-G solvation = -4125.997879 eV No. of doubly occupied orbitals = 66 Molecular weight (most abundant/longest-lived isotopes) = 340.220 amu Computer time = 1.71 seconds Orbital eigenvalues (eV) -40.55115 -38.24761 -37.64697 -36.73035 -36.21649 -34.68382 -31.96002 -31.26348 -30.85318 -30.63467 -27.93578 -26.55743 -25.39787 -24.02240 -23.32479 -22.83304 -22.60427 -21.91661 -20.86880 -19.67696 -18.42630 -16.87333 -16.70292 -16.46750 -16.26654 -16.13442 -15.74141 -15.31572 -15.16793 -14.81913 -14.71989 -14.42859 -14.01385 -13.68680 -13.59150 -13.36494 -13.17232 -12.86680 -12.71460 -12.41536 -12.39016 -12.13896 -11.97211 -11.82104 -11.70883 -11.60631 -11.50701 -11.18285 -11.08302 -10.78313 -10.70595 -10.66549 -10.51798 -10.26937 -10.17835 -9.93285 -9.65843 -9.61410 -9.49496 -9.45249 -9.07550 -8.46363 -8.29936 -6.79140 -5.83673 -3.21100 2.66584 3.05078 3.36518 3.41245 3.52096 3.54441 3.79311 4.41396 4.57522 4.63407 4.66874 4.78121 4.81493 5.08167 5.10944 5.17422 5.24696 5.36290 5.39159 5.40926 5.42560 5.45797 5.50149 5.55654 5.59661 5.64691 5.69266 5.70765 5.75848 5.84240 5.88637 5.89068 5.93044 5.94698 5.95860 6.09442 6.13007 6.17448 6.19770 6.34958 6.45378 6.57932 6.74902 6.80466 7.00647 7.05559 7.17769 7.57090 8.14267 8.29130 8.46130 8.53323 9.08011 9.86631 10.20841 11.32767 Molecular weight = 340.22amu Principal moments of inertia in cm(-1) A = 0.011712 B = 0.003854 C = 0.003500 Principal moments of inertia in units of 10**(-40)*gram-cm**2 A = 2390.052960 B = 7263.781623 C = 7998.389262 Net atomic charges, atomic populations, and dipole contributions calculated with CM2 Atom NO. Type Charge No. of electrons 1 C 0.033 3.967 2 O -0.392 6.392 3 C 0.064 3.936 4 C -0.075 4.075 5 C -0.152 4.152 6 C 0.551 3.449 7 O -0.544 6.544 8 N -0.598 5.598 9 C 0.113 3.887 10 C -0.141 4.141 11 C 0.181 3.819 12 N -0.593 5.593 13 C -0.244 4.244 14 H 0.077 0.923 15 C 0.043 3.957 16 N -0.897 5.897 17 Si 0.870 3.130 18 H -0.288 1.288 19 H -0.306 1.306 20 H -0.273 1.273 21 C 0.129 3.871 22 C 0.063 3.937 23 O -0.533 6.533 24 C -0.138 4.138 25 C 0.107 3.893 26 C -0.146 4.146 27 C -0.157 4.157 28 H 0.045 0.955 29 H 0.086 0.914 30 H 0.038 0.962 31 H 0.050 0.950 32 H 0.052 0.948 33 H 0.085 0.915 34 H 0.081 0.919 35 H 0.061 0.939 36 H 0.085 0.915 37 H 0.067 0.933 38 H 0.076 0.924 39 H 0.050 0.950 40 H 0.259 0.741 41 H 0.064 0.936 42 H 0.050 0.950 43 H 0.020 0.980 44 H 0.030 0.970 45 H 0.352 0.648 46 H 0.075 0.925 47 H 0.070 0.930 48 H 0.062 0.938 49 H 0.084 0.916 50 H 0.101 0.899 51 H 0.103 0.897 52 H 0.105 0.895 53 H 0.099 0.901 For the charges calculated by CM2: Dipole moment (debyes) X Y Z Total from point charges 11.191 2.553 14.127 18.203 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.064 4.064 2 O -0.311 6.311 3 C -0.013 4.013 4 C -0.113 4.113 5 C -0.156 4.156 6 C 0.338 3.662 7 O -0.422 6.422 8 N -0.331 5.331 9 C -0.010 4.010 10 C -0.181 4.181 11 C 0.073 3.927 12 N -0.314 5.314 13 C -0.373 4.373 14 H 0.094 0.906 15 C -0.157 4.157 16 N -0.539 5.539 17 Si 0.701 3.299 18 H -0.217 1.217 19 H -0.234 1.234 20 H -0.198 1.198 21 C 0.003 3.997 22 C -0.017 4.017 23 O -0.334 6.334 24 C -0.178 4.178 25 C -0.016 4.016 26 C -0.183 4.183 27 C -0.194 4.194 28 H 0.064 0.936 29 H 0.104 0.896 30 H 0.056 0.944 31 H 0.069 0.931 32 H 0.071 0.929 33 H 0.103 0.897 34 H 0.100 0.900 35 H 0.079 0.921 36 H 0.103 0.897 37 H 0.085 0.915 38 H 0.095 0.905 39 H 0.068 0.932 40 H 0.069 0.931 41 H 0.082 0.918 42 H 0.068 0.932 43 H 0.039 0.961 44 H 0.049 0.951 45 H 0.195 0.805 46 H 0.093 0.907 47 H 0.088 0.912 48 H 0.080 0.920 49 H 0.102 0.898 50 H 0.119 0.881 51 H 0.121 0.879 52 H 0.123 0.877 53 H 0.117 0.883 Dipole moment (debyes) X Y Z Total from point charges 10.188 3.204 14.199 17.767 hybrid contribution 0.275 0.588 -1.061 1.244 sum 10.463 3.792 13.138 17.218 Atomic orbital electron populations 1.22731 0.80775 1.02695 1.00153 1.87809 1.19971 1.27850 1.95499 1.22434 0.93322 0.84627 1.00867 1.21036 0.93670 0.92434 1.04148 1.21286 1.00853 0.96768 0.96649 1.19454 0.80482 0.84015 0.82212 1.90670 1.32724 1.38271 1.80488 1.48057 1.49255 1.24637 1.11136 1.21704 1.00473 0.89755 0.89073 1.22377 1.01641 0.99216 0.94881 1.20304 0.94346 0.91067 0.86956 1.44104 1.23032 1.61362 1.02875 1.18472 0.93748 1.20966 1.04068 0.90594 1.25261 0.96205 0.96730 0.97484 1.69893 1.24735 1.30156 1.29071 0.87133 0.75278 0.86160 0.81298 1.21657 1.23402 1.19815 1.21396 0.93898 0.93282 0.91146 1.22868 0.94456 0.87849 0.96575 1.86969 1.35080 1.15033 1.96288 1.22317 1.02358 0.99248 0.93897 1.21932 0.93000 0.88277 0.98412 1.23043 1.01022 1.02306 0.91920 1.23112 0.89729 1.03795 1.02774 0.93612 0.89598 0.94361 0.93142 0.92942 0.89698 0.90031 0.92121 0.89665 0.91471 0.90500 0.93218 0.93144 0.91809 0.93186 0.96145 0.95149 0.80485 0.90661 0.91163 0.91980 0.89780 0.88103 0.87854 0.87660 0.88270 Number of geometries 1 Number of calculations of the screened coulomb radii 1 The total number of SCF iterations 17. 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.40 11.01 101.05 1.11 1.52 16 2 O -0.39 -5.79 10.70 -35.23 -0.38 -6.16 16 3 C 0.06 0.90 5.61 37.16 0.21 1.10 16 4 C -0.08 -0.89 5.27 -26.73 -0.14 -1.03 16 5 C -0.15 -1.76 0.84 -155.66 -0.13 -1.89 16 6 C 0.55 8.28 5.89 -10.98 -0.06 8.22 16 7 O -0.54 -9.99 16.40 5.56 0.09 -9.90 16 8 N -0.60 -8.63 2.97 -172.78 -0.51 -9.15 16 9 C 0.11 1.70 6.44 -3.72 -0.02 1.67 16 10 C -0.14 -2.54 5.35 -26.61 -0.14 -2.68 16 11 C 0.18 3.25 2.59 -67.88 -0.18 3.07 16 12 N -0.59 -13.17 2.04 -165.61 -0.34 -13.51 16 13 C -0.24 -6.58 8.29 -15.06 -0.12 -6.71 16 14 H 0.08 2.26 6.90 -52.48 -0.36 1.90 16 15 C 0.04 1.24 5.29 -17.43 -0.09 1.15 16 16 N -0.90 -27.77 13.97 55.49 0.78 -26.99 16 17 Si 0.87 22.37 25.18 -169.99 -4.28 18.09 16 18 H -0.29 -7.28 6.49 56.52 0.37 -6.91 16 19 H -0.31 -8.10 7.11 56.52 0.40 -7.69 16 20 H -0.27 -6.76 7.11 56.52 0.40 -6.36 16 21 C 0.13 2.69 4.91 -4.04 -0.02 2.67 16 22 C 0.06 1.15 7.39 37.15 0.27 1.43 16 23 O -0.53 -11.16 13.13 -35.23 -0.46 -11.62 16 24 C -0.14 -2.26 5.25 -26.61 -0.14 -2.40 16 25 C 0.11 1.36 5.99 -3.71 -0.02 1.34 16 26 C -0.15 -1.40 9.82 -26.69 -0.26 -1.67 16 27 C -0.16 -1.36 9.25 -26.69 -0.25 -1.60 16 28 H 0.05 0.59 8.14 -51.93 -0.42 0.17 16 29 H 0.09 0.90 8.14 -51.93 -0.42 0.48 16 30 H 0.04 0.40 8.14 -51.93 -0.42 -0.02 16 31 H 0.05 0.67 8.14 -51.93 -0.42 0.24 16 32 H 0.05 0.86 7.60 -51.93 -0.39 0.46 16 33 H 0.08 1.05 6.39 -51.93 -0.33 0.72 16 34 H 0.08 0.82 8.09 -51.93 -0.42 0.40 16 35 H 0.06 0.79 8.14 -51.93 -0.42 0.37 16 36 H 0.08 1.38 7.01 -51.93 -0.36 1.01 16 37 H 0.07 1.25 8.14 -51.93 -0.42 0.82 16 38 H 0.08 1.63 7.77 -51.93 -0.40 1.23 16 39 H 0.05 0.80 7.82 -51.93 -0.41 0.39 16 40 H 0.26 7.69 8.31 -40.82 -0.34 7.35 16 41 H 0.06 1.53 2.97 -51.93 -0.15 1.38 16 42 H 0.05 0.90 7.94 -51.93 -0.41 0.49 16 43 H 0.02 0.32 8.14 -51.93 -0.42 -0.10 16 44 H 0.03 0.50 8.14 -51.93 -0.42 0.08 16 45 H 0.35 5.39 9.12 45.56 0.42 5.81 16 46 H 0.07 1.48 6.95 -51.93 -0.36 1.11 16 47 H 0.07 1.14 8.14 -51.93 -0.42 0.72 16 48 H 0.06 0.68 8.14 -51.93 -0.42 0.26 16 49 H 0.08 0.95 4.24 -51.93 -0.22 0.73 16 50 H 0.10 1.05 8.14 -51.93 -0.42 0.63 16 51 H 0.10 0.83 8.14 -51.93 -0.42 0.41 16 52 H 0.11 0.72 8.14 -51.93 -0.42 0.30 16 53 H 0.10 0.85 8.14 -51.93 -0.42 0.42 16 LS Contribution 409.37 15.07 6.17 6.17 Total: -1.00 -34.65 409.37 -7.43 -42.08 By element: Atomic # 1 Polarization: 15.31 SS G_CDS: -8.50 Total: 6.81 kcal Atomic # 6 Polarization: 4.18 SS G_CDS: 0.01 Total: 4.19 kcal Atomic # 7 Polarization: -49.58 SS G_CDS: -0.08 Total: -49.65 kcal Atomic # 8 Polarization: -26.94 SS G_CDS: -0.75 Total: -27.69 kcal Atomic # 14 Polarization: 22.37 SS G_CDS: -4.28 Total: 18.09 kcal Total LS contribution 6.17 Total: 6.17 kcal Total: -34.65 -7.43 -42.08 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. ZINC001098387896.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 -27.927 kcal (2) G-P(sol) polarization free energy of solvation -34.648 kcal (3) G-ENP(sol) elect.-nuc.-pol. free energy of system -62.575 kcal (4) G-CDS(sol) cavity-dispersion-solvent structure free energy -7.429 kcal (5) G-P-CDS(sol) = G-P(sol) + G-CDS(sol) = (2) + (4) -42.078 kcal (6) G-S(sol) free energy of system = (1) + (5) -70.004 kcal 1SCF WAS SPECIFIED, SO GEOMETRY OPTIMIZATION WAS NOT USED Total computer time = 1.71 seconds