Wall clock time and date at job start Tue Mar 31 2020 13:56:37 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 C 1.53002 * 109.47115 * 2 1 4 4 C 1.50696 * 109.47270 * 120.00123 * 2 1 3 5 5 O 1.21279 * 120.00150 * 264.99269 * 4 2 1 6 6 N 1.34777 * 119.99788 * 85.00335 * 4 2 1 7 7 C 1.46495 * 120.00118 * 179.97438 * 6 4 2 8 8 C 1.53039 * 109.46476 * 204.99626 * 7 6 4 9 9 C 1.53192 * 109.31624 * 178.61195 * 8 7 6 10 10 N 1.46925 * 108.89163 * 54.56131 * 9 8 7 11 11 C 1.36940 * 120.62821 * 126.46045 * 10 9 8 12 12 H 1.08000 * 119.99751 * 331.92858 * 11 10 9 13 13 C 1.38810 * 120.00078 * 151.92771 * 11 10 9 14 14 N 1.31619 * 119.99689 * 343.17980 * 13 11 10 15 15 Si 1.86795 * 119.99988 * 163.18147 * 13 11 10 16 16 H 1.48509 * 109.47304 * 59.99932 * 15 13 11 17 17 H 1.48501 * 109.47440 * 179.97438 * 15 13 11 18 18 H 1.48501 * 109.47535 * 299.99868 * 15 13 11 19 19 C 1.46954 * 118.71946 * 306.64559 * 10 9 8 20 20 C 1.53044 * 109.46278 * 84.92375 * 7 6 4 21 21 N 1.46900 * 109.46921 * 240.00710 * 2 1 3 22 22 C 1.46850 * 110.99956 * 184.20686 * 21 2 1 23 23 C 1.53041 * 109.52029 * 174.15286 * 22 21 2 24 24 C 1.53118 * 109.27478 * 59.18500 * 23 22 21 25 25 C 1.53125 * 109.16051 * 302.33927 * 24 23 22 26 26 C 1.46849 * 111.00309 * 59.99192 * 21 2 1 27 27 H 1.09004 * 109.46602 * 299.99717 * 1 2 3 28 28 H 1.08997 * 109.47194 * 59.99631 * 1 2 3 29 29 H 1.09001 * 109.46953 * 180.02562 * 1 2 3 30 30 H 1.08994 * 109.47360 * 180.02562 * 3 2 1 31 31 H 1.09003 * 109.47143 * 300.00558 * 3 2 1 32 32 H 1.09004 * 109.47034 * 60.00078 * 3 2 1 33 33 H 0.96996 * 119.99854 * 359.97438 * 6 4 2 34 34 H 1.09009 * 109.46016 * 324.97414 * 7 6 4 35 35 H 1.09005 * 109.49298 * 58.66180 * 8 7 6 36 36 H 1.08998 * 109.49763 * 298.56780 * 8 7 6 37 37 H 1.08997 * 109.58838 * 174.42439 * 9 8 7 38 38 H 1.09002 * 109.58672 * 294.71002 * 9 8 7 39 39 H 0.97000 * 119.99637 * 180.02562 * 14 13 11 40 40 H 1.09000 * 109.47125 * 293.56762 * 19 10 9 41 41 H 1.09000 * 109.71403 * 173.28134 * 19 10 9 42 42 H 1.09000 * 109.50000 * 301.34377 * 20 7 6 43 43 H 1.09010 * 109.49012 * 61.40617 * 20 7 6 44 44 H 1.08988 * 109.46100 * 294.17163 * 22 21 2 45 45 H 1.09004 * 109.46033 * 54.14623 * 22 21 2 46 46 H 1.09001 * 109.50725 * 299.24329 * 23 22 21 47 47 H 1.09001 * 109.49962 * 179.12072 * 23 22 21 48 48 H 1.09004 * 109.52467 * 182.43714 * 24 23 22 49 49 H 1.09001 * 109.52875 * 62.24303 * 24 23 22 50 50 H 1.08998 * 109.50283 * 297.71659 * 25 24 23 51 51 H 1.08996 * 109.50200 * 177.59733 * 25 24 23 52 52 H 1.08995 * 109.46066 * 305.85956 * 26 21 2 53 53 H 1.09003 * 109.45516 * 65.82728 * 26 21 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.5300 0.0000 0.0000 3 6 2.0400 1.4425 0.0000 4 6 2.0324 -0.7104 1.2304 5 8 2.3210 -0.0749 2.2222 6 7 2.1612 -2.0520 1.2285 7 6 2.6500 -2.7426 2.4244 8 6 3.2851 -4.0750 2.0200 9 6 3.8260 -4.7783 3.2689 10 7 2.7419 -4.9007 4.2530 11 6 2.4253 -6.1181 4.7944 12 1 2.6075 -7.0225 4.2330 13 6 1.8703 -6.1895 6.0647 14 7 1.9487 -5.1562 6.8761 15 14 1.0123 -7.7522 6.6225 16 1 1.9728 -8.8843 6.5864 17 1 0.5103 -7.5736 8.0086 18 1 -0.1279 -8.0415 5.7161 19 6 2.0014 -3.6924 4.6418 20 6 1.4810 -3.0013 3.3776 21 7 2.0196 -0.6924 -1.1995 22 6 3.4841 -0.7997 -1.1843 23 6 3.9469 -1.6499 -2.3696 24 6 3.4779 -0.9978 -3.6732 25 6 1.9537 -0.8535 -3.6429 26 6 1.5486 -0.0265 -2.4207 27 1 -0.3633 0.5138 0.8901 28 1 -0.3633 0.5139 -0.8899 29 1 -0.3633 -1.0277 0.0005 30 1 3.1300 1.4426 -0.0005 31 1 1.6766 1.9564 -0.8900 32 1 1.6767 1.9564 0.8900 33 1 1.9301 -2.5603 0.4353 34 1 3.3945 -2.1221 2.9233 35 1 4.1028 -3.8915 1.3229 36 1 2.5348 -4.7059 1.5434 37 1 4.1901 -5.7701 3.0009 38 1 4.6409 -4.1924 3.6941 39 1 1.5612 -5.2063 7.7639 40 1 2.6658 -3.0167 5.1804 41 1 1.1617 -3.9655 5.2809 42 1 1.0145 -2.0536 3.6468 43 1 0.7473 -3.6419 2.8880 44 1 3.9220 0.1957 -1.2584 45 1 3.8042 -1.2689 -0.2539 46 1 3.5204 -2.6498 -2.2889 47 1 5.0349 -1.7164 -2.3669 48 1 3.7679 -1.6225 -4.5180 49 1 3.9350 -0.0133 -3.7731 50 1 1.4959 -1.8408 -3.5822 51 1 1.6177 -0.3511 -4.5499 52 1 0.4631 0.0660 -2.3883 53 1 1.9957 0.9653 -2.4889 RHF calculation, no. of doubly occupied orbitals= 63 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 ZINC001002814521.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 13:56:37 Heat of formation + Delta-G solvation = -21.298317 kcal Electronic energy + Delta-G solvation = -29725.535407 eV Core-core repulsion = 26018.050438 eV Total energy + Delta-G solvation = -3707.484969 eV No. of doubly occupied orbitals = 63 Molecular weight (most abundant/longest-lived isotopes) = 323.231 amu Computer time = 0.91 seconds Orbital eigenvalues (eV) -40.77716 -38.72065 -37.29158 -35.12968 -34.35911 -32.06972 -31.24159 -30.54058 -29.98427 -29.40123 -27.49284 -25.49352 -24.55072 -23.49135 -22.52277 -22.16990 -21.24950 -20.08715 -19.46283 -18.81213 -17.48795 -16.74127 -16.02867 -15.60074 -15.06485 -14.81671 -14.71310 -14.57082 -14.28321 -14.13343 -13.83248 -13.74910 -13.46172 -13.39809 -12.98756 -12.83556 -12.57111 -12.45359 -12.31858 -12.13048 -12.10334 -11.91679 -11.57736 -11.43447 -11.20014 -11.08434 -10.98651 -10.88501 -10.69788 -10.67881 -10.52708 -10.34400 -10.22303 -10.01602 -9.72458 -9.70855 -9.33043 -8.60966 -8.57875 -8.36863 -6.64821 -5.75541 -3.18617 2.79905 3.39565 3.43622 3.46316 3.47361 3.54515 4.48369 4.50845 4.51995 4.63211 4.70145 4.87300 4.96248 5.00446 5.09792 5.18998 5.31642 5.35625 5.38996 5.41296 5.44661 5.50121 5.51439 5.55591 5.60008 5.65605 5.68884 5.71012 5.75096 5.87006 5.88069 5.89341 5.98982 6.03162 6.13019 6.14182 6.29249 6.43524 6.50259 6.66049 6.69298 6.71704 6.81058 6.82359 6.97052 7.03135 7.41606 7.71229 7.82101 7.87444 8.48009 8.52237 9.24094 9.87489 10.44514 11.33065 Molecular weight = 323.23amu Principal moments of inertia in cm(-1) A = 0.020821 B = 0.003029 C = 0.002862 Principal moments of inertia in units of 10**(-40)*gram-cm**2 A = 1344.458528 B = 9242.741485 C = 9781.610806 Net atomic charges, atomic populations, and dipole contributions calculated with CM2 Atom NO. Type Charge No. of electrons 1 C -0.137 4.137 2 C 0.090 3.910 3 C -0.166 4.166 4 C 0.525 3.475 5 O -0.543 6.543 6 N -0.712 5.712 7 C 0.149 3.851 8 C -0.137 4.137 9 C 0.149 3.851 10 N -0.592 5.592 11 C -0.255 4.255 12 H 0.065 0.935 13 C 0.003 3.997 14 N -0.902 5.902 15 Si 0.903 3.097 16 H -0.283 1.283 17 H -0.291 1.291 18 H -0.284 1.284 19 C 0.174 3.826 20 C -0.145 4.145 21 N -0.524 5.524 22 C 0.043 3.957 23 C -0.115 4.115 24 C -0.121 4.121 25 C -0.116 4.116 26 C 0.051 3.949 27 H 0.067 0.933 28 H 0.067 0.933 29 H 0.070 0.930 30 H 0.059 0.941 31 H 0.062 0.938 32 H 0.074 0.926 33 H 0.403 0.597 34 H 0.078 0.922 35 H 0.054 0.946 36 H 0.063 0.937 37 H 0.057 0.943 38 H 0.019 0.981 39 H 0.254 0.746 40 H 0.019 0.981 41 H 0.088 0.912 42 H 0.051 0.949 43 H 0.054 0.946 44 H 0.037 0.963 45 H 0.082 0.918 46 H 0.073 0.927 47 H 0.066 0.934 48 H 0.063 0.937 49 H 0.056 0.944 50 H 0.072 0.928 51 H 0.064 0.936 52 H 0.080 0.920 53 H 0.037 0.963 For the charges calculated by CM2: Dipole moment (debyes) X Y Z Total from point charges -0.208 11.209 -22.218 24.886 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.195 4.195 2 C -0.003 4.003 3 C -0.225 4.225 4 C 0.309 3.691 5 O -0.419 6.419 6 N -0.362 5.362 7 C 0.046 3.954 8 C -0.177 4.177 9 C 0.022 3.978 10 N -0.315 5.315 11 C -0.385 4.385 12 H 0.083 0.917 13 C -0.196 4.196 14 N -0.546 5.546 15 Si 0.732 3.268 16 H -0.209 1.209 17 H -0.217 1.217 18 H -0.210 1.210 19 C 0.048 3.952 20 C -0.185 4.185 21 N -0.250 5.250 22 C -0.085 4.085 23 C -0.152 4.152 24 C -0.159 4.159 25 C -0.154 4.154 26 C -0.077 4.077 27 H 0.086 0.914 28 H 0.086 0.914 29 H 0.089 0.911 30 H 0.078 0.922 31 H 0.081 0.919 32 H 0.093 0.907 33 H 0.238 0.762 34 H 0.097 0.903 35 H 0.073 0.927 36 H 0.081 0.919 37 H 0.076 0.924 38 H 0.037 0.963 39 H 0.064 0.936 40 H 0.037 0.963 41 H 0.106 0.894 42 H 0.070 0.930 43 H 0.073 0.927 44 H 0.055 0.945 45 H 0.100 0.900 46 H 0.092 0.908 47 H 0.085 0.915 48 H 0.081 0.919 49 H 0.075 0.925 50 H 0.091 0.909 51 H 0.083 0.917 52 H 0.098 0.902 53 H 0.056 0.944 Dipole moment (debyes) X Y Z Total from point charges -0.236 11.171 -21.282 24.037 hybrid contribution -0.243 -0.135 -0.041 0.281 sum -0.479 11.036 -21.322 24.014 Atomic orbital electron populations 1.21847 0.91552 1.03413 1.02639 1.21303 0.97121 0.92939 0.88936 1.22246 1.01797 0.95574 1.02880 1.20155 0.77578 0.82838 0.88540 1.90716 1.49960 1.66380 1.34876 1.46347 1.66048 1.06615 1.17196 1.21033 0.96711 0.92702 0.84986 1.22308 0.99776 0.95662 0.99924 1.20828 0.90691 0.97709 0.88583 1.44074 1.34569 1.05226 1.47622 1.19132 1.37815 0.84966 0.96537 0.91727 1.25057 1.01912 0.97314 0.95299 1.70028 1.47015 1.30470 1.07053 0.86371 0.79884 0.83044 0.77456 1.20898 1.21654 1.21041 1.20741 0.96597 0.86472 0.91377 1.22545 0.97166 0.99849 0.98984 1.62168 1.10654 1.51072 1.01085 1.22280 0.87063 0.99156 0.99975 1.21566 1.00283 0.99820 0.93577 1.21585 0.95928 1.00622 0.97767 1.21583 0.96025 1.01096 0.96660 1.22189 1.00630 0.95626 0.89252 0.91362 0.91439 0.91128 0.92215 0.91944 0.90702 0.76171 0.90345 0.92736 0.91867 0.92434 0.96330 0.93552 0.96262 0.89375 0.92988 0.92710 0.94508 0.89991 0.90836 0.91535 0.91870 0.92481 0.90926 0.91694 0.90180 0.94450 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.14 -1.35 8.23 37.16 0.31 -1.05 16 2 C 0.09 1.00 0.64 -132.75 -0.08 0.92 16 3 C -0.17 -1.79 7.32 37.16 0.27 -1.51 16 4 C 0.52 7.90 5.31 -10.99 -0.06 7.85 16 5 O -0.54 -10.41 15.76 5.56 0.09 -10.33 16 6 N -0.71 -10.38 3.97 -53.77 -0.21 -10.59 16 7 C 0.15 2.64 2.41 -67.89 -0.16 2.48 16 8 C -0.14 -2.31 5.55 -26.61 -0.15 -2.46 16 9 C 0.15 3.10 6.51 -3.71 -0.02 3.08 16 10 N -0.59 -14.95 3.00 -172.30 -0.52 -15.47 16 11 C -0.25 -7.02 9.58 -15.06 -0.14 -7.17 16 12 H 0.07 1.73 7.88 -52.49 -0.41 1.32 16 13 C 0.00 0.08 4.92 -17.43 -0.09 0.00 16 14 N -0.90 -26.58 11.20 55.49 0.62 -25.96 16 15 Si 0.90 21.81 29.59 -169.99 -5.03 16.78 16 16 H -0.28 -6.73 7.11 56.52 0.40 -6.33 16 17 H -0.29 -7.08 7.11 56.52 0.40 -6.68 16 18 H -0.28 -6.78 7.11 56.52 0.40 -6.38 16 19 C 0.17 4.31 5.30 -3.70 -0.02 4.29 16 20 C -0.15 -2.96 5.65 -26.61 -0.15 -3.11 16 21 N -0.52 -4.95 3.76 -220.35 -0.83 -5.78 16 22 C 0.04 0.42 4.61 -3.83 -0.02 0.41 16 23 C -0.11 -0.93 6.09 -26.65 -0.16 -1.10 16 24 C -0.12 -0.84 6.00 -26.61 -0.16 -1.00 16 25 C -0.12 -0.79 6.09 -26.65 -0.16 -0.95 16 26 C 0.05 0.39 4.64 -3.83 -0.02 0.37 16 27 H 0.07 0.72 8.14 -51.93 -0.42 0.30 16 28 H 0.07 0.52 6.28 -51.93 -0.33 0.20 16 29 H 0.07 0.72 8.14 -51.93 -0.42 0.30 16 30 H 0.06 0.69 6.63 -51.93 -0.34 0.34 16 31 H 0.06 0.54 6.54 -51.93 -0.34 0.20 16 32 H 0.07 0.89 7.87 -51.93 -0.41 0.49 16 33 H 0.40 5.17 7.90 -40.82 -0.32 4.85 16 34 H 0.08 1.57 7.58 -51.92 -0.39 1.18 16 35 H 0.05 0.79 8.14 -51.93 -0.42 0.36 16 36 H 0.06 1.07 8.14 -51.93 -0.42 0.64 16 37 H 0.06 1.16 7.93 -51.93 -0.41 0.74 16 38 H 0.02 0.40 8.14 -51.93 -0.42 -0.02 16 39 H 0.25 7.27 8.31 -40.82 -0.34 6.93 16 40 H 0.02 0.51 8.14 -51.93 -0.42 0.09 16 41 H 0.09 2.43 6.07 -51.93 -0.32 2.11 16 42 H 0.05 1.07 8.14 -51.93 -0.42 0.65 16 43 H 0.05 1.11 8.14 -51.92 -0.42 0.69 16 44 H 0.04 0.37 6.71 -51.94 -0.35 0.02 16 45 H 0.08 0.96 6.08 -51.93 -0.32 0.65 16 46 H 0.07 0.63 8.14 -51.93 -0.42 0.21 16 47 H 0.07 0.51 8.14 -51.93 -0.42 0.09 16 48 H 0.06 0.40 8.14 -51.93 -0.42 -0.02 16 49 H 0.06 0.40 8.14 -51.93 -0.42 -0.03 16 50 H 0.07 0.53 8.14 -51.93 -0.42 0.11 16 51 H 0.06 0.37 8.14 -51.93 -0.42 -0.06 16 52 H 0.08 0.57 6.00 -51.93 -0.31 0.26 16 53 H 0.04 0.29 6.63 -51.93 -0.34 -0.05 16 LS Contribution 389.87 15.07 5.88 5.88 Total: -1.00 -30.80 389.87 -10.47 -41.27 By element: Atomic # 1 Polarization: 12.80 SS G_CDS: -9.65 Total: 3.15 kcal Atomic # 6 Polarization: 1.87 SS G_CDS: -0.82 Total: 1.05 kcal Atomic # 7 Polarization: -56.86 SS G_CDS: -0.94 Total: -57.80 kcal Atomic # 8 Polarization: -10.41 SS G_CDS: 0.09 Total: -10.33 kcal Atomic # 14 Polarization: 21.81 SS G_CDS: -5.03 Total: 16.78 kcal Total LS contribution 5.88 Total: 5.88 kcal Total: -30.80 -10.47 -41.27 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. ZINC001002814521.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 19.971 kcal (2) G-P(sol) polarization free energy of solvation -30.796 kcal (3) G-ENP(sol) elect.-nuc.-pol. free energy of system -10.825 kcal (4) G-CDS(sol) cavity-dispersion-solvent structure free energy -10.474 kcal (5) G-P-CDS(sol) = G-P(sol) + G-CDS(sol) = (2) + (4) -41.269 kcal (6) G-S(sol) free energy of system = (1) + (5) -21.298 kcal 1SCF WAS SPECIFIED, SO GEOMETRY OPTIMIZATION WAS NOT USED Total computer time = 0.92 seconds