Abstract

A series of Zn-substituted hydroxylapatites [(ZnxCa1−x)5(PO4)3OH, Zn-Ca-HA] with the Zn/(Zn + Ca) molar ratio () of 0~0.16 was prepared and characterized, and then the dissolution of the synthesized solids in aqueous solution was investigated by batch experiment. The results indicated that the aqueous zinc, calcium, and phosphate concentrations greatly depended on the Zn/(Zn + Ca) molar ratio of the Zn-Ca-HA solids (). For the Zn-Ca-HA dissolution at 25°C with an initial pH of 2.00, the final solution pH increased, while the final solution calcium and phosphate concentrations decreased with the increasing . The final solution zinc concentrations increased with the increasing when and decreased with the increasing when = 0.08~0.16. The mean values for (ZnxCa1−x)5(PO4)3OH at 25°C decreased from 10−57.75 to 10−58.59 with the increasing from 0.00 to 0.08 and then increased from 10–58.59 to 10–56.63 with the increasing from 0.08 to 0.16. This tendency was consistent with the dependency of the lattice parameter on . The corresponding free energies of formation increased lineally from −6310.45 kJ/mol to −5979.39 kJ/mol with the increasing from 0.00 to 0.16.

1. Introduction

Phosphate apatites form an enormous mineral group due to their huge isomorphic capacity [1, 2], which play an important role in many research areas such as biomaterials and environmental science [37].

As the main inorganic constituent of bone and dental enamel of vertebrates, calcium hydroxylapatite (HA) has been broadly used in osteoinductive coatings, bone replacement and repair, dental orthopaedics, and so forth [3, 6, 811]. The substitution of trace ions in hydroxylapatite can affect not only its lattice parameters and , crystallinity, and morphology, but also its dissolution mechanism and other physicochemical properties [1012]. Zinc is one of the most important essential trace elements for the growth of humans and its incorporation in Ca-hydroxylapatite can significantly improve the bioactivity of Ca-hydroxylapatite [810]. The slow release of zinc substituted in an implant Ca-hydroxylapatite material can also promote bone metabolism and growth around the implant. Thus, the zinc-substituted hydroxylapatite can be a novel biomaterial for bone tissue engineering [10, 13].

Calcium hydroxylapatite (HA) can also be applied to immobilize dangerous metallic compounds in metal-contaminated soils and industrial wastewaters due to its huge ion substitution capacity, which can considerably decrease the mobility and bioavailability of Zn2+, Pb2+, Cd2+, Cu2+, Ni2+, and U2+ by transforming these toxic metal ions into some new forms having low solubility and high geochemical stability [1, 5, 1418]. Heavy metal cations can easily substitute for Ca2+ in the hydroxylapatite structure and form zinc-calcium hydroxylapatite (Zn-Ca-HA), lead-calcium hydroxylapatite (Pb-Ca-HA), or cadmium-calcium hydroxylapatite (Cd-Ca-HA) through dissolution-precipitation, ion-exchange, or adsorption process [15]. Therefore, a fundamental knowledge of the apatite physicochemical properties, especially the solubility, stability, and water-mineral interaction, is required to understand mineral evolution and natural phenomenon or to optimize the industrial processes concerning apatite [5, 19, 20].

However, the thermodynamic data for Zn-substituted hydroxylapatites are now lacking, regardless of the fact that its dissolution and elemental release from solid to aqueous solution exert a great influence on the cycling of zinc, calcium, and phosphate. So far, no experiment on the dissolution and stability of the Zn-substituted hydroxylapatite [(ZnxCa1−x)5(PO4)3OH] has been carried out, for which little information has been reported in literatures. Hence, no thermodynamic data can be obtained to assess the bioactivity and bioavailability of an implant Zn-Ca-hydroxylapatite biomaterial or the environmental risk of zinc concerning the Zn-substituted hydroxylapatite. Additionally, the previous data and results about the effects of the Zn substitution for the Ca sites on the apatite structure and properties are still ambiguous and rather inconsistent [811, 21, 22].

In this work, calcium hydroxylapatite [Ca5(PO4)3OH, Ca-HA] and Zn-substituted hydroxylapatites [(ZnxCa1−x)5(PO4)3OH, Zn-Ca-HA] with various Zn/(Zn + Ca) atomic ratios were prepared and the influences of zinc replacement on the hydroxylapatite properties were investigated with XRD, FT-IR, FE-SEM, and FE-TEM instruments. Then, the dissolution of the synthesized solids and the release of components (Zn2+, Ca2+, and ) were studied, and the solubility product () and the corresponding free energy of formation () of the Zn-Ca hydroxylapatites were determined.

2. Experimental Methods

2.1. Solid Preparation and Characterization
2.1.1. Synthesis

The synthesis of the Ca-HA and Zn-Ca-HA solids was carried out by the precipitation method after the following precipitation reaction: 5M2+ + + OH- = M5(PO4)3OH, where M = Ca for Ca-HA and (Zn + Ca) for Zn-Ca-HA. An aqueous solution with [P] = 0.12 mol/L was first prepared by dissolving NH4H2PO4 into ultrapure water, and a series of the mixed aqueous solutions with [Zn + Ca] = 0.4 mol/L were then prepared by dissolving Zn(CH3COO)2·2H2O and Ca(CH3COO)2·H2O into ultrapure water. The moles of Zn(CH3COO)2·2H2O and Ca(CH3COO)2·H2O were varied in each preparation to get the mixed aqueous solutions with different [Zn]/[Zn + Ca] molar ratios of 0, 0.02, 0.04, 0.06, 0.08, 0.10, 0.12, 0.14, 0.16, 0.18, and 0.20. 250 mL of the Zn2+ and Ca2+ mixture solution was added to 250 mL of 4.4 mol/L CH3COONH4 buffer solution and then 500 mL of the NH4H2PO4 solution was also added with vigorous stirring at °C, which resulted in the forming of white suspension (Table 1). The NH4OH solution was used to adjust the pH of the resulting suspension to 7.5. The suspension was aged at 100°C for 48 h and then subjected to suction filtration. Finally, the white precipitates obtained were cleaned cautiously with ultrapure water and dried at 70°C for 16 h.

2.1.2. Characterization

10 mg of each synthetic solid was digested in 20 mL of 1 mol/L HNO3 solution and then diluted to 100 mL with ultrapure water. It was measured for zinc, calcium, and phosphate using an inductively coupled plasma-optical emission spectrometer (ICP-OES, Perkin-Elmer Optima 7000 DV) to calculate the solid compositions. The solids were measured using a powder X-ray diffractometer (XRD, X’Pert PRO) that was set to 40 kV and 40 mA with a Cu Kα radiation at a scan speed of 0.2°/min. Phase identifications were made by comparing the recorded XRD patterns of the solids with the reference code 00-024-0033 for calcium hydroxylapatite and the reference code 00-020-1427 for ammonium zinc phosphate (NH4ZnPO4) from the ICDD standards. A field emission scanning electron microscope (Hitachi FE-SEM S-4800) was used to observe the morphology of each solid.

2.2. Dissolution Experiments

Each dissolution was designed by adding 2.0 g of the synthesized Ca-HA or Zn-Ca-HA to 100 mL of HNO3 solution (pH = 2) or ultrapure water (pH = 5.6) or NaOH solution (pH = 9) in a series of 150 mL polypropylene bottles, which were then capped and soaked in the water bath of 25°C. From the bottles, 5 mL of the aqueous solutions was sampled at 23 intervals (1, 3, 6, 12, 24, 48, 72, 120, 240, 480, 720, 1080, 1440, 1800, 2160, 2880, 3600, 4320, 5040, 5760, 7200, 7920, and 8640 hours), filtered through a 0.22-μm-pore-size filter and stabilized using 0.2% HNO3 solution. After each sampling, 5 mL of the corresponding HNO3 solution, ultrapure water, or NaOH solution was added. The dissolved zinc, calcium, and phosphate in the aqueous solutions were measured using ICP-OES. After 8640 h dissolution, the solids were sampled, cleaned, dried, and then examined using XRD, FT-IR, FE-SEM, and FE-TEM.

2.3. Thermodynamic Calculations

The PHREEQC program (version 3) [23] was used to calculate the aqueous activities of Zn2+, Ca2+, , and OH-, and then the ion activity product (IAP) for Ca-HA and Zn-Ca-HA was determined after the mass-action expression. The minteq.v4.dat and llnl.dat databases were chosen in the simulation [24]. The following solution species were considered in the calculation: Zn2+, ZnOH+, , Zn, Zn,  , ZnHPO4, and ZnH2 for zinc; Ca2+, CaOH+, Ca, CaHPO4, and CaH2 for calcium; and ,  , H2, H3, Ca, CaHPO4, and CaH2 for phosphate.

3. Results and Discussion

3.1. Solid Characterizations
3.1.1. Chemical Component

The component of the precipitate was dependent on the [Zn]/[Zn + Ca] molar fraction in the starting solution (Table 1). The Zn/(Zn + Ca) molar ratio () and the (Zn + Ca)/P molar ratio of each solid sample were equal to the designed composition of the Zn-Ca hydroxylapatites [(ZnxCa1−x)5(PO4)3OH] when .

3.1.2. XRD

The XRD results proved that Ca-HA and Zn-Ca-HA before dissolution (Figure 1(a)) and after dissolution (Figure 1(b)) were the apatite group minerals that belong to the hexagonal crystal system P63/m. The precipitate of was identified to be Ca-HA (ICDD reference code 00-024-0033). The XRD patterns of the Zn-Ca-HA precipitates of differed from each other only in their peak location, peak intensity, and peak width. The (002), (211), (102), and (210) reflection peaks of the solid samples shift regularly and slightly to the high-angle direction with the increasing due to the replacement of Ca2+ (0.099 nm) by Zn2+ (0.074 nm), which indicated that the Zn-Ca-HA solids were a continuous solid solution when [25]. When , the characteristic diffraction peaks for ammonium zinc phosphate (NH4ZnPO4) were also observed and the peaks for Zn-Ca-HA weakened, which showed that NH4ZnPO4 gradually became the main product; when , the diffraction peaks of HA disappeared in our preliminary experiment. The XRD examination showed that the characters of the Ca-HA and Zn-Ca-HA samples before and after dissolution were not obviously distinguishable (Figure 1). No secondary solid phases formed in the Ca-HA and Zn-Ca-HA dissolution.

The continuous Zn-Ca-HA solid solution could be formed within limited [8, 10, 25]. The solids prepared can be examined for their compositional homogeneity by considering the broadening of the powder XRD peaks of the major reflections [26]. The XRD peak width significantly increased with the increasing from 0.00 to 0.16, which indicated that the crystallinity of Zn-Ca-HA considerably decreased with the increasing . On the other hand, the NH4ZnPO4 phase formed when . The peak intensity of NH4ZnPO4 increased and the peak intensity of apatite decreased with the increasing from 0.16 to 0.20. No parascholzite (CaZn2(PO4)2·2H2O) phase was observed when = 0.16~0.20 [27].

The cell parameter decreased with the increasing from 0.00 to 0.08, increased with the increasing from 0.08 to 0.12, and then decreased with the increasing from 0.12 to 0.16 (Figure 2), which had also been confirmed by some previous researchers [8, 10, 25]. The cell parameter decreased up to and began to increase over [25]. The cell parameter decreased with the increasing up to 0.10 and increased over [8, 10]. The cell parameter decreased with a lower zinc substitution in the Ca-HA lattice ( = 0.00~0.08) because the ion radius of Ca2+ (0.099 nm) is larger than that of Zn2+ (0.074 nm). The increase in the cell parameter for higher (0.08~0.12) was attributed to the increasing amount of lattice H2O that could incorporate in OH sites in the apatite structure [25].

3.1.3. FE-SEM

Figure 3 shows the FE-SEM images of the solids with various to 0.19. The Ca-HA solid () was an aggregate of fine rod-like particles with 20~50 nm in width and 100~150 nm in length. The apatite particle size decreased with the increasing up to 0.19, which also indicated that the crystallinity of the apatite solids decreased with the increasing , as showed also in the XRD diffraction (Figure 1).

3.2. Dissolution
3.2.1. Change of Solution pH and Elemental Concentrations with Time

Evolution trends of the solution pH and concentrations of Zn2+, Ca2+, and with time for the Zn-Ca-HA dissolution in HNO3 solution (pH = 2.00) or ultrapure water (pH = 5.60) or NaOH solution (pH = 9.00) were showed in Figures 4(a), 4(b), and 4(c).

For the Zn-Ca-HA dissolution at 25°C with an initial pH of 2.00 (Figure 4(a)), the solution pH increased from 2.00 to 4.42~4.90 in 1 h dissolution and became stable (pH = 4.88~6.43) after 5040~5760 h dissolution. Commonly, the solution pH increased with the increasing Zn/(Zn + Ca) molar ratios of the Zn-Ca apatites () (Figure 5). The solution Zn2+, Ca2+, and concentrations were greatly affected by (Figure 5). The solids with lower (0.00, 0.02, 0.04, 0.06, and 0.08) showed a different dissolution process from the solids with higher (0.10, 0.12, 0.14, 0.16, 0.17, and 0.19) (Figures 4(a) and 5).

For the Zn-Ca apatites with lower or higher , the solution Ca2+ concentrations increased quickly with time and reached the highest values after 24~480 h dissolution. Then the concentrations decreased slowly and were stable after 1800~2160 h dissolution. The solution Zn2+ concentrations increased rapidly with time and reached the highest values in 1 h dissolution and then decreased progressively and attained a stable state after 4320~5040 h dissolution. For the Zn-Ca apatites with higher or lower , the solution Ca2+ concentrations increased steadily with increasing time and reached the highest values after 720 h dissolution and then decreased and became stable after 5760 h dissolution. The solution Zn2+ concentrations increased quickly with increasing time and achieved the highest values after 1~24 h and then decreased slowly and became stable after 4320~5040 h dissolution. The solution phosphate concentrations had an evolution trend similar to the solution Ca2+ concentrations. Generally, the final aqueous Ca2+ and phosphate concentrations decreased with the increasing of the Zn-Ca apatites (Figure 5). The final aqueous Zn2+ concentrations increased with the increasing when and decreased with the increasing when = 0.10~0.19 (Figure 5).

For the Zn-Ca-HA dissolution at 25°C with an initial pH of 5.60 and 9.00, the solution pH, zinc, and phosphate concentrations became stable after 5040–5760 h (Figures 4(b) and 4(c)). The solution zinc and phosphate concentrations were significantly lower than those for the Zn-Ca-HA dissolution at 25°C with an initial pH of 2.00; that is, the solubility of the Zn-Ca apatites at pH 5.60 or 9.00 was considerably smaller than that at pH 2.00 (Figure 4).

The Zn-Ca apatites dissolved in the acidic solution stoichiometrically during the early stages and then nonstoichiometrically to the end of dissolution. Commonly, the solution [Zn]/[Zn + Ca] molar ratios () decreased with increasing time and were not equal to the stoichiometric Zn/(Zn + Ca) atomic ratios of the Zn-Ca apatites () (Figure 6). During the early stage of dissolution, the solution [Zn]/[Zn + Ca] molar ratios () were nearly equal to the stoichiometric for the Zn-Ca apatites of . The solution [Zn]/[Zn + Ca] molar ratios () decreased steadily with time. The higher the , the higher the solution [Zn]/[Zn + Ca] molar ratios (). For the dissolution of the Zn-Ca apatites with = 0.10~0.19, the solution [Zn + Ca]/[P] molar ratios increased slowly with time and reached the highest values in 3~48 h; after that, the solution [Zn + Ca]/[P] molar ratios decreased steadily and became stable to the end of dissolution. For the dissolution of the Zn-Ca apatites with , the solution [Zn + Ca]/[P] molar ratios increased steadily to 2.47~3.59 after 3~48 h dissolution and reached 1.61~2.93 at the end of dissolution. The solution [Zn + Ca]/[P] molar ratios increased with the increasing when and decreased with the increasing when = 0.10~0.19. The solution [Ca]/[P] molar ratios changed with time in much the same manner as the solution [Zn + Ca]/[P] molar ratios because of the very low Zn concentrations in comparison to the Ca concentrations.

3.2.2. Dissolution Mechanism

For the (ZnxCa1−x)5(PO4)3OH dissolution in acidic solution (pH 2.00 or 5.60), the H+ consumption indicated that the adsorption of H+ ions to the phosphate groups on the (ZnxCa1−x)5(PO4)3OH surface could transform the phosphate groups from to and enhanced the dissolution [28]. Besides, the coexisting replacement of H+ for metallic cations on the solid surface could also cause a H+ consumption in the (ZnxCa1−x)5(PO4)3OH dissolution. In order to describe the H+ depletion in the (ZnxCa1−x)5(PO4)3OH dissolution comprehensively, many processes should be considered: the stoichiometric dissolution of (ZnxCa1−x)5(PO4)3OH during the early stage, the substitution of 2H+ for Zn2+ or Ca2+ on the solid surface, and the adsorption/desorption of H+ on the solid surface [29]. Additionally, the experimental conditions could significantly affect the apatite dissolution [30]. Various dissolution models for apatite have been proposed in literatures, but they consider only some specific dissolution aspects and cannot describe the dissolution process comprehensively [30].

Derived from the results of the present experiment and some previous works [30], the following coexisting steps or processes are considered in the (ZnxCa1−x)5(PO4)3OH dissolution in acidic solution:(A)Diffusion of H+ to the solution-solid interface and adsorption of H+ onto the Zn-Ca-HA surface.(B)Transformation of to on the Zn-Ca-HA surface in acidic solution.(C)Desorption of Zn2+, Ca2+, and ions from the Zn-Ca-HA surface and ion complexation.(D)Readsorption of Zn2+, Ca2+, and/or ions from solution back onto the Zn-Ca-HA surface.

In Steps (A) and (B), the solution pH was increased from 2.00 to 4.42~4.90 in 1 h due to the diffusion and adsorption of H+ ions onto the Zn-Ca-HA surface for the dissolution at 25°C with an initial pH of 2.00. In Step (C), Zn2+, Ca2+, and ions were removed from the Zn-Ca-HA surface to water solution. Many chemical reactions could happen in the apatite dissolution because of the structural complexity [7]. In comparison to Zn2+ cations, Ca2+ cations could be preferentially removed from the Zn-Ca-HA surface. Reaction (1) for the Zn-Ca-HA dissolution in aqueous acidic media could be significantly affected by the solution pH and together with the reactions (2)~(5) of protonation and complexation, which had caused an increase in the solution pH. As (ZnxCa1−x)5(PO4)3OH dissolved in water, Zn2+ cations were transformed into ZnOH+, Zn, Zn, Zn,  , ZnHPO4, and ZnH2. Ca2+ cations were transformed into CaOH+, , CaHPO4, and CaH2, and   cations were transformed into , H2, H3, Ca, CaHPO4, and CaH2, and only a small portion of zinc, calcium, and phosphate existed in the dissociation forms such as Zn2+, Ca2+, and .

In Step (D), Zn2+ and Ca2+ cations were partly readsorbed from solution onto the Zn-Ca-HA surface as an initial portion of Zn-Ca-HA dissolved and the solution Zn2+ and Ca2+ concentrations decreased as the dissolution progressed. In comparison to Ca2+ ions, Zn2+ ions were preferentially readsorbed from solution onto the Zn-Ca-HA surface, which resulted in an obvious decrease in the solution [Zn]/[Zn + Ca] molar ratios () with time. Finally, adsorption and desorption of Zn2+ and Ca2+ reached a stable state. The solution Zn2+, Ca2+, and phosphate concentrations were nearly invariable from 7200 h to 8640 h for the (ZnxCa1−x)5(PO4)3OH dissolution at 25°C with an initial pH of 2.00.

3.3. Determination of Solubility

The dissolution experiments had been carried out until the analytical uncertainty for the ion activity products calculated from the last two or three samples was less than ±0.25 log units [31]. To obtain the solubility products of the (ZnxCa1−x)5(PO4)3OH solids, the aqueous activities of the zinc, calcium, and phosphate species for the last two or three solution samples (7200 h, 7920 h, and 8640 h) were considered in the calculation. The PHREEQC simulation results indicated that the final equilibrated solutions were unsaturated with respect to any potential secondary minerals (e.g., portlandite [Ca(OH)2], lime [CaO], CaHPO4·2H2O, CaHPO4, Ca4H(PO4)3·3H2O, Ca3(PO4)2(beta); Zn(OH)2, Zn(OH)2(am), and Zn3(PO4)2·4H2O).

The (ZnxCa1−x)5(PO4)3OH dissolution and the release of Zn2+, Ca2+, and can be expressed using reaction (1). The equilibrium constant () for reaction (1) can be expressed as follows:where is the activity of the solution species.

The standard free energy of reaction () can be calculated from its byFor the dissolution reaction (1),or

Table 2 lists the calculated solubility products () for Ca-HA and Zn-Ca-HA, together with the solution pH, zinc, calcium, and phosphate analyses for the dissolution at 25°C with an initial pH of 2.00. By using the Gibbs free energies of formation for Zn2+, Ca2+, , and OH- from literatures [32], that is, , , , and , the Gibbs free energies of formation for (ZnxCa1−x)5(PO4)3OH, [(ZnxCa1−x)5(PO4)3OH], were also calculated (Table 2).

For calcium hydroxylapatite [Ca5(PO4)3OH, Ca-HA], the average value was determined to be 10–57.75±0.12 at 25°C and the Gibbs free energy of formation () was calculated to be −6310.45 kJ/mol in the present work, which were consistent with the results of many previous researches. The value for Ca5(PO4)3OH has been reported to be 10–57.65 [14], 10–57 [32], 10-58±1 [33], 10–59 [34], and 10–58.3 [35].

No reports on the solubility properties of the zinc-substituted hydroxylapatites have been found in literatures. The mean values for (ZnxCa1−x)5(PO4)3OH at 25°C decreased from 10–57.75±0.12 to 10–58.60±0.18 with the increasing from 0.00 to 0.06 and then increased from 10–58.60±0.18 (10−58.42–10−58.69) to 10–56.63±0.19 with the increasing from 0.06 to 0.16 (Table 2 and Figure 7). This tendency was consistent with the dependency of the lattice parameter on . On the other hand, the corresponding free energies of formation () increased lineally from −6310.45 kJ/mol to −5979.39 kJ/mol with the increasing from 0.00 to 0.16 (Table 2 and Figure 7).

4. Summary

Examination by using XRD and FE-SEM confirmed that no obvious variation of the Zn-Ca apatites was observed in the dissolution. The cell parameter of the Zn-Ca apatites decreased with the increasing from 0.00 to 0.08, increased with the increasing from 0.08 to 0.12, and then decreased with the increasing from 0.12 to 0.16. When , ammonium zinc phosphate (NH4ZnPO4) was also observed in the precipitates. No apatite phase formed when the [Zn]/[Zn + Ca] molar ratio in the mixed aqueous solution was greater than 0.20.

The solution concentrations of zinc, calcium, and phosphate were greatly correlated to the Zn/(Zn + Ca) molar ratios of the Zn-Ca apatites (). For the dissolution at 25°C with an initial pH of 2.00, the solids of showed a different dissolution process from the solids of = 0.10~0.16. The final solution pH values increased, and the final solution Ca2+ and phosphate concentrations decreased with the increasing . The final solution Zn2+ concentrations increased with the increasing when and decreased with the increasing when = 0.08~0.16. For the dissolution at 25°C with an initial pH of 5.60 and 9.00, the solution zinc and phosphate concentrations were significantly lower than those for the dissolution at 25°C with an initial pH of 2.00. Generally, the solution [Zn]/[Zn + Ca] molar ratios () were lower than the values of the corresponding solids. Ca2+ ions were preferentially removed from solid to solution in comparison to Zn2+ ions, while Zn2+ ions were preferentially readsorbed from solution onto the Zn-Ca-HA surface, which resulted in a significant decrease of the solution [Zn]/[Zn + Ca] molar ratios () with time.

The (ZnxCa1−x)5(PO4)3OH dissolution is considered to include four coexisting steps: diffusion and adsorption of H+ onto the solid surface, transformation of to on the solid surface, desorption and ion complexation of Zn2+, Ca2+, and , and readsorption of these ions from solution back onto the solid surface.

The mean values for (ZnxCa1−x)5(PO4)3OH at 25°C decreased from 10–57.75 to 10–58.59 with the increasing from 0.00 to 0.08 and then increased from 10–58.59 to 10–56.63 with the increasing from 0.08 to 0.16. This tendency was consistent with the dependency of the lattice parameter on . The corresponding free energies of formation () increased lineally from −6310.45 kJ/mol to −5979.39 kJ/mol with the increasing from 0.00 to 0.16.

Conflicts of Interest

The authors declare that there are no conflicts of interest regarding the publication of this paper.

Acknowledgments

The authors thank the Guangxi Collaborative Innovation Centre for Water Pollution Control and Water Safety in Karst Area for the research assistance and the financial supports from the National Natural Science Foundation of China (51638006, 41263009), the Guangxi Science and Technology Development Project (GuiKeGong14124004-3-3), and the Provincial Natural Science Foundation of Guangxi (2014GXNSFBA118054, 2012GXNSFDA053022).