{"id":"f75bd2f8-f056-4b19-b6f5-f82080826825","arxiv_id":"1908.11160","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":4.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":7,"one_line_summary":"In an S4 x Zn flavor-symmetric SO(10) SUSY model with type-II seesaw, the mu-to-e-gamma branching ratio constrains CMSSM, NUHM and NUSM parameter space, with MEG-II projected to probe most remaining regions.","lead":"This paper calculates how often a muon should decay into an electron and a photon inside three supersymmetric models, using a specific family-symmetry grand unified theory framework. It maps which parts of each model's parameter space are already ruled out or could be seen by the future MEG-II experiment.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The central mu->e gamma boundaries are not independently checkable: Table I's delta inputs (Eqs. 5-7) are given without the f_nu matrix or intermediate calculation, and the delta_31 entry as rendered is unphysical.","rationale":"The paper's central contribution is not a new model but a set of phenomenological exclusion contours: the same S4 x Z_n type-II seesaw Dirac Yukawa structure is used to forecast MEG-II reach in three SUSY-breaking schemes. The one input that is both specific to the model and not shown is f_nu, and it enters all LFV rates multiplicatively. If f_nu differs - whether because the cited construction has free parameters that were not fixed in the text or because the table entries are mis-typed - the TeV-scale boundaries shift by large factors, since BR scales roughly as f_nu^4 / m_SUSY^8. The reader's weak-assumption note about the leading-log factor-of-10 caveat is real, but it is secondary: the paper says it uses SuSeFLAV's full two-loop RGE for the scans, so the leading-log equations are illustrations unless the paper explicitly substitutes them for the numerics. The unavailability of the f_nu matrix and the unphysical-looking delta_31 entry are what prevent an independent check of the headline numbers. I therefore keep the CONDITIONAL verdict: the qualitative ordering among CMSSM, NUHM and NUSM is credible, but the precise exclusion limits need to be re-derived from a documented f_nu input before they can be used as predictions.","tokens_in":16091,"tokens_out":8526,"duration_ms":80661,"concrete_test":"Take the S4 x Z_n construction of [28], write out f_nu explicitly, and recompute (f_nu*)_{13}(f_nu)_{23} and (f_nu*)_{23}(f_nu)_{33} with M_R1 = 1e13 GeV, M_R2 = 1e14 GeV, M_R3 = 1e16 GeV, M_X = 2e16 GeV; compare to the delta_12 and delta_23 rows of Table I. Then rerun SuSeFLAV at one representative CMSSM point on the claimed exclusion edge, e.g. m0 = 6 TeV, M_1/2 = 3 TeV, tan beta = 5, A0 = 0, and verify that the resulting BR(mu -> e gamma) matches the plotted contour. If the reconstructed delta values or the BR differ by more than a factor of about 2, the numerical boundaries in the strongest claim do not follow from the model as presented.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The strongest claim's quantitative boundaries all pass through the delta_ij inputs of Table I, which are supposed to enter Eqs. (5)-(7). The paper never gives the S4 x Z_n Dirac Yukawa matrix f_nu from [28] nor the intermediate calculation that yields delta_12 = 0.6672e-4, delta_23 = 1.5634e-4, and delta_31 = 0.7377e43. As rendered, the last entry is unphysical, so at least one table value is not a reliable input. For the mu->e channel that drives the headline, Eq. (5) depends on (f_nu*)_{13}(f_nu)_{23} ln(M_X/M_R3), and no numerical value for this product is supplied at the stated M_R1 = 1e13 GeV, M_R2 = 1e14 GeV, M_R3 = 1e16 GeV. Without these numbers a reader cannot tell whether the TeV boundaries in Figs. 3-8 reflect the S4 x Z_n model or just the internal choices of the SuSeFLAV scan. The paper's own caveat after Eq. (13) - that at M_1/2 ~ 1 TeV the leading-log BR can differ from full RGE running by up to a factor of 10 - also applies to the NUSM lower boundary, and the paper does not document which plots use the leading-log formulas and which use the claimed full two-loop running. The qualitative hierarchy (CMSSM heavy, NUHM light via cancellations, NUSM wide) is plausible and consistent with prior literature, but the specific mass limits in the central claim are not independently verifiable from the text.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"This paper studies the charged-lepton-flavor-violating decay mu -> e gamma in supersymmetric SO(10) models with an S4 x Z_n flavor symmetry and type-II seesaw neutrino masses. Using the Dirac neutrino Yukawa texture of Ref. [28], the author scans the soft SUSY-breaking parameter space of CMSSM, NUHM, and NUSM with the public code SuSeFLAV and applies the MEG 2016 bound BR(mu -> e gamma) < 4.2 x 10^-13 and the projected MEG-II sensitivity 6 x 10^-14. The main results are allowed regions in (m0, M1/2, tan beta, A0) for each model: CMSSM survives only for heavy spectra, NUHM allows much lighter spectra because of cancellations involving m_Hu, and NUSM leaves a wide range of M1/2 up to about 6 TeV. The paper also comments on the reach of HE/HL-LHC for the surviving spectra.","tokens_in":16470,"tokens_out":9087,"duration_ms":89004,"significance":"If the numerical results were fully supported, the paper would provide a useful model-discrimination statement for MEG-II and HL-LHC: the three SUSY boundary conditions produce qualitatively different allowed regions, and the projected MEG-II sensitivity would probe essentially all of NUHM and much of NUSM while leaving CMSSM largely untouched. The use of the public SuSeFLAV package and of an externally published S4 x Z_n Yukawa texture are appropriate, and the comparison with the current MEG limit is concrete. However, the quantitative boundaries are not currently reproducible from the text: the key delta_ij inputs are not derived, one table entry is unphysical, and the relationship between the leading-log equations and the claimed full two-loop running is not resolved. The qualitative hierarchy is plausible and consistent with the earlier literature, but the specific mass limits should be treated with caution until the inputs are documented.","major_comments":[{"comment":"The numerical inputs that drive the LFV predictions are stated without derivation: the Dirac neutrino Yukawa matrix f_nu from [28] is not displayed, not even the (1,3), (2,3), and (3,3) combinations that enter Eqs. (5)-(7), and no intermediate calculation is shown that yields delta_12 = 0.6672 x 10^-4, delta_23 = 1.5634 x 10^-4, or delta_31 = 0.7377 x 10^43. The last value is impossible as a dimensionless mass insertion and indicates at least a typographical error. Because delta_ij is defined in Sec. II.A as Delta_ij / m_tilde_l^2, fixed table values also presuppose a fixed slepton scale that is not specified. Please provide f_nu or a step-by-step derivation, correct the table, and state how the table entries are used in the SuSeFLAV runs. This is load-bearing because the excluded/allowed boundaries in Figs. 3-8 ultimately rest on these numbers.","section":"Sec. III, Table I and Eqs. (5)-(7)"},{"comment":"The paper's own caveat after Eq. (13) states that for M1/2 around 1 TeV the branching fractions can differ by up to a factor of 10 from the full RGE running, while the abstract and Sec. IV claim that the numerical analysis includes full two-loop RGE running. Since the scans extend to M1/2 of 4.5-6 TeV and the NUSM lower boundary sits at M1/2 around 1 TeV, which is exactly where the stated factor-of-10 discrepancy applies, the text must state explicitly which figures are produced by SuSeFLAV's full running and which by the leading-log equations (3) and (13). Without this specification, a factor-of-10 uncertainty applies to the lower boundaries and the sharp quantitative limits quoted in Sec. IV are not supported as stated.","section":"Sec. II.A (after Eq. (13)) and Sec. IV"},{"comment":"Several central quantitative statements are read from scatter plots without documented acceptance or rejection criteria or scan density; examples include Sec. IV.A's statement that m0 lies between 4.5 TeV and 8 TeV and the ranges compiled in Tables II and III. In addition, the column headers of Tables II and III are inconsistent, with Table III's first column labeled 'CMSSM' although the section and table title concern NUSM. Please provide the SuSeFLAV input files or benchmark points, and state the number of scan points and the criterion for 'allowed'. Without this, the printed boundaries are not reproducible and the reader cannot tell whether they reflect the S4 x Z_n texture or internal choices of the scan.","section":"Sec. IV and Tables II-III"},{"comment":"The abstract states that regions excluded by LHC searches are specified, but the results sections apply only the MEG bound, the Higgs-mass window, and projected sensitivities; no explicit LHC sparticle-search exclusion, such as gluino or squark mass limits, is referenced or overlaid in the figures. Either add the LHC constraints actually used, or soften the claim in the abstract and conclusion to match what is presented.","section":"Abstract and Sec. IV"}],"minor_comments":[{"comment":"The title contains a typo, 'L ight' for 'Light', and the abstract has 'for the the above mentioned'; these should be corrected.","section":"Title and Abstract"},{"comment":"There are repeated language and typographical errors, including 'Feynmann' for 'Feynman', 'paramater' for 'parameter', and inconsistent formatting such as 'SuSeFL A V'; a careful editorial pass is needed.","section":"Throughout"},{"comment":"The sentence 'the parameter space M1/2 >= 10 GeV is permitted by present MEG bounds' is clearly missing a factor of 10^3 and should read at least 1 TeV, consistent with the surrounding discussion and Table II.","section":"Sec. IV.A"},{"comment":"After Eq. (9), the neutrino masses are listed as 'm_nu3 = 0.05 eV, m_nu3 = 0.01 eV, and m_nu1 = 0.005 eV'; the second occurrence of m_nu3 should be m_nu2.","section":"Sec. II.B"},{"comment":"References [29] and [39] appear to be the same paper and are redundant; please check and consolidate them.","section":"References"}],"recommendation":"major_revision","confidential_remarks":"The manuscript needs substantial language editing and a proper derivation or appendix for the delta_ij inputs. The self-citation [31] is contextual rather than a circularity problem; the main concern is technical reproducibility, not citation ethics."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Kavita,\n\nQuick take: this is a parameter scan, not an analytic derivation. The paper updates an earlier type-I seesaw LFV analysis to a type-II seesaw in the S4 x Zn SO(10) model, and it compares three SUSY-breaking schemes against the current MEG bound and MEG-II projection. The qualitative punchline—CMSSM is pushed to heavy spectra, NUHM can stay light because of cancellations among m0, A0, and mHu, and NUSM has a wide allowed region—is plausible and consistent with older literature. That part is worth something.\n\nWhat is genuinely new is modest: the texture is imported from [28], the SUSY models are standard, and the formulas are the usual leading-log mass insertions. But applying them to the updated MEG limit and the MEG-II reach in this specific SO(10) flavor model is a legitimate exercise, and the direct comparison across boundary conditions is useful for the LFV phenomenology community.\n\nThe soft spots are real. Table 1 lists delta31 as 0.7377 x 10^43, which is unphysical; that is almost certainly a typo, but it is not the only problem. The delta12 and delta23 values are given without the f_nu matrix or any intermediate step, so a reader cannot reproduce the numbers that feed Eqs. (5)-(7). The paper says the analysis includes 'full 2 loop RG running' in the abstract, but the LFV rates are computed with leading-log approximations, and the caveat after Eq. (13)—that at M1/2 ~ 1 TeV the branching fraction can differ from full RGE by up to a factor of 10—applies to most of the plotted range since M1/2 goes up to 4.5-6 TeV. No code or input files are provided. There are also copy-paste slips in Tables II and III. None of this kills the qualitative conclusion, but it does mean the specific mass boundaries in the figures should not be taken at face value.\n\nWho is this for? A phenomenologist working on LFV in SO(10) or scanning SUSY parameter space might cite it after the table issues are fixed. It deserves a serious referee, but one who asks the author to show the f_nu computation, correct the typo, and clarify which regions use the leading-log approximation rather than the full running.\n\nMy recommendation: send it back for a major revision, not a desk rejection. The kernel is sound; the presentation is not.","headline":"A useful but under-verified parameter scan: the qualitative CMSSM/NUHM/NUSM hierarchy under MEG-II is credible, but the paper's own leading-log caveat and an unphysical table entry keep the quantitative boundaries from being trusted.","tokens_in":17084,"tokens_out":2913,"would_cite":false,"duration_ms":26515,"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":"In the $S_4 \\times Z_n$ type-II seesaw SO(10) framework, the $\\mu \\to e \\gamma$ rate is computed to sharply separate CMSSM, NUHM, and NUSM supersymmetry-breaking patterns, with MEG-II projected to probe nearly all of NUHM and much of NUSM…","keywords":["lepton flavor violation","mu to e gamma","supersymmetry","SO(10) grand unified theory","S4 flavor symmetry","type-II seesaw","non-universal Higgs model","MEG experiment"],"falsifier":"Compute $\\text{BR}(\\mu\\to e\\gamma)$ by full two-loop RGE running of the $S_4\\times Z_n$ Dirac neutrino Yukawa matrix at representative points (for instance $m_0 = 8$ TeV, $M_{1/2} = 3$ TeV, $\\tan\\beta = 5$) and compare with the leading-log formula; an order-of-magnitude discrepancy would shift the claimed allowed regions. Alternatively, a null result from MEG-II would exclude every predicted NUHM point whose central rate exceeds $6\\times 10^{-14}$, directly contradicting the paper's claim that almost all of NUHM lies within MEG-II reach.","tokens_in":15832,"feed_emoji":"⚛️","tokens_out":10592,"duration_ms":86699,"temperature":0.7,"pith_summary":"The paper argues that in a supersymmetric SO(10) grand unified theory with an $S_4 \\times Z_n$ flavor symmetry and a type-II seesaw for neutrino masses, the decay $\\mu \\to e \\gamma$ becomes a sharp discriminator among three supersymmetry-breaking patterns. It claims that the current MEG upper bound $\\text{BR}(\\mu\\to e\\gamma) < 4.2 \\times 10^{-13}$ and the projected MEG-II sensitivity of $6 \\times 10^{-14}$ leave almost disjoint surviving regions: CMSSM must have very heavy scalar and gaugino masses, NUHM can be light only because of cancellations in the off-diagonal slepton-mass entries, and NUSM allows a wide range of masses. If these predictions are right, the next round of muon lepton-flavor-violation experiments would effectively survey the NUHM and NUSM parameter spaces while leaving CMSSM mostly untouched, offering a way to distinguish the models at the HE/HL-LHC.","feed_headline":"MEG-II would probe nearly all NUHM and much of NUSM SUSY space","feed_subtitle":"In this seesaw SO(10) model, current MEG bounds force CMSSM to multi-TeV masses, leaving light NUHM and NUSM regions for MEG-II.","key_machinery":"The central machinery is the leading-log mass-insertion approximation for the off-diagonal left-handed slepton mass matrix. In the $S_4\\times Z_n$ type-II seesaw SO(10) model, the Dirac neutrino Yukawa matrix $f_\\nu$ is fixed by the flavor symmetry, and the entries $(\\delta_{LL})_{ij}$ are proportional to $(f_\\nu^\\dagger)_{ik}(f_\\nu)_{jk}\\log(M_X/M_{R_k})$. For CMSSM the prefactor is $(-3m_0^2 + A_0^2)/(8\\pi^2)$, while for NUHM it becomes $(-2m_0^2 + A_0^2 + m_{H_u}^2)/(8\\pi^2)$, and the relative sign between $m_{H_u}^2$ and $m_0^2$ is what permits cancellations. The numerical scans use the paper's chosen spectrum and LFV computation code, with full two-loop RGE running of the Yukawa couplings.","core_discovery":"Using the $S_4 \\times Z_n$ constrained type-II seesaw framework, the paper computes the lepton-flavor-violating mass insertions $(\\delta_{LL})_{ij}$ from the Dirac neutrino Yukawa couplings at the GUT scale and evaluates $\\text{BR}(\\mu\\to e\\gamma)$ under CMSSM, NUHM, and NUSM boundary conditions. The central finding is that the current bound of $4.2\\times 10^{-13}$ and the future reach of $6\\times 10^{-14}$ carve out qualitatively different allowed regions: CMSSM requires $m_0$ roughly 4.5--8 TeV with $M_{1/2}$ above about 3 TeV and a narrow $\\tan\\beta$ band; NUHM permits spectra as light as about 1 TeV in $M_{1/2}$ because negative $A_0$ and the Higgs soft mass $m_{H_u}$ partially cancel against $m_0^2$ in the off-diagonal slepton mass; and NUSM with multi-TeV first-two-generation scalars survives over $M_{1/2}$ from about 1 to 6 TeV and $m_0$ up to 16 TeV, with $\\tan\\beta$ restricted to 5--47 and, for $m_h \\simeq 125.9$ GeV, to 15--30. The paper reads these differing fates as a way to distinguish the three supersymmetry-breaking patterns at MEG-II and the HE/HL-LHC.","pith_inferences":["If the paper's admitted factor-of-10 uncertainty in the leading-log formula propagates up to $M_{1/2}$ of several TeV, the true exclusion boundaries could shift by thousands of GeV; a dedicated full-RGE benchmark scan on the $S_4\\times Z_n$ $f_\\nu$ matrix would settle this.","Because the $S_4\\times Z_n$ model fixes the relative sizes of $\\delta_{12}$, $\\delta_{23}$, and $\\delta_{31}$, a future measurement of $\\tau\\to\\mu\\gamma$ and $\\tau\\to e\\gamma$ alongside $\\mu\\to e\\gamma$ would test the flavor-symmetry structure itself, not just the individual rates.","The paper does not derive the tabulated $\\delta$ values from the model's vacuum alignments; deriving them would turn the phenomenological scan into a first-principles test of the flavor symmetry, and any inconsistency would point to corrections to the leading-log or type-II seesaw assumptions."],"forward_implications":["MEG-II at $6\\times 10^{-14}$ would probe essentially all of the NUHM parameter space that survives the 2016 MEG bound, because the cancellation mechanism keeps rates above the future sensitivity.","In CMSSM, only very heavy spectra ($m_0 \\sim 4.5$--$8$ TeV, $M_{1/2} \\gtrsim 3$ TeV) remain, so a null MEG-II result would not further constrain CMSSM, while a positive signal would strongly disfavor it.","In NUSM, MEG-II would restrict $\\tan\\beta$ to below about 20, and the surviving points predict low LFV rates, making MEG-II and HE/HL-LHC complementary probes.","The model-by-model allowed regions listed in the summary tables provide direct target lists for HE/HL-LHC sparticle searches, since each model corresponds to a distinct $m_0$--$M_{1/2}$--$A_0$ pattern."],"supporting_citations":[{"why":"Supplies the Dirac neutrino Yukawa matrix $f_\\nu$ from the $S_4 \\times Z_n$ type-II seesaw model, the input that fixes the flavor-violating mass insertions.","marker":"[28]"},{"why":"Sets the current experimental upper bound $\\text{BR}(\\mu\\to e\\gamma) < 4.2 \\times 10^{-13}$ used to define allowed regions.","marker":"[16]"},{"why":"Defines the projected MEG-II sensitivity of $6 \\times 10^{-14}$ that carves out the future probe regions.","marker":"[17]"},{"why":"Provides the SUSY spectrum and the LFV branching-ratio computation used for the scans.","marker":"[50]"},{"why":"Introduces the leading-log mass-insertion framework for SUSY SO(10) seesaw and the earlier type-I seesaw comparison.","marker":"[11]"},{"why":"Defines the non-universal scalar mass model (NUSM) whose boundary conditions are scanned.","marker":"[39]"},{"why":"Defines the non-universal Higgs mass (NUHM) parametrization used in the NUHM scans.","marker":"[40]"}],"fun_headline_variants":["MEG-II to distinguish SUSY breaking via mu->e gamma","mu->e gamma could tell CMSSM, NUHM, NUSM apart","S4 seesaw: MEG-II probes NUHM, NUSM, not CMSSM","LFV in S4 SO(10): MEG-II to test non-universal SUSY"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that the leading-log mass-insertion formulas with the tabulated $\\delta$ values give the right branching fractions, even though the paper concedes that for $M_{1/2}$ around 1 TeV the result may differ from a full RGE evaluation by up to a factor of 10.","fun_headline_variants_meta":{"raw":{"variants":["MEG-II to distinguish SUSY breaking via mu->e gamma","mu->e gamma could tell CMSSM, NUHM, NUSM apart","S4 seesaw: MEG-II probes NUHM, NUSM, not CMSSM","LFV in S4 SO(10): MEG-II to test non-universal SUSY"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000312,"raw_usage":{"total_tokens":1876,"prompt_tokens":1148,"completion_tokens":728,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":764,"completion_tokens_details":{"reasoning_tokens":635}},"tokens_in":764,"tokens_out":728,"duration_ms":7276,"temperature":1.0,"reasoning_tokens":635,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T10:22:39.667770+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Compute $\\text{BR}(\\mu\\to e\\gamma)$ by full two-loop RGE running of the $S_4\\times Z_n$ Dirac neutrino Yukawa matrix at representative points (for instance $m_0 = 8$ TeV, $M_{1/2} = 3$ TeV, $\\tan\\beta = 5$) and compare with the leading-log formula; an order-of-magnitude discrepancy would shift the claimed allowed regions. Alternatively, a null result from MEG-II would exclude every predicted NUHM point whose central rate exceeds $6\\times 10^{-14}$, directly contradicting the paper's claim that almost all of NUHM lies within MEG-II reach.","supporting_citations":[{"cited_title":"The Scale-Invariant NMSSM and the 126 GeV Higgs Boson","cited_arxiv_id":"1212.5243","evidence_quote":"Provides the SUSY spectrum and the LFV branching-ratio computation used for the scans."},{"cited_title":"Witten, Nucl","cited_arxiv_id":null,"evidence_quote":"Introduces the leading-log mass-insertion framework for SUSY SO(10) seesaw and the earlier type-I seesaw comparison."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Defines the non-universal scalar mass model (NUSM) whose boundary conditions are scanned."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Defines the non-universal Higgs mass (NUHM) parametrization used in the NUHM scans."}],"review_version":1}