Abstract
This study systematically investigates the high-temperature oxidation behaviour of (Zr,Ti,W)C-Me
x
B
y
multiphase ceramics fabricated by spark plasma sintering, focusing on their performance in static air at 1000°C
Introduction
Ultra-high-temperature ceramics (UHTCs) are in urgent demand for advanced equipment in aerospace, nuclear energy systems and hypersonic vehicles, yet their oxidation resistance in extreme environments still faces severe challenges.1–3 Refractory metal carbides (e.g., ZrC, TiC, WC), with their high melting points, excellent high-temperature strength and good wear resistance, are considered promising candidates for UHTCs, with broad application prospects in aerospace thermal protection, nuclear fusion reactors, cutting tools and other fields.4–8 However, these materials are prone to severe oxidation in high-temperature oxygen environments, leading to structural degradation and performance failure. Particularly, single-component carbides exhibit poor high-temperature oxidation resistance,9–12 significantly limiting their engineering application prospects.9,10,13–15
To enhance the oxidation resistance of refractory carbides, a synergistic strengthening strategy commonly employed involves multiple solid solutions combined with second-phase incorporation. Our previous study utilised a cation-anion dual solid solution (Zr,Ti)(C,N) and introduced SiC through an in-situ reaction, which significantly improved the oxidation resistance of the composite material at 850°C–950°C. 16 Kang et al. 17 enhanced the material's oxidation resistance by employing (Ti,W,Mo,Ta)(C,N) multi-component solid solution and second-phase composite effects.
Borides such as ZrB2 exhibit certain oxidation resistance and thermal shock resistance at high temperatures.18,19 ZrB2 reacts with O2 at temperatures above 650°C to form ZrO2 and B2O3. 20 During oxidation, a ZrO2 framework and B2O3 glass phase are formed, with B2O3 effectively filling cracks and pores within the oxide layer, thereby hindering further oxygen diffusion.1,20–22 Consequently, ZrB2 is considered an effective means to enhance the oxidation resistance of refractory metal carbides and is often used synergistically with SiC to improve the oxidation resistance of ZrC ceramics. Studies on the high-temperature oxidation behaviour of the ZrC-ZrB2-SiC system23,24 indicate that oxidation below 1200°C is reaction-controlled, while above 1300°C it becomes diffusion-controlled. At 1500°C, a ZrO2 oxide layer with large pores forms, leading to a decline in oxidation resistance.
In our previous study, 25 (Zr,Ti,W)C-Me x B y multiphase ceramics were prepared by adding varying contents of ZrB2 into a (Ti,W)C matrix, achieving synergistic enhancement in both strength and toughness of the material. On the one hand, the (Zr,Ti,W)C solid solution integrates the advantages of each component, improving material strength while maintaining a high melting point. On the other hand, ZrB2 can react in-situ with the matrix to form platelet-like toughening phases such as TiB2 and W2B₅. However, there remains a lack of systematic understanding regarding the microstructural evolution, oxide layer formation mechanisms and kinetic behaviour of (Zr,Ti,W)C-Me x B y multiphase ceramics during oxidation at 1000°C–1100°C (the primary operating temperature range for divertor components in nuclear fusion devices). Particularly, the influence of ZrB2 content on the phase composition, structural stability and protective performance of the oxide layer under coupled vapourisation-diffusion effects remains unclear, warranting further investigation. Notably, B2O3 exhibits significantly reduced protective capability due to its decreased viscosity and vapourisation at temperatures around 1000°C, 26 while WO3 generated from the oxidation of W in the material also undergoes intense vapourisation within this temperature range, 27 collectively leading to a fundamental transformation in the oxidation mechanism. Therefore, studying the oxidation behaviour in this critical 1000°C–1100°C temperature range is decisive for elucidating the oxidation mechanisms of (Zr,Ti,W)C-Me x B y multiphase ceramics.
In this study, two characteristic temperature points of 1000°C and 1100°C were selected to systematically investigate the high-temperature oxidation behaviour of (Zr,Ti,W)C-Me x B y ceramics in static air. The phase evolution and microstructural characteristics of the oxide layer were revealed. Combined with thermogravimetric analysis and oxidation kinetic fitting, their oxidation mechanisms and variations in apparent activation energy were elucidated. This research aims to provide theoretical foundations and experimental support for the design and application of components such as nuclear fusion device divertors or oxidation-resistant structural parts.
Experimental procedure
(Zr,Ti,W)C-Me x B y multiphase ceramics were fabricated by Spark Plasma Sintering (SPS) according to our previous study. 25 The experiment used (Ti0.602W0.193)C powder (purity > 99.2 wt%, average particle size 2.05 μm) and ZrB2 powder (purity > 99.9%, average particle size 1 μm) as raw materials. The powders were placed in a polyethylene ball milling jar according to the designed ratio, with silicon nitride (Si3N4) as the grinding medium and a ball-to-powder mass ratio of 5:1, and milled at 150 rpm for 24 h. After ball milling, the powder was sieved through a 120-mesh screen. Subsequently, an SPS system (SPS-10T-10-1 V, China) was used for sintering: the powder was first cold-pressed, then heated to 1800°C at a heating rate of 100°C/min under a vacuum environment of <0.1 Pa, held for 10 min and subjected to a uniaxial pressure of 50 MPa throughout the sintering process. After sintering, the samples were rapidly cooled to room temperature with the furnace. The relative densities of the sintered ceramics, as reported in our previous study, 25 were 88.8%, 96.8%, 98.7% and 98.3% for TW0ZB, TW40ZB, TW50ZB and TW60ZB, respectively. Table 1 presents the composition design, sintering parameters and relative density of the (Zr,Ti,W)C-Me x B y multiphase ceramics used in this study.
The compositional design and sintering schedule and relative density for (Zr,Ti,W)C-Me x B y multiphase ceramics.
The (Zr,Ti,W)C-based multiphase ceramic samples were cut into rectangular test blocks measuring 10 mm × 10 mm × 4 mm. The surfaces were sequentially ground and polished using SiC sandpaper with varying grit sizes. The pre-oxidation masses were precisely measured using an electronic balance (FA1004B, accuracy 0.1 mg). The samples were placed vertically. A large crucible with an inner diameter of 12 mm was used as the outer container, inside which a small crucible with an inner diameter of 7.5 mm was placed. The sample was then vertically placed into the large crucible. This configuration minimises the contact area between the sample and the crucible, thereby effectively reducing experimental errors. The mass measurements were performed as follows: an initial mass was recorded before oxidation; during the oxidation process, samples were removed after 1, 3, 5, 7, 9, 11, 13 and 15 h of exposure, cooled to room temperature and then weighed. To ensure data reliability, each sample was measured at least three times at each time point, and the most accurate value was taken as the final result. X-ray diffraction (XRD; X'Pert, Netherlands) was employed to analyze the phase composition of the sample surface layers and detached oxide powders, with the following test parameters: diffraction angle range 10°–90°, scanning speed 10°/min. Following oxidation, the samples were embedded in resin, and their cross-sections were sequentially ground and polished for microscopic examination using a scanning electron microscope (Hitachi SU5000, Japan) to observe the oxide layer microstructure, accompanied by energy-dispersive spectroscopy (EDS) for elemental composition analysis. Thermodynamic calculations to explore reaction energy changes were performed using HSC Chemistry 6.0 software.
Results and discussion
Oxidation kinetics
Figure 1 presents the thermogravimetric curves of the ceramics after oxidation at 1000°C and 1100°C. The oxidation weight gain of ceramic samples with varying ZrB2 contents shows significant differences. At 1000°C, the order of oxidation weight gain is: TW > TW40ZB > TW50ZB > TW60ZB, while at 1100°C, the sequence changes to: TW40ZB > TW50ZB > TW > TW60ZB. The TW60ZB sample exhibits the smallest oxidation weight gain at both temperatures. Combined with the previous analysis of oxide layer depth, TW60ZB demonstrates the shallowest total oxide layer depth, the thickest dense inner oxide layer and the minimal mass change, collectively indicating its superior oxidation resistance.

Evolution of weight gain per unit area (mg/mm2) over time during oxidation tests in static air for ceramics TW0ZB, TW40ZB, TW50ZB and TW60ZB at 1000°C and 1100°C.
The weight gain trend was successfully fitted using the kinetic equation (1):
where Δm represents the weight gain, S is the surface area of the ceramic, t denotes the oxidation time, kn is the oxidation rate constant, n signifies the oxidation exponent and c corresponds to the constant oxidation before reaching isothermal conditions.
The value of n was determined by fitting experimental data using equation (1). As shown in Table 2, the exponent constants of the fitting function can be categorised into two types: the first involves n fluctuating around 1 (linear), while the second shows n fluctuating around 0.5 (parabolic). Linear and parabolic kinetic models are commonly used simplified models to describe the oxidation behaviour of pure ceramics. During pure metal oxidation, an oxide layer typically forms covering the surface; however, in ceramic oxidation processes, the presence of pores, cracks or volatile products within the oxide layer often leads to the formation of non-protective oxide layers.
Values of oxidation exponent (n) and coefficient of determination (R2) obtained by fitting curves using equation (1) for ceramics oxidised at 1000°C and 1100°C.
At both temperatures, the TW samples exhibit oxidation kinetic fitting exponents n close to 1, indicating their oxidation behaviour follows linear kinetics. This is attributed to the single-phase (Ti,W)O2 oxidation product of TW, which forms a porous oxide layer structure under the combined action of CO, CO2 and O2, primarily due to WO3 volatilisation. Although the TiO2 therein provides some protective effect to the substrate, the overall oxidation resistance remains limited. In contrast, TW40ZB follows parabolic kinetics (n ≈ 0.5) at 1000°C, significantly differing from the linear behaviour of TW, demonstrating that ZrB2 addition markedly alters the oxidation mechanism and oxidation resistance at moderate temperatures. Parabolic kinetics typically indicates the formation of protective oxide layers, where the growth rate is inversely proportional to the oxide layer thickness, gradually slowing with oxidation progress.
At 1000°C, TW40ZB, TW50ZB and TW60ZB all exhibit n values close to 0.5, confirming their oxidation weight gain follows parabolic law and the process is diffusion-controlled. However, deviations exist between the actual system and ideal models: oxygen diffusion through grain boundaries, presence of defects and CO/CO2 gas generation from carbon oxidation in the material can all lead to porous oxide layers, accelerating oxidation and causing kinetic behaviour to deviate from the ideal parabola. Notably, at 1100°C, all ZrB2-containing samples show n values close to 1, with oxidation behaviour transitioning to linear kinetics. This suggests the protective oxidation mechanism dominant at 1000°C weakens with rising temperature. The cause may be aggravated WO3 volatilisation at high temperatures, forming numerous gas diffusion channels that render the oxide layer structure porous, degrading oxidation resistance.
The oxidation rate constant (Kp) for each ceramic sample at different temperatures was calculated from the relationship between the square of oxidation weight gain and oxidation time, with the corresponding values listed in Table 3. The results show that at the same oxidation temperature, the oxidation rate constant of the samples is significantly correlated with ZrB2 content, confirming that ZrB2 addition effectively enhances the oxidation resistance. For example, at 1000°C, the Kp value of TW ceramic is 16 times higher than that of TW60ZB. With increasing ZrB2 content, the Kp value gradually decreases, further indicating enhanced oxidation resistance of the material.
Values of parabolic rate constant (Kp) obtained from slope of linear fitting and coefficient of determination (r) for ceramics oxidised at 1000°C and 1100°C.
The apparent activation energy (Ea) of the oxidation process for each system was determined using the Arrhenius equation (equation (2)) and the coefficient of determination (r) of ceramic samples, including TW, TW40ZB, TW50ZB and TW60ZB. The apparent activation energy is a key parameter reflecting the high-temperature oxidation resistance of ceramic materials. In this study, the multiphase oxide layer formed by (Zr,Ti)O2 and (Ti,Zr)O2 solid solutions doped with ZrO2 and TiO2 effectively inhibited the diffusion of oxygen ions and matrix elements in the ceramic, reduced the overall oxidation rate and thereby significantly altered the apparent activation energy of the oxidation process.
According to the Arrhenius equation, the oxidation rate constant Kp at thermodynamic temperature T can be expressed as:
where A is the pre-exponential or frequency factor; Ea is the apparent activation energy; R is the universal gas constant; and T is the absolute temperature.
The apparent activation energy was calculated by plotting the relationship curve between lnKp and 1/T (Figure 2) and determining its slope, with the results summarised in Table 4. The results demonstrate that ZrB2 addition significantly enhances the apparent activation energy of oxidation for TW-based ceramics. Within the 1000°C–1100°C temperature range, the oxidation resistance of ceramics improves with increasing activation energy. Although ZrB2-containing ceramics exhibit similar activation energy values, comprehensive analysis of total oxide layer depth, inner oxide layer thickness and oxidation weight gain confirms TW60ZB possesses the optimal overall oxidation resistance. Compared to unmodified TW, its oxidation resistance improves by 355%.

Arrhenius plots showing the temperature dependence of oxidation rate constants for four ceramic materials at oxidation temperatures of 1000°C and 1100°C.
Apparent activation energy (Ea) of TW0ZB, TW40ZB, TW50ZB and TW60ZB ceramics oxidised at 1000°C–1100°C.
Oxidation substance analysis
Figure 3 shows the XRD patterns of (Zr,Ti,W)C-based multiphase ceramics with varying ZrB2 contents before oxidation. The reaction between (Ti,W)C and ZrB2 generated (Ti,Zr)B2, (Ti,W,Zr)C and (Zr,Ti,W)C solid solutions. In the 40ZB sample, in addition to the diffraction peaks of these solid solutions, residual (Ti,W)C diffraction peaks were observed; in the 50ZB sample, (Ti,W)C reacted completely with ZrB2, exhibiting diffraction peaks of WB, W2B₅, W2B and WB₄; while in the 60ZB sample, residual ZrB2 diffraction peaks were detected.

X-ray diffraction (XRD) pattern of (Zr,Ti,W)C-Me x B y ceramics before oxidation.
Figure 4 presents backscattered electron images of the as-sintered ceramics. As shown, the black phase is identified as (Ti, Zr)B2, while the grey phase corresponds to (Zr,Ti,W)C. However, further analysis reveals that the grey phase actually consists of two distinct phases: (Zr,Ti,W)C formed via in-situ reaction, and (Ti,W,Zr)C resulting from the solid solution of Zr into the residual (Ti,W)C. These two phases are distributed among plate-like (Ti, Zr)B2 grains. 25

Backscattered electron (BSE) images of (a) TW40ZB, (b) TW50ZB and (c–d) TW60ZB ceramics before oxidation.
Figure 5 shows the macroscopic morphologies of (Zr,Ti,W)C-based multiphase ceramics after oxidation at 1000°C and 1100°C for durations of 1, 3, 5, 7, 9, 11, 13 and 15 h, respectively. During oxidation, with prolonged oxidation time, the sample without ZrB2 addition (0ZB) exhibited cracking at the edges of its prismatic structure after oxidation, while ZrB2-added samples maintained overall structural integrity. Additionally, the surface oxidation products gradually turned yellow with increasing oxidation time, due to the oxidation of TiC to rutile-type TiO2 and the oxidation of ZrC to monoclinic ZrO2 (m-ZrO2). 28 This is confirmed by the strong diffraction peaks observed near 27° and 36° from TiO2 in the temperature range of 800°C–1200°C, 29 as well as the oxidation of ZrC to ZrO2 and CO2 above 900°C. 30 With increasing ZrB2 content, cracking and spalling of the oxide layer were significantly reduced, with the 60ZB sample showing the lowest degree of cracking. From the macroscopic morphology, it can be observed that even under the highest experimental temperature and duration conditions, the TW60ZB effectively resists severe oxidation damage. As shown in Figure 5(b), after 15-h oxidation at 1100°C, compared to other samples, the TW60ZB sample maintained the higher structural integrity, demonstrating its superior high-temperature oxidation resistance.

The morphology of (Zr,Ti,W)C-Me x B y ceramics after oxidation at (a) 1000°C and (b) 1100°C.
To determine the phase composition of the oxide layer in (Zr,Ti,W)C-based multiphase ceramics, the detached oxide layer was ground into powder and subjected to XRD analysis, with the results shown in Figure 6. Figure 6(a) reveals that the primary phases of all detached oxide layers are tetragonal TiO2 (t-TiO2, PDF#99-000-3236) and monoclinic ZrO2 (m-ZrO2, PDF#04-005-7378). Further analysis of Figure 6(b) and (d) demonstrates that in ZrB2-added ceramic samples, with Zr atoms dissolving into the TiO2 unit cell, the lattice parameter of TiO2 increases, causing its diffraction peaks to shift to lower angles; Simultaneously, Ti atoms dissolving into the ZrO2 unit cell reduce the lattice parameter of ZrO2, inducing its diffraction peaks to shift to higher angles. In contrast, for ZrB2-free pure TW ceramic samples, W atoms dissolving into the TiO2 unit cell decrease the lattice parameter of TiO2, resulting in high-angle shifts of its diffraction peaks.

(a) X-ray diffraction (XRD) patterns and (b) expanded patterns of detached oxide layers of (Zr,Ti,W)C-Me x B y ceramics after oxidation at 1000°C for 15 h; (c) XRD patterns and (d) expanded patterns of detached oxide layers of (Zr,Ti,W)C-Me x B y ceramics after oxidation at 1100°C for 15 h.
After 15-h oxidation at 1000°C, XRD analysis revealed that the sample without ZrB2 addition was oxidised to form (Ti,W)O x solid solution. When the oxidation temperature was increased to 1100°C and maintained for 15 h, this sample underwent further oxidation, producing TiO2 and WO3. Due to the high-temperature volatility of WO3, only TiO2 was detected in the XRD pattern. For TW60ZB ceramics, the initial phases consisted of (Ti,Zr)B2, (Ti,W,Zr)C, (Zr,Ti,W)C, as well as WB, W2B5, W2B, WB₄ and unreacted residual ZrB2 . After 15-h high-temperature oxidation, these phases ultimately transformed into (Ti,Zr)O2, B2O3, WO3, ZrO2 and (Zr,Ti)O2. Among these, WO3 and B2O3 partially volatilised at high temperatures, while the remaining oxides collectively formed a multiphase oxide layer. Compared to a single-phase oxide layer, this multiphase oxide layer exhibited a denser structure, superior oxidation resistance and provided more effective protection to the underlying substrate.31–33
To evaluate the oxidation tendencies of each phase, Figure 7 shows the standard Gibbs free energy curves of reactions between TiC, ZrC, WC, ZrB2, TiB2, WB, W2B5, W2B and O2 as a function of temperature. Thermodynamic analysis of the oxidation products indicates that the reaction forming CO2 is more favourable in terms of Gibbs free energy compared to those forming CO or free carbon.

Curves of Gibbs free energy change with temperature for different reactions.
The (Zr,W,Ti)C solid solution in this study has a crystal structure similar to ZrC, therefore its oxidation process can be described by referring to the oxidation model of ZrC, with the reaction equation expressed as equation (3):
To simplify the thermodynamic expression, the influences of trace Ti and W elements can be neglected, and the competitive reaction pathways between Zr and the oxidation processes of TiC, WC and tungsten borides (W x B y ) can also be ignored. Similarly, the (Ti,W,Zr)C solid solution, which has a crystal structure analogous to TiC, can be treated using the oxidation model of TiC.
The Gibbs free energy variation trend clearly indicates that the oxidation reaction of W2B₅ to form WO3 and B2O3 proceeds most readily, with the reaction equation expressed as equation (4):
Among all the phases considered, WC exhibits the highest thermodynamic stability,
34
with its oxidation reaction expressed as equation (5):
Therefore, based on the thermodynamic results, within the temperature range of 1000°C–1100°C, the ease of oxidation reactions for each phase in this ceramic system follows the order: W2B5> ZrB2 > W2B > TiB2 > ZrC > TiC > WB > WC.
Oxidation microstructure characterisation
Figure 8 presents the cross-sectional microstructures and elemental distributions of the TW60ZB ceramic after oxidation at 1000°C and 1100°C for 15 h. At 1000°C (Figure 8(a) and (b)), irregularly distributed pores are observed within the oxide layer, likely originating from diffusion channels formed by the volatilisation of WO3 and B2O3. 35 Energy-dispersive spectroscopy mapping reveals partial overlap of Zr and Ti between the inner and outer oxide layers, indicating the formation of Zr-rich (Zr,Ti)O2 and Ti-rich (Ti,Zr)O2 solid solutions. The W content in the outer layer is significantly lower than that in the inner layer, which is attributed to the volatilisation of WO3 – formed by the reaction of W x B y with O2 – in the temperature range of 800°C–900°C 36 ; the dense inner oxide layer suppresses WO3 volatilisation, retaining more WO3 within, whereas the vigorous reaction in the outer layer leads to substantial WO3 loss. At 1100°C (Figure 8(c) and (d)), the overlap of Zr and Ti distributions persists, confirming the presence of both Zr-rich (Zr,Ti)O2 and Ti-rich (Ti,Zr)O2 phases in the inner and outer layers. Compared with the results at 1000°C, the W content further decreases due to the intensified oxidation reaction at elevated temperatures, which promotes more severe WO3 volatilisation. During this process, a volatilisation-induced migration mechanism may transport some TiO2 to the oxide surface via physical or transient chemical adsorption37–40; as WO3 continuously volatilises into the atmosphere, the transported TiO2 remains on the surface, resulting in significant TiO2 enrichment in the outermost oxide layer. A comparison of Zr content between the inner and outer layers shows that the outer layer is richer in Zr. This phenomenon may be related to the oxidation of ZrB2 to form B2O3: B2O3 fills the voids in ZrO2 and, during volatilisation at high temperatures, carries ZrO2 toward the outer oxide layer, thereby enriching Zr in the outer region.41,42

Cross-sectional scanning electron microscope (SEM) and energy-dispersive spectroscopy (EDS) images of the oxidation outer layer (a) and oxidation inner layer (b) of TW60ZB multiphase ceramic after 15-h oxidation at 1000°C; the oxidation outer layer (c) and oxidation inner layer (d) after 15-h oxidation at 1100°C.
It should be clarified that B2O3 exists as an amorphous glassy phase in the temperature range of 1000°C–1100°C. According to the literature, B2O3 has a melting point of approximately 450°C, 43 above which it transforms into a low-viscosity glassy phase with good fluidity and filling ability, allowing it to effectively infiltrate pores and microcracks within the oxide layer to form a dense protective layer.43,44 Although the volatilisation rate of B2O3 increases with temperature in this range, it predominantly remains in the glassy state within the oxide layer until complete volatilisation. 45 Therefore, the conclusion that B2O3 functions as a glassy phase for filling and densification within the experimental temperature range is widely accepted in the field of oxidation research on ZrB2-based ultra-high temperature ceramics.46,47
Figure 9 further identifies the phase composition of the oxide layer via point EDS. At 1000°C (Figure 9(a) and (b)), based on XRD and EDS analyses, the dark grey regions are identified as (Ti,Zr)O2, the light grey regions as (Zr,Ti)O2, the bright grey regions as WO3 and the pale grey regions as ZrO2. Additionally, a (Ti,W)O x phase is detected at this temperature. As the temperature increases to 1100°C (Figure 9(c) and (d)), the (Ti,W)O x phase disappears with the intensified volatilisation of WO3, and TiO2 is formed.

Cross-sectional scanning electron microscope (SEM) images and energy-dispersive spectroscopy (EDS) analysis results of the oxidation inner layer (a) and outer layer (b) of (Zr,Ti,W)C-based multiphase ceramic after 15-h oxidation at 1000°C; the oxidation inner layer (c) and outer layer (d) after 15-h oxidation at 1100°C.
Collectively, Figure 8 and Figure 9 reveal a distinct temperature-dependent evolution of the oxide layer structure in the TW60ZB ceramic. At 1000°C, the oxide layer exhibits a bilayer structure comprising a porous outer layer and a relatively dense inner layer. The inner layer benefits from the filling of pores and microcracks by B2O3, which enhances the density and oxygen-barrier performance of the oxide layer. Meanwhile, the formation of (Ti,Zr)O2 and (Zr,Ti)O2 solid solutions further contributes to the structural integrity of the oxide layer. The Zr-rich phase is uniformly distributed as large blocks, whereas the Ti-rich phase exists as small agglomerated particles interspersed with ZrO2 and (Ti,W)O x ; such a multiphase composite structure effectively improves the density and oxidation resistance of the oxide layer. 31 At 1100°C, the volatilisation of WO3 and B2O3 intensifies, leading to increased porosity and cracking in the outer layer and a reduction in the protective efficiency of the oxide layer, although the inner layer retains a relatively dense structure.
Based on the above microstructural observations, the multiphase mixed oxide layer enhances oxidation resistance primarily through the following mechanisms. First, the oxygen ion detouring mechanism: oxygen ions must bypass ZrO2 regions with low diffusion rates during transport, substantially lengthening the diffusion path and markedly reducing the overall oxygen diffusion rate, thereby inhibiting the oxidation reaction. This principle is analogous to the barrier network formed by second-phase particles (graphene) in the coatings studied by Wei et al., 48 which forces corrosive media to detour, prolonging the diffusion path and reducing the overall diffusion rate. Second, the grain boundary pinning mechanism: phase boundaries between different phases can effectively pin grain boundaries, suppressing grain growth of the oxide layer at high temperatures. Fine grains and tortuous grain boundaries collectively hinder the short-circuit diffusion of oxygen ions and metal ions, maintaining the stability and density of the oxide layer even under prolonged high-temperature exposure. This observation aligns with the conclusion of Mao et al., 49 who reported that phase boundaries (TiB2-ZrO2) in multilayer structures pin grain boundaries and improve oxidation resistance. Additionally, the mismatch in thermal expansion coefficients (CTE) between TiO2 and ZrO2 generates a micro-stress field near the phase interfaces, which helps deflect or terminate microcracks generated during oxidation, effectively inhibiting crack propagation and oxide layer spallation, thereby preserving the structural integrity of the oxide layer. 50
Other studies on the high-temperature oxidation behaviour of refractory carbides also indicate that mixed oxides often exhibit superior protective performance. For instance, during the high-temperature oxidation of ZrC-20TaSi2, mixed oxides such as TaZr2.75O8 formed a thick, dense oxide layer with isolated pores, effectively hindering oxygen diffusion into the substrate. 11 Similarly, in the present study, the (Zr,Ti)O2 solid solution forms a dense and continuous outer layer that protects the underlying porous ZrO2 and TiO2 layers, reducing the oxygen diffusion rate and progressively depleting the oxygen partial pressure within the composite, thereby effectively suppressing further oxidation of the substrate. 51 Under identical oxidation conditions, the multiphase oxide layer composed of (Ti,Zr)O2 and (Zr,Ti)O2 solid solutions interspersed with TiO2 and ZrO2 exhibits superior oxidation resistance compared to a single-phase solid solution oxide layer. This conclusion is further supported by extensive literature: single-phase WC develops an oxide layer depth of approximately 1.35 mm after oxidation in air at 800°C for 2 h, 52 whereas under the same conditions, the TWT sample investigated by Zheng et al. exhibits a reduced oxide layer depth of 812 μm, 53 indicating that the incorporation of MeC (Me = Ti, Ta) partially suppresses oxygen penetration through the porous WO3 layer into the substrate, indirectly confirming the beneficial role of multiphase structures in enhancing oxide layer protectiveness 54 ; Song et al. introduced reinforcing phases such as YAG (yttrium aluminium garnet) and Al2O3 to form a multiphase structure, significantly improving the oxidation resistance of ZrB2 ceramics – compared with pure ZrB2 ceramics, the resulting multiphase ceramics exhibited markedly thinner oxide layers and higher oxide layer densification 55 ; furthermore, Xia et al. reported that multiphase oxide layers typically contain glass phases (SiO2) or dense oxide phases (Al18B4O33), which can infiltrate and fill pores and cracks within the oxide layer at elevated temperatures, forming a dense protective film that substantially reduces oxygen diffusion rates and effectively blocks oxygen penetration into the substrate. 56 Collectively, these findings demonstrate that the superior oxidation resistance of multiphase oxide layers compared to single-phase counterparts is well established in high-temperature ceramics such as ZrB2 and SiC, as well as in metal matrix composites.
Based on the above microstructural analysis and literature corroboration, the TW60ZB ceramic exhibits outstanding oxidation resistance. In the TW60ZB composition, residual ZrB2 reacts with oxygen to form ZrO2 and B2O3. ZrO2, which possesses excellent high-temperature chemical stability, serves as an oxygen barrier layer that reduces oxygen diffusion into the substrate, further enhancing the overall oxidation resistance of the material. 57 Considering the oxide layer thickness and activation energy analysis, TW60ZB demonstrates the best oxidation resistance among all tested samples.
Figure 10 shows the cross-sectional morphologies of the oxide layer in (Zr,Ti,W)C-based multiphase ceramics after 15-h oxidation at 1000°C and 1100°C. The oxidised interface exhibits a distinct layered structure, sequentially comprising the outer oxide layer, inner oxide layer and substrate, with a continuous oxidation gradient formed between them. The outer oxide layer contains numerous voids and cracks, while the inner oxide layer is relatively dense. The higher-density inner layer effectively hinders inward oxygen diffusion, thereby providing excellent protection to the substrate.

Scanning electron microscope (SEM) images of oxidised cross-section of TW40ZB, TW50ZB and TW60ZB after oxidation at 1000°C and 1100°C for 15 h.
Quantitative analysis of the oxide layer thickness reveals that after oxidation at 1000°C, the TW60ZB sample exhibits an inner oxide layer thickness of 112.59 μm and a total oxide layer depth of 280.9 μm. When the temperature increases to 1100°C, the inner oxide layer thickness rises to 159.47 μm, with the total oxide layer depth reaching 440.61 μm. For comparison, existing studies show that pure ZrC can achieve an oxide layer depth of 550 μm after 1-h oxidation in air at 900°C, 9 while ZrC-TiC composites reduce this to 200 μm under identical conditions. 31 Remarkably, our study obtains shallower oxide layer depths under more severe conditions (higher temperature at 1100°C and longer duration of 15 h), a phenomenon analogous to that observed in ZrC + SiC + ZrB2 composites after oxidation. 58 The experimental results demonstrate that ZrB2 addition and the formation of (Zr,Ti)O2 solid solution significantly enhance the material's oxidation resistance.
Numerous pores and microcracks are observed in the outer oxide layers of TW50ZB and TW60ZB samples. Their formation may be attributed to the high-temperature oxidation of volatile WO3 and B2O3 generated from W x B y compounds, with these gaseous products escaping outward from the oxide layer interior, creating diffusion channels. Particularly in TW60ZB samples after 15-h oxidation at 1100°C, cracks between the oxide layer and substrate are visible, likely caused by CTE mismatch between them. During cooling, if the average linear CTE of the oxide layer (e.g., αZrC = 6.7 × 10−6 K−1, 11 αZrO2=10.5 × 10−6 K−1,59,60), thermal mismatch stress will induce tensile stress on the outer surface while the substrate experiences compressive stress, 59 leading to interfacial separation.
To visually compare the oxidation resistance of different samples, Figure 11 presents a line chart of oxide layer depth versus ZrB2 content. The results show that with increasing ZrB2 content, the inner oxide layer thickness gradually increases while the total oxide layer thickness progressively decreases. Among them, TW60ZB samples exhibit the thickest inner oxide layer and thinnest total oxide layer, demonstrating superior oxidation resistance.

Total oxidation layer depth and oxidation inner layer thickness of (Zr,Ti,W)C-based multiphase ceramics with varying ZrB2 content after 15-h oxidation at 1000°C and 1100°C.
The effect of ZrB2 content on the oxidation behaviour of the composites exhibits a distinct temperature dependence, which can be attributed to the varying roles of oxidation products in densification and volatilisation at different temperatures.
At 1000°C, the B2O3 glassy phase derived from ZrB2 oxidation possesses relatively high viscosity and low vapour pressure, enabling it to effectively fill pores and microcracks within the oxide layer. Consequently, the inner layer undergoes significant densification, which substantially impedes inward oxygen diffusion. As a result, the total oxide layer thickness decreases markedly with increasing ZrB2 content. 61 In this regime, the increase in inner layer thickness is negatively correlated with the reduction in total layer thickness,61,62 indicating that inner layer densification serves as the dominant mechanism governing oxidation resistance.
At 1100°C, the oxidation mechanism shifts fundamentally due to enhanced volatilisation of key oxidation products. On one hand, WO3 generated from the oxidation of tungsten begins to volatilise above 800°C–900°C 63 and undergoes severe volatilisation at 1100°C, leaving numerous gas escape channels in the outer layer. On the other hand, the viscosity of B2O3 decreases with increasing temperature, accompanied by an accelerated volatilisation rate, which compromises its ability to continuously fill voids, particularly in the outer layer. These combined effects result in a loose and porous outer layer structure, allowing rapid oxygen diffusion toward the inner layer. Although the inner layer thickens with increasing ZrB2 content, it fails to suppress the growth of the total oxide layer thickness. 64
Collectively, the densification-dominated mechanism that governs oxidation resistance at 1000°C is rendered ineffective at 1100°C due to the pronounced volatilisation of WO3 and B2O3. This transition from ‘inner-layer densification dominance’ to ‘outer-layer volatilisation dominance’ accounts for the temperature-dependent influence of ZrB2 content on the oxidation behaviour of the composites.
Oxidation mechanism
Based on the above analysis, the oxidation mechanism schematic of TW60ZB after 15-h oxidation at 1000°C and 1100°C is illustrated in Figure 12. The diagram distinguishes different oxide layer regions using yellow and dark grey colours, visually presenting the oxidation process and product distribution.

Schematic diagram showing oxidation mechanism of TW60ZB oxidised at 1000°C and 1100°C. (a) and (d) show the unoxidised state of the samples at 1000°C and 1100°C; (b) and (e) The samples during oxidation at 1000°C and 1100°C; (c) and (f) The samples after oxidation at 1000°C and 1100°C.
Figure 12 presents schematic diagrams illustrating the oxidation mechanism of TW60ZB at 1000°C and 1100°C. Figure 12(a) shows the initial state of TW60ZB before oxidation at 1000°C (unoxidised). Figure 12(b) illustrates the oxidation process of TW60ZB during 15 h of oxidation at 1000°C, depicting the gradual formation of the oxide layer, the filling of pores by the B2O3 glass phase and the formation of diffusion channels resulting from WO3 volatilisation. Figure 12(c) presents the final oxidation mechanism of TW60ZB after 15 h of oxidation at 1000°C, revealing a complete multi-layered oxide structure: a porous volatile outer layer, a dense inner layer formed by B2O3 consolidation and the formation of (Zr,Ti)O2 and (Ti,Zr)O2 solid solution layers at the interface. Figure 12(d) shows the initial state of TW60ZB before oxidation at 1100°C (unoxidised). Figure 12(e) illustrates the oxidation process of TW60ZB during 15 h of oxidation at 1100°C, highlighting significantly enhanced volatilisation of WO3 and B2O3 increased porosity and cracking in the outer layer, and an increased thickness of the dense inner layer accompanied by weakened protective capability. Figure 12(f) depicts the final oxidation mechanism of TW60ZB after 15 h of oxidation at 1100°C, showing the resulting oxide layer structure characterised by severe porosity in the outer layer and a relatively dense inner layer still containing volatile escape channels.
Under 1000°C oxidation conditions, the outer oxide layer primarily consists of monoclinic m-(Zr,Ti)O2, (Ti,Zr)O2 solid solution, (Ti,W)O x and incompletely volatilised WO3. Within the m-(Zr,Ti)O2 and (Ti,Zr)O2 composite oxide layer, ZrO2 and TiO2 particles are dispersed. These particles likely originate from WO3 and B2O3 being transported to the surface through physical adsorption or micro-area chemical reactions during high-temperature outward volatilisation. The inner oxide layer contains significant amounts of unvolatilised WO3 and B2O3 in addition to the aforementioned oxides. B2O3 fills pores in ZrO2 and oxidation microcracks, thereby enhancing the density of the inner layer. The dense inner oxide layer effectively inhibits outward carbon diffusion and inward oxygen penetration, significantly improving the material's overall oxidation resistance. 65
When the oxidation temperature increases to 1100°C, substantial thickening of the entire oxide layer is observed due to intensified oxidation reactions. At this temperature, WO3 becomes undetectable in the outer layer owing to aggravated high-temperature volatilisation, while TiO2 and ZrO2 contents significantly increase, attributable to migration effects induced by WO3 and B2O3 volatilisation. Furthermore, the (Ti,W)O x phase is no longer observed in either the inner or outer oxide layers, indicating complete oxidation of W elements into WO3 through vigorous volatilisation at elevated temperatures. Under these conditions, rapid evaporation of B2O3 reduces the effectiveness of the oxidation barrier.20,66
Conclusions
This study systematically investigated the high-temperature oxidation behaviour of (Zr,Ti,W)C-Me x B y multiphase ceramics prepared by SPS with varying ZrB2 additions, in static air at 1000°C–1100°C (the primary operating temperature range for divertor components in nuclear fusion devices). The oxide layer primarily consists of rutile-type t-TiO2, monoclinic m-ZrO2 and (Ti,Zr)O2 and (Zr,Ti)O2 solid solutions. Volatilisation of W and B elements in the forms of WO3 and B2O3, respectively, leads to porous structure formation in the outer oxide layer, while B2O3 effectively fills pores and microcracks in the inner layer, significantly enhancing density and thereby inhibiting inward oxygen diffusion and outward carbon migration. This multiphase oxide structure substantially reduces oxygen diffusion rate by impeding oxygen ion migration and pinning grain boundaries, thereby strengthening protective effects. At 1000°C, ZrB2-containing ceramics follow parabolic oxidation kinetics, with the rate constant Kp decreasing markedly with increasing ZrB2 addition. At 1100°C, vigorous WO3 and B2O3 volatilisation render the oxide layer porous, transitioning oxidation behaviour to linear kinetics. ZrB2 addition elevates oxidation activation energy from 13.3641 kJ/mol for (Ti,W)C to 47.4178 kJ/mol. When ZrB2 content reaches 60 mol% (TW60ZB), the material exhibits optimal oxidation resistance, with 355% improvement over (Ti,W)C. This research provides theoretical foundations and experimental support of (Zr,Ti,W)C-Me x B y multiphase ceramics for the design and application in components such as nuclear fusion device divertors or oxidation-resistant structural parts.
Footnotes
Acknowledgements
This work was financially supported by the National Natural Science Foundation of China (Nos. 52002098, 52002003 and 52402074), the Natural Science Foundation of Heilongjiang Province (LH2023E080) and the Natural Science Foundation of Hubei Province (Nos. 2025AFD102 and 2025AFD040).
Funding
The authors disclosed receipt of the following financial support for the research, authorship and/or publication of this article: This work was supported by the Natural Science Foundation of Hubei Province, National Natural Science Foundation of China, Natural Science Foundation of Heilongjiang Province (grant numbers 2025AFD040, 2025AFD102, 52002003, 52002098, 52402074, LH2023E080).
Declaration of conflicting interests
The authors declared no potential conflicts of interest with respect to the research, authorship and/or publication of this article.
