
Physical Review B 107, 024103 (2023)
Finite-dimensional signature of spinodal instability in an athermal hysteretic
transition
Anurag Banerjee1and Tapas Bar2, ∗
1Department of Physics, Ben-Gurion University of the Negev, Beer-Sheva 84105, Israel
2ICN2-Institut Catal`a de Nanoci`encia i Nanotecnologia (CERCA-BIST-CSIC),
Campus Universitat Aut`onoma de Barcelona, 08193 Bellaterra, Barcelona, Spain
(Dated: January 16, 2023)
We study the off-equilibrium critical phenomena across a hysteretic first-order transition in dis-
ordered athermal systems. The study focuses on the zero temperature random field Ising model
(ZTRFIM) above the critical disorder for spatial dimensions d= 2,3,and 4. We use Monte Carlo
simulations to show that disorder suppresses critical slowing down in phase ordering time for finite-
dimensional systems. The dynamic hysteresis scaling, the measure of explicit finite-time scaling, is
used to subsequently quantify the critical slowing down. The scaling exponents in all dimensions
increase with disorder strength and finally reach a stable value where the transformation is no longer
critical. The associated critical behavior in the mean-field limit is very different, where the exponent
values for various disorders in all dimensions are similar. The non-mean-field exponents asymptot-
ically approach the mean-field value (Υ ≈2/3) with increase in dimensions. The results suggest
that the critical features in the hysteretic metastable phase are controlled by inherent mean-field
spinodal instability that gets blurred by disorder in low-dimension athermal systems.
I. INTRODUCTION
The critical-like features in abrupt hysteretic tran-
sition have recently been observed in various materi-
als including transition metal oxide [1–4], metal alloys
[5, 6], martensitic transformation [7–9], functional mate-
rials [10], amorphous solids [11–13], microbiology, and
social, economic, climate, and other complex systems
[14, 15]. Such “surprising” [7] behavior is not normal
in terms of typical first-order phase transition formal-
ism. Some of such transitions have been explained in
terms of classical spinodal instability, a limiting point
of metastability (Fig. 1), where the system behaves
like a mean-field [3, 4, 16, 17]. The stability of the
metastable phase depends on the competition of dis-
order, thermal fluctuation, and activation barriers sep-
arating the two phases [10]. Any fluctuations, linked
with disorder or thermal, in the abrupt transition ini-
tiate nucleations before the extreme limit of metastabil-
ity [18]. In the long-range interacting system, thermal
fluctuations are suppressed [16, 19], and the metastable
phase of the system approaches the spinodal point af-
ter multiple cycling of the materials (training) across the
transition [3, 20]. The divergence of correlation length
and relaxation time scale (spinodal slowing down) signals
the instability in experiments [2, 4, 6]. The mean-field
spinodal universality in disorder material might be ex-
plained in terms of training-induced self-organized criti-
cality [21, 22]. However, the critical exponents often vary
widely from mean-field predictions [23] [see Table II] and
therefore remains unexplained. In general, the training
cannot tune the quenched disorders such as domain walls,
friction, defects due to an underlying heterogeneous sub-
∗tapas.bar@icn2.cat
strate, pinning defects, and kinetically arrested hetero-
geneity. Therefore, the correlation length of the system
would be bounded by the local disorder points, and het-
erogeneous nucleation sites start to emerge before ap-
proaching the spinodal [24–28]. As a result, a suppressed
spinodal slowing down associated with a mild finite-size
effect is expected to be observed [6, 13, 29] that may
explain such non-mean-field critical exponents. In this
article, we investigate spinodal instability using a ran-
dom field Ising model (RFIM) in the presence of quench
disorder and under athermal conditions. The athermal
(zero temperature) model mimics the fluctuationless ki-
netics associated with long-ranged potential, whereas the
short-ranged Ising model only deals with the interplay of
disorder and metastable barrier.
In RFIM, the critical signature in hysteretic tran-
sition has generally been observed in two distinct as-
pects: steady-state (slow-driven or quasistatic) and off-
equilibrium (highly-driven). The steady-state studies are
limited to the avalanche distribution and can explain the
disorder-induced critical transition near the critical dis-
order [30, 31]. Away from the critical point, the power-
law behavior of avalanche distribution is not adequately
understood [31]. One study attempts to explain such
phenomena at a low disorder regime in the context of
spinodal instability [13]. However, most of the hysteretic
transitions in materials are not single-step processes; in-
stead they show a broad transition accompanied by re-
turn point memory indicating the disorder in the system
is greater than the critical disorder [7, 30, 32]. There-
fore, further investigations are required above the crit-
ical point. On the other hand, the off-equilibrium as-
pect of critical phenomena such as dynamic hysteresis
scaling and phase ordering dynamics are comparatively
easy to measure in experiments. Not surprisingly, nu-
merous assessments have been reported for different ma-
terials [3, 6, 33–50]. In theory, several attempts have also
arXiv:2210.04057v2 [cond-mat.stat-mech] 12 Jan 2023