Wall clock time and date at job start Sat Apr 11 2020 06:13:02 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.52994 * 1 3 3 N 1.46903 * 109.47096 * 2 1 4 4 C 1.46897 * 110.99716 * 66.05086 * 3 2 1 5 5 C 1.46901 * 111.00257 * 190.00345 * 3 2 1 6 6 C 1.52994 * 109.47401 * 170.00295 * 5 3 2 7 7 C 1.50703 * 109.47095 * 180.02562 * 6 5 3 8 8 O 1.21277 * 120.00298 * 359.97438 * 7 6 5 9 9 N 1.34780 * 119.99785 * 180.02562 * 7 6 5 10 10 C 1.47428 * 125.64780 * 359.96987 * 9 7 6 11 11 C 1.54913 * 104.83199 * 204.16650 * 10 9 7 12 12 C 1.55157 * 101.58039 * 322.99121 * 11 10 9 13 13 H 1.08994 * 110.72470 * 276.96160 * 12 11 10 14 14 N 1.46898 * 110.71938 * 153.86222 * 12 11 10 15 15 C 1.46906 * 110.99888 * 86.45374 * 14 12 11 16 16 C 1.50699 * 109.47140 * 180.02562 * 15 14 12 17 17 O 1.21297 * 119.99980 * 359.97438 * 16 15 14 18 18 N 1.34779 * 120.00037 * 180.02562 * 16 15 14 19 19 C 1.37097 * 119.99944 * 180.02562 * 18 16 15 20 20 H 1.08003 * 120.00055 * 359.97438 * 19 18 16 21 21 C 1.38948 * 119.99682 * 179.97438 * 19 18 16 22 22 N 1.31408 * 120.00492 * 359.97438 * 21 19 18 23 23 Si 1.86803 * 119.99897 * 180.02562 * 21 19 18 24 24 H 1.48501 * 109.46948 * 359.97438 * 23 21 19 25 25 H 1.48497 * 109.47378 * 120.00054 * 23 21 19 26 26 H 1.48501 * 109.47174 * 239.99969 * 23 21 19 27 27 C 1.47021 * 125.64787 * 180.02562 * 9 7 6 28 28 H 1.08996 * 109.88469 * 59.48561 * 27 9 7 29 29 C 1.53003 * 109.88134 * 298.31903 * 27 9 7 30 30 H 1.08997 * 109.47432 * 59.99935 * 1 2 3 31 31 H 1.09005 * 109.47311 * 180.02562 * 1 2 3 32 32 H 1.09003 * 109.47172 * 299.99495 * 1 2 3 33 33 H 1.09003 * 109.47172 * 240.00505 * 2 1 3 34 34 H 1.08997 * 109.47432 * 120.00065 * 2 1 3 35 35 H 1.09008 * 109.47068 * 183.95224 * 4 3 2 36 36 H 1.08997 * 109.47061 * 303.95311 * 4 3 2 37 37 H 1.09000 * 109.47329 * 63.95436 * 4 3 2 38 38 H 1.08999 * 109.47077 * 49.99747 * 5 3 2 39 39 H 1.09003 * 109.46974 * 290.00439 * 5 3 2 40 40 H 1.08999 * 109.47453 * 59.99682 * 6 5 3 41 41 H 1.09003 * 109.47174 * 300.00021 * 6 5 3 42 42 H 1.09002 * 110.36632 * 85.33078 * 10 9 7 43 43 H 1.08993 * 110.36844 * 323.00320 * 10 9 7 44 44 H 1.08997 * 111.00332 * 204.91969 * 11 10 9 45 45 H 1.09000 * 111.00128 * 81.06531 * 11 10 9 46 46 H 1.00902 * 111.00533 * 322.49780 * 14 12 11 47 47 H 1.09000 * 109.47033 * 300.00113 * 15 14 12 48 48 H 1.09001 * 109.46544 * 60.00058 * 15 14 12 49 49 H 0.96996 * 119.99844 * 359.97438 * 18 16 15 50 50 H 0.97001 * 120.00876 * 179.97438 * 22 21 19 51 51 H 1.09003 * 109.47086 * 302.19341 * 29 27 9 52 52 H 1.09001 * 109.46964 * 62.18985 * 29 27 9 53 53 H 1.08998 * 109.47155 * 182.19313 * 29 27 9 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.5299 0.0000 0.0000 3 7 2.0196 1.3850 0.0000 4 6 1.6702 2.0668 1.2534 5 6 3.4684 1.4312 -0.2382 6 6 3.8998 2.8792 -0.4791 7 6 5.3860 2.9265 -0.7241 8 8 6.0363 1.9029 -0.7072 9 7 5.9927 4.1064 -0.9610 10 6 5.3265 5.4207 -1.0092 11 6 6.2504 6.3016 -1.8868 12 6 7.6557 5.8030 -1.4578 13 1 7.9784 6.2995 -0.5428 14 7 8.6330 6.0143 -2.5340 15 6 9.2088 7.3640 -2.4646 16 6 10.2040 7.5436 -3.5819 17 8 10.4134 6.6363 -4.3592 18 7 10.8622 8.7120 -3.7168 19 6 11.7672 8.8755 -4.7335 20 1 11.9530 8.0680 -5.4262 21 6 12.4453 10.0803 -4.8729 22 7 12.2188 11.0631 -4.0305 23 14 13.6778 10.3034 -6.2587 24 1 13.7493 9.0572 -7.0630 25 1 13.2485 11.4294 -7.1265 26 1 15.0161 10.6009 -5.6878 27 6 7.4284 4.2960 -1.2151 28 1 7.7267 3.7292 -2.0971 29 6 8.2369 3.8369 0.0000 30 1 -0.3634 0.5138 -0.8899 31 1 -0.3634 -1.0277 0.0005 32 1 -0.3634 0.5138 0.8901 33 1 1.8933 -0.5138 -0.8901 34 1 1.8933 -0.5138 0.8899 35 1 1.9824 3.1098 1.1989 36 1 0.5920 2.0178 1.4053 37 1 2.1770 1.5789 2.0859 38 1 3.7121 0.8291 -1.1135 39 1 3.9925 1.0359 0.6320 40 1 3.6561 3.4814 0.3962 41 1 3.3757 3.2745 -1.3493 42 1 5.2425 5.8405 -0.0068 43 1 4.3417 5.3302 -1.4673 44 1 6.1220 7.3581 -1.6512 45 1 6.0764 6.1153 -2.9466 46 1 8.2181 5.8442 -3.4379 47 1 8.4142 8.1035 -2.5639 48 1 9.7104 7.4961 -1.5059 49 1 10.6950 9.4374 -3.0949 50 1 12.6924 11.9040 -4.1275 51 1 7.9183 4.3952 0.8804 52 1 8.0715 2.7723 0.1653 53 1 9.2969 4.0164 -0.1801 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 ZINC001486093324.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:13:02 Heat of formation + Delta-G solvation = -65.728289 kcal Electronic energy + Delta-G solvation = -29534.698090 eV Core-core repulsion = 25441.178390 eV Total energy + Delta-G solvation = -4093.519700 eV No. of doubly occupied orbitals = 66 Molecular weight (most abundant/longest-lived isotopes) = 340.226 amu Computer time = 0.88 seconds Orbital eigenvalues (eV) -40.53120 -38.75818 -38.16352 -37.55353 -34.76649 -33.81533 -32.88224 -31.47685 -31.26510 -30.71814 -29.51535 -27.16013 -26.86821 -24.99699 -23.96606 -23.47134 -22.41884 -21.25684 -20.75468 -19.35612 -18.93386 -17.91704 -17.46937 -16.79299 -16.20008 -16.10274 -15.81480 -15.45781 -15.09250 -14.76796 -14.56160 -14.53112 -14.12438 -13.94030 -13.86630 -13.67731 -13.41804 -13.20270 -13.02662 -12.70286 -12.51353 -12.51211 -12.28235 -12.24668 -12.14548 -11.98443 -11.77079 -11.64078 -11.57941 -11.48008 -11.28917 -11.23658 -11.13096 -11.02024 -10.83252 -10.35207 -10.23980 -9.90901 -9.53684 -8.82081 -8.77722 -8.49953 -8.10387 -7.97834 -6.43520 -3.84280 2.48188 2.68276 3.15454 3.26387 3.33738 3.51412 3.88163 4.00920 4.20935 4.27671 4.40606 4.48494 4.60002 4.67748 4.77219 4.80112 4.86710 5.00604 5.02308 5.07353 5.08682 5.10122 5.14455 5.17763 5.23845 5.29019 5.38595 5.41351 5.48682 5.50346 5.54681 5.58306 5.62372 5.63247 5.72327 5.78820 5.81096 5.83597 5.94082 6.06871 6.23187 6.25780 6.30082 6.60770 6.62073 7.14184 7.35967 7.54035 7.56253 8.01556 8.04741 8.59200 9.16406 9.52234 10.03516 10.74207 Molecular weight = 340.23amu Principal moments of inertia in cm(-1) A = 0.034433 B = 0.001584 C = 0.001574 Principal moments of inertia in units of 10**(-40)*gram-cm**2 A = 812.967240 B =17676.604181 C =17781.757796 Net atomic charges, atomic populations, and dipole contributions calculated with CM2 Atom NO. Type Charge No. of electrons 1 C -0.149 4.149 2 C 0.050 3.950 3 N -0.545 5.545 4 C 0.017 3.983 5 C 0.070 3.930 6 C -0.132 4.132 7 C 0.530 3.470 8 O -0.539 6.539 9 N -0.616 5.616 10 C 0.103 3.897 11 C -0.170 4.170 12 C 0.084 3.916 13 H 0.095 0.905 14 N -0.657 5.657 15 C 0.053 3.947 16 C 0.440 3.560 17 O -0.582 6.582 18 N -0.563 5.563 19 C -0.299 4.299 20 H 0.104 0.896 21 C 0.023 3.977 22 N -0.902 5.902 23 Si 0.929 3.071 24 H -0.268 1.268 25 H -0.282 1.282 26 H -0.281 1.281 27 C 0.142 3.858 28 H 0.086 0.914 29 C -0.155 4.155 30 H 0.063 0.937 31 H 0.058 0.942 32 H 0.054 0.946 33 H 0.076 0.924 34 H 0.033 0.967 35 H 0.069 0.931 36 H 0.071 0.929 37 H 0.026 0.974 38 H 0.083 0.917 39 H 0.043 0.957 40 H 0.091 0.909 41 H 0.100 0.900 42 H 0.070 0.930 43 H 0.077 0.923 44 H 0.086 0.914 45 H 0.080 0.920 46 H 0.353 0.647 47 H 0.045 0.955 48 H 0.087 0.913 49 H 0.392 0.608 50 H 0.271 0.729 51 H 0.050 0.950 52 H 0.065 0.935 53 H 0.072 0.928 For the charges calculated by CM2: Dipole moment (debyes) X Y Z Total from point charges -23.247 -12.530 12.870 29.378 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.206 4.206 2 C -0.078 4.078 3 N -0.270 5.270 4 C -0.132 4.132 5 C -0.058 4.058 6 C -0.172 4.172 7 C 0.318 3.682 8 O -0.416 6.416 9 N -0.351 5.351 10 C -0.021 4.021 11 C -0.210 4.210 12 C -0.028 4.028 13 H 0.113 0.887 14 N -0.288 5.288 15 C -0.080 4.080 16 C 0.225 3.775 17 O -0.466 6.466 18 N -0.200 5.200 19 C -0.429 4.429 20 H 0.122 0.878 21 C -0.179 4.179 22 N -0.543 5.543 23 Si 0.756 3.244 24 H -0.192 1.192 25 H -0.208 1.208 26 H -0.207 1.207 27 C 0.039 3.961 28 H 0.104 0.896 29 C -0.214 4.214 30 H 0.082 0.918 31 H 0.077 0.923 32 H 0.073 0.927 33 H 0.094 0.906 34 H 0.052 0.948 35 H 0.088 0.912 36 H 0.090 0.910 37 H 0.044 0.956 38 H 0.101 0.899 39 H 0.061 0.939 40 H 0.109 0.891 41 H 0.119 0.881 42 H 0.088 0.912 43 H 0.095 0.905 44 H 0.104 0.896 45 H 0.099 0.901 46 H 0.174 0.826 47 H 0.063 0.937 48 H 0.105 0.895 49 H 0.225 0.775 50 H 0.081 0.919 51 H 0.069 0.931 52 H 0.084 0.916 53 H 0.091 0.909 Dipole moment (debyes) X Y Z Total from point charges -22.397 -13.631 12.313 28.966 hybrid contribution 0.150 1.041 0.139 1.061 sum -22.246 -12.590 12.452 28.434 Atomic orbital electron populations 1.21807 0.94743 1.01615 1.02484 1.22084 0.96117 0.88543 1.01040 1.62000 1.13858 1.16244 1.34852 1.22503 1.01564 0.99330 0.89762 1.21818 0.86496 0.95468 1.01988 1.21605 0.92091 0.99181 1.04368 1.20242 0.89677 0.81917 0.76335 1.90708 1.66839 1.32749 1.51349 1.48623 1.09216 1.06452 1.70789 1.22271 0.96760 0.83013 1.00043 1.23616 0.96644 1.00155 1.00574 1.22692 0.89561 0.94115 0.96444 0.88658 1.59958 1.35154 1.27208 1.06516 1.21682 0.97176 0.89247 0.99884 1.19794 0.85700 0.84094 0.87905 1.90538 1.68337 1.41406 1.46353 1.41860 1.36843 1.12691 1.28596 1.19762 1.16253 0.99826 1.07034 0.87845 1.25253 0.98816 0.95198 0.98618 1.69613 1.43072 1.04630 1.36991 0.85806 0.79700 0.78459 0.80426 1.19196 1.20774 1.20737 1.21503 0.79334 0.96018 0.99229 0.89589 1.22053 1.00549 1.02153 0.96596 0.91768 0.92266 0.92690 0.90601 0.94835 0.91220 0.91047 0.95575 0.89862 0.93895 0.89069 0.88132 0.91173 0.90468 0.89555 0.90117 0.82624 0.93675 0.89492 0.77475 0.91923 0.93112 0.91572 0.90873 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 -0.80 9.84 37.15 0.37 -0.44 16 2 C 0.05 0.33 5.68 -3.83 -0.02 0.31 16 3 N -0.55 -4.10 4.76 -242.92 -1.16 -5.26 16 4 C 0.02 0.10 9.13 60.06 0.55 0.65 16 5 C 0.07 0.61 4.25 -3.83 -0.02 0.60 16 6 C -0.13 -1.06 4.63 -27.88 -0.13 -1.18 16 7 C 0.53 6.05 7.64 -10.98 -0.08 5.97 16 8 O -0.54 -8.17 15.03 5.56 0.08 -8.09 16 9 N -0.62 -6.47 3.06 -163.01 -0.50 -6.97 16 10 C 0.10 0.78 6.69 -2.53 -0.02 0.76 16 11 C -0.17 -1.50 6.63 -24.58 -0.16 -1.67 16 12 C 0.08 0.95 2.76 -65.88 -0.18 0.77 16 13 H 0.10 1.02 7.61 -51.93 -0.40 0.63 16 14 N -0.66 -9.97 5.91 -106.05 -0.63 -10.60 16 15 C 0.05 0.87 6.02 -4.97 -0.03 0.84 16 16 C 0.44 10.12 7.82 -10.99 -0.09 10.03 16 17 O -0.58 -15.63 15.78 5.54 0.09 -15.54 16 18 N -0.56 -14.00 5.29 -10.96 -0.06 -14.05 16 19 C -0.30 -8.27 9.68 -14.93 -0.14 -8.42 16 20 H 0.10 3.02 7.16 -52.48 -0.38 2.64 16 21 C 0.02 0.60 5.50 -17.47 -0.10 0.51 16 22 N -0.90 -24.23 13.05 55.50 0.72 -23.50 16 23 Si 0.93 20.91 29.55 -169.99 -5.02 15.89 16 24 H -0.27 -6.05 7.11 56.52 0.40 -5.65 16 25 H -0.28 -6.23 7.11 56.52 0.40 -5.83 16 26 H -0.28 -6.25 7.11 56.52 0.40 -5.85 16 27 C 0.14 1.74 3.21 -66.95 -0.22 1.53 16 28 H 0.09 1.26 8.03 -51.93 -0.42 0.85 16 29 C -0.16 -1.90 9.10 37.16 0.34 -1.57 16 30 H 0.06 0.38 8.14 -51.93 -0.42 -0.05 16 31 H 0.06 0.26 8.14 -51.93 -0.42 -0.16 16 32 H 0.05 0.28 6.47 -51.93 -0.34 -0.06 16 33 H 0.08 0.51 7.93 -51.93 -0.41 0.10 16 34 H 0.03 0.23 8.12 -51.93 -0.42 -0.20 16 35 H 0.07 0.37 6.47 -51.92 -0.34 0.04 16 36 H 0.07 0.37 6.30 -51.93 -0.33 0.05 16 37 H 0.03 0.16 8.14 -51.93 -0.42 -0.26 16 38 H 0.08 0.81 7.89 -51.93 -0.41 0.40 16 39 H 0.04 0.42 8.08 -51.93 -0.42 0.00 16 40 H 0.09 0.59 6.55 -51.93 -0.34 0.25 16 41 H 0.10 0.72 7.85 -51.93 -0.41 0.31 16 42 H 0.07 0.44 8.14 -51.93 -0.42 0.02 16 43 H 0.08 0.46 7.48 -51.93 -0.39 0.07 16 44 H 0.09 0.67 8.14 -51.93 -0.42 0.24 16 45 H 0.08 0.77 7.81 -51.93 -0.41 0.36 16 46 H 0.35 5.63 8.00 -39.29 -0.31 5.31 16 47 H 0.05 0.65 8.13 -51.93 -0.42 0.23 16 48 H 0.09 1.32 8.08 -51.93 -0.42 0.90 16 49 H 0.39 9.26 7.90 -40.82 -0.32 8.94 16 50 H 0.27 6.88 8.31 -40.82 -0.34 6.54 16 51 H 0.05 0.55 8.14 -51.93 -0.42 0.13 16 52 H 0.07 0.91 6.90 -51.93 -0.36 0.55 16 53 H 0.07 0.93 8.14 -51.93 -0.42 0.51 16 LS Contribution 420.41 15.07 6.34 6.34 Total: -1.00 -32.68 420.41 -9.39 -42.07 By element: Atomic # 1 Polarization: 20.36 SS G_CDS: -9.32 Total: 11.03 kcal Atomic # 6 Polarization: 8.62 SS G_CDS: 0.07 Total: 8.69 kcal Atomic # 7 Polarization: -58.77 SS G_CDS: -1.62 Total: -60.39 kcal Atomic # 8 Polarization: -23.80 SS G_CDS: 0.17 Total: -23.63 kcal Atomic # 14 Polarization: 20.91 SS G_CDS: -5.02 Total: 15.89 kcal Total LS contribution 6.34 Total: 6.34 kcal Total: -32.68 -9.39 -42.07 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. ZINC001486093324.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 -23.662 kcal (2) G-P(sol) polarization free energy of solvation -32.678 kcal (3) G-ENP(sol) elect.-nuc.-pol. free energy of system -56.339 kcal (4) G-CDS(sol) cavity-dispersion-solvent structure free energy -9.389 kcal (5) G-P-CDS(sol) = G-P(sol) + G-CDS(sol) = (2) + (4) -42.066 kcal (6) G-S(sol) free energy of system = (1) + (5) -65.728 kcal 1SCF WAS SPECIFIED, SO GEOMETRY OPTIMIZATION WAS NOT USED Total computer time = 0.88 seconds