Wall clock time and date at job start Sat Apr 11 2020 06:12:20 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.42901 * 1 3 3 C 1.42898 * 114.00025 * 2 1 4 4 C 1.52996 * 109.47136 * 179.97438 * 3 2 1 5 5 C 1.50703 * 109.54576 * 268.09815 * 4 3 2 6 6 O 1.21280 * 120.00066 * 235.00460 * 5 4 3 7 7 N 1.34777 * 119.99793 * 55.00652 * 5 4 3 8 8 C 1.47419 * 125.64786 * 5.79580 * 7 5 4 9 9 C 1.55135 * 104.72239 * 155.81727 * 8 7 5 10 10 C 1.55161 * 101.48237 * 37.29843 * 9 8 7 11 11 C 1.47023 * 125.64557 * 186.10502 * 7 5 4 12 12 H 1.08994 * 109.88545 * 300.43960 * 11 7 5 13 13 C 1.52995 * 109.88413 * 61.45112 * 11 7 5 14 14 N 1.46501 * 109.47427 * 303.38351 * 13 11 7 15 15 C 1.46503 * 119.99556 * 278.70243 * 14 13 11 16 16 C 1.36935 * 120.00672 * 98.69853 * 14 13 11 17 17 H 1.08007 * 119.99773 * 318.82063 * 16 14 13 18 18 C 1.38807 * 120.00361 * 138.81562 * 16 14 13 19 19 N 1.31618 * 119.99773 * 357.53951 * 18 16 14 20 20 Si 1.86799 * 120.00243 * 177.53595 * 18 16 14 21 21 H 1.48498 * 109.47097 * 90.00359 * 20 18 16 22 22 H 1.48497 * 109.47428 * 210.00447 * 20 18 16 23 23 H 1.48507 * 109.47105 * 330.00209 * 20 18 16 24 24 C 1.53176 * 109.54559 * 147.91998 * 4 3 2 25 25 C 1.53093 * 109.16810 * 183.15295 * 24 4 3 26 26 O 1.42897 * 109.41270 * 57.60922 * 25 24 4 27 27 C 1.42906 * 114.10127 * 298.84293 * 26 25 24 28 28 C 1.53099 * 109.41254 * 61.15498 * 27 26 25 29 29 H 1.09003 * 109.46901 * 59.99439 * 1 2 3 30 30 H 1.09000 * 109.46967 * 179.97438 * 1 2 3 31 31 H 1.08998 * 109.46894 * 299.99966 * 1 2 3 32 32 H 1.09005 * 109.46912 * 60.00266 * 3 2 1 33 33 H 1.08997 * 109.47593 * 300.00493 * 3 2 1 34 34 H 1.09000 * 110.37005 * 274.74744 * 8 7 5 35 35 H 1.08997 * 110.36787 * 37.07368 * 8 7 5 36 36 H 1.08999 * 111.19535 * 279.19334 * 9 8 7 37 37 H 1.09003 * 110.93781 * 155.28474 * 9 8 7 38 38 H 1.09009 * 110.71633 * 205.97200 * 10 9 8 39 39 H 1.08995 * 110.87626 * 82.76577 * 10 9 8 40 40 H 1.09004 * 109.47182 * 63.37809 * 13 11 7 41 41 H 1.09001 * 109.47351 * 183.37789 * 13 11 7 42 42 H 1.09004 * 109.47268 * 84.26064 * 15 14 13 43 43 H 1.08995 * 109.47002 * 204.25808 * 15 14 13 44 44 H 1.08996 * 109.47134 * 324.25665 * 15 14 13 45 45 H 0.97005 * 120.00033 * 179.97438 * 19 18 16 46 46 H 1.09003 * 109.52462 * 63.24766 * 24 4 3 47 47 H 1.09003 * 109.57311 * 303.09361 * 24 4 3 48 48 H 1.08996 * 109.48073 * 177.58997 * 25 24 4 49 49 H 1.09003 * 109.47892 * 297.62766 * 25 24 4 50 50 H 1.09000 * 109.47948 * 301.18099 * 27 26 25 51 51 H 1.08995 * 109.49519 * 181.14873 * 27 26 25 52 52 H 1.09003 * 109.51902 * 182.48462 * 28 27 26 53 53 H 1.08998 * 109.51923 * 62.30150 * 28 27 26 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.0102 1.3054 0.0000 4 6 3.5354 1.1846 0.0006 5 6 4.0412 1.1921 1.4203 6 8 4.8780 2.0008 1.7616 7 7 3.5635 0.2999 2.3104 8 6 2.6408 -0.8086 2.0052 9 6 2.8908 -1.8500 3.1276 10 6 3.1239 -0.9307 4.3556 11 6 3.8713 0.2750 3.7479 12 1 4.9449 0.1582 3.8952 13 6 3.3943 1.5719 4.4047 14 7 4.0124 2.7166 3.7311 15 6 5.3783 3.1191 4.0756 16 6 3.3137 3.4104 2.7795 17 1 2.7081 2.8793 2.0599 18 6 3.3834 4.7961 2.7398 19 7 4.0823 5.4461 3.6461 20 14 2.4857 5.7385 1.3999 21 1 1.1135 6.0689 1.8615 22 1 3.2218 6.9925 1.0984 23 1 2.4080 4.9033 0.1743 24 6 4.1432 2.3633 -0.7660 25 6 5.6633 2.1925 -0.8277 26 8 5.9761 0.9440 -1.4484 27 6 5.4740 -0.1970 -0.7497 28 6 3.9462 -0.1222 -0.6847 29 1 -0.3633 0.5139 -0.8900 30 1 -0.3633 -1.0277 -0.0005 31 1 -0.3633 0.5138 0.8900 32 1 1.6886 1.8463 0.8901 33 1 1.6886 1.8464 -0.8899 34 1 1.6082 -0.4606 2.0306 35 1 2.8734 -1.2384 1.0310 36 1 3.7747 -2.4520 2.9168 37 1 2.0154 -2.4827 3.2739 38 1 3.7419 -1.4318 5.1007 39 1 2.1754 -0.6167 4.7913 40 1 2.3098 1.6437 4.3217 41 1 3.6789 1.5723 5.4569 42 1 6.0884 2.5192 3.5062 43 1 5.5181 4.1729 3.8351 44 1 5.5447 2.9645 5.1416 45 1 4.1314 6.4145 3.6181 46 1 3.9018 3.2945 -0.2534 47 1 3.7380 2.3888 -1.7776 48 1 6.0971 3.0062 -1.4090 49 1 6.0722 2.2089 0.1825 50 1 5.8804 -0.2119 0.2616 51 1 5.7715 -1.1048 -1.2744 52 1 3.5640 -0.9686 -0.1141 53 1 3.5370 -0.1480 -1.6947 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 ZINC001281940670.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:20 Heat of formation + Delta-G solvation = -54.112279 kcal Electronic energy + Delta-G solvation = -35422.786023 eV Core-core repulsion = 31297.477279 eV Total energy + Delta-G solvation = -4125.308743 eV No. of doubly occupied orbitals = 66 Molecular weight (most abundant/longest-lived isotopes) = 340.220 amu Computer time = 0.82 seconds Orbital eigenvalues (eV) -40.71726 -38.41875 -37.60463 -36.43292 -34.95515 -33.62347 -31.40049 -31.21290 -30.39190 -28.94034 -27.96446 -26.07501 -25.63038 -24.90223 -23.65402 -23.13669 -22.31960 -21.37112 -20.69725 -19.59318 -17.93948 -16.94833 -16.34357 -16.22863 -15.90738 -15.83170 -15.55978 -15.41701 -14.81156 -14.53550 -14.44086 -14.01350 -13.64176 -13.47695 -13.24271 -13.12869 -12.77697 -12.64253 -12.33206 -12.16554 -11.98365 -11.94024 -11.75718 -11.52884 -11.34939 -11.22539 -11.16493 -11.01750 -10.84827 -10.64058 -10.59745 -10.45868 -10.22401 -10.10337 -9.97395 -9.84306 -9.79013 -9.65178 -9.43146 -9.12115 -8.41158 -8.30920 -7.64777 -6.07779 -5.44411 -3.16616 3.35992 3.40177 3.45255 3.52296 3.78676 4.00389 4.24327 4.56508 4.82265 4.89173 4.96206 5.05122 5.06665 5.22686 5.30025 5.37388 5.43652 5.54193 5.57867 5.59023 5.64195 5.67698 5.70399 5.72951 5.77477 5.84628 5.88137 5.94066 6.00871 6.06547 6.11320 6.21773 6.24564 6.26711 6.27975 6.44162 6.46551 6.55600 6.55953 6.61619 6.68297 6.98086 7.17444 7.20434 7.28488 7.57571 7.65656 7.87443 8.16935 8.51078 8.82520 9.34446 9.46276 10.06409 10.66099 11.45468 Molecular weight = 340.22amu Principal moments of inertia in cm(-1) A = 0.009360 B = 0.007146 C = 0.005997 Principal moments of inertia in units of 10**(-40)*gram-cm**2 A = 2990.653016 B = 3917.554638 C = 4667.514744 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.383 6.383 3 C 0.070 3.930 4 C -0.066 4.066 5 C 0.526 3.474 6 O -0.449 6.449 7 N -0.609 5.609 8 C 0.102 3.898 9 C -0.137 4.137 10 C -0.126 4.126 11 C 0.117 3.883 12 H 0.057 0.943 13 C 0.141 3.859 14 N -0.572 5.572 15 C 0.131 3.869 16 C -0.201 4.201 17 H 0.048 0.952 18 C 0.003 3.997 19 N -0.900 5.900 20 Si 0.884 3.116 21 H -0.291 1.291 22 H -0.291 1.291 23 H -0.278 1.278 24 C -0.129 4.129 25 C 0.052 3.948 26 O -0.391 6.391 27 C 0.061 3.939 28 C -0.141 4.141 29 H 0.037 0.963 30 H 0.083 0.917 31 H 0.041 0.959 32 H 0.066 0.934 33 H 0.051 0.949 34 H 0.065 0.935 35 H 0.086 0.914 36 H 0.066 0.934 37 H 0.071 0.929 38 H 0.069 0.931 39 H 0.071 0.929 40 H 0.046 0.954 41 H 0.032 0.968 42 H -0.008 1.008 43 H 0.098 0.902 44 H -0.011 1.011 45 H 0.249 0.751 46 H 0.091 0.909 47 H 0.064 0.936 48 H 0.086 0.914 49 H 0.088 0.912 50 H 0.048 0.952 51 H 0.087 0.913 52 H 0.093 0.907 53 H 0.073 0.927 For the charges calculated by CM2: Dipole moment (debyes) X Y Z Total from point charges -4.496 -14.783 -6.013 16.581 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.302 6.302 3 C -0.006 4.006 4 C -0.069 4.069 5 C 0.312 3.688 6 O -0.318 6.318 7 N -0.347 5.347 8 C -0.021 4.021 9 C -0.175 4.175 10 C -0.165 4.165 11 C 0.011 3.989 12 H 0.076 0.924 13 C 0.012 3.988 14 N -0.292 5.292 15 C -0.014 4.014 16 C -0.332 4.332 17 H 0.066 0.934 18 C -0.194 4.194 19 N -0.546 5.546 20 Si 0.715 3.285 21 H -0.218 1.218 22 H -0.217 1.217 23 H -0.204 1.204 24 C -0.167 4.167 25 C -0.024 4.024 26 O -0.310 6.310 27 C -0.016 4.016 28 C -0.178 4.178 29 H 0.055 0.945 30 H 0.101 0.899 31 H 0.060 0.940 32 H 0.084 0.916 33 H 0.069 0.931 34 H 0.083 0.917 35 H 0.104 0.896 36 H 0.085 0.915 37 H 0.090 0.910 38 H 0.088 0.912 39 H 0.090 0.910 40 H 0.064 0.936 41 H 0.050 0.950 42 H 0.010 0.990 43 H 0.116 0.884 44 H 0.008 0.992 45 H 0.060 0.940 46 H 0.109 0.891 47 H 0.083 0.917 48 H 0.104 0.896 49 H 0.106 0.894 50 H 0.066 0.934 51 H 0.105 0.895 52 H 0.111 0.889 53 H 0.092 0.908 Dipole moment (debyes) X Y Z Total from point charges -3.705 -14.186 -6.041 15.857 hybrid contribution -1.179 0.844 -0.454 1.519 sum -4.883 -13.341 -6.495 15.621 Atomic orbital electron populations 1.22679 0.81781 1.02244 1.00049 1.87964 1.19068 1.27879 1.95259 1.22597 0.92163 0.84511 1.01345 1.20810 0.93957 0.97383 0.94769 1.19291 0.83939 0.80525 0.85044 1.90797 1.29995 1.34069 1.76962 1.48128 1.48520 1.30353 1.07723 1.22753 0.94331 0.87012 0.98009 1.22849 1.01388 0.97346 0.95940 1.22634 1.00870 0.95812 0.97198 1.22366 0.98600 0.96634 0.81286 0.92442 1.21179 0.96400 0.87446 0.93738 1.43409 1.11316 1.31192 1.43318 1.21178 0.84460 1.03197 0.92555 1.18798 1.21465 0.87263 1.05676 0.93441 1.25256 0.99799 0.96991 0.97395 1.69708 1.43098 1.06496 1.35290 0.86607 0.79735 0.80133 0.82072 1.21830 1.21684 1.20387 1.22050 0.96993 0.98539 0.99121 1.22868 0.94208 0.84910 1.00366 1.87826 1.77089 1.14482 1.51595 1.22528 0.93022 0.89020 0.97005 1.22322 0.96559 0.97409 1.01525 0.94456 0.89860 0.94004 0.91580 0.93104 0.91737 0.89644 0.91506 0.91017 0.91234 0.91026 0.93561 0.94966 0.98963 0.88418 0.99234 0.94025 0.89072 0.91732 0.89633 0.89427 0.93370 0.89510 0.88877 0.90824 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.03 0.30 10.58 101.05 1.07 1.37 16 2 O -0.38 -4.98 5.33 -35.23 -0.19 -5.17 16 3 C 0.07 1.11 4.90 37.16 0.18 1.29 16 4 C -0.07 -1.10 0.49 -155.48 -0.08 -1.18 16 5 C 0.53 10.20 2.53 -10.98 -0.03 10.17 16 6 O -0.45 -10.43 7.51 5.56 0.04 -10.39 16 7 N -0.61 -10.61 3.04 -162.88 -0.49 -11.10 16 8 C 0.10 1.39 4.66 -2.41 -0.01 1.38 16 9 C -0.14 -1.58 7.10 -24.46 -0.17 -1.75 16 10 C -0.13 -1.66 6.48 -24.89 -0.16 -1.82 16 11 C 0.12 2.07 3.28 -66.95 -0.22 1.85 16 12 H 0.06 1.12 8.14 -51.93 -0.42 0.70 16 13 C 0.14 2.99 5.34 -4.04 -0.02 2.97 16 14 N -0.57 -14.47 1.62 -181.12 -0.29 -14.76 16 15 C 0.13 3.43 9.38 59.85 0.56 3.99 16 16 C -0.20 -5.36 6.39 -15.06 -0.10 -5.46 16 17 H 0.05 1.24 4.22 -52.48 -0.22 1.02 16 18 C 0.00 0.07 5.30 -17.43 -0.09 -0.02 16 19 N -0.90 -26.67 10.75 55.49 0.60 -26.07 16 20 Si 0.88 21.74 29.46 -169.99 -5.01 16.73 16 21 H -0.29 -7.19 7.11 56.52 0.40 -6.78 16 22 H -0.29 -7.31 7.11 56.52 0.40 -6.91 16 23 H -0.28 -6.45 6.50 56.52 0.37 -6.09 16 24 C -0.13 -2.30 4.42 -26.59 -0.12 -2.42 16 25 C 0.05 0.94 6.55 37.21 0.24 1.18 16 26 O -0.39 -6.66 10.77 -35.23 -0.38 -7.04 16 27 C 0.06 0.91 6.62 37.21 0.25 1.16 16 28 C -0.14 -2.01 4.09 -26.59 -0.11 -2.11 16 29 H 0.04 0.36 8.14 -51.93 -0.42 -0.06 16 30 H 0.08 0.68 8.14 -51.93 -0.42 0.25 16 31 H 0.04 0.45 8.14 -51.93 -0.42 0.03 16 32 H 0.07 1.22 6.06 -51.93 -0.31 0.91 16 33 H 0.05 0.76 7.75 -51.93 -0.40 0.36 16 34 H 0.06 0.91 6.35 -51.93 -0.33 0.58 16 35 H 0.09 1.12 3.69 -51.93 -0.19 0.93 16 36 H 0.07 0.74 8.14 -51.93 -0.42 0.32 16 37 H 0.07 0.68 8.14 -51.93 -0.42 0.26 16 38 H 0.07 0.82 8.14 -51.92 -0.42 0.39 16 39 H 0.07 0.92 7.80 -51.93 -0.41 0.52 16 40 H 0.05 1.00 7.90 -51.93 -0.41 0.59 16 41 H 0.03 0.67 8.06 -51.93 -0.42 0.25 16 42 H -0.01 -0.20 6.84 -51.93 -0.36 -0.56 16 43 H 0.10 2.88 5.14 -51.93 -0.27 2.61 16 44 H -0.01 -0.28 8.05 -51.93 -0.42 -0.70 16 45 H 0.25 7.20 8.31 -40.82 -0.34 6.87 16 46 H 0.09 1.90 7.39 -51.93 -0.38 1.52 16 47 H 0.06 1.03 8.06 -51.93 -0.42 0.61 16 48 H 0.09 1.40 8.14 -51.93 -0.42 0.98 16 49 H 0.09 1.80 5.64 -51.93 -0.29 1.51 16 50 H 0.05 0.79 7.38 -51.93 -0.38 0.41 16 51 H 0.09 1.07 8.14 -51.93 -0.42 0.65 16 52 H 0.09 1.24 4.22 -51.93 -0.22 1.02 16 53 H 0.07 1.00 8.14 -51.93 -0.42 0.58 16 LS Contribution 371.56 15.07 5.60 5.60 Total: -1.00 -31.08 371.56 -7.76 -38.84 By element: Atomic # 1 Polarization: 11.59 SS G_CDS: -8.83 Total: 2.76 kcal Atomic # 6 Polarization: 9.41 SS G_CDS: 1.20 Total: 10.61 kcal Atomic # 7 Polarization: -51.74 SS G_CDS: -0.19 Total: -51.94 kcal Atomic # 8 Polarization: -22.08 SS G_CDS: -0.53 Total: -22.60 kcal Atomic # 14 Polarization: 21.74 SS G_CDS: -5.01 Total: 16.73 kcal Total LS contribution 5.60 Total: 5.60 kcal Total: -31.08 -7.76 -38.84 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. ZINC001281940670.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 -15.275 kcal (2) G-P(sol) polarization free energy of solvation -31.079 kcal (3) G-ENP(sol) elect.-nuc.-pol. free energy of system -46.355 kcal (4) G-CDS(sol) cavity-dispersion-solvent structure free energy -7.757 kcal (5) G-P-CDS(sol) = G-P(sol) + G-CDS(sol) = (2) + (4) -38.837 kcal (6) G-S(sol) free energy of system = (1) + (5) -54.112 kcal 1SCF WAS SPECIFIED, SO GEOMETRY OPTIMIZATION WAS NOT USED Total computer time = 0.82 seconds