Charmonium radiative transitions to dileptons from lattice QCD: The case of $h_c \to η_c \ell^+\ell^-$ and $χ_{c1} \to J/ψ\,\ell^+\ell^-$

Fuente: arXiv
Saved in:
Bibliographic Details
Main Authors: Bečirević, D., Di Palma, R., Frezzotti, R., Gagliardi, G., Lubicz, V., Sanfilippo, F., Tantalo, N.
Format: Preprint
Published: 2026
Subjects:
Online Access:
Tags: Add Tag
No Tags, Be the first to tag this record!
_version_ 1866913038501675008
author Bečirević, D.
Di Palma, R.
Frezzotti, R.
Gagliardi, G.
Lubicz, V.
Sanfilippo, F.
Tantalo, N.
author_facet Bečirević, D.
Di Palma, R.
Frezzotti, R.
Gagliardi, G.
Lubicz, V.
Sanfilippo, F.
Tantalo, N.
contents We present a lattice QCD study of dilepton production in charmonium transitions, specifically focusing on the $1^{+-} \to 0^{-+}$ and $1^{++} \to 1^{--}$ processes: $h_c \to η_c \ell^+ \ell^-$ and $χ_{c1} \to J/ψ\ell^+ \ell^-$, where $\ell = e, μ$. The relevant hadronic matrix elements are computed using gauge field configurations generated by the Extended Twisted Mass Collaboration with $N_f = 2+1+1$ dynamical Wilson--Clover twisted-mass fermions at four lattice spacings. Simulations are performed at physical dynamical $u$, $d$, $s$, and $c$ quark masses, except for the coarsest lattice, where the lightest sea quark mass corresponds to a slightly heavier pion mass. A controlled continuum extrapolation is carried out. In the continuum limit for the $h_c$ decays, we obtain $Γ(h_c \to η_c e^+ e^-) = 5.45(19)~\mathrm{keV}$, and $Γ(h_c \to η_c μ^+ μ^-) = 0.635(22)~\mathrm{keV}$. For the $χ_{c1}$ decays, we find: $Γ(χ_{c1} \to J/ψe^+ e^-)= 2.869(90)~\mathrm{keV}$, and $Γ(χ_{c1} \to J/ψμ^+ μ^-) = 0.1993(72)~\mathrm{keV}$. Our results for the $χ_{c1}$ decays show good compatibility with experimental data. However, our prediction for the $h_c \to η_c e^+ e^- $ decay rate is approximately $3σ$ larger than the BESIII result. We also present predictions for the differential decay widths as functions of the dilepton invariant mass, $q^2$, and for angular observables sensitive to longitudinal transition form factors, which are inaccessible in radiative decays with real photon emission. These results constitute the first fully dynamical lattice QCD predictions for dilepton decay rates in $h_c$ and $χ_{c1}$ charmonium transitions, including their differential distributions and angular observables. They provide benchmark predictions for future experimental studies.
format Preprint
id arxiv_https___arxiv_org_abs_2604_15256
institution arXiv
publishDate 2026
record_format arxiv
spellingShingle Charmonium radiative transitions to dileptons from lattice QCD: The case of $h_c \to η_c \ell^+\ell^-$ and $χ_{c1} \to J/ψ\,\ell^+\ell^-$
Bečirević, D.
Di Palma, R.
Frezzotti, R.
Gagliardi, G.
Lubicz, V.
Sanfilippo, F.
Tantalo, N.
High Energy Physics - Lattice
High Energy Physics - Phenomenology
We present a lattice QCD study of dilepton production in charmonium transitions, specifically focusing on the $1^{+-} \to 0^{-+}$ and $1^{++} \to 1^{--}$ processes: $h_c \to η_c \ell^+ \ell^-$ and $χ_{c1} \to J/ψ\ell^+ \ell^-$, where $\ell = e, μ$. The relevant hadronic matrix elements are computed using gauge field configurations generated by the Extended Twisted Mass Collaboration with $N_f = 2+1+1$ dynamical Wilson--Clover twisted-mass fermions at four lattice spacings. Simulations are performed at physical dynamical $u$, $d$, $s$, and $c$ quark masses, except for the coarsest lattice, where the lightest sea quark mass corresponds to a slightly heavier pion mass. A controlled continuum extrapolation is carried out. In the continuum limit for the $h_c$ decays, we obtain $Γ(h_c \to η_c e^+ e^-) = 5.45(19)~\mathrm{keV}$, and $Γ(h_c \to η_c μ^+ μ^-) = 0.635(22)~\mathrm{keV}$. For the $χ_{c1}$ decays, we find: $Γ(χ_{c1} \to J/ψe^+ e^-)= 2.869(90)~\mathrm{keV}$, and $Γ(χ_{c1} \to J/ψμ^+ μ^-) = 0.1993(72)~\mathrm{keV}$. Our results for the $χ_{c1}$ decays show good compatibility with experimental data. However, our prediction for the $h_c \to η_c e^+ e^- $ decay rate is approximately $3σ$ larger than the BESIII result. We also present predictions for the differential decay widths as functions of the dilepton invariant mass, $q^2$, and for angular observables sensitive to longitudinal transition form factors, which are inaccessible in radiative decays with real photon emission. These results constitute the first fully dynamical lattice QCD predictions for dilepton decay rates in $h_c$ and $χ_{c1}$ charmonium transitions, including their differential distributions and angular observables. They provide benchmark predictions for future experimental studies.
title Charmonium radiative transitions to dileptons from lattice QCD: The case of $h_c \to η_c \ell^+\ell^-$ and $χ_{c1} \to J/ψ\,\ell^+\ell^-$
topic High Energy Physics - Lattice
High Energy Physics - Phenomenology
url https://arxiv.org/abs/2604.15256