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

紫外-可見吸收光譜結合高斯多峰擬合技術測定甲基紅酸離解常數

2012-12-12 02:42:34張建華陳玉苗劉兆清徐常威
物理化學學報 2012年5期

張建華 劉 瓊 陳玉苗 劉兆清 徐常威

(廣州大學化學化工學院,廣州510006)

1 Introduction

The acid-base ionization equilibrium exists in aqueous solution of many organic dyes.Acid dissociation constant is a very important parameter to indicate degree of ionization at the different pH values in acidic organic dye solutions.1-3The acidbase ionization equilibrium of organic dyes and the determination of acid dissociation constant in aqueous solution are of great significance for many practical applications and scientific research areas such as acid-base titration,complex formation, solvent extractions,capillary electrophoresis,4chromatographic retention characteristics,5,6potentiometric titration,fabric dyeing and finishing,7environmental monitoring and protection,8-12drug research and development.In particular,drug synthesis, production,purification,formulation,dissolution,absorption, distribution and metabolism processes are closely related with the pKa.13,14Many determination methods for acid dissociation constants of organic dyes have been developed.3,15-17The most commonly used methods are potentiometric titration,5,18-21conductivity,22,23capillary electrophoresis,15capillary electrophoresis mass spectrometry,24nuclear magnetic resonance spectroscopy,25-28liquid chromatography,6,16,29infrared spectroscopy,30,31Raman scattering,32fluorescence spectroscopy,33,34UV-visible spectrophotometry,35-37theoretical calculation38and so on.Among them,spectrophotometric method has high precision and accuracy.The linear relationship between absorbance and concentration of color solution with a certain thickness is determined by Lambert-Beer?s law,which gives a theoretical foundation of the spectrophotometric determination of the pKaof an acidbase indicator.

In 1958,Tobey39determined the pKaof methyl red using singlewavelength spectrophotometry and the pKais 5.02 at 300.35 K.Since then,a lot of work to determine the pKaof organic dyes using single-wavelength spectrophotometry has been reported.3,35-37Single-wavelength spectrophotometry is very suitable to determine the position of absorption peaks of acid and base species which are separated.However,in most cases,the maximum absorption peaks of acid and base species are association with each other in organic dye aqueous solutions.This situation leads to difficulties in data processing and determination of the pKa.3The multi-wavelengths spectrophotometric method has been adopted to determine the value of pKa.4,40,41Target factor analysis,4rank annihilation factor analysis42,43and other methods6,44have been developed to deduce the pKavalues from the multi-wavelengths spectrophotometric data obtained at different pH values.

UV-visible absorption spectra of methyl red were measured at different pH values regulated by a series of acetic acid (HAc)-sodium acetate(NaAc)buffer solutions with different concentrations and characterized by an overlap of a principal peak of acidic specie of methyl red(HMR)at 520-550 nm and a shoulder peak of basic specie of methyl red(MR-)at 425-460 nm.A multi-peaks Gaussian fitting method based on Origin from Microcal Company was used to interpret the spectra in this study.The multi-peaks Gaussian fitting calculation45-47on the overlap peaks gave the integrated absorbance ratioA1/A2,then the pKaof methyl red was obtained.

The studies on organic dye-surfactant interactions in aqueous buffered systems are of great importance in analytical chemistry,pesticide efficiency,pharmaceutical development, fabric dyeing and so on.Therefore,a new procedure based on multi-peaks Gaussian fitting method was firstly performed in dye-surfactant interactions of methyl red with sodium dodecyl sulfate(SDS),cetylammonium bromide(CTAB),which will enrich research methods in this research area.

2 Experimental

2.1 Reagents

Methyl red,anhydrous NaAc,CTAB,SDS were of analytical grade and purchased from Sigma-Aldrich.Ethanol,HAc, HCl were of analytical grade and purchased from Guangzhou chemical reagent factory.All solutions were prepared with distilled water.

(1)All the solutions of methyl red were prepared according to the literature37and the pH values were achieved using a certain concentration of NaAc-HAc buffer solution.

(2)A series of methyl red solutions with SDS were prepared by a concentration range of SDS from 0.001 to 0.018 mol·L-1and the concentration of methyl red was fixed.

(3)Instead of SDS,a series of methyl red solutions with CTAB were prepared by the concentration of CTAB at 0.0001 mol·L-1(below the critical micelle concentration(CMC)of CTAB)and 0.001 mol·L-1(above the CMC of CTAB)and the concentration of methyl red was fixed.

2.2 Apparatus

UV-visible absorption spectra were recorded on Shimadzu UV2550-UV-visible spectrophotometer(Japan)equipped with 10 mm path length quartz cell.Distilled water was used as reference solution.All the spectra were obtained between 320 and 750 nm and the sampling interval was 0.5 nm.The pH values of the solutions were measured by PHSJ-4A-type PH meter furnished with a combined glass electrode(Shanghai Precision Division-Shanghai Lei magnetic)which was pre-calibrated with at least two buffer solutions at pH 4.00 and 10.00.Each pH value was obtained from the average of three measurements.All pH values and spectra were measured at a constant temperature which was controlled by a super-heated water circulating thermostat bath.The measurement data were imported into PC with Microcal Origin 7.0 for data processing.

2.3 Multi-peaks Gaussian fitting method

The UV-visible spectrum data were imported and plotted with Origin 7.5 software.When menu command of Analysis/ Fit Multi-Peaks/Gaussian was selected,a dialog box with the number and the initial half-width estimated default values of the peaks appeared in the current graphics window.After doubleclicking the mouse at 425 and 520 nm of the spectra respectively,Origin 7.5 automatically completed a multi-peaks Gaussian fitting procedure on certain spectra and gave the line-shape parameters of the UV-visible spectra in the result window.

3 Results and discussion

3.1 Principle of multi-peaks Gaussian fitting method

Ionization equilibrium of methyl red in aqueous solution is given as the following equation

The pH value range of color change of methyl red in aqueous solution is well known as 4.4-6.2.When pH values are 4.63, 4.93,5.39,5.68,the UV-visible absorption spectra and their multi-peaks Gaussian fitting results are shown in Fig.1.It is shown that integrated absorbanceA1of base MR-peaks increases and integrated absorbance A2of acid HMR peaks decreases with the increase of pH value.The increase and decrease of the integrated absorbance of the MR-and HMR are objectively due to the change of the relative concentrations of MR-and HMR.

Existence simultaneously of HMR and MR-in methyl red solution results in spectra with two peaks at a certain pH value. Multi-peaks Gaussian fitting on the spectra with two peaks satisfies following equation

where,y0is baseline,λmax1and λmax2are the maximum absorption wavelengths,w1and w2are half peak widths,A1and A2are the integrated absorbances of the two peaks for MR-and HMR.

Multi-peaks Gaussian fitting method based on the spectragives the λmax1,λmax2,w1,w2,A1,A2and these data are listed in Table 1.

Table 1 Results of the multi-peaks Gaussian fitting method on the spectra of methyl red solution at different pH values and 298.15 K

3.2 Relationship between pKaand the relative integrated absorbance

Acid dissociation constant of methyl red is given as following equation

where,ε1and ε2are the molar absorption coefficients of MR-and HMR,respectively.Then,the methyl red absorption spectra(Fig.2)were measured in base(pH=9)and acid(pH=2)conditions respectively.

When the value of pH is 9,there is only MR-in methyl red solution and the[MR-]=C,then

When the value of pH is 2,there is only HMR in methyl red so-lution and[HMR]=C,then

Fig.2 UV-visible absorption spectra of methyl red solution measured in basic(pH=9)and acid(pH=2)solutions

Setting ε=ε1/ε2and A1or A2are obtained with the same concentration of methyl red,then

The value of ε is obtained as 0.56.

From Eqs.(4)-(8),pKacan be calculated from Eq.(10)

The value of pKais obtained from Eq.(10)at different pH values.The Multi-peaks Gaussian Fitting results are shown in Table 2.

3.3 Error analysis for determination of pKa

It is very clear that the values of pKalisted in Table 2 are slightly lower than the values(4.90±0.20)in literature.39,42,48However,the results are within the error range for the values of pKaand prove the reliability of multi-peaks Gaussian fittingmethod for determination of pKa.The concentrations of MR-and HMR in solutions are determined by the curves for absorbance of maximum absorption wavelength at 425 nm for MR-and that at 520 nm for HMR.However the concentrations of MR-and HMR are determined by the methyl red color range which are unreliable because of blue shift of MR-and red shift of HMR absorption peak.18The absorption peak of MR-bluely shifts from 456 to 436 nm and the absorption peak of HMR redly shifts from 531 to 546 nm when the pH value increases from 4.63 to 5.68 in Table 1.The shifts of absorption peaks make the relationship of the concentration and absorbance diverge obviously from the standard curve and cause about ±(3%-5%)systematic error for pKameasurement which is higher than true value.

Table 2 Relationship between Kaand temperature

A multi-peaks Gaussian fitting based on the absorption spectra of methyl red within the color change interval has been made and gives the relative integrated absorbance of MR-and HMR absorption peaks which is used to determine the relative concentration of MR-and HMR instead of using the unreliable standard curves obtained from the acid and base conditions in this study.So the pKameasurement results are more reliable and repeatable.

3.4 Thermodynamics on ionization equilibrium of methyl red

The pKavalues of methyl red were determined by UV-visible absorption spectroscopy with multi-peaks Gaussian fitting method at different temperatures and shown in Table 2.The values of pKadecrease considerably with increase of temperature,which is consistent with the literature.49Generally it is believed that the increase of temperature promotes ionization of organic dyes and leads to decrease of pKa.

A few of thermodynamic models have been developed for acid dissociation equilibrium of oganic dyes.49Here,the‘density’model has been selected.49The standard Gibbs free energy of reaction(ΔrG?)for the ionization equilibrium of methyl red satisfies the following Eq.(12)

where pKa=-lgKa,ρwis water density(kg·m-3),T is the thermodynamic temperature(K),and a-g are model parameters.The relationship between the lgKaand 1/T is shown in Fig.3.The ΔrG?obtained from Eq.(12)with parameters which were obtained from non-linear curve fitting on data of Fig.3 with change of temperature is shown in Fig.4.The ΔrG?increases with increase of temperature and is consistent with the literature.49

3.5 Effect of surfactants on pKa

3.5.1 SDS-methyl red system

The variety of CMC of surfactants can be measured accurately by UV-visible absorption spectrum combined with multipeaks Gaussian fitting method in order to study their aggregation behavior.45To study the effects of aggregation behavior of surfactant on pKa,SDS-methyl red and CTAB-methyl red aque-ous solution systems were selected.Fig.5 shows the relationships between pKaof methyl red and the SDS concentration at different temperatures.The values of pKadecrease considerably with the increase of temperature,indicating that the ionization equilibrium of methyl red moves to the right in the existence of anionic surfactant SDS.When the SDS concentration is lower than the CMC,no significant change of pKavalues is observed.When the SDS concentration is higher than the CMC,pKavalues increase with the increase of SDS concentration.It is undoubted that the pKais sensitive to the CMC and used for determination of the CMC of SDS,for example the CMC of SDS are 8.54×10-3,8.68×10-3,8.82×10-3,9.66×10-3mol·L-1at 303.15,308.15,313.15,318.15 K,respectively, which are consistent with the literature.50

Fig.3 Relationship between lgKaand 1/T

Fig.4 Relationship between ΔrG?and temperature

Fig.5 Relationships between pKaof methyl red and the SDS concentration at different temperatures

Table 3 Spectrum line-shape parameters of methyl red with different concentrations of SDS at 298.15 K and pH=5.40

In order to investigate further effect of SDS on ionization equilibrium of methyl red in solution,a multi-peaks Gaussian fitting method on UV-visible spectra of the series of solutions (pH=5.40)containing different concentrations of SDS from 0.001 to 0.018 mol·L-1gives spectral line-shape parameters such as A1and A2,λmax1and λmax2,w1and w2of MR-and HMR at 298.15 K in Table 3.

The spectral line-shape parameters listed in Table 3 are sensitive to the SDS concentration and changes suddenly at CMC of SDS which are shown in Fig.6 and Fig.7.The red shift of λmax1and blue shift of λmax2with the increase of the SDS concentration are observed and more obvious after formation of the SDS micelle.There is a sudden change at 0.0083 mol·L-1in Fig.6 and Fig.7.At first,the λmax1of MR-absorption peak increases slowly and then increases quickly with the increase of SDS concentration.Sudden change occurs at the CMC of SDS.The change of λmax2is exactly the opposite with λmax1.The w1of MR-absorption peak increases then decreases quickly with the increase of SDS concentration.Sudden change occurs at the CMC of SDS.The change of w2is exactly the opposite with w1.Generally,relationship between spectral line-shape parameters such as A1,A2,λmax1,λmax2,w1,w2and concentration of surfac-tant has been used by us for CMC determination.45

Fig.6 Dependence of the maximum absorption wavelengths of MR-and HMR on the concentration of SDS at 298.15 K

Table 4 Spectrum line-shape parameters and effect of aggregation behavior of CTAB on pKaof methyl red with different concentra

3.5.2 CTAB-methyl red system

CTAB with methyl red system is much more complicated than SDS with methyl red system.The UV-visible absorption spectrum of CTAB with methyl red solution is shown in Fig.8 when CTAB concentration is bellow and above the CMC.Table 4 gives the spectrum line-shape parameters and the effect of aggregation behavior of CTAB on the pKaof methyl red solutions at 303.15 K.

When the CTAB concentration is 0.0001 mol·L-1which is bellow the CMC of CTAB,the UV-visible absorption spectrum of CTAB with methyl red solution is similar with that of methyl red solution.However,the pKavalue of methyl red with 0.0001 mol·L-1CTAB obtained by multi-peaks Gaussian fitting method is 4.67 which is 0.09 lower than that of the methyl red(4.76)at the same temperature.Furthermore,pKavalue decreases slowly with the increase of CTAB concentration.When the CTAB concentration is 0.001 mol·L-1,the pKavalue decreases rapidly to below 3.32 and the color of solution system changes completely from red to yellow at the same time,which indicates that there is mainly MR-species in the solution. These results illustrate that formation of the CTAB micelle can changes the range of changed color of methyl red from 4.4-6.2 to 3.0-5.0 and makes the pKavalue of methyl red decrease remarkably.Generally the significant decrease of pKavalue is attributed to electrical double layer of CTAB micelle which absorbs selectively MR-and makes the ionization equilibrium of methyl red move obviously to the right.

Fig.7 Dependence of the half-widths of MR-and HMR on the concentration of SDS at 298.15 K

Fig.8 UV-visible absorption spectra of methyl red solution with different concentrations of CTAB at 303.15 K

4 Conclusions

In this study,a multi-peaks Gaussian fitting method based on the UV-visible absorption spectra is firstly used to determine the pKavalue of organic dyes such as methyl red.The reliability of the method is adequately proved by an excellent agreement of the measurement results with literature.There are several advantages,such as the easy operational procedure,explicit physical meaning,and the accurate measurement results for the method.

The relative concentration of the MR-and HMR in solution has been determined by the relative integrated absorbance calculated from multi-peaks Gaussian fitting method based on the UV-visible spectra.This method avoids successfully the systematic error of 3%-5%of pKameasurement from the standard curves which is established from the absorbance of MR-maximum absorption wavelength at 425 nm and HMR maximum absorption wavelength at 520 nm versus concentration of MR-and HMR at the methyl red color range,respectively. Moreover,the effect of SDS and CTAB on ionization equilibrium of methyl red also has been studied and some reliable results with the method have been obtained.The spectral lineshape parameters such as A1,A2,λmax1,λmax2,w1and w2of MR-and HMR that obtained from multi-peaks Gaussian fitting method on the UV-visible spectra of SDS-methyl red and CTAB-methyl red solutions have been firstly discovered to be sensitive to aggregation behavior of surfactants SDS and CTAB.The CMC values can be determined by the dependence of the three sets of parameters on the surfactant concentration which should support each other.

(1) Kara,D.;Alkan,M.Spectrochim.Acta A 2000,56,2753.

(2) Niyazi,A.;Yazdanipour,A.;Ramezani,M.Chin.Chem.Lett. 2007,18,989.

(3) Babi?,S.;Horvat,A.J.M.;Pavlovi?,D.M.;Ka?telan-Macan, M.Trends Anal.Chem.2007,26,1043.

(4)Allen,R.I.;Box,K.J.;Comer,J.E.A.;Peake,C.;Tam,K.Y. J.Pharmaceut.Biomed.Anal.1998,17,699.

(5) Beltran,J.L.;Sanli,N.;Fonrodona,G.;Barron,D.;Ozkanb,G.; Barbosa,J.Anal.Chim.Acta 2003,484,253.

(6) Erdemgil,F.Z.;Sanli,S.;Sanli,N.;Ozkan,G.;Barbosa,J.; Guiteras,J.;Beltran,J.L.Talanta 2007,72,489.

(7)Tang,R.C.;Tang,H.;Yang,C.Ind.Eng.Chem.Res.2010,49, 8894.

(8) Simon,E.W.;Beevers,H.New Phytol.1952,51,163.

(9) Adam,R.S.,Jr.Res.Rev.1973,47,1.

(10) Weber,J.B.Adv.Chem.Ser.1972,111,55.

(11) Halling-S?rensen,B.;Nielsen,S.N.;Lanzky,P.F.;Ingerslev,F.; Lutzhoft,H.C.H.;J?rgensen,S.E.Chemosphere 1998,36, 357.

(12) Burns,D.C.;Ellis,D.A.;Li,H.X.;McMurdo,C.J.;Webster, E.Environ.Sci.Technol.2008,42,9283.

(13) Lin,J.H.;Lu,A.Y.Pharmacol.Rev.1997,49,403.

(14) Frey,P.A.;Kokesh,F.O.;Westheimer,F.H.J.Am.Chem.Soc. 1971,93,7266.

(15) Poole,S.K.;Patel,S.;Dehring,K.;Workman,H.;Poole,C.F. J.Chromatogr.A 2004,1037,445.

(16) Hardcastle,J.E.;Jano,I.J.Chromatogr.B 1998,717,39.

(17) Ye,L.;Zhu,Q.Q.;Wu,S.K.Acta Phys.-Chim.Sin.1987,3, 272. [葉 玲,朱琴琴,吳世康.物理化學學報,1987,3,272.]

(18) Qiang,Z.;Adams,C.Water Res.2004,38,2874.

(19) Zhang,W.M.;Yang,Z.D.;Liu,J.;Sun,Z.X.Acta Phys.-Chim. Sin.2010,26,2109.[張衛民,楊振東,劉 嘉,孫中溪.物理化學學報,2010,26,2109.]

(20) Li,L.F.;Hou,W.G.;Jiao,Y.N.;Liu,C.X.Acta Phys.-Chim. Sin.,2004,20,459.[李麗芳,侯萬國,焦燕妮,劉春霞.物理化學學報,2004,20,459.]

(21) Zhang,X.D.;Liu,Y.;Sun,J.Y.;Liu,Q.T.Acta Phys.-Chim. Sin.2000,16,351.[張向東,劉 巖,孫錦玉,劉祁濤.物理化學學報,2000,16,351.]

(22) Kolthoff,I.M.;Chantooni,M.K.J.Phys.Chem.1966,70,856.

(23) Nag,S.;Datta,D.Indian J.Chem.2007,46A,1263.

(24)Wan,H.;Holmen,A.G.;Wang,Y.;Lindberg,W.;Englund,M.; Nagard,M.B.;Thompson,R.A.Rapid Commun.Mass Sp. 2003,17,2639.

(25) Szakacs,Z.;Hagele,G.Talanta 2004,62,819.

(26) Rabenstein,D.L.;Hari,S.P.;Kaerner,A.Anal.Chem.1997,69, 4310.

(27) Rabenstein,D.L.;Sayer,T.L.Anal.Chem.1976,48,1141.

(28) Wang,J.;Rabenstein,D.L.Anal.Chem.2007,79,6799.

(29) Oumada,F.Z.;Rafols,C.;Roses,M.;Bosch,E.J.Pharm.Sci. 2002,91,991.

(30) Lebrón-Paler,A.;Pemberton,J.E.Anal.Chem.2006,78,7649.

(31) Lachenwitzer,A.;Li,N.;Lipkowski,J.J.Electroanal.Chem. 2002,532,85.

(32) Edwards,H.G.M.Spectrochim.Acta A 1989,45,715.

(33) Cagigal,E.;Gonzalez,L.;Alonso,R.M.;Jimenez,R.M. J.Pharmaceut.Biomed.Anal.2001,26,477.

(34) Ferrari,V.;Cutler,D.J.J.Pharmaceut.Sci.1987,76,554.

(35) Cessna,A.J.;Grover,R.J.Agric.Food Chem.1978,26,289.

(36) Foulon,C.;Duhal,N.;Lacroix-Callens,B.;Vaccher,C.;Bonte, J.P.;Goossens,J.F.Eur.J.Pharm.Sci.2007,31,165.

(37) Barbosa,J.;Barron,D.;Jimenez-Lozano,E.;Sanz-Nebot,V. Anal.Chim.Acta 2001,437,309.

(38) Jang,Y.H.;Hwang,S.G.;Chang,S.B.;Ku,J.;Chung,D.S. J.Phys.Chem.A 2009,113,13036.

(39) Tobey,S.W.J.Chem.Educ.1958,35,514.

(40) Tam,K.Y.;Takacs-Novak,T.Anal.Chim.Acta 2001,434,157.

(41) Tam,K.Y.;Hadley,M.;Patterson,W.Talanta 1999,49,539.

(42) Khalafi,L.;Rohani,M.;Afkhami,A.J.Chem.Eng.Data 2008, 53,2389.

(43) Zarei,K.;Atabati,M.;Abdinasab,E.E.J.Anal.Chem.2009,4, 314.

(44) Jimenez-Lozano,E.;Marques,I.;Barron,D.;Beltran,J.L.; Barbosa,J.Anal.Chim.Acta 2002,464,37.

(45) Li,X.G.;Zhang,J.H.;Liu,Z.Q.;Chen,S.;Su,Y.Z.;Xu,C. W.Global J.Phys.Chem.2011,2,34.

(46) Zhang,J.H.;Kong,K.Q.;He,Z.L.;Liu,Z.L.Spectroscopy and Spectral Analysis 2007,27,1412. [張建華,孔凱清,何爭玲,劉自立.光譜學與光譜分析,2007,27,1412.]

(47) Feng,W.S.;Fang,Y.;Xu,J.X.;Fang,C.H.;Jia,Q.J.;Wang, H.H.;Jiang,X.M.Acta Phys.-Chim.Sin.2008,24,497.[馮望生,房 艷,徐繼香,房春暉,賈全杰,王煥華,姜曉明.物理化學學報,2008,24,497.]

(48) Patterson,G.S.J.Chem.Educ.1999,76,395.

(49) Ehlerova,J.;Trevani,L.;Sedlbauer,J.J.Sol.Chem.2008,37, 857.

(50) Mukerjee,P.;Mysels,K.National Standards Reference,Data Series;National Bureau of Standards.US.Government Printing Office:Washington,D.C.,U.S.,1971;Vo1.36,pp 8-71.

主站蜘蛛池模板: 视频二区亚洲精品| 成人毛片在线播放| 日韩小视频在线播放| 日本欧美在线观看| 三区在线视频| AV不卡国产在线观看| 在线看AV天堂| 丁香婷婷在线视频| 亚洲天堂网在线观看视频| 欧美五月婷婷| 国产高清国内精品福利| 国产一区亚洲一区| 精品福利一区二区免费视频| 全部免费特黄特色大片视频| 欧美日韩在线第一页| 精品一区二区三区无码视频无码| 久久久久免费精品国产| 国产精品大白天新婚身材| 真实国产乱子伦高清| 亚洲浓毛av| 亚洲精品第一页不卡| 99精品这里只有精品高清视频| 欧美国产日产一区二区| 精品午夜国产福利观看| 国内精自视频品线一二区| 欧美性天天| 狠狠色成人综合首页| 久久久成年黄色视频| 久久99热这里只有精品免费看| 免费在线色| 免费啪啪网址| 成人字幕网视频在线观看| 狠狠干欧美| 日韩第一页在线| 亚洲综合一区国产精品| 亚洲综合第一区| 国产成人凹凸视频在线| 亚洲一区二区三区麻豆| 久久成人免费| 青青国产在线| 国内精自线i品一区202| 亚洲热线99精品视频| 四虎成人免费毛片| 无码一区18禁| 园内精品自拍视频在线播放| 波多野结衣久久高清免费| 国产精品美女免费视频大全| 亚洲久悠悠色悠在线播放| 黄色污网站在线观看| 欧美视频二区| 婷婷亚洲综合五月天在线| 国产精品视频系列专区| 国产精品专区第一页在线观看| 在线欧美日韩| 午夜欧美理论2019理论| 国产福利大秀91| 日韩精品少妇无码受不了| 精品无码一区二区三区在线视频| 国产亚洲精品在天天在线麻豆 | 亚洲熟女中文字幕男人总站| 91在线播放免费不卡无毒| 日本三级精品| 亚洲动漫h| 久草国产在线观看| 欧美不卡视频一区发布| 一级成人a做片免费| 国产成人精品高清不卡在线| 成人免费网站久久久| 制服无码网站| 久久精品这里只有国产中文精品| 亚洲第七页| 色综合天天娱乐综合网| 国产在线第二页| 日韩亚洲综合在线| 99久久国产综合精品女同| 免费jjzz在在线播放国产| 日韩精品久久久久久久电影蜜臀| 91精品国产丝袜| 亚洲国产成人久久77| 久久精品国产电影| 国产精品欧美亚洲韩国日本不卡| 国产自在线拍|