2026-03-16 | | Total: 10
(See paper for full abstract) This study reports on the benchmarking of a protocol for the adjustment of the parameters of a local leading-order (LO) T=1 (like-particle) pairing EDF that consists in adjusting the density-dependence of the 1S0 pairing gap at the chemical potential in infinite nuclear matter (INM). When using a suitably chosen reference calculation, this protocol leads to consistent results for the odd-even staggering of masses of spherical and heavy deformed nuclei and also for the rotational moments of inertia calculated in a time-reversal-breaking cranked HFB approach. The implementation of the HFB solver for infinite matter at arbitrary isospin asymmetry used for this study is sketched in appendices. Additional points that are discussed concern (i) the illustration that the gaps at the chemical potential are not necessarily sufficient to completely characterise the pairing interaction in infinite matter, and that adjusting LO pairing EDF to reproduce gaps obtained from finite-range pairing interactions in HFB or from more microscopic calculations can lead to unrealistic predictions for finite nuclei; (ii) the finding that some regions of the parameter space for the density dependence of the T=1 LO pairing EDF lead to a spurious transition to a Bose-Einstein condensate of di-nucleons in spite of producing realistic pairing correlations for well-bound nuclei; (iii) the correlation between effective mass and the parameters of the LO pairing EDF that control the form factor of the density dependence when reproducing pairing gaps in infinite matter, (iv) the significant impact on the odd-even staggering of masses made by either including or not the spin-gradient terms in the particle-hole part of the Skyrme EDF; (v) the sizeable impact of keeping or not the contribution from the density-dependent LO pairing EDF to the mean fields on the odd-even staggering of radii.
Reaction-rate calculations relevant to r-process nucleosynthesis depend sensitively on the nuclear $γ$-strength function ($γ$SF). Here we investigate the impact of low-lying pygmy dipole strength (PDS) in $(n,γ)$ and $(γ,n)$ reactions using $γ$SF based on relativistic nuclear energy density functional theory and propagate these strengths into Hauser--Feshbach statistical model calculations of the reaction rates. We show that considerable reaction-rate enhancements at temperatures relevant for r-process nucleosynthesis are governed by the alignment of the pygmy dipole strength energy with the neutron separation threshold $S_n$ rather than by the total low-energy strength. Consequently, nuclei such as $^{68}$Ni and $^{132}$Sn, where the PDS energy-$S_n$ alignment occurs, exhibit the strongest effects on reaction-rate enhancements. These results demonstrate that modeling reliable reaction rates in r-process nucleosynthesis necessitates accurate microscopic descriptions of low-energy dipole strength, in close synergy with experimental investigations in the vicinity of neutron threshold.
The most accurate approach to determine the compressibility of nuclear matter remains the one based on microscopic Energy Density Functionals (EDFs). Recent analyses yield a value for nuclear incompressibility modulus $K_\sat=240\pm 20$~MeV, defined in nuclear matter as the second derivative of the energy per particle at saturation density. However, we demonstrate that the compressibility modulus can be reduced to values shifted by four times the suggested uncertainty, i.e., $K_\sat\approx 160$~MeV, by providing examples based on models where the second derivative ($K_\sat$) and third derivative ($Q_\sat$) of the energy per particle at saturation density can be independently varied, while the experimental binding energies, charge radii, and ISGMR data in $^{120}$Sn and $^{208}$Pb are enforced. The present work suggests a new methodology to access the compressibility of nuclear matter from nuclear experiments, still based on microscopic models, but using EDFs containing more flexibility than the ones employed up to now. Consequences of our results for nuclear matter at supra-saturation density are also discussed by exploring the quarkyonic cross-over. We predict that, for our models with low values for $K_\sat$, the quark onset density has to be low for neutron stars to exist.
Measurements of the event-by-event correlation between elliptic flow ($v_2$) and the mean transverse momentum ($[p_{\rm T}]$) using the modified Pearson correlation coefficient $ρ(v_2^2,[p_{\rm T}])$ are reported in pp collisions at $\sqrt{s} = 13$ TeV, and in p-Pb and Pb-Pb collisions at a center-of-mass energy per nucleon pair $\sqrt{s_{\rm NN}} = 5.02$ TeV. This analysis is based on the full LHC Run 2 dataset recorded by ALICE and is performed for the first time in small collision systems with the ALICE detector. In Pb-Pb collisions, the $ρ(v_2^2,[p_{\rm T}])$ measurement shows a non-monotonic dependence on charged particle multiplicity ($N_{\rm ch}$); it first decreases and then increases with an increase in multiplicity. The decreasing trends of $ρ(v_2^2,[p_{\rm T}])$ with increasing multiplicity are also observed in p-Pb and pp collisions for the presented multiplicity range. All three systems show consistent values of $ρ(v_2^2,[p_{\rm T}])$ for $N_{\rm ch} \lesssim 80$. These measurements are also compared with theoretical model calculations, including PYTHIA, where no collectivity is generated, as well as AMPT and IP-Glasma + MUSIC + UrQMD, which produce collective effects in small systems. These comparisons offer unique insights into the origin of collectivity in small systems. They improve the understanding of the initial geometry, size, and their correlations. The comparison also allows an investigation of the role of initial momentum correlations predicted by the Color Glass Condensate framework. The new measurements could not be explained by current state-of-the-art models, offering insights into the initial stage of collisions in small systems and also imposing strong constraints on the existing theoretical models. This will significantly advance our understanding of the collective phenomena observed in small systems at the LHC.
This paper presents the measurement of $π^{\pm}$, K$^{\pm}$, and p($\bar{\text{p}}$) production in high-multiplicity proton-proton collisions at $\sqrt{s} = 13$ TeV at midrapidity ($|y| < 0.5$) using the ALICE detector at the LHC. The transverse-momentum ($p_{\rm T}$) spectra of these particles are reported for three high-multiplicity classes. The results show a mass- and multiplicity-dependent hardening of the $p_{\rm T}$ spectra and an enhancement of the p/$π$ ratio at intermediate $p_{\rm T}$. These features are similar to those observed in heavy-ion collisions, where quark-gluon plasma formation is expected. The new measurements have extended the highest average charged-particle multiplicity density per unit of pseudorapidity achieved in pp collisions, roughly a factor five higher than that in average inelastic pp collisions, thereby reducing the multiplicity gap between small and large collision systems. In addition, the results are further compared with previously published measurements and with model calculations obtained using distinct tunes of the PYTHIA 8 Monte Carlo generator, as well as with predictions from the EPOS4. The comparison of the $p_{\rm T}$-integrated K/$π$ and p/$π$ ratios across different collision systems and energies suggests that particle production scales with charged-particle multiplicity, rather than with collision energy or system size. While the PYTHIA 8 tunes and the EPOS4 model are able to reproduce some of these measurements, either quantitatively or qualitatively, none of them consistently describes all observed features of the data.
We investigate the restoration patterns of chiral and $U(1)$ axial symmetries at finite temperature using a soft-wall holographic QCD model. The study employs two distinct parameter sets (Case I and Case II), both calibrated to reproduce a pseudocritical temperature $T_{\rm pc} \sim 155$ MeV and the physical pion mass. The temperature dependence of the light and strange quark condensates confirms a smooth chiral crossover transition, with pseudocritical temperatures of $T_{\rm pc}=0.157$ GeV and $T_{\rm pc}=0.154$ GeV for Cases I and II, respectively. The screening masses of chiral partner mesons ($π$-$σ$ and $η$-$a_0$) become degenerate near $T_{\rm pc}$, providing a clear signature of chiral symmetry restoration. Analysis of the corresponding meson susceptibilities further supports this conclusion. However, the indicator for $U(1)$ axial symmetry restoration, $χ_π- χ_{a_0}$, vanishes at a temperature $T \sim 0.190 $ GeV, which indicates a distinct restoration scale with chiral symmetry restoration scale within the present holographic framework. The temperature-dependent topological susceptibility $χ_{\rm top}^{1/4}$ is also computed, showing a sharp drop near $T_{\rm pc}$ and a subsequent slight decrease. While the model qualitatively captures established features of the chiral transition, the results highlight a limitation in the qualitative description of the $U(1)$ axial anomaly compared to LQCD in our work.
We develop a fully relativistic framework to study the rotational response of gravitationally coupled two-fluid neutron stars within the slow-rotation approximation. Treating the two components as independently conserved perfect fluids interacting only through spacetime curvature, we derive the coupled equilibrium and frame-dragging equations and exploit their linear structure to construct a basis decomposition of the rotational response. This formulation leads to a natural definition of the effective total moment of inertia, which generalizes the single-fluid concept and depends solely on the equilibrium background. It further reveals that the coupled system admits two intrinsic collective rotational eigenmodes, characterized by distinct eigen-moments of inertia, even in the absence of relative rotation between the fluids. Applying this framework to neutron stars containing dark matter, we explore how the presence of an additional gravitationally bound component modifies the global rotational response and its relation to tidal deformability. Our results demonstrate that the persistence or breakdown of rotational-tidal universality in two-fluid neutron stars is governed by dark-sector microphysics rather than by the mere presence of an additional component, and establish a unified framework for interpreting rotational observables, intrinsic mode structure, and universal relations in multi-component relativistic stars.
Strong magnetic fields modify the force that confines quarks inside hadrons and make it direction-dependent. Using quark-antiquark potentials obtained from lattice simulations as inputs to a quark potential model, we investigate how the anisotropic confinement affects the mass spectrum of quarkonium. In the strong-field regime, we find downward mass shifts induced by a softening of the confining potential along the field direction. In particular, the mass shifts of radially excited states are more significant than that of the ground state. For the longitudinal spin eigenstates, the excited-state spectrum strongly depends on the magnetic-field strength, in contrast to the spectrum with conventional isotropic confinement, which is insensitive to the field strength. This provides a clean probe of magnetically induced confinement anisotropy that can be confirmed in future lattice simulations.
We present a new determination of the $\bar b \bar b u d$ ($J^P=1^+$, $I=0$) tetraquark binding energy using lattice QCD with domain-wall light quarks and a nonperturbatively tuned three-parameter anisotropic-clover ``relativistic'' action for the $b$ quarks. We also perform a direct comparison with a reanalysis of data generated in prior work using a lattice-NRQCD action for the $b$ quarks and otherwise identical parameters. Using the new data with relativistic $b$ quarks from seven different ensembles with multiple lattice spacings and pion masses, we perform combined chiral and continuum extrapolations and obtain $(m_{T_{bb}}-m_B-m_{B^*})_{\rm RHQ}=(-79 \pm 23)$ MeV. For the NRQCD data from five ensembles, we perform chiral-only extrapolations and obtain $(m_{T_{bb}}-m_B-m_{B^*})_{\rm NRQCD}=(-74 \pm 17 \pm 10)$ MeV. The lower magnitude of the results obtained here, compared to the original analysis in Phys. Rev. D 100, 014503 (2019), is due to the use of the symmetric parts of the correlation matrices with local four-quark operators only.
The surprising manifestation of collectivity in small collision systems, such as nucleon-nucleon and nucleon-nucleus collisions, is perhaps even more striking when discussed at higher momenta. In larger systems, high-$p_T$ elliptic anisotropy is understood as a selection bias effect due to the smaller energy loss experienced along the shorter direction that aligns with the event plane. However, in small systems the amount of energy loss appears insufficient to reproduce the sizable angular anisotropy observed experimentally. In this work, we explore a new mechanism generating preferred orientations for energetic particles without the need of energy loss. We exploit a simple model that is based on two basic although inalienable ingredients: geometry and quantum mechanics. Our findings suggest that this sum-over-paths mechanism can provide a relevant contribution to so-called flow coefficients of energetic particles traversing deconfined media of any size.