Wall clock time and date at job start Tue Mar 31 2020 06:32:56 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.52996 * 1 3 3 C 1.53006 * 109.46799 * 2 1 4 4 C 1.53006 * 109.47139 * 119.99503 * 2 1 3 5 5 H 1.09003 * 109.46527 * 60.00441 * 4 2 1 6 6 N 1.46491 * 109.47294 * 300.00379 * 4 2 1 7 7 C 1.34776 * 120.00311 * 84.99666 * 6 4 2 8 8 O 1.21287 * 119.99661 * 0.02562 * 7 6 4 9 9 C 1.50699 * 120.00433 * 179.97438 * 7 6 4 10 10 S 1.81002 * 109.47437 * 180.02562 * 9 7 6 11 11 C 1.76199 * 99.99873 * 180.02562 * 10 9 7 12 12 H 1.07991 * 120.00360 * 359.97438 * 11 10 9 13 13 C 1.38801 * 119.99612 * 180.02562 * 11 10 9 14 14 N 1.31633 * 119.99815 * 179.97438 * 13 11 10 15 15 Si 1.86799 * 120.00325 * 359.97438 * 13 11 10 16 16 H 1.48501 * 109.47294 * 89.99662 * 15 13 11 17 17 H 1.48500 * 109.47152 * 210.00197 * 15 13 11 18 18 H 1.48501 * 109.46788 * 329.99926 * 15 13 11 19 19 C 1.50696 * 109.47078 * 179.97438 * 4 2 1 20 20 C 1.38230 * 120.00328 * 260.00371 * 19 4 2 21 21 C 1.38235 * 120.00326 * 179.97438 * 20 19 4 22 22 C 1.38234 * 119.99814 * 359.97438 * 21 20 19 23 23 C 1.50706 * 120.00196 * 179.97438 * 22 21 20 24 24 C 1.53000 * 109.46964 * 60.00061 * 23 22 21 25 25 C 1.52999 * 109.46977 * 179.97438 * 23 22 21 26 26 C 1.52998 * 109.46922 * 300.00453 * 23 22 21 27 27 C 1.38240 * 120.00061 * 0.07178 * 22 21 20 28 28 C 1.38234 * 119.99689 * 359.90912 * 27 22 21 29 29 H 1.09006 * 109.46995 * 180.02562 * 1 2 3 30 30 H 1.09005 * 109.47418 * 299.99907 * 1 2 3 31 31 H 1.08993 * 109.47517 * 60.00411 * 1 2 3 32 32 H 1.08996 * 109.47031 * 239.99667 * 2 1 3 33 33 H 1.09001 * 109.46843 * 179.97438 * 3 2 1 34 34 H 1.09002 * 109.46754 * 299.99847 * 3 2 1 35 35 H 1.08995 * 109.46952 * 60.00395 * 3 2 1 36 36 H 0.97004 * 120.00167 * 265.00172 * 6 4 2 37 37 H 1.08997 * 109.47019 * 300.00321 * 9 7 6 38 38 H 1.08998 * 109.46703 * 59.99953 * 9 7 6 39 39 H 0.97003 * 119.99878 * 179.97438 * 14 13 11 40 40 H 1.08000 * 120.00376 * 0.02562 * 20 19 4 41 41 H 1.08000 * 120.00616 * 179.72578 * 21 20 19 42 42 H 1.09009 * 109.46612 * 59.99991 * 24 23 22 43 43 H 1.08991 * 109.47096 * 179.97438 * 24 23 22 44 44 H 1.08992 * 109.47275 * 300.00500 * 24 23 22 45 45 H 1.08995 * 109.47063 * 60.00845 * 25 23 22 46 46 H 1.09001 * 109.46742 * 179.97438 * 25 23 22 47 47 H 1.08993 * 109.47496 * 300.00240 * 25 23 22 48 48 H 1.09001 * 109.47153 * 59.99823 * 26 23 22 49 49 H 1.09001 * 109.46779 * 179.97438 * 26 23 22 50 50 H 1.09003 * 109.46942 * 299.99559 * 26 23 22 51 51 H 1.07997 * 119.99852 * 179.97438 * 27 22 21 52 52 H 1.07999 * 119.99751 * 180.02562 * 28 27 22 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.0399 1.4426 0.0000 4 6 2.0400 -0.7212 1.2493 5 1 1.6767 -1.7489 1.2493 6 7 1.5516 -0.0306 2.4454 7 6 0.3297 -0.3164 2.9369 8 8 -0.3655 -1.1457 2.3892 9 6 -0.1725 0.3936 4.1677 10 16 -1.8316 -0.2047 4.5748 11 6 -2.1632 0.7647 6.0083 12 1 -1.4253 1.4661 6.3684 13 6 -3.3725 0.6260 6.6752 14 7 -3.6201 1.3498 7.7465 15 14 -4.6487 -0.5878 6.0528 16 1 -4.4282 -1.9096 6.6927 17 1 -6.0078 -0.0944 6.3914 18 1 -4.5246 -0.7243 4.5793 19 6 3.5469 -0.7217 1.2490 20 6 4.2377 -1.8469 0.8396 21 6 5.6200 -1.8472 0.8388 22 6 6.3116 -0.7221 1.2469 23 6 7.8187 -0.7225 1.2466 24 6 8.3283 -0.9746 -0.1738 25 6 8.3292 0.6336 1.7378 26 6 8.3285 -1.8268 2.1748 27 6 5.6209 0.4029 1.6573 28 6 4.2386 0.4033 1.6577 29 1 -0.3633 -1.0277 0.0005 30 1 -0.3634 0.5138 0.8900 31 1 -0.3634 0.5137 -0.8899 32 1 1.8933 -0.5139 -0.8899 33 1 3.1299 1.4426 0.0005 34 1 1.6765 1.9563 -0.8900 35 1 1.6766 1.9563 0.8900 36 1 2.1077 0.6325 2.8837 37 1 0.5001 0.1944 5.0019 38 1 -0.2090 1.4665 3.9787 39 1 -4.4654 1.2531 8.2124 40 1 3.6974 -2.7256 0.5196 41 1 6.1597 -2.7278 0.5232 42 1 7.9650 -0.1879 -0.8351 43 1 9.4182 -0.9745 -0.1741 44 1 7.9646 -1.9405 -0.5242 45 1 7.9662 0.8136 2.7497 46 1 9.4192 0.6331 1.7381 47 1 7.9659 1.4199 1.0763 48 1 7.9648 -2.7927 1.8244 49 1 9.4185 -1.8274 2.1743 50 1 7.9654 -1.6471 3.1868 51 1 6.1613 1.2818 1.9765 52 1 3.6989 1.2823 1.9778 RHF calculation, no. of doubly occupied orbitals= 64 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). S: (AM1): M.J.S. DEWAR, Y.-C. YUAN, INORGANIC CHEM., 29, 589 (1990) 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 ZINC000040088104.mol2 52 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 06:32:56 Heat of formation + Delta-G solvation = -37.209986 kcal Electronic energy + Delta-G solvation = -28615.008718 eV Core-core repulsion = 24897.697033 eV Total energy + Delta-G solvation = -3717.311685 eV No. of doubly occupied orbitals = 64 Molecular weight (most abundant/longest-lived isotopes) = 349.182 amu Computer time = 1.07 seconds Orbital eigenvalues (eV) -40.72170 -39.48542 -37.39903 -35.54546 -33.04256 -32.42920 -31.44404 -30.71926 -29.13846 -27.82632 -27.62652 -27.00096 -25.93747 -25.54210 -22.73043 -22.31545 -20.73287 -19.65763 -19.41276 -18.00978 -17.11015 -16.51611 -15.96728 -15.69230 -15.26148 -15.18690 -14.96222 -14.65380 -14.28293 -14.13431 -14.04973 -13.88074 -13.71220 -13.49024 -13.24127 -12.95780 -12.78910 -12.77676 -12.37947 -12.34755 -12.23277 -12.06717 -11.93028 -11.77447 -11.49614 -11.41253 -11.25879 -11.24274 -11.04584 -11.00338 -10.81054 -10.78232 -10.67019 -10.51644 -9.85267 -9.47035 -9.10877 -9.00209 -8.92194 -8.68360 -8.48895 -6.90279 -6.28232 -3.99540 1.07325 1.18968 2.78741 3.00859 3.30983 3.34537 3.50673 3.52046 3.71562 4.19938 4.33677 4.35711 4.47100 4.61326 4.74177 4.75248 4.76377 4.88015 4.91419 4.99896 5.10865 5.12945 5.16337 5.21670 5.23467 5.30812 5.36469 5.37344 5.49081 5.49762 5.57856 5.61356 5.66895 5.67423 5.74164 5.76602 5.78863 5.83137 5.89298 5.91689 5.92954 6.00403 6.03147 6.11809 6.22985 6.27684 6.33415 6.39320 6.58609 6.66904 6.76403 7.64701 7.94890 8.42579 8.53593 9.18010 10.81025 Molecular weight = 349.18amu Principal moments of inertia in cm(-1) A = 0.017416 B = 0.002211 C = 0.002094 Principal moments of inertia in units of 10**(-40)*gram-cm**2 A = 1607.333961 B =12662.783209 C =13369.676803 Net atomic charges, atomic populations, and dipole contributions calculated with CM2 Atom NO. Type Charge No. of electrons 1 C -0.141 4.141 2 C -0.099 4.099 3 C -0.152 4.152 4 C 0.192 3.808 5 H 0.095 0.905 6 N -0.728 5.728 7 C 0.532 3.468 8 O -0.514 6.514 9 C -0.125 4.125 10 S -0.139 6.139 11 C -0.521 4.521 12 H 0.088 0.912 13 C 0.095 3.905 14 N -0.881 5.881 15 Si 0.857 3.143 16 H -0.267 1.267 17 H -0.275 1.275 18 H -0.235 1.235 19 C -0.088 4.088 20 C -0.105 4.105 21 C -0.113 4.113 22 C -0.071 4.071 23 C -0.030 4.030 24 C -0.138 4.138 25 C -0.145 4.145 26 C -0.139 4.139 27 C -0.120 4.120 28 C -0.108 4.108 29 H 0.063 0.937 30 H 0.060 0.940 31 H 0.045 0.955 32 H 0.074 0.926 33 H 0.059 0.941 34 H 0.052 0.948 35 H 0.056 0.944 36 H 0.397 0.603 37 H 0.095 0.905 38 H 0.095 0.905 39 H 0.270 0.730 40 H 0.123 0.877 41 H 0.120 0.880 42 H 0.056 0.944 43 H 0.057 0.943 44 H 0.057 0.943 45 H 0.059 0.941 46 H 0.057 0.943 47 H 0.059 0.941 48 H 0.057 0.943 49 H 0.057 0.943 50 H 0.057 0.943 51 H 0.126 0.874 52 H 0.123 0.877 For the charges calculated by CM2: Dipole moment (debyes) X Y Z Total from point charges 25.484 -2.477 -17.174 30.830 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.198 4.198 2 C -0.118 4.118 3 C -0.209 4.209 4 C 0.087 3.913 5 H 0.113 0.887 6 N -0.382 5.382 7 C 0.319 3.681 8 O -0.388 6.388 9 C -0.276 4.276 10 S 0.087 5.913 11 C -0.669 4.669 12 H 0.106 0.894 13 C -0.114 4.114 14 N -0.517 5.517 15 Si 0.682 3.318 16 H -0.192 1.192 17 H -0.200 1.200 18 H -0.159 1.159 19 C -0.089 4.089 20 C -0.123 4.123 21 C -0.131 4.131 22 C -0.072 4.072 23 C -0.030 4.030 24 C -0.195 4.195 25 C -0.202 4.202 26 C -0.196 4.196 27 C -0.137 4.137 28 C -0.127 4.127 29 H 0.082 0.918 30 H 0.079 0.921 31 H 0.064 0.936 32 H 0.092 0.908 33 H 0.078 0.922 34 H 0.071 0.929 35 H 0.075 0.925 36 H 0.232 0.768 37 H 0.114 0.886 38 H 0.114 0.886 39 H 0.080 0.920 40 H 0.141 0.859 41 H 0.138 0.862 42 H 0.075 0.925 43 H 0.076 0.924 44 H 0.076 0.924 45 H 0.078 0.922 46 H 0.076 0.924 47 H 0.078 0.922 48 H 0.076 0.924 49 H 0.076 0.924 50 H 0.076 0.924 51 H 0.144 0.856 52 H 0.141 0.859 Dipole moment (debyes) X Y Z Total from point charges 24.180 -3.813 -17.931 30.344 hybrid contribution 0.393 0.928 0.907 1.356 sum 24.573 -2.885 -17.024 30.033 Atomic orbital electron populations 1.21893 0.93633 1.02831 1.01465 1.21470 0.96990 0.96517 0.96847 1.21991 1.01041 0.96010 1.01885 1.19813 0.88536 0.97836 0.85141 0.88706 1.46069 1.18211 1.44303 1.29637 1.20180 0.80681 0.82489 0.84797 1.90703 1.55535 1.32872 1.59684 1.23288 1.03222 1.03373 0.97672 1.84759 1.12652 1.56659 1.37253 1.22324 1.06541 1.26583 1.11423 0.89434 1.24720 0.97421 0.95348 0.93878 1.69942 1.24762 1.39290 1.17732 0.86731 0.81116 0.81673 0.82299 1.19237 1.19982 1.15914 1.19906 0.94743 0.94308 0.99919 1.20975 0.94097 0.97590 0.99684 1.20914 0.93804 0.97944 1.00461 1.19804 0.94668 0.93931 0.98753 1.19507 0.90900 0.96414 0.96199 1.21887 1.00729 1.01902 0.95004 1.21966 1.00927 0.95942 1.01339 1.21894 1.00724 0.97707 0.99289 1.21146 0.93861 0.98316 1.00413 1.20994 0.94144 0.97760 0.99755 0.91775 0.92052 0.93628 0.90800 0.92186 0.92873 0.92494 0.76804 0.88647 0.88634 0.91968 0.85909 0.86197 0.92491 0.92392 0.92432 0.92153 0.92402 0.92235 0.92420 0.92381 0.92372 0.85606 0.85939 Number of geometries 1 Number of calculations of the screened coulomb radii 1 The total number of SCF iterations 14. 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.91 7.72 37.16 0.29 -1.63 16 2 C -0.10 -1.09 2.79 -90.62 -0.25 -1.34 16 3 C -0.15 -1.49 7.96 37.16 0.30 -1.19 16 4 C 0.19 2.22 2.28 -69.09 -0.16 2.06 16 5 H 0.09 1.26 7.44 -51.93 -0.39 0.88 16 6 N -0.73 -9.28 4.15 -54.94 -0.23 -9.50 16 7 C 0.53 9.15 6.20 -10.99 -0.07 9.08 16 8 O -0.51 -10.42 12.99 5.55 0.07 -10.34 16 9 C -0.12 -2.40 6.16 36.01 0.22 -2.17 16 10 S -0.14 -3.51 18.55 -107.50 -1.99 -5.50 16 11 C -0.52 -15.01 10.04 30.94 0.31 -14.70 16 12 H 0.09 2.63 8.03 -52.49 -0.42 2.21 16 13 C 0.09 2.78 5.50 -17.43 -0.10 2.69 16 14 N -0.88 -26.63 13.97 55.49 0.78 -25.86 16 15 Si 0.86 21.84 27.74 -169.99 -4.71 17.12 16 16 H -0.27 -6.74 7.11 56.52 0.40 -6.34 16 17 H -0.27 -6.74 7.11 56.52 0.40 -6.34 16 18 H -0.24 -5.97 6.74 56.52 0.38 -5.59 16 19 C -0.09 -0.87 4.71 -104.63 -0.49 -1.37 16 20 C -0.10 -0.98 9.67 -39.59 -0.38 -1.36 16 21 C -0.11 -0.95 8.71 -39.58 -0.34 -1.30 16 22 C -0.07 -0.58 4.76 -104.62 -0.50 -1.08 16 23 C -0.03 -0.19 0.85 -155.66 -0.13 -0.33 16 24 C -0.14 -0.80 8.29 37.16 0.31 -0.49 16 25 C -0.14 -0.83 7.13 37.16 0.27 -0.57 16 26 C -0.14 -0.86 8.29 37.16 0.31 -0.55 16 27 C -0.12 -1.00 8.98 -39.58 -0.36 -1.35 16 28 C -0.11 -0.97 6.84 -39.58 -0.27 -1.24 16 29 H 0.06 0.96 7.53 -51.93 -0.39 0.57 16 30 H 0.06 0.97 5.78 -51.93 -0.30 0.67 16 31 H 0.04 0.56 8.14 -51.93 -0.42 0.13 16 32 H 0.07 0.75 8.14 -51.93 -0.42 0.33 16 33 H 0.06 0.54 5.60 -51.93 -0.29 0.24 16 34 H 0.05 0.50 8.14 -51.93 -0.42 0.08 16 35 H 0.06 0.58 7.60 -51.93 -0.39 0.19 16 36 H 0.40 4.20 6.87 -40.82 -0.28 3.92 16 37 H 0.10 1.73 8.14 -51.92 -0.42 1.31 16 38 H 0.10 1.71 8.14 -51.92 -0.42 1.29 16 39 H 0.27 7.60 8.31 -40.82 -0.34 7.26 16 40 H 0.12 1.03 8.04 -52.49 -0.42 0.61 16 41 H 0.12 0.83 7.10 -52.49 -0.37 0.46 16 42 H 0.06 0.35 8.14 -51.92 -0.42 -0.07 16 43 H 0.06 0.28 8.14 -51.93 -0.42 -0.14 16 44 H 0.06 0.32 7.39 -51.93 -0.38 -0.06 16 45 H 0.06 0.36 7.51 -51.93 -0.39 -0.03 16 46 H 0.06 0.29 8.14 -51.93 -0.42 -0.13 16 47 H 0.06 0.33 7.52 -51.93 -0.39 -0.06 16 48 H 0.06 0.35 7.38 -51.93 -0.38 -0.03 16 49 H 0.06 0.30 8.14 -51.93 -0.42 -0.12 16 50 H 0.06 0.40 8.14 -51.93 -0.42 -0.02 16 51 H 0.13 0.87 5.90 -52.49 -0.31 0.56 16 52 H 0.12 0.98 6.09 -52.49 -0.32 0.66 16 LS Contribution 410.71 15.07 6.19 6.19 Total: -1.00 -32.53 410.71 -9.77 -42.30 By element: Atomic # 1 Polarization: 11.25 SS G_CDS: -8.82 Total: 2.44 kcal Atomic # 6 Polarization: -15.78 SS G_CDS: -1.06 Total: -16.84 kcal Atomic # 7 Polarization: -35.91 SS G_CDS: 0.55 Total: -35.36 kcal Atomic # 8 Polarization: -10.42 SS G_CDS: 0.07 Total: -10.34 kcal Atomic # 14 Polarization: 21.84 SS G_CDS: -4.71 Total: 17.12 kcal Atomic # 16 Polarization: -3.51 SS G_CDS: -1.99 Total: -5.50 kcal Total LS contribution 6.19 Total: 6.19 kcal Total: -32.53 -9.77 -42.30 kcal The number of atoms in the molecule is 52 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. ZINC000040088104.mol2 52 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 5.092 kcal (2) G-P(sol) polarization free energy of solvation -32.528 kcal (3) G-ENP(sol) elect.-nuc.-pol. free energy of system -27.435 kcal (4) G-CDS(sol) cavity-dispersion-solvent structure free energy -9.774 kcal (5) G-P-CDS(sol) = G-P(sol) + G-CDS(sol) = (2) + (4) -42.302 kcal (6) G-S(sol) free energy of system = (1) + (5) -37.210 kcal 1SCF WAS SPECIFIED, SO GEOMETRY OPTIMIZATION WAS NOT USED Total computer time = 1.07 seconds