Wall clock time and date at job start Sat Apr 11 2020 06:10:28 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.53000 * 1 3 3 C 1.52998 * 109.49884 * 2 1 4 4 C 1.53049 * 109.49906 * 239.90374 * 2 1 3 5 5 C 1.53036 * 109.53858 * 178.68097 * 4 2 1 6 6 C 1.52769 * 109.38409 * 61.24125 * 5 4 2 7 7 H 1.08996 * 109.43074 * 65.42234 * 6 5 4 8 8 C 1.52991 * 109.50831 * 185.70564 * 6 5 4 9 9 N 1.46498 * 109.46914 * 179.84857 * 8 6 5 10 10 C 1.34785 * 120.00080 * 180.02562 * 9 8 6 11 11 O 1.21276 * 119.99847 * 359.97438 * 10 9 8 12 12 C 1.50706 * 119.99852 * 180.02562 * 10 9 8 13 13 C 1.53001 * 109.46650 * 179.97438 * 12 10 9 14 14 O 1.42973 * 109.47211 * 294.99758 * 13 12 10 15 15 C 1.43089 * 113.77038 * 181.20958 * 14 13 12 16 16 C 1.53261 * 108.94248 * 57.96878 * 15 14 13 17 17 C 1.53259 * 108.84484 * 305.47607 * 16 15 14 18 18 O 1.42965 * 109.46806 * 175.00163 * 13 12 10 19 19 N 1.46930 * 109.00604 * 305.58685 * 6 5 4 20 20 C 1.36930 * 120.63009 * 233.58328 * 19 6 5 21 21 H 1.08007 * 119.99974 * 29.59770 * 20 19 6 22 22 C 1.38816 * 120.00192 * 209.59743 * 20 19 6 23 23 N 1.31624 * 119.99990 * 197.18300 * 22 20 19 24 24 Si 1.86793 * 119.99857 * 17.17866 * 22 20 19 25 25 H 1.48496 * 109.47251 * 265.52201 * 24 22 20 26 26 H 1.48508 * 109.47238 * 25.52046 * 24 22 20 27 27 H 1.48500 * 109.47105 * 145.51843 * 24 22 20 28 28 C 1.46928 * 118.74065 * 53.30841 * 19 6 5 29 29 H 1.08998 * 109.46690 * 299.93333 * 1 2 3 30 30 H 1.09000 * 109.46844 * 59.93677 * 1 2 3 31 31 H 1.08998 * 109.46526 * 179.97438 * 1 2 3 32 32 H 1.09002 * 109.47628 * 180.02562 * 3 2 1 33 33 H 1.09004 * 109.46840 * 300.06248 * 3 2 1 34 34 H 1.08999 * 109.47525 * 60.06139 * 3 2 1 35 35 H 1.09002 * 109.45291 * 298.70329 * 4 2 1 36 36 H 1.08997 * 109.45950 * 58.65735 * 4 2 1 37 37 H 1.08996 * 109.49582 * 181.31449 * 5 4 2 38 38 H 1.08997 * 109.50193 * 301.41398 * 5 4 2 39 39 H 1.09001 * 109.47373 * 299.84452 * 8 6 5 40 40 H 1.08994 * 109.47569 * 59.84959 * 8 6 5 41 41 H 0.97000 * 120.00668 * 0.02562 * 9 8 6 42 42 H 1.08997 * 109.46982 * 299.99301 * 12 10 9 43 43 H 1.08998 * 109.46711 * 59.99897 * 12 10 9 44 44 H 1.08999 * 109.47157 * 54.99941 * 13 12 10 45 45 H 1.09001 * 109.55536 * 298.13425 * 15 14 13 46 46 H 1.09001 * 109.55807 * 177.80637 * 15 14 13 47 47 H 1.09004 * 109.57268 * 185.66672 * 16 15 14 48 48 H 1.08996 * 109.68124 * 65.34922 * 16 15 14 49 49 H 1.09000 * 109.56056 * 294.68524 * 17 16 15 50 50 H 1.09007 * 109.64670 * 174.41851 * 17 16 15 51 51 H 0.97001 * 119.99839 * 359.97438 * 23 22 20 52 52 H 1.09000 * 109.58575 * 66.42541 * 28 19 6 53 53 H 1.09000 * 109.58779 * 186.71532 * 28 19 6 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.0407 1.4422 0.0000 4 6 2.0409 -0.7235 -1.2482 5 6 3.5709 -0.6953 -1.2659 6 6 4.1056 -1.4216 -0.0328 7 1 3.8414 -2.4774 -0.0909 8 6 5.6271 -1.2732 0.0277 9 7 6.1386 -1.9669 1.2123 10 6 7.4610 -1.9701 1.4729 11 8 8.2275 -1.3992 0.7263 12 6 7.9872 -2.6842 2.6912 13 6 9.5088 -2.5357 2.7518 14 8 10.0980 -3.2076 1.6358 15 6 11.5276 -3.1495 1.6167 16 6 12.0677 -3.7608 2.9143 17 6 11.4185 -3.0474 4.1053 18 8 9.9958 -3.1120 3.9662 19 7 3.5047 -0.8288 1.1698 20 6 4.2936 -0.3778 2.1941 21 1 5.2862 -0.0053 1.9875 22 6 3.8179 -0.3995 3.4980 23 7 4.4075 0.3244 4.4258 24 14 2.3407 -1.4557 3.9357 25 1 1.1058 -0.6342 3.8630 26 1 2.2399 -2.5868 2.9787 27 1 2.5022 -1.9859 5.3134 28 6 2.0414 -0.7234 1.2498 29 1 -0.3632 0.5128 0.8906 30 1 -0.3633 0.5148 -0.8894 31 1 -0.3632 -1.0277 -0.0005 32 1 3.1307 1.4418 -0.0005 33 1 1.6767 1.9565 -0.8895 34 1 1.6785 1.9560 0.8905 35 1 1.6967 -1.7576 -1.2348 36 1 1.6582 -0.2252 -2.1389 37 1 3.9334 -1.1903 -2.1668 38 1 3.9160 0.3385 -1.2543 39 1 6.0706 -1.7085 -0.8679 40 1 5.8872 -0.2163 0.0850 41 1 5.5257 -2.4239 1.8093 42 1 7.5435 -2.2497 3.5870 43 1 7.7276 -3.7412 2.6332 44 1 9.7723 -1.4785 2.7200 45 1 11.8497 -2.1109 1.5410 46 1 11.9044 -3.7133 0.7633 47 1 13.1492 -3.6315 2.9564 48 1 11.8253 -4.8229 2.9484 49 1 11.7364 -2.0050 4.1238 50 1 11.7171 -3.5364 5.0326 51 1 5.1748 0.8726 4.1985 52 1 1.6053 -1.7212 1.2985 53 1 1.7638 -0.1594 2.1403 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 ZINC001099043865.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:10:28 Heat of formation + Delta-G solvation = -117.840833 kcal Electronic energy + Delta-G solvation = -32171.913811 eV Core-core repulsion = 28043.841590 eV Total energy + Delta-G solvation = -4128.072221 eV No. of doubly occupied orbitals = 66 Molecular weight (most abundant/longest-lived isotopes) = 340.220 amu Computer time = 0.88 seconds Orbital eigenvalues (eV) -40.05380 -39.30456 -38.43246 -35.59761 -35.18824 -33.23858 -33.12306 -32.26944 -30.68071 -30.09990 -27.49448 -26.78687 -26.00296 -25.28305 -24.18537 -22.74695 -22.42251 -21.16720 -19.95399 -19.53540 -18.51755 -17.56031 -17.05529 -16.54749 -16.37995 -15.65675 -15.41338 -15.36171 -14.81499 -14.75254 -14.67642 -14.33037 -14.12826 -13.61206 -13.46487 -13.29630 -13.22849 -12.89196 -12.82447 -12.58731 -12.16043 -12.09742 -11.94578 -11.60554 -11.58178 -11.45339 -11.40640 -11.11529 -11.04396 -10.99552 -10.87480 -10.58717 -10.46500 -10.18021 -10.01389 -9.90873 -9.69858 -9.61599 -9.43563 -9.26593 -9.02605 -8.61197 -8.21756 -6.99602 -6.28016 -3.35823 3.24875 3.35827 3.38565 3.48628 3.63787 3.73095 4.16207 4.55097 4.61364 4.62248 4.83654 4.94316 4.99835 5.15309 5.15371 5.25670 5.26647 5.31452 5.32316 5.45107 5.47699 5.58476 5.62242 5.64877 5.73502 5.79232 5.84839 5.93764 5.95615 5.98173 6.01533 6.05331 6.14424 6.21686 6.29524 6.31579 6.33109 6.37052 6.48050 6.51152 6.65481 6.72686 6.83928 6.87268 6.87810 7.08692 7.18869 7.24443 8.03673 8.25498 8.41296 8.89762 9.03850 9.56758 10.17741 11.11922 Molecular weight = 340.22amu Principal moments of inertia in cm(-1) A = 0.016629 B = 0.003228 C = 0.002855 Principal moments of inertia in units of 10**(-40)*gram-cm**2 A = 1683.379884 B = 8673.195227 C = 9804.615077 Net atomic charges, atomic populations, and dipole contributions calculated with CM2 Atom NO. Type Charge No. of electrons 1 C -0.139 4.139 2 C -0.086 4.086 3 C -0.126 4.126 4 C -0.096 4.096 5 C -0.132 4.132 6 C 0.163 3.837 7 H 0.034 0.966 8 C 0.140 3.860 9 N -0.693 5.693 10 C 0.518 3.482 11 O -0.557 6.557 12 C -0.128 4.128 13 C 0.268 3.732 14 O -0.377 6.377 15 C 0.058 3.942 16 C -0.193 4.193 17 C 0.058 3.942 18 O -0.394 6.394 19 N -0.599 5.599 20 C -0.341 4.341 21 H 0.077 0.923 22 C -0.011 4.011 23 N -0.919 5.919 24 Si 1.056 2.944 25 H -0.297 1.297 26 H -0.283 1.283 27 H -0.303 1.303 28 C 0.138 3.862 29 H 0.052 0.948 30 H 0.048 0.952 31 H 0.049 0.951 32 H 0.068 0.932 33 H 0.035 0.965 34 H 0.047 0.953 35 H 0.054 0.946 36 H 0.050 0.950 37 H 0.054 0.946 38 H 0.067 0.933 39 H 0.055 0.945 40 H 0.064 0.936 41 H 0.412 0.588 42 H 0.120 0.880 43 H 0.104 0.896 44 H 0.064 0.936 45 H 0.053 0.947 46 H 0.095 0.905 47 H 0.082 0.918 48 H 0.080 0.920 49 H 0.053 0.947 50 H 0.098 0.902 51 H 0.263 0.737 52 H 0.027 0.973 53 H 0.070 0.930 For the charges calculated by CM2: Dipole moment (debyes) X Y Z Total from point charges 5.668 -6.432 -6.269 10.620 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.197 4.197 2 C -0.088 4.088 3 C -0.184 4.184 4 C -0.133 4.133 5 C -0.171 4.171 6 C 0.056 3.944 7 H 0.052 0.948 8 C 0.015 3.985 9 N -0.343 5.343 10 C 0.304 3.696 11 O -0.436 6.436 12 C -0.167 4.167 13 C 0.172 3.828 14 O -0.297 6.297 15 C -0.019 4.019 16 C -0.231 4.231 17 C -0.018 4.018 18 O -0.315 6.315 19 N -0.324 5.324 20 C -0.469 4.469 21 H 0.094 0.906 22 C -0.213 4.213 23 N -0.558 5.558 24 Si 0.889 3.111 25 H -0.225 1.225 26 H -0.209 1.209 27 H -0.232 1.232 28 C 0.012 3.988 29 H 0.071 0.929 30 H 0.067 0.933 31 H 0.068 0.932 32 H 0.087 0.913 33 H 0.054 0.946 34 H 0.066 0.934 35 H 0.073 0.927 36 H 0.069 0.931 37 H 0.072 0.928 38 H 0.086 0.914 39 H 0.073 0.927 40 H 0.082 0.918 41 H 0.249 0.751 42 H 0.138 0.862 43 H 0.122 0.878 44 H 0.082 0.918 45 H 0.071 0.929 46 H 0.114 0.886 47 H 0.100 0.900 48 H 0.099 0.901 49 H 0.071 0.929 50 H 0.116 0.884 51 H 0.073 0.927 52 H 0.044 0.956 53 H 0.088 0.912 Dipole moment (debyes) X Y Z Total from point charges 5.007 -6.650 -5.544 10.001 hybrid contribution 0.920 0.590 -0.885 1.406 sum 5.927 -6.060 -6.429 10.639 Atomic orbital electron populations 1.21831 0.95696 1.01068 1.01125 1.21062 0.94003 0.96615 0.97115 1.21715 1.02440 0.94551 0.99662 1.21410 0.95466 0.99638 0.96819 1.21936 0.95879 1.01597 0.97714 1.20640 0.90675 0.95149 0.87950 0.94804 1.21199 0.93178 0.97459 0.86618 1.45871 1.06658 1.54255 1.27564 1.19973 0.82719 0.80763 0.86162 1.90649 1.56923 1.46074 1.49971 1.21439 0.92053 1.03709 0.99548 1.22904 0.91802 0.93057 0.75008 1.88182 1.18987 1.76824 1.45712 1.22575 0.81519 1.00310 0.97463 1.22809 1.01329 1.02120 0.96818 1.22610 0.80828 0.99839 0.98540 1.88088 1.21374 1.82046 1.39953 1.44381 1.00564 1.70360 1.17124 1.20337 1.00661 1.39886 0.86036 0.90573 1.26107 1.03606 1.00016 0.91541 1.70884 1.22522 1.26692 1.35734 0.83903 0.75561 0.76868 0.74795 1.22538 1.20870 1.23160 1.21457 0.84771 0.96778 0.95770 0.92884 0.93288 0.93155 0.91310 0.94556 0.93430 0.92742 0.93096 0.92755 0.91442 0.92700 0.91755 0.75142 0.86238 0.87810 0.91817 0.92920 0.88637 0.89984 0.90120 0.92868 0.88408 0.92693 0.95564 0.91152 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.14 -1.99 8.48 37.16 0.32 -1.67 16 2 C -0.09 -1.43 0.77 -154.38 -0.12 -1.55 16 3 C -0.13 -2.17 7.99 37.16 0.30 -1.88 16 4 C -0.10 -1.42 4.71 -26.69 -0.13 -1.55 16 5 C -0.13 -2.06 5.00 -26.83 -0.13 -2.19 16 6 C 0.16 2.93 2.99 -67.82 -0.20 2.72 16 7 H 0.03 0.61 8.14 -51.93 -0.42 0.19 16 8 C 0.14 2.60 4.42 -4.05 -0.02 2.58 16 9 N -0.69 -13.13 3.82 -60.30 -0.23 -13.36 16 10 C 0.52 10.01 7.71 -10.98 -0.08 9.93 16 11 O -0.56 -12.22 15.42 -8.85 -0.14 -12.36 16 12 C -0.13 -2.12 5.78 -27.88 -0.16 -2.29 16 13 C 0.27 4.36 2.80 37.16 0.10 4.47 16 14 O -0.38 -6.17 9.77 -64.76 -0.63 -6.80 16 15 C 0.06 0.70 6.96 37.29 0.26 0.96 16 16 C -0.19 -1.91 6.35 -26.46 -0.17 -2.07 16 17 C 0.06 0.65 6.96 37.29 0.26 0.91 16 18 O -0.39 -5.96 10.54 -59.03 -0.62 -6.58 16 19 N -0.60 -12.58 2.18 -165.02 -0.36 -12.94 16 20 C -0.34 -8.31 7.36 -15.06 -0.11 -8.42 16 21 H 0.08 2.05 5.97 -52.48 -0.31 1.73 16 22 C -0.01 -0.28 5.28 -17.43 -0.09 -0.38 16 23 N -0.92 -26.61 13.71 55.49 0.76 -25.85 16 24 Si 1.06 24.64 24.70 -169.99 -4.20 20.44 16 25 H -0.30 -6.79 6.74 56.52 0.38 -6.41 16 26 H -0.28 -6.16 6.32 56.52 0.36 -5.80 16 27 H -0.30 -7.36 7.11 56.52 0.40 -6.96 16 28 C 0.14 2.75 2.86 -3.71 -0.01 2.74 16 29 H 0.05 0.77 8.14 -51.93 -0.42 0.35 16 30 H 0.05 0.62 8.14 -51.93 -0.42 0.19 16 31 H 0.05 0.67 8.14 -51.93 -0.42 0.25 16 32 H 0.07 1.32 5.97 -51.93 -0.31 1.01 16 33 H 0.04 0.54 8.14 -51.93 -0.42 0.12 16 34 H 0.05 0.86 8.14 -51.93 -0.42 0.44 16 35 H 0.05 0.79 8.14 -51.93 -0.42 0.36 16 36 H 0.05 0.68 8.14 -51.93 -0.42 0.26 16 37 H 0.05 0.74 8.14 -51.93 -0.42 0.32 16 38 H 0.07 1.12 6.34 -51.93 -0.33 0.79 16 39 H 0.05 0.96 8.14 -51.93 -0.42 0.53 16 40 H 0.06 1.35 7.42 -51.93 -0.39 0.97 16 41 H 0.41 7.87 7.40 -40.82 -0.30 7.56 16 42 H 0.12 2.01 8.14 -51.93 -0.42 1.58 16 43 H 0.10 1.56 8.14 -51.93 -0.42 1.13 16 44 H 0.06 1.09 7.96 -51.93 -0.41 0.68 16 45 H 0.05 0.62 8.14 -51.93 -0.42 0.19 16 46 H 0.10 0.98 8.14 -51.93 -0.42 0.56 16 47 H 0.08 0.59 8.14 -51.93 -0.42 0.17 16 48 H 0.08 0.82 8.14 -51.93 -0.42 0.39 16 49 H 0.05 0.55 8.14 -51.93 -0.42 0.13 16 50 H 0.10 0.89 8.14 -51.92 -0.42 0.47 16 51 H 0.26 7.70 8.73 -40.82 -0.36 7.35 16 52 H 0.03 0.52 6.85 -51.93 -0.36 0.17 16 53 H 0.07 1.57 4.62 -51.93 -0.24 1.33 16 LS Contribution 394.54 15.07 5.95 5.95 Total: -1.00 -30.17 394.54 -8.94 -39.12 By element: Atomic # 1 Polarization: 19.54 SS G_CDS: -9.48 Total: 10.07 kcal Atomic # 6 Polarization: 2.30 SS G_CDS: 0.01 Total: 2.31 kcal Atomic # 7 Polarization: -52.32 SS G_CDS: 0.17 Total: -52.15 kcal Atomic # 8 Polarization: -24.34 SS G_CDS: -1.39 Total: -25.74 kcal Atomic # 14 Polarization: 24.64 SS G_CDS: -4.20 Total: 20.44 kcal Total LS contribution 5.95 Total: 5.95 kcal Total: -30.17 -8.94 -39.12 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. ZINC001099043865.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 -78.723 kcal (2) G-P(sol) polarization free energy of solvation -30.174 kcal (3) G-ENP(sol) elect.-nuc.-pol. free energy of system -108.898 kcal (4) G-CDS(sol) cavity-dispersion-solvent structure free energy -8.943 kcal (5) G-P-CDS(sol) = G-P(sol) + G-CDS(sol) = (2) + (4) -39.117 kcal (6) G-S(sol) free energy of system = (1) + (5) -117.841 kcal 1SCF WAS SPECIFIED, SO GEOMETRY OPTIMIZATION WAS NOT USED Total computer time = 0.88 seconds