• <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.

    国产在视频线精品| 欧美人与性动交α欧美精品济南到| 俄罗斯特黄特色一大片| 国产成人精品无人区| 19禁男女啪啪无遮挡网站| 女人被躁到高潮嗷嗷叫费观| 国产精品亚洲一级av第二区| 一级a爱视频在线免费观看| 少妇 在线观看| 激情视频va一区二区三区| 久久婷婷成人综合色麻豆| 亚洲伊人色综图| 久久狼人影院| 天天操日日干夜夜撸| 成人国语在线视频| 精品卡一卡二卡四卡免费| 女人被躁到高潮嗷嗷叫费观| 亚洲av日韩精品久久久久久密| 亚洲熟女精品中文字幕| 免费一级毛片在线播放高清视频 | 波多野结衣一区麻豆| 搡老乐熟女国产| av免费在线观看网站| 亚洲国产毛片av蜜桃av| 午夜福利免费观看在线| e午夜精品久久久久久久| 正在播放国产对白刺激| 欧美乱妇无乱码| 亚洲第一av免费看| 亚洲av成人不卡在线观看播放网| 亚洲精品成人av观看孕妇| 亚洲色图av天堂| 一区二区三区激情视频| 999久久久精品免费观看国产| 超碰97精品在线观看| 亚洲 国产 在线| 久热这里只有精品99| 巨乳人妻的诱惑在线观看| 精品国产亚洲在线| 亚洲 国产 在线| 国产伦理片在线播放av一区| 在线看a的网站| 欧美亚洲 丝袜 人妻 在线| 成年女人毛片免费观看观看9 | 人人妻人人澡人人看| 亚洲国产av影院在线观看| 欧美久久黑人一区二区| 高清欧美精品videossex| 999精品在线视频| 免费久久久久久久精品成人欧美视频| 两人在一起打扑克的视频| 黄色丝袜av网址大全| 高清在线国产一区| 一级,二级,三级黄色视频| 肉色欧美久久久久久久蜜桃| 久久久水蜜桃国产精品网| 亚洲精品美女久久av网站| 久久久国产精品麻豆| 国产不卡一卡二| 国产精品国产av在线观看| 80岁老熟妇乱子伦牲交| 国产在线视频一区二区| 国产黄色免费在线视频| 精品熟女少妇八av免费久了| 色婷婷av一区二区三区视频| 久久久久精品国产欧美久久久| 丁香六月欧美| 超碰成人久久| 黑人巨大精品欧美一区二区蜜桃| 日韩人妻精品一区2区三区| 亚洲熟女毛片儿| 欧美性长视频在线观看| 一区二区三区乱码不卡18| 高潮久久久久久久久久久不卡| 国产aⅴ精品一区二区三区波| 亚洲一卡2卡3卡4卡5卡精品中文| 黑丝袜美女国产一区| 国产男女超爽视频在线观看| 国产午夜精品久久久久久| 国产一卡二卡三卡精品| 俄罗斯特黄特色一大片| 人人妻人人澡人人爽人人夜夜| 两个人免费观看高清视频| 精品亚洲成a人片在线观看| 韩国精品一区二区三区| 久久人人97超碰香蕉20202| 婷婷成人精品国产| 欧美 日韩 精品 国产| 婷婷成人精品国产| 老熟妇乱子伦视频在线观看| 国产成+人综合+亚洲专区| 免费一级毛片在线播放高清视频 | 麻豆国产av国片精品| 如日韩欧美国产精品一区二区三区| 桃红色精品国产亚洲av| 麻豆av在线久日| 啪啪无遮挡十八禁网站| 精品少妇黑人巨大在线播放| 欧美黑人精品巨大| 午夜福利欧美成人| 久久精品熟女亚洲av麻豆精品| 久久人人97超碰香蕉20202| 国产精品久久久久久精品古装| 一边摸一边做爽爽视频免费| 成年版毛片免费区| 午夜日韩欧美国产| 精品久久蜜臀av无| 日韩人妻精品一区2区三区| 亚洲视频免费观看视频| 大片免费播放器 马上看| a在线观看视频网站| avwww免费| 精品国产一区二区三区四区第35| 老熟妇仑乱视频hdxx| 超碰97精品在线观看| 超碰97精品在线观看| 精品国产亚洲在线| 一边摸一边做爽爽视频免费| 久久国产精品影院| 美女福利国产在线| 亚洲色图综合在线观看| 9色porny在线观看| 欧美成人免费av一区二区三区 | 日韩欧美一区视频在线观看| 少妇精品久久久久久久| 亚洲中文字幕日韩| 国产精品电影一区二区三区 | 一进一出抽搐动态| 高清黄色对白视频在线免费看| 电影成人av| 亚洲色图 男人天堂 中文字幕| 香蕉国产在线看| 桃花免费在线播放| www日本在线高清视频| 人妻 亚洲 视频| 国产成人影院久久av| 久久热在线av| 大型黄色视频在线免费观看| 国产成人精品无人区| www.熟女人妻精品国产| 午夜精品国产一区二区电影| 亚洲精品中文字幕一二三四区 | 国产精品.久久久| 久久久精品国产亚洲av高清涩受| 国产精品久久久久久精品电影小说| 日韩欧美一区视频在线观看| 99香蕉大伊视频| 黑人操中国人逼视频| 丁香欧美五月| a在线观看视频网站| 人人妻人人添人人爽欧美一区卜| 免费黄频网站在线观看国产| 最新的欧美精品一区二区| 高潮久久久久久久久久久不卡| 日韩欧美国产一区二区入口| 交换朋友夫妻互换小说| 久久午夜综合久久蜜桃| 丝袜在线中文字幕| 满18在线观看网站| 看免费av毛片| 丁香六月天网| 国产色视频综合| 久久 成人 亚洲| 91成人精品电影| 亚洲精华国产精华精| 91老司机精品| 亚洲成人免费av在线播放| 国产福利在线免费观看视频| 淫妇啪啪啪对白视频| 黑人巨大精品欧美一区二区蜜桃| 母亲3免费完整高清在线观看| 91老司机精品| tube8黄色片| 国产福利在线免费观看视频| 久久 成人 亚洲| 一个人免费看片子| 久久性视频一级片| 久久精品亚洲熟妇少妇任你| 亚洲精品美女久久av网站| 亚洲国产av影院在线观看| 曰老女人黄片| 黄色视频不卡| 欧美日韩亚洲国产一区二区在线观看 | 男女高潮啪啪啪动态图| 成人国产一区最新在线观看| 老鸭窝网址在线观看| 色视频在线一区二区三区| 女人高潮潮喷娇喘18禁视频| 最新的欧美精品一区二区| 日日爽夜夜爽网站| 80岁老熟妇乱子伦牲交| 十分钟在线观看高清视频www| 欧美精品亚洲一区二区| 亚洲中文av在线| 午夜激情久久久久久久| 国产高清视频在线播放一区| 老鸭窝网址在线观看| 日日摸夜夜添夜夜添小说| 欧美精品高潮呻吟av久久| 免费看十八禁软件| 中文字幕精品免费在线观看视频| 日本精品一区二区三区蜜桃| 91麻豆av在线| 91字幕亚洲| 亚洲av日韩在线播放| 国产在线精品亚洲第一网站| 午夜福利欧美成人| 久久久精品国产亚洲av高清涩受| 国产日韩欧美视频二区| 日本a在线网址| 久久国产精品人妻蜜桃| tube8黄色片| 国产精品成人在线| 国产精品国产高清国产av | 国产av又大| 香蕉久久夜色| 精品人妻熟女毛片av久久网站| 亚洲精品乱久久久久久| 每晚都被弄得嗷嗷叫到高潮| 欧美精品人与动牲交sv欧美| 亚洲专区字幕在线| 日本黄色日本黄色录像| 国产男女内射视频| 91成人精品电影| 91国产中文字幕| 91九色精品人成在线观看| 久热这里只有精品99| 国产精品美女特级片免费视频播放器 | videosex国产| 国产人伦9x9x在线观看| 国产精品麻豆人妻色哟哟久久| 国产色视频综合| 99国产精品一区二区三区| 高清视频免费观看一区二区| 亚洲欧洲日产国产| 国产日韩欧美在线精品| 亚洲伊人久久精品综合| 国产一区二区三区视频了| 麻豆乱淫一区二区| 亚洲精品中文字幕一二三四区 | 欧美日韩亚洲综合一区二区三区_| av一本久久久久| 成人国产一区最新在线观看| 国产黄频视频在线观看| 欧美日韩精品网址| 51午夜福利影视在线观看| 久久热在线av| 国产高清激情床上av| 三级毛片av免费| 少妇 在线观看| 天天添夜夜摸| 一级毛片精品| 宅男免费午夜| 老司机福利观看| 少妇粗大呻吟视频| 十八禁人妻一区二区| 女人爽到高潮嗷嗷叫在线视频| 蜜桃在线观看..| 国产成人免费观看mmmm| 老熟妇仑乱视频hdxx| 亚洲av欧美aⅴ国产| 亚洲色图综合在线观看| 日本精品一区二区三区蜜桃| 色尼玛亚洲综合影院| 嫩草影视91久久| 国产精品偷伦视频观看了| 无限看片的www在线观看| 人人妻人人爽人人添夜夜欢视频| 精品一区二区三区四区五区乱码| 亚洲自偷自拍图片 自拍| 亚洲av美国av| 三级毛片av免费| 国产在视频线精品| 91成年电影在线观看| 最近最新中文字幕大全电影3 | 十八禁网站网址无遮挡| 制服人妻中文乱码| 美女视频免费永久观看网站| 亚洲人成77777在线视频| av又黄又爽大尺度在线免费看| 国产视频一区二区在线看| 91精品三级在线观看| 99久久精品国产亚洲精品| 亚洲精品成人av观看孕妇| 久久久水蜜桃国产精品网| 欧美黄色淫秽网站| 亚洲欧美日韩另类电影网站| 啦啦啦 在线观看视频| 成年女人毛片免费观看观看9 | 在线观看舔阴道视频| 亚洲少妇的诱惑av| 男女之事视频高清在线观看| 狂野欧美激情性xxxx| 成年人午夜在线观看视频| 亚洲国产看品久久| 最近最新免费中文字幕在线| 黄色a级毛片大全视频| 欧美+亚洲+日韩+国产| 狠狠婷婷综合久久久久久88av| 亚洲欧洲精品一区二区精品久久久| 91成人精品电影| 69精品国产乱码久久久| 亚洲三区欧美一区| 久久久久国内视频| 国产亚洲av高清不卡| 国产伦人伦偷精品视频| 亚洲欧美一区二区三区久久| 日韩视频一区二区在线观看| 亚洲成人手机| 国产不卡av网站在线观看| 操出白浆在线播放| 午夜福利视频精品| 欧美成人免费av一区二区三区 | 日本vs欧美在线观看视频| 777久久人妻少妇嫩草av网站| 国产欧美日韩一区二区三| 操出白浆在线播放| 咕卡用的链子| 天天影视国产精品| av有码第一页| 精品欧美一区二区三区在线| 国产日韩欧美视频二区| 成人免费观看视频高清| 亚洲一卡2卡3卡4卡5卡精品中文| 日日摸夜夜添夜夜添小说| 国产精品香港三级国产av潘金莲| 欧美日韩成人在线一区二区| 啦啦啦视频在线资源免费观看| 日韩成人在线观看一区二区三区| 亚洲精华国产精华精| 亚洲精品久久午夜乱码| 亚洲 国产 在线| 免费观看av网站的网址| 精品国内亚洲2022精品成人 | 国产一区二区激情短视频| 久久国产精品人妻蜜桃| 男女床上黄色一级片免费看| 青青草视频在线视频观看| 国产欧美日韩一区二区三| 免费久久久久久久精品成人欧美视频| 日韩视频在线欧美| 老司机福利观看| 色在线成人网| 成人国语在线视频| 久久国产精品男人的天堂亚洲| 女人精品久久久久毛片| 国产在线观看jvid| 中文字幕色久视频| 99re在线观看精品视频| 成人国语在线视频| 极品少妇高潮喷水抽搐| 十分钟在线观看高清视频www| 老熟女久久久| 久久久精品国产亚洲av高清涩受| 免费看a级黄色片| 午夜老司机福利片| 97在线人人人人妻| 欧美黄色淫秽网站| 亚洲欧美色中文字幕在线| 在线av久久热| 日韩成人在线观看一区二区三区| 人妻久久中文字幕网| 午夜福利视频在线观看免费| 夜夜夜夜夜久久久久| 乱人伦中国视频| 亚洲一卡2卡3卡4卡5卡精品中文| 女人被躁到高潮嗷嗷叫费观| 精品一区二区三区四区五区乱码| 免费看a级黄色片| 女人精品久久久久毛片| 欧美精品高潮呻吟av久久| 中文字幕色久视频| 国产不卡一卡二| 欧美成人免费av一区二区三区 | 午夜福利,免费看| www.熟女人妻精品国产| 久久中文看片网| 亚洲情色 制服丝袜| 亚洲精品一二三| 国产日韩一区二区三区精品不卡| 99精品欧美一区二区三区四区| 黄片播放在线免费| 人妻一区二区av| 欧美久久黑人一区二区| 亚洲精品粉嫩美女一区| 午夜福利视频精品| 国产一区二区三区在线臀色熟女 | 精品国产乱子伦一区二区三区| 精品乱码久久久久久99久播| 色视频在线一区二区三区| 国产一区二区三区综合在线观看| 99精品在免费线老司机午夜| 老司机影院毛片| 99久久99久久久精品蜜桃| 亚洲九九香蕉| 在线 av 中文字幕| 久久免费观看电影| 久久99热这里只频精品6学生| 久久av网站| 久久99热这里只频精品6学生| 国产一区二区三区综合在线观看| 中文字幕av电影在线播放| 女性生殖器流出的白浆| 欧美日韩国产mv在线观看视频| 肉色欧美久久久久久久蜜桃| a级毛片黄视频| 亚洲精华国产精华精| 欧美日韩精品网址| 国产色视频综合| 每晚都被弄得嗷嗷叫到高潮| 日本vs欧美在线观看视频| 免费黄频网站在线观看国产| 午夜成年电影在线免费观看| 日本精品一区二区三区蜜桃| 精品高清国产在线一区| 日韩欧美免费精品| 国产精品1区2区在线观看. | 久久这里只有精品19| 91成人精品电影| 久久久国产一区二区| 国产淫语在线视频| 国产高清videossex| 午夜福利乱码中文字幕| 久久午夜综合久久蜜桃| 国产三级黄色录像| 一本久久精品| 久久天躁狠狠躁夜夜2o2o| 黄色成人免费大全| 久久久国产成人免费| 18禁观看日本| 夜夜夜夜夜久久久久| 99re在线观看精品视频| 欧美中文综合在线视频| 久久av网站| 欧美久久黑人一区二区| 国产免费av片在线观看野外av| 欧美黄色淫秽网站| av在线播放免费不卡| 少妇的丰满在线观看| 亚洲av成人不卡在线观看播放网| 18禁黄网站禁片午夜丰满| 乱人伦中国视频| 精品欧美一区二区三区在线| 他把我摸到了高潮在线观看 | 久久中文字幕人妻熟女| 亚洲午夜精品一区,二区,三区| 99riav亚洲国产免费| 999久久久精品免费观看国产| 免费不卡黄色视频| 在线观看一区二区三区激情| 高清在线国产一区| 国产欧美日韩精品亚洲av| 亚洲午夜理论影院| 国产福利在线免费观看视频| 在线观看www视频免费| 久久久国产成人免费| 久久精品国产a三级三级三级| 老司机午夜福利在线观看视频 | 色婷婷av一区二区三区视频| 久久亚洲精品不卡| 国产在线精品亚洲第一网站| 黑丝袜美女国产一区| 免费日韩欧美在线观看| 桃红色精品国产亚洲av| 亚洲专区中文字幕在线| 大码成人一级视频| 久久天堂一区二区三区四区| 国产成人免费无遮挡视频| 叶爱在线成人免费视频播放| 侵犯人妻中文字幕一二三四区| 两性夫妻黄色片| 亚洲 国产 在线| 91麻豆av在线| 国产伦理片在线播放av一区| 黑人操中国人逼视频| 最近最新中文字幕大全免费视频| 另类精品久久| 国产精品秋霞免费鲁丝片| 久久国产精品男人的天堂亚洲| 一级毛片电影观看| 亚洲,欧美精品.| 91精品国产国语对白视频| 国内毛片毛片毛片毛片毛片| 国产视频一区二区在线看| 中文字幕人妻熟女乱码| 国产又爽黄色视频| 国产一区二区 视频在线| 亚洲av电影在线进入| 女性被躁到高潮视频| 免费人妻精品一区二区三区视频| 日韩制服丝袜自拍偷拍| 自拍欧美九色日韩亚洲蝌蚪91| 国产伦理片在线播放av一区| 精品卡一卡二卡四卡免费| 天天操日日干夜夜撸| 国产精品久久久久久精品古装| 极品人妻少妇av视频| 久久精品熟女亚洲av麻豆精品| 黄色怎么调成土黄色| 精品久久久久久久毛片微露脸| 美女国产高潮福利片在线看| 国产av又大| 少妇精品久久久久久久| 制服人妻中文乱码| 俄罗斯特黄特色一大片| 久久国产精品影院| 亚洲国产av新网站| 国产高清激情床上av| 国精品久久久久久国模美| 中文字幕另类日韩欧美亚洲嫩草| 女同久久另类99精品国产91| 一区二区av电影网| 在线永久观看黄色视频| 啦啦啦视频在线资源免费观看| 久久99一区二区三区| 亚洲三区欧美一区| 亚洲国产av影院在线观看| 久热这里只有精品99| xxxhd国产人妻xxx| 国产福利在线免费观看视频| 蜜桃国产av成人99| 久久久久视频综合| 高清视频免费观看一区二区| 欧美 亚洲 国产 日韩一| 久久精品成人免费网站| 亚洲精品在线观看二区| 久久久精品94久久精品| 91成年电影在线观看| 9热在线视频观看99| 欧美日韩福利视频一区二区| 日韩 欧美 亚洲 中文字幕| 看免费av毛片| 久久精品国产亚洲av高清一级| 蜜桃国产av成人99| 亚洲第一欧美日韩一区二区三区 | 色尼玛亚洲综合影院| 亚洲精品国产精品久久久不卡| 免费日韩欧美在线观看| 国产av又大| 美女午夜性视频免费| 午夜免费鲁丝| 大型av网站在线播放| 免费在线观看影片大全网站| 高清av免费在线| 99在线人妻在线中文字幕 | 久热这里只有精品99| 黄色毛片三级朝国网站| 亚洲国产欧美一区二区综合| 狠狠精品人妻久久久久久综合| 国产精品久久电影中文字幕 | 男女下面插进去视频免费观看| 亚洲精品中文字幕在线视频| 精品久久蜜臀av无| 欧美在线黄色| 国产精品久久电影中文字幕 | 免费观看av网站的网址| 窝窝影院91人妻| 女性生殖器流出的白浆| 咕卡用的链子| 国产成人精品无人区| 9热在线视频观看99| 日韩大码丰满熟妇| 自线自在国产av| 桃花免费在线播放| 日韩视频在线欧美| 丰满人妻熟妇乱又伦精品不卡| 两个人免费观看高清视频| 黄色a级毛片大全视频| 久久久国产精品麻豆| 91九色精品人成在线观看| 99国产综合亚洲精品| 欧美 亚洲 国产 日韩一| 69av精品久久久久久 | 国产伦人伦偷精品视频| 一边摸一边抽搐一进一出视频| 夜夜爽天天搞| 久热爱精品视频在线9| 国产精品成人在线| 每晚都被弄得嗷嗷叫到高潮| 免费在线观看黄色视频的| 女人高潮潮喷娇喘18禁视频| 可以免费在线观看a视频的电影网站| 精品一区二区三卡| 成人精品一区二区免费| a在线观看视频网站| 丝袜喷水一区| 怎么达到女性高潮| 欧美在线一区亚洲| 久久精品国产a三级三级三级| 国产一区二区 视频在线| 黄频高清免费视频| 老熟妇仑乱视频hdxx| 伦理电影免费视频| 欧美激情高清一区二区三区| 久久久久久久大尺度免费视频| 黑人猛操日本美女一级片| 久久国产精品大桥未久av| 欧美日韩黄片免| 国产亚洲午夜精品一区二区久久| 日韩视频一区二区在线观看| 婷婷成人精品国产| 成人18禁在线播放| 久久中文字幕一级| 国产高清视频在线播放一区| 中文字幕最新亚洲高清| 男女免费视频国产| 亚洲国产成人一精品久久久| 搡老乐熟女国产| 国产成人精品无人区| 久久久国产欧美日韩av| 91老司机精品| 久久免费观看电影| 亚洲欧美色中文字幕在线| 精品高清国产在线一区| 免费看十八禁软件|