2026-10-20 | | Total: 6
The finite density equation of state of QCD is a fundamental quantity for understanding the phase structure of the theory, and also for the realistic modeling of the QCD medium in simulations. While direct lattice simulations at finite density are extremely challenging, several expansion schemes have been devised to extrapolate into the QCD phase diagram. All of these schemes are based on susceptibilities of conserved charges, i.e. derivatives of the pressure with respect to the chemical potential(s). In this work we set out to determine which of these schemes possesses the best convergence properties. To this end, we employ extreme-statistics lattice simulations on a $16^3\times8$ lattice, which allow us to obtain susceptibilities up to order eight with remarkably small errors. We compare in particular two methods, based on contours of constant (normalized) baryon density and of constant entropy density. We find that the latter converges more quickly, effectively requiring susceptibilities only up to order four to cover the chemical potential range up to at least $μ_B \approx 500$ MeV. We also extend the comparison to the three dimensional space of chemical potentials, and find that the entropy-based method consistently shows superior convergence.
Unraveling the three-dimensional structure of hadrons is a central objective in modern nuclear and particle physics. Transverse Momentum Dependent parton distributions (TMDs) provide crucial insights into the interplay between partonic transverse motion and spin, and are indispensable for describing a wide range of high-energy processes at facilities such as JLab, COMPASS, the LHC, and the upcoming Electron-Ion Collider. In recent years, significant progress has been made in accessing TMDs from lattice QCD using novel approaches. Methods based on ratios of Mellin moments, large-momentum effective theory (LaMET), and Coulomb-gauge techniques have enabled first-principles investigations of the Collins-Soper kernel, intrinsic soft functions, TMD wave functions, and TMDPDFs. Recent lattice results for both unpolarized and spin-dependent distributions, including the Boer-Mulders function, demonstrate the feasibility of extracting transverse-momentum--dependent information directly from QCD. This paper reviews recent advances in lattice-QCD studies of TMDs, outlines remaining challenges such as continuum and infinite-momentum extrapolations, and discusses future directions toward precision determinations of hadron structure in three dimensions.
We extend the perturbative investigation of the renormalization and mixing patterns of the supercurrent operator in lattice ${\cal N}=1$ supersymmetric Yang-Mills (SYM) theory, using stout-smeared gauge links. We employ the SLiNC action for gluinos and a family of Symanzik improved gauge actions for gluons. Building upon previous studies employing clover-improved fermions and the Wilson gauge action, we compute the one-loop mixing matrix of the supercurrent operator among gauge-invariant and gauge-variant operators. Particular emphasis is placed on the structure of operator mixing in the presence of stout smearing parameters. Our results contribute toward refining fine-tuning strategies necessary for the restoration of supersymmetry in the continuum limit, and provide a stepping stone for future perturbative and nonperturbative investigations in supersymmetric theories on the lattice.
We present a quantum simulation of the SU(2) lattice gauge theory in 1+1 dimensions with one flavor of fermions, performed using the trapped-ion Quantinuum System Model H2. This work is based on the loop-string-hadron (LSH) formulation of the gauge theory, and is distinguished by using qubits to represent both bosonic-gauge and fermionic-matter degrees of freedom. Furthermore, we access the $j=1$ electric-flux sector beyond the $j_{\max}=1/2$ truncation commonly used in SU(2) quantum simulations. To that end, we digitize an electric-field degree of freedom using a multiqubit register for representing the bosonic quantum number. Our simulation therefore synthesizes several key ingredients relevant to quantum simulations in higher spatial dimensions, in particular explicit local gauge-field degrees of freedom represented here by a multiqubit bosonic register, while relying solely on circuit-construction optimizations that are expected to scale polynomially with lattice size. This work demonstrates that modern quantum hardware is increasingly able to accommodate the features of complicated but physically important lattice gauge theories while successfully reproducing the expected Trotterized dynamics.
The statistical long-range Ising (LRI) model in one dimension provides a family of some of the simplest non-exactly solvable critical points. These critical points are described by long-range conformal field theories (CFTs) that are well understood analytically with perturbative expansions and some exact results. Guided by this, we revisit the Monte Carlo (MC) approach to the corresponding lattice model which has many favorable features compared to other systems studied with MC methods. In particular, with very modest computational cost, we are able to estimate scaling dimensions of leading and subleading operators as well as the associated three-point functions, and even certain four-point functions. These results help highlight the fact that 1d LRI CFTs are a phenomenal playground to more generally test non-perturbative methods and ideas for strongly coupled CFTs.
We define string breaking by the midpoint electric field falling to half a flux unit. In the massive Schwinger model, an energy balance between intact string and screened configuration predicts this separation using only vacuum observables and meson masses as input. We observe numerically that this half flux point corresponds to equal weights of the two configurations thus explaining the observed peak of the entanglement entropy as maximal uncertainty of the binary flux. The same binary uncertainty explains the entanglement peak at the vacuum transition. Near breaking, screening organizes part of the spatial entanglement between the string interior and exterior into a nearly maximally entangled bit. Within the selected flux sectors, this approximate encoded Bell pair records whether the flux is screened and nearly separates from the remaining correlations. The vacuum transition exhibits the same structure. The string and the $θ$-vacuum reveal a matching structure of the capacity of entanglement minimum at their transition point, linking their entanglement spectra. The transition rearranges degrees of freedom coherently as revealed by the entanglement asymmetry and Robustness of Magic of disjoint subregions. Finally, we construct the $θ$-vacuum and string wavefunctions with high fidelity using a restricted Boltzmann machine ansatz.