Wall clock time and date at job start Sat Mar 6 2021 16:04:50 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= 0 * 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.45202 * 1 3 3 C 1.34233 * 117.00346 * 2 1 4 4 O 1.20828 * 119.99577 * 0.02562 * 3 2 1 5 5 C 1.50699 * 120.00172 * 179.97438 * 3 2 1 6 6 H 1.09004 * 112.84479 * 48.80723 * 5 3 2 7 7 C 1.53778 * 113.61086 * 179.97438 * 5 3 2 8 8 C 1.53787 * 87.08035 * 140.06661 * 7 5 3 9 9 H 1.08994 * 113.61447 * 220.02014 * 8 7 5 10 10 N 1.46496 * 113.61210 * 89.11438 * 8 7 5 11 11 C 1.34774 * 120.00099 * 155.00440 * 10 8 7 12 12 O 1.21560 * 119.99832 * 359.97438 * 11 10 8 13 13 C 1.47848 * 120.00646 * 180.02562 * 11 10 8 14 14 C 1.39549 * 120.14732 * 179.97438 * 13 11 10 15 15 C 1.37980 * 119.85208 * 179.97438 * 14 13 11 16 16 C 1.50707 * 119.92618 * 179.97438 * 15 14 13 17 17 O 1.42900 * 109.47056 * 90.00425 * 16 15 14 18 18 C 1.35896 * 116.99804 * 180.02562 * 17 16 15 19 19 C 1.38713 * 120.05333 * 180.02562 * 18 17 16 20 20 C 1.38162 * 119.94496 * 180.02562 * 19 18 17 21 21 C 1.50702 * 119.97274 * 180.02562 * 20 19 18 22 22 O 1.42895 * 109.46990 * 299.99430 * 21 20 19 23 23 Si 1.86301 * 109.46980 * 179.97438 * 21 20 19 24 24 C 1.38242 * 120.05451 * 359.97438 * 20 19 18 25 25 C 1.38283 * 120.11303 * 0.02562 * 24 20 19 26 26 C 1.38129 * 120.05612 * 359.67742 * 25 24 20 27 27 C 1.38339 * 120.14653 * 0.02562 * 15 14 13 28 28 C 1.38381 * 120.28786 * 359.97438 * 27 15 14 29 29 C 1.37935 * 120.14669 * 359.73918 * 28 27 15 30 30 C 1.53778 * 113.69101 * 277.53681 * 5 3 2 31 31 H 1.08996 * 109.47063 * 59.99778 * 1 2 3 32 32 H 1.08998 * 109.46979 * 179.97438 * 1 2 3 33 33 H 1.09001 * 109.46670 * 299.99840 * 1 2 3 34 34 H 1.08998 * 113.61879 * 25.52144 * 7 5 3 35 35 H 1.08999 * 113.61887 * 254.61251 * 7 5 3 36 36 H 0.97002 * 120.00490 * 335.00104 * 10 8 7 37 37 H 1.08002 * 120.07239 * 359.95349 * 14 13 11 38 38 H 1.09003 * 109.46651 * 209.99953 * 16 15 14 39 39 H 1.08997 * 109.46870 * 330.00212 * 16 15 14 40 40 H 1.08001 * 120.02926 * 359.97390 * 19 18 17 41 41 H 1.08995 * 109.47090 * 59.99844 * 21 20 19 42 42 H 0.96704 * 114.00079 * 180.02562 * 22 21 20 43 43 H 1.48500 * 110.00233 * 59.99671 * 23 21 20 44 44 H 1.48506 * 109.99875 * 181.31889 * 23 21 20 45 45 H 1.08000 * 119.94394 * 180.02562 * 24 20 19 46 46 H 1.08007 * 119.97570 * 179.97438 * 25 24 20 47 47 H 1.07998 * 120.02612 * 180.31809 * 26 25 24 48 48 H 1.07995 * 119.85153 * 179.97438 * 27 15 14 49 49 H 1.08002 * 119.92272 * 180.02562 * 28 27 15 50 50 H 1.08004 * 120.06902 * 180.27884 * 29 28 27 51 51 H 1.08993 * 113.61229 * 105.46294 * 30 5 3 52 52 H 1.09000 * 113.61263 * 334.55944 * 30 5 3 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.4520 0.0000 0.0000 3 6 2.0615 1.1960 0.0000 4 8 1.4034 2.2093 0.0005 5 6 3.5664 1.2748 0.0006 6 1 4.0239 0.6361 -0.7550 7 6 4.1078 2.7141 -0.0001 8 6 5.1821 2.2269 0.9865 9 1 5.3979 2.9384 1.7835 10 7 6.3908 1.7179 0.3337 11 6 7.5677 1.7310 0.9902 12 8 7.6264 2.1657 2.1239 13 6 8.7877 1.2179 0.3313 14 6 10.0080 1.2306 1.0081 15 6 11.1426 0.7498 0.3875 16 6 12.4615 0.7633 1.1168 17 8 12.6225 -0.4629 1.8327 18 6 13.7674 -0.6132 2.5492 19 6 13.9893 -1.7789 3.2676 20 6 15.1539 -1.9291 3.9957 21 6 15.3959 -3.1954 4.7760 22 8 14.3589 -3.3637 5.7445 23 14 17.0424 -3.0868 5.6409 24 6 16.0978 -0.9191 4.0081 25 6 15.8789 0.2439 3.2928 26 6 14.7190 0.3964 2.5584 27 6 11.0725 0.2562 -0.9029 28 6 9.8656 0.2408 -1.5796 29 6 8.7227 0.7131 -0.9686 30 6 4.1926 1.1230 1.3969 31 1 -0.3633 0.5138 -0.8899 32 1 -0.3633 -1.0277 -0.0005 33 1 -0.3633 0.5138 0.8900 34 1 3.4178 3.4404 0.4294 35 1 4.5100 3.0290 -0.9630 36 1 6.3440 1.3713 -0.5711 37 1 10.0642 1.6164 2.0154 38 1 13.2731 0.8733 0.3975 39 1 12.4808 1.5982 1.8171 40 1 13.2521 -2.5680 3.2586 41 1 15.3990 -4.0467 4.0953 42 1 14.4508 -4.1604 6.2849 43 1 17.0518 -1.9299 6.5719 44 1 17.3074 -4.3424 6.3883 45 1 17.0077 -1.0386 4.5775 46 1 16.6177 1.0317 3.3044 47 1 14.5509 1.3028 1.9957 48 1 11.9637 -0.1186 -1.3842 49 1 9.8175 -0.1462 -2.5868 50 1 7.7805 0.6963 -1.4963 51 1 4.6631 0.1537 1.5613 52 1 3.5257 1.4140 2.2084 There are 73 doubly filled levels RHF calculation, no. of doubly occupied orbitals= 73 No. of singly occupied orbitals= 1 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=0 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 s_22_12121848_4549630.mol2 52 1SCF WAS SPECIFIED, SO GEOMETRY OPTIMIZATION WAS NOT USED SCF FIELD WAS ACHIEVED Wall clock time and date at job start Sat Mar 6 2021 16:04:50 Heat of formation + Delta-G solvation = -97.057969 kcal Electronic energy + Delta-G solvation = -35032.133189 eV Core-core repulsion = 30123.017111 eV Total energy + Delta-G solvation = -4909.116078 eV No. of doubly occupied orbitals = 73 and no. of open levels = 1 Molecular weight (most abundant/longest-lived isotopes) = 398.165 amu Computer time = 1.38 seconds Orbital eigenvalues (eV) -42.08738 -41.44355 -41.06749 -40.06625 -39.10795 -38.63754 -37.36674 -37.32956 -34.72138 -33.30331 -32.53885 -32.23190 -30.99521 -29.45572 -27.82918 -27.13974 -26.96233 -26.56745 -25.07285 -23.58215 -23.18398 -22.54287 -22.09527 -21.43035 -20.50225 -19.12497 -18.91921 -18.60690 -18.50162 -17.89678 -17.44377 -17.36828 -17.24713 -17.13899 -17.06805 -16.80581 -16.10964 -16.06242 -15.81196 -15.61538 -15.52449 -15.22633 -14.99799 -14.75831 -14.52424 -14.46561 -14.27787 -14.05048 -13.90224 -13.71011 -13.57529 -13.52385 -13.43117 -13.22381 -13.18557 -12.88566 -12.80966 -12.51764 -12.18659 -12.02299 -11.94928 -11.91096 -11.84499 -11.79540 -11.64291 -11.33833 -11.02680 -10.54493 -10.29006 -10.10850 -10.04107 -9.82424 -9.17218 -4.95182 -0.31915 0.10140 0.27612 0.48741 1.13769 1.27726 1.60869 1.66858 1.69096 1.90351 1.97007 2.26453 2.69389 2.82742 3.12326 3.21696 3.33578 3.44677 3.47357 3.59205 3.65401 3.66552 3.80582 3.83436 3.91622 3.93863 3.99235 4.04433 4.13087 4.14891 4.18378 4.22721 4.29140 4.30470 4.36521 4.38272 4.51593 4.56398 4.57725 4.58900 4.65120 4.68068 4.69704 4.76040 4.80750 4.87793 5.03924 5.10609 5.18085 5.31215 5.31988 5.37159 5.39061 5.43416 5.53290 5.67541 5.73260 6.32797 6.42694 6.44706 6.71473 7.09671 Net atomic charges, atomic populations, and dipole contributions calculated with CM2 Atom NO. Type Charge No. of electrons 1 C 0.034 3.966 2 O -0.358 6.358 3 C 0.470 3.530 4 O -0.503 6.503 5 C -0.113 4.113 6 H 0.125 0.875 7 C -0.124 4.124 8 C 0.122 3.878 9 H 0.116 0.884 10 N -0.711 5.711 11 C 0.565 3.435 12 O -0.527 6.527 13 C -0.140 4.140 14 C -0.033 4.033 15 C -0.104 4.104 16 C 0.101 3.899 17 O -0.308 6.308 18 C 0.131 3.869 19 C -0.123 4.123 20 C -0.034 4.034 21 C -0.078 4.078 22 O -0.535 6.535 23 Si 0.686 3.314 24 C -0.153 4.153 25 C -0.072 4.072 26 C -0.205 4.205 27 C -0.066 4.066 28 C -0.128 4.128 29 C -0.075 4.075 30 C -0.132 4.132 31 H 0.064 0.936 32 H 0.103 0.897 33 H 0.064 0.936 34 H 0.114 0.886 35 H 0.086 0.914 36 H 0.401 0.599 37 H 0.139 0.861 38 H 0.076 0.924 39 H 0.076 0.924 40 H 0.132 0.868 41 H 0.079 0.921 42 H 0.380 0.620 43 H -0.255 1.255 44 H -0.261 1.261 45 H 0.129 0.871 46 H 0.132 0.868 47 H 0.131 0.869 48 H 0.132 0.868 49 H 0.135 0.865 50 H 0.128 0.872 51 H 0.088 0.912 52 H 0.100 0.900 For the charges calculated by CM2: Dipole moment (debyes) X Y Z Total from point charges 3.120 -2.520 -5.945 7.172 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.060 4.060 2 O -0.274 6.274 3 C 0.312 3.688 4 O -0.391 6.391 5 C -0.132 4.132 6 H 0.142 0.858 7 C -0.162 4.162 8 C 0.017 3.983 9 H 0.134 0.866 10 N -0.363 5.363 11 C 0.352 3.648 12 O -0.404 6.404 13 C -0.143 4.143 14 C -0.052 4.052 15 C -0.104 4.104 16 C 0.026 3.974 17 O -0.222 6.222 18 C 0.086 3.914 19 C -0.142 4.142 20 C -0.035 4.035 21 C -0.187 4.187 22 O -0.342 6.342 23 Si 0.587 3.413 24 C -0.170 4.170 25 C -0.090 4.090 26 C -0.224 4.224 27 C -0.084 4.084 28 C -0.146 4.146 29 C -0.094 4.094 30 C -0.171 4.171 31 H 0.083 0.917 32 H 0.122 0.878 33 H 0.083 0.917 34 H 0.132 0.868 35 H 0.105 0.895 36 H 0.237 0.763 37 H 0.157 0.843 38 H 0.095 0.905 39 H 0.094 0.906 40 H 0.150 0.850 41 H 0.096 0.904 42 H 0.228 0.772 43 H -0.179 1.179 44 H -0.187 1.187 45 H 0.146 0.854 46 H 0.149 0.851 47 H 0.148 0.852 48 H 0.150 0.850 49 H 0.153 0.847 50 H 0.146 0.854 51 H 0.107 0.893 52 H 0.118 0.882 Dipole moment (debyes) X Y Z Total from point charges 2.961 -1.480 -4.641 5.701 hybrid contribution -1.216 -0.712 -0.056 1.411 sum 1.745 -2.192 -4.697 5.469 Atomic orbital electron populations 1.23326 0.75924 1.04020 1.02732 1.86521 1.23752 1.33348 1.83789 1.23555 0.91650 0.81440 0.72126 1.90690 1.66958 1.36420 1.44992 1.22657 0.87695 1.01028 1.01845 0.85760 1.23371 0.96965 0.94828 1.01050 1.21815 0.83408 0.96437 0.96676 0.86586 1.45571 1.07388 1.64779 1.18577 1.18249 0.82954 0.77945 0.85612 1.90818 1.85829 1.45837 1.17901 1.19830 0.93874 1.05012 0.95551 1.21192 0.90201 0.94171 0.99588 1.19429 0.94130 1.02036 0.94777 1.21597 0.94797 0.86123 0.94835 1.86159 1.36785 1.38590 1.60627 1.18974 0.85675 0.94443 0.92356 1.20222 0.98263 0.96954 0.98757 1.19021 0.94152 0.94376 0.95922 1.24846 1.01748 0.97695 0.94456 1.86603 1.47573 1.48806 1.51199 0.95017 0.80632 0.80461 0.85209 1.20854 0.99578 0.94018 1.02594 1.20885 0.96441 0.96964 0.94745 1.21313 0.96073 0.99809 1.05204 1.21342 0.98527 0.94885 0.93639 1.21136 0.92156 1.00587 1.00676 1.21364 0.98964 0.95609 0.93436 1.23413 0.97813 1.01259 0.94602 0.91735 0.87848 0.91717 0.86808 0.89524 0.76272 0.84270 0.90543 0.90621 0.84986 0.90350 0.77174 1.17912 1.18651 0.85371 0.85059 0.85158 0.85008 0.84706 0.85385 0.89342 0.88188 Number of geometries 1 Number of calculations of the screened coulomb radii 1 The total number of SCF iterations 12. 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.02 10.56 101.05 1.07 1.05 16 2 O -0.36 -0.46 10.57 -39.26 -0.41 -0.88 16 3 C 0.47 0.71 7.49 36.01 0.27 0.98 16 4 O -0.50 -2.19 15.65 -18.75 -0.29 -2.48 16 5 C -0.11 0.18 4.13 -90.96 -0.38 -0.20 16 6 H 0.12 -0.46 8.14 -51.93 -0.42 -0.88 16 7 C -0.12 0.22 7.15 -25.92 -0.19 0.03 16 8 C 0.12 -0.11 4.33 -67.09 -0.29 -0.40 16 9 H 0.12 0.08 7.46 -51.93 -0.39 -0.31 16 10 N -0.71 0.90 5.21 -53.94 -0.28 0.62 16 11 C 0.56 1.07 7.77 -12.33 -0.10 0.98 16 12 O -0.53 -3.07 15.95 5.31 0.08 -2.98 16 13 C -0.14 -0.10 5.87 -104.97 -0.62 -0.71 16 14 C -0.03 -0.05 9.24 -39.19 -0.36 -0.41 16 15 C -0.10 -0.02 5.50 -104.68 -0.58 -0.60 16 16 C 0.10 -0.02 4.86 36.01 0.18 0.15 16 17 O -0.31 -0.79 10.12 -19.41 -0.20 -0.98 16 18 C 0.13 0.33 6.69 -39.22 -0.26 0.07 16 19 C -0.12 -0.37 9.55 -39.43 -0.38 -0.75 16 20 C -0.03 -0.07 5.42 -104.65 -0.57 -0.64 16 21 C -0.08 -0.04 2.81 36.01 0.10 0.06 16 22 O -0.53 -1.05 12.74 -35.23 -0.45 -1.50 16 23 Si 0.69 -0.41 29.79 -169.99 -5.06 -5.47 16 24 C -0.15 -0.14 8.67 -39.56 -0.34 -0.48 16 25 C -0.07 -0.02 10.04 -39.60 -0.40 -0.42 16 26 C -0.21 -0.17 9.04 -39.43 -0.36 -0.53 16 27 C -0.07 0.08 9.69 -39.49 -0.38 -0.31 16 28 C -0.13 0.26 10.03 -39.64 -0.40 -0.14 16 29 C -0.08 0.14 9.54 -39.19 -0.37 -0.24 16 30 C -0.13 0.16 7.43 -25.90 -0.19 -0.03 16 31 H 0.06 -0.04 8.13 -51.93 -0.42 -0.46 16 32 H 0.10 -0.30 8.14 -51.93 -0.42 -0.72 16 33 H 0.06 -0.03 8.13 -51.93 -0.42 -0.45 16 34 H 0.11 -0.07 7.44 -51.93 -0.39 -0.45 16 35 H 0.09 -0.25 8.14 -51.93 -0.42 -0.68 16 36 H 0.40 -1.64 6.74 -40.82 -0.27 -1.91 16 37 H 0.14 0.30 7.65 -52.49 -0.40 -0.10 16 38 H 0.08 -0.16 7.59 -51.93 -0.39 -0.56 16 39 H 0.08 -0.07 7.61 -51.93 -0.40 -0.47 16 40 H 0.13 0.41 8.06 -52.49 -0.42 -0.02 16 41 H 0.08 -0.09 7.96 -51.93 -0.41 -0.51 16 42 H 0.38 -0.83 8.74 45.56 0.40 -0.43 16 43 H -0.25 -0.67 7.11 56.52 0.40 -0.27 16 44 H -0.26 -0.34 7.11 56.52 0.40 0.06 16 45 H 0.13 -0.13 4.87 -52.49 -0.26 -0.39 16 46 H 0.13 -0.28 8.06 -52.48 -0.42 -0.70 16 47 H 0.13 -0.17 6.30 -52.49 -0.33 -0.50 16 48 H 0.13 -0.38 8.06 -52.49 -0.42 -0.80 16 49 H 0.14 -0.57 8.06 -52.49 -0.42 -0.99 16 50 H 0.13 -0.62 6.41 -52.48 -0.34 -0.95 16 51 H 0.09 -0.10 8.14 -51.93 -0.42 -0.52 16 52 H 0.10 -0.10 8.10 -51.93 -0.42 -0.52 16 LS Contribution 438.00 15.07 6.60 6.60 Total: 0.00 -11.57 438.00 -11.58 -23.14 By element: Atomic # 1 Polarization: -6.51 SS G_CDS: -7.02 Total: -13.53 kcal Atomic # 6 Polarization: 2.00 SS G_CDS: -4.54 Total: -2.54 kcal Atomic # 7 Polarization: 0.90 SS G_CDS: -0.28 Total: 0.62 kcal Atomic # 8 Polarization: -7.55 SS G_CDS: -1.27 Total: -8.82 kcal Atomic # 14 Polarization: -0.41 SS G_CDS: -5.06 Total: -5.47 kcal Total LS contribution 6.60 Total: 6.60 kcal Total: -11.57 -11.58 -23.14 kcal The number of atoms in the molecule is 52 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. s_22_12121848_4549630.mol2 52 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 -73.914 kcal (2) G-P(sol) polarization free energy of solvation -11.569 kcal (3) G-ENP(sol) elect.-nuc.-pol. free energy of system -85.483 kcal (4) G-CDS(sol) cavity-dispersion-solvent structure free energy -11.575 kcal (5) G-P-CDS(sol) = G-P(sol) + G-CDS(sol) = (2) + (4) -23.144 kcal (6) G-S(sol) free energy of system = (1) + (5) -97.058 kcal 1SCF WAS SPECIFIED, SO GEOMETRY OPTIMIZATION WAS NOT USED Total computer time = 1.38 seconds