Next Article in Journal
Coexistence of Long-Range Magnetic Order and Magnetic Frustration of a Novel Two-Dimensional S = 1/2 Structure: Na2Cu3(SeO3)4
Next Article in Special Issue
Unraveling the Magnetic Properties of NiO Nanoparticles: From Synthesis to Nanostructure
Previous Article in Journal
Changes in Material Behavior according to the Amount of Recycled Magnetic Materials in Polymer-Bonded Magnets Based on Thermoplastics
Previous Article in Special Issue
Reflection and Refraction of a Spin at the Edge of a Quasi-Two-Dimensional Semiconductor Layer (Quantum Well) and a Topological Insulator
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Phase Diagram Mapping out the Complex Magnetic Structure of Single Crystals of (Gd, Er)B4 Solid Solutions

by
Sueli H. Masunaga
1,2,*,
Vagner B. Barbeta
1,
Fábio Abud
2,3,
Milton S. Torikachvili
4 and
Renato F. Jardim
2
1
Departamento de Física, Centro Universitário FEI, São Bernardo do Campo 09850-901, SP, Brazil
2
Instituto de Física, Universidade de São Paulo, São Paulo 05315-970, SP, Brazil
3
Departamento de Engenharia de Materiais, Escola de Engenharia de Lorena, Universidade de São Paulo, Lorena 12612-550, SP, Brazil
4
Department of Physics, San Diego State University, San Diego, CA 92182, USA
*
Author to whom correspondence should be addressed.
Magnetism 2024, 4(1), 24-34; https://doi.org/10.3390/magnetism4010002
Submission received: 30 November 2023 / Revised: 18 January 2024 / Accepted: 25 January 2024 / Published: 4 February 2024

Abstract

:
Measurements of specific heat and magnetization in single crystals were used to map out the magnetic phase diagram of Gd1−xErxB4 (x = 0.2 and 0.4) solid solutions along the c-axis. While GdB4 orders antiferromagnetically (AF) at 41.7 K, with the easy plane of magnetization oriented perpendicularly to the c-axis, ErB4 displays AF ordering below 15.4 K, with the easy axis along c. Therefore, in solid solutions, the competition between the different spin anisotropies, as well as frustration, lead to a complex spin configuration. These measurements reveal that a 40% substitution of Er for Gd is sufficient for generating a phase diagram similar to the one for the ErB4 system, characterized by the occurrence of plateau phases and other exotic features attributed to the interplay of competing magnetic anisotropies.

Graphical Abstract

1. Introduction

The simultaneous occurrence of geometric frustration, competing interactions, and external magnetic fields can result in the appearance of a rich variety of magnetic states with different spin configurations [1,2,3]. A very well-known geometrical frustrated system is the Shastry–Sutherland lattice, which can be realized in a group of layered compounds [3]. The Shastry–Sutherland model has also been realized in systems displaying magnetization plateaus in the presence of a magnetic field [3,4,5].
Rare-earth tetraborides with the general composition RB4 (R = La-Lu, Pm and Eu excepted) are isostructural, and their crystal structures are tetragonal (space group P4/mbm) [6,7,8]. Most RB4 have an AF ground state [9], and they may be considered as an experimental realization of the Shastry–Sutherland frustrated system model [10,11]. The geometrical frustration in them is due to the presence of two antiferromagnetic exchange interactions: J1 between the nearest neighboring (NN) R atoms and J2 between the next-nearest neighbor (NNN) [11,12]. Although many frustrated compounds are insulating, rare-earth tetraborides are known to be good metals, and their magnetic behavior is determined by Ruderman–Kittel–Kasuya–Yosida (RKKY) interactions [13].
The AF ordering temperature (TN) in the RB4 series ranges from TN∼7 K for R = Ho to 44 K for R = Tb [8]. Some compounds have extra transitions; e.g., R = Nd, Sm, Tb, Dy, Ho, and Tm, and in some cases, these extra transitions are associated with quadrupolar orbital fluctuation [14,15]. ErB4 undergoes an AF transition at TN∼15.4 K, and it has a strong Ising-like anisotropy with the easy axis along the [001] direction. GdB4 displays an AF transition at TN∼42 K, the easy axis is located on the ab-plane, and the system is of Heisenberg-type. Therefore, the magnetic structure in Gd1−xErxB4 mixed compounds should be informed by the competing anisotropies. For magnetic fields parallel or perpendicular to the [001] direction of ErB4 below TN, the M × H curves exhibit a plateau for M = MS/2, where MS is the saturation magnetization [11,16]. On the other hand, the M × H curves for GdB4 below TN do not show plateaus, perhaps due to the small magnetic anisotropy when compared to the other compounds of the RB4 family [17].
The origin of the magnetization plateaus is still not fully understood, though they are thought to result from the interplay between anisotropy, Zeeman energy, and frustration, just as other complex magnetic behaviors. For example, in the case of TmB4, the applied magnetic field induces an ordered phase with stripes reminiscent of domain boundaries [18]. The origin of the “plateau phase” in ErB4 is attributed to the competition between the Zeeman energy and the strong anisotropy [11]. Therefore, the change in the delicate balance of these energies can lead to the appearance of new phases, making the magnetic phase diagram of these compounds rich and complex. Given the sensitivity of the overall magnetic behavior to the competing interactions, a study of Gd1−xErxB4 solid solutions is in order, since it could help to identify the energetics involved in the various magnetic interactions.
In the case of the Gd1−xErxB4 compounds, in addition to the competing anisotropies, which can induce non-conventional spin alignments, the presence of frustration adds a new complication which can result in a complex phase diagram. Systems with competing anisotropies are a matter of intense research [19,20,21]. It has been shown that solid solutions of isostructural materials with orthogonal anisotropies, subjected to a small concentration of one species, can have spins aligned parallel or perpendicular to the magnetic field, depending on the dominant ion and orientation of the magnetic field. On the other hand, in the intermediate range of concentrations, a new mixed phase with spins pointing neither parallel nor perpendicular to the original easy axis can be induced [22].
The goal of this study was to map out the magnetic phases of Gd1−xErxB4 single crystals resulting from the convolution of magnetic frustration and competing anisotropies. High-quality single crystals of Gd1−xErxB4 were grown using the flux method and characterized for single-phase composition, crystal quality, and orientation through Laue X-ray and X-ray powder diffraction methods. These single-crystal samples were studied by means of measurements of DC magnetization (M) and specific heat (Cp) as a function of temperature (T) and magnetic field (H). These experiments were performed using a commercial Physical Property Measurement System (PPMS) from Quantum Design. The phase boundary lines were identified in the H × T phase diagram, and different spin configurations are proposed.

2. Materials and Methods

Single crystals of Gd1−xErxB4 (x = 0.2 and 0.4) were grown in molten aluminum flux using the methodology described in Ref. [23]. The resulting crystals were shaped like rectangular platelets with in-plane (ab-plane) dimensions in the 1–2 mm range and thicknesses in the 0.5–1.0 mm range. X-ray powder (XRD) analyses and X-ray Laue scans indicated the single-phase nature and high crystallinity of these samples. The characterization of Gd1−xErxB4 single crystals has been presented elsewhere [24,25] and is also summarized in the Supplementary Materials.
Specific heat [Cp(T, H)] and magnetization [M(T, H)] measurements, as a function of temperature, applied magnetic field, and orientation, were performed using the dedicated options of the Physical Property Measurement System (PPMS) from Quantum Design, equipped with a 9 T magnet. The calorimeter relies on a relaxation technique. The samples sat on a sapphire stage that had a minute temperature sensor and a heater underneath. A heat pulse was applied, the temperature of the whole assembly as a function of time was monitored, and the heat capacity was extracted from the relaxation rate [26]. The magnetization was measured using a vibrating sample magnetometer (VSM).
Cp(T, H) data were measured as a function of temperature, ranging from 2 K to 100 K, with magnetic fields up to 9 T applied parallel and perpendicular to the c-axis. The M(T) curves were taken under zero-field-cooling (ZFC) and field-cooling (FC) conditions, with magnetic fields ranging from ∼5 mT to 9 T, applied parallel and perpendicular to the c-axis. For the ZFC measurements, a sample was first cooled down to 2 K in zero magnetic field; the magnetic field was subsequently applied, and the M × T data were recorded upon warming up to ∼300 K. The FC data were collected upon cooling to 2 K following the ZFC cycle. M(H) curves were taken at several temperatures ranging from 4 to 30 K and magnetic fields between −9 and +9 T after cooling the crystals in zero applied magnetic field.

3. Experimental Results

Figure 1 shows the specific heat Cp(T, H) for the Gd1−xErxB4 (x = 0.2 and 0.4) compositions. For H = 0 and under the partial substitution of Er3+ for Gd3+, the Néel temperature (TN) decreases from 41.7 K (GdB4) to ~32 K (x = 0.2) and to ~27 K (x = 0.4). For magnetic fields up to 9.0 T, applied perpendicular to the c-axis, TN shifts slightly to lower temperatures for the x = 0.2 and 0.4 compositions, while the overall behavior of Cp(T, H) remains the same, as shown in the lower panels of Figure 1a,b. However, a second shallow maximum centered near 12 K (Tm) develops below TN in the Cp data for H c , in the 2.9–4.5 T field range for x = 0.2, and in the 4.4–6.0 T range x = 0.4, as shown in Figure 1a and Figure 1b, respectively. Small bumps in the Cp data also developed near Tm for the H ab-plane. In the pristine GdB4 compound, this characteristic was attributed to the Schottky anomaly [27], which is independent of the applied magnetic field [28], and to a hidden transition below which an unusually large magnetoresistance is observed [28]. In the case of Gd1−xErxB4, this feature exhibits a slight dependence on the magnetic field for both compositions; as this feature is more responsive to the magnetic field in the x = 0.4 samples, it is possible that it arises from competing anisotropies.
Isothermal magnetization curves M vs. H below TN for H c for the x = 0.2 and x = 0.4 samples are shown in Figure 2 and Figure S1 in the Supplementary Materials. Unlike ErB4, for which the M(H) curves exhibit a magnetization plateau near MS/2 at temperatures as low as 1.5 K [11], this plateau could not be identified in the Gd0.8Er0.2B4 at temperatures as low as 2 K (Figure S1a); only a spin realignment field-induced paramagnetic (FIP) transition could be identified near HFIP ~4 T, where the Gd3+ moments flip away from the ab-plane. The derivative dM/dH (Figure 2) can be used to monitor the FIP transition and to confirm that a plateau phase is not present. The magnetic phase in fields exceeding HFIP probably corresponds to a fully saturated moment J = 15/2 for Er3+. Increasing temperatures softened the FIP transitions, meaning that the corresponding features in the dM/dH curves could no longer be observed for T ≥ 30 K.
On the other hand, the emergence of the plateau phase (PP) can be clearly seen in the x = 0.4 samples for H c between ~4.4 and ~5.5 T and T = 4 K, as shown in Figure S1b. Two features in dM/dH are consistent with the onset and outset of the PP phase, respectively. The transition at ~5.5 T is nearly independent of temperature, and the corresponding feature in dM/dH can be identified up to ~12 K. Higher temperatures result in a less robust PP transition, meaning that above ≈ 14 K, only the FIP transition can be identified. The maxima in the dM/dH curves characterizing the FIP transitions are clearly identifiable in the isotherms for T ≤ 24 K. The minor features observed in the dM/dH curves, close to H = 0 in both samples, are attributed to changes in slope arising from the presence of remanent magnetization and the coercive field within the magnetic hysteresis loops (Figure S2).
Along the c-axis, the magnetization above HFIP increases linearly with the magnetic field for the x = 0.4 samples, and no indication of saturation could be found up to 9 T, the highest field of this experiment. Considering that the Er3+ saturation moment ( μ S ) is 9 μ B and that only 40% of the total magnetization is due to the Er3+, it is expected that Er3+ contributes ≈ 3.6 μ B to the total magnetization. Therefore, the jump observed in magnetization as it transitions to the PP region corresponds to approximately half of the Er3+ ions changing orientation and aligning with the external magnetic field.
The magnetization curves obtained for H c [25] exhibited just a monotonic increase in magnetization with the applied field for the both x = 0.2 and 0.4 samples (Figure S1). This behavior is consistent with the specific heat measurements shown in Figure 1 (lower panels a and b), for which no field-induced phase transition was observed. The reorientation of the spins in GdB4 [17] and the onset of metamagnetic transitions in ErB4 [11] take place under applied magnetic fields perpendicular to the c-axis, typically exceeding 12 T. Given that the highest field in our experiment was 9 T, if a similar behavior took place in the Gd1−xErxB4 compounds, we would not have been able to observe it.
The behavior of M/H × T for the Gd0.8Er0.2B4 and Gd0.6Er0.4B4 single crystals obtained under FC conditions is shown in Figure 3. For low field with H a b -plane and H c , we observed the development of a weak ferromagnetic (wFM) phase, which persisted up to 200 mT. This weak ferromagnetism is quite anisotropic, a behavior reminiscent of RB4 (R = Sm, Gd and Er) when diluted with a nonmagnetic ion [15,24]. The M/H × T data collected under FC and ZFC conditions show significant irreversibility up to ~ 100 mT (Figure S3).
The deconvolution of the components of the wFM can be accomplished by considering that to a first approximation, the magnetization at low fields can be described as follows [24]:
M T = χ 0 ( T ) H + M R ( T ) ,
where χ 0 ( T ) is the magnetic susceptibility of the AF sublattice, and M R ( T ) is the remnant magnetization corresponding to the wFM phase. χ 0 ( T ) was selected from the curves obtained within the 2–3 T range, assuming a negligible contribution from the other phases. The M/H × T curves collapse across the entire temperature range, suggesting the presence of a stable magnetic state. The M R ( T ) curves we obtained are shown in Figure 4.
The M R ( T ) data for both samples (Figure 4) show much higher values along the c-axis, which corresponds to the easy axis in ErB4. The maximum M R values along the c-axis amount to 0.19% and 0.35% of the saturated moments ( μ S ) for the x = 0.2 and 0.4 samples, respectively. In contrast, for the H a b -plane, these values were reduced to 0.05% and 0.007% of μ S for the x = 0.2 and 0.4 samples, respectively. These values are consistent with the findings for the wFM induced by Y substitution in RB4 compounds [15]. In the case of RB4 (R = Sm, Gd, Tb, Dy, and Er) compounds, the highest values of M R were found for ≈ 30% of Y partial substitution in RxYxB4 [15]. For Gd3+ and Er3+, the M R values correspond to 0.42% and 0.21% of the saturated moments, respectively.
The onset of the wFM for H c is defined by the inflexion points taking place at ~32 K and ~27 K for x = 0.2 and 0.4, respectively, both below TN, as shown in the upper panels of Figure 4. This phase emerges at the slightly higher temperatures of ~39 K (x = 0.2) and ~40 K (x = 0.4) for the H a b -plane, values just slightly below TN~42 K for GdB4. The data for H a b -plane suggest that the transitions to the wFM phase take place in two steps, as shown in the bottom panels of Figure 4.
The M × H data for H c in the Gd1−xErxB4 (x = 0.2 and 0.4) compounds reveal remanence and coercivity, as shown in the insets of Figure S1. The magnetic hysteresis for Gd0.8Er0.2B4 persists up to 32 K and 35 K for the H c and H a b -plane, respectively. For Gd0.6Er0.4B4, the magnetic hysteresis can be observed up to 30 K for H c ; however, M R and H C are negligibly small for the H a b -plane.
The Gd0.8Er0.2B4 sample for the field range from 0.5 to 3.0 T ( H c ) as well as Gd0.6Er0.4B4 in the 0.5–4.3 T field range both exhibit M/H × T curves reminiscent of AF systems, displaying a peak near TN. Examples of such AF-like curves are shown in the upper panels of Figure 3 and Figure S4 for both samples. Conversely, when a magnetic field higher than 0.5 T was applied perpendicular to the c-axis, these characteristics of AF systems were observed in both samples, as shown in the lower panels of Figure 3 and Figure S4.
The metamagnetic transitions observed in the M × H curves (Figure 2) can also be identified in the M/H × T curves for H c (upper panels of Figure 3 and Figure S4) within the 3.5–4.5 T range (x = 0.20) and 4.3–6.0 T range (x = 0.4). An increase in magnetization with decreasing temperature can be observed below TN, with the inflection point being identified as the transition temperature. For the x = 0.2 samples, this increase cannot be clearly seen given the logarithmic scale of the graph; however, it is indicated in Figure S4a. In the case of the x = 0.4 samples, the inflection points at Tm, as depicted in the upper panels of Figure 3b and Figure S4b, correspond to the second peak observed in the Cp × T curves, as shown in the upper panel of Figure 1b.
Upon further increasing H c , the M/H × T curves exhibit near saturation at low T, suggesting the complete polarization of Er3+ moments for μ 0 H 5 T and 6 T for x = 0.2 and 0.4, respectively, as depicted in the upper panels of Figure 3 and Figure S4. This finding is consistent with the trend observed in the M × H curves shown in Figure 2, where dM/dH displays a peak just below 5 T and 6 T for x = 0.2 and 0.4, respectively.
For H c exceeding 5 T (x = 0.2) and 6 T (x = 0.4), the magnetization M(H) exhibits linear increases at rates of 0.120 and 0.121 μ B / T , respectively. Remarkably, these values correspond closely to the 0.121 μ B / T slope for the GdB4 measured in Ref. [17]. If this linear trend is maintained, saturation is expected to occur at approximately 50 T and 42 T for x = 0.2 and 0.4, respectively. These values are slightly below the saturation field of 54 T observed in the GdB4 compound, suggesting that the presence of Er3+ in the GdB4 lattice affects the saturation field. It is conceivable that the aligned Er spins might provide an offset field, therefore reducing the applied field necessary to fully saturate the Gd sublattice.
Considering the results described above, the effects of the partial substitution of Er for Gd in Gd1−xErxB4 solid solutions can be summarized as follows: (i) an increase in the magnitude of magnetization, as shown in Figure S4, and in the experimentally determined µeff in the PM state (see, for instance, Reference [25]); (ii) a slight decrease in TN for the stabilization of the AF state; (iii) an appreciable anisotropy of the wFM phase along the c-axis; and (iv) the stabilization of the plateau state at ½ MS for the x = 0.4 crystals. However, the transition to the FIP state, certainly related to the presence of Er3+ ions, occurs at lower applied magnetic fields (HFIP~4 T) for crystals with x = 0.2, whereas it is HFIP~5.5 T for x = 0.4, following the emergence of the plateau state at 4.4 T.
We also emphasize that increasing Er concentration in Gd1−xErxB4 solid solutions is accompanied by a systematic decrease in the lattice parameters due to the smaller ionic radius of Er [25]. As far as this point is concerned, we might recall that the frustration in the Shastry–Sutherland lattice model is closely related to the competition between the Zeeman effect and interactions J1 and J2. An ordered antiferromagnetic ground state has been found for J1/J2 < 1.43 [11]. In our solid solutions, due to the partial substitution of smaller Er ions, such a J1/J2 ratio is altered and maybe a source of the appearance of the ½ MS plateau in the M × H curves.
We also point out that the pristine ErB4 compound exhibits a strong Ising-type single-ion anisotropy with an easy magnetization axis along the c-axis [4]. On the other hand, the single-ion anisotropy of Gd is considered small, and GdB4 can be described as a noncollinear Heisenberg antiferromagnet with the easy axis along the (110) direction [29]. Therefore, the partial substitution of Er by Gd has an impact on the magnetic anisotropy of the Gd1−xErxB4 solid solution. Such a little but important variation in ionic radii, along with the difference in the magnetic moments between Gd and Er, may, in turn, influence the geometrically frustrated exchange interaction and the exchange couplings, with these effects playing fundamental roles in stabilizing distinct magnetic states [4]. This point is currently under investigation.

4. Discussion

Informed by the detailed magnetization data on the Gd0.8Er0.2B4 and Gd0.6Er0.4B4 single crystals, we were able to assemble field-temperature (H × T) phase diagrams for H c , as shown in Figure 5. The value of TN decreases slightly for magnetic field H a b -plane. Given that 9 T is the highest field of this set of data, the phase diagram misses the possible field-induced phase transition at higher fields. Consequently, the H × T phase diagram for both samples in the ab-plane only includes vertical boundaries delineating paramagnetic (PM), wFM, and AF regions.
Regions I and II of the phase diagrams in Figure 5 correspond to the paramagnetic PM and wFM phases for Gd0.8Er0.2B4 and Gd0.6Er0.4B4. However, region II for x = 0.2 and 0.4 are clearly different. For the x = 0.2 samples, a minor anomaly in the Cp × T data up to μ 0 H = 9 T is observed just above TN. This anomaly coincides with the onset of the wFM phase, as indicated by the vertical borderline at ~ 34 K. On the other hand, for the x = 0.4 samples, there is no corresponding anomaly in the Cp(H, T) data associated with the wFM onset. Our Cp(H, T) and M(H, T) data lack the resolution required to definitively identify the origin of the observed differences in region II for the two Er compositions.
The AF phase was observed to occur within regions III and IV of the phase diagrams for both x = 0.2 and 0.4 compositions. The phase boundary between regions I/II and III/IV are clearly correlated with anomalies in the Cp × T data in applied magnetic fields up to 9 T. In contrast, the ordering temperature is suppressed with increasing magnetic field, as shown in the M/H × T data, in which saturation in M was observed at low temperatures. Specifically, TN could be reliably, unambiguously determined from the M/H × T curves in magnetic fields up to ~4.5 T and ~6.0 T for the x = 0.2 and 0.4 compositions, respectively. Region III represents the AF phase of the Gd1−xErxB4 system, whereas region IV denotes the coexistence of the AF phase, associated with the Gd3+ sublattice and the FIP phase of the Er3+ ions.
The phase boundary between regions III and IV was well defined by the observation of well-defined magnetization jumps observed in the M × H isotherms for H c . Measurements of Cp × T and M/H × T provide limited insight into this boundary due to the small variations in this boundary with H. Phase V (PP), resulting from the plateau at ½ MS of the Er+3 ions, was clearly observed only for the x = 0.4 composition. The transitions between phase III and V, as well as between phase V and IV, are characterized by the magnetization discontinuities observed in the M × H isotherms. The boundary between phases III and V was also clearly delineated through Cp × T characterization.
Establishing whether the presence of phase V for the x = 0.2 composition results from the possible realignment of some the Er3+ spins is quite difficult. In this sample, it is anticipated that the Er sublattice contributes 1.8 μ B to the overall magnetization of the system. This magnitude is consistent with the observed jump in the M × H curves, a feature which is very pronounced at T = 4 K. As a result, the system appears to transition directly from the AF phase to the mixed phase IV without an intermediate phase V (PP). The temperature of the triple point, in this case, which delineates the boundary lines of phases II, III, and IV, appears to be near 32 K and 3 T.
The phase diagram highlights the role of the partial Er substitution in the GdB4 lattice and its impact on the boundary lines within the system. For instance, a partial substitution of only 40% of Er for Gd reveals a striking resemblance between the resulting phase diagram (Figure 5b) and that of the ErB4 compound [11]. Analogous to ErB4, the x = 0.4 composition displays the presence of two discernible metamagnetic phase transitions, albeit at higher fields. In the ErB4 phase diagram, derived from neutron diffraction data, three distinct regions identified as being antiferromagnetic, ferrimagnetic, and ferromagnetic merge to a triple point at 11.7 K and 2.1 T [30]. In the Gd0.6Er0.4B4 composition, this triple point is located at ~12.8 K, nearly the same temperature as that of ErB4, but at a higher magnetic field of ~5.4 T, where the phases in region III, IV, and V coalesce, as shown in Figure 5b. Additionally, another triple point can be defined at ~24 K and ~6 T, where phases I, II, and IV converge.
Remarkably, a partial substitution of just 20% of Er is sufficient for a metamagnetic transition from AF to FIP, as shown in Figure 5a. However, the FIP phase corresponds solely to the orientation of the Er3+ ions, as the magnetization corresponds to only 26% and 44% of the theoretical magnetization of μ S = 7.4 and 7.8 μ B expected for the x = 0.2 and 0.4 compositions, respectively. The complete saturation of the Gd1−xErxB4 compounds can be achieved only at saturation fields close to the values found for GdB4 (i.e., 50 T and 42 T, respectively). In the case of GdB4, this metamagnetic transition occurs only at the very high fields of ~52 T and ~54 T, applied along the a and c axes, respectively [17]. It is important to emphasize that both the Gd0.8Er0.2B4 and Gd0.6Er0.4B4 compositions exhibit pronounced anisotropies, while the preferred orientation for ErB4 is the c-axis.
Our combined analysis of the magnetization and specific heat data suggest that the Er and Gd sublattices keep their original easy directions (i.e., the ab-plane for Gd and the c-axis for Er). The possibility of a small canting angle for both easy directions, which would also explain the appearance of the wFM, cannot be completely ruled out.

5. Conclusions

In summary, we relied on the magnetization and specific heat data for Gd0.8Er0.2B4 and Gd0.6Er0.4B4 to construct a magnetic phase diagram detailing a set of complex magnetic structures. Notably, for both samples, the Néel temperature drops slightly for the H a b -plane, while the occurrence of a field-induced plateau phase in fields below 9 T ( H c ) is specific to the x = 0.4 compositions. Therefore, the phase diagram focuses primarily on the emergence of phases for H c . We have identified five different regions in the phase diagrams. Regions I to IV for x = 0.2 comprise paramagnetic, weak ferromagnetic, antiferromagnetic, and mixed field-induced paramagnetic and antiferromagnetic phases. A plateau phase composed of the ½ MS plateau state (PP) was found for the x = 0.4 compositions.
The phase diagrams displayed in Figure 5 summarize the effect of the partial substitution of Er for Gd in the GdB4 lattice, and the boundary lines were experimentally determined along the c-axis using the M × T, M × H, and Cp × T curves. In the Gd1−xErxB4 compounds, the exchange couplings J1 and J2, which play an essential role in stabilizing distinct magnetic states, were altered due to the differences between Gd and Er ionic radii, as well as differences in magnetic moment and single-ion anisotropy. Contrary to classical competing anisotropy systems, such as Fe1−xCoxCl2 [31], FexNi1−xF2 [32], and others, it is important to point out that a tiny (~20%) partial Er substitution is sufficient for inducing a metamagnetic transition from the antiferromagnetic state to the field-induced paramagnetic state. Also, we observed that a 40% Er substitution results in a phase diagram closely resembling that of ErB4. Both compositions studied exhibit strong anisotropies favoring the c-axis, the easy axis for ErB4. These findings provide insights into the interplay of the magnetic phases as a c-axis easy-axis ion partially replaces an ab-plane easy-orientation ion.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/magnetism4010002/s1, magnetism–Supplementary Materials: Figure S1. M × H curves under H c and H c for single crystals of Gd1−xErxB4 with (a) x = 0.2 at T = 2 K; (b) x = 0.4 at T = 4 K. The arrows indicate the sweeping direction of the applied magnetic field. The insets show an expanded view of the low-field limits of the curves; Figure S2. M × H under H c for single crystals of Gd1−xErxB4 with (a) x = 0.2 at T = 20 K; (b) x = 0.4 at T = 14 K. The corresponding derivatives, dM/dH, for both compositions are presented at the bottom of the plots; Figure S3. M/H × T for Gd1−xErxB4 samples with (a) x = 0.2 and (b) x = 0.4 measured under magnetic fields up to 0.1 T. The closed symbols represent zero-field-cooled (ZFC) curves, while the open symbols depict field-cooled (FC) curves, as indicated by the arrows; Figure S4. M/H × T for Gd1−xErxB4 samples with (a) x = 0.2 and (b) x = 0.4 measured under magnetic fields ranging from 0.5 T to 9.0 T in the FC process. The shallow maxima observed at Tm in the Cp(T) curves correspond to the inflexion points in the M/H × T curves (marked by the arrows in the figure).

Author Contributions

Conceptualization, R.F.J.; formal analysis, S.H.M. and V.B.B.; investigation, S.H.M., F.A., R.F.J. and M.S.T.; resources, R.F.J.; writing—original draft preparation, S.H.M. and V.B.B.; writing—review and editing, F.A., R.F.J. and M.S.T.; supervision, R.F.J.; funding acquisition, R.F.J. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by FAPESP (Grants No. 2013/07296–2, 2019/26141–6, 2022/02691–0, and 2022/10874–7), CNPq (Grant No. 301463/2019–0), and CAPES (Finance Code 001).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data presented in this study are available upon request.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Seabra, L.; Momoi, T.; Sindzingre, P.; Shannon, N. Phase diagram of the classical Heisenberg antiferromagnet on a triangular lattice in an applied magnetic field. Phys. Rev. B 2011, 84, 214418. [Google Scholar] [CrossRef]
  2. Gabani, S.; Flachbart, K.; Siemensmeyer, K.; Mori, T. Magnetism and superconductivity of rare earth borides. J. Alloys Compd. 2020, 821, 153201. [Google Scholar] [CrossRef]
  3. Brunt, D.C.; Balakrishnan, G.; Mayoh, D.; Lees, M.R.; Gorbunov, D.; Qureshi, N.; Petrenko, O. Magnetisation process in the rare earth tetraborides, NdB4 and HoB4. Sci. Rep. 2018, 8, 232. [Google Scholar]
  4. Farkašovský, P.; Regeciová, L. Formation of Magnetization Plateaus in Rare Earth Tetraborides: Exact Diagonalization and Quantum Monte Carlo Studies. J. Supercond. Nov. Magn. 2020, 33, 3463–3467. [Google Scholar] [CrossRef]
  5. Pristáš, G.; Bačkai, J.; Orendáč, M.; Gabáni, S.; Gažo, E.; Shitsevalova, N.; Dukhnenko, A.; Siemensmeyer, K.; Flachbart, K. Magnetic Phase Diagram of Geometrically Frustrated Magnetic System ErB4 under High Pressure. In Proceedings of the International Conference on Strongly Correlated Electron Systems (SCES2019), Okayama, Japan, 23–28 September 2019. [Google Scholar]
  6. Fisk, Z.; Maple, M.B.; Johnston, D.C.; Woolf, L.D. Multiple phase transitions in rare earth tetraborides at low temperature. Solid State Commun. 1981, 39, 1189–1192. [Google Scholar] [CrossRef]
  7. Etourneau, J.; Mercurio, J.P.; Berrada, A.; Hagenmuller, P.; Georges, R.; Bourezg, R.; Gianduzzo, J.C. The magnetic and electrical properties of some rare earth tetraborides. J. Common Met. 1979, 67, 531–539. [Google Scholar] [CrossRef]
  8. Yin, Z.P.; Pickett, W.E. Rare-earth–boron bonding and 4f state trends in RB4 tetraborides. Phys. Rev. B 2008, 77, 035135. [Google Scholar] [CrossRef]
  9. Buschow, K.H.J. Magnetic Properties of Borides. In Boron and Refractory Borides; Matkovich, V.I., Ed.; Springer: Berlin/Heidelberg, Germany, 1977; pp. 494–515. ISBN 978-3-642-66620-9. [Google Scholar]
  10. Watanuki, R.; Sato, G.; Suzuki, K.; Ishihara, M.; Yanagisawa, T.; Nemoto, Y.; Goto, T. Geometrical Quadrupolar Frustration in DyB4. J. Phys. Soc. Jpn. 2005, 74, 2169–2172. [Google Scholar] [CrossRef]
  11. Michimura, S.; Shigekawa, A.; Iga, F.; Sera, M.; Takabatake, T.; Ohoyama, K.; Okabe, Y. Magnetic frustrations in the Shastry–Sutherland system ErB4. Phys. B Condens. Matter 2006, 378–380, 596–597. [Google Scholar] [CrossRef]
  12. Inami, T.; Ohwada, K.; Matsuda, H.Y.; Wen, O.Z.; Nojiri, H.; Matsumura, T.; Okuyama, D.; Murakami, Y. Resonant Magnetic X-ray Diffraction Study on Successive Metamagnetic Transitions in TbB4. J. Phys. Soc. Jpn. 2009, 78, 033707. [Google Scholar] [CrossRef]
  13. Schäfer, W.; Will, G.; Buschow, K.H.J. The magnetic structure of the rare earth tetraborides ErB4 and DyB4. J. Chem. Phys. 1976, 64, 1994–1997. [Google Scholar] [CrossRef]
  14. Okuyama, D.; Matsumura, T.; Nakao, H.; Murakami, Y. Quadrupolar Frustration in Shastry–Sutherland Lattice of DyB4 Studied by Resonant X-ray Scattering. J. Phys. Soc. Jpn. 2005, 74, 2434–2437. [Google Scholar] [CrossRef]
  15. Kang, B.Y.; Lee, S.S.; Song, M.S.; Cho, K.K.; Han, S.H.; Cho, B.K. Universality of weak ferromagnetism in the magnetically frustrated system R1−xYxB4. Curr. Appl. Phys. 2016, 16, 1001–1004. [Google Scholar] [CrossRef]
  16. Mat’aš, S.; Siemensmeyer, K.; Wheeler, E.; Wulf, E.; Beyer, R.; Hermannsdörfer, T.; Ignatchik, O.; Uhlarz, M.; Flachbart, K.; Gabáni, S.; et al. Magnetism of rare earth tetraborides. J. Phys. Conf. Ser. 2010, 200, 032041. [Google Scholar] [CrossRef]
  17. Kikkawa, A.; Katsumata, K.; Narumi, Y.; Suga, K.; Fukui, T.; Sugaya, T.; Kindo, K.; Iga, F.; Michimura, S. Magnetization Process in GdB4. J. Phys. Soc. Jpn. 2007, 76, 024711. [Google Scholar] [CrossRef]
  18. Siemensmeyer, K.; Wulf, E.; Mikeska, H.-J.; Flachbart, K.; Gabáni, S.; Mat’aš, S.; Priputen, P.; Efdokimova, A.; Shitsevalova, N. Fractional Magnetization Plateaus and Magnetic Order in the Shastry-Sutherland Magnet TmB4. Phys. Rev. Lett. 2008, 101, 177201. [Google Scholar] [CrossRef]
  19. Katsumata, K.; Tawaraya, T. Tetracritical behaviours of random mixtures with competing spin anisotropies. J. Magn. Magn. Mater. 1980, 15–18, 253–254. [Google Scholar] [CrossRef]
  20. Mano, H. Pair-Spin Theory of a Random Mixture with Competing Ising and XY Anisotropies. J. Phys. Soc. Jpn. 1986, 55, 908–919. [Google Scholar] [CrossRef]
  21. Katsumata, K.; Shapiro, S.M.; Matsuda, M.; Shirane, G.; Tuchendler, J. Simultaneous ordering of orthogonal spin components in a random magnet with competing anisotropies. Phys. Rev. B 1992, 46, 14906–14908. [Google Scholar] [CrossRef]
  22. Mano, H. Possible Phase Diagrams of Systems with Competing Anisotropies. Prog. Theor. Phys. Suppl. 1990, 101, 597–610. [Google Scholar] [CrossRef]
  23. Fisk, Z.; Cooper, A.S.; Schmidt, P.H.; Castellano, R.N. Preparation and lattice parameters of the rare earth tetraborides. Mater. Res. Bull. 1972, 7, 285–288. [Google Scholar] [CrossRef]
  24. Masunaga, S.H.; Barbeta, V.B.; Jardim, R.F.; Becerra, C.C.; Torikachvili, M.S.; Rosa, P.F.S.; Fisk, Z. Anomalous remnant magnetization in dilute antiferromagnetic Gd1−xYxB4. Phys. Rev. Mater. 2018, 2, 084415. [Google Scholar] [CrossRef]
  25. Masunaga, S.H.; Barbeta, V.B.; Abud, F.; Torikachvili, M.S.; Jardim, R.F. Anomalous Ferromagnetic Phase in the Gd1−xErxB4 Series: Crystal Growth, Thermal, and Magnetic Properties. Crystals 2023, 13, 1137. [Google Scholar] [CrossRef]
  26. Hwang, J.S.; Lin, K.J.; Tien, C. Measurement of heat capacity by fitting the whole temperature response of a heat-pulse calorimeter. Rev. Sci. Instrum. 1997, 68, 94–101. [Google Scholar] [CrossRef]
  27. Novikov, V.V.; Mitroshenkov, N.V.; Morozov, A.V.; Matovnikov, A.V.; Avdashchenko, D.V. Heat capacity and thermal expansion of gadolinium tetraboride at low temperatures. J. Appl. Phys. 2012, 111, 063907. [Google Scholar] [CrossRef]
  28. Cho, B.K.; Rhyee, J.-S.; Kim, J.Y.; Emilia, M.; Canfield, P.C. Anomalous magnetoresistance at low temperatures (T ≤ 10K) in a single crystal of GdB4. J. Appl. Phys. 2005, 97, A923. [Google Scholar] [CrossRef]
  29. Johnston, D.C. Magnetic Susceptibility of Collinear and Noncollinear Heisenberg Antiferromagnets. Phys. Rev. Lett. 2012, 109, 077201. [Google Scholar] [CrossRef] [PubMed]
  30. Pfeiffer, F.; Schäfer, W.; Will, G.; Etourneau, J.; Georges, R. The magnetic phase diagram of ErB4. J. Magn. Magn. Mater. 1979, 14, 306–308. [Google Scholar] [CrossRef]
  31. Wong, P. Specific-Heat Study of Random-Field and Competing-Anisotropy Effects in Fe1−xCoxCl2. Phys. Rev. B 1986, 34, 1864–1879. [Google Scholar] [CrossRef]
  32. Perez, F.A.; Borisov, P.; Johnson, T.A.; Stanescu, T.D.; Trappen, R.; Holcomb, M.B.; Lederman, D.; Fitzsimmons, M.R.; Aczel, A.A.; Hong, T. Phase Diagram of a Three-Dimensional Antiferromagnet with Random Magnetic Anisotropy. Phys. Rev. Lett. 2015, 114, 097201. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Cp × T curves under different applied magnetic fields H for Gd1−xErxB4 single crystals. (a) x = 0.2; (b) x = 0.4. Besides the peak at TN associated with the AF transition, a shallow maximum is also observed at lower temperatures (Tm ~12 K) along the c-axis.
Figure 1. Cp × T curves under different applied magnetic fields H for Gd1−xErxB4 single crystals. (a) x = 0.2; (b) x = 0.4. Besides the peak at TN associated with the AF transition, a shallow maximum is also observed at lower temperatures (Tm ~12 K) along the c-axis.
Magnetism 04 00002 g001
Figure 2. M × H curves for single crystals of Gd1−xErxB4 with (a) x = 0.2 and temperatures T = 4, 10, 12, 15, 20, 25, and 30 K and (b) x = 0.4 at T = 4, 6, 8, 10, 12, 14, 16, 18, and 20 K. In both cases, the magnetic field H was applied along the c direction. The derivatives dM/dH for two compositions, ranging from −9 T to 9 T, are shown at the bottom of the plots.
Figure 2. M × H curves for single crystals of Gd1−xErxB4 with (a) x = 0.2 and temperatures T = 4, 10, 12, 15, 20, 25, and 30 K and (b) x = 0.4 at T = 4, 6, 8, 10, 12, 14, 16, 18, and 20 K. In both cases, the magnetic field H was applied along the c direction. The derivatives dM/dH for two compositions, ranging from −9 T to 9 T, are shown at the bottom of the plots.
Magnetism 04 00002 g002
Figure 3. M/H × T for Gd1−xErxB4 samples with (a) x = 0.2 and (b) x = 0.4 for magnetic fields up to 9 T. The vertical scales are logarithmic to emphasize the field dependence of M/H.
Figure 3. M/H × T for Gd1−xErxB4 samples with (a) x = 0.2 and (b) x = 0.4 for magnetic fields up to 9 T. The vertical scales are logarithmic to emphasize the field dependence of M/H.
Magnetism 04 00002 g003
Figure 4. MR × Τ for single crystals of Gd1−xErxB4 with (a) x = 0.2 and (b) x = 0.4 under different external magnetic fields H. The upper panels show the results when the magnetic field H is parallel to the c-axis, and the lower panels show the results for when the magnetic field H was perpendicular to the c-axis.
Figure 4. MR × Τ for single crystals of Gd1−xErxB4 with (a) x = 0.2 and (b) x = 0.4 under different external magnetic fields H. The upper panels show the results when the magnetic field H is parallel to the c-axis, and the lower panels show the results for when the magnetic field H was perpendicular to the c-axis.
Magnetism 04 00002 g004
Figure 5. Magnetic phase diagrams for (a) Gd0.8Er0.2B4 and (b) Gd0.6Er0.4B4 for H c . The characteristics of the different phases are discussed in the text. Regions I, II, III, IV, and V correspond to PM, wFM, AF, (FIP + AF), and PP, respectively.
Figure 5. Magnetic phase diagrams for (a) Gd0.8Er0.2B4 and (b) Gd0.6Er0.4B4 for H c . The characteristics of the different phases are discussed in the text. Regions I, II, III, IV, and V correspond to PM, wFM, AF, (FIP + AF), and PP, respectively.
Magnetism 04 00002 g005
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Masunaga, S.H.; Barbeta, V.B.; Abud, F.; Torikachvili, M.S.; Jardim, R.F. Phase Diagram Mapping out the Complex Magnetic Structure of Single Crystals of (Gd, Er)B4 Solid Solutions. Magnetism 2024, 4, 24-34. https://doi.org/10.3390/magnetism4010002

AMA Style

Masunaga SH, Barbeta VB, Abud F, Torikachvili MS, Jardim RF. Phase Diagram Mapping out the Complex Magnetic Structure of Single Crystals of (Gd, Er)B4 Solid Solutions. Magnetism. 2024; 4(1):24-34. https://doi.org/10.3390/magnetism4010002

Chicago/Turabian Style

Masunaga, Sueli H., Vagner B. Barbeta, Fábio Abud, Milton S. Torikachvili, and Renato F. Jardim. 2024. "Phase Diagram Mapping out the Complex Magnetic Structure of Single Crystals of (Gd, Er)B4 Solid Solutions" Magnetism 4, no. 1: 24-34. https://doi.org/10.3390/magnetism4010002

APA Style

Masunaga, S. H., Barbeta, V. B., Abud, F., Torikachvili, M. S., & Jardim, R. F. (2024). Phase Diagram Mapping out the Complex Magnetic Structure of Single Crystals of (Gd, Er)B4 Solid Solutions. Magnetism, 4(1), 24-34. https://doi.org/10.3390/magnetism4010002

Article Metrics

Back to TopTop