{"id":"a67a465e-7ca3-4165-9a9b-104b212cc737","arxiv_id":"2412.02348","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":7.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":5,"one_line_summary":"A first-principles lattice calculation determines the mass of the gluino-glue bound state in N=1 supersymmetric Yang-Mills theory in the large-N limit.","lead":"Researchers ran large-scale lattice simulations of a supersymmetric gauge theory and measured the mass of its lightest bound state, the gluino-glue. This gives the first non-perturbative large-N value and a concrete prediction for the mass gap of the theory.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Chiral-continuum fit (3.11) is the weakest link: the paper's own mixed-term check shifts w0 M from 1.21 to 1.61, so the quoted 0.11 error is understated.","rationale":"I agree with the reader's weakest-assumption identification. The analysis has independent support in several places: the two mass extraction methods agree (Sec. 2.2), the finite-volume extrapolation (3.8) has reasonable chi2 and is stable under cuts, and the final SU(infinity) value lies on a smooth 1/N^2 curve with the SU(2) and SU(3) results (Sec. 3.3). The paper is also honest about the mixed-term check, which is what makes the concern concrete rather than speculative. The problem is that the quoted error 0.11 comes from the two-term fit, while the three-term fit moves the central value by 0.40. In lattice QCD, a systematic shift this large relative to the statistical error normally gets added in quadrature or the expanded model becomes the central result. Merely noting that the two values overlap after enlarging the error is not sufficient, because the default result and its error remain the narrow ones. The bias is in the direction of a smaller mass; if the mixed term is real at the level suggested by the fit, the large-N gluino-glue mass is closer to 1.6 w0, and the comparison with N=2 and N=3 data changes quantitatively (though the trend of increasing mass with N would survive). A finer lattice spacing is the cleanest discriminator because the mixed term is O(a), so it should decrease as a shrinks if the naive continuum extrapolation is valid. I therefore keep the CONDITIONAL verdict and recommend the new-ensemble check.","tokens_in":17046,"tokens_out":13237,"duration_ms":154229,"concrete_test":"Generate a new TEK ensemble at b = 0.355 (finer lattice spacing) for the same kappa values, perform the thermodynamic N=infinity extrapolation as in Sec. 3.1, and repeat the chiral-continuum fit of Sec. 3.2 both with Eq. (3.11) and with the mixed term c3 (a/sqrt(8t1))(8t1 m_pi^2) added. If the two fits still differ by ~0.4 in the SUSY-limit intercept (as they do with the present three b values), Eq. (3.11) is not supported and the final value must be quoted with a systematic error covering the mixed-term shift; if the two fits converge within combined errors, the simpler ansatz and the value 1.21(11) are confirmed.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central result (3.12) is the intercept at a=0, m_pi=0 of the fit in Eq. (3.11), which assumes a linear O(a) lattice artifact and a linear term in 8t1 m_pi^2 with no cross term. The data in Table 2 have the property that, within each b value, smaller m_pi is achieved at smaller a, so the points sit on a tilted line in the (a, 8t1 m_pi^2) plane. The SUSY-limit intercept is therefore an extrapolation with high leverage. The authors' own check at the end of Sec. 3.2 adds a mixed term c3 (a/sqrt(8t1))(8t1 m_pi^2); they find c3 = 0.21(20) and the intercept moves from 1.21 to 1.61(40). A shift of 0.40 is about 3.6 times the quoted error on (3.12). Since c3 is only 1 sigma from zero, the data do not exclude the mixed term, yet it changes the central value substantially. The statement that the two are 'compatible' is reached only by expanding the error by roughly a factor of four; it does not justify quoting 1.21(11) as the headline. The same sensitivity propagates to M/Lambda_NSVZ = 8.94(1.01) in Eq. (4.1). The concern is not that the calculation is wrong in some hidden way; the N=infinity limit, mass extraction, and comparison with finite-N data all look sound. It is that the final SUSY-limit value is not robust to a plausible, controlled extension of the fitting ansatz, and the quoted uncertainty does not reflect this.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"This paper computes the mass of the gluino-glue bound state in large-N N=1 supersymmetric Yang-Mills theory using twisted Eguchi-Kawai volume reduction and dynamical adjoint Wilson fermions. The authors use previously generated ensembles from Ref. [27] for N=169, 289, 361 at three values of the inverse 't Hooft coupling b and several gluino masses. They extract the mass from a GEVP-improved time correlator in the reduced model, verify it against an effective-mass plateau from the GEVP eigenvalues, and extrapolate to the thermodynamic limit with an exponential finite-volume ansatz. They then perform a combined chiral and continuum extrapolation, Eq. (3.11), to reach the SUSY limit and quote w0 M_gluino-glue = 1.21(11) at N=infinity, equivalently M/Lambda_NSVZ = 8.94(1.01). The large-N result is compared with SU(2) and SU(3) determinations and found to be consistent with a mild 1/N^2 growth.","tokens_in":17389,"tokens_out":9292,"duration_ms":95469,"significance":"If the quoted central value survives a more careful treatment of the extrapolation uncertainty, this is the first non-perturbative large-N value of the gluino-glue mass and a strong result for the large-N SUSY theory; the comparison with finite-N lattice data provides a useful test of 1/N^2 scaling. The paper has notable internal strengths: two consistent mass-determination procedures, a stable thermodynamic extrapolation across volume cuts, and a complete propagation of the result into scale-invariant units. The main weakness is the model dependence of the chiral-continuum extrapolation, which is currently not reflected in the quoted error.","major_comments":[{"comment":"The central SUSY-limit value is not robust to the inclusion of the mixed term c3 (a/sqrt(8t1))(8t1 m_pi^2) that the authors themselves test in the paragraph following Eq. (3.12). They find c3 = 0.21(20), which is only about one standard deviation from zero, yet the intercept moves from w0 M = 1.21 to 1.61, a shift of 0.40, about 3.6 times the quoted error 0.11. The statement that the two results are compatible is achieved only by accepting a four-times-larger uncertainty; it does not justify quoting 1.21(11) as the headline. Because within each b value the data points are strongly correlated, with smaller m_pi occurring at smaller a, the intercept of fit (3.11) has high leverage and the omission of the mixed term is a load-bearing modeling assumption. Please either include the c3 term in the central fit and report the corresponding value, or provide a systematic error that covers the variation between fit ansatze, and propagate the resulting uncertainty to Eqs. (3.13), (3.16)/(4.1) and the comparison in Sec. 3.3.","section":"Sec. 3.2, Eq. (3.11)-(3.12)"},{"comment":"The uncertainty quoted in Eq. (4.1), M/Lambda_NSVZ = 8.94 +/- 1.01, includes only the statistical errors from the mass fit and the inputs of Refs. [27,28]. The dominant systematic effect identified in Sec. 3.2, namely the change of the chiral-continuum fit ansatz, is not represented in this number. Without a quantitative systematic error from the extrapolation, the error budget is incomplete; the revised version should either state explicitly that Eq. (4.1) is conditional on the linear ansatz (3.11) or enlarge the error accordingly.","section":"Sec. 4, Eq. (4.1)"}],"minor_comments":[{"comment":"The phrase 'SUSY-restorting limit' contains a typo and should read 'SUSY-restoring limit'.","section":"Sec. 3, opening paragraph"},{"comment":"The b=0.345 rows and some b=0.350 rows have no entries in the N=169 column; the caption should state whether those ensembles were not generated or whether the values were omitted.","section":"Table 1"},{"comment":"For the 1/N^2 fit with three data points and two parameters, please report the number of degrees of freedom and the p-value; the statement that the chi-squared is 'very small' is not informative without these details.","section":"Eq. (3.19)"},{"comment":"Please clarify the mismatch between the SUSY-limit value of w0 M_gluino-glue reported in Table 1 and the value shown in Fig. 1 of Ref. [37]; a brief explanation would help the reader assess the reliability of the quoted SU(3) value.","section":"Sec. 3.3"},{"comment":"The large-N points are continuum-subtracted, while the SU(2) and SU(3) points are SUSY-limit values at a=0; using distinct symbols and an explicit legend entry for these two categories would improve readability.","section":"Fig. 3"}],"recommendation":"major_revision","confidential_remarks":"The central number may change after the requested robustness analysis, but this is a corrigible issue rather than a fundamental flaw. The fit-model sensitivity should be the main focus of the revision. I do not see a citation-practice concern beyond the expected reliance on the authors' own previous scale-setting and gluino-condensate determinations."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Colleague,\n\nThe headline result is a genuine first: the gluino-glue mass in the large-N limit of N=1 SYM, obtained from TEK volume reduction. The mass extraction is careful, the N=∞ thermodynamic extrapolation is stable across volume cuts, and the result is consistent with the finite-N trend from SU(2) and SU(3). That is real progress for lattice SUSY and large-N physics.\n\nThe soft spot is exactly where the stress-test note puts it: the SUSY-limit intercept from the fit in Eq. (3.11). The data points in Table 2 are tilted in the (a, m_pi^2) plane, so the intercept is high-leverage, and the authors' own check with the mixed term c3 moves the central value from 1.21 to 1.61(40). Saying these are 'compatible' only after inflating the error by a factor of four is technically true but not reassuring. The quoted 1.21(11) understates the systematic uncertainty from the fit form. This is not a hidden error in the simulation; it's a modeling assumption about the approach to the chiral-continuum limit. The paper should either adopt the more conservative value 1.61(40) as the headline, or justify the two-term fit more strongly and add a systematic error to 1.21.\n\nOther concerns are minor. The reliance on the authors' own scale setting and gluino condensate is not circular; those are independent published determinations. The lack of raw data and code makes reproduction a full lattice campaign, but that's normal for this field.\n\nThe paper is written clearly and the authors are transparent about the fit check. It deserves serious peer review. My recommendation: accept with revision, and make the fit-systematic explicit.","headline":"First large-N gluino-glue mass is a real result, but the quoted SUSY-limit error misses a fit-shape systematic that shifts the central value by roughly 0.4.","tokens_in":18014,"tokens_out":1966,"would_cite":true,"duration_ms":21961,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["11.15.Ha","12.60.Jv"],"model":"deepseek-v4-flash","headline":"The paper reports the first non-perturbative determination of the gluino-glue mass in large-$N$ $\\mathcal{N}=1$ supersymmetric Yang-Mills theory, obtaining $w_0M_{\\tilde g g}=1.21(11)$ in the SUSY limit.","keywords":["gluino-glue bound state","large-N limit","N=1 supersymmetric Yang-Mills theory","lattice gauge theory","twisted Eguchi-Kawai reduction","volume independence","chiral-continuum limit","gluino mass"],"falsifier":"Take the same $N=169,289,361$ ensembles and redo the chiral-continuum fit including the mixed term $c_3(a/\\sqrt{8t_1})(8t_1 m_\\pi^2)$ with higher statistics; the reported coefficient $c_3=0.21(20)$ is consistent with zero but shifts the central value from 1.21 to 1.61, so a determination of $c_3$ at the two-$\\sigma$ level would settle which SUSY-limit mass is right. Alternatively, data at a new, smaller lattice spacing would test whether the linear ansatz still describes the approach to the continuum limit.","tokens_in":16824,"feed_emoji":"⚛️","tokens_out":10566,"duration_ms":101502,"temperature":0.7,"pith_summary":"The paper establishes the first non-perturbative value of the gluino-glue mass in large-$N$ $\\mathcal{N}=1$ supersymmetric Yang-Mills theory, using numerical Monte Carlo simulations of a twisted volume-reduced lattice model. After taking the thermodynamic limit and then the chiral-continuum limit that restores supersymmetry, the authors obtain $w_0M_{\\tilde g g}=1.21(11)$, equivalently $M/\\Lambda_{\\mathrm{NSVZ}}=8.94\\pm1.01$. This matters because the gluino-glue state is expected to belong to the lightest supermultiplet, so the number fixes the mass gap of the infinite-color theory and provides a non-perturbative target for analytic approaches. The result is larger than the $SU(2)$ and $SU(3)$ values and is consistent with a mild growth in $1/N^2$.","feed_headline":"Infinite-color gluino-glue mass: 1.21(11)","feed_subtitle":"First non-perturbative SUSY-limit value; it fixes the large-N mass gap and tops SU(2), SU(3).","key_machinery":"The argument runs through large-$N$ twisted volume reduction: the twisted Eguchi-Kawai (TEK) model, a single-site matrix model with twisted boundary conditions, reproduces the infinite-volume, infinite-color lattice theory when center symmetry is preserved. Adjoint Wilson fermions are added through the TEK Wilson-Dirac operator, and the gluino-glue mass is extracted from the exponential decay of the correlator built from a clover-discretized field strength and the gluino propagator, with a generalized eigenvalue problem and stout smearing used to isolate the lightest state. Finite-$N$ effects in the TEK model act as finite-volume effects, so the $N=169,289,361$ data are extrapolated with an exponential finite-volume ansatz; the SUSY limit is then reached with a combined linear fit in $a/\\sqrt{8t_1}$ and $8t_1 m_\\pi^2$, justified by the absence of chiral logarithms in partially quenched chiral perturbation theory.","core_discovery":"The paper's central claim is that the mass of the lightest gluino-glue bound state in $\\mathcal{N}=1$ supersymmetric Yang-Mills theory at infinite number of colors is $w_0M_{\\tilde g g}=1.21(11)$ in the supersymmetric (chiral-continuum) limit, with $w_0$ the gradient-flow hadronic scale defined in Eq. (3.5). Expressed in the NSVZ scheme, this is $M_{\\tilde g g}/\\Lambda_{\\mathrm{NSVZ}}=8.94\\pm1.01$. The value comes from lattice Monte Carlo simulations of the twisted Eguchi-Kawai reduced model at $N=169,289,361$, extrapolated first to the thermodynamic limit and then to zero lattice spacing and zero adjoint-pion mass. The authors state that the large-$N$ result is close to but larger than the known $SU(2)$ and $SU(3)$ values, confirming the trend that the mass grows mildly with $N$ and is consistent with a naive $1/N^2$ extrapolation of the finite-$N$ data.","pith_inferences":["The obvious next check is to resolve the mixed-term coefficient $c_3$: with current errors $c_3=0.21(20)$ is consistent with zero, but the central value moves from 1.21 to 1.61 when it is included, so a higher-statistics determination of $c_3$ would tell whether the quoted SUSY-limit mass is stable.","The paper's proposed route to the other supermultiplet members is a spatially reduced lattice with temporal extent $N_t>1$; if such a simulation confirmed threefold degeneracy at large $N$, it would strengthen the supersymmetry-restoration picture and turn the gluino-glue mass into a genuine mass-gap prediction.","One could also test the $1/N^2$ trend with an independent standard-lattice $SU(4)$ or $SU(5)$ calculation; agreement with the fit would support large-$N$ scaling at surprisingly small color numbers, while disagreement would signal corrections beyond the naive leading term."],"forward_implications":["Under the standard supersymmetry argument that the gluino-glue state is degenerate with the lightest scalar and pseudoscalar states, the result fixes the mass gap of large-$N$ $\\mathcal{N}=1$ supersymmetric Yang-Mills at $M/\\Lambda_{\\mathrm{NSVZ}}=8.94\\pm1.01$.","The infinite-$N$ value $w_0M_{\\tilde g g}=1.21(11)$ is larger than the $N=2$ value $0.823(56)$ and the $N=3$ value $1.042(46)$; combined with those it supports a mild $1/N^2$ approach to the large-$N$ limit, with the fitted leading coefficient $1.212(71)$ matching the direct determination.","The result gives a concrete non-perturbative target that analytic, holographic, or semiclassical approaches to large-$N$ supersymmetric gauge theories must reproduce.","Because the same twisted-reduction setup already produced the large-$N$ gluino condensate and scale setting, the framework is now able to deliver renormalized large-$N$ SUSY observables independent of standard finite-volume lattice spectroscopy."],"supporting_citations":[{"why":"Supplies the gauge ensembles, the scale setting, and the adjoint-pion masses used for all extrapolations.","marker":"[27]"},{"why":"Provides the non-perturbative large-N value $\\sqrt{8t_1}\\Lambda_{\\mathrm{NSVZ}}=0.231(15)$ used to convert the gluino-glue mass into units of $\\Lambda$.","marker":"[28]"},{"why":"Establishes the twisted Eguchi-Kawai reduction as a matrix-model equivalent of large-N lattice gauge theory.","marker":"[12]"},{"why":"Supplies the symmetric-twist choice and the conditions needed for center symmetry to hold, justifying the simulation setup.","marker":"[17]"},{"why":"Derives the Wilson-Dirac operator for adjoint fermions in the TEK model used in the calculation.","marker":"[25]"},{"why":"Provides the partially quenched chiral perturbation theory result that there is no chiral logarithm in the adjoint-pion mass dependence, justifying the fit form.","marker":"[33]"},{"why":"Provides the SU(3) SUSY-limit value $w_0M_{\\tilde g g}=1.042(46)$ used for comparison.","marker":"[37]"},{"why":"Provides the SU(2) variational analysis result used for comparison.","marker":"[38]"},{"why":"Provides the ratio $w_0(c=0.1)/w_0(c=0.15)$ needed to convert the SU(2) result to the same scale definition.","marker":"[67]"}],"fun_headline_variants":["First non-perturbative large-N gluino-glue mass: 1.21(11)","Gluino-glue mass at N=∞: 1.21(11) from lattice","Twisted reduction gives infinite-color gluino-glue mass: 1.21(11)","Large-N gluino-glue mass matches SU(2), SU(3) trend: 1.21(11)","Non-perturbative large-N gluino-glue mass: 1.21(11) in NSVZ"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"Everything hinges on the assumption that the data approach the supersymmetric limit along a straight line in the lattice spacing $a$ and in the squared adjoint-pion mass $m_\\pi^2$, with no extra mixed correction; the paper itself reports that including such a mixed term changes the central value from 1.21 to 1.61, so the final number is only as solid as this assumption.","fun_headline_variants_meta":{"raw":{"variants":["First non-perturbative large-N gluino-glue mass: 1.21(11)","Gluino-glue mass at N=∞: 1.21(11) from lattice","Twisted reduction gives infinite-color gluino-glue mass: 1.21(11)","Large-N gluino-glue mass matches SU(2), SU(3) trend: 1.21(11)","Non-perturbative large-N gluino-glue mass: 1.21(11) in NSVZ"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.001268,"raw_usage":{"total_tokens":5157,"prompt_tokens":880,"completion_tokens":4277,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":496,"completion_tokens_details":{"reasoning_tokens":4145}},"tokens_in":496,"tokens_out":4277,"duration_ms":30638,"temperature":1.0,"reasoning_tokens":4145,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-11T23:34:38.548595+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Take the same $N=169,289,361$ ensembles and redo the chiral-continuum fit including the mixed term $c_3(a/\\sqrt{8t_1})(8t_1 m_\\pi^2)$ with higher statistics; the reported coefficient $c_3=0.21(20)$ is consistent with zero but shifts the central value from 1.21 to 1.61, so a determination of $c_3$ at the two-$\\sigma$ level would settle which SUSY-limit mass is right. Alternatively, data at a new, smaller lattice spacing would test whether the linear ansatz still describes the approach to the continuum limit.","supporting_citations":[{"cited_title":"Variational analysis of low-lying states in supersymmetric Yang-Mills theory","cited_arxiv_id":"1901.02416","evidence_quote":"Provides the SU(2) variational analysis result used for comparison."},{"cited_title":"Influence of topology on the scale setting","cited_arxiv_id":"1411.6995","evidence_quote":"Provides the ratio $w_0(c=0.1)/w_0(c=0.15)$ needed to convert the SU(2) result to the same scale definition."}],"review_version":1}