This paper was converted on www.awesomepapers.org from LaTeX by an anonymous user.
Want to know more? Visit the Converter page.

Anomalous depinning of magnetic domain walls within the ferromagnetic phase of the Weyl semimetal Co3Sn2S2

Zihao Shen Department of Physics and Astronomy, University of California, Davis, California 95616, USA    X.D. Zhu Department of Physics and Astronomy, University of California, Davis, California 95616, USA Department of Optical Sciences and Engineering, Fudan University, Shanghai 200045, China [email protected]    Rahim R. Ullah Department of Physics and Astronomy, University of California, Davis, California 95616, USA    Peter Klavins Department of Physics and Astronomy, University of California, Davis, California 95616, USA    Valentin Taufour Department of Physics and Astronomy, University of California, Davis, California 95616, USA [email protected]
Abstract

We report bulk magnetization measurements and spatially resolved measurements of magnetic domains in Co3Sn2S2\mathrm{Co}_{3}\mathrm{Sn}_{2}\mathrm{S}_{2} single crystals. The results indicate that a previously reported magnetic anomaly around 130130 K is due to an anomalous domain wall depinning upon cooling. Our measurements also reveal a hysteresis between field-cooled-cooling (FCC) and field-cooled-warming (FCW) magnetization curves acquired under a constant magnetic field below 300 Oe. This observation rules out the possibility that the anomaly stems from a second-order phase transition. Our results further suggest that changes in the shape of hysteresis loops from 55 K to 170170 K are caused by an unusual temperature-dependent domain nucleation field that changes sign around 130130 K. The Kerr rotation images of the magnetic domains confirm that the domain walls depin between 120120 K and 140140 K.

\ioptwocol

1 Introduction

Weyl Semimetals (WSM) are material systems with non-degenerate bands that cross each other near the Fermi level when inversion symmetry or time-reversal symmetry is broken [1, 2, 3]. As one of the most interesting examples of WSM [4], Co3Sn2S2 is a ferromagnetic material with a Curie temperature of TC=174T_{\mathrm{C}}=174 K [5]. The topological features of this material lead to a large anomalous Hall conductivity [6], an anomalous Nernst Hall effect [7], a negative magnetoresistance due to a chiral anomaly [4], and an unusual breaking-down of Ohm’s law [8]. In addition, an unusually large magneto-optic effect has also been observed [9].

Co3Sn2S2 crystallizes in a hexagonal lattice with triangular S and Sn layers interspersed between the kagome lattice planes of Co. The magnetic structure of this compound is still under investigation due to an anomaly in its magnetization occurring around 130130 K upon cooling [10, 11, 12, 13, 14, 15, 16]. Various interpretations have been proposed to explain the anomaly. Wu et al. proposed a hidden magnetic phase in Co3Sn2S2 as revealed in AC susceptibility measurements [17]. Guguchia et al. suggested an antiferromagnetic phase (AFM) in this compound based on two different muon precessions below TCT_{\mathrm{C}} [18]. Meanwhile, other groups interpreted the anomaly as a spin glass transition [19, 20], a domain wall transition from linear domain walls to elliptical domain walls [21], and an enhanced domain wall motion observed by magnetic force microscopy [22]. Some of these proposals, however, have been challenged by experiments. For example, no strong evidence for an AFM phase was found in a neutron scattering experiment [13].

In this context, we report bulk magnetization measurements on Co3Sn2S2 single crystals in external magnetic fields applied along the cc axis, where we find a depinning of magnetic domain walls around 130130 K upon cooling and a hysteresis between field-cooled-cooling (FCC) and field-cooled-warming (FCW) magnetization at low fields. We also present magneto-optic Kerr effect (MOKE) images of magnetic domains in Co3Sn2S2 from 120120 to 140140 K obtained under an external magnetic field of 9494 Oe applied along the cc axis. From the MOKE images, we extract fractions of single domains and domain walls near the magnetic anomaly around 130130 K. While we find no strong evidence of sudden domain reorganization or domain size changes, we do find a change in domain fractions consistent with the partial depinning of magnetic domains. Our results indicate that the magnetic anomaly can be explained by a domain wall depinning effect. The origin of the domain wall depinning effect remains to be further investigated.

2 Methods

2.1 Single Crystal Growth

Single crystals of Co3Sn2S2 were synthesized by solution growth [8, 23, 24]. A ternary mixture with initial composition Co12S8Sn80 was first heated to 400400°C over two hours and held there for another two hours. It was then heated to 10501050°C over six hours and held there for 1010 hours, followed by a slow, 9090 hour cool down to 740740°C. The remaining molten flux was removed by centrifugation. Shiny hexagonal crystals were obtained. The powder x-ray diffraction data refines to a hexagonal unit cell with a=5.3641(8)a=5.3641(8) Å and c=13.1724c=13.1724 Å, which are consistent with previously reported values [25].

2.2 Magnetization measurements and magnetic domain imaging

Magnetization measurements of Co3Sn2S2 single crystals were performed in a Quantum Design MPMS SQUID up to 77 T with the magnetic field parallel to the cc axis. Magnetization as a function of temperature was measured under various applied fields in both field-cooled-cooling (FCC) and field-cooled-warming (FCW) conditions. Hysteresis loops were measured by first cooling down to the desired temperature in zero applied field and then sweeping the field. Before each measurement, we followed a demagnetization procedure to minimize the remanent field in the magnet and in the sample chamber by systematically oscillating the magnet to zero field at room temperature.

Polar MOKE images are acquired from as-grown Co3Sn2S2(0001) single crystal samples, using a normal-incidence Sagnac interferometric scanning microscope  [26]. The microscope measures Kerr rotation as opposed to Kerr ellipticity [9]. It has a spatial resolution of 0.85 μ\mum and a sensitivity of 4.44.4μ\murad. The single crystal sample is mounted on a cold finger inside a cryostat. The sample is optically accessible from outside through an optical window. The temperature of the sample is variable from 300300 K to 55 K. We use a permanent magnet placed outside of the cryostat to produce a variable magnetic field of up to 100 Oe (0.01 T) along the cc axis (the surface normal) of the sample. For this experiment, the field is held constant at 94 Oe, measured by a Hall effect sensor.

3 Results and Discussion

3.1 Magnetization measurements

The FCC curve with a magnetic field of 10001000 Oe applied parallel to the cc axis is shown in Fig. 1, and illustrates the ferromagnetic-paramagnetic transition. The inverse susceptibility is fit to the Curie-Weiss law and we obtain an effective moment of 1.141.14μB\mu_{B} per Co and a Curie-Weiss temperature of θCW=175\theta_{\mathrm{CW}}=175 K, which are consistent with previously published work [23, 25, 27, 4]. There is no difference in the FCC and FCW magnetization at fields above 300300 Oe. If the applied field is below 300300 Oe, however, two anomalies emerge. First, the FCC and FCW curves start to deviate near TCT_{\mathrm{C}}. Second, a local maximum arises at TAT_{\mathrm{A}} \approx 130 K in the FCC data, which is also reported in other studies [18, 10, 11, 12, 13, 14, 15, 16]. The local maximum disappears above 300 Oe and below 50 Oe, while a small hysteresis between FCC and FCW still remains at 50 Oe. The field dependence of TAT_{\mathrm{A}} is relatively small, as shown in the inset of Fig. 1(b). The fact that the anomaly disappears above 300300 Oe, where the magnetization is saturated, suggests its origin is related to magnetic domains. A local minimum also emerges in the FCW magnetization albeit at higher temperature. The hysteresis between FCC and FCW magnetizations indicates it is unlikely that the local maximum and minimum in the FCC and FCW magnetization are due to a second order phase transition, which should not depend on magnetic history.

Refer to caption
Figure 1: a) Magnetization as a function of temperature measured by FCC in a 10001000 Oe applied field parallel to the cc axis. The blue curve shows a Curie-Weiss fitting with an effective moment 1.141.14μB\mu_{B} per Co. The inset shows a Co3Sn2S2 single crystal grown from solution growth. b) FCC and FCW curves show a hysteresis between 100100 K and TC=174T_{\mathrm{C}}=174 K. TAT_{A} is the temperature of the magnetic anomaly where the FCC magnetization reaches a local maximum followed by a downward kink. The inset shows the field dependence of TAT_{A}.

We will compare the fraction of magnetic domains assuming no pinning with the fraction extracted from our experimental magnetization where pinning effects are inevitable. In the absence of domain wall pinning, the demagnetising field Hd=NcMH_{d}=N_{c}M perfectly counteracts the applied field, HappliedH_{\mathrm{applied}}, as long as Happlied<NcMsH_{\mathrm{applied}}<N_{c}M_{\mathrm{s}}. NcN_{c} is the demagnetization factor along the cc-direction. Because the cc axis is the easy magnetization axis, NcN_{c} can be measured experimentally from a plot of M2M^{2} versus Happlied/MH_{\mathrm{applied}}/M for the M(H)M(H) data along the easy axis of magnetization as shown in Fig. 2(c). The x-intercept gives the value of the demagnetization factor along the easy magnetization axis [28, 29]. Due to the existence of domain pinning, the curves at low temperatures tend to shift to the right. Therefore, we used the x-intercept at 170170 K, where domain pinning is negligible, as the demagnetization factor and obtained Nc=0.608N_{c}=0.608. MsM_{\mathrm{s}} is the spontaneous magnetization, which is obtained from the extrapolation to zero field of M(H)M(H) curves, as shown in Fig. 4 (c) and (d). The magnetization of the sample can be related to the fraction of up domains nn_{\uparrow}, where M=(2n1)MsM=(2n_{\uparrow}-1)M_{\mathrm{s}}. From Happlied=NcMH_{\mathrm{applied}}=N_{c}M, we can determine nn_{\uparrow} as follows:

n=12MsHappliedNc+12n_{\uparrow}=\frac{1}{2M_{\mathrm{s}}}\frac{H_{\mathrm{applied}}}{N_{c}}+\frac{1}{2} (1)

The applied field, HappliedH_{\mathrm{applied}}, and the demagnetization factor, NcN_{c}, are temperature independent. Equation 1 gives us the result that, without pinning, the fraction of up domains should decrease as temperature decreases because MsM_{\mathrm{s}} increases.

Refer to caption
Figure 2: a) Temperature dependence of the fraction of up domains extracted from FCC (blue curve) and FCW (red curve) magnetization data at H=100H=100 Oe compared to our model based on demagnetization theory with Nc=0.608N_{c}=0.608 (black curve). b) Comparison between FCC and FCW magnetization at 100 Oe and δMs\delta M_{\mathrm{s}} with δ=0.31\delta=0.31. c) M2M^{2} versus H/MH/M at various temperatures. The intercept on the x-axis (orange line) reveals a demagnetization factor NcN_{\mathrm{c}} of 0.6080.608.

In Fig. 2(a) we compare the calculated fraction of magnetically aligned (up) domains with the value deduced from the experimental magnetization. The blue and red curves are the temperature dependence of the up-domain fraction extracted from the FCC and FCW magnetization respectively at 100100 Oe, with the following formula:

n=M2Ms+12.n_{\uparrow}=\frac{M}{2M_{\mathrm{s}}}+\frac{1}{2}. (2)

The black curve is the fraction of up-domain calculated from a demagnetization theory [30] in Eq. 1. We can see that domain wall pinning keeps the up-domain fraction above the theoretical value for all but the temperatures near TCT_{C}. While cooling down below TCT_{C}, the system encounters a domain-wall pinning effect where the blue curve begins to deviate from the black curve. This domain-wall pinning causes the magnetization to increase nearly proportionally to MsM_{s} as shown by the blue curve in Fig. 2(b). Below around 132132 K, the system enters a partial depinning region where the up-domain fraction decreases towards its theoretical no pinning value. The magnetization consequently decreases, which causes a local maximum near 132132 K. Upon further cooling, there is a competition between this anomalous depinning effect and the reduced thermal energy which tends to favor pinning. The reduced thermal energy eventually dominates around 120120 K and domain walls tend to pin again. Therefore, the magnetization tends to be proportional to MsM_{\mathrm{s}}, i.e. M=δMsM=\delta M_{\mathrm{s}} with δ=0.31\delta=0.31. The repinning is consistent with the fact that domain walls in Co3Sn2S2 rarely move below 119 K as reported by a Lorentz microscopy study [31]. Upon warming up, as shown by the red curve in Fig. 2, since domain walls are pinned, the domain fraction remains unchanged until around 155 K, where the increased thermal energy allows more up domains to form. Due to the presence of domain-wall pinning, the total magnetization of the sample when cooled down and warmed up follows two different pathways, resulting in the hysteresis between the FCC and FCW magnetization curves.

Refer to caption
Figure 3: a) Hysteresis loop with applied field parallel to the cc axis up to 100100 Oe. The arrow shows the coercive field HcH_{\mathrm{c}}. Nucleation field HnucH_{\mathrm{nuc}}^{\downarrow} is obtained by the intersection of red lines. b) HcH_{\mathrm{c}} and HnucH_{\mathrm{nuc}}^{\downarrow}, with maximum field up to 100100 Oe, as a function of temperature.
Refer to caption
Figure 4: Hysteresis loops with the applied field up to 77 T parallel to the cc axis at various temperatures. Only a selection of temperatures are shown. The hysteresis gradually changes from rectangle to bi-triangular shape. Spontaneous magnetization MsM_{\mathrm{s}} is obtained as shown in (c) and (d).

For evidence of domain-wall depinning effects, we measured hysteresis loops with the applied field only up to 100100 Oe at various temperatures. An example is shown in Fig. 3a. We note that domain structures depend on defects and the shape of a sample [32, 33]. For this reason, we use the same sample for all the bulk magnetization measurements presented in this article. The virgin curve at 130130 K remains flat at first, then increases with a constant slope, which is the hallmark of a domain-wall pinning-type virgin curve[30]. The magnetization of the field-increasing part of the hysteresis loop is smaller than that of the virgin curve, indicating that the domains encounter some pinning effect. Asymmetry between positive and negative coercive field emerges up to 160160 K. The asymmetry can be attributed to the Exchange Bias (EB), which has been observed by other groups [19, 34], or the remanent field during the measurement. At 130130 K, the negative coercive field is Hc=34.3(5)H_{\mathrm{c}}^{-}=-34.3(5) Oe, while the positive coercive field is Hc+=42.1(3)H_{\mathrm{c}}^{+}=42.1(3) Oe, and the actual coercive field, HcH_{\mathrm{c}}, is taken to be the average of the absolute value of the two. The nucleation field HnucH_{\mathrm{nuc}}^{\downarrow}, below which the down domains start to form, was obtained as shown in Fig. 3(a). If there were no domain-wall pinning, the M(H)M(H) curve should be a straight line. Therefore, the coercive field should always be zero and the nucleation field HnucH_{\mathrm{nuc}}^{\downarrow} should be the highest field of the hysteresis loop, i.e., 100100 Oe. However, Fig. 3(b) reveals a local maximum of HcH_{\mathrm{c}} and a local minimum of HnucH_{\mathrm{nuc}}^{\downarrow} at 132132 K, below which HcH_{\mathrm{c}} starts to decrease and HnucH_{\mathrm{nuc}}^{\downarrow} starts to increase, indicating a partial domain-wall depinning. The domain wall depinning around 130 K is consistent with the high mobility of domain walls observed near 130 K with a magnetic force microscope [22]. Further below 120120 K, domain-wall pinning strength increases again. HcH_{\mathrm{c}} consequently increases and HnucH_{\mathrm{nuc}}^{\downarrow} decreases. We note that the local maximum disappears at fields above 300300 Oe. As shown in Fig 5, 300300 Oe is the largest value of the nucleation field when this sample is fully magnetized. Therefore, down domains cannot form above 300300 Oe and the local maximum disappears. The actual value of the nucleation field depends on the shape of the sample and defects [32, 33]. Although the field above which the local maximum in the FCC curve disappears varies among different samples, the depinning around 130130 K is a general observation in this compound [18, 10, 11, 12, 13, 14, 15, 16]. While domain-wall pinning upon cooling is very common in many systems, depinning upon cooling is anomalous. The origin of the depinning, however, is still under investigation. It is possible that the depinning is enabled by the recently observed local structure distortion [20] which can reduce the local anisotropy energy and thus allows domains to flip, causing the depinning around 130130 K.

Since the sample is single domain at high fields, one can wonder if there is still evidence for the magnetic anomaly in the hysteresis loops up to 7 T. As shown in Fig. 4, the spontaneous magnetization, MsM_{\mathrm{s}}, is about 0.3 μB\mu_{B}/Co at 55 K [6, 25]. We also found that increasing temperature causes the shape of the hysteresis loop to change from a rectangular shape (5 K), to a triangular-tailed rectangle (100 K), to a bi-triangular shape (132 K), until the hysteresis disappears completely (160 K). The change from triangular-tailed-rectangle to rectangular shape of hysteresis loops can be explained with a domain-wall motion model. As shown in Fig. 5, the magnitude of the nucleation field increases as cooling down, until, around 8585 K, it exceeds the magnitude of the demagnetization field, causing the hysteresis loop changes from a triangular-tailed shape to a rectangular shape [35]. While there is no minimum in the nucleation field, we find that the nucleation field changes its sign around 130130 K, causing the hysteresis loop changes from a bi-triangular shape to a triangular-tailed shape.

Refer to caption
Figure 5: Temperature dependence of nucleation field and maximum demagnetization field. The red dashed line shows the 300300 Oe field.

3.2 Magnetic domain structures from 140 K to 120 K

Except for a recent MOKE study reported by Lee and coworkers  [21], domain structures of Co3Sn2S2(0001) single crystal samples are largely unknown. By measuring spatially resolved magnetization through the Kerr effect as a function of temperature, we can examine whether the magnetic anomaly is intrinsic to the physics in a single domain, or, as we suggested previously, is due to the presence of multiple domains and domain walls. Such detailed information is elusive from bulk measurements averaged over the entire sample.

To examine the domain structures near T=130T=130 K, we acquired a set of MOKE images at 140140 K, 136136 K, 132132 K, 128128 K, 124124 K, and 120120 K after the sample is field cooled in a magnetic field of H=+94{}_{\parallel}=+94 Oe. These images are shown in Fig. 6. Other than the disappearance of a small embedded red domain in the upper left corner (between 136136 K and 132132 K) and that of an embedded blue domain to its immediate right (between 128128 K and 124124 K), the domain structure remains essentially unchanged between 140140 K and 120120 K. Limited by the spatial resolution of the scanning microscope (0.85μ0.85\,\mum), from these MOKE images alone, we cannot rule out transitions that involve domain walls, such as those proposed in [21].

Refer to caption
Figure 6: Polar Kerr rotation images of Co3Sn2S2(0001) after FCC in HH_{\parallel} = + 9494 Oe (pointing out of paper along the c-axis) from room temperature acquired at (a) 140 K; (b) 136 K; (c) 132 K; (d) 128 K; (e) 124 K; (f) 120 K. Red regions are domains with magnetization pointing out of paper along the direction of the applied magnetic field; blue regions are domains with magnetization pointing in the opposite direction. The unit of Kerr rotation is μ\murad. The scale bar is 50 μ\mum. The image size is 240 μ\mum x 240 μ\mum.

We extract the fraction of up domains, down domains, and domain walls from these images as follows. First, we measure the average absolute Kerr rotation of up and down domains. Then, we count points above 75%75\% of the average value as up or down domains, and points below that value as domain walls. Fig. 7 shows the evolution of these fractions as a function of temperature. As the sample is cooled down, the fraction of up domains shrinks while that of down domains grows due to the demagnetization effect. The fact that the temperature variation of the domain fractions follows the same trend as the evolution expected from the case with no pinning confirms that the magnetic domain walls are indeed able to move between 120K and 140K.

Refer to caption
Figure 7: Evolution of the fraction of up and down domains and domain walls while cooling down extracted from MOKE images. The fraction of up domains shrinks while the fraction of down domains grows.

4 Conclusion

We find that the hysteresis between FCC and FCW magnetization in low fields is due to a domain-wall pinning effect and attribute the magnetic anomaly TAT_{\mathrm{A}} near 130 K in the FCC magnetization curve in Co3Sn2S2 to an anomalous domain wall depinning effect upon cooling. We also find that the nucleation field changes its sign around the anomaly temperature TA=130T_{\mathrm{A}}=130 K, causing the shape of the hysteresis loops to change from triangular-tailed rectangles into double-triangles. The MOKE images also reveal a magnetic domain walls depinning between 140140 K and 120120 K upon cooling. Our findings corroborate some other domain studies in this compound [31, 22]. However, the origin the domain-wall depinning still requires further investigation.

5 Acknowledgments

V.T., Z.S. and R.R.U. acknowledge support from the UC Lab Fees Research Program (LFR-20-653926) and UC Davis Startup funds. X.D.Z. acknowledges a Visiting Lecture Professorship from Fudan University in support of this work. We acknowledge support from the Physics Liquid Helium Laboratory fund.

6 Reference

References

  • [1] Wan X, Turner A M, Vishwanath A and Savrasov S Y 2011 Phys. Rev. B 83(20) 205101
  • [2] Burkov A A and Balents L 2011 Phys. Rev. Lett. 107(12) 127205
  • [3] Armitage N P, Mele E J and Vishwanath A 2018 Rev. Mod. Phys. 90(1) 015001
  • [4] Liu E, Sun Y, Kumar N, Muechler L, Sun A, Jiao L, Yang S Y, Liu D, Liang A, Xu Q et al. 2018 Nature physics 14 1125–1131
  • [5] Yan W, Zhang X, Shi Q, Yu X, Zhang Z, Wang Q, Li S and Lei H 2018 Solid State Communications 281 57–61 ISSN 0038-1098
  • [6] Wang Q, Xu Y, Lou R, Liu Z, Li M, Huang Y, Shen D, Weng H, Wang S and Lei H 2018 Nature communications 9 1–8
  • [7] Yang H, You W, Wang J, Huang J, Xi C, Xu X, Cao C, Tian M, Xu Z A, Dai J and Li Y 2020 Phys. Rev. Materials 4(2) 024202
  • [8] Nagpal V and Patnaik S 2020 J. Phys.: Condens. Matter 32(40) 405602
  • [9] Okamura Y, Minami S, Kato Y, Fujishiro Y, Kaneko Y, Ikeda J, Muramoto J, Kaneko R, Ueda K, Kocsis V et al. 2020 Nature communications 11 1–8
  • [10] Kassem M A, Tabata Y, Waki T and Nakamura H 2017 Phys. Rev. B 96(1) 014429
  • [11] Kassem M A, Tabata Y, Waki T and Nakamura H 2020 Journal of Physics: Condensed Matter 33 015801
  • [12] Zhang Q, Okamoto S, Samolyuk G D, Stone M B, Kolesnikov A I, Xue R, Yan J, McGuire M A, Mandrus D and Tennant D A 2021 Phys. Rev. Lett. 127(11) 117201
  • [13] Soh J R, Yi C, Zivkovic I, Qureshi N, Stunault A, Ouladdiaf B, Rodríguez-Velamazán J A, Shi Y, Rønnow H M and Boothroyd A T 2022 Phys. Rev. B 105(9) 094435
  • [14] Shin D H, Jun J H, Lee S E and Jung M H 2021 arXiv preprint arXiv:2105.03892
  • [15] Hu J, Kan X, Chen Z, Zheng G and Ma Y Journal of the American Ceramic Society
  • [16] Shen J, Zeng Q, Zhang S, Tong W, Ling L, Xi C, Wang Z, Liu E, Wang W, Wu G et al. 2019 Applied Physics Letters 115 212403
  • [17] Wu H, Sun P, Hsieh D, Chen H, Kakarla D C, Deng L, Chu C and Yang H 2020 Materials Today Physics 12 100189
  • [18] Guguchia Z, Verezhak J, Gawryluk D, Tsirkin S, Yin J X, Belopolski I, Zhou H, Simutis G, Zhang S S, Cochran T et al. 2020 Nature communications 11 1–9
  • [19] Lachman E, Murphy R A, Maksimovic N, Kealhofer R, Haley S, McDonald R D, Long J R and Analytis J G 2020 Nature communications 11 1–8
  • [20] Zhang Q, Zhang Y, Matsuda M, Garlea V O, Yan J, McGuire M A, Tennant D A and Okamoto S 2022 Journal of the American Chemical Society
  • [21] Lee C, Vir P, Manna K, Shekhar C, Moore J, Kastner M, Felser C and Orenstein J 2022 Nature Communications 13 1–6
  • [22] Howlader S, Ramachandran R, Singh Y, Sheet G et al. 2020 Journal of Physics: Condensed Matter 33 075801
  • [23] Lin X, Bud’ko S L and Canfield P C 2012 Philosophical Magazine 92 2436–2447
  • [24] Kassem M A, Tabata Y, Waki T and Nakamura H 2015 Journal of Crystal Growth 426 208–213
  • [25] Vaqueiro P and Sobany G G 2009 Solid state sciences 11 513–518
  • [26] Zhu X 2021 AIP Advances 11 085214
  • [27] Schnelle W, Leithe-Jasper A, Rosner H, Schappacher F M, Pöttgen R, Pielnhofer F and Weihrich R 2013 Phys. Rev. B 88(14) 144404
  • [28] Lamichhane T N, Taufour V, Thimmaiah S, Parker D S, Bud’ko S L and Canfield P C 2016 Journal of Magnetism and Magnetic Materials 401 525–531
  • [29] Arrott A 1957 Phys. Rev. 108(6) 1394–1396
  • [30] du Tremolet de Lacheisserie E, Gignoux D and Schlenker M (eds) 2004 Magnetism Fundamentals (Springer) ISBN 0387229671
  • [31] Sugawara A, Akashi T, Kassem M A, Tabata Y, Waki T and Nakamura H 2019 Phys. Rev. Materials 3(10) 104421
  • [32] Vial F, Joly F, Nevalainen E, Sagawa M, Hiraga K and Park K 2002 Journal of magnetism and magnetic materials 242 1329–1334
  • [33] Zeper W, Van Kesteren H, Jacobs B, Spruit J and Carcia P 1991 Journal of applied physics 70 2264–2271
  • [34] Noah A, Toric F, Feld T D, Zissman G, Gutfreund A, Tsruya D, Devidas T R, Alpern H, Vakahi A, Steinberg H, Huber M E, Analytis J G, Gazit S, Lachman E and Anahory Y 2022 Phys. Rev. B 105(14) 144423
  • [35] Tsymbal L, Kakazei G and Bazaliy Y B 2009 Physical Review B 79 092414