We investigate generative diffusion models as denoising samplers for two- and three-dimensional lattice $\phi^4$ theory across the symmetric, near-critical, and broken phases. Validated against ensembles generated by Fourier-accelerated HMC combined with Wolff cluster updates, the reverse-SDE sampler reproduces scalar observables and the momentum-space propagator $G(|k|)$, with residual bias concentrated in the zero-mode and, in three dimensions, the action density. We introduce two local diagnostics and an HMC-referenced effective sample size (ESS), which probe the learned drift directly, through a Metropolis-adjusted Langevin acceptance rate, and through observable-level bias and variance. Exploiting a fully convolutional architecture with weights shared across different volumes ($V=L^D$), we show that cross-volume training transfers to unseen sizes, matching or slightly improving in-distribution training in the two-dimensional symmetric and broken phases. A three-dimensional model trained on $L \in \{4, 8, 16, 32\}$ reproduces the propagator and most scalar observables at the unseen lattice size $L = 64$ across the phase diagram, with the residual susceptibility excess in the broken phase as the main exception, and improves several critical observables relative to in-distribution $L = 64$ training. This establishes cross-volume generalization as a viable mechanism for large-volume sampling, and the score learned from many cheap small-lattice configurations transfers to the target volume without retraining.

We develop in this paper the Gram/Wishart/Stiefel formulation of the \(N=2\), large--\(d\) planar endpoint theory of the BFSS/BMN matrix quantum mechanics on the lattice, obtained in our previous work. In this formulation, the endpoint degrees of freedom are reorganized into rank--two Wishart eigenvalues and relative Stiefel angular variables. This allows the holonomy invariants \(A\), \(B\), and \(R^2=A^2+B^2\) to be analyzed directly in terms of radial and angular Gram data. A central point is the large-\(R\) aligned asymptotics of the holonomy potential. Its universal linear contribution \(-A\) is absorbed into the Gaussian sector, producing the shifted mass parameter \((\alpha_\Lambda)_{\rm eff}=\alpha_\Lambda-1/2\). In the Gram/Wishart/Stiefel variables, the exact \(O(2)\) angular integral encodes this shifted sector in a rank--two Bessel kernel. The pure \(-A\) theory, which is exactly solvable in Cartesian variables, then fixes the leading Bessel/HCIZ structure: its exponential part selects the aligned configuration, while its prefactor removes the spurious doubled Wishart entropy. We then apply this structure to the transverse \(B\)-type expansion and its non-polynomial toy completion. Finite polynomial truncations lead to an apparent large--\(d\) perturbativity bound incompatible with the continuum limit, but this bound is shown to be an artifact of truncation. After summing the local transverse completion and balancing the compensating \(+A\) term, the Wishart saddle is recovered with the physical shifted mass. The resulting continuum behavior reproduces the universal \(-2d\) contribution of the \(D_\Lambda\)-channel, while the genuinely anisotropic \(\beta_\Lambda\)-channel lies outside the scope of a pure transverse \(B\)-type description.

Dilaton effective field theory (dEFT) can be employed to analyze lattice data in gauge theories that lie in close proximity of the lower edge of the conformal window. Under special conditions, we show that it can be used as a diagnostic tool to distinguish near-conformal, yet confining, theories from infrared conformal ones. We demonstrate this efficacy by analyzing two sets of lattice measurements taken from the literature. For the $SU(3)$ theory coupled to $N_f=8$ Dirac fermions transforming in the fundamental representation, our analysis favors confinement. For the $SU(2)$ theory with $N_f=1$ adjoint fermion, our fits favor infrared conformal behavior. We discuss future lattice measurements, and analysis refinements, that can further test this framework.

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.

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.

We study the effects of uncorrelated quenched disorder to the phase diagram and continuous transitions of three-dimensional lattice ${\mathbb Z}_2$ gauge Higgs models. For this purpose, we consider two types of quenched disorder, associated with the sites and plaquettes of the cubic lattice. In both cases, for sufficiently weak disorder, the phase diagram remains similar to that of the pure system, showing two different phases (one of them being a topologically ordered phase), separated by two different continuous transition lines. However, the quenched disorder changes the universality classes of the critical behaviors along some of the transition lines. The random-plaquette disorder turns out to be relevant along the topological ${\mathbb Z}_2$ gauge transition line, so the critical behaviors belong to the different random-plaquette $\mathbb{Z}_2$ gauge (RP${\mathbb Z}_2$G) universality class with length-scale exponent $\nu=\nu_{\rm rp}\approx 0.82$; on the other hand, it turns out to be irrelevant along the other Ising$^\times$ transition line (a variant of the Ising transitions with a gauge-dependent order parameter), leaving unchanged its asymptotic critical behaviors with $\nu=\nu_{\cal I}\approx 0.63$. The random-site disorder leads to a substantially different scenario: it destabilizes the Ising$^\times$ critical behaviors of the pure model, changing them into those of the randomly-dilute Ising$^{\times}$ (RDI$^{\times}$) universality class with $\nu=\nu_{\rm rdi}\approx 0.68$, while the critical behaviors along the other ${\mathbb Z}_2$ gauge topological transition line remains stable, with $\nu=\nu_{\cal I}\approx 0.63$.