2026-05-29 | | Total: 8
We present a lattice QCD determination of the Mellin moments of the unpolarized gluon parton distribution function in the proton. The analysis is based on matrix elements of nonlocal gluon operators coupled to momentum-boosted proton states. The calculation relies on an $N_f=2+1+1$ ensemble of maximally twisted mass fermions with clover improvement and the Iwasaki-improved gauge action, at a pion mass of approximately 260 MeV. Working within the short-distance operator product expansion (OPE) of the reduced gluon Ioffe-time distribution, we extract ratios of higher-order gluon moments, $\langle x^n\rangle$ with $n>1$, to the gluon momentum fraction, $\langle x\rangle$. We investigate systematic effects associated with the truncation of the order of moment in the OPE, the choice of minimum and maximum Wilson-line separations entering the analysis, and the treatment of mixing with the quark-singlet under perturbative matching. The stability of the extracted moments is further studied under scale evolution using DGLAP equations, allowing us to assess uncertainties related to perturbative truncation by varying the scale. Our work provides a determination of the ratio $\langle x^3\rangle_g/\langle x\rangle_g$ at a scale of 2 GeV, with uncertainties that account for both statistical and the dominant theoretical systematic uncertainties.
We present results on the Mellin moments of the unpolarized parton distribution function (PDF) of the pion and kaon up to the fourth order. The computation is done using one $N_f=2+1+1$ gauge ensemble of twisted mass fermions with quark masses tuned to approximately their physical values. We reconstruct the valence pion and kaon PDFs using the connected contributions to the three Mellin moments. We compare our results on the Mellin moments and the reconstructed PDFs with other lattice QCD and phenomenological determinations.
We present a lattice QCD calculation of the charged-neutral pion mass splitting $M_{π^+} - M_{π^0}$ at $\mathcal{O}(α_\mathrm{em})$ using a recently proposed framework based on a Pauli-Villars (PV) regulated photon propagator defined in the continuum and infinite-volume limit, with $Λ$ acting as an additional UV cutoff scale. The use of this propagator avoids power-law finite-volume effects, allowing for a straightforward treatment of the infinite-volume limit. We perform the calculation using CLS ensembles, studying finite-volume effects, the continuum limit and the extrapolation to the physical point for several values of the scale $Λ$. By means of the Cottingham formula, we further decompose the result into elastic and inelastic contributions at fixed $Λ$. Our final result, after removing the cutoff scale $Λ$, is $M_{π^+} - M_{π^0} = 4.56(22)$ MeV, in good agreement with the experimental measurement. This calculation serves as a validation of the formalism in a well-controlled setting and offers useful insights into the application of electromagnetic corrections to other observables.
This is a brief note on 'A better conditioned Domain Wall Operator', which provides a more detailed explanation of the domain wall-to-overlap transformation with the inclusion of the alpha parameter.
As a first step towards machine identification of confining objects in thermalized lattice gauge configurations, we present our 2dVoId model for center vortex identification on pure SU(2) lattices in $D = 2$ dimensions. We create a training set by inserting thin Z2 vortices at various locations on a zero action lattice, and then distort those configurations by applying random SU(2) gauge transformations, noise, and by thickening the vortices via cooling. For moderate vortex visibility, our model is able to reliably identify the location of center vortices. We additionally demonstrate scalability through tiling strategies, which will enable generalization to higher dimensions while reducing training costs.
Various physical models can be expressed in terms of matrices. A valuable tool for analysing matrix models is numerical simulations, often the Metropolis algorithm with various improvements. The downside of this approach is that the simulation may become stuck in a vacuous state, and probing the relevant parts of the configuration space might be difficult. Here, we propose an algorithm that moves around a cluster of eigenvalues and show that it converges to the true vacuum state.
We consider three-dimensional (3D) lattice Abelian Higgs models, with compact U(1) gauge variables coupled to a doubly-charged $N$-component complex scalar field (CLAH). We focus on their phase transitions between the disordered-confined (DC) and ordered-deconfined (OD) phases. When they are continuous they belong to the 3D Abelian Higgs (AH) universality class associated with the stable charged fixed point (CFP) of the renormalization-group flow of the 3D AH field theory, or scalar electrodynamics, describing $N$-component complex scalar fields minimally coupled to a U(1) gauge field. This CFP exists only for a sufficiently large number of components, i.e., $N \ge N_d^*$, where the integer $N_d^*$ depends on the spatial dimension $d$ (for example $N_4^*=183$). To estimate $N_3^*$, we look for the minimum number $N_{\rm cL}$ of scalar components of 3D doubly-charged CLAH models developing continuous transitions along their DC-OD transition line. For this purpose, we present finite-size scaling analyses of Monte Carlo simulations for $N\in[4,10]$, up to lattice sizes $L\approx 100$. The results provide evidence of continuous DC-OD transitions for $N=10$, and weak first-order transitions for $N\le 7$. They are not conclusive for $N=8,\,9$. Therefore, we estimate $N_{\rm cL}=9(1)$.
Finite-size effects limit the accuracy with which conformal data can be extracted from lattice simulations of critical systems. While action improvement suppresses some corrections to scaling, it does not address operator-dependent effects arising from imperfect lattice representations of continuum conformal fields. In this work, we propose a data-driven method for improving lattice operators themselves, constructing estimators with enhanced overlap with the corresponding primary operators of the continuum conformal field theory. We identify improved lattice representations of leading spin and energy operators in three two-dimensional critical systems: the Ising model, the q = 3 Potts model, and the dilute q = 3 Potts model. In all cases, the resulting operators exhibit reduced corrections to scaling and yield more accurate estimates of scaling dimensions compared to conventional lattice choices. The code and analysis workflows used to produce these results are made available in an accompanying GitHub repository.