999精品在线视频,手机成人午夜在线视频,久久不卡国产精品无码,中日无码在线观看,成人av手机在线观看,日韩精品亚洲一区中文字幕,亚洲av无码人妻,四虎国产在线观看 ?

以甲氧基苯甲酸為配體的兩個雙核銅配合物的合成、晶體結(jié)構(gòu)和理論計(jì)算

2015-06-01 10:45:19于良民夏樹偉閆星辰倪春花
關(guān)鍵詞:海洋大學(xué)理論

張 琦 于良民,2 夏樹偉 李 霞*,,2 閆星辰 倪春花

(1中國海洋大學(xué),海洋化學(xué)理論與工程技術(shù)教育部重點(diǎn)實(shí)驗(yàn)室,青島266100)

(2中國海洋大學(xué),海洋科學(xué)與技術(shù)青島協(xié)同創(chuàng)新中心,青島266100)

以甲氧基苯甲酸為配體的兩個雙核銅配合物的合成、晶體結(jié)構(gòu)和理論計(jì)算

張 琦1于良民1,2夏樹偉1李 霞*,1,2閆星辰1倪春花1

(1中國海洋大學(xué),海洋化學(xué)理論與工程技術(shù)教育部重點(diǎn)實(shí)驗(yàn)室,青島266100)

(2中國海洋大學(xué),海洋科學(xué)與技術(shù)青島協(xié)同創(chuàng)新中心,青島266100)

以2-甲氧基苯甲酸(HL1)、2,3-二甲氧基苯甲酸(HL2)及甲醇為配體,合成了配合物Cu2(L1)4(CH3OH)2(1)和Cu2(L2)4(CH3OH)2(2),并通過紅外、元素分析、X-射線粉末和單晶衍射等研究手段表征了其結(jié)構(gòu)。配合物1屬單斜晶系,空間群P21/n;配合物2屬三斜晶系,空間群P1。2個配合物都具有雙核銅結(jié)構(gòu),由2個銅離子、4個L配體分子和2個甲醇配體分子組成,其中配體L通過雙齒配位模式與銅離子配合。研究了2個配合物的熱穩(wěn)定性,并通過Gaussian 09軟件密度泛函理論B3LYP方法進(jìn)行了理論研究。

甲氧基苯甲酸;雙核銅配合物;晶體結(jié)構(gòu);量子化學(xué)計(jì)算

0 Introduction

In recent years,the design and synthesis of metal complexes have attracted widely interest because of their versatile networktopologies[1-3]and potential applicationsinelectrochemicalsensor,catalysis, magnetism,luminescence,antibacterial activity,biological activity,mental-base drugs[4-11]and so on.As anionic O-donor ligands,carboxylates have strong coordination ability to metal ions,which involve different coordination modes and their complexes exhibit diverse network structures[12-16].Copper is often selected to design and synthesize new complexes with novel structure because of its high probability in forming a variety of mononuclear or multinuclear complexes[17-18].The carboxylates can coordinate with copper cations by monodentate or bidentate bridging coordination mode.

In this paper,two dinuclear copper complexes with methoxybenzoic acids were synthesized and the singlecrystalstructureswereobtainedand characterized.DFT[19-21]studies of the complexes were performed,and frontier molecular orbital energies and components were calculated also.These results would provide theoretical support for further study of the chemical properties and applications.

1 Experimental

1.1 Materials and physical measurement

All reagents used in this work were of analytical grade.The X-ray diffraction data were collected on a Bruker Smart Apex CCD X-ray single-crystal diffractometer.X-ray powder diffraction was determined by a Bruker D8 powder diffractometer at room temperature.FT-IR spectra were recorded on an Avatar 360 FT-IR spectrometer in dry KBr pellets in the range of 4 000~400 cm-1.Elemental analyses were carried out with a model 2400 Perkin-Elmer analyzer.The quantum chemistry calculation was carried out by DFT B3LYP method with 6-31G(d,p)bases for C,H,O atoms and LANL2DZ for copper on Work Station by using Gaussian 09 program.Atom coordinates used in the calculation were derived from crystallographic data.

1.2 Syntheses of the complexes

2-methoxybenzoic acid(0.152 g,1.0 mmol)and KOH(0.056 g,1.0 mmol)were dissolved in 15 mL methanol,and then the solution was stirred for 15 min on a magnetic stirrer.After that Cu(CF3SO3)2·6H2O (0.235 g,0.5 mmol)was added and the solution was stirred for another 15 min.The obtained solution was filtered,and deep green block crystals of complex 1 were gained by slow evaporation of the filtrate at room temperature in the air for two weeks.The crystals were preserved in mother liquor for further detection (Yield:50.5%,based on Cu).Anal.Calcd.for C34H36Cu2O14(%):C 51.32,H 4.56;Found(%):C 51.15, H 4.52.IR(KBr,cm-1):3 488,1 606,1 396,1 247, 1 019,759,669.

The preparation of complex 2 was performed in the same way as complex 1,and the dose of 2,3-dimethoxybenzioc acid was 0.182 g(1.0 mmol).Green block crystals were gained after about two weeks.The crystals were preserved in mother liquor for further detection(Yield:52.4%,based on Cu).Anal.Calcd. for C38H44Cu2O18(%):C 49.84,H 4.84;Found(%):C 49.67,H4.81.IR(KBr,cm-1):3451,1624,1397,1260, 1 061,1 000,780,757.

1.3 Crystallographic data collection and structure determination

A suitable crystal of complex 1 with dimensions of 0.45 mm×0.40 mm×0.37 mm was mounted on a Bruker Smart Apex CCD X-ray single-crystal diffractometer with a Mo Kα radiation(λ=0.071 073 nm)at 298(2)K by using a φ-ω scan mode.A total of 8 472 reflections were collected in the range of 2.66°≤θ≤25.02°(7≤h≤8,13≤k≤13,-25≤l≤17),and 3031 reflections were independent withRint=0.028 8,of which 2 366 reflections were observed with I>2σ(I).

A suitable crystal of complex 2 with dimensions of 0.43 mm×0.31 mm×0.16 mm was mounted on the diffractometer mentioned above.A total of 4 979 reflections were collected in the range of 2.53°≤θ≤25.02°(10≤h≤9,10≤k≤12,-14≤l≤13),and 3422 reflections were independent with Rint=0.019 5,of which 2 787 reflections were observed with I>2σ(I).

Empirical absorption corrections were appliedusing SADABS program[22].The structures of compoundsweresolvedbydirectmethodsusingthe SHELXS-97[23].Positional and thermal parameters of all non-hydrogen atoms were refined by full-matrix least-squaresmethodtoconvergenceusingthe SHELXTL software package[24].The methyl hydrogen atoms were positioned geometrically and refined using a riding model,with dC-H=0.096 nm and Uiso(H)=1.5 Ueq(C).All the other hydrogen atoms were fixed geometrically at calculated distances and allowed to ride on the parent non-hydrogen atoms.Crystal data and structure refinement of the complexes were listed in Table 1.

CCDC:1007777,1;1013765,2.

2 Results and discussion

2.1 Description of crystal structures

Complex1crystallizesinmonoclinic,space group P21/n,and complex 2 crystallizes in triclinic, space group P1.The crystal structures are shown in Fig.1 and Fig.2 ,and the selected bond lengths and bond angles are listed in Table 2.

Complex 1 has dinuclear structure,consisting of two Cucations,four L1ligands and two methanol Ligands molecules,in which the L1ligands are coordinated with copper cations by bidentate bridging coordination mode.The structure of complex 1 is centrosymmetric.The structure of complex 2 is similar as complex 1,and they are shown in Fig.1 (a)and Fig. 2(a).Each copper atom is coordinated with five oxygen atoms derived from four L1or L2ligands and one methanol,Cu1 and Cu1iare connected through the metallic bond.Namely,each Cuis six-coordinated, forming an octahedral configuration.Taking Cu1 of complex 1 as an example,the O5i-Cu1-O2i(88.69(10)°), O4-Cu1-O2i(89.82(10)°),O5i-Cu1-O1(91.31(10)°), O4-Cu1-O1(88.22(10)°)bond angles are all near 90°, indicating that the four atoms(O1,O2i,O4 and O5i) are approximately coplanar.In addition,the O5i-Cu1-O4(169.23(9)°),O2i-Cu1-O1(169.51(9)°),O7-Cu1-Cu1i(175.22(6)°)bond angles are deviated from 180°, indicating that the octahedron is slight distorted.TheCu1 of complex 2 forms a similarly distorted octahedral configuration as complex 1.

Table1 Crystallographic data and structure refinement for complexes 1 and 2

Fig.1 (a)Crystal structure of complex 1(Hydrogen atoms are omitted for clarity;Displacement ellipsoids are drawn at the 30%probability level),(b)1D chain structure of complex 1,(c)Packing diagram of complex 1

The Cu1-Cu1ibond length in complex 1 and complex 2 are 0.257 16(6)nm and 0.261 53(7)nm respectively,which are both shorter than 0.30 nm, indicating the existence of intense coupling effect[25-26]between Cu1 and Cu1i.In complex 1,the Cu1-O1 (0.197 23(19)nm),Cu1-O2i(0.196 72(19)nm),Cu1-O4(0.194 9(2)nm)and Cu1-O5i(0.194 6(2)nm)bond lengths are all shorter than that of Cu1-O7(0.213 87(19) nm),and in the complex 2,Cu1-O1(0.196 9(2)nm), Cu1-O2i(0.195 7(2)nm),Cu1-O5(0.196 7(2)nm), Cu1-O6i(0.196 6(2)nm)bond lengths are all shorter than that of Cu1-O9(0.219 7(2)nm),indicating that the oxygen atoms of L1and L2ligands have better coordinationabilitythantheoxygenatomsfrom methanol molecules.This might result from sterichindrance.In addition,the Cu1-O9 bond of complex 2 is longer than the Cu1-O7 bond of complex 1, suggesting that the steric hindrance effect of L2ligand is stronger than L1ligand.

Table2 Selected bond lengths(nm)and angles(°)for complexes 1 and 2

Hydrogen bonds exist in complex 1 and complex 2 leading to a one-dimensional structure,as shown in Fig.1 (b)and Fig.2 (b).In complex 1,there exist O7-H7…O1iiand O7-H7…O3iihydrogen bonds,and the symmetry code is-x,1-y,1-z.The distances between the hydrogen atom and oxygen atom are 0.220 7 and 0.255 6 nm,and the angles are 158.08°and 128.85° respectively.While in complex 2,there only exists O9-H9…O3iihydrogen bond(Symmetry code:2-x,1-y,1-z).The distance between the hydrogen atom and oxygen atom is 0.209 4 nm,and the angle is 169.68°. The differences are derived from the steric hindrance effect of methoxy group in meta-position.The steric hindrance effect of L2Ligand also leads to the deflectionofbenzenering,sothecomplex2 crystallizes in triclinic rather than monoclinic system.

There exist C-H…π interactions[27-28]in complex 1,and the distance between the hydrogen atom and the center of the phenyl ring is 0.346 7 nm,and the angle is 134.89°.A 3D network constructed by the CH…π interactions and hydrogen bond is shown in Fig.1 (c).In complex 2,there exist the C-H…O bond interaction.The distances between the hydrogen atom and oxygen atom are 0.270 3 and 0.261 4 nm,and the angles are 116.70°,136.31°respectively.A 3D network formed by the O-H…O and C-H…O hydrogen bond interactions is shown in Fig.2 (c).

2.2 IR spectroscopy

In the spectrum of complex 1,a strong wide absorption peak is observed at 3 488 cm-1,which can be attributed to the stretching vibration band of OH in CH3OH.The stretching vibration band of C=O in HL1at 1 693 cm-1disappears and appears at 1 606 cm-1after coordination,and the red shift of νC=Oindicates that the carbonyl group of L1ligand coordinates with Cucations.

In complex 2,a strong wide absorption peak is observed at 3 451 cm-1,which can be attributed to the stretching vibration band of OH in CH3OH.The stretching vibration band of C=O in HL2at 1 685 cm-1disappear and appears at 1 624 cm-1after coordination, and the red shift of νC=Oindicates that the carbonyl group of L2ligand coordinates with Cucations.The resultsareconsistentwiththecrystalstructure analysis.

2.3 X-ray powder diffraction

The simulated and experimental X-ray powder diffraction patterns of complex 1 and complex 2 are shown in Fig.3 .All peaks presented in the measured curves approximately match the simulated curves generated from single-crystal diffraction data,which confirms the phase purity of the complexes.

Fig.3 Comparison of simulated X-ray powder diffraction patterns(upper)and the measured values(lower)of complex 1(a)and complex 2(b)

2.4 Thermogravimetric analysis

Toinvestigatethermalstabilityofthetwo complexes,TG analyses were carried out with a Netzsch STA449 F3 instrument at a heating rate of 10℃·min-1from 30 to 800℃under nitrogen(Flow rate= 50 mL·min-1).

As shown in Fig.4 ,complex 1 displays three weight loss stages.The weight loss of 2.8%(Calcd. 8.1%)in the temperature range of 120~220℃corresponds to partial release of the coordinated methanol molecule.The quick weight loss of 34.2% (Calcd.43.3%)during the temperature range of 220~ 300℃couldbeattributedtothereleaseof coordinated methanol molecule and partial of L1ligand.The third step occurs between 300~800℃, which corresponds to the decomposition of residual ligands,and the final weight loss is 78.8%(Calcd. 84.0%).

Fig.4 TG curves of complex 1 and 2

In comparison with complex 1,complex 2 shows two weight loss stages.The weight loss of 1.6% (Calcd.7.0%)in the temperature range of 120~200℃might be caused by partial release of the coordinated methanol molecule.The second step occurs between 200~800℃,corresponding to release of residual methanol and all the ligands,and the final weight loss is 81.7%(Calcd.86.1%).

2.5 Quantum chemistry calculation

The natural atomic charges of the two complexes are listed in Table 3.The natural charges of Cuwere decreased to 0.77 or 0.78,suggesting that part of theelectrondensityhastransferredfromthe coordinated O atoms to the copper ions.The negative charges mostly concentrate on coordinated O atoms, which indicate O atoms are prone to coordinate with Cucations.The calculation results are in accordancewith the determination.

Some frontier molecular orbital energies and components of the two complexes are listed in Table 4,and views of frontier molecular orbital are shown in Fig.5 and Fig.6 .The total energy of complex 1 and complex 2 are-2 761.443 10 and-3 220.144 25 hartree,respectively.The HOMO-1,HOMO,LUMO, LUMO+1 orbital energies of the two complexes are -0.204 83,-0.166 89,-0.112 59,-0.013 31 a.u.and -0.201 09,-0.175 21,-0.123 61,-0.014 34 a.u..The energy gap between HOMO and LUMO of complex 1 and 2 are 0.054 3 and 0.051 6 a.u.respectively.

The orbital components of HOMO and LUMO are mainly distributed on central ion Cuand ligating O atom of L ligands,and show quite well symmetry as shown in Fig.5 and Fig.6 .The HOMO and LUMO components of Cucations in complex 1 are about 49.5%,42.6%,and in complex 2 are about 51.3%, 49.7%.Apart from the cations,the frontier molecular orbital evenly distributes on the ligands.The HOMO and LUMO components of coordinated O atoms from L1ligand are about 5.3%and 3.6%averagely,and that of L2ligand are about 5.2%and 4.6%averagely, while the orbital components of coordinated O atoms from methanol in the two complexes are both 0%. This indicates that although the methanol exists in the unit cell of the coordination compounds,it does not donate electrons to central ion Cuand just acts as an interstitial molecule.

Table4 Selected frontier molecular orbital energy (a.u.)and component(%)of complexes 1 and 2

Fig.5 View of frontier molecular orbital of complex 1

Fig.6 View of frontier molecular orbital of complex 2

The dipole moments of complex 1 and complex 2 are 0.011 4D,0.012 4D respectively,which due to the high symmetry of complexes.The calculation results identifycentrosymmetricstructuresofthetwo complexes.

3 Conclusions

Reference:

[1]Guo F,Wang F,Yang H,et al.Inorg.Chem.,2012,51(18): 9677-9682

[2]Xu J K,Sun X C,Ju C X,et al.J.Coord.Chem.,2013,66 (14):2541-2548

[3]Lin J D,Jia C C,Li Z H,et al.Inorg.Chem.Commun., 2009,12(6):558-562

[4]TANG Si-Ping(唐斯萍),XU Zhi-Feng(許志鋒),FENG Yong-Lan(馮泳蘭),et al.Chinese J.Inorg.Chem.(無機(jī)化學(xué)學(xué)報(bào)),2013,29(12):2683-2687

[5]Wang X L,Lin H Y,Liu G C,et al.J.Organomet.Chem., 2008,693(16):2767-2774

[6]Qu H,Qiu L,Leng X K,et al.Inorg.Chem.Commun.,2011, 14(9):1347-1351

[7]Qin L,Liu L W,Du X,et al.Transition Met.Chem.,2013, 38(1):85-91

[8]Ma L F,Han M L,Qin J H,et al.Inorg.Chem.,2012,51 (17):9431-9442

[9]LI Hai-Hua(李海華),ZHOU Xiao-Xia(周曉霞),YOU Zhong-Lu(由忠錄).Chinese J.Inorg.Chem.(無機(jī)化學(xué)學(xué)報(bào)), 2013,29(3):649-653

[10]Mostafa M H K,Eglal R S,Eman H I,et al.Chinese J. Inorg.Chem.(無機(jī)化學(xué)學(xué)報(bào)),2013,29(9):1969-1978

[11]Huxford R C,Rocca J D,Lin W.Curr.Opin.Chem.Biol., 2010,14(2):262-268

[12]Haribabu P,Patil Y P,Reddy K H,et al.Transition Met. Chem.,2011,36(8):867-874

[13]WU Gang(吳剛),WANG Xiao-Feng(王小鋒),XIAN Hua(鮮華),et al.Chinese J.Inorg.Chem.(無機(jī)化學(xué)學(xué)報(bào)),2010,26 (7):1315-1318

[14]Ke C H,Lin G R,Kuo B C,et al.Cryst.Growth Des.,2012, 12(7):3758-3765

[15]TANG Si-Ping(唐斯萍),ZHANG Shao-Hua(張少華),YANG Ying-Qun(楊穎群).Chinese J.Inorg.Chem.(無機(jī)化學(xué)學(xué)報(bào)),2013,29(11):2465-2469

[16]Yang W B,Lin X,Blake A J,et al.Inorg.Chem.,2009,48 (23):11067-11078

[17]Soltani B,Sadr M H,Engle J T,et al.Transition Met.Chem., 2012,37(8):687-694

[18]Lin H Y,Mu B,Wang X L,et al.J.Organomet.Chem., 2012,702:36-44

[19]SHI Zhi-Qiang(石志強(qiáng)),JI Ning-Ning(季寧寧),HE Guo-Fang(何國芳),et al.Chinese J.Inorg.Chem.(無機(jī)化學(xué)學(xué)報(bào)),2012,28(6):1279-1285

[20]He H S,Zhong Y H,Si L P,et al.Inorg.Chim.Acta,2011, 378(1):30-35

[21]Jana M S,Pramanik A K,Sarkar D,et al.Polyhedron,2014, 81:66-73

[22]Sheldrick G M.SADABS,Program for Empirical Absorption Correction of Area Detector Data,University of G?ttingen, Germany,1996.

[23]Sheldrick G M.SHELXS-97,Program for the Solution of Crystal Structures,University of G?ttingen,Germany,1997.

[24]Sheldrick G M.SHELXL-97,Program for the Refinement of Crystal Structures,University of G?ttingen,Germany,1997.

[25]Zhao L L,Luo X Z,Xu L,et al.Inorg.Chem.Commun., 2010,13(4):554-557

[26]Yang Q,Zhang X F,Zhao J P,et al.Cryst.Growth Des., 2011,11(7):2839-2845

[27]Kamiński R,Kowalski J,Mames I,et al.Eur.J.Inorg. Chem.,2011,2011(4):479-488

[28]Santana A M,Ferreira J G,Moro A C,et al.Inorg.Chem. Commun.,2011,14(1):83-86

Syntheses,Crystal Structures and Theoretical Calculation of Two Dinuclear CopperComplexes with Methoxybenzoic Acids Ligands

ZHANG Qi1YU Liang-Min1,2XIA Shu-Wei1LI Xia*,1,2YAN Xing-Chen1NI Chun-Hua1

(1Key Laboratory of Marine Chemistry Theory and Technology,Ministry of Education, Ocean University of China,Qingdao,Shandong 266100,China)
(2Qingdao colobrative innovation center of marine science and technology, Ocean University of China,Qingdao,Shandong 266100,China)

Two complexes,Cu2(L1)4(CH3OH)2(1)and Cu2(L2)4(CH3OH)2(2)(HL1=2-methoxybenzoic acid,HL2=2,3-dimethoxybenzoic acid)were synthesized and characterized by IR spectorscopy,elemental analysis,X-ray powder and X-ray single-crystal diffraction.Complex 1 crystallizes in monoclinic,space group P21/n;Complex 2 crystallizes in triclinic,space group P1.Both of the complexes have dinuclear structure,consisting of two Cucations,four L ligands and two methanol ligands molecules,in which the L ligands are coordinated with copper cations by bidentate bridging coordination mode.The thermal stability of the complexes was investigated,and theoretical study of the complexes was carried out by Density Functional Theory(DFT)B3LYP method using Gaussian 09 program.CCDC:1007777,1;1013765,2.

methoxybenzoic acid;dinuclear coppercomplex;crystal structure;quantum chemistry calculation

O614.121

A

1001-4861(2015)03-0585-09

10.11862/CJIC.2015.086

2014-10-13。收修改稿日期:2014-12-23。

國家自然科學(xué)基金(No.51003099、51102219),中國海洋大學(xué)青年教師專項(xiàng)基金(No.201013017、201113024),國家科技支撐計(jì)劃(No.2012BAB15B02)資助項(xiàng)目。

*通訊聯(lián)系人。E-mail:xiali@ouc.edu.cn

猜你喜歡
海洋大學(xué)理論
堅(jiān)持理論創(chuàng)新
神秘的混沌理論
理論創(chuàng)新 引領(lǐng)百年
中國海洋大學(xué)作品選登
相關(guān)于撓理論的Baer模
中國海洋大學(xué) 自主招生,讓我同時(shí)被兩所211大學(xué)錄取
Multilingual Mix in Women Fashion Industry Advertising A Comparison between Hong Kong and China
?? ??? ???? ????
理論宣講如何答疑解惑
La communication sino-fran?aise
主站蜘蛛池模板: 男女性午夜福利网站| 国产欧美又粗又猛又爽老| 亚洲欧美日韩中文字幕在线一区| 五月天丁香婷婷综合久久| 久久影院一区二区h| 亚洲男人天堂网址| 亚洲最猛黑人xxxx黑人猛交| 九色免费视频| 2020国产精品视频| 日韩毛片基地| 久久久久国产一区二区| 成人午夜网址| 日本欧美视频在线观看| AV在线天堂进入| 99精品热视频这里只有精品7| 蝴蝶伊人久久中文娱乐网| 国产精品原创不卡在线| 免费在线视频a| 97视频在线观看免费视频| 午夜成人在线视频| 免费不卡在线观看av| 亚洲毛片网站| 67194在线午夜亚洲 | 亚洲国产日韩在线成人蜜芽| 二级特黄绝大片免费视频大片| 国产又爽又黄无遮挡免费观看 | 国产精品第一区在线观看| 久久久久人妻精品一区三寸蜜桃| 91精品人妻互换| AV不卡无码免费一区二区三区| 狠狠色噜噜狠狠狠狠色综合久| 亚洲综合18p| 这里只有精品在线播放| 国产69精品久久久久孕妇大杂乱 | 秋霞一区二区三区| 另类专区亚洲| 一级毛片在线播放免费观看| 日韩美毛片| 亚洲第一成人在线| 久久久亚洲色| 国产美女精品一区二区| 在线色国产| 99re视频在线| 国产精品美乳| 亚洲AⅤ永久无码精品毛片| 99re精彩视频| 国产精品hd在线播放| a毛片在线免费观看| 亚洲欧美日韩精品专区| 女人18毛片一级毛片在线| 一级黄色网站在线免费看| 国产原创演绎剧情有字幕的| 香蕉综合在线视频91| 国产精品刺激对白在线| 久久国产黑丝袜视频| 刘亦菲一区二区在线观看| 亚洲人成在线免费观看| 国产特级毛片| 日韩av在线直播| 国产成人高清在线精品| 亚洲国产欧美目韩成人综合| 在线精品视频成人网| 久久国产精品影院| 亚洲日韩精品伊甸| 超碰精品无码一区二区| 亚洲欧美日韩动漫| 亚洲AV无码久久精品色欲| 嫩草在线视频| 国产特一级毛片| 亚洲精品老司机| 久久久久久久久亚洲精品| 在线毛片网站| 久久久久久久久亚洲精品| 午夜在线不卡| 国产午夜人做人免费视频| a毛片在线免费观看| 中国国产A一级毛片| 国产丝袜啪啪| 久久黄色影院| 欧美一级特黄aaaaaa在线看片| 尤物成AV人片在线观看| 国产主播一区二区三区|