Wall clock time and date at job start Tue Mar 31 2020 21:01:56 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.53006 * 1 3 3 H 1.09005 * 109.47280 * 2 1 4 4 C 1.53094 * 109.47893 * 120.02700 * 2 1 3 5 5 C 1.53181 * 109.16067 * 177.59228 * 4 2 1 6 6 H 1.09004 * 109.61413 * 62.93710 * 5 4 2 7 7 C 1.50697 * 109.54048 * 183.15346 * 5 4 2 8 8 O 1.21275 * 120.00102 * 359.56305 * 7 5 4 9 9 N 1.34782 * 119.99723 * 179.55617 * 7 5 4 10 10 C 1.46500 * 119.99930 * 180.27906 * 9 7 5 11 11 C 1.52999 * 109.47216 * 179.72693 * 10 9 7 12 12 H 1.09010 * 109.43439 * 59.65915 * 11 10 9 13 13 C 1.53146 * 109.43601 * 179.52368 * 11 10 9 14 14 C 1.53695 * 108.87852 * 66.17718 * 13 11 10 15 15 H 1.09001 * 112.58805 * 187.93007 * 14 13 11 16 16 C 1.54614 * 107.83191 * 312.22368 * 14 13 11 17 17 C 1.54315 * 104.74605 * 88.66112 * 16 14 13 18 18 C 1.54605 * 104.75203 * 0.02562 * 17 16 14 19 19 H 1.09002 * 112.69304 * 147.07504 * 18 17 16 20 20 C 1.53143 * 109.43610 * 299.75046 * 11 10 9 21 21 N 1.48289 * 104.26365 * 25.13068 * 18 17 16 22 22 C 1.36937 * 128.20812 * 138.38111 * 21 18 17 23 23 H 1.08003 * 119.99821 * 157.87535 * 22 21 18 24 24 C 1.38810 * 120.00315 * 337.87909 * 22 21 18 25 25 N 1.31621 * 120.00304 * 159.24055 * 24 22 21 26 26 Si 1.86796 * 120.00073 * 339.24100 * 24 22 21 27 27 H 1.48502 * 109.47214 * 92.78443 * 26 24 22 28 28 H 1.48500 * 109.47572 * 212.78890 * 26 24 22 29 29 H 1.48501 * 109.46999 * 332.78528 * 26 24 22 30 30 C 1.53177 * 109.03840 * 303.02732 * 5 4 2 31 31 C 1.53083 * 109.16372 * 56.97863 * 30 5 4 32 32 O 1.42894 * 109.48416 * 239.97686 * 2 1 3 33 33 H 1.09001 * 109.46380 * 180.02562 * 1 2 3 34 34 H 1.08994 * 109.47355 * 300.02576 * 1 2 3 35 35 H 1.08996 * 109.46926 * 60.02861 * 1 2 3 36 36 H 1.08999 * 109.51602 * 57.69436 * 4 2 1 37 37 H 1.08999 * 109.52597 * 297.50070 * 4 2 1 38 38 H 0.97002 * 120.00301 * 0.02562 * 9 7 5 39 39 H 1.09003 * 109.47162 * 299.72488 * 10 9 7 40 40 H 1.08994 * 109.47474 * 59.72303 * 10 9 7 41 41 H 1.09003 * 109.56580 * 185.99551 * 13 11 10 42 42 H 1.08997 * 109.57128 * 306.35473 * 13 11 10 43 43 H 1.09002 * 110.38519 * 207.48107 * 16 14 13 44 44 H 1.08994 * 110.37997 * 329.84720 * 16 14 13 45 45 H 1.08996 * 110.37952 * 241.18340 * 17 16 14 46 46 H 1.09008 * 110.37498 * 118.82602 * 17 16 14 47 47 H 1.09002 * 109.57354 * 53.63641 * 20 11 10 48 48 H 1.08997 * 109.57452 * 173.88489 * 20 11 10 49 49 H 0.96994 * 120.00249 * 180.02562 * 25 24 22 50 50 H 1.08998 * 109.52256 * 176.88534 * 30 5 4 51 51 H 1.09003 * 109.51890 * 297.00833 * 30 5 4 52 52 H 1.08995 * 109.48503 * 62.37522 * 31 30 5 53 53 H 1.09001 * 109.48226 * 182.39811 * 31 30 5 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.5301 0.0000 0.0000 3 1 1.8934 1.0277 0.0000 4 6 2.0406 -0.7222 1.2496 5 6 3.5711 -0.7708 1.2122 6 1 3.9703 0.2431 1.2417 7 6 4.0800 -1.5475 2.3991 8 8 3.2973 -2.0101 3.2018 9 7 5.4048 -1.7282 2.5685 10 6 5.8992 -2.4885 3.7190 11 6 7.4272 -2.5501 3.6712 12 1 7.8320 -1.5382 3.6964 13 6 7.9437 -3.3349 4.8806 14 6 7.4918 -4.7984 4.7537 15 1 7.7077 -5.3823 5.6484 16 6 5.9856 -4.8002 4.4047 17 6 5.9395 -4.7438 2.8632 18 6 7.4222 -4.7138 2.4259 19 1 7.5913 -5.2243 1.4778 20 6 7.8690 -3.2435 2.3791 21 7 8.1469 -5.3592 3.5471 22 6 9.1769 -6.2594 3.4834 23 1 9.8409 -6.3807 4.3265 24 6 9.3676 -7.0158 2.3353 25 7 10.5410 -7.5609 2.0936 26 14 7.9590 -7.2602 1.1331 27 1 8.0176 -6.2126 0.0822 28 1 8.0725 -8.6002 0.5032 29 1 6.6668 -7.1621 1.8582 30 6 4.0177 -1.4599 -0.0809 31 6 3.4308 -0.7091 -1.2789 32 8 2.0067 -0.6740 -1.1664 33 1 -0.3632 -1.0277 0.0005 34 1 -0.3634 0.5142 0.8897 35 1 -0.3633 0.5134 -0.8902 36 1 1.6444 -1.7374 1.2708 37 1 1.7148 -0.1843 2.1399 38 1 6.0308 -1.3580 1.9266 39 1 5.5841 -1.9990 4.6406 40 1 5.4937 -3.4997 3.6878 41 1 9.0325 -3.2899 4.9070 42 1 7.5366 -2.9050 5.7957 43 1 5.5134 -5.7140 4.7654 44 1 5.4960 -3.9249 4.8313 45 1 5.4272 -3.8407 2.5316 46 1 5.4445 -5.6299 2.4656 47 1 7.4088 -2.7491 1.5236 48 1 8.9541 -3.1926 2.2901 49 1 10.6744 -8.0891 1.2911 50 1 5.1059 -1.4495 -0.1421 51 1 3.6624 -2.4904 -0.0859 52 1 3.8193 0.3092 -1.2957 53 1 3.7103 -1.2198 -2.2004 RHF calculation, no. of doubly occupied orbitals= 65 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 ZINC001048577523.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 21:01:56 Heat of formation + Delta-G solvation = 29.739161 kcal Electronic energy + Delta-G solvation = -31813.884386 eV Core-core repulsion = 27882.667766 eV Total energy + Delta-G solvation = -3931.216620 eV No. of doubly occupied orbitals = 65 Molecular weight (most abundant/longest-lived isotopes) = 336.220 amu Computer time = 1.41 seconds Orbital eigenvalues (eV) -41.64039 -40.21241 -38.60108 -34.96555 -33.98469 -33.15092 -32.73502 -31.71177 -30.79917 -30.61042 -28.23706 -27.08461 -24.98477 -23.70670 -23.36335 -22.89659 -21.15863 -20.98134 -20.13247 -20.03584 -18.80325 -17.33068 -17.17689 -16.65496 -16.33497 -16.14090 -15.59182 -15.32245 -14.90411 -14.72477 -14.35327 -13.77758 -13.59565 -13.47165 -13.26900 -13.14316 -12.87516 -12.82472 -12.61877 -12.41037 -12.24759 -12.13068 -11.91104 -11.69786 -11.63900 -11.56909 -11.43184 -11.23100 -11.04391 -10.75719 -10.71471 -10.51968 -10.32483 -10.08623 -10.02373 -9.77107 -9.75679 -9.55324 -9.48868 -9.28550 -8.64178 -8.32391 -7.09304 -5.97580 -3.28248 2.63135 3.20796 3.26857 3.45663 3.48259 3.48651 4.31441 4.44359 4.48119 4.50684 4.60879 4.77509 4.83449 4.88884 5.04515 5.06894 5.09967 5.15053 5.19634 5.24654 5.31216 5.33639 5.36988 5.46393 5.48804 5.57994 5.59724 5.61645 5.69062 5.72391 5.76792 5.81139 5.87586 5.90800 5.97724 6.04015 6.08781 6.18817 6.22255 6.31176 6.38855 6.50126 6.58744 6.61008 6.67425 6.75741 6.81595 6.94612 6.95611 7.54651 7.65381 8.25379 8.30384 8.77649 9.64657 10.24292 11.19040 Molecular weight = 336.22amu Principal moments of inertia in cm(-1) A = 0.016739 B = 0.003383 C = 0.003061 Principal moments of inertia in units of 10**(-40)*gram-cm**2 A = 1672.295862 B = 8274.453404 C = 9146.241812 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.086 3.914 3 H 0.057 0.943 4 C -0.129 4.129 5 C -0.110 4.110 6 H 0.095 0.905 7 C 0.513 3.487 8 O -0.542 6.542 9 N -0.722 5.722 10 C 0.001 3.999 11 C -0.082 4.082 12 H 0.058 0.942 13 C -0.121 4.121 14 C 0.153 3.847 15 H 0.078 0.922 16 C -0.163 4.163 17 C -0.169 4.169 18 C 0.140 3.860 19 H 0.073 0.927 20 C -0.124 4.124 21 N -0.596 5.596 22 C -0.233 4.233 23 H 0.092 0.908 24 C 0.022 3.978 25 N -0.910 5.910 26 Si 0.877 3.123 27 H -0.289 1.289 28 H -0.304 1.304 29 H -0.264 1.264 30 C -0.150 4.150 31 C 0.065 3.935 32 O -0.380 6.380 33 H 0.065 0.935 34 H 0.071 0.929 35 H 0.062 0.938 36 H 0.081 0.919 37 H 0.079 0.921 38 H 0.394 0.606 39 H 0.051 0.949 40 H 0.193 0.807 41 H 0.070 0.930 42 H 0.054 0.946 43 H 0.061 0.939 44 H 0.094 0.906 45 H 0.099 0.901 46 H 0.058 0.942 47 H 0.046 0.954 48 H 0.070 0.930 49 H 0.263 0.737 50 H 0.080 0.920 51 H 0.080 0.920 52 H 0.051 0.949 53 H 0.098 0.902 For the charges calculated by CM2: Dipole moment (debyes) X Y Z Total from point charges -14.825 17.859 -3.039 23.409 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.202 4.202 2 C 0.028 3.972 3 H 0.075 0.925 4 C -0.167 4.167 5 C -0.132 4.132 6 H 0.113 0.887 7 C 0.299 3.701 8 O -0.419 6.419 9 N -0.373 5.373 10 C -0.118 4.118 11 C -0.101 4.101 12 H 0.077 0.923 13 C -0.161 4.161 14 C 0.045 3.955 15 H 0.096 0.904 16 C -0.201 4.201 17 C -0.207 4.207 18 C 0.034 3.966 19 H 0.090 0.910 20 C -0.164 4.164 21 N -0.318 5.318 22 C -0.360 4.360 23 H 0.109 0.891 24 C -0.176 4.176 25 N -0.553 5.553 26 Si 0.709 3.291 27 H -0.217 1.217 28 H -0.233 1.233 29 H -0.190 1.190 30 C -0.188 4.188 31 C -0.012 4.012 32 O -0.299 6.299 33 H 0.084 0.916 34 H 0.090 0.910 35 H 0.081 0.919 36 H 0.100 0.900 37 H 0.098 0.902 38 H 0.228 0.772 39 H 0.070 0.930 40 H 0.208 0.792 41 H 0.089 0.911 42 H 0.073 0.927 43 H 0.080 0.920 44 H 0.112 0.888 45 H 0.117 0.883 46 H 0.077 0.923 47 H 0.065 0.935 48 H 0.089 0.911 49 H 0.073 0.927 50 H 0.099 0.901 51 H 0.099 0.901 52 H 0.070 0.930 53 H 0.116 0.884 Dipole moment (debyes) X Y Z Total from point charges -15.260 17.206 -2.659 23.152 hybrid contribution 0.099 0.267 -0.902 0.946 sum -15.161 17.473 -3.562 23.406 Atomic orbital electron populations 1.21755 0.92878 1.03194 1.02360 1.22136 0.95307 0.94803 0.84958 0.92500 1.22051 0.93041 1.03236 0.98335 1.21395 0.97773 1.00179 0.93823 0.88736 1.20439 0.82448 0.80565 0.86684 1.90688 1.55674 1.50037 1.45486 1.46347 1.06793 1.55761 1.28404 1.25183 0.96084 1.03589 0.86961 1.21178 0.93325 0.98577 0.97007 0.92332 1.22100 1.01381 0.97084 0.95500 1.21619 0.89654 0.92278 0.91936 0.90365 1.23383 0.99437 1.01746 0.95545 1.23525 0.99160 1.02106 0.95948 1.21733 0.90095 0.91562 0.93251 0.90997 1.22165 1.01642 0.97386 0.95170 1.45475 1.38570 1.48264 0.99478 1.18628 1.01760 1.12686 1.02904 0.89074 1.25036 0.96610 0.98059 0.97918 1.69810 1.20934 1.37334 1.27256 0.86960 0.85161 0.75910 0.81101 1.21747 1.23253 1.18950 1.22233 1.01661 1.01272 0.93681 1.22541 0.81045 0.99455 0.98142 1.87947 1.20891 1.78033 1.43023 0.91615 0.90994 0.91862 0.89983 0.90239 0.77165 0.93036 0.79226 0.91120 0.92740 0.92007 0.88794 0.88320 0.92309 0.93520 0.91116 0.92735 0.90132 0.90088 0.93027 0.88379 Number of geometries 1 Number of calculations of the screened coulomb radii 1 The total number of SCF iterations 17. 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.57 37.16 0.36 -0.72 16 2 C 0.09 0.74 3.26 -26.68 -0.09 0.65 16 3 H 0.06 0.43 8.14 -51.93 -0.42 0.00 16 4 C -0.13 -1.27 4.01 -26.59 -0.11 -1.38 16 5 C -0.11 -1.06 2.61 -91.59 -0.24 -1.30 16 6 H 0.09 0.73 8.14 -51.93 -0.42 0.30 16 7 C 0.51 6.57 6.91 -10.98 -0.08 6.50 16 8 O -0.54 -8.93 15.97 5.56 0.09 -8.84 16 9 N -0.72 -8.66 4.30 -60.30 -0.26 -8.92 16 10 C 0.00 0.01 2.68 -4.04 -0.01 0.00 16 11 C -0.08 -1.19 2.76 -90.47 -0.25 -1.44 16 12 H 0.06 0.76 8.14 -51.92 -0.42 0.34 16 13 C -0.12 -1.91 5.37 -26.29 -0.14 -2.05 16 14 C 0.15 2.83 4.32 -65.73 -0.28 2.55 16 15 H 0.08 1.45 8.14 -51.93 -0.42 1.03 16 16 C -0.16 -2.75 5.72 -25.19 -0.14 -2.89 16 17 C -0.17 -3.06 5.55 -25.20 -0.14 -3.20 16 18 C 0.14 2.87 2.88 -65.71 -0.19 2.68 16 19 H 0.07 1.65 3.12 -51.93 -0.16 1.49 16 20 C -0.12 -2.08 5.28 -26.27 -0.14 -2.22 16 21 N -0.60 -13.60 2.57 -152.73 -0.39 -13.99 16 22 C -0.23 -6.38 10.74 -15.06 -0.16 -6.54 16 23 H 0.09 2.66 8.06 -52.48 -0.42 2.24 16 24 C 0.02 0.64 5.47 -17.43 -0.10 0.54 16 25 N -0.91 -27.71 13.97 55.49 0.78 -26.94 16 26 Si 0.88 22.18 24.71 -169.99 -4.20 17.98 16 27 H -0.29 -7.09 6.67 56.52 0.38 -6.71 16 28 H -0.30 -7.90 7.11 56.52 0.40 -7.50 16 29 H -0.26 -6.49 5.62 56.52 0.32 -6.17 16 30 C -0.15 -1.34 5.29 -26.60 -0.14 -1.48 16 31 C 0.06 0.53 6.84 37.20 0.25 0.78 16 32 O -0.38 -3.89 10.41 -35.23 -0.37 -4.25 16 33 H 0.06 0.53 8.14 -51.93 -0.42 0.11 16 34 H 0.07 0.44 8.14 -51.93 -0.42 0.01 16 35 H 0.06 0.44 8.14 -51.93 -0.42 0.01 16 36 H 0.08 0.97 8.03 -51.93 -0.42 0.55 16 37 H 0.08 0.77 8.14 -51.93 -0.42 0.35 16 38 H 0.39 3.91 6.82 -40.82 -0.28 3.63 16 39 H 0.05 0.73 7.99 -51.93 -0.42 0.32 16 40 H 0.19 2.95 1.12 -165.90 -0.19 2.76 16 41 H 0.07 1.21 8.14 -51.93 -0.42 0.79 16 42 H 0.05 0.74 8.14 -51.93 -0.42 0.31 16 43 H 0.06 1.05 8.14 -51.93 -0.42 0.62 16 44 H 0.09 1.42 4.84 -51.93 -0.25 1.16 16 45 H 0.10 1.60 3.95 -51.93 -0.21 1.40 16 46 H 0.06 1.14 6.99 -51.92 -0.36 0.78 16 47 H 0.05 0.66 6.73 -51.93 -0.35 0.31 16 48 H 0.07 1.32 8.14 -51.93 -0.42 0.89 16 49 H 0.26 7.63 8.31 -40.82 -0.34 7.29 16 50 H 0.08 0.61 7.67 -51.93 -0.40 0.21 16 51 H 0.08 0.86 8.14 -51.93 -0.42 0.44 16 52 H 0.05 0.34 8.14 -51.93 -0.42 -0.08 16 53 H 0.10 0.70 8.14 -51.93 -0.42 0.28 16 LS Contribution 376.35 15.07 5.67 5.67 Total: -1.00 -32.32 376.35 -9.31 -41.63 By element: Atomic # 1 Polarization: 16.21 SS G_CDS: -9.03 Total: 7.18 kcal Atomic # 6 Polarization: -7.92 SS G_CDS: -1.59 Total: -9.51 kcal Atomic # 7 Polarization: -49.97 SS G_CDS: 0.12 Total: -49.84 kcal Atomic # 8 Polarization: -12.82 SS G_CDS: -0.28 Total: -13.10 kcal Atomic # 14 Polarization: 22.18 SS G_CDS: -4.20 Total: 17.98 kcal Total LS contribution 5.67 Total: 5.67 kcal Total: -32.32 -9.31 -41.63 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. ZINC001048577523.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 71.368 kcal (2) G-P(sol) polarization free energy of solvation -32.318 kcal (3) G-ENP(sol) elect.-nuc.-pol. free energy of system 39.050 kcal (4) G-CDS(sol) cavity-dispersion-solvent structure free energy -9.311 kcal (5) G-P-CDS(sol) = G-P(sol) + G-CDS(sol) = (2) + (4) -41.629 kcal (6) G-S(sol) free energy of system = (1) + (5) 29.739 kcal 1SCF WAS SPECIFIED, SO GEOMETRY OPTIMIZATION WAS NOT USED Total computer time = 1.42 seconds