I derive the regularized formula for glue scalar density (gluon condensate) in terms of Dirac spectral density [arXiv:2509.03509], and elaborate on its uses and meaning. Particular attention is given to understanding of what this new formula reveals about the relation between glue and quark scalar densities, how it relates to IR phase, how it clarifies the distinction between anomalous and spontaneous ways of breaking symmetries, and what it says about the relation between UV and IR in QCD.

We expand the treatment of the problem of the extraction of smeared spectral densities from Euclidean correlators introduced in [Phys. Rev. D 99, 094508], providing an alternative which does not rely on the Backus-Gilbert regularization. This is possible due to the observation that the solution can be decomposed into a sum of terms, in the spirit of the singular value decomposition, where those with the largest contribution to the statistical noise happen to contribute the least to the central value of the smeared spectral density. The analysis of the systematics of the inverse problem is then shifted to finding the optimal truncation of such summation, so that the signal is saturated before the noise explodes. We scrutinise the performance and systematics of this approach either as a standalone procedure, or to complement the stability analysis required to extrapolate the unbiased result in the Backus-Gilbert regulated version of the solution.

We apply the Worldvolume Hybrid Monte Carlo (WV-HMC) method [arXiv:2012.08468] to the two-dimensional Hubbard model, which is known to suffer from a severe sign problem when the system is doped (away from half filling). We show that the method predicts physical observables with controlled statistical errors on an $8 \times 8$ lattice at temperature $T/t = 1/6.4 \approx 0.156$ and interaction strength $U/t = 8.0$ ($t$ is the hopping amplitude), for which the standard determinant quantum Monte Carlo fails.

Quantifying quantum resources for simulating the fundamental forces of Nature is sensitive to the mapping of gauge fields onto finite quantum computational architectures. When locally truncating lattice gauge theories in the irreducible representation basis, it has been proposed to further deform the theory via quantum groups. The purpose of this deformation is (1) to provide an infinite tower of finite-dimensional ($d = k+1$) groups systematically approximating the infinite-dimensional gauge links and (2) to restore the physical unitarity of a plaquette operator diagonalization procedure analytically derived from the field continuum by recontracting vertex pairs. For the SU(2)$_k$ Yang-Mills pure-gauge theory, we provide a constructive strategy of gauge-variant completions to extend this unitarity to the entire computational Hilbert space, leading to well-defined time evolution unitaries as targets for optimized circuit synthesis. Leveraging basic circuit decompositions and symmetries of the diagonalized plaquette operator, we report resource upper-bounds on the generalized-controlled-X two-qudit gates for arbitrary local truncation $d$, reducing estimates and scaling relative to the non-deformed theory by three polynomial powers from $O(d^8)$ to $O(d^5)$. Examining the stronger q-deformed gauge constraint, which softens the total flux at vertices, we show that the physical Hilbert space dimension of the deformed plaquette operator scales equivalently to its non-deformed counterpart with a constant factor $0.2563(5)$. Thus, despite affecting interactions at all scales as exemplified by the observed flux hierarchy inversion symmetry, q-deformation continues to pass scrutiny as a reliable truncation offering advantages in quantum circuit synthesis.

In this work we investigate the influence of weak acceleration on the confinement-deconfinement phase transition in gluodynamics. Our study is carried out within lattice simulation in the comoving reference frame of accelerated observer which is parameterized by the Rindler coordinates. We find that finite temperature confinement-deconfinement phase transition turns into spatial crossover in the Rindler spacetime. In other words, spatially separated confinement and deconfinement phases can coexist in the Rindler spacetime within certain intervals of temperature and acceleration. We determine the position of the boundary between the phases as a function of temperature for several accelerations and find that it can be described by the Tolman-Ehrenfest law with rather good accuracy although a minor deviation takes place. Moreover, the critical temperature of the system in the weak acceleration regime is found to remain unchanged as that of the standard homogeneous gluodynamics. Our results imply that the spatial confinement-deconfinement transition might take place in the vicinity of the Schwarzschild black hole horizon.

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.

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 Preskill (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 Worldvolume Hybrid Monte Carlo (WV-HMC) method [arXiv:2012.08468] is an efficient algorithm for addressing the numerical sign problem at moderate computational cost. It mitigates the sign problem while avoiding the ergodicity issues inherent in approaches based on Lefschetz thimbles. In this study, we apply WV-HMC to the two-dimensional Hubbard model doped away from half filling, which is known to suffer from a severe sign problem. We compute the number density and the energy density on lattices of size $6 \times 6$ and $8 \times 8$ at temperature $T/t = 1/6.4 \simeq 0.156$ and interaction strength $U/t = 8.0$, using Trotter number $N_t = 20$ (Trotter step $\epsilon = 0.32$). Our results demonstrate that WV-HMC remains effective even in parameter regimes where standard (non-thimble) determinant quantum Monte Carlo methods fail. In this work, fermion matrix inversions are performed using direct solvers, leading to a computational cost of $O(N^3)$, where $N$ denotes the number of degrees of freedom and is proportional to the spacetime lattice volume. An alternative algorithm employing pseudofermions and iterative solvers, which reduces the cost to $O(N^2)$ at the expense of careful parameter tuning, will be discussed in a separate publication.