Wall clock time and date at job start Sat Apr 11 2020 06:13:14 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.42895 * 1 3 3 C 1.42898 * 114.00062 * 2 1 4 4 H 1.09005 * 109.47090 * 30.00443 * 3 2 1 5 5 C 1.50697 * 109.47428 * 270.00514 * 3 2 1 6 6 O 1.21277 * 119.99976 * 354.99255 * 5 3 2 7 7 N 1.34779 * 119.99906 * 175.00462 * 5 3 2 8 8 C 1.46505 * 119.99700 * 179.97438 * 7 5 3 9 9 H 1.08993 * 109.47130 * 324.99787 * 8 7 5 10 10 C 1.53006 * 109.47023 * 84.99299 * 8 7 5 11 11 C 1.53000 * 109.46846 * 179.97438 * 10 8 7 12 12 C 1.52993 * 109.47345 * 300.00008 * 11 10 8 13 13 H 1.09002 * 109.47493 * 300.00298 * 12 11 10 14 14 C 1.53000 * 109.47345 * 180.02562 * 12 11 10 15 15 N 1.46499 * 109.46840 * 65.00000 * 14 12 11 16 16 C 1.36936 * 120.00002 * 180.02562 * 15 14 12 17 17 H 1.08004 * 120.00140 * 0.02562 * 16 15 14 18 18 C 1.38811 * 120.00367 * 180.02562 * 16 15 14 19 19 N 1.31622 * 120.00562 * 180.02562 * 18 16 15 20 20 Si 1.86804 * 119.99699 * 359.97438 * 18 16 15 21 21 H 1.48497 * 109.47093 * 90.00457 * 20 18 16 22 22 H 1.48504 * 109.47189 * 210.00360 * 20 18 16 23 23 H 1.48500 * 109.47392 * 329.99935 * 20 18 16 24 24 C 1.53006 * 109.47155 * 60.00181 * 12 11 10 25 25 C 1.52993 * 109.47015 * 204.99213 * 8 7 5 26 26 C 1.53003 * 109.47303 * 149.99910 * 3 2 1 27 27 H 1.09002 * 110.80076 * 60.70035 * 26 3 2 28 28 C 1.54805 * 110.90883 * 297.07903 * 26 3 2 29 29 C 1.54100 * 104.18912 * 220.91868 * 28 26 3 30 30 O 1.43742 * 107.29228 * 357.96649 * 29 28 26 31 31 C 1.44417 * 107.00023 * 26.49221 * 30 29 28 32 32 H 1.08999 * 109.47042 * 180.02562 * 1 2 3 33 33 H 1.09001 * 109.46854 * 299.99516 * 1 2 3 34 34 H 1.08995 * 109.47485 * 59.99904 * 1 2 3 35 35 H 0.96996 * 119.99948 * 359.97438 * 7 5 3 36 36 H 1.09001 * 109.47241 * 300.00223 * 10 8 7 37 37 H 1.09002 * 109.46824 * 60.00121 * 10 8 7 38 38 H 1.09002 * 109.46642 * 180.02562 * 11 10 8 39 39 H 1.08993 * 109.46977 * 60.00690 * 11 10 8 40 40 H 1.09002 * 109.47398 * 185.00327 * 14 12 11 41 41 H 1.09001 * 109.47194 * 305.00507 * 14 12 11 42 42 H 0.97000 * 120.00420 * 0.02562 * 15 14 12 43 43 H 0.96993 * 120.00222 * 179.97438 * 19 18 16 44 44 H 1.09002 * 109.46824 * 59.99708 * 24 12 11 45 45 H 1.09001 * 109.47241 * 179.97438 * 24 12 11 46 46 H 1.09002 * 109.47179 * 60.00402 * 25 8 7 47 47 H 1.09002 * 109.47493 * 300.00284 * 25 8 7 48 48 H 1.09002 * 110.49037 * 339.59461 * 28 26 3 49 49 H 1.09006 * 110.49055 * 102.10775 * 28 26 3 50 50 H 1.09003 * 109.88092 * 117.39499 * 29 28 26 51 51 H 1.08998 * 109.88185 * 238.54702 * 29 28 26 52 52 H 1.09002 * 110.65579 * 201.48302 * 31 30 29 53 53 H 1.09001 * 110.60699 * 78.46914 * 31 30 29 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.3054 0.0000 4 1 1.3449 1.9994 0.5139 5 6 2.2144 1.7644 -1.4208 6 8 1.9791 1.0107 -2.3413 7 7 2.6588 3.0125 -1.6684 8 6 2.8579 3.4585 -3.0496 9 1 3.2089 2.6225 -3.6544 10 6 1.5333 3.9757 -3.6145 11 6 1.7414 4.4421 -5.0568 12 6 2.7807 5.5644 -5.0849 13 1 2.4294 6.4003 -4.4799 14 6 2.9883 6.0312 -6.5271 15 7 1.7513 6.6359 -7.0275 16 6 1.6976 7.1287 -8.3040 17 1 2.5672 7.0713 -8.9420 18 6 0.5257 7.7021 -8.7780 19 7 0.4740 8.1754 -10.0051 20 14 -0.9785 7.8008 -7.6747 21 1 -0.9535 9.0775 -6.9168 22 1 -2.2072 7.7429 -8.5066 23 1 -0.9705 6.6603 -6.7237 24 6 4.1053 5.0473 -4.5201 25 6 3.8974 4.5807 -3.0784 26 6 3.3589 1.2633 0.7213 27 1 4.0463 0.5838 0.2173 28 6 3.1789 0.8656 2.2065 29 6 4.1512 1.7894 2.9654 30 8 4.8394 2.5835 1.9846 31 6 3.9694 2.6833 0.8361 32 1 -0.3633 -1.0277 0.0005 33 1 -0.3633 0.5138 0.8901 34 1 -0.3634 0.5138 -0.8899 35 1 2.8468 3.6154 -0.9322 36 1 0.7927 3.1762 -3.5940 37 1 1.1822 4.8118 -3.0097 38 1 0.7978 4.8109 -5.4590 39 1 2.0926 3.6063 -5.6617 40 1 3.7914 6.7676 -6.5587 41 1 3.2535 5.1775 -7.1508 42 1 0.9702 6.6871 -6.4547 43 1 -0.3450 8.5757 -10.3363 44 1 4.4566 4.2114 -5.1251 45 1 4.8459 5.8468 -4.5406 46 1 4.8410 4.2123 -2.6760 47 1 3.5463 5.4169 -2.4736 48 1 2.1529 1.0451 2.5279 49 1 3.4473 -0.1801 2.3568 50 1 3.5947 2.4372 3.6427 51 1 4.8685 1.1912 3.5273 52 1 4.5409 2.9288 -0.0590 53 1 3.1887 3.4239 1.0095 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 ZINC001724749821.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:13:14 Heat of formation + Delta-G solvation = -129.492739 kcal Electronic energy + Delta-G solvation = -30888.334003 eV Core-core repulsion = 26759.756518 eV Total energy + Delta-G solvation = -4128.577485 eV No. of doubly occupied orbitals = 66 Molecular weight (most abundant/longest-lived isotopes) = 340.220 amu Computer time = 1.13 seconds Orbital eigenvalues (eV) -41.04063 -39.61646 -38.70699 -37.53883 -35.39857 -34.08704 -33.25325 -31.54155 -30.81339 -30.43145 -29.27747 -27.84279 -26.31166 -25.03243 -23.39067 -22.84874 -22.27097 -21.45718 -20.83010 -20.72258 -19.74749 -18.12923 -17.59837 -17.38360 -17.01942 -16.58888 -15.99911 -15.67594 -15.29758 -15.02386 -14.84146 -14.59885 -14.36974 -14.29462 -14.07328 -13.99237 -13.72323 -13.41668 -12.99026 -12.80033 -12.52214 -12.36228 -12.28320 -12.18482 -12.03126 -11.96372 -11.87765 -11.75059 -11.31683 -11.22884 -11.12496 -11.08210 -10.93082 -10.82386 -10.37875 -10.27591 -10.21532 -10.15317 -9.99826 -9.59415 -9.52260 -9.22522 -8.86641 -7.19093 -5.74447 -3.09140 2.22991 2.64306 3.13640 3.37170 3.37780 3.40120 3.45692 3.92201 4.13654 4.26744 4.37090 4.45700 4.52930 4.57389 4.60761 4.64563 4.70249 4.75542 4.89136 4.94102 5.04966 5.07954 5.17023 5.17651 5.18621 5.24155 5.31878 5.39942 5.41896 5.43795 5.51405 5.53009 5.56434 5.57560 5.80002 5.86517 5.93734 6.02310 6.12541 6.17978 6.24335 6.26248 6.31981 6.44935 6.50446 6.56787 6.74939 6.86668 7.16841 7.90551 8.02818 8.44999 9.24591 9.97061 10.17383 11.44233 Molecular weight = 340.22amu Principal moments of inertia in cm(-1) A = 0.018372 B = 0.002252 C = 0.002115 Principal moments of inertia in units of 10**(-40)*gram-cm**2 A = 1523.731351 B =12428.753595 C =13236.644529 Net atomic charges, atomic populations, and dipole contributions calculated with CM2 Atom NO. Type Charge No. of electrons 1 C 0.028 3.972 2 O -0.353 6.353 3 C 0.084 3.916 4 H 0.088 0.912 5 C 0.488 3.512 6 O -0.516 6.516 7 N -0.728 5.728 8 C 0.143 3.857 9 H 0.090 0.910 10 C -0.127 4.127 11 C -0.109 4.109 12 C -0.099 4.099 13 H 0.069 0.931 14 C 0.161 3.839 15 N -0.756 5.756 16 C -0.241 4.241 17 H 0.090 0.910 18 C 0.002 3.998 19 N -0.935 5.935 20 Si 0.884 3.116 21 H -0.285 1.285 22 H -0.291 1.291 23 H -0.277 1.277 24 C -0.107 4.107 25 C -0.124 4.124 26 C -0.116 4.116 27 H 0.102 0.898 28 C -0.143 4.143 29 C 0.043 3.957 30 O -0.369 6.369 31 C 0.046 3.954 32 H 0.092 0.908 33 H 0.042 0.958 34 H 0.047 0.953 35 H 0.396 0.604 36 H 0.062 0.938 37 H 0.061 0.939 38 H 0.070 0.930 39 H 0.065 0.935 40 H 0.026 0.974 41 H 0.026 0.974 42 H 0.360 0.640 43 H 0.259 0.741 44 H 0.066 0.934 45 H 0.070 0.930 46 H 0.067 0.933 47 H 0.067 0.933 48 H 0.078 0.922 49 H 0.085 0.915 50 H 0.061 0.939 51 H 0.087 0.913 52 H 0.103 0.897 53 H 0.065 0.935 For the charges calculated by CM2: Dipole moment (debyes) X Y Z Total from point charges 5.436 -14.694 30.574 34.355 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.068 4.068 2 O -0.270 6.270 3 C 0.022 3.978 4 H 0.105 0.895 5 C 0.272 3.728 6 O -0.390 6.390 7 N -0.383 5.383 8 C 0.040 3.960 9 H 0.108 0.892 10 C -0.166 4.166 11 C -0.147 4.147 12 C -0.118 4.118 13 H 0.088 0.912 14 C 0.034 3.966 15 N -0.401 5.401 16 C -0.371 4.371 17 H 0.107 0.893 18 C -0.194 4.194 19 N -0.584 5.584 20 Si 0.715 3.285 21 H -0.212 1.212 22 H -0.218 1.218 23 H -0.204 1.204 24 C -0.145 4.145 25 C -0.162 4.162 26 C -0.135 4.135 27 H 0.121 0.879 28 C -0.180 4.180 29 C -0.034 4.034 30 O -0.288 6.288 31 C -0.031 4.031 32 H 0.111 0.889 33 H 0.061 0.939 34 H 0.066 0.934 35 H 0.232 0.768 36 H 0.081 0.919 37 H 0.080 0.920 38 H 0.089 0.911 39 H 0.084 0.916 40 H 0.044 0.956 41 H 0.044 0.956 42 H 0.193 0.807 43 H 0.069 0.931 44 H 0.085 0.915 45 H 0.088 0.912 46 H 0.085 0.915 47 H 0.086 0.914 48 H 0.096 0.904 49 H 0.104 0.896 50 H 0.080 0.920 51 H 0.105 0.895 52 H 0.121 0.879 53 H 0.083 0.917 Dipole moment (debyes) X Y Z Total from point charges 5.643 -15.296 29.700 33.880 hybrid contribution -1.675 0.557 0.939 1.999 sum 3.968 -14.739 30.639 34.230 Atomic orbital electron populations 1.22907 0.80190 1.03552 1.00183 1.87991 1.20340 1.24305 1.94382 1.21978 0.95861 0.87438 0.92495 0.89461 1.21282 0.77544 0.81994 0.91968 1.90748 1.48785 1.53669 1.45803 1.46428 1.68155 1.14191 1.09487 1.21123 0.95909 0.97423 0.81580 0.89185 1.21929 0.96717 1.00553 0.97368 1.21515 1.00115 0.98318 0.94765 1.21410 0.94774 0.97660 0.97924 0.91243 1.20662 0.89622 0.95537 0.90771 1.42734 1.10256 1.73001 1.14146 1.19081 0.98937 1.30577 0.88522 0.89261 1.24677 0.97566 1.00883 0.96228 1.69700 1.27483 1.45663 1.15531 0.86549 0.83952 0.77756 0.80280 1.21152 1.21771 1.20360 1.21447 0.97723 1.01135 0.94213 1.21796 0.98157 0.97492 0.98730 1.22764 0.95268 0.98983 0.96532 0.87942 1.22912 1.00583 1.01442 0.93102 1.23247 0.95338 0.92999 0.91769 1.88943 1.40929 1.72483 1.26440 1.23746 0.94143 0.94727 0.90440 0.88923 0.93893 0.93395 0.76847 0.91893 0.92043 0.91087 0.91644 0.95573 0.95559 0.80699 0.93117 0.91485 0.91155 0.91472 0.91438 0.90356 0.89593 0.92040 0.89472 0.87867 0.91715 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 C 0.03 0.21 10.81 101.05 1.09 1.30 16 2 O -0.35 -3.39 9.66 -55.09 -0.53 -3.93 16 3 C 0.08 0.65 2.00 -27.88 -0.06 0.60 16 4 H 0.09 0.47 7.78 -51.93 -0.40 0.07 16 5 C 0.49 5.17 7.06 -10.99 -0.08 5.09 16 6 O -0.52 -7.88 16.17 -14.30 -0.23 -8.11 16 7 N -0.73 -6.21 4.33 -53.81 -0.23 -6.44 16 8 C 0.14 1.51 2.45 -67.93 -0.17 1.34 16 9 H 0.09 1.14 7.58 -51.93 -0.39 0.74 16 10 C -0.13 -1.56 5.50 -26.73 -0.15 -1.70 16 11 C -0.11 -1.60 5.08 -26.73 -0.14 -1.74 16 12 C -0.10 -1.44 2.22 -90.62 -0.20 -1.64 16 13 H 0.07 1.02 8.14 -51.93 -0.42 0.59 16 14 C 0.16 3.09 5.65 -4.04 -0.02 3.06 16 15 N -0.76 -18.08 5.13 -59.91 -0.31 -18.38 16 16 C -0.24 -6.93 10.47 -15.06 -0.16 -7.09 16 17 H 0.09 2.75 8.06 -52.48 -0.42 2.32 16 18 C 0.00 0.06 5.50 -17.43 -0.10 -0.04 16 19 N -0.94 -29.78 13.97 55.49 0.78 -29.00 16 20 Si 0.88 22.35 27.73 -169.99 -4.71 17.64 16 21 H -0.28 -7.12 7.11 56.52 0.40 -6.72 16 22 H -0.29 -7.47 7.11 56.52 0.40 -7.06 16 23 H -0.28 -6.57 6.52 56.52 0.37 -6.21 16 24 C -0.11 -1.20 5.20 -26.75 -0.14 -1.34 16 25 C -0.12 -1.12 5.39 -26.76 -0.14 -1.26 16 26 C -0.12 -0.76 2.88 -88.61 -0.25 -1.02 16 27 H 0.10 0.81 8.14 -51.93 -0.42 0.39 16 28 C -0.14 -0.72 6.38 -25.21 -0.16 -0.88 16 29 C 0.04 0.23 7.70 37.73 0.29 0.52 16 30 O -0.37 -2.88 11.52 -35.23 -0.41 -3.28 16 31 C 0.05 0.27 5.30 38.21 0.20 0.47 16 32 H 0.09 0.57 8.14 -51.93 -0.42 0.15 16 33 H 0.04 0.22 8.09 -51.93 -0.42 -0.20 16 34 H 0.05 0.41 8.14 -51.93 -0.42 -0.01 16 35 H 0.40 2.29 5.83 -40.82 -0.24 2.05 16 36 H 0.06 0.85 8.14 -51.93 -0.42 0.43 16 37 H 0.06 0.70 8.14 -51.93 -0.42 0.28 16 38 H 0.07 1.21 7.18 -51.93 -0.37 0.84 16 39 H 0.06 0.99 8.14 -51.93 -0.42 0.57 16 40 H 0.03 0.51 8.14 -51.93 -0.42 0.09 16 41 H 0.03 0.52 8.14 -51.93 -0.42 0.10 16 42 H 0.36 8.38 4.71 -40.82 -0.19 8.19 16 43 H 0.26 7.83 8.31 -40.82 -0.34 7.49 16 44 H 0.07 0.77 8.14 -51.93 -0.42 0.34 16 45 H 0.07 0.75 8.14 -51.93 -0.42 0.33 16 46 H 0.07 0.52 8.14 -51.93 -0.42 0.10 16 47 H 0.07 0.55 8.14 -51.93 -0.42 0.13 16 48 H 0.08 0.33 7.86 -51.93 -0.41 -0.08 16 49 H 0.09 0.43 8.14 -51.93 -0.42 0.00 16 50 H 0.06 0.26 8.14 -51.93 -0.42 -0.17 16 51 H 0.09 0.38 8.14 -51.93 -0.42 -0.04 16 52 H 0.10 0.58 7.31 -51.93 -0.38 0.20 16 53 H 0.06 0.26 7.37 -51.93 -0.38 -0.12 16 LS Contribution 409.19 15.07 6.17 6.17 Total: -1.00 -35.68 409.19 -9.20 -44.89 By element: Atomic # 1 Polarization: 14.34 SS G_CDS: -9.55 Total: 4.79 kcal Atomic # 6 Polarization: -4.16 SS G_CDS: -0.17 Total: -4.33 kcal Atomic # 7 Polarization: -54.06 SS G_CDS: 0.23 Total: -53.83 kcal Atomic # 8 Polarization: -14.15 SS G_CDS: -1.17 Total: -15.32 kcal Atomic # 14 Polarization: 22.35 SS G_CDS: -4.71 Total: 17.64 kcal Total LS contribution 6.17 Total: 6.17 kcal Total: -35.68 -9.20 -44.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. ZINC001724749821.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 -84.607 kcal (2) G-P(sol) polarization free energy of solvation -35.682 kcal (3) G-ENP(sol) elect.-nuc.-pol. free energy of system -120.289 kcal (4) G-CDS(sol) cavity-dispersion-solvent structure free energy -9.204 kcal (5) G-P-CDS(sol) = G-P(sol) + G-CDS(sol) = (2) + (4) -44.886 kcal (6) G-S(sol) free energy of system = (1) + (5) -129.493 kcal 1SCF WAS SPECIFIED, SO GEOMETRY OPTIMIZATION WAS NOT USED Total computer time = 1.13 seconds