PaperPanorama

Nuclear Theory·nucl-th

Tuesday·December 30, 2025

20 papers12 primary·8 cross-listed

  1. 01

    Bidirectional Neural Networks for Global Nucleon-Nucleus Optical Model Calculations

    Jin Lei

    Modern nuclear data evaluation increasingly requires not only accurate scattering calculations, but also efficient methods for uncertainty quantification and parameter optimization, tasks that benefit from differentiable solvers amenable to gradient-based algorithms. I present a neural network emulator based on Bidirectional Liquid Neural Networks (BiLNN) that provides a fully differentiable mapping from optical potential parameters to scattering wave functions. The key innovation enabling generalization across the parameter space is the use of phase-space coordinates that normalize the oscillation wavelength regardless of projectile energy, allowing a single network to span 1 to 200~MeV. Trained on Numerov solutions for twelve target nuclei (\nuc{12}{C} to \nuc{208}{Pb}), both protons and neutrons, and partial waves up to , the network achieves an overall relative error of 1.2\%. The predicted wave functions yield accurate -matrix elements and elastic scattering cross sections, reproducing diffraction patterns spanning four orders of magnitude. Importantly, the model extrapolates successfully to nuclei not included in training (\nuc{24}{Mg}, \nuc{63}{Cu}, \nuc{184}{W}) with comparable accuracy, demonstrating that it has learned the physics of the optical model rather than memorizing specific targets. The differentiable nature of the trained model opens the door to gradient-based optimization of optical model parameters and efficient uncertainty quantification.

    nucl-thPRC(2026)·3 citations
  2. 02

    Exploring Interplays Between Neutron Superfluid Vortices and Proton Fluxtubes in the Outer Core of Neutron Stars

    Tatsuhiro Hattori · Kazuyuki Sekizawa

    In the outer core of neutron stars, P superfluid neutrons and S superconducting protons are deemed to exist, forming quantum vortices and magnetic fluxtubes, respectively. Those quantum vortices and fluxtubes play an important role in explaining observed sudden changes of rotational frequency, known as pulsar ``glitches.'' While the most of conventional glitch models rely on pinning/unpinning dynamics of neutron superfluid vortices in the inner crust, contributions of the outer core have not been ruled out. However, the latter possibility has been less explored so far and further thorough investigations are desired. In this study, we are thus developing a microscopic model based on spin-2 Gross-Pitaevskii equation (GPE) for neutron P superfluid vortices coupled with Ginzburg-Landau equation (GLE) for fluxtubes associated with superconducting protons. In this contribution, we outline our theoretical framework and report tentative results showing how shape of quantum vortices could be affected by the presence of a proton fluxtube.

    nucl-thastro-ph.HEcond-mat.quant-gasEPJ Web Conf.(2026)·3 citations
  3. 03

    Bound States of Baryons in Light Nuclei

    Igor Filikhin🇺🇸 · Roman Ya. Kezerashvili🇺🇸 · Branislav Vlahovic🇺🇸

    We investigate bound states of light -clusters (), motivated by the potential recently developed by the HAL QCD collaboration. To regularize this potential, we remove the deeply attractive core at and parametrize the long-range component () using a two-range Gaussian form. This procedure preserves the relevant two-body bound state energy while having a negligible effect on the and systems. An effective potential is then constructed by fitting a two-range Gaussian function to the long-range component of the folding potential, enabling calculations of the bound state energies of the , , and systems. The regularization procedure leads to a substantial reduction in bound state energies compared to those obtained with the original potential. We further extend the analysis to -cluster systems by introducing an interaction, derived by comparing the existing and potentials. Our results suggest that several parametrizations predict bound states in -containing clusters. Finally, the interaction is described using a contact-like potential approach, motivated by the effective field theory.

    nucl-thPRC(2026)·2 citations
  4. 04

    Elucidating the role of the surface energy in density functional theory

    Md Jafrul Islam · Athul Kunjipurayil · J. Piekarewicz · A. Volya

    The saturation of symmetric nuclear matter -- reflected in the nearly constant interior density of heavy nuclei -- is a defining property of nuclear matter. Modern relativistic energy density functionals (EDFs) calibrated exclusively to the properties of finite nuclei, make robust predictions with quantified uncertainties about the bulk properties of symmetric nuclear matter in the vicinity of the saturation density. Following the same fitting protocol, nonrelativistic Skyrme EDFs systematically predict higher saturation densities than their relativistic counterparts. To investigate this tension in the bulk limit, we study the ground-state properties of hypothetical symmetric macroscopic nuclei containing thousands of nucleons. Using both relativistic and non-relativistic EDF frameworks, we extract the corresponding liquid-drop parameters. We find a clear correlation between the volume and surface energy coefficients: Skyrme models, which saturate at higher densities, develop softer and more diffuse surfaces with lower surface energies, whereas relativistic EDFs, which saturate at lower densities, produce more defined and less diffuse surfaces with higher surface energies. This compensating behavior allows both classes of models to reproduce empirical nuclear radii despite their distinct saturation properties. Our analysis suggests that the apparent disparity in saturation densities arises from the intrinsic balance among saturation density, bulk binding energy, and surface tension, rather than from the fitting protocol.

    nucl-thPRC(2026)·0 citations
  5. 05

    Kinematic and dynamical origins of mean- fluctuations in heavy-ion collisions

    Lipei Du🇺🇸

    Event-by-event fluctuations of the mean transverse momentum (mean-) provide a sensitive probe of collective dynamics beyond single-particle spectra and anisotropic flow. We present a systematic study of mean- fluctuation observables using a Bayesian-calibrated multistage hydrodynamic framework, including quantitative comparisons to RHIC measurements and model-based investigations of beam-energy and kinematic-acceptance effects. The experimental definitions employed by the STAR and ALICE Collaborations are implemented explicitly and found to yield consistent results within controlled limits. We study the centrality and beam-energy dependence of the observable, its sensitivity to key soft-sector ingredients, and the impact of the kinematic acceptance. By introducing scaled- cuts, we demonstrate that a part of the apparent energy dependence arises from kinematic projection effects, while the remaining trends reflect genuine collective dynamics. Our results establish mean- fluctuations as a nontrivial and independent validation of calibrated hydrodynamic descriptions of the quark--gluon plasma.

    nucl-thnucl-exPLB(2026)·3 citations
  6. 06

    Factorized distorted wave calculations for electron, neutrino, and BSM processes

    Alexis Nikolakopoulos🇺🇸 · Raúl González-Jiménez🇪🇸

    We point out that, under certain conditions, the nuclear currents that couple to vector bosons can be written as the trace of the product of two matrices. One contains `nucleon dynamics', e.g. form factors, the other contains the overlaps of nuclear wavefunctions. This factorized form may always be obtained and viewed as a `local' approximation, in which particle four-momenta that enter in the transition operator are fixed to their asymptotic values. We write the overlap matrix in a general form in terms of Dirac matrices. The current with arbitrary couplings to the nucleon can then be evaluated using standard trace identities. We show that this encompasses non-relativistic models as well. We have tabulated overlaps obtained in the relativistic distorted-wave impulse approximation, which can be used to compute the single-nucleon knockout cross section. We give a self-contained overview of the formulae. We discuss some properties of the overlap matrices which may be derived from general principles, and highlight differences with the commonly used plane-wave impulse approximation. The factorized form is attractive for (neutrino) event generators: it abstracts away the nuclear model and allows to easily modify couplings to the nucleon. This allows to consistently treat electron, neutrino and beyond Standard Model (BSM) processes.

    nucl-th2 citations
  7. 07

    The effect of inversion of and orbits on halo formation in heavy sodium isotopes

    Jagjit Singh · J. Casal · L. Fortunato · N. R. Walet

    The role of the inversion of the and shell-model orbits in the emergence of halo structures in the ground states of neutron-rich Na is investigated. Families of two- and three-body models are constructed with effective core-neutron interactions, with parameter choices based on a combination of the available experimental data and systematic trends, as well as the GPT - interaction and a phenomenological three-body force. Our results indicate a possible one-neutron halo in Na, while Na exhibit features of Borromean halos. The halo formation is driven by the weakening of the shell gap and inversion of the and orbits expected to occur somewhere near these masses. We further show that the electric dipole response provides a clear and sensitive probe of halo structure in these isotopes.

    nucl-thnucl-ex2 citations
  8. 08

    Basis truncation, statistical errors, and systematic uncertainties in relativistic approaches to nuclear response

    A.V. Afanasjev · E. Litvinova · B. Osei

    Although there exists a clear and, in principle, exact theoretical formulation for the equation of motion for the response of a correlated fermionic system, its numerical implementations for atomic nuclei require feasible approximations. One of the widely accepted approximations is a truncated harmonic oscillator (HO) basis, whose wave functions are used to expand the solutions obtained with realistic interactions. In this work, we extend previously employed HO basis truncated at = 20 fermionic shells to = 50 and perform a systematic study of the effects of such basis increase on nuclear resonances. The relativistic random phase approximation (RRPA) and its extension by the particle-vibration coupling dubbed as relativistic time-blocking approximation (RTBA) are applied to the description of the monopole, dipole, quadrupole, and octupole resonances in Ca, Ni, and Sn, and the RRPA studies are extended to Ca and Pb. A considerable sensitivity of the strength distributions to the HO basis size is found, especially for low-spin resonances in the light neutron-rich nuclei. The effects of the HO basis extension to = 50 are analyzed and linked to the involvement of proton and neutron continuum states and proton quasi-bound states in the strength formation. The obtained results point to the importance of the HO basis completeness and continuum effects in the nuclear response calculations and evaluation of the associated parameters of the nuclear equation of state. Statistical errors and systematic uncertainties in the RRPA strength functions are analyzed. They are found to be substantial for the monopole response, but significantly smaller for the dipole, quadrupole, and octupole ones. Neither of them shows a pronounced mass dependence, and statistical errors are generally smaller than systematic uncertainties.

    nucl-thPRC(2026)·1 citation
  9. 09

    Neutron Star Equation of State with Nucleon Short-Range Correlations: A Concise Review and Open Issues

    Bao-Jun Cai🇨🇳 · Bao-An Li🇺🇸 · Yu-Gang Ma🇨🇳

    Nucleon short-range correlations (SRCs) and the associated high-momentum tail (HMT) in its momentum distribution represent a universal feature of strongly interacting Fermi systems. In nuclear matter, SRCs arise primarily from the spin-isospin dependence of the tensor and short-range components of the nucleon-nucleon interaction, leading to a substantial depletion of its Fermi sea and a characteristic tail populated predominantly by isosinglet neutron-proton pairs. These microscopic structures modify both the kinetic and interaction contributions to the Equation of State (EOS) of dense matter and thereby influence a broad range of neutron-star (NS) properties. This short review provides a streamlined overview of how SRC-induced changes in reshape the kinetic EOS, including its symmetry energy part and how these effects propagate into macroscopic NS observables, including mass-radius relations, tidal deformabilities, direct Urca thresholds and core-crust transition. We summarize key existing results, highlight current observational constraints relevant for testing SRC-HMT effects, and outline open questions for future theoretical, experimental, and multimessenger studies of dense nucleonic matter.

    nucl-thastro-ph.HEhep-phhep-th+1Modern Physics Letters A, 2630005 (2026)·1 citation
  10. 10

    Poles from the conserved kinetic equation: The emerging gradient structure and causality riddle of relativistic hydrodynamics

    Sukanya Mitra🇮🇳

    In this work, the poles and the resulting dispersion spectra from the relativistic kinetic equation have been analyzed with the help of a proposed collision kernel that conserves both the energy-momentum tensor and particle current by construction. The dispersion relations, which originally come out in the form of logarithmic divergences, in the long wavelength limit exhibit the systematic gradient structure of the relativistic hydrodynamics. The key result is that, in the derivative expansion series, the spatial gradients appear in perfect unison with the temporal gradients in the non-local relaxation operator like forms. It is then shown that this dispersion structure, including non-local temporal derivatives, is essential for the preservation of causality of the theory truncated at any desired order.

    nucl-thgr-qchep-phhep-thPLB(2026)·0 citations
  11. 11

    Compact star and compact star matter properties from a baryonic extended linear sigma model with explicit chiral symmetry breaking

    Yao Ma🇨🇳 · Yong-Liang Ma🇨🇳 · Lu-Qi Zhang🇨🇳

    Based on a baryonic extended linear sigma model including explicit chiral symmetry breaking effect, the structure of neutron stars with the emergence of hyperons is investigated using the relativistic mean field approximation. It is found that, except for the lightest scalar meson whose structure is not well understood so far, the vacuum mass spectra of relevant hadrons and nuclear matter properties around saturation density can be well reproduced. Nevertheless, based on the present model and the applied relativistic mean field approach, we found that, to have a realistic mass-radius relation of neutron stars, the sigma term that denotes the contribution of explicit symmetry breaking should deviate from its empirical values at vacuum. Specifically, MeV, rather than at vacuum. With an appropriate choice of and , our framework can give a more observationally favored mass-radius relation of neutron stars with the emergence of hyperons, suggesting a possible density dependence of the low energy constants, at least within the present leading order framework with the relativistic mean field approach. The present result provides a new perspective on the relation between microscopic explicit chiral symmetry breaking in dense matter and macroscopic structure of compact stars and calls for more systematic treatments beyond leading order relativistic mean field calculation.

    nucl-thastro-ph.HEPRD(2026)·0 citations
  12. 12

    Recent advancements in the strongly coupled many-body theory for nuclear spectral computation

    Elena Litvinova

    Some recent advancements of the nuclear many-body theory and selected results on nuclear giant and pygmy resonances are presented. The theory is compactly reviewed, with a special focus on the emergent scale of the quasiparticle-vibration coupling (qPVC), which carries the order parameter associated with the qPVC vertex, and an efficient treatment of the nuclear many-body problem organized around the qPVC hierarchy. Self-consistent numerical solutions of the relativistic Bethe-Salpeter-Dyson equation for the nuclear response function in medium-heavy nuclei are discussed. The presented update on the pygmy dipole resonance focuses on establishing the formation of its two-component structure as a result of the fragmentation of the low-energy dipole mode due to the qPVC and its mixing with the similarly fragmented giant dipole resonance. The centroid of the isoscalar giant monopole resonance is also linked to qPVC effects, particularly to its sensitivity to the coupling of the collective breathing mode to the lowest quadrupole vibrations, which is enhanced by quadrupole collectivity. The resolution of the long-standing "fluffiness" puzzle regarding the compressibility of open-shell tin isotopes is summarized. The recently developed thermal variant of the superfluid response theory is briefly introduced and continuously linked to the description of the isoscalar monopole response at finite temperature with the prospect of refining the temperature-dependent nuclear equation of state.

    nucl-thEPJ Web Conf.(2026)·1 citation
  13. 13

    Anomalous Transport In Low Dimension Materials

    Leonardo Lopes

    This dissertation presents a systematic theoretical investigation into realizing a condensed matter analogue of the Chiral Magnetic Effect (CME) in a quasi-planar, 2+1D system. The research establishes a conceptual bridge between the anomalous transport phenomena of high-energy physics and the emergent electronic properties of engineered honeycomb lattices. The central objective is the formulation of a low-energy effective Hamiltonian that incorporates the necessary ingredients for a CME-like effect. This is achieved by moving beyond pristine graphene, whose inherent sublattice symmetry precludes the formation of a mass gap necessary for defining robust pseudo-chiral states. The core of this work is a model based on a honeycomb lattice with explicitly broken sublattice symmetry, which introduces a band gap and endows the quasi-particles system with a well-defined pseudo-chirality based on sublattice polarization. A time-reversal symmetry-breaking parameter is introduced to asymmetrically modify the valley gaps, creating a controllable non-equilibrium imbalance analogous to the chiral chemical potential in relativistic systems. A key finding is the validation of the physical model consistency; through commutator calculations, the total angular momentum - comprising both orbital and an emergent lattice spin component - is shown to be a conserved quantity. This research successfully transforms the abstract possibility of a 2D CME into a concrete, self-consistent theoretical framework, detailing the precise symmetry conditions required for its manifestation.

    cond-mat.othercond-mat.mes-hallhep-phhep-th+10 citations
  14. 14

    Orbital angular momentum in the pion and kaon: rest-frame and light-front

    Y.-Y. Xiao🇨🇳 · Z.-N. Xu🇪🇸 · Z.-Q. Yao🇩🇪 · C. D. Roberts🇨🇳 · J. Rodríguez-Quintero🇪🇸

    Orbital angular momentum (OAM) is not a Poincaré invariant quantity; so, its value is observer dependent. Notwithstanding that, in quantum chromodynamics, a Poincaré-invariant theory, OAM is part of every hadron wave function. Using continuum Schwinger function methods, we elucidate both the subjective character of in-hadron OAM and expose some of its impacts on pion and kaon structure and observables. For instance, working with light-front projections of their Bethe-Salpeter wave functions, it is found that the pion is a roughly 50/50 mix of light-front OAM zero and one components and the kaon is a 60/40 system. The overall picture is that (near) Nambu-Goldstone modes are complex bound states, each with significant intrinsic OAM, independent of the observer's reference frame. This feature must be accounted for in the calculation of observables. Inductively, the same is true for all hadrons.

    hep-phhep-exhep-latnucl-ex+1PLB(2026)·4 citations
  15. 15

    Nonequilibrium QCD in heavy-ion collisions: Kinetic theory and jet modifications during the initial stages

    Florian Lindenbauer🇦🇹

    This thesis focuses on how jets are modified by the nonequilibrium quark-gluon plasma during the initial stages in heavy-ion collisions. Its influence on their propagation is typically encoded in a single medium function, the dipole cross section. Its small distance behavior is characterized by the jet quenching parameter , and we obtain its numerical value throughout the pre-equilibrium stage, finding values comparable in magnitude to the earlier Glasma stage. We also compute the more general elastic collision kernel, obtained by Fourier transforming the dipole cross section. This constitutes an important step to facilitate the understanding of jet-medium interactions during the initial stages in heavy-ion collisions. Additionally, we improve QCD kinetic theory simulations by employing a more realistic (HTL) screening mechanism to incorporate medium effects, which we compare with simpler screening mechanisms. An expanding plasma realized in the initial stages of heavy-ion collisions exhibits a significantly reduced maximum anisotropy and reduced specific shear viscosity when using the improved screening prescription. Moreover, we investigate the gluon splitting rates, which are typically obtained using an isotropic model for the collision kernel. Going beyond that approximation, we find that the splitting rates obtained from the nonequilibrium anisotropic kernel differ significantly both in magnitude and in their qualitative time evolution. We further identify a novel type of weak-coupling attractor, which can be observed in the ratio of the jet quenching parameter and pressure ratio, and is obtained by extrapolating to vanishing coupling. This improved kinetic theory description and novel limiting attractors contribute towards a more realistic modeling of the nonequilibrium QCD plasma and its equilibration and hydrodynamization process during the initial stages.

    hep-phnucl-th0 citations
  16. 16

    Quark Coalescence: Formation of Mesons Including Excited States

    R. J. Fries · P. Virupapuram · J. Purcell · H. Anconetani · W. Lippincott · S. Robicheaux · M. Kordell · C. M. Ko

    We discuss the quantum mechanics of coalescence of quark-antiquark pairs into mesons using a non-relativistic quark model. We derive the coalescence probabilities assuming a harmonic oscillator potential and generic Gaussian wave packet shapes for the initial quarks and antiquarks. Our particular emphasis is on modeling excited states of the meson spectrum consistently. We provide the formalism to systematically include excited states from the Particle Data Book, and many more predicted by the quark model, up to and masses of about 2.2 GeV for light flavors. We provide estimates of masses and decay branching ratios for unconfirmed states. We use a phase space picture which is appropriate for the quasi-classical nature of the information typically available for the quarks and antiquarks in applications like Monte Carlo simulations. We demonstrate that for typical parton configurations expected in jets, excited meson states are populated abundantly.

    hep-phnucl-exnucl-th1 citation
  17. 17

    Fine-tuning final state interactions model in NuWro Monte Carlo event generator

    Hemant Prasad🇵🇱 · Jan T. Sobczyk🇵🇱 · Rwik Dharmapal Banerjee🇵🇱 · J. Luis Bonilla🇵🇱 · Krzysztof M. Graczyk🇵🇱 · Beata E. Kowal🇵🇱 · Artur M. Ankowski🇵🇱

    Recent experimental data from MINERvA on transverse kinematics observables across four different nuclear targets - carbon, oxygen, iron, and lead - have been utilized to refine the modeling of final state interaction effects in the NuWro Monte Carlo neutrino event generator. For this purpose, we have developed an event reweighting tool for future applications to adjust the strength of final-state interactions. This study highlights the requirement for stronger nucleon reinteractions than previously assumed, but it still falls within the uncertainty range observed in a study comparing proton transparency measurements. This conclusion has significant implications for both experimental and theoretical work involving NuWro.

    hep-phnucl-thPRD(2026)·1 citation
  18. 18

    Perturbative results for the matrix elements of the vector current and the role of different infrared regulators

    Alessio Carmelo Alvaro🇮🇹 · Ignacio Castelli🇺🇸 · Cédric Lorcé🇫🇷 · Andreas Metz🇺🇸 · Barbara Pasquini🇮🇹 · Simone Rodini🇮🇹

    We investigate the twist-2 unpolarized generalized parton distributions (GPDs) of quarks for an on-shell gluon target in quantum chromodynamics. These GPDs parametrize the leading-twist matrix elements of the nonlocal light-like flavor-singlet vector current. We compute them at one-loop accuracy in perturbation theory using a quark mass and dimensional regularization as infrared regulators. In particular, we discuss the limit of vanishing momentum transfer. The present work extends our previous related study on the axial current.

    hep-phnucl-thPLB(2026)·0 citations
  19. 19

    Limits on dark matter existence in neutron stars from recent astrophysical observations and mass correlation analysis

    Jing Fu Hu · Hang Lu · Bao Yuan Sun

    Dark matter admixed neutron stars (DANSs) serve as a specific astrophysical laboratory for probing the features of dark matter (DM) and have emerged as a promising candidate for interpreting recent astrophysical observations (e.g., by NICER and LIGO/Virgo). Accurately constraining the internal DM content of DANSs is therefore of critical importance. In this work, we construct the equations of state (EoS) for DANS matter by employing twelve nuclear matter (NM) models within the covariant density functional (CDF) theory and a self-interacting fermionic model for DM. Using these EoSs as input, we solve the two-fluid Tolman-Oppenheimer-Volkov (TOV) equations to systematically investigate the influence of DM on the global properties of neutron stars (NSs). By incorporating recent observational constraints on NS properties, the maximum DM mass fraction in DANSs is determined for each NM EoS model. Our analysis reveals a strong linear correlation (Pearson coefficient ) between and the maximum mass of a pure NS, , described by . Leveraging this correlation and the observed NS maximum mass distribution, , we derive the probability distribution function (PDF) for the maximum DM mass, , in DANSs. We find that at the 68\% confidence level, . This quantitative constraint on the DM mass provides a critical prior for interpreting potential observational signatures of DANSs, such as anomalous tidal deformabilities and distinctive gravitational-wave signals.

    astro-ph.HEnucl-thAstrophys.Space Sci.(2026)·1 citation
  20. 20

    The Effect of Hadronic Matter on Parton Energy Loss

    Ritoban Datta🇺🇸 · Abhijit Majumder🇺🇸

    Modified thermal distributions (dispersion relations) are introduced within both the MATTER and LBT event generators used to describe jet modification in a heavy-ion collision, within the JETSCAPE framework. Hard partons, propagating through dense matter, scatter off the partonic substructure of the medium, leading to stimulated emission, accompanied by recoiling medium partons. We introduce a simple modification, a multiplicative correction to the dispersion relation of quarks and gluons (equivalent to an effective fugacity). This leads to calculated transport coefficients (e.g. ) showing the expected behavior of depreciating at lower temperatures, including within the hot hadronic gas. This simple modification recovers the light-like dispersion relations at high temperatures, and introduces an excess depreciation factor for parton populations at lower temperatures, allowing partonic energy loss and recoil calculations to be extended into the hadronic phase. This modified distribution, in combination with initial state cold nuclear matter effects (shadowing), is used to simultaneously describe the nuclear modification factor and elliptic anisotropy of jets and leading hadrons, over multiple centralities and collision energies.

    hep-phhep-exnucl-exnucl-th2 citations

Affiliations

first authorsco-authorsvia INSPIRE