We present a study of tuning of the anisotropic highly improved staggered quark (aHISQ) action on pure gauge ensembles with the renormalized anisotropy ranging from 1 to 8. We discuss multiple gradient flow schemes for tuning the gauge anisotropy and comment on what scheme may be optimal for anisotropic simulations. Next, we compare tuning of the fermion anisotropy for the naive staggered and aHISQ actions. Finally, we study the dependence of the staggered pion taste mass splittings on anisotropy for the two actions and develop an empirical model that captures the main features of the aHISQ spectrum. We observe qualitatively different behavior of the naive and aHISQ taste spectrum with anisotropy.

We report the first complete lattice-QCD determination of the leading-twist light-cone distribution amplitudes (LCDAs) of the $\Lambda$ baryon, obtained as full two-dimensional functions of the valence-quark momentum fractions. The calculation employs large-momentum effective theory to relate the light-cone amplitudes to equal-time nonlocal three-quark matrix elements of boosted $\Lambda$ baryons. Controlled physical extrapolations to the continuum, physical pion mass, and infinite momentum, together with hybrid renormalization, large-$\lambda$ extrapolation, and perturbative matching, yield the three leading-twist LCDAs $V$, $A$, and $T$. Using the lattice-determined LCDAs in place of the asymptotic form, we find an $\mathcal{O}(10\%)$ shift in the $\Lambda$ electromagnetic form factor at perturbative scales, demonstrating that the full two-dimensional LCDAs, rather than only their asymptotic shapes or lowest moments, are required for precision baryonic phenomenology. This work, together with the companion paper [1] detailing the baryon-LaMET framework, provides the first complete multi-dimensional $x$-dependent baryon LCDAs from first principles and establishes a benchmark for lattice access to multi-dimensional baryon structure.

This paper presents the first non-perturbative lattice determination of the Yang--Mills topological susceptibility slope $\chi^\prime$ in the large-$N$ limit. This quantity represents the $\mathcal{O}(p^2)$ term of the momentum expansion of the topological charge density two-point correlator, and has important theoretical and phenomenological implications for strong interactions. This calculation is based on a novel algorithm that avoids topological freezing at large $N$ on fine lattices, and on a novel method to reliably compute $\chi^\prime$ on the lattice. The results of this study are relevant for the description of the proton spin in deep inelastic scattering experiments via the Shore--Veneziano formula.

Strong magnetic fields produced in relativistic heavy-ion collisions can modify fluctuations of conserved charges and, consequently, their associated chemical potentials. We present first-principles $(2+1)$-flavor lattice-QCD results for isospin-driven splittings of conserved-charge chemical potentials between the isobar systems $^{96}_{44}\mathrm{Ru}+^{96}_{44}\mathrm{Ru}$ and $^{96}_{40}\mathrm{Zr}+^{96}_{40}\mathrm{Zr}$ in the QCD crossover region, both at vanishing and nonzero magnetic fields along the pseudo-critical line $T_{pc}(eB)$. We outline a framework that, under strangeness neutrality and charge-to-baryon ratio $r\equiv n_{\rm Q}/n_{\rm B}$, maps the isospin difference between two nuclei, as encoded in $r_{\rm Zr}$ and $r_{\rm Ru}$, onto splitting ratios $\Delta\mu_{\rm Q}/\Delta\mu_{\rm B}$, $\Delta\mu_{\rm S}/\Delta\mu_{\rm B}$, and $\Delta\mu_{\rm S}/\Delta\mu_{\rm Q}$ as functions of $\mu_{\rm B}(r_{\rm Ru})/\Delta\mu_{\rm B}$. Using continuum-estimated lattice results for the leading-order coefficients $q_1\equiv(\mu_{\rm Q}/\mu_{\rm B})_{\rm LO}$ and $s_1\equiv(\mu_{\rm S}/\mu_{\rm B})_{\rm LO}$, we find that, at vanishing magnetic field, the splitting ratios are of similar magnitude to recent Bayesian extractions from STAR isobar data and yield $\Delta\mu_{\rm Q}<0$ and $\Delta\mu_{\rm S}>0$, with the electric-charge sector dominating. At nonzero magnetic fields, the splitting ratios show only moderate $eB$ dependence. We therefore further examine Ru--Zr differences in the normalized magnetic-field response of chemical-potential ratios, particularly those involving $\mu_{\rm Q}/\mu_{\rm B}$, which display a pronounced enhancement in lattice QCD. We also present hadron resonance gas (HRG) results and experimentally motivated proxy observables with kinematic cuts to facilitate contact with experiment.

Baryon light-cone distribution amplitudes (LCDAs) are inherently multidimensional objects parametrized by two independent longitudinal momentum fractions, making their first-principles determination substantially more challenging than that of meson LCDAs. We present a systematic large-momentum effective theory (LaMET) framework for determining baryon leading-twist LCDAs from lattice QCD. The framework covers the complete path from equal-time three-quark quasi-distribution amplitudes to physical baryon LCDAs. We formulate the leading-twist $V$, $A$, and $T$ quasi-DAs and analyze their spin-flavor and coordinate-space symmetries, including antisymmetric amplitudes with vanishing local limits. We develop a hybrid renormalization prescription on the $(z_1,z_2)$ plane, introduce a newly developed large-$\lambda$ extrapolation strategy based on the asymptotic large-distance behavior of Euclidean correlators, and derive the corresponding one-loop LaMET matching relation in the hybrid renormalization scheme. As a demonstration, we apply the complete analysis pipeline to the $\Lambda$-baryon $A$-structure quasi-DAs using seven $2+1$--flavor lattice ensembles, and use this amplitude to examine the impact of large-distance extrapolation, perturbative matching, and extrapolation to the continuum, physical-pion-mass, and infinite-momentum limits, together with the associated systematic uncertainties. This work provides the formalism, renormalization, extrapolation, and matching infrastructure for first-principles determinations of $x$-dependent baryon LCDAs.

Low-energy $\alpha$-$\alpha$ scattering underlies stellar helium burning and sharply tests nuclear forces in the reaction regime. We present its first calculation using the high-fidelity N3LO chiral NLEFT interaction, incorporated through wave function matching, on a fine lattice, using the adiabatic projection method. On the fine lattice, the two-cluster norm matrix becomes severely ill-conditioned, and its direct inversion is unstable. We address this with Tikhonov regularization, extrapolating the regulator to zero, and confirm the result with an independent truncated singular-value decomposition. The S- and D-wave phase shifts agree with empirical analyses, extending the validation of this interaction from bound states and charge radii to scattering and providing a practical route to ab initio nuclear reactions on fine lattices

We present a formulation for understanding color confinement on the basis of quantum reference frames (QRFs) and relational observables. In the QRF approach to color confinement, colored quantities are not defined as isolated local fields, but rather as relational observables with respect to a color frame or a dressing field. By the Gauss law, local color charge is excluded from the physical bulk algebra, whereas semi-local data such as boundary fluxes and Wilson lines may remain. Color confinement is characterized by the absence of a globally well-defined long-distance color QRF capable of supporting isolated non-singlet relational observables. This formulation preserves the insight of the Kugo-Ojima type picture, while avoiding dependence on a particular covariant gauge, an unbroken global BRST symmetry, and a specific infrared confinement criterion. As concrete examples, we consider (1+1)-dim. Yang-Mills theory, (1+1)-dim. U(1) gauge-Higgs model, and the two-dim. U(1) gauge-Higgs model on $\mathbb{H}^2$ ($AdS_2$) and three-dim. SU(2) gauge-Higgs model on $\mathbb{H}^3$ ($AdS_3$) obtained by dimensional reduction of four-dim. SU(2) Yang-Mills theory restricted to symmetric-instanton sectors. Through explicit calculations in these examples and in controlled sectors, we provide nontrivial consistency checks for the validity of the present formulation. We also discuss prospects for four-dim. Yang-Mills theory and gauge-Higgs theories. QRF-based color confinement provides a relational formulation of why isolated colored asymptotic sectors are absent. At the same time, it clarifies the role played by topological defects and shows that other confinement criteria -- the Wilson-loop area law, the preservation of generalized symmetry, namely center one-form symmetry, and the restoration of residual gauge symmetry -- can be organized as manifestations of a common QRF structure.

We study the scattering of the ground states of doubly charmed baryons ($\Xi_{cc}^{++},\Xi_{cc}^{+},\Omega_{cc}^{+}$) and light-flavor pseudoscalar mesons ($\pi,K,\eta$) up to the next-to-leading order within chiral effective theory. We perform the unitarization of the $S$-wave scattering amplitudes in order to study the excited doubly charmed baryons. The unknown next-to-leading order low energy constants are determined through the fits to recent lattice data in the elastic scattering processes based on the CLQCD ensembles. Following the chiral extrapolation to physical quark masses, we predict resonance, virtual and bound doubly-charmed-baryon states arising from the single- and coupled-channel scattering of $\Xi_{cc}^{++},\Xi_{cc}^{+},\Omega_{cc}^{+}$ with $\pi,K,\eta$. Furthermore, we also calculate the corresponding scattering lengths, effective ranges, phase shifts and inelasticities at physical quark masses, which could shed light on future experimental searches and lattice simulations.

We investigate the impact of including a dynamical charm quark on the properties of light hadrons. Our study compares the calculations performed on 2+1+1 flavor (HISQ fermion) ensembles at four lattice spacings to those on 2+1 flavor (clover fermion) ensembles at six lattice spacings, with both sets of ensembles employing the identical Symanzik gauge action. For the light, strange and charm flavor observables, we employ the same tadpole-improved clover fermion action. From numerical results for light and strange quark masses, pion and kaon decay constants, and $\Omega$ and $\Omega_{ccc}$ baryon masses, we find that the values obtained after continuum, chiral, and infinite-volume extrapolations are consistent within uncertainties. Even though the mixed action setup can introduce additional discretization effects, our calculation shows evidences that those effects can cancel with the discretization error in the unitary setup, resulting in better convergence in the continuum extrapolation.

We investigate the quantum phase diagram of an ``odd'' variant of the two-dimensional Ising Fradkin--Shenker model, characterized by a uniform background of static $e$ and $m$ charges. Using large-scale tensor network and exact diagonalization methods, we determine the topology of the phase diagram, identifying an ``odd'' deconfined phase, confinement- and Higgs-induced valence bond solids (VBS), and a trivial paramagnet. Most notably, we uncover an exotic multicritical point along the self-dual line, where electric and magnetic excitations are related by an enriched $\mathbb{Z}_2$ duality. This transition is marked by the simultaneous onset of confinement, Higgs condensation, translational symmetry breaking, and spontaneous duality symmetry breaking. Within our numerical accuracy, the transition appears continuous, involving the softening of excitation gaps for $e$ and $m$ anyons at finite momentum. At intermediate couplings, we further identify VBS phases with enlarged unit cells, potentially indicating frustration-induced crystalline order beyond commensurate limits.

Understanding the phases of quantum chromodynamics (QCD) at finite temperature and baryon density is crucial for describing matter in heavy-ion collisions and the interior of neutron stars. While the transition from the confined phase to the deconfined phase has been extensively studied at vanishing chemical potential, recent theoretical work suggests the existence of intermediate phases with partial deconfinement, where only a subset of the color degrees of freedom is deconfined. In this paper, we study the generalization of partial deconfinement in QCD in the Veneziano large-$N_{\rm c}$ limit ($N_{\rm f}/N_{\rm c}$ fixed) to finite baryon chemical potential. We find that the partially deconfined phase has finer structure, and hence that there are four phases in QCD at finite temperature and moderately large chemical potential: complete confinement, complete deconfinement, and two kinds of partial deconfinement. A key ingredient is the refined understanding of the meaning of the 'Gross-Witten-Wadia point', i.e., the opening of a gap in the distribution of Polyakov line phases: either string condensation or baryon condensation causes the GWW transition. As a by-product, we observe the emergence of the QCD critical point from the interplay between baryon condensation and partial deconfinement. While our approach is necessarily qualitative due to the fermion sign problem, it provides a unified theoretical framework for understanding the rich phase structure of dense QCD matter and offers new perspectives on the location and nature of the QCD critical point.

The profile of the pion valence quark distribution function (DF) remains controversial. Working from the concepts of QCD effective charges and generalised parton distributions, we show that since the pion elastic electromagnetic form factor is well approximated by a monopole, then, at large light-front momentum fraction, the pion valence quark DF is a convex function described by a large-$x$ power law that is practically consistent with expectations based on quantum chromodynamics.