Wall clock time and date at job start Wed Apr 1 2020 00:20:51 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.52997 * 1 3 3 H 1.09001 * 109.52374 * 2 1 4 4 C 1.53100 * 109.52190 * 239.80900 * 2 1 3 5 5 N 1.47026 * 109.41899 * 182.54042 * 4 2 1 6 6 C 1.46893 * 111.00090 * 182.98831 * 5 4 2 7 7 C 1.50695 * 109.47294 * 66.04478 * 6 5 4 8 8 H 1.08000 * 120.00373 * 0.02562 * 7 6 5 9 9 C 1.37883 * 120.00368 * 180.02562 * 7 6 5 10 10 N 1.31511 * 120.00323 * 180.02562 * 9 7 6 11 11 Si 1.86800 * 119.99775 * 0.02562 * 9 7 6 12 12 H 1.48496 * 109.47368 * 90.00001 * 11 9 7 13 13 H 1.48503 * 109.47076 * 210.00175 * 11 9 7 14 14 H 1.48508 * 109.46976 * 329.99955 * 11 9 7 15 15 C 1.47178 * 110.63374 * 59.29369 * 5 4 2 16 16 C 1.52168 * 109.35489 * 302.88530 * 15 5 4 17 17 C 1.54384 * 110.46237 * 291.05209 * 16 15 5 18 18 C 1.54458 * 102.88112 * 207.45464 * 17 16 15 19 19 N 1.47077 * 107.28961 * 21.51131 * 18 17 16 20 20 C 1.34780 * 125.62532 * 180.96693 * 19 18 17 21 21 O 1.21513 * 119.99577 * 180.02562 * 20 19 18 22 22 O 1.34628 * 120.00286 * 0.02562 * 20 19 18 23 23 C 1.45202 * 117.00407 * 179.97438 * 22 20 19 24 24 C 1.53006 * 109.46867 * 59.99955 * 23 22 20 25 25 C 1.52999 * 109.47019 * 300.00782 * 23 22 20 26 26 C 1.52997 * 109.47464 * 180.02562 * 23 22 20 27 27 C 1.47573 * 108.74291 * 1.23215 * 19 18 17 28 28 O 1.41958 * 111.47459 * 54.48517 * 16 15 5 29 29 H 1.09005 * 109.47222 * 299.86870 * 1 2 3 30 30 H 1.08994 * 109.47506 * 59.87279 * 1 2 3 31 31 H 1.09006 * 109.47144 * 179.87287 * 1 2 3 32 32 H 1.09000 * 109.48226 * 302.52447 * 4 2 1 33 33 H 1.08998 * 109.47980 * 62.55081 * 4 2 1 34 34 H 1.09000 * 109.47604 * 186.05093 * 6 5 4 35 35 H 1.09009 * 109.46935 * 306.04388 * 6 5 4 36 36 H 0.96996 * 120.00351 * 179.97438 * 10 9 7 37 37 H 1.09000 * 109.49011 * 182.91709 * 15 5 4 38 38 H 1.09004 * 109.48695 * 62.84933 * 15 5 4 39 39 H 1.09004 * 110.73308 * 89.10179 * 17 16 15 40 40 H 1.09001 * 110.88523 * 325.87137 * 17 16 15 41 41 H 1.08998 * 109.87538 * 262.08555 * 18 17 16 42 42 H 1.09004 * 109.87933 * 140.93535 * 18 17 16 43 43 H 1.09001 * 109.46801 * 300.00397 * 24 23 22 44 44 H 1.09001 * 109.46998 * 59.99993 * 24 23 22 45 45 H 1.09000 * 109.47050 * 180.02562 * 24 23 22 46 46 H 1.09000 * 109.47160 * 179.97438 * 25 23 22 47 47 H 1.09002 * 109.46719 * 299.99410 * 25 23 22 48 48 H 1.08993 * 109.47839 * 59.99251 * 25 23 22 49 49 H 1.09001 * 109.47375 * 300.00043 * 26 23 22 50 50 H 1.09002 * 109.46777 * 60.00599 * 26 23 22 51 51 H 1.09000 * 109.47140 * 180.02562 * 26 23 22 52 52 H 1.08999 * 110.37805 * 95.25729 * 27 19 18 53 53 H 1.09004 * 110.37220 * 217.60393 * 27 19 18 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.8942 1.0273 0.0000 4 6 2.0416 -0.7257 -1.2473 5 7 3.5106 -0.7777 -1.2145 6 6 4.0424 -1.4052 -2.4315 7 6 3.7411 -0.5301 -3.6207 8 1 3.2082 0.3995 -3.4860 9 6 4.1458 -0.9162 -4.8810 10 7 3.8833 -0.1523 -5.9188 11 14 5.0669 -2.5245 -5.1143 12 1 4.0954 -3.6200 -5.3614 13 1 5.9846 -2.4054 -6.2758 14 1 5.8531 -2.8266 -3.8912 15 6 3.9778 -1.4828 -0.0100 16 6 3.4224 -0.7913 1.2265 17 6 4.0730 0.5972 1.4064 18 6 3.9897 0.8212 2.9324 19 7 3.8457 -0.5006 3.5611 20 6 3.7219 -0.7318 4.8831 21 8 3.6091 -1.8702 5.2928 22 8 3.7234 0.2984 5.7498 23 6 3.5856 -0.0250 7.1586 24 6 4.7435 -0.9271 7.5904 25 6 2.2594 -0.7530 7.3874 26 6 3.6080 1.2640 7.9825 27 6 3.8654 -1.5322 2.5060 28 8 2.0086 -0.6774 1.1673 29 1 -0.3634 0.5118 0.8912 30 1 -0.3634 0.5158 -0.8888 31 1 -0.3634 -1.0277 -0.0023 32 1 1.6423 -1.7397 -1.2673 33 1 1.7166 -0.1893 -2.1388 34 1 5.1211 -1.5286 -2.3351 35 1 3.5762 -2.3807 -2.5707 36 1 4.1683 -0.4237 -6.8053 37 1 5.0672 -1.4654 0.0219 38 1 3.6308 -2.5158 -0.0358 39 1 3.5069 1.3607 0.8728 40 1 5.1111 0.5859 1.0742 41 1 3.1247 1.4407 3.1693 42 1 4.9012 1.3028 3.2866 43 1 4.7270 -1.8454 7.0034 44 1 5.6883 -0.4086 7.4274 45 1 4.6398 -1.1703 8.6479 46 1 2.1561 -0.9961 8.4449 47 1 1.4345 -0.1103 7.0798 48 1 2.2429 -1.6712 6.8005 49 1 4.5526 1.7830 7.8196 50 1 2.7828 1.9062 7.6748 51 1 3.5050 1.0212 9.0400 52 1 4.8718 -1.9332 2.3861 53 1 3.1633 -2.3314 2.7436 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 ZINC000880221844.mol2 53 1SCF WAS SPECIFIED, SO GEOMETRY OPTIMIZATION WAS NOT USED SCF FIELD WAS ACHIEVED Wall clock time and date at job start Wed Apr 1 2020 00:20:51 Heat of formation + Delta-G solvation = -110.165025 kcal Electronic energy + Delta-G solvation = -31792.831851 eV Core-core repulsion = 27665.092478 eV Total energy + Delta-G solvation = -4127.739373 eV No. of doubly occupied orbitals = 66 Molecular weight (most abundant/longest-lived isotopes) = 340.220 amu Computer time = 1.00 seconds Orbital eigenvalues (eV) -41.03085 -40.14549 -37.79537 -36.34306 -35.37091 -34.24436 -32.66953 -32.45157 -30.53698 -29.41484 -27.84097 -27.65464 -27.59774 -26.29826 -24.45043 -23.05984 -22.18662 -21.43424 -20.57736 -19.68413 -19.02242 -17.99068 -17.14015 -16.72854 -16.52784 -16.15549 -15.60938 -15.26042 -15.04807 -14.84434 -14.49803 -14.28621 -14.04014 -13.96753 -13.52203 -13.35345 -13.10708 -12.89013 -12.74539 -12.60518 -12.43187 -12.36495 -12.23342 -12.12670 -12.10443 -12.08464 -11.70451 -11.63883 -11.49679 -11.45433 -11.19229 -11.07409 -10.85226 -10.75228 -10.66867 -10.54130 -10.14888 -10.02539 -9.89774 -9.42468 -9.13491 -8.80039 -8.64876 -7.57639 -6.12576 -4.10781 1.96403 2.77372 3.23248 3.31497 3.32487 3.61210 3.78408 4.32955 4.35700 4.47879 4.51851 4.52680 4.66499 4.78106 4.87122 5.13238 5.19646 5.20922 5.21904 5.27851 5.28955 5.37632 5.41370 5.43449 5.48931 5.51583 5.56541 5.58695 5.62324 5.66363 5.71061 5.72677 5.88037 5.92006 5.95462 5.98691 6.05835 6.08560 6.14583 6.18877 6.27022 6.33619 6.36949 6.50430 6.72847 6.83280 7.39645 7.42794 7.48610 7.71614 7.86038 8.25313 8.54775 8.94370 9.48607 11.00224 Molecular weight = 340.22amu Principal moments of inertia in cm(-1) A = 0.022270 B = 0.002514 C = 0.002406 Principal moments of inertia in units of 10**(-40)*gram-cm**2 A = 1256.994910 B =11136.639305 C =11632.791397 Net atomic charges, atomic populations, and dipole contributions calculated with CM2 Atom NO. Type Charge No. of electrons 1 C -0.141 4.141 2 C 0.086 3.914 3 H 0.068 0.932 4 C 0.032 3.968 5 N -0.520 5.520 6 C 0.163 3.837 7 C -0.511 4.511 8 H 0.081 0.919 9 C 0.070 3.930 10 N -0.898 5.898 11 Si 0.875 3.125 12 H -0.278 1.278 13 H -0.285 1.285 14 H -0.267 1.267 15 C 0.034 3.966 16 C 0.114 3.886 17 C -0.144 4.144 18 C 0.109 3.891 19 N -0.598 5.598 20 C 0.670 3.330 21 O -0.571 6.571 22 O -0.361 6.361 23 C 0.135 3.865 24 C -0.183 4.183 25 C -0.183 4.183 26 C -0.140 4.140 27 C 0.142 3.858 28 O -0.367 6.367 29 H 0.059 0.941 30 H 0.074 0.926 31 H 0.060 0.940 32 H 0.031 0.969 33 H 0.100 0.900 34 H 0.031 0.969 35 H -0.014 1.014 36 H 0.264 0.736 37 H 0.084 0.916 38 H 0.036 0.964 39 H 0.095 0.905 40 H 0.090 0.910 41 H 0.069 0.931 42 H 0.069 0.931 43 H 0.088 0.912 44 H 0.059 0.941 45 H 0.069 0.931 46 H 0.069 0.931 47 H 0.059 0.941 48 H 0.088 0.912 49 H 0.066 0.934 50 H 0.066 0.934 51 H 0.079 0.921 52 H 0.088 0.912 53 H 0.090 0.910 For the charges calculated by CM2: Dipole moment (debyes) X Y Z Total from point charges -1.338 1.787 24.599 24.700 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.199 4.199 2 C 0.029 3.971 3 H 0.086 0.914 4 C -0.096 4.096 5 N -0.244 5.244 6 C 0.037 3.963 7 C -0.549 4.549 8 H 0.099 0.901 9 C -0.132 4.132 10 N -0.537 5.537 11 Si 0.702 3.298 12 H -0.203 1.203 13 H -0.211 1.211 14 H -0.192 1.192 15 C -0.096 4.096 16 C 0.073 3.927 17 C -0.182 4.182 18 C -0.014 4.014 19 N -0.339 5.339 20 C 0.425 3.575 21 O -0.461 6.461 22 O -0.276 6.276 23 C 0.098 3.902 24 C -0.240 4.240 25 C -0.240 4.240 26 C -0.197 4.197 27 C 0.020 3.980 28 O -0.286 6.286 29 H 0.078 0.922 30 H 0.093 0.907 31 H 0.079 0.921 32 H 0.049 0.951 33 H 0.118 0.882 34 H 0.049 0.951 35 H 0.003 0.997 36 H 0.074 0.926 37 H 0.102 0.898 38 H 0.055 0.945 39 H 0.113 0.887 40 H 0.109 0.891 41 H 0.088 0.912 42 H 0.088 0.912 43 H 0.106 0.894 44 H 0.078 0.922 45 H 0.088 0.912 46 H 0.087 0.913 47 H 0.078 0.922 48 H 0.107 0.893 49 H 0.085 0.915 50 H 0.085 0.915 51 H 0.098 0.902 52 H 0.106 0.894 53 H 0.108 0.892 Dipole moment (debyes) X Y Z Total from point charges -1.302 1.878 24.848 24.952 hybrid contribution 0.384 -1.589 -0.546 1.723 sum -0.918 0.289 24.302 24.321 Atomic orbital electron populations 1.21714 0.93140 1.02717 1.02282 1.22333 0.94331 0.95987 0.84433 0.91426 1.22708 0.88833 0.99465 0.98601 1.62292 1.07857 1.52376 1.01898 1.20207 0.95156 0.93219 0.87758 1.20969 1.37359 1.07441 0.89175 0.90096 1.24889 0.95634 0.96792 0.95912 1.69880 1.46972 1.37656 0.99161 0.86872 0.80382 0.84433 0.78095 1.20339 1.21099 1.19183 1.22451 1.00661 0.95822 0.90654 1.21823 0.82517 0.94816 0.93591 1.23064 1.03111 0.98524 0.93472 1.22096 1.01690 0.84182 0.93387 1.47101 1.74861 1.09650 1.02281 1.17419 0.76323 0.82822 0.80926 1.90890 1.56148 1.21900 1.77158 1.86312 1.86164 1.36650 1.18469 1.22057 0.96328 0.96108 0.75661 1.22516 0.97240 1.01690 1.02569 1.22513 0.95708 1.02725 1.03094 1.21850 1.03383 0.95656 0.98815 1.22096 1.02006 0.90429 0.83435 1.87789 1.18994 1.78562 1.43282 0.92191 0.90677 0.92101 0.95052 0.88213 0.95121 0.99674 0.92574 0.89774 0.94544 0.88663 0.89122 0.91212 0.91248 0.89355 0.92221 0.91235 0.91254 0.92164 0.89346 0.91499 0.91457 0.90210 0.89423 0.89160 Number of geometries 1 Number of calculations of the screened coulomb radii 1 The total number of SCF iterations 12. 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.14 -1.79 9.57 37.16 0.36 -1.43 16 2 C 0.09 1.30 3.20 -26.68 -0.09 1.22 16 3 H 0.07 1.01 6.49 -51.93 -0.34 0.67 16 4 C 0.03 0.59 4.72 -3.70 -0.02 0.57 16 5 N -0.52 -10.08 4.59 -236.20 -1.08 -11.16 16 6 C 0.16 3.74 3.79 -4.98 -0.02 3.72 16 7 C -0.51 -13.83 9.03 -34.44 -0.31 -14.15 16 8 H 0.08 2.37 7.36 -52.49 -0.39 1.99 16 9 C 0.07 1.95 5.47 -17.82 -0.10 1.85 16 10 N -0.90 -26.68 13.96 55.43 0.77 -25.90 16 11 Si 0.87 20.89 27.47 -169.99 -4.67 16.22 16 12 H -0.28 -6.57 7.11 56.52 0.40 -6.17 16 13 H -0.29 -6.78 7.11 56.52 0.40 -6.38 16 14 H -0.27 -6.11 6.54 56.52 0.37 -5.74 16 15 C 0.03 0.53 4.29 -4.10 -0.02 0.52 16 16 C 0.11 1.60 0.84 -89.61 -0.08 1.52 16 17 C -0.14 -1.79 5.81 -25.24 -0.15 -1.94 16 18 C 0.11 1.26 6.90 -2.96 -0.02 1.24 16 19 N -0.60 -7.96 3.33 -176.40 -0.59 -8.55 16 20 C 0.67 9.51 7.93 54.06 0.43 9.94 16 21 O -0.57 -9.05 12.15 -19.28 -0.23 -9.28 16 22 O -0.36 -4.54 9.49 -40.29 -0.38 -4.92 16 23 C 0.14 1.35 1.13 -90.62 -0.10 1.25 16 24 C -0.18 -1.74 8.37 37.16 0.31 -1.43 16 25 C -0.18 -1.81 8.37 37.16 0.31 -1.50 16 26 C -0.14 -1.00 8.85 37.16 0.33 -0.67 16 27 C 0.14 1.86 6.09 -2.75 -0.02 1.85 16 28 O -0.37 -5.68 9.99 -35.23 -0.35 -6.03 16 29 H 0.06 0.70 8.14 -51.93 -0.42 0.27 16 30 H 0.07 0.88 8.14 -51.93 -0.42 0.46 16 31 H 0.06 0.78 8.14 -51.93 -0.42 0.35 16 32 H 0.03 0.58 8.14 -51.93 -0.42 0.16 16 33 H 0.10 2.00 6.79 -51.93 -0.35 1.64 16 34 H 0.03 0.71 7.13 -51.93 -0.37 0.34 16 35 H -0.01 -0.34 7.55 -51.92 -0.39 -0.73 16 36 H 0.26 7.38 8.31 -40.82 -0.34 7.04 16 37 H 0.08 1.24 7.98 -51.93 -0.41 0.83 16 38 H 0.04 0.59 8.14 -51.93 -0.42 0.17 16 39 H 0.09 1.24 6.31 -51.93 -0.33 0.92 16 40 H 0.09 1.09 7.91 -51.93 -0.41 0.68 16 41 H 0.07 0.81 8.14 -51.93 -0.42 0.38 16 42 H 0.07 0.71 8.14 -51.93 -0.42 0.29 16 43 H 0.09 1.06 5.88 -51.93 -0.31 0.75 16 44 H 0.06 0.54 8.14 -51.93 -0.42 0.12 16 45 H 0.07 0.51 8.14 -51.93 -0.42 0.09 16 46 H 0.07 0.53 8.14 -51.93 -0.42 0.11 16 47 H 0.06 0.58 8.14 -51.93 -0.42 0.16 16 48 H 0.09 1.10 5.88 -51.93 -0.31 0.80 16 49 H 0.07 0.46 8.14 -51.93 -0.42 0.04 16 50 H 0.07 0.47 8.14 -51.93 -0.42 0.05 16 51 H 0.08 0.41 8.14 -51.93 -0.42 -0.01 16 52 H 0.09 1.02 8.01 -51.93 -0.42 0.60 16 53 H 0.09 1.27 8.05 -51.93 -0.42 0.86 16 LS Contribution 403.73 15.07 6.08 6.08 Total: -1.00 -31.10 403.73 -9.15 -40.25 By element: Atomic # 1 Polarization: 10.26 SS G_CDS: -9.52 Total: 0.74 kcal Atomic # 6 Polarization: 1.73 SS G_CDS: 0.83 Total: 2.56 kcal Atomic # 7 Polarization: -44.72 SS G_CDS: -0.90 Total: -45.61 kcal Atomic # 8 Polarization: -19.26 SS G_CDS: -0.97 Total: -20.23 kcal Atomic # 14 Polarization: 20.89 SS G_CDS: -4.67 Total: 16.22 kcal Total LS contribution 6.08 Total: 6.08 kcal Total: -31.10 -9.15 -40.25 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. ZINC000880221844.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 -69.917 kcal (2) G-P(sol) polarization free energy of solvation -31.102 kcal (3) G-ENP(sol) elect.-nuc.-pol. free energy of system -101.019 kcal (4) G-CDS(sol) cavity-dispersion-solvent structure free energy -9.146 kcal (5) G-P-CDS(sol) = G-P(sol) + G-CDS(sol) = (2) + (4) -40.248 kcal (6) G-S(sol) free energy of system = (1) + (5) -110.165 kcal 1SCF WAS SPECIFIED, SO GEOMETRY OPTIMIZATION WAS NOT USED Total computer time = 1.00 seconds