Wall clock time and date at job start Wed Apr 1 2020 00:29:13 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.52995 * 1 3 3 H 1.09002 * 109.49928 * 2 1 4 4 C 1.53192 * 109.49831 * 239.93190 * 2 1 3 5 5 N 1.46927 * 108.77419 * 185.37947 * 4 2 1 6 6 C 1.34777 * 120.62966 * 233.58380 * 5 4 2 7 7 O 1.21274 * 120.00145 * 0.02562 * 6 5 4 8 8 C 1.50696 * 119.99735 * 180.02562 * 6 5 4 9 9 C 1.52998 * 109.46973 * 180.02562 * 8 6 5 10 10 C 1.50706 * 109.47348 * 180.02562 * 9 8 6 11 11 H 1.08005 * 119.99758 * 359.97438 * 10 9 8 12 12 C 1.37871 * 120.00102 * 179.97438 * 10 9 8 13 13 N 1.31516 * 120.00217 * 0.02562 * 12 10 9 14 14 Si 1.86800 * 119.99766 * 179.97438 * 12 10 9 15 15 H 1.48501 * 109.47057 * 179.97438 * 14 12 10 16 16 H 1.48500 * 109.47389 * 300.00068 * 14 12 10 17 17 H 1.48504 * 109.47102 * 60.00417 * 14 12 10 18 18 C 1.46922 * 118.73621 * 53.86712 * 5 4 2 19 19 C 1.53629 * 108.54335 * 306.08668 * 18 5 4 20 20 C 1.53044 * 109.23737 * 54.86124 * 19 18 5 21 21 H 1.08998 * 109.46325 * 58.52941 * 20 19 18 22 22 N 1.46499 * 109.44482 * 178.49888 * 20 19 18 23 23 C 1.34779 * 120.00011 * 275.05787 * 22 20 19 24 24 O 1.21511 * 119.99753 * 0.02562 * 23 22 20 25 25 O 1.34633 * 120.00293 * 180.02562 * 23 22 20 26 26 C 1.45198 * 117.00064 * 180.02562 * 25 23 22 27 27 C 1.52992 * 109.47386 * 59.99644 * 26 25 23 28 28 C 1.53001 * 109.47088 * 179.97438 * 26 25 23 29 29 C 1.53003 * 109.47109 * 299.99478 * 26 25 23 30 30 H 1.08999 * 109.47189 * 180.02562 * 1 2 3 31 31 H 1.09005 * 109.47489 * 300.06702 * 1 2 3 32 32 H 1.09000 * 109.47114 * 60.06483 * 1 2 3 33 33 H 1.09002 * 109.58112 * 65.59310 * 4 2 1 34 34 H 1.09000 * 109.58855 * 305.17416 * 4 2 1 35 35 H 1.09000 * 109.47333 * 60.00571 * 8 6 5 36 36 H 1.09002 * 109.47152 * 300.00044 * 8 6 5 37 37 H 1.09001 * 109.47671 * 60.00215 * 9 8 6 38 38 H 1.09006 * 109.47297 * 299.99934 * 9 8 6 39 39 H 0.97001 * 119.99646 * 359.97438 * 13 12 10 40 40 H 1.08998 * 109.58902 * 65.92698 * 18 5 4 41 41 H 1.09003 * 109.58673 * 186.34728 * 18 5 4 42 42 H 1.08992 * 109.63384 * 294.86996 * 19 18 5 43 43 H 1.09001 * 109.45497 * 174.74363 * 19 18 5 44 44 H 0.96998 * 119.99985 * 95.05995 * 22 20 19 45 45 H 1.09000 * 109.47483 * 60.00135 * 27 26 25 46 46 H 1.09007 * 109.47330 * 180.02562 * 27 26 25 47 47 H 1.09000 * 109.47558 * 299.99786 * 27 26 25 48 48 H 1.08994 * 109.47173 * 59.99696 * 28 26 25 49 49 H 1.08997 * 109.46993 * 179.97438 * 28 26 25 50 50 H 1.09000 * 109.46795 * 299.99261 * 28 26 25 51 51 H 1.09001 * 109.47110 * 59.99989 * 29 26 25 52 52 H 1.08997 * 109.47551 * 180.02562 * 29 26 25 53 53 H 1.09006 * 109.47219 * 300.00763 * 29 26 25 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 1 1.8938 1.0275 0.0000 4 6 2.0413 -0.7235 -1.2497 5 7 3.5047 -0.8281 -1.1701 6 6 4.2807 -0.3850 -2.1790 7 8 3.7722 0.1061 -3.1645 8 6 5.7802 -0.4999 -2.0826 9 6 6.4180 0.0729 -3.3499 10 6 7.9176 -0.0425 -3.2538 11 1 8.3704 -0.4802 -2.3763 12 6 8.7114 0.4107 -4.2859 13 7 8.1601 0.9441 -5.3542 14 14 10.5701 0.2683 -4.1664 15 1 11.1892 0.8230 -5.3970 16 1 11.0499 1.0285 -2.9844 17 1 10.9489 -1.1605 -4.0232 18 6 4.1056 -1.4306 0.0276 19 6 3.5714 -0.6942 1.2655 20 6 2.0414 -0.7230 1.2482 21 1 1.6977 -1.7573 1.2351 22 7 1.5274 -0.0529 2.4453 23 6 1.4000 -0.7382 3.5988 24 8 1.7118 -1.9117 3.6460 25 8 0.9271 -0.1227 4.6989 26 6 0.8172 -0.9223 5.9059 27 6 2.1987 -1.4531 6.2936 28 6 0.2663 -0.0579 7.0418 29 6 -0.1288 -2.0982 5.6543 30 1 -0.3633 -1.0276 0.0005 31 1 -0.3634 0.5149 0.8894 32 1 -0.3633 0.5128 -0.8906 33 1 1.7632 -0.1577 -2.1389 34 1 1.6056 -1.7215 -1.2982 35 1 6.0583 -1.5485 -1.9776 36 1 6.1322 0.0581 -1.2150 37 1 6.1400 1.1216 -3.4546 38 1 6.0660 -0.4847 -4.2178 39 1 7.1949 1.0184 -5.4160 40 1 3.8350 -2.4851 0.0812 41 1 5.1902 -1.3316 -0.0185 42 1 3.9160 0.3398 1.2535 43 1 3.9345 -1.1883 2.1667 44 1 1.2788 0.8839 2.4077 45 1 2.8727 -0.6154 6.4729 46 1 2.1162 -2.0537 7.1995 47 1 2.5917 -2.0688 5.4845 48 1 -0.7178 0.3206 6.7659 49 1 0.1834 -0.6583 7.9477 50 1 0.9406 0.7795 7.2210 51 1 0.2641 -2.7138 4.8451 52 1 -0.2117 -2.6982 6.5604 53 1 -1.1131 -1.7200 5.3780 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 ZINC000496642011.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 00:29:13 Heat of formation + Delta-G solvation = -142.372547 kcal Electronic energy + Delta-G solvation = -30746.450086 eV Core-core repulsion = 26617.314090 eV Total energy + Delta-G solvation = -4129.135996 eV No. of doubly occupied orbitals = 66 Molecular weight (most abundant/longest-lived isotopes) = 340.220 amu Computer time = 1.10 seconds Orbital eigenvalues (eV) -41.18683 -40.11585 -38.41321 -36.94914 -34.97000 -34.29476 -33.86960 -33.08683 -31.06905 -30.78597 -27.79491 -27.73473 -27.21583 -25.84351 -24.49168 -22.88033 -22.31781 -21.72887 -20.80781 -20.11106 -18.77463 -18.27177 -17.33020 -16.93338 -16.43345 -16.21306 -15.78887 -15.62833 -15.28610 -15.02660 -14.73814 -14.42792 -13.96659 -13.86269 -13.73535 -13.61825 -13.52636 -13.07145 -12.87624 -12.60768 -12.48478 -12.30624 -12.26740 -12.25326 -12.16693 -12.04926 -11.88899 -11.79800 -11.64941 -11.39971 -11.31312 -11.23753 -10.82224 -10.74764 -10.68922 -10.32818 -10.12758 -10.11534 -9.89417 -9.45412 -9.25390 -8.55594 -8.52676 -8.27845 -6.31291 -4.01426 1.85092 2.57053 3.08350 3.15686 3.45661 3.48048 3.66851 3.79521 4.23773 4.33507 4.37828 4.48923 4.52360 4.57437 4.60613 4.85625 5.00227 5.03158 5.08243 5.10067 5.15993 5.22090 5.28257 5.30747 5.37255 5.40071 5.44884 5.46317 5.48530 5.53994 5.57626 5.60578 5.68420 5.77262 5.84910 5.88652 5.89894 5.92719 6.01468 6.07728 6.09771 6.14544 6.24098 6.60908 6.75403 7.12763 7.51850 7.57164 7.61846 8.10629 8.19774 8.49856 8.61386 8.92574 9.72164 11.09593 Molecular weight = 340.22amu Principal moments of inertia in cm(-1) A = 0.017868 B = 0.002111 C = 0.001948 Principal moments of inertia in units of 10**(-40)*gram-cm**2 A = 1566.647699 B =13262.621272 C =14369.004863 Net atomic charges, atomic populations, and dipole contributions calculated with CM2 Atom NO. Type Charge No. of electrons 1 C -0.145 4.145 2 C -0.113 4.113 3 H 0.083 0.917 4 C 0.120 3.880 5 N -0.626 5.626 6 C 0.514 3.486 7 O -0.536 6.536 8 C -0.122 4.122 9 C 0.024 3.976 10 C -0.608 4.608 11 H 0.049 0.951 12 C 0.015 3.985 13 N -0.910 5.910 14 Si 1.070 2.930 15 H -0.282 1.282 16 H -0.288 1.288 17 H -0.288 1.288 18 C 0.112 3.888 19 C -0.147 4.147 20 C 0.170 3.830 21 H 0.094 0.906 22 N -0.705 5.705 23 C 0.651 3.349 24 O -0.567 6.567 25 O -0.367 6.367 26 C 0.137 3.863 27 C -0.184 4.184 28 C -0.141 4.141 29 C -0.183 4.183 30 H 0.058 0.942 31 H 0.058 0.942 32 H 0.063 0.937 33 H 0.098 0.902 34 H 0.066 0.934 35 H 0.081 0.919 36 H 0.081 0.919 37 H 0.014 0.986 38 H 0.014 0.986 39 H 0.272 0.728 40 H 0.070 0.930 41 H 0.097 0.903 42 H 0.072 0.928 43 H 0.081 0.919 44 H 0.405 0.595 45 H 0.059 0.941 46 H 0.070 0.930 47 H 0.088 0.912 48 H 0.066 0.934 49 H 0.080 0.920 50 H 0.067 0.933 51 H 0.088 0.912 52 H 0.070 0.930 53 H 0.058 0.942 For the charges calculated by CM2: Dipole moment (debyes) X Y Z Total from point charges -19.297 -4.107 25.275 32.064 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.203 4.203 2 C -0.133 4.133 3 H 0.101 0.899 4 C -0.003 4.003 5 N -0.362 5.362 6 C 0.305 3.695 7 O -0.412 6.412 8 C -0.162 4.162 9 C -0.013 4.013 10 C -0.646 4.646 11 H 0.067 0.933 12 C -0.190 4.190 13 N -0.546 5.546 14 Si 0.898 3.102 15 H -0.208 1.208 16 H -0.215 1.215 17 H -0.215 1.215 18 C -0.011 4.011 19 C -0.187 4.187 20 C 0.066 3.934 21 H 0.112 0.888 22 N -0.363 5.363 23 C 0.404 3.596 24 O -0.456 6.456 25 O -0.282 6.282 26 C 0.100 3.900 27 C -0.241 4.241 28 C -0.198 4.198 29 C -0.240 4.240 30 H 0.077 0.923 31 H 0.077 0.923 32 H 0.082 0.918 33 H 0.116 0.884 34 H 0.084 0.916 35 H 0.100 0.900 36 H 0.099 0.901 37 H 0.032 0.968 38 H 0.033 0.967 39 H 0.082 0.918 40 H 0.088 0.912 41 H 0.115 0.885 42 H 0.091 0.909 43 H 0.100 0.900 44 H 0.240 0.760 45 H 0.079 0.921 46 H 0.088 0.912 47 H 0.107 0.893 48 H 0.085 0.915 49 H 0.099 0.901 50 H 0.086 0.914 51 H 0.107 0.893 52 H 0.088 0.912 53 H 0.077 0.923 Dipole moment (debyes) X Y Z Total from point charges -18.015 -4.116 24.018 30.304 hybrid contribution -1.582 0.028 1.126 1.942 sum -19.597 -4.088 25.144 32.140 Atomic orbital electron populations 1.21768 0.93793 1.02436 1.02265 1.21825 0.96588 1.00852 0.94061 0.89886 1.21503 0.80589 1.00454 0.97801 1.48326 1.06039 1.62642 1.19188 1.21573 0.90715 0.75860 0.81385 1.90646 1.75771 1.43834 1.30924 1.21744 0.88143 1.02766 1.03522 1.19026 0.95319 0.96002 0.90984 1.22024 0.91340 1.45056 1.06224 0.93286 1.25848 1.06744 0.93690 0.92689 1.70738 1.14032 1.41234 1.28612 0.83540 0.74456 0.76276 0.75921 1.20780 1.21482 1.21466 1.21760 0.99579 0.96673 0.83135 1.22455 0.95019 1.02288 0.98922 1.21027 0.93323 0.97524 0.81500 0.88780 1.44872 1.68803 1.14806 1.07790 1.18096 0.76709 0.83820 0.80991 1.90905 1.54866 1.15165 1.84640 1.86301 1.78647 1.44362 1.18856 1.22062 0.96498 0.90601 0.80838 1.22518 0.95691 1.02316 1.03578 1.21865 1.01845 1.00176 0.95942 1.22516 0.99041 0.98274 1.04207 0.92270 0.92299 0.91819 0.88390 0.91599 0.90042 0.90087 0.96773 0.96734 0.91752 0.91167 0.88503 0.90879 0.90010 0.76048 0.92144 0.91151 0.89281 0.91515 0.90124 0.91443 0.89311 0.91156 0.92255 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.14 -1.08 9.02 37.16 0.34 -0.75 16 2 C -0.11 -1.03 2.46 -90.50 -0.22 -1.26 16 3 H 0.08 0.81 8.14 -51.93 -0.42 0.39 16 4 C 0.12 1.42 5.71 -3.71 -0.02 1.40 16 5 N -0.63 -8.87 2.97 -172.50 -0.51 -9.38 16 6 C 0.51 9.86 7.70 -10.99 -0.08 9.77 16 7 O -0.54 -12.38 15.39 5.56 0.09 -12.30 16 8 C -0.12 -2.49 3.98 -27.89 -0.11 -2.61 16 9 C 0.02 0.65 3.57 -27.88 -0.10 0.55 16 10 C -0.61 -17.50 9.11 -34.44 -0.31 -17.81 16 11 H 0.05 1.34 7.74 -52.48 -0.41 0.94 16 12 C 0.01 0.43 5.46 -17.82 -0.10 0.33 16 13 N -0.91 -28.84 13.65 55.43 0.76 -28.08 16 14 Si 1.07 25.51 29.90 -169.99 -5.08 20.43 16 15 H -0.28 -6.76 7.11 56.52 0.40 -6.36 16 16 H -0.29 -6.51 7.11 56.52 0.40 -6.10 16 17 H -0.29 -6.51 7.11 56.52 0.40 -6.11 16 18 C 0.11 1.27 6.37 -3.48 -0.02 1.25 16 19 C -0.15 -1.47 5.67 -26.38 -0.15 -1.62 16 20 C 0.17 1.56 1.88 -67.89 -0.13 1.43 16 21 H 0.09 0.96 7.58 -51.93 -0.39 0.57 16 22 N -0.70 -6.05 4.37 -58.69 -0.26 -6.31 16 23 C 0.65 6.77 8.02 54.06 0.43 7.20 16 24 O -0.57 -7.03 11.54 -19.28 -0.22 -7.26 16 25 O -0.37 -3.51 9.94 -40.31 -0.40 -3.91 16 26 C 0.14 1.08 1.13 -90.62 -0.10 0.98 16 27 C -0.18 -1.47 8.37 37.15 0.31 -1.15 16 28 C -0.14 -0.78 8.85 37.16 0.33 -0.45 16 29 C -0.18 -1.41 8.37 37.16 0.31 -1.10 16 30 H 0.06 0.44 8.14 -51.93 -0.42 0.01 16 31 H 0.06 0.39 7.37 -51.93 -0.38 0.01 16 32 H 0.06 0.47 8.14 -51.93 -0.42 0.05 16 33 H 0.10 1.35 7.01 -51.93 -0.36 0.98 16 34 H 0.07 0.69 8.14 -51.93 -0.42 0.27 16 35 H 0.08 1.54 7.93 -51.93 -0.41 1.13 16 36 H 0.08 1.53 7.70 -51.93 -0.40 1.13 16 37 H 0.01 0.41 7.97 -51.93 -0.41 -0.01 16 38 H 0.01 0.42 7.97 -51.93 -0.41 0.00 16 39 H 0.27 8.74 7.22 -40.82 -0.29 8.45 16 40 H 0.07 0.71 8.14 -51.93 -0.42 0.29 16 41 H 0.10 1.12 5.93 -51.93 -0.31 0.81 16 42 H 0.07 0.76 8.14 -51.93 -0.42 0.34 16 43 H 0.08 0.83 8.14 -51.93 -0.42 0.41 16 44 H 0.40 2.82 8.42 -40.82 -0.34 2.47 16 45 H 0.06 0.47 8.14 -51.93 -0.42 0.05 16 46 H 0.07 0.44 8.14 -51.92 -0.42 0.01 16 47 H 0.09 0.91 5.88 -51.93 -0.31 0.61 16 48 H 0.07 0.35 8.14 -51.93 -0.42 -0.07 16 49 H 0.08 0.31 8.14 -51.93 -0.42 -0.11 16 50 H 0.07 0.37 8.14 -51.93 -0.42 -0.06 16 51 H 0.09 0.89 5.88 -51.93 -0.31 0.58 16 52 H 0.07 0.42 8.14 -51.93 -0.42 0.00 16 53 H 0.06 0.43 8.14 -51.93 -0.42 0.01 16 LS Contribution 413.35 15.07 6.23 6.23 Total: -1.00 -35.22 413.35 -8.49 -43.71 By element: Atomic # 1 Polarization: 10.15 SS G_CDS: -9.46 Total: 0.69 kcal Atomic # 6 Polarization: -4.19 SS G_CDS: 0.37 Total: -3.82 kcal Atomic # 7 Polarization: -43.76 SS G_CDS: -0.01 Total: -43.77 kcal Atomic # 8 Polarization: -22.93 SS G_CDS: -0.54 Total: -23.47 kcal Atomic # 14 Polarization: 25.51 SS G_CDS: -5.08 Total: 20.43 kcal Total LS contribution 6.23 Total: 6.23 kcal Total: -35.22 -8.49 -43.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. ZINC000496642011.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 -98.661 kcal (2) G-P(sol) polarization free energy of solvation -35.220 kcal (3) G-ENP(sol) elect.-nuc.-pol. free energy of system -133.881 kcal (4) G-CDS(sol) cavity-dispersion-solvent structure free energy -8.491 kcal (5) G-P-CDS(sol) = G-P(sol) + G-CDS(sol) = (2) + (4) -43.711 kcal (6) G-S(sol) free energy of system = (1) + (5) -142.373 kcal 1SCF WAS SPECIFIED, SO GEOMETRY OPTIMIZATION WAS NOT USED Total computer time = 1.11 seconds