Wall clock time and date at job start Sun Mar 7 2021 06:28:05 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.42903 * 1 3 3 C 1.35891 * 117.00396 * 2 1 4 4 C 1.38606 * 120.03114 * 180.02562 * 3 2 1 5 5 C 1.38662 * 119.89593 * 179.97438 * 4 3 2 6 6 O 1.35891 * 120.03210 * 180.02562 * 5 4 3 7 7 C 1.42904 * 116.99845 * 359.97438 * 6 5 4 8 8 C 1.50697 * 109.47160 * 180.02562 * 7 6 5 9 9 O 1.21279 * 120.00180 * 0.02562 * 8 7 6 10 10 N 1.34775 * 119.99837 * 180.02562 * 8 7 6 11 11 C 1.46506 * 119.99641 * 180.02562 * 10 8 7 12 12 H 1.08995 * 112.85102 * 334.66645 * 11 10 8 13 13 C 1.53780 * 113.61264 * 105.86012 * 11 10 8 14 14 C 1.53957 * 87.01936 * 139.95116 * 13 11 10 15 15 H 1.09001 * 113.54020 * 89.18544 * 14 13 11 16 16 C 1.52998 * 113.74395 * 220.05631 * 14 13 11 17 17 N 1.46898 * 109.47436 * 272.59688 * 16 14 13 18 18 H 1.00906 * 109.46779 * 305.00395 * 17 16 14 19 19 C 1.46900 * 109.47202 * 64.99876 * 17 16 14 20 20 C 1.46895 * 109.47392 * 184.99827 * 17 16 14 21 21 C 1.53781 * 113.61075 * 203.46657 * 11 10 8 22 22 C 1.38753 * 119.94175 * 359.97438 * 5 4 3 23 23 C 1.50700 * 119.97354 * 180.02562 * 22 5 4 24 24 O 1.42896 * 109.47219 * 74.99672 * 23 22 5 25 25 Si 1.86303 * 109.46836 * 195.00095 * 23 22 5 26 26 C 1.38201 * 120.05496 * 0.02562 * 22 5 4 27 27 C 1.38160 * 120.11026 * 359.70189 * 26 22 5 28 28 H 1.09002 * 109.47382 * 179.97438 * 1 2 3 29 29 H 1.09001 * 109.46938 * 59.99709 * 1 2 3 30 30 H 1.09004 * 109.46926 * 299.99830 * 1 2 3 31 31 H 1.07998 * 120.05393 * 359.95736 * 4 3 2 32 32 H 1.08997 * 109.47765 * 300.00979 * 7 6 5 33 33 H 1.09006 * 109.46513 * 60.00515 * 7 6 5 34 34 H 0.96997 * 120.00645 * 359.97438 * 10 8 7 35 35 H 1.08999 * 113.61855 * 254.52845 * 13 11 10 36 36 H 1.09000 * 113.62001 * 25.43079 * 13 11 10 37 37 H 1.09004 * 109.47200 * 32.59651 * 16 14 13 38 38 H 1.08997 * 109.47540 * 152.59460 * 16 14 13 39 39 H 1.09007 * 109.47298 * 300.00280 * 19 17 16 40 40 H 1.09004 * 109.47192 * 59.99962 * 19 17 16 41 41 H 1.09004 * 109.47321 * 180.02562 * 19 17 16 42 42 H 1.08997 * 109.46985 * 60.00427 * 20 17 16 43 43 H 1.08999 * 109.47181 * 179.97438 * 20 17 16 44 44 H 1.08999 * 109.47641 * 299.99821 * 20 17 16 45 45 H 1.09000 * 113.61806 * 334.56857 * 21 11 10 46 46 H 1.08999 * 113.61118 * 105.39044 * 21 11 10 47 47 H 1.08997 * 109.46704 * 314.99398 * 23 22 5 48 48 H 0.96702 * 113.99986 * 60.00564 * 24 23 22 49 49 H 1.48504 * 109.99416 * 179.97438 * 25 23 22 50 50 H 1.48498 * 109.99788 * 301.31160 * 25 23 22 51 51 H 1.08003 * 119.94806 * 179.97438 * 26 22 5 52 52 H 1.07995 * 119.97710 * 180.29266 * 27 26 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 8 1.4290 0.0000 0.0000 3 6 2.0460 1.2108 0.0000 4 6 3.4302 1.2839 -0.0005 5 6 4.0569 2.5209 -0.0011 6 8 5.4137 2.5959 -0.0021 7 6 6.1317 1.3603 -0.0021 8 6 7.6125 1.6399 -0.0027 9 8 8.0136 2.7844 -0.0026 10 7 8.4912 0.6179 -0.0032 11 6 9.9308 0.8897 -0.0033 12 1 10.1783 1.8605 -0.4326 13 6 10.6077 0.6250 1.3519 14 6 11.7847 0.0307 0.5571 15 1 12.5400 0.7668 0.2819 16 6 12.3836 -1.2266 1.1908 17 7 13.4165 -0.8412 2.1616 18 1 14.1008 -0.2536 1.7092 19 6 12.7957 -0.1050 3.2709 20 6 14.0740 -2.0474 2.6818 21 6 10.7758 -0.2658 -0.5650 22 6 3.2973 3.6821 -0.0016 23 6 3.9776 5.0267 -0.0028 24 8 4.5034 5.2934 -1.3044 25 14 2.7442 6.3471 0.4514 26 6 1.9174 3.6065 -0.0011 27 6 1.2906 2.3752 -0.0062 28 1 -0.3634 -1.0277 -0.0005 29 1 -0.3633 0.5139 -0.8900 30 1 -0.3633 0.5138 0.8900 31 1 4.0197 0.3790 -0.0003 32 1 5.8706 0.7883 0.8882 33 1 5.8693 0.7877 -0.8917 34 1 8.1705 -0.2975 -0.0029 35 1 10.8626 1.5334 1.8976 36 1 10.0827 -0.1060 1.9669 37 1 11.5985 -1.7867 1.6987 38 1 12.8289 -1.8482 0.4140 39 1 12.3077 0.7901 2.8849 40 1 12.0564 -0.7397 3.7596 41 1 13.5622 0.1814 3.9910 42 1 14.5347 -2.5937 1.8588 43 1 14.8407 -1.7614 3.4018 44 1 13.3349 -2.6822 3.1705 45 1 10.2969 -1.2409 -0.4752 46 1 11.1650 -0.0760 -1.5653 47 1 4.7902 5.0233 0.7237 48 1 3.8346 5.3115 -2.0027 49 1 3.4030 7.6780 0.4544 50 1 1.6153 6.3408 -0.5133 51 1 1.3278 4.5114 -0.0011 52 1 0.2122 2.3187 -0.0106 There are 66 doubly filled levels RHF calculation, no. of doubly occupied orbitals= 66 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=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 s_22_16608394_4974312.mol2 52 1SCF WAS SPECIFIED, SO GEOMETRY OPTIMIZATION WAS NOT USED SCF FIELD WAS ACHIEVED Wall clock time and date at job start Sun Mar 7 2021 06:28:05 Heat of formation + Delta-G solvation = 54.427669 kcal Electronic energy + Delta-G solvation = -31383.319546 eV Core-core repulsion = 27050.428385 eV Total energy + Delta-G solvation = -4332.891161 eV No. of doubly occupied orbitals = 66 and no. of open levels = 1 Molecular weight (most abundant/longest-lived isotopes) = 352.200 amu Computer time = 0.86 seconds Orbital eigenvalues (eV) -44.55784 -42.60138 -42.09828 -40.96289 -40.22473 -38.88532 -37.52969 -36.01569 -32.93885 -32.64634 -32.03400 -31.92216 -30.26030 -29.09678 -28.56575 -28.05719 -26.68027 -25.87455 -24.30382 -23.13517 -22.88880 -22.63597 -21.49065 -20.66081 -20.29547 -19.80715 -19.42554 -18.79419 -18.66475 -18.36272 -18.09926 -17.98979 -17.69717 -17.38914 -17.08210 -16.95766 -16.58909 -16.52672 -16.30685 -16.13494 -15.92311 -15.89833 -15.78163 -15.61427 -15.56362 -15.30707 -15.19043 -15.05653 -14.74909 -14.32617 -14.20844 -14.15701 -14.03308 -13.76525 -13.59537 -13.31425 -13.25807 -12.86994 -12.59132 -12.32670 -12.11760 -11.88851 -11.40332 -10.98801 -10.34593 -9.71153 -5.43461 -1.48406 -0.41181 -0.30396 0.35565 0.44045 0.50572 0.54085 0.70346 0.82831 1.09485 1.24884 1.34316 1.50089 1.50956 1.57936 1.64176 1.68690 1.76625 1.83403 1.91721 1.95819 1.98135 2.08006 2.15714 2.24671 2.27914 2.33017 2.51562 2.85509 2.89134 2.94201 2.98158 2.98889 3.05365 3.07913 3.12436 3.22022 3.28538 3.40612 3.44659 3.56151 3.60473 3.63244 3.73730 3.78944 3.98050 4.11393 4.13355 4.41467 4.59800 4.68046 4.81096 5.14347 5.46592 5.70082 6.01818 6.03954 Molecular weight = 352.20amu Principal moments of inertia in cm(-1) A = 0.012936 B = 0.002146 C = 0.001934 Principal moments of inertia in units of 10**(-40)*gram-cm**2 A = 2164.063605 B =13045.986829 C =14475.429272 Net atomic charges, atomic populations, and dipole contributions calculated with CM2 Atom NO. Type Charge No. of electrons 1 C 0.024 3.976 2 O -0.315 6.315 3 C 0.175 3.825 4 C -0.248 4.248 5 C 0.189 3.811 6 O -0.273 6.273 7 C 0.036 3.964 8 C 0.519 3.481 9 O -0.500 6.500 10 N -0.724 5.724 11 C 0.132 3.868 12 H 0.124 0.876 13 C -0.161 4.161 14 C -0.143 4.143 15 H 0.118 0.882 16 C 0.003 3.997 17 N -0.386 5.386 18 H 0.427 0.573 19 C -0.056 4.056 20 C -0.048 4.048 21 C -0.132 4.132 22 C -0.180 4.180 23 C -0.057 4.057 24 O -0.523 6.523 25 Si 0.683 3.317 26 C -0.041 4.041 27 C -0.236 4.236 28 H 0.104 0.896 29 H 0.060 0.940 30 H 0.061 0.939 31 H 0.136 0.864 32 H 0.097 0.903 33 H 0.097 0.903 34 H 0.405 0.595 35 H 0.088 0.912 36 H 0.097 0.903 37 H 0.130 0.870 38 H 0.133 0.867 39 H 0.127 0.873 40 H 0.121 0.879 41 H 0.124 0.876 42 H 0.123 0.877 43 H 0.125 0.875 44 H 0.121 0.879 45 H 0.099 0.901 46 H 0.094 0.906 47 H 0.116 0.884 48 H 0.371 0.629 49 H -0.256 1.256 50 H -0.259 1.259 51 H 0.140 0.860 52 H 0.136 0.864 For the charges calculated by CM2: Dipole moment (debyes) X Y Z Total from point charges 28.081 -18.325 9.950 34.976 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.070 4.070 2 O -0.228 6.228 3 C 0.129 3.871 4 C -0.267 4.267 5 C 0.143 3.857 6 O -0.185 6.185 7 C -0.043 4.043 8 C 0.304 3.696 9 O -0.372 6.372 10 N -0.378 5.378 11 C 0.026 3.974 12 H 0.142 0.858 13 C -0.200 4.200 14 C -0.163 4.163 15 H 0.136 0.864 16 C -0.119 4.119 17 N 0.038 4.962 18 H 0.263 0.737 19 C -0.198 4.198 20 C -0.190 4.190 21 C -0.170 4.170 22 C -0.181 4.181 23 C -0.167 4.167 24 O -0.329 6.329 25 Si 0.585 3.415 26 C -0.059 4.059 27 C -0.254 4.254 28 H 0.123 0.877 29 H 0.078 0.922 30 H 0.080 0.920 31 H 0.154 0.846 32 H 0.114 0.886 33 H 0.115 0.885 34 H 0.240 0.760 35 H 0.106 0.894 36 H 0.116 0.884 37 H 0.148 0.852 38 H 0.151 0.849 39 H 0.145 0.855 40 H 0.140 0.860 41 H 0.142 0.858 42 H 0.142 0.858 43 H 0.143 0.857 44 H 0.140 0.860 45 H 0.118 0.882 46 H 0.112 0.888 47 H 0.134 0.866 48 H 0.218 0.782 49 H -0.181 1.181 50 H -0.184 1.184 51 H 0.157 0.843 52 H 0.154 0.846 Dipole moment (debyes) X Y Z Total from point charges 28.517 -16.886 10.162 34.664 hybrid contribution -0.194 -0.832 -1.337 1.586 sum 28.323 -17.719 8.826 34.555 Atomic orbital electron populations 1.23183 0.77284 1.04624 1.01949 1.86062 1.22268 1.26582 1.87896 1.18700 0.91458 0.84499 0.92451 1.20330 0.92954 0.97085 1.16337 1.18339 0.81579 0.92497 0.93330 1.86252 1.15067 1.29040 1.88107 1.22469 0.89770 0.87644 1.04372 1.20163 0.87663 0.85500 0.76234 1.90832 1.79689 1.20790 1.45871 1.45722 1.05518 1.09429 1.77165 1.22323 0.79172 1.00420 0.95469 0.85808 1.23943 0.95345 1.02325 0.98388 1.22966 0.98123 0.95618 0.99550 0.86398 1.24063 0.92760 1.01235 0.93821 1.45081 1.17274 1.17183 1.16615 0.73749 1.24325 1.03876 1.01327 0.90247 1.24140 1.02450 0.86965 1.05422 1.23360 0.95030 0.99914 0.98739 1.19514 0.92746 0.94347 1.11522 1.24779 1.05841 0.97762 0.88329 1.86673 1.37090 1.92835 1.16269 0.95171 0.77692 0.78246 0.90394 1.20838 0.92713 0.98036 0.94303 1.21187 1.00684 0.91745 1.11820 0.87723 0.92157 0.92027 0.84585 0.88564 0.88502 0.76016 0.89396 0.88428 0.85165 0.84861 0.85482 0.86015 0.85768 0.85820 0.85681 0.86014 0.88248 0.88780 0.86587 0.78175 1.18063 1.18439 0.84265 0.84575 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.02 -0.22 9.81 101.05 0.99 0.77 16 2 O -0.32 2.01 10.49 -19.40 -0.20 1.80 16 3 C 0.18 -1.10 6.68 -39.24 -0.26 -1.36 16 4 C -0.25 1.77 8.97 -39.29 -0.35 1.42 16 5 C 0.19 -1.16 6.66 -39.23 -0.26 -1.42 16 6 O -0.27 1.68 8.99 -39.31 -0.35 1.33 16 7 C 0.04 -0.41 5.22 36.00 0.19 -0.22 16 8 C 0.52 -5.69 7.91 -10.99 -0.09 -5.78 16 9 O -0.50 2.79 16.12 -14.35 -0.23 2.56 16 10 N -0.72 11.94 5.36 -53.02 -0.28 11.65 16 11 C 0.13 -2.47 4.21 -67.12 -0.28 -2.75 16 12 H 0.12 -1.93 7.50 -51.93 -0.39 -2.32 16 13 C -0.16 3.75 5.70 -25.83 -0.15 3.60 16 14 C -0.14 3.83 3.30 -89.71 -0.30 3.54 16 15 H 0.12 -3.22 8.09 -51.93 -0.42 -3.64 16 16 C 0.00 -0.08 4.57 -3.83 -0.02 -0.10 16 17 N -0.39 12.69 0.66 -251.26 -0.17 12.52 16 18 H 0.43 -14.72 8.41 -39.29 -0.33 -15.05 16 19 C -0.06 1.74 8.57 60.06 0.51 2.25 16 20 C -0.05 1.53 9.19 60.06 0.55 2.09 16 21 C -0.13 2.99 7.39 -25.92 -0.19 2.80 16 22 C -0.18 1.06 5.27 -104.45 -0.55 0.51 16 23 C -0.06 0.31 2.37 36.01 0.09 0.39 16 24 O -0.52 1.88 12.88 -35.23 -0.45 1.43 16 25 Si 0.68 -4.43 30.24 -169.99 -5.14 -9.57 16 26 C -0.04 0.28 8.65 -39.62 -0.34 -0.06 16 27 C -0.24 1.62 9.04 -39.40 -0.36 1.26 16 28 H 0.10 -1.08 8.14 -51.93 -0.42 -1.50 16 29 H 0.06 -0.57 7.65 -51.93 -0.40 -0.97 16 30 H 0.06 -0.60 7.67 -51.93 -0.40 -1.00 16 31 H 0.14 -1.31 6.30 -52.49 -0.33 -1.65 16 32 H 0.10 -1.33 7.66 -51.93 -0.40 -1.73 16 33 H 0.10 -1.25 7.66 -51.93 -0.40 -1.65 16 34 H 0.40 -7.88 8.71 -40.82 -0.36 -8.24 16 35 H 0.09 -1.96 7.19 -51.93 -0.37 -2.33 16 36 H 0.10 -2.35 7.91 -51.93 -0.41 -2.76 16 37 H 0.13 -4.02 7.80 -51.93 -0.41 -4.43 16 38 H 0.13 -4.15 8.09 -51.93 -0.42 -4.57 16 39 H 0.13 -3.79 5.13 -51.93 -0.27 -4.05 16 40 H 0.12 -3.75 8.11 -51.93 -0.42 -4.17 16 41 H 0.12 -3.86 8.14 -51.93 -0.42 -4.28 16 42 H 0.12 -3.93 8.11 -51.93 -0.42 -4.36 16 43 H 0.12 -3.96 8.14 -51.93 -0.42 -4.38 16 44 H 0.12 -3.83 8.14 -51.93 -0.42 -4.25 16 45 H 0.10 -2.35 8.07 -51.93 -0.42 -2.76 16 46 H 0.09 -1.98 8.14 -51.93 -0.42 -2.41 16 47 H 0.12 -0.59 7.82 -51.93 -0.41 -0.99 16 48 H 0.37 -2.47 8.70 45.56 0.40 -2.07 16 49 H -0.26 0.73 7.11 56.52 0.40 1.13 16 50 H -0.26 1.32 6.68 56.52 0.38 1.70 16 51 H 0.14 -1.12 4.85 -52.48 -0.25 -1.38 16 52 H 0.14 -1.17 6.29 -52.49 -0.33 -1.50 16 LS Contribution 410.52 15.07 6.19 6.19 Total: 1.00 -40.83 410.52 -9.95 -50.78 By element: Atomic # 1 Polarization: -77.14 SS G_CDS: -8.48 Total: -85.62 kcal Atomic # 6 Polarization: 7.75 SS G_CDS: -0.82 Total: 6.93 kcal Atomic # 7 Polarization: 24.63 SS G_CDS: -0.45 Total: 24.18 kcal Atomic # 8 Polarization: 8.36 SS G_CDS: -1.24 Total: 7.12 kcal Atomic # 14 Polarization: -4.43 SS G_CDS: -5.14 Total: -9.57 kcal Total LS contribution 6.19 Total: 6.19 kcal Total: -40.83 -9.95 -50.78 kcal The number of atoms in the molecule is 52 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_16608394_4974312.mol2 52 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 105.208 kcal (2) G-P(sol) polarization free energy of solvation -40.834 kcal (3) G-ENP(sol) elect.-nuc.-pol. free energy of system 64.374 kcal (4) G-CDS(sol) cavity-dispersion-solvent structure free energy -9.947 kcal (5) G-P-CDS(sol) = G-P(sol) + G-CDS(sol) = (2) + (4) -50.780 kcal (6) G-S(sol) free energy of system = (1) + (5) 54.428 kcal 1SCF WAS SPECIFIED, SO GEOMETRY OPTIMIZATION WAS NOT USED Total computer time = 0.86 seconds