hep-lat

Trained lattice samplers are usually judged by the ensembles they generate. Here we instead analyze the trained field-space function itself: a flow-matching velocity, a diffusion score, or a normalizing-flow action residual. We project these functions onto operator bases fixed before the fit, chosen from symmetry, exact Gaussian path limits, finite-volume modes, and gauge covariance. For two-dimensional lattice \(\phi^4\), a trained straight-flow teacher is not described by a local force basis alone. After the local transport basis, the residual separates into a zero-mode Binder component and a lowest-shell finite-\(k\) correlator component. The deflated zero-mode polynomial \(P_5(M;t)\) reduces the dominant Binder-tail component, while \(\phi^\perp_{|n|^2=1}\) reduces the finite-\(k\) correlator component; wrong-parity, off-zero-mode, and random controls do not produce the same reductions. The same projection distinguishes other sampler classes. Diffusion follows the force-resolvent ordering predicted by the free theory, reverse-KL normalizing-flow collapse appears as a forbidden odd zero-mode residual, and gauge-equivariant teachers are resolved by Wilson-loop-force tangent directions. The operator basis is model- and symmetry-dependent, but the test is common: project the trained field-space function and retain sectors that lower held-out residuals and pass the available controls.

We present a method to quantize free fermions which eliminates the doublers when implemented on the lattice in any number of dimensions and in the $m=0$ limit. The elimination of doublers is achieved by combining a second-order description of fermions, with the tools associated with non-Hermitian Hamiltonians. We identify a new Pseudo-Hermitian symmetry of the second-order fermion equation, and we identify the associated $U(1)$ symmetry which will become charge when shifted to a local gauge theory. We validated the methods numerically.

We demonstrate that conditional Masked Autoregressive Flows constitute a flexible interpolation tool for lattice QCD observables, conditioned on bare lattice parameters. As a benchmark, we use the chiral phase structure of QCD with five degenerate light quark flavours, which on coarse lattices exhibits a region of first-order chiral transitions terminating in a critical quark mass. The method successfully reproduces standard reweighting in the gauge coupling, and naturally extends to interpolation in quark mass and spatial volume, for which reweighting is computationally prohibitive or inapplicable, respectively. Once trained, the model generates samples across the full parameter space in minutes, which can be used to obtain consistent first estimates of the critical quark mass without simulating all intermediate parameter values. This offers a concrete reduction in the number of lattice ensembles required. Precision on the critical mass from learned distributions is so far prohibited by the mode-covering effect inherent to maximum-likelihood-based training, which introduces a systematic bias near first-order transitions. At the current stage, the method is well-suited for a range of practical applications: localising phase boundaries, identifying the universal scaling axes at a critical point, and accelerating informed determinations of parameter values ahead of high-precision Monte Carlo campaigns.

We study the temperature dependence of meson screening masses in two-flavour lattice QCD using dynamical truncated overlap fermions (TOF), a type of lattice chiral fermions. The screening masses for the $\pi$, $\rho$, $a_1$, $a_0$, and the sigma $(f_0)$ mesons are extracted by computing spatial correlation functions. Above the pseudocritical temperature $T_{\rm pc}$, the $\pi$ and $f_0$ screening masses become degenerate, consistent with chiral restoration. The $(\pi,f_0)$ and $(\rho,a_1)$ pairs also show the expected degeneracy. Decomposition of the $f_0$ propagator reveals that the connected contribution dominates above $T_{\rm pc}$, while the disconnected part becomes significant below $T_{\rm pc}$, explaining the reduced statistical clarity observed at low $T$. These results demonstrate that dynamical TOF simulations can capture the qualitative thermal behaviour of the scalar sector.

We present a new set of lattice QCD ensembles with four flavors of smeared M\"obius 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_{\pi^{\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.

We study the signatures of the Chiral Magnetic Effect (CME) in $SU(2)$ gauge theory with $N_f = 2$ flavours of dynamical fermions at finite temperature $T$, quark chemical potential $\mu$ and a weak external magnetic field $e B$. We consider both the correlator of the axial density and the vector current, which gives direct access to the out-of-equilibrium CME, and the correlator of two vector currents, which probes the CME indirectly via the enhancement of the longitudinal electric conductivity (Negative Magnetoresistance, NMR). We find that the CME response extracted from the vector-axial correlator exhibits a rather weak dependence on temperature and density in the quark-gluon plasma regime, and is very close to the universal value for free massless quarks. The CME is mildly suppressed at low temperatures and large densities in the hadronic phase. In contrast, the NMR behaves in a qualitatively different way across the phase diagram, and is strongly suppressed at either large densities or temperatures. The magnitude of the NMR response appears to be considerably smaller than the prediction based on the lowest Landau level calculation for free quarks. Our findings suggest that for relatively small magnetic field strengths $e B \lesssim m_{\pi}^2$ the relation between the CME and NMR might not be as direct as expected. We also do not find statistically significant indications for an enhancement of the CME strength in the vicinity of the crossover or second-order phase transition lines in the $(\mu, T)$ phase diagram.

Using forward matrix elements of local leading-twist operators, we present a determination of the isovector third Mellin moments $\left< x^2 \right>$ of nucleon unpolarized, polarized, and transversity parton distribution functions. Two lattice QCD ensembles at the physical pion mass are used, which were generated using a tree-level Symanzik-improved gauge action and 2+1 flavor tree-level improved Wilson Clover fermions coupling via 2-level HEX-smearing. Leveraging a wide set of operators, two extraction methods for the matrix elements, and the automatic inclusion of model uncertainties via bootstrapped model averages, we extract values of the third Mellin moments. This is the first direct calculation of these observables performed at the physical pion mass.

We present the strange electromagnetic form factors of the nucleon using lattice QCD with $N_f=2+1+1$ twisted mass clover-improved fermions and quark masses tuned to their physical values. Using four ensembles with lattice spacings of $a=0.080$ fm, $0.068$ fm, $0.057$ fm and $0.049$ fm, and similar physical volume, we obtain the continuum limit directly at the physical pion mass. The disconnected strange contributions are computed using high statistics two-point functions combined with stochastic noise mitigation techniques, such as spin-color dilution and hierarchical probing in the estimation of the quark loop. From the momentum dependence of the form factors, we provide the strange electric and magnetic radii, as well as the strange magnetic moment in the continuum limit.

Isospin-breaking corrections to the HVP are among the leading sources of uncertainty in the Standard Model prediction of the muon $g-2$. In recent work by the RC$^{\star}$ collaboration, we compute the intermediate window contribution for a flavour non-singlet current using two strategies to include isospin-breaking corrections: the RM123 approach and a fully non-perturbative dynamical QCD+QED simulation. In both computations, we use $C$-periodic spatial boundary conditions to ensure that locality, gauge invariance, and translational invariance are preserved throughout the calculation. At fixed lattice spacing and volume with $N_f =1+2+1$ dynamical fermions, and fully including sea-quark effects in both computations, we find that simulating the full QCD+QED distribution yields smaller uncertainties for a fixed statistics. We summarize the comparison of the two approaches and discuss the implications for future lattice QCD+QED computations.

A significant component of the cost of making predictions from lattice QCD stems from the computation of correlation functions on a given ensemble of gauge fields. This cost depends on the observable of interest and the details of its representation, including any approximation needed to estimate it. Moreover, the variance of such estimators may depend strongly on physical and kinematical parameters such as the lattice spacing, volume or separation, which gives an important insight into the costs of reaching the relevant physical limits. In these proceedings, I review some observables involving quark propagators, including both quark-line connected and disconnected Wick contractions, and discuss variance-reduction schemes based on decompositions of the quark propagators. Such strategies have already proven useful for precision physics observables and in future may help reduce the computational cost of reaching large volumes.

In this chapter we provide a pedagogical introduction to the main theoretical aspects related to topology and $\theta$-dependence in Quantum Chromo-Dynamics (QCD), and to their phenomenological relevance in the Standard Model ($\eta^\prime$ physics, neutron electric dipole moment) and beyond (strong CP problem and the axion solution). We then provide an overview of the main analytic predictions for $\theta$-dependence obtained using several different approaches (chiral effective theories, large-$N$ arguments, semiclassical methods) and their regimes of validity, as well as a selection of the most recent numerical results about QCD topology obtained via Monte Carlo simulations of the lattice-discretized theory.

This chapter, to appear in the section on QCD under extreme conditions within the Encyclopedia of Nuclear Physics, aims to provide a pedagogical introduction to the physics of quarks and gluons in the presence of high temperature, nonzero (isospin) density and strong background electromagnetic fields. Extreme conditions of these types are relevant for the description of high-energy heavy-ion collisions, neutron stars and their mergers, as well as the evolution of the early Universe in its first microsecond. Most of the existing results on this topic have been obtained by means of first-principles simulations of the discretized theory of the strong interactions, lattice Quantum Chromodynamics (QCD). This lays the focus of this review chapter, although various calculations within effective theories of QCD -- most notably chiral perturbation theory -- are also discussed. Furthermore, we provide an outlook concerning open questions and yet uncharted parameter regions within this fascinating system.

Heavy meson HQET light-cone distribution amplitudes (LCDAs) are critical for precision predictions of $B$ meson weak decays, but currently are one of dominant theoretical uncertainties that obscure interpretations of $B$ anomalies and CP-violating measurements. Building on the established HQLaMET framework, supplemented by lattice QCD calculations of the OPE moments, we present a precise lattice QCD calculation of HQET LCDAs by employing multi-ensemble simulations for continuum and physical pion mass extrapolation, quantifying comprehensive systematic errors, and validating results through OPE moment cross-validation. Details of the lattice calculations are provided in a companion paper \cite{HeavymesonDA_long_paper}. Our final results for key inverse moments (at $\mu=1$ GeV) are $\lambda_B=0.340(20)$ GeV and $\sigma_B^{(1)}=1.685(63)$, with the total uncertainty reduced by a factor of three relative to the previous analysis. These results can greatly reduce the uncertainty in the $B \to K^*$ form factors in the large-recoil region. This work resolves the long-standing bottleneck in first-principles predictions of heavy meson LCDAs, advancing precision flavor physics to new frontiers.

We discuss the calculation of the inclusive semileptonic decay for the process $B_s \to X_c \, l \nu_l$ using lattice QCD. This calculation could be decisive in understanding the long-standing tension between inclusive and exclusive determinations of the CKM matrix element, $|V_{cb}|$. In this talk, we investigate the main sources of systematic uncertainty in these decays, including the impact of Jacobi smearing at the source and sink, variations in source-sink separation, and the intrinsic uncertainties of the inclusive reconstruction method itself. In addition, we explain how we can restrict the reconstruction of the inclusive decay rate to just the excited-state contributions. This is achieved by treating the ground-state contributions as an exclusive decay with well-controlled conventional techniques. Systematic effects from the reconstruction then only affect excited-state contributions. Where these are sub-dominant, a suppression of systematic effects is expected. We show results based on Chebyshev reconstruction, which are part of a larger effort towards a first phenomenologically relevant computation of the inclusive decay rate in the continuum and infinite-volume limits.

We present a first-principles determination of heavy meson light-cone distribution amplitudes (LCDAs) from lattice QCD in the continuum limit, improving substantially on our previous pioneering study. Within the heavy-quark large-momentum effective theory (HQLaMET) framework, supplemented by lattice QCD calculations of the OPE moments, we analyze six ensembles with lattice spacings ranging from $a=0.0519-0.1053$\,fm and pion masses from $m_\pi=135.5-317.2$\,MeV, thereby enabling controlled continuum, chiral, and infinite-momentum extrapolations to the physical point. Momentum-smeared sources, hypercubic-smeared Wilson lines, and optimized interpolating operators are adopted to significantly improved signals for the nonlocal correlators. Within a unified framework, we determine both QCD LCDAs and HQET LCDAs. Our resulting QCD LCDAs of $D$ meson peak at $y\approx 0.2-0.3$, with total uncertainties below $30\%$ for $0.1<y<0.9$. The leading-twist HQET LCDA is constructed using a peak-and-tail factorization, in which the nonperturbative peak region is obtained from lattice QCD and the perturbative tail is incorporated from HQET, with the two regions combined through a model-independent Laguerre-polynomial parametrization. At $\mu=1$\,GeV, we obtain the inverse moment of HQET LCDA $\lambda_B=0.340(20)$\,GeV and first inverse-logarithmic moment $\sigma_B^{(1)}=1.685(63)$, consistent with experimental constraints and phenomenological determinations. Direct lattice calculations based on operator product expansion provide a nontrivial cross-check of the LaMET results. Final results and phenomenological impact of these results are presented in a companion paper~\cite{HeavymesonDA_short_paper}. Our results remove the single-lattice-spacing limitation of the previous study, and provide a robust determinations of heavy meson LCDAs in both QCD and HQET for next-generation heavy flavor physics.

The type IIB matrix model has been proposed as a nonperturbative formulation of superstring theory. While numerical simulations of this model are essential for probing nonperturbative effects, such as the emergence of time and an expanding 3--dimensional space, they are hindered by the sign problem. We address this using the Complex Langevin Method (CLM). Furthermore, to suppress spurious numerical artifacts that originate from large Lorentz boosts due to the Lorentz symmetry of the model, we nonperturbatively fix the Lorentz symmetry using the Faddeev--Popov procedure. We then study this model to investigate the impact of supersymmetry on the dynamical generation of (3+1)--dimensional spacetime.

We propose and analyze a new method of detecting center vortices and monopoles in lattice Yang-Mills theory. This procedure is sensitive to the intrinsic degeneracy of the center charges, which play a crucial role in how these topological objects interact and correlate with one another. Our approach is based on fixing the Maximal Abelian gauge (MAG) and decomposing the link configuration in a suitable way to look for so-called Cartan fluxes. Our discussion is general for $SU(N)$ gauge theory, but we focus our applications on the $SU(3)$ case. For the $SU(2)$ case, our proposed parametrization is equivalent to the usual one.

Quantum simulation offers a promising framework for quantum field theory calculations. Obtaining reliable results, however, requires careful characterization of systematic uncertainties. One important source is the boson truncation error, which arises from representing infinite-dimensional local Hilbert spaces with finite-dimensional ones. Previous studies have examined this problem from several perspectives. In particular, Jordan, Lee, and Prekill (arXiv:1111.3633) derived an energy-based bound applicable to generic low-energy states across a broad class of field theories. However, this approach often yields overly conservative bounds, especially at large volumes. In this work, we introduce a new methodology that significantly tightens the energy-based boson truncation bound through two complementary advances: an improved analytic derivation and a Monte Carlo-based numerical procedure. We demonstrate the method in (1+1)-dimensional scalar field theory and (2+1)-dimensional U(1) gauge theory in the dual formalism. Our approach substantially mitigates the volume dependence of the required truncation cutoff, achieving reductions nearly proportional to the volume in some cases and to the square root of the volume in others.

The sign problem obstructs the determination of the QCD phase diagram in the temperature-baryon chemical potential plane using lattice QCD. We review the sign problem in QCD and related field theories, including applications to real-time dynamics. We focus on approaches where the sign problem can potentially be solved or controlled, irrespective of its severeness. These include holomorphic extensions -- Lefschetz thimbles, holomorphic flow, contour deformations, and complex Langevin dynamics --, and the introduction of new degrees of freedom -- dual variables and the tensor renormalisation group. We also highlight directions in which machine learning approaches have shown promise. Since many methods are first tested in simpler models, we provide an outlook on their feasibility for lattice systems.

The sign problem obstructs the determination of the QCD phase diagram in the temperature-baryon chemical potential plane using lattice QCD. We review the sign problem in QCD and related field theories, including applications to real-time dynamics. We focus on approaches where the sign problem can potentially be solved or controlled, irrespective of its severeness. These include holomorphic extensions -- Lefschetz thimbles, holomorphic flow, contour deformations, and complex Langevin dynamics --, and the introduction of new degrees of freedom -- dual variables and the tensor renormalisation group. We also highlight directions in which machine learning approaches have shown promise. Since many methods are first tested in simpler models, we provide an outlook on their feasibility for lattice systems.

We calculate second- and selected fourth-order conserved-charge fluctuations in $(2+1)$-flavor QCD using M\"obius domain-wall fermions (MDWF) along a line of constant physics. Gauge ensembles were generated for two light-to-strange quark-mass ratios, $m_l/m_s=1/10$ and $1/27.4$, corresponding to heavier-than-physical and physical pion masses, respectively. For $m_l/m_s=1/10$, calculations were carried out on lattices with temporal extents $N_\tau=12$ and $16$, enabling an assessment of lattice-spacing effects at heavier pion mass. For $m_l/m_s=1/27.4$, calculations were performed at $N_\tau=12$, allowing us to study the light-quark-mass dependence down to the physical point. Below the pseudocritical temperature, second-order electric-charge, strangeness, and off-diagonal conserved-charge fluctuations are consistent with QMHRG2020 hadron resonance gas calculations. Across the crossover region, these observables rise rapidly and tend toward their Stefan--Boltzmann limits. Selected fourth-order cumulants were also computed at the physical pion mass. Although these observables are statistically more demanding, several channels with controlled uncertainties permit a first comparison with hadron resonance gas calculations.

We present preliminary results from a lattice-QCD study of the hadronic contributions to the running of the electromagnetic coupling, $\Delta\alpha(Q^2)$, and the electroweak mixing angle, $\Delta\sin^2\theta_{W}(Q^2)$. Using $N_f = 2+1+1$ HISQ ensembles at physical quark masses, we discuss the challenges posed by strong statistical correlations in the time-momentum representation and propose a spectral-reconstruction strategy to obtain controlled continuum-extrapolated results across the full energy range.

We present results for the inclusive semileptonic $D_s \to X_s \ell\bar\nu$ decay rate from lattice QCD. Chiral and continuum extrapolations are performed using gauge ensembles generated with 2+1 flavours of M\"obius domain-wall fermions. Systematic errors are fully addressed including those from the integral over all possible final states. Our results are in agreement with currently available experimental data, with an error at the few-percent level.

We propose a reconstruction of the full $(x, \xi, t)$ dependence of unpolarized isovector proton generalized parton distributions (GPDs) $H^{u-d}$ and $E^{u-d}$ from lattice QCD data in the pseudo-distribution formalism. For the first time, we extract double distributions (DDs) directly from lattice data, enforcing therefore an important property of GPDs linked to Lorentz symmetry. We use the flexible framework of multidimensional Gaussian process regression to regularize the inverse problem and present an assessment of the impact of model dependence on the systematic uncertainty. Our lattice ensemble corresponds to a pion mass $m_\pi = 358$~MeV and a lattice spacing $a = 0.094$~fm. We use larger hadron momenta, up to 2.7~GeV, and kinematic coverage compared to our previous computations and extract additional skewness-dependent moments of the GPD.

We present recent results from the FASTSUM collaboration, using anisotropic lattice QCD to study spectral properties of heavy quarkonia and open heavy flavour systems at high temperature. For heavy quarkonium, our results using a number of different methods suggest a small but significant and robust negative mass shift as well as an increasing thermal width. We present the first lattice results for masses and spectral functions of B mesons at high temperature, and preliminary results for a high-precision calculation of the static quark potential.

We present the results of continuum-extrapolated lattice simulations of quantum chromodynamics (QCD) above the crossover temperature and for unprecedentedly high baryon densities at the physical point, employing the complex Langevin equation. In particular, we determine the QCD equation of state by computing the baryon density as well as the pressure as functions of the baryon chemical potential and the temperature. Potential issues with wrong convergence of complex Langevin dynamics are under control and we indeed find agreement with previous lattice studies working at smaller chemical potentials, as well as with perturbative hard-thermal-loop calculations at high temperatures.

Motivated by recent experimental observations of the flavour-exotic $T^*_{cs0}(2870)^0$ and $T^*_{c\bar{s}0}(2900)$, we present the first lattice QCD study of coupled-channel scattering of a charm meson with a light meson in the flavour-exotic sectors at the $SU(3)_f$ flavour symmetric point. Utilising five volumes with $m_\pi \approx 700$ MeV and employing large bases of meson-meson operators, finite-volume spectra are extracted and used to constrain infinite-volume scattering amplitudes with $J^P = \{0, 1, 2, 3, 4\}^+$ via the L\"uscher formalism. In the flavour $\mathbf{6}$ sector, each $S$-wave channel considered is found to be attractive with the scattering amplitudes having an associated pole singularity on an unphysical sheet below threshold, giving six flavour-exotic poles in the energy region constrained. In $J^P = 0^+$ there is a virtual bound state and a resonance. The latter is identified with the $T^*_{cs0}(2870)^0$ and $T^*_{c\bar{s}0}(2900)$, appearing as one state in the $SU(3)_f$ flavour symmetric limit, and suggests the existence of an isospin-$\frac{1}{2}$ partner. In $J^P =1^+$ there are three poles, one of which is identified as a $J^P =1^+$ partner of the $T^*_{cs0}(2870)^0$ and $T^*_{c\bar{s}0}(2900)$, and $J^P =2^+$ contains one pole which is identified as their $J^P =2^+$ partner. Only mild interactions and no poles are seen in the $J^P = \{3, 4\}^+$ scattering amplitudes. In the flavour $\overline{\mathbf{15}}$ sector, weak interactions are observed in $J^P = \{0, 1, 2, 3, 4\}^+$ with no well-determined poles in the energy region constrained.

We present an update on the ongoing computation of the isospin-breaking effects in the Pion Decay Constant from the BMW Collaboration. The calculation is carried out with N$_f$=2+1+1 staggered quarks with a near-physical pion mass and QED$_{\text{L}}$. We give an update on the isosymmetric value and the current determination used to compute the gradient-flow scale $w_{0}$, then we present some preliminary results of the valence-valence contribution to the axial-pseudoscalar correlator for different volumes and lattice spacings. We also discuss the next steps and plans.

In this paper, we apply the deterministic quantum imaginary time evolution (QITE) algorithm to obtain the ground state of a $2+1$-dimensional pure $\mathbb{Z}_2$ lattice gauge theory. We first construct the set of Pauli operators commuting with Gauss's law constraints, generalizing a previous result. This makes the deterministic QITE gauge-invariant and reduces both the measurement and gate costs significantly without adding extra algorithm errors in the QITE. Then, the classical numerical simulation of the deterministic QITE using tensor networks is performed, and the results are compared with the density matrix renormalization group (DMRG) to evaluate the accuracy of the algorithm. Specifically, we investigate the coupling and system size dependence, and find that the deterministic QITE can achieve a relative error of less than $0.1\%$ up to a twelve-plaquette system and coupling values in a regime that we study. Furthermore, the error dependence on the number of time steps is studied and discussed.

We numerically investigate whether the center-symmetric confined phase of large-$N$ $SU(N)$ gauge theory with one adjoint Dirac fermion persists under spatial compactification on $\mathbb{R}^3 \times S^1$. To this end, we employ a partially reduced twisted Eguchi-Kawai (TEK) model on a $1^3 \times L_4$ lattice with an adjoint Wilson fermion, and measure both the Polyakov loop around $S^1$ and order parameters for volume independence in the reduced directions. For $N=36$, $L_4=2$, $b=0.30\text{-}0.46$, and $\kappa=0.03\text{-}0.16$, we find that, with periodic boundary conditions, the Polyakov loop remains near zero in the light-fermion regime as the circle size is reduced. For the modified twist, the volume-independence order parameters are also consistent with zero in the explored region, supporting the validity of the partially reduced description. These results provide numerical evidence, within the reduced-model setup and parameter range studied, for an adiabatic-continuity scenario in which the confined phase is smoothly connected between large and small circles. By contrast, with antiperiodic boundary conditions, the Polyakov loop exhibits a clear deconfinement transition. We also discuss how this scenario is compatible with the anomaly constraints of the underlying four-dimensional theory. The symmetric twist is examined as a useful comparison, although its volume-independence properties appear less robust at the present value of $N$.

We present a determination of the scalar and tensor $\Lambda\to p$ transition form factors using lattice QCD. These form factors are relevant for semileptonic hyperon decays in the presence of extensions of the Standard Model that include scalar and tensor interactions. The calculation is carried out using a gauge ensemble of twisted mass fermions at the physical pion mass, following the same strategy as our recent study on vector and axial form factors for the same transition. We provide the complete set of form factors as functions of $q^2$ employing a model-independent parametrization. We examine their impact on searches for non-standard charged-current interactions via the muon-to-electron decay-rate ratio $R^{\mu e}=\Gamma(\Lambda\to p\mu\bar\nu_\mu)/\Gamma(\Lambda\to pe\bar\nu_e)$, where scalar and tensor contributions enter linearly and are helicity-enhanced relative to the electron channel. We compare this first-principles prediction for the decay-rate ratio with recent experimental measurements, thereby enabling improved constraints on non-standard charged-current interactions.

We discuss the status and progress of recent efforts to modernize the International Lattice Data Grid(ILDG).This includes activities of the metadata and middleware workinggroups concerning deployment and operation of crucial services (user management, metadata catalogues, file catalogues) and extensions of the metadata format, which have been tailored according to the needs of the large collaborations. We also report on developments and extensions that are planned to be addressed in the foreseeable future.

The Hamiltonian formulation of lattice gauge theories plays a central role in quantum simulations of gauge theories, and understanding their spectrum and other properties is expected to become crucial in the upcoming years. The relevant Hamiltonians in this framework possess local symmetry at each lattice site and may exhibit higher-form symmetries. There are then an exponentially large number of dynamically disconnected symmetry sectors, most of which are not translation-invariant. An exponential number of dynamically disconnected sectors, i.e., Hilbert space fragmentation, can also occur in systems in which no such symmetries have been identified. In this contribution, we describe an emergent gauge symmetry that is valid only in a subset of sectors of the fragmented $S=1$ dipole-conserving spin chain. These non-invertible symmetries can label exponentially many of the model's sectors. Simulating this Hamiltonian, which is not gauge-invariant, yields an exact quantum simulation of a gauge theory.

We present a lattice QCD study of dilepton production in charmonium transitions, specifically focusing on the $1^{+-} \to 0^{-+}$ and $1^{++} \to 1^{--}$ processes: $h_c \to \eta_c \ell^+ \ell^-$ and $\chi_{c1} \to J/\psi \ell^+ \ell^-$, where $\ell = e, \mu$. The relevant hadronic matrix elements are computed using gauge field configurations generated by the Extended Twisted Mass Collaboration with $N_f = 2+1+1$ dynamical Wilson--Clover twisted-mass fermions at four lattice spacings. Simulations are performed at physical dynamical $u$, $d$, $s$, and $c$ quark masses, except for the coarsest lattice, where the lightest sea quark mass corresponds to a slightly heavier pion mass. A controlled continuum extrapolation is carried out. In the continuum limit for the $h_c$ decays, we obtain $\Gamma(h_c \to \eta_c e^+ e^-) = 5.45(19)~\mathrm{keV}$, and $\Gamma(h_c \to \eta_c \mu^+ \mu^-) = 0.635(22)~\mathrm{keV}$. For the $\chi_{c1}$ decays, we find: $\Gamma(\chi_{c1} \to J/\psi e^+ e^-)= 2.869(90)~\mathrm{keV}$, and $\Gamma(\chi_{c1} \to J/\psi \mu^+ \mu^-) = 0.1993(72)~\mathrm{keV}$. Our results for the $\chi_{c1}$ decays show good compatibility with experimental data. However, our prediction for the $h_c \to \eta_c e^+ e^- $ decay rate is approximately $3\sigma$ larger than the BESIII result. We also present predictions for the differential decay widths as functions of the dilepton invariant mass, $q^2$, and for angular observables sensitive to longitudinal transition form factors, which are inaccessible in radiative decays with real photon emission. These results constitute the first fully dynamical lattice QCD predictions for dilepton decay rates in $h_c$ and $\chi_{c1}$ charmonium transitions, including their differential distributions and angular observables. They provide benchmark predictions for future experimental studies.

We present a minimal implementation of SU($N$) pure Yang-Mills theory in $3+1$ dimensions for digital quantum simulation, designed to enable quantum advantage. Building on the orbifold lattice simulation protocol with logarithmic scaling in the local Hilbert-space truncation, we introduce further simplified Hamiltonians. Furthermore, we test simple methods that improve the convergence to the infinite mass limit, thereby removing the requirement of a large scalar mass to obtain the Kogut-Susskind Hamiltonian. For the SU(2) theory, we can cut the resource requirement further by utilizing the embedding of $\mathrm{SU}(2)\cong\mathrm{S}^3$ into $\mathbb{R}^4$. Monte Carlo simulations of the Euclidean path integral were used to benchmark the accuracy of these new analytical improvements to the theory. These results provide further support for the noncompact-variable-based approach as a practical framework for quantum simulation of non-Abelian gauge theories.

We propose an effective Hamiltonian formulation of quantum field theories using a Daubechies wavelet basis in position space. Combined with flow-equation methods of the similarity renormalization group (SRG), this approach provides an efficient framework for analyzing quantum field theories by reducing the dimensionality of the Hamiltonian and systematically decoupling degrees of freedom across scales. As an application, the free scalar field theory has been reformulated within this framework to calculate the low-lying energy spectrum of the theory. These basis elements are known to transform the free scalar field theory into a theory of coupled localized oscillators, each of which is labeled by a location and a resolution index. In this representation, the Hamiltonian is naturally organized into fixed-resolution blocks, alongside blocks associated with the interactions between different resolutions. To decouple the different resolution modes and obtain a block diagonalized Hamiltonian with each block associated with a fixed resolution, the flow equation approach of SRG is applied. Finally, we demonstrate that with increasing resolution, the low-energy spectrum can be extracted from the effective lowest-resolution block of the Hamiltonian, leading to a significant reduction in computational cost.

Quantum Chromodynamics (QCD) at finite density is most often formulated on the lattice as a grand canonical ensemble. Since lattice QCD has a complex action problem at finite baryo-chemical potential ($\mu_B$), its results at finite density are indirect: e.g. in the form of a set of expansion coefficients. In contrast, the canonical formulation offers direct results for integer-valued net-baryon number. In this work we present for the first time results in the canonical formulation with physical quark masses. To this end we use a high statistics finite-volume lattice ($16^3\times8$) data set that we generated at $\mu_B=0$ with our 4HEX staggered action. We extend the canonical ensemble to non-integer net-baryon number and connect the results back to the grand canonical ensemble. Unlike reweighing to real $\mu_B$, this method can also be used with rooted staggered quarks. For densities where the sign problem can be overcome by brute force computing power, this scheme provides lattice QCD results (e.g. for pressure, baryon density) directly, without relying on any extrapolation in the baryo-chemical potential. In this work we chart the phase diagram by studying bulk thermodynamic observables, which we show to be feasible up to $\mu_B\approx500$~MeV.

We present a lattice QCD calculation of the electric polarizability of the charged kaon using a four-point function approach, which is the Euclidean analog of low-energy Compton scattering. In the case of the charged kaon, the polarizability is separated into an elastic (Born) term, determined from the charge radius extracted via the kaon electromagnetic form factor, and an inelastic (non-Born) term obtained from the time-integrated difference of four-point correlation functions. Our study employs 500 configurations of Wilson quenched $24^3\times 48$ lattices, and we compute connected diagrams as a proof of principle. From this analysis, we obtain values for the charged kaon electric polarizability of $\alpha_E = (0.988 \pm 0.534) \times 10^{-4}\;\mathrm{fm}^3$ as well as $\langle r_E^2\rangle =0.3303\pm 0.0028\;\mathrm{fm}^2$ for the squared kaon charge radius, after extrapolation to the physical pion mass. The study demonstrates the applicability of the four-point function framework to strange mesons, extends previous four-point function polarizability studies, and provides a foundation for future calculations with increased statistics, dynamical fermions, and improved control of systematic uncertainties.

We present a lattice QCD study of heavy baryons containing charm and bottom quarks, with particular emphasis on the relativistic treatment of all valence quarks. We use $N_f=2+1+1$ HISQ ensembles at the physical point to compute ground-state energies of spin-$3/2^+$ baryons, including singly-, doubly-, and triply-heavy charmed and bottom baryons. This work represents the first investigation of heavy baryons using fully relativistic bottom quarks.

We introduce the \emph{flavoured lattice Schwinger model}, a $(1{+}1)$-dimensional $U(1)$ lattice gauge theory in which the fermion doubling problem is resolved by staggering a $\mathbb{Z}_{2}$ flavour degree of freedom rather than staggering chirality. Unlike all standard approaches, the flavoured construction preserves an exact axial $U(1)$ symmetry at finite lattice spacing. We derive the continuum limit, showing the model reduces to two copies of the massless Schwinger model labelled by $\alpha\in\{0,1\}$. The central result is that the flavoured construction admits a well-defined, regularized, gauge-invariant lattice axial charge $Q_{G}^{A}$ with chiral anomaly equation $\langle dQ_{G}^{A}/dt\rangle = -(2g/\pi)\int dx\,\langle E(x)\rangle$ in the continuum limit, derived as a direct dynamical consequence of minimal gauge coupling at finite lattice spacing. Restricting to the $\alpha=0$ sector recovers the standard single-flavour result. We further show that spatial separation of the flavour sectors can be realised as a helical edge states living on the boundaries of a ribbon shaped $(2{+}1)$-dimensional Bernevig--Hughes--Zhang topological insulator. This provides a bulk-boundary picture solution to fermion doubling and allows the chiral anomaly to be put on the lattice for a single flavour.

We propose a variant of the model-independent method for determining meson charge radii from spatial moments of correlation functions on the lattice. Traditional determinations based on fits to the momentum transfer squared dependence of form factors are subject to systematic uncertainties arising from the choice of fit ansatz. By contrast, model-independent methods based on spatial moments provide a useful framework for determining the slope of the form factor without assuming its functional form. Recently, Feng et al. proposed a model-independent method, which drastically suppresses the finite-volume effect in the charge radius coming from higher-order contributions of the expansion of the form factor with respect to the momentum transfer squared. In this work, we introduce an auxiliary function of the momentum transfer squared and reformulate the method in terms of its product with the form factor, rather than the form factor itself, thereby further suppressing higher-order contributions, notably in cases of small volume and large radius. In particular, we investigate quadratic and logarithmic forms as practical choices for this auxiliary function. Applying this method to mock data based on a monopole form factor, as well as to actual lattice QCD data using $N_f=2+1$ gauge ensembles at $m_\pi \simeq 0.5$ and $0.3$ GeV, we find that it reduces residual finite-volume effects and provides an effective framework for meson charge radius determinations.

Hilbert space fragmentation refers to exponential growth in the number of dynamically disconnected Krylov sectors with system size. It is taken as evidence of ergodicity breaking, since conventional symmetries generate at most a polynomial number of sectors. However, we demonstrate that generalized symmetries can fragment the Hilbert space. Models with higher-form, subsystem, and gauge symmetries can have exponentially many symmetry sectors. We further prove that non-invertible symmetries can induce additional fragmentation within individual symmetry sectors. Fragmentation in several known models arises from generalized symmetries, and the presence of exponentially many Krylov sectors therefore does not by itself imply ergodicity breaking. Finally, we show that disorder free localization arises naturally from Krylov-restricted thermalization when sectors lack translation invariance, requiring neither ergodicity breaking nor gauge symmetry.

In our recent work [1] on lattice QCD calculation of the baryon leading-twist LCDAs within the framework of LaMET, a novel hybrid renormalization scheme is implemented for octet baryon quasi-DAs, yielding reliable results across both perturbative and non-perturbative regions. The numerical simulations are performed using CLQCD ensembles with $N_f = 2+1$ stout-smeared clover fermions and a Symanzik-improved gauge action. Calculations are carried out at three lattice spacings, $a = {0.052, 0.077, 0.105}$ fm. After renormalization, the linear divergences inherent in quasi-DAs are effectively removed, leading to smooth and well-behaved continuum coordinate-space distributions. These results demonstrate the viability of hybrid renormalization frameworks for light-baryon quasi-DAs and provide a robust foundation for future LaMET-based determinations of baryon LCDAs.

In this review I summarize how machine learning can be used in lattice gauge theory simulations and what ap\-proaches are currently available to improve the sampling of gauge field configurations, with a focus on applications in four-dimensional SU(3) gauge theories. These include approaches based on generative machine-learning models such as (stochastic) normalizing flows and diffusion processes, and an approach based on renormalization group (RG) transformations, more specifically the machine learning of RG-improved gauge actions using gauge-equivariant convolutional neural networks. In particular, I present scaling results for a machine-learned fixed-point action in four-dimensional SU(3) gauge theory towards the continuum limit. The results include observables based on the classically perfect gradient-flow scales, which are free of tree-level lattice artefacts to all orders, and quantities related to the static potential and the deconfinement transition.

The complex Langevin approach is a promising method for the numerical treatment of systems with a sign problem, for which conventional lattice field theory techniques based on importance sampling cannot be applied. However, complex Langevin dynamics may fail to converge in some cases and converge to a wrong limit in others, motivating the development of various diagnostic tools over the years to assess the correctness of given simulation results. This work aims at providing a systematic comparison between the most prominent such correctness criteria. In particular, the main goal is to contrast their applicability, ease of use, and - most importantly - their predictive power. To this end, four simple but nontrivial models are considered and the criteria applied to each of them. The obtained conclusions are expected to carry over to more realistic theories as well.

We determine the interaction potential of a solitonic dipole in the singlet state, modeled as an SU(2) field, using improved lattice simulations of two stationary solitons at varying separations. The potential is extracted from the energy of two-soliton configurations as a function of distance. At large separations, the interaction reproduces the classical Coulomb potential quantitatively up to an energy shift $\delta E_\infty\approx 9\;\text{keV}$ of the fitted asymptotic constant relative to $2m_ec_0^2$, assumed to be related to limited numerical precision on the lattice. At shorter distances, deviations from the Coloumb potential of point-like charges appear, that are in qualitative agreement with the asymptotic formula of perturbative Quantum Electrodynamics, reflecting the running of the fine-structure constant, with the inverse fine-structure constant ($\alpha^{-1} \approx 137$) reproduced.

The $U(1)_A$ symmetry of the massless QCD Lagrangian is explicitly broken by the axial anomaly, but it may be effectively restored at finite temperature. Determining the temperature at which this occurs is important for understanding the chiral transition and the structure of the QCD phase diagram. A commonly used probe of effective $U(1)_A$ restoration is the degeneracy of flavour non-singlet pseudoscalar and scalar susceptibilities. Using anisotropic lattice QCD ensembles with Wilson-clover fermions generated by the \textsc{Fastsum} collaboration, we study this degeneracy through hadronic correlation functions over a wide range of temperatures. The fine temporal resolution of our Generation 3 ensembles allows us to determine the temperature at which the pseudoscalar and scalar channels become degenerate. We find evidence for the effective restoration of $U(1)_A$ symmetry at $T_{U(1)_A}=319(22)$ MeV, well above the chiral crossover temperature.

Sampling lattice field theories near criticality is severely hindered by critical slowing down, which makes standard Markov chain methods increasingly inefficient at large lattice volumes. We introduce a multiscale generative sampler, inspired by renormalization-group ideas, that models the Boltzmann distribution through a coarse-to-fine hierarchy across length scales. At each level, a conditional Gaussian mixture model captures the main local dependence of newly introduced variables on the already-sampled coarse field, while a masked continuous normalizing flow refines the remaining conditional structure. Coarse levels encode the dominant long-wavelength modes, and finer levels progressively add short-distance fluctuations. In addition, because the architecture preserves coarse fields exactly during refinement, it provides exact restriction maps at no additional computational cost and directly enables unbiased Multilevel Monte Carlo (MLMC) variance reduction. For the two-dimensional scalar $\phi^4$ theory at criticality, the method achieves integrated autocorrelation times orders of magnitude smaller than Hybrid Monte Carlo (HMC) on large volumes, maintains high importance-sampling efficiency relative to other generative baselines, and reproduces unbiased physical observables in statistical agreement with long HMC simulations.

Relativistic wavefunctions of nucleon excitations are scrutinised to understand their node structure and the underlying role of local interpolating fields in generating the nucleon spectrum. In addressing quark model perspectives, approximately 4000 propagators are employed on the heaviest PACS-CS ensemble at $m_\pi \simeq$ 702 MeV. We examine the ground and four lowest-lying excited states at zero momentum for both positive- and negative-parity spectra, where the proton's d-quark wavefunction is calculated about the two u quarks at the origin. This is achieved using two local interpolating fields that each carry the quantum numbers of the nucleon but with differing spin-flavour structures, one of which vanishes in the nonrelativistic limit. We find that two distinct types of wavefunction nodes are manifest: "superposition nodes" formed through a linear combination of interpolating fields, and novel "built-in nodes" that are fundamentally built in to the s-wave Dirac components of an individual interpolating field. These are investigated qualitatively through visualisations in the form of both volume and surface renderings, and quantitatively by the calculation of radial wavefunctions. Combined, these findings build a comprehensive picture of the single-particle nucleon spectrum and how its properties derive from fundamental lattice operators.

We study the duality of lattice Maxwell theory in the modified Villain formulation, employing an ultra-local action with a theta term. Although this action is known to become non ultra-local through the Poisson resummation formula, we show that this non ultra-locality can be removed by incorporating a non-local transformation procedure into the definition of the S-transformation. As a result, the ultra-local action with a theta term exhibits an exact SL(2,Z)-duality. We further analyze the SL(2,Z)-structure of Wilson and 't Hooft loops, demonstrating that they transform properly up to a nontrivial phase factor arising from the nontrivial self-linking of the loops. This effect originates from the non-local transformation procedure in the S-transformation. Remarkably, the resulting SL(2,Z)-structure closely resembles that of non-spin Maxwell theory.

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.

We comment on a recent manuscript "Ether of Orbifolds" by Henry Lamm. In the first version, it was mistakenly claimed that the orbifold lattice Hamiltonian is not gauge invariant, and a quantity $\epsilon_g$, which has nothing to do with a non-existent "gauge violation", was introduced. The scaling of this $\epsilon_g$ was used to claim a huge simulation cost. In fact, $\epsilon_g$ characterizes the shift of the effective lattice spacing -- because, in the orbifold lattice formulation, the lattice is generated dynamically from the vacuum expectation value of the complex matrices. In the second version, the claim about the gauge symmetry was partially corrected, based on our comments. However, $\epsilon_g$ is still mistakenly interpreted as a measure of "departure from SU($N$)", inconsistently with the foundational results by Kaplan, Katz, and \"{U}nsal, and also by Arkani-Hamed, Cohen, and Georgi. This interpretation plays a central role in sustaining the argument introduced in the first version.

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 consider the problem of identifying the onset of the conformal window for QCD with $N_f$ massless flavors in the fundamental representation, and propose a new effective method to determine it from lattice simulations. This method is based on the investigation of the so-called Roberge-Weiss transition temperature $T_{RW}$, which is encountered at specific values of the imaginary baryon chemical potential, and can also be interpreted as the inverse of the critical spatial size at which charge conjugation is spontaneously broken in a finite box. Since $T_{RW}$ corresponds to a genuine phase transition for any value of the quark masses, it is a well-defined quantity; we argue that the critical $N_f$ at which $T_{RW}$ vanishes in the chiral limit coincides with the onset of the conformal window. We implement our proposal by investigating QCD with $N_f = 8$ flavors, discretized via stout improved staggered fermions and the tree-level improved Symanzik pure gauge action, at Euclidean temporal extents $N_t = 8, 10, 12, 16, 24$. In this case, we find evidence that $T_{RW}$ already vanishes in the chiral limit, indicating that $N_f = 8$ is already in the conformal window.

We report our new lattice results for the sphaleron rate calculated within a thermal effective field theory of soft SU(N) gluons whose momenta are below the magnetic scale, where $N=2,3$, for a wide range of temperatures spanning from $0.6$-$10^{15}$ GeV at sufficiently large volumes. Comparing these results with sphaleron rates in a non-thermal SU(N) plasma where the infrared gluons are over-occupied, we estimate the typical thermalization time for these ultra-soft soft gluons during the early stages of reheating after inflation. We also calculate the thermal production rate of relativistic axions due to these non-perturbatively interacting soft gluons which shows a significant deviation from its perturbative estimate even at the electroweak scale.

We report our new lattice results for the sphaleron rate calculated within a thermal effective field theory of soft SU(N) gluons whose momenta are below the magnetic scale, where $N=2,3$, for a wide range of temperatures spanning from $0.6$-$10^{15}$ GeV at sufficiently large volumes. Comparing these results with sphaleron rates in a non-thermal SU(N) plasma where the infrared gluons are over-occupied, we estimate the typical thermalization time for these ultra-soft soft gluons during the early stages of reheating after inflation. We also calculate the thermal production rate of relativistic axions due to these non-perturbatively interacting soft gluons which shows a significant deviation from its perturbative estimate even at the electroweak scale.

The extraction of the static quark-antiquark potential from lattice QCD suffers from the poor signal-to-noise ratio of Wilson loops at large Euclidean times. To overcome this, smearing methods or the Coulomb gauge are used to improve the ground-state overlap with respect to the straight Wilson line trial state within the Wilson loop. To find excited states, complicated shapes are introduced to generate specific quantum numbers. Here, we introduce a neural-network parametrization of trial states, constructed with gauge-equivariant layers and optimized with a loss function that favors ground and excited states. In the quenched theory, we automatically obtain the interpolators for the ground and excited states.

Lattice QCD (LQCD) calculations play a key role in the establishment of flavor anomalies. One of the most recent advancements in LQCD to this end has been the publication of several calculations of the $B\to D^\ast\ell\nu$ form factors, but despite all the anticipation, the LQCD results have been unable to give a final answer to the questions it was destined to answer. In this work I briefly review what is the current status of heavy-to-heavy and heavy-to-light semileptonic decays calculations in LQCD, and what we can expect for the near and not-so-near future.

We present the status of calculations of the form factors of the most relevant heavy-to-heavy and heavy-to-light decay channels. Using seven $N_f = 2+1+1$ HISQ ensembles, with lattice spacings ranging from 0.15 fm down to 0.06 fm, we calculate the form factors of the decays, including correlations among them. More than half of our ensembles feature physical pion masses, and the heavy quarks are simulated at their physical masses using the Wilson-clover action with the Fermilab interpretation. Even though we have recently seen huge qualitative and quantitative leaps in the characterization of heavy-to-heavy decays, these advances have failed to translate into improvements for the inclusive vs exclusive question, or the matter of the Lepton Flavor Universality ratios. In particular, in the $B\to D^{\ast}\ell\nu$ channel, the current situation of the lattice-QCD form factors is far from clear. Further, the latest lattice-QCD results on the heavy-to-light form factors display unexplained tensions that must urgently be resolved. The work presented here is an attempt to address these issues.

We present a first real-time study of hadronic scattering in a $(1+1)$-dimensional SU(2) lattice gauge theory with fundamental fermions using tensor-network techniques. Working in the gaugeless Hamiltonian formulation, we investigate scattering processes across sectors of fixed global baryon number $B = 0, 1, 2$, corresponding respectively to meson--meson, meson--baryon, and baryon--baryon collisions. At strong coupling, the $B = 0$ and $B = 2$ channels exhibit predominantly elastic dynamics closely resembling the U(1) Schwinger model. The mixed $B = 1$ sector displays qualitatively new behavior: meson and baryon wavepackets become entangled during the collision, with the slower state becoming spatially delocalized while the faster one propagates ballistically. We characterize these processes through local observables, entanglement entropy, and the information lattice.

Non-equilibrium properties of strongly interacting gauge theories are often intractable with classical simulation methods. Due to recent developments of quantum simulations, studies of their properties in two spatial dimensions are becoming accessible. By demonstrating the existence of an approximate spectrum-generating algebra for a pure gauge plaquette ladder, we predict and verify the existence of Quantum Many-Body Scars in spin-1 Quantum Link Models. The analysis of the model is facilitated by a dualization process that maps the original gauge theory to a constrained spin chain. Was it not for the constraint, the system would have an exact spectrum-generating algebra. We propose a set of observables for diagnosing an approximate spectrum-generating algebra, which is expected to guide quantum simulators toward interesting physical regimes.

Simulating lattice gauge theories on quantum computers presents unique challenges that drive the development of novel theoretical frameworks. The orbifold lattice approach offers a scalable method for simulating SU($N$) gauge theories in arbitrary dimensions. In this work, we present three improvements: (i) two new simplified Hamiltonians, (ii) an encoding of the SU(2) theory with smaller number of qubits, and (iii) a reduction in the requirement for large scalar masses to reach the Kogut-Susskind limit, achieved via the inclusion of an additional term in the Hamiltonian. These advancements significantly reduce circuit depth and qubit requirements for quantum simulations. We benchmarked these improvements using Monte Carlo simulations of SU(2) in (2+1) dimensions. Preliminary results demonstrate the effectiveness of these developments and further validate the use of noncompact variables as a promising framework for scalable quantum simulations of gauge theories.

High-precision determinations of Cabibbo-Kobayashi-Maskawa (CKM) matrix elements are essential probes of physics Beyond the Standard Model (BSM). Current precision tests show a deficit in the first row unitarity relation. At the current level of precision, the only relevant CKM matrix elements that contribute to this test are $|V_{ud}|$ and $|V_{us}|$. Without resorting to nuclear inputs, they can be extracted from the combination of the experimental decay width of kaon and pion leptonic decays, along with the theoretical calculation of their decay constants; combined with the decay width of semileptonic kaon decays, with the computation of the corresponding form factor at zero momentum transfer. We review current efforts by the Fermilab Lattice and MILC collaborations towards a correlated analysis of the lattice inputs needed for this test using Highly Improved Staggered Quarks (HISQ) on the $N_f=2+1+1$ MILC configurations along with Staggered Chiral Perturbation Theory (SChPT) as a functional form for the chiral-continuum limit.

We study chimera baryons, fermion bound states composed of two (hyper)quarks transforming in the fundamental and one in the antisymmetric representation of a non-Abelian gauge group. While in QCD they coincide with ordinary baryons, in composite Higgs models (CHMs) with top partial compositeness, spin-1/2 chimera baryons serve as partners of the top quark and are responsible for its large mass. We perform non-perturbative lattice calculations of the low-lying spectrum of the chimera baryons, in a specific realization of CHMs based on a Sp(4) gauge theory. In the quenched approximation, we present the numerical results in the continuum and massless limits. Then, for dynamical fermions, we measure the spectrum and matrix elements by employing a newly developed spectral density analysis for several choices of the lattice parameters.

Signatures of symmetric mass generation (SMG) have recently been reported in lattice QCD calculations employing staggered fermions. We discuss the general criteria for SMG, and demonstrate that these conditions are indeed met by the staggered fermion action, in particular by the continuum action derived by Lee and Sharpe. We propose possible RG flow around the SMG transition, guided by the numerical results. We also point out that the Goldstone tetraquark meson states provide a phenomenological signature of the "type-II" SMG phase.

We investigate the symmetry structure of the $3+1$ D staggered fermion Hamiltonian and its implications for anomalies. Since the spin and flavor degrees of freedom of Dirac fermions are distributed over the lattice, in addition to the standard on-site mass term, the staggered fermion system also admits one-, two-, and three-link bilinear terms within a unit cube as local, charge conserving mass terms with different spin and flavor dependence. We identify the spin flavor structures of all those bilinear mass terms and determine the symmetries preserved by each of them. Among them, one of the one-link mass terms preserves a larger residual symmetry associated with conserved charges that generate the Onsager algebra. Motivated by this structure, we consider a kink profile of the one-link mass and analyze the resulting domain-wall system. In the low-energy limit, the $3+1$ D bulk becomes gapped, while two-flavor massless Dirac fermions appear as localized modes on the $2+1$ D domain wall. We show that the bulk conserved charges act on the wall as generators of a flavor $\mathrm{SU}(2)$ symmetry, and that no symmetric mass gap is allowed for the boundary theory when this $\mathrm{SU}(2)$ symmetry and space reflection symmetry are both imposed. This realizes the parity anomaly of the boundary theory and shows that the boundary flavor symmetry and anomaly descend from the ultraviolet staggered-fermion Hamiltonian rather than emerging only in the infrared.

In theories with topological sectors, such as lattice QCD and four-dimensional SU(N) gauge theories with periodic boundary conditions, conventional update algorithms suffer from topological freezing due to large action barriers separating distinct sectors. With appropriately constructed bias potentials, Metadynamics and related enhanced sampling techniques can mitigate this problem and significantly reduce the integrated autocorrelation times of the topological charge and associated observables. We test strategies to accelerate the buildup of bias potentials and the possibility of extrapolating potentials from small to large volumes. We also investigate the effectiveness of orthogonal algorithmic improvements, such as longer HMC trajectories and HMC variants, which may benefit conventional simulations as well.

We consider a quenched SU(2)$\times$U(1) gauge Higgs theory on the lattice, coupled to a static vector-like fermion which, in this case, is in the same gauge group representation as the Higgs field. Physical (i.e. locally gauge invariant) electrically charged and electrically neutral states of matter particles in the electroweak theory were described decades ago, but those constructions do not exhaust all the possibilities, and new types of electrically charged/neutral states, orthogonal to former constructions, are described here. The difference has to do with how the static source, which by itself does not create a physical state, is dressed by dynamical fields. We find that, unsurprisingly, the neutral static fermion is much lighter than any of the charged fermion states. But a lattice study of the propagation of the charged fermion states indicates the existence of (at least) two particle states with different masses in charged particle spectrum.

The static energy is an excellent observable for extracting the strong coupling $\alpha_s$ on the lattice. For short distances, the static energy can be calculated both on the lattice using Wilson line correlators, and with perturbation theory up to three loop accuracy with leading ultrasoft log resummation. Comparing the perturbative expression and lattice data allows for precise determination of $\alpha_s$. We present early results for 1-loop lattice perturbation theory improvement of the Wilson loop and show how it improves the $\alpha_s$ extraction. We present a preliminary reanalysis of the TUMQCD (2+1)-flavor QCD data.

Reliable approximations for correlation functions at intermediate and strong coupling remain hard to obtain for general quantum field theories. Perturbative expansions are often asymptotic or have a finite radius of convergence, which limits their applicability beyond weak coupling. Here we combine weak- and strong-coupling expansions and propose to use two-point Pad\'e schemes to construct approximants. For lattice $\phi^4$ theory, we show that this two-point interpolation strategy yields accurate global approximations to the two-point correlation function across broad coupling regimes and compares favorably with standard one-point resummation methods. We also provide heuristic explanations for the observed convergence behavior and discuss the practical range of validity of the approach.

We present a lattice quantum chromodynamics calculation of the $x$-dependent isovector quark helicity parton distribution function (PDF) of the proton in the large momentum effective theory (LaMET) framework. Through operator product expansion (OPE) we also extract the $\tilde{d}_2$ moment of the twist-3 PDF $g_T(x)$ for the first time in the $\overline{\rm MS}$ scheme, which is proportional to the average color Lorentz force experienced by the quark in the proton. This calculation is performed on a lattice of spacing $a$ = 0.076 fm at physical quark masses. The quasi-PDF matrix elements are measured in proton states boosted to momenta $P_z=\{0, 0.25, 1.02, 1.53\}$ GeV. We first extract the lowest few helicity PDF moments from the renormalization-group (RG) invariant ratios of the matrix elements with OPE. Combined with the matrix elements relevant for $g_T(x)$, we obtain $\tilde{d}_2^{u-d}(2\ {\rm GeV})=0.0024(46)$ at next-to-leading order in $\overline{\rm MS}$. Then, the helicity quasi-PDF matrix elements are renormalized in the hybrid scheme with linear renormalon resummation and Fourier transformed to the $x$-space after an asymptotic extrapolation. The quasi-PDF is perturbatively matched to the $\overline{\rm MS}$ PDF with RG and threshold resummations at next-to-leading power and next-to-next-to-leading logarithmic accuracies. After resummations, we determine the PDF in the region $x\in[0.25,0.75]$ with controlled systematic uncertainties. The end-point regions are then parameterized, combined with the LaMET prediction at moderate $x$, and fitted to the short-distance matrix elements in coordinate space.

The orbifold lattice has been proposed as a route to practical quantum simulation of Yang--Mills theory, with claims of exponential speedup over all known approaches. Through analytical derivations, Monte Carlo simulation, and explicit circuit construction, we identify compounding costs entirely absent in Kogut--Susskind formulations: a mass-dependent Trotter overhead that scales as $m^4$, non-singlet contamination that grows as $m^2$ and worsens with penalty terms, and a mandatory mass extrapolation. Monte Carlo simulations of SU(3) establish a universal scaling: the continuum limit forces $m^2 \propto 1/a$, binding the Trotter step to the lattice spacing through a cost unique to orbifolds. For a fiducial $10^3$ calculation, the orbifold is $10^4$--$10^{10}$ times more expensive than every published alternative. These results indicate that the claimed computational advantages do not at present survive quantitative scrutiny.

Hadron structure calculations using lattice Quantum Chromodynamics (QCD) have advanced significantly in recent years. Results for charges, form factors, and lower Mellin moments can be obtained to high precision, generalized parton distributions can now be computed either directly or reconstructed from moments, and transverse-momentum-dependent distributions can be accessed through direct lattice calculations. Together, these quantities provide detailed and complementary insights into the internal structure of hadrons. These theoretical developments are highly relevant to the experimental program of the Electron-Ion Collider (EIC) and of other facilities. We review the most pertinent lattice QCD results for hadron structure that inform the EIC scientific agenda, with particular emphasis on the pion, kaon, and nucleon.