Wall clock time and date at job start Sat Apr 11 2020 06:10:11 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.42899 * 1 3 3 C 1.42906 * 114.00376 * 2 1 4 4 C 1.52998 * 109.47482 * 179.97438 * 3 2 1 5 5 C 1.50698 * 109.54558 * 300.18316 * 4 3 2 6 6 O 1.21283 * 119.99965 * 239.82464 * 5 4 3 7 7 N 1.34783 * 120.00232 * 59.81539 * 5 4 3 8 8 C 1.46498 * 119.99863 * 180.02562 * 7 5 4 9 9 C 1.53002 * 109.46932 * 180.02562 * 8 7 5 10 10 H 1.09006 * 110.88741 * 51.11705 * 9 8 7 11 11 C 1.54906 * 111.00302 * 287.20080 * 9 8 7 12 12 N 1.47430 * 104.83403 * 204.91701 * 11 9 8 13 13 C 1.36928 * 125.65400 * 204.10958 * 12 11 9 14 14 H 1.08007 * 119.99519 * 359.97438 * 13 12 11 15 15 C 1.38813 * 120.00359 * 179.97438 * 13 12 11 16 16 N 1.31611 * 119.99938 * 359.97438 * 15 13 12 17 17 Si 1.86798 * 119.99950 * 179.97438 * 15 13 12 18 18 H 1.48495 * 109.47528 * 359.97438 * 17 15 13 19 19 H 1.48499 * 109.47206 * 120.00055 * 17 15 13 20 20 H 1.48504 * 109.46869 * 239.99758 * 17 15 13 21 21 C 1.47027 * 108.70056 * 24.12443 * 12 11 9 22 22 C 1.54330 * 107.26902 * 358.94279 * 21 12 11 23 23 H 1.08998 * 110.86756 * 96.24433 * 22 21 12 24 24 C 1.53001 * 110.71777 * 219.44044 * 22 21 12 25 25 C 1.53177 * 109.61441 * 60.38966 * 4 3 2 26 26 C 1.53092 * 109.16443 * 176.87824 * 25 4 3 27 27 O 1.42905 * 109.42428 * 302.41109 * 26 25 4 28 28 C 1.42901 * 114.09770 * 61.14988 * 27 26 25 29 29 C 1.53095 * 109.41489 * 298.84274 * 28 27 26 30 30 H 1.08991 * 109.46979 * 179.97438 * 1 2 3 31 31 H 1.08999 * 109.47226 * 300.00053 * 1 2 3 32 32 H 1.09004 * 109.46785 * 60.00007 * 1 2 3 33 33 H 1.09007 * 109.46596 * 60.00040 * 3 2 1 34 34 H 1.08995 * 109.47101 * 300.00127 * 3 2 1 35 35 H 0.97005 * 119.99507 * 0.02562 * 7 5 4 36 36 H 1.08992 * 109.47380 * 300.00398 * 8 7 5 37 37 H 1.08995 * 109.46824 * 60.00549 * 8 7 5 38 38 H 1.08997 * 110.36842 * 323.74908 * 11 9 8 39 39 H 1.08997 * 110.37446 * 86.07643 * 11 9 8 40 40 H 0.97004 * 120.00200 * 180.02562 * 16 15 13 41 41 H 1.08999 * 109.88217 * 239.53014 * 21 12 11 42 42 H 1.08992 * 109.88354 * 118.36115 * 21 12 11 43 43 H 1.09005 * 109.46827 * 173.52062 * 24 22 21 44 44 H 1.08999 * 109.47809 * 293.52590 * 24 22 21 45 45 H 1.09005 * 109.46929 * 53.52780 * 24 22 21 46 46 H 1.09004 * 109.51365 * 56.98141 * 25 4 3 47 47 H 1.08998 * 109.53811 * 296.79585 * 25 4 3 48 48 H 1.09002 * 109.48343 * 182.42096 * 26 25 4 49 49 H 1.08999 * 109.47546 * 62.39689 * 26 25 4 50 50 H 1.09002 * 109.48282 * 58.82206 * 28 27 26 51 51 H 1.08997 * 109.48587 * 178.84728 * 28 27 26 52 52 H 1.09001 * 109.52141 * 177.51212 * 29 28 27 53 53 H 1.09000 * 109.52256 * 297.64026 * 29 28 27 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.0103 1.3055 0.0000 4 6 3.5355 1.1847 0.0006 5 6 3.9812 0.4333 1.2286 6 8 4.7255 0.9627 2.0266 7 7 3.5527 -0.8270 1.4401 8 6 3.9856 -1.5571 2.6341 9 6 3.3612 -2.9539 2.6320 10 1 3.5362 -3.4564 1.6807 11 6 1.8476 -2.8859 2.9547 12 7 1.5209 -4.2105 3.5133 13 6 0.2813 -4.7922 3.5154 14 1 -0.5599 -4.2721 3.0813 15 6 0.1058 -6.0503 4.0752 16 7 1.1310 -6.6842 4.6036 17 14 -1.5854 -6.8434 4.0785 18 1 -2.5612 -5.9315 3.4295 19 1 -1.5324 -8.1255 3.3311 20 1 -2.0073 -7.1037 5.4784 21 6 2.7465 -4.8035 4.0683 22 6 3.8943 -3.8031 3.8160 23 1 4.0647 -3.1784 4.6928 24 6 5.1763 -4.5387 3.4205 25 6 3.9923 0.4331 -1.2535 26 6 5.5223 0.3862 -1.2819 27 8 6.0393 1.7181 -1.2489 28 6 5.6838 2.4530 -0.0761 29 6 4.1603 2.5832 0.0007 30 1 -0.3633 -1.0276 -0.0005 31 1 -0.3633 0.5138 0.8900 32 1 -0.3633 0.5139 -0.8900 33 1 1.6887 1.8464 0.8901 34 1 1.6887 1.8464 -0.8899 35 1 2.9577 -1.2504 0.8017 36 1 5.0720 -1.6452 2.6322 37 1 3.6664 -1.0174 3.5256 38 1 1.6511 -2.1054 3.6898 39 1 1.2736 -2.7054 2.0459 40 1 1.0085 -7.5635 4.9945 41 1 2.9610 -5.7477 3.5678 42 1 2.6268 -4.9679 5.1391 43 1 5.9385 -3.8132 3.1363 44 1 5.5326 -5.1280 4.2655 45 1 4.9710 -5.1991 2.5779 46 1 3.6273 0.9501 -2.1410 47 1 3.5971 -0.5825 -1.2341 48 1 5.8539 -0.1099 -2.1941 49 1 5.8851 -0.1662 -0.4152 50 1 6.0486 1.9277 0.8066 51 1 6.1319 3.4457 -0.1193 52 1 3.8850 3.1051 0.9172 53 1 3.7989 3.1446 -0.8609 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 ZINC001102027882.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:11 Heat of formation + Delta-G solvation = -114.063350 kcal Electronic energy + Delta-G solvation = -32280.890192 eV Core-core repulsion = 28152.981775 eV Total energy + Delta-G solvation = -4127.908417 eV No. of doubly occupied orbitals = 66 Molecular weight (most abundant/longest-lived isotopes) = 340.220 amu Computer time = 0.75 seconds Orbital eigenvalues (eV) -41.07830 -38.56208 -37.96317 -37.70858 -34.63350 -34.36223 -32.89825 -31.28865 -30.71074 -29.29431 -28.36240 -26.27683 -25.97567 -25.19166 -23.82885 -23.56415 -22.71817 -21.77487 -20.03320 -19.80303 -18.88210 -17.14346 -16.80222 -16.68936 -16.51244 -16.09478 -15.85013 -15.63071 -15.38680 -14.90239 -14.66166 -14.57756 -14.24436 -13.78598 -13.70644 -13.62813 -13.19391 -12.89463 -12.65751 -12.54579 -12.45163 -12.35514 -12.15201 -11.79038 -11.67031 -11.61885 -11.45592 -11.35714 -11.18287 -10.87309 -10.81523 -10.76677 -10.68121 -10.13008 -10.01823 -9.96238 -9.91000 -9.86026 -9.72258 -9.58348 -9.32954 -8.80736 -8.60120 -6.81499 -5.70692 -3.03781 2.62631 3.23480 3.38708 3.42189 3.43810 3.46567 3.74152 4.29947 4.39385 4.49185 4.58002 4.79442 4.82505 5.00886 5.03238 5.10904 5.15805 5.21301 5.25798 5.27034 5.35880 5.37266 5.40646 5.44367 5.53701 5.55248 5.58784 5.60587 5.64530 5.79511 5.79972 5.82406 5.88185 5.95992 6.05801 6.16122 6.21701 6.32632 6.50751 6.58567 6.62566 6.66130 6.80532 6.83331 6.98526 7.04219 7.57578 7.66143 7.70587 7.84728 8.26900 8.34106 9.22857 9.83468 10.50237 11.40964 Molecular weight = 340.22amu Principal moments of inertia in cm(-1) A = 0.015601 B = 0.003084 C = 0.002729 Principal moments of inertia in units of 10**(-40)*gram-cm**2 A = 1794.285376 B = 9077.004649 C =10257.006323 Net atomic charges, atomic populations, and dipole contributions calculated with CM2 Atom NO. Type Charge No. of electrons 1 C 0.035 3.965 2 O -0.389 6.389 3 C 0.074 3.926 4 C -0.078 4.078 5 C 0.522 3.478 6 O -0.548 6.548 7 N -0.708 5.708 8 C 0.131 3.869 9 C -0.111 4.111 10 H 0.079 0.921 11 C 0.145 3.855 12 N -0.607 5.607 13 C -0.226 4.226 14 H 0.077 0.923 15 C -0.009 4.009 16 N -0.921 5.921 17 Si 0.910 3.090 18 H -0.281 1.281 19 H -0.291 1.291 20 H -0.291 1.291 21 C 0.184 3.816 22 C -0.111 4.111 23 H 0.066 0.934 24 C -0.134 4.134 25 C -0.139 4.139 26 C 0.058 3.942 27 O -0.385 6.385 28 C 0.051 3.949 29 C -0.127 4.127 30 H 0.093 0.907 31 H 0.043 0.957 32 H 0.041 0.959 33 H 0.056 0.944 34 H 0.061 0.939 35 H 0.410 0.590 36 H 0.063 0.937 37 H 0.065 0.935 38 H 0.016 0.984 39 H 0.037 0.963 40 H 0.251 0.749 41 H 0.049 0.951 42 H 0.040 0.960 43 H 0.044 0.956 44 H 0.053 0.947 45 H 0.054 0.946 46 H 0.081 0.919 47 H 0.087 0.913 48 H 0.097 0.903 49 H 0.051 0.949 50 H 0.079 0.921 51 H 0.092 0.908 52 H 0.082 0.918 53 H 0.077 0.923 For the charges calculated by CM2: Dipole moment (debyes) X Y Z Total from point charges 5.824 16.686 -12.129 21.435 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.061 4.061 2 O -0.308 6.308 3 C -0.003 4.003 4 C -0.082 4.082 5 C 0.308 3.692 6 O -0.427 6.427 7 N -0.360 5.360 8 C 0.008 3.992 9 C -0.130 4.130 10 H 0.097 0.903 11 C 0.019 3.981 12 N -0.333 5.333 13 C -0.354 4.354 14 H 0.095 0.905 15 C -0.205 4.205 16 N -0.569 5.569 17 Si 0.741 3.259 18 H -0.206 1.206 19 H -0.218 1.218 20 H -0.218 1.218 21 C 0.060 3.940 22 C -0.131 4.131 23 H 0.085 0.915 24 C -0.191 4.191 25 C -0.177 4.177 26 C -0.019 4.019 27 O -0.304 6.304 28 C -0.025 4.025 29 C -0.165 4.165 30 H 0.112 0.888 31 H 0.062 0.938 32 H 0.059 0.941 33 H 0.074 0.926 34 H 0.080 0.920 35 H 0.248 0.752 36 H 0.081 0.919 37 H 0.084 0.916 38 H 0.034 0.966 39 H 0.055 0.945 40 H 0.061 0.939 41 H 0.067 0.933 42 H 0.059 0.941 43 H 0.063 0.937 44 H 0.072 0.928 45 H 0.073 0.927 46 H 0.100 0.900 47 H 0.105 0.895 48 H 0.115 0.885 49 H 0.069 0.931 50 H 0.097 0.903 51 H 0.110 0.890 52 H 0.101 0.899 53 H 0.095 0.905 Dipole moment (debyes) X Y Z Total from point charges 6.447 16.734 -11.617 21.367 hybrid contribution -1.765 -0.587 0.055 1.861 sum 4.682 16.148 -11.562 20.405 Atomic orbital electron populations 1.22777 0.79860 1.03622 0.99888 1.87720 1.20609 1.27067 1.95405 1.22476 0.91010 0.85629 1.01167 1.20959 0.95006 0.97203 0.94986 1.20183 0.80403 0.81989 0.86656 1.90632 1.36770 1.68159 1.47151 1.46175 1.49709 1.14268 1.25880 1.21388 0.98758 0.91733 0.87367 1.22642 0.97473 0.93263 0.99661 0.90287 1.21540 0.91558 0.89119 0.95905 1.44479 1.05744 1.14782 1.68267 1.19082 0.84950 0.98813 1.32604 0.90527 1.24793 0.97531 0.97319 1.00887 1.69971 1.28276 1.14492 1.44117 0.86067 0.83390 0.79184 0.77248 1.20588 1.21794 1.21814 1.21209 0.82038 0.92660 0.98133 1.22443 0.96599 0.97753 0.96301 0.91525 1.21692 0.96518 0.99622 1.01305 1.22207 0.97131 1.01924 0.96405 1.22596 0.92422 0.84440 1.02429 1.87894 1.74611 1.25531 1.42322 1.22812 0.93435 0.97169 0.89112 1.21962 0.97278 0.93915 1.03319 0.88841 0.93829 0.94084 0.92582 0.92049 0.75221 0.91862 0.91633 0.96596 0.94483 0.93850 0.93325 0.94141 0.93717 0.92807 0.92723 0.89993 0.89478 0.88505 0.93102 0.90293 0.88997 0.89928 0.90477 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.03 0.31 11.01 101.05 1.11 1.42 16 2 O -0.39 -4.07 7.61 -35.23 -0.27 -4.34 16 3 C 0.07 0.69 4.81 37.16 0.18 0.86 16 4 C -0.08 -0.80 0.53 -155.48 -0.08 -0.89 16 5 C 0.52 6.99 5.14 -10.98 -0.06 6.93 16 6 O -0.55 -8.84 13.84 5.55 0.08 -8.76 16 7 N -0.71 -9.51 5.02 -60.30 -0.30 -9.82 16 8 C 0.13 1.96 5.13 -4.04 -0.02 1.94 16 9 C -0.11 -1.86 2.50 -88.48 -0.22 -2.08 16 10 H 0.08 1.33 8.05 -51.93 -0.42 0.91 16 11 C 0.15 2.93 6.65 -2.53 -0.02 2.91 16 12 N -0.61 -15.29 3.39 -170.07 -0.58 -15.87 16 13 C -0.23 -6.24 10.28 -15.06 -0.15 -6.39 16 14 H 0.08 2.04 7.83 -52.48 -0.41 1.63 16 15 C -0.01 -0.26 5.49 -17.43 -0.10 -0.35 16 16 N -0.92 -26.91 10.31 55.49 0.57 -26.34 16 17 Si 0.91 21.94 29.55 -169.99 -5.02 16.92 16 18 H -0.28 -6.53 7.11 56.52 0.40 -6.12 16 19 H -0.29 -7.05 7.11 56.52 0.40 -6.65 16 20 H -0.29 -7.09 7.11 56.52 0.40 -6.69 16 21 C 0.18 4.42 4.97 -3.06 -0.02 4.40 16 22 C -0.11 -2.03 3.02 -88.79 -0.27 -2.30 16 23 H 0.07 1.19 8.04 -51.93 -0.42 0.77 16 24 C -0.13 -2.10 9.16 37.16 0.34 -1.76 16 25 C -0.14 -1.28 4.67 -26.59 -0.12 -1.41 16 26 C 0.06 0.59 6.62 37.20 0.25 0.84 16 27 O -0.38 -4.46 10.77 -35.23 -0.38 -4.84 16 28 C 0.05 0.55 6.55 37.21 0.24 0.79 16 29 C -0.13 -1.20 4.56 -26.59 -0.12 -1.33 16 30 H 0.09 0.86 8.14 -51.93 -0.42 0.44 16 31 H 0.04 0.38 8.14 -51.93 -0.42 -0.04 16 32 H 0.04 0.29 8.14 -51.93 -0.42 -0.13 16 33 H 0.06 0.53 8.14 -51.93 -0.42 0.11 16 34 H 0.06 0.45 8.14 -51.93 -0.42 0.03 16 35 H 0.41 5.28 5.84 -40.82 -0.24 5.05 16 36 H 0.06 0.93 8.10 -51.93 -0.42 0.51 16 37 H 0.07 1.06 7.96 -51.93 -0.41 0.65 16 38 H 0.02 0.34 7.91 -51.93 -0.41 -0.07 16 39 H 0.04 0.72 8.14 -51.93 -0.42 0.30 16 40 H 0.25 7.13 8.31 -40.82 -0.34 6.79 16 41 H 0.05 1.29 7.00 -51.93 -0.36 0.92 16 42 H 0.04 1.09 7.51 -51.93 -0.39 0.70 16 43 H 0.04 0.61 8.10 -51.93 -0.42 0.19 16 44 H 0.05 0.82 8.14 -51.93 -0.42 0.40 16 45 H 0.05 0.88 8.06 -51.93 -0.42 0.46 16 46 H 0.08 0.67 8.14 -51.93 -0.42 0.25 16 47 H 0.09 0.81 7.18 -51.93 -0.37 0.44 16 48 H 0.10 0.81 8.14 -51.93 -0.42 0.39 16 49 H 0.05 0.59 7.38 -51.93 -0.38 0.20 16 50 H 0.08 1.02 5.80 -51.93 -0.30 0.72 16 51 H 0.09 0.82 8.14 -51.93 -0.42 0.40 16 52 H 0.08 0.80 8.03 -51.93 -0.42 0.38 16 53 H 0.08 0.60 8.14 -51.93 -0.42 0.18 16 LS Contribution 403.56 15.07 6.08 6.08 Total: -1.00 -31.82 403.56 -8.46 -40.28 By element: Atomic # 1 Polarization: 12.68 SS G_CDS: -9.58 Total: 3.10 kcal Atomic # 6 Polarization: 2.64 SS G_CDS: 0.95 Total: 3.59 kcal Atomic # 7 Polarization: -51.72 SS G_CDS: -0.31 Total: -52.03 kcal Atomic # 8 Polarization: -17.37 SS G_CDS: -0.57 Total: -17.94 kcal Atomic # 14 Polarization: 21.94 SS G_CDS: -5.02 Total: 16.92 kcal Total LS contribution 6.08 Total: 6.08 kcal Total: -31.82 -8.46 -40.28 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. ZINC001102027882.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 -73.783 kcal (2) G-P(sol) polarization free energy of solvation -31.823 kcal (3) G-ENP(sol) elect.-nuc.-pol. free energy of system -105.607 kcal (4) G-CDS(sol) cavity-dispersion-solvent structure free energy -8.457 kcal (5) G-P-CDS(sol) = G-P(sol) + G-CDS(sol) = (2) + (4) -40.280 kcal (6) G-S(sol) free energy of system = (1) + (5) -114.063 kcal 1SCF WAS SPECIFIED, SO GEOMETRY OPTIMIZATION WAS NOT USED Total computer time = 0.75 seconds