亚洲免费av电影一区二区三区,日韩爱爱视频,51精品视频一区二区三区,91视频爱爱,日韩欧美在线播放视频,中文字幕少妇AV,亚洲电影中文字幕,久久久久亚洲av成人网址,久久综合视频网站,国产在线不卡免费播放

        ?

        紫外-可見(jiàn)吸收光譜結(jié)合高斯多峰擬合技術(shù)測(cè)定甲基紅酸離解常數(shù)

        2012-12-12 02:42:34張建華陳玉苗劉兆清徐常威
        物理化學(xué)學(xué)報(bào) 2012年5期
        關(guān)鍵詞:張建華化工學(xué)院物理化學(xué)

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

        (廣州大學(xué)化學(xué)化工學(xué)院,廣州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. [葉 玲,朱琴琴,吳世康.物理化學(xué)學(xué)報(bào),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.[張衛(wèi)民,楊振東,劉 嘉,孫中溪.物理化學(xué)學(xué)報(bào),2010,26,2109.]

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

        (21) Zhang,X.D.;Liu,Y.;Sun,J.Y.;Liu,Q.T.Acta Phys.-Chim. Sin.2000,16,351.[張向東,劉 巖,孫錦玉,劉祁濤.物理化學(xué)學(xué)報(bào),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. [張建華,孔凱清,何爭(zhēng)玲,劉自立.光譜學(xué)與光譜分析,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.[馮望生,房 艷,徐繼香,房春暉,賈全杰,王煥華,姜曉明.物理化學(xué)學(xué)報(bào),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.

        猜你喜歡
        張建華化工學(xué)院物理化學(xué)
        使固態(tài)化學(xué)反應(yīng)100%完成的方法
        張建華靜物素描作品選
        Square grid pattern with direction-selective surface discharges in dielectric barrier discharge
        物理化學(xué)課程教學(xué)改革探索
        云南化工(2021年9期)2021-12-21 07:44:16
        物理化學(xué)課堂教學(xué)改進(jìn)的探索
        云南化工(2021年6期)2021-12-21 07:31:42
        國(guó)家開(kāi)放大學(xué)石油和化工學(xué)院學(xué)習(xí)中心列表
        【鏈接】國(guó)家開(kāi)放大學(xué)石油和化工學(xué)院學(xué)習(xí)中心(第四批)名單
        張建華·書(shū)法篆刻作品欣賞
        Chemical Concepts from Density Functional Theory
        張建華:打造健康產(chǎn)業(yè)新模式
        商周刊(2017年6期)2017-08-22 03:42:47
        亚洲女厕偷拍一区二区| 日韩男女av中文字幕| 国产真实夫妇视频| 亚洲国产精品久久无人区| 少妇bbwbbw高潮| 免费特级毛片| 人妖一区二区三区四区| 久久精品国产亚洲一级二级| 欧美日韩视频无码一区二区三| 免费av网站大全亚洲一区| 激情亚洲综合熟女婷婷| 中文字幕一区二区三区精华液| 91九色老熟女免费资源| 日韩精品免费一区二区中文字幕| 激情偷乱人伦小说视频在线| 天堂中文а√在线| 一区二区三区观看视频在线| 亚洲va在线va天堂va四虎| 好看的欧美熟妇www在线| 黑人老外3p爽粗大免费看视频| 国产特黄1区2区3区4区| 麻豆国产人妻欲求不满| 在线精品无码字幕无码av| 日本免费一二三区在线| 东风日产系列全部车型| 艳妇乳肉豪妇荡乳av无码福利| 国产国产人免费人成免费视频 | 欧美巨大xxxx做受中文字幕| 亚洲国产激情一区二区三区| 欧洲熟妇色xxxx欧美老妇性| 少妇精品揄拍高潮少妇桃花岛| 无码AV午夜福利一区| 国产69精品久久久久9999| 亚洲第一最快av网站| 一区二区三区中文字幕脱狱者 | 久久久久亚洲av片无码| 尤物精品国产亚洲亚洲av麻豆| 丝袜美腿爆炒国产在线观看| 亚洲精品成人网久久久久久| 欧洲熟妇色 欧美| 伊人久久精品亚洲午夜|