Wall clock time and date at job start Sat Apr 11 2020 06:12:36 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.52993 * 1 3 3 O 1.42893 * 109.47415 * 2 1 4 4 C 1.42901 * 114.00429 * 179.97438 * 3 2 1 5 5 H 1.09006 * 109.47133 * 29.99850 * 4 3 2 6 6 C 1.50696 * 109.47348 * 270.00180 * 4 3 2 7 7 O 1.21280 * 120.00118 * 355.00035 * 6 4 3 8 8 N 1.34773 * 120.00309 * 174.99688 * 6 4 3 9 9 C 1.46503 * 120.00431 * 179.97438 * 8 6 4 10 10 C 1.52998 * 109.47056 * 179.97438 * 9 8 6 11 11 H 1.08996 * 110.88725 * 54.30546 * 10 9 8 12 12 C 1.55162 * 111.00337 * 178.18007 * 10 9 8 13 13 C 1.54319 * 102.94492 * 153.57805 * 12 10 9 14 14 N 1.47022 * 107.26875 * 337.80890 * 13 12 10 15 15 C 1.36943 * 125.64777 * 178.91115 * 14 13 12 16 16 H 1.08001 * 120.00043 * 174.88049 * 15 14 13 17 17 C 1.38813 * 119.99679 * 354.87050 * 15 14 13 18 18 N 1.31618 * 120.00328 * 7.07364 * 17 15 14 19 19 Si 1.86800 * 119.99592 * 187.06507 * 17 15 14 20 20 H 1.48504 * 109.47220 * 94.99426 * 19 17 15 21 21 H 1.48509 * 109.47157 * 214.98962 * 19 17 15 22 22 H 1.48497 * 109.47417 * 334.98973 * 19 17 15 23 23 C 1.47470 * 108.72528 * 358.97450 * 14 13 12 24 24 C 1.53002 * 109.47262 * 150.00038 * 4 3 2 25 25 C 1.53171 * 109.54306 * 59.55634 * 24 4 3 26 26 C 1.53091 * 109.16519 * 183.15144 * 25 24 4 27 27 O 1.42899 * 109.41499 * 57.61011 * 26 25 24 28 28 C 1.42912 * 114.09767 * 298.84658 * 27 26 25 29 29 C 1.53084 * 109.41569 * 61.15340 * 28 27 26 30 30 H 1.09006 * 109.47418 * 179.97438 * 1 2 3 31 31 H 1.08996 * 109.47699 * 299.99923 * 1 2 3 32 32 H 1.09004 * 109.46705 * 60.00312 * 1 2 3 33 33 H 1.09005 * 109.47254 * 239.99720 * 2 1 3 34 34 H 1.09004 * 109.46923 * 120.00193 * 2 1 3 35 35 H 0.97005 * 120.00331 * 0.02562 * 8 6 4 36 36 H 1.09000 * 109.47264 * 300.00517 * 9 8 6 37 37 H 1.09007 * 109.46683 * 60.00110 * 9 8 6 38 38 H 1.08998 * 110.72294 * 271.94426 * 12 10 9 39 39 H 1.09007 * 110.71801 * 35.03357 * 12 10 9 40 40 H 1.08997 * 109.88194 * 97.22421 * 13 12 10 41 41 H 1.09000 * 109.88294 * 218.39285 * 13 12 10 42 42 H 0.96995 * 119.99976 * 180.02562 * 18 17 15 43 43 H 1.08998 * 110.34696 * 265.25414 * 23 14 13 44 44 H 1.08998 * 110.40776 * 142.91642 * 23 14 13 45 45 H 1.09002 * 109.54359 * 179.82358 * 24 4 3 46 46 H 1.08994 * 109.52814 * 303.06222 * 25 24 4 47 47 H 1.09007 * 109.52069 * 63.24749 * 25 24 4 48 48 H 1.08997 * 109.48076 * 297.62837 * 26 25 24 49 49 H 1.08998 * 109.48483 * 177.59316 * 26 25 24 50 50 H 1.09000 * 109.47813 * 301.17123 * 28 27 26 51 51 H 1.09000 * 109.48982 * 181.14717 * 28 27 26 52 52 H 1.09002 * 109.52636 * 182.48527 * 29 28 27 53 53 H 1.09003 * 109.52582 * 62.29860 * 29 28 27 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.5299 0.0000 0.0000 3 8 2.0063 1.3472 0.0000 4 6 3.4309 1.4601 0.0006 5 1 3.8636 0.6014 -0.5129 6 6 3.9311 1.4999 1.4215 7 8 3.1416 1.5293 2.3418 8 7 5.2557 1.5028 1.6699 9 6 5.7420 1.5421 3.0513 10 6 7.2720 1.5363 3.0536 11 1 7.6646 2.3439 2.4357 12 6 7.8262 1.6234 4.5003 13 6 9.2051 0.9404 4.3839 14 7 9.1635 0.0722 3.1981 15 6 10.1987 -0.6945 2.7335 16 1 10.0435 -1.3633 1.8998 17 6 11.4475 -0.6118 3.3340 18 7 11.6037 0.0950 4.4332 19 14 12.9062 -1.5040 2.5821 20 1 13.6586 -0.5712 1.7051 21 1 13.7979 -1.9942 3.6638 22 1 12.4220 -2.6559 1.7797 23 6 7.8229 0.1646 2.5906 24 6 3.8409 2.7456 -0.7207 25 6 3.2246 3.9525 -0.0068 26 6 3.5749 5.2265 -0.7801 27 8 3.1018 5.1120 -2.1236 28 6 3.6794 4.0288 -2.8553 29 6 3.3333 2.7068 -2.1654 30 1 -0.3634 -1.0277 -0.0005 31 1 -0.3634 0.5138 0.8899 32 1 -0.3633 0.5138 -0.8901 33 1 1.8933 -0.5139 -0.8900 34 1 1.8932 -0.5139 0.8900 35 1 5.8872 1.4796 0.9339 36 1 5.3809 2.4481 3.5380 37 1 5.3749 0.6689 3.5906 38 1 7.9356 2.6628 4.8096 39 1 7.1819 1.0821 5.1933 40 1 9.9832 1.6942 4.2643 41 1 9.3990 0.3425 5.2744 42 1 12.4764 0.1531 4.8526 43 1 7.1860 -0.6421 2.9534 44 1 7.8952 0.1323 1.5035 45 1 4.9273 2.8344 -0.7182 46 1 2.1416 3.8376 0.0342 47 1 3.6233 4.0196 1.0055 48 1 4.6564 5.3624 -0.7859 49 1 3.1032 6.0842 -0.3005 50 1 4.7622 4.1491 -2.8883 51 1 3.2831 4.0237 -3.8707 52 1 3.8096 1.8829 -2.6969 53 1 2.2524 2.5666 -2.1680 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 ZINC001714632116.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:12:36 Heat of formation + Delta-G solvation = -119.979703 kcal Electronic energy + Delta-G solvation = -31089.896114 eV Core-core repulsion = 26961.731145 eV Total energy + Delta-G solvation = -4128.164969 eV No. of doubly occupied orbitals = 66 Molecular weight (most abundant/longest-lived isotopes) = 340.220 amu Computer time = 1.53 seconds Orbital eigenvalues (eV) -40.64913 -39.15748 -38.05445 -37.32775 -34.83635 -34.41194 -32.51588 -31.40244 -31.30045 -30.01100 -28.77057 -27.66528 -25.53644 -24.19820 -23.84786 -23.12967 -22.99224 -20.85547 -20.28995 -19.97522 -19.71133 -17.63849 -17.20768 -16.78992 -16.61020 -16.14718 -15.79817 -15.44400 -15.36852 -15.03599 -14.81216 -14.67732 -14.28621 -14.13624 -13.98470 -13.66646 -13.34454 -13.10596 -12.67998 -12.56375 -12.39276 -12.24894 -12.12052 -12.00748 -11.78436 -11.55373 -11.52208 -11.42379 -11.34862 -11.02358 -10.93983 -10.84464 -10.53205 -10.51051 -10.38157 -10.15949 -10.02770 -9.90586 -9.82978 -9.71791 -9.26507 -9.01053 -8.61836 -6.87526 -5.76752 -3.09331 2.53936 2.98086 3.36414 3.41582 3.43915 3.46928 3.62999 4.26335 4.30116 4.45268 4.50817 4.68971 4.71137 4.80349 4.88234 4.97739 5.00537 5.08489 5.12243 5.16558 5.19251 5.21943 5.34837 5.35948 5.38742 5.42050 5.47320 5.49264 5.59292 5.64342 5.70789 5.76457 5.84648 5.88023 5.90452 6.00030 6.01971 6.12900 6.29472 6.38997 6.44378 6.49270 6.63156 6.82210 6.85913 6.93781 7.48916 7.51481 7.59649 7.78422 8.19035 8.21674 9.17423 9.78090 10.44462 11.35388 Molecular weight = 340.22amu Principal moments of inertia in cm(-1) A = 0.012877 B = 0.002463 C = 0.002347 Principal moments of inertia in units of 10**(-40)*gram-cm**2 A = 2173.942573 B =11367.210783 C =11929.112485 Net atomic charges, atomic populations, and dipole contributions calculated with CM2 Atom NO. Type Charge No. of electrons 1 C -0.148 4.148 2 C 0.054 3.946 3 O -0.350 6.350 4 C 0.080 3.920 5 H 0.086 0.914 6 C 0.489 3.511 7 O -0.528 6.528 8 N -0.723 5.723 9 C 0.133 3.867 10 C -0.117 4.117 11 H 0.079 0.921 12 C -0.135 4.135 13 C 0.177 3.823 14 N -0.606 5.606 15 C -0.231 4.231 16 H 0.077 0.923 17 C -0.007 4.007 18 N -0.919 5.919 19 Si 0.910 3.090 20 H -0.289 1.289 21 H -0.291 1.291 22 H -0.281 1.281 23 C 0.143 3.857 24 C -0.080 4.080 25 C -0.148 4.148 26 C 0.061 3.939 27 O -0.383 6.383 28 C 0.061 3.939 29 C -0.148 4.148 30 H 0.070 0.930 31 H 0.065 0.935 32 H 0.061 0.939 33 H 0.048 0.952 34 H 0.051 0.949 35 H 0.400 0.600 36 H 0.065 0.935 37 H 0.067 0.933 38 H 0.061 0.939 39 H 0.059 0.941 40 H 0.048 0.952 41 H 0.041 0.959 42 H 0.253 0.747 43 H 0.018 0.982 44 H 0.035 0.965 45 H 0.088 0.912 46 H 0.085 0.915 47 H 0.073 0.927 48 H 0.050 0.950 49 H 0.094 0.906 50 H 0.051 0.949 51 H 0.094 0.906 52 H 0.073 0.927 53 H 0.083 0.917 For the charges calculated by CM2: Dipole moment (debyes) X Y Z Total from point charges -18.900 6.117 -12.703 23.580 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.205 4.205 2 C -0.022 4.022 3 O -0.268 6.268 4 C 0.018 3.982 5 H 0.104 0.896 6 C 0.274 3.726 7 O -0.403 6.403 8 N -0.377 5.377 9 C 0.009 3.991 10 C -0.137 4.137 11 H 0.098 0.902 12 C -0.174 4.174 13 C 0.052 3.948 14 N -0.333 5.333 15 C -0.360 4.360 16 H 0.095 0.905 17 C -0.203 4.203 18 N -0.566 5.566 19 Si 0.740 3.260 20 H -0.216 1.216 21 H -0.217 1.217 22 H -0.206 1.206 23 C 0.016 3.984 24 C -0.099 4.099 25 C -0.185 4.185 26 C -0.015 4.015 27 O -0.302 6.302 28 C -0.016 4.016 29 C -0.186 4.186 30 H 0.089 0.911 31 H 0.084 0.916 32 H 0.080 0.920 33 H 0.066 0.934 34 H 0.070 0.930 35 H 0.234 0.766 36 H 0.083 0.917 37 H 0.086 0.914 38 H 0.080 0.920 39 H 0.078 0.922 40 H 0.066 0.934 41 H 0.059 0.941 42 H 0.063 0.937 43 H 0.036 0.964 44 H 0.053 0.947 45 H 0.106 0.894 46 H 0.103 0.897 47 H 0.091 0.909 48 H 0.069 0.931 49 H 0.112 0.888 50 H 0.069 0.931 51 H 0.112 0.888 52 H 0.092 0.908 53 H 0.102 0.898 Dipole moment (debyes) X Y Z Total from point charges -19.526 6.364 -12.035 23.803 hybrid contribution 1.640 -1.057 -0.405 1.992 sum -17.886 5.307 -12.440 22.424 Atomic orbital electron populations 1.21697 0.93449 1.02202 1.03158 1.22582 0.95952 0.83581 1.00098 1.87979 1.17868 1.26549 1.94363 1.22013 0.83883 0.99723 0.92533 0.89630 1.21249 0.82636 0.76728 0.92000 1.90709 1.55154 1.48481 1.45996 1.46146 1.06724 1.74989 1.09790 1.21399 0.93896 1.01939 0.81834 1.22611 0.91935 0.99951 0.99210 0.90229 1.22852 0.99689 1.00537 0.94272 1.21153 0.90547 0.92956 0.90159 1.44496 1.08105 1.51868 1.28801 1.19049 0.87627 1.16664 1.12612 0.90488 1.24827 0.97621 1.00249 0.97620 1.69976 1.22596 1.44839 1.19232 0.86128 0.82767 0.78858 0.78199 1.21598 1.21735 1.20639 1.21572 0.87943 0.92831 0.96010 1.21305 1.00344 0.93029 0.95186 1.22177 1.03200 0.93941 0.99213 1.22596 0.98738 0.97749 0.82441 1.87935 1.67560 1.50143 1.24554 1.22568 0.97662 0.85871 0.95479 1.22208 1.03422 0.98451 0.94511 0.91142 0.91606 0.91967 0.93405 0.93028 0.76558 0.91696 0.91442 0.91981 0.92232 0.93393 0.94075 0.93710 0.96441 0.94673 0.89413 0.89677 0.90859 0.93127 0.88778 0.93117 0.88821 0.90798 0.89798 Number of geometries 1 Number of calculations of the screened coulomb radii 1 The total number of SCF iterations 15. 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 -1.21 10.28 37.15 0.38 -0.83 16 2 C 0.05 0.52 6.46 37.15 0.24 0.76 16 3 O -0.35 -3.98 8.07 -55.09 -0.44 -4.42 16 4 C 0.08 0.81 2.26 -27.88 -0.06 0.75 16 5 H 0.09 0.70 7.79 -51.93 -0.40 0.30 16 6 C 0.49 6.34 5.96 -10.99 -0.07 6.27 16 7 O -0.53 -8.96 16.51 -14.30 -0.24 -9.20 16 8 N -0.72 -8.66 5.41 -60.29 -0.33 -8.99 16 9 C 0.13 1.85 5.12 -4.04 -0.02 1.83 16 10 C -0.12 -1.87 3.07 -88.47 -0.27 -2.14 16 11 H 0.08 1.24 8.14 -51.93 -0.42 0.81 16 12 C -0.13 -2.40 6.42 -24.90 -0.16 -2.56 16 13 C 0.18 4.19 5.38 -3.07 -0.02 4.17 16 14 N -0.61 -15.10 3.40 -170.04 -0.58 -15.67 16 15 C -0.23 -6.34 10.26 -15.06 -0.15 -6.49 16 16 H 0.08 2.03 7.84 -52.49 -0.41 1.62 16 17 C -0.01 -0.19 5.48 -17.43 -0.10 -0.28 16 18 N -0.92 -26.74 10.36 55.49 0.57 -26.16 16 19 Si 0.91 21.75 29.56 -169.99 -5.02 16.73 16 20 H -0.29 -6.88 7.11 56.52 0.40 -6.48 16 21 H -0.29 -7.04 7.11 56.52 0.40 -6.64 16 22 H -0.28 -6.52 7.11 56.52 0.40 -6.12 16 23 C 0.14 2.77 6.56 -2.51 -0.02 2.75 16 24 C -0.08 -0.69 2.13 -90.44 -0.19 -0.88 16 25 C -0.15 -1.46 5.01 -26.60 -0.13 -1.59 16 26 C 0.06 0.55 6.84 37.20 0.25 0.80 16 27 O -0.38 -3.89 10.77 -35.23 -0.38 -4.27 16 28 C 0.06 0.47 6.84 37.20 0.25 0.72 16 29 C -0.15 -1.11 5.32 -26.60 -0.14 -1.26 16 30 H 0.07 0.45 8.14 -51.93 -0.42 0.03 16 31 H 0.06 0.61 8.14 -51.93 -0.42 0.19 16 32 H 0.06 0.47 8.14 -51.93 -0.42 0.05 16 33 H 0.05 0.37 8.09 -51.93 -0.42 -0.05 16 34 H 0.05 0.56 8.14 -51.93 -0.42 0.14 16 35 H 0.40 4.01 8.22 -40.82 -0.34 3.67 16 36 H 0.06 0.90 8.14 -51.93 -0.42 0.48 16 37 H 0.07 1.03 8.03 -51.93 -0.42 0.61 16 38 H 0.06 0.99 8.14 -51.93 -0.42 0.57 16 39 H 0.06 1.01 8.04 -51.93 -0.42 0.59 16 40 H 0.05 1.26 7.30 -51.93 -0.38 0.88 16 41 H 0.04 1.09 7.61 -51.93 -0.40 0.69 16 42 H 0.25 7.13 8.31 -40.82 -0.34 6.79 16 43 H 0.02 0.35 7.96 -51.93 -0.41 -0.06 16 44 H 0.04 0.64 8.14 -51.93 -0.42 0.22 16 45 H 0.09 0.66 7.92 -51.93 -0.41 0.25 16 46 H 0.08 0.98 7.59 -51.93 -0.39 0.59 16 47 H 0.07 0.80 7.36 -51.93 -0.38 0.42 16 48 H 0.05 0.41 8.14 -51.93 -0.42 -0.01 16 49 H 0.09 0.76 8.14 -51.93 -0.42 0.34 16 50 H 0.05 0.35 8.14 -51.93 -0.42 -0.08 16 51 H 0.09 0.59 8.14 -51.93 -0.42 0.17 16 52 H 0.07 0.44 8.14 -51.93 -0.42 0.02 16 53 H 0.08 0.73 7.60 -51.93 -0.39 0.34 16 LS Contribution 414.28 15.07 6.24 6.24 Total: -1.00 -33.18 414.28 -10.18 -43.36 By element: Atomic # 1 Polarization: 10.16 SS G_CDS: -9.81 Total: 0.35 kcal Atomic # 6 Polarization: 2.24 SS G_CDS: -0.20 Total: 2.04 kcal Atomic # 7 Polarization: -50.49 SS G_CDS: -0.33 Total: -50.82 kcal Atomic # 8 Polarization: -16.83 SS G_CDS: -1.06 Total: -17.89 kcal Atomic # 14 Polarization: 21.75 SS G_CDS: -5.02 Total: 16.73 kcal Total LS contribution 6.24 Total: 6.24 kcal Total: -33.18 -10.18 -43.36 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. ZINC001714632116.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 -76.624 kcal (2) G-P(sol) polarization free energy of solvation -33.180 kcal (3) G-ENP(sol) elect.-nuc.-pol. free energy of system -109.804 kcal (4) G-CDS(sol) cavity-dispersion-solvent structure free energy -10.176 kcal (5) G-P-CDS(sol) = G-P(sol) + G-CDS(sol) = (2) + (4) -43.356 kcal (6) G-S(sol) free energy of system = (1) + (5) -119.980 kcal 1SCF WAS SPECIFIED, SO GEOMETRY OPTIMIZATION WAS NOT USED Total computer time = 1.53 seconds