• <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
    (口歐)!鷹笛
    3wmmmm亚洲av在线观看| 夜夜爽天天搞| 日日夜夜操网爽| 成人国产麻豆网| 欧美又色又爽又黄视频| 中文在线观看免费www的网站| 香蕉av资源在线| 亚洲午夜理论影院| 日本与韩国留学比较| 少妇丰满av| 国产男靠女视频免费网站| 男女啪啪激烈高潮av片| 我的女老师完整版在线观看| 久久精品夜夜夜夜夜久久蜜豆| 国产精品国产三级国产av玫瑰| 神马国产精品三级电影在线观看| 十八禁国产超污无遮挡网站| 国产精品永久免费网站| 老司机午夜福利在线观看视频| 99在线视频只有这里精品首页| 久久久久久久久大av| 亚洲av美国av| 精品国产三级普通话版| 免费av毛片视频| 国内精品久久久久精免费| 久久久久久久久久黄片| 一个人免费在线观看电影| 国产精品野战在线观看| or卡值多少钱| 久久人人爽人人爽人人片va| www.色视频.com| 欧美性感艳星| 欧美精品国产亚洲| av视频在线观看入口| 亚洲精品粉嫩美女一区| 欧美成人a在线观看| 三级毛片av免费| 欧美另类亚洲清纯唯美| 又爽又黄a免费视频| 91麻豆av在线| 午夜爱爱视频在线播放| 国产人妻一区二区三区在| 在现免费观看毛片| 国产成人影院久久av| 日本精品一区二区三区蜜桃| 国产伦精品一区二区三区视频9| 嫩草影视91久久| 深夜a级毛片| av国产免费在线观看| 亚洲国产精品合色在线| 成年版毛片免费区| 内射极品少妇av片p| 日本一本二区三区精品| 色尼玛亚洲综合影院| 久久天躁狠狠躁夜夜2o2o| 国产欧美日韩一区二区精品| 成年人黄色毛片网站| 国产成人aa在线观看| 亚洲欧美精品综合久久99| 能在线免费观看的黄片| 亚洲男人的天堂狠狠| 国产淫片久久久久久久久| 悠悠久久av| 日韩一本色道免费dvd| 日韩中字成人| 亚洲一级一片aⅴ在线观看| 十八禁网站免费在线| 性插视频无遮挡在线免费观看| 精品人妻熟女av久视频| 麻豆av噜噜一区二区三区| 色视频www国产| 九九爱精品视频在线观看| 久久久久久久久中文| 在线观看av片永久免费下载| 国产亚洲91精品色在线| 色精品久久人妻99蜜桃| 日本在线视频免费播放| 国产精品亚洲美女久久久| 女同久久另类99精品国产91| 少妇被粗大猛烈的视频| 91久久精品国产一区二区三区| 亚洲五月天丁香| 成人性生交大片免费视频hd| 亚洲久久久久久中文字幕| 在线免费观看不下载黄p国产 | 狂野欧美激情性xxxx在线观看| 国产精品人妻久久久影院| 麻豆国产97在线/欧美| 黄色视频,在线免费观看| 亚洲av一区综合| 99精品久久久久人妻精品| 成人国产麻豆网| 日本与韩国留学比较| 亚洲国产欧洲综合997久久,| 日本欧美国产在线视频| 国产大屁股一区二区在线视频| 成人综合一区亚洲| 一个人看视频在线观看www免费| 国产乱人视频| 日日撸夜夜添| 久久久精品欧美日韩精品| 久久九九热精品免费| 亚洲成人久久爱视频| 欧美极品一区二区三区四区| 91在线精品国自产拍蜜月| 欧美日韩乱码在线| 高清日韩中文字幕在线| 嫩草影院精品99| 亚洲av免费高清在线观看| 啪啪无遮挡十八禁网站| 人人妻人人澡欧美一区二区| 久久精品国产自在天天线| 88av欧美| 亚洲av二区三区四区| 亚洲欧美日韩东京热| 91狼人影院| 国产精品一区二区性色av| 三级国产精品欧美在线观看| 精品午夜福利在线看| 国产淫片久久久久久久久| 午夜免费男女啪啪视频观看 | 久久精品夜夜夜夜夜久久蜜豆| www.www免费av| 观看免费一级毛片| 久久久午夜欧美精品| 国内毛片毛片毛片毛片毛片| 一本精品99久久精品77| 男人舔奶头视频| 少妇裸体淫交视频免费看高清| 免费看av在线观看网站| 一级黄色大片毛片| 九色成人免费人妻av| 精品无人区乱码1区二区| 内射极品少妇av片p| 久久这里只有精品中国| 在线免费十八禁| 美女被艹到高潮喷水动态| 自拍偷自拍亚洲精品老妇| 日日干狠狠操夜夜爽| 欧美日韩亚洲国产一区二区在线观看| 精品一区二区三区人妻视频| av在线天堂中文字幕| 老司机午夜福利在线观看视频| .国产精品久久| 长腿黑丝高跟| 国产精品伦人一区二区| 国产一区二区三区视频了| 亚洲av日韩精品久久久久久密| 免费看av在线观看网站| 不卡一级毛片| 国产高清有码在线观看视频| 中文字幕人妻熟人妻熟丝袜美| 国产免费一级a男人的天堂| 99九九线精品视频在线观看视频| 一本精品99久久精品77| 国产三级中文精品| 亚洲精品久久国产高清桃花| 欧美高清成人免费视频www| 五月伊人婷婷丁香| 国产美女午夜福利| 婷婷精品国产亚洲av| 国产精品国产三级国产av玫瑰| 亚洲人与动物交配视频| 亚洲av五月六月丁香网| 亚洲熟妇中文字幕五十中出| 国产精品久久久久久久久免| 精品不卡国产一区二区三区| 亚洲无线在线观看| 尤物成人国产欧美一区二区三区| 自拍偷自拍亚洲精品老妇| 亚洲五月天丁香| 啦啦啦韩国在线观看视频| 99国产极品粉嫩在线观看| 国产黄a三级三级三级人| 国产精品女同一区二区软件 | 中文资源天堂在线| 午夜福利在线观看吧| 69av精品久久久久久| a在线观看视频网站| 两人在一起打扑克的视频| 亚洲av中文字字幕乱码综合| 亚洲成人久久性| 我的女老师完整版在线观看| 天美传媒精品一区二区| www.www免费av| 三级毛片av免费| 国产亚洲av嫩草精品影院| 午夜亚洲福利在线播放| 日韩高清综合在线| 在线免费观看不下载黄p国产 | 国产午夜精品久久久久久一区二区三区 | 啦啦啦韩国在线观看视频| 色综合站精品国产| 女同久久另类99精品国产91| 在线天堂最新版资源| 日韩人妻高清精品专区| 午夜影院日韩av| 精品一区二区三区人妻视频| 99精品在免费线老司机午夜| 日韩中文字幕欧美一区二区| 精品久久久久久成人av| 久久草成人影院| 国产成人a区在线观看| 欧美不卡视频在线免费观看| 我的老师免费观看完整版| 亚洲美女视频黄频| 一进一出好大好爽视频| 综合色av麻豆| 舔av片在线| 一本一本综合久久| 日韩欧美在线乱码| 亚洲av美国av| 国产精品伦人一区二区| 中文字幕av成人在线电影| 天天躁日日操中文字幕| 国产成人一区二区在线| 亚洲国产欧美人成| 精品久久久久久久久久免费视频| 国产一区二区在线观看日韩| 亚洲成a人片在线一区二区| 精品不卡国产一区二区三区| 亚洲av熟女| 啦啦啦观看免费观看视频高清| 无遮挡黄片免费观看| 男插女下体视频免费在线播放| av在线亚洲专区| 精品乱码久久久久久99久播| 女的被弄到高潮叫床怎么办 | 亚洲精品影视一区二区三区av| 可以在线观看毛片的网站| 一区福利在线观看| 97超级碰碰碰精品色视频在线观看| 在线a可以看的网站| 国产精品日韩av在线免费观看| 成人亚洲精品av一区二区| 99久久久亚洲精品蜜臀av| 欧美成人性av电影在线观看| 亚洲成av人片在线播放无| 91麻豆精品激情在线观看国产| 一进一出抽搐动态| 亚洲中文字幕日韩| 亚洲五月天丁香| 韩国av一区二区三区四区| 国产91精品成人一区二区三区| 美女高潮的动态| 国产午夜福利久久久久久| 最近中文字幕高清免费大全6 | 丰满的人妻完整版| 欧美另类亚洲清纯唯美| 精品人妻视频免费看| 1024手机看黄色片| 亚洲av二区三区四区| 18禁在线播放成人免费| 国产在线男女| 亚洲熟妇中文字幕五十中出| 国产又黄又爽又无遮挡在线| 美女免费视频网站| xxxwww97欧美| 亚洲成人久久性| 亚洲精华国产精华精| 91在线观看av| 久久精品综合一区二区三区| 麻豆av噜噜一区二区三区| 波野结衣二区三区在线| 一进一出好大好爽视频| 嫩草影院入口| 国产精品久久电影中文字幕| 99久久九九国产精品国产免费| 免费不卡的大黄色大毛片视频在线观看 | 国产精品98久久久久久宅男小说| 日韩欧美一区二区三区在线观看| 少妇裸体淫交视频免费看高清| 亚洲国产精品久久男人天堂| 欧美精品啪啪一区二区三区| 亚洲美女视频黄频| 春色校园在线视频观看| av在线老鸭窝| 大又大粗又爽又黄少妇毛片口| 国产精品综合久久久久久久免费| 国产主播在线观看一区二区| 性插视频无遮挡在线免费观看| 国产真实乱freesex| 99精品在免费线老司机午夜| 亚洲国产日韩欧美精品在线观看| 国产人妻一区二区三区在| 一边摸一边抽搐一进一小说| 在现免费观看毛片| 天堂网av新在线| 99久久久亚洲精品蜜臀av| 美女黄网站色视频| 天美传媒精品一区二区| 精品99又大又爽又粗少妇毛片 | 男女视频在线观看网站免费| 午夜精品在线福利| 亚洲,欧美,日韩| 少妇的逼好多水| 亚洲人成网站在线播| 亚洲成人久久性| 校园人妻丝袜中文字幕| 亚洲avbb在线观看| 熟女电影av网| 在线国产一区二区在线| 97碰自拍视频| 国产精品久久久久久精品电影| 丰满乱子伦码专区| 国产美女午夜福利| 欧美一级a爱片免费观看看| 日本成人三级电影网站| 在线观看一区二区三区| 免费高清视频大片| 欧美激情久久久久久爽电影| 日韩精品青青久久久久久| 两个人视频免费观看高清| 国产综合懂色| 国产三级中文精品| 夜夜爽天天搞| 成人综合一区亚洲| 老熟妇仑乱视频hdxx| 亚洲成人中文字幕在线播放| 国产精品一区二区三区四区免费观看 | 午夜福利成人在线免费观看| 亚洲中文字幕一区二区三区有码在线看| 麻豆国产97在线/欧美| 内地一区二区视频在线| 日本免费一区二区三区高清不卡| 亚洲av第一区精品v没综合| av在线老鸭窝| 婷婷精品国产亚洲av在线| 中亚洲国语对白在线视频| 又爽又黄无遮挡网站| 欧美不卡视频在线免费观看| 一区福利在线观看| 久久国内精品自在自线图片| 特级一级黄色大片| 国产成人福利小说| 国产精品久久电影中文字幕| 久久欧美精品欧美久久欧美| 亚洲性夜色夜夜综合| 日本黄色片子视频| 国内少妇人妻偷人精品xxx网站| 一边摸一边抽搐一进一小说| 国产探花极品一区二区| 久久人人爽人人爽人人片va| 国产久久久一区二区三区| 身体一侧抽搐| 淫秽高清视频在线观看| 观看免费一级毛片| 亚洲精品在线观看二区| 免费人成视频x8x8入口观看| 国产成人a区在线观看| 日日啪夜夜撸| 少妇裸体淫交视频免费看高清| 久久久久精品国产欧美久久久| 久久人人精品亚洲av| 久久久色成人| 老司机深夜福利视频在线观看| 在线免费观看的www视频| 能在线免费观看的黄片| 欧美日韩瑟瑟在线播放| 久久精品影院6| 精品日产1卡2卡| 白带黄色成豆腐渣| 美女大奶头视频| 免费搜索国产男女视频| 国产亚洲精品综合一区在线观看| 九九久久精品国产亚洲av麻豆| 男女视频在线观看网站免费| 精品一区二区三区av网在线观看| 中文字幕av成人在线电影| 色播亚洲综合网| 日韩在线高清观看一区二区三区 | 性插视频无遮挡在线免费观看| 成人一区二区视频在线观看| 午夜福利欧美成人| av女优亚洲男人天堂| 九九在线视频观看精品| 日本爱情动作片www.在线观看 | 日日摸夜夜添夜夜添小说| 97人妻精品一区二区三区麻豆| 春色校园在线视频观看| 99久久九九国产精品国产免费| 久久久精品欧美日韩精品| 中文在线观看免费www的网站| 不卡一级毛片| 99国产精品一区二区蜜桃av| 免费看日本二区| 免费看光身美女| 成年女人永久免费观看视频| x7x7x7水蜜桃| 国内精品一区二区在线观看| 一a级毛片在线观看| 性插视频无遮挡在线免费观看| 国产精品不卡视频一区二区| 久久精品夜夜夜夜夜久久蜜豆| 国产熟女欧美一区二区| 18禁黄网站禁片免费观看直播| 亚洲av第一区精品v没综合| 特大巨黑吊av在线直播| 国产精品亚洲一级av第二区| 国产精品无大码| 欧美xxxx黑人xx丫x性爽| 亚洲欧美日韩高清在线视频| 国产一区二区三区视频了| 午夜亚洲福利在线播放| 国产亚洲欧美98| 一区二区三区激情视频| 日本免费a在线| 午夜福利成人在线免费观看| 国产午夜福利久久久久久| 国产精品精品国产色婷婷| 成年免费大片在线观看| 亚洲国产高清在线一区二区三| 精品久久久久久久末码| 中文亚洲av片在线观看爽| 日本一二三区视频观看| 22中文网久久字幕| 麻豆av噜噜一区二区三区| 亚洲精品成人久久久久久| 欧美另类亚洲清纯唯美| 国产精品免费一区二区三区在线| 国内精品久久久久精免费| 久久精品综合一区二区三区| 熟妇人妻久久中文字幕3abv| 色精品久久人妻99蜜桃| 国产精品一区二区性色av| 欧美xxxx性猛交bbbb| 两个人的视频大全免费| 三级毛片av免费| 国产精品无大码| 亚洲avbb在线观看| 欧美黑人欧美精品刺激| 国产精品人妻久久久久久| 我要搜黄色片| 亚洲专区国产一区二区| 国产欧美日韩一区二区精品| 22中文网久久字幕| 国产精品伦人一区二区| 美女免费视频网站| 97超级碰碰碰精品色视频在线观看| ponron亚洲| 亚洲精品亚洲一区二区| 久久九九热精品免费| 夜夜夜夜夜久久久久| 成人av一区二区三区在线看| 欧美成人a在线观看| av在线亚洲专区| 九色成人免费人妻av| 精品久久久久久久久久免费视频| 精品一区二区免费观看| 91精品国产九色| 亚洲精品一区av在线观看| 国产欧美日韩一区二区精品| 久久精品91蜜桃| 最近最新免费中文字幕在线| 免费在线观看日本一区| 美女黄网站色视频| 亚洲国产精品久久男人天堂| 日本色播在线视频| 亚洲一区高清亚洲精品| 久久久午夜欧美精品| 中文字幕精品亚洲无线码一区| 精品不卡国产一区二区三区| 亚洲aⅴ乱码一区二区在线播放| netflix在线观看网站| 精品久久国产蜜桃| 精品国内亚洲2022精品成人| 亚洲精品色激情综合| 亚洲av电影不卡..在线观看| eeuss影院久久| 中文字幕免费在线视频6| 久久久成人免费电影| 观看免费一级毛片| 亚洲熟妇中文字幕五十中出| 大又大粗又爽又黄少妇毛片口| 国产一区二区亚洲精品在线观看| 久久久久久久久中文| 成人特级av手机在线观看| 真人做人爱边吃奶动态| 国产精品三级大全| 99热精品在线国产| 五月玫瑰六月丁香| 日本一本二区三区精品| 香蕉av资源在线| 精品人妻一区二区三区麻豆 | 久久久色成人| 又黄又爽又免费观看的视频| 欧美在线一区亚洲| 国产伦一二天堂av在线观看| 舔av片在线| 亚洲精华国产精华液的使用体验 | 精品人妻偷拍中文字幕| 国产视频内射| 亚洲精品日韩av片在线观看| 久久精品国产99精品国产亚洲性色| 男女啪啪激烈高潮av片| 国产视频内射| 欧美日本亚洲视频在线播放| 欧美中文日本在线观看视频| 99国产极品粉嫩在线观看| 国产视频内射| 久久国产乱子免费精品| 一卡2卡三卡四卡精品乱码亚洲| 男女做爰动态图高潮gif福利片| 国产又黄又爽又无遮挡在线| 天堂√8在线中文| 国产av麻豆久久久久久久| 亚洲无线在线观看| 欧美日韩国产亚洲二区| 一个人免费在线观看电影| 伦理电影大哥的女人| 日韩高清综合在线| 欧美最新免费一区二区三区| 日本欧美国产在线视频| 又粗又爽又猛毛片免费看| 桃色一区二区三区在线观看| 婷婷色综合大香蕉| 高清在线国产一区| 在线免费观看的www视频| 欧美最黄视频在线播放免费| 国产精品三级大全| 国产伦在线观看视频一区| 一个人观看的视频www高清免费观看| 熟女电影av网| 99在线视频只有这里精品首页| 亚洲专区中文字幕在线| 在线免费观看不下载黄p国产 | 一卡2卡三卡四卡精品乱码亚洲| 久久久久久伊人网av| 国产综合懂色| 国产精品人妻久久久影院| 久久久久精品国产欧美久久久| 国产高清视频在线观看网站| 成人一区二区视频在线观看| 欧美一区二区国产精品久久精品| 精品久久久久久久人妻蜜臀av| 国产av不卡久久| 中文字幕免费在线视频6| 久久99热这里只有精品18| 国产大屁股一区二区在线视频| 国产成年人精品一区二区| 人妻丰满熟妇av一区二区三区| 国产一区二区激情短视频| 午夜免费男女啪啪视频观看 | 日本一二三区视频观看| 身体一侧抽搐| 国产麻豆成人av免费视频| 日韩精品青青久久久久久| 久久99热6这里只有精品| 麻豆久久精品国产亚洲av| 精品久久久久久久末码| 国产视频一区二区在线看| 国产亚洲精品久久久com| 国内精品美女久久久久久| 狂野欧美白嫩少妇大欣赏| 午夜影院日韩av| 女的被弄到高潮叫床怎么办 | 日本在线视频免费播放| 99热只有精品国产| 国产精品野战在线观看| a级毛片免费高清观看在线播放| 亚洲成av人片在线播放无| 最近中文字幕高清免费大全6 | 精品乱码久久久久久99久播| 日日夜夜操网爽| 欧美性猛交╳xxx乱大交人| 欧美激情国产日韩精品一区| 夜夜爽天天搞| 999久久久精品免费观看国产| 乱系列少妇在线播放| 色在线成人网| 日本三级黄在线观看| 久久久久久久久久黄片| 成人欧美大片| 少妇裸体淫交视频免费看高清| 国产乱人伦免费视频| 国产精华一区二区三区| 久久热精品热| 蜜桃亚洲精品一区二区三区| 一进一出好大好爽视频| 啦啦啦韩国在线观看视频| 一级毛片久久久久久久久女| 免费无遮挡裸体视频| 美女cb高潮喷水在线观看| 亚洲av中文av极速乱 | 精品一区二区免费观看| 国产在线精品亚洲第一网站| 国国产精品蜜臀av免费| 久久精品国产亚洲av涩爱 | 黄片wwwwww| 亚洲人成网站高清观看| 中文字幕熟女人妻在线| 网址你懂的国产日韩在线| 久久久色成人| 午夜福利18| 老师上课跳d突然被开到最大视频| 有码 亚洲区| 亚洲最大成人中文| 高清毛片免费观看视频网站| 亚洲 国产 在线| 婷婷丁香在线五月| 88av欧美| 尾随美女入室| 97热精品久久久久久| ponron亚洲| 成人二区视频| 香蕉av资源在线| 麻豆一二三区av精品| 日本a在线网址| 九九在线视频观看精品| 午夜爱爱视频在线播放| 国产精品福利在线免费观看| 亚洲不卡免费看| 亚洲熟妇中文字幕五十中出| 欧美日韩国产亚洲二区| 欧美xxxx黑人xx丫x性爽|