Wall clock time and date at job start Wed Apr 1 2020 03:28: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= -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.53001 * 1 3 3 C 1.52990 * 109.46990 * 2 1 4 4 C 1.53010 * 109.47145 * 180.02562 * 3 2 1 5 5 H 1.09000 * 109.46510 * 305.92644 * 4 3 2 6 6 C 1.50695 * 109.47465 * 185.92699 * 4 3 2 7 7 O 1.21288 * 119.99935 * 247.05144 * 6 4 3 8 8 N 1.34770 * 120.00381 * 67.04310 * 6 4 3 9 9 C 1.46498 * 120.00255 * 185.49874 * 8 6 4 10 10 C 1.46500 * 119.99782 * 5.50061 * 8 6 4 11 11 C 1.53042 * 109.46184 * 234.38192 * 10 8 6 12 12 C 1.53196 * 109.31279 * 178.60334 * 11 10 8 13 13 N 1.46926 * 108.76784 * 54.65682 * 12 11 10 14 14 C 1.36929 * 120.63358 * 126.40788 * 13 12 11 15 15 H 1.08008 * 119.99515 * 328.43273 * 14 13 12 16 16 C 1.38807 * 120.00403 * 148.43928 * 14 13 12 17 17 N 1.31624 * 119.99879 * 349.09380 * 16 14 13 18 18 Si 1.86801 * 120.00288 * 169.10219 * 16 14 13 19 19 H 1.48500 * 109.47360 * 89.99370 * 18 16 14 20 20 H 1.48502 * 109.47139 * 209.99457 * 18 16 14 21 21 H 1.48504 * 109.46934 * 329.99375 * 18 16 14 22 22 C 1.46933 * 118.73630 * 306.38298 * 13 12 11 23 23 C 1.53039 * 109.46037 * 114.32797 * 10 8 6 24 24 C 1.50699 * 109.47045 * 65.91956 * 4 3 2 25 25 C 1.38227 * 119.61522 * 85.00392 * 24 4 3 26 26 C 1.38627 * 119.16288 * 180.02562 * 25 24 4 27 27 C 1.38679 * 118.41981 * 359.97438 * 26 25 24 28 28 C 1.38187 * 119.16207 * 0.27807 * 27 26 25 29 29 N 1.31848 * 119.62183 * 265.02284 * 24 4 3 30 30 H 1.08997 * 109.47568 * 180.02562 * 1 2 3 31 31 H 1.09002 * 109.47288 * 300.00061 * 1 2 3 32 32 H 1.08991 * 109.47356 * 59.99783 * 1 2 3 33 33 H 1.09002 * 109.46548 * 239.99861 * 2 1 3 34 34 H 1.09000 * 109.47191 * 120.00361 * 2 1 3 35 35 H 1.09002 * 109.47588 * 60.00499 * 3 2 1 36 36 H 1.08998 * 109.47498 * 299.99659 * 3 2 1 37 37 H 1.09006 * 109.47199 * 85.26075 * 9 8 6 38 38 H 1.09006 * 109.47202 * 205.25789 * 9 8 6 39 39 H 1.08993 * 109.47313 * 325.25908 * 9 8 6 40 40 H 1.08996 * 109.46219 * 354.35603 * 10 8 6 41 41 H 1.08996 * 109.49548 * 298.56559 * 11 10 8 42 42 H 1.08994 * 109.49515 * 58.66563 * 11 10 8 43 43 H 1.08998 * 109.59560 * 294.86755 * 12 11 10 44 44 H 1.09000 * 109.69907 * 174.50566 * 12 11 10 45 45 H 0.97000 * 119.99928 * 180.02562 * 17 16 14 46 46 H 1.09000 * 109.58366 * 293.82335 * 22 13 12 47 47 H 1.08994 * 109.70176 * 173.46830 * 22 13 12 48 48 H 1.08998 * 109.50155 * 301.33583 * 23 10 8 49 49 H 1.09000 * 109.50022 * 61.40143 * 23 10 8 50 50 H 1.08004 * 120.42080 * 0.04105 * 25 24 4 51 51 H 1.07999 * 120.79123 * 180.02562 * 26 25 24 52 52 H 1.07994 * 120.42035 * 180.02562 * 27 26 25 53 53 H 1.07994 * 119.62162 * 179.70121 * 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.5300 0.0000 0.0000 3 6 2.0399 1.4424 0.0000 4 6 3.5700 1.4425 -0.0006 5 1 3.9336 0.8393 0.8313 6 6 4.0725 2.8556 0.1461 7 8 4.6604 3.3885 -0.7713 8 7 3.8684 3.5274 1.2964 9 6 4.2486 4.9385 1.3980 10 6 3.2665 2.8464 2.4454 11 6 2.0425 3.6328 2.9204 12 6 1.3955 2.8994 4.0997 13 7 2.4130 2.6827 5.1372 14 6 2.1949 3.0937 6.4250 15 1 1.5919 3.9689 6.6176 16 6 2.7488 2.3859 7.4828 17 7 3.2988 1.2090 7.2709 18 14 2.7135 3.1035 9.2071 19 1 1.4699 2.6750 9.8963 20 1 3.8932 2.6200 9.9687 21 1 2.7495 4.5859 9.1272 22 6 3.6688 2.0158 4.7668 23 6 4.2901 2.7589 3.5798 24 6 4.0718 0.8629 -1.2981 25 6 4.2004 1.6791 -2.4063 26 6 4.6623 1.1360 -3.5952 27 6 4.9797 -0.2136 -3.6250 28 6 4.8205 -0.9683 -2.4784 29 7 4.3828 -0.4167 -1.3632 30 1 -0.3634 -1.0276 0.0005 31 1 -0.3634 0.5138 0.8900 32 1 -0.3633 0.5138 -0.8899 33 1 1.8932 -0.5139 -0.8900 34 1 1.8934 -0.5139 0.8899 35 1 1.6767 1.9563 0.8900 36 1 1.6767 1.9563 -0.8900 37 1 3.4389 5.5613 1.0175 38 1 4.4400 5.1891 2.4415 39 1 5.1496 5.1153 0.8107 40 1 2.9616 1.8416 2.1532 41 1 2.3506 4.6294 3.2364 42 1 1.3241 3.7149 2.1050 43 1 1.0061 1.9386 3.7633 44 1 0.5825 3.5015 4.5053 45 1 3.6855 0.7142 8.0101 46 1 3.4646 0.9829 4.4848 47 1 4.3578 2.0348 5.6110 48 1 5.1715 2.2190 3.2338 49 1 4.5772 3.7638 3.8892 50 1 3.9443 2.7266 -2.3467 51 1 4.7732 1.7496 -4.4770 52 1 5.3428 -0.6700 -4.5340 53 1 5.0609 -2.0210 -2.4946 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 ZINC001265195686.mol2 53 1SCF WAS SPECIFIED, SO GEOMETRY OPTIMIZATION WAS NOT USED SCF FIELD WAS ACHIEVED Wall clock time and date at job start Wed Apr 1 2020 03:28:50 Heat of formation + Delta-G solvation = 20.742523 kcal Electronic energy + Delta-G solvation = -32101.647953 eV Core-core repulsion = 28166.845209 eV Total energy + Delta-G solvation = -3934.802744 eV No. of doubly occupied orbitals = 66 Molecular weight (most abundant/longest-lived isotopes) = 345.216 amu Computer time = 1.28 seconds Orbital eigenvalues (eV) -41.00502 -39.57694 -37.68527 -35.84003 -35.08353 -32.52645 -32.34514 -32.06886 -30.23275 -29.67973 -29.53540 -26.94340 -25.27893 -24.69972 -23.56631 -22.57652 -22.27522 -22.03766 -20.93163 -19.41731 -18.91031 -17.27999 -17.03392 -16.63604 -16.11442 -15.85886 -15.43033 -15.29844 -14.88044 -14.77821 -14.23526 -13.93635 -13.71041 -13.58923 -13.54316 -13.32996 -13.06372 -12.95940 -12.79514 -12.60870 -12.58512 -12.31276 -12.12811 -11.93862 -11.55327 -11.45756 -11.14408 -11.13915 -11.04857 -10.82996 -10.77559 -10.58837 -10.39535 -10.28350 -10.14787 -9.93832 -9.83429 -9.73270 -9.47778 -9.28295 -8.93004 -8.59527 -8.36168 -6.69028 -5.83949 -3.32065 1.09088 1.31198 2.61999 3.05170 3.34062 3.41287 3.42442 3.79897 4.15776 4.38507 4.44279 4.78889 4.82444 4.89741 4.96034 5.05719 5.08191 5.14850 5.20837 5.33137 5.36783 5.37602 5.42916 5.49924 5.56547 5.59002 5.64026 5.65919 5.69395 5.82567 5.83432 5.89460 5.91520 5.97658 6.08476 6.12778 6.14913 6.20982 6.24410 6.31505 6.33800 6.40027 6.46478 6.56033 6.67574 6.70192 6.73522 6.88331 7.28821 7.50251 7.62516 7.74485 7.97439 8.21309 8.42181 9.16584 9.77702 10.35251 11.24421 Molecular weight = 345.22amu Principal moments of inertia in cm(-1) A = 0.015985 B = 0.002919 C = 0.002792 Principal moments of inertia in units of 10**(-40)*gram-cm**2 A = 1751.270411 B = 9590.123102 C =10026.468081 Net atomic charges, atomic populations, and dipole contributions calculated with CM2 Atom NO. Type Charge No. of electrons 1 C -0.151 4.151 2 C -0.119 4.119 3 C -0.110 4.110 4 C -0.060 4.060 5 H 0.116 0.884 6 C 0.522 3.478 7 O -0.547 6.547 8 N -0.601 5.601 9 C 0.086 3.914 10 C 0.134 3.866 11 C -0.141 4.141 12 C 0.148 3.852 13 N -0.591 5.591 14 C -0.255 4.255 15 H 0.064 0.936 16 C 0.006 3.994 17 N -0.903 5.903 18 Si 0.903 3.097 19 H -0.286 1.286 20 H -0.290 1.290 21 H -0.279 1.279 22 C 0.175 3.825 23 C -0.150 4.150 24 C 0.139 3.861 25 C -0.135 4.135 26 C -0.089 4.089 27 C -0.177 4.177 28 C 0.088 3.912 29 N -0.475 5.475 30 H 0.054 0.946 31 H 0.058 0.942 32 H 0.049 0.951 33 H 0.063 0.937 34 H 0.069 0.931 35 H 0.086 0.914 36 H 0.066 0.934 37 H 0.052 0.948 38 H 0.075 0.925 39 H 0.069 0.931 40 H 0.093 0.907 41 H 0.064 0.936 42 H 0.054 0.946 43 H 0.020 0.980 44 H 0.058 0.942 45 H 0.255 0.745 46 H 0.021 0.979 47 H 0.091 0.909 48 H 0.051 0.949 49 H 0.060 0.940 50 H 0.157 0.843 51 H 0.130 0.870 52 H 0.129 0.871 53 H 0.153 0.847 For the charges calculated by CM2: Dipole moment (debyes) X Y Z Total from point charges -0.412 -0.258 -19.351 19.357 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.209 4.209 2 C -0.157 4.157 3 C -0.148 4.148 4 C -0.082 4.082 5 H 0.134 0.866 6 C 0.309 3.691 7 O -0.426 6.426 8 N -0.331 5.331 9 C -0.056 4.056 10 C 0.031 3.969 11 C -0.181 4.181 12 C 0.020 3.980 13 N -0.314 5.314 14 C -0.385 4.385 15 H 0.082 0.918 16 C -0.193 4.193 17 N -0.547 5.547 18 Si 0.733 3.267 19 H -0.213 1.213 20 H -0.216 1.216 21 H -0.205 1.205 22 C 0.048 3.952 23 C -0.191 4.191 24 C 0.001 3.999 25 C -0.155 4.155 26 C -0.115 4.115 27 C -0.197 4.197 28 C -0.069 4.069 29 N -0.183 5.183 30 H 0.073 0.927 31 H 0.077 0.923 32 H 0.068 0.932 33 H 0.082 0.918 34 H 0.088 0.912 35 H 0.104 0.896 36 H 0.085 0.915 37 H 0.070 0.930 38 H 0.094 0.906 39 H 0.088 0.912 40 H 0.111 0.889 41 H 0.083 0.917 42 H 0.073 0.927 43 H 0.038 0.962 44 H 0.076 0.924 45 H 0.065 0.935 46 H 0.039 0.961 47 H 0.110 0.890 48 H 0.070 0.930 49 H 0.078 0.922 50 H 0.174 0.826 51 H 0.148 0.852 52 H 0.147 0.853 53 H 0.170 0.830 Dipole moment (debyes) X Y Z Total from point charges -0.116 -0.588 -18.923 18.932 hybrid contribution 0.293 0.775 0.175 0.847 sum 0.177 0.187 -18.747 18.749 Atomic orbital electron populations 1.21761 0.95724 1.01640 1.01759 1.21488 0.95663 0.95591 1.02941 1.21734 0.90658 0.98828 1.03544 1.20145 0.97731 0.93778 0.96504 0.86649 1.20291 0.78906 0.87752 0.82125 1.90583 1.42483 1.70277 1.39251 1.48120 1.60820 1.11670 1.12527 1.21572 1.00481 0.80845 1.02727 1.21530 0.94043 0.96366 0.84920 1.22419 0.96198 1.00353 0.99103 1.20880 0.91877 0.96550 0.88657 1.44022 1.12718 1.69394 1.05244 1.19136 1.31039 1.06348 0.81950 0.91834 1.25164 1.01029 0.95069 0.98019 1.69961 1.42143 1.11730 1.30882 0.86290 0.77285 0.78617 0.84500 1.21312 1.21567 1.20478 1.20764 0.88466 0.92162 0.93760 1.22643 0.97223 1.02134 0.97066 1.20503 0.95762 0.89457 0.94139 1.22601 0.98200 1.02442 0.92233 1.21631 0.96707 0.94987 0.98149 1.22111 1.02307 0.95229 1.00017 1.22625 0.96076 1.00089 0.88125 1.67380 1.10495 1.12302 1.28163 0.92733 0.92293 0.93173 0.91831 0.91187 0.89592 0.91523 0.92956 0.90593 0.91195 0.88911 0.91698 0.92704 0.96198 0.92390 0.93513 0.96098 0.89048 0.92998 0.92155 0.82562 0.85235 0.85277 0.82991 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.15 -1.56 9.91 37.16 0.37 -1.19 16 2 C -0.12 -1.38 5.73 -26.74 -0.15 -1.53 16 3 C -0.11 -1.30 4.20 -26.73 -0.11 -1.41 16 4 C -0.06 -0.81 1.65 -92.92 -0.15 -0.96 16 5 H 0.12 1.65 7.26 -51.93 -0.38 1.27 16 6 C 0.52 7.71 6.84 -10.99 -0.08 7.64 16 7 O -0.55 -9.02 12.49 5.55 0.07 -8.95 16 8 N -0.60 -8.56 2.76 -173.89 -0.48 -9.04 16 9 C 0.09 1.12 9.35 59.85 0.56 1.68 16 10 C 0.13 2.04 1.49 -67.89 -0.10 1.93 16 11 C -0.14 -2.17 4.95 -26.61 -0.13 -2.30 16 12 C 0.15 2.80 6.52 -3.71 -0.02 2.77 16 13 N -0.59 -13.76 2.98 -172.31 -0.51 -14.27 16 14 C -0.26 -6.65 9.58 -15.06 -0.14 -6.79 16 15 H 0.06 1.63 7.85 -52.48 -0.41 1.22 16 16 C 0.01 0.16 5.04 -17.43 -0.09 0.07 16 17 N -0.90 -25.22 11.08 55.49 0.61 -24.60 16 18 Si 0.90 21.02 29.57 -169.99 -5.03 15.99 16 19 H -0.29 -6.68 7.11 56.52 0.40 -6.28 16 20 H -0.29 -6.81 7.11 56.52 0.40 -6.41 16 21 H -0.28 -6.38 7.11 56.52 0.40 -5.98 16 22 C 0.17 3.89 5.25 -3.70 -0.02 3.87 16 23 C -0.15 -2.67 5.41 -26.61 -0.14 -2.82 16 24 C 0.14 1.93 4.96 -82.59 -0.41 1.52 16 25 C -0.14 -1.77 8.34 -39.44 -0.33 -2.10 16 26 C -0.09 -0.96 10.18 -39.28 -0.40 -1.36 16 27 C -0.18 -1.74 10.11 -39.44 -0.40 -2.14 16 28 C 0.09 0.98 11.18 -17.55 -0.20 0.78 16 29 N -0.47 -6.48 9.62 -9.47 -0.09 -6.57 16 30 H 0.05 0.53 8.14 -51.93 -0.42 0.10 16 31 H 0.06 0.61 8.14 -51.93 -0.42 0.19 16 32 H 0.05 0.48 8.14 -51.93 -0.42 0.05 16 33 H 0.06 0.78 6.92 -51.93 -0.36 0.42 16 34 H 0.07 0.86 8.14 -51.93 -0.42 0.44 16 35 H 0.09 1.04 5.04 -51.93 -0.26 0.78 16 36 H 0.07 0.76 8.14 -51.93 -0.42 0.34 16 37 H 0.05 0.62 8.14 -51.93 -0.42 0.20 16 38 H 0.08 0.97 6.59 -51.93 -0.34 0.63 16 39 H 0.07 0.92 7.50 -51.93 -0.39 0.53 16 40 H 0.09 1.40 4.51 -51.93 -0.23 1.17 16 41 H 0.06 1.04 7.46 -51.93 -0.39 0.65 16 42 H 0.05 0.74 7.52 -51.93 -0.39 0.35 16 43 H 0.02 0.38 8.14 -51.93 -0.42 -0.04 16 44 H 0.06 1.09 7.93 -51.93 -0.41 0.68 16 45 H 0.25 6.93 8.31 -40.82 -0.34 6.59 16 46 H 0.02 0.49 8.14 -51.93 -0.42 0.07 16 47 H 0.09 2.33 6.13 -51.93 -0.32 2.01 16 48 H 0.05 0.89 8.14 -51.93 -0.42 0.46 16 49 H 0.06 1.10 6.97 -51.93 -0.36 0.73 16 50 H 0.16 2.16 5.69 -52.48 -0.30 1.86 16 51 H 0.13 1.08 8.06 -52.49 -0.42 0.66 16 52 H 0.13 0.89 8.06 -52.49 -0.42 0.46 16 53 H 0.15 1.27 8.06 -52.49 -0.42 0.84 16 LS Contribution 403.63 15.07 6.08 6.08 Total: -1.00 -29.66 403.63 -10.05 -39.71 By element: Atomic # 1 Polarization: 12.74 SS G_CDS: -8.75 Total: 3.99 kcal Atomic # 6 Polarization: -0.39 SS G_CDS: -1.95 Total: -2.34 kcal Atomic # 7 Polarization: -54.02 SS G_CDS: -0.47 Total: -54.49 kcal Atomic # 8 Polarization: -9.02 SS G_CDS: 0.07 Total: -8.95 kcal Atomic # 14 Polarization: 21.02 SS G_CDS: -5.03 Total: 15.99 kcal Total LS contribution 6.08 Total: 6.08 kcal Total: -29.66 -10.05 -39.71 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. ZINC001265195686.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 60.456 kcal (2) G-P(sol) polarization free energy of solvation -29.664 kcal (3) G-ENP(sol) elect.-nuc.-pol. free energy of system 30.792 kcal (4) G-CDS(sol) cavity-dispersion-solvent structure free energy -10.049 kcal (5) G-P-CDS(sol) = G-P(sol) + G-CDS(sol) = (2) + (4) -39.714 kcal (6) G-S(sol) free energy of system = (1) + (5) 20.743 kcal 1SCF WAS SPECIFIED, SO GEOMETRY OPTIMIZATION WAS NOT USED Total computer time = 1.28 seconds