Diffusivities and Atomic Mobilities in BCC Ti-Fe-Cr Alloys

In this research, the diffusion behaviors within the Ti-Fe-Cr ternary system were examined at the temperatures of 1273 K and 1373 K through the diffusion couple technique. This study led to the determination of both ternary inter-diffusion and impurity diffusion coefficients in the body-centered cubic (bcc) phase for the Ti-Fe-Cr alloy, utilizing the Whittle–Green and Hall methods. The statistics show that the average diffusion coefficients D˜FeFeTi and D˜CrCrTi measured at 1273 K were 1.34 × 10−12 and 3.66 × 10−13, respectively. At 1373 K, the average values of D˜FeFeTi and D˜CrCrTi were 4.89 × 10−12 and 1.43 × 10−12. By adopting the CALPHAD method, a self-consistent database for atomic mobility in the bcc phase of the Ti-Fe-Cr system was established. This database underwent refinement by comparing the newly acquired diffusion coefficients with data from the existing literature. Diffusion simulations for the diffusion couples were performed, drawing on the established database. The error between the simulated diffusion coefficient and the experimental measurement data is within 15%, and the simulated data of the component distance distribution and diffusion path are in good agreement with the experimental data. The simulations generated results that aligned well with the observed experimental diffusion characteristics, thereby affirming the reliability and accuracy of the database.


Introduction
Over the past few decades, the excellent performance of titanium and its alloys, such as their specific strength, relatively low density, superior resistance to corrosion, and favorable biocompatibility, have positioned them as materials of choice for both the aerospace and biomedical fields [1][2][3].One of the main challenges of their broader application, compared to steel and aluminum, is the high cost of production [4][5][6].This issue largely arises from the reliance on expensive β-type elements, like V, Mo, and Nb, instead of more affordable ones, such as Fe, Cr, and Mn.Opting for cost-effective β-stabilizers is a strategic method to decrease alloy expenses [7].Notably, Cr has been identified as an efficient eutectoid element [8], and Fe is considered an economical option for achieving high-strength, lowdensity alloys.Consequently, several economical Ti-Fe-Cr alloys have been developed, such as Ti-4.3Fe-7.1Cr[9], Ti-4.2Fe-6.9Cr[10], and Ti-3.0Al-8.5Cr-2.0Fe[11].
The properties of titanium alloys are closely tied to the microstructures that develop from their processing and subsequent heat treatments.A considerable body of research has been focused on exploring the diffusion mechanisms of alloying elements within titanium alloys, particularly within binary and ternary compositions, to understand the underlying processes that govern microstructure formation in β-titanium alloys [12][13][14][15].
Materials 2024, 17, 1927 3 of 18 within the bcc solid solution region of the Ti-Fe-Cr system at temperatures of 1273 K and 1373 K, as illustrated in Table 1.Alloys comprising solely Ti, binary mixtures of Ti-Fe, Ti-Cr, and the ternary Ti-Fe-Cr, were fabricated using arc melting in an argon environment.The melting stock included granules of Ti, Fe, and Cr, each with a purity of 99.99 wt%.The materials used in this experiment were provided by Shanghai Chunming Metal Materials Co., Ltd. in Shanghai, China.To achieve homogeneity of the compositions across the alloys, the melting procedures of all alloys were carried out six times.Following this, the ingots were subjected to vacuum annealing at 1473 K for 30 h before being quenched in ice water.This specific regimen led to the formation of alloys with a grain size averaging several millimeters, thereby significantly reducing the role of grain boundary diffusion.
Table 1.Terminal compositions of the diffusion couples.
Following the annealing process, the diffusion couples were precisely divided along the diffusion path with the assistance of WEDM.This step was succeeded by their preparation through a series of standard metallography techniques, which included mounting, grinding, and polishing to achieve a suitable surface for examination.To further dissect the compositional gradients within these couples, electron probe microanalysis (EPMA, using the JXA-8230 model by JEOL, Tokyo, Japan) was employed.This analysis was conducted under specific conditions: a 15 kV operating voltage, a 20 nA beam current, and a 40-degree take-off angle, ensuring the detailed and accurate measurement of the composition-distance profiles.

Extraction of Inter-Diffusion Coefficients
In the context of one-dimensional diffusion, the diffusion coefficients within multicomponent metallic solid solutions from the composition-distance profiles were calculated using the Matano-Kirkaldy method [44].Kirkaldy [45] introduced the Boltzmann parameter λ = z/ √ t, where z represents the spatial position and t denotes the time, initially applying this parameter at the boundary condition when t = 0, The diffusion flux of an element i, symbolized by J i , correlates with its concentration, c i , where c iL stands for the element's concentration at the diffusion couple's left boundary.The parameter z reflects the distance from this boundary, z M denotes the position of the Matano plane, and t represents the elapsed time of diffusion.The symbol D k ij refers to the diffusion coefficient of element i in medium k, influenced by the concentration gradient of element j.Within a ternary alloy system, it is essential to identify four independent diffusion coefficients: D 3  11 , D 3 22 , D 3 12 , and D 3 21 .To accurately determine these coefficients, employing two diffusion couples that converge at a singular intersection point is crucial.
Identifying the Matano plane, which is established as the zero-distance point, is the initial step in this method, yet this task is known for being both time-intensive and prone to inaccuracies.To address these challenges, Whittle and Green [46] proposed the use of a normalized concentration variable: In the study of diffusion processes, the concentration levels of species at the diffusion couple's boundaries, denoted as c iL and c iR for the left and right ends, respectively, play a critical role.To refine the analysis of inter-diffusion fluxes, represented by J i , a novel approach involves the utilization of a newly introduced variable, Y i .This methodology facilitates the recalibration of the fluxes, enhancing the precision of the diffusion analysis.
where Y * i is the Y i value at specific position z * .Combining Equations ( 1) and (3), we can obtain in the Ti-Fe-Cr system: Given the scarcity of dependable information regarding how the composition influences the molar volume within Ti-Fe-Cr alloys, this research proceeded under the premise of a uniform molar volume across the board.The expectation is that the potential deviations stemming from this premise are sufficiently minor, thereby residing within the acceptable precision range for diffusion coefficients obtained via the Whittle-Green method.

Extraction of Impurity Diffusion Coefficients
Utilizing the Hall method [47], this study computed the impurity diffusion coefficients for chromium in the titanium-iron alloy and iron in the titanium-chromium alloy based on the final compositions of the diffusion pairs.This approach is notably effective for calculating inter-diffusion coefficients at the extremities of composition-distance profiles, where the Matano-Kirkaldy method encounters limitations due to the slope of the curve approaching zero, rendering the coefficients infinite.Initially, the composition profiles underwent a transformation to achieve a near-linear relation, depicted as µ = hλ + k.This step involved incorporating µ = er f −1 (2Y − 1) and the Boltzmann parameter λ = z/ √ t.Following this, the analysis proceeded to extract parameters h 1 , k 1 for the left and h 2 , k 2 for the right sides of the diffusion pairs through linear fitting.Conclusively, this study utilized specific equations to accurately calculate the inter-diffusion coefficients at the endpoints of the diffusion couples.
D(x′) specifically refers to the inter-diffusion coefficient at the left or right end of the profiles obtained when Y(x) tends toward 0 or 1.If the composition of an element is 0, the impurity diffusion coefficient of that element is as follows:

Atomic Mobility and Diffusion Coefficient
In the context of multi-component substitutional solutions, the diffusion flux experienced by any given species is directly affected by the thermodynamic forces at play.These forces manifest as gradients in chemical potential across the solution.Employing a synthesis of Fick's first law with Onsager's relational expressions enables the quantification of species-specific diffusion flux, symbolized by J i .This approach elucidates the intricate relationship between thermodynamic gradients and the resultant diffusion behaviors.
By incorporating the inter-diffusion coefficient D n ij , it becomes feasible to articulate the diffusion flux, J i , leveraging the phenomenological coefficient L ij alongside the chemical potential of species j, µ j , and the parameter z, representing distance.
In scenarios where diffusion primarily occurs through mono-vacancy exchange mechanisms, the Einstein relation serves to clarify the correlation between the species-specific tracer diffusion coefficient, denoted as D * i , and the corresponding atomic mobility, represented by M i .This pivotal relation sheds light on the influence of atomic mobility on the diffusion coefficient within contexts where vacancy exchange is the key mode of atom movement.
Within the context where mono-vacancy exchange predominates as the chief diffusion process, the Einstein relation provides a clarifying linkage between the tracer diffusion coefficient, D * i for a given species i, and its atomic mobility, M i .This relation demystifies how the mobility of atoms influences the diffusion coefficient in materials where vacancymediated diffusion is the primary mechanism.
In a substitutional solution phase with a constant volume, the interdiffusion coefficient linked to atomic mobility is articulated through the following formula, where R represents the gas constant, and T the temperature in kelvins: The Kronecker delta δ ki equals 1 when i = k, and is 0 otherwise; x i , x j , x k , and x n denote the mole fractions of species i, j, k, and n, respectively; M k represents the atomic mobility of species k.Initially introduced by Andersson and Ågren [48] and further elaborated by Jönsson [49], the formula for the atomic mobility M i of species i in relation to temperature within disordered solid solutions is as follows: The frequency factor for species i is denoted by M 0 i , and the activation energy is represented by Q i .The term mg Ω , indicative of the ferromagnetic contribution, is assumed to be 1, reflecting the absence of ferromagnetic transformation in the Ti-Fe-Cr alloys.Given that ∆G 12) can be reformulated as: Using the CALPHAD approach, the parameter ∆G φ i is acknowledged to vary with the composition, and it is formulated using the Redlich-Kister polynomial [50] as follows: For species designated as p, q, and v, the corresponding mole fractions are labeled as x p , x q , and x v .Furthermore, the symbol Φ p i denotes the mobility of species i when it is in its pure form, labeled as p.The parameters that describe interactions in mixtures with two and three components are respectively notated as (r) Φ p,q i for binary interactions and (s) Φ p,q,v i for ternary interactions.It is pertinent to note that within the framework of binary systems, the equation labeled as Equation ( 14) disregards the ternary interaction term, treating its value effectively as zero.

Composition-Distance Profiles and Diffusion Paths
Figure 2a showcases the microstructural view of the A1 diffusion couple, captured via backscattered electron microscopy, after a 44 h thermal treatment at 1273 K. Figure 2b illustrates the composition-distance profiles for the same couple, obtained through EPMA analysis.All diffusion couples within the Ti-Fe-Cr system are observed in a consistent bcc phase region, with the A1 microstructure serving as a prime example.The EPMA examination along a predetermined direction yielded accurate composition profiles, revealing that iron's penetration depth significantly exceeded that of chromium, approximately reaching 2900 µm.This finding highlights iron's substantially greater diffusion velocity compared to chromium in the β-phase Ti-Fe-Cr alloys.
For species designated as , , and , the corresponding mole fractions are labeled as  ,  , and  .Furthermore, the symbol  denotes the mobility of species  when it is in its pure form, labeled as .The parameters that describe interactions in mixtures with two and three components are respectively notated as  , ( ) for binary interactions and  , , ( ) for ternary interactions.It is pertinent to note that within the framework of binary systems, the equation labeled as Equation ( 14) disregards the ternary interaction term, treating its value effectively as zero.

Composition-Distance Profiles and Diffusion Paths
Figure 2a showcases the microstructural view of the A1 diffusion couple, captured via backscattered electron microscopy, after a 44 h thermal treatment at 1273 K. Figure 2b illustrates the composition-distance profiles for the same couple, obtained through EPMA analysis.All diffusion couples within the Ti-Fe-Cr system are observed in a consistent bcc phase region, with the A1 microstructure serving as a prime example.The EPMA examination along a predetermined direction yielded accurate composition profiles, revealing that iron's penetration depth significantly exceeded that of chromium, approximately reaching 2900 µm.This finding highlights iron's substantially greater diffusion velocity compared to chromium in the β-phase Ti-Fe-Cr alloys.Due to the discreteness of the experimental data points measured using EPMA, errors introduced in the fitting or smoothing process should be avoided.The compositiondistance profiles obtained by EPMA are expressed analytically via using error function expansion (ERFEX) [37,51]: At the specified distance labeled , the mixture's composition is characterized by the parameters  ,  ,  , and  each serving as variables for adjustment.Figure 3 provides a visualization of the composition versus distance profiles for diffusion couples, specifically for A1 and B2, which were maintained at a temperature of 1273 K for a duration of 44 h, and for C2 and D2, subjected to 1373 K for 22 h.These profiles are described utilizing the ERFEX model.In the execution of this study's methodology, the selection of Due to the discreteness of the experimental data points measured using EPMA, errors introduced in the fitting or smoothing process should be avoided.The composition-distance profiles obtained by EPMA are expressed analytically via using error function expansion (ERFEX) [37,51]: At the specified distance labeled z, the mixture's composition is characterized by the parameters a i , b i , c i , and d i each serving as variables for adjustment.Figure 3 provides a visualization of the composition versus distance profiles for diffusion couples, specifically for A1 and B2, which were maintained at a temperature of 1273 K for a duration of 44 h, and for C2 and D2, subjected to 1373 K for 22 h.These profiles are described utilizing the ERFEX model.In the execution of this study's methodology, the selection of i = 4 for the ERFEX model was guided by the detailed examination of the diffusion couples' composition-distance profiles.The employment of the error function in the analysis of diffusion processes offers a framework that is not only methodical but also imbues the findings with a more profound understanding of the physical phenomena being studied.
Materials 2024, 17, x FOR PEER REVIEW 8 of 19 i = 4 for the ERFEX model was guided by the detailed examination of the diffusion couples' composition-distance profiles.The employment of the error function in the analysis of diffusion processes offers a framework that is not only methodical but also imbues the findings with a more profound understanding of the physical phenomena being studied.In the study presented in Figure 4, diffusion paths for various couples were meticulously charted following their exposure to annealing at two distinct temperatures: 1273 K and 1373 K.The resulting diffusion trajectories prominently exhibited an S-shaped curve, indicative of the differential diffusion rates exhibited by Fe and Cr.Close examination of the graphical data revealed a total of 20 intersection points within the diffusion couples annealed at 1273 K, and 15 points at 1373 K. Employing the Whittle-Green method, the derivation of the ternary interdiffusion coefficients at the precise compositions coincide with these intersection points.In addition, the Hall method was used to calculate the impurity diffusion coefficients of chromium in titanium iron alloy and iron in titanium chromium alloy at the end of the A1-A5 and C1-C5 diffusion couples.

Diffusion Coefficients
Tables 2-4 detail the ternary diffusion coefficients in bcc Ti-Fe-Cr alloys, calculated at 1273 K and 1373 K using the Whittle-Green method.The standard deviations were determined from four independent calculations upon two independent measurements.This study further confirms that these coefficients, as derived through the Whittle-Green methodology, are in complete accord with the thermodynamic constraints outlined in reference [45].These constraints are established as follows: In the study presented in Figure 4, diffusion paths for various couples were meticulously charted following their exposure to annealing at two distinct temperatures: 1273 K and 1373 K.The resulting diffusion trajectories prominently exhibited an S-shaped curve, indicative of the differential diffusion rates exhibited by Fe and Cr.Close examination of the graphical data revealed a total of 20 intersection points within the diffusion couples annealed at 1273 K, and 15 points at 1373 K. Employing the Whittle-Green method, the derivation of the ternary interdiffusion coefficients at the precise compositions coincide with these intersection points.In addition, the Hall method was used to calculate the impurity diffusion coefficients of chromium in titanium iron alloy and iron in titanium chromium alloy at the end of the A1-A5 and C1-C5 diffusion couples.
Materials 2024, 17, x FOR PEER REVIEW 8 of 19 i = 4 for the ERFEX model was guided by the detailed examination of the diffusion couples' composition-distance profiles.The employment of the error function in the analysis of diffusion processes offers a framework that is not only methodical but also imbues the findings with a more profound understanding of the physical phenomena being studied.In the study presented in Figure 4, diffusion paths for various couples were meticulously charted following their exposure to annealing at two distinct temperatures: 1273 K and 1373 K.The resulting diffusion trajectories prominently exhibited an S-shaped curve, indicative of the differential diffusion rates exhibited by Fe and Cr.Close examination of the graphical data revealed a total of 20 intersection points within the diffusion couples annealed at 1273 K, and 15 points at 1373 K. Employing the Whittle-Green method, the derivation of the ternary interdiffusion coefficients at the precise compositions coincide with these intersection points.In addition, the Hall method was used to calculate the impurity diffusion coefficients of chromium in titanium iron alloy and iron in titanium chromium alloy at the end of the A1-A5 and C1-C5 diffusion couples.

Diffusion Coefficients
Tables 2-4 detail the ternary diffusion coefficients in bcc Ti-Fe-Cr alloys, calculated at 1273 K and 1373 K using the Whittle-Green method.The standard deviations were determined from four independent calculations upon two independent measurements.This study further confirms that these coefficients, as derived through the Whittle-Green methodology, are in complete accord with the thermodynamic constraints outlined in reference [45].These constraints are established as follows:

Diffusion Coefficients
Tables 2-4 detail the ternary diffusion coefficients in bcc Ti-Fe-Cr alloys, calculated at 1273 K and 1373 K using the Whittle-Green method.The standard deviations were determined from four independent calculations upon two independent measurements.This study further confirms that these coefficients, as derived through the Whittle-Green methodology, are in complete accord with the thermodynamic constraints outlined in reference [45].These constraints are established as follows: 1.678 ± 0.33 In the datasets presented in Tables 2 and 3, the investigation clearly highlights that the predominant interdiffusion coefficients, specifically D Ti FeFe and D Ti CrCr , significantly exceed the values of their cross-interdiffusion counterparts, D Ti FeCr and D Ti CrFe .Notably, the primary coefficients are observed to be larger by one to two orders of magnitude when compared to the cross coefficients.This pattern underscores the dominant role of the intrinsic concentration gradients of Fe and Cr in steering their diffusion behaviors, with minimal influence from the concentration gradients of the other element.Furthermore, the derivations of the cross-interdiffusion coefficients, as facilitated by the Whittle-Green method, encounter substantial variability and oscillations.This outcome manifests a reduced level of precision and highlights the intrinsic variability associated with these calculations.
In Table 2, the value of the main interdiffusion coefficient D Ti FeFe is between 0.79 × 10 −12 m 2 /s and 2.33 × 10 −12 m 2 /s.In Table 3, the value of D Ti FeFe is between 2.71 × 10 −12 m 2 /s and 6.33  5. Comparing the average values in Table 5, the following relationships are obtained: It is shown that the average interdiffusion coefficient and average impurity diffusion coefficient of Fe elements in bcc-Ti-Fe-Cr alloy are higher than those of Cr elements.Then, the average diffusion coefficients of the same element at different temperatures were compared to analyze the effect of increasing temperature on the diffusion behavior of Fe and Cr elements.Specifically, they are as follows: The above equations show that the diffusion coefficient increases by about three times when the temperature increases by 100 K.
In Figures 5 and 6, the dependencies of the primary interdiffusion coefficients, namely D Ti FeFe and D Ti CrCr , on the concentrations of Fe and Cr are depicted for temperatures of 1273 K and 1373 K, respectively.These figures illustrate how the diffusion coefficients adjust as a function of varying elemental compositions at the specified thermal conditions.The different colored markers in the figure represent different diffusion couples.It can be found in Figures 5 and 6 that the main interdiffusion coefficient D Ti FeFe decreases with the increase in Fe and Cr elements.However, the main interdiffusion coefficient, D Ti CrCr , exhibits a completely different phenomenon at 1000 • C and 1100 • C. As shown in Figures 5 and 6 It is shown that the average interdiffusion coefficient and average impurity diffusion coefficient of Fe elements in bcc-Ti-Fe-Cr alloy are higher than those of Cr elements.Then, the average diffusion coefficients of the same element at different temperatures were compared to analyze the effect of increasing temperature on the diffusion behavior of Fe and Cr elements.Specifically, they are as follows: The above equations show that the diffusion coefficient increases by about three times when the temperature increases by 100 K.
In Figures 5 and 6, the dependencies of the primary interdiffusion coefficients, namely  and  , on the concentrations of Fe and Cr are depicted for temperatures of 1273 K and 1373 K, respectively.These figures illustrate how the diffusion coefficients adjust as a function of varying elemental compositions at the specified thermal conditions.The different colored markers in the figure represent different diffusion couples.It can be found in Figures 5 and 6 that the main interdiffusion coefficient  decreases with the increase in Fe and Cr elements.However, the main interdiffusion coefficient,  , exhibits a completely different phenomenon at 1000 °C and 1100 °C.As shown in Figures 5 and 6,  slowly decreases with the increase in Fe and Cr elements at 1000 °C, but increases with the increase in Fe and Cr elements at 1100 °C.

Atomic Mobility Optimization
Employing Thermo-Calc (Educational 2023b) software, which is seamlessly integrated with the DICTRA module, facilitated the refinement of the atomic mobility parameters within the ternary bcc Ti-Fe-Cr alloy system.The results of this optimization process are systematically documented in Table 6, as found in the referenced literature.Based on the thermodynamic parameters of binary system and the ternary diffusion coefficient calculated at 1273 K and 1373 K, the ternary atomic mobility parameters were optimized.

Atomic Mobility Optimization
Employing Thermo-Calc (Educational 2023b) software, which is seamlessly integrated with the DICTRA module, facilitated the refinement of the atomic mobility parameters within the ternary bcc Ti-Fe-Cr alloy system.The results of this optimization process are systematically documented in Table 6, as found in the referenced literature.Based on the thermodynamic parameters of binary system and the ternary diffusion coefficient calculated at 1273 K and 1373 K, the ternary atomic mobility parameters were optimized.
The refined atomic mobility parameters were then applied to calculate the ternary diffusion coefficients for the Ti-Fe-Cr system, ensuring a rigorous comparison with the experimental data at the aforementioned temperatures.A comparison between the calculated and experimental values of the diffusion coefficient is shown in Figure 7.All calculated diffusion coefficients are within acceptable deviations, so the data points are distributed near the x = y line.In the graph, it can be seen that the calculated values of the main diffusion coefficients D Ti FeFe and D Ti CrCr at 1373 K have very small errors compared to the measured values, indicating that the data points are closely aligned with the straight line.The distribution of data points at 1273 K is relatively discrete; except for a few points that deviate from the straight line, most points are distributed on both sides of the line.The above results indicate a strong consistency between the calculated diffusion coefficient values and the experimental measurement values, thus confirming the reliability and accuracy of the optimization results in Table 6.
A diffusion model was developed to mimic the atomic diffusion behaviors under identical initial conditions and thermal treatments.Figures 8 and 9 depict the component distance curves for nine diffusion couple pairs annealed at 1273 K and eight pairs at 1373 K, as determined through experiments and simulations.The figures demonstrate that the simulation curves for the majority of the diffusion couples align closely with the experimental data.However, discrepancies are noted in some simulations, such as for diffusion couples D1 and D3, where the component gradient in the diffusion zone is underestimated, resulting in a simulated diffusion coefficient larger than the experimental findings.This discrepancy might stem from the differential diffusion rates of Fe and Cr, leading to fitting errors during the composition distance curve analysis via electron probe microscopy, as it may not accurately represent the true terminal positions of Fe and Cr diffusion.Figure 10 showcases the comparison of the simulated diffusion paths against the experimental paths for 17 diffusion couple pairs.The substantial concordance between the simulation outcomes and the experimental data underscores the high precision of the derived atomic mobility parameters.The refined atomic mobility parameters were then applied to calculate diffusion coefficients for the Ti-Fe-Cr system, ensuring a rigorous comparis experimental data at the aforementioned temperatures.A comparison betwee lated and experimental values of the diffusion coefficient is shown in Figure lated diffusion coefficients are within acceptable deviations, so the data points uted near the x = y line.In the graph, it can be seen that the calculated values diffusion coefficients  and  at 1373 K have very small errors comp measured values, indicating that the data points are closely aligned with the s The distribution of data points at 1273 K is relatively discrete; except for a few deviate from the straight line, most points are distributed on both sides of th above results indicate a strong consistency between the calculated diffusion values and the experimental measurement values, thus confirming the reliabi curacy of the optimization results in Table 6.

Conclusions
In this investigation, the researchers fabricated two groups of bcc Ti-Fe-Cr diffusion couples, which underwent annealing processes at temperatures of 1273 K and 1373 K for durations of 44 h and 22 h, respectively.This study harnessed the methodologies devised by Whittle-Green and Hall to extract the ternary interdiffusion and impurity diffusion coefficients from the collected data.The atomic mobility parameters of the bcc-Ti-Fe-Cr system were optimized and evaluated using the CALPHAD method.The following conclusions were obtained: (1) Observations from the study elucidate that Fe manifests a higher diffusion velocity than Cr in the bcc Ti-Fe-Cr alloy when subjected to the same thermal conditions.The diffusion rates of both elements in the alloy are positively correlated with the temperature; as the temperature increases, the diffusion rate also increases.(2) At 1273 K and 1373 K, the main diffusion coefficient  decreases with the increase in Fe and Cr elements.The main diffusion coefficient  slowly decreases with the increase in Fe and Cr elements at 1273 K, but increases with the increase in Fe and Cr elements at 1373 K.
(3) Using the atomic mobility parameters, the diffusion process of alloy in diffusion couple was simulated in DICTRA software.A comparative analysis of the simulation outcomes with the experimental data highlighted minor discrepancies, all of which fell within the bounds of acceptable error margins.This congruence attests to the reliability and accuracy of the atomic mobility parameters, validating the comprehensive approach employed in matching the experimental findings with the theoretical models.
-diffusion coefficients at the point where the diffusion paths of two elements crossdetermined by resolving Equation (4a,b).

Figure 2 .
Figure 2. Depictions of (a) the diffusion zone through a SEM backscattered electron image and (b) the composition-distance profiles of the B1 diffusion couple after undergoing 44 h of diffusion at 1273 K.

Figure 2 .
Figure 2. Depictions of (a) the diffusion zone through a SEM backscattered electron image and (b) the composition-distance profiles of the B1 diffusion couple after undergoing 44 h of diffusion at 1273 K.

Figure 3 .
Figure 3. ERFEX depiction of diffusion profiles for (a) the A1-B2 diffusion couple after 44 h at 1273 K and (b) the C2-D2 diffusion couple following 22 h at 1373 K.

Figure 4 .
Figure 4. Diffusion paths of diffusion couples annealed at (a) 1273 K for 44 h and (b) 1373 K for 22 h.

Figure 3 .
Figure 3. ERFEX depiction of diffusion profiles for (a) the A1-B2 diffusion couple after 44 h at 1273 K and (b) the C2-D2 diffusion couple following 22 h at 1373 K.

Figure 3 .
Figure 3. ERFEX depiction of diffusion profiles for (a) the A1-B2 diffusion couple after 44 h at 1273 K and (b) the C2-D2 diffusion couple following 22 h at 1373 K.

Figure 4 .
Figure 4. Diffusion paths of diffusion couples annealed at (a) 1273 K for 44 h and (b) 1373 K for 22 h.

Figure 4 .
Figure 4. Diffusion paths of diffusion couples annealed at (a) 1273 K for 44 h and (b) 1373 K for 22 h.
× 10 −12 m 2 /s.The average values of D Ti FeFe at 1273 K and 1373 K are 1.34 × 10 −12 m 2 /s and 4.89 × 10 −12 m 2 /s, respectively.The average values and ranges of D Ti FeFe , D Ti CrCr , D * Fe(Ti−Cr) , and D * Cr(Ti−Fe)) are shown in Table , D TiCrCr slowly decreases with the increase in Fe and Cr elements at 1000 • C, but increases with the increase in Fe and Cr elements at 1100 • C. Materials 2024, 17, x FOR PEER REVIEW 11 of 19

Figure 5 .
Figure 5.The variations in ternary inter-diffusion coefficients with the compositions: (a)  with Fe, (b)  with Cr, (c)  with Fe, and (d)  with Cr at 1273 K.

Figure 5 .
Figure 5.The variations in ternary inter-diffusion coefficients with the compositions: (a) D Ti FeFe with Fe, (b) D Ti FeFe with Cr, (c) D Ti CrCr with Fe, and (d) D Ti CrCr with Cr at 1273 K.

Figure 6 .
Figure 6.The variations in the ternary inter-diffusion coefficients with the compositions: (a)  with Fe, (b)  with Cr, (c)  with Fe, and (d)  with Cr at 1373 K.

Figure 6 .
Figure 6.The variations in the ternary inter-diffusion coefficients with the compositions: (a) D TiFeFe

Figure 7 .
Figure 7. Analysis of main inter-diffusion coefficients in bcc Ti-Fe-Cr alloys, showcasin ison of logarithmically expressed calculated values against those obtained experiment

Figure 7 .
Figure 7. Analysis of main inter-diffusion coefficients in bcc Ti-Fe-Cr alloys, showcasing a comparison of logarithmically expressed calculated values against those obtained experimentally.

Materials 2024 ,
17, x FOR PEER REVIEW 15 of 19 K, (f) B1 at 1273 K, (g) B2 at 1273 K, (h) B3 at 1273 K and (i) B4 at 1273 K, where the simulated results is illustrated by the solid line.

Figure 10 .
Figure 10.Discrepancy analysis between the modeled and observed diffusion channel values in the Ti-Fe-Cr system post-annealing at (a) 1273 K for 44 h and (b) 1373 K for 22 h.

Table 2 .
Experimental interdiffusion coefficients in bcc Ti-Fe-Cr alloys at 1273 K.

Table 3 .
Experimental interdiffusion coefficients in bcc Ti-Fe-Cr alloys at 1373 K.

Table 4 .
Experimental impurity diffusion coefficients in bcc Ti-Fe-Cr alloys.

Table 5 .
Statistical results of diffusion coefficients obtained experimentally.

Table 6 .
Atomic mobility parameters for the bcc phase of the Ti-Fe-Cr system.

Table 6 .
Atomic mobility parameters for the bcc phase of the Ti-Fe-Cr system.