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

    [Co3]-Cluster Based Metal-Organic Framework Enables “Two Birds with One Stone” in Efficient Transformation of CO2 to Oxazolidinones

    2023-12-08 03:15:08ZhuohaoJiaoXinyuanZhaoJianZhaoYaoXieShengliHouBinZhao
    物理化學(xué)學(xué)報(bào) 2023年11期

    Zhuohao Jiao,Xinyuan Zhao,Jian Zhao,Yao Xie,Shengli Hou ,Bin Zhao

    College of Chemistry,Key Laboratory of Advanced Energy Material Chemistry,Ministry of Education,Renewable Energy Conversion and Storage Center,Nankai University,Tianjin 300071,China.

    Abstract: CO2,as a greenhouse gas,has excessive emissions that lead to many environmental problems and is a rich and cheap C1 resource.Effective utilization and transformation of CO2 has become an important means of achieving carbon neutrality.Oxazolidinones are important intermediates in pharmaceutical chemistry that can be synthesized by carboxylation cyclization of CO2 with propargyl amines or cycloaddition of CO2 with aziridines.Owing to CO2’s high stability,these reactions typically require harsh conditions,such as high temperatures or pressures.It is desirable,but challenging,to find a catalyst that can catalyze these two types of reactions under relatively mild conditions.Metal-organic frameworks (MOFs) are an emerging class of heterogeneous catalysts that with great potential in the catalytic conversion of CO2 to value-added products because of their attractive features,such as abundant metal active sites,inherent porosity,and easy functionalities.Herein,a unique three-dimensional (3D) MOF,{(CH3NH2CH3)2[Co3(BCP)2]·6H2O·4DMF}n (1) (H4BCP: 5-(2,6-bis(4-carboxyphenyl)pyridin-4-yl) isophthalic acid; DMF: N,N'-dimethylformamide),was synthesized using carboxylic acid ligands and cobalt salts via a solvothermal method.According to structural analysis,[Co3] clusters as secondary building units (SBU) are bridged by BCP4- ligands,forming an anion framework with flu topology,and dimethylamine cations act as counter ions in the pores.The framework has rectangular channels of approximately 0.4 nm × 0.9 nm along the a-axis direction,exhibiting its porous property.Infrared spectroscopy (IR) and X-ray photoelectron spectroscopy (XPS) characterizations proved the coordination interaction between the carboxyl groups in the ligands and the metal ions.The powder X-Ray diffraction (PXRD) test further confirmed the phase purity of the synthesized samples.PXRD and thermogravimetry (TG) analyses indicated that 1 possessed good solvent and thermal stabilities.The catalytic experiments revealed that 1 could effectively catalyze CO2 with aziridines or propargyl amines to prepare oxazolidinones.In the cycloaddition of CO2 with aziridines,1 can facilitate the reaction under relatively mild conditions compared to other reported MOF-based catalysts.It shows excellent universality for substrates with various substitutions on the N atom or benzene ring.Investigation of the mechanism indicated that the coordination interaction of cobalt metal sites with the nitrogen atoms of aziridines can activate the substrates.For the carboxylative cyclization of CO2 with propargylic amines,this catalyst also has a broad substrate scope.Control experiments and nuclear magnetic resonance (NMR) tests suggest that Lewis acid metal sites are responsible for the high catalytic efficiency achieved by activating the alkyne groups.Moreover,1 showed good reusability in both reactions.Compound 1 represents a new catalyst that enables “two birds with one stone” in the catalytic synthesis of oxazolidinones using CO2.

    Key Words: Metal-organic framework; Bifunctional catalyst; CO2 utilization; Oxazolidinone; Reusability

    1 Introduction

    With the proposal of carbon neutrality,the emission reduction and effective utilization of CO2have become the focus of researches.CO2as an inexpensive and readily available C1 feedstock,can be converted into a variety of high value-added products,including urea,formic acid,methanol,cyclic carbonates,oxazolidinones,and polymers through photocatalysis,electro-catalysis,or thermos-catalysis1–11.Among them,oxazolidinones as heterocyclic compounds have been widely used in synthesizing antibacterial drugs such as tedizolid,radezolid,and linezolid12.And they also play an important role in serving as intermediates in organic synthesis and chiral auxiliary reagents13,14.Heretofore,diverse methods have been proposed for the preparation of oxazolidinones,among which two of the most common and efficient approaches are the cycloaddition of CO2with aziridines or carboxylation cyclization of CO2with propargyl amines due to their high atom economy and simple procedure15.Some catalysts have been applied in preparation of oxazolidinones by these two approaches,including amino acids16,haloids17,noble metal complexes18,etc.However,high cost or complicated postprocessing procedures limit their large-scale industrial applications,and these studies only focus on one kind of reaction on account of the high thermodynamic stability and kinetic inertness of CO2.Hence,it is important to construct a heterogeneous catalyst that can catalyze these two kinds of reactions.

    Metal-organic frameworks (MOFs),as an emerging class of heterogeneous catalysts,have gained increasing attention due to their abundant metal active sites19–22,inherent porosity23–25,and easy functionalities26–29.In recent years,researches focused on MOFs for catalyzing the formation of oxazolidinones from CO2with aziridines or propargyl amines,and many exciting works have been reported successively30,31.However,these catalysts are usually used to catalyze only one type of these reactions because these two reactions are Lewis-acid and Lewis-base preferred reactions,respectively.In addition,due to the high stability of CO2,these reactions usually require harsh conditions such as high temperature or pressure.Hence,the construction of one MOF that can simultaneously catalyze these two types of reactions under mild conditions is desirable but challenging.

    The Co2+ion has been used to construct a variety of MOFs structures due to its diverse coordination modes,and its Lewis acidity also allows it to apply in numerous catalytic reactions32–37.In this contribution,we synthesized a unique 3D MOF{(CH3NH2CH3)2[Co3(BCP)2]·6H2O·4DMF}n(1) by combining Co2+and a carboxylate ligand 5-(2,6-bis(4-carboxyphenyl)pyridin-4-yl) isophthalic acid (H4BCP) through solvothermal method,which has a large 1D channel with a size of 0.4 nm ×0.9 nm.Experiment results reveal that 1 shows excellent catalytic performance for promoting the conversion of both aziridines and propargylic amines with CO2into oxazolidinones,and it can be reused for five cycles and maintains high catalytic activity.In particular,1 enables “two birds with one stone” in the synthesis of oxazolidinones by CO2,which represents a rarely reported catalyst for preparing oxazolidinones through two pathways.

    2 Experimental and computational section

    2.1 Preparation of {( CH3NH2CH3)2[Co3(BCP)2]·6H2O·4DMF}n (1)

    Co(NO3)2·6H2O (99.99%,15 mg) and H4BCP (20 mg) were mixed withN,N'-dimethylformamide (DMF,99.5%,2.5 mL)and deionized water (0.5 mL).Then 200 μL of nitric acid(65.0%–68.0%) was dropped in the mixture to adjust the pH.Subsequently,the mixture was transferred to a 25 mL Teflonlined stainless-steel vessel and heated at 145 °C for 72 h under autogenous pressure.After cooling to room temperature naturally,dark purple lamellar crystals were collected.The resultant crystals were washed with DMF and H2O several times to remove the residual metal salts or ligands.Elemental analysis of crystal 1,calculated results: C 51.33%,H 4.83%,N 6.84%.Tested results: C 51.93%,H 4.79%,N 6.89%.CCDC number:2160038.

    2.2 Catalyst characterization

    Powder X-ray diffraction (PXRD) patterns were measured on an Ultima IV diffractometer (Rigaku,Japan) using CuKαas radiation sources with a 5 (°)·min-1scanning speed.X-ray photoelectron Spectrometer (XPS) analysis was carried out on an ESCALAB 250Xi spectrometer (ThermoFisher,USA) using Mono AlKαradiation,and the internal standard carbon is 284.8 eV.1H NMR chemical shifts were measured on 400 MHz spectrometer (ASCEND 400,Bruker,Germany) using CDCl3as internal standard.13C NMR chemical shifts were carried out on a AVANCE III HD 400 spectrometer (101 MHz,Bruker,Germany).The metal content was tested by Aglient 5110 (USA)inductively coupled plasma optical emission spectrometer (ICPOES).Typically,100 μL of reaction mixture obtained by centrifugation is decomposed by aqua regia and then diluted to 5 mL by deionized water.The resultant sample was measured by ICP-OES.Infrared spectra (IR) were measured by a Bruker Tensor 27 (Germany) spectrophotometer.The thermogravimetric analysis (TGA) data were measured on thermal gravimetric analyzer (NETZSCH TG209,Germany) under air atmosphere.Single-crystal X-ray diffraction was performed on a SuperNova Single Crystal diffractometer (Rigaku,Japan) with MoKαradiation (λ= 0.71073 ?,1 ? = 0.1 nm).The obtained data was analyzed by Olex2 software to get the structure38.Due to the highly disordered state of solvent molecules in frameworks,we use PLATON/SQUEEZE programs to remove the disordered solvent molecules39,and the number of free solvent molecules was determined by elemental analysis and TGA.Detail of single crystal data and structure refinement for 1 were recorded in Table S1 in Supporting Information.

    2.3 Catalytic experiments

    Cycloaddition of CO2with aziridines: Compound 1 (10 mg),aziridine substrates (0.5 mmol) and co-catalyst TBAB (32.4 mg,0.1 mmol) were sealed in an autoclave.After purging with 0.5 MPa CO2,the autoclave was placed in water bath of 30 °C and stirred for 9 h.When the reaction was completed,the catalyst was separated by centrifugation and filtration.The filtrate was dissolved in 0.5 mL CH2Cl2(99.9%),and then the yields of corresponding products (2x and 3x) were analyzed by1H NMR with 1,3,5-trimethoxybenzene (99.9%) as an internal standard.

    Carboxylation cyclization of CO2with propargylic amines:Compound 1 (10 mg),propargylic amines (0.5 mmol) and additive 1,8-diazabicyclo-[5.4.0]undec-7-ene (DBU,98%,19 mg,0.125 mmol) were mixed in CH3CN (99.7%,0.5 mL).After adding the mixture in a 15 mL Schlenk tube and purging 0.1 MPa of CO2,the tube was placed in oil bath at 80 °C and kept stirred for 10 h.The mixture was separated by centrifugation and filtration when the reaction was completed,and the filtrate was analyzed by1H NMR with 1,3,5-trimethoxybenzene as an internal standard.

    Cycle experiments: The filter residue was washed with fresh CH3CN three times and dried under vacuum at 60 °C.The recovered catalyst was used for the next catalytic experiments.

    3 Results and discussion

    3.1 Characterization

    According to single-crystal X-ray diffraction analysis,compound 1 crystallizes in a triclinic system withP-1 space group.The asymmetric unit contains one and a half Co2+ions,one independent BCP4-ligand,and one dimethylamine cation([(CH3)2NH2]+).There are two types of Co ions (denoted as Co1 and Co2),and both of them are coordinated with six oxygen atoms from the carboxyl group of ligands (O1,O1A,O3,O3A,O7,O7A for Co1,and O3,O5,O9 O12,O13,O16 for Co2) (Fig.S1).One Co1 and two Co2 ions are bridged by the carboxyl group to form a [Co3] cluster (Fig.1a).Each ligand coordinates with four [Co3] clusters,expanding into a 3D anion framework(Fig.1b,c).[(CH3)2NH2]+acts as counterion to balance the charge which may generate from the decomposition of DMF and its existence was confirmed by structural analysis40.By considering the [Co3] clusters and H4BCP ligands as nodes,the structure can be rationalized to a 4,8-connected net withflutopology (Fig.1d).Along thea-axis direction,a rectangular channel about 0.4 nm × 0.9 nm is present in the framework,and the total potential pore volume of the structure is estimated to be 39.6% by PLATON software (Fig.1e).

    Fig.1 (a) The 8-connected [Co3] cluster SBU; (b) coordination mode of H4BCP; (c) 3D structure of the framework along the a-axis; (d) flu topology of 1; (e) the rectangular channel of 1 along the a-axis direction.

    According to the analysis of infrared spectroscopy (IR),the ―OH vibration peak of the carboxyl group in the ligand is located at 2940–3300 cm-1,while it disappears in the spectrum of 141.Meanwhile,the peak at 1672.8 cm-1which belongs to the stretching vibration of C=O shifted to 1656.4 cm-1in that of 1,which demonstrates the coordination of carboxyl groups to the Co ions (Fig.S2)42.

    XPS analysis shows the fitting peak at the binding energy of 531.78 eV in O 1sspectrum corresponds to ―C―O―Co in the framework,conforming the coordination interaction between carboxyl groups and Co ions (Fig.S3c)43.In Co 2pspectrum,the peaks located at 781.7 and 797.4 eV belong to Co 2p3/2and Co 2p1/2of Co2+,while the binding energy of 786.2 and 802.0 eV is in accordance with the satellite peaks of Co2+(Fig.S3d)44.The PXRD peaks of 1 are identical to those of the simulated one,demonstrating the satisfactory phase purity of synthesized crystals (Fig.S4a).After immersing the crystals in several common solvents (DMF,N,N'-dimethylacetamide,methanol,ethanol,acetonitrile,dichloromethane,dioxane,Nmethylpyrrolidone) for 24 h,the patterns still agree with that of synthesized samples and standard patterns,which confirms the excellent solvent stability of 1 (Fig.S4b).In addition,the thermal stability was explored by performing the thermogravimetric analysis (TGA) at 30–700 °C (Fig.S5).As the temperature rises,there is a weight loss about 25.05% before 330 °C,which is assigned to the escape of free solvent molecules(H2O and DMF,calculated 25%).Then the framework begins to gradually collapse and decompose.Meanwhile,PXRD test was performed under different temperatures (Fig.S6).Combined with TGA and PXRD results,1 can maintain its structure integrity under 200 °C,which indicates that this compound has good thermal stability.The favorable thermal and solvent stabilities provide a guarantee for compound 1 to be used for thermal catalysis.

    3.2 Catalytic conversion of CO2

    3.2.1 Cycloaddition with aziridines

    Benefiting from the porous nature of the MOF material,the catalyst may enrich the substrates and CO2.Hence,we tested the CO2adsorption capacity of the activated samples at 298 K (Fig.S7).The adsorption-desorption isotherms display that the CO2absorption is 10.68 cm3·g-1at 298 K,which indicates that 1 possesses the ability to adsorb CO2in thermal catalysis.Owing to the attractive stability and adsorption ability of this compound,we conducted to study the MOF-catalysed conversion of CO2to oxazolidinones.Initially,we explored the performance of 1 in catalytic reaction of CO2with aziridines,and 1-ethyl-2-phenylaziridine (1a) was chosen as an example to explore the optimal reaction conditions.The effect of pressure was screened at first and the target product 3-ethyl-5-phenyloxazolidin-2-one (2a) was obtained with 19% NMR yield under atmospheric pressure.By increasing pressure to 0.5 MPa,the substrate could be quantitatively converted to the desired compound (entries 1–3,Table 1).As the temperature decreased,the substrate could be completely converted at 40 or 30 °C,but when the temperature dropped to room temperature (25 °C),the yield could only reach 70% (entries 4–6,Table 1).With increasing reaction time,the yield gradually increased,and 99% yield could be obtained in 9 h (entries 7 and 8,Table 1).Therefore,the optimal reaction condition is considered as 10 mg of catalyst 1 with 0.1 mmol of TBAB as co-catalyst,0.5 mmol of substrates,and 0.5 MPa CO2for 9 h at 30 °C in solvent-free conditions.Compared to the reported MOF-based catalysts,1 can facilitate the reaction under relatively mild conditions.(Table S2).

    Table 1 Catalytic reaction of CO2 with 1-ethyl-2-phenylaziridine under different conditions a.

    The universality of 1 with different aziridines was discussed on the basis of the optimal reaction conditions.For the alkyl substituents from methyl ton-pentyl on the N atom,the yields of corresponding products decreased from 99%–79%,which is mainly attributed to the increase of size for R1 group hindering the proceeding of the reaction (entries 1,2,4–6,Table 2).In addition,when benzyl substituent with larger steric hindrance was introduced on the N atom,the reaction yield decreased to 83%,which further proves the steric hindrance effect for the reaction (entry 7,Table 2).Besides,isopropyl substituent led to a relatively lower yield (70%),probably because the larger bulk of the N-atom-connected propyl group hinders the coordination of metal sites with N atoms (entry 3,Table 2)45.As for the substituents at theparaposition of the benzene ring,the electron-donating substituent ―CH3had almost no effect on the yield and still gave excellent yields.In contrast,the electronwithdrawing substituents such as chloro (―Cl) and bromo (―Br)reduced the reaction efficiency to some extent,resulting in the corresponding product 5-(4-chlorophenyl)-3-ethyloxazolidin-2-one (2j) and 5-(4-bromophenyl)-3-ethyloxazolidin-2-one (2k) of 78% and 75% yields,respectively (entries 8–10,Table 2).To sum up,compound 1 can successfully prepare oxazolidinonesviathe pathway of cycloaddition of aziridines and CO2with broad substrate scope.

    Table 2 Exploration of substrate universality of catalyst 1 under optimal conditions a.

    Recyclability is of crucial importance for heterogeneous catalysts in practical applications.By simple centrifugation and filtration,the catalyst was recovered and applied to the next catalytic experiment.1 preserved high catalytic activity and only a slight decrease in yield was observed after five continuous reaction cycles (Fig.2).The catalyst recovered after five cycles was analysed by PXRD.The pattern agreed well with that of the original sample,indicating that 1 can maintain stable during the reaction (Fig.S8,The ICP-OES tests showed that only a small amount of Co ions leached into the reaction mixture (calculated to be 1.39% leakage of Co2+),which could be caused by the destruction of trace amounts of MOF during vigorous stirring(Table S3)46.All the results above show that 1 maintains its efficiency and robustness during cyclic reactions.

    Fig.2 The cycling performance of compound 1 in catalytic reactions.

    To explore the possible mechanism,we conducted control experiments with different conditions (Table S4).It is apparent that the ligand H4BCP could not make the experiment proceed,and the cobalt nitrate and co-catalyst TBAB could only obtain 30% or 48% yield,respectively.The mechanical mixture of these three components also only yielded 68% of desired products,which is much lower than that of 1.The results indicate that the binding effect between free Co2+and Br-may inhibit the nucleophilicity of Br-,thereby weakening the ability of Br-to attack the aziridines.Combining the above findings and literature reports30,47,a possible mechanism is proposed for the cycloaddition reaction of the aziridines with CO2.Firstly,the porous structure of 1 enriches CO2and substrates,and cobalt metal sites coordinate with nitrogen atoms of the aziridines to activate the substrates; secondly,the nucleophilic attack of Br–at the C atoms of aziridines leads to ring-opening by two pathways which are presented as I and II; thirdly,CO2is inserted into the N―Co bond to form a carbamate salt intermediate;finally,the desired oxazolidinone product is obtained by an intramolecular nucleophilic substitution ring-closure reaction with the regeneration of 1 and TBAB (Scheme 1).As for regioselectivity,path I is the preferred one,mainly because Br–tends to attack the aziridines from the more substituted carbon with the generation of a more stable carbamate salt than path II.Thus,the 5-substituted oxazolidinones can be obtained in excellent yields and regioselectivity when 1 is used as the catalyst.

    Scheme 1 Mechanistic proposal for cycloaddition reaction of CO2 and aziridines.

    3.2.2 Carboxylative cyclization with propargylic amines

    Then we turned to investigate whether 1 could act as a bifunctional catalyst to promote the carboxylative cyclization of propargylic amines with CO2.At the outset,we choseN-benzyl-2-propynylamine (A1) as the model compound and the reaction was carried out using 1 (10 mg) as catalyst and DBU (0.25 mmol) as base in CH3CN at 40 °C under atmospheric CO2pressure.Under given conditions,the target product 3-benzyl-5-methyleneoxazolidin-2-one was detected with an NMR yield of 19% (entry 1,Table S5).Subsequently,the influence of reaction temperature (40,60 and 80 °C) was explored and the results indicated that the rise of temperature significantly promoted the reaction,with a higher yield of 76% at 80 °C (entries 2 and 3,Table S5).Then the influence of time was examined,and when the reaction time was extended to 10 h,the substrate was almost quantitatively converted to the desired product (entry 4,Table S5).Under the above-optimized conditions,the reaction still proceeds smoothly with an excellent yield (97%) after reducing the amount of DBU to 0.125 mmol (entry 5,Table S5).Noticeably,when the reaction performed without solvent,it still achieved a satisfactory yield up to 85% (entry 6,Table S5).Therefore,we choose 10 mg of catalyst 1 with 0.125 mmol of DBU as base at 80 °C and 1 atm (1 atm = 101325 Pa) of CO2in 0.5 mL CH3CN for 10 h as the optimal reaction conditions.In addition,control experiments show that compound 1 exhibits superior catalytic performance compared with the mixture of Co(NO3)2and H4BCP,suggesting that the framework of 1 plays an important role in this catalysis (Table S6).

    Similarly,the catalytic activity of 1 was examined with various propargyl amines (Table 3).The substrates bearing electron-donating substituent methoxy (―OCH3) or electron-withdrawing groups like chloro (―Cl) and bromo (―Br) at thepara-position of the benzyl group could smoothly transform to corresponding oxazolidinone products with excellent yields (B2 93%,C2 99%,D2 91%,entries 2–4,Table 3).However,whentert-butyl was introduced on the benzene ring,there was a significant decrease in yield to 64%,probably caused by the large steric hindrance effect oftert-butyl group (entry 5,Table 3).In addition,when there are alkyl substituents such as isopropyl,n-butyl,and cyclohexyl on the N atom of propargyl amines,the corresponding products were detected with 92%–95% yield(entries 6–8,Table 3).Thus,1 also exhibits satisfactory substrate universality for the carboxylative cyclization of propargylic amines with CO2.

    Table 3 Preparation of various oxazolidinones with CO2 and propargylic amines a.

    Using the same recycling procedures above,we investigated the recycling performance of 1 for this approach.1 could be recovered after 5 cycles of reactions and maintained its catalytic activity for catalysing the reaction in good yields (Fig.2).The PXRD of the catalyst still matched with the simulated data,illustrating the integrity of the framework (Fig.S7).Furthermore,ICP tests showed that only about 1.16% Co leaked into the reaction mixture after five continuous catalytic runs(Table S3).These results demonstrate that 1 can act as a heterogeneous catalyst for carbonylative cyclization of CO2and propargyl amines with good recyclability and stability.

    To investigate the reaction mechanism,we attempted to monitor the activation of the propargyl amines by 1 and DBU through1H NMR and13C NMR spectroscopy.The N―H signal of A1 at 2.63 disappeared when DBU was added,which indicates that DBU can activate the hydrogen proton of the N―H group in substrates (Fig.S9)48.At the same time,the strength of the13C signals atδ= 73.7 and 82.8 assigned to the alkyne group was obviously weakened after the addition of 1,suggesting that the C≡C bond of propargyl amine can be activated by 1 (Fig.S10).Based on the NMR results and literature research49,50,a possible mechanism is proposed for the carboxylative cyclization of propargylic amines with CO2.Initially,the imino and alkyne groups of propargyl amine A1 are activated by DBU and 1,respectively; Then the activated substrate further reacts with CO2to generate the carbamate intermediate; Subsequently,the negatively charged oxygen ion undergoes an intramolecular ring closure reaction through nucleophilic attack; Finally,proto-demetallation occurs with the help the proton carrier to obtain the desired product A2,accompanied by the regeneration of DBU and 1 (Scheme S1).

    4 Conclusions

    In summary,a unique 3D Co-based MOF 1 with 1D channels was successfully synthesized.Compound 1 features good solvent and thermal stability.The results of catalytic experiments indicated that 1 can serve as an effective heterogeneous catalyst for promoting the cycloaddition of CO2with aziridines or carboxylative cyclization of CO2with propargylic amines to prepare oxazolidinones.The catalyst presents broad substrate universality for both reactions.Besides,it can be reused for five cycles without significant decrease in catalytic efficiency.Compound 1 represents a rare catalyst that enables “two birds with one stone” in catalytic synthesis of oxazolidinones by CO2,which may expand the application of MOFs in the field of catalysis and the potential of MOFs to catalyse multiple reactions simultaneously.

    Supporting Information: Materials,control experiments details,and IR,TGA,PXRD,CO2adsorption-desorption,1H NMR and13C NMR spectra (PDF); X-ray crystallographic data of 1 (CIF) available free of chargeviathe internet at http://www.whxb.pku.edu.cn.

    欧美日韩精品成人综合77777| 啦啦啦在线免费观看视频4| 久久人人爽av亚洲精品天堂| 我要看黄色一级片免费的| 久久精品久久久久久久性| 亚洲精品久久成人aⅴ小说| 国产av精品麻豆| 午夜日本视频在线| 在线观看国产h片| 亚洲综合精品二区| 精品99又大又爽又粗少妇毛片| 国产不卡av网站在线观看| 中文字幕人妻丝袜制服| 一本—道久久a久久精品蜜桃钙片| 日韩成人av中文字幕在线观看| 97精品久久久久久久久久精品| 久久国产精品男人的天堂亚洲| 欧美最新免费一区二区三区| 99热网站在线观看| 丰满少妇做爰视频| 久久精品国产a三级三级三级| 久久久久精品人妻al黑| 看非洲黑人一级黄片| 中国国产av一级| 在现免费观看毛片| 中文字幕亚洲精品专区| 亚洲在久久综合| 亚洲激情五月婷婷啪啪| 国产亚洲午夜精品一区二区久久| 亚洲欧美中文字幕日韩二区| 80岁老熟妇乱子伦牲交| 午夜福利网站1000一区二区三区| 成年女人毛片免费观看观看9 | 大香蕉久久网| 最新中文字幕久久久久| 色播在线永久视频| 99久久人妻综合| 成年人免费黄色播放视频| 青草久久国产| xxxhd国产人妻xxx| 丝袜在线中文字幕| 亚洲经典国产精华液单| 亚洲国产av影院在线观看| 看免费av毛片| 精品卡一卡二卡四卡免费| 最近中文字幕2019免费版| 天天躁狠狠躁夜夜躁狠狠躁| 国产精品久久久久久久久免| 自线自在国产av| 免费黄频网站在线观看国产| 秋霞伦理黄片| 人体艺术视频欧美日本| 亚洲一区二区三区欧美精品| 美女主播在线视频| 各种免费的搞黄视频| 一级毛片电影观看| 国产av国产精品国产| 亚洲av免费高清在线观看| 激情视频va一区二区三区| 精品99又大又爽又粗少妇毛片| 午夜av观看不卡| 视频区图区小说| av在线播放精品| 999精品在线视频| 久久av网站| 亚洲人成电影观看| 亚洲国产精品一区三区| 亚洲美女搞黄在线观看| 中文字幕另类日韩欧美亚洲嫩草| 免费久久久久久久精品成人欧美视频| 一边摸一边做爽爽视频免费| 欧美激情极品国产一区二区三区| 国产成人免费观看mmmm| 久久精品国产综合久久久| 国产深夜福利视频在线观看| 成年美女黄网站色视频大全免费| 国产精品人妻久久久影院| 久久99热这里只频精品6学生| 久久久久精品性色| videos熟女内射| 日本黄色日本黄色录像| 久久久久精品久久久久真实原创| 国产成人免费无遮挡视频| 亚洲精品成人av观看孕妇| videos熟女内射| 日本-黄色视频高清免费观看| 26uuu在线亚洲综合色| 精品一区在线观看国产| 国产精品蜜桃在线观看| 亚洲精品美女久久av网站| 日韩中字成人| 国产精品女同一区二区软件| 男人舔女人的私密视频| 亚洲欧美精品自产自拍| 国产一区亚洲一区在线观看| 国产日韩欧美视频二区| 久久精品国产亚洲av涩爱| 国产成人精品久久二区二区91 | 热re99久久精品国产66热6| 在线观看一区二区三区激情| 这个男人来自地球电影免费观看 | 久久午夜福利片| 乱人伦中国视频| 欧美精品av麻豆av| videossex国产| 涩涩av久久男人的天堂| 久久国产精品大桥未久av| av免费在线看不卡| 性色avwww在线观看| 欧美av亚洲av综合av国产av | 精品国产乱码久久久久久男人| 波多野结衣一区麻豆| 国精品久久久久久国模美| 久久久精品区二区三区| 欧美日韩综合久久久久久| 五月伊人婷婷丁香| av国产久精品久网站免费入址| 中文字幕制服av| 久久国产精品大桥未久av| 69精品国产乱码久久久| 婷婷色综合www| 久久精品国产亚洲av天美| 免费看不卡的av| 大香蕉久久成人网| 久久精品久久精品一区二区三区| 波多野结衣av一区二区av| 日韩,欧美,国产一区二区三区| 最黄视频免费看| 9色porny在线观看| 婷婷成人精品国产| 1024香蕉在线观看| 国产精品一国产av| 美女福利国产在线| 国产精品久久久久久久久免| 美女脱内裤让男人舔精品视频| 99热国产这里只有精品6| 好男人视频免费观看在线| 国产精品蜜桃在线观看| 欧美日韩综合久久久久久| 菩萨蛮人人尽说江南好唐韦庄| 亚洲内射少妇av| 九色亚洲精品在线播放| 久久这里有精品视频免费| 亚洲精品视频女| 亚洲内射少妇av| 在线观看www视频免费| 蜜桃在线观看..| 日本色播在线视频| 大片免费播放器 马上看| 亚洲欧美色中文字幕在线| 国产高清不卡午夜福利| 中文字幕色久视频| 一边摸一边做爽爽视频免费| 欧美日韩综合久久久久久| 七月丁香在线播放| 亚洲三区欧美一区| 久久久久国产网址| 国产精品免费大片| 国产 精品1| 人人妻人人爽人人添夜夜欢视频| 哪个播放器可以免费观看大片| 精品少妇内射三级| 国产精品麻豆人妻色哟哟久久| 久久免费观看电影| 久久影院123| 18禁观看日本| 国产精品久久久av美女十八| 婷婷成人精品国产| 久久国产精品男人的天堂亚洲| 天天躁狠狠躁夜夜躁狠狠躁| 成人午夜精彩视频在线观看| 亚洲人成网站在线观看播放| av.在线天堂| 汤姆久久久久久久影院中文字幕| 99香蕉大伊视频| 国产一区二区 视频在线| 午夜91福利影院| 亚洲美女搞黄在线观看| 少妇精品久久久久久久| 不卡视频在线观看欧美| 国产精品人妻久久久影院| a级毛片黄视频| 秋霞伦理黄片| 两性夫妻黄色片| 一级毛片我不卡| 亚洲伊人久久精品综合| 曰老女人黄片| 尾随美女入室| av网站在线播放免费| 亚洲精品中文字幕在线视频| 久久精品人人爽人人爽视色| 七月丁香在线播放| 黄色 视频免费看| 女人久久www免费人成看片| 国产一区亚洲一区在线观看| 亚洲三区欧美一区| 午夜福利在线观看免费完整高清在| 亚洲成国产人片在线观看| 侵犯人妻中文字幕一二三四区| 国产成人午夜福利电影在线观看| 国产一区亚洲一区在线观看| 超碰97精品在线观看| 亚洲成人av在线免费| 亚洲一级一片aⅴ在线观看| 99久久综合免费| 久久精品熟女亚洲av麻豆精品| 精品福利永久在线观看| 久久午夜福利片| 新久久久久国产一级毛片| 午夜日本视频在线| 激情五月婷婷亚洲| 亚洲欧美色中文字幕在线| 满18在线观看网站| 嫩草影院入口| 成年女人在线观看亚洲视频| 久久精品久久久久久噜噜老黄| 日日摸夜夜添夜夜爱| 十八禁高潮呻吟视频| 美女高潮到喷水免费观看| 午夜激情久久久久久久| 欧美国产精品va在线观看不卡| 久久精品熟女亚洲av麻豆精品| 久久久亚洲精品成人影院| 综合色丁香网| 国产亚洲精品第一综合不卡| 婷婷色综合大香蕉| 欧美 亚洲 国产 日韩一| 在线亚洲精品国产二区图片欧美| 婷婷色麻豆天堂久久| 亚洲精品国产色婷婷电影| 久久精品国产a三级三级三级| 久久国产精品大桥未久av| av福利片在线| 亚洲av免费高清在线观看| 18+在线观看网站| 久久99蜜桃精品久久| 亚洲欧美清纯卡通| 日韩在线高清观看一区二区三区| 亚洲av.av天堂| 国产老妇伦熟女老妇高清| 亚洲国产日韩一区二区| 日韩三级伦理在线观看| 伊人久久国产一区二区| 亚洲视频免费观看视频| 亚洲男人天堂网一区| 在线观看www视频免费| 久久精品夜色国产| 一本久久精品| 亚洲一区中文字幕在线| 自线自在国产av| 欧美激情高清一区二区三区 | 亚洲欧美成人精品一区二区| 人人妻人人爽人人添夜夜欢视频| 黄频高清免费视频| 久久久久久伊人网av| 在线观看国产h片| 免费观看无遮挡的男女| h视频一区二区三区| 国产亚洲精品第一综合不卡| 亚洲第一av免费看| 精品少妇黑人巨大在线播放| 成人毛片a级毛片在线播放| 99热网站在线观看| 我要看黄色一级片免费的| 欧美成人精品欧美一级黄| 国产一区有黄有色的免费视频| 波野结衣二区三区在线| 亚洲美女黄色视频免费看| 久久精品久久精品一区二区三区| 97精品久久久久久久久久精品| 午夜福利视频在线观看免费| 18禁国产床啪视频网站| 日韩三级伦理在线观看| 成人二区视频| 精品人妻熟女毛片av久久网站| 国产精品免费大片| 精品亚洲乱码少妇综合久久| 国产精品 国内视频| 亚洲情色 制服丝袜| 欧美在线黄色| 美女主播在线视频| 久久久精品区二区三区| 一级毛片黄色毛片免费观看视频| 丝袜美足系列| 18+在线观看网站| 啦啦啦在线观看免费高清www| 99热全是精品| 午夜福利在线免费观看网站| 久久 成人 亚洲| 亚洲第一av免费看| 七月丁香在线播放| 国产成人欧美| 亚洲一区二区三区欧美精品| 久久午夜福利片| 亚洲国产日韩一区二区| √禁漫天堂资源中文www| 丁香六月天网| av福利片在线| 国产 精品1| kizo精华| 亚洲国产精品999| 亚洲av福利一区| 一级毛片电影观看| 亚洲欧美中文字幕日韩二区| 伊人久久国产一区二区| 91精品三级在线观看| 99热国产这里只有精品6| 综合色丁香网| 亚洲精品成人av观看孕妇| 久久精品熟女亚洲av麻豆精品| 亚洲av中文av极速乱| 成人国语在线视频| videossex国产| 中文字幕制服av| 久久鲁丝午夜福利片| av片东京热男人的天堂| av电影中文网址| 国产不卡av网站在线观看| 亚洲av成人精品一二三区| 亚洲成人一二三区av| 老女人水多毛片| 黄色一级大片看看| 亚洲综合色网址| 国产精品嫩草影院av在线观看| 国产成人午夜福利电影在线观看| 成人黄色视频免费在线看| 久久综合国产亚洲精品| 亚洲国产精品国产精品| 你懂的网址亚洲精品在线观看| 最新中文字幕久久久久| 美女国产高潮福利片在线看| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 亚洲国产毛片av蜜桃av| 国产女主播在线喷水免费视频网站| 少妇被粗大猛烈的视频| av国产久精品久网站免费入址| 国产在线一区二区三区精| 成人免费观看视频高清| 视频区图区小说| 丰满迷人的少妇在线观看| 亚洲三区欧美一区| 亚洲,欧美,日韩| 久久99精品国语久久久| 黄色毛片三级朝国网站| 日韩一区二区视频免费看| 久久精品国产亚洲av高清一级| 各种免费的搞黄视频| 夜夜骑夜夜射夜夜干| 精品亚洲成国产av| 亚洲色图 男人天堂 中文字幕| 久久久久久久久免费视频了| 中文字幕亚洲精品专区| 在线看a的网站| av不卡在线播放| 岛国毛片在线播放| 亚洲情色 制服丝袜| 久久久精品免费免费高清| 国产麻豆69| 九色亚洲精品在线播放| 天堂中文最新版在线下载| 午夜老司机福利剧场| 蜜桃在线观看..| 精品久久蜜臀av无| av线在线观看网站| 日韩av不卡免费在线播放| 97人妻天天添夜夜摸| 国产av精品麻豆| 国产人伦9x9x在线观看 | 性色av一级| 亚洲一区二区三区欧美精品| 在线 av 中文字幕| 久久久亚洲精品成人影院| 午夜免费鲁丝| 美女高潮到喷水免费观看| 亚洲精品av麻豆狂野| 热99久久久久精品小说推荐| 成人漫画全彩无遮挡| 精品国产一区二区三区四区第35| 中文精品一卡2卡3卡4更新| 在线亚洲精品国产二区图片欧美| 亚洲欧美成人综合另类久久久| 日韩中文字幕欧美一区二区 | 久久久国产欧美日韩av| 精品久久蜜臀av无| 色婷婷久久久亚洲欧美| 久久精品人人爽人人爽视色| 亚洲成人手机| 国产精品二区激情视频| 国产精品人妻久久久影院| 亚洲第一av免费看| 日韩伦理黄色片| 99热全是精品| 成人毛片60女人毛片免费| 久久这里只有精品19| 亚洲三区欧美一区| 亚洲视频免费观看视频| 日韩成人av中文字幕在线观看| 国产色婷婷99| 男女边吃奶边做爰视频| 青春草视频在线免费观看| 一边亲一边摸免费视频| 久久久国产欧美日韩av| 又大又黄又爽视频免费| 99久久人妻综合| 欧美精品一区二区免费开放| 一级a爱视频在线免费观看| 一级,二级,三级黄色视频| 2021少妇久久久久久久久久久| 国产一区二区激情短视频 | 亚洲av综合色区一区| 最近中文字幕2019免费版| 午夜免费鲁丝| 欧美老熟妇乱子伦牲交| 伦理电影免费视频| 国产毛片在线视频| 国产精品三级大全| 国产野战对白在线观看| 亚洲综合色惰| 精品99又大又爽又粗少妇毛片| 欧美亚洲 丝袜 人妻 在线| 婷婷色综合www| 午夜福利在线观看免费完整高清在| 卡戴珊不雅视频在线播放| 尾随美女入室| 日韩成人av中文字幕在线观看| 日韩在线高清观看一区二区三区| 欧美日韩亚洲国产一区二区在线观看 | 国语对白做爰xxxⅹ性视频网站| 宅男免费午夜| 久久青草综合色| 亚洲精华国产精华液的使用体验| 日本猛色少妇xxxxx猛交久久| 97在线人人人人妻| 亚洲国产欧美在线一区| 国产日韩欧美亚洲二区| 五月伊人婷婷丁香| 欧美成人午夜精品| 69精品国产乱码久久久| 如何舔出高潮| 久久精品亚洲av国产电影网| 天堂俺去俺来也www色官网| 日韩一本色道免费dvd| 我要看黄色一级片免费的| 久久国产精品大桥未久av| 性色avwww在线观看| 国产一区二区三区综合在线观看| 欧美日韩视频精品一区| 中文字幕制服av| 美女视频免费永久观看网站| 男女无遮挡免费网站观看| kizo精华| 哪个播放器可以免费观看大片| 色视频在线一区二区三区| 欧美变态另类bdsm刘玥| 亚洲国产av新网站| 黄网站色视频无遮挡免费观看| 人体艺术视频欧美日本| 欧美bdsm另类| 纯流量卡能插随身wifi吗| 美女大奶头黄色视频| 2021少妇久久久久久久久久久| av网站免费在线观看视频| 久久国内精品自在自线图片| 婷婷色麻豆天堂久久| 国产精品久久久久久久久免| 街头女战士在线观看网站| 性色av一级| 日韩人妻精品一区2区三区| 国产无遮挡羞羞视频在线观看| 最新中文字幕久久久久| 在线亚洲精品国产二区图片欧美| 欧美日韩视频精品一区| 在线观看人妻少妇| 桃花免费在线播放| 久久精品aⅴ一区二区三区四区 | 亚洲欧美一区二区三区久久| 美女脱内裤让男人舔精品视频| 精品久久久精品久久久| 你懂的网址亚洲精品在线观看| 狂野欧美激情性bbbbbb| 男女免费视频国产| 少妇熟女欧美另类| 国产精品成人在线| 人妻少妇偷人精品九色| 亚洲美女黄色视频免费看| 日韩一区二区三区影片| 精品国产一区二区久久| 精品亚洲乱码少妇综合久久| 午夜久久久在线观看| 亚洲三级黄色毛片| 欧美日韩亚洲国产一区二区在线观看 | 精品酒店卫生间| 人人妻人人添人人爽欧美一区卜| 人体艺术视频欧美日本| 国产1区2区3区精品| 欧美日韩一区二区视频在线观看视频在线| 最近2019中文字幕mv第一页| av网站免费在线观看视频| 老汉色∧v一级毛片| 69精品国产乱码久久久| 国产又爽黄色视频| 午夜激情久久久久久久| 精品人妻熟女毛片av久久网站| 国产在视频线精品| 99国产综合亚洲精品| 免费黄色在线免费观看| 成人手机av| 精品亚洲乱码少妇综合久久| 最近最新中文字幕免费大全7| 久久久精品免费免费高清| 一边亲一边摸免费视频| 久久久久国产网址| 亚洲精品av麻豆狂野| 久久久久久久久久久免费av| 黄色视频在线播放观看不卡| 久久久亚洲精品成人影院| 狂野欧美激情性bbbbbb| 精品国产一区二区久久| tube8黄色片| 高清av免费在线| 不卡视频在线观看欧美| 久久韩国三级中文字幕| 人体艺术视频欧美日本| 青春草国产在线视频| 哪个播放器可以免费观看大片| 飞空精品影院首页| 亚洲伊人久久精品综合| 激情视频va一区二区三区| 日韩熟女老妇一区二区性免费视频| 精品国产乱码久久久久久男人| √禁漫天堂资源中文www| 黄色 视频免费看| freevideosex欧美| 亚洲 欧美一区二区三区| 男人舔女人的私密视频| 日韩人妻精品一区2区三区| 中文字幕av电影在线播放| 亚洲第一av免费看| av免费在线看不卡| 天天躁狠狠躁夜夜躁狠狠躁| 美女主播在线视频| 熟妇人妻不卡中文字幕| 成人手机av| 久久国内精品自在自线图片| 国产免费福利视频在线观看| 各种免费的搞黄视频| 国产成人精品婷婷| 欧美日韩av久久| 人人妻人人爽人人添夜夜欢视频| 亚洲伊人久久精品综合| 亚洲欧美一区二区三区国产| 免费观看av网站的网址| 久久久久久久久久人人人人人人| 亚洲精品av麻豆狂野| 黄色 视频免费看| 香蕉国产在线看| 亚洲av综合色区一区| 免费女性裸体啪啪无遮挡网站| 久久韩国三级中文字幕| 亚洲成av片中文字幕在线观看 | 精品一区二区三区四区五区乱码 | 咕卡用的链子| 亚洲国产精品一区二区三区在线| 国产无遮挡羞羞视频在线观看| 日韩精品有码人妻一区| 女的被弄到高潮叫床怎么办| av又黄又爽大尺度在线免费看| 成人午夜精彩视频在线观看| 国产人伦9x9x在线观看 | 2022亚洲国产成人精品| 亚洲成人一二三区av| 一本—道久久a久久精品蜜桃钙片| 两个人免费观看高清视频| 国产精品女同一区二区软件| 亚洲一码二码三码区别大吗| 黄网站色视频无遮挡免费观看| 亚洲天堂av无毛| 精品一区二区三卡| 搡老乐熟女国产| 日本黄色日本黄色录像| 国产老妇伦熟女老妇高清| 大香蕉久久成人网| 久久久久久人人人人人| 亚洲国产av影院在线观看| 免费在线观看完整版高清| av卡一久久| 欧美日韩成人在线一区二区| 999久久久国产精品视频| 亚洲国产av新网站| 高清视频免费观看一区二区| 中文字幕av电影在线播放| 美女大奶头黄色视频| 国产无遮挡羞羞视频在线观看| 精品一区二区三区四区五区乱码 | 亚洲一区二区三区欧美精品| 丁香六月天网| 2018国产大陆天天弄谢| 国产又色又爽无遮挡免| 亚洲精品国产av成人精品| 日韩av在线免费看完整版不卡| 麻豆精品久久久久久蜜桃| 精品久久久久久电影网| 日韩一区二区视频免费看| av国产久精品久网站免费入址| 成人毛片a级毛片在线播放| 伊人久久国产一区二区| 国产精品免费大片| 国产熟女欧美一区二区| 97人妻天天添夜夜摸| 最近最新中文字幕免费大全7| 考比视频在线观看| 亚洲欧美清纯卡通| 久久鲁丝午夜福利片| 少妇 在线观看| 美女国产视频在线观看|