{"id":"4d9e3d2f-7012-4040-98ef-72726e6dc07f","arxiv_id":"2411.09650","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":4,"one_line_summary":"In a realistic general gauge mediation model, including gravitino decays of all sparticles weakens ATLAS mono-photon exclusion limits on the gluino-NLSP mass plane, especially for negative mu and compressed spectra.","lead":"This paper revisits ATLAS's 13 TeV mono-photon search for gauge-mediated supersymmetry and finds that the published limits are too strong in certain regions. In a more complete general gauge mediation model, the gluino mass limit drops to about 2.3 TeV instead of 2.4 TeV near the best-case neutralino mass, and the flat 2.2 TeV limit does not hold for large neutralino masses.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The quantitative limit shift rests on manually added gravitino decay widths for non-NLSP states, but the validation only checks the ATLAS scenario in which those widths are switched off; without an independent cross-check of the modified branching ratios, the central numerical claim is not yet…","rationale":"The central claim requires that non-NLSP gravitino decays can dominate in relevant parameter regions. The analytic widths are standard, but the implementation is not validated in the region where they matter. The reader identified this as the weakest assumption; I agree. The comparison with ATLAS in the ATLAS scenario is good support for the detector and recasting framework, but not for the modified decay tables. Because the result is quantitative—specific GeV limits and specific validity ranges—this is load-bearing. The fixed parameter slice (M2 = 3 TeV, tan beta = 1.5, decoupled sfermions, |M1| ~ |mu|) is a secondary generality limitation, but the central conclusion is explicitly about that slice, so it does not supersede the decay-table concern. This is not a disagreement with the physics motivation; it is a verification gap in the numerical engine behind the headline numbers. The paper remains CONDITIONAL pending independent validation or release of the modified decay tables.","tokens_in":19754,"tokens_out":5614,"duration_ms":61986,"concrete_test":"Take a benchmark grid covering the negative-mu and compressed regions (e.g., M_g ~ 2.0–2.5 TeV, M_chi10 ~ 1.3–1.8 TeV on the mu<0 lines of Fig. 2) and produce the SLHA decay blocks used in Section 3.3. Recompute the decays of g, chi_20, chi_30 and chi_1± with an independent implementation: SPheno cascade widths plus Eqs. 2.4–2.6 for gravitino widths, using the same mixing matrices, and compare every non-zero branching ratio to the submitted SLHA blocks. Then rerun the recasting for the benchmark points where any branching ratio differs by more than 5 percentage points and check whether the 95% CL exclusion boundary in Fig. 13 moves by more than ~50 GeV. If it does, the quantitative claim is not robust; if it does not, the concern is resolved.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The paper's quantitative conclusion—that the ATLAS limit drops from 2.4 to 2.3 TeV and that the 2.2 TeV bound holds only for restricted neutralino masses—is produced by allowing gravitino decays of the gluino, chi_20, chi_30 and chi_1+ in addition to cascade decays. In Section 3.3 the authors state that they 'modify the decay tables of the MSSM particles after calculating the gravitino partial widths using the analytical expressions provided in Section 2.1' (Eqs. 2.4–2.6). The cascade widths come from SPheno. The validation shown in the left panels of Figs. 7–12 compares acceptance×efficiency with ATLAS in the simplified scenario where, by construction, non-NLSP gravitino decays are absent. That validates the event generation, detector simulation and selection, but not the very branching ratios that are responsible for the right-panel suppression. The branching ratios BR(chi_20→X+G), BR(chi_1±→W+G) and BR(g→g+G) in the negative-mu and compressed regions are the difference between the left and right panels; if the SPheno three-body cascade widths are inaccurate in the compressed regime, or if the manual insertion changes total widths or branching-ratio normalization incorrectly, the regions with BR_gravitino > 50% (Figs. 3–4) and hence the new exclusion boundaries (Fig. 13) would shift. Since the claimed shift is only O(100) GeV, an unvalidated O(10–20%) BR error could alter or erase the conclusion. No code, SLHA grids, or modified decay tables are released.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"This paper reinterprets the ATLAS mono-photon + missing transverse energy search at 139 fb^-1 for simplified General Gauge Mediation, relaxing the assumption that only the lightest neutralino decays into a gravitino. The authors compute gravitino partial widths from analytic formulas in Eqs. (2.4)-(2.6), insert them into SPheno decay tables for the gluino, chi_20, chi_30, and chi_1+, and recast the ATLAS signal regions on the gluino-NLSP mass plane for both signs of mu. They find that in compressed spectra and especially for negative mu, non-NLSP gravitino decays suppress the mono-photon rate, so the ATLAS exclusion is locally weakened: the most stringent gluino mass limit changes from about 2.4 TeV to about 2.3 TeV for neutralino masses near 1.3-1.4 TeV, and the 2.2 TeV overall limit holds only in restricted neutralino-mass windows.","tokens_in":20052,"tokens_out":8509,"duration_ms":87141,"significance":"The paper identifies a physically plausible and useful effect: phase-space suppression of cascade decays in compressed GGM spectra allows direct gravitino decays of the gluino and of heavier electroweak-inos to compete, reducing the mono-photon yield on which the ATLAS limit is based. The analytic decay-width formulas and the decay phase diagrams in Figures 1-4 are a clear strength, and the simulation chain is standard and is validated against ATLAS acceptance x efficiency in the simplified scenario where non-NLSP gravitino decays are absent. The work is not circular: no model parameter is fitted to data, and the recast uses ATLAS background predictions and observed event counts. If the modified branching ratios are independently verified, the conclusion that the ATLAS simplified-GGM bounds are locally overestimated is a credible, modest correction to the experimental limits.","major_comments":[{"comment":"The quantitative conclusion rests on inserting gravitino partial widths for the gluino, chi_20, chi_30, and chi_1+ into the SPheno decay tables, but the validation shown in the left panels of Figs. 7-12 only tests the ATLAS scenario, in which those gravitino decay channels are absent by construction. Section 3.3 states that the default MSSM model file has no gravitino decay channels, so the acceptance x efficiency agreement in the simplified scenario validates event generation, detector simulation, and event selection, but does not test the modified branching ratios that produce the right-panel suppression and the new exclusion boundaries in Fig. 13. Because the reported shift is only about 0.1 TeV in the most stringent limit and up to about 0.3 TeV in parts of the plane, an unvalidated error of order 10-20% in BR(chi_20 -> X + G), BR(chi_1+ -> W + G), or BR(g -> g + G) could alter the claimed limits. Please add an independent cross-check of the modified decay tables, for example by comparing the SPheno-plus-manual-insertion branching ratios with an independent decay code or with the analytic partial widths on representative grid points, and verify the total-width normalization after insertion.","section":"Section 3.3, Eqs. (2.4)-(2.6), Figs. 7-12"},{"comment":"The scan is confined to the ATLAS benchmark slice M2 = 3 TeV, tan(beta) = 1.5, sfermions decoupled at 5 TeV, heavy Higgs states decoupled at 2 TeV, |M1| ~ |mu|, and m_gravitino ~ 1 eV. This is a model point inherited from the experimental simplified setup, not the full GGM parameter space, and the phrase \"realistic GGM scenario\" overstates the coverage of the scan. The competition between gravitino and cascade decays of chi_20, chi_30, and chi_1+ depends on the wino and higgsino admixtures set by M2, mu, and tan(beta); a sensitivity scan over at least M2 and tan(beta), even on a coarse grid, is needed to establish that the displaced exclusion boundaries in Fig. 13 are not artifacts of this benchmark slice. If the authors prefer to present only the ATLAS benchmark, the title and conclusions should be scoped accordingly.","section":"Section 2 and Section 3.3"}],"minor_comments":[{"comment":"The sentence \"the overall lower limit of 2.2 TeV only holds for M_chi10 < 1500 (1300) GeV and M_chi10 < 150 GeV for the positive (negative) mu scenario\" is internally confusing; the summary in Section 4 states instead that the limit holds only below 1500 GeV for positive mu and below 1300 GeV for negative mu. Please clarify the intended statement.","section":"Section 3.3.1"},{"comment":"The validation comparison with ATLAS is presented for only a handful of grid points, with the ATLAS values in brackets; a table listing all validation points and the relative differences in acceptance x efficiency would make the validation quantitative and reproducible.","section":"Figs. 7-12"},{"comment":"The statements that mu > 0 gives equal gamma and Z branching fractions while mu < 0 gives equal gamma and h branching fractions are not derived; a short expression for the relevant neutralino-mixing combinations would help readers reproduce the phase diagram in Fig. 2.","section":"Section 2.1"},{"comment":"No modified SLHA decay tables, spectrum files, or analysis scripts are released. For a recast paper of this type, providing at least representative SLHA files with the modified decay tables would greatly facilitate independent verification of the central step.","section":"Section 3.3"}],"recommendation":"major_revision","confidential_remarks":"The main concern is validation of the manual decay-table modification that drives the central result. If the authors can supply an independent branching-ratio cross-check and a sensitivity scan over at least M2 and tan(beta), the paper is publishable as a major revision; without such a check, the quantitative claims should be substantially qualified. I see no circularity or apparent misuse of the ATLAS data; the issue is reproducibility and validation of a key technical step."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Short version: this paper identifies a real, under-appreciated flaw in the ATLAS 139/fb mono-photon GGM interpretation—only the NLSP is allowed to decay to the gravitino—and shows where that assumption breaks down. The qualitative point is solid and the quantitative shift is plausible, though the headline numbers rest on a decay-table modification that is not independently validated.\n\nWhat's new: earlier GGM recasts don't address this specific assumption. The authors compute gravitino partial widths analytically, add them to SPheno decay tables, and map where non-NLSP gravitino decays dominate: compressed gluino-NLSP spectra for both signs of mu, and large parts of the negative-mu plane. The phase diagrams are clear and the physics is well explained—phase-space suppression of three-body cascade decays lets two-body gravitino decays win. The validation against ATLAS acceptance times efficiency in the simplified scenario is a legitimate check of the event generation and selection pipeline.\n\nSoft spots, in order. (1) The load-bearing ingredient is the manual insertion of gravitino widths into SPheno decay tables. The validation only exercises the ATLAS scenario, where those widths are switched off. So the agreement with ATLAS tells you nothing about whether the modified branching ratios are right where they matter. The formulas are standard and the implementation is probably fine, but 'probably' isn't the standard you want for a claimed 100–300 GeV shift in exclusion limits. No code, SLHA grids, or modified decay tables are released, so a referee can't check this directly. (2) The fixed slice (M2=3 TeV, tan beta=1.5, sfermions at 5 TeV) is taken from ATLAS, but the conclusion is framed as the 'realistic GGM' result. That overstates generality; other slices could give different quantitative limits, even if the qualitative effect persists. (3) Minor: the text says 'all possible decay modes' but the treatment covers the states in the gluino cascade; that's fine, but the wording overshoots.\n\nBottom line: the central physical argument holds up and the paper deserves a serious referee. The referee should require the decay-table code and a cross-check of the modified branching ratios with an independent tool, and should push back on the 'realistic GGM' framing. If those checks come back clean, the corrected limits are a useful contribution to SUSY recasting.","headline":"A solid recast paper with a real point: the ATLAS GGM mono-photon limit overstates exclusions where non-NLSP gravitino decays win, but the headline numbers need an independent check of the modified decay tables.","tokens_in":20662,"tokens_out":1978,"would_cite":true,"duration_ms":19763,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"The ATLAS mono-photon bound on GGM supersymmetry is weakened once direct gravitino decays of the gluino and heavier electroweakinos are included.","keywords":["Gauge-mediated supersymmetry breaking","General gauge mediation","gravitino LSP","mono-photon search","ATLAS constraints","gluino mass limits","neutralino NLSP","branching-ratio recast"],"falsifier":"Take a benchmark point in the disputed region, for example $m_{\\tilde g}=2.3$ TeV and $m_{\\tilde\\chi_1^0}=1.35$ TeV with $\\mu<0$, and compute the mono-photon signal yield with and without the added gravitino decays; if the suppression is smaller than the 95% CL uncertainty quoted by ATLAS, the claimed relaxation of the bound does not hold.","tokens_in":19487,"feed_emoji":"📉","tokens_out":13543,"duration_ms":111425,"temperature":0.7,"pith_summary":"Supersymmetry with gauge mediation predicts a very light gravitino, and collider searches usually assume that only the lightest neutralino decays into it. This paper argues that in a realistic General Gauge Mediation model, the gluino, the next two neutralinos, and the lightest chargino can also decay directly to a gravitino, especially when the spectrum is compressed or the Higgsino mass parameter $\\mu$ is negative. Because those direct decays bypass the lightest neutralino and often produce no energetic photon, the mono-photon final state used by the ATLAS search is depleted, so the published exclusion region is too aggressive. Recasting the $139~\\text{fb}^{-1}$ ATLAS search with all gravitino decays included lowers the most stringent gluino bound from about 2.4 TeV to about 2.3 TeV for negative $\\mu$ and shrinks the region where the overall 2.2 TeV bound applies. If correct, the paper changes what counts as excluded on the gluino-NLSP mass plane and points future searches toward boosted $W$, $Z$, Higgs, and top-quark signatures.","feed_headline":"ATLAS gluino limits shrink when all gravitino decays are included","feed_subtitle":"In compressed and negative-mu regions, extra gravitino decays cut the photon yield and relax the 2.4 TeV bound.","key_machinery":"The machinery is the set of gravitino partial widths for neutralinos, charginos, and the gluino, computed analytically in terms of the neutralino mixing matrix and masses. These widths are added to the spectrum generator's decay tables, where they compete with the standard cascade decays. The branching fractions are set by phase space: compressed spectra and negative $\\mu$ suppress the cascade channels, so the direct-to-gravitino modes dominate. This competition is what removes photons from the final state and relaxes the exclusion limits.","core_discovery":"The paper's central claim is that the assumption behind the ATLAS mono-photon search—that the gravitino LSP is reached only through the lightest neutralino—fails over a substantial part of the GGM parameter space. Using gravitino partial widths computed from the analytic formulas and inserted into the spectrum generator's decay tables, the authors find that for positive $\\mu$ the assumption mostly holds, while for negative $\\mu$ the second and third neutralinos and the lightest chargino decay to gravitinos with large branching fractions; in quasi-degenerate gluino-NLSP regions the gluino itself decays directly to a gravitino plus a gluon. The consequence is a photon-deficient final state. Recasting the ATLAS signal regions with these decays, the authors find that the published 2.4 TeV most-stringent gluino bound does not survive for negative $\\mu$: the updated limit is roughly 2.3 TeV for neutralino masses between 1.3 and 1.4 TeV, and the overall 2.2 TeV bound holds only for neutralino masses below 1500 GeV (positive $\\mu$) or 1300 GeV (negative $\\mu$), plus the below-150 GeV neutralino region, with discrepancies of up to 300 GeV compared with the original limit.","pith_inferences":["Beyond the paper's scan, the same phase-space competition should affect other GMSB searches—diphoton, multilepton, and long-lived particle searches—so simplified-model limits in those channels may also need an all-particles gravitino treatment.","The sign of $\\mu$ becomes a testable handle: a mono-photon deficit in the negative-$\\mu$ plane relative to the positive-$\\mu$ plane would support this mechanism over an overall cross-section suppression.","Quantitative shifts could change outside the paper's fixed parameter slice, but the qualitative conclusion—that ignoring non-NLSP gravitino decays overstates the excluded region—is likely robust to those choices.","A natural next step, not pursued here, is a recast of electroweakino pair production with the same modified decay tables, where direct gravitino decays of $\\tilde\\chi_2^0$ and $\\tilde\\chi_1^+$ should similarly weaken the limits."],"forward_implications":["For $\\mu < 0$, the published ATLAS most-stringent gluino mass limit of 2.4 TeV is replaced by roughly 2.3 TeV for neutralino masses between 1.3 and 1.4 TeV.","The overall gluino mass bound of 2.2 TeV survives only for neutralino masses below 1500 GeV (positive $\\mu$) or 1300 GeV (negative $\\mu$), plus the region below 150 GeV.","In quasi-degenerate gluino-NLSP spectra, the gluino's direct gravitino decay suppresses the mono-photon rate, so the compressed region is largely unconstrained by the ATLAS search.","For negative $\\mu$, the heavier neutralinos and the lightest chargino often decay straight to gravitinos, so the cascade never reaches the NLSP and no hard photon is produced.","In the regions where the bounds weaken, the final states contain highly boosted $W$, $Z$, Higgs, or top jets, so a fat-jet search could recover sensitivity."],"supporting_citations":[{"why":"The ATLAS mono-photon search at 139 fb^-1 whose simplified-GGM limits are the target of the reinterpretation.","marker":"[1]"},{"why":"Supplies the gravitino partial decay widths for neutralinos, charginos, and gluino used to modify the decay tables.","marker":"[37]"},{"why":"Defines General Gauge Mediation, the model framework whose mass plane is scanned.","marker":"[35]"},{"why":"Provides the three-body cascade decay modes of the heavier neutralinos and chargino that compete with the gravitino widths.","marker":"[84]"},{"why":"Gives the gluino cascade decay modes whose phase-space suppression drives the direct-gravitino dominance in compressed spectra.","marker":"[85]"},{"why":"The spectrum generator whose decay tables receive the manually added gravitino partial widths.","marker":"[93]"},{"why":"Provides the gluino pair production cross sections used to set the recast limits.","marker":"[95]"},{"why":"The statistical routine used to convert signal yields into 95% CL expected and observed exclusion limits.","marker":"[98]"}],"fun_headline_variants":["Gravitino decays cut ATLAS gluino limit by up to 300 GeV","Negative mu region weakens ATLAS mono-photon bounds","Direct gluino-to-gravitino decay weakens LHC reach","ATLAS gluino bound relaxes to 2.3 TeV with full decays"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The result stands or falls on the gravitino branching fractions for the gluino and the heavier neutralinos and charginos, which are computed by hand and added to the spectrum generator without being independently validated in the regime where they dominate.","fun_headline_variants_meta":{"raw":{"variants":["Gravitino decays cut ATLAS gluino limit by up to 300 GeV","Negative mu region weakens ATLAS mono-photon bounds","Direct gluino-to-gravitino decay weakens LHC reach","ATLAS gluino bound relaxes to 2.3 TeV with full decays"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000609,"raw_usage":{"total_tokens":2883,"prompt_tokens":1043,"completion_tokens":1840,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":659,"completion_tokens_details":{"reasoning_tokens":1758}},"tokens_in":659,"tokens_out":1840,"duration_ms":13996,"temperature":1.0,"reasoning_tokens":1758,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-12T20:25:22.277743+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Take a benchmark point in the disputed region, for example $m_{\\tilde g}=2.3$ TeV and $m_{\\tilde\\chi_1^0}=1.35$ TeV with $\\mu<0$, and compute the mono-photon signal yield with and without the added gravitino decays; if the suppression is smaller than the 95% CL uncertainty quoted by ATLAS, the claimed relaxation of the bound does not hold.","supporting_citations":[{"cited_title":"Chargino and Neutralino Decays Revisited","cited_arxiv_id":"hep-ph/0104115","evidence_quote":"Provides the three-body cascade decay modes of the heavier neutralinos and chargino that compete with the gravitino widths."},{"cited_title":"Barbieri, G","cited_arxiv_id":null,"evidence_quote":"Gives the gluino cascade decay modes whose phase-space suppression drives the direct-gravitino dominance in compressed spectra."},{"cited_title":"SUSYCrossSections","cited_arxiv_id":null,"evidence_quote":"Provides the gluino pair production cross sections used to set the recast limits."}],"review_version":1}