An Analytical Approach for Computing the Coefﬁcient of Refrigeration Performance in Giant Inverse Magnetocaloric Materials

: An analytical approach for computing the coefﬁcient of refrigeration performance (CRP) was described for materials that exhibited a giant inverse magnetocaloric effect (MCE), and their governing thermodynamics were reviewed. The approach deﬁnes the magnetic work input using thermodynamic relationships rather than isothermal magnetization data discretized from the literature. The CRP was computed for only cyclically reversible temperature and entropy changes in materials that exhibited thermal hysteresis by placing a limit on their operating temperature in a thermodynamic cycle. The analytical CRP serves to link meaningful material properties in ﬁrst-order MCE refrigerants to their potential work and efﬁciency and can be employed as a metric to compare the behaviors of dissimilar alloy compositions or for materials design. We found that an optimum in the CRP may exist that depends on the applied ﬁeld level and Clausius–Clapeyron (CC) slope. Moreover, through a large literature review of NiMn-based materials, we note that NiMn(In/Sn) alloys offer the most promising materials properties for applications within the bounds of the developed framework.


Introduction
Giant inverse magnetocaloric materials have gained recent attention in the scientific community due to their high magnetic-to-thermal energy conversion efficiencies and potential for eliminating environmentally harmful chemicals typically used in conventional vapor-compression refrigeration units. Residential and commercial heating and cooling appliances account for nearly 40% of the United States' energy consumption, totaling 20.26 quadrillion British thermal units (BTU) annually. As global populations continue to grow, the amount of energy allocated for heating and cooling applications is expected to increase by 84% by the year 2050, compared to energy levels from 2010 [1,2]. The environmentally damaging fluids used in vapor compression systems, such as chlorofluorocarbons (CFCs) or hydrofluorocarbons (HCFCs) [3,4], are well-known to be leading contributors to greenhouse gas emissions and can destroy stratospheric ozone [5,6].
Solid-state refrigeration devices, such as magnetocaloric and thermoelastic regenerators, have the potential to mitigate greenhouse gas emissions [7], as they do not require HCFCs in thermodynamic cooling cycles. In a report submitted to the U.S. Department of Energy [8], magnetocaloric and thermoelastic technologies are projected to reduce energy demands for HVAC systems in both commercial and residential sectors by at least 20%. However, to succeed in that goal, there must be further development of caloric materials. Thus, the authors aim to define an analytical approach for computing a performance metric useful for comparing the energy conversion efficiency in novel solid-state refrigerant materials to facilitate high-throughput performance screening.
In previous works, Wood and Potter developed a thermodynamic metric, namely the refrigeration capacity (RC), that could be applied to materials that exhibited secondorder magnetic transitions without hysteretic loss [9]. Nearly 20 years following their work, Pecharsky and Gsneidner conceived a similar metric, namely the relative cooling power, RCP, which was applied to the first-order magnetostructural transition observed in GdSiGe [10,11]. They also employed the RCP to compare materials that exhibited a second-order magnetic transition, such as pure Gd [10]. However, using the RCP to compare materials that exhibited hysteretic losses across magnetostructural transitions with those that did not (e.g., across purely magnetic second-order transitions) led to misleading comparisons. More recently, other metrics applicable to materials exhibiting first-order magnetostructural transitions have been derived by the present author [12], Moya [13,14], and Brück [15]. For instance, Brück et al. defined a coefficient of performance, COP, which accounted for hysteretic losses through its implementation and was applicable to the physics across first-order structural phase transformations. It was mentioned that, in quantifying the applied magnetic-work required to generate the cooling in the magnetocaloric material (i.e., energy input), the isothermal magnetization loops needed to be numerically integrated. Data for complete magnetization loops across magnetostructural transitions were limited in the literature and were cumbersome to digitize.
Here, we aim to analytically compute the metric so that commonly measured and reported materials properties can be used to quickly compare the performance of novel magnetocaloric materials. In this way, solid refrigerants that exhibit a first-order magnetostructural transition with transformation hysteresis can be quickly assessed. First, the thermodynamics governing the conventional (second-order) caloric effects in solids are briefly described, which can be defined for any externally applied stimuli (thermodynamic force) that influences a material's free energy. The thermodynamics are then incorporated into Wood and Potter's RC parameter. Next, the thermodynamics governing first-order structural phase transformations are described. The entropy change and adiabatic temperature change are derived from the Clausius-Clapeyron relationship and a method developed by Porcari et al., [16,17], respectively. These are later applied to Wood and Potter's RC parameter, and the disadvantages of employing this metric to first-order structural transitions are revealed.
Finally, we present an analytical approach to compute the coefficient of refrigeration performance (CRP) that is directly applicable to first-order structural phase transitions. The analytical approach accounts for hysteretic losses and employs basic and commonly reported materials properties in the literature, rather than cumbersome digitized magnetization data. The development of the analytical CRP is followed with a discussion and examples used for comparing the caloric effects in first-order magnetostructural materials.

Thermodynamic Quantities of Interest
A caloric, or thermal, effect is defined by a material's isothermal entropy (∆S iso ) or adiabatic temperature (∆T ad ) change in response to some external stimulus (force). For the elastocaloric effect (ECE), this stimulus is mechanical load, and for the magnetocaloric effect (MCE), a magnetic field. Reversible heating and cooling effects can be generated in a solid material by a change in the material's free energy, which might include latent heat effects from a structural transition. Heating can also be generated through any irreversible processes that occur at the time of applying the external field. Magnetic and structural hystereses serve as indicators that thermodynamically irreversible processes take place, and through the second law of thermodynamics, the entropy produced (and thus heat) from these irreversible processes will always be greater than zero.
The internal energy of a substance changes with an applied thermodynamic "force", such as stress, temperature, and magnetic field. These forces do not depend on the material's volume and are referred to, herein, as intensive thermodynamic variables [18]. Since we aim to determine the caloric effects resulting from a change in internal energy, we consider the Gibbs free energy, G, which implicitly assumes intensive variables are applied to the substance in question during experiments [18]. Here, the G describes the free energy of a single structural phase. The change in free energy, dG, of that phase is defined as [18] where V is specific volume (m 3 ·kg −1 ), P is hydrostatic pressure (Pa), S is the entropy (J·kg −1 ·K −1 ), T is temperature (K), X i are extensive (volume dependent) material properties, such as bulk magnetization M (emu·g −1 ) or specific strain ε (m 3 ·kg −1 ), and y i are their thermodynamic intensive-force conjugates, such as applied magnetic field, µ 0 H (Tesla), or uniaxial stress, σ (Pa), respectively. Note that µ 0 is the permeability of free space in Henry·m −1 , and H is the applied magnetic field in A·m −1 . The isothermal entropy change defining the caloric effect can then be quantified using Equation (1). This is performed by employing Maxwell relations, first assuming that the application of the external stimuli is performed at atmospheric pressure and that all other applied forces are constant, i.e., dP = dy i = 0. The entropy, S, is expressed as The forces that are assumed constant in Equation (2) are denoted by subscripts. Similarly, other extensive thermodynamic quantities can be solved for from Equation (1) using the same approach. In the case, below, pressure and temperature are assumed constant.
Here, X i can be expressed as where P, T, and all other y i are constant intensive variables. Next, Equations (2) and (3) can be differentiated with respect to the other's independent variable. Equation (2) is differentiated with respect to y i resulting in and Equation (3) with respect to T, resulting in Assuming that the second partial derivative of G is smooth and continuous, the righthand side of Equations (4) and (5) are identical, and therefore, the left-hand sides are also mathematically equivalent. The incremental isothermal entropy change and "caloric effect", dS iso , from stimulus y i can then be denoted as which simplifies to In Equation (7), the isothermal entropy change has been derived for any single phase exposed to the intensive force y i in G. It is important to note that, in experiments where strain drives the caloric effect and is thus the independent variable, the Helmholtz free energy expression should be employed, instead of G, to insure consistency between the measured and predicted caloric behaviors [19].
If a magnetic field, y i = µ 0 H, is applied to a magnetocaloric refrigerant, the resulting entropy change is defined as (7), where X i has been replaced with M and is the energetic conjugate of µ 0 H. The MCE, or magnetic field-induced isothermal entropy change across a second-order magnetic transition, is illustrated in Figure 1a,b.
In Figure 1a, the entropy versus temperature diagram of a ferromagnetic material is depicted around its T Curie . Under zero magnetic field (top curve), the entropy is shown to exhibit a cusp at T Curie [32]. According to Equation (2), the entropy curve is defined by the partial derivative of the G with respect to temperature. Clearly, the second temperature derivative of the free energy, i.e., dS/dT, exhibits a discontinuity at T Curie , thus making the change in magnetization at T Curie a second-order magnetic transition. Interestingly, the assumption of commutativity, described above, no longer holds exactly at T Curie due to the inability to quantify the derivative of G. However, at temperatures above and below T Curie , Equation (7) is applicable for determining magnetic field-induced ∆S iso .
Magnetism 2022, 1, FOR PEER REVIEW 4 energy expression should be employed, instead of , to insure consistency between the measured and predicted caloric behaviors [19].
If a magnetic field, = 0 , is applied to a magnetocaloric refrigerant, the resulting entropy change is defined as ∫ ( ) ( 0 ) 0 0 by Equation (7), where has been replaced with and is the energetic conjugate of 0 . The MCE, or magnetic field-induced isothermal entropy change across a second-order magnetic transition, is illustrated in Figure 1a,b.
In Figure 1a, the entropy versus temperature diagram of a ferromagnetic material is depicted around its . Under zero magnetic field (top curve), the entropy is shown to exhibit a cusp at . [32] According to Equation (2), the entropy curve is defined by the partial derivative of the with respect to temperature. Clearly, the second temperature derivative of the free energy, i.e., / , exhibits a discontinuity at , thus making the change in magnetization at a second-order magnetic transition. Interestingly, the assumption of commutativity, described above, no longer holds exactly at due to the inability to quantify the derivative of . However, at temperatures above and below , Equation (7) is applicable for determining magnetic field-induced ∆ .   Upon applying µ 0 H to a ferromagnetic material around its T Curie , the total entropy (magnetic disorder) decreases, as depicted in Figure 1a,b. This can manifest an increase in temperature. Consequently, the MCE in most ferromagnetic materials that exhibit secondorder magnetic transitions is achieved through adiabatic or isentropic demagnetization, i.e., removing the applied magnetic field.
The adiabatic temperature change corresponding to a decrease in entropy is illustrated by a red horizontal arrow labeled as ∆T ad in Figure 1a. The magnitude of ∆T ad in a single structural phase material (described with the above Gibbs free energy) can be derived using the following thermodynamic assumptions and analyses. It is assumed that the total change in entropy, dS total , is a function of all other thermodynamic quantities as posited by the Truesdell's theory of equipresence and the structure of thermodynamic state functions [33]. As such, dS total is defined as According to Figure 1a, the ∆T ad is produced when dS total = 0. Assuming isobaric (dP = 0) and isentropic (dS total = 0) conditions in Equation (8) where by the entropy change generated from applying y i can be moved to the left side of the equality. Thus, where dS dT = C p T per the second law of thermodynamics, and C p is the isobaric heat capacity [18]. Therefore, Finally, the terms in Equation (12) can be separated and integrated, leading to where the adiabatic temperature change can be computed for any isochoric driving force, y i , if the isobaric heat capacity, C p , and extensive X i histories are known. Oftentimes, C p is assumed to be independent of the driving force, y i , and ∆T ad is approximated as [34] ∆T It is important to note that Equation (14) has been developed using G for a single structural phase. Furthermore, this is simply an approximation, assuming only small fluctuations occur in C p by applying the external stimuli. Moreover, T in Equation (14) is nominally assumed to equal the starting temperature at which y i is applied during the adiabatic temperature change.

Wood and Potters Refrigeration Capacity for Magnetic Cooling
The refrigeration capacity (RC) defines the suitability of a magnetic material for cooling applications and was originally conceived by Wood and Potter [9] approximately forty years following the seminal work of Weiss and Piccard [35]. Wood and Potter proposed that a measure of how well a magnetic refrigerant will operate should be based on its capability to perform thermal work through field cycling used in magnetic regenerative heat cycles. They considered that the entropy change at the cold end of a heat cycle should not exceed that at the hot end to satisfy the second law of thermodynamics. As such, the "thermal work" that was performed by a magnetocaloric material was defined as the entropy change achieved at the cold temperature reservoir of a heat cycle, ∆S c , multiplied by the temperature gradient across which the heat was moved, ∆T. Originally, the ∆T was defined as the difference between the hot (T h ) and cold (T c ) reservoir temperature of a regenerative cycle, i.e., ∆T = T h − T c . Due to the ambiguity of the proposed hot and cold temperature reservoirs, in practice, these temperatures were assumed to be the full-width-half-maximum (T FW HM ) of the entropy change versus temperature diagram, as shown in Figure 1b [32]. Figure 1b illustrates the RC computed for conventional second-order magnetocaloric refrigerants from the above entropy change versus temperature diagram. To generate the entropy change versus temperature curves, the difference was computed between the zero field (µ 0 H = 0) and applied field (µ 0 H > 0) curves from Figure 1a. As shown in Figure 1b, an entropy change at some cold temperature, T c , is defined as ∆S c , whereas an entropy change of ∆S h is achievable at temperature T h . Thus, the RC, was defined as [9] Clearly, the goal of much present-day research is increasing the RC in state-of-the-art MCE materials [7,36], because it corresponds to the amount of reversible thermal work performed by the refrigerant.
In magnetocaloric refrigerants exhibiting a second-order magnetic transition, the largest entropy change occurs around the T Curie , but this temperature does not change as a function of applied field. As such, applying greater magnetic fields to second-order magnetocaloric compounds simply increases the ∆S c and only marginally broadens the entropy change peak (i.e., T FW HM in Equation (15)). To compare materials using Wood and Potter's RC parameter, Equations (14) and (7) were substituted into Equation (15), so RC = ∆S c ∆T ad . Using ∆T ad in place of T FW HM leading to a thermodynamic performance metric for second-order phase transforming materials in heat cycles around their T Curie [10]. It is important to note, however, that ∆S c and ∆T ad are both dependent on the applied field level and operating temperature. In the following section, we describe the magnetocaloric effects that are generated across first-order magnetostructural phase transitions and how they differ from those in the above sections. Performance criteria are analytically developed for reversible behavior in first-order materials so they may be properly compared with reversible second-order systems.

Thermodynamics of Giant Caloric Effects across First-Order Meta-Magnetic Transitions
First-order magnetostructural transitions are observed in many multifunctional alloys [37]. For simplicity, we focus on the observed behaviors in a well-known NiMnIn meta-magnetic shape memory alloy (MMSMA) [38,39]. MMSMAs are rare-earth-free intermetallic compounds that exhibit a reversible first-order magnetostructural diffusionless phase transformation between low temperature martensite (M) and high temperature austenite (A) and are candidates for cooling applications [40][41][42]. In our analyses, isothermal magnetization data measured in Ni 48 Mn 38 In 14 (at.%) polycrystals were used to describe the meta-magnetic behavior across first-order magnetostructural transitions and the caloric effect. The NiMnIn polycrystals, herein, were subjected to a homogenization treatment at 1123 K for 24 h in a protective argon atmosphere and then water quenched.
At temperatures below a critical transformation temperature, T 0 , NiMnIn alloys are composed of a non-magnetic M-phase (paramagnetic, superparamagnetic, etc.) [43][44][45][46][47]. This is illustrated by the thermomagnetic curves in Figure 2a. On heating the M-phase in a small µ 0 H, the MMSMA begins to spontaneously magnetize (see green curve) as it transforms into the high temperature A-phase and is accompanied by an endothermic (cooling) reaction from the latent heat of the martensitic transition. On heating, critical M-to-A transformation temperatures are referred to as the austenite start (A s ) temperature and austenite finish (A f ) temperatures. Conversely, cooling causes an exothermic A-to-M transformation at the critical temperatures denoted as martensite start (M s ) and martensite finish (M f ). The difference in A f and M f has been used to describe energy barriers to the martensitic transformation, including elastic and irreversible parts [48][49][50][51][52], and is the main reason for developing an analytical performance metric for only the reversible behavior in first-order materials, as second-order materials do not exhibit these same energy losses or barriers. Moreover, second-order materials exhibit the conventional magnetocaloric effect as described by Equation (7), but first-order materials exhibit magneto-structural transformations which are accompanied by a latent heat component and therefore exhibit "giant" magnetocaloric effects.
Magnetism 2022, 1, FOR PEER REVIEW 7 At temperatures below a critical transformation temperature, 0 , NiMnIn alloys are composed of a non-magnetic M-phase (paramagnetic, superparamagnetic, etc.) [43][44][45][46][47]. This is illustrated by the thermomagnetic curves in Figure 2a. On heating the M-phase in a small 0 , the MMSMA begins to spontaneously magnetize (see green curve) as it transforms into the high temperature A-phase and is accompanied by an endothermic (cooling) reaction from the latent heat of the martensitic transition. On heating, critical M-to-A transformation temperatures are referred to as the austenite start ( ) temperature and austenite finish ( ) temperatures. Conversely, cooling causes an exothermic A-to-M transformation at the critical temperatures denoted as martensite start ( ) and martensite finish ( ). The difference in and has been used to describe energy barriers to the martensitic transformation, including elastic and irreversible parts [48][49][50][51][52], and is the main reason for developing an analytical performance metric for only the reversible behavior in first-order materials, as second-order materials do not exhibit these same energy losses or barriers. Moreover, second-order materials exhibit the conventional magnetocaloric effect as described by Equation (7), but first-order materials exhibit magneto-structural transformations which are accompanied by a latent heat component and therefore exhibit "giant" magnetocaloric effects. Upon applying a constant magnetic field, * , to the MMSMA and measuring magnetization across the reversible martensitic transition, the , , , and temperatures experience a decrease (see blue curve in Figure 2a), indicating that the magnetic field stabilized the high-temperature ferromagnetic A-phase. Indeed, each critical temperature exhibits a unique sensitivity to the applied field, also referred to as the Clausius-Clapeyron (CC) slopes, / ( 0 ), / ( 0 ), / ( 0 ), and / ( 0 ). Often, in the literature, these slopes are assumed to be equivalent to / ( 0 ); however, this can be an oversimplification. In inverse giant MCE materials, the CC slopes are negative in sign.
If the MMSMA is held at constant temperature under zero magnetic field, 0 =0 , then there exists a field level ( 0 = * ) at which the M-phase will completely transform to the A-phase, as shown in Figure 2a Upon applying a constant magnetic field, H * , to the MMSMA and measuring magnetization across the reversible martensitic transition, the M f , M s , A s , and A f temperatures experience a decrease (see blue curve in Figure 2a), indicating that the magnetic field stabilized the high-temperature ferromagnetic A-phase. Indeed, each critical temperature exhibits a unique sensitivity to the applied field, also referred to as the Clausius-Clapeyron (CC) slopes, dM f /d(µ 0 H), dM s /d(µ 0 H), dA s /d(µ 0 H), and dA f /d(µ 0 H). Often, in the literature, these slopes are assumed to be equivalent to dM s /d(µ 0 H); however, this can be an oversimplification. In inverse giant MCE materials, the CC slopes are negative in sign.
If the MMSMA is held at constant temperature M f under zero magnetic field, M . Therefore, we can define H * as the field level needed achieve a complete and reversible isother- [52].
In this expression, d(µ 0 H)/A f is assumed to be a positive quantity that simply denotes the sensitivity of the A f temperature to the applied µ 0 H. Similarly, the magnetic field required to initiate the meta-magnetic transition at M µ 0 H=0 f can be defined by In this context of isothermal meta-magnetic magnetic field-driven transformations, the volume fraction of stabilized A-phase, ξ, can be computed as Clearly, if the MMSMA is constrained to a constant temperature equal to M µ 0 H=0 f and then H * is applied, the complete M-to-A transition will occur. Moreover, subsequent removal of the field will produce the complete A-to-M transformation, making it cyclically reversible. If the same process was performed above M µ 0 H=0 f , only partial or no transformation will occur on subsequent field cycling.
The magnetization and demagnetization response across the complete magnetic field induced transformation at M µ 0 H=0 f in NiMnIn is depicted (blue) in Figure 2b up to H * . In application, the endothermic reaction resulting from the M-to-A transition during adiabatic loading contributes to the temperature change of the MMSMA. Therefore, if the MMSMA was originally exposed to a magnetic field large enough to initiate transformation (see H req in the figure) at M µ 0 H=0 f , the adiabatic temperature change must be overcome by applying an even greater magnetic field to progress the transformation to completion. For instance, if the complete M-to-A transformation produced ∆T ad of −5 K, then the temperature barrier that must be overcome to complete the transition could be described by A In terms of quantifying the magnetic field needed to complete the adiabatic transition, we assume the starting temperature is less than M To illustrate the adiabatic magnetization curve up to H * ad , we have depicted an isothermal magnetization curve of our NiMnIn alloy at a temperature lower than M µ 0 H=0 f (green) and then found the average magnetization response (see dash-dot-dot curve) between the two temperatures. The adiabatic magnetization curve is essential for later development of the analytical CRP, as it defines the input energy required by the MMSMA to achieve the giant MCE.
First-order phase transitions are associated with a change in crystal symmetry (structure), which implies the G is discontinuous at the thermodynamic equilibrium temperature, T 0 . The free energy above and below T 0 describe that of two separate structural phases with independent physical properties including C p , M, V, etc. At T 0 , thermodynamic equilibrium between the two phases leads to the well-known CC relations [18,52,53]. According to Equation (1), the G in either phase of the bulk material is influenced by external stimuli including hydrostatic pressure, temperature, and y i (stress or magnetic field).
In the case of M-to-A first-order structural phase transitions [49,[53][54][55] the G of each phase is equivalent at the point of the transition; thus, so is their differential, Since each phase is exposed to the same driving forces at the point of transition, Assuming the phase transition occurs in a constant pressure atmosphere, dP = 0, Equation (17) can be simplified to Or where ∆S M→A is the entropy change, or difference, between the two phases at T 0 , ∆X M→A i is the measured change in the i th extensive property across the first-order transition, and dy i /dT 0 is the inverse CC slope [53]. For partial transformations, i.e., ξ < 1, only a fraction of the ∆X M→A i and ∆S M→A are achieved. Equation (19) is applicable to magnetocaloric materials that exhibit first-order magnetostructural phase transitions, such as NiMnIn, where y i = µ 0 H and ∆X i = ∆M M→A ; thus, ∆S M→A = −∆M M→A ξ·d(µ 0 H)/dT 0 and ξ = 1 for a complete transformation. Equation (19) appears to be similar to Equation (7); however, these two expressions were derived under dissimilar thermodynamic conditions and therefore provide the ∆S iso for fundamentally different processes. There has been much debate in the literature whether Equation (7) and its discretized form [56] are applicable to quantify ∆S M→A across first-order transitions [57]. While we recognize the large body of contributing discussions, this topic is out of the scope of the present work.
Using the described physics, we aim to develop an analytical form for the CRP that accounts for the cyclic irreversibility in first-order materials. We thus turn our attention to the entropy (S) versus temperature (T) diagrams, as was done with second-order materials (in Figure 1). The isofield total S vs. T curves are depicted in Figure 3 for our NiMnIn alloy [52]. Since S is an extensive thermodynamic property, its behavior parallels the thermomagnetic response depicted in Figure 2a. On heating the M-phase under no applied magnetic field, the MMSMA undergoes a structural transition starting at the A µ 0 H=0 s temperature and finishing at A µ 0 H=0 f , as indicated by the abrupt increase in S. On subsequent cooling, a thermal hysteresis is clearly observed when the material transforms back to the M-phase, and the S decreases abruptly.
Similar to the thermomagnetic response in Figure 2a, the S vs. T response in Figure 3a depicts a decrease in transformation temperatures under H * . It can be seen clearly in Figure 3a that if the MMSMA was originally at the M µ 0 H=0 f temperature and H * was applied, then ∆S iso will equal approximately the entropy difference between A and M. Applying this field results in a decrease in critical transformation temperatures.
Using the applied-and zero-field entropy versus temperature (S − T) data in Figure 3, the entropy change, ∆S, has been plotted below by computing the difference between the curves. The entropy change versus temperature depicted by the ∆S − T curve in Figure 3 and the corresponding labels are representative of the physics in all inverse NiMn-based first-order magnetostructural transitions. The cyan curve describing the entropy change represents those which are most commonly reported in giant magnetocaloric studies of MMSMAs [12]. It is crucial to note that the features of the ∆S − T curve align with critical transformation temperatures from the S − T diagram above, so they have been labeled explicitly.
Typically, in studies on giant MCE materials, the RC is computed using the commonly reported cyan ∆S − T curve by multiplying ∆S iso by T FW HM , as was done in Equation (15). Note the T FW HM encompasses the thermal hysteresis of the S − T diagram above. If field cycling was performed at a temperature above M µ 0 H=0 f , the MMSMA would only exhibit a partial entropy change on the first field ramping and zero on subsequent, and therefore would be cyclically irreversible-a large difference in the behavior that was described for second-order materials. Similar to the thermomagnetic response in Figure 2a, the vs. response in Figure  3a depicts a decrease in transformation temperatures under * . It can be seen clearly in Figure 3a that if the MMSMA was originally at the 0 =0 temperature and * was applied, then ∆ will equal approximately the entropy difference between A and M. Applying this field results in a decrease in critical transformation temperatures.
Using the applied-and zero-field entropy versus temperature ( − ) data in Figure  3, the entropy change, ∆ , has been plotted below by computing the difference between the curves. The entropy change versus temperature depicted by the ∆ − curve in Figure 3 and the corresponding labels are representative of the physics in all inverse NiMnbased first-order magnetostructural transitions. The cyan curve describing the entropy change represents those which are most commonly reported in giant magnetocaloric studies of MMSMAs [12]. It is crucial to note that the features of the ∆ − curve align with critical transformation temperatures from the − diagram above, so they have been labeled explicitly.
Typically, in studies on giant MCE materials, the RC is computed using the commonly reported cyan ∆ − curve by multiplying ∆ by , as was done in Equation (15). Note the encompasses the thermal hysteresis of the − diagram above. If field cycling was performed at a temperature above 0 =0 , the MMSMA would only exhibit a partial entropy change on the first field ramping and zero on subsequent,  As mentioned in Section 2.2, RC = ∆S c ∆T, where ∆T was assumed to equal T FW HM , and more recently, ∆T = ∆T ad [10,15]. In first-order materials, we propose ∆S c = ∆S iso , which is nearly equivalent to the entropy change across the transition (see Equation (19)). It is interesting to note that T FW HM appears to be related to critical transformation temperatures in Figure 3.
, the RC parameter will yield erroneously high values. Since it is assumed that T h = M µ 0 H=0 f , ∆S iso will be reversible with field cycling. Moreover, ∆T ad is expected to be cyclically reversible under these constraints.
Oftentimes, in the literature, ∆T ad is approximated using Equations (8)- (14) [58] for first-order phase transitions. However, through the author's experience, this method of computing ∆T ad often leads to approximations that nearly double that of measured values. In published works, similar discrepancies are typically attributed to non-ideal experimental conditions; however, the author posits that they are the result of neglecting the original assumptions imposed in the development of Equation (8). That is, the expression for dS total is only valid for one of the two structural phases that exist across first-order phase transitions, and Equation (8) is not applicable when two or more phases co-exist. In the simplification of Equation (14), we assumed dS/dT = C p /T from the second law of thermodynamics. However, A and M phases in first-order materials are well-known to exhibit differences in C p as large as 50 J·kg -1 ·K −1 [59]. Conversely, accurate predictions for ∆T ad have been developed by Cugini and Porcari et al. [16,17], whereby the S − T diagrams (see Figure 3) have been analyzed across first-order phase transitions with a discontinuous change in G.
In their work, ∆T ad was derived across a first-order phase transition in NiMn-based materials through graphical analysis of S − T diagrams, such as those shown in Figure 3. The ∆T ad was determined by assuming isentropic conditions (dS total = 0, as was performed in deriving Equation (9) where C M p is the isobaric heat capacity of the M phase depicted in Figure 3, and ∆S is the entropy change produced by applying y i , i.e., ∆S (M→A) . It is important to note that, since Equation (20) is empirical, dT 0 /dy i was assumed to be a positive quantity used in defining the length of one side of the similar triangles. Similar to H * , the sign of the CC slope is neglected here and simply represents the sensitivity of the transformation temperatures to an applied field. When comparing ∆T ad computed using Equations (14) and (20) with measured values across first-order transitions in the literature, it is clear that predictions computed from Equation (20) more closely match experiments [52]. Thus, for developing the analytical CRP, herein, Equation (20) was used to compute ∆T ad for materials demonstrating first-order magnetostructural transitions.

Refrigeration Capacity (RC) in First-Order Phase Transforming Materials
In a process similar to that applied to second-order materials, Equations (19) and (20) are substituted into Equation (15) to compute the RC in first-order systems. Assuming that the A-phase in the first-order material is ferromagnetic and the M-phase is non-magnetic, such as those depicted in Figure 2, applying a sufficient magnetic field to the M phase is expected to produce the M-to-A transition. On removing the applied field, the material will return to the M-phase if the temperature was at or below M µ 0 H=0 f . As such, cycling a sufficient magnitude of magnetic field at temperatures at or below M µ 0 H=0 f will produce reversible and repeatable entropy or temperature changes in the magnetocaloric solid. Such reversibility is analogous to that of second-order systems. Therefore, from Equations (19) and (20), assuming the only applied force is µ 0 H under isothermal conditions and that the magnetocaloric material was originally in the M-phase at the M µ 0 H=0 f temperature, the RC can be defined as which simplifies to where dT 0 /d(µ 0 H) = dA f /d(µ 0 H), and µ 0 H is, in fact, the isothermally applied magnetic field consistent with Equations (19) and (20) and bound between H req and H * . Within these constraints, if µ 0 H < H req , then RC = 0, and if µ 0 H > H * , the RC will saturate at a maximum. Per Equation (22), the RC in first-order materials can be approximated using the magnetization change across the magnetostructural transformation, ∆M M→A , the A f and A s temperature field sensitivity, the field-free transformation temperatures, and the heat capacity of the M-phase, C M p , at M µ 0 H=0 f .

Coefficient of Performance (CRP) in First-Order Phase Transforming Materials
Recent works [15,34,37,[60][61][62][63][64][65] have described the problems in using the RC as a metric to compare the MCE performance of first-order materials because it does not account for the energy input needed to drive the magnetostructural transformation. Moreover, the thermal gradient across which heat is moved is assumed to be the T FW HM (the bigger cyan curve in Figure 3b). A metric that addresses these shortcomings is the coefficient of (refrigeration) performance, CRP, defined as η = Q is the reversible thermal work that can be achieved using the solid refrigerant, and . W is the work (in this case magnetic work) required to excite the system. The coefficient of performance applied to a first-order materials was described as [15,37,62] where the numerator was the RC in Equation (15), and the denominator was the magnetic work applied to the first-order material at T Curie , where B = µ 0 H. The magnetic field applied to the first-order material in the present case is at M µ 0 H=0 f . As shown for completely reversible first-order transitions with field cycling, T Curie in Equation (23) must equal M µ 0 H=0 f . In the literature, the CRP was tabulated with Equation (23) for various materials using different quantities for ∆T for the RC in the numerator [37]. ∆T has been assumed to be T FW HM of the ∆S − T curve, ∆T ad , or ∆T cyclic , where ∆T cyclic was described to be the cyclically reversible ∆T ad at T Curie . Since no standard is used in the literature, CRPs greater than unity were found when assuming ∆T = T FW HM , because the effect of hysteresis was not removed from the calculation [37]. In these past works, it is still unclear whether (1) ∆T ad was computed using Equations (14) or (20) or was from direct measurement, and (2) if T Curie in Equation (23) corresponded to a thermodynamic equilibrium temperature, T 0 . In the NiMn-based material mentioned above, the ferromagnetic A-phase exhibited its own ferromagnetic T Curie , which is often above A

Analytical Approach for Computing the CRP in First-Order Phase Transforming Materials
In the present work, RC = ∆S M→A |∆T ad | as imposed by the bounds defining the RC-, i.e., the first-order material is only cyclically reversible, and thus comparable to secondorder materials when excited at or below M µ 0 H=0 f . As such, the numerator in Equation (23) has been defined by Equation (22). Furthermore, we have defined T h = M  Figure 3.
In application, the ∆T ad is achieved by ramping the magnetic field on the first-order material in thermally insulated conditions, and therefore, the magnetic field required to achieve the full temperature and entropy change at T h = M µ 0 H=0 f was defined by H * ad (see Section 3) and is depicted in Figure 2b. The magnetization change across the complete martensitic transition was defined as ∆M M→A , and therefore, the denominator of Equation (23) has been approximated as the area of a triangle (see magnetization response in Figure 2b), or Substituting Equations (24) and (22) in to Equation (23), and assuming H/dA f , per the above descriptions, results in the analytical form of the CRP for first-order materials. The analytical CRP reduces to (25) and is confirmed to be a dimensionless parameter, with limits from zero to unity, that can be quantified with basic and commonly reported materials properties in giant inverse MCE materials. The dA f /d(µ 0 H) in Equation (25) is assumed to be a positive quantity as a result of employing Equation (20) to empirically quantify ∆T ad . In a previous work [15], the community was encouraged to publish complete datasets, so that the CRP could be computed with numerically integrated digitized magnetization loops; however, here we find that the CRP can be approximated with , and C M p from only a few basic thermal and thermomagnetic experiments. Moreover, binding T h to M µ 0 H=0 f ensures that the CRP will be computed for only the reversible part of first-order transformations, which then lends the ability to employ the metric to compare the performance of first-order materials that undergo magnetostructural transformations and hysteresis with those that do not, i.e., purely ferromagnetic coolants.

Comparison of CRP in NiMn-Based Meta-Magnetic SMA
The analytical form of the CRP was given in Equation (25), which was defined by physically meaningful and fundamental materials properties. The metric can be used to compare the performance of magnetic refrigerants that exhibit inverse first-order magnetostructural transformations. Moreover, it allows for a uniform comparison in performance in both second-order and first-order materials. Thus, it can serve as a materials design tool to determine the optimal candidate for a solid-state cooling application with a specific T h and applied field level. In the present case, T h is bound to M ) is too large in comparison with dA s /d(µ 0 H) for any meaningful phase change to take place [66].
To demonstrate the use of the analytic CRP, data for inverse MCE materials were tabulated in Table 1 for over forty MMSMAs with various compositions. We found that C M p and dA s /d(µ 0 H) were scarce in the literature; therefore, Thermocalc was employed to approximate C M p at M µ 0 H=0 f for a few of the compositions and the dA s /d(µ 0 H), which was needed to compute H req , were obtained through digitization of commonly reported isofield thermomagnetization curves. The CRP was computed for each alloy, assuming 5 T was applied at the M µ 0 H=0 f temperature.
Ref. As shown in Table 1, nearly half of the listed MMSMA exhibited a non-zero CRP. The alloys that exhibited a CRP of zero required H req > 5 T. Since no M-to-A transformation would be generated by applying 5 T, for these alloys, the RC, i.e., the numerator of the CRP, was also zero. A larger applied field would result in more alloys exhibiting a nonzero CRP, because a M-to-A transformation may begin to take place. By subjecting all the MMSMAs listed in Table 1 to the same applied field, the alloys most suitable for a given MCE application could be identified as those with the largest CRP. Interestingly, the alloys containing indium offered the most prevalent non-zero CRP for all the NiMn-based MMSMAs with an applied field of 5 T.
Upon inspecting values in Table 1, the CRP was non-intuitive and could not be simply predicted using a single material property, such as dA f /d(µ 0 H). A parametric analysis of each material property would reveal its effect on CRP and would provide a clearer understanding of the important factors dictating performance. Prior to performing a parametric analysis, the CRP was computed for an architype MMSMA to use as a reference. Ni 50 Mn 35 In 15 at.% was selected as the archetypal alloy, which exhibited the properties listed in Table 2. Using the characterization parameters in Table 2 Table 2. Upper and lower bounds in ∆ → and from Table 1 are also listed in Table 2. Figure 4b depicts the CRP vs. 0 for various . Since the is present in only one term of the ∆ in Equation (20), and the ratio of the highest to lowest values as given in Table 2 is low (a factor of two times, compared to nearly 17 times for ∆ → ), the effect to the CRP is minor. Nonetheless, lower values would obviously produce greater ∆ and a higher RC. Higher ∆ would, in turn, increase * , which might decrease the CRP, because greater magnetic work would need to be applied to the material to generate the cooling effect. Thus, the CRP of an alloy may be increased by decreasing , although the change is relatively small.
Finally, Figure 5a-c depict the CRP vs. inverse CC slopes for an alloy with the same properties as the architype NiMnIn. The inverse CC slopes appear to have the greatest influence on the CRP out of all the material properties in Equation (25). This is not surprising, because the entropy change, or ∆ → in Equation (19), depends on the inverse CC slope. In Figure 5a, the CRP is depicted as a function of ( 0 )/ for 0 = 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, and 10.5 Tesla. For the Ni50Mn35In15 alloy, * = 10.5 T. As mentioned earlier, / ( 0 ) was used to compute * , * , and the RC in Equation (22); therefore, the influence of ( 0 )/ on the CRP is highly nonlinear. Interestingly, the results in Figure 5a suggest an optimum may exist in the CRP, which appears to be dependent on 0 . For example, if the cycled magnetic field was 0 = 2 T, a CRP of 0.  Table 1 that exhibit a CRP of zero also exhibit a / ( 0 ) with a magnitude less than −0.85 K·T −1 .   Table 2. Upper and lower bounds in ∆M M→A and C M p from Table 1 are also listed in Table 2. Figure 4b depicts the CRP vs. µ 0 H for various C M p . Since the C M p is present in only one term of the ∆T ad in Equation (20), and the ratio of the highest to lowest values as given in Table 2 is low (a factor of two times, compared to nearly 17 times for ∆M M→A ), the effect to the CRP is minor. Nonetheless, lower values would obviously produce greater ∆T ad and a higher RC. Higher ∆T ad would, in turn, increase H * ad , which might decrease the CRP, because greater magnetic work would need to be applied to the material to generate the cooling effect. Thus, the CRP of an alloy may be increased by decreasing C M p , although the change is relatively small.
Finally, Figure 5a-c depict the CRP vs. inverse CC slopes for an alloy with the same properties as the architype NiMnIn. The inverse CC slopes appear to have the greatest influence on the CRP out of all the material properties in Equation (25). This is not surprising, because the entropy change, or ∆S M→A in Equation (19), depends on the inverse CC slope. In Figure 5a, the CRP is depicted as a function of d(µ 0 H)/dA f for µ 0 H = 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, and 10.5 Tesla. For the Ni 50 Mn 35 In 15 alloy, H * = 10.5 T. As mentioned earlier, dA f /d(µ 0 H) was used to compute H * ad , H * , and the RC in Equation (22); therefore, the influence of d(µ 0 H)/dA f on the CRP is highly nonlinear. Interestingly, the results in Figure 5a suggest an optimum may exist in the CRP, which appears to be dependent on µ 0 H.  Table 1 that exhibit a CRP of zero also exhibit a dA f /d(µ 0 H) with a magnitude less than −0.85 K·T −1 .  Table 2.
Similarly, the CRP vs. ( 0 )/ was plotted in Figure 5b for an alloy with materials properties matching those of Ni50Mn35In15 at.%. The / ( 0 ) parameter was only used in quantifying , therefore affecting for a given 0 . The / ( 0 ) parameter does not appear to affect the maximum attainable CRP; however, the results in Figure 5b can be used to determine acceptable 0 needed achieve a non-zero CRP.
In most literature, only one CC slope is reported, i.e., the thermodynamic equilibrium point of the forward (A-to-M) martensitic transformation, but in most MMSMAs, there are four separate CC slopes. In the case of the inverse MCE, we focus on the reverse (Mto-A) martensitic transformation, and therefore, the / ( 0 ) and / ( 0 ) slopes are the most important. Assuming the thermodynamic equilibrium from M-to-A at the 0 =0 and 0 =0 exhibit the same sensitivity to 0 , the CRP vs. the inverse CC slope would exhibit features from both Figure 5a,b, as depicted in Figure 5c. In Figure 5c, non-zero CRP is attained at specific inverse CC slopes and 0 that also correspond to some CRP maximum.

Conclusions and Future Work
In the present work, an analytical approach for computing the coefficient of refrigeration performance (CRP) was described in the context of the giant inverse magnetocaloric effect (MCE), whereby discretized isothermal magnetization data were not required to compute the input work. The framework developed for computing the CRP from firstorder magnetostructural and meta-magnetic transformations was intended to mirror the cyclic reversibility attainable in second-order and strictly magnetic transitions. While computing the refrigeration capacity (RC), we found that cyclic reversibility was achieved if the "hot temperature reservoir", namely ℎ , was bound to the zero-field martensite finish temperature, 0 =0 . Due to this limitation, not all first-order meta-magnetic transformations were capable of exhibiting non-zero RC or CRP. A parametric analysis of the analytical CRP was performed, revealing that optimum Clausius-Clapeyron (CC) slopes may exist for a given magnitude of an applied magnetic field. Moreover, we found that the CRP can be maximized by ensuring a magnetization change across the martensitic transition, ∆ → , is greater than 27 emu·g −1 , which suggests non-magnetic martensite-to-paramagnetic austenite transformations are not useful for giant MCE applications.
Finally, due to the highly nonlinear behavior of the CRP with respect to the CC slopes, NiMnIn and NiMnSn alloys were found to exhibit the most favorable properties for giant inverse MCE applications when compared to over forty meta-magnetic shape memory alloys. This was attributed to the favorable CC slopes, ∆ → , and . In future work, we aim to employ the given framework to compare the CRP in first-order MCE  Table 2.
Similarly, the CRP vs. d(µ 0 H)/dA s was plotted in Figure 5b for an alloy with materials properties matching those of Ni 50 Mn 35 In 15 at.%. The dA s /d(µ 0 H) parameter was only used in quantifying H req , therefore affecting ξ for a given µ 0 H. The dA s /d(µ 0 H) parameter does not appear to affect the maximum attainable CRP; however, the results in Figure 5b can be used to determine acceptable µ 0 H needed achieve a non-zero CRP.
In most literature, only one CC slope is reported, i.e., the thermodynamic equilibrium point of the forward (A-to-M) martensitic transformation, but in most MMSMAs, there are four separate CC slopes. In the case of the inverse MCE, we focus on the reverse (M-to-A) martensitic transformation, and therefore, the dA s /d(µ 0 H) and dA f /d(µ 0 H) slopes are the most important. Assuming the thermodynamic equilibrium from M-to-A at the A µ 0 H=0 f and A µ 0 H=0 s exhibit the same sensitivity to µ 0 H, the CRP vs. the inverse CC slope would exhibit features from both Figure 5a,b, as depicted in Figure 5c. In Figure 5c, non-zero CRP is attained at specific inverse CC slopes and µ 0 H that also correspond to some CRP maximum.

Conclusions and Future Work
In the present work, an analytical approach for computing the coefficient of refrigeration performance (CRP) was described in the context of the giant inverse magnetocaloric effect (MCE), whereby discretized isothermal magnetization data were not required to compute the input work. The framework developed for computing the CRP from first-order magnetostructural and meta-magnetic transformations was intended to mirror the cyclic reversibility attainable in second-order and strictly magnetic transitions. While computing the refrigeration capacity (RC), we found that cyclic reversibility was achieved if the "hot temperature reservoir", namely T h , was bound to the zero-field martensite finish temperature, M µ 0 H=0 f . Due to this limitation, not all first-order meta-magnetic transformations were capable of exhibiting non-zero RC or CRP.
A parametric analysis of the analytical CRP was performed, revealing that optimum Clausius-Clapeyron (CC) slopes may exist for a given magnitude of an applied magnetic field. Moreover, we found that the CRP can be maximized by ensuring a magnetization change across the martensitic transition, ∆M M→A , is greater than 27 emu·g −1 , which suggests non-magnetic martensite-to-paramagnetic austenite transformations are not useful for giant MCE applications.
Finally, due to the highly nonlinear behavior of the CRP with respect to the CC slopes, NiMnIn and NiMnSn alloys were found to exhibit the most favorable properties for giant inverse MCE applications when compared to over forty meta-magnetic shape memory alloys. This was attributed to the favorable CC slopes, ∆M M→A , and C M p . In future work, we aim to employ the given framework to compare the CRP in first-order MCE materials with that in second-order materials at their respective T h (or M µ 0 H=0 f ) to identify the most favorable compositions for a specific application and operating temperature interval.

Conflicts of Interest:
The authors declare no conflict of interest.