Design of high-strength refractory complex solid-solution alloys

Nickel-based superalloys and near-equiatomic high-entropy alloys containing Molybdenum are known for higher temperature strength and corrosion resistance. Yet, complex solid-solution alloys offer a huge design space to tune for optimal properties at slightly reduced entropy. For refractory Mo-W-Ta-Ti-Zr, we showcase KKR electronic-structure methods via the coherent-potential approximation to identify alloys over 5-dimensional design space with improved mechanical properties and necessary global (formation enthalpy) and local (short-range order) stability. Deformation is modeled with classical molecular dynamic simulations, validated from our first-principles data. We predict complex solid-solution alloys of improved stability with greatly enhanced modulus of elasticity ($3\times$ at 300 K) over near-equiatomic cases, as validated experimentally, and with higher moduli above 500~K over commercial alloys ($2.3\times$ at 2000 K). We also show that optimal complex solid-solution alloys are not described well by classical potentials due to critical electronic effects.


INTRODUCTION
Nickel-based superalloys exhibit high-temperature strength, toughness, and oxidation resistance in harsh environments. 1 Improving existing single-crystal alloys is unlikely as melting is near 1350 • C, and, heat treatment lowers this to ∼1270 • C. In high-speed turbines, melting reduces below 1250 • C at the zone between the bond coat (e.g., NiAl) and the singlecrystal blade. 2 As such, the engine efficiency and thrustto-weight ratio can be improved by a guided search for new materials. High-entropy alloys based on refractory elements may achieve higher temperature operation with superior creep strength. 3 Typical refractory high-entropy alloys exhibit a yield strength of 500-700 MPa at 1200 • C, surpassing Ni-based superalloys. 4 Indeed, at elevated temperatures Mo-based alloys show good thermal (higher conductivities with lower strains 5 ) and mechanical (machinability) 6 properties, making them promising candidates.
Almost all high-entropy alloys for which mechanical properties have been reported are based on Cr-Fe-Co-Ni with other elements added, e.g., Al, 7 Mn, 8 Mo, 9 and Ti. 10 CoCrFeNi exhibits very high compression strength at 300 K, often exceeding 1500 MPa. Strains in as-cast condition do not often exceed 5−7%, albeit a few exhibit 25−33%. 10,11 Annealing does improve ductility of as-cast alloys. 11 As with conventional alloys, a rapid decrease in strength (i.e., Young's modulus, E) occurs above 0.6 of the melting temperature T m , and the strength of alloys approaches 100 MPa at 1273 K. 7 High-entropy alloys consist of N (≥5) elements in near equiatomic compositions (c α ∼1/N ), giving maximal point (mixing) entropy (S pt =− N α=1 c α ln c α max − − → ln N ), that may better form solid solutions due to a compromise between the large S pt and a formation (mixing) energy ∆E f orm that is not too positive (strongly clustering) nor too negative (strongly ordering). 7 As for binary solid solutions, Hume-Rothery's rules 12 for atomic size difference (δ), crystal structure, valence electron concentration (VEC), and electronegativity difference (∆χ) play a similar role in high-entropy alloy formation. The production of several single or multi-phase alloys with face-centered cubic (FCC or A1), body-centered cubic (BCC or A2), hexagonal closepacked (HCP or A3), or cubic diamond (A4) structures exhibiting enhanced high-temperature strength, ductility, fracture and creep resistance to corrosion, [7][8][9][10]13,14 and thermal stability 15 validates the concept of HEAs. 7 MoWVNbTa, e.g., with a density of 12.2 g/cm 3 , has a reported usable strength up to 1873 K. 4 Nonetheless, from an alloy design perspective, complex solid-solution alloys (CSAs) offer a huge design space to tune properties, especially considering the strong effects alloying has on electronic properties ("band" filling, hybridization, Fermi-surface nesting, ...), phase stability, and structure. The CSAs comprised of whole composition (Gibbs) space, however, high-entropy alloys are a subset of it. Optimized CSAs offer a slightly reduced entropy with a singlephase region, or two-phase region for enhanced mechanical properties, existing in a desired operational temperature range. [16][17][18] Here we narrow the design of high-strength, refractory (Mo-W)-Ta-(Ti-Zr) alloys via KKR electronic-structure methods within density-functional theory (DFT) using the coherentpotential approximation (CPA) to handle chemical disorder and thermodynamic averaging. [16][17][18][19] The well-established KKR-CPA predicts structural properties [e.g., Young's (E) or bulk modulus (B)], and phase stability (∆E f orm vs. {c α }), as well as short-range order (SRO) via thermodynamic linear response, 18,20-23 a method, in particular, which revealed the origin for Hume-Rothery's size-effect rule. 21,24 Notably, global stability (∆E f orm ) and local instability (SRO) should be jointly assessed: While segregation is expected for ∆E f orm > 0, SRO can be segregating from local compositional instabilities even if ∆E f orm < 0. To predict mechanical behavior (e.g., E) versus temperature (T) rapidly, we performed extensive molecular dynamics (MD) simulations based on semi-empirical potentials, validated in part by first-principles results (and also highlighting limitations of such methods). The tuned and proposed refractory quinary alloys and their properties are placed in context to Hume-Rothery-type design targets and compared to experiments.

RESULTS AND DISCUSSION
A. Hume-Rothery Design Targets: High-entropy alloys contain elements with c α ∼35-12 at.% (N =3-8). Trial-and-error has led to alloys with simple crystal structures, and a few with extraordinary properties, 25 e.g., formability using size disparate elements for confusion by design. 26 For CSAs design of phase stability and of electronic and mechanical behavior, targets for DFTbased KKR-CPA are limited by extending Hume-Rothery 12 criteria: 1. Size: Solute and host atomic radii (in elemental solid) must differ by < 15%. 12,21 For CSAs, withr= N i=1 ciri, size limit in terms of standard deviation is sensible: A few comments are warranted. In #1, the 6-6.8% achieves 15%-rule for CSAs with ∼50% confidence level, an inequality also found empirically. 27,28 Extending #3 via electronic density of states concepts, A2 forms for 4<VEC<6 as stability increases when bonding d-states fill, and is maximal when half-filled (VEC≈6); antibonding states fill with VEC>6 (above a pseudogap, see results) and stability decreases. Indeed, A2 CSAs are observed when VEC is 5±1. 29 For 6.8 ≤ VEC ≤ 8 other phases compete, e.g., FeCr has VEC=7 (like Mn) and constituent's structure are both A2, yet the CSA is unstable to the σ-phase, as often appears. 27 Again, from band filling, A1 becomes more stable for VEC>8. 30 CSAs are indeed observed to form within these rules. 31 ∆E f orm lower limit in #5 is set by −T a S pt , where annealing temperature (needed for kinetics) is T a ∼0.55T m (∼1000-1650 K for refractories); upper limit is set such that miscibility gap T M G c <T a (where 158 K∼1 mRy). For δ > 5%, CSAs with ∆E f orm > 5 mRy form complex phases, but tend to form metastable metallic glasses for ∆E f orm < −11 mRy. 32 Considering binaries and supercells, these limits for CSA formation are supported. 33 As we have shown, transition temperatures from α (CSA) to β (ordered or segregated) phases are well estimated from calculated ∆E f orm 's: 34,35 For segregating CSAs (∆E f orm >0), T M G c ≈ ∆E f orm /S pt ; 34 and, for ordering CSAs with ∆E α→β f orm = ∆E α f orm − ∆E β f orm > 0, the orderdisorder transition is T od c ≈∆E α→β f orm . 35 Notably, #4 also reveals if vibrations are important, as vibrational entropy in binaries correlates as ∆S vib =−∆χ/3 (±0.06∆χ). 36  , which reproduces measured trends without phonon calculations. 34 These estimates are within 5 − 10%. 34,35 As a predictive guide, we use KKR-CPA results to tune (∆E f orm , δ, VEC, and ∆χ) versus {c α } to find (Mo-W)-Ta-(Ti-Zr) alloys in 5D space with better stability and mechanical properties. Results identify the stability of competing phases, possible multiphase regions, electronic properties, and practical design limits. We use the above criteria to restrict search space for mechanical simulations. Here, via the KKR-CPA, we search all CSAs without restrictions on {c α }, or the need for large supercells, as A1,A2 (A3) have only 1 (2) atoms per cell. In this quinary, atomic size of Zr (1.60Å) is largest, followed by Ti, Ta, and W, Mo (1.46, 1.43, and 1.37, 1.36Å), where bandwidths (inversely related to atomic size) and alloy hybridization determine the effect of size. 21 For χ (or ∆χ's), reflecting solubility and vibrational entropy, 34 First, we exemplify in Fig. 1, top-panel, our accuracy for ∆E f orm vs.
x in Ta 1−x W x , which agrees well with measured values (within 5%), and ordering enthalpies are low (<310 K) compared to melting. Also, we show results for specific ordered cells, which are compared to and agree well with other reliable bandstructure methods [e.g., VASP pseudo-potential 38 and fullpotential linear-augmented plane-wave (FLAPW) 39 ]. For Mo x (WTaTiZr) 1−x in Fig. 1 (bottom-panel), we find that A2 is favored over A1 or A3, and that increasing %Mo (larger x) helps stabilize A2. So, in this work, we focus on A2-phase of (Mo-W)-Ta-(Ti-Zr) alloy. For x 0.4 A4 phase competes with A2, and Frank-Kasper phases, like C15-Mo 2 (Ti-Zr), may be anticipated. Usually elements from group IIA/IVA of the periodic table, e.g., Al/Si, are added to stabilize or change, e.g., oxidation resistance. However, we find that adding Al stabilizes A2-phase up to 20%Mo, and similar behavior of Al-addition has been seen in other alloys too. 19 On the other hand, Si-addition comes out to be energetically less favorable than Al-addition. High-Throughput Assessments: For (Mo-W)-Ta-(Ti-Zr) results are most easily presented in a cut through 5D {c α } space to visualize with only two parameters (x, y) along lines or planes (Fig. 2), changing {c α } in obvious ways. For fast screening of (Mo-W)-Ta-(Ti-Zr) 5D composition design space, we used an estimate for lattice constant to perform "high-throughput" [40][41][42][43][44][45][46] calculations to discover the best alloys in terms of phase-stability and/or mechanical behavior. Specifically, we estimated alloy lattice constants via Vegard's rule, which is the concentration-weighted sum of volume optimized elemental lattice constants in the parent alloys (A2) phase, or, simply, a alloy = i c X i a X i , where i=1, 5 and X= Mo, W, Ta, Ti, Zr. The estimated lattice constants are within 1-3% with respect to the optimized lattice constants for all considered compositions. To downselect regions of interest, we perform the calculation over the entire design space and chose increments in {c α } every 5% to sweep whole 5D space, see Fig. 2. For selected alloys, we perform full lattice-optimization to determine ∆E f orm and B, and detail the electronic-structure (dispersion and density of states), and the thermodynamic short-range order (incipient ordering) for alloy design. Assessment & Validation: We now assess CSAs that best satisfy design criteria, and local stability. Along with other targets, KKR-CPA ∆E f orm vs. {c α } for (Mo-W) x Ta y (Ti-Zr) 1−x−y are shown in Fig. 2. Clearly, ∆E f orm for equiatomic case is too positive (+12.7 mRy), and decomposition is expected (with T M G c =1244 K from estimates in Hume-Rothery section). Our calculated SRO also indicates phase decomposition at spinodal T sp =1240 K, agreeing with T M G c , in the near-equiatomic and Ti-Zr-rich alloys (see Fig. 5). This predicted segregation is corroborated by our X-ray diffraction experiments, Fig. 3(a), that indicate presence of two (major/minor) phases. The phases were indexed as a disordered A2 phase with Im3m space group, and a minor phase of F d3m space group. The A2 lattice parameter was measured as 3.1713Å (std. dev: 0.0002Å). Figure 3(b) shows the SEM micrograph of the alloy with a two-phase alloy evident. The major phase (A2) is Mo, W and Ta rich, with small amounts of Ti and Zr incorporated in it. Given the higher melting temperatures of Mo, W and Ta, the major phase is likely to be the primary solidifying phase during the final step of casting. As this phase forms during casting, Ti and Zr are rejected into the surrounding liquid, which subsequently freezes. Hence, the minor phase is Ti and Zr rich and incorporates small amounts of the refractory metals. To visualize key alloying effects for these CSAs, we plot in Fig. 4 the electronic dispersion and projected total density of states (TDOS), referenced to each alloy's Fermi energy, E F . With disorder, dispersion exhibits broadening in E and k, showing that k is a "good" (on the scale of the Brillouin zone) but not an exact quantum number (as for zero-width, ordered bands); the width dk ∼ l −1 e (l e is the electron scattering length) and gives rise to increased residual resistivity, as may be calculated. 47 Guided by such details, we can improve CSA properties. The equiatomic alloy has a TDOS with E F not yet in the pseudogap between bonding and antibonding states (top Fig. 4), so this alloy does not satisfy the design criteria. The VEC (the average electrons per atom outside the closed shells of the component atoms) is a dominant factor in controlling the phase stability of the alloys. The electronic states present on/near the E F are chemically most active, which affect the chemical property of the alloy, i.e., more states at E F destabilizes the alloy. This means, for such cases, adding or removing electron we can manipulate the electronic properties very quickly. By integrating states from E F to the pseudogap for equiatomi case, 0.2 electrons are needed to fill bonding states and improved stability. More at.%Mo-W (VEC=6) adds electrons, moving E F up (Fig. 4), and ∆E f orm reduces to stabilize CSAs (Fig. 2). Adding small %Ta helps in altering states near E F : the flat bands near Γ in Fig. 4 are moves from E F for C7, lowering ∆E f orm , shown in supplementary figure 2 & 3. We show in Fig. 5 that the SRO changes from clustering in C6 to ordering in C7, while ∆E f orm ∼0. Here, ∆E f orm reduces quickly for MoW with a small %Ta and %TiZr, while bulk modulus (B) increases quickly (inset Fig. 2). Notably, the dispersion of A2-metals is canonical when scaled by bandwidth (inverse atomic size), and so the behavior of the alloy dispersion is fairly generic and predominantly determined by relative composition and size, hybridization, and band filling. To promote oxide-scale formation for protection, and lightweighting, Al is often added. In Fig. 4, 5%Al added at the expense of Ta to C6 (whose ∆E f orm =+0.70 mRy) increases disorder broadening (from Al sp-d hybridization) and causes d-state around Γ (predominantly TiZr) to again straddle at E F . This Fermi surface feature energetically destabilizes the alloy making ∆E f orm much more positive (+6.8 mRy, with T M G c =828 K), which is also visible in SRO showing strong clustering behavior at T sp =780 K (see Supplementary Figure 4). The most significant Al-Mo pair suggest that Al will segregate to surfaces due its faster kinetics, as needed for oxide formation, i.e., adding Al at the expense of Ta or TiZr decreases VEC and drops E F into localized d-states, reduces stability; so a balance must be struck by keeping some Ta and TiZr and making VEC high enough to be near ∆E f orm ∼0 but with a large B (Fig. 2). Unlike in other systems, Al is not generically a good A2 stabilizer, as it leads to larger electron scattering for reduced stability, increased resistivity, and decreased thermal transport, see [18] and references therein.

C. Chemical SRO:
From KKR-CPA linear response (see Methods), we predict (Fig. 5) Warren-Cowley SRO (or atomic pair correlations) α µν (k;T), whose largest peak at wavevector k 0 reveals the unstable (Fourier) modes to ordering, or clustering at Γ=(000). As an alloying guide, SRO identifies pairs driving the instability, and predicts the spinodal T sp , where α −1 µν (k 0 ; T sp )=0 signifying the absolute instability to this chemical fluctuation. 18,19 In real-space, pair probabilities are P µν ij = c i µ c j ν (1 − α ij µν ), with α i =j µν =0 for no SRO, and α<0 (α>0) indicates ordering (clustering) with bounds of −[min(c µ , c ν )] 2 (c µ c ν ) −1 ≤ α i =j µν ≤ 1. Near-equiatomic alloys in Fig. 5 have maximal SRO peaks in α µν (k 0 =Γ;T>T sp ) signaling spinodal (infinite wavelength) decomposition in specified pairs at T sp of 1240 K for C1, and at 500 K for C6. This 60% drop in T sp is unsurprising given that ∆E f orm reduces with Mo-W and Ta addition (Fig. 2). For C7, where ∆E f orm has become slightly negative due to movement of bands present at Γ away from E F , a weak incommensurate (long-period) ordering is found with SRO peak (Fig. 5) at 70% along N-H at k 0 =(0.85, 0.15, 0). This SRO arises from Fermi-surface nesting, 57 with contributions at a radius of |k 0 |∼0.86, as confirmed along Γ-H (Fig. 5). [SRO is B2 type if it peaks at k 0 =H={100}, commensurate with A2 lattice.] For theory and detailed examples, see [18] and [19]. We also show, in supplementary figure 4, that 5%Al addition to the C6 alloy instigates a clustering instability. The Al-Mo pair drives spinodal decomposition at T sp of 780 K, which shows the tendency of Al to phase separate from Mo, an indication that Al's clustering tendency might be helpful in promoting stable oxide-layer at high temperatures. These results indicate that alloying may improve oxidation behavior, just as for Fe-Cr with a narrow window for chromia formation. Clearly, KKR-CPA methods address profound electronic and alloying effects not possible from effective potentials, or methods that approximate disorder by ordered configurations.

D. Deformation Analysis:
Mechanical properties in CSAs have been studied at macro-and micro-scopic levels, 48,49 but deformation analysis is key to establish high-T structural candidates. We perform quasi-static uniaxial loading via MD simulations (see Methods) by deforming an ideal single-crystal alloy in small but finite steps and equilibrating after each step. For equiatomic case, <100> compression (Fig. 6a) reveals a smooth stress-strain curve signaling simple plastic flow. In contrast, C3 (0.425at.%Mo) ideal crystal has stress drops and strain-hardening triggered by <111> dislocations; a stress drop at 0.065 strain marks the initiation of dislocation with A2 Burgers vector, b = 1 2 <111>, from 77-2000 K. Snapshots of the evolution show that dislocations (edge and screw type) triggered these instabilities (Fig. 6b). The defect mobility is affected by local distortions caused by the different sizes and modulus of the solutes. The rise and drop in stress with increasing strain in the ideal crystal corresponds to the defect evolution where new dislocations occur after every major stress drop followed by strain-hardening due to dislocation interactions and drag. An investigation of the local structural environment (Fig. 6b) reveals deviation from perfect A2, as yielding occurs for 300 K. Shear bands (black) are promoted, denoting deformed regions with higher compression. At very high strain the interplay between edge and screw dislocations can be visualized via the band dynamics (Fig. 6b-4). For engineering, Young's modulus E =3(1-2ν)B is pertinent, so Poisson's ratio ν is also key. Small-deformation MD simulations determined E and ν at 300 K, and E vs. T was found from the elastic stress-strain curve (Fig. 6a). The KKR-CPA energy versus a (A2 lattice constants) at {c α } determines the equilibrium a 0 , ∆E f orm , and B (Fig. 6d), all used in Fig. 2. We compare temperature dependence of E, in Fig. 6c, calculated from MD, KKR-CPA (using Grüneisen approximation at low T), and experiments for commercial Mo-rich TZM alloy, 3 which again validate theory results. As MD is performed on ideal crystals, an "ideal" yield strength is obtained, with a qualitative relative change versus temperature. To confirm our predicted E in equiatomic C1, we performed indentation on samples prepared by arc-melting (see Methods), in Fig. 6(c), and show that value from measurement 104±12 GPa at 300 K compared very well to our predictions 115 GPa from KKR-CPA and 120 GPa from MD. As B changes slowly for binary Mo-W (Fig. 6d), the Poisson effect (variation of ν) controls strength, which requires W-or Mo-rich alloys for larger E values (Fig. 6d). Similarly, for quinary (Mo z W 1−z ) 0.85 Ta 0.10 (TiZr) 0.05 , we find that C10 (z = 0.50) has strength similar to C4 (Fig. 6c). Whereas we predict a region around C10 (z = 0.05, highlighted in Fig. 6d, which is perpendicular to plane in Fig. 2 at C10) that shows both enhanced stability and E (Fig. 6c). For these CSAs, we find 3× larger E than high-entropy alloys at 300 K, and alloys like C10 have a much larger, less temperature-dependent modulus (Fig. 6c) above 500 K (2.3× at 2000 K) over existing commercial TZM alloys, and lie midway between pure Mo and W, unlike TZM alloys. Finally, one notable point, albeit not surprising, the classical MD simulations fail to represent properly the alloys that crossover from ∆E f orm positive to negative (e.g., C4-to-C7 or C3-to-C10) in which electronic dispersion (not addressed by semi-empirical potentials) is controlling the materials physics. Hence, we plot only C10 KKR-CPA values in Fig. 6c, as MD values of E vs. T for C10 and C10 are similar to C4, whereas C10 values from first-principles increases over C10, as expected from Fig. 6(d). From a design perspective, in general, alloys in complex solid-solution alloys have superior properties over near-   Fig. 2) shows enhanced E (above Mo) and better T-dependent slope (Fig. 6c) from a favorable Poisson effect (green area in Fig. 6d). Above 500 K, C10 has larger, less T-dependent E than TZM (2.3× at 2000 K).
equiatomic alloys (so-called high-entropy alloys), although the design space becomes enormous. Using a first-principles KKR-CPA, we predicted the relative phase stability, dispersion, short-range order (i.e., incipient long-range order, including T c ) and its electronic origin, and mechanical properties over all compositions as a design guide. Using electronic alloy design concepts and criteria, we identified higher strength refractory (Mo-W)-Ta-(Ti-Zr) alloys from materials physics and engineering perspectives. Temperaturedependent deformation (most relevant the elastic behavior) in selected set of alloys was modeled using classical MD simulations, validated from first-principles data; we also identi-fied failures in classical potentials that arose from dispersion effects.
Based on our calculation, we designed a Mo-rich region of improved stability with enhanced Young's moduli over high-entropy alloys, as we confirmed experimentally, and an improved temperature-dependence above 500 K (2.3× at 2000 K) over existing commercial alloys. Our electronicstructure approaches and analysis of alloying and stability (formation energies, dispersion, short-range ordering) highlights how instructive these details are in guiding design. The techniques are quite general for assessing any arbitrary complex solid-solution alloys, where alloying and non-trivial electronic effects play a key role.

E. DFT Methods:
KKR electronic-structure is used with the coherent-potential approximation (CPA) to handle chemical disorder; 16,17 screened-CPA addresses Friedel screening from charge-correlations. 17 Scalar-relativistic effects are included (no spin-orbit). Generalized gradient approximation to exchange-correlation was included through use of libXC libraries. 51 CSAs require only 1-atom (2atom) cells for A1,A2 (A3). Brillouin zone (BZ) integrations were performed with Monkhorst-Pack k-point method, 52 with 12×12×12(6) for A1,A2 (A3) meshes. We used 300 k -points in the irreducible-BZ to visualize dispersion along symmetry lines. Each scatterer's radii were defined by neutral "atoms-in-cell", with interstitial divided proportionally to each scatterer, to improve radial density representation near saddle-points in the electronic density. 53,54 We chose Lmax=3 spherical-harmonic basis to include s, p, d and f-orbital symmetries. Shallow core states were included in the valence in all calculations. A variational potential zero v0 was used to yield kinetic energies nearing those of full-potential methods. 55 For self-consistent densities, complexenergy contour integration 56 used 20-point Gauss-Legendre semicircular contour.

F. Chemical SRO:
From KKR-CPA linear-response, we calculate SRO parameters, αµν (k;T), for µ-ν pairs, 18,19 as detailed elsewhere. [20][21][22][23]57 Dominant pairs driving SRO are identified from pair-interchange energies, S (i,j) µν (T ), or curvature (concentration 2 nd -variation) of the KKR-CPA grand potential, yielding energy cost for concomitant fluctuations of c i µ , c j ν at atomic sites i, j. S (2) µν (k; T ) reveals the unstable (Fourier) modes with ordering wavevector k0 (or clustering at (000)), identifies the origin for phase transitions, and dictates the SRO: µν (k; T )]. The spinodal temperature, where α −1 µν (k0; Tsp)=0, signifies an absolute instability to this fluctuation and provides an estimate for T M G c or T od c . 18,19,58 For N>2, pairs driving ordering (clustering) will not necessarily be the same pairs that peak in the SRO due to the matrix inversion that relates them (Fig. 5).

G. MD Simulations:
Deformation is evaluated using Large-scale Atomistic/Molecular Massively Parallel Simulator (LAMMPS) package. 59 The KKR-CPA structural parameters are used to validate potentials for finite-T modeling. The force-field parameters for the quinary are established from available ternary EAM potentials. 60 We verified similar hybrid potential parameter combinatorial technique for high-entropy aloy, like Al10CrCoFeNi. 19 In A2 lattice, with dimensions (30×30×30)a (54,000 atoms), we distributed 5 elements via composition to form (Mo-W)-Ta-(Ti-Zr) solid solutions. Initially the lattice was melted at 4,000 K for 90 ps, followed by a quench to 300 K within 10 ns. Uniaxial deformation was performed after equilibrating and relaxing the structure at high strain, as detailed in supplementary video 1. 19

H. Synthesis & Characterization:
The equiatomic MoWTaTiZr was synthesized by arc-melting pellets of elemental powder blends (Alfa Aesar, purity ≥ 99.9%) in an ultra-high purity argon atmosphere on water-cooled copper hearth. Powders were used to reduce the large macro-segregation that occurs during casting when using elemental chips. With the significant difference in melting temperatures (3695 K for W vs. 1941 K for Ti) a three-step melting process was adopted.
Step 1: W and Ta powders were mixed thoroughly in a SPEX 8000 mill, and pressed using a Carver hydraulic press; and the pellet was then arc-melted.
Step 2: elemental blends of Mo, Ti and Zr were similarly mixed, pressed and arc-melted.
Step 3: both arc-melted buttons were re-melted together for a total of four times to ensure better homogeneity.
Phase analyses were carried out using a Philips PANalytical X-Ray Diffractometer (XRD), in a Bragg-Brentano geometry using Cu-Kα radiation. Microstructure and phase compositions were analyzed using a FEI Helios NanoLab G3UC Scanning Electron Microscope (SEM), equipped with Oxford Energy Dispersive Spectroscopy (EDS) system. Accelerating voltages of 10-15 kV were employed for imaging and compositional analyses. Compositions were measured at 7 different locations for each phase, with the average composition and X-ray diffraction shown in Fig. 3. The diffraction pattern indicated the presence of two phases, indexed as a bcc (A2) phase with lattice parameter 3.1713(2)Å and a minor phase with F d3m space group, like for B32 (NaTl prototype) or C15 (MgCu2 prototype) structures, with lattice constant 7.6148(9)Å.

I. Nanoindentation:
Nanoindentation utilized a tribo-indenter HYSITRON TI-900 with a Berkovich (3 µm) tip. The indenter control module applies a trapezoidal load on the sample for 10 s, followed by 5 s rest, and unloads in 10 s. To calibrate the sample measurements, which also determines the best applied load and optimum contact depth, the alloy was scanned (15 measurements) on an arbitrarily chosen sample location to optimize for force vs. displacement. For a minimum of 200 nm of contact depth, a 6000 µN load was found to suffice. With these set, we indented the alloy at 30 manually chosen locations to measure the sample's elastic response; for the equiatomic case, the mean values were: Young's modulus of 103.73 ± 11.49 GPa, hardness of 4.6 ± 0.34 GPa, and contact depth of 231.18 ± 8.74 nm.