Wall clock time and date at job start Tue Mar 31 2020 05:50: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 * SM5.42R - SM5.42R CALCULATIONS WILL BE PERFORMED * - CHARGE MODEL 2 WILL BE USED (CM2A) * - THE SOLVENT IS WATER ******************************************************************************* 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.52992 * 1 3 3 H 1.09002 * 109.47202 * 2 1 4 4 C 1.50701 * 109.47010 * 119.99780 * 2 1 3 5 5 O 1.21283 * 120.00130 * 92.80407 * 4 2 1 6 6 N 1.34774 * 120.00352 * 272.80075 * 4 2 1 7 7 C 1.46504 * 119.99936 * 5.06712 * 6 4 2 8 8 C 1.50699 * 109.47194 * 270.00223 * 7 6 4 9 9 H 1.08002 * 119.99922 * 359.97438 * 8 7 6 10 10 C 1.37877 * 120.00137 * 180.02562 * 8 7 6 11 11 N 1.31514 * 120.00030 * 359.97438 * 10 8 7 12 12 Si 1.86805 * 119.99958 * 179.97438 * 10 8 7 13 13 H 1.48500 * 109.47067 * 359.97438 * 12 10 8 14 14 H 1.48501 * 109.47417 * 119.99793 * 12 10 8 15 15 H 1.48503 * 109.47024 * 239.99764 * 12 10 8 16 16 C 1.39923 * 120.00656 * 185.06601 * 6 4 2 17 17 C 1.38868 * 120.07462 * 173.01983 * 16 6 4 18 18 C 1.38118 * 119.93142 * 179.97438 * 17 16 6 19 19 C 1.38275 * 120.06600 * 359.97438 * 18 17 16 20 20 C 1.38282 * 120.14681 * 0.04436 * 19 18 17 21 21 C 1.38108 * 120.06687 * 359.95885 * 20 19 18 22 22 N 1.46898 * 109.47441 * 239.99936 * 2 1 3 23 23 C 1.46852 * 110.99928 * 182.38728 * 22 2 1 24 24 C 1.53039 * 109.52378 * 174.15866 * 23 22 2 25 25 C 1.53122 * 109.27571 * 59.17911 * 24 23 22 26 26 N 1.46896 * 109.57450 * 182.40947 * 25 24 23 27 27 C 1.46901 * 111.00143 * 294.71425 * 26 25 24 28 28 C 1.46903 * 110.99970 * 58.66992 * 26 25 24 29 29 C 1.53122 * 109.16234 * 302.34800 * 25 24 23 30 30 C 1.46848 * 111.00029 * 58.17977 * 22 2 1 31 31 H 1.09000 * 109.46770 * 53.74549 * 1 2 3 32 32 H 1.09006 * 109.47157 * 173.73976 * 1 2 3 33 33 H 1.08998 * 109.47632 * 293.74402 * 1 2 3 34 34 H 1.09001 * 109.47133 * 29.99385 * 7 6 4 35 35 H 1.08997 * 109.47706 * 149.99582 * 7 6 4 36 36 H 0.96997 * 119.99624 * 179.97438 * 11 10 8 37 37 H 1.08001 * 120.03462 * 359.69641 * 17 16 6 38 38 H 1.07995 * 119.96135 * 179.97438 * 18 17 16 39 39 H 1.08005 * 119.93078 * 179.97438 * 19 18 17 40 40 H 1.07999 * 119.97019 * 179.97438 * 20 19 18 41 41 H 1.08011 * 120.03561 * 180.02562 * 21 20 19 42 42 H 1.08995 * 109.46583 * 294.17680 * 23 22 2 43 43 H 1.09006 * 109.46101 * 54.14917 * 23 22 2 44 44 H 1.08998 * 109.50200 * 299.24322 * 24 23 22 45 45 H 1.08996 * 109.50291 * 179.11981 * 24 23 22 46 46 H 1.08997 * 109.52155 * 62.25367 * 25 24 23 47 47 H 1.09005 * 109.46553 * 60.00312 * 27 26 25 48 48 H 1.09003 * 109.47006 * 179.97438 * 27 26 25 49 49 H 1.08991 * 109.47235 * 300.00505 * 27 26 25 50 50 H 1.08998 * 109.46915 * 60.00502 * 28 26 25 51 51 H 1.08999 * 109.46962 * 180.02562 * 28 26 25 52 52 H 1.08994 * 109.47101 * 300.00118 * 28 26 25 53 53 H 1.09000 * 109.50094 * 297.71330 * 29 25 24 54 54 H 1.08996 * 109.50248 * 177.59896 * 29 25 24 55 55 H 1.08994 * 109.45835 * 305.85657 * 30 22 2 56 56 H 1.09000 * 109.46203 * 65.82700 * 30 22 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 6 1.5299 0.0000 0.0000 3 1 1.8933 1.0277 0.0000 4 6 2.0322 -0.7104 1.2305 5 8 2.2828 -1.8962 1.1863 6 7 2.2031 -0.0280 2.3801 7 6 1.7965 1.3769 2.4653 8 6 2.9507 2.2591 2.0645 9 1 3.8930 1.8175 1.7756 10 6 2.8041 3.6301 2.0676 11 7 1.6568 4.1679 2.4199 12 14 4.2350 4.7236 1.5715 13 1 5.4055 3.8790 1.2222 14 1 3.8472 5.5434 0.3956 15 1 4.5902 5.6177 2.7027 16 6 2.7688 -0.6612 3.4923 17 6 3.0792 0.0742 4.6286 18 6 3.6370 -0.5539 5.7249 19 6 3.8858 -1.9137 5.6920 20 6 3.5782 -2.6488 4.5619 21 6 3.0204 -2.0266 3.4623 22 7 2.0197 -0.6925 -1.1994 23 6 3.4865 -0.7617 -1.2050 24 6 3.9539 -1.6124 -2.3882 25 6 3.4488 -0.9868 -3.6913 26 7 3.8447 -1.8294 -4.8276 27 6 5.3056 -1.8563 -4.9793 28 6 3.3047 -3.1885 -4.6885 29 6 1.9220 -0.8817 -3.6397 30 6 1.5136 -0.0522 -2.4202 31 1 -0.3633 0.6077 -0.8287 32 1 -0.3634 -1.0216 -0.1121 33 1 -0.3634 0.4138 0.9406 34 1 0.9550 1.5530 1.7952 35 1 1.5007 1.6082 3.4885 36 1 1.5538 5.1324 2.4224 37 1 2.8888 1.1370 4.6534 38 1 3.8790 0.0176 6.6087 39 1 4.3225 -2.4026 6.5504 40 1 3.7742 -3.7106 4.5391 41 1 2.7799 -2.6014 2.5800 42 1 3.8972 0.2438 -1.2960 43 1 3.8324 -1.2123 -0.2747 44 1 3.5548 -2.6220 -2.2906 45 1 5.0431 -1.6506 -2.4012 46 1 3.8786 0.0079 -3.8084 47 1 5.7576 -2.2578 -4.0722 48 1 5.5713 -2.4876 -5.8273 49 1 5.6718 -0.8441 -5.1504 50 1 2.2164 -3.1447 -4.6485 51 1 3.6133 -3.7912 -5.5427 52 1 3.6843 -3.6381 -3.7710 53 1 1.4910 -1.8799 -3.5618 54 1 1.5597 -0.3981 -4.5468 55 1 0.4267 0.0125 -2.3725 56 1 1.9336 0.9500 -2.5056 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=WATER ZINC000095430158.mol2 56 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 05:50:24 Heat of formation + Delta-G solvation = 11.372391 kcal Electronic energy + Delta-G solvation = -34471.195952 eV Core-core repulsion = 30380.449270 eV Total energy + Delta-G solvation = -4090.746682 eV No. of doubly occupied orbitals = 69 Molecular weight (most abundant/longest-lived isotopes) = 359.231 amu Computer time = 1.54 seconds Orbital eigenvalues (eV) -42.05634 -40.62142 -39.31159 -37.68903 -36.97507 -34.67111 -34.10500 -32.06232 -31.73939 -31.16349 -30.07066 -29.12495 -28.52018 -26.53217 -25.43799 -25.01180 -23.88544 -23.53675 -22.82906 -21.98250 -20.31701 -20.27874 -18.82644 -18.49471 -17.99396 -17.45529 -16.90898 -16.59537 -16.44045 -16.26822 -15.94182 -15.71710 -15.51821 -15.37289 -15.10190 -14.83059 -14.73426 -14.47838 -14.36575 -14.24682 -14.08735 -13.80450 -13.53865 -13.31331 -13.18844 -13.10904 -12.97072 -12.84012 -12.78111 -12.63598 -12.60478 -12.59927 -12.34811 -12.27543 -12.02820 -11.96891 -11.92500 -11.83114 -11.57849 -11.25168 -10.90471 -10.70262 -10.11041 -9.77653 -9.33933 -9.21464 -8.78957 -8.27203 -6.43721 0.41618 0.57404 1.32660 1.36921 1.45778 1.63354 1.99809 2.33124 2.81934 3.00447 3.01615 3.04149 3.43370 3.74244 3.81685 4.01727 4.06687 4.13922 4.19181 4.20771 4.25808 4.32831 4.35238 4.41349 4.44544 4.49756 4.56846 4.59970 4.62981 4.64548 4.65519 4.69063 4.71563 4.75368 4.83422 4.90992 4.93555 4.95972 4.99619 5.01800 5.06407 5.07247 5.17916 5.23657 5.26219 5.36058 5.44255 5.59126 5.61993 5.70201 5.80193 5.89413 5.90674 6.05864 6.10890 6.13336 6.32775 6.42144 6.68823 7.16741 7.28374 8.76543 Molecular weight = 359.23amu Principal moments of inertia in cm(-1) A = 0.008081 B = 0.003767 C = 0.003032 Principal moments of inertia in units of 10**(-40)*gram-cm**2 A = 3464.016621 B = 7431.216657 C = 9233.233338 Net atomic charges, atomic populations, and dipole contributions calculated with CM2 Atom NO. Type Charge No. of electrons 1 C -0.131 4.131 2 C 0.070 3.930 3 H 0.063 0.937 4 C 0.520 3.480 5 O -0.537 6.537 6 N -0.543 5.543 7 C 0.244 3.756 8 C -0.542 4.542 9 H 0.043 0.957 10 C 0.033 3.967 11 N -0.910 5.910 12 Si 0.949 3.051 13 H -0.278 1.278 14 H -0.273 1.273 15 H -0.275 1.275 16 C 0.139 3.861 17 C -0.168 4.168 18 C -0.102 4.102 19 C -0.151 4.151 20 C -0.105 4.105 21 C -0.198 4.198 22 N -0.514 5.514 23 C 0.062 3.938 24 C -0.131 4.131 25 C 0.086 3.914 26 N -0.543 5.543 27 C 0.019 3.981 28 C 0.017 3.983 29 C -0.096 4.096 30 C 0.051 3.949 31 H 0.089 0.911 32 H 0.073 0.927 33 H 0.057 0.943 34 H 0.046 0.954 35 H 0.028 0.972 36 H 0.271 0.729 37 H 0.115 0.885 38 H 0.144 0.856 39 H 0.152 0.848 40 H 0.144 0.856 41 H 0.166 0.834 42 H 0.018 0.982 43 H 0.050 0.950 44 H 0.071 0.929 45 H 0.066 0.934 46 H 0.073 0.927 47 H 0.024 0.976 48 H 0.080 0.920 49 H 0.067 0.933 50 H 0.074 0.926 51 H 0.078 0.922 52 H 0.024 0.976 53 H 0.075 0.925 54 H 0.089 0.911 55 H 0.097 0.903 56 H 0.029 0.971 For the charges calculated by CM2: Dipole moment (debyes) X Y Z Total from point charges 1.666 -10.544 -8.492 13.640 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.189 4.189 2 C -0.041 4.041 3 H 0.081 0.919 4 C 0.311 3.689 5 O -0.417 6.417 6 N -0.263 5.263 7 C 0.124 3.876 8 C -0.581 4.581 9 H 0.061 0.939 10 C -0.168 4.168 11 N -0.548 5.548 12 Si 0.774 3.226 13 H -0.202 1.202 14 H -0.198 1.198 15 H -0.201 1.201 16 C 0.044 3.956 17 C -0.188 4.188 18 C -0.120 4.120 19 C -0.169 4.169 20 C -0.124 4.124 21 C -0.218 4.218 22 N -0.240 5.240 23 C -0.067 4.067 24 C -0.170 4.170 25 C -0.023 4.023 26 N -0.265 5.265 27 C -0.129 4.129 28 C -0.131 4.131 29 C -0.134 4.134 30 C -0.077 4.077 31 H 0.108 0.892 32 H 0.092 0.908 33 H 0.076 0.924 34 H 0.063 0.937 35 H 0.046 0.954 36 H 0.080 0.920 37 H 0.133 0.867 38 H 0.162 0.838 39 H 0.169 0.831 40 H 0.162 0.838 41 H 0.182 0.818 42 H 0.037 0.963 43 H 0.069 0.931 44 H 0.089 0.911 45 H 0.084 0.916 46 H 0.091 0.909 47 H 0.043 0.957 48 H 0.099 0.901 49 H 0.086 0.914 50 H 0.093 0.907 51 H 0.096 0.904 52 H 0.042 0.958 53 H 0.093 0.907 54 H 0.108 0.892 55 H 0.116 0.884 56 H 0.047 0.953 Dipole moment (debyes) X Y Z Total from point charges 1.272 -11.194 -9.015 14.429 hybrid contribution 1.463 1.714 0.997 2.464 sum 2.735 -9.479 -8.018 12.714 Atomic orbital electron populations 1.21761 0.90427 1.03492 1.03221 1.21539 0.98400 0.96052 0.88090 0.91913 1.20917 0.76756 0.86640 0.84562 1.90227 1.51479 1.15967 1.84015 1.46010 1.56377 1.13888 1.10006 1.19381 0.95876 0.74248 0.98112 1.21404 0.98096 0.91570 1.46998 0.93901 1.25487 0.97970 0.98770 0.94601 1.69983 1.25010 1.12426 1.47422 0.85644 0.79609 0.78866 0.78447 1.20249 1.19796 1.20054 1.18416 0.96750 0.92654 0.87774 1.21111 1.03834 1.00055 0.93826 1.21216 0.97753 0.95525 0.97555 1.21409 1.02429 0.93859 0.99244 1.21401 0.97599 1.00756 0.92599 1.22747 1.02969 0.93399 1.02643 1.61827 1.08889 1.52443 1.00845 1.21941 0.87755 0.97808 0.99199 1.22093 1.00156 1.00063 0.94693 1.21752 0.96593 0.96597 0.87366 1.62317 1.11684 1.12242 1.40281 1.22452 0.86964 1.02653 1.00817 1.22509 1.01167 0.88854 1.00547 1.21574 0.93895 1.01198 0.96696 1.22292 1.01466 0.95441 0.88485 0.89184 0.90800 0.92418 0.93653 0.95422 0.92005 0.86711 0.83781 0.83078 0.83779 0.81774 0.96342 0.93101 0.91087 0.91593 0.90895 0.95717 0.90082 0.91407 0.90706 0.90371 0.95765 0.90660 0.89231 0.88447 0.95254 Number of geometries 1 Number of calculations of the screened coulomb radii 1 The total number of SCF iterations 15. 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.13 -4.10 8.41 71.98 0.61 -3.49 16 2 C 0.07 2.61 1.63 44.17 0.07 2.68 16 3 H 0.06 2.67 6.40 -2.39 -0.02 2.65 16 4 C 0.52 22.03 5.71 87.66 0.50 22.53 16 5 O -0.54 -23.55 11.78 -3.04 -0.04 -23.58 16 6 N -0.54 -24.37 2.68 -635.31 -1.70 -26.07 16 7 C 0.24 11.73 4.10 85.63 0.35 12.08 16 8 C -0.54 -28.09 8.51 25.60 0.22 -27.87 16 9 H 0.04 2.34 7.81 -2.91 -0.02 2.32 16 10 C 0.03 1.72 5.46 84.65 0.46 2.18 16 11 N -0.91 -48.44 13.29 21.45 0.28 -48.15 16 12 Si 0.95 41.44 29.54 68.60 2.03 43.47 16 13 H -0.28 -12.28 7.11 99.48 0.71 -11.58 16 14 H -0.27 -11.51 7.11 99.48 0.71 -10.80 16 15 H -0.28 -11.75 7.11 99.48 0.71 -11.05 16 16 C 0.14 6.02 6.49 38.15 0.25 6.27 16 17 C -0.17 -6.83 9.15 22.39 0.20 -6.62 16 18 C -0.10 -3.31 10.04 22.25 0.22 -3.09 16 19 C -0.15 -4.33 10.04 22.29 0.22 -4.11 16 20 C -0.11 -3.33 10.04 22.25 0.22 -3.11 16 21 C -0.20 -7.81 9.05 22.39 0.20 -7.61 16 22 N -0.51 -16.71 4.26 -875.52 -3.73 -20.44 16 23 C 0.06 2.10 5.12 86.33 0.44 2.54 16 24 C -0.13 -3.53 4.69 30.65 0.14 -3.38 16 25 C 0.09 1.99 2.29 45.00 0.10 2.10 16 26 N -0.54 -11.11 5.01 -922.81 -4.62 -15.73 16 27 C 0.02 0.34 9.32 127.69 1.19 1.53 16 28 C 0.02 0.31 8.82 127.69 1.13 1.44 16 29 C -0.10 -2.10 5.24 30.65 0.16 -1.94 16 30 C 0.05 1.32 5.18 86.32 0.45 1.77 16 31 H 0.09 2.35 6.61 -2.39 -0.02 2.33 16 32 H 0.07 2.19 8.14 -2.38 -0.02 2.17 16 33 H 0.06 1.92 6.81 -2.39 -0.02 1.91 16 34 H 0.05 2.19 5.29 -2.39 -0.01 2.18 16 35 H 0.03 1.42 6.53 -2.39 -0.02 1.41 16 36 H 0.27 13.60 8.31 -92.71 -0.77 12.83 16 37 H 0.11 5.11 6.27 -2.91 -0.02 5.09 16 38 H 0.14 3.79 8.06 -2.91 -0.02 3.77 16 39 H 0.15 3.14 8.06 -2.91 -0.02 3.12 16 40 H 0.14 3.62 8.06 -2.91 -0.02 3.60 16 41 H 0.17 6.94 4.54 -61.68 -0.28 6.66 16 42 H 0.02 0.66 8.14 -2.39 -0.02 0.65 16 43 H 0.05 2.03 6.04 -2.38 -0.01 2.02 16 44 H 0.07 2.00 6.27 -2.39 -0.01 1.98 16 45 H 0.07 1.71 6.66 -2.39 -0.02 1.69 16 46 H 0.07 1.80 8.13 -2.39 -0.02 1.78 16 47 H 0.02 0.46 6.51 -2.38 -0.02 0.45 16 48 H 0.08 1.19 8.12 -2.39 -0.02 1.18 16 49 H 0.07 1.19 8.14 -2.39 -0.02 1.17 16 50 H 0.07 1.33 6.28 -2.39 -0.02 1.32 16 51 H 0.08 1.18 8.12 -2.39 -0.02 1.16 16 52 H 0.02 0.48 6.23 -2.39 -0.01 0.47 16 53 H 0.07 1.66 6.32 -2.39 -0.02 1.64 16 54 H 0.09 1.71 8.14 -2.39 -0.02 1.69 16 55 H 0.10 2.24 6.09 -2.39 -0.01 2.23 16 56 H 0.03 0.81 8.14 -2.39 -0.02 0.79 16 Total: -1.00 -59.76 415.35 -0.02 -59.78 By element: Atomic # 1 Polarization: 36.21 SS G_CDS: 0.61 Total: 36.82 kcal Atomic # 6 Polarization: -13.25 SS G_CDS: 7.15 Total: -6.10 kcal Atomic # 7 Polarization: -100.62 SS G_CDS: -9.77 Total: -110.39 kcal Atomic # 8 Polarization: -23.55 SS G_CDS: -0.04 Total: -23.58 kcal Atomic # 14 Polarization: 41.44 SS G_CDS: 2.03 Total: 43.47 kcal Total: -59.76 -0.02 -59.78 kcal The number of atoms in the molecule is 56 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. ZINC000095430158.mol2 56 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 71.155 kcal (2) G-P(sol) polarization free energy of solvation -59.759 kcal (3) G-ENP(sol) elect.-nuc.-pol. free energy of system 11.397 kcal (4) G-CDS(sol) cavity-dispersion-solvent structure free energy -0.024 kcal (5) G-P-CDS(sol) = G-P(sol) + G-CDS(sol) = (2) + (4) -59.783 kcal (6) G-S(sol) free energy of system = (1) + (5) 11.372 kcal 1SCF WAS SPECIFIED, SO GEOMETRY OPTIMIZATION WAS NOT USED Total computer time = 1.54 seconds