Wall clock time and date at job start Sat Apr 11 2020 06:11:15 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 O 1.42900 * 1 3 3 C 1.42907 * 114.00067 * 2 1 4 4 H 1.08997 * 109.57467 * 29.84560 * 3 2 1 5 5 C 1.53093 * 109.52270 * 269.68371 * 3 2 1 6 6 O 1.42902 * 109.41620 * 297.70650 * 5 3 2 7 7 C 1.42908 * 114.09665 * 298.83798 * 6 5 3 8 8 C 1.53111 * 109.48615 * 61.11187 * 7 6 5 9 9 C 1.53177 * 109.15629 * 302.52808 * 8 7 6 10 10 H 1.08993 * 109.54367 * 296.92994 * 9 8 7 11 11 C 1.52999 * 109.54256 * 176.75486 * 9 8 7 12 12 C 1.50703 * 109.46793 * 65.17681 * 11 9 8 13 13 O 1.21275 * 120.00385 * 359.97438 * 12 11 9 14 14 N 1.34787 * 119.99972 * 180.02562 * 12 11 9 15 15 C 1.46493 * 119.99823 * 180.02562 * 14 12 11 16 16 C 1.53006 * 109.47083 * 179.97438 * 15 14 12 17 17 H 1.09007 * 109.49500 * 304.93351 * 16 15 14 18 18 C 1.53034 * 109.49754 * 184.83603 * 16 15 14 19 19 C 1.53047 * 109.53794 * 178.68227 * 18 16 15 20 20 C 1.53190 * 109.31207 * 61.36243 * 19 18 16 21 21 N 1.46932 * 108.77700 * 305.36609 * 20 19 18 22 22 C 1.36931 * 120.63228 * 233.62896 * 21 20 19 23 23 H 1.08001 * 120.00275 * 151.81096 * 22 21 20 24 24 C 1.38812 * 119.99860 * 331.81371 * 22 21 20 25 25 N 1.31614 * 120.00428 * 343.07522 * 24 22 21 26 26 Si 1.86806 * 119.99729 * 163.07818 * 24 22 21 27 27 H 1.48501 * 109.47241 * 0.02562 * 26 24 22 28 28 H 1.48508 * 109.47030 * 119.99810 * 26 24 22 29 29 H 1.48497 * 109.46972 * 239.99616 * 26 24 22 30 30 C 1.46924 * 118.73553 * 53.61183 * 21 20 19 31 31 H 1.08995 * 109.46711 * 300.00178 * 1 2 3 32 32 H 1.08994 * 109.46780 * 59.99747 * 1 2 3 33 33 H 1.08998 * 109.47084 * 180.02562 * 1 2 3 34 34 H 1.09002 * 109.47491 * 57.68505 * 5 3 2 35 35 H 1.09000 * 109.48057 * 177.72172 * 5 3 2 36 36 H 1.08995 * 109.47989 * 181.14278 * 7 6 5 37 37 H 1.08996 * 109.48194 * 301.17064 * 7 6 5 38 38 H 1.08993 * 109.59735 * 182.57309 * 8 7 6 39 39 H 1.09008 * 109.44687 * 62.38299 * 8 7 6 40 40 H 1.09001 * 109.47301 * 305.17697 * 11 9 8 41 41 H 1.08996 * 109.47273 * 185.17500 * 11 9 8 42 42 H 0.97001 * 119.99448 * 359.97438 * 14 12 11 43 43 H 1.09001 * 109.47460 * 59.99465 * 15 14 12 44 44 H 1.09005 * 109.47105 * 299.99370 * 15 14 12 45 45 H 1.09004 * 109.46139 * 298.69920 * 18 16 15 46 46 H 1.09001 * 109.46380 * 58.66622 * 18 16 15 47 47 H 1.08992 * 109.49774 * 181.31914 * 19 18 16 48 48 H 1.09005 * 109.49553 * 301.38123 * 19 18 16 49 49 H 1.08999 * 109.58911 * 65.15332 * 20 19 18 50 50 H 1.08993 * 109.58743 * 185.57882 * 20 19 18 51 51 H 0.96997 * 120.00514 * 180.02562 * 25 24 22 52 52 H 1.08998 * 109.58575 * 186.60623 * 30 21 20 53 53 H 1.09001 * 109.59233 * 66.17854 * 30 21 20 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 8 1.4290 0.0000 0.0000 3 6 2.0103 1.3055 0.0000 4 1 1.3450 2.0014 0.5111 5 6 2.2256 1.7696 -1.4429 6 8 3.1433 0.8883 -2.0934 7 6 4.4381 0.8452 -1.4902 8 6 4.3135 0.3452 -0.0484 9 6 3.3607 1.2646 0.7219 10 1 3.7806 2.2693 0.7671 11 6 3.1663 0.7268 2.1410 12 6 4.4730 0.7980 2.8885 13 8 5.4623 1.2353 2.3400 14 7 4.5423 0.3747 4.1663 15 6 5.8124 0.4445 4.8930 16 6 5.6181 -0.0934 6.3122 17 1 4.8092 0.4507 6.7999 18 6 6.9114 0.0902 7.1095 19 6 6.7079 -0.4175 8.5389 20 6 6.3700 -1.9112 8.5027 21 7 5.2127 -2.1131 7.6201 22 6 4.1039 -2.7818 8.0655 23 1 3.1338 -2.5575 7.6473 24 6 4.2267 -3.7479 9.0546 25 7 5.4068 -4.2614 9.3301 26 14 2.7169 -4.3194 9.9946 27 1 1.5205 -3.5933 9.4981 28 1 2.5298 -5.7785 9.7908 29 1 2.9003 -4.0413 11.4417 30 6 5.2678 -1.5835 6.2508 31 1 -0.3632 0.5138 0.8899 32 1 -0.3633 0.5139 -0.8899 33 1 -0.3633 -1.0276 0.0005 34 1 1.2737 1.7597 -1.9738 35 1 2.6314 2.7813 -1.4426 36 1 5.0800 0.1694 -2.0552 37 1 4.8728 1.8447 -1.4906 38 1 5.2937 0.3549 0.4280 39 1 3.9180 -0.6706 -0.0495 40 1 2.8304 -0.3091 2.0942 41 1 2.4187 1.3282 2.6582 42 1 3.7511 0.0247 4.6048 43 1 6.1477 1.4806 4.9402 44 1 6.5605 -0.1563 4.3757 45 1 7.7140 -0.4742 6.6345 46 1 7.1761 1.1473 7.1339 47 1 7.6219 -0.2650 9.1127 48 1 5.8894 0.1300 9.0063 49 1 7.2238 -2.4699 8.1193 50 1 6.1272 -2.2554 9.5080 51 1 5.4926 -4.9368 10.0210 52 1 4.2979 -1.7113 5.7703 53 1 6.0318 -2.1153 5.6837 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 ZINC001650833607.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:11:15 Heat of formation + Delta-G solvation = -124.669681 kcal Electronic energy + Delta-G solvation = -30827.329301 eV Core-core repulsion = 26698.960959 eV Total energy + Delta-G solvation = -4128.368342 eV No. of doubly occupied orbitals = 66 Molecular weight (most abundant/longest-lived isotopes) = 340.220 amu Computer time = 0.92 seconds Orbital eigenvalues (eV) -40.49712 -39.34564 -37.81271 -37.21179 -34.96340 -34.21248 -33.01971 -30.81132 -30.59109 -30.44887 -28.94141 -26.86251 -25.57040 -25.26202 -23.64993 -23.45590 -22.11442 -20.98459 -20.43876 -20.03033 -19.04799 -17.23740 -16.88438 -16.71353 -16.57580 -16.06234 -15.84354 -15.56398 -15.24233 -15.12369 -14.82554 -14.42112 -14.16826 -13.92638 -13.78772 -13.46189 -13.22321 -13.01193 -12.80491 -12.68501 -12.47937 -12.34358 -12.24236 -12.18935 -11.72390 -11.59894 -11.30755 -11.20724 -11.04083 -10.96154 -10.74440 -10.55194 -10.50421 -10.27355 -10.23351 -10.20390 -9.96657 -9.80856 -9.57275 -9.44007 -9.29091 -8.89546 -8.49746 -6.69853 -5.81849 -3.25047 2.70986 3.27390 3.34132 3.38627 3.42008 3.43727 3.77512 4.36377 4.43122 4.45424 4.56885 4.65457 4.92519 4.95873 5.03507 5.08015 5.11845 5.12933 5.19854 5.26217 5.28189 5.32537 5.42032 5.43503 5.45072 5.46119 5.55975 5.60537 5.74124 5.80983 5.82763 5.94035 5.98550 6.03156 6.21409 6.24106 6.31768 6.35718 6.53967 6.57158 6.68474 6.75815 6.81117 6.83747 6.94560 7.07678 7.37501 7.60840 7.65035 7.75934 8.30606 8.41706 9.17931 9.81427 10.38427 11.26468 Molecular weight = 340.22amu Principal moments of inertia in cm(-1) A = 0.014033 B = 0.002383 C = 0.002124 Principal moments of inertia in units of 10**(-40)*gram-cm**2 A = 1994.853974 B =11746.201174 C =13180.612791 Net atomic charges, atomic populations, and dipole contributions calculated with CM2 Atom NO. Type Charge No. of electrons 1 C 0.029 3.971 2 O -0.367 6.367 3 C 0.057 3.943 4 H 0.070 0.930 5 C 0.042 3.958 6 O -0.375 6.375 7 C 0.062 3.938 8 C -0.152 4.152 9 C -0.066 4.066 10 H 0.089 0.911 11 C -0.134 4.134 12 C 0.512 3.488 13 O -0.547 6.547 14 N -0.724 5.724 15 C 0.127 3.873 16 C -0.118 4.118 17 H 0.072 0.928 18 C -0.100 4.100 19 C -0.142 4.142 20 C 0.169 3.831 21 N -0.592 5.592 22 C -0.254 4.254 23 H 0.065 0.935 24 C 0.003 3.997 25 N -0.901 5.901 26 Si 0.901 3.099 27 H -0.278 1.278 28 H -0.290 1.290 29 H -0.289 1.289 30 C 0.155 3.845 31 H 0.041 0.959 32 H 0.038 0.962 33 H 0.087 0.913 34 H 0.092 0.908 35 H 0.058 0.942 36 H 0.093 0.907 37 H 0.048 0.952 38 H 0.094 0.906 39 H 0.074 0.926 40 H 0.103 0.897 41 H 0.089 0.911 42 H 0.403 0.597 43 H 0.060 0.940 44 H 0.065 0.935 45 H 0.056 0.944 46 H 0.046 0.954 47 H 0.049 0.951 48 H 0.053 0.947 49 H 0.017 0.983 50 H 0.086 0.914 51 H 0.253 0.747 52 H 0.050 0.950 53 H 0.018 0.982 For the charges calculated by CM2: Dipole moment (debyes) X Y Z Total from point charges -6.099 13.842 -17.237 22.933 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.068 4.068 2 O -0.286 6.286 3 C -0.002 4.002 4 H 0.087 0.913 5 C -0.035 4.035 6 O -0.293 6.293 7 C -0.015 4.015 8 C -0.189 4.189 9 C -0.085 4.085 10 H 0.107 0.893 11 C -0.174 4.174 12 C 0.298 3.702 13 O -0.426 6.426 14 N -0.375 5.375 15 C 0.003 3.997 16 C -0.138 4.138 17 H 0.090 0.910 18 C -0.138 4.138 19 C -0.181 4.181 20 C 0.043 3.957 21 N -0.314 5.314 22 C -0.383 4.383 23 H 0.083 0.917 24 C -0.196 4.196 25 N -0.545 5.545 26 Si 0.731 3.269 27 H -0.203 1.203 28 H -0.217 1.217 29 H -0.216 1.216 30 C 0.028 3.972 31 H 0.059 0.941 32 H 0.057 0.943 33 H 0.105 0.895 34 H 0.110 0.890 35 H 0.076 0.924 36 H 0.111 0.889 37 H 0.066 0.934 38 H 0.112 0.888 39 H 0.093 0.907 40 H 0.121 0.879 41 H 0.108 0.892 42 H 0.238 0.762 43 H 0.079 0.921 44 H 0.083 0.917 45 H 0.075 0.925 46 H 0.065 0.935 47 H 0.068 0.932 48 H 0.072 0.928 49 H 0.036 0.964 50 H 0.104 0.896 51 H 0.064 0.936 52 H 0.068 0.932 53 H 0.036 0.964 Dipole moment (debyes) X Y Z Total from point charges -5.289 13.721 -17.653 22.976 hybrid contribution -1.308 -0.153 1.126 1.733 sum -6.597 13.568 -16.527 22.378 Atomic orbital electron populations 1.22814 0.81378 1.02972 0.99642 1.88007 1.19160 1.26515 1.94892 1.22878 0.95225 0.85897 0.96174 0.91252 1.23017 0.94278 0.93078 0.93123 1.87944 1.22422 1.65318 1.53639 1.22613 0.84420 1.00613 0.93851 1.22298 1.01472 0.99871 0.95287 1.21027 0.93623 1.00904 0.92972 0.89277 1.21653 0.95658 1.04760 0.95345 1.20627 0.90666 0.77144 0.81726 1.90637 1.37957 1.46335 1.67636 1.46128 1.13689 1.64588 1.13137 1.21441 0.85291 1.00668 0.92256 1.21953 0.99072 1.00035 0.92760 0.90956 1.21341 0.96324 0.99382 0.96726 1.22172 0.99865 1.01131 0.94920 1.20797 0.88692 0.89047 0.97144 1.44044 1.16502 1.63903 1.06949 1.19190 0.86358 1.13708 1.19061 0.91727 1.25052 0.96535 0.96342 1.01683 1.70041 1.22270 1.28847 1.33384 0.86245 0.82583 0.78165 0.79890 1.20269 1.21712 1.21557 1.20891 0.98358 0.89820 0.88160 0.94078 0.94291 0.89491 0.89002 0.92388 0.88870 0.93361 0.88755 0.90696 0.87866 0.89231 0.76225 0.92133 0.91664 0.92489 0.93495 0.93192 0.92786 0.96438 0.89581 0.93650 0.93171 0.96386 Number of geometries 1 Number of calculations of the screened coulomb radii 1 The total number of SCF iterations 16. 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.03 0.20 10.82 101.05 1.09 1.30 16 2 O -0.37 -3.66 9.08 -42.72 -0.39 -4.05 16 3 C 0.06 0.48 2.43 -26.59 -0.06 0.42 16 4 H 0.07 0.46 7.79 -51.93 -0.40 0.05 16 5 C 0.04 0.36 6.49 37.21 0.24 0.60 16 6 O -0.37 -4.20 10.70 -42.72 -0.46 -4.65 16 7 C 0.06 0.62 6.84 37.21 0.25 0.87 16 8 C -0.15 -1.68 4.97 -26.58 -0.13 -1.81 16 9 C -0.07 -0.65 1.75 -90.44 -0.16 -0.81 16 10 H 0.09 0.89 7.96 -51.93 -0.41 0.47 16 11 C -0.13 -1.32 4.93 -27.88 -0.14 -1.46 16 12 C 0.51 6.56 7.57 -10.98 -0.08 6.48 16 13 O -0.55 -8.59 13.99 5.56 0.08 -8.51 16 14 N -0.72 -9.22 5.33 -60.29 -0.32 -9.54 16 15 C 0.13 1.81 5.07 -4.04 -0.02 1.79 16 16 C -0.12 -1.89 2.38 -90.50 -0.21 -2.10 16 17 H 0.07 1.19 8.14 -51.92 -0.42 0.77 16 18 C -0.10 -1.61 5.17 -26.69 -0.14 -1.75 16 19 C -0.14 -2.63 6.08 -26.61 -0.16 -2.79 16 20 C 0.17 3.96 5.30 -3.71 -0.02 3.94 16 21 N -0.59 -14.37 2.99 -172.32 -0.52 -14.88 16 22 C -0.25 -6.77 9.58 -15.06 -0.14 -6.91 16 23 H 0.07 1.68 7.88 -52.49 -0.41 1.27 16 24 C 0.00 0.07 4.92 -17.43 -0.09 -0.01 16 25 N -0.90 -25.88 11.20 55.49 0.62 -25.26 16 26 Si 0.90 21.28 29.59 -169.99 -5.03 16.25 16 27 H -0.28 -6.36 7.11 56.52 0.40 -5.96 16 28 H -0.29 -6.94 7.11 56.52 0.40 -6.53 16 29 H -0.29 -6.85 7.11 56.52 0.40 -6.45 16 30 C 0.15 3.06 5.60 -3.71 -0.02 3.04 16 31 H 0.04 0.24 8.09 -51.93 -0.42 -0.18 16 32 H 0.04 0.23 7.86 -51.93 -0.41 -0.18 16 33 H 0.09 0.57 8.14 -51.93 -0.42 0.15 16 34 H 0.09 0.66 7.86 -51.93 -0.41 0.25 16 35 H 0.06 0.41 8.14 -51.93 -0.42 -0.01 16 36 H 0.09 0.85 8.14 -51.93 -0.42 0.43 16 37 H 0.05 0.44 8.14 -51.93 -0.42 0.02 16 38 H 0.09 1.20 6.05 -51.93 -0.31 0.88 16 39 H 0.07 0.91 7.84 -51.92 -0.41 0.51 16 40 H 0.10 1.04 7.77 -51.93 -0.40 0.64 16 41 H 0.09 0.71 8.14 -51.93 -0.42 0.28 16 42 H 0.40 4.64 7.70 -40.82 -0.31 4.32 16 43 H 0.06 0.85 8.14 -51.93 -0.42 0.43 16 44 H 0.07 1.00 8.14 -51.93 -0.42 0.57 16 45 H 0.06 0.93 8.14 -51.93 -0.42 0.51 16 46 H 0.05 0.67 8.14 -51.93 -0.42 0.24 16 47 H 0.05 0.84 8.14 -51.93 -0.42 0.42 16 48 H 0.05 1.02 8.14 -51.93 -0.42 0.60 16 49 H 0.02 0.43 8.14 -51.93 -0.42 0.01 16 50 H 0.09 2.27 6.06 -51.93 -0.31 1.95 16 51 H 0.25 7.02 8.31 -40.82 -0.34 6.69 16 52 H 0.05 0.95 7.13 -51.93 -0.37 0.58 16 53 H 0.02 0.38 8.14 -51.93 -0.42 -0.04 16 LS Contribution 406.38 15.07 6.12 6.12 Total: -1.00 -31.74 406.38 -9.33 -41.06 By element: Atomic # 1 Polarization: 12.33 SS G_CDS: -9.65 Total: 2.69 kcal Atomic # 6 Polarization: 0.57 SS G_CDS: 0.21 Total: 0.78 kcal Atomic # 7 Polarization: -49.47 SS G_CDS: -0.22 Total: -49.69 kcal Atomic # 8 Polarization: -16.44 SS G_CDS: -0.77 Total: -17.21 kcal Atomic # 14 Polarization: 21.28 SS G_CDS: -5.03 Total: 16.25 kcal Total LS contribution 6.12 Total: 6.12 kcal Total: -31.74 -9.33 -41.06 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. ZINC001650833607.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 -83.609 kcal (2) G-P(sol) polarization free energy of solvation -31.735 kcal (3) G-ENP(sol) elect.-nuc.-pol. free energy of system -115.344 kcal (4) G-CDS(sol) cavity-dispersion-solvent structure free energy -9.325 kcal (5) G-P-CDS(sol) = G-P(sol) + G-CDS(sol) = (2) + (4) -41.060 kcal (6) G-S(sol) free energy of system = (1) + (5) -124.670 kcal 1SCF WAS SPECIFIED, SO GEOMETRY OPTIMIZATION WAS NOT USED Total computer time = 0.92 seconds