{"id":"994e79e4-f6e6-4e37-8c4f-6d0c951e4714","arxiv_id":"2502.00331","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"A first-order confinement transition is observed at am=0.4 in the one-flavor SU(4) Hyper Stealth Dark Matter theory.","lead":"Lattice simulations of an SU(4) dark gauge theory with one fermion flavor find evidence for a first-order confinement transition at heavy dark quark mass. If confirmed, this composite dark matter model could generate a detectable gravitational wave background from the early universe.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The first-order claim at am=0.4 rests on a two-volume FSS fit: 16^3 is excluded for deviating, with M5=1.8 vs 1.5, leaving only 24^3 and 32^3. A matched fourth volume is needed to test whether b=1.02(19) is asymptotic.","rationale":"The reader's weakest_assumption identifies exactly the issue I find most load-bearing: the am = 0.4 first-order claim depends on a finite-size scaling exponent fit to only two volumes, with the third volume excluded for reasons tied to an action-parameter mismatch (M5 = 1.8 vs 1.5). My independent reading of Sec. 4 and Fig. 7 agrees that neither the two-point exponent nor the single-volume double-peaked histogram alone is sufficient to establish asymptotic first-order behavior. I do not see a more serious problem: the paper is careful about thermalization (hot/cold streams), controls residual chiral symmetry breaking to sub-percent level, and explicitly presents the result as preliminary. The correct disposition is therefore to keep the existing CONDITIONAL verdict: accept the result as a valuable lattice data point, but do not treat the first-order nature at am = 0.4 as established until a matched third or fourth volume confirms the scaling and, ideally, a continuum extrapolation is attempted. This does not change the reader's verdict.","tokens_in":8350,"tokens_out":4259,"duration_ms":43675,"concrete_test":"Run a matched 40^3 x 8 ensemble at am = 0.4, Nt = 8, M5 = 1.5, tuning beta to the Polyakov-loop susceptibility peak (e.g., near beta = 11.035) with statistics comparable to the 32^3 run. Refit chi_max = c N_s^{3b} using 24^3, 32^3, and 40^3, with and without the 16^3 point. If the new point lies on the 24/32 power law and b remains near 1, the first-order inference is supported; if it deviates by more than the fit uncertainty or b drops materially below 1, the two-volume exponent was not asymptotic and the first-order conclusion is not established.","verdict_should_be":"UNCHANGED","load_bearing_attack":"In Sec. 4, the transition order is decided by fitting the Polyakov-loop susceptibility peak to chi_max proportional to N_s^{3b}, with b = 1 expected for a first-order transition. At am = 0.4 the quoted b = 1.02(19) is obtained from only the 24^3 x 8 and 32^3 x 8 ensembles; the 16^3 x 8 point is not used because it deviates significantly and, per footnote 1, was generated with M5 = 1.8 instead of 1.5. The remaining fit has zero degrees of freedom: any two points define a power law, so the reported error reflects only the uncertainty of the individual Gaussian peak fits, not the stability of the scaling form. The same two-point situation applies to the crossover exponents b ~ 0.46(10) at am = 0.2 and 0.49(19) at am = 0.3. The double-peaked |PL| histogram (Fig. 7) is shown on 24^3 only; a two-peak structure at a single spatial volume can appear in a sufficiently sharp crossover and does not independently establish a true first-order transition unless the separation persists or sharpens with increasing volume. The paper is appropriately labeled preliminary, but the central claim is load-bearing for the gravitational-wave motivation and is currently supported by one two-point exponent plus one single-volume histogram.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","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.","tokens_in":8584,"tokens_out":6465,"duration_ms":54552,"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":[{"comment":"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.","section":"Sec. 4 (FSS ansatz)"},{"comment":"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.","section":"Sec. 4 (crossover exponents)"},{"comment":"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.","section":"Sec. 4, Fig. 7"}],"minor_comments":[{"comment":"The phrase 'over at least in some finite range' contains a grammatical error; rephrase to 'over at least a finite range.'","section":"Abstract"},{"comment":"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.","section":"Sec. 4"},{"comment":"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.","section":"Footnote 1"},{"comment":"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.","section":"Figure 5 caption"},{"comment":"The phrase 'coincidental peak locations' should be 'coincident peak locations' or 'coinciding peak locations.'","section":"Sec. 5"},{"comment":"The observation that the topological susceptibility shows no peak structure around beta_crit is not discussed; a sentence explaining its interpretation would be useful.","section":"Sec. 2 / Fig. 3"},{"comment":"Reference [1] is cited as a preprint (2412.14540); if it has been published by the time of final submission, please update the citation.","section":"References"}],"recommendation":"major_revision","confidential_remarks":"The paper is a proceedings-style contribution, and the collaboration is transparent about the preliminary nature of the analysis. The load-bearing issue is the two-volume FSS fit; this is a common limitation for such projects, and the requested revision (added volumes or a matched 16^3) is feasible within the project's scope. I do not see any sign of circular reasoning in the quenched benchmark or the residual-mass control. If the authors are unable to add volumes, the central claim should be downgraded to a 'candidate first-order' statement rather than a definitive finding."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Quick take: this is a solid preliminary proceedings report, not a finished result. The new thing is the first N_t=8 scan of transition order versus bare quark mass in one-flavor SU(4), with finite-size scaling exponents b for am=0.2, 0.3, 0.4 and quenched. The quenched benchmark b=1.009(23) is a good sanity check, and the paper is honest about its preliminary status: it reports two-start thermalization, residual mass control, topological charge tunneling, and does not oversell the chiral observables.\n\nThe load-bearing claim for the gravitational-wave motivation is that the transition is first-order at am=0.4. I think that is plausible but not yet established. The exponent b=1.02(19) comes from the 24^3 and 32^3 volumes only; the 16^3 point is excluded because it was generated with M5=1.8 instead of 1.5 and deviates from the scaling ansatz. Two points define a power law, so the quoted error reflects the uncertainty in the individual Gaussian peak fits, not the stability of the scaling form. The same two-volume situation applies to the crossover exponents at am=0.2 and 0.3. The double-peaked |PL| histogram is shown on 24^3 only; a sufficiently sharp crossover can also produce two peaks at a single volume. Together, b near 1 plus the two-peak structure is suggestive, but it does not independently confirm a true first-order transition. What is needed is a matched fourth volume (or the corrected 16^3 ensembles), ideally larger volumes, and eventually a continuum extrapolation. The paper is appropriately labeled preliminary, and these limitations are normal for a proceedings.\n\nCredit where it is due: the residual mass is sub-percent for the parameters used, the authors correctly note that the one-flavor theory has no chiral transition to look for, and the quenched comparison shows the method reproduces a known first-order transition. The citation pattern is proper: it builds on the LSD collaboration's earlier report and on one-flavor QCD studies, with no obvious omissions. The paper is a genuine new numerical result, even if it is an incremental step within an ongoing program.\n\nThis is for readers working on composite dark matter, gravitational-wave phenomenology, and lattice finite-temperature methods. It deserves a serious referee rather than a desk reject, precisely because the central claim is important and the current evidence, while suggestive, is not yet conclusive. My own verdict would be conditional: accept as a preliminary proceedings contribution, with the understanding that the first-order determination needs more volumes and, eventually, a continuum limit.","headline":"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.","tokens_in":9198,"tokens_out":1940,"would_cite":true,"duration_ms":23098,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"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.","keywords":["lattice gauge theory","SU(4) gauge theory","composite dark matter","Möbius domain wall fermions","finite-temperature phase transition","Polyakov loop","gravitational waves","first-order transition"],"falsifier":"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.","tokens_in":8080,"feed_emoji":"🌌","tokens_out":6250,"duration_ms":58025,"temperature":0.7,"pith_summary":"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.","feed_headline":"SU(4) dark matter transition is first-order at moderate quark mass","feed_subtitle":"Lattice simulation shows the Polyakov-loop signal growing with volume, a signature that could seed gravitational waves.","key_machinery":"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.","core_discovery":"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.","pith_inferences":["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."],"forward_implications":["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)."],"supporting_citations":[{"why":"Proposes the Hyper Stealth Dark Matter model and identifies the composite baryon as the dark matter candidate.","marker":"[1]"},{"why":"Previous report on this theory that the present study continues, providing the baseline ensembles and context.","marker":"[6]"},{"why":"Introduces the domain-wall fermion discretization underlying the Möbius variant used in the simulations.","marker":"[7]"},{"why":"Defines the Möbius domain wall fermion algorithm that the paper uses for the finite-temperature ensembles.","marker":"[8]"},{"why":"Supplies the exact one-flavor pseudofermion action used to generate gauge configurations.","marker":"[9]"},{"why":"One-flavor QCD analogue that motivates the expectation of first-order transitions at heavy masses and crossover at light masses.","marker":"[11]"},{"why":"Supplies the finite-size scaling ansatz for first-order phase transitions, $\\chi_{\\max}\\propto N_s^{3b}$ with $b=1$.","marker":"[15]"},{"why":"Provides the $b=1$ finite-size scaling signature for first-order deconfinement in a lattice gauge theory.","marker":"[16]"}],"fun_headline_variants":["First-order dark phase transition could seed gravity waves","SU(4) dark matter deconfinement found first-order on lattice","Lattice shows first-order transition in stealth dark matter","One-flavor SU(4) transition first-order, potential GW source","Möbius fermions reveal first-order dark matter phase change"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"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.","fun_headline_variants_meta":{"raw":{"variants":["First-order dark phase transition could seed gravity waves","SU(4) dark matter deconfinement found first-order on lattice","Lattice shows first-order transition in stealth dark matter","One-flavor SU(4) transition first-order, potential GW source","Möbius fermions reveal first-order dark matter phase change"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000172,"raw_usage":{"total_tokens":1285,"prompt_tokens":967,"completion_tokens":318,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":583,"completion_tokens_details":{"reasoning_tokens":231}},"tokens_in":583,"tokens_out":318,"duration_ms":3726,"temperature":1.0,"reasoning_tokens":231,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-09T19:23:42.766378+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"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.","supporting_citations":[{"cited_title":"Exploring Composite Dark Matter with an SU(4) gauge theory with 1 fermion flavor","cited_arxiv_id":"2402.07362","evidence_quote":"Previous report on this theory that the present study continues, providing the baseline ensembles and context."},{"cited_title":"Binder and D","cited_arxiv_id":null,"evidence_quote":"Supplies the finite-size scaling ansatz for first-order phase transitions, $\\chi_{\\max}\\propto N_s^{3b}$ with $b=1$."},{"cited_title":"Fukugita, M","cited_arxiv_id":null,"evidence_quote":"Provides the $b=1$ finite-size scaling signature for first-order deconfinement in a lattice gauge theory."}],"review_version":1}