Wall clock time and date at job start Sat Apr 11 2020 06:12:22 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 N 1.46504 * 1 3 3 C 1.46503 * 119.99848 * 2 1 4 4 C 1.52996 * 109.47248 * 89.99806 * 3 2 1 5 5 H 1.08998 * 109.46876 * 304.99609 * 4 3 2 6 6 O 1.42894 * 109.47458 * 64.99821 * 4 3 2 7 7 C 1.53006 * 109.47159 * 185.00104 * 4 3 2 8 8 N 1.46494 * 109.47167 * 179.97438 * 7 4 3 9 9 C 1.36934 * 120.00016 * 179.97438 * 8 7 4 10 10 H 1.08008 * 120.00173 * 0.02562 * 9 8 7 11 11 C 1.38810 * 120.00547 * 180.02562 * 9 8 7 12 12 N 1.31627 * 119.99737 * 0.02562 * 11 9 8 13 13 Si 1.86795 * 120.00547 * 180.02562 * 11 9 8 14 14 H 1.48505 * 109.47182 * 90.00037 * 13 11 9 15 15 H 1.48502 * 109.47369 * 209.99773 * 13 11 9 16 16 H 1.48504 * 109.47173 * 329.99873 * 13 11 9 17 17 C 1.34780 * 120.00081 * 180.26850 * 2 1 3 18 18 O 1.21280 * 119.99867 * 185.14742 * 17 2 1 19 19 C 1.50700 * 119.99875 * 5.14262 * 17 2 1 20 20 C 1.53005 * 112.84765 * 293.75723 * 19 17 2 21 21 H 1.09004 * 109.48068 * 64.81687 * 20 19 17 22 22 C 1.53092 * 109.48337 * 304.78565 * 20 19 17 23 23 C 1.53177 * 109.16141 * 182.41076 * 22 20 19 24 24 C 1.53173 * 109.03913 * 56.98035 * 23 22 20 25 25 C 1.53090 * 109.16386 * 303.01789 * 24 23 22 26 26 O 1.42900 * 109.47802 * 184.84150 * 20 19 17 27 27 C 1.53773 * 113.61386 * 162.45425 * 19 17 2 28 28 C 1.53787 * 87.07808 * 220.01789 * 27 19 17 29 29 C 1.53771 * 87.08223 * 25.43326 * 28 27 19 30 30 H 1.08992 * 109.47265 * 264.56868 * 1 2 3 31 31 H 1.09002 * 109.46942 * 24.57344 * 1 2 3 32 32 H 1.08995 * 109.46711 * 144.56354 * 1 2 3 33 33 H 1.09005 * 109.46796 * 209.99463 * 3 2 1 34 34 H 1.08994 * 109.46884 * 329.99497 * 3 2 1 35 35 H 0.96700 * 114.00419 * 60.00176 * 6 4 3 36 36 H 1.09002 * 109.46722 * 300.00028 * 7 4 3 37 37 H 1.08995 * 109.47292 * 59.99177 * 7 4 3 38 38 H 0.97000 * 119.99851 * 359.97438 * 8 7 4 39 39 H 0.96999 * 120.00078 * 180.02562 * 12 11 9 40 40 H 1.08995 * 109.52013 * 302.31330 * 22 20 19 41 41 H 1.08996 * 109.51983 * 62.50578 * 22 20 19 42 42 H 1.08998 * 109.53965 * 176.84976 * 23 22 20 43 43 H 1.08997 * 109.54256 * 297.11145 * 23 22 20 44 44 H 1.09002 * 109.52093 * 62.92072 * 24 23 22 45 45 H 1.08998 * 109.52338 * 183.11202 * 24 23 22 46 46 H 1.09000 * 109.48279 * 177.59132 * 25 24 23 47 47 H 1.08997 * 109.47995 * 297.61663 * 25 24 23 48 48 H 1.08998 * 113.61654 * 334.56375 * 27 19 17 49 49 H 1.09004 * 113.61749 * 105.47318 * 27 19 17 50 50 H 1.08997 * 113.61103 * 270.87752 * 28 27 19 51 51 H 1.09009 * 113.61222 * 139.98011 * 28 27 19 52 52 H 1.08998 * 113.61422 * 220.02017 * 29 28 27 53 53 H 1.09007 * 113.69677 * 89.12550 * 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 7 1.4650 0.0000 0.0000 3 6 2.1975 1.2688 0.0000 4 6 2.4525 1.7105 1.4425 5 1 1.5061 1.7551 1.9814 6 8 3.3208 0.7727 2.0814 7 6 3.1060 3.0940 1.4443 8 7 3.3495 3.5172 2.8254 9 6 3.9314 4.7309 3.0771 10 1 4.2103 5.3765 2.2573 11 6 4.1627 5.1318 4.3858 12 7 3.8232 4.3450 5.3848 13 14 4.9557 6.7879 4.7291 14 1 3.9004 7.8243 4.8607 15 1 5.7347 6.7116 5.9911 16 1 5.8626 7.1443 3.6084 17 6 2.1390 -1.1672 -0.0055 18 8 3.3481 -1.1697 0.0888 19 6 1.3900 -2.4691 -0.1286 20 6 0.5036 -2.7615 1.0838 21 1 -0.2844 -2.0111 1.1474 22 6 1.3515 -2.7209 2.3579 23 6 0.4700 -3.0736 3.5599 24 6 -0.1649 -4.4481 3.3277 25 6 -0.9542 -4.4220 2.0162 26 8 -0.0823 -4.0574 0.9443 27 6 2.2779 -3.6536 -0.5448 28 6 1.1135 -4.1103 -1.4395 29 6 0.6674 -2.6392 -1.4753 30 1 -0.3633 -0.0973 1.0230 31 1 -0.3633 0.9346 -0.4274 32 1 -0.3632 -0.8373 -0.5958 33 1 3.1500 1.1384 -0.5138 34 1 1.6085 2.0283 -0.5139 35 1 4.1836 0.6812 1.6544 36 1 4.0523 3.0494 0.9052 37 1 2.4436 3.8093 0.9568 38 1 3.0986 2.9376 3.5616 39 1 3.9845 4.6252 6.2993 40 1 2.1640 -3.4429 2.2778 41 1 1.7635 -1.7204 2.4894 42 1 1.0793 -3.1013 4.4632 43 1 -0.3136 -2.3239 3.6699 44 1 0.6171 -5.2050 3.2672 45 1 -0.8373 -4.6829 4.1529 46 1 -1.3729 -5.4099 1.8242 47 1 -1.7615 -3.6936 2.0922 48 1 3.1669 -3.3582 -1.1020 49 1 2.4978 -4.3423 0.2710 50 1 0.4070 -4.7688 -0.9343 51 1 1.4321 -4.4937 -2.4090 52 1 -0.4129 -2.5088 -1.4107 53 1 1.1152 -2.0652 -2.2866 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 ZINC001486289112.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:22 Heat of formation + Delta-G solvation = -97.978319 kcal Electronic energy + Delta-G solvation = -32314.089595 eV Core-core repulsion = 28186.878677 eV Total energy + Delta-G solvation = -4127.210917 eV No. of doubly occupied orbitals = 66 Molecular weight (most abundant/longest-lived isotopes) = 340.220 amu Computer time = 0.73 seconds Orbital eigenvalues (eV) -41.37369 -39.21303 -37.81620 -36.80328 -35.27603 -33.92752 -33.78905 -31.56978 -31.19589 -29.59375 -27.30938 -27.03603 -25.83968 -25.18740 -24.20367 -23.38616 -22.57538 -21.13319 -20.42350 -20.11607 -18.78273 -17.34785 -17.12102 -16.74029 -16.21965 -16.10435 -15.80290 -15.33985 -15.08258 -14.75505 -14.63383 -14.39709 -14.34147 -14.11202 -13.85896 -13.46796 -13.13413 -12.92681 -12.77055 -12.68719 -12.56582 -12.45817 -12.26970 -12.18002 -11.90566 -11.57613 -11.42874 -11.39229 -11.20437 -11.08683 -10.90544 -10.80191 -10.68733 -10.60085 -10.47195 -10.16286 -10.12732 -9.77982 -9.70024 -9.54835 -9.43892 -8.86859 -8.45002 -6.88704 -5.72137 -3.01357 2.57473 2.95927 3.45089 3.48673 3.50064 3.52164 4.13371 4.32398 4.45616 4.56943 4.60695 4.66952 4.69901 4.72891 4.83312 4.91755 4.93436 4.99014 5.03754 5.08422 5.19235 5.25839 5.36398 5.39848 5.42186 5.47924 5.56398 5.56754 5.59822 5.66157 5.67179 5.69050 5.73732 5.77270 5.79358 5.91876 5.96549 6.03031 6.05800 6.07019 6.21857 6.28326 6.38621 6.49949 6.56266 6.78484 6.98203 7.36765 7.71090 8.21132 8.28898 8.46435 9.57489 10.16179 10.52066 11.50236 Molecular weight = 340.22amu Principal moments of inertia in cm(-1) A = 0.016747 B = 0.002764 C = 0.002525 Principal moments of inertia in units of 10**(-40)*gram-cm**2 A = 1671.561636 B =10129.190625 C =11085.478236 Net atomic charges, atomic populations, and dipole contributions calculated with CM2 Atom NO. Type Charge No. of electrons 1 C 0.068 3.932 2 N -0.608 5.608 3 C 0.085 3.915 4 C 0.088 3.912 5 H 0.093 0.907 6 O -0.545 6.545 7 C 0.125 3.875 8 N -0.725 5.725 9 C -0.242 4.242 10 H 0.069 0.931 11 C -0.013 4.013 12 N -0.930 5.930 13 Si 0.907 3.093 14 H -0.290 1.290 15 H -0.291 1.291 16 H -0.282 1.282 17 C 0.539 3.461 18 O -0.523 6.523 19 C -0.099 4.099 20 C 0.119 3.881 21 H 0.068 0.932 22 C -0.151 4.151 23 C -0.121 4.121 24 C -0.156 4.156 25 C 0.060 3.940 26 O -0.370 6.370 27 C -0.104 4.104 28 C -0.134 4.134 29 C -0.122 4.122 30 H 0.056 0.944 31 H 0.074 0.926 32 H 0.075 0.925 33 H 0.085 0.915 34 H 0.084 0.916 35 H 0.374 0.626 36 H 0.020 0.980 37 H 0.023 0.977 38 H 0.389 0.611 39 H 0.255 0.745 40 H 0.073 0.927 41 H 0.084 0.916 42 H 0.070 0.930 43 H 0.062 0.938 44 H 0.070 0.930 45 H 0.071 0.929 46 H 0.092 0.908 47 H 0.048 0.952 48 H 0.086 0.914 49 H 0.074 0.926 50 H 0.101 0.899 51 H 0.062 0.938 52 H 0.075 0.925 53 H 0.081 0.919 For the charges calculated by CM2: Dipole moment (debyes) X Y Z Total from point charges -13.099 -19.974 -13.591 27.482 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.075 4.075 2 N -0.341 5.341 3 C -0.038 4.038 4 C 0.027 3.973 5 H 0.111 0.889 6 O -0.350 6.350 7 C -0.004 4.004 8 N -0.369 5.369 9 C -0.373 4.373 10 H 0.087 0.913 11 C -0.210 4.210 12 N -0.578 5.578 13 Si 0.738 3.262 14 H -0.218 1.218 15 H -0.217 1.217 16 H -0.208 1.208 17 C 0.326 3.674 18 O -0.398 6.398 19 C -0.102 4.102 20 C 0.061 3.939 21 H 0.085 0.915 22 C -0.189 4.189 23 C -0.159 4.159 24 C -0.194 4.194 25 C -0.016 4.016 26 O -0.289 6.289 27 C -0.142 4.142 28 C -0.171 4.171 29 C -0.159 4.159 30 H 0.074 0.926 31 H 0.093 0.907 32 H 0.094 0.906 33 H 0.103 0.897 34 H 0.102 0.898 35 H 0.221 0.779 36 H 0.038 0.962 37 H 0.041 0.959 38 H 0.223 0.777 39 H 0.065 0.935 40 H 0.092 0.908 41 H 0.103 0.897 42 H 0.089 0.911 43 H 0.081 0.919 44 H 0.089 0.911 45 H 0.089 0.911 46 H 0.111 0.889 47 H 0.067 0.933 48 H 0.104 0.896 49 H 0.092 0.908 50 H 0.119 0.881 51 H 0.081 0.919 52 H 0.093 0.907 53 H 0.100 0.900 Dipole moment (debyes) X Y Z Total from point charges -12.835 -19.990 -13.316 27.233 hybrid contribution 0.616 1.253 0.483 1.478 sum -12.219 -18.737 -12.833 25.788 Atomic orbital electron populations 1.22072 0.79994 1.03461 1.01948 1.47672 1.06609 1.05922 1.73854 1.22299 0.98142 0.84594 0.98772 1.22444 0.94443 0.89022 0.91396 0.88929 1.86529 1.29171 1.62389 1.56891 1.21238 0.96958 0.92492 0.89666 1.42235 1.73787 1.20105 1.00820 1.19221 1.32101 0.94177 0.91790 0.91325 1.24758 1.02042 0.98568 0.95599 1.69688 1.53125 1.36791 0.98203 0.86189 0.78865 0.83748 0.77367 1.21756 1.21748 1.20791 1.19979 0.87595 0.84190 0.75591 1.90669 1.14027 1.86778 1.48313 1.22050 0.98095 0.94152 0.95925 1.21730 0.93992 0.85581 0.92608 0.91451 1.22182 0.98777 1.03995 0.93905 1.21566 0.99331 0.96423 0.98550 1.22193 1.00650 0.98411 0.98141 1.22569 0.91721 1.00559 0.86786 1.87930 1.62926 1.33815 1.44238 1.22703 0.96767 0.94990 0.99698 1.22904 0.98603 0.96032 0.99587 1.22924 1.00088 0.98119 0.94792 0.92557 0.90729 0.90624 0.89679 0.89788 0.77872 0.96248 0.95932 0.77670 0.93542 0.90840 0.89741 0.91150 0.91877 0.91087 0.91064 0.88940 0.93321 0.89583 0.90776 0.88077 0.91883 0.90672 0.90038 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.07 0.60 8.98 59.85 0.54 1.13 16 2 N -0.61 -7.24 2.86 -181.57 -0.52 -7.76 16 3 C 0.09 1.12 5.14 -4.04 -0.02 1.10 16 4 C 0.09 1.48 3.30 -26.73 -0.09 1.39 16 5 H 0.09 1.62 8.14 -51.93 -0.42 1.20 16 6 O -0.54 -10.38 11.82 -48.02 -0.57 -10.95 16 7 C 0.13 2.53 5.92 -4.04 -0.02 2.51 16 8 N -0.72 -19.05 5.49 -59.91 -0.33 -19.38 16 9 C -0.24 -6.96 9.99 -15.06 -0.15 -7.11 16 10 H 0.07 1.91 7.83 -52.48 -0.41 1.49 16 11 C -0.01 -0.40 5.50 -17.43 -0.10 -0.50 16 12 N -0.93 -29.00 13.06 55.49 0.72 -28.28 16 13 Si 0.91 22.86 29.55 -169.99 -5.02 17.83 16 14 H -0.29 -7.23 7.11 56.52 0.40 -6.82 16 15 H -0.29 -7.39 7.11 56.52 0.40 -6.98 16 16 H -0.28 -6.91 7.11 56.52 0.40 -6.51 16 17 C 0.54 7.21 5.17 -10.98 -0.06 7.15 16 18 O -0.52 -8.77 14.46 -7.23 -0.10 -8.87 16 19 C -0.10 -1.06 0.30 -154.85 -0.05 -1.10 16 20 C 0.12 1.20 1.98 -26.68 -0.05 1.15 16 21 H 0.07 0.64 5.40 -51.93 -0.28 0.36 16 22 C -0.15 -1.78 5.04 -26.59 -0.13 -1.92 16 23 C -0.12 -1.22 6.03 -26.55 -0.16 -1.38 16 24 C -0.16 -1.26 6.14 -26.60 -0.16 -1.42 16 25 C 0.06 0.45 6.84 37.20 0.25 0.71 16 26 O -0.37 -3.45 7.29 -35.23 -0.26 -3.70 16 27 C -0.10 -1.12 6.52 -25.92 -0.17 -1.29 16 28 C -0.13 -1.17 7.53 -25.92 -0.20 -1.37 16 29 C -0.12 -1.05 5.65 -25.92 -0.15 -1.20 16 30 H 0.06 0.50 7.19 -51.93 -0.37 0.13 16 31 H 0.07 0.54 7.83 -51.93 -0.41 0.13 16 32 H 0.08 0.55 4.37 -51.93 -0.23 0.32 16 33 H 0.08 1.19 7.49 -51.93 -0.39 0.80 16 34 H 0.08 0.91 7.86 -51.93 -0.41 0.50 16 35 H 0.37 6.81 8.63 45.56 0.39 7.21 16 36 H 0.02 0.40 8.14 -51.93 -0.42 -0.02 16 37 H 0.02 0.45 8.14 -51.93 -0.42 0.02 16 38 H 0.39 11.03 7.54 -40.82 -0.31 10.72 16 39 H 0.25 7.66 8.31 -40.82 -0.34 7.32 16 40 H 0.07 0.93 7.44 -51.93 -0.39 0.55 16 41 H 0.08 1.24 7.38 -51.93 -0.38 0.86 16 42 H 0.07 0.74 8.14 -51.93 -0.42 0.32 16 43 H 0.06 0.62 8.14 -51.93 -0.42 0.20 16 44 H 0.07 0.60 8.14 -51.93 -0.42 0.18 16 45 H 0.07 0.48 8.14 -51.93 -0.42 0.06 16 46 H 0.09 0.55 8.14 -51.93 -0.42 0.12 16 47 H 0.05 0.33 8.14 -51.93 -0.42 -0.09 16 48 H 0.09 1.02 7.73 -51.93 -0.40 0.62 16 49 H 0.07 0.84 7.48 -51.93 -0.39 0.45 16 50 H 0.10 0.90 6.17 -51.93 -0.32 0.58 16 51 H 0.06 0.49 8.14 -51.92 -0.42 0.07 16 52 H 0.07 0.56 7.08 -51.93 -0.37 0.19 16 53 H 0.08 0.71 8.09 -51.93 -0.42 0.29 16 LS Contribution 401.14 15.07 6.05 6.05 Total: -1.00 -33.77 401.14 -9.18 -42.95 By element: Atomic # 1 Polarization: 22.70 SS G_CDS: -8.44 Total: 14.26 kcal Atomic # 6 Polarization: -1.44 SS G_CDS: -0.71 Total: -2.15 kcal Atomic # 7 Polarization: -55.29 SS G_CDS: -0.12 Total: -55.41 kcal Atomic # 8 Polarization: -22.60 SS G_CDS: -0.93 Total: -23.53 kcal Atomic # 14 Polarization: 22.86 SS G_CDS: -5.02 Total: 17.83 kcal Total LS contribution 6.05 Total: 6.05 kcal Total: -33.77 -9.18 -42.95 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. ZINC001486289112.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 -55.024 kcal (2) G-P(sol) polarization free energy of solvation -33.772 kcal (3) G-ENP(sol) elect.-nuc.-pol. free energy of system -88.796 kcal (4) G-CDS(sol) cavity-dispersion-solvent structure free energy -9.182 kcal (5) G-P-CDS(sol) = G-P(sol) + G-CDS(sol) = (2) + (4) -42.954 kcal (6) G-S(sol) free energy of system = (1) + (5) -97.978 kcal 1SCF WAS SPECIFIED, SO GEOMETRY OPTIMIZATION WAS NOT USED Total computer time = 0.73 seconds