Fnglei Tong ,Xize Chen ,Shnghi Wei,? ,Jenny Mlmstr?m,b ,Joseph Vell ,Wei Go

aDepartment of Chemical &Materials Engineering,Faculty of Engineering,The University of Auckland,New Zealand

b MacDiarmid Institute for Advanced Materials and Nanotechnology,Wellington,New Zealand

cSchool of Chemical Science,Faculty of Science,The University of Auckland,New Zealand

Abstract Four Mg-xZn-ySn(x=2,4 and y=1,3 wt.%)alloys are investigated as anode materials for magnesium-air(Mg-air)battery.The self-corrosion and battery discharge behavior of these four Mg-Zn-Sn alloys are analyzed by electrochemical measurements and Mg-air battery tests.The results show that addition of Sn stimulates the electrochemical activity and significantly improves the anodic efficiency and specific capacity of Mg-Zn alloy anodes.Among the four alloy anodes,Mg-2Zn-3Sn(ZT23)shows the best battery discharge performance at low current densities(≤5 mA cm?2),achieving high energy density of 1367 mWh g?1 at 2 mA cm?2.After battery discharging,the surface morphology and electrochemical measurement results illustrate that a ZnO and SnO/SnO2 mixed film on alloy anode surface decreases self-corrosion and improves anodic efficiency during discharging.The excessive intermetallic phases lead to the failure of passivation films,acting as micro-cathodes to accelerate self-corrosion.

Keywords: Magnesium alloys;Alloy anode;Self-corrosion;Magnesium-air battery;Discharge performance.

1.Introduction

Metal-air batteries contain a metal anode,an air cathode and electrolyte,have been regarded as promising electricity storage systems for future energy storage because of their low cost and high specific density [1].Additionally,the cathode(O2)comes from the air,which reduces the total weight of the battery [2,3].The Mg-air battery is a type of metal-air battery,which is being considered as next-generation energy storage technology for future energy demand.

Mg anodes have the advantages of a very negative potential(?2.37 V vs.standard hydrogen electrode,SHE),a high Faradic capacity(2.2 Ah g?1),and a low density(1.74 g cm?3)[4,5].Mg-air batteries have a high theoretical voltage(3.09 V)and energy density(6.8 kWh kg?1)[6].Although the Mg-air battery is a primary battery,the consumed Mg anode can be replaced by a fresh Mg anode to make Mg-air battery "refuelable" [7,8].However,the metal anode of Mg-air batteries suffers from low anodic efficiency due to its high self-corrosion rate,and the discharge products tend to accumulate on the electrode,suppressing the electrochemical reaction kinetics between Mg anode and aqueous electrolyte.Much effort has been attempted to reduce the corrosion rate by changing the composition of Mg alloy anode [9–11],and by adding corrosion inhibitors in the electrolyte [7,12–15].

Introduction of alloying elements into Mg is one of the most efficient methods to improve its battery performance.Alloying elements including Al,Zn,Pb,Sn,In,Li,Ca and RE(rare earth elements)influence the microstructure and chemical properties of Mg alloy,which could modify the kinetics of anodic or cathodic [16].Both Zn and Sn are eco-friendly elements with high solid solubility in Mg,6.2 and 14.5 wt.%,respectively.A small amount of Zn(<5 wt.%)can refine the grain size,while high Zn concentration in Mg forms MgxZnyphases distributed as a network structure along the dendrite and grain boundaries [17].The MgxZnyphases provide an age-hardening response [18],but large MgxZnyphases can act as the cathode to the Mg matrix and increase the corrosion rate [19,20].

Previous studies have shown that a small amount of Zn can improve the corrosion resistance of Mg alloy and form a protective film [21,22].During battery discharging,Zn can reduce the pH value of the electrolyte close to the surface,and accelerate the dissolution of discharge products Mg(OH)2[23].In our previous research,the environmental-friendly and low-cost Mg-Zn binary alloys have been studied as anode materials for a Mg-air battery [24].Mg-2Zn alloy anode shows a high utilization efficiency and specific capacity of 54.42%and 1185.50 mAh g?1at 10 mA cm?2,respectively [24].Mg-Zn based ternary alloys have also shown interesting anodic performance,for example,Chen et al.[25]reported that Mg-6Zn-1Y(ZW61)has a high discharge capacity of 1162.8 mAh g?1and anodic efficiency of 55.1% at the current density of 40 mA cm?2.

Sn has been used as a high-performance anode material for Mg-ion rechargeable batteries due to its high theoretical specific capacity(903 mAh g?1)and low standard H2potential(?0.1375 V)[26,27].We recently reviewed the performance of Mg alloys as anodes for Mg-air batteries with aqueous electrolyte system [16].The secondary phase Mg2Sn accelerates the dissolution of Mg anode and breaks the passive film,significantly improving the discharge performance of the Mg-air battery [28].Gu et al.reported that Mg-4Sn-1Zn-1Ag(TZQ411)had good anodic efficiency of 69.53%at 120 mA cm?2[11].Although these studies reported that addition of Sn improves the discharge performance of Mg alloy anodes in the Mg-air battery [16,28],the electrochemical behavior and battery properties of ternary Mg-Zn-Sn alloys have not been systematically studied.In this work,the microstructure,electrochemical properties,and battery discharge performance of four ZT alloys have been studied as the anode materials for Mg-air aqueous batteries.The effect of MgxZnyand Mg2Sn secondary phases on Mg-air battery performance are discussed.

2.Experimental

2.1.Materials

Four Mg-xZn-ySn alloys(x=2,4 andy=1,3 wt.%)were prepared in a vacuum induction furnace.High purity raw Mg,Zn and Sn(>99.9%)materials were melted in a mild steel crucible under argon atmosphere at 720 °C.The melt was stirred well and homogenized before pouring into a steel mould,preheated to 250°C.For comparison,Mg-2Zn and Mg-4Zn alloys were also prepared via the same casting procedures.The nominal and analyzed compositions of ascast Mg-Zn and Mg-Zn-Sn(ZT)alloys are presented in Table 1,showing that the alloys have very low concentrations of other elements,Fe,Ni and Cu.

Table 1 Chemical compositions of Mg-Zn and ZT alloys(analyzed by Agilent 7700 ICP-MS).

Table 2 Electrochemical parameters of Mg-Zn and ZT alloys obtained from polarization curves.

2.2.Microstructure characterization

The microstructure of ZT alloys were examined using optical microscopy(OM,Nikon Eclipse LV100ND),scanning electron microscopy(SEM,Hitachi SU70),and energy dispersive X-ray spectroscopy(EDS,Thermo Scientific Noran System 7).Metallographic specimens of as-cast ZT alloys were ground on SiC papers of grit sizes up to 1200# and cleaned with distilled water.In the final polishing stage,all specimens were carefully polished by 0.05 μm silica in ethanol to obtain high-quality surfaces for microstructure observation.The samples were etched with a solution consisting of 5 mL acetic acid,6 g picric acid,10 mL H2O,and 100 mL ethanol.

2.3.Electrochemical measurements

All Mg alloy samples were ground with abrasive papers up to 1200 grit before electrochemical measurements and battery testing.A typical three electrodes system was used for electrochemical analysis by an electrochemical workstation(BioLogic SP300):Pt plate as a counter electrode,saturated calomel electrode(SCE)as a reference electrode and Mg alloy samples as working electrodes in 3.5 wt.%NaCl solution.Potentiodynamic polarization analysis was conducted at a scan rate of 1 mV s?1after 10 min immersion in the electrolyte.The open-circuit potential(OCP)was recorded in a 3.5 wt.%NaCl solution for 30 min waiting time to achieve a stable surface.Electrochemical impedance spectroscopy(EIS)was conducted at OCP with frequencies of 100 kHz to 0.05 Hz and 5 mV voltage amplitude.The potentiodynamic polarization and EIS tests were performed more than two times to ensure the repeatability of the data.In addition,Atomic Force Microscope(MFP-3D Origin AFM)combined with Scanning Kelvin Probe force microscopy(SKPFM)was used to study the volta potential difference between intermetallic compounds and Mg matrix.The silicon probe is coated with a Cr/Pt conductive coating(Budget Sensors Electritap300-G,resonant frequency 300 kHz).

2.4.Battery discharge and performance calculation

The battery discharge properties were tested by a LAND battery test system(CT3001A)in a Mg-air battery test kit at four current densities of 1,2,5,and 10 mA cm?2as shown in Fig.1.Mg alloys worked as anodes with an exposed area of approximately 1.77 cm2(1.5 cm diameter)in a 3.5 wt.%NaCl solution electrolyte.The cathode was a commercial air cathode with 0.3 mg cm?240%platinum nanoparticles as catalysts(Fuel Cell Store.Ltd).After battery discharge testing,the reaction products on the anode surface were removed by a 200 g L?1chromic acid solution.Then,the surface morphologies were characterized by SEM(SE,10 kV).The specific capacity,anodic efficiency and specific energy density were calculated as follows [6,24,29]:

Fig.1.Schematic of the Mg-air cell for battery discharge test.

Where?W(g)is the weight loss of samples during the discharge,Wtheo(g)is the theoretical weight loss,Uis the cell voltage,Iis the discharge current(A),tis discharge time(h),Fis the Faraday constant(26.8 Ah mol?1),andxi,ni,mirepresent the mass fraction,number of exchanged electrons and atomic weight,respectively.

3.Results and discussion

3.1.Microstructure analysis

The optical micrographs of as-cast ZT21,ZT23,ZT41 and ZT43 alloys are shown in Fig.2(a–d).These as-cast ZT alloys exhibit a typical dendritic microstructure,which consists ofα-Mg dendrites and intermetallic phases at dendritic regions and grain boundaries.Backscattered electron SEM(BSE-SEM)technique has been used for characterising these four alloys,and results are shown in Fig.2(e–h).With increasing Zn and Sn content,the contrast betweenα-Mg dendrites and interdendritic regions increases,indicating that solute elements Zn and Sn segregate around the interdendritic regions and grain boundaries.The intermetallic phases in ZT alloys are distributed along dendritic regions and grain boundaries.

Wei et al.[30]have conducted a comprehensive microstructure analysis on Mg-Zn-Sn alloys and reported that the microstructure of as-cast ZT43 alloy consists ofα-Mg,Mg2Sn phase,eutectic phase(α-Mg+Mg4Zn7),and globularshaped phase((α-Mg +Mg4Zn7)+Mg2Sn).High magnification SEM images and EDS results in Fig.2(i,j)confirm that ZT alloys have three types of intermetallic phases:eutectic phaseA(α-Mg+Mg4Zn7),Mg2Sn phaseBand a combination of two structures particleC.Instead of globular-shape,particleCdisplays as dolphin-like morphology with a hybrid structure of eutectic phase(α-Mg+Mg4Zn7)on the right side and Mg2Sn on the left side.Thus,the microstructure of as-cast ZT alloys consists ofα-Mg,Mg2Sn phase,eutectic phase(α-Mg+Mg4Zn7),and the hybrid phase((α-Mg+Mg4Zn7)+Mg2Sn).

3.2.Electrochemical analysis

In Mg-air aqueous batteries,a high self-corrosion rate is a critical problem for Mg metal anodes.The corrosion performances of ZT alloys were evaluated by potentiodynamic polarization and electrochemical impedance spectroscopy(EIS).For comparison purpose,the potentiodynamic polarization and EIS results of Mg-2Zn and Mg-4Zn were also presented.Fig.3 shows the Tafel curves of Mg-Zn and ZT alloys,and the fitting value of Ecorrand jcorrare presented in Table 2.The corrosion potential Ecorrof Mg-Zn alloys shows a decrease trend by adding Sn alloy element,indicating that Sn can reduce the cathodic kinetics.The cathode polarization curve indicates the hydrogen evolution rate,showing that ZT alloys have a low hydrogen evolution rate than Mg-Zn alloys.In other words,addition of Sn improves corrosion resistance of Mg-Zn alloys.

Fig.2.The microstructure of as-cast Mg-Zn-Sn(ZT)alloys:(a)and(e)ZT21,(b)and(f)ZT23,(c)and(g)ZT41 and(d)and(h)ZT43.(i)and(j)shows SEM-EDS elemental analysis of the intermetallic phases in ZT43 alloy.

Fig.3.Potentiodynamic polarization curves of Mg-Zn and ZT alloys in 3.5 wt.% NaCl electrolyte solution.

However,corrosion current density jcorrof Mg-Zn alloys varied with the amount of Sn content.It is known that Mg has a“negative difference effect”in which the rate of partial cathodic reaction increases with applied anodic potential.Addition of Sn and Zn has affected the anodic reaction of Mg alloys in NaCl solution.At the anode part of the Tafel curve,ZT21 and ZT41 alloys show apparent passive points(insert image in Fig.3),indicating that MgxZn-1Sn alloys can produce passive films on the surface of the metal anode to decrease their corrosion progress.This agrees with the report by Jiang et al.[31]that solute Sn in Mg matrix may form SnO2at alloy surface in an aqueous solution,improving alloy’s corrosion resistance.However,when the concentration of Sn is higher than 2 wt.%,a large amount of Mg2Sn increases corrosion by galvanic effect [31].

Fig.4.Electrochemical impedance spectra of Mg-Zn and ZT alloys in 3.5 wt.% NaCl solution:(a)Nyquist plots,(b)Bode plots of impedance modulus vs.frequency,and(c)Bode plots of phase angle vs.frequency.

The EIS of Mg-Zn and ZT alloys were investigated at open circuit potential(OCP).The Nyquist plot and their fitting circuit are shown in Fig.4(a),in whichRsis the solution resistance,Rtis the charge transfer resistance that equals to the diameters of capacitive semi-circles,and CPE represents the electric double-layer capacitance at the interface between the anode surface and the NaCl solution.In general,relatively high CPE values indicate large active corrosion areas and lowRtvalues,representing poor corrosion performance [4,32].

It can be seen from Table 3 that ZT alloys have larger CPE values than Mg-Zn alloys,implying that a larger active area for the discharge reaction.ZT alloys also show much lowerRtthan that of Mg-Zn binary alloys,which means that ZT alloys have lower charge transfer resistance than those of Mg-Zn alloys.Thus,the protection of initial oxide films of ZT alloys is not better than Mg-Zn alloys.Fig.4(b)shows the Bode plot in terms of the |Z| vs frequency.For all materials tested,the impedance increases with decreasing frequency and then maintains a stable value at low frequency.The impedance modulus is generally employed to evaluate the corrosion resistance [33,34].Materials with higher impedance modulus usually mean higher corrosion resistance.In this work,the impedance modulus was found to decrease in the order of Mg-2 Zn>ZT21>Mg-4 Zn>ZT41>ZT23>ZT43.

Volta potential differences of metal microstructure have been used to predict corrosion behavior [35,36].This method can measure the potential in a small area to investigate localized corrosion and micro-galvanic activities [37].In this research,we have conducted volta potential measurements on secondary phases of ZT43 alloy,Mg-4Zn and Mg-3Sn alloys for comparison,as shown in Fig.5.The results indicate that all the intermetallic compounds(bright particles)have higher Volta potential than the surrounding Mg matrix.The volta potentials of Mg2Sn and Mg4Zn7phases are around 120 and 55 mV,respectively.Therefore,Mg2Sn phase acts as stronger cathodic sites than Mg4Zn7for corrosion,thus increases the electrochemical activity of Mg-Zn alloy,constant with the results from Tafel curves and EIS.

Fig.5.Volta potential map and potential curves of:(a)Mg-4Zn-3Sn(ZT43),(b)Mg-4 Zn and(c)Mg-3Sn alloys by SKPFM.

3.3.Discharge behavior

The battery discharge performances of ZT alloy anodes were studied in Mg-air electrochemical cells,and the results are shown in Fig.6.Discharge properties of Mg-2Zn and Mg-4Zn alloys are presented for comparison [24].ZT alloys display improved discharge behaviors.Fig.6 shows that the operating voltages of ZT alloy anodes are higher than Mg-Zn alloy.Fig.6(e)shows an interesting phenomenon,jagged fluctuations,observed in the discharge curves,which represents the release of hydrogen as a side reaction.Compared to Mg-Zn alloys,ZT alloys have minor fluctuations,especially ZT21 alloy shows a very smooth discharge reaction at low current densities.

Fig.6.Discharge curves of Mg-Zn and ZT alloy anodes at four current densities:(a)1 mA cm?2;(b)and(e)2 mA cm?2;(c)5 mA cm?2,and(d)10 mA cm?2 in 3.5 wt.% NaCl solution.The cathode is a Pt catalyst carbon cloth [24].

At 1 and 2 mA cm?2,all Mg alloy anodes can be discharged up to 20 h.When the current density increases to 5 and 10 mA cm?2,the discharge voltage fluctuate severely.However,ZT21 and ZT23 anodes can still discharge smoothly for 20 h.The deteriorated battery performance of anode is directly related to the broken balance between the reaction products deposition and detachment [4].Intermetallic compounds Mg2Sn or solid soluble Sn may cause easy detachment or dissolution of discharge products.Compared to ZT21 and ZT23 alloys,the turbulent voltages of ZT41 and ZT43 are caused by the larger number of intermetallic compounds than that in ZT21 and ZT23.In addition,Zn and Sn segregated on interdendritic regions and grain boundaries in Mg alloys could accelerate the self-corrosion,and affect the accumulation of reaction products on the anode surface [24].

The battery properties of ZT alloys anodes in Mg-air batteries are shown in Fig.7(a,b).The battery discharge performance of recast Mg,Mg-Zn alloy are included for comparison [24].ZT alloys have much higher discharge voltage,anodic efficiency,specific capacity,and specific energy density than Mg-Zn and pure Mg at the current density of 1 to 5 mA cm?2.ZT23 alloy shows the highest utilization efficiency of 47.2%,specific capacity of 1009 mAh g?1and specific energy density of 1367 mWh g?1at 2 mA cm?2.ZT21 and Mg-2Zn alloy anodes show the similar anodic efficiency and specific capacity,which are higher than other anode materials at 10 mA cm?2.Although the properties of ZT41 and ZT43 alloys are not better than those of ZT21 and ZT23 alloys,they still much better than Mg-4Zn at the current densities of 1,2,5 mA cm?2.With increasing current density from 1 to 10 mA cm?2,the anodic efficiencies of pure Mg changes from 20.84 to 50.60%,while the anodic efficiencies of ZT21 alloy anode increases from 41.28 to 53.41%.Compared to recast Mg and Mg-Zn alloys,ZT alloy anodes demonstrate a much stable discharge performance.This may be related to the formation,accumulation and detachment of discharge product on the anode surface.

Fig.7.Discharge properties of Mg-Zn and ZT alloy anodes at different current densities for 20 h:(a)average discharge voltage and specific capacity vs current density,(b)anodic efficiency and specific energy density vs current density [24].

Fig.8.Battery discharge performance of ZT anodes comparing with pure Mg and commercial alloys anodes at the current density 5 mA cm?2[16,24,25,38–40].

Fig.9.Intermittent discharge of Mg-air battery with recast Mg,Mg-Zn and ZT anodes at 2 mA cm?2 in 3.5 wt.% NaCl solution.

Gu et al.have studied the discharge performance of Mg-Sn-Zn-Ag alloy anode with three different Sn concentrations(2,4 and 8 wt.%),and reported that Mg-4Sn-1Zn-1Ag(TZQ411)has high anodic efficiency of 69.53% at 120 mA cm?2and 45.30% at 10 mA cm?2[11].Compared to TZQ411,ZT21 and ZT23 show better anodic efficiency at 53.41% and 50.90% at 10 mA cm?2,respectively.Fig.8 compares the battery performance of ZT21 and ZT23 anodes with commercial Mg anodes such as AZ31 and AM60 at 5 mA cm?2[16,24,25,38–40].The specific energy density and anodic efficiency of ZT alloy are much better than pure Mg,AZ31 and AM60 alloys.

Since the discharge in practical application runs periodically,the intermittent discharge is a useful test to simulate the practical application.ZT21 and ZT23 alloys have good discharge performance during 20 h time,so they are tested for intermittent discharging.Intermittent discharge testing was also performed on recast Mg and Mg-2Zn alloy anodes for comparison purpose,which was carried out for 10 h at a constant current density 2 mA cm?2and stop for 5 h as one cycle.The test is terminated at the fifth cycle as these anodes showed a significant voltage drop.Fig.9 shows that ZT,Mg-2Zn and recast Mg anodes exhibited very similar behavior in the first 4 cycles with discharge voltages of ～1.4 V.When the discharge stopped,the OCV increased to 1.8 V.However,at the following cycles,the voltage drop of pure Mg and Mg-Zn were more apparent than ZT alloys.The percentage of voltage drop are calculated after 30 min discharge for each cycle,and results are listed in Table 4,indicating that ZT21 and ZT23 have more stable discharge voltage than Mg and Mg-Zn.

Fig.10.Surface morphologies of anodes after discharge:(a)ZT21,(b)ZT23,(c)ZT41,(d)and(e)ZT43 at 2 mA cm?2 for 20 h in 3.5 wt.% NaCl solution.

Table 4 The voltage drop of 5 cycles(%),compared with the previous cycle.

At the fifth discharge cycle,the voltages of ZT21 and ZT23 are around 1.3 V,while the Mg and Mg-2Zn’s voltages were around 1.2 V.The voltage drops percentages of Mg,Mg-2Zn,ZT21 and ZT23 are 12.39%,10.64%,5.67%,and 4.28%,respectively.As the discharge products accumulated on the anode and electrolyte interface,which cause a decrease of operating voltage.This result indicates that Sn could minimize the discharge product accumulation,improving the discharge stability.

3.4.Anode surface analysis after battery discharging

Fig.10 shows the surface morphology of ZT anode after discharge for 20 h at 2 mA cm?2.ZT anodes have relatively fine corrosion holes uniformly distributed inside grains(pointed by arrowA),and obvious microcracks extend along the grain boundaries(pointed by arrowB).The surface of ZT21 and ZT23 looks relatively flat,while ZT41 and ZT43 have obviously deep and large corrosion pits.Since the standard potential of Sn(?0.1375 V)and Zn(?0.7618 V)are much positive than that of Mg(?2.372 V)[27],the Zn and Sn-rich site near the dendritic and grain boundaries may act as a cathode to accelerate micro-galvanic corrosion of the surrounding matrix.The cracks at the grain boundaries shown in Fig 10(e)may come from the Zn and Sn-rich regions.According to the SKPFM result in Section 3.2,the intermetallics are more positive than the Mg matrix.The corrosion pits in areaAcame from the intermetallic compounds and distributed evenly inside the grains,which accelerate the corrosion of Mg matrix and the discharge performance.However,the deep pits make the corrosion products not easy to fall off.

Fig.11(a)and(b)show the cross-section of Mg-2Zn and ZT21 anodes after discharge for 20 h at 2 mA cm?2.EDS line scan in Fig.11(a)shows that the discharge reaction interface of Mg-2 Zn and ZT21 alloy are enriched with Zn and Zn/Sn elements,respectively,implying that the anode surfaces are protected by ZnO and ZnO/SnOx(x=1 or 2),consistent with the corrosion in NaCl solution [41–43].The EIS test at OCV after discharge at 2 mA cm?2for 30 min is shown in Fig.11(c)and(d).The charge transfer resistances(Rt)is related to the corrosion products film.TheRtof Mg-2Zn-ySn are lower than those of Mg-4Zn-ySn(y=0,1 and 3),so the alloys with high Zn content exist thick discharge product layers.This result is consistent with the previous result that Mg alloys with high Zn content have a large amount of intermetallic compounds to facilitate self-corrosion and discharge products accumulation[24].

Fig.11.Cross-section of(a)ZT21 and(b)Mg-2Zn alloy anodes discharge at 2 mA cm?2 for 20 h with discharge products and their EDS line scan,(c)EIS and the fitting circuit of the investigated alloy anodes after discharge at 2 mA cm?2 for 30 min,and(d)Charge transfer resistance of the EIS results.

TheRtof Mg alloys follows the order of Mg-xZn

4.Conclusion

The effect of Sn on the microstructure,electrochemical and discharge properties of Mg-Zn alloy anodes for Mg-air battery have been investigated.The addition of Sn affects the microstructure of Mg-Zn alloys,forming SnO or SnO2film on the anodes during discharge and significantly improves the battery performance of Mg-Zn alloy at low current densities(<10 mA cm?2).Sn can also stimulate the electrochemical activity of Mg matrix during discharge.Mg2Sn phase has stronger cathodic effect than MgxZnyphase.Mg air battery with Mg-Zn-Sn alloy anode has significantly higher operating voltage than Mg-Zn anode.Compared to Mg-2Zn,the anodic efficiency of ZT23 increase from 41.6% to 50.3% at 5 mA cm?2.Among the ZT anodes,ZT23 has the best battery performance at low current density(<5 mA cm?2)and shows the high energy density of 1367 mWh g?1at 2 mA cm?2.

Acknowledgments

This work is partially supported by the Marsden Fund managed by the Royal Society of New Zealand Te Apˉarangi(Fast-Start Marsden Grant project No.UOA1817).Fanglei Tong is thankful for the scholarship from China Scholarship Council(No.201808060410).