Wall clock time and date at job start Tue Mar 9 2021 13:22:06 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 N 1.46502 * 109.46882 * 2 1 4 4 C 1.34772 * 119.99933 * 180.02562 * 3 2 1 5 5 O 1.21859 * 120.00562 * 0.02562 * 4 3 2 6 6 C 1.40194 * 119.99696 * 179.97438 * 4 3 2 7 7 S 1.55360 * 119.99567 * 180.02562 * 6 4 3 8 8 O 1.52204 * 102.99970 * 286.87161 * 7 6 4 9 9 C 1.76206 * 102.99683 * 180.02562 * 7 6 4 10 10 C 1.38259 * 119.96611 * 29.97578 * 9 7 6 11 11 C 1.38304 * 120.14671 * 180.02562 * 10 9 7 12 12 C 1.38083 * 120.07263 * 359.72268 * 11 10 9 13 13 C 1.38923 * 119.92494 * 0.57742 * 12 11 10 14 14 N 1.39832 * 120.07165 * 179.70062 * 13 12 11 15 15 C 1.34774 * 119.99948 * 211.91724 * 14 13 12 16 16 O 1.21705 * 119.99909 * 5.19518 * 15 14 13 17 17 C 1.46648 * 120.00184 * 185.19901 * 15 14 13 18 18 C 1.37891 * 126.23421 * 179.97438 * 17 15 14 19 19 C 1.50699 * 126.46006 * 359.97438 * 18 17 15 20 20 C 1.39829 * 107.08330 * 179.97438 * 18 17 15 21 21 C 1.50702 * 126.06321 * 179.97438 * 20 18 17 22 22 O 1.42892 * 109.47523 * 308.37743 * 21 20 18 23 23 Si 1.86306 * 109.47134 * 68.38048 * 21 20 18 24 24 H 1.48501 * 109.46984 * 65.29582 * 23 21 20 25 25 H 1.48494 * 109.47014 * 185.29876 * 23 21 20 26 26 H 1.48501 * 109.46678 * 305.29683 * 23 21 20 27 27 C 1.35456 * 107.86906 * 0.02562 * 20 18 17 28 28 C 1.50696 * 125.56129 * 179.97438 * 27 20 18 29 29 N 1.35543 * 108.87047 * 0.24138 * 27 20 18 30 30 C 1.38131 * 119.96168 * 210.00181 * 9 7 6 31 31 H 1.09007 * 109.46985 * 299.99394 * 1 2 3 32 32 H 1.08993 * 109.47569 * 59.99607 * 1 2 3 33 33 H 1.09005 * 109.46999 * 180.02562 * 1 2 3 34 34 H 1.08993 * 109.47092 * 240.00660 * 2 1 3 35 35 H 1.08997 * 109.47326 * 119.99698 * 2 1 3 36 36 H 0.97004 * 120.00382 * 0.02562 * 3 2 1 37 37 H 1.08000 * 120.00383 * 0.02562 * 6 4 3 38 38 H 1.08001 * 119.92254 * 359.97438 * 10 9 7 39 39 H 1.08003 * 119.96301 * 180.02562 * 11 10 9 40 40 H 1.07996 * 120.03943 * 180.27691 * 12 11 10 41 41 H 0.97001 * 120.00113 * 31.91294 * 14 13 12 42 42 H 1.08997 * 109.47009 * 268.14716 * 19 18 17 43 43 H 1.09006 * 109.47018 * 28.13782 * 19 18 17 44 44 H 1.08997 * 109.47519 * 148.14444 * 19 18 17 45 45 H 1.09004 * 109.46810 * 188.37335 * 21 20 18 46 46 H 0.96698 * 114.00040 * 59.99929 * 22 21 20 47 47 H 1.09004 * 109.47205 * 84.55367 * 28 27 20 48 48 H 1.08991 * 109.47598 * 324.55199 * 28 27 20 49 49 H 1.09001 * 109.46827 * 204.54542 * 28 27 20 50 50 H 0.96999 * 125.67334 * 179.73069 * 29 27 20 51 51 H 1.08001 * 120.03150 * 0.02562 * 30 9 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 6 1.5300 0.0000 0.0000 3 7 2.0183 1.3813 0.0000 4 6 3.3433 1.6276 -0.0005 5 8 4.1353 0.7014 -0.0015 6 6 3.8105 2.9494 -0.0011 7 16 5.3380 3.2333 -0.0023 8 8 5.7543 2.8729 1.4167 9 6 5.4138 4.9937 -0.0018 10 6 4.4075 5.7281 0.5977 11 6 4.4630 7.1101 0.5998 12 6 5.5273 7.7607 0.0077 13 6 6.5370 7.0259 -0.6010 14 7 7.6130 7.6818 -1.2070 15 6 8.2122 7.1391 -2.2854 16 8 7.7588 6.1281 -2.7891 17 6 9.4105 7.7666 -2.8519 18 6 10.1229 7.3232 -3.9460 19 6 9.8171 6.1062 -4.7806 20 6 11.1827 8.2155 -4.1353 21 6 12.2302 8.1326 -5.2156 22 8 12.8079 6.8256 -5.2197 23 14 11.4435 8.4691 -6.8705 24 1 10.9564 9.8713 -6.9110 25 1 12.4459 8.2619 -7.9463 26 1 10.3011 7.5421 -7.0731 27 6 11.1015 9.1670 -3.1746 28 6 12.0492 10.3267 -3.0073 29 7 10.0321 8.9024 -2.3848 30 6 6.4784 5.6385 -0.6009 31 1 -0.3633 0.5138 0.8901 32 1 -0.3634 0.5138 -0.8899 33 1 -0.3633 -1.0277 0.0005 34 1 1.8933 -0.5137 -0.8900 35 1 1.8934 -0.5138 0.8900 36 1 1.3880 2.1186 0.0004 37 1 3.1087 3.7702 -0.0006 38 1 3.5758 5.2215 1.0647 39 1 3.6755 7.6808 1.0693 40 1 5.5724 8.8397 0.0135 41 1 7.9282 8.5274 -0.8513 42 1 10.3606 5.2483 -4.3848 43 1 8.7464 5.9039 -4.7493 44 1 10.1230 6.2861 -5.8112 45 1 13.0074 8.8728 -5.0256 46 1 12.1725 6.1153 -5.3832 47 1 11.7404 11.1450 -3.6579 48 1 13.0584 10.0124 -3.2729 49 1 12.0339 10.6616 -1.9701 50 1 9.7506 9.4275 -1.6194 51 1 7.2638 5.0641 -1.0697 RHF calculation, no. of doubly occupied orbitals= 73 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). S: (AM1): M.J.S. DEWAR, Y.-C. YUAN, INORGANIC CHEM., 29, 589 (1990) 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_9704022_12297754.mol2 51 1SCF WAS SPECIFIED, SO GEOMETRY OPTIMIZATION WAS NOT USED SCF FIELD WAS ACHIEVED Wall clock time and date at job start Tue Mar 9 2021 13:22:06 Heat of formation + Delta-G solvation = -119.750973 kcal Electronic energy + Delta-G solvation = -35380.742105 eV Core-core repulsion = 30564.775504 eV Total energy + Delta-G solvation = -4815.966601 eV No. of doubly occupied orbitals = 73 Molecular weight (most abundant/longest-lived isotopes) = 406.144 amu Computer time = 1.09 seconds Orbital eigenvalues (eV) -41.96750 -40.02958 -37.91792 -37.65417 -37.08178 -34.79979 -33.84710 -32.97474 -32.35258 -31.43111 -30.45041 -29.90232 -29.43014 -28.64284 -27.12246 -26.49918 -25.59561 -22.97875 -22.38293 -21.94407 -21.05857 -20.03591 -19.55049 -19.23504 -18.55568 -17.59486 -17.11279 -16.01646 -15.97155 -15.71620 -15.67054 -15.12124 -15.07755 -14.99613 -14.71497 -14.43976 -14.34356 -14.13302 -13.85075 -13.73434 -13.52732 -13.39472 -13.21802 -13.11166 -12.95052 -12.76304 -12.64039 -12.53785 -12.42675 -12.19625 -11.99188 -11.76036 -11.66969 -11.43546 -11.39712 -11.19740 -11.05213 -10.98555 -10.79496 -10.69678 -10.33304 -9.99591 -9.92594 -9.77314 -8.77093 -8.44918 -8.30287 -8.06995 -7.98585 -7.74261 -7.68296 -6.83439 -5.40603 1.09524 1.70354 1.78638 1.89643 2.11754 2.16158 2.38778 2.66052 2.73547 2.93927 3.21902 3.82597 4.22545 4.26461 4.38201 4.57281 4.68877 4.70753 4.86082 4.99076 5.06689 5.07030 5.15869 5.24117 5.28992 5.32031 5.41709 5.45218 5.59012 5.60572 5.68277 5.71786 5.77704 5.83970 5.89050 5.95199 6.04628 6.12864 6.18374 6.20884 6.40843 6.52682 6.57648 6.60332 6.63604 6.68907 6.70400 6.82862 7.02158 7.14245 7.30625 7.36611 7.47049 7.68392 7.94483 8.05695 8.08966 9.22916 10.35305 Molecular weight = 406.14amu Principal moments of inertia in cm(-1) A = 0.012369 B = 0.001689 C = 0.001569 Principal moments of inertia in units of 10**(-40)*gram-cm**2 A = 2263.150212 B =16569.546394 C =17846.444205 Net atomic charges, atomic populations, and dipole contributions calculated with CM2 Atom NO. Type Charge No. of electrons 1 C -0.159 4.159 2 C 0.124 3.876 3 N -0.762 5.762 4 C 0.609 3.391 5 O -0.642 6.642 6 C -1.081 5.081 7 S 1.464 4.536 8 O -0.935 6.935 9 C -0.389 4.389 10 C -0.074 4.074 11 C -0.124 4.124 12 C -0.140 4.140 13 C 0.119 3.881 14 N -0.672 5.672 15 C 0.584 3.416 16 O -0.495 6.495 17 C -0.052 4.052 18 C -0.057 4.057 19 C -0.114 4.114 20 C -0.211 4.211 21 C -0.037 4.037 22 O -0.542 6.542 23 Si 0.935 3.065 24 H -0.254 1.254 25 H -0.254 1.254 26 H -0.239 1.239 27 C 0.074 3.926 28 C -0.095 4.095 29 N -0.536 5.536 30 C -0.058 4.058 31 H 0.053 0.947 32 H 0.055 0.945 33 H 0.061 0.939 34 H 0.050 0.950 35 H 0.048 0.952 36 H 0.382 0.618 37 H 0.101 0.899 38 H 0.123 0.877 39 H 0.117 0.883 40 H 0.115 0.885 41 H 0.398 0.602 42 H 0.078 0.922 43 H 0.092 0.908 44 H 0.065 0.935 45 H 0.110 0.890 46 H 0.375 0.625 47 H 0.079 0.921 48 H 0.088 0.912 49 H 0.069 0.931 50 H 0.416 0.584 51 H 0.140 0.860 For the charges calculated by CM2: Dipole moment (debyes) X Y Z Total from point charges 12.302 19.783 -12.727 26.545 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.217 4.217 2 C -0.002 4.002 3 N -0.415 5.415 4 C 0.417 3.583 5 O -0.536 6.536 6 C -1.213 5.213 7 S 1.665 4.335 8 O -0.930 6.930 9 C -0.480 4.480 10 C -0.093 4.093 11 C -0.143 4.143 12 C -0.160 4.160 13 C 0.025 3.975 14 N -0.316 5.316 15 C 0.373 3.627 16 O -0.369 6.369 17 C -0.164 4.164 18 C -0.065 4.065 19 C -0.171 4.171 20 C -0.217 4.217 21 C -0.148 4.148 22 O -0.348 6.348 23 Si 0.757 3.243 24 H -0.178 1.178 25 H -0.178 1.178 26 H -0.163 1.163 27 C -0.044 4.044 28 C -0.153 4.153 29 N -0.134 5.134 30 C -0.079 4.079 31 H 0.072 0.928 32 H 0.074 0.926 33 H 0.080 0.920 34 H 0.068 0.932 35 H 0.066 0.934 36 H 0.213 0.787 37 H 0.120 0.880 38 H 0.140 0.860 39 H 0.135 0.865 40 H 0.133 0.867 41 H 0.232 0.768 42 H 0.097 0.903 43 H 0.111 0.889 44 H 0.083 0.917 45 H 0.128 0.872 46 H 0.223 0.777 47 H 0.098 0.902 48 H 0.106 0.894 49 H 0.088 0.912 50 H 0.253 0.747 51 H 0.157 0.843 Dipole moment (debyes) X Y Z Total from point charges 13.321 18.370 -13.642 26.476 hybrid contribution -2.884 2.205 2.758 4.559 sum 10.437 20.574 -10.884 25.509 Atomic orbital electron populations 1.21976 0.96271 1.01751 1.01687 1.20898 0.92273 0.86241 1.00748 1.45818 1.08129 1.05476 1.82094 1.15178 0.80226 0.88239 0.74620 1.90603 1.54899 1.47576 1.60488 1.26603 1.20070 1.00627 1.74032 1.66923 0.85080 0.75663 1.05827 1.94470 1.83829 1.77884 1.36848 1.28865 0.96432 1.23946 0.98801 1.21279 0.98420 0.93284 0.96305 1.20994 0.97726 0.95049 1.00575 1.20670 0.94153 0.98502 1.02697 1.17702 0.88705 0.92248 0.98836 1.43261 1.33143 1.25096 1.30091 1.17064 0.81939 0.84275 0.79416 1.90883 1.63202 1.23490 1.59321 1.20019 1.00014 0.93861 1.02486 1.19700 0.93868 0.97544 0.95436 1.20735 1.05114 0.93907 0.97343 1.19470 1.01074 0.99091 1.02097 1.24204 0.97736 0.88461 1.04418 1.86469 1.38203 1.14941 1.95198 0.86014 0.80273 0.81736 0.76296 1.17804 1.17767 1.16256 1.20198 0.91382 0.99505 0.93332 1.20493 0.97717 0.94119 1.02921 1.41839 1.23221 1.22470 1.25843 1.20755 0.98166 0.91668 0.97276 0.92757 0.92629 0.91964 0.93153 0.93373 0.78671 0.88042 0.85962 0.86528 0.86717 0.76797 0.90306 0.88894 0.91670 0.87214 0.77712 0.90198 0.89363 0.91243 0.74687 0.84264 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.16 -2.40 10.03 37.16 0.37 -2.02 16 2 C 0.12 2.56 6.39 -4.04 -0.03 2.53 16 3 N -0.76 -18.71 5.54 -63.64 -0.35 -19.06 16 4 C 0.61 18.27 7.50 -15.53 -0.12 18.15 16 5 O -0.64 -21.13 16.28 5.05 0.08 -21.05 16 6 C -1.08 -31.39 8.32 31.47 0.26 -31.13 16 7 S 1.46 39.82 12.78 -107.50 -1.37 38.44 16 8 O -0.93 -31.97 18.25 -35.23 -0.64 -32.62 16 9 C -0.39 -9.08 5.67 -39.61 -0.22 -9.31 16 10 C -0.07 -1.63 8.93 -39.55 -0.35 -1.98 16 11 C -0.12 -2.20 10.04 -39.61 -0.40 -2.59 16 12 C -0.14 -2.18 9.96 -39.38 -0.39 -2.57 16 13 C 0.12 1.99 6.30 -83.72 -0.53 1.46 16 14 N -0.67 -8.52 5.34 -11.05 -0.06 -8.58 16 15 C 0.58 8.13 7.53 -12.87 -0.10 8.03 16 16 O -0.50 -8.50 12.36 5.18 0.06 -8.44 16 17 C -0.05 -0.59 7.05 -81.99 -0.58 -1.16 16 18 C -0.06 -0.63 5.87 -104.15 -0.61 -1.24 16 19 C -0.11 -1.19 8.17 36.01 0.29 -0.89 16 20 C -0.21 -1.91 4.44 -105.02 -0.47 -2.37 16 21 C -0.04 -0.28 2.85 36.01 0.10 -0.18 16 22 O -0.54 -4.81 12.87 -35.23 -0.45 -5.27 16 23 Si 0.93 6.83 26.86 -169.99 -4.57 2.27 16 24 H -0.25 -1.98 7.11 56.52 0.40 -1.58 16 25 H -0.25 -1.98 7.11 56.52 0.40 -1.58 16 26 H -0.24 -1.95 5.58 56.52 0.32 -1.64 16 27 C 0.07 0.55 7.13 -81.95 -0.58 -0.03 16 28 C -0.10 -0.38 10.00 36.00 0.36 -0.02 16 29 N -0.54 -4.13 5.99 -5.70 -0.03 -4.17 16 30 C -0.06 -1.16 8.36 -39.39 -0.33 -1.49 16 31 H 0.05 0.79 8.14 -51.93 -0.42 0.37 16 32 H 0.05 0.78 8.14 -51.93 -0.42 0.36 16 33 H 0.06 0.79 8.14 -51.93 -0.42 0.36 16 34 H 0.05 1.05 8.14 -51.93 -0.42 0.63 16 35 H 0.05 1.06 8.14 -51.93 -0.42 0.64 16 36 H 0.38 8.40 8.42 -40.82 -0.34 8.06 16 37 H 0.10 2.83 6.09 -45.47 -0.28 2.55 16 38 H 0.12 2.89 6.48 -52.49 -0.34 2.54 16 39 H 0.12 1.71 8.06 -52.48 -0.42 1.29 16 40 H 0.11 1.34 8.06 -52.49 -0.42 0.92 16 41 H 0.40 3.38 7.96 -40.82 -0.32 3.05 16 42 H 0.08 0.75 8.01 -51.93 -0.42 0.33 16 43 H 0.09 1.15 6.38 -51.93 -0.33 0.82 16 44 H 0.06 0.55 4.40 -51.93 -0.23 0.32 16 45 H 0.11 0.62 6.66 -51.93 -0.35 0.28 16 46 H 0.37 2.61 6.29 45.56 0.29 2.90 16 47 H 0.08 0.23 8.14 -51.93 -0.42 -0.20 16 48 H 0.09 0.30 7.23 -51.93 -0.38 -0.08 16 49 H 0.07 0.19 8.14 -51.93 -0.42 -0.23 16 50 H 0.42 1.91 8.20 -40.82 -0.33 1.57 16 51 H 0.14 2.73 6.32 -52.49 -0.33 2.40 16 LS Contribution 426.17 15.07 6.42 6.42 Total: -1.00 -44.51 426.17 -10.28 -54.78 By element: Atomic # 1 Polarization: 30.14 SS G_CDS: -6.05 Total: 24.09 kcal Atomic # 6 Polarization: -23.51 SS G_CDS: -3.31 Total: -26.82 kcal Atomic # 7 Polarization: -31.36 SS G_CDS: -0.45 Total: -31.81 kcal Atomic # 8 Polarization: -66.42 SS G_CDS: -0.95 Total: -67.37 kcal Atomic # 14 Polarization: 6.83 SS G_CDS: -4.57 Total: 2.27 kcal Atomic # 16 Polarization: 39.82 SS G_CDS: -1.37 Total: 38.44 kcal Total LS contribution 6.42 Total: 6.42 kcal Total: -44.51 -10.28 -54.78 kcal The number of atoms in the molecule is 51 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_9704022_12297754.mol2 51 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 -64.970 kcal (2) G-P(sol) polarization free energy of solvation -44.505 kcal (3) G-ENP(sol) elect.-nuc.-pol. free energy of system -109.475 kcal (4) G-CDS(sol) cavity-dispersion-solvent structure free energy -10.276 kcal (5) G-P-CDS(sol) = G-P(sol) + G-CDS(sol) = (2) + (4) -54.781 kcal (6) G-S(sol) free energy of system = (1) + (5) -119.751 kcal 1SCF WAS SPECIFIED, SO GEOMETRY OPTIMIZATION WAS NOT USED Total computer time = 1.09 seconds