Wall clock time and date at job start Sat Apr 11 2020 06:09:49 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.52999 * 1 3 3 H 1.08997 * 109.47281 * 2 1 4 4 O 1.42902 * 109.46995 * 239.99784 * 2 1 3 5 5 C 1.53006 * 109.47134 * 119.99976 * 2 1 3 6 6 C 1.52997 * 109.47239 * 186.31443 * 5 2 1 7 7 H 1.09002 * 109.70079 * 58.69695 * 6 5 2 8 8 C 1.52759 * 109.50931 * 298.57636 * 6 5 2 9 9 C 1.53041 * 109.38662 * 185.70401 * 8 6 5 10 10 C 1.53043 * 109.53513 * 61.24149 * 9 8 6 11 11 C 1.53191 * 109.31160 * 298.63867 * 10 9 8 12 12 N 1.46925 * 109.58522 * 179.05007 * 6 5 2 13 13 C 1.34767 * 120.63360 * 353.42137 * 12 6 5 14 14 O 1.21356 * 120.00101 * 7.16548 * 13 12 6 15 15 C 1.49453 * 120.00253 * 187.17011 * 13 12 6 16 16 O 1.21342 * 119.99978 * 6.30975 * 15 13 12 17 17 N 1.34774 * 119.99744 * 186.31564 * 15 13 12 18 18 C 1.46507 * 119.99668 * 179.97438 * 17 15 13 19 19 C 1.53001 * 109.46995 * 180.02562 * 18 17 15 20 20 C 1.52997 * 109.47086 * 60.00456 * 19 18 17 21 21 C 1.53001 * 109.47257 * 299.99888 * 19 18 17 22 22 C 1.50704 * 109.46929 * 180.02562 * 19 18 17 23 23 H 1.08003 * 120.00046 * 3.22539 * 22 19 18 24 24 C 1.37875 * 119.99857 * 183.23497 * 22 19 18 25 25 N 1.31510 * 120.00123 * 4.32060 * 24 22 19 26 26 Si 1.86797 * 119.99992 * 184.32164 * 24 22 19 27 27 H 1.48498 * 109.47542 * 265.00063 * 26 24 22 28 28 H 1.48503 * 109.47314 * 144.99754 * 26 24 22 29 29 H 1.48507 * 109.46973 * 25.00064 * 26 24 22 30 30 H 1.09002 * 109.47213 * 299.99813 * 1 2 3 31 31 H 1.08996 * 109.47053 * 60.00161 * 1 2 3 32 32 H 1.08997 * 109.47281 * 179.97438 * 1 2 3 33 33 H 0.96704 * 113.99997 * 299.99480 * 4 2 1 34 34 H 1.08997 * 109.47153 * 66.31098 * 5 2 1 35 35 H 1.08997 * 109.46626 * 306.31478 * 5 2 1 36 36 H 1.09003 * 109.56500 * 65.67162 * 8 6 5 37 37 H 1.08998 * 109.39184 * 305.63668 * 8 6 5 38 38 H 1.09001 * 109.45583 * 301.22564 * 9 8 6 39 39 H 1.08996 * 109.45217 * 181.25812 * 9 8 6 40 40 H 1.09000 * 109.49692 * 178.68701 * 10 9 8 41 41 H 1.09001 * 109.49645 * 58.58383 * 10 9 8 42 42 H 1.09008 * 109.58624 * 294.84394 * 11 10 9 43 43 H 1.09002 * 109.59040 * 174.42429 * 11 10 9 44 44 H 0.96996 * 120.00582 * 359.97438 * 17 15 13 45 45 H 1.09001 * 109.46537 * 60.00226 * 18 17 15 46 46 H 1.08995 * 109.47281 * 299.99921 * 18 17 15 47 47 H 1.09002 * 109.46803 * 175.35796 * 20 19 18 48 48 H 1.08998 * 109.46699 * 295.35132 * 20 19 18 49 49 H 1.08998 * 109.47311 * 55.35621 * 20 19 18 50 50 H 1.08996 * 109.47164 * 300.00642 * 21 19 18 51 51 H 1.09001 * 109.46770 * 60.00119 * 21 19 18 52 52 H 1.08995 * 109.47421 * 180.02562 * 21 19 18 53 53 H 0.97003 * 119.99974 * 179.97438 * 25 24 22 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 8 2.0063 -0.6737 -1.1668 5 6 2.0400 -0.7213 1.2493 6 6 3.5617 -0.5854 1.3312 7 1 3.8353 0.4693 1.3628 8 6 4.1989 -1.2393 0.1065 9 6 5.7217 -1.2283 0.2592 10 6 6.1174 -2.0330 1.4994 11 6 5.5019 -1.3831 2.7427 12 7 4.0518 -1.2578 2.5421 13 6 3.1877 -1.7552 3.4488 14 8 1.9979 -1.7671 3.2102 15 6 3.6946 -2.2920 4.7482 16 8 4.8678 -2.1766 5.0360 17 7 2.8488 -2.9043 5.6002 18 6 3.3458 -3.4311 6.8738 19 6 2.1911 -4.0724 7.6459 20 6 1.5966 -5.2142 6.8191 21 6 1.1118 -3.0218 7.9150 22 6 2.7023 -4.6136 8.9562 23 1 3.7348 -4.4648 9.2360 24 6 1.8519 -5.2990 9.7977 25 7 0.6182 -5.5544 9.4204 26 14 2.4521 -5.8646 11.4738 27 1 2.9222 -7.2703 11.3828 28 1 1.3395 -5.7773 12.4534 29 1 3.5718 -4.9960 11.9177 30 1 -0.3634 0.5138 0.8900 31 1 -0.3633 0.5138 -0.8900 32 1 -0.3634 -1.0276 -0.0005 33 1 1.7210 -1.5953 -1.2328 34 1 1.7719 -1.7763 1.1946 35 1 1.5873 -0.2768 2.1356 36 1 3.9193 -0.6857 -0.7899 37 1 3.8492 -2.2682 0.0225 38 1 6.0683 -0.2005 0.3674 39 1 6.1783 -1.6745 -0.6241 40 1 7.2032 -2.0443 1.5946 41 1 5.7494 -3.0544 1.4030 42 1 5.9370 -0.3948 2.8917 43 1 5.6976 -2.0061 3.6155 44 1 1.9109 -2.9963 5.3703 45 1 3.7698 -2.6176 7.4626 46 1 4.1144 -4.1798 6.6821 47 1 0.8294 -5.7237 7.4021 48 1 1.1531 -4.8108 5.9087 49 1 2.3837 -5.9214 6.5574 50 1 1.5354 -2.2084 8.5040 51 1 0.7420 -2.6299 6.9675 52 1 0.2893 -3.4785 8.4654 53 1 0.0201 -6.0369 10.0124 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 ZINC000422488447.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:09:49 Heat of formation + Delta-G solvation = -109.648258 kcal Electronic energy + Delta-G solvation = -31686.359212 eV Core-core repulsion = 27558.642248 eV Total energy + Delta-G solvation = -4127.716964 eV No. of doubly occupied orbitals = 66 Molecular weight (most abundant/longest-lived isotopes) = 340.220 amu Computer time = 0.71 seconds Orbital eigenvalues (eV) -41.11747 -39.19242 -37.84511 -37.31270 -36.48019 -33.59908 -32.96525 -32.66204 -30.93602 -30.86243 -27.51742 -26.26149 -25.74905 -25.36761 -24.72881 -23.62996 -22.27266 -21.52206 -20.69493 -19.62813 -18.80269 -17.52382 -16.92127 -16.58028 -15.99488 -15.84928 -15.64969 -15.37132 -15.02420 -14.84540 -14.57791 -14.45805 -14.31872 -13.77394 -13.64994 -13.28703 -13.16290 -12.96274 -12.79842 -12.70525 -12.41556 -12.26002 -12.00993 -11.81266 -11.54382 -11.45978 -11.21531 -11.19744 -11.12050 -11.06065 -10.90835 -10.83300 -10.72110 -10.44951 -10.19209 -10.03737 -9.95372 -9.68721 -9.44267 -8.87485 -8.79996 -8.60408 -8.55221 -8.39620 -6.04052 -4.16731 1.62878 2.50518 3.25565 3.33451 3.37085 3.41469 3.71688 4.13481 4.35913 4.45447 4.65000 4.70179 4.79012 4.82556 4.92945 4.95726 5.04332 5.16794 5.20540 5.30012 5.39879 5.47210 5.57670 5.66385 5.75709 5.79219 5.85671 5.90554 5.98838 6.01719 6.07610 6.11474 6.15355 6.19983 6.25426 6.43761 6.49406 6.66939 6.81569 6.84472 6.92466 6.97787 7.14012 7.30109 7.40176 7.50838 7.60363 7.66712 7.92884 7.98859 8.14205 8.32340 8.39946 8.83778 9.49743 10.95158 Molecular weight = 340.22amu Principal moments of inertia in cm(-1) A = 0.018021 B = 0.002601 C = 0.002354 Principal moments of inertia in units of 10**(-40)*gram-cm**2 A = 1553.393268 B =10763.358808 C =11892.899155 Net atomic charges, atomic populations, and dipole contributions calculated with CM2 Atom NO. Type Charge No. of electrons 1 C -0.181 4.181 2 C 0.093 3.907 3 H 0.094 0.906 4 O -0.563 6.563 5 C -0.141 4.141 6 C 0.147 3.853 7 H 0.071 0.929 8 C -0.136 4.136 9 C -0.123 4.123 10 C -0.136 4.136 11 C 0.065 3.935 12 N -0.580 5.580 13 C 0.508 3.492 14 O -0.520 6.520 15 C 0.494 3.506 16 O -0.528 6.528 17 N -0.677 5.677 18 C 0.127 3.873 19 C 0.056 3.944 20 C -0.140 4.140 21 C -0.130 4.130 22 C -0.507 4.507 23 H 0.058 0.942 24 C 0.067 3.933 25 N -0.884 5.884 26 Si 0.894 3.106 27 H -0.282 1.282 28 H -0.286 1.286 29 H -0.274 1.274 30 H 0.066 0.934 31 H 0.059 0.941 32 H 0.059 0.941 33 H 0.373 0.627 34 H 0.082 0.918 35 H 0.095 0.905 36 H 0.105 0.895 37 H 0.061 0.939 38 H 0.060 0.940 39 H 0.062 0.938 40 H 0.073 0.927 41 H 0.064 0.936 42 H 0.064 0.936 43 H 0.148 0.852 44 H 0.411 0.589 45 H 0.061 0.939 46 H 0.062 0.938 47 H 0.101 0.899 48 H 0.005 0.995 49 H 0.024 0.976 50 H 0.029 0.971 51 H 0.011 0.989 52 H 0.077 0.923 53 H 0.260 0.740 For the charges calculated by CM2: Dipole moment (debyes) X Y Z Total from point charges 4.218 9.474 -19.327 21.934 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.238 4.238 2 C 0.034 3.966 3 H 0.111 0.889 4 O -0.369 6.369 5 C -0.179 4.179 6 C 0.045 3.955 7 H 0.089 0.911 8 C -0.174 4.174 9 C -0.161 4.161 10 C -0.174 4.174 11 C -0.056 4.056 12 N -0.314 5.314 13 C 0.291 3.709 14 O -0.397 6.397 15 C 0.275 3.725 16 O -0.408 6.408 17 N -0.328 5.328 18 C 0.005 3.995 19 C 0.055 3.945 20 C -0.198 4.198 21 C -0.187 4.187 22 C -0.546 4.546 23 H 0.077 0.923 24 C -0.135 4.135 25 N -0.522 5.522 26 Si 0.722 3.278 27 H -0.208 1.208 28 H -0.212 1.212 29 H -0.198 1.198 30 H 0.085 0.915 31 H 0.078 0.922 32 H 0.078 0.922 33 H 0.221 0.779 34 H 0.100 0.900 35 H 0.114 0.886 36 H 0.123 0.877 37 H 0.079 0.921 38 H 0.079 0.921 39 H 0.081 0.919 40 H 0.092 0.908 41 H 0.082 0.918 42 H 0.083 0.917 43 H 0.164 0.836 44 H 0.248 0.752 45 H 0.080 0.920 46 H 0.081 0.919 47 H 0.120 0.880 48 H 0.024 0.976 49 H 0.043 0.957 50 H 0.048 0.952 51 H 0.030 0.970 52 H 0.095 0.905 53 H 0.071 0.929 Dipole moment (debyes) X Y Z Total from point charges 4.495 10.028 -19.238 22.156 hybrid contribution -0.314 -0.824 0.938 1.288 sum 4.181 9.204 -18.300 20.906 Atomic orbital electron populations 1.22419 0.96259 1.02699 1.02449 1.22113 0.93093 0.96568 0.84791 0.88852 1.86533 1.76326 1.29299 1.44780 1.22910 0.93594 1.02969 0.98388 1.21310 0.94588 0.95776 0.83829 0.91080 1.22218 0.94648 1.01885 0.98632 1.21593 0.96360 1.00280 0.97878 1.21955 1.00703 0.99809 0.94912 1.22982 0.77784 1.03291 1.01589 1.48043 1.07525 1.55454 1.20341 1.20363 0.86279 0.79061 0.85216 1.90772 1.15965 1.55082 1.77833 1.20166 0.86563 0.80247 0.85564 1.90405 1.16239 1.58587 1.75539 1.46247 1.12999 1.54133 1.19387 1.21969 0.98991 0.97472 0.81036 1.18380 0.90234 0.92186 0.93653 1.21981 1.01818 0.99225 0.96784 1.21804 1.01205 0.97846 0.97888 1.20962 0.96617 1.35169 1.01825 0.92334 1.24880 0.94921 0.95394 0.98275 1.69763 1.05747 1.40709 1.35978 0.86547 0.78462 0.78938 0.83872 1.20776 1.21183 1.19829 0.91502 0.92174 0.92217 0.77903 0.89977 0.88634 0.87742 0.92065 0.92108 0.91947 0.90813 0.91770 0.91726 0.83576 0.75190 0.92019 0.91922 0.88024 0.97571 0.95706 0.95200 0.97027 0.90457 0.92924 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.18 -1.45 9.53 37.16 0.35 -1.09 16 2 C 0.09 0.79 3.50 -26.73 -0.09 0.70 16 3 H 0.09 0.77 8.14 -51.93 -0.42 0.34 16 4 O -0.56 -4.61 9.94 -35.23 -0.35 -4.96 16 5 C -0.14 -1.48 2.93 -26.73 -0.08 -1.56 16 6 C 0.15 1.55 2.99 -67.82 -0.20 1.34 16 7 H 0.07 0.71 8.14 -51.93 -0.42 0.29 16 8 C -0.14 -1.23 4.32 -26.83 -0.12 -1.35 16 9 C -0.12 -1.06 5.89 -26.69 -0.16 -1.21 16 10 C -0.14 -1.31 6.08 -26.61 -0.16 -1.47 16 11 C 0.06 0.76 6.02 -3.71 -0.02 0.74 16 12 N -0.58 -7.33 2.98 -165.97 -0.49 -7.82 16 13 C 0.51 7.68 7.66 -11.59 -0.09 7.59 16 14 O -0.52 -8.24 11.91 5.49 0.07 -8.17 16 15 C 0.49 8.43 7.76 -11.59 -0.09 8.34 16 16 O -0.53 -9.55 13.14 5.50 0.07 -9.48 16 17 N -0.68 -12.20 4.98 -60.73 -0.30 -12.51 16 18 C 0.13 2.59 4.50 -4.04 -0.02 2.57 16 19 C 0.06 1.28 0.50 -155.66 -0.08 1.20 16 20 C -0.14 -3.19 8.04 37.16 0.30 -2.89 16 21 C -0.13 -2.90 8.11 37.16 0.30 -2.59 16 22 C -0.51 -12.98 8.01 -34.44 -0.28 -13.26 16 23 H 0.06 1.50 7.75 -52.48 -0.41 1.09 16 24 C 0.07 1.73 5.14 -17.82 -0.09 1.63 16 25 N -0.88 -23.29 9.36 55.43 0.52 -22.78 16 26 Si 0.89 20.37 29.54 -169.99 -5.02 15.35 16 27 H -0.28 -6.42 7.11 56.52 0.40 -6.02 16 28 H -0.29 -6.52 7.11 56.52 0.40 -6.12 16 29 H -0.27 -6.22 7.11 56.52 0.40 -5.82 16 30 H 0.07 0.55 8.14 -51.93 -0.42 0.12 16 31 H 0.06 0.45 8.14 -51.93 -0.42 0.03 16 32 H 0.06 0.48 8.14 -51.93 -0.42 0.06 16 33 H 0.37 2.21 9.05 45.56 0.41 2.62 16 34 H 0.08 0.98 7.05 -51.93 -0.37 0.62 16 35 H 0.10 1.19 6.43 -51.93 -0.33 0.86 16 36 H 0.10 0.92 5.71 -51.93 -0.30 0.62 16 37 H 0.06 0.60 8.14 -51.93 -0.42 0.18 16 38 H 0.06 0.50 8.14 -51.93 -0.42 0.08 16 39 H 0.06 0.48 8.14 -51.93 -0.42 0.06 16 40 H 0.07 0.62 8.14 -51.93 -0.42 0.20 16 41 H 0.06 0.67 8.14 -51.93 -0.42 0.25 16 42 H 0.06 0.71 8.14 -51.92 -0.42 0.28 16 43 H 0.15 2.06 4.22 -55.26 -0.23 1.83 16 44 H 0.41 7.30 5.68 -40.82 -0.23 7.07 16 45 H 0.06 1.29 8.14 -51.93 -0.42 0.87 16 46 H 0.06 1.31 8.14 -51.93 -0.42 0.89 16 47 H 0.10 2.57 5.53 -51.93 -0.29 2.28 16 48 H 0.01 0.11 7.00 -51.93 -0.36 -0.25 16 49 H 0.02 0.55 8.14 -51.93 -0.42 0.12 16 50 H 0.03 0.65 8.14 -51.93 -0.42 0.23 16 51 H 0.01 0.22 6.86 -51.93 -0.36 -0.14 16 52 H 0.08 1.90 6.50 -51.93 -0.34 1.56 16 53 H 0.26 6.63 8.31 -40.82 -0.34 6.29 16 LS Contribution 396.38 15.07 5.97 5.97 Total: -1.00 -26.90 396.38 -8.34 -35.23 By element: Atomic # 1 Polarization: 18.75 SS G_CDS: -8.28 Total: 10.48 kcal Atomic # 6 Polarization: -0.79 SS G_CDS: -0.52 Total: -1.31 kcal Atomic # 7 Polarization: -42.83 SS G_CDS: -0.28 Total: -43.11 kcal Atomic # 8 Polarization: -22.40 SS G_CDS: -0.21 Total: -22.62 kcal Atomic # 14 Polarization: 20.37 SS G_CDS: -5.02 Total: 15.35 kcal Total LS contribution 5.97 Total: 5.97 kcal Total: -26.90 -8.34 -35.23 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. ZINC000422488447.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 -74.416 kcal (2) G-P(sol) polarization free energy of solvation -26.896 kcal (3) G-ENP(sol) elect.-nuc.-pol. free energy of system -101.312 kcal (4) G-CDS(sol) cavity-dispersion-solvent structure free energy -8.336 kcal (5) G-P-CDS(sol) = G-P(sol) + G-CDS(sol) = (2) + (4) -35.232 kcal (6) G-S(sol) free energy of system = (1) + (5) -109.648 kcal 1SCF WAS SPECIFIED, SO GEOMETRY OPTIMIZATION WAS NOT USED Total computer time = 0.71 seconds