Near-Zero Thermal Expansion and Phase Transitions in HfMg_1−xZn_xMo₃O₁₂

Introduction

Large difference in coefficients of thermal expansion (CTE) of materials could lead to performance deterioration and even failure of devices due to thermal stress when temperature changes abruptly or frequently. Since most materials expand on heating and contract on cooling, materials with opposite property, namely negative thermal expansion (NTE), are particularly desired for tailoring CTEs. The rediscovery of NTE in ZrW₂O₈ in a wide temperature range (Evans et al., 1996, 1997a) triggered continuous efforts on understanding the NTE phenomenon and searching for more NTE materials (Yang et al., 2007; Chen et al., 2013, 2015; Tallentire et al., 2013; Lama et al., 2014; Liu et al., 2014; Peng et al., 2014; Xiao et al., 2014; Hu et al., 2015). To date, different families of NTE materials based on various mechanisms, such as the phonon effect (Pryde et al., 1996; Wang et al., 2011; Bridges et al., 2014; Cheng et al., 2016a; Ge et al., 2016a), magnetovolume effect (Takenaka and Takagi, 2005; Qu et al., 2012; Yan et al., 2014), spontaneous ferroelectric polarization (Chen et al., 2013; Peng et al., 2016), and charge transfer (Long et al., 2009; Azuma et al., 2011; Yamada et al., 2016) have been reported. Among the materials, the family of A₂M₃O₁₂ (A = transition metal or a mixture of tetravalent and bivalent cations, M = W, Mo) have been particularly attractive, because whose NTEs go over a wide temperature range and can be tuned from low positive to large negative values due to chemical flexibility (Evans et al., 1997b; Suzuki and Omote, 2006; Wu et al., 2009, 2012, 2014; Li et al., 2011; Das et al., 2013; Miller et al., 2013; Song et al., 2014a; Liu et al., 2015; Chen et al., 2016; Cheng et al., 2016a).

In recent years, a number of novel NTE materials have been designed based on the basic structure of A₂M₃O₁₂ family, including those with a general formula ABM₃O₁₂ where A is tetravalent Hf⁴⁺ or Zr⁴⁺ and B is bivalent cation Mg²⁺ or Mn²⁺, and M is W or Mo or a combination of them (Suzuki and Omote, 2004; Baiz et al., 2008; Gindhart et al., 2008; Marinkovic et al., 2008; Song et al., 2013; Li et al., 2014, 2016, 2017; Ge et al., 2016a; Liu et al., 2018) and those with a formula ABM₂XO₁₂ where A and M are the same as in ABM₃O₁₂, B is a trivalent cation and X is P⁵⁺ or V⁵⁺ (Chen et al., 2016; Cheng et al., 2016b, 2017; Ge et al., 2016b,c). The most distinct characteristics of the materials with formula ABM₂XO₁₂ are that they exhibit NTE over a wide temperature range and intense photoluminescence in the visible range. Nearly an order higher ionic conductivity was observed for HfMgW₃O₁₂ with respect to the family A₂M₃O₁₂ (Omote et al., 2011). HfMgMo₃O₁₂ with a linear CTE of 1.02 × 10⁻⁶ K⁻¹ from 298 to 1013 K was reported by Marinkovic et al. (2008). It crystallizes in orthorhombic symmetry with space group Pnma(62) or Pna2₁ (33) and transforms to monoclinic structure at 175 K (Miller et al., 2012).

In this paper, we investigate the effects of Zn²⁺ incorporation on the structure, phase transition, thermal expansion, and vibrational properties of HfMgMo₃O₁₂. It is shown that single phase solid solution of HfMg_1−xZn_xMo₃O₁₂ can be achieved only for the compositions of x ≤ 0.5, otherwise, additional phases of HfMo₂O₈ and ZnMoO₄ appear. The monoclinic to orthorhombic phase transition temperature increases with the content of Zn²⁺ for x ≤ 0.5 so that HfMg_0.5Zn_0.5Mo₃O₁₂ crystallizes in monoclinic phase and all other samples (x ≤ 0.4) adopt orthorhombic structure at room temperature (RT). The incorporation of Zn²⁺ alters the axial CTE differently for each axis and finally results in near-zero thermal expansion (ZTE) property over wide temperature ranges with smaller thermal expansion anisotropy with respect to HfMgMo₃O₁₂. The mechanisms of Zn²⁺ incorporation on the phase transition, thermal expansion and vibrational properties are discussed.

Experimental

Analytic grade reagents of HfO₂, MgO, ZnO, and MoO₃ were mixed with stoichiometric ratios for HfMg_1−xZn_xMo₃O₁₂ with x = 0.0, 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, and 1.0. The mixtures were ground in an agate mortar for 2 h, then, pressed under 325 MPa into cylinders with diameter of 10 mm and height of 6 mm using a uniaxial tablet machine. The cylinders were sintered at 1,073 K for 5 h in a muffle furnace in air and cooled down to 300 K naturally.

The as-prepared samples were analyzed by XRD with a PANalytical X'Pert PRO X-ray Diffractometer to identify the crystalline phase. Variable-temperature X-ray powder data were collected on a Rigaku (Japan, SmartLab 3KW) diffractometer with Cu Kα (λ = 0.15405 nm) radiation. Diffraction data were collected with a step size of 0.01° in the 2θ range of 10°–120°. The sample was heated at a rate of 10 K/min and remained at each measurement temperature for 5 min before measurement. Unit cell dimensions above the phase transition temperature were determined with software of PowderX. Variable-temperature/RT Raman spectra were recorded with A LabRAM HR Evolution Raman spectrometer (France HORIBA JobinYvon S.A.A.) equipped with a Linkam THMS600 Heating and Freezing Stage (Japan Hightech) (an accuracy of ±0.1 K). The excitation wavelength is 633 nm and low excitation laser power is necessary to avoid local heating by the laser. The microstructures and energy dispersive spectra of the samples were examined with a scanning electron microscope (SEM, Model Quanta 250). The relative length changes were measured with LINSEIS DIL L75 dilatometer at the heating and cooling rates of 5 K/min.

Results and Discussion

Figure 1A shows the XRD patterns of the solid solutions of HfMg_1−xZn_xMo₃O₁₂. When x = 0.0, the diffraction peaks are corresponding to HfMgMo₃O₁₂, which adopts an orthorhombic structure with space group Pnma or Pna2₁ (Marinkovic et al., 2008). No obvious changes in the XRD patterns could be observed with increasing the content of Zn²⁺ till x = 0.4. It is reasonable to conclude that HfMg_1−xZn_xMo₃O₁₂ for x ≤ 0.4 crystallized in an orthorhombic structure. Nevertheless, some subtle changes are observed for x = 0.5, such as the weak peak appearing at about 25.6° which is characteristic for a monoclinic structure (Song et al., 2014b; Ge et al., 2016a) of ABMo₃O₁₂. HfMg_0.5Zn_0.5Mo₃O₁₂ at RT is thus identified as a monoclinic structure. The XRD patterns change obviously with further increasing content of Zn²⁺. Detailed analyses show that the newly appeared peaks correspond well to HfMo₂O₈ and ZnMoO₄ (Reichelt et al., 2000; Allen et al., 2004), respectively.

Figure 1. (A) X-ray diffraction patterns of the solid solutions of HfMg_1−xZn_xMo₃O₁₂; (B) Raman spectra of the solid solutions of HfMg_1−xZn_xMo₃O₁₂.

The above analysis is supported by Raman spectroscopic analysis (Figure 1B). The Raman spectra are consistent with each other for x ≤ 0.4 while the Raman band at 988 cm⁻¹ splits into two bands at 980 and 998 cm⁻¹ for x = 0.5 (as indicated by the black arrowhead), which is characteristic for a phase transition from higher orthorhombic symmetry to lower monoclinic symmetry (Li et al., 2011; Ge et al., 2016a) for the ABMo₃O₁₂ family. Distinct changes of the Raman spectra occur for higher content of Zn²⁺. The new Raman bands at about 175, 328, 362, 772, 943, and 1,008 cm⁻¹ correspond well to HfMo₂O₈ (Liang et al., 2008b) and that around 968 cm⁻¹ arises from ZnMoO₄ (Ahsaine et al., 2016). Both XRD and Raman analyses demonstrate that a single phase solid solution of HfMg_1−xZn_xMo₃O₁₂ is only possible for x ≤ 0.5 and additional phases of HfMo₂O₈ and ZnMoO₄ form for x ≥ 0.6. At RT, HfMg_1−xZn_xMo₃O₁₂ for x ≤ 0.4 adopt an orthorhombic structure while HfMg_0.5Zn_0.5Mo₃O₁₂ crystallizes in a monoclinic structure.

Raman spectroscopy is very sensitive to the monoclinic-to-orthorhombic phase transition (Li et al., 2011, 2016; Ge et al., 2016a). In order to get some insights into the influence of Zn²⁺ on the phase transition, we carried out temperature-dependent Raman spectral observation of HfMg_1−xZn_xMo₃O₁₂ (x ≤ 0.5) as shown in Figure 2. The XRD analyses suggest that HfMg_1−xZn_xMo₃O₁₂ (x ≤ 0.5) have similar open framework structure as HfMgMo₃O₁₂. In the orthorhombic phase, there are four molecular formulas in a unit cell, in which each MoO₄ tetrahedron sharing its four vortexes with HfO₆/MgO₆ octahedra and each HfO₆/MgO₆ octahedron shares its corners with six MoO₄ tetrahedra. Hf and Mg are alternatively aligned in the [010] direction forming a quasi-layered structure (Omote et al., 2011). The Raman modes from 1,050 to 900 cm⁻¹, from 900 to 750 cm⁻¹, from 400 to 320 cm⁻¹, and from 320 to 280 cm⁻¹ are identified as symmetric stretching (ν₁), asymmetric stretching (ν₃), asymmetric bending (ν₄), and symmetric bending (ν₂) modes in the MoO₄ tetrahedra, respectively (Liang et al., 2008a; Li et al., 2011). Figure 2A shows the temperature dependent Raman spectra of HfMgMo₃O₁₂. The most distinctive change of the Raman spectra is the disappearance of the band at about 1,001 cm⁻¹ with temperature increase from 168 to 178 K, which can be regarded as characteristic of the phase transition from low temperature monoclinic to high temperature orthorhombic structure (Li et al., 2011; Ge et al., 2016a). The phase transition temperature agrees well with the result derived from XRD analysis (Miller et al., 2012). In Figures 2B–F we present the temperature dependent Raman spectra for Zn²⁺-containing samples. It is shown that the vanishing of the characteristic Raman band for the monoclinic structure occurs in the ranges of 168–178, 203–213, 223–233, 258–268, 283–293, and 318–328 K for x = 0.0, 0.1, 0.2, 0.3, 0.4, and 0.5, respectively, demonstrating that the phase transition temperature increases with the content of Zn²⁺.

Figure 2. Temperature-dependent Raman spectra of HfMg_1−xZn_xMo₃O₁₂ with x = 0.0 (A), 0.1 (B), 0.2 (C), 0.3 (D), 0.4 (E), and 0.5 (F).

In the orthorhombic structure of HfMg_1−xZn_xMo₃O₁₂ (x ≤ 0.5), the four vortexes of each MoO₄ tetrahedron are shared with two HfO₆ and two MgO₆/ZnO₆ octahedra while each HfO₆/MgO₆/ZnO₆ octahedron shares its corners with six MoO₄ tetrahedra. Statistically, each MoO₄ tetrahedron links to 0.0, 0.2, 0.4, 0.8, and 1.0 ZnO₆ octahedron for x = 0.0, 0.1, 0.2, 0.3, 0.4, and 0.5. Since the ionic radius of Zn²⁺ is 74 pm which is slightly larger than that of Mg²⁺ (72 pm), large lattice distortion and increase in phase transition temperature is not expected if only the ionic radius is considered. The experimentally observed obvious increase in phase transition temperature is therefore attributed to the difference in electronegativity between Zn²⁺ (1.65 Pauling) and Mg²⁺ (1.31 Pauling). Replacement of Mg²⁺ by Zn²⁺ causes an increase in electronegativity at theZn²⁺-cation side and a decrease in the effective negative charge on oxygen, and hence a decrease in the oxygen-oxygen repulsion. With increasing the content of Zn²⁺, oxygen-oxygen attractive forces increase, causing the network collapse transition to occur at higher temperatures (Evans et al., 1997b).

Figure 3 shows the relative length changes of sintered cylinders with increasing temperature measured by dilatometry. All the samples for x ≤ 0.6 exhibit abrupt length increase around the temperature of monoclinic to orthorhombic phase transition. The phase transition temperature increases with increasing the content of Zn²⁺ except the one for x = 0.6 whose phase transition temperature is lower than that of x ≤ 0.5 due to the generation of HfMo₂O₈ and ZnMoO₄. In this case, it can be deduced that the real content of Zn²⁺ in HfMg_1−xZn_xMo₃O₁₂ is lower than HfMg_0.5Zn_0.5Mo₃O₁₂. These results comply well with the above Raman spectroscopic analyses. The CTEs are calculated from the relative length change and shown in the Table 1. It indicates that all the single phase samples present excellent near-ZTE property above the phase transition temperature. It is interesting to notice that even for the multi-phase samples for x = 0.8 and 1.0, a near-ZTE property in a wide temperature range are realized. However, in this paper we focus on the effect of Zn²⁺ incorporation on the structure and properties of the single phase HfMg_1−xZn_xMo₃O₁₂.

Figure 3. Variations of relative length of the solid solutions of HfMg_1−xZn_xMo₃O₁₂ with increasing temperature from 170 to 673 K.

Table 1. Values of CTEs of HfMg_1−xZn_xMo₃O₁₂.

In order to get insight into the axial thermal expansion property, we carried out temperature-dependent powder XRD measurements for the samples of x = 0.2 and 0.3. For comparison, variable-temperature XRD data of HfMgMo₃O₁₂ were also collected. Figure 4A shows the selected temperature-dependent XRD patterns for HfMg_0.7Zn_0.3Mo₃O₁₂ at different temperatures. It is obvious that its XRD pattern changes distinctively around 225 K, which is attributed to the phase transformation from lower temperature monoclinic to higher temperature orthorhombic structure. Lattice constants and cell volume at each temperature are calculated and given in Figure 4B. It is evident that the a- and c-axes contract while the b-axis expands with increasing temperature. The CTEs for the a-, b-, and c-axes and volume are calculated to be ${\alpha }_a=-2.70\times 10^{-6}{K}^{-1}$, ${\alpha }_b=5.30\times 10^{-6}{K}^{-1}$, ${\alpha }_c=-1.72\times 10^{-6}{K}^{-1}$, ${\alpha }_V=0.86\times 10^{-6}{K}^{-1}(350-573K)$, respectively. This gives rise to a linear CTE ${\alpha }_l=0.29\times 10^{-6}{K}^{-1}$. Similar axial thermal expansion behaviors are obtained for HfMg_0.8Zn_0.2Mo₃O₁₂ from temperature dependent XRD measurements (not shown here). The changes of its lattice constants and volume with temperature are given in Figure 4C. The CTEs for the a-, b-, and c-axes and volume are calculated to be ${\alpha }_a=-2.26\times 10^{-6}{K}^{-1}$, ${\alpha }_b=5.21\times 10^{-6}{K}^{-1}$, ${\alpha }_c=-1.80\times 10^{-6}{K}^{-1}$, ${\alpha }_V=1.12\times 10^{-6}{K}^{-1}$, respectively, corresponding to a linear CTE ${\alpha }_l=0.37\times 10^{-6}{K}^{-1}(350-573K)$. These results are consistent with the values measured by dilatometry, confirming HfMg_0.7Zn_0.3Mo₃O₁₂ and HfMg_0.8Zn_0.2Mo₃O₁₂ being intrinsically ZTE materials. Figure 4D shows the changes of lattice constants and volume of HfMgMo₃O₁₂ with temperature.

Figure 4. (A) Variable-temperature XRD patterns of HfMg_0.7Zn_0.3Mo₃O₁₂ (The XRD patterns were measured for 2θ = 10°-120°, here only a small range is shown for clarity); Changes of lattice constants and volume of HfMg_1−xZn_xMo₃O₁₂ with temperature: (B) for x = 0.3, (C) for x = 0.2, (D) for HfMgMo₃O₁₂.

Considering the fact that the intrinsic linear CTE of HfMgMo₃O₁₂ is 1.02 × 10⁻⁶ K⁻¹ (Marinkovic et al., 2008), it is reasonable to conclude that the incorporation of Zn²⁺ reduces the linear CTE and results in near ZTE of HfMg_1−xZn_xMo₃O₁₂. A comparison of the axial CTEs for HfMgMo₃O₁₂, HfMg_0.8Zn_0.2Mo₃O₁₂, and HfMg_0.7Zn_0.3Mo₃O₁₂ are given in Table 2. It is found that partial substitution of Mg²⁺ by Zn²⁺ leads to a significant reduction of the CTE in the b-axis and an increase of the NTE in the c-axis, resulting in a near ZTE and lower anisotropy of thermal expansion in the Zn-containing compounds with respect to HfMgMo₃O₁₂ (see Table 2). The anisotropy of thermal expansion is defined as the maximum difference in the axial thermal expansion coefficients (Srikanth et al., 1992; Miller et al., 2013). The near zero linear thermal expansion and lower anisotropy property of the Zn-containing compounds suggest that they could withstand higher thermal shock resistance.

Table 2. Intrinsic thermal expansion coefficients (α) for HfMg_0.7Zn_0.3Mo₃O₁₂ and HfMg_0.8Zn_0.2Mo₃O₁₂ as obtained from variable-temperature XRD and literature and experimental values for HfMgMo₃O₁₂.

The difference in the linear CTEs measured by XRD and dilatometry could be understood by the microstructural effects. In contrast to XRD measurement which gives the thermal expansion property of cell lattice, dilatometry reveals the bulk thermal expansion property, including both intrinsic (thermal expansion of a material arising from the lattice dynamics) and extrinsic (thermal expansion related to microstructures such as texture, grain size, grain boundaries, poses, and microcracks) effects. The difference measured by the two methods reflects the extrinsic effect in the sintered bulk, which, on heating, can add a small negative component to the intrinsic linear expansion coefficients. Generally speaking, a smaller difference suggests a better sintered quality of the bulk material which is desired for most applications. The absolute differences for the Zn-containing compounds (Δ${\alpha }_l=0.42\times 10^{-6}{K}^{-1}$ for x = 0.2 and Δ${\alpha }_l=0.38\times 10^{-6}{K}^{-1}$ for x = 0.3) are obviously smaller than that for HfMgMo₃O₁₂ (Δ${\alpha }_l=1.17\times 10^{-6}{K}^{-1}$). It means that partial substitution of Mg²⁺by Zn²⁺ in HfMgMo₃O₁₂ could improve the sintering performance of the material and minimize the possible contributions of extrinsic effects. The analysis is supported by microstructural observation.

Figure 5 shows the SEM images of HfMg_1−xZn_xMo₃O₁₂ ceramics with x = 0.0 (Figure 5A), 0.1 (Figure 5B), 0.2 (Figure 5C), 0.3 (Figure 5D), 0.4 (Figure 5E), and 0.5 (Figure 5F). The micro morphology of the sample for x = 0.1 is dominated by well crystallized polyhdra or truncated polyhedra and the average particle size is obviously smaller than that of HfMgMo₃O₁₂. With increasing the content of Zn²⁺, the polyhedra become more rounded. Compared to HfMgMo₃O₁₂, the incorporation of Zn²⁺ seems to lead to less pores in the ceramic bodies and pores can hardly be found in the solid solutions for x = 0.1–0.3. It illustrates that proper amount incorporation of Zn²⁺ favors the formation of uniform distribution of particles and efficient reduction of porosity in the sintered body.

Figure 5. SEM images of HfMg_1−xZn_xMo₃O₁₂ with x = 0.0 (A), 0.1 (B), 0.2 (C), 0.3 (D), 0.4 (E), and 0.5 (F).

The schematic diagram of HfMgMo₃O₁₂ is given in Figure 6A to help us understand the mechanism of the phenomenons. In HfMg_1−xZn_xMo₃O₁₂ (x ≤ 0.5), Zn²⁺ is expected to substitute for Mg²⁺ due to the same valence and similar cation radius and each ZnO₆ octahedron shares all its corners with six MoO₄ tetrahedra. In order to see the bond strength changes induced by local electronic environment upon substitution of Zn²⁺for Mg²⁺, we show in Figure 6B the Raman spectra of the stretching region for x = 0.0, 0.1, 0.2, 0.3, 0.4, and 0.5. It is obvious that all the stretching modes shift successively to lower wavenumbers with increasing the content of Zn²⁺, indicating a softening of the Mo-O bonds upon incorporation of Zn²⁺. Once an Mg²⁺ is replaced by Zn²⁺, the local electronic equilibrium around the MoO₄ tetrahedron is broken. Zn²⁺ has obviously a higher electronegativity and ability to drag electrons to the ZnO₆ octahedron from its connected six MoO₄ tetrahedra than Mg²⁺, resulting in the weakening of the Mo–O bonds. The differences in ionic radius and electronegativity could cause a slight rotation of the connected polyhedra and hence the M-O-M linkages. This is probably the reason that the positive expansion of the b-axis is pronouncedly reduced and the NTE in the c-axes become more negative, resulting hence in a lower anisotropy in thermal expansion and near-zero CTEs of the Zn-containing compounds. Due to the large difference in electronegativity between Zn²⁺ (1.65 Pauling) and Mg²⁺ (1.31 Pauling), the more of Zn²⁺ is incorporated, the more of the MoO₄ tetrahedra get distorted as revealed by Raman spectroscopy, resulting in larger distortion and instability of the lattice. When the forming energy for the single phase exceeds that for the multi-phases, then the multi-phases form.

Figure 6. (A) Schematic diagram of HfMgMo₃O₁₂ building block (red sphere indicates oxygen atom); (B) Raman spectra of HfMg_1−xZn_xMo₃O₁₂ in the high wavenumber regions; (C) Temperature-dependent Raman spectra of HfMg_0.9Zn_0.1Mo₃O₁₂ in the low wavenumber region.

Figure 6C shows the temperature-dependent Raman spectra of HfMg_0.9Zn_0.1Mo₃O₁₂ in the low wavenumber region. Obvious change of the Raman spectra occur between 203 and 213 K, corresponding to the monoclinic to orthorhombic phase transition, such as the appearance of new Raman modes at about 27,44, 86, 159, and 182 cm^–1. The modes at about 44 and 86 cm^–1 are split into two or three modes in the low temperature phase and become degenerated in the high temperature phase. The low wavenumber modes arise from the external librational and translational vibrations of the connected octahedra–tetrahedra, or the librational and translational motions of metal ions in the Hf(Mg/Zn)–O–Mo linkages, which can also be regarded as the transverse vibrations of the bridging oxygen from the point of view of relative movement. Such an harmonic vibrations along with the distortion of the polyhedra are believed to be the origin of the NTE in the open frame work structure since they bring the two end atoms closer upon heating (Evans, 1999; Ding et al., 2008; Marinkovic et al., 2009; Wang et al., 2013).

Conclusion

Solid solutions of HfMg_1−xZn_xMo₃O₁₂ with near-ZTE are successfully synthesized by solid state reaction and the effects of Zn²⁺ incorporation on the phase formation, thermal expansion, phase transition, and vibrational properties and micro-morphologies are investigated by XRD, dilatometry, Raman spectroscopy, and SEM. It is shown that (i) single phase formation is only possible for x ≤ 0.5, otherwise, additional phases of HfMo₂O₈ and ZnMoO₄ generate; (ii) HfMg_1−xZn_xMo₃O₁₂ crystallize in an orthorhombic structure for x ≤ 0.4 and in a monoclinic structure for x = 0.5 at RT; (iii) The phase transition temperature from monoclinic to orthorhombic structure increases with the content of Zn²⁺, which occurs within 168–178, 203–213, 223–233, 258–268, 283–293, and 318–328 K for x = 0.0, 0.1, 0.2, 0.3, 0.4, and 0.5, respectively; (iv) The incorporation of Zn²⁺ leads to an pronounced reduction in the positive expansion of the b-axis and an enhanced NTE in c-axes, making the Zn-containing materials exhibit near-ZTE over a wide temperature range and lower anisotropy in thermal expansion in the orthorhombic phase; (v) Replacement of Mg²⁺ by Zn²⁺ breaks the local electronic equilibrium around the MoO₄ tetrahedron and weakens the Mo–O bonds, leading to obvious red shifts of all the Mo–O stretching modes with increasing the content of Zn²⁺ due to obviously higher electronegativity of Zn²⁺ than Mg²⁺; (vi) The incorporation of Zn²⁺ improves sintering property of samples, minimizing the possible contributions of extrinsic effects such as pores, which is preferred for most applications.

Author Contributions

EL conceived the idea and supervised the research. SL and RS are in charge of the synthesis and part measurements of the materials. XG, HY, and DC are in charge of the thermal expansion and Raman characterization. JG and MC are in charge of the XRD characterization and structural analyses. SL, XG, and EL are in charge of the manuscript preparation.

Conflict of Interest Statement

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Acknowledgments

This work was supported by the National Natural Science Foundation of China (No. 11574276; No. 51302249).

References