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

    Effect of solute atom adsorption on heterogeneous nucleation by in-situ MgO particles: Experimental and theoretical studies

    2023-12-27 12:42:52LilingMoHengbinLioXiongZhouYuJunZhoJunDu
    Journal of Magnesium and Alloys 2023年10期

    Liling Mo ,Hengbin Lio ,Xiong Zhou ,Yu-Jun Zho ,Jun Du,?

    a Department of Metallic Materials, School of Materials Science and Engineering, South China University of Technology, 381 Wushan Rd, Guangzhou,Guangdong 510640, China

    b Department of Physics, South China University of Technology, 381 Wushan Rd, Guangzhou, Guangdong 510640, China

    Abstract The introduction of oxide inclusions during the smelting process has a clear promising heterogeneous nucleation potency on Mg-based alloys,but the mechanism has not been explored clearly yet.In the present work,the grain refinement mechanism of MgO in pure Mg,Mg-3Al and AZ31 (Mg-2.9Al-0.9Zn-0.3Mn) alloys is investigated by combing first-principles calculations and experiments.The theoretical results show that solute atoms adsorption will affect the nucleation and the subsequent growth process,which is an important factor affecting the refinement efficiency.A contradiction between the experimental results and the grain growth restriction factor (GRF) theory is observed,that is the refinement ratio of AZ31 is worse than Mg-3Al alloy.This is explained by an adsorption model which reveals that Al promotes the adsorption of Mg on MgO surface so as to stimulate more particles available as nucleating sites.Meanwhile,Fe and Mn also have favorable effects on the adsorption of Mg,Zn may play the opposite role compared to Al.The theoretical analyzes provide a strong support to the experiments that Al benefits the initial nucleation of α-Mg on MgO to promote the grain refining effect of Mg-3Al prior to AZ31 alloy.

    Keywords: Magnesium alloy;First-principles calculation;MgO;Grain refinement;Adsorption.

    1.Introduction

    Magnesium (Mg) and its alloys have attracted extensive attention for their excellent properties,and are considered to have broad application prospects in industry [1–3].However,some undesirable properties such as low mechanical properties and poor workability limit their wider applications and developments.Grain refinement can simultaneously improve the mechanical properties and processing abilities of Mg-based alloys and further broaden their applications effectively [4].

    Over the past decades,some grain refinement techniques have been developed [5–10].Nevertheless,many researchers and organizations are still striving for an efficient and stable refining method.Among them,most of the attempts are concentrated on adding nonmetallic particles or making themin-situform in the Mg melt [4,11,12].Noting that Mg has a high tendency to oxidize and will produce oxide inclusions during smelting process.Various studies have shown that oxides can act as inoculation refiners and reduce the grain size of Mg alloys effectively.Fu et al.[13] studied the refining effect of ZnO by experiments and the edge-to-edge matching(E2EM)model.They found the amounts of ZnO would be reduced by the reaction of ZnO+Mg→MgO+Zn and the grain refinement efficiency was limited.Men et al.[14] used the intensive melt shearing method to investigate the grain refinement mechanism of AZ91 alloys,where the dispersive MgO particles in the melt were thought to be potent nucleation sites.Wang et al.[15] carried out transmission electron microscope (TEM) studies,revealing that there are interfacial orientation relationships betweenα-Mg and MgO during solidification of [20]Mg(0001)Mg||[011]MgO(11)MgOand[23]Mg(101)Mg||[110]MgO(1)MgO.However,there still exist some ambiguous questions about the mechanism of oxide inoculation process,which is difficult to be observed by experiments.In recent years,theoretical calculation has attracted more and more attention and is used to assist experimental research.The first-principles calculation method has been applied to study surface and interface properties between heterogeneous nuclei and substrates [16–20].Several research groups have investigated the interface properties theoretically,such as ZnO/Mg [21],MgO/Mg [22],and so on.

    Most of experiments and theoretical calculations have been carried out from the perspective of interface properties so far.Besides,there are no stable refiners forα-Mg till now.This may be actually critical to the high activity of Mg melt and the effect of solute atoms adsorption on the heterogeneous nucleation particles.The influence of solute atoms on heterogeneous nucleation process was discussed by the interdependence theory in experiments frequently [23–26].Furthermore,the adsorption of solute atoms on nuclei surface during the initial nucleation stage is worthy to be thought profoundly.Cantor and Kim [27,28] developed a model for atomic adsorption on heterogeneous nucleation.They pointed out that the adsorption during the solidification process will affect the efficiency of nuclei as well as the subsequent grain growth process.Therefore,the common alloying and impurity elements may affect the structural stability and performance of nucleating catalyst,including Ca,Al,Zn,Mn,Fe,etc.Yao et al.[29] and Wei et al.[30] studied the influence of Mn and Fe adsorption on Al4C3(0001) surface through first-principles calculations.Yang et al.used density functional theory (DFT) method to analyze the influence of Ca on Al4C3[31],Al2MgC2[32] and MgAl2O4[33].According to these simulations,adsorbed atoms will change the surface configuration of nuclei,which in turn affect their heterogeneous nucleation performance.

    Following studies in the past years,it is clear that MgO can be the potent heterogeneous nucleation forα-Mg grain[14,15,22,34].Besides,it is easy to form MgO during the smelting process.Al/Zn/Mn are the main additive elements in commercial Mg alloys,such as AZ and AM series alloys.While Fe is the inevitable impurity element during the smelting process [35].It is believed that the existence of the solute elements has an indispensable influence on the grain refinement efficiency.Unfortunately,the adsorption details of the solute atoms on MgO surface are not well studied,leading to a blurry picture on the initial heterogeneous nucleation process.

    Here,we aim to investigate the influence of solute atoms(including Al,Zn,Mn and Fe)on inoculation potency of MgO in Mg-based alloys by combining experiments and theoretical calculations.The adsorption and co-adsorption of Al,Zn,Mn,Fe and Mg on MgO surface are simulated.We demonstrate that solute atoms adsorption will affect the nucleation and the subsequent growth process,consequently affecting the grain refinement efficiency.The current work provides a theoretical insight to better understand the grain refinement mechanism of Mg-based alloys inoculated by oxide inclusions.

    2.Experiments and theoretical calculations

    2.1.Experimental procedures

    Pure Mg,Mg-3Al,and commercial AZ31 (Mg-2.9Al-0.9Zn-0.3Mn) (wt.%,the same below) alloys were chosen as experimental materials in this work.Here,pure Mg (99.95%)and AZ31 were provided by Shanxi Regal Advanced Materials Co.,Ltd.,China.High purity aluminum (99.99%) was provided by ZhongNuo Advanced Material (Beijing) Technology Co.,Ltd.,China.Mg-3Al was synthesized by using pure Mg and high purity Al.Two series of experiments were designed,i.e.,pure Mg,Mg-3Al and AZ31 alloys without or with inoculation.For alloys with inoculation,the melt was treated byin-situoxidation process,in which the mixed gas of O2and Ar was introduced to make the oxide inclusionsin-situformed.According to reference[36],the gas flow rates of O2and Ar were set to be 10 and 100 ml/min,the rotor flowmeter was used to control the flow.

    Firstly,the raw materials with different compositions(250 g in total)were placed into the MgO crucible and melted in a resistance furnace under a mixed protective atmosphere of 99.5 vol.% N2and 0.5 vol.% SF6.Then the melt was untreated or inoculated for 10 min,and finally poured into a cylindrical low carbon steel mold(? 40 mm × 60 mm) preheated to 500 °C.Here,the casting temperature was maintained at 760 °C.Afterwards,all the ingots were sectioned,ground,polished,and etched by the picric acid solution (picric acid 2.1 g,acetic acid 5 ml,ethanol 40 ml and water 5 ml).

    Grain morphologies were observed by optical camera.Then the microstructure of the samples and the heterogeneous nuclei were characterized via scanning electron microscope(SEM,Quanta 200) equipped with energy-dispersive X-ray spectroscopy (EDS).The chemical concentrations of samples were determined via the inductively coupled plasma spectrometry (ICP,Agilent 5100).According to ASTM E112-88,average grain size was evaluated by the linear intercept method.For pure Mg,the grain size was evaluated by rings taken along different diameters of samples.The columnar crystal width was used to characterize the grain size.For Mg-3Al and AZ31,the area for measurement was the inner region of the samples with a diameter of ?30 mm.Grain refining ratio(δ) was obtained by Eq.(1):

    2.2.Computational details

    Fig.1.MgO (11) surfaces saturated by the pseudo-hydrogen: (a) Mg-terminated and (b) O-terminated MgO polar surfaces;(c) structure of the pseudomolecule (MgH4 and OH4).

    Based on density functional theory (DFT),all the firstprinciples calculations were conducted by the Viennaab initiosimulation package(VASP)[37].The exchange-correlation functional was described by the generalized gradient approximation and Perdew-Burke-Ernzerhof (GGA-PBE) [38].The interaction between the ions and the electrons was described by the standard frozen-core potentials,which were characterized by 3svalence state for Mg atom,2s,2pfor O atom,3s,3pfor Al atom and 3d,4sfor Zn/Mn/Fe atoms,respectively.Spin polarization was considered when the magnetic atoms of Fe or Mn were involved.The kinetic energy of the plane wave was set to 520 eV based on convergence tests.Brillouin zone was sampled with 9×9×6 and 9×9×1 Monkhorst-Pack k-point meshes for the bulk and slab calculations.All the atomic positions were fully optimized until the force was converged to 0.03 eV/?A for each atom,and the total energy change criterion was chosen to 1 × 10-5eV/atom.

    MgO exhibits the NaCl type simple cubic structure,with the space group of Fmm (No.225).The relaxed lattice parameters area=b=c=4.247,which is consistent with experimental data (fora=b=c=4.217 ?A).According to the previous studies,{100} and {111} MgO facets were observed in Mg alloys.The {111} facet is relatively stable in thermodynamics and kinetics,and its research is still inadequate.According to the TEM results,the(11) surface in {111} was selected for the current research[15,22,39].A MgO(11)-(2×2) multilayer slab model with periodic boundaries was adopted for surface simulation.MgO(11) surface consists of Mg and O layers alternatively aligned along the z axis.Mg-terminated and O-terminated MgO(11) polar surfaces were taken into consideration,as illustrated in Fig.1(a,b).The residual internal electric field in the MgO slabs caused by the excessive atomic characteristic was weaken through the adsorption of charged substances,like pseudo-hydrogens [17].1/2 and 3/2 artificial hydrogen atoms were added to the surface layer atoms for both O-termination and Mg-termination,respectively (as shown in Fig.1).The vacuum thickness was set to 15 ?A.It is known that the polarized ideal surfaces are likely unstable and often undergo a surface reconstruction [40].Here,the reconstruction is not considered as the polarized surfaces appear as part of the interface in the final models.

    3.Results and discussion

    3.1.Grain refinement of Mg alloys

    Fig.2 shows the macro grain morphologies of pure Mg,Mg-3Al and AZ31 alloys without or with inoculation.Fig.3 is the corresponding average grain size and grain refining ratio.Pure Mg possesses coarse columnar structure,as shown in Fig.2(a).Noting that only the average width of columnar structure was counted to analyze the effect of oxide inoculation,similar method has been applied in the literature[26].The Mg-3Al and AZ31 alloys exhibit an equiaxed crystal morphology,with the average grain sizes of 1140 μm and 1612 μm,respectively.All of the alloys show an effective refinement after inoculation,as shown in Fig.2 and Fig.3.

    Fig.2.Grain morphologies of as-cast pure Mg,Mg-3Al and AZ31 alloys without and with inoculation by in-situ oxidation process.

    Fig.3.Average grain size and corresponding grain refining ratio of pure Mg,Mg-3Al and AZ31 alloys without and with inoculation by in-situ oxidation process.

    The relative grain refining ratios are about 45.0%,83.2% and 44.1%,respectively.

    Fig.4.Typical SEM micrographs and EDS results of as-cast pure Mg,Mg-3Al and AZ31 alloys before and after inoculation.

    Fig.4 illustrates the microstructures of pure Mg,Mg-3Al and AZ31 alloys before and after inoculation,as well as the EDS results.Pure Mg exhibits a single-phase ofα-Mg matrix,as shown in Fig.4(a).For the Mg-3Al alloy,there are two phases that detected to beα-Mg matrix and Mg17Al12(Fig.4(d)).Three phases are found in the commercial AZ31 alloys,which areα-Mg,Mg17Al12,and Al8Mn5,as shown in Fig.4(g).Besides,there are some Zn-rich Mg-Al particles,which could be the Mg-Al-Zn ternary phase according to reference [41] or may be due to Zn atoms adsorption on Mg17Al12surface.It is clear that MgO particles arein-situformed after inoculation (see Fig.4).Besides,some Al-rich MgO particles were produced in Al-containing Mg alloys after inoculation (see Fig.4(e,f,h,i)).The existence of Al solute elements may lead to the formation of MgAl2O4according to reference[36],but is not considered here for its relatively low amount.In addition,Al2O3may be formed during the melting process.However,it is unstable in the Mg melt and finally converted into MgO or MgAl2O4,similarly to other oxidation products such as MnO2and Fe2O3.In addition,the competition relationship of MgO,Al2O3,MgAl2O4is analyzed from the perspective of chemical potential,as shown in Fig.5.The results showed that in the case of Mg-rich and Al-poor (the current experimental condition),MgO is the host phase,with the chemical potential of Al increases,MgAl2O4and Al2O3form gradually.Due to the strong oxygen affinity of Mg,in the case of introducing oxygen into the Mg melt,MgO is easier to form than other oxide inclusions from kinetic point.Moreover,most of the oxide products observed in this experiment are MgO particles.Therefore,MgO is the primary oxide in pure Mg and Mg-based alloys with inoculation,which is confirmed to be potent heterogeneous nucleation substrate here.

    According to the GRF theory [42,43],nucleation of pure Mg cannot occur by constitutional undercooling.The initial thermal undercooling generated by the mould provides the condition for the formation of grains,which further grow into columnar grains.Then the grain morphology is changed from columnar to equiaxed with alloying elements addition.In the case of being inoculated,MgO particles arein-situformedand become the heterogeneous nucleation sites forα-Mg.All of the pure Mg,Mg-3Al and AZ31 show an effective refinement.It is obvious that the grain refining ratios of these alloys are different,and the ratio of Mg-3Al is significantly greater than AZ31 alloy.Generally,the grain growth restriction factorQis greater with the solute elements increasing.Alloys with higherQvalue always exhibit better refinement.Qis calculated by the formula ofQ=mc0(k-1),herec0is the alloy’s composition,m(k-1) can be known from reference [44].The chemical compositions of pure Mg,Mg-3Al and AZ31 alloys before and after inoculation can be seen in Table 1.TheQvalue of inoculated AZ31 is 0.183,which is larger than inoculated Mg-3Al of 0.135.However,the refinement ratio of AZ31 alloy is lower than Mg-3Al alloy.This is contradictory to the GRF theory.

    Table 1 Chemical compositions (wt.%) of pure Mg,Mg-3Al and AZ31 alloys before and after inoculation.

    Fig.5.Chemical potential domain calculation of Mg-Al-O system.The allowed region for ΔμMg, ΔμAl,and ΔμO at (ΔμMg, ΔμAl) plane.All the heat of formations for the compounds are from DFT calculated results.

    The difference of solute elements in these samples is the most obvious change observed.Considering the solute elements are easily segregated in front of the S/L interface,thus it will adsorb on the nucleus surface and influence the adsorption of Mg during the initial stage of solidification.Besides,at the beginning of the solidification,the content of solidified Al8Mn5is far less thanα-Mg at the same temperature [45].Therefore,the effect of Al8Mn5on the adsorption performance is not consider here.Mg17Al12is also excluded for the low precipitating temperature.Limited by the experimental and analytical methods,it is difficult to analyze the influence mechanism of solute atoms on the adsorption and growth process on nuclei.Therefore,the adsorption model is constructed by first-principles calculation to explore whether the solute elements Al,Zn,Mn and Fe affect the surface structure and heterogeneous nucleation of MgO.

    3.2.Adsorption simulation

    3.2.1.Surface properties

    It is important to make sure that both Mg-terminated and O-terminated MgO(11) surface slabs are thick enough to stabilize the surface characteristics.Therefore,the change of interlayer spacings and corresponding surface energies of different layers were performed to determine the minimum thickness of the slab models.Surface energies of different models saturated by pseudo-H were calculated as follows [46–49]:

    This is derived from the following formulas:

    Fig.6.The changes of interlayer spacing of MgO polar surfaces saturated with pseudo-H: (a) O-terminated and (b) Mg-terminated MgO (11) surface.

    Table 2 The surface energies of Mg-terminated and O-terminated MgO (11) surface slabs.

    Table 2 The surface energies of Mg-terminated and O-terminated MgO (11) surface slabs.

    The symbol “?” indicates that layer meets the convergence criterion and was adopted for the latter calculations.

    The calculated changes of interlayer spacings are shown in Fig.6 and the surface energies of optimized MgO(11)surface slabs with different thickness are listed in Table 2.It shows that as the layer thickness increases,the change of the interlayer spacing for both Mg-terminated and Oterminated surfaces decreases gradually.Meanwhile,corresponding surface energies of these different slabs are well converged(as listed in Table 2).The differences of surface energies between the 10th-layer and 12th-layer are 0.01 eV/?A2(Δσ10-12=0.2 J/m2both for Mg-rich and O-rich condition)for O-termination and 0.004 eV/ ?A2(Δσ10-12=0.06 J/m2for Mg-rich andΔσ10-12=0.07 J/m2for O-rich situation) for Mgtermination,respectively.Besides,the relaxation of atomic positions is also converged according to the change of interlayer spacing in Fig.6.Therefore,12 atomic layers Mg-termination and O-termination were selected for the subsequent investigation.

    Furthermore,surface energy versus chemical potentialΔμMgwas calculated to discuss the stability of O-terminated and Mg-terminated MgO surfaces,as shown in Fig.7.It is clear that O-termination is energetically favored than Mgtermination under O-rich condition,while Mg-terminated surface is more stable under Mg-rich condition.Surface energies of these configurations have an intersection atΔμMg=-3.74 eV,indicating that both O-terminated and Mgterminated surfaces may exist and should be considered.Considering that the actual smelting and cooling process is in the Mg melt environment,the subsequent calculations and analyzes mainly focus on the Mg-termination.

    3.2.2.Stability for solute atoms adsorption on MgO(11)surface

    Fig.7.Surface energies versus ΔμMg for O-terminated and Mg-terminated MgO (11) surfaces.

    Adsorption of Al,Zn,Mn,Fe and Mg on MgO were discussed to describe the influence of these atoms on initial nucleation process ofα-Mg.The 12-layers O(Mg)-terminated MgO(11) surfaces were adopted according to previous calculation.Four high symmetry adsorption sites i.e.,H1 (Hollow site 1),H2 (Hollow site 2),T (Top site) and B (Bridge site),were taken into consideration.Adsorption sites of H1,H2 and T are shown in Fig.8,while configuration of B is not shown for its instability of adsorbing Mg or solute atoms.For Mg adsorption,coverages of 0.25–1 ML were considered,only 0.25 ML and 0.5 ML for other solute atoms because the concentration of these atoms is generally low in the Mg melt.Adsorption energies of Mg or solute atoms adsorption on MgO(11) surfaces can be given by:

    hereEis the energy of MgO(11) clean surface,niis the number of Mg or solute atoms,respectively.Eiis the energy of the isolated adatom,andEis the total energy of MgO(11) surface with adatoms.The calculated individual adsorption energies of Al,Zn,Mn,Fe are listed in Table 3,while the adsorption energies of Mg on different conditions are shown in Fig.9.

    As listed in Table 3,the adsorption energies at site H2 are higher than other sites under 0.25 ML and 0.5 ML coverages for Al,Zn,Mn and Fe on O-termination.The initial adsorption of these atoms on T site will be relaxed to H2 spontaneously,and their x and y coordinates were fixed here to calculate the energies for comparison.With the coverage increasing from 0.25 ML to 0.5 ML,the difference in adsorption energy among H1,H2 and T sites decreases.For Mg-terminated surface,the adsorption energies at site H1 are mostly higher than other sites under 0.25 ML,while sites H1 and H2 are almost the same as the coverage increases to 0.5 ML.Similar phenomena can also be found for Mg adsorption on O-terminated and Mg-terminated MgO surfaces.With the coverage increasing,the adsorption energy for Mg decreases on O-termination while it increases on Mg-termination(as shown in Fig.8Fig.9).

    Comparing the energies of Al,Zn,Mn,Fe adsorption on MgO surface,the adsorption energies follow the sequence ofEad(Al)>Ead(Fe)>Ead(Mn)>Ead(Mg)>Ead(Zn)both for O-terminated and Mg-terminated surfaces.It is clear that the adsorption energy of Al is greater than Mg cases of 0.25–1 ML and other elements on both O-termination and Mg-termination at the same adsorption sites and coverages.This indicates that Al is more easily to adsorb on MgO(11)surfaces,the same as Mn and Fe.As for Zn adsorption,the adsorption energy is less than Mg and other solute atoms.According to Eq.(8),the larger positive value of adsorption energy means the more stable structure.The adsorption energies of Zn are positive under different conditions,indicating that Zn is likely to have a tendency to adsorb on the MgO.

    To study the effect of solute atoms on the initial nucleation process ofα-Mg,co-adsorption of Mg and Al/Mn/Fe on H1 and H2 sites of Mg-terminated MgO surfaces were considered.The adsorption energy at the condition of Mg-X(X=Al,Mn,Fe) co-adsorption on the MgO(11) surfaces can be calculated by:

    Eis the total energy of the slab that single solute atom adsorption,andE(X+Mg)/MgO(1?1)is the total energy of Mg and solute atom co-adsorption.The actual coverage of Al or Mg depends on their chemical potentials in the experiment.Here the chemical potential of Mg is about zero with respect to bulk Mg as Mg melt exists in the sample,while the chemical potential of Al will be much smaller as the content of Al is relatively small.This is corresponding to a Mg-rich and Al-poor condition in the experiment.Therefore,only the low coverage of Al is considered in this work following the practical condition,similarly to other solute atoms.Based on above discussion and the computational cost of DFT calculations,the coverage ratio of Mg and solute atoms were selected to be 3:1 in the model construction to simulate the impact of solute atoms on the adsorption of Mg.

    Fig.8.The diagrammatic drawing of adsorption sites and co-adsorption on the 12-layers of (a) O-terminated and (b) Mg-terminated MgO (11) surfaces.

    Table 3 The adsorption energies of Al,Zn,Mn and Fe atoms on O(Mg)-terminated MgO (11) surfaces.(F) in the Table stands for fixing the x and y coordinate of adatoms for their instability.

    Table 3 The adsorption energies of Al,Zn,Mn and Fe atoms on O(Mg)-terminated MgO (11) surfaces.(F) in the Table stands for fixing the x and y coordinate of adatoms for their instability.

    The calculated co-adsorption energies of Mg/Al,Mg/Mn and Mg/Fe on Mg-termination are shown in Fig.10.There are some discrepancies on energies between different solute atoms and Mg co-adsorption.The co-adsorption energies follow the order ofEco-ad(MgAl)>Eco-ad(MgFe)>Eco-ad(MgMn)>Ead(1 ML Mg) for Mg-termination.It is obvious that all the energies of Mg/Al,Mg/Fe,Mg/Mn co-adsorption are greater than 1 ML Mg case on Mg-terminated MgO surfaces.This means that Al,Fe,Mn atoms will adsorb on MgO surface prior to Mg,then the existence of them promotes the adsorption of Mg atoms on MgO.Therefore,moreα-Mg grains are nucleating on oxide particles acting as nucleation catalyst.This will lead to the fact that the grain refining ratio of Mg-3Al is greater than pure Mg,so as to AZ31 alloy.Unfortunately,the effect of Zn on Mg adsorption is not clear from the theoretical point of view,which requires further investigation.

    Fig.9.The adsorption energies of Mg on (a) O-terminated and (b) Mg-terminated MgO (11) surfaces.

    Fig.10.The effects of Al,Fe and Mn on the adsorption energy of Mg on the Mg terminated MgO (11) surface.

    3.2.3.Electronic structure and bonding characteristics

    The structural evolution and performance of the nuclei after adsorption are closely associated with the bonding between the adatoms and surface atoms.Therefore,to have some insight into the surface properties more intuitively from a microscopic perspective,the charge density differences (CDD)for Mg and Al/Zn/Mn/Fe adsorption and co-adsorption at H1 for Mg-termination MgO surface were calculated here.The CDDΔρcan be calculated as:

    whereρis the electron density of the total i/MgO(11) adsorption system,ρiandρa(bǔ)re the unperturbed electron densities of isolated adatoms and MgO(11) surfaces,respectively.

    The calculated CDDs are shown in Fig.11.Here,the blue areas indicate charge accumulation while yellow and red areas represent charge depletion.It is worth noting that the rearrangement of the charge is mainly distributed to the adatoms and the top layer atoms on the surface.Taking Al adsorption and Mg-Al co-adsorption as an example,the CDDs of Mg and Al adsorption/co-adsorption in Fig.11 show that the charge mainly transfers from the adatoms to the middle regions,and there are charge accumulation between Mg-Mg and Mg-Al atoms,implying that polar covalent bonds are formed between metal-metal atoms.Similar phenomena can be found for other solute atoms (Zn,Mn,Fe) adsorption on Mg-terminated surface.

    Bader charge analysis was carried out to verify the charge transfer between adatoms and the top layer atoms,as shown in Fig.12.The charges of Mg atoms on the top layers of Al/Mn/Fe adsorption models are more than those in clean surface and Mg adsorption case,and the charge of the adsorbed metallic atoms are changed for single and co-adsorption case.For the case of Mg and Al co-adsorption,the charge of adsorbed Al gains 2.47 e,which is more than single Al adsorption of 1.21 e.Besides,adatoms seem to have stronger interaction with the surface atoms compared with single adsorption case.It indicates that solute atoms adsorption preferentially would influence the bonding on the surface,then strengthen the stability of surface structure and affect the refinement efficiency of heterogeneous nuclei.The calculated results of Bader charge are consistent with the CDDs.

    The density of states(DOS)and corresponding partial density of states (PDOS) were carried out to further explore the surface bonding after adsorption.Fig.13(a)indicates that bulk MgO is an insulator with a band gap of 3.32 eV (much smaller than the experimental gap of 6.36 eV due to the well-known underestimation of GGA calculations).The surface become metallic as the band gap disappears,mainly due to the charge residual from the Mg atoms on the surface for the unstable polarized surface.After single Al adsorption,the characteristic of metallic is degraded.The resonance of PDOS with different adsorbates is observed at around -20~-15 eV and -8~8 eV Besides,it can be seen that theselectron of Al atoms are slightly hybridized with s/pelectron of Mg around -5 eV for single alloying elements adsorption and co-adsorption,which shows the formation of covalent bonds between the adatoms and Mg atoms.Besides,the value of the PDOS at Fermi level for Mn,Fe and Mg adsorption and coadsorption in the Mg-terminated system are larger than clean surface and single element adsorption,reflecting the stronger metallic feature after Mg and Mn/Fe co-adsorption on surface.In consequence,the bonding characteristic between the adsorbed metallic atoms and the surface Mg-layer atoms on the Mg-termination is mainly covalent bonds.

    There are significant differences in PDOS results of different adsorption slabs,demonstrating that there are some discrepancies about the interaction between various adatoms and Mg atoms,and bond strength will be different for these atoms.

    Fig.12.Bader charge analysis of selected atoms for Al/Mn/Fe and Mg adsorption and co-adsorption on Mg-terminated MgO (11) surfaces (positive means get charge while negative is lost charge).Here,Mg 1–4 represent the adsorbed Mg atoms in the surface model.Mg 11–14 are Mg atoms on the top layer of MgO surface.

    This indicates that adsorption of solute atoms will significantly influence the surface properties of the grain refiners,consequently influencing their inoculation performance in the initial nucleation process,in line with the experiment results.

    3.3.Analysis of grain refinement mechanism

    Fig.13.(a) Density of states (DOS) of MgO bulk and corresponding projected density of states (PDOS) of Mg,O atoms;(b)-(f) Total density of states and PDOS for single Mg,Al,Zn,Mn,Fe adsorption and (g)-(h) 1 ML Mg,Mg/Al,Mg/Mn,Mg/Fe co-adsorption on Mg-terminated MgO surface,respectively.Here,dashed lines denote the Fermi level.

    Fig.14.(a–c) The schematic diagram of initial nucleation progress of pure Mg,Mg-3Al and AZ31 alloys inoculated by in-situ MgO;(d) adsorption process of MgO surface in Mg-3Al alloy.The orange dashed lines in (a)-(c) denote the molten Mg.

    MgO was observed forin-situoxidation inoculated pure Mg,Mg-3Al and AZ31 alloys,and is considered to be the heterogeneous nucleation sites in the present work.According to the interdependence theory,constitutional supercooling(CS) induced by solute elements offers an extra undercooling for heterogeneous substrates and plays an important role in connecting the nucleation and grain growth [23].Actually,growth restriction factorQis a calculation for development rate of the constitutional zone [43].The addition of solute elements with a largerQvalue not only increase the restricting effect for growth,but also create a stronger CS for stimulating heterogeneous nucleation,thereby increasing the grain refinement efficiency.However,we confirmed that AZ31 has a lower refinement ratio,though theQvalue is larger than Mg-3Al alloy.From the perspective of the nucleation and growth process,the initial nucleation is affected by the adsorption of solute elements.

    The schematic diagram of nucleation process for inoculated Mg-based alloys is shown in Fig.14.After being inoculated,MgO particles are formed in the melt.For the inoculated pure Mg,part of the particles can be the heterogeneous sites at the beginning of solidification.As for the Mg-3Al alloy,the dissociative Al atom in the melt has a strong adsorption capacity on MgO surface in the initial nucleation process.Besides,it can be observed from Table 3 and Fig.9 that the adsorption energy of Al is higher than Mg.For the commercial AZ31 alloy,the adsorption of Al and Mn will also increase the adsorption energy of Mg on MgO.These results suggest that Al is easier to adsorb on MgO surface than Mg.In the case of Al adsorption on Mg-terminated MgO(11)surface,the energies of Mg adsorption increase from 1.65 eV to 2.10 eV and 1.65 eV to 2.14 eV for H1 and H2 sites,respectively.The influence of increasing Al concentration on the calculation results of adsorption energy was considered to verified the accuracy.The co-adsorption energy was calculated when the coverage of Al is 0.5 ML.The results showed that with the coverage increases from 0.25 ML to 0.5 ML,the co-adsorption energies on Mg-termination increase from 2.10 eV to 2.25 eV and 2.14 eV to 2.33 eV for H1 and H2 sites,respectively.It is still higher than 1 ML Mg adsorption case.Since the experiment is in the Mg melt environment,the coverage of Al is not considered greater than 0.5 ML.These results provide important evidence that Al promotes the adsorption of Mg on MgO surface,and Fe and Mn seem to have weakly favorable role on the adsorption of Mg.For Zn adsorption case,because the single adsorption energy is lower than Mg,it is difficult to determine the influence of Zn on the adsorption of Mg from the view of theoretical point.According to the experimental results,it is deduced that Zn may be unfavorable to the adsorption of Mg on MgO surface,which deserved further exploration.

    The above analysis implies that the grain refinement efficiency of nucleating particles is related to the effect of solute atoms on the adsorption of Mg to MgO particles surface,for Al/Mn/Fe are positive.This means that Al adsorption on MgO would significantly promote the nucleating ofα-Mg and more particles available as nucleating sites,thus strengthen the ability of MgO for acting as heterogeneous nucleation substrate to a great extent,as shown in Fig.14(d).The theoretical results provide a powerful supporting for the experimental phenomena that Al has a positive effect on the adsorption and nucleation ofα-Mg grains on MgO particles,it is one of the reasons that the refinement efficiency of Mg-3Al is higher than AZ31 alloy.

    4.Conclusion

    Through a combination of experiments and first-principles calculations,the potent nucleant of pure Mg,Mg-3Al and AZ31 alloys inoculated byin-situoxidation process has been determined.Adsorption energies,electronic structures and bonding of Mg and Al/Zn/Mn/Fe adsorption on MgO(11)surface were studied.The conclusions are summarized as following:

    (1) MgO is proved to have a great grain refinement effect on pure Mg,Mg-3Al,and AZ31 alloys.The refining ratio of Mg-3Al is greater than AZ31 alloy,which is contradictory to the GRF theory.

    (2) Al is preferentially adsorbed on MgO surface than Mg and other solute atoms at Mg-termination and O-termination.Besides,the energies of Mg/X(X=Al/Mn/Fe) co-adsorption demonstrate that Al promotes the adsorption of Mg on MgO surface,Fe and Mn also have weakly favorable role on the adsorption of Mg.

    (3) The electronic structure and bonding characteristics of the adsorption structure show that polar covalent bonds are formed between adatoms and top layers atoms,which will influence the structure and performance of the nuclei.

    (4) The adsorption energy of Mg on MgO surface will change significantly after solute atoms adsorption.It concludes that solute atoms adsorb on MgO surface will affect the subsequent stacking and growth ofα-Mg on MgO.Therefore,the addition of solute atoms will influence the grain refining efficiency of primary Mg grains by oxide inoculation.

    Declaration of competing interest

    The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper.

    Acknowledgment

    This work was supported by the National Natural Science Foundation of China (No.51871100).Computing resources from national supercomputer center in Guangzhou are gratefully acknowledged.

    精品人妻熟女av久视频| 婷婷亚洲欧美| 91精品一卡2卡3卡4卡| 蜜桃久久精品国产亚洲av| 亚洲成人久久爱视频| 国产黄色小视频在线观看| 一级毛片电影观看 | 国产在线精品亚洲第一网站| 色尼玛亚洲综合影院| 免费看日本二区| 亚洲一区高清亚洲精品| 色视频www国产| 日韩视频在线欧美| 岛国毛片在线播放| 黄片无遮挡物在线观看| 国国产精品蜜臀av免费| 国产精品久久久久久精品电影小说 | 蜜桃亚洲精品一区二区三区| 99精品在免费线老司机午夜| 身体一侧抽搐| 亚洲av免费在线观看| 久久6这里有精品| 中文字幕精品亚洲无线码一区| 乱码一卡2卡4卡精品| 人人妻人人看人人澡| av天堂中文字幕网| 一区福利在线观看| 岛国在线免费视频观看| 一级av片app| 国产探花在线观看一区二区| 身体一侧抽搐| 我的女老师完整版在线观看| 18禁裸乳无遮挡免费网站照片| 极品教师在线视频| 久久久久久久久久成人| 亚洲欧美日韩无卡精品| 一个人免费在线观看电影| 日韩 亚洲 欧美在线| 国产成人精品一,二区 | 秋霞在线观看毛片| 欧美潮喷喷水| ponron亚洲| 高清日韩中文字幕在线| 一本久久中文字幕| 国产高清三级在线| 成人特级黄色片久久久久久久| 一区福利在线观看| 18禁在线无遮挡免费观看视频| 成人鲁丝片一二三区免费| 成人一区二区视频在线观看| 日本免费一区二区三区高清不卡| 麻豆国产97在线/欧美| 成人午夜高清在线视频| 成人二区视频| 国产午夜精品论理片| 亚洲激情五月婷婷啪啪| 九九热线精品视视频播放| 一本一本综合久久| 一级毛片aaaaaa免费看小| 亚洲中文字幕日韩| 亚洲自偷自拍三级| 亚洲欧美中文字幕日韩二区| 免费人成在线观看视频色| 高清午夜精品一区二区三区 | 丝袜喷水一区| 亚洲欧美精品自产自拍| 国产一区亚洲一区在线观看| 精品一区二区三区视频在线| 亚洲丝袜综合中文字幕| 欧美+亚洲+日韩+国产| 国产精品一区二区三区四区久久| 国产淫片久久久久久久久| 精品人妻一区二区三区麻豆| 国产精品久久久久久久电影| 日韩欧美三级三区| 久久久色成人| 99久久精品一区二区三区| 国产成人精品久久久久久| 国产精品久久久久久精品电影| 嫩草影院新地址| 少妇熟女欧美另类| 成人永久免费在线观看视频| 91aial.com中文字幕在线观看| 亚洲最大成人av| 可以在线观看的亚洲视频| 日韩大尺度精品在线看网址| 91久久精品电影网| 麻豆国产97在线/欧美| 好男人在线观看高清免费视频| 国产 一区 欧美 日韩| 黄色日韩在线| 亚州av有码| 一本久久精品| 国产精品永久免费网站| 熟女电影av网| av又黄又爽大尺度在线免费看 | 一级黄片播放器| av专区在线播放| 日本黄大片高清| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 国产成人一区二区在线| 国产一区二区三区av在线 | 在线a可以看的网站| 亚洲欧洲国产日韩| 亚洲婷婷狠狠爱综合网| 男女那种视频在线观看| 三级男女做爰猛烈吃奶摸视频| 国产高清激情床上av| 少妇猛男粗大的猛烈进出视频 | 亚洲国产欧美在线一区| 亚洲国产色片| 91精品一卡2卡3卡4卡| 日本五十路高清| 激情 狠狠 欧美| 人妻系列 视频| 看黄色毛片网站| 国产老妇女一区| 男女啪啪激烈高潮av片| 亚洲精品乱码久久久久久按摩| 久久久午夜欧美精品| 欧美bdsm另类| 日日撸夜夜添| 国产男人的电影天堂91| 最近视频中文字幕2019在线8| 中国美女看黄片| 国产亚洲欧美98| 日韩强制内射视频| 激情 狠狠 欧美| 成人性生交大片免费视频hd| 中文字幕免费在线视频6| 天堂av国产一区二区熟女人妻| 亚洲成人av在线免费| 3wmmmm亚洲av在线观看| 激情 狠狠 欧美| 禁无遮挡网站| 高清毛片免费看| kizo精华| 身体一侧抽搐| 男女下面进入的视频免费午夜| 黄片无遮挡物在线观看| 婷婷六月久久综合丁香| 久久久久九九精品影院| 天堂av国产一区二区熟女人妻| 波多野结衣高清无吗| 淫秽高清视频在线观看| 亚洲av不卡在线观看| 婷婷亚洲欧美| 波多野结衣高清作品| 国产亚洲av嫩草精品影院| 蜜桃亚洲精品一区二区三区| 秋霞在线观看毛片| 最后的刺客免费高清国语| 一进一出抽搐动态| 亚洲美女搞黄在线观看| 天天一区二区日本电影三级| 自拍偷自拍亚洲精品老妇| 久久婷婷人人爽人人干人人爱| 天美传媒精品一区二区| av在线播放精品| 国产单亲对白刺激| 欧美日韩国产亚洲二区| 午夜福利高清视频| 色吧在线观看| 插阴视频在线观看视频| 欧洲精品卡2卡3卡4卡5卡区| 国产精品野战在线观看| 亚洲电影在线观看av| 中文欧美无线码| 爱豆传媒免费全集在线观看| 中文字幕熟女人妻在线| 欧美3d第一页| 变态另类成人亚洲欧美熟女| 国产成人福利小说| 久久久久性生活片| 色综合色国产| 国产av不卡久久| 国产麻豆成人av免费视频| 18+在线观看网站| 国产真实伦视频高清在线观看| 床上黄色一级片| 午夜激情欧美在线| 免费av观看视频| 国产精品无大码| 搡老妇女老女人老熟妇| 精品久久久久久久久久久久久| 色综合色国产| 欧美bdsm另类| 人妻少妇偷人精品九色| 亚洲丝袜综合中文字幕| 午夜a级毛片| 午夜福利视频1000在线观看| 国产av一区在线观看免费| 日韩大尺度精品在线看网址| 少妇人妻精品综合一区二区 | 亚洲第一电影网av| 中文字幕人妻熟人妻熟丝袜美| 久久亚洲国产成人精品v| 99视频精品全部免费 在线| 国产av在哪里看| 日韩三级伦理在线观看| 少妇裸体淫交视频免费看高清| 禁无遮挡网站| 亚洲第一电影网av| 国产伦一二天堂av在线观看| 两个人的视频大全免费| 3wmmmm亚洲av在线观看| 97在线视频观看| 男人的好看免费观看在线视频| 久久久久久久久久久丰满| 一个人看的www免费观看视频| 亚洲精品色激情综合| 国产老妇伦熟女老妇高清| 午夜视频国产福利| 黄片无遮挡物在线观看| 亚洲av第一区精品v没综合| 99久国产av精品| 国产亚洲精品久久久久久毛片| 九色成人免费人妻av| 国产69精品久久久久777片| 国语自产精品视频在线第100页| 亚洲国产欧美在线一区| 国产精品美女特级片免费视频播放器| 中国美白少妇内射xxxbb| 哪个播放器可以免费观看大片| 欧美成人a在线观看| 青春草视频在线免费观看| 日本与韩国留学比较| 18禁在线无遮挡免费观看视频| 欧美潮喷喷水| 亚洲欧美成人综合另类久久久 | 大又大粗又爽又黄少妇毛片口| 日产精品乱码卡一卡2卡三| 偷拍熟女少妇极品色| 精品人妻熟女av久视频| 久久久久国产网址| 国产一区二区激情短视频| 亚洲av不卡在线观看| h日本视频在线播放| 久久精品久久久久久噜噜老黄 | 久久精品久久久久久久性| 国产老妇伦熟女老妇高清| 国产高清激情床上av| 丰满乱子伦码专区| 99久久精品国产国产毛片| 亚洲精品456在线播放app| 国产午夜精品一二区理论片| 全区人妻精品视频| 欧美bdsm另类| 国产高清有码在线观看视频| 人妻系列 视频| 亚洲丝袜综合中文字幕| 女同久久另类99精品国产91| 久久亚洲精品不卡| 91在线精品国自产拍蜜月| 久久人人爽人人片av| 亚洲最大成人中文| 亚洲国产欧美人成| 精品人妻一区二区三区麻豆| 国产 一区精品| 一本久久中文字幕| 黄色一级大片看看| 成人特级av手机在线观看| 亚洲美女搞黄在线观看| 日本与韩国留学比较| 欧美性猛交黑人性爽| 悠悠久久av| 欧美日韩在线观看h| 一级毛片电影观看 | 亚洲自偷自拍三级| 中国美女看黄片| 中文字幕制服av| 91狼人影院| 亚洲欧美精品专区久久| 欧美一区二区国产精品久久精品| 国产老妇伦熟女老妇高清| 久久久久久久亚洲中文字幕| 九九在线视频观看精品| 草草在线视频免费看| 亚洲精品影视一区二区三区av| 别揉我奶头 嗯啊视频| 九九在线视频观看精品| 爱豆传媒免费全集在线观看| 亚洲精品国产av成人精品| 国产精品伦人一区二区| 免费观看精品视频网站| 国产熟女欧美一区二区| 久久久成人免费电影| 一个人观看的视频www高清免费观看| 国产乱人偷精品视频| av天堂中文字幕网| 丰满的人妻完整版| 亚洲七黄色美女视频| 美女高潮的动态| 麻豆一二三区av精品| 午夜福利高清视频| 长腿黑丝高跟| 狂野欧美白嫩少妇大欣赏| 国产探花在线观看一区二区| 天堂中文最新版在线下载 | 狂野欧美白嫩少妇大欣赏| 亚洲久久久久久中文字幕| 别揉我奶头 嗯啊视频| 老熟妇乱子伦视频在线观看| 国产熟女欧美一区二区| 秋霞在线观看毛片| 成人二区视频| 老女人水多毛片| www.av在线官网国产| 日本五十路高清| 色综合站精品国产| 久久国内精品自在自线图片| 男女那种视频在线观看| 永久网站在线| 我的女老师完整版在线观看| 狂野欧美白嫩少妇大欣赏| 午夜福利成人在线免费观看| 久久亚洲国产成人精品v| 国产女主播在线喷水免费视频网站 | 菩萨蛮人人尽说江南好唐韦庄 | 国产精品国产高清国产av| 欧美日韩精品成人综合77777| 一个人观看的视频www高清免费观看| 中文字幕久久专区| 国产乱人视频| 18禁在线播放成人免费| 老司机福利观看| 国产精品.久久久| 精品一区二区免费观看| 亚洲欧美日韩高清专用| 国产在线精品亚洲第一网站| 成年女人看的毛片在线观看| 黑人高潮一二区| 精品一区二区三区视频在线| 国产伦精品一区二区三区四那| 波多野结衣高清无吗| 午夜亚洲福利在线播放| 老司机福利观看| 亚洲av二区三区四区| 日韩欧美三级三区| 人妻少妇偷人精品九色| 精品少妇黑人巨大在线播放 | 黑人高潮一二区| 99久国产av精品国产电影| 美女被艹到高潮喷水动态| 日本五十路高清| 欧美在线一区亚洲| 国产综合懂色| 亚洲国产精品sss在线观看| 久久国内精品自在自线图片| 人人妻人人看人人澡| www日本黄色视频网| 男人的好看免费观看在线视频| av在线亚洲专区| 久久久久久久久久成人| 日韩制服骚丝袜av| 久久午夜亚洲精品久久| 免费av观看视频| 深夜a级毛片| 亚洲精品国产av成人精品| 丝袜喷水一区| 久久精品人妻少妇| 亚洲七黄色美女视频| 日本一本二区三区精品| 大型黄色视频在线免费观看| 日日摸夜夜添夜夜爱| 日韩一区二区视频免费看| 日韩中字成人| 三级毛片av免费| 国产黄片美女视频| 亚洲国产欧洲综合997久久,| 99热这里只有精品一区| 国产高清不卡午夜福利| 国产一区二区激情短视频| 国产单亲对白刺激| 国产午夜精品论理片| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 亚洲国产欧美在线一区| 久久人人爽人人片av| 三级经典国产精品| 亚洲av.av天堂| 久久精品国产亚洲网站| 婷婷亚洲欧美| 国产私拍福利视频在线观看| 成年免费大片在线观看| av.在线天堂| 18禁在线播放成人免费| 免费看日本二区| 成人毛片a级毛片在线播放| 久久久久久大精品| 乱码一卡2卡4卡精品| 神马国产精品三级电影在线观看| 欧美成人a在线观看| 国产乱人偷精品视频| 国产免费男女视频| 久久婷婷人人爽人人干人人爱| 亚洲人与动物交配视频| 黄色欧美视频在线观看| 日韩av不卡免费在线播放| 欧美zozozo另类| 99精品在免费线老司机午夜| 在线观看一区二区三区| 久久6这里有精品| 一本精品99久久精品77| 99久久久亚洲精品蜜臀av| a级毛片免费高清观看在线播放| 久久久久性生活片| 亚洲成人av在线免费| 综合色av麻豆| 欧洲精品卡2卡3卡4卡5卡区| 欧美精品一区二区大全| 亚洲av第一区精品v没综合| 亚洲在线观看片| 黄片无遮挡物在线观看| 最近手机中文字幕大全| 尾随美女入室| 黑人高潮一二区| 婷婷色av中文字幕| 全区人妻精品视频| 欧美一区二区亚洲| 最近最新中文字幕大全电影3| 午夜福利在线在线| 国产成人影院久久av| 欧美+日韩+精品| 深夜精品福利| 国产日本99.免费观看| 午夜福利高清视频| 精品无人区乱码1区二区| 91午夜精品亚洲一区二区三区| 亚洲国产精品成人久久小说 | 久99久视频精品免费| 国产精品综合久久久久久久免费| 日本一二三区视频观看| 国产免费一级a男人的天堂| 久久久色成人| 国产乱人视频| 麻豆成人午夜福利视频| 日韩av不卡免费在线播放| 欧美xxxx性猛交bbbb| 日韩一本色道免费dvd| 国产黄色视频一区二区在线观看 | 丰满人妻一区二区三区视频av| 亚洲成a人片在线一区二区| 午夜福利高清视频| 午夜爱爱视频在线播放| 国产精品一区二区在线观看99 | 国产真实伦视频高清在线观看| 国产单亲对白刺激| 在线观看一区二区三区| av黄色大香蕉| 国产男人的电影天堂91| 看片在线看免费视频| 97在线视频观看| 欧美高清成人免费视频www| 欧美性猛交╳xxx乱大交人| 91狼人影院| 最近中文字幕高清免费大全6| 日韩国内少妇激情av| 成人午夜高清在线视频| 欧美最新免费一区二区三区| 性色avwww在线观看| 亚洲人成网站在线播放欧美日韩| 毛片女人毛片| 国产高潮美女av| 最近2019中文字幕mv第一页| 欧美日本亚洲视频在线播放| 老司机福利观看| 九九久久精品国产亚洲av麻豆| 久久久久久久久久久丰满| 高清日韩中文字幕在线| 亚洲三级黄色毛片| 色综合站精品国产| 成人永久免费在线观看视频| 卡戴珊不雅视频在线播放| 91av网一区二区| 成人性生交大片免费视频hd| 国产三级中文精品| 男人狂女人下面高潮的视频| 搡女人真爽免费视频火全软件| 日本三级黄在线观看| 美女国产视频在线观看| 国产探花在线观看一区二区| 久久久精品欧美日韩精品| 一个人看视频在线观看www免费| 欧美zozozo另类| 一级黄色大片毛片| 国产精华一区二区三区| 男女啪啪激烈高潮av片| 一级av片app| 热99在线观看视频| 国产成人精品一,二区 | 成年版毛片免费区| 午夜福利视频1000在线观看| 男女啪啪激烈高潮av片| 嘟嘟电影网在线观看| 热99在线观看视频| 天堂网av新在线| 九九热线精品视视频播放| 欧美+日韩+精品| 99久久无色码亚洲精品果冻| 欧美三级亚洲精品| 51国产日韩欧美| 自拍偷自拍亚洲精品老妇| 欧美成人免费av一区二区三区| 亚洲欧洲国产日韩| 国内精品久久久久精免费| 综合色av麻豆| 黄色一级大片看看| 禁无遮挡网站| 神马国产精品三级电影在线观看| 国产精品蜜桃在线观看 | 晚上一个人看的免费电影| 精品一区二区免费观看| 免费人成在线观看视频色| 中国国产av一级| 国产精品一区二区三区四区免费观看| 一本久久中文字幕| 91午夜精品亚洲一区二区三区| 欧美3d第一页| 免费av毛片视频| 青春草视频在线免费观看| 最近视频中文字幕2019在线8| 亚洲国产精品合色在线| 亚洲精品国产av成人精品| 日本免费一区二区三区高清不卡| 免费看日本二区| 亚洲av男天堂| 欧美三级亚洲精品| 亚洲人与动物交配视频| www日本黄色视频网| 乱码一卡2卡4卡精品| 亚洲欧洲日产国产| 伊人久久精品亚洲午夜| 亚洲精品国产av成人精品| 中文字幕av成人在线电影| 蜜桃亚洲精品一区二区三区| 免费电影在线观看免费观看| 波多野结衣巨乳人妻| 狠狠狠狠99中文字幕| 亚洲国产高清在线一区二区三| 国产精品一区www在线观看| 欧美日本亚洲视频在线播放| 色尼玛亚洲综合影院| 亚洲内射少妇av| 伦精品一区二区三区| 国模一区二区三区四区视频| 国产精品综合久久久久久久免费| 国产一级毛片在线| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 91aial.com中文字幕在线观看| 日本在线视频免费播放| 久久精品影院6| 特级一级黄色大片| 一本久久精品| 99国产极品粉嫩在线观看| 亚洲美女搞黄在线观看| 亚洲内射少妇av| 国产精品乱码一区二三区的特点| 男人的好看免费观看在线视频| 波多野结衣高清无吗| 国产精品久久电影中文字幕| 长腿黑丝高跟| 日韩国内少妇激情av| 禁无遮挡网站| 久久99蜜桃精品久久| 国产在视频线在精品| 91久久精品国产一区二区三区| 少妇高潮的动态图| av女优亚洲男人天堂| 亚洲精品日韩av片在线观看| 观看免费一级毛片| 中国国产av一级| 三级男女做爰猛烈吃奶摸视频| 一进一出抽搐动态| 少妇丰满av| 一级毛片aaaaaa免费看小| 69人妻影院| 国内精品美女久久久久久| 国产69精品久久久久777片| 国产精品,欧美在线| 97超碰精品成人国产| 国模一区二区三区四区视频| 亚洲av中文av极速乱| 亚洲精品乱码久久久v下载方式| 精品久久国产蜜桃| 老熟妇乱子伦视频在线观看| 国产人妻一区二区三区在| 国产不卡一卡二| 99精品在免费线老司机午夜| 久久6这里有精品| 国产女主播在线喷水免费视频网站 | 成人漫画全彩无遮挡| 人体艺术视频欧美日本| 99热全是精品| 91麻豆精品激情在线观看国产| 久久久久久久久中文| 69人妻影院| 99久久久亚洲精品蜜臀av| 一级黄片播放器| 亚洲av.av天堂| 中国美白少妇内射xxxbb| 午夜免费激情av| 日韩人妻高清精品专区| 成人二区视频| 免费不卡的大黄色大毛片视频在线观看 | а√天堂www在线а√下载| 亚洲一区二区三区色噜噜| 国产v大片淫在线免费观看| 长腿黑丝高跟| 婷婷亚洲欧美| 在线免费十八禁| av在线观看视频网站免费| 久久精品国产鲁丝片午夜精品| 嫩草影院精品99| 国内精品一区二区在线观看| 国产精品久久视频播放|