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

    Formation of carbon and oxygen rich surface layer on high purity magnesium by atmospheric carbon dioxide plasma

    2023-03-24 13:19:10GyoungJangSinhulYeomJongKeumMinaYoonHarryIIIMeyerYiFengSuJiheonJun
    Journal of Magnesium and Alloys 2023年1期

    Gyoung G.Jang ,Sinhul Yeom ,Jong K.Keum ,Mina Yoon ,Harry III Meyer ,Yi-Feng Su ,Jiheon Jun,?

    a Manufacturing Science Division,Oak Ridge National Laboratory (ORNL),One Bethel Valley Road,PO Box 2008,MS6054,Oak Ridge,TN 37831,United States

    b Material Science and Technologies Division,Oak Ridge National Laboratory (ORNL),One Bethel Valley Road,Oak Ridge,TN 37831,United States

    c Neutron Scattering Division &Center for Nanophase Materials Sciences,ORNL,One Bethel Valley Road,Oak Ridge,TN 37831,United States

    d Chemical Sciences Division,ORNL,One Bethel Valley Road,Oak Ridge,TN 37831,United States

    Abstract Carbon and oxygen-rich corrosion barrier layer formed on Mg by a simple and scalable CO2 atmospheric plasma (CO2-AP) process.The reactive CO2-AP interacts with the Mg surface and forms a unique layered structure with the top MgCO3/MgO-intermixed particulates pillars and the bottom dense layer.The surface features were simultaneously formed on the nano-/micro-structured MgO layer by carbonate molecules,plasma-active CO2 molecules,and/or other volatile organic compounds on the nano-/micro-structured MgO particle layer.The resulting surfaces after CO2-AP were either hydrophobic or hydrophilic and exhibited lower anodic current or high resistance for Mg corrosion.For the hydrophobic surfaces of CO2-AP treated Mg,molecular dynamic simulations were performed to understand the origin of hydrophobicity and identifie that the amorphous carbon layers formed on the Mg surface are the source.The environmentally benign abundant-gas-based process enables the cost reduction associated with waste treatment,generation of by-product,and supply of raw material.

    Keywords: Atmospheric plasma;Magnesium;Corrosion protection;Surface modification

    1.Introduction

    Poor corrosion resistance of Mg is one of the major bottlenecks to limit the wide application of Mg alloys in lightweight vehicle and aircraft components,bioengineering implants,and battery and fuel cell electrodes [1–6].In humid air and water,the surface of Mg naturally forms Mg(OH)2and MgO layers[7] which are often porous and unprotective,especially in the presence of NaCl.Various alloying strategies have demonstrated a moderate increase of Mg corrosion resistance;however,protective coatings are still considered to be the most cost-effective way for improving the corrosion resistance of Mg and Mg alloys.

    Conventional coating techniques such as chemical and electrochemical conversion coatings [8],anodizing [9,10],electroplating [11],thermal powder spray coatings [12] or polymer coatings have been developed to fabricate corrosionresistant coatings and thus provide protection for practical applications.However,reliable Mg protection typically requires multilayer coatings via a series of processes [2,13].Also,these wet-chemical-based coating approaches must follow the standard of US Environmental Protection Agency regulations regarding volatile organic compound generation.These conditions raise the issues of the complexity of processes,high cost,and an inadequate degree of corrosion protection.

    Recently,superhydrophobic coatings have been considered as a promising technology to significantl improve the corrosion resistance of metallic substrates because this method can effectively reduce the contact area from aqueous corrosive medium [14–16].As is well-known,the superhydrophobic surfaces on metallic substrates are usually fabricated by a two-step route: (1) creating nano/micro surface structures on a substrate with particle deposition,etching,and conversion;(2) modifying the surface with low-surface-energy compounds such as F and C.Still,the technical issue is increasing adhesion between nano/microstructures and metal surfaces to achieve durable superhydrophobicity.Hydrothermal treatment and chemical conversion in a salt solution have been studied to growin-situMgO/Mg(OH)2-corroded product on the Mg surface for the strong nano/microstructure as a single entity.The structured surfaces are then modifie with low-surface-energy compounds such as steric acid and fluorosilan [15–18].

    A new coating approach recently reported that a corrosionresistant MgCO3layer on the surface of the micro Mg structure was formed via direct reaction of an already-existing surface oxide fil with excited CO2under electron-beam irradiation [19].Magnesium carbonate is a stable compound with a very low solubility of 20 mg per 100 mL in water;it can stay inert in water for a long time and grow in a dense and compact layer on the Mg metal surface [19].Additionally,research on film formed on Mg substrates under atmospheric conditions containing CO2shows that layers containing carbonates/hydroxy carbonates on the surface of Mg alloys can suppress corrosion by blocking the anodic and/or cathodic sites [20].

    Inspired by the reported MgCO3layer formation under CO2-driven,electron-beam/plasma irradiation and improved corrosion resistance [19],the objective of this study is to demonstrate the facile formation of MgCO3for an anticorrosion layer on Mg via an atmospheric room temperature CO2chemical plasma technique (Fig.A1).The atmospheric plasma(AP)is generated by a high-frequency electromagnetic fiel under high voltage.As compressed gas fl ws across the field the atoms become ionized with free electrons and atomic and molecular radicals in excited levels,producing a “cold”plasma [21].This highly reactive environment is maintained close to the room temperature.

    Depending on the principle of generation,these APs can be employed directly in the production line because no need exists for high-cost vacuum chambers.The cold plasma has a high chemical energy but a low thermal energy,so it should not damage temperature-sensitive materials [22].Currently,the AP technique has been successfully industrialized in largescale deployments such as an automated in-line continuous system for rapid,low-cost,open-air surface cleaning,activating,and cutting.By simply using CO2instead of O,N,and Ar as the source feedstock gas,and ideally with CO2as a cover gas,our approach provides improved corrosion protection for Mg.

    In the present study,we report a new corrosion mitigation method using the CO2-AP treatment.Among the CO2plasma-treated Mg samples,some surfaces showed hydrophobicity which remarkably reduced Mg corrosion.To mechanistically demonstrate the hydrophobic Mg surfaces,molecular dynamics (MD) simulations were also performed.

    2.Experimental

    2.1.Materials

    The cast commercial purity Mg,approximately 99.98%Mg[23],was obtained from Magnesium Elektron North America,Madison,IL,USA.Test samples,25 mm×25 mm in size and ~1 to 1.5 mm thick,were obtained for all samples by electric discharge machining.The test samples were wet ground with SiC paper 120 to 600 grit finis for rough and smooth references,cleaned with deionized water,and dried with an air stream.The samples were stored in a desiccator before AP treatment.

    2.2.Preparation of the CO2-AP–treated Mg substrate

    The CO2-AP treatment was carried out in a customdesigned CO2-enriched cover gas chamber to increase the MgCO3formation.The cold-plasma treatment was performed using a PlasmaFlux system manufactured by Atmospheric Plasma Solutions [22].The process was performed in a working chamber with a 3D robotic surface equipped with an exhaust and filtratio system.The plasma generated a highvelocity fl w of chemically active CO2radicals that exited the nozzle tip as a plasma plume,as shown in Fig.A1 in the appendix.The robotic stage and plasma pen were controlled by three stepper motors and three rack and pinion drives that have the ability to move in the x (forward),y (side),and z(height) directions.Various plasma operation conditions were tested to uniformly treat on surface.For example,the nozzlegenerating CO2-AP by 100% CO2gas passed through the sample surface with 0.38 cm height,762 cm/min and 0.06 cm spacing at 70 SLPM of fl w rate.The swing number,which is the overlapping of plasma treatments,increased to intensify the CO2-AP effect for uniform and multiple layers.

    2.3.Characterization of materials

    Surface treated samples were analyzed by surface scanning electron microscopy (SEM) with energy dispersive Xray analysis (EDS),X-ray diffraction (XRD),X-ray photoelectron spectroscopy (XPS),and cross-section scanning transmission electron microscopy (STEM) in dark fiel (DF)and high angle annular dark fiel (HAADF) modes.Water contact angles (WCA) were measured using an optical tensiometer (OneAttension,Biolin Scientific) The XPS analysis used a monochromated,micro-focusing,Al KαX-ray source(1486.6 eV) with an ~400 μm X-ray spot size for maximum signal and to obtain an average surface composition over the largest possible area [24].Depth profilin analyses were conducted with a Thermo Scientifi EX06 argon ion gun operated at 1000 eV and rastered over a 2 mm×4 mm area.The specimens for STEM analysis were prepared via the focus ion beam (FIB) milling technique by thein-situlift-out method using a Hitachi NB5000 FIB-SEM.A tungsten and/or carbon and tungsten overlayer was deposited to protect the exposed top surface of the fil during ion milling.Further details of the sample preparation and analysis methodology are provided elsewhere [24].

    2.4.Corrosion evaluation

    Potentiodynamic polarization measurements were conducted for the 600-grit finishe and CO2-AP-treated Mg samples as working electrodes,with the exposed area of 0.98 cm2(7×14 mm2) or 0.833 cm2(?=10.3 mm) prepared by insulation tape masking.Reference saturated calomel electrode (SCE) and Pt plate counter electrode with ~5 cm2of surface area were used in 3.5 wt.% NaCl solution.In the beginning,the Mg samples were resting at open-circuit potential(OCP) for 50–60 min.After this open-circuit potential delay,potentiodynamic polarization was performed with a rate of 1 mV·s?1in two different potential windows: for the exposed area of 0.98 cm2,?1.8 to ?1.2 VSCE,and for 0.833 cm2,?2.25 to ?0.95 VSCE.

    Electrochemical impedance spectroscopy (EIS) was also performed to assess the corrosion behavior with increasing immersion time in 3.5 wt.% NaCl solution.OCP was measured for 1 h before the firs EIS measurement for the impedance data with 1 h immersion.The next EIS measurements were also performed for the same samples at longer immersion times.The frequency window was 200 kHz to 7 mHz,and the amplitude was ±10 mV with respect to the last OCP value.The impedance data were analyzed using a commercial computer software(Zview,Scribner)to determine the key resistance values for corrosion evaluation.All electrochemical measurements were conducted at room temperature and open to air.

    H2gas evolution measurements during immersion were also carried out in 3.5 wt.% NaCl solution at ambient temperature for untreated and CO2-AP-treated Mg samples.The edges of all samples were masked with a polymer tape to expose a fi ed area of 7×14 mm2(0.98 cm2) to NaCl solution.Each sample was positioned upright in a glass beaker containing 500 mL of 3.5 wt.% NaCl.The solution was open to air for the duration of the experiment.An inverted funnel was placed over the sample,and a burette fille with 3.5 wt.%NaCl solution was placed over the tip of the inverted funnel.This apparatus collected H2gas evolved during Mg corrosion[25].After exposure,the samples were removed from the solution,rinsed with deionized water and then ethanol,dried using a stream of air,and photographed.

    3.Computational modeling

    3.1.Modeling magnesite (MgCO3) and its interaction with water droplets

    Magnesium carbonate or magnesite is an inorganic salt of trigonal crystal (symmetry group of) with experimental lattice parameters of a=b=4.635,c=15.019,α=β=90°,andγ=120° [26].We employed a forcefiel potential [27,28] that was established for alkaline-earth carbonate and water systems and leads to the lattice parameters a=b=4.611,c=15.069,α=β=90°,andγ=120°.At 300 K,these values are slightly modifie to a=b=4.615,and c=15.165[26].The interaction between C and O atoms of carbonate is described by the 12–6 Lenard-Jones potential,,whereεC–O=0.1 eV andσC-O=2.98.These parameters were determined to eliminate any internal pressure at 300 K by performing canonical MD simulations as implemented in the Large-scale Molecular Massively Parallel Simulator (LAMMPS) program [29].

    We then modeled the magnesite surface as a rectangularlike slab consisting of a 117.3×119.7×13.03size supercell (containing total 20,800 atoms) with large vacuum layers of 100,20,and 20each along the z,x,and y axes,with periodic boundary conditions,where ± (104) crystallographic orientation along the z axis andalong the x and y axes,respectively,were exposed.Then,a water droplet of ~6 nm diameter (consisting of 1000 water molecules) was introduced on the (104) surface and the configuratio was fully optimized until the maximum force became less than 10?8eV/,using the conjugate gradient scheme implemented in LAMMPS [28].After that,we began MD simulations using the optimized configurations

    3.2.Modeling graphite and amorphous carbon and their interaction with water droplets

    We compared the wetting behavior of water molecules on two different carbon surfaces,graphite and amorphous carbon.Our graphite is modeled as 6 graphene layers with a surface area of 25.6×398.52containing 22,842 carbon atoms,where the atomic distance between the nearest carbon atoms is 1.42and the distance between the layers is 3.35.We then used the adaptive intermolecular reactive empirical bond order—Morse (AIREBO-M) force-fiel potential [30] to generate an amorphous carbon structure according to the scheme in the literature [31].The initial structure is a simple cubic cell with 27,900 carbon atoms of the same atomic density of graphite,2.267 g/cm3.We then performed MD simulations to treat the structure with heating,annealing,and cooling cycles[30].The structure was heated to 8,000 K for 10 ps,cooled to 2,000 K for 100 ps,annealed at 4,000 K for 100 ps,and cooled to 300 K for 10 ps,resulting in an amorphous carbon slab with dimensions of 31.0×384.0×20.63.We applied periodic boundary conditions in x and y directions for all the supercell calculations.

    As shown in Fig.1a,the amorphous carbon completely lost the essential features of the initial configuratio (“graphite,”six-layered graphenes),such as the network of sp2hexagonal rings and their arrangement in two-dimensional layers,and it became completely amorphous.It consisted of a randomly distributed mixture of hexagons and polygons interconnected to form a three-dimensional sp3network of carbon atoms.This structural feature is confirme in the radial distribution functions (RDFs) comparing the atomic distance with their distribution.Fig.1b compares the RDFs of the graphite and the amorphous carbon.Because of the well-define structure of graphite,its RDF shows discrete peaks at r~1.4,2.5,2.8,and so on,corresponding to the nth-nearest-neighbor atomic distances.Meanwhile,the RDF of the amorphous carbon shows relatively smooth and broad peaks with the firs two peaks located at r~1.4 and 2.5.

    Fig.1.(a) Side view of the amorphous C slab with magnificatio of the structure in the yellow region.(b) Radial distribution functions of “graphite” (our six-layer graphene system) and amorphous C.

    To describe the interaction between carbon and water system,we employed a force fiel [32] which combines point charge interactions for electrostatic interactions and a Lenard-Jones potential for the chemical interactions.Our parameters of partial charges for Coulomb interactions are O=?0.8476 e and H=+0.4238 e,the Lenard-Jones parameters for O–O areεO-OO=0.00674 eV andσO-O=3.16,for C–O areεC-O=0.00412 eV andσC-O=3.19with a cut-off distance of 10.A water droplet of ~7 nm diameter (consisting of 3,969 water molecules) was relaxed and annealed on the carbon surface at 300 K for 10 ps by canonical MD simulations with the coordinates of the surface fi ed.Next,we performed 0.9 ns of MD simulations for sampling.After that,we collected the configuration of the last 10 snapshots at every 10 ps and analyzed the contact angles between a water droplet on each surface for the further analysis.

    4.Results and discussion

    4.1.Characteristics of the CO2 plasma-treated surface

    CO2-AP treated Mg samples resulted in either hydrophobic or hydrophilic surfaces where some examples are shown in Fig.2a and b.The untreated Mg surface (600 grit SiC finish was hydrophilic as exemplifie in Table 1.Despite our repeated reviews on plasma processing parameters,it was impossible to defin the cause (es) for hydrophobic or hydrophilic treated Mg surfaces.As Mg is highly reactive with any environmental factor (e.g.humidity),the affinit to water on treated Mg samples may have been impacted by the conditions in which the samples were stored during the transportation between the plasma treatment facility in Cary,North Carolina and corrosion testing and microscopic characterization labs in Oak Ridge,Tennessee.For the treated samples with hydrophobic surface,the value of water contact angle increased with increasing the swing number of CO2-AP as presented in Table 1.

    Table 1 Water contact angle measurement on CO2-AP–treated Mg surfaces and corresponding surface morphologies of the treated surfaces.The hydrophobicity increased with the plasma swing number,correlated with surface morphology changes.

    Here,we focus on characterizing a treated hydrophobic Mg specimen with 5-swing of CO2-AP,and for the sample,two SEM images are presented in Fig.2c and d.In plan view,the treated surface exhibited agglomerated sub-micron particles which resulted in a rough profile In cross-sectional view(Fig.2d),duplex surface structure,with ~1 μm thick inner dense layer and outer porous pillars containing sub-micron particles,is observed on Mg substrate.For the observed sample crosssection,there was no apparent defect or gap between plasmadriven duplex layer structures and Mg substrate (see Fig.2d).The surface nanoparticles appeared interconnected each other,which implies that the particles would not easily be detached upon external friction or other stresses.

    4.2.Corrosion mitigation effect of the CO2 plasma-treated surface

    In a simple immersion test in 3.5 wt.% NaCl,a CO2-AP treated Mg sample showed reduced corrosion attack with respect to an untreated Mg sample as observed in Fig.A2 in the appendix.Typical potentiodynamic curves are presented for hydrophobic and hydrophilic CO2-AP treated (with 5-swing)Mg samples as well as untreated Mg samples in Fig.3.Hydrophobic treated surfaces showed notably low current compared to the untreated Mg samples (Fig.3a),which is highly likely due to minimal wetting between the hydrophobic surface and NaCl solution.However,a hydrophilic Mg sample showed cathodic and anodic currents that were slightly larger than the untreated sample (Fig.3b),implying no corrosion mitigation if CO2-AP treated surfaces are hydrophilic.

    Fig.2.SEM images showing a CO2 plasma-driven protective layer that consisted of particle deposition on a dense layer with ~1.5 μm thickness.(a) Photoimage of a CO2-AP–treated superhydrophobic Mg surface repelling water droplets (b) Photo-image of hydrophilic Mg surface after CO2-AP-treatment,(c,d)Plane and cross-sectional SEM imges of 5-swing CO2-AP treated Mg sample.

    To further investigate the corrosion behavior of ‘hydrophilic’ CO2-AP-treated surfaces,EIS measurements were conducted,and some impedance spectra are plotted in Fig.4 for 1 h and longer immersion times in 3.5 wt.% NaCl.For 1 h immersion,the treated Mg sample showed two impedance arcs,and the total impedance scale was greater in the treated sample than the untreated Mg.In longer immersion times,the impedance scale decreased significantl for both the treated and untreated Mg samples,but the treated sample still exhibited a greater impedance scale.The resistance values,R2or R2+R2’,determined from impedance data fittin are summarized in Table 2.The resistance values suggest that even hydrophilic CO2-AP-treated Mg surface could exhibit slight corrosion mitigation.

    Table 2 OCP values and EIS fittin results of the impedance spectra presented in Fig.4.

    H2collection measurements were also performed for untreated and CO2-AP-treated Mg samples,and the results are presented in Fig.5.Note that the two separate comparisons are from the different exposed area (0.98 vs.0.833 cm2) and baseline (untreated) surface (120 vs.600 grit SiC finish)Compared to the untreated Mg surfaces,both hydrophobic and hydrophilic CO2-AP-treated Mg surfaces produced lower H2vol during the immersion in 3.5 wt.% NaCl,indicating that CO2-AP-treatment of Mg would result in long-term corrosion mitigation regardless of the surface wetting property(hydrophobic or hydrophilic).

    4.3.Protective carbonized MgO/MgCO3 surface fil

    Fig.3.Polarization curves of (a) untreated and hydrophobic CO2-AP treated Mg surfaces with 0.98 cm2 area and (b) untreated and hydrophilic CO2-AP treated Mg surfaces with 0.833 cm2 area.3.5 wt.% NaCl was used.

    Fig.4.Impance plots of (a) untreated and hydrophilic CO2-AP Mg surfaces witrh different immersion times and (b) two impedance plots in (a) shown in smaller impedance axes.3.5 wt.% NaCl was used.The impedance values at low frequencies,often in inductive domain,were noisy to be precisley analyzed so that they were not plotted.Although more preciase corrosion rates can be derived by analyzing low frequency inductive impedances [33],fair comparison of corrosion resistance was previously made without considering the low frequency impedances [34,35].

    Fig.5.H2 collection volume vs.immersion time plots: (a) untreated and hydrophobic CO2-AP treated Mg surfaces with 0.98 cm2 area and (b) untreated and hydrophilic CO2-AP treated Mg surfaces with 0.833 cm2 area.

    Fig.6.SEM-Energy-dispersive X-ray spectroscopy images of CO2-AP–treated Mg surface.The Mg surface was treated with 1 swing.(a) Top view of SEM SE image and corresponding elemental maps (b–d).

    The corrosion mitigation effect is presumably originated from the chemical composition and crystalline changes (the MgCO3/MgO mixture) resulting from the CO2-AP treatment.Hydrophobicity,observed from some CO2-AP treated samples,should also contribute significantl to corrosion mitigation and could be originated from surface roughness associated with the nanoparticle (i.e.,the MgO/MgCO3intermixed)pillars (see Fig.2d) We suggest that chemically active CO2plasma molecules adsorbed on the MgO nanoparticle surfaces,then formed carbonated layers.In principle,MgCO3can be formed via the chemical reaction,MgO+CO2→MgCO3.At ambient temperature and pressure,Gibbs free energy change of this chemical reaction is ?21.4 KJ mol?1(i.e.,the exothermic and thermodynamically favorable reaction even at room temperature).Also,the MgO particles have diverse basic sites where the acidic CO2molecules be chemically/physically adsorbed on the surface,resulting in various carbonate surface species,including bicarbonate,bidentate,and unidentate carbonates [36].The theoretical CO2uptake capacity of MgO is 109 wt.% (1.09 of CO2/1 g of MgO sorbent).Nano-scale mesoporous MgO can show an adsorption capacity of ~8–10 wt.% at ~25–30 °C [37,38].

    We also consider that heat generation could be another factor to accelerate CO2adsorption at the surface.It was reported that the substrate temperature under atmospheric plasma could reach to 200 °C during the treatment with the plasma [24].A formation of Mg hydroxycarbonate may be considered as another route for MgCO3.Because of the affinit between CO2and Mg(OH)2,hydromagnesite [i.e.,Mg5(CO3)4(OH2)·4H2O]was formed in humid air with CO2-rich conditions [39].The hydromagnesite becomes MgCO3in dry conditions.The CO2-driven Mg hydroxycarbonate is less soluble in water and could provide a protective and electrically insulating layer on the surface,possibly suppressing corrosion by blocking the anodic and/or cathodic sites [22].

    The CO2-AP–treated surface exhibited a higher carbon content.Fig.6 shows the top surface treated with 1-swing CO2-AP.The treated surface had more O and C contents,suggesting that the MgO adsorbed CO2or molecular C (e.g.,volatile organic carbons) from CO2-rich environment.Surface fil formed on 1-swing CO2-AP treatment was further characterized by cross-sectional DF-and HAADF-STEM with energy-dispersive X-ray spectroscopy analyses,as presented in Fig.7.It shows that agglomerated nanoparticles in pillar structures and underlying dense layer on Mg substrate,which is similar to the previous SEM image in Fig.2d.The elemental mapping of the top layer indicates that three major elements (Mg,O,and C) are well-dispersed throughout the layer,verifying the coexistence of MgO and MgCO3.

    X-ray diffraction(XRD)measurements in Fig.8 confirme that the phases of the resultant layer were primarily MgO and MgCO3.The peaks at 32.0,43.0,46.5,and 52.0° could be assigned to [104],[113],[202] and [024] planes in the hexagonal phase of MgCO3,whereas other peaks at 37.2 and 43.0° could be assigned to MgO.The surface carbonization was analyzed by XPS spectra of Mg,O,and C (Fig.9).The CO2-AP sample,as shown in Fig.9a,has a distinct Ccarbonate peak at 292 eV and a broad C–C peak fraction at 285 eV,whereas the untreated sample has a relatively low Ccarbonate peak,which is very similar to the reported pure Mg results [36].The untreated referential surface was generated by scratching on the same sample before the XPS examination to simulate the referential corrosion test sample and identify the carbon adsorption.Notably,all carbon-related peaks in the reference scratched area,as shown in Fig.9e,disappear after 15 s of Ar-sputtering.

    Fig.7.High-magnificatio cross-sectional DF-and HAADF-STEM images of the CO2-AP–treated Mg along with corresponding energy-dispersive X-ray spectroscopy elemental mapping showing the distribution of Mg,C,and O in the top treated layer.In the Yellow dotted box area of DF4 image,the compostions of C,Mg and O were ~15,~40 and 35 wt.%,respectively.

    Fig.8.XRD pattern for the CO2-AP–treated Mg (fi e-times swing).XRD patterns of untreated Mg were measured at two angles (0° and 90°) to avoid the crystalline directional effect.

    A pure Mg fil was reported to form a ~4 nm oxide/hydroxide film after exposure to the ambient atmosphere within 1 day,the oxide fil thickness increased to 10 nm over three weeks,and the air-formed oxide/hydroxide film on Mg adsorbed a significan amount of carbon from the ambient air to form MgCO3[40].The same study showed that upon exposure to the ambient atmosphere,the thin oxide/hydroxide fil Mg surface adsorbed volatile organic compounds to form predominantly bidentate-bound carboxylate salts,which cause a significan surface change from hydrophilic to hydrophobic.Also,another reported XPS study on carbonization of MgO particles confirme that the formation of carbonate species on MgO particles increases significantl with the increase in the CO2concentration [41].

    Fig.9b shows two O 1s peaks are easily distinguished in the XPS spectra of two comparable samples.The literature has shown that the peaks related to MgO have binding energies at ~530.0–531.0 eV,MgOH have binding energies at ~530.0–533.2 eV,and MgCO3present peaks at ~533.2–533.5 eV [42,43].Similarly,a strong O–C peak in the O 1s spectra of the treated sample is present after the long-time sputtering (Fig.9b),while the O–C peak from the untreated sample significantl reduced after the short-time sputtering(Fig.9f).The Mg 2p and Mg Auger indicate that the CO2-AP treatment formed thicker MgCO3/MgO film than the untreated sample as the peaks of respective species were higher in the treated sample (Fig.9c,d,g and h).Additionally,XPS depth profile in Fig.10 show that the C–C/?CO3amount of the treated surface is significantl higher than the untreated one.

    Fig.9.XPS core energy spectra for (a–d) a CO2-AP treated Mg (with 5-swing) compared with (e–h) untreated Mg.

    Fig.10.XPS depth profilin of carbon component(C–C/?CO3)for the CO2-AP–treated Mg.

    5.Molecular dynamics simulations

    Due to analytical limitation of surface characterization at the atomic scale surface layer,we performed classical MD simulations based on empirical potentials [26,27,31] to understand the origin of the hydrophobicity at the atomic level for the hydrophobic treated samples.In the beginning,we investigated the interaction of water molecules on the pristine MgCO3(magnesite) surface using a force fiel [26,27] established for alkaline-earth carbonate and water systems.Our model system (Fig.11a) consists of 1,000 water molecules on a MgCO3slab of 20,800 atoms and a 13.0thickness with (104) facet,which corresponds to the strongest magnesite peak of our XRD analysis,as shown in Fig.8.The thermally equilibrated structure shows the contact angle between water molecules and the surface of ~20°,which confirm the hydrophilicity of magnesite;the experimental contact angle reported in the literature is ~10.4° [44].

    We hypothesize that the surface roughness of a hydrophobic surface is the origin of the superhydrophobicity of our system [45].Unlike the fla magnesite or MgO,which are hydrophilic,hydrophobic carbon can add surface roughness to the magnesite surface when the carbon molecules bind on the surface during the CO2-AP treatment.The formation of a carbon layer is confirme (Fig.10) but characterizing its structural properties is challenging.Therefore,we model the crystalline and amorphous carbon structures by performing MD simulations: a simple cubic carbon structure is annealed at a high temperature(~8,000 K),followed by a rapid cooling to room temperature,resulting in an amorphous carbon with the density of graphite (2.267 g/cm3) [30] (see the section: 3.Computational Modeling).

    Fig.11.(a) MD simulation modeling a system of a water droplet (1,000 water molecules) on a (104) MgCO3 slab;(b) a water droplet (~4,000 water molecules) on 6 layers of graphene layers,and (c) on an amorphous carbon slab.Here,the red,white,gray,and green balls represent O,H,C,and Mg atoms,respectively.The atomic simulations confir the hydrophilicity of the pristine magnesium carbonate surface and hydrophobicity or superhydrophobicity on the graphene surface and the amorphous carbon surface.

    Next,our MD simulations compare the contact angles of a water droplet on graphite and amorphous carbon with same atomic densities.We obtain a contact angle of the water nanodroplet on graphite of 94.3 ± 2.5° (Fig.11b),which corresponds well with experimental values in the range of 90° to 95° [46].On the amorphous carbon surface,the contact angle increases to 138.3 ± 4.0° (Fig.11c).This angle is slightly less than the superhydrophobic angle,150°,but considering our slab is a fla surface,a roughened amorphous carbon surface would show superhydrophobicity.Indeed,the CO2-AP surfaces show rough surfaces consisting of nanoparticle accumulation (Figs.2d and 7a).

    The reason for the increase of the contact angle may be explained by the relatively weakened cohesive force between the water molecules and carbon atoms.The amorphous carbon slab would have relatively less surface atomic density than that of graphite due to defects,and this would decrease the overall cohesive force between water molecules and carbon atoms.Notably,for the system-size independency,we adapted an approach from literature that is free of a systemsize effect by introducing thin slab water droplets instead of dome-shaped water droplets [29].In the literature,the authors found more than 2,000 water molecules show size independency,and more than 4 layers of graphene showed graphitelike behavior;therefore,we choose 4,000 water molecules on 6 layers of graphene and a similar thickness of the amorphous carbon slab.

    6.Conclusion

    A microstructured MgCO3/MgO based surface was formed on commercially pure Mg substrates via simple CO2-AP treatment.The treated Mg samples exhibited either hydrophobic or hydrophilic surfaces.In hydrophobic treated surfaces,reduction of Mg corrosion was remarkable as confirme from significantl low anodic and cathodic currents during polarization measurements.In hydrophilic treated surfaces,moderated increase of corrosion resistance was confirme by EIS measurements.The results from H2collection measurements showed that both hydrophobic and hydrophilic treated surfaces can exhibit long-term corrosion mitigation.

    Electron-microscopic and X-ray based characterization suggest that the treated surfaces had a unique structure,composed of top nano pillars and dense bottom layer and two major phases,MgCO3and MgO,with excessive carbon.Presence of carbon-carbon bonds,possibly a cause of hydrophobicity,were confirme on the very top surface of a CO2-AP treated Mg.The MD simulations demonstrated that superhydrophobicity could originate from the amorphous carbon layers formed on the treated Mg surface.

    Declaration of Competing InterestThe authors declare that they have no known competing financia interests or personal relationships that could have appeared to influenc the work reported in this paper.

    Acknowledgments

    This work was supported by the US Department of Energy’s (DOE’s) Offic of Energy Efficien y and Renewable Energy,Vehicle Technology Office Powertrain Materials Core and Light Metals Core Programs.Sample characterization is also partially supported by the Technology Commercialization Fund Fiscal Year 2020 of DOE’s Offic of Technology Transitions and by the Creative Materials Discovery Program through the National Research Foundation of Korea,with computational modeling of amorphous carbons funded by the Ministry of Science,ICT and Future Planning(NRF-2016M3D1A1919181).The research was conducted at ORNL,which is managed by UT Battelle LLC for DOE under contract DE-AC05-00OR22725.Computational modeling and part of the materials characterization (SEM and XRD)were performed at the Center for Nanophase Materials Sciences (proposal ID: CNMS 2021-A-00626),which is sponsored at ORNL by DOE’s Scientifi User Facilities Division.This research used resources of the Oak Ridge Leadership Computing Facility and the National Energy Research Scientifi Computing Center,which are DOE Offic of Science user facilities.Peter Yancey (Atmospheric Plasma Solutions,11301 Penny Road,Cary,NC 27518,United States) provided significan technical support and supervised the production of plasma-treated Mg samples.

    Appendix

    Fig.A1.(a) Plasma treatment inside a CO2-enriched chamber with a 3-inch wide observation window,(b) inline automated CO2 atmospheric plasma treatment on an epoxy-mounted Mg disk specimen (?≈0.5 inch),(c) the schematic diagram shows a proposed mechanism;CO2 plasma converts the Mg surface to an MgO/MgCO3 intermixed-particulates layer.

    Fig.A2.Exposure surface photos of the Mg specimens (7×14 mm2 or 0.98 cm2) after 3.5 wt.% NaCl solution immersion for 20~24 h: (a) untreated (120 and 600 grit) and (b) hydrophobic CO2-AP treated Mg surfaces (5-swings).The CO2-AP treated (hydrophobic) surface has much smaller areas of corrosion attack,compared with the untreated Mg.(c) CO2-AP swing on Mg sample.

    Fig.A3.Two equivanent circuit models with (a) single and (b) double R//CPE used for fittin the impedance data presented in Fig.4.R1 is solution resistance,R2 corresponds to charge transfer resistance at the interface of Mg alloy and electrolyte (3.5 wt.% NaCl) and R2’ represents resistance of porous surface layer.CPE and CPE’ correspond to non-ideal capacitances of double layer and porous surface layer,respectively.

    久久久久久久午夜电影| 法律面前人人平等表现在哪些方面| 美女大奶头视频| 亚洲av成人不卡在线观看播放网| www.999成人在线观看| 超碰av人人做人人爽久久 | 特大巨黑吊av在线直播| 国产精品乱码一区二三区的特点| 又紧又爽又黄一区二区| 成人高潮视频无遮挡免费网站| 啪啪无遮挡十八禁网站| 国产真实伦视频高清在线观看 | 欧美黄色片欧美黄色片| 我要搜黄色片| 久久久久性生活片| 少妇的丰满在线观看| 五月玫瑰六月丁香| 欧美日韩黄片免| 一本综合久久免费| 国产精品国产高清国产av| 久久久久久久久中文| 国产精品国产高清国产av| 美女大奶头视频| 免费看美女性在线毛片视频| 色综合欧美亚洲国产小说| 搡老熟女国产l中国老女人| 国产乱人视频| av在线天堂中文字幕| 少妇人妻精品综合一区二区 | 高潮久久久久久久久久久不卡| 99视频精品全部免费 在线| 在线观看一区二区三区| 久久精品91无色码中文字幕| 天堂影院成人在线观看| 老汉色∧v一级毛片| 久久中文看片网| 欧美xxxx黑人xx丫x性爽| 蜜桃久久精品国产亚洲av| 少妇熟女aⅴ在线视频| 无遮挡黄片免费观看| 亚洲性夜色夜夜综合| 午夜激情欧美在线| 无限看片的www在线观看| 免费看a级黄色片| 久久久精品欧美日韩精品| 午夜激情欧美在线| 日本熟妇午夜| 婷婷六月久久综合丁香| 老鸭窝网址在线观看| 首页视频小说图片口味搜索| 久久久久国产精品人妻aⅴ院| 18禁黄网站禁片免费观看直播| 免费大片18禁| 色av中文字幕| 午夜两性在线视频| 国产91精品成人一区二区三区| 国模一区二区三区四区视频| 国产色婷婷99| 神马国产精品三级电影在线观看| 69人妻影院| 亚洲欧美精品综合久久99| 日韩国内少妇激情av| 欧美av亚洲av综合av国产av| 欧美黑人巨大hd| 亚洲男人的天堂狠狠| 国产精品野战在线观看| 波多野结衣巨乳人妻| 国产成人啪精品午夜网站| xxxwww97欧美| 久99久视频精品免费| 免费在线观看成人毛片| 99在线视频只有这里精品首页| 国产精品亚洲av一区麻豆| 国产成+人综合+亚洲专区| 久久国产乱子伦精品免费另类| 我要搜黄色片| 欧美另类亚洲清纯唯美| 亚洲欧美日韩高清在线视频| 色综合欧美亚洲国产小说| 免费观看人在逋| 亚洲,欧美精品.| 亚洲片人在线观看| 亚洲精品国产精品久久久不卡| 精品欧美国产一区二区三| 国内精品久久久久精免费| 人妻夜夜爽99麻豆av| 1000部很黄的大片| 亚洲18禁久久av| 久久久色成人| 国产伦一二天堂av在线观看| 日韩人妻高清精品专区| 国产精品精品国产色婷婷| 亚洲一区高清亚洲精品| 亚洲av二区三区四区| 国产成人系列免费观看| 国产一区二区三区在线臀色熟女| 久久久久久久精品吃奶| 在线免费观看的www视频| 悠悠久久av| 国语自产精品视频在线第100页| 国产视频一区二区在线看| 国产av在哪里看| 欧美午夜高清在线| 午夜免费观看网址| 成年女人永久免费观看视频| aaaaa片日本免费| 免费av观看视频| 麻豆久久精品国产亚洲av| 日韩人妻高清精品专区| 久久人妻av系列| 看片在线看免费视频| 白带黄色成豆腐渣| 黄色日韩在线| 又粗又爽又猛毛片免费看| 日韩精品青青久久久久久| 在线免费观看的www视频| 欧美性猛交黑人性爽| 久久亚洲真实| 法律面前人人平等表现在哪些方面| 精品久久久久久久人妻蜜臀av| 在线观看免费视频日本深夜| 午夜免费成人在线视频| 香蕉av资源在线| 国产精品综合久久久久久久免费| 最近最新中文字幕大全电影3| 欧美另类亚洲清纯唯美| 99久久精品一区二区三区| 国产精品99久久久久久久久| 欧美区成人在线视频| 中文字幕av成人在线电影| 久久这里只有精品中国| 哪里可以看免费的av片| 成人特级av手机在线观看| 一边摸一边抽搐一进一小说| 国产精品精品国产色婷婷| 亚洲精品亚洲一区二区| 18禁黄网站禁片免费观看直播| 最新中文字幕久久久久| 亚洲内射少妇av| 性色avwww在线观看| 成人18禁在线播放| 一本综合久久免费| 又粗又爽又猛毛片免费看| 国产精品香港三级国产av潘金莲| 精品久久久久久久毛片微露脸| 嫁个100分男人电影在线观看| 午夜久久久久精精品| 国产精品电影一区二区三区| 韩国av一区二区三区四区| 亚洲精品在线观看二区| 床上黄色一级片| 最近最新中文字幕大全免费视频| 99精品久久久久人妻精品| 美女被艹到高潮喷水动态| 精品久久久久久久人妻蜜臀av| 99视频精品全部免费 在线| 欧美黄色片欧美黄色片| av女优亚洲男人天堂| 在线播放国产精品三级| 成人高潮视频无遮挡免费网站| 日韩人妻高清精品专区| 国产一区二区在线av高清观看| 免费大片18禁| 亚洲久久久久久中文字幕| 国产亚洲精品一区二区www| 亚洲乱码一区二区免费版| 国产黄色小视频在线观看| 免费在线观看成人毛片| 免费看光身美女| 亚洲国产精品sss在线观看| 噜噜噜噜噜久久久久久91| 亚洲国产色片| 男插女下体视频免费在线播放| 午夜福利在线观看免费完整高清在 | 中亚洲国语对白在线视频| 午夜老司机福利剧场| 国产视频内射| 99精品在免费线老司机午夜| 午夜视频国产福利| 国产午夜精品论理片| 亚洲一区高清亚洲精品| 国产在线精品亚洲第一网站| 手机成人av网站| 搡女人真爽免费视频火全软件 | 18禁在线播放成人免费| 淫秽高清视频在线观看| 熟女人妻精品中文字幕| 人妻丰满熟妇av一区二区三区| 色综合婷婷激情| 成人av一区二区三区在线看| 日本 欧美在线| xxx96com| 国产精品美女特级片免费视频播放器| 亚洲熟妇熟女久久| 免费观看人在逋| 国产一区二区激情短视频| 成人特级av手机在线观看| 一个人观看的视频www高清免费观看| 日韩有码中文字幕| av黄色大香蕉| 91在线观看av| 国产一区二区三区视频了| 99热只有精品国产| 真人一进一出gif抽搐免费| 夜夜躁狠狠躁天天躁| 在线十欧美十亚洲十日本专区| av天堂中文字幕网| 啦啦啦免费观看视频1| 国产欧美日韩精品亚洲av| 少妇裸体淫交视频免费看高清| 99久久无色码亚洲精品果冻| 亚洲最大成人手机在线| 中出人妻视频一区二区| 在线播放无遮挡| av天堂中文字幕网| 国产免费av片在线观看野外av| av女优亚洲男人天堂| 国产熟女xx| www日本黄色视频网| 69av精品久久久久久| 一边摸一边抽搐一进一小说| 日韩亚洲欧美综合| 高潮久久久久久久久久久不卡| 我的老师免费观看完整版| 欧美黄色片欧美黄色片| 99riav亚洲国产免费| 怎么达到女性高潮| 精品国产亚洲在线| 久久精品夜夜夜夜夜久久蜜豆| 免费高清视频大片| 日日干狠狠操夜夜爽| 亚洲国产色片| 亚洲精品乱码久久久v下载方式 | 久久久久免费精品人妻一区二区| 欧美三级亚洲精品| 啦啦啦免费观看视频1| 欧美在线一区亚洲| 国产单亲对白刺激| 国产午夜福利久久久久久| 日日干狠狠操夜夜爽| 美女被艹到高潮喷水动态| 国产伦精品一区二区三区视频9 | 悠悠久久av| 色播亚洲综合网| 高潮久久久久久久久久久不卡| 亚洲国产精品sss在线观看| 麻豆国产av国片精品| 内射极品少妇av片p| 国产免费av片在线观看野外av| 十八禁人妻一区二区| 日本撒尿小便嘘嘘汇集6| 亚洲国产中文字幕在线视频| 国产熟女xx| 久久久久久久亚洲中文字幕 | 国产伦精品一区二区三区四那| 久久久久国产精品人妻aⅴ院| 国内精品久久久久精免费| 特级一级黄色大片| 国产午夜精品论理片| 久久欧美精品欧美久久欧美| 亚洲熟妇熟女久久| 在线观看免费午夜福利视频| 亚洲国产中文字幕在线视频| 麻豆成人av在线观看| 一级a爱片免费观看的视频| 18禁国产床啪视频网站| 国产免费av片在线观看野外av| 丰满乱子伦码专区| 中国美女看黄片| 国产在线精品亚洲第一网站| 日本免费一区二区三区高清不卡| 天天添夜夜摸| 成年女人毛片免费观看观看9| 亚洲av中文字字幕乱码综合| 网址你懂的国产日韩在线| 又粗又爽又猛毛片免费看| 日韩欧美一区二区三区在线观看| 欧美日韩乱码在线| 精品无人区乱码1区二区| 中文字幕高清在线视频| 婷婷精品国产亚洲av在线| 露出奶头的视频| 成人性生交大片免费视频hd| xxxwww97欧美| 日日摸夜夜添夜夜添小说| 成人特级黄色片久久久久久久| av欧美777| 丝袜美腿在线中文| 国产aⅴ精品一区二区三区波| 老熟妇仑乱视频hdxx| 日韩欧美精品免费久久 | 黑人欧美特级aaaaaa片| 国产精品98久久久久久宅男小说| 精品一区二区三区视频在线观看免费| 久久精品国产自在天天线| 午夜福利免费观看在线| 精品熟女少妇八av免费久了| 国产aⅴ精品一区二区三区波| 少妇熟女aⅴ在线视频| 国产蜜桃级精品一区二区三区| 亚洲成av人片免费观看| 天堂网av新在线| 99riav亚洲国产免费| 91久久精品国产一区二区成人 | 欧洲精品卡2卡3卡4卡5卡区| 看免费av毛片| 欧美绝顶高潮抽搐喷水| 精品一区二区三区视频在线观看免费| 国产精品,欧美在线| 日韩欧美精品免费久久 | 国产午夜精品论理片| 一区福利在线观看| 淫秽高清视频在线观看| 欧美成人一区二区免费高清观看| 欧美黄色淫秽网站| 精品午夜福利视频在线观看一区| 成人永久免费在线观看视频| 免费看a级黄色片| www.999成人在线观看| 国产一区二区三区在线臀色熟女| 欧美丝袜亚洲另类 | 欧美成人免费av一区二区三区| 国语自产精品视频在线第100页| 午夜福利高清视频| 久久精品国产清高在天天线| 有码 亚洲区| 天美传媒精品一区二区| 久久草成人影院| 国产精品乱码一区二三区的特点| 热99在线观看视频| 午夜免费激情av| 精品久久久久久久人妻蜜臀av| 18禁黄网站禁片免费观看直播| 国产国拍精品亚洲av在线观看 | 亚洲av免费高清在线观看| 又爽又黄无遮挡网站| 网址你懂的国产日韩在线| 国产主播在线观看一区二区| 亚洲欧美日韩高清专用| 国产伦在线观看视频一区| 亚洲精品日韩av片在线观看 | 免费看美女性在线毛片视频| 日本一二三区视频观看| 内地一区二区视频在线| 欧美最新免费一区二区三区 | 少妇高潮的动态图| 国产一区二区三区视频了| 亚洲专区国产一区二区| 色在线成人网| 99国产综合亚洲精品| 国产成人aa在线观看| 搞女人的毛片| 久久中文看片网| 51国产日韩欧美| 精品国内亚洲2022精品成人| 一本综合久久免费| 成人av一区二区三区在线看| 久久精品国产亚洲av香蕉五月| 亚洲精品久久国产高清桃花| 欧美一级毛片孕妇| 女同久久另类99精品国产91| 波多野结衣高清作品| 美女 人体艺术 gogo| 美女被艹到高潮喷水动态| 亚洲电影在线观看av| 亚洲美女黄片视频| 男女视频在线观看网站免费| 观看免费一级毛片| 久久久久亚洲av毛片大全| 夜夜爽天天搞| 麻豆一二三区av精品| 少妇丰满av| 嫩草影院精品99| 女警被强在线播放| 淫妇啪啪啪对白视频| 1000部很黄的大片| 免费大片18禁| 观看免费一级毛片| 亚洲国产精品sss在线观看| 国产亚洲精品久久久com| 国产亚洲精品一区二区www| 国内揄拍国产精品人妻在线| 亚洲性夜色夜夜综合| 成年女人毛片免费观看观看9| 99久久99久久久精品蜜桃| 亚洲精品成人久久久久久| 一边摸一边抽搐一进一小说| 成人精品一区二区免费| www.色视频.com| 欧美色欧美亚洲另类二区| 白带黄色成豆腐渣| 在线视频色国产色| 最近最新中文字幕大全免费视频| 欧美日韩亚洲国产一区二区在线观看| 少妇高潮的动态图| 欧美午夜高清在线| 热99re8久久精品国产| 国语自产精品视频在线第100页| 亚洲熟妇熟女久久| 小蜜桃在线观看免费完整版高清| 国产精品女同一区二区软件 | 国产精品日韩av在线免费观看| 少妇的丰满在线观看| 网址你懂的国产日韩在线| 男女床上黄色一级片免费看| 麻豆成人午夜福利视频| 一区二区三区国产精品乱码| 国产不卡一卡二| 亚洲成人久久爱视频| h日本视频在线播放| 亚洲激情在线av| 免费观看的影片在线观看| 看免费av毛片| 精品一区二区三区视频在线观看免费| 国产精品精品国产色婷婷| 校园春色视频在线观看| 91久久精品电影网| 在线观看午夜福利视频| 国产熟女xx| 亚洲内射少妇av| 99热这里只有精品一区| 久久国产精品影院| 亚洲精品粉嫩美女一区| 日韩亚洲欧美综合| 一本综合久久免费| 18禁国产床啪视频网站| av女优亚洲男人天堂| av福利片在线观看| 成人特级av手机在线观看| 最近最新中文字幕大全电影3| 国产三级在线视频| 欧美成人一区二区免费高清观看| 天堂影院成人在线观看| 国产免费一级a男人的天堂| 18禁在线播放成人免费| 国产成人av激情在线播放| 夜夜看夜夜爽夜夜摸| 亚洲七黄色美女视频| 天堂av国产一区二区熟女人妻| 色综合亚洲欧美另类图片| 毛片女人毛片| 欧美日本视频| 国产91精品成人一区二区三区| 美女 人体艺术 gogo| 岛国在线观看网站| 日本精品一区二区三区蜜桃| 欧美黄色淫秽网站| 熟女电影av网| 亚洲av熟女| 成人精品一区二区免费| 久久精品国产亚洲av涩爱 | 91在线精品国自产拍蜜月 | 女生性感内裤真人,穿戴方法视频| 精品99又大又爽又粗少妇毛片 | tocl精华| 久久精品亚洲精品国产色婷小说| 精品乱码久久久久久99久播| 狂野欧美激情性xxxx| 亚洲精品在线观看二区| 特大巨黑吊av在线直播| 高清日韩中文字幕在线| 中文字幕久久专区| 亚洲精品美女久久久久99蜜臀| 中出人妻视频一区二区| 精品人妻偷拍中文字幕| 欧美成人a在线观看| 91九色精品人成在线观看| 香蕉av资源在线| 亚洲人与动物交配视频| 国产精品99久久99久久久不卡| 国产精品女同一区二区软件 | 宅男免费午夜| 怎么达到女性高潮| 欧美性感艳星| 欧美最新免费一区二区三区 | 老司机午夜福利在线观看视频| 婷婷六月久久综合丁香| 精品国内亚洲2022精品成人| 免费观看人在逋| 又紧又爽又黄一区二区| 美女cb高潮喷水在线观看| 精品不卡国产一区二区三区| 日韩av在线大香蕉| 韩国av一区二区三区四区| 变态另类丝袜制服| 九色国产91popny在线| 搡老岳熟女国产| 可以在线观看的亚洲视频| 中文亚洲av片在线观看爽| 成人午夜高清在线视频| 国产亚洲精品一区二区www| 听说在线观看完整版免费高清| 99久久精品一区二区三区| 国产精品一及| 国产成人啪精品午夜网站| 国产高清激情床上av| 亚洲av电影不卡..在线观看| 男女床上黄色一级片免费看| 亚洲一区二区三区不卡视频| 黄色女人牲交| 黄色片一级片一级黄色片| 国产一级毛片七仙女欲春2| 欧美一区二区国产精品久久精品| 日本黄色片子视频| 亚洲成人久久性| 日韩高清综合在线| 亚洲精品一卡2卡三卡4卡5卡| av视频在线观看入口| 熟妇人妻久久中文字幕3abv| 麻豆国产97在线/欧美| www.熟女人妻精品国产| 日韩 欧美 亚洲 中文字幕| 校园春色视频在线观看| 亚洲片人在线观看| 欧美日韩亚洲国产一区二区在线观看| 日本与韩国留学比较| 狂野欧美白嫩少妇大欣赏| 在线观看66精品国产| 91在线精品国自产拍蜜月 | av国产免费在线观看| 两人在一起打扑克的视频| 男人舔奶头视频| 亚洲欧美激情综合另类| 国产精品影院久久| 亚洲avbb在线观看| 哪里可以看免费的av片| 精品久久久久久久人妻蜜臀av| 午夜福利免费观看在线| 一级毛片高清免费大全| 国产精品自产拍在线观看55亚洲| 一本综合久久免费| 69人妻影院| 欧美日韩精品网址| 中亚洲国语对白在线视频| 国产亚洲精品久久久久久毛片| 一边摸一边抽搐一进一小说| 熟女电影av网| 十八禁人妻一区二区| 麻豆成人av在线观看| 精品免费久久久久久久清纯| 国产成+人综合+亚洲专区| 黄色视频,在线免费观看| 小蜜桃在线观看免费完整版高清| 亚洲av成人av| 久久伊人香网站| 免费av毛片视频| 天堂影院成人在线观看| 国产91精品成人一区二区三区| 最新美女视频免费是黄的| 亚洲熟妇中文字幕五十中出| 在线国产一区二区在线| 俄罗斯特黄特色一大片| 韩国av一区二区三区四区| 淫秽高清视频在线观看| 国产av麻豆久久久久久久| 国产国拍精品亚洲av在线观看 | 一级毛片高清免费大全| 亚洲av电影在线进入| 好看av亚洲va欧美ⅴa在| 国产国拍精品亚洲av在线观看 | 日韩欧美国产一区二区入口| e午夜精品久久久久久久| 丁香欧美五月| 亚洲久久久久久中文字幕| www日本黄色视频网| 欧美一级毛片孕妇| 女人高潮潮喷娇喘18禁视频| 很黄的视频免费| 久久久久久九九精品二区国产| 嫁个100分男人电影在线观看| 日韩中文字幕欧美一区二区| 欧美乱妇无乱码| 国产99白浆流出| 麻豆国产av国片精品| 精品电影一区二区在线| 久久久久性生活片| 亚洲,欧美精品.| 亚洲最大成人中文| 国产精品美女特级片免费视频播放器| 色综合亚洲欧美另类图片| 脱女人内裤的视频| 90打野战视频偷拍视频| 亚洲国产欧美网| 亚洲性夜色夜夜综合| 在线观看舔阴道视频| 舔av片在线| 中文资源天堂在线| 九九热线精品视视频播放| 99精品在免费线老司机午夜| 国产精品三级大全| 国产亚洲精品综合一区在线观看| 99视频精品全部免费 在线| 十八禁人妻一区二区| 日韩免费av在线播放| www国产在线视频色| 黄色女人牲交| 亚洲欧美日韩高清在线视频| 亚洲18禁久久av| 真人做人爱边吃奶动态| 亚洲狠狠婷婷综合久久图片| 一本综合久久免费| 久久6这里有精品| 免费无遮挡裸体视频| 一区二区三区免费毛片| 日本成人三级电影网站| 1024手机看黄色片| 91字幕亚洲| 99久久综合精品五月天人人| 色吧在线观看| 别揉我奶头~嗯~啊~动态视频| 久久久久国产精品人妻aⅴ院| 一本久久中文字幕| 日韩精品青青久久久久久| 亚洲av不卡在线观看| 色吧在线观看|