2026-04-15 | | Total: 10
We propose a variant of the model-independent method for determining meson charge radii from spatial moments of correlation functions on the lattice. Traditional determinations based on fits to the momentum transfer squared dependence of form factors are subject to systematic uncertainties arising from the choice of fit ansatz. By contrast, model-independent methods based on spatial moments provide a useful framework for determining the slope of the form factor without assuming its functional form. Recently, Feng et al. proposed a model-independent method, which drastically suppresses the finite-volume effect in the charge radius coming from higher-order contributions of the expansion of the form factor with respect to the momentum transfer squared. In this work, we introduce an auxiliary function of the momentum transfer squared and reformulate the method in terms of its product with the form factor, rather than the form factor itself, thereby further suppressing higher-order contributions, notably in cases of small volume and large radius. In particular, we investigate quadratic and logarithmic forms as practical choices for this auxiliary function. Applying this method to mock data based on a monopole form factor, as well as to actual lattice QCD data using $N_f=2+1$ gauge ensembles at $m_π\simeq 0.5$ and $0.3$ GeV, we find that it reduces residual finite-volume effects and provides an effective framework for meson charge radius determinations.
Hilbert space fragmentation refers to exponential growth in the number of dynamically disconnected Krylov sectors with system size. It is taken as evidence of ergodicity breaking, since conventional symmetries generate at most a polynomial number of sectors. However, we demonstrate that generalized symmetries can fragment the Hilbert space. Models with higher-form, subsystem, and gauge symmetries can have exponentially many symmetry sectors. We further prove that non-invertible symmetries can induce additional fragmentation within individual symmetry sectors. Fragmentation in several known models arises from generalized symmetries, and the presence of exponentially many Krylov sectors therefore does not by itself imply ergodicity breaking. Finally, we show that disorder free localization arises naturally from Krylov-restricted thermalization when sectors lack translation invariance, requiring neither ergodicity breaking nor gauge symmetry.
In our recent work [1] on lattice QCD calculation of the baryon leading-twist LCDAs within the framework of LaMET, a novel hybrid renormalization scheme is implemented for octet baryon quasi-DAs, yielding reliable results across both perturbative and non-perturbative regions. The numerical simulations are performed using CLQCD ensembles with $N_f = 2+1$ stout-smeared clover fermions and a Symanzik-improved gauge action. Calculations are carried out at three lattice spacings, $a = {0.052, 0.077, 0.105}$ fm. After renormalization, the linear divergences inherent in quasi-DAs are effectively removed, leading to smooth and well-behaved continuum coordinate-space distributions. These results demonstrate the viability of hybrid renormalization frameworks for light-baryon quasi-DAs and provide a robust foundation for future LaMET-based determinations of baryon LCDAs.
In this review I summarize how machine learning can be used in lattice gauge theory simulations and what ap\-proaches are currently available to improve the sampling of gauge field configurations, with a focus on applications in four-dimensional SU(3) gauge theories. These include approaches based on generative machine-learning models such as (stochastic) normalizing flows and diffusion processes, and an approach based on renormalization group (RG) transformations, more specifically the machine learning of RG-improved gauge actions using gauge-equivariant convolutional neural networks. In particular, I present scaling results for a machine-learned fixed-point action in four-dimensional SU(3) gauge theory towards the continuum limit. The results include observables based on the classically perfect gradient-flow scales, which are free of tree-level lattice artefacts to all orders, and quantities related to the static potential and the deconfinement transition.
The complex Langevin approach is a promising method for the numerical treatment of systems with a sign problem, for which conventional lattice field theory techniques based on importance sampling cannot be applied. However, complex Langevin dynamics may fail to converge in some cases and converge to a wrong limit in others, motivating the development of various diagnostic tools over the years to assess the correctness of given simulation results. This work aims at providing a systematic comparison between the most prominent such correctness criteria. In particular, the main goal is to contrast their applicability, ease of use, and - most importantly - their predictive power. To this end, four simple but nontrivial models are considered and the criteria applied to each of them. The obtained conclusions are expected to carry over to more realistic theories as well.
We determine the interaction potential of a solitonic dipole in the singlet state, modeled as an SU(2) field, using improved lattice simulations of two stationary solitons at varying separations. The potential is extracted from the energy of two-soliton configurations as a function of distance. At large separations, the interaction reproduces the classical Coulomb potential quantitatively up to an energy shift $δE_\infty\approx 9\;\text{keV}$ of the fitted asymptotic constant relative to $2m_ec_0^2$, assumed to be related to limited numerical precision on the lattice. At shorter distances, deviations from the Coloumb potential of point-like charges appear, that are in qualitative agreement with the asymptotic formula of perturbative Quantum Electrodynamics, reflecting the running of the fine-structure constant, with the inverse fine-structure constant ($α^{-1} \approx 137$) reproduced.
The $U(1)_A$ symmetry of the massless QCD Lagrangian is explicitly broken by the axial anomaly, but it may be effectively restored at finite temperature. Determining the temperature at which this occurs is important for understanding the chiral transition and the structure of the QCD phase diagram. A commonly used probe of effective $U(1)_A$ restoration is the degeneracy of flavour non-singlet pseudoscalar and scalar susceptibilities. Using anisotropic lattice QCD ensembles with Wilson-clover fermions generated by the \textsc{Fastsum} collaboration, we study this degeneracy through hadronic correlation functions over a wide range of temperatures. The fine temporal resolution of our Generation 3 ensembles allows us to determine the temperature at which the pseudoscalar and scalar channels become degenerate. We find evidence for the effective restoration of $U(1)_A$ symmetry at $T_{U(1)_A}=319(22)$ MeV, well above the chiral crossover temperature.
The impact of open-flavor thresholds on the quarkonium spectrum has been a subject of study since the introduction of the Cornell potential and has been quantified through various phenomenological approaches, most notably the $^3P_0$ model. We revisit this problem using the Born--Oppenheimer effective field theory (BOEFT), an effective field theory systematically derived from QCD by exploiting hierarchies of energy scales and symmetries. Within the BOEFT, open-flavor threshold effects emerge from the mixing between quarkonium and tetraquark static potentials sharing the same Born--Oppenheimer quantum numbers. The shapes of the static potentials are constrained by lattice QCD calculations. Furthermore, we account for the distinctive behavior of the BOEFT tetraquark static potentials at short and large distances: at short distances they are repulsive, reflecting the color-octet configuration of the heavy quark-antiquark pair, while at large distances they asymptotically approach heavy-light meson-antimeson thresholds. To quantify threshold effects on the quarkonium spectrum below threshold, we solve a set of coupled Schrödinger equations dictated by the BOEFT, whose only free parameter, the adjoint meson mass, is fixed to the mass of the $χ_{c1}(3872)$ state. These coupled equations are solved both in the spin-isospin averaged threshold limit and, for the first time, including the spin splittings of the physical thresholds. We validate our results by computing the same threshold effects as self-energy corrections to the quarkonium propagator. We compare our predictions with existing experimental data and previous literature. Finally, we provide a field-theoretical interpretation of the pair-creation constant $γ$ appearing in the $^3P_0$ model.
The subprocess $nn\to ppe^-e^-$ is a key ingredient in the interpretation of nuclear neutrinoless double-beta decay. Intermediate $Δ$ resonances may provide additional enhancements to this transition. We take a first step toward a $Δ$-full description of $nn\to ppe^-e^-$ by investigating the neutrinoless double-beta decay $Δ^- \to p e^-e^-$ in the framework of chiral effective field theory. We systematically derive the long-range contribution from light-Majorana-neutrino exchange through loop diagrams and incorporate the short-range part through counterterms required by renormalization. We predict the pion-mass dependence of the decay amplitude in the kinematic configuration with collinear electrons. Furthermore, to facilitate lattice-QCD matching, we calculate the decay amplitude in the degenerate $Δ$-nucleon mass limit and provide the corresponding long-range prediction.
We calculate electromagnetic multipole moments of $Σ$-type strange hidden-charm pentaquarks $P^Σ_{ψs}$ (isospin triplet $Σ^+,Σ^0,Σ^-$) using QCD light-cone sum rules, with six (spin-1/2) and seven (spin-3/2) interpolating currents built from diquark-diquark-antiquark operators. We compute magnetic dipole $μ$ for all channels and, for spin-3/2, electric quadrupole ${\cal Q}$ and magnetic octupole ${\cal O}$ moments (first computation), and give the first quark-flavor decomposition. Scalar diquark currents yield charm-dominated, flavor-insensitive moments ($μ\in[-1.92,-1.21]μ_N$ for spin-1/2, $|μ|\lesssim1.2μ_N$ for spin-3/2), consistent with heavy-quark spin symmetry. Axial-vector diquark currents produce larger, flavor-sensitive moments with sign reversals governed by $e_u/e_d=-2$. For ${\cal Q}$, scalar-diquark currents give oblate deformations ($Q_0\approx-2.0\times10^{-2}{\rm fm}^2$) dominated by charm, while two-axial-vector-diquark currents predict prolate values up to $Q_0=+8.0\times10^{-2}{\rm fm}^2$, with sign reversal for $[su][uc]\bar{c}$ in two currents. Currents with scalar antiquark coupling yield a topology-independent octupole ${\cal O}\approx-0.25\times10^{-3}{\rm fm}^3$, a lattice QCD benchmark. Comparison with constituent quark models identifies four discriminants: $|μ|\gtrsim3μ_N$ in spin-1/2; sign of $μ$ for $[su][uc]\bar{c}$ in spin-3/2; non-zero ${\cal Q}$ (vanishes in $S$-wave molecules); and the ${\cal Q}$-${\cal O}$ sign correlation, probing $1/m_q$ weighting.