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

    Improving the electrochemical stability of AZ31 Mg alloy in a 3.5wt.%NaCl solution via the surface functionalization of plasma electrolytic oxidation coating

    2022-07-13 03:09:54MosKseemTehseenZehrBurkDikiiAliDfliHeWoongYngYoungGunKo
    Journal of Magnesium and Alloys 2022年5期

    Mos Kseem, Tehseen Zehr, Burk Dikii, Ali Dfli, He Woong Yng,Young Gun Ko

    a Department of Nanotechnology and Advanced Materials Engineering, Sejong University, Seoul 05006, Republic of Korea

    b Department of Metallurgical and Materials Engineering, Ataturk University, Erzurum 25240, Turkey

    c Faculty of Sciences Oujda, Department of Chemistry, Mohamed the First University, Oujda, Morocco

    d Pohang Institute of Metal Industry Advancement, Pohang 37666, Republic of Korea

    e Materials Electrochemistry Group, School of Materials Science and Engineering, Yeungnam University, Gyeongsan 38541, Republic of Korea

    Abstract

    Keywords: Hybrid composite; Coordination bond; Albumin; Fatty acid; Porous layer; Corrosion; Theoretical approach.

    1.Introduction

    In recent years,growing interests have been poured into the development of magnesium and its alloys in the automobile and electronic fields due to their lightweight, high strength,and good machinability [1,2].However, the high corrosion rate of Mg and its alloys would limit the industrial applications of these materials [3].Therefore, alloying [4]and surface treatments methods namely chemical conversion coatings, sol-gel, organic/polymer coatings, chemical/physical vapor deposition, anodizing, and plasma electrolytic oxidation(PEO) were suggested by several research groups to improve the protection properties of Mg-based materials in the corrosive environments [5–10].Among them, PEO has been considered as an environmentally friendly treatment method that leads to form protective oxide layers on the metallic substrates.Briefly, the high voltage applied in an alkaline electrolyte could create numerous micro-sized plasma discharges through a complex process involving thermal, electrochemical, and plasma reactions [11–18].Since PEO is a surfaceconversion coating technique, the surface of the metallic alloy is transformed into ceramicviaelectrochemical plasma reactions without damaging the metal alloys at its core.Structures with complex geometries can be coated without distorting the bulk microstructure because the whole substrate does not go through an elevated temperature during PEO coating [4,12].The plasma discharge activity during PEO would create a porous structure allowing the corrosive species to attack the underlying substrate.Thus, additional treatment would be required to delay or suppress the movement of corrosive species through the defective layer (IL) producedviaPEO [12].

    The application of hydrophobic coatings was proven as an effective approach in inhibiting the corrosion of Mg and its alloys [19].For example, highly protective hydrophobic coatings can be obtained by using fatty acids that bind to the metallic oxide surfaces through hard acid-hard base interactions in which the metallic cations from the oxide layer acts as hard acid while carboxylate anions from fatty acid behave like hard base [20,21].Moreover, the long aliphatic chain in these acids would help the development of a highly protective hydrophobic coating [22].According to Cui et al.[23], the immersion of PEO-coated AZ31 Mg alloy in an ethanolic solution of octadecanoic acid for 10 h would lead to sealing the majority of structural defects of PEO coating where a hydrophobic layer with a contact angle ~151 ° was obtained through the bonding of CH3(CH2)16COO-anions to the surface of PEO coating which resulted in significant improvements in the protective properties of the metallic alloy.However,unwanted side reactions would occur during the forming of such hydrophobic coatings because the high electrochemical activity of Mg alloys with many acids, alkalis,and saline solutions would lead to unwanted side reactions[24].Therefore, the hydrophobic coatings on Mg alloys are not popular when compared to other inert metals, such as Cu,Al, and Ti whose tendency to form passive films on their surfaces [25–27].Moreover, the chemical activity of Mg alloys is heterogeneous due to the presence of more than two phases with the structure of these materials [18].

    On the other hand, novel hybrid materials with various structures could be discovered by hybridization by exploiting the tendency of proteins to interact with a variety of nanomaterials [28].Thanks to its availability and nontoxicity, albumin (ALB) was utilized to tailor the surface characteristics of IL prepared on Mg alloys.For instance, Kaseem and Ko[29]reported the growth of flowers and rods-like structures on the IL as a result of immersion of IL in solutions containing different concentrations of albumin (ALB) and sodium phosphate which led to forming novel protective composites on the AZ31 Mg alloy.The ALB molecules were used to modify the layered double hydroxide (LDHs) films grown on IL for smart corrosion protection of AZ31 Mg alloy [30].Although the structural defects of IL can be disadvantages for the corrosion property, they have potential benefits.For example, the large surface area provided by the defect structure would facilitate the self-assembly of organic and inorganic compounds where the local anchoring of these compounds with simple structures around these defects was reported [31].

    Based on the discussion above, an issue can be raised of whether it is possible to exploit the competitive interactions between fatty acids and proteins to fabricate an organic layer(OL) on the top of the defective IL to improve the electrochemical stability of the AZ31 Mg alloy substrate.Here, we reported for the first time the formation of novel organic structures on top of the porous surface of IL to fabricate a novel hybrid composite composed of HA-ALB-IL as an individual layer on the surface of AZ31 Mg alloy substrate.The hybrid composite is fabricated by immersing the IL made on Mg alloyviaPEO in an ethanolic solution containing ALB and hexadecanoic acid (HA).So far, the use of such a defective layer as an immobilization platform for hybrid organic-inorganic flowery structures remains undocumented when two organic compounds were used.The porous surface of IL would play a critical role in providing nucleation sites for self-assembly of ALB-HA which will drive the self-assembly of flowery-like structure.The primary bonding mechanism responsible for flowery structure is proposed based upon DFT calculations of HA and ALB.The effectiveness of the prepared hybrid materials in protecting Mg alloy substrate against corrosion attack was studied and compared with other methods available in the literature.

    2.Experimental procedures

    2.1.Fabrication and functionalization of the defective inorganic layer

    The material used as a substrate in this study was AZ31 magnesium alloy sheets (30 × 25 × 4 mm) with a chemical composition (in wt.%) of 2.89 Al, 0.96 Zn, 0.31 Mn,0.15 Fe, 0.12 Si, and balance Mg.The AZ31 Mg alloy samples were treated in an electrolyte comprising of 4 g/L KOH,2 g/L WO3, and 10 g/L of Na3PO4.The particle size and zeta potential of WO3were measured to be ~300 nm and~-39 mV, respectively.Both values were determined by Zetapotentiometer(Nano ZS, ZEN360).The obtained inorganic layer is called IL.The processing parameters, namely oxidation time, current density, and frequency were controlled to be 600 s, 100 mA/cm2, and 60 Hz, respectively.To form an OL-IL hybrid composite, the IL samples were functionalized by immersion them in an ethanolic solution (70% ethanol)containing a mixture of 3 g/L albumin (ALB) and (2.5 g/L)hexadecanoic acid (HA) under continuous stirring for 10 and 24 h at 60 °C.The resultant samples were called IL-OL-10 h and IL-OL-24 h, respectively.Based on the preliminary experiments, the immersion periods were selected to be 10 and 24 h because the immersion times below 12 h did not show significant improvement in the corrosion resistance of IL.Whereas the immersion times higher than 24 h caused inferior improvements in the corrosion resistance of the IL sample as compared to the IL sample treated at 24 h.So,we thought the immersion of the IL sample for 24 h in the ALB/HA solution was enough to modify its porous structure.

    2.2.Microstructural and compositional analysis

    The structural morphologies of the examined samples were observed with scanning electron microscopy (SEMHITACHI, PT-S2800) equipped with dispersive X-ray spectroscopy (EDS).The thickness of the obtained coating was acquired from>10 different points of at least 5 different cross-sectional morphologies of the samples.The phase compositional analysis of the samples was carried out using X-ray diffraction (XRD; RIGAKU, D-MAX 2500) and Xray photoelectron spectroscopy (XPS, VG Microtech, ESCA 2000).The Scan range during XRD measurements was from 10 ° to 90 °.The functional groups that appeared due to the post-treatment were analyzed with FTIR-ATR spectroscopy(PerkinElmer Spectrum 100) from a wave-number range of 400–4000 cm-1.

    2.3.Electrochemical analysis

    Electrochemical measurements were conducted in an aerated 3.5 wt.%NaCl solution with a potentiostat(Gamry Interface 1000) in a three-electrode configuration, where platinum and Ag/AgCl were set as counter and reference electrodes,respectively.Before potentiodynamic polarization (PDP) measurements, the samples were left in the corrosive solution for 5 h to stabilize the open circuit potential (OCP).PDP curves were scanned from - 0.2 to + 0.4 V versus the OCP at a scanning rate of 1 mV?s-1.The Electrochemical impedance spectroscopy (EIS) measurements were performed with AC amplitude of 10 mV over the OCP in the frequency range from 100 kHz to 0.01 Hz after immersion periods of 5, 12,24, 96, 168, and 336 in the same solution.All electrochemical tests were performed at least four times to ensure data repeatability.

    2.4.Theoretical calculations

    Calculations were performed under the principles of density functional theory (DFT) and the geometry of hexadecanoic acid (HA) and albumin (ALB) molecules was optimizedviathe B3LYP/6–311G(d,p) basis set using the Gaussian (09) software [32–34].The reactivity of ALB and HA molecules was examined by calculations of the energy of the highest occupied molecular orbital (EHOMO) as well as the energy of the highest occupied molecular orbital (EHOMO).The values of energy gap(ΔEgap), dipole moment (μ), ionization potential (I), electron affinity (A), Mulliken’s electronegativity (χ), absolute hardness (η) and absolute softness (S) were estimated by following equations [31]:

    The fraction of emitted electrons (ΔN) was calculated using the above formula [31]:

    Here,χMgandηMgrepresent the theoretical value of Mulliken’s electronegativity and hardness value of Mg, respectively, whileχmolandηmolare the calculated electronegativities and hardness values of organic molecule respectively.The theoretical values ofχMgandηMgfor Mg are found to be 3.75 eV mol-1and 3.90 eV mol-1, respectively [35].

    3.Results and discussion

    3.1..Microstructure of hybrid composites

    Fig.1(a–c) shows the surface morphologies of IL, IL-OL-10 h, and IL-OL-24 h samples before and after the posttreatment using an organic solution containing HA and ALB.As shown in Fig.1a,the surface of IL reveals micropores with a multi-level layout.The hyper-level pores show an irregular,pit-like shape, while their hypo-level counterparts show a circular shape.These multi-level micropores arise due to the plasma discharge activity during PEO [12, 36–38].Cracks highlighted by arrows in Fig.1a are attributed to the residual thermal stress generated within the IL as a result of the rapid solidification of molten materials by the cold electrolyte [39].In contrast,multi-level pores observed on the surface of the IL sample became less dominant on the surface of IL-OL-10 h and IL-OL-24 h samples.After the post-treatment process in a solution of HA and ALB, the nucleation of flowery flakes on top of the IL surface was started as shown in Fig.1b for the IL-OL-10 h sample.The large surface area provided by IL due to the presence of micropores would responsible for such nucleation of flowery flakes like agglomerates of HAALB (OL) on the surface of IL.Accordingly, the structural defects that exist on the IL surface before the post-treatment are sealed due to the growth flowery-flake structure as shown in Fig.1 (b–e).Such effective sealing would be attributed to the unique interactions between the ALB and HA on the surface of IL which lead to an anisotropic growth of primary clusters by self-assembly [12,13,21].The active adsorption of OL (HA-ALB) operating the growth process is affected by the porous structure of IL and the bonding compatibility between IL and OL by electron donor-acceptor reactions, as we will be approved later by DFT calculations.Thus, ALB and/or HA molecules would provide electrons to the solution and the porous IL surface leading to the formation of the hybrid organic-inorganic film.EDS mappings on the surface of OL-IL-10 h indicated the existence of Mg, O, P, W, C,and N elements which tended to distribute homogeneously throughout the surface of the hybrid film.The presence of C and N would confirm the active adsorption of the OL on the surface of IL.The EDS elemental composition on the surface of IL, IL-OL-10 h, and IL-OL-24 h samples (Table 1) shows the contents of the elements that appeared in EDS mappings shown in Fig.1(f–n).

    Table 1EDS analysis on the surface of IL, IL-OL-10 h, and IL-OL-24 h samples.

    Fig.1.(a–e) SEM images showing the surface morphologies of the samples (a) IL, (b and d) IL-OL-10 h, (c and e) IL-OL-24 h.(f–n) EDS mapping analysis on the surface of IL-OL-24 h sample shown in Fig.1c where the Mg, O, P, W, C, and N elements were detected.

    Fig.2 (a–c) displays the corresponding cross-sectional images of the IL,IL-OL-10 h,and IL-OL-24 h.The thickness of the hybrid composites was found to be ~23.2 and ~24.3 μm for IL-OL-10 h and IL-OL-24 h samples, respectively while the thickness of IL was calculated to be ~20.1 μm.As indicated by yellow circles, the formation of flowery flakes led to the sealing of micro-pores in IL-OL-10 h and IL-OL-24 hsamples.The EDS elemental mappings, shown in Fig.2 (d–j), confirmed the existence of Mg, O, P, W, N, and C in the sealed porous structure.

    Fig.2.(a–c) SEM images showing the cross-sectional morphologies of (a) IL, (b) IL-OL-10 h, (c) IL-OL-24 h.(d–j) The corresponding EDS mapping of the IL-OL-24 h sample showing the presence of Mg, O, P, W, C, and N elements (For interpretation of the references to color in this figure, the reader is referred to the web version of this article).

    3.2.Compositional analysis of hybrid composites

    Fig.3(a) displays the FT-IR spectra of IL-OL-10 h, and Il-OL-24 h to prove the interaction of PA and ALB with IL.The peaks evident at ~1520 cm-1and ~1573 cm-1are attributed to –COO-asymmetric and symmetric stretching vibrations,respectively; though according to Liu et al.[22], the symmetric stretching peak of–COO-should be at ~1468 cm-1which suggested that a combination of different interactions would occur on the surface of IL.Similar findings were reported by Cui et al.[23]who examined the sealing effect of AZ31 Mg alloy by stearic acid.The peaks at ~1022 cm-1would be attributed to the ethanol adsorbed on the surface [23].The peaks at about 2867 and 2943 cm-1were contributed to the C–H asymmetric and symmetric stretching vibration, respectively [40].The position of the N–CO–NH2group in ALB is commonly used to discover the conformational changes on the ALB molecule during the adsorption on the metallic sur-face [41,42].In this work, the N–CO–NH2group vibrations were noticed at 1686 and 1645 cm-1which is slightly different from the values of the amide group reported of pure ALB which usually occurs at 1629–1748 cm-1[41].This change in band position suggested that interaction with ALB with IL and HA surface would induce strong changes in ALB conformational state in favor of less-ordered conformations.The band at 1520 cm-1would be related to an interaction between NH2and positively charged sites (Mg sites) on the IL surface[21].

    Fig.3.(a) FTIR spectra of IL, IL-OL-10 h, and IL-OL-24 h samples.(b) XRD analysis of the IL, IL-OL-10 h, and IL-OL-24 h samples.(c–f) The high-resolution scan for (c) Mg 1 s, (d) W 4f, (e) C 1 s, and (f) N 1 s for IL-OL-24 h sample.

    The compositional analyses of IL,IL-OL-10 h,and IL-OL-24 h samples were investigated by XRD, as shown in Fig.3b.MgO, Mg3(PO4)2, WO3, W18O49, and Mg were detected in the IL sample.MgO can be formed by the Mg2+generated from AZ31 Mg alloy substrate and O2–anions produced from the electrolyte during PEO (Eq.(7)):

    PO43-anions produced from the ionization of Na3PO4are moved toward the anode under the effect of the high electric field and react with Mg2+cations in the vicinity of the anode(Eq.(8)):

    The peaks of Mg were detected due to the penetration of X-rays into the AZ31 Mg alloy substrate through the porous thin layer.The presence of WO3and W18O49is ascribed to the successful incorporation of the WO3additive, and its transformation to W18O49under plasma conditions[31].Since WO3nanoparticles possessed a negative zeta potential, their electrophoretic movement towards the anode would ensure the delivery of the particles at the anode-electrolyte interface [11,12].Thus, the incorporation of these nanoparticles during the PEO process could occur by the electrophoresis process, then the high local temperature of plasma discharges melted and sintered the WO3nanoparticles into the coating[12].Here, two assumptions were considered to be important during the incorporation behavior of WO3nanoparticlesviathe PEO process.First, the extent of oxygen evolution during the initial stages of the PEO process could be helpful to determine the preservation of these particles that are already embedded in the solution.Second, the breakdown phenomena of the coating and generation micro-pores and cracks could decide the migration of WO3nanoparticles at the later stages of the PEO process.The high plasma temperature during PEO would facilitate the phase transformation of WO3to W18O49[31].On the other hand, the XRD patterns of OL-IL-10 h and OL-IL-24 h samples showed the additional peaks at 20.01, 20.76, 21.50, 23.57, and 39.41° along with their observed peaks of IL,suggesting the successful formation of OL on the top of IL.The complete growth of the flowery flakes which seal the entire porous structure in the case of the ILOL-24 h sample would explain the high intensity of these peaks in comparison to the case of the IL-OL-10 h sample.

    The chemical states of the flowery-flake structures are investigated in detail by XPS analysis.Fig.3 (c–f) shows the high-resolution spectra of Mg 1s, W 4f, C 1s, and N 1s for the IL-OL-24 h sample.The deconvolution of the Mg 1s spectrum led to three peaks at binding energies of ~1303.28 eV,~1303.88 eV, and ~1304.28 eV corresponded to Mg, MgO,and Mg3(PO4)2respectively [31].The spectrum W 4f was decomposed into peaks of W 4f7/2,and W 4f5/2for IL,where thebinding energy at ~34.38 eV and ~36.5 eV associated with W18O49, and the binding energies at 35.48 eV, and 37.75 eV corresponded to WO3[31].The deconvolution of the N 1s spectrum shows three deconvoluted peaks at ~398.4 eV representing C–N bond,and at ~399.2 eV representing N–C=O group, and at ~399.8 eV describing the N-Mg bond, suggesting the binding of Mg with ALB to form coordinating complexes [31].As such, the deconvolution of C 1s spectrum showed five peaks at ~284 eV, ~285.7 eV, ~287.3 eV,~289.4 eV, and ~290.4 eV related to (C–C/C–H), (C–N),and (N–C = O), (O–C = O), andπ-πrespectively [19, 43–45].The detection of the C 1s and N 1s peaks, arising from the post-treatment in a solution composed of ALB and HA,would suggest the adsorption of ALB and HA on the surface of IL, forming flower-like organic structures.

    Fig.4.(a) PDP curves scanned from - 0.25 to + 0.40 V concerning the OCP at a scan rate of 1 mV/s for AZ31 Mg alloy, IL, IL-OL-10 h, and IL-OL-24 h samples.(b) Ecorr vs.icorr for AZ31 Mg alloy, IL, IL-OL-10 h, and IL-OL-24 h samples.

    3.3.Electrochemical performance of hybrid composites

    Fig.4a shows PDP curves of the AZ31 Mg alloy substrate,IL, IL-OL-10 h, and IL-OL-24 h samples with a direct comparison with the IL before the post-treatments in a solution of ALB and HA.The parameters of corrosion current density (icorr) and corrosion potential extracted from PDP curves are plotted in Fig.4b.As shown in Fig.4b, the AZ31 Mg alloy substrate had a higher corrosion susceptibility as compared to other samples which was reflected by the largest value oficorrand the most negative value ofEcorr.Although the PEO coating (IL sample) led to improve the corrosion resistance of the AZ31 Mg alloy, the formation of flowery structures on the surface of IL-OL-10 h and IL-OL-24 h samples resulted in significant changes in terms of theEcorrandicorrwhich correlate with the thermodynamics and kinetics of the corrosion reactions, respectively.As shown in Fig.4a,the PDP curve of the IL sample was shifted to the nobler potential after deposition of the OL layer where the most positive increase in the potential was observed in the IL-OL-24 h sample.This is attributed to the sealing of structural defects of IL by the OL, limiting the activity of the whole surface.Therefore, its surface becomes less prone to chemical processes, as indicated by the noble shift ofEcorr.The PDP curves of IL-OL-10 h and IL-OL-24 h samples showed a shift to the lower region of theicorr, which indicates slower anodic reactions.This result suggests that the flowery-flakes structure is an effective approach to improve the stability of IL by blocking the porous structure and improving the barrier properties.In other words, the flowery-like structures can be considered as a potential candidate to improve the durability of anti-corrosion layers by sealing the structural defects and strengthening the barrier properties.This would,in turn,block the additional inward flow of the corrosive ions (Cl-) into the substrate.Moreover, such structures would provide additional shielding that prevented the flow of the corrosive ions towards the multi-level micropores.Therefore, it is believed that the high stability in the corrosive solution of the IL-OL-24 h sample is attributed to the complete and uniform sealing of micro-pores due to the deposition of OL acting as a barrier film.

    Fig.5(a–c) shows the cross-sectional images of IL, ILOL-10 h, and IL-OL-24 h samples after 5 h immersion in a 3.5 wt.% NaCl solution.It can be seen from Fig.5a that IL had suffered severe corrosion since many corroded areas were found on the surface of this sample.In contrast, the flake-like structure in IL-OL-10 h and IL-OL-24 h samples film protected the Mg alloy substrate to some extent where some degraded areas are rarely observed, as shown in Fig.5 (b and c).Interestingly, the cross-sectional image of the IL-OL-24 h sample (Fig.5c) was almost unaffected by the corrosion, confirming that the effectiveness of the flower-flakes structure in improving the protection properties of AZ31 Mg alloy substrate.

    To approve the effectiveness of the approach proposed in the present work, theicorrandEcorrvalues derived from PDP tests were compared with data in the literature as displayedin Fig.5d.For the examples shown in Fig.5d, the immersion times before PDP tests were varied from 10 min to 4 h.As it can be seen from Fig.5d, the formation of flowery flakes structures on the surface of IL, through the unique interactions between ALB and HA molecules,is an effective method to improve the protection properties of Mg substrates in comparison to different PEO coatings, conversion coatings, LDH coating,anodizing, hybrid coating,and anodizing+LDH coatings [46–59].Although the value ofEcorrof the IL-OL-24 h sample was more negative than the other cases, the smallest value oficorr(2.47 × 10-11A/cm2) for this sample would imply a superior corrosion resistance of the AZ31 Mg alloy in comparison to the literature data.This was attributed to the stable barrier protection achieved by flowery flakes structures formed due to the unique interactions between ALB and HA on the surface of IL.Even a longer immersion period of 5 h was used in the present work, better corrosion resistance was obtained in comparison to other studies shown in Fig.5d in which the samples were immersed in the corrosive solutions for shorter periods (from10 min to 4 h).Since all materials employed in the present works were almost less toxic and eco-friendly, the present approach can be applied to fabricate functional hybrid composites with potential biomedical applications.

    Fig.5.(a–c) Morphological features of the samples after 5 h immersion in 3.5 wt.% NaCl solution where (a) IL sample, (b) IL-OL-10 h sample, and (c)IL-OL-24 h sample.(d) Comparison of Ecorr and icorr values obtained in the present work with the counterparts values collected for Mg alloy substrates treated by couples of surface treatments [46–59].The immersion periods in the corrosive solutions were varied from 10 min to 5 h.

    3.4.EIS evaluation for long time immersion

    For the assessment of the corrosion performance of the coatings with prolonged exposition to 3.5 wt.% of NaCl solution, EIS measurements were performed at different immersion times up to 336 h.The results obtained for the IL, ILOL-10 h, IL-OL-24 h samples are presented in Fig.6, in the form of Nyquist and Bode plots.Regardless of the immersion time, the IL sample showed lower corrosion resistance than IL-OL-10 h and IL-OL-24 h samples as deduced from the smaller diameter of the capacitive loops in the Nyquist plots and the lower values of modulus of impedance at low frequencies.As the immersion time increase, IL exposes to thecorrosive media and form loose Mg(OH)2[60], which has no blocking effect on Cl-, thus facilitating a further degradation during 12–168 h immersion.This degradation became more pronounced when the immersion time reaches 336 h where the diameter of the capacitance loop of IL diminished significantly.The significant shrink of capacitive loops and drop of the total impedance of IL sample with increasing immersion time confirm the less protection of the coating and the active dissolution of the IL sample.The defects in the IL sample(micropores and microcracks) acted as the active sites, providing pathways for quick diffusion of corrosive ions to the substrate with a higher corrosion rate.With increasing immersion time up to 368 h, the large drop of corrosion resistance would imply that the substrate had been corroded seriously because of the insufficient protection.Although the diameters of the capacitive loops in the Nyquist plots for IL-OL-10 h and IL-OL-24 h samples were reduced with increasing the immersion times, the corrosion resistance of these samples was higher than that for the IL sample, which can be considered as sealing effect, and shows a slow-corrosion process during long-term immersion.

    Fig.6.EIS results obtained for the IL, IL-OL-10 h, and IL-OL-24 h samples at different immersion times, where (a–c) IL sample, (d-f) IL-OL-10 h sample,and (g–i) IL-OL-24 h sample.The solid lines represent the fitting curves obtained with the equivalent circuits.

    In general, the EIS spectrum in a broad frequency range may exhibit one, two, or more time constants depending on whether the coating has a single-layered or multi-layered structure, whether it is porous or compact, whether the spacecharge region that puts restrictions on the charge transfer forms at the potential used, or the charge transfer process is effected by diffusion [61].The impedance response for the three samples at different immersion periods indicates the presence of several relaxation processes, as revealed by phase plots in Fig.6 (c,f,i).At least two-time constants for IL, ILOL-10 h, and IL-OL-24 h samples can be identified from the phase angle plots shown in Fig.6(c,f,i).It is also evident that the IL-OL-10 h and IL-OL-24 h samples had higher phase angle values corresponding to the higher frequency (1 MHz) as compared to the IL, suggesting higher corrosion protective ability.Moreover, it can be seen some deviation of the phase maximum values from 90owhich indicates a non-ideal capacitor behavior.Lopez-Ortega et al.[62], attributed this dispersion to inhomogeneities in the dielectric material,i.e.passive layer, such as porosity, mass transport, and the relaxation effect.The time constants, however, are hardly observed due to a strong overlying [63,64].Similar findings were reported by Zhang et al.[65]and Kaseem et al.[13]for the PEO coating made on AZ31 Mg alloy followed by several types of post-treatments.The protective properties can be measured from the value of |Z|f→0in the modulus Bode plots [66,67].As shown in Fig.6 (b, e, h), the high values of|Z|0.1Hzcanclearly present the beneficial effect of immersion in an organic solution composed of ALB and HA for all immersion times.Interestingly, the values of |Z|0.1Hzfor the IL-OL-24 h sample were significantly higher than other cases which would confirm that the flowery flake structures could significantly suppress the movement of corrosion species towards metallic alloy substrate, which was consistent with the PDP results.

    Fig.7.(a) The equivalent circuit model used for fitting the EIS results of the IL sample at immersion periods up to 96 h.(b) The equivalent circuit model used for fitting the EIS results of IL sample at immersion periods of 168 and 336 h.(c) The equivalent circuit model used for fitting the EIS results of IL-OL-10 h, and IL-OL-24 h samples.

    To explore the different electrochemical processes that the hybrid composite films undergo, the impedance responses are fitted to equivalent circuit (EC) models shown in Fig.7 by considering the physical and chemical features.Due to the dual-layered structure of the IL sample, two constant phase elements, such asCPEoandCPEicorresponded to the outer layer and the inner layer, respectively were used.Here,Rodescribes the resistance of the outer layer whileRiindicates the resistance of the inner layer.Therefore, the EC model shown in Fig.7a was used for fitting EIS results of the IL sample at immersion times up to 96 h while the EC model in Fig.7b was used for fitting the results at longer immersion periods of 168 and 336 h.Here,CPEdlandRctdescribe the double-layer capacitance and charge transfer resistance at the Mg alloy/solution interface.As for IL-OL-10 h and IL-OL-24 h samples,the EC model in Fig.7c was used for fitting the EIS results.In this model,CPEfandRfrepresented the constant phase element and the resistance of the flowery flakes grown on the surface of IL.The pure capacitor was replaced with aCPEowing to the inhomogeneous surface of the coatings obtainedviaPEO [68].The chi-squared (χ2) error value was 10-3order of magnitude for all cases, implying that the EC models in Fig.7 were adequate for fitting Fig.6.The parameter values drawn from these equivalent circuit models are available in Tables 2–4.In all samples, the values ofRiwere higher than values ofRodue to the compact structure of the coatings producedviaPEO [12].In addition to the high values ofRoandRiin comparison to that of the IL sample, the higher protective ability provided by IL-OL-10 h and IL-OL-24 h samples would be reflected by the presence ofRf.It can be noted that the IL-OL-24 h samples showed the highest values of theRfandRofor all times of immersion.The evolution of theRifor the three samples after different immersion periods in a 3.5 wt.% NaCl solution revealed that the values ofRiwere almost identical for samples at short immersion periods (5 and 12 h) in a 3.5 wt.% NaCl solution.Although the values ofRishowed a sharp decrease with increasing immersion times, the values ofRifor IL-OL-24 h were almost unchanged by increasing the immersion time, indicating considerable corrosion protection of the inner layer of the IL sample was obtained by the deposition of OL which has flowery flake structure.

    Table 2Electrochemical parameters obtained from fitting IL sample.The unit of Y is (S· sn· cm-2) while the unit of R is (Ω cm2).

    Table 3Electrochemical parameters obtained from fitting IL-OL-10 h sample.The unit of Y is (S· sn· cm-2) while the unit of R is (Ω?cm2).

    Table 4Electrochemical parameters obtained from fitting IL-OL-24 h sample.The unit of Y is (S· sn· cm-2) while the unit of R is (Ω·cm2).

    Table 5Calculated descriptors of ALB and HA molecules using the DFT method.

    The physical meaning ofCPE’sparameters would be correlated with the morphology of the IL, IL-OL-10 h, and ILOL-24 h samples.The high value ofCPE-Yvalue for the IL sample would suggest that a larger surface area would be exposed to the 3.5 wt.% NaCl solution, facilitating the charge transfer of Cl-species towards the AZ31Mg alloy substrate[69].In contrast, the low values of theCPEffor the IL-OL-24 h sample would be an indicator of the effective suppression of the infiltration of corrosive ions due to the formation of flowery flakes structure, covering the entire porous surface of IL.On the other hand, the values ofCPE-nwould describe the uniformity of each sample where the high values ofCPE-nfin the IL-OL-24 h sample would reflect the high uniformity of this sample which led to reduce the exposed interfacial surface, impeding anodic dissolution of AZ31 Mg alloy [12,70].

    3.5.Formation mechanism of hybrid films

    DFT was utilized to evaluate the active sites responsible for the adsorption of ALB and HA molecules on the surface of IL.Figs.8 and 9 show the computed FMO structures of both ALB and HA molecules.The chemical reactivity of a reacting molecule would be described by the interaction between its HOMO and LUMO levels [71–74].The ability of the molecule to donate electrons to electrophilic metal orbitals is described by the high values ofEHOMOwhile the low value ofELUMOimplies the tendency of a molecule to accept electrons from stronger donating molecules[21].Accordingly,the high binding ability of the molecule and ions to the metal surface would be described by the higher values ofEHOMOand lower values ofELUMO[13,31,39].Based on the DFT results listed in Table 5, the ALB molecule is more reactive than the HA molecule as it possesses a higherEHOMOvalue of (-7.5094 eV) and a lowerELUMOvalue of (-1.0216 eV).This would also be described by the low value of the energy gapΔE= (ELUMO-EHOMO) for the ALB molecule (~6.49 eV) in comparison to that for HA molecule (~7.78 eV)[31].Besides, it is well known that the higher reactivity of a molecule would also be described by increasing the softness value (S) or by decreasing the hardnessηvalue [75].According to the results tabulated in Table 5, ALB would be rigorously reactive than HA.In addition, the interaction of the system metal-surface/organic molecules is subjected to the variation of their electronegativity [33].This difference in electronegativity trains generally the electron transfer from the organic molecules to the metal surface.The obtained absolute electronegativity values for both tested molecules show thatχof ALB is greater than that of HA.Such a result has a significant impact on-trend to attract electrons in the bonding formed between these molecules and the metal surface.This suggestion is in good agreement with the experimental results obtained above.As the high values ofΔNwould describe the stability of donor-acceptor, respecting the condition thatΔNis less than 3.6, both molecules show a high electron-donating ability, which leads to good coverage of the metal surface, confirming thus, the electron transfer from the tested organic molecules to the surface.Moreover, from Table 5, the ALB molecule had a higher electron-donating ability than HA which was attributed to the higher number of lone pairs and denticity in ALB due to the presence of thepeptide group.On the other hand,it is well-known that the adsorption capacity of organic molecules over the metal surface would affect the corrosion inhibition property of the organic molecules [76,77].In general, the adsorption phenomenon of organic molecules would rely on their physical and chemical characteristics due to several factors related to the functional groups, molecular structure, and electron density, present at the donor atom [78].First, the –COOH group present in HA would form a secondary bond with the IL surface in a process of physisorption [79].Second, chemisorption may also be involved for the better adsorption capacity where the dissociated –COO-anions form bonds with Mg atoms of the metal surface, leading to the formation of the magnesium carboxylate complex [79].ALB is involved in the anticorrosion activity as it carries more lone pairs which impede the dissolution of Mg and creating a uniform barrier layer on the IL surface.Here, competitive adsorption between ALB and HA molecules would occur on the surface of IL.Since ALB molecules have higher reactivity than the HA molecules, as approved by DFT calculations, the surface of IL would be occupied by ALB molecules which hindered to some extent the adsorption of HA molecules.ALB contains more functional groups than that in HA which indicates its strong capability to form a stable barrier monolayer on the IL surface through physical and chemical interactions.This would be occurred by directing the ALB molecules to be adsorbedviathe lone pairs of electrons of the nitrogen atom of the peptide group (CONH2), and the oxygen atom of the carbonyl group(-COO), which would form coordinate bonds with IL surface,protecting it from the corrosion[80].Moreover,with the presence of ALB and HA together, a strong barrier layer could be formed, and the corrosion protection increased because of the uniform, deposition of the OL, which may prevent water molecules and corrosive anions (Cl-) from penetrating the IL surface through its hydrophobic nature due to the presence of hydrophobic long-chain in HA.The possible bonding be-tween ALB and HA molecules would also contribute to form a stable OL on the IL, leading to form a stable barrier layer effectively protecting the AZ31 Mg alloy substrate from corrosion.It is interesting to recognize that the better protective properties observed in the case of OL-IL-24 h sample were ascribed to their polarity and the adsorption of OL on the IL through a lone pair of electrons of the (CONH2), and the oxygen atom of the carbonyl group (-CO) present in solution in higher amounts.

    Fig.8.Optimized structure, the electron density distribution of the HOMO and LUMO orbitals, and Mulliken atomic charge distribution of ALB.

    Fig.9.Optimized structure, the electron density distribution of the HOMO and LUMO orbitals, and Mulliken atomic charge distribution of HA.

    4.Conclusions

    In summary, a flowery-like multidimensional structured ALB-HA nanocomposite was successfully made on the defective surface of IL by the chemical doping in the mixed solvent of ethanol and water containing ALB and HA molecules.The flowery flake structures sealed the entire surface of IL helped to suppress the corrosion rate of AZ31 Mg alloy which was ascribed not only to reducing the porosity of IL but also to the strong orbital interactions between ALB and HA which resulted in the formation of a stable organic barrier layer on the surface of IL.Theoretical calculation shows efficient interactions between tested molecules (ALB & HA) and the metal surface thanks to the free-electrons pairs of functional groups existing in these molecules.The results of the theoretical study are in good agreement with those obtained experimentally.The current inexpensive method of producing flowery flakes like organic compounds with a large surface area could find potential catalytic and photocatalytic applications.

    Declaration of Competing Interest

    The authors declare that no conflict of interests.

    Acknowledgment

    This work was supported by the National Research Foundation of Korea (NRF) grant funded by the Korea government (MSIT) (No.NRF-2019R1G1A1099335).This work is supported also by the Mid-Level Researcher National Project of the National Research Foundation (NRF) funded by the Ministry of Science and ICT, Republic of Korea (NRF-2020R1A2C2004192).This research is supported partly by Basic Research Program through the National Research Foundation, Republic of Korea (NRF-2019R1FA1062702).

    亚洲av成人不卡在线观看播放网| 午夜两性在线视频| tocl精华| 真人一进一出gif抽搐免费| 两个人视频免费观看高清| 草草在线视频免费看| 老司机在亚洲福利影院| 国产精品乱码一区二三区的特点| 中文字幕av在线有码专区| 国产av麻豆久久久久久久| 欧美高清成人免费视频www| 看免费av毛片| 一级毛片高清免费大全| 国产精品98久久久久久宅男小说| eeuss影院久久| 淫秽高清视频在线观看| e午夜精品久久久久久久| 99精品在免费线老司机午夜| 日韩欧美 国产精品| 亚洲熟妇中文字幕五十中出| 黄色片一级片一级黄色片| 亚洲欧美激情综合另类| 色播亚洲综合网| 免费在线观看影片大全网站| 日本五十路高清| 精品一区二区三区av网在线观看| 中文字幕熟女人妻在线| 国产三级黄色录像| 日本 av在线| 最近视频中文字幕2019在线8| 18禁裸乳无遮挡免费网站照片| av片东京热男人的天堂| 中文字幕人妻熟人妻熟丝袜美 | 身体一侧抽搐| 精品久久久久久,| 国产精品免费一区二区三区在线| 变态另类丝袜制服| 老汉色av国产亚洲站长工具| www日本黄色视频网| 亚洲精品在线观看二区| 成人三级黄色视频| 18禁国产床啪视频网站| АⅤ资源中文在线天堂| 丰满乱子伦码专区| 极品教师在线免费播放| 综合色av麻豆| 亚洲av一区综合| 男人的好看免费观看在线视频| 最近在线观看免费完整版| 热99re8久久精品国产| 首页视频小说图片口味搜索| svipshipincom国产片| 淫妇啪啪啪对白视频| 亚洲无线观看免费| 亚洲 欧美 日韩 在线 免费| 日本免费一区二区三区高清不卡| 精品久久久久久,| 草草在线视频免费看| 天堂√8在线中文| 欧美av亚洲av综合av国产av| 亚洲国产中文字幕在线视频| 99久久成人亚洲精品观看| 可以在线观看毛片的网站| av天堂中文字幕网| 2021天堂中文幕一二区在线观| 日韩成人在线观看一区二区三区| 精品国产美女av久久久久小说| 88av欧美| 国产精品三级大全| 欧美激情久久久久久爽电影| 噜噜噜噜噜久久久久久91| 老汉色av国产亚洲站长工具| 久久久久九九精品影院| 国产v大片淫在线免费观看| 中文资源天堂在线| 12—13女人毛片做爰片一| 一级a爱片免费观看的视频| 亚洲国产欧美人成| 中亚洲国语对白在线视频| 亚洲一区高清亚洲精品| 午夜视频国产福利| 全区人妻精品视频| 亚洲成人中文字幕在线播放| 国产精品久久久久久亚洲av鲁大| 哪里可以看免费的av片| 国产高清视频在线观看网站| 久久精品国产亚洲av香蕉五月| 色av中文字幕| 日本 欧美在线| 99久久精品国产亚洲精品| 色综合婷婷激情| 午夜福利在线在线| 90打野战视频偷拍视频| 狂野欧美白嫩少妇大欣赏| 亚洲av中文字字幕乱码综合| 悠悠久久av| 麻豆国产97在线/欧美| 日韩有码中文字幕| 天堂av国产一区二区熟女人妻| 久久久久久九九精品二区国产| 床上黄色一级片| 成人av在线播放网站| 脱女人内裤的视频| 亚洲一区二区三区色噜噜| 五月玫瑰六月丁香| 国产精品嫩草影院av在线观看 | 香蕉av资源在线| 久久精品91无色码中文字幕| 国内久久婷婷六月综合欲色啪| 99久久九九国产精品国产免费| 精华霜和精华液先用哪个| 精品免费久久久久久久清纯| 亚洲精品456在线播放app | 精品午夜福利视频在线观看一区| 久久久久亚洲av毛片大全| 五月玫瑰六月丁香| 久久久色成人| 两人在一起打扑克的视频| 看黄色毛片网站| h日本视频在线播放| 人人妻人人澡欧美一区二区| 少妇裸体淫交视频免费看高清| 免费在线观看成人毛片| 欧美乱码精品一区二区三区| 男女午夜视频在线观看| 亚洲av第一区精品v没综合| 人妻久久中文字幕网| 99久久九九国产精品国产免费| 婷婷六月久久综合丁香| 桃红色精品国产亚洲av| 日韩大尺度精品在线看网址| 日韩欧美 国产精品| 久久久国产精品麻豆| 日本撒尿小便嘘嘘汇集6| 久久精品亚洲精品国产色婷小说| 国产成人系列免费观看| 99国产精品一区二区蜜桃av| 国内久久婷婷六月综合欲色啪| 午夜两性在线视频| 免费高清视频大片| 女警被强在线播放| 欧美日韩国产亚洲二区| 精华霜和精华液先用哪个| 国产精品 国内视频| 级片在线观看| 亚洲一区二区三区色噜噜| 亚洲国产欧美人成| 成年免费大片在线观看| 搞女人的毛片| 日本 av在线| 美女cb高潮喷水在线观看| 国产激情欧美一区二区| 久久久久亚洲av毛片大全| 男女之事视频高清在线观看| 欧美日韩亚洲国产一区二区在线观看| av天堂中文字幕网| 日本精品一区二区三区蜜桃| 日韩有码中文字幕| 香蕉久久夜色| 色噜噜av男人的天堂激情| 在线观看午夜福利视频| 欧美一级a爱片免费观看看| ponron亚洲| 国内毛片毛片毛片毛片毛片| 亚洲av电影不卡..在线观看| 波野结衣二区三区在线 | 99久久精品国产亚洲精品| 成人高潮视频无遮挡免费网站| 国产精品久久久久久精品电影| 久久午夜亚洲精品久久| 成人精品一区二区免费| 日韩欧美精品v在线| 国产高清有码在线观看视频| 日本三级黄在线观看| 国产精品久久久久久久电影 | 国产亚洲精品av在线| 美女大奶头视频| 精品久久久久久久末码| 免费观看精品视频网站| 日本五十路高清| 99国产精品一区二区蜜桃av| 九色国产91popny在线| 精品电影一区二区在线| 少妇裸体淫交视频免费看高清| 亚洲不卡免费看| 动漫黄色视频在线观看| 露出奶头的视频| 精品日产1卡2卡| 日本黄色视频三级网站网址| 床上黄色一级片| 欧美中文日本在线观看视频| 国语自产精品视频在线第100页| 久久中文看片网| 亚洲人成电影免费在线| 国产免费男女视频| 一进一出好大好爽视频| 99在线视频只有这里精品首页| 中文资源天堂在线| 草草在线视频免费看| 久9热在线精品视频| 在线观看舔阴道视频| 亚洲美女视频黄频| 1024手机看黄色片| 久久久久免费精品人妻一区二区| 国产不卡一卡二| 久久性视频一级片| 国产亚洲精品综合一区在线观看| 国产真人三级小视频在线观看| 性色avwww在线观看| 国产激情偷乱视频一区二区| 91av网一区二区| 18禁裸乳无遮挡免费网站照片| 精品国内亚洲2022精品成人| 久久人人精品亚洲av| 免费看十八禁软件| 欧美成人a在线观看| 一本精品99久久精品77| 91字幕亚洲| 亚洲欧美日韩卡通动漫| 国产精品一区二区免费欧美| 人人妻,人人澡人人爽秒播| 亚洲国产精品成人综合色| 一本精品99久久精品77| 日韩精品青青久久久久久| av视频在线观看入口| 欧美一级a爱片免费观看看| 久久人人精品亚洲av| 99riav亚洲国产免费| 亚洲国产欧美网| 国产精品香港三级国产av潘金莲| 欧美乱妇无乱码| 搞女人的毛片| 88av欧美| 亚洲成a人片在线一区二区| 久久精品国产综合久久久| 色噜噜av男人的天堂激情| 狠狠狠狠99中文字幕| 亚洲久久久久久中文字幕| 一本综合久久免费| 国产精品日韩av在线免费观看| 国产亚洲欧美在线一区二区| 在线免费观看不下载黄p国产 | 欧美性感艳星| av片东京热男人的天堂| 日韩精品中文字幕看吧| 九九久久精品国产亚洲av麻豆| 性色av乱码一区二区三区2| 日本三级黄在线观看| 身体一侧抽搐| 嫩草影院精品99| 成年免费大片在线观看| 一级毛片高清免费大全| 激情在线观看视频在线高清| 亚洲精品一区av在线观看| 亚洲av成人av| 亚洲欧美日韩高清专用| 黄色视频,在线免费观看| 欧美日本亚洲视频在线播放| 国产精品免费一区二区三区在线| www.www免费av| 99热这里只有精品一区| 婷婷精品国产亚洲av| 国产91精品成人一区二区三区| www.www免费av| 脱女人内裤的视频| 午夜福利在线观看免费完整高清在 | 国产中年淑女户外野战色| 日韩精品中文字幕看吧| av国产免费在线观看| 狠狠狠狠99中文字幕| 国产精品久久电影中文字幕| 亚洲av二区三区四区| 国产精品久久久久久精品电影| 国产精华一区二区三区| 一级毛片女人18水好多| 99久久久亚洲精品蜜臀av| 有码 亚洲区| 老司机午夜十八禁免费视频| 欧美乱色亚洲激情| 最新中文字幕久久久久| 亚洲精品久久国产高清桃花| 国产高清videossex| 日本成人三级电影网站| 免费电影在线观看免费观看| 一本一本综合久久| 成人欧美大片| 人人妻人人澡欧美一区二区| 成人午夜高清在线视频| 亚洲久久久久久中文字幕| 在线a可以看的网站| 亚洲精品色激情综合| 1024手机看黄色片| 麻豆一二三区av精品| 婷婷精品国产亚洲av在线| or卡值多少钱| 欧美一级毛片孕妇| 国产精品1区2区在线观看.| 宅男免费午夜| 国产美女午夜福利| 国产精品99久久久久久久久| 久久精品国产自在天天线| 色播亚洲综合网| 真实男女啪啪啪动态图| 国产爱豆传媒在线观看| 国内毛片毛片毛片毛片毛片| 色av中文字幕| 一区二区三区激情视频| 哪里可以看免费的av片| e午夜精品久久久久久久| 又黄又粗又硬又大视频| 精品国产三级普通话版| 国产精华一区二区三区| 午夜视频国产福利| 啦啦啦免费观看视频1| 男人的好看免费观看在线视频| 丁香六月欧美| 97人妻精品一区二区三区麻豆| 亚洲精华国产精华精| 人妻丰满熟妇av一区二区三区| 伊人久久大香线蕉亚洲五| а√天堂www在线а√下载| 国产又黄又爽又无遮挡在线| 精品电影一区二区在线| 51午夜福利影视在线观看| 亚洲第一电影网av| 色哟哟哟哟哟哟| 舔av片在线| 亚洲人成网站在线播| 日韩中文字幕欧美一区二区| 国产一区二区激情短视频| 国产v大片淫在线免费观看| 国产精品亚洲一级av第二区| 九色成人免费人妻av| 97超级碰碰碰精品色视频在线观看| 在线十欧美十亚洲十日本专区| 国产91精品成人一区二区三区| 国产美女午夜福利| 成人高潮视频无遮挡免费网站| 在线免费观看的www视频| 51国产日韩欧美| 欧美日韩瑟瑟在线播放| 亚洲精品色激情综合| 日本黄色片子视频| 在线观看免费视频日本深夜| 天堂√8在线中文| 性欧美人与动物交配| 国产老妇女一区| 女生性感内裤真人,穿戴方法视频| 国产色爽女视频免费观看| 国产一级毛片七仙女欲春2| 听说在线观看完整版免费高清| 国产精品久久电影中文字幕| 我的老师免费观看完整版| 制服丝袜大香蕉在线| 欧美大码av| 国产三级中文精品| 啪啪无遮挡十八禁网站| 国产精品野战在线观看| www日本黄色视频网| 午夜精品在线福利| 好看av亚洲va欧美ⅴa在| 国产精品久久久人人做人人爽| 一进一出抽搐动态| 色老头精品视频在线观看| 88av欧美| 成人鲁丝片一二三区免费| 日本成人三级电影网站| 中文字幕熟女人妻在线| 色尼玛亚洲综合影院| 欧美+亚洲+日韩+国产| 欧美成人性av电影在线观看| 亚洲人成网站高清观看| 国产真实伦视频高清在线观看 | 大型黄色视频在线免费观看| 国产精品久久久久久人妻精品电影| 中文字幕久久专区| 欧美av亚洲av综合av国产av| 禁无遮挡网站| 亚洲一区二区三区不卡视频| 男女床上黄色一级片免费看| 亚洲一区高清亚洲精品| 亚洲电影在线观看av| 91在线精品国自产拍蜜月 | 十八禁网站免费在线| 老司机福利观看| 日韩免费av在线播放| 啪啪无遮挡十八禁网站| 日日摸夜夜添夜夜添小说| av欧美777| 18禁黄网站禁片午夜丰满| h日本视频在线播放| 嫩草影院精品99| 99热这里只有是精品50| 美女免费视频网站| 成人av在线播放网站| 亚洲精品一卡2卡三卡4卡5卡| 亚洲无线在线观看| 国产色爽女视频免费观看| av国产免费在线观看| 久9热在线精品视频| 成年版毛片免费区| 最近视频中文字幕2019在线8| 欧美大码av| 在线视频色国产色| 国产淫片久久久久久久久 | xxx96com| 91av网一区二区| 在线观看日韩欧美| 免费av不卡在线播放| 欧美一区二区亚洲| 亚洲欧美日韩高清在线视频| 人人妻人人澡欧美一区二区| 国产精华一区二区三区| 国产精品乱码一区二三区的特点| 变态另类丝袜制服| 成人18禁在线播放| 日本与韩国留学比较| 男人和女人高潮做爰伦理| 日韩精品中文字幕看吧| 亚洲无线在线观看| 最新美女视频免费是黄的| 国产真实伦视频高清在线观看 | 91av网一区二区| bbb黄色大片| 亚洲成人免费电影在线观看| 精品免费久久久久久久清纯| 成人欧美大片| 91麻豆精品激情在线观看国产| xxxwww97欧美| 1000部很黄的大片| 老鸭窝网址在线观看| 在线播放无遮挡| 亚洲精品色激情综合| 欧美另类亚洲清纯唯美| 国产av在哪里看| 成人av在线播放网站| 又紧又爽又黄一区二区| 性色av乱码一区二区三区2| 看片在线看免费视频| 在线免费观看的www视频| 精品午夜福利视频在线观看一区| 一进一出抽搐gif免费好疼| a在线观看视频网站| 成人高潮视频无遮挡免费网站| 噜噜噜噜噜久久久久久91| 欧美一级a爱片免费观看看| 啦啦啦免费观看视频1| 秋霞在线观看毛片| 成年人午夜在线观看视频 | 国产视频首页在线观看| av在线播放精品| 亚洲婷婷狠狠爱综合网| 极品少妇高潮喷水抽搐| 欧美不卡视频在线免费观看| 十八禁网站网址无遮挡 | 色哟哟·www| 蜜桃亚洲精品一区二区三区| 精品一区二区三区人妻视频| 啦啦啦中文免费视频观看日本| 一级毛片aaaaaa免费看小| 婷婷色综合大香蕉| 精品一区二区三区视频在线| 内射极品少妇av片p| 日韩精品有码人妻一区| 人妻少妇偷人精品九色| 成年女人在线观看亚洲视频 | 国产成人精品福利久久| 欧美bdsm另类| 性插视频无遮挡在线免费观看| 中文字幕久久专区| 欧美精品一区二区大全| 国产国拍精品亚洲av在线观看| 成人鲁丝片一二三区免费| 肉色欧美久久久久久久蜜桃 | 国产亚洲最大av| videos熟女内射| 国产亚洲91精品色在线| 日韩在线高清观看一区二区三区| 久久久久久久午夜电影| 亚洲精品国产av蜜桃| 91午夜精品亚洲一区二区三区| 亚洲国产最新在线播放| 精品久久久噜噜| 国产一区二区在线观看日韩| 久久久久性生活片| 一级毛片黄色毛片免费观看视频| 一级毛片久久久久久久久女| 啦啦啦韩国在线观看视频| 亚洲欧洲国产日韩| 97精品久久久久久久久久精品| 久久精品国产亚洲av涩爱| 国产国拍精品亚洲av在线观看| 男人狂女人下面高潮的视频| 国内精品美女久久久久久| 国产av国产精品国产| 国产在线男女| 久久这里只有精品中国| 成人二区视频| 边亲边吃奶的免费视频| 69av精品久久久久久| 成人亚洲欧美一区二区av| 熟女人妻精品中文字幕| 97人妻精品一区二区三区麻豆| 精品亚洲乱码少妇综合久久| 麻豆国产97在线/欧美| 欧美人与善性xxx| 国产精品久久视频播放| 成人毛片a级毛片在线播放| 久久精品久久久久久久性| 精品国产三级普通话版| 中文字幕久久专区| 视频中文字幕在线观看| 九九在线视频观看精品| 免费黄色在线免费观看| 久久精品久久久久久噜噜老黄| 国产一区有黄有色的免费视频 | 国产成人a区在线观看| 国产视频首页在线观看| 三级男女做爰猛烈吃奶摸视频| 国产淫片久久久久久久久| 街头女战士在线观看网站| 亚洲国产精品国产精品| 国产一区二区亚洲精品在线观看| 天天一区二区日本电影三级| 美女cb高潮喷水在线观看| 神马国产精品三级电影在线观看| 熟妇人妻不卡中文字幕| 亚洲内射少妇av| 国产成人freesex在线| 国产在视频线精品| 日产精品乱码卡一卡2卡三| 欧美高清性xxxxhd video| 国产国拍精品亚洲av在线观看| 亚洲欧美日韩东京热| 国产亚洲av嫩草精品影院| 国产av在哪里看| 国产在线一区二区三区精| 国产精品日韩av在线免费观看| 久久鲁丝午夜福利片| 2021少妇久久久久久久久久久| 亚洲成人精品中文字幕电影| 国产亚洲av片在线观看秒播厂 | 午夜免费观看性视频| 一二三四中文在线观看免费高清| 国产精品一区二区性色av| 干丝袜人妻中文字幕| 我的女老师完整版在线观看| 亚洲精品一二三| videossex国产| 亚洲电影在线观看av| 色播亚洲综合网| 伦精品一区二区三区| 久久久精品94久久精品| 国产成人a区在线观看| 精品久久久久久电影网| 女人被狂操c到高潮| 蜜臀久久99精品久久宅男| 在现免费观看毛片| 99re6热这里在线精品视频| 久久精品综合一区二区三区| 亚洲在线观看片| 亚洲一级一片aⅴ在线观看| 日日啪夜夜爽| 婷婷六月久久综合丁香| 成人亚洲精品av一区二区| av播播在线观看一区| 少妇熟女欧美另类| 寂寞人妻少妇视频99o| videossex国产| 免费无遮挡裸体视频| 国产av国产精品国产| 国产伦理片在线播放av一区| 亚洲国产欧美人成| 乱码一卡2卡4卡精品| 菩萨蛮人人尽说江南好唐韦庄| 蜜桃久久精品国产亚洲av| 精品国产三级普通话版| 高清毛片免费看| 少妇被粗大猛烈的视频| 国产精品伦人一区二区| 国产乱人视频| 水蜜桃什么品种好| 午夜日本视频在线| 一本久久精品| 成年免费大片在线观看| 国产午夜精品久久久久久一区二区三区| 久久99精品国语久久久| 免费看av在线观看网站| 亚洲成人精品中文字幕电影| 女的被弄到高潮叫床怎么办| 日本wwww免费看| 日本三级黄在线观看| 国产精品一及| 久久久精品免费免费高清| 久久久久久久国产电影| 精品人妻偷拍中文字幕| av女优亚洲男人天堂| 亚洲精品亚洲一区二区| 婷婷色av中文字幕| 亚洲美女视频黄频| 在线观看人妻少妇| 日本猛色少妇xxxxx猛交久久| 18+在线观看网站| 搡老妇女老女人老熟妇| 亚洲成色77777| 国产亚洲av嫩草精品影院| 国产乱人视频| 久久草成人影院| 久久久久精品久久久久真实原创| 国产在视频线在精品| 我的女老师完整版在线观看| 午夜免费观看性视频| 日韩成人av中文字幕在线观看| 久久精品国产亚洲av涩爱| 免费观看性生交大片5|