2026-03-23 | | Total: 13
Construction of creation operators which can properly sample the underlying energy eigenstates remains a fundamental first step in lattice QCD spectroscopy calculations, particularly when the spectrum includes states with different composition such as mesons, glueballs, multi-particle states, etc. We tackle this issue in the study of the scalar glueball and charmonium mixing, where we use improved operators for both types of states to resolve the low-lying spectrum and identify the dominant composition of each state in a mass regime where the glueball is stable. We include derivative-based meson operators combined with distillation profiles, as well as glueball operators built from the chromo-magnetic field and its derivatives which retain angular momentum information from their continuum counterparts. We comment on the advantages of these operators, particularly on the construction and implementation of the glueball ones, thanks to which we identify the lightest iso-scalar state as glueball-dominated $0^{++}$.
When the number of massless fermions exceeds a critical value $N_f^*$, QCD enters the conformal window and becomes chirally symmetric already in the vacuum. Determining $N_f^*$ from lattice simulations is challenging, since calculations are performed at finite lattice spacing, quark mass, and temporal lattice size, where both a thermal transition and an unphysical bulk transition obscure the conformal behaviour. In this work, we present results on the chiral phase boundaries in the bare lattice parameter space $(N_τ,\;β,\;am,\;N_f)$ of unimproved staggered fermions. Our analysis indicates that the chiral transition in continuum QCD is of second order for all $N_f$ up to the onset of the conformal window. By systematically studying the thermal chiral transition and its interplay with the bulk transition, we obtain a coherent picture of the lattice phase structure and suggest how the onset of the conformal window can be identified from simulations performed away from the chiral and continuum limits.
We report our ongoing lattice QCD study of radiative leptonic decays of the charged pseudoscalar mesons $D$, $D_s$, $B$, and $B_c \to \ell ν_\ell γ$. We carry out our analysis on a single JLQCD ensemble with lattice spacing $a=0.044~\text{fm}$. This work is a step towards a complete QCD+QED lattice calculation of these modes, aimed at reducing theoretical uncertainties in the extraction of $|V_{cd}|$ and $|V_{cs}|$ and providing first-principles estimates of the corresponding form factors in the $B$ sector.
We report on our lattice QCD study of coupled $DD_s^* - D^*D_s$ scattering in the $J^P=1^+$ channel and elastic $DD_s$ scattering in the $J^P=0^+$ channel, aimed at investigating the possible existence of $cc\bar{u}\bar{s}$ tetraquarks near threshold. The calculation uses CLS ensembles with $m_π\approx 280$ MeV, lattice spacing $a \approx 0.09$ fm, and spatial extents $L/a=24, 32$. Finite-volume spectra are obtained from a variational analysis of two-point correlation matrices constructed from two-meson operator bases using distillation. The $l=0$ partial-wave scattering amplitudes are determined from the lattice spectra in multiple moving frames using Lüscher's formalism as well as a finite-volume implementation of the Lippmann-Schwinger equation. In both channels we observe small but nonzero shifts relative to the noninteracting spectrum, indicating weak meson interactions. The extracted physically plausible $S$-wave amplitudes show no pole structures near threshold.
We report new results obtained in our lattice studies of the $Sp(4)$ gauge theory coupled to two fundamental Dirac fermions. This theory provides a candidate for the dynamical origin of dark matter models within the strongly interacting massive particle paradigm. We employ Lüscher's formalism to analyse finite-volume energy levels and study the scattering amplitude of two pseudoscalar states in the spin-1 channel. We present our preliminary findings for a set of ensembles generated within a broad range of (Wilson) fermion masses.
We develop a one-class, deep-learning framework to detect the phase transition and recover critical behavior of the 3D Ising model. A 3D convolutional neural network autoencoder (CAE) is trained on ground-state configurations only, without prior knowledge of the critical temperature, the Hamiltonian, or the order parameter. After training, the model is applied to Monte Carlo configurations across a wide temperature range and different lattice sizes. The mean-square reconstruction error is shown to be sensitive to the transition. Finite-size scaling of the peak location for the reconstruction error susceptibility yields the critical temperature $T_c=4.5128(58)$ and the correlation-length critical exponent $ν=0.63(27)$, consistent with results from the literature. Our results show that a one-class CAE, trained on zero-temperature configurations only, can recover nontrivial critical behavior of the 3D Ising model.
Recent progress in experimental facilities, together with larger data samples and more refined analysis strategies has enabled the observation of many exotic hadronic states, adding new members to the hadron spectrum. Each newly reported signal encourages further experimental searches and simultaneously motivates theoretical studies aimed at uncovering additional nonconventional states. Motivated by this perspective and by the increasing interest in systems containing multiple heavy quarks, we present a spectroscopic study of fully heavy pentaquark candidates with spin-parity quantum numbers $J^{P}=\frac{1}{2}^{-}$ and quark contents $QQQ'Q\bar{Q'}$, $QQQ'Q'\bar{Q}$, and $Q'Q'QQ\bar{Q}$, where $Q(Q')$ represents either $c(b)$ or $b(c)$ quarks. We employ the QCD sum rule approach with three different types of interpolating currents to obtain the corresponding masses and current coupling constants of the considered states. The following masses for the states containing three $c$ and two $b$ quarks are predicted: $m_{(3c2b)}=14479.30\pm75.06~\mathrm{MeV}$ using the current $J_1$, $\tilde{m}_{(3c2b)}=14276.80\pm76.29~\mathrm{MeV}$ using $J_2$, and $\bar{m}_{(3c2b)}=14276.80\pm76.29~\mathrm{MeV}$ using $J_3$. The corresponding predictions for the states containing three $b$ and two $c$ quarks are as $m_{(3b2c)}=17458.90\pm130.11~\mathrm{MeV}$ with $J_1$, $\tilde{m}_{(3b2c)}=17202.70\pm132.37~\mathrm{MeV}$ with $J_2$, and $\bar{m}_{(3b2c)}=17250.80\pm131.98~\mathrm{MeV}$ with $J_3$, respectively. In addition, we provide the corresponding current coupling constants, which can serve as useful inputs for analyses of decay properties and interaction mechanisms of these fully heavy pentaquark candidates.
We derive the factorization theorem for the quasi-transverse-momentum-dependent (quasi-TMD) correlator, including kinematic power corrections to all orders. The resulting expression involves only twist-two TMD distributions and is frame invariant. As in the leading-power approximation, the reduced soft factor factorizes multiplicatively; however, in contrast, the TMD evolution factor is not multiplicative but enters as a convolution with the nonperturbative TMD distribution. We present a numerical estimate of the difference between the derived expression and the leading-power approximation, finding it to be of the order of tens of percent for current lattice simulations. We also demonstrate that the inclusion of these corrections improves the agreement between phenomenological extractions of the Collins-Soper kernel and lattice results.
We introduce homotopy lattice gauge fields (HLGFs), a version of gauge fields over a discretized base, based on a notion of higher parallel transport that enriches the usual parallel transport along paths on a lattice to also consider higher dimensional paths. Higher dimensional data keeps information about the parallel transport along homotopies of curves. With this data, a HLGF on a base space of dimension two or three determines a principal bundle over the base manifold. This data is also responsible for our formulas for the topological charge on two-dimensional bases. Our framework is an application of a nonabelian algebraic topology framework developed to solve the local to global problem in higher dimensional homotopy. No previous knowledge of higher category theory is assumed. The second part will be devoted to the space of fields as an arena for doing Quantum Field Theory, and to give the first examples of how our framework refines standard lattice gauge theory.
This paper investigates the crucial role of parton distribution functions (PDFs) in high-energy physics, particularly their impact on the extraction of generalized parton distributions (GPDs) at zero skewness. To this aim, we perform six global analyses of GPDs using different modern PDF sets (\texttt{NNPDF40}, \texttt{CT18}, and \texttt{MSHT20}) at three specific factorization scales ($ μ= 2 $, $ 1.3 $, and $ 1 $ GeV) and different perturbative orders. A wide range of elastic electron scattering data is included in the analysis to constrain GPDs in a broad interval in the momentum transfer squared $ t $. The analyses reveal that the best overall description of experimental data is achieved using \texttt{NNPDF40} PDFs at the next-to-leading order (NLO), with moderate sensitivity to the choice of PDF input, especially in the region of low $ t $. The dependence on the perturbative order is relatively mild, indicating stability in the extraction procedure. We also show that different GPD sets become more consistent at larger $ |t| $ values, with down-quark GPDs experiencing greater suppression than up-quark GPDs. The six extracted GPD sets, available at different scales and perturbative orders, provide valuable resources for future theoretical and phenomenological studies, offering flexibility for researchers in exploring the proton's internal structure.
We demonstrate the renormalisability of quantum field theories in four dimensions with elementary self-interacting Dirac fermions and to leading order in the limit of many fermion flavours $N_{\rm f}$. Starting from the underlying divergence structure and using Gross-Neveu-type interactions as a template, we explain why extended four-fermion theories including higher-derivative interactions are well-defined, renormalisable, and predictive with only a few free parameters. We also provide the exact large-$N_{\rm f}$ leading beta functions of couplings and discuss quantum scaling dimensions, universality, $1/N_{\rm f}$ corrections, and extensions to other types of four fermion interactions. Implications for effective theory and model building are indicated.
The infrared structure of gauge theories with chiral fermions remains largely unexplored. In this work we investigate the Bars-Yankielowicz class using the functional renormalisation group, building on recent developments in gauge-fermion systems that provide clear criteria for confinement and dynamical symmetry breaking. We show that two distinct phases arise: one exhibiting both confinement and symmetry breaking at small numbers of colours, and another characterised by confinement without symmetry breaking in the large-colour limit. The latter realises a novel regime, opening the possibility of exotic spectra and phenomena that can now be studied within a systematic framework.
Advances in LLMs have produced agents with knowledge and operational capabilities comparable to human scientists, suggesting potential to assist, accelerate, and automate research. However, existing studies mainly evaluate such systems on well-defined benchmarks or general tasks like literature retrieval, limiting their end-to-end problem-solving ability in open scientific scenarios. This is particularly true in physics, which is abstract, mathematically intensive, and requires integrating analytical reasoning with code-based computation. To address this, we propose PhysMaster, an LLM-based agent functioning as an autonomous theoretical and computational physicist. PhysMaster couples absract reasoning with numerical computation and leverages LANDAU, the Layered Academic Data Universe, which preserves retrieved literature, curated prior knowledge, and validated methodological traces, enhancing decision reliability and stability. It also employs an adaptive exploration strategy balancing efficiency and open-ended exploration, enabling robust performance in ultra-long-horizon tasks. We evaluate PhysMaster on problems from high-energy theory, condensed matter theory to astrophysics, including: (i) acceleration, compressing labor-intensive research from months to hours; (ii) automation, autonomously executing hypothesis-driven loops ; and (iii) autonomous discovery, independently exploring open problems.