Wall clock time and date at job start Sat Apr 11 2020 06:10:46 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.53003 * 1 3 3 C 1.52993 * 109.47152 * 2 1 4 4 H 1.08998 * 109.47525 * 56.41164 * 3 2 1 5 5 C 1.50704 * 109.47105 * 176.41469 * 3 2 1 6 6 H 1.07997 * 119.99787 * 120.00085 * 5 3 2 7 7 C 1.37876 * 119.99801 * 299.99825 * 5 3 2 8 8 N 1.31516 * 119.99854 * 180.02562 * 7 5 3 9 9 Si 1.86799 * 120.00231 * 359.97438 * 7 5 3 10 10 H 1.48498 * 109.47002 * 89.99934 * 9 7 5 11 11 H 1.48495 * 109.47231 * 210.00137 * 9 7 5 12 12 H 1.48500 * 109.47129 * 330.00312 * 9 7 5 13 13 N 1.46902 * 109.47073 * 296.41357 * 3 2 1 14 14 C 1.46847 * 111.00309 * 174.38804 * 13 3 2 15 15 C 1.53080 * 109.49530 * 186.12755 * 14 13 3 16 16 C 1.52887 * 109.26699 * 301.28270 * 15 14 13 17 17 C 1.54171 * 109.49804 * 177.27606 * 16 15 14 18 18 C 1.55356 * 102.18838 * 98.54877 * 17 16 15 19 19 H 1.09006 * 111.11755 * 276.02534 * 18 17 16 20 20 O 1.42902 * 111.15359 * 151.89596 * 18 17 16 21 21 C 1.42901 * 114.00000 * 156.94472 * 20 18 17 22 22 C 1.50700 * 109.46999 * 179.97438 * 21 20 18 23 23 O 1.21276 * 120.00215 * 359.97438 * 22 21 20 24 24 N 1.34777 * 119.99840 * 179.97438 * 22 21 20 25 25 C 1.46500 * 120.00030 * 179.97438 * 24 22 21 26 26 C 1.54957 * 100.97094 * 33.95175 * 18 17 16 27 27 O 1.43014 * 110.06648 * 295.64912 * 16 15 14 28 28 C 1.52887 * 109.84215 * 56.99845 * 16 15 14 29 29 C 1.46843 * 110.99815 * 298.69648 * 13 3 2 30 30 H 1.08995 * 109.47116 * 305.07703 * 1 2 3 31 31 H 1.09004 * 109.46658 * 65.07824 * 1 2 3 32 32 H 1.08996 * 109.46999 * 185.07379 * 1 2 3 33 33 H 1.09000 * 109.47499 * 239.99277 * 2 1 3 34 34 H 1.09007 * 109.46806 * 119.99649 * 2 1 3 35 35 H 0.96999 * 119.99932 * 180.02562 * 8 7 5 36 36 H 1.08999 * 109.45959 * 306.13755 * 14 13 3 37 37 H 1.08999 * 109.46339 * 66.11597 * 14 13 3 38 38 H 1.09010 * 109.50575 * 181.34369 * 15 14 13 39 39 H 1.09000 * 109.54336 * 61.24640 * 15 14 13 40 40 H 1.08999 * 110.85699 * 216.72203 * 17 16 15 41 41 H 1.08995 * 110.85730 * 340.37318 * 17 16 15 42 42 H 1.09004 * 109.46877 * 299.99460 * 21 20 18 43 43 H 1.08995 * 109.47531 * 59.99382 * 21 20 18 44 44 H 0.97007 * 119.99683 * 0.02562 * 24 22 21 45 45 H 1.09000 * 109.47017 * 180.02562 * 25 24 22 46 46 H 1.08993 * 109.47190 * 299.99808 * 25 24 22 47 47 H 1.08995 * 109.46978 * 60.00072 * 25 24 22 48 48 H 1.09002 * 110.31144 * 204.75759 * 26 18 17 49 49 H 1.09002 * 110.39298 * 82.62666 * 26 18 17 50 50 H 1.09010 * 109.50087 * 183.07032 * 28 16 15 51 51 H 1.08991 * 109.50928 * 62.93905 * 28 16 15 52 52 H 1.09003 * 109.46703 * 293.87504 * 29 13 3 53 53 H 1.09001 * 109.46527 * 53.86885 * 29 13 3 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 6 2.0400 1.4424 0.0000 4 1 1.6252 1.9745 0.8560 5 6 3.5444 1.4433 0.0889 6 1 4.1315 1.8943 -0.6975 7 6 4.1712 0.8684 1.1740 8 7 5.4841 0.8697 1.2518 9 14 3.1557 0.0891 2.5344 10 1 2.9266 -1.3444 2.2218 11 1 3.8837 0.2044 3.8236 12 1 1.8496 0.7877 2.6393 13 7 1.6224 2.1095 -1.2404 14 6 1.9790 3.5338 -1.2179 15 6 1.3961 4.2260 -2.4526 16 6 1.9425 3.5543 -3.7126 17 6 1.3062 4.2024 -4.9584 18 6 2.4070 5.2034 -5.4057 19 1 2.3332 6.1406 -4.8539 20 8 2.3527 5.4337 -6.8150 21 6 2.9380 6.6713 -7.2247 22 6 2.8262 6.8088 -8.7212 23 8 2.2978 5.9330 -9.3728 24 7 3.3133 7.9050 -9.3355 25 6 3.2052 8.0385 -10.7904 26 6 3.6891 4.4254 -5.0159 27 8 3.3577 3.7433 -3.7951 28 6 1.6129 2.0615 -3.6895 29 6 2.1893 1.4368 -2.4162 30 1 -0.3633 0.5905 0.8410 31 1 -0.3633 0.4331 -0.9320 32 1 -0.3633 -1.0236 0.0909 33 1 1.8934 -0.5139 -0.8899 34 1 1.8933 -0.5138 0.8901 35 1 5.9251 0.4649 2.0151 36 1 1.5740 3.9954 -0.3174 37 1 3.0641 3.6366 -1.2221 38 1 1.6808 5.2782 -2.4484 39 1 0.3093 4.1431 -2.4380 40 1 1.1217 3.4572 -5.7321 41 1 0.3876 4.7273 -4.6964 42 1 2.4146 7.4975 -6.7434 43 1 3.9888 6.6895 -6.9360 44 1 3.7363 8.6054 -8.8142 45 1 3.6515 8.9825 -11.1028 46 1 3.7295 7.2126 -11.2709 47 1 2.1545 8.0195 -11.0798 48 1 4.5156 5.1166 -4.8506 49 1 3.9479 3.7060 -5.7929 50 1 2.0516 1.5784 -4.5626 51 1 0.5314 1.9273 -3.7034 52 1 3.2731 1.5522 -2.4133 53 1 1.9367 0.3769 -2.3853 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 ZINC001255167933.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:10:46 Heat of formation + Delta-G solvation = -112.102788 kcal Electronic energy + Delta-G solvation = -31242.854153 eV Core-core repulsion = 27115.030753 eV Total energy + Delta-G solvation = -4127.823400 eV No. of doubly occupied orbitals = 66 Molecular weight (most abundant/longest-lived isotopes) = 340.220 amu Computer time = 1.02 seconds Orbital eigenvalues (eV) -40.47315 -39.91385 -38.03400 -36.94667 -34.99651 -34.06727 -32.01167 -31.68278 -30.63243 -30.14575 -29.03156 -27.32429 -25.97382 -25.60360 -22.93801 -22.79528 -21.78597 -21.36787 -21.02865 -19.70091 -18.87712 -17.89302 -17.51209 -16.68690 -16.36688 -16.08894 -15.96276 -15.30480 -14.90655 -14.77463 -14.64491 -14.30011 -14.17046 -13.95755 -13.56913 -13.41266 -13.12752 -13.04573 -12.78609 -12.68361 -12.67837 -12.43876 -12.07173 -11.91035 -11.84838 -11.48505 -11.42726 -11.36654 -11.19084 -10.82734 -10.65124 -10.48880 -10.41885 -10.20886 -10.03990 -9.89477 -9.70717 -9.60422 -9.51958 -9.47745 -9.14291 -9.05691 -8.50002 -7.03716 -6.04462 -4.04132 2.21247 2.38144 3.29828 3.40487 3.42777 3.47104 4.01581 4.25886 4.38146 4.43055 4.65322 4.70289 4.76207 4.78120 4.85427 5.00118 5.02845 5.05953 5.13667 5.26633 5.33313 5.34453 5.57086 5.67367 5.74859 5.80736 5.86078 5.88420 6.00704 6.02975 6.06887 6.08525 6.16928 6.24927 6.31037 6.36761 6.52106 6.60908 6.71049 6.73648 6.84459 6.93133 6.98041 7.00487 7.06425 7.17165 7.30337 7.32449 7.65282 7.79425 8.20733 8.26971 8.54029 9.01349 9.59292 11.08477 Molecular weight = 340.22amu Principal moments of inertia in cm(-1) A = 0.026150 B = 0.002243 C = 0.002149 Principal moments of inertia in units of 10**(-40)*gram-cm**2 A = 1070.498201 B =12479.469550 C =13023.244274 Net atomic charges, atomic populations, and dipole contributions calculated with CM2 Atom NO. Type Charge No. of electrons 1 C -0.141 4.141 2 C -0.123 4.123 3 C 0.181 3.819 4 H 0.037 0.963 5 C -0.539 4.539 6 H 0.079 0.921 7 C 0.072 3.928 8 N -0.899 5.899 9 Si 0.886 3.114 10 H -0.279 1.279 11 H -0.300 1.300 12 H -0.264 1.264 13 N -0.526 5.526 14 C 0.066 3.934 15 C -0.113 4.113 16 C 0.113 3.887 17 C -0.139 4.139 18 C 0.053 3.947 19 H 0.082 0.918 20 O -0.346 6.346 21 C 0.041 3.959 22 C 0.512 3.488 23 O -0.515 6.515 24 N -0.737 5.737 25 C 0.089 3.911 26 C 0.029 3.971 27 O -0.351 6.351 28 C -0.116 4.116 29 C 0.062 3.938 30 H 0.042 0.958 31 H 0.060 0.940 32 H 0.031 0.969 33 H 0.059 0.941 34 H 0.045 0.955 35 H 0.262 0.738 36 H 0.061 0.939 37 H 0.050 0.950 38 H 0.059 0.941 39 H 0.074 0.926 40 H 0.089 0.911 41 H 0.087 0.913 42 H 0.082 0.918 43 H 0.082 0.918 44 H 0.403 0.597 45 H 0.077 0.923 46 H 0.061 0.939 47 H 0.060 0.940 48 H 0.082 0.918 49 H 0.068 0.932 50 H 0.062 0.938 51 H 0.073 0.927 52 H 0.051 0.949 53 H 0.068 0.932 For the charges calculated by CM2: Dipole moment (debyes) X Y Z Total from point charges -6.425 15.591 -17.870 24.570 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.199 4.199 2 C -0.160 4.160 3 C 0.074 3.926 4 H 0.055 0.945 5 C -0.578 4.578 6 H 0.097 0.903 7 C -0.130 4.130 8 N -0.537 5.537 9 Si 0.717 3.283 10 H -0.207 1.207 11 H -0.227 1.227 12 H -0.189 1.189 13 N -0.250 5.250 14 C -0.063 4.063 15 C -0.151 4.151 16 C 0.074 3.926 17 C -0.176 4.176 18 C -0.006 4.006 19 H 0.099 0.901 20 O -0.264 6.264 21 C -0.039 4.039 22 C 0.297 3.703 23 O -0.388 6.388 24 N -0.392 5.392 25 C -0.055 4.055 26 C -0.048 4.048 27 O -0.270 6.270 28 C -0.154 4.154 29 C -0.066 4.066 30 H 0.061 0.939 31 H 0.079 0.921 32 H 0.050 0.950 33 H 0.078 0.922 34 H 0.063 0.937 35 H 0.072 0.928 36 H 0.080 0.920 37 H 0.069 0.931 38 H 0.078 0.922 39 H 0.092 0.908 40 H 0.107 0.893 41 H 0.106 0.894 42 H 0.100 0.900 43 H 0.100 0.900 44 H 0.236 0.764 45 H 0.096 0.904 46 H 0.079 0.921 47 H 0.079 0.921 48 H 0.100 0.900 49 H 0.086 0.914 50 H 0.080 0.920 51 H 0.091 0.909 52 H 0.070 0.930 53 H 0.086 0.914 Dipole moment (debyes) X Y Z Total from point charges -7.129 14.639 -18.253 24.460 hybrid contribution 0.872 0.230 1.098 1.421 sum -6.257 14.870 -17.155 23.549 Atomic orbital electron populations 1.21711 0.97387 0.99082 1.01670 1.21957 0.94015 0.99476 1.00593 1.19695 0.94376 0.90294 0.88275 0.94545 1.21675 0.90527 1.38082 1.07512 0.90314 1.24998 0.94677 0.95530 0.97844 1.69927 1.08074 1.45137 1.30583 0.86871 0.80477 0.79861 0.81095 1.20654 1.22748 1.18945 1.62322 1.58160 1.01479 1.03070 1.21769 0.99563 0.89417 0.95533 1.21906 1.00392 0.98387 0.94460 1.21915 0.80844 0.95996 0.93882 1.23227 0.98019 0.98445 0.97944 1.23979 0.96409 0.98977 0.81203 0.90055 1.88052 1.79751 1.36385 1.22186 1.22250 0.99853 0.89201 0.92565 1.20349 0.78382 0.81903 0.89646 1.90758 1.41650 1.39725 1.66695 1.45680 1.62306 1.22159 1.09024 1.21546 1.02442 1.02086 0.79426 1.23802 0.96791 0.95762 0.88490 1.88615 1.22190 1.77909 1.38259 1.21927 1.01185 0.93385 0.98904 1.21874 0.98548 0.97598 0.88622 0.93893 0.92140 0.94968 0.92194 0.93659 0.92819 0.92043 0.93146 0.92200 0.90775 0.89269 0.89436 0.89996 0.89980 0.76363 0.90431 0.92061 0.92089 0.90027 0.91406 0.91976 0.90874 0.93029 0.91387 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.14 -2.58 9.72 37.16 0.36 -2.22 16 2 C -0.12 -2.68 4.71 -26.73 -0.13 -2.81 16 3 C 0.18 4.24 1.35 -68.86 -0.09 4.15 16 4 H 0.04 0.90 7.12 -51.93 -0.37 0.53 16 5 C -0.54 -14.90 7.78 -34.44 -0.27 -15.17 16 6 H 0.08 2.41 6.09 -52.49 -0.32 2.09 16 7 C 0.07 2.07 5.45 -17.82 -0.10 1.97 16 8 N -0.90 -27.63 13.96 55.43 0.77 -26.86 16 9 Si 0.89 22.14 24.40 -169.99 -4.15 18.00 16 10 H -0.28 -6.85 6.84 56.52 0.39 -6.47 16 11 H -0.30 -7.54 7.11 56.52 0.40 -7.14 16 12 H -0.26 -6.34 6.53 56.52 0.37 -5.97 16 13 N -0.53 -10.64 4.11 -228.02 -0.94 -11.58 16 14 C 0.07 1.25 5.01 -3.81 -0.02 1.24 16 15 C -0.11 -1.68 4.94 -26.75 -0.13 -1.82 16 16 C 0.11 1.66 0.79 -90.12 -0.07 1.59 16 17 C -0.14 -1.67 5.82 -24.87 -0.14 -1.81 16 18 C 0.05 0.58 3.54 -24.45 -0.09 0.49 16 19 H 0.08 0.75 7.91 -51.93 -0.41 0.33 16 20 O -0.35 -4.09 9.70 -55.14 -0.53 -4.63 16 21 C 0.04 0.33 6.30 36.01 0.23 0.55 16 22 C 0.51 4.97 7.91 -10.99 -0.09 4.88 16 23 O -0.51 -7.62 16.50 -14.35 -0.24 -7.86 16 24 N -0.74 -4.21 5.71 -67.08 -0.38 -4.59 16 25 C 0.09 0.45 10.70 59.85 0.64 1.09 16 26 C 0.03 0.37 7.52 38.19 0.29 0.66 16 27 O -0.35 -5.83 9.79 -35.23 -0.35 -6.18 16 28 C -0.12 -1.83 4.95 -26.74 -0.13 -1.96 16 29 C 0.06 1.21 4.69 -3.81 -0.02 1.19 16 30 H 0.04 0.78 8.14 -51.93 -0.42 0.35 16 31 H 0.06 1.05 7.86 -51.93 -0.41 0.64 16 32 H 0.03 0.52 8.14 -51.93 -0.42 0.10 16 33 H 0.06 1.30 6.13 -51.93 -0.32 0.98 16 34 H 0.05 1.06 3.22 -51.92 -0.17 0.89 16 35 H 0.26 7.60 8.31 -40.82 -0.34 7.26 16 36 H 0.06 1.16 8.06 -51.93 -0.42 0.74 16 37 H 0.05 1.13 6.72 -51.93 -0.35 0.78 16 38 H 0.06 0.79 8.14 -51.92 -0.42 0.37 16 39 H 0.07 1.01 8.12 -51.93 -0.42 0.59 16 40 H 0.09 1.11 7.83 -51.93 -0.41 0.71 16 41 H 0.09 0.92 7.78 -51.93 -0.40 0.51 16 42 H 0.08 0.43 8.11 -51.93 -0.42 0.01 16 43 H 0.08 0.47 8.14 -51.93 -0.42 0.05 16 44 H 0.40 0.89 8.83 -40.82 -0.36 0.53 16 45 H 0.08 0.15 8.14 -51.93 -0.42 -0.28 16 46 H 0.06 0.38 8.14 -51.93 -0.42 -0.05 16 47 H 0.06 0.36 8.14 -51.93 -0.42 -0.06 16 48 H 0.08 0.90 8.14 -51.93 -0.42 0.48 16 49 H 0.07 0.90 8.14 -51.93 -0.42 0.48 16 50 H 0.06 0.94 8.14 -51.92 -0.42 0.52 16 51 H 0.07 1.06 8.14 -51.93 -0.42 0.64 16 52 H 0.05 1.19 6.32 -51.93 -0.33 0.87 16 53 H 0.07 1.33 5.96 -51.93 -0.31 1.02 16 LS Contribution 399.79 15.07 6.02 6.02 Total: -1.00 -35.33 399.79 -8.79 -44.12 By element: Atomic # 1 Polarization: 10.77 SS G_CDS: -9.25 Total: 1.53 kcal Atomic # 6 Polarization: -8.22 SS G_CDS: 0.24 Total: -7.97 kcal Atomic # 7 Polarization: -42.48 SS G_CDS: -0.55 Total: -43.03 kcal Atomic # 8 Polarization: -17.55 SS G_CDS: -1.12 Total: -18.66 kcal Atomic # 14 Polarization: 22.14 SS G_CDS: -4.15 Total: 18.00 kcal Total LS contribution 6.02 Total: 6.02 kcal Total: -35.33 -8.79 -44.12 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. ZINC001255167933.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 -67.978 kcal (2) G-P(sol) polarization free energy of solvation -35.331 kcal (3) G-ENP(sol) elect.-nuc.-pol. free energy of system -103.309 kcal (4) G-CDS(sol) cavity-dispersion-solvent structure free energy -8.793 kcal (5) G-P-CDS(sol) = G-P(sol) + G-CDS(sol) = (2) + (4) -44.124 kcal (6) G-S(sol) free energy of system = (1) + (5) -112.103 kcal 1SCF WAS SPECIFIED, SO GEOMETRY OPTIMIZATION WAS NOT USED Total computer time = 1.02 seconds