Wall clock time and date at job start Tue Mar 31 2020 21:12:54 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 N 2 2 C 1.31619 * 1 3 3 Si 1.86800 * 119.99953 * 2 1 4 4 H 1.48498 * 109.47091 * 240.00055 * 3 2 1 5 5 H 1.48509 * 109.47096 * 0.02562 * 3 2 1 6 6 H 1.48497 * 109.47112 * 119.99679 * 3 2 1 7 7 C 1.38812 * 120.00108 * 179.97438 * 2 1 3 8 8 H 1.07997 * 120.00167 * 179.97438 * 7 2 1 9 9 N 1.36936 * 119.99717 * 0.02562 * 7 2 1 10 10 C 1.47427 * 125.65054 * 359.94606 * 9 7 2 11 11 C 1.54913 * 104.82971 * 204.17546 * 10 9 7 12 12 C 1.55161 * 101.58474 * 322.98963 * 11 10 9 13 13 H 1.08993 * 110.71606 * 277.14271 * 12 11 10 14 14 C 1.52998 * 110.72216 * 154.04378 * 12 11 10 15 15 H 1.08994 * 109.58258 * 56.28373 * 14 12 11 16 16 C 1.53197 * 109.58413 * 295.86100 * 14 12 11 17 17 C 1.53042 * 109.30994 * 185.57786 * 16 14 12 18 18 C 1.53046 * 109.53665 * 61.36378 * 17 16 14 19 19 C 1.53186 * 109.31374 * 298.63400 * 18 17 16 20 20 N 1.46924 * 108.77880 * 54.63579 * 19 18 17 21 21 C 1.34780 * 120.63320 * 126.40892 * 20 19 18 22 22 O 1.21285 * 119.99769 * 184.36265 * 21 20 19 23 23 C 1.50700 * 119.99601 * 4.36746 * 21 20 19 24 24 C 1.52998 * 115.55035 * 284.28308 * 23 21 20 25 25 C 1.53003 * 109.47071 * 303.73160 * 24 23 21 26 26 F 1.39901 * 109.47378 * 185.09066 * 25 24 23 27 27 C 1.52998 * 117.49662 * 138.61621 * 23 21 20 28 28 C 1.52999 * 117.50098 * 69.95823 * 23 21 20 29 29 C 1.47024 * 125.64829 * 180.02562 * 9 7 2 30 30 H 0.97004 * 120.00466 * 359.97438 * 1 2 3 31 31 H 1.08994 * 110.36547 * 85.33862 * 10 9 7 32 32 H 1.08997 * 110.36548 * 323.01004 * 10 9 7 33 33 H 1.09002 * 111.00231 * 204.91598 * 11 10 9 34 34 H 1.08998 * 111.00648 * 81.06304 * 11 10 9 35 35 H 1.08998 * 109.49504 * 305.52872 * 16 14 12 36 36 H 1.08997 * 109.49869 * 65.62823 * 16 14 12 37 37 H 1.09004 * 109.45959 * 301.34796 * 17 16 14 38 38 H 1.09000 * 109.46292 * 181.38368 * 17 16 14 39 39 H 1.09004 * 109.49150 * 178.68177 * 18 17 16 40 40 H 1.09000 * 109.49784 * 58.59097 * 18 17 16 41 41 H 1.09001 * 109.59062 * 294.85228 * 19 18 17 42 42 H 1.08998 * 109.70411 * 174.49506 * 19 18 17 43 43 H 1.09000 * 109.47005 * 63.72610 * 24 23 21 44 44 H 1.08993 * 109.47113 * 183.73280 * 24 23 21 45 45 H 1.09004 * 109.47147 * 305.08825 * 25 24 23 46 46 H 1.09003 * 109.47431 * 65.09310 * 25 24 23 47 47 H 1.09003 * 117.49972 * 0.02562 * 27 23 21 48 48 H 1.08999 * 117.49901 * 145.01344 * 27 23 21 49 49 H 1.09000 * 117.49814 * 214.97957 * 28 23 21 50 50 H 1.09002 * 117.49853 * 359.96970 * 28 23 21 51 51 H 1.09001 * 109.88204 * 59.46875 * 29 9 7 52 52 H 1.09001 * 110.00934 * 298.38629 * 29 9 7 Note: An asterisk (*) indicates this parameter was or will be optimized. Cartesian coordinates (angstroms) NO. Atom X Y Z 1 7 0.0000 0.0000 0.0000 2 6 1.3162 0.0000 0.0000 3 14 2.2502 1.6177 0.0000 4 1 3.1039 1.6964 -1.2125 5 1 1.2851 2.7465 0.0006 6 1 3.1038 1.6964 1.2125 7 6 2.0103 -1.2021 0.0005 8 1 3.0902 -1.2021 0.0001 9 7 1.3256 -2.3880 0.0005 10 6 -0.1416 -2.5331 0.0011 11 6 -0.4014 -3.9319 0.6142 12 6 0.7785 -4.7525 0.0293 13 1 0.5455 -5.0987 -0.9776 14 6 1.1207 -5.9332 0.9401 15 1 0.2340 -6.5496 1.0878 16 6 1.6181 -5.4116 2.2920 17 6 2.0893 -6.5903 3.1469 18 6 3.2504 -7.2993 2.4457 19 6 2.7714 -7.8422 1.0957 20 7 2.1837 -6.7388 0.3238 21 6 2.6112 -6.4708 -0.9260 22 8 2.1749 -5.5096 -1.5231 23 6 3.6314 -7.3658 -1.5814 24 6 5.0705 -7.1961 -1.0904 25 6 5.4986 -5.7371 -1.2602 26 9 6.8493 -5.6048 -0.9208 27 6 3.4246 -7.7098 -3.0578 28 6 3.1831 -8.7724 -1.9828 29 6 1.9317 -3.7275 0.0016 30 1 -0.4851 0.8400 0.0004 31 1 -0.5276 -2.4876 -1.0171 32 1 -0.5991 -1.7584 0.6163 33 1 -1.3589 -4.3326 0.2813 34 1 -0.3513 -3.8969 1.7024 35 1 2.4477 -4.7228 2.1331 36 1 0.8068 -4.8932 2.8028 37 1 1.2659 -7.2912 3.2845 38 1 2.4209 -6.2245 4.1186 39 1 3.6007 -8.1243 3.0660 40 1 4.0647 -6.5926 2.2855 41 1 2.0202 -8.6150 1.2583 42 1 3.6152 -8.2622 0.5482 43 1 5.1300 -7.4720 -0.0375 44 1 5.7310 -7.8387 -1.6723 45 1 5.3533 -5.4341 -2.2972 46 1 4.8966 -5.1035 -0.6088 47 1 2.5509 -7.2885 -3.5552 48 1 4.3138 -7.8156 -3.6792 49 1 3.9135 -9.5770 -1.8974 50 1 2.1504 -9.0508 -1.7725 51 1 2.5583 -3.8501 0.8851 52 1 2.5298 -3.8669 -0.8989 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) F: (AM1): M.J.S. DEWAR AND E. G. ZOEBISCH, THEOCHEM, 180, 1 (1988). 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 ZINC001670924057.mol2 52 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 21:12:54 Heat of formation + Delta-G solvation = -7.093416 kcal Electronic energy + Delta-G solvation = -32626.504927 eV Core-core repulsion = 28542.919343 eV Total energy + Delta-G solvation = -4083.585584 eV No. of doubly occupied orbitals = 65 Molecular weight (most abundant/longest-lived isotopes) = 338.211 amu Computer time = 0.90 seconds Orbital eigenvalues (eV) -47.59971 -41.02834 -38.43203 -36.92731 -36.39211 -32.98327 -32.72648 -31.45953 -30.84014 -28.90829 -27.23573 -25.48350 -24.46423 -24.26627 -23.17215 -22.75564 -22.07935 -20.67288 -19.92877 -19.46954 -17.98804 -16.95811 -16.60346 -16.24862 -16.09149 -15.73882 -15.31509 -15.07519 -14.88288 -14.59547 -14.39292 -14.29230 -14.10410 -13.71029 -13.55310 -13.11672 -13.02842 -12.84692 -12.39012 -12.26247 -11.85769 -11.77742 -11.64270 -11.57289 -11.54283 -11.18859 -10.98421 -10.90277 -10.81640 -10.72015 -10.69939 -10.54164 -10.40040 -10.21987 -10.01246 -9.90239 -9.73764 -9.52897 -9.34776 -8.73336 -8.42953 -8.23837 -6.56269 -5.48923 -2.87664 2.68356 2.96489 3.41019 3.49003 3.52718 3.54326 4.16272 4.36279 4.51613 4.55321 4.63199 4.68933 4.82177 4.86786 4.88128 4.98434 5.10709 5.14260 5.23340 5.28584 5.31082 5.34671 5.45248 5.52246 5.60676 5.63170 5.65337 5.67925 5.73610 5.85027 5.87626 5.89480 5.91251 5.93118 5.94880 6.10168 6.14193 6.24454 6.40800 6.48134 6.57892 6.72889 6.89641 6.99193 7.17507 7.43314 7.55341 7.71818 7.79907 8.02984 8.43162 8.55750 9.40155 10.02643 10.70471 11.54797 Molecular weight = 338.21amu Principal moments of inertia in cm(-1) A = 0.011345 B = 0.004279 C = 0.003785 Principal moments of inertia in units of 10**(-40)*gram-cm**2 A = 2467.372944 B = 6542.558055 C = 7395.375888 Net atomic charges, atomic populations, and dipole contributions calculated with CM2 Atom NO. Type Charge No. of electrons 1 N -0.925 5.925 2 C -0.015 4.015 3 Si 0.910 3.090 4 H -0.287 1.287 5 H -0.291 1.291 6 H -0.287 1.287 7 C -0.212 4.212 8 H 0.072 0.928 9 N -0.604 5.604 10 C 0.157 3.843 11 C -0.150 4.150 12 C -0.094 4.094 13 H 0.102 0.898 14 C 0.161 3.839 15 H 0.067 0.933 16 C -0.141 4.141 17 C -0.118 4.118 18 C -0.136 4.136 19 C 0.107 3.893 20 N -0.602 5.602 21 C 0.530 3.470 22 O -0.490 6.490 23 C -0.156 4.156 24 C -0.107 4.107 25 C 0.060 3.940 26 F -0.220 7.220 27 C -0.143 4.143 28 C -0.159 4.159 29 C 0.144 3.856 30 H 0.250 0.750 31 H 0.007 0.993 32 H 0.102 0.898 33 H 0.051 0.949 34 H 0.057 0.943 35 H 0.083 0.917 36 H 0.077 0.923 37 H 0.059 0.941 38 H 0.063 0.937 39 H 0.071 0.929 40 H 0.067 0.933 41 H 0.064 0.936 42 H 0.086 0.914 43 H 0.082 0.918 44 H 0.082 0.918 45 H 0.080 0.920 46 H 0.078 0.922 47 H 0.105 0.895 48 H 0.102 0.898 49 H 0.102 0.898 50 H 0.099 0.901 51 H 0.014 0.986 52 H 0.047 0.953 For the charges calculated by CM2: Dipole moment (debyes) X Y Z Total from point charges 6.057 -22.846 1.658 23.694 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 N -0.574 5.574 2 C -0.210 4.210 3 Si 0.741 3.259 4 H -0.214 1.214 5 H -0.217 1.217 6 H -0.214 1.214 7 C -0.341 4.341 8 H 0.090 0.910 9 N -0.330 5.330 10 C 0.032 3.968 11 C -0.189 4.189 12 C -0.113 4.113 13 H 0.120 0.880 14 C 0.058 3.942 15 H 0.085 0.915 16 C -0.178 4.178 17 C -0.156 4.156 18 C -0.174 4.174 19 C -0.016 4.016 20 N -0.336 5.336 21 C 0.317 3.683 22 O -0.365 6.365 23 C -0.159 4.159 24 C -0.144 4.144 25 C 0.005 3.995 26 F -0.201 7.201 27 C -0.180 4.180 28 C -0.196 4.196 29 C 0.019 3.981 30 H 0.060 0.940 31 H 0.025 0.975 32 H 0.119 0.881 33 H 0.070 0.930 34 H 0.075 0.925 35 H 0.101 0.899 36 H 0.096 0.904 37 H 0.078 0.922 38 H 0.081 0.919 39 H 0.089 0.911 40 H 0.086 0.914 41 H 0.082 0.918 42 H 0.105 0.895 43 H 0.100 0.900 44 H 0.101 0.899 45 H 0.098 0.902 46 H 0.097 0.903 47 H 0.123 0.877 48 H 0.121 0.879 49 H 0.120 0.880 50 H 0.117 0.883 51 H 0.031 0.969 52 H 0.064 0.936 Dipole moment (debyes) X Y Z Total from point charges 5.811 -22.082 1.333 22.873 hybrid contribution 0.549 0.778 0.128 0.961 sum 6.360 -21.304 1.462 22.281 Atomic orbital electron populations 1.69804 1.04786 1.24448 1.58403 1.24733 0.95655 0.97437 1.03164 0.86216 0.78402 0.83516 0.77791 1.21392 1.21673 1.21373 1.19017 0.90305 0.82812 1.41966 0.90985 1.44527 1.01918 1.00059 1.86477 1.21818 0.84964 0.92229 0.97764 1.23225 0.96336 0.99628 0.99731 1.22693 0.96871 0.92996 0.98765 0.88005 1.21286 0.87933 0.92033 0.92913 0.91471 1.22196 1.02945 0.97799 0.94898 1.21524 0.97660 0.97681 0.98719 1.21932 0.98923 1.00599 0.95977 1.21763 0.98424 0.88910 0.92493 1.48000 1.38508 1.32933 1.14169 1.19427 0.82358 0.84668 0.81805 1.90826 1.55488 1.26119 1.64071 1.21353 0.95354 0.99132 1.00102 1.21708 0.92683 0.97255 1.02795 1.21975 0.77370 0.97182 1.02930 1.92751 1.35905 1.97072 1.94383 1.23048 1.04243 1.00900 0.89827 1.23094 1.02774 0.90753 1.03012 1.21843 0.93905 0.84691 0.97668 0.93974 0.97541 0.88064 0.93024 0.92480 0.89871 0.90418 0.92199 0.91854 0.91073 0.91431 0.91767 0.89548 0.89994 0.89890 0.90158 0.90324 0.87693 0.87941 0.87977 0.88272 0.96857 0.93608 Number of geometries 1 Number of calculations of the screened coulomb radii 1 The total number of SCF iterations 15. 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 N -0.92 -28.24 10.09 55.49 0.56 -27.68 16 2 C -0.02 -0.45 5.49 -17.43 -0.10 -0.54 16 3 Si 0.91 23.16 29.55 -169.99 -5.02 18.14 16 4 H -0.29 -7.25 7.11 56.52 0.40 -6.85 16 5 H -0.29 -7.40 7.11 56.52 0.40 -7.00 16 6 H -0.29 -7.17 7.11 56.52 0.40 -6.76 16 7 C -0.21 -6.28 10.27 -15.06 -0.15 -6.43 16 8 H 0.07 2.08 7.83 -52.49 -0.41 1.67 16 9 N -0.60 -16.54 3.39 -170.07 -0.58 -17.12 16 10 C 0.16 4.10 5.89 -2.53 -0.01 4.09 16 11 C -0.15 -3.03 6.43 -24.58 -0.16 -3.19 16 12 C -0.09 -1.79 2.06 -88.79 -0.18 -1.97 16 13 H 0.10 2.11 5.90 -52.31 -0.31 1.80 16 14 C 0.16 2.33 2.76 -67.60 -0.19 2.15 16 15 H 0.07 0.90 8.14 -51.93 -0.42 0.48 16 16 C -0.14 -1.97 4.79 -26.61 -0.13 -2.09 16 17 C -0.12 -1.22 5.88 -26.69 -0.16 -1.38 16 18 C -0.14 -1.22 6.08 -26.61 -0.16 -1.38 16 19 C 0.11 0.96 5.94 -3.71 -0.02 0.94 16 20 N -0.60 -7.52 2.96 -165.22 -0.49 -8.01 16 21 C 0.53 7.36 5.62 -10.98 -0.06 7.30 16 22 O -0.49 -8.96 9.92 5.55 0.06 -8.91 16 23 C -0.16 -1.59 0.83 -155.66 -0.13 -1.72 16 24 C -0.11 -1.01 5.32 -26.73 -0.14 -1.15 16 25 C 0.06 0.75 7.23 37.16 0.27 1.02 16 26 F -0.22 -3.04 17.36 2.25 0.04 -3.00 16 27 C -0.14 -1.33 9.81 -26.70 -0.26 -1.59 16 28 C -0.16 -1.22 9.20 -26.70 -0.25 -1.46 16 29 C 0.14 3.26 5.63 -3.07 -0.02 3.25 16 30 H 0.25 7.40 8.31 -40.82 -0.34 7.06 16 31 H 0.01 0.19 8.14 -51.93 -0.42 -0.24 16 32 H 0.10 2.92 5.87 -51.93 -0.30 2.62 16 33 H 0.05 0.94 8.14 -51.93 -0.42 0.52 16 34 H 0.06 1.10 6.26 -51.93 -0.33 0.78 16 35 H 0.08 1.39 5.51 -51.93 -0.29 1.10 16 36 H 0.08 1.14 6.57 -51.93 -0.34 0.80 16 37 H 0.06 0.56 8.14 -51.93 -0.42 0.14 16 38 H 0.06 0.63 8.14 -51.93 -0.42 0.21 16 39 H 0.07 0.49 8.14 -51.93 -0.42 0.07 16 40 H 0.07 0.70 8.14 -51.93 -0.42 0.27 16 41 H 0.06 0.50 8.14 -51.93 -0.42 0.08 16 42 H 0.09 0.66 4.25 -51.93 -0.22 0.43 16 43 H 0.08 0.71 6.41 -51.93 -0.33 0.38 16 44 H 0.08 0.65 8.10 -51.93 -0.42 0.22 16 45 H 0.08 0.99 8.14 -51.93 -0.42 0.56 16 46 H 0.08 1.14 7.66 -51.93 -0.40 0.75 16 47 H 0.10 1.14 8.14 -51.93 -0.42 0.72 16 48 H 0.10 0.79 8.14 -51.93 -0.42 0.37 16 49 H 0.10 0.58 8.14 -51.93 -0.42 0.15 16 50 H 0.10 0.78 8.14 -51.93 -0.42 0.36 16 51 H 0.01 0.31 5.12 -51.93 -0.27 0.05 16 52 H 0.05 1.13 5.64 -51.96 -0.29 0.84 16 LS Contribution 383.12 15.07 5.77 5.77 Total: -1.00 -33.39 383.12 -10.05 -43.45 By element: Atomic # 1 Polarization: 10.09 SS G_CDS: -8.54 Total: 1.55 kcal Atomic # 6 Polarization: -2.33 SS G_CDS: -1.85 Total: -4.18 kcal Atomic # 7 Polarization: -52.30 SS G_CDS: -0.51 Total: -52.81 kcal Atomic # 8 Polarization: -8.96 SS G_CDS: 0.06 Total: -8.91 kcal Atomic # 9 Polarization: -3.04 SS G_CDS: 0.04 Total: -3.00 kcal Atomic # 14 Polarization: 23.16 SS G_CDS: -5.02 Total: 18.14 kcal Total LS contribution 5.77 Total: 5.77 kcal Total: -33.39 -10.05 -43.45 kcal The number of atoms in the molecule is 52 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. ZINC001670924057.mol2 52 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 36.352 kcal (2) G-P(sol) polarization free energy of solvation -33.393 kcal (3) G-ENP(sol) elect.-nuc.-pol. free energy of system 2.959 kcal (4) G-CDS(sol) cavity-dispersion-solvent structure free energy -10.052 kcal (5) G-P-CDS(sol) = G-P(sol) + G-CDS(sol) = (2) + (4) -43.445 kcal (6) G-S(sol) free energy of system = (1) + (5) -7.093 kcal 1SCF WAS SPECIFIED, SO GEOMETRY OPTIMIZATION WAS NOT USED Total computer time = 0.91 seconds