2026-10-06 | | Total: 12
We investigate the longitudinal momentum distributions (LMDs) of the deuteron using a relativistic deuteron wave function derived from a realistic nucleon-nucleon (NN) interaction within the light-front framework. By convoluting the deuteron LMDs with the nucleon parton distribution functions (PDFs), we analyze the unpolarized, longitudinally polarized, and tensor-polarized structure functions of the deuteron. In addition, we also study the EMC effect in the weakly bound deuteron through the structure function ratio $F_{2}^{D}/(F_{2}^{p}+F_{2}^{n})$, which deviates from unity due to nuclear binding and Fermi motion effects.
We present an overview of the relativistic distorted-wave impulse approximation (RDWIA) for lepton-nucleus scattering, with emphasis on its use as an elementary-vertex model in event generators. In the intermediate-energy regime, neutrino-nucleus event generators typically follow a two-step scheme: an elementary neutrino-hadron interaction is first generated, and the outgoing hadrons are then propagated through the nuclear medium using an intranuclear cascade (INC) model. Using exclusive and inclusive $(e,e'p)$ and $(e,e')$ data for several nuclei, we assess how nuclear effects should be distributed between these two stages. In particular, we identify which effects should be included in the elementary-vertex description and which can be treated more effectively by the INC, while addressing the consistency required between both components. We also discuss nuclear transparency in quasielastic scattering. We show that it is directly connected to $(e,e'p)$ data below the two-nucleon knockout threshold; and depends on the full nuclear density distribution and the initial position of the nucleon, therefore, it varies substantially with the shell from which the nucleon is removed. We demonstrate that this shell dependence can be incorporated straightforwardly into the INC component of event generators.
We consider monoenergetic particles propagating in three spatial dimensions and undergoing elastic scattering off an external medium. For an arbitrary rotationally invariant scattering kernel, we derive an exact expression for the frequency- and wave-number-dependent diffusivity, thereby resumming the hydrodynamic gradient expansion to all orders. The result is an infinite continued fraction, with each successive level determined by the relaxation rate of a higher angular harmonic of the collision operator. Thus, while ordinary diffusion probes only the first angular harmonic, the resummed constitutive relation retains information about the full scattering kernel. We show how this microscopic information manifests itself in the analytic structure of the diffusivity and in the response to a localized external source. In particular, deviations from Fick's law decay exponentially with distance, at a rate determined by the nearest complex-wavenumber singularity of the inverse diffusivity. Thus, exponentially small departures from hydrodynamics retain detailed information about the underlying scattering process.
Self-consistent mean-field calculations with semi-realistic interactions have predicted that $N=58$ becomes a magic number near Ni, where the $n2s_{1/2}$ orbit is the highest occupied state. As the last neutrons occupy the shallowly bound $s$-orbit, $^{86}$Ni can be highly distinctive, being a doubly-magic nucleus and possessing a neutron halo simultaneously. With minimal ambiguity in the many-body wave function, it may provide a unique opportunity to investigate, \textit{e.g.}, reaction mechanisms in halo nuclei. As an example, we present differential cross sections for proton elastic scattering predicted via a folding potential without local approximations, and compare them with those obtained using an empirical potential.
Observed massive neutron stars may contain heavier baryons such as $Δ$ isobars and hyperons in their interiors. In the presence of $Δ$ isobars, the neutrino emitting nucleonic Urca processes can be accompanied by $Δ$-direct-Urca processes involving nucleons and hyperons. Over the past three decades, significant theoretical effort has been dedicated to understanding the thermal evolution of neutron stars with the nucleonic and hyperonic Urca processes emerging as dominant cooling mechanisms. While recent studies have proposed the presence of $Δ$ isobars in the dense core of neutron stars, the specific contribution of the $Δ$-direct-Urca process to neutrino emissivity has not been explicitly quantified. In this work, we address this gap by deriving the neutrino emissivity associated with $Δ$-direct-Urca processes in dense matter inside neutron stars. We show that once $Δ$-direct-Urca processes become kinematically allowed, their emissivities are comparable to those of nucleonic-direct-Urca processes.
We systematically investigate the correlations between nuclear fission properties and the surface energy coefficient $a_{\rm surf}$ extracted from semi-infinite nuclear matter. To this end, we employ the covariant density functional theory based on a point-coupling energy density functional (EDF). Taking $^{240}$Pu and $^{236}$U as representative examples, we employ Gaussian-process emulators to construct the joint posterior distributions of $a_{\rm surf}$ and selected fission observables, constrained by the ground state properties of the two nuclei. The results reveal a strong correlation between $a_{\rm surf}$ and the fission barrier heights. A sensitivity analysis further shows that the height of the outer fission barrier is governed predominantly by the isoscalar-vector derivative coupling constant $δ_V$ in the point-coupling EDF, whereas that of the inner barrier is by the combined effects of $δ_V$ and several mutually correlated coupling constants.
We develop a new model that provides a unified description of the equilibration and subsequent transport of conserved charges, including baryon number, electric charge, and strangeness, as well as the dynamical formation and evolution of the quark-gluon plasma medium in high-energy nuclear collisions. This is achieved by extending the state-of-the-art model of dynamical core-corona initialization to incorporate phenomenological source terms for conserved charges within a full (3+1)-dimensional hydrodynamic framework; the updated model is referred to as DCCI3. Consequently, the model simultaneously describes the core (equilibrium) and corona (nonequilibrium), including their mutual interactions. We present a detailed construction of the DCCI3 model and employ it in numerical simulations of Pb+Pb collisions at $\sqrt{s_{NN}}=2.76~\mathrm{TeV}$ to investigate the resulting dynamics of matter and conserved charges, with particular emphasis on the longitudinal structure. We find that each of the chemical potentials exhibits sizable fluctuations around midrapidity, which challenges the conventional assumption of vanishing values as an oversimplification. Moreover, the baryon chemical potential increases sharply with rapidity and becomes comparable to the temperature at forward rapidities. This suggests that the Rapidity Scan (RS) approach in high-energy nuclear collisions can serve as a complementary probe to the Beam Energy Scan (BES) in exploring the QCD phase diagram.
The mass table in the deformed relativistic Hartree-Bogoliubov theory in continuum (DRHBc) with the PC-PK1 density functional has been established for nuclei with $8 \leq Z \leq 120$, extended from the previous works for even-even nuclei [Zhang et al. (DRHBc mass table collaboration), At. Data Nucl. Data Tables 144, 101488 (2022)] and for even-$Z$ nuclei [Guo et al. (DRHBc mass table collaboration), At. Data Nucl. Data Tables 158, 101661 (2024)]. The calculated binding energies, two- and one-nucleon separation energies, root-mean-square (rms) radii of neutron, proton, matter, and charge distributions, quadrupole deformations, neutron and proton Fermi surfaces, and the blocked neutron (proton) orbitals of odd-$N$ ($Z$) nuclei are tabulated and compared with the available experimental data. A total of 9495 nuclei are predicted to be bound, with an rms deviation of 1.444 MeV from the 2380 mass data. Good agreement with the available experimental pairing gaps, $α$ decay energies, and charge radii is also achieved. The accuracies of the calculated nuclear masses and nucleon separation energies as well as the prediction for drip lines are compared with those obtained by other relativistic and nonrelativistic density functional calculations. It turns out that the DRHBc theory with PC-PK1 provides one of the best microscopic descriptions for nuclear masses. The systematics of nucleon separation energies, pairing gaps, pairing energies, two-nucleon gaps, $α$ decay energies, rms radii, quadrupole deformations, potential energy curves, neutron density distributions, and neutron mean-field potentials are discussed.
Nucleon scattering off strongly deformed nuclei such as the actinides requires coupled-channels (CC) calculations in which several members of the ground-state rotational band are coupled explicitly, which makes Bayesian uncertainty quantification of deformed global optical potentials costly. We reformulate the direct boundary matching method (DBMM) by imposing its outgoing-wave boundary condition through a generalized Bloch operator; the resulting Bloch-DBMM solver converges as fast as the modern R-matrix method and extends to charged projectiles on targets with deformed charge distributions. On top of it we build a reduced-basis emulator, CCB-ROM, whose basis vectors span all channels of a $J^π$ block simultaneously, and generalize the coupled-channels empirical interpolation method to arbitrary non-affine deformed interactions. For $n,p+{}^{232}$Th inelastic scattering with the ground-state band coupled through the $8^+$ state, at 4 and 35 MeV, the Bloch-DBMM solver is benchmarked against Frescox, and the emulator replaces a 2700-dimensional system by one of dimension $N_b\le200$, with a cost of $\mathcal{O}(N_b^3)$ independent of the number of channels, reproducing the exact solver calculations for all five band members with speed-ups of over two orders of magnitude. CCB-ROM thus delivers differential cross sections, transmission coefficients and reaction cross sections from one calculation, at a cost compatible with Bayesian calibration of deformed optical potentials for neutrons and protons alike, and extends naturally to closed channels and to compound-nucleus reaction codes.
We present two complementary neural routes for microscopic neutron-star EOS inference: truncated sequential neural posterior estimation with mixture importance sampling (TSNPE+MIS), which targets one observation, and amortized neural posterior estimation with importance sampling (A-NET+IS), which is reused across nuclear inputs and NICER sources. Both are corrected against the original prior and full likelihood, with the effective sample size as a diagnostic. We test nucleonic and hyperonic relativistic mean-field models with nuclear data, NICER, GW170817, and pQCD. For the fiducial analyses, both agree with UltraNest in the parameter posteriors and 90% mass-radius and tidal-deformability bands. TSNPE+MIS band-edge differences are only 0.024-0.033 km and 1.02-1.03%, while A-NET+IS differences are 0.019-0.050 km and 0.86-1.13%. TSNPE+MIS reduces the computing time by factors of 1.6-2.3; after training, A-NET+IS is 15-240 times faster per configuration. The frozen A-NET also analyzes three NICER sources absent from training, including the newly reported high-mass, compact PSR J1614-2230 posterior. For this source, the inferred $R_{1.4}$ agrees with UltraNest within 0.010 km for nucleonic matter and 0.009 km for hyperonic matter. We also introduce a Green-function Evidence Network (EN), a regression-based variant of Evidence Networks and, to our knowledge, the first Evidence Network applied to neutron-star EOS inference. After a one-time cost, its network queries return absolute evidences in about 0.03 s for new nuclear data and 3.5 s for a new NICER source; all 16 differences from UltraNest are below $1σ$, and independent importance-sampling evidences agree within the quoted EN uncertainties. Under the adopted likelihood, both evidence calculations favor nucleonic matter by Bayes factors of about 22:1-27:1.
We show that a collisionless point-vortex gas behaves as a genuine Vlasov system. Starting from the non-canonical Hamiltonian structure of point-vortex dynamics, we derive a reduced Klimontovich equation and the corresponding Vlasov description directly on the physical plane, without introducing auxiliary momentum variables. The theory predicts a rich spectrum of collective phenomena, including vortex Rossby-like waves, global oscillations, resonant critical layers, and a dielectric response operator analogous to that of plasmas and self-gravitating systems. For smooth vortex distributions, differential rotation gives rise to a vortex analogue of Landau damping, whereby collective modes transfer energy and angular momentum to resonant vortex annuli in the absence of collisions or dissipation. Remarkably, differential rotation acts as a stringent mode-selection mechanism: for a Gaussian vortex cloud, only a single quadrupolar quasi-mode survives, while all higher multipoles are destroyed by phase mixing. We further predict a vortex counterpart of the beam--plasma instability, in which a thin vortex ring orbiting an Onsager patch selectively excites Kelvin modes with an azimuthal symmetry controlled by the ring population. These results establish kinetic theory as a natural framework for collective vortex dynamics and suggest new mechanisms for vortex transport, vortex avalanches, and rotational glitches in neutron-star superfluids.
Relativistic heavy quarks equilibrate even when their microscopic drag and momentum diffusion violate the Einstein relation. We investigate how equilibration proceeds as a function of the coupling and compared to simplified Fokker-Planck models. To do this, we present a detailed calculation of the momentum-transfer kernel in weakly coupled non-Abelian plasmas through strict $\mathcal{O}(g^4)$, and study the resulting dynamics at leading order in the inverse heavy-quark mass. Beyond leading logarithm, its non-Gaussian structure plays a crucial role in equilibration. We examine closely the analytic structure of the evolution kernel, determine the large-order asymptotics of the momentum-transfer cumulants, and show how the kernel's bounded analytic domain fixes asymmetric exponential tails at the level of the momentum transfer probability. Strongly coupled $\mathcal{N}=4$ SYM also features an approximately Gaussian core and exponential tails, but differs in the nature of the kernel's singularities and the velocity dependence of the transport coefficients and tail exponents. These differences have distinct dynamical consequences: for the parameters studied, the persistence of the relativistic population in steeply falling spectra is largely explained by drag alone at very weak coupling, whereas fluctuation-induced survivor bias produces a large relative enhancement at strong coupling. Nevertheless, the shared qualitative structure between weak and strong coupling points the way towards a practical phenomenological description of heavy quarks anchored in first-principles field theory.