2026-10-28 | | Total: 20
In this work, we focus on investigating the impact of the short-range correlations on the properties of quarkyonic matter. These correlations induce the scattering of nucleons to momenta exceeding the Fermi momentum, thereby altering the momentum-space occupation distribution that underlies the quarkyonic framework. In particular, we examine to what extent the emergence of a high-momentum tail, due to the short-range nucleon-nucleon correlations, in the momentum distribution, for momenta exceeding the Fermi momentum, affects the fundamental properties of quarkyonic matter, such as the energy density, pressure, and speed of sound, among others. In fact, at high densities, the behavior of the momentum distribution has been established, both experimentally and theoretically, to follow a characteristic $\sim 1/k^4$ dependence at large momenta. To the best of our knowledge, this effect on quarkyonic matter has not yet been systematically investigated. Moreover, this experimentally confirmed feature of the momentum distribution suggests that heavy-ion collisions may probe quarkyonic matter, underscoring the importance of incorporating short-range correlations into a complete quarkyonic model.
We study the covariant Wigner function for the spin-$3/2$ Rarita-Schwinger field and its relation to the spin density matrix and the mean spin polarization vector. We obtain a formula relating the covariant Wigner function to the mean spin vector for a weakly interacting field. We then calculate the mean spin vector induced by thermal vorticity, thermal shear, and spin potential at local thermal equilibrium in the Belinfante, canonical and GLW pseudo-gauges. The resulting formula can be used to predict the spin polarization of spin-$3/2$ hadrons, such as the $Ω^-$, in heavy-ion collisions. Differences in the thermal-shear induced polarization in the canonical and Belinfante pseudo-gauges make it possible to use local spin polarization measurements of the $Ω^-$ to discriminate between the two pseudo-gauges.
We analyze the on-shell content of the local Lorentz covariants of the five-dimensional spin-zero Duffin-Kemmer-Petiau algebra in elastic spin-zero-spin-$1/2$ scattering. The scalar, vector, and tensor structures are reduced explicitly to the two independent amplitudes allowed by the physical parity-conserving process. The two scalar structures are proportional on shell, the genuine nonminimal vector $i[P,β^μ]$ is longitudinal and vanishes after contraction with the conserved fermion vector current, and the projected vectors $Pβ^μ$ and $β^μP$ each reduce to one half of the minimal-vector contribution. The tensor covariants introduce no additional on-shell directions. We discuss the consequences for the four invariant choices used in meson-nucleus optical-potential calculations and show that different nuclear equations can arise only after on-shell-equivalent representatives are promoted to optical operators and iterated beyond the order at which the elementary impulse approximation is fixed. These results provide a systematic basis for distinguishing genuinely independent elementary amplitudes from representation-dependent effects generated by the optical embedding.
The process of $β$-delayed fission ($β$DF) consists of a mother nucleus $β$-decaying to different states of a daughter isotope, followed by fission of the latter with a probability that strongly depends on the energy of the states populated in the $β$ decay. Studying this process experimentally can be challenging due to the typically small $β$DF probabilities. The theoretical framework presented in this work was developed to globally characterise $β$DF probabilities. This approach uses energy density functional based models to obtain the necessary inputs for the final probabilities, which are obtained through the nuclear reaction code TALYS. Calculations were performed for all isotopes with an experimentally measured value of $β$DF probabilities. Following the optimisation of the theoretical framework to better reproduce known data, a global RMS deviation of about two orders of magnitude was achieved, an accuracy that can be used to identify promising cases for future experimental $β$DF studies.
Antineutrinos produce hyperons without kaons through Cabibbo-suppressed quasi-elastic scattering on a single nucleon. A two-body current acting on a correlated nucleon pair, nubar [NN] -> l+ Y N', gives the same hyperons plus a second nucleon. A companion paper constructs this current without fixing its size. StrangeMC, the only generator with the channel, set it to a flat 15% of the quasi-elastic hyperon cross section. We determine the ratio f = sigma_2p2h/sigma_QE of the two-body to the quasi-elastic hyperon cross section on carbon versus antineutrino energy, and the detector signatures of two-body events. The strange two-body current and its non-strange counterpart are computed on the same nucleon pair. Their ratio is converted into f with the Valencia non-strange 2p2h and quasi-elastic cross sections for carbon, assuming that nuclear effects missing from the pair are common to both. Under this assumption f rises from 0.14 at 0.6 GeV to 0.20-0.22 at 2-5 GeV. Averaged over the Booster antineutrino spectra it is 0.15-0.17, close to the 15% placeholder, which StrangeMC now replaces with f(E). The largest evaluated uncertainty comes from separating the two-body current from one-body production followed by rescattering: across the separation prescriptions explored, f at 1 GeV ranges from 0.12 to 0.22. Interference with one-body production is not calculated and could be as large as the two-body rate. In a simulated SBND antineutrino-mode exposure on argon with the carbon f, two-body events are about 11% of events that produce a strange hadron. They supply more than a quarter of the kaon-less Sigma+ events used to tag hyperon charge exchange, partly as Sigma+ made at the weak vertex. They reconstruct 12-30% low in energy under the quasi-elastic hypothesis. After final-state interactions half of them have a proton above 40 MeV at the vertex, against a fifth of quasi-elastic events.
We present a microscopic description of the α-decay processes from axially-symmetric superfluid nuclei in terms of proton (π) and neutron (ν) single-particle (sp) degrees of freedom within a standard Woods-Saxon mean field (WSMF) and a density-dependent clustering potential centered around the average Mott point of the nuclear density, where the α-particle is formed. In this way the clustering process is automatically enhanced on the nuclear surface. As an illustration of the method, a systematics of α-decay from even-even emitters is presented.
The equation of state of strongly interacting matter plays a %is one of the central quantities characterizing central role in characterizing quark--gluon plasma created in relativistic heavy-ion collisions. It has recently been recognized that the QCD speed of sound is responsible for the linear correlations between the mean transverse momentum and charged-particle multiplicity in ultra-central collisions. In realistic collisions, however, this thermodynamic correlation is inevitably biased by quantum-initiated event-by-event fluctuations. In this work, we present a novel statistical framework to correct this bias by leveraging the Gaussianity condition of the quantum fluctuations. Furthermore, we demonstrate that this approach can also serve as a probe of thermalization. Validated by the event-by-event hydrodynamic simulations, the extracted speed of sound in our framework is successfully consistent with the predictions from lattice QCD, from large to small collision systems.
Spin-orbit coupling makes the mean position of a spinning particle transform anomalously under a Lorentz boost, with an additional shift known as a side jump. In this work, we develop a trace method to derive an exact but compact transformation kernel for the probability center of a general three-dimensional Dirac wave packet under an arbitrary finite Lorentz boost. The expression includes the physics of finite momentum width, spin-momentum correlations, and coherent phase structure without a semi-classical expansion. We identify texture-dependent side jumps even at fixed initial mean momentum, spin, and orbital angular momentum, as well as a non-monotonic mass dependence for Gaussian packets at fixed momentum width. The explicit kernels further resolve the spin and orbital angular-momentum contributions for massive and massless fermions, providing a kinematic basis for covariant particle-based spin transport that retains wave-packet information.
We construct a density-dependent covariant energy density functional (CE-ME1) from full-Dirac-space relativistic Brueckner--Hartree--Fock (RBHF) self-energies and the corresponding equation of state. The space-like component $U_V$ is incorporated at the Fermi surface through reduced scalar and time-like self-energies, which are mapped onto local Hartree $σ$, $ω$, $δ$, and $ρ$ channels. The density-dependent couplings are represented by cubic Chebyshev polynomials over the density range $0.08\leqρ\leq0.48~{\rm fm}^{-3}$. A constrained thermodynamic calibration uses 199 RBHF energy-per-particle values from published microscopic computations at six isospin asymmetries, while maintaining positive, monotonic couplings. The Bonn-A RBHF saturation point, incompressibility $K_\infty=258$ MeV, symmetry energy $J=33.1$ MeV, and slope $L=65.2$ MeV are reproduced. The predicted curvature $K_{\rm sym}=33$ MeV differs from the published RBHF value of $-67$ MeV, illustrating the derivative sensitivity of the local reduction. The combined energy RMS is $0.613$ MeV, measuring model discrepancy rather than statistical uncertainty. Uniform matter remains mechanically stable and causal above saturation within the calibrated range. In cold charge-neutral $npeμ$ matter, CE-ME1 reasonably reproduces the interpolated RBHF pressure and proton fraction over most of the interval, although the pressure discrepancy increases near its upper edge. The DD-ME2 benchmark is stiffer and predicts a lower proton fraction at high density.
The nuclear equation of state shapes both the structure of rapidly rotating neutron stars and the orbital motion in their surrounding spacetime. Its uncertainty affects the mass and radius at which the innermost stable circular orbit reaches the stellar surface, a condition that defines the transition between disks whose inner edge is set by the surface and those truncated by orbital instability. We quantify how measurements of the electric dipole polarizability of neutron-rich nuclei reduce the uncertainty in this transition and in the orbital properties at the disk inner edge. We compute uniformly rotating neutron stars in general relativity for a family of equations of state that systematically vary the symmetry energy, derived from relativistic energy density functionals with density-dependent point-coupling interactions. We compare the results for the full family with those for the subset consistent with the measured dipole polarizabilities of finite nuclei. The dipole polarizability constraint reduces the spread in transition mass and equatorial radius by factors of about seven to ten across the spin frequencies considered. Along this transition, it also confines the orbital frequency at the disk inner edge to a band narrower than $20$~Hz at each spin. Combined with a measured neutron-star spin, a kilohertz quasi-periodic oscillation interpreted as disk orbital motion further restricts the mass--radius region, allowing an independent mass estimate to probe the symmetry-energy constraint. The dipole response of finite nuclei thus informs both the structure of rotating neutron stars and the orbital properties of their inner accretion disks.
We investigate how sensitive the effects of neutrino trapping on hot hybrid stars are to uncertainties in the hadronic and quark equations of state and in the thermodynamic construction of the hadron--quark transition. The hadronic phase is described by finite-temperature covariant density functional equations of state, and the quark phase by a three-flavor Nambu--Jona-Lasinio (NJL) model with repulsive vector interactions and two-flavor color superconductivity. Because confinement is absent from the NJL model, an effective pressure shift $B^*$ is used phenomenologically to locate the deconfinement transition. We construct isentropic stellar equations of state with and without trapped electron neutrinos, vary the quark vector coupling and the nuclear symmetry-energy slope, and compare global with local conservation of electron-lepton number within phase coexistence. Neutrino trapping generically shifts deconfinement toward higher density throughout the parameter range considered. Increasing the vector interaction stiffens quark matter and raises the maximum stellar mass, but only moderately for the couplings studied. We also examine stronger diquark pairing and RG-consistent regularization, instead of sharp momentum cut-off: both affect the pairing gap, phase boundaries, and quark-core composition more strongly than they alter the qualitative trends of the stellar mass--radius sequences. Conversely, trapping reduces the sensitivity to the hadronic symmetry-energy slope because the enhanced proton fraction drives matter toward greater isospin symmetry. Global and local lepton-number prescriptions produce rather similar mass--radius relations but substantially different internal phase structures: global conservation leads to a considerably more extended mixed phase.
Forward-backward multiplicity correlations probe the longitudinal structure of particle production, but fluctuations in the number of emitting sources complicate their interpretation. Strongly intensive observables remove the contribution of source-number fluctuations at second order. Here we extend the construction to all cumulant orders. For independent identical sources whose emission law does not depend on the source count, we build a normalized hierarchy of strongly intensive observables with $n$ new algebraically independent quantities at order $n$. Every observable is normalized to unity for statistically independent particles, vanishes when both window multiplicities are fixed, and remains finite at equal positive window means. We give explicit expressions through fourth order and an all-order recursion. Analytical particle production models and PYTHIA calculations for $pp$ collisions at $\sqrt{s}=0.9$ and 13 TeV demonstrate the sensitivity to window geometry, event selection, and connected multiparticle correlations. Fits constrained by the mean multiplicities and second-order observables provide conditional higher-order benchmarks. Deviations of these benchmarks from the PYTHIA results test the assumed correlation hierarchy, subject to the quality of the lower-order fit. Preliminary ALICE measurements at $\sqrt{s}=13$ TeV are used to illustrate the corresponding model-dependent scenarios.
This work investigates the properties of spin polarized strange quark matter SQM and the resulting structure of strange quark stars SQSs within a thermodynamically self consistent quasiparticle framework. By developing a dynamic bag constant intrinsically coupled to the medium dependent reduction of quark effective masses we ensure strict internal thermodynamic consistency. This coupling gives rise to a distinct non monotonic behavior in the speed of sound significantly stiffening the Equation of State EoS at high densities. Consequently our model predicts maximum stable gravitational masses up to 2.19 Mo. These results demonstrate that hot magnetized SQM provides a robust and viable mechanism to explain ultra massive pulsars such as PSR J0740 6620 and PSR J0952 0607 without requiring unrealistic thermal or structural assumptions.
Multiparticle correlations of event-plane angles, known as symmetry plane correlators (SPCs), and correlations of anisotropic flow amplitudes of different harmonics, known as symmetric cumulants (SCs), help constrain the initial conditions and the transport properties of the system formed in heavy-ion collisions. Such multiparticle correlations are robust against nonflow biases in heavy-ion collisions; however, in small systems the nonflow contribution is sizeable and can distort SPCs and SCs at low multiplicity, mimicking a flow-like signal. Evidence of nuclear-geometry-induced collectivity in light-ion collisions, such as OO and Ne--Ne, has been established at the LHC. Measurements of SPCs and SCs in these systems, and their comparison to heavy-ion collisions, can therefore clarify the system-size dependence of these correlations and the strength of nonlinear mode coupling. In this work, an improved nonflow subtraction framework based on multiplicity reweighting ($\czero$) and $v_1$-estimator-based ($\cone$) scaling is validated for selected SPCs and SCs using a controlled toy model and the \textsc{Pythia~8~Angantyr} model, where collective flow is absent by construction. We find that $\cone$ provides the more reliable nonflow control and apply it to predict the SPCs and normalized SCs in OO collisions at $\sqrt{s_{\rm NN}}=5.36$~TeV with the \textsc{ampt} string-melting model. The results are compared with published Pb--Pb measurements from ALICE and from \textsc{ampt} as a function of midrapidity charged-particle multiplicity $\langle N_{\rm ch}\rangle(|η|<0.5)$, placing light- and heavy-ion systems on a common system-size axis.
Initial stages of relativistic ion and proton-proton collisions can be described within the framework of color flux tube formation between the receding color charges. In the approximation of fully penetrable color string volumes, regions of high string density yield overlaps of randomly oriented color fields and create string clusters. The paper discusses two ways of characterizing the medium formed by the resulting field configurations. One is based on the chromo-electric field strength that scales as the square root of the number of color sources in the cluster. The other is defined by the energy density of the confined color field and thus scales linearly with that number. We compare the effect of the two approaches on the extracted effective temperature of the medium, which is the input for calculations of collisional and radiative energy loss of a probe traversing it.
While the body of experimental nuclear data continues to grow at a rapid pace, significant gaps persist. The lack of evaluated experimental data is far more pronounced in charged-particle reactions which are essential to a large portion of the isotope production community for vital medical, industrial, and national security interests. In the absence of experimental and evaluated data, isotope production relies on theoretical models that may result in low yield and isotopic impurities. This work investigates model and parameter adjustments in TALYS 2.0 to establish uncertainty bands on proton-induced cross sections, highlight influential parameters in calculations, and constrain adjustments within physically defensible bounds. The EXFOR dataset for proton-induced reactions on natural copper, comprised of 24 independent cross sections, was used for iterative Bayesian model updating. Upon convergence, parameters and their uncertainties were compared against 23 cumulative cross sections from EXFOR, highlighting the limitations of generalization. Various machine learning models were leveraged to predict TALYS performance, with a small neural network showing useful ranking capabilities across several generations, suggesting the use of a surrogate approach to screening for future evaluations.
We study how magnetic fields affect the onset of quark saturation, which limits the domain of applicability of a purely baryonic description. We consider matter composed of electrically neutral $ddu$ baryons, including neutron- and $Δ^0$-like states. Using a simple constituent-quark model, we compute the baryon wave functions and spectra in magnetic fields. The magnetic field polarizes the spin-flavor structure, compresses the internal wave function in the transverse direction, and strongly suppresses the transverse kinetic energy. From the baryon wave functions, we construct the one-body reduced density matrix and derive the corresponding single-quark distributions in longitudinal momentum and Landau-level space. We then use these distributions to evaluate quark occupation probabilities in an ideal baryon gas and determine the onset of quark saturation. At finite magnetic field in the $+z$-direction, the $d_\downarrow$ quark in the lowest Landau level reaches saturation first. The resulting saturation density depends nonmonotonically on the magnetic-field strength: it decreases at weak fields and increases at stronger fields. This behavior reflects the competition between spin-flavor polarization, which promotes saturation, and transverse spatial squeezing, which delays saturation by broadening of the quark transverse momentum distribution.
We report a benchmark study of the electric-quadrupole hyperfine interaction in the electronic ground state of the RaF molecule, including finite-nuclear-size effects. The calculations include relativistic effects using the Dirac-Coulomb-Breit Hamiltonian and electron correlation up to perturbative quadruple excitations within the coupled cluster framework. At the achieved level of accuracy, penetration of the electronic density into the finite nuclear quadrupole distribution produces a non-negligible correction to the usual electric field gradient (EFG) description at the center of the nucleus and must be considered. At a fixed rms charge radius, the EFG is much less sensitive to the shape of the monopole nuclear charge distribution: electron wavefunctions computed using the Gaussian and Fermi nuclear charge models yield nearly identical EFGs. The recommended effective EFG parameter for the $X\,{}^{2}Σ_{1/2}$ electronic ground state is $\langle q\rangle_0=-5.178(21)$ a.u., with an estimated uncertainty reduced by an order of magnitude compared to previous studies. Combined with future measurements of electric-quadrupole hyperfine coupling constants in RaF, this result can be used to determine the nuclear quadrupole moments of Ra isotopes with record precision.
The beta decay of neutron-rich nuclei probes the microscopic nature of weak interactions in the nuclear medium under extreme conditions. Such weak decays drive the synthesis of heavy elements in astrophysical environments, but remain poorly understood. Far from stability, experimental data are scarce, placing strong demands on ab initio nuclear theory. Two-body weak currents have been proposed as the dominant mechanism responsible for the universal quenching of theoretically predicted Gamow-Teller decay strength. However, direct experimental validation in neutron-rich nuclei has been lacking. Here we measure the beta decay of 54K, a nucleus at the experimental frontier of radioactive ion-beam facilities. In this work, the doubly magic nature of the daughter nucleus 54Ca results in a remarkably simple decay pattern, dominated by a single transition to a neutron-unbound state at 6.2 MeV with a comparative partial half-life logft=4.4(1) or Gamow-Teller decay strength B(GT)=0.26 (+0.08) (-0.04). Our measurement is reproduced by new state-of-the-art calculations of the decay, providing direct evidence for the essential role of two-body currents in nuclear beta decay for systems with a significant neutron excess. We establish a key benchmark for theoretical predictions of weak processes in exotic nuclei.
We investigate axion-like particle (ALP) production in the interior of Betelgeuse via nuclear de-excitation of thermally populated excited states of $^{23}$Na and $^{55}$Mn. As these ALPs propagate to Earth, they can produce photons either through their conversion in the Galactic magnetic field or through their decay. Using INTEGRAL-ISGRI and NuSTAR observations, we search for these photon signals, find no significant excess, and derive 95% CL exclusion limits on the ALP couplings. For magnetic-field conversion of ALPs, the $^{23}\mathrm{Na}$ ($^{55}\mathrm{Mn}$) line observed with INTEGRAL-ISGRI constrains $|g_{aγ}g_{ap}| < 4.1\times10^{-21}$ ($6.9\times10^{-19}$) $\mathrm{GeV}^{-1}$ for $m_a < 10^{-10} \mathrm{eV}$. For ALP decays, the $^{23}\mathrm{Na}$ ($^{55}\mathrm{Mn}$) line observed with INTEGRAL-ISGRI (NuSTAR) constrains $|g_{aγ}g_{ap}|$ down to $3.0\times10^{-22}$ ($2.9\times10^{-19}$) $\mathrm{GeV}^{-1}$ at $m_a \simeq 400$ ($83$) keV, weakening in both cases as $m_a^{-3/2}$ for lower ALP masses, and as $m_a^{-2}$ for $^{23}$Na below $m_a \sim 1$ keV, with sensitivity extending to $m_a \simeq 440.2$ ($124$) keV. Although our $^{23}\mathrm{Na}$-based limits are the strongest for both channels, their larger transition energy makes them more sensitive to the precise evolutionary stage of Betelgeuse, whereas the $^{55}\mathrm{Mn}$-based limits are considerably more robust against such uncertainties.