On the Strong Composition Dependence of the Martensitic Transformation Temperature and Heat in Shape Memory Alloys

General derivation of the well-known Ren–Otsuka relationship, 1αdTodx=−αβ (where To, x, α and β(>0) are the transformation temperature and composition, as well as the composition and temperature coefficient of the critical shear constant, c′, respectively) for shape memory alloys, SMAs, is provided based on the similarity of interatomic potentials in the framework of dimensional analysis. A new dimensionless variable, tox=ToxTmx, describing the phonon softening (where Tm is the melting point) is introduced. The dimensionless values of the heat of transformation, ΔH, and entropy, ΔS, as well as the elastic constants c′, c44, and A=c44c′ are universal functions of to(x) and have the same constant values at to(0) within sub-classes of host SMAs having the same type of crystal symmetry change during martensitic transformation. The ratio of dtodx and α has the same constant value for all members of a given sub-class, and relative increase in c′ with increasing composition should be compensated by the same decrease in to. In the generalized Ren–Otsuka relationship, the anisotropy factor A appears instead of c′, and α as well as β are the differences between the corresponding coefficients for the c44 and c′ elastic constants. The obtained linear relationship between h and to rationalizes the observed empirical linear relationships between the heat of transformation measured by differential scanning calorimetry (DSC) (QA⟶M) and the martensite start temperature, Ms.


Introduction
It is well known that the martensitic transformation temperature has strong composition dependence in shape memory alloys (SMAs) [1][2][3][4][5] and, for instance, 1 at% composition change can alter the transition temperature by about 100 K.In addition, the change of the content of defects (e.g., the change of concentration of vacancies by quench) can also have a similar effect [1].The explanation of the above effect, even after the publication of a semiquantitative derivation in [1], is still under discussion in the literature (see, e.g., [4,[6][7][8][9][10][11][12][13]).In addition, it was found (see, e.g., Figure 5 in [3] and Figures 5b and 5d in [4]) that the transformation heat was a common linear function of the martensite start temperature, M s , in some different SMAs.We will also discuss this relationship and will derive it.
It was shown in [1], from a Landau-type model of first-order phase transformations, that the martensitic transformation (MT) occurs at a critical elastic (basal plane shear) constant, c ′ (=(c 11 − c 12 )/2) and the value of it is constant at the transformation temperature, T o .Considering that the elastic constant has strong composition, x, and temperature dependence, it was concluded in [1] that the constancy of c ′ at T o demands that the transformation temperature must exhibit an opposite effect, i.e., if c ′ increases, T o should decrease with increasing composition.Thus, they arrived at the requirement Assuming that the composition and temperature dependence of c ′ can be expressed as where c ′ o is the value of c ′ at T o , the following relationship was obtained: or, using the definition of α, It was also shown in [3] that, also taking into account that, typically, β = 4% 100K = 4 × 10 −4 1 K [1-3,5,13] and α = 4 − 10 % at% = (4 − 10) [1,11], dT o dx has a negative sign and can be given as This means that the transformation temperature is strongly affected by even a small change in composition or by quench.Investigating the validity of predictions ( 3) and ( 4) in [1], experimental data obtained in Cu-based shape memory alloys [2,3] as well as in Ti 50 Ni 30 Cu 20 alloys [11] were used, implicitly suggesting that (3) and ( 4) can be general predictions for all SMAs, although according to their derivation c ′ and dT o dx can have different constant values at T o (0) in different alloy systems.Experimental data supported this expectation: the values of c ′ and the dT o dx slopes were slightly but definitely different in different SMAs [2,3,7,[14][15][16][17].For instance, it was shown in [14] that the value of c ′ at T o was 30% smaller in Ti 50 Ni 30 Cu 20 alloy than in binary NiTi.Furthermore, it was demonstrated that in Ni 2 MnGa dT o dx was even positive (and α was negative) [7,15] and, e.g., the value of dT o dx in TiPd-based shape memory alloys [17] with different alloying elements varied by about a factor of four (i.e., it changed between −15 K/at% and −60K/at%) by changing the type of the third alloying element.Since the softening of the corresponding elastic moduli is a key characteristic for martensitic transformations in SMAs, and not only c ′ , but c 44 (belonging to non-basal plane shear) can also show softening, the above Ren-Otsuka relationships are expected to be valid only if c ′ has phonon-softening behavior and c 44 is practically independent of the temperature [11][12][13]18].This latter assumption is, in a good approximation, valid for Cu-based alloys [5,18] or for Ti 50 Ni 30 Cu 20 [14] but, e.g., in binary NiTi, both above moduli have phonon-softening-related temperature and composition dependence [11][12][13]18] and it was concluded (see, e.g., [14]) that the transformation temperature is more sensitive to the variation of c 44 than to that of c ′ .
From the state-of-the-art, summarized above, it is clear that there is a need to generalize the semi-quantitative derivation of (3) and to get a definite answer on the following question: whether the predicted constants at T o are different in different alloys or there exists a classification rule providing groups of host SMAs, within which they are the same.
In this paper, we provide a different, general derivation of relationships of type (3) and (4), based on the law of corresponding states (LCSs) for metals with phonon softening.It will be shown that the general forms of (3) or (4), which contain the composition and temperature dependence of both c ′ and c 44 , are just the consequence of the similarity of interatomic potentials [19].For the derivation of it, one can avoid the use of phrasing like "critical value at the transition temperature" [1,4,5], or "criticality of the austenite", which can be typical formulation for second-order phase transformations [20].This is in line with the Ren-Otsuka approach, in which also no complete mode softening (if e.g., c ′ −→ 0) is required [18] for first order phase transformations.In addition, since the explanation of the strong composition dependence of the heat of transformation, ∆H, and the linear relationship between ∆H and M s are still the question under debate [3,4,[6][7][8]21] (e.g., in [21], it was concluded that composition dependence of ∆H "remains to be rationalized"), these will also be discussed.
The organization of the paper is as follows.Basic relationships for the dependence of the transformation heat (∆H), entropy (∆S), the shear constants c ′ , c 44 , and the anisotropy constant, A = c 44 c ′ on t o = T o T m (T m is the meting point) are given in Section 2. The validity of the derived linearized relationships between h = ∆H kT m and t o is demonstrated on the examples of binary NiTi alloy (where the concentration, x, denotes the deviation from the stoichiometric (50/50 at%) composition on the Ni-rich side), in Cu 74.08 Al 13.13 Be 2.79 alloys (where x shows the increase in the Be content from x Be = 2.79 at%), in Ti 50−x Ni 40+x Cu 10 (where x changes between 0 at% and 1.2 at%) as well as in Ni 2 MnGa alloys (where x is the Ni excess in at% from the stoichiometric composition).In Section 3, our predictions will be compared with other experimental data.Section 4 contains the conclusions.

Law of Corresponding States for Phonon-Softening Systems
The LCS is the consequence of similarity of interatomic potentials [19], which can be written, in general, as: where r i is the position vector of the particle i (N is the number of atoms).Its form along a given direction can be represented by a periodic function with a period of a o and with an energy parameter of ε = f max − f min (e.g., similar to a sinus-type function with wavelength λ = a o and amplitude A = ε 2 , where a o is the nearest neighbor distance).This shape is similar for all solids of the same bonding type (e.g., for metals and metallic alloys) and crystal structure, forming so-called similarity classes, i.e., f is the same function of its arguments within a similarity class.
In order to derive useful relationships between different physical quantities, one has to start from the fundamental theorem of dimensional analysis ("Pi theorem") [22].This is based on the dimensional homogeneity.It is well known that, in physics, there are fundamental and derived quantities.The fundamental quantities are dimensionally independent, and their number, n, is finite.It is easy to show that, e.g., the mass, m, the energy, ε, the length, a o , and the Boltzmann constant, k, can be taken as dimensionally independent quantities (the physical dimension of none of them can be combined from the others), forming the basis of the dimensional analysis [19,22].According to the fundamental theorem, if Q denotes a (derived) physical quantity, then it can be given in the following form [19]: Here, q and its variables q 1 ,. . ., q g−n are dimensionless.This also means that q = Q/m a ε b a o c k d is dimensionless, i.e., the exponents a, b, c, and d should be chosen in such a way that the above dimensional combination of the fundamental quantities should give the dimension of Q.Furthermore, the number of independent dimensionless variables q i , in principle, is equal to g − n, where g is the number of variables present in the physically meaningful equation under investigation [22].It can be shown that, neglecting quantum effects, and considering macroscopic (thermodynamic) quantities, in most of the cases, the only plausible variables are the dimensionless pressure and temperature [19, 23,24]: where the dimensionless temperature and pressure are given by t = kT/ε and p = pa 3 /ε.Furthermore, in accordance with Equation ( 6), q should be the same function for all members of the similarity class in question.It was shown that, using scaling parameters kT m (∝ ε), Ω (∝ a 3 o ), m and k, the relationships derived from (8) provided nice agreement with experimental data at p ∼ = 0 for most metals and alloys (see, for instance, relationships for the diffusion and point defect properties [19, [23][24][25]), even for binary alloys using compositiondependent melting points, T m (x), molar volume, Ω(x), and mass m(x).As it follows from (8), for the activation energy of diffusion at p ∼ = 0 which is the well-known rule of thumb for self-diffusion in normal metals [23,25] (Q D is independent of T, according to the well-known Arrhenius-type T-dependence of the diffusion coefficient).On the other hand, for bcc metals showing phonon-softening behavior in the form of a curved Arrhenius function (anomalous behavior), introduction of one new dimensionless parameter, ξ, in the argument of Q D was needed [23,26].According to [26,27], the curved Arrhenius plot can be described by the following temperature dependence of Q D : where is the constant obtained in normal metals and T ′ o was called, in [26], "hypothetical critical temperature".In addition, it was shown in [26] that at a fixed (low) temperature the Q D T m ratio showed a linear dependence on the phonon-softening parameter, , where ξ is the ratio of the transversal and longitudinal sound velocities along the < 111 > direction, i.e., ξ ∼ T ′ o .Since phonon softening is a key characteristic for martensitic transformations in SMAs, let us also formally introduce this parameter, denoted by ξ too, for shape memory alloys.Thus, we will use (8) in the form at atmospheric pressure (t = T T m , p r ∼ = 0).Accordingly, the equilibrium transformation temperature, T o , can be given as i.e., t o = ϑ(t o , ξ).
Equation (13) means that a universal relationship should exist between ξ and t o .Thus, T m can be considered a good measure of ξ at t = t o and, in the following, ξ will be replaced by t o = T o T m .It is well known that T o has a strong composition dependence and thus t o should also have similar behavior.

Dependence of the Reduced Transformation Heat and Entropy on the Transformation Temperature
Let us first consider, in general, the t o -dependence of the heat of transformation, h = ∆H kT m , and transformation entropy and Here, we considered that ∆H and ∆S are usually independent of the temperature, t.Equation ( 14) means that h should depend universally on t o only and should have the same constant value at t o (0).Consequently, for instance, the η h = 1 h dh dt o derivative should also have the same constant value for all SMAs at t o (0).Thus, for shape memory alloys, we can plot |∆H| kT m versus t o = T o T m , for NiTi (from [21], with T m ∼ = 1583 K), CuAlBe alloys (from [3], with T m ∼ = 1353 K), and Ni 2 MnGa (from [15], with T m ∼ = 1403 K), as well as for Ti 50−x Ni 40+x Cu 10 (from [4], with T m ∼ = 1550 K) as it is shown in Figure 1.(In these plots |Q| = −Q A−→M ∼ = ∆H and T o ∼ = M s +A f 2 assumptions were implicitly assumed, where Q A−→M (< 0) is the transformation heat measured by DSC, and M s and A f are the martensite start and austenite finish temperatures, respectively).[7,15,16,31,33,34] B2/18R CuAlBe [2,3] 1353 0.20 CuAlNi data for elastic constants and their T-dependence are averages of those published in [31] and [37], which deviate from the given average values by about ±18%.It is worth mentioning that in a recent paper [31] the temperature dependence of c′ and c44 was investigated in the very near vicinity of To in NiTi, Ni2MnGa, and CuAlNi and, except for the result for NiTi, their data provides about an order of magnitude larger values for βc′ than those given in Table 1 above.Data with upper index * are from theoretical papers also cited in the first column.Data for αc′ indexed by ** are estimated from the empirical relationship proposed by Veringen and Delaey [38] (see also [41]).1).
From the results shown in Figure 1, it can be seen that the h(t o ), in the investigated parameter ranges, can be well approximated by straight lines.The corresponding slopes, as well as the values of η h at t o (0) (see also Table 1), although it would be expected that h has the same dependence on t o for all SMAs, are characteristically different.The main difference between the above four alloys is that they have different types of symmetry change during the martensitic transformation: (B 2 (bcc)/B19 ′ (monclinic), in NiTi, B 2 (bcc)/18R(rombohedral), in CuAlBe, L2 1 (bcc)/tetragonal, in Ni 2 MnGa (where the structure of the tetragonal phase can be complex, being non-modulated or modulated [7]), and B 2 (bcc)/B19(orhorombic) in Ti 50−x Ni 40+x Cu 10 , respectively.Thus, plausibly, we have to make a distinction between SMAs on the basis of the type of symmetry change during the martensitic transformation; groups of alloys having the same symmetry change form different similarity sub-classes.
Table 1.Experimental dimensionless parameters at t o (0) in different sub-classes of SMAs (the references are given in the first column).For calculation of the dimensionless values, the melting points were estimated from the corresponding phase diagrams (see also Table 2) and we assumed that the atomic volume is the same for all alloys (Ω = Ω NiTi = 8.4 × 10 −6 m 3 mol [28]).In the third and fifth columns, values of η h , calculated from the t o -dependence of h (from Equation ( 14) and Figure 1), as well as from Equation (21), respectively, are shown for comparison.The above classification is also supported by the following arguments since, by definition, the dimensionless transformation entropy is given by

Sub-Class/Alloy
and the transformation entropy, as it is well known [3], is also different for different subclasses, i.e., it depends on the structure (symmetry) of the martensite.Furthermore, the fact that ∆H kT m is approximately a linear function of T o T m would dictate that the reduced entropy, s, should be independent of t o (and thus from x) within a sub-class, and has different constant values for different sub-classes.In addition, the function h(t o ) should go through the origin.
It can be seen that it is quite well fulfilled for CuAlBe and Ti = 0.30.It can be seen that the DSC data and the above slopes indeed agree very well (the slopes are 0.16 and 0.30, respectively).In addition, the experimental data of [3,4] confirm that s is indeed constant, i.e., independent of t o , in these alloys, and s e ∼ = s (see also the Appendix A).In the case of NiTi and Ni 2 MnGa, the slopes provide different values than those of the experimental values of the entropies ) at t o (0) (see also the data in Table 1).In addition, the s e =

∆S exp k
values show approximately a linear dependence on t o , as illustrated in  Since h has composition dependence via its universal to-dependence, a general relationship should exist between its composition derivative and , as follows: According to this, the universal constants  can also be calculated from the ratio of the slopes, giving the composition dependence of h and to (see also Table 1).
It is worth emphasizing that the found approximately linear universal dependence of ℎ on to can also explain that, even if both ∆ and To have non-linear dependence on Since h has composition dependence via its universal t o -dependence, a general relationship should exist between its composition derivative and dt o dx , as follows: According to this, the universal constants η h can also be calculated from the ratio of the slopes, giving the composition dependence of h and t o (see also Table 1).
It is worth emphasizing that the found approximately linear universal dependence of h on t o can also explain that, even if both ∆H and T o have non-linear dependence on x, ∆H versus T o can be linear (see, e.g., Figure 2a,b in [4] for binary NiTi, where the x-dependence of both M s and ∆H is non-linear, but ∆H versus M s is approximately a linear function.(Regarding the validity of assumption M s ∼ = T o , see our comments in Section 3).
Before considering the corresponding expressions for the reduced (dimensionless) c ′ , c 44 , and A, it is worth recalling that besides the Boltzmann constant, the other three scaling parameters can also have composition dependence.But it can be easily accepted that the composition dependence of m, Ω, and T m is week as compared to the composition dependence of T o in SMAs.Indeed, in shape memory host alloys, the value of the atomic mass and volume is expected to have only a few percent change upon alloying.Furthermore, for the composition dependence of T m , we can take, as examples, the binary NiTi system as well as the CuAlBe systems, in which 1 . Thus, for the sake of simplicity, in the following, only T o will be considered composition dependent.

Derivation of General Relationship for the Composition Dependence of the Transformation Temperature
Let us consider the elastic constant, c ′ , given in the following form: where γ c ′ is a universal function of t and t o and, at According to Equations ( 18) and ( 19), the dimensionless c ′ is a universal function of the reduced temperature, t, and the phonon-softening parameter, t o , within a given sub-class of SMAs.Equation (20) corresponds to Equation (1) in dimensionless form, and the t o -dependence is the consequence of the phonon softening, with an additional conclusion that the constant is expected to be the same within a sub-class but different in different sub-classes.From (19) (considering that 1 According to (19), in the γ c ′ = const.Condition, the constant value is the limit of the bilinear function of γ c ′ (t, x) taken at t = t o (0).Introducing t ′ = t − t o as a new variable in the vicinity of t o , the t o -dependence stems from the t'-dependence, and it is given by Combination of (21) and ( 22) leads to This is just Equation (3) and, by the same arguments that were used in [1], Equation (4) is also valid, i.e., dT o dc ′ = − In addition, dt o should have the same constant value at t o (0) within a sub-class, i.e., Thus, the value of η c ′ can be calculated from experimental data on β and T m (see also Table 3).In addition, from ( 23) and, since the right-hand side is the same negative constant value within a sub-class, we get that the ratio of slopes expressing the composition dependence of t o and c ′ is constant, i.e., a relative increase in c ′ with increasing composition should be compensated by the same relative decrease in t o (β > 0).This is more definite statement than simply saying the lower α is, the higher dT o dx is (see [1]).
Table 3.Estimated parameters related to the composition dependence of the transformation temperature.The last column shows the experimental values of dT o dx for comparison.Since, in Table 2, except for the binary NiTi alloy, in all cases, β c 44 < 0 (i.e., c 44 did not show phonon softening) the original Ren-Otsuka relationship (Equation ( 3)) with β c ′ and α c ′ was used.For NiTi, the generalized relationship (Equation ( 29)) with β = β c4 − β c ′ and α = α c4 − α c ′ was used ((Equation (3) would even lead to a positive value for − α β , i.e., the predicted value for dT o dx would be wrong: − α β = 2533).As we summarized in the introduction, Equations (4) and thus ( 24) can be valid only if c ′ shows phonon softening while c 44 does not.The same results as above can be obtained by replacing c ′ with c 44 (with β now describing the phonon-softening-caused temperature dependence of c 44 ); in this case the softening of c ′ should be neglected.This conclusion is in line with the comments of [12], where it was mentioned that in (24) ". . .c can be either c 44 or c ′ . . .".

Sub-Class
Following the same procedure as above for c ′ we can write for the dimensionless anisotropy factor, A and respectively; (η A = 1 A dA dt o ) and η A and A = γ A are constants at t o (0).Thus, (see also (21) and ( 22)) finally, is obtained, where and Equation ( 29) is the generalized Ren-Otsuka relationship with generalized α and β.If the composition and temperature dependence of c 44 can be neglected (or normal, i.e., it does not show phonon softening and its contribution to the martensitic transformation can be neglected) then one gets back the original Equation ( 3).For phonon-softening systems, the temperature coefficient contains two contributions.The positive one describes the direct phonon-softening contribution while the small, usually negligible, negative one is related to the usual (normal) anharmonicity-related softening of the crystal, and the first one dominates in the vicinity of T o .It can be noted that the anharmonicity-related temperature and composition dependence was also neglected in [1].Furthermore, if both c ′ and c 44 have phonon softening in the above difference, the "normal" anharmonicity-related contributions approximately cancel out.
It is worth mentioning that the above approach, namely that the constancy of the anisotropy constant is the best starting point for finding a generalized relationship for the composition dependence of T o , can also be confirmed from the general form of the Landau expansion of the free energy, F, as used in [18].This paper, instead of using only two strains (basal plane shear, e 1 , basal plane shuffle, η in Equation ( 1) of [1]), contained three ones: e 1 , η, and the {001} < 110 > non-basal plane shear, e 2 (and c ′ , ω 2 η , and c 44 are the corresponding energy terms, respectively).Now it is easy to show that by minimizing F with respect to both e 2 and η strains, one can get, besides condition (1) (and ω 2 η (T, x) = conts.: see Equations (7a) and (7b) in [1]), that c 44 (T, x) has to also be constant at T=T o .Now, the ratio of Equation ( 1) and c 44 (T, x)=const.shows that A = γ A = const.at t o (0) (see Equation ( 27)).Of course, if there is no phonon softening in c 44 (i.e., if its t o -dependence can be neglected) then only the phonon softening of c ′ occurs, and the original Ren-Otsuka relationship can be approximately valid (like the case of the Ti 50−x Ni 40+x Cu 10 or CuAlBe alloy with a large and increasing anisotropy by approaching T o ).On the other hand, the NiTi alloys represent the other limit, when both c ′ and c 44 show phonon softening with a small value of A and a decreasing tendency of A with approaching T o [18] (i.e., η A > 0, see also Table 1 and the discussion below).

Comparison with Experimental Data
Some general features of the t o -dependence of the reduced characteristic quantities were already analyzed in the previous chapter, and it led to the conclusion that the host SMAs can be divided into sub-classes (having the same type of symmetry change during MT) within which the above quantities have the same constant values at t o (0).In this chapter, in addition to summarizing these, we also consider other SMAs to support the conclusions based on the data analyzed in Section 2.2.Furthermore, the reduced values of c ′ , c 44 , and A will be collected and compared with the data available in the literature.
It has to be noted that, in many publications, the composition dependence of T o and M s is taken to be the same, and similarly to [1], we can also assume it here.M s can be given as −∆s c , where d o and e o denote the first derivatives of the dissipative and elastic energies per unit volume, during the cooling process at the beginning of the transformation, and ∆s c is the entropy change per unit volume during cooling [42].The second term in M s , in fact, determines the dissipation and elastic energy accumulation, i.e., the above assumption means that we neglect the composition dependence of these terms, although, in a more refined treatment, this should be necessary to take into account, since, e.g., according to [8,9], the dissipative energy (the integral of d o ) also shows a composition dependence.It is also worth mentioning that, in general, the x-dependence of T o (x) and the transformation heat, ∆H(x), are not strictly linear, but have a small downward curvature [8,21,30], but for the sake of simplicity, we neglect this moderate x-dependence of the slopes.Furthermore, it is also worth emphasizing that, in light of the results obtained in the previous chapter, it is indeed not expected that the values of the slopes dt o dx = 1 T m dT o dx should have the same value even within a given sub-class, while, according to Equation ( 26 3

below).
Table 1 contains the values of those dimensionless constants at t o (0), which are predicted to be the same within a given subclass; h, η h , s, and 1 T m dT o dx ) used in the calculation of the data given in Table 1.It has to be noted that most of the experimental data suffer from relatively large errors for the elastic parameters (typically between 15% and 25%), and there is a lack of reliable data, especially for the composition dependence of c ′ and c 44 .Table 3 shows the estimated parameters related to the composition dependence of the transformation temperature.
Table 1 contains the summary of the parameters predicted to be the same within the five sub-classes, represented by the NiTi, Ti 45−x Ni 50+x Cu 5 , Ti 50−x Ni 40+x Cu 10 , Ni 2 MnGa, Cu-Al-Be, CuZn, CuZnAl, and Cu 68 Al 28 Ni 4 alloys.It can indeed be seen that the estimated values are characteristically different for the sub-classes.Furthermore, data for Cu-based alloys with B 2 /18R transformation (6th, 7th, and 8th rows) are rather similar, demonstrating that these quantities have the same constant values within a certain sub-class, as predicted.Regarding the CuAlNi (with B 2 /2H transformation), the constants are also not much different from the values of the above Cu-based alloys, suggesting that these two subclasses behave similarly.On the other hand, the observations concluded in [3] support that the CuAlNi belongs to a different sub-class.In Figure 5 of [3], where the transformation heat ( ∆H A→M ) was plotted versus the M s temperature, the slopes of the straight lines were slightly but definitely different for CuAlNi (1.59 J/molK) from the common slope belonging to the fitted line on the data of CuAlBe and CuZnAl (1.30 J/molK; see also Figure 3

below).
From the above slopes, dM s , 5.2 and 4.8, respectively.It can be seen, as expected, that 5.2 is in a very good agreement with the value given in Table 1 for CuAlBe from Equation (17) (4.9) while, for CuAlNi, the agreement is still acceptable (from Equation (17) 3.9 was obtained).While this example also illustrates the experimental scatter (which is still in the range of the differences between 5.2 and 4.9 as well as 4.8 and 3.9), since the difference of the above slopes obtained on the basis of large number of experimental data collected in [3] was definite, one can confirm that the CuAlNi belongs to a different sub-class.
Interestingly, in Ni 2 MnGa, the signs of both α and 1 T m dT o dx are even negative (see Table 2), and thus Equations ( 4) and (5) (as it was also mentioned in [7]) remain valid, since the sign of β is still the same as for other phonon-softening alloys (i.e., it is positive).
One additional comment supporting that bcc metals (and alloys) with phonon softening behave differently than the "normal" metals can be made.According to the above results, βT m should be a universal constant (at t = 1, i.e., at T m ) for all "normal" metals, while its value can be different for phonon-softening systems and, in addition, it should be different for different sub-classes of SMAs.The most salient result is that, indeed, β > 0, belonging to phonon-softening elastic constants in SMAs.On the other hand, it is negative, e.g., for Ag, Au, and Cu [43] (and βT m is approximately constant for temperatures larger than the Debye temperature for all "normal" metals [43]: ~−0.56).versus Ms plots are more commonly used in the analysis of experimental data (see, e.g., [3,4]), Figure 3a  βT m (columns 4th and 5th), as calculated at T o (0) from the experimental data given in Table 2.It can be seen that the signs in all cases are correct.Note (as it is also mentioned in the caption of Table 3) that for NiTi alloys, which are the most frequently used shape memory materials, only the generalized relationship provides the correct sign.Furthermore, as it can be seen from columns 5 and 6, the agreement between the value of − α β and the experimental data for dT o dx is also satisfactory, taking into account the uncertainties of the experimental values of α and β at present.As an example, we can mention the case of CuAlBe alloys.Here, the temperature dependence of the elastic constants is well known (as it is also shown in Table 2), and this slope does not change with the composition [2].On the other hand, the composition dependence of c ′ and c 44 at room temperature and at the transformation temperature (see Figure 2 in [3]) is remarkably different (α c ′ is about 10 at room temperature and about 0.7 at T o , or the composition dependence of A is given by η A = −9.3 as well as −1.0, respectively).In the Tables above, the room temperature value of α c ′ was taken, since the reported value at T o would lead to about an order of magnitude smaller value, although a value between 0.7 and 10 (and closer to 10) would lead to better agreement (α c ′ ∼ = 6.7 would lead to exact agreement between − α β and the experimental dT o dx value).Finally, it is also worth adding that the constancy of 1 αT m dT o dx provides an explanation of the conjecture already proposed in 1988 by Verlinder and Delaey [44]: "the M s temperatures of all the alloys can be correlated with an expression similar to that given for the composition dependence of c ′ . . ." i.e., 1   T m dT o dx ∼ α."In addition, they expressed that "similar calculations and conclusions as those presented in this paper for the two observations concerning the composition dependence of c ′ and M s could be made for the other alloy systems, providing the necessary experimental data are available".
Finally, since, on the |Q| kT m versus t o plots, the intercept I also depends on the position of the fitted t o -interval for systems in which the entropy has a linear t o -dependence, and thus T m plots.It can be seen in the compiled plots in Figure 3b that the only difference between the straight lines is that their slopes are different for different sub-classes of SMAs.Thus, this is a nice illustration of our prediction that the slopes of the  The latter two results will be important in understanding and classification of experimental results that will be obtained from new measurements on different SMAs.
as s e versus t o plots.In these alloys, there is a well-defined linear dependence of s e on t o , which can be attributed to other than vibrational contributions (which are fully related to the symmetry change during the martensitic transformation [3,15]) and, for instance, in [15], the above dependence was interpreted by magnetic contribution to s e in Ni 2 MnGa.Thus, we can write for such an intrinsic dependence:  A1, and the linear fit is made in the certain t o interval with t o median value (t o ∼ = 0.20 and 0.21 for NiTi as well as Ni 2 MnGa, respectively).Thus, the intercepts of the linearized |Q| kT m versus t o have no physical meaning (and depend on the position of the fitted range).In contrast, the slopes of the |Q| kT m versus t o as well as s e versus t o plots are relevant parameters, and, as one can check it from the data collected in Table A1, they form a consistent set, in accordance with Equations (A8) and (A9).It can be noted that the position of the minimum in Figure A1 (form (A9)) is given by t omin = 1 2 t oo − s oo β s , i.e., it is determined by the fitting parameters of the linear s e versus t o function.

Figure 1 .
Figure 1.Graphs showing: |Q| kT m versus T o T m in binary NiTi (a), in CuAlBe (b), in Ni 2 MnGa alloys (c), and in Ti 50−x Ni 40+x Cu 10 (d) (on the basis of data published in [3,4,15,21], respectively (see also the text).The slopes, are 0.83, 0.16, 1.77, and 0.30, respectively.Blank squares are the experimental points and the dotted line is the extrapolation of the straight line fitted.It can be seen that these are linear functions, and the slopes of |Q| kT m versus T o T m are 0.83, 0.16, 1.77, and 0.30 for NiTi, CuAlBe, Ni 2 MnGa, and Ti 50−x Ni 50+x Cu 10 , respectively.

Furthermore, η h (t o ( 0
NiTi, 5.2 for the CuAlBe, 59 for Ni 2 MnGa, and 4.2 for Ti 50−x Ni 40+x Cu 10 systems, respectively (see also Table 50−x Ni 40+x Cu 10 , where the |Q| kT m versus T o T m function goes through the origin, while the linear extrapolation of the fitted straight lines have definite intercept values for binary NiTi and binary Ni 2 MnGa alloys.In the case of CuAlBe and Ti 50−x Ni 40+x Cu 10 alloys, it also means that the slopes in Figure 1 should be equal to the (constant) entropy, as calculated from the DSC data at t o (0): s e = ∆S exp k = 0.15 and ∆S exp k

Figure 2 .
Figure 2. Thus, both the linear dependence of s e and the different slopes of the |Q| kT m versus t o plots from s e (t o (0)) can be related to (i); the s e (t o ) function is not constant, but the linear function of t o (s e = s eo + β s (t o − t o (0) ) and/or (ii), the approximations used when the |Q| kT m versus t o is plotted instead of ∆H kT m versus t o .As discussed in Appendix A, (i) has the dominating effect and |Q| kT m is a quadratic function of t o , and the slope of the fitted linear relationship between |Q| kT m and t o , in a certain interval, is given approximately by s e (t o (0)) + t o (0)β s = s eo + t o (0)β s .Materials 2024, 17, x FOR PEER REVIEW 8 of 18

Figure 2 .
Figure 2. Graphs showing: ∆S k versus T o T m in binary NiTi alloys [21] (a) and in Ni 2 MnGa alloys (b) [15].The slopes, d ∆S k dt o , are 3.0 and 6.4, respectively (see also the text).Blank squares are the experimental points and the dotted line is the extrapolation of the straight line fitted.
have the same value (see also Table

Figure 3 .
Figure 3. Graphs showing | | versus (a) as well as | | −  versus (b) plots (where I denotes the values of the intercepts in Figure 1 (see also the text).See also Figures 1 and 2. Finally, since, on the | | versus to plots, the intercept I also depends on the position of the fitted to-interval for systems in which the entropy has a linear to-dependence, and thus | | versus to is the quadratic function (see the Appendix A), it would be worthwhile to compile these plots into a common plot of | | −  versus to.Furthermore, since the | | ,b shows the | | versus to as well as | | −  versus plots.It can be seen in the compiled plots in Figure 3b that the only difference between the straight lines is that their slopes are different for different sub-classes of SMAs.Thus, this is a nice illustration of our prediction that the slopes of the | | versus to plots should be different for SMAs different symmetry changes during MTs.It is so, even if one takes into account that: (i) In those systems where the entropy has an intrinsic to-dependence, the slopes differ from the  ( (0)) = ∆ ≅ | | values (shown in the fifth column ofTable 1); and (ii) the slopes of | | −  versus are obviously slightly different from those of the | | −  versus to plots (for instance the slopes are 1.11 as well 0.83 in NiTi or 0.160 and 0.156 in CuAlBe, respectively).

Figure 3 .
Figure 3. Graphs showing |Q| kT m versus M s T m (a) as well as |Q| kT m − I versus M s T m (b) plots (where I denotes the values of the intercepts in Figure 1 (see also the text).See also Figures 1 and 2.Table 3 contains the comparison of the predicted values of η −1 A = 1 αT m dT odx as well as η −1 A = − 1 βT m (columns 4th and 5th), as calculated at T o (0) from the experimental data given in Table2.It can be seen that the signs in all cases are correct.Note (as it is also mentioned in the caption of Table3) that for NiTi alloys, which are the most frequently used shape memory materials, only the generalized relationship provides the correct sign.Furthermore, as it can be seen from columns 5 and 6, the agreement between the value of − α β and the experimental data for dT o dx is also satisfactory, taking into account the uncertainties of the experimental values of α and β at present.As an example, we can mention the case of CuAlBe alloys.Here, the temperature dependence of the elastic constants is well known (as it is also shown in Table2), and this slope does not change with the composition[2].On the other hand, the composition dependence of c ′ and c 44 at room temperature and at the transformation temperature (see Figure2in[3]) is remarkably different (α c ′ is about 10 at room temperature and about 0.7 at T o , or the composition dependence of A is given by η A = −9.3 as well as −1.0, respectively).In the Tables above, the room temperature value of α c ′ was taken, since the reported value at T o would lead to about an order of magnitude smaller value, although a value between 0.7 and 10 (and closer to 10) would lead to better agreement (α c ′ ∼ = 6.7 would lead to exact agreement between − α β and the experimental

|Q|
kT m versus t o is the quadratic function (see the Appendix A), it would be worthwhile to compile these plots into a common plot of |Q| kT m − I versus t o .Furthermore, since the |Q| kT m versus M s plots are more commonly used in the analysis of experimental data (see, e.g.,[3,4]), Figure3a,b shows the |Q| kT m versus t o as well as |Q| kT m − I versus M s that the application of the law of corresponding states for martensitic transformations of shape memory alloys with phonon softening requires the introduction of a new dimensionless phonon-softening parameter, which is proportional to t o = T o T m .-Both the dimensionless heat and entropy of transformation (h = ∆H kT m and s = ∆S k ) are universal functions of t o , and the composition dependences of them are determined by the composition dependence of t o (or T o , since the composition dependence of T m can be neglected).-The slopes of the linearized h versus t o plots were different for SMAs with different symmetry changes during martensitic transformation, forming sub-classes.-Within a given sub-class, the normalized parameters, like the c ′ elastic constant or the anisotropy constant (γ = c ′ Ω kT m and A = c 44 c ′ ), are the same constants at T o T m = t o (0).-From the above property of A, the generalized Ren-Otsuka relationship is obtained with generalized α and ββ parameters (α = α c4 − α c ′ = as β = β c4 − β c ′ = , where these are different from zero only for parameters showing phonon softening).-It is shown that 1 α dt o dx is the same constant within a given sub-class.-The obtained a linear relationship between ∆H and T o rationalizes the observed empirical linear relationships between the heat of transformation measured by DSC (Q A−→M ) and the martensite start temperature, M s .-

|Q|
kT m versus t o as well as s e versus t o plots for NiTi and Ni 2 MnGa alloys.Materials 2024, 17, x FOR PEER REVIEW 16 of 18 linearized | | versus to have no physical meaning (and depend on the position of the fitted range).In contrast, the slopes of the | | versus to as well as  versus to plots are relevant parameters, and, as one can check it from the data collected in Table A1, they form a consistent set, in accordance with Equations (A8) and (A9).It can be noted that the position of the minimum in Figure A1 (form (A9)) is given by  =  − , i.e., it is determined by the fitting parameters of the linear se versus to function.

Table A1 .Figure A1 .
Figure A1.Graphs showing | | versus to plots for NiTi (a) and Ni2MnGa (b).The fitted curves show the fit with Equation (A9) and both fitting parameters are in very good agreement with data given in Table A1 (obtained from the linear of se versus to):  −   = −0.15 and  −   = −0.75 as well as  = 2.9 and  = 6.4,for NiTi and Ni2MnGa systems, respectively.The same is valid for the positions of the minima.Blank squares are the experimental points and the dotted line demonstrate the fitted quadratic function.

Figure A1 .
Figure A1.Graphs showing |Q| kT m versus t o plots for NiTi (a) and Ni 2 MnGa (b).The fitted curves show the fit with Equation (A9) and both fitting parameters are in very good agreement with data given in Table A1 (obtained from the linear of s e versus t o ): s eo − β s t oo = −0.15 and s eo − β s t oo = −0.75 as well as β s = 2.9 and β s = 6.4,for NiTi and Ni 2 MnGa systems, respectively.The same is valid for the positions of the minima.Blank squares are the experimental points and the dotted line demonstrate the fitted quadratic function.

β c ' T m β c 44 T m α c ' α c '
|Q| kT m versus t o plots should be different for SMAs with different symmetry changes during MTs.It is so, even if one takes into account that: (i) In those systems where the entropy has an intrinsic t o -dependence, the slopes differ from the s e (t o (0)) = − I versus M s T m are obviously slightly different from those of the |Q| kT m − I versus t o plots (for instance the slopes are 1.11 as well 0.83 in NiTi or 0.160 and 0.156 in CuAlBe, respectively).
s eo + β s (t o − t oo ) = s eo − β s t oo + β s t o , (A8) where s eo = s e (t o (0)) and t oo = t o (0) are constants.Putting (A8) into (A5) we get |Q| kT m = t o [s eo + β s (t o − t oo )] = (s eo − β s t oo )t o + β s t 2 = s eo − β s t oo + 2β s t o = s eo + β s (2t o − t oo ), (A10) where t o is the median value of the interval where the linear fit was made to the |Q| kT m versus t o plot.It can be seen from (A9) that |Q| kT m is a quadratic function of t o , as is illustrated in Figure

Table A1 .
Slope and intercepts of the