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

    The adsorption of acidic gaseous pollutants on metal and nonmetallic surface studied by first-principles calculation: A review

    2020-01-14 07:54:38XueYeShengguiXiJingZhishnYngWenjuJingHulinWng
    Chinese Chemical Letters 2019年12期

    Xue Ye,Shenggui M,,Xi Jing,,*,Zhishn Yng,,Wenju Jing,,Hulin Wng

    a College of Architecture and Environment, Sichuan University, Chengdu 610065, China

    b National Engineering Research Center for Flue Gas Desulfurization, Sichuan University, Chengdu 610065, China

    Keywords:

    Acidic gases

    Adsorption

    Metals

    Carbon-based materials

    DFT

    ABSTRACT

    The acidic gases such as SO2,NOx,H2S and CO2 are typical harmful pollutants and greenhouse gases in the atmosphere,which are also the main sources of PM2.5.The most widely used method of treating these gas molecules is to capture them with different adsorption materials, i.e., metal and nonmetallic materials such as MnO2, MoS2 and carbon-based materials.And doping transition metal atoms in adsorption materials are beneficial to the gas adsorption process.The first-principles calculation is a powerful tool for studying the adsorption properties of contaminant molecules on different materials at the molecular and atomic levels to understand surface adsorption reactions, adsorption reactivity, and structureactivity relationships which can provide theoretical guidance for laboratory researches and industrial applications.This review introduces the adsorption models and surface properties of these gas molecules on metal and nonmetallic surfaces by first-principles calculation in recent years.The purpose of this review is to provide the theoretical guidance for experimental research and industrial application,and to inspire scientists to benefit from first-principles calculation for applying similar methods in future work.

    1.Introduction

    In recent years, with the rapid development of the global economy,the air haze pollution problems have become increasingly prominent in China, which principal substance is PM2.5[1].PM2.5refers to atmospheric fine particles with an aerodynamics diameter less than 2.5 μm.According to the Global Burden of Disease,environmentalpollutionfromPM2.5accountedfor 3.1milliondeaths worldwide as of 2010[2].Moreover,the transportation distance of haze isverylong,whichwill reducethe visibilityof a large area of the city and greatly reduce the effective visual distance[3-5].The main components of PM2.5are secondary particles such as sulfate,nitrate and ammonium salt,and their contribution is as high as 50%,which are mainly produced by the chemical transformation of acidic gases such as sulfur dioxide (SO2), nitrogen oxide (NOx) and hydrogen sulfide (H2S) in the atmosphere [6,7].In addition, global warming caused by massive emissions of greenhouse gases such as carbon dioxide (CO2) is another major environmental problem.It causes temperature rise, sea level rise, land desertification, increased disease and insect pests and so on[8].

    These acidic gases such as SO2, NOx, H2S and CO2, are mainly emitted from anthropogenic sources.The total emissions of SO2in China accounts for about 26%of the global emissions,which is over 19 million tons per year for more than ten years, and about 85%comes from the combustion of coal-based fossil fuels[9].The NOxare mainly derived from fossil fuel combustion, metal smelting,municipal solid waste incineration and so on,and NOxemissions of China has increased to 25 million tons in 2012 [10].The main sources of H2S are natural gas processing, petroleum refinement,gas manufacturing and wastewater treatment.The CO2derived mainly from fossil fuel combustion and automobile exhaust emissions, and it is estimated that the total CO2emission is responsible for 2/3 of the growth of all greenhouse gas emissions[11].Thus, it is significantly important to take effective measures for the control of the emission of these gases into atmosphere.

    At present,the commonly used purification methods of gaseous pollutants are absorption, adsorption, catalytic and biological methods.Among them, adsorption method is one of the most widely used methods, because of the high selectivity, low cost,good separation effect and high purification efficiency [12].Adsorption is divided into physical adsorption and chemical adsorption, based on the intermolecular force and chemical bond of adsorbents and adsorbates, respectively.The adsorbents commonly used can be divided into metals and their compounds such as MnO2and nonmetallic materials such as carbon-based materials.Supported metal compounds can be used as the adsorbents for the removal of SO2, NOx, H2S and CO2, i.e., Fe, Co,Ni,Cu,Cr,Ce[13-17].Jin et al.[16]used pyrolusite slurry to absorb SO2from industrial waste gases and found that the removal rate of SO2can reach more than 90%.Ren et al.[17]found that the removal rate of SO2was nearly 100% and the leaching rate of manganese was close to 100% using a self-developed jet bubbling reactor.On the other hand, there are mainly two kinds of nonmetallic acidic gases adsorbents, including activated carbon [18-20]and zeolite molecular sieve [21-26].Mochida et al.[24]believed that SO2,O2and H2O molecules are adsorbed firstly on the carbon surface,while the formation of H2SO4can be achieved,when there is a close enough distance and the appropriate spatial configuration between them.Zeolite molecular sieve can also adsorb SO2effectively because of its good porosity and permeability [25,26].However, from the experimental operation, the various reactions can only be judged from the macroscopic change of pollutants substances, while the specific mechanism of adsorption reactions cannot be explained very well.

    First-principles calculation is valuable in understanding surface adsorption reactions, adsorption reactivity, and structure-activity relationships [27-31].The first-principles calculation is an algorithm for solving the Schr?dinger equation directly after some approximate treatment according to the principle of the interaction between nucleus and electron,and its basic motion law,using the principle of quantum mechanics, starting from the specific requirements [32].In recent years, with the vast development of various algorithms and high-performance computing equipment,first-principles calculation have been widely used in various surface chemical reactions,including surface adsorption reactions[33,34],which provides reliable theoretical guidance for practical applications and makes up for some deficiencies in experimental operation.Especially in the field of adsorption, it is possible to understand the energy change, charge density and magnetic properties of pollutant molecules adsorbed on different materials by the method of first-principles calculation [35,36], which can help us to explain the experimental results more effectively in terms of mechanism.

    Thus, this review first introduced the theoretical basis behind first-principles calculation, then some commonly software, and application examples of acidic gases (i.e., SO2, NOx, H2S and CO2)adsorbed on metal and nonmetallic surfaces studied by firstprinciples calculation are summarized.This review will provide theoretical basis for further modification of catalysts for removing these polluting gases, so as to further improve their removal efficiency,reduce energy consumption and prolong the lifetime of catalysts.It also helps to encourage researchers in this field to consider first-principles calculations when studying surface reactions at the atomic scale.

    2.Theory section

    2.1.Overview

    Quantum mechanics is a physical theory that describes microscopic matter.The quantum mechanics of matrix expression are established by W.K.Heisenberg in 1925 [37].The following year,the classical wave equation is proposed by E.Schr?dinger and it is discussed that the equation is essentially unified with the expression of quantum mechanics matrix [38].The generation of quantum mechanics provides a way to describe the motion of electrons and nuclei under arbitrary conditions [39,40].All properties of the material can be obtained by solving the Schr?dinger equation in principle.In the microscopic field, the movement of electrons determines most of the properties of the material, but the interaction between electrons in the material is highly complex.Moreover,in multi-particle research systems such as gas phase molecules, clusters, and crystals, the interaction between particles is extremely complicated.Therefore, it is difficult to understand the properties of multibody materials strictly based on the solution of Schr?dinger equation.The use of density functional theory(DFT)can solve this problem well and it is more effective to investigate the properties of the materials.

    2.2.Density functional theory

    In the 1960s, the DFT was proposed by P.Hohenberg, W.Kohn and Sham [41,42].DFT is a quantum theory based on Thomas-Fermi's theory, in which the electron density distribution n(r) is used as the basic variable to study the electronic structure of atoms and molecules.In the traditional quantum theory system, wave function is regarded as the basic physical quantity.However, the electron density is regarded as the basic variable in DFTcalculation.Hence, the solution of Schr?dinger equation is transformed into the solution of the optimal electron density based on the DFT.Meanwhile, the multi-electron problem is successfully transformed into the single-electron equation [43,44].In recent years,the combination of DFT and molecular dynamics has made a lot of progress in material calculation, and has gradually become the core technology in materials science and related fields.

    Hohenberg and Kohn put forward the basic theoretical basis of DFT, which is called Hohenberg-Kohn theorem [45,46].Although the Hohenberg-Kohn theorem establishes the framework of DFT,it has not explained how to establish its concrete form.Therefore,it cannot be directly used to solve the problem.In 1965,Kohn worked with Sham to introduce a multi-electron system with no interaction in the potential field.The charge density of the system is the same as that of the original system,and the kinetic energy of the non-interacting system can be simply written as the sum of each electron kinetic energy.With these basic theoretical frameworks and dealing with the exchange correlation energy,the Schr?dinger equation of the ground state of the complex multielectron system can be transformed into a simple single-electron equation, which brings great convenience to scientific research.However,they do not accurately express the exchange correlation energy,so it is still difficult to get practical application.In order to solve this problem, the local density approximation (LDA) was proposed after the Kohn-Sham equation was put forward in the same year[46].LDA assumes that the functional of charge density of the system changes slowly in space,so the charge density in each volume element is affected only by its spatial position.This approximate method can accurately find the exchange correlation functional and is suitable for calculating the uniform distribution of charge density in space.However,the LDA approximation has a great error in the calculation of transition metals in systems where the charge density varies more sharply.For the sake of obtain more accurate results in different calculation processes, the inhomogeneity of charge distribution must also be considered.To solve this problem,Becke introduced the gradient of charge density in 1988 to modify the result of LDA calculation, which resulted in the generalized gradient approximation (GGA) [47].The method assumes that the exchange correlation energy is related to the local charge density and the nearby charge density and can deal with the non-uniform system, so the calculated results are more accurate.At present, the GGA method has developed a variety of models, such as the common Perdew-Burke-Emzetho (PBE) [48],Perdew Wang 91 (PW91) [49], Becke and BLYP exchange correlation potential.At present,in order to improve the accuracy of theoretical calculation, the hybrid functional method is introduced, that is, the exchange correlation functional of the system is obtained by using the exchange energy of Hartree-Fock and the exchange energy in DFT as a linear combination.Such as the commonly used HSE06 hybrid functionals and B3LYP functionals in Gaussian calculation method [50].

    2.3.Software packages based on DFT

    With the above theoretical foundation, DFT has been continuously developed and widely used in practice.In order to satisfy the huge computational requirements of first-principles calculation,scientists have developed various computing programs and software after combining with the modern computer technology through continuous research and experiments.This software can be used on a variety of large computers, servers and personal computers with a wider range of applications and better compatibility.And these programs and software have different features and can be applied to different computing systems.Therefore, researchers can choose the software according to their actual demand in practical application.Here, we will introduce several kinds of software commonly used by researchers: VASP,Material Studio and Gaussian.

    The VASP(Vienna Ab-initio Simulation Package)is a commercial software package to describe the electronic-core interaction used in quantum mechanical and molecular dynamics calculation,which comes from the Hafner team of the University of Vienna[51,52].It is one of the most popular software in the related fields of science[53-55].VASP used plane wave pseudopotential method to carry out ab initio simulation, and used periodic boundary conditions to deal with atoms, molecules, clusters, crystals, thin films, solids and surface systems, etc.The electronic structure,optical and magnetic properties of the crystal can be calculated more exactly.Similarly, some basic physical and chemical properties of the material surface can also be calculated accurately.VASP has been improved step by step as the continuous development of computer science and technology.Now VASP has been widely used in the optimization of complex systems and the adsorption of molecules at the surface,and can be used for high efficiency parallel computation of mainframe computers.

    Material Studio is a powerful platform for material science researchers,and it is also a powerful modeling tool for configuration optimization, property prediction, X-ray diffraction analysis and the calculation of complex dynamic simulation.Material Studio contains several modules to perform varies calculations,including Materials Visualizer, Discover, COMPASS, Amorphous Cell, Reflex, DMol3, CASTEP and other modules.Among them,DMol3and CASTEP are widely applied in various areas.Unlike other quantum mechanical programs,DMol3uses numerical functions to describe atomic orbitals.The method describing the electronic state of the system is the linear combination of atomic orbitals,which results in extremely high computational accuracy and efficiency.DMol3can also simulate the chemical reaction processes of gaseous, liquid, and solid surfaces, which is one of its major characteristics.Therefore,it has been widely applied in the fields of physics,chemistry, biology, chemical engineering, in particular to the field of catalytic reaction and other molecular reactions.CASTEP uses the plane wave pseudopotentials to carry out the firstprinciples calculation.It is often used to study metals, semiconductors,ceramics,and so on.The surface chemical and optical properties of crystals and the wave functions of the system can be obtained.Now it has been widely used in solid state physics,material science and other fields.

    Gaussian is also a comprehensive software package widely used in the field of quantum chemistry.It can be applied not only to density functional method,but also to molecular dynamics,semiempirical method and so on.Gaussian is often used to calculate molecular orbitals,vibrational frequencies,thermodynamic properties,FT-IR and Raman spectra,multiple moment,reaction paths,etc.The calculation can be performed on the ground state or excited state of the system and predict the energy, structure and molecular orbital of the periodic system.It is suitable for the study of chemical reaction systems [56,57].

    3.Applications

    The adsorption properties of several acidic gaseous pollutants(SO2,NOx,H2S and CO2)on different materials are discussed in this section.Adsorption materials we discussed mainly include metals and their compounds, and non-metallic materials, such as graphene and nanotubes.DFT method is used to calculate the adsorption process and properties of these substances at different interfaces.Some representative examples are selected to discuss the application of first-principles calculation to the adsorption of these pollutants.

    3.1.Adsorption of acidic gas molecules on metal surface

    Transition metals and precious metals play a significant role in the catalytic oxidation of some gas pollutant molecules such as SO2and NOx, and often exist as catalysts with excellent performance[58-66].In general,the experimental study is to make metals into large single crystals.After cutting the single crystals,the required crystal faces can be obtained, which provides a material basis for the adsorption of the latter molecules.The lattice of the metal is replaced by different forms of crystal planes, which are represented by Miller indices.For example,the metal Fe has(100),(110)and(111) plane.At the same time, certain metal oxides may have some properties which are not found in single metal, such as copper oxides have higher sulfur capacity and higher affinity to SO2molecules as desulfurizer [67].From the last century, some scholars have studied the adsorption behavior of gas molecules on some metal surfaces by DFT theories[68-72].Therefore,the study of the adsorption or decomposition mechanism of these pollutants on the surface of metal and its compounds is of great significance for improving the existing catalysts and producing new catalysts.

    The metal monolayer-stanene material is a newly discovered zero-gap semiconductors with buckled honeycomb structure,which has attracted extensive attention in material science and other disciplines [73,74].Chen Xianping et al.[75]studied the adsorption properties of H2S,NO,NO2and other gas molecules on the surface of stanene by first-principles calculations.There are three absorbable positions in stanene, which are the top of Sn atom,the center of a Sn hexagonal and the center of Sn-Sn bond.It is found that the interaction of H2S and stanene is weak and belongs to physical adsorption.There is a strong interaction between NO and stanene,and the adsorption energy is-0.801 eV.However, the adsorption energy of NO2on stannene is above-1.12 eV,and its configuration is shown in Fig.1.At the same time,it is found from density of states (DOS) that the Sn p and O p orbitals of the NO2/stanene structure have strong orbital hybridization, which shows that the adsorption structures for NO2adsorbed on stanine are quite stable.And the charge density difference (CDD) also proved this result.The results indicate that stanene has a potential application in the fields of gas sensor and high-performance catalyst for the detection and adsorption of NO and NO2.

    Fig.1.Top and side view of the lowest energy structure of NO2 adsorbed on stanene.Reproduced with permission [75].Copyright 2016, The Journal of Physical Chemistry C.

    In recent years,some oxygen-functionalized MXenes materials,i.e., M2CO2(M=Sc, Hf, Zr, Ti), have been found to have a broad prospect in gas molecular sensing or capture owing to their semiconductor characters [76-78].Ma Shuhong et al.[79]have calculated certain properties of SO2gas molecules adsorption on M2CO2(M=Sc,Hf,Zr,Ti).By finding the favorable configuration of M2CO2and the most stable structure of SO2adsorbed on the surface of the material,it is found that Sc2CO2is the most preferred for the adsorption of SO2molecules.Its adsorption energy is-0.646 eV.As seen from Fig.2, the significant orbital hybridizations between SO2and Sc2CO2are localized near the Fermi level,so that SO2adsorbs on Sc2CO2.Therefore, the charge transfer between the SO2and Sc2CO2increases multiply, which is consistent with the CDD diagram on the right of Fig.2, indicating that the adsorption is relative stable.Moreover, the method of improving the adsorption capacity of SO2by applying tensile strains or negative electric fields has been proposed.All these results show that Sc2CO2can be considered as a gas sensor or adsorption material for SO2.

    MoS2is a new type of transition metal 2D material with excellent structure and properties, such as high activity, high selectivity and high stability, now it has attracted extensive attention in the fields of chemistry [80-82].On the basis of the first-principle calculations, Wei Huangli et al.[83]studied the adsorption properties of the original MoS2monolayer and the molecular layer doped with Ni metal atoms to the H2S and SO2gas molecules,and compared the adsorbed structural parameters,the adsorption energy and the amount of charge transfer.The adsorption energy of H2S and SO2reaches -1.319 eV and-1.382 eV, respectively, indicating that the interaction between the two gases and Ni-doped MoS2monolayers is pretty strong.The doping of metal atoms provides a substantial number of free electrons for MoS2, and the strong orbital hybridization with adsorbed molecules improves the adsorption capacity of the MoS2molecular layer.Accordingly,Ni-doped MoS2might be a promising gas absorbent material to remove H2S and SO2gas pollutants.

    Fig.2.PDOS of SO2 molecule (green shadow), Sc2CO2 monolayers with (red shadow)and without(blue solid line)SO2 adsorption,respectively.The side views of charge density difference are given in the upright panels.Reproduced with permission [79].Copyright 2017, The Journal of Physical Chemistry C.

    Table 1 summarizes the adsorption examples of acidic gas pollutant molecules on different metals and their compounds,and gives the software, methods and results of adsorption[67,75,79,83-92].According to the reported results,gas molecules may be adsorbed on the surface of materials in different states,and there are many adsorption sites on the surface of materials.Therefore,the adsorption configurations of these gas molecules on different surfaces of materials are not one way but many ways.In order to obtain the desired results briefly, the most stable configuration of adsorption and the corresponding results are selected.It can be seen from Table 1 that the adsorption energy of the gas molecules is relatively large,which belongs to the chemical adsorption process.Moreover, we found that the adsorbent materials are mainly transition metal elements.In addition,doping another transition metal element into metal compounds can significantly improve the adsorption properties of the materials [83,86].These results indicate that these acidic gas pollutant molecules could undergo a great adsorption reaction on the transition metal and its compounds.

    3.2.Adsorption of acidic gas molecules on non-metallic surface

    3.2.1.Adsorption of gas molecules on graphene surface

    Graphene(GP)is a kind of two-dimensional monatomic carbon simple material with perfect hexagonal structure [93].It is also a multifunctional catalyst carrier with high conductivity, high structure stability and high mechanical strength.Its structure is a two-dimensional crystal structure composed of hexagonal crystals.There are two kinds of carbon atoms in each unit cell.After oxidation with strong oxidant, various graphene oxides can be obtained,and it is easy to be modified.At the same time,it can form a stable structure with other atoms to carry out a series of reaction processes.For example, the chemical activity of the system can be enhanced by decoration of graphene with adatoms or dopants,and the bonding strength of gases on graphene can be improved at the same time[94-98].Therefore,graphene has been widely used in chemical catalysis, adsorption and other related fields.

    Based on DFT,Diego Cortés-Arriagada et al.[99]calculated the adsorption behavior of CO2,SO2and H2S and other pollutant gases on graphene nanosheets doped with Fe atoms.The adsorption capacity of the original graphene and the doped graphene for gas molecules was compared.The adsorption of H2S onto pristine graphene and Fe-doped graphene is shown in Fig.3.The results show that the adsorption energy of CO2, SO2and H2S gas on graphene are about-0.11,-0.28 and-0.15 eV,respectively,while the adsorption energy on Fe-GP is about-0.54,-1.80 and-1.19 eV.

    The chemical adsorption is still strong even at room temperature.It is demonstrated that graphene doped with Fe atom is a feasible catalyst or absorbent.Chen Chi et al.[100]reported the adsorption of CO2,NO2,SO2on metal(Li,Al,Ti)modified graphene oxide (GO) by VASP by using the first-principles calculations method.Hydroxyl and epoxy groups were placed on the hexatomic ring to obtain the most stable configurations of metal decorated GO.It is found that the binding energies of CO2,NO2and SO2on Li,Al, Ti decorated GO surface are about -0.19 to -0.91 eV, -1.25 to-3.60 eV and -0.44 to -2.15 eV, respectively.However, other scholars have studied that the binding energies of CO2, NO2and SO2on graphene surface are-0.40[101],0.07[102]and 0.01[103]eV,respectively.In addition,the adsorption distance between gas molecules and GO surface and the results of Bader charge analysis are consistent with the results of binding energies, as shown in Fig.4,an example is given.The results show that metal decorated GO has a significant effect on the capture of CO2,NO2and SO2gases.

    Ali Shokuhi Rad[104]also used the first-principles calculation method to study the adsorption properties of NO2on Al-doped graphene.In the pristine graphene sheets, one C atom was substituted by one Al atom to form a stable structure.The results show that the adsorption energy of NO2molecules on origin graphene is about -0.002 eV, while the NO2have the lager adsorption energy on the Al-doped graphene of -0.645 eV,demonstrating that the presence of Al increases the capability ofgraphene to adsorb these NO2molecules.The results of charge transfer and other properties also illustrate that doping metal atoms can significantly improve the interaction between adsorbates and substrate.Therefore, Al-doped graphene material is expected to be an excellent gas sensor material or adsorption material for NO2.

    Table 1 Summary of adsorption of acidic gas molecules on metal surface.

    Fig.3.The adsorption of H2S onto pristine graphene (G) and Fe-doped graphene(FeG).Reproduced with permission [99].Copyright 2018,Applied Surface Science.

    The adsorption properties of SO2on the surface of graphene doped with three groups was calculated by Zhang Huijuan et al.[105]with van der Waals force modified DFT method.The first was doped with hydroxyl group, the second was doped with epoxy group, and the third was doped with both hydroxyl and epoxy group.Three kinds of graphene oxides (GO) were obtained.The adsorption structure,adsorption energy and charge transfer were calculated and compared with the adsorption of SO2on graphene.The results reveal that the metallicity of GO decreases with the addition of oxygen functional groups,which leads to the decrease of SO2adsorption by GO.However,the adsorption energy of SO2on graphene surface doped with hydroxyl group is larger than pristine graphene because of the formed hydrogen bond between hydroxyl group and SO2.Through the Bader charge analysis and differential charge density analysis, it can be found that the existence of hydroxyl group can promote the charge transfer, while the existence of epoxy group will hinder the charge transfer.It is concluded that GO doped with hydroxyl groups can enhance the adsorption of SO2.

    Table 2 summarizes the adsorption examples of these acidic gas pollutant molecules on different graphene,and gives the software,methods and results of adsorption [99,100,104-111].The most stable configuration of adsorption and its corresponding calculation results are still selected.We can see from the results in Table 2 that the pure graphene is different from the structurally modified graphene.Although the gas adsorbed is different, they all have a relatively similar rule: (i) The adsorption energy between pure graphene and gas molecules is small and belongs to weak physical adsorption, while the adsorption distance between the molecule and the surface of graphene is longer and the charge transfer amount is less.Therefore, their adsorption is extremely unstable,while the structure modified graphene often adsorbs these gas molecules by chemisorption.(ii) The adsorption energy of modified graphene is larger than that of the pristine graphene and the adsorption distance is closer.The charge transfer is also changed to a certain extent.Therefore,they are easy to form stable adsorption structure.According to the results, we found that the adsorption performance of graphene doped with transition metal atoms on gas molecules was greatly improved in the modified graphene [106,108,111].This indicates that doped groups or transition metal atoms in graphene is favorable to gas adsorption compared with the original structure of graphene,which is vastly different from the adsorption properties on the metal surface.

    Fig.4.The bond angles(green)of acidic gases(a is CO2;b is NO2;c is SO2),the distances(black)from gases to metals and the electron transfers(blue)between gases and metal decorated GO (a-c is Li decorations).Reproduced with permission [100].Copyright 2014, Physical Chemistry Chemical Physics.

    Table 2 Summary of the adsorption of acidic gas molecules on graphene surface.

    3.2.2.Adsorption of gas molecules on nanotubes surface

    Since the discovery of carbon nanotubes (CNTs) in 1991 [112],scientists have found that CNTs have many special properties,such as great strength and elasticity, nanoscale space.CNTs are essentially graphene sheets rolled into seamless tubular monolayers, and composed of carbon atoms on the surface [113].According to the different surface shapes of carbon atoms,they can be divided into armchair type,sawtooth type and spiral type.They have good bonding properties and have greater reaction activity than other graphite.The pore structure on the surface and inside of carbon nanotubes determines that it can be used as a good gas adsorbent and has a strong adsorption potential energy for some gases.The properties of BN are similar to that of C, and carbon nanotubes and fullerenes may be formed when there are layered structures in bulk structures in theory.Therefore,the BN layer was curled up like C layer to form BN nanotubes (BNNTs), and the structural properties of BNNTs were studied.Consequently, using carbon nanotubes and BNNTs as a new gas adsorbent candidate and studying their adsorption mechanism theoretically are of great significance for breaking through the traditional adsorption theory and methods.

    Based on DFT calculations, the sensing and adsorption properties of NO2gas molecules on single-wall carbon nanotubes(SWCNTs) doped with metal atoms were calculated by Icell M.Sharafeldin et al.[114].The(8,0)SWCNT was selected as the matrix material,and a single carbon atom was replaced with a Cu,Pt or Ti to form a stable nanotube structure.The adsorption energy of NO2on their surface come to reach -1.827 eV, -2.262 eV and-2.653 eV, respectively.However, it has been reported that the adsorption energy of NO2on the surface of pristine CNTs is only about -0.267 eV [115], which indicates that there are strong interactions between metal doped SWCNTs and NO2gas molecules.These conclusions prove that SWCNTs doped with Cu,Pt or Ti materials has high sensitivity to NO2molecules.Jiao Yan et al.[116].calculated the adsorption behavior of(7,7)SWCNTs for CO2under several nitrogen doping conditions, including graphite nitrogen,pyridine nitrogen and tripyridine nitrogen.The results show that the adsorption energy of CO2molecules by all structures is less than 0.3 eV.However, when four extra electrons are injected into CNT doped with pyridine nitrogen,CO2molecules can be adsorbed spontaneously by forming C-N bonds,and the adsorption capacity is significantly improved.Fig.5 shows the adsorption configurations with 0,2,and 4 electrons injection.According to this,the feasibility of electrocatalytic in gas adsorption and separation is considered, and a new thought for the preparation of gas adsorbents is provided.

    Fig.5.Adsorption configurations of CO2 on pyridine-CNT with 0,2,and 4 electrons injection.Reproduced with permission [116].Copyright 2014, ChemSusChem.

    Fig.6.SO2 adsorbed on MnN-BNNT site.(a) Difference spin charge density of side view.(b)PDOS.Reproduced with permission[117].Copyright 2015,Applied Surface Science.

    The adsorption behavior of SO2molecules on BNNTs surface of original and doped with Mn atom was studied by Deng Zunyi et al.[117].The simulated model uses zigzag (8, 0) BNNT to form two different Mn-BNNTconfigurations by replace one B atom or N atom with a Mn atom.The formation energy of SO2molecule by pristine BNNT is only about-0.056 eV,which indicates that SO2undergoes weakly physical adsorption.However, the formation energy of B and N sites of SO2on Mn-BNNT is about-0.993 eV and-1.569 eV,which is much stronger than the original ones.As can be seen from Fig.6, there is a strong covalent interaction and charge transfer between SO2and MnN-BNNT.Therefore, the doping of Mn improves the adsorption capacity of BNNT remarkably, and provides a theoretical basis for the further application of BNNTs.

    Table 3 summarizes the adsorption examples of common acidic gaseous pollutant molecules on the surface of nanotubes,and gives the software,methods and results of adsorption[36,114-121].The most stable configuration of adsorption and its corresponding calculation results are still selected.According to the results of literature, we find that the adsorption ability of gas molecules on the surface of carbon nanotubes is weak and belongs to physical adsorption.However, on the surface of carbon nanotubes doped with metal atoms, the adsorption ability of gas molecules is relative stronger,and the adsorption structure is more stable.This is similar to the adsorption results of graphene.Similarly, we can conclude that the doping of metal atoms facilitates the adsorption of gas molecules on carbon nanotubes.It also shows that carbon nanotubes are theoretically feasible as adsorbents for acidic gaseous pollutant molecules.

    4.Conclusions and prospects

    The first-principles calculation, especially the DFT method, is extremely useful for studying the adsorption process of gases on different materials.In this paper, we summarize the adsorption properties of several common acidic pollutants in the atmosphere on the surface of several kinds of materials.It was discussed on the calculated results of adsorption energy,the change of bond length and the configuration before and after adsorption,and the charge transfer of adsorption.The results suggested that SO2,NOx,H2S and CO2can be adsorbed on metal and its compounds, graphene and nanotubes.Among them, metals and its compounds can achieve good adsorption effect without doping.The adsorption effect of the original graphene and carbon nanotubes is weak, while the adsorption performance can be greatly improved by doping the functional groups or metal atoms, especially the transition metal atoms.Accordingly, it is concluded that doping transition metal atoms in these adsorption materials is beneficial to the gas adsorption process.In addition,the first principles can investigate more properties of molecular adsorption, such as bond angle,energy band,density of state,magnetic properties,which have not been summarized in this paper.Through the calculation of these properties,we can obtain a series of physical and chemical changebetween them,which can explain some of the results obtained in the experimental research and provide theoretical support for the actual research results.Therefore,this is a significant and valuable work.

    Table 3 Summary of the adsorption of acidic gas molecules on nanotubes surface.

    Although the first-principles calculation has been widely used in adsorption in recent years,we still put forward some prospects and suggestions: i) Metal and its compounds and graphene have been used to adsorb SO2,NOxand other gas molecules,which have been widely used in experimental research and some applications in practical engineering.However,the application of nanotubes in gas adsorption is still very few, and most of them are used in gas sensors.Hence,we hope that this review will give some inspiration to the researchers in the study of new adsorption materials.ii)According to the summary of this paper,it can be found that noble metal elements have very good effect on the adsorption of several gas molecules,but the cost of using these precious metal materials in the actual experiment is relatively high.Therefore,we hope that this review can be used as a theoretical reference for researchers to study the adsorption of related gas molecules on related materials,so as to develop more excellent adsorbents with good performance, and even develop more economical and green adsorption materials that can replace precious metal elements.iii) The firstprinciples calculation often has ideal conditions for simulating the molecular and material surfaces, and the adsorption materials used in the actual experiments are not pure substances.In addition, the experimental conditions are changeable and cannot be simulated,so in some cases,the first-principles calculation can only represent a trend.In consequence, more optimization is needed in the future in order to achieve results more consistent with actual experiments.

    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.

    Acknowledgments

    This work is supported by the National Nature Science Foundation of China (No.51778383) and Sichuan Science and Technology Program (No.2019YFS0500).

    中文乱码字字幕精品一区二区三区 | 亚洲av电影在线观看一区二区三区 | 亚洲欧美日韩卡通动漫| 国产老妇伦熟女老妇高清| 建设人人有责人人尽责人人享有的 | 国产精品日韩av在线免费观看| 啦啦啦中文免费视频观看日本| 99久国产av精品国产电影| 国产精品蜜桃在线观看| 十八禁网站网址无遮挡 | 国产亚洲一区二区精品| 亚洲精华国产精华液的使用体验| 久久久久久九九精品二区国产| 高清日韩中文字幕在线| 男插女下体视频免费在线播放| 非洲黑人性xxxx精品又粗又长| 蜜桃久久精品国产亚洲av| 成年女人在线观看亚洲视频 | 久久久久久久久久久免费av| 熟女电影av网| 亚洲乱码一区二区免费版| 亚洲18禁久久av| 久久久精品欧美日韩精品| 国产精品久久视频播放| 蜜桃亚洲精品一区二区三区| 天堂av国产一区二区熟女人妻| 久久久久国产网址| 亚洲欧美精品专区久久| 亚洲内射少妇av| 国产免费视频播放在线视频 | 国产av在哪里看| 精品一区二区免费观看| 嫩草影院精品99| 欧美成人a在线观看| 99热全是精品| 秋霞伦理黄片| 国产有黄有色有爽视频| 久久久久精品久久久久真实原创| 观看美女的网站| 午夜亚洲福利在线播放| 精品一区在线观看国产| 国产精品日韩av在线免费观看| 日韩中字成人| 亚洲伊人久久精品综合| 综合色丁香网| 中文乱码字字幕精品一区二区三区 | 久久精品国产亚洲av涩爱| 国产免费福利视频在线观看| 精品一区二区免费观看| 91久久精品国产一区二区成人| 国产亚洲av嫩草精品影院| 国产亚洲一区二区精品| 精品国内亚洲2022精品成人| 国产精品一区二区在线观看99 | 亚洲精品国产成人久久av| 成人欧美大片| .国产精品久久| 日韩精品青青久久久久久| 99热网站在线观看| 97热精品久久久久久| 熟妇人妻不卡中文字幕| 日韩国内少妇激情av| 亚洲av中文av极速乱| 2018国产大陆天天弄谢| 在现免费观看毛片| 大又大粗又爽又黄少妇毛片口| av福利片在线观看| 精品人妻视频免费看| 日韩强制内射视频| 亚洲人与动物交配视频| 真实男女啪啪啪动态图| 两个人的视频大全免费| 国产高清国产精品国产三级 | 久久精品久久精品一区二区三区| 午夜福利成人在线免费观看| 最近中文字幕2019免费版| 日韩中字成人| 久久99热这里只频精品6学生| 午夜激情福利司机影院| kizo精华| 26uuu在线亚洲综合色| 美女xxoo啪啪120秒动态图| 最近的中文字幕免费完整| 美女被艹到高潮喷水动态| 亚洲av日韩在线播放| 免费不卡的大黄色大毛片视频在线观看 | 80岁老熟妇乱子伦牲交| 亚洲精品aⅴ在线观看| 亚州av有码| 一级毛片我不卡| av线在线观看网站| av福利片在线观看| 国产精品一及| 男人爽女人下面视频在线观看| 韩国av在线不卡| 亚洲av不卡在线观看| 精品人妻一区二区三区麻豆| 亚洲av成人精品一区久久| 精品国内亚洲2022精品成人| 欧美不卡视频在线免费观看| 一级毛片黄色毛片免费观看视频| 国产熟女欧美一区二区| 国产精品综合久久久久久久免费| 少妇人妻一区二区三区视频| av又黄又爽大尺度在线免费看| 国产亚洲91精品色在线| 国产亚洲最大av| 精品久久久久久久末码| 亚洲欧美中文字幕日韩二区| 久久精品国产鲁丝片午夜精品| 亚洲人成网站高清观看| 国产精品不卡视频一区二区| 熟女电影av网| 午夜日本视频在线| 国产成人91sexporn| 国产高清国产精品国产三级 | 男女下面进入的视频免费午夜| 狂野欧美激情性xxxx在线观看| av卡一久久| 国产精品综合久久久久久久免费| 久久精品久久久久久噜噜老黄| 午夜福利视频精品| 成人欧美大片| 亚洲熟女精品中文字幕| 春色校园在线视频观看| 小蜜桃在线观看免费完整版高清| eeuss影院久久| 91精品伊人久久大香线蕉| 中文字幕人妻熟人妻熟丝袜美| 免费av不卡在线播放| 亚洲av免费在线观看| 99久国产av精品国产电影| 天堂av国产一区二区熟女人妻| 91久久精品电影网| 国产精品久久久久久精品电影| 91在线精品国自产拍蜜月| 午夜福利在线观看吧| 黄色欧美视频在线观看| 国产成人精品久久久久久| 久久国产乱子免费精品| 男插女下体视频免费在线播放| 国产精品久久久久久久久免| 免费av观看视频| 国产熟女欧美一区二区| 一区二区三区高清视频在线| 午夜激情久久久久久久| 久久午夜福利片| 成人性生交大片免费视频hd| 国产免费又黄又爽又色| 欧美最新免费一区二区三区| 91av网一区二区| 菩萨蛮人人尽说江南好唐韦庄| 久久久久精品久久久久真实原创| 91久久精品电影网| 亚洲av二区三区四区| 精品久久久久久成人av| 嫩草影院精品99| 一边亲一边摸免费视频| 亚洲真实伦在线观看| 中文天堂在线官网| 内射极品少妇av片p| 亚洲av电影不卡..在线观看| 国产老妇伦熟女老妇高清| 最后的刺客免费高清国语| 超碰av人人做人人爽久久| 国产高清有码在线观看视频| 国产精品爽爽va在线观看网站| 亚洲欧美精品自产自拍| 边亲边吃奶的免费视频| 日韩av在线大香蕉| 欧美 日韩 精品 国产| 国产淫语在线视频| 亚洲欧美一区二区三区黑人 | 欧美成人精品欧美一级黄| 三级经典国产精品| 日韩制服骚丝袜av| 色综合站精品国产| 亚洲精品国产av蜜桃| 国产伦在线观看视频一区| 久久热精品热| 一区二区三区四区激情视频| 亚洲精品国产av蜜桃| 免费大片黄手机在线观看| 中文欧美无线码| 国产视频首页在线观看| 国产一级毛片在线| 国产伦一二天堂av在线观看| 免费观看a级毛片全部| 亚洲国产精品成人综合色| 干丝袜人妻中文字幕| 国产又色又爽无遮挡免| 亚洲欧美精品专区久久| 一边亲一边摸免费视频| 日韩 亚洲 欧美在线| 亚洲av成人精品一区久久| 亚洲色图av天堂| 汤姆久久久久久久影院中文字幕 | 亚洲av在线观看美女高潮| 日韩一区二区三区影片| 18+在线观看网站| 久久久久久久久久人人人人人人| 久久精品夜夜夜夜夜久久蜜豆| 亚洲aⅴ乱码一区二区在线播放| 国产亚洲一区二区精品| 久久久久免费精品人妻一区二区| 亚洲精品日韩在线中文字幕| 国产在视频线精品| 欧美97在线视频| 欧美日韩在线观看h| 一级二级三级毛片免费看| 国产91av在线免费观看| 白带黄色成豆腐渣| 精品一区二区免费观看| 国产视频内射| 看免费成人av毛片| eeuss影院久久| 久久久久精品性色| 亚洲av中文字字幕乱码综合| 午夜福利网站1000一区二区三区| 校园人妻丝袜中文字幕| 久热久热在线精品观看| 国内精品一区二区在线观看| 亚洲欧洲日产国产| 日韩 亚洲 欧美在线| 白带黄色成豆腐渣| 2018国产大陆天天弄谢| 精品人妻熟女av久视频| 久久久久免费精品人妻一区二区| 亚洲欧美精品专区久久| 亚洲aⅴ乱码一区二区在线播放| 成人综合一区亚洲| av在线观看视频网站免费| 丰满乱子伦码专区| av在线老鸭窝| 国产午夜福利久久久久久| 国产男人的电影天堂91| 欧美高清性xxxxhd video| 国产男女超爽视频在线观看| 国内揄拍国产精品人妻在线| 三级国产精品片| 尾随美女入室| 99久久精品热视频| 国产在视频线精品| 午夜免费男女啪啪视频观看| 我的女老师完整版在线观看| 简卡轻食公司| 你懂的网址亚洲精品在线观看| av又黄又爽大尺度在线免费看| 91在线精品国自产拍蜜月| 国产成人a区在线观看| 国产高潮美女av| 少妇的逼好多水| 一本一本综合久久| 午夜精品国产一区二区电影 | 高清av免费在线| 日本爱情动作片www.在线观看| 国产亚洲精品久久久com| 韩国av在线不卡| 国产精品蜜桃在线观看| 成人美女网站在线观看视频| 人妻制服诱惑在线中文字幕| 五月伊人婷婷丁香| 高清日韩中文字幕在线| 搡老乐熟女国产| 一级毛片久久久久久久久女| 午夜激情欧美在线| 2021少妇久久久久久久久久久| a级毛色黄片| 日韩av在线大香蕉| 男人和女人高潮做爰伦理| 在线免费十八禁| 日本一二三区视频观看| 亚洲国产精品成人久久小说| 亚洲国产成人一精品久久久| 欧美潮喷喷水| 老司机影院毛片| 亚洲欧洲国产日韩| 久久久久久国产a免费观看| 一二三四中文在线观看免费高清| 最后的刺客免费高清国语| 欧美丝袜亚洲另类| 在现免费观看毛片| 色视频www国产| 国产高清三级在线| 欧美日韩一区二区视频在线观看视频在线 | 成人午夜精彩视频在线观看| 色尼玛亚洲综合影院| 91狼人影院| 亚洲电影在线观看av| 天堂中文最新版在线下载 | 国产在线男女| 女人久久www免费人成看片| 国内精品美女久久久久久| 国产真实伦视频高清在线观看| 美女国产视频在线观看| 国产黄片视频在线免费观看| 男插女下体视频免费在线播放| 乱系列少妇在线播放| 97超碰精品成人国产| 国产毛片a区久久久久| 美女脱内裤让男人舔精品视频| 日韩av不卡免费在线播放| 亚洲精品色激情综合| 五月玫瑰六月丁香| 亚洲婷婷狠狠爱综合网| 国产麻豆成人av免费视频| 成年av动漫网址| 一个人看视频在线观看www免费| 亚洲电影在线观看av| 亚洲欧美日韩无卡精品| 亚洲av.av天堂| 噜噜噜噜噜久久久久久91| 一本久久精品| 欧美变态另类bdsm刘玥| 26uuu在线亚洲综合色| 老司机影院成人| 九九在线视频观看精品| 欧美性猛交╳xxx乱大交人| 精品国内亚洲2022精品成人| 久久99蜜桃精品久久| 精品人妻熟女av久视频| 一区二区三区高清视频在线| 啦啦啦中文免费视频观看日本| 日日摸夜夜添夜夜添av毛片| 一级毛片 在线播放| 成人亚洲精品av一区二区| 国产女主播在线喷水免费视频网站 | 亚洲成人一二三区av| 国产男女超爽视频在线观看| 我要看日韩黄色一级片| 国产伦理片在线播放av一区| 日韩三级伦理在线观看| 亚洲国产精品国产精品| 日韩精品青青久久久久久| 国产欧美另类精品又又久久亚洲欧美| 天堂俺去俺来也www色官网 | 午夜精品国产一区二区电影 | 高清av免费在线| 春色校园在线视频观看| 日本黄色片子视频| 久久鲁丝午夜福利片| 最新中文字幕久久久久| 激情五月婷婷亚洲| 中文乱码字字幕精品一区二区三区 | 午夜免费观看性视频| 亚洲欧美精品自产自拍| 亚洲av成人精品一二三区| 日本免费a在线| 国产老妇伦熟女老妇高清| 国产精品麻豆人妻色哟哟久久 | 嫩草影院精品99| 国产黄色视频一区二区在线观看| 日韩电影二区| 色网站视频免费| 亚洲欧美精品专区久久| 免费观看无遮挡的男女| 青春草国产在线视频| 国产免费又黄又爽又色| 看免费成人av毛片| 麻豆久久精品国产亚洲av| 亚洲精品成人久久久久久| 777米奇影视久久| 建设人人有责人人尽责人人享有的 | 97热精品久久久久久| 精品人妻一区二区三区麻豆| 亚洲精品日韩在线中文字幕| 淫秽高清视频在线观看| 国产中年淑女户外野战色| 免费人成在线观看视频色| 波野结衣二区三区在线| 精品一区二区三区人妻视频| 精品午夜福利在线看| 亚洲精品影视一区二区三区av| 青春草国产在线视频| 亚洲丝袜综合中文字幕| 午夜福利在线观看免费完整高清在| 亚洲国产欧美人成| 成人av在线播放网站| 寂寞人妻少妇视频99o| 69av精品久久久久久| 乱码一卡2卡4卡精品| 一级毛片 在线播放| 国产69精品久久久久777片| 一个人看的www免费观看视频| 日韩av在线大香蕉| 久久久欧美国产精品| 看非洲黑人一级黄片| 五月天丁香电影| 日韩av在线大香蕉| 午夜精品国产一区二区电影 | 日韩av在线大香蕉| av专区在线播放| 亚洲国产av新网站| 一级av片app| 欧美高清成人免费视频www| 国产白丝娇喘喷水9色精品| 免费播放大片免费观看视频在线观看| 日日干狠狠操夜夜爽| 国产精品一二三区在线看| 三级经典国产精品| 26uuu在线亚洲综合色| 一区二区三区乱码不卡18| 免费看日本二区| av黄色大香蕉| 久久久久久久久久人人人人人人| 国产一区二区在线观看日韩| 老司机影院成人| 久久精品国产亚洲网站| 亚洲色图av天堂| 婷婷色麻豆天堂久久| 久久人人爽人人爽人人片va| av线在线观看网站| 在线a可以看的网站| 国产真实伦视频高清在线观看| 国内精品宾馆在线| 国产高清三级在线| 日韩在线高清观看一区二区三区| 免费看不卡的av| 欧美日韩在线观看h| 九色成人免费人妻av| 人妻少妇偷人精品九色| 国产淫语在线视频| 国产黄色免费在线视频| 日日撸夜夜添| 国产精品国产三级国产专区5o| 免费看不卡的av| 国产黄频视频在线观看| 舔av片在线| 免费黄网站久久成人精品| 亚洲最大成人中文| 亚洲精品乱久久久久久| 中文字幕免费在线视频6| 亚洲av二区三区四区| 成人漫画全彩无遮挡| 国产久久久一区二区三区| 国产永久视频网站| 亚洲四区av| 国产精品伦人一区二区| 亚洲av国产av综合av卡| 嫩草影院新地址| 成人综合一区亚洲| 国产亚洲5aaaaa淫片| 少妇裸体淫交视频免费看高清| 可以在线观看毛片的网站| 国产一区二区三区av在线| 国产精品久久视频播放| 成人亚洲欧美一区二区av| av在线播放精品| 日日摸夜夜添夜夜添av毛片| 日日摸夜夜添夜夜爱| 亚洲欧美成人精品一区二区| 七月丁香在线播放| 天堂影院成人在线观看| 日本色播在线视频| 最近最新中文字幕大全电影3| 亚洲精品日韩av片在线观看| 中文字幕人妻熟人妻熟丝袜美| 国产精品一区二区三区四区久久| 色尼玛亚洲综合影院| 777米奇影视久久| 国产欧美日韩精品一区二区| 九色成人免费人妻av| 国产在视频线在精品| 日日撸夜夜添| 熟妇人妻不卡中文字幕| 亚洲av成人av| 99热这里只有是精品50| 淫秽高清视频在线观看| 啦啦啦韩国在线观看视频| 国产精品爽爽va在线观看网站| 精品国内亚洲2022精品成人| 日本免费a在线| 亚洲国产日韩欧美精品在线观看| 成人漫画全彩无遮挡| 精品久久久久久久久久久久久| 午夜福利网站1000一区二区三区| 18禁裸乳无遮挡免费网站照片| 久久久久免费精品人妻一区二区| 97在线视频观看| 中文字幕亚洲精品专区| 亚洲av中文av极速乱| 久久久久九九精品影院| 成人午夜精彩视频在线观看| 女人十人毛片免费观看3o分钟| 老女人水多毛片| 免费高清在线观看视频在线观看| 大片免费播放器 马上看| av在线老鸭窝| 日韩欧美三级三区| 嘟嘟电影网在线观看| 伊人久久精品亚洲午夜| 日韩欧美 国产精品| 国产av国产精品国产| 夜夜看夜夜爽夜夜摸| 一级毛片我不卡| 91午夜精品亚洲一区二区三区| 女人被狂操c到高潮| 日韩av在线免费看完整版不卡| 久久久久久久大尺度免费视频| 国产av在哪里看| 亚洲图色成人| 婷婷六月久久综合丁香| 精品久久久噜噜| 在线观看美女被高潮喷水网站| 日日摸夜夜添夜夜添av毛片| 毛片一级片免费看久久久久| 国产高清不卡午夜福利| 久久韩国三级中文字幕| 国产精品久久久久久久久免| 久久鲁丝午夜福利片| 99久久中文字幕三级久久日本| 国产91av在线免费观看| 国产真实伦视频高清在线观看| 建设人人有责人人尽责人人享有的 | 欧美激情国产日韩精品一区| 国产成人91sexporn| 国产美女午夜福利| 亚洲在久久综合| 久久综合国产亚洲精品| 中文字幕人妻熟人妻熟丝袜美| 婷婷六月久久综合丁香| 国产片特级美女逼逼视频| av福利片在线观看| 欧美一区二区亚洲| 亚洲无线观看免费| 欧美成人一区二区免费高清观看| 乱码一卡2卡4卡精品| 国产成人精品婷婷| 中文字幕免费在线视频6| 在线观看av片永久免费下载| 免费观看性生交大片5| 亚洲不卡免费看| av网站免费在线观看视频 | 日本色播在线视频| 日韩亚洲欧美综合| 久久这里有精品视频免费| 少妇被粗大猛烈的视频| 日韩亚洲欧美综合| h日本视频在线播放| 久久久久久国产a免费观看| 在线 av 中文字幕| 国产精品综合久久久久久久免费| 久久久久免费精品人妻一区二区| 在线观看美女被高潮喷水网站| 不卡视频在线观看欧美| 日韩成人av中文字幕在线观看| 久久久精品免费免费高清| 夫妻午夜视频| 国产91av在线免费观看| 精品国产三级普通话版| 看黄色毛片网站| 一级毛片电影观看| 欧美日韩综合久久久久久| 久久久久九九精品影院| 天堂中文最新版在线下载 | 亚洲av中文字字幕乱码综合| a级毛片免费高清观看在线播放| 精品亚洲乱码少妇综合久久| 中文欧美无线码| 成年免费大片在线观看| 永久网站在线| 午夜激情久久久久久久| 赤兔流量卡办理| 综合色av麻豆| 精品人妻一区二区三区麻豆| av又黄又爽大尺度在线免费看| 又爽又黄无遮挡网站| 亚洲av日韩在线播放| 校园人妻丝袜中文字幕| 久久久精品免费免费高清| 亚洲av成人精品一区久久| 91av网一区二区| 亚洲综合精品二区| 联通29元200g的流量卡| 一级二级三级毛片免费看| 国产色爽女视频免费观看| 久久久久久久久大av| 乱码一卡2卡4卡精品| 黄片无遮挡物在线观看| 国产中年淑女户外野战色| 日韩av在线免费看完整版不卡| 国产在视频线精品| 日日摸夜夜添夜夜添av毛片| 精品国产三级普通话版| 中文在线观看免费www的网站| 午夜日本视频在线| 一区二区三区免费毛片| 尤物成人国产欧美一区二区三区| 日韩人妻高清精品专区| 欧美成人一区二区免费高清观看| 99热这里只有是精品在线观看| 亚洲精品日韩av片在线观看| 午夜激情欧美在线| av在线蜜桃| 99热6这里只有精品| 午夜福利网站1000一区二区三区| 国产人妻一区二区三区在| 熟女电影av网| 国产黄片视频在线免费观看| 亚洲成色77777| 一本久久精品| 黄片wwwwww| 亚洲精品第二区| 我的女老师完整版在线观看| 国产一区二区三区综合在线观看 | 最新中文字幕久久久久| 天美传媒精品一区二区| 国产69精品久久久久777片| av.在线天堂| a级一级毛片免费在线观看| 欧美一级a爱片免费观看看| 国产精品.久久久| 又黄又爽又刺激的免费视频.| kizo精华| 亚洲一级一片aⅴ在线观看|