2026-05-08 | | Total: 35
This work introduces a novel, nonparametric pixel-based framework for the Bayesian inference and imaging of transverse momentum dependent (TMD) parton distributions. The methodology is built upon a fully differentiable framework that integrates TMD evolution with the Collins-Soper-Sterman formalism, enabling the simultaneous extraction of partonic distributions and the nonperturbative evolution kernel. To achieve efficient and exact sampling of the high-dimensional posterior, we leverage generative AI through a hybrid normalizing flow-driven Metropolis-Hastings approach. The framework is validated through multi-scale closure tests of increasing complexity, ranging from basic functional models to convoluted structure functions. Using singular value decomposition (SVD), we rigorously characterize the uncertainty of the reconstructed distributions and reveal the existence of null TMDs, which are functional components in the null space of the kernel that remain unconstrained by observables. The new framework provides the first integration of pixel-based discretization, generative AI, and SVD within a Bayesian context to solve the TMD inverse problem. This synergy between machine learning and multi-scale data removes inherent degeneracies and enables unbiased 3D partonic imaging.
In this paper we revisit, from a different perspective, a long-standing question: ``Is the systematic deficit observed in all branching ratios mediated by a $b\to s μμ$ transition due to New Physics, or to a hypothetical constant unknown universal hadronic contribution that mimics New Physics?'' The key observation that allows us to distinguish between these two possibilities is that non-perturbative contributions associated with $c\bar c$ loops affect inclusive $B\to X_s \ell\ell$ and exclusive $B\to K^{(*)}\ell\ell$ modes differently. In inclusive decays, factorizable contributions are exactly determined from data on $e^+e^-\to \mathrm{hadrons}$, while non-factorizable corrections are described by resolved-photon contributions at low $q^2$ and by local power corrections at high $q^2$. In exclusive decays, by contrast, hypothetical charming-penguin effects, beyond those already included in current uncertainty estimates, could appear, in a worst-case scenario, as a constant, universal contribution that it seems, in principle, indistinguishable from genuine New Physics. We identify two observables, constructed from ratios of exclusive to inclusive modes, that can discriminate between a New Physics contribution and a constant hadronic contribution. Moreover, these ratios can be measured directly by LHCb, as they do not require any normalisation to $J/ψK^{(*)}$ branching fractions from B factories. A preliminary evaluation of these observables with present data shows some preference for the New Physics interpretation. In a complementary test, a comparison between inclusive measurements and the corresponding sum of exclusive modes at high $q^2$ similarly disfavours an explanation based on a constant hadronic contribution. Finally, we provide projections for the new observables based on expected LHCb and Belle II measurements in the near future.
We develop a quantum spin-density-matrix framework for heavy-quark spin dynamics in relativistic heavy-ion collisions. Starting from an initial polarization induced along the magnetic-field direction, we derive the evolution equation for spin polarization within this framework and obtain analytic solutions. The evolved polarization is connected to open heavy-flavor observables via a fragmentation-based hadronization prescription. For vector mesons, the spin-alignment parameter $ρ_{00}$ is constructed by coupling the heavy-quark spin to that of the light antiquark produced during fragmentation. We confront our results with recent ALICE measurements of prompt $D^{*+}$ spin alignment in Pb--Pb collisions at $\sqrt{s_{\rm NN}}=5.02~{\rm TeV}$ and extract an effective depolarization strength that determines the spin-relaxation time scale. Using this fitted parameter, we provide benchmark estimates for $Λ_c^+$ and $\barΛ_c^-$ polarization, up to an overall spin-transfer normalization. We further estimate the recently proposed elliptic polarization harmonic arising from path-length-dependent depolarization in an anisotropic fireball.
The leading Generalized Transverse Momentum Dependent parton distributions (GTMDs) are studied in the bag model. The model description is shown to be theoretically consistent. The orbital angular momentum is studied in terms of the GTMD $F_{1,4}^q$ and Ji sum rule. Analytical proofs of the associated sum rules are given. A deeper relationship between orbital angular momentum and the pretzelosity TMD is established in this model.
The interpretation of the SN1987A neutrino data continues to be hindered by significant absolute timing uncertainties and a well-known tension in the angular distributions. We perform a quantitative, high-precision alignment of the Kamiokande-II and Baksan data with the IMB clock, using for the first time the relative time offset (RTO) and reducing the temporal uncertainty by two orders of magnitude to the sub-second level. The chi-square analysis shows that Baksan's absolute timestamps require an advancement of 30.4 s, while those of Kamiokande-II require a delay of about 6.4 s. The knowledge of the unified timeline provides a necessary and rigorous basis for testing whether the first Kamiokande-II event arose from the neutronization burst - a possibility motivated by its angular properties. Our analysis favors an accretion phase electron anti-neutrino origin over a neutronization-burst origin, with a likelihood ratio of 3-6, depending on the specific MSW oscillation scenario. This result corroborates the standard interpretation - that only inverse beta decay events were detected. Our framework yields the most stringent constraints to date on the SN1987A chronology and establishes a precision benchmark for future Galactic supernova observations.
In this work, we perform a dedicated analysis of fermion mass fits in the minimal $SO(10)$ grand unified theory (GUT), going beyond the tree-level approximation by incorporating one-loop finite corrections to the neutrino mass matrix. We show that parameter regions that successfully reproduce all fermion masses and mixings at tree level can lead to significant deviations in neutrino masses and leptonic mixing parameters once the radiative corrections are included. These results expose a limitation of conventional tree-level fitting procedures and highlight the sensitivity of neutrino observables to loop effects. Since in the minimal $SO(10)$ GUT the same set of Yukawa parameters simultaneously governs quark masses, charged lepton masses, and neutrino properties, these radiative corrections propagate across all fermion sectors, reshaping the viable parameter space in a highly non-trivial and correlated manner. We find that the largest corrections to the masses and mixing angles are of order $\mathcal{O}(30\%)$-$\mathcal{O}(40\%)$, therefore, cannot be neglected. In light of the current precision of neutrino oscillation measurements, and the expected improvements from ongoing and future experiments, we demonstrate that the inclusion of one-loop effects is essential for a consistent and reliable exploration of the parameter space, with important implications for the predictivity of $SO(10)$ GUTs.
We study the origin of the non-monotonic magnetic-field dependence of the lowest charged pion excitation observed in lattice QCD. In a magnetic field, the charged pion mixes with the longitudinally polarized charged rho meson, which shares the same quantum numbers. Within the SU(2)$_f$ Nambu--Jona-Lasinio model supplemented by a gauge invariant tree-level $π-ρ$ mixing operator constrained by the experimental $ρ^\pm\rightarrowπ^\pmγ$ decay width, we compare four mass-extraction schemes: rest-mass reconstruction, local bosonization, direct determinant solving with Landau projection, and near-pole expansion. The rest-mass scheme cannot reproduce the lattice-type turnover, while in the local derivative-expansion scheme the turnover presence but is weak which occurs at large magnetic field. By contrast, the direct determinant and near-pole schemes both retain a robust non-monotonic lowest mode. The former is most faithful to the Landau-level kinematics of the charged excitation, while the latter most clearly shows that residue suppression enhances the effective mixing after canonical normalization. Our results indicate that the lattice behavior is a genuine quasiparticle mixing effect, but one whose robustness depends crucially on how the charged-meson pole structure is extracted in a magnetic field.
In the present work we explore for the first time the general structure and properties of the nonperturbative quark-gluon vertex in the complex plane. Specifically, we focus on the transversely-projected quark-gluon vertex that emerges from a recently developed symmetry-preserving approach for the study of meson properties beyond the rainbow-ladder approximation. The analysis focuses on the so-called "soft-gluon" limit, which reduces the momentum-dependence of the corresponding vertex form factors to a single momentum variable. The complexification of this variable inside the defining integrals furnishes unambiguously all eight vertex form factors within a concrete domain of the complex variable, delimited by a characteristic parabola. The extent of this reliable domain is determined by the appearance of the first singularity in the integrands of the vertex integrals, where the standard Wick rotation must be duly supplemented by additional crucial contributions. This primary analytic region may be extended considerably by resorting to standard extrapolation methods, which remain valid up until the appearance of complex structures associated with the onset of physical processes. The generalization of the method to arbitrary gluon momenta, and its relevance for the determination of the quark propagator in the complex plane, are briefly discussed.
We investigate freeze-in production of light dark matter through the quark flavor-changing operator $(\bar{s}γ_μd)(\barχγ^μχ)$ in a low-reheating cosmology. For reheating temperatures below the QCD crossover, kaon decays and scatterings generate the dark matter abundance through $K\toπχ\barχ$ and $Kπ\toχ\barχ$. The same interaction induces the rare kaon decays $K^+\toπ^+χ\barχ$ and $K_L\toπ^0χ\barχ$. This links the freeze-in relic abundance to searches at NA62, KOTO, and KOTO II. We find that lower reheating temperatures require larger couplings to compensate for the Boltzmann-suppressed kaon abundance, making kaon-driven freeze-in dark matter testable at rare kaon decay experiments.
We have studied the interaction of $π^0 (η) f_1(1285)$ assuming the $f_1(1285)$ to be a molecular state of $K^* \bar K - \bar K^* K$. We use a framework in which a $π^0 (η) f_1(1285)$ optical potential is obtained, which is later used as the kernel of the Lippmann-Schwinger equation, following the standard method for the interaction of particles with nuclei. The optical potential is obtained using the fixed center approximation to the Faddeev equations, where a cluster, here the $f_1(1285)$, remains unchanged during the interaction, appropriate to the situation that one has here. We have obtained the scattering matrix for this system, the scattering length and effective range, plus the correlation functions. The framework used has been previously tested in the study of the $p f_1(1285)$ interaction and has been shown to give results in agreement with the recent experimental measurement of the $p f_1(1285)$ correlation function. On the other hand, from this interaction we do not obtain clear signals for the $π_1(1400)$ or $π_1(1600)$, nor for the $η_1(1855)$ resonances, which in other approaches have been claimed to arise from the same dynamics. We, however, obtain a structure in the $π^0 f_1(1285)$ amplitude around $1500-1600$ MeV and a strong cusp at the $ηf_1(1285)$ threshold of $1833$ MeV.
Under certain conditions, the production of a Higgs boson in association with a top-anti-top pair at hadron colliders can be described via a factorisation theorem using perturbative fragmentation functions. The latter describe the nearly collinear emission of a Higgs boson from a top-quark and reproduce the leading mass dependence of the exact next-to-leading-order (NLO) calculation. Although the NLO fragmentation functions have been calculated a few years ago, it has not been possible up to now to demonstrate the applicability of the approximation in a realistic setup. At NLO, we analyse two different ways of treating the top-quark mass, called the zero-mass-top-quark (ZMTQ) and the hybrid prescription. We show that the method yields reliable results at LHC center-of-mass (cms) energies in the hybrid prescription. In the ZMTQ prescription, the results at LHC cms energies are only reliable in the quark-anti-quark channel, but become viable for the full $pp \rightarrow t\bar{t}H$ process at a 100 TeV hadron collider. In addition, we discuss some subtleties and complications arising when extending the formalism to next-to-next-to-leading-order (NNLO) and beyond.
We discuss the properties of the extended three-flavor quark-meson diquark (EQMD) model as a renormalizable low-energy effective model for QCD. The effective degrees of freedom are quarks, scalar- and pseudoscalar mesons, diquarks, vector- and axial-vector mesons. We calculate the equation of state (EoS) in the mean-field approximation at $T=0$ imposing charge neutrality for electric and color charges. We match the EoS with a low-density nuclear equation of state. We discuss how the choice of parameters in the model affects the EoS and thereby the mass-radius for hybrid stars. We show that it is possible to construct hybrid stars whose masses and radii are in agreement with recent astrophysical observations and perturbative QCD (pQCD). The addition of vector and axial vector mesons to the quark-meson diquark is essential, since it makes the EoS sufficiently stiff for intermediate densities. Our results suggest that stars with a mass larger than $M\sim2M_{\odot}$ have a quark core with a central density $n_B\geq 3.9n_{\rm sat}$, where $n_{\rm sat}\approx0.165$fm$^{-3}$ is the saturation density. The speed of sound has a double-peak structure and relaxes to the conformal limit from above for large baryon chemical potentials $μ_B$. This structure is caused by the decrease in the mass of the $s$ quark as $μ_B$ increases.
We systematically study the mass spectra of $S$-wave singly-heavy tetraquark states $Qq\bar{q}\bar{q}$ ($Q=c,b$; $q=u,d,s$) in a mass splitting model. We adopt the assumption that the $X(4140)$ is the lowest $J^{PC}=1^{++}$ $cs\bar{c}\bar{s}$ tetraquark and use this state as a reference to determine the mass splittings. According to the obtained results, we also estimate the rearrangement decay widths of the tetraquarks within a simple scheme. We find that the recently observed states $T^a_{c\bar{s}0}(2900)^{++/0}$ and $T_{cs0}^*(2870)^0$ by the LHCb Collaboration can be consistently interpreted as the second highest $I(J^P)=1(0^+)$ $cn\bar{s}\bar{n}$ ($n=u,d$) and the higher $I(J^P)=0(0^+)$ $cs\bar{n}\bar{n}$ tetraquark states, respectively. We predict several narrow tetraquark candidates: the lowest $cn\bar{s}\bar{n}$ and $cs\bar{n}\bar{n}$ with $I(J^P)=0(0^+)$ and $0(1^+)$, and their bottom counterparts. The obtained information from mass spectrum and rearrangement decay properties will help search for the new singly-heavy tetraquark states.
This study constrains the range of in-medium mass modification through a comparison of theoretical calculations with experimental transverse momentum spectra and the yield ratio {\bar{p}}/p of protons and antiprotons. Based on the constrained range and a Gaussian source model with radial ow, the theoretical predictions for the fermion back-to-back correlation (fBBC) of p{\bar{p}} pairs at RHIC energies are presented. The results reveal a strong sensitivity of the fBBC signal to the assumed source time distribution: a Lorentzian form generates a pronounced high-momentum signal, whereas an α-stable Lévy form leads to a marked low-momentum signal. Moreover, the in-medium mass modification is shown to enhance the yield ratio {\bar{p}}/p. Therefore, events characterized by a larger {\bar{p}}/p ratio are predicted to have a significantly higher probability of exhibiting a detectable fBBC signal. This study may propose a promising new direction for the experimental observation of this phenomenon
We extract a target-side Cu transverse-momentum profile from fixed-target $p$--Cu Drell--Yan data by holding a momentum-space proton reference fixed and training only an asymmetric Cu kernel in the small-$q_T$ region. In the supported window, $0.15 \le x_{Cu} \le 0.46$ and $7.5 \le Q_M \le 15.75$ GeV, the nuclear modification is not a universal width increase. It appears as $Q$-dependent redistribution: an $\mathcal{O}(1~{\rm GeV})$ shoulder and compensating probability flow between shoulder and resolved-tail regions, beyond one-parameter broadening.
We analytically calculate the all-loop bare perturbative part of the four-quark Bethe-Salpeter kernel using modern scattering amplitude methods. We work to subleading order in the large number of quark flavors approximation of massless Quantum Chromodynamics, which simultaneously makes an all-loop calculation feasible, is systematically improvable, and preserves asymptotic freedom. It also allows for avoiding the ambiguity of choosing a truncation scheme in Dyson-Schwinger equations. We exploit state-of-the-art methods in Integration-By-Parts reduction of Lorentz scalar Feynman integrals into a minimal Master Integral basis, and direct integration into Generalized Polylogarithms. As a byproduct of our calculation, we also provide the result for the gluon and quark propagators. We discuss a path towards nonperturbative formulation and potential future phenomenological applications.
We present MAPFF1.0_Lambda, the first global analysis at next-to-next-to-leading order in perturbative QCD of the collinear unpolarised fragmentation functions of Lambda hyperons. The fit is based on data from single-inclusive electron-positron annihilation, and from both neutral-current and -- for the first time -- charged-current semi-inclusive deep-inelastic scattering. We have adopted a statistical framework based on Monte Carlo sampling and parametrised fragmentation functions in terms of a neural network. The fragmentation function set comprises a total of seven independent parton flavours, allowing for the first independent determination of valence-quark distributions. Our analysis offers new insights into the hadronisation mechanism of strange baryons and establishes a baseline for future phenomenological and experimental investigations.
The scotogenic model defines a framework for radiative neutrino masses and provides a viable dark matter candidate. Since the scotogenic dark matter is leptophilic, indirect searches appear as an especially interesting possibility. Here we propose a simple variation of the model with a very distinct phenomenology. The scotogenic fermion singlets are naturally grouped into pseudo-Dirac pairs of mass of 0.1-1 TeV. We show that the lightest one constitutes a dark matter candidate that near threshold annihilates with a 90% branching ratio into neutrino pairs. The model gives the observed relic abundance consistently with the bounds from direct searches and with all neutrino and charged lepton data. We also show that, for a sub-MeV dark matter particle, the model suggests a scenario that could address the lithium problem.
The annihilation of self-interacting dark matter with long-range interactions can be significantly enhanced at low velocities through the Sommerfeld effect. At special points in parameter space, where near-zero-energy resonances exist in the spectrum of the theory, the standard calculation of this enhancement appears to violate unitarity. Recently, several approaches have been proposed to regulate this behavior, some introducing explicit ultraviolet (UV) scales and others not, raising the question of whether these prescriptions are consistent. In this paper, we compare these approaches and show that even in nominally cutoff-dependent methods, the regulated Sommerfeld-enhanced cross sections are independent of the UV regulator to a good approximation, and that when the unitarity-preserving corrections are large, the regulation schemes coincide to leading order. We use these insights to write down a regulator-independent prescription for unitarization applicable to multi-state systems, where the modified enhancement can be written solely in terms of the standard enhancement factor, the hard annihilation amplitude, and the $S$-matrix for scattering in the long-range potential.
We explore the coincidence of scales between cosmic inflation and right-handed neutrinos in seesaw models. We show that inflation models, which will be tested by next-generation CMB experiments, can produce right-handed neutrinos in sufficient abundance to explain the observed baryon asymmetry of the universe. The model can be tested by gravitational wave signatures from cosmic inflation and particle production.
We introduce a new strategy for compositional neural surrogates for radiation-matter interactions, a key task spanning domains from particle physics through nuclear and space engineering to medical physics. Exploiting the locality and the Markov nature of particle interactions, we create a \emph{next-particle prediction} kernel using hybrid discrete-continuous transformer models based on Riemannian Flow Matching on product manifolds. The model generates variable-sized typed sets of particles and radiation side effects that are the result of the interaction of an incident particle with a material volume. The resulting kernel can be composed to simulate unseen large-scale material distributions in a zero-shot manner. Unlike mechanistic simulators, our model is designed to be differentiable, provides tractable likelihoods for future downstream applications. A significant computational speed-up on GPU compared to CPU-bound mechanistic simulation is observed for single-kernel execution. We evaluate the model at the kernel level and demonstrate predictive stability over multi-round autoregressive rollouts. We additionally release a novel 20M-event radiation-matter interaction dataset for further research.
We present a new set of lattice QCD ensembles with four flavors of smeared Möbius 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_{π^{\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.
The Hagedorn regime of string theory is usually understood as an equilibrium limiting phenomenon: the exponential growth of the density of states makes the canonical partition function singular at the Hagedorn temperature, while in the microcanonical description additional energy is absorbed predominantly by highly excited long-string configurations. In this work we revisit this regime from a nonequilibrium perspective using Steepest-Entropy-Ascent Quantum Thermodynamics (SEAQT), where thermodynamic evolution is formulated directly on the state manifold and does not require a globally well-defined canonical ensemble. The inverse temperature is treated as an instantaneous, state-dependent quantity, and we derive its exact scalar evolution equation. In the commuting limit, this dynamics is controlled by higher-order fluctuation moments, showing that the Hagedorn regime may act as a dynamical bottleneck for the response of the effective intensive variable. We then extend the construction to an open-system setting through a system--reservoir splitting of the SEAQT metric and show that reservoir coupling can drive the subsystem toward effective Hagedorn slowing-down. A diagonal Hagedorn evaluation further shows that the strength of this bottleneck depends not only on the exponential density of states, but also on its algebraic prefactor. These results provide a nonequilibrium interpretation of Hagedorn behavior and suggest a connection between long-string dominance, thermodynamic slowing-down, and the breakdown of effective descriptions in quantum gravity.
The density dependence of the nuclear symmetry energy remains one of the key uncertainties in contemporary nuclear physics, with significant implications for the structure of exotic nuclei, the dynamics of heavy-ion collisions, and the properties of astrophysical objects such as neutron stars and core-collapse supernovae. However, extracting robust constraints requires observables that are minimally affected by final-state interactions and are reliably predicted by transport models. This review synthesizes recent theoretical and experimental advancements in constraining the symmetry energy by leveraging isospin diffusion in heavy-ion reactions within the Fermi energy domain. Recent results from the INDRA-FAZIA collaboration, including isospin transport ratio data, and Boltzmann-Uehling-Uhlenbeck (BUU) transport model calculations are highlighted. Confidence regions for the symmetry energy are extracted from isospin transport ratios and isospin diffusion currents by utilizing state-of-the-art nuclear functionals, including both ab initio and phenomenological approaches, with a particular focus on the density regions probed by these experiments. The resulting constraints will aid future Bayesian studies of the nuclear equation of state and contribute to a more unified understanding of dense matter in both terrestrial experiments and astrophysical environments.