Wall clock time and date at job start Sat Apr 11 2020 06:10:02 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.53000 * 1 3 3 H 1.09005 * 109.46960 * 2 1 4 4 C 1.52997 * 109.47497 * 240.00503 * 2 1 3 5 5 N 1.46503 * 109.47420 * 64.99871 * 4 2 1 6 6 C 1.46500 * 119.99744 * 270.00490 * 5 4 2 7 7 C 1.34769 * 120.00082 * 89.99966 * 5 4 2 8 8 O 1.21583 * 120.00289 * 3.83399 * 7 5 4 9 9 N 1.34783 * 119.99936 * 183.83903 * 7 5 4 10 10 C 1.46497 * 119.99807 * 185.44543 * 9 7 5 11 11 C 1.53004 * 109.47191 * 185.11363 * 10 9 7 12 12 C 1.50697 * 115.54786 * 57.06919 * 11 10 9 13 13 H 1.08002 * 119.99936 * 270.64950 * 12 11 10 14 14 C 1.37882 * 120.00032 * 90.37480 * 12 11 10 15 15 N 1.31512 * 119.99934 * 172.35786 * 14 12 11 16 16 Si 1.86799 * 119.99786 * 352.35242 * 14 12 11 17 17 H 1.48504 * 109.46870 * 93.56532 * 16 14 12 18 18 H 1.48499 * 109.47270 * 213.56104 * 16 14 12 19 19 H 1.48496 * 109.47387 * 333.56916 * 16 14 12 20 20 C 1.52995 * 117.50097 * 202.73932 * 11 10 9 21 21 C 1.52998 * 117.49881 * 271.40253 * 11 10 9 22 22 C 1.50697 * 109.47128 * 119.99915 * 2 1 3 23 23 O 1.20831 * 119.99779 * 120.00553 * 22 2 1 24 24 O 1.34222 * 120.00220 * 300.00387 * 22 2 1 25 25 C 1.45200 * 117.00189 * 179.97438 * 24 22 2 26 26 C 1.53003 * 109.47388 * 59.99883 * 25 24 22 27 27 C 1.52994 * 109.47482 * 180.02562 * 25 24 22 28 28 C 1.53001 * 109.47152 * 299.99961 * 25 24 22 29 29 H 1.09005 * 109.46610 * 179.97438 * 1 2 3 30 30 H 1.08992 * 109.47356 * 299.99464 * 1 2 3 31 31 H 1.09004 * 109.47026 * 59.99832 * 1 2 3 32 32 H 1.08996 * 109.46748 * 304.99824 * 4 2 1 33 33 H 1.09000 * 109.47463 * 185.00041 * 4 2 1 34 34 H 1.08999 * 109.47094 * 264.59815 * 6 5 4 35 35 H 1.08997 * 109.47370 * 24.59713 * 6 5 4 36 36 H 1.08997 * 109.47096 * 144.59755 * 6 5 4 37 37 H 0.96995 * 120.00094 * 5.44453 * 9 7 5 38 38 H 1.09008 * 109.47316 * 65.11735 * 10 9 7 39 39 H 1.08994 * 109.47504 * 305.11638 * 10 9 7 40 40 H 0.96995 * 119.99903 * 180.02562 * 15 14 12 41 41 H 1.08992 * 117.50043 * 215.00215 * 20 11 10 42 42 H 1.09004 * 117.50100 * 0.02599 * 20 11 10 43 43 H 1.09006 * 117.50037 * 359.97181 * 21 11 10 44 44 H 1.08999 * 117.49788 * 144.97000 * 21 11 10 45 45 H 1.09001 * 109.46781 * 60.00132 * 26 25 24 46 46 H 1.08993 * 109.47347 * 180.02562 * 26 25 24 47 47 H 1.08996 * 109.46998 * 300.00433 * 26 25 24 48 48 H 1.09006 * 109.47471 * 59.99976 * 27 25 24 49 49 H 1.09002 * 109.47622 * 179.97438 * 27 25 24 50 50 H 1.09001 * 109.47065 * 300.00099 * 27 25 24 51 51 H 1.09004 * 109.46952 * 60.00122 * 28 25 24 52 52 H 1.08996 * 109.47353 * 179.97438 * 28 25 24 53 53 H 1.08999 * 109.47146 * 300.00214 * 28 25 24 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.0277 0.0000 4 6 2.0401 -0.7211 -1.2493 5 7 1.6526 0.0355 -2.4424 6 6 0.3747 -0.2422 -3.1028 7 6 2.4716 0.9871 -2.9324 8 8 3.5620 1.1700 -2.4267 9 7 2.0818 1.7360 -3.9830 10 6 3.0155 2.6871 -4.5911 11 6 2.2913 3.4931 -5.6713 12 6 1.0702 4.2467 -5.2107 13 1 0.1013 3.7702 -5.2330 14 6 1.1908 5.5427 -4.7556 15 7 0.1138 6.2514 -4.4962 16 14 2.8840 6.2916 -4.5074 17 1 3.2839 7.0283 -5.7332 18 1 2.8467 7.2292 -3.3564 19 1 3.8675 5.2127 -4.2355 20 6 3.1458 4.1343 -6.7665 21 6 2.2737 2.9184 -7.0891 22 6 2.0323 -0.7104 1.2304 23 8 2.7271 -1.6927 1.1196 24 8 1.7080 -0.2519 2.4495 25 6 2.2267 -0.9864 3.5896 26 6 3.7561 -0.9860 3.5437 27 6 1.7547 -0.3200 4.8833 28 6 1.7155 -2.4277 3.5431 29 1 -0.3633 -1.0277 -0.0005 30 1 -0.3633 0.5137 0.8900 31 1 -0.3633 0.5139 -0.8900 32 1 1.6055 -1.7196 -1.2952 33 1 3.1264 -0.7993 -1.2056 34 1 0.5343 -0.9411 -3.9238 35 1 -0.3186 -0.6785 -2.3838 36 1 -0.0424 0.6865 -3.4922 37 1 1.1831 1.6453 -4.3365 38 1 3.8471 2.1435 -5.0398 39 1 3.3948 3.3639 -3.8256 40 1 0.1987 7.1632 -4.1764 41 1 2.8248 5.1012 -7.1539 42 1 4.2240 3.9816 -6.7178 43 1 2.7785 1.9662 -7.2525 44 1 1.3787 3.0850 -7.6886 45 1 4.1202 0.0408 3.5768 46 1 4.1454 -1.5377 4.3992 47 1 4.0920 -1.4610 2.6220 48 1 0.6651 -0.3207 4.9159 49 1 2.1437 -0.8716 5.7392 50 1 2.1185 0.7070 4.9167 51 1 2.0515 -2.9028 2.6214 52 1 2.1052 -2.9793 4.3986 53 1 0.6260 -2.4281 3.5758 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 ZINC000458746046.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:10:02 Heat of formation + Delta-G solvation = -90.690779 kcal Electronic energy + Delta-G solvation = -31126.056530 eV Core-core repulsion = 26999.161624 eV Total energy + Delta-G solvation = -4126.894906 eV No. of doubly occupied orbitals = 66 Molecular weight (most abundant/longest-lived isotopes) = 340.220 amu Computer time = 0.72 seconds Orbital eigenvalues (eV) -41.10692 -39.82462 -37.61326 -37.09217 -35.81517 -34.37989 -33.72181 -32.53455 -32.12458 -27.96954 -27.66744 -27.64602 -26.91578 -25.58696 -25.01764 -24.06212 -21.72135 -21.33784 -19.62122 -19.51327 -19.06003 -17.63757 -17.41108 -16.77709 -16.46673 -16.15974 -15.67284 -15.55441 -14.97874 -14.79732 -14.32690 -14.21245 -14.09759 -13.81702 -13.50262 -13.43308 -13.10380 -12.86916 -12.74767 -12.65658 -12.61258 -12.38305 -12.28720 -12.20433 -12.15863 -12.13809 -11.95119 -11.66907 -11.47207 -11.23577 -11.16916 -11.08330 -10.69754 -10.58491 -10.51682 -10.15018 -10.04360 -9.79784 -9.61409 -9.17044 -8.69234 -8.62461 -8.55515 -8.09023 -6.17771 -4.20939 2.13022 2.40730 3.07188 3.28962 3.39645 3.46195 3.63235 4.34335 4.40104 4.42290 4.48755 4.50606 4.63876 4.72471 4.77197 4.94433 5.11890 5.14284 5.15786 5.19887 5.26292 5.37233 5.41055 5.45165 5.50124 5.56975 5.60812 5.65771 5.68478 5.72805 5.73846 5.77063 5.93859 5.97613 6.03203 6.03963 6.06659 6.08114 6.26309 6.39503 6.43268 6.45680 6.65031 6.72245 7.00493 7.09163 7.18962 7.31018 7.68080 7.71522 8.38952 8.47346 8.81773 9.06504 9.46623 10.93434 Molecular weight = 340.22amu Principal moments of inertia in cm(-1) A = 0.015727 B = 0.002611 C = 0.002482 Principal moments of inertia in units of 10**(-40)*gram-cm**2 A = 1779.986211 B =10719.552026 C =11279.081458 Net atomic charges, atomic populations, and dipole contributions calculated with CM2 Atom NO. Type Charge No. of electrons 1 C -0.145 4.145 2 C -0.106 4.106 3 H 0.114 0.886 4 C 0.154 3.846 5 N -0.635 5.635 6 C 0.074 3.926 7 C 0.701 3.299 8 O -0.600 6.600 9 N -0.715 5.715 10 C 0.171 3.829 11 C -0.016 4.016 12 C -0.490 4.490 13 H 0.077 0.923 14 C 0.082 3.918 15 N -0.894 5.894 16 Si 0.876 3.124 17 H -0.279 1.279 18 H -0.293 1.293 19 H -0.255 1.255 20 C -0.192 4.192 21 C -0.179 4.179 22 C 0.478 3.522 23 O -0.505 6.505 24 O -0.348 6.348 25 C 0.115 3.885 26 C -0.184 4.184 27 C -0.141 4.141 28 C -0.182 4.182 29 H 0.058 0.942 30 H 0.068 0.932 31 H 0.078 0.922 32 H 0.082 0.918 33 H 0.094 0.906 34 H 0.054 0.946 35 H 0.075 0.925 36 H 0.073 0.927 37 H 0.401 0.599 38 H 0.039 0.961 39 H 0.065 0.935 40 H 0.265 0.735 41 H 0.090 0.910 42 H 0.074 0.926 43 H 0.072 0.928 44 H 0.087 0.913 45 H 0.062 0.938 46 H 0.069 0.931 47 H 0.084 0.916 48 H 0.066 0.934 49 H 0.077 0.923 50 H 0.067 0.933 51 H 0.084 0.916 52 H 0.069 0.931 53 H 0.061 0.939 For the charges calculated by CM2: Dipole moment (debyes) X Y Z Total from point charges -0.168 -16.370 12.816 20.791 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.202 4.202 2 C -0.127 4.127 3 H 0.132 0.868 4 C 0.030 3.970 5 N -0.373 5.373 6 C -0.069 4.069 7 C 0.402 3.598 8 O -0.479 6.479 9 N -0.370 5.370 10 C 0.046 3.954 11 C -0.017 4.017 12 C -0.527 4.527 13 H 0.095 0.905 14 C -0.121 4.121 15 N -0.531 5.531 16 Si 0.704 3.296 17 H -0.205 1.205 18 H -0.220 1.220 19 H -0.179 1.179 20 C -0.229 4.229 21 C -0.217 4.217 22 C 0.321 3.679 23 O -0.393 6.393 24 O -0.265 6.265 25 C 0.079 3.921 26 C -0.241 4.241 27 C -0.198 4.198 28 C -0.239 4.239 29 H 0.077 0.923 30 H 0.087 0.913 31 H 0.097 0.903 32 H 0.100 0.900 33 H 0.112 0.888 34 H 0.073 0.927 35 H 0.094 0.906 36 H 0.092 0.908 37 H 0.237 0.763 38 H 0.057 0.943 39 H 0.083 0.917 40 H 0.075 0.925 41 H 0.109 0.891 42 H 0.093 0.907 43 H 0.091 0.909 44 H 0.106 0.894 45 H 0.081 0.919 46 H 0.088 0.912 47 H 0.103 0.897 48 H 0.085 0.915 49 H 0.095 0.905 50 H 0.086 0.914 51 H 0.102 0.898 52 H 0.088 0.912 53 H 0.080 0.920 Dipole moment (debyes) X Y Z Total from point charges 0.775 -16.344 13.213 21.031 hybrid contribution 0.403 0.570 0.105 0.706 sum 1.179 -15.774 13.318 20.678 Atomic orbital electron populations 1.21868 0.91905 1.02472 1.03994 1.21601 0.98632 1.02372 0.90066 0.86782 1.20993 1.02202 0.92858 0.80963 1.46912 1.20784 1.39664 1.29953 1.21863 0.85230 1.01946 0.97869 1.15957 0.83717 0.79387 0.80715 1.90873 1.21964 1.72264 1.62775 1.45103 1.17641 1.41083 1.33128 1.21034 0.89485 0.89412 0.95434 1.19613 0.97399 0.98433 0.86268 1.20118 0.91664 0.97220 1.43701 0.90516 1.24943 0.96501 0.96678 0.94013 1.69896 1.28523 1.09874 1.44758 0.86903 0.85617 0.79514 0.77590 1.20481 1.22016 1.17911 1.22917 0.95918 0.97628 1.06483 1.22869 1.01476 0.95572 1.01733 1.23514 0.77050 0.82630 0.84726 1.90564 1.35464 1.25996 1.87280 1.86570 1.69603 1.55626 1.14744 1.22181 0.94137 0.91662 0.84118 1.22495 0.93945 1.03113 1.04559 1.21850 1.02363 1.01536 0.94044 1.22482 1.01566 0.95295 1.04546 0.92257 0.91290 0.90299 0.89985 0.88766 0.92704 0.90619 0.90841 0.76261 0.94265 0.91656 0.92507 0.89135 0.90739 0.90924 0.89428 0.91888 0.91162 0.89680 0.91480 0.90452 0.91388 0.89757 0.91170 0.92042 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.15 -1.45 8.11 37.16 0.30 -1.15 16 2 C -0.11 -1.31 3.03 -91.77 -0.28 -1.59 16 3 H 0.11 1.64 8.14 -51.93 -0.42 1.22 16 4 C 0.15 2.01 4.25 -4.04 -0.02 2.00 16 5 N -0.64 -9.00 2.70 -187.36 -0.51 -9.51 16 6 C 0.07 0.79 9.10 59.85 0.54 1.33 16 7 C 0.70 13.10 8.27 -86.92 -0.72 12.38 16 8 O -0.60 -13.37 15.88 5.29 0.08 -13.29 16 9 N -0.72 -13.52 5.25 -65.23 -0.34 -13.86 16 10 C 0.17 3.61 5.34 -4.04 -0.02 3.58 16 11 C -0.02 -0.34 0.93 -155.66 -0.15 -0.49 16 12 C -0.49 -12.62 8.92 -34.44 -0.31 -12.93 16 13 H 0.08 2.11 8.06 -52.49 -0.42 1.69 16 14 C 0.08 2.24 5.46 -17.82 -0.10 2.15 16 15 N -0.89 -26.09 13.96 55.43 0.77 -25.31 16 16 Si 0.88 21.38 25.04 -169.99 -4.26 17.12 16 17 H -0.28 -6.65 7.11 56.52 0.40 -6.25 16 18 H -0.29 -7.30 7.11 56.52 0.40 -6.90 16 19 H -0.25 -6.04 5.52 56.52 0.31 -5.72 16 20 C -0.19 -3.91 7.70 -26.69 -0.21 -4.12 16 21 C -0.18 -3.40 9.94 -26.69 -0.27 -3.67 16 22 C 0.48 6.07 6.57 36.00 0.24 6.31 16 23 O -0.50 -7.02 11.42 -18.75 -0.21 -7.24 16 24 O -0.35 -4.17 9.79 -39.24 -0.38 -4.55 16 25 C 0.11 1.13 1.13 -90.62 -0.10 1.03 16 26 C -0.18 -1.83 8.37 37.16 0.31 -1.52 16 27 C -0.14 -1.06 8.85 37.15 0.33 -0.73 16 28 C -0.18 -1.62 8.37 37.16 0.31 -1.31 16 29 H 0.06 0.50 8.14 -51.93 -0.42 0.08 16 30 H 0.07 0.67 8.11 -51.93 -0.42 0.25 16 31 H 0.08 0.81 5.99 -51.93 -0.31 0.49 16 32 H 0.08 0.89 8.07 -51.93 -0.42 0.47 16 33 H 0.09 1.47 7.18 -51.93 -0.37 1.09 16 34 H 0.05 0.52 8.14 -51.93 -0.42 0.10 16 35 H 0.08 0.62 6.23 -51.93 -0.32 0.30 16 36 H 0.07 0.83 6.43 -51.93 -0.33 0.50 16 37 H 0.40 6.99 5.97 -40.82 -0.24 6.75 16 38 H 0.04 0.79 8.11 -51.92 -0.42 0.37 16 39 H 0.07 1.55 6.29 -51.93 -0.33 1.22 16 40 H 0.27 7.35 8.31 -40.82 -0.34 7.02 16 41 H 0.09 1.95 7.29 -51.93 -0.38 1.57 16 42 H 0.07 1.39 8.14 -51.93 -0.42 0.97 16 43 H 0.07 1.23 8.14 -51.93 -0.42 0.81 16 44 H 0.09 1.67 8.14 -51.93 -0.42 1.24 16 45 H 0.06 0.63 8.14 -51.93 -0.42 0.21 16 46 H 0.07 0.55 8.14 -51.93 -0.42 0.13 16 47 H 0.08 1.02 5.89 -51.93 -0.31 0.72 16 48 H 0.07 0.48 8.14 -51.93 -0.42 0.05 16 49 H 0.08 0.44 8.14 -51.93 -0.42 0.01 16 50 H 0.07 0.52 8.14 -51.93 -0.42 0.10 16 51 H 0.08 0.91 5.89 -51.93 -0.31 0.60 16 52 H 0.07 0.50 8.14 -51.93 -0.42 0.07 16 53 H 0.06 0.51 8.14 -51.93 -0.42 0.08 16 LS Contribution 411.78 15.07 6.21 6.21 Total: -1.00 -31.82 411.78 -8.07 -39.89 By element: Atomic # 1 Polarization: 18.56 SS G_CDS: -9.31 Total: 9.25 kcal Atomic # 6 Polarization: 1.42 SS G_CDS: -0.12 Total: 1.29 kcal Atomic # 7 Polarization: -48.61 SS G_CDS: -0.07 Total: -48.68 kcal Atomic # 8 Polarization: -24.56 SS G_CDS: -0.51 Total: -25.08 kcal Atomic # 14 Polarization: 21.38 SS G_CDS: -4.26 Total: 17.12 kcal Total LS contribution 6.21 Total: 6.21 kcal Total: -31.82 -8.07 -39.89 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. ZINC000458746046.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 -50.797 kcal (2) G-P(sol) polarization free energy of solvation -31.823 kcal (3) G-ENP(sol) elect.-nuc.-pol. free energy of system -82.620 kcal (4) G-CDS(sol) cavity-dispersion-solvent structure free energy -8.071 kcal (5) G-P-CDS(sol) = G-P(sol) + G-CDS(sol) = (2) + (4) -39.894 kcal (6) G-S(sol) free energy of system = (1) + (5) -90.691 kcal 1SCF WAS SPECIFIED, SO GEOMETRY OPTIMIZATION WAS NOT USED Total computer time = 0.72 seconds