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

芳香性與環伸縮振動拉曼光譜相關性的理論研究

2014-11-07 04:45:23付景琳等
科技資訊 2014年7期
關鍵詞:相關性

付景琳等

摘 要:本研究采用密度泛函理論B3LYP方法對芳香性分子LnHn(L=C、Si、Ge,n=3~8)進行幾何優化及NICS和拉曼光譜的計算,進而計算其環伸縮振動拉曼光譜頻率與NICS數值間的相關系數,并以C6H6為例進行了環伸縮振動拉曼光譜理論計算值與實驗測定值間一致性的比對研究。研究表明,C6H6的環伸縮振動拉曼光譜理論計算值與實驗測定值相對誤差在1.8%以內,LnHn(L=C、Si、Ge,n=3~8)系列芳香性分子的芳香性大小判別指標NICS理論計算值與其環伸縮振動拉曼光譜頻率理論計算值間存在高度正相關,即芳香性越大,其環伸縮振動頻率也越大。預測可通過測定芳香性物質的環伸縮振動A1′或A1g拉曼光譜頻率的大小來實現對其芳香性大小的實驗定量測定。

關鍵詞:環伸縮振動頻率 芳香性 相關性

中圖分類號:G64 文獻標識碼:A 文章編號:1672-3791(2014)03(a)-0081-03

早在20世紀90年代中期,Schleyer等在反復計算驗證的基礎上,提出核獨立化學位移(Nucleus-Independent Chemical Shifts,NICS)可作為判別物質芳香性大小的一個簡單而有效的探針和判別標準[1~2]。此標準不僅能用作有機物芳香性大小的判據[1],而且能作為無機芳香物[2]、全金屬簇芳香合物體系[3~5]和三明治芳香配合物[6~9]芳香性大小的判據。但是,利用NICS的大小判別物質芳香性的大小,目前僅僅停留在理論計算階段,在實驗上至今還沒有找到測定物質芳香性的有效辦法。研究發現環狀的芳香性有機分子具有高對稱的環伸縮振動拉曼光譜[10]。研究也發現物質芳香性與其共軛性(即環狀共軛π和σ鍵)有關,而這種共軛性也由其芳香性分子中各相鄰原子間的大小相等的共軛作用力大小所決定,相鄰原子間的大小相等的這種作用力在實驗上表現為存在高對稱的環伸縮振動的A1g/A1′拉曼光譜。那么我們完全有理由推測物質的芳香性大小也一定與物質的高對稱環伸縮振動光譜拉曼光譜大小相關,完全可以通過測定芳香性物質的環伸縮振A1g/A1′拉曼光譜頻率的大小來實現對其芳香性大小的實驗定量測定,解決實驗上無法測定物質芳香性大小的這一難題。為了尋找芳香性分子的芳香性大小(NICS)與其高對稱的環伸縮振動A1g/A1′光譜拉曼光譜頻率大小間的相關性,我們在前人對C6H6、C5H5-、C7H7+[11]、Si6H6、Ge6H6[2]等芳香性大小理論計算的基礎上,用密度泛函理論(DFT)對LnHn(L=C、Si、Ge,n=3~8)系列芳香性分子的芳香性大小NICS、環伸縮振動光譜拉曼光譜頻率大小及其二者間的相關性進行理論計算,以期為實驗尋求測試物質芳香性大小有效辦法提供理論參考。

1 計算和實驗測定方法

利用Gaussian09[12]量子化學計算程序,用高精度的量子化學計算方法B3LYP/6-311+g(d,p)對芳香性分子—CnHn(n=3~8)、SinHn(n=3,5~8)、GenHn(n=3,5~8)進行幾何優化、頻率和拉曼光譜的計算(如表1所示),并利用B3LYP-GIAO/6-311+g(d,p)方法對所研究系列分子的NICS值進行了計算,計算結果如表1所示。采用NBO3.1程序對各基態結構的鍵級進行自然鍵軌道(NBO)分析[13]。文中所涉及的分子幾何結構圖及分子軌道圖均采用Guassview5.0繪制而成。我們對苯的NICS(0.0)計算值與Schleyer等的計算結果完全一致[14]。同時,在常溫常壓下,用NX RFT-RAMAN MODULE的拉曼光譜儀,對測定了分析純苯的拉曼光譜值,實驗測定環伸縮振動(A1g)頻率值993.1 cm-1與文獻報道實驗值996.7 cm-1相近[15],理論計算的環伸縮振動(A1g)頻率值為1011.8 cm-1,與實驗測定值的相對誤差在1.8%以內,說明用我們的方法計算芳香性物質拉曼光譜值是非常準確的。

2 結果

為了尋找所研究系列分子(CnHn(n=3~8)、SinHn(n=3,5~8)、GenHn(n=3,5~8))芳香性大小與其環伸縮振動拉曼光譜的相關性,我們對其進行幾何優化、頻率、鍵級、NICS和拉曼光譜的計算(如表1所示)。

從表1可以看出我所計算的分子均為Dnh,NICS的絕對值、環伸縮振動(A1g/A1′)拉曼光譜與MBO鍵級均隨著n值的增大而減小。例如CnHn(n=3~8)的NICS(min)從|-28.6︱、|-13.6︱、|-12.6︳、|-10.5︱、|-9.5︱、︱-9.0︱依次減小,freq(A1g/A1′)從1664.3、1274.1、1143.7、1011.8、874.4、752.1依次減小,MBO鍵級從1.65、1.76、1.28、1.11、0.83、0.70依次減小。從表1還可以看出CnHn(n=3~8)除C4H42-有4個虛頻外,其余均為穩定結構沒有虛頻,而從SinHn、GenHn(n=3,5~8)的NIF可以看出除Si3H3+、Ge3H3+沒有虛頻,其余均有大于或等于1個虛頻,為不穩定結構。由此看出,這三個系列分子的NICS理論計算值與其環伸縮振動(A1g/A1′)拉曼光譜頻率理論計算值大小均呈高度正相關,分子的NICS絕對值越大,其環伸縮振動頻率也越大,預測可通過測定芳香性物質的環伸縮振動(A1g/A1′)拉曼光譜頻率的大小來實現對其芳香性大小的實驗定量測定。

3 討論與分析

C6H6的環伸縮振動A1g正則振動形式如圖1所示,發現其為高對稱環伸縮振動模式,該振動模式對應的拉曼光譜頻率值的大小與其各相鄰碳原子間的大小相等的作用力呈正相關,即C6H6的各相鄰碳原子間的作用力越大,則其A1g正則振動方式大小的拉曼光譜頻率值也越大。苯的芳香性與碳原子間的共軛性(即存在環狀共軛π和σ鍵)有關,為說明苯的共軛性,我們用Guassview5.0繪制其基態構型分子軌道中存在環電流的幾個重要特征的前線分子軌道圖,如圖2所示,從圖2看出C6H6的HOMO-4特征分子軌道形成π環流,使苯具有π芳香性,而HOMO-9、HOMO-14軌道形成σ環流,使苯具有σ芳香性,而這種共軛性大小也由苯分子中各相鄰碳原子間的大小相等的這種作用力大小(主要由MBO表示)所決定,在實驗上必表現為其共軛性越大其高對稱的環伸縮振動A1g拉曼光譜頻率越高,即苯的芳香性大小與其高對稱的環伸縮振動A1g拉曼光譜頻率大小也呈正相關關系,即苯的芳香性越大,其環伸縮振動頻率也越大。endprint

因為C、Si、Ge屬于同一主族,最外層都具有4個電子,形成化合物的性質十分相似,我們經過對三個系列分子的數據進行整理分析,發現CnHn(n=3~8)、SinHn(n=3,5~8)、GenHn(n=3,5~8)符合C6H6所討論的結果,它們都具有環狀的共軛π和σ鍵,芳香性越大,鍵級越大,環伸縮頻率越大,如SinHn(n=3,5~8)從Si8H82+到Si3H3+,︱NICS(min)︳由︱-7.7︳增加到︱-20.9︳,MBO鍵級由0.72增加到1.56,環伸縮振動頻率由277.9增加到604.7。說明我們找到了芳香性分子的芳香性大小(NICS)與其高對稱的環伸縮振動A1g/A1′光譜拉曼光譜頻率大小間相關性的有力依據。

4 結論

CnHn(n=3~8)、SinHn(n=3,5~8)、GenHn(n=3,5~8)三個系列分子的NICS理論計算值大小與其環伸縮振動(A1g/A1′)拉曼光譜頻率理論計算值大小間呈高度正相關,分子的芳香性越大,其環伸縮振動頻率也越大,預測可通過測定芳香性物質的環伸縮振動(A1g/A1′)拉曼光譜頻率的大小來實現對其芳香性大小的實驗定量測定。

參考文獻

[1] Schleyer,P.v R.;Maerker,C.; Dransfeld,A.;Jiao,H. J.;Hommes, N. J.R.V.E.J.Am.Chem.Soc[Z].1996,118:6317-6318.

[2] Schleyer,P.v R.;Jiao,H.-J.; Hommes,N.J.R. v.E.;Malkin,V. G.; Malkina,O.L.J.Am.Chem.Soc[Z].1997,119:12669-12670.

[3] Srinivasu.K.,Swapan K.Ghosh,Das.R,The Royal Society of Chemistry[Z].2012,2:2914-2922.

[4] Jose M.Mercero,Jesus M.Ugalde,J. Am.Chem.Soc[Z].2004,126:3380-3381.

[5] Li M.Y,Yi H.D,Chem.Eur.J.2007, 13:2546-2555.

[6] Liu Z. Z.,Tian W.Q.,Feng J.,K., Zhang G.,Li.W.,Q.,J.Phys.Chem. A[Z],2005,109:5645-5655.

[7] Liu Z.Z.,Tian W.Q.,Feng J.,K.,Zhang G.,Li.W.Q J.Mol.Struc.:THEOCHEM[Z],2006,758(2-3):127-138.

[8] Liu Z.Z.,Tian W.Q.,Feng J.,K.,Zhang G.,Li.W.,Q.Eur.J.Inorg.Chem[Z].2006(14):2808-2818.

[9] Liu Z.Z.Tian W.Q.Feng J.,K., Zhang G.Li.W.Q J.Mol.Struct.: THEOCHEM[Z].2007,809(1-3):171-179.

[10] Chen,Z.F.Wannere,C.S.;Corminboeuf,C.;Puchta, R.;Schleyer, P.v R.,Chem.Rev[Z].2005:105,3842-3888.

[11] Patrick W.Fowler.;Erich Steiner.J.Phys.Chem.1997,101:1409-1413.

[12] M.J.Frisch,G.W.Trucks,H.B.Schlegel,G.E.Scuseria,M.A.Robb,J.R.Cheeseman,G.Scalmani,V.Barone,B.Mennucci,G.A.Petersson,H. Nakatsuji,M.Caricato,X.Li,H.P.Hratchian,A.F.Izmaylov,J.Bloino,G.Zheng,J.L.Sonnenberg,M.Hada,M.Ehara,K.Toyota,R.Fukuda,J. Hasegawa,M.Ishida,T.Nakajima,Y. Honda,O.Kitao, H.Nakai,T.Vreven, J.A.Montgomery,Jr.,J.E.Peralta, F.Ogliaro,M.Bearpark, J. J.Heyd, E.Brothers,K.N.Kudin,V.N. Staroverov,R.Kobayashi,J. Normand, K.Raghavachari,A.Rendell,J.C. Burant, S. S. Iyengar, J. Tomasi, M. Cossi, N. Rega, J. M. Millam, M. Klene, J. E. Knox, J. B. Cross, V. Bakken, C. Adamo, J. Jaramillo, R. Gomperts,R.E.Stratmann,O. Yazyev,A.J.Austin,R.Cammi,C. Pomelli,J.W.Ochterski, R. L.Martin,K.Morokuma,V.G.Zakrzewski,G.A. Voth,P.Salvador,J.J.Dannenberg,S. Dapprich,A.D.Daniels,O.Farkas,J. Foresman,J.V.Ortiz,J.Cioslowski, and D.J.Fox,Gaussian,Inc.,Wallingford CT[Z].2009.

[13] Glendening,E.D.;Carpenter,J.E.;Weinhold,F.NBO,Version3.1,Theoretical Chemistry Institute and Department of Chemistry,University of Wisconsin[Z].1990.

[14] Schleyer,P.vR.;Manoharan,M.;Wang,Z.-X;Kiran,;Jiao,H,;Puchta,R,;Hommes,N.J.R.v.E.Org.Lett[Z].2001,16:2465-2468.

[15] Ma J.Lin.S.L.Hang.R. [J].Physics and Engineering,2013,23(1):29-32.(馬靖,林斯樂,黃蓉[J].物理與工程,2013,23(1):29-32.endprint

因為C、Si、Ge屬于同一主族,最外層都具有4個電子,形成化合物的性質十分相似,我們經過對三個系列分子的數據進行整理分析,發現CnHn(n=3~8)、SinHn(n=3,5~8)、GenHn(n=3,5~8)符合C6H6所討論的結果,它們都具有環狀的共軛π和σ鍵,芳香性越大,鍵級越大,環伸縮頻率越大,如SinHn(n=3,5~8)從Si8H82+到Si3H3+,︱NICS(min)︳由︱-7.7︳增加到︱-20.9︳,MBO鍵級由0.72增加到1.56,環伸縮振動頻率由277.9增加到604.7。說明我們找到了芳香性分子的芳香性大小(NICS)與其高對稱的環伸縮振動A1g/A1′光譜拉曼光譜頻率大小間相關性的有力依據。

4 結論

CnHn(n=3~8)、SinHn(n=3,5~8)、GenHn(n=3,5~8)三個系列分子的NICS理論計算值大小與其環伸縮振動(A1g/A1′)拉曼光譜頻率理論計算值大小間呈高度正相關,分子的芳香性越大,其環伸縮振動頻率也越大,預測可通過測定芳香性物質的環伸縮振動(A1g/A1′)拉曼光譜頻率的大小來實現對其芳香性大小的實驗定量測定。

參考文獻

[1] Schleyer,P.v R.;Maerker,C.; Dransfeld,A.;Jiao,H. J.;Hommes, N. J.R.V.E.J.Am.Chem.Soc[Z].1996,118:6317-6318.

[2] Schleyer,P.v R.;Jiao,H.-J.; Hommes,N.J.R. v.E.;Malkin,V. G.; Malkina,O.L.J.Am.Chem.Soc[Z].1997,119:12669-12670.

[3] Srinivasu.K.,Swapan K.Ghosh,Das.R,The Royal Society of Chemistry[Z].2012,2:2914-2922.

[4] Jose M.Mercero,Jesus M.Ugalde,J. Am.Chem.Soc[Z].2004,126:3380-3381.

[5] Li M.Y,Yi H.D,Chem.Eur.J.2007, 13:2546-2555.

[6] Liu Z. Z.,Tian W.Q.,Feng J.,K., Zhang G.,Li.W.,Q.,J.Phys.Chem. A[Z],2005,109:5645-5655.

[7] Liu Z.Z.,Tian W.Q.,Feng J.,K.,Zhang G.,Li.W.Q J.Mol.Struc.:THEOCHEM[Z],2006,758(2-3):127-138.

[8] Liu Z.Z.,Tian W.Q.,Feng J.,K.,Zhang G.,Li.W.,Q.Eur.J.Inorg.Chem[Z].2006(14):2808-2818.

[9] Liu Z.Z.Tian W.Q.Feng J.,K., Zhang G.Li.W.Q J.Mol.Struct.: THEOCHEM[Z].2007,809(1-3):171-179.

[10] Chen,Z.F.Wannere,C.S.;Corminboeuf,C.;Puchta, R.;Schleyer, P.v R.,Chem.Rev[Z].2005:105,3842-3888.

[11] Patrick W.Fowler.;Erich Steiner.J.Phys.Chem.1997,101:1409-1413.

[12] M.J.Frisch,G.W.Trucks,H.B.Schlegel,G.E.Scuseria,M.A.Robb,J.R.Cheeseman,G.Scalmani,V.Barone,B.Mennucci,G.A.Petersson,H. Nakatsuji,M.Caricato,X.Li,H.P.Hratchian,A.F.Izmaylov,J.Bloino,G.Zheng,J.L.Sonnenberg,M.Hada,M.Ehara,K.Toyota,R.Fukuda,J. Hasegawa,M.Ishida,T.Nakajima,Y. Honda,O.Kitao, H.Nakai,T.Vreven, J.A.Montgomery,Jr.,J.E.Peralta, F.Ogliaro,M.Bearpark, J. J.Heyd, E.Brothers,K.N.Kudin,V.N. Staroverov,R.Kobayashi,J. Normand, K.Raghavachari,A.Rendell,J.C. Burant, S. S. Iyengar, J. Tomasi, M. Cossi, N. Rega, J. M. Millam, M. Klene, J. E. Knox, J. B. Cross, V. Bakken, C. Adamo, J. Jaramillo, R. Gomperts,R.E.Stratmann,O. Yazyev,A.J.Austin,R.Cammi,C. Pomelli,J.W.Ochterski, R. L.Martin,K.Morokuma,V.G.Zakrzewski,G.A. Voth,P.Salvador,J.J.Dannenberg,S. Dapprich,A.D.Daniels,O.Farkas,J. Foresman,J.V.Ortiz,J.Cioslowski, and D.J.Fox,Gaussian,Inc.,Wallingford CT[Z].2009.

[13] Glendening,E.D.;Carpenter,J.E.;Weinhold,F.NBO,Version3.1,Theoretical Chemistry Institute and Department of Chemistry,University of Wisconsin[Z].1990.

[14] Schleyer,P.vR.;Manoharan,M.;Wang,Z.-X;Kiran,;Jiao,H,;Puchta,R,;Hommes,N.J.R.v.E.Org.Lett[Z].2001,16:2465-2468.

[15] Ma J.Lin.S.L.Hang.R. [J].Physics and Engineering,2013,23(1):29-32.(馬靖,林斯樂,黃蓉[J].物理與工程,2013,23(1):29-32.endprint

因為C、Si、Ge屬于同一主族,最外層都具有4個電子,形成化合物的性質十分相似,我們經過對三個系列分子的數據進行整理分析,發現CnHn(n=3~8)、SinHn(n=3,5~8)、GenHn(n=3,5~8)符合C6H6所討論的結果,它們都具有環狀的共軛π和σ鍵,芳香性越大,鍵級越大,環伸縮頻率越大,如SinHn(n=3,5~8)從Si8H82+到Si3H3+,︱NICS(min)︳由︱-7.7︳增加到︱-20.9︳,MBO鍵級由0.72增加到1.56,環伸縮振動頻率由277.9增加到604.7。說明我們找到了芳香性分子的芳香性大小(NICS)與其高對稱的環伸縮振動A1g/A1′光譜拉曼光譜頻率大小間相關性的有力依據。

4 結論

CnHn(n=3~8)、SinHn(n=3,5~8)、GenHn(n=3,5~8)三個系列分子的NICS理論計算值大小與其環伸縮振動(A1g/A1′)拉曼光譜頻率理論計算值大小間呈高度正相關,分子的芳香性越大,其環伸縮振動頻率也越大,預測可通過測定芳香性物質的環伸縮振動(A1g/A1′)拉曼光譜頻率的大小來實現對其芳香性大小的實驗定量測定。

參考文獻

[1] Schleyer,P.v R.;Maerker,C.; Dransfeld,A.;Jiao,H. J.;Hommes, N. J.R.V.E.J.Am.Chem.Soc[Z].1996,118:6317-6318.

[2] Schleyer,P.v R.;Jiao,H.-J.; Hommes,N.J.R. v.E.;Malkin,V. G.; Malkina,O.L.J.Am.Chem.Soc[Z].1997,119:12669-12670.

[3] Srinivasu.K.,Swapan K.Ghosh,Das.R,The Royal Society of Chemistry[Z].2012,2:2914-2922.

[4] Jose M.Mercero,Jesus M.Ugalde,J. Am.Chem.Soc[Z].2004,126:3380-3381.

[5] Li M.Y,Yi H.D,Chem.Eur.J.2007, 13:2546-2555.

[6] Liu Z. Z.,Tian W.Q.,Feng J.,K., Zhang G.,Li.W.,Q.,J.Phys.Chem. A[Z],2005,109:5645-5655.

[7] Liu Z.Z.,Tian W.Q.,Feng J.,K.,Zhang G.,Li.W.Q J.Mol.Struc.:THEOCHEM[Z],2006,758(2-3):127-138.

[8] Liu Z.Z.,Tian W.Q.,Feng J.,K.,Zhang G.,Li.W.,Q.Eur.J.Inorg.Chem[Z].2006(14):2808-2818.

[9] Liu Z.Z.Tian W.Q.Feng J.,K., Zhang G.Li.W.Q J.Mol.Struct.: THEOCHEM[Z].2007,809(1-3):171-179.

[10] Chen,Z.F.Wannere,C.S.;Corminboeuf,C.;Puchta, R.;Schleyer, P.v R.,Chem.Rev[Z].2005:105,3842-3888.

[11] Patrick W.Fowler.;Erich Steiner.J.Phys.Chem.1997,101:1409-1413.

[12] M.J.Frisch,G.W.Trucks,H.B.Schlegel,G.E.Scuseria,M.A.Robb,J.R.Cheeseman,G.Scalmani,V.Barone,B.Mennucci,G.A.Petersson,H. Nakatsuji,M.Caricato,X.Li,H.P.Hratchian,A.F.Izmaylov,J.Bloino,G.Zheng,J.L.Sonnenberg,M.Hada,M.Ehara,K.Toyota,R.Fukuda,J. Hasegawa,M.Ishida,T.Nakajima,Y. Honda,O.Kitao, H.Nakai,T.Vreven, J.A.Montgomery,Jr.,J.E.Peralta, F.Ogliaro,M.Bearpark, J. J.Heyd, E.Brothers,K.N.Kudin,V.N. Staroverov,R.Kobayashi,J. Normand, K.Raghavachari,A.Rendell,J.C. Burant, S. S. Iyengar, J. Tomasi, M. Cossi, N. Rega, J. M. Millam, M. Klene, J. E. Knox, J. B. Cross, V. Bakken, C. Adamo, J. Jaramillo, R. Gomperts,R.E.Stratmann,O. Yazyev,A.J.Austin,R.Cammi,C. Pomelli,J.W.Ochterski, R. L.Martin,K.Morokuma,V.G.Zakrzewski,G.A. Voth,P.Salvador,J.J.Dannenberg,S. Dapprich,A.D.Daniels,O.Farkas,J. Foresman,J.V.Ortiz,J.Cioslowski, and D.J.Fox,Gaussian,Inc.,Wallingford CT[Z].2009.

[13] Glendening,E.D.;Carpenter,J.E.;Weinhold,F.NBO,Version3.1,Theoretical Chemistry Institute and Department of Chemistry,University of Wisconsin[Z].1990.

[14] Schleyer,P.vR.;Manoharan,M.;Wang,Z.-X;Kiran,;Jiao,H,;Puchta,R,;Hommes,N.J.R.v.E.Org.Lett[Z].2001,16:2465-2468.

[15] Ma J.Lin.S.L.Hang.R. [J].Physics and Engineering,2013,23(1):29-32.(馬靖,林斯樂,黃蓉[J].物理與工程,2013,23(1):29-32.endprint

猜你喜歡
相關性
我國創意產業集群與區域經濟發展研究
商情(2016年42期)2016-12-23 14:25:52
淺析財務管理與稅收籌劃的相關性
商情(2016年42期)2016-12-23 13:35:35
醫學教育中基于蛋白質驅動的miRNA與疾病相關性研究
東方教育(2016年4期)2016-12-14 22:15:13
財務管理和稅收籌劃的相關性分析
基于Kronecker信道的MIMO系統通信性能分析
科技視界(2016年21期)2016-10-17 17:37:34
小兒支氣管哮喘與小兒肺炎支原體感染相關性分析
腦梗死與高同型半胱氨酸的相關性研究(2)
腦梗死與高同型半胱氨酸的相關性研究
會計信息質量可靠性與相關性的矛盾與協調
主站蜘蛛池模板: 色视频国产| 亚洲欧美成人网| 国产原创第一页在线观看| 国产成人综合亚洲欧洲色就色| 99视频全部免费| 亚洲香蕉在线| 日本精品视频| 最新国语自产精品视频在| 欧美区国产区| 免费看黄片一区二区三区| 国产JIZzJIzz视频全部免费| 久久久久亚洲av成人网人人软件| 久久久久久久97| 亚洲午夜国产精品无卡| 国产福利在线免费| 久久这里只有精品8| 综合社区亚洲熟妇p| 亚洲国产理论片在线播放| 午夜不卡视频| 国产人成午夜免费看| 国产性生大片免费观看性欧美| 91网址在线播放| 国产亚洲日韩av在线| 美女毛片在线| 午夜免费视频网站| 无码免费视频| 亚欧美国产综合| 色国产视频| 精品国产女同疯狂摩擦2| 一本无码在线观看| 日本一区中文字幕最新在线| 国产极品美女在线观看| 尤物精品视频一区二区三区| 国产午夜精品一区二区三区软件| 四虎成人免费毛片| www.91在线播放| 中文字幕无线码一区| 色哟哟色院91精品网站| 国产打屁股免费区网站| 欧美成人怡春院在线激情| 福利视频99| 国产精品男人的天堂| 国产乱人视频免费观看| 亚洲人成网站日本片| 亚洲乱亚洲乱妇24p| 黄色片中文字幕| 精品撒尿视频一区二区三区| 国产精品亚洲欧美日韩久久| 亚洲欧美色中文字幕| 久久久久青草线综合超碰| 无码久看视频| 日韩不卡高清视频| 在线精品欧美日韩| 一区二区日韩国产精久久| 天天综合网站| 激情综合图区| 好久久免费视频高清| 伊人久久精品无码麻豆精品| 三上悠亚一区二区| 亚洲欧美日韩动漫| 亚洲无码精品在线播放| 伊人五月丁香综合AⅤ| 欧美性精品| 亚洲最大情网站在线观看| 国产成人91精品| 国产精品久久精品| 久久精品视频亚洲| 国产精品原创不卡在线| 真实国产乱子伦高清| 成人在线第一页| 亚洲av无码成人专区| 久久天天躁狠狠躁夜夜躁| 国产一区二区三区在线精品专区| 国产流白浆视频| 久久亚洲精少妇毛片午夜无码| 日韩一级二级三级| 婷婷色一二三区波多野衣| 国产精品任我爽爆在线播放6080| 香蕉蕉亚亚洲aav综合| 毛片一区二区在线看| 热热久久狠狠偷偷色男同| 一本大道AV人久久综合|