Setting the Basis for the Interpretation of Temperature First Order Reversal Curve (TFORC) Distributions of Magnetocaloric Materials

: First Order Reversal Curve (FORC) distributions of magnetic materials are a well-known tool to extract information about hysteresis sources and magnetic interactions, or to ﬁngerprint them. Recently, a temperature variant of this analysis technique (Temperature-FORC, TFORC) has been used for the analysis of the thermal hysteresis associated with ﬁrst-order magnetocaloric materials. However, the theory supporting the interpretation of the diagrams is still lacking, limiting TFORC to a ﬁngerprinting technique so far. This work is a ﬁrst approach to correlate the modeling of ﬁrst-order phase transitions, using the Bean–Rodbell model combined with a phenomenological transformation mechanism, with the features observed in experimental TFORC distributions of magnetocaloric materials. The di ﬀ erent characteristics of the transformations, e.g., transition temperatures, symmetry, temperature range, etc., are correlated to distinct features of the distributions. We show a catalogue of characteristic TFORC distributions for magnetocaloric materials that exhibit some of the features observed experimentally.


Introduction
Magnetocaloric (MC) materials [1], those that experience a significant temperature change when submitted to magnetic field variations in adiabatic conditions, deserve the attention of the magnetic research community due to their implementation on a new, green and energy efficient solid-state cooling technology at room temperature [2,3]. The effect is associated with significant magnetization changes driven by temperature and magnetic field, classifying the MC materials according to the order of their thermomagnetic phase transition: first-and second-order type [4]. Since the first experiments with pure Gd in 1976 [5], which undergoes a Curie transition (second-order type) from the ferromagnetic (FM) to the paramagnetic (PM) state, the interest has moved nowadays to those magnetocaloric materials which undergo magnetoelastic or magnetostructural transitions (first-order type) such as GdSiGe [6,7], LaFeSi [8,9], MnFePAs [10,11] or Heusler [12,13]-based alloys. Typically, these materials present higher MC response than those with second order transitions, although they have associated undesired thermal and magnetic hysteresis, which significantly reduces its response in cyclical conditions, limiting their applicability in refrigeration devices [14]. Therefore, the understanding of the different sources and mechanisms of hysteresis is a relevant topic in the MC research field [15].
First Order Reversal Curve (FORC) distributions of magnetic materials are a well-recognized tool to extract information about magnetic interactions or to fingerprint magnetic materials [16][17][18]. Recently, the possibility to use this technique to analyze the thermal hysteresis associated with the first-order magnetocaloric materials has been shown [19,20]. This so called TFORC technique (which is Metals 2020, 10 formally defined analogously to its more conventional field counterpart) is based on the analysis of the minor loops of the magnetization (M) in the hysteretic region, defining a certain set of reversal curves, in an analogous way to what was previously done for spin crossover materials [21]. To measure these curves, first the system is heated up (cooled down) from saturation to a certain temperature (reversal temperature, T r ) and, next, the system is cooled down (heated up) to the initial saturated value while the magnetization values are measured. The set of magnetization curves for different T r are processed to obtain a TFORC distribution in the case of heating and cooling according to: where the subindex C and H refer to the cooling and heating curves, respectively. These distributions are explained in terms of hysteron distributions although a direct connection between the FORC and TFORC distributions and the Preisach model [22] is not usually found. The hysteron is a fundamental hysteresis loop, which is characterized by two switching fields between two stable states called "up" and "down" (the fundamental hysteresis unit) and is defined by two parameters for the thermal hysterons: (1) the hysteron width (T h ) and (2) the central position of the hysteron (T u ). The distributions shown in Equations (1) and (2) can be obtained as a function of these two temperatures by a simple variable change: T h = |T − T r |/2 and T u = (T + T r )/2. In the case of conventional FORC, the extensive modelling efforts performed over decades have allowed the extraction of different characteristics of the material from the features found in the distribution [18,[23][24][25][26][27]. However, the recent appearance of TFORC in magnetocaloric materials means that there is no fundamental basis that connects TFORC diagrams and the underlying physics of the thermal hysteresis, limiting TFORC to a fingerprinting technique so far. In this work, we set the basis for the interpretation of TFORC diagrams by performing a systematic study based on modeling of first-order phase transitions using the Bean-Rodbell model [28] extended with a phenomenological transformation mechanism and correlating the obtained diagrams with the characteristics of the transition. The use of skewed normal distributions for modeling the transformation process between FM and PM phases allows the calculation of the reversal magnetization curves, unlike with the pure Bean and Rodbell model. The transformation characteristics explored are: transition temperatures, temperature ranges of the transformations and symmetry/asymmetry of cooling and heating processes. Each of them produce different shapes in the distributions, such as circular, ellipsoidal, semicircular or triangular. This first catalogue of images obtained in this work, analogous to those catalogues available for conventional FORC, can help the interpretation of experimental TFORC distributions.

Methods
The Bean-Rodbell model was proposed in 1962, to describe the first-order magnetoelastic transition observed in the MnAs compound [28]. Apart from this material, the model has been successfully applied to other materials that undergo a magnetoelastic first-order transition such as LaFeSi [29,30] and MnFePAs [31,32] alloys, which are of great technological interest. The model assumes that magnetic interactions are affected by such volume changes as: where T t is the transition temperature, T 0 the transition temperature in absence of deformations (or Curie temperature), β an introduced parameter and ω the relative volume change.
Introducing the assumed dependence shown in Equation (3) into the expression for the Gibbs energy (which considers the Zeeman, exchange, distortion, external pressure and entropy terms in the framework of mean field theories [28]) and minimizing with respect to ω and M variables, the following expression for M can be obtained: where M S is the saturation magnetization, B J the Brillouin function, g the Landé factor, µ B the Bohr magneton, J the total quantum angular momentum, µ 0 the vacuum permeability, H the magnetic field and η a parameter of the model, defined as: where k is the compressibility. The η parameter controls the order of the transition, being first-order if η > 1, second-order if η < 1 and η = 1 corresponds to the critical point (where the order of the transition changes from second to first). T has a single-valued solution as a function of M and H in Equation (4). For η > 1, metastable and instable regions appear in the solution. The stable solution is obtained when the energy of FM and PM phases are the same (being T t,ST , the transition temperature of the stable branch), which is calculated according to the so-called equal-area construction from the single-valued solution. The existence of the metastable branches leads to two different solutions when the material is under heating or cooling conditions (T H t,MS and T C t,MS ), which gives rise to associated hysteresis. For simulating the Bean-Rodbell model, the following parameters have been used in this work: g = 2, J = 7/2, M S = 1.96 × 10 6 A·m −1 , T 0 = 300 K, k = 2.5·× 10 −11 Pa −1 and β = 13.54 (η = 1.5). These parameters are arbitrarily chosen to represent a first-order phase transition material with a transition temperature close to room temperature. Figure 1 shows the temperature dependence of the Bean-Rodbell solution for a magnetic field of 0.01 T.
Metals 2020, 10, x FOR PEER REVIEW 3 of 16 the framework of mean field theories [28]) and minimizing with respect to and variables, the following expression for can be obtained: where is the saturation magnetization, the Brillouin function, the Landé factor, the Bohr magneton, the total quantum angular momentum, the vacuum permeability, the magnetic field and a parameter of the model, defined as: where is the compressibility. The parameter controls the order of the transition, being firstorder if > 1, second-order if < 1 and = 1 corresponds to the critical point (where the order of the transition changes from second to first).
has a single-valued solution as a function of and in Equation (4). For > 1, metastable and instable regions appear in the solution. The stable solution is obtained when the energy of FM and PM phases are the same (being , , the transition temperature of the stable branch), which is calculated according to the so-called equal-area construction from the single-valued solution. The existence of the metastable branches leads to two different solutions when the material is under heating or cooling conditions ( , and , ), which gives rise to associated hysteresis. For simulating the Bean-Rodbell model, the following parameters have been used in this work: = 2, = 7/2, = 1.96 × 10 6 A·m −1 , = 300 K, = 2.5·× 10 −11 Pa −1 and = 13.54 ( = 1.5). These parameters are arbitrarily chosen to represent a first-order phase transition material with a transition temperature close to room temperature. Figure 1 shows the temperature dependence of the Bean-Rodbell solution for a magnetic field of 0.01 T. In the framework of the pure Bean-Rodbell model, the transformation between both FM and PM phases is described as just an abrupt step in temperature (or magnetic field) between the values In the framework of the pure Bean-Rodbell model, the transformation between both FM and PM phases is described as just an abrupt step in temperature (or magnetic field) between the values corresponding to both phases. This feature does not agree with the observations of first-order transitions in which coexistence between phases is intrinsic to the nature of the phase transition (and crucial for understanding FORC distributions). To reproduce the transformed fraction ( f ) through the first-order transition, we propose the modeling through the cumulative distribution function of a skewed normal distribution [33] as: where erf is the error function and O T is the Owen's T function. The different parameters are ξ which is the location, σ the scale and α the shape or skewness factor (which can be either positive or negative). For this distribution, the mean value is ξ + σ α √ 1+α 2 2 π and the variance is σ 2 (1 − 2α 2 π(1+α 2 ) ). It is observed that for α = 0, the normal distribution function is recovered (as O T (α = 0 ) = 0). This phenomenological skewed normal distribution function can be linked to different transformation characteristics of the materials, such as nucleation-growth processes, transition temperature distributions (caused by compositional inhomogeneities) or internal stress effects, among others. Some of these details of the transformation can be directly identified (as will be mentioned later); for others (e.g., to distinguish between compositional inhomogeneities or kinetic processes), a detailed analysis of the samples using complementary techniques would be required. To be able to obtain a hysteretic behavior, the parameters of the transformed function (i.e., ξ, σ and α) while heating (FM to PM transformation) or cooling (PM to FM transformation) should be different. Therefore, two transformed functions, f H and f C , are defined for the heating and cooling processes, respectively: Then by knowing the magnetization curves of both FM and PM phases, M FM and M PM , respectively, which are obtained from the Bean-Rodbell model (shown in Figure 1), the magnetization for the heating and cooling processes are expressed as: According to the Bean-Rodbell model, we assume a transformation process only in the metastable temperature region. Therefore, the transformed function for any set of parameters should fulfil that , the transformation cannot be extended further than the metastable region). An example of the modeled magnetization curves that can be obtained is included in Figure 1 using the normal distribution with T H t,MS = 307 K, T C t,MS = 305 K and σ H = σ C = 0.5 K. This kind of approximation in which the Bean-Rodbell model is combined with a transformation process has been applied successfully to the description of the hysteresis in MC materials [34].
Once the transformed fractions together with the magnetization values are obtained for both heating and cooling processes, the different reversal magnetization curves can be derived as a function of the temperature and the reversal temperature (according to the measurement protocol explained in the introduction) as: As an example, minor hysteresis loops to obtain the first order reversal magnetization curves while cooling and heating are shown in Figure 2. With these reversal curves, the TFORC distributions are calculated using Equations (1) and (2). while cooling and heating are shown in Figure 2. With these reversal curves, the TFORC distributions are calculated using Equations (1) and (2).

Symmetric Transformations
First, we tested the simplest cases for which both heating and cooling process are symmetric ( = = 0, the resulting distribution being the usual normal distribution) and both have the same transition temperature range ( = ), i.e., the FM to PM transformation is equivalent to the PM to FM one. The only difference is the location parameter of both processes ( and ). This symmetric transformation is typically ascribed to nucleation-growth phenomena or to transition temperature distributions.

Symmetric Transformations
First, we tested the simplest cases for which both heating and cooling process are symmetric (α H = α C = 0, the resulting distribution being the usual normal distribution) and both have the same transition temperature range (σ H = σ C ), i.e., the FM to PM transformation is equivalent to the PM to FM one. The only difference is the location parameter of both processes (ξ H and ξ C ). This symmetric transformation is typically ascribed to nucleation-growth phenomena or to transition temperature distributions. Figure 3a'-d' shows the obtained cooling TFORC diagrams (which are accompanied by their corresponding first order reversal magnetization curves Figure 3a-d) under different conditions (different transition temperatures and different ranges of the transformation). As a general feature, we observe a circular TFORC distribution around a maximum at certain coordinates (T u , T h ). In addition to that, and only for temperatures far from the maximum of the distribution, deviations from the circular shape can be observed (further discussed below). We can observe that the maximum distribution values are only dependent on the values of the transition temperatures while the radius distributions depend only on the σ values, increasing as the σ values increase. For the normal distribution, just the σ parameter is related to the transformation range (e.g., as the variance is only dependent on σ). The coordinates of the maxima coincide with the center of the loop (T center ) in the T h axis and with the half width of the hysteresis loop (∆T loop ) in the T h axis. On the one hand, T center can be experimentally calculated as the average value of the temperatures that make maximum ∂M/∂T for both heating and cooling processes. On the other hand, ∆T loop can be calculated as the half of the difference between the temperatures that make maximum ∂M/∂T for both heating and cooling processes. The maximum value of the distribution is highly dependent on the transformation range, increasing as the transformation range decreases.
Metals 2020, 10, x FOR PEER REVIEW 6 of 16 we observe a circular TFORC distribution around a maximum at certain coordinates ( , ). In addition to that, and only for temperatures far from the maximum of the distribution, deviations from the circular shape can be observed (further discussed below). We can observe that the maximum distribution values are only dependent on the values of the transition temperatures while the radius distributions depend only on the values, increasing as the values increase. For the normal distribution, just the parameter is related to the transformation range (e.g., as the variance is only dependent on ). The coordinates of the maxima coincide with the center of the loop ( ) in the axis and with the half width of the hysteresis loop (∆ ) in the axis. On the one hand, can be experimentally calculated as the average value of the temperatures that make maximum / for both heating and cooling processes. On the other hand, ∆ can be calculated as the half of the difference between the temperatures that make maximum / for both heating and cooling processes. The maximum value of the distribution is highly dependent on the transformation range, increasing as the transformation range decreases.  Cooling first order reversal magnetization curves and associated Temperature First Order Reversal Curve (TFORC) distributions with different transformation characteristics: However, as mentioned above, for low ρ C . values (far from the maximum), a tail of small negative values located along the 135 • line with respect the maximum can be observed for all the diagrams. When analyzing the heating TFORC in comparison to the cooling one (shown in Figure 4 a',b'), only positive values are obtained and now a small tail with positive values is located along the 225 • line with respect to the maximum (this difference will be justified below). Apart from these features at low ρ values, both cooling and heating TFORC distributions are equivalent (it has also been proven for all the other cases studied in this work, even if not explicitly mentioned again). Slightly higher values of ρ. are obtained for the cooling compared to the heating one, being ascribed to the larger magnetization jump at the transition for the cooling curve. The origin of the observed tails is related to the different temperature evolution of the FM and PM phases, and the transformation processes have no influence on it, as can be mathematically demonstrated for the cooling case (heating would be analogous). Starting from Equation (11), ρ C (T, T r ) can be analytically calculated from Equation (1) leading to: Metals 2020, 10, x FOR PEER REVIEW 7 of 16 However, as mentioned above, for low values (far from the maximum), a tail of small negative values located along the 135° line with respect the maximum can be observed for all the diagrams. When analyzing the heating TFORC in comparison to the cooling one (shown in Figure 4 a',b'), only positive values are obtained and now a small tail with positive values is located along the 225° line with respect to the maximum (this difference will be justified below). Apart from these features at low values, both cooling and heating TFORC distributions are equivalent (it has also been proven for all the other cases studied in this work, even if not explicitly mentioned again). Slightly higher values of are obtained for the cooling compared to the heating one, being ascribed to the larger magnetization jump at the transition for the cooling curve. The origin of the observed tails is related to the different temperature evolution of the FM and PM phases, and the transformation processes have no influence on it, as can be mathematically demonstrated for the cooling case (heating would be analogous). Starting from Equation (11), ( , ) can be analytically calculated from Equation (1) leading to: is always positive (its minimum value is 0) as ⁄ is negative and ∆ positive. On the other hand, the sign of the first term depends on the sign of ∆ ⁄ , as (1 − ) is always positive. Close to the heating branch of the hysteresis loop, this term is negative (as can be easily inferred from Figure  1) while ⁄ ≅ 0 and then, the distribution takes negative values. This fits with the observations, as negative values are obtained for > and < ∆ . The same analysis can be done for the heating distribution (from Equations (2) and (12)), the result being in agreement with the observation of Figure 4b'. As we conclude that the observed tails come from the temperature dependence on FM and PM magnetization, it is possible to check this feature if we just reduce As ∆M = M FM − M PM . As ∂ f H /∂T r is always positive, the sign of the distribution is given by the sign of the bracket. On the one hand, we can see that the second term inside the bracket (− ∂ f C ∂T ∆M) is always positive (its minimum value is 0) as ∂ f C /∂T is negative and ∆M positive. On the other hand, the sign of the first term depends on the sign of ∂∆M/∂T, as (1 − f C ) is always positive. Close to the heating branch of the hysteresis loop, this term is negative (as can be easily inferred from Figure 1) while ∂ f C /∂T 0 and then, the distribution takes negative values. This fits with the observations, as negative values are obtained for T u > T center and T h < ∆T loop . The same analysis can be done for the heating distribution (from Equations (2) and (12)), the result being in agreement with the observation of Figure 4b'. As we conclude that the observed tails come from the temperature dependence on FM and PM magnetization, it is possible to check this feature if we just reduce M FM and M PM to 1 and 0, respectively. Making this, the resulting TFORC is just dependent on the transformation characteristics. Figure 4c' shows the cooling TFORC distribution without magnetization contribution for the cooling process, and it is clearly observed that previous tails do not appear (in agreement with the previous arguments). In fact, despite those tails, both distributions are quite similar to the ones using magnetization data, illustrating that TFORC is a useful tool to extract information on the transformation processes. At this point, it has to be mentioned that the general features of the transformation (e.g., asymmetry or width) could also be extracted from magnetization curves. However, M(T) curves along the transformation are considerably affected by the temperature dependence of the magnetic moment of the pure phases (mainly the FM) jeopardizing the extraction of transformed functions from magnetization data. Moreover, TFORC would make it possible to discern the case of several transformation processes during the transformation (possibly emerging from coexisting phases), which would not be evident from the simple observation of the M(T) curves, analogously to the advantages of conventional FORC vs. M(H) curves. These various transformation processes have been evidenced by TFORC in [19].
Next, we checked the cases for which both heating and cooling processes remain symmetric (α H = α C = 0) but with different transition ranges (σ H σ C ). Figure 5a'-c' shows the cooling TFORC distribution for different σ (for simplification, the transition temperatures are ξ H = T t,ST + 2 K and ξ C = T t,ST − 2 K for all the cases). It can be observed that as the difference between σ H and σ C appears, the circular behavior (shown in Figure 3a') changes to a rotated ellipse for which its major symmetry axis is located along the +45 • line with respect the maximum if σ H > σ C . As the difference between σ H and σ C increases, the ellipse becomes narrower. When σ C > σ H , the major symmetry axis of the ellipse is rotated 90 • . This observation highlights that not only the shape of the distribution is important but the orientation also gives us useful information. It can be noted that, although the shape is not circular, the ellipse keeps its center at T u = T center and T H = ∆T loop .
Metals 2020, 10, x FOR PEER REVIEW 8 of 16 and to 1 and 0, respectively. Making this, the resulting TFORC is just dependent on the transformation characteristics. Figure 4c' shows the cooling TFORC distribution without magnetization contribution for the cooling process, and it is clearly observed that previous tails do not appear (in agreement with the previous arguments). In fact, despite those tails, both distributions are quite similar to the ones using magnetization data, illustrating that TFORC is a useful tool to extract information on the transformation processes. At this point, it has to be mentioned that the general features of the transformation (e.g., asymmetry or width) could also be extracted from magnetization curves. However, M(T) curves along the transformation are considerably affected by the temperature dependence of the magnetic moment of the pure phases (mainly the FM) jeopardizing the extraction of transformed functions from magnetization data. Moreover, TFORC would make it possible to discern the case of several transformation processes during the transformation (possibly emerging from coexisting phases), which would not be evident from the simple observation of the M(T) curves, analogously to the advantages of conventional FORC vs ( ) curves. These various transformation processes have been evidenced by TFORC in [19].
Next, we checked the cases for which both heating and cooling processes remain symmetric ( = = 0) but with different transition ranges ( ≠ ). Figure 5a'-c' shows the cooling TFORC distribution for different (for simplification, the transition temperatures are = , + 2 K and = , − 2 K for all the cases). It can be observed that as the difference between and appears, the circular behavior (shown in Figure 3a') changes to a rotated ellipse for which its major symmetry axis is located along the +45° line with respect the maximum if > . As the difference between and increases, the ellipse becomes narrower. When > , the major symmetry axis of the ellipse is rotated 90°. This observation highlights that not only the shape of the distribution is important but the orientation also gives us useful information. It can be noted that, although the shape is not circular, the ellipse keeps its center at = and = ∆ .

Asymmetric Transformations
We can introduce asymmetry in the transformation processes by using non-null α H and α C factors in the model. This asymmetry of the transformations has been observed in experimental data and is typically linked to internal stress/constrains during the transformations [35,36]. At the first stages of the transformation, the emerging nuclei can grow without interaction (contact) with others, but growth is hindered at final stages (reducing the rate) by the effect of the nearby transformed fractions (potentially leading to the existence of kinetic arrest). With the proposed model, we are able to reproduce: (1) transitions for which the transformation rate is faster at initial stages and (2) transformations with slower rates at initial stages. These two features can be reproduced using positive and negative α parameters, respectively. To illustrate them, the temperature dependence of the magnetization at the transition and its derivative with respect the temperature using different α parameters are shown in Figure 6. When the skewness is not considered (as was the case up to this section), the curves are symmetric (i.e., the maximum transformation rate occurs at the middle of the transformation), as can be observed in the derivative. When skewness is introduced, the distribution becomes asymmetric, with the maximum being shifted either to low transformed fractions when α is positive (i.e., the process is faster at the beginning of the transformation, before f = 0.5) or to high transformed fractions when α is negative (i.e., the process is faster at the end of the transformation, after f = 0.5). In this simple way, the possible transformation cases can be qualitatively modeled without the necessity of introducing complicated kinetic theories. factors in the model. This asymmetry of the transformations has been observed in experimental data and is typically linked to internal stress/constrains during the transformations [35,36]. At the first stages of the transformation, the emerging nuclei can grow without interaction (contact) with others, but growth is hindered at final stages (reducing the rate) by the effect of the nearby transformed fractions (potentially leading to the existence of kinetic arrest). With the proposed model, we are able to reproduce: 1) transitions for which the transformation rate is faster at initial stages and 2) transformations with slower rates at initial stages. These two features can be reproduced using positive and negative parameters, respectively. To illustrate them, the temperature dependence of the magnetization at the transition and its derivative with respect the temperature using different parameters are shown in Figure 6. When the skewness is not considered (as was the case up to this section), the curves are symmetric (i.e., the maximum transformation rate occurs at the middle of the transformation), as can be observed in the derivative. When skewness is introduced, the distribution becomes asymmetric, with the maximum being shifted either to low transformed fractions when is positive (i.e., the process is faster at the beginning of the transformation, before = 0.5) or to high transformed fractions when is negative (i.e., the process is faster at the end of the transformation, after = 0.5). In this simple way, the possible transformation cases can be qualitatively modeled without the necessity of introducing complicated kinetic theories.   Figure 7 a'-d' shows the TFORC distributions when considering that one process is asymmetric while the other remains symmetric, using different combination of the skewness factors (keeping the other parameters fixed to ξ H = T t,ST + 2 K, ξ C = T t,ST − 2 K, σ H = σ C = 0.5 K). It can be observed that the TFORC distribution has a semicircular-like shape. The semicircular shape is more evident as the asymmetry factor increases (Figure 7a',b'). Moreover, the orientation of the semicircle depends on the different combinations of parameters. The orientation is denoted with respect to the symmetry axis (which is the reflection axis). When the asymmetric transformation is the cooling one, the semicircle is orientated along 135 • or 315 • for positive and negative α values, respectively (Figure 7c',d'). Similarly, when the asymmetric transformation is the heating one, the semicircle is orientated along 45 • or 225 • for positive and negative α values, respectively. It should be noted that when introducing non-negative α values, the temperature range is no longer dependent only on σ (although it is the main factor) and therefore, the heating and cooling processes have slightly different transformation ranges.
the TFORC distribution has a semicircular-like shape. The semicircular shape is more evident as the asymmetry factor increases (Figure 7a',b'). Moreover, the orientation of the semicircle depends on the different combinations of parameters. The orientation is denoted with respect to the symmetry axis (which is the reflection axis). When the asymmetric transformation is the cooling one, the semicircle is orientated along 135° or 315° for positive and negative values, respectively ( Figure  7c',d'). Similarly, when the asymmetric transformation is the heating one, the semicircle is orientated along 45° or 225° for positive and negative values, respectively. It should be noted that when introducing non-negative values, the temperature range is no longer dependent only on (although it is the main factor) and therefore, the heating and cooling processes have slightly different transformation ranges.
Different situations when both transformation processes are considered asymmetric (same α factors) are illustrated in Figure 8a'-d', for which ξ H = T t,ST + 2 K, ξ C = T t,ST − 2 K, σ H = σ C = 0.5 K have been kept constant. It is inferred from the plots that TFORC distributions have isosceles triangle-like shape. It should be noted that as σ and α heating and cooling factors have the same absolute value, the transformation range is the same. It is observed that one vertex of the triangle is always pointing to the central coordinate (i.e., T u = T center and T h = ∆T loop ) but its reflection symmetry axis is oriented with different angles, 0, 90, 180 and 270 • for α H = α C > 0, α H = −α C > 0, α H = −α C < 0 and α H = α C < 0, respectively. It is worth noting that the orientations observed for the triangular distributions are not the same of those semicircular ones. If both skewness factors are different, the isosceles triangle is modified to a scalene one as the triangle does not have any symmetry axis.
pointing to the central coordinate (i.e., = and = ∆ ) but its reflection symmetry axis is oriented with different angles, 0, 90, 180 and 270° for = > 0, = − > 0, = − < 0 and = < 0, respectively. It is worth noting that the orientations observed for the triangular distributions are not the same of those semicircular ones. If both skewness factors are different, the isosceles triangle is modified to a scalene one as the triangle does not have any symmetry axis. Another possible effect in addition to both asymmetric processes is that both transformations also present different parameters (as mentioned before, with being the main factor affecting the transformation range). We take as an example the previous case with = = 4 (shown in Another possible effect in addition to both asymmetric processes is that both transformations also present different σ parameters (as mentioned before, with σ being the main factor affecting the transformation range). We take as an example the previous case with α H = α C = 4 (shown in Figure 8a') with different σ values when cooling and heating, illustrated in Figure 9a'-c'. When σ H = σ C the distribution has a triangular-like shape, as previously observed. However, when σ H σ C the triangle shape is distorted. The distortion produced is quite similar to the one studied previously in Figure 5. It can be observed that as the difference between σ H and σ C increases, the triangular-like shape is distorted, becoming elongated along the 45 • line if σ H > σ C (Figure 9a'). If σ H < σ C , the elongation is produced along the −45 • line, this being more significant as the differences increase (Figure 9b',c'). The isosceles triangle is modified to a scalene one (similar to the case when α are different) if the differences between σ H and σ C are not significant, or to a skewed ellipse if the differences are quite relevant. the triangle shape is distorted. The distortion produced is quite similar to the one studied previously in Figure 5. It can be observed that as the difference between and increases, the triangular-like shape is distorted, becoming elongated along the 45° line if > (Figure 9a'). If < , the elongation is produced along the −45° line, this being more significant as the differences increase (Figure 9b',c'). The isosceles triangle is modified to a scalene one (similar to the case when are different) if the differences between and are not significant, or to a skewed ellipse if the differences are quite relevant. In order to summarize the main results of the work, the different characteristic TFORC distributions as a function of the transformation parameters ( and ) are tabulated in Table 1. It is worth mentioning that the published TFORC distributions of magnetocaloric materials [19] include some of the features predicted in this work, like triangular shaped maxima and multiple peaks corresponding to multiple transitions. The combination of the FORC distribution catalogue with the details of the transformation function will make it possible to interpret the features giving rise to the hysteresis in magnetocaloric materials and eventually lead to designing strategies to minimize the irreversibility of the transformation. At the same time, the knowledge of the minor loops allows the implementation of more accurate models for the simulation of the performance of magnetocaloric devices.
In order to summarize the main results of the work, the different characteristic TFORC distributions as a function of the transformation parameters (σ and α) are tabulated in Table 1. It is worth mentioning that the published TFORC distributions of magnetocaloric materials [19] include some of the features predicted in this work, like triangular shaped maxima and multiple peaks corresponding to multiple transitions. The combination of the FORC distribution catalogue with the details of the transformation function will make it possible to interpret the features giving rise to the hysteresis in magnetocaloric materials and eventually lead to designing strategies to minimize the irreversibility of the transformation. At the same time, the knowledge of the minor loops allows the implementation of more accurate models for the simulation of the performance of magnetocaloric devices.

Conclusions
Through the modeling of first-order thermomagnetic transitions using the Bean-Rodbell model extended with phase transformation processes by skew normal distributions, it has been possible to extract information of TFORC distributions for magnetocaloric materials. It is relevant to highlight that a TFORC distribution (measured either upon cooling or upon heating) contains information of both cooling and heating transformations. In the simplest case, for which both heating and cooling transitions are symmetric and equivalent, the TFORC distribution has a circular shape around a maximum value. When both transitions are symmetric but with different temperature ranges, the TFORC distribution has an ellipsoidal shape with the major symmetry axis along the 45 • (if the heating process is broader than the cooling one) or −45 • lines (if the cooling process is broader than the heating one). If one transformation is asymmetric, the TFORC distribution has a semicircular-like shape, the reflection axis oriented along 45, 135, 225 and 315 • if α H > 0 and α C = 0, α H = 0 and α C < 0, α H < 0 and α C = 0 and α H = 0 and α C > 0, respectively. When it is considered that both transformations are asymmetric (the maximum transformation rate occurs before or after the center of the transformation temperature span), the TFORC distribution has an isosceles triangular-like shape with its reflection axis along 0, 90, 180 and 270 • orientations for α H = α C > 0, α H = −α C > 0, α H = −α C < 0 and α H = α C < 0, respectively. Finally, it is observed that when the transformations are asymmetric and occur with different transformation ranges, the triangular-like shape becomes distorted along the 45 • (if the heating process is broader than the cooling one) or 315 • lines (if the cooling process is broader than the heating one), obtaining scalene triangle or skewed ellipse shapes. This information can be useful when analyzing TFORC distributions of first-order MC materials and sets the basis for correlating the observed TFORC distributions with characteristics of the transformation.