Wall clock time and date at job start Tue Mar 31 2020 20:58:42 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.46504 * 109.47032 * 2 1 4 4 C 1.34775 * 119.99934 * 274.59999 * 3 2 1 5 5 O 1.21352 * 120.00052 * 186.96296 * 4 3 2 6 6 C 1.49455 * 120.00040 * 6.96869 * 4 3 2 7 7 O 1.21347 * 119.99992 * 7.95904 * 6 4 3 8 8 N 1.34772 * 119.99747 * 187.96007 * 6 4 3 9 9 C 1.46923 * 120.63526 * 7.55216 * 8 6 4 10 10 C 1.52766 * 109.01037 * 233.58504 * 9 8 6 11 11 C 1.53039 * 109.37868 * 305.58918 * 10 9 8 12 12 C 1.50708 * 109.44639 * 181.25502 * 11 10 9 13 13 H 1.07999 * 119.99815 * 239.93787 * 12 11 10 14 14 C 1.37874 * 120.00000 * 59.93376 * 12 11 10 15 15 N 1.31518 * 120.00150 * 0.02562 * 14 12 11 16 16 Si 1.86802 * 119.99833 * 180.02562 * 14 12 11 17 17 H 1.48500 * 109.47216 * 60.00159 * 16 14 12 18 18 H 1.48503 * 109.47048 * 180.02562 * 16 14 12 19 19 H 1.48498 * 109.47007 * 299.99955 * 16 14 12 20 20 C 1.53032 * 109.54408 * 61.24063 * 11 10 9 21 21 C 1.46926 * 120.62822 * 187.83243 * 8 6 4 22 22 C 1.46505 * 119.99871 * 94.60218 * 3 2 1 23 23 C 1.53784 * 113.61504 * 28.32042 * 22 3 2 24 24 C 1.53779 * 87.07846 * 89.11437 * 23 22 3 25 25 C 1.53002 * 113.61319 * 139.98196 * 24 23 22 26 26 C 1.52999 * 113.61622 * 270.79827 * 24 23 22 27 27 C 1.53779 * 87.17860 * 25.39083 * 24 23 22 28 28 H 1.09000 * 109.47063 * 60.00151 * 1 2 3 29 29 H 1.08993 * 109.47119 * 180.02562 * 1 2 3 30 30 H 1.09000 * 109.47047 * 300.00628 * 1 2 3 31 31 H 1.08995 * 109.47352 * 239.99708 * 2 1 3 32 32 H 1.09002 * 109.46597 * 119.99627 * 2 1 3 33 33 H 1.09001 * 109.58526 * 353.41862 * 9 8 6 34 34 H 1.08993 * 109.58580 * 113.84085 * 9 8 6 35 35 H 1.08994 * 109.56201 * 185.54853 * 10 9 8 36 36 H 1.09007 * 109.35800 * 65.49461 * 10 9 8 37 37 H 1.08997 * 109.46109 * 301.21170 * 11 10 9 38 38 H 0.96992 * 119.99815 * 180.02562 * 15 14 12 39 39 H 1.08997 * 109.50105 * 178.68326 * 20 11 10 40 40 H 1.09005 * 109.52747 * 58.59617 * 20 11 10 41 41 H 1.09001 * 109.58767 * 246.15841 * 21 8 6 42 42 H 1.09004 * 109.58672 * 6.43747 * 21 8 6 43 43 H 1.09003 * 112.84728 * 257.12216 * 22 3 2 44 44 H 1.08997 * 113.61442 * 334.56911 * 23 22 3 45 45 H 1.08998 * 113.68769 * 203.67071 * 23 22 3 46 46 H 1.08999 * 109.47212 * 181.52717 * 25 24 23 47 47 H 1.09003 * 109.47315 * 301.52615 * 25 24 23 48 48 H 1.09001 * 109.47430 * 61.52726 * 25 24 23 49 49 H 1.09004 * 109.46957 * 285.64791 * 26 24 23 50 50 H 1.08994 * 109.47159 * 45.64656 * 26 24 23 51 51 H 1.09005 * 109.47412 * 165.64676 * 26 24 23 52 52 H 1.09001 * 113.61491 * 89.15336 * 27 24 23 53 53 H 1.09003 * 113.61890 * 220.06225 * 27 24 23 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 2.1547 2.0478 -1.1634 5 8 2.6764 3.1434 -1.1775 6 6 1.6633 1.4377 -2.4362 7 8 1.0147 0.4125 -2.4068 8 7 1.9406 2.0261 -3.6166 9 6 2.8592 3.1707 -3.6850 10 6 2.1637 4.3247 -4.4050 11 6 1.6574 3.8423 -5.7663 12 6 0.9981 4.9876 -6.4909 13 1 -0.0405 4.9181 -6.7787 14 6 1.7208 6.1239 -6.7862 15 7 2.9857 6.2084 -6.4360 16 14 0.9034 7.5437 -7.6837 17 1 -0.2403 8.0447 -6.8799 18 1 1.8881 8.6378 -7.8800 19 1 0.4127 7.0770 -9.0053 20 6 0.6426 2.7147 -5.5648 21 6 1.3252 1.5390 -4.8587 22 6 2.3585 2.0379 1.2647 23 6 2.8186 1.0626 2.3610 24 6 4.2453 1.2495 1.8187 25 6 5.3101 1.3412 2.9136 26 6 4.6149 0.2599 0.7119 27 6 3.7841 2.6187 1.2919 28 1 -0.3633 0.5138 -0.8900 29 1 -0.3633 -1.0276 0.0005 30 1 -0.3633 0.5139 0.8899 31 1 1.8934 -0.5138 -0.8899 32 1 1.8933 -0.5138 0.8901 33 1 3.1309 3.4820 -2.6763 34 1 3.7558 2.8872 -4.2360 35 1 2.8688 5.1434 -4.5482 36 1 1.3216 4.6698 -3.8050 37 1 2.4967 3.4745 -6.3564 38 1 3.4940 7.0080 -6.6435 39 1 0.2650 2.3876 -6.5336 40 1 -0.1855 3.0746 -4.9542 41 1 2.0942 1.1197 -5.5074 42 1 0.5869 0.7721 -4.6244 43 1 1.5926 2.7388 1.5970 44 1 2.4328 0.0509 2.2354 45 1 2.6824 1.4450 3.3726 46 1 6.2874 1.4995 2.4577 47 1 5.0790 2.1751 3.5763 48 1 5.3225 0.4142 3.4869 49 1 4.7927 -0.7232 1.1481 50 1 3.7976 0.1963 -0.0064 51 1 5.5178 0.6014 0.2057 52 1 3.9120 3.4285 2.0102 53 1 4.1766 2.8614 0.3044 RHF calculation, no. of doubly occupied orbitals= 65 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 ZINC001261368346.mol2 53 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 20:58:42 Heat of formation + Delta-G solvation = -19.664737 kcal Electronic energy + Delta-G solvation = -31646.981583 eV Core-core repulsion = 27713.622649 eV Total energy + Delta-G solvation = -3933.358934 eV No. of doubly occupied orbitals = 65 Molecular weight (most abundant/longest-lived isotopes) = 336.220 amu Computer time = 0.88 seconds Orbital eigenvalues (eV) -41.45283 -39.80483 -37.90657 -37.20743 -35.30053 -33.38203 -31.24877 -31.16301 -30.02749 -29.73928 -27.62723 -26.60822 -26.10966 -24.86291 -23.31257 -22.80056 -22.32021 -21.09777 -19.69081 -18.63498 -18.10441 -17.22517 -16.61766 -16.45829 -16.12600 -15.94851 -15.10604 -14.94125 -14.51065 -14.33849 -14.15696 -14.05139 -13.78821 -13.52472 -13.11243 -12.93111 -12.88709 -12.80011 -12.72355 -12.61947 -12.32737 -12.05764 -11.98333 -11.77786 -11.59998 -11.38034 -11.21612 -10.98160 -10.91595 -10.81618 -10.71261 -10.52543 -10.49189 -10.30185 -10.18258 -10.01333 -9.57684 -9.17754 -9.11180 -8.73197 -8.44276 -8.25084 -8.02976 -5.97859 -4.04360 1.65994 2.43724 3.32623 3.37384 3.39956 3.41927 3.60689 4.17717 4.27143 4.43815 4.44452 4.58948 4.62796 4.86791 5.02004 5.26170 5.31701 5.32983 5.34767 5.39430 5.44349 5.53620 5.62418 5.64040 5.66966 5.69569 5.78575 5.84620 5.86251 5.89992 5.95104 6.00339 6.07510 6.10098 6.14999 6.21042 6.26057 6.35784 6.47962 6.52033 6.65165 6.75256 6.81874 7.01097 7.14549 7.21784 7.34426 7.69999 7.82278 7.93864 7.99928 8.12341 8.31266 8.45263 8.98497 9.60236 11.07501 Molecular weight = 336.22amu Principal moments of inertia in cm(-1) A = 0.017548 B = 0.003084 C = 0.002815 Principal moments of inertia in units of 10**(-40)*gram-cm**2 A = 1595.267355 B = 9075.686460 C = 9945.204462 Net atomic charges, atomic populations, and dipole contributions calculated with CM2 Atom NO. Type Charge No. of electrons 1 C -0.162 4.162 2 C 0.102 3.898 3 N -0.570 5.570 4 C 0.510 3.490 5 O -0.496 6.496 6 C 0.497 3.503 7 O -0.535 6.535 8 N -0.560 5.560 9 C 0.048 3.952 10 C -0.116 4.116 11 C 0.046 3.954 12 C -0.514 4.514 13 H 0.053 0.947 14 C 0.062 3.938 15 N -0.890 5.890 16 Si 0.890 3.110 17 H -0.278 1.278 18 H -0.287 1.287 19 H -0.278 1.278 20 C -0.123 4.123 21 C 0.109 3.891 22 C 0.121 3.879 23 C -0.140 4.140 24 C -0.088 4.088 25 C -0.134 4.134 26 C -0.139 4.139 27 C -0.137 4.137 28 H 0.074 0.926 29 H 0.061 0.939 30 H 0.045 0.955 31 H 0.097 0.903 32 H 0.078 0.922 33 H 0.138 0.862 34 H 0.055 0.945 35 H 0.104 0.896 36 H 0.036 0.964 37 H 0.026 0.974 38 H 0.261 0.739 39 H 0.062 0.938 40 H 0.052 0.948 41 H 0.065 0.935 42 H 0.079 0.921 43 H 0.082 0.918 44 H 0.080 0.920 45 H 0.094 0.906 46 H 0.055 0.945 47 H 0.055 0.945 48 H 0.053 0.947 49 H 0.054 0.946 50 H 0.067 0.933 51 H 0.057 0.943 52 H 0.089 0.911 53 H 0.091 0.909 For the charges calculated by CM2: Dipole moment (debyes) X Y Z Total from point charges 2.028 -12.990 18.325 22.554 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.220 4.220 2 C -0.020 4.020 3 N -0.304 5.304 4 C 0.294 3.706 5 O -0.372 6.372 6 C 0.279 3.721 7 O -0.415 6.415 8 N -0.291 5.291 9 C -0.072 4.072 10 C -0.155 4.155 11 C 0.028 3.972 12 C -0.552 4.552 13 H 0.072 0.928 14 C -0.139 4.139 15 N -0.529 5.529 16 Si 0.719 3.281 17 H -0.204 1.204 18 H -0.213 1.213 19 H -0.204 1.204 20 C -0.161 4.161 21 C -0.013 4.013 22 C 0.017 3.983 23 C -0.177 4.177 24 C -0.088 4.088 25 C -0.191 4.191 26 C -0.196 4.196 27 C -0.175 4.175 28 H 0.093 0.907 29 H 0.080 0.920 30 H 0.064 0.936 31 H 0.115 0.885 32 H 0.096 0.904 33 H 0.154 0.846 34 H 0.073 0.927 35 H 0.122 0.878 36 H 0.055 0.945 37 H 0.044 0.956 38 H 0.071 0.929 39 H 0.080 0.920 40 H 0.070 0.930 41 H 0.083 0.917 42 H 0.098 0.902 43 H 0.100 0.900 44 H 0.098 0.902 45 H 0.113 0.887 46 H 0.074 0.926 47 H 0.074 0.926 48 H 0.072 0.928 49 H 0.074 0.926 50 H 0.086 0.914 51 H 0.076 0.924 52 H 0.107 0.893 53 H 0.109 0.891 Dipole moment (debyes) X Y Z Total from point charges 2.137 -12.804 18.294 22.431 hybrid contribution -0.510 0.920 -0.599 1.210 sum 1.627 -11.884 17.695 21.377 Atomic orbital electron populations 1.22158 0.95266 1.02174 1.02384 1.21940 0.95422 0.81163 1.03436 1.47614 1.63476 1.14575 1.04784 1.20426 0.79572 0.85582 0.84997 1.90401 1.41739 1.20614 1.84443 1.19856 0.82261 0.84959 0.85004 1.90654 1.41165 1.23084 1.86602 1.48225 1.43342 1.31974 1.05588 1.23311 0.91112 0.85980 1.06817 1.21835 1.00354 0.97050 0.96236 1.18677 0.93510 0.91958 0.93069 1.21151 0.97113 1.00509 1.36395 0.92845 1.24957 0.94840 0.98360 0.95772 1.69837 1.07605 1.28511 1.46911 0.86685 0.79166 0.82033 0.80258 1.20358 1.21308 1.20366 1.22013 0.99559 0.95803 0.98688 1.21864 0.97470 0.97791 0.84194 1.21948 0.97277 0.95443 0.83594 1.23458 0.96165 1.00098 0.97965 1.21824 0.94222 0.96610 0.96106 1.21709 0.98421 1.01392 0.97607 1.21794 1.02252 0.97807 0.97782 1.23440 0.94625 0.96036 1.03405 0.90681 0.91971 0.93625 0.88543 0.90422 0.84607 0.92699 0.87829 0.94510 0.95553 0.92936 0.91954 0.92967 0.91681 0.90241 0.89986 0.90208 0.88740 0.92581 0.92565 0.92758 0.92646 0.91418 0.92429 0.89263 0.89059 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.16 -1.79 9.62 37.16 0.36 -1.43 16 2 C 0.10 1.10 3.82 -4.04 -0.02 1.08 16 3 N -0.57 -7.06 1.38 -165.57 -0.23 -7.29 16 4 C 0.51 8.06 7.24 -11.58 -0.08 7.98 16 5 O -0.50 -8.78 11.05 5.49 0.06 -8.72 16 6 C 0.50 8.57 7.47 -11.58 -0.09 8.48 16 7 O -0.54 -9.12 10.71 5.50 0.06 -9.06 16 8 N -0.56 -10.42 2.98 -173.54 -0.52 -10.94 16 9 C 0.05 0.96 6.23 -3.93 -0.02 0.93 16 10 C -0.12 -2.75 5.14 -26.83 -0.14 -2.89 16 11 C 0.05 1.13 2.07 -91.73 -0.19 0.94 16 12 C -0.51 -13.99 8.56 -34.44 -0.29 -14.29 16 13 H 0.05 1.44 7.74 -52.49 -0.41 1.04 16 14 C 0.06 1.70 5.30 -17.82 -0.09 1.61 16 15 N -0.89 -25.10 11.31 55.43 0.63 -24.47 16 16 Si 0.89 21.29 29.54 -169.99 -5.02 16.27 16 17 H -0.28 -6.62 7.11 56.52 0.40 -6.22 16 18 H -0.29 -6.82 7.11 56.52 0.40 -6.42 16 19 H -0.28 -6.58 7.11 56.52 0.40 -6.17 16 20 C -0.12 -2.53 5.30 -26.61 -0.14 -2.67 16 21 C 0.11 1.97 6.41 -3.71 -0.02 1.94 16 22 C 0.12 1.28 4.34 -67.02 -0.29 0.99 16 23 C -0.14 -1.10 6.64 -25.92 -0.17 -1.27 16 24 C -0.09 -0.77 1.39 -153.70 -0.21 -0.99 16 25 C -0.13 -1.00 9.18 37.16 0.34 -0.66 16 26 C -0.14 -1.33 8.85 37.16 0.33 -1.00 16 27 C -0.14 -1.47 6.45 -25.83 -0.17 -1.64 16 28 H 0.07 1.05 6.13 -51.93 -0.32 0.73 16 29 H 0.06 0.60 8.14 -51.93 -0.42 0.18 16 30 H 0.04 0.46 8.14 -51.93 -0.42 0.04 16 31 H 0.10 1.18 6.03 -51.93 -0.31 0.87 16 32 H 0.08 0.66 5.08 -51.93 -0.26 0.39 16 33 H 0.14 2.70 3.83 -56.51 -0.22 2.48 16 34 H 0.05 1.09 8.14 -51.93 -0.42 0.67 16 35 H 0.10 2.73 6.03 -51.93 -0.31 2.42 16 36 H 0.04 0.90 8.14 -51.92 -0.42 0.48 16 37 H 0.03 0.70 8.14 -51.93 -0.42 0.27 16 38 H 0.26 7.07 8.31 -40.82 -0.34 6.73 16 39 H 0.06 1.23 8.14 -51.93 -0.42 0.80 16 40 H 0.05 1.11 8.14 -51.93 -0.42 0.69 16 41 H 0.06 1.12 8.14 -51.93 -0.42 0.70 16 42 H 0.08 1.36 7.06 -51.93 -0.37 0.99 16 43 H 0.08 0.88 8.14 -51.93 -0.42 0.46 16 44 H 0.08 0.59 5.48 -51.93 -0.28 0.30 16 45 H 0.09 0.61 8.10 -51.93 -0.42 0.19 16 46 H 0.06 0.43 8.14 -51.93 -0.42 0.01 16 47 H 0.06 0.39 7.98 -51.93 -0.41 -0.03 16 48 H 0.05 0.36 8.14 -51.93 -0.42 -0.06 16 49 H 0.05 0.45 8.14 -51.93 -0.42 0.03 16 50 H 0.07 0.73 4.24 -51.93 -0.22 0.51 16 51 H 0.06 0.58 8.14 -51.93 -0.42 0.16 16 52 H 0.09 0.86 7.99 -51.93 -0.41 0.44 16 53 H 0.09 1.24 6.14 -51.93 -0.32 0.92 16 LS Contribution 388.28 15.07 5.85 5.85 Total: -1.00 -28.66 388.28 -8.98 -37.64 By element: Atomic # 1 Polarization: 12.51 SS G_CDS: -8.90 Total: 3.60 kcal Atomic # 6 Polarization: -1.98 SS G_CDS: -0.91 Total: -2.88 kcal Atomic # 7 Polarization: -42.58 SS G_CDS: -0.12 Total: -42.70 kcal Atomic # 8 Polarization: -17.89 SS G_CDS: 0.12 Total: -17.77 kcal Atomic # 14 Polarization: 21.29 SS G_CDS: -5.02 Total: 16.27 kcal Total LS contribution 5.85 Total: 5.85 kcal Total: -28.66 -8.98 -37.64 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. ZINC001261368346.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 17.973 kcal (2) G-P(sol) polarization free energy of solvation -28.659 kcal (3) G-ENP(sol) elect.-nuc.-pol. free energy of system -10.686 kcal (4) G-CDS(sol) cavity-dispersion-solvent structure free energy -8.979 kcal (5) G-P-CDS(sol) = G-P(sol) + G-CDS(sol) = (2) + (4) -37.638 kcal (6) G-S(sol) free energy of system = (1) + (5) -19.665 kcal 1SCF WAS SPECIFIED, SO GEOMETRY OPTIMIZATION WAS NOT USED Total computer time = 0.88 seconds