Wall clock time and date at job start Thu Mar 4 2021 12:25:40 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= 0 * AM1 - THE AM1 HAMILTONIAN TO BE USED * SM5.42R - SM5.42R CALCULATIONS WILL BE PERFORMED * - CHARGE MODEL 2 WILL BE USED (CM2A) * - THE SOLVENT IS WATER ******************************************************************************* 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.53004 * 1 3 3 C 1.52998 * 109.46966 * 2 1 4 4 C 1.52993 * 109.46697 * 240.00276 * 2 1 3 5 5 O 1.45195 * 109.47008 * 120.00086 * 2 1 3 6 6 C 1.34637 * 116.99983 * 60.00151 * 5 2 1 7 7 O 1.21507 * 120.00381 * 359.97438 * 6 5 2 8 8 N 1.34774 * 119.99966 * 179.97438 * 6 5 2 9 9 C 1.46504 * 119.99976 * 179.97438 * 8 6 5 10 10 H 1.08991 * 112.84732 * 336.19812 * 9 8 6 11 11 C 1.53791 * 113.68879 * 107.46622 * 9 8 6 12 12 C 1.53785 * 87.07932 * 139.99230 * 11 9 8 13 13 H 1.09000 * 113.69063 * 89.12541 * 12 11 9 14 14 C 1.53002 * 113.61304 * 220.01691 * 12 11 9 15 15 C 1.52997 * 109.46906 * 272.49831 * 14 12 11 16 16 N 1.46505 * 109.46890 * 179.97438 * 15 14 12 17 17 C 1.34776 * 119.99773 * 180.02562 * 16 15 14 18 18 O 1.21549 * 119.99928 * 359.97438 * 17 16 15 19 19 C 1.47844 * 119.99630 * 180.02562 * 17 16 15 20 20 C 1.39549 * 120.14465 * 359.97438 * 19 17 16 21 21 C 1.38064 * 119.84544 * 179.97438 * 20 19 17 22 22 Br 1.89105 * 119.93873 * 179.97438 * 21 20 19 23 23 C 1.38426 * 120.12372 * 0.02562 * 21 20 19 24 24 C 1.38357 * 120.28441 * 359.97355 * 23 21 20 25 25 C 1.50699 * 119.92221 * 180.02562 * 24 23 21 26 26 O 1.42903 * 109.47152 * 120.00572 * 25 24 23 27 27 Si 1.86293 * 109.47431 * 0.02562 * 25 24 23 28 28 C 1.37945 * 120.15763 * 0.30103 * 24 23 21 29 29 C 1.53780 * 113.61410 * 204.99859 * 9 8 6 30 30 H 1.09008 * 109.46861 * 299.99964 * 1 2 3 31 31 H 1.08998 * 109.47435 * 59.99254 * 1 2 3 32 32 H 1.08998 * 109.47026 * 179.97438 * 1 2 3 33 33 H 1.08997 * 109.47002 * 179.97438 * 3 2 1 34 34 H 1.08998 * 109.47653 * 300.00613 * 3 2 1 35 35 H 1.09008 * 109.47301 * 60.00306 * 3 2 1 36 36 H 1.08998 * 109.47184 * 299.99660 * 4 2 1 37 37 H 1.09002 * 109.47607 * 60.00435 * 4 2 1 38 38 H 1.09005 * 109.46789 * 179.97438 * 4 2 1 39 39 H 0.97005 * 120.00636 * 359.97438 * 8 6 5 40 40 H 1.08998 * 113.61383 * 25.44095 * 11 9 8 41 41 H 1.08995 * 113.60881 * 254.54168 * 11 9 8 42 42 H 1.08995 * 109.47244 * 152.49692 * 14 12 11 43 43 H 1.09003 * 109.46605 * 32.49642 * 14 12 11 44 44 H 1.09000 * 109.47196 * 60.00372 * 15 14 12 45 45 H 1.08991 * 109.47104 * 299.99998 * 15 14 12 46 46 H 0.97002 * 119.99860 * 0.02562 * 16 15 14 47 47 H 1.08001 * 120.07154 * 0.02562 * 20 19 17 48 48 H 1.08004 * 119.85649 * 179.97438 * 23 21 20 49 49 H 1.09008 * 109.47393 * 239.99967 * 25 24 23 50 50 H 0.96700 * 114.00311 * 180.02562 * 26 25 24 51 51 H 1.48506 * 109.99973 * 59.99460 * 27 25 24 52 52 H 1.48493 * 110.00375 * 181.31561 * 27 25 24 53 53 H 1.07997 * 120.06717 * 179.70289 * 28 24 23 54 54 H 1.09000 * 113.61522 * 105.38688 * 29 9 8 55 55 H 1.09001 * 113.61663 * 334.48214 * 29 9 8 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 6 2.0400 1.4425 0.0000 4 6 2.0399 -0.7212 -1.2492 5 8 2.0140 -0.6845 1.1855 6 6 1.6522 -1.9723 1.3382 7 8 0.9509 -2.5095 0.5040 8 7 2.0680 -2.6657 2.4165 9 6 1.6739 -4.0669 2.5830 10 1 0.7592 -4.3160 2.0453 11 6 2.8202 -5.0638 2.3436 12 6 2.2715 -5.8609 3.5388 13 1 1.5201 -6.6010 3.2637 14 6 3.3621 -6.4311 4.4478 15 6 3.7985 -7.8006 3.9235 16 7 4.8433 -8.3463 4.7936 17 6 5.3831 -9.5515 4.5242 18 8 5.0027 -10.1858 3.5596 19 6 6.4370 -10.1024 5.4026 20 6 6.8755 -9.3769 6.5110 21 6 7.8596 -9.8975 7.3275 22 35 8.4563 -8.9162 8.8299 23 6 8.4112 -11.1359 7.0477 24 6 7.9804 -11.8597 5.9501 25 6 8.5861 -13.2074 5.6536 26 8 7.5632 -14.2050 5.6773 27 14 9.8724 -13.6070 6.9405 28 6 7.0017 -11.3482 5.1234 29 6 1.6755 -4.5375 4.0470 30 1 -0.3633 0.5139 0.8901 31 1 -0.3634 0.5139 -0.8899 32 1 -0.3633 -1.0276 -0.0005 33 1 3.1300 1.4425 0.0005 34 1 1.6766 1.9564 -0.8899 35 1 1.6767 1.9564 0.8901 36 1 1.6766 -1.7488 -1.2493 37 1 1.6766 -0.2073 -2.1393 38 1 3.1300 -0.7208 -1.2495 39 1 2.6282 -2.2371 3.0825 40 1 3.8098 -4.6413 2.5178 41 1 2.7464 -5.5941 1.3942 42 1 2.9727 -6.5369 5.4603 43 1 4.2177 -5.7558 4.4563 44 1 4.1881 -7.6947 2.9110 45 1 2.9430 -8.4758 3.9149 46 1 5.1465 -7.8402 5.5636 47 1 6.4453 -8.4108 6.7301 48 1 9.1811 -11.5388 7.6891 49 1 9.0514 -13.1877 4.6680 50 1 7.8818 -15.0996 5.4947 51 1 9.2501 -13.6369 8.2885 52 1 10.4997 -14.9193 6.6414 53 1 6.6700 -11.9124 4.2643 54 1 0.6777 -4.6430 4.4729 55 1 2.3522 -3.9712 4.6870 There are 74 doubly filled levels RHF calculation, no. of doubly occupied orbitals= 74 No. of singly occupied orbitals= 1 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). Br: (AM1): M.J.S. DEWAR AND E. G. ZOEBISCH, THEOCHEM, 180, 1 (1988). 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=0 AM1 1SCF TLIMIT=15 GEO-OK SM5.42R & SOLVNT=WATER m_22_16659234_12335870.mol2 55 1SCF WAS SPECIFIED, SO GEOMETRY OPTIMIZATION WAS NOT USED SCF FIELD WAS ACHIEVED Wall clock time and date at job start Thu Mar 4 2021 12:25:40 Heat of formation + Delta-G solvation = -97.806294 kcal Electronic energy + Delta-G solvation = -35369.704237 eV Core-core repulsion = 30420.285692 eV Total energy + Delta-G solvation = -4949.418544 eV No. of doubly occupied orbitals = 74 and no. of open levels = 1 Molecular weight (most abundant/longest-lived isotopes) = 455.119 amu Computer time = 1.85 seconds Orbital eigenvalues (eV) -41.90186 -41.55984 -40.44214 -39.88734 -38.03557 -37.81635 -36.45311 -35.39092 -34.34894 -33.67355 -33.36797 -31.26237 -30.46214 -28.39235 -28.32017 -27.10198 -26.74289 -26.50130 -24.87848 -23.59355 -22.70926 -22.18514 -21.50170 -21.28380 -20.39382 -19.10638 -18.96714 -18.64608 -17.79974 -17.59182 -17.44832 -17.19290 -16.78842 -16.61168 -16.43971 -16.22224 -16.05045 -15.58668 -15.44701 -15.42577 -15.16614 -14.90112 -14.84157 -14.51170 -14.27203 -14.02667 -13.80076 -13.72809 -13.59311 -13.38088 -13.36724 -13.17388 -12.93696 -12.88465 -12.84335 -12.83442 -12.62272 -12.42767 -12.36781 -12.29337 -12.01791 -11.87222 -11.83774 -11.68826 -11.44194 -11.40887 -11.22419 -11.13853 -11.09134 -10.64967 -10.37428 -10.12608 -10.09288 -9.89024 -4.98639 -0.51663 -0.15431 0.42396 1.18309 1.19546 1.52032 1.55841 1.73136 1.89406 1.99398 2.45194 3.14048 3.16013 3.39230 3.49399 3.52158 3.57071 3.72457 3.76987 3.78328 3.83995 3.85786 3.87968 3.92566 3.98864 4.09691 4.17194 4.27423 4.33033 4.36137 4.49373 4.52246 4.58744 4.63784 4.67206 4.74121 4.75862 4.77615 4.79583 4.83535 4.85656 4.89010 4.93235 4.97114 4.99929 5.04719 5.09758 5.13506 5.21626 5.24660 5.27360 5.30602 5.34519 5.38964 5.56113 5.61765 6.36739 6.41373 6.52251 7.03862 7.34441 Net atomic charges, atomic populations, and dipole contributions calculated with CM2 Atom NO. Type Charge No. of electrons 1 C -0.180 4.180 2 C 0.128 3.872 3 C -0.138 4.138 4 C -0.180 4.180 5 O -0.370 6.370 6 C 0.641 3.359 7 O -0.597 6.597 8 N -0.688 5.688 9 C 0.128 3.872 10 H 0.097 0.903 11 C -0.149 4.149 12 C -0.098 4.098 13 H 0.096 0.904 14 C -0.108 4.108 15 C 0.117 3.883 16 N -0.709 5.709 17 C 0.557 3.443 18 O -0.576 6.576 19 C -0.137 4.137 20 C -0.044 4.044 21 C -0.099 4.099 22 Br -0.025 7.025 23 C -0.063 4.063 24 C -0.079 4.079 25 C -0.088 4.088 26 O -0.557 6.557 27 Si 0.758 3.242 28 C -0.032 4.032 29 C -0.124 4.124 30 H 0.062 0.938 31 H 0.086 0.914 32 H 0.062 0.938 33 H 0.070 0.930 34 H 0.098 0.902 35 H 0.070 0.930 36 H 0.062 0.938 37 H 0.086 0.914 38 H 0.062 0.938 39 H 0.418 0.582 40 H 0.097 0.903 41 H 0.059 0.941 42 H 0.091 0.909 43 H 0.088 0.912 44 H 0.059 0.941 45 H 0.059 0.941 46 H 0.416 0.584 47 H 0.167 0.833 48 H 0.155 0.845 49 H 0.102 0.898 50 H 0.398 0.602 51 H -0.290 1.290 52 H -0.279 1.279 53 H 0.130 0.870 54 H 0.083 0.917 55 H 0.110 0.890 For the charges calculated by CM2: Dipole moment (debyes) X Y Z Total from point charges 4.799 4.867 4.430 8.145 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.238 4.238 2 C 0.091 3.909 3 C -0.195 4.195 4 C -0.237 4.237 5 O -0.284 6.284 6 C 0.394 3.606 7 O -0.488 6.488 8 N -0.347 5.347 9 C 0.024 3.976 10 H 0.115 0.885 11 C -0.187 4.187 12 C -0.117 4.117 13 H 0.114 0.886 14 C -0.145 4.145 15 C -0.006 4.006 16 N -0.362 5.362 17 C 0.345 3.655 18 O -0.456 6.456 19 C -0.141 4.141 20 C -0.064 4.064 21 C -0.178 4.178 22 Br 0.061 6.939 23 C -0.083 4.083 24 C -0.080 4.080 25 C -0.196 4.196 26 O -0.368 6.368 27 Si 0.662 3.338 28 C -0.051 4.051 29 C -0.162 4.162 30 H 0.081 0.919 31 H 0.105 0.895 32 H 0.081 0.919 33 H 0.089 0.911 34 H 0.116 0.884 35 H 0.089 0.911 36 H 0.081 0.919 37 H 0.105 0.895 38 H 0.081 0.919 39 H 0.255 0.745 40 H 0.115 0.885 41 H 0.078 0.922 42 H 0.109 0.891 43 H 0.107 0.893 44 H 0.077 0.923 45 H 0.078 0.922 46 H 0.254 0.746 47 H 0.184 0.816 48 H 0.172 0.828 49 H 0.120 0.880 50 H 0.249 0.751 51 H -0.217 1.217 52 H -0.205 1.205 53 H 0.148 0.852 54 H 0.102 0.898 55 H 0.128 0.872 Dipole moment (debyes) X Y Z Total from point charges 4.261 4.622 4.598 7.788 hybrid contribution -0.688 -1.224 -0.635 1.541 sum 3.573 3.398 3.963 6.326 Atomic orbital electron populations 1.22478 0.94118 1.03633 1.03528 1.22331 0.94409 0.92610 0.81526 1.21891 1.02474 0.90831 1.04266 1.22472 1.01630 1.02197 0.97446 1.86309 1.71545 1.23378 1.47195 1.18487 0.78620 0.80699 0.82749 1.90901 1.42469 1.71095 1.44375 1.44805 1.54670 1.12365 1.22899 1.22660 0.98767 0.77044 0.99094 0.88477 1.23636 0.99119 0.96939 0.99049 1.22195 0.98003 0.95995 0.95526 0.88565 1.21696 0.98922 0.93845 1.00069 1.21465 0.89740 0.95556 0.93818 1.45822 1.40070 1.18770 1.31505 1.18331 0.81553 0.81575 0.84044 1.90799 1.63361 1.60610 1.30859 1.19893 0.99756 0.96796 0.97668 1.21471 0.93714 0.99946 0.91267 1.21705 1.01663 0.97048 0.97424 1.96848 1.88666 1.72006 1.36407 1.21322 0.98450 0.93007 0.95479 1.19357 0.97886 0.95615 0.95158 1.25049 0.97110 0.86801 1.10688 1.86542 1.35429 1.20671 1.94143 0.94289 0.79949 0.84779 0.74810 1.21259 0.92464 0.93187 0.98203 1.23393 1.03095 0.95555 0.94134 0.91859 0.89511 0.91895 0.91142 0.88367 0.91100 0.91931 0.89509 0.91940 0.74501 0.88491 0.92244 0.89096 0.89301 0.92306 0.92239 0.74598 0.81612 0.82804 0.87989 0.75056 1.21742 1.20538 0.85193 0.89848 0.87207 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.18 -0.66 8.37 71.98 0.60 -0.06 16 2 C 0.13 0.42 1.13 -10.80 -0.01 0.40 16 3 C -0.14 0.28 8.85 71.98 0.64 0.91 16 4 C -0.18 -0.70 8.37 71.98 0.60 -0.09 16 5 O -0.37 -2.05 9.94 -55.46 -0.55 -2.60 16 6 C 0.64 4.39 8.06 129.79 1.05 5.44 16 7 O -0.60 -7.63 11.48 19.82 0.23 -7.40 16 8 N -0.69 0.23 5.36 -567.60 -3.04 -2.82 16 9 C 0.13 -0.25 4.24 45.52 0.19 -0.06 16 10 H 0.10 0.27 7.46 -2.39 -0.02 0.25 16 11 C -0.15 0.56 7.43 31.13 0.23 0.79 16 12 C -0.10 0.70 3.65 -10.26 -0.04 0.66 16 13 H 0.10 -0.57 8.13 -2.39 -0.02 -0.59 16 14 C -0.11 1.01 5.17 30.59 0.16 1.16 16 15 C 0.12 -0.33 5.66 86.38 0.49 0.16 16 16 N -0.71 1.80 5.51 -465.98 -2.57 -0.76 16 17 C 0.56 2.92 7.77 86.79 0.67 3.59 16 18 O -0.58 -8.26 16.35 -3.87 -0.06 -8.32 16 19 C -0.14 -0.34 5.87 -20.09 -0.12 -0.46 16 20 C -0.04 0.11 9.28 22.54 0.21 0.32 16 21 C -0.10 0.15 6.22 22.27 0.14 0.28 16 22 Br -0.02 0.06 33.69 -20.37 -0.69 -0.62 16 23 C -0.06 0.05 8.41 22.34 0.19 0.24 16 24 C -0.08 -0.12 5.42 -19.91 -0.11 -0.22 16 25 C -0.09 0.11 2.81 71.24 0.20 0.31 16 26 O -0.56 -1.14 12.74 -148.98 -1.90 -3.04 16 27 Si 0.76 -2.52 29.78 68.60 2.04 -0.48 16 28 C -0.03 -0.15 9.16 22.53 0.21 0.06 16 29 C -0.12 1.00 7.39 31.12 0.23 1.23 16 30 H 0.06 0.17 8.14 -2.38 -0.02 0.15 16 31 H 0.09 0.02 8.14 -2.39 -0.02 0.00 16 32 H 0.06 0.58 5.88 -2.39 -0.01 0.56 16 33 H 0.07 -0.16 8.14 -2.39 -0.02 -0.18 16 34 H 0.10 -0.51 8.14 -2.39 -0.02 -0.53 16 35 H 0.07 -0.17 8.14 -2.38 -0.02 -0.19 16 36 H 0.06 0.59 5.88 -2.39 -0.01 0.58 16 37 H 0.09 0.04 8.14 -2.39 -0.02 0.02 16 38 H 0.06 0.19 8.14 -2.38 -0.02 0.17 16 39 H 0.42 -2.30 8.84 -92.70 -0.82 -3.12 16 40 H 0.10 -0.52 7.90 -2.39 -0.02 -0.54 16 41 H 0.06 0.00 8.14 -2.39 -0.02 -0.02 16 42 H 0.09 -1.05 8.14 -2.39 -0.02 -1.07 16 43 H 0.09 -1.03 7.84 -2.39 -0.02 -1.05 16 44 H 0.06 0.03 8.14 -2.39 -0.02 0.01 16 45 H 0.06 0.01 8.14 -2.39 -0.02 -0.01 16 46 H 0.42 -3.86 6.51 -92.71 -0.60 -4.46 16 47 H 0.17 -1.52 6.40 -2.91 -0.02 -1.54 16 48 H 0.15 -0.61 4.86 -2.91 -0.01 -0.63 16 49 H 0.10 -0.52 7.96 -2.38 -0.02 -0.54 16 50 H 0.40 -2.43 8.74 -74.06 -0.65 -3.08 16 51 H -0.29 -0.99 7.11 99.48 0.71 -0.28 16 52 H -0.28 -0.18 7.11 99.48 0.71 0.52 16 53 H 0.13 0.97 7.65 -2.91 -0.02 0.95 16 54 H 0.08 -0.66 8.14 -2.39 -0.02 -0.67 16 55 H 0.11 -1.26 8.08 -2.39 -0.02 -1.28 16 Total: 0.00 -25.80 462.16 -2.09 -27.90 By element: Atomic # 1 Polarization: -15.45 SS G_CDS: -1.09 Total: -16.54 kcal Atomic # 6 Polarization: 9.15 SS G_CDS: 5.53 Total: 14.68 kcal Atomic # 7 Polarization: 2.03 SS G_CDS: -5.61 Total: -3.58 kcal Atomic # 8 Polarization: -19.07 SS G_CDS: -2.29 Total: -21.36 kcal Atomic # 14 Polarization: -2.52 SS G_CDS: 2.04 Total: -0.48 kcal Atomic # 35 Polarization: 0.06 SS G_CDS: -0.69 Total: -0.62 kcal Total: -25.80 -2.09 -27.90 kcal The number of atoms in the molecule is 55 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. m_22_16659234_12335870.mol2 55 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 -69.911 kcal (2) G-P(sol) polarization free energy of solvation -25.802 kcal (3) G-ENP(sol) elect.-nuc.-pol. free energy of system -95.713 kcal (4) G-CDS(sol) cavity-dispersion-solvent structure free energy -2.093 kcal (5) G-P-CDS(sol) = G-P(sol) + G-CDS(sol) = (2) + (4) -27.895 kcal (6) G-S(sol) free energy of system = (1) + (5) -97.806 kcal 1SCF WAS SPECIFIED, SO GEOMETRY OPTIMIZATION WAS NOT USED Total computer time = 1.85 seconds