Arxiv:1902.01645V1 [Cond-Mat.Mtrl-Sci] 5 Feb 2019 Magnetic Phases of Fe-Cr Alloys Are Fairly Complex

The Dynamics of Magnetism in Fe-Cr Alloys with Cr Clustering

Jacob B. J. Chapman,∗ Pui-Wai Ma, and S. L. Dudarev UK Atomic Energy Authority, Culham Science Centre, Oxfordshire, OX14 3DB, United Kingdom

The dynamics of magnetic moments in iron-chromium alloys with different levels of Cr clustering show unusual features resulting from the fact that even in a perfect body-centred cubic structure, magnetic moments experience geometric magnetic frustration resembling that of a spin glass. Due to the long range exchange coupling and configuration randomness, magnetic moments of Cr solutes remain non-collinear at all temperatures. To characterise magnetic properties of Fe-Cr alloys, we explore the temperature dependence of magnetisation, susceptibility, Curie temperature and spin- spin correlations with spatial resolution. The static and dynamic magnetic properties are correlated with the microstructure of Fe-Cr, where magnetisation and susceptibility are determined by the size of Cr precipitates at nominal Cr concentrations. The Curie temperature is always maximised when the solute concentration of Cr in the α phase is close to 5 to 6 at.%, and the susceptibility of Fe atoms is always enhanced at the boundary between a precipitate and solid solution. Interaction between Cr and Fe stimulates magnetic disorder, lowering the effective Curie temperature. Dynamic simulation of evolution of magnetic correlations shows that the spin-spin relaxation time in Fe-Cr alloys is in the 20 to 40 ps range.

I. INTRODUCTION and 84 at.% Cr, Fe magnetic moments prevent long range AFM-like correlations, causing the moments to freeze 15 Iron and chromium are the main components of steels. into a spin glass state . For higher concentrations of Fe, Fe-Cr systems attain FM ordering. The Curie tempera- They are both magnetic transition metals. Fe-Cr alloys exp adopt various structural and magnetic phases depending ture TC is found to be maximised at TC ≈ 1049K with 16 on the chemical composition and temperature. Theoret- Cr concentration near 6 at.% . An ordered structure of ical understanding of the correlation between magnetic the alloy at 6.25 at.% Cr is stabilised by a minimum of 17 properties and microstructure provides valuable insight the spin-down density of states at the Fermi energy . for tuning and characterising the mechanical properties Fe-Cr alloys generally adopt a BCC structure across of Fe-Cr alloys and steels. Understanding the coupled the whole chemical composition range, except at low Cr mechanical and magnetic properties that evolve as the concentration and at high temperature, where they ex- underlying microstructure changes is particularly signifi- hibit an FCC γ-loop in the phase diagram18. Despite cant for steels development, including reduced activation having similar lattice constants, where aF e=2.87A˚ and ferritic martensitic (RAFM) steels for advanced fission aCr=2.91A,˚ they do not always form a solid solution. and fusion reactors1. Between room temperature and 700K, a miscibility gap Magnetic ordering strongly affects structural and me- is observed starting from 5-10 at.% Cr and terminat- chanical properties of Fe-Cr alloys. For example, the to- ing at 90-95 at.% Cr. The precise boundaries remain 18–20 tal free energy resulting from both magnon and phonon contested . Within the miscibility gap, there is a do- excitations stabilises the ferromagnetic (FM) body- main of spinodal decomposition. At temperatures above centred cubic (BCC) phase of Fe in comparison with the 1100K, iron and chromium are fully miscible. Fe-Cr al- hexagonal close-packed phase structure stabilised by the loys exhibit anomalous formation enthalpy as a function 21,22 s and d orbital occupation alone2,3. Self-interstitial atom of chemical composition and temperature . At low (SIA) defects in all the non-magnetic transition metals temperature, when Cr concentration is below approxi- form h111i crowdion configurations4–6. On the other mately 10 at.%, short-range ordering of the Cr magnetic hand, SIAs in magnetic Fe and Fe-Cr alloys prefer h110i moments gives rise to a negative formation enthalpy, dumbbell configurations with antiferromagnetic (AFM) which favours complete dissolution of Cr. ordering at the core of the defect7,8. On the macroscopic The 475◦C embrittlement in steels23 is a well known scale, the softening of elastic constant C0 near the Curie phenomenon associated with spinodal decomposition. 9,10 temperature TC results from magnetic excitations . Thermally aged Fe-Cr alloys undergo slow coherent phase arXiv:1902.01645v1 [cond-mat.mtrl-sci] 5 Feb 2019 Magnetic phases of Fe-Cr alloys are fairly complex. separation into Cr-rich α0 and Fe-rich α precipitates if The magnetic ground state of pure BCC Cr is an in- the chemical composition is in the miscibility gap23. Ion commensurate spin density wave (SDW) with the Néel irradiation experiments24 have shown that irradiation- exp 11 accelerated Cr precipitation is closely related to the hard- temperature TN ≈ 310K . Magnetic moments are in antiparallel alignment with magnitudes varying si- ening of Fe-Cr alloys. The first wall of fusion reactors is nusoidally and being incommensurate with the lattice expected to sustain high irradiation doses approaching periodicity11–13. The introduction of isolated Fe atoms 10 to 20 displacements per atom (dpa) per full operation into the Cr matrix suppress the ordering by destabilising year1,25. the correlation between electron-hole pairs, which makes In general, precipitation during thermal ageing occurs the SDW commensurate at 96-98 at.% Cr14. Between 81 across a time scale exceeding 104 hours26. Under irradi- 2 ation, this can be accelerated by 6 to 7 orders of mag- both the direction and magnitude of an magnetic moment nitude due to a substantial increase in point defect con- (or atomic spin) as a three dimensional vector3,40. Trans- centrations enabling fast mass diffusion27. Despite dif- verse and longitudinal fluctuations are treated within a ferent time scales, the growth morphology of Cr clusters unified framework, corresponding to the localised and is indistinguishable between irradiation and annealing27, itinerant property of hybridised s and d electrons41. with the exception in the vicinity of sinks, such as cav- The temperature of a dynamic spin system can be con- ities and grain boundaries, where precipitate free zones trolled by Langevin thermostat40,42–44. Considering an form28,29. arbitrary spin Hamiltonian H, the Langevin equation of 40 The kinetics of Cr clustering in Fe-Cr alloys has been motion of an atomic spin Si can be written as : studied extensively through experimental26,30–34 and the- 27,34,35 dSi 1 oretical approaches . However, the effect of Cr pre- = Si × Hi + γsHi + ξi (1) cipitation on the magnetic properties of Fe-Cr alloys has dt ¯h largely been ignored from quantitative analysis. This mo- where the Hi = −∂H/∂Si is the effective field for spin i, tivates a systematic investigation on the correlation be- γs is a damping parameter and ξi is δ-correlated random 0 tween the magnetic properties of Fe-Cr alloys and their noise, which satisfies hξi(t)i = 0 and hξiα(t)ξjβ(t )i = 0 microstructure. µsδijδαβδ(t − t ). Greek alphabet subscripts repre- In this paper we quantitatively explore the relation- sent the Cartesian coordinates. According to the ship of magnetic properties of Fe-Cr alloys with the un- fluctuation-dissipation theorem (FDT)45,46, one can ob- derlying atomic ordering. We show that the Curie tem- tain a relationship between the fluctuation and dissi- 40,47 perature is not only dependent upon the nominal con- pative forces , namely µs = 2γskBT . The mag- centration of solute atoms but also on the degree of Cr netic moment and atomic spin are simply related by clustering. Furthermore, using correlation functions we Mi = −gµBSi, where g is the electronic g-factor and show that the local Curie and Néel temperatures change µB is the Bohr magneton. in the vicinity of the cluster interfaces. We shall under- We perform SD simulations using our code stand these phenomena from the scientific perspective, SPILADY48. The source code has been generalised showing that some of the peculiar magnetic properties to handle multiple elemental species and complex stem from geometric magnetic frustration. The random- Hamiltonian forms. In all the calculations, the equations ness of exchange coupling between magnetic moments in of motion of spins are integrated using the parallelised Fe-Cr alloys essentially form spin glasses. Knowing the symplectic Suzuki-Trotter decomposition algorithm49 relationship between magnetic properties and precipita- with a time step of 1 fs. tion in a quantitative manner may help in the development of a non-invasive engineering technique to estimate the level of Cr precipitation by probing local magnetic B. Magnetic Cluster Expansion Hamiltonian properties. Our paper is structured as follows. In Section II, the We adopt a Heisenberg-Landau type magnetic cluster method of Langevin spin dynamics, the Hamiltonian and expansion (MCE) Hamiltonian developed for Fe1−xCrx the quantitative properties of interest are introduced. In alloys, which has been used in various Monte Carlo (MC) Section III A, magnetisation of solid solutions of Fe-Cr simulations50–53. Terms in the MCE Hamiltonian, con- alloys are investigated, which are compared to alloys con- tributing to the effective field, can roughly be decom- sisting of ideal and impure Cr clusters in Sections III B posed into inter-site interaction Heisenberg terms HH and III D, respectively. We examine the change of Curie and on-site Landau terms HL: temperature TC as a function of the size and chemical composition of clusters. The dynamics of magnetisation H = HH + HL, (2) characterised using correlation functions of the idealised where clusters is examined in Section III C and for (001) Fe/Cr interfaces in Section III E. Our results are summarised in X (0) (1) (2) Section IV. HH = Jij +Jij (σi +σj)+Jij σiσj Mi ·Mj, (3) i6=j and II. METHODOLOGY X (0) (1) X (2) 2 HL = A + A σi + Aij σiσj Mi A. Spin Dynamics i j6=i X X (2) + B(0) + B(1)σ + B σ σ M4. (4) Spin Dynamics (SD) is an indispensable tool for study- i ij i j i ing magnetic excitations. Originating from the Landau- i j6=i Lifshitz and Gilbert equations36,37, SD has vastly ex- The magnetic energy of an alloy configuration depends tended our understanding of solitons and magnons in on atomic site occupational variables σi and the non- 3,38,39 magnetic media . Recent developments consider collinear magnetic moments of each atom Mi. The 3 discrete occupational site variable σi takes a value of the trajectory in a time period of τ, assuming ergodicity: th σi = +1 if the i atomic site is occupied by an Fe atom,   or σi = −1 if it is an Cr atom. Nη 1 Z τ 1 X Since the magnetic interactions depend on the atomic hMηi =  Mi(t) dt. (5) 54 τ N configuration, both the magnetic order and magnitudes 0 η i evolve according to the local atomic environment via the Landau coefficients Ai, Bi and the Heisenberg exchange The summation spans all atoms Nη of the species. The coupling parameter Jij. The Jij are parameterised up to averaged system magnetic moment is the contribution th (2) (2) the 5 nearest neighbour. Aij and Bij are functions from all element types. Magnetisation is the magnitude st nd of atomic species in the 1 and 2 coordination shells. of the averaged magnetic moment i.e. M¯ η = |hMηi|. The numerical values of the parameters for the Fe1−xCrx The susceptibility of an observable with respect to its Hamiltonian in Eq. 2 to 4 are given in Table 1 of Sup- conjugate field can be calculated from the fluctuations plementary Material. Details on the parametrisation are of the observable within a canonical ensemble. This fa- given in Refs. 50 and 52. cilitates the calculation of the magnetic susceptibility χ Magnetocrystalline anisotropic energy is not included through the spin dynamics calculations: in the MCE Hamiltonian, but can in principle be in- corporated phenomenologically55. It would require ex- hM2i − hMi2 tra fitting of parameters from ab initio calculations that χ = µ0 , (6) kBT consider spin-orbit coupling. As we are working on perfect bcc configurations, it is expected that the anisotropic 56,57 where µ0 is the vacuum permeability and M is the aver- energy is in the µeV/atom level . As we will see in aged magnetic moment within a defined volume. the following sections, the origin of magnetic frustration One can estimate the fundamental timescales of mag- is due to the long range exchange coupling which is in netic excitations using the correlation of fluctuating mag- meV level. We may safely ignore the magnetocrystalline netic moments. The correlation of a magnetic moment anisotropy effects in this study. at site i over a time interval of τ from a starting time t We verified our implementation by performing time av- can be calculated using a time-displaced spin-spin auto- eraged benchmarking comparisons for pure Fe, Fe-Cr al- correlation function. loys and (001) Fe/Cr superlattices. They show perfect agreement with MC calculations performed by Lavren- N 50–52 1 X tiev et al. . The results and details can be found in c (r, τ) = e (r, t + τ) · e (r, t) , (7) 0 N i i Supplementary Materials. n A fundamental difference between Langevin SD (Eq. 1) and MC methods is the process of attaining equilib- where ei = Mi/Mi is the unit vector of magnetic mo- rium. Whilst both approaches should give the same ther- ment. The autocorrelation function is calculated on-the- modynamically averaged quantities in the limit of time fly59 at thermal equilibrium and averaged using multiple going to infinity, dynamics simulations show if a system statistically independent starting times. can really achieve equilibrium within realistic timescales. The time-displaced correlated behaviour between the MC methods, on the other hand, do not contain an ob- ith magnetic moment and the magnetic moments con- jective concept of time. It provides no information on stituting its 1st coordination shell give insights into the the dynamics or time relaxation of events. magnetic ordering and timescale of its decoherence: The current model is also applicable to the study of spin glasses. Fe-Cr adopts a spin glass magnetic state N 8 at high Cr concentrations15. A common theoretical ap- 1 X X 0 c1(r, τ) = ej(r , t + τ) · ei(r, t) , (8) proach is to consider a phenomenological Ising model. In 8N n j=1 these cases the Heisenberg parameters are drawn from a Gaussian distribution about a non-zero mean58. On where atom j is in the 1st nearest neighbour shell of atom the other hand, the MCE Hamiltonian (Eq. 2-4) encodes i. At cryogenic temperatures, a FM monodomain would configurational disorder using physically meaningful ex- yield c = 1 whereas a layered AFM domain such as in change parameters that were fitted to ab initio data52. 1 L10 phase FeNi would return c1 = −1. In the following, when we consider a spherical precipitate, we measure the averaged correlation functions of C. Simulation Metrics magnetic moments displaced from a chosen origin r0, i.e. the geometric centre of a cluster. We call c0(|r − r0|) An important macroscopic order parameter is the time the radially resolved time-displaced spin-spin autocorre- st averaged magnetic moment. The contribution by species lation functions (RRACF), and c1(|r−r0|) the 1 nearest η = Fe or Cr is calculated using the magnetic moments in neighbour radially resolved time-displaced spin-spin cor- thermal equilibrium at a particular temperature T over relation functions (1nRRCF). 4

III. RESULTS AND DISCUSSION

A. Disordered Solid Solutions

Our investigation starts from disordered solid solutions of Fe1−xCrx, where 0 ≤ x ≤ 0.25. We study the change of Curie temeprature TC through the change of magnetic susceptibility χ as functions of Cr concentration and temperature. SD simulations were performed using a cubic simulation cell with periodic boundary conditions containing 31250 atoms. Collinear ordering of magnetic moments along [001] was chosen for the initial configuration. We define the magnetisation axis parallel to the z direction. Statistics were collected over a 5ps production run following a 1ps thermal equilibration period. Each configuration was generated by random substitution of Cr atoms into the BCC Fe matrix to satisfy the concentration requirements. As we use a relatively large simulation cell, we assume sufficient statistical sampling of local structure. To check this assumption, for x = 15 and 25 at.%, we simulated three random configurations and found comparable results for each case. The magnetisation of the disordered solid solutions are shown in Fig. 1(a) over a temperature range of 100K ≤ T ≤ 1500K. At low temperatures, the magnetisation decreases as the Cr concentration increases. The magnetisation profiles are similar to pure Fe as the temperature is raised. For Cr concentrations of 5 and 9 at.%, the magnetisation is found to remain finite to temperatures higher than the TC of pure Fe. We determine TC through the turning point in the bulk 2 susceptibility (∂T χ = 0) as shown in Fig. 1(b). TC is plotted against the Cr concentration in Fig. 1(c). We note the consolute temperature is approximately 900K. Therefore, disordered solid solutions are appropriate sta- ble configurations near the TC for aged Fe-Cr alloys across the full compositional range. Data from experi- FIG. 1. Magnetic properties of disordered Fe1−xCrx solid solutions with different Cr concentrations. (a) Magnetisa- ments and MC simulations16 are also shown for compar- 16 tion versus temperature. (b) Magnetic susceptibility versus ison. In the MC simulations, Lavrentiev et al. deter- temperature. (Inset: Susceptibility around the Curie temper- mined the TC through the observation of peaks in the ature TC .) The points of inflection for x = 0, 5 and 25 at.% specific heat instead. Excellent agreement is found be- Cr are indicated by arrows. (c) Curie temperature TC versus tween our SD calculations and the MC results, where Cr concentration as determined from the susceptibility. Data both show the same trend as experiments. from Monte Carlo simulations and experiment16 are given for TC increases for small Cr concentrations up to 5 comparison. The determination of TC presented in Ref. 16 at.%Cr. Further increases in the Cr concentration de- is via observation of peak in the specific heat. Error bars are set as the temperature interval in sampling, which is 20K. crease the TC . This is in agreement with experimental observations that TC is maximised for a Cr concentration around 6 at.%16,60. The physical mechanism for this phenomenon is known to occur due to the strengthening of the Fe magnetic moment induced by the strong antiferromagnetic coupling to Cr atoms on nearest-neighbour lattice sites16,21,52. If an Fe atom has multiple Cr near- support perfect AFM ordering of all Cr magnetic mo- est neighbours, it results in magnetic frustration, because ments with respect to their neighbouring Fe atoms. This antiparallel ordering with respect to the Fe atom cannot is achieved when the Cr atoms take the corner sites of a be held for all Cr neighbouring moments simultaneously. repeating unit constructed of 2 × 2 × 2 unit cells. This is A concentration of x = 6.25 at.% Cr is an idealised reflected with a pseudo-gap minimum in the spin-resolved case. It allows for the highest concentration of Cr which density of states17. 5

B. Fe-Cr Alloys with Idealised Cr Clusters

Irradiation alters the microstructure of steels, which influences their mechanical properties24. We expect different microstructures to exhibit different magnetic properties, even when the nominal chemical composition of the overall system remains the same. If verified it sug- gests magnetic properties may be probed externally with- out requiring the destructive extraction of test samples from components. This may help develop a non-invasive means of assessing microstructure in steels. Towards this end, in this section we investigate how the formation of pure Cr α0 precipitates might change the magnetic properties with respect to the disordered solid solutions studied in previous section. Atom probe experiments show precipitates form during ageing and irradiation. They starts forming as elongated coagulated clusters and becomes roughly spherical over time26. We model spherical α0 precipitates by substitut- ing Fe atoms with Cr atoms within a sphere of radius rc = |r − r0|. Simulation cells containing 31,250 atoms in BCC structures are created. To achieve a particular nominal composition for the Fe-Cr alloy, we distribute extra Cr substitutions randomly throughout the remaining Fe matrix. The system size corresponds to an α0 precipitate number density of 2.7×1023 m−3 which is comparable to the lower limits observed experimentally in neutron irradi- FIG. 2. The averaged magnetisation of Fe Cr alloys is 61 1−x x ated Fe-Cr steels . Clusters are generated with radii plotted as a function of temperature and cluster size with from 6 to 26A˚ for nominal concentrations of 9, 12, 15, different chemical composition. In each simulation cell, we 18, 22 and 25 at.%Cr. The cluster sizes are chosen to created a pure Cr α0 spherical precipitate. Extra Cr atoms is be comparable to precipitates aged from 10 to 10,000 randomly distributed in Fe matrix to attain the required over- 26,61 hours . all Cr concentration x. The radius of precipitates rc ranges ˚ ˚ The averaged magnetisation for Cr concentrations of from 6 to 26A. Disordered solid solution is denoted as rc = 0A. 9% and 25%, with different precipitate sizes and at dif- Only x = 9% and 25% are shown here. Figures of x =12%, ferent temperatures, are shown in Fig. 2. Figures re- 15%, 18% and 22% are presented in Supplementary Material. garding Cr concentrations of 12%, 15%, 18% and 22% are presented in Supplementary Material for complete- atomic configurations with clustering reduces the over- ness. The magnetisation of bulk Fe (|hMFei|), the scaled bulk value proportional to the nominal concentration all frustration in the system. This reinforces collinearity and enables the constructive superposition of the local ((1 − x)|hMFei|) and the magnetisation of the disordered moments. solid solution (rc = 0) are also plotted. The scaled bulk Fe curve gives an approximate limit obtained upon In Fig. 3(a)-(c) we present a snapshot of magnetic complete phase segregation as could theoretically occur moments of a system with a precipitate of radius 18A.˚ within the spinodal decomposition. The Fe-9Cr system has 31,250 atoms. Within the α0 In general, the disordered solid solutions have the precipitate, Cr moments are antiferromagnetic and per- smallest magnetisation. Both the Fe and Cr moments pendicular to Fe moments in the α phase. Cr moments prefer being antiparallel to neighbouring Cr moments. situated on the cluster interface are aligned antiparallel With disordered configurations (when x > 0.06) this can to the adjacent Fe moments. They gradually rotate to be rarely be satisfied. In order to minimize the energy, Cr perpendicular with the Fe moments when they are well moments try to be as antiparallel as they can with respect separated from the interface by 5-8A.˚ Throughout this to their neighbours, so they tilt. This results in geometric tilting from the Fr moments, the Cr magnetic moments magnetic frustration and non-collinearity of the Cr mag- attempt to maintain an antiparallel ordering to their Cr netic moments (see supporting Figure 2(e) in the Sup- neighbours. plementary Material). Such ordering reduces the total An alternative representation is the magnetic moment magnetisation. On the other hand, for the same nominal space as shown in Fig. 3(d) and (e). Each dot represents concentration x, the averaged magnetisation increases as a magnetic moment vector, where all moment vectors Cr is grouped into larger cluster sizes. We find that local are centred at the origin. In (d) the points are coloured 6

FIG. 3. Snapshot of magnetic moments of a system containing a spherical pure Cr precipitate with radius of 18A.˚ The simulation cell contains 31,250 atoms with 9 at.% Cr. The magnetic system is thermalized to 100K. (a) Magnetic moments of Fe is in orange, where Cr is in white. (b) Magnetic moments colourised according to double conical space. (c) Focused view near the interface between precipitate and bulk. (d) Magnetic moment space map of all moments. (e) Magnetic moment space P map colourised according to the sum of the occupational site variables of the nearest neighbours, e.g. σj = 8 means all neighbours are Fe, 0 means 50% is Fe and 50% is Cr, and -8 means all are Cr). according to the elemental species. In (e), the colour quite general. They relax for low and high index in- corresponds to the number of Fe atoms in the first coor- terfaces, and at elevated temperatures, as a method of dination shell of the atom producing the magnetic mo- minimising free energy subjected to magnetic frustration. ment. Red represents 8 Fe nearest neighbours whereas This highlights the need to consider non-collinear mag- blue represents 8 Cr. Again, Fe moments in the α are netism in the study of steels. observed to be ferromagnetic. Fe magnetic moments on Knowledge that we just obtained can be applied to the interface have less deviance from the magnetisation realistic materials such as EUROFER-97. This is the axis. The antiferromagnetism of the Cr atoms in the α0 is European reference steel for the first wall of a DEMO shown by the clustering points on the left and right of the fusion reactor which has a Cr concentration of 9 at.%. As moment space representation with a void in the centre. we can ascertain from Fig. 2, the change of magnetisation The arcing of these points show the moments which are is delimited by curves representing the disordered solid tilted in the vicinity of the interface. The bulge at the solution (0%) and complete phase separation ((1 − x)| < base reveals the magnitude of Cr magnetic moments are MFe > |). strengthened when surrounded by Fe atoms, i.e. those in The first wall operating temperature of a fusion reactor the α phase. is expected to be in the range of 350◦C to 550◦C63 (i.e. Non-collinear arrangements of Cr and Fe moments 623 to 823K). We calculated the magnetisation for dif- across the interface of the precipitate are similar to ob- ferent Cr concentrations at 600K and 800K. Details can servations of (011) Fe-Cr interfaces16,62. Our simulations be found in the Supplementary Material. The difference on clusters show that these magnetic configurations are between the magnetisation of the solid solution and com- 7

FIG. 4. The Curie temperature versus the nominal Cr concentration for different α0 precipitate sizes with radius of 0, 6, 10, 14, 18 and 22A.˚ The Curie temperature is maximised when the effective Cr concentration in the Fe matrix is about 5-6 at.%. plete phase separation at 9 at.%Cr is 10.4% at 600K and 10.3% at 800K. Our analysis shows that there are mea- surable changes to the macroscopic magnetic properties, which arise due to Cr clustering. In Fig. 4, we plot the TC as a function of nominal composition with different cluster sizes. TC is measured by the turning point in the magnetic susceptibility. Our calculations show that the TC of Fe1−xCrx can vary significantly for a given nominal concentration x. For x = 18at.% Cr, we find an 80K difference in the TC between Cr dispersed randomly compared to cases where most atoms are in a large cluster. FIG. 5. Magnetic properties of a spherical pure Cr precipi- As a metric to quantify the conditions which maximise tate with radius rc = 4 nm in pure Fe. (a) Radially averaged the Curie temperature, we consider the effective concen- magnetisation and (b) susceptibility evaluated as a function 0 tration x . This is defined as the concentration of the Cr of distance |r−r0| from the centre of the precipitate. (c) Time remaining the α phase (see supplementary material for averaged magnetic configurations for |r − r0| = (i) 10.1A,˚ (ii) further details). Indeed for each concentration, we iden- 35.1A,˚ (iii) 40.0A˚ and (iv) 69.97A˚ at T =100K. tify that TC is maximised when enough Cr conglomerates to produce an effective concentration of approximately 6 at.%Cr. The effective concentration x0 in the systems higher temperatures. with nominal concentrations x of 5, 12 and 18 at.%Cr Our simulations show the alloy acts as a composite are shown explicitly in Figure 4 where TC is maximised system. The α phase, with an effective concentration 0 to highlight this trend. of x , behaves as a disordered solid solution Fe1−x0 Crx0 . 0 0 The increase of TC when x ≈ 6 at.%Cr is not coinci- Cr within the α precipitate are bulk-like with an anti- dental. It can be attributed to the same mechanism that ferromagnetic ordering below the Néeltemperature. In drives TC in solid solutions to be maximised at x = 6.25 the succeeding sections we consider the dynamics at the at.% Cr17. When more Cr atoms are in the precipitate, cluster interface and how Fe substitutions within the α0 the relative concentration of Cr within the α phase de- change the magnetic properties. creases. This increases the average number of Fe nearest neighbours and decreases the number of nearest and next-nearest neighbouring Cr atoms. This reduces the C. Spatial and Time Resolved Magnetic Properties 0 geometric magnetic frustration of the Cr magnetic mo- of α Clusters ments in the α phase. They can then align antiferromag- netically with the adjacent Fe atoms which is energeti- For larger cluster sizes, we can study the magnetic cally favourable. This in turn strengthens the magnetic properties with spatial resolution. We demonstrate this moment of the Fe atoms which can remain correlated to by examining a spherical Cr cluster, of radius 4 nm with 8

FIG. 6. Radially resolved autocorrelation and first nearest neighbour correlation functions at T =100K, 500K, 800K and 1200K. They are taken as averaged values for magnetic moments located at 13.1, 35.1, 37.0, 40.1, 45.0 and 70.0A˚ away from the centre of a Cr precipitate in pure Fe with a radius of 40A.˚ no impurities, situated within pure Fe. We used a simu- haviours are due to the long range nature of the exchange lation cell containing 250,000 atoms, that is 50 × 50 × 50 coupling parameters. unit cells. This corresponds to a number density of We repeated the calculation three times with slightly 23 −3 3.2×10 m which can be observed experimentally, different initial magnitudes of the magnetic moments such as experiments on neutron irradiated Fe-12Cr at such that different phase space trajectories were explored. ◦ 61 320 C by Reese et al. . We did this in order to confirm the oscillations are phys- In Fig. 5 the radially averaged magnetisation and sus- ical, and not caused by noise or poor convergence. It ceptibility for temperatures from 100K to 1500K are plot- also improves our statistical sampling. The average is ted as a function of distance from the centre of the cluster plotted as the solid line for each temperature with the shaded region limiting the maximum and minimum val- |r − r0|. The magnetisation in the pure Fe surrounding a Cr cluster are found to have bulk values. When ues of the repeated runs. We confirm the oscillations ˚ ˚ approaching the Cr-Fe interface from the α phase, the at 32A < |r − r0| < 40A are physical by noting only magnetisation drops when the Fe coordination shells be- small variances between the peaks of all the calculations. come populated by Cr atoms. The magnetisation drops Furthermore, the peaks occur at fixed distances across sharply within the α0. Both the susceptibility and aver- different temperatures. aged magnetic moment are observed to oscillate within To visualise the local magnetic configurations, the time a penetration depth of 1 to 2nm in the α0. These be- averaged magnetic moments within a representative con- 9 ventional unit cell centred at each radius is shown in Fig. of 20 to 40ps is generally observed. 5(c). Four radii of interest were selected (i) 10.1A˚ (ii) At T = 100K the RRACF is positive for all radii. 35.1A˚ (iii) 40.0A˚ and (iv) 69.97A.˚ Cr magnetic moments In the α phase away from the cluster, the RRACF re- situated deep within the cluster are arranged antiferro- mains close to unity indicating the magnetic moments magnetically. They are perpendicular to the ferromag- are strongly correlated over time. The Fe moments on netic Fe moments in the α. The magnitudes of Cr mo- the interface (40.1A)˚ have a lower RRACF due to the ments on the interface are enhanced. At the interface magnetic frustration reducing the correlation. This coin- their orientations are locked antiparallel to the interfacial cides with the sharp increase in magnetic susceptibility Fe moments. Cr atoms further away from the interface along the interfacial Fe layers. Near the centre of the pre- rotated in plane. A region of frustration exists within a cipitate (13.1A),˚ the Cr magnetic moments have bulk-like penetration depth of 10A˚ from the interface were the Cr ordering and are strongly correlated. moments are less correlated. These properties agree with Near the interface the magnetisation within the clus- the smaller scale calculations in the previous section. ter fluctuates. As one can see in Fig. 5(a), the radially The susceptibility within a cluster is particularly high averaged magnetisation peaks at 37.0A˚ and troughs at when the temperature is below the Néeltemperature 35.1A.˚ The RRACF shows that at 100K, the moments TN ≈ 350K. This is due to the magnetic frustration of are more correlated at 37A˚ where a peak is observed in the Cr near the interface. A peak in susceptibility is ob- the magnetic moment (see Fig. 5a). This contrasts with served within the Fe matrix even at a temperatures much 35.1 A˚ which corresponds to a radius of reduced mag- lower than TC . From Eq. 6, it is apparent that the Fe netisation. Further for 35.1A,˚ the RRACF is weakest moments near the interface fluctuate more readily than indicating localised frustration. in the bulk Fe. We can deduce that the disorder of the The 1nRRCF provides a measure of correlation be- configuration at the Fe-Cr interface leads to a reduction tween a magnetic moment at |r − r0| with its nearest of TC locally. This is evidenced further in the analysis of neighbours. At 100K, it shows that the Fe moments in correlation functions below. the α phase remain extremely collinear with their neigh- We note that the magnetisation level in the centre of bours since c1 ≈ 1. Collinearity is decreased slightly for a pure Cr cluster is not unique. This indicates there the Fe magnetic moments near the Cr cluster (40.1A).˚ Cr is metastability enabling the magnetic moments to or- magnetic order is bulk-like near the centre of the precip- der according to initial conditions akin to a spin glass. itate at 13.1A˚ where c1 ≈ −1. Towards the Fe interface The system also meets the criteria of frustration due to the Cr moments become less correlated with their nearest the local configuration dependence of the exchange pa- neighbours over time as indicated by the increased value rameters. Since Fe-Cr is known to become a spin glass of the 1nRRCF. This indicates the Cr moments are devi- for 81-84%Cr15, it would be interesting to determine if ating from AFM ordering. Indeed this assessment agrees Cr rich α0 precipitates form spin glasses locally for lower with Fig. 5c(ii) where magnetic moments are seen to be nominal concentrations. highly non-collinear with their neighbours. An advantage of SD compared to MC is its capability At 500K, Cr moments in the centre of the precipitate of simulating the dynamic time resolved magnetic prop- (13.1A˚ and 35.1A)˚ become uncorrelated, which is consis- erties. In Fig. 6, we plot the radially resolved autocorre- tent with TN ≈ 350K. These moments dephase within lation function (RRACF, Eq. 7) and the first nearest 20-30ps. Moments at 37.0A˚ from the precipitate centre neighbour radially resolved correlation function (1nR- are still weakly correlated with their neighbouring atoms, RCF, Eq. 8) of magnetic moments at temperatures of and the 1nRRCF becomes positive. This implies the lo- 100K, 500K, 800K and 1200K. They are measured at cal TN is enhanced for Cr near the interface in contrast different distances from the centre of the cluster. We se- to the reduction in TC for interfacial Fe. lected six radii to compare. Three of them are within At 800K, Fe magnetic moments remain strongly corre- the precipitate, 13.1A,˚ 35.1A˚ and 37.0A.˚ The next is lated as shown in RRACF and 1nRRCF albeit collinear- at 40.1A,˚ which is at the first Fe layer surrounding the ity has decreased due to thermal fluctuations. Fe mo- Cr cluster. The remaining two are at 45.0A˚ and 70.0A,˚ ments near the α0 interface at 40.1A˚ remain less corre- which are in the α phase. lated than those in the Fe matrix and take a longer time The profile of the correlation functions c0 and c1 con- to dephase. These interfacial Fe moments have a higher sist of a decaying envelope function indicative of the char- oscillation frequency indicating that their effective fields F e−Cr acteristic relaxation or dephasing time, where its asymp- are stronger because of the strong J1nn interaction. totic value represents the strength of correlation. The The larger fluctuations give rise to the enhanced suscep- oscillations in the RRACF are due to the precession of tibility (Eq. 6). This provides further evidence that the 64 magnetic moments . Fe moments precess at a higher effective TC is reduced around the cluster. Consequently frequency than the Cr moments. If the rotational fre- the interface acts as a nucleation site for the ferromag- quency ωi = −|Hi|/¯h is considered, we immediately re- netic to paramagnetic phase transition. alise the effective field is greatest on the Fe moments. At 1200K, all magnetic moments become incoherent, Correspondingly the effective field is weak where there is as shown by both the RRACF and the 1nRRCF dropping reduced correlation such as at 35.1A.˚ A relaxation time to zero (Figure 6c). This is consistent with the loss of the 10 net magnetisation (Fig. 5a). The dephasing time of the To visualise the magnetic ordering in the system and bulk Fe moments with their neighbours is approximately the change with temperature, the magnetic moment ten times greater than Cr moments within the α0 cluster. space for each precipitate concentration is shown in Fig. 7(c). The moments are averaged over every timestep during the 5 ps SD simulations. For xc = 1 the Cr magnetic D. Impure Cr Clusters moments in the centre of the cluster arrange in antiferromagnetic order and perpendicular to the magnetisation In previous sections, we considered idealised spheri- axis of the Fe as was shown in Fig. 3. As the temper- cal α0 clusters. In reality, a Cr cluster is rich in Cr ature increases, the moments became more disordered. but impure61. Recent studies show Cr clusters formed Nonetheless, at 500K, which is above TN , Cr moments in neutron, ion and electron irradiated alloys, as well as retain some ordering and remain mostly parallel to the Fe in thermally aged materials, have sizes between 10 and magnetisation axis as indicated by the prolated distribu- 40A˚ and have local Cr concentrations between 50 and tion. At 1500K, when temperature is well above TN and 95%61,65–67. We investigate whether the same phenom- TC , the moments of Fe and Cr are spherically distributed ena exhibited by pure Cr clusters can also be observed indicating a paramagnetic state. with impure clusters as the α0 phase is diluted by Fe For xc = 0.8, 0.6 and 0.4, the Fe moments dispersed 0 substitutions. within the α precipitate experience significant frustra- 0 tion. This non-collinearity is shown in the magnetic mo- Spherical α clusters with Cr concentration xc and radius 40A˚ are created in a pure Fe matrix. Three ment space by the egg-shaped / ovoid distribution. The Fe moments maintain weak ordering with a net mag- configurations are generated for concentrations of xc = 0.8, 0.6, 0.4, 0.2 in a supercell with 250,000 atoms. Since netisation and with maxima along the major axis of the irradiated steels have been experimentally observed to ovoid. The Cr moments are similarly distributed albeit have α0 clusters with concentrations between 50 and 95% with a weak total magnetic moment orientated oppo- we report the x = 0.8 and 0.6 cases in Fig. 7. Figures site to that of Fe moments. Surprisingly at 500K, the c collinearity increases. The thermal agitation enables Fe for xc =0.4 and 0.2 are in Supplementary Material. In Fig. 7, column (i) shows the radially averaged magneti- moments perpendicular to the magnetisation axis to re- orient. We note that experiments have confirmed spin sation and column (ii) the magnetic susceptibility as a 15 function of distance from the centre of the cluster. The glass state 81-84 at.% Cr at low temperatures . This region between the smallest and largest values measured spin glass state transits into ferromagnetic state as tem- for each simulation is shaded and the average of the three perature increases. Our results suggest Cr clusters share runs is plotted as a solid line. glassy characteristics. In summary, the idealised pure Cr precipitates have At low temperatures (T < TC ), the magnetisation in the Fe matrix recovers the bulk values typically within an ordered non-collinear AFM structure which becomes 10A˚ from the α0 interface. When the temperature in- more disordered with Fe substitutions. Despite the reduction of magnetic ordering, the overall properties re- creases, the recovery length increases for all xc. At 0 main similar. The effective TC of Fe surrounding the α 1100K, the recovery length is 20-30A˚ for all xc, even where there is no Cr in the Fe matrix. This indicates the cluster is reduced whereas TN of interfacial Cr atoms is interactions and associated reduced ferromagnetic order- increased. This results in a spatial dependence of time- ing proximate to the interface have long range effects, resolved properties such as the susceptibility which was which reduces their correlated reinforcement. used to identify interfacial Fe to be more susceptible than For higher concentrations of Fe within the α0, the to- Fe in the α phase. tal magnetisation within the cluster steadily increases. A discontinuous change in the radially averaged magnetisation is observed across the interface. The total magneti- E. Ideal (001) Interface of Fe-Cr Superlattices sation is enhanced within 10A˚ of the precipitate interface. The enhanced magnetic susceptibility of the Fe atoms For large clusters, we may approximate the interface by near the interface discussed in the previous section is ideal flat interfaces with different orientations. Indeed, it also observed for the non-ideal clusters (see arrow in has been shown that when non-collinear magnetic config- Fig.7(b)(ii)). When the local Cr concentration decreases, urations are permitted, the formation energies of (001), the enhancement decreases. Nonetheless, even for the (011) and (111) Fe/Cr interfaces predicted using density 54 xc = 0.2 case the susceptibility of Fe on the interface is functional theory are nearly degenerate . As such, in still increased relative to the rest of the α phase. This this section, we limit our investigation to the study of indicates that α0 acts as a nucleation site for the phase the magnetic structure of a (001) Fe/Cr interface. transitions even for clusters with dilute Cr concentration. Spin dynamics calculations are performed on a cell con- Within the cluster, the susceptibility is reduced near taining 32000 atoms, where 20 × 20 × 20 Cr unit cells the interface. This corresponds to an increased effective stacked upon 20×20×20 Fe unit cells. Periodic boundary TN originating from the interfacial Cr moments coupling conditions are enforced. Essentially, we create an Fe/Cr to Fe moments. superlattice with two interfaces in the principle cell. The 11

FIG. 7. (i) Radially averaged magnetisation and (ii) magnetic susceptibility as a function of radius from the centre of a spherical α0 cluster with radius of 40A.˚ The Cr precipitate has a concentration of (a) 80 at.%Cr, (b) 60 at.%Cr. Concentrations of 40 at.%Cr and 20 at.%Cr are provided in Supplementary Material. Three samples are simulated for each case. The averaged value is shown as the solid line. The region between the smallest and largest values measured for each simulation is shaded. (c) Time averaged magnetic moment space diagrams for the Cr clusters of different local concentrations at 100K, 500K, 900K and 1500K. Each dot represent a magnetic moment centred at the origin. Fe moments are in red, where Cr moments are in blue. Each map is shown with respect to a scale of 1µB . Magnetic moments were averaged over every time step during the 5 ps production run. 12 time averaged magnetic properties of such a confirgura- tion has been studied using MC MCE simulations in Ref. 54. In the Supplementary Material we provide the averaged magnetisation of the superlattice system and the layer resolved magnetic moments which we compare to these earlier calculations by Lavrentiev et al.54. The magnetisation almost equals to the scaled value for bulk Fe implying the Fe magnetic moments in the Fe layers remain bulk like at low and intermediate temperatures. It is similar to the case of large clusters with high Cr concentration. At higher temperatures the magnetisation of the interfacial system reduces faster than the bulk Fe (Supplementary Fig. 3). We can infer the TC is reduced by approximately 100K. From the layer averaged calculations, the magnetisation of the middle Cr layers, the averaged magnetisation vanishes between 300 and 400K corresponding to the Néeltemperature (TN ≈ 350K). A finite averaged moment in Cr is observed within a penetration depth of about 20A˚ for all temperatures < TC . The middle layers of the Fe recover bulk magnetisation values at low temperatures. For Fe in layers close to the interface, the magnetisation is sup- pressed. The bulk magnetisation recovery length from the interface increases with temperature. There is excellent agreement between our SD calculations with the MC simulations for these time averaged properties. Using SD, FIG. 8. (a) Schematic view of an Fe/Cr superlattice with we can also study the continuous evolution of magnetic (001) interfaces, and a snap-shot of the magnetic moments moments as a function of time. We analyse the time cor- at the interface at 100K. Fe moments are in orange. Cr mo- relation between spins by considering correlation func- ments are in grey. (b) Magnetic susceptibility as a function of temperature for each (001) layer in a Fe/Cr superlattice. tions across a simple (001) interface. These results are compared to the dynamics observed for more complex interfaces of clusters. This allows us to conclude whether the behaviour identified in Section III C is general or is a consequence of the topology of a cluster. to be resolved per (001) layer in Figure 9. Each layer is In Fig. 8(a), we show a snapshot of the magnetic mo- labelled by the elemental type occupying the layer (either ments near the interface at 100K. The Fe moments re- Fe or Cr) and a number indexing the layers proximity to main collinear whereas the orientation of the Cr moments the interface, such as Fe(1) or Cr(4). Index 1 represents are strongly dependent upon their proximity to the in- interfacial layers (sample indices are shown in Figure 8a). terface where they have a larger proportion of Fe in their first five coordinate shells. The first interfacial Cr layer In Fig. 9, the autocorrelation function c0(z, τ) shows is orientated antiparallel to the Fe moment alignment. the 1st and 3rd Fe layers (Fe(1,3)) at 100K are less cor- The deeper into the Cr slab, the magnetic moments re- related than the 2nd layer. This corresponds to the dip orientate gradually until perpendicular to the Fe slab. recorded in the susceptibility below 800K (see Fig. 8b). Since each Cr sublayer is attempting to minimise energy Above 800K, the 2nd Fe layer becomes increasingly un- by aligning antiparallel to one another, there is substan- correlated with respect to the interface layer Fe(1). tial frustration near the interface. This phenomenon was likewise observed near the Cr cluster interface. The interfacial Cr layers Cr(1), remain strongly corre- In Fig. 8(b), the susceptibility resolved for each (001) lated in comparison to the other Cr layers. A minimum is layer in the simulation cell is shown between 100 K and observed in the susceptibility at the interface for all tem- 1500 K with 100 K increments. Remarkable similarity is peratures correspondingly. At 800K, almost twice TN , 0 observed with respect to our results for an ideal α cluster c1(z, τ) identifies the moments in Cr(1) are still corre- (Fig. 5(b)). Again, enhanced susceptibility is identified lated to the moments of the nearest neighbouring atoms. for the Fe moments near the Fe/Cr interface. It becomes This behaviour of the susceptibility and correlation is clear that presence of Cr destabilises the Fe ordering near observed for both the (001) interface and the high in- an interface irrespective to the index of the plane. dex interfaces of spherical α0 clusters. This shows simple Additional evidence can be observed from the time cor- (001) interface are good approximations to more complex relation functions (Eq. 7 and 8) which are now defined interfaces. 13

IV. CONCLUSION

Spin dynamics simulations were performed to investigate the magnetic properties in Fe1−xCrx alloys. Differ- ent Cr concentrations and degrees of Cr clustering were considered. We adopted a magnetic cluster expansion Hamiltonian that considers both longitudinal and trans- verse fluctuations of magnetic moments within the framework of Langevin spin dynamics. We observed a strong dependence between the magnetic properties of the FeCr alloys with the local microstructure. Using the SPILADY code48 we compared the change of Curie temperature between disordered solutions and systems with precipitation. In both cases the Curie temperature was determined from the magnetic susceptibility. For disordered solid solutions, the Curie temperature was maximised for a Cr concentration of approximately 5 at.%. When Cr clustering was considered, the total magnetisation for any nominal Cr concentration increased with the size of the Cr precipitate at low temperatures and was minimum for solid solutions. Configurational disorder in the solid solution created geometric magnetic frustration that reduced the magnetisation. This means that the higher is the level of Cr precipitation, the higher is the total magnetisation. When we observed the change in the Curie temperature for nominal concentrations of Cr with precipitates of varying size, we found a dependence upon the effective concentration. The effective concentration is the Cr concentration outside the α0 cluster. The transition temperature was found to be highest for configurations with an effective α concentration of 5-6 at.%. For a particular case of Fe-9Cr alloys, which has the same Cr concentration as EUROFER-97, the magnetisation can vary by 10% over the operational temperature ranges of the first wall of fusion power plant, for example from 600K to 800K. The change of magnetisation is due to the varying level of Cr precipitation. This potentially has significant impact on the integrity of the mechanical properties of steels, which is coupled to its magnetic properties. Cr located at the α-α0 interface were observed to have moments that remained strongly correlated over time and with respect to the neighbouring moments. The Cr moments were aligned anti-parallel to the moments of adjacent Fe atoms. This reduced their susceptibility relative to Cr in the core of the cluster. For pure Cr precipitates, the moments were identified to rotate from being anti-parallel with Fe at the interface, FIG. 9. Time-displaced spin-spin autocorrelation and near- to being perpendicular in the centre. The Cr moments est neighbour correlation functions of each (001) layer in a were in an antiferromagnetic configuration with respect Fe/Cr superlattice at 100K, 800K and 1200K. Fe/Cr inter- to their nearest neighbours. Cr moments near the interfaces are located between layers 80 & 1, and 40 & 41. Curves face, where the non-collinear configuration was observed, of interest are labelled according to the atoms in the layer and their proximity to the interface. The Fe interfacial lay- were determined to be the least correlated and therefore ers (1,40) are labelled as Fe(1). The next layers (2,39) are frustrated. For impure clusters of experimentally consis- labelled Fe(2), Cr interface layers (41,80) are labelled Cr(1), tent Fe concentrations (50-95%Cr), the ordered antifer- and so on. romagnetism in the centre of the cluster was destroyed. Instead, disordered non-collinear configurations were ob- 14 served at all temperatures. A net magnetisation was still analysis of dynamic magnetism as a function of both mi- supported. crostructure and an external magnetic field would be in- When temperature was increased, the Fe moments sightful future work for reactor material design. close to the interfaces became more disordered than in In summary, using spin dynamics we have investigated the bulk. The Cr cluster interface acts as a nucleation the magnetisation, the susceptibility and the Curie site with a reduced effective Curie temperature. Con- temperature of Fe1−xCrx alloys with different Cr con- versely, the Cr moments on the interface remain statis- centration and different degrees of Cr precipitation. This tically correlated at temperatures far exceeding the Néel enhances our knowledge of how magnetic properties of point. steels change upon ageing and under irradiation damage. We may understand the underlying physics through an This may help to find a non-invasive way of determining investigation into dynamic properties which can be calcu- the level of Cr precipitation through the observation of lated from spin dynamics calculations in contrast to the magnetic properties of alloys. Monte Carlo method. We studied the spin-spin autocorrelation and 1st nearest neighbour correlation functions in Cr clusters and (001) Fe-Cr superlattices. Contrary to expectation, the high index of the cluster interface did not significantly change the observed correlation of ACKNOWLEDGMENTS Fe and Cr across the interface with respect to the simple (001) case. Fe magnetic moments near the interface We would like to express our gratitude to Mikhail Yu. fluctuate more readily than in the bulk of the α. This Lavrentiev for valuable discussions and providing data resulted in enhanced susceptibility of the Fe in the inter- from his work. This work has been carried out within facial layers. The dynamic properties observed for the the framework of the EUROfusion Consortium and has impure clusters remained consistent with the pure clus- received funding from the Euratom research and training ters and (001) interface. programme 2014-2018 and 2019-2020 under grant agree- The effect of external magnetic field, which was not ment No. 633053 and No. 755039 and from the RCUK considered here, may indeed worth attention from the Energy Programme [grant number EP/P012450/1]. SLD technological perspective. In the design of ITER68, acknowledges support from the Centre for Non-Linear the maximum toroidal magnetic field will be 12T and Studies at LANL. To obtain further information on the poloidal field will be 6T. According to a spin-lattice dy- data and models underlying this paper please contact namics study on pure iron69, they may increase the tran- [email protected]. The views and opinions sition temperature by approximately 1K T−1. Further expressed herein do not necessarily reflect those of the verification is required for Fe-Cr alloys. A quantitative European Commission.

∗ Corresponding author: [email protected] structure in body-centred cubic metals,” Phys. Rev. Mater. 1 D. Stork, P. Agostini, J.-L. Boutard, D. Buckthorpe, (2019). E. Diegele, S. L. Dudarev, C. English, G. Federici, M. R. 7 C. C. Fu, F. Willaime, and P. Ordejón,“Stability and Gilbert, S. Gonzalez, A. Ibarra, C. Linsmeier, A. Li Puma, mobility of mono- and di-intersistials in α-fe,” Physical G. Marbach, L. W. Packer, B. Raj, M. Rieth, M. Q. Review Letters 92, 175503 (2004). Tran, D. J. Ward, and S. J. Zinkle, “Materials r&d for a 8 P. Olsson, C. Domain, and J. Wallenius, “Ab initio study timely demo: Key findings and recommendations of the eu of cr interactions with point defects in bcc fe,” Physical roadmap materials assessment group,” Fusion Engineering Review B 75, 014110 (2007). and Design 89, 1586 (2014). 9 D. J. Dever, “Temperature dependence of the elastic con- 2 H. Hasegawa and D. G. Pettifor, “Microscopic theory of stants in α-iron single crystals: relationship to spin order the temperature-pressure phase diagram of iron,” Physical and diffusion anomalies*,” Journal of Applied Physics 43, Review Letters 50, 130 (1983). 3293 (1972). 3 P.-W. Ma, S. L. Dudarev, and J. S. Wróbel, “Dynamic 10 H. Hasegawa, M. W. Finnis, and D. G. Pettifor, “A cal- simulation of structural phase transitions in magnetic culation of elastic constants of ferromagnetic iron at finite iron,” Physical Review B 96, 094418 (2017). temperatures,” Journal of Physics F: Metal Physics 15, 19 4 D. Nguyen-Manh, A. P. Horsfield, and S. L. Dudarev, (1985). “Self-interstital atom defects in bcc transition metals: 11 E. Fawcett, “Spin-density-wave antiferromagnetism in Group-specific trends,” Physical Review B 73, 020101 chromium,” Reviews of Modern Physics 60, 209 (1988). (2006). 12 A. W. Overhauser, “Spin density waves in an electron gas,” 5 P. M. Derlet, D. Nguyen-Manh, and S. L. Dudarev, “Mul- Physical Review 128, 1437 (1962). tiscale modelling of crowdion and vacancy defects in body- 13 R. Hafner, D. Spiˇsák,R. Lorenz, and J. Hafner, “Magnetic centred-cubic transition metals,” Physical Review B 76, ground state of cr in density-functional theory,” Phys. Rev. 054107 (2007). B 65, 184432 (2002). 6 Pui-Wai Ma and S. L. Dudarev, “University of point defect 14 S. K. Burke and B. D. Rainford, “The evolution of mag- 15

netic order in crfe alloys: I. antiferromagnetic alloys close (1992). to the critical concentration,” Journal of Physics F: Metal 32 S. Dubiel and J. Zukrowski, “Change of cr atoms distribu- Physics 13, 441 (1983). tion in fe85-cr15 alloy caused by 250 kev he+ ion irradia- 15 S. K. Burke, R. Cywinski, J. R. Davis, and B. D. Rainford, tion to different doses,” Journal of Alloys and Compounds “The evolution of magnetic order in crfe alloys .2. onset of 624, 165 (2015). ferromagnetism,” Journal of Physics F: Metal Physics 13, 33 M. Miller, J. Hyde, M. Hetherington, A. Cerezo, G. Smith, 451 (1983). and C. Elliott, “Spinodal decomposition in fe-cr alloys: ex- 16 M. Y. Lavrentiev, K. Mergia, M. Gjoka, D. Nguyen-Manh, perimental study at the atomic level and comparison with G. Apostolopoulos, and S. L. Dudarev, “Magnetic clus- computer models i. introduction and methodology,” Acta ter expansion simulation and experimental study of high Metallurgica et Materialia 43, 3385 (1995). temperature magnetic properties of fe-cr alloys,” Journal 34 A. Kuronen, S. Granroth, M. H. Heinonen, R. E. Perälä, of Physics: Condensed Matter 24, 326001 (2012). T. Kilpi, P. Laukkanen, J. L˚ang,J. Dahl, M. P. J. Punkki- 17 D. Nguyen-Manh, M. Y. Lavrentiev, and S. L. Dudarev, nen, K. Kokko, M. Ropo, B. Johansson, and L. Vitos, “Ab initio and monte carlo modeling in fecr system: Mag- “Segregation, precipitation, and α − α0 phase separation netic origin of anomalous thermodynamic and kinetic prop- in fe-cr alloys,” Physical Review B 92, 214113 (2015). erties,” Computational Materials Science 44, 1 (2008). 35 I. Dopico, P. Castrillo, and I. Martin-Bragado, “Quasi- 18 J. O. Andersson and B. Sundman, “Thermodynamic prop- atomistic modeling of the microstructure evolution in bi- ertis of the cr-fe system,” CALPHAD 11, 83 (1987). nary alloys and its application to the fecr case,” Acta Ma- 19 W. Xiong, M. Selleby, Q. Chen, J. Odqvist, and Y. Du, terialia 95, 324 (2015). “Phase equilibria and thermodynamic properties in the fe- 36 L. D. Landau and L. M. Lifshitz, “On the theory of the dis- cr system,” Crtical Reviews in Solid State and Materials persion of magnetic permeability in ferromagnetic bodies,” Science 35, 125 (2010). Physikalische Zeitschrift der Sowjetunion 8, 153 (1935). 20 G. Bonny, D. Terentyev, and L. Malerba, “On the α − α0 37 T. L. Gilbert, “A phenomenological theory of damping in miscibility gap of fe-cr alloys,” Scripta Materialia 59, 1193 ferromagnetic materials,” IEEE Transactions on Magnetics (2008). 40, 3443 (2004). 21 T. P. C. Klaver, R. Drautz, and M. W. Finnis, “Magnetism 38 X. Tao, D. P. Landau, T. C. Schulthess, and G. M. Stocks, and thermodynamics of defect-free fe-cr alloys,” Physical “Spin waves in paramagnetic bcc iron: spin dynamics sim- Review B 74, 094435 (2006). ulations,” Physical Review Letters 95, 087207 (2005). 22 P. Olsson, I. A. Abrikosov, L. Vitos, and J. Wallenius, 39 P. W. Ma, C. H. Woo, and S. L. Dudarev, “High- “Ab initio formation energies of fe-cr alloys,” Journal of temperature dynamics of surface magnetism in iron thin Nuclear Materials 321, 84 (2003). films,” Philosophical Magazine 89, 2921 (2009). 23 P. J. Grobner, “The 885◦ f (475◦ c) embrittlement of fer- 40 P. W. Ma and S. L. Dudarev, “Longitudinal magnetic fluc- ritic stainless steels,” Metallurgical Transactions 4, 251– tuations in langevin spin dynamics,” Physical Review B 260 (1973). 86, 054416 (2012). 24 C. D. Hardie, C. A. Williams, S. Xu, and S. G. Roberts, 41 J. Hubbard, “Magnetism of iron,” Physical Review B 19, “Effects of irradiation temperature and dose rate on the 2626 (1979). mechanical properties of self-ion implanted fe and fe-cr al- 42 W. F. Brown Jr., “Thermal fluctuations of a single-domain loys,” Journal of Nuclear Materials 439, 33 (2013). particle,” Physical Review 130, 1677 (1963). 25 D. Stork, P. Agostini, J. L. Boutard, D. Buckthorpe, 43 J. L. Garcia-Palacios and F. J. Lázaro, “Langevin- E. Diegele, S. L. Dudarev, C. English, G. Federici, M. R. dynamics study of the dynamical properties of small mag- Gilbert, S. Gonzalez, A. Ibarra, Ch. Linsmeier, A. Li. netic particles,” Physical Review B 58, 14937 (1998). Puma, G. Marbach, P.F. Morris, L. W. Packer, B. Raj, 44 P. W. Ma, S. L. Dudarev, A. A. Semenov, and C. H. Woo, M. Rieth, Tran M. Q., D. J. Ward, and S. J. Zinkle, “Temperature for a dynamic spin ensemble,” Physical Re- “Developing structural, high-heat flux and plasma facing view E 82, 031111 (2010). materials for a near-term demo fusion power plant: The eu 45 S. Chandrasekhar, “Stochastic problems in physics and as- assessment,” Journal of Nuclear Materials 455, 277 (2014). tronomy,” Reviews of Modern Physics 15, 1 (1943). 26 S. Novy, P. Pareige, and C. Pareige, “Atomic scale analysis 46 R. Kubo, “The fluctuation-dissipation theorem,” Reports and phase separation understanding in a thermally aged on Progress in Physics 29, 255 (1966). fe-20 at.% cr alloy,” Journal of Nuclear Materials 384, 96 47 P. W. Ma and S. L. Dudarev, “Langevin spin dynamics,” (2009). Physical Review B 83, 134418 (2011). 27 F. Soisson and T. Jourdan, “Radiation-accelerated precip- 48 P.-W. Ma, S. L. Dudarev, and C. H. Woo, “Spilady: A itation in fe-cr alloys,” Acta Materialia 103, 870 (2016). parallel cpu and gpu code for spinlattice magnetic molec- 28 E. Little and D. Stow, “Void-swelling in irons and ferritic ular dynamics simulations,” Computer Physics Communi- steels: Ii. an experimental survey of materials irradiated in cations 207, 350 (2016). a fast reactor,” Journal of Nuclear Materials 87, 25 (1979). 49 P. W. Ma and C. H. Woo, “Parallel algorithm for spin 29 M. Bachhav, G. R. Odette, and E. A. Marquis, “Mi- and spin-lattice dynamics,” Physical Review E 79, 046703 crostructural changes in a neutron irradiated fe-15at.%cr (2009). alloy,” Journal of Nuclear Materials 454, 381 (2014). 50 M. Y. Lavrentiev, S. L. Dudarev, and D. Nguyen-Manh, 30 M. Furusaka, Y. Ishikawa, S. Yamaguchi, and Y. Fujino, “Magnetic cluster expansion simulations of fecr alloys,” “Phase separation process in fe-cr alloys studied by neu- Journal of Nuclear Materials 386, 22 (2009). tron small angle scattering,” Journal of the Physical Soci- 51 M. Y. Lavrentiev, D. Nguyen-Manh, and S. L. Dudarev, ety of Japan 55, 2253 (1986). “Cluster expansion models for fe-cr alloys, the prototype 31 F. Bley, “Neutron small-angle scattering study of unmixing materials for a fusion power plant,” Computational Mate- in fecr alloys,” Acta Metallurgica et Materialia 40, 1505 rials Science 49, S199 (2010). 16

52 M. Yu. Lavrentiev, D. Nguyen-Manh, and S. L. Dudarev, Odette, and E. A. Marquis, “On α0 precipitate compo- “Magnetic cluster expansion model for bcc-fcc transitions sition in thermally annealed and neutron-irradiated fe- 9- in fe and fe-cr alloys,” Physical Review B 81, 184202 18cr alloys,” Journal of Nuclear Materials 500, 192 (2018). (2010). 62 H. Fritzsche, S. Bonn, J. Hauschild, J. Klenke, K. Prokes, 53 D. Nguyen-Manh, M. Y. Lavrentiev, M. Muzyk, and S. L. and G. J. McIntyre, “Antiferromagnetic order of thin epi- Dudarev, “First-principles models for phase stability and taxial cr layers in an fe/cr(110) multilayer,” Physical Re- radiation defects in structural materials for future fusion view B 65, 144408 (2002). power-plant applications,” Journal of Materials Science 47, 63 M. Gasparotto, R. Andreani, L. V. Boccaccini, 7685 (2012). A. Cardella, G. Federici, L. Giancarli, G. Le Marois, 54 M. Yu. Lavrentiev, R. Soulairol, C. C. Fu, D. Nguyen- D. Maisonnier, S. Malang, A. Moeslang, Y. Poitevin, Manh, and S. L. Dudarev, “Noncollinear magnetism at B. van der Schaaf, and M. Victoria, “Study of in-vessel interfaces in iron-chromium alloys: The ground states and candidate materials for fusion power plants- the european finite-temperature configurations,” Physical Review B 84, materials r&d programme,” Fusion Engineering and 144203 (2011). Design 66, 129 (2003). 55 D. Perera, M. Eisenbach, D. M. Nicholson, G. M. Stocks, 64 P.-W. Ma, C. H. Woo, and S. L. Dudarev, “Large-scale and D. P. Landau, “Reinventing atomistic magnetic sim- simulation of the spon-lattice dynamics in ferromagnetic ulations with spin-orbit coupling,” Physical Review B 93, iron,” Physical Review B 78, 024434 (2008). 060402 (2016). 65 O. Tissot, C. Pareige, E. Meslin, B. Decamps, and 56 S. V. Halilov, A. Ya. Perlov, P. M. Oppeneer, A. N. J. Henry, “Kinetics of alpha’ precipitation in an electron Yaresko, and V. N. Antonov, “Magnetiocrystalline irradiated fe15cr alloy,” Scripta Materialia 122, 31 (2016). anisotropy energy in cubic fe, co, and ni: Applicability 66 S. V Rogozhkin, O. A. Korchuganova, and A. A. Aleev, of local-spin-density theory reexamined,” Physical Review “Microstructural evolution of fe22%cr model alloy under B 57, 9557 (1998). thermal ageing and ion irradiation conditions studied by 57 Y. Xie and J. A. Blackman, “Magnetocrystalline atom probe tomography,” Inorganic Materials: Applied anisotropy and orbital polarization in ferromagnetic tran- Research 7, 210 (2016). sition metals,” Physical Review B 69, 172407 (2004). 67 Wei-Ying Chen, Yinbin Miao, Yaqiao Wu, Kun Tomchik, 58 D Sherringdon and B. W. Southern, “Spin glass verses fer- Carolyn A.and Mo, Jian Gan, Maria A. Okuniewski, Stu- romagnet,” Journal of Physics F: Metal Physics 5, L49 art A. Maloy, and James F. Stubbins, “Atom probe study (1975). of irradiation-enhanced alpha’ precipitation in neutron- 59 K. Chen and D. P. Landau, “Spin-dynamics study of the irradiated fe-cr model alloys,” Journal of Nuclear Materials critical behavior of the three-dimensional classical heisen- 462, 242–249 (2015). berg ferromagnet,” Physical Review B 49, 3266 (1994). 68 INTERNATIONAL ATOMIC ENERGY AGENCY, Sum- 60 W. Xiong, P. Hedström,M. Selleby, J. Odqvist, M. Thu- mary of the ITER Final Design Report, ITER EDA Doc- vander, and Q. Chen, “An improved thermodynamic mod- umentation Series No. 22 (INTERNATIONAL ATOMIC eling of the fe-cr system down to zero kelvin coupled with ENERGY AGENCY, Vienna, 2001). experiments,” CALPHAD 35, 355 (2011). 69 C. P. Chui and Y. Zhou, “Spin-lattice dynamics simulation 61 E. R. Reese, M. Bachhav, P. Wells, T. Yamamoto, G. R. of external field effect on magnetic order of ferromagnetic iron,” AIP Advances 4, 037110 (2014).