Wall clock time and date at job start Sat Apr 11 2020 06:13:57 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.42904 * 1 3 3 C 1.42895 * 114.00172 * 2 1 4 4 H 1.09002 * 109.52258 * 329.87151 * 3 2 1 5 5 C 1.53096 * 109.52388 * 89.99751 * 3 2 1 6 6 O 1.42900 * 109.41624 * 62.35817 * 5 3 2 7 7 C 1.42905 * 114.09687 * 61.15982 * 6 5 3 8 8 C 1.53096 * 109.41308 * 298.84159 * 7 6 5 9 9 C 1.53175 * 109.21377 * 57.57874 * 8 7 6 10 10 H 1.09004 * 109.54889 * 62.91724 * 9 8 7 11 11 C 1.53006 * 109.54231 * 183.18353 * 9 8 7 12 12 C 1.50702 * 109.47031 * 294.99545 * 11 9 8 13 13 O 1.21289 * 119.99779 * 0.02562 * 12 11 9 14 14 N 1.34770 * 120.00167 * 179.97438 * 12 11 9 15 15 C 1.46502 * 120.00338 * 174.55649 * 14 12 11 16 16 C 1.52997 * 110.43534 * 303.00891 * 15 14 12 17 17 N 1.46500 * 109.47097 * 307.76108 * 16 15 14 18 18 C 1.36941 * 120.00316 * 174.95670 * 17 16 15 19 19 H 1.07994 * 119.99402 * 0.02562 * 18 17 16 20 20 C 1.38807 * 120.00209 * 179.97438 * 18 17 16 21 21 N 1.31627 * 119.99886 * 0.02562 * 20 18 17 22 22 Si 1.86795 * 120.00018 * 180.02562 * 20 18 17 23 23 H 1.48500 * 109.47497 * 90.00446 * 22 20 18 24 24 H 1.48507 * 109.47312 * 210.00452 * 22 20 18 25 25 H 1.48507 * 109.47159 * 330.00338 * 22 20 18 26 26 C 1.54883 * 110.48701 * 180.46461 * 15 14 12 27 27 C 1.55237 * 102.58797 * 279.18945 * 26 15 14 28 28 C 1.54269 * 104.23318 * 322.08590 * 27 26 15 29 29 C 1.53873 * 106.59532 * 23.60791 * 28 27 26 30 30 H 1.09006 * 109.46632 * 59.99920 * 1 2 3 31 31 H 1.08999 * 109.47399 * 179.97438 * 1 2 3 32 32 H 1.08993 * 109.47147 * 300.00455 * 1 2 3 33 33 H 1.09008 * 109.47505 * 302.37376 * 5 3 2 34 34 H 1.08999 * 109.48309 * 182.34435 * 5 3 2 35 35 H 1.08998 * 109.48063 * 178.86260 * 7 6 5 36 36 H 1.08998 * 109.48197 * 58.82059 * 7 6 5 37 37 H 1.09002 * 109.51874 * 177.51375 * 8 7 6 38 38 H 1.08995 * 109.52353 * 297.63940 * 8 7 6 39 39 H 1.08997 * 109.47103 * 54.99669 * 11 9 8 40 40 H 1.09006 * 109.47016 * 175.00159 * 11 9 8 41 41 H 0.97002 * 119.99692 * 354.56571 * 14 12 11 42 42 H 1.08997 * 109.47452 * 67.75665 * 16 15 14 43 43 H 1.08997 * 109.47041 * 187.76075 * 16 15 14 44 44 H 0.97002 * 119.99923 * 354.95490 * 17 16 15 45 45 H 0.97000 * 119.99374 * 179.97438 * 21 20 18 46 46 H 1.08999 * 110.75992 * 37.58260 * 26 15 14 47 47 H 1.09005 * 110.75776 * 160.95340 * 26 15 14 48 48 H 1.09002 * 110.46951 * 203.39438 * 27 26 15 49 49 H 1.08995 * 110.58383 * 80.83309 * 27 26 15 50 50 H 1.09001 * 110.03300 * 264.33068 * 28 27 26 51 51 H 1.08999 * 110.07345 * 142.92101 * 28 27 26 52 52 H 1.09000 * 110.03353 * 240.71393 * 29 28 27 53 53 H 1.09002 * 110.06979 * 119.31177 * 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 8 1.4290 0.0000 0.0000 3 6 2.0103 1.3054 0.0000 4 1 1.3467 1.9996 -0.5157 5 6 2.2184 1.7728 1.4429 6 8 3.1339 0.8939 2.0997 7 6 4.4315 0.8508 1.5026 8 6 4.3124 0.3463 0.0621 9 6 3.3637 1.2631 -0.7162 10 1 3.7831 2.2683 -0.7612 11 6 3.1769 0.7207 -2.1347 12 6 4.4884 0.7855 -2.8741 13 8 5.4755 1.2221 -2.3208 14 7 4.5648 0.3574 -4.1497 15 6 5.8069 0.5297 -4.9072 16 6 6.2085 2.0054 -4.9460 17 7 5.0765 2.8074 -5.4168 18 6 5.1663 4.1737 -5.4362 19 1 6.0718 4.6600 -5.1048 20 6 4.0939 4.9335 -5.8827 21 7 2.9904 4.3407 -6.2870 22 14 4.2167 6.7972 -5.9098 23 1 4.7552 7.2389 -7.2213 24 1 2.8685 7.3852 -5.7049 25 1 5.1240 7.2494 -4.8245 26 6 5.6473 -0.0268 -6.3438 27 6 5.8549 -1.5518 -6.1409 28 6 6.9674 -1.6448 -5.0762 29 6 6.9368 -0.3214 -4.2918 30 1 -0.3633 0.5139 -0.8900 31 1 -0.3634 -1.0276 -0.0005 32 1 -0.3633 0.5139 0.8899 33 1 1.2638 1.7631 1.9693 34 1 2.6230 2.7849 1.4424 35 1 5.0714 0.1767 2.0720 36 1 4.8651 1.8508 1.5028 37 1 5.2952 0.3546 -0.4092 38 1 3.9176 -0.6696 0.0634 39 1 2.8378 -0.3141 -2.0862 40 1 2.4342 1.3223 -2.6588 41 1 3.7972 -0.0647 -4.5663 42 1 6.4932 2.3310 -3.9456 43 1 7.0522 2.1346 -5.6240 44 1 4.2633 2.3705 -5.7148 45 1 2.2409 4.8717 -6.5987 46 1 4.6497 0.1818 -6.7303 47 1 6.4120 0.3829 -7.0039 48 1 6.1770 -2.0197 -7.0711 49 1 4.9392 -2.0182 -5.7775 50 1 7.9366 -1.7704 -5.5589 51 1 6.7737 -2.4814 -4.4049 52 1 7.8912 0.1952 -4.3938 53 1 6.7280 -0.5149 -3.2397 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 ZINC001568001089.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:13:57 Heat of formation + Delta-G solvation = -129.153188 kcal Electronic energy + Delta-G solvation = -32328.380840 eV Core-core repulsion = 28199.818079 eV Total energy + Delta-G solvation = -4128.562761 eV No. of doubly occupied orbitals = 66 Molecular weight (most abundant/longest-lived isotopes) = 340.220 amu Computer time = 1.00 seconds Orbital eigenvalues (eV) -40.34257 -39.55390 -37.49452 -36.98655 -34.32948 -33.93699 -32.96027 -31.06617 -30.49855 -29.06451 -28.52843 -27.68302 -26.66331 -24.50568 -23.61779 -23.33224 -22.16435 -20.65889 -20.25064 -19.73704 -18.87570 -17.06593 -16.91125 -16.58534 -16.41900 -15.88566 -15.61423 -15.50469 -14.99173 -14.92690 -14.57037 -14.47062 -14.16907 -13.91474 -13.70143 -13.22399 -13.06364 -12.82834 -12.70836 -12.57309 -12.43222 -12.18804 -11.96852 -11.67444 -11.38859 -11.29956 -11.14693 -11.03904 -10.89985 -10.87124 -10.82120 -10.73412 -10.63698 -10.33728 -10.29094 -10.25236 -10.04503 -9.99610 -9.67109 -9.27632 -9.01750 -8.74581 -8.41336 -6.83259 -5.73469 -2.99090 3.19811 3.44407 3.50022 3.51765 3.60251 3.65138 3.99768 4.60861 4.64665 4.73090 4.77751 4.85383 4.89261 4.94641 5.04124 5.17423 5.25300 5.29687 5.38610 5.42362 5.47085 5.54450 5.56003 5.58615 5.62353 5.64528 5.65463 5.70909 5.75481 5.79667 5.93793 5.98887 5.99253 6.02226 6.05942 6.11247 6.23358 6.30472 6.35119 6.41592 6.46234 6.50071 6.60913 6.65665 6.97537 7.02948 7.21689 7.29724 7.74967 8.04647 8.46546 8.80076 9.57991 10.16985 10.53205 11.49787 Molecular weight = 340.22amu Principal moments of inertia in cm(-1) A = 0.007987 B = 0.004285 C = 0.003241 Principal moments of inertia in units of 10**(-40)*gram-cm**2 A = 3504.691555 B = 6533.512468 C = 8638.425435 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.369 6.369 3 C 0.057 3.943 4 H 0.071 0.929 5 C 0.043 3.957 6 O -0.377 6.377 7 C 0.061 3.939 8 C -0.148 4.148 9 C -0.064 4.064 10 H 0.092 0.908 11 C -0.131 4.131 12 C 0.527 3.473 13 O -0.536 6.536 14 N -0.710 5.710 15 C 0.170 3.830 16 C 0.146 3.854 17 N -0.724 5.724 18 C -0.227 4.227 19 H 0.067 0.933 20 C -0.016 4.016 21 N -0.937 5.937 22 Si 0.903 3.097 23 H -0.291 1.291 24 H -0.293 1.293 25 H -0.283 1.283 26 C -0.119 4.119 27 C -0.125 4.125 28 C -0.122 4.122 29 C -0.123 4.123 30 H 0.041 0.959 31 H 0.083 0.917 32 H 0.038 0.962 33 H 0.091 0.909 34 H 0.059 0.941 35 H 0.090 0.910 36 H 0.049 0.951 37 H 0.094 0.906 38 H 0.069 0.931 39 H 0.093 0.907 40 H 0.093 0.907 41 H 0.394 0.606 42 H 0.045 0.955 43 H 0.018 0.982 44 H 0.382 0.618 45 H 0.252 0.748 46 H 0.075 0.925 47 H 0.082 0.918 48 H 0.062 0.938 49 H 0.062 0.938 50 H 0.059 0.941 51 H 0.057 0.943 52 H 0.069 0.931 53 H 0.073 0.927 For the charges calculated by CM2: Dipole moment (debyes) X Y Z Total from point charges -1.543 -13.510 9.875 16.805 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.288 6.288 3 C -0.001 4.001 4 H 0.089 0.911 5 C -0.035 4.035 6 O -0.296 6.296 7 C -0.015 4.015 8 C -0.186 4.186 9 C -0.083 4.083 10 H 0.110 0.890 11 C -0.171 4.171 12 C 0.314 3.686 13 O -0.413 6.413 14 N -0.360 5.360 15 C 0.084 3.916 16 C 0.017 3.983 17 N -0.368 5.368 18 C -0.358 4.358 19 H 0.085 0.915 20 C -0.212 4.212 21 N -0.587 5.587 22 Si 0.735 3.265 23 H -0.219 1.219 24 H -0.219 1.219 25 H -0.209 1.209 26 C -0.157 4.157 27 C -0.162 4.162 28 C -0.160 4.160 29 C -0.162 4.162 30 H 0.060 0.940 31 H 0.102 0.898 32 H 0.056 0.944 33 H 0.109 0.891 34 H 0.078 0.922 35 H 0.108 0.892 36 H 0.067 0.933 37 H 0.112 0.888 38 H 0.088 0.912 39 H 0.111 0.889 40 H 0.111 0.889 41 H 0.228 0.772 42 H 0.063 0.937 43 H 0.036 0.964 44 H 0.215 0.785 45 H 0.062 0.938 46 H 0.093 0.907 47 H 0.100 0.900 48 H 0.081 0.919 49 H 0.080 0.920 50 H 0.078 0.922 51 H 0.076 0.924 52 H 0.087 0.913 53 H 0.092 0.908 Dipole moment (debyes) X Y Z Total from point charges -0.688 -13.283 10.459 16.920 hybrid contribution -0.945 1.246 -0.754 1.736 sum -1.634 -12.037 9.705 15.548 Atomic orbital electron populations 1.22791 0.81676 1.02560 0.99782 1.87961 1.18939 1.26941 1.94929 1.22897 0.95571 0.85600 0.96069 0.91078 1.23019 0.94125 0.93242 0.93079 1.87906 1.22317 1.65570 1.53766 1.22609 0.84316 1.00646 0.93966 1.22263 1.01616 0.99496 0.95182 1.21000 0.93168 1.01201 0.92963 0.88960 1.21599 0.96405 1.03651 0.95482 1.20317 0.90255 0.76457 0.81565 1.90597 1.37935 1.45472 1.67305 1.45818 1.15254 1.61998 1.12967 1.21321 0.83979 0.96039 0.90240 1.21052 0.91340 0.87662 0.98200 1.42196 1.14142 1.02886 1.77602 1.19238 0.98900 0.82858 1.34851 0.91464 1.24654 0.96330 0.97617 1.02554 1.69653 1.04584 1.33921 1.50540 0.86232 0.76564 0.86063 0.77647 1.21857 1.21921 1.20882 1.22789 1.03101 0.94246 0.95599 1.22451 0.98587 0.96673 0.98527 1.22128 0.98685 0.97731 0.97448 1.22433 0.98591 0.94855 1.00300 0.93990 0.89815 0.94383 0.89132 0.92237 0.89166 0.93261 0.88786 0.91185 0.88911 0.88914 0.77239 0.93700 0.96430 0.78481 0.93790 0.90675 0.89956 0.91918 0.91972 0.92170 0.92426 0.91253 0.90797 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.03 0.26 10.82 101.05 1.09 1.36 16 2 O -0.37 -4.53 9.07 -42.66 -0.39 -4.92 16 3 C 0.06 0.69 2.43 -26.59 -0.06 0.62 16 4 H 0.07 0.81 7.79 -51.93 -0.40 0.40 16 5 C 0.04 0.50 6.49 37.21 0.24 0.74 16 6 O -0.38 -5.09 10.70 -42.66 -0.46 -5.54 16 7 C 0.06 0.76 6.84 37.21 0.25 1.02 16 8 C -0.15 -1.98 4.96 -26.59 -0.13 -2.12 16 9 C -0.06 -0.89 1.76 -90.43 -0.16 -1.05 16 10 H 0.09 1.44 7.97 -51.93 -0.41 1.03 16 11 C -0.13 -1.84 4.94 -27.88 -0.14 -1.98 16 12 C 0.53 8.75 6.81 -10.99 -0.07 8.68 16 13 O -0.54 -10.29 10.94 5.55 0.06 -10.23 16 14 N -0.71 -10.98 4.54 -46.08 -0.21 -11.19 16 15 C 0.17 2.77 0.39 -130.16 -0.05 2.72 16 16 C 0.15 3.07 4.72 -4.04 -0.02 3.05 16 17 N -0.72 -18.15 5.30 -59.91 -0.32 -18.46 16 18 C -0.23 -6.42 9.99 -15.06 -0.15 -6.57 16 19 H 0.07 1.90 7.83 -52.49 -0.41 1.49 16 20 C -0.02 -0.47 5.50 -17.43 -0.10 -0.57 16 21 N -0.94 -27.85 13.06 55.49 0.72 -27.13 16 22 Si 0.90 22.67 29.55 -169.99 -5.02 17.65 16 23 H -0.29 -7.27 7.11 56.52 0.40 -6.87 16 24 H -0.29 -7.37 7.11 56.52 0.40 -6.97 16 25 H -0.28 -7.00 7.11 56.52 0.40 -6.60 16 26 C -0.12 -1.72 5.67 -24.55 -0.14 -1.86 16 27 C -0.12 -1.44 6.80 -24.88 -0.17 -1.61 16 28 C -0.12 -1.44 6.80 -25.61 -0.17 -1.61 16 29 C -0.12 -1.78 5.68 -25.61 -0.15 -1.93 16 30 H 0.04 0.36 8.09 -51.93 -0.42 -0.06 16 31 H 0.08 0.66 8.14 -51.93 -0.42 0.24 16 32 H 0.04 0.31 7.83 -51.93 -0.41 -0.09 16 33 H 0.09 0.91 7.83 -51.92 -0.41 0.50 16 34 H 0.06 0.65 8.14 -51.93 -0.42 0.22 16 35 H 0.09 0.99 8.14 -51.93 -0.42 0.57 16 36 H 0.05 0.60 8.14 -51.93 -0.42 0.18 16 37 H 0.09 1.40 6.03 -51.93 -0.31 1.08 16 38 H 0.07 0.96 7.84 -51.93 -0.41 0.55 16 39 H 0.09 1.16 7.77 -51.93 -0.40 0.76 16 40 H 0.09 1.34 8.14 -51.93 -0.42 0.92 16 41 H 0.39 5.30 8.50 -40.82 -0.35 4.96 16 42 H 0.05 1.06 6.49 -51.93 -0.34 0.72 16 43 H 0.02 0.38 8.02 -51.93 -0.42 -0.03 16 44 H 0.38 9.81 7.27 -40.82 -0.30 9.51 16 45 H 0.25 7.30 8.31 -40.82 -0.34 6.96 16 46 H 0.07 1.14 8.10 -51.93 -0.42 0.72 16 47 H 0.08 1.24 7.97 -51.93 -0.41 0.82 16 48 H 0.06 0.65 8.14 -51.93 -0.42 0.23 16 49 H 0.06 0.67 8.14 -51.93 -0.42 0.25 16 50 H 0.06 0.65 8.14 -51.93 -0.42 0.23 16 51 H 0.06 0.63 8.14 -51.93 -0.42 0.21 16 52 H 0.07 1.02 7.84 -51.93 -0.41 0.61 16 53 H 0.07 1.16 6.63 -51.93 -0.34 0.81 16 LS Contribution 406.46 15.07 6.13 6.13 Total: -1.00 -32.54 406.46 -8.91 -41.46 By element: Atomic # 1 Polarization: 22.85 SS G_CDS: -9.51 Total: 13.34 kcal Atomic # 6 Polarization: -1.17 SS G_CDS: 0.08 Total: -1.10 kcal Atomic # 7 Polarization: -56.98 SS G_CDS: 0.20 Total: -56.79 kcal Atomic # 8 Polarization: -19.91 SS G_CDS: -0.78 Total: -20.69 kcal Atomic # 14 Polarization: 22.67 SS G_CDS: -5.02 Total: 17.65 kcal Total LS contribution 6.13 Total: 6.13 kcal Total: -32.54 -8.91 -41.46 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. ZINC001568001089.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 -87.697 kcal (2) G-P(sol) polarization free energy of solvation -32.541 kcal (3) G-ENP(sol) elect.-nuc.-pol. free energy of system -120.239 kcal (4) G-CDS(sol) cavity-dispersion-solvent structure free energy -8.914 kcal (5) G-P-CDS(sol) = G-P(sol) + G-CDS(sol) = (2) + (4) -41.456 kcal (6) G-S(sol) free energy of system = (1) + (5) -129.153 kcal 1SCF WAS SPECIFIED, SO GEOMETRY OPTIMIZATION WAS NOT USED Total computer time = 1.00 seconds