Wall clock time and date at job start Sat Apr 11 2020 06:11:04 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.53007 * 1 3 3 C 1.53003 * 109.46799 * 2 1 4 4 C 1.53001 * 109.47042 * 120.00083 * 2 1 3 5 5 N 1.46500 * 109.46746 * 174.99884 * 4 2 1 6 6 C 1.34782 * 119.99674 * 179.97438 * 5 4 2 7 7 O 1.21586 * 119.99888 * 5.07783 * 6 5 4 8 8 N 1.34779 * 119.99677 * 185.07410 * 6 5 4 9 9 C 1.46928 * 120.63092 * 354.96784 * 8 6 5 10 10 C 1.53194 * 108.77288 * 126.43473 * 9 8 6 11 11 C 1.53045 * 109.31360 * 54.63355 * 10 9 8 12 12 C 1.53045 * 109.53677 * 298.63473 * 11 10 9 13 13 H 1.09003 * 109.49626 * 301.38327 * 12 11 10 14 14 C 1.53000 * 109.49464 * 181.31449 * 12 11 10 15 15 C 1.50697 * 109.46710 * 175.06419 * 14 12 11 16 16 O 1.20826 * 120.00429 * 0.02562 * 15 14 12 17 17 O 1.34233 * 119.99870 * 180.02562 * 15 14 12 18 18 C 1.45202 * 117.00167 * 179.97438 * 17 15 14 19 19 H 1.09001 * 109.46954 * 30.00372 * 18 17 15 20 20 C 1.53007 * 109.47109 * 270.00270 * 18 17 15 21 21 C 1.50697 * 109.46991 * 150.00009 * 18 17 15 22 22 H 1.07995 * 120.00158 * 359.72570 * 21 18 17 23 23 C 1.37883 * 120.00043 * 179.97438 * 21 18 17 24 24 N 1.31514 * 119.99487 * 0.02562 * 23 21 18 25 25 Si 1.86791 * 120.00293 * 180.02562 * 23 21 18 26 26 H 1.48503 * 109.47111 * 359.97438 * 25 23 21 27 27 H 1.48498 * 109.47319 * 120.00172 * 25 23 21 28 28 H 1.48501 * 109.47197 * 240.00273 * 25 23 21 29 29 C 1.46934 * 120.59916 * 174.99191 * 8 6 5 30 30 H 1.08996 * 109.47066 * 300.00114 * 1 2 3 31 31 H 1.08995 * 109.46898 * 59.99879 * 1 2 3 32 32 H 1.08998 * 109.47072 * 179.97438 * 1 2 3 33 33 H 1.08995 * 109.46888 * 240.00409 * 2 1 3 34 34 H 1.08999 * 109.47005 * 299.99744 * 3 2 1 35 35 H 1.08995 * 109.47066 * 60.00518 * 3 2 1 36 36 H 1.09004 * 109.46986 * 180.02562 * 3 2 1 37 37 H 1.08993 * 109.46922 * 54.99886 * 4 2 1 38 38 H 1.08996 * 109.46722 * 294.99985 * 4 2 1 39 39 H 0.96997 * 120.00309 * 0.02562 * 5 4 2 40 40 H 1.08997 * 109.58739 * 6.64659 * 9 8 6 41 41 H 1.08998 * 109.58523 * 246.22133 * 9 8 6 42 42 H 1.09001 * 109.49720 * 294.68182 * 10 9 8 43 43 H 1.08995 * 109.49736 * 174.58845 * 10 9 8 44 44 H 1.09000 * 109.46140 * 178.61869 * 11 10 9 45 45 H 1.08997 * 109.45829 * 58.65617 * 11 10 9 46 46 H 1.08999 * 109.47132 * 55.06565 * 14 12 11 47 47 H 1.08992 * 109.47127 * 295.06421 * 14 12 11 48 48 H 1.09003 * 109.46823 * 180.02562 * 20 18 17 49 49 H 1.08996 * 109.46753 * 299.99432 * 20 18 17 50 50 H 1.08992 * 109.47207 * 59.99836 * 20 18 17 51 51 H 0.97007 * 119.99718 * 179.97438 * 24 23 21 52 52 H 1.08996 * 109.55901 * 113.76794 * 29 8 6 53 53 H 1.08992 * 109.61561 * 353.33390 * 29 8 6 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 6 2.0400 1.4426 0.0000 4 6 2.0401 -0.7213 1.2492 5 7 3.5001 -0.8264 1.1905 6 6 4.1732 -1.4283 2.1911 7 8 3.5645 -1.9651 3.0965 8 7 5.5209 -1.4352 2.1880 9 6 6.2697 -0.8886 1.0482 10 6 7.2687 0.1504 1.5672 11 6 8.1454 -0.4862 2.6482 12 6 7.2682 -0.9391 3.8176 13 1 6.7268 -0.0826 4.2193 14 6 8.1490 -1.5456 4.9117 15 6 7.3007 -1.8803 6.1114 16 8 6.1134 -1.6571 6.0931 17 8 7.8642 -2.4263 7.2005 18 6 6.9866 -2.7219 8.3189 19 1 5.9882 -2.9511 7.9466 20 6 6.9183 -1.5082 9.2482 21 6 7.5235 -3.9071 9.0792 22 1 8.4309 -4.3927 8.7518 23 6 6.8535 -4.3747 10.1899 24 7 5.7457 -3.7867 10.5854 25 14 7.5194 -5.8434 11.1327 26 1 8.7692 -6.3202 10.4879 27 1 7.8060 -5.4436 12.5338 28 1 6.5119 -6.9344 11.1254 29 6 6.2686 -1.9903 3.3246 30 1 -0.3633 0.5138 0.8899 31 1 -0.3633 0.5138 -0.8899 32 1 -0.3633 -1.0276 -0.0005 33 1 1.8933 -0.5138 -0.8900 34 1 1.6767 1.9563 -0.8900 35 1 1.6767 1.9563 0.8900 36 1 3.1300 1.4426 -0.0005 37 1 1.6054 -1.7197 1.2951 38 1 1.7517 -0.1583 2.1369 39 1 3.9825 -0.4622 0.4319 40 1 5.5791 -0.4149 0.3506 41 1 6.8069 -1.6923 0.5447 42 1 6.7269 0.9966 1.9899 43 1 7.8965 0.4938 0.7451 44 1 8.8734 0.2439 3.0019 45 1 8.6679 -1.3472 2.2314 46 1 8.9175 -0.8280 5.1989 47 1 8.6215 -2.4532 4.5362 48 1 6.2598 -1.7303 10.0880 49 1 7.9167 -1.2791 9.6205 50 1 6.5300 -0.6510 8.6983 51 1 5.2741 -4.1159 11.3666 52 1 6.8054 -2.8841 3.0068 53 1 5.5777 -2.2429 4.1288 RHF calculation, no. of doubly occupied orbitals= 66 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 ZINC000122572867.mol2 53 1SCF WAS SPECIFIED, SO GEOMETRY OPTIMIZATION WAS NOT USED SCF FIELD WAS ACHIEVED Wall clock time and date at job start Sat Apr 11 2020 06:11:04 Heat of formation + Delta-G solvation = -137.608571 kcal Electronic energy + Delta-G solvation = -30733.679619 eV Core-core repulsion = 26604.750205 eV Total energy + Delta-G solvation = -4128.929414 eV No. of doubly occupied orbitals = 66 Molecular weight (most abundant/longest-lived isotopes) = 340.220 amu Computer time = 0.78 seconds Orbital eigenvalues (eV) -40.57897 -39.51373 -37.99208 -36.71842 -35.69562 -34.02961 -33.20144 -32.90204 -31.91208 -29.63576 -27.80442 -27.56885 -26.62272 -25.26186 -24.75921 -23.11822 -21.63141 -20.79636 -20.58278 -19.55198 -19.34902 -18.14581 -16.84446 -16.58206 -16.47081 -15.85011 -15.57303 -15.33256 -15.00721 -14.88307 -14.58956 -14.50324 -14.28484 -14.05640 -13.41494 -13.40322 -13.27224 -13.08448 -12.76076 -12.66259 -12.52300 -12.32602 -12.09883 -11.92704 -11.59836 -11.50027 -11.31284 -11.28358 -11.17350 -11.09981 -11.04296 -10.94308 -10.82409 -10.76343 -10.56156 -10.28032 -10.22618 -9.99529 -9.50437 -9.36340 -9.16842 -8.82328 -8.52241 -8.43557 -6.12603 -4.19319 2.59216 3.09838 3.27390 3.27474 3.33725 3.38688 3.69741 3.79037 4.37254 4.38230 4.50010 4.51130 4.65033 4.66304 4.75242 4.87457 4.99611 5.15162 5.19002 5.20476 5.25465 5.35301 5.37024 5.42117 5.44560 5.49514 5.51976 5.61567 5.63545 5.67010 5.72898 5.74618 5.86336 5.87985 5.91783 5.95737 5.98582 6.27381 6.33085 6.42166 6.53624 6.71973 6.91839 6.94993 6.99142 7.27040 7.66758 7.80492 7.92180 8.02499 8.21586 8.65442 8.72141 8.96429 9.45945 10.94886 Molecular weight = 340.22amu Principal moments of inertia in cm(-1) A = 0.013843 B = 0.002399 C = 0.002158 Principal moments of inertia in units of 10**(-40)*gram-cm**2 A = 2022.232124 B =11668.079842 C =12972.388258 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.108 4.108 3 C -0.148 4.148 4 C 0.140 3.860 5 N -0.750 5.750 6 C 0.704 3.296 7 O -0.575 6.575 8 N -0.627 5.627 9 C 0.112 3.888 10 C -0.134 4.134 11 C -0.114 4.114 12 C -0.090 4.090 13 H 0.084 0.916 14 C -0.110 4.110 15 C 0.445 3.555 16 O -0.512 6.512 17 O -0.359 6.359 18 C 0.214 3.786 19 H 0.047 0.953 20 C -0.154 4.154 21 C -0.519 4.519 22 H 0.069 0.931 23 C 0.061 3.939 24 N -0.888 5.888 25 Si 0.911 3.089 26 H -0.271 1.271 27 H -0.282 1.282 28 H -0.282 1.282 29 C 0.125 3.875 30 H 0.059 0.941 31 H 0.059 0.941 32 H 0.059 0.941 33 H 0.074 0.926 34 H 0.056 0.944 35 H 0.059 0.941 36 H 0.053 0.947 37 H 0.066 0.934 38 H 0.066 0.934 39 H 0.392 0.608 40 H 0.072 0.928 41 H 0.069 0.931 42 H 0.068 0.932 43 H 0.075 0.925 44 H 0.069 0.931 45 H 0.067 0.933 46 H 0.102 0.898 47 H 0.100 0.900 48 H 0.098 0.902 49 H 0.035 0.965 50 H 0.015 0.985 51 H 0.265 0.735 52 H 0.060 0.940 53 H 0.118 0.882 For the charges calculated by CM2: Dipole moment (debyes) X Y Z Total from point charges -0.634 10.764 -25.912 28.066 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.203 4.203 2 C -0.127 4.127 3 C -0.206 4.206 4 C 0.016 3.984 5 N -0.409 5.409 6 C 0.405 3.595 7 O -0.451 6.451 8 N -0.364 5.364 9 C -0.012 4.012 10 C -0.172 4.172 11 C -0.152 4.152 12 C -0.109 4.109 13 H 0.102 0.898 14 C -0.148 4.148 15 C 0.288 3.712 16 O -0.402 6.402 17 O -0.275 6.275 18 C 0.162 3.838 19 H 0.065 0.935 20 C -0.211 4.211 21 C -0.557 4.557 22 H 0.087 0.913 23 C -0.141 4.141 24 N -0.526 5.526 25 Si 0.738 3.262 26 H -0.195 1.195 27 H -0.208 1.208 28 H -0.208 1.208 29 C 0.003 3.997 30 H 0.078 0.922 31 H 0.078 0.922 32 H 0.078 0.922 33 H 0.093 0.907 34 H 0.075 0.925 35 H 0.078 0.922 36 H 0.072 0.928 37 H 0.084 0.916 38 H 0.084 0.916 39 H 0.228 0.772 40 H 0.090 0.910 41 H 0.087 0.913 42 H 0.086 0.914 43 H 0.093 0.907 44 H 0.088 0.912 45 H 0.086 0.914 46 H 0.120 0.880 47 H 0.119 0.881 48 H 0.117 0.883 49 H 0.054 0.946 50 H 0.034 0.966 51 H 0.075 0.925 52 H 0.078 0.922 53 H 0.136 0.864 Dipole moment (debyes) X Y Z Total from point charges -1.618 10.394 -25.021 27.143 hybrid contribution 0.484 -0.443 0.207 0.688 sum -1.133 9.952 -24.814 26.759 Atomic orbital electron populations 1.21758 0.93930 1.02327 1.02273 1.21360 0.96802 0.96624 0.97880 1.21867 1.00895 0.95630 1.02167 1.21274 0.78692 1.00668 0.97764 1.45303 1.07611 1.64019 1.24005 1.15863 0.81535 0.79567 0.82554 1.90998 1.67036 1.49577 1.37452 1.47260 1.03702 1.64879 1.20600 1.21611 0.94627 0.96683 0.88304 1.21907 0.97879 0.98339 0.99067 1.21572 0.97558 1.00846 0.95191 1.21345 0.95853 0.99620 0.94113 0.89784 1.21350 0.98369 1.03461 0.91663 1.23368 0.87656 0.75128 0.85098 1.90636 1.16975 1.50095 1.82484 1.86472 1.45168 1.71092 1.24808 1.20618 0.89994 0.93968 0.79229 0.93494 1.22147 1.01052 0.96021 1.01918 1.21066 1.10859 1.12624 1.11175 0.91321 1.24903 0.94790 0.97550 0.96879 1.69920 1.16453 1.40597 1.25666 0.86228 0.78952 0.81630 0.79362 1.19503 1.20763 1.20776 1.21812 0.94919 0.94117 0.88881 0.92163 0.92222 0.92236 0.90717 0.92469 0.92229 0.92812 0.91561 0.91581 0.77195 0.91026 0.91309 0.91360 0.90674 0.91248 0.91397 0.87969 0.88123 0.88316 0.94557 0.96569 0.92499 0.92165 0.86426 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.15 -0.87 9.19 37.16 0.34 -0.53 16 2 C -0.11 -0.71 2.96 -90.61 -0.27 -0.98 16 3 C -0.15 -0.89 9.07 37.16 0.34 -0.55 16 4 C 0.14 1.41 5.07 -4.04 -0.02 1.39 16 5 N -0.75 -8.03 5.20 -65.23 -0.34 -8.37 16 6 C 0.70 10.74 8.29 -86.91 -0.72 10.02 16 7 O -0.57 -11.96 15.78 5.29 0.08 -11.87 16 8 N -0.63 -8.09 2.97 -178.49 -0.53 -8.62 16 9 C 0.11 0.89 6.43 -3.71 -0.02 0.87 16 10 C -0.13 -0.97 6.08 -26.61 -0.16 -1.13 16 11 C -0.11 -0.98 5.17 -26.68 -0.14 -1.12 16 12 C -0.09 -1.15 1.98 -90.50 -0.18 -1.33 16 13 H 0.08 1.26 7.96 -51.93 -0.41 0.84 16 14 C -0.11 -1.54 5.29 -27.88 -0.15 -1.68 16 15 C 0.44 9.29 7.68 36.00 0.28 9.57 16 16 O -0.51 -12.36 13.21 -18.75 -0.25 -12.61 16 17 O -0.36 -8.70 9.86 -39.27 -0.39 -9.09 16 18 C 0.21 5.96 2.56 -27.88 -0.07 5.88 16 19 H 0.05 1.45 7.16 -51.93 -0.37 1.08 16 20 C -0.15 -4.04 9.33 37.16 0.35 -3.69 16 21 C -0.52 -14.98 8.90 -34.44 -0.31 -15.29 16 22 H 0.07 1.95 7.81 -52.49 -0.41 1.54 16 23 C 0.06 1.72 5.30 -17.82 -0.09 1.63 16 24 N -0.89 -26.07 11.32 55.43 0.63 -25.44 16 25 Si 0.91 21.61 29.54 -169.99 -5.02 16.59 16 26 H -0.27 -6.28 7.11 56.52 0.40 -5.88 16 27 H -0.28 -6.55 7.11 56.52 0.40 -6.15 16 28 H -0.28 -6.63 7.11 56.52 0.40 -6.23 16 29 C 0.13 1.87 5.37 -3.70 -0.02 1.85 16 30 H 0.06 0.38 8.14 -51.93 -0.42 -0.04 16 31 H 0.06 0.30 8.14 -51.93 -0.42 -0.13 16 32 H 0.06 0.36 8.14 -51.93 -0.42 -0.06 16 33 H 0.07 0.44 8.14 -51.93 -0.42 0.01 16 34 H 0.06 0.28 8.14 -51.93 -0.42 -0.14 16 35 H 0.06 0.40 8.14 -51.93 -0.42 -0.02 16 36 H 0.05 0.31 7.18 -51.93 -0.37 -0.07 16 37 H 0.07 0.74 8.14 -51.93 -0.42 0.32 16 38 H 0.07 0.77 8.14 -51.93 -0.42 0.34 16 39 H 0.39 2.97 5.11 -40.82 -0.21 2.76 16 40 H 0.07 0.42 5.57 -51.93 -0.29 0.13 16 41 H 0.07 0.48 8.14 -51.93 -0.42 0.06 16 42 H 0.07 0.56 8.14 -51.93 -0.42 0.14 16 43 H 0.07 0.38 8.14 -51.93 -0.42 -0.05 16 44 H 0.07 0.53 8.14 -51.93 -0.42 0.11 16 45 H 0.07 0.52 8.14 -51.93 -0.42 0.10 16 46 H 0.10 1.20 8.14 -51.93 -0.42 0.78 16 47 H 0.10 1.22 8.14 -51.93 -0.42 0.79 16 48 H 0.10 2.77 6.07 -51.93 -0.32 2.45 16 49 H 0.04 0.89 8.14 -51.93 -0.42 0.47 16 50 H 0.02 0.38 8.14 -51.93 -0.42 -0.04 16 51 H 0.27 7.30 8.31 -40.82 -0.34 6.96 16 52 H 0.06 0.81 8.14 -51.93 -0.42 0.39 16 53 H 0.12 2.37 5.00 -51.93 -0.26 2.11 16 LS Contribution 414.57 15.07 6.25 6.25 Total: -1.00 -35.87 414.57 -9.80 -45.67 By element: Atomic # 1 Polarization: 11.97 SS G_CDS: -9.38 Total: 2.59 kcal Atomic # 6 Polarization: 5.75 SS G_CDS: -0.85 Total: 4.90 kcal Atomic # 7 Polarization: -42.19 SS G_CDS: -0.24 Total: -42.43 kcal Atomic # 8 Polarization: -33.02 SS G_CDS: -0.55 Total: -33.57 kcal Atomic # 14 Polarization: 21.61 SS G_CDS: -5.02 Total: 16.59 kcal Total LS contribution 6.25 Total: 6.25 kcal Total: -35.87 -9.80 -45.67 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. ZINC000122572867.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 -91.935 kcal (2) G-P(sol) polarization free energy of solvation -35.871 kcal (3) G-ENP(sol) elect.-nuc.-pol. free energy of system -127.806 kcal (4) G-CDS(sol) cavity-dispersion-solvent structure free energy -9.802 kcal (5) G-P-CDS(sol) = G-P(sol) + G-CDS(sol) = (2) + (4) -45.673 kcal (6) G-S(sol) free energy of system = (1) + (5) -137.609 kcal 1SCF WAS SPECIFIED, SO GEOMETRY OPTIMIZATION WAS NOT USED Total computer time = 0.78 seconds