Wall clock time and date at job start Tue Mar 31 2020 23:34:49 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.53004 * 1 3 3 H 1.08997 * 110.67634 * 2 1 4 4 C 1.54997 * 110.61638 * 237.03848 * 2 1 3 5 5 C 1.54815 * 102.09311 * 204.47396 * 4 2 1 6 6 C 1.54093 * 104.18595 * 20.88073 * 5 4 2 7 7 H 1.08996 * 109.87999 * 242.45739 * 6 5 4 8 8 C 1.53006 * 109.87918 * 121.45377 * 6 5 4 9 9 C 1.52999 * 109.46997 * 175.00095 * 8 6 5 10 10 C 1.50706 * 109.46573 * 180.02562 * 9 8 6 11 11 O 1.21274 * 120.00042 * 359.72068 * 10 9 8 12 12 N 1.34777 * 119.99766 * 179.97438 * 10 9 8 13 13 C 1.46929 * 120.63091 * 359.72591 * 12 10 9 14 14 C 1.52765 * 109.00690 * 233.58011 * 13 12 10 15 15 C 1.53032 * 109.38265 * 305.58660 * 14 13 12 16 16 O 1.42894 * 109.46125 * 301.21882 * 15 14 13 17 17 C 1.53001 * 109.44971 * 181.25365 * 15 14 13 18 18 N 1.46492 * 109.47193 * 59.93832 * 17 15 14 19 19 C 1.36939 * 120.00194 * 180.02562 * 18 17 15 20 20 H 1.07996 * 119.99667 * 359.97438 * 19 18 17 21 21 C 1.38813 * 120.00016 * 180.02562 * 19 18 17 22 22 N 1.31615 * 120.00140 * 359.97438 * 21 19 18 23 23 Si 1.86802 * 120.00057 * 179.97438 * 21 19 18 24 24 H 1.48508 * 109.46797 * 0.02562 * 23 21 19 25 25 H 1.48500 * 109.47108 * 119.99582 * 23 21 19 26 26 H 1.48500 * 109.47425 * 240.00133 * 23 21 19 27 27 C 1.53041 * 109.53897 * 61.24115 * 15 14 13 28 28 C 1.46923 * 120.63281 * 179.97438 * 12 10 9 29 29 O 1.43748 * 107.29498 * 2.03042 * 6 5 4 30 30 H 1.09000 * 109.46860 * 182.96193 * 1 2 3 31 31 H 1.09000 * 109.46791 * 302.96060 * 1 2 3 32 32 H 1.08991 * 109.47055 * 62.95995 * 1 2 3 33 33 H 1.09010 * 110.90716 * 86.27556 * 4 2 1 34 34 H 1.09001 * 110.90870 * 322.67791 * 4 2 1 35 35 H 1.08996 * 110.48449 * 262.20354 * 5 4 2 36 36 H 1.08993 * 110.48430 * 139.55788 * 5 4 2 37 37 H 1.09000 * 109.46866 * 54.99883 * 8 6 5 38 38 H 1.09001 * 109.46611 * 295.00000 * 8 6 5 39 39 H 1.09000 * 109.47486 * 60.00167 * 9 8 6 40 40 H 1.09001 * 109.47362 * 299.99445 * 9 8 6 41 41 H 1.09004 * 109.58660 * 353.40881 * 13 12 10 42 42 H 1.09001 * 109.59243 * 113.83779 * 13 12 10 43 43 H 1.09001 * 109.56725 * 185.54838 * 14 13 12 44 44 H 1.09006 * 109.38446 * 65.51568 * 14 13 12 45 45 H 0.96700 * 113.99983 * 300.06726 * 16 15 14 46 46 H 1.09006 * 109.47187 * 179.97438 * 17 15 14 47 47 H 1.09006 * 109.46949 * 299.93298 * 17 15 14 48 48 H 0.97007 * 120.00253 * 359.97438 * 18 17 15 49 49 H 0.97004 * 120.00209 * 180.02562 * 22 21 19 50 50 H 1.09008 * 109.49661 * 178.68278 * 27 15 14 51 51 H 1.09000 * 109.52837 * 58.60129 * 27 15 14 52 52 H 1.08998 * 109.58754 * 246.16009 * 28 12 10 53 53 H 1.08998 * 109.58450 * 6.58545 * 28 12 10 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 1 1.9149 1.0198 0.0000 4 6 2.0758 -0.7893 -1.2172 5 6 3.4795 -1.2167 -0.7236 6 6 3.3840 -1.1255 0.8116 7 1 4.0873 -0.3802 1.1831 8 6 3.6936 -2.4885 1.4341 9 6 3.7188 -2.3591 2.9584 10 6 4.0243 -3.7016 3.5713 11 8 4.1992 -4.6676 2.8594 12 7 4.0982 -3.8295 4.9110 13 6 3.8891 -2.6683 5.7866 14 6 5.0827 -2.5359 6.7308 15 6 5.2896 -3.8563 7.4761 16 8 4.1059 -4.1805 8.2078 17 6 6.4653 -3.7156 8.4451 18 7 6.1795 -2.6492 9.4080 19 6 7.0945 -2.3399 10.3788 20 1 8.0266 -2.8825 10.4347 21 6 6.8234 -1.3298 11.2915 22 7 5.6872 -0.6690 11.2236 23 14 8.0712 -0.9084 12.6162 24 1 9.2579 -1.7905 12.4770 25 1 8.4889 0.5091 12.4699 26 1 7.4593 -1.1086 13.9544 27 6 5.5885 -4.9704 6.4703 28 6 4.3946 -5.1317 5.5234 29 8 2.0385 -0.7354 1.1336 30 1 -0.3633 -1.0263 0.0531 31 1 -0.3633 0.5591 0.8623 32 1 -0.3633 0.4672 -0.9153 33 1 2.1509 -0.1475 -2.0951 34 1 1.4549 -1.6605 -1.4260 35 1 4.2406 -0.5331 -1.0999 36 1 3.6974 -2.2390 -1.0323 37 1 4.6651 -2.8361 1.0828 38 1 2.9246 -3.2039 1.1427 39 1 2.7472 -2.0115 3.3097 40 1 4.4878 -1.6437 3.2497 41 1 3.8016 -1.7657 5.1818 42 1 2.9792 -2.8094 6.3699 43 1 4.8931 -1.7384 7.4492 44 1 5.9766 -2.3001 6.1531 45 1 3.8557 -3.5169 8.8652 46 1 6.6127 -4.6558 8.9765 47 1 7.3686 -3.4688 7.8870 48 1 5.3424 -2.1616 9.3576 49 1 5.4979 0.0371 11.8612 50 1 5.7589 -5.9058 7.0035 51 1 6.4775 -4.7127 5.8946 52 1 3.5267 -5.4775 6.0849 53 1 4.6404 -5.8540 4.7451 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 ZINC001700308518.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 23:34:49 Heat of formation + Delta-G solvation = -136.406751 kcal Electronic energy + Delta-G solvation = -30590.622836 eV Core-core repulsion = 26461.745537 eV Total energy + Delta-G solvation = -4128.877299 eV No. of doubly occupied orbitals = 66 Molecular weight (most abundant/longest-lived isotopes) = 340.220 amu Computer time = 1.08 seconds Orbital eigenvalues (eV) -40.41813 -39.83978 -38.46551 -35.49641 -34.18904 -34.00883 -32.56013 -32.06482 -31.41371 -30.60248 -30.30484 -27.28695 -26.24342 -24.09933 -23.32274 -23.06631 -21.75630 -21.49980 -21.15280 -20.12664 -18.81457 -17.55044 -17.16579 -16.89546 -16.37835 -15.93421 -15.63619 -15.28444 -15.17178 -14.94671 -14.35546 -14.22983 -14.12471 -13.99040 -13.68468 -13.34427 -13.26565 -13.09990 -12.91253 -12.70072 -12.56134 -12.43836 -12.26438 -12.06761 -11.89256 -11.85103 -11.52723 -11.34593 -11.16476 -11.13591 -11.05519 -11.03840 -10.97167 -10.87092 -10.86113 -10.28803 -10.12419 -9.91901 -9.73624 -9.67156 -9.39616 -8.97078 -8.50055 -7.02922 -5.85909 -3.14791 2.68265 3.05042 3.38727 3.45090 3.46825 3.47091 4.21207 4.27266 4.46828 4.51086 4.55751 4.57033 4.68198 4.80565 4.85706 4.94199 4.99847 5.06927 5.09733 5.17994 5.19781 5.22269 5.25552 5.37492 5.41845 5.44410 5.49631 5.50668 5.52072 5.64844 5.66616 5.67469 5.74211 5.80554 5.85502 5.91245 5.91580 6.00115 6.04959 6.06715 6.13468 6.34882 6.39212 6.60617 6.67749 6.69210 6.71205 7.48165 7.63902 8.08692 8.22068 8.33317 9.43975 9.97682 10.31860 11.36939 Molecular weight = 340.22amu Principal moments of inertia in cm(-1) A = 0.015338 B = 0.002136 C = 0.001979 Principal moments of inertia in units of 10**(-40)*gram-cm**2 A = 1825.041706 B =13104.859827 C =14146.830873 Net atomic charges, atomic populations, and dipole contributions calculated with CM2 Atom NO. Type Charge No. of electrons 1 C -0.147 4.147 2 C 0.070 3.930 3 H 0.069 0.931 4 C -0.151 4.151 5 C -0.153 4.153 6 C 0.076 3.924 7 H 0.070 0.930 8 C -0.099 4.099 9 C -0.138 4.138 10 C 0.514 3.486 11 O -0.553 6.553 12 N -0.607 5.607 13 C 0.111 3.889 14 C -0.158 4.158 15 C 0.130 3.870 16 O -0.560 6.560 17 C 0.129 3.871 18 N -0.729 5.729 19 C -0.240 4.240 20 H 0.071 0.929 21 C -0.013 4.013 22 N -0.934 5.934 23 Si 0.908 3.092 24 H -0.280 1.280 25 H -0.291 1.291 26 H -0.291 1.291 27 C -0.116 4.116 28 C 0.118 3.882 29 O -0.369 6.369 30 H 0.064 0.936 31 H 0.063 0.937 32 H 0.069 0.931 33 H 0.077 0.923 34 H 0.077 0.923 35 H 0.075 0.925 36 H 0.080 0.920 37 H 0.081 0.919 38 H 0.076 0.924 39 H 0.099 0.901 40 H 0.092 0.908 41 H 0.082 0.918 42 H 0.075 0.925 43 H 0.092 0.908 44 H 0.072 0.928 45 H 0.382 0.618 46 H 0.028 0.972 47 H 0.025 0.975 48 H 0.380 0.620 49 H 0.256 0.744 50 H 0.081 0.919 51 H 0.076 0.924 52 H 0.074 0.926 53 H 0.087 0.913 For the charges calculated by CM2: Dipole moment (debyes) X Y Z Total from point charges -8.130 -0.069 -24.220 25.549 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.204 4.204 2 C 0.013 3.987 3 H 0.086 0.914 4 C -0.189 4.189 5 C -0.191 4.191 6 C 0.018 3.982 7 H 0.087 0.913 8 C -0.137 4.137 9 C -0.178 4.178 10 C 0.303 3.697 11 O -0.432 6.432 12 N -0.340 5.340 13 C -0.012 4.012 14 C -0.196 4.196 15 C 0.089 3.911 16 O -0.368 6.368 17 C 0.001 3.999 18 N -0.375 5.375 19 C -0.371 4.371 20 H 0.089 0.911 21 C -0.209 4.209 22 N -0.582 5.582 23 Si 0.739 3.261 24 H -0.206 1.206 25 H -0.218 1.218 26 H -0.218 1.218 27 C -0.155 4.155 28 C -0.005 4.005 29 O -0.288 6.288 30 H 0.083 0.917 31 H 0.082 0.918 32 H 0.087 0.913 33 H 0.096 0.904 34 H 0.096 0.904 35 H 0.094 0.906 36 H 0.099 0.901 37 H 0.099 0.901 38 H 0.095 0.905 39 H 0.117 0.883 40 H 0.110 0.890 41 H 0.101 0.899 42 H 0.094 0.906 43 H 0.110 0.890 44 H 0.091 0.909 45 H 0.232 0.768 46 H 0.046 0.954 47 H 0.043 0.957 48 H 0.215 0.785 49 H 0.065 0.935 50 H 0.099 0.901 51 H 0.094 0.906 52 H 0.093 0.907 53 H 0.105 0.895 Dipole moment (debyes) X Y Z Total from point charges -7.979 -1.166 -24.668 25.953 hybrid contribution 0.664 1.188 1.166 1.791 sum -7.316 0.022 -23.503 24.615 Atomic orbital electron populations 1.21771 0.93205 1.03059 1.02365 1.23113 0.95204 0.95258 0.85124 0.91356 1.23184 0.96394 1.00835 0.98476 1.22993 0.97681 1.03092 0.95314 1.22789 0.84878 0.96481 0.94044 0.91255 1.21534 1.03239 0.95525 0.93401 1.21669 1.04310 0.95185 0.96617 1.20713 0.76241 0.91004 0.81720 1.90653 1.51034 1.40869 1.60624 1.48167 1.69847 1.10524 1.05473 1.21669 0.99136 0.87695 0.92650 1.22472 0.98520 0.99382 0.99249 1.21279 0.86084 0.92905 0.90834 1.86374 1.43876 1.51331 1.55187 1.21317 0.94737 0.92832 0.91050 1.42290 1.21564 1.39755 1.33890 1.19251 0.98090 1.13954 1.05766 0.91092 1.24731 0.97492 0.98943 0.99684 1.69637 1.23559 1.26058 1.38987 0.86102 0.80430 0.77830 0.81720 1.20565 1.21840 1.21803 1.21852 0.98165 0.99240 0.96238 1.21509 0.98105 0.84804 0.96034 1.88913 1.27702 1.73371 1.38852 0.91716 0.91779 0.91259 0.90410 0.90416 0.90638 0.90094 0.90076 0.90480 0.88254 0.88995 0.89932 0.90649 0.89005 0.90918 0.76834 0.95372 0.95711 0.78535 0.93463 0.90057 0.90561 0.90742 0.89486 Number of geometries 1 Number of calculations of the screened coulomb radii 1 The total number of SCF iterations 13. 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.95 9.66 37.16 0.36 -0.59 16 2 C 0.07 0.53 3.98 -25.67 -0.10 0.43 16 3 H 0.07 0.47 8.14 -51.93 -0.42 0.05 16 4 C -0.15 -0.96 6.59 -24.72 -0.16 -1.13 16 5 C -0.15 -1.05 6.64 -25.20 -0.17 -1.22 16 6 C 0.08 0.65 3.41 -26.16 -0.09 0.56 16 7 H 0.07 0.56 8.14 -51.93 -0.42 0.13 16 8 C -0.10 -0.98 3.99 -26.73 -0.11 -1.09 16 9 C -0.14 -1.45 3.96 -27.88 -0.11 -1.56 16 10 C 0.51 7.02 7.70 -10.98 -0.08 6.93 16 11 O -0.55 -9.17 15.39 5.56 0.09 -9.08 16 12 N -0.61 -8.08 2.98 -173.00 -0.52 -8.59 16 13 C 0.11 1.37 6.08 -3.93 -0.02 1.35 16 14 C -0.16 -2.28 4.91 -26.84 -0.13 -2.42 16 15 C 0.13 2.05 0.71 -90.58 -0.06 1.99 16 16 O -0.56 -9.22 11.71 -35.23 -0.41 -9.64 16 17 C 0.13 2.50 5.40 -4.05 -0.02 2.48 16 18 N -0.73 -17.27 5.07 -59.91 -0.30 -17.57 16 19 C -0.24 -6.49 9.99 -15.06 -0.15 -6.64 16 20 H 0.07 1.90 7.83 -52.49 -0.41 1.49 16 21 C -0.01 -0.35 5.50 -17.43 -0.10 -0.45 16 22 N -0.93 -27.01 13.06 55.49 0.72 -26.28 16 23 Si 0.91 21.96 29.55 -169.99 -5.02 16.94 16 24 H -0.28 -6.61 7.11 56.52 0.40 -6.21 16 25 H -0.29 -7.08 7.11 56.52 0.40 -6.68 16 26 H -0.29 -7.09 7.11 56.52 0.40 -6.69 16 27 C -0.12 -1.62 5.11 -26.59 -0.14 -1.76 16 28 C 0.12 1.61 6.20 -3.69 -0.02 1.59 16 29 O -0.37 -3.72 10.54 -35.23 -0.37 -4.10 16 30 H 0.06 0.45 8.14 -51.93 -0.42 0.03 16 31 H 0.06 0.43 8.14 -51.93 -0.42 0.01 16 32 H 0.07 0.33 8.14 -51.93 -0.42 -0.10 16 33 H 0.08 0.38 8.14 -51.92 -0.42 -0.05 16 34 H 0.08 0.51 8.06 -51.93 -0.42 0.09 16 35 H 0.07 0.43 8.14 -51.93 -0.42 0.01 16 36 H 0.08 0.55 7.70 -51.93 -0.40 0.15 16 37 H 0.08 0.80 8.10 -51.93 -0.42 0.38 16 38 H 0.08 0.87 8.10 -51.93 -0.42 0.45 16 39 H 0.10 1.03 7.72 -51.93 -0.40 0.63 16 40 H 0.09 0.85 7.68 -51.93 -0.40 0.45 16 41 H 0.08 0.86 5.93 -51.93 -0.31 0.55 16 42 H 0.08 0.95 7.81 -51.93 -0.41 0.54 16 43 H 0.09 1.48 6.82 -51.93 -0.35 1.13 16 44 H 0.07 1.06 8.14 -51.93 -0.42 0.64 16 45 H 0.38 6.30 7.85 45.56 0.36 6.65 16 46 H 0.03 0.58 8.14 -51.93 -0.42 0.16 16 47 H 0.03 0.49 8.14 -51.93 -0.42 0.07 16 48 H 0.38 9.12 5.58 -40.82 -0.23 8.89 16 49 H 0.26 7.18 8.31 -40.82 -0.34 6.84 16 50 H 0.08 1.10 8.14 -51.92 -0.42 0.68 16 51 H 0.08 1.03 8.14 -51.93 -0.42 0.61 16 52 H 0.07 1.05 7.81 -51.93 -0.41 0.64 16 53 H 0.09 1.22 7.01 -51.93 -0.36 0.86 16 LS Contribution 409.44 15.07 6.17 6.17 Total: -1.00 -31.72 409.44 -9.54 -41.26 By element: Atomic # 1 Polarization: 21.20 SS G_CDS: -8.78 Total: 12.41 kcal Atomic # 6 Polarization: -0.41 SS G_CDS: -1.11 Total: -1.52 kcal Atomic # 7 Polarization: -52.36 SS G_CDS: -0.09 Total: -52.45 kcal Atomic # 8 Polarization: -22.12 SS G_CDS: -0.70 Total: -22.82 kcal Atomic # 14 Polarization: 21.96 SS G_CDS: -5.02 Total: 16.94 kcal Total LS contribution 6.17 Total: 6.17 kcal Total: -31.72 -9.54 -41.26 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. ZINC001700308518.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 -95.145 kcal (2) G-P(sol) polarization free energy of solvation -31.720 kcal (3) G-ENP(sol) elect.-nuc.-pol. free energy of system -126.865 kcal (4) G-CDS(sol) cavity-dispersion-solvent structure free energy -9.542 kcal (5) G-P-CDS(sol) = G-P(sol) + G-CDS(sol) = (2) + (4) -41.262 kcal (6) G-S(sol) free energy of system = (1) + (5) -136.407 kcal 1SCF WAS SPECIFIED, SO GEOMETRY OPTIMIZATION WAS NOT USED Total computer time = 1.08 seconds