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.

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.

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.

We present a lattice QCD calculation of the charged-neutral pion mass splitting $M_{\pi^+} - M_{\pi^0}$ at $\mathcal{O}(\alpha_\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 $\Lambda$ 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 $\Lambda$. By means of the Cottingham formula, we further decompose the result into elastic and inelastic contributions at fixed $\Lambda$. Our final result, after removing the cutoff scale $\Lambda$, is $M_{\pi^+} - M_{\pi^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.

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 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.

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.

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)$.

Parton distributions encode the momentum-space structure and, in their generalizations, the spatial tomography of quarks and gluons inside hadrons, the building blocks of visible matter. We present a unified neural-network approach that learns these distributions directly from matrix elements calculated via numerical simulations of quantum chromodynamics (QCD) on the lattice by fitting two complementary inputs simultaneously: data matched to physical quantities via known momentum-space and coordinate-space formalisms. Utilizing data from both methods stabilizes the extraction and mitigates biases that can arise when either is used alone. We validate the method on controlled mock data and apply it to lattice-QCD matrix elements to extract parton distribution functions (PDFs). We show benefits of such an approach for determining the physical quantities. We further extend the framework to zero-skewness generalized parton distributions and demonstrate nucleon tomography within the same neural-network parameterization. Our results provide an adaptable and systematically improvable approach for extracting partonic distributions from Euclidean correlators. It can incorporate polarization, additional channels, and future experimental constraints from current and future facilities, such as the Electron-Ion Collider.

We present the first lattice QCD determination of the $\Xi_b \to \Xi$ vector, axial-vector, and tensor form factors, which are relevant for the theory of rare decays including $\Xi_b \to \Xi \ell^+\ell^-$ and $\Xi_b \to \Xi \gamma$. The calculation is performed with 2+1 flavors of domain-wall fermions at three different lattice spacings and pion masses in the range from approximately 430 to 230 MeV. The bottom quark is implemented using an anisotropic clover action. Three-point functions with a wide range of source-sink separations and model averaging are used to extract the ground-state contributions. We fit the dependence of the form factors on the momentum transfer, the pion mass, and the lattice spacing using modified $z$ expansions that account for subthreshold branch cuts, and apply dispersive bounds and asymptotic-behavior constraints to achieve controlled uncertainties in the full semileptonic kinematic region. Using our form factor results, we present Standard-Model predictions for the $\Xi_b^- \to \Xi^- \gamma$ and $\Xi_b^- \to \Xi^- \mu^+\mu^-$ branching fractions and two angular observables.

We develop Microcanonical Hamiltonian Monte Carlo (MCHMC), a class of models which follow a fixed energy Hamiltonian dynamics, in contrast to Hamiltonian Monte Carlo (HMC), which follows canonical distribution with different energy levels. MCHMC tunes the Hamiltonian function such that the marginal of the uniform distribution on the constant-energy-surface over the momentum variables gives the desired target distribution. We show that MCHMC requires occasional energy conserving billiard-like momentum bounces for ergodicity, analogous to momentum resampling in HMC. We generalize the concept of bounces to a continuous version with partial direction preserving bounces at every step, which gives an energy conserving underdamped Langevin-like dynamics with non-Gaussian noise (MCLMC). MCHMC and MCLMC exhibit favorable scalings with condition number and dimensionality. We develop an efficient hyperparameter tuning scheme that achieves high performance and consistently outperforms NUTS HMC on several standard benchmark problems, in some cases by more than an order of magnitude.