Wall clock time and date at job start Sat Apr 11 2020 06:10:14 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.30996 * 1 3 3 C 1.50700 * 119.99831 * 2 1 4 4 C 1.50698 * 120.00288 * 180.02562 * 2 1 3 5 5 H 1.09001 * 109.47271 * 4.99244 * 4 2 1 6 6 C 1.53005 * 109.46791 * 244.99926 * 4 2 1 7 7 C 1.50703 * 109.46957 * 174.01714 * 6 4 2 8 8 O 1.21282 * 119.99742 * 4.57627 * 7 6 4 9 9 N 1.34771 * 120.00073 * 184.57745 * 7 6 4 10 10 C 1.46496 * 120.00155 * 185.10739 * 9 7 6 11 11 C 1.53006 * 109.54642 * 300.77082 * 10 9 7 12 12 N 1.46505 * 109.46498 * 305.41564 * 11 10 9 13 13 C 1.36938 * 119.99609 * 185.01671 * 12 11 10 14 14 H 1.08008 * 119.99843 * 4.98270 * 13 12 11 15 15 C 1.38811 * 119.99979 * 184.97952 * 13 12 11 16 16 N 1.31625 * 119.99861 * 184.81246 * 15 13 12 17 17 Si 1.86798 * 120.00167 * 4.81064 * 15 13 12 18 18 H 1.48502 * 109.47086 * 84.86864 * 17 15 13 19 19 H 1.48496 * 109.47299 * 204.87424 * 17 15 13 20 20 H 1.48502 * 109.46878 * 324.87144 * 17 15 13 21 21 C 1.53180 * 109.54644 * 180.59351 * 10 9 7 22 22 C 1.53287 * 109.05935 * 297.03157 * 21 10 9 23 23 O 1.42902 * 109.45961 * 302.44084 * 22 21 10 24 24 C 1.42904 * 114.09844 * 61.16366 * 23 22 21 25 25 C 1.53094 * 109.41271 * 298.84056 * 24 23 22 26 26 O 1.42898 * 109.47222 * 124.99997 * 4 2 1 27 27 C 1.42899 * 113.99684 * 266.22237 * 26 4 2 28 28 C 1.53000 * 109.47124 * 184.93517 * 27 26 4 29 29 H 1.07994 * 120.00064 * 359.97438 * 1 2 3 30 30 H 1.08005 * 119.99789 * 179.73006 * 1 2 3 31 31 H 1.08997 * 109.47570 * 90.00301 * 3 2 1 32 32 H 1.09003 * 109.47141 * 210.00149 * 3 2 1 33 33 H 1.09002 * 109.47235 * 329.99922 * 3 2 1 34 34 H 1.09000 * 109.47376 * 294.01071 * 6 4 2 35 35 H 1.08993 * 109.47431 * 54.01706 * 6 4 2 36 36 H 0.97005 * 120.00033 * 5.11292 * 9 7 6 37 37 H 1.08999 * 109.47133 * 65.41122 * 11 10 9 38 38 H 1.08995 * 109.47167 * 185.41891 * 11 10 9 39 39 H 0.96999 * 120.00368 * 5.01502 * 12 11 10 40 40 H 0.97002 * 119.99716 * 180.02562 * 16 15 13 41 41 H 1.08996 * 109.52350 * 57.01557 * 21 10 9 42 42 H 1.09001 * 109.52487 * 177.21124 * 21 10 9 43 43 H 1.09003 * 109.45656 * 182.45120 * 22 21 10 44 44 H 1.08994 * 109.46837 * 62.44890 * 22 21 10 45 45 H 1.08996 * 109.48545 * 178.86094 * 24 23 22 46 46 H 1.08997 * 109.49272 * 58.82347 * 24 23 22 47 47 H 1.08999 * 109.52096 * 297.69995 * 25 24 23 48 48 H 1.09004 * 109.57540 * 177.54410 * 25 24 23 49 49 H 1.08998 * 109.47750 * 304.94101 * 27 26 4 50 50 H 1.09004 * 109.46938 * 64.93551 * 27 26 4 51 51 H 1.09005 * 109.46832 * 175.02865 * 28 27 26 52 52 H 1.09000 * 109.47374 * 295.02066 * 28 27 26 53 53 H 1.08992 * 109.47281 * 55.02335 * 28 27 26 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.3100 0.0000 0.0000 3 6 2.0634 1.3051 0.0000 4 6 2.0635 -1.3050 -0.0006 5 1 1.3586 -2.1316 -0.0904 6 6 2.8465 -1.4424 1.3068 7 6 3.4811 -2.8077 1.3711 8 8 3.3926 -3.5627 0.4261 9 7 4.1471 -3.1898 2.4788 10 6 4.8585 -4.4702 2.5007 11 6 3.8715 -5.6089 2.2355 12 7 2.7665 -5.5303 3.1942 13 6 1.8093 -6.5093 3.2174 14 1 1.9114 -7.3793 2.5856 15 6 0.7079 -6.3835 4.0528 16 7 -0.1630 -7.3659 4.1460 17 14 0.4618 -4.8199 5.0449 18 1 -0.2318 -3.8073 4.2090 19 1 -0.3606 -5.1191 6.2445 20 1 1.7826 -4.2903 5.4698 21 6 5.5140 -4.6689 3.8708 22 6 6.5568 -3.5673 4.0914 23 8 7.5192 -3.6046 3.0357 24 6 6.9690 -3.3786 1.7363 25 6 5.9445 -4.4715 1.4204 26 8 2.9709 -1.3307 -1.1042 27 6 2.4422 -1.9528 -2.2770 28 6 3.5328 -2.0341 -3.3470 29 1 -0.5400 0.9353 0.0004 30 1 -0.5400 -0.9354 -0.0044 31 1 2.2452 1.6198 1.0276 32 1 3.0159 1.1748 -0.5139 33 1 1.4743 2.0648 -0.5138 34 1 3.6231 -0.6786 1.3472 35 1 2.1688 -1.3166 2.1511 36 1 4.1552 -2.6169 3.2616 37 1 3.4801 -5.5217 1.2220 38 1 4.3819 -6.5654 2.3470 39 1 2.7129 -4.7848 3.8126 40 1 -0.9329 -7.2779 4.7295 41 1 4.7534 -4.6146 4.6497 42 1 6.0014 -5.6433 3.9036 43 1 7.0581 -3.7277 5.0460 44 1 6.0629 -2.5958 4.0990 45 1 7.7670 -3.4030 0.9943 46 1 6.4803 -2.4046 1.7127 47 1 6.4406 -5.4419 1.4037 48 1 5.4914 -4.2773 0.4482 49 1 2.0981 -2.9575 -2.0317 50 1 1.6058 -1.3645 -2.6546 51 1 3.1565 -2.5867 -4.2080 52 1 3.8143 -1.0275 -3.6561 53 1 4.4043 -2.5460 -2.9392 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 ZINC001502626342.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:14 Heat of formation + Delta-G solvation = -100.078232 kcal Electronic energy + Delta-G solvation = -32619.935110 eV Core-core repulsion = 28492.633133 eV Total energy + Delta-G solvation = -4127.301976 eV No. of doubly occupied orbitals = 66 Molecular weight (most abundant/longest-lived isotopes) = 340.220 amu Computer time = 0.72 seconds Orbital eigenvalues (eV) -40.61875 -39.07766 -37.87092 -35.81705 -35.21733 -33.73584 -32.18280 -31.36884 -31.12607 -30.76621 -29.71339 -27.49553 -26.85751 -24.29956 -23.71710 -23.20201 -22.68770 -21.31125 -20.65415 -19.17615 -18.84561 -17.71865 -17.08635 -16.85661 -16.15378 -15.98273 -15.62308 -15.45023 -15.00069 -14.71927 -14.54572 -14.47390 -14.07878 -13.85092 -13.63110 -13.44999 -13.29108 -13.16366 -12.81331 -12.67886 -12.48855 -12.39629 -12.23199 -12.03453 -11.80149 -11.76892 -11.69178 -11.59981 -11.24797 -11.15792 -11.04428 -10.90067 -10.80290 -10.59543 -10.32129 -10.20785 -9.63341 -9.57376 -9.47539 -9.39835 -9.22869 -8.80594 -8.62969 -7.14862 -5.74982 -3.08302 2.00736 2.98145 3.28367 3.33303 3.35752 3.41142 3.76359 3.89099 4.55859 4.58065 4.71428 4.81968 4.91418 4.91737 4.99510 5.05614 5.08386 5.15324 5.20001 5.28908 5.33153 5.39061 5.42212 5.45136 5.50883 5.54863 5.56953 5.67601 5.68320 5.73828 5.75456 5.79394 5.83905 5.85884 5.95096 5.99108 6.01515 6.11232 6.14067 6.22618 6.31598 6.42595 6.46009 6.62439 6.82318 7.01281 7.19510 7.31823 7.77186 8.03851 8.44141 8.56660 9.24428 9.96223 10.17702 11.44016 Molecular weight = 340.22amu Principal moments of inertia in cm(-1) A = 0.007758 B = 0.005051 C = 0.003759 Principal moments of inertia in units of 10**(-40)*gram-cm**2 A = 3608.358467 B = 5542.206344 C = 7446.107094 Net atomic charges, atomic populations, and dipole contributions calculated with CM2 Atom NO. Type Charge No. of electrons 1 C -0.184 4.184 2 C -0.147 4.147 3 C -0.114 4.114 4 C 0.144 3.856 5 H 0.076 0.924 6 C -0.139 4.139 7 C 0.522 3.478 8 O -0.514 6.514 9 N -0.724 5.724 10 C 0.172 3.828 11 C 0.148 3.852 12 N -0.748 5.748 13 C -0.233 4.233 14 H 0.088 0.912 15 C 0.002 3.998 16 N -0.935 5.935 17 Si 0.884 3.116 18 H -0.285 1.285 19 H -0.292 1.292 20 H -0.283 1.283 21 C -0.155 4.155 22 C 0.061 3.939 23 O -0.385 6.385 24 C 0.067 3.933 25 C -0.161 4.161 26 O -0.367 6.367 27 C 0.064 3.936 28 C -0.146 4.146 29 H 0.100 0.900 30 H 0.100 0.900 31 H 0.064 0.936 32 H 0.070 0.930 33 H 0.059 0.941 34 H 0.099 0.901 35 H 0.100 0.900 36 H 0.392 0.608 37 H 0.054 0.946 38 H 0.031 0.969 39 H 0.357 0.643 40 H 0.258 0.742 41 H 0.084 0.916 42 H 0.089 0.911 43 H 0.093 0.907 44 H 0.045 0.955 45 H 0.089 0.911 46 H 0.047 0.953 47 H 0.078 0.922 48 H 0.093 0.907 49 H 0.056 0.944 50 H 0.039 0.961 51 H 0.066 0.934 52 H 0.055 0.945 53 H 0.064 0.936 For the charges calculated by CM2: Dipole moment (debyes) X Y Z Total from point charges 10.507 17.456 -7.436 21.689 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.221 4.221 2 C -0.148 4.148 3 C -0.170 4.170 4 C 0.086 3.914 5 H 0.094 0.906 6 C -0.179 4.179 7 C 0.308 3.692 8 O -0.389 6.389 9 N -0.377 5.377 10 C 0.086 3.914 11 C 0.020 3.980 12 N -0.394 5.394 13 C -0.363 4.363 14 H 0.106 0.894 15 C -0.193 4.193 16 N -0.584 5.584 17 Si 0.715 3.285 18 H -0.212 1.212 19 H -0.219 1.219 20 H -0.209 1.209 21 C -0.194 4.194 22 C -0.015 4.015 23 O -0.304 6.304 24 C -0.010 4.010 25 C -0.200 4.200 26 O -0.286 6.286 27 C -0.013 4.013 28 C -0.204 4.204 29 H 0.118 0.882 30 H 0.118 0.882 31 H 0.083 0.917 32 H 0.089 0.911 33 H 0.078 0.922 34 H 0.117 0.883 35 H 0.119 0.881 36 H 0.226 0.774 37 H 0.072 0.928 38 H 0.049 0.951 39 H 0.190 0.810 40 H 0.068 0.932 41 H 0.103 0.897 42 H 0.107 0.893 43 H 0.111 0.889 44 H 0.064 0.936 45 H 0.108 0.892 46 H 0.066 0.934 47 H 0.097 0.903 48 H 0.112 0.888 49 H 0.074 0.926 50 H 0.057 0.943 51 H 0.085 0.915 52 H 0.074 0.926 53 H 0.083 0.917 Dipole moment (debyes) X Y Z Total from point charges 10.716 16.727 -8.175 21.481 hybrid contribution -1.086 1.256 1.295 2.105 sum 9.630 17.982 -6.880 21.527 Atomic orbital electron populations 1.23820 0.95175 1.02229 1.00891 1.22255 0.95465 0.96599 1.00446 1.20750 1.00713 0.93789 1.01789 1.20691 0.90845 0.93800 0.86109 0.90594 1.21821 1.01366 0.96171 0.98536 1.20248 0.78023 0.88315 0.82570 1.90651 1.58122 1.53544 1.36585 1.46046 1.54077 1.21596 1.16004 1.20625 0.91638 0.84336 0.94789 1.21084 0.88263 0.93529 0.95128 1.42545 1.27913 1.28098 1.40841 1.19066 0.98205 0.98865 1.20206 0.89375 1.24664 0.97213 0.97417 1.00026 1.69704 1.14733 1.31512 1.42442 0.86585 0.77200 0.84712 0.79997 1.21232 1.21869 1.20948 1.22485 0.98857 1.00891 0.97118 1.22621 0.90775 0.97303 0.90848 1.87903 1.34207 1.92962 1.15295 1.22507 0.96654 0.97283 0.84546 1.22498 0.96264 1.00375 1.00816 1.87845 1.40118 1.79711 1.20898 1.22413 0.96106 0.96988 0.85794 1.21691 0.99411 1.01960 0.97290 0.88159 0.88210 0.91685 0.91110 0.92242 0.88284 0.88146 0.77391 0.92780 0.95085 0.80997 0.93180 0.89726 0.89264 0.88864 0.93635 0.89249 0.93428 0.90303 0.88824 0.92571 0.94272 0.91494 0.92576 0.91698 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.18 -1.97 13.91 29.02 0.40 -1.57 16 2 C -0.15 -1.64 5.08 -101.07 -0.51 -2.15 16 3 C -0.11 -1.04 9.75 36.01 0.35 -0.69 16 4 C 0.14 1.90 1.88 -27.88 -0.05 1.84 16 5 H 0.08 1.16 7.17 -51.93 -0.37 0.79 16 6 C -0.14 -1.71 5.24 -27.88 -0.15 -1.86 16 7 C 0.52 8.01 7.13 -10.99 -0.08 7.93 16 8 O -0.51 -9.50 11.72 -9.76 -0.11 -9.61 16 9 N -0.72 -10.19 3.59 -47.43 -0.17 -10.36 16 10 C 0.17 2.63 0.32 -131.63 -0.04 2.58 16 11 C 0.15 2.98 4.84 -4.04 -0.02 2.96 16 12 N -0.75 -17.67 4.94 -59.90 -0.30 -17.97 16 13 C -0.23 -6.69 10.47 -15.06 -0.16 -6.85 16 14 H 0.09 2.75 8.06 -52.48 -0.42 2.33 16 15 C 0.00 0.06 5.50 -17.43 -0.10 -0.04 16 16 N -0.94 -29.67 13.97 55.49 0.78 -28.90 16 17 Si 0.88 22.04 27.76 -169.99 -4.72 17.32 16 18 H -0.29 -7.07 7.11 56.52 0.40 -6.67 16 19 H -0.29 -7.41 7.11 56.52 0.40 -7.01 16 20 H -0.28 -6.46 6.51 56.52 0.37 -6.09 16 21 C -0.16 -2.14 4.83 -26.49 -0.13 -2.27 16 22 C 0.06 0.70 6.67 37.31 0.25 0.95 16 23 O -0.38 -4.99 10.77 -35.23 -0.38 -5.37 16 24 C 0.07 0.81 6.68 37.21 0.25 1.06 16 25 C -0.16 -2.23 4.77 -26.59 -0.13 -2.35 16 26 O -0.37 -5.02 8.71 -50.55 -0.44 -5.46 16 27 C 0.06 0.80 6.53 37.16 0.24 1.05 16 28 C -0.15 -1.56 10.28 37.16 0.38 -1.18 16 29 H 0.10 0.90 8.04 -52.49 -0.42 0.48 16 30 H 0.10 1.12 7.84 -52.48 -0.41 0.71 16 31 H 0.06 0.53 8.14 -51.93 -0.42 0.10 16 32 H 0.07 0.68 7.72 -51.93 -0.40 0.28 16 33 H 0.06 0.49 8.06 -51.93 -0.42 0.07 16 34 H 0.10 1.01 8.07 -51.93 -0.42 0.59 16 35 H 0.10 1.15 8.14 -51.93 -0.42 0.73 16 36 H 0.39 4.61 7.51 -40.82 -0.31 4.30 16 37 H 0.05 1.22 6.15 -51.93 -0.32 0.90 16 38 H 0.03 0.64 8.14 -51.93 -0.42 0.22 16 39 H 0.36 8.09 4.58 -40.82 -0.19 7.90 16 40 H 0.26 7.80 8.31 -40.82 -0.34 7.46 16 41 H 0.08 1.25 7.43 -51.93 -0.39 0.87 16 42 H 0.09 1.26 8.14 -51.93 -0.42 0.84 16 43 H 0.09 0.93 8.14 -51.93 -0.42 0.51 16 44 H 0.05 0.46 7.08 -51.93 -0.37 0.09 16 45 H 0.09 0.96 8.14 -51.93 -0.42 0.54 16 46 H 0.05 0.56 7.75 -51.93 -0.40 0.16 16 47 H 0.08 1.08 8.14 -51.93 -0.42 0.66 16 48 H 0.09 1.44 6.36 -51.93 -0.33 1.11 16 49 H 0.06 0.80 7.98 -51.93 -0.41 0.39 16 50 H 0.04 0.45 8.14 -51.93 -0.42 0.03 16 51 H 0.07 0.61 8.14 -51.93 -0.42 0.19 16 52 H 0.06 0.55 8.14 -51.93 -0.42 0.12 16 53 H 0.06 0.76 8.14 -51.93 -0.42 0.33 16 LS Contribution 413.71 15.07 6.23 6.23 Total: -1.00 -33.79 413.71 -7.99 -41.78 By element: Atomic # 1 Polarization: 22.34 SS G_CDS: -9.40 Total: 12.94 kcal Atomic # 6 Polarization: -1.11 SS G_CDS: 0.52 Total: -0.59 kcal Atomic # 7 Polarization: -57.54 SS G_CDS: 0.31 Total: -57.23 kcal Atomic # 8 Polarization: -19.52 SS G_CDS: -0.93 Total: -20.45 kcal Atomic # 14 Polarization: 22.04 SS G_CDS: -4.72 Total: 17.32 kcal Total LS contribution 6.23 Total: 6.23 kcal Total: -33.79 -7.99 -41.78 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. ZINC001502626342.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 -58.301 kcal (2) G-P(sol) polarization free energy of solvation -33.786 kcal (3) G-ENP(sol) elect.-nuc.-pol. free energy of system -92.086 kcal (4) G-CDS(sol) cavity-dispersion-solvent structure free energy -7.992 kcal (5) G-P-CDS(sol) = G-P(sol) + G-CDS(sol) = (2) + (4) -41.778 kcal (6) G-S(sol) free energy of system = (1) + (5) -100.078 kcal 1SCF WAS SPECIFIED, SO GEOMETRY OPTIMIZATION WAS NOT USED Total computer time = 0.72 seconds