Wall clock time and date at job start Tue Mar 31 2020 19:41:33 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.52998 * 1 3 3 H 1.08997 * 109.46905 * 2 1 4 4 C 1.53003 * 109.47353 * 119.99712 * 2 1 3 5 5 C 1.50700 * 109.47203 * 295.00092 * 4 2 1 6 6 O 1.21273 * 119.99715 * 359.97438 * 5 4 2 7 7 N 1.34778 * 119.99464 * 179.97438 * 5 4 2 8 8 C 1.46505 * 119.99544 * 180.02562 * 7 5 4 9 9 C 1.53946 * 113.58087 * 154.99969 * 8 7 5 10 10 C 1.52905 * 86.78200 * 89.86899 * 9 8 7 11 11 C 1.52913 * 121.13639 * 150.61618 * 10 9 8 12 12 C 1.53945 * 86.77956 * 100.82368 * 11 10 9 13 13 N 1.46506 * 113.57811 * 139.34092 * 12 11 10 14 14 C 1.36931 * 119.99939 * 154.99803 * 13 12 11 15 15 H 1.08000 * 120.00644 * 359.97438 * 14 13 12 16 16 C 1.38815 * 119.99715 * 179.97438 * 14 13 12 17 17 N 1.31611 * 120.00107 * 0.02562 * 16 14 13 18 18 Si 1.86805 * 119.99746 * 179.72150 * 16 14 13 19 19 H 1.48497 * 109.47323 * 0.02562 * 18 16 14 20 20 H 1.48506 * 109.46790 * 120.00304 * 18 16 14 21 21 H 1.48501 * 109.47364 * 240.00332 * 18 16 14 22 22 C 1.52906 * 121.13836 * 259.17847 * 10 9 8 23 23 C 1.52908 * 88.04323 * 24.89970 * 10 9 8 24 24 C 1.52995 * 109.47242 * 239.99928 * 2 1 3 25 25 H 1.08997 * 109.82807 * 185.56423 * 24 2 1 26 26 C 1.54387 * 109.82909 * 64.58275 * 24 2 1 27 27 C 1.55310 * 102.25188 * 219.12816 * 26 24 2 28 28 C 1.54943 * 100.94020 * 34.13773 * 27 26 24 29 29 O 1.43246 * 109.82924 * 306.45239 * 24 2 1 30 30 H 1.09006 * 109.47000 * 180.02562 * 1 2 3 31 31 H 1.09007 * 109.46739 * 299.99953 * 1 2 3 32 32 H 1.08994 * 109.47166 * 59.99882 * 1 2 3 33 33 H 1.09002 * 109.47241 * 174.99620 * 4 2 1 34 34 H 1.09001 * 109.47047 * 55.00113 * 4 2 1 35 35 H 0.96994 * 120.00569 * 359.97438 * 7 5 4 36 36 H 1.09005 * 112.87852 * 23.85657 * 8 7 5 37 37 H 1.08999 * 113.66236 * 335.39637 * 9 8 7 38 38 H 1.08997 * 113.66296 * 204.34609 * 9 8 7 39 39 H 1.09006 * 113.65906 * 346.34780 * 11 10 9 40 40 H 1.09001 * 113.74510 * 215.31106 * 11 10 9 41 41 H 1.09004 * 113.57989 * 270.05207 * 12 11 10 42 42 H 0.97000 * 119.99307 * 334.99597 * 13 12 11 43 43 H 0.97000 * 120.00330 * 179.97438 * 17 16 14 44 44 H 1.09003 * 113.66036 * 144.70622 * 22 10 9 45 45 H 1.08997 * 113.66322 * 13.65622 * 22 10 9 46 46 H 1.08999 * 113.66416 * 89.57307 * 23 10 9 47 47 H 1.08999 * 113.66154 * 220.62741 * 23 10 9 48 48 H 1.08990 * 110.84903 * 337.32253 * 26 24 2 49 49 H 1.09006 * 110.84551 * 100.94542 * 26 24 2 50 50 H 1.08987 * 111.12840 * 276.21475 * 27 26 24 51 51 H 1.09007 * 111.15551 * 152.06773 * 27 26 24 52 52 H 1.08997 * 110.30783 * 204.63239 * 28 27 26 53 53 H 1.09008 * 110.31216 * 82.46897 * 28 27 26 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.8933 1.0276 0.0000 4 6 2.0401 -0.7212 1.2493 5 6 1.6414 0.0571 2.4766 6 8 1.0118 1.0874 2.3632 7 7 1.9846 -0.3919 3.7002 8 6 1.5966 0.3645 4.8934 9 6 1.4955 -0.4989 6.1639 10 6 2.9726 -0.1786 6.3951 11 6 3.5576 -0.1209 7.8068 12 6 4.1146 -1.5310 7.5399 13 7 5.4729 -1.7309 8.0512 14 6 5.9094 -2.9914 8.3605 15 1 5.2523 -3.8382 8.2280 16 6 7.1966 -3.1809 8.8444 17 7 7.9973 -2.1490 9.0064 18 14 7.7894 -4.8996 9.2738 19 1 6.7041 -5.8781 9.0092 20 1 8.9729 -5.2370 8.4427 21 1 8.1626 -4.9501 10.7103 22 6 4.0655 -1.1873 6.0401 23 6 2.7720 1.0930 5.5699 24 6 2.0400 -0.7212 -1.2492 25 1 3.1254 -0.8115 -1.2067 26 6 1.6269 0.0640 -2.5126 27 6 1.2539 -1.0745 -3.5010 28 6 0.6556 -2.1229 -2.5296 29 8 1.4471 -2.0224 -1.3351 30 1 -0.3633 -1.0277 0.0005 31 1 -0.3633 0.5139 0.8901 32 1 -0.3633 0.5138 -0.8899 33 1 3.1264 -0.7995 1.2056 34 1 1.6055 -1.7198 1.2952 35 1 2.4879 -1.2161 3.7908 36 1 0.7268 1.0002 4.7278 37 1 1.3027 -1.5526 5.9626 38 1 0.8335 -0.0834 6.9237 39 1 2.8025 -0.0901 8.5923 40 1 4.3298 0.6383 7.9312 41 1 3.4325 -2.3310 7.8278 42 1 6.0630 -0.9704 8.1705 43 1 8.8966 -2.2813 9.3449 44 1 4.9804 -0.7277 5.6663 45 1 3.7162 -2.0086 5.4143 46 1 2.4689 1.9560 6.1627 47 1 3.5989 1.3109 4.8941 48 1 0.7650 0.6995 -2.3098 49 1 2.4620 0.6512 -2.8947 50 1 2.1391 -1.4644 -4.0032 51 1 0.5103 -0.7397 -4.2243 52 1 0.7379 -3.1231 -2.9550 53 1 -0.3866 -1.8874 -2.3138 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 ZINC001122113278.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 19:41:33 Heat of formation + Delta-G solvation = -42.598658 kcal Electronic energy + Delta-G solvation = -29320.180256 eV Core-core repulsion = 25385.826832 eV Total energy + Delta-G solvation = -3934.353424 eV No. of doubly occupied orbitals = 65 Molecular weight (most abundant/longest-lived isotopes) = 336.220 amu Computer time = 0.86 seconds Orbital eigenvalues (eV) -40.76954 -40.16541 -39.01770 -35.90014 -34.63368 -33.48821 -32.47735 -32.17656 -29.66830 -28.91412 -27.86092 -25.74674 -25.20600 -24.37403 -23.78261 -22.45527 -21.95659 -20.80618 -19.83696 -19.74597 -17.85976 -17.47467 -17.07279 -17.00007 -16.18844 -15.94726 -15.62076 -15.12880 -14.90873 -14.68139 -14.36318 -14.21474 -14.02089 -13.96119 -13.46492 -13.13435 -13.04279 -12.62193 -12.46925 -12.33344 -12.15893 -12.12726 -11.98837 -11.72416 -11.64288 -11.58527 -11.22750 -11.20365 -10.98927 -10.96921 -10.87285 -10.81992 -10.78030 -10.20975 -10.01543 -9.91766 -9.78971 -9.75176 -9.59916 -9.32061 -8.95707 -8.84960 -6.95884 -5.78745 -3.07214 2.61031 3.11454 3.41822 3.47594 3.49791 3.50593 4.10883 4.33743 4.49134 4.55452 4.58287 4.64590 4.69369 4.75291 4.75842 4.83476 4.90039 4.96484 5.02685 5.18788 5.30199 5.34346 5.36899 5.37782 5.44884 5.46471 5.48504 5.55559 5.63526 5.68546 5.77369 5.79720 5.84892 5.85992 5.90899 5.97359 5.99261 6.05693 6.10543 6.22352 6.25669 6.29224 6.34351 6.41605 6.54405 6.67764 6.73206 6.94129 7.13945 7.51516 7.69333 8.22142 8.39997 9.51073 10.10100 10.44284 11.43752 Molecular weight = 336.22amu Principal moments of inertia in cm(-1) A = 0.013735 B = 0.002099 C = 0.001892 Principal moments of inertia in units of 10**(-40)*gram-cm**2 A = 2038.148866 B =13336.345516 C =14792.328970 Net atomic charges, atomic populations, and dipole contributions calculated with CM2 Atom NO. Type Charge No. of electrons 1 C -0.151 4.151 2 C -0.082 4.082 3 H 0.086 0.914 4 C -0.135 4.135 5 C 0.516 3.484 6 O -0.544 6.544 7 N -0.709 5.709 8 C 0.122 3.878 9 C -0.122 4.122 10 C -0.104 4.104 11 C -0.139 4.139 12 C 0.169 3.831 13 N -0.713 5.713 14 C -0.242 4.242 15 H 0.071 0.929 16 C -0.015 4.015 17 N -0.932 5.932 18 Si 0.906 3.094 19 H -0.280 1.280 20 H -0.292 1.292 21 H -0.292 1.292 22 C -0.140 4.140 23 C -0.137 4.137 24 C 0.083 3.917 25 H 0.081 0.919 26 C -0.154 4.154 27 C -0.157 4.157 28 C 0.050 3.950 29 O -0.374 6.374 30 H 0.054 0.946 31 H 0.075 0.925 32 H 0.049 0.951 33 H 0.093 0.907 34 H 0.099 0.901 35 H 0.403 0.597 36 H 0.097 0.903 37 H 0.072 0.928 38 H 0.093 0.907 39 H 0.059 0.941 40 H 0.078 0.922 41 H 0.074 0.926 42 H 0.385 0.615 43 H 0.255 0.745 44 H 0.079 0.921 45 H 0.051 0.949 46 H 0.086 0.914 47 H 0.068 0.932 48 H 0.083 0.917 49 H 0.081 0.919 50 H 0.077 0.923 51 H 0.081 0.919 52 H 0.082 0.918 53 H 0.059 0.941 For the charges calculated by CM2: Dipole moment (debyes) X Y Z Total from point charges -16.681 5.003 -22.512 28.462 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.208 4.208 2 C -0.101 4.101 3 H 0.104 0.896 4 C -0.175 4.175 5 C 0.303 3.697 6 O -0.422 6.422 7 N -0.360 5.360 8 C 0.016 3.984 9 C -0.160 4.160 10 C -0.105 4.105 11 C -0.178 4.178 12 C 0.057 3.943 13 N -0.351 5.351 14 C -0.373 4.373 15 H 0.089 0.911 16 C -0.211 4.211 17 N -0.581 5.581 18 Si 0.737 3.263 19 H -0.205 1.205 20 H -0.219 1.219 21 H -0.219 1.219 22 C -0.181 4.181 23 C -0.176 4.176 24 C 0.025 3.975 25 H 0.099 0.901 26 C -0.191 4.191 27 C -0.195 4.195 28 C -0.026 4.026 29 O -0.293 6.293 30 H 0.073 0.927 31 H 0.094 0.906 32 H 0.068 0.932 33 H 0.111 0.889 34 H 0.117 0.883 35 H 0.238 0.762 36 H 0.115 0.885 37 H 0.091 0.909 38 H 0.111 0.889 39 H 0.077 0.923 40 H 0.097 0.903 41 H 0.092 0.908 42 H 0.219 0.781 43 H 0.065 0.935 44 H 0.098 0.902 45 H 0.070 0.930 46 H 0.104 0.896 47 H 0.086 0.914 48 H 0.102 0.898 49 H 0.099 0.901 50 H 0.096 0.904 51 H 0.099 0.901 52 H 0.101 0.899 53 H 0.077 0.923 Dipole moment (debyes) X Y Z Total from point charges -17.058 5.368 -22.369 28.639 hybrid contribution 0.959 -0.885 0.184 1.318 sum -16.099 4.483 -22.185 27.775 Atomic orbital electron populations 1.21926 0.94365 1.01877 1.02612 1.21233 0.97397 0.99227 0.92262 0.89614 1.21697 1.02915 1.00691 0.92210 1.20723 0.78743 0.85566 0.84665 1.90639 1.40190 1.24703 1.86640 1.45783 1.57403 1.27576 1.05261 1.21835 0.98139 0.94324 0.84053 1.23033 0.96283 1.00226 0.96434 1.22345 0.94283 0.96291 0.97593 1.23276 0.99732 1.00250 0.94547 1.21534 0.88080 0.93767 0.90956 1.42435 1.09696 1.06937 1.76059 1.19350 0.96523 0.85723 1.35697 0.91076 1.24769 0.96829 0.97276 1.02176 1.69667 1.07812 1.28513 1.52070 0.86128 0.78466 0.83932 0.77730 1.20536 1.21875 1.21874 1.23586 0.96397 0.96466 1.01661 1.23224 0.98401 0.97296 0.98655 1.22839 0.97393 0.82923 0.94326 0.90127 1.23238 1.02780 0.98429 0.94691 1.23257 1.01745 0.96118 0.98353 1.23363 0.95012 0.99486 0.84780 1.88429 1.66144 1.31990 1.42728 0.92694 0.90640 0.93228 0.88915 0.88262 0.76152 0.88525 0.90930 0.88872 0.92261 0.90307 0.90847 0.78129 0.93537 0.90218 0.93034 0.89553 0.91361 0.89841 0.90068 0.90438 0.90076 0.89935 0.92292 Number of geometries 1 Number of calculations of the screened coulomb radii 1 The total number of SCF iterations 13. 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.15 -1.43 8.47 37.16 0.31 -1.12 16 2 C -0.08 -0.74 1.86 -90.62 -0.17 -0.91 16 3 H 0.09 0.83 7.95 -51.93 -0.41 0.42 16 4 C -0.14 -1.26 5.04 -27.88 -0.14 -1.40 16 5 C 0.52 6.13 7.55 -10.98 -0.08 6.04 16 6 O -0.54 -8.03 13.57 5.56 0.08 -7.95 16 7 N -0.71 -8.00 5.26 -52.85 -0.28 -8.28 16 8 C 0.12 1.45 4.31 -66.95 -0.29 1.16 16 9 C -0.12 -1.41 7.59 -26.29 -0.20 -1.61 16 10 C -0.10 -1.43 1.02 -154.71 -0.16 -1.59 16 11 C -0.14 -2.25 7.59 -26.28 -0.20 -2.45 16 12 C 0.17 3.26 4.48 -66.95 -0.30 2.96 16 13 N -0.71 -17.71 5.40 -52.47 -0.28 -17.99 16 14 C -0.24 -6.65 9.97 -15.06 -0.15 -6.80 16 15 H 0.07 1.90 7.83 -52.49 -0.41 1.49 16 16 C -0.01 -0.42 5.50 -17.43 -0.10 -0.51 16 17 N -0.93 -28.09 13.06 55.49 0.72 -27.36 16 18 Si 0.91 22.26 29.55 -169.99 -5.02 17.24 16 19 H -0.28 -6.64 7.11 56.52 0.40 -6.24 16 20 H -0.29 -7.20 7.11 56.52 0.40 -6.80 16 21 H -0.29 -7.22 7.11 56.52 0.40 -6.81 16 22 C -0.14 -2.28 7.18 -26.29 -0.19 -2.47 16 23 C -0.14 -1.74 7.64 -26.28 -0.20 -1.94 16 24 C 0.08 0.66 3.21 -26.01 -0.08 0.57 16 25 H 0.08 0.61 8.14 -51.93 -0.42 0.19 16 26 C -0.15 -0.96 6.36 -24.78 -0.16 -1.11 16 27 C -0.16 -0.85 7.16 -24.48 -0.18 -1.03 16 28 C 0.05 0.34 7.71 38.18 0.29 0.63 16 29 O -0.37 -3.44 9.82 -35.23 -0.35 -3.78 16 30 H 0.05 0.53 7.54 -51.93 -0.39 0.14 16 31 H 0.07 0.88 5.95 -51.92 -0.31 0.57 16 32 H 0.05 0.41 6.34 -51.93 -0.33 0.08 16 33 H 0.09 0.77 8.14 -51.93 -0.42 0.34 16 34 H 0.10 0.90 8.01 -51.93 -0.42 0.48 16 35 H 0.40 4.14 7.57 -40.82 -0.31 3.83 16 36 H 0.10 1.22 7.47 -51.93 -0.39 0.83 16 37 H 0.07 0.82 8.14 -51.93 -0.42 0.40 16 38 H 0.09 0.98 8.14 -51.93 -0.42 0.56 16 39 H 0.06 0.90 8.14 -51.93 -0.42 0.48 16 40 H 0.08 1.35 8.12 -51.93 -0.42 0.93 16 41 H 0.07 1.42 8.10 -51.93 -0.42 1.00 16 42 H 0.38 10.20 7.92 -40.82 -0.32 9.88 16 43 H 0.25 7.42 8.31 -40.82 -0.34 7.08 16 44 H 0.08 1.43 8.14 -51.93 -0.42 1.01 16 45 H 0.05 0.77 7.47 -51.93 -0.39 0.38 16 46 H 0.09 1.03 8.14 -51.93 -0.42 0.61 16 47 H 0.07 0.93 8.14 -51.93 -0.42 0.51 16 48 H 0.08 0.52 6.07 -51.93 -0.32 0.21 16 49 H 0.08 0.43 8.14 -51.93 -0.42 0.01 16 50 H 0.08 0.40 8.14 -51.94 -0.42 -0.02 16 51 H 0.08 0.34 8.14 -51.93 -0.42 -0.08 16 52 H 0.08 0.48 8.14 -51.93 -0.42 0.05 16 53 H 0.06 0.38 8.14 -51.92 -0.42 -0.05 16 LS Contribution 411.10 15.07 6.20 6.20 Total: -1.00 -31.68 411.10 -10.38 -42.06 By element: Atomic # 1 Polarization: 20.93 SS G_CDS: -9.47 Total: 11.47 kcal Atomic # 6 Polarization: -9.61 SS G_CDS: -1.98 Total: -11.59 kcal Atomic # 7 Polarization: -53.79 SS G_CDS: 0.16 Total: -53.63 kcal Atomic # 8 Polarization: -11.47 SS G_CDS: -0.27 Total: -11.74 kcal Atomic # 14 Polarization: 22.26 SS G_CDS: -5.02 Total: 17.24 kcal Total LS contribution 6.20 Total: 6.20 kcal Total: -31.68 -10.38 -42.06 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. ZINC001122113278.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 -0.540 kcal (2) G-P(sol) polarization free energy of solvation -31.677 kcal (3) G-ENP(sol) elect.-nuc.-pol. free energy of system -32.217 kcal (4) G-CDS(sol) cavity-dispersion-solvent structure free energy -10.382 kcal (5) G-P-CDS(sol) = G-P(sol) + G-CDS(sol) = (2) + (4) -42.059 kcal (6) G-S(sol) free energy of system = (1) + (5) -42.599 kcal 1SCF WAS SPECIFIED, SO GEOMETRY OPTIMIZATION WAS NOT USED Total computer time = 0.86 seconds