Wall clock time and date at job start Sat Apr 11 2020 06:11:40 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 O 1.42894 * 1 3 3 C 1.42903 * 113.99971 * 2 1 4 4 C 1.53000 * 109.47267 * 179.97438 * 3 2 1 5 5 O 1.42906 * 109.46972 * 65.00304 * 4 3 2 6 6 C 1.42895 * 114.00394 * 180.02562 * 5 4 3 7 7 N 1.46904 * 109.47495 * 79.99394 * 6 5 4 8 8 C 1.48047 * 111.00084 * 205.00638 * 7 6 5 9 9 C 1.53456 * 86.01993 * 222.73218 * 8 7 6 10 10 C 1.48048 * 111.00005 * 105.54230 * 7 6 5 11 11 C 1.52632 * 113.68952 * 268.90770 * 9 8 7 12 12 C 1.53187 * 109.35212 * 286.02589 * 11 9 8 13 13 C 1.53140 * 109.35167 * 300.91225 * 12 11 9 14 14 N 1.46899 * 109.08168 * 56.44783 * 13 12 11 15 15 C 1.34778 * 120.68711 * 124.67792 * 14 13 12 16 16 O 1.21277 * 119.99917 * 179.97438 * 15 14 13 17 17 C 1.50700 * 120.00097 * 359.97223 * 15 14 13 18 18 C 1.53001 * 109.47053 * 180.02562 * 17 15 14 19 19 C 1.50700 * 109.46887 * 180.02562 * 18 17 15 20 20 H 1.07994 * 120.00655 * 0.02562 * 19 18 17 21 21 C 1.37885 * 119.99983 * 179.97438 * 19 18 17 22 22 N 1.31514 * 119.99848 * 0.02562 * 21 19 18 23 23 Si 1.86796 * 120.00149 * 180.02562 * 21 19 18 24 24 H 1.48492 * 109.47129 * 89.99955 * 23 21 19 25 25 H 1.48498 * 109.47197 * 210.00628 * 23 21 19 26 26 H 1.48505 * 109.46831 * 330.00042 * 23 21 19 27 27 C 1.47117 * 118.62327 * 304.65290 * 14 13 12 28 28 H 1.08992 * 109.47170 * 299.99936 * 1 2 3 29 29 H 1.09006 * 109.47165 * 60.00365 * 1 2 3 30 30 H 1.09005 * 109.46914 * 179.97438 * 1 2 3 31 31 H 1.09005 * 109.47042 * 59.99831 * 3 2 1 32 32 H 1.08999 * 109.47064 * 299.99575 * 3 2 1 33 33 H 1.08996 * 109.47383 * 305.00349 * 4 3 2 34 34 H 1.08996 * 109.47437 * 184.99576 * 4 3 2 35 35 H 1.09009 * 109.46724 * 319.99687 * 6 5 4 36 36 H 1.09000 * 109.47474 * 199.99989 * 6 5 4 37 37 H 1.09005 * 113.78524 * 108.45883 * 8 7 6 38 38 H 1.08992 * 113.78636 * 337.00742 * 8 7 6 39 39 H 1.09003 * 113.78187 * 22.99038 * 10 7 6 40 40 H 1.08999 * 113.78712 * 251.55389 * 10 7 6 41 41 H 1.09007 * 109.49006 * 45.98113 * 11 9 8 42 42 H 1.08994 * 109.49333 * 166.05526 * 11 9 8 43 43 H 1.08995 * 109.48647 * 180.95056 * 12 11 9 44 44 H 1.08996 * 109.49089 * 60.87916 * 12 11 9 45 45 H 1.09005 * 109.53997 * 296.56592 * 13 12 11 46 46 H 1.09000 * 109.53313 * 176.32607 * 13 12 11 47 47 H 1.08995 * 109.47731 * 300.00456 * 17 15 14 48 48 H 1.09006 * 109.46850 * 60.00890 * 17 15 14 49 49 H 1.09006 * 109.47202 * 299.99899 * 18 17 15 50 50 H 1.08994 * 109.47491 * 60.00485 * 18 17 15 51 51 H 0.96991 * 120.00228 * 179.97438 * 22 21 19 52 52 H 1.08999 * 109.63882 * 291.72101 * 27 14 13 53 53 H 1.09001 * 109.80362 * 171.20766 * 27 14 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 8 1.4289 0.0000 0.0000 3 6 2.0102 1.3055 0.0000 4 6 3.5354 1.1847 0.0006 5 8 3.9658 0.5793 -1.2202 6 6 5.3821 0.4188 -1.3217 7 7 5.7960 -0.7690 -0.5628 8 6 7.2090 -0.6594 -0.1349 9 6 7.3373 -2.1311 -0.5500 10 6 6.1348 -1.8877 -1.4713 11 6 7.0540 -3.1146 0.5823 12 6 8.2346 -3.1187 1.5584 13 6 9.5095 -3.5205 0.8111 14 7 9.7134 -2.6033 -0.3182 15 6 10.8688 -1.9223 -0.4511 16 8 11.0241 -1.1687 -1.3886 17 6 11.9678 -2.0971 0.5651 18 6 13.1601 -1.2175 0.1833 19 6 14.2593 -1.3929 1.1992 20 1 14.1216 -2.0648 2.0334 21 6 15.4414 -0.6962 1.0631 22 7 15.6094 0.1217 0.0470 23 14 16.8042 -0.9141 2.3219 24 1 17.6971 -2.0222 1.8979 25 1 17.5889 0.3425 2.4244 26 1 16.2098 -1.2357 3.6443 27 6 8.6305 -2.4441 -1.3012 28 1 -0.3633 0.5138 0.8899 29 1 -0.3634 0.5138 -0.8901 30 1 -0.3633 -1.0277 -0.0005 31 1 1.6885 1.8464 0.8900 32 1 1.6886 1.8464 -0.8900 33 1 3.8508 0.5686 0.8426 34 1 3.9784 2.1766 0.0901 35 1 5.8777 1.3005 -0.9154 36 1 5.6597 0.2960 -2.3686 37 1 7.3327 -0.4920 0.9351 38 1 7.8014 0.0226 -0.7448 39 1 6.4089 -1.5695 -2.4772 40 1 5.4044 -2.6967 -1.4638 41 1 6.1483 -2.8124 1.1084 42 1 6.9192 -4.1146 0.1703 43 1 8.0420 -3.8328 2.3590 44 1 8.3599 -2.1222 1.9818 45 1 9.4079 -4.5400 0.4389 46 1 10.3623 -3.4613 1.4873 47 1 11.6022 -1.8056 1.5496 48 1 12.2790 -3.1416 0.5876 49 1 13.5260 -1.5087 -0.8013 50 1 12.8488 -0.1732 0.1611 51 1 16.4407 0.6120 -0.0485 52 1 8.5110 -3.3678 -1.8674 53 1 8.8665 -1.6255 -1.9812 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 ZINC001271252441.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:40 Heat of formation + Delta-G solvation = -94.472427 kcal Electronic energy + Delta-G solvation = -30309.409657 eV Core-core repulsion = 26182.350767 eV Total energy + Delta-G solvation = -4127.058890 eV No. of doubly occupied orbitals = 66 Molecular weight (most abundant/longest-lived isotopes) = 340.220 amu Computer time = 0.76 seconds Orbital eigenvalues (eV) -41.14770 -39.19739 -38.49224 -37.31409 -36.06841 -33.50504 -32.40981 -31.50037 -30.86560 -30.06759 -26.93463 -26.58886 -26.35967 -25.87630 -24.85590 -23.07982 -22.48954 -22.05663 -21.12730 -18.67184 -18.32359 -17.31159 -17.06311 -16.99148 -16.38510 -16.23004 -15.90847 -15.73756 -15.37370 -15.20932 -14.92875 -14.75218 -14.42151 -14.14032 -14.07883 -13.56036 -13.35702 -13.25810 -13.06227 -12.76712 -12.63809 -12.40612 -12.17996 -12.02032 -11.91109 -11.76014 -11.47494 -11.21357 -11.12706 -11.03711 -10.94405 -10.90534 -10.62643 -10.55141 -10.46637 -10.04573 -10.00045 -9.95769 -9.60125 -9.17500 -9.09261 -8.95646 -8.36301 -8.28896 -5.96169 -4.00473 3.20453 3.26797 3.32001 3.33071 3.39365 3.41894 3.45527 3.77568 4.18805 4.43901 4.47181 4.60715 4.72876 4.78560 4.85902 4.89735 4.97726 4.98469 5.00826 5.05429 5.11553 5.13845 5.17619 5.21854 5.27662 5.34124 5.38467 5.38809 5.47834 5.52662 5.54706 5.62305 5.70807 5.74528 5.77528 5.80711 5.95839 6.01682 6.17870 6.19615 6.22044 6.23764 6.69268 6.75476 6.81708 6.91724 7.00767 7.74336 7.77634 7.86345 8.07154 8.38635 8.75439 9.08029 9.63005 11.14782 Molecular weight = 340.22amu Principal moments of inertia in cm(-1) A = 0.019253 B = 0.001949 C = 0.001894 Principal moments of inertia in units of 10**(-40)*gram-cm**2 A = 1453.974454 B =14361.182065 C =14782.004878 Net atomic charges, atomic populations, and dipole contributions calculated with CM2 Atom NO. Type Charge No. of electrons 1 C 0.031 3.969 2 O -0.380 6.380 3 C 0.064 3.936 4 C 0.066 3.934 5 O -0.371 6.371 6 C 0.189 3.811 7 N -0.499 5.499 8 C 0.020 3.980 9 C -0.097 4.097 10 C 0.022 3.978 11 C -0.087 4.087 12 C -0.135 4.135 13 C 0.108 3.892 14 N -0.622 5.622 15 C 0.514 3.486 16 O -0.537 6.537 17 C -0.128 4.128 18 C 0.037 3.963 19 C -0.526 4.526 20 H 0.056 0.944 21 C 0.058 3.942 22 N -0.890 5.890 23 Si 0.892 3.108 24 H -0.281 1.281 25 H -0.287 1.287 26 H -0.274 1.274 27 C 0.125 3.875 28 H 0.039 0.961 29 H 0.040 0.960 30 H 0.085 0.915 31 H 0.062 0.938 32 H 0.055 0.945 33 H 0.066 0.934 34 H 0.059 0.941 35 H 0.060 0.940 36 H 0.084 0.916 37 H 0.090 0.910 38 H 0.074 0.926 39 H 0.073 0.927 40 H 0.091 0.909 41 H 0.075 0.925 42 H 0.065 0.935 43 H 0.075 0.925 44 H 0.068 0.932 45 H 0.064 0.936 46 H 0.091 0.909 47 H 0.079 0.921 48 H 0.078 0.922 49 H 0.023 0.977 50 H 0.023 0.977 51 H 0.261 0.739 52 H 0.064 0.936 53 H 0.090 0.910 For the charges calculated by CM2: Dipole moment (debyes) X Y Z Total from point charges -28.880 -3.282 0.899 29.080 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.066 4.066 2 O -0.299 6.299 3 C -0.013 4.013 4 C -0.011 4.011 5 O -0.286 6.286 6 C 0.023 3.977 7 N -0.228 5.228 8 C -0.108 4.108 9 C -0.099 4.099 10 C -0.107 4.107 11 C -0.125 4.125 12 C -0.173 4.173 13 C -0.015 4.015 14 N -0.358 5.358 15 C 0.304 3.696 16 O -0.412 6.412 17 C -0.168 4.168 18 C 0.000 4.000 19 C -0.565 4.565 20 H 0.074 0.926 21 C -0.144 4.144 22 N -0.528 5.528 23 Si 0.720 3.280 24 H -0.207 1.207 25 H -0.212 1.212 26 H -0.199 1.199 27 C 0.001 3.999 28 H 0.057 0.943 29 H 0.058 0.942 30 H 0.104 0.896 31 H 0.080 0.920 32 H 0.073 0.927 33 H 0.085 0.915 34 H 0.077 0.923 35 H 0.078 0.922 36 H 0.102 0.898 37 H 0.108 0.892 38 H 0.092 0.908 39 H 0.091 0.909 40 H 0.110 0.890 41 H 0.094 0.906 42 H 0.084 0.916 43 H 0.094 0.906 44 H 0.087 0.913 45 H 0.082 0.918 46 H 0.110 0.890 47 H 0.097 0.903 48 H 0.097 0.903 49 H 0.041 0.959 50 H 0.041 0.959 51 H 0.070 0.930 52 H 0.082 0.918 53 H 0.108 0.892 Dipole moment (debyes) X Y Z Total from point charges -29.224 -2.968 0.625 29.381 hybrid contribution 1.239 0.480 0.238 1.350 sum -27.985 -2.488 0.864 28.109 Atomic orbital electron populations 1.22740 0.81336 1.02806 0.99716 1.88007 1.19309 1.27225 1.95314 1.22567 0.92729 0.85455 1.00539 1.22495 0.92788 0.97061 0.88755 1.87994 1.21748 1.81738 1.37164 1.23806 0.82427 0.93531 0.97907 1.68176 1.20297 1.05732 1.28644 1.23931 0.88167 0.96243 1.02466 1.22324 0.94756 0.95842 0.96956 1.23979 0.96211 0.95558 0.94962 1.21188 0.98457 0.97923 0.94932 1.21930 0.94736 1.02598 0.98054 1.21675 0.98686 0.92711 0.88454 1.48279 1.15543 1.45439 1.26510 1.21567 0.82495 0.78400 0.87114 1.90648 1.80040 1.37013 1.33515 1.21874 0.95208 1.03064 0.96644 1.18903 0.89572 0.94530 0.97038 1.21220 0.99727 1.25655 1.09876 0.92624 1.25025 0.98507 0.95233 0.95613 1.69994 1.25846 1.27457 1.29538 0.86595 0.82125 0.77500 0.81751 1.20674 1.21243 1.19895 1.21446 0.84992 1.02199 0.91230 0.94267 0.94176 0.89615 0.92027 0.92714 0.91539 0.92256 0.92164 0.89792 0.89172 0.90825 0.90903 0.89027 0.90633 0.91591 0.90646 0.91326 0.91772 0.89031 0.90314 0.90319 0.95877 0.95901 0.92972 0.91826 0.89212 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.03 0.21 11.01 101.05 1.11 1.32 16 2 O -0.38 -4.00 10.67 -40.07 -0.43 -4.43 16 3 C 0.06 0.53 6.32 37.16 0.23 0.77 16 4 C 0.07 0.62 6.19 37.16 0.23 0.85 16 5 O -0.37 -4.31 10.54 -40.07 -0.42 -4.73 16 6 C 0.19 1.94 6.05 60.07 0.36 2.31 16 7 N -0.50 -5.33 6.21 -246.35 -1.53 -6.86 16 8 C 0.02 0.22 6.80 -2.96 -0.02 0.20 16 9 C -0.10 -1.00 0.87 -154.41 -0.13 -1.13 16 10 C 0.02 0.21 7.10 -2.97 -0.02 0.19 16 11 C -0.09 -0.80 5.02 -26.82 -0.13 -0.94 16 12 C -0.13 -1.26 6.00 -26.56 -0.16 -1.42 16 13 C 0.11 1.20 6.29 -3.75 -0.02 1.17 16 14 N -0.62 -8.95 2.98 -172.91 -0.52 -9.46 16 15 C 0.51 9.91 7.70 -10.98 -0.08 9.83 16 16 O -0.54 -12.54 15.42 5.56 0.09 -12.46 16 17 C -0.13 -2.60 3.98 -27.88 -0.11 -2.71 16 18 C 0.04 1.00 4.26 -27.88 -0.12 0.88 16 19 C -0.53 -14.89 9.11 -34.44 -0.31 -15.20 16 20 H 0.06 1.48 7.74 -52.49 -0.41 1.07 16 21 C 0.06 1.68 5.46 -17.82 -0.10 1.58 16 22 N -0.89 -27.40 13.29 55.43 0.74 -26.67 16 23 Si 0.89 21.62 29.54 -169.99 -5.02 16.60 16 24 H -0.28 -6.69 7.11 56.52 0.40 -6.29 16 25 H -0.29 -6.90 7.11 56.52 0.40 -6.50 16 26 H -0.27 -6.43 7.11 56.52 0.40 -6.03 16 27 C 0.12 1.55 5.60 -3.81 -0.02 1.53 16 28 H 0.04 0.21 8.14 -51.93 -0.42 -0.22 16 29 H 0.04 0.23 8.14 -51.93 -0.42 -0.19 16 30 H 0.09 0.50 8.14 -51.93 -0.42 0.07 16 31 H 0.06 0.38 8.14 -51.93 -0.42 -0.05 16 32 H 0.05 0.42 8.14 -51.93 -0.42 0.00 16 33 H 0.07 0.69 8.14 -51.93 -0.42 0.26 16 34 H 0.06 0.44 8.09 -51.93 -0.42 0.02 16 35 H 0.06 0.57 7.64 -51.92 -0.40 0.17 16 36 H 0.08 0.80 7.16 -51.93 -0.37 0.43 16 37 H 0.09 1.02 7.36 -51.93 -0.38 0.64 16 38 H 0.07 0.89 7.75 -51.93 -0.40 0.49 16 39 H 0.07 0.66 7.33 -51.93 -0.38 0.28 16 40 H 0.09 0.86 8.09 -51.93 -0.42 0.44 16 41 H 0.08 0.73 8.02 -51.93 -0.42 0.31 16 42 H 0.07 0.55 8.14 -51.93 -0.42 0.13 16 43 H 0.07 0.58 8.14 -51.93 -0.42 0.15 16 44 H 0.07 0.72 7.36 -51.93 -0.38 0.33 16 45 H 0.06 0.62 8.14 -51.93 -0.42 0.20 16 46 H 0.09 1.04 5.92 -51.93 -0.31 0.73 16 47 H 0.08 1.46 7.73 -51.93 -0.40 1.06 16 48 H 0.08 1.47 7.91 -51.93 -0.41 1.06 16 49 H 0.02 0.68 8.10 -51.93 -0.42 0.26 16 50 H 0.02 0.67 8.10 -51.93 -0.42 0.25 16 51 H 0.26 7.54 8.31 -40.82 -0.34 7.20 16 52 H 0.06 0.72 8.14 -51.93 -0.42 0.29 16 53 H 0.09 1.30 6.89 -51.93 -0.36 0.94 16 LS Contribution 418.65 15.07 6.31 6.31 Total: -1.00 -35.20 418.65 -9.74 -44.95 By element: Atomic # 1 Polarization: 7.19 SS G_CDS: -9.66 Total: -2.47 kcal Atomic # 6 Polarization: -1.48 SS G_CDS: 0.70 Total: -0.78 kcal Atomic # 7 Polarization: -41.68 SS G_CDS: -1.31 Total: -42.99 kcal Atomic # 8 Polarization: -20.85 SS G_CDS: -0.76 Total: -21.61 kcal Atomic # 14 Polarization: 21.62 SS G_CDS: -5.02 Total: 16.60 kcal Total LS contribution 6.31 Total: 6.31 kcal Total: -35.20 -9.74 -44.95 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. ZINC001271252441.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 -49.526 kcal (2) G-P(sol) polarization free energy of solvation -35.203 kcal (3) G-ENP(sol) elect.-nuc.-pol. free energy of system -84.729 kcal (4) G-CDS(sol) cavity-dispersion-solvent structure free energy -9.743 kcal (5) G-P-CDS(sol) = G-P(sol) + G-CDS(sol) = (2) + (4) -44.946 kcal (6) G-S(sol) free energy of system = (1) + (5) -94.472 kcal 1SCF WAS SPECIFIED, SO GEOMETRY OPTIMIZATION WAS NOT USED Total computer time = 0.76 seconds