Wall clock time and date at job start Tue Mar 31 2020 16:58:12 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.52993 * 1 3 3 C 1.53003 * 109.47180 * 2 1 4 4 C 1.53002 * 109.47558 * 239.99865 * 2 1 3 5 5 C 1.47197 * 109.47054 * 119.99955 * 2 1 3 6 6 C 2.64599 * 109.47057 * 119.99864 * 2 1 3 7 7 C 4.07803 * 109.46971 * 119.99810 * 2 1 3 8 8 O 1.21610 * 119.99967 * 8.73989 * 7 6 5 9 9 N 1.34771 * 119.99903 * 188.73824 * 7 6 5 10 10 C 1.46930 * 120.63306 * 6.11231 * 9 7 6 11 11 C 1.53188 * 108.77477 * 233.58710 * 10 9 7 12 12 C 1.53036 * 109.31561 * 305.36105 * 11 10 9 13 13 C 1.53039 * 109.53803 * 61.36703 * 12 11 10 14 14 C 1.46923 * 120.63241 * 186.08916 * 9 7 6 15 15 H 1.08999 * 109.58619 * 246.34457 * 14 9 7 16 16 C 1.53005 * 109.58687 * 6.62356 * 14 9 7 17 17 H 1.08997 * 110.88742 * 304.40044 * 16 14 9 18 18 C 1.55161 * 111.00602 * 180.51635 * 16 14 9 19 19 C 1.54325 * 102.94077 * 206.42650 * 18 16 14 20 20 N 1.47015 * 107.27248 * 22.19481 * 19 18 16 21 21 C 1.36936 * 125.65092 * 181.02195 * 20 19 18 22 22 H 1.07992 * 120.00421 * 180.02562 * 21 20 19 23 23 C 1.38814 * 119.99621 * 359.75176 * 21 20 19 24 24 N 1.31619 * 119.99946 * 0.27210 * 23 21 20 25 25 Si 1.86799 * 119.99839 * 180.27608 * 23 21 20 26 26 H 1.48499 * 109.47193 * 359.97438 * 25 23 21 27 27 H 1.48508 * 109.46997 * 119.99969 * 25 23 21 28 28 H 1.48500 * 109.46859 * 239.99323 * 25 23 21 29 29 C 1.47431 * 108.70362 * 1.28579 * 20 19 18 30 30 H 1.08994 * 109.47572 * 300.00170 * 1 2 3 31 31 H 1.09005 * 109.46978 * 60.00595 * 1 2 3 32 32 H 1.08996 * 109.47430 * 179.97438 * 1 2 3 33 33 H 1.09006 * 109.46926 * 179.97438 * 3 2 1 34 34 H 1.09001 * 109.47147 * 299.99664 * 3 2 1 35 35 H 1.09000 * 109.47268 * 59.99787 * 3 2 1 36 36 H 1.09006 * 109.46752 * 300.00102 * 4 2 1 37 37 H 1.09001 * 109.47055 * 60.00160 * 4 2 1 38 38 H 1.08994 * 109.47248 * 180.02562 * 4 2 1 39 39 H 1.09005 * 109.58455 * 353.36861 * 10 9 7 40 40 H 1.09000 * 109.58597 * 113.79851 * 10 9 7 41 41 H 1.09003 * 109.50127 * 65.31462 * 11 10 9 42 42 H 1.09005 * 109.49673 * 185.37831 * 11 10 9 43 43 H 1.08995 * 109.46408 * 181.38847 * 12 11 10 44 44 H 1.08997 * 109.46010 * 301.36332 * 12 11 10 45 45 H 1.08998 * 109.50061 * 178.68602 * 13 12 11 46 46 H 1.09000 * 109.50122 * 58.58394 * 13 12 11 47 47 H 1.09001 * 110.71489 * 88.06156 * 18 16 14 48 48 H 1.09003 * 110.71974 * 324.79305 * 18 16 14 49 49 H 1.08998 * 109.88149 * 262.77474 * 19 18 16 50 50 H 1.08995 * 110.01014 * 141.68490 * 19 18 16 51 51 H 0.97005 * 120.00212 * 180.02562 * 24 23 21 52 52 H 1.08996 * 110.36683 * 94.48518 * 29 20 19 53 53 H 1.09001 * 110.36614 * 216.80266 * 29 20 19 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.5299 0.0000 0.0000 3 6 2.0400 1.4425 0.0000 4 6 2.0400 -0.7213 -1.2492 5 6 2.0206 -0.6939 1.2019 6 6 2.4119 -1.2473 2.1605 7 6 2.8892 -1.9223 3.3298 8 8 2.3173 -2.9166 3.7339 9 7 3.9721 -1.4557 3.9824 10 6 4.5937 -0.1841 3.5879 11 6 6.0851 -0.4228 3.3322 12 6 6.7128 -1.0725 4.5675 13 6 6.0454 -2.4245 4.8296 14 6 4.5534 -2.2086 5.1021 15 1 4.4265 -1.6450 6.0264 16 6 3.8535 -3.5634 5.2275 17 1 2.7761 -3.4360 5.3322 18 6 4.4343 -4.3821 6.4107 19 6 4.1633 -5.8412 5.9873 20 7 4.0060 -5.8496 4.5256 21 6 3.7335 -6.9507 3.7585 22 1 3.6369 -6.8522 2.6875 23 6 3.5872 -8.1946 4.3571 24 7 3.7051 -8.3147 5.6625 25 14 3.2232 -9.6980 3.3099 26 1 3.1390 -9.2981 1.8822 27 1 4.3096 -10.6958 3.4818 28 1 1.9312 -10.2951 3.7334 29 6 4.1966 -4.4772 4.0218 30 1 -0.3634 0.5138 0.8899 31 1 -0.3633 0.5138 -0.8901 32 1 -0.3634 -1.0276 -0.0005 33 1 3.1300 1.4425 0.0005 34 1 1.6767 1.9563 -0.8900 35 1 1.6766 1.9564 0.8900 36 1 1.6767 -1.7490 -1.2492 37 1 1.6768 -0.2074 -2.1392 38 1 3.1300 -0.7216 -1.2489 39 1 4.1217 0.1864 2.6779 40 1 4.4723 0.5457 4.3883 41 1 6.2060 -1.0826 2.4730 42 1 6.5767 0.5292 3.1319 43 1 7.7791 -1.2209 4.3975 44 1 6.5696 -0.4236 5.4314 45 1 6.5087 -2.8973 5.6955 46 1 6.1665 -3.0655 3.9563 47 1 3.9107 -4.1456 7.3370 48 1 5.5045 -4.2032 6.5149 49 1 3.2497 -6.2009 6.4606 50 1 5.0029 -6.4747 6.2735 51 1 3.6033 -9.1839 6.0808 52 1 5.2305 -4.3273 3.7112 53 1 3.5180 -4.2802 3.1918 RHF calculation, no. of doubly occupied orbitals= 64 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 ZINC001716706465.mol2 53 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 16:58:12 Heat of formation + Delta-G solvation = 57.408627 kcal Electronic energy + Delta-G solvation = -30147.980158 eV Core-core repulsion = 26408.419427 eV Total energy + Delta-G solvation = -3739.560731 eV No. of doubly occupied orbitals = 64 Molecular weight (most abundant/longest-lived isotopes) = 332.221 amu Computer time = 1.47 seconds Orbital eigenvalues (eV) -40.68002 -37.83427 -37.37221 -36.75621 -33.35430 -32.78415 -31.91394 -30.81171 -28.85504 -27.64977 -27.60281 -27.21011 -25.55156 -24.28622 -23.90940 -22.50023 -21.10895 -20.39928 -19.68020 -19.05861 -17.39845 -16.70512 -16.05423 -15.80859 -15.26570 -15.08437 -14.60255 -14.58993 -14.47087 -14.37979 -14.15105 -13.81733 -13.57139 -13.16044 -12.95606 -12.77299 -12.72899 -12.31288 -12.05513 -12.02827 -11.98305 -11.74525 -11.61294 -11.55933 -11.50632 -11.20425 -11.01320 -10.85602 -10.76351 -10.69504 -10.48122 -10.31911 -10.29590 -9.93781 -9.88416 -9.75629 -9.55531 -9.43345 -8.75615 -8.40457 -8.34053 -6.54903 -5.51976 -2.82109 1.69131 2.31390 3.01164 3.50686 3.55716 3.58618 3.65009 3.90088 4.59949 4.60801 4.61791 4.67213 4.68389 4.81261 4.84743 5.05521 5.17636 5.25184 5.27044 5.29456 5.33066 5.41418 5.44095 5.50925 5.53456 5.55932 5.67107 5.69008 5.70123 5.75084 5.82354 5.83506 5.90442 5.94398 5.95375 5.95899 6.05464 6.12831 6.20469 6.37983 6.45219 6.56588 6.74522 6.81490 6.94709 7.10602 7.25292 7.40748 7.51707 7.76604 7.85209 8.02325 8.36792 8.47912 9.43515 10.04937 10.73370 11.60272 Molecular weight = 332.22amu Principal moments of inertia in cm(-1) A = 0.010329 B = 0.003785 C = 0.003009 Principal moments of inertia in units of 10**(-40)*gram-cm**2 A = 2710.078989 B = 7396.467540 C = 9301.921625 Net atomic charges, atomic populations, and dipole contributions calculated with CM2 Atom NO. Type Charge No. of electrons 1 C -0.142 4.142 2 C 0.040 3.960 3 C -0.140 4.140 4 C -0.142 4.142 5 C -0.093 4.093 6 C -0.203 4.203 7 C 0.616 3.384 8 O -0.479 6.479 9 N -0.599 5.599 10 C 0.110 3.890 11 C -0.135 4.135 12 C -0.119 4.119 13 C -0.141 4.141 14 C 0.164 3.836 15 H 0.071 0.929 16 C -0.108 4.108 17 H 0.116 0.884 18 C -0.147 4.147 19 C 0.182 3.818 20 N -0.603 5.603 21 C -0.214 4.214 22 H 0.074 0.926 23 C -0.014 4.014 24 N -0.925 5.925 25 Si 0.908 3.092 26 H -0.281 1.281 27 H -0.292 1.292 28 H -0.292 1.292 29 C 0.143 3.857 30 H 0.061 0.939 31 H 0.068 0.932 32 H 0.063 0.937 33 H 0.061 0.939 34 H 0.068 0.932 35 H 0.060 0.940 36 H 0.064 0.936 37 H 0.068 0.932 38 H 0.061 0.939 39 H 0.081 0.919 40 H 0.068 0.932 41 H 0.068 0.932 42 H 0.071 0.929 43 H 0.064 0.936 44 H 0.061 0.939 45 H 0.080 0.920 46 H 0.082 0.918 47 H 0.055 0.945 48 H 0.054 0.946 49 H 0.042 0.958 50 H 0.038 0.962 51 H 0.248 0.752 52 H 0.011 0.989 53 H 0.050 0.950 For the charges calculated by CM2: Dipole moment (debyes) X Y Z Total from point charges 1.529 24.955 -8.029 26.260 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.200 4.200 2 C 0.040 3.960 3 C -0.197 4.197 4 C -0.199 4.199 5 C -0.095 4.095 6 C -0.206 4.206 7 C 0.404 3.596 8 O -0.354 6.354 9 N -0.333 5.333 10 C -0.013 4.013 11 C -0.173 4.173 12 C -0.156 4.156 13 C -0.179 4.179 14 C 0.062 3.938 15 H 0.089 0.911 16 C -0.127 4.127 17 H 0.133 0.867 18 C -0.186 4.186 19 C 0.057 3.943 20 N -0.329 5.329 21 C -0.343 4.343 22 H 0.092 0.908 23 C -0.209 4.209 24 N -0.574 5.574 25 Si 0.739 3.261 26 H -0.207 1.207 27 H -0.220 1.220 28 H -0.219 1.219 29 C 0.017 3.983 30 H 0.080 0.920 31 H 0.087 0.913 32 H 0.083 0.917 33 H 0.080 0.920 34 H 0.087 0.913 35 H 0.079 0.921 36 H 0.083 0.917 37 H 0.087 0.913 38 H 0.080 0.920 39 H 0.099 0.901 40 H 0.087 0.913 41 H 0.086 0.914 42 H 0.090 0.910 43 H 0.083 0.917 44 H 0.080 0.920 45 H 0.098 0.902 46 H 0.100 0.900 47 H 0.074 0.926 48 H 0.073 0.927 49 H 0.060 0.940 50 H 0.056 0.944 51 H 0.059 0.941 52 H 0.028 0.972 53 H 0.068 0.932 Dipole moment (debyes) X Y Z Total from point charges 1.185 24.173 -7.769 25.418 hybrid contribution -0.065 -0.850 -0.474 0.975 sum 1.120 23.323 -8.243 24.762 Atomic orbital electron populations 1.21907 0.92450 1.03079 1.02520 1.17870 0.95926 0.94146 0.88091 1.21897 1.01627 0.93571 1.02630 1.21903 1.01725 1.00914 0.95388 1.22424 0.97165 0.96332 0.93593 1.21046 1.02411 1.00235 0.96906 1.14882 0.80977 0.83030 0.80678 1.90971 1.56505 1.23402 1.64528 1.47894 1.28746 1.22680 1.33996 1.21694 0.93586 0.85324 1.00675 1.21921 0.94942 1.01239 0.99204 1.21532 0.99649 0.96951 0.97501 1.22197 0.94631 0.97156 1.03904 1.21284 0.92573 0.93026 0.86927 0.91094 1.23077 1.01822 0.90677 0.97099 0.86672 1.23078 0.99148 0.99345 0.97046 1.21051 1.00019 0.89057 0.84181 1.44441 1.83196 1.02468 1.02752 1.19108 1.39816 0.85193 0.90164 0.90795 1.24755 1.02869 0.97403 0.95907 1.69988 1.57254 1.22420 1.07728 0.86082 0.77456 0.83323 0.79225 1.20673 1.21969 1.21916 1.21997 0.96116 0.85152 0.95025 0.92030 0.91342 0.91740 0.91994 0.91301 0.92050 0.91736 0.91287 0.91961 0.90091 0.91330 0.91359 0.90996 0.91744 0.92030 0.90166 0.89994 0.92644 0.92750 0.94002 0.94379 0.94124 0.97209 0.93223 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.14 -1.12 8.43 37.15 0.31 -0.81 16 2 C 0.04 0.34 0.95 -157.30 -0.15 0.19 16 3 C -0.14 -0.89 8.43 37.16 0.31 -0.58 16 4 C -0.14 -1.13 8.43 37.16 0.31 -0.81 16 5 C -0.09 -1.15 11.71 -37.12 -0.43 -1.58 16 6 C -0.20 -3.03 11.58 -38.83 -0.45 -3.48 16 7 C 0.62 10.56 6.71 -14.34 -0.10 10.46 16 8 O -0.48 -10.27 10.75 5.27 0.06 -10.21 16 9 N -0.60 -8.53 2.96 -168.07 -0.50 -9.03 16 10 C 0.11 1.13 6.34 -3.71 -0.02 1.10 16 11 C -0.13 -1.25 6.08 -26.62 -0.16 -1.41 16 12 C -0.12 -1.19 5.88 -26.69 -0.16 -1.35 16 13 C -0.14 -1.91 4.72 -26.61 -0.13 -2.03 16 14 C 0.16 2.47 2.80 -67.60 -0.19 2.28 16 15 H 0.07 0.96 8.14 -51.93 -0.42 0.54 16 16 C -0.11 -2.14 2.13 -88.35 -0.19 -2.32 16 17 H 0.12 2.58 5.74 -52.89 -0.30 2.28 16 18 C -0.15 -2.99 6.32 -24.89 -0.16 -3.15 16 19 C 0.18 4.82 5.37 -3.07 -0.02 4.81 16 20 N -0.60 -17.04 3.39 -169.94 -0.58 -17.61 16 21 C -0.21 -6.53 10.28 -15.06 -0.15 -6.69 16 22 H 0.07 2.22 7.83 -52.49 -0.41 1.80 16 23 C -0.01 -0.43 5.49 -17.43 -0.10 -0.52 16 24 N -0.92 -29.00 10.31 55.49 0.57 -28.42 16 25 Si 0.91 23.66 29.55 -169.99 -5.02 18.64 16 26 H -0.28 -7.13 7.11 56.52 0.40 -6.73 16 27 H -0.29 -7.57 7.11 56.52 0.40 -7.17 16 28 H -0.29 -7.63 7.11 56.52 0.40 -7.23 16 29 C 0.14 3.36 5.67 -2.41 -0.01 3.35 16 30 H 0.06 0.49 8.14 -51.93 -0.42 0.06 16 31 H 0.07 0.39 8.14 -51.93 -0.42 -0.03 16 32 H 0.06 0.58 8.14 -51.93 -0.42 0.15 16 33 H 0.06 0.38 8.14 -51.93 -0.42 -0.04 16 34 H 0.07 0.32 8.14 -51.93 -0.42 -0.10 16 35 H 0.06 0.39 8.14 -51.93 -0.42 -0.03 16 36 H 0.06 0.59 8.14 -51.93 -0.42 0.16 16 37 H 0.07 0.40 8.14 -51.93 -0.42 -0.02 16 38 H 0.06 0.49 8.14 -51.93 -0.42 0.07 16 39 H 0.08 0.78 6.31 -51.93 -0.33 0.45 16 40 H 0.07 0.60 8.14 -51.93 -0.42 0.17 16 41 H 0.07 0.72 8.14 -51.93 -0.42 0.30 16 42 H 0.07 0.50 8.14 -51.93 -0.42 0.08 16 43 H 0.06 0.60 8.14 -51.93 -0.42 0.17 16 44 H 0.06 0.55 8.14 -51.93 -0.42 0.13 16 45 H 0.08 1.08 6.48 -51.93 -0.34 0.75 16 46 H 0.08 1.30 5.67 -51.93 -0.29 1.01 16 47 H 0.05 1.03 8.14 -51.93 -0.42 0.60 16 48 H 0.05 0.99 6.10 -51.93 -0.32 0.68 16 49 H 0.04 1.23 7.39 -51.93 -0.38 0.85 16 50 H 0.04 1.09 7.49 -51.93 -0.39 0.71 16 51 H 0.25 7.52 8.31 -40.82 -0.34 7.18 16 52 H 0.01 0.24 5.38 -51.93 -0.28 -0.04 16 53 H 0.05 1.26 5.64 -51.93 -0.29 0.97 16 LS Contribution 398.23 15.07 6.00 6.00 Total: -1.00 -35.28 398.23 -10.18 -45.45 By element: Atomic # 1 Polarization: 6.96 SS G_CDS: -9.23 Total: -2.27 kcal Atomic # 6 Polarization: -1.07 SS G_CDS: -1.47 Total: -2.54 kcal Atomic # 7 Polarization: -54.56 SS G_CDS: -0.50 Total: -55.06 kcal Atomic # 8 Polarization: -10.27 SS G_CDS: 0.06 Total: -10.21 kcal Atomic # 14 Polarization: 23.66 SS G_CDS: -5.02 Total: 18.64 kcal Total LS contribution 6.00 Total: 6.00 kcal Total: -35.28 -10.18 -45.45 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. ZINC001716706465.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 102.862 kcal (2) G-P(sol) polarization free energy of solvation -35.278 kcal (3) G-ENP(sol) elect.-nuc.-pol. free energy of system 67.584 kcal (4) G-CDS(sol) cavity-dispersion-solvent structure free energy -10.175 kcal (5) G-P-CDS(sol) = G-P(sol) + G-CDS(sol) = (2) + (4) -45.453 kcal (6) G-S(sol) free energy of system = (1) + (5) 57.409 kcal 1SCF WAS SPECIFIED, SO GEOMETRY OPTIMIZATION WAS NOT USED Total computer time = 1.47 seconds