Low-Temperature Crystal Structure and Mean-Field Modeling of Er x Dy 1 − x Al 2 Intermetallics

: Low-temperature crystal structure of the Er x Dy 1 − x Al 2 alloys with x = 0.45, 0.67, 0.90 was examined using temperature-dependent powder X-ray di ﬀ raction. The Er-rich sample, Er 0.9 Dy 0.1 Al 2 , exhibits a rhombohedral distortion associated with the magnetic ordering that occurs around 20 K. The rhombohedral distortion is suppressed in Er 0.67 Dy 0.33 Al 2 , while a weak low-temperature tetragonal distortion is observed in Er 0.45 Dy 0.55 Al 2 . The mean-ﬁeld theory supports the correlation between the type of structural distortion and the variable easy magnetization axis in Er x Dy 1 − x Al 2 intermetallics.


Introduction
The interactions between two or more elements containing 4f electrons is a fascinating topic [1][2][3] with potential practical importance for novel technologies such as magnetocaloric cooling [4][5][6]. The sheer number of known and yet to be discovered intermetallic compounds containing 4f elements is enormous, and this broad family becomes nearly infinite when considering the ability to partially substitute one lanthanide element with another. Further, intra rare earth substitutions involving twelve magnetic lanthanide ions (excluding non-magnetic Lu and La, as well as unstable Pm) lead to captivating basic science rooted in the steeply increased complexity of magnetic interactions. It is known that magnetic exchange between the localized, spatially separated 4f orbitals is mediated by the conduction electrons (Ruderman-Kittel-Kasuya-Yosida, RKKY-type interactions) [7][8][9][10]. In magnetic lanthanides with a non-zero orbital quantum number, L, the presence of spin-orbit coupling and site-dependent crystalline electric fields lead to the splitting of the 4f energy levels and their population by electrons creating unusually complex magnetic structures [10][11][12][13][14]. Further, the intimate coupling between magnetic and crystallographic sublattices often translates the magnetic complexity into crystallographic one (and vice versa) leading to a variety of magnetostructural phenomena [15][16][17]. Giant magnetostriction in (Tb 1−x Dy x Fe 2 ) [18,19] and giant room temperature magnetocaloric effect in (Gd 5 Si 2 Ge 2 ) [20,21] among others, take origin in the magnetoelastic coupling, signifying its practical and fundamental significance. Among broadly known model systems routinely tapped to study magnetostructural phenomena is the family of intermetallic Laves phases that adopt MgCu 2 -type crystal structure crystallizing in cubic space group Fd3m; in particular compounds exhibiting strong magnetocrystalline anisotropy in their magnetically ordered states [3,10,19,[22][23][24][25][26][27]. The RM 2 Laves phase compounds, where R is a rare earth element and M is a dor p-element, for example, Mn, Fe, Co, Ni or Al, show clear correlations between magnetic anisotropy and the type of crystallographic distortions that occur in parallel with magnetic ordering [10,28]. In such compounds, the easy magnetization axis (EMA) determines whether the material undergoes cubic to tetragonal (EMA <100>), rhombohedral (EMA <111>), or orthorhombic (EMA <110>) distortions. Consequently, many magnetic binary RM 2 compounds adopt low-temperature crystal structures related to their EMAs. In rare cases, e.g., in HoCo 2 , a compound may exhibit more than one crystallographic transformation reflecting spin-reorientation transitions that change EMA [29].
The structural transitions become harder to predict and model in pseudobinary systems containing two rare earth elements whose corresponding binary parents adopt different EMAs in the magnetically ordered states, and, consequently, different low-temperature crystal structures. Here, additional complexities associated with intra-lanthanide substitutions, such as lattice disorder (e.g., when different rare-earth atoms randomly occupy the same atomic site), and exchange interactions between 4f orbitals of different elements modified by crystalline electric field splitting come into play, producing unexpected results. For example, a recent study of Er x Dy 1−x Co 2 compounds shows that at x = 0.75 the compound is mimicking the behavior of HoCo 2 , including the presence of a second crystallographic transition [3]. However, earlier temperature-dependent X-ray powder diffraction study of the Er 0.75 Dy 0.25 Al 2 alloy [30], where behaviors similar to Er 0.75 Dy 0.25 Co 2 may be expected due to competition between Er and Dy ions showed no signs of any structural transformation below its Curie temperature, T C , despite the first-order nature of spin reorientation transition observed below T C clearly evidenced by heat capacity data.
Both ErCo 2 and ErAl 2 Laves phase compounds adopt <111> as EMA, and both undergo a rhombohedral distortion [24,26], while DyCo 2 and DyAl 2 are both reported to be tetragonal at low temperature [23,26]. In Er x Dy 1−x Co 2 the M element, Co, carries magnetic moment, and the onsets of magnetic ordering in the Co sublattice are responsible for the discontinuous volume changes observed during the first-order transitions in DyCo 2 and ErCo 2 . The non-magnetic aluminum, on the other hand, does not bring itinerant magnetism to bear, and Er x Dy 1−x Al 2 compounds should, in principle, exhibit magnetostructural behavior that reflects fundamental interactions between the Er and Dy 4f orbitals. With this in mind, we performed a temperature-dependent crystallographic study of the Er x Dy 1−x Al 2 compounds with x = 0.45, 0.67, and 0.90 in order to understand how the interactions between two magnetic rare-earth sublattices influence the low-temperature crystallography. The physical behaviors (heat capacity and magnetization) of these materials are known and reported [31,32], yet data about their low-temperature crystal structures are lacking. In this work, we also use previously published X-ray powder diffraction data of ErAl 2 , Er 0.75 Dy 0.25 Al 2 , and DyAl 2 compounds [23,24,30] for comparison. Further, we use the mean-field theory tested earlier on similar R'R"Al 2 pseudobinary systems [33,34] to explain how the low-temperature crystallographic behavior evolves with a composition by modeling how the EMAs of Er and Dy sublattices change as functions of temperature (T) and of Er concentration (x).

Experimental
The samples used in this study were the same specimens that were prepared and thoroughly investigated in [32]. The samples are stable at ambient conditions and the sample surface retains metallic luster showing no visible traces of oxidation or corrosion after storage in a standard laboratory climate-controlled atmosphere for~10 years. The samples were ground into fine powders, screened to eliminate particles greater than 25 µm, and mixed with GE varnish; the obtained paste was placed and solidified in a custom-made copper sample holder and the sample surface was polished flat. The temperature-dependent X-ray powder diffraction data were obtained on a Rigaku TTRAX rotating anode diffractometer (Rigaku, Japan) equipped with a continuous helium-flow cryostat and a superconducting magnet; the sample space is pumped down to 10 −6 Torr to ensure temperature stability The first terms in (1) and (2) represent exchange interactions for a given concentration of x and 1 − x of Er and Dy, respectively. The g and J are the Landé factors (g Er = 6 5 and g Dy = 4 3 ) and total angular momentum quantum numbers (J Er = J Dy = 15/2), respectively, µ B is the Bohr magneton, and λ Er = 0.07 meV, λ Dy = 0.2621 meV and λ Er,Dy = 0.21 meV are the exchange parameters for the Er-Er, Dy-Dy and Er-Dy interactions, respectively [33]. The H Er,Dy CEF terms represent the crystalline electrical field (CEF) Hamiltonian for a cubic symmetry in Lea, Leask, and Wolf (LLW) notation [41,42]. Using cubic symmetry for the low-temperature structures is a valid approximation considering that the studied samples, as shown experimentally below, either develop only minor lattice distortions or do not exhibit measurable distortions at all. The CEF parameters used for each rare-earth sublattice are: (i) Er sublattice: F 4 = 60, F 6 = 13, 860, X = −0.2620 and W = − 0.0252 meV; (ii) Dy sublattice: F 4 = 60, F 6 = 13, 860, X = 0.3 and W = − 0.011 meV (taken from [10,33]). The F 4 and F 6 parameters depend only on J and are, therefore, identical for Er and Dy.
The eigenvalues (ε n ) and eigenvectors ( |ε n ) of the Hamiltonians (1) and (2) are evaluated through a self-consistent procedure to obtain the magnetization components, M R i (i = x, y, z and R = Er, Dy), of each sublattice, given by: where β = 1/k B T and k B is the Boltzmann constant. Taking the crystallographic z-axis as a reference, the angle (ϕ R ) between the magnetic moments and principal crystallographic axes can then be calculated from the magnetization components (3) using the following equation: The magnetic entropy (S mag ) is given by: where R is the gas constant. The magnetic heat capacity is then calculated as C = T ∂S ∂T . The lattice vibrational (Debye model approximation) and electronic contributions to total heat capacity are given by: where Θ D (x, T) and γ(x) are Debye temperature and Sommerfeld coefficient, respectively. They were calculated considering the relative position of a rare-earth element in the lanthanides series LuAl 2 14 and the concentration of each sublattice (Θ D , [43]. The coefficients γ LaAl 2 , γ LuAl 2 Θ LaAl 2 D and Θ LuAl 2 D were taken from the literature [43,44].

Results
The Er 1−x Dy x Al 2 compounds form a continuous solid solution at room temperature, adopting cubic MgCu 2 structure type in the paramagnetic state [32]. Close examination of Bragg peaks in the range of 2θ ≥ 50 • (Mo Kα radiation) confirms that all of the examined samples, including binaries reported in earlier studies (x = 0, 1) [23,24], are cubic (as follows from the absence of splitting) and isostructural at room temperature ( Figure 1a). The lattice parameters at 295 K shown in Table 1 confirm the lattice contraction when Er substitutes Dy.
where Θ , and are Debye temperature and Sommerfeld coefficient, respectively. They were calculated considering the relative position of a rare-earth element in the lanthanides series Θ , , • , and the concentration of each sublattice Θ , Θ , 1 Θ , [43]. The coefficients , Θ and Θ were taken from the literature [43,44].

Results
The Er1−xDyxAl2 compounds form a continuous solid solution at room temperature, adopting cubic MgCu2 structure type in the paramagnetic state [32]. Close examination of Bragg peaks in the range of 2θ ≥ 50° (Mo Kα radiation) confirms that all of the examined samples, including binaries reported in earlier studies (x = 0, 1) [23,24], are cubic (as follows from the absence of splitting) and isostructural at room temperature ( Figure 1a). The lattice parameters at 295 K shown in Table 1 confirm the lattice contraction when Er substitutes Dy. Table 1. Critical temperatures of the heat capacity anomalies [31] and crystallographic parameters of the ErxDy1−xAl2 compounds with x = 0.45, 0.67, and 0.90. T * is the temperature of the second anomaly observed in the heat capacity data. The compositional dependence of crystal structure is more complex at 10 K as shown in Figure  1b, where PXRD patterns are plotted in the same range of Bragg angles. The splitting of different Bragg peaks seen in the binary compounds, ErAl2 and DyAl2, reflects different structural distortions reported earlier. Mixing Er and Dy markedly suppresses those distortions. Thus, Bragg peaks of Er0.67Dy0.33Al2 do not show visible splitting compared to the room temperature pattern, which is in agreement with our earlier study of the neighboring Er0.75Dy0.25Al2 [30], where the cubic structure is preserved down to 5 K. Close examination of the PXRD patterns of the other two samples, Er0.45Dy0.55Al2 and Er0.9Dy0.1Al2, indicates the possibility of much weaker distortions compared to the binaries, seen as a minor but noticeable broadening of the corresponding Bragg peaks.   Table 1. Critical temperatures of the heat capacity anomalies [31] and crystallographic parameters of the Er x Dy 1−x Al 2 compounds with x = 0.45, 0.67, and 0.90. T * is the temperature of the second anomaly observed in the heat capacity data. The compositional dependence of crystal structure is more complex at 10 K as shown in Figure 1b do not show visible splitting compared to the room temperature pattern, which is in agreement with our earlier study of the neighboring Er 0.75 Dy 0.25 Al 2 [30], where the cubic structure is preserved down to 5 K. Close examination of the PXRD patterns of the other two samples, Er 0.45 Dy 0.55 Al 2 and Er 0.9 Dy 0.1 Al 2 , indicates the possibility of much weaker distortions compared to the binaries, seen as a minor but noticeable broadening of the corresponding Bragg peaks.
It is worth noting that the Rietveld refinement of the low-temperature PXRD patterns of all pseudobinary samples studied here can be performed satisfactorily using the cubic symmetry with only a minuscule increase in residuals. For example, for x = 0.9 (Figure 2a), the refinement using the distorted (rhombohedral) structure leads to profile residual, R p = 8.4%, which is only slightly lower than R p = 8.7% for the cubic symmetry ( Figure 2b); the Bragg residuals, R B , are nearly indistinguishable (3.2% vs. 3.1%). Yet, both the visual examination of the Bragg peaks as well as the ability to converge the least-squares refinements using the lower symmetry without imposing constraints on the lattice parameters indicate that the corresponding distortions (rhombohedral for Er 0.9 Dy 0.1 Al 2 and tetragonal for Er 0.45 Dy 0.55 Al 2 , see Tables 1 and 2) do indeed occur in these samples at low temperatures. On the contrary, the refinements of the room temperature or even 100 K data using non-cubic structural models do not converge without constraining lattice parameters. The data with higher resolution, such as x-ray synchrotron radiation, can examine these distortions with a higher level of accuracy but are unlikely to change the main conclusion of this work, which is the much-suppressed distortions due to Er/Dy substitutions.
Metals 2020, 10, x FOR PEER REVIEW 5 of 12 It is worth noting that the Rietveld refinement of the low-temperature PXRD patterns of all pseudobinary samples studied here can be performed satisfactorily using the cubic symmetry with only a minuscule increase in residuals. For example, for x = 0.9 (Figure 2a), the refinement using the distorted (rhombohedral) structure leads to profile residual, Rp = 8.4%, which is only slightly lower than Rp = 8.7% for the cubic symmetry ( Figure 2b); the Bragg residuals, RB, are nearly indistinguishable (3.2% vs. 3.1%). Yet, both the visual examination of the Bragg peaks as well as the ability to converge the least-squares refinements using the lower symmetry without imposing constraints on the lattice parameters indicate that the corresponding distortions (rhombohedral for Er0.9Dy0.1Al2 and tetragonal for Er0.45Dy0.55Al2, see Tables 1 and 2) do indeed occur in these samples at low temperatures. On the contrary, the refinements of the room temperature or even 100 K data using non-cubic structural models do not converge without constraining lattice parameters. The data with higher resolution, such as x-ray synchrotron radiation, can examine these distortions with a higher level of accuracy but are unlikely to change the main conclusion of this work, which is the muchsuppressed distortions due to Er/Dy substitutions.  Some of the examined patterns show strong peak shape anisotropy, which, in our opinion, is not purely instrumental and is likely sample-related. Even the use of multi-parameter axial divergence model, employed by FullProf [38,45], did not allow to appropriately treat the low-angle shape anisotropy for x = 0.45 and x = 0.67 while the peak shape anisotropy of Er0.9Dy0.1Al2 could be accurately refined even by using simple Howard's model [46]. The unaccounted peak shape anisotropy is partially responsible for the elevated values of residuals in the samples with higher Dy content. Additionally, the lack of sample spinning leads to random errors in peak intensities. Below we provide a brief description of the low-temperature structural behavior for each studied sample.
The Er0.9Dy0.1Al2 (Figure 2) shows the lowest Bragg residuals among the three samples. We note  Some of the examined patterns show strong peak shape anisotropy, which, in our opinion, is not purely instrumental and is likely sample-related. Even the use of multi-parameter axial divergence model, employed by FullProf [38,45], did not allow to appropriately treat the low-angle shape anisotropy for x = 0.45 and x = 0.67 while the peak shape anisotropy of Er 0.9 Dy 0.1 Al 2 could be accurately refined even by using simple Howard's model [46]. The unaccounted peak shape anisotropy is partially responsible for the elevated values of residuals in the samples with higher Dy content. Additionally, the lack of sample spinning leads to random errors in peak intensities. Below we provide a brief description of the low-temperature structural behavior for each studied sample. The Er 0.9 Dy 0.1 Al 2 (Figure 2) shows the lowest Bragg residuals among the three samples. We note that the sample may be refined using cubic symmetry down to the lowest temperature (5 K), but below 20 K a rhombohedral distortion (space group R3m) model provides a stable and reproducible solution indicating that it is likely a stable ground state structure for this sample. The temperature dependence of the lattice parameters ( Figure 3) corroborates this suggestion and the degree of the distortion clearly increases on cooling.
Metals 2020, 10, x FOR PEER REVIEW 6 of 12 below ~20 K a rhombohedral distortion (space group 3 ) model provides a stable and reproducible solution indicating that it is likely a stable ground state structure for this sample. The temperature dependence of the lattice parameters ( Figure 3) corroborates this suggestion and the degree of the distortion clearly increases on cooling.

Er0.45Dy0.55Al2
The low-temperature structural behavior of Er0.45Dy0.55Al2 indicates a very weak tetragonal distortion. The X-ray powder diffraction patterns of this compound may be refined in the cubic symmetry at any temperature, but the refinement shows reproducible convergence using a tetragonal (space group I41/amd) model below 30 K (Figure 4a). The tetragonal splitting is noticeable at 10 K pattern (right panel of Figure 4b), and the splitting increases with cooling. The low angle asymmetry in the Bragg peaks of this compound significantly affects the quality of the refinement. The temperature dependence of the lattice parameters calculated using the tetragonal model is shown in Figure 5. The Rietveld refinement of the X-ray powder diffraction pattern of the Er0.45Dy0.55Al2 measured at 10 K using tetragonal model (Rp = 10.7%, RB = 4.8%). A weak Bragg peak at ~11 deg 2θ is external to the sample (acknowledging that it is impossible to identify a phase from a single Bragg peak, its location coincides with the strongest Bragg peak of the low-temperature cubic ice polymorph [47], and this peak appears below 150 K at low pressures, matching the conditions reported in [47]).

Er 0.45 Dy 0.55 Al 2
The low-temperature structural behavior of Er 0.45 Dy 0.55 Al 2 indicates a very weak tetragonal distortion. The X-ray powder diffraction patterns of this compound may be refined in the cubic symmetry at any temperature, but the refinement shows reproducible convergence using a tetragonal (space group I4 1 /amd) model below 30 K (Figure 4a). The tetragonal splitting is noticeable at 10 K pattern (right panel of Figure 4b), and the splitting increases with cooling. The low angle asymmetry in the Bragg peaks of this compound significantly affects the quality of the refinement. The temperature dependence of the lattice parameters calculated using the tetragonal model is shown in Figure 5.
Metals 2020, 10, x FOR PEER REVIEW 6 of 12 below ~20 K a rhombohedral distortion (space group 3 ) model provides a stable and reproducible solution indicating that it is likely a stable ground state structure for this sample. The temperature dependence of the lattice parameters ( Figure 3) corroborates this suggestion and the degree of the distortion clearly increases on cooling.

Er0.45Dy0.55Al2
The low-temperature structural behavior of Er0.45Dy0.55Al2 indicates a very weak tetragonal distortion. The X-ray powder diffraction patterns of this compound may be refined in the cubic symmetry at any temperature, but the refinement shows reproducible convergence using a tetragonal (space group I41/amd) model below 30 K (Figure 4a). The tetragonal splitting is noticeable at 10 K pattern (right panel of Figure 4b), and the splitting increases with cooling. The low angle asymmetry in the Bragg peaks of this compound significantly affects the quality of the refinement. The temperature dependence of the lattice parameters calculated using the tetragonal model is shown in Figure 5. The Rietveld refinement of the X-ray powder diffraction pattern of the Er0.45Dy0.55Al2 measured at 10 K using tetragonal model (Rp = 10.7%, RB = 4.8%). A weak Bragg peak at ~11 deg 2θ is external to the sample (acknowledging that it is impossible to identify a phase from a single Bragg peak, its location coincides with the strongest Bragg peak of the low-temperature cubic ice polymorph [47], and this peak appears below 150 K at low pressures, matching the conditions reported in [47]). (b) The tetragonal splitting of the (008) cubic Bragg reflection. The Rietveld refinement of the X-ray powder diffraction pattern of the Er 0.45 Dy 0.55 Al 2 measured at 10 K using tetragonal model (R p = 10.7%, R B = 4.8%). A weak Bragg peak at~11 deg 2θ is external to the sample (acknowledging that it is impossible to identify a phase from a single Bragg peak, its location coincides with the strongest Bragg peak of the low-temperature cubic ice polymorph [47], and this peak appears below 150 K at low pressures, matching the conditions reported in [47]). (b) The tetragonal splitting of the (008) cubic Bragg reflection.

Er0.67Dy0.33Al2
The behavior of the Er0.67Dy0.33Al2 sample is nearly identical to that reported for the Er0.75Dy0.25Al2 [30]. No indication of a structural distortion was observed in our measurements (Figures 1 and 6). The lattice parameter a steadily decreases on cooling but shows a minor anomaly around TC ~25 K (Figure 7).  The behavior of the Er 0.67 Dy 0.33 Al 2 sample is nearly identical to that reported for the Er 0.75 Dy 0.25 Al 2 [30]. No indication of a structural distortion was observed in our measurements (Figures 1 and 6). The lattice parameter a steadily decreases on cooling but shows a minor anomaly around T C~2 5 K (Figure 7).

Er0.67Dy0.33Al2
The behavior of the Er0.67Dy0.33Al2 sample is nearly identical to that reported for the Er0.75Dy0.25Al2 [30]. No indication of a structural distortion was observed in our measurements (Figures 1 and 6). The lattice parameter a steadily decreases on cooling but shows a minor anomaly around TC ~25 K (Figure 7).

Easy Magnetization Axis as a Function of Composition
Given the experimentally well-established connection between the EMA and the type of structural distortion, one can reasonably conclude that the structural behaviors observed in this study indicate considerable changes in the magnetic anisotropy of the ErxDy1−xAl2 compounds with x. However, due to competition between the EMAs of Er and Dy sublattices, there is no longer a clearly defined easy magnetization axis in these alloys. Further, even in the binary DyAl2, the <111> direction is switched to <100> as temperature varies (our calculations do not differentiate between [100], [010], and [001] directions of the cubic lattice). The evolution of the magnetization angle with x (for 0 1), obtained by the mean-field theory analysis, suggests that for low Er concentrations, i.e., when x < 0.47, the easy direction for both sublattices is <100> ( 90°) in agreement with the DyAl2 ground state EMA (Figure 8). When x increases, the angle decreases for both sublattices changing the easy magnetization direction until the moments align along <111> ( 54.7°), coinciding with ErAl2 easy direction ( Figure 8). Figure 8 shows this behavior for two different temperatures in the ordered region, T = 3 K and 10 K, black and blue curves, respectively. At 3 K, the <100> is the EMA until 0.47, the intermediate state occurs between 0.47 0.76, and for 0.76 the easy magnetization direction is <111>. At 10 K the <100> remains the EMA until 0.7. However, the transition from the <100> towards the <111> direction happens in a narrower concentration range 0.7 0.75. Figure 9 shows as a function of temperature for Er (solid lines) and Dy (dashed lines) sublattices calculated for Er0.45Dy0.55Al2 (a), Er0.67Dy0.33Al2 (b) and Er0.90Dy0.10Al2 (c) compounds. For x = 0.45 the EMA is [100] in almost all temperature range, with a minor (~5 deg) deviation below 3 K. For x = 0.67 (transition region in Figure 8) the sublattices are not aligned in the same direction and there is no clear EMA. The lack of well-defined EMA correlates with the presence of heat capacity anomalies when x = 0.67, 0.75, and 0.82 [32]. Our modeling of heat capacity data supports this hypothesis ( Figure 10). The discontinuity in the entropy in the magnetically ordered region is observed in x = 0.67 but not in x = 0.45 and 0.90 samples. For x = 0.90 the high Er concentration is responsible for the EMA of the compound to be <111>.

Easy Magnetization Axis as a Function of Composition
Given the experimentally well-established connection between the EMA and the type of structural distortion, one can reasonably conclude that the structural behaviors observed in this study indicate considerable changes in the magnetic anisotropy of the Er x Dy 1−x Al 2 compounds with x. However, due to competition between the EMAs of Er and Dy sublattices, there is no longer a clearly defined easy magnetization axis in these alloys. Further, even in the binary DyAl 2 , the <111> direction is switched to <100> as temperature varies (our calculations do not differentiate between [100], [010], and [001] directions of the cubic lattice). The evolution of the magnetization angle with x (for 0 ≤ x ≤ 1), obtained by the mean-field theory analysis, suggests that for low Er concentrations, i.e., when x < 0.47, the easy direction for both sublattices is <100> (ϕ = 90 • ) in agreement with the DyAl 2 ground state EMA (Figure 8). When x increases, the angle ϕ decreases for both sublattices changing the easy magnetization direction until the moments align along <111> (ϕ ≈ 54.7 • ), coinciding with ErAl 2 easy direction ( Figure 8). Figure 8 shows this behavior for two different temperatures in the ordered region, T = 3 K and 10 K, black and blue curves, respectively. At 3 K, the <100> is the EMA until x ≈ 0.47, the intermediate state occurs between 0.47 < x < 0.76, and for x > 0.76 the easy magnetization direction is <111>. At 10 K the <100> remains the EMA until x = 0.7. However, the transition from the <100> towards the <111> direction happens in a narrower concentration range 0.7 < x < 0.75. Figure 9 shows ϕ as a function of temperature for Er (solid lines) and Dy (dashed lines) sublattices calculated for Er 0. 45    , and x = 0.90 (c). The symbols are the experimental measurements reported in [32] and the solid lines represent the theoretical results calculated using the model described above.
The calculation results agree well with the XRD data for the x = 0.45 and 0.90 samples since the calculated EMA matches the type of the observed distortion: rhombohedral for <111> and tetragonal for <100>. The calculations also explain the lack of structural distortion when x = 0.67 since there is no clearly established EMA. At the same time, the magnetic moments in both sublattices are mainly oriented in the xy plane for T = 10 K. We speculate that using high-resolution synchrotron X-ray    , and x = 0.90 (c). The symbols are the experimental measurements reported in [32] and the solid lines represent the theoretical results calculated using the model described above.
The calculation results agree well with the XRD data for the x = 0.45 and 0.90 samples since the calculated EMA matches the type of the observed distortion: rhombohedral for <111> and tetragonal for <100>. The calculations also explain the lack of structural distortion when x = 0.67 since there is no clearly established EMA. At the same time, the magnetic moments in both sublattices are mainly oriented in the xy plane for T = 10 K. We speculate that using high-resolution synchrotron X-ray   , and x = 0.90 (c). The symbols are the experimental measurements reported in [32] and the solid lines represent the theoretical results calculated using the model described above.
The calculation results agree well with the XRD data for the x = 0.45 and 0.90 samples since the calculated EMA matches the type of the observed distortion: rhombohedral for <111> and tetragonal for <100>. The calculations also explain the lack of structural distortion when x = 0.67 since there is no clearly established EMA. At the same time, the magnetic moments in both sublattices are mainly oriented in the xy plane for T = 10 K. We speculate that using high-resolution synchrotron X-ray The symbols are the experimental measurements reported in [32] and the solid lines represent the theoretical results calculated using the model described above.
The calculation results agree well with the XRD data for the x = 0.45 and 0.90 samples since the calculated EMA matches the type of the observed distortion: rhombohedral for <111> and tetragonal for <100>. The calculations also explain the lack of structural distortion when x = 0.67 since there is no clearly established EMA. At the same time, the magnetic moments in both sublattices are mainly oriented in the xy plane for T = 10 K. We speculate that using high-resolution synchrotron X-ray diffraction one may be able to observe a low-symmetry distortion in the ab plane in the x = 0.67 (this work) and x = 0.75 [30] materials.

Conclusions
In the Er x Dy 1−x Al 2 pseudobinary system the structural distortions, clearly observed in the corresponding binary parents, are strongly suppressed. Minor distortions can be detected in the x = 0.45 and 0.9 samples, but at the x = 0.67 concentration no evidence of a structural distortion was observed. The experimental results agree with the mean-field theory modeling that confirms the EMA [100] for x = 0.45 and EMA <111> for x = 0.9 but indicates the absence of a clearly defined EMA (intermediate state with different orientations of sublattices) for x = 0.67. This behavior is contrasting with the structural properties of the Er 1−x Dy x Co 2 alloys, where mixing Er and Dy produces clear and strong low-temperature lattice distortions similar to those observed in the HoCo 2 compound [3]. Given the limited resolution of laboratory PXRD, we acknowledge that the accuracy of our results may be improved by using the synchrotron radiation, but the conclusions presented here will likely stand.