REVIEW 3 major objections 5 minor 29 references
Magnetic polarisability of octet baryons near the physical quark-mass point
T0 review · 3 major / 5 minor · reviewed 2026-08-11 · deepseek-v4-flash
Pith's one-line read The paper reports the first lattice QCD calculation of octet-baryon magnetic polarisabilities at a pion mass of 156 MeV, made possible by a new algorithm that identifies and removes exceptional configurations caused by electro-quenching…
desk verdict First near-physical-mass background-field lattice QCD calculation of octet-baryon magnetic polarisabilities, with a genuinely new exceptional-configuration cleaning algorithm, but the proton result rests on a subjective stopping rule and should be treated as indicative rather than final. 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 ratio $R(B,t) = G_{\uparrow\uparrow}(B,t)\,G_{\uparrow\downarrow}(B,t)/G(0,t)^2$ of spin-aligned and anti-aligned two-point correlation functions, whose effective-energy difference strips away the baryon mass and magnetic-moment terms and leaves the magnetic polarisability energy shift $\delta E_\beta(B,t) = |q_B eB|/(2m) - 2\pi\beta B^2 + O(B^3)$. The new machinery is the exceptional-configuration identification algorithm: it slides a window of $M=20$ correlation functions in configuration time, computes the uncertainty $\Delta E_i(t)$ of the sub-ensemble effective energy at early time slices, and flags as exceptional any configuration whose entry causes $\Delta E_i$ to jump above a threshold $(\Delta E)_{\max} = \overline{\Delta E} + n_\sigma \sigma(\Delta E)$ for at least $M$ consecutive indices. The threshold $n_\sigma$ is lowered from 20 to 2 in steps of 0.1--0.25 over several randomised passes, and removal stops when the $\delta E_\beta$ plateau visibly resolves.
What would settle it
Generate new dynamical ensembles with the background magnetic field coupled to the sea-quark charges (removing electro-quenching) and recompute the proton and neutron polarisabilities at the same lattice spacing, volume, and pion mass; if the plateau appears without any removal or the extracted values move by more than the quoted errors, the visual stopping criterion is biasing the central values.
Extended reading notes
Core claim
The central discovery is that the exceptional configuration problem in background-field lattice QCD at near-physical quark masses originates from electro-quenching: the dynamical gauge-field generation neglects the electric charges of the sea quarks, so the fermion determinant does not suppress gauge fields that become improbable once the background magnetic field is introduced. The paper establishes an algorithm that detects these configurations as step-function jumps in the uncertainty of the effective energy of small sub-ensembles of correlation functions, removes them, and repeats with a decreasing threshold until the energy-shift plateau resolves. After removal, the magnetic polarisability energy shift of the outer octet baryons ($p$, $n$, $\Sigma^+$, $\Sigma^-$, $\Xi^0$, $\Xi^-$) exhibits stable plateaus, and the extracted polarisabilities are in accord with chiral perturbation theory predictions at the $1\sigma$ level.
Load-bearing premise
The stopping criterion for the removal algorithm is chosen by visually inspecting animations of the effective-energy shift and stopping as soon as the plateau resolves; the paper notes that further removal continues to change the central value, with the magnetic polarisability tending to increase.
Editorial extensions
If this is right
- At mπ = 156 MeV, octet-baryon magnetic polarisabilities can be extracted with good precision once exceptional configurations are removed, with values agreeing with chiral perturbation theory at the 1σ level.
- The proton requires removing about 1.4% of the sampled configurations, showing that the exceptional-configuration problem is finite and manageable under the algorithm's conservative threshold sweep.
- A single strange quark sharply reduces susceptibility to exceptional configurations, while a doubly represented up quark increases it, consistent with the up quark's larger electric charge.
- Including the new light-quark point has minimal impact on the physical-point chiral extrapolations from heavier masses, so the previous predictions remain stable.
- Higher statistics alone did not cure the proton signal; configuration removal rather than sample size was the decisive step.
Reading between the lines
- Inference: The algorithm's sensitivity to the up-quark charge suggests that the field-strength quantum could be re-optimised for observables dominated by the up quark, where a finer field grid might reduce the exceptional-configuration rate.
- Inference: Because the stopping criterion is visual, the reported central values are best read as provisional lower bounds on the polarisability if further removal continues to increase the value; a blinded or automated stopping rule would settle the direction of the residual systematic.
- Inference: The same outlier-identification machinery could be applied to other electro-quenched background-field observables, such as magnetic moments or form factors at light quark masses, where similar outliers may currently inflate uncertainties without being recognised.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. This paper presents the first lattice QCD calculation of the magnetic polarisabilities of the outer octet baryons at a near-physical pion mass, m_pi = 156 MeV, using the background field method on electro-quenched PACS-CS gauge configurations. The authors identify a new 'exceptional configuration' problem caused by electro-quenching, in which improbable gauge fields lead to large additive mass renormalisations and outlier correlation functions. They develop an algorithm to identify and remove such configurations, and after removal they extract magnetic polarisabilities for the proton, neutron, Sigma^+, Sigma^-, Xi^0, and Xi^- using a weighted fit to the magnetic-field-dependent energy shift. The results are compared with chiral perturbation theory through a finite-volume and electro-quenching-corrected chiral extrapolation based on the authors' previous heavier-mass results. The paper finds agreement typically at the 1-sigma level, while noting that the new light-quark points have minimal influence on the extrapolation.
Significance. If the extraction is robust, this is the first calculation of octet baryon magnetic polarisabilities at a near-physical pion mass, and it identifies an important systematic effect of electro-quenching that will affect future background-field lattice QCD calculations. The paper is transparent about the removal procedure and provides supplementary animations illustrating the evolution of the effective energy as configurations are removed. However, the quantitative claim of agreement with chiral perturbation theory is currently weakened by the subjective stopping criterion for configuration removal and by the absence of any systematic uncertainty assigned to that choice. The work is nevertheless significant as a first step toward controlling this exceptional-configuration problem.
major comments (3)
- [Sec. IV.D and V.A] The stopping criterion for exceptional configuration removal is chosen by visual inspection of the supplementary animations to identify a 'step change' in the effective-energy plateau. The paper concedes in Sec. IV.D that 'it is not easy to define a consistent stopping criterion,' and in Sec. V.A that the magnetic polarisability is 'sensitive to continued configuration removal, tending to increase to larger values.' No systematic error is assigned for this choice. The proton result exists only after removing 773 of approximately 50,000 configurations, and the hyperon values shift by more than their statistical errors (e.g., Sigma+ from 1.7(10) to 2.17(51); Xi0 from 2.14(32) to 2.32(20)). The claimed 1-sigma agreement with chiral perturbation theory is therefore not robust; the authors should estimate the systematic uncertainty by scanning over removal thresholds or n_sigma,min values and quoting the spread, or provide an objective a priori stopping rule.
- [Sec. II and Table II] The fit windows are very short for the noisiest channels. For the neutron, the fit range [tmin,tmax] = [20,22] admits a single window of length three, which cannot establish a plateau and provides no lever arm against excited-state contamination. The AIC weighting in Eq. (9) cannot suppress excited states when only one window is eligible. The authors should demonstrate that the extracted value is stable under small changes to tmax (for example by fitting at tmax = 23 or 24 when signal permits), or include a systematic error that reflects the short-window sensitivity.
- [Sec. VI] The chiral extrapolation and the correction scheme are taken from the authors' own Ref. [5], and the new light-quark-mass points have 'minimal impact on the result of the extrapolation to the physical point' (as stated in Sec. VI). Consequently, the agreement with chiral perturbation theory is a consistency check with a fit dominated by the four heavier masses, not a strong quantitative test of the new results. The paper should state this limitation clearly when claiming 1-sigma accord.
minor comments (5)
- [Abstract] The sentence 'the dynamical-fermion gauge-field generation algorithm the electric charges of the quarks' appears to be missing a verb; please revise.
- [Sec. VII] 'choral effective field theory' should be 'chiral effective field theory'.
- [Sec. III.G] '250 K unique correlation functions' is ambiguous; please write '250,000' or define K explicitly.
- [Captions of Figs. 4-6 and Supplemental] The red bar indicating the number of removed configurations uses a nonlinear scale; please state explicitly in the captions that the scale is logarithmic.
- [Supplemental Material, Fig. S-6] The caption for Fig. S-6 says 'for the proton' but it should refer to the Xi^- baryon, consistent with the '72 exceptional configurations' and the bottom row of Fig. 6; this is a copy-paste error.
Circularity Check
No significant circularity: the lattice polarisabilities are fit coefficients of a defined energy-shift relation, and the chiral comparison is not forced by the new data.
full rationale
The derivation of the magnetic polarisabilities is a direct lattice measurement: Eq. (6) defines beta as the coefficient of k_d^2 in the energy shift, and the values in Table II are obtained by fitting that relation to correlation-function plateaus. Nothing in this step is defined in terms of the final beta. The exceptional-configuration removal uses a statistical outlier criterion (Eq. 13) on effective-energy uncertainties at early times, not on the final polarisability; while the stopping point is chosen by inspecting animations of the energy-shift plateau (Secs. IV.D, V), the paper explicitly concedes this is difficult and that beta is sensitive to continued removal ('the magnetic polarisability is sensitive to continued configuration removal, tending to increase to larger values', Sec. V.A). That is a systematic-error limitation, not a circular reduction: the plateau is the measured quantity, and choosing where it resolves is a fitting-window choice, not an equation that forces beta. The chiral comparison in Sec. VI uses the authors' prior Ref. [5] for the heavier points and the finite-volume/electro-quench corrections, but the paper states the new light point has minimal impact on the fit ('the inclusion of values at the lightest quark mass in the chiral fit has minimal impact on the result of the extrapolation'), so the agreement with chiral EFT is an independent consistency check rather than a self-citation-forced result. No step reduces by construction to its inputs.
Assumptions & free parameters
free parameters (3)
- Number of exceptional configurations removed =
p: 773, n: 60, Σ+: 680, Σ−: 225, Ξ0: 921, Ξ−: 72
- Fit window end tmax per baryon =
p:25, n:22, Σ+:26, Σ−:25, Ξ0:29, Ξ−:28
- nσ,min threshold for exceptional config removal =
Various values such as 3.5, 4.0, 4.5, 5.0, 6.0, 6.25, 6.75, 7.25, 7.5, 9.25, 9.5
assumptions (5)
- domain assumption Wilson-clover fermion action with additive mass renormalisation
- domain assumption Background field method: energy shift formula Eq. (2) with Landau term and polarisability term
- domain assumption Single-state dominance for t ≥ tmin=20
- domain assumption Electro-quenched approximation (sea quarks blind to the background magnetic field)
- domain assumption Chiral perturbation theory as benchmark
Cite this review
Pith. "Pith review of Magnetic polarisability of octet baryons near the physical quark-mass point." pith.science (2026). https://pith.science/paper/ZJN5XMU7
@misc{pith2026241208960,
author = {Pith},
title = {Pith review of: Magnetic polarisability of octet baryons near the physical quark-mass point},
year = {2026},
howpublished = {\url{https://pith.science/paper/ZJN5XMU7}},
note = {Machine review of arXiv:2412.08960}
}
read the original abstract
The magnetic polarisabilities of octet baryons are calculated close to the physical quark-mass point using the background field method in lattice QCD. This first calculation draws on the identification and elimination of exceptional configurations that have hindered previous attempts. The origin of the exceptional configuration problem lies in the use of a Wilson-type fermion action on electro-quenched gauge field configurations, where the dynamical-fermion gauge-field generation algorithm the electric charges of the quarks. Changes in the fermion determinant that would suppress some gauge fields in the background magnetic field are neglected, leaving improbable gauge fields that generate large additive mass renormalisations which manifest as significant outliers in correlation-function distributions. An algorithm for the systematic identification and removal of these exceptional configurations is described. We find the light up and down quarks to be problematic, particularly the up quark with its larger electric charge. The heavier mass of the strange quark protects the hyperon correlation functions to some extent. However, these also benefit from the removal of exceptional configurations. In many cases, the magnetic polarisability is calculated with good precision. We find our results to be in accord with the behaviour anticipated by chiral perturbation theory.
Figures
Figures from the paper (4 more)
Reference graph
Works this paper leans on
-
[5]
Neutron magnetic polarisability with Landau mode operators
R. Bignell, J. Hall, W. Kamleh, D. Leinweber, and M. Burkardt, Phys. Rev. D 98, 034504 (2018), arXiv:1804.06574 [hep-lat]
work page Pith review arXiv 2018
-
[1]
Construct a periodic list of correlation functions of length N + M − 1
-
[2]
Given the ith subset of correlation functions {Ci, ..., Ci+M −1} determine the effective energy Ei(t) and uncertainty ∆ Ei(t)
-
[3]
Repeat for i ∈ [1, N]. We emphasise that the ith elements of {E1, · · ·EN } and {∆E1, · · ·∆EN } are the effective energy and associated uncertainty of the corresponding subset of correlation functions in the ith subset. We’ll refer to evolution of the index i as evolution in configuration time. It is also important to emphasise that these quantities are ...
work page 2009
-
[6]
R. Bignell, W. Kamleh, and D. Leinweber, Phys. Rev. D 101, 094502 (2020), arXiv:2002.07915 [hep-lat]
arXiv 2020
-
[7]
R. Bignell, W. Kamleh, and D. Leinweber, Phys. Rev. D 100, 114518 (2019), arXiv:1910.14244 [hep-lat]
arXiv 2019
-
[8]
Magnetic Polarisability of Octet Baryons via Lattice QCD
T. Kabelitz, R. Bignell, W. Kamleh, and D. Leinweber, (2024), arXiv:2407.01064 [hep-lat]
work page Pith review arXiv 2024
-
[9]
S. Aoki et al. (PACS-CS), Phys. Rev. D 79, 034503 (2009), arXiv:0807.1661 [hep-lat]
arXiv 2009
Show all 29 references
-
[10]
Martinelli, G
G. Martinelli, G. Parisi, R. Petronzio, and F. Rapuano, Phys. Lett. B 116, 434 (1982)
1982
-
[11]
Primer, W
T. Primer, W. Kamleh, D. Leinweber, and M. Burkardt, Phys. Rev. D 89, 034508 (2014), arXiv:1307.1509 [hep- lat]
2014 arXiv
-
[12]
Akaike, IEEE Trans
H. Akaike, IEEE Trans. Automatic Control 19, 716 (1974)
1974
- [13]
-
[14]
B. J. Menadue, A Study of theΛ(1405) in Lattice QCD, Ph.D. thesis, Adelaide U. (2018)
2018
-
[15]
Di Renzo, in 40th International Symposium on Lattice Field Theory (2024) arXiv:2401.14752 [hep-lat]
F. Di Renzo, in 40th International Symposium on Lattice Field Theory (2024) arXiv:2401.14752 [hep-lat]
2024 arXiv
-
[16]
B. J. Menadue, W. Kamleh, D. B. Leinweber, and M. S. Mahbub, Phys. Rev. Lett. 108, 112001 (2012), arXiv:1109.6716 [hep-lat]
2012 arXiv
-
[17]
D. B. Leinweber, R. M. Woloshyn, and T. Draper, Phys. Rev. D 43, 1659 (1991)
1991
-
[18]
Deshmukh and B
A. Deshmukh and B. C. Tiburzi, Phys. Rev. D97, 014006 (2018), arXiv:1709.04997 [hep-ph]
2018 arXiv
-
[19]
Bruckmann, G
F. Bruckmann, G. Endrodi, M. Giordano, S. D. Katz, T. G. Kovacs, F. Pittler, and J. Wellnhofer, Phys. Rev. D 96, 074506 (2017), arXiv:1705.10210 [hep-lat]
2017 arXiv
-
[20]
Chang, W
E. Chang, W. Detmold, K. Orginos, A. Parreno, M. J. Savage, B. C. Tiburzi, and S. R. Beane (NPLQCD), Phys. Rev. D 92, 114502 (2015), arXiv:1506.05518 [hep- lat]
2015 arXiv
-
[21]
Morningstar and M
C. Morningstar and M. J. Peardon, Phys. Rev. D 69, 054501 (2004), arXiv:hep-lat/0311018
2004 arXiv
-
[22]
Gusken, Nucl
S. Gusken, Nucl. Phys. B Proc. Suppl. 17, 361 (1990)
1990
-
[23]
D. S. Roberts, P. O. Bowman, W. Kamleh, and D. B. Leinweber, Phys. Rev. D 83, 094504 (2011), arXiv:1011.1975 [hep-lat]
2011 arXiv
-
[24]
B. C. Tiburzi and S. O. Vayl, Phys. Rev. D 87, 054507 (2013), arXiv:1210.4464 [hep-lat]
2013 arXiv
-
[25]
Boinepalli, W
S. Boinepalli, W. Kamleh, D. B. Leinweber, A. G. Williams, and J. M. Zanotti, Phys. Lett. B 616, 196 (2005), arXiv:hep-lat/0405026
2005 arXiv
-
[26]
J. M. Zanotti, S. O. Bilson-Thompson, F. D. R. Bonnet, P. D. Coddington, D. B. Leinweber, A. G. Williams, J. B. Zhang, W. Melnitchouk, and F. X. Lee (CSSM Lattice), Phys. Rev. D 65, 074507 (2002), arXiv:hep-lat/0110216
2002 arXiv
-
[27]
J. M. M. Hall, D. B. Leinweber, and R. D. Young, Phys. Rev. D 89, 054511 (2014), arXiv:1312.5781 [hep-lat]
2014 arXiv
-
[28]
Kamleh, PoS LA TTICE2022, 339 (2023), arXiv:2302.00850 [hep-lat]
W. Kamleh, PoS LA TTICE2022, 339 (2023), arXiv:2302.00850 [hep-lat]
2023 arXiv
-
[29]
https://www2.jldg.org/
-
[30]
J. A. Mickley, W. Kamleh, and D. B. Leinweber, Phys. Rev. D 110, 034516 (2024), arXiv:2405.10670 [hep-lat]. Supplemental Material: Magnetic polarisability energy shift animations Thomas Kabelitz, Waseem Kamleh, and Derek Leinweber Special Research Centre for the Subatomic Stru...
2024 arXiv
Reviewed August 11, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.