Wall clock time and date at job start Tue Mar 31 2020 19:25:01 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.52997 * 1 3 3 C 1.52993 * 109.47057 * 2 1 4 4 C 1.53001 * 109.47140 * 119.99698 * 2 1 3 5 5 C 1.53002 * 109.47365 * 175.00161 * 4 2 1 6 6 H 1.09000 * 115.55766 * 29.29474 * 5 4 2 7 7 C 1.53002 * 117.49773 * 243.61158 * 5 4 2 8 8 C 1.53000 * 60.02940 * 252.51785 * 7 5 4 9 9 H 1.09008 * 117.49673 * 107.47884 * 8 7 5 10 10 C 1.50690 * 117.50063 * 252.49714 * 8 7 5 11 11 O 1.21283 * 120.00042 * 68.61154 * 10 8 7 12 12 N 1.34777 * 120.00347 * 248.60878 * 10 8 7 13 13 C 1.46501 * 120.00401 * 179.97438 * 12 10 8 14 14 C 1.52995 * 109.47437 * 180.02562 * 13 12 10 15 15 H 1.09004 * 110.26969 * 60.61580 * 14 13 12 16 16 C 1.55002 * 110.30410 * 182.79789 * 14 13 12 17 17 C 1.55456 * 101.03112 * 155.01048 * 16 14 13 18 18 C 1.53773 * 102.01483 * 326.40497 * 17 16 14 19 19 C 1.53265 * 114.73066 * 151.73839 * 18 17 16 20 20 N 1.47713 * 85.91537 * 92.61985 * 19 18 17 21 21 C 1.36936 * 134.43981 * 204.20291 * 20 19 18 22 22 H 1.07996 * 119.99459 * 179.97438 * 21 20 19 23 23 C 1.38810 * 119.99930 * 359.93451 * 21 20 19 24 24 N 1.31617 * 119.99818 * 0.02562 * 23 21 20 25 25 Si 1.86800 * 119.99838 * 179.97438 * 23 21 20 26 26 H 1.48499 * 109.47029 * 0.02562 * 25 23 21 27 27 H 1.48500 * 109.46968 * 120.00454 * 25 23 21 28 28 H 1.48495 * 109.47209 * 240.00409 * 25 23 21 29 29 C 1.47709 * 91.12298 * 24.15628 * 20 19 18 30 30 O 1.42534 * 108.37624 * 20.95664 * 18 17 16 31 31 H 1.09002 * 109.47363 * 299.99831 * 1 2 3 32 32 H 1.08997 * 109.47688 * 60.00513 * 1 2 3 33 33 H 1.09006 * 109.47017 * 180.02562 * 1 2 3 34 34 H 1.08997 * 109.47091 * 239.99293 * 2 1 3 35 35 H 1.08999 * 109.47519 * 179.97438 * 3 2 1 36 36 H 1.08999 * 109.47047 * 299.99787 * 3 2 1 37 37 H 1.09005 * 109.47503 * 59.99451 * 3 2 1 38 38 H 1.09007 * 109.47338 * 55.00337 * 4 2 1 39 39 H 1.08996 * 109.47269 * 295.00408 * 4 2 1 40 40 H 1.09002 * 117.49918 * 145.01819 * 7 5 4 41 41 H 1.09002 * 117.49754 * 0.02562 * 7 5 4 42 42 H 0.97000 * 120.00030 * 0.02562 * 12 10 8 43 43 H 1.09006 * 109.47359 * 300.00196 * 13 12 10 44 44 H 1.08997 * 109.47077 * 60.00008 * 13 12 10 45 45 H 1.09003 * 111.10849 * 272.94839 * 16 14 13 46 46 H 1.08996 * 111.11455 * 37.06523 * 16 14 13 47 47 H 1.08999 * 110.92190 * 208.23503 * 17 16 14 48 48 H 1.09000 * 110.92759 * 84.58205 * 17 16 14 49 49 H 1.08996 * 113.80281 * 206.87151 * 19 18 17 50 50 H 1.09003 * 113.80045 * 338.37090 * 19 18 17 51 51 H 0.97009 * 119.99790 * 180.02562 * 24 23 21 52 52 H 1.08997 * 113.80508 * 221.59645 * 29 20 19 53 53 H 1.08999 * 113.90921 * 90.10305 * 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.5300 0.0000 0.0000 3 6 2.0399 1.4424 0.0000 4 6 2.0400 -0.7212 1.2493 5 6 3.5648 -0.8310 1.1881 6 1 4.0720 -0.0356 0.6420 7 6 4.1697 -2.2345 1.1158 8 6 4.2804 -1.4200 2.4063 9 1 3.6646 -1.7291 3.2509 10 6 5.6325 -0.8562 2.7594 11 8 6.1139 0.0312 2.0872 12 7 6.3069 -1.3378 3.8222 13 6 7.6217 -0.7901 4.1652 14 6 8.1657 -1.5065 5.4027 15 1 8.2749 -2.5728 5.2043 16 6 9.5255 -0.8937 5.8246 17 6 9.5545 -1.2478 7.3381 18 6 8.0664 -1.1610 7.7157 19 6 7.6471 -2.0631 8.8817 20 7 7.9164 -0.8946 9.7442 21 6 8.2381 -0.7318 11.0652 22 1 8.3886 0.2596 11.4661 23 6 8.3699 -1.8411 11.8892 24 7 8.1862 -3.0493 11.4006 25 14 8.8095 -1.6191 13.6911 26 1 8.9621 -0.1723 13.9888 27 1 7.7297 -2.1871 14.5376 28 1 10.0857 -2.3215 13.9795 29 6 7.7008 0.0281 8.6111 30 8 7.2925 -1.2947 6.5262 31 1 -0.3634 0.5138 0.8900 32 1 -0.3634 0.5137 -0.8900 33 1 -0.3633 -1.0277 0.0005 34 1 1.8933 -0.5139 -0.8899 35 1 3.1299 1.4425 0.0005 36 1 1.6766 1.9562 -0.8900 37 1 1.6765 1.9564 0.8899 38 1 1.6054 -1.7198 1.2953 39 1 1.7516 -0.1583 2.1370 40 1 5.0749 -2.3625 0.5222 41 1 3.4813 -3.0796 1.1115 42 1 5.9221 -2.0477 4.3597 43 1 8.3062 -0.9381 3.3299 44 1 7.5281 0.2757 4.3737 45 1 10.3513 -1.3683 5.2946 46 1 9.5335 0.1857 5.6738 47 1 10.1409 -0.5187 7.8972 48 1 9.9377 -2.2561 7.4948 49 1 6.5981 -2.3576 8.8548 50 1 8.3219 -2.9017 9.0538 51 1 8.2779 -3.8246 11.9765 52 1 8.4180 0.8477 8.5685 53 1 6.6683 0.3599 8.5018 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 ZINC001053665713.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:25:01 Heat of formation + Delta-G solvation = -23.352211 kcal Electronic energy + Delta-G solvation = -28961.516981 eV Core-core repulsion = 25027.998146 eV Total energy + Delta-G solvation = -3933.518835 eV No. of doubly occupied orbitals = 65 Molecular weight (most abundant/longest-lived isotopes) = 336.220 amu Computer time = 0.93 seconds Orbital eigenvalues (eV) -40.26183 -39.72121 -37.87637 -36.55568 -35.25391 -34.44140 -32.09019 -31.86762 -29.44365 -27.92483 -27.45179 -26.95799 -24.93951 -24.67436 -23.07004 -21.86136 -21.54658 -21.20618 -20.11568 -19.95464 -18.75668 -17.46308 -16.90342 -16.53073 -16.04523 -15.63308 -15.15883 -15.12538 -14.83630 -14.62104 -14.35565 -14.07739 -13.82107 -13.66077 -13.33042 -13.08637 -12.94616 -12.75951 -12.52706 -12.44204 -12.33228 -12.13283 -11.94242 -11.77410 -11.41541 -11.36404 -11.21632 -11.12220 -10.98808 -10.89418 -10.84130 -10.69069 -10.51505 -10.46987 -10.34674 -10.25281 -10.01400 -9.82733 -9.59008 -9.09111 -8.83207 -8.73213 -6.91937 -5.80393 -3.12649 2.42391 3.04526 3.36756 3.42171 3.44818 3.60057 3.98759 4.30429 4.45519 4.46695 4.53803 4.68570 4.73771 4.75180 4.83301 4.87533 5.00480 5.08638 5.10557 5.17065 5.28369 5.31815 5.32227 5.39736 5.41535 5.46250 5.47765 5.53595 5.56101 5.60465 5.67145 5.70026 5.71532 5.74487 5.85689 5.91324 5.93972 6.07849 6.08872 6.15951 6.34392 6.34809 6.41502 6.49151 6.58767 6.71965 7.01922 7.21187 7.48511 7.49995 7.69641 7.94076 8.18925 9.01376 9.67060 10.55885 11.37121 Molecular weight = 336.22amu Principal moments of inertia in cm(-1) A = 0.018367 B = 0.001901 C = 0.001807 Principal moments of inertia in units of 10**(-40)*gram-cm**2 A = 1524.083082 B =14726.017154 C =15489.813339 Net atomic charges, atomic populations, and dipole contributions calculated with CM2 Atom NO. Type Charge No. of electrons 1 C -0.148 4.148 2 C -0.094 4.094 3 C -0.151 4.151 4 C -0.087 4.087 5 C -0.118 4.118 6 H 0.120 0.880 7 C -0.150 4.150 8 C -0.191 4.191 9 H 0.124 0.876 10 C 0.542 3.458 11 O -0.549 6.549 12 N -0.716 5.716 13 C 0.124 3.876 14 C 0.067 3.933 15 H 0.074 0.926 16 C -0.148 4.148 17 C -0.112 4.112 18 C 0.056 3.944 19 C 0.165 3.835 20 N -0.605 5.605 21 C -0.219 4.219 22 H 0.082 0.918 23 C -0.008 4.008 24 N -0.921 5.921 25 Si 0.913 3.087 26 H -0.279 1.279 27 H -0.290 1.290 28 H -0.290 1.290 29 C 0.118 3.882 30 O -0.365 6.365 31 H 0.054 0.946 32 H 0.054 0.946 33 H 0.053 0.947 34 H 0.071 0.929 35 H 0.060 0.940 36 H 0.052 0.948 37 H 0.054 0.946 38 H 0.066 0.934 39 H 0.069 0.931 40 H 0.097 0.903 41 H 0.103 0.897 42 H 0.407 0.593 43 H 0.076 0.924 44 H 0.073 0.927 45 H 0.080 0.920 46 H 0.078 0.922 47 H 0.084 0.916 48 H 0.078 0.922 49 H 0.032 0.968 50 H 0.076 0.924 51 H 0.253 0.747 52 H 0.045 0.955 53 H 0.040 0.960 For the charges calculated by CM2: Dipole moment (debyes) X Y Z Total from point charges -8.210 1.601 -23.556 24.997 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.206 4.206 2 C -0.112 4.112 3 C -0.209 4.209 4 C -0.124 4.124 5 C -0.136 4.136 6 H 0.138 0.862 7 C -0.187 4.187 8 C -0.212 4.212 9 H 0.142 0.858 10 C 0.328 3.672 11 O -0.427 6.427 12 N -0.367 5.367 13 C 0.001 3.999 14 C 0.008 3.992 15 H 0.092 0.908 16 C -0.186 4.186 17 C -0.150 4.150 18 C 0.016 3.984 19 C 0.039 3.961 20 N -0.333 5.333 21 C -0.347 4.347 22 H 0.099 0.901 23 C -0.205 4.205 24 N -0.568 5.568 25 Si 0.743 3.257 26 H -0.204 1.204 27 H -0.217 1.217 28 H -0.217 1.217 29 C -0.008 4.008 30 O -0.283 6.283 31 H 0.074 0.926 32 H 0.073 0.927 33 H 0.072 0.928 34 H 0.089 0.911 35 H 0.079 0.921 36 H 0.071 0.929 37 H 0.073 0.927 38 H 0.084 0.916 39 H 0.088 0.912 40 H 0.116 0.884 41 H 0.122 0.878 42 H 0.242 0.758 43 H 0.094 0.906 44 H 0.091 0.909 45 H 0.098 0.902 46 H 0.097 0.903 47 H 0.103 0.897 48 H 0.097 0.903 49 H 0.050 0.950 50 H 0.094 0.906 51 H 0.064 0.936 52 H 0.063 0.937 53 H 0.057 0.943 Dipole moment (debyes) X Y Z Total from point charges -7.996 2.154 -23.704 25.109 hybrid contribution 0.156 -0.352 1.543 1.591 sum -7.840 1.802 -22.161 23.575 Atomic orbital electron populations 1.21849 0.95251 1.01643 1.01813 1.21064 0.95458 0.96504 0.98190 1.21916 1.01603 0.95740 1.01601 1.21298 0.91964 1.00767 0.98363 1.22517 0.95105 0.97318 0.98679 0.86208 1.23098 1.03369 0.95033 0.97202 1.22316 0.96609 1.07587 0.94650 0.85795 1.19644 0.85545 0.81943 0.80020 1.90677 1.69564 1.30637 1.51804 1.45835 1.20081 1.40552 1.30243 1.21570 0.83885 0.98835 0.95647 1.23168 0.90164 0.98741 0.87095 0.90793 1.23267 0.96501 1.01706 0.97166 1.22614 0.95631 1.03134 0.93616 1.24304 0.90268 0.99796 0.84018 1.22506 0.98386 0.92177 0.83001 1.45659 1.80419 1.01903 1.05301 1.19024 1.39851 0.93139 0.82669 0.90060 1.24794 1.02160 0.95200 0.98317 1.69832 1.55306 0.96993 1.34638 0.86029 0.77669 0.77375 0.84595 1.20425 1.21720 1.21673 1.22808 0.98170 0.93242 0.86629 1.88349 1.32004 1.92750 1.15219 0.92644 0.92731 0.92756 0.91094 0.92111 0.92934 0.92658 0.91563 0.91226 0.88417 0.87848 0.75783 0.90631 0.90857 0.90155 0.90346 0.89727 0.90344 0.94993 0.90616 0.93646 0.93678 0.94270 Number of geometries 1 Number of calculations of the screened coulomb radii 1 The total number of SCF iterations 14. 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.09 -0.63 2.92 -90.62 -0.26 -0.89 16 3 C -0.15 -1.10 9.01 37.15 0.33 -0.77 16 4 C -0.09 -0.62 4.96 -26.73 -0.13 -0.75 16 5 C -0.12 -1.03 5.04 -90.58 -0.46 -1.49 16 6 H 0.12 1.25 5.91 -51.93 -0.31 0.94 16 7 C -0.15 -1.17 9.88 -26.73 -0.26 -1.44 16 8 C -0.19 -1.81 5.91 -91.74 -0.54 -2.36 16 9 H 0.12 1.02 8.14 -51.92 -0.42 0.60 16 10 C 0.54 6.87 7.61 -10.99 -0.08 6.78 16 11 O -0.55 -8.79 16.26 5.56 0.09 -8.70 16 12 N -0.72 -8.47 5.48 -60.30 -0.33 -8.80 16 13 C 0.12 1.52 5.40 -4.04 -0.02 1.50 16 14 C 0.07 0.87 3.52 -25.68 -0.09 0.78 16 15 H 0.07 0.89 8.14 -51.93 -0.42 0.47 16 16 C -0.15 -1.84 6.58 -24.37 -0.16 -2.01 16 17 C -0.11 -1.83 6.32 -25.02 -0.16 -1.99 16 18 C 0.06 1.08 1.68 -90.05 -0.15 0.92 16 19 C 0.16 3.93 7.01 -3.25 -0.02 3.91 16 20 N -0.60 -15.57 3.67 -181.29 -0.67 -16.23 16 21 C -0.22 -6.14 10.70 -15.06 -0.16 -6.30 16 22 H 0.08 2.20 7.83 -52.49 -0.41 1.79 16 23 C -0.01 -0.23 5.50 -17.43 -0.10 -0.33 16 24 N -0.92 -27.08 11.39 55.49 0.63 -26.45 16 25 Si 0.91 22.01 29.55 -169.99 -5.02 16.98 16 26 H -0.28 -6.51 7.11 56.52 0.40 -6.11 16 27 H -0.29 -7.06 7.11 56.52 0.40 -6.66 16 28 H -0.29 -6.98 7.11 56.52 0.40 -6.58 16 29 C 0.12 2.50 8.37 -3.24 -0.03 2.48 16 30 O -0.36 -6.18 10.67 -35.23 -0.38 -6.56 16 31 H 0.05 0.33 8.14 -51.93 -0.42 -0.09 16 32 H 0.05 0.30 8.14 -51.93 -0.42 -0.12 16 33 H 0.05 0.29 8.14 -51.93 -0.42 -0.13 16 34 H 0.07 0.47 8.14 -51.93 -0.42 0.05 16 35 H 0.06 0.50 6.41 -51.93 -0.33 0.17 16 36 H 0.05 0.34 8.14 -51.93 -0.42 -0.08 16 37 H 0.05 0.41 8.14 -51.93 -0.42 -0.01 16 38 H 0.07 0.42 7.88 -51.93 -0.41 0.01 16 39 H 0.07 0.53 8.14 -51.93 -0.42 0.11 16 40 H 0.10 0.81 8.14 -51.93 -0.42 0.39 16 41 H 0.10 0.63 8.03 -51.93 -0.42 0.21 16 42 H 0.41 4.32 8.41 -40.82 -0.34 3.97 16 43 H 0.08 0.85 8.14 -51.93 -0.42 0.43 16 44 H 0.07 1.00 8.14 -51.93 -0.42 0.58 16 45 H 0.08 0.82 8.14 -51.93 -0.42 0.40 16 46 H 0.08 0.94 8.05 -51.93 -0.42 0.52 16 47 H 0.08 1.44 7.83 -51.93 -0.41 1.03 16 48 H 0.08 1.29 8.01 -51.93 -0.42 0.88 16 49 H 0.03 0.80 8.14 -51.93 -0.42 0.38 16 50 H 0.08 1.95 7.05 -51.93 -0.37 1.58 16 51 H 0.25 7.20 8.31 -40.82 -0.34 6.86 16 52 H 0.05 0.90 7.98 -51.93 -0.41 0.48 16 53 H 0.04 0.84 8.14 -51.93 -0.42 0.42 16 LS Contribution 421.79 15.07 6.36 6.36 Total: -1.00 -32.42 421.79 -10.99 -43.41 By element: Atomic # 1 Polarization: 12.19 SS G_CDS: -9.72 Total: 2.47 kcal Atomic # 6 Polarization: -0.52 SS G_CDS: -1.96 Total: -2.48 kcal Atomic # 7 Polarization: -51.13 SS G_CDS: -0.36 Total: -51.49 kcal Atomic # 8 Polarization: -14.97 SS G_CDS: -0.29 Total: -15.26 kcal Atomic # 14 Polarization: 22.01 SS G_CDS: -5.02 Total: 16.98 kcal Total LS contribution 6.36 Total: 6.36 kcal Total: -32.42 -10.99 -43.41 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. ZINC001053665713.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 20.062 kcal (2) G-P(sol) polarization free energy of solvation -32.424 kcal (3) G-ENP(sol) elect.-nuc.-pol. free energy of system -12.361 kcal (4) G-CDS(sol) cavity-dispersion-solvent structure free energy -10.991 kcal (5) G-P-CDS(sol) = G-P(sol) + G-CDS(sol) = (2) + (4) -43.414 kcal (6) G-S(sol) free energy of system = (1) + (5) -23.352 kcal 1SCF WAS SPECIFIED, SO GEOMETRY OPTIMIZATION WAS NOT USED Total computer time = 0.93 seconds