Jump to ContentJump to Main Navigation
Show Summary Details
More options …

Open Physics

formerly Central European Journal of Physics

Editor-in-Chief: Seidel, Sally

Managing Editor: Lesna-Szreter, Paulina

IMPACT FACTOR 2018: 1.005

CiteScore 2018: 1.01

SCImago Journal Rank (SJR) 2018: 0.237
Source Normalized Impact per Paper (SNIP) 2018: 0.541

ICV 2017: 162.45

Open Access
See all formats and pricing
More options …
Volume 16, Issue 1


Volume 13 (2015)

Modeling hysteresis curves of La(FeCoSi)13 compound near the transition point with the GRUCAD model

Roman Gozdur
  • Corresponding author
  • Łódź University of Technology, Department of Semiconductor and Optoelectronics Devices, Wólczańska 211/215, 90-924, Łódź, Poland
  • Email
  • Other articles by this author:
  • De Gruyter OnlineGoogle Scholar
/ Piotr Gębara / Krzysztof Chwastek
Published Online: 2018-05-24 | DOI: https://doi.org/10.1515/phys-2018-0038


The paper presents a new approach to modeling of magnetic properties and the influence of temperature on hysteresis loops in La(FeCoSi)13 magnetocaloric compound. The GRUCAD hysteresis model has been used for analysis of magnetic behavior in the vicinity of the magnetic phase transition. The modeling results are in a qualitative agreement with experiment.

Keywords: LaFeCoSi; magnetocaloric materials; GRUCAD model; magnetic power losses; specific power loss measurements

PACS: 75.60.Ej; 75.30.Sg

1 Introduction

Modern cooling based on magnetocaloric effect is an environment-friendly technique used to lower the temperature. The well-known magnetocaloric materials are pure gadolinium and its alloys, especially Gd5Ge2Si2 discovered by Pecharsky and Gschneidner Jr. in 1997 [1]. Due to the fact that Gd is lanthanide, its price is relatively high. Relatively cheap pseudo-binary La(FeCoSi)13 compounds have been the subject of considerable interest of engineering community as the magnetocaloric effect occurs in these materials near room temperature [2, 3]. This group of alloys is based on face centered cubic NaZn13 structure and it is responsible for excellent magnetic properties. The La atoms occupy 8a (in Wyckoff notation) positions and Fe atoms two non-equivalent positions (8b and 96i). The La atoms may be substituted by other lanthanides (Ce [4], Pr or Ho [5]). Moreover, the Fe atoms placed in position 96i may be substituted by i.e. Co [6], Ni [7] or Al [8]. Such modifications of chemical composition have an important influence on the value of the Curie point and change in magnetic entropy.

Most of the papers concerning the La(Fe,Si)13 alloys deliver information about their processing, microstructure and magnetocaloric properties. From the scientific point of view, it is also important to investigate various effects in the vicinity of phase transition observed in the material. An interesting technique used to study the response of magnetic material near the phase transition is modeling of hysteresis curves [9].

The present paper considers the possibility to model the hysteresis curves of La(FeCoSi)13 compound near the transition temperature with the use of a phenomenological description. It is an extension of the previous work [10], where the variations of shape of the anhysteretic (hysteresis-free) curves were considered.

The paper is structured as follows. Section 2 provides basic pieces of information on contemporary hysteresis models used in studies of magnetocaloric effect and the GRUCAD description applied for this purpose in this paper. Section 3 discusses in detail the measurement methods and results. Section 4 is devoted to modeling. Finally, Section 5 comments on the obtained results.

2 Phenomenological hysteresis models

Hysteresis modeling may be carried out with an appropriate description like the bottom-up approach by Preisach, extended by Mayergoyz [11] or the phenomenological proposal advanced by Jiles and Atherton [12]. The latter model has attracted much attention of the scientific community in the last thirty years due to the uncomplicated structure of its equations and physical meaning attributed to its parameters. The Jiles-Atherton (JA) model still remains one of the most used descriptions for modeling magneto-structural first order transitions in magneto-caloric materials as well [13]. However, some authors have pointed out a number of its drawbacks, to mention the necessity to “patch” the model equations with an artificially introduced term suppressing the irreversible magnetization term after a sudden field reversal [14], non-closure of minor hysteresis loops if the values of model parameters are kept constant [15], the ambiguities due to the use of modified Langevin function as the description of the anhysteretic curve [16, 17] (“modified” means that the argument of the function includes the mean field term).

The model developed by the GRUCAD research group [18, 19] is a valuable alternative to the original JA formulation. In the GRUCAD description the reversible and irreversible processes are decoupled. The irreversibility is implemented by offsetting the loop branches along the H-axis from the anhysteretic (truly reversible in thermodynamic sense) curve, given with Langevin function. There is no need to introduce additional conditions, which are aimed at the modification of model equations after a field reversal. This description is similar in spirit to some models consistent with the laws of irreversible thermodynamics [20, 21, 22, 23]. The equations that define the model are as follows:






Here δ = sign (dB/dt) = ± 1, whereas α, a, γ, HHS and Ms are model parameters. Their values may be obtained using any robust optimization routine. Two first equations refer to the reversible, whereas the third and the fourth one – to the irreversible magnetization processes.

3 Material and measurements

Magnetocaloric bulk samples of La(FeCoSi)13 applied for the tests are made of reactively sintered fine-grained powders manufactured in Vacuumschmelze Company (Figure 1) [24]. Curie temperature Tc and peak entropy changes ΔSM(T) can be tailored through introduction of cobalt into structure La(FexCoySi1−xy)13 [3, 24]. The study and modeling were carried out with composition tuned-up for the Curie point of 300 K (Figure 2).

SEM micrograph of a fracture of sintered La(FeCoSi)13 samples
Figure 1

SEM micrograph of a fracture of sintered La(FeCoSi)13 samples

Influence of temperature and magnetic field strength on magnetic entropy changes of La(Fe10.8Co1.1Si1.1) composition
Figure 2

Influence of temperature and magnetic field strength on magnetic entropy changes of La(Fe10.8Co1.1Si1.1) composition

The final-form of plates with overall dimension 18 × 36 × 5 mm has been obtained in machining process [24]. Nominal parameters of examined composition are comprised in Table 1. Temperature and magnetic dependencies of magnetic entropy changes |ΔSM(T)| are depicted in Figure 2.

Table 1

Parameters of the tested magnetocaloric material

The experimental study of magnetic properties concerns hysteresis loops, magnetic polarization and power losses in the vicinity of Curie temperature. The presented results of measurements have been obtained with application of IEC 60404-6:2004 method. The reference measurements of power losses have been performed simultaneously with the use of Unbalanced Bridge Method (UBM) [25]. Instrumentation and primal uncertainties of the measurement methods are presented in Table 2.

Table 2

Accuracy of the measurement system

Measurements of magnetic waveforms H(t), J(t), uB(t) and specific power losses PS have been carried out in a closed magnetic core. The core was made by setting of four homogeneous LaFe10.8Co1.1Si1.1 (total mass = 0.094 kg) plates in non-magnetic yoke. Each component of the core was equipped with a section of pickup coil W2 (4 × 150 turns) and exciting coil W1 (4 × 72 turns) (see Figure 4). Comprehensive description of the measurement setup was presented previously in Ref. [26], whereas the most significant uncertainty components are collected in Table 2.

Investigations of thermal effects on magnetic properties and power losses were evaluated through analysis of the magnetic hysteresis loops J(H) measured at the frequency of 0.1 Hz (quasi-static conditions - cf. Figure 3).

Influence of temperature on the hysteresis loops at given peak field strength 1900 A/m near the transition point (Tc = 300K).
Figure 3

Influence of temperature on the hysteresis loops at given peak field strength 1900 A/m near the transition point (Tc = 300K).

Transitional stage of magnetic coercive field HC, remanent magnetization Jr and saturation Jmax of La(FeCoSi)13 core are shown in Figure 4. The measurements confirm a strong influence of the temperature on hysteresis loops of La(FeCoSi)13 compound.

Magnetic properties of LaFe10.8Co1.1Si1.1 near the transition point (Tc = 300 K).
Figure 4

Magnetic properties of LaFe10.8Co1.1Si1.1 near the transition point (Tc = 300 K).

The measured hysteresis loops (Figure 3) and the magnetization curve Jmax(T) in Figure 4 clearly indicate a transition from ferromagnetic to paramagnetic regime close to 300 K.

4 Modeling

In the previous work [10], we have focused on the description of the anhysteretic curve near the transition temperature. The shape parameter in the Langevin function used in modeling was found to be directly related to temperature. This fact allowed us to obtain a unified description for anhysteretic curves for temperatures both below and above the transition point. For temperatures above 300 K the hysteresis loops practically disappeared and the slope of the anhysteretic curve changed abruptly. In the present paper we focus on modeling hysteresis curves for temperatures below 300 K.

The values of GRUCAD model parameters were determined for successive temperatures, keeping the values of fixed to those determined previously [10]. The results of modeling are presented in Figure 5.

The modeled hysteresis curves for temperatures below phase transition
Figure 5

The modeled hysteresis curves for temperatures below phase transition

In order to check how well the GRUCAD model is able to describe the hysteresis loops of La(FeCoSi)13 compound near the transition temperature, we have compared the measured and the modeled values of coercive field strength and remanence induction for quasi-static magnetization conditions (f = 0.1 Hz). Moreover we have integrated numerically the area of the modeled hysteresis loops in order to find out the approximate value of power loss density. The results for chosen temperatures are summarized in Table 3.

Table 3

A comparison of measured and modeled quantities

It can be stated, that the measured and modeled values of coercive field strength and power loss density are in qualitative agreement. The maximum relative error understood as the absolute value of the difference between the measured and the modeled values of the quantity referred to the measured value does not exceed 17.1%. However it can be noticed that the model overestimated significantly the remanence value at 295.5 K. Figure 6 depicts the hysteresis loops for temperatures considered in Table 3 in a more detailed way.

The hysteresis curves for temperatures 295.5 K and 299.5 K. Dots – measurement data
Figure 6

The hysteresis curves for temperatures 295.5 K and 299.5 K. Dots – measurement data

5 Conclusions

The presented modeling results clearly indicate that the GRUCAD model might be a valuable tool supporting prediction of the effects of temperature variations in La(FeCoSi)13 magnetocaloric compound near the transition temperature. The obtained results are in qualitative agreement with the experiment and for coercive field strength and power loss density there is also an excellent quantitative agreement.

Future work shall be aimed at fine-tuning of the GRUCAD description. In the present contribution it was assumed that the compound may be described with a single set of model parameters. It should be born in mind that in reality magneto-caloric materials exhibit a highly complicated internal structure due to the presence of multiple phases. Therefore a more adequate description should include distributions of model parameters to take this fact into account.


This research was partly supported by the grants No. 6370/B/T02/2011/40 and 4889/B/T02/2011/40, from the National Science Centre of Poland.


  • [1]

    Tishin A.M., Spichkin Y.I., Themagnetocaloric effect and its application, Institute of Physics Publishing, London, 2003 Google Scholar

  • [2]

    Liu J., Moore J. D., Skokov K. P., Krautz M., Löwe K., Barcza A., et al., Exploring La(Fe,Si)13-based magnetic refrigerants towards application, Scripta Mat., 2012, 67, 584-589. CrossrefGoogle Scholar

  • [3]

    M’nassri R., Cheikhrouhou A., Magnetocaloric effect in LaFe10.7Co0.8Si1.5 compound near room temperature, J. Supercond. Nov. Mater., 2014, 27, 1059-1064. CrossrefGoogle Scholar

  • [4]

    Fujita A., Fujieda S., Fukamichi K., Realtive cooling power of La(FexSi1−x)13 after controlling the Curie temperature by hydrogenation and partial substitution of Ce, J. Magn. Magn. Mat., 2007, 310, e1006-e1007. Google Scholar

  • [5]

    Gębara P., Kovac J., Magnetocaloric effect of the LaFe11.2Co0.7Si1.1 modified by partial substitutiton of La by Pr or Ho, J. Mat. Des., 2017, 129, 111-115. Google Scholar

  • [6]

    Yan A., Muller K.H., Gutfleisch O., Magnetocaloric effect in the LaFe11.8−xCoxSi1.2 melt-spun ribbons, J. All. Compd., 2008, 450, 18-21. CrossrefGoogle Scholar

  • [7]

    Gębara P., Pawlik P., Hasiak M., Alteration of negative lattice expansion of the La(Fe,Si)13-type phase in LaFe11.14−xCo0.66NixSi1.2 alloys, J. Magn. Magn. Mat., 2017, 422, 61-65. CrossrefGoogle Scholar

  • [8]

    Gębara P., Pawlik P., Michalski B., Wysłocki J.J., Measurements of magnetocaloric effect in LaFe11.14Co0.66Si1.2-xAlx (x = 0.1, 0.2, 0.3) alloys, Acta Phys. Pol. A, 2015, 127, 576-578. CrossrefGoogle Scholar

  • [9]

    Basso V., Kuepferling M., Sasso C. P., LoBue M., Modeling hysteresis of first-order magneto-structural phase transformations, IEEE Trans. Magn., 2008, 44, 3177-3180. CrossrefWeb of ScienceGoogle Scholar

  • [10]

    Gozdur R., Chwastek K., Najgebauer M., Lebioda M., Bernacki Ł., Wodzynski A., Scaling of anhysteretic curves for LaFeCoSi alloy near the transition point, Acta Phys. Pol. A, 2017, 131, 801-803. CrossrefWeb of ScienceGoogle Scholar

  • [11]

    Mayergoyz I.D., Mathematical models of hysteresis and their applications, 2nd Ed., Elsevier Academic Press, Amsterdam, 2003 Google Scholar

  • [12]

    Jiles D.C., Atherton D.L., Theory of ferromagnetic hysteresis, J. Magn. Magn. Mater., 1986, 61, 48-60. CrossrefGoogle Scholar

  • [13]

    Melikhov Y., Hadimani R.L., Raghunathan A., Phenomenological modelling of first order transitions in magnetic systems, J. Appl. Phys., 2014, 115, 183902 Web of ScienceCrossrefGoogle Scholar

  • [14]

    Steentjes S., Chwastek K., Petrun M., Dolinar D., Hameyer K., Sensitivity analysis and modeling of symmetric minor hysteresis loops using the GRUCAD description, IEEE Trans. Magn., 2014, 50, 7300804 Web of ScienceGoogle Scholar

  • [15]

    Benabou A., Leite J. V., Clénet S., Simão C., Sadowski N., Minor loops modelling with a modified Jiles-Atherton model and comparison with the Preisach Model, J. Magn. Magn. Mater., 2008, 320, e1034-e1038. Web of ScienceGoogle Scholar

  • [16]

    Zirka S. E., Moroz Yu. I., Harrison R. G., Chwastek K., On physical aspects of Jiles-Atherton models, J. Appl. Phys., 2012, 112, 043916 CrossrefWeb of ScienceGoogle Scholar

  • [17]

    Qingyou Liu, Xu Luo, Haiyan Zhu, Yiwei Han, Modified magnetomechanical model in the constant and low intensity magnetic field based on J-A theory, Chinese Phys. B, 2017, 26, 077502 CrossrefGoogle Scholar

  • [18]

    Koltermann P. I., Righi L. A., Bastos J. P. A., Carlson R., Sadowski N., Batistela N. J., A modified Jiles method for hysteresis computation including minor loops, Physica B, 2000, 275, 233-237. CrossrefGoogle Scholar

  • [19]

    Righi L. A., Sadowski N., Carlson R., Bastos J. P. A., Batistela N. J., A new approach for iron losses calculation in voltage fed time stepping Finite Elements, IEEE Trans. Magn., 2001, 37, 3353-3356. CrossrefGoogle Scholar

  • [20]

    Bergqvist A., Magnetic vector hysteresis model with dry friction-like pinning, Physica B, 1997, 233, 342-347. CrossrefGoogle Scholar

  • [21]

    Henrotte F., Nicolet A., Hameyer K., An energy-based vector hysteresis model for ferromagnetic materials, COMPEL, 2006, 25, 71-80. CrossrefGoogle Scholar

  • [22]

    Tavakoli H., Bormann D., Ribberfjärd D., Engdahl G., Comparison of a simple and a detailed model of magnetic hysteresis with measurements on electrical steel, COMPEL, 2009, 28, 700-710. CrossrefWeb of ScienceGoogle Scholar

  • [23]

    François-Lavet V., Henrotte F., Stainer L., Noels L. , Geuzaine C., An energy-based variational model of ferromagnetic hysteresis for finite element computations, J. Comp. Appl. Math., 2013, 246, 243-250. CrossrefGoogle Scholar

  • [24]

    Katter M., Zellmann V., Reppel G.W., Uestuener K., Magnetocaloric Properties of La(Fe,Co,Si)13 Bulk Material Prepared by Powder Metallurgy, IEEE Trans. Magn., 2008, 44, 3044-3047. Web of ScienceCrossrefGoogle Scholar

  • [25]

    Gozdur R., Majocha A., Power losses measurements of nanocrystalline and amorphous magnetic cores, Electr. Machin. - Trans. J., 2013, 100, 175-179. Google Scholar

  • [26]

    Gozdur R., Lebioda M., Bernacki Ł., Power losses in LaFexCoySi1.1 intermetallics near the magnetic phase transition, Acta Phys. Pol. A., 2015, 128, 98-103. CrossrefGoogle Scholar

About the article

Received: 2017-11-01

Accepted: 2017-11-30

Published Online: 2018-05-24

Citation Information: Open Physics, Volume 16, Issue 1, Pages 266–270, ISSN (Online) 2391-5471, DOI: https://doi.org/10.1515/phys-2018-0038.

Export Citation

© 2018 Roman Gozdur et al., published by De Gruyter. This work is licensed under the Creative Commons Attribution-NonCommercial-NoDerivatives 4.0 License. BY-NC-ND 4.0

Comments (0)

Please log in or register to comment.
Log in