Wall clock time and date at job start Sat Apr 11 2020 06:08:24 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.42902 * 1 3 3 C 1.42903 * 113.99413 * 2 1 4 4 C 1.52998 * 109.46957 * 304.38064 * 3 2 1 5 5 N 1.46500 * 109.47551 * 300.57108 * 4 3 2 6 6 C 1.34783 * 120.00047 * 185.27765 * 5 4 3 7 7 O 1.21342 * 119.99856 * 0.02562 * 6 5 4 8 8 C 1.49450 * 119.99680 * 179.97438 * 6 5 4 9 9 O 1.21341 * 119.99998 * 0.02562 * 8 6 5 10 10 N 1.34773 * 119.99667 * 179.97438 * 8 6 5 11 11 C 1.46496 * 120.00197 * 180.02562 * 10 8 6 12 12 C 1.52998 * 109.47536 * 180.02562 * 11 10 8 13 13 C 1.52999 * 109.46967 * 180.02562 * 12 11 10 14 14 C 1.50709 * 109.46953 * 180.02562 * 13 12 11 15 15 H 1.08003 * 119.99766 * 359.97438 * 14 13 12 16 16 C 1.37876 * 120.00174 * 179.97438 * 14 13 12 17 17 N 1.31515 * 119.99552 * 359.97438 * 16 14 13 18 18 Si 1.86796 * 119.99914 * 180.02562 * 16 14 13 19 19 H 1.48491 * 109.47158 * 359.97438 * 18 16 14 20 20 H 1.48504 * 109.47158 * 120.00378 * 18 16 14 21 21 H 1.48507 * 109.46811 * 239.99525 * 18 16 14 22 22 C 1.52999 * 109.47304 * 64.38548 * 3 2 1 23 23 C 1.52998 * 109.47035 * 60.00106 * 22 3 2 24 24 C 1.52998 * 109.47401 * 59.99886 * 23 22 3 25 25 C 1.53004 * 109.47037 * 59.99811 * 24 23 22 26 26 C 1.52999 * 109.47119 * 300.00063 * 24 23 22 27 27 C 1.53001 * 109.46874 * 184.38895 * 3 2 1 28 28 H 1.09006 * 109.46997 * 184.79057 * 1 2 3 29 29 H 1.08991 * 109.47596 * 304.79405 * 1 2 3 30 30 H 1.09007 * 109.47058 * 64.79821 * 1 2 3 31 31 H 1.09001 * 109.47087 * 60.57151 * 4 3 2 32 32 H 1.08999 * 109.46962 * 180.56660 * 4 3 2 33 33 H 0.97001 * 119.99704 * 5.27310 * 5 4 3 34 34 H 0.97006 * 120.00204 * 0.02562 * 10 8 6 35 35 H 1.09004 * 109.47166 * 300.00999 * 11 10 8 36 36 H 1.09002 * 109.46919 * 60.00136 * 11 10 8 37 37 H 1.08995 * 109.47150 * 300.00325 * 12 11 10 38 38 H 1.09002 * 109.46915 * 60.00316 * 12 11 10 39 39 H 1.09001 * 109.47427 * 299.99322 * 13 12 11 40 40 H 1.08992 * 109.47744 * 60.00403 * 13 12 11 41 41 H 0.97010 * 119.99683 * 179.97438 * 17 16 14 42 42 H 1.09002 * 109.46720 * 179.97438 * 22 3 2 43 43 H 1.08998 * 109.46885 * 299.99358 * 22 3 2 44 44 H 1.08999 * 109.46613 * 180.02562 * 23 22 3 45 45 H 1.09002 * 109.47170 * 299.99484 * 23 22 3 46 46 H 1.08996 * 109.47026 * 179.97438 * 24 23 22 47 47 H 1.08996 * 109.47123 * 299.99830 * 25 24 23 48 48 H 1.08996 * 109.47122 * 60.00279 * 25 24 23 49 49 H 1.09000 * 109.46763 * 180.02562 * 25 24 23 50 50 H 1.09002 * 109.46720 * 299.99954 * 26 24 23 51 51 H 1.09000 * 109.47332 * 180.02562 * 26 24 23 52 52 H 1.08999 * 109.47680 * 59.99366 * 27 3 2 53 53 H 1.09002 * 109.47146 * 180.02562 * 27 3 2 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.0101 1.3055 0.0000 4 6 1.4733 2.1027 -1.1905 5 7 1.8085 1.4067 -2.4353 6 6 1.5286 1.9784 -3.6234 7 8 0.9994 3.0697 -3.6619 8 6 1.8710 1.2686 -4.8932 9 8 2.4006 0.1775 -4.8546 10 7 1.5907 1.8400 -6.0812 11 6 1.9258 1.1441 -7.3260 12 6 1.5061 2.0039 -8.5199 13 6 1.8556 1.2767 -9.8199 14 6 1.4427 2.1239 -10.9960 15 1 0.9772 3.0846 -10.8325 16 6 1.6592 1.6724 -12.2806 17 7 2.2265 0.5027 -12.4797 18 14 1.1481 2.7228 -13.7383 19 1 0.5371 3.9847 -13.2490 20 1 2.3414 3.0377 -14.5642 21 1 0.1590 1.9806 -14.5606 22 6 1.6478 2.0251 1.3007 23 6 2.1847 1.2279 2.4911 24 6 3.7061 1.1089 2.3808 25 6 4.3283 2.5067 2.3807 26 6 4.0684 0.3894 1.0801 27 6 3.5315 1.1865 -0.1104 28 1 -0.3633 -1.0241 0.0858 29 1 -0.3634 0.5863 0.8438 30 1 -0.3633 0.4376 -0.9299 31 1 0.3906 2.1961 -1.1067 32 1 1.9242 3.0950 -1.1962 33 1 2.2319 0.5346 -2.4045 34 1 1.1670 2.7121 -6.1121 35 1 3.0008 0.9677 -7.3656 36 1 1.3997 0.1902 -7.3619 37 1 0.4313 2.1805 -8.4805 38 1 2.0326 2.9576 -8.4838 39 1 2.9304 1.0998 -9.8597 40 1 1.3288 0.3232 -9.8559 41 1 2.3786 0.1849 -13.3835 42 1 2.0915 3.0207 1.3008 43 1 0.5640 2.1099 1.3792 44 1 1.9262 1.7404 3.4177 45 1 1.7414 0.2321 2.4912 46 1 4.0885 0.5414 3.2291 47 1 3.9459 3.0746 1.5326 48 1 4.0702 3.0194 3.3074 49 1 5.4121 2.4219 2.3017 50 1 3.6251 -0.6064 1.0802 51 1 5.1522 0.3042 1.0017 52 1 3.7897 0.6739 -1.0370 53 1 3.9748 2.1823 -0.1099 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 ZINC000343333276.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:08:24 Heat of formation + Delta-G solvation = -148.100092 kcal Electronic energy + Delta-G solvation = -30025.430239 eV Core-core repulsion = 25896.045878 eV Total energy + Delta-G solvation = -4129.384361 eV No. of doubly occupied orbitals = 66 Molecular weight (most abundant/longest-lived isotopes) = 340.220 amu Computer time = 1.19 seconds Orbital eigenvalues (eV) -41.13196 -40.24898 -38.44378 -37.09355 -35.78852 -34.86716 -33.32106 -31.81675 -31.38621 -31.24465 -28.41911 -27.27309 -27.09478 -25.84418 -24.16043 -23.22103 -22.37265 -21.17562 -20.77096 -20.12827 -19.20041 -18.47469 -16.95383 -16.65878 -16.37252 -16.24160 -16.09888 -15.83810 -15.64926 -15.33845 -14.82249 -14.67502 -14.57527 -14.43110 -14.11093 -13.71665 -13.68404 -13.50558 -13.11861 -12.84076 -12.73417 -12.54739 -12.32116 -12.16510 -11.88878 -11.79439 -11.74774 -11.68855 -11.65897 -11.35775 -10.93290 -10.79140 -10.76452 -10.73537 -10.58405 -10.38267 -9.97718 -9.90821 -9.75412 -9.49789 -9.37133 -9.12049 -8.97578 -8.38574 -6.03106 -4.06026 1.41963 2.30749 3.22119 3.29198 3.29893 3.37516 3.41443 3.69914 4.01469 4.10919 4.26446 4.43498 4.51468 4.55392 4.57597 4.67953 4.86391 4.90495 4.99071 5.01637 5.07084 5.10471 5.14273 5.27022 5.31221 5.32987 5.37585 5.41276 5.46170 5.56629 5.59490 5.60409 5.69200 5.74080 5.81314 5.86937 5.92406 5.95401 6.08091 6.08184 6.27428 6.53944 6.69131 6.76056 6.91969 6.94995 7.40910 7.61715 7.71942 7.73156 7.95386 8.05637 8.36624 9.00688 9.58398 11.08173 Molecular weight = 340.22amu Principal moments of inertia in cm(-1) A = 0.029214 B = 0.001627 C = 0.001620 Principal moments of inertia in units of 10**(-40)*gram-cm**2 A = 958.199152 B =17206.876242 C =17282.537122 Net atomic charges, atomic populations, and dipole contributions calculated with CM2 Atom NO. Type Charge No. of electrons 1 C 0.021 3.979 2 O -0.367 6.367 3 C 0.109 3.891 4 C 0.091 3.909 5 N -0.696 5.696 6 C 0.508 3.492 7 O -0.522 6.522 8 C 0.500 3.500 9 O -0.524 6.524 10 N -0.692 5.692 11 C 0.114 3.886 12 C -0.123 4.123 13 C 0.024 3.976 14 C -0.522 4.522 15 H 0.057 0.943 16 C 0.056 3.944 17 N -0.894 5.894 18 Si 0.893 3.107 19 H -0.273 1.273 20 H -0.284 1.284 21 H -0.284 1.284 22 C -0.142 4.142 23 C -0.107 4.107 24 C -0.092 4.092 25 C -0.153 4.153 26 C -0.111 4.111 27 C -0.108 4.108 28 H 0.087 0.913 29 H 0.045 0.955 30 H 0.044 0.956 31 H 0.075 0.925 32 H 0.084 0.916 33 H 0.417 0.583 34 H 0.411 0.589 35 H 0.071 0.929 36 H 0.071 0.929 37 H 0.057 0.943 38 H 0.057 0.943 39 H 0.020 0.980 40 H 0.020 0.980 41 H 0.261 0.739 42 H 0.079 0.921 43 H 0.069 0.931 44 H 0.062 0.938 45 H 0.069 0.931 46 H 0.069 0.931 47 H 0.059 0.941 48 H 0.052 0.948 49 H 0.054 0.946 50 H 0.073 0.927 51 H 0.063 0.937 52 H 0.073 0.927 53 H 0.081 0.919 For the charges calculated by CM2: Dipole moment (debyes) X Y Z Total from point charges 0.262 1.783 34.652 34.699 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.075 4.075 2 O -0.286 6.286 3 C 0.069 3.931 4 C -0.033 4.033 5 N -0.351 5.351 6 C 0.288 3.712 7 O -0.398 6.398 8 C 0.279 3.721 9 O -0.401 6.401 10 N -0.343 5.343 11 C -0.008 4.008 12 C -0.161 4.161 13 C -0.014 4.014 14 C -0.560 4.560 15 H 0.075 0.925 16 C -0.145 4.145 17 N -0.533 5.533 18 Si 0.721 3.279 19 H -0.197 1.197 20 H -0.210 1.210 21 H -0.210 1.210 22 C -0.180 4.180 23 C -0.145 4.145 24 C -0.111 4.111 25 C -0.210 4.210 26 C -0.148 4.148 27 C -0.146 4.146 28 H 0.105 0.895 29 H 0.064 0.936 30 H 0.063 0.937 31 H 0.093 0.907 32 H 0.102 0.898 33 H 0.255 0.745 34 H 0.247 0.753 35 H 0.089 0.911 36 H 0.089 0.911 37 H 0.076 0.924 38 H 0.076 0.924 39 H 0.038 0.962 40 H 0.038 0.962 41 H 0.070 0.930 42 H 0.098 0.902 43 H 0.087 0.913 44 H 0.080 0.920 45 H 0.088 0.912 46 H 0.088 0.912 47 H 0.078 0.922 48 H 0.071 0.929 49 H 0.073 0.927 50 H 0.091 0.909 51 H 0.082 0.918 52 H 0.091 0.909 53 H 0.100 0.900 Dipole moment (debyes) X Y Z Total from point charges 0.430 1.497 34.489 34.525 hybrid contribution -0.559 0.775 -0.923 1.328 sum -0.130 2.271 33.567 33.644 Atomic orbital electron populations 1.22841 0.81340 1.02981 1.00324 1.88001 1.19555 1.26272 1.94767 1.21757 0.94681 0.83814 0.92799 1.22092 1.01313 0.96921 0.82979 1.45645 1.61609 1.23491 1.04363 1.20377 0.79556 0.85910 0.85320 1.90756 1.41419 1.20106 1.87534 1.19980 0.80797 0.86201 0.85080 1.90766 1.41970 1.20037 1.87336 1.46067 1.59753 1.23086 1.05428 1.21845 1.01146 0.96512 0.81330 1.21926 1.00697 0.98789 0.94646 1.19042 0.97276 0.95707 0.89375 1.21243 1.39844 1.04079 0.90881 0.92515 1.25030 0.95219 0.95120 0.99121 1.69979 1.43170 1.18805 1.21312 0.86516 0.78272 0.80478 0.82608 1.19739 1.20996 1.21003 1.22077 1.01500 1.00041 0.94407 1.21441 0.95370 1.00719 0.96944 1.21102 0.95322 0.96600 0.98050 1.21965 1.00601 0.96667 1.01785 1.21491 1.00438 0.99577 0.93311 1.21608 0.91382 1.01849 0.99720 0.89454 0.93616 0.93714 0.90699 0.89771 0.74491 0.75344 0.91084 0.91092 0.92444 0.92429 0.96172 0.96173 0.92981 0.90228 0.91287 0.91967 0.91240 0.91223 0.92170 0.92881 0.92722 0.90867 0.91828 0.90869 0.90044 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.02 0.14 10.13 101.05 1.02 1.16 16 2 O -0.37 -2.87 8.65 -35.23 -0.30 -3.17 16 3 C 0.11 0.73 0.36 -90.62 -0.03 0.70 16 4 C 0.09 0.74 4.06 -4.04 -0.02 0.72 16 5 N -0.70 -7.19 5.08 -60.73 -0.31 -7.50 16 6 C 0.51 6.68 7.89 -11.58 -0.09 6.58 16 7 O -0.52 -7.75 16.47 5.50 0.09 -7.66 16 8 C 0.50 7.54 7.89 -11.59 -0.09 7.45 16 9 O -0.52 -9.01 16.48 5.50 0.09 -8.92 16 10 N -0.69 -10.44 5.55 -60.73 -0.34 -10.78 16 11 C 0.11 1.96 5.67 -4.04 -0.02 1.94 16 12 C -0.12 -2.39 5.27 -26.73 -0.14 -2.53 16 13 C 0.02 0.59 4.97 -27.88 -0.14 0.46 16 14 C -0.52 -14.56 9.11 -34.44 -0.31 -14.87 16 15 H 0.06 1.53 7.74 -52.48 -0.41 1.13 16 16 C 0.06 1.61 5.46 -17.82 -0.10 1.51 16 17 N -0.89 -26.68 13.29 55.43 0.74 -25.94 16 18 Si 0.89 21.51 29.54 -169.99 -5.02 16.49 16 19 H -0.27 -6.42 7.11 56.52 0.40 -6.02 16 20 H -0.28 -6.77 7.11 56.52 0.40 -6.37 16 21 H -0.28 -6.77 7.11 56.52 0.40 -6.37 16 22 C -0.14 -0.76 4.21 -26.73 -0.11 -0.87 16 23 C -0.11 -0.56 5.30 -26.73 -0.14 -0.70 16 24 C -0.09 -0.49 3.20 -90.62 -0.29 -0.78 16 25 C -0.15 -0.75 8.50 37.16 0.32 -0.43 16 26 C -0.11 -0.69 5.31 -26.73 -0.14 -0.83 16 27 C -0.11 -0.72 4.51 -26.73 -0.12 -0.84 16 28 H 0.09 0.50 8.14 -51.93 -0.42 0.08 16 29 H 0.05 0.24 6.32 -51.93 -0.33 -0.09 16 30 H 0.04 0.32 6.64 -51.93 -0.34 -0.02 16 31 H 0.07 0.61 6.63 -51.93 -0.34 0.27 16 32 H 0.08 0.68 8.13 -51.93 -0.42 0.26 16 33 H 0.42 4.29 6.80 -40.82 -0.28 4.01 16 34 H 0.41 5.66 7.95 -40.82 -0.32 5.34 16 35 H 0.07 1.26 8.14 -51.93 -0.42 0.84 16 36 H 0.07 1.26 8.14 -51.93 -0.42 0.84 16 37 H 0.06 1.08 8.14 -51.93 -0.42 0.66 16 38 H 0.06 1.09 8.14 -51.93 -0.42 0.66 16 39 H 0.02 0.53 8.14 -51.93 -0.42 0.11 16 40 H 0.02 0.53 8.14 -51.93 -0.42 0.11 16 41 H 0.26 7.32 8.31 -40.82 -0.34 6.98 16 42 H 0.08 0.39 6.43 -51.93 -0.33 0.05 16 43 H 0.07 0.34 6.53 -51.93 -0.34 0.00 16 44 H 0.06 0.29 8.14 -51.93 -0.42 -0.14 16 45 H 0.07 0.42 7.68 -51.93 -0.40 0.02 16 46 H 0.07 0.36 8.14 -51.93 -0.42 -0.06 16 47 H 0.06 0.30 4.81 -51.93 -0.25 0.05 16 48 H 0.05 0.24 8.14 -51.93 -0.42 -0.18 16 49 H 0.05 0.27 8.14 -51.93 -0.42 -0.16 16 50 H 0.07 0.53 7.68 -51.93 -0.40 0.13 16 51 H 0.06 0.37 8.14 -51.93 -0.42 -0.05 16 52 H 0.07 0.56 7.01 -51.93 -0.36 0.20 16 53 H 0.08 0.49 6.43 -51.93 -0.33 0.16 16 LS Contribution 411.00 15.07 6.19 6.19 Total: -1.00 -31.83 411.00 -8.35 -40.18 By element: Atomic # 1 Polarization: 11.52 SS G_CDS: -9.07 Total: 2.44 kcal Atomic # 6 Polarization: -0.92 SS G_CDS: -0.41 Total: -1.33 kcal Atomic # 7 Polarization: -44.31 SS G_CDS: 0.09 Total: -44.22 kcal Atomic # 8 Polarization: -19.63 SS G_CDS: -0.12 Total: -19.75 kcal Atomic # 14 Polarization: 21.51 SS G_CDS: -5.02 Total: 16.49 kcal Total LS contribution 6.19 Total: 6.19 kcal Total: -31.83 -8.35 -40.18 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. ZINC000343333276.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 -107.922 kcal (2) G-P(sol) polarization free energy of solvation -31.833 kcal (3) G-ENP(sol) elect.-nuc.-pol. free energy of system -139.755 kcal (4) G-CDS(sol) cavity-dispersion-solvent structure free energy -8.345 kcal (5) G-P-CDS(sol) = G-P(sol) + G-CDS(sol) = (2) + (4) -40.178 kcal (6) G-S(sol) free energy of system = (1) + (5) -148.100 kcal 1SCF WAS SPECIFIED, SO GEOMETRY OPTIMIZATION WAS NOT USED Total computer time = 1.19 seconds