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

    Promoting wetting of Mg on the SiC surfaces by addition of Al,Zn and Zr elements:A study via first-principl calculations

    2022-07-26 11:41:36WeibingGuoWenshnBinHitoXueXiomingZhng
    Journal of Magnesium and Alloys 2022年6期

    Weibing Guo ,Wenshn Bin ,Hito Xue ,Xioming Zhng

    aSchool of Materials Science and Engineering,Hebei University of Technology,Tianjin 300401,China

    b State Key Laboratory of Advanced Welding and Joining,Harbin Institute of Technology,Harbin 150001,China

    Abstract When preparing SiC/Mg composites,alloy elements play key roles in the nucleation and interfacial wetting.In this paper,effects of Al,Zn and Zr additions on the Mg/SiC interfacial bonding properties were investigated comprehensively via method of first-principl calculations.Mg(0001)/SiC(0001) interfaces with different terminations and stacking sequences were built and C Ⅱ-T C Ⅱtop interface has the largest work of adhesion (Wad).Zn dopants can not improve the Wad for both C-T and Si-T interfaces.Al atom can only strengthen the C-T interface.Zr addition can greatly improve the Wad for both C-T and Si-T interfaces.Wad of C-T interfaces can reach up to 11.55 J/m2 and 12.55 J/m2 after doping 1 monolayer (ML) of Al and Zr atoms.Larger Wad can lead to lower contact angles of Mg on SiC surfaces,which can improve wetting and nucleation in SiC/Mg composites.Analysis of electronic structure shows that Al-C and Zr-C bonds have more covalent composition than the Mg-C bond,which is responsible for the improvement of interfacial bonding strength.Experimental results in references were also analyzed,which are in well agreement with our calculation results.

    Keywords: Mg alloy;SiC;Wetting;Electronic structure;First-principle calculations.

    1.Introduction

    Mg alloys have excellent properties of high specifi strength and low density,which make them as potential candidates in automotive,electronic,and aerospace industries[1,2].Recently,Mg and its alloys are also found to be potential materials for energy storages [3,4] due to high gravimetric and volumetric storage capacities.However,Mg alloys also have shortcomings of low strength,poor room temperature toughness and ductility,which have limited the wide applications of magnesium.In conventional processes,alloy elements of Al,Zn,Y,Zr and Ga are always added into Mg alloys to improve their mechanical properties [5].In order to improve specifi strength,creep resistance,thermal shock resistance and corrosion resistance further,some ceramic particles (such as SiC [6],TiC [7],ZrB2[8],Al2O3[9]) are added when preparing Mg matrix composites [10].Among these ceramics,SiC is common used for its properties of high strength,elastic modulus and hardness [11].SiC/pure Mg [12],SiC/AZ80[13],SiC/AZ61 [14,15] SiC/AZ91 [16,17] and SiCw/ZK60[18] composites were prepared by different processes.

    A challenge to prepare SiC/Mg composites is the poor wettability between Mg melt and SiC ceramics.Additions of alloy elements to the SiC/Mg interfaces play key role in wetting and nucleation of the composites,which will influenc the microstructure and performance of composites further.First-principle calculation is based on density functional theory.Results of work of adhesion [19] and electronic structures [20] can be obtained through first-principl calculations,which can help us to understand the interface from atomic or electronic level.Adhesion properties and electronic structures of Mg/TiB2[21],Mg/ZrB2[22],Mg/Al4C3[23] and Mg/Al2CO [24] interfaces were all studied through first-principl calculations.

    Fig.1.The relationships between cut-off energy and total energies of (a) SiC and Mg crystals,(b) surface and interface structures.

    However,little study focused on the influenc mechanism of alloy elements on the SiC/Mg interface,which makes it difficul to design the composition of Mg alloy matrixs.Al[25,26],Zn [27,28] and Zr [29] are the common used alloy elements in Mg alloys and composites.In this work,through method of first-principl calculations,we examine that whether Al,Zn and Zr alloy elements can promote interfacial wetting between Mg alloy and SiC particles.The interfacial bonding properties and electronic structures are analyzed and discussed.The results can help us to optimize alloy compositions for SiC/Mg composites.

    2.Calculation methods

    Code of CASTEP (Cambridge Serial Total Energy Package),which is in framework of DFT (density functional theory),was used to study the energies and properties of the interfaces [30].The generalized gradient approximation (GGA)in Perdew-BurkeErnzerhof (PBE) was used as the exchangecorrelation functional.Method of ultrasoft pseudopotential was used to describe the interactions between ionic core and valence electrons.The valence electronic configuration were 2p63s2for Mg,2s23p2for C,3s23p2for Si,3s23p1for Al,4s24p64d25s2for Zr,and 3d104s2for Zn,respectively.Selection of cut-off energy has influenc on calculation accuracy and computing time.The relationship among cut-off energy and total energies of crystals (Mg and SiC),surfaces (Mg(0001) and SiC(0001)) and interface (Mg/SiC)structures are shown in Fig.1.When the cut-off energy is larger than 400 eV,energies of different structures all change little.Therefore,in the following calculation process,the cutoff energies of the plane waves were all set to be 400 eV.The Brillouin zone was sampled with a Monkhorst-Packk-point mesh.Thek-points meshing were 6×6×6 for the bulks and 6×6×1 for the interface and surface slabs.The tolerance was 1.0×10?5eV/atom when calculating energy.In the relaxation process,the convergence tolerances of stress,displacement and force were 0.05 GPa,0.001 °A and 0.03 eV/°A,respectively.All the calculation parameters were selected based on convergence test to balance the calculation accuracy and computing time.

    Fig.2.Crystal structures (a) Mg,(b) 4H-SiC (0001).

    3.Results and discussion

    3.1.Bulk properties of Mg and 4H-SiC

    Fig.2 shows the structures of Mg and 4H-SiC crystals.Mg crystal has HCP structure (P63/MMC).The lattice constants of Mg are:a=b=3.2094 °A,c=5.2105 °A andα=β=90°,γ=120°.The 4H-SiC crystal has space group of P63MC.The lattice constants of 4H-SiC are:a=b=3.0780 °A,c=10.0460 °A andα=β=90°,γ=120°.

    Selection of exchange-correlation functional can influenc the accuracy of calculation results.Therefore,the lattice constants and bulk modulus (B0) of Mg and SiC crystals were calculated using different exchangecorrelation functionals.Computed lattice constants andB0are listed in Table 1.Lattice constants of Mg using GGA-PBE functional area=b=3.2222 °A,c=5.1707 °A.The calculated results are very close to the experimental values,and also in well agreement with the values in Ref.[22] (a=b=3.222 °A,c=5.170 °A).The calculatedB0is 35.49 GPa,which is slightly lower than that in Ref.[22] (B0=36.3 GPa).Lattice constants using GGA-PW91 are close to that using GGA-PBE,but value ofB0is much lower.Lattice constants using LDA-CAPZ functional are lower (a=b=3.1597 °A,c=5.0766 °A) than that using GGA.The value ofB0(30.87 GPa) is also lower.Lattice constants(a=b=2.9390 °A,c=4.5208 °A) andB0(21.61 GPa) using HSE functional are all the lowest,and also much lower than the experimental values.As to the result of 4H-SiC,the calculated results with GGA-PBE and GGA-PW91 are close to the results in experiments and Ref.[31].Value of lattice constants (a=b=3.0553 °A,c=9.9649 °A) are much lower using HSE functional.In contrasting,value ofB0(224.13 GPa) is the highest,which is also higher than that in Ref.[21].Lattice constants (a=b=3.0414 °A,c=9.9546 °A)andB0(195.17 GPa) using LDA-CAPZ functional are all the lowest.Therefore,functional of GGA-PBE is used in the following calculations to ensure the accuracy.

    Table 1 Computed values of lattice constants and bulk modulus (B0) using different exchange-correlation functional.

    3.2.Surface properties of Mg (0001) and 4H-SiC (0001)

    4H-SiC (0001) polar surfaces are adopted for they can be either C-terminated (C-T) or Si-terminated (Si-T).The lattice constants of the SiC (0001) surface are:a=b=3.0780 °A andγ=120°.Mg (0001) surface is the closed-packed surface,and the lattice constants of Mg (0001) surface slab are:a=b=3.2094 °A andγ=120°.If Mg (0001) and SiC (0001)surfaces are matched together to form an interface,the misfi of interface is as low as 4.09%.

    Firstly,Mg(0001)and SiC(0001)surface slabs were built.For the purpose of building thick enough slab to present bulklike properties,surface slabs with different thickness experienced fully relaxation.3~15 atomic layers of Mg(0001) and SiC(0001) surface slabs were built.Vacuum layer with thickness of 16 °A was added for each surface slab to eliminate the interactions between the surface atoms.For the SiC(0001)surfaces,both C-T and Si-T surfaces were considered.Changes of spacing (Δij) between neighboring layers were monitored,which can be calculated by Eq.(1):

    Where dijand dij0are the interlayer distances of neighboring layers after and before geometry optimization,respectively.

    Fig.3.Surface configuration (a) Mg(0001),(b) C-T SiC(0001),(c) Si-T SiC(0001).

    Results ofΔijafter fully relaxation are listed in Table.2.In general,value ofΔijdecreases with increase of thickness.For the Mg(0001) slabs,large values ofΔijare mainly found between the upper three layers.Δijbecomes relatively low when the slab is thicker than 5 atomic layers.For the CT SiC(0001) surface slab,value ofΔijbetween the second and third layers is always the largest,and it becomes very low when the slabs are thicker than 11 atomic layer.As to the Si-T SiC(0001) surface slab,much larger values ofΔijare found between the upper four atomic layers,and they decrease to lower values when the slab is thicker than 11 atomic layers.To ensure the calculation accuracy,7 atomiclayer Mg(0001) and 13 atomic-layer SiC (0001) surface slabs were built for both C-T and Si-T surfaces.The 13 atomiclayer 4H-SiC (0001) surfaces are non-stoichiometric,which can eliminate the dipole effect on the surfaces.C and Si atoms in the SiC slab have different positions.Therefore,C-T slabs have both C Ⅰ-T and C Ⅱ-T surfaces,Si-T slabs have both SiⅠ-T and Si Ⅱ-T surfaces,which are shown in Fig.3a),b) and c).

    Table 2 Structure relaxation of Mg(0001) and SiC (0001) surfaces.

    Surface energies of Mg(0001) and SiC (0001) surfaces were calculated.For Mg(0001) slab,surface energy (γ) can be calculated by Eq.(2) [32]:

    WhereEslabandΔEare the energies of surface slab and single one atomic layer,Nis the atomic layers andArepresents the area of surface.Calculatedγof Mg(0001) is 0.52 J/m2.Surface energy of SiC (0001) is calculated considering chemical potential (μ) of both Si and C atoms.γcan be calculated by using Eq.(3) [33]:

    Where,Eslabis the total energy of the surface slab,Niandμiare numbers and chemical potentials of the corresponding atoms.Ais the area of surface slab.In general,the terms of PV and TS can be neglected at 0 K.

    The total energy of the bulk SiCcan be given by Eqs.(4) and (5) [34]:

    Where,are the chemical potentials of atoms in the corresponding slab and bulk.ΔHforis the formation enthalpy of the SiC.Value ofΔHforcan be calculated by Eq.(6):

    The chemical potential of atom in the slab is always lower than that in the bulk,therefore:

    Fig.4.Relationship between surface energy of SiC(0001)slabs and chemical potential of C.

    Surface energy of SiC can be calculated by Eq.(9):

    The calculated formation enthalpy of the SiC is 0.68 eV,which is very close to the result (0.66 eV) in Ref.[35].Fig.4 shows the relationship between surface energy of SiC(0001)surface and chemical potential of C.The surface energy is not a constant for the surface slab is either Si-rich or Crich.As increase of chemical potential of C,surface energy of C-T surface decreases from 8.59 J/m2to 7.92 J/m2,and surface energy of Si-T surface increases from 2.86 J/m2to 3.52 J/m2.It is found that surface energy of C-T slab is always larger than that of Si-T slab.The result indicates that C-T surface is more active.The termination of SiC(0001) surface can influenc the Mg/SiC interfacial bonding.

    Fig.5.Interface configurations (a)C Ⅰ-T C Ⅰtop,(b) C Ⅰ-T C Ⅱtop,(c) C Ⅱ-T C Ⅱtop,(d) C Ⅱ-T Si Ⅱtop,(e) Si Ⅰ-T Si Ⅱtop,(f) Si Ⅰ-T Si Ⅰtop,(g) SiⅡ-T C Ⅱtop,(h) Si Ⅱ-T Si Ⅱtop.

    3.3.Mg(0001)/SiC(0001) interfaces

    Mg(0001)/SiC(0001) interfaces with different terminations and stacking sequences were constructed.When the SiC(0001) surfaces are terminated with C Ⅰ,C Ⅱ,Si Ⅰor Si Ⅱatom,the interfaces are named as CⅠ-T,CⅡ-T,SiⅠ-T or Si Ⅱ-T interfaces.For each interface,two kinds of stacking sequences were considered.When the Mg atom was placed on the top site of C Ⅰ,C Ⅱ,Si Ⅰor Si Ⅱatom directly,the interfaces are named as C Ⅰtop,C Ⅱtop,Si Ⅰtop or SiⅡtop interfaces.Therefore,8 kinds of interfaces were built,as shown in Fig.5a)~h).When building interface structure,optimized lattice parameters were used.The softer Mg(0001)slab was strained to match the SiC(0001) slab.For all the interfaces,a 16 °A-thick vacuum layer was added to separate the two free surfaces.

    Work of adhesion (Wad) is used to evaluate the interfacial bonding strength quantitatively.Value of Wadcan be calculated by using Eq.(10) [36]:

    Fig.6.Relationship between Wad and interfacial separation.

    Interfacial distance (d) represents vertical distance between the two slabs.In the calculating process,value of d was adjusted manually from 0 °A to 6 °A,and the corresponding Wadwas calculated.Results of Wadwith various d are shown in Fig.6.The curves all have the same trends.Peak value of Wadcan be found with the equilibrium interfacial distance(d0).In general,Wadof C-T interface is larger than that of Si-T interface.d0of C-T interface is smaller.Smaller value of d0always indicates shorter interfacial bonding length and stronger bonding.For the Si-T interfaces,Si Ⅱ-T C Ⅱtop interface has lowest Wadof 3.70 J/m2and d0of 2.4 °A.As to the C-T interfaces,Wadof C Ⅱ-T C Ⅱtop interface is 7.21 J/m2,which is much larger than that of other two interfaces.d0of C Ⅱ-T C Ⅱtop interface is only 1.3 °A,which is also the smallest in all interfaces.Results of Wadand d0imply that Mg-C bonds are stronger than Mg-Si bonds at the interface.Si Ⅱ-T C Ⅱtop and C Ⅱ-T C Ⅱtop interfaces with largest Wadare used to study the influenc of alloy elements in the following.

    Fig.7.Side views of clean and doped interfaces after relaxation.(a) clean Si-T interface,(b) Si-T interface with Al,(c) Si-T interface with Zn,(d) Si-T interface with Zr,(e) clean C-T interface,(f) C-T interface with Al,(g) C-T interface with Zn,(h) C-T interface with Zr.

    3.4.Mg(0001)/SiC(0001) interfaces with alloy elements

    Al,Zn and Zr atoms are the common used alloy elements in AZ and ZK series Mg alloys and composites.Next,Al,Zn and Zr atoms were introduced to the clean Mg/SiC interfaces by replacing Mg atoms in-situ.In order to investigate the influenc of alloy elements coverage,2×2 supercells of Mg/SiC interfaces were built.Firstly,the selected clean Si-T and C-T interfaces were relaxed.Fig.7a) and e) show the Si-T and C-T interfaces after relaxation.Wadof Si-T and C-T interfaces are 3.76 J/m2and 7.65 J/m2after relaxation,which are both slightly higher than that before geometry optimization.Lengths of Mg-Si and Mg-C bonds at interfaces are 2.926 °A and 2.166 °A.Then one Al (or Zn,Zr) atom was doped into interface,and the interfaces also experienced relaxation.The relaxed Si-T and C-T interfaces with alloy elements are shown in Fig.7b)~d) and f)~h).During the relaxation process,Al dopant moved towards SiC slab at the Si-T interface.Bond length of Al-Si decreases to 2.898 °A after relaxation.While,Zn atom moved away from the SiC slab,and Zn-Si bond has length of 2.994 °A.It is also found that Zr and Si atoms moved towards each other and Zr-Si bond has length of 2.787 °A,which is much smaller than other bonds at interface.As to the doped C-T interfaces,Al,Zn and Zr atoms all moved closer to the C atoms.Bond lengths of Al-C,Zn-C and Zr-C bonds are shortened to 2.006 °A,2.005 °A and 2.082 °A,respectively.

    Fig.8.Wad of clean and doped Mg/SiC interfaces.

    Fig.8 shows Wadvalues of clean and doped Mg/SiC interfaces.The clean Si-T interface has Wadof 3.76 J/m2.After doping Al and Zn atoms,values of Waddecrease to 3.63 J/m2and 3.66 J/m2.The results indicate that Al and Zn dopants would weaken the Si-T interfaces.When the interface is doped with Zr atom,Wadincreases to 4.97 J/m2.The large Wadand short Zr-Si bond length imply that Zr atom can form strong bond with Si atom,which can strengthen the Si-T interface.Wadof clean C-T interface is 7.65 J/m2.Zn dopant has little influenc on the bonding strength of C-T interface.Wadof Zn-doped C-T interface is 7.63 J/m2,which is almost the same as that of clean C-T interface.Al and Zr dopants can both improve Wadof C-T interfaces.Values of Wadincrease to 8.39 J/m2and 9.65 J/m2with Al and Zr atoms.Therefore,we can conclude that Zn atom can not improve interfacial bonding for both C-T and Si-T interfaces.Al atom can only strengthen the C-T interface.Zr atom can strengthen both C-T and Si-T interfaces.The strengthening effect of Zr atom is better than that of Al atom.

    Fig.9.Values of Wad with different coverage of dopant atoms.

    Different numbers of Al and Zr atoms were introduced into the C-T interfaces.Therefore,we can study the influenc of dopant numbers on the interfacial bonding.In the 2×2 C-T interface,each layer contains four atoms.Therefore,1,2,3 and 4 dopant atoms correspond to 0.25,0.5,0.75 and 1 monolayer (ML) of atom coverage.Fig.9 shows the values of Wadwith different ML of Al and Zr atoms.It is found that Wadincreases continually with coverage of dopant atoms.Wadof Zr-doped interface is always larger than that of Al-doped interface.When the interfaces are doped with four atoms,the interface will be saturated with Al-C and Zr-C bonds.Values of Wadincrease to 11.55 J/m2and 12.55 J/m2with 1 ML of Al and Zr atoms,which are much larger than that of clean Mg/SiC interfaces.

    Heat of segregation (ΔGseg) is calculated to evaluate the stability of dopant atoms at the interfaces.The value ofΔGsegcan be obtained by Eq.(11) [37]:

    WhereEMg/SiC:nAandEMg/SiCrepresent total energies of the interfaces with dopant atoms and clean Mg/SiC interface.μand n are chemical potential and number of Al or Zr atoms.It is reported that more negative value ofΔGsegrepresents more stability in thermodynamics [38].

    Fig.10 showsΔGsegvalues with different coverage of dopant atoms.When C-T interfaces are doped with one Al and Zr atom,ΔGsegare ?2.02 eV and ?3.98 eV.The negative values indicate that doping Al and Zr atoms to the CT Mg/SiC interface is stable in thermodynamics.The lowerΔGsegimplies that doping Zr atom is more stable.ΔGsegdecreases with increase of Al atom coverage,and increases with Zr atom coverage.Values ofΔGsegare all negative with different atom coverage.When the C-T interfaces are doped with 1 ML of Al and Zr atoms,ΔGsegare ?2.22 eV and?3.17 eV.The results indicate that the interfaces are still stable with four doped Al and Zr atoms.

    Fig.10.Relationships between ΔGseg and coverage of dopant atoms.

    Fig.11.Charge density difference of Mg/SiC interfaces (a) clean interface,(b) Al-doped interface,(c) Zr-doped interface.

    3.5.Electronic structures

    To have a deeper insight on the bonding mechanism of the doped interfaces,electronic structures are analyzed through methods of charge density difference and partial density of states(PDOS).Results of charge density difference are shown in Fig.11.In the plots,the charge enrichment area is labeled as red color and charge dissipation area is labeled as blue color.Fig.11a) shows plot of charge density difference for clean interface.Red and blue areas are found around C and Mg atoms,which indicates that some electrons are transferred from Mg to C atoms.Mg and C atoms become cation and anion.They bond with each other through electrostatic attraction,which implies ionic bond.Fig.11b) and c) show charge density difference plots of Al and Zr-doped interfaces.Charge transfers are all found from the Al and Zr atoms to the neighboring C atoms.Therefore,Mg-C,Al-C and Zr-C bonds at the interface all have compositions of ionic bonds.Al-C and Zr-C bonds are stronger than the Mg-C bond for more transferred charges.

    Fig.12.(a) totals DOS results of SiC and Mg/SiC interface,PDOS results of (b) clean Mg/SiC interface,(c) Al-doped interface,(d) Zr-doped interface.

    Total DOS results of SiC and Mg/SiC interface are shown in Fig.12a).PDOS results of selected atoms in clean and doped Mg/SiC interface are shown in Fig.12b)~d).In the clean interface (Fig.12b),Mg-spand C-sporbitals overlap from ?9 eV to 1 eV.Overlapping of Mg-spand C-sporbital is also found within the energy range of ?12~?9.7 eV.These peaks are narrow,which represents localization of valence electrons.This kind of overlapping indicates covalent bond.Therefore,Mg-C bond at the interface is a mix of covalent and ionic bond.Bond population analysis is also performed to evaluate the degree of covalent.Bond population of Mg-C is 0.20,indicating that Mg-C bond is mainly ionic.Fig.12c)and d) show PDOS of atoms at Al and Zr-doped interfaces.The similar overlapping between Al(or Zr)and neighboring C atom is found,which indicates that Al-C and Zr-C bonds also have both ionic and covalent compositions.Bond populations of Al-C and Zr-C bonds are 0.93 and 1.02,which suggests that these two bonds are mainly covalent.

    3.6.Discussion

    SiC exists as form of solid particles in SiC/Mg composites for SiC has high melting point and high stability.Mg and 4H-SiC crystals all have hexagonal structure,and the lattice mismatch of Mg(0001)/SiC(0001) is only 4.09%.Therefore,SiC can be heterogeneity nucleus for Mg.When preparing SiC/Mg composites,wetting of Mg melt on SiC particles is the key factor to influenc the quality and properties of composites.Doping of Al and Zr atoms to the interface can raise Wadvalue,which also has influenc on wetting.The relationship between Wadand contact angle (θ) can be expressed as[39]:

    WhereγLVis liquid-vapor interfacial tension.

    When preparing SiC/Mg composites,γLVof Mg-Al (or Mg-Zr) alloy/SiC interface is a constant.Al and Zr atoms would diffuse to the interface and lead to larger Wad.The contact angle will be smaller with larger Wad.Therefore,segregation of Al and Zr to the interfaces can promote wetting and nucleation.

    The calculation results are compared with the experimental results in the references.In Ref.[40],SiCw/AZ91 composites were prepared by squeeze cast at 800 °C.Massive Mg17Al12intermetallic compounds (IMC) are found at the interface between Mg alloy and SiC whiskers.The results confir segregation of Al to the Mg/SiC interface.Al would stay at the interface as form of Mg17Al12IMC.It also suggests that the IMC preferentially nucleates at the matrix-whisker interface for the strong bonds between Al and C atoms.In Ref.[41],SiC-reinforced magnesium alloy metal matrix composites are produced by rheocasting process.STEM micrograph as well as STEM/EDX maps are carried out to analyze the interface of SiC and Mg alloys.A lot of particles were found close to the SiC/metal interface,and the particles were rich in Al and C.Composition analysis shows the stoichiometry corresponds to AlC2.The results confir the strong bonds between Al and C,and Al-C compounds can form under certain conditions.In Ref.[42],SiCw/ZK51A magnesium matrix composite was fabricated by two-step squeeze casting.Through the TEM observation,SiCw/Mg interface was found to be quite clean,indicating no IMC formed during the fabrication process of the composite.It is also mentioned the SiC whiskers are very stable and carbides are difficul to form at the interface.From the results,we notice that although Zr can form strong bonds with C atoms,no Zr-C compounds formed at the interface.Mg alloys containing Zr can be excellent candidates for the matrix of composites.

    4.Conclusions

    In this paper,effects of Al,Zn and Zr atoms on the Mg(0001)/SiC(0001) interfacial bonding properties were investigated via methods of first-principl calculations.The C-T interfaces were found to have larger Wadthan the Si-T interfaces.C Ⅱ-T C Ⅱtop interface has largest Wadof 7.21 J/m2.Then Al,Zn and Zr atoms were doped into the Mg/SiC interfaces.Zn atom can not improve interfacial bonding for both C-T and Si-T interfaces.Al atom can only strengthen the C-T interface.Zr atom can strengthen both C-T and Si-T interfaces.Wadof Zr-doped C-T interface is improved to 9.65 J/m2.In the Al and Zr-doped C-T interfaces,Wadincreases as the ML of atom coverage.Wadcan reach up to 11.55 J/m2and 12.55 J/m2with 1 ML of Al and Zr atoms.Increase of Wadcan lead to lower contact angles,which can improve wetting and nucleation.NegativeΔGsegvalues imply that doping Al and Zr to the C-T interface is stable in thermodynamics.Al-C and Zr-C bonds formed at interfaces have composition of both ionic and covalent.Al-C and Zr-C bonds have more covalent composition than the Mg-C bond,which is responsible for strong bonding strength.Experimental results in references were also analyzed,which can confir the calculation results.

    Conflic of Interest

    None.

    Acknowledgments

    This research was funded by China Postdoctoral Science Foundation (No.2020M670670),Foundation Strengthening Program (No:2019-JCJQ-00),Natural Science Foundation of Hebei Province (E2019202407),Science and Technology Research Project of Hebei Province Colleges and Universities(QN2019028) and“Yuan Guang Scholar”Plan of Hebei University of Technology.

    99国产综合亚洲精品| 啦啦啦中文免费视频观看日本| 春色校园在线视频观看| 丁香六月天网| 波野结衣二区三区在线| 国产成人免费观看mmmm| 好男人视频免费观看在线| 在线亚洲精品国产二区图片欧美| 亚洲中文av在线| 街头女战士在线观看网站| 九色成人免费人妻av| 久久精品国产综合久久久 | 日韩伦理黄色片| 国产亚洲最大av| 男人爽女人下面视频在线观看| 一个人免费看片子| 精品熟女少妇av免费看| 成人国产av品久久久| 菩萨蛮人人尽说江南好唐韦庄| 看免费av毛片| 久久人人97超碰香蕉20202| 久久av网站| 亚洲av国产av综合av卡| 国产综合精华液| 一边摸一边做爽爽视频免费| 考比视频在线观看| 1024视频免费在线观看| 夫妻午夜视频| 国产精品久久久久久av不卡| 日本欧美国产在线视频| 久久久国产精品麻豆| 国产女主播在线喷水免费视频网站| 日本色播在线视频| 欧美精品国产亚洲| 高清不卡的av网站| av女优亚洲男人天堂| 婷婷色av中文字幕| 国产精品麻豆人妻色哟哟久久| 久久精品国产鲁丝片午夜精品| tube8黄色片| 日韩视频在线欧美| 少妇被粗大的猛进出69影院 | 成人毛片60女人毛片免费| 在线观看免费日韩欧美大片| 一区在线观看完整版| 哪个播放器可以免费观看大片| 成人国语在线视频| 搡老乐熟女国产| 精品国产国语对白av| 曰老女人黄片| 精品99又大又爽又粗少妇毛片| 久热这里只有精品99| 午夜福利网站1000一区二区三区| 夫妻午夜视频| 激情视频va一区二区三区| 日韩成人伦理影院| 欧美日韩综合久久久久久| 日本av免费视频播放| 黄片播放在线免费| 三上悠亚av全集在线观看| 在线天堂中文资源库| 精品国产乱码久久久久久小说| 一区二区三区四区激情视频| 丝瓜视频免费看黄片| 久久亚洲国产成人精品v| 成人免费观看视频高清| 视频区图区小说| 少妇的丰满在线观看| 精品久久久精品久久久| 亚洲国产毛片av蜜桃av| 天天躁夜夜躁狠狠久久av| 免费av中文字幕在线| 久久久久国产精品人妻一区二区| 久久精品国产亚洲av涩爱| 国产精品一国产av| 免费观看性生交大片5| 女人被躁到高潮嗷嗷叫费观| 久久久久精品人妻al黑| 日韩制服骚丝袜av| 精品一区二区三区视频在线| 色5月婷婷丁香| 超色免费av| 青青草视频在线视频观看| 超碰97精品在线观看| 国产探花极品一区二区| 少妇人妻久久综合中文| 国产精品久久久久久精品电影小说| 亚洲丝袜综合中文字幕| 国产一区二区三区av在线| 最近中文字幕高清免费大全6| 美女视频免费永久观看网站| 国产精品.久久久| 精品久久久久久电影网| 亚洲av免费高清在线观看| 美女中出高潮动态图| 国产精品 国内视频| 男女高潮啪啪啪动态图| 欧美 日韩 精品 国产| 久久精品人人爽人人爽视色| 日日撸夜夜添| 最近的中文字幕免费完整| 亚洲欧美成人精品一区二区| 五月伊人婷婷丁香| 亚洲精品日本国产第一区| 日韩欧美一区视频在线观看| 国产精品国产三级专区第一集| 9色porny在线观看| 伊人久久国产一区二区| 一级片免费观看大全| 欧美精品av麻豆av| 丝瓜视频免费看黄片| 肉色欧美久久久久久久蜜桃| 久久久久人妻精品一区果冻| 老熟女久久久| 1024视频免费在线观看| 丝袜美足系列| 国产xxxxx性猛交| 日韩人妻精品一区2区三区| 搡老乐熟女国产| 香蕉精品网在线| 成人亚洲欧美一区二区av| 18禁裸乳无遮挡动漫免费视频| 丝袜脚勾引网站| 亚洲欧美成人精品一区二区| 久久久久久伊人网av| 免费日韩欧美在线观看| 午夜免费男女啪啪视频观看| 一区二区av电影网| 在线免费观看不下载黄p国产| 人人妻人人爽人人添夜夜欢视频| 久久久久视频综合| 欧美丝袜亚洲另类| 永久网站在线| 视频中文字幕在线观看| 免费少妇av软件| 亚洲av中文av极速乱| 人人妻人人添人人爽欧美一区卜| 国产高清三级在线| 亚洲精品,欧美精品| 精品久久国产蜜桃| 亚洲欧美中文字幕日韩二区| 97超碰精品成人国产| 中文字幕另类日韩欧美亚洲嫩草| 日本猛色少妇xxxxx猛交久久| 麻豆精品久久久久久蜜桃| 岛国毛片在线播放| 有码 亚洲区| 妹子高潮喷水视频| 亚洲欧美精品自产自拍| 久久青草综合色| 极品少妇高潮喷水抽搐| 中文欧美无线码| 少妇猛男粗大的猛烈进出视频| 少妇精品久久久久久久| av免费观看日本| 午夜福利影视在线免费观看| 亚洲精品一区蜜桃| 免费高清在线观看日韩| 妹子高潮喷水视频| 麻豆精品久久久久久蜜桃| 九九爱精品视频在线观看| 亚洲成色77777| av片东京热男人的天堂| 最近手机中文字幕大全| 亚洲图色成人| 人人妻人人澡人人爽人人夜夜| 亚洲美女视频黄频| 精品久久久精品久久久| 久久久久久久亚洲中文字幕| 男女下面插进去视频免费观看 | 国产精品嫩草影院av在线观看| 巨乳人妻的诱惑在线观看| 天天操日日干夜夜撸| 国产国拍精品亚洲av在线观看| 亚洲综合精品二区| 亚洲国产毛片av蜜桃av| 亚洲精品aⅴ在线观看| 少妇的逼好多水| a 毛片基地| h视频一区二区三区| 寂寞人妻少妇视频99o| 国产男人的电影天堂91| 中文乱码字字幕精品一区二区三区| 国产亚洲精品第一综合不卡 | 精品国产国语对白av| 高清毛片免费看| 性高湖久久久久久久久免费观看| 亚洲天堂av无毛| 国产精品三级大全| 欧美xxxx性猛交bbbb| 精品人妻在线不人妻| 母亲3免费完整高清在线观看 | 亚洲国产av新网站| 一级爰片在线观看| 美女脱内裤让男人舔精品视频| 亚洲av成人精品一二三区| 99久久综合免费| 韩国高清视频一区二区三区| 国产免费又黄又爽又色| 男女免费视频国产| 各种免费的搞黄视频| 亚洲精品久久午夜乱码| 下体分泌物呈黄色| 午夜福利视频在线观看免费| 十八禁网站网址无遮挡| 亚洲色图 男人天堂 中文字幕 | 国产一级毛片在线| 在线精品无人区一区二区三| 哪个播放器可以免费观看大片| 国产成人av激情在线播放| 免费看av在线观看网站| 久久精品国产亚洲av涩爱| 亚洲天堂av无毛| 国产av国产精品国产| 18禁动态无遮挡网站| 国产在线视频一区二区| 又大又黄又爽视频免费| 中文乱码字字幕精品一区二区三区| 精品人妻熟女毛片av久久网站| 国产精品国产三级国产av玫瑰| 建设人人有责人人尽责人人享有的| 久久久久久久久久人人人人人人| 精品第一国产精品| 一级爰片在线观看| 少妇人妻久久综合中文| 亚洲国产成人一精品久久久| 大香蕉久久成人网| 又黄又粗又硬又大视频| 波野结衣二区三区在线| 一区二区三区四区激情视频| 免费在线观看完整版高清| 一本—道久久a久久精品蜜桃钙片| 久久精品久久久久久久性| 国产免费福利视频在线观看| 久久 成人 亚洲| av视频免费观看在线观看| 菩萨蛮人人尽说江南好唐韦庄| 秋霞在线观看毛片| 国产福利在线免费观看视频| 最近手机中文字幕大全| 免费高清在线观看日韩| 久久 成人 亚洲| 国产又色又爽无遮挡免| 亚洲内射少妇av| 久久精品国产亚洲av天美| 这个男人来自地球电影免费观看 | 亚洲欧美日韩卡通动漫| 蜜桃在线观看..| 国产精品国产三级国产av玫瑰| 亚洲精品456在线播放app| 成人毛片a级毛片在线播放| 精品卡一卡二卡四卡免费| 熟女av电影| 国产成人精品婷婷| 国产一区有黄有色的免费视频| 欧美人与性动交α欧美精品济南到 | 亚洲综合色网址| 狠狠精品人妻久久久久久综合| 亚洲激情五月婷婷啪啪| 国产精品久久久久久精品电影小说| 90打野战视频偷拍视频| 天美传媒精品一区二区| 国产av一区二区精品久久| 99热网站在线观看| 熟妇人妻不卡中文字幕| 日本午夜av视频| 久久久久国产精品人妻一区二区| 亚洲av欧美aⅴ国产| av国产久精品久网站免费入址| 午夜福利,免费看| 一边亲一边摸免费视频| 一个人免费看片子| av免费在线看不卡| 极品少妇高潮喷水抽搐| 亚洲国产av影院在线观看| 丰满饥渴人妻一区二区三| 亚洲av日韩在线播放| av卡一久久| 国产毛片在线视频| 国产精品一区二区在线不卡| 欧美变态另类bdsm刘玥| 日本猛色少妇xxxxx猛交久久| 乱码一卡2卡4卡精品| 久久精品国产自在天天线| 一边摸一边做爽爽视频免费| 男女边摸边吃奶| 欧美 亚洲 国产 日韩一| 国产成人精品福利久久| 国精品久久久久久国模美| 日韩成人伦理影院| 成人18禁高潮啪啪吃奶动态图| 日本欧美视频一区| 五月玫瑰六月丁香| 国产一区有黄有色的免费视频| 日韩不卡一区二区三区视频在线| 18禁观看日本| 美女主播在线视频| 18禁动态无遮挡网站| 超色免费av| 国产精品熟女久久久久浪| 人妻一区二区av| 久久久久国产网址| 深夜精品福利| 国产欧美日韩综合在线一区二区| 久久韩国三级中文字幕| 久久久久久人妻| 欧美日韩综合久久久久久| 国产男女内射视频| 久久久精品区二区三区| 久久久国产精品麻豆| 午夜激情久久久久久久| 超碰97精品在线观看| 在线天堂最新版资源| 在线观看国产h片| 精品99又大又爽又粗少妇毛片| 国产免费福利视频在线观看| 成人漫画全彩无遮挡| 欧美成人午夜免费资源| 久久99一区二区三区| 午夜免费鲁丝| 精品视频人人做人人爽| 精品久久国产蜜桃| 赤兔流量卡办理| 亚洲内射少妇av| 国产免费又黄又爽又色| 一级片'在线观看视频| av在线老鸭窝| av又黄又爽大尺度在线免费看| 高清在线视频一区二区三区| 极品人妻少妇av视频| 中文字幕免费在线视频6| 黄片无遮挡物在线观看| 建设人人有责人人尽责人人享有的| 男女啪啪激烈高潮av片| 久久久久国产精品人妻一区二区| 这个男人来自地球电影免费观看 | 久久人人爽av亚洲精品天堂| 国产成人一区二区在线| 视频中文字幕在线观看| 国产一区有黄有色的免费视频| 插逼视频在线观看| 久久毛片免费看一区二区三区| 亚洲av欧美aⅴ国产| 成人国产av品久久久| 亚洲精品国产av蜜桃| 国产成人精品无人区| 人人妻人人添人人爽欧美一区卜| 在线观看一区二区三区激情| 日本与韩国留学比较| 色视频在线一区二区三区| 亚洲精品456在线播放app| 精品酒店卫生间| 精品久久久久久电影网| 美国免费a级毛片| 国产毛片在线视频| 少妇的丰满在线观看| 悠悠久久av| 电影成人av| 日韩免费高清中文字幕av| 午夜福利一区二区在线看| 在线观看日韩欧美| 欧美精品亚洲一区二区| 久久久久精品人妻al黑| 久久久久国产一级毛片高清牌| 久久久久久免费高清国产稀缺| 露出奶头的视频| 久久人人爽av亚洲精品天堂| 精品一品国产午夜福利视频| 日本wwww免费看| 黄色成人免费大全| 亚洲av第一区精品v没综合| 啦啦啦在线免费观看视频4| 美女午夜性视频免费| 欧美日韩一级在线毛片| 亚洲精品一卡2卡三卡4卡5卡| 国产精品自产拍在线观看55亚洲 | 精品视频人人做人人爽| 国产97色在线日韩免费| 午夜两性在线视频| 91av网站免费观看| 国产精品一区二区免费欧美| 亚洲全国av大片| 成人手机av| 精品卡一卡二卡四卡免费| 久久亚洲精品不卡| 欧美日韩福利视频一区二区| 欧美精品亚洲一区二区| 女性生殖器流出的白浆| 国内久久婷婷六月综合欲色啪| 亚洲全国av大片| 别揉我奶头~嗯~啊~动态视频| 亚洲色图av天堂| 亚洲av日韩精品久久久久久密| 法律面前人人平等表现在哪些方面| 成人国产一区最新在线观看| 最新的欧美精品一区二区| 好看av亚洲va欧美ⅴa在| 国产一卡二卡三卡精品| 在线av久久热| 中文字幕制服av| 亚洲精华国产精华精| 韩国av一区二区三区四区| 老汉色av国产亚洲站长工具| 国产激情久久老熟女| 国产日韩欧美亚洲二区| tocl精华| 国产精品1区2区在线观看. | 老司机午夜福利在线观看视频| 18禁观看日本| 亚洲国产欧美一区二区综合| 天天躁日日躁夜夜躁夜夜| 亚洲精品美女久久久久99蜜臀| 亚洲第一欧美日韩一区二区三区| 亚洲国产欧美日韩在线播放| 久久九九热精品免费| 精品福利观看| 中文字幕人妻丝袜制服| 男男h啪啪无遮挡| 国产一卡二卡三卡精品| 一级片'在线观看视频| 精品人妻1区二区| 在线观看日韩欧美| 欧美亚洲日本最大视频资源| 亚洲av日韩在线播放| 亚洲 欧美一区二区三区| 国产又色又爽无遮挡免费看| 一进一出抽搐动态| 黄色成人免费大全| 亚洲熟妇熟女久久| 欧美日韩亚洲国产一区二区在线观看 | 悠悠久久av| 国产av又大| 在线观看免费视频日本深夜| 露出奶头的视频| 我的亚洲天堂| 日本精品一区二区三区蜜桃| 成人18禁在线播放| 女人高潮潮喷娇喘18禁视频| 一区在线观看完整版| 成年版毛片免费区| 午夜福利在线免费观看网站| 国产成人av激情在线播放| 久久中文看片网| 久久热在线av| 国产高清激情床上av| 18禁黄网站禁片午夜丰满| 丰满饥渴人妻一区二区三| 天堂中文最新版在线下载| 看免费av毛片| 欧美激情 高清一区二区三区| 亚洲精品国产一区二区精华液| 韩国精品一区二区三区| 桃红色精品国产亚洲av| 精品久久久精品久久久| 五月开心婷婷网| 国产aⅴ精品一区二区三区波| 一级黄色大片毛片| 狠狠婷婷综合久久久久久88av| 精品卡一卡二卡四卡免费| 在线天堂中文资源库| 久久影院123| 国产精品香港三级国产av潘金莲| 精品午夜福利视频在线观看一区| 久久国产精品男人的天堂亚洲| 99久久国产精品久久久| 欧美日韩黄片免| 国产精品偷伦视频观看了| 老司机亚洲免费影院| av有码第一页| 久久人妻av系列| 久久狼人影院| 亚洲精品国产一区二区精华液| 十八禁高潮呻吟视频| 人人妻人人澡人人看| 欧美亚洲日本最大视频资源| 超色免费av| 国产精品免费视频内射| 亚洲五月天丁香| 又紧又爽又黄一区二区| 一级,二级,三级黄色视频| 777久久人妻少妇嫩草av网站| 91av网站免费观看| 91九色精品人成在线观看| 人妻一区二区av| 欧美日韩亚洲国产一区二区在线观看 | 在线观看66精品国产| 男女下面插进去视频免费观看| 欧美人与性动交α欧美软件| 国产成人精品在线电影| 一级毛片高清免费大全| 亚洲精品国产精品久久久不卡| 看片在线看免费视频| 久久久久久久午夜电影 | 亚洲人成电影免费在线| 国产一区二区三区视频了| 在线观看免费日韩欧美大片| 水蜜桃什么品种好| a级毛片在线看网站| 亚洲va日本ⅴa欧美va伊人久久| 丝袜美腿诱惑在线| 丝袜美足系列| 亚洲av日韩精品久久久久久密| 午夜91福利影院| 99精国产麻豆久久婷婷| 王馨瑶露胸无遮挡在线观看| 国产欧美亚洲国产| 热re99久久国产66热| 久久国产精品大桥未久av| 少妇猛男粗大的猛烈进出视频| 19禁男女啪啪无遮挡网站| 精品一区二区三区av网在线观看| 亚洲精品一二三| 校园春色视频在线观看| 欧美在线黄色| 亚洲第一欧美日韩一区二区三区| 亚洲少妇的诱惑av| 99热网站在线观看| 午夜免费成人在线视频| 18禁黄网站禁片午夜丰满| 一进一出抽搐gif免费好疼 | 亚洲国产欧美一区二区综合| 成年人午夜在线观看视频| 亚洲九九香蕉| 真人做人爱边吃奶动态| 亚洲av日韩精品久久久久久密| 色综合婷婷激情| 大陆偷拍与自拍| 在线观看舔阴道视频| 少妇裸体淫交视频免费看高清 | 99精品久久久久人妻精品| 999久久久精品免费观看国产| 这个男人来自地球电影免费观看| 水蜜桃什么品种好| 久久久久久久久久久久大奶| 一区二区三区国产精品乱码| 久久香蕉国产精品| 国产成人精品久久二区二区免费| 国产三级黄色录像| 亚洲国产欧美网| 亚洲国产精品一区二区三区在线| 手机成人av网站| 满18在线观看网站| 久久久久国产精品人妻aⅴ院 | 亚洲一区中文字幕在线| 18在线观看网站| 国产精品综合久久久久久久免费 | 欧美亚洲日本最大视频资源| 久久九九热精品免费| 精品国产超薄肉色丝袜足j| 成人手机av| 国产成人欧美| 久久久国产成人免费| 亚洲全国av大片| 精品久久蜜臀av无| 亚洲精品久久成人aⅴ小说| 午夜久久久在线观看| 国产免费av片在线观看野外av| 亚洲精品国产色婷婷电影| 免费在线观看完整版高清| 50天的宝宝边吃奶边哭怎么回事| 99在线人妻在线中文字幕 | 啦啦啦 在线观看视频| 丝袜美足系列| 俄罗斯特黄特色一大片| 一本综合久久免费| 欧美激情高清一区二区三区| av不卡在线播放| 亚洲欧美一区二区三区黑人| 国产精品国产高清国产av | 欧美人与性动交α欧美软件| 中文字幕人妻丝袜制服| 女性被躁到高潮视频| 高清黄色对白视频在线免费看| 国产精品99久久99久久久不卡| 中文欧美无线码| 女人高潮潮喷娇喘18禁视频| 高清av免费在线| 国产欧美日韩一区二区三| 亚洲 欧美一区二区三区| 亚洲人成77777在线视频| 亚洲va日本ⅴa欧美va伊人久久| 欧美精品人与动牲交sv欧美| 热re99久久精品国产66热6| av福利片在线| 免费看a级黄色片| 精品欧美一区二区三区在线| 一本一本久久a久久精品综合妖精| 啦啦啦在线免费观看视频4| 中文字幕人妻丝袜制服| av免费在线观看网站| 久久性视频一级片| 伊人久久大香线蕉亚洲五| 久久久国产成人精品二区 | 国产男靠女视频免费网站| 99久久综合精品五月天人人| 18禁国产床啪视频网站| 窝窝影院91人妻| 国产精品av久久久久免费| 在线观看66精品国产| 99精品久久久久人妻精品| 老熟女久久久| 亚洲综合色网址| 18禁国产床啪视频网站| xxxhd国产人妻xxx| 露出奶头的视频| 老司机亚洲免费影院| 亚洲精华国产精华精| 久久婷婷成人综合色麻豆| 少妇的丰满在线观看| 久久精品国产清高在天天线| 窝窝影院91人妻| 性色av乱码一区二区三区2| 一级黄色大片毛片| 日本一区二区免费在线视频| 国产精品香港三级国产av潘金莲|