Abstract

The hot deformation behaviors of homogenized direct-chill (DC) casting 6061 aluminum alloys and Mn/Cr-containing aluminum alloys denoted as WQ1 were studied systematically by uniaxial compression tests at various deformation temperatures and strain rates. Hot deformation behavior of WQ1 alloy was remarkably changed compared to that of 6061 alloy with the presence of α-Al(MnCr)Si dispersoids. The hyperbolic-sine constitutive equation was employed to determine the materials constants and activation energies of both studied alloys. The evolution of the activation energies of two alloys was investigated on a revised Sellars’ constitutive equation. The processing maps and activation energy maps of both alloys were also constructed to reveal deformation stable domains and optimize deformation parameters, respectively. Under the influence of α dispersoids, WQ1 alloy presented a higher activation energy, around 40 kJ/mol greater than 6061 alloy’s at the same deformation conditions. Dynamic recrystallization (DRX) is main dynamic softening mechanism in safe processing domain of 6061 alloy, while dynamic recovery (DRV) was main dynamic softening mechanism in WQ1 alloy due to pinning effect of α-Al(MnCr)Si dispersoids. α dispersoids can not only resist DRX but also increase power required for deformation of WQ1 alloy. The microstructure analysis revealed that the flow instability was attributed to the void formation and intermetallic cracking during hot deformation of both alloys.

1. Introduction

Al-Mg-Si-Cu alloys (6xxx series alloys) have been widely applied in auto industry for their high specific strength, excellent formability, and corrosion resistance [1, 2]. Adding of transition elements, such as Mn and/or Cr, will significantly enhance mechanical properties of 6xxx alloys at both ambient and elevated temperatures due to formation of α-Al(FeMnCr)Si dispersoids via homogenization heat treatments [35]. The α-Al(FeMnCr)Si dispersoids can also effectively pin migration of dislocation and substructures and therefore restrain recovery and recrystallization. In these hot processings, the workability of alloys directly relates to their deformation conditions (i.e., temperature and strain as well as strain rate) [6] and their chemical composition [7] and microstructure evolution. The processing maps based on dynamic material modeling (DMM) have been proved to be an effective approach to distinguish optimum processing domains with high power dissipation efficiency from instability domains associated with microstructural defects for hot deformed alloys [8, 9]. Wang et al. [10] predicted the optimum processing parameters for 7050 Al alloy as 653–693 K and 0.001–0.18 s−1 for temperatures and strain rates, respectively. Kai et al. [11] proposed that the optimum processing window for 6X82 was 465–535°C and 0.09–1.2 s−1. Owing to alloying elements addition and dispersion strengthening, the plastic characteristics of 6xxx alloys with α-Al(FeMnCr)Si dispersoids are quite different from those of traditional 6xxx alloys and the related study is rare. Given current status, it is necessary to clearly understand effect of α-Al(MnCr)Si dispersoids on hot deformation behaviors and workability of 6xxx alloys.

Arrhenius type constitutive equations established by Jonas et al. [12] are also widely used to characterize flow behaviors and interpret activation energy (Q) of alloys deformed over a wide range of deformation temperatures and stain rates. Q can be expressed as necessary free energy for dislocation to slip on slip planes and hence indicate deformation resistance of alloys [13]. Q was used to be treated as a constant under different hot deformation conditions, which is contrary to recent researches. The results [13] revealed that the activation energy value of AA7150 Al alloy during hot deformation was not constant but decreased with increasing of deformation temperature and strain rate. Wang et al. [10] proved that Q of 7050 Al alloy during hot deformation were greatly affected by activation of dynamic recovery mechanism. Qin et al. [14] demonstrated that the activation energy map of the composites containing Sc and Zr was divided into two regions associated with dynamic precipitation. Qian et al. [3] showed that activation energy increased with increasing Mn contents in 6082 Al alloy. However, the influences of α-Al(MnCr)Si dispersoids on activation energy evolution of corresponding alloys during hot deformation process remain rarely reported.

The hot deformation behaviors of alloys are directly affected by the dynamic balance between work hardening and softening mechanisms [1517]. Therefore, the microstructural evolution and flow characteristics are the key to understand plastic behaviors of Al alloys during deformation process. It is commonly accepted that DRV and DRX are the main softening mechanisms for alloys deformed at high temperature. But how α-Al(MnCr)Si dispersoids affect microstructural evolution and softening mechanics of hot deformed 6xxx is still not well clarified.

In the present study, a new type of 6xxx alloy with addition of Mn and Cr (donated as WQ1) was designed aiming to form fine α-Al(MnCr)Si dispersoids. The hot deformation behaviors of 6061 and WQ1 Al alloys were compared as a function of deformation temperatures and strain rates. The materials constants and activation energies of the two studied alloys were derived by application of hyperbolic-sine constitutive equation. The activation energies of the two alloys as a function of deformation conditions were constructed and compared to illustrate the impacts of deformation parameters and α-Al(MnCr)Si dispersoids in detail. The processing maps that present processing safe domains and flow instability domains of both alloys were also established. Effects of α dispersoids on dynamic softening mechanisms and working windows of studied alloys were also analyzed. The workability of both alloys was discussed in combination with processing maps and activation energy maps. Microstructure evolution of both alloys was also observed to reveal the mechanism by which dispersoids affect softening mechanisms and verify workability of the alloys during hot deformation.

2. Experimental Procedure

The samples used in the present work were from ingots with 150 mm in diameter of 6061 and Mn/Cr-containing aluminum alloys denoted as WQ1 which were fabricated by DC casting. The ingots of 6061 and WQ1 were homogenized at 555°C for 6 h with a heating rate of 111°C/h and 530°C for 6 h with a heating rate of 133°C/h, respectively, followed by air cooling to ambient temperature. Chemical compositions of these two alloys were examined by SPECTROLAB optical emission spectroscopy and corresponding results are listed in Table 1. The cylindrical hot compression samples with 10 mm in diameter and 15 mm in height were extracted from at half radius to center position of homogenized ingots as shown in Figure 1. To characterize the distribution of α-Al(MnCr)Si dispersoids in studied alloys on as-homo condition, the samples were etched in 0.5 wt.% HF solution for 15 s at ambient temperature. To determine average grain size of two alloys, the homogenized samples were anode-coated at 18 V in a solution of 38 wt.% H2SO4 + 43 wt.% HNO3 + 19 wt.% deionized water for 1.5 min.

Isothermal compression deformation tests were performed on a Gleeble-3500 thermomechanical simulator with deformation temperature ranging from 300 to 500°C and strain rate range of 0.01–10 s−1 to strain of 1.2. To assure a uniform deformation process, graphite foils were placed at both sides of the sample and the anvil. Thermocouples were welded on cylindrical surface of each sample to record and control deformation temperature. The samples are heated to the targeted temperature with a heating rate of 5°C/s and soaked for 2 min to ensure thermal uniformity.

All deformed samples were immediately quenched into water at room temperature after hot deformation to preserve compressed microstructures. The specimens for microstructural observation were selected from the central part of the deformed samples along compression direction and metallographically prepared. Prior to electron-back scattering diffraction (EBSD) analysis, the prepared samples were electropolished in a solution (10 wt.% perchloric acid and 90 wt.% ethanol) at the voltage of 20 V for 40 s to produce strain-free surface. All EBSD data was processed by commercial software HKL Channel 5. In EBSD inverse pole figures (IPF), high-angle boundaries (HABs) expressed as black lines were defined as a misorientation larger than 15°, while low-angle boundaries (LABs) shown as white lines are set as a misorientation between 1° and 15°. Boundaries with misorientation less than 1° were taken as noises. Microstructure evolution was studied by JEOL JEM-2100F field emission transmission electron microscope (TEM) operated at 200 kV. The TEM samples were prepared by twin-jet electropolishing method at −20°C in a solution containing 75 vol.% of methanol and 25 vol.% HNO3. Convergent beam electron diffraction (CBED) patterns were used to measure the thickness of the observed TEM foils. The volume fraction of dispersoids of 6061 and WQ1 alloys were determined, using equation proposed by Li and Arnberg [18].

3. Results and Discussion

3.1. Microstructure after Homogenization Annealing

Figure 2 shows the optical microscope images of 6061 and WQ1 alloys on homogenization condition after anode-coating. The microstructures of both alloys are composed of equiaxed grains with average of 95.75 ± 6.34 μm and 216.76 ± 13 μm for 6061 and WQ1, respectively. According to previous study [19], grain size variation plays an insignificant role in affecting flow stress during high temperature deformation. However, grain size in the range among commercial alloys has limit effect on flow stress.

It can be seen from dark field optical microscope images in Figure 3 that few and unevenly size of dispersoids precipitate in 6061 (Figure 3(a)) due to lack of Mn and Cr elements. Meanwhile fine and dense dispersoids generated in WQ1 (Figure 3(b)) during homogenization heat treatment. The distribution of dispersoids in 6061 and WQ1 alloys is inhomogeneous, especially in former, and large area of dispersoid free zone (DFZ) is observed. The nonuniform distribution of dispersoids is caused by initial segregation of Mn/Cr at interdendritic regions during solidification [2] and low diffusion coefficients of Mn/Cr elements.

TEM images of as-homo 6061 and WQ1 are shown in Figure 4. It is obvious that the number density of dispersoids in WQ1 is significantly higher than that of 6061. Based on quantitative analysis of TEM image of dispersoids in both alloys, the average equivalent diameters are determined as 124 ± 6 nm and 129 ± 23 nm for 6061 and WQ1, respectively. According to the method proposed by Li and Arnberg [18], the volume fraction of dispersoids in 6061 and WQ1 was calculated as 0.1% and 1.1%, respectively. The chemical composition of the dispersoids was also examined by energy disperse spectroscopy (EDS) as shown in Figure 4(c). Based on the morphology [5] and composition, the dispersoids were identified as α-Al(MnCr)Si dispersoids [20].

3.2. Flow Stress Behavior

The true strain-true stress curves of the isothermal compression deformation of studied alloys are drawn in Figure 5. It is obvious that both deformation temperature and strain rate are of great importance to the flow behaviors of the alloys. Flow stresses of both alloys increase with decreasing deformation temperatures as well as increasing strain rates. Overall, the true stress ascends rapidly at initial stage of deformation and subsequently reaches a plateau or slightly increases throughout compression process, indicating a dynamic equilibrium between work hardening and dynamic softening. All results show the flow stresses of WQ1 (Figure 5(b)) are obviously higher than those of 6061 (Figure 5(a)) under the same deformation conditions.

At initial stage, compression deformation triggers generation of massive microscopic defects (vacancies and dislocations). The complex interaction within these overwhelming microdefects causes dislocation multiplication and tangling as well as vacancy-dislocation kinking and jogging throughout the material, leading to conspicuous work hardening effect. With increasing of strain, two mutually counteracting phenomena impact the microstructure simultaneously. On one side, newly formed movable defects trapped into the preexisting dislocations and vacancies causing further strengthen work hardening. On the other side, a critical dislocation density value was reached; the dislocation cross-slip was activated along favorable slip system, resulting in dislocation annihilation and recombination as well as reduction in dislocation density. Therefore, the dynamic balance between work hardening and dynamic softening eventually developed into a steady state of flow stress during further compression. As the steady state of flow stress always presents at the strain where it reaches its peak stress, hence, peak stress is applied for the following constitutive analysis.

Figure 6 illustrates the peak flow stresses evolution of 6061 and WQ1 as function of deformation temperatures and strain rates. There is a strong connection between peak flow stresses and deformation conditions (deformation temperature and strain rates). The peak flow stresses are set as the maximum stress value of each flow curve, which increase associated with increasing strain rates and decreasing deformation temperatures. Under the same deformation condition, the peak stresses of 6061 (Figure 6(a)) are lower than those of WQ1 (Figure 6(b)), which means α-Al(MnCr)Si dispersoids considerably increase deformation resistance of WQ1. The mobility of dislocations and grain boundaries is thermally activated. Restrained by deformation temperature, lower temperature will limit DRV and DRX, therefore causing a higher flow stress. Lower strain rates can provide sufficient time for activation energy of dynamic softening mechanism leading to lower flow stress.

3.3. Constitutive Analysis

Constitutive equations are often applied to reveal the intrinsic relationships between flow stress of alloys and deformation conditions (deformation temperatures and strain rates) and predict flow stress behaviors of alloys during hot deformation. Among numerous equations, Arrhenius constitutive equation proposed by Sellars and Tegart is mostly utilized set [21]:

is strain rate (s−1), A1, A2, A, n1, β, n, and α (α = β/n1) are constants of the material, σ is the flow stress (MPa), Q is deformation activation energy (kJ/mol), R is the gas constant (8.314 J·mol−1·K−1), T is the absolute temperature (K), and Z is Zener–Hollomon parameter.

Based on previous study, a power relationship equation (2) is applied to low flow stress situation. On the contrary, an exponential relationship equation (3) is applied to high flow stress situation. The hyperbolic sine law equation (4) is suitable for wide range of flow stress situation. Applying natural logarithms to both sides of equations (2)–(4), the equations for constants of the materials can be obtained:

Subsequently, the values of n1 and β can be determined by the slope of the plots for and , respectively. According to α = β/n1, the value of α can be obtained by substituting the average values of n1 and β.

After undergoing partial differentiation on equation (7), Q for hot deformation of the studied alloys can be calculated as

By substituting experimental data of 6061 and WQ1 into equation (8), the corresponding material constants (A, β, α, and n) and Q can be obtained and calculated values are presented in Table 2. The average Q value of 6061 alloy is 148.73 kJ/mol; meanwhile Q value of WQ1 reaches 189.12 kJ/mol. Apparently, Q of WQ1 is higher than that of 6061, which is caused by synthetic effect of higher solute level of Mn and α-Al(MnCr)Si dispersoids in WQ1. Solute atoms can work as obstacles to restrict dislocation motions during hot deformation, causing ascending activation energy of alloys [19, 22]. α-Al(MnCr)Si dispersoids are able to effectively pin dislocation migration as well as subgrain boundaries migration. Therefore, extra energy was required for dislocations or subgrain boundaries for unpinning from α-Al(MnCr)Si dispersoids, therefore resulting in a greater Q for WQ1 and higher hot deformation resistance.

Recent works [10, 14, 23] demonstrate that activation energy for hot deformation of a given material varies with deformation parameters rather than remaining a constant. Shi et al. [13] proposed a revised Sellars’ constitutive equation (equation (9)) to present the relation between activation energy and deformation parameters (T and ) (equation (12)): represents the slope of versus plotted at different deformation temperatures and is the slope of versus 1/T plotted at different strain rates. The value of Q at various deformation conditions can be obtained according to M and N; hence, the corresponding Q maps for 6061 and WQ1 can be illustrated (Figure 7).

For 6061 alloys (Figure 7(a)), Q is 170.82 kJ/mol at 300°C with strain rate of 0.01/s, while it degrades to 132.33 kJ/mol at 550°C with strain rate of 10/s. For WQ1 alloys (Figure 7(b)), Q is about 208.12 kJ/mol at 300°C with strain rate of 0.01/s, while it decreases to 152.03 kJ/mol at 550°C with strain rate of 0.01/s. Q of WQ1 is significantly increased compared to that of 6061 under the same deformation conditions. It is reported that Q is directly connected with density and movement of dislocations [10]. The presence of α-Al(MnCr)Si dispersoids in WQ1 alloy effectively retards movement of dislocations and (sub)grain boundaries, resulting in higher Q value than 6061’s.

As shown in Figure 7, it is obvious that Q is not constant but is influenced by thermomechanical conditions, that is, deformation temperatures as well as strain rates. Q of 6061 and WQ1 alloy decreases with increasing deformation temperatures and increasing strain rates. On one hand, dynamic restoration process is expedited when the alloy is deformed at higher strain rates [24], leading to promoted dislocation movement. On the other hand, increasing strain rates leads to increasing resolved shear stress along the dislocation slipping direction [25], resulting in advanced initiation of dislocation movement. With increasing deformation temperature, DRV and DRX are greatly boosted, including formation, growth, and polygonization of well-defined subgrains as well as formation of recrystallized grains, respectively. Therefore, the activation energy decreases with increasing deformation temperature.

The activation energy can reflect workability of alloys under various deformation conditions. Higher Q domains in activation map are often related to inferior workability [26], while lower Q domains indicate less processing difficulty. There is a correlation between Q and hot deformation mechanism as well as microstructural evolution, which will be discussed in the following section.

3.4. Processing Maps

Based on the dynamic material model, the processing map is established by superposition of a flow instability map over a power dissipation efficiency map. The efficiency of power dissipation (η) is employed to assess the power dissipation capacity of the material, which is represented as [27, 28]where m is strain rate sensitivity.

Basically, higher efficiency of power dissipation indicates that more power is consumed by microstructural evolution, meaning better hot workability. The flow instability of the material, including microcracks, flow localization, and adiabatic shear bands, can be expressed by a parameter (ξ) [8, 29, 30]:

The instability map is composed of two parts according to the value of ξ. One is the processing zone, which is considered as safe domain when ξ > 0, while the remaining part (ξ < 0) is taken as flow instability domain.

Figure 8 shows the processing maps at the true strain of 1.2 for 6061 and WQ1 alloys as a function of deformation temperatures and strain rates. As shown in Figure 8, the contour numbers represent the value of power dissipation efficiency and the shadow area stands for instability regions. A safe processing domain is always with a high-power dissipation efficiency and absence from instability region. The selected safe domains are marked in red boxes in Figure 8 and the corresponding processing window values of deformation temperatures and strain rates are listed in Table 3.

As shown in Figure 8(a) and Table 3, there are three safe domains for 6061 alloys; domain I locates at temperatures ranging from 485 to 550°C with strain rates of 0.01–0.08 s−1 where power dissipation efficiency is 0.34–0.40 and domain II locates at temperatures of 440–480°C and strain rates of 0.05–0.6 s−1 with power dissipation efficiency of 0.34, while optimal domain III covers the region with higher temperatures and higher strain rates (530–550°C, 6–10 s−1). For WQ1 alloys (Figure 8(b)), there is mainly one safe domain, which locates at temperature ranging from 485 to 550°C with strain rates of 0.01–0.1 s−1 whose power dissipation efficiency is 0.26–0.29. Comparing processing maps of 6061 and WQ1 alloys, the safe domain of both alloys locates at high deformation temperatures with low strain rates but 6061 alloy is accompanied with higher power dissipation efficiency. It has been reported [31, 32] that the deformation mechanism is DRV if power dissipation efficiency value during hot deformation is less than 0.3; otherwise, higher value is associated with DRX, which indicates that the dynamic softening mechanism of WQ1 is inhibited by α-Al(MnCr)Si dispersoids. This point will be verified in the examination of microstructure evolution part. Overall, safe processing domain and low activation energy should be key determinants for the optimum workability and elimination of deformation defects.

3.5. Microstructure Evolution
3.5.1. Interpretation of Safe Domains

Figures 9(a) and 9(b) show the microstructures of deformed 6061 and WQ1 alloys at safe domains with low Q (550°C and 10 s−1). The microstructures of both alloys mainly consisted of elongated grains and a certain area fraction of recrystallized grains, indicating a typical DRV process. Moreover, different degrees of recrystallization occur for 6061 and WQ1 alloys. It can be seen from Figure 9(a) that recrystallization happens extensively along deformed grain boundaries (marked by black box) and the area fraction of recrystallized grains is as much as 50.12% for 6061 alloy. Figure 9(b) shows that partial recrystallization mainly initiates around triple junction grain boundaries (noted by black box) and the area fraction of recrystallized grains is only 9.78% for WQ1 alloy. Particle-stimulated nucleation (PSN) of recrystallization has been identified as an active nucleation mechanism of DRX for hot deformed alloys, which is strongly dependent on the size of second phase particles. Normally, particles can induce PSN effect and enhance the DRX when they are larger than 1 μm, while they can, on the other hand, suppress the DRX for the strong pining effect when particles are smaller than 1 μm [33, 34]. Hence, recrystallization in WQ1 alloy is severely restrained by the presence of large amount of nanoscale (∼120 nm) α-Al(MnCr)Si dispersoids. Misorientation distributions based on Figures 9(a) and 9(b) are presented in Figures 9(c) and 9(d), respectively. It is obvious that WQ1 alloys hold more percentage of LABs and lower average misorientation angle (θave) compared to 6061 alloys, which is consistent with results of Figures 9(a) and 9(b). The substructure of 6061 (Figure 9(e)) is composed of abundant equiaxed subgrains with LABs and recrystallized grains with size of 10 μm. It can be concluded that the main mechanism for DRX is continuous dynamic recrystallization [35]. The recrystallized grains formed by transformation of LABs to HABs (pointed out by black box). For WQ1 alloys (Figure 9(f)), the microstructure is composed of subgrains and recrystallized grains with size around 5 μm. It should be noticed that some subgrains in WQ1 alloys are elongated rather than equiaxed (marked by black arrows), indicating that subgrain rotation is strongly suppressed. Owing to pinning effect of α-Al(MnCr)Si dispersoids, subgrain rotation and (sub)grain boundaries migration are inhibited. Hence, DRV was main dynamic softening mechanism for WQ1 alloy at its safe processing domain. Figure 10(a) demonstrated the interaction of α-Al(MnCr)Si dispersoids with dislocation and dislocation walls when WQ1 was deformed at 450°C and 0.1 s−1. It can be seen that dislocation walls movement is retarded and dislocations bow when encountered with dispersoids, showing a strong pinning effect of α-Al(MnCr)Si dispersoids on dislocation migration. Moreover, dispersoids also inhibit grain boundaries (Figure 10(b)), even though WQ1 was compressed at 550°C. Hence, DRV in WQ1 is significantly restrained due to pinning effects of α-Al(MnCr)Si dispersoids on dislocation motions and grain rotations [19, 24].

3.5.2. Interpretation of Instable Domains

The flow instability is often related to generation of microstructural defects, including flow localization, adiabatic shear bands, micro/macro cracks, and intermetallic/particle cracking and debonding. Micro/macro damages are closely related to shear bands and flow localization. Based on research of Ramanathan et al. [36], flow localization refers to discontinuous and inhomogeneous deformation around adjacent areas resulting from composition differences, interphase boundaries, and loading conditions [37].

The optical micrographs of 6061 deformed at 300°C with 0.1 s−1 to a strain of 1.2 were shown in Figure 11. Micro cracks and voids (marked by arrows) are observed around intermetallic particles indicating the formation of micro defects related to intermetallic particles.

In flow instable domain of 6061 alloys, cracking of primary Mg2Si phases (pointed in Figure 12(a)) along compression direction and interfacial debonding between Mg2Si and Al matrix was observed. Those sites are preferential places for micro crack formation and propagation. Flow localization happens around brittle Fe-bearing particles (Figure 12(b)) and cracked Fe bearing intermetallic causes matrix cracks initiation and propagate perpendicular to compression direction.

As the compression proceeds, both primary Mg2Si and Fe-containing particles are prone to cracking and debonding rather than deforming in WQ1 alloys (Figure 13). Thus inhomogeneous deformation leads to severe stress concentration and enables cracking at the interface between the second phase particles and Al matrix under the action of external loading. Visible voids and cracks formed due to flow localization (Figures 13(a) and 13(b)) and then the pores grew and finally formed a fracture. Multiple voids and micro cracks tend to expand by linking adjacent voids.

4. Conclusions

(1)A great amount of α dispersoids formed in WQ1 alloy, which is with addition of Mn and Cr, while limited α dispersoids were observed in 6061 alloy. Flow stresses of both alloys increased with decreasing temperatures and increasing strain rates. Furthermore, flow stresses of WQ1 alloy were higher than those of 6061 alloy under all deformation conditions owing to dispersion strengthening caused by α-Al(MnCr)Si dispersoids.(2)The materials constants in constitutive equation and activation energies of WQ1 and 6061 alloys during hot deformation were calculated by applying constitutive analysis. Due to presence of α-Al(MnCr)Si dispersoids in WQ1, activation energy of it for hot deformation was 40.4 kJ/mol higher than that of 6061.(3)The activation energy maps and processing maps of the two alloys during hot deformation were constructed. The activation energy of both alloys is variant and is greatly influenced by deformation temperatures and strain rates. Optimized hot deformation parameters for WQ1 and 6061 were determined with high dissipation efficiency and low activation energy within safe processing domains.(4)The microstructure analysis revealed that main dynamic softening mechanism for 6061 under optimized hot deformation parameters was DRV and extensive DRX, while DRV was predominant softening mechanism for WQ1 due to strong pinning effect of α-Al(MnCr)Si dispersoids on subgrain rotation and subgrain boundaries migration.(5)For both alloys, flow instability was associated with cracking and debonding of intermetallic particles as well as void formation during hot deformation.

Data Availability

All experimental data are available.

Disclosure

Part of the manuscript was presented in Sub-Session: Materials Forming and Modeling, but the manuscript was not submitted in the session. Furthermore, the content of the manuscript was revised.

Conflicts of Interest

The authors declare no conflicts of interest.

Acknowledgments

The authors acknowledge the financial support provided by Jincheng Science and Technology Plan Projects of Shanxi Province (no. 201702014) and National Natural Science Foundation of China (Grant no. U1864209).