Wall clock time and date at job start Sat Apr 11 2020 06:09:44 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.42896 * 1 3 3 C 1.42903 * 114.00046 * 2 1 4 4 C 1.53002 * 109.47235 * 179.97438 * 3 2 1 5 5 C 1.52998 * 109.45699 * 299.99509 * 4 3 2 6 6 N 1.46505 * 109.47327 * 60.03141 * 5 4 3 7 7 C 1.34777 * 120.00094 * 179.72619 * 6 5 4 8 8 O 1.21278 * 120.00008 * 0.26795 * 7 6 5 9 9 C 1.50704 * 119.99908 * 180.27726 * 7 6 5 10 10 H 1.08995 * 110.33095 * 55.69374 * 9 7 6 11 11 C 1.54479 * 110.36810 * 293.44345 * 9 7 6 12 12 C 1.54510 * 104.06325 * 217.19082 * 11 9 7 13 13 H 1.08998 * 110.50962 * 241.35821 * 12 11 9 14 14 C 1.53002 * 110.51082 * 118.78374 * 12 11 9 15 15 C 1.54478 * 104.06069 * 0.02562 * 12 11 9 16 16 O 1.44424 * 110.37295 * 177.94803 * 9 7 6 17 17 C 1.53039 * 109.45108 * 59.96103 * 4 3 2 18 18 C 1.53193 * 109.31237 * 181.37774 * 17 4 3 19 19 N 1.46924 * 108.77427 * 305.36637 * 18 17 4 20 20 C 1.36932 * 120.62687 * 233.58266 * 19 18 17 21 21 H 1.08000 * 119.99658 * 152.30961 * 20 19 18 22 22 C 1.38814 * 120.00648 * 332.31082 * 20 19 18 23 23 N 1.31621 * 119.99633 * 342.60253 * 22 20 19 24 24 Si 1.86795 * 120.00416 * 162.59760 * 22 20 19 25 25 H 1.48500 * 109.47103 * 94.34825 * 24 22 20 26 26 H 1.48502 * 109.47316 * 214.35256 * 24 22 20 27 27 H 1.48498 * 109.47050 * 334.35476 * 24 22 20 28 28 C 1.46918 * 118.73814 * 53.61183 * 19 18 17 29 29 C 1.53037 * 109.45819 * 180.02562 * 4 3 2 30 30 H 1.08995 * 109.47367 * 299.99858 * 1 2 3 31 31 H 1.09001 * 109.46815 * 60.00970 * 1 2 3 32 32 H 1.09007 * 109.46747 * 180.02562 * 1 2 3 33 33 H 1.08999 * 109.47227 * 59.99515 * 3 2 1 34 34 H 1.08999 * 109.46562 * 299.99521 * 3 2 1 35 35 H 1.09002 * 109.47466 * 300.02823 * 5 4 3 36 36 H 1.09000 * 109.47212 * 180.02562 * 5 4 3 37 37 H 0.97005 * 119.99545 * 359.97438 * 6 5 4 38 38 H 1.08999 * 110.51184 * 335.83822 * 11 9 7 39 39 H 1.09002 * 110.50994 * 98.41130 * 11 9 7 40 40 H 1.09002 * 109.47168 * 59.99983 * 14 12 11 41 41 H 1.09000 * 109.46894 * 180.02562 * 14 12 11 42 42 H 1.08996 * 109.46974 * 300.00227 * 14 12 11 43 43 H 1.08997 * 110.37272 * 142.80861 * 15 12 11 44 44 H 1.09007 * 110.37123 * 265.03526 * 15 12 11 45 45 H 1.09001 * 109.49775 * 301.33361 * 17 4 3 46 46 H 1.09002 * 109.49657 * 61.42323 * 17 4 3 47 47 H 1.08998 * 109.58099 * 185.57542 * 18 17 4 48 48 H 1.09005 * 109.58286 * 65.15816 * 18 17 4 49 49 H 0.96995 * 119.99771 * 179.97438 * 23 22 20 50 50 H 1.08998 * 109.58757 * 186.60137 * 28 19 18 51 51 H 1.08998 * 109.70846 * 66.24657 * 28 19 18 52 52 H 1.08998 * 109.49492 * 298.57433 * 29 4 3 53 53 H 1.09000 * 109.50044 * 58.63916 * 29 4 3 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.4290 0.0000 0.0000 3 6 2.0102 1.3055 0.0000 4 6 3.5355 1.1847 0.0006 5 6 3.9861 0.4257 1.2504 6 7 3.3906 -0.9129 1.2500 7 6 3.6535 -1.7676 2.2584 8 8 4.3791 -1.4269 3.1686 9 6 3.0460 -3.1468 2.2549 10 1 3.2985 -3.6694 1.3323 11 6 1.5127 -3.0593 2.4211 12 6 1.1616 -4.2728 3.3107 13 1 0.4938 -4.9532 2.7824 14 6 0.5247 -3.8096 4.6225 15 6 2.5245 -4.9495 3.5772 16 8 3.5013 -3.8997 3.4001 17 6 3.9870 0.4241 -1.2483 18 6 5.5115 0.2750 -1.2271 19 7 6.1149 1.6066 -1.0804 20 6 7.0612 2.0452 -1.9676 21 1 7.1780 3.1026 -2.1536 22 6 7.8704 1.1309 -2.6280 23 7 7.9800 -0.0995 -2.1734 24 14 8.8036 1.6464 -4.1618 25 1 10.1753 2.0739 -3.7864 26 1 8.8823 0.4998 -5.1022 27 1 8.0951 2.7755 -4.8161 28 6 5.6840 2.4582 0.0366 29 6 4.1576 2.5829 0.0013 30 1 -0.3634 0.5138 0.8899 31 1 -0.3633 0.5137 -0.8901 32 1 -0.3633 -1.0278 0.0005 33 1 1.6885 1.8464 0.8899 34 1 1.6886 1.8463 -0.8900 35 1 3.6635 0.9666 2.1400 36 1 5.0727 0.3399 1.2513 37 1 2.8101 -1.1852 0.5222 38 1 1.2347 -2.1292 2.9168 39 1 1.0185 -3.1373 1.4527 40 1 -0.3874 -3.2531 4.4069 41 1 0.2847 -4.6778 5.2363 42 1 1.2236 -3.1678 5.1589 43 1 2.5631 -5.3371 4.5952 44 1 2.6965 -5.7519 2.8597 45 1 3.5249 -0.5631 -1.2599 46 1 3.6877 0.9768 -2.1388 47 1 5.8472 -0.1777 -2.1601 48 1 5.8054 -0.3541 -0.3868 49 1 8.5452 -0.7385 -2.6350 50 1 6.1337 3.4461 -0.0626 51 1 5.9923 2.0096 0.9809 52 1 3.8568 3.1153 -0.9010 53 1 3.8166 3.1325 0.8786 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 ZINC001162366515.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:44 Heat of formation + Delta-G solvation = -106.074811 kcal Electronic energy + Delta-G solvation = -32477.012030 eV Core-core repulsion = 28349.450023 eV Total energy + Delta-G solvation = -4127.562008 eV No. of doubly occupied orbitals = 66 Molecular weight (most abundant/longest-lived isotopes) = 340.220 amu Computer time = 0.84 seconds Orbital eigenvalues (eV) -40.38525 -39.49265 -37.96272 -37.10224 -34.48366 -33.61307 -33.00033 -31.04258 -30.21089 -29.93319 -27.68190 -26.69888 -26.25940 -25.64342 -23.85427 -23.40190 -22.25003 -20.88476 -20.61842 -20.00407 -18.46618 -17.34723 -16.83508 -16.41614 -16.17774 -15.84886 -15.77763 -15.35105 -15.05413 -14.71744 -14.59935 -14.45562 -14.01464 -13.80086 -13.57856 -13.36389 -13.16994 -12.95245 -12.82244 -12.74125 -12.62450 -12.24042 -11.92910 -11.83834 -11.56304 -11.43246 -11.40118 -11.27113 -11.07225 -10.99144 -10.90588 -10.69384 -10.53414 -10.24206 -10.17609 -9.96159 -9.87076 -9.74443 -9.51637 -9.38926 -9.29575 -8.55108 -8.39990 -6.62883 -5.75330 -3.17374 2.92967 3.16123 3.38928 3.46286 3.48172 3.69023 4.05247 4.33639 4.51943 4.54085 4.68519 4.90868 4.95078 4.97131 5.02743 5.14407 5.21730 5.24387 5.28477 5.35801 5.40371 5.43148 5.45728 5.49383 5.52464 5.53347 5.57296 5.65627 5.68524 5.71965 5.72700 5.80674 5.99660 6.03634 6.07157 6.24759 6.38112 6.40708 6.61109 6.63464 6.65293 6.71946 6.82621 6.90319 7.00051 7.09453 7.40238 7.71368 7.82595 8.00166 8.47905 8.64176 9.24390 9.88033 10.44636 11.32886 Molecular weight = 340.22amu Principal moments of inertia in cm(-1) A = 0.014307 B = 0.003258 C = 0.002919 Principal moments of inertia in units of 10**(-40)*gram-cm**2 A = 1956.671753 B = 8591.462612 C = 9590.085304 Net atomic charges, atomic populations, and dipole contributions calculated with CM2 Atom NO. Type Charge No. of electrons 1 C 0.034 3.966 2 O -0.391 6.391 3 C 0.073 3.927 4 C -0.059 4.059 5 C 0.133 3.867 6 N -0.705 5.705 7 C 0.513 3.487 8 O -0.532 6.532 9 C 0.058 3.942 10 H 0.092 0.908 11 C -0.137 4.137 12 C -0.126 4.126 13 H 0.085 0.915 14 C -0.144 4.144 15 C 0.044 3.956 16 O -0.342 6.342 17 C -0.125 4.125 18 C 0.170 3.830 19 N -0.591 5.591 20 C -0.253 4.253 21 H 0.065 0.935 22 C 0.002 3.998 23 N -0.900 5.900 24 Si 0.900 3.100 25 H -0.288 1.288 26 H -0.290 1.290 27 H -0.280 1.280 28 C 0.143 3.857 29 C -0.119 4.119 30 H 0.037 0.963 31 H 0.040 0.960 32 H 0.082 0.918 33 H 0.046 0.954 34 H 0.051 0.949 35 H 0.055 0.945 36 H 0.084 0.916 37 H 0.408 0.592 38 H 0.081 0.919 39 H 0.085 0.915 40 H 0.056 0.944 41 H 0.056 0.944 42 H 0.062 0.938 43 H 0.094 0.906 44 H 0.061 0.939 45 H 0.054 0.946 46 H 0.056 0.944 47 H 0.086 0.914 48 H 0.022 0.978 49 H 0.253 0.747 50 H 0.056 0.944 51 H 0.022 0.978 52 H 0.066 0.934 53 H 0.054 0.946 For the charges calculated by CM2: Dipole moment (debyes) X Y Z Total from point charges -20.440 -5.517 7.119 22.337 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.063 4.063 2 O -0.310 6.310 3 C -0.004 4.004 4 C -0.059 4.059 5 C 0.010 3.990 6 N -0.363 5.363 7 C 0.297 3.703 8 O -0.408 6.408 9 C -0.003 4.003 10 H 0.110 0.890 11 C -0.175 4.175 12 C -0.144 4.144 13 H 0.103 0.897 14 C -0.201 4.201 15 C -0.033 4.033 16 O -0.260 6.260 17 C -0.164 4.164 18 C 0.044 3.956 19 N -0.314 5.314 20 C -0.382 4.382 21 H 0.083 0.917 22 C -0.196 4.196 23 N -0.544 5.544 24 Si 0.730 3.270 25 H -0.215 1.215 26 H -0.216 1.216 27 H -0.205 1.205 28 C 0.016 3.984 29 C -0.159 4.159 30 H 0.056 0.944 31 H 0.058 0.942 32 H 0.100 0.900 33 H 0.064 0.936 34 H 0.069 0.931 35 H 0.073 0.927 36 H 0.102 0.898 37 H 0.249 0.751 38 H 0.099 0.901 39 H 0.103 0.897 40 H 0.076 0.924 41 H 0.076 0.924 42 H 0.081 0.919 43 H 0.112 0.888 44 H 0.080 0.920 45 H 0.072 0.928 46 H 0.075 0.925 47 H 0.104 0.896 48 H 0.040 0.960 49 H 0.063 0.937 50 H 0.074 0.926 51 H 0.040 0.960 52 H 0.084 0.916 53 H 0.073 0.927 Dipole moment (debyes) X Y Z Total from point charges -19.498 -5.519 7.820 21.721 hybrid contribution -0.236 0.225 -1.320 1.360 sum -19.734 -5.294 6.500 21.441 Atomic orbital electron populations 1.22694 0.81486 1.02335 0.99808 1.87560 1.19200 1.28868 1.95325 1.22768 0.93932 0.83867 0.99821 1.20807 0.91393 0.97289 0.96451 1.21448 0.99543 0.81656 0.96339 1.46113 1.50359 1.14604 1.25181 1.20057 0.81583 0.85952 0.82741 1.90707 1.36746 1.77367 1.35957 1.22846 0.93787 0.92002 0.91662 0.88995 1.22974 0.93987 0.99506 1.01006 1.22530 0.97707 0.97029 0.97173 0.89651 1.21832 1.00612 1.01403 0.96228 1.23571 0.86792 0.90680 1.02286 1.89308 1.46427 1.42882 1.47391 1.22179 1.00892 0.98348 0.95026 1.20814 0.86343 0.86732 1.01696 1.44043 1.44711 1.14780 1.27817 1.19159 1.10419 0.92250 1.16407 0.91735 1.25034 0.97575 0.96651 1.00380 1.70047 1.40696 1.04569 1.39093 0.86319 0.79819 0.77510 0.83340 1.21510 1.21624 1.20518 1.20926 0.90576 0.95470 0.91413 1.21986 0.99350 0.93513 1.01018 0.94394 0.94169 0.89980 0.93598 0.93104 0.92668 0.89756 0.75052 0.90052 0.89663 0.92443 0.92448 0.91929 0.88787 0.92048 0.92760 0.92532 0.89573 0.96006 0.93678 0.92586 0.95988 0.91563 0.92714 Number of geometries 1 Number of calculations of the screened coulomb radii 1 The total number of SCF iterations 17. 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.27 11.01 101.05 1.11 1.39 16 2 O -0.39 -4.49 6.73 -35.23 -0.24 -4.73 16 3 C 0.07 0.90 4.80 37.16 0.18 1.08 16 4 C -0.06 -0.91 0.51 -154.47 -0.08 -0.99 16 5 C 0.13 2.28 3.80 -4.04 -0.02 2.26 16 6 N -0.70 -11.00 4.62 -60.29 -0.28 -11.28 16 7 C 0.51 8.82 7.07 -10.98 -0.08 8.74 16 8 O -0.53 -11.59 16.55 -13.65 -0.23 -11.82 16 9 C 0.06 0.79 3.71 -27.11 -0.10 0.69 16 10 H 0.09 1.10 8.14 -51.93 -0.42 0.68 16 11 C -0.14 -1.33 6.27 -25.16 -0.16 -1.49 16 12 C -0.13 -1.04 3.71 -89.05 -0.33 -1.37 16 13 H 0.09 0.54 8.14 -51.93 -0.42 0.12 16 14 C -0.14 -1.13 8.95 37.16 0.33 -0.80 16 15 C 0.04 0.44 7.18 37.93 0.27 0.71 16 16 O -0.34 -5.18 10.92 -54.44 -0.59 -5.78 16 17 C -0.13 -2.26 4.59 -26.61 -0.12 -2.38 16 18 C 0.17 3.97 4.84 -3.71 -0.02 3.95 16 19 N -0.59 -14.51 2.99 -172.33 -0.52 -15.02 16 20 C -0.25 -6.88 9.58 -15.06 -0.14 -7.02 16 21 H 0.06 1.72 7.89 -52.49 -0.41 1.31 16 22 C 0.00 0.06 4.91 -17.43 -0.09 -0.02 16 23 N -0.90 -26.28 11.19 55.49 0.62 -25.66 16 24 Si 0.90 21.69 29.59 -169.99 -5.03 16.67 16 25 H -0.29 -6.99 7.11 56.52 0.40 -6.59 16 26 H -0.29 -7.03 7.11 56.52 0.40 -6.63 16 27 H -0.28 -6.59 7.11 56.52 0.40 -6.19 16 28 C 0.14 2.92 6.04 -3.71 -0.02 2.89 16 29 C -0.12 -1.92 4.91 -26.61 -0.13 -2.05 16 30 H 0.04 0.29 8.14 -51.93 -0.42 -0.14 16 31 H 0.04 0.31 8.14 -51.93 -0.42 -0.12 16 32 H 0.08 0.53 8.14 -51.92 -0.42 0.11 16 33 H 0.05 0.53 8.14 -51.93 -0.42 0.11 16 34 H 0.05 0.60 8.14 -51.93 -0.42 0.18 16 35 H 0.06 0.92 8.14 -51.93 -0.42 0.50 16 36 H 0.08 1.76 4.98 -51.93 -0.26 1.50 16 37 H 0.41 5.67 4.71 -40.82 -0.19 5.48 16 38 H 0.08 0.87 7.46 -51.93 -0.39 0.48 16 39 H 0.08 0.65 8.14 -51.93 -0.42 0.22 16 40 H 0.06 0.40 8.14 -51.93 -0.42 -0.02 16 41 H 0.06 0.37 8.14 -51.93 -0.42 -0.05 16 42 H 0.06 0.61 8.14 -51.93 -0.42 0.19 16 43 H 0.09 0.82 7.87 -51.93 -0.41 0.42 16 44 H 0.06 0.52 8.14 -51.93 -0.42 0.10 16 45 H 0.05 0.91 6.50 -51.93 -0.34 0.57 16 46 H 0.06 1.05 8.14 -51.93 -0.42 0.62 16 47 H 0.09 2.29 6.04 -51.93 -0.31 1.97 16 48 H 0.02 0.56 6.69 -51.93 -0.35 0.21 16 49 H 0.25 7.17 8.31 -40.82 -0.34 6.84 16 50 H 0.06 1.12 7.93 -51.93 -0.41 0.71 16 51 H 0.02 0.47 6.70 -51.93 -0.35 0.12 16 52 H 0.07 1.07 8.14 -51.93 -0.42 0.65 16 53 H 0.05 0.75 8.14 -51.93 -0.42 0.33 16 LS Contribution 401.16 15.07 6.05 6.05 Total: -1.00 -33.39 401.16 -8.92 -42.31 By element: Atomic # 1 Polarization: 13.01 SS G_CDS: -9.32 Total: 3.69 kcal Atomic # 6 Polarization: 4.96 SS G_CDS: 0.61 Total: 5.58 kcal Atomic # 7 Polarization: -51.78 SS G_CDS: -0.17 Total: -51.96 kcal Atomic # 8 Polarization: -21.27 SS G_CDS: -1.06 Total: -22.32 kcal Atomic # 14 Polarization: 21.69 SS G_CDS: -5.03 Total: 16.67 kcal Total LS contribution 6.05 Total: 6.05 kcal Total: -33.39 -8.92 -42.31 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. ZINC001162366515.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 -63.768 kcal (2) G-P(sol) polarization free energy of solvation -33.385 kcal (3) G-ENP(sol) elect.-nuc.-pol. free energy of system -97.153 kcal (4) G-CDS(sol) cavity-dispersion-solvent structure free energy -8.922 kcal (5) G-P-CDS(sol) = G-P(sol) + G-CDS(sol) = (2) + (4) -42.307 kcal (6) G-S(sol) free energy of system = (1) + (5) -106.075 kcal 1SCF WAS SPECIFIED, SO GEOMETRY OPTIMIZATION WAS NOT USED Total computer time = 0.85 seconds