2026-11-10 | | Total: 16
We generalize the Akaike Information Criterion (AIC) to allow model and data subset selection when regularized fits are performed on correlated data. While the standard generalization to regularized fits is inherently dependent on the type of regularization applied to the data covariance matrix, our definition, the RIAIC, is, on average, regularization-independent. The regularization-independence is achieved by defining the new Information Criterion as an estimator of the model-dependent contribution to the correlated Kullback--Leibler divergence. The Bayesian extension is derived as well. We test the new criterion on lattice QCD data by performing regularized fits and show that it reproduces the selection properties of the standard correlated AIC, with fluctuations which are compatible with the expected statistical uncertainties. These results enable the application of subset and model selection techniques to cases where regularizations are applied to the covariance matrix of the fitted dataset.
We compute the pion quark momentum fraction in the large-$N_c$ limit using the Twisted Eguchi-Kawai (TEK) reduced model. Gauge configurations for different values of the lattice spacing are generated with the pure-gauge TEK action, while the pion correlators are computed for different quark masses with maximally twisted mass fermions. The bare matrix element corresponding to the quark momentum fraction is obtained using the Feynman-Hellmann method, and we use perturbative renormalization factors to connect to the $\overline{\mathrm{MS}}$ scheme at $2\,\mathrm{GeV}$. A global analysis at several lattice spacings, pion masses, and values of $N_c$ yields $\langle x\rangle_q^{\overline{\mathrm{MS}}}(2\,\mathrm{GeV}) = 0.267(105)$ at the physical pion mass in the large-$N_c$ limit. This value is significantly smaller than existing $N_c=3$ determinations from phenomenology and lattice QCD, indicating that, in the large-$N_c$ limit, the pion's light-cone momentum is carried predominantly by gluons.
We present a lattice quantum chromodynamics (QCD) calculation of the proton isovector ($u-d$) unpolarized parton distribution function (PDF) at physical sea and valence quark masses, using Coulomb-gauge-fixed quark bilinears within the large-momentum effective theory framework. The calculation is performed on a (2+1)-flavor ensemble with a lattice spacing of $a=0.076$ fm and proton momenta up to $2.3$ GeV. We combine kinematically enhanced nucleon interpolating operators, the Coulomb-gauge method, and high-statistics measurements to improve the signal quality at large momenta and quark separations. We use the static regularization-independent momentum subtraction and hybrid schemes to renormalize the bare matrix elements. The light-cone PDF in the $\overline{\rm MS}$ scheme is obtained from the quasi-PDF through next-to-leading-order (NLO) perturbative matching supplemented by next-to-leading-logarithmic renormalization-group resummation. With reduced uncertainties, the resulting PDF is consistent with our previous calculation on a similar lattice ensemble with $a=0.06$~fm and a valence pion mass of 300 MeV. Notably, despite the reduced uncertainties, the agreement with the phenomenological PDF persists over the longitudinal momentum fraction range $0.25\lesssim|x|\lesssim0.75$, where the lattice calculation is expected to be reliable. This agreement is further supported by a comparison in coordinate space, where our short-distance lattice matrix elements are consistent with those obtained from the phenomenological PDF through NLO matching. In particular, we observe close agreement in the imaginary part for the first time. Together, these results represent an encouraging step toward precision benchmarking of lattice QCD calculations of PDFs.
We propose a framework to construct the Symanzik effective Hamiltonian for lattice quantum field theories directly based on dimensional analysis and symmetries in the Hamiltonian formulation. This allows to classify all operators which can contribute at a given order in the spatial lattice spacing. This list can be reduced using the equations of motion and other identities. We exemplify the framework for $φ^4$ and pure gauge theory, obtaining results consistent with earlier publications for the case of classical improvement. It provides the means to develop a non-perturbative improvement programme in the Hamiltonian formulation of lattice quantum field theories in the future.
The calculation and analysis of nucleon correlation functions is essential for understanding the observed spectra and structure of these ubiquitous hadrons. Two challenges arise in their analysis: at early Euclidean times, they are dominated by excited state contributions and at later Euclidean times, the signal-to-noise ratio degrades exponentially. This leads to a common challenge in identifying a reasonable window on which to extract the ground state energy. By modulating the contribution of each Laplacian eigenmode to the correlator via distillation profiles, we build optimised operators which improve the overlap with the nucleon ground state. In these proceedings, we present our progress on constructing optimised distillation profiles for nucleons and their excited state suppression in two- and three-point correlation functions.
The degeneracy of hadron correlators in symmetry-related channels signals symmetry restoration in high-temperature QCD. RG-invariant probes of this degeneracy need care when the regularization breaks chiral symmetry, or when the symmetry is not exact in the massless action, as for the emergent chiral-spin symmetry $SU(2)_{CS}$. We renormalize the thermal and screening meson correlators of every flavor sector by the thermal correlator of the same channel in the $\bar ud$ sector at a common reference point. The symmetry ratio $κ_{AB}$ then depends on the renormalization only through the ratio $R_{AB}$. For exact symmetries, $R_{AB}=1$ and $κ_{AB}$ is RG invariant and scheme independent in mass-independent schemes. For emergent symmetries, $R_{AB}\neq1$ carries the relative logarithmic cutoff dependence of the two channels, which cancels in $κ_{AB}$, and the residual dependence on the reference scale is physical. The framework applies to Wilson, twisted-mass and staggered fermions, where $R_{AB}$ absorbs the channel-dependent renormalization, and the tensor pair $(T_k,X_k)$ keeps $R_{T_kX_k}=1$ on isotropic Wilson-type lattices. Lattice artifacts are a small correction to the chiral-spin probes $κ_{CS}$ and $κ_{fade}$, but as large as the signal of the exact-symmetry probes. For free quarks, we derive $κ_{fade}=2/3$ for every mass and separation, and $κ_{CS}\le1/3$. We apply the renormalized probes to $N_f=2+1+1$ QCD with optimal domain-wall quarks at the physical point, with $a=0.064$ fm and $T=192$-$1539$ MeV. The hierarchical restoration of the chiral symmetries from the light to the charm sectors is confirmed, and is invariant under renormalization. The $\bar uc$ and $\bar sc$ sectors have the widest chiral-spin windows, also at a tighter threshold. Thus spin-1 mesons with one charm quark are the most attractive channels for the chiral-spin regime.
Generative samplers for lattice field theory are commonly evaluated by their effective sample size (ESS). We show that a small global ESS need not imply inaccurate long-distance observables: errors spread over many short-distance modes can dominate this diagnostic. We introduce a scale-resolved diagnostic based on the overlap between coarse-grained model and target distributions. At the lattice cutoff, this overlap determines the asymptotic ESS fraction; lower cutoffs isolate long-wavelength errors. In the exactly solvable free theory, two samplers with comparable ESS yield a $32\%$ bias and no bias in the susceptibility, respectively. In two-dimensional $φ^4$ theory near criticality, trained diffusion samplers show the same separation of scales: their path-weight ESS fraction falls rapidly over the volumes studied, while low moments of the magnetization remain close to their target values. Varying the reference parameters substantially changes infrared observables, with little effect on the log-weight variance. Given the target magnetization distribution, a marginal Metropolis test maintains an acceptance rate of order one in our tests. Combined with exact local updates, it reduces the remaining mismatch, although residual bias remains at finite relaxation length. In the free field, uniformly bounded relative errors in the proposal mode variances give a volume-independent bound on the number of updates needed for a prescribed accuracy of the chain's stationary mode variances at fixed correlation length. The target marginal must be supplied or computed separately; estimating it from the samplers' own weights can inherit the poor global overlap.
We determine the nucleon axial and tensor charges using four ensembles of twisted-mass clover-improved fermions with approximately equal physical volumes of about 5.5~fm and four lattice spacings, $ a=0.05$~fm, 0.06~fm, 0.07~fm, and 0.08~fm. The mass of the degenerate up and down, strange and charm quarks are tuned to approximately their physical values. We compute both connected and disconnected contributions to determine the isovector and isoscalar charges and their flavour decomposition. For the nucleon isovector axial charge we find $g_A^{u-d}=1.263(34)(26)$, in agreement with the experimental value. For the isoscalar axial charge we find $g_A^{u+d}=0.458(16)(1)$, while for the corresponding tensor charges we find $g_T^{u-d}=0.951(21)(5)$ and $g_T^{u+d}=0.556(16)(5)$.
Background-field lattice calculations extract weak-field response coefficients from correlators computed at several field strengths. We treat this extraction as linear estimation in field space and compare fixed cancellation stencils with covariance-optimized, generalized-least-squares estimators of the same coefficient, using archival nucleon double ratios on about 1000 quenched configurations at six quark masses and decuplet correlators for which exact answers follow from the lattice contractions. The fixed $O(B^4)$-$O(B^8)$ stencils are nearly coincident and cost almost nothing in precision. Weight-free null-space tests reject the two-term response model in most nucleon ensembles; whether they reject the three-term model depends on the temporal fits and on the noise of the estimated temporal covariance, which raises the null-contrast error by a factor of 1.7-2.4. With that noise included, covariance optimization gains 1-2% in precision at the five lighter masses (3-30% with two-state time fits, which also lower moments by 3-10%) and nothing on independent data halves, while its weights carry three to four orders of magnitude more leverage on omitted response modes than the fixed stencils and move the moments by up to 1.8 marginal standard errors. The decuplet controls show what such leverage does under misspecification: the first-order responses obey $μ_{Δ^+}/μ_{Δ^{++}}=1/2$ and $μ_{Δ^0}=0$ exactly. The lowest-order estimate misses the first by 10-25 standard errors. The fixed $O(B^4)$-$O(B^8)$ stencils satisfy it within errors in the production windows, whereas the correlated linear fit in $B^2$, with the same constraints as the $O(B^4)$ stencil, misses it by 3-17%; the null-space test rejects that fit's model. The results distinguish statistical precision from robustness to omitted response modes; they are not physical-point magnetic moments.
Assessing the reliability of spectral reconstruction from finite, noisy lattice correlators requires understanding the constraints supplied by the data relative to prior assumptions. We develop an information-theoretic framework using Kullback--Leibler (KL) divergences to quantify the information supplied by the data through changes from prior to posterior probability distributions. Our analysis accounts for both changes in the probability distributions of spectral strength at individual energies and constraints encoded in correlations across the spectrum. We smear these information measures over neighboring energies using Gaussian weights to obtain a smoothed assessment at a specified energy scale. We demonstrate this framework using a hierarchical Bayesian reconstruction method with posterior sampling, applied to mock correlator data generated from a three-peak spectrum and a near-physical lattice proton correlator. These calculations show that the information measures help explain the properties of reconstructed spectra. They also provide a semiquantitative criterion for assessing how strongly reconstructed spectral structures are constrained by the data and how reliable they are in light of prior effects.
We compute the magnetic moments of the $Σ_c$ and $Σ_b$ isotriplets with external-field QCD sum rules. The heavy-quark propagator is kept with its full mass dependence and with its couplings to the photon, to soft gluons and to both simultaneously, and the magnetic sum rule is normalized by the mass sum rule of the same interpolating field, which removes the residue and gives the baryon masses within a few percent. The quark-condensate magnetic susceptibility $χ$, which otherwise dominates the uncertainty, is fixed by the proton magnetic moment computed with the same interpolating field; the neutron moment then follows to about $5\%$. Uncertainties from the QCD inputs are evaluated with a Monte Carlo analysis, and those from the interpolating field, the Borel window, the scale of the heavy-quark mass, the calibration and the normalization are added. With $\overline{\rm MS}$ heavy-quark masses we obtain $μ_{Σ_c^{++}}=2.37\pm0.65$, $μ_{Σ_c^{+}}=0.50\pm0.15$, $μ_{Σ_c^{0}}=-1.39\pm0.37$, $μ_{Σ_b^{+}}=2.16\pm0.65$, $μ_{Σ_b^{0}}=0.55\pm0.16$ and $μ_{Σ_b^{-}}=-1.06\pm0.32$ in nuclear magnetons; at this order the results depend sizably on the heavy-quark mass scheme. The moments are carried by the light diquark; the heavy-quark contribution, whose separation depends on the interpolating field, is $-0.13\pm0.02\,μ_N$ for $Σ_c$ and $+0.013\pm0.003\,μ_N$ for $Σ_b$.
We study negative-parity hidden-flavor pentaquarks with three to five heavy quarks, $QQ\bar{Q}qq$, $QQQ\bar{Q}q$, and $QQQQ\bar{Q}$ with $Q=c,b$ and $q=u,d$, in QCD sum rules, using the complete set of S-wave baryon--meson currents with $J^P=1/2^-$, $3/2^-$, and $5/2^-$ and the operator product expansion up to dimension 10, 9, and 6. We focus on the splittings between states of the same baryon--meson configuration, which, unlike the absolute masses, are only weakly affected by the choice of the heavy-quark mass scheme and by a finite width. Configurations whose light degrees of freedom have a unique spin are nearly degenerate or show splittings that decrease from charm to bottom by a factor compatible with $\bar m_c/\bar m_b$, whereas configurations with two light quarks in different hadrons show splittings that are not reduced, as expected from heavy-quark spin symmetry. A decomposition of the correlators into factorizable two-hadron and connected quark-exchange contributions shows that the correlators are dominated by the two-hadron part and that this pattern holds separately for both parts: the quark-exchange contribution decreases from charm to bottom when a heavy quark is exchanged and is not reduced when a light quark is exchanged. In the fully heavy systems the $1/2^-$ and $3/2^-$ states are degenerate within about $20$~MeV, while the $5/2^-$ state lies $126$--$139$~MeV ($cccc\bar c$) and $183$--$228$~MeV ($bbbb\bar b$) above the $3/2^-$ state. This splitting is generated by the exchange terms between the identical heavy quarks and does not decrease from charm to bottom, which points to the Pauli principle rather than to a hyperfine effect. Currents whose baryon and meson share no quark flavor are exactly factorizable at this order and cannot distinguish bound states from two-hadron states.
The construction of three-dimensional conformal field theories (3D CFTs) using interacting fermions in the lowest Landau level on a sphere has recently proven remarkably efficient, while it differs from that of conventional field theories on the fuzzy sphere. Focusing on the 3D Ising CFT, we consider a class of models that are defined by deformed Hamiltonians specified by a function which controls how the interaction is damped at high angular momenta, and this class includes the original model. We show that each model in the class can be reformulated in a basis-independent manner as a matrix--vector quantum-mechanical model, and, upon integrating out the fermionic vectors, as a matrix model. We then focus on a Gaussian choice, for which the bosonic matrices acquire a heat-kernel propagator whose width is controlled by a parameter $ξ$. Through exact diagonalization with the number of fermions up to $16$, we find that the low-energy spectra of the deformed Hamiltonians at criticality agree with those in the 3D Ising CFT for a range of $ξ$. We thus obtain matrix-model realizations of the 3D Ising CFT that, in contrast to conventional fuzzy-sphere field theories, show no sign of UV/IR mixing in the low-lying spectrum.
Gauge fixing in Yang-Mills theories encounters a non-perturbative obstruction: distinct gauge-equivalent configurations can satisfy the same gauge condition. These are the so-called Gribov copies and their existence challenges the global validity of the Faddeev-Popov prescription. Addressing this ambiguity motivates a framework that connects the geometry of gauge-field space to infrared dynamics. This mini-review introduces the Gribov problem in Euclidean Yang-Mills theories and traces the construction of the Gribov-Zwanziger action and its refinement through the introduction of dimension-two condensates. We discuss the resulting infrared propagators, their comparison to lattice results, and the role of BRST symmetry, before surveying selected applications and open questions.
The question of pairing between quarks inside baryons is addressed by calculating the fragmentation functions that describe production of nucleons (protons and neutrons) in high-energy collisions and comparing attendant predictions with phenomenological global fits to available data and directly with data on hadron multiplicities. The fragmentation functions are obtained from nucleon parton distribution functions using the Drell-Levy-Yan relation and subsequently solving a coupled set of pion, kaon, and nucleon jet cascade equations. Strong empirical evidence is found for the presence of diquark correlations. In comparison with global fits, $d$ (down) quark to proton fragmentation ($d\to p$) exhibits the clearest signals. The strongest evidence is found in comparisons with the momentum-dependent proton multiplicity distribution: FFs built from a proton wave function that contains diquark correlations deliver a $z$-profile that is fully consistent with available data, whereas that obtained from a no-pairing wave function is disfavoured by data at the $9σ$ level.
We compute the electromagnetic form factors of the pion and the kaon in a symmetry-preserving contact interaction whose effective coupling depends on the number of active light flavors Nf , so that dynamical chiral symmetry breaking is screened as Nf grows toward the chiral transition. The calculation uses the rainbow-ladder truncation of the Schwinger-Dyson and Bethe-Salpeter equations with proper-time regularization, and it includes the dressed quark-photon vertex obtained from its own inhomogeneous Bethe-Salpeter equation in the same truncation. At Nf = 2 the model reproduces the established contact interaction values of the meson masses, decay constants and condensates, and the charge radii rpi = 0.447 fm and rK = 0.417 fm, which undershoot experiment by about one third, as is known for this class of interactions. We therefore restrict every comparison to Q2 less then or equal to 0.6 GeV2, confront it with the NA7 data and with one-loop chiral perturbation theory, and focus on the response to Nf . Both radii grow monotonically with Nf , and the growth survives every variation of the screening function and of the infrared confinement scale that we test. The growth is carried by the quark core rather than by the vector-meson dressing: between Nf = 2 and 6 the bare-vertex contribution to <r2pi> rises from 0.097 to 0.226 fm2, while the dressing contribution moves from 0.102 to 0.127 fm2 even though the rho pole drops from 0.93 to 0.70 GeV. The finite strange mass slows the kaon response only moderately.