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

    Role of biochar surface characteristics in the adsorption of aromatic compounds:Pore structure and functional groups

    2021-12-29 05:46:42XueFeiTnShiShuZhuRuPengWngYiDiChenPuLokeShowFengZhngShihHsinHo
    Chinese Chemical Letters 2021年10期

    Xue-Fei Tn,Shi-Shu Zhu,Ru-Peng Wng,Yi-Di Chen,e,Pu-Loke Show,Feng-F Zhng,Shih-Hsin Ho,?

    a College of Materials and Chemical Engineering,Heilongjiang Institute of Technology,Harbin 150050,China

    b State Key Laboratory of Urban Water Resource and Environment,School of Environment,Harbin Institute of Technology,Harbin 150090,China

    c School of Environmental Science and Engineering,Sun Yat-sen University,Guangzhou 510275,China

    d Guangdong Provincial Key Laboratory of Environmental Pollution Control and Remediation Technology,Sun Yat-sen University,Guangzhou 510275,China

    e State Key Laboratory of Urban Water Resource and Environment,School of Civil and Environmental Engineering,Harbin Institute of Technology(Shenzhen),Shenzhen 518055,China

    f Department of Chemical and Environmental Engineering,Faculty of Science and Engineering,University of Nottingham Malaysia,43500 Semenyih,Selangor,Malaysia

    Keywords:Biochar Surface characteristics Aromatics Adsorption Surface functional groups

    ABSTRACT Biochar (BC) are widely used as highly efficient adsorbents to alleviate aromatics-based contaminants due to their ease of preparation,wide availability,and high sustainability.The surface properties of BCs usually vary greatly due to their complex chemical constituents and different preparation processes and are reflected in the values of parameters such as the specific surface area (SSA),pore volume/size,and surface functional groups (SFGs).The effects of SSA and pore volume/size on the adsorption of aromatics have been widely reported.However,the corresponding mechanisms of BC SFGs towards aromatics adsorption remains unclear as the compositions of the SFGs are usually complex and hard to determine.To address in this gap in the literature,this review introduces a new perspective on the adsorption mechanisms of aromatics.Through collecting previously-reported results,the parameters logP (logarithm of the Kow),polar surface area,and the positive/negative charges were carefully calculated using ChemDraw 3D,which allowed the hydrophobicity/hydrophilicity properties,electron donor-acceptor interactions,Hbonding,and electrostatic interactions between SFGs and aromatics-based contaminates to be inferred intuitively.These predictions were consistent with the reported results and showed that tailor-made BCs can be designed according to the molecular weights,chemical structures,and polarities of the target aromatics.Overall,this review provides new insight into predicting the physicochemical properties of BCs through revealing the relationship between SFGs and adsorbates,which may provide useful guidance for the preparing of highly-efficient,functional BCs for the adsorption of aromatics.

    1.Introductio n

    With the increasing development of industry,large amounts of aromatics are discharged into ecological surroundings as a result of environmental migration and transformation,and exhibit high toxicity towards both humans and animals [1,2].Thus,the deployment of highly-efficient technology to eliminate aromatics present in the environment is an urgent need.Over the decades,adsorption has been regarded as the preferred method for removing organic contaminants in wastewater or groundwater [3],and various benchmark carbonaceous adsorbents such as activated carbon(AC) [4],carbon nanotubes (CNT) [5],graphene oxides (GO) [6],and biochar (BC) [7,8]have been widely used with satisfactory removal ability.Of these,BC in particular has been identified as a promising candidate for alleviating aromatics-based contaminates,owing to its relative ease of preparation,wide availability,and functional structure [9,10].Recent studies have identified many aromatics,which,along with their derivatives,can be rapidly adsorbed by BC,such as polycyclic aromatic hydrocarbons (PAHs),chlorinated pesticide,antibiotics,and polychlorinated biphenyl [11,12].However,as the adsorption mechanisms of aromatics on BC are still unclear,the production of tailor-made BCs designed to target specific aromatics is currently lacking,which hinders its development as a highly efficient adsorbent.

    Fig.1.(a-c) The adsorption mechanisms between biochar and aromatic compounds.

    Typically,the surface structural properties of BCs,such as surface functional groups (SFGs),specific surface area (SSA),and pore volume/size,are mainly responsible for the adsorption of aromatics from the aqueous phase,as opposed to the molecular weights,chemical structures,and polarities of the target aromatics [13].In particular,for BCs with higher pyrolysis temperatures,surface adsorption usually emerges as the dominant domain rather than the partitioning of matrix diffusion [14].Unlike in diffusion-limiting adsorption,the surface properties of BCs regulate the interfacial bonding by means of covalent bonding,hydrogen (H-) bonding,or electron donor-acceptor (EDA) interactions [15,16].Moreover,the pyrolysis temperature can alter the surface structure and phase of BCs.As the pyrolysis temperature increases from 100°C to 800°C,the surface pore sizes of BCs evolve from micropores to macropores,while moieties as functional groups gradually emerge on the defective edges of the BCs between 100 and 500°C and are eliminated above 600°C [17,18].In addition,the types of functional groups (e.g.,oxygen-containing) appearing on the BC surface can directly affect its adsorptive properties [19].Thus,gaining an indepth understanding of these surface properties,which are key factors affecting the adsorption of aromatics on BCs,is of great significance.

    The impact of the nature and characteristics of BC surfaces have been widely investigated through exploring the collinear relationship between surface properties and the adsorption of aromatics[20,21].The underlying mechanisms of the interactions between aromatics and surface adsorption sites are summarized in Fig.1,which suggests that the hydrophobic effect (weak hydrophobic interaction between solute and adsorbent surface),H-bonding (noncovalent specific interaction involving the aromaticπ-system),and EDA interactions (between the dissociated SFGs and the carbonaceous layers of BCs with electronegative atoms) as shown in Fig.1a,as well as coulombic forces,correspond to the surface functional groups and pore size or volume.Most studies focused on the large specific surface area of biochar,which can provide more active sites for adsorption (Fig.1b).And the larger pore volume can only act as a channel,while the smaller pore volume will prevent molecules from entering the pore channel to reduce the probability of pollutants interacting with biochar.Only a suitable pore volume can promote the occurrence of adsorption reactions (Fig.1c).However,to date,a comprehensive review focused on discussing and comparing the relationships between the surface properties of BCs and the adsorption of aromatics,as well as the underlying mechanisms,is still lacking.Therefore,to address this gap in the literature,the current paper aims to further enhance the applicability of BCs towards aromatics removal through critically reviewing,discussing,and comparing the relationships between the surface properties of BCs and the adsorption of aromatics.Additionally,and more importantly,the potential mechanisms by which the surface structures of BCs facilitate aromatics removal are also comprehensively revealed and summarized based on molecular computing analysis of data obtained from the literature using the chemical structure analysis tool ChemDraw 3D.

    2.Specific surface area and pore volume

    The porosity of BCs has been identified as a crucial factor in determining the adsorption of aromatics.It has also been found that the toluene adsorption capacity is linearly correlated to the BC specific surface area (SSA),with a correlation coefficient of 0.897 [22],suggesting that higher surface areas may provide more available sites for aromatics [23,24].Thus,improving the SSA is regarded as a feasible strategy for further increasing aromatics adsorption on BCs.Several hole-making methods have been used as activation approaches for enlarging the SSA [25].In particular,physical activation tend to result in microporous structures and larger SSAs.For instance,increasing the pyrolysis temperature could accelerate the development of porosity [26].As the temperature increases,the unstable components in biomass such as starch,hemicellulose,and lignin tend to decompose into volatile compounds such as acetic acid,tar,methanol,CO2,CO and H2.The overflow of volatile compounds can turn on the blocked pore channels to form more porous structures and produce new active sites [26,27].Moreover,condensation of the graphene layer distributed on the BC would occur,leading to the emergence of microcrystalline structures.The process of creating the random stacking graphite clusters also contributes to the enlargement of pore sizes and SSA [28].

    In addition to physical activation,chemical activation methods have also been widely employed to increase the porosity of BCs.The co-pyrolysis of a dehydrating agent or a corrosive activator with biomass causes the carbon layer to undergo a series of crosslinking condensation reactions or carving operations to form microporous structures [29].Liet al.[29]demonstrated that impregnating biomass with chemical activators could significantly develop the porosity and SSA and pore volume of BCs.Strong acids,alkalis,and salts,including H2SO4[29],H3PO4[30],KOH [31]and NaOH [32],have been identified as efficient hole-making activators;however,different chemical activation methods may result in different degrees of SSA and pore volume of BCs.For instance,it has been shown that KOH as an activator is beneficial for forming micropores,whereas H3PO4is more favorable towards the creation of microporous or mesoporous,as shown in Table S1 (Supporting information) [31,33].Moreover,apart from the type of species used for chemical activation,the incubation ratio and reaction period also have a significant impact on the development of SSA and pore structure as well as the on the ability of aromatics to adsorb on BCs.

    It is also important to establish the relationships between the SSA and pore volume of BCs and the adsorption of aromatic contaminants.Zhuet al.[34]utilized ZnCl2to enlarge the pore volume of BC from 0.08 cm3/g to 1.68 cm3/g,thereby improving its adsorption capacity towards toluene.To demonstrate this effect,they divided the pore volumes of BC into four size ranges:0.2-1.6 cm3/g (Fig.2a),0.2-0.9 cm3/g (Fig.2b),0.2-0.6 cm3/g (Fig.2c),0.08 to 0.2 cm3/g (Fig.2d).Their analysis showed that pore volumes ranging from 0.2 cm3/g to 0.9 cm3/g were best correlated with adsorption capacity,clearly indicating that the development of appropriate pore volumes can enhance the adsorption of aromatics on BCs.This phenomenon may result from the intragranular diffusion between suitable pore volumes and toluene during the adsorption process through van der Waals forces.Simultaneously,capillary condensation occurs in the mesopores,which is also an important factor affecting the adsorption efficiency [25].

    Fig.2.The effect of pores volume on toluene adsorption.Pore volume is (a) 0.2-1.6 cm3/g,(b) 0.2-0.9 cm3/g,(c) 0.2-0.6 cm3/g,(d) 0.06-0.2 cm3/g,respectively.Modified with permission [34].Copyright 2018,Elsevier.

    Fig.3.The effect of pore volume/size on adsorption.

    Furthermore,the pore sizes of BCs also have noticeable impacts on their adsorption ability.The pores of BCs can be generally divided into macropores (>50 nm),mesopores (2-50 nm),micropores (<2 nm),and narrow micropores (<1 nm) [35].By comparison,the molecular sizes of most aromatics are much smaller than the pore sizes of most BCs,where the varied pore size leads to different adsorption mechanisms.Crespoet al.[36]illustrated that pore sizes of 12 and 16.8 nm on carbonaceous materials exhibited distinct adsorption capacities towards benzene and thiophene.Stronger adsorption was observed for the smaller carbonaceous.However,the presence of pore sizes far larger than the adsorbate molecular diameters would cause the interactive forces between the adsorbent and the adsorbate to be largely eliminated,resulting in the pores acting mainly as diffusion channels.Wanget al.[37]synthesized mesoporous carbon materials for benzene and cyclohexane adsorption,demonstrating that the diffusion coefficient of the synthesized mesoporous carbon with 5.0 nm pores was twice that of the material with 1.8 nm pores.This finding further confirms that pores which are slightly larger in size than the adsorbate can provide adsorption sites,whereas oversized pores likely enhance the diffusion of aromatics.The adsorption mechanisms of aromatics on BCs with different pore sizes are illustrated in Fig.3.

    3.Surface functional groups

    It is worth noting that although the functional groups on the BC surface largely affect the chemical reactions between BCs and aromatics,it is still not straightforward to clearly identify and distinguish all of the SFGs.To date,however,several technologies,such as FT-IR,X-ray photoelectron spectroscopy (XPS),and nearedge X-ray absorption fine structure spectroscopy (NEXAFS) have been applied towards investigating these functional groups.Based on the biomass composition,and the different pretreatment,activation,and pyrolysis conditions,a given BC would possess various functional groups distributed on its surface,which greatly affect the adsorption properties.

    The SFGs existing on BC are commonly classified as either oxygen-,nitrogen-,or sulfur-containing,respectively [38].Of these,the oxygen-containing functional groups are considered the most common,and can be further divided into acidic,neutral,and basic groups based on their properties.The acidic groups normally include the carboxyl,lactonic and phenolic groups.On the other hand,the pyrone-and chromene-type groups are accepted as two common basic active sites,whose properties can be attributed to the basic nature of the carbon surface,which is more pronounced in carbon atoms without oxygen due to the presence of delocalized electrons [39].Many studies have reported that the adsorption capacity of biochar adsorbents has increased after oxidation modification [40,41].The introduction of nitrogen-containing functional groups can significantly enhance the polarity of the carbon surface,thereby increasing the specific interaction between pollutants and polar adsorbents [42].Whereas,sulfur-containing functional groups can adjust the adsorption capacity of biochar by changing its specific surface area and pore volume [43,44].

    During the adsorption of aromatics,and in particular ionic aromatic pollutants such as herbicides and antibiotics,various forces occur between the adsorbate and SFGs [45].The SFGs with different physical and chemical properties interact with the adsorbent and the liquid phase depending on their typical structural characteristics.It is generally accepted that the adsorption process of the aromatics mainly depends on their chemical interactions with the multifarious SFGs of BC adsorbents [45].The distribution,polarity,and hydrophobicity of SFGs greatly influence the adsorption capacity of aromaticsviaH-bonding,surface charges,and other interactions [45].However,it is difficult to be intuitively observed through experiments.According to the spatial structure of SFGs,parameters such as logP(logarithm of theKow),PSA (polar surface area),positive and negative charges are calculated by ChemDraw 3D to measure their hydrophobicity/hydrophilicity,H-bonding,and electron donating/accepting properties,as shown in Table 1.Therefore,new insights into the relationship between SFGs and adsorbates should be provided to enable better understanding of the properties of the functional groups,paving the way for the design of“smart”surfaces with enhanced surface chemistry for achieving the goal of highly-selective adsorption of target contaminants.

    Table 1Structure,charge,logP and polar surface area of SFGs on biochar.

    3.1.Hydrophobicity/hydrophilicity of biochar

    The mechanism of biochar adsorption of aromatic hydrocarbons usually depends on the surface properties of the adsorbent,the structure of the organic adsorbate and environmental conditions.Ohet al.[46]confirmed that the hydration of SFGs could greatly hamper the process of phenol and nitrobenzene adsorption by blocking the inherent active sites.In the case of BC,the hydrophobicity/hydrophilicity has also been used to predict the relationship between the adsorption efficiency and SFGs [47].The polarity index (i.e.,the ratio (O+N)/C or O/C) is commonly recognized as an indicator of adsorbability [48],and it is often found that the hydrophobic interactions play an important role in the adsorptionof aromatics onto BC [49].However,using the polarity index exclusively to predict the hydrophobicity/hydrophilicity of BC by the bulk elements may be inadequate,as the adsorption process usually occurs between the organic adsorbate and the adsorbent SFGs.Thus,it is necessary to determine the surface polarity,which has been employed as an alternative predictor,using XPS.

    Recently,the impact of the hydrophobicity/hydrophilicity of SFGs on adsorption was also proposed as way to predict its adsorption capacity through calculating its organic solvent-water partitioning coefficient (Kow) for octanol.However,althoughKowis related to the surface polarity of the adsorbent,it is hard to specify the polarity of each SFG on BC due to interference from surrounding factors such as ionic strength and temperature [50].For BCs in particular,various atoms (e.g.,C,H,O,N,P and S) contributing to the BC structure can form different combinations of SFGsviathe formation of hydroxyl,epoxy,carboxyl,acyl,carbonyl,ether,ester,sulfonic,pyrrole,pyridine,amine,imine,acylamide,nitroso and nitro groups on the surface [51,52].Thus,calculating the pKa(basic dissociation constants) values and concentrations of SFGs by acidbase titration is likely more effective in enabling the polarity ofthe BC surface as well as its contribution towards the adsorption capacity of aromatics to be predicted [51].Hence,in this review,we demonstrate how to predict the hydrophobicity/hydrophilicity of the BC surface by calculating the partition coefficient (logP) using ChemDraw 3D,based on the molecular orbital theory and the default molecular mechanics (MM2) (Table 1).

    The logPparameter represents logarithm of theKowof neutral species [21].Since the presence of SFGs on BC would modulate the adsorption mechanisms of aromatics due to the varied hydrophobicity,using logPcould enable prediction of the corresponding adsorption variability information.LogPis driven by the free energy of the transfer from aqueous phase to nonpolar phase of the aromatics [21,53].In this vein,several well-established mathematical correlations between chemical structure and logPhave been developed [54,55];however,there is barely any literature to date on the relationship between logPand the adsorption of aromatics on BC.According to previously-reported results [51],the common SFGs on the BC surface are shown in Table 1.The authors of the current paper hypothesized that the SFGs with different properties and pKavalues would significantly impact the logPvalues of the BC surface.The largest differences in logPwere seen in the oxygen-containing groups with varied carbon configurations and spatial distributions,such as structures 1,2,3,4 and 10 in Table 1.The merging of nitrogen and sulfur atoms with oxygencontaining groups caused the hydrophilicity of the BC to increase,with resulting logPvalues of -0.28,-0.99 and -0.87 for the nitrogen and sulfur-containing structures 5,8,and 9,respectively,and 0.49 and 0.52 for the phosphorus-containing structures 6 and 7,respectively.As a result,using logPas an indicator of hydrophobicity/hydrophilicity to reliably predict the adsorption of aromatics on BC seems a promising approach,and bears further investigation.

    3.2.Electron donor-acceptor interactions

    Electron donor-acceptor (EDA) interactions are recognized as an important type of noncovalent specific interaction involving the aromaticπ-system.Due to the graphene units often co-existing with carbonaceous adsorbents,π-EDA interactions on BC are likely to occur in the presence of ionic aromatic adsorbates [50].Generally,the hydroxyl groups on the surface of carbonaceous adsorbents can act as electron donors,while the carboxyl functional group may act as acceptors [56].

    Xieet al.demonstrated thatπ-πEDA interactions occurred between sulfamethoxazole and the aromaticity of BCs [57].A structural analysis of the sulfamethoxazole showed that it contains an amino functional group,an N-heteroaromatic ring,and an unprotonated sulfonamide group with high electron withdrawing ability,all of which can act as strongπ-electron acceptors [58].Moreover,it has been demonstrated thatπ-πEDA interaction acted as an important role between a series of heteraromatic amines and graphite was also reported [47,57].

    In addition,theπ-πand n-πEDA interactions of anion contaminants may be increased due to the dissociation of -OH groups which are stronger electron donors compared to neutral species[59].For instance,Wuet al.[60]found that the removal efficiency of naphthol using carbon nanotubes (CNTs) and graphite as adsorbents was improved,with enhanced EDA interactions,at higher pH values.In general,compared to cations and zwitterions,the adsorption capacities of anions on carbonaceous adsorbents are lower,while theirπ-πand n-πEDA interaction intensities are also relatively weak.Additionally,the negative repulsions occurring between the adsorbate and adsorbent due to acid groups [50].Many aromatic amines are positively charged in a specific pH range,enabling the charges of cations located on the aromatic ring to interact with the electron-rich carbonaceous materialvia π+-πEDA interactions [61].Recently,ten aromatic cations (including triazine herbicides,heterocyclic aromatic amines,substituted anilines,and other amines) were reacted with pyrogenic carbonaceous materials(graphite and BCs),which clarified the function of theπ+-πEDA interactions,thus confirming that theπ-electron acceptors provided by the quinoline cations could interact with graphite [62].Further,based on the density function theory (DFT),electrons can be transferred to the ring system during adsorption,which may contribute to the enhancement of cation-πinteractions [62].

    On the other hand,the presence of zwitterions in amphoteric pollutants may make it difficult to obtain a detailed understanding of the different types of EDA interactions occurring and their interplay during the adsorption process.For instance,bothπ-πand cation-πEDA interactions may occur between tetracycline and the graphene surface [63].It has been suggested that the varied EDA interactions play different roles between the carbonaceous adsorbents and the SFGs of amphoteric compounds (e.g.,quinolone groups in antibiotics) [56,64].

    SFGs as the electron donor/acceptor part greatly affect the electrochemical behavior of the material,the charges of the SFGs were calculated and are shown in Fig.4 and Table 1.By comparison,the corresponding oxygen-,phosphorous-,sulfur-and nitrogencontaining functional groups possess obvious charge characteristics.For instance,the charges of O(2) in structure 2,O(5) in structure 3,and O(4) in structure 4 (Fig.4) are -0.288,-0.414,-0.531,respectively,resulting in their negative charges enhancing the donor strength of the electrons in the oxygen-containing functional groups.Klüpfelet al.also confirmed that the phenolic-OH and -C-O groups can be considered as potential electron donating capacity groups [65].However,carboxyl -COOH,phenolic-OH and electron-rich ketones groups (-C-O),generally can be used as an electron donating group [66,67].The establishment of N-containing and S-containing functional groups can increase the electroactive sites of BC,which is an effective method to promote the adsorption reaction.As we all know,N atoms are rich in electrons,which can adjust the electronic properties of BC and help improve the conductivity of the carbon matrix [68,69].Because of this,the introduction of nitrogen-containing functional groups cannot only increase the adsorption rate of BC,but also increase the adsorption capacity of aromatic compounds,especially pyridine N,pyrrole N and graphite N [70].Additionally,as the nitrogen atoms can serve as both electronic donors and receptors,the charge of N(3) in structure 5 is 0.236 (i.e.,positive) while other N atoms have negative charges,which suggests that N atoms located at different sites in nitrogenous functional groups play different roles in the adsorption process.Moreover,due to sulfur atoms exhibiting similar electronic receptor properties to phosphorus,positive charges were observed in structures 7,8 and 9.Furthermore,Du et al.found that the introduction of S-containing functional groups can accelerate the reaction [71].It has been reported that sulfurcontaining functional groups,especially thiophene-S,can help induce adjacent C sites to have a higher positive charge density [72-74].Thus,these results clearly demonstrate that the atom distribution and structure of SFGs can greatly affect the electron donating/accepting properties of carbonaceous adsorbents.

    3.3.Hydrogen bonding

    Another mechanism by which SFGs may have a significant impact on the adsorption of aromatics isviaH-bonding between the dissociated SFGs and the carbonaceous layers of BCs with electronegative atoms (e.g.,F,N,O) in organic compounds [75,76].These H-bonds could largely contribute to the adsorption of polar aromatics (e.g.,phenol and bisphenol),despite the existence of columbic forces with electrostatic repulsion [77].However,although H-bonds have been found to occur between on series ofadsorbents [49],intuitive evidence for the corresponding mechanism still remains elusive.

    Fig.4.The charge data of common SFGs on biochar.

    It has been reported that the polar surface area (PSA),simply defined as that part of a molecular surface corresponding to oxygen and nitrogen atoms (including the attached hydrogens),can be used as an index for measuring the H-bonding capacity [78].Another study also demonstrated a high correlation between PSA and the numbers of H-bond donors and acceptors [79].Thus,PSA has been applied as a useful indicator to describe the accessible area for H-bonding,and is positively correlated with the water solubility of the species as well [80].Table 1 shows that the oxygencontaining functional groups which also contain phosphorus and sulfur (6,7,8,9 in Fig.4) have a larger PSA,which suggests that such groups could provide more H-bonding interactions between the carbonaceous adsorbent surface and aromatics.In contrast,the basic groups containing C-C and C=C on the surface of carbon material are unable to provide hydrogen bonds as their PSA is close to zero.Additionally,another type of H-bonding formed by charge strengthening (a.k.a.charge assisted hydrogen bonds -CAHBs) also largely influences the interaction properties between the ionic organic pollutants and SFGs (N,O,and S,including ester,carbonyl,nitrile,and thiol groups) during the adsorption process.

    CAHBs usually have much higher bond energy than ordinary hydrogen bonds [81].Recent studies have pointed out that(-)CAHB can be used to explain the anomalous adsorption between surface functionalized adsorbents and ionic organic compounds[82,83].The following (-)CAHB mechanism has been proposed for carboxyl groups on biochar and CNTs with increasing pH (Eq.1)[50]:

    When the acid constant (pKa) of the adsorbent containing O or N SFGs is close to that of the adsorbate,the effect of (-)CAHB on carbon nanotubes is stronger [81].

    During the adsorption of ionized sulfadimidine by BC,Leeet al.found that CAHBs are an alternative route to van der Waals forces and hydrophobicity for promoting adsorption [84].It has also been reported that the anions produced by the dissociation of sulfamethoxazole can interact with weak acid groups (such as acidic oxyl and carboxyl groups) in BCs to form CAHB,and the CAHB interactions increased with increasing pKafor pollutants closer to the BC surface [85].Furthermore,under neutral and basic conditions,the (-)CAHB interactions can be strengthened to further enhance adsorption by dissociation of pollutants [84].Additionally,as shown in Fig.4 and Table 1,the O atom seems to display higher electronegativity than other atoms within the oxygen-containing functional groups.In structures 6,7,8 and 9 (Fig.4) in particular,the electronegativity of O atoms was noticeably enhanced when P and S atoms were present,suggesting that these atoms are electropositive within oxygen-containing functional groups.

    3.4.Electrostatic interaction

    Electrostatic interactions require adsorbents with positive or negative charges.Anion exchange by electrostatic attraction has recently been reported as an important contributor to the adsorption of aromatics,particularly for functionalized carbonaceous materials with amino groups [86].In the data analysis presented in this pa-per,it was found that the N atoms at different sites of the nitrogenous functional groups usually owned negative charges (with only structure 5 in Fig.4 showing a positive charge),demonstrating that most of the nitrogen-containing groups would be unlikely to promote electrostatic attraction by anion exchange.On the other hand,it was found that the P atoms in the phosphorus-oxygen groups and the S atoms in the thio-oxygen groups owned some positive charges,which were 2.08 and 2.09 in structures 6 and 7 and 2.63 and 2.59 in structures 8 and 9,respectively (Fig.4 and Table 1).Based on these finds,it can be speculated that the phosphorusoxygen and thio-oxygen functional groups play the dominant role in strengthening the anion exchange;however,as most carbonaceous adsorbents are rich in negative charges with plentiful acid functional groups,the function of anion exchange is limited.Thus,cation exchange is likely to contribute more electrostatic interactions than anion exchange on most of the carbonaceous adsorbents.

    The carbonaceous adsorbents with various SFGs usually display a superior cation exchange capacity towards positively-charged pollutants when their pKa>PZC (the point of zero charge).Additionally,the cation exchange may be further improved by the addition of acidic groups on the adsorbent surface.For instance,Maet al.confirmed that the carbon nanotubes loaded with acidic groups can increase the cation exchange capacity and show better adsorption capacity [87].BC is also usually recognized as an excellent adsorbent due to its abundant SFGs capable of improving the cation exchange with positively-charged pollutants [88].Additionally,all the O atoms in the oxygen-containing functional groups can carry more negative charges (Fig.4 and Table 1),clearly suggesting that these groups greatly contribute to the electrostatic attraction by cation exchange.Moreover,in addition to electrostatic attractions,the adsorption of pollutants also contains electrostatic repulsions.However,although few studies to date have highlighted the function of electrostatic repulsion or mentioned the positive correlation between the repulsion force and the degree of oxidation of carbonaceous material [89,90],most of the electrostatic repulsions which occur between carbonaceous materials and aromatic pollutants are negligible when compared with the electrostatic attractions.

    4.Conclusions

    The main conclusions derived from the literature review and data analysis presented in this paper are as follows:(1) In interactions between BC adsorbents and aromatic adsorbates,the presence of adsorbent pore sizes which are only slightly larger than the molecular diameters of the adsorbates contributes to better adsorption capacity,as much larger pore sizes would act as diffusion channels due to weaker adsorptive forces.(2) Through calculation of certain parameters (such as logP,PSA and positive/negative charges),the mechanisms of different kinds of SFGs in the adsorption process may be analyzed.In particular,the value of logPcan act as a predictor of the hydrophobicity/hydrophilicity on the BC surface,while the introduction of PSA allowed the H-bonding accessible area to be positively correlated with the H-bonding capacity.(3) Based on the atomic distributions and structures of SFGs,the positive/negative charges present could also be useful indicators to reflect the electron donating/accepting properties of BC.Additionally,the amount of positive/negative charges are likely responsible for the strength of electrostatic interactions.Through using these parameters calculated in this paper,the interactions between the aromatic adsorbents and the SFGs of BCs can be intuitively revealed,thus enabling both the preparation of tailor-made BCs as well as a better understanding of the physicochemical properties of aromatics.

    It is hoped that this review will bring new insight into an indepth understanding of the roles of the main surface characteristics of BCs towards the adsorption of aromatic compounds,thereby providing new information to guide the directional design of BCs,further enhancing their adsorption efficiency.

    Declaration of competing interest

    The authors declare no conflict of interest.

    Acknowledgment

    This study was funded by the National Natural Science Foundation of China (No.51961165104),Project of Thousand Youth Talents (No.AUGA2160100917),University Nursing Program for Young Scholars with Creative Talents in Heilongjiang Province (No.UNPYSCT-2020050) and Provincial Leading Talent Echelon Cultivation Project of Heilongjiang Institute of Technology (No.2020LJ04).

    Appendix A.Supplementary data

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

    国产成人欧美在线观看| 又大又爽又粗| 91老司机精品| av超薄肉色丝袜交足视频| 欧美激情 高清一区二区三区| 美女 人体艺术 gogo| 国产三级黄色录像| 国产又爽黄色视频| 美女扒开内裤让男人捅视频| 老汉色av国产亚洲站长工具| 老司机靠b影院| 亚洲精品国产区一区二| 可以在线观看的亚洲视频| 91成人精品电影| 精品国产超薄肉色丝袜足j| netflix在线观看网站| 亚洲欧美精品综合久久99| 日日干狠狠操夜夜爽| 日韩欧美一区视频在线观看| 久久精品91蜜桃| 亚洲最大成人中文| 国产精品永久免费网站| tocl精华| 欧美zozozo另类| 亚洲男人天堂网一区| 在线视频色国产色| 嫩草影视91久久| 好看av亚洲va欧美ⅴa在| 成人免费观看视频高清| 国产av在哪里看| 国产精品亚洲一级av第二区| 伦理电影免费视频| 亚洲精品久久国产高清桃花| 国产精品影院久久| 黄网站色视频无遮挡免费观看| 亚洲一区中文字幕在线| xxx96com| 国内精品久久久久久久电影| a在线观看视频网站| 满18在线观看网站| 久久国产精品男人的天堂亚洲| 亚洲国产日韩欧美精品在线观看 | 老熟妇乱子伦视频在线观看| 亚洲第一欧美日韩一区二区三区| 亚洲国产精品成人综合色| 亚洲欧美日韩无卡精品| 成人国产综合亚洲| 91大片在线观看| 美女高潮到喷水免费观看| 日韩免费av在线播放| 91老司机精品| 欧美成人免费av一区二区三区| 欧美丝袜亚洲另类 | 亚洲五月婷婷丁香| 亚洲国产欧洲综合997久久, | 在线av久久热| 麻豆成人午夜福利视频| 丝袜美腿诱惑在线| 国产午夜精品久久久久久| 哪里可以看免费的av片| 啦啦啦免费观看视频1| 又黄又粗又硬又大视频| 啦啦啦韩国在线观看视频| 色综合欧美亚洲国产小说| 一级片免费观看大全| 美女午夜性视频免费| 老司机靠b影院| 亚洲国产欧美日韩在线播放| 国产一区二区三区视频了| 中文字幕高清在线视频| 国产激情久久老熟女| 可以在线观看的亚洲视频| 一本久久中文字幕| 成人亚洲精品一区在线观看| 欧美黄色片欧美黄色片| 美女免费视频网站| 每晚都被弄得嗷嗷叫到高潮| 欧美不卡视频在线免费观看 | 麻豆国产av国片精品| 最近在线观看免费完整版| 麻豆av在线久日| 美女午夜性视频免费| 国产av在哪里看| 亚洲成人精品中文字幕电影| 女人高潮潮喷娇喘18禁视频| 国产亚洲精品一区二区www| 久久久国产成人精品二区| 99在线视频只有这里精品首页| 亚洲色图 男人天堂 中文字幕| 亚洲av成人不卡在线观看播放网| 免费在线观看完整版高清| 久久婷婷人人爽人人干人人爱| 男女做爰动态图高潮gif福利片| 一进一出好大好爽视频| 久热这里只有精品99| 激情在线观看视频在线高清| 国产精品亚洲美女久久久| 国产亚洲精品av在线| 在线观看舔阴道视频| 欧美黄色淫秽网站| 精品熟女少妇八av免费久了| tocl精华| 后天国语完整版免费观看| 国产一区二区三区在线臀色熟女| 免费电影在线观看免费观看| 国产成人系列免费观看| 午夜福利高清视频| 夜夜看夜夜爽夜夜摸| 在线十欧美十亚洲十日本专区| 一区二区三区国产精品乱码| 欧美大码av| 99热这里只有精品一区 | 亚洲第一电影网av| 人妻久久中文字幕网| x7x7x7水蜜桃| 久久香蕉激情| 十八禁人妻一区二区| 成人午夜高清在线视频 | 精品国产乱子伦一区二区三区| 日韩欧美在线二视频| 久久精品国产99精品国产亚洲性色| av视频在线观看入口| 国产三级黄色录像| 久久久久久久午夜电影| 久久人妻福利社区极品人妻图片| 丁香欧美五月| 久久人人精品亚洲av| 精品国产亚洲在线| 国产爱豆传媒在线观看 | 欧美成人午夜精品| 精品不卡国产一区二区三区| 久久婷婷成人综合色麻豆| 天天躁夜夜躁狠狠躁躁| 2021天堂中文幕一二区在线观 | 此物有八面人人有两片| 一边摸一边做爽爽视频免费| 欧美乱色亚洲激情| 国产aⅴ精品一区二区三区波| 欧美日韩亚洲国产一区二区在线观看| 国产精品久久电影中文字幕| 亚洲精品粉嫩美女一区| 亚洲激情在线av| 国产一卡二卡三卡精品| 女生性感内裤真人,穿戴方法视频| 国产伦人伦偷精品视频| 亚洲精品中文字幕一二三四区| 啪啪无遮挡十八禁网站| 国产精品美女特级片免费视频播放器 | 99热只有精品国产| 亚洲国产精品999在线| 精品人妻1区二区| 神马国产精品三级电影在线观看 | 欧美性长视频在线观看| 搡老熟女国产l中国老女人| 欧美 亚洲 国产 日韩一| 啦啦啦观看免费观看视频高清| 国产亚洲精品久久久久5区| 日韩欧美免费精品| 色精品久久人妻99蜜桃| 成人手机av| 欧美在线一区亚洲| 亚洲av电影在线进入| 丰满的人妻完整版| 一区二区日韩欧美中文字幕| 中文字幕人妻丝袜一区二区| 欧美性猛交黑人性爽| 日韩视频一区二区在线观看| 两个人视频免费观看高清| 亚洲一码二码三码区别大吗| 麻豆av在线久日| 久久香蕉激情| 久久精品91无色码中文字幕| 视频在线观看一区二区三区| 午夜免费激情av| 久久这里只有精品19| 欧美色视频一区免费| 视频在线观看一区二区三区| 久久久久久久精品吃奶| 午夜久久久久精精品| 中文字幕最新亚洲高清| 久久国产精品人妻蜜桃| 在线播放国产精品三级| 国产人伦9x9x在线观看| 日日夜夜操网爽| 黄色毛片三级朝国网站| 中文资源天堂在线| 老司机在亚洲福利影院| 久久久久久人人人人人| 麻豆久久精品国产亚洲av| 欧美亚洲日本最大视频资源| 国产伦在线观看视频一区| 亚洲熟女毛片儿| 欧美又色又爽又黄视频| 免费电影在线观看免费观看| 男女床上黄色一级片免费看| 久久狼人影院| 国产三级黄色录像| 搡老妇女老女人老熟妇| 亚洲五月天丁香| 香蕉丝袜av| 日韩欧美三级三区| 哪里可以看免费的av片| 大型av网站在线播放| 免费人成视频x8x8入口观看| 国产一区二区激情短视频| 亚洲av美国av| 成人亚洲精品一区在线观看| 这个男人来自地球电影免费观看| 人人妻人人澡人人看| 国产视频内射| 国产成人av激情在线播放| 一级毛片高清免费大全| 国产精品乱码一区二三区的特点| 国产一级毛片七仙女欲春2 | 亚洲精品美女久久av网站| 久久天躁狠狠躁夜夜2o2o| 露出奶头的视频| 18禁裸乳无遮挡免费网站照片 | 欧美色视频一区免费| 人人妻人人澡人人看| 国产亚洲欧美精品永久| 看免费av毛片| 搞女人的毛片| 国产伦一二天堂av在线观看| 91在线观看av| 欧美三级亚洲精品| 国产视频一区二区在线看| 黑人巨大精品欧美一区二区mp4| 精品一区二区三区视频在线观看免费| www日本在线高清视频| 老司机靠b影院| 亚洲精品久久成人aⅴ小说| 男人舔女人下体高潮全视频| 最好的美女福利视频网| 精品国产一区二区三区四区第35| 美女免费视频网站| 99在线人妻在线中文字幕| 午夜福利在线观看吧| 亚洲人成伊人成综合网2020| 美女 人体艺术 gogo| 亚洲国产中文字幕在线视频| netflix在线观看网站| 久久久久久免费高清国产稀缺| 国产主播在线观看一区二区| 国产野战对白在线观看| 欧美午夜高清在线| 午夜福利免费观看在线| 亚洲五月天丁香| 亚洲免费av在线视频| 国产又黄又爽又无遮挡在线| netflix在线观看网站| 少妇 在线观看| 自线自在国产av| 国产真人三级小视频在线观看| 久久中文看片网| 俄罗斯特黄特色一大片| www.www免费av| 国产视频内射| 日韩欧美国产在线观看| 超碰成人久久| 国产精品国产高清国产av| 悠悠久久av| 草草在线视频免费看| 国产伦人伦偷精品视频| 国产精品自产拍在线观看55亚洲| 精品久久久久久久末码| 午夜精品久久久久久毛片777| 91大片在线观看| 国语自产精品视频在线第100页| 女人高潮潮喷娇喘18禁视频| 国产私拍福利视频在线观看| 视频区欧美日本亚洲| 成人手机av| 国产片内射在线| 欧美不卡视频在线免费观看 | 国产精品久久电影中文字幕| 一级毛片精品| 777久久人妻少妇嫩草av网站| 女同久久另类99精品国产91| 亚洲欧美精品综合一区二区三区| ponron亚洲| 久久婷婷人人爽人人干人人爱| 夜夜爽天天搞| 久久性视频一级片| 丝袜人妻中文字幕| 精品午夜福利视频在线观看一区| 看免费av毛片| 国产精品99久久99久久久不卡| 好男人电影高清在线观看| 亚洲av美国av| a级毛片在线看网站| 欧美成人免费av一区二区三区| 成人永久免费在线观看视频| 久久久久久人人人人人| 国产精品免费一区二区三区在线| 视频在线观看一区二区三区| 日本在线视频免费播放| 国产精品国产高清国产av| 国产三级在线视频| 美女午夜性视频免费| 亚洲久久久国产精品| 麻豆av在线久日| www.熟女人妻精品国产| 在线观看免费视频日本深夜| 午夜福利欧美成人| 久久久久久国产a免费观看| 国产精品 国内视频| 色综合欧美亚洲国产小说| 麻豆一二三区av精品| 老熟妇仑乱视频hdxx| 国产亚洲欧美在线一区二区| 亚洲成a人片在线一区二区| 亚洲国产日韩欧美精品在线观看 | 可以在线观看毛片的网站| 精品国产超薄肉色丝袜足j| 精品国产一区二区三区四区第35| 一二三四社区在线视频社区8| 日韩欧美免费精品| 国产国语露脸激情在线看| 精品久久久久久久久久免费视频| 国产亚洲精品av在线| 精品国产亚洲在线| 亚洲国产精品sss在线观看| 久热这里只有精品99| 国产黄色小视频在线观看| 中文字幕另类日韩欧美亚洲嫩草| 丝袜人妻中文字幕| 在线av久久热| 国产激情欧美一区二区| 国产精品98久久久久久宅男小说| 亚洲专区中文字幕在线| 国产一区在线观看成人免费| 亚洲中文av在线| 在线观看免费日韩欧美大片| 丰满人妻熟妇乱又伦精品不卡| 精品久久久久久久人妻蜜臀av| 熟女电影av网| 男女下面进入的视频免费午夜 | 哪里可以看免费的av片| 亚洲精品久久国产高清桃花| 久久草成人影院| 亚洲专区字幕在线| 国产精品综合久久久久久久免费| 亚洲av电影在线进入| 欧洲精品卡2卡3卡4卡5卡区| 在线观看66精品国产| www日本黄色视频网| 19禁男女啪啪无遮挡网站| x7x7x7水蜜桃| 搡老熟女国产l中国老女人| 色在线成人网| 日韩欧美免费精品| 成人亚洲精品av一区二区| 国产乱人伦免费视频| 国产一区在线观看成人免费| 黑丝袜美女国产一区| 成年免费大片在线观看| av电影中文网址| 无遮挡黄片免费观看| 天天添夜夜摸| 可以在线观看的亚洲视频| 欧美一级a爱片免费观看看 | 两性午夜刺激爽爽歪歪视频在线观看 | 免费在线观看亚洲国产| 91在线观看av| 亚洲真实伦在线观看| 女性生殖器流出的白浆| 国产区一区二久久| 久久久久免费精品人妻一区二区 | 国产高清激情床上av| 亚洲狠狠婷婷综合久久图片| av片东京热男人的天堂| www日本在线高清视频| 久久99热这里只有精品18| 中亚洲国语对白在线视频| 一卡2卡三卡四卡精品乱码亚洲| 侵犯人妻中文字幕一二三四区| 国产一卡二卡三卡精品| 777久久人妻少妇嫩草av网站| 国产野战对白在线观看| 满18在线观看网站| 国产高清videossex| 久热爱精品视频在线9| 他把我摸到了高潮在线观看| 精品久久久久久久末码| 级片在线观看| 久久久久久免费高清国产稀缺| 男女床上黄色一级片免费看| 日本一区二区免费在线视频| 好看av亚洲va欧美ⅴa在| 在线观看一区二区三区| 视频在线观看一区二区三区| 级片在线观看| 999久久久精品免费观看国产| 国产一区二区激情短视频| 欧美成人免费av一区二区三区| 色婷婷久久久亚洲欧美| 亚洲国产毛片av蜜桃av| 精品免费久久久久久久清纯| 国产97色在线日韩免费| 欧美精品亚洲一区二区| 亚洲熟妇中文字幕五十中出| 午夜激情福利司机影院| 美女国产高潮福利片在线看| 国产av一区二区精品久久| 国产精华一区二区三区| 看片在线看免费视频| 午夜影院日韩av| 一边摸一边抽搐一进一小说| 成人欧美大片| 色在线成人网| 亚洲av熟女| 午夜久久久久精精品| 999精品在线视频| 成人亚洲精品一区在线观看| 成年免费大片在线观看| 欧美激情久久久久久爽电影| 免费av毛片视频| 真人一进一出gif抽搐免费| 成人一区二区视频在线观看| 一边摸一边做爽爽视频免费| 亚洲第一青青草原| 亚洲五月婷婷丁香| 啦啦啦韩国在线观看视频| 亚洲av成人不卡在线观看播放网| 免费观看人在逋| 人妻久久中文字幕网| 熟女电影av网| 精品久久久久久久毛片微露脸| 久久精品国产99精品国产亚洲性色| 在线观看66精品国产| 色综合站精品国产| 国产精品 欧美亚洲| 久久久久久人人人人人| 欧美成人性av电影在线观看| 色av中文字幕| 精品一区二区三区视频在线观看免费| 欧美+亚洲+日韩+国产| 国产91精品成人一区二区三区| 日韩一卡2卡3卡4卡2021年| 一二三四在线观看免费中文在| 国产一区二区在线av高清观看| av片东京热男人的天堂| 日本熟妇午夜| 免费av毛片视频| 在线观看免费午夜福利视频| 日韩国内少妇激情av| 天堂√8在线中文| 国产单亲对白刺激| 亚洲精品在线美女| 91成年电影在线观看| 日本三级黄在线观看| 国内少妇人妻偷人精品xxx网站 | 久久久久久久久免费视频了| 一级a爱片免费观看的视频| 国产爱豆传媒在线观看 | 琪琪午夜伦伦电影理论片6080| 一本精品99久久精品77| 美女大奶头视频| 无人区码免费观看不卡| 亚洲熟女毛片儿| 90打野战视频偷拍视频| 美女扒开内裤让男人捅视频| 亚洲精品在线观看二区| 亚洲黑人精品在线| 国产在线精品亚洲第一网站| 国产成人系列免费观看| 午夜精品久久久久久毛片777| 久久久久国产精品人妻aⅴ院| 国产午夜精品久久久久久| 高清在线国产一区| 亚洲真实伦在线观看| a在线观看视频网站| 精品久久久久久久人妻蜜臀av| 一级a爱片免费观看的视频| 一级作爱视频免费观看| 不卡av一区二区三区| 黄色视频不卡| 中文字幕最新亚洲高清| 巨乳人妻的诱惑在线观看| 国产av又大| 国产伦一二天堂av在线观看| 国产精品,欧美在线| 日本撒尿小便嘘嘘汇集6| 熟女电影av网| 久久精品国产亚洲av高清一级| 看片在线看免费视频| 黄片小视频在线播放| 美女扒开内裤让男人捅视频| 此物有八面人人有两片| 午夜福利在线观看吧| 久久久久免费精品人妻一区二区 | 久久精品国产亚洲av高清一级| 丁香六月欧美| 亚洲,欧美精品.| 午夜成年电影在线免费观看| 久久久久亚洲av毛片大全| 高清毛片免费观看视频网站| 亚洲午夜理论影院| 国产高清激情床上av| 免费在线观看成人毛片| 国产一区二区三区视频了| 国产精品 国内视频| 欧美精品啪啪一区二区三区| 国产一区二区激情短视频| 精品久久久久久成人av| 亚洲精品中文字幕在线视频| 大型av网站在线播放| 精品久久久久久久人妻蜜臀av| 久久婷婷人人爽人人干人人爱| 免费看日本二区| 亚洲全国av大片| 18禁裸乳无遮挡免费网站照片 | 不卡一级毛片| 国产精品久久久久久人妻精品电影| 亚洲av中文字字幕乱码综合 | 精品欧美一区二区三区在线| 国产爱豆传媒在线观看 | 亚洲一码二码三码区别大吗| 波多野结衣巨乳人妻| 国产私拍福利视频在线观看| 精品久久蜜臀av无| 91字幕亚洲| 久久性视频一级片| 婷婷亚洲欧美| 久久中文字幕一级| 国产精品电影一区二区三区| bbb黄色大片| 免费无遮挡裸体视频| 香蕉丝袜av| 国产激情偷乱视频一区二区| 丰满的人妻完整版| 午夜久久久在线观看| 搡老岳熟女国产| 亚洲免费av在线视频| 此物有八面人人有两片| 91国产中文字幕| 色婷婷久久久亚洲欧美| 亚洲无线在线观看| 亚洲 欧美一区二区三区| 亚洲成人国产一区在线观看| 丰满人妻熟妇乱又伦精品不卡| 男女下面进入的视频免费午夜 | 真人做人爱边吃奶动态| 国产精品一区二区三区四区久久 | 国产99久久九九免费精品| 亚洲男人天堂网一区| 久久午夜综合久久蜜桃| 亚洲精品久久成人aⅴ小说| 少妇的丰满在线观看| 女人高潮潮喷娇喘18禁视频| 真人一进一出gif抽搐免费| 黑人巨大精品欧美一区二区mp4| av视频在线观看入口| 色综合婷婷激情| 亚洲激情在线av| 国产单亲对白刺激| 亚洲自拍偷在线| 亚洲一区高清亚洲精品| 午夜福利在线在线| 99在线视频只有这里精品首页| 国产视频一区二区在线看| 亚洲av熟女| 成人永久免费在线观看视频| 成在线人永久免费视频| 亚洲av成人一区二区三| 免费高清在线观看日韩| 91在线观看av| 国产精品精品国产色婷婷| 91九色精品人成在线观看| 欧美性猛交黑人性爽| 夜夜躁狠狠躁天天躁| 欧美黑人欧美精品刺激| 一本精品99久久精品77| 国产私拍福利视频在线观看| 级片在线观看| 国产精品久久久久久精品电影 | 国内少妇人妻偷人精品xxx网站 | 欧美日韩黄片免| 青草久久国产| 亚洲精品中文字幕在线视频| 给我免费播放毛片高清在线观看| 精品不卡国产一区二区三区| 色精品久久人妻99蜜桃| av免费在线观看网站| 欧美成人一区二区免费高清观看 | 亚洲av熟女| 亚洲色图av天堂| 1024视频免费在线观看| 一边摸一边做爽爽视频免费| 禁无遮挡网站| 久久婷婷成人综合色麻豆| 麻豆国产av国片精品| 每晚都被弄得嗷嗷叫到高潮| 欧美色视频一区免费| 999精品在线视频| 亚洲色图 男人天堂 中文字幕| 亚洲精品中文字幕在线视频| 最近在线观看免费完整版| 一区二区三区高清视频在线| 欧美xxxx性猛交bbbb| 午夜福利在线在线| 欧美3d第一页| 日韩中字成人| 欧美另类亚洲清纯唯美| 久久久成人免费电影| 乱码一卡2卡4卡精品| 亚州av有码| 联通29元200g的流量卡| 女生性感内裤真人,穿戴方法视频| 成人精品一区二区免费| 淫秽高清视频在线观看| 亚洲,欧美,日韩| 99久久无色码亚洲精品果冻| 一个人看的www免费观看视频|