Wall clock time and date at job start Tue Mar 31 2020 19:33:58 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 C 1.53002 * 1 3 3 C 1.52995 * 109.47102 * 2 1 4 4 C 1.52995 * 109.47536 * 120.00088 * 2 1 3 5 5 C 1.52997 * 109.47291 * 174.79820 * 4 2 1 6 6 C 1.53004 * 110.54295 * 323.21245 * 5 4 2 7 7 C 1.54226 * 110.60892 * 200.36634 * 5 4 2 8 8 C 1.51310 * 103.55917 * 241.48313 * 7 5 4 9 9 O 1.21228 * 125.40219 * 180.02562 * 8 7 5 10 10 N 1.33879 * 109.20585 * 0.02562 * 8 7 5 11 11 C 1.46495 * 122.76482 * 179.97438 * 10 8 7 12 12 N 1.46901 * 109.46918 * 98.12270 * 11 10 8 13 13 C 1.46898 * 111.00054 * 275.88159 * 12 11 10 14 14 C 1.50704 * 109.47451 * 169.99911 * 13 12 11 15 15 H 1.07997 * 119.99562 * 359.73334 * 14 13 12 16 16 C 1.37880 * 120.00092 * 180.02562 * 14 13 12 17 17 N 1.31516 * 119.99895 * 359.72878 * 16 14 13 18 18 Si 1.86799 * 119.99801 * 179.72838 * 16 14 13 19 19 H 1.48497 * 109.47085 * 359.97438 * 18 16 14 20 20 H 1.48494 * 109.47140 * 120.00187 * 18 16 14 21 21 H 1.48502 * 109.46998 * 240.00233 * 18 16 14 22 22 C 1.46904 * 110.99879 * 151.93189 * 12 11 10 23 23 C 1.53001 * 109.47175 * 287.91525 * 22 12 11 24 24 C 1.52999 * 117.49873 * 243.80958 * 23 22 12 25 25 C 1.52996 * 117.49791 * 175.18506 * 23 22 12 26 26 C 1.33877 * 114.47301 * 0.26351 * 10 8 7 27 27 O 1.21232 * 125.39729 * 179.70863 * 26 10 8 28 28 H 1.08998 * 109.46772 * 299.99724 * 1 2 3 29 29 H 1.08992 * 109.47312 * 59.99573 * 1 2 3 30 30 H 1.08996 * 109.46832 * 180.02562 * 1 2 3 31 31 H 1.09007 * 109.46999 * 240.00232 * 2 1 3 32 32 H 1.09001 * 109.47259 * 196.72948 * 3 2 1 33 33 H 1.09003 * 109.47449 * 316.73168 * 3 2 1 34 34 H 1.08999 * 109.47605 * 76.73089 * 3 2 1 35 35 H 1.09006 * 109.46877 * 54.80067 * 4 2 1 36 36 H 1.08994 * 109.47199 * 294.80208 * 4 2 1 37 37 H 1.08995 * 109.47185 * 297.85426 * 6 5 4 38 38 H 1.09003 * 109.46790 * 57.85412 * 6 5 4 39 39 H 1.09000 * 109.47446 * 177.84954 * 6 5 4 40 40 H 1.08998 * 110.60797 * 122.96538 * 7 5 4 41 41 H 1.09005 * 110.73828 * 0.06393 * 7 5 4 42 42 H 1.09005 * 109.47279 * 338.12282 * 11 10 8 43 43 H 1.08997 * 109.47550 * 218.11848 * 11 10 8 44 44 H 1.09005 * 109.46904 * 50.00410 * 13 12 11 45 45 H 1.09000 * 109.47451 * 290.00161 * 13 12 11 46 46 H 0.97001 * 119.99460 * 179.97438 * 17 16 14 47 47 H 1.08995 * 109.47285 * 47.91848 * 22 12 11 48 48 H 1.09007 * 109.46987 * 167.91732 * 22 12 11 49 49 H 1.08997 * 115.54677 * 29.51364 * 23 22 12 50 50 H 1.09006 * 117.49909 * 359.97438 * 24 23 22 51 51 H 1.09000 * 117.49782 * 145.01629 * 24 23 22 52 52 H 1.09008 * 117.49538 * 214.97771 * 25 23 22 53 53 H 1.08998 * 117.50247 * 359.97438 * 25 23 22 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 6 1.5300 0.0000 0.0000 3 6 2.0400 1.4424 0.0000 4 6 2.0401 -0.7212 1.2492 5 6 3.5645 -0.8354 1.1854 6 6 4.0245 -1.0850 -0.2523 7 6 4.0642 -1.9578 2.1177 8 6 4.9453 -1.2295 3.1091 9 8 5.5452 -1.7459 4.0272 10 7 4.9574 0.0758 2.8116 11 6 5.7199 1.0583 3.5857 12 7 6.9944 1.3327 2.9088 13 6 8.0067 0.3299 3.2658 14 6 9.2120 0.4939 2.3761 15 1 9.2229 1.2631 1.6181 16 6 10.2983 -0.3419 2.5258 17 7 10.2847 -1.2789 3.4486 18 14 11.7891 -0.1435 1.4179 19 1 11.5553 0.9746 0.4690 20 1 12.9856 0.1516 2.2463 21 1 12.0112 -1.3994 0.6572 22 6 7.4689 2.6900 3.2099 23 6 6.5939 3.7101 2.4786 24 6 7.2632 4.5789 1.4118 25 6 6.8568 5.1921 2.7534 26 6 4.2105 0.4101 1.7520 27 8 4.0756 1.5324 1.3139 28 1 -0.3633 0.5138 0.8900 29 1 -0.3633 0.5139 -0.8899 30 1 -0.3633 -1.0276 0.0005 31 1 1.8934 -0.5138 -0.8901 32 1 3.0890 1.4570 -0.2958 33 1 1.4556 2.0345 -0.7044 34 1 1.9388 1.8637 1.0002 35 1 1.6027 -1.7185 1.2969 36 1 1.7548 -0.1565 2.1367 37 1 3.6156 -2.0316 -0.6055 38 1 3.6725 -0.2758 -0.8921 39 1 5.1134 -1.1254 -0.2836 40 1 4.6428 -2.6925 1.5578 41 1 3.2283 -2.4371 2.6273 42 1 5.9157 0.6628 4.5825 43 1 5.1458 1.9813 3.6674 44 1 7.5909 -0.6690 3.1337 45 1 8.3019 0.4660 4.3062 46 1 11.0489 -1.8671 3.5536 47 1 7.4107 2.8653 4.2841 48 1 8.5025 2.7951 2.8800 49 1 5.5618 3.4057 2.3053 50 1 8.3270 4.4241 1.2312 51 1 6.6714 4.8462 0.5363 52 1 5.9975 5.8628 2.7605 53 1 7.6527 5.4404 3.4554 RHF calculation, no. of doubly occupied orbitals= 65 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 ZINC000077164573.mol2 53 1SCF WAS SPECIFIED, SO GEOMETRY OPTIMIZATION WAS NOT USED SCF FIELD WAS ACHIEVED Wall clock time and date at job start Tue Mar 31 2020 19:33:58 Heat of formation + Delta-G solvation = -35.947831 kcal Electronic energy + Delta-G solvation = -31745.768750 eV Core-core repulsion = 27811.703728 eV Total energy + Delta-G solvation = -3934.065022 eV No. of doubly occupied orbitals = 65 Molecular weight (most abundant/longest-lived isotopes) = 336.220 amu Computer time = 0.86 seconds Orbital eigenvalues (eV) -41.89849 -38.60547 -37.81512 -36.91367 -35.69038 -35.13594 -33.21006 -32.26930 -28.69728 -27.55823 -27.26145 -26.90765 -26.44296 -24.34753 -23.63369 -22.00223 -21.35818 -21.07829 -20.09505 -19.59135 -18.59353 -16.87833 -16.69314 -16.21375 -16.00935 -15.47032 -15.23847 -15.08323 -14.75725 -14.22621 -14.14193 -14.02171 -13.87661 -13.83327 -13.59527 -13.01247 -12.77545 -12.49116 -12.40520 -12.32509 -12.28136 -12.00382 -11.97596 -11.89837 -11.56845 -11.37620 -11.20604 -11.13785 -10.87490 -10.79720 -10.72371 -10.68414 -10.48400 -10.41659 -10.11667 -9.90969 -9.75097 -9.54897 -9.41907 -9.32203 -9.10249 -8.44250 -7.41580 -6.01721 -4.04331 1.89429 2.80937 3.15408 3.29857 3.37481 3.42280 4.08638 4.19651 4.33155 4.38393 4.42109 4.61266 4.83002 4.88010 4.90138 5.06268 5.23085 5.23849 5.33048 5.35006 5.39840 5.40406 5.49922 5.54125 5.54949 5.59935 5.64420 5.71297 5.72384 5.79840 5.83885 5.89970 5.91348 5.93442 5.99798 6.02186 6.13002 6.17704 6.20341 6.32086 6.34146 6.34764 6.48827 6.66848 6.87087 7.18267 7.24415 7.33771 7.73140 7.76529 7.98479 8.12390 8.24578 8.46799 9.03503 9.59617 11.08440 Molecular weight = 336.22amu Principal moments of inertia in cm(-1) A = 0.010663 B = 0.004230 C = 0.003557 Principal moments of inertia in units of 10**(-40)*gram-cm**2 A = 2625.263638 B = 6617.968352 C = 7868.841147 Net atomic charges, atomic populations, and dipole contributions calculated with CM2 Atom NO. Type Charge No. of electrons 1 C -0.146 4.146 2 C -0.102 4.102 3 C -0.151 4.151 4 C -0.092 4.092 5 C -0.095 4.095 6 C -0.138 4.138 7 C -0.152 4.152 8 C 0.525 3.475 9 O -0.473 6.473 10 N -0.558 5.558 11 C 0.263 3.737 12 N -0.541 5.541 13 C 0.193 3.807 14 C -0.510 4.510 15 H 0.064 0.936 16 C 0.057 3.943 17 N -0.892 5.892 18 Si 0.895 3.105 19 H -0.272 1.272 20 H -0.285 1.285 21 H -0.284 1.284 22 C 0.099 3.901 23 C -0.140 4.140 24 C -0.172 4.172 25 C -0.190 4.190 26 C 0.541 3.459 27 O -0.467 6.467 28 H 0.051 0.949 29 H 0.053 0.947 30 H 0.051 0.949 31 H 0.078 0.922 32 H 0.084 0.916 33 H 0.038 0.962 34 H 0.054 0.946 35 H 0.076 0.924 36 H 0.066 0.934 37 H 0.071 0.929 38 H 0.078 0.922 39 H 0.066 0.934 40 H 0.119 0.881 41 H 0.114 0.886 42 H 0.066 0.934 43 H 0.089 0.911 44 H 0.036 0.964 45 H -0.007 1.007 46 H 0.261 0.739 47 H 0.027 0.973 48 H 0.093 0.907 49 H 0.109 0.891 50 H 0.093 0.907 51 H 0.086 0.914 52 H 0.086 0.914 53 H 0.089 0.911 For the charges calculated by CM2: Dipole moment (debyes) X Y Z Total from point charges -18.476 2.782 -5.063 19.358 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.204 4.204 2 C -0.120 4.120 3 C -0.208 4.208 4 C -0.130 4.130 5 C -0.098 4.098 6 C -0.195 4.195 7 C -0.192 4.192 8 C 0.313 3.687 9 O -0.344 6.344 10 N -0.301 5.301 11 C 0.052 3.948 12 N -0.267 5.267 13 C 0.069 3.931 14 C -0.548 4.548 15 H 0.082 0.918 16 C -0.145 4.145 17 N -0.532 5.532 18 Si 0.723 3.277 19 H -0.197 1.197 20 H -0.211 1.211 21 H -0.210 1.210 22 C -0.029 4.029 23 C -0.159 4.159 24 C -0.209 4.209 25 C -0.228 4.228 26 C 0.329 3.671 27 O -0.338 6.338 28 H 0.071 0.929 29 H 0.072 0.928 30 H 0.070 0.930 31 H 0.096 0.904 32 H 0.103 0.897 33 H 0.057 0.943 34 H 0.073 0.927 35 H 0.095 0.905 36 H 0.084 0.916 37 H 0.090 0.910 38 H 0.096 0.904 39 H 0.085 0.915 40 H 0.137 0.863 41 H 0.133 0.867 42 H 0.084 0.916 43 H 0.107 0.893 44 H 0.054 0.946 45 H 0.011 0.989 46 H 0.070 0.930 47 H 0.045 0.955 48 H 0.111 0.889 49 H 0.127 0.873 50 H 0.111 0.889 51 H 0.104 0.896 52 H 0.105 0.895 53 H 0.108 0.892 Dipole moment (debyes) X Y Z Total from point charges -18.150 2.956 -5.184 19.106 hybrid contribution 1.142 -0.018 0.370 1.201 sum -17.008 2.938 -4.814 17.919 Atomic orbital electron populations 1.21853 0.95557 1.01284 1.01666 1.21299 0.95202 0.97614 0.97921 1.22300 1.02720 0.94187 1.01611 1.21627 0.90856 1.01463 0.99043 1.21471 0.97474 0.92984 0.97902 1.21917 1.02305 1.02874 0.92438 1.22713 1.00920 0.98713 0.96867 1.20717 0.80769 0.82218 0.85044 1.90816 1.42622 1.70613 1.30303 1.46827 1.42906 1.10450 1.29889 1.22058 0.90344 0.91131 0.91310 1.61719 1.03394 1.03630 1.57970 1.19831 0.88759 0.89799 0.94729 1.21160 0.99211 1.16749 1.17713 0.91794 1.24988 0.98649 0.95229 0.95612 1.69999 1.30777 1.19327 1.33056 0.86490 0.82610 0.78044 0.80590 1.19651 1.21086 1.21026 1.21521 0.99212 0.86309 0.95881 1.22448 1.00896 0.90475 1.02062 1.22892 1.01202 1.04478 0.92357 1.22928 1.03128 0.96484 1.00265 1.20461 0.77414 0.90496 0.78689 1.90696 1.59679 1.23595 1.59868 0.92937 0.92769 0.92983 0.90421 0.89741 0.94285 0.92685 0.90501 0.91554 0.91035 0.90365 0.91486 0.86278 0.86738 0.91571 0.89340 0.94623 0.98913 0.92974 0.95530 0.88899 0.87314 0.88895 0.89558 0.89549 0.89230 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.15 -1.11 9.13 37.16 0.34 -0.77 16 2 C -0.10 -0.90 2.39 -90.62 -0.22 -1.12 16 3 C -0.15 -1.68 7.22 37.15 0.27 -1.41 16 4 C -0.09 -0.82 4.10 -26.74 -0.11 -0.93 16 5 C -0.09 -1.06 0.71 -154.72 -0.11 -1.17 16 6 C -0.14 -1.44 7.28 37.16 0.27 -1.17 16 7 C -0.15 -1.76 6.27 -26.94 -0.17 -1.93 16 8 C 0.52 8.93 7.47 -11.09 -0.08 8.85 16 9 O -0.47 -9.99 16.54 5.60 0.09 -9.90 16 10 N -0.56 -9.50 2.80 -118.62 -0.33 -9.83 16 11 C 0.26 4.89 5.09 -63.32 -0.32 4.57 16 12 N -0.54 -11.66 4.55 -239.97 -1.09 -12.76 16 13 C 0.19 5.10 4.56 -4.98 -0.02 5.08 16 14 C -0.51 -14.36 8.28 -34.44 -0.29 -14.64 16 15 H 0.06 1.80 7.27 -52.49 -0.38 1.41 16 16 C 0.06 1.63 5.46 -17.82 -0.10 1.53 16 17 N -0.89 -26.95 13.29 55.43 0.74 -26.21 16 18 Si 0.89 21.80 29.54 -169.99 -5.02 16.78 16 19 H -0.27 -6.50 7.11 56.52 0.40 -6.10 16 20 H -0.28 -6.86 7.11 56.52 0.40 -6.46 16 21 H -0.28 -6.88 7.11 56.52 0.40 -6.48 16 22 C 0.10 1.92 5.04 -3.82 -0.02 1.90 16 23 C -0.14 -2.34 5.51 -90.62 -0.50 -2.84 16 24 C -0.17 -2.67 10.18 -26.73 -0.27 -2.94 16 25 C -0.19 -2.66 10.25 -26.73 -0.27 -2.94 16 26 C 0.54 8.14 6.45 -11.08 -0.07 8.07 16 27 O -0.47 -7.64 12.00 5.60 0.07 -7.57 16 28 H 0.05 0.41 8.14 -51.93 -0.42 -0.02 16 29 H 0.05 0.39 8.14 -51.93 -0.42 -0.03 16 30 H 0.05 0.34 8.14 -51.93 -0.42 -0.08 16 31 H 0.08 0.64 6.11 -51.93 -0.32 0.33 16 32 H 0.08 1.10 4.47 -51.93 -0.23 0.87 16 33 H 0.04 0.40 8.03 -51.93 -0.42 -0.02 16 34 H 0.05 0.68 7.05 -51.93 -0.37 0.31 16 35 H 0.08 0.53 7.84 -51.93 -0.41 0.13 16 36 H 0.07 0.63 7.82 -51.93 -0.41 0.22 16 37 H 0.07 0.57 8.14 -51.93 -0.42 0.14 16 38 H 0.08 0.81 4.91 -51.93 -0.25 0.55 16 39 H 0.07 0.84 8.14 -51.93 -0.42 0.42 16 40 H 0.12 1.29 7.74 -51.93 -0.40 0.89 16 41 H 0.11 1.00 7.59 -51.93 -0.39 0.60 16 42 H 0.07 1.30 7.38 -51.93 -0.38 0.91 16 43 H 0.09 1.48 6.40 -51.93 -0.33 1.15 16 44 H 0.04 1.03 6.56 -51.93 -0.34 0.69 16 45 H -0.01 -0.20 8.12 -51.93 -0.42 -0.62 16 46 H 0.26 7.43 8.31 -40.82 -0.34 7.09 16 47 H 0.03 0.50 8.05 -51.93 -0.42 0.08 16 48 H 0.09 2.00 6.33 -51.92 -0.33 1.67 16 49 H 0.11 1.90 6.48 -51.93 -0.34 1.56 16 50 H 0.09 1.50 8.03 -51.93 -0.42 1.08 16 51 H 0.09 1.25 8.14 -51.93 -0.42 0.83 16 52 H 0.09 1.07 8.14 -51.92 -0.42 0.65 16 53 H 0.09 1.20 8.14 -51.93 -0.42 0.78 16 LS Contribution 405.04 15.07 6.10 6.10 Total: -1.00 -32.50 405.04 -10.19 -42.69 By element: Atomic # 1 Polarization: 11.64 SS G_CDS: -9.07 Total: 2.56 kcal Atomic # 6 Polarization: -0.19 SS G_CDS: -1.67 Total: -1.87 kcal Atomic # 7 Polarization: -48.11 SS G_CDS: -0.69 Total: -48.80 kcal Atomic # 8 Polarization: -17.63 SS G_CDS: 0.16 Total: -17.47 kcal Atomic # 14 Polarization: 21.80 SS G_CDS: -5.02 Total: 16.78 kcal Total LS contribution 6.10 Total: 6.10 kcal Total: -32.50 -10.19 -42.69 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. ZINC000077164573.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 6.742 kcal (2) G-P(sol) polarization free energy of solvation -32.499 kcal (3) G-ENP(sol) elect.-nuc.-pol. free energy of system -25.756 kcal (4) G-CDS(sol) cavity-dispersion-solvent structure free energy -10.192 kcal (5) G-P-CDS(sol) = G-P(sol) + G-CDS(sol) = (2) + (4) -42.690 kcal (6) G-S(sol) free energy of system = (1) + (5) -35.948 kcal 1SCF WAS SPECIFIED, SO GEOMETRY OPTIMIZATION WAS NOT USED Total computer time = 0.87 seconds