Introduction

β-Titanium (β-Ti) alloys that are free from vanadium (V), nickel (Ni) and aluminium (Al) are one of the most promising groups of metallic biomaterials for orthopaedic applications. By definition, β-Ti alloys refer to those which contain sufficient β-stabilizing element(s) to enable the β-phase to be retained in a metastable condition by water quenching.1 They are generally divided into metastable (age-hardenable) and stable β-Ti alloys (not age-hardenable).2 The critical minimum amount of β-stabilizers needed to retain 100% β-phase on water quenching is denoted as βc, while the minimum content of β-stabilizers needed to form stable β-Ti alloys is referred to as βs. Metastable β-Ti alloys are those that fall in between βc and βs, within which α can precipitate upon ageing. As the values of βc and βs depend on alloy chemistry, a parameter known as the Mo equivalent (Mo-Eq) has been defined and commonly used as a general guideline to assess the β-phase stability of a Ti alloy for both alloy design and processing.3,4,5 It is therefore an important design parameter for Ti alloys. This paper first examines all existing Mo-Eq expressions and then focuses on the development of a new Mo-Eq equation proposed specifically for the Ti-Nb-Zr system. This is followed by a detailed experimental assessment and comparison with Morinaga’s d-electron theory for alloy design.

Existing Mo-Eq Expressions and Their Deficiencies

The concept of the Mo equivalence was proposed by Molchanova [Phase diagrams of titanium alloys (Figs. 193 and 194, p. 158)],3 which was a translation from Atlas Diagramm Sostoyaniya Titanovykh Splavov (Israel Program for Scientific Translations, Jerusalem, 1965). Boyer et al.4 presented the data shown in Fig. 193 of Ref. 3 as the Mo equivalence of a titanium alloy in the well-known Materials properties handbook: titanium alloys (1994) in the following form:

$$ \begin{aligned} \left( {\text{Mo-Eq}} \right)^{{{\text{Molchanova}}\, [ 3 ]}} & = 1.0\left( {\% {\text{Mo}}} \right) + 0.2\left( {\% {\text{Ta}}} \right) + 0.28\left( {\% {\text{Nb}}} \right) \\ & \quad + \,0.4\left( {\% {\text{W}}} \right) + 0.67\left( {\% {\text{V}}} \right) + 1.25\left( {\% {\text{Cr}}} \right) \\ & \quad + \,1.25\left( {\% {\text{Ni}}} \right) + 1.7\left( {\% {\text{Mn}}} \right) + 1.7\left( {\% {\text{Co}}} \right) \, \\ & \quad + \,2.5\left( {\% {\text{Fe}}} \right) \\ \end{aligned} $$
(1)

Bania5 then further discussed the Mo-Eq concept and expressed it in the form of Eq. 2 in which the prefix coefficient of each β-stabilizer represents the ratio of the βc value of Mo (i.e., 10 wt.%) to that of respective β-stabilizer. As an α-stabilizer, Al is also included to show its opposing effect when it is present in a β-Ti alloy.

$$ \begin{aligned} \left( {\text{Mo-Eq}} \right)^{\text{Bania [5]}} & = \, 1.0\left( {\% {\text{Mo}}} \right) + 0.67\left( {\% {\text{V}}} \right) + 0.44\left( {\% {\text{W}}} \right) \\ & \quad + \,0.28\,\left( {\% {\text{Nb}}} \right) + 0.22\left( {\% {\text{Ta}}} \right) + 2.9\left( {\% {\text{Fe}}} \right) \\ & \quad + \,1.6\left( {\% {\text{Cr}}} \right) + 0.77\left( {\% {\text{Cu}}} \right) + 1.11\left( {\% {\text{Ni}}} \right) \\ & \quad + \,1.43\left( {\% {\text{Co}}} \right) + 1.54\left( {\% {\text{Mn}}} \right) - 1.0\left( {\% {\text{Al}}} \right) \\ \end{aligned} $$
(2)

Equation 2 overrides the original Eq. 1 as a commonly used expression to estimate the Mo-Eq of a Ti alloy. Apart from Eq. 2, three other expressions for Mo-Eq have also been proposed,6,7,8 namely:

$$ \begin{aligned} \left( {\text{Mo-Eq}} \right)^{{{\text{Zhou}}\, [ 6 ]}} & = 1.0\left( {\% {\text{Mo}}} \right) + 0.74\left( {\% {\text{V}}} \right) + 0.5\left( {\% {\text{W}}} \right) + 0.39\left( {\% {\text{Nb}}} \right) \\ & \quad + \,0.28\left( {\% {\text{Ta}}} \right) + 2.2\left( {\% {\text{Fe}}} \right) + 1.69\left( {\% {\text{Cr}}} \right) + 0.85\left( {\% {\text{Cu}}} \right) \\ & \quad + \,1.22\left( {\% {\text{Ni}}} \right) + 1.57\left( {\% {\text{Co}}} \right) + 1.69\left( {\% {\text{Mn}}} \right) \\ \end{aligned} $$
(3)
$$ \begin{aligned} \left( {\text{Mo-Eq}} \right)^{{{\text{Mo}}\, [ 7 ]}} & = 1.0\left( {\% {\text{Mo}}} \right) + 0.73\left( {\% {\text{V}}} \right) + 0.5\left( {\% {\text{W}}} \right) + 0.31\left( {\% {\text{Nb}}} \right) \\ & \quad + \,0.24\left( {\% {\text{Ta}}} \right) + 2\left( {\% {\text{Fe}}} \right) + 1.69\left( {\% {\text{Cr}}} \right) \\ & \quad + \,1.29\left( {\% {\text{Ni}}} \right) + 1.16\left( {\% {\text{Co}}} \right) + 1.69\left( {\% {\text{Mn}}} \right) \\ \end{aligned} $$
(4)

and

$$ \begin{aligned} \left( {\text{Mo-Eq}} \right)^{{{\text{Wang}}\,[8]}} & = 1.0\left( {\% {\text{Mo}}} \right) + 1.25\left( {\% {\text{V}}} \right) + 0.59\left( {\% {\text{W}}} \right) + 0.28\left( {\% {\text{Nb}}} \right) \\ & \quad + \,0.22\left( {\% {\text{Ta}}} \right) + 1.93\left( {\% {\text{Fe}}} \right) + 1.84\left( {\% {\text{Cr}}} \right) + 1.5\left( {\% {\text{Cu}}} \right) \\ & \quad + \,2.46\left( {\% {\text{Ni}}} \right) + 2.67\left( {\% {\text{Co}}} \right) + 2.26\left( {\% {\text{Mn}}} \right) + 3.01\left( {\% {\text{Si}}} \right) \\ & \quad + \,0.3\left( {\% {\text{Sn}}} \right) + 0.47\left( {\% {\text{Zr}}} \right) -1.47\left( {\% {\text{Al}}} \right) \\ \end{aligned} $$
(5)

Equations 24 are similar and were all formulated on the basis of experimental data, where Zr was treated as a neutral element and therefore excluded. It should be pointed out that no information was given about the sample size used in the water-quenching experiments, which affects the resultant microstructure. However, Eq. 3 considered that the βc value of Mo should be 11 (rather than 10), which was also mentioned in Ref. 3. The slight difference between Eqs. 2, 3 and 4 resulted from experimental and analytical errors. In contrast, Eq. 5 combines both experimental and phase diagram data. For example, the coefficients for Nb and Ta were taken from Eq. 2, while the coefficient for Zr was derived from the Ti-Zr phase diagram.8 As a result, Zr was treated as a β-stabilizer that is more potent than Nb (0.28) and Ta (0.22), but no experimental support was given. Limited evidence suggests9,10 that the β-stabilizing effect of Zr depends on the quantity of other β-stabilizers present in the alloy.

Table I summarizes nine different Ti-Nb-Zr alloys reported in the literature with identified phase constituents in the water-quenched state.9,11 They provide an experimental basis for a detailed assessment of Eqs. 2 and 5. As can be seen from Table I, despite the high Mo-Eq values predicted by Eq. 5, e.g., 24.12 for Ti-12Nb-30Zr (at.%), α″-martensite still forms on water quenching. In other words, Eq. 5 significantly overestimates the β-stabilizing effect of the alloying elements in these alloys. With regards to Eq. 2, it substantially underestimates the overall β-stabilizing effect of the alloying elements in Ti-Nb-Zr alloys. For example, Eq. 2 predicts Mo-Eq values of 3.64 and 5.43 for Ti-10Nb-45Zr (at.%) and Ti-14Nb-30Zr (at.%), respectively, which are much lower than 10. However, the β-phase was fully retained in both alloys upon water quenching. Equations 2 and 4 offer similar predictions, while Eq. 3 also shows a serious underestimation of the β-stabilizing effect. For instance, the Mo-Eq of Ti-14Nb-30Zr (at.%) is predicted to be 7.57 but full β phase is retained. So is the case for Ti-10Nb-45Zr (at.%), where Mo-Eq = 5.07 but full β phase is retained. On the other hand, Eq. 3 overestimates the β-stabilizing effect in Ti-20.9Nb-9Zr (at.%) (clear formation of α″, although Mo-Eq = 12.36). In summary, none of the existing Mo-Eq equations offers predictions that are consistent with experimental observations for Ti-Nb-Zr alloys. In addition, when all these predictions based on Eqs. 25 are compared with the d-electron theory diagram, there is also substantial inconsistency.

Table I Metastable β-region for Ti-Nb-Zr alloy systems and their Mo-Eq9,11

New Mo-Eq Expression for Ti-Nb-Zr Alloys

We first determine the Mo-Eq coefficient of Nb. Three different values were used previously in Eqs. 25, namely, 0.28%, 0.31% and 0.39%. Recent experimental work has defined that the βc value for Ti-Nb binary alloys is 42 wt.%.12 Consequently, the Mo-Eq coefficient of Nb can be defined as

$$ \psi = \left( {10/42} \right) = 0.238 $$
(6)

We then need to properly consider the β-stabilizing effect of Zr. As can be seen from Table I, a small increase (e.g., 1–3 at.%) in the content of Zr or Nb + Zr can change the phase constituents from α″ + β to single β under water-quenching conditions. For instance, a single β-phase was retained in the Ti-22Nb-4Zr (at.%) alloy but not in the Ti-22Nb-2Zr (at.%) alloy. Similarly, a single β-phase was retained in the Ti-14Nb-30Zr (at.%) alloy but not in the Ti-12Nb-30Zr (at.%) alloy. Considering that Nb and Zr are generally weak β-stabilizers, it can be assumed that the Mo-Eq values of the Ti-Nb-Zr alloys listed in Table I are all in the vicinity of the critical value, i.e., 10 wt.%. Consequently, the β-stabilizing effect of Zr in these Ti-Nb-Zr alloys can be estimated by the quantity, Δ = 10 − ψ × Nb (wt.%). Figure 1 shows the dependence of the quantity, Δ = 10 − ψ × Nb (wt.%), on the Zr content for all the Ti-Nb-Zr alloys listed in Table I. The resultant significant linear relationship identifies the contribution of Zr in stabilizing the β phase in Ti-Nb-Zr alloys considered in Table I. Accordingly, a new Mo-Eq for Ti-Nb-Zr alloys can be defined from Fig. 1 as

$$ \left( {\text{Mo-Eq}} \right)_{\text{Ti -Nb -Zr}} = 0.238{\text{Nb}}\,\left( {{\text{wt}} .\% } \right) + 0.11{\text{Zr}}\left( {{\text{wt}} .\% } \right) + 0.97 $$
(7)

Equation 7 was applied to all the Ti-Nb-Zr alloys listed in Table I. Good consistency was observed with the experimental observations, where, when the new Mo-Eq parameter is close to or above 10, the alloy contains β-phase as the primary phase (along with some other secondary phases) after water quenching, while, when it is clearly below 10, martensite is the dominant phase.

Fig. 1
figure 1

Dependence of the quantity of Δ = 10 − ψ × Nb (wt.%) on Zr content in Ti-Nb-Zr alloys

Applying the New Mo-Eq to the Design of Ti-Nb-Zr Alloys and Validation

Table II summarizes the Mo-Eq values predicted by Eqs. 25 and 7 for six newly designed Ti-Nb-Zr alloys in this study (the first six alloys in Table II) and those available in the literature. Ingot samples of these alloys were made by arc-melting of a mixture of pure metals (purity > 99.9 wt.%) in a Ti-gettered high-purity argon atmosphere. The ingots were re-melted four times to ensure chemical homogeneity. Then, the melted alloys were cast into a copper mold with a cylindrical cavity of 10 mm in diameter and 60 mm in depth. The as-cast samples with a diameter of 10 mm and thickness of 2 mm were homogenized at 1000°C for 4 h and then solution-treated at 900°C for 0.5 h, followed by water quenching. Samples were ground with SiC grinding papers of up to 4000 grit for microstructural characterization. Final polish was performed with a mixture of colloidal silica and H2O2. The samples were subsequently cleaned in an ultrasonic bath using ethanol and etched with Kroll’s reagent. Phase identification was conducted using a Bruker D8 Advance x-ray diffractometer (XRD) with a Cu Kα radiation source at room temperature, and the resultant XRD data was analyzed by MAUD software.

Table II Newly designed metastable Ti-Nb-Zr alloys and corresponding Mo-Eq values

Figure 2 shows the XRD data of six water-quenched Ti-Nb-Zr alloys and optical micrographs of two of them. Ti-10Nb-26Zr (at.%) was predicted to have martensite due to its low Mo-Eq (8.6, Table II) while the rest should essentially be single β-phase in the water-quenched state. Both the XRD data and microstructural observations have confirmed these predictions, where the β-phase is predominant in each of the other five alloys (e.g., Figure 2c) plus a small amount of secondary phases (α″, ω) while Ti-10Nb-26Zr (at.%) contains significant α″-martensite (Fig. 2d). The newly defined Mo-Eq expression by Eq. 7 is in good agreement with the experimental observations of each alloy.

Fig. 2
figure 2

XRD results of (a) two Zr-based Ti-Nb-Zr alloys and (b) four Ti-based Ti-Nb-Zr alloys. (c) Water-quenched Ti-15Nb-40Zr (at.%), nearly single β. (d) Water-quenched Ti-10Nb-26Zr (at.%), noticeable α″-martensite

With regard to the alloys from the literature (Table II), they were all reported to be nearly single β-phase alloys (i.e., may contain a small amount of α″ and/or ω).13,14,15 By the definition of the Mo-Eq, their Mo-Eq values should be equal to or above 10 wt.%. As can be seen from Table II, the new Mo-Eq (Eq. 7) also shows excellent consistency with experimental observations.

The Newly Proposed Mo-Eq and the d-Electron Theory

The d-electron theory describes the potential relationship between phase stability and elastic properties of Ti alloys through two electronic properties, bond order, Bo and d-orbital energy level, Md.12,16 It was developed by Morinaga16 on the basis of the DV-Xα molecular orbital calculations of electronic structures. The theory has become an essential fundamental tool for the design of β-Ti alloys.1,12 In this approach, Bo denotes the covalent bond strength between Ti and a specific alloying element and Md correlates the electronegativity and metallic radius of that element. Bo and Md values for all the alloying elements of Ti can be found in Ref. 16. Using these Bo and Md values, the compositional average of Bo and Md can be calculated based on alloy composition through the following equations.16

$$ \overline{{B_{\text{o}} }} = \sum X_{i} \left( {B_{\text{o}} } \right)_{i} $$
(8)
$$ \overline{{M_{d} }} = \sum X_{i} \left( {M_{d} } \right)_{i} $$
(9)

where \( X_{i} \), (Bo)i and \( (M_{d} )_{i} \) are the atomic fraction, bond order and d-orbital energy level of the given element, respectively. Then, a \( \overline{{B_{\text{o}} }} \)\( \overline{{M_{d} }} \) diagram can be plotted by taking the \( \overline{{B_{\text{o}} }} \) and \( \overline{{M_{d} }} \) values of different Ti alloys. Several sets of Ti alloys are used to define the three regions, α, α + β, β, of this \( \overline{{B_{\text{o}} }} \)\( \overline{{M_{d} }} \) diagram. Along the (α + β)/β boundary in the \( \overline{{B_{\text{o}} }} \)\( \overline{{M_{d} }} \) diagram, the β-phase is considered metastable.

Figure 3 shows the \( \overline{{B_{\text{o}} }} \)\( \overline{{M_{d} }} \) diagram plotted for this study on the basis of Refs. 12 and 14. A triangle is drawn by connecting the Bo and Md values for pure Ti, Nb and Zr, which define the regime of the Ti-Nb-Zr alloys, similar to ternary phase diagrams. Each side of this triangle indicates the compositional average values of Bo and Md of the Ti-Nb, Ti-Zr and Zr-Nb binary alloys. The blue dotted line represents the boundary of critical compositions for metastable β alloys according to the inventor of the d-electron theory.12 This boundary may shift upwards or downwards depending on the exact alloying compositions. In the Ti-Nb-Zr system, adding Zr was found to shift this boundary downwards,14 as shown by the black dashed line in Fig. 3. The Mo-Eq values predicted by Eq. 7 are superimposed onto Fig. 3, which indicates that the newly developed Mo-Eq expression is also in good agreement with the d-electron theory in predicting the β-phase stability. For example, Ti-Nb-Zr alloys having Mo-Eq values higher than 10 wt.% are situated above the black dashed line and vice versa. The new Mo-Eq is thus consistent with the d-electron theory for Ti-Nb-Zr alloys, which was not the case for all the other four Mo-Eq expressions.

Fig. 3
figure 3

Ti-Nb-Zr alloys in the \( \overline{{B_{\text{o}} }} \varvec{ } \)\( \overline{{M_{\text{d}} }} \) map13,14,15,16

Summary

A new Mo-Eq expression has been formulated for the Ti-Nb-Zr system, i.e., (Mo-Eq)Ti-Nb-Zr = 0.238Nb (wt.%) + 0.11Zr (wt.%) + 0.97, which shows good consistency with experimental observations as well as literature data. In addition, it is highly consistent with the well-established d-electron theory in predicting the β-phase stability of various Ti-Nb-Zr alloys. In contrast, none of the existing four Mo-Eq expressions is consistent with the d-electron theory, as well as experimental observations. However, caution should be exercised when the alloy composition goes significantly beyond the bounds of the chemistries discussed in this study. It is proposed that the approach developed in this work can also be applied to the assessment of the β-phase stability in other Zr-containing Ti alloys.