Pith. sign in

REVIEW 3 major objections 7 minor 16 references

Finite Temperature Transition in Hyper Stealth Dark Matter using M\"{o}bius Domain Wall Fermions

T0 review · 3 major / 7 minor · reviewed 2026-08-09 · deepseek-v4-flash

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗

Pith's one-line read The paper claims that the deconfinement transition of one-flavor SU(4) dark gauge theory is first-order at $am=0.4$, a result that would make the theory a candidate source of gravitational waves.

desk verdict A clean, honestly-labeled preliminary lattice scan of the HSDM transition; the first-order claim at am=0.4 is plausible but rests on a two-volume FSS fit and a single-volume histogram. read the letter →

arxiv 2502.00331 v1 pith:QT73EJWE submitted 2025-02-01 hep-lat

Venkitesh Ayyar , Nobuyuki Matsumoto , Aaron S. Meyer , Sungwoo Park This is my paper · ORCID
classification hep-lat
keywords latticegaugetheorySU(4)compositedarkmatterMöbiusdomainwallfermionsfinite-temperaturephasetransitionPolyakovloopgravitationalwavesfirst-order
verification ladder T0 review T1 audit T2 compute T3 formal

The pith

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The reading

This paper investigates the finite-temperature deconfinement transition of one-flavor SU(4) gauge theory, a candidate dark matter model in which the dark matter particle is a composite bosonic baryon of four dark quarks. Using Möbius domain wall fermions on $N_t=8$ lattices with three spatial volumes, the authors measure the Polyakov loop, chiral condensate, and topological charge across a range of dark-quark masses from $am=0.01$ to the quenched limit. The paper's central claim is that the transition is first-order at the moderately heavy quark mass $am=0.4$, evidenced by a Polyakov-loop susceptibility peak growing as $N_s^3$ with exponent $b=1.02(19)$ and by a double-peaked $|P_L|$ histogram, while lighter masses show crossover behavior. This matters because a first-order transition in this theory could generate a detectable stochastic gravitational-wave background in the early universe.

What carries the argument

The central machinery is the finite-size scaling of the Polyakov loop susceptibility, $\chi_{|P_L|}\propto N_s^{3b}$, with a first-order transition predicted to have $b=1$. The argument uses Möbius domain wall fermions, a lattice discretization of chiral fermions, together with an exact one-flavor pseudofermion action, and it examines the Polyakov loop's complex-plane distribution and histogram to separate the first-order region from the crossover region. The residual mass from finite fifth-dimension extent is checked to be sub-percent, making the chiral symmetry breaking controlled.

What would settle it

Measure the Polyakov-loop susceptibility at $am=0.4$ on at least a third, larger volume such as $40^3\times 8$ and check whether the peak height grows as $N_s^3$; if the exponent falls toward the crossover value $b\simeq 0.5$ or the $|P_L|$ histogram becomes single-peaked, the first-order claim is contradicted. A demonstration that the two-volume scaling can be reproduced by a crossover model would also falsify the interpretation.

Watch

Extended reading notes

Core claim

The authors claim that one-flavor SU(4) with Möbius domain wall fermions has a first-order deconfinement transition at bare quark mass $am=0.4$, with the second-order endpoint located between $am=0.3$ and $am=0.4$. At $am=0.4$, the peak of the Polyakov-loop susceptibility scales linearly with spatial volume, $b=1.02(19)$ in $\chi_{\max}\propto N_s^{3b}$, and the histogram of $|P_L|$ at $\beta_{\rm crit}$ is double-peaked, both signatures of a first-order transition. At lighter masses the exponent is $b=0.46(10)$ for $am=0.2$ and $b=0.49(19)$ for $am=0.3$, consistent with crossover, while the quenched theory gives $b=1.009(23)$, consistent with first order. The residual chiral symmetry breaking from finite fifth-dimension extent is sub-percent, and no peak structure is observed in the topological susceptibility around the transition, consistent with the absence of a chiral phase transition in a one-flavor theory.

Load-bearing premise

The claim that the transition is first-order at $am=0.4$ rests on the assumption that the growth of the susceptibility peak from the $24^3$ lattice to the $32^3$ lattice already captures the asymptotic infinite-volume behavior, even though only two volumes are used.

Editorial extensions

If this is right

  • At $am=0.4$ the theory realizes a genuine first-order deconfinement transition, providing a concrete candidate for a gravitational-wave source from a strongly coupled dark sector.
  • The second-order endpoint of the transition sits between $am=0.3$ and $am=0.4$, defining the boundary between crossover and first-order behavior in this one-flavor theory.
  • The steep change of the chiral condensate near $\beta_{\rm crit}$ and the coincident peaks of chiral and Polyakov-loop susceptibilities suggest that chiral and deconfinement transitions are tied in this model.
  • Zero-temperature generation at $\beta_{\rm crit}$ is needed to set the scale and determine whether the composite baryon mass is indeed a few GeV.
  • The quenched limit is first-order with the center symmetry broken spontaneously, as expected for pure SU(4).

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

  • If the first-order line continues past $am=0.4$, the theory could produce stronger gravitational-wave signals; a run at $am\ge 0.5$ would test this directly.
  • The apparent coincidence of the Polyakov and chiral susceptibility peaks, if confirmed with improved statistics, would be evidence that deconfinement drives the chiral response in a one-flavor theory despite the axial anomaly.
  • A three-volume analysis at $am=0.4$ that includes a $40^3$ lattice would either confirm asymptotic scaling or reveal that the two-volume exponent is a crossover artifact.
  • The sub-percent residual mass suggests that extrapolating to infinite fifth-dimension extent is unlikely to change the order of the transition, making the first-order claim robust to the chiral discretization.
Share X Bluesky LinkedIn Reddit HN

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Referee Report

3 major / 7 minor

Summary. The paper reports a lattice study of the finite-temperature confinement transition in the one-flavor SU(4) 'Hyper Stealth Dark Matter' theory using Möbius domain-wall fermions at N_t=8. The authors measure the Polyakov loop, chiral condensate, and topological charge for bare quark masses am = 0.01-0.4 and in the quenched limit, using spatial volumes 16^3, 24^3, and 32^3. They determine the transition coupling beta_crit from the Polyakov-loop susceptibility peak and classify the transition order from the finite-volume scaling exponent b in chi_max proportional to N_s^{3b}, with b=1 expected for first-order. They report b=1.02(19) at am=0.4 (from 24^3 and 32^3), b=0.46(10) at am=0.2, b=0.49(19) at am=0.3, and b=1.009(23) in the quenched theory, and they show a double-peaked |PL| histogram at am=0.4 on 24^3. The central claim is that the transition is first-order at a 'large but moderate' quark mass am=0.4, with a second-order point somewhere between am=0.3 and 0.4, making the model a potential source of early-universe gravitational waves.

Significance. If correct, the claim that one-flavor SU(4) has a first-order deconfinement transition at moderate fermion mass is phenomenologically important: it would identify a concrete composite dark matter model whose confinement transition can generate a detectable stochastic gravitational-wave background. The paper's methodological strengths are its benchmark against the quenched theory, where the known first-order behavior is reproduced with b=1.009(23), and its careful control of the residual chiral symmetry breaking (sub-percent m_res/m). The work is honestly labeled as preliminary, and the text does not over-interpret the data beyond the central claim. However, the central classification at am=0.4 currently rests on a two-volume FSS fit and a single-volume histogram, which is a correct but fragile support base.

major comments (3)
  1. [Sec. 4 (FSS ansatz)] The first-order claim at am=0.4 is supported by the exponent b=1.02(19) obtained from a fit of chi_max proportional to N_s^{3b} to only the 24^3 x 8 and 32^3 x 8 ensembles. With two points the fit has zero degrees of freedom, so the quoted uncertainty reflects only the Gaussian peak-height errors and not the stability of the power-law form. The 16^3 x 8 point is excluded because, per footnote 1, it was generated with M5=1.8 rather than 1.5; this parameter mismatch confounds the exclusion with a genuine finite-volume effect. Please add a matched 16^3 dataset, add at least a fourth volume, or explicitly state that the two-point fit is a preliminary indicator rather than a determination.
  2. [Sec. 4 (crossover exponents)] The same two-volume limitation applies to the crossover classifications: b=0.46(10) at am=0.2 and b=0.49(19) at am=0.3 are quoted after the 16^3 point is removed 'as it deviates significantly from the scaling ansatz.' Excluding a point because it deviates from the ansatz under test is circular, and the M5 mismatch (footnote 1) provides an alternative explanation for the deviation. Consequently the claimed second-order endpoint between am=0.3 and 0.4 is not yet robustly located.
  3. [Sec. 4, Fig. 7] The double-peaked |PL| histogram at am=0.4 is shown only for the 24^3 x 8 lattice. At a single volume, a sufficiently sharp crossover can produce a bimodal distribution, so this plot does not independently confirm a first-order transition. A volume-dependence plot (e.g., 16^3 vs 24^3 vs 32^3, with matched M5) showing that the two peaks separate and the barrier grows would turn this into supporting evidence.
minor comments (7)
  1. [Abstract] The phrase 'over at least in some finite range' contains a grammatical error; rephrase to 'over at least a finite range.'
  2. [Sec. 4] In the sentence 'we do not use the results of 163 x 8 in estimating b for am=0.2, 0.3 as it is deviates significantly from the scaling ansatz,' the phrase 'is deviates' is a typo; also clarify that the deviation is from the expected scaling behavior rather than from the ansatz itself.
  3. [Footnote 1] The expectation that the M5=1.8 vs 1.5 difference is 'small' should be quantified (e.g., by comparing m_res or |PL| on a small test ensemble), because the parameter mismatch is the stated reason for excluding a whole volume from the FSS analysis.
  4. [Figure 5 caption] The caption says 'Green, blue, and red points show the results' but the plotted points are colored by volume (16^3, 24^3, 32^3) with a legend that may not correspond to those colors; please check that the caption matches the actual plot colors.
  5. [Sec. 5] The phrase 'coincidental peak locations' should be 'coincident peak locations' or 'coinciding peak locations.'
  6. [Sec. 2 / Fig. 3] The observation that the topological susceptibility shows no peak structure around beta_crit is not discussed; a sentence explaining its interpretation would be useful.
  7. [References] Reference [1] is cited as a preprint (2412.14540); if it has been published by the time of final submission, please update the citation.

Circularity Check

0 steps flagged · score 0.0 of 10

No significant circularity: the transition-order claim is a direct lattice measurement benchmarked against the quenched theory.

full rationale

The paper is a numerical lattice study whose central claim, that the SU(4) one-flavor theory has a first-order deconfinement transition at am = 0.4, is obtained by direct measurement of the Polyakov loop susceptibility and histogram. The finite-size scaling ansatz chi_max ∝ N_s^{3b} with b = 1 for first-order transitions is taken from standard literature [15,16], not from the authors' own prior work. The method is independently validated on the quenched theory, where b = 1.009(23) reproduces the known first-order transition. The residual-mass check uses the axial Ward identity and only establishes that MDWF chiral symmetry breaking is controlled; it does not enter the transition-order determination. No fitted parameter is renamed as a prediction, and no uniqueness theorem or load-bearing result is imported from a self-citation. The acknowledged limitation that the b = 1.02(19) estimate at am = 0.4 relies on only two spatial volumes (24^3 and 32^3), with 16^3 excluded due to the M5 = 1.8 versus 1.5 mismatch, is a statistical robustness concern, not a circularity: the two volumes do not by construction determine the quoted exponent, and the conclusion is additionally supported by the double-peaked |PL| histogram. The self-citations to Refs. [1] and [6] provide model motivation and prior context only, and the central numerical result is self-contained against the quenched benchmark. Therefore no circular step is present.

Assumptions & free parameters 0 free parameters · 4 assumptions · 0 invented entities

No free parameters are fitted to produce the claimed result; the bare quark masses and beta are simulation inputs. The analysis relies on the standard finite-size scaling ansatz for first-order transitions, the assumption that the Polyakov loop is a useful order parameter with one dynamical flavor, the validity of the MDWF approximation demonstrated by the small residual mass, and the representation of the heavy-mass limit by the quenched theory. No new entities are introduced.

assumptions (4)
  • domain assumption Finite-size scaling ansatz chi_max proportional to N_s^{3b}, with b=1 for a first-order transition (Refs. 15,16).
    Invoked in Sec. 4 to translate the volume dependence of the Polyakov loop susceptibility peak into a statement about the order of the transition.
  • domain assumption The Polyakov loop is a valid order parameter for deconfinement in the presence of a dynamical fermion that explicitly breaks the Z4 center symmetry.
    Stated in Sec. 4; the authors argue it should still work for large quark masses, but this is an assumption for the unquenched cases.
  • domain assumption The Möbius domain-wall fermion action with L_s=16 and the chosen M5 values gives a negligible residual mass, so the simulated theory is close to the target one-flavor theory.
    Checked in Sec. 2 with the eigenvalue and residual mass analysis; sub-percent residual mass, but it remains a discretization assumption.
  • domain assumption The quenched theory (am=infinity) correctly represents the heavy-mass limit of the one-flavor theory.
    Used in Sec. 4 as the infinite-mass reference point for comparing the FSS exponent.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Finite Temperature Transition in Hyper Stealth Dark Matter using M\"{o}bius Domain Wall Fermions." pith.science (2026). https://pith.science/paper/QT73EJWE

@misc{pith2026250200331,
  author       = {Pith},
  title        = {Pith review of: Finite Temperature Transition in Hyper Stealth Dark Matter using M\"obius Domain Wall Fermions},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/QT73EJWE}},
  note         = {Machine review of arXiv:2502.00331}
}
abstract

The first-order confinement transition of a strongly coupled composite dark matter theory can provide a possible source of gravitational waves in the early universe. In this work, on behalf of the Lattice Strong Dynamics (LSD) Collaboration, we present our recent investigation on the finite temperature confinement transition of the one-flavor SU(4) dark gauge theory named Hyper Stealth Dark Matter (HSDM). The dark matter candidate in this theory is a composite bosonic baryon and can have a remarkably low mass of a few GeV. We expect the finite temperature transition to be first-order over at least in some finite range of fermionic masses and to be a potential source of observable gravitational radiation. The finite temperature simulation of one-flavor SU(4) is done by using M\"{o}bius Domain wall fermions. The order of the transition and its fermionic mass dependence are explored by monitoring the Polyakov loop, chiral condensate and topological charge using three lattice volumes at $N_t=8$.

Figures

Figures reproduced from arXiv: 2502.00331 by the authors.

Figure 1
Figure 1. The ten smallest magnitude eigenvalues of the four-dimensional hermitian Dirac operator: 𝐻4 (𝑀5) = 𝛾5𝐷̸ 𝑊 (−𝑀5), calculated with ten thermalized configurations for 𝛽 = 10.6, 10.8, 11.0 (from left to right) with the volume 243 ×8 and the bare quark mass 𝑎𝑚 = 0.1. The red vertical dotted line indicates 𝑎𝑀5 = 1.8, which is the parameter used in the HMC for 𝑎𝑚 = 0.1 and 0.4. 0 5 10 15 20 x 0.00000 0.00008 0.00016 0.0002… view at source ↗
Figure 2
Figure 2. (Left) The effective residual mass 𝑎𝑚eff res(𝑥) = ⟨𝐽5𝑞 (𝑥)𝑃(0)⟩/⟨𝑃(𝑥)𝑃(0)⟩ as a function of the spatial separation 𝑥, evaluated at 𝛽crit for each quark mass. (Right) The residual mass 𝑚res determined from the midpoint of 𝑎𝑚eff res(𝑥), plotted as a function of 𝛽. Both plots are calculated with the 243 × 8 ensembles. given in Sec. 4. The red vertical line in the figure marks 𝑀5 = 1.8, which is used in the HMC to gener… view at source ↗
Figure 3
Figure 3. Monte Carlo time history of the Wilson-flowed topological charge 𝑄𝑊, and the correspond￾ing histogram calculated at 𝛽crit for the masses: (left) 𝑎𝑚 = 0.1 and (middle) 𝑎𝑚 = 0.4. In the right panel, we show the topological susceptibility as a function of 𝛽 for the six different quark masses: 𝑚 = 0.01, 0.05, 0.1, 0.2, 0.3, 0.4. Vertical dotted lines are drawn at 𝛽 = 𝛽crit to indicate the transition point for each mass … view at source ↗
Figures from the paper (5 more)
Figure 4
Figure 4. Figure 4: Polyakov Loop absolute value versus 𝛽 on 243 × 8 ensembles. for 𝑎𝑚 = 0.4, and 𝑏 = 1.009(23) for the quenched theory. We do not use the results of 163 × 8 in estimating 𝑏 for 𝑎𝑚 = 0.2, 0.3 as it is deviates significantly from the scaling ansatz, in relation to which we …
Figure 5
Figure 5. Figure 5: The Polyakov loop susceptibility normalized by the spatial volume, 𝜒|PL|/𝑉, where 𝑉 = 𝑁 3 𝑠 . Green, blue, and red points show the results for 163 × 8, 243 × 8, and 323 × 8, respectively. The ensembles with hot and cold starts are combined in this plot. As a preliminar…
Figure 6
Figure 6. Figure 6: The scattered plot of the spatially averaged Polyakov loop in the complex plane for the masses 𝑎𝑚 = 0.2, 0.4, and the quenched theory (from top to the bottom rows) on the 243 ×8 lattice. We choose three 𝛽 values for each 𝑎𝑚: below, on top of, and above the coupling 𝛽cr…
Figure 7
Figure 7. Figure 7: The Monte Carlo time history of the absolute value of the Polyakov loop, |𝑃𝐿|, and its histogram at 𝛽crit for the masses below and above the second-order transition point: 𝑎𝑚 = 0.2, 0.4 on the 243 ×8 lattice. The red and blue colors represent the hot and cold start str…
Figure 8
Figure 8. Figure 8: (Top) The chiral condensate and the Polyakov loop as a function of 𝛽. (Bottom) The disconnected chiral susceptibility and the Polyakov loop susceptibility as a function of 𝛽. The lattice size is fixed to 243 ×8 in this figure, while the dark-quark mass is varied as 𝑎𝑚 …

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

16 extracted references · 4 canonical work pages

  1. [1]

    G. T. Fleming, G. D. Kribs, E. T. Neil, D. Schaich, and P. M. Vranas,Hyper Stealth Dark Matter and Long-Lived Particles, 2412.14540

    arXiv

  2. [2]

    115(2015), no

    P.Schwaller, GravitationalWavesfromaDarkPhaseTransition ,Phys.Rev.Lett. 115(2015), no. 18 181101, [1504.07263]

    arXiv 2015

  3. [3]

    Capriniet al., Detecting gravitational waves from cosmological phase transitions with LISA: an update, JCAP03 (2020) 024, [1910.13125]

    C. Capriniet al., Detecting gravitational waves from cosmological phase transitions with LISA: an update, JCAP03 (2020) 024, [1910.13125]

    arXiv 2020

  4. [4]

    Auclairet al.,Cosmology with the Laser Interferometer Space Antenna, Living Rev

    LISA Cosmology Working GroupCollaboration, P. Auclairet al.,Cosmology with the Laser Interferometer Space Antenna, Living Rev. Rel.26(2023), no. 1 5, [2204.05434]. 8 Finite Temperature Transition in HSDM using MDWF Sungwoo Park

    arXiv 2023

  5. [5]

    Escrivà, R

    A. Escrivà, R. Inui, Y. Tada, and C.-M. Yoo,LISA forecast on a smooth crossover beyond the standard model through the scalar-induced gravitational waves,Phys. Rev. D111(2025), no. 2 023528, [2404.12591]

    arXiv 2025

  6. [6]

    Exploring Composite Dark Matter with an SU(4) gauge theory with 1 fermion flavor

    LSDCollaboration, V. Ayyar,Exploring Composite Dark Matter with an SU(4) gauge theory with 1 fermion flavor, PoSLATTICE2023(2024) 102, [2402.07362]

    work page Pith review arXiv 2024

  7. [7]

    D. B. Kaplan,A Method for simulating chiral fermions on the lattice,Phys. Lett. B288 (1992) 342–347, [hep-lat/9206013]

    arXiv 1992

  8. [8]

    R. C. Brower, H. Neff, and K. Orginos,The Möbius domain wall fermion algorithm,Comput. Phys. Commun.220 (2017) 1–19, [1206.5214]

    arXiv 2017

Show all 16 references
  1. [9]

    Chen and T.-W

    TWQCD Collaboration, Y.-C. Chen and T.-W. Chiu,Exact Pseudofermion Action for Monte Carlo Simulation of Domain-Wall Fermion,Phys. Lett. B738 (2014) 55–60, [1403.1683]

    2014 arXiv

  2. [10]

    Boyle, A

    P. Boyle, A. Yamaguchi, G. Cossu, and A. Portelli,Grid: A next generation data parallel C++ QCD library, 1512.03487

    arXiv

  3. [11]

    Alexandrou, A

    C. Alexandrou, A. Borici, A. Feo, P. de Forcrand, A. Galli, F. Jegerlehner, and T. Takaishi, The Deconfinement phase transition in one flavor QCD,Phys. Rev. D60(1999) 034504, [hep-lat/9811028]

    1999 arXiv

  4. [12]

    R. D. Pisarski and F. Wilczek,Remarks on the Chiral Phase Transition in Chromodynamics, Phys. Rev. D29(1984) 338–341

    1984

  5. [13]

    RBC, UKQCDCollaboration, D. J. Antonioet al.,Localization and chiral symmetry in three flavor domain wall QCD, Phys. Rev. D77(2008) 014509, [0705.2340]

    2008 arXiv

  6. [14]

    CP-PACSCollaboration, A. A. Khanet al., Chiral properties of domain wall quarks in quenched QCD,Phys. Rev. D63(2001) 114504, [hep-lat/0007014]

    2001 arXiv

  7. [15]

    Binder and D

    K. Binder and D. P. Landau,Finite-size scaling at first-order phase transitions, Phys. Rev. B 30(1984), no. 3 1477

    1984

  8. [16]

    Fukugita, M

    M. Fukugita, M. Okawa, and A. Ukawa,Order of the Deconfining Phase Transition in SU(3) Lattice Gauge Theory, Phys. Rev. Lett.63(1989) 1768. 9

    1989

Pith tools

Reviewed August 9, 2026 · model on record in the stance chip above.