Wall clock time and date at job start Tue Mar 31 2020 20:43:08 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.52996 * 1 3 3 H 1.09002 * 109.46863 * 2 1 4 4 O 1.42899 * 109.47526 * 239.99865 * 2 1 3 5 5 C 1.42896 * 114.00254 * 89.99790 * 4 2 1 6 6 C 1.53002 * 109.47074 * 180.02562 * 5 4 2 7 7 C 1.53002 * 117.49786 * 243.74213 * 6 5 4 8 8 C 1.53003 * 117.49665 * 175.08709 * 6 5 4 9 9 C 1.50702 * 109.47447 * 119.99242 * 2 1 3 10 10 O 1.21285 * 119.99932 * 114.99843 * 9 2 1 11 11 N 1.34776 * 119.99885 * 295.01166 * 9 2 1 12 12 C 1.46500 * 119.99720 * 179.97438 * 11 9 2 13 13 C 1.53004 * 115.54953 * 179.97438 * 12 11 9 14 14 C 1.53039 * 109.45791 * 34.41719 * 13 12 11 15 15 C 1.53194 * 109.31279 * 178.60906 * 14 13 12 16 16 N 1.46925 * 108.77112 * 54.64111 * 15 14 13 17 17 C 1.36927 * 120.63534 * 126.41528 * 16 15 14 18 18 H 1.08001 * 120.00618 * 207.94201 * 17 16 15 19 19 C 1.38820 * 120.00081 * 27.94124 * 17 16 15 20 20 N 1.31615 * 119.99968 * 17.51930 * 19 17 16 21 21 Si 1.86799 * 119.99918 * 197.52506 * 19 17 16 22 22 H 1.48501 * 109.47488 * 94.06659 * 21 19 17 23 23 H 1.48502 * 109.46661 * 214.06718 * 21 19 17 24 24 H 1.48495 * 109.47274 * 334.06079 * 21 19 17 25 25 C 1.46931 * 118.73572 * 306.38465 * 16 15 14 26 26 C 1.53039 * 109.46178 * 274.34890 * 13 12 11 27 27 C 1.53001 * 117.50119 * 325.66470 * 12 11 9 28 28 C 1.53002 * 117.50181 * 34.32552 * 12 11 9 29 29 H 1.09001 * 109.47181 * 60.00659 * 1 2 3 30 30 H 1.09002 * 109.46863 * 179.97438 * 1 2 3 31 31 H 1.09002 * 109.47457 * 300.00461 * 1 2 3 32 32 H 1.09005 * 109.46971 * 300.00325 * 5 4 2 33 33 H 1.09000 * 109.46897 * 60.00199 * 5 4 2 34 34 H 1.08997 * 115.54870 * 29.42006 * 6 5 4 35 35 H 1.08995 * 117.49553 * 0.02562 * 7 6 5 36 36 H 1.09003 * 117.49577 * 145.01845 * 7 6 5 37 37 H 1.08997 * 117.49604 * 214.98029 * 8 6 5 38 38 H 1.09003 * 117.49842 * 359.97438 * 8 6 5 39 39 H 0.97001 * 120.00227 * 359.97438 * 11 9 2 40 40 H 1.09006 * 109.44674 * 154.37866 * 13 12 11 41 41 H 1.08998 * 109.50115 * 298.56564 * 14 13 12 42 42 H 1.09005 * 109.49582 * 58.66190 * 14 13 12 43 43 H 1.08999 * 109.58540 * 174.41777 * 15 14 13 44 44 H 1.08998 * 109.59002 * 294.85034 * 15 14 13 45 45 H 0.96996 * 120.00010 * 180.02562 * 20 19 17 46 46 H 1.09001 * 109.58574 * 293.82691 * 25 16 15 47 47 H 1.08997 * 109.58825 * 173.39908 * 25 16 15 48 48 H 1.09002 * 109.49328 * 301.33591 * 26 13 12 49 49 H 1.08997 * 109.49711 * 61.40210 * 26 13 12 50 50 H 1.08998 * 117.49823 * 215.00697 * 27 12 11 51 51 H 1.09003 * 117.49358 * 0.02562 * 27 12 11 52 52 H 1.09002 * 117.49655 * 359.97438 * 28 12 11 53 53 H 1.09002 * 117.50210 * 144.97643 * 28 12 11 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 8 2.0064 -0.6736 -1.1667 5 6 2.2001 0.1828 -2.2940 6 6 2.7080 -0.6410 -3.4791 7 6 4.1061 -0.3228 -4.0131 8 6 2.8687 0.0763 -4.8210 9 6 2.0324 -0.7102 1.2305 10 8 2.6530 -1.7465 1.1217 11 7 1.7918 -0.1940 2.4520 12 6 2.2797 -0.8847 3.6483 13 6 1.9500 -0.2092 4.9809 14 6 0.5810 0.4679 4.8844 15 6 0.2389 1.1134 6.2309 16 7 1.3329 2.0162 6.6139 17 6 1.0796 3.3239 6.9316 18 1 1.8394 4.0750 6.7736 19 6 -0.1539 3.6842 7.4567 20 7 -0.9448 2.7654 7.9690 21 14 -0.6912 5.4732 7.4360 22 1 -1.5078 5.7308 6.2227 23 1 -1.5003 5.7609 8.6475 24 1 0.5073 6.3498 7.4215 25 6 2.7039 1.4891 6.6503 26 6 3.0184 0.8394 5.2991 27 6 2.3073 -2.4142 3.6191 28 6 3.6156 -1.6231 3.5423 29 1 -0.3633 0.5137 -0.8900 30 1 -0.3633 -1.0277 -0.0005 31 1 -0.3634 0.5139 0.8899 32 1 1.2538 0.6545 -2.5588 33 1 2.9320 0.9512 -2.0450 34 1 2.3860 -1.6821 -3.5003 35 1 4.6641 0.4730 -3.5200 36 1 4.7037 -1.1547 -4.3860 37 1 2.6527 -0.4930 -5.7250 38 1 2.6122 1.1350 -4.8596 39 1 1.2955 0.6350 2.5391 40 1 1.9295 -0.9585 5.7724 41 1 0.6080 1.2344 4.1099 42 1 -0.1759 -0.2758 4.6347 43 1 -0.6876 1.6803 6.1395 44 1 0.1223 0.3385 6.9885 45 1 -1.8069 3.0171 8.3355 46 1 2.7883 0.7448 7.4422 47 1 3.4039 2.3038 6.8356 48 1 3.9963 0.3601 5.3446 49 1 3.0228 1.6020 4.5204 50 1 2.0903 -2.9393 4.5492 51 1 1.9628 -2.9064 2.7096 52 1 4.1312 -1.5953 2.5823 53 1 4.2594 -1.6276 4.4218 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) 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 ZINC001284844282.mol2 53 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 20:43:08 Heat of formation + Delta-G solvation = -10.255935 kcal Electronic energy + Delta-G solvation = -30168.476754 eV Core-core repulsion = 26235.525817 eV Total energy + Delta-G solvation = -3932.950938 eV No. of doubly occupied orbitals = 65 Molecular weight (most abundant/longest-lived isotopes) = 336.220 amu Computer time = 1.04 seconds Orbital eigenvalues (eV) -40.59400 -38.92889 -37.89792 -36.81958 -35.79606 -33.40271 -32.76908 -31.11271 -30.43024 -29.65855 -27.41308 -25.64495 -24.53025 -22.82275 -22.74367 -22.53063 -21.67869 -21.18082 -20.63859 -19.79293 -18.86299 -17.66488 -17.26908 -16.76388 -16.09109 -15.71939 -15.29475 -14.89015 -14.78527 -14.61851 -14.17589 -13.94523 -13.67358 -13.53514 -13.38646 -13.24901 -13.03602 -12.90215 -12.68609 -12.56925 -12.30564 -11.93949 -11.82735 -11.51235 -11.32387 -11.20567 -11.10667 -10.96478 -10.82301 -10.76940 -10.71599 -10.62886 -10.39264 -10.26417 -10.11875 -9.95754 -9.90525 -9.59416 -9.45305 -9.24463 -8.52060 -8.40315 -6.69961 -5.81892 -3.23901 2.54951 2.67480 3.34223 3.36282 3.41973 3.43232 3.67965 4.12427 4.40348 4.46794 4.49114 4.52640 4.64904 4.69064 4.81626 4.86734 4.95399 5.08792 5.10980 5.23241 5.25455 5.28437 5.33580 5.35667 5.41121 5.46778 5.52741 5.61329 5.69307 5.73383 5.76562 5.84535 5.86257 5.92298 5.99053 6.08972 6.20250 6.29719 6.49750 6.50757 6.54780 6.56714 6.71679 6.73250 6.80670 6.85701 6.98702 7.31236 7.56303 7.61546 7.77068 8.24672 8.42209 9.18715 9.82275 10.39330 11.26562 Molecular weight = 336.22amu Principal moments of inertia in cm(-1) A = 0.013756 B = 0.002586 C = 0.002337 Principal moments of inertia in units of 10**(-40)*gram-cm**2 A = 2034.950365 B =10824.572775 C =11980.113505 Net atomic charges, atomic populations, and dipole contributions calculated with CM2 Atom NO. Type Charge No. of electrons 1 C -0.171 4.171 2 C 0.071 3.929 3 H 0.085 0.915 4 O -0.353 6.353 5 C 0.090 3.910 6 C -0.151 4.151 7 C -0.177 4.177 8 C -0.174 4.174 9 C 0.513 3.487 10 O -0.518 6.518 11 N -0.687 5.687 12 C 0.110 3.890 13 C -0.042 4.042 14 C -0.133 4.133 15 C 0.171 3.829 16 N -0.592 5.592 17 C -0.252 4.252 18 H 0.065 0.935 19 C 0.002 3.998 20 N -0.902 5.902 21 Si 0.901 3.099 22 H -0.288 1.288 23 H -0.292 1.292 24 H -0.279 1.279 25 C 0.144 3.856 26 C -0.124 4.124 27 C -0.179 4.179 28 C -0.177 4.177 29 H 0.069 0.931 30 H 0.071 0.929 31 H 0.075 0.925 32 H 0.055 0.945 33 H 0.051 0.949 34 H 0.115 0.885 35 H 0.095 0.905 36 H 0.096 0.904 37 H 0.098 0.902 38 H 0.095 0.905 39 H 0.407 0.593 40 H 0.079 0.921 41 H 0.051 0.949 42 H 0.052 0.948 43 H 0.087 0.913 44 H 0.019 0.981 45 H 0.253 0.747 46 H 0.017 0.983 47 H 0.053 0.947 48 H 0.052 0.948 49 H 0.057 0.943 50 H 0.091 0.909 51 H 0.104 0.896 52 H 0.104 0.896 53 H 0.092 0.908 For the charges calculated by CM2: Dipole moment (debyes) X Y Z Total from point charges 5.214 -5.978 -19.771 21.303 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.229 4.229 2 C 0.010 3.990 3 H 0.103 0.897 4 O -0.271 6.271 5 C 0.014 3.986 6 C -0.170 4.170 7 C -0.214 4.214 8 C -0.211 4.211 9 C 0.298 3.702 10 O -0.393 6.393 11 N -0.339 5.339 12 C 0.023 3.977 13 C -0.061 4.061 14 C -0.173 4.173 15 C 0.045 3.955 16 N -0.315 5.315 17 C -0.381 4.381 18 H 0.083 0.917 19 C -0.197 4.197 20 N -0.547 5.547 21 Si 0.731 3.269 22 H -0.215 1.215 23 H -0.218 1.218 24 H -0.204 1.204 25 C 0.018 3.982 26 C -0.164 4.164 27 C -0.217 4.217 28 C -0.216 4.216 29 H 0.088 0.912 30 H 0.090 0.910 31 H 0.094 0.906 32 H 0.073 0.927 33 H 0.070 0.930 34 H 0.133 0.867 35 H 0.114 0.886 36 H 0.115 0.885 37 H 0.116 0.884 38 H 0.114 0.886 39 H 0.244 0.756 40 H 0.097 0.903 41 H 0.069 0.931 42 H 0.070 0.930 43 H 0.105 0.895 44 H 0.037 0.963 45 H 0.064 0.936 46 H 0.035 0.965 47 H 0.072 0.928 48 H 0.071 0.929 49 H 0.075 0.925 50 H 0.110 0.890 51 H 0.123 0.877 52 H 0.122 0.878 53 H 0.110 0.890 Dipole moment (debyes) X Y Z Total from point charges 5.537 -6.865 -19.775 21.653 hybrid contribution -0.659 1.563 0.002 1.696 sum 4.878 -5.302 -19.773 21.045 Atomic orbital electron populations 1.22244 0.93355 1.03947 1.03355 1.21886 0.95211 0.96241 0.85679 0.89730 1.88031 1.86488 1.33907 1.18689 1.22119 0.99267 0.91965 0.85283 1.22099 1.00497 1.01091 0.93279 1.23084 0.90850 1.03626 1.03806 1.22994 1.02570 1.01106 0.94478 1.20311 0.80065 0.85737 0.84108 1.90669 1.39274 1.22943 1.86401 1.45362 1.56992 1.27474 1.04067 1.20934 0.99715 0.90916 0.86125 1.20778 0.96850 0.97117 0.91327 1.22209 0.94702 0.98544 1.01817 1.20765 0.93071 0.92238 0.89398 1.44065 1.00844 1.08396 1.78192 1.19198 0.99646 0.83642 1.35634 0.91711 1.25040 0.99273 0.97346 0.98031 1.70035 1.11283 1.32111 1.41228 0.86295 0.78470 0.85013 0.77121 1.21526 1.21802 1.20428 1.20920 0.88272 0.96142 0.92898 1.22055 0.96611 0.97478 1.00210 1.23437 1.01518 0.93535 1.03253 1.23352 0.91169 1.03674 1.03372 0.91241 0.90997 0.90572 0.92660 0.93035 0.86715 0.88633 0.88544 0.88374 0.88648 0.75627 0.90289 0.93065 0.92967 0.89455 0.96327 0.93649 0.96490 0.92817 0.92910 0.92470 0.89037 0.87744 0.87772 0.88955 Number of geometries 1 Number of calculations of the screened coulomb radii 1 The total number of SCF iterations 16. 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.17 -1.57 9.48 37.16 0.35 -1.22 16 2 C 0.07 0.78 3.13 -27.88 -0.09 0.70 16 3 H 0.08 0.76 7.79 -51.93 -0.40 0.36 16 4 O -0.35 -4.37 9.53 -55.09 -0.52 -4.89 16 5 C 0.09 0.83 6.03 37.16 0.22 1.05 16 6 C -0.15 -1.31 6.06 -90.62 -0.55 -1.86 16 7 C -0.18 -1.34 10.18 -26.69 -0.27 -1.61 16 8 C -0.17 -1.15 10.25 -26.69 -0.27 -1.42 16 9 C 0.51 7.11 6.88 -10.98 -0.08 7.04 16 10 O -0.52 -8.86 11.99 -14.31 -0.17 -9.03 16 11 N -0.69 -9.10 5.35 -47.60 -0.25 -9.35 16 12 C 0.11 1.59 2.48 -131.82 -0.33 1.26 16 13 C -0.04 -0.66 2.28 -90.58 -0.21 -0.86 16 14 C -0.13 -2.31 4.83 -26.61 -0.13 -2.44 16 15 C 0.17 3.86 5.31 -3.71 -0.02 3.84 16 16 N -0.59 -14.17 2.99 -172.32 -0.52 -14.68 16 17 C -0.25 -6.67 9.58 -15.06 -0.14 -6.82 16 18 H 0.07 1.70 7.89 -52.49 -0.41 1.29 16 19 C 0.00 0.05 4.90 -17.43 -0.09 -0.04 16 20 N -0.90 -25.53 11.22 55.49 0.62 -24.91 16 21 Si 0.90 21.12 29.59 -169.99 -5.03 16.09 16 22 H -0.29 -6.71 7.11 56.52 0.40 -6.31 16 23 H -0.29 -6.95 7.11 56.52 0.40 -6.55 16 24 H -0.28 -6.42 7.11 56.52 0.40 -6.02 16 25 C 0.14 2.90 6.52 -3.71 -0.02 2.88 16 26 C -0.12 -2.05 5.22 -26.61 -0.14 -2.19 16 27 C -0.18 -2.60 9.37 -26.69 -0.25 -2.85 16 28 C -0.18 -2.56 8.91 -26.69 -0.24 -2.80 16 29 H 0.07 0.52 7.87 -51.93 -0.41 0.12 16 30 H 0.07 0.73 8.14 -51.93 -0.42 0.31 16 31 H 0.08 0.64 7.92 -51.93 -0.41 0.23 16 32 H 0.06 0.42 7.87 -51.93 -0.41 0.01 16 33 H 0.05 0.43 7.82 -51.93 -0.41 0.03 16 34 H 0.11 1.08 8.14 -51.93 -0.42 0.65 16 35 H 0.10 0.70 8.03 -51.93 -0.42 0.28 16 36 H 0.10 0.68 8.14 -51.93 -0.42 0.25 16 37 H 0.10 0.55 8.14 -51.93 -0.42 0.13 16 38 H 0.10 0.54 8.14 -51.93 -0.42 0.12 16 39 H 0.41 4.91 6.44 -40.82 -0.26 4.65 16 40 H 0.08 1.28 7.90 -51.93 -0.41 0.87 16 41 H 0.05 0.86 6.50 -51.93 -0.34 0.52 16 42 H 0.05 0.84 8.14 -51.93 -0.42 0.42 16 43 H 0.09 2.22 6.06 -51.93 -0.31 1.90 16 44 H 0.02 0.45 8.14 -51.93 -0.42 0.02 16 45 H 0.25 6.95 8.31 -40.82 -0.34 6.61 16 46 H 0.02 0.35 8.14 -51.93 -0.42 -0.07 16 47 H 0.05 1.07 7.93 -51.93 -0.41 0.66 16 48 H 0.05 0.78 7.24 -51.93 -0.38 0.41 16 49 H 0.06 0.97 8.14 -51.93 -0.42 0.55 16 50 H 0.09 1.22 8.05 -51.93 -0.42 0.80 16 51 H 0.10 1.67 6.38 -51.93 -0.33 1.34 16 52 H 0.10 1.66 6.37 -51.93 -0.33 1.33 16 53 H 0.09 1.23 7.60 -51.93 -0.39 0.83 16 LS Contribution 410.64 15.07 6.19 6.19 Total: -1.00 -30.86 410.64 -11.33 -42.19 By element: Atomic # 1 Polarization: 15.14 SS G_CDS: -9.40 Total: 5.74 kcal Atomic # 6 Polarization: -5.09 SS G_CDS: -2.24 Total: -7.34 kcal Atomic # 7 Polarization: -48.80 SS G_CDS: -0.15 Total: -48.95 kcal Atomic # 8 Polarization: -13.23 SS G_CDS: -0.70 Total: -13.93 kcal Atomic # 14 Polarization: 21.12 SS G_CDS: -5.03 Total: 16.09 kcal Total LS contribution 6.19 Total: 6.19 kcal Total: -30.86 -11.33 -42.19 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. ZINC001284844282.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 31.931 kcal (2) G-P(sol) polarization free energy of solvation -30.861 kcal (3) G-ENP(sol) elect.-nuc.-pol. free energy of system 1.071 kcal (4) G-CDS(sol) cavity-dispersion-solvent structure free energy -11.327 kcal (5) G-P-CDS(sol) = G-P(sol) + G-CDS(sol) = (2) + (4) -42.187 kcal (6) G-S(sol) free energy of system = (1) + (5) -10.256 kcal 1SCF WAS SPECIFIED, SO GEOMETRY OPTIMIZATION WAS NOT USED Total computer time = 1.04 seconds