REVIEW 3 major objections 5 minor 37 references
Effects of random fields on the reversal of magnetisation of Ising ferromagnet
T0 review · 3 major / 5 minor · reviewed 2026-08-14 · deepseek-v4-flash
Pith's one-line read Adding quenched random fields to a uniform field reduces the magnetisation reversal time in the 2D Ising ferromagnet, and bimodal disorder does so most effectively.
desk verdict The main observation—quenched random fields shorten magnetization reversal time in 2D Ising—is solid and supported by multiple figures, but the cross-distribution comparison is undercut by an inconsistent Gaussian width convention. read the letter →
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
What carries the argument
The central object is the two-dimensional Ising ferromagnet with Hamiltonian $H = -J \sum_{\langle i,j\rangle} S_i S_j - \sum_i h_i S_i$, where $h_i = h_0 + h_r$ is a uniform field plus a quenched random field drawn from one of three zero-mean distributions (bimodal, uniform, or Gaussian). Dynamics is Metropolis single-spin-flip Monte Carlo, and the reversal time $\tau$ is defined as the first time the magnetisation becomes negative starting from a fully ordered state. The argument is carried by the Becker-Doring classical nucleation relation $\log \tau \sim 1/|h|^{d-1}$ in the single-droplet regime and $\sim 1/((d+1)|h|^{d-1})$ in the coalescence regime, together with the Avrami law for the metastable volume fraction. These supply the diagnostics: slopes in the $\log \tau$ versus $1/|h|$ plot identify the strong-field, coalescence, and nucleation regimes, and spin snapshots at reversal distinguish single-droplet from multidroplet growth.
What would settle it
Re-run the reversal-time comparison with all three distributions normalised to the same standard deviation (for example, set the Gaussian $\sigma$ equal to the standard deviation of the uniform and bimodal distributions at each $w$), and measure $\log \tau$ versus $1/|h_0|$ at $L=100$, $T=1.6$, $|h_0|=0.125$ for $w$ from $0.45$ to $0.5$; if the single-droplet regime still survives for uniform and Gaussian disorder at $w=0.5$, or if the Gaussian ranking reverses once its width is matched, the reported thresholds and ranking would not hold.
Extended reading notes
Core claim
In the two-dimensional Ising ferromagnet below the Curie temperature, adding a site-dependent but time-independent random field of zero mean to the applied uniform field systematically reduces the reversal time for all three distributions studied, and the reduction is largest for the bimodal distribution. Mean and most-probable reversal times both follow $\exp(-B w)$ in the disorder width $w$ and $\exp(A/T)$ in inverse temperature. The decay of the metastable volume fraction obeys the Avrami $t^3$ law, just as in the pure case. Becker-Doring analysis shows that the strong-field and coalescence regimes survive the addition of random fields, but the single-droplet nucleation regime is progressively destroyed as $w$ grows; for bimodal disorder it disappears between $w=0.3$ and $w=0.35$ at $|h_0|=0.125$ and $T=1.6$, while for uniform and Gaussian disorder it survives up to at least $w=0.45$.
Load-bearing premise
The comparison across disorder types assumes that $w$ means the same thing for all three distributions, but the Gaussian is described in the text as having standard deviation $\sigma=2w$ while the Fig. 5 caption says $\sigma=w/5$; if the text value is correct, the Gaussian disorder at a given $w$ is much broader than the bimodal and uniform disorder, so the effectiveness ranking and the reported loss of the nucleation regime are not controlled comparisons.
Editorial extensions
If this is right
- Adding quenched disorder of any of the three zero-mean forms shortens the reversal time, so disorder width can serve as a control knob for switching speed rather than only a source of noise.
- The exponential dependence $\tau \sim \exp(-B w)$ means modest increases in disorder width produce large reductions in reversal time; the fitted decay constants are roughly twice as large for bimodal disorder as for uniform or Gaussian disorder at $T=1.0$ and $|h_0|=0.5$.
- Bimodal disorder is the most effective of the three distributions studied, and the snapshots show it changes reversal morphology from single-droplet to multidroplet growth already by $w=0.35$ at $|h_0|=0.125$.
- For sufficiently large random-field width the single-droplet nucleation regime disappears, so the Becker-Doring single-droplet scaling no longer describes reversal; the coalescence description still applies.
- At higher temperatures the reversal-time distributions for the three disorder types converge, indicating that thermal fluctuations dominate over the details of the random-field distribution.
Reading between the lines
- A controlled comparison that matches the standard deviations of all three distributions would likely change the reported ranking and the critical width for loss of the nucleation regime; the manuscript is internally inconsistent about the Gaussian width (text says $\sigma=2w$, Fig. 5 caption says $\sigma=w/5$).
- If the equivalence between disorder width and temperature is taken seriously, engineered binary disorder could lower switching times without raising temperature, a design route for magnetic storage that the paper only gestures at.
- The disappearance of the single-droplet regime for bimodal disorder suggests a testable scaling: the critical width beyond which nucleation disappears should scale with $|h_0|$, since the mechanism is local fields opposing droplet growth; the paper does not derive this scaling.
- The finite-size check of the exponential $\tau(w)$ decay is performed only for bimodal disorder; repeating it for uniform and Gaussian disorder at $L=400$ would show whether the reported thresholds are size-dependent.
Signed reviews
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. This manuscript reports Monte Carlo simulations of magnetization reversal in a two-dimensional Ising ferromagnet driven by a uniform field plus quenched random fields drawn from bimodal, uniform, and Gaussian distributions. The central observations are that any zero-mean random field reduces the reversal time relative to the uniform-field-only case, that the bimodal distribution appears most effective among the three considered, that reversal times decay roughly exponentially with temperature and with disorder width, that the metastable volume fraction follows an Avrami-type law, and that for sufficiently large disorder width the single-droplet nucleation regime disappears at weak uniform field. The paper compares the measured slopes of log-reversal-time versus inverse field with Becker-Doring predictions for the strong-field, coalescence, and nucleation regimes.
Significance. If the quantitative claims are established, the paper would be a useful contribution to the study of metastable decay in disordered Ising systems: the qualitative effect of quenched random fields in shortening reversal time is robust and is supported by several independent figures (Figs. 7, 9, 13, 16, and 18), and the finite-size check in Fig. 12 is a commendable feature. The explicit fits in Figs. 10 and 11 provide concrete parameterizations of the temperature and width dependence. However, the cross-distribution comparisons are not currently controlled because the Gaussian width is defined inconsistently, and the regime classification and slope estimates lack objective criteria and error bars. These issues directly affect the quantitative ranking of the disorder distributions and the reported threshold for loss of the nucleation regime, so the manuscript needs major revision before the quantitative claims can be accepted.
major comments (3)
- [Section III, Eq. (12c) and Fig. 5 caption] The definition of the Gaussian disorder width is inconsistent. The text after Eq. (12c) states that sigma = 2w is 'considered throughout our work', while the Fig. 5 caption describes the plotted Gaussian random field as having sigma = w/5 and adds parenthetically that sigma = 2w was used in the work. The sentence following Eq. (12c) also says that the random-field values are distributed between -w and w, which is not true for an unbounded Gaussian. Because Figs. 7, 11, and 13-19 use w as the common control parameter, the comparison of bimodal, uniform, and Gaussian disorder is not controlled: at w = 0.45, the Gaussian standard deviation is either 0.9 or 0.09, leading to opposite interpretations of why the nucleation regime persists in Fig. 18. Please specify the exact sigma used in every simulation, and either replot the cross-distribution comparisons using a common measure of disorder strength (such as the standard deviation or a common interquartile range) or explicitly restrict the Gaussian ranking claims to the chosen convention.
- [Section IV, Figs. 13, 15, 16, and 18] The claim that the nucleation regime disappears for bimodal disorder between w = 0.30 and w = 0.35 but persists for uniform and Gaussian disorder at w = 0.45 is a central quantitative result, yet it is supported only by visual inspection of the log tau versus 1/|h| curves and by a small number of snapshots. No criterion is given for classifying a curve as possessing a nucleation regime, and the slopes in the purported SFR, CR, and NR regions are not fitted with error bars for the disordered cases. Please define an objective classifier (for example, a slope ratio close to 1/3 between the CR and NR regions, a curvature threshold, or a droplet-count measurement) and report the uncertainties in the fitted slopes and in the inferred critical width.
- [Section IV, Figs. 3, 13, 16, and 18] The comparison with Becker-Doring theory rests on the fitted slopes in the three regimes, but the pure-system slopes are quoted without uncertainties (Fig. 3: 0.75, 0.37, 1.08) and the disordered-system curves in Figs. 13, 16, and 18 are not fitted at all. Consequently, the statement that the coalescence-regime slope is approximately one third of the nucleation-regime slope, and the assertion that this relation survives in the presence of random fields, cannot be quantitatively assessed from the presented data. Please provide least-squares fits with confidence intervals for each regime and state how the regime boundaries were chosen.
minor comments (5)
- [Section III, after Eq. (12c)] The phrase 'distributed between -w to w' is inaccurate for the Gaussian distribution; please rephrase to describe the standard deviation or a quantile range.
- [Fig. 4 caption] The caption writes '1/h0 = -1.5' and similarly for the other regimes, but the axis labels use positive 1/|h0|; the signs are inconsistent and should be corrected to 1/|h0| = 1.5, 4.5, and 8.0.
- [Fig. 12 and Section IV] The decay exponent for L = 100 (b = 4.21) differs strongly from the values for L = 200-400 (b between 2.42 and 2.62), so the statement that L = 300 is 'free from any finite size effect' is stronger than the data support; please discuss the L = 100 deviation or restrict the claim to larger system sizes.
- [Section IV, near Fig. 6] There are typographical errors in the text, for example 'Acase for randomly distributed field thefter applying negative field'; the manuscript should be proofread for such issues.
- [Section IV, Fig. 8 and surrounding text] The Avrami-law statement refers to the 'third power of time' and the figure label is (t/tau)^3; the text should consistently write t^3 and specify the functional form ln(N1/N) proportional to -(t/tau)^3.
Circularity Check
No circularity: reversal-time claims are direct Monte Carlo measurements checked against independent Becker-Doring and Avrami benchmarks; the Gaussian width inconsistency is a non-circular control concern.
full rationale
The paper's central quantitative results are not derived from fitted inputs. Reversal times are measured directly from Monte Carlo metastable decays (Fig. 2), and the Becker-Doring check (Fig. 3 and Eqs. 7-9) compares measured slopes (about 0.37 versus the predicted 1/3, and 1.08 in the nucleation regime) against an independent classical nucleation formula. The temperature and width dependences in Figs. 10-11 are descriptive exponential fits to simulation data, not predictions generated from those fits, so no fitted parameter is renamed as a predicted result. The Avrami-law check in Fig. 8 is an external benchmark, not an input to the simulation. The only self-citation, Ref. [5] in the coalescence-regime discussion, merely cross-references a prior paper for a formula that is re-derived in Eqs. (8)-(9) within the present paper; it is not load-bearing. There is a genuine internal inconsistency in the Gaussian-width convention: Eq. (12c) states sigma = 2w, while the Fig. 5 caption states sigma = w/5. That ambiguity undermines the controlled comparison across distributions and the location of the nucleation-regime crossover, but it is a correctness or reproducibility concern rather than a circular derivation: no claimed prediction is defined in terms of the data it purports to explain, and no input quantity is equivalent by construction to the output claim. Accordingly, no specific circular step can be exhibited, and the appropriate circularity score is 0.
Assumptions & free parameters
free parameters (5)
- Arrhenius prefactor A and exponent B in tau ~ A exp(B/T) =
A ~ 0.4-0.7, B ~ 6.3-7.1 per distribution (Fig. 10)
- Prefactor C and exponent D in tau ~ C exp(-D w) =
C ~ 600-760, D ~ 0.94-2.53 per distribution (Fig. 11)
- Slopes of log tau versus 1/|h0| in SFR, CR, and NR =
0.75, 0.37, 1.08 (Fig. 3)
- Linear reference slope and intercept for Avrami plots =
-0.712, -0.1 and -0.524, -0.202 (Fig. 8)
- Gaussian standard deviation sigma relative to w =
sigma=2w in text, sigma=w/5 in Fig. 5 caption
assumptions (5)
- domain assumption The 2D Ising Hamiltonian with nearest-neighbor coupling J and local fields h_i = h0 + h_r describes the ferromagnet, and Metropolis dynamics samples the equilibrium distribution.
- domain assumption Becker-Doring classical nucleation theory: droplet free energy E_l = -2hl + sigma l^((d-1)/d), critical droplet size, nucleation rate I proportional to exp(-E_c/(kBT)), and reversal-time relations (7a) and (9a).
- domain assumption Avrami law: metastable volume fraction decays as exp(-const * t^(d+1)).
- ad hoc to paper For weak random fields, log tau versus 1/h0 remains linear with the theoretical slopes, and linearity is interpreted as validity of Becker-Doring theory in the disordered system.
- ad hoc to paper The width parameter w is a common measure across the three distributions.
Cite this review
Pith. "Pith review of Effects of random fields on the reversal of magnetisation of Ising ferromagnet." pith.science (2026). https://pith.science/paper/5S3AMSOV
@misc{pith2026190809268,
author = {Pith},
title = {Pith review of: Effects of random fields on the reversal of magnetisation of Ising ferromagnet},
year = {2026},
howpublished = {\url{https://pith.science/paper/5S3AMSOV}},
note = {Machine review of arXiv:1908.09268}
}
read the original abstract
We have studied the reversal time of the magnetisation in two dimensional Ising ferromagnet in the presence of externally applied uniform magnetic field using Monte Carlo simulation based on Metropolis single spin flip algorithm. Then we have investigated the change in reversal time due to the presence of quenched random field in addition to the uniform magnetic field. We report the results of statistical distribution of reversal times in the presence of three different types of the distributions (namely, uniform, bimodal and normal) of random fields and compared the results with those obtained for uniform field only. We have observed that the reversal time decreases due to the presence of any kind (of distribution) of the random fields. The metastable volume fraction is observed to follow the Avrami's law. Dependence of reversal times on temperature and different widths of the distributions of random fields are also reported. We have also checked whether the system obeyed Becker-Doring theory of classical nucleation in presence of additional random field and tried to investigate the range of the width of the distribution of random field. For larger width of the distribution of random field, the system fails to show the reversal via the nucleation of a single droplet (for small values of uniform field only). The possible reason is analysed.
Figures
Figures from the paper (16 more)
Reference graph
Works this paper leans on
- [1]
-
[2]
S. N. Piramanayagam and T. C. Chong, Development in Data S torage: Material Perspective, Wiley- IEEE Press, 2011
work page 2011
- [3]
-
[4]
P. A. Rikvold, H. Tomita, S. Miyashita and S. W. Sides, Phy s. Rev. E 49, 5080 (1994)
work page 1994
- [5]
- [6]
- [7]
- [8]
Show all 37 references
-
[9]
A. N. Kolmogorov, Bull. Acad. Sci. USSR Ser. Math. 3 (1937 ) 355
1937
-
[10]
A Johnson and P
W. A Johnson and P. A. Mehl, Trans. Am. Inst. Min. Metall. Eng. 135 (1939) 416
1939
-
[11]
Avrami, J
M. Avrami, J. Chem. Phys. 7 (1939) 1103; 8 (1940) 212; 9 (1 941) 177
1939
-
[12]
P. A. Rikvold, G. Brown, S. J. Mitchell and M. A. Novotny, Chapter 10, Nanostructured Magnetic Materials and Their Applications, Springer, Berlin, Heide lberg, 2002
2002
-
[13]
Hinzke, U
D. Hinzke, U. Nowak, phys. stat. sol. (a) 189 (2002) 475
2002
-
[14]
Beckmann, U
B. Beckmann, U. Nowak and K. D. Usadel, Phys. Rev. Lett. 9 1 (2003) 187201
2003
-
[15]
Brendel, G
K. Brendel, G. T. Barkema and H. V. Beijeren, Phys. Rev. E 71 (2005) 031601
2005
-
[16]
Machado , G
E. Machado , G. M. Buendia and P. A. Rikvold, Phys. Rev. E 7 1 (2005) 031603
2005
-
[17]
Acharyya, Physica Scripta, 82 (2010) 065703
M. Acharyya, Physica Scripta, 82 (2010) 065703
2010
-
[18]
W. R. Deskins, G. Brown, S. H. Thompson and P. A. Rikvold, Phys. Rev. B 84, 094431 (2011)
2011
-
[19]
Acharyya, Physica A 403, 94-99 (2014)
M. Acharyya, Physica A 403, 94-99 (2014)
2014
-
[20]
M.O.A Ellis and R.W Chantrell, Appl. Phys. Lett. 106, 16 2407 (2015)
2015
-
[21]
Dhar and M
A. Dhar and M. Acharyya, Commun. Theor. Phys., 66 (2016) 563
2016
-
[22]
Dutta, M
R. Dutta, M. Acharyya and A Dhar, Heliyon 4, e00892 (2018 )
2018
-
[23]
Imry and S
Y. Imry and S. -K. Ma, Phys. Rev. Lett. 35 (1975) 1399
1975
-
[24]
D. S. Fisher, J. Frolich and T. Spencer, J. Stat. Phys. 34 (1984) 863
1984
-
[25]
Aharony, Y
A. Aharony, Y. Imry and S.-K. Ma, Phys. Rev. Lett. 37 (197 6) 1364
-
[26]
N. G. Fytas, V. M. Mayor, M. Picco and N. Sourlas, J. Stat. Phys, 172 (2018) 665
2018
-
[27]
N. G. Fytas, V. M. Mayor, M. Picco, N. Sourlas, Phys. Rev. Lett. 116 (2016) 227201
2016
-
[28]
N. G. Fytas, V. M. Mayor, G. Parisi, M. Picco, N. Sourlas, Phys. Rev. Lett. 122 (2019) 240603. 11
2019
-
[29]
Kolesik, H
M. Kolesik, H. L. Richards, M. A. Novotny, P. A. Rikvold a nd P. A. Lindgard, J. Appl. Phys. 81 (1997) 5600
1997
-
[30]
Schiefele, I
B. Schiefele, I. S. Voivod, R. K. Bowles and P. H. Poole, P hys. Rev. E 87 (2013) 042407
2013
-
[31]
J. D. Gunton, M. Droz, Introduction to theory of Metasta ble and Unstable states, springer- verlag Berlin, 1983
1983
-
[32]
Vehkam¨ aki, Classical Nucleation Theory in Multico mponent Systems, Springer (2006)
H. Vehkam¨ aki, Classical Nucleation Theory in Multico mponent Systems, Springer (2006)
2006
-
[33]
Becker and W
R. Becker and W. D¨ oring, Ann. Phys. (Leipzig) 416 (1935 ) 719
1935
-
[34]
Binder and D
K. Binder and D. W. Heermann, Monte Carlo Simulation in S tatistical physics, Second edition, Springer-Verlag (1992), Berlin
1992
-
[35]
Metropolis, A
N. Metropolis, A. W. Rosenbluth, M. N. Rosenbluth, A. H. Teller, J. Chem. Phys. 21 (1953) 1087
1953
-
[36]
G. E. P. Box and M. E. Muller, The Annals of Mathematical S tatistics, 29 (1958) 610
1958
-
[37]
R. A. Ramos, P. A. Rikvold and M. A. Novotny, Phys. Rev. B 5 9 (1999) 9053. 12 Figure 1: Schematic variation of the free energy (for drople t formation) with size of the droplet -1 -0.8 -0.6 -0.4 -0.2 0 0.2 0.4 0.6 0.8 1 0 2000 4000 6000 8000 10000 L= 100 T= 1.6 |h0|= 0.14 <--...
1999
Reviewed August 14, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.