2026-04-17 | | Total: 3
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.
We present a minimal implementation of SU($N$) pure Yang-Mills theory in $3+1$ dimensions for digital quantum simulation, designed to enable quantum advantage. Building on the orbifold lattice simulation protocol with logarithmic scaling in the local Hilbert-space truncation, we introduce further simplified Hamiltonians. Furthermore, we test simple methods that improve the convergence to the infinite mass limit, thereby removing the requirement of a large scalar mass to obtain the Kogut-Susskind Hamiltonian. For the SU(2) theory, we can cut the resource requirement further by utilizing the embedding of $\mathrm{SU}(2)\cong\mathrm{S}^3$ into $\mathbb{R}^4$. Monte Carlo simulations of the Euclidean path integral were used to benchmark the accuracy of these new analytical improvements to the theory. These results provide further support for the noncompact-variable-based approach as a practical framework for quantum simulation of non-Abelian gauge theories.
We propose an effective Hamiltonian formulation of quantum field theories using a Daubechies wavelet basis in position space. Combined with flow-equation methods of the similarity renormalization group (SRG), this approach provides an efficient framework for analyzing quantum field theories by reducing the dimensionality of the Hamiltonian and systematically decoupling degrees of freedom across scales. As an application, the free scalar field theory has been reformulated within this framework to calculate the low-lying energy spectrum of the theory. These basis elements are known to transform the free scalar field theory into a theory of coupled localized oscillators, each of which is labeled by a location and a resolution index. In this representation, the Hamiltonian is naturally organized into fixed-resolution blocks, alongside blocks associated with the interactions between different resolutions. To decouple the different resolution modes and obtain a block diagonalized Hamiltonian with each block associated with a fixed resolution, the flow equation approach of SRG is applied. Finally, we demonstrate that with increasing resolution, the low-energy spectrum can be extracted from the effective lowest-resolution block of the Hamiltonian, leading to a significant reduction in computational cost.