Low-temperature hysteresis broadening emerging from domain-wall creep dynamics in a two-phase competing system

Hysteretic behaviour accompanies any first-order phase transition, forming a basis for many applications. However, its quantitative understanding remains challenging, and even a qualitative understanding of pronounced hysteresis broadening at low temperature, which is often observed in magnetic-field-induced first-order phase transition materials, is unclear. Here, we show that such pronounced hysteresis broadening emerges if the phase-front velocity during the first-order phase transition exhibits an activated behaviour as a function of both temperature and magnetic field. This is demonstrated by using real-space magnetic imaging techniques, for the magnetic-field-induced first-order phase transition between antiferromagnetic and ferrimagnetic phases in (Fe$_{0.95}$Zn$_{0.05}$)$_{2}$Mo$_{3}$O$_{8}$. When combined with the Kolmogorov-Avrami-Ishibashi model, the observed activated temperature- and field-dependences of the growth velocity of the emerging antiferromagnetic domain quantitatively reproduce the pronounced hysteresis broadening. Furthermore, the same approach also reproduces the field-sweep-rate dependence of the transition field observed in the experiment. Our findings thus provide a quantitative and comprehensive understanding of pronounced hysteresis broadening from the microscopic perspective of domain growth.


Introduction
Functionalities carried by ferroic materials are closely related to their hysteretic behaviour that accompanies the first-order phase transition (FOT) associated with the order parameter reversal.Ferromagnets or ferroelectric materials with a high coercive magnetic or electric field facilitate a remnant order parameter; these materials have been applied to permanent magnets 1 , piezo elements 2 and recording media 3,4 .Alternatively, those materials with a low coercive field have superior sensitivity and switching ability 5 have been applied to motors, capacitors, transformers, and power supplies.In addition to these commercial products, the hysteretic behaviour accompanying a field-induced FOT between competing phases with different symmetries plays a key role in exotic functionalities, such as nonvolatile resistance control by a magnetic field 6,7 and magneto-or electrocaloric refrigeration 8 .From a microscopic viewpoint, hysteretic behaviour accompanying an FOT originates from the nonequilibrium evolution dynamics of the emerging phase, such as nucleation and growth; thus, the hysteretic behaviour tends to vary with the sweep rate of relevant parameters, also dictating the performance under high-speed operation 9 .
For a field-induced FOT in a single-crystalline sample, the transition profile is usually sharp, and accordingly, two characteristic transition fields, H * high and H * low, can be defined: one is higher than the field of the equilibrium FOT, Hc, and the other is lower.Thus, by plotting H * high and H * low at various temperatures, two hysteresis lines, H * high(T) and H * low(T), sandwiching the equilibrium FOT line, Hc(T), can be drawn in the phase diagram.However, in some materials exhibiting a magnetic-field-induced FOT, the difference between H * high(T) and H * low(T), i.e., the hysteresis width, starts to increase below a certain temperature and eventually shows a large hysteresis width of more than several Tesla at low temperatures.For example, the hysteresis width between ferrimagnetic up and down states in multiferroic LuFe2O4 reaches ~18 T at 4.2 K 10 .Although it is rare for ferroic materials to exhibit such high coercivity, there are many examples of field-induced FOTs between competing phases of different symmetries that exhibit pronounced hysteresis broadening toward low temperatures (Fig. 1a), for instance, the colossal magnetoresistive manganite 6 , Gd5Ge4 11 , LaFe12B6 12 , doped CeFe2 alloy 13 , and doped Mn2Sb 14,15 .Similar hysteresis broadening has also been observed for a structural FOT in martensitic materials 16,17 and aqueous solutions 18 , implying that common physics, which is not limited to magnetic materials, underlies these observations.However, different equations have been used for different systems to analyze hysteresis lines 10,16,17,[19][20][21] , and thus, the universal relationship between microscopic domain-wall motion and the broadening of hysteresis lines at low temperature remains unclear.Such hysteresis broadening is also potentially relevant to thermal-quenching-induced metastable states in quantum materials 22 .
For example, IrTe2 exhibits a metastable superconducting state when thermal quenching is applied 23 .This behaviour can be understood by assuming the hysteresis broadening accompanying the chemical-doping-induced FOT, although it is not possible to demonstrate directly the hypothetical broadening in this system by continuously changing the composition at low temperatures.
Among materials exhibiting low-temperature hysteresis broadening, we targeted (Fe0.95Zn0.05)2Mo3O8(Fig. 1b) because the hysteretic magnetic-field-induced FOT is accessible with a moderate magnetic field 24,25 , as shown in the phase diagram (Figs.1c,d); the FOT nature is confirmed by the fact that the extensive variables such as entropy and magnetization are discontinuous between the antiferromagnetic (AFM) and ferrimagnetic (FRI) phases and no critical phenomena are observed 32 .In (Fe1-yZny)2Mo3O8, the magnetism originates from two crystallographically inequivalent Fe 2+ (I) and Fe 2+ (II) sites (Fig. 1b), and Mo 4+ layers have no contribution to the magnetism 26 .The magnetic moments of Fe(II) sites are slightly greater than those of Fe(I) sites 24 .In Fe2Mo3O8, a metamagnetic transition occurs under magnetic fields along the c-axis 25,27 .In the low-field region, the two Fe(I) sites are antiferromagnetically ordered, and the same is true for the two Fe(II) sites, thus resulting in an antiferromagnetic (AFM) phase with no net macroscopic magnetization.When a sufficiently high magnetic field (~7 T) is applied, by contrast, the ferromagnetically aligned Fe(I) sites and the ferromagnetically aligned Fe(II) sites are antiferromagnetically ordered, thus resulting in a ferrimagnetic (FRI) phase with macroscopic magnetization.Doped Zn ions selectively occupy the Fe(I)-sites, which decreases the critical field of the metamagnetic transition 28 and causes the broadening of the hysteresis region at low temperatures 25 .For compositions above y=0.05, the FRI phase is so stable that the FRI to AFM transition can be observed only in a very narrow temperature-field region; on the other hand, for compositions below y=0.05, the AFM phase is so stable that the magnetic-field-induced FOT is difficult to observe.For these reasons, the present study focused on the composition of y=0.05, in which the transition field is well accessible and thus experiments are feasible.
In this Article, we show that the pronounced hysteresis broadening originates from an activated behaviour of domain-wall dynamics as a function of both temperature and magnetic field, i.e., the creep motion of a domain wall.This is evidenced by real-space magnetic imaging experiments and numerical simulations.Furthermore, we also find that a sweep-rate dependent magnetic hysteresis loop, which is important for understanding the performance under highspeed operation, is also well reproduced within our numerical simulations.The present study unveils a strong correlation between microscopic domain-wall motion and the macroscopic transition accompanied by low-temperature hysteresis broadening.

Results and discussion
Observation of an antiferromagnetic-domain growth.
To reveal the microscopic origin of the hysteresis broadening in (Fe0.95Zn0.05)2Mo3O8from the perspective of the nonequilibrium phase-evolution dynamics, we performed magnetic force microscopy (MFM) and magneto-optical Kerr effect (MOKE) imaging.In particular, we focused on AFM phase evolution in the matrix of the FRI phase, close to the lower-field boundary of the hysteresis region.To this end, a single-phase FRI state was first prepared by applying a magnetic field of 7 T at high temperature (~25 K), followed by field-cooling to a target temperature under 7 T.Then, the magnetic field was decreased to a target field, and the time evolution of the emergent AFM domain was measured while holding the temperature and magnetic field constant.Figures 1e-g show the time evolution of the MFM images at 15 K and 1.8 T, noting that this comes from a magnetic signal (Supplementary Note 1).In the initial state, a small disk-shaped AFM domain was observed in the matrix of the FRI phase, and it gradually and continuously expanded over ~10 4 s (Figs.1f,g).A similar AFM domain was also observed at a different position.While MFM has a high spatial resolution of < 50 nm, image acquisition takes 10-20 minutes, and thus, rapid domain growth cannot be tracked.To track the faster motion of the AFM-domain growth, we utilized MOKE imaging, which has a lower spatial resolution, ~1 m, but a faster image acquisition time, i.e., 1 s.In MOKE imaging, the AFM domains were observed as similar multiple circular spots with dark boundaries in a much wider field of view (Fig. 1h).Note that images shown in panels Figs.1h and 1i-k are divided by the raw MOKE image at 0 s to remove the influence of the position-dependent intensity of the incident light.Figures 1i-k display the AFM domain growth at 15 K and 0.9 T, demonstrating that compared with the data at a higher magnetic field (Figs.1e-g), the growth speed of the AFM domains is pronouncedly faster.
We found that the radius of the circular AFM domain changes almost linearly with respect to elapsed time, as shown in Figs.1l and 1m, and from this behaviour, a constant growth velocity, v, of the AFM domain was defined as the slope of this plot; the time evolution of the circular domain appears monotonic, suggesting that the spatial distribution of impurities can be regarded as homogenous as far as the present domain-wall dynamics are concerned.The same analysis was performed for various temperatures, T, and magnetic fields, H, and thus we were successful in deriving v(H, T) for T ≤ 20 K: Above this temperature, the image contrast significantly decreased, although the reason is not clear; therefore, the analysis was not possible.
The results are summarized as a contour plot of the growth velocity in Fig. 1n (see also Supplementary Note 2 for the reproducibility at different positions), and to gain insight into the relationship between v(H, T) and the hysteresis line H * low(T), the contour plot is superimposed on the phase diagram.Here, two important points can be highlighted.On the one hand, the curvature of the contour line is similar to that of the H * low(T) line, suggesting that the profile of the hysteresis line is closely related to the microscopic growth process of the AFM domain.On the other hand, although the growth speed near the H * low(T) line is expected to have a particularly important impact on the macroscopic phase evolution under the continuous fieldsweeping experiments, the information was not obtained from the above analysis because the growth speed is far beyond the detection limit, ~1 m s -1 .

Creep motion of a domain wall between the antiferromagnetic and ferrimagnetic phases.
To obtain v(H, T) near the H * low(T) line, the microscopic mechanism should be identified together with the corresponding formula.To this end, below we consider the application of the so-called modified Merz's law that describes the creep motion of an elastic interface, such as a domain wall (DW), in a disordered medium 29,30,31 .The characteristic of the modified Merz's law is that the velocity of a moving DW obeys activated behaviour with respect to both temperature, T, and the driving force exerted on the DW, F. The equation is given as where  is a so-called dynamical exponent and A and  represent the limit of high speed and temperature-force-composite activation barrier, respectively.When applying this equation to (Fe1-yZny)2Mo3O8, the driving force should be measured with reference to Hc(T) (for more details, see Supplementary Note 3), which has recently been determined by making full use of the thermodynamic relations 32 as indicated by white circles in Fig. 1d.Thus, in the present case, the modified Merz's law is given as where H(T) = |H(T) -Hc(T)|.Equation ( 1) describes the average velocity of the creep motion, rather than the instantaneous velocity that may be affected by the details of the spatial distribution of impurities.
Figure 2a displays the growth velocities at various temperatures, plotted against 1/(H)  .We find that the observed velocities are well reproduced by equation ( 1) (solid lines, Fig. 2a), indicating that the observed DW motion should be dictated by the creep motion.We note, however, that it is not so obvious to what extent the creep motion describes the actual DW dynamics outside the examined parameter region; in fact, previous theoretical studies predict that the DW dynamics change from the creep to flow regimes as the driving force increases 30 .We will discuss this issue later.The optimized adjustable parameters are A=3806.5 [m s -1 ], which is quite close to the expected velocity of the acoustic phonon in Fe2Mo3O8, 2.7-6.610 3 m s -1 calculated from the elastic constants 33 , and μ=0.4.The similarity between the value of A and the phonon velocity may imply that structural change accompanies the metamagnetic transition and thus affects the dynamics of DW.The value of  is nearly constant (~500) within the measured temperature range (Fig. 2b), and the residual weak temperature dependence may be attributed to that of the order parameter.In ferroic materials, it is known that the value of μ depends on both the dimensionality of the DW and the type of disorder, which is generally classified into random-bond or random-field types 30,34,35 .However, for the case of DWs separating distinct symmetry phases, such as AFM and FRI, such classification of disorder is not straightforward (for details, see Supplementary Note 4).Thus, although μ = 0.4 appears to be close to μ = 0.5-0.6,which is a value reported for two-dimensional DWs under the influence of random-bond-type disorder 31,36 , its implications are not clear.
Origin of the pronounced hysteresis broadening.
Nevertheless, having established that equation ( 1) reproduces the observed (, ) well, we can supplement the data near the H * low(T) line by the extrapolation according to equation ( 1), enabling us to quantitatively investigate the relationship between the creep-like growth velocity and the pronounced hysteresis broadening.To this end, we use the Kolmogorov-Avrami-Ishibashi (KAI) model 37,38 , which can calculate the volume fraction of a growing domain for a given DW velocity while avoiding the complexity caused by domain coalescence (Fig. 3a).
Given that in the present observation, the nucleation sites are always at the same positions, we applied the formalism of the KAI model with the so-called one-step nucleation mechanism (for details, see Supplementary Note 5): where c(t), N, and S(t) are the time-dependent volume fraction of the AFM domain, the density of nucleation sites, and the time (t)-dependent extended area caused by a single nucleation site, respectively.When calculating S(t), the DW velocity should be integrated with respect to t: In the present study, because a phase transition under continuous field-sweeping experiments is considered, the DW velocity is time dependent via the time-dependent driving force, H(t).
The calculation of the phase evolution was performed by using the experimentally obtained parameter set, A=3806.5, μ=0.4 and =500.We define the transition field Hc(T) as the field at which the phase evolution is half complete.The experimental transition profile at 20 K and the corresponding calculated results are shown in Figs.3b and 3c, respectively.The details of the transition profile, such as the value of Hc(T) and the sharpness of the transition, depend on N and n (Fig. 3c).However, we find that the shape of the hysteresis line is not sensitive to N and n (Fig. 3d and Supplementary Figs.4a,b).Therefore, upon comparing the hysteresis line profiles of the calculation and experiments, the normalized temperature dependence of Hc(T)/Hc(12 K) is the quantity that should be considered.The comparison is displayed in Fig. 3d, showing good agreement with the experiment.Furthermore, the profile of the hysteresis broadening is less pronounced for a higher μ value, such as μ = 1 (Supplementary Note 5), which is a value reported for DWs in ferroic materials with randomfield type disorder 31,39,40 .These observations indicate that the peculiar T and H dependences of the DW creep velocity, underpinned by a low μ value, are at the heart of the pronounced hysteresis broadening.

Sweep-rate dependence.
Finally, to further corroborate the relationship between the creep dynamics of the DWs and the pronounced hysteresis broadening, it would be important to test another aspect of the hysteresis line, that is, the sweep-rate dependence.If the H * low(T) line is dictated by the creep dynamics, the H * low(T) line should depend on the sweep rate as a natural consequence of the kinetic aspect.
Figure 4a shows the sweep rate dependences of the isothermal magnetization curves at 20 K, and it can be seen that the FRI-AFM transition field (white circles) appreciably decreases as the sweep rate increases.This feature is conspicuous especially below 30 K (Fig. 4b), highlighting the kinetic aspect of the FOT process.Using the KAI model, the field-induced transition at various field-sweeping rates at 20 K can be calculated and compared with experimental results (Figs. 4c and 4d).As in the case of the temperature dependence of the hysteresis line, we find that the normalized sweep-rate dependence [Hc(r)-Hc(rmin)]/[Hc(rmax)-Hc(rmin)] depends on N and n only weakly.We therefore compare this value between the experiments and calculations at selected temperatures, as shown in Fig. 4e (see also Supplementary Figs.4d-f), confirming that the sweep-rate dependence of the transition field is successfully captured by the calculation at a quantitative level.

Conclusion
The present analysis relies on extrapolating the (, ) data towards Hlow within the framework of the modified Merz's law (equation ( 1)), and thus, the validity of the model near Hlow may not be necessarily clear.However, we have demonstrated that the two main characteristics of the hysteresis line at low temperatures, that is, the pronounced broadening and the pronounced sweep-rate dependence are both quantitatively explained by considering the creep motion of the DW obeying an activate form with respect to temperature and magnetic field.These observations indicate that the modified Merz's law still captures the DW dynamics near Hlow.Given that not all FOT materials exhibit pronounced hysteresis broadening, the DW creep dynamics, including the μ value, are likely to vary widely from material to material.So far, DW dynamics has been investigated exclusively in ferroic systems because of their importance in applications.Toward exotic applications based on an FOT beyond ferroic systems, our findings suggest that control of μ in a phase-competing system is likely a key ingredient.

Methods
Crystal growth.Single crystals of (Fe0.95Zn0.05)2Mo3O8were grown by chemical vapor transport reaction method following the literature 24,25 , 41 .Powders of MoO2, Fe, Fe2O3, and ZnO were mixed and ground in the stoichiometric molecular ratio.The powder was sealed in the inner quartz tube (diameter 15 mm, length 130 mm) together with 100 mg of TeCl4 as the transport agent.This ampoule was double-sealed with the outer quartz tube (diameter 20 mm, length 200 mm), and placed in a three-zone furnace.After 24 hours of back transport (845 ℃ on the precursor-loaded side; 980 ℃ on the growth side), the temperature gradient was reverted to be 980 °C for the load side and 845 °C for the growth side, and maintained for 11 days.
Magnetization measurement.Magnetization was measured by reciprocating sample option (RSO) and vibrating sample magnetometer (VSM) modes in a commercial superconducting quantum interference device magnetometer (MPMS-XL and MPMS-3; Quantum Design).Magnetic sweeping rates from 2  10 -4 T s -1 to 2  10 -2 T s -1 were controlled in MPMS-3 and the VSM measurement was performed in the continuous sampling mode.
Magnetic force microscopy measurement.Frequency-modulated MFM was performed under high vacuum conditions with a commercially available scanning probe microscope (attocube AFM/MFM I).We used the MFMR tip (supplied by NANOSENSORS) with a tip radius of ~50 nm.The typical sample dimension was approximately 3.0×2.0×1.0 mm 3 with a mass of 31.95 mg.We prepared a sample with a naturally grown surface that is normal to the c-axis.Gold with a thickness of 50-100 nm was deposited on the surface to avoid the charge up problem during the measurement.The MFM measurement was performed in noncontact mode with a lift height of 100-200 nm.The amplitude of the cantilever oscillation was approximately 10-20 nm.The resonant frequency was ≈63.6 kHz, and the Q-factor was ≈3.5×10 4 under the measurement conditions.

Magneto-optical Kerr effect (MOKE) measurement. MOKE measurement setup
follows the literature 42 .A commercial polarizing microscope (BXFM, Olympus) was fixed to the aluminum frame self-constructed above the physical property measurement system (PPMS; PPMS-14T, Quantum Design).An infinity-corrected objective lens (PLN4XP, Olympus) attached to a fiberglass reinforced plastic tube (diameter 16 mm, length 1m) was inserted into the sample chamber of PPMS.Light from a LED source (625 nm, M624L4-C1, Thorlabs) passed through a polarizer and coaxially illuminated onto the sample.Reflected light through the objective lens, an analyzer, and a television (TV) lens (magnification of 5) was captured by a CMOS camera (PCO).The polarization is shifted about 5 degrees from the orthogonal configuration 42 .The temperature and magnetic field were controlled by the PPMS control system.The sample on the sapphire plate was placed on the PPMS sample puck (P101/3A) with varnish and installed into the sample chamber of PPMS.The measurement was performed with the naturally grown surface of the same sample used in the MFM measurement.wall is on nonmagnetic impurity sites.Thus, doped Zn ions should work as RB-type impurities in the present system.Second, given that a delicate energy balance between the AFM and FRI phases results from the magnitudes of the interlayer Fe(I)-Fe(II) and interlayer Fe(I)-Fe(I) interactions in this system, doped nonmagnetic Zn ions are expected to stabilize the FRI phase by weakening Fe(I)-Fe(I) magnetic interactions.In fact, as the amount of Zn doping increases, the FRI phase tends to be realized at relatively low magnetic fields, indicating that Zn doping plays a role in changing the energy balance between the two competing phases.These two mechanisms are expected to work simultaneously, and therefore, a doped Zn ion can be considered to have characteristics of both RF-type and RB-type impurities.
For a system that exhibits competition between phases with different symmetries, any impurity can have an asymmetric impact on the two competing phases, and it can also interact with the domain wall when the domain wall is on the impurity.It thus appears that any impurity can invariably have the characteristics of RB-and RF-type impurities, and it is therefore nontrivial what value of  should be obtained by the introduction of impurities in a phasecompeting system, such as (Fe1-yZny)2Mo3O8.

Fig. 1 |
Fig. 1 | Low-temperature hysteresis broadening and growth of the AFM domain in (Fe1-yZny)2Mo3O8.a, Archetypal phase diagram of a field-induced first-order phase transition (FOT) between two competing phases (labelled by A and B) with different symmetries.The solid lines in a represent magnetic fields at which the phase transition is observed during magnetic field sweeps at a given rate.Such hysteresis broadening is not always but is frequently observed in a field-induced FOT between competing phases.b The crystal structure of (Fe1-yZny)2Mo3O8.Fe 2+ ions at the Fe(I)-sites (brown) and Fe(II)-sites (blue) are surrounded by oxygen tetrahedra and octahedra, respectively.Mo 4+ ions (gray) form the nonmagnetic spin-trimer.c Isothermal magnetization curve of (Fe0.95Zn0.05)2Mo3O8at 20 K. Blue (red) open circles indicate the transition fields from the FRI to AFM (AFM to FRI) phases, where AFM and FRI represent antiferromagnetic and ferrimagnetic, respectively.These transition fields are defined as the inflection points of M-H curves.d Temperature-magnetic-field phase diagram of (Fe0.95Zn0.05)2Mo3O8.The transitions from the AFM to FRI phases and from the FRI to the AFM phases are indicated by red and blue open circles, respectively.Black open circles with error bars represent the equilibrium FOT phase boundary, taken from ref. 32 .Orange and blue arrows in the schematic represent spins at Fe(I) and Fe(II) sites, respectively.e-g Magnetic force microscopy (MFM) images of a timeevolving AFM domain at 15 K and 1.8 T at (e) 0 s, (f) 4.810 3 s and (g) 8.810 3 s.h Magneto-optical-Kerr-effect-microscopy (MOKE) image of a wider area at (15 K, 0.9 T). i-k MOKE images of a time-evolving AFM domain at (15 K, 0.9 T) indicated by the dotted rectangle in panel h at (i) 50 s, (j) 100 s and (k) 150 s.Images shown in panels h-k are divided by the raw MOKE image at 0 s to remove the influence of the positiondependent intensity of the incident light.l,m Time dependences of a radius at (l) (15 K, 1.8 T) and (m) (15 K, 0.9T).n Contour map of the growth velocity, superimposed on the phase diagram.The blue solid line with blue open circles indicates the transitions from the FRI to AFM phases on a field-decreasing process.The black dashed lines on the contour map represent contour lines drawn at 1 and 3 series from 1.010 -10 to 1.010 -6 in logarithmic representation.Gray cross and plus marks represent the points measured with MFM and MOKE, respectively.

Fig. 2 |
Fig. 2 | The fitting of growth velocities of an antiferromagnetic domain with Merz's law. a Growth velocities plotted against the inverse of the driving magnetic field.The data indicated by filled circles and open triangles are obtained from the MFM and MOKE experiments, respectively.The solid lines represent the fitting with

Fig. 3 |
Fig. 3 | Calculation of the low-temperature hysteresis broadening based on the Kolmogorov-Avrami-Ishibashi (KAI) model combined with modified Merz's law. a Schematic illustration of the KAI model for a two-dimensional case.Gray regions represent nuclei of a new phase, which are isotropically growing.b,c The experimental (b) and corresponding calculated (c) transition profiles at 20 K.The vertical axis in panel b represents the volume fraction that is calculated as [M(7T)-M(H)]/[M(7T)-M(0T)], where M(H) is magnetization.The horizontal axis is H(T) = |H(T) -Hc(T)|.The calculation was performed following equation (3) at selected values of (n,N), where n and N represent the system dimension and the density of the nucleation sites, respectively.d Comparison of the transition fields from the FRI to AFM phases, obtained from the calculation and experiments.The transition fields are normalized by the transition field at 12 K.The experimental and simulated data are represented by the filled circles and curves, respectively.

Fig. 4 |
Fig. 4 | Sweep-rate-dependent behaviour of the hysteresis lines.a Isothermal magnetization curves at 20 K measured at various sweep rates, obtained in the experiments.b The transition lines determined by various sweep rates.Open circles represent the sweep-rate dependent transition fields from FRI (AFM) to AFM (FRI) phases, determined by isothermal magnetization-curve measurements.c,d The experimental (c) and corresponding calculated (d) sweep-rate-dependent transition profiles at 20 K.The definitions of the vertical and horizontal axes are the same as in Figs.3b and 3c.e Comparison of the sweep-rate dependent transition fields from the FRI to the AFM phases at selected temperatures, obtained from the calculation and the experiments.The transition fields are normalized as [Hc(r)-Hc(rmin)]/[Hc(rmax)-Hc(rmin)] at the lowest and highest sweep rates, i.e., rmin = 210 -4 T s -1 and rmax = 210 - 2 T s -1 .The experimental and simulated data are represented by the filled circles and curves, respectively.Each curve is shifted by 0.2 for visibility.