Wall clock time and date at job start Tue Mar 31 2020 06:57: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 N 1.46494 * 1 3 3 C 1.36880 * 126.40012 * 2 1 4 4 C 1.34673 * 106.81894 * 180.25587 * 3 2 1 5 5 N 1.34126 * 108.01962 * 359.73950 * 4 3 2 6 6 C 1.30503 * 109.25575 * 0.02562 * 5 4 3 7 7 C 1.50697 * 125.64552 * 180.02562 * 6 5 4 8 8 N 1.46902 * 109.47188 * 90.02235 * 7 6 5 9 9 C 1.46264 * 110.99970 * 169.99420 * 8 7 6 10 10 C 1.47542 * 114.37618 * 148.55081 * 9 8 7 11 11 C 1.54947 * 109.72903 * 62.22255 * 10 9 8 12 12 N 1.44678 * 109.52612 * 290.52191 * 11 10 9 13 13 C 1.34772 * 118.54713 * 269.05903 * 12 11 10 14 14 O 1.21553 * 120.00020 * 172.22894 * 13 12 11 15 15 C 1.47886 * 120.00373 * 352.22964 * 13 12 11 16 16 C 1.39917 * 120.02113 * 27.41326 * 15 13 12 17 17 C 1.38330 * 120.35384 * 179.97438 * 16 15 13 18 18 C 1.37753 * 120.41494 * 0.02562 * 17 16 15 19 19 C 1.40745 * 120.04045 * 359.97438 * 18 17 16 20 20 C 1.41429 * 120.18176 * 180.02562 * 19 18 17 21 21 H 1.08005 * 120.00169 * 24.77712 * 20 19 18 22 22 C 1.39937 * 120.00467 * 204.77597 * 20 19 18 23 23 N 1.30824 * 119.99893 * 17.70621 * 22 20 19 24 24 Si 1.86796 * 120.00092 * 197.70811 * 22 20 19 25 25 H 1.48507 * 109.47112 * 90.00113 * 24 22 20 26 26 H 1.48498 * 109.47285 * 209.99844 * 24 22 20 27 27 H 1.48497 * 109.47531 * 329.99996 * 24 22 20 28 28 C 1.38998 * 120.02471 * 207.14210 * 15 13 12 29 29 C 1.47784 * 122.90420 * 89.07813 * 12 11 10 30 30 C 1.44761 * 110.99970 * 46.51149 * 8 7 6 31 31 H 1.09001 * 109.47370 * 270.18596 * 1 2 3 32 32 H 1.08998 * 109.47353 * 30.19430 * 1 2 3 33 33 H 1.09007 * 109.46916 * 150.18757 * 1 2 3 34 34 H 1.08001 * 126.58972 * 359.97155 * 3 2 1 35 35 H 1.08000 * 125.98561 * 179.72345 * 4 3 2 36 36 H 1.09005 * 109.46731 * 210.01907 * 7 6 5 37 37 H 1.09000 * 109.47607 * 330.02341 * 7 6 5 38 38 H 1.08997 * 108.50740 * 269.82439 * 9 8 7 39 39 H 1.09006 * 108.50261 * 27.44621 * 9 8 7 40 40 H 1.08998 * 109.43397 * 182.30195 * 10 9 8 41 41 H 1.09004 * 109.42945 * 302.14703 * 10 9 8 42 42 H 1.08996 * 109.46218 * 50.54210 * 11 10 9 43 43 H 1.08999 * 109.46235 * 170.52227 * 11 10 9 44 44 H 1.07994 * 119.82498 * 0.02562 * 16 15 13 45 45 H 1.08008 * 119.79324 * 180.02562 * 17 16 15 46 46 H 1.07998 * 119.97883 * 179.97438 * 18 17 16 47 47 H 0.96997 * 119.99990 * 179.97438 * 23 22 20 48 48 H 1.08001 * 120.20403 * 359.97438 * 28 15 13 49 49 H 1.08998 * 109.06373 * 78.64499 * 29 12 11 50 50 H 1.09001 * 109.06058 * 197.64831 * 29 12 11 51 51 H 1.09000 * 109.43376 * 337.91732 * 30 8 7 52 52 H 1.08997 * 109.44076 * 97.79152 * 30 8 7 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 7 1.4649 0.0000 0.0000 3 6 2.2772 1.1017 0.0000 4 6 3.5460 0.6504 -0.0058 5 7 3.5095 -0.6904 -0.0034 6 6 2.2662 -1.0870 0.0031 7 6 1.8018 -2.5206 0.0082 8 7 1.6567 -2.9834 1.3948 9 6 0.9866 -4.2825 1.4464 10 6 1.4092 -5.1416 2.5690 11 6 1.0608 -4.4545 3.9134 12 7 1.9111 -3.2990 4.1006 13 6 3.0792 -3.4581 4.7538 14 8 3.8772 -2.5425 4.8036 15 6 3.3900 -4.7464 5.4100 16 6 2.8297 -5.9327 4.9240 17 6 3.1162 -7.1417 5.5322 18 6 3.9556 -7.1918 6.6233 19 6 4.5287 -6.0088 7.1263 20 6 5.3914 -6.0568 8.2460 21 1 5.9293 -6.9650 8.4747 22 6 5.5481 -4.9277 9.0576 23 7 4.6688 -3.9602 9.0092 24 14 7.0126 -4.8020 10.2102 25 1 8.1425 -4.1451 9.5051 26 1 6.6378 -3.9978 11.4010 27 1 7.4221 -6.1629 10.6409 28 6 4.2435 -4.7802 6.5066 29 6 1.5553 -1.9521 3.6075 30 6 0.9442 -2.0098 2.1948 31 1 -0.3634 0.0033 1.0276 32 1 -0.3634 0.8882 -0.5168 33 1 -0.3633 -0.8917 -0.5110 34 1 1.9613 2.1345 -0.0004 35 1 4.4368 1.2611 -0.0116 36 1 0.8410 -2.5930 -0.5015 37 1 2.5343 -3.1418 -0.5073 38 1 -0.0876 -4.1149 1.5243 39 1 1.1861 -4.8119 0.5146 40 1 0.8929 -6.0995 2.5070 41 1 2.4853 -5.3064 2.5143 42 1 0.0174 -4.1393 3.9004 43 1 1.2178 -5.1571 4.7318 44 1 2.1706 -5.9046 4.0690 45 1 2.6789 -8.0520 5.1493 46 1 4.1742 -8.1394 7.0929 47 1 4.7772 -3.1777 9.5720 48 1 4.6821 -3.8675 6.8822 49 1 0.8336 -1.5025 4.2894 50 1 2.4528 -1.3341 3.5812 51 1 1.0214 -1.0288 1.7259 52 1 -0.1052 -2.2960 2.2650 RHF calculation, no. of doubly occupied orbitals= 69 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 ZINC000069895585.mol2 52 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 06:57:36 Heat of formation + Delta-G solvation = 120.777424 kcal Electronic energy + Delta-G solvation = -34378.304038 eV Core-core repulsion = 30153.272420 eV Total energy + Delta-G solvation = -4225.031618 eV No. of doubly occupied orbitals = 69 Molecular weight (most abundant/longest-lived isotopes) = 368.195 amu Computer time = 1.55 seconds Orbital eigenvalues (eV) -42.44055 -40.18184 -38.17148 -36.87966 -35.33237 -34.30053 -33.64251 -32.81708 -31.14633 -29.77467 -29.42735 -28.44154 -27.28779 -26.38704 -25.28786 -24.84366 -23.04966 -21.60798 -21.29257 -20.93133 -20.62273 -19.63280 -18.54823 -17.54002 -17.26219 -16.93283 -16.46379 -15.75691 -15.38430 -15.30397 -15.10574 -14.87019 -14.77206 -14.38002 -13.98260 -13.84159 -13.45941 -13.31661 -13.28034 -13.21357 -13.06140 -12.94301 -12.71177 -12.38347 -12.32694 -12.19670 -11.82409 -11.55992 -11.43226 -11.19464 -11.13126 -11.04788 -10.88470 -10.69340 -10.39623 -10.15138 -10.03133 -9.90347 -9.78384 -9.43715 -9.12440 -8.78284 -8.62036 -8.36395 -8.04732 -7.49912 -7.38253 -6.74463 -4.33989 1.78172 2.38305 2.48299 2.87757 2.98116 3.12021 3.14998 3.21540 3.57182 3.86757 3.92450 3.97328 4.15617 4.27575 4.49447 4.62613 4.65230 4.76636 4.81794 4.88178 4.94048 5.01897 5.13433 5.29044 5.31200 5.36204 5.41049 5.54753 5.56560 5.64434 5.66770 5.72894 5.84105 5.89607 6.02991 6.18030 6.19020 6.26542 6.39200 6.48887 6.60179 6.79334 6.81647 6.92995 7.03127 7.06356 7.11655 7.20174 7.28957 7.47573 7.77134 7.87660 8.00288 8.00845 8.14252 8.26616 8.58615 8.87829 8.91917 9.09556 10.20850 Molecular weight = 368.20amu Principal moments of inertia in cm(-1) A = 0.013739 B = 0.002571 C = 0.002307 Principal moments of inertia in units of 10**(-40)*gram-cm**2 A = 2037.552738 B =10888.791697 C =12133.314639 Net atomic charges, atomic populations, and dipole contributions calculated with CM2 Atom NO. Type Charge No. of electrons 1 C 0.066 3.934 2 N -0.464 5.464 3 C -0.059 4.059 4 C -0.024 4.024 5 N -0.464 5.464 6 C 0.209 3.791 7 C 0.146 3.854 8 N -0.518 5.518 9 C 0.068 3.932 10 C -0.146 4.146 11 C 0.084 3.916 12 N -0.596 5.596 13 C 0.566 3.434 14 O -0.547 6.547 15 C -0.121 4.121 16 C -0.256 4.256 17 C -0.101 4.101 18 C -0.213 4.213 19 C 0.100 3.900 20 C -0.461 4.461 21 H 0.080 0.920 22 C 0.074 3.926 23 N -0.789 5.789 24 Si 0.910 3.090 25 H -0.266 1.266 26 H -0.278 1.278 27 H -0.262 1.262 28 C -0.096 4.096 29 C 0.125 3.875 30 C 0.033 3.967 31 H 0.080 0.920 32 H 0.088 0.912 33 H 0.085 0.915 34 H 0.174 0.826 35 H 0.170 0.830 36 H 0.045 0.955 37 H 0.106 0.894 38 H 0.034 0.966 39 H 0.066 0.934 40 H 0.082 0.918 41 H 0.088 0.912 42 H 0.056 0.944 43 H 0.128 0.872 44 H 0.138 0.862 45 H 0.091 0.909 46 H 0.093 0.907 47 H 0.283 0.717 48 H 0.128 0.872 49 H 0.054 0.946 50 H 0.089 0.911 51 H 0.080 0.920 52 H 0.041 0.959 For the charges calculated by CM2: Dipole moment (debyes) X Y Z Total from point charges -14.222 6.500 -17.358 23.363 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.077 4.077 2 N -0.148 5.148 3 C -0.191 4.191 4 C -0.172 4.172 5 N -0.182 5.182 6 C -0.053 4.053 7 C 0.017 3.983 8 N -0.243 5.243 9 C -0.059 4.059 10 C -0.185 4.185 11 C -0.040 4.040 12 N -0.329 5.329 13 C 0.358 3.642 14 O -0.425 6.425 15 C -0.124 4.124 16 C -0.274 4.274 17 C -0.120 4.120 18 C -0.233 4.233 19 C 0.098 3.902 20 C -0.490 4.490 21 H 0.098 0.902 22 C -0.144 4.144 23 N -0.413 5.413 24 Si 0.734 3.266 25 H -0.191 1.191 26 H -0.203 1.203 27 H -0.186 1.186 28 C -0.116 4.116 29 C 0.002 3.998 30 C -0.097 4.097 31 H 0.098 0.902 32 H 0.106 0.894 33 H 0.103 0.897 34 H 0.191 0.809 35 H 0.188 0.812 36 H 0.063 0.937 37 H 0.124 0.876 38 H 0.052 0.948 39 H 0.085 0.915 40 H 0.101 0.899 41 H 0.106 0.894 42 H 0.074 0.926 43 H 0.145 0.855 44 H 0.155 0.845 45 H 0.110 0.890 46 H 0.111 0.889 47 H 0.094 0.906 48 H 0.146 0.854 49 H 0.072 0.928 50 H 0.107 0.893 51 H 0.098 0.902 52 H 0.059 0.941 Dipole moment (debyes) X Y Z Total from point charges -12.890 6.669 -17.041 22.383 hybrid contribution -1.008 0.950 0.249 1.408 sum -13.898 7.619 -16.792 23.091 Atomic orbital electron populations 1.22129 0.77114 1.04641 1.03768 1.45448 1.07309 1.06847 1.55241 1.22693 0.91422 0.92369 1.12617 1.22377 0.96477 0.91219 1.07134 1.70302 1.18216 1.11846 1.17869 1.22950 0.85893 0.90302 1.06179 1.20504 0.99838 0.87333 0.90583 1.60889 1.52460 1.02627 1.08355 1.20971 0.95408 0.89780 0.99767 1.21569 1.04103 0.97932 0.94890 1.22574 0.93174 0.89350 0.98881 1.47284 1.20802 1.07988 1.56806 1.18932 0.80134 0.88353 0.76808 1.90670 1.47466 1.43591 1.60798 1.19462 1.02579 0.91290 0.99069 1.21769 1.08150 0.92108 1.05402 1.20495 0.98848 0.95841 0.96775 1.21111 1.03491 0.97680 1.00998 1.17575 0.90326 0.90667 0.91656 1.19578 1.22876 0.98306 1.08281 0.90176 1.25933 0.96260 0.94805 0.97407 1.69789 1.30936 1.07088 1.33525 0.86289 0.81764 0.79344 0.79208 1.19129 1.20317 1.18619 1.21540 0.97284 0.99435 0.93357 1.21673 1.00833 0.85185 0.92109 1.22738 0.95145 0.96636 0.95153 0.90199 0.89376 0.89669 0.80854 0.81240 0.93681 0.87551 0.94761 0.91535 0.89950 0.89383 0.92557 0.85466 0.84504 0.89047 0.88869 0.90643 0.85431 0.92780 0.89267 0.90159 0.94083 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.07 0.34 9.85 59.84 0.59 0.93 16 2 N -0.46 -4.39 2.40 -121.28 -0.29 -4.68 16 3 C -0.06 -0.57 11.58 -16.56 -0.19 -0.76 16 4 C -0.02 -0.29 11.87 -17.81 -0.21 -0.51 16 5 N -0.46 -7.14 11.25 -10.82 -0.12 -7.27 16 6 C 0.21 2.59 6.01 -155.41 -0.93 1.66 16 7 C 0.15 1.56 5.31 -4.98 -0.03 1.53 16 8 N -0.52 -6.46 4.30 -241.04 -1.04 -7.50 16 9 C 0.07 0.71 4.92 -6.80 -0.03 0.67 16 10 C -0.15 -1.80 4.36 -28.34 -0.12 -1.92 16 11 C 0.08 1.22 4.92 -3.99 -0.02 1.20 16 12 N -0.60 -10.80 2.54 -173.31 -0.44 -11.23 16 13 C 0.57 13.25 7.44 -12.31 -0.09 13.16 16 14 O -0.55 -14.79 15.91 5.32 0.08 -14.70 16 15 C -0.12 -2.94 5.59 -105.04 -0.59 -3.53 16 16 C -0.26 -5.33 7.13 -38.92 -0.28 -5.61 16 17 C -0.10 -2.15 10.01 -39.72 -0.40 -2.55 16 18 C -0.21 -5.12 9.78 -38.81 -0.38 -5.50 16 19 C 0.10 2.71 5.57 -106.87 -0.59 2.11 16 20 C -0.46 -12.42 9.05 -37.73 -0.34 -12.77 16 21 H 0.08 2.06 7.91 -52.48 -0.41 1.64 16 22 C 0.07 1.93 5.14 -17.34 -0.09 1.84 16 23 N -0.79 -21.26 10.99 55.55 0.61 -20.65 16 24 Si 0.91 19.37 29.59 -169.99 -5.03 14.33 16 25 H -0.27 -5.69 7.11 56.52 0.40 -5.29 16 26 H -0.28 -5.76 7.11 56.52 0.40 -5.36 16 27 H -0.26 -5.43 7.11 56.52 0.40 -5.03 16 28 C -0.10 -2.64 8.23 -38.44 -0.32 -2.96 16 29 C 0.13 2.00 6.50 -2.80 -0.02 1.98 16 30 C 0.03 0.39 4.83 -4.43 -0.02 0.37 16 31 H 0.08 0.38 7.03 -51.93 -0.36 0.02 16 32 H 0.09 0.27 8.14 -51.93 -0.42 -0.15 16 33 H 0.08 0.34 7.04 -51.93 -0.37 -0.02 16 34 H 0.17 1.07 8.06 -52.49 -0.42 0.64 16 35 H 0.17 1.76 8.06 -52.49 -0.42 1.34 16 36 H 0.05 0.35 6.99 -51.93 -0.36 -0.02 16 37 H 0.11 1.18 8.14 -51.93 -0.42 0.76 16 38 H 0.03 0.30 7.18 -51.93 -0.37 -0.07 16 39 H 0.07 0.59 7.68 -51.93 -0.40 0.19 16 40 H 0.08 0.80 8.13 -51.93 -0.42 0.38 16 41 H 0.09 1.36 7.14 -51.93 -0.37 0.99 16 42 H 0.06 0.68 8.09 -51.93 -0.42 0.26 16 43 H 0.13 2.10 4.23 -52.37 -0.22 1.88 16 44 H 0.14 2.34 2.11 -52.49 -0.11 2.23 16 45 H 0.09 1.62 8.06 -52.48 -0.42 1.20 16 46 H 0.09 2.09 8.06 -52.49 -0.42 1.67 16 47 H 0.28 7.01 8.30 -40.82 -0.34 6.67 16 48 H 0.13 3.82 6.01 -52.49 -0.32 3.51 16 49 H 0.05 0.78 8.14 -51.93 -0.42 0.35 16 50 H 0.09 1.73 6.91 -51.93 -0.36 1.37 16 51 H 0.08 0.89 2.87 -51.93 -0.15 0.75 16 52 H 0.04 0.40 7.14 -51.93 -0.37 0.03 16 LS Contribution 397.82 15.07 6.00 6.00 Total: -1.00 -35.01 397.82 -11.41 -46.42 By element: Atomic # 1 Polarization: 17.06 SS G_CDS: -7.11 Total: 9.95 kcal Atomic # 6 Polarization: -6.59 SS G_CDS: -4.07 Total: -10.65 kcal Atomic # 7 Polarization: -50.06 SS G_CDS: -1.28 Total: -51.33 kcal Atomic # 8 Polarization: -14.79 SS G_CDS: 0.08 Total: -14.70 kcal Atomic # 14 Polarization: 19.37 SS G_CDS: -5.03 Total: 14.33 kcal Total LS contribution 6.00 Total: 6.00 kcal Total: -35.01 -11.41 -46.42 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. ZINC000069895585.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 167.193 kcal (2) G-P(sol) polarization free energy of solvation -35.006 kcal (3) G-ENP(sol) elect.-nuc.-pol. free energy of system 132.188 kcal (4) G-CDS(sol) cavity-dispersion-solvent structure free energy -11.410 kcal (5) G-P-CDS(sol) = G-P(sol) + G-CDS(sol) = (2) + (4) -46.416 kcal (6) G-S(sol) free energy of system = (1) + (5) 120.777 kcal 1SCF WAS SPECIFIED, SO GEOMETRY OPTIMIZATION WAS NOT USED Total computer time = 1.56 seconds