• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    Simple Ligand Modifications to Modulate the Activity of RutheniumCatalysts for CO2 Hydrogenation: Trans Influence of Boryl Ligands and Nature of Ru―H Bond

    2018-10-19 08:00:44LIUTianLIJunLIUWeijiaZHUYudanLUXiaohua
    物理化學(xué)學(xué)報(bào) 2018年10期

    LIU Tian , LI Jun ,*, LIU Weijia , ZHU Yudan ,*, LU Xiaohua

    1 College of Chemical Engineering, State Key Laboratory of Materials-Oriented Chemical Engineering, Nanjing Tech University,Nanjing 210009, Jiangsu Province, P. R. China.

    2 Nanjing Boiler and Pressure Vessel Inspection Institute, Nanjing 210019, Jiangsu Province, P. R. China.

    Abstract: The development of efficient catalysts for the hydrogenation of CO2 to formic acid (FA) or formate has attracted significant interest as it can address the increasingly severe energy crisis and environmental problems. One of the most efficient methods to transform CO2 to FA is catalytic homogeneous hydrogenation using noble metal catalysts based on Ir, Ru, and Rh. In our previous work, we demonstrated that the activity of CO2 hydrogenation via direct addition of hydride to CO2 on Ir(III) and Ru(II) complexes was determined by the nature of the metal-hydride bond. These complexes could react with the highly stable CO2 molecule because they contain the same distinct metal-hydride bond formed from the mixing of the sd2 hybrid orbital of metal with the 1s orbital of H, and evidently, this property can be influenced by the trans ligand. Since boryl ligands exhibit a strong trans influence, we proposed that introducing such ligands may enhance the activity of the Ru―H bond by weakening it as a result of the trans influence. In this work, we designed two potential catalysts, namely, Ru-PNP-HBcat and Ru-PNPHBpin, which were based on the Ru(PNP)(CO)H2 (PNP = 2,6-bis(dialkylphosphinomethyl)pyridine) complex, and computationally investigated their reactivity toward CO2 hydrogenation. Bcat and Bpin (cat = catecholate, pin = pinacolate)are among the most popular boryl ligands in transition metal boryl complexes and have been widely applied in catalytic reactions. Our optimization results revealed that the complexes modified by boryl ligands possessed a longer Ru―H bond.Similarly, natural bond orbital (NBO) charge analysis indicated that the nucleophilic character of the hydride in Ru-PNPHBcat and Ru-PNP-HBpin was higher as compared to that in Ru-PNP-H2. NBO analysis of the nature of Ru―H bond indicated that these complexes also followed the law of the bonding of Ru―H bond proved in the previous works (Bull.Chem. Soc. Jpn. 2011, 84 (10), 1039; Bull. Chem. Soc. Jpn. 2016, 89 (8), 905), and the d orbital contribution of the Ru atom in Ru-PNP-HBcat and Ru-PNP-HBpin was smaller than that in Ru-PNP-H2. Consequently, the Ru-PNP-HBcat and Ru-PNP-HBpin complexes were more active than Ru-PNP-H2 for the direct hydride addition to CO2 because of the lower activation energy barrier, i.e., from 29.3 kJ·mol?1 down to 24.7 and 23.4 kJ·mol?1, respectively. In order to further verify our proposed catalyst-design strategy for CO2 hydrogenation, the free energy barriers of the complete pathway for the hydrogenation of CO2 to formate catalyzed by complexes Ru-PNP-H2, Ru-PNP-HBcat, and Ru-PNP-HBpin were calculated to be 76.2, 67.8, and 54.4 kJ·mol?1, respectively, indicating the highest activity of Ru-PNP-HBpin. Thus, the reactivity of Ru catalysts for CO2 hydrogenation could be tailored by the strong trans influence of the boryl ligands and the nature of the Ru―H bond.

    Key Words: CO2 hydrogenation; Ru complex; Boryl ligand; Trans influence; Ru―H bond

    1 Introduction

    The carbon dioxide emissions have increased owing to rising energy consumption. The catalytic hydrogenation of CO2into fuels and useful chemicals, especially for formic acid (FA) or formate, holds potential to address this problem. On one hand,hydrogenation of CO2to produce FA1–7could mitigate the CO2in atmosphere. On the other hand, FA can serve as the most promising hydrogen storage media and can be directly used for formic acid fuel cell1,8. In addition, FA and its salts are widely used in the textile, leather, dye industry and other industrial areas9.It is well known that CO2 is an attractive C1 feedstock with the advantages of abundance, non-flammability and nontoxicity, but its inherent stability still poses a big challenge to chemists10,which determines the key issue in the process of CO2transformation is activation of CO2.

    Homogeneous hydrogenation of CO2to formic acid by noble metal catalysts has become one of the most efficient pathways11–16.In 2009, Tanaka et al.14reported an impressive pincer-supported iridium trihydride catalyst Ir(III)-PNP for CO2hydrogenation to formate, the turnover frequencies (TOFs) and turnover number(TON) achieved 150000 h?1and 3500000 at 200 °C and 5000000 Pa. Compared to iridium catalyst, ruthenium catalysts are more common in CO2 hydrogenation. More importantly, the ruthenium element is much less expensive. Currently, a significant breakthrough was achieved by Filonenko et al. by combination of complex Ru(PNP)(CO)H2and strong DBU base15,17. It allows TOFs as high as 1892000 h?1at the lower temperature and pressure (132 °C and 4000000 Pa). Previously,we found that they are so active to catalyze highly stable CO2because they contain same distinct metal-hydride bond formed from the mixing of the sd2hybrid orbital of metal with the 1s orbital of H18. In addition, this unique nature of metal-hydride bond also plays vital role in CO2activation in other ruthenium catalysts, such as trans-Ru(dmpe)2H2, cis-Ru(dmpe)2H2, (dmpe =Me2PCH2CH2PMe2), and cis-Ru(PMe3)3(H2O)H219. A distinct feature of them is that the hydrido complex is capable of transferring a hydride to CO2, which makes the CO2 activation is possible. Even though these studies have achieved extraordinary efficiency, the harsh reaction conditions have to be carried out. Thus, the design and development of high-efficiency catalysts that can operate under milder conditions is highly desirable.

    Developing a highly efficient catalyst need to accelerate the process of hydride transferring, since this process is postulated as key issue in CO2activation. The comparison of trans-Ru(dmpe)2H2and cis-Ru(dmpe)2H2may give us some enlightenment. The similarity between the two is that they are in the same ligand environment. Differently, the trans ligand of hydride of trans-Ru(dmpe)2H2 is another hydride, the counterpart of cis-Ru(dmpe)2H2is phosphine ligand. The difference of Ru―H bond properties and activation barriers of two isomers might relate to the ligand trans to the hydride. Our earlier study also showed that the ligand trans to hydride is an important factor to influence the bond distance and nature of Ru―H bond of Ru(PNP)(CO)H219. Such a phenomenon can be attributed to trans influence of ligand opposite to hydride, which means that a ligand can influence some properties of its trans ligand, such as bond distance, the vibrational frequency, force constant and a battery of other parameters20. Therefore, the activity of the hydride on a Ru(II) complex might be adjusted by the change of the ligand trans to the hydride. Schmeier and Langer’s group21,22also consciously controlled this factor to optimize catalytic performance and designed active pincer complexes for CO2hydrogenation containing two hydrides trans to each other. However, they have already thought hydride is the most active trans ligand and have not explored the nature of metal-hydride bond. These studies inspire us to explore stronger trans influence than hydride ligand trans to reacting hydride to weaken Ru―H bond and apply it to CO2hydrogenation reaction.

    Boryl ligands are known to its strong trans influence in transition-metal boryl complexes23,24. Besides, boryl ligands have shown to alter the electronic property of transition metals markedly25. Lin and Zhu have investigated the trans influence of boryl ligands via density functional theory (DFT) calculations on square-planar platinum(II) complexes, and presented that trans influence of boryl ligands are stronger than hydride group20.In addition, bis-phosphino-boryl (PBP) pincer transition-metal complexes were synthesized26, and were applied in the field of catalysis. For instance, Peters et al.27have synthesized cobalt(I)-N2complex supported by PBP ligand and explored its catalytic performance on catalytic olefin hydrogenation as well as amineborane dehydrogenation. He also presented the synthesis of boryl-Ni complex (tBuPBP)NiH, which reacts with CO2due to trans influence engendered by boryl ligand to give formate complex28.

    In continuation of our work on ruthenium catalysts and combined with strong trans influence of boryl ligands, we would like to design more efficient catalysts based on Ru(II) complex Ru(PNP)(CO)H2for CO2hydrogenation. Designing more efficient ruthenium catalysts, in other words in this system,designing more efficient hydride that attack to CO2. Given strong trans influence of boryl ligands, we would like to take an easy coordinating form between ruthenium and boryl ligands into consideration, and position the boryl ligands trans to hydride in ruthenium complexes so as to exert its strong trans influence. Therefore, it will be interesting to study the role of boryl ligand for changing character of its trans hydride, which in turn might be helpful for CO2hydrogenation reaction. To this end, two boryl ligands are employed, and the corresponding complexes are presented in Fig. 1 through substituting one of hydrides in Ru-PNP-H2with Bcat and Bpin ligands, named Ru-PNP-HBcat and Ru-PNP-HBpin. Bcat and Bpin ligands are among the most popular boryl ligands in transition metal boryl complexes and widely applied in catalysis reaction29,30. In this paper, we have computationally examined the activities of two complexes for the hydrogenation of CO2using DFT, and compared the reactivity with their parent complex Ru-PNP-H2.Since the key step in hydrogenation of CO2 is insertion of CO2 into the active metal-hydride bond31,32, so as to Ru-PNP-H2.Furthermore, the step of formato ligand isomerization is also vital to CO2transformation since it may be more difficult than hydride transfer18, and the barriers of CO2transformation to formate were also calculated in this work.

    2 Computational details

    DFT calculations were implemented using the Gaussian 09 suite of programs33. The geometry optimization has employed the B3LYP exchange-correlation functional34. The choice of the B3LYP functional is based on its successful application in describing the ruthenium complexes13,35,36. The Stuttgart-Dresden basis set-relativistic effective core potential (RECP)combination, which was supplemented with two sets of f functions and a set of g functions37, was used for Ru atom. The Dunning cc-pVDZ basis set was used for P, O, N, C, B and H atom38. Same basis sets were employed for these atoms in our earlier works18,19,39and others40–42.

    In order to maintain consistency with the previous work19, all geometry optimizations were carried out in vacuum. To verify the validity of geometry optimization of complexes, the comparison with related complexes is necessary, and found that they are consistent with experimentally reported Ru complexes43–45.For the sake of identifying all the stationary points as minima(zero imaginary frequency) or transition states (one imaginary frequency) and obtaining ZPE corrections, frequency calculation were carried out on all the optimized structures at the same level of theory used for the optimization step. Intrinsic reaction coordinate (IRC) calculations were adopted to confirm the connection of the transition states and two relevant minima. The properties and WBI of Ru―H bond were characterized by natural bond orbital (NBO) analysis46, and carried out on the structures of some intercepted stationary points. For all ruthenium complexes, the tert-butyl groups of P atom were replaced by methyl groups, this simplification also carried out in others paper to decrease computational cost21. Gibbs free energies of all species were computed at 298.15 K and 101300 Pa.

    3 Results and discussion

    Fig. 1 Optimized geometries and schematic drawings of the structures of Ru(II) complexes generated from Ru-PNP-H2 with modified ligands.Bond distances are in nanometer (nm). H atoms on methyl group are omitted for clarity. Green balls: Ru; violet: P; red: O; pink: B; blue: N; grey:C; white: H; orange: methyl group, color online.

    In order to present more straightforward trans influence of boryl ligand, a series of parameters of hydride ligand are considered as the metric of trans influence and catalytic efficiency. The hydride ligand is in the same position as boryl ligands and commonly used in pincer complex for homogeneous CO2hydrogenation. In this work, we would explore boryl ligands how to exert its strong trans influence in key structural parameters, the nature of Ru―H bond and modulate the activity toward CO2hydrogenation compared to hydride ligand at an atomic level. These Ru complexes (Fig. 1) are modeled by changing the ligand trans to reacting hydride. Ru-PNP-H2 have two hydrides trans to each other, one of the axial hydride ligands of complex Ru-PNP-H2is substituted by a Bcat and Bpin ligand to model complex Ru-PNP-HBcat and Ru-PNP-HBpin. The reserved hydride ligand is considered as reacting hydride to form ruthenium hydrido complex, and what discussed are based on these complexes.

    3.1 Atomic and electronic structure of Ru(II) complex

    The optimized geometries and schematic drawings of the structures of Ru-PNP-H2, Ru-PNP-HBcat and Ru-PNP-HBpin are shown in Fig. 1. Some key structural parameters are listed in Table 1. First of all, distance of metal-hydride bond play an important role in hydride transferring process. Ru―H bond distance of Ru-PNP-H2 is 0.1697 nm, and the Ru―H bonddistances of Ru-PNP-HBcat and Ru-PNP-HBpin are 0.1714 and 0.1716 nm, respectively. The distances of Ru―H bonds on the two complexes Ru-PNP-HBcat and Ru-PNP-HBpin are larger than that of Ru-PNP-H2, illustrate that the boryl ligands exert strong trans influence to stretch the Ru―H bond. Furthermore,the bond length of Ru―B are 0.2140 and 0.2170 nm in Ru-PNPHBcat and Ru-PNP-HBpin, indicate that the bond length of Ru―B is influenced by its substituents in the B atom. The tendency of distance of Ru―B bonds is similar to Pt―B bonds of square-planar platinum(II) boryl complexes47. The angle of three atoms in trans position of three complexes are 175.7°,176.0°, 177.7°, respectively. This trend implies that boryl ligands are more prone to locate trans position of hydride, and the magnitude of the trend is in the order of Bpin > Bcat > H. Finally,the relative position of boryl ligands is appealing. The dihedral angle of P1-P2-O1-O2 in Ru-PNP-HBcat and Ru-PNP-HBpin are 62.0° and 35.2°, and are not ideal octahedral geometry. The smaller dihedral angle in Ru-PNP-HBpin might result from the bulky nature of Bpin ligand. The NBO charge analysis is done to understand the charge variation when varied the trans ligand of hydride. Table 2 listed NBO charges of various ligands, the stretching frequency and WBI of Ru―H bond. When we changed the trans ligand of hydride, the PNP and CO ligands areacting as donators. In addition, introducing boryl ligands also increases electron density at the metal center, and the electron density of metal center is one of influencing factor related to the barrier for hydride transfer48. More importantly, the NBO charge on the reacting hydride becomes more negative as we change boryl ligands with different trans influence. The nucleophilicity of the hydride can be estimated by NBO charge,the more negative charge, and the more nucleophilic. From this perspective, the natural charge on hydride of Ru-PNP-H2 is?0.245e, the natural charge on hydride are ?0.270e and ?0.276e for Ru-PNP-HBcat and Ru-PNP-HBpin, respectively. Boryl ligands make the natural charge on hydride more negative,indicates that the character of nucleophilicity of hydride are increased in Ru-PNP-HBcat and Ru-PNP-HBpin, we attribute a strongly donating trans ligand of Bcat and Bpin20. Similar to bond distance of Ru―H, stretching frequencies of Ru―H bonds in complexes modified by Bcat and Bpin (1658, 1644 cm?1) are significantly smaller than that of Ru-PNP-H2(1853 cm?1),demonstrate that a weaker Ru―H bond in Ru-PNP-HBcat and Ru-PNP-HBpin than that of Ru-PNP-H2. Moreover, Wiberg bond indices (WBI) also support the relatively weak interaction between hydride and the central Ru atom in Ru-PNP-HBcat and Ru-PNP-HBpin. The WBI value for the Ru―H bond of 0.8620 in Ru-PNP-H2is higher than those of Ru-PNP-HBcat and Ru-PNP-HBpin with 0.8372, 0.8370, respectively. In the present case, the bond distance, NBO charge on hydride, stretching frequency and WBI of Ru―H bond all indicate that the hydride in Ru-PNP-HBcat and Ru-PNP-HBpin should be more favorable for hydrogenation of CO2, and more detailed theoretical study of CO2hydrogenation mediated by them is required to confirm our suppose.

    Table 1 Structural comparison among three Ru(II) complexes

    Table 2 Summed up NBO charges on the ligands and metal, stretching frequency and WBI of Ru―H bond of three complexes.

    3.2 The nature of Ru―H bond of Ru(II) complex

    A longer, weaker Ru―H bond, higher electron density of Ru,and more nucleophilic hydride of Ru-PNP-HBcat and Ru-PNPHBpin indicate that the hydride of them are more prone to nucleophilic attack to CO2than Ru-PNP-H2. We turn to calculate the activation barrier for the hydrogenation of CO2via direct addition of hydride to CO2. And the nature of Ru―H bond was explored, since it is the nature of the Ru―H bond that determines the activity of hydride addition to CO2. New properties of Ru―H bonds were investigated based on NBO analysis, and the bonding parameters of these Ru―H bonds were collected in Table 3. We remark that the property of Ru―H bond of Ru-PNP-H2is slightly different from our earlier paper because tert-butyl groups of P atoms were replaced by methylgroups in this work. Compared to hydride ligand, trans influence of boryl ligands on the nature of Ru―H bond is consistent to their influence of Ru―H bond distance, stretching frequency,and natural charge of hydride. The properties of Ru–H bonds in Ru-PNP-HBcat, Ru-PNP-HBpin are formed via the mixing of the sd1.83hybrid orbital of Ru with the 1s orbital of H and the sd1.80hybrid orbital of Ru with the 1s orbital of H, respectively.The d orbital contribution of Ru atom in Ru-PNP-HBcat and Ru-PNP-HBpin are less than that of Ru-PNP-H2 formed from sd2.11hybrid orbital of Ru with the 1s orbital of H. For activity of CO2 activation, Ru-PNP-HBcat and Ru-PNP-HBpin make the activation energies for hydride addition to CO2from 29.3 kJ·mol?1drop to 24.7 and 23.4 kJ·mol?1, respectively, and the Gibbs free energy change of direct hydride transfer (1 → 2) to the C atoms of the CO2 are calculated from 29.7 kJ·mol?1decline to 22.6 and 18.8 kJ·mol?1, respectively. Thus, if the trans ligand alters the property of Ru―H bond with less d orbital components than hydride, its activation energy and Gibbs free energy change for hydride direct addition to CO2would decrease to some extent. The longer Ru―H bond, the less d orbital component, and the lower activation barrier, as shown in Table 3, which is consistent with our previous results19. The complexes Ru-PNP-HBcat and Ru-PNP-HBpin would be more active than complex Ru-PNP-H2for the step of direct hydride addition to CO2due to lower activation energy and larger Ru―H bond distance. Combined with the properties of Ru―H bond and activation barriers, the enhancement of activity of direct hydride addition to CO2 for Ru-PNP-HBcat and Ru-PNP-HBpin might originates from the less d orbital component of Ru atom, which leads to the ease of transferring hydride to CO2. This conjecture is supported by the WBI value of Ru―H bonds, the lower d orbital components are involved in the formation of the Ru―H bond, corresponding to the relatively weak interaction between hydride and Ru.

    Table 3 Bond distance, bonding parameter of Ru―H bond, activation energy for hydride (E a) addition to CO2 on three complexes.

    3.3 Hydrogenation process of CO2 to formate on Ru(II) complex

    Besides the step of hydride direct addition to CO2, we also have investigated the detailed whole pathway for the hydrogenation of CO2 to formate on the new systems, Ru-PNPHBcat and Ru-PNP-HBpin, and compared the activity with their parent complex Ru-PNP-H2. The whole pathway for the hydrogenation of CO2to formate are composed of below several steps: 1) CO2interacts with the Ru(II) complex; 2) hydride transfer from the Ru hydrido complex to CO2; 3) the rotation of the formato ligand, Scheme 1. In earlier work, our attention was focused on the activation energy barrier of hydride addition to CO2due to this process is determined by the nature of Ru―H bond. In the present case, our work presented that not only the nature of Ru―H bond correlates with the hydride transferring process, but also the trans ligand of hydride relates to catalytic activity in the whole pathway for the hydrogenation of CO2 to formate. Furthermore, some study49demonstrated that the process of CO2insertion into Ru―H bond via direct hydride abstraction pathway is a sequentially uphill step. Therefore, the whole process of hydrogenation of CO2 to formate should be paid attention and the overall barriers are used to measure the catalytic efficiency of these complexes.

    Scheme 1 The hydrogenation of CO2 to formate on three complexes (X = H, Bcat, Bpin).

    The calculated free energy profile is presented in Fig. 2 and the optimized structures of key intermediates and transition states are shown in Fig. 3. All catalysts and CO2were considered as the starting point of energy profiles. Firstly, complexes Ru-PNP-H2 and Ru-PNP-HBcat interact with CO2 and form an octahedral precursor complex 1, 1-Bcat. This process is an endergonic step, which are slightly uphill by 14.7 and 18.0 kJ·mol?1, respectively. In the precursor complex 1 and 1-Bcat,the Ru―H bonds distances maintain to 0.1695 and 0.1714 nm,and the distance of H and CO2 are far to 0.2629 and 0.2690 nm,respectively. This behavior resembles to that found in the case of Ir(III) trihydride complex that did present a CO2directly attacked to the nucleophilic hydride18. Similarly, the bond distances of Ru―H bond and Ru―B bond positioned trans to reacting hydride also slightly shorten to 0.1694 and 0.2133 nm,and this tendency is similar to other study50. Furthermore, CO2 is weakly activated because its bond angle changes from linear to 176.5°, 176.8°, respectively. The angle of CO2, distance between C and H atoms, and distance between Ru and H atoms are listed in Table 4 as the hydride addition to CO2proceeds.Next step is hydride migrate to CO2to form H-bound formate complexes (2, 2-Bcat), the formation of H-bound formate complex requires a transition state TS1-2. This step is key to whole process, the calculated free energies barriers are 29.3 and 24.7 kJ·mol?1for Ru-PNP-H2and Ru-PNP-HBcat, respectively.Our calculation results indicate that Ru-PNP-HBcat catalyst has the lower activation barrier for hydride transfer to CO2. In formate complexes (2, 2-Bcat), the distance of Ru―H bond becomes longer are accompany with the approach of C atom and H atom. The distance between C atom and H atom are 0.1280,0.1266 nm, and the distances of the Ru―H bond are 0.1857 and 0.1901 nm, respectively, indicate that the breaking of Ru―H bond and the formation of C―H bond. The CO2molecules also change from linear to bent fashion and reach to 137.3° in 2-Bcat,which is greater deformation than counterpart complex 2.Subsequently, we investigated the next step of catalytic cycle,namely, the rotation of the formato ligand to form the O-bound formate complex 3 and 3-Bcat owing to H-bound formate complex is not stable. This step is more difficult than forward step, and the step of formation of formate is a successive uphill process. This step is exergonic by 17.2 kJ·mol–1with an activation energy barrier of 31.8 kJ·mol–1for Ru-PNP-H2. For Ru-PNPHBcat, it is more exergonic by 1.7 kJ·mol?1with a less activation energy of 27.2 kJ·mol–1. The series of pivotal thermodynamic and kinetic data (Gibbs free energy change, activation barrier)illustrate that the model complex Ru-PNP-HBcat is more efficient than the extraordinarily active catalyst Ru-PNP-H2. The optimized transition state structures (TS2-3, TS2-3-Bcat) of rotation of the formato ligand are shown in Fig. S1, the broken Ru―H bond continuously stretch to 0.2413 and 0.2501 nm,respectively. In addition, the two C―H bonds both shorten to as much as 0.1160 nm. The ultimate O-bound formate complexes 3, 3-Bcat are stable, the C―H bonds both shorten to 0.1122 nm in formato ligand, and the distance between Ru and O atoms are 0.2230 and 0.2260 nm, respectively.

    Fig. 2 Free energy profiles for the key steps of CO2 hydrogenation by three Ru(II) complexes.Schematic drawings of the structures for the species involved are provided for clarity.

    Fig. 3 Optimized geometries of the species involved in the direct addition of hydride to CO2 on the three Ru(II) complexes.Bond distances are in nanometer (nm) and bond angles are in degree (°). H atoms on methyl group are omitted for clarity. Green balls: Ru; violet: P; red: O;pink: B; blue: N; gray: C; white: H; orange: methyl groups, color online.

    Table 4 Angle of CO2, distance between C and H atoms, distance between Ru and H atoms for hydride addition to CO2 on three complexes.

    For complex Ru-PNP-HBpin, its reaction pathway is similar to the case for Ru-PNP-HBcat. First of all, CO2interacts with hydride of Ru-PNP-HBpin to form association complex 1-Bpin.It is worth noting that this association complex 1-Bpin is less endergonic than complex 1 by 0.4 kJ·mol?1, while in 1-Bcat is more endergonic by 3.3 kJ·mol?1. In the association complex 1-Bpin, the Ru―H bond elongated to 0.1717 nm, and the distance between hydride and C atom of CO2is 0.2674 nm. The angel of CO2 is 176.6°, which is weakly activated comparable to complex 1 and 1-Bcat. Next step is hydride transfer to electron-deficient C atom to form H-bound formate complex (1-Bpin), transition state is found in this process. The calculated activation barrier is 23.4 kJ·mol?1, which have a moderate reduction in comparison to Ru-PNP-H2(29.3 kJ·mol?1). In addition, its Gibbs free energy change also have a moderate decreases by 11.3 and 7.5 kJ·mol?1compared to Ru-PNP-H2 and Ru-PNP-HBcat, respectively. The distance between Ru atom and H atom of formed H-bound formate complex 2-Bpin has a longer with 0.1931 nm compared to complex 2 (0.1857 nm), 2-Bcat (0.1901 nm). Besides, the angel of CO2in H-bound formate complex (2-Bpin) is 136.3°,indicates that CO2is polarized completely due to electrons have been transferred into CO218. For the rotation of the formato ligand of the H-bound formate complex (2-Bpin) to form O-bound formate complex (3-Bpin), this step is exothermic by 24.7 kJ·mol?1and has an activation barrier of 21.4 kJ·mol?1. In transition state structures of TS2-3-Bpin, the distance between Ru and H atoms continuously elongates to 0.2562 nm, which is longer than that of TS2-3 and TS2-3-Bcat (Fig. S1).Furthermore, the distance between Ru and O atoms in the O-bound formate complexes is 0.2281 nm due to strong trans influence of Bpin ligand, which is larger than those of 3 and 3-Bcat. Surprisingly, in this model catalyst, rotation of the formato ligand is easier than hydride transfer to CO2. The overall free energy barrier is 54.4 kJ·mol?1, which leads to decreases of 21.8 kJ·mol?1in the overall barrier compared to Ru-PNP-H2. And the results illustrate that a simple ligand modification imparts changes in the overall activation barrier for the CO2hydrogenation on complex Ru-PNP-HBpin.

    The process of hydrogenation of CO2to formate is just half catalytic cycle of production of FA, and we focus on this process because CO2addition to hydride is the rate-determining step calculated by Filonenko et al.17. They predicted the total activation energy barrier is 67 kJ·mol?1, which is in good agreement with our result. But they ignored the step of formato ligand isomerization, and this step also studied in other ruthenium complexes13,35,51. We can see from the energy profile in Fig. 2 shown above: all the complexes proceed via direct abstraction hydride mechanism for hydrogenation of CO2 to formate. The overall activation barriers for CO2 hydrogenation to formate on complexes Ru-PNP-H2, Ru-PNP-HBcat, and Ru-PNP-HBpin are 76.2, 67.8, and 54.4 kJ·mol?1, respectively. Ru-PNP-HBcat and Ru-PNP-HBpin that contain boryl ligand exerted strong trans influence both have lower activation energy of whole pathway for the hydrogenation of CO2to formate than Ru-PNP-H2. The complex Ru-PNP-HBpin could be a very efficient catalyst for the hydrogenation of CO2under mild condition. The free energy barriers of hydrogenation of CO2are influenced by ligand trans to reacting hydride significantly, and an appropriate adjustment of the functional groups trans to reacting hydride could greatly improve the catalytic activity.

    4 Conclusions

    In conclusion, our current calculations performed by detailed DFT calculations verify our proposed catalyst-design strategy for CO2hydrogenation. Ru-PNP-HBcat and Ru-PNP-HBpin are proposed based on strong trans influence of boryl ligands and examined the activities of two complexes for the hydrogenation of CO2. From perspective of related structural, electronic parameters and the nature of Ru―H bond, the introduction of boryl ligand increase the nucleophilic character of hydride, and make the hydride involved in a Ru―H bond have less d orbital components than that of Ru-PNP-H2. As a consequence, the step of direct addition of hydride to CO2 for Ru-PNP-HBcat and Ru-PNP-HBpin have lower activation barriers compared to Ru-PNP-H2. Moreover, introducing strongly trans-influence boryl ligands trans to reacting hydride is found to speed up the overall process of hydrogenation of CO2to formate. The overall activation barriers for CO2hydrogenation on Ru-PNP-HBcat,Ru-PNP-HBpin are 67.8, 54.4 kJ·mol–1, which are 8.4 and 21.8 kJ·mol–1lower than the reaction catalyzed by Ru-PNP-H2,respectively. Therefore, the reactivity of ruthenium complexes can be fine-tuned by changing trans ligand of active hydride to influence the nature of Ru―H bond. Further design of non-noble metal catalysts based on strong trans influence of boryl ligands and more in depth reasons for their excellent performance in CO2hydrogenation are underway.

    Acknowledgment: We are grateful to the High Performance Computing Center of Nanjing Tech University for supporting the computational resources.

    Supporting Information: available free of charge via the internet at http://www.whxb.pku.edu.cn.

    亚洲欧美精品综合久久99| 最近最新免费中文字幕在线| 精品99又大又爽又粗少妇毛片 | 欧美一区二区亚洲| 亚洲成a人片在线一区二区| 高清在线国产一区| 国产色爽女视频免费观看| 亚洲精品乱码久久久v下载方式| 最新在线观看一区二区三区| 国产成人aa在线观看| 热99re8久久精品国产| 别揉我奶头~嗯~啊~动态视频| 亚洲欧美日韩高清专用| 禁无遮挡网站| 美女xxoo啪啪120秒动态图| 天美传媒精品一区二区| 国国产精品蜜臀av免费| 欧美性猛交黑人性爽| 性色avwww在线观看| 蜜桃久久精品国产亚洲av| 熟妇人妻久久中文字幕3abv| 大型黄色视频在线免费观看| 国产国拍精品亚洲av在线观看| 嫁个100分男人电影在线观看| 如何舔出高潮| 欧美不卡视频在线免费观看| 国产蜜桃级精品一区二区三区| 欧美成人性av电影在线观看| 精品99又大又爽又粗少妇毛片 | 国产精品久久久久久av不卡| 国产大屁股一区二区在线视频| 色综合色国产| 九九久久精品国产亚洲av麻豆| 日日干狠狠操夜夜爽| 麻豆精品久久久久久蜜桃| a级一级毛片免费在线观看| 悠悠久久av| 亚洲18禁久久av| 最后的刺客免费高清国语| 欧洲精品卡2卡3卡4卡5卡区| 亚洲精品一卡2卡三卡4卡5卡| 久久精品久久久久久噜噜老黄 | 日本黄色视频三级网站网址| 人人妻,人人澡人人爽秒播| 成年女人永久免费观看视频| 亚洲成人久久爱视频| 午夜精品一区二区三区免费看| 亚洲成人免费电影在线观看| 国产精品国产三级国产av玫瑰| 午夜福利在线观看吧| 精品久久久久久成人av| 国产高清不卡午夜福利| 久久婷婷人人爽人人干人人爱| 国产精品自产拍在线观看55亚洲| 在线免费十八禁| 观看免费一级毛片| 国产视频一区二区在线看| 黄色女人牲交| 久久九九热精品免费| 精品久久久久久久久久免费视频| 成人毛片a级毛片在线播放| 亚洲美女视频黄频| 久久久久久国产a免费观看| 嫩草影院精品99| 国产精品自产拍在线观看55亚洲| 日韩欧美精品v在线| 99在线人妻在线中文字幕| 精品久久久久久久末码| 日本黄色片子视频| 欧美中文日本在线观看视频| 99热这里只有是精品50| 黄色一级大片看看| 成人av一区二区三区在线看| 精华霜和精华液先用哪个| 日本精品一区二区三区蜜桃| 日韩av在线大香蕉| av在线老鸭窝| 久久热精品热| 亚洲精品国产成人久久av| 无人区码免费观看不卡| 日韩,欧美,国产一区二区三区 | 国产69精品久久久久777片| 欧美人与善性xxx| 国产 一区精品| 一区二区三区四区激情视频 | 国产视频内射| 特级一级黄色大片| 久久久久久久久大av| 午夜老司机福利剧场| 亚洲久久久久久中文字幕| 99久久九九国产精品国产免费| 久久精品夜夜夜夜夜久久蜜豆| 亚洲黑人精品在线| 日本一本二区三区精品| 亚洲av中文av极速乱 | 白带黄色成豆腐渣| 免费不卡的大黄色大毛片视频在线观看 | 在线观看一区二区三区| 18禁裸乳无遮挡免费网站照片| 观看免费一级毛片| 波多野结衣高清无吗| 亚洲狠狠婷婷综合久久图片| 亚洲精品456在线播放app | 禁无遮挡网站| 18+在线观看网站| 神马国产精品三级电影在线观看| 女人十人毛片免费观看3o分钟| 麻豆精品久久久久久蜜桃| 国产精品99久久久久久久久| 久久精品国产亚洲av香蕉五月| 国模一区二区三区四区视频| 免费人成在线观看视频色| 又粗又爽又猛毛片免费看| a级毛片a级免费在线| 可以在线观看的亚洲视频| 中文字幕高清在线视频| 精品久久久久久久久亚洲 | 国产 一区 欧美 日韩| 国产成人一区二区在线| 色噜噜av男人的天堂激情| 欧美成人一区二区免费高清观看| 男女那种视频在线观看| 2021天堂中文幕一二区在线观| 国产黄色小视频在线观看| 别揉我奶头 嗯啊视频| 亚洲中文字幕日韩| or卡值多少钱| 蜜桃久久精品国产亚洲av| 国产精品亚洲一级av第二区| 色综合站精品国产| 欧美日韩亚洲国产一区二区在线观看| 亚洲精品一卡2卡三卡4卡5卡| 香蕉av资源在线| 毛片女人毛片| 在线观看舔阴道视频| 91av网一区二区| 日本 欧美在线| 老女人水多毛片| 亚洲中文字幕一区二区三区有码在线看| 人妻夜夜爽99麻豆av| 亚洲熟妇熟女久久| 国产精品久久久久久久久免| 中文字幕av在线有码专区| 久久久久国产精品人妻aⅴ院| 色哟哟哟哟哟哟| 精华霜和精华液先用哪个| 亚洲最大成人av| 熟女人妻精品中文字幕| 美女被艹到高潮喷水动态| 欧美日韩黄片免| 精品人妻视频免费看| 女生性感内裤真人,穿戴方法视频| 精品人妻熟女av久视频| 日韩一区二区视频免费看| 亚洲专区中文字幕在线| 久久中文看片网| 久久久久国内视频| 国产又黄又爽又无遮挡在线| 夜夜看夜夜爽夜夜摸| 美女高潮喷水抽搐中文字幕| 国产高清视频在线观看网站| 久久久久九九精品影院| 免费看日本二区| 别揉我奶头 嗯啊视频| 亚洲国产精品合色在线| 亚洲 国产 在线| 久久人人爽人人爽人人片va| 亚洲av熟女| 国产美女午夜福利| 真实男女啪啪啪动态图| 看免费成人av毛片| 日本免费一区二区三区高清不卡| 99久久精品热视频| 欧美一区二区亚洲| 熟女人妻精品中文字幕| 国内久久婷婷六月综合欲色啪| av在线亚洲专区| 国产高清视频在线播放一区| 国产男人的电影天堂91| 国产亚洲av嫩草精品影院| 亚洲精华国产精华液的使用体验 | 国产精品综合久久久久久久免费| 亚洲av第一区精品v没综合| 亚洲天堂国产精品一区在线| 久久国产乱子免费精品| 国产成人aa在线观看| 国产精品一区二区三区四区免费观看 | 亚洲自偷自拍三级| 色av中文字幕| 国产高潮美女av| 91在线精品国自产拍蜜月| 亚洲av成人精品一区久久| 国产精品福利在线免费观看| 日韩亚洲欧美综合| 在线观看66精品国产| 久久久久久久午夜电影| 国产大屁股一区二区在线视频| 久久久午夜欧美精品| 精品久久久久久久末码| 国产探花极品一区二区| 久久中文看片网| 国产亚洲欧美98| 中文亚洲av片在线观看爽| 狂野欧美激情性xxxx在线观看| 日韩欧美在线二视频| 成人毛片a级毛片在线播放| av视频在线观看入口| 成人综合一区亚洲| 99久久精品热视频| 老司机午夜福利在线观看视频| 动漫黄色视频在线观看| 夜夜爽天天搞| 国产在线男女| 亚洲av免费在线观看| 久久99热6这里只有精品| 又紧又爽又黄一区二区| 日日干狠狠操夜夜爽| 啦啦啦韩国在线观看视频| 色噜噜av男人的天堂激情| 国产人妻一区二区三区在| 欧美精品啪啪一区二区三区| a级毛片a级免费在线| 又爽又黄a免费视频| 国产av在哪里看| 男人的好看免费观看在线视频| 久久久久免费精品人妻一区二区| 有码 亚洲区| 免费黄网站久久成人精品| 国产乱人视频| 春色校园在线视频观看| 亚洲av五月六月丁香网| 免费看美女性在线毛片视频| 欧美潮喷喷水| 草草在线视频免费看| 午夜精品久久久久久毛片777| netflix在线观看网站| 欧美日韩综合久久久久久 | 国产大屁股一区二区在线视频| 日本熟妇午夜| 亚洲精品亚洲一区二区| 波多野结衣高清作品| 亚洲精品影视一区二区三区av| 中文在线观看免费www的网站| 成人特级黄色片久久久久久久| 久久婷婷人人爽人人干人人爱| 成人午夜高清在线视频| 99久国产av精品| 免费观看在线日韩| 亚洲av日韩精品久久久久久密| 久久99热6这里只有精品| 两性午夜刺激爽爽歪歪视频在线观看| 五月伊人婷婷丁香| 97热精品久久久久久| 久久国产精品人妻蜜桃| 一区二区三区高清视频在线| 久久久精品欧美日韩精品| 成人av在线播放网站| 亚洲三级黄色毛片| 一区二区三区高清视频在线| 制服丝袜大香蕉在线| 如何舔出高潮| av天堂中文字幕网| 欧美精品国产亚洲| 很黄的视频免费| 亚洲,欧美,日韩| 久久亚洲精品不卡| 亚洲va在线va天堂va国产| 亚洲在线自拍视频| 成人av一区二区三区在线看| 91久久精品国产一区二区三区| 国内精品久久久久久久电影| 色5月婷婷丁香| 午夜a级毛片| 久久久久九九精品影院| 免费av不卡在线播放| 一级毛片久久久久久久久女| 大又大粗又爽又黄少妇毛片口| 国产主播在线观看一区二区| 国内精品久久久久久久电影| 久久精品91蜜桃| 岛国在线免费视频观看| 欧美+日韩+精品| 99在线人妻在线中文字幕| 亚洲第一电影网av| 成人国产麻豆网| 成人鲁丝片一二三区免费| 1000部很黄的大片| 一本一本综合久久| 婷婷精品国产亚洲av| 校园人妻丝袜中文字幕| 久久久久久久久久成人| 欧美一区二区亚洲| 午夜日韩欧美国产| 有码 亚洲区| 中文字幕人妻熟人妻熟丝袜美| 亚洲精品成人久久久久久| 乱系列少妇在线播放| 熟女电影av网| 免费人成视频x8x8入口观看| 欧美日韩瑟瑟在线播放| 少妇裸体淫交视频免费看高清| 天堂av国产一区二区熟女人妻| 久久久久久伊人网av| 午夜福利成人在线免费观看| 亚洲欧美日韩无卡精品| 日本黄色片子视频| 久久久午夜欧美精品| 国产精品嫩草影院av在线观看 | or卡值多少钱| 亚洲成av人片在线播放无| 十八禁网站免费在线| 国产aⅴ精品一区二区三区波| 久久亚洲精品不卡| 国产久久久一区二区三区| 国产午夜精品久久久久久一区二区三区 | 嫩草影院入口| 伦理电影大哥的女人| 国产黄色小视频在线观看| netflix在线观看网站| 麻豆国产97在线/欧美| 欧美一级a爱片免费观看看| 国产精品一及| 久久人妻av系列| 久久亚洲精品不卡| 别揉我奶头 嗯啊视频| 久久国产乱子免费精品| 国产亚洲精品久久久com| 观看美女的网站| 午夜激情欧美在线| 成年版毛片免费区| 美女高潮喷水抽搐中文字幕| 日韩亚洲欧美综合| 联通29元200g的流量卡| 国产精品1区2区在线观看.| 男女做爰动态图高潮gif福利片| 欧美xxxx黑人xx丫x性爽| 日本一本二区三区精品| 又粗又爽又猛毛片免费看| 听说在线观看完整版免费高清| 老司机午夜福利在线观看视频| 中文亚洲av片在线观看爽| 欧美精品国产亚洲| 男人狂女人下面高潮的视频| 成人毛片a级毛片在线播放| 一本久久中文字幕| 国产精品久久久久久久电影| 美女cb高潮喷水在线观看| 女同久久另类99精品国产91| 性色avwww在线观看| 又紧又爽又黄一区二区| 俺也久久电影网| 国产精品98久久久久久宅男小说| 亚洲精品影视一区二区三区av| 白带黄色成豆腐渣| 中国美女看黄片| 伦理电影大哥的女人| 一区二区三区四区激情视频 | 91精品国产九色| 美女大奶头视频| 最新在线观看一区二区三区| 久久国产乱子免费精品| 欧美日韩亚洲国产一区二区在线观看| 亚洲美女视频黄频| 午夜久久久久精精品| 亚洲综合色惰| 国产在线精品亚洲第一网站| 亚洲av免费在线观看| 又黄又爽又刺激的免费视频.| 91麻豆精品激情在线观看国产| 一本精品99久久精品77| 88av欧美| av专区在线播放| 有码 亚洲区| a级一级毛片免费在线观看| 国产一区二区在线av高清观看| 联通29元200g的流量卡| 91久久精品国产一区二区三区| av专区在线播放| 亚洲内射少妇av| 国产精品野战在线观看| 我要看日韩黄色一级片| 免费看光身美女| 女人十人毛片免费观看3o分钟| 亚洲无线观看免费| 亚洲美女视频黄频| 村上凉子中文字幕在线| 国产精品久久久久久av不卡| 久久久精品大字幕| 午夜福利欧美成人| 中文字幕精品亚洲无线码一区| 欧美最新免费一区二区三区| 亚洲欧美日韩卡通动漫| .国产精品久久| 亚洲狠狠婷婷综合久久图片| 日本一本二区三区精品| 五月伊人婷婷丁香| 国产亚洲精品久久久com| 亚洲午夜理论影院| 成人三级黄色视频| 欧美xxxx黑人xx丫x性爽| 中文资源天堂在线| 国内揄拍国产精品人妻在线| 真人一进一出gif抽搐免费| 看十八女毛片水多多多| 久久久久久久久久黄片| 别揉我奶头~嗯~啊~动态视频| 亚洲av二区三区四区| 精品午夜福利在线看| 欧美另类亚洲清纯唯美| 日韩精品青青久久久久久| 午夜福利在线观看吧| 国产精品99久久久久久久久| 国产精品三级大全| 天堂影院成人在线观看| 日韩欧美精品v在线| 99精品久久久久人妻精品| 黄色日韩在线| 男女之事视频高清在线观看| 窝窝影院91人妻| 日韩欧美一区二区三区在线观看| 日本五十路高清| 一区二区三区免费毛片| 久久久久久久久久成人| 99在线人妻在线中文字幕| 色视频www国产| 国产极品精品免费视频能看的| 狂野欧美激情性xxxx在线观看| 欧美另类亚洲清纯唯美| 日本免费一区二区三区高清不卡| 男插女下体视频免费在线播放| 日本熟妇午夜| 狂野欧美激情性xxxx在线观看| 国内精品久久久久精免费| 丰满的人妻完整版| 极品教师在线视频| 午夜福利在线在线| 国产真实乱freesex| 动漫黄色视频在线观看| 最近最新免费中文字幕在线| 一个人看视频在线观看www免费| 综合色av麻豆| 亚洲三级黄色毛片| 国产私拍福利视频在线观看| 精品久久久久久成人av| 美女免费视频网站| 成年免费大片在线观看| 国产精品嫩草影院av在线观看 | 欧美黑人欧美精品刺激| 少妇的逼水好多| 制服丝袜大香蕉在线| 国产白丝娇喘喷水9色精品| 高清毛片免费观看视频网站| 又粗又爽又猛毛片免费看| 久久精品国产99精品国产亚洲性色| 久久亚洲精品不卡| 免费观看精品视频网站| 少妇裸体淫交视频免费看高清| 国产精品98久久久久久宅男小说| 国产欧美日韩一区二区精品| 村上凉子中文字幕在线| 国产av一区在线观看免费| 欧美黑人欧美精品刺激| 在线观看舔阴道视频| 日韩av在线大香蕉| 麻豆久久精品国产亚洲av| x7x7x7水蜜桃| 热99在线观看视频| 成人综合一区亚洲| 日韩欧美免费精品| 国产老妇女一区| 乱系列少妇在线播放| 亚洲国产精品sss在线观看| 日韩欧美在线二视频| 欧美精品国产亚洲| 久久久久久久久久成人| 一本久久中文字幕| 婷婷精品国产亚洲av在线| 精品99又大又爽又粗少妇毛片 | 国产伦在线观看视频一区| 亚洲精华国产精华精| 悠悠久久av| 免费看av在线观看网站| 香蕉av资源在线| 看片在线看免费视频| 两个人视频免费观看高清| 久久精品91蜜桃| 一级黄片播放器| 日韩精品有码人妻一区| 成人综合一区亚洲| 日本爱情动作片www.在线观看 | 亚洲av免费在线观看| 亚洲精华国产精华液的使用体验 | 午夜亚洲福利在线播放| 中文亚洲av片在线观看爽| 美女 人体艺术 gogo| 久久精品久久久久久噜噜老黄 | 国产伦人伦偷精品视频| 国产视频一区二区在线看| 狂野欧美白嫩少妇大欣赏| 少妇被粗大猛烈的视频| 亚洲欧美日韩高清在线视频| 国产黄色小视频在线观看| 全区人妻精品视频| xxxwww97欧美| 久久久久久伊人网av| 小蜜桃在线观看免费完整版高清| 丝袜美腿在线中文| 国产爱豆传媒在线观看| 精品午夜福利视频在线观看一区| 亚州av有码| 久久久久免费精品人妻一区二区| www.www免费av| 亚洲av中文av极速乱 | 国产成人福利小说| 国产精品一区二区三区四区免费观看 | 熟女人妻精品中文字幕| 最新中文字幕久久久久| 国产av麻豆久久久久久久| 欧美激情在线99| 成人av一区二区三区在线看| 91精品国产九色| 天堂动漫精品| 看免费成人av毛片| 在线免费十八禁| 国产不卡一卡二| 免费人成在线观看视频色| 欧美成人性av电影在线观看| 12—13女人毛片做爰片一| 99久国产av精品| 久久精品影院6| 看免费成人av毛片| 成人av在线播放网站| 欧美bdsm另类| 日日摸夜夜添夜夜添小说| 亚洲国产精品成人综合色| 国内少妇人妻偷人精品xxx网站| 一卡2卡三卡四卡精品乱码亚洲| 天堂√8在线中文| 又黄又爽又免费观看的视频| 欧美色视频一区免费| 深夜a级毛片| 午夜精品一区二区三区免费看| 一夜夜www| 精品一区二区三区人妻视频| 欧美极品一区二区三区四区| 精品午夜福利在线看| 国产精品乱码一区二三区的特点| av女优亚洲男人天堂| 男女下面进入的视频免费午夜| 国产精品电影一区二区三区| 男女视频在线观看网站免费| 嫩草影视91久久| av黄色大香蕉| 免费看美女性在线毛片视频| 亚洲一区二区三区色噜噜| 亚洲精品456在线播放app | 欧美绝顶高潮抽搐喷水| 好男人在线观看高清免费视频| 3wmmmm亚洲av在线观看| 联通29元200g的流量卡| 啪啪无遮挡十八禁网站| 久久精品人妻少妇| 午夜久久久久精精品| 在线观看美女被高潮喷水网站| 成年女人看的毛片在线观看| 欧美成人免费av一区二区三区| 麻豆国产97在线/欧美| 婷婷精品国产亚洲av| 日本撒尿小便嘘嘘汇集6| 97超视频在线观看视频| 天天一区二区日本电影三级| 熟女人妻精品中文字幕| 久久久久久久久大av| 欧美人与善性xxx| 99在线视频只有这里精品首页| 欧美bdsm另类| 免费av观看视频| 99久久精品国产国产毛片| 亚洲成人中文字幕在线播放| 美女被艹到高潮喷水动态| 久久天躁狠狠躁夜夜2o2o| 亚洲乱码一区二区免费版| 成人综合一区亚洲| 欧美成人a在线观看| 日本五十路高清| 午夜精品一区二区三区免费看| 免费人成视频x8x8入口观看| 精品人妻1区二区| 淫妇啪啪啪对白视频| 久久久久久久精品吃奶| 日日啪夜夜撸| 久久久精品大字幕| 男女视频在线观看网站免费| av在线观看视频网站免费| 日韩大尺度精品在线看网址| 亚洲精品影视一区二区三区av| 欧美日本亚洲视频在线播放| 中文字幕久久专区| 国产精品久久久久久av不卡| 欧美日本亚洲视频在线播放| 免费人成在线观看视频色| 97热精品久久久久久| 搡老岳熟女国产| 国产在线男女| 毛片女人毛片| 久久亚洲精品不卡| 欧美日韩国产亚洲二区| 简卡轻食公司| 免费电影在线观看免费观看| 看免费成人av毛片| 韩国av在线不卡| 国产av不卡久久| 99视频精品全部免费 在线| 一级av片app| 久久久久久大精品|