Microstructure and Mechanical Properties of High-Specific-Strength (TiVCrZr)100−xWx (x = 5, 10, 15 and 20) Refractory High-Entropy Alloys

With the increasing demand for high-specific-strength materials for high-temperature applications, particularly in the aerospace field, novel (TiVCrZr)100−xWx (x = 5, 10, 15 and 20) refractory high-entropy alloys (RHEAs) were developed. The phase formation, microstructure, and mechanical properties were studied. The (TiVCrZr)100−xWx RHEAs exhibit a relatively high specific strength and low density compared with the W-containing RHEAs and most of the W-free RHEAs. In (TiVCrZr)100−xWx RHEAs, Laves, BCC and Ti-rich phases are formed, where the Laves phase is the major phase, and the volume fraction of the BCC phase increases with increasing W content. (TiVCrZr)100−xWx RHEAs exhibit dendrite structures, where W is enriched in the dendrite region, and increasing W-rich precipitations corresponding to the BCC phase are observed. The improvement of the strength and hardness of RHEAs is mainly attributed to the evolution of the microstructure and corresponding strengthening effect of W. The empirical parameters and calculated phase diagram were investigated, which further explain and verify the formation and variation of phases. The present findings give more insights into the formation of multi phases in (TiVCrZr)100−xWx RHEAs, and explore their application potential in the aerospace industry and nuclear reactors due to their high specific strength and low-activation constituent elements.


Introduction
High-entropy alloys (HEAs), firstly proposed in 2004, have attracted extensive attention due to their outstanding mechanical properties, corrosion resistance, irradiation resistance and thermal stability [1][2][3][4][5]. HEAs are generally composed of at least five kinds of elements with an atomic percentage between 5% and 35% in equal mole or near equal mole ratios [6,7]. Differing from conventional alloys, HEAs tend to form simple solid solution structures, such as face-centered cubic (FCC), body-centered cubic (BCC) and close-packed hexagonal (HCP), rather than intermetallic compounds or amorphous phases due to their higher mix entropy [8][9][10][11]. However, the existing HEAs also contain several intermetallic compounds, which can strengthen the mechanical properties [12,13]. With the urgent demand for structural materials for application in high-temperature conditions, refractory high-entropy alloys (RHEAs), generally composed of refractory elements, were proposed in 2010 [14]. Senkov et al. [15] studied the microstructure and mechanical properties of WNbMoTa and WNbMoTaV RHEAs; the results showed that these two RHEAs have a typical BCC structure, and exhibit a higher strength and hardness especially at elevated temperature. RHEAs are therefore considered as promising next-generation high-temperature materials, which show great application potential in the nuclear fusion and aerospace fields [16,17].
High-specific-strength materials play a critical role in engineering applications, especially for the aerospace and automobile industries, as they can reduce the weight of

Materials and Methods
(TiVCrZr) 100−x W x (x = 5, 10, 15 and 20, denoted as W5, W10, W15 and W20 respectively) RHEAs were prepared by vacuum arc melting of the raw metals with a purity higher than 99.95% under a Ti-getter argon atmosphere. Each alloy ingot was remelted 12 times to ensure its homogeneity. The phase formation of RHEAs was characterized by a PANalytical X'Pert PRO MPD X-ray Diffractometer. The microstructures of RHEAs were inspected by a field emission scanning electron microscope (FESEM, Hitachi SU8020, Japan), where samples were polished by abrasive papers (Matador, Germany) with different grits from 400 to 5000. The attached energy dispersive spectrometer (EDS) was employed to inspect the elemental content and distribution of RHEAs. The Vickers hardness values of the RHEAs were measured by a Vickers hardness tester (WeiYee, HV-1000SA, China) under 1000 gf loads lasting 10 s, and the hardness values of 10 different positions for each sample were collected to obtain the average value. The samples were cut from a button-shaped ingot by a wire electrical discharge machine. The sample surfaces were polished by abrasive papers (Matador, Germany) with different grits from 400 to 2000. The typical W10 RHEA was annealed at 1000 • C for 10 h under argon atmosphere, where the heating rate was 10 • C/min. Compression tests of RHEAs were conducted by an electronic universal material testing machine (WanCe, DNS100, China) at an initial strain rate of 1 × 10 −3 s −1 at room temperature, where the samples' dimensions were 3 mm × 3 mm × 6 mm [28]. The fracture morphology of samples was also observed by the FESEM.

Phase Characterization of (TiVCrZr) 100−x W x RHEAs
According to the X-ray diffraction (XRD) patterns (Figure 1), (TiVCrZr) 100−x W x RHEAs have three different phases including the Laves (cF24-Cu 2 Mg, SG No. 194) [29], BCC and Ti-rich phases. Particularly, the W5 RHEA only has the Laves and Ti-rich phases. It is emphasized that the main peaks of the BCC and Laves phase are close to coincidence. The Laves phase is the major phase based on the diffraction intensity of the peak. It is worth noting that there were two peaks in the main peak of the Laves phase (Figure 1b), and with the increase of W content, the intensity of the left peak gradually increased, while the intensity of the right peak gradually decreased. For W15 and W20 RHEAs, the right peak was almost invisible. It is speculated that there are two Laves phases with slightly different compositions, and a Laves phase gradually disappears while another Laves phase gradually increases with increasing W content. A detailed discussion is presented in a later section. Moreover, the BCC increased based on the intensity of the corresponding diffraction peak. For the W20 RHEA, the main peak shifted obviously right, which may be influenced by the formation of more BCC phase.
Entropy 2022, 24, x FOR PEER REVIEW 3 of 19 of 1 × 10 −3 s −1 at room temperature, where the samples' dimensions were 3 mm × 3 mm × 6 mm [28]. The fracture morphology of samples was also observed by the FESEM.

Phase Characterization of (TiVCrZr)100−xWx RHEAs
According to the X-ray diffraction (XRD) patterns ( Figure 1), (TiVCrZr)100−xWx RHEAs have three different phases including the Laves (cF24-Cu2Mg, SG No. 194) [29], BCC and Ti-rich phases. Particularly, the W5 RHEA only has the Laves and Ti-rich phases. It is emphasized that the main peaks of the BCC and Laves phase are close to coincidence. The Laves phase is the major phase based on the diffraction intensity of the peak. It is worth noting that there were two peaks in the main peak of the Laves phase (Figure 1b), and with the increase of W content, the intensity of the left peak gradually increased, while the intensity of the right peak gradually decreased. For W15 and W20 RHEAs, the right peak was almost invisible. It is speculated that there are two Laves phases with slightly different compositions, and a Laves phase gradually disappears while another Laves phase gradually increases with increasing W content. A detailed discussion is presented in a later section. Moreover, the BCC increased based on the intensity of the corresponding diffraction peak. For the W20 RHEA, the main peak shifted obviously right, which may be influenced by the formation of more BCC phase.

Microstructure of (TiVCrZr)100−xWx RHEAs
The microstructure of (TiVCrZr)100−xWx RHEAs is shown in Figure 2. (TiVCrZr)100−xWx RHEAs exhibit a typical dendrite structure, which widely exists in HEAs due to the element segregation [30,31]. For the W5 RHEA, the microstructure comprised three distinct constituents, having bright, gray and dark regions. As shown in Figure 2a, the gray region was mainly distributed around the bright region, for example, the marked region between the blue and yellow dotted line. The fine gray region was dispersed in the interdendrite region. With the increase of W, the gray region around the bright region (W5) turned to a bright region, corresponding to the disappearance of the gray region around the bright region (W10, W15 and W20). However, a small amount of fine gray region dispersed in interdendrite regions was still observed. It is worth noting that the white region appeared inside the dendrite in W10, W15 and W20 RHEAs, and the proportion of the white region increased with the addition of W, corresponding to the appearance and enhancement of the BCC phase diffraction peak. Based on the XRD patterns and SEM images, the bright region, corresponding to the Laves phase, occupied most area of the microstructure, and

Microstructure of (TiVCrZr) 100−x W x RHEAs
The microstructure of (TiVCrZr) 100−x W x RHEAs is shown in Figure 2. (TiVCrZr) 100−x W x RHEAs exhibit a typical dendrite structure, which widely exists in HEAs due to the element segregation [30,31]. For the W5 RHEA, the microstructure comprised three distinct constituents, having bright, gray and dark regions. As shown in Figure 2a, the gray region was mainly distributed around the bright region, for example, the marked region between the blue and yellow dotted line. The fine gray region was dispersed in the interdendrite region. With the increase of W, the gray region around the bright region (W5) turned to a bright region, corresponding to the disappearance of the gray region around the bright region (W10, W15 and W20). However, a small amount of fine gray region dispersed in interdendrite regions was still observed. It is worth noting that the white region appeared inside the dendrite in W10, W15 and W20 RHEAs, and the proportion of the white region increased with the addition of W, corresponding to the appearance and enhancement of the BCC phase diffraction peak. Based on the XRD patterns and SEM images, the bright region, corresponding to the Laves phase, occupied most area of the microstructure, and the white region corresponds to the BCC phase, while the dark region corresponds to the Ti-rich phase. The gray region will be analyzed in detail in the subsequent section. the white region corresponds to the BCC phase, while the dark region corresponds to the Ti-rich phase. The gray region will be analyzed in detail in the subsequent section. EDS analysis of different regions was conducted to study the elemental distribution and further verify the types of each phase. As shown in Figure 3, V was distributed uniformly in the whole region, while other elements occur to different degrees of segregation due to their different melting points [32]. W was enriched in the dendrite region, particularly the white region, which is similar to previous studies [26,27,33]. This may result from its high melting point and limited solubility. The difference of W concentration between the bright and gray region drove the diffusion of W from the bright to gray region, resulting in the disappearance of the gray region around the bright region. With the increase of W content, W precipitated when the concentration of W in the dendrite region reached a certain limit, corresponding to the white region in Figure 2. In the solidification process of alloys, high melting point elements like W, Ta, Mo and Nb favor nucleate preferentially, resulting in the enrichment of elements in distinct regions [16,28,32]. However, Ti was mainly enriched in the interdendrite region. This may be attributed to the first nucleation of high melting point elements, leading to the latter solidified liquid with more Ti content, and further formation of a Ti-rich interdendrite region. Moreover, the interdendrite region was depleted of Cr, while the white region was depleted of Zr. EDS analysis of different regions was conducted to study the elemental distribution and further verify the types of each phase. As shown in Figure 3, V was distributed uniformly in the whole region, while other elements occur to different degrees of segregation due to their different melting points [32]. W was enriched in the dendrite region, particularly the white region, which is similar to previous studies [26,27,33]. This may result from its high melting point and limited solubility. The difference of W concentration between the bright and gray region drove the diffusion of W from the bright to gray region, resulting in the disappearance of the gray region around the bright region. With the increase of W content, W precipitated when the concentration of W in the dendrite region reached a certain limit, corresponding to the white region in Figure 2. In the solidification process of alloys, high melting point elements like W, Ta, Mo and Nb favor nucleate preferentially, resulting in the enrichment of elements in distinct regions [16,28,32]. However, Ti was mainly enriched in the interdendrite region. This may be attributed to the first nucleation of high melting point elements, leading to the latter solidified liquid with more Ti content, and further formation of a Ti-rich interdendrite region. Moreover, the interdendrite region was depleted of Cr, while the white region was depleted of Zr.
In order to obtain the specific element content of distinct regions, EDS point tests were carried out, and three to five points of each region were inspected to obtain their average values. The results are shown in Table 1. As mentioned above, the phase corresponding to the gray region still needed to be analyzed. According to the elemental content of the gray region, it can be found that the content of V, Cr, Zr satisfied the formula (relation) of Zr(V, Cr) 2.04 , which is almost consistent with the form of pseudo-binary system Zr(V 1−x Cr x ) 2 Laves intermetallic compound [34,35]. In general, the Laves phase, as a topologically close-packed (TCP) phase, is an intermetallic compound with a chemical formula of the main AB 2 type, where the ratio of the atomic radius of elements A and B ranges from 1.09 to 1.55 [29]. For instance, ZrV 2 and ZrCr 2 Laves phases existing in the HEAs have been reported [36][37][38]. The appearance of the pseudo-binary system Zr(V 1−x Cr x ) 2 Laves phase may be attributed to the mutual substitution (replacement) of Cr and V. As mentioned earlier, the bright region also corresponded to the Laves phase (cF24-Cu 2 Mg, SG No. 194). The EDS results of the bright region showed that the content of V, Cr, Zr and W satisfied the formula (relation) of Zr(V, Cr, W) 2.05 , a quaternary Laves, which may result from the replacement of V or Cr by W. Additionally, Senkov et al. studied the microstruc- ture and phase formation of CrNbTiVZr RHEA, and the experimental results showed that the Laves phase is also not entirely the binary ZrCr 2 Laves phase, but also contains other compositional elements. The actual composition of the Laves phase fulfills the formula (Cr, Nb, V) 2 (Ti, Zr), which indicates that part of Cr was substituted by V and Nb, and part of Zr was substituted with Ti. As a result, both the gray and bright regions belong to the Laves phase, and the difference between them was the W content. As shown in Figure 1b (B-spline curve fitting), there were two diffraction peaks in the main peak of the Laves phase, corresponding to two Laves phases with slightly different compositions. With the increase of W, the gray region disappeared due to the diffusion of W, and accordingly the right diffraction peak also disappeared, as shown in Figure 1b   In order to obtain the specific element content of distinct regions, EDS point tests were carried out, and three to five points of each region were inspected to obtain their average values. The results are shown in Table 1. As mentioned above, the phase corresponding to the gray region still needed to be analyzed. According to the elemental content of the gray region, it can be found that the content of V, Cr, Zr satisfied the formula (relation) of Zr(V, Cr)2.04, which is almost consistent with the form of pseudo-binary system Zr(V1−xCrx)2 Laves intermetallic compound [34,35]. In general, the Laves phase, as a topologically close-packed (TCP) phase, is an intermetallic compound with a chemical formula of the main AB2 type, where the ratio of the atomic radius of elements A and B ranges from 1.09 to 1.55 [29]. For instance, ZrV2 and ZrCr2 Laves phases existing in the Moreover, the difference of atomic sizes of V, Cr and W was less than 15%, and V, Cr and W had similar BCC crystalline structures, which accorded with the replacement condition. Additionally, considering the enthalpy of mixing between pairs of atoms, the calculated results demonstrate that the enthalpy of mixing of Zr-W (−9 kJ/mol) ranked only second to that of Zr-Cr (−12 kJ/mol) in alloys, as shown in Figure 4. Large negative enthalpy of mixing usually represents strong attraction, i.e., combination tendency between atoms [39]. It is possible that the attraction of Zr to W leads to the substitution of Cr and V by W due to the relatively large negative enthalpy of mixing of Zr-W. However, the white precipitation was enriched with W while depleted of Zr, which may be attributed to the large difference in melting point.  Overall, the content of W in the dendrite region gradually increased, while the elemental content of interdendrite almost remained the same. The variation of microstructure in the dendrite region has a critical effect on the mechanical properties of (TiVCrZr)100−xWx RHEAs.

Mechanical Properties of (TiVCrZr)100−xWx RHEAs
The compressive stress-strain curves of (TiVCrZr)100−xWx RHEAs at room temperature are shown in Figure 5. The corresponding fracture strength (σb) and strain (εf) are listed in Table 2. With the increase of W content, the fracture strength (σb) and strain (εf) were gradually enhanced except for the W20 RHEA. The strength of the alloy system increased from an initial 1574 MPa (W5) to a maximum of 1894 MPa (W15), while the strength of W20 (1754 MPa) decreased with further increase of W content. The high strength and low ductility of the alloy system are mainly attributed to the formation of the Laves phase, which exhibited a high strength and brittleness at room temperature. The fitting results (Figure 5c) of hardness and x value indicate that the hardness of the alloy system increased nearly linearly with the increase of x value, where the goodness of fit R 2 = 0.958 was close to 1. This implies a highly linear relationship between the hardness and x value [40]. With the increase of W, the W content of the dendrite region gradually in- Overall, the content of W in the dendrite region gradually increased, while the elemental content of interdendrite almost remained the same. The variation of microstructure in the dendrite region has a critical effect on the mechanical properties of (TiVCrZr) 100−x W x RHEAs.

Mechanical Properties of (TiVCrZr) 100−x W x RHEAs
The compressive stress-strain curves of (TiVCrZr) 100−x W x RHEAs at room temperature are shown in Figure 5. The corresponding fracture strength (σ b ) and strain (ε f ) are listed in Table 2. With the increase of W content, the fracture strength (σ b ) and strain (ε f ) were gradually enhanced except for the W20 RHEA. The strength of the alloy system increased from an initial 1574 MPa (W5) to a maximum of 1894 MPa (W15), while the strength of W20 (1754 MPa) decreased with further increase of W content. The high strength and low ductility of the alloy system are mainly attributed to the formation of the Laves phase, which exhibited a high strength and brittleness at room temperature. The fitting results (Figure 5c) of hardness and x value indicate that the hardness of the alloy system increased nearly linearly with the increase of x value, where the goodness of fit R 2 = 0.958 was close to 1. This implies a highly linear relationship between the hardness and x value [40]. With the increase of W, the W content of the dendrite region gradually increased, particularly with increased W content of the Laves phase, resulting in the increase of hardness. For the W20 RHEA, although a large amount of W-rich phase led to a decrease of the Laves phase, the grain size was reduced and the high hardness of the W-rich phase still improved the hardness of the alloys. The addition of W can therefore significantly improve the strength and hardness of alloys over a certain range.  The above-mentioned mechanical properties results show that (TiVCrZr)100−xWx RHEAs exhibit high strength and brittleness due to the formation of the Laves phase. As shown in Figure 6, (TiVCrZr)100−xWx RHEAs exhibit a typical brittle fracture morphology, which is similar to the fracture morphology observed in CrNbTiVZr RHEA [24]. Some particles were observed on the fracture surface, particularly for W5 and W10 RHEAs, which may be attributed to the catastrophic fracture of the brittle Laves phase. It is worth noting that some spherical particles and a cluster of spherical particles were also observed. This fracture characteristic of alloys resulted from the instantaneous release of elastic energy, which led to the melting of the local region. According to the stress-strain curves of (TiVCrZr)100−xWx RHEAs, there were a large amount of accumulated elastic energy in specimens during the compressive process, and light emission appeared at the moment of fracture.  The above-mentioned mechanical properties results show that (TiVCrZr) 100−x W x RHEAs exhibit high strength and brittleness due to the formation of the Laves phase. As shown in Figure 6, (TiVCrZr) 100−x W x RHEAs exhibit a typical brittle fracture morphology, which is similar to the fracture morphology observed in CrNbTiVZr RHEA [24]. Some particles were observed on the fracture surface, particularly for W5 and W10 RHEAs, which may be attributed to the catastrophic fracture of the brittle Laves phase. It is worth noting that some spherical particles and a cluster of spherical particles were also observed. This fracture characteristic of alloys resulted from the instantaneous release of elastic energy, which led to the melting of the local region. According to the stress-strain curves of (TiVCrZr) 100−x W x RHEAs, there were a large amount of accumulated elastic energy in specimens during the compressive process, and light emission appeared at the moment of fracture. Moreover, after annealing at 1000 °C for 10 h, the W10 RHEA slightly demonstrated increased hardness and strength compared with the as-prepared specimens. As shown in Table 3, the hardness increased from 639 ± 23 HV to 647 ± 13 HV, while the strength increased from 1867 ± 52.5 MPa to 1930 ± 45 MPa. This indicates that the (TiVCrZr)100−xWx RHEAs may have better mechanical properties at high temperature. As shown in Figure  7a, the microstructure of the annealed specimens still exhibited four different contrasts including white, bright, gray and dark regions, which is similar to the microstructure of as-prepared specimens (Figure 7b). Thus, the (TiVCrZr)100−xWx RHEAs may possess higher thermostability and have potential for high temperature applications.  Moreover, after annealing at 1000 • C for 10 h, the W10 RHEA slightly demonstrated increased hardness and strength compared with the as-prepared specimens. As shown in Table 3, the hardness increased from 639 ± 23 HV to 647 ± 13 HV, while the strength increased from 1867 ± 52.5 MPa to 1930 ± 45 MPa. This indicates that the (TiVCrZr) 100−x W x RHEAs may have better mechanical properties at high temperature. As shown in Figure 7a, the microstructure of the annealed specimens still exhibited four different contrasts including white, bright, gray and dark regions, which is similar to the microstructure of as-prepared specimens (Figure 7b). Thus, the (TiVCrZr) 100−x W x RHEAs may possess higher thermostability and have potential for high temperature applications.

Phase Structure
The experimental results indicate that (TiVCrZr)100−xWx RHEAs exhibit a complex phase formation rather than the simple solid solution as some HEAs exhibit. Empirical parameters and calculated phase diagram were employed for further understanding of the phase formation. Generally, HEAs tend to form solid solution phases due to the high entropy of mixing, while partial existing HEAs also tend to form Laves, σ, μ and χ phases of intermetallic compounds [41]. Particularly, the Laves phase is usually formed in Crcontaining and Zr-containing HEAs [23,42,43]. The phase type has a critical effect on the properties of HEAs, which are related to the thermodynamic, electronic parameters and atomic size of the alloy system [29,44]. The corresponding parameter values of compositional elements are listed in Table 4. Empirical parameters including the entropy of mixing (ΔSmix), enthalpy of mixing (ΔHmix) and atomic size difference (δ) were proposed, which can be calculated by the following formulas [46]: where R is the gas constant; n is the number of elements; ci and cj are the atomic percentage of the ith and jth elements; Δ ij mix is the enthalpy of mixing between the ith and jth elements, which can be obtained by Ref. [47]; ri is the atomic radius of the ith element; and

Phase Structure
The experimental results indicate that (TiVCrZr) 100−x W x RHEAs exhibit a complex phase formation rather than the simple solid solution as some HEAs exhibit. Empirical parameters and calculated phase diagram were employed for further understanding of the phase formation. Generally, HEAs tend to form solid solution phases due to the high entropy of mixing, while partial existing HEAs also tend to form Laves, σ, µ and χ phases of intermetallic compounds [41]. Particularly, the Laves phase is usually formed in Cr-containing and Zr-containing HEAs [23,42,43]. The phase type has a critical effect on the properties of HEAs, which are related to the thermodynamic, electronic parameters and atomic size of the alloy system [29,44]. The corresponding parameter values of compositional elements are listed in Table 4. Empirical parameters including the entropy of mixing (∆S mix ), enthalpy of mixing (∆H mix ) and atomic size difference (δ) were proposed, which can be calculated by the following formulas [46]: where R is the gas constant; n is the number of elements; c i and c j are the atomic percentage of the ith and jth elements; ∆H mix ij is the enthalpy of mixing between the ith and jth elements, which can be obtained by Ref. [47]; r i is the atomic radius of the ith element; and r a = ∑ n i=1 c i r i is the average atomic radius. Then, Yang et al. [48] optimized the aforementioned rules and proposed a novel parameter Ω, which can be calculated by the following formulas: where T m is the average melting point of the alloy and (T m ) i is the melting point of the ith element. The result shows that the formed phase types can be divided into four categories including solid solutions; intermetallics; bulk metallic glasses; and solid solutions and intermetallics in terms of the relationship of Ω and δ. The corresponding calculated values of (TiVCrZr) 100−x W x RHEAs are distributed in the region of formation of solid solutions and intermetallics. However, these parameters cannot predict the specific type of phases. Guo et al. [49] studied the effect of valence electron concentration (VEC) on the phase type of HEAs, which can be calculated by the following formula: where (VEC) i is the valence electron concentration of the ith element, referring to Ref. [45].
The results indicate that HEAs tend to form a BCC phase or FCC phase with VEC< 6.87 or VEC ≥ 8 respectively. When 6.87 ≤ VEC < 8.0, BCC and FCC phase coexist in HEAs. Yurchenko et al. [50] studied the effect of Allen electronegativity difference (∆χ Allen ) and δ on the phase type of HEAs, which can be calculated by the following formula: where χ Allen i is the Allen electronegativity of the ith element and χ a = ∑ n i=1 c i χ Allen i is the average Allen electronegativity. The results show that the Laves phase tends to form in HEAs when δ > 5.0% and ∆χ Allen > 7.0% [50]. Therefore, the VEC, δ and ∆χ Allen values shown in Table 5 demonstrate that (TiVCrZr) 100−x W x RHEAs have the Laves and BCC phases, which is consistent with the experimental results. Beyond empirical criteria, the calculated phase diagram is also widely applied in the analysis of the phase formation of HEAs. In this work, thermo-Calc (v. 2019a, Thermo-Calc Software AB, Stockholm, Sweden) in combination with the HEAs thermodynamic database (v. 2.1.1) were used to simulate the equilibrium solidification process. The calculated phase diagram also indicates that (TiVCrZr) 100−x W x RHEAs exhibit a complex phase formation, including the Laves, BCC and HCP phases. Although the calculated phase diagram cannot entirely be consistent with the experimental results, some of the results still have reference value. As shown in Figure 8, (TiVCrZr) 100−x W x RHEAs exhibit the ordered intermetallic compound phase in the as-cast state despite the high entropy of mixing. Nevertheless, (TiVCrZr) 100−x W x RHEAs exhibit a mixture of single BCC solid solution and liquid phases at a higher temperature, which can be explained from the viewpoint of Gibbs free energy.
Gibbs free energy has a critical effect on the formation of phases, which can be expressed by the following formula [12]: where H is the enthalpy; T is the thermodynamic temperature; and S is the entropy. According to Formula (8), the effect of entropy on Gibbs free energy becomes more significant under high temperature conditions, and can effectively reduce the Gibbs free energy of the system [51]. Lower Gibbs free energy can stabilize the solid solution. For (TiVCrZr) 100−x W x RHEAs, the single BCC solid solution only exists in a narrow temperature range. However, the temperature range of a single BCC solid solution increases with increasing mixing entropy, resulting from the addition of W as shown in Table 5. This can be explained by the effect of the entropy increase on the reduction of Gibbs free energy. With the decrease of temperature, the Laves phase appears and compresses the single BCC solid solution range. The Laves phase becomes the major compositional phase until complete solidification. The equilibrium phase diagram shows that the fraction of the Laves phase decreases and the fraction of BCC phase increases with the increase of W content, which is consistent with the evolution of microstructure, i.e., the gradual precipitation of the white region corresponding to the BCC phase. It is worth noting that the Laves phase exists over a larger temperature range and the stable temperature range of the Laves phase increases with increasing W content. Additionally, the temperature at which the Laves phase starts to dissolve and the temperature at which the Laves phase completely dissolves also gradually increase. As mentioned earlier, W substitutes Cr and V to form a similar (Cr, V, W) 2 Zr-type Laves phase, resulting in an improved melting point of the Laves phase. The calculated phase diagram can predict the existence and variation trend of the Laves and BCC phases more clearly, while the HCP phase was not observed according to the XRD patterns. Nevertheless, the XRD results and microstructure show the existence of the Ti-rich phase. Indeed, Ti exhibits a HCP crystal structure at low temperature, and the equilibrium phase diagram simulates the formed phases at a slow cooling rate [32], which may result in the existence of a HCP phase as shown in Figure 8.

Specific Strength
Materials with a high specific strength have attracted extensive attention due to their huge advantages in engineering applications. (TiVCrZr)100−xWx RHEAs have high specific strength as shown in Table 6, exhibiting great application potential. In order to better reveal the behavior of high specific strength, the following contents will be discussed. On

Specific Strength
Materials with a high specific strength have attracted extensive attention due to their huge advantages in engineering applications. (TiVCrZr) 100−x W x RHEAs have high specific strength as shown in Table 6, exhibiting great application potential. In order to better reveal the behavior of high specific strength, the following contents will be discussed. On the one hand, (TiVCrZr) 100−x W x RHEAs comprise four low-density refractory elements, which make the alloy system have a lower density. The density of (TiVCrZr) 100−x W x RHEAs was calculated by the rule of mixture, which can be expressed by the following formula [16]: where c i is the atomic percentage of the ith element; n is the number of elements; A i is the atomic weight of the ith element; and ρ i is the density of the ith element. The theoretical density and actually measured density results are almost consistent with previous research [23,32], and thus the theoretical density was employed in this work. As shown in Formula (9), the low density of components in the alloy system can reduce the whole density of the alloys. On the other hand, W addition significantly improves the strength of the alloy system while obtaining low density. With the increase of W content, the W content in the dendrite region increases, and the initial gray region (Zr(V 1−x Cr x ) 2 Laves phase) around the bright region transforms into a bright region (Zr(V, Cr, W) 2 Laves phase). It is known that HEAs have a cocktail effect, which implies that HEAs can sufficiently utilize the comprehensive properties of compositional elements. W possesses high strength and hardness [52,53], which has a positive effect on the mechanical properties of the alloy system. As a result, the bright region composed of the Zr(V, Cr, W) 2 Laves phase has a higher strength and hardness than that of the gray region composed of the Zr(V 1−x Cr x ) 2 Laves phase. The enhancement of the major phase (Laves phase), benefiting from W addition, improves the hardness and strength of the alloy system. For the W20 RHEA, a larger number of the W-rich BCC phase (brittle phase) precipitates in the dendrite region, leading to increasing brittleness at room temperature and reduced volume fraction of the Laves phase. These result in the decrease of the exhibited strength of the W20 RHEA. Although W addition also increases the density of the alloy system, the strengthening effect of W on the strength exceeds that of W on the density of the alloy system over a certain range. Therefore, with the increase of W, the specific strength of the alloy system first increases and then decreases, where W10 RHEAs can reach a maximum specific strength of 256.65 MPa·cm 3 ·g −1 .
To further understand the specific strength of (TiVCrZr) 100−x W x RHEAs, some RHEAs were selected as comparative materials, as these RHEAs meet the following conditions. The involved RHEAs were prepared by vacuum arc melting. The selected RHEAs comprised W-containing RHEAs and W-free RHEAs, where W-free RHEAs consisted of nine refractory elements and contained any three of the four elements (Ti, V, Cr, Zr). As shown in Figure 9, (TiVCrZr) 100−x W x RHEAs exhibited a higher specific strength and lower density, particularly for W10 RHEAs. Overall, the specific strength of (TiVCrZr) 100−x W x RHEAs exceeded almost all W-containing RHEAs prepared by vacuum arc melting. The specific strength of the W10 RHEA was only inferior to the two W-free RHEAs among the selected RHEAs, i.e., Mo 0.5 VNbTiCr 1.5 and Mo 0.5 VNbTi Cr 2.0 . As a result, (TiVCrZr) 100−x W x RHEAs show great application potential in the aerospace industry due to their high specific strength. Additionally, the compositional elements of the alloys are low activation elements, and tungsten and its alloys possess great application prospects [40,54,55], implying that (TiVCrZr) 100−x W x RHEAs may have excellent irradiation resistance. (TiVCrZr) 100−x W x RHEAs not only have great application potential for turbine blades, but also outstanding application value in structural materials for nuclear reactors. Improving the plasticity of the (TiVCrZr) 100−x W x RHEAs could expand their engineering applications. However, improving the plasticity is challenging due to the strength-ductility trade-off. Many ways have been proposed to further improve the plasticity by maintaining the strength [28,56,57]; for example, Sun et al. [56] tuned elemental composition to minimize the Gibbs free energy, which effectively suppressed the formation of intermetallic compounds. This effectively improved plasticity. For the present (TiVCrZr) 100−x W x RHEAs, due to the larger negative mixing enthalpy between Zr and other compositional elements, it was more likely to result in the formation of the Zr-containing Laves phase. Therefore, reducing Zr content could be effective to suppress the formation of the Laves phase, further improving the plasticity of the (TiVCrZr) 100−x W x RHEAs.

Conclusions
In this work, (TiVCrZr)100−xWx (x = 5, 10, 15 and 20) RHEAs with a high specific strength were developed and investigated. The phase formation, microstructure and mechanical properties were examined. The main findings are summarized as follows: (1) (TiVCrZr)100−xWx RHEAs have the Laves, BCC and Ti-rich phases, where the Laves phase is the major phase, and the volume fraction of the BCC phase increases with increasing W content. Dendrite structures were observed. The EDS results show that W is enriched, and then gradually diffuses, and finally precipitates in the dendrite region, corresponding to the transformation from the gray region (Zr(Cr, V)2-type

Conclusions
In this work, (TiVCrZr) 100−x W x (x = 5, 10, 15 and 20) RHEAs with a high specific strength were developed and investigated. The phase formation, microstructure and mechanical properties were examined. The main findings are summarized as follows: (1) (TiVCrZr) 100−x W x RHEAs have the Laves, BCC and Ti-rich phases, where the Laves phase is the major phase, and the volume fraction of the BCC phase increases with increasing W content. Dendrite structures were observed. The EDS results show that W is enriched, and then gradually diffuses, and finally precipitates in the dendrite region, corresponding to the transformation from the gray region (Zr(Cr, V) 2 -type Laves phase) to bright region (Zr(Cr, V, W) 2 -type Laves phase) and the increased white region. (2) Empirical parameters and the calculated phase diagram (Thermo-Calc Software AB) are beneficial for understanding the mechanisms of phase formation and variation.
The results indicate that the calculated phase type and variation are consistent with the experimental results. Moreover, the calculated phase diagram confirmed the substitution of W for Cr and V to some extent. (3) W addition significantly enhances the strength and hardness of (TiVCrZr) 100−x W x RHEAs. The strength of the alloy system improved from 1574 MPa (W5) to 1894 MPa (W15), and the hardness improved from 588 HV to 694 HV. Generally, the higher specific strength benefits those alloys to be envisaged as structure materials. Compared with the selected RHEAs and Inconel 718, the specific strength of (TiVCrZr) 100−x W x RHEAs exceed that of almost all W-containing RHEAs, and the highest specific strength was found to be 256.65 MPa·cm 3 ·g −1 of W10 RHEA, which was just below that of the Mo 0.5 VNbTiCr 1.5 and Mo 0.5 VNbTiCr 2.0 W-free RHEAs. Moreover, the specific strength of (TiVCrZr) 100−x W x RHEAs also exceed that of Inconel 718 (about 148 MPa·cm 3 ·g −1 ), which is widely applied in the aerospace industry. The improvement of strength and hardness resulted from the strengthening effect of W on the Laves phase.