{"id":"b28c3a88-b10e-4dd0-9083-a1047fd18fcd","arxiv_id":"1909.02029","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":8.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":4,"one_line_summary":"The complete leading-order two-loop Majoron couplings to all Standard Model gauge bosons and to flavor-changing quarks are derived, and they change the expected photon coupling and rare-decay rates compared with earlier one-loop estimates.","lead":"This paper computes, for the first time, the two-loop couplings of the Majoron (a hypothetical particle associated with neutrino mass) to photons, gluons, W/Z bosons, and flavor-changing quarks. These couplings determine the main observable signals for Majoron dark matter and rare decays, so the calculation sharpens experimental searches in both muon and tau sectors.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Central completeness claim is contradicted by Sec. III F: the two-loop J–hZ and J–hγ couplings are acknowledged to exist but are not computed, so 'all Majoron couplings to SM particles' is not established.","rationale":"The reader's weakest_assumption focused on the free-quark treatment for low Majoron masses, which is real but acknowledged and primarily affects phenomenological limits, not the existence of the central two-loop calculation. The more load-bearing issue is the paper's own admission in Sec. III F that J–hZ and J–hγ couplings arise at two loops but are not presented. Since the abstract and conclusion claim completeness, this is an internal inconsistency in the central claim rather than a numerical approximation. I do not challenge the correctness of the presented formulas: the calculation is detailed, cross-checked against limits, and published. The concern is that the scope of the claim exceeds the content. The reader's CONDITIONAL verdict already accommodates this issue, so no verdict change is needed; the manuscript should either include the missing amplitudes or explicitly restrict the completeness claim to the couplings actually presented.","tokens_in":21657,"tokens_out":7237,"duration_ms":86327,"concrete_test":"Re-run the authors' FeynArts/FeynRules implementation with the expansion-by-regions code used in Sec. III to evaluate the two-loop J→hZ and J→hγ amplitudes at leading order in the seesaw expansion, with m_J small and m_h fixed. If either amplitude is nonzero at O(1/f), the claim to have derived 'all' two-loop Majoron couplings to SM particles is false as stated; the paper would need to present those amplitudes or explicitly restrict its claim. If both amplitudes vanish at this order, the completeness claim can be restored by adding a proof that they vanish.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The paper's central claim—'complete the program ... derive all Majoron couplings to SM particles' (Sec. I and the Conclusion)—is not supported by the paper's own Sec. III F. There the authors state that CP invariance allows couplings of J to hZ and hγ, that the dominant contributions to both couplings arise at two-loop order, and that 'we will not, however, present the results here.' These are couplings to SM particles at the same claimed order in the seesaw expansion. Unless those amplitudes vanish identically at leading order—which the text does not show—the completeness claim is overbroad. The acknowledged hadronic replacement for m_J around the MeV scale (Fig. 6) is a numerical-uncertainty issue in derived constraints, but the missing hZ/hγ amplitudes are an internal inconsistency in the headline claim. The two-loop formulas that are presented may well be correct; the problem is the scope of the claim.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper studies the minimal singlet Majoron model with spontaneous lepton-number breaking. In the seesaw limit MD << MR, it computes the two-loop couplings of the Majoron to gluons, photons, Z-gamma, ZZ, WW, and to flavor-changing quark pairs dd' and uu' at leading order in the seesaw expansion. The couplings are given as explicit expressions in terms of one-loop integral functions; the authors verify UV finiteness, independence of the gamma5 treatment, and the expected low-energy and soft limits. They then use these couplings to discuss constraints from rare kaon/B decays, supernova cooling, lepton-flavor-violating decays, and Majoron dark matter, and compare with standard seesaw observables.","tokens_in":21812,"tokens_out":7991,"duration_ms":84750,"significance":"The calculation is technically demanding, involving on the order of one hundred two-loop diagrams, and is performed with multiple cross-checks, including UV-finiteness, gamma5-scheme independence, and agreement with known limits. The main physics findings—an mJ^2-suppressed and potentially cancellation-prone J-gamma-gamma coupling, non-vanishing J-Z-gamma, J-ZZ, and J-WW couplings, and log(MR/mW)-enhanced flavor-changing quark couplings—are novel and phenomenologically relevant. The explicit formulas constitute a useful resource for experimental searches. If the missing hZ/h-gamma amplitudes are supplied or the scope claim is revised, this would be a valuable contribution appropriate for the journal.","major_comments":[{"comment":"The paper's headline claim of deriving 'all Majoron couplings to SM particles' is not supported by its own text: Sec. III F states that CP invariance permits couplings of the Majoron to hZ and h-gamma, that the dominant contributions to both arise at two-loop order, and that 'we will not, however, present the results here.' Unless these couplings vanish identically at leading order in the seesaw expansion—which is not shown—the completeness claim is false as stated. The authors should either compute these amplitudes or explicitly revise the abstract, introduction, and conclusion to state that the paper covers the phenomenologically important gauge-boson and flavor-changing-quark couplings, listing hZ and h-gamma as omitted.","section":"Sec. III F (with Sec. I and Sec. V)"},{"comment":"The low-mass phenomenological constraints (mJ below roughly 1 GeV) rely on the free-quark expression for gJ-gamma-gamma with the ad hoc replacements mu=md=m_pi and ms=m_K, a procedure the text itself labels 'very naive' and cautions 'should not be taken too seriously.' Because Fig. 6 is presented as the central constraint summary, the exclusion regions between about 1 MeV and 1 GeV should be clearly marked as model-dependent estimates, and the text should state that robust quantitative bounds in this mass range require a hadronic treatment rather than free-quark thresholds.","section":"Sec. IV A and Fig. 6"},{"comment":"The subleading contribution Lsub_Jdd' in Eq. (37) is neglected in the subsequent phenomenological analysis with the justification that log(MR/mW) >~ 1. Since log(MR/mW) is only O(1) when MR is near the TeV scale, and no quantitative comparison between the leading-log and subleading terms is provided, the dominance of the leading-log coupling in the K -> pi J and B -> K J constraints should be demonstrated, or Lsub should be included in the numerical analysis.","section":"Sec. III G"}],"minor_comments":[{"comment":"The rendering of the Z-boson coupling term appears garbled: 'ni/Z' seems to be missing a gamma matrix (presumably gamma^mu Z_mu or a slashed Z). Please check the typeset version.","section":"Eq. (3)"},{"comment":"The caption describes the blue dashed line as (MDMD†)ee=(MDMD†)mu-mu=(100 GeV)^2, (MDMD†)tau-tau=0 and simultaneously states that this configuration corresponds to gJee=gJmu-mu=0. These statements are inconsistent with the definitions in Eqs. (10) and (11), since nonzero diagonal entries of MDMD† do not make the corresponding diagonal fermion couplings vanish. The caption should be corrected.","section":"Fig. 4 caption"},{"comment":"The paper relies on an in-house Mathematica implementation of expansion by regions and two-loop tensor reduction, but no code or detailed algebraic steps are provided. An ancillary file or an appendix with additional intermediate steps would substantially improve reproducibility.","section":"Sec. III"},{"comment":"The one-argument loop functions h(x) and g(x) are defined in Eqs. (16) and (27), while a two-argument function g(x,y) is introduced later in Eq. (34). Collecting all loop-function definitions in a single table or appendix would improve readability.","section":"Sec. III B and III E"}],"recommendation":"major_revision","confidential_remarks":"The technical content appears sound and the calculations are carefully validated. The main issue is the overbroad completeness claim, which conflicts with the explicit omission of the hZ and h-gamma two-loop amplitudes. This is fixable either by computing those amplitudes or by narrowing the claim, so I do not see grounds for rejection. The hadronic limitation is honestly acknowledged but should be more prominently reflected in the presentation of the low-mass constraints."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The important physics here is the two-loop, leading-seesaw Majoron couplings to gauge bosons and flavor-changing quarks. The genuinely new piece is the Set I 1PI contribution to Jγγ and the full set of JZγ, JZZ, JWW, and Jgg amplitudes, plus the log-enhanced Jdd′ and Juu′ couplings. Before this paper only the J–Z mixing (Set II) piece of Jγγ had been estimated, so this fills a real gap in the 40-year-old Majoron program.\n\nThe calculation is careful and transparent. The authors give explicit formulas, check UV finiteness using the identities of Eq. (5), verify γ5-scheme independence, and confirm the expected low-mass and soft limits. The FeynArts/FeynRules pipeline is validated on known tree- and one-loop results. The physics points, like the mJ^2 suppression of Jgg and the possible cancellations in Jγγ from the two diagram sets, are credible and interesting. The phenomenology section is honest about its approximations.\n\nThe soft spots are real but not damaging. The abstract and introduction claim to derive “all Majoron couplings to SM particles,” yet Sec. III F explicitly states that the two-loop J–hZ and J–hγ couplings exist but are not presented. That is a load-bearing claim in the framing only; the authors acknowledge the omission and give a reasonable justification (low phenomenological relevance). Still, “complete” is overbroad and should be qualified. Second, the low-mass hadronic treatment is acknowledged to be crude: setting mu=md=mπ and ms=mK in the Jγγ loop is a naive substitution, and the authors themselves say the resulting limit “should not be taken too seriously.” That weakens the low-mass exclusion plot but does not touch the two-loop formulas in the free-quark regime.\n\nThe citation pattern is appropriate, including the reuse of the Set II formula from the authors’ earlier work, which is properly credited. I have not re-run the integrals, but the internal consistency checks and the published record give me reasonable confidence in the central results.\n\nThis paper is for anyone working on Majoron models, axion-like particles, rare B and K decays, or Majoron dark matter. It deserves serious peer review. The main revision I would ask for is to temper the completeness claim and make clear that hZ/hγ are not computed here. I would bring it to a reading group and would cite it in related work.","headline":"Solid two-loop calculation that completes most of the Majoron coupling program; the completeness claim is slightly overbroad because the hZ/hγ couplings are omitted, but that is a framing issue, not a physics flaw.","tokens_in":22421,"tokens_out":1457,"would_cite":true,"duration_ms":17534,"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 paper claims that the Majoron's two-loop couplings to gauge bosons and flavor-changing quarks complete the leading-order description of how the Majoron talks to every Standard Model particle.","keywords":["Majoron","two-loop couplings","seesaw mechanism","lepton number violation","axion-like particle","flavor-changing neutral currents","Majoron dark matter","rare kaon decays"],"falsifier":"An independent numerical evaluation of the two-loop amplitudes at a benchmark point (say $m_J=1\\,\\mathrm{GeV}$, $f=M_R=1\\,\\mathrm{TeV}$, diagonal $M_DM_D^\\dagger=(100\\,\\mathrm{GeV})^2$) would settle the formal calculation, since the paper's expressions are explicit enough to compare numerically. For the phenomenological layer, a hadronic computation of $J\\to\\gamma\\gamma$ between 1 and 100 MeV would decide whether the quoted exclusion curves survive.","tokens_in":21395,"feed_emoji":"⚛️","tokens_out":12259,"duration_ms":126411,"temperature":0.7,"pith_summary":"The paper completes a program that began with the original 1981 Majoron proposal: it derives the couplings of the Majoron, the Goldstone boson of spontaneously broken lepton number, to every Standard Model particle, at leading order in the seesaw expansion. The missing pieces were the two-loop couplings to gauge bosons and to quarks of different generations. If the calculation is right, the Majoron is a fully predictive axion-like particle whose couplings are fixed by the seesaw parameters, and rare processes such as $\\tau\\to\\ell J$, $K\\to\\pi J$, and $B\\to K J$ become direct probes of those parameters. The result also reshapes Majoron dark matter expectations, because the two-loop diphoton coupling can be suppressed by cancellations, leaving neutrinos as the dominant decay signal.","feed_headline":"Majoron two-loop couplings expose rare decays as seesaw probes","feed_subtitle":"Known Majoron couplings make kaon and tau rare decays direct probes of the seesaw.","key_machinery":"The load-bearing object is the two-loop effective-vertex calculation in the minimal singlet Majoron model, organized by the seesaw expansion in the small ratio $v/f$ of electroweak scale to lepton-number-breaking scale. The diagrams split into Set I (genuine two-loop 1PI diagrams) and Set II (reducible diagrams dominated by $J$-$Z$ mixing), and the computation uses standard identities among the seesaw mixing matrices $C$, $B$, and $M_n$ to ensure the amplitudes are ultraviolet finite. The quoted outputs are coefficients $g_{JVV'}$ of on-shell decay amplitudes $M(J \\to VV') = -g_{JVV'} \\varepsilon^{\\mu\\nu\\rho\\sigma} \\varepsilon^*_\\mu(k_1) \\varepsilon^*_\\nu(k_2) k_{1\\rho} k_{2\\sigma}$, with the loop functions $h(x)$, $g(x)$, and the Passarino-Veltman $C_0$ carrying all mass dependence; after reduction, no two-loop master integrals survive.","core_discovery":"At leading order in the seesaw expansion, the Majoron's couplings to two gluons, two photons, a photon and a $Z$, two $Z$s, two $W$s, and to quark pairs of different generations all arise for the first time at two-loop order, and this paper gives them in closed form. The amplitudes are written as on-shell decay amplitudes $J \\to VV'$, they are ultraviolet finite by virtue of the neutrino-mixing identities of Eq. (5), and after tensor reduction they collapse to sums of elementary one-loop functions and rational terms. Two important structural results follow: $J \\to gg$ and $J \\to \\gamma\\gamma$ vanish like $m_J^2$ for small Majoron masses, matching onto derivative operators such as $(\\partial^2 J) F \\tilde F$ rather than $J F \\tilde F$, while $J \\to Z\\gamma$, $J \\to ZZ$, and $J \\to WW$ remain nonzero in the zero-momentum limit. The flavor-changing quark couplings are of minimal-flavor-violating form and are enhanced by $\\log(M_R/m_W)$, making $s \\to dJ$ the largest quark-level transition.","pith_inferences":["Implicit in the paper's hierarchy plots: because $g_{J\\gamma\\gamma}$ can be suppressed to near zero for particular hierarchies of the diagonal entries of $M_DM_D^\\dagger$, monochromatic gamma-ray line searches for Majoron dark matter may be far less sensitive than previously estimated, making the neutrino-line mode the more robust discovery channel.","A direct extension of the calculation would apply the same two-loop machinery to other spontaneously broken lepton numbers, such as familons or lepton-flavored majorons, where the flavor-changing quark couplings would inherit the same $\\log(M_R/m_W)$ structure.","The paper leaves the hadronic replacement as future work; because the free-quark amplitudes are the only input in the low-mass regime, a chiral or lattice calculation of the relevant matrix elements would settle the reliability of the low-mass exclusion regions.","The paper does not pursue it, but the explicit dependence of the flavor-changing couplings on $\\mathrm{tr}(M_D \\log(M_R/m_W) M_D^\\dagger)$ means that a single measured rare-decay rate, if combined with neutrino oscillation data, could begin to overconstrain the seesaw parameter space."],"forward_implications":["For a strictly massless Majoron, the new photon and gluon couplings vanish as $m_J^2$, and the paper concludes that the two-loop couplings are phenomenologically irrelevant there.","For nonzero Majoron masses, the flavor-changing quark couplings provide new rare-decay probes: $K\\to\\pi J$ and $B\\to K J$ can beat stellar-cooling limits when the $\\log(M_R/m_W)$ enhancement is active.","The full $J\\to\\gamma\\gamma$ coupling can be significantly smaller than the earlier Set II estimate when the diagonal entries of $M_DM_D^\\dagger$ are hierarchical, so dark-matter diphoton limits weaken and the neutrino-line channel becomes relatively more important.","In low-scale seesaw realizations, the special texture of $M_D$ makes all Majoron and non-Majoron coupling matrices share one flavor structure, so $\\tau\\to\\ell J$ and $\\mu\\to e J$ can be observable even where $\\tau\\to\\ell\\gamma$ and $\\mu\\to e\\gamma$ are not.","Because Majoron operators scale as $M_DM_D^\\dagger/f$ while ordinary seesaw operators scale as $M_D M_R^{-2}M_D^\\dagger$, the ratio $f/M_R$ controls whether Majoron couplings dominate or are negligible relative to standard seesaw observables."],"supporting_citations":[{"why":"Original definition of the singlet Majoron model and its tree-level neutrino couplings; this paper completes that program.","marker":"[3]"},{"why":"Computed the earlier Set II contribution to the diphoton coupling and established the Majoron dark matter phenomenology that this paper extends.","marker":"[10]"},{"why":"Gave the one-loop charged-lepton couplings and lepton-flavor-violating decay rates that feed the charged-lepton flavor-violating constraints used here.","marker":"[18]"},{"why":"Showed that Goldstone couplings to chiral gauge bosons are not fixed by anomalies, motivating the full perturbative two-loop calculation.","marker":"[40]"},{"why":"Supplied the expansion-by-regions method used for the double asymptotic expansion of the two-loop vertex integrals.","marker":"[46]"},{"why":"Supplied the tensor reduction of two-loop vacuum integrals used to reduce the amplitudes to one-loop functions and rational terms.","marker":"[48]"},{"why":"Provided the low-mass effective-field-theory form of the pseudoscalar-photon coupling used as a consistency check.","marker":"[52]"},{"why":"Gave the effective non-Majoron seesaw observables and limits used for comparison in the phenomenology section.","marker":"[17]"},{"why":"Supplied the current and future kaon and B-meson rare-decay constraints used in the limit plots.","marker":"[61]"}],"fun_headline_variants":["Two-loop Majoron couplings turn rare decays into seesaw probes","Flavor-changing decays expose two-loop Majoron couplings","Majoron rare decays at two loops probe seesaw scale","Two-loop Majoron: rare decays become seesaw probes"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is the free-quark description of the two-loop amplitudes: the paper's low-mass phenomenology, such as MeV-to-100 MeV $J\\to\\gamma\\gamma$ limits and $K\\to\\pi J$, assumes that replacing quarks by hadrons in the loops does not change the results by order-one factors, and that replacement is explicitly not done.","fun_headline_variants_meta":{"raw":{"variants":["Two-loop Majoron couplings turn rare decays into seesaw probes","Flavor-changing decays expose two-loop Majoron couplings","Majoron rare decays at two loops probe seesaw scale","Two-loop Majoron: rare decays become seesaw probes"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.001055,"raw_usage":{"total_tokens":4383,"prompt_tokens":857,"completion_tokens":3526,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":473,"completion_tokens_details":{"reasoning_tokens":3459}},"tokens_in":473,"tokens_out":3526,"duration_ms":26939,"temperature":1.0,"reasoning_tokens":3459,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T05:03:16.232042+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"An independent numerical evaluation of the two-loop amplitudes at a benchmark point (say $m_J=1\\,\\mathrm{GeV}$, $f=M_R=1\\,\\mathrm{TeV}$, diagonal $M_DM_D^\\dagger=(100\\,\\mathrm{GeV})^2$) would settle the formal calculation, since the paper's expressions are explicit enough to compare numerically. For the phenomenological layer, a hadronic computation of $J\\to\\gamma\\gamma$ between 1 and 100 MeV would decide whether the quoted exclusion curves survive.","supporting_citations":[{"cited_title":"Neutrino Mixing and the Axion-Gluon Vertex","cited_arxiv_id":"1203.3886","evidence_quote":"Showed that Goldstone couplings to chiral gauge bosons are not fixed by anomalies, motivating the full perturbative two-loop calculation."},{"cited_title":"Conformal Standard Model with an extended scalar sector","cited_arxiv_id":"1507.01755","evidence_quote":"Provided the low-mass effective-field-theory form of the pseudoscalar-photon coupling used as a consistency check."}],"review_version":1}