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

        ?

        一種評(píng)估陽(yáng)離子交換吸附活化能的方法

        2017-11-01 18:11:07李嫣婷劉新敏丁武泉修維寧唐玲玲
        物理化學(xué)學(xué)報(bào) 2017年10期
        關(guān)鍵詞:重慶市特異性重慶

        李嫣婷 劉新敏,* 田 銳 丁武泉 修維寧唐玲玲 張 靜 李 航

        (1西南大學(xué)資源環(huán)境學(xué)院,土壤多尺度界面過(guò)程與調(diào)控重慶市重點(diǎn)實(shí)驗(yàn)室,重慶 400715;2農(nóng)業(yè)工程研究所,重慶市農(nóng)業(yè)科學(xué)院,重慶 401329)

        一種評(píng)估陽(yáng)離子交換吸附活化能的方法

        李嫣婷1劉新敏1,*田 銳1丁武泉1修維寧2唐玲玲1張 靜1李 航1

        (1西南大學(xué)資源環(huán)境學(xué)院,土壤多尺度界面過(guò)程與調(diào)控重慶市重點(diǎn)實(shí)驗(yàn)室,重慶 400715;2農(nóng)業(yè)工程研究所,重慶市農(nóng)業(yè)科學(xué)院,重慶 401329)

        離子和帶電粒子之間的相互作用受納米/膠體顆粒表面上離子吸附控制,同時(shí)它決定著顆粒的微觀性質(zhì)、界面反應(yīng)過(guò)程以及顆粒間相互作用?,F(xiàn)有的離子界面反應(yīng)理論沒(méi)有考慮到離子特異性對(duì)離子-表面相互作用的影響。本研究建立了帶電顆粒表面離子吸附模型。結(jié)果表明蒙脫石膠體表面離子吸附平衡存在著強(qiáng)烈的離子特異性效應(yīng)。離子和表面之間的吸附能隨離子濃度的增加而降低,隨著粒子電荷密度(絕對(duì)值)升高而增加。根據(jù)建立的模型還可預(yù)測(cè)吸附反離子在擴(kuò)散層中的位置,進(jìn)而可計(jì)算出離子的吸附能與活化能?;罨苤苯記Q定了離子的吸附飽和度, 且不同的離子遵循相同的規(guī)律。本研究表明建立的離子-表面相互作用模型在固/液界面反應(yīng)中具有普適性的應(yīng)用。

        離子特異性效應(yīng);離子吸附;帶電表面;活化能

        1 Introduction

        The interactions between ion and charged particle determine the surface properties of the particles, ion adsorption,interactions between particles, colloidal suspension stability and settlement etc. Ion adsorption obeys Boltzmann distribution law at the particle/water interface, suggesting the interaction between ion and the surface is restricted by concentration gradient and electric potential gradient1. The Boltzmann equation is a basic equation in electrical double layer theory, and it is therefore a very important approach to study the ion-surface interactions.

        The classical double layer theory only considered the Coulomb interaction between ion and surface. Thus the interaction between particles and ions with equal valence causes the same interfacial reaction phenomena. This obviously does not match a large number of experimental phenomena2?6.Some studies indicated that ion volume7, hydration8?10and dispersion11etc. also strongly affect the ion-surface interactions.These factors may be the origin of ion specificity or Hofmeister effects, which were discovered through protein aggregation experiments during the 1880s and 1890 s12. Hofmeister effects affect the ion-particle interactions, and then influence many of the micro processes and macro phenomena, such as the adsorption/desorption of ions on charged surface, surface properties of charged particles, particle-particle interactions, the stability of soil aggregates, the structure of the electrical double layer, and so on. But the physical mechanism of the Hofmeister effects is not clear, for more than one hundred years, many researchers demonstrated that this difference is caused only by the radius of the ion or the hydration radius of the ion13. Ion dispersion force or quantum fluctuation is one of origins of Hofmeister effects at high electrolyte concentration or for low charge density of colloidal particles. The changes of Hofmeister effects with different electrolyte concentrations5,14,15, however, cannot be explained by the ion dispersion. On the deep discussion of the driving force of ion-surface interactions, we further find that a strong electric field near the surface exists due to surface charges of nano/colloidal particles16?26. Induced ionic dipole moments in such a strong electric field are different from the softness of ionic electron cloud, and this makes the Hofmeister effects on ion exchange adsorption5.

        The activation energy has been considered as a crucial parameter to characterize the ion-surface interactions27.Generally, the activation energy of ion exchange reaction was estimated using the classical Derjaguin-Landau-Verwey-Overbeek (DLVO) theory, which has not taken Hofmeister effects into account, although Hofmeister effects may play a considerable role28. In this study, a model of ion-surface interactions considering the Hofmeister effects was established in 2 : 1 electrolyte (MgCl2, CaCl2, CuCl2) solutions. Based on the model, the activation energy of ion exchange adsorption at the solid/liquid interface would be estimated. Ion exchange adsorption is a useful approach to study interfacial reactions,and the study of Hofmeister effects on ion exchange adsorption is also helpful to reveal the mechanism of ion-surface interactions.

        2 Theory

        The adsorption capacity of ions in the diffuse layer can be expressed as:

        where S is the specific surface area of the particles; κ is Debye parameter, 1/κ is the thickness of the electrical double layer; l is the distance between surface and the diffusion position of exchangeable counterion. l = 1/κ indicates that the ion does not diffuse into the electrical double layer, and the adsorption capacity is 0. On the contrary, l = 0 indicates that the ion diffuses to the surface of the particles. At this point, the adsorption capacity reached the maximum, i.e. saturated adsorption capacity. The fi(x) is the concentration distribution in the diffuse layer, which obeys the Boltzmann distribution considering the ion-ion interactions in bulk solution29:

        where aiis the activity of ion i species in bulk solution; wi(x) is the adsorption energy of ion i species at the x position in the diffuse layer, R is gas constant, T is absolute temperature.

        Taking the additional energies arising from the Hofmeister effects into account, the ion-surface interaction energy can be expressed as30:

        where Ziis the valence of i ion species; F is Faraday constant;φ(x) is the potential at the x position in the diffuse layer; wiadd(x)is the additional energy arising from the Hofmeister effects.

        Since the Poisson-Boltzmann equation with Hofmeister effects is considered to have no analytical solution, and we consider the correction of the ion specificity on the analytical solution of the classical nonlinear Poisson-Boltzmann equation.In 2 : 1 electrolyte system (2 : 1 electrolyte refers to the two valence cations and monovalent anions form a chemical electrolyte, e.g. CaCl2), the solution of the Poisson-Boltzmann equation taking ion specificity into consideration is31:

        where m is a function of mean interaction energy between ion and surface in 2 : 1 electrolyte. When x = l, the actual

        Equations (5)?(8) show that the activation energy is a function of Niand NiTwhich depends on the values of m and l.The Ni/NiTrepresents the adsorption saturation.

        3 Materials and methods

        Purified montmorillonite was purchased from Wuhuatianbao Mineral Resources Co., Lid, in Inner Mongolia in China. The average particle size and specific surface area are determined to be 1480 nm and 731 m2·g?1, respectively.

        Pretreatment of samples: the Na-saturated montmorillonite was prepared by equilibrating with 0.5 mol·L?1NaCl (AR)solution for five times (mass ratio of montmorillonite and water is 1 : 10), the adsorbed cations were saturated by Na+. After repeatedly washing the Na+-equilibrated samples with pure water for two times, and then the Na-saturated montmorillonite was dried at a temperature of 298 K, and sieved through a 0.25 mm sieve.

        Exchange adsorption: approximate 1 g Na-saturated montmorillonite was weighed and transferred into the triangular bottles, 30 mL 0.05, 0.1, 0.2 and 0.3 mol·L?1MgCl2(AR) solutions were added into each group, respectively (the low concentration was selected to test the ion adsorption saturation). The suspensions in the triangular bottles were allowed to equilibrate for 12 h at 298 K, and then centrifugal separation, the supernatant was then collected in a 200 mL volumetric flask. To make the cation exchange completely, the experiment was repeated four times, and each supernatant was preserved in the same corresponding volumetric flask, and then diluted with deionized water to 200 mL.

        The same method was applied to the Ca/Na and Cu/Na exchange adsorption experiment. The concentration of Na+was determined by flame photometer. The experimental method was proposed by Low32.

        4 Results and discussion

        In the experiment, the supernatant was collected to make the exchange adsorption more fully through adding adsorption ion repeatedly (Mg2+, Ca2+and Cu2+). Since the adsorbed ions were added repeatedly, the adsorption capacity could not be calculated by the difference between the added and remaining amount. Considering montmorillonite is constantly charged mineral, the chemical adsorption could be neglected. According to the charge balance principle, the adsorption capacity of Na-montmorillonite was calculated by the desorption of Na+.

        The data in Table 1 showed that the adsorption capacity of the three ions increased with increasing ion concentrations.Obviously, there was a significant difference among the adsorption capacities for the three species of bivalent ions,suggesting Hofmeister effects exert influence on ion exchange adsorption equilibrium. Due to the constant charge of the montmorillonite surface, the chemical adsorption could be ignored, and the ion adsorption specificity is equal to the ion exchange specificity. In addition, the adsorption in montmorillonite is not controlled by the Stern layer. The experimental results could not be explained by the classical theory which only takes Coulomb interaction into account.

        Based on the experimental process, the true equilibrium adsorption of montmorillonite should be a constant in different electrolyte solutions, and this value is the total surface charge number. By comparing the equilibrium adsorption capacity of montmorillonite under four different concentrations (Table 1),the adsorption capacity of the three ion species can be found:Mg2+< Ca2+< Cu2+, indicating strong Hofmeister effects on ion adsorption at montmorillonite surface exist on the one hand.On the other hand, the experimental adsorption capacity suggests that the achieved ion exchange equilibrium is not in a real state of equilibrium, but a metastable equilibrium. At this time, the exchange reaction rate is so slow that it cannot be determined during the experiment. This may be result from the fact that the adsorption energy and the exchange energy are different and the former is higher than the latter. Here we suppose the adsorption capacity of Cu2+at 0.3 mol·L?1reach approximately the real equilibrium, thus the surface charge of montmorillonite in this study is 1.2 mmol·g?1, this value is approximately equal to the experimental result in the previous study30. The maximum adsorption capacity of the bivalent ions is 0.6 mmol·g?1.

        Table 1 Ionic adsorption capacity on montmorillonite surface in different counterion solutions.

        As long as the surface charge number of particle was obtained, the relationship between the maximum adsorption energy (negative value) and the ion concentration could be described by Eq.(7) (Fig.1). Only when the adsorption energy reaches the maximum value, the counterion adsorption could be saturation.

        Fig.1 shows that the adsorption energy between ion and surface decreases with increasing ion concentration. The larger the surface charge density is (e.g. 0.3168 C·m?2), the higher the adsorption energy will be (the dot line in Fig.1). Obtaining the total saturated adsorption energy, the relationship between the adsorption saturation of ions and the diffusion distance in the double layer can be quantitatively calculated by Eq.(5), as shown in Fig.2. When κl is 0, the adsorbed ions reach the surface and the ion adsorption reaches saturation (100%); when κl is 1, the ions are not adsorbed and the ion adsorption is zero.Simultaneously, ion concentration also affects the adsorption saturation. The higher the ion concentration is, the larger the ion adsorption will be. According to the Eq.(5) or Fig.2, the position could be estimated when the ions are adsorbed in the diffuse layer.

        The activation energy, the difference between the actual adsorption energy and the saturated adsorption energy, is the determining factor of the ion adsorption saturation. Generally,higher activation energy indicates lower adsorption saturation.When the activation energy is zero, the ion can get closer to the particle surface (κl = 0). After determining the position of the ion in the diffuse layer, the actual adsorption energy of the ion and the surface can be calculated.

        Fig.2 Relationship between ion adsorption saturation and diffusion distance.κl is dimensionless distance.

        Fig.1 Maximum adsorption energy as a function of ionic activity in bulk solution for different surface charge densities.

        The metastable equilibrium of ion adsorption depends on the diffusion distance of ions in the diffuse layer. The actual adsorption energy is not the saturation adsorption energy, but the adsorption energy at the diffusion position l. If the diffusion distance of counterions from bulk phase (Gouy plane) equal 1/κ,the actual adsorption energy is equal to saturation adsorption,and the activation energy is approximately zero. Fig.3 shows that the activation energy decreases with increasing ion concentration. The activation energies for different cations exhibit Hofmeister effects. According to Eq.(8), the higher activation energy represents the lower adsorption energy.

        The adsorption energy of ions decreases with the increase of the ionic concentration in the diffuse layer, the series for ions is Cu2+> Ca2+> Mg2+. Although activation energies for the above bivalent cations were different (Fig.3), the adsorption saturation could be expressed as the same function of activation energy(Fig.4). In other words, the model of the interactions between ion and surface considering the Hofmeister effects is universal.

        Fig.3 Activation energies of cation exchange adsorption in different ion concentrations.

        The existing theories reveal that the ionic dispersion force and the volume effect are the main factors of the Hofmeister effects11,33. However, these factors only work at high concentrations, but the Hofmeister effects vanish slowly at low concentrations. In this study, the strong Hofmeister effects are observed from low concentration (0.05 mol·L?1) to high concentration (0.3 mol·L?1). The surface electric field of montmorillonite is 2.2 × 108V·m?1in this experiment, and the distribution of the electric field in the diffuse layer increases with ion concentration decreasing. Such a strong electric field can inevitably change the structure of electron clouds of the adsorbed ions. This effect is called the coupling effect between electric field and quantum fluctuation5. Recently, Parsons et al.34hold that surface charge transfer was a new origin of Hofmeister effects at low salt concentration. Therefore, we predict that the combined effects among ion dispersion forces,ionic volume, charge transfer and the coupling effect between electric field and quantum fluctuation exist in the interaction between ion and particle, but the main factor changes with different conditions. The coupling effect between the electric field and the quantum fluctuation for higher surface charge density and charge transfer for very low surface charge density are the main factor under low concentration, while the dispersion force and volume effect are enhanced at very high concentration.

        Fig.4 Relationship between the adsorption saturation and activation energy.

        5 Summary

        Ion-surface interactions strongly affect the ion adsorption. A new model of ion adsorption with Hofmeister effects was established, and the activation energy of ion exchange adsorption was estimated in this study. Ionic adsorption saturation at the montmorillonite surface obeys the sequence Cu2+> Ca2+> Mg2+. Hofmeister effects result in the difference of activation energy of ion exchange adsorption, and then determine the ion adsorption saturation. The calculated activation energies for the tested bivalent cations based on the new model was a same function of adsorption saturation. For the constant charge surfaces, the model of the interaction between ions and the surface considering the Hofmeister effects is universal, which must be included in the ion interface reactions.

        (1) Li, H.; Wu, L. S. Soil Sci. Soc. Am. J. 2007, 71, 1694.doi: 10.2136/sssaj2007.0004

        (3) Chen, X.; Yang, T.; Kataoka, S.; Cremer, P. S. J. Am. Chem.Soc. 2007, 129, 12272. doi: 10.1021/ja073869r

        (4) Borah, J. M.; Mahiuddin, S.; Sarma, N.; Parsons, D. F.;Ninham, B. W. Langmuir 2011, 27, 8710. doi:10.1021/la2006277

        (5) Liu, X.; Li, H.; Li, R.; Xie, D.; Ni, J.; Wu, L. Sci. Rep. 2014, 4,5047. doi: 10.1038/srep05047

        (6) Vereda, F.; Martin-Molina, A.; Hidalgo-Alvarez, R.;Quesada-Perez, M. Phys. Chem. Chem. Phys. 2015, 17, 17069.doi: 10.1039/C5CP01011J

        (7) Borukhov, I.; Andelman, D.; Orland, H. Phys. Rev. Lett. 1997,79, 435. doi: 10.1103/PhysRevLett.79.435

        (8) Ben-Yaakov, D.; Andelman, D.; Podgornik, R.; Harries, D.Curr. Opin. Colloid Interface Sci. 2011, 16, 542.doi: 10.1016/j.cocis.2011.04.012

        (9) Brown, M. A.; Bossa, G. V.; May, S. Langmuir 2015, 31,11477. doi: 10.1021/acs.langmuir.5b02389

        (10) Brown, M. A.; Abbas, Z.; Kleibert, A.; Green, R. G.; Goel, A.;May, S.; Squires, T. M. Physical Review X 2016, 6, 011007.doi: 10.1103/PhysRevX.6.011007

        (11) Parsons, D. F.; Bostr? m, M.; Nostro, P. L.; Ninham, B. W.Phys. Chem. Chem. Phys. 2011, 13, 12352.doi: 10.1039/C1CP20538B

        (12) Kunz, W.; Henle, J.; Ninham, B. Curr. Opin. Colloid Interface Sci. 2004, 9, 19. doi: 10.1016/j.cocis.2004.05.005

        (13) Jungwirth, P.; Cremer, P. S. Nature Chemistry 2014, 6, 261.doi: 10.1038/nchem.1899

        (14) Tian, R.; Yang, G.; Li, H.; Gao, X.; Liu, X.; Zhu, H.; Tang, Y.Phys. Chem. Chem. Phys. 2014, 16, 8828.doi: 10.1039/c3cp54813a

        (15) Xu, C.; Li, H.; Hu, F.; Li, S.; Liu, X.; Li, Y. Eur. J. Soil Sci.2015, 66, 615. doi: 10.1111/ejss.12252

        (16) Bolt, G. H. Soil Sci. 1955, 79, 267.doi: 10.1097/00010694-195504000-00004

        (17) Barber, J. Biochim. Biophys. Acta Bioenerg. 1980, 594, 253.doi: 10.1016/0304-4173(80)90003-8

        (18) Smit, W.; Holten, C. L. M. J. Colloid Interface Sci. 1980, 78,1. doi: 10.1016/0021-9797(80)90489-0

        (19) Bowen, W. R.; Mukhtar, H. Colloids Surf. A 1993, 81, 93.doi: 10.1016/0927-7757(93)80237-9

        (20) Mao, M.; Fornasiero, D.; Ralston, J.; Smart, R. S. C.; Sobieraj,S. Colloids Surf. A 1994, 85, 37.doi: 10.1016/0927-7757(94)02745-5

        (21) Mullet, M.; Fievet, P.; Reggiani, J. C.; Pagetti, J. J. Membr.Sci. 1997, 123, 255. doi: 10.1016/S0376-7388(96)00220-7

        (22) Boroudjerdi, H.; Kim, Y. W.; Naji, A.; Netz, R. R.;Schlagberger, X.; Serr, A. Phys. Rep. 2005, 416, 129.doi: 10.1016/j.physrep.2005.06.006

        (23) Dobrzynska, I.; Skrzydlewska, E.; Figaszewski, Z.Bioelectrochemistry 2006, 69, 142.doi: 10.1016/j.bioelechem.2006.01.004

        (24) Kinraide, T. B.; Wang, P. J. Exp. Bot. 2010, 61, 2507.doi: 10.1093/jxb/erq082

        (25) Li, H.; Hou, J.; Liu, X.; Wu, L. Soil Sci. Soc. Am. J. 2011, 75,2128. doi: 10.2136/sssaj2010.0301

        (26) Liu, X.; Li, H.; Li, R.; Tian, R.; Xu, C. Analyst 2013, 138,1122. doi: 10.1039/c2an36069a

        (27) Rowe, A. J. Biophys. Chem. 2001, 93, 93.doi: 10.1016/S0301-4622(01)00214-9

        (28) Bostr? m, M.; Williams, D. R. M.; Ninham, B. W. Phys. Rev.Lett. 2001, 87, 168103. doi: 10.1103/PhysRevLett.87.168103

        (29) Liu, X.; Li, H.; Li, R.; Tian, R.; Xu, C. Commun. Theor. Phys.2012, 58, 437. doi: 10.1088/0253-6102/58/3/20

        (30) Liu, X.; Li, H.; Du, W.; Tian, R.; Li, R.; Jiang, X. J. Phys.Chem. C 2013, 117, 6245. doi: 10.1021/jp312682u

        (31) Li, H.; Qing, C. L.; Wei, S. Q.; Jiang, X. J. J. Colloid Interface Sci. 2004, 275, 172. doi: 10.1016/j.jcis.2003.12.055

        (32) Low, P. F. Soil Sci. Soc. Am. J. 1980, 44, 667.doi: 10.2136/sssaj1980.03615995004400040001x

        (33) Lima, E. R. A.; Tavares, F. W.; Biscaia Jr, E. C. Phys. Chem.Chem. Phys. 2007, 9, 3174. doi: 10.1039/B701170A

        (34) Parsons, D. F.; Salis, A. Curr. Opin. Colloid Interface Sci.2016, 23, 41. doi: 10.1016/j.cocis.2016.05.005

        An Approach to Estimate the Activation Energy of Cation Exchange Adsorption

        LI Yan-Ting1LIU Xin-Min1,* TIAN Rui1DING Wu-Quan1XIU Wei-Ning2
        TANG Ling-Ling1ZHANG Jing1LI Hang1
        (1Chongqing Key Laboratory of Soil Multi-Scale Interfacial Process, College of Resources and Environment, Southwest University,Chongqing 400715, P. R. China;2Institute of Agricultural Engineering, Chongqing Academy of Agricultural Sciences,Chongqing 401329, P. R. China)

        Ion exchange adsorption is an important physicochemical process at solid/liquid interfaces. In this study, an approach was established to estimate the activation energy of cation exchange reaction on the charged surface, considering Hofmeister effects. The experimental results showed that Hofmeister effects strongly affect the ionic adsorption equilibrium on the charged particle surface. The position of the adsorbed counterion in the diffuse layer was predicted according to the established model, and the ion exchange activation energies for different bivalent cations were estimated via the cation exchange experiments. The activation energy decreases with increasing ion concentration, and the adsorption saturation of cations is a function of the activation energy at different concentrations. The established model of cation exchange adsorption in the present study has universal applicability in solid/liquid

        interface reactions.

        Specific ion effect; Ion adsorption; Charged surface; Activation energy

        December 29, 2016; Revised: May 8, 2017; Published online: May 18, 2017.*Corresponding author. Email: lucimir@163.com; Tel: +86-23-68250674.

        O641;O648;S153

        10.1351/pac200678081611

        [Article]

        10.3866/PKU.WHXB201705181 www.whxb.pku.edu.cn

        The project was supported by the National Natural Science Foundation of China (41530855, 41501240,41501241), Natural Science Foundation Project of CQ CSTC, China (cstc2015jcyjA80015), Young Talent Project of China Association for Science and Technology (2015?2017), Undergraduate Innovation Fund of Southwest University, China (20152402009) and “Guangjiong” Projects of Southwest University, China (2016005).

        國(guó)家自然科學(xué)基金(41530855, 41501240, 41501241), 重慶市自然科學(xué)基金(cstc2015jcyA80015), 中國(guó)科協(xié)青年人才托舉工程(2015?2017), 西南大學(xué)本科生創(chuàng)新基金(20152402009)和西南大學(xué)“光炯”創(chuàng)新基金(2016005)資助項(xiàng)目

        ? Editorial office of Acta Physico-Chimica Sinica

        猜你喜歡
        重慶市特異性重慶
        重慶客APP
        重慶市光海養(yǎng)蜂場(chǎng)
        蜜蜂雜志(2022年5期)2022-07-20 09:54:06
        重慶市光海養(yǎng)蜂場(chǎng)
        蜜蜂雜志(2022年3期)2022-05-23 05:11:42
        重慶人為什么愛吃花
        鑄造輝煌
        ——重慶市大足區(qū)老年大學(xué)校歌
        “逗樂(lè)坊”:徜徉相聲里的重慶味
        在這里看重慶
        今日重慶(2017年5期)2017-07-05 12:52:25
        精確制導(dǎo) 特異性溶栓
        BOPIM-dma作為BSA Site Ⅰ特異性探針的研究及其應(yīng)用
        重復(fù)周圍磁刺激治療慢性非特異性下腰痛的臨床效果
        久久久精品免费国产四虎| av免费一区二区久久| 在线观看国产成人av天堂野外| 日日碰狠狠添天天爽五月婷| 国语对白做受xxxxx在线中国 | 国产精品女丝袜白丝袜| 午夜免费观看一区二区三区| 亚洲av成人av三上悠亚| 亚洲性无码一区二区三区| 麻豆国产原创视频在线播放| 牛鞭伸入女人下身的真视频| 丝袜美腿网站一区二区| 亚洲精品久久久中文字| 亚洲中文字幕日韩综合| 亚洲图片日本视频免费| 亚洲精品你懂的在线观看| 免费人人av看| 日韩一区二区三区精品视频| 极品少妇xxxx精品少妇偷拍| 日本午夜精品一区二区三区电影| 富婆如狼似虎找黑人老外| 中文岛国精品亚洲一区| 免费毛片一区二区三区女同| 日本一区二区三区综合视频| 日日噜噜夜夜狠狠va视频v| 超碰97人人做人人爱少妇| 中文亚洲AV片在线观看无码| 国产熟妇一区二区三区网站| 国产精品国产三级国产a | 三级全黄裸体| 欧美丰满大屁股ass| 深夜福利国产| 国产白色视频在线观看| 久久精品国产色蜜蜜麻豆| 蜜桃精品免费久久久久影院| 亚洲性爱区免费视频一区| 亚洲一区二区三区2021| 97人妻人人做人碰人人爽| 日韩一线无码av毛片免费| 青青草国内视频在线观看| 最新在线观看免费的a站国产|