A major goal of the quantum simulation of high-energy physics (HEP) is to probe real-time nonperturbative far-from-equilibrium quantum processes underlying phenomena such as hadronization in quantum chromodynamics (QCD). The quantum simulation of the dynamics of confining strings and glueballs, both essential aspects of quark confinement, in a controllable first-principles way is an important step towards this goal. Here, we realize a $\mathbb{Z}_2$ lattice gauge theory in $2+1$D with a tunable plaquette term on a \texttt{Quantinuum System Model H2} trapped-ion quantum computer. We implement a shallow depth-6 Trotter circuit on a $6 \times 5$ matter-site square lattice utilizing all $56$ available qubits to execute over $1000$ entangling gates. We prepare far-from-equilibrium initial string configurations that we quench across a range of parameters to observe rich dynamical phenomena, such as the formation of gauge-invariant closed-loop excitations reminiscent of glueballs in QCD and multi-order string breaking accompanied by spontaneous matter creation. We further demonstrate experimentally that the system displays genuine $2+1$D dynamics, as evidenced by string snapshots over time that cannot be trivially mapped to $1+1$D physics. Our results demonstrate digital quantum simulations of nonequilibrium dynamics in a higher-dimensional lattice gauge theory and provide an experimentally accessible setting for phenomena related to confinement physics.

We present the first lattice QCD calculation of the next-to-leading order (NLO) hadronic vacuum polarization (HVP) contribution to the muon anomalous magnetic moment with sub-percent precision. We employ the time-momentum representation combined with the spatially summed vector correlator computed on CLS ensembles with $N_{\mathrm{f}}=2+1$ flavors of $\mathrm{O}(a)$-improved Wilson fermions, spanning six lattice spacings ($0.039$-$0.097\,$fm) and a range of pion masses including the physical value. After accounting for finite-size corrections and isospin-breaking effects, we obtain in the continuum limit $a_\mu^{\mathrm{hvp,\,nlo}} = (-101.57 \pm 0.26_{\rm stat} \pm 0.54_{\rm syst}) \times 10^{-11}$, corresponding to a total relative error of 0.6$\%$. Our result lies 1.4$\sigma$ below the estimate of the 2025 White Paper update and is two times more precise. It also shows a tension of $4.6\sigma$ with data-driven evaluations based on hadronic cross section measurements prior to the CMD-3 result.

Quantum simulations of high-energy physics in $2+1$D can probe dynamical phenomena nonexistent in one spatial dimension and access regimes that are challenging for existing classical simulation methods. For string dynamics -- relevant to hadronization -- a plaquette term is required to realize genuine $2+1$D behavior, as it endows the gauge field with dynamics and enables the propagation of photon-like excitations. Here, we realize a U$(1)$ quantum link model of quantum electrodynamics in two spatial dimensions with a tunable plaquette term on a \texttt{Quantinuum System Model H2} quantum computer. We implement, to our knowledge, the largest quantum simulation of string-breaking dynamics reported to date, on a $5 \times 4$ matter-site square lattice using $51$ qubits. The simulation uses a shallow circuit design with a two-qubit gate depth of $28$ per Trotter step and up to $1540$ entangling gates. Starting from far-from-equilibrium string configurations, we measure the probability for the string to propagate within the lattice plane and find signatures of genuine $2+1$D dynamics only when the plaquette term is present. In a resonant regime, we observe the annihilation of string segments accompanied by the production of electron--positron pairs that screen them. We further find that, only with a nonzero plaquette term, matter creation extends across the lattice plane rather than remaining confined to the initial string path. These results experimentally realize string breaking and demonstrate the emergence of dynamical gauge fields in two spatial dimensions, establishing a route to photon-like propagation in programmable quantum simulators of gauge theories.

A first-order, confinement/deconfinement phase transition appears in the finite temperature behavior of many non-Abelian gauge theories. These theories play an important role in proposals for completion of the Standard Model of particle physics, hence the phase transition might have occurred in the early stages of evolution of our universe, leaving behind a detectable relic stochastic background of gravitational waves. Lattice field theory studies implementing the density of states method have the potential to provide detailed information about the phase transition, and measure the parameters determining the gravitational-wave power spectrum, by overcoming some the challenges faced with importance-sampling methods. We assess this potential for a representative choice of Yang-Mills theory with $Sp(4)$ gauge group. We characterize its finite-temperature, first-order phase transition, in the thermodynamic (infinite volume) limit, for two different choices of number of sites in the compact time direction, hence taking the first steps towards the continuum limit extrapolation. We demonstrate the persistence of non-perturbative phenomena associated to the first-order phase transition: coexistence of states, metastability, latent heat, surface tension. We find consistency between several different strategies for the extraction of the volume-dependent critical coupling, hence assessing the size of systematic effects. We also determine the minimum choice of ratio between spatial and time extent of the lattice that allows to identify the contribution of the surface tension to the free energy. We observe that this ratio scales non-trivially with the time extent of the lattice, and comment on the implications for future high-precision numerical studies.

We develop a systematic Hamiltonian formulation of minimally doubled lattice fermions in (3+1) dimensions, derive their nodal structures (structures of zeros), and classify their symmetry patterns for both four-component Dirac and two-component Weyl constructions. Motivated by recent single-Weyl proposals based on Bogoliubov-de Gennes (BdG) representation, we argue that the corresponding single-Weyl Hamiltonians are obtained from the minimal-doubling Hamiltonians supplemented by an appropriate species-splitting mass term, and we re-examine the non-onsite symmetry protecting the physical Weyl node in terms of a Ginsparg-Wilson-type relation. We then construct a one-parameter family of deformations that preserves all the symmetries and demonstrate that, once the parameter exceeds a critical value, additional Weyl nodes emerge and the system exits the single-node regime. This indicates that in interacting theories radiative corrections can generate symmetry-allowed counterterms, so maintaining the desired single-Weyl phase generically requires "moderate" parameter tuning.