Wall clock time and date at job start Tue Mar 31 2020 20:39:55 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.53003 * 1 3 3 H 1.09007 * 109.47247 * 2 1 4 4 N 1.46506 * 109.46985 * 119.99216 * 2 1 3 5 5 C 1.34774 * 120.00110 * 275.00504 * 4 2 1 6 6 O 1.21277 * 120.00259 * 359.97438 * 5 4 2 7 7 C 1.50704 * 119.99842 * 180.02562 * 5 4 2 8 8 C 1.53001 * 109.46786 * 299.99794 * 7 5 4 9 9 C 1.52999 * 109.47216 * 180.02562 * 8 7 5 10 10 C 1.52996 * 109.47389 * 59.99504 * 9 8 7 11 11 C 1.53006 * 109.47034 * 300.00414 * 10 9 8 12 12 C 1.52993 * 109.47228 * 59.99863 * 11 10 9 13 13 C 1.50699 * 109.47368 * 239.99382 * 2 1 3 14 14 O 1.21277 * 120.00243 * 264.93556 * 13 2 1 15 15 N 1.34775 * 119.99472 * 84.93431 * 13 2 1 16 16 C 1.46922 * 120.63457 * 354.81584 * 15 13 2 17 17 C 1.53192 * 108.77739 * 126.38782 * 16 15 13 18 18 C 1.53039 * 109.31499 * 54.63470 * 17 16 15 19 19 C 1.53041 * 109.53580 * 298.63537 * 18 17 16 20 20 H 1.09002 * 109.49948 * 301.38187 * 19 18 17 21 21 C 1.50704 * 109.49775 * 181.31628 * 19 18 17 22 22 H 1.08005 * 119.99962 * 240.06259 * 21 19 18 23 23 C 1.37875 * 119.99753 * 59.78845 * 21 19 18 24 24 N 1.31508 * 120.00257 * 0.27897 * 23 21 19 25 25 Si 1.86804 * 119.99705 * 180.27656 * 23 21 19 26 26 H 1.48496 * 109.47283 * 60.00181 * 25 23 21 27 27 H 1.48498 * 109.46911 * 180.02562 * 25 23 21 28 28 H 1.48502 * 109.47145 * 299.99785 * 25 23 21 29 29 C 1.46935 * 120.59929 * 174.78767 * 15 13 2 30 30 H 1.09005 * 109.47207 * 305.00815 * 1 2 3 31 31 H 1.08999 * 109.47102 * 65.00487 * 1 2 3 32 32 H 1.08994 * 109.47120 * 185.00557 * 1 2 3 33 33 H 0.97000 * 119.99549 * 95.00265 * 4 2 1 34 34 H 1.08997 * 109.47487 * 59.99098 * 7 5 4 35 35 H 1.09008 * 109.47199 * 300.00005 * 8 7 5 36 36 H 1.09001 * 109.47408 * 60.00024 * 8 7 5 37 37 H 1.09011 * 109.46698 * 299.99296 * 9 8 7 38 38 H 1.09000 * 109.47538 * 179.97438 * 9 8 7 39 39 H 1.08997 * 109.47165 * 180.02562 * 10 9 8 40 40 H 1.08995 * 109.47015 * 60.00369 * 10 9 8 41 41 H 1.09005 * 109.46672 * 299.99693 * 11 10 9 42 42 H 1.08994 * 109.46789 * 180.02562 * 11 10 9 43 43 H 1.09000 * 109.47396 * 179.97438 * 12 11 10 44 44 H 1.09003 * 109.47338 * 59.99759 * 12 11 10 45 45 H 1.08999 * 109.59058 * 6.59787 * 16 15 13 46 46 H 1.09004 * 109.58601 * 246.17395 * 16 15 13 47 47 H 1.08996 * 109.50000 * 294.68192 * 17 16 15 48 48 H 1.09004 * 109.52675 * 174.60298 * 17 16 15 49 49 H 1.08998 * 109.46014 * 178.61702 * 18 17 16 50 50 H 1.08997 * 109.46377 * 58.65089 * 18 17 16 51 51 H 0.97003 * 119.99927 * 180.02562 * 24 23 21 52 52 H 1.08999 * 109.56252 * 113.81082 * 29 15 13 53 53 H 1.08996 * 109.61701 * 353.38143 * 29 15 13 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.8934 1.0277 0.0000 4 7 2.0183 -0.6905 1.1963 5 6 2.1470 -0.0180 2.3572 6 8 1.8589 1.1588 2.4115 7 6 2.6498 -0.7282 3.5876 8 6 1.7063 -1.8825 3.9316 9 6 2.2163 -2.6031 5.1811 10 6 2.2673 -1.6197 6.3521 11 6 3.2109 -0.4655 6.0080 12 6 2.7008 0.2552 4.7586 13 6 2.0324 -0.7105 -1.2304 14 8 2.3220 -1.8869 -1.1755 15 7 2.1599 -0.0385 -2.3917 16 6 1.7075 1.3548 -2.5039 17 6 2.8643 2.2076 -3.0343 18 6 3.3774 1.6060 -4.3445 19 6 3.8991 0.1903 -4.0879 20 1 4.6960 0.2253 -3.3450 21 6 4.4338 -0.3892 -5.3722 22 1 4.0023 -1.2924 -5.7780 23 6 5.4701 0.2366 -6.0320 24 7 5.9956 1.3362 -5.5379 25 14 6.1277 -0.4775 -7.6280 26 1 5.0377 -0.5250 -8.6353 27 1 7.2329 0.3757 -8.1337 28 1 6.6359 -1.8513 -7.3833 29 6 2.7546 -0.6859 -3.5690 30 1 -0.3634 0.5896 0.8418 31 1 -0.3633 0.4342 -0.9314 32 1 -0.3633 -1.0237 0.0897 33 1 2.2490 -1.6316 1.1528 34 1 3.6491 -1.1201 3.3981 35 1 0.7068 -1.4906 4.1211 36 1 1.6697 -2.5829 3.0972 37 1 3.2158 -2.9950 4.9919 38 1 1.5439 -3.4252 5.4265 39 1 2.6311 -2.1331 7.2421 40 1 1.2679 -1.2280 6.5412 41 1 4.2103 -0.8574 5.8189 42 1 3.2475 0.2348 6.8424 43 1 3.3727 1.0778 4.5137 44 1 1.7013 0.6470 4.9477 45 1 1.4037 1.7206 -1.5231 46 1 0.8656 1.4112 -3.1940 47 1 3.6705 2.2236 -2.3010 48 1 2.5150 3.2246 -3.2131 49 1 4.1844 2.2248 -4.7368 50 1 2.5644 1.5674 -5.0695 51 1 6.7244 1.7767 -6.0023 52 1 1.9980 -0.7965 -4.3459 53 1 3.1419 -1.6661 -3.2911 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 ZINC000597911920.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 20:39:55 Heat of formation + Delta-G solvation = -108.377026 kcal Electronic energy + Delta-G solvation = -30311.304718 eV Core-core repulsion = 26374.098930 eV Total energy + Delta-G solvation = -3937.205788 eV No. of doubly occupied orbitals = 65 Molecular weight (most abundant/longest-lived isotopes) = 336.220 amu Computer time = 0.81 seconds Orbital eigenvalues (eV) -40.64962 -39.64651 -37.82454 -36.92202 -34.97399 -33.60368 -32.04405 -31.38071 -31.00339 -30.51350 -29.68377 -25.57690 -25.11499 -24.32187 -23.05897 -22.65810 -21.75768 -20.63718 -19.94778 -19.82869 -18.58784 -17.61950 -16.66450 -16.41944 -15.97468 -15.62890 -15.28912 -14.92688 -14.71993 -14.51521 -14.38654 -14.17024 -13.79459 -13.67981 -13.37814 -13.21484 -12.97109 -12.77432 -12.60290 -12.52776 -12.43647 -12.24540 -12.11042 -11.96571 -11.65962 -11.52123 -11.41768 -10.99553 -10.92766 -10.88845 -10.79868 -10.71229 -10.39948 -10.36780 -9.99804 -9.89801 -9.77407 -9.57743 -9.41848 -9.11882 -8.87706 -8.42833 -8.01036 -6.05430 -4.13093 2.53337 2.73556 3.03217 3.28923 3.33683 3.36032 3.52898 4.34500 4.38808 4.39323 4.44191 4.57338 4.88196 4.89159 4.94217 5.05583 5.07674 5.11289 5.17603 5.22668 5.23514 5.31678 5.33287 5.50193 5.54272 5.57329 5.66517 5.72953 5.76085 5.77784 5.81017 5.87478 5.94258 5.99754 6.03767 6.05282 6.15854 6.27846 6.29795 6.36500 6.42747 6.59279 6.73608 6.85873 7.11142 7.27533 7.45240 7.64544 7.81408 7.90568 8.02955 8.15419 8.30687 8.44519 8.90089 9.52033 11.00006 Molecular weight = 336.22amu Principal moments of inertia in cm(-1) A = 0.018888 B = 0.002382 C = 0.002326 Principal moments of inertia in units of 10**(-40)*gram-cm**2 A = 1482.085473 B =11750.851204 C =12037.057392 Net atomic charges, atomic populations, and dipole contributions calculated with CM2 Atom NO. Type Charge No. of electrons 1 C -0.155 4.155 2 C 0.122 3.878 3 H 0.110 0.890 4 N -0.712 5.712 5 C 0.520 3.480 6 O -0.542 6.542 7 C -0.105 4.105 8 C -0.111 4.111 9 C -0.118 4.118 10 C -0.122 4.122 11 C -0.120 4.120 12 C -0.098 4.098 13 C 0.508 3.492 14 O -0.550 6.550 15 N -0.595 5.595 16 C 0.103 3.897 17 C -0.136 4.136 18 C -0.113 4.113 19 C 0.040 3.960 20 H 0.029 0.971 21 C -0.516 4.516 22 H 0.057 0.943 23 C 0.064 3.936 24 N -0.888 5.888 25 Si 0.894 3.106 26 H -0.277 1.277 27 H -0.287 1.287 28 H -0.277 1.277 29 C 0.117 3.883 30 H 0.071 0.929 31 H 0.071 0.929 32 H 0.064 0.936 33 H 0.409 0.591 34 H 0.098 0.902 35 H 0.068 0.932 36 H 0.069 0.931 37 H 0.065 0.935 38 H 0.065 0.935 39 H 0.061 0.939 40 H 0.062 0.938 41 H 0.063 0.937 42 H 0.064 0.936 43 H 0.070 0.930 44 H 0.069 0.931 45 H 0.074 0.926 46 H 0.068 0.932 47 H 0.058 0.942 48 H 0.059 0.941 49 H 0.097 0.903 50 H 0.036 0.964 51 H 0.262 0.738 52 H 0.057 0.943 53 H 0.079 0.921 For the charges calculated by CM2: Dipole moment (debyes) X Y Z Total from point charges -10.044 -2.474 19.786 22.326 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.214 4.214 2 C 0.013 3.987 3 H 0.128 0.872 4 N -0.366 5.366 5 C 0.307 3.693 6 O -0.419 6.419 7 C -0.127 4.127 8 C -0.149 4.149 9 C -0.156 4.156 10 C -0.159 4.159 11 C -0.158 4.158 12 C -0.136 4.136 13 C 0.294 3.706 14 O -0.429 6.429 15 N -0.326 5.326 16 C -0.020 4.020 17 C -0.174 4.174 18 C -0.151 4.151 19 C 0.021 3.979 20 H 0.047 0.953 21 C -0.554 4.554 22 H 0.075 0.925 23 C -0.138 4.138 24 N -0.527 5.527 25 Si 0.722 3.278 26 H -0.202 1.202 27 H -0.213 1.213 28 H -0.202 1.202 29 C -0.003 4.003 30 H 0.090 0.910 31 H 0.090 0.910 32 H 0.083 0.917 33 H 0.246 0.754 34 H 0.116 0.884 35 H 0.087 0.913 36 H 0.088 0.912 37 H 0.084 0.916 38 H 0.084 0.916 39 H 0.080 0.920 40 H 0.081 0.919 41 H 0.081 0.919 42 H 0.083 0.917 43 H 0.089 0.911 44 H 0.088 0.912 45 H 0.092 0.908 46 H 0.086 0.914 47 H 0.077 0.923 48 H 0.078 0.922 49 H 0.115 0.885 50 H 0.055 0.945 51 H 0.071 0.929 52 H 0.075 0.925 53 H 0.097 0.903 Dipole moment (debyes) X Y Z Total from point charges -9.950 -2.109 19.774 22.236 hybrid contribution 0.554 -0.648 -1.094 1.387 sum -9.396 -2.757 18.680 21.091 Atomic orbital electron populations 1.22080 0.92925 1.03141 1.03269 1.20849 0.95729 0.98096 0.83980 0.87231 1.45864 1.67243 1.14869 1.08584 1.20698 0.77149 0.87171 0.84313 1.90685 1.50090 1.15972 1.85120 1.21355 1.01168 0.97121 0.93012 1.21541 0.98756 0.96080 0.98525 1.21530 1.01651 0.97844 0.94541 1.21567 1.00043 0.96497 0.97806 1.21566 0.99187 0.97450 0.97589 1.21337 1.02692 0.96955 0.92642 1.20392 0.78432 0.87204 0.84572 1.90644 1.51723 1.15708 1.84809 1.48276 1.59397 1.13325 1.11617 1.21896 0.97810 0.81236 1.01039 1.22012 0.97009 0.97760 1.00612 1.21642 1.00422 0.99987 0.93057 1.18977 0.91531 0.91140 0.96221 0.95293 1.21019 1.17160 1.10119 1.07140 0.92499 1.24936 0.96272 0.94890 0.97684 1.69858 1.25354 1.18133 1.39347 0.86630 0.78525 0.79095 0.83595 1.20229 1.21291 1.20221 1.22062 0.96458 0.97772 0.84045 0.91012 0.90973 0.91688 0.75367 0.88354 0.91341 0.91239 0.91623 0.91631 0.92013 0.91932 0.91857 0.91689 0.91148 0.91198 0.90804 0.91418 0.92325 0.92221 0.88522 0.94482 0.92856 0.92457 0.90256 Number of geometries 1 Number of calculations of the screened coulomb radii 1 The total number of SCF iterations 12. 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.16 -1.81 9.24 37.16 0.34 -1.46 16 2 C 0.12 1.71 2.11 -69.08 -0.15 1.56 16 3 H 0.11 1.61 5.93 -51.93 -0.31 1.30 16 4 N -0.71 -9.18 4.95 -54.93 -0.27 -9.45 16 5 C 0.52 6.64 7.12 -10.98 -0.08 6.56 16 6 O -0.54 -8.55 15.64 5.56 0.09 -8.47 16 7 C -0.11 -1.00 2.53 -91.77 -0.23 -1.23 16 8 C -0.11 -0.87 5.10 -26.73 -0.14 -1.01 16 9 C -0.12 -0.78 5.92 -26.73 -0.16 -0.94 16 10 C -0.12 -0.79 5.92 -26.73 -0.16 -0.95 16 11 C -0.12 -0.86 5.92 -26.73 -0.16 -1.02 16 12 C -0.10 -0.91 4.49 -26.74 -0.12 -1.03 16 13 C 0.51 8.77 6.90 -10.99 -0.08 8.69 16 14 O -0.55 -10.90 15.44 5.56 0.09 -10.82 16 15 N -0.59 -10.66 2.97 -172.76 -0.51 -11.17 16 16 C 0.10 1.57 5.81 -3.71 -0.02 1.55 16 17 C -0.14 -2.40 6.08 -26.61 -0.16 -2.56 16 18 C -0.11 -2.49 4.96 -26.69 -0.13 -2.63 16 19 C 0.04 0.98 2.25 -91.65 -0.21 0.77 16 20 H 0.03 0.79 8.14 -51.93 -0.42 0.37 16 21 C -0.52 -13.75 8.57 -34.44 -0.30 -14.04 16 22 H 0.06 1.49 7.74 -52.48 -0.41 1.09 16 23 C 0.06 1.72 5.30 -17.82 -0.09 1.63 16 24 N -0.89 -24.73 11.31 55.43 0.63 -24.10 16 25 Si 0.89 20.89 29.54 -169.99 -5.02 15.87 16 26 H -0.28 -6.37 7.11 56.52 0.40 -5.97 16 27 H -0.29 -6.68 7.11 56.52 0.40 -6.28 16 28 H -0.28 -6.46 7.11 56.52 0.40 -6.06 16 29 C 0.12 2.47 5.64 -3.70 -0.02 2.45 16 30 H 0.07 0.81 8.14 -51.93 -0.42 0.39 16 31 H 0.07 0.77 7.40 -51.93 -0.38 0.39 16 32 H 0.06 0.74 8.14 -51.93 -0.42 0.32 16 33 H 0.41 5.03 7.20 -40.82 -0.29 4.74 16 34 H 0.10 0.87 8.14 -51.93 -0.42 0.45 16 35 H 0.07 0.56 8.14 -51.92 -0.42 0.14 16 36 H 0.07 0.53 7.56 -51.93 -0.39 0.13 16 37 H 0.07 0.42 8.14 -51.92 -0.42 0.00 16 38 H 0.06 0.39 8.14 -51.93 -0.42 -0.04 16 39 H 0.06 0.35 8.14 -51.93 -0.42 -0.07 16 40 H 0.06 0.41 8.14 -51.93 -0.42 -0.01 16 41 H 0.06 0.44 8.14 -51.93 -0.42 0.02 16 42 H 0.06 0.43 8.14 -51.93 -0.42 0.01 16 43 H 0.07 0.73 8.10 -51.93 -0.42 0.31 16 44 H 0.07 0.70 8.10 -51.93 -0.42 0.28 16 45 H 0.07 0.95 5.14 -51.93 -0.27 0.69 16 46 H 0.07 0.94 8.14 -51.93 -0.42 0.52 16 47 H 0.06 1.10 8.14 -51.93 -0.42 0.68 16 48 H 0.06 0.92 8.14 -51.93 -0.42 0.50 16 49 H 0.10 2.43 6.03 -51.93 -0.31 2.12 16 50 H 0.04 0.80 8.14 -51.93 -0.42 0.37 16 51 H 0.26 6.99 8.31 -40.82 -0.34 6.65 16 52 H 0.06 1.17 8.14 -51.93 -0.42 0.75 16 53 H 0.08 1.78 7.04 -51.93 -0.37 1.42 16 LS Contribution 403.85 15.07 6.09 6.09 Total: -1.00 -30.31 403.85 -10.24 -40.55 By element: Atomic # 1 Polarization: 14.64 SS G_CDS: -9.47 Total: 5.17 kcal Atomic # 6 Polarization: -1.81 SS G_CDS: -1.85 Total: -3.67 kcal Atomic # 7 Polarization: -44.57 SS G_CDS: -0.16 Total: -44.73 kcal Atomic # 8 Polarization: -19.46 SS G_CDS: 0.17 Total: -19.28 kcal Atomic # 14 Polarization: 20.89 SS G_CDS: -5.02 Total: 15.87 kcal Total LS contribution 6.09 Total: 6.09 kcal Total: -30.31 -10.24 -40.55 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. ZINC000597911920.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 -67.827 kcal (2) G-P(sol) polarization free energy of solvation -30.306 kcal (3) G-ENP(sol) elect.-nuc.-pol. free energy of system -98.133 kcal (4) G-CDS(sol) cavity-dispersion-solvent structure free energy -10.244 kcal (5) G-P-CDS(sol) = G-P(sol) + G-CDS(sol) = (2) + (4) -40.550 kcal (6) G-S(sol) free energy of system = (1) + (5) -108.377 kcal 1SCF WAS SPECIFIED, SO GEOMETRY OPTIMIZATION WAS NOT USED Total computer time = 0.81 seconds