{"id":"81d272c0-5497-4ad7-8977-4d51ab757056","arxiv_id":"2411.09124","paper_version":2,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":4.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"The invisible branching fractions of J/psi, psi(2S), and Upsilon(nS) to neutrino pairs are predicted at the 10^-8 to 10^-5 level, offering a way to measure the weak mixing angle at quarkonia mass scales.","lead":"This paper calculates how often the heavy quark-antiquark particles J/psi, psi prime, and Upsilon(nS) decay invisibly into neutrino pairs, finding tiny Standard Model rates. It argues that measuring these rates with future collider data could pin down the weak mixing angle, a key electroweak parameter, at a previously unmeasured energy scale.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Use of alpha(0) instead of the running alpha(M^2) in Eq. (8) introduces a 7-11% systematic shift in the predicted invisible branching fractions, comparable to the assumed experimental precision.","rationale":"I verified that Table I follows from Eq. (8): for Upsilon(1S) with t3 = -1/2, q_b = -1/3, sin^2 = 0.233, g_V = -0.3447, the formula gives 1.01e-5, matching. So the arithmetic is not the issue. The reader's weakest_assumption, that the same f_V governs photon and Z channels, is less problematic than it appears: both amplitudes use the same local vector current, so QCD corrections to the annihilation factorize and cancel in the ratio at leading twist. The concrete, self-identified weak point is the replacement of the running coupling alpha(M^2) by alpha(0). This is not a minor convention: because Eq. (8) contains 1/alpha^2 and the photon virtuality is M^2, the choice changes the central predictions by 7-11%, the same order as the 10% measurement accuracy assumed for the sin^2(theta_W) projections. The extracted mixing angle is thereby biased by an amount comparable to the quoted uncertainties. This does not invalidate the method, since using the correct alpha restores the predictions, but it means the paper's precision claims are incomplete until the running-alpha systematic is included and folded into Table I. I therefore keep the CONDITIONAL verdict. The suggested check, recomputing with alpha(m_c) and alpha(m_b), would settle the size of the bias in one step.","tokens_in":8598,"tokens_out":15771,"duration_ms":203279,"concrete_test":"Recompute Table I and the Delta sin^2(theta_W) columns using alpha(m_c) ~ 1/132 and alpha(m_b) ~ 1/130 (or, better, the full time-like running alpha including hadronic vacuum polarization) instead of 1/137.036. If the resulting BR(V->nu nubar) values decrease by more than 5% and the inferred sin^2(theta_W) shifts by more than the quoted statistical projections, the paper must state the running-alpha choice as a systematic uncertainty and revise the precision claims.","verdict_should_be":"UNCHANGED","load_bearing_attack":"Eq. (8) is the central prediction, and the paper's extraction proposal relies on its normalization to the measured BR(V->e+e-). The derivation cancels f_V between Eqs. (2) and (7), but the alpha appearing in Eq. (8) is the coupling of the virtual photon in the normalization channel, whose virtuality is M^2, not zero. The paper states 'We have used alpha = 1/137.036 ... instead of the running alpha', but does not quantify the effect. With alpha(m_c) ~ 1/132 and alpha(m_b) ~ 1/130, [alpha(0)/alpha(M)]^2 ~ 0.93 for J/psi and ~0.90 for Upsilon(1S), so the tabulated branching fractions are high by ~7% and ~11%. This is the same size as the assumed 10% measurement uncertainty, and when propagated into the extraction it moves sin^2(theta_W) by roughly +0.005 (J/psi) and +0.03 (Upsilon(1S)), i.e. at or beyond the quoted +/-0.007 and +/-0.027. The correct time-like alpha(M^2), including hadronic vacuum polarization and final-state QED corrections, is needed before the projected uncertainties in Table I can be taken at face value.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript computes Standard Model branching fractions for invisible decays of vector quarkonia (J/psi, psi(2S), Upsilon(1S), Upsilon(4S)) into neutrino pairs. Using vector-meson dominance, the authors obtain the compact formula of Eq. (8), in which the quarkonium decay constant cancels against the measured e+e- branching ratio. The predicted values are BR(J/psi -> nu nubar) = 2.04e-8 and BR(Upsilon(1S) -> nu nubar) = 1.02e-5. The authors propose that a 5-10% measurement of these invisible branching fractions would determine sin^2 theta_W at the quarkonia mass scale with uncertainties of roughly 0.007 (charmonia) and 0.027 (bottomonia). They also discuss gamma-Z interference in the charged-lepton modes and argue that non-standard neutrino couplings could, in principle, distinguish Dirac and Majorana neutrinos.","tokens_in":8858,"tokens_out":16967,"duration_ms":176211,"significance":"If the central formula is correct, the paper offers a clean, hadron-uncertainty-free ratio method for a low-energy determination of the weak mixing angle, complementary to NuTeV and coherent neutrino scattering. The cancellation of the decay constant in Eq. (8) is elegant, the numerical predictions are falsifiable, and the comparison with the existing BaBar upper limits is informative. The significance is, however, conditional on correcting an unquantified running-alpha systematic, on including the electroweak rho corrections in the extraction, and on a more realistic assessment of experimental backgrounds.","major_comments":[{"comment":"The numerical predictions use alpha = 1/137.036 instead of the time-like running fine-structure constant at the quarkonium mass scale. Since the ratio in Eq. (8) is built from the theoretical leptonic width in Eq. (2) and then multiplied by the measured BR(V_Q -> e+e-), the alpha that appears is the coupling of the virtual photon at virtuality M^2, not at zero momentum. With alpha(m_c) ~ 1/132 and alpha(m_b) ~ 1/130, the tabulated branching fractions are systematically high by about 7% for charmonia and 11% for bottomonia. For the assumed 10% measurement uncertainty, this shift is comparable to the quoted precision; propagated into the extracted sin^2 theta_W, it moves the central value by roughly +0.005 (J/psi) and +0.03 (Upsilon(1S)), at or beyond the quoted +/-0.007 and +/-0.027. Please replace alpha(0) by the appropriate scale-dependent alpha(M) (including hadronic vacuum polarization and final-state QED corrections) and recompute Table I and the projected sensitivities.","section":"Section 2, Eq. (8) and Table I"},{"comment":"The numerical analysis in Table I uses the tree-level relation g_V^Q = t_3 - 2 q_Q sin^2 theta_hat, while footnote 1 correctly states that electroweak corrections modify this to sqrt(rho_f) (t_3 - 2 q_f sin^2 theta_hat). The rho_f factor shifts g_V^Q by about a percent, which translates into a shift in the extracted sin^2 theta_W of roughly 0.001 for J/psi and 0.003 for Upsilon(1S). These shifts are smaller than the 10% statistical projections but become a non-negligible fraction of the 5% scenario and should be either implemented in the formula or explicitly bounded before the precision statement in the conclusions can be taken at face value.","section":"Section 2, footnote 1 and Eq. (8)"},{"comment":"The paper states that the predicted branching fractions 'could be within the reach of current and future data sets at BES-III and Belle-II' and that with 10^12 J/psi and a 1-10% detection efficiency the weak mixing angle can be determined. No estimate of backgrounds is provided for the invisible decay signature, such as e+e- -> gamma gamma with both photons escaping detection, radiative Bhabha with undetected tracks, or other two-photon processes. Since the signal rates are at the 1e-8 to 1e-5 level and the extraction requires a few-percent measurement of a completely invisible final state, the feasibility claim is not yet supported. Please either add a quantitative background estimate or soften the reach claim accordingly.","section":"Introduction and Conclusions, experimental reach"}],"minor_comments":[{"comment":"There is a typo: 'botommonium' should be 'bottomonium'.","section":"Introduction"},{"comment":"The section title contains 'DECA Y' and should read 'DECAY'.","section":"Section on Majorana neutrinos"},{"comment":"The notation V_Q is used both for the quarkonium state and, in Eq. (10), for the factor multiplying the electromagnetic amplitude. This is confusing and should be clarified, for example by renaming the factor in Eq. (10) to R_V.","section":"Eqs. (9)-(11)"},{"comment":"The statement that the predicted Upsilon(1S) rate is 'only one order of magnitude below' the BaBar limit is imprecise: the limit 3.0e-4 is a factor of about 30 above the prediction 1.02e-5, i.e., about 1.5 orders of magnitude.","section":"Table I and text after it"},{"comment":"The sentence 'Precise measurements of decay fractions would indeed provide a competitive determination at those scales, but even a not-so-precise measurement could give us for the first time indications of the central value' is vague; it would be useful to state the minimum precision required for a meaningful determination.","section":"Conclusions"}],"recommendation":"major_revision","confidential_remarks":"The paper is a compact theory letter with a clean central formula, but the unquantified running-alpha systematic directly affects the headline precision numbers and the experimental reach discussion is optimistic. The self-citation to Ref. [24] is appropriate because the amplitude template originates there. The manuscript is better suited to a journal that accepts short phenomenological letters; after the alpha and rho corrections are implemented and the reach claim is tempered, it could be suitable for publication."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Worth knowing: the paper's real contribution is not the invisible-width formula—that goes back to Refs [5-11]—but the explicit use of quarkonium invisible decays as a weak-mixing-angle probe, with concrete precision projections. The derivation is clear, the f_V cancellation is neat, and I checked the Table I numbers: they follow from Eq. (8). The experimental context is sensible: the BaBar upper limit on Upsilon(1S)->invisible is only an order of magnitude above the SM prediction, and BES-III/Belle-II event counts make the measurement not absurd.\n\nCredit where due: this is a legitimate, honest SM calculation. It cites the earlier invisible-width estimates rather than hiding them, it keeps the Z propagator mass instead of dropping it, and the Majorana neutrino discussion is correct—the Dirac-Majorana difference is neutrino-mass-suppressed, and only non-standard couplings could evade the confusion theorem. The amplitude template borrowed from Kim's general Z->nu nu analysis is used in a genuinely different context, not fitted to the prediction, so I don't see a circularity problem.\n\nSoft spots, in proportion:\n\nThe alpha issue is real and the paper flags it without quantifying it. Eq. (8) uses alpha = 1/137.036 while the normalization channel has virtuality M^2. Using alpha(M^2) reduces the predicted branching fractions by roughly 7% for J/psi and 11% for Upsilon(1S). That is the same order as the assumed 10% measurement uncertainty, and it shifts the extracted sin^2(theta_W) by about +0.005 and +0.03—at or beyond the quoted +/-0.007 and +/-0.027. The Table I projections cannot be taken at face value until the running time-like QED coupling, including hadronic vacuum polarization, is used.\n\nSecond, the same f_V is assumed to govern the Z and photon matrix elements. That is the standard vector-meson-dominance assumption and likely fine at the few-percent level once f_V cancels in the ratio, but the QCD and electroweak corrections to the weak vertex are not estimated. This is a minor-to-moderate gap, not a fatal one.\n\nThird, the novelty claim in the conclusions, that invisible decays for this purpose have not been considered in any literature before, is too strong as written. The paper does cite Refs [5-11], and the specific sin^2(theta_W) extraction goal is genuinely not in those references, so this is a wording issue, not a substantive one.\n\nThe reader's conditional verdict is fair: the central idea holds up, but the precision projections need revision. This paper deserves a serious referee. Send it to peer review, and make sure the referee demands a corrected running-alpha treatment plus an estimate of the vertex corrections before the numbers in Table I are used.","headline":"A clean, checkable SM calculation with a modestly new angle-extraction proposal, but the alpha(0) choice biases the central numbers at the level of the claimed precision.","tokens_in":9381,"tokens_out":2920,"would_cite":true,"duration_ms":41386,"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":"Invisible decays of vector quarkonia could determine the weak mixing angle at quarkonium mass scales.","keywords":["invisible decays","vector quarkonia","weak mixing angle","neutrino pairs","J/psi","Upsilon(1S)","Dirac-Majorana confusion","Standard Model predictions"],"falsifier":"A high-statistics $e^+e^-$ sample that limits $\\mathrm{BR}(\\Upsilon(1S)\\to\\text{invisible})$ below about $1\\times10^{-5}$ at 90% C.L., with no new invisible particles, would falsify the prediction; equally, a measurement that breaks the predicted scaling across $J/\\psi$, $\\psi(2S)$, and $\\Upsilon(1S)$ would falsify the ratio formula.","tokens_in":8389,"feed_emoji":"⚛️","tokens_out":12481,"duration_ms":120518,"temperature":0.7,"pith_summary":"This paper sets out to show that the otherwise-unobserved decays of vector quarkonium states into neutrino pairs are clean Standard Model rates that change predictably with the weak vector coupling of the heavy quark, and that measuring them would give the first determination of the weak mixing angle at charmonium and bottomonium mass scales. The key formula, Eq. (8), relates each invisible branching fraction to the measured branching fraction into electron-positron pairs, with the quarkonium annihilation constant cancelling out, so the prediction is essentially free of strong-interaction uncertainty. The paper's numerical predictions are $\\mathrm{BR}(J/\\psi\\to\\nu\\bar\\nu)=2.04\\times10^{-8}$ and $\\mathrm{BR}(\\Upsilon(1S)\\to\\nu\\bar\\nu)=1.02\\times10^{-5}$, which would require a few-percent measurement to compete with existing low-energy determinations of $\\sin^2\\theta_W$. It also shows that within the Standard Model the Dirac and Majorana neutrino rates coincide in the massless limit, and that only new neutrino couplings could separate them, at a precision beyond current facilities.","feed_headline":"Rare invisible quarkonium decays could fix the weak mixing angle","feed_subtitle":"Predicted branching fractions of about 2×10⁻⁸ and 1×10⁻⁵ put the weak mixing angle within reach.","key_machinery":"The load-bearing object is the vector quarkonium annihilation constant $f_{V_Q}$, defined by $\\langle 0|\\bar Q\\gamma_\\mu Q|V_Q\\rangle = M^2 f_{V_Q}^{-1}\\epsilon_\\mu$. It is the same hadronic matrix element appearing in $V_Q\\to\\ell^+\\ell^-$ and $V_Q\\to\\nu\\bar\\nu$, and because Eq. (8) is the ratio of the two widths, $f_{V_Q}$ cancels exactly. What is left is the weak vector current, so only $g_V^Q$ matters and the axial coupling drops out. The Z propagator is kept with finite quarkonium mass $M$, producing the factor $m_Z^2/(m_Z^2-M^2)$ in the amplitude.","core_discovery":"The central claim is Eq. (8): the invisible branching fraction of a vector quarkonium state is fixed by its measured electromagnetic annihilation rate, $$\\mathrm{BR}(V_Q\\to\\nu\\bar\\nu)=N_\\nu\\left[\\frac{G_F $M^{2}$}{4\\pi\\$\\alpha$}\\left(\\frac{g_V^Q}{q_Q}\\right)\\left(\\frac{$m_Z^{2}$}{$m_Z^{2}$-$M^{2}$}\\right)\\right]^2 \\mathrm{BR}(V_Q\\to e^+e^-),$$ with $N_\\nu=3$ light neutrinos. Because the same decay constant $f_{V_Q}$ enters the photon- and Z-mediated amplitudes, taking the ratio cancels it; the remaining dependence is the heavy-quark weak vector charge $g_V^Q=t_3^Q-2q_Q\\sin^2\\theta_W$, so an invisible-rate measurement reads $\\sin^2\\theta_W$ at the quarkonium mass scale. The numerical table gives $2.04\\times10^{-8}$, $5.52\\times10^{-9}$, $1.02\\times10^{-5}$, and $1.05\\times10^{-8}$ for $J/\\psi$, $\\psi(2S)$, $\\Upsilon(1S)$, and $\\Upsilon(4S)$, respectively. The paper also derives the Majorana version and finds the Standard Model Dirac-Majorana difference is proportional to $m_\\nu^2$, vanishing for massless neutrinos; with nonstandard couplings the difference scales with $\\epsilon_V-\\epsilon_A$ and is too small for current experiments.","pith_inferences":["Beyond the paper, the same ratio formula should apply, with obvious changes of mass and charge, to other $1^{--}$ vector states such as $\\Upsilon(2S)$ and $\\Upsilon(3S)$ for which $e^+e^-$ branching fractions are measured, giving cross-checks of the extracted $\\sin^2\\theta_W$ at several nearby scales.","Beyond the paper, if QCD or electroweak corrections to the Z-mediated vertex differ from the photon-mediated one, Eq. (8) would carry a correction not visible in the ratio; comparing the extracted $\\sin^2\\theta_W$ across charmonium and bottomonium states could expose such an effect.","A natural testable extension is to push the search for $\\Upsilon(1S)\\to\\nu\\bar\\nu$ with tagged radiative or dipion transitions at high-luminosity $e^+e^-$ colliders, since an improved upper limit near $10^{-5}$ would either confirm the prediction or require new invisible decay products."],"forward_implications":["A 10% measurement of $\\mathrm{BR}(J/\\psi\\to\\nu\\bar\\nu)$ would determine $\\sin^2\\theta_W$ at the charmonium scale to about $\\pm0.007$; a 5% measurement would give $\\pm0.003$, comparable to existing low-energy determinations.","For $\\Upsilon(1S)$, the predicted $1.02\\times10^{-5}$ sits one order of magnitude below the current experimental upper limit, so a moderate increase in $e^+e^-$ statistics could either confirm the Standard Model prediction or constrain new invisible decay modes.","The cancellation of $f_{V_Q}$ removes the dominant hadronic uncertainty, making these decays a relatively clean electroweak precision observable at low energy scales, complementary to neutrino-scattering determinations.","Within the Standard Model, the invisible width of quarkonia cannot distinguish Dirac from Majorana neutrinos; only nonstandard neutrino couplings create a difference, and that difference is smaller than the expected experimental precision."],"supporting_citations":[{"why":"Supplies the measured $V_Q\\to e^+e^-$ branching fractions used as inputs in Eq. (8) and the reference running of $\\sin^2\\theta_W$.","marker":"[4]"},{"why":"Provides the 90% C.L. upper limits on $J/\\psi$ and $\\psi(2S)$ invisible branching fractions that the predictions are compared against.","marker":"[15]"},{"why":"Reports the current upper limit on $\\Upsilon(1S)\\to\\text{invisible}$, the benchmark closest to the predicted Standard Model rate.","marker":"[16]"},{"why":"Supplies the MS-scheme running values of $\\sin^2\\hat\\theta_W(\\mu)$ used at each quarkonium mass scale in the numerical table.","marker":"[20–22]"},{"why":"Gives the general Lorentz-, CP-, and CPT-invariant $Z\\to\\nu\\bar\\nu$ amplitude form that the quarkonium amplitude is mapped onto.","marker":"[24]"},{"why":"Provides the companion general parametrization of the $Z\\to\\nu\\bar\\nu$ amplitude with possibly modified couplings.","marker":"[25]"},{"why":"Establishes the Dirac-Majorana confusion theorem, which the paper invokes to state that massless-neutrino Standard Model rates do not depend on neutrino nature.","marker":"[26]"},{"why":"Supplies the electroweak precision data and the measured invisible width of the Z used to constrain nonstandard neutrino couplings in the Majorana analysis.","marker":"[12]"}],"fun_headline_variants":["Invisible quarkonium decays probe weak mixing angle","Quarkonia neutrinos: a new weak angle ruler","Rare quarkonium decays set sight on sin²θW","Vector quarkonia decays map weak mixing angle","Neutrino-pair decays of quarkonia weigh weak angle"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that a quarkonium state annihilates into a Z boson exactly as it annihilates into a photon, so that no QCD or electroweak correction shifts the weak process relative to the electromagnetic one.","fun_headline_variants_meta":{"raw":{"variants":["Invisible quarkonium decays probe weak mixing angle","Quarkonia neutrinos: a new weak angle ruler","Rare quarkonium decays set sight on sin²θW","Vector quarkonia decays map weak mixing angle","Neutrino-pair decays of quarkonia weigh weak angle"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000149,"raw_usage":{"total_tokens":1203,"prompt_tokens":965,"completion_tokens":238,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":581,"completion_tokens_details":{"reasoning_tokens":160}},"tokens_in":581,"tokens_out":238,"duration_ms":2962,"temperature":1.0,"reasoning_tokens":160,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-12T21:01:42.552727+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A high-statistics $e^+e^-$ sample that limits $\\mathrm{BR}(\\Upsilon(1S)\\to\\text{invisible})$ below about $1\\times10^{-5}$ at 90% C.L., with no new invisible particles, would falsify the prediction; equally, a measurement that breaks the predicted scaling across $J/\\psi$, $\\psi(2S)$, and $\\Upsilon(1S)$ would falsify the ratio formula.","supporting_citations":[{"cited_title":"A Search for Invisible Decays of the Upsilon(1S)","cited_arxiv_id":"0908.2840","evidence_quote":"Reports the current upper limit on $\\Upsilon(1S)\\to\\text{invisible}$, the benchmark closest to the predicted Standard Model rate."},{"cited_title":"Probing the non-standard neutrino interactions using quantum statistics","cited_arxiv_id":"2209.10110","evidence_quote":"Provides the companion general parametrization of the $Z\\to\\nu\\bar\\nu$ amplitude with possibly modified couplings."},{"cited_title":"Kayser, Phys","cited_arxiv_id":null,"evidence_quote":"Establishes the Dirac-Majorana confusion theorem, which the paper invokes to state that massless-neutrino Standard Model rates do not depend on neutrino nature."}],"review_version":1}