Wall clock time and date at job start Tue Mar 31 2020 06:13:37 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 * SM5.42R - SM5.42R CALCULATIONS WILL BE PERFORMED * - CHARGE MODEL 2 WILL BE USED (CM2A) * - THE SOLVENT IS WATER ******************************************************************************* 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 N 1.46496 * 1 3 3 C 1.46507 * 120.00140 * 2 1 4 4 C 1.50694 * 109.46960 * 270.00184 * 3 2 1 5 5 H 1.08007 * 120.00305 * 0.02562 * 4 3 2 6 6 C 1.37880 * 120.00116 * 179.97438 * 4 3 2 7 7 N 1.31519 * 119.99625 * 359.97438 * 6 4 3 8 8 Si 1.86797 * 120.00302 * 180.02562 * 6 4 3 9 9 H 1.48498 * 109.47115 * 89.99880 * 8 6 4 10 10 H 1.48500 * 109.47062 * 209.99874 * 8 6 4 11 11 H 1.48496 * 109.47160 * 329.99831 * 8 6 4 12 12 C 1.39679 * 120.00034 * 179.97438 * 2 1 3 13 13 C 1.38904 * 120.07158 * 0.02562 * 12 2 1 14 14 C 1.38108 * 119.92469 * 179.97438 * 13 12 2 15 15 C 1.38287 * 120.07052 * 359.97438 * 14 13 12 16 16 C 1.50703 * 119.92391 * 180.02562 * 15 14 13 17 17 N 1.46902 * 109.47272 * 89.99722 * 16 15 14 18 18 C 1.46895 * 110.99793 * 179.97438 * 17 16 15 19 19 C 1.53007 * 109.47103 * 180.02562 * 18 17 16 20 20 H 1.08998 * 109.46674 * 59.99801 * 19 18 17 21 21 O 1.42898 * 109.46891 * 300.00235 * 19 18 17 22 22 C 1.52997 * 109.47128 * 180.02562 * 19 18 17 23 23 N 1.46908 * 109.46719 * 184.99648 * 22 19 18 24 24 C 1.47024 * 110.99261 * 66.32173 * 23 22 19 25 25 C 1.53103 * 109.24985 * 182.34870 * 24 23 22 26 26 O 1.43012 * 109.29295 * 303.21723 * 25 24 23 27 27 C 1.43011 * 113.74271 * 58.65008 * 26 25 24 28 28 C 1.47018 * 110.99661 * 189.99712 * 23 22 19 29 29 C 1.38283 * 120.15337 * 0.27387 * 15 14 13 30 30 C 1.38109 * 120.07182 * 359.45125 * 29 15 14 31 31 H 1.08995 * 109.47079 * 270.00048 * 1 2 3 32 32 H 1.08998 * 109.47236 * 29.99637 * 1 2 3 33 33 H 1.09000 * 109.47410 * 150.00029 * 1 2 3 34 34 H 1.09001 * 109.46681 * 30.00126 * 3 2 1 35 35 H 1.08995 * 109.47260 * 150.00036 * 3 2 1 36 36 H 0.96998 * 119.99719 * 179.97438 * 7 6 4 37 37 H 1.08000 * 120.04140 * 359.95258 * 13 12 2 38 38 H 1.07993 * 119.96821 * 179.97438 * 14 13 12 39 39 H 1.08994 * 109.46957 * 329.99908 * 16 15 14 40 40 H 1.09001 * 109.46999 * 210.00090 * 16 15 14 41 41 H 1.00903 * 110.99813 * 56.04663 * 17 16 15 42 42 H 1.08998 * 109.47088 * 60.00284 * 18 17 16 43 43 H 1.09006 * 109.47632 * 299.99814 * 18 17 16 44 44 H 0.96701 * 113.99894 * 300.00253 * 21 19 18 45 45 H 1.09001 * 109.47059 * 65.00273 * 22 19 18 46 46 H 1.08995 * 109.46804 * 304.99511 * 22 19 18 47 47 H 1.08998 * 109.47369 * 62.43595 * 24 23 22 48 48 H 1.08998 * 109.54483 * 302.29802 * 24 23 22 49 49 H 1.09002 * 109.50635 * 183.25881 * 25 24 23 50 50 H 1.09003 * 109.50698 * 63.17493 * 25 24 23 51 51 H 1.08997 * 109.50806 * 61.26171 * 27 26 25 52 52 H 1.08997 * 109.50924 * 181.35636 * 27 26 25 53 53 H 1.08994 * 109.50209 * 57.67478 * 28 23 22 54 54 H 1.08997 * 109.50514 * 297.54292 * 28 23 22 55 55 H 1.07995 * 119.96775 * 179.70690 * 29 15 14 56 56 H 1.08002 * 120.03683 * 180.27478 * 30 29 15 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 7 1.4650 0.0000 0.0000 3 6 2.1975 1.2688 0.0000 4 6 2.4486 1.7038 -1.4208 5 1 2.0973 1.0962 -2.2417 6 6 3.1260 2.8780 -1.6728 7 7 3.5525 3.6186 -0.6731 8 14 3.4366 3.4176 -3.4339 9 1 2.3062 4.2602 -3.9002 10 1 4.6967 4.2009 -3.4966 11 1 3.5540 2.2207 -4.3050 12 6 2.1634 -1.2097 0.0005 13 6 1.4703 -2.4135 0.0005 14 6 2.1640 -3.6077 0.0005 15 6 3.5468 -3.6059 0.0000 16 6 4.3004 -4.9110 -0.0006 17 7 4.5446 -5.3348 -1.3858 18 6 5.2794 -6.6062 -1.4254 19 6 5.5150 -7.0139 -2.8812 20 1 6.0969 -6.2420 -3.3850 21 8 4.2574 -7.1696 -3.5417 22 6 6.2797 -8.3384 -2.9226 23 7 6.6105 -8.6695 -4.3151 24 6 7.5667 -7.7024 -4.8736 25 6 7.8378 -8.0499 -6.3398 26 8 8.3144 -9.3955 -6.4262 27 6 7.4070 -10.3632 -5.8921 28 6 7.1322 -10.0399 -4.4210 29 6 4.2398 -2.4093 0.0052 30 6 3.5524 -1.2114 -0.0005 31 1 -0.3633 0.0000 1.0276 32 1 -0.3633 0.8900 -0.5138 33 1 -0.3634 -0.8900 -0.5138 34 1 1.6084 2.0284 0.5139 35 1 3.1499 1.1384 0.5138 36 1 4.0293 4.4445 -0.8504 37 1 0.3904 -2.4156 0.0009 38 1 1.6259 -4.5441 0.0010 39 1 3.7115 -5.6707 0.5133 40 1 5.2530 -4.7808 0.5129 41 1 5.0333 -4.6188 -1.9021 42 1 4.6977 -7.3779 -0.9213 43 1 6.2389 -6.4873 -0.9220 44 1 3.6859 -7.8418 -3.1460 45 1 5.6605 -9.1288 -2.4984 46 1 7.1982 -8.2458 -2.3431 47 1 8.4990 -7.7467 -4.3106 48 1 7.1490 -6.6977 -4.8090 49 1 8.5898 -7.3709 -6.7418 50 1 6.9160 -7.9528 -6.9134 51 1 6.4722 -10.3336 -6.4519 52 1 7.8470 -11.3573 -5.9708 53 1 6.3972 -10.7407 -4.0254 54 1 8.0578 -10.1220 -3.8513 55 1 5.3198 -2.4113 0.0095 56 1 4.0943 -0.2772 -0.0010 RHF calculation, no. of doubly occupied orbitals= 70 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=WATER ZINC000023624943.mol2 56 1SCF WAS SPECIFIED, SO GEOMETRY OPTIMIZATION WAS NOT USED SCF FIELD WAS ACHIEVED Wall clock time and date at job start Tue Mar 31 2020 06:13:37 Heat of formation + Delta-G solvation = -73.619419 kcal Electronic energy + Delta-G solvation = -32075.140487 eV Core-core repulsion = 27790.252232 eV Total energy + Delta-G solvation = -4284.888255 eV No. of doubly occupied orbitals = 70 Molecular weight (most abundant/longest-lived isotopes) = 363.231 amu Computer time = 1.20 seconds Orbital eigenvalues (eV) -40.68432 -40.43540 -38.93742 -37.65187 -36.34846 -35.68069 -34.69488 -32.24271 -31.97168 -31.87617 -30.58943 -29.60989 -28.28470 -26.92331 -25.69941 -25.07852 -24.29267 -23.41123 -22.91506 -22.34954 -21.35564 -20.99026 -18.63918 -18.22158 -17.86871 -17.43543 -17.12020 -16.86188 -16.76823 -16.40727 -16.28710 -16.04373 -15.79011 -15.70147 -15.44177 -15.35823 -15.22492 -14.88833 -14.76313 -14.40930 -14.27055 -14.05943 -13.97627 -13.44107 -13.28546 -13.18620 -13.05965 -13.03090 -12.92236 -12.69462 -12.60995 -12.53742 -12.44469 -12.30540 -12.19872 -11.95932 -11.74812 -11.66773 -11.54423 -11.39822 -10.92337 -10.64803 -10.57878 -10.38842 -9.65699 -9.36341 -9.23363 -8.36575 -7.79699 -6.37934 0.71886 0.85192 1.36116 1.37135 1.48538 2.33854 2.61229 2.64242 2.91749 3.04789 3.05763 3.44995 3.69167 3.86035 3.90493 3.97787 4.05184 4.10998 4.15317 4.17259 4.23045 4.26868 4.30073 4.36221 4.42411 4.47181 4.55044 4.59689 4.64756 4.67070 4.70842 4.71965 4.74402 4.77674 4.78645 4.87666 4.88597 4.93091 4.99383 5.02156 5.13765 5.15137 5.20469 5.34145 5.59448 5.67773 5.70365 5.89791 5.93067 5.94582 6.06932 6.14621 6.15476 6.24982 6.43623 6.50343 6.82688 6.94287 7.07065 7.33883 8.82477 Molecular weight = 363.23amu Principal moments of inertia in cm(-1) A = 0.010559 B = 0.001711 C = 0.001537 Principal moments of inertia in units of 10**(-40)*gram-cm**2 A = 2651.051900 B =16363.761906 C =18218.469774 Net atomic charges, atomic populations, and dipole contributions calculated with CM2 Atom NO. Type Charge No. of electrons 1 C 0.098 3.902 2 N -0.632 5.632 3 C 0.267 3.733 4 C -0.541 4.541 5 H 0.048 0.952 6 C 0.012 3.988 7 N -0.926 5.926 8 Si 0.961 3.039 9 H -0.270 1.270 10 H -0.278 1.278 11 H -0.273 1.273 12 C 0.212 3.788 13 C -0.219 4.219 14 C -0.048 4.048 15 C -0.159 4.159 16 C 0.122 3.878 17 N -0.685 5.685 18 C 0.032 3.968 19 C 0.121 3.879 20 H 0.080 0.920 21 O -0.556 6.556 22 C 0.030 3.970 23 N -0.538 5.538 24 C 0.016 3.984 25 C 0.062 3.938 26 O -0.407 6.407 27 C 0.059 3.941 28 C 0.019 3.981 29 C -0.069 4.069 30 C -0.217 4.217 31 H 0.053 0.947 32 H 0.056 0.944 33 H 0.077 0.923 34 H 0.009 0.991 35 H 0.011 0.989 36 H 0.267 0.733 37 H 0.132 0.868 38 H 0.141 0.859 39 H 0.087 0.913 40 H 0.036 0.964 41 H 0.337 0.663 42 H 0.088 0.912 43 H 0.042 0.958 44 H 0.381 0.619 45 H 0.096 0.904 46 H 0.047 0.953 47 H 0.039 0.961 48 H 0.086 0.914 49 H 0.108 0.892 50 H 0.056 0.944 51 H 0.058 0.942 52 H 0.115 0.885 53 H 0.102 0.898 54 H 0.039 0.961 55 H 0.123 0.877 56 H 0.095 0.905 For the charges calculated by CM2: Dipole moment (debyes) X Y Z Total from point charges 5.742 -30.312 -3.027 30.999 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.045 4.045 2 N -0.340 5.340 3 C 0.145 3.855 4 C -0.580 4.580 5 H 0.066 0.934 6 C -0.187 4.187 7 N -0.565 5.565 8 Si 0.786 3.214 9 H -0.196 1.196 10 H -0.203 1.203 11 H -0.198 1.198 12 C 0.107 3.893 13 C -0.242 4.242 14 C -0.067 4.067 15 C -0.162 4.162 16 C -0.006 4.006 17 N -0.316 5.316 18 C -0.098 4.098 19 C 0.061 3.939 20 H 0.098 0.902 21 O -0.362 6.362 22 C -0.098 4.098 23 N -0.263 5.263 24 C -0.112 4.112 25 C -0.015 4.015 26 O -0.326 6.326 27 C -0.018 4.018 28 C -0.108 4.108 29 C -0.089 4.089 30 C -0.240 4.240 31 H 0.071 0.929 32 H 0.074 0.926 33 H 0.096 0.904 34 H 0.027 0.973 35 H 0.029 0.971 36 H 0.076 0.924 37 H 0.149 0.851 38 H 0.158 0.842 39 H 0.105 0.895 40 H 0.055 0.945 41 H 0.154 0.846 42 H 0.107 0.893 43 H 0.060 0.940 44 H 0.229 0.771 45 H 0.115 0.885 46 H 0.065 0.935 47 H 0.057 0.943 48 H 0.105 0.895 49 H 0.126 0.874 50 H 0.074 0.926 51 H 0.076 0.924 52 H 0.133 0.867 53 H 0.120 0.880 54 H 0.057 0.943 55 H 0.140 0.860 56 H 0.113 0.887 Dipole moment (debyes) X Y Z Total from point charges 5.084 -29.812 -3.774 30.477 hybrid contribution 1.490 1.205 0.506 1.982 sum 6.574 -28.608 -3.268 29.534 Atomic orbital electron populations 1.21596 0.79450 1.02491 1.00959 1.46960 1.05605 1.03714 1.77709 1.18767 0.91073 0.80276 0.95387 1.21778 1.36405 1.06487 0.93346 0.93376 1.25779 0.95103 0.96322 1.01487 1.70037 1.39009 1.14741 1.32747 0.85508 0.77966 0.78230 0.79709 1.19576 1.20321 1.19769 1.18236 0.91679 0.86516 0.92894 1.20823 1.00066 0.90803 1.12478 1.20883 0.93552 0.97032 0.95275 1.19534 0.93952 0.94073 1.08638 1.20630 0.97683 0.93021 0.89283 1.60437 1.52534 1.05409 1.13183 1.22653 0.97190 0.92997 0.96940 1.21783 0.85247 0.95484 0.91400 0.90199 1.86622 1.28384 1.63013 1.58155 1.22804 0.99959 0.98680 0.88366 1.62125 1.42496 1.02338 1.19380 1.22815 0.94336 0.97514 0.96508 1.22723 1.00369 0.83396 0.94965 1.87954 1.45506 1.16676 1.82458 1.22823 0.93889 0.92743 0.92319 1.22829 0.99548 0.90580 0.97887 1.20839 0.99302 0.91357 0.97360 1.20763 0.93213 0.98545 1.11495 0.92879 0.92559 0.90410 0.97290 0.97078 0.92396 0.85061 0.84170 0.89457 0.94539 0.84596 0.89326 0.94009 0.77115 0.88550 0.93502 0.94261 0.89533 0.87353 0.92594 0.92410 0.86723 0.87983 0.94312 0.85951 0.88723 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.10 3.47 10.19 127.77 1.30 4.77 16 2 N -0.63 -28.15 2.99 -726.58 -2.17 -30.32 16 3 C 0.27 14.21 4.91 85.63 0.42 14.63 16 4 C -0.54 -30.87 9.39 25.60 0.24 -30.63 16 5 H 0.05 2.79 7.81 -2.91 -0.02 2.77 16 6 C 0.01 0.70 5.46 84.65 0.46 1.16 16 7 N -0.93 -54.67 13.29 21.45 0.28 -54.39 16 8 Si 0.96 44.62 29.54 68.60 2.03 46.65 16 9 H -0.27 -12.08 7.11 99.48 0.71 -11.37 16 10 H -0.28 -12.49 7.11 99.48 0.71 -11.79 16 11 H -0.27 -12.58 7.11 99.48 0.71 -11.87 16 12 C 0.21 8.78 6.55 38.21 0.25 9.03 16 13 C -0.22 -7.38 9.21 22.40 0.21 -7.17 16 14 C -0.05 -1.39 9.69 22.25 0.22 -1.18 16 15 C -0.16 -4.51 5.40 -19.84 -0.11 -4.62 16 16 C 0.12 2.52 5.81 85.56 0.50 3.02 16 17 N -0.68 -13.77 6.36 -517.92 -3.29 -17.06 16 18 C 0.03 0.44 5.35 86.31 0.46 0.90 16 19 C 0.12 1.76 3.14 30.60 0.10 1.86 16 20 H 0.08 1.40 6.30 -2.39 -0.02 1.38 16 21 O -0.56 -10.42 13.15 -148.98 -1.96 -12.38 16 22 C 0.03 0.34 4.59 86.30 0.40 0.74 16 23 N -0.54 -7.54 4.67 -895.17 -4.18 -11.73 16 24 C 0.02 0.25 5.44 86.31 0.47 0.71 16 25 C 0.06 0.95 7.01 72.02 0.50 1.46 16 26 O -0.41 -7.15 10.83 -148.98 -1.61 -8.76 16 27 C 0.06 0.73 7.01 72.02 0.50 1.23 16 28 C 0.02 0.22 5.65 86.31 0.49 0.71 16 29 C -0.07 -2.33 9.69 22.25 0.22 -2.12 16 30 C -0.22 -9.17 9.21 22.40 0.21 -8.97 16 31 H 0.05 1.66 8.14 -2.39 -0.02 1.64 16 32 H 0.06 2.07 8.07 -2.39 -0.02 2.05 16 33 H 0.08 2.36 6.20 -2.39 -0.01 2.34 16 34 H 0.01 0.51 8.07 -2.39 -0.02 0.49 16 35 H 0.01 0.64 6.20 -2.39 -0.01 0.63 16 36 H 0.27 14.83 8.31 -92.71 -0.77 14.06 16 37 H 0.13 3.86 6.14 -2.91 -0.02 3.84 16 38 H 0.14 3.19 8.06 -2.91 -0.02 3.17 16 39 H 0.09 1.55 8.08 -2.39 -0.02 1.53 16 40 H 0.04 0.71 8.08 -2.39 -0.02 0.69 16 41 H 0.34 7.63 8.72 -89.69 -0.78 6.85 16 42 H 0.09 0.98 8.14 -2.39 -0.02 0.97 16 43 H 0.04 0.52 8.14 -2.38 -0.02 0.50 16 44 H 0.38 4.73 9.04 -74.06 -0.67 4.06 16 45 H 0.10 0.78 8.04 -2.39 -0.02 0.76 16 46 H 0.05 0.50 8.12 -2.39 -0.02 0.48 16 47 H 0.04 0.59 8.14 -2.39 -0.02 0.57 16 48 H 0.09 1.35 6.29 -2.39 -0.02 1.33 16 49 H 0.11 1.35 8.14 -2.39 -0.02 1.33 16 50 H 0.06 0.91 8.14 -2.39 -0.02 0.89 16 51 H 0.06 0.74 8.14 -2.39 -0.02 0.72 16 52 H 0.11 0.94 8.14 -2.39 -0.02 0.92 16 53 H 0.10 0.80 8.09 -2.39 -0.02 0.78 16 54 H 0.04 0.46 8.14 -2.39 -0.02 0.44 16 55 H 0.12 3.60 8.06 -2.91 -0.02 3.58 16 56 H 0.09 4.60 6.14 -2.91 -0.02 4.59 16 Total: -1.00 -69.47 442.93 -4.66 -74.13 By element: Atomic # 1 Polarization: 28.91 SS G_CDS: -0.58 Total: 28.33 kcal Atomic # 6 Polarization: -21.30 SS G_CDS: 6.83 Total: -14.47 kcal Atomic # 7 Polarization: -104.14 SS G_CDS: -9.36 Total: -113.50 kcal Atomic # 8 Polarization: -17.56 SS G_CDS: -3.57 Total: -21.14 kcal Atomic # 14 Polarization: 44.62 SS G_CDS: 2.03 Total: 46.65 kcal Total: -69.47 -4.66 -74.13 kcal The number of atoms in the molecule is 56 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. ZINC000023624943.mol2 56 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 0.510 kcal (2) G-P(sol) polarization free energy of solvation -69.471 kcal (3) G-ENP(sol) elect.-nuc.-pol. free energy of system -68.961 kcal (4) G-CDS(sol) cavity-dispersion-solvent structure free energy -4.658 kcal (5) G-P-CDS(sol) = G-P(sol) + G-CDS(sol) = (2) + (4) -74.130 kcal (6) G-S(sol) free energy of system = (1) + (5) -73.619 kcal 1SCF WAS SPECIFIED, SO GEOMETRY OPTIMIZATION WAS NOT USED Total computer time = 1.21 seconds