Wall clock time and date at job start Tue Mar 31 2020 22:51: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.53002 * 1 3 3 C 1.53005 * 109.46982 * 2 1 4 4 O 1.42905 * 109.47079 * 119.99508 * 2 1 3 5 5 C 1.42898 * 113.99615 * 150.00204 * 4 2 1 6 6 H 1.09000 * 112.85228 * 339.66415 * 5 4 2 7 7 C 1.53788 * 113.61515 * 110.86630 * 5 4 2 8 8 C 1.54053 * 86.98447 * 89.15877 * 7 5 4 9 9 H 1.09002 * 113.61560 * 139.98084 * 8 7 5 10 10 N 1.46500 * 113.68900 * 270.87174 * 8 7 5 11 11 C 1.34771 * 120.00282 * 252.53392 * 10 8 7 12 12 O 1.21347 * 119.99820 * 359.97438 * 11 10 8 13 13 C 1.49451 * 120.00302 * 180.02562 * 11 10 8 14 14 O 1.21356 * 119.99696 * 0.02562 * 13 11 10 15 15 N 1.34771 * 120.00302 * 179.97438 * 13 11 10 16 16 C 1.46501 * 119.99911 * 180.02562 * 15 13 11 17 17 C 1.53004 * 109.46897 * 179.97438 * 16 15 13 18 18 C 1.53000 * 109.47179 * 179.97438 * 17 16 15 19 19 C 1.50697 * 109.47401 * 179.97438 * 18 17 16 20 20 H 1.08006 * 120.00048 * 0.02562 * 19 18 17 21 21 C 1.37878 * 120.00221 * 180.02562 * 19 18 17 22 22 N 1.31509 * 119.99804 * 0.02562 * 21 19 18 23 23 Si 1.86803 * 120.00011 * 180.02562 * 21 19 18 24 24 H 1.48499 * 109.47168 * 89.99741 * 23 21 19 25 25 H 1.48503 * 109.47024 * 209.99784 * 23 21 19 26 26 H 1.48493 * 109.46989 * 329.99938 * 23 21 19 27 27 C 1.53782 * 87.07838 * 25.43341 * 8 7 5 28 28 C 1.53004 * 113.61269 * 220.02242 * 27 8 7 29 29 C 1.52994 * 113.61490 * 89.11807 * 27 8 7 30 30 H 1.09004 * 109.47059 * 299.99867 * 1 2 3 31 31 H 1.08992 * 109.47060 * 60.00351 * 1 2 3 32 32 H 1.09004 * 109.46930 * 180.02562 * 1 2 3 33 33 H 1.09000 * 109.47312 * 239.99678 * 2 1 3 34 34 H 1.09001 * 109.46812 * 180.02562 * 3 2 1 35 35 H 1.09002 * 109.46935 * 299.99733 * 3 2 1 36 36 H 1.08992 * 109.47061 * 59.99708 * 3 2 1 37 37 H 1.09001 * 113.61210 * 203.66615 * 7 5 4 38 38 H 1.08999 * 113.60897 * 334.56467 * 7 5 4 39 39 H 0.97008 * 119.99903 * 72.53426 * 10 8 7 40 40 H 0.97000 * 120.00059 * 359.97438 * 15 13 11 41 41 H 1.09001 * 109.46728 * 299.99353 * 16 15 13 42 42 H 1.08995 * 109.47480 * 60.00052 * 16 15 13 43 43 H 1.09002 * 109.47232 * 299.99936 * 17 16 15 44 44 H 1.09009 * 109.47008 * 60.00002 * 17 16 15 45 45 H 1.09002 * 109.47320 * 300.00290 * 18 17 16 46 46 H 1.09002 * 109.47269 * 59.99810 * 18 17 16 47 47 H 0.96999 * 120.00179 * 180.02562 * 22 21 19 48 48 H 1.09003 * 109.47094 * 177.00111 * 28 27 8 49 49 H 1.08994 * 109.46979 * 297.00327 * 28 27 8 50 50 H 1.09002 * 109.46810 * 57.00095 * 28 27 8 51 51 H 1.09000 * 109.47253 * 302.99547 * 29 27 8 52 52 H 1.08999 * 109.47149 * 62.99585 * 29 27 8 53 53 H 1.09004 * 109.47163 * 182.99397 * 29 27 8 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.4426 0.0000 4 8 2.0064 -0.6736 1.1669 5 6 3.2660 -1.3244 0.9887 6 1 3.5006 -1.5131 -0.0589 7 6 4.4159 -0.6648 1.7684 8 6 4.0431 -1.5909 2.9417 9 1 4.9064 -1.9943 3.4709 10 7 3.0375 -1.0231 3.8432 11 6 3.3965 -0.5969 5.0703 12 8 4.5526 -0.6839 5.4287 13 6 2.3706 -0.0182 5.9901 14 8 1.2142 0.0675 5.6322 15 7 2.7294 0.4076 7.2174 16 6 1.7236 0.9743 8.1193 17 6 2.3886 1.3750 9.4377 18 6 1.3384 1.9675 10.3793 19 6 1.9932 2.3616 11.6781 20 1 3.0528 2.2055 11.8178 21 6 1.2399 2.9215 12.6880 22 7 -0.0503 3.1111 12.5181 23 14 2.0518 3.4106 14.2977 24 1 2.5226 4.8161 14.2076 25 1 1.0699 3.2875 15.4049 26 1 3.2085 2.5171 14.5602 27 6 3.4583 -2.5527 1.8939 28 6 2.1445 -3.2061 2.3275 29 6 4.4742 -3.5602 1.3521 30 1 -0.3633 0.5138 0.8900 31 1 -0.3633 0.5137 -0.8899 32 1 -0.3633 -1.0277 0.0005 33 1 1.8934 -0.5139 -0.8899 34 1 3.1300 1.4425 -0.0005 35 1 1.6767 1.9563 -0.8900 36 1 1.6766 1.9563 0.8899 37 1 5.4010 -0.8675 1.3481 38 1 4.2525 0.3926 1.9759 39 1 2.1133 -0.9539 3.5568 40 1 3.6535 0.3380 7.5039 41 1 1.2748 1.8533 7.6565 42 1 0.9502 0.2313 8.3140 43 1 2.8372 0.4958 9.9003 44 1 3.1623 2.1180 9.2433 45 1 0.8901 2.8467 9.9165 46 1 0.5644 1.2250 10.5737 47 1 -0.5803 3.5053 13.2285 48 1 1.7643 -3.8295 1.5181 49 1 2.3189 -3.8222 3.2095 50 1 1.4143 -2.4320 2.5636 51 1 5.3240 -3.0263 0.9268 52 1 4.8175 -4.2020 2.1635 53 1 4.0049 -4.1701 0.5802 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 ZINC000343328236.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 22:51:46 Heat of formation + Delta-G solvation = -118.865489 kcal Electronic energy + Delta-G solvation = -30182.974711 eV Core-core repulsion = 26054.858058 eV Total energy + Delta-G solvation = -4128.116653 eV No. of doubly occupied orbitals = 66 Molecular weight (most abundant/longest-lived isotopes) = 340.220 amu Computer time = 1.04 seconds Orbital eigenvalues (eV) -41.32012 -40.04678 -38.31386 -37.42882 -35.63703 -33.77692 -33.17956 -32.62194 -31.21221 -30.81640 -27.78760 -27.56992 -27.23695 -25.69258 -23.97377 -23.03334 -22.22872 -21.72144 -20.73108 -20.07405 -19.32713 -18.13981 -16.86804 -16.79218 -16.33058 -15.98350 -15.80145 -15.69699 -15.53134 -15.06952 -14.70896 -14.45231 -14.33363 -14.05330 -13.83946 -13.62045 -13.52032 -13.26552 -12.84356 -12.67938 -12.59482 -12.51692 -12.34364 -12.29179 -12.19189 -11.86439 -11.76987 -11.68517 -11.57084 -11.49011 -11.03486 -10.94418 -10.90022 -10.73313 -10.59149 -10.43383 -9.87737 -9.79319 -9.72858 -9.33485 -9.27539 -9.11114 -8.92186 -8.36564 -6.01284 -4.04160 1.49762 2.38546 3.06383 3.29057 3.32128 3.38386 3.41680 3.72518 4.00227 4.10726 4.38339 4.45175 4.51904 4.56282 4.60596 4.68063 4.74719 4.79692 4.86681 5.13588 5.19762 5.24594 5.26018 5.34646 5.36905 5.41668 5.44400 5.47403 5.50910 5.56678 5.63300 5.63698 5.67457 5.72781 5.82695 5.86198 5.89537 6.04607 6.11533 6.23231 6.32748 6.57152 6.73912 6.90839 6.97740 7.08635 7.51788 7.63578 7.75472 7.83711 7.97308 8.12908 8.38673 9.02367 9.60050 11.10019 Molecular weight = 340.22amu Principal moments of inertia in cm(-1) A = 0.018004 B = 0.001887 C = 0.001790 Principal moments of inertia in units of 10**(-40)*gram-cm**2 A = 1554.797873 B =14833.813410 C =15637.079330 Net atomic charges, atomic populations, and dipole contributions calculated with CM2 Atom NO. Type Charge No. of electrons 1 C -0.152 4.152 2 C 0.085 3.915 3 C -0.176 4.176 4 O -0.356 6.356 5 C 0.073 3.927 6 H 0.076 0.924 7 C -0.174 4.174 8 C 0.133 3.867 9 H 0.106 0.894 10 N -0.677 5.677 11 C 0.510 3.490 12 O -0.524 6.524 13 C 0.500 3.500 14 O -0.524 6.524 15 N -0.692 5.692 16 C 0.114 3.886 17 C -0.122 4.122 18 C 0.024 3.976 19 C -0.522 4.522 20 H 0.056 0.944 21 C 0.056 3.944 22 N -0.894 5.894 23 Si 0.894 3.106 24 H -0.280 1.280 25 H -0.286 1.286 26 H -0.274 1.274 27 C -0.093 4.093 28 C -0.135 4.135 29 C -0.125 4.125 30 H 0.067 0.933 31 H 0.068 0.932 32 H 0.061 0.939 33 H 0.059 0.941 34 H 0.056 0.944 35 H 0.063 0.937 36 H 0.066 0.934 37 H 0.100 0.900 38 H 0.081 0.919 39 H 0.419 0.581 40 H 0.411 0.589 41 H 0.070 0.930 42 H 0.070 0.930 43 H 0.057 0.943 44 H 0.057 0.943 45 H 0.020 0.980 46 H 0.020 0.980 47 H 0.261 0.739 48 H 0.057 0.943 49 H 0.057 0.943 50 H 0.069 0.931 51 H 0.061 0.939 52 H 0.067 0.933 53 H 0.064 0.936 For the charges calculated by CM2: Dipole moment (debyes) X Y Z Total from point charges 6.924 -11.915 -31.196 34.104 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.210 4.210 2 C 0.028 3.972 3 C -0.234 4.234 4 O -0.275 6.275 5 C 0.014 3.986 6 H 0.093 0.907 7 C -0.213 4.213 8 C 0.027 3.973 9 H 0.124 0.876 10 N -0.331 5.331 11 C 0.291 3.709 12 O -0.400 6.400 13 C 0.280 3.720 14 O -0.401 6.401 15 N -0.344 5.344 16 C -0.008 4.008 17 C -0.160 4.160 18 C -0.013 4.013 19 C -0.560 4.560 20 H 0.074 0.926 21 C -0.145 4.145 22 N -0.533 5.533 23 Si 0.721 3.279 24 H -0.207 1.207 25 H -0.212 1.212 26 H -0.199 1.199 27 C -0.094 4.094 28 C -0.192 4.192 29 C -0.183 4.183 30 H 0.086 0.914 31 H 0.087 0.913 32 H 0.080 0.920 33 H 0.077 0.923 34 H 0.075 0.925 35 H 0.082 0.918 36 H 0.085 0.915 37 H 0.118 0.882 38 H 0.100 0.900 39 H 0.259 0.741 40 H 0.246 0.754 41 H 0.089 0.911 42 H 0.089 0.911 43 H 0.076 0.924 44 H 0.075 0.925 45 H 0.038 0.962 46 H 0.038 0.962 47 H 0.071 0.929 48 H 0.076 0.924 49 H 0.076 0.924 50 H 0.088 0.912 51 H 0.081 0.919 52 H 0.086 0.914 53 H 0.083 0.917 Dipole moment (debyes) X Y Z Total from point charges 6.625 -11.832 -30.899 33.744 hybrid contribution 0.523 0.467 0.682 0.979 sum 7.148 -11.365 -30.217 33.065 Atomic orbital electron populations 1.21888 0.93422 1.02800 1.02840 1.22190 0.95055 0.92487 0.87496 1.22346 1.01490 0.96780 1.02778 1.87857 1.34894 1.69054 1.35703 1.23263 0.84528 0.92745 0.98029 0.90654 1.24017 0.98949 1.01091 0.97241 1.21970 0.91864 0.93938 0.89553 0.87637 1.45429 1.13554 1.60684 1.13396 1.20356 0.88679 0.78994 0.82902 1.90742 1.19413 1.54389 1.75473 1.20011 0.88382 0.80262 0.83395 1.90764 1.19138 1.54924 1.75238 1.46056 1.14217 1.60236 1.13906 1.21831 0.90591 0.98320 0.90059 1.21926 1.00520 1.00585 0.93017 1.19032 0.93301 0.96042 0.92960 1.21204 0.94950 1.40743 0.99126 0.92574 1.25030 0.94849 0.95836 0.98786 1.69962 1.05577 1.44285 1.33491 0.86570 0.79153 0.78895 0.83240 1.20659 1.21221 1.19869 1.22510 0.97268 0.94189 0.95420 1.21564 0.95790 1.00672 1.01163 1.21550 0.97553 0.97822 1.01349 0.91392 0.91294 0.91954 0.92262 0.92495 0.91779 0.91521 0.88194 0.90023 0.74147 0.75364 0.91146 0.91131 0.92436 0.92459 0.96212 0.96190 0.92933 0.92370 0.92356 0.91188 0.91950 0.91369 0.91682 Number of geometries 1 Number of calculations of the screened coulomb radii 1 The total number of SCF iterations 12. 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.15 -0.98 9.56 37.16 0.36 -0.62 16 2 C 0.09 0.57 3.34 -26.73 -0.09 0.48 16 3 C -0.18 -1.19 9.53 37.16 0.35 -0.84 16 4 O -0.36 -2.97 7.28 -35.23 -0.26 -3.23 16 5 C 0.07 0.50 3.74 -25.91 -0.10 0.40 16 6 H 0.08 0.41 7.42 -51.93 -0.39 0.02 16 7 C -0.17 -1.32 7.38 -25.77 -0.19 -1.51 16 8 C 0.13 1.22 3.78 -66.97 -0.25 0.96 16 9 H 0.11 1.05 7.42 -51.93 -0.39 0.67 16 10 N -0.68 -7.81 4.66 -53.31 -0.25 -8.06 16 11 C 0.51 7.19 7.89 -11.59 -0.09 7.10 16 12 O -0.52 -8.18 16.08 5.50 0.09 -8.09 16 13 C 0.50 7.91 7.89 -11.59 -0.09 7.82 16 14 O -0.52 -9.37 16.48 5.49 0.09 -9.28 16 15 N -0.69 -10.84 5.55 -60.73 -0.34 -11.18 16 16 C 0.11 2.02 5.67 -4.04 -0.02 1.99 16 17 C -0.12 -2.43 5.27 -26.73 -0.14 -2.57 16 18 C 0.02 0.62 4.97 -27.88 -0.14 0.48 16 19 C -0.52 -14.65 9.11 -34.44 -0.31 -14.97 16 20 H 0.06 1.52 7.74 -52.48 -0.41 1.12 16 21 C 0.06 1.61 5.46 -17.82 -0.10 1.52 16 22 N -0.89 -26.86 13.29 55.43 0.74 -26.12 16 23 Si 0.89 21.65 29.54 -169.99 -5.02 16.62 16 24 H -0.28 -6.70 7.11 56.52 0.40 -6.30 16 25 H -0.29 -6.88 7.11 56.52 0.40 -6.48 16 26 H -0.27 -6.49 7.11 56.52 0.40 -6.09 16 27 C -0.09 -0.65 1.34 -153.69 -0.21 -0.85 16 28 C -0.13 -1.03 8.45 37.16 0.31 -0.72 16 29 C -0.13 -0.62 9.07 37.15 0.34 -0.29 16 30 H 0.07 0.52 8.14 -51.93 -0.42 0.10 16 31 H 0.07 0.34 8.14 -51.93 -0.42 -0.08 16 32 H 0.06 0.39 8.14 -51.93 -0.42 -0.03 16 33 H 0.06 0.32 7.42 -51.93 -0.39 -0.06 16 34 H 0.06 0.38 8.14 -51.93 -0.42 -0.04 16 35 H 0.06 0.35 8.14 -51.93 -0.42 -0.07 16 36 H 0.07 0.55 8.14 -51.93 -0.42 0.12 16 37 H 0.10 0.57 7.63 -51.93 -0.40 0.18 16 38 H 0.08 0.75 8.14 -51.93 -0.42 0.33 16 39 H 0.42 4.93 5.51 -40.82 -0.22 4.70 16 40 H 0.41 5.89 7.94 -40.82 -0.32 5.56 16 41 H 0.07 1.28 8.14 -51.93 -0.42 0.86 16 42 H 0.07 1.28 8.14 -51.93 -0.42 0.86 16 43 H 0.06 1.10 8.14 -51.93 -0.42 0.68 16 44 H 0.06 1.10 8.14 -51.92 -0.42 0.68 16 45 H 0.02 0.53 8.14 -51.93 -0.42 0.10 16 46 H 0.02 0.53 8.14 -51.93 -0.42 0.11 16 47 H 0.26 7.40 8.31 -40.82 -0.34 7.06 16 48 H 0.06 0.38 8.14 -51.93 -0.42 -0.04 16 49 H 0.06 0.47 8.14 -51.93 -0.42 0.05 16 50 H 0.07 0.63 5.64 -51.93 -0.29 0.34 16 51 H 0.06 0.25 7.32 -51.93 -0.38 -0.13 16 52 H 0.07 0.35 8.14 -51.93 -0.42 -0.08 16 53 H 0.06 0.27 8.14 -51.93 -0.42 -0.15 16 LS Contribution 427.41 15.07 6.44 6.44 Total: -1.00 -32.12 427.41 -8.38 -40.50 By element: Atomic # 1 Polarization: 13.50 SS G_CDS: -9.50 Total: 4.00 kcal Atomic # 6 Polarization: -1.22 SS G_CDS: -0.37 Total: -1.60 kcal Atomic # 7 Polarization: -45.51 SS G_CDS: 0.15 Total: -45.36 kcal Atomic # 8 Polarization: -20.53 SS G_CDS: -0.08 Total: -20.61 kcal Atomic # 14 Polarization: 21.65 SS G_CDS: -5.02 Total: 16.62 kcal Total LS contribution 6.44 Total: 6.44 kcal Total: -32.12 -8.38 -40.50 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. ZINC000343328236.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 -78.365 kcal (2) G-P(sol) polarization free energy of solvation -32.120 kcal (3) G-ENP(sol) elect.-nuc.-pol. free energy of system -110.485 kcal (4) G-CDS(sol) cavity-dispersion-solvent structure free energy -8.380 kcal (5) G-P-CDS(sol) = G-P(sol) + G-CDS(sol) = (2) + (4) -40.500 kcal (6) G-S(sol) free energy of system = (1) + (5) -118.865 kcal 1SCF WAS SPECIFIED, SO GEOMETRY OPTIMIZATION WAS NOT USED Total computer time = 1.04 seconds