Effect of Stress on Spinodal Decomposition in Binary Alloys: Atomistic Modeling and Atom Probe Tomography

Self-organizing nanostructure evolution through spinodal decomposition is a critical phenomenon determining the properties of many materials. Here, we study the influence of stress on the morphology of the nanostructure in binary alloys using atomistic modeling and atom probe tomography. The atomistic modeling is based on the quasi-particle approach, and it is compared to quantitative three-dimensional (3-D) atom mapping results. It is found that the magnitude of the stress and the crystallographic direction of the applied stress directly affect the development of spinodal decomposition and the nanostructure morphology. The modulated nanostructure of the binary bcc alloy system is quantified by a characteristic wavelength, λ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ \lambda $$\end{document}. From modeling the tensile stress effect on the A-35 at. pct B system, we find that λ001<λ111<λ101<λ112\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ \lambda _{001}< \, \lambda _{111}< \, \lambda _{101}< \, \lambda _{112}$$\end{document} and the same trend are observed in the experimental measurements on an Fe-35 at. pct Cr alloy. Furthermore, the effect of applied compressive and shear stress states differs from the effect of the applied tensile stress regarding morphological anisotropy. 

Self-organizing nanostructure evolution through spinodal decomposition is a critical phenomenon determining the properties of many materials. Here, we study the influence of stress on the morphology of the nanostructure in binary alloys using atomistic modeling and atom probe tomography. The atomistic modeling is based on the quasi-particle approach, and it is compared to quantitative three-dimensional (3-D) atom mapping results. It is found that the magnitude of the stress and the crystallographic direction of the applied stress directly affect the development of spinodal decomposition and the nanostructure morphology. The modulated nanostructure of the binary bcc alloy system is quantified by a characteristic wavelength, k. From modeling the tensile stress effect on the A-35 at. pct B system, we find that k 001 < k 111 < k 101 < k 112 and the same trend are observed in the experimental measurements on an Fe-35 at. pct Cr alloy. Furthermore, the effect of applied compressive and shear stress states differs from the effect of the applied tensile stress regarding morphological anisotropy.
https://doi.org/10.1007/s11661-021-06467-3 Ó The Author(s) 2021 THE role of computational materials design is consistently increasing for the development of novel materials. The design process frequently adopts a hierarchic engineering approach since properties of materials at the macroscale are related to material structures on both the micro-and nanoscale. The bridging of the hierarchic scales is still challenging, for example, connecting first principles calculations to the arrangement of atoms at the nanoscale. However, there are theories for structural evolution that have been proven valid at the nanoscale. One such technique is phase-field crystal (PFC) modeling, which provides a mean to bridge the mentioned gap. [1,2] A more general formulation of the PFC approach is the atomic density function (ADF) theory, which has been constructed by Khachaturyan. [3] The continuum approach of the ADF theory adopts a pseudo-particle approach where atoms are built up by atomic fragments called fratons. This method was developed by Lavrskyi et al. [4] to be able to model complex patterning and treat different types of phase transformations including diffusionless transformations, which involve short-range displacements of atoms.
Spinodal decomposition is a phase separation phenomenon found, for example, in Ag-Cu nanoparticles, Ti-Al-N thin films, Fe-Mn, and Fe-Cr steels. [5][6][7][8][9] This phenomenon is responsible for the self-organization of a nanostructure isostructurally. The phenomenon is beneficial in cutting tools with Ti-Al-N coatings, where the increased hardness gives an increased lifetime of the cutting tools, [10] while, for example, in Fe-Cr alloys, spinodal decomposition restricts dislocation motion and causes severe embrittlement of the alloys. [11,12] The effect of stress is an important consideration in spinodal decomposition. The Ti-Al-N coatings produced by physical vapor deposition have large residual stresses and are also exposed to large stresses during cutting at high temperatures where decomposition occurs. In the case of stainless steel weldments, it is known that the existing residual stresses can enhance the decomposition kinetics. [13] The kinetics and morphology of spinodal decomposition, when exposed to applied tension, have been studied by Cahn. [14] Although Cahn did not perform any simulations, he could deduce from his theoretical treatment that uniaxial tension or compression should result in a directionality of the decomposed morphology that was normal to the axis of the applied tension or parallel to the axis of the applied compression. In a more recent study by Thompson and Voorhees, [15] elastic energy was assumed to be a function of the expansion difference between two coexisting phases constituting different solute concentrations. Their work considered stress in an inhomogeneous system as a function of composition-dependent elastic constants and an inhomogeneous strain field. Principal stress states were investigated to elucidate the relationship between the applied tension and the wavevector of the composition modulations. It was shown that the uniaxial loading both yields directional effects on the developing structure and can alter the limit of instability; that is, the temperature at which the decomposition process starts can be either suppressed or elevated dependent on the elastic properties of the system. The directional effect exhibited alignments of the nanostructure morphology perpendicular or parallel to the applied traction; this effect was found to be coupled to the sign of the applied strain, the solute expansion, and the difference between C 11 -C 12 elastic constants.
Despite the many significant results obtained in the phase-field (PF) modeling of the spinodal decomposition in the Fe-Cr system under external load, the understanding of the influence of atomic rearrangement at the atomic scale in the Fe-Cr alloy under applied stress in the zone of spinodal decomposition is still an open question. In this work, we aim to elucidate the effect of the anisotropic elastic properties of the crystal structure in response to the applied stress through three-dimensional (3-D) simulations and direct experimental comparison. The binary system Fe-Cr is a favorable model system to study the effect of stress on spinodal decomposition due to the small lattice mismatch between Fe and Cr and the relatively slow decomposition kinetics. In addition, Fe-Cr is a technically important alloy used, for example, in reactor pressure vessels in the nuclear industry, [16,17] where thermal embrittlement is of major concern. [18] Furthermore, Fe-Cr alloys are also important structural materials, for example, in the aviation industry, where complex load states may influence the spinodal decomposition and the related fatigue performance. [19] Thus, understanding spinodal decomposition at an atomic scale beyond traditional thermodynamics is important. [20] We apply an extended version of the continuum ADF theory, called the quasiparticle approach (QA), to model the spinodal decomposition of an elastically strained binary alloy at the atomic scale. [4] The benefit of using QA is that we operate at a time and length scale that bridges the gap between DFT and MD simulations to mesoscopic techniques such as PF. This new approach opens a way to answer numerous outstanding questions concerning the atomistic mechanisms of the displacive phase transformation, [21] formation of defects, dislocations, grain boundaries, [22,23] vacancies, [24] and crystallization, the formation of polymers due to the aggregation of monomers in their solution, [25] pattern formation, [26] etc. This list can be significantly extended. In other words, the QA allows description of the physical phenomena at the atomistic scale, keeping the diffusion time scale. Another benefit of the QA is that we do not have to introduce a complex expression for the elastic energy. As in the MD calculations, the interaction potential contains all of the information about the thermodynamic properties of the system and elastic interaction is integrated self-consistently in the model. The result is that we do not assume anything about the phase transformation itself, such as developing composition modulations and material-dependent interaction constants, which have major effects on the developing structure. Instead, we use a simple thermodynamic model that can capture the major effects found in the experiment.
This article is organized as follows. First, an overview of the QA with application to the binary system is presented, with special emphasis on the choice of model parameters. The QA is then applied to the model spinodal decomposition in the binary system under external stress. The effect of the tensile load along the [101], [111], and [112] crystallographic directions is examined. Finally, the simulation results are compared to the 3-D atom probe results. As will be shown in this article, the simulation results are in good agreement with the experimental observations.
To describe the spinodal decomposition dynamics in binary alloys, we use the recently proposed quasi-particle approach (QA). [4] In this approach, the size of the simulation grid is smaller than the distance between neighboring atoms, and we assume that each atom is a sphere comprised of a number of finite elements called fratons. The atomic displacement can be associated with the creation and annihilation of fratons at site r. The choice of model Hamiltonian should describe the interaction of fratons that results in both their ''condensation'' into atomic spheres and the redistribution of these spheres into the desired equilibrium atomic configuration. Then the main variable describing the atomic configurations in the QA is the fraton density probability function, n a (r, t), which is the occupation probability of finding a fraton of the kind a (a = 1,2,…, m) at the site r, where m is the number of components in the system. In other words, n a (r, t) is the probability that a given simulation grid point belongs to the atomic sphere of the kind a. The temporal evolution of the density function of the fratons n a (r, t) is described by the microscopic diffusion equation: [3] dn a *r; t Here, the summation is carried out over all points, r¢, of the computational grid approximating the continuum space. L ab ðr À r 0 Þ is the matrix of the kinetic coefficients (i.e., exchange probabilities per unit time which in this work are: L AA = L BB = 1 and L AB = L BA = -0.5), and F is the free-energy functional. To guarantee the conservation of the number of atoms during the temporal evolution of the system, the condition P *r L ab r * ¼ 0 should be satisfied. The free energy function can be written in the following form: where the first term describes the internal energy with an approximation of pair interaction and the second term represents the entropy term. V ab (r-r') is the potential of the interaction of a pair of fratons, a and b. In the QA, V ab (r-r') is constructed as a sum of two terms corresponding to a short-range and a long-range part of the interaction potential. The short-range interaction potential accounts for the spontaneous formation of atoms and prevents the overlapping of atoms, and the long-range interaction assures the auto assembly of atomic spheres in the desired crystal symmetry. The Fourier transformation (FT) of the interaction potential isṼ where the summation is carried out over all sites of the computational grid, and the wave vector, k, is defined at all quasi-continuum points, k, in the first Brillouin zone of the computational grid, that is, at all points (N 0 ) in the k space permitted by the periodic boundary conditions.
In accordance with Lavrskyi et al., [4] the short-range interaction potential was chosen as the step function where -1 accounts for the strength of the i-i fraton agglomeration into atomic spheres, the constant R is associated with the atomic radius (i.e., at r < R, fratons agglomerate due to the attractive forces), Dr is the width of the contact repulsion part of the short-range interaction between atoms, and f is the height of the repulsive barrier ensuring structural integrity and repulsion between atoms.
For the long-range interaction, an isotropic potential is considered in this work by using a Gaussian distribution function. The number of wells and their position determines the crystal symmetry of the system. In this work, the binary bcc system is considered; it was previously shown that only one minimum potential is needed to reproduce this crystal symmetry. [4] Then the isotropic long-range potential for the bcc crystal can be written as The minimum of the Gaussian function is situated at , where a bcc is the lattice parameter of the bcc crystal.
The elastic properties are essentially defined by the width of the Gaussian r. Thus, the total interaction potential V ð abÞ normalized by the energy density change determines the thermodynamic and mechanical properties of the system, which, in the binary case with a and b fratons, is where k Sr and k Lr are weight constants for short-and long-range interaction potentials, respectively. In these simulations, these constants were set to k Sr ¼ 1 and k Lr ¼ 0:5. [23,27] To assure appropriate selection of the interaction potential parameters, for any given system, the elastic constant C mn can be calculated by C ijkl ¼ @ 2 F=@e ij @e kl .
Strain, e, is introduced by an infinitesimal displacement of the atoms through a deformation matrix, c I is the identity matrix and e * is the strain vector directly defined in real space. Thus, the free energy of a deformed crystal can be evaluated by the change in the interaction potential from V r . [27] Then the elastic constants are calculated through the change in free energy of the system DFðeÞ due to the infinitesimal deformation of the crystal.
For a cubic crystal, C ijkl is reduced to only three independent values: C 11 , C 12 , and C 44 . To calculate these constants, three characteristic deformations are defined in compression/expansion: [24] hydrostatic deformation ðx; y; zÞ Þy; zÞ for C 12 , and monoclinic deformation ( ðx; y; zÞ ! x þ ey; y; z ð Þ ) for C 11 . Thus, the free energy associated with the corresponding deformations is given by where F is the free energy of the unstrained system given by F in Eq. [1]. Usually, to quantify the relation between the elastic constants of a cubic crystal, the Zener anisotropy ratio (A Zn ¼ C 44 = C 11 À C 12 ð Þ ) is used. A Zn ¼ 1 corresponds to a fully isotropic material. It will be shown that the Gauss interaction potential used in our model can be used to reproduce the anisotropic properties of a bcc crystal. To calculate the elastic constants for the bcc a-Fe crystal, the following set of parameters have been used: r = 0.3, k 0 = -0.005, r i = 6.4, Dr = 6.4, and f= 4. The computations are carried out for a 128Dx 3 computational cell with the lattice parameter a 0 = 16Dx and the fraton density q A = 0.125. First, to estimate the free energy of the unconstrained bcc state F 0 , small periodic fluctuations were introduced in the simulation domain and then the system was relaxed to an equilibrium bcc state using the microscopic diffusion Eq. [1]. To evaluate the elastic constants, hydrostatic, orthorhombic, and monoclinic deformations were applied to the equilibrium bcc state by changing the mesh spacing, as mentioned previously. Then the free energy curves as a function of strain were calculated and the values of elastic constants were estimated using the second derivative of the free energy. The obtained results are presented in Table I; included in the table are also ab-initio and experimental data [18,19] for comparison.
The work of Mu¨ller et al. [28] is based on DFT calculations; however, in the article by Mendelev et al., [29] the embedded-atom method potential was used. Hence, there are differences between the predicted and the experimentally measured values of the elastic constants in Table I. Using the ADF approach, proper elastic constants of bcc Fe can be obtained through an atomic interaction defined by either a Gauss potential (as in this work) or by fitting to a structure factor determined by X-ray diffraction, as in Kapikranian et al. [24] The reproduction of Cr elastic constants is cumbersome as Cr elastic is a brittle material with negative Cauchy pressure that requires a negative curvature of the embedding function. Hence, Pasianot et al. [30] only managed to reproduce Cr by a polynomial pair potential, something which was not done in this work due to the complexity of the potential. The anisotropy is of high importance as the influence of crystallographic orientation is the main objective in this work. Thus, it is more important to obtain robust values for Fe C 44 and C', which will dictate the elastic properties of this system.
For infinitesimal displacements where the material exhibits linear elastic behavior, Hooke's law (b r ¼ b Cb e) is generally assumed valid. In this work, the stress tensor b r and strain tensor b e are both field tensors of rank 2 coupled by a material tensor b C of rank 4 characterizing the elastic compliance of the material. Then, it is possible to use the stiffness tensor b S to relate stress, stiffness, and strain (b e ¼ b Sb r). The position of an atom, initially situated at point r, can after deformation be defined using the three-by-three deformation matrix b D. The atom position vector after deformation is given by b r 0 ¼ b Db r. Similarly, in Fourier space, the relation between the vectors k and k', corresponding to the undeformed and deformed structures, can be written using the following expression: In the simulations,kr T ¼k 0r0T = 2p, where k ! is a row vector and k ! T is a colon vector. From this expression, we can relate the displacement of the wave vector in reciprocal space to the deformation matrix in . Thus, the deformation matrix becomesD  [29] ab initio 116 49 178 2.37 In this case, uniaxial deformation in the [001] direction is considered; however, any 3 9 3 deformation matrix can be introduced, and through the introduction of a 3-D rotation matrix, any crystal can be simulated, rotated, and deformed. In this work, the rotation matrix is Hence, the strain is introduced by an infinitesimal displacement of the atoms sitting on the lattice; the displacement is inherently dependent on the discretization of the simulation grid. In this work, the morphology is of importance; thus, a trade-off between computational cost and discretization is made. In this work, the compromise was a 0 = 8Dx.
To relate the simulated and real deformation, the length ratio L/L 0 (current length L over initial length L 0 ) is in this case described by a Green-Lagrange strain tensor, as it accounts for first-and second-order elastic energy contributions given by the mean square displacement. Thus, the displacement of atoms represents strain that leads to a deformation of the simulated volume. The Green-Lagrange description of the deformed simulation volume provides a better scaling to the experiments, and it more easily defines the minimum of the second-order equation, Eq. [7], that defines the elastic limit. The experimental macroscopic deformation and strain were obtained through tensile testing of the Fe-35 at. pct Cr alloy.
As previously stated for the definition of the elastic limit in our simulations, we assume that the linear elastic behavior with the relation DF / e 2 holds. Thus, the limit of elasticity is here defined by the point after which DF starts to oscillate by higher-order terms. It is known that the anisotropic properties of the crystal correlate with anisotropic plastic deformation. [31] Thus, the DF ratio as a function of the uniaxial tensile and shear deformations is obtained by DF e k ð Þ ¼ D mono: e k ð Þ=D otho: e k ð Þ, where e k is strain or, more accurately, displacement introduced in k space; the result is shown in Figure 1. This stress state is related to the Zener anisotropy ratio and, therefore, related to the effective shear. The deformation in our simulations is introduced in reciprocal space; thus, for e k >0, the system is compressed, and if e k <0, the system is elongated (Figure 1). For bcc Fe, the elastic limit is easily defined in compression by the minimum in Figure 1 at e k = 0.61 as it thereafter starts to oscillate. The equivalent volume dilatation in tension is given by e k ¼ À0:213. However, the minima and subsequent deviation from linear elastic behavior in tension begins just after e k ¼ À0:195, which corresponds to 27.4 pct volumetric deformation. This shows a difference in response upon tensile or compression loading of the system. Similar asymmetric stress behavior has been observed by Huang et al., [32] where the difference was attributed to the change in Poisson's ratio. Thus, a discontinuous behavior may be recognized as a change in the crystallography due to deformation. Compression, therefore, is used to define the elastic limit and normalize the load scale as it is well defined and continuous.
The effect of stress on spinodal decomposition was also studied experimentally for a binary Fe-35 at. pct Cr alloy. The model alloy was received in hot-rolled and air-cooled condition. Tensile specimens were produced according to ASTM E-8, and the tensile axis was aligned parallel to the transverse rolling direction. The samples were cut to 2-mm sample thickness using a diamond wire saw. The samples were aged at 525°C for 100 hours under an applied uniaxial tensile load of r = 235 MPa in an MTS Landmark tensile rig. The yield strength of the specimens was r 0:2 = 310 MPa at room temperature. To prepare APT needles along certain crystallographic directions, the gage section of the tensile specimen was analyzed by electron backscatter diffraction after aging. Specific grain orientations were selected, and site-specific APT needles were prepared by FIB. At least three APT needles per orientation were analyzed per condition, and this forms the basis for the results presented here. Statistical analysis was performed on representative subvolumes for comparison. The APT analyses were performed on a CAMECAÒLEAP 4000 HR instrument at 50 K, pulse fraction 20 pct, repetition rate 200 kHz,  and evaporation rate 0.3 pct. The 3-D reconstructions of the tips were performed using IVASÒ3.8 according to the reconstruction scheme by Moody et al. [33] The influence of an applied load on Fe-Cr morphology was previously studied using various simulations and experiments. [7,13,34] In this work, the 3-D experimental atomic-scale comparison is key; thus, we solely work with uniaxial tensile loads to simplify the experiments. The strength measured at room temperature, furthermore, is expected to be reduced at elevated temperatures. To assure that experiments are performed purely within the spinodal regime, experimental aging is made at 525°C, where one can expect~23 pct reduction in yield strength. [35] The simulations are performed at a slightly higher temperature than the experiment; since the simulation kinetics is very sensitive to the selection of temperature, lowering the temperature would dramatically increase the simulation time. Thus, the loads of 235 MPa in experiments and 231 MPa (or e k ¼ À0:171) in simulations are used for the direct comparison. Furthermore, we assume that the Taylor model holds, meaning all grains experience the same level of strain in the polycrystalline alloy. [36] Using the simulations, we start by applying different stress levels in the [001] direction using a pure tensile load. Then selecting the load equivalent to the one used in the experiment, the crystal is rotated in the simulation box so the different principal directions of the bcc crystal are aligned in the z-axis of the simulation box. In addition, the response of the system to different stress states is simulated as well. Figure 2 shows the microstructural evolution of the binary A-35 at. pct B system as a function of tensile load applied along the [001] direction. The minimum strain to trigger a measurable response on the system is in this case equivalent to~19 pct of the elastic limit. The lack of elastic energy contribution at lower strain can be attributed to softening since these simulations are performed at elevated temperatures to favor the kinetics. The critical simulation temperature of instability (T critical Þ is defined in the simulations by the potential k 0 = 0.00; then the undercooling is dependent on the negative depth of the potential and the resulting simulation temperature is easily compared to the critical value. The undercooling in this work is 0.9912 9T critical or close to 565°C in experimental values. This information is used to calibrate the scale for the load level, that is, assuming we have a 19 pct reduction of the elastic limit for the A-35 at. pct B system at 565°C in comparison to room temperature. Thus, with an increased load to~98 pct of the elastic limit of 256 MPa at 565°C, a pronounced alignment of the B-rich region appears. Given the system's response to the external load given in Figure 1, 90 pct of the defined elastic limit, or e k ¼ À0:192, is a more suitable choice for the maximum allowed stress in tension as the abrupt change in behavior may be attributed to a deformation of the crystal.
In addition, Figure 2 shows an alignment of the decomposed structure directly related to the level of the applied load. This alignment naturally affects the wavelength of the structure in the load direction. The morphological effect is expected as the perpendicular directions would contain lower elastic energy during high loads.
Based on the initial results, we expand our simulations to investigate the influence of uniaxial compression along the [001] direction (Figure 3(b)) and shear of the {011} planes (Figure 3(c)) on the morphology; the applied strain was again fixed to 231 MPa equivalent. The simulation results are shown in Figure 3  applied to the [101] direction, a more pronounced alignment of the structure perpendicular to the load is observed. The {101} plane is a close-packed plane; thus, a greater orientation effect in relation to the [001] direction can be expected. However, when the load is applied along the {111} and {112} directions, the alignment of the microstructure starts to rotate in relation to the applied load. To quantitatively characterize this effect, the structure factor s(k) was calculated for the principal directions in their equilibrium state. The result is shown in Figure 4, where s(k) is obtained through the FT of the density probability function n(k,t), also used to describe the primitive unit cell. In Figure 4, we can clearly distinguish the effect of the rotated crystal symmetry and applied load on the peak position. Thus, the first peak of s(k) determines the nearest neighboring atom in the direction of the static concentration wave that is the equilibrium solution to our simulation. As an example, for the 001 orientation, we apply 231 MPa through -0.192 Å displacement of the atomic position in reciprocal space. In Figure 4, the first peak is located at 1. orientations, we can see that their s(k) tends to align themselves in a similar way as the 101 orientation. This gives rise to an approximately 30 deg rotation of the structure in real space, as seen in Figure 3.
In this section, simulating the response to the tensile load was mainly performed in the elastically soft [001] direction for load scale calibration and to potentially maximize the effect of the applied stress. [14,37] This is a known effect on Fe-Cr morphological anisotropy, previously modeled by phase-field methods, [7,38,39] where morphological effects only were observed at very high load levels. The advantage of the QA approach compared to using the phase-field methods is the ability to consider externally applied strains where elastic properties are independent from the developing composition modulations.  Focusing the analysis on the 235 MPa tensile load presented in Figure 5, there is an apparent difference in the volume fraction of a 0 represented by the green atoms for QA and APT. The volume fraction a 0 obtained through the lever rule in QA is roughly~0.40, while in APT, it is closer to~0.1. This difference is due to the construction of the miscibility gap in the model and partially because we run the model to equilibrium, which is not reached in experiments. In APT, a 0 near equilibrium composition consists of~80 at. pct Cr and a contains~14 at. pct Cr, [40] while in QA, we model atoms that are considered pure elements with an associated lattice site occupation probability. It should also be noted that the reason we accept a higher volume fraction in QA is that the simulated volume is much smaller; if the volume fraction is much smaller, we are not be able to quantify the wavelength of the spinodal structure in real space. In Figure 6(b), the QA diffraction pattern is shown as the probability distribution intensity I(k) = IðkÞ ¼ñ k,t ð Þ Áñ k,t ð Þ Ã . In this figure, we can clearly distinguish the fundamental diffraction spots of the bcc lattice and diffuse scattering around some main spots. This pattern is used in comparison with the APT detector hit density to make sure the right crystallographic orientation is analyzed; calculated stereographic projections in the Appendix were used for APT reconstruction calibration.
In Figure 5(a), there is a tendency of a 0 to elongate perpendicular to the applied stress seen experimentally; in that case, the load was applied in the [001] direction, which resulted in a 0 alignment in the [010] direction. The experimental aging time was 100 hours; still, the nanostructural comparison of the results shows good agreement in the morphological alignment between APT and QA. The main morphological discrepancy is a higher degree of anisotropic morphological alignment in QA, presumably attributed to the difference in volume fraction a 0 between the modeling and experiment. Figure 5 also shows the simulated diffraction pattern and detector hit density map of the [001] direction, in which the external load was applied. Naturally, in the experiment, there are a few degrees deviation from the  true calculated projections, which might slightly influence the results. Still, without a perfect APT/QA comparison based on the morphology and behavior dependent on the applied load, we obtain good reproduction of the nanostructure between the experiment and modeling. Note in the comparison that the elastic properties of A and B elements in the A-35 at. pct B system are only fitted to a-Fe for the main component A. Thus, a-Fe is the base element of the simulated system and of the Fe-35 at. pct Cr alloy, which is used for comparison. Hence, no difference in elastic behavior between A and B atoms is modeled, even though that is not true for the real Fe-Cr systems. Figure 5 illustrates the visual correlation between our model and the experimental work; only tensile load is considered for this comparison. High a 0 morphological anisotropy for very high tensile load levels has been seen before by APT in duplex Fe-Cr weld alloys, [13] where the question regarding the thermal effect on the local stress state due to the difference in thermal expansion was unresolved. This issue is not of concern in this analysis as these alloys are fully ferritic binary Fe-Cr alloys with conservative load levels. The conservative load level is presumably the reason no significant morphological anisotropy could be quantified in these experiments. However, there is good agreement in the system's response to the applied load regarding its effect on the characteristic wavelength. Experimental investigations in this case show accelerated fluctuations in Cr composition amplitudes dependent on crystallographic orientation, which can be seen in Figure 6.
The bulk normalized concentration at zero distance is evaluated by the radial distribution function (RDF), which, in this case, relates to the amplitude of Cr fluctuations. [41] In Figure 6, one can clearly evaluate [111] < [001] < [112] < [101] even though the differences are small. The Cr amplitude is known to be related to the degradation of mechanical properties. [11] The difference in morphology of the nanostructure after applying tensile load along the [001], [101], [111], and [112] directions is very small in both APT and QA. However, in the trend of the characteristic wavelength k obtained through autocorrelation in real space, we find from the modeling that k 001 <k 111 <k 101 <k 112 , and in the experiments, we obtain the same results from RDF: k 001 < k 111 <k 101 <k 112 . These results are presented in Figure 7.
In APT, the k 001 and k 111 are relatively similar, while the close-packed slip plane of [101] and the [112] slip plane both produce longer wavelengths k 101 and k 112 , respectively. A similar trend is observed by QA. However, in QA, the supposedly soft orientation [001] produced a much shorter wavelength, even though we expect the opposite given s(k) in Figure 4. Still the results in Figure 7(a) confirm the results shown in Figures 3 and 5. Thus, this result may be attributed to a visualization effect in 3-D and the subsequent voxelization of real space to calculate the wavelength k by autocorrelation. By visual inspection of the QA volumes in Figures 3(b) and 2, one would expect the [001] wavelengths to be longer than the others. Thus, one should be careful to draw definitive conclusions given the relatively small difference even though the trend seems obvious.
In phase-field modeling of phase separation in anisotropic elastic bodies, anisotropic alignment of the a 0 phase is derived from the Cahn-Hilliard theory. [38,39] Zhou et al. [7] attributed the alignment of a 0 to the ratio of the shear modules between the precipitates and matrix, given that a 0 precipitates in Fe-Cr alloys are soft. A soft precipitate will align itself perpendicular to the applied tensile strain. Li et al. [34] predicted directional growth through effective eigenstrain; a direction with lower eigenstrain is generally a favorable direction for decomposition to proceed. In the case of Fe-Cr alloys, e xx < 0 and e yy > 0; thus, eigenstrain is compressive in the x-axis and tensile in the y-axis. Therefore, elongation occurs in the y-axis perpendicular to an applied load in the z-direction. This is in line with our observations in QA and APT ( Figure 5); the merit of using QA is to maintain true atomic resolution and crystal structure in the analysis of nanostructural segregation on the continuum time scale without assuming concentration-dependent elastic properties.
Treatment of anisotropy of drift diffusion due to elastic strain fields in cubic crystals requires consideration when comparing results from different orientations. [42] In the presence of an external force field, the diffusion current density needs to consider drift current, which is proportional to the gradient of the interaction energy. The symmetry of the host lattice can strongly affect the symmetry of the elementary-jump mechanism of the defect (e.g., vacancy) itself. The addition of a continuous strain field alters the saddle point in the activation energy barrier and jump rates in different directions. In addition, the diffusion constant itself is a hyperbolic strain-dependent coefficient. [43] The migration energy of diffusion in a force field is proportional to the gradient energy between the equilibrium energy and the saddle point; the external field affects the saddle point and the equilibrium energy remains unchanged. Under normal conditions, the external strain is small so that it can be linearly approximated. Thus, the expected effect of uniaxial loading in different crystallographic directions on the diffusion constant is only a few percent.
1. The effect of applied strain can be modeled at the atomic scale using the QA approach. For the first time, the real strain in the system was quantified and directly compared with experimental data. Thus, different strain tensors and rotation matrices are introduced to alter the stress states. This approach has shown promising results for tensile, compression, and shear stress in addition to the rotation of the crystal structure. Introducing an infinitesimal displacement of an atom triggers a free energy response of the system to alter its morphology in good agreement with APT. In the case of Fe-Cr, high-load levels are required to have a noticeable morphological anisotropy effect near the elastic limit. 2. The added external load does affect the diffusion-controlled decomposition process in the experiment. The addition of elastic load in the Fe-35 at.
pct Cr system only has minor effects dependent on the differences in the crystallographic orientation. Still, the close-packed [101] orientation favors decomposition in these experiments, indicating that the ease of diffusion is more important than the difference in elastic properties in the direction of the applied load between the different orientations.
3. The exact stress state in the polycrystalline materials is still unknown. However, it is unlikely to have a dramatic effect on these results. It should be mentioned that there is still room for improvement in the treatment of the polycrystalline materials.

ACKNOWLEDGMENTS
The financial support from the Centre National de la Recherche Scientifique (CNRS), Region-Normandie, EIT Raw Materials project ENDUREIT, and Carl Tryggers Research Foundation is gratefully acknowledged. This work was partly carried out owing to the experimental GENESIS platform. GENESIS 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. The authors thank Drs. A Vaugeois, C. Keller, and A. Guillet for their support and discussion of this work.

CONFLICT OF INTEREST
The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this article.

DATA AVAILABILITY
All datasets, code, and results occupy several gigabytes of data, which are available from the authors upon serious request.

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/.