Precision Thermal Treatments, Atom Probe Characterization, and Modeling to Describe the Fe-Cr Metastable Miscibility Gap

The Fe-Cr metastable miscibility gap has been studied by high-precision thermal treatments, Vickers micro-hardness (HV) measurements, and atom probe tomography (APT). Thermodynamic modeling further supplements the experimental work. The results obtained show that recent thermodynamic descriptions of the metastable miscibility gap found in literature generally overestimates the consolute temperature. We can show that the source of ambiguity in previous studies is most likely a lack of clear distinction between Cr-Cr clustering and α′\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ \alpha^{\prime} $$\end{document} formation. This distinction is here made by APT results, and it leads to a determined consolute temperature of 580 ± 1 °C for Fe0.50Cr0.50. The revised thermodynamic modeling of the metastable miscibility gap captures the experimental results and is consistent with the overall picture from the Fe-Cr data in the literature.


I. INTRODUCTION
TECHNICAL development of Fe-Cr-based alloys such as stainless steel is driven by their application in industry and consumer items. [1,2] The widespread use of Fe-Cr alloys is due to favorable mechanical properties, corrosion resistance, and flexibility of materials design. [2] Their microstructure and thus innate mechanical properties can be broadly divided into four classes: austenitic, ferritic, martensitic, and austenitic-ferritic (duplex). In particular, ferritic Fe-Cr alloys are of high interest in the nuclear industry where they are considered as structural materials for the next generation reactors. [3][4][5][6][7] Their resistance to irradiation-induced swelling and thermal shock is especially attractive in this application. However, Fe-Cr alloys with ferrite or martensite as the main phase, even austenitic stainless steel welds that contain high-temperature d-ferrite, are all susceptible to the ''475°C embrittlement.'' This phenomenon is of critical concern for highly corrosion-resistant ferritic stainless steels at intermediate temperature, which are otherwise known to have good thermal stability. [8,9] The ''475°C embrittlement'' is caused by the formation of a modulated nanostructure comprised of Fe-rich a-ferrite and Cr-rich a 0 -ferrite, which leads to a dramatic loss of toughness. The demixing of Fe and Cr is due to the miscibility gap that is present in the Fe-Cr phase diagram. [10,11] The a toa 0 phase separation mechanisms that govern this demixing phenomenon is either spinodal decomposition or nucleation and growth. [12] In addition to the technical relevance, binary Fe-Cr alloys are favorable for fundamental studies of phase separation kinetics because of the wide miscibility gap, relatively slow decomposition process, and small a toa 0 lattice mismatch.
There is a vast amount of literature available on Fe-Cr, including an extensive review by Xiong et al., [13] that led to the most recent thermodynamic description of the binary Fe-Cr system. [14] In their work, [14] the thermodynamic modeling was supported by atom probe tomography (APT) experiments focused on locating the spinodal line within the miscibility gap on the Fe-rich side to supplement prior experiments. [15] However, recent works [16,17] have indicated that the miscibility gap and in particular the consolute temperature is not conclusively established. In the present work, we, therefore, investigate the upper limit of the miscibility gap by distinguishing between phase separation and Cr-Cr clustering. This is critical to be able to accurately determine the undercooling, which drives the temporal evolution of phase separation. The upper limit of the miscibility gap is located by APT results and Vickers micro-hardness (HV) measurements on samples exposed to high-precision thermal treatments. Thermodynamic modeling of the metastable miscibility gap is performed using the new experimental data and literature data.

A. Materials and Heat Treatments
The investigated samples were fully ferritic Fe-Cr model alloys, supplied by OCAS NV. Three 33 kg ingots were cast in a vacuum induction furnace, a 125 9 125 9 70 mm 3 block was cut, reheated to 1200°C, and then hot-rolled to 15 mm thickness. The realized chemical composition upon delivery is provided in Table I. Before thermal aging, a 10 9 15 9 50 mm 3 bar was cut and homogenized in air at 1100°C for 2 hours followed by water quenching to room temperature. Thermal aging was performed in a Jofra 601 dry-block calibrator furnace, the temperature accuracy of the furnace during heat treatment was ± 0.9°C verified by a K-type thermocouple. Two types of samples were prepared, 0.3 mm sheets that were later cut into comb-shaped samples for atom probe analysis (illustrated in Figure 1) and square samples for hardness measurements (15 9 10 9 3 mm 3 ). After isothermal aging, the samples were quenched in water to room temperature.
The experimental matrix used in this study is provided in Table II. Literature data and the current thermodynamic description were used to select the starting point of heat treatments within the miscibility gap. When the limit was located within a temperature range, the interval was systematically reduced to determine the location of the miscibility gap boundary. Micro-hardness measurements were first used to study the presence of phase separation, [18] then APT was used to verify the demixing.

B. Hardness Measurements
Vickers micro-hardness (HV) was tested on samples after polishing with a final 4000 grit SiC paper. Twelve measurements were made on each sample with a load of 200 g (Fe 0.80 Cr 0.20 , Fe 0.65 Cr 0.35 )-300 g (Fe 0.50 Cr 0.50 ) and a dwell time of 10 seconds. The different loads were used to maintain similar indent sizes to assure equivalent measurement accuracy for all alloys. The reported HV value is the average after removing the lowest and highest two values. The standard deviation was calculated and presented as error bars. The measurements were performed in accordance with the ISO 6507-1 standard. The DHV measurement is the measured relative change of hardness in as-quenched (AQ) state and after isothermal aging (i.e., DHV = HV aged À HV AQ ).

Instrument and reconstruction parameters
The APT measurements were conducted using a CAMECA Ò LEAP 4000HR instrument. The final preparation of the blanks with dimensions 0.3 9 0.3 9 12 mm 3 ( Figure 1) was made using the classic two-step electropolishing method. [19] For this purpose, a microloop was used with acidic/perchloric acid in 75/25-vol pct as the first step and in the second step, a 98/2-vol pct 2-butoxy-ethanol/perchloric acid solution was used. The analysis conditions for the LEAP instrument were as follows: operating temperature 50 K, pulse fraction 20 pct, and detection rate 0.20 pct. Reconstructions were performed in IVAS 3.8, utilizing spatial distribution maps for depth calibration. [20,21] Shank angle, image compression,, and evaporation field were adjusted in accordance with Vurpillot et al. [22] In addition, the material-dependent evaporation field of the Fe-Cr alloys was corrected in IVAS linearly depending on nominal composition and temperature, [23] to account for the atomic volume under the assumption that Vegard's law holds.

Statistical quantification of Cr segregation
The reconstructed volumes were analyzed using frequency distribution in representative subvolumes for direct quantitative comparison of samples. [21] The at. pctCr frequency distribution diagrams were in this work generally taken from 25 9 25 9 75 nm 3 blocks with a fixed bin volume of 50 ions. To quantify the composition amplitude in relation to a theoretically random distribution, the variation parameter (V), as defined by Blavette et al., [24] was used, where P Ã (n) is the observed experimental frequency distribution and BðnÞ is the theoretical binomial frequency distribution. Thus, V in this work is the normalized area difference, which generates a quantitative value 0 £ V £ 2. If V = 0 then P Ã ðnÞ ¼ BðnÞ, and if V = 2, there is no overlap of the P Ã ðnÞ and BðnÞ frequency distributions.
3. The analysis of short-range order The radial distribution function (RDF) using APT data has previously been defined by Marquis [25] and the bulk-normalized RDF can be defined as follows: [26] RDF r ð Þ ¼ 1 q n RDF ðrÞ ð4=3Þpððr þ Dr=2Þ 3 À ðr À Dr=2Þ 3 Þ ½2 where n RDF ðr) is the number of atoms in the shell at distance r around each atom, q is the average ionic density of the volume: ð4=3Þpððr þ Dr=2Þ 3 À ðr À Dr=2Þ 3 Þ is the volume of each shell, and the distance Dr is the thickness of each shell, in this work Dr¼ 0:25nm: The Warren-Cowley (a WC ) parameter is a way of describing the chemical dependency of Cr-Cr shortrange ordering. [27,28] The a WC parameters of APT data have been defined by Ceguerra et al. [29] according to, where pðBjBÞ is the probability of finding B-B neighboring atoms and x B is the nominal mole fraction of element B. The probability of pðBjBÞ is obtained from the first RDF peak as by Moody et al. [30] In this work, the first peak RDFð0Þ is obtained as in Zhou et al. [31] The coordination number of the BCC structure is 8 so the probability of element i and its reference state j becomes, where f is the experimental detection efficiency and C 0 is the nominal concentration, a relative short-range order parameter has previously been used by Ceguerra et al. [32] Thus, the relative DSRO parameter is in this work is defined as, where a j WC is the WC parameter of the nominal composition and a i WC corresponds to the WC parameter within the first RDF shell. Thus, if DSRO = 0 there will be no short-range order present, if DSRO <0 there is negative co-segregation and if DSRO >0 there is positive co-segregation.
Similarly, short-range order has been investigated by APT using the nearest neighbor (NN) distributions by Mukherjee et al. [33] In Table III, the classic variation parameter is given by V as defined by Eq. [1], but the same quantification method has also been used to quantify the variations of the 1st -and 5th-order NN distributions, to analyze short-range interaction, [34,35] with their respective abbreviation, which are V (1st-NN) and V (5th-NN).

Analysis of wavelength
Concentration profiles were extracted from a 120-nm long cylinder with 5 nm radius to maximize the signal-to-noise ratio, with a fixed bin size of 0.6 nm in the analysis direction. The autocorrelation function was used to extract the mean characteristic wavelength k m of the modulated structure by, where R k is the correlation factor, C i is the nominal concentration of element i and C i is the local concentration of element i in a specific binning volume. The selection of bin size is essential as it affects the local composition uncertainty as defined by Danoix et al. [36] Uncertainty becomes significant in small volumes, while over-averaging might hide important features.

Visualization of a 0 morphology
A simple assumption is that the nominal composition (Cr 0 ) is comprised of randomly distributed atoms in the volume. Then during phase separation, Cr atoms migrate out from the bulk to create Cr-enriched a 0 , which results in Cr depletion of the matrix, given the assumption conservation of mass. Therefore, proxigrams over the Cr 0 iso-concentration surface were used to estimate the at. pctCr in a 0 , [37] which is more easily done for clearly defined precipitates as in the Fe 0.80 Cr 0.20 alloy.
The selection of a 0 composition when the amplitude of Cr fluctuations is low or not fully developed means greater uncertainty. Hence, to visualize spinodal decomposition, we utilize the proxigram to plot the Cr sample count over the Cr 0 surface and selected the a 0 threshold as the maximum height of the distribution curve. As an example in the Fe 0.50 Cr 0.50 alloy at 570°C after 120 hours, we obtain a maximum at a position x = 1.113 nm relative to the Cr 0 iso-concentration surface at x = 0.0 seen in Figure 2(a), this position is then equivalent to 58.72 at. pctCr in Figure 2(b) which is how we obtain the iso-concentration threshold used for a 0 visualization.

D. Thermodynamic Modeling
CALculation of PHAse Diagrams (CALPHAD) is a thermodynamic modeling approach with widespread use in materials science, [38] and it is an invaluable tool in materials design. [39] Based on the minimization of Gibbs energy of the system, it is possible to predict its equilibrium state. In the present work, the body-centered cubic (BCC) a phase is described using a substitutional model assuming random mixing between Fe and Cr. The molar Gibbs energy is given by, where o G i is the Gibbs energy of a given constituent and E G m is the excess contribution to Gibbs energy. The excess Gibbs energy is expressed using a Redlich-Kister polynomial, where the interaction parameters, k L, can be expressed as L CrFe ¼ a þ bT. In the present work, a sub-sub regular solution description was used. A drawback of the model is that it does not explicitly take short-range ordering (SRO) into account. Instead, SRO can be taken into account through the excess energy term, [38] and the interaction parameters were here re-optimized to highlight the effect of SRO.

III. RESULTS AND DISCUSSION
A. The Fe 0.65 Cr 0.35 Reference State Investigation Solution-treated Fe 0.65 Cr 0.35 samples were investigated to ensure a homogenous Cr distribution in the as-quenched reference state. Energy-dispersive spectroscopy (EDS) in an electron microscope (Zeiss XB 540) was used first. The EDS analysis consisted of three 1 mm homogeneous line scans on the surface of a HV-prepared sample. Thus, three line scans in the thickness directions were also made which consistently gave a 0.2 at. pct Cr deficit at the surface, and the Cr depleted zone was~60 lm deep. Hence, 6 quantitative spot analyses were made in the center in the thickness direction stepping 5 lm in the thickness direction, and the spot analyses were consistent with a standard deviation of 0.2 at. pctCr.
Subsequently, the solution-treated Fe 0.65 Cr 0.35 alloy was investigated by APT, a resulting 20 9 20 9 60 nm 3 volume of Fe and Cr ions is showed in Figure 3(a). To assess the Cr homogeneity V analysis is used as a measurement. In Figure 3(b), the Cr frequency distributions are presented by the experimentally observed (Cr-Obs.) and the binomial (Cr-Bin.) distributions, and their difference is given as experimental minus random (EMR). The EMR curve reveals that there are random fluctuations of Cr within the solution-treated sample. The result shows that V is highly sensitive to any fluctuations that deviate from the theoretical randomness. In this case, the as-quenched solution-treated state gives V = 0.026, in agreement with previous reports. [18,40] Statistical analysis of the APT results is essential to determine whether the elemental distributions are homogeneous or not. Furthermore, minor Cr segregation can easily be detected in the solution-treated state by the Cr concentration frequency distributions. In this case, both EDS and APT revealed a satisfactory level of homogeneity in the as-quenched reference state. Thus, the solution treatment 2 hours at 1100°C was applied to all alloys, followed by APT homogeneity verification.

B. Hardness Measurements
It is well known that the HV of Fe-Cr-based alloys are greatly affected by the Cr fluctuations present in the ferrite phase. [18,41] The general understanding is that the a 0 formation in these alloys generates coherency strain that restricts dislocation mobility, causing significant loss of toughness. Thus, it is possible to detect phase separation and relate its kinetics to nucleation and growth or spinodal decomposition through the evolution of mechanical properties relative to the as-quenched reference state by DHV.   The DHV of all investigated alloys is presented in Figure 4. If significant hardening develops during the first 120 hours, it indicates that phase separation is present during aging. On the other hand, if the DHV values stay close to zero, it means that the alloys remain close to the as-quenched solution-treated state without significant at. pctCr amplitudes. Hence, the temperature window of the DHV transition region is interesting to further study directly in 3D by APT to determine if phase separation is present.
The different phase separation characteristics of these alloys are largely attributed to the nature of their decomposition mechanism and its thermodynamic properties. Thus, note that the Fe 0.80 Cr 0.20 alloy (indicated in Figure 4 by red circles) is expected to decompose through nucleation and growth at these temperatures, [42] i.e., there is an a 0 activation energy barrier that restricts the onset of phase separation. While Fe 0.65 Cr 0.35 and Fe 0.50 Cr 0.50 alloys (shown Figure 4 by stars and squares) are both known to be unstable, i.e., they decompose continuously without an activation energy barrier. Therefore, it is expected that the DHV transition over the top of the miscibility gap at Fe 0.50 Cr 0.50 is more gradual as the driving force is reduced, whereas, in relation to the nucleation process, in Fe 0.80 Cr 0.20 , DHV exhibits more of a stepwise behavior. The result promotes a more careful APT investigation of the Fe 0.80 Cr 0.20 alloy in the interval from 537°C to 540°C , investigation of the Fe 0.65 Cr 0.35 alloy in the interval from 560°C to 570°C, and an investigation of the Fe 0.50 Cr 0.50 alloy in the interval from 565°C to 580°C.

C. Atom Probe Tomography (APT)
APT allows for 3D investigations of all elements with subnanometer resolution. In this work, mass-resolution was high enough to resolve all isotopes of Fe and Cr except the 54 Cr and 54 Fe isobars, which was not accounted for, the peak at 27 Da was attributed to Fe 2+ . The background noise was in the range of 25 to 50 ppm, which gives a very low detection limit as well. [43] Thus, APT allows for highly accurate investigation of a 0 and short-range Cr fluctuations in the vicinity of the limit of the miscibility gap.
1. Analysis of the Fe 0.80 Cr 0.20 alloy Based on the DHV results, samples of the Fe 0.80 Cr 0.20 alloy were prepared for APT after 120 hours of aging at 537°C and 540°C; the result is shown in Figure 5.
It is evident from Figure 5 that Fe 0.80 Cr 0.20 contains dispersed precipitation of a 0 phase present at 537°C, while at 540°C, there are no such precipitates. However, during the early stages of phase separation, a 0 is not always this obvious. Thus, there is a need for statistical analysis to consistently quantify the strength of Cr-Cr clustering and phase separation in-between samples. Hence, the statistical analysis of the Fe 0.80 Cr 0.20 phase separation, compiled in Table III, serves as a good reference for further comparisons. In addition, in the case of the Fe 0.80 Cr 0.20 alloy, the cluster extraction procedure proposed by Stephenson et al. [44] was used to isolate a 0 precipitates and obtain their composition through radial concentration profiles [45] and the gyration radius has been used to obtain the Guinier radius (r a 0 G ) of a 0 . [46] Figure 5(b) displays a relatively homogenous sample without the presence of any a 0 precipitates, but with a crystallographic low-density Fe-pole, which gives rise to the high-Cr region. It is an experimental artifact that would introduce a sample-biased error, which is not representative of the overall Cr segregation and therefore is excluded in the statistical analysis. Similarly, the less visible pole in Figure 5(a) is also excluded from the analysis presented in Table III. In Figure 5(a), the morphology consists of close to spherical Cr-enriched precipitates. At this point, the selection of at. pctCr threshold is crucial to highlight the morphology qualitatively given the high density of precipitates. Hence, the a 0 threshold selection is made by proxigrams as described in Section II-B-E.
The morphology in Figure 5(a) clearly shows a 0 precipitates which is a result of phase separation while in Figure 5(b), there are no precipitates which is why we believe we have located the limit of the miscibility gap even though the a 0 precipitates do not reach the equilibrium composition. In experiments, equilibrium conditions are rarely met for numerous reasons, e.g., in Fe-Cr-reduced kinetics of transformation is expected as equilibrium is approached, [47] in addition to the experimental uncertainties of precipitate composition measurements inherent by the field evaporation technique. However, these effects do not account for the lack of equilibrium; instead, this deviation is expected due to the precipitate size and the contribution of the surface energy in the early stages; this is previously well known as non-classical nucleation and growth. [48] 2. Analysis of the Fe 0.65 Cr 0.35 alloy During spinodal decomposition, as in the Fe 0.65 Cr 0.35 alloy, the characterization of the a 0 phase morphology is more sensitive to the selected iso-surfaces thresholds. The major difficulty is due to the interconnected structure that develops continuously. The result of spinodal decomposition is more of a gradual a 0 ! a morphological transition in comparison to nucleation and growth, that is when in this case the aging time is restricted to 120 hrs. Even though a wider transition range would suggest a greater undercooling from the MG limit to give greater Cr-Cr interaction, it is seen by the DSRO parameter that Cr-Cr interaction in a 0 in the Fe 0.65 Cr 0.35 is weaker than in the Fe 0.80 Cr 0.20 alloy at equivalent undercooling from the suggested miscibility gap during this transition. Yet, despite the lack of equilibrium, once again, it is still possible to follow the presence of phase separation through the a 0 Cr measurements.
In Figure 6, there is an evolution of the Fe 0.65 Cr 0.35 morphology from 560°C to 568°C where there is a transition into a much more homogeneous state at 570°C . Still, at 570°C, the sample is far from as homogeneous as in the as-quenched state. It should be noted that a small region of nucleation and growth in-between 568°C and 570°C is not ruled out, since the driving force should be very low and could lead to heterogeneous precipitation outside of the probed volume. The statistical analysis of the Fe 0.65 Cr 0.35 alloy is compiled in Table IV.
The V-based analysis reveals a rapid decrease of the Cr composition amplitudes caused by phase separation when approaching 570°C. The Cr composition measurement binning volume was 50 ions or approximately 1.6 nm 3 . There is at the same time as V decreases one observes a rapid increase of the characteristic wavelength, obtained by autocorrelation in the analysis direction. In relation to the solution-treated state, the fluctuations at 570°C may be attributed to its proximity to the miscibility gap even though the fluctuations are not large enough in amplitude to form a 0 . The iso-surface threshold selected according to Section II-B-E. At 540°C a high at. pctCr pillar is visible (also seen at 537°C), which is a characteristic feature of the low hit-density Fe h110i pole.
The DSRO and DNN are local short-range quantification measurements that only consider Cr-Cr interaction on a subnanometer scale. The DNN includes a distance dependency, and thus, it makes it obvious that distribution analysis is dependent on binning volumes. Hence, comparing V and DNN between Fe 0.65 Cr 0.35 and Fe 0.80 Cr 0.20 , we see similar trends and that there is approximately twice the amount of clustering in the alloy Fe 0.65 Cr 0.35 vs. the Fe 0.80 Cr 0.20 alloy for those samples that are considered as being outside of the miscibility gap (i.e., 570°C vs. 540°C). In addition, the DSRO is closely related to the strength of Cr-Cr interactions, within a fixed spherical volume of maximum 0.25 nm. Thus, considering DSRO, we see that the strength of Cr-Cr interaction is equivalent outside of the miscibility gap but much stronger in the Fe 0.80 Cr 0.20 alloy just inside the miscibility gap.
3. Analysis of the Fe 0.50 Cr 0.50 alloy In the Fe 0.50 Cr 0.50 alloy, visualization of a 0 is very sensitive to the threshold selection because of the high content of solute atoms and relatively low Cr amplitudes. The result is that small Cr fluctuations close to the miscibility gap limit can easily be missed or misinterpreted. Hence, a systematic sample-specific approach is necessary, as the one presented in Figure 2. The issue in the Fe 0.50 Cr 0.50 alloy is that a 0 is restricted by the limits of the miscibility gap, which approaches a single point at the consolute temperature. Meaning lower temperatures are required to allow for greater Cr fluctuation amplitudes. Thus, the consolute temperature of the miscibility gap is more difficult to determine than the limits of the lower Cr alloys. The evolution of a 0 morphology of the Fe 0.50 Cr 0.50 alloy can be seen in Figure 7.
In Figures 7(c) and (d), there are still Cr-rich regions present whether these are just the top of the interconnected structure or just fluctuations present due to the proximity of the miscibility gap is not clear only from visual inspection. Therefore, kinetic measurements were made and analyzed by APT. Thus, comparing phase separation after 24 and 120 hours to quantify the evolution of Cr segregation at nanoscale near the indicated limit, the results are presented in Table V. At 565°C and 570°C, one can see a clearly interconnected structure of a 0 (Figures 7(a) and (b)) and in Table V, statistical analysis shows that there is phase separation present. However, the extent of phase separation is rapidly decreased when increasing the temperature from 565°C to 570°C, which is an indication of a decreased driving force for phase separation. Thus, due to the above-mentioned reasons in this case, one needs to rely upon the statistical measurements and the expectation of a reduction in driving force to locate the limit of the miscibility gap. Furthermore, it is expected that there will still be significant clustering present just outside of the miscibility gap and in this case phase separation may be retarded by a sluggish sigma phase formation. [15] From 570°C to 578°C, the V parameter has reduced to less than half and is becoming very close to the same value as for the Fe 0.65 Cr 0.35 alloy at 570°C. Thus, given the DHV indication (Figure 4), the kinetic evolution has been investigated by APT to see how the potential phase separation evolves with time at 578°C and 580°C. In Table V, one can see that the 578°C alloy has a structure that slowly but clearly evolves in V with time (24 vs. 120 hours), this evolution is small but measurable considering the sensitivity of the method in the reference state and taken as an indication of phase separation. On the other hand, there is no evolution of the decomposition at 580°C between 24 and 120 hours, indicating that 580°C is outside of the miscibility gap. In addition, at short range after 120 hours, DSRO and DNN are greater in Fe 0.50 Cr 0.50 at 578°C 120 hours than at 580°C 120 hours which is almost identical to 570°C 120 hours Fe 0.65 Cr 0.35 .

D. Thermodynamic Description
The limit of Fe-Cr phase separation and its temporal evolution is determined by the thermodynamic descriptions of the system. The metastable miscibility gap has been re-optimized using the Thermo-Calc software with a second-order Redlich-Kister interaction parameter in the excess energy term. The new description is based on literature data, compiled in Table VI, and the experimental work presented in this paper, using the thermodynamic model in Section II-D. The thermodynamic Fe-Cr description presented by Andersson and Sundman [49] has been used as it is currently implemented in the Thermo-Calc Software TCFE10 Steels/Fe-alloys database.
Experimentally, Fe-Cr phase separation by a 0 formation is generally explored at temperatures around 500°C , motivated by the rate of transformation. The solubility of Cr in a and a 0 phases determines the limits   of the miscibility gap and thus are used as input for the new re-optimization. Therefore, long-term heat treatments are sought out, and the near-equilibrium Cr content in a 0 is used as miscibility gap estimates. [42,[50][51][52][53] The upper limit of the miscibility gap has been investigated on several occasions, e.g., see Reference 55, 57, 15, but not with the same accuracy as in this work. Still, some of the prior work (e.g., Williams, [57] ) also suggest that the miscibility gap should be lower than the current thermodynamic descriptions. [14,49,58] A reason for these non-conclusive results in literature could be the use of indirect experimental techniques that lacks to distinguish between phase separation and Cr-Cr clustering. Thus, the results obtained in this work by APT in combination with the literature data are presented in Table VI are compiled in Figure 8, presented with the corresponding re-optimization. The re-optimization of the miscibility gap results in an addition of +500 to the regular term 0 L, and a subregular term, 2 L = +2500, to partially account for clustering using the excess energy term. The addition of these parameters improves the description of the consolute temperature of the miscibility gap and reveals a ''Nishizawa horn'' that appears due to the magnetic transition. Still, this re-optimized miscibility gap is by no means final but only shows what can be done using the excess term. In addition, the Cr solubility in a on the Fe side is well established but it could not accurately be accounted for when trying to include SRO. Thus, in a revised thermodynamic description, one should include the new magnetic description as by Xiong et al. [59] and the 3rd-generation unary descriptions of Fe and Cr. Still, an accurate consolute temperature of the miscibility gap is important to determine the undercooling which drives the temporal evolution of phase separation.

IV. CONCLUSIONS
-The limit of a 0 formation is in this work defining the limit of the metastable miscibility gap, and Cr-Cr clustering is separated from phase separation. This definition is supported by the fact that Cr-Cr clustering outside the miscibility gap is not large enough to generate a significant HV impact. The Cr-Cr clustering should be regarded as a different temperature-dependent phenomenon in need of further investigation. -The APT result in this work is in good agreement with a selection of the most cited experimental studies of Fe-Cr alloys (seen in Figure 8). However, the consolute temperature of the miscibility gap of 580°C ± 1°C at Fe 0.50 Cr 0.50 is set slightly higher than obtained by, e.g., Williams. [57] This deviation is probably due to high-precision thermal treatment and state-of-the-art APT characterization used in this work.  [51] Fe 0.60 Cr 0.40 20 SANS Novy et al. [42] Fe 0.80 Cr 0.20 1067 APT Chandra [52] Fe 0.66 Cr 0.24 1738 MB Bergner [53] Fe 0.91 Cr 0.09 Na* SANS Gou et al. [54] Fe 0.75 Cr 0.24 Ni 0.01 10,000 APT Fe 0.70 Cr 0.26 Ni 0. 04 10,000 APT Miller et al. [55] Fe Bal. Cr 0.17,0.19,0.32, 0.45 500 FIM-APT Kuwano [56] Fe 1-X Cr X 50 MB Williams and Paxton [57] Fe 1-X Cr X 1000 HV / RES This Work Fe Bal. Cr 0.20,0.35,0.50 120 APT MB Mo¨ssbauer spectroscopy, SANS small-angle neutron scattering, APT atom probe tomography, HV Vickers hardness, RES resistivity measurements.
*Irradiated samples 0.6 and 1.5 dPa. Fig. 8-CALPHAD re-optimization of the Fe-Cr miscibility gap, using the Thermo-Calc software. [39] -One of the findings in this work is that the phase separation through APT Cr composition distribution analysis confirms that the limit of miscibility gap at Fe 0.50 Cr 0.50 is higher than at Fe 0.65 Cr 0.35 , i.e., the miscibility gap is not flat at the top. A reason as to why the critical temperature of the Fe 0.50 Cr 0.50 alloy might be interpreted as located below the Fe 0.65 Cr 0.35 alloy by DHV is because there is a significant drop of DHV around/just before 570°C in the Fe 0.50 Cr 0.50 alloy, seen in Figure 4. However, in the Fe 0.50 Cr 0.50 alloy, there is a long gradual transition over the top of the MG from the first significant DHV drop which occurs before the first noticeable DHV drop in the Fe 0.65 Cr 0.35 alloy, hence a probable source of confusion in indirect measurements. -Apart from the DSRO analysis, the RDF has been used to investigate the presence of a 2nd maximum on its curve to determine the presence of segregation periodicity in the alloys where a 0 could not clearly be determined. Thus, it exists a 2nd maximum in the Fe 0.50 Cr 0.50 RDF at 565°C, 570°C, and 578°C after 120 hours while it is absent at 580°C. The disappearance of the 2nd RDF peak correlates well with the disappearance of DHV increases as well. This supports the conclusion that 580°C is the upper limit for a 0 at Fe 0.50 Cr 0.50 apart from DSRO and DNN. Using the RDFs, it is possible to see that there are no 2 nd maxima at 570°C Fe 0.65 Cr 0.35 nor at 540°C Fe 0.80 Cr 0.20 , which means there is no repeating periodicity present. Therefore, these temperatures are identified as the upper limit of a 0 formation and thus the limit of the miscibility gap. These RDF curves are included in Appendix A.

ACKNOWLEDGMENTS
The financial support from Centre National de la Recherche Scientifique (CNRS), Region-Normandie, and Carl Tryggers Research Foundation is gratefully acknowledged. This work was partly carried out owing to the experimental GENISIS platform. GENISIS is supported by the region Haute Normandie, the Metropole Rouen Normandie, the CNRS via LABEX EMC, and the French National Research Agency as part of the program ''Investissement d'avenir'' ARN-11-EQPX-0020. This work was also partly performed within the Competence Center Hero-m 2i financed by VINNOVA, the Swedish Governmental Agency for Innovation Systems, Swedish Industry, and KTH Royal Institute of Technology. One of the authors (MS) would like to thank Prof. Wei Xiong for useful discussions.

FUNDING
Open access funding provided by Royal Institute of Technology..

OPEN ACCESS
This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creat ivecommons.org/licenses/by/4.0/.

APPENDIX A: THE NORMALIZED RDF'S FROM 25 3 25 3 25NM VOLUMES
The slope of the A1 Fe 0.80 Cr 0.20 120 hours at 540°C sample is due to minor crystallographic poles that strongly affect Cr segregation but not large enough to be visible in the reconstruction. A2 shows Cr-enriched region with a radius smaller than 4 nm. In A3 Fe 0.50 Cr 0.50 120 hours at 578°C, there is an increase in Cr again from 8 nm which indicates a periodic structure. A4 shows no periodic Cr enrichment at this length scale.