Wall clock time and date at job start Sat Apr 11 2020 06:10:31 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 O 1.42903 * 1 3 3 C 1.42908 * 113.99750 * 2 1 4 4 C 1.52993 * 109.47148 * 179.97438 * 3 2 1 5 5 N 1.46497 * 109.47354 * 65.00184 * 4 3 2 6 6 C 1.37099 * 120.00420 * 265.88236 * 5 4 3 7 7 H 1.07994 * 120.00716 * 138.14939 * 6 5 4 8 8 C 1.38954 * 119.99551 * 318.14686 * 6 5 4 9 9 N 1.31410 * 119.99989 * 341.74524 * 8 6 5 10 10 Si 1.86800 * 119.99900 * 161.74030 * 8 6 5 11 11 H 1.48500 * 109.46918 * 60.00362 * 10 8 6 12 12 H 1.48501 * 109.46832 * 180.02562 * 10 8 6 13 13 H 1.48494 * 109.47266 * 300.00343 * 10 8 6 14 14 C 1.34775 * 119.99775 * 85.87388 * 5 4 3 15 15 O 1.21290 * 120.00072 * 5.51208 * 14 5 4 16 16 C 1.50702 * 120.00244 * 185.50441 * 14 5 4 17 17 H 1.09005 * 109.46666 * 37.97355 * 16 14 5 18 18 C 1.53000 * 109.47347 * 277.98285 * 16 14 5 19 19 N 1.46498 * 109.47311 * 157.97696 * 16 14 5 20 20 C 1.34780 * 119.99813 * 205.00358 * 19 16 14 21 21 O 1.21282 * 120.00003 * 359.97438 * 20 19 16 22 22 C 1.50704 * 119.99870 * 179.97438 * 20 19 16 23 23 C 1.52998 * 109.46961 * 179.97438 * 22 20 19 24 24 C 1.53005 * 109.47472 * 295.00311 * 23 22 20 25 25 C 1.52995 * 109.47129 * 180.02562 * 24 23 22 26 26 C 1.52996 * 109.47366 * 59.99801 * 25 24 23 27 27 C 1.52997 * 109.47207 * 299.99996 * 26 25 24 28 28 C 1.52996 * 109.46865 * 175.00209 * 23 22 20 29 29 H 1.08993 * 109.47159 * 299.99710 * 1 2 3 30 30 H 1.09003 * 109.47162 * 59.99959 * 1 2 3 31 31 H 1.09002 * 109.46576 * 179.97438 * 1 2 3 32 32 H 1.09000 * 109.46632 * 59.99965 * 3 2 1 33 33 H 1.08998 * 109.47518 * 300.00342 * 3 2 1 34 34 H 1.09006 * 109.47399 * 305.00511 * 4 3 2 35 35 H 1.09004 * 109.47241 * 185.00497 * 4 3 2 36 36 H 0.97003 * 119.99823 * 179.97438 * 9 8 6 37 37 H 1.09001 * 109.47405 * 61.68487 * 18 16 14 38 38 H 1.09000 * 109.47175 * 181.69207 * 18 16 14 39 39 H 1.09007 * 109.46576 * 301.68746 * 18 16 14 40 40 H 0.96999 * 120.00129 * 25.00407 * 19 16 14 41 41 H 1.08992 * 109.46899 * 300.00395 * 22 20 19 42 42 H 1.09006 * 109.46579 * 60.00117 * 22 20 19 43 43 H 1.09002 * 109.47404 * 54.99962 * 23 22 20 44 44 H 1.09002 * 109.46521 * 300.00191 * 24 23 22 45 45 H 1.08997 * 109.46985 * 59.99382 * 24 23 22 46 46 H 1.09002 * 109.47245 * 299.99344 * 25 24 23 47 47 H 1.09006 * 109.46866 * 180.02562 * 25 24 23 48 48 H 1.08999 * 109.46982 * 179.97438 * 26 25 24 49 49 H 1.09002 * 109.47222 * 60.00133 * 26 25 24 50 50 H 1.09002 * 109.47141 * 300.00186 * 27 26 25 51 51 H 1.08999 * 109.47490 * 180.02562 * 27 26 25 52 52 H 1.08997 * 109.47502 * 299.99879 * 28 23 22 53 53 H 1.09002 * 109.47222 * 60.00246 * 28 23 22 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.4290 0.0000 0.0000 3 6 2.0102 1.3056 0.0000 4 6 3.5354 1.1848 0.0006 5 7 3.9767 0.5642 -1.2508 6 6 4.1794 -0.7904 -1.3101 7 1 3.8471 -1.3490 -2.1726 8 6 4.8125 -1.4447 -0.2604 9 7 5.5272 -0.7608 0.6046 10 14 4.6474 -3.2965 -0.0788 11 1 5.1842 -3.9627 -1.2926 12 1 5.4111 -3.7470 1.1124 13 1 3.2153 -3.6518 0.0885 14 6 4.1833 1.3250 -2.3440 15 8 4.1013 2.5326 -2.2668 16 6 4.5163 0.6751 -3.6623 17 1 3.9407 -0.2442 -3.7699 18 6 6.0106 0.3493 -3.7073 19 7 4.1831 1.5909 -4.7562 20 6 3.9082 1.1040 -5.9825 21 8 3.9372 -0.0922 -6.1806 22 6 3.5648 2.0460 -7.1076 23 6 3.2993 1.2425 -8.3823 24 6 4.5865 0.5409 -8.8204 25 6 4.3210 -0.2620 -10.0954 26 6 3.8520 0.6820 -11.2043 27 6 2.5649 1.3836 -10.7663 28 6 2.8300 2.1862 -9.4914 29 1 -0.3633 0.5138 0.8899 30 1 -0.3633 0.5138 -0.8900 31 1 -0.3632 -1.0277 -0.0005 32 1 1.6885 1.8464 0.8900 33 1 1.6886 1.8465 -0.8899 34 1 3.8508 0.5686 0.8427 35 1 3.9784 2.1768 0.0903 36 1 5.9688 -1.2175 1.3376 37 1 6.2529 -0.3555 -2.9118 38 1 6.2570 -0.0933 -4.6724 39 1 6.5863 1.2646 -3.5696 40 1 4.1602 2.5476 -4.5978 41 1 2.6739 2.6162 -6.8449 42 1 4.3971 2.7294 -7.2761 43 1 2.5275 0.4977 -8.1883 44 1 5.3582 1.2858 -9.0145 45 1 4.9208 -0.1313 -8.0302 46 1 3.5491 -1.0069 -9.9017 47 1 5.2380 -0.7622 -10.4073 48 1 3.6634 0.1100 -12.1128 49 1 4.6239 1.4268 -11.3981 50 1 1.7930 0.6388 -10.5725 51 1 2.2305 2.0559 -11.5564 52 1 1.9130 2.6860 -9.1794 53 1 3.6018 2.9310 -9.6855 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 ZINC000457184826.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:31 Heat of formation + Delta-G solvation = -143.963120 kcal Electronic energy + Delta-G solvation = -31310.630941 eV Core-core repulsion = 27181.425973 eV Total energy + Delta-G solvation = -4129.204968 eV No. of doubly occupied orbitals = 66 Molecular weight (most abundant/longest-lived isotopes) = 340.220 amu Computer time = 1.04 seconds Orbital eigenvalues (eV) -39.95882 -38.90660 -37.91804 -37.01273 -34.90288 -34.18226 -32.71929 -31.63566 -30.74457 -30.52059 -28.86155 -26.93997 -25.86263 -24.85509 -23.69312 -22.65902 -22.09412 -21.14475 -19.81104 -19.70863 -19.14410 -17.56240 -16.69335 -15.92221 -15.56536 -15.46431 -15.37720 -15.20309 -14.94053 -14.42004 -14.33163 -14.24162 -14.16100 -14.10096 -13.84164 -13.34580 -13.14907 -12.76409 -12.46636 -12.21186 -12.14278 -12.02162 -11.84670 -11.60211 -11.55605 -11.38627 -11.17447 -11.13786 -11.06859 -11.01669 -10.97838 -10.89830 -10.58553 -10.47300 -10.36069 -10.27094 -10.07875 -10.05877 -9.28608 -9.16597 -9.01428 -8.69824 -8.37033 -7.50867 -6.35249 -4.07308 3.01560 3.20715 3.27037 3.32549 3.50670 3.80634 4.19289 4.46423 4.59648 4.67441 4.68076 4.76390 5.01547 5.06370 5.07755 5.20012 5.33727 5.45574 5.46654 5.53750 5.58464 5.62096 5.63819 5.68536 5.72695 5.80773 5.85944 5.93832 5.96356 5.99124 6.01426 6.03864 6.06740 6.14365 6.15432 6.16711 6.25021 6.26191 6.43894 6.49546 6.52419 6.59493 6.65775 6.87138 7.05410 7.18934 7.48560 7.78088 7.94685 8.02346 8.34117 8.65929 9.08096 9.41382 9.83944 10.68668 Molecular weight = 340.22amu Principal moments of inertia in cm(-1) A = 0.013606 B = 0.002821 C = 0.002652 Principal moments of inertia in units of 10**(-40)*gram-cm**2 A = 2057.396582 B = 9922.151167 C =10556.570823 Net atomic charges, atomic populations, and dipole contributions calculated with CM2 Atom NO. Type Charge No. of electrons 1 C 0.031 3.969 2 O -0.364 6.364 3 C 0.052 3.948 4 C 0.119 3.881 5 N -0.419 5.419 6 C -0.331 4.331 7 H 0.087 0.913 8 C 0.054 3.946 9 N -0.863 5.863 10 Si 0.913 3.087 11 H -0.272 1.272 12 H -0.286 1.286 13 H -0.271 1.271 14 C 0.470 3.530 15 O -0.611 6.611 16 C 0.142 3.858 17 H 0.104 0.896 18 C -0.144 4.144 19 N -0.712 5.712 20 C 0.510 3.490 21 O -0.543 6.543 22 C -0.137 4.137 23 C -0.069 4.069 24 C -0.120 4.120 25 C -0.117 4.117 26 C -0.120 4.120 27 C -0.118 4.118 28 C -0.115 4.115 29 H 0.030 0.970 30 H 0.031 0.969 31 H 0.078 0.922 32 H 0.040 0.960 33 H 0.039 0.961 34 H 0.108 0.892 35 H 0.072 0.928 36 H 0.270 0.730 37 H 0.087 0.913 38 H 0.053 0.947 39 H 0.049 0.951 40 H 0.404 0.596 41 H 0.090 0.910 42 H 0.090 0.910 43 H 0.078 0.922 44 H 0.053 0.947 45 H 0.085 0.915 46 H 0.061 0.939 47 H 0.057 0.943 48 H 0.056 0.944 49 H 0.059 0.941 50 H 0.061 0.939 51 H 0.058 0.942 52 H 0.060 0.940 53 H 0.061 0.939 For the charges calculated by CM2: Dipole moment (debyes) X Y Z Total from point charges -4.846 4.514 -13.716 15.231 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.066 4.066 2 O -0.282 6.282 3 C -0.025 4.025 4 C -0.005 4.005 5 N -0.136 5.136 6 C -0.459 4.459 7 H 0.105 0.895 8 C -0.152 4.152 9 N -0.499 5.499 10 Si 0.740 3.260 11 H -0.197 1.197 12 H -0.211 1.211 13 H -0.196 1.196 14 C 0.258 3.742 15 O -0.498 6.498 16 C 0.035 3.965 17 H 0.122 0.878 18 C -0.202 4.202 19 N -0.365 5.365 20 C 0.296 3.704 21 O -0.421 6.421 22 C -0.177 4.177 23 C -0.088 4.088 24 C -0.158 4.158 25 C -0.155 4.155 26 C -0.157 4.157 27 C -0.156 4.156 28 C -0.153 4.153 29 H 0.049 0.951 30 H 0.050 0.950 31 H 0.096 0.904 32 H 0.058 0.942 33 H 0.057 0.943 34 H 0.126 0.874 35 H 0.091 0.909 36 H 0.080 0.920 37 H 0.106 0.894 38 H 0.072 0.928 39 H 0.068 0.932 40 H 0.240 0.760 41 H 0.108 0.892 42 H 0.108 0.892 43 H 0.096 0.904 44 H 0.071 0.929 45 H 0.104 0.896 46 H 0.080 0.920 47 H 0.076 0.924 48 H 0.075 0.925 49 H 0.077 0.923 50 H 0.080 0.920 51 H 0.077 0.923 52 H 0.079 0.921 53 H 0.080 0.920 Dipole moment (debyes) X Y Z Total from point charges -4.631 4.023 -13.572 14.894 hybrid contribution -0.382 -0.197 0.429 0.607 sum -5.013 3.825 -13.143 14.577 Atomic orbital electron populations 1.22649 0.83358 1.01821 0.98804 1.88151 1.17424 1.27574 1.95027 1.23040 0.96462 0.84432 0.98612 1.21396 0.89105 0.99117 0.90910 1.43184 1.55643 1.07072 1.07671 1.19228 1.37888 0.81577 1.07191 0.89506 1.25290 0.93538 0.99865 0.96522 1.69911 1.27586 1.36990 1.15458 0.86222 0.78787 0.83156 0.77884 1.19710 1.21103 1.19581 1.20372 0.80414 0.86639 0.86808 1.90491 1.58499 1.14505 1.86316 1.20558 0.97839 0.93292 0.84788 0.87831 1.21893 0.93700 1.02033 1.02609 1.46004 1.69951 1.11333 1.09183 1.20598 0.76654 0.89557 0.83561 1.90650 1.52425 1.16185 1.82830 1.21745 1.03428 0.98567 0.93991 1.20890 0.98853 0.96901 0.92161 1.21715 0.96550 0.99580 0.97933 1.21500 1.00918 0.97946 0.95121 1.21551 0.97458 0.98698 0.98009 1.21499 0.97770 0.99166 0.97125 1.21598 1.00951 0.97822 0.94894 0.95126 0.95025 0.90360 0.94157 0.94289 0.87382 0.90943 0.92038 0.89447 0.92763 0.93213 0.76030 0.89166 0.89174 0.90405 0.92873 0.89629 0.92029 0.92380 0.92507 0.92262 0.92036 0.92340 0.92091 0.92030 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.03 0.49 11.01 101.05 1.11 1.60 16 2 O -0.36 -7.64 10.37 -35.23 -0.37 -8.01 16 3 C 0.05 1.06 6.26 37.15 0.23 1.29 16 4 C 0.12 2.75 5.35 -4.05 -0.02 2.73 16 5 N -0.42 -10.25 2.45 -114.81 -0.28 -10.53 16 6 C -0.33 -8.31 7.34 -14.93 -0.11 -8.42 16 7 H 0.09 2.18 6.64 -52.49 -0.35 1.83 16 8 C 0.05 1.36 4.96 -17.47 -0.09 1.27 16 9 N -0.86 -22.43 11.65 55.50 0.65 -21.78 16 10 Si 0.91 20.00 29.59 -169.99 -5.03 14.97 16 11 H -0.27 -5.88 7.11 56.52 0.40 -5.48 16 12 H -0.29 -6.21 7.11 56.52 0.40 -5.81 16 13 H -0.27 -6.00 7.11 56.52 0.40 -5.60 16 14 C 0.47 11.01 6.85 -10.99 -0.08 10.93 16 15 O -0.61 -14.94 15.57 5.55 0.09 -14.86 16 16 C 0.14 2.95 2.27 -69.08 -0.16 2.79 16 17 H 0.10 2.32 6.16 -51.93 -0.32 2.00 16 18 C -0.14 -2.83 9.38 37.16 0.35 -2.48 16 19 N -0.71 -12.45 5.05 -54.93 -0.28 -12.72 16 20 C 0.51 8.26 7.51 -10.98 -0.08 8.18 16 21 O -0.54 -10.43 13.63 5.56 0.08 -10.35 16 22 C -0.14 -1.61 5.17 -27.88 -0.14 -1.75 16 23 C -0.07 -0.79 1.84 -90.62 -0.17 -0.95 16 24 C -0.12 -1.43 4.95 -26.73 -0.13 -1.57 16 25 C -0.12 -1.22 5.92 -26.74 -0.16 -1.38 16 26 C -0.12 -1.05 5.92 -26.73 -0.16 -1.21 16 27 C -0.12 -1.01 5.92 -26.75 -0.16 -1.17 16 28 C -0.12 -1.04 5.20 -26.76 -0.14 -1.18 16 29 H 0.03 0.42 8.14 -51.93 -0.42 0.00 16 30 H 0.03 0.47 8.14 -51.93 -0.42 0.05 16 31 H 0.08 1.12 8.14 -51.93 -0.42 0.70 16 32 H 0.04 0.71 8.14 -51.93 -0.42 0.29 16 33 H 0.04 0.80 8.14 -51.93 -0.42 0.38 16 34 H 0.11 2.70 5.79 -51.93 -0.30 2.39 16 35 H 0.07 1.68 7.48 -51.93 -0.39 1.29 16 36 H 0.27 6.69 8.31 -40.82 -0.34 6.35 16 37 H 0.09 1.88 7.23 -51.93 -0.38 1.50 16 38 H 0.05 1.00 8.14 -51.93 -0.42 0.58 16 39 H 0.05 0.95 8.14 -51.92 -0.42 0.53 16 40 H 0.40 6.65 7.75 -40.82 -0.32 6.34 16 41 H 0.09 0.93 8.14 -51.93 -0.42 0.50 16 42 H 0.09 0.91 8.14 -51.93 -0.42 0.48 16 43 H 0.08 1.02 7.95 -51.93 -0.41 0.61 16 44 H 0.05 0.59 8.14 -51.93 -0.42 0.17 16 45 H 0.09 1.26 5.95 -51.93 -0.31 0.96 16 46 H 0.06 0.68 8.14 -51.93 -0.42 0.26 16 47 H 0.06 0.58 8.14 -51.93 -0.42 0.15 16 48 H 0.06 0.45 8.14 -51.93 -0.42 0.03 16 49 H 0.06 0.50 8.14 -51.93 -0.42 0.08 16 50 H 0.06 0.55 8.14 -51.93 -0.42 0.12 16 51 H 0.06 0.44 8.14 -51.93 -0.42 0.01 16 52 H 0.06 0.51 8.14 -51.93 -0.42 0.09 16 53 H 0.06 0.52 8.14 -51.93 -0.42 0.10 16 LS Contribution 415.31 15.07 6.26 6.26 Total: -1.00 -29.16 415.31 -8.30 -37.46 By element: Atomic # 1 Polarization: 20.41 SS G_CDS: -9.52 Total: 10.89 kcal Atomic # 6 Polarization: 8.57 SS G_CDS: 0.10 Total: 8.67 kcal Atomic # 7 Polarization: -45.12 SS G_CDS: 0.09 Total: -45.04 kcal Atomic # 8 Polarization: -33.01 SS G_CDS: -0.20 Total: -33.22 kcal Atomic # 14 Polarization: 20.00 SS G_CDS: -5.03 Total: 14.97 kcal Total LS contribution 6.26 Total: 6.26 kcal Total: -29.16 -8.30 -37.46 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. ZINC000457184826.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 -106.500 kcal (2) G-P(sol) polarization free energy of solvation -29.164 kcal (3) G-ENP(sol) elect.-nuc.-pol. free energy of system -135.664 kcal (4) G-CDS(sol) cavity-dispersion-solvent structure free energy -8.300 kcal (5) G-P-CDS(sol) = G-P(sol) + G-CDS(sol) = (2) + (4) -37.463 kcal (6) G-S(sol) free energy of system = (1) + (5) -143.963 kcal 1SCF WAS SPECIFIED, SO GEOMETRY OPTIMIZATION WAS NOT USED Total computer time = 1.04 seconds