2026-05-08 | | Total: 7
We present a new set of lattice QCD ensembles with four flavors of smeared Möbius Domain Wall Fermions with good chiral symmetry and small fifth-dimensional extent. A modest amount of computing resources was sufficient to generate 30 publicly-available ensembles spanning five lattice spacings and a broad range of pion masses down to physical. To scrutinize our action we determine $F_{K^{\pm}}/F_{π^{\pm}} = 1.1962(34)$, a key quantity for precision CKM unitarity tests, heralding a future of inexpensive high-precision calculations of hadronic observables with chiral fermions.
The singlet axial form factor of the nucleon provides essential input for a complete understanding of the nucleon axial structure. Together with the isovector and isoscalar octet channels, in the forward limit it forms the basis for a full flavor decomposition of the proton spin. In this work we present a lattice QCD determination of the singlet axial form factor $G^{u+d+s}_A(Q^2)$ and related strange contribution $G^{s}_A(Q^2)$ using a set of $N_f = 2 + 1$ CLS gauge ensembles with $O(a)$-improved Wilson fermions, with a full error budget for the extrapolation to the chiral, continuum and infinite-volume limits. Particular focus is placed on the treatment of the disconnected contributions, which constitute the crucial element for the extraction of the strange component. Together with determinations of the isovector and isoscalar octet axial form factors, this work provides a comprehensive lattice QCD determination of the nucleon axial structure across different flavor channels.
Implicit score matching provides a computationally efficient approach for training diffusion models and generating high-quality samples from complex distributions. In this work, we develop a score-matching framework for SU(N) lattice gauge theories, which can be extended to other Lie groups. We apply the method to SU(3) gauge configurations with the Wilson gauge action in two and four dimensions and assess the quality of the generated samples by comparison with Hybrid Monte Carlo (HMC) simulations. We show that the diffusion models can be successfully trained and applied for sampling the Wilson gauge action. For large values of inverse coupling, accurate reverse-time integration requires predictor-corrector schemes, for which we introduce a corrector based on Hamiltonian molecular dynamics. While the corrector significantly improves sampling quality, it also increases the computational cost. We outline several strategies for improving sampling efficiency.
We propose a novel, machine-learning-based framework for constructing lattice fermions using Physics-Informed Neural Networks (PINNs). Our approach treats the formulation of the Dirac operator as an optimization problem guided by physical requirements, such as symmetries, locality and doubler-decoupling conditions. We first demonstrate that, when trained to satisfy the Ginsparg-Wilson (GW) relation as a soft constraint, a neural network reproduces the overlap fermion operator to high numerical accuracy and learns an effective sign-function mapping without explicitly using a prescribed polynomial or rational approximation. Secondly, we extend the framework from operator construction to machine-assisted algebraic discovery. Within a generalized polynomial ansatz, the network autonomously drives higher-order terms to zero and recovers the standard Ginsparg-Wilson relation. Remarkably, by changing the initial search bias, the same framework also finds a distinct solution corresponding to a Fujikawa-type generalized GW relation.
In the present work we explore for the first time the general structure and properties of the nonperturbative quark-gluon vertex in the complex plane. Specifically, we focus on the transversely-projected quark-gluon vertex that emerges from a recently developed symmetry-preserving approach for the study of meson properties beyond the rainbow-ladder approximation. The analysis focuses on the so-called "soft-gluon" limit, which reduces the momentum-dependence of the corresponding vertex form factors to a single momentum variable. The complexification of this variable inside the defining integrals furnishes unambiguously all eight vertex form factors within a concrete domain of the complex variable, delimited by a characteristic parabola. The extent of this reliable domain is determined by the appearance of the first singularity in the integrands of the vertex integrals, where the standard Wick rotation must be duly supplemented by additional crucial contributions. This primary analytic region may be extended considerably by resorting to standard extrapolation methods, which remain valid up until the appearance of complex structures associated with the onset of physical processes. The generalization of the method to arbitrary gluon momenta, and its relevance for the determination of the quark propagator in the complex plane, are briefly discussed.
We analytically calculate the all-loop bare perturbative part of the four-quark Bethe-Salpeter kernel using modern scattering amplitude methods. We work to subleading order in the large number of quark flavors approximation of massless Quantum Chromodynamics, which simultaneously makes an all-loop calculation feasible, is systematically improvable, and preserves asymptotic freedom. It also allows for avoiding the ambiguity of choosing a truncation scheme in Dyson-Schwinger equations. We exploit state-of-the-art methods in Integration-By-Parts reduction of Lorentz scalar Feynman integrals into a minimal Master Integral basis, and direct integration into Generalized Polylogarithms. As a byproduct of our calculation, we also provide the result for the gluon and quark propagators. We discuss a path towards nonperturbative formulation and potential future phenomenological applications.
We study the crystalline phase of the $O(2N)$ Gross--Neveu model with a chemical potential for $a \leq N-2$ of the fermions. We analyze the problem in three independent ways: using perturbative QFT methods, a semiclassical large $N$ analysis, and integrability techniques (both at finite and large $N$). The resulting picture is consistent across all three approaches: at sufficiently large chemical potential $h$, an inhomogeneous phase emerges in which $a$-particle bound states condense and which, at large $N$, corresponds to a periodically oscillating chiral condensate. In this phase, the usual dynamically generated scale $Λ$ is replaced by two new dynamically generated scales $Λ_{\rm n}$ and $Λ_{\rm c}$. These two scales govern the multiple nonperturbative effects in the theory, corresponding in particular to the mass gaps of neutral and charged excitations on top of the inhomogeneous vacuum, respectively. They also control the nonperturbative corrections to observables such as the free energy and provide the parameters characterizing the oscillatory profile of the mean field at large $N$. In this paper, we provide the necessary details of each of the three methods, thereby complementing the results announced in a previous, shorter publication.