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

    The mechanism study of low-pressure air plasma cleaning on large-aperture optical surface unraveled by experiment and reactive molecular dynamics simulation

    2022-07-13 00:36:58YuhaiLI李玉海QingshunBAI白清順YuhengGUAN關(guān)宇珩HaoLIU劉昊PengZHANG張鵬BuerlikeBATELIBIEKE巴特勒別克RongqiSHEN沈榮琦LihuaLU盧禮華XiaodongYUAN袁曉東XinxiangMIAO苗心向WeiHAN韓偉andCaizhenYAO姚彩珍
    Plasma Science and Technology 2022年6期
    關(guān)鍵詞:韓偉張鵬巴特勒

    Yuhai LI (李玉海),Qingshun BAI (白清順),Yuheng GUAN (關(guān)宇珩),Hao LIU (劉昊),Peng ZHANG (張鵬),3,?,Buerlike BATELIBIEKE (巴特勒·別克),Rongqi SHEN (沈榮琦),Lihua LU (盧禮華),Xiaodong YUAN (袁曉東),?,Xinxiang MIAO (苗心向),Wei HAN (韓偉) and Caizhen YAO (姚彩珍)

    1 School of Mechatronics Engineering,Harbin Institute of Technology,Harbin 150000,People’s Republic of China

    2 Research Center of Laser Fusion,China Academy of Engineering Physics,Mianyang 621900,People’s Republic of China

    3 Chongqing Research Institute,Harbin Institute of Technology,Chongqing 401135,People’s Republic of China

    Abstract Low-pressure air plasma cleaning is an effective method for removing organic contaminants on large-aperture optical components in situ in the inertial confinement fusion facility.Chemical reactions play a significant role in plasma cleaning,which is a complex process involving abundant bond cleavage and species generation.In this work,experiments and reactive molecular dynamics simulations were carried out to unravel the reaction mechanism between the benchmark organic contaminants of dibutyl phthalate and air plasma.The optical emission spectroscopy was used to study the overall evolution behaviors of excited molecular species and radical signals from air plasma as a reference to simulations.Detailed reaction pathways were revealed and characterized,and specific intermediate radicals and products were analyzed during experiments and simulation.The reactive species in the air plasma,such as O,HO2 and O3 radicals,played a crucial role in cleaving organic molecular structures.Together,our findings provide an atomic-level understanding of complex reaction processes of low-pressure air plasma cleaning mechanisms and are essential for its application in industrial plasma cleaning.

    Keywords: organic contaminants,large-aperture optical components,low-pressure air plasma,plasma cleaning,reactive species,reactive molecular dynamics

    1.Introduction

    Over the past few decades,inertial confinement fusion (ICF)as a controlled fusion is expected to bring possible solutions to the energy shortages and environmental pollution on the Earth [1].Clean and sustainable energy can be obtained by driving the fusion through multiple large-aperture optical components (>400 mm) in ICF facilities,which transfer and focus high-energy laser beams [2].At present,the contamination damage of large-aperture optical components is the bottleneck limiting the output energy improvement of the ICF facility [3].The surface quality and optical properties of the contaminated optical components are seriously degraded[4].Organic contaminants are an essential component of surface contaminants of optical components,mainly from vacuum pump oil and rubber components.According to the detection of contaminants in the ICF facility [3],organic contaminants were adsorbed on the surface of large-aperture optical components after serving in a vacuum environment for a long time.

    Currently,ex situ wet cleaning (deionized water,acid treatments,ultrasonic,megabound cleaning,etc) of organic contaminants is generally used for large-aperture optical components after manufacturing and coating [3].Due to the limitation of the disassembly and installation of the largeaperture optical components in the ICF facility,it is challenging to remove organic contaminants in situ using existing wetting cleaning methods.At the same time,in the long-time vacuum operation environment of the ICF facility,the surface of optical components will be repeatedly contaminated.Thus,it is urgent to explore new technology to clean large-aperture optical components in situ in the ICF facility.Low-pressure air plasma cleaning (LPAPC) is a promising alternative to removing organic contaminants.Traditional dry cleaning technology (such as laser [5],UV/O3[6],CO2/SNOW [7],etc) is not considered for in situ cleaning of optical components due to low cleaning efficiency and uncontrolled cleaning process.With the development of plasma technology,low-pressure plasma can realize large-area uniform discharge,high-efficiency cleaning,no damage and secondary contamination to the optical surface in the vacuum environment[8].Low-pressure plasma is considered an ideal in situ cleaning technology for large-aperture optical components,which has been initially applied in ITER[9,10],synchrotron beamline optics [11],extreme ultraviolet (EUV) [12],microelectronics[13]facilities.Contaminants such as carbon,Sn and Si produced a reduction reaction with hydrogen plasma to generate CHx,SnH4and SiH in synchrotron beamline optics,EUV facilities.Inert gases (such as helium,argon,and neon) plasmas removed beryllium and tungsten particles from the surface of optical components in ITER by physical bombardment [14].

    Although plasma cleaning is a mature commercial technology,few studies elucidate the interaction between plasma and organic contaminants.To understand the plasma cleaning mechanism,including reaction frequency factor,activation energy,cleaning rate,the dissociation evolution characteristics of organic contaminants were studied using reaction kinetics [15] and finite element analysis [16].Although Arrhenius’s theorem in reaction kinetics described the coupling of the activation reaction and energy of various reactive species in the plasma quantitatively,it cannot provide pathway on cleavage and forming chemical bonds in the detailed reaction process at an extremely short period [17].It is still challenging to fully understand the fundamental chemical reactions between the plasma and organic contaminants and how reactive species in air plasma affect cleavage behaviors by experimental detections.Recently,reactive molecular dynamics (RMD) based on a reactive force field (ReaxFF)simulation have been performed to clarify the mechanisms of interactions between plasma-generated species and various biomolecules[18,19].It could provide insights into the entire interwoven reaction pathway of the plasma cleaning at the atomic level compared with the traditional theoretical calculation that researches a single reaction [20].

    This study conducted an in situ low-pressure air plasma cleaning experiment on the large-aperture optical components to remove the organic contaminant.Reactive species and excitation signals in low-pressure air plasma were studied using optical emission spectroscopy (OES).Then,the chemical composition and components content of the surface organic contaminants were also analyzed by attenuated total reflection flourier transformed infrared spectroscopy (FTIR)and X-ray photoelectron spectroscopy (XPS) to provide a reference to simulation results.We aim to understand the plasma cleaning mechanism by RMD simulation,which provides a possible reaction pathway for the dissociation of organic contaminants by reactive species in the plasma.

    2.Experimental methods and simulation details

    2.1.Experimental methods

    SiO2sol-gel antireflective (AR) film is widely coated on the surface of large-aperture optical components to improve the transmittance and laser-induced damage threshold(LIDT)[2].In this experiment,the fused quartz samples (50×50×5 mm3)were prepared with the same surface quality as the largeaperture optical components in the ICF facility.According to previous research [21],the sol-gel AR film was deposited on the fused quartz substrates using the dip-coating method at 355 nm.The fused quartz was immersed in the sol-gel liquid and then pulled at a constant speed(~80 mm min-1).After the post-treatment of hexamethyldisilazane and ammonia,the dense film was obtained,consistent with the actual coating method of large-aperture optical components [22].

    The samples were placed in a vacuum chamber to contaminate for 5 h with a high vacuum(10-5Pa)to simulate the organic contamination environment in the ICF facility.Dibutyl phthalate(DBP)has the benzene ring,ester group and carbon chain,representing the chemical properties of most organic contaminants (alkanes,aromatic hydrocarbons,esters),which are found in the high energy facilities [23].Therefore,DBP was chosen as a typical contaminant in our research.

    Figure 1 outlines the experimental equipment designed to verify in situ plasma cleaning of large-aperture optical components,which consisted of a power system,laser path,vacuum system and four electrodes(made of copper plates of 20×5×430 mm3).A large area of uniform plasma was generated by discharging between electrodes and the cavity wall in the vacuum chamber.The air was excited by an audio frequency power supply and the plasma was dispersed on the optical surface (430×430 mm2).The cleaning parameters are listed in table 1.The discharges use frequencies usually in the radio frequency band,such as 13.56 MHz.The plasma sheath voltage can reach hundreds or even thousands of volts in a capacitively coupled discharge plasma reactor.When the ions pass through the sheath,they will be accelerated by the sheath voltage to bombard the surface of the optical component with high energy,which is easy to cause surface damage.According to previous research on grating cleaning[24],oxygen plasma is overaggressive and challenging to control accurately.Thus,we choose more mild air as the plasma gas source to avoid surface damage to large-aperture optical components.At the same time,the audio-frequency(20 kHz) discharge structure is straightforward and suitable for mounting on the cavity wall of the laser path with limited space.

    Figure 1.Experiment setup for low-pressure air plasma cleaning.

    Figure 2.Molecular model of a chemical reaction between plasma and organic contaminants.

    Figure 3.(a) Transmittance and (b) LIDT of optical components.

    Figure 4.The excited species in air plasma.(a) The full optical emission spectra of air plasma,(b) an emission line of O I.

    Figure 5.FTIR spectra of sol-gel AR film with different treatments.

    Table 1.Low-pressure plasma cleaning parameters.

    The transmittance of optical components coated with the sol-gel AR film is measured by the spectrophotometer(Perkin Eimer Lambda 950).Single-pulse Nd:YAG laser(wavelength 355 nm) measures the LIDT of optical components by the 1-on-1 method.Optical emission spectra (Maya 2000 pro)were used to measure reactive species and their concentration during air plasma cleaning.The measurement range of spectrum signal was from 165 to 1100 nm and the optical resolution can reach approximately 0.035 nm.The atomic force microscope (PSIA XE-100) was used to measure the surface morphology and roughness of optical components.A total reflection accessory of FTIR(Thermo Nicolet 5700)was conducted to examine the formation and cleavage bonding information after plasma cleaning.It is challenging to detect infrared absorption spectra of tiny organic contaminants on the surface of SiO2sol-gel AR film.Thus,the sample surface was coated with a thin layer of DBP.An XPS(Thermo Fisher ESCALAB 250Xi)was conducted to analyze surface element content and functional groups quantitatively.XPS spectra were obtained using an Al Kα (1486.7 eV) X-ray source at 300 W.High-resolution spectra (0.1 eV) were calibrated by C(1s) signal at 285.0 eV [25].

    2.2.Simulation details

    Plasma cleaning simulation is based on RMD,a bond orderbased potential proposed by van Duin [26].ReaxFF uses a general relationship between the bond distance and bond order,bond energy,leading to proper dissociation of bonds to separated atoms.Although RMD does not explicitly simulate electrochemical reactions involving electrons (e-) transfer,it is beneficial to construct the entire interwoven chemical reactions in the plasma cleaning and the mass transfer in the femtosecond time range [20].In RMD simulations,the general system energy is divided into various partial energy contributions.Each expression of the bond order function is demonstrated in equation (1).

    Here,Esystemrepresents the system energy.The Ebond,Eover,Eunder,Eval,Epen,Etors,Econj,EvdWaals,ECoulombdenote the bond energy,over-/under-coordination energy,valence angle energy,penalty energy,torsion energy,conjugation energy,non-bonded van der Waals interactions and Coulomb interactions,respectively.A ReaxFF potential comprising C/H/O parameters is employed to accurately describe the plasma cleaning model.

    In this study,plasma concentration refers to the concentration of reactive species involved in chemical reactions.Plasma temperature refers to the gas temperature,which is mainly determined by the average kinetic energy of the ions during elastic collisions.According to previous experimental results,plasma composition is mainly considered reactive oxygen species (e.g.,O,OH,O3,H2O2,and HO2) compared to reactive nitrogen species [27].O2is regarded as the main component of the gas source for comparing destruction capability with other exciting species.Therefore,five DBP molecules and two hundred reactive species were randomly placed in a sealed box,as shown in figure 2.The plasma concentration was much higher than the actual plasma concentration,but this setting was widely accepted to explain the chemical reaction mechanism between plasma and organic contaminants from a microscopic perspective [18].Periodic boundary conditions were applied to simulate an infinite system.Thus,energy minimization was performed by a smart minimizer to make the model more reasonable.After that,the model was simulated in the canonical ensemble (NVT) at room temperature for 1 ns and the time step was set at 0.25 fs.

    3.Results analysis

    3.1.Optical properties of sol-gel AR films

    The ultimate goal of low-pressure air plasma in situ cleaning is to restore the optical properties of large-aperture optical components.The transmittance of the sol-gel AR film decreased rapidly from about 99.80% to 99.38% after contamination for 10 h,as shown in figure 3(a),which was unacceptable for high-energy laser transmission.It was noted that the transmittance of the optical component recovered after 5 min of plasma cleaning,which indicated that organic contaminants were removed by the plasma from the surface.The LIDT of optical components in the high-energy laser facility is an essential factor to limit the improvement of the load capacity ability of the facility [2].The LIDT of the optical component decreased by nearly one-third(from ~13.6 to ~10.4 J cm-2)after 10 h of contamination,with significant ablations in the damaged area,as shown in figure 3(b).This is because organic contaminants are adsorbed in the pores of sol-gel AR films and carbonized under high-energy laser irradiation.After 5 min of plasma cleaning,the LIDT was recovered and ablation regions caused by organic contaminants were eliminated.Although the optical properties of the components were degraded by organic contamination in the high vacuum environment,the optical properties of the components recovered completely after plasma cleaning.

    3.2.Optical emission spectra of air plasma

    OES analysis of air plasma is conducted in situ in different regions within the window to obtain reactive species information,as shown in figure 4.OES is an effective method to analyze excited species,vibrational temperature,electron density and temperature during plasma excitation,which has no disturbance to the plasma compared to the Langmuir probe technique in the in situ diagnosis[28,29].It can be found that there are mainly sharp emission peaks (like H,O,N2,CO,N2+,C2,C3) and wide emission bands (like N2,) by measuring the full spectrum of low-pressure air plasma,which is in good accordance with previous results [30-32].Among them,the signal of O I was found in spectra,which was the main reactive species driving chemical reaction and induced the production of other species (e.g.,OH,O3,H2O2,and HO2).

    In addition,we found that air plasma concentration distribution in the in situ plasma cleaning was different in the large-aperture window according to purple intensity in figure 4(a).So the window was divided into nine regions to measure the emission spectra of O I.The plasma concentration was positively related to the spectral intensity,according to equation (2) [32].

    where nAis the concentration of reactive species A,Iλis the light emission intensity integrated against line width,and αλAis the coefficient related to electron velocity and emission wavelength.

    Among them,the most typical emission line of O I was 777.9 nm (3p5P→3s5S),which was a direct excitation from the ground state[33].The line 844.6 nm(3p3P→3s3S)of O I was not used to evaluate plasma excitation due to a weak signal.The plasma concentration distribution showed a high concentration of O radical around the window and a low concentration in the middle,consistent with the diffusion law of charged particles in plasma [32].The highest O I radical spectral intensity in region 9 was ~4540.5,twice that in region 5,representing the higher excitation degree of air in region 9.The plasma species spread toward the middle region of low density through the capacitive coupling charge from the four sides of the window.It is necessary to study the influence of the contributions of air plasma concentration on the plasma cleaning effect in future research.

    3.3.Chemical composition evolutions on optics surfaces

    To understand the composition of chemical elements on the surfaces of sol-gel AR film,the surface chemical bond information and element contents were detected by FTIR and XPS,as shown in figures 5 and 6.Figure 5 shows FTIR spectra of cleavage and formation bond information before and after plasma cleaning.It can be noted that abundant C-H bonds were found on the surface of sol-gel AR film after organic contamination.In DBP molecular structure,1078 cm-1and 1123 cm-1correspond to the in-plane swing deformation vibration of the benzene ring and the stretching vibration of the C-O single bond near the C=O double bond.1281 cm-1refers to the in-plane swing vibration of the benzene ring and the stretching vibration of the C-C bond between the C atom in the benzene ring and C=O.2859 cm-1and 2924 cm-1represent the stretching vibrations of the C-H bond on CH3and CH2,dropping below the untreated level after LPAPC.The results showed that plasma could remove most contamination from the surface in situ.Besides,2360 cm-1and 3500-3950 cm-1correspond to CO2and H2O(gaseous) in the background environment,which increased significantly relative to the untreated surface.The plasma is widely used to clean surfaces containing trace amounts of organic contaminants,and it is still difficult to completely remove contaminants coated with a layer of DBP.

    Figure 6.High-resolution XPS spectra of the film before and after low-pressure air plasma cleaning.(a) The full spectrum of XPS after plasma cleaning for 5 min,(b)element contents,(c)fitting of the C(1s)signal after plasma cleaning for the 30 s,(d)fitting of the O(1s)signal after plasma cleaning for 5 min.

    Furthermore,both C-C (1281 cm-1) and C-O(1123 cm-1) bonds were easily oxidized to C=O and C=C double bonds after plasma cleaning[34],indicating that DBP molecular structure was destroyed.A wide range of C=C was stretching vibration frequency at 1370-1700 cm-1,mainly the π bond in the benzene ring and olefin,with a vibration frequency of 1370-1610 cm-1and 1620-1700 cm-1,respectively.These results indicated that unsaturated bonds were produced during the cleaning process.Different types of C=O bonds (such as lipid group,aldehyde group) were generated in the peak positions of 1650-1800 cm-1,which was due to the intermediate oxidation products of the carbon chain,which had been proved by Bonova et al [35].Therefore,plasma removed organic contaminants by destroying CH and C-C bonds and oxidizing to generate small molecular groups such as C-O,C=O,C=C,O-H,CO2and H2O.

    A high-resolution XPS analysis of samples was performed to understand the changes in element contents and physicochemical properties of sol-gel AR film before and after plasma cleaning.Figure 6(a) shows the XPS signals C(1s),O(1s),Si(2 s),Si(2p),Al(2p),Fe(2p) and Cu(2p),which contain the components contents and binding energy.Since the electrode(copper plate),the cavity wall(aluminum alloy)and the screw(stainless steel) are made of metal.The elements of Al(2p),Fe(2p),and Cu(2p) were not found in all XPS signals after plasma cleaning for 5 min in figure 6(a),which indicated that the surface of sol-gel AR film was not contaminated by the metal during plasma cleaning.However,it is inevitable that the electrode sputter metal particles after multiple cleaning cycles [9].Since the discharge in this device is between the electrode and the cavity wall,the electrode sputtering can be blocked by coating the electrode surface or adding borosilicate glass.In this work,it is proposed for the first time to use audiofrequency capacitive coupled air plasma to remove organic contaminants on the surface of large-aperture optical components in situ.Future research aims to reduce and avoid metal contamination and damage on the surface of optical components after multiple cycles of cleaning.

    Figure 7.Schematic diagram of low-pressure air plasma cleaning mechanism.

    The main elements of the sol-gel AR film surface contaminated by DBP were C(1s)and O(1s).Among them,C(1s)content increased rapidly with the increase of contamination time because the main component of DBP was C(1s),as shown in figure 6(b).As the surface of sol-gel AR films consisted of O(1s) and Si,the content was reduced after contamination because of forming a layer of oil film on the surface.The organic contaminants were not completely removed due to a slight C(1s)decrease after plasma cleaning for the 30 s.After cleaning for 5 min,C(1s)content was lower than that of the uncontaminated surface and the content of O(1s)and Si reached the uncontaminated level.Furthermore,the C(1s)and O(1s) element signals were analyzed to quantitatively evaluate the chemical composition changes of the products during the plasma cleaning.Figure 6(c)shows that the C(1s)represents the middle product during the chemical reaction after plasma cleaning for 30 s,characterized by having three well-defined components,284.8 eV,286 eV,and 289 eV,respectively [36].These signals indicated that large amounts of C-O,C=O bonds were produced during the plasma cleaning due to plasma oxidation,consistent with the FTIR results.Similarly,it was possible to assign three components to the O(1s) signal placed at 531.9 eV,532.9 eV and 533 eV,corresponding to oxygen atoms in different chemical states after plasma cleaning for 5 min O(1s)had a characteristic signal of 531.9 eV,assigned to oxygen atoms bonded with hydrogen forming hydroxide groups(O-H).At the same time,there were still more C=O bonds left after plasma cleaning,which indicated that C=O was not easy to decompose and formed the final residuals.

    Figure 8.The interaction characteristics between reactive species and DBP in the reaction process.(a) Total energy evolution,(b) mean square displacement of reactive species.

    4.Discussion

    4.1.Analysis of low-pressure air plasma cleaning mechanism

    When an audio-frequency alternating current electric field was applied to the discharge electrodes,free electrons were accelerated in the electric field and gained kinetic energy.Electrons with certain energy collided with atoms and molecules in the gas in the process of movement.In inelastic collisions,electrons transferred energy to atoms and molecules to ionize them and produce abundant reactive species,as shown in figure 7.When these particles accelerated through the sheath to bombard the solid surface,they transferred energy to organic contaminants and substrates,facilitating the desorption of organic contaminants on the surface.Among them,secondary electron emission has a significant influence on the generation and concentration of reactive species,but it has no influence on the dissociation mechanism of organic contaminants by reactive species.Previous studies revealed that oxidation and ultraviolet reaction,ion and electrons collision,magnetic fields attraction mainly occurred during plasma treatment,chemical reaction played a crucial role in plasma cleaning[27].In plasma cleaning,ions energy is usually reduced to avoid surface damage,so the contribution of physical bombardment to plasma cleaning can be ignored.

    A layer of oil film was formed on the sol-gel AR film after contamination,as the AFM detection shows in figure 7.Reactive species in low-pressure air plasma combined with organic contaminants to form new unstable radicals.The FTIR and XPS test results proved that the organic contaminants on the surface were removed with significant destruction of C-H and C-C bonds after plasma cleaning.Further,those chemical bonds were oxidized to C=C,C=O and O-H bonds during plasma cleaning.New free radicals were generated on the surface,along with the formation of small molecules.These chemical reactions continued until the hydrocarbons were decomposed entirely into small molecules such as H2O and CO2,which have been examined in previous research [37].Then,these volatile small molecules were extracted easily from the surface by a vacuum pump due to weak van der Waals and Coulomb forces between small molecules and the surface.After plasma cleaning,the surface morphology and roughness of sol-gel AR film were restored and showed the loose porous structure without damage,as shown in figure 7.Finally,organic contaminants on the surface of large-aperture optical components were removed by in situ plasma cleaning.

    Although the reactive species and the intermediate products of plasma cleaning can be found through various experimental detection methods,it is still challenging to obtain the evolution behaviors and formation mechanisms of chemical bonds.Especially,some functional groups,such as C=O(carboxylic acid,ester,aldehyde,ketone,in 1650-1800 cm-1) and C=C (olefins,aromatic hydrocarbons,in 1370-1700 cm-1),have overlapping wavelengths making it impossible to distinguish the product species.Thus,RMD was conducted to study the evolution of cleavage and the formation of chemical bonds during plasma cleaning at the atomic level.

    4.2.Reaction data analysis

    The reaction statistics between air plasma and organic contaminants were analyzed further to elucidate the chemical reaction characteristic at the atomic level.Figure 8(a)shows the total energy change process of a chemical reaction within 1 ns during RMD simulation.Although the total energy of the system had not reached equilibrium,the primary chemical reaction was complete in plasma cleaning,including richer reaction information.When reactive species were involved in plasma cleaning,a cascade of bond cleavage and formation events within approximately 0.5 ns reached the reaction scope.This phenomenon was that the plasma concentration in the system decreased and the reaction gradually stopped as the plasma reacted with organic contaminants.O2molecules did not destroy the DBP structure because they exhibited weak,non-bond interactions during plasma cleaning.Notably,O,O3,and H2O radicals reacted with the DBP like H2O2and OH radicals,with the only difference being that O,O3,and H2O radicals reacted violently with DBP in 0.2 ns,resulting in a significantly decreased total energy.

    Figure 9.Temporal snapshots of hydrogen bond formation between atoms.The dashed blue lines are hydrogen bonds.

    Figure 10.The cleaning capability of reactive species was evaluated.(a)The cleaning capability of different reactive species in air plasma is considered within 1 ns,which defines the cleavage of C-C bonds compared to total C-C bonds,(b)statistics of C-H bond dissociation and C=O,C=C,and C-O bond formation upon interaction with reactive species and organic contaminants.The percentages are calculated by breaking and forming bonds divided by the total bonds present in the model.

    The translational mobility of reactive species during plasma cleaning was examined by the mean squared displacement (MSD).This work calculated the MSD of reactive species based on each species’ center of mass,as shown in figure 8(b).The slope of the MSD curve is related to the mobility of reactive species,which means greater mobility of reactive species with a higher slope [21].The dynamic properties of these reactive species were related to the chemical characteristics and density of reactive species.Oxygen atoms abstracted hydrogen atoms from contaminants to form OH radicals,which tended to form long chains due to the formation of hydrogen bonds,as shown in figure 9.The concentration of OH generated in the oxygen radicals reaction model was less than that of other models,which indicated that oxygen radicals had the most significant mobility characteristics.O2,O3,and HO2had similar mobility characteristics due to similar molecular densities after dissociation.Many hydrogen bonds were formed between reaction products of H2O2,which resulted in the strong interaction between molecules.Therefore,O radical had better mobility characteristics during chemical reactions and reacted more sufficiently with organic contaminants during long-time plasma cleaning.

    Furthermore,statistics regarding hydrogen abstraction reactions and carbon chains oxidation were conducted after the chemical reaction.Figure 10 shows that the total number of broken and formed bonds is calculated based on a simulation time of 1 ns.Figure 10(a)evaluates the evolution of the cleaning capability of reactive species in the air plasma along with reaction time,which can be divided into three stages.The reaction between the reaction species and the organic contaminants was the most intense in stage I(within 0.05 ns).The reactive species had a similar strong cleaning capability except for the H2O2radicals.With the consumption of reactive species in the system,the reaction rate gradually decreased from 0.05 ns to 0.35 ns in stage II.Eventually,the chemical reaction basically stopped at stage III after the plasma concentration dropped to a certain level.The reaction between OH and DBP continued throughout the whole reaction process,while the cleaning capability of OH was slightly weaker than those of O,O3,and HO2.O radicals reacted most strongly with organic contaminants and had the most vital cleaning capability,followed by O3and HO2.High plasma concentration was adopted in the simulation to shorten the reaction time,which resulted in the inconsistency between the reaction time in the experimental time and the simulation,as shown in figures 6(b) and 10(a).

    Figure 11.Real-time dynamics of air plasma cleaning.(a)Snapshots of the partial cleaning system extracted from the RMD simulation of air plasma cleaning.The red color represents new generated molecules and bonds.The gray color represents the molecules that were not involved in the reaction.(b) Time dependences of the numbers of oxygen and main molecular species in real-time RMD simulation.

    Figure 12.Temporal snapshots DBP dissociation mechanism induced by O radicals.(a)Dissociation of DBP,(b)O-H bonds formation,(c)and (d) C-C bonds cleavage,(e) DBP dissociated products,(f) main reaction pathways in the plasma cleaning.

    In addition,the effects of reactive species on chemical bond breaking and formation after plasma cleaning were evaluated.HO2radicals triggered more hydrogen abstraction reactions (~39%) than other reactive species,indicating that the HO2radicals played a more crucial role in the bond cleavage of C-H,followed by O and O3radicals.Meanwhile,this characteristic was also observed in the capability of HO2atoms to form C=O bonds (~16%).The number of C-O bonds was little after plasma cleaning because C-O bonds were unstable and easily oxidized to C=O because of high plasma concentration.However,more C=C bonds were found in the FTIR and XPS results due to low plasma concentrations in the actual reaction process.Therefore,abundant C-O,C=C,and C=O bonds were generated in air plasma cleaning,in which O,O3,and HO2played a leading role compared with OH and H2O2radicals.

    4.3.Possible reaction pathway

    The plasma cleaning reaction has a similar pathway through breaking the molecular structure into small groups between various reactive species.As the representation of chemical reactions,figure 11 shows the time-dependent progression of O radicals and the main molecular species during the RMD simulation.After 1 ns,about 120 O were consumed,and about 45 C-H and 25 C-C were destroyed,which indicated that the DBP molecule structure had already been destructed except benzene ring.At the same time,about 110 O-H,20 C=C,and 30 C=O were generated after the chemical reaction.DBP molecules were destructed gradually slightly with the O consumed quickly,especially the chemical reaction almost stopped after 0.4 ns.Because oxygen concentration progressively reduced along with chemical reaction,the carbon chain had been oxidated into C=C and C=O.It was challenging to have further oxidation of benzene rings.Therefore,the final products of plasma cleaning were benzene ring,C=C,and C=O,consistent with the XPS results.

    The primary reaction was hydrogen abstraction from the carbon chain,which triggered the hemolytic cleavage of important bonds (i.e.,C-C and C-H bonds) and formed double bonds(i.e.,C-O,C=O,and C=C bonds).Therefore,cleavage and formation of mechanisms of DBP are discussed with statistics of O radicals.Plasma cleaning is a complex chemical reaction that involves abundant reaction pathways and products during the process.An O radicals-based DBP dissociation main pathway is illustrated in figure 12.Most of the reaction species in each figure were removed to analyze the products more clearly during the chemical reaction.Due to the high concentration of O atoms and the intense oxidation,the C-O bond on carbon 1 and carbon 5 was broken firstly,which indicated that C-O-C in the ester structure of the DBP molecule was the most easily destroyed.At the same time,two hydrogen atoms on carbon were taken away by O to form O-H bonds,which further abstracted hydrogen from the carbon chain to become an H2O molecular (see green lines,blue dash circles in figures 12(a)and(b)).Due to the high concentration of oxygen radicals,OH radicals were further oxidized into groups such as HO2and H2O2.Subsequently,carbon 1 and carbon 5 connected two oxygen atoms,one of which formed a C=O double bond.The hydrogen atoms were abstracted away from the carbon chain,resulting in C≡O(shè) and C=O bonds forming.Further,the carbon chain length was broken down (green line,blue dashed circle in figures 12(c) and (d)).Therefore,these changes led to the complete oxidation and destruction of four CO2molecules of a carbon chain (blue dashed circle in figure 12(e)).Finally,the two DBP branched chains were oxidized into CO2and OH radicals except for the benzene ring,which completely dissociated organic contaminants by plasma.The main reaction pathways of plasma cleaning are provided in figure 12(f).In the actual reaction process,there is a large amount of intermediate residue on the surface of the reaction due to the limitation of plasma concentration and cleaning time.The benzene ring is difficult to destroy in the plasma cleaning process,but the removal efficiency can be improved with increasing temperature and cleaning time [38,39].

    5.Conclusion

    Experiments demonstrated that the organic contaminants on the surface of large-aperture components (430×430 mm2)were removed by low-pressure air plasma in situ without secondary contamination.Abundant reactive species were detected in optical emission spectra during plasma cleaning.Unsaturated bonds,such as C-O,C=C,C=O,and O-H bonds,were generated during the cleaning process and C=O was the final production.Besides,the cleaning mechanism of the organic contaminants with various reactive species in air plasma was illustrated by experiments and the reactive molecular dynamics simulation.In the process of plasma cleaning,O,HO2,and O3radicals had a more robust cleaning capacity in dissociating organic contaminants than H2O2and OH.O radicals had the most cleaning ability due to significant mobility and reacted more sufficiently with organic contaminants during the reaction process.The rate of chemical reaction slowed down gradually with the decrease of the plasma concentration.Reactive molecular dynamics simulations revealed the main pathway that the DBP molecular structure was oxidized into small groups by oxygen radicals and elucidated the dissociation and formation mechanism of the chemical bonds.The primary reaction pathway was hydrogen abstraction from the carbon chain,which triggered the hemolytic cleavage of important bonds (i.e.,C-C and C-H bonds) and formed unsaturated groups (i.e.,C-O,C=C and C=O),consistent with the experimental results.Aromatic hydrocarbons are more difficult to dissociate than alkanes and remain on the surface after plasma cleaning.

    Acknowledgments

    This work was supported by the Joint Funds of National Natural Science Foundation of China and China Academy of Engineering Physics (NSAF) (No.U2030109) and National Natural Science Foundation of China (No.52075129).

    猜你喜歡
    韓偉張鵬巴特勒
    Quantitative analysis of the main components in ceramic raw materials based on the desktop LIBS analyzer
    Competitive effect between roughness and mask pattern on charging phenomena during plasma etching
    基于PLC控制的平移式自動(dòng)門設(shè)計(jì)
    理發(fā)風(fēng)波
    “律法之前”與“刑具之后”:巴特勒律法述行思想的文學(xué)闡釋
    事不關(guān)己
    NBA特刊(2017年17期)2017-11-11 07:44:07
    乖乖蛋,不搗亂
    特別的外宿
    UAV Velocity Measurement for Ground Moving Target
    (口歐)!鷹笛
    精品一区二区三区视频在线观看免费 | 80岁老熟妇乱子伦牲交| 九色亚洲精品在线播放| 人人澡人人妻人| 性少妇av在线| 极品人妻少妇av视频| 欧美日韩福利视频一区二区| 90打野战视频偷拍视频| 长腿黑丝高跟| 亚洲,欧美精品.| 日韩 欧美 亚洲 中文字幕| 免费看a级黄色片| 女人被躁到高潮嗷嗷叫费观| 在线观看www视频免费| 9热在线视频观看99| 中文字幕色久视频| 99热国产这里只有精品6| 黄频高清免费视频| 91成人精品电影| 悠悠久久av| 精品少妇一区二区三区视频日本电影| 咕卡用的链子| 高清在线国产一区| 国产成人av教育| 中亚洲国语对白在线视频| 亚洲av熟女| 在线看a的网站| 又紧又爽又黄一区二区| 成熟少妇高潮喷水视频| 精品久久蜜臀av无| 婷婷丁香在线五月| 操出白浆在线播放| 午夜精品在线福利| 一级片免费观看大全| 久久久国产成人免费| 电影成人av| 啪啪无遮挡十八禁网站| 一个人观看的视频www高清免费观看 | 黄频高清免费视频| 性色av乱码一区二区三区2| 午夜影院日韩av| 亚洲欧美激情在线| 在线观看舔阴道视频| 免费在线观看日本一区| 国产成人一区二区三区免费视频网站| 人成视频在线观看免费观看| 90打野战视频偷拍视频| 欧美黄色片欧美黄色片| 男女下面进入的视频免费午夜 | 一进一出抽搐gif免费好疼 | 欧美激情 高清一区二区三区| 新久久久久国产一级毛片| 国产精品国产高清国产av| 一级毛片精品| 在线播放国产精品三级| 午夜精品久久久久久毛片777| 中国美女看黄片| 精品久久久久久久毛片微露脸| 久久亚洲精品不卡| 91字幕亚洲| 激情视频va一区二区三区| 精品久久久精品久久久| 人人妻人人添人人爽欧美一区卜| 精品欧美一区二区三区在线| 国产精品美女特级片免费视频播放器 | 日韩人妻精品一区2区三区| 国产精品爽爽va在线观看网站 | 一夜夜www| 国产av精品麻豆| 久久天堂一区二区三区四区| 国产深夜福利视频在线观看| 亚洲五月色婷婷综合| 欧美精品啪啪一区二区三区| 亚洲精品av麻豆狂野| www国产在线视频色| 亚洲人成电影免费在线| 国产精品 国内视频| 久久中文字幕一级| 免费av毛片视频| 在线永久观看黄色视频| 国产深夜福利视频在线观看| 99久久人妻综合| 12—13女人毛片做爰片一| 久久国产亚洲av麻豆专区| 亚洲欧美精品综合久久99| 欧美人与性动交α欧美精品济南到| 日日干狠狠操夜夜爽| 国产精品国产av在线观看| 国产xxxxx性猛交| 天天影视国产精品| 亚洲全国av大片| 不卡一级毛片| 天堂俺去俺来也www色官网| 午夜91福利影院| 女人爽到高潮嗷嗷叫在线视频| 在线国产一区二区在线| 啦啦啦在线免费观看视频4| 人人妻人人添人人爽欧美一区卜| 欧美久久黑人一区二区| 亚洲精品av麻豆狂野| 国产1区2区3区精品| 首页视频小说图片口味搜索| 一边摸一边抽搐一进一小说| 欧美成人性av电影在线观看| 国产黄a三级三级三级人| 亚洲熟妇中文字幕五十中出 | 国产国语露脸激情在线看| 国产精品美女特级片免费视频播放器 | 黄片大片在线免费观看| 十八禁网站免费在线| 88av欧美| 亚洲美女黄片视频| 身体一侧抽搐| 夜夜看夜夜爽夜夜摸 | 日韩免费av在线播放| 午夜91福利影院| 黑人操中国人逼视频| 在线十欧美十亚洲十日本专区| 成年版毛片免费区| 伊人久久大香线蕉亚洲五| 亚洲成人精品中文字幕电影 | 国产成人精品无人区| 两性夫妻黄色片| 99热国产这里只有精品6| 久久人妻福利社区极品人妻图片| 丰满的人妻完整版| 欧美黑人精品巨大| 国产真人三级小视频在线观看| 成人三级黄色视频| 深夜精品福利| 在线天堂中文资源库| 男女床上黄色一级片免费看| 午夜免费观看网址| 亚洲国产毛片av蜜桃av| 亚洲自拍偷在线| 日本五十路高清| 男人舔女人的私密视频| 老熟妇仑乱视频hdxx| 黄色 视频免费看| 欧美日韩黄片免| 午夜两性在线视频| 国产熟女xx| 看免费av毛片| 国产高清国产精品国产三级| svipshipincom国产片| 美女扒开内裤让男人捅视频| 校园春色视频在线观看| 高潮久久久久久久久久久不卡| av国产精品久久久久影院| 中文字幕av电影在线播放| 麻豆一二三区av精品| 精品国内亚洲2022精品成人| 美女福利国产在线| 在线国产一区二区在线| 不卡av一区二区三区| 曰老女人黄片| 麻豆成人av在线观看| 亚洲三区欧美一区| 欧美日韩黄片免| ponron亚洲| 新久久久久国产一级毛片| 日本欧美视频一区| 国产又色又爽无遮挡免费看| 午夜福利一区二区在线看| 日本 av在线| 久久精品国产清高在天天线| 久久中文看片网| 国产黄a三级三级三级人| 国产av精品麻豆| 久久精品国产亚洲av高清一级| 国产精品久久久久久人妻精品电影| 母亲3免费完整高清在线观看| 日韩人妻精品一区2区三区| 国产色视频综合| 日本wwww免费看| 精品久久久久久成人av| 精品一区二区三区av网在线观看| xxx96com| 麻豆一二三区av精品| 色播在线永久视频| 欧美精品啪啪一区二区三区| 亚洲成人免费电影在线观看| 亚洲专区中文字幕在线| 久久人人精品亚洲av| 男女做爰动态图高潮gif福利片 | 一本综合久久免费| 无限看片的www在线观看| 亚洲色图av天堂| 中文亚洲av片在线观看爽| 涩涩av久久男人的天堂| 两个人免费观看高清视频| 国产精品日韩av在线免费观看 | 亚洲成av片中文字幕在线观看| 日本黄色日本黄色录像| 久久久久久久精品吃奶| 在线观看午夜福利视频| 一进一出抽搐动态| 老司机福利观看| 麻豆一二三区av精品| 无限看片的www在线观看| 又紧又爽又黄一区二区| 韩国av一区二区三区四区| 成人18禁高潮啪啪吃奶动态图| 热re99久久国产66热| 1024视频免费在线观看| av视频免费观看在线观看| 免费在线观看完整版高清| 一个人免费在线观看的高清视频| 中亚洲国语对白在线视频| 在线观看免费午夜福利视频| 成人影院久久| 免费在线观看视频国产中文字幕亚洲| 在线观看66精品国产| 精品国产一区二区三区四区第35| 成人亚洲精品一区在线观看| 成年人黄色毛片网站| 国产亚洲欧美精品永久| 欧美日韩国产mv在线观看视频| 激情在线观看视频在线高清| 丰满饥渴人妻一区二区三| 美女午夜性视频免费| 嫩草影院精品99| 别揉我奶头~嗯~啊~动态视频| 亚洲男人的天堂狠狠| 精品一区二区三区av网在线观看| 99久久99久久久精品蜜桃| 亚洲av片天天在线观看| 国产精品乱码一区二三区的特点 | 99精品在免费线老司机午夜| 中文亚洲av片在线观看爽| 欧美日韩亚洲高清精品| 长腿黑丝高跟| 电影成人av| 一级,二级,三级黄色视频| 久99久视频精品免费| 国产亚洲欧美精品永久| 99精品在免费线老司机午夜| 亚洲精品成人av观看孕妇| 波多野结衣av一区二区av| 欧美黄色片欧美黄色片| 亚洲精华国产精华精| 亚洲第一欧美日韩一区二区三区| 男男h啪啪无遮挡| 18禁观看日本| 欧美老熟妇乱子伦牲交| 人成视频在线观看免费观看| 欧美一级毛片孕妇| 女性被躁到高潮视频| 亚洲第一av免费看| 露出奶头的视频| 嫁个100分男人电影在线观看| 性欧美人与动物交配| 91九色精品人成在线观看| 波多野结衣高清无吗| 国产欧美日韩精品亚洲av| 麻豆国产av国片精品| 正在播放国产对白刺激| 国产乱人伦免费视频| 亚洲精品在线观看二区| 亚洲熟妇中文字幕五十中出 | 久久久久久亚洲精品国产蜜桃av| 在线国产一区二区在线| 悠悠久久av| 久久精品影院6| 男人的好看免费观看在线视频 | 91大片在线观看| 国产精品99久久99久久久不卡| 欧美av亚洲av综合av国产av| 国产av精品麻豆| 视频区图区小说| av国产精品久久久久影院| 999精品在线视频| 黄色视频不卡| 视频区图区小说| 两个人看的免费小视频| 在线观看66精品国产| av视频免费观看在线观看| www.www免费av| x7x7x7水蜜桃| 黄色怎么调成土黄色| 亚洲精品在线美女| 亚洲成人国产一区在线观看| 中亚洲国语对白在线视频| 日本三级黄在线观看| 麻豆成人av在线观看| 男女下面进入的视频免费午夜 | 色婷婷av一区二区三区视频| 久久青草综合色| 中文亚洲av片在线观看爽| 80岁老熟妇乱子伦牲交| 男女之事视频高清在线观看| 午夜老司机福利片| 热re99久久国产66热| 女人高潮潮喷娇喘18禁视频| 成在线人永久免费视频| bbb黄色大片| 国产成人免费无遮挡视频| 纯流量卡能插随身wifi吗| 久久久久久大精品| 精品无人区乱码1区二区| 色播在线永久视频| 99国产极品粉嫩在线观看| 欧美日韩中文字幕国产精品一区二区三区 | 亚洲精品久久午夜乱码| 亚洲五月色婷婷综合| 1024香蕉在线观看| 久久久久久久久久久久大奶| 亚洲精品国产一区二区精华液| 亚洲色图av天堂| 精品国产美女av久久久久小说| 亚洲精品美女久久av网站| 黄片大片在线免费观看| 制服人妻中文乱码| 999久久久精品免费观看国产| 99热只有精品国产| 69精品国产乱码久久久| 99国产精品99久久久久| 久久久久久人人人人人| 午夜福利在线免费观看网站| 高清毛片免费观看视频网站 | 国产精品野战在线观看 | 午夜精品在线福利| 亚洲在线自拍视频| 不卡一级毛片| 久久久久久久午夜电影 | 欧美日韩亚洲高清精品| 亚洲精品粉嫩美女一区| 免费在线观看黄色视频的| 黑人操中国人逼视频| 亚洲精品一区av在线观看| 免费女性裸体啪啪无遮挡网站| 色综合欧美亚洲国产小说| 变态另类成人亚洲欧美熟女 | 麻豆久久精品国产亚洲av | 人人澡人人妻人| 丰满饥渴人妻一区二区三| 亚洲久久久国产精品| 久久性视频一级片| 国产激情欧美一区二区| 亚洲伊人色综图| 国产亚洲精品久久久久5区| 可以免费在线观看a视频的电影网站| 两性午夜刺激爽爽歪歪视频在线观看 | 中文字幕精品免费在线观看视频| 国产xxxxx性猛交| 欧美日韩乱码在线| 多毛熟女@视频| 手机成人av网站| 欧美老熟妇乱子伦牲交| 久久精品成人免费网站| 久久影院123| 激情视频va一区二区三区| 免费高清视频大片| 黄色片一级片一级黄色片| 午夜日韩欧美国产| 亚洲成人精品中文字幕电影 | 在线观看日韩欧美| 啦啦啦在线免费观看视频4| 亚洲专区中文字幕在线| videosex国产| 亚洲一区高清亚洲精品| 欧美 亚洲 国产 日韩一| 一级a爱片免费观看的视频| 满18在线观看网站| 国内久久婷婷六月综合欲色啪| 夜夜看夜夜爽夜夜摸 | 黄片大片在线免费观看| 50天的宝宝边吃奶边哭怎么回事| 女人被躁到高潮嗷嗷叫费观| 久久久久久大精品| 午夜老司机福利片| 咕卡用的链子| 丁香欧美五月| 麻豆av在线久日| 最新美女视频免费是黄的| 色哟哟哟哟哟哟| 国产一区二区激情短视频| 在线观看免费午夜福利视频| 黑人操中国人逼视频| 一级a爱片免费观看的视频| 视频在线观看一区二区三区| 成人18禁高潮啪啪吃奶动态图| 亚洲国产欧美日韩在线播放| 淫秽高清视频在线观看| 少妇的丰满在线观看| 麻豆成人av在线观看| 日韩三级视频一区二区三区| 国产成人av教育| 色老头精品视频在线观看| 久久精品国产清高在天天线| 久久九九热精品免费| av国产精品久久久久影院| 久久久久亚洲av毛片大全| 电影成人av| 一级毛片高清免费大全| 男人舔女人下体高潮全视频| 亚洲色图 男人天堂 中文字幕| 日韩国内少妇激情av| 一进一出好大好爽视频| 精品一区二区三区四区五区乱码| 欧美日韩瑟瑟在线播放| 日韩欧美三级三区| 国产精品秋霞免费鲁丝片| 日韩精品中文字幕看吧| 国产不卡一卡二| 极品教师在线免费播放| 久久精品91无色码中文字幕| 人人妻人人添人人爽欧美一区卜| 亚洲伊人色综图| 黄色视频,在线免费观看| 亚洲欧洲精品一区二区精品久久久| 自线自在国产av| cao死你这个sao货| 夜夜躁狠狠躁天天躁| 欧美一区二区精品小视频在线| 成人亚洲精品一区在线观看| 老司机午夜福利在线观看视频| 女性被躁到高潮视频| 老汉色av国产亚洲站长工具| 久99久视频精品免费| 久久久久国产精品人妻aⅴ院| 12—13女人毛片做爰片一| 亚洲,欧美精品.| 在线观看免费日韩欧美大片| 亚洲国产中文字幕在线视频| 久久久久久人人人人人| 俄罗斯特黄特色一大片| 丰满人妻熟妇乱又伦精品不卡| 麻豆久久精品国产亚洲av | 国产成人av激情在线播放| 纯流量卡能插随身wifi吗| 久久人妻熟女aⅴ| 正在播放国产对白刺激| 午夜精品国产一区二区电影| 最近最新中文字幕大全免费视频| 50天的宝宝边吃奶边哭怎么回事| 亚洲欧美日韩无卡精品| 成人手机av| 免费在线观看视频国产中文字幕亚洲| 黄片小视频在线播放| 97超级碰碰碰精品色视频在线观看| 亚洲精品中文字幕在线视频| 天堂影院成人在线观看| 成年女人毛片免费观看观看9| 亚洲精品国产区一区二| 国产精品自产拍在线观看55亚洲| 国产日韩一区二区三区精品不卡| 在线视频色国产色| 后天国语完整版免费观看| 欧美精品啪啪一区二区三区| 长腿黑丝高跟| 成熟少妇高潮喷水视频| 欧美亚洲日本最大视频资源| 亚洲欧美一区二区三区黑人| 国产精品电影一区二区三区| 嫩草影视91久久| 又大又爽又粗| av电影中文网址| 99在线人妻在线中文字幕| 欧美中文综合在线视频| 午夜影院日韩av| 亚洲人成电影免费在线| 亚洲三区欧美一区| 久久久久久久久久久久大奶| 91在线观看av| 88av欧美| 天堂影院成人在线观看| 啦啦啦免费观看视频1| 国产欧美日韩综合在线一区二区| 国产亚洲精品久久久久久毛片| 最近最新中文字幕大全免费视频| 一区二区三区国产精品乱码| 丰满的人妻完整版| 99久久精品国产亚洲精品| 亚洲一区二区三区色噜噜 | 成年人黄色毛片网站| 欧美日韩av久久| 国产亚洲欧美在线一区二区| 国产人伦9x9x在线观看| 国产精品久久久av美女十八| 日本一区二区免费在线视频| 黑人巨大精品欧美一区二区mp4| 国内毛片毛片毛片毛片毛片| 怎么达到女性高潮| 国产精品久久久人人做人人爽| 十分钟在线观看高清视频www| 国产深夜福利视频在线观看| 黄色片一级片一级黄色片| 中文字幕精品免费在线观看视频| 久久亚洲精品不卡| 亚洲成国产人片在线观看| 久热爱精品视频在线9| 长腿黑丝高跟| 成人三级做爰电影| 动漫黄色视频在线观看| 免费搜索国产男女视频| 午夜免费鲁丝| 久久人人精品亚洲av| 国产精品偷伦视频观看了| 看黄色毛片网站| 99热只有精品国产| 日本vs欧美在线观看视频| 精品福利观看| 日韩视频一区二区在线观看| 一级片'在线观看视频| 国产一区二区三区视频了| 香蕉丝袜av| 最新美女视频免费是黄的| 久久精品aⅴ一区二区三区四区| 中文欧美无线码| 咕卡用的链子| 免费高清视频大片| 亚洲第一欧美日韩一区二区三区| 日本黄色视频三级网站网址| 少妇裸体淫交视频免费看高清 | 一个人观看的视频www高清免费观看 | 国产国语露脸激情在线看| 色老头精品视频在线观看| 少妇裸体淫交视频免费看高清 | 欧美色视频一区免费| 一级黄色大片毛片| cao死你这个sao货| 超碰成人久久| av在线天堂中文字幕 | 麻豆久久精品国产亚洲av | 成人特级黄色片久久久久久久| 最近最新中文字幕大全免费视频| 在线观看舔阴道视频| videosex国产| 国产激情欧美一区二区| 一二三四社区在线视频社区8| 久久午夜综合久久蜜桃| 国产蜜桃级精品一区二区三区| 国产区一区二久久| 97碰自拍视频| 女警被强在线播放| 免费在线观看影片大全网站| 中文欧美无线码| 不卡一级毛片| 欧美色视频一区免费| 热re99久久精品国产66热6| 久久精品亚洲av国产电影网| 国产伦一二天堂av在线观看| 黄频高清免费视频| 免费在线观看视频国产中文字幕亚洲| 91av网站免费观看| 超碰成人久久| cao死你这个sao货| 老熟妇乱子伦视频在线观看| 亚洲国产欧美一区二区综合| 一级黄色大片毛片| 午夜免费激情av| 精品一区二区三区视频在线观看免费 | 亚洲精华国产精华精| 国产精品一区二区精品视频观看| 国产精品 国内视频| 日日爽夜夜爽网站| 成年人免费黄色播放视频| 欧美久久黑人一区二区| 国产精品久久久人人做人人爽| 欧美另类亚洲清纯唯美| 巨乳人妻的诱惑在线观看| 久久久久久亚洲精品国产蜜桃av| 午夜免费成人在线视频| 无人区码免费观看不卡| 国产高清视频在线播放一区| 制服诱惑二区| 亚洲精品国产区一区二| 青草久久国产| 成人特级黄色片久久久久久久| 日韩高清综合在线| 一级黄色大片毛片| 久久精品91无色码中文字幕| 黄色成人免费大全| 久久精品91无色码中文字幕| 亚洲专区国产一区二区| 超色免费av| 激情视频va一区二区三区| 中亚洲国语对白在线视频| 777久久人妻少妇嫩草av网站| 久久天躁狠狠躁夜夜2o2o| 少妇裸体淫交视频免费看高清 | 97超级碰碰碰精品色视频在线观看| 新久久久久国产一级毛片| av片东京热男人的天堂| 村上凉子中文字幕在线| 久久精品国产99精品国产亚洲性色 | av福利片在线| 国产精品综合久久久久久久免费 | 99国产精品免费福利视频| 亚洲欧洲精品一区二区精品久久久| 国产精品久久久人人做人人爽| 日韩欧美在线二视频| 少妇被粗大的猛进出69影院| 国产高清videossex| 亚洲av美国av| 人人妻人人添人人爽欧美一区卜| 老汉色av国产亚洲站长工具| 可以在线观看毛片的网站| videosex国产| 日韩大尺度精品在线看网址 | ponron亚洲| 国产精品 国内视频| 国产成+人综合+亚洲专区| 精品免费久久久久久久清纯| a级毛片黄视频| 视频在线观看一区二区三区| 久久香蕉精品热| 波多野结衣一区麻豆| 伊人久久大香线蕉亚洲五| 操美女的视频在线观看| 在线观看一区二区三区激情| 欧美另类亚洲清纯唯美| 国产精品久久久人人做人人爽| 在线播放国产精品三级|