{"id":"4fb78d66-0414-40b4-81b2-45bfff319511","arxiv_id":"2606.24189","paper_version":1,"verdict":"UNVERDICTED","confidence":"LOW","novelty_score":4.0,"correctness_risk":"unknown","formal_verification":"none","parameter_count":0,"one_line_summary":"BS model calculations of Z boson radiative decays to charmonium yield higher branching fraction for J/ψ than prior theoretical results.","lead":"The paper calculates branching fractions for rare radiative decays of the Z boson into various S- and P-wave charmonium states using the Bethe-Salpeter model in instantaneous approximation. A smart generalist might read it to see how theoretical predictions for these decays compare to earlier work and what that implies for future collider searches.","discovery_kind":"new_application","skeptic_critique":{"model":"grok-4.3","headline":"Instantaneous BS approximation + LO diagrams may not reliably establish that prior J/ψ results were underestimated","rationale":"The reader's weakest_assumption directly identifies the load-bearing step. Because the full text was not examined here, no additional internal inconsistency or parameter-tuning artifact could be checked, so the existing flag remains the dominant concern and the provisional UNVERDICTED status is unchanged.","tokens_in":1602,"tokens_out":316,"duration_ms":17058,"concrete_test":"Recompute the J/ψ branching fraction after replacing the instantaneous kernel with a covariant ladder kernel (or adding the leading retardation correction) while keeping the same quark mass and coupling parameters; if the result moves by more than the difference quoted versus prior literature, the claim that earlier results were underestimated weakens.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim (higher J/ψ branching fraction than prior work, hence more promising) rests on decay amplitudes computed from BS wave functions in the instantaneous approximation using only leading-order Feynman diagrams. This approximation sets the relative energy component to zero and truncates the kernel, which can alter the momentum-space overlap integrals that enter the Z → charmonium + γ matrix element. Because the Z mass sets a hard scale and the photon carries away significant momentum, retardation and higher-order QCD/electroweak corrections omitted by construction could shift the rate by an amount comparable to the reported difference from earlier calculations. No sensitivity study or comparison to a covariant treatment is supplied to show the effect is under control.","agreement_with_reader":"agree"},"referee_report":{"model":"grok-4.3","summary":"The manuscript employs the Bethe-Salpeter model in the instantaneous approximation to compute branching fractions for the radiative decays Z → charmonium + γ. It treats S-wave states (ψ(nS), η_c(nS) for n=1,2,3) and P-wave states (χ_c0(mP), χ_c1(mP), h_c(mP) for m=1,2) via leading-order Feynman diagrams and concludes that the J/ψ branching fraction is larger than prior theoretical estimates, making it more promising for experiment.","tokens_in":1750,"tokens_out":495,"duration_ms":25302,"significance":"If the numerical results hold under the stated approximations, the work supplies updated predictions for rare Z decays that could guide searches at future e+e- colliders. The relativistic bound-state treatment via the BS equation is a standard tool in this area and the explicit listing of multiple channels provides a systematic survey.","major_comments":[{"comment":"Abstract and results section: the central claim that prior J/ψ results were underestimated rests on decay amplitudes obtained in the instantaneous BS approximation with only LO diagrams. Because the Z mass sets a hard scale and the photon carries O(M_Z) momentum, the neglect of retardation and higher-order corrections can shift the overlap integrals by an amount comparable to the reported difference; no sensitivity study or comparison to a covariant kernel is provided to demonstrate control.","section":"Abstract and results section"},{"comment":"Method of amplitude construction: the instantaneous approximation sets the relative energy component to zero by construction, which directly affects the momentum-space wave-function overlap that enters the Z → charmonium + γ matrix element. Without an explicit check (e.g., variation of the kernel or comparison to a non-instantaneous treatment) it is unclear whether the enhancement relative to earlier work is physical or an artifact of the truncation.","section":"Method of amplitude construction"}],"minor_comments":[{"comment":"Notation for the P-wave states should be clarified (e.g., explicit definition of the polarization tensors or radial wave functions used for χ_c1 and h_c).","section":"Notation"},{"comment":"A table comparing the new branching fractions directly with the numerical values from the cited prior works would make the “underestimated” statement quantitative rather than qualitative.","section":"Results"}],"recommendation":"major_revision","confidential_remarks":null},"author_rebuttal":{"model":"grok-4.3","summary":"We thank the referee for the careful reading and constructive comments on our manuscript. We address each major comment below, clarifying our use of the instantaneous Bethe-Salpeter framework while acknowledging its limitations.","responses":[{"response":"The instantaneous approximation is the standard reduction of the BS equation employed in the majority of charmonium decay calculations in the literature. Our enhancement for the J/ψ channel originates from the relativistic momentum-space wave functions solved with the chosen kernel, which differ systematically from the non-relativistic or quark-model inputs used in earlier estimates. We agree that a dedicated sensitivity study or covariant comparison would strengthen the robustness claim; we will therefore add an explicit paragraph in the revised manuscript discussing the expected size of retardation effects and the model dependence of the result.","revision_made":"partial","referee_comment":"[Abstract and results section] Abstract and results section: the central claim that prior J/ψ results were underestimated rests on decay amplitudes obtained in the instantaneous BS approximation with only LO diagrams. Because the Z mass sets a hard scale and the photon carries O(M_Z) momentum, the neglect of retardation and higher-order corrections can shift the overlap integrals by an amount comparable to the reported difference; no sensitivity study or comparison to a covariant kernel is provided to demonstrate control."},{"response":"Setting the relative energy to zero is the defining step of the instantaneous approximation and is applied uniformly to both the bound-state wave functions and the decay amplitude. This truncation is the same one used in our prior BS studies of charmonium transitions, where it reproduces known decay constants and widths to within typical hadronic uncertainties. The reported enhancement is therefore a direct consequence of the relativistic overlap integrals obtained inside this consistent framework rather than an uncontrolled artifact. A full non-instantaneous calculation lies outside the scope of the present work.","revision_made":"no","referee_comment":"[Method of amplitude construction] Method of amplitude construction: the instantaneous approximation sets the relative energy component to zero by construction, which directly affects the momentum-space wave-function overlap that enters the Z → charmonium + γ matrix element. Without an explicit check (e.g., variation of the kernel or comparison to a non-instantaneous treatment) it is unclear whether the enhancement relative to earlier work is physical or an artifact of the truncation."}],"tokens_in":1313,"tokens_out":530,"duration_ms":28588,"standing_objections":["A quantitative comparison to a covariant (non-instantaneous) Bethe-Salpeter kernel or an explicit retardation study would require a separate, substantially larger computational project."]},"desk_editor":{"model":"grok-4.3","letter":"The paper applies the instantaneous Bethe-Salpeter model to radiative Z decays into S- and P-wave charmonium states plus a photon, using leading-order diagrams to produce branching fraction predictions for several channels including psi(nS), eta_c(nS), chi_c0/1, and h_c.\n\nIt does the straightforward job of extending an established framework to these Z channels and reports specific numbers. The J/psi result is presented as larger than earlier estimates, which the authors say makes it more accessible to experiment.\n\nThe soft spot is the reliance on the instantaneous approximation together with only leading-order diagrams. Setting the relative energy component to zero and truncating the kernel can change the overlap integrals that control the decay rate, especially when the photon carries substantial momentum. Without a sensitivity study or comparison to a covariant treatment, it is hard to know whether the reported difference from prior work is robust or an artifact of the truncation. The abstract gives no indication that such checks were done.\n\nThis is a specialized phenomenology calculation. Readers who already use Bethe-Salpeter wave functions for charmonium might want the numbers for comparison, but the central claim about underestimated priors does not yet look solid enough to shift the literature.\n\nIt deserves a serious referee to examine the numerical implementation and the size of the omitted corrections.","headline":"BS model gives new branching fraction numbers for Z to charmonium gamma, but the claim that prior J/psi results were underestimated rests on untested approximations.","tokens_in":2217,"tokens_out":343,"would_cite":false,"duration_ms":8875,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"grok-4.3","headline":"Bethe-Salpeter calculations find a higher branching fraction for Z to J/ψ gamma than prior theory.","keywords":["Bethe-Salpeter equation","Z boson radiative decays","charmonium","branching fractions","J/psi","instantaneous approximation","S-wave states","P-wave states"],"falsifier":"A measured branching fraction for Z to J/ψ gamma that lies well below the Bethe-Salpeter prediction would show the calculation overestimates the rate.","tokens_in":2518,"feed_emoji":"⚛️","tokens_out":626,"duration_ms":36135,"temperature":0.7,"pith_summary":"The paper applies the Bethe-Salpeter model in its instantaneous approximation to compute branching fractions for radiative decays of the Z boson into S-wave charmonium states psi(nS) and eta_c(nS) as well as P-wave states chi_c(mP) and h_c(mP). The authors evaluate all amplitudes from leading-order Feynman diagrams and obtain numerical results for n and m up to 3. Their key result is that the J/ψ channel branching fraction exceeds earlier theoretical estimates, which would make observation more feasible at future experiments. A reader would care because these rare processes test how bound-state dynamics in heavy quarks couple to electroweak currents.","feed_headline":"Bethe-Salpeter model raises predicted Z to J/ψ gamma rate","feed_subtitle":"New calculations indicate earlier estimates were too low and improve prospects for detection.","key_machinery":"Instantaneous approximation to the Bethe-Salpeter equation, supplying the charmonium wave functions that enter the leading-order decay amplitudes.","core_discovery":"Within the instantaneous approximation of the Bethe-Salpeter equation the paper calculates the decay amplitudes for Z boson radiative transitions to charmonium using leading-order Feynman diagrams. For the J/ψ channel the resulting branching fraction is larger than previous theoretical results, indicating that those earlier estimates were likely too low.","pith_inferences":["The result implies that relativistic effects retained in the Bethe-Salpeter approach matter for electroweak decay rates of charmonium.","Similar calculations could be performed for radiative decays of the Higgs boson into the same states.","Tension with other quarkonium models would isolate the role of the instantaneous kernel choice."],"forward_implications":["Branching fractions are now available for all listed S- and P-wave channels as direct predictions.","The larger J/ψ rate raises the expected number of events at a future Z factory.","The same framework supplies consistent rates for the radial excitations n=2 and n=3.","Comparison with data on multiple channels can test the model uniformly across states."],"fun_headline_variants":["Bethe-Salpeter model boosts Z to J/ψ rate","BS calc finds higher branching for Z J/ψ gamma","Updated estimate increases Z radiative decay to J/ψ","Z to J/ψ gamma branching larger than prior theories","BS framework yields elevated Z J/ψ branching fraction"],"cache_read_input_tokens":2112,"weakest_assumption_plain":"The instantaneous approximation of the Bethe-Salpeter equation together with leading-order Feynman diagrams suffice to compute the decay amplitudes for all listed channels.","fun_headline_variants_meta":{"raw":{"variants":["Bethe-Salpeter model boosts Z to J/ψ rate","BS calc finds higher branching for Z J/ψ gamma","Updated estimate increases Z radiative decay to J/ψ","Z to J/ψ gamma branching larger than prior theories","BS framework yields elevated Z J/ψ branching fraction"]},"model":"grok-4.3","cost_usd":0.003919,"raw_usage":{"total_tokens":1956,"prompt_tokens":562,"num_sources_used":0,"completion_tokens":79,"cost_in_usd_ticks":39187000,"prompt_tokens_details":{"text_tokens":562,"audio_tokens":0,"image_tokens":0,"cached_tokens":256},"completion_tokens_details":{"audio_tokens":0,"reasoning_tokens":1315,"accepted_prediction_tokens":0,"rejected_prediction_tokens":0}},"tokens_in":562,"tokens_out":79,"duration_ms":7884,"temperature":1.0,"reasoning_tokens":1315,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-06-26T00:02:19.525435+00:00","model_set":{"reader":"grok-4.3"},"falsifier":"A measured branching fraction for Z to J/ψ gamma that lies well below the Bethe-Salpeter prediction would show the calculation overestimates the rate.","supporting_citations":[],"review_version":1}