• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    Derivative-extremum analysis of current-potential curves showing electrochemical kinetics in the full reversibility range

    2023-03-14 06:52:50FengjunYinHongLiu
    Chinese Chemical Letters 2023年1期

    Fengjun Yin,Hong Liu

    a Chongqing Institute of Green and Intelligent Technology,Chinese Academy of Sciences,Chongqing 400714,China

    b Key Laboratory of Reservoir Aquatic Environment,Chinese Academy of Sciences,Chongqing 400714,China

    c China University of Chinese Academy of Sciences,Beijing 100049,China

    Keywords:Reversibility Electrochemical kinetics Half-wave potential Derivative-extremum analysis Parameter determinations

    ABSTRACT Derivative-extremum analysis (DEA) of j-E curves is a newly proposed method of half wave potential(E1/2) and activation feature extraction from steady-state voltammetry.Here,the DEA is demonstrated to be valid in the full range of reversibility using numerical simulations with a derived universal electrode equation,providing a novel perspective of electrochemical kinetics in the reversibility domain.The results reveal that E1/2 is a better choice of the reference potential instead of equilibrium potential (Eeq) in electrode equations,especially since Eeq is meaningless in an irreversible case.The equations referenced with standard potential, E1/2 and Eeq,are summarized in three tables,and their applications in parameter determinations are specified.Finally,reversibility is proved to be a relative measure between kinetic slowness and mass transport of electroactive species,and the reversibility classifications are proposed according to the DEA feature in the reversibility domain.This work,based on the DEA principle,refines the electrode equation forms and generalizes their applicability in the full range of reversibility.

    Steady-state current (j)-potential (E) voltammograms provide an easy-to-follow technique for extracting the kinetic information of an electrode reaction because of the simplicity of steady-state kinetic theory.In mathematics,steady-state voltammetry avoids charging current interference [1] and possesses time-independent mass transport processes between the electrode surface and the bulk solution [2].

    Reversibility is a basic scale in the recognition of the reactivity of a reaction,and the electrode equations are often divided into reversible,irreversible andquasi-reversible types [3].The electrode equations and their applicability are dependent on the reference potential against which the potential is quantified,e.g.,standard(or formal) potential,equilibrium potential (Eeq) and half wave potential (E1/2).The electrode equations are multifarious in terms of reversibility and reference potential.If possible,it is strongly desired to unify all the electrode equations and elaborately differentiate their applicability in steady-state voltammetry.

    “Reversible” and “irreversible” are two extreme cases; the former is only dependent on thermodynamics,while the latter is fully governed by kinetics.The variations in the kinetic features embodied in the steady-state voltammograms between the two cases have long been absent from a quantitative measure.E1/2is an important indicator of reaction occurrence,whereas its physical significance in the full range of reversibility has yet to be clarified.

    Recently,we found that the peak point of the derivative of aj-Ecurve in irreversible cases corresponded to a special kinetic state,where the peak potential was demonstrated to beE1/2and the peak value designated the activation feature of a reaction driven by potential [4].Additionally,this derivative extremum analysis (DEA)can be easily demonstrated to be applicable in reversible cases; although,the peak potential and peak value have different expressions.In this work,DEA was further proven to be valid in the full range of reversibility through numerical simulations using a derived universal electrode equation,and it provided a new perspective in understanding the kinetic feature evolution in the reversibility domain.We reorganized the equations referenced with formal potential,EeqandE1/2in the full reversibility scale and elucidated their applicability.Reversibility classifications were proposed according to the variation trends of DEA in the reversibility domain.

    To investigate the DEA applicability,we first derived a universal electrode equation in the full range of reversibility with three premises as follows: the uniformly accessible electrode treating the mass transport process,a large excess of supporting electrolyte to eliminate the migration effect,and a Nernst diffusion layer with stirring maintaining the bulk concentration constant.The cathodic current is defined as positive,and the anodic current is negative.

    The most basic equation describing an elementary reaction O+ne-=R is an original Butler-Volmer form (Eq.1) [5]:

    wherek0is the standard rate constant,E0is the standard potential and can be replaced by formal potential in practice,c*Oandc*Rare the oxidant and reductant concentrations on the electrode surface,respectively,αandβare the electron transfer coefficients (ETCs),α+β=1,Fis the Faraday constant,andf=F/RT,where R is the ideal gas constant,andTis the temperature in Kelvin.All the symbols are defined in Table S1 (Supporting information).

    Eq.1 is difficult to use sincec*Oandc*Rcan hardly be defined and measured.However,c*Oandc*Rcan be replaced by introducing the mass transport relationships of reactants (O and R).In the steady state,the reaction rate equals the mass transport fluxes of O and R between the electrode surface and bulk solution (Eqs.2)and 3:

    wherecbOandcbRare the concentrations in bulk solution,mis the mass transport coefficient,Dis the diffusion coefficient,andδis the diffusion layer thickness.The maximum transport fluxes of O and R are expressed as the cathodic and anodic limiting current densities (jl,candjl,a),which are obtained atc*Oandc*Rreaching zero (Eq.4):

    Introducing Eq.4 into Eqs.2 and 3 obtains (Eq.5):

    Finally,introducing Eq.5 into Eq.1 derives the following equation with the potential referenced againstE0(Eq.6):

    This equation is a universal form in consideration of the concentration polarization in the full range of reversibility,which excludes any simplification of kinetic factors.The electrode geometric size information can be introduced to the equation by using the analytical expressions of the diffusion layer thickness [6].All other equation forms,such as reversible and irreversible,in the presence of only O or R,and with equal mass transport coefficients,can be simplified from this equation.For example,in the presence of only O,jl,a=0,Eq.6 can be simplified to the following well-known equation proposed by Mirkin and Bard (Eq.7) [2]:

    Under the premise of equal mass transport coefficients(mO=mR),Eq.6 can be transformed to the following form as reported by Molinaet al.(Eq.8) [6]:

    whereη=nf(E-E0),=k0/mOandμ=cbR/cbO.

    If specifying the electrode geometric size inmOandmR,more detailed electrode equations hold for hemispherical and disc microelectrodes [3,7].Here,the electrode equations referenced withE0in the full range of reversibility are summarized in Table S2(Supporting information).

    The DEA derivation in reversible case is shown as follows.To obtain the reversible equations,Eq.6 is transformed to the following form (Eq.9):

    Under reversible conditions,mO/k0→0 andmR/k0→0.We proved thatmo/k0exp(αnf(E-E0))in the first term andmR/k0exp(-βnf(E-E0))in the second term are negligible; the detailed derivation is shown in Supporting information.Thus,a reversible equation form can be simplified from the above equation:

    By solving for the extremal solution of the derivative of Eq.10,we derived the expressions of the peak potential (EP) and peak vale(PV) in the dj/dEcurve and demonstrated thatEPwas identical toE1/2.The detailed derivation process is shown in the Supporting information.

    The DEA derivation in irreversible case is shown here.To date,the reported irreversible equations have basically applied the forms usingEeqas the reference potential.Eq.6 can be transformed to aEeq-referenced form.The exchange current density(j0) is an equilibrium current at the initial state wherec*O=cbOandc*R=cbRand has the following expressions (Eq.13):

    In combination with the expressions ofj0,jl,candjl,a,the cathodic and anodic terms in Eq.6 can be rearranged as follows (Eqs.14 and 15):

    Please note that the above transformations are only valid atcbOandcbRare not zero.Based on the above relationships,Eq.6 can be transformed to the following form referenced withEeq(Eq.16):

    In the irreversible case,j0is so small that a very high overpotential is required to drive the reaction occurrence.Thus,the term exp(βnf(E-Eeq))for the cathode reaction is negligible,and exp(-αnf(E-Eeq))for the anode reaction is negligible.Eq.16 can be simplified to the following cathodic and anodic irreversible equations (Eqs.17 and 18):

    Cathode:

    Fig.1.Calculation results of Eq.6 with lgk0 varying from 0 to–9, cbO=cbR=10 mol/m3,and mO=mR=5×10–5 m/s, n=1, E0=0 V, α=0.5 and T = 25 °C: (a) j-E curves,the dashed lines were calculated from the reversible equation (Eq.10) with lgk0=0 and irreversible equations (Eqs.17 and 18) with lgk0=–9,(b) associated dj/dE curves,(c)peak potential (EP) and potentials at jl,c/2 and jl,a/2 plotted vs. lgk0,and (d) PV plotted vs. lgk0,the lines were calculated from the reversible and irreversible expressions of PV.

    Anode:

    These equations have similar forms to those derived from the multielectron and multistep irreversible process with one sole rate-limiting step [4].The DEA applicability in irreversible cases has been demonstrated by solving for the extremal solution of the derivatives of Eqs.17 and 18 in a previous work [4].The expressions ofE1/2and PV in irreversible cathodic reactions are shown as follows (Eqs.19 and 20):

    In both reversible and irreversible cases,Epis identical toE1/2,and PV has a unit of A m-2V-1and can serve to characterize the kinetic slowness of a reaction driven by potential.However,E1/2and PV have different expressions in the two cases.

    Furthermore,the DEA simulations in the full range of reversibility using Eq.6 are calculated to show the DEA in the transition range between reversible and irreversible cases,because mathematically solving for the derivatives of Eq.6 is too difficult.In a given scenario,a set ofj-Ecurves withk0values varying from 0 m/s to 10-10m/s are calculated first,then the dj/dEcurves are calculated,and theEpand PV values are extracted from the peak points of dj/dEcurves.Finally,the plots ofEpand PVvs.lgk0are obtained.For all the calculations,n=1,α=β=0.5,E0=0 Vvs.SHE andT= 25 °C.The other parameter values (mO,mR,cbOandcbR) are presented in the figure titles.The calculations were performed using MATLAB 2012a.

    The calculation results in the presence of both O and R are shown in Fig.1.Figs.1a clearly presents the variation trend of thej-Ecurves from the reversible to irreversible case: in the condition ofk0>10-4m/s,the curves approached the calculation result of the reversible equation (Eq.10) withk0=1 m/s; oncek0was small enough,namely,k0=10-9m/s,the curves were consistent with the results of irreversible equations (Eqs.17 and 18).In Fig.1b,the associated dj/dEcurves show that the peak point gradually shifts from one point in the reversible case to two peak points.The isolation degree between the cathode and anode reactions can be characterized by the distance between the two points and is intensified ask0decreases.

    To further investigate the relationship betweenEpandE1/2,theEpand the potentials ofjl,c/2 andjl,a/2 were plotted in Fig.1c.Epchanged suddenly from the reversibleE1/2feature,namely,the potential of (jl,c+jl,a)/2,to the irreversibleE1/2feature,namely,the potentials ofjl,c/2 andjl,a/2.In Fig.1d,the PV gradually shifts from the reversible case,as designated by the line calculated from Eq.12 to the irreversible case designated by the line calculated from Eq.20.These results confirmed that DEA is applicable in most cases in the presence of both O and R,except for a narrow transition range between reversible and irreversible cases where the exactE1/2can hardly be defined.In light of the variation features ofEpand PV,the exact reversible and irreversible ranges can be defined as shown by the colored regions.

    Fig.2.Calculation results of Eq.6 with lgk0 varying from 0 to–10, cbO = 10 mol/m3, cbR = 0, mO = mR=5×10–5 m/s, n=1, E0=0 V, α=0.5 and T = 25 °C: (a) j-E curves,(b) associated dj/dE curves,(c) Ep vs. lgk0 and (d) PV vs. lgk0,the lines were calculated from the reversible and irreversible expressions of Ep and PV.

    Table 1 Electrode equation forms referenced with E1/2 in the full range of reversibility.

    Table 2 Electrode equation forms referenced with Eeq in the full range of reversibility.

    The DEA simulation results in the presence of only O are shown in Fig.2.Figs.2a and b demonstrate that theEpin the dj/dEcurves corresponded to the potential ofjl,c/2 in thej-Ecurves in the full range of reversibility.Thus,Epis invariably identical toE1/2in this case.TheEpvs.lgk0plot in Fig.2c also presented a sudden change from the reversibleE1/2feature to the irreversibleE1/2feature.In Fig.2d,PV gradually shifts from the reversible case to the irreversible case.

    Finally,the above DEA simulation results confirmed thatE1/2is a kinetically special potential corresponding to the positions of the half limiting current in thej-Ecurves and the peak point in the dj/dEcurve.DEA is a general method to extractE1/2as well as PV in the full range of reversibility.The variation trends ofE1/2and PV along thek0axis can give a classification criterion ofquasi-reversible andquasi-irreversible regions,namely,(i)quasi-reversible:E1/2is approximate to the reversibleE1/2value;(ii)quasi-irreversible:E1/2is approximate to the irreversibleE1/2value.

    The above calculations indicate thatE1/2is a good choice of reference potential because it can be steadily determined in the full range of reversibility.Here,the electrode equations referenced withE1/2in the full range of reversibility are summarized in Table 1.The reversible equation form (Eq.21) can be obtained by combining Eqs.10 and 11:

    TheE1/2-referenced equations of irreversible cathode and anode reactions can be derived by introducing the irreversible expression ofE1/2(Eq.19) into Eqs.17 and 18,respectively (Eqs.22 and 23):

    Cathode:

    Anode:

    In irreversible cases,it is not necessary to distinguish the equations in the presence of both O and R and in the presence of only O because the anode and cathode reactions are evidently separated from each other.

    Inquasi-reversible/irreversible cases,E1/2-referenced equations are not available because the exact expression ofE1/2and PV can hardly be derived.Nonetheless,DEA is still applicable.

    Furthermore,the electrode equations referenced withEeqin the full range of reversibility are summarized in Table 2.The reversible equation is derived as follows.Simultaneously,dividing byj0in the numerator and denominator of Eq.16,it becomes (Eq.24):

    In the reversible case,j0is so large that 1/j0is negligible.The above equation can be transformed to the following reversible form (Eq.25):

    Although irreversibleEeq-referenced equations are derived and are as given in Table 2,users should be highly cautious about their applications.The applications ofEeqandj0are introduced in the following part.TheEeqused as a reference potential in kinetics is the open circuit potential of an electrode reaction,which satisfies the Nernst relationship (Eq.26):

    At this potential,the cathodic and anodic rates reach balance and are expressed as exchange current density (j0) (Eq.27):

    It is obvious thatEeqandj0are relevant only in the presence of both O and R.OncecbOorcbRapproaches zero,Eeqandj0in Eqs.26 and 27 become an infinite problem without explicit value.Moreover,Eeqis experimentally indeterminable in irreversible cases because the zero current is not a point but a wide gap between the cathode and anode reactions (Fig.1).Thus,Eeqandj0are applicable for the kinetically fast reactions in the presence of both O and R.

    Despite this,Eeqandj0are still frequently applied in irreversible cases.Here,we demonstrate thatEeqandj0are an interdependent pair in irreversible cases.Choosing any potential in the zero current stage as a formal equilibrium potential (E*eq),the deviation ofE*eqfrom the theoreticalEeqcan be described byEeq=E*eq+ΔE.Introducing it to the irreversible equation (Eq.17)obtains Eq.28:

    It revealed that the application ofE*eqonly results in the variation ofj0.In irreversible cases,Eeqis indeterminable and is usually determined as the onset potential,the lowest potential where the faradaic current is observed [8].Therefore,thej0value has little physical significance in irreversible cases.If selecting the sameEeqvalue,j0is still useful to compare the activity of a reaction between different conditions or catalysts.

    In addition,the activation overpotential (ηact) in the irreversible case is also meaningless because it is dependent on thej0andEeqvalues.Nonetheless,the overpotential analysis of the mass transport effect is feasible,viz.,the concentration overpotential assigned to reactant mass transport [9] and the pH overpotential assigned to H+/OH-transport [10],because the overpotential partE-Eeq-ηacthas excluded theEeqvalue.

    Currently,theEeq-referenced equations are typically used to fit thej-Ecurves for parameter determinations,such as the Tafel fitting method and a recently reported nonlinear fitting method [11].The Tafel method either directly ignores the mass transport effect by fitting a narrowj-Ecurve range [12,13] or corrects the mass transport current part to improve the accuracy of parameter values [14–16].The nonlinear fitting method is performed based on Eqs.17 or 18 to fit the entirej-Ecurve range,which overcomes the disadvantages of the Tafel method to obtain reliable ETC,j0andjlvalues.

    TheE1/2-referenced equations are more significant for quantifying the activity of irreversible reactions.The information ofEeqandj0is involved inE1/2(see Eq.19) to reveal the priority of reaction occurrence.The nonlinear fitting method can also be performed using the irreversibleE1/2-referenced equations and obtain the values of ETC,E1/2andjl.In Fig.S1 (Supporting information),the nonlinear fitting results of the two types of equations (Eqs.17 and 22)showed that the fitting lines were totally identical and could well fit the entirej-Ecurve.All the parameter values (Table S3 in Supporting information) determined from Eq.22 have a narrow 95%confidence interval,and the intervals of ETC,E1/2andjl,aare less than 1.10%,0.06% and 0.25%,respectively.This confirms that the nonlinear fits based on the irreversibleE1/2-referenced equations are a reliable parameter determination method of irreversible reactions.

    Moreover,the irreversibleE1/2-referenced equations are also applicable in conditions deviating from the irreversible case,wherein the ETC reveals reversibility variations.The reversibleE1/2-referenced equations are the special form of irreversibleE1/2-referenced equations with ETC equal to 1.With the aid of micro- and nanoscale electrode techniques to measure steady-state voltammograms [17],the practical application scope of irreversibleE1/2-referenced equations can be greatly extended.

    Importantly,the variation trends ofE1/2and PV in DEA in the reversibility domain provide classification criteria of reversible,irreversible,quasi-reversible andquasi-irreversible cases.Here,the influences of transport coefficients (mOandmR) and concentrations (cbOandcbR) on the reversibility classifications are investigated.In the presence of only O,the effect ofmOwith ak0scale is shown in Figs.3a–c.The reversible and irreversible ranges presented by the colored regions gradually shift toward the decreasing direction of lgk0asmOdecreases from 10-4m/s to 10-6m/s.The transition region can be divided intoquasi-reversible andquasi-irreversible parts bounded by lgmOas indicated by the dashed lines in Figs.3a–c.

    Fig.3. Ep and PV vs. lgk0 plots in the presence of only O calculated from Eq.6 with cbO=10 mol/m3, n=1, α=0.5, E0=0 V and T = 25 °C: (a) lgmO=–4,(b) lgmO=–5,and (c) lgmO=–6; (d) plots of Ep in Fig.(a–c) vs. lg(k0/mO).

    TheEpvs.lgk0plots in Figs.3a–c were redrawn using a logarithmick0/mOscale as shown in Fig.3d,and the results at differentmOconditions were the same.It revealed thatk0/mOwas a dimensionless parameter and the most basic scale for reversibility classification.Hence,reversibility was a relative measure between reaction kinetics and mass transport of electroactive species.Roughly,

    k0/mO?1: reversible

    k0/mO?1: irreversible

    In theory,the full “reversible” and “irreversible” are two extreme cases that can hardly be reached,and all the reactions may be regarded asquasi-states [3].Here,a detailed classification of the effectivelyquasi-reversible andquasi-irreversible cases was proposed,as listed in Table 3,namely,

    100>k0/mO>1:quasi-reversible,whereE1/2approaches the reversibleE1/2value

    0.01<k0/mO<1:quasi-irreversible,whereE1/2approaches the irreversibleE1/2value

    Thek0scale for the reversibility classification is dependent on themOand/ormRvalues.In practice,the diffusion coefficient ormass transport coefficient values of most solutes rarely differ by more than one order of magnitude.The diffusion coefficient usually has a magnitude of 10-9m2/s [18,19],and the mass transport coefficient has a magnitude of 10-5m/s [20].The reversibility classification in terms of thek0scale withmO=10-5m/s (see Fig.3b)is listed in Table 3.

    Table 3 Reversibility classifications in terms of k0/mO and k0 in the presence of only O.

    In addition,j0,in comparison to the transport limiting current(jl,corjl,a),is also frequently applied to the quantification of reversibility in electrochemistry [21].The ratio betweenj0andjl,ccan be derived from Eqs.4 and 27 as follows (Eq.29):

    This indicates thatj0/jl,cis an extended form ofk0/mOthat involves an extra concentration term.Obviously,j0/jl,cis applicable in the presence of both O and R and is identical tok0/mOin the condition ofcbO=cbR.The DEA calculation result under the condition ofcbO=cbRis shown in Fig.S2 (Supporting information),confirming that the reversibility scale is independent of the concentration extent.Under the condition ofcbO/=cbR(Fig.S3 in Supporting information),the reversible and irreversible regions can be well identified but thequasi-state regions for anode and cathode reactions become very complicated.

    In this work,numerical simulations of DEA using a derived universal electrode equation successfully verified the validity of DEA in the full range of reversibility and gained panoramic insight into electrochemical kinetics.The equations referenced withE0,E1/2andEeqin conjunction with their applicability were reorganized in three tables in terms of the reversibility scale.Reversibility classifications in terms of the variation features of DEA in the reversibility domain were proposed.

    The outcomes of this work are especially important for reactions suffering from mass transport limitations.One can apply the DEA method to extractE1/2and PV from the steady statej-Ecurves.E1/2is a potential feature parameter pronouncing the priority of reaction occurrence on a potential scale.PV is an activation feature parameter revealing the kinetic slowness driven by potential.Alternatively,the nonlinear fitting method ofj-Ecurves using the irreversibleE1/2-referenced equations can obtain reliable ETC,E1/2andjlvalues.In summary,theE1/2-referenced equations have wide application scope in the presence of only O,while theEeq-referenced equations are primarily useful for the kinetically fast reactions in the presence of both O and R.Users should be highly cautious about the applications ofEeqandj0in irreversible cases because they are an interdependent pair and have little physical significance.

    Declaration of competing interest

    The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper.

    Acknowledgment

    This work was financially supported by the National Natural Science Foundation of China (Nos.52131003,52170059,51808526,51727812).

    Supplementary materials

    Supplementary material associated with this article can be found,in the online version,at doi:10.1016/j.cclet.2022.01.078.

    夜夜夜夜夜久久久久| 嫩草影视91久久| 久久午夜亚洲精品久久| 久久国产精品影院| 国产黄色小视频在线观看| 国产色婷婷99| 哪里可以看免费的av片| 日韩欧美国产在线观看| 成人三级黄色视频| 99久久精品国产亚洲精品| 波多野结衣高清无吗| 午夜精品在线福利| 悠悠久久av| 亚洲成人久久爱视频| 色哟哟哟哟哟哟| 国产精品嫩草影院av在线观看 | 国产av不卡久久| av专区在线播放| 国产一区二区亚洲精品在线观看| 欧美黑人欧美精品刺激| 两个人的视频大全免费| 国产熟女xx| 色av中文字幕| 真人一进一出gif抽搐免费| 午夜免费成人在线视频| 极品教师在线视频| 久久久成人免费电影| 精品一区二区三区av网在线观看| 有码 亚洲区| 国产伦在线观看视频一区| 狠狠狠狠99中文字幕| 757午夜福利合集在线观看| 18美女黄网站色大片免费观看| 91麻豆精品激情在线观看国产| 男人舔女人下体高潮全视频| 亚洲一区二区三区色噜噜| 在线观看舔阴道视频| 久久精品国产亚洲av天美| 亚洲人成电影免费在线| 极品教师在线免费播放| 欧美日韩综合久久久久久 | av在线蜜桃| x7x7x7水蜜桃| 久久精品久久久久久噜噜老黄 | 国产精品1区2区在线观看.| 国产精华一区二区三区| 亚洲专区中文字幕在线| 九色国产91popny在线| 啦啦啦韩国在线观看视频| 久久久精品大字幕| 18+在线观看网站| 两个人视频免费观看高清| 中文亚洲av片在线观看爽| 久99久视频精品免费| 国产野战对白在线观看| 免费黄网站久久成人精品 | 极品教师在线视频| 亚洲av免费在线观看| 亚洲av电影在线进入| 成人欧美大片| 日韩欧美免费精品| 网址你懂的国产日韩在线| 哪里可以看免费的av片| 99精品在免费线老司机午夜| 国内毛片毛片毛片毛片毛片| 国产精品1区2区在线观看.| 精品熟女少妇八av免费久了| 少妇人妻一区二区三区视频| 国产av一区在线观看免费| 亚洲内射少妇av| 成人高潮视频无遮挡免费网站| a在线观看视频网站| 欧美+亚洲+日韩+国产| 精品99又大又爽又粗少妇毛片 | 亚洲真实伦在线观看| 色尼玛亚洲综合影院| 精品久久久久久久久av| 人妻制服诱惑在线中文字幕| 国产亚洲av嫩草精品影院| 国产爱豆传媒在线观看| 国产单亲对白刺激| 精品午夜福利在线看| 欧美+日韩+精品| 性色avwww在线观看| 国产精品影院久久| 国产aⅴ精品一区二区三区波| 国产一区二区亚洲精品在线观看| 老熟妇乱子伦视频在线观看| 国产真实乱freesex| 久久精品久久久久久噜噜老黄 | 国产免费一级a男人的天堂| 又爽又黄无遮挡网站| 亚洲自偷自拍三级| 伦理电影大哥的女人| 久久久国产成人免费| 久久精品国产自在天天线| 免费人成在线观看视频色| 欧美成人a在线观看| 国产三级中文精品| 乱码一卡2卡4卡精品| 国产伦人伦偷精品视频| 三级国产精品欧美在线观看| 老司机午夜福利在线观看视频| 五月玫瑰六月丁香| 国产精品一区二区三区四区免费观看 | 国产成人福利小说| h日本视频在线播放| 国产日本99.免费观看| 五月玫瑰六月丁香| 天堂√8在线中文| 国产蜜桃级精品一区二区三区| 搡老岳熟女国产| 亚洲精品一区av在线观看| 午夜福利高清视频| 一个人看的www免费观看视频| 蜜桃亚洲精品一区二区三区| 偷拍熟女少妇极品色| 国产爱豆传媒在线观看| 亚洲在线观看片| 久久久久九九精品影院| 久久精品夜夜夜夜夜久久蜜豆| 在线播放无遮挡| 国产精品美女特级片免费视频播放器| 国产精品一及| 婷婷精品国产亚洲av| 成人国产一区最新在线观看| 日韩成人在线观看一区二区三区| 午夜福利欧美成人| 日本 av在线| 精品人妻视频免费看| 国产精品综合久久久久久久免费| 在线天堂最新版资源| 又黄又爽又免费观看的视频| 久久久久久久久久成人| 99久久无色码亚洲精品果冻| 国产精品一区二区三区四区免费观看 | 久久6这里有精品| 亚洲国产精品sss在线观看| 欧美日韩国产亚洲二区| aaaaa片日本免费| 久久精品综合一区二区三区| 国内精品美女久久久久久| 日本一本二区三区精品| 十八禁人妻一区二区| 91麻豆精品激情在线观看国产| 黄色日韩在线| 国产成人啪精品午夜网站| 免费人成视频x8x8入口观看| 丁香欧美五月| 搞女人的毛片| 欧美3d第一页| 国内毛片毛片毛片毛片毛片| av在线观看视频网站免费| 国产av麻豆久久久久久久| 国产真实乱freesex| 国产一区二区在线av高清观看| 亚洲一区二区三区色噜噜| 久久久精品大字幕| 99在线视频只有这里精品首页| 五月玫瑰六月丁香| 我的老师免费观看完整版| 日韩大尺度精品在线看网址| 午夜精品久久久久久毛片777| 一二三四社区在线视频社区8| 动漫黄色视频在线观看| 日本一本二区三区精品| 欧美成狂野欧美在线观看| 欧美黄色片欧美黄色片| 夜夜看夜夜爽夜夜摸| 十八禁国产超污无遮挡网站| 亚洲欧美清纯卡通| 天堂√8在线中文| ponron亚洲| 人人妻人人看人人澡| 久久久久性生活片| 俄罗斯特黄特色一大片| 少妇人妻一区二区三区视频| 欧美日韩瑟瑟在线播放| 欧美国产日韩亚洲一区| 午夜老司机福利剧场| 美女高潮喷水抽搐中文字幕| 日本黄大片高清| 成人一区二区视频在线观看| 三级毛片av免费| 噜噜噜噜噜久久久久久91| 自拍偷自拍亚洲精品老妇| 欧美乱色亚洲激情| 成人亚洲精品av一区二区| 亚洲五月天丁香| 亚洲av美国av| 欧美+亚洲+日韩+国产| 亚洲,欧美,日韩| 日韩欧美在线二视频| a在线观看视频网站| 国产欧美日韩精品亚洲av| 中文字幕av成人在线电影| 久久久久久久久久成人| 亚洲国产精品久久男人天堂| 一本久久中文字幕| 哪里可以看免费的av片| 亚洲真实伦在线观看| 少妇高潮的动态图| 十八禁国产超污无遮挡网站| 久99久视频精品免费| 日本免费a在线| 亚洲激情在线av| 国内少妇人妻偷人精品xxx网站| 国产亚洲精品av在线| 亚洲男人的天堂狠狠| 禁无遮挡网站| 成人特级av手机在线观看| 亚洲欧美日韩无卡精品| 日本成人三级电影网站| 亚洲三级黄色毛片| 两人在一起打扑克的视频| 一个人观看的视频www高清免费观看| 国产精品,欧美在线| 午夜免费成人在线视频| 国内精品一区二区在线观看| 久久久久精品国产欧美久久久| 国产欧美日韩一区二区精品| 99国产综合亚洲精品| 最后的刺客免费高清国语| 一级作爱视频免费观看| 久久99热6这里只有精品| 91麻豆av在线| 国产大屁股一区二区在线视频| 欧美潮喷喷水| 亚洲美女搞黄在线观看 | 欧美乱色亚洲激情| 夜夜躁狠狠躁天天躁| 精品日产1卡2卡| 久9热在线精品视频| 日日干狠狠操夜夜爽| 国产av不卡久久| 男女下面进入的视频免费午夜| 成人av一区二区三区在线看| 国产探花极品一区二区| 国产欧美日韩精品一区二区| 欧美精品啪啪一区二区三区| 亚洲 欧美 日韩 在线 免费| 久久久久亚洲av毛片大全| 婷婷色综合大香蕉| 88av欧美| 国产伦人伦偷精品视频| 精品熟女少妇八av免费久了| 97热精品久久久久久| 国产伦人伦偷精品视频| 亚洲精品粉嫩美女一区| 国产主播在线观看一区二区| 91字幕亚洲| 婷婷精品国产亚洲av| 91狼人影院| 精品国内亚洲2022精品成人| 亚洲男人的天堂狠狠| 国产视频内射| 欧美乱妇无乱码| 亚洲黑人精品在线| 亚洲av电影在线进入| 久久久久久久久久黄片| 男女下面进入的视频免费午夜| 嫩草影视91久久| 国产午夜福利久久久久久| 伦理电影大哥的女人| 国产亚洲精品av在线| 亚洲欧美日韩高清在线视频| 久久精品国产亚洲av天美| 欧美3d第一页| 国产成人av教育| 欧美成狂野欧美在线观看| 毛片一级片免费看久久久久 | 在线十欧美十亚洲十日本专区| 麻豆一二三区av精品| 国产aⅴ精品一区二区三区波| 91麻豆av在线| aaaaa片日本免费| 亚洲人成电影免费在线| 久久久久九九精品影院| 欧美在线一区亚洲| 黄色视频,在线免费观看| 精品久久国产蜜桃| 久久久久久久久久黄片| 草草在线视频免费看| 无遮挡黄片免费观看| 成人国产一区最新在线观看| 亚洲中文字幕一区二区三区有码在线看| 国产伦一二天堂av在线观看| 成人欧美大片| 99热精品在线国产| 欧美高清性xxxxhd video| 亚洲一区高清亚洲精品| x7x7x7水蜜桃| bbb黄色大片| 成人av在线播放网站| 男人和女人高潮做爰伦理| 桃红色精品国产亚洲av| 在线观看午夜福利视频| 久久久久久久精品吃奶| 制服丝袜大香蕉在线| 99国产综合亚洲精品| 亚洲狠狠婷婷综合久久图片| 一二三四社区在线视频社区8| 黄片小视频在线播放| 国产免费av片在线观看野外av| 国产一区二区在线观看日韩| 日韩成人在线观看一区二区三区| av中文乱码字幕在线| 亚洲精品在线观看二区| 91午夜精品亚洲一区二区三区 | 久久久久久久久久成人| 亚洲精品亚洲一区二区| 成年女人看的毛片在线观看| 亚洲中文字幕日韩| 夜夜看夜夜爽夜夜摸| 国产精品久久久久久久电影| 免费av不卡在线播放| 亚洲av五月六月丁香网| 伦理电影大哥的女人| 国产白丝娇喘喷水9色精品| 欧美精品啪啪一区二区三区| 伊人久久精品亚洲午夜| 少妇熟女aⅴ在线视频| 内射极品少妇av片p| 欧美+日韩+精品| 成人高潮视频无遮挡免费网站| 国产精品乱码一区二三区的特点| 欧美绝顶高潮抽搐喷水| 亚洲欧美日韩东京热| 老熟妇乱子伦视频在线观看| h日本视频在线播放| 国产精品1区2区在线观看.| 欧美最新免费一区二区三区 | 国产免费一级a男人的天堂| 悠悠久久av| 日韩精品青青久久久久久| 波多野结衣高清无吗| 最好的美女福利视频网| 亚洲欧美日韩无卡精品| 最新中文字幕久久久久| 亚洲综合色惰| 一二三四社区在线视频社区8| 精品久久久久久久久av| 国产一区二区三区在线臀色熟女| 18禁在线播放成人免费| 91在线精品国自产拍蜜月| 毛片一级片免费看久久久久 | 最好的美女福利视频网| 欧美性猛交╳xxx乱大交人| 国产精品美女特级片免费视频播放器| 国产高清三级在线| 又紧又爽又黄一区二区| 最新中文字幕久久久久| 久久6这里有精品| 欧美色欧美亚洲另类二区| 内地一区二区视频在线| 啦啦啦韩国在线观看视频| 国产伦一二天堂av在线观看| 亚洲av二区三区四区| 欧美日韩黄片免| 亚洲真实伦在线观看| 最新中文字幕久久久久| 99热只有精品国产| a级毛片免费高清观看在线播放| 搡老妇女老女人老熟妇| 午夜老司机福利剧场| 国产成人aa在线观看| 91麻豆av在线| 99热这里只有是精品50| 精品国产亚洲在线| 欧美日韩中文字幕国产精品一区二区三区| 午夜精品久久久久久毛片777| 99热精品在线国产| 免费在线观看影片大全网站| 十八禁网站免费在线| 国产亚洲av嫩草精品影院| 亚洲av成人精品一区久久| 日韩欧美 国产精品| 午夜a级毛片| 亚洲精品粉嫩美女一区| 亚洲欧美日韩卡通动漫| 99在线视频只有这里精品首页| 成人鲁丝片一二三区免费| 中文字幕人妻熟人妻熟丝袜美| 少妇的逼水好多| 国产精品1区2区在线观看.| 精品人妻偷拍中文字幕| 少妇被粗大猛烈的视频| 免费高清视频大片| 少妇裸体淫交视频免费看高清| 桃红色精品国产亚洲av| 日本撒尿小便嘘嘘汇集6| 午夜久久久久精精品| 中文资源天堂在线| 小说图片视频综合网站| 久久精品人妻少妇| 一进一出好大好爽视频| av在线天堂中文字幕| 久久久久九九精品影院| 国产69精品久久久久777片| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 免费在线观看日本一区| 搞女人的毛片| 成人av一区二区三区在线看| 午夜福利欧美成人| 欧美一区二区国产精品久久精品| 老鸭窝网址在线观看| 无人区码免费观看不卡| av黄色大香蕉| 日日摸夜夜添夜夜添av毛片 | 99久久99久久久精品蜜桃| 国产亚洲精品久久久久久毛片| 亚洲自拍偷在线| 长腿黑丝高跟| 色精品久久人妻99蜜桃| 无人区码免费观看不卡| 国产aⅴ精品一区二区三区波| 啦啦啦韩国在线观看视频| 亚洲国产精品合色在线| 动漫黄色视频在线观看| eeuss影院久久| a级毛片a级免费在线| 美女cb高潮喷水在线观看| 99久国产av精品| 小蜜桃在线观看免费完整版高清| 欧美区成人在线视频| 成人三级黄色视频| 欧美日韩中文字幕国产精品一区二区三区| 欧美xxxx黑人xx丫x性爽| 久久久国产成人精品二区| 日日摸夜夜添夜夜添小说| 99视频精品全部免费 在线| 97人妻精品一区二区三区麻豆| 亚洲熟妇熟女久久| 亚洲五月天丁香| 99久国产av精品| 99久久99久久久精品蜜桃| 亚洲一区二区三区色噜噜| www.色视频.com| 亚洲av熟女| 又爽又黄无遮挡网站| 久久久久久久久中文| 1000部很黄的大片| 美女黄网站色视频| 夜夜爽天天搞| 蜜桃亚洲精品一区二区三区| 18+在线观看网站| 老熟妇仑乱视频hdxx| 亚州av有码| 一本久久中文字幕| 国产久久久一区二区三区| 欧美绝顶高潮抽搐喷水| 国产av在哪里看| 久久久久久大精品| 能在线免费观看的黄片| 国产精品一区二区免费欧美| .国产精品久久| 日本 av在线| 成人午夜高清在线视频| 97超级碰碰碰精品色视频在线观看| 日日摸夜夜添夜夜添小说| 三级国产精品欧美在线观看| 亚洲激情在线av| 亚洲精品久久国产高清桃花| 精品国产亚洲在线| 丰满的人妻完整版| 性插视频无遮挡在线免费观看| 精品日产1卡2卡| 黄色配什么色好看| 美女免费视频网站| 狠狠狠狠99中文字幕| 亚洲精品456在线播放app | 人妻久久中文字幕网| 99国产精品一区二区三区| 亚洲成人精品中文字幕电影| 国产精品自产拍在线观看55亚洲| 非洲黑人性xxxx精品又粗又长| 赤兔流量卡办理| 日本五十路高清| 99国产综合亚洲精品| 美女高潮的动态| 亚洲国产日韩欧美精品在线观看| 好男人电影高清在线观看| 国产黄a三级三级三级人| 精品人妻熟女av久视频| 可以在线观看毛片的网站| 国产高清视频在线播放一区| 少妇的逼好多水| 精品熟女少妇八av免费久了| 欧美最黄视频在线播放免费| av在线蜜桃| 丰满乱子伦码专区| 国产成人av教育| 久久久久国产精品人妻aⅴ院| 一个人免费在线观看的高清视频| 中文字幕人妻熟人妻熟丝袜美| 久久人人精品亚洲av| 亚洲欧美日韩高清专用| 成人高潮视频无遮挡免费网站| 少妇的逼好多水| 色尼玛亚洲综合影院| 2021天堂中文幕一二区在线观| 最近最新免费中文字幕在线| 久久人妻av系列| 大型黄色视频在线免费观看| 欧美一区二区亚洲| 国产成人啪精品午夜网站| 中文亚洲av片在线观看爽| 少妇的逼好多水| 最新在线观看一区二区三区| 黄色一级大片看看| 国产91精品成人一区二区三区| 最近在线观看免费完整版| 欧美高清成人免费视频www| 丁香欧美五月| 日韩欧美精品免费久久 | 欧美bdsm另类| 美女cb高潮喷水在线观看| 麻豆一二三区av精品| 我的女老师完整版在线观看| 欧美国产日韩亚洲一区| 在线观看一区二区三区| 亚洲美女搞黄在线观看 | 99在线视频只有这里精品首页| 波野结衣二区三区在线| 嫩草影院新地址| 18+在线观看网站| 国产欧美日韩精品亚洲av| 久久精品国产亚洲av天美| 亚洲乱码一区二区免费版| 亚洲国产色片| 在线观看66精品国产| av视频在线观看入口| 久久久久久久久久成人| 国产私拍福利视频在线观看| 中文字幕精品亚洲无线码一区| 亚洲欧美日韩东京热| 日本黄色片子视频| 麻豆成人av在线观看| 女生性感内裤真人,穿戴方法视频| 夜夜夜夜夜久久久久| 又紧又爽又黄一区二区| 九九热线精品视视频播放| 婷婷丁香在线五月| 身体一侧抽搐| 一夜夜www| 啦啦啦观看免费观看视频高清| 最好的美女福利视频网| 成人无遮挡网站| 在线国产一区二区在线| 男人舔女人下体高潮全视频| 亚洲国产色片| 成年版毛片免费区| 国产三级黄色录像| 欧美成人免费av一区二区三区| 欧美日韩黄片免| 一区二区三区四区激情视频 | 狠狠狠狠99中文字幕| 亚洲av成人精品一区久久| 美女大奶头视频| 国产淫片久久久久久久久 | 成年版毛片免费区| 国产成人福利小说| 日本三级黄在线观看| 亚洲中文字幕日韩| 偷拍熟女少妇极品色| 99热这里只有是精品50| 97人妻精品一区二区三区麻豆| 男人和女人高潮做爰伦理| 欧美绝顶高潮抽搐喷水| 欧美黑人欧美精品刺激| 久久精品国产99精品国产亚洲性色| 不卡一级毛片| 岛国在线免费视频观看| a级一级毛片免费在线观看| 国产欧美日韩一区二区三| 首页视频小说图片口味搜索| 我的老师免费观看完整版| 久久亚洲真实| 成人永久免费在线观看视频| 亚洲av免费在线观看| 色av中文字幕| 欧美又色又爽又黄视频| 亚洲最大成人手机在线| 免费无遮挡裸体视频| 国产精品精品国产色婷婷| 亚洲在线观看片| 国产高清激情床上av| 亚洲一区高清亚洲精品| 久久热精品热| 给我免费播放毛片高清在线观看| 日本熟妇午夜| 久久久精品欧美日韩精品| 人妻丰满熟妇av一区二区三区| 在线观看一区二区三区| 国产亚洲欧美98| 中文字幕av成人在线电影| 免费观看的影片在线观看| 一个人看的www免费观看视频| 在线观看美女被高潮喷水网站 | 伊人久久精品亚洲午夜| 成人特级av手机在线观看| 天堂网av新在线| 国产精品,欧美在线| 久久久久久久久久久丰满| 欧美日韩精品成人综合77777| 麻豆乱淫一区二区| 男的添女的下面高潮视频| 国产亚洲精品久久久com| av天堂中文字幕网| 亚洲欧美日韩卡通动漫| 91狼人影院| 狂野欧美激情性xxxx在线观看| 国产探花在线观看一区二区| 亚洲一级一片aⅴ在线观看| 99久久人妻综合|