2026-06-05 | | Total: 13
We study charge-dependent modifications of identical-pion and identical-kaon femtoscopic correlation functions from two sources: the residual Coulomb field of the charged source and isospin-related hadronic dynamics. The residual Coulomb effect is modeled with a modified Retiere--Lisa blast-wave source, where the same emitted particles are propagated with positive and negative charge signs through an effective residual field. The residual Coulomb field produces a small but systematic positive-to-negative splitting of the correlation functions, strongest at low $k_T$ and sensitive to the effective charge, spatial width, and expansion velocity of the residual source. It also modifies the height and shape of the correlation function, complicating interpretation of the fitted radii. UrQMD 3.4 cascade calculations for Au+Au collisions at $\sqrt{s_{NN}}=7.7~\mathrm{GeV}$ show that charge-dependent splittings can also appear even without a residual Coulomb field. For pions, this difference is mainly driven by the initial isospin composition, while for kaons it is strongly affected by different $K^+$ and $K^-$ production mechanisms and subsequent hadronic evolution. These results show that residual Coulomb and isospin-related effects can compete, and neither can be interpreted reliably without constraining the other within the same model framework.
Whether cold dense QCD matter undergoes a strong first-order phase transition remains an open question. In nature, neutron stars provide the most direct probe of cold dense QCD matter. Theoretically, chiral effective field theory constrains the equation of state of dense matter near nuclear saturation density, while perturbative QCD calculations constrain it at densities well beyond stable neutron-star interiors. We perform Bayesian inference with non-parametric Gaussian-process equation of state for $β$-equilibrated neutron-star matter under the assumption with and without a strong first-order phase transition, using the tidal deformability from GW170817, the NICER mass--radius measurements of PSR~J0740$+$6620, PSR~J0030$+$0451, PSR~J0437$-$4715, PSR~J0614$-$3329, chiral effective field theory, and perturbative QCD. Our results favor a strong first-order phase transition, with its onset most likely lying \emph{above} the central density of the most massive neutron star. Such an onset reconciles the stiffness required to support massive neutron stars with the softening favored by perturbative QCD from asymptotically high density.
A consistent transformation from momentum space to energy space is performed for the diffusion equation within kinetic theory, with particle-number conservation explicitly preserved. Under the assumption of a constant single-particle level density, the equation reduces to a one-dimensional diffusion equation in energy space with consistent kinetic coefficients. The equivalence of the definitions of the equilibrium temperature in momentum space and in energy space is demonstrated. It is established that the inclusion of the energy dependence of the kinetic coefficients leads to an energy-dependent equilibrium temperature and a modification of the distribution function. The obtained results may be used to analyze relaxation processes in atomic nuclei and nonequilibrium dynamics of Fermi systems.
The present Reply addresses the Comment as posted on arXiv (arXiv:2606.04137 [nucl-th], June 2026).
We investigate physics-informed neural networks (PINNs) for computing the $s$-wave bound states of the hypernucleus $_Λ^{208}$Pb, modeled as a two-body system composed of a $Λ$ hyperon and a $^{207}$Pb core. The interaction is described by a Woods--Saxon potential without spin--orbit coupling. In the PINN formulation, the radial bound-state wave function is represented by an artificial neural network, and the eigenenergies are obtained from the Rayleigh--Ritz variational quotient. Because PINN eigenvalue calculations can depend on the residual-loss formulation and random-seed initialization, two residual losses are compared across four independent random seeds. Their performance is assessed using eigenvalue accuracy, coefficient of variation, signal-to-noise ratio, bias--variance decomposition, and Hermitian spectral-ordering consistency. The normalized residual loss gives the most stable and physically consistent bound-state spectrum for the four seeds considered. With this loss, the computed bound-state energies and root-mean-square radii are in very good agreement with the corresponding theoretical values. The resulting wave functions are used to construct the Jost functions through the Wronskian with incoming and outgoing Riccati--Hankel functions. From these Jost functions, the residue of the partial-wave $S$-matrix at the bound-state pole, the Asymptotic Normalization and the Nuclear Vertex Constants are extracted. These quantities show reasonable agreement with theoretical values, although with larger standard deviations than those obtained for the eigenenergies and radii. The results indicate that the Rayleigh--Ritz formulation combined with a seed-robust normalized residual loss provides a stable PINN framework for bound-state and Jost-function-based calculations in hypernuclear two-body systems.
In this work, we discuss the $γW$-exchange contributions in neutron $β$ decay with an elastic intermediate state, beyond the forward-angle limit (FAL). By decomposing the one-$W$-exchange and $γW$-exchange amplitudes in terms of $16$ independent Pauli-spinor structures, we calculate the $γW$-exchange corrections to the relevant coefficients. Our numerical results show that the relative corrections to the Fermi Born term $C_{\text{Born}}^{\text{F}}$ and the Gamow-Teller Born term $C_{\text{Born}}^{\text{GT}}$ are enhanced by about 8\% and 18\%, respectively. In particular, we find a non-zero contribution to $C_{\text{Born}}^{\text{GT}}$ from the axial-vector current, which is identically zero in the FAL. The corresponding effect on the extracted $V_{ud}$ from the neutron lifetime is also analyzed, and we find the correction to be at the $10^{-4}$ level.
We derive dissipative relativistic hydrodynamic equations in the presence of multiple conserved charges, i.e., baryon number ($B$), electric charge ($Q$), and strangeness ($S$), using the Chapman-Enskog (CE) method within the kinetic theory approach. The relativistic Boltzmann equation is solved within the relaxation-time approximation with a momentum-independent relaxation time in the collision term. We derive both first-order (Navier-Stokes limit) and second-order dissipative hydrodynamic equations. Within the kinetic theory framework, using the Boltzmann's H-theorem, and by demanding that for a dissipative system, the entropy must be produced, we find different transport coefficients at the first-order and second-order gradient expansion of the out-of-equilibrium distribution function around the local equilibrium. Apart from the well-known transport coefficients, the shear ($η$) and the bulk ($ζ$) viscosities , we also find the diffusion matrix elements ($κ_{qq^{\prime}}$) for the conserved charges $B$, $Q$ and $S$. The diffusion matrix elements ($κ_{qq^{\prime}}$) are important to model the multi-component diffusion dynamics sourced by inhomogeneous baryon stopping in the initial state of heavy-ion collisions. We estimate the temperature ($T$) and chemical potential dependence of diagonal and off-diagonal elements of the diffusion matrix elements for the (2+1) flavor quark-gluon plasma. We further estimate the ratio $κ_{qq^{\prime}}T/η$ for a wide range of temperature and chemical potentials to show the relative importance of the diffusion matrix elements compared to other transport coefficients.
The spin-orbit correlation (SOC) provides a unique probe into the internal spin structure of hadrons. Defined via the parity-odd (P-odd) energy-momentum tensor (EMT), this observable can remain non-vanishing even in systems where the total angular momentum is zero. In this study, we connect the formal field-theoretical definition of the SOC to a non-perturbative quantum many-body framework utilizing light-front dynamics. Furthermore, we show how the SOC can be extracted directly from the hadronic matrix elements of the P-odd EMT, establishing a pathway to access this observable within the partonic picture. As a practical application, we compute the transverse and longitudinal SOC distributions for charmonium and $B_c$ mesons. While our findings align with rough estimates based on the Clebsch-Gordan decomposition, we demonstrate that these observables yield rich, non-trivial information regarding partonic dynamics.
We report the first observation of a large directed flow ($v_1$) of $φ$ mesons in mid-central Au+Au collisions at $\sqrt{s_{\mathrm{NN}}}=3.0$-4.5 GeV, together with the $\sqrt{s_{\mathrm{NN}}}=7.7$ GeV collider-energy result for the collision-energy dependence, measured by the STAR experiment at RHIC. Despite its mesonic nature, the $φ$-meson $v_{1}$ exhibits a magnitude comparable to that of protons and $Λ$ baryons, and is significantly larger (by a factor of $\gtrsim3$) than that of $K^0_S$ mesons at the lowest energies. The midrapidity slope $dv_1/dy$ shows a pronounced energy dependence similar to that observed for protons and $Λ$ baryons, in contrast to lighter mesons, in the high-baryon-density region. Comparisons with hadronic transport calculations indicate that the observed $φ$-meson directed flow is consistent with production via high-mass baryon resonances and the collective motion of baryons influenced by baryonic mean-field interactions. Due to its long lifetime and relatively weak effective $φ$-$N$ interaction, the $φ$ meson can retain sensitivity to collective dynamics established during the dense baryonic stage of the collision. These results provide new experimental constraints on baryon-resonance production and strange-hadron transport dynamics in dense QCD matter.
The proposal to formulate chiral gauge theories using domain wall fermions on $2n+1$ dimensional Euclidean lattice with a slab geometry involves $2n$ dimensional dynamical gauge fields residing on one of the domain walls. The gauge fields are extended into the extra dimension using gradient flow decoupling the mirror fermions on the anti-wall. We implement this construction on the lattice for $n=1$ in the presence of $2n$ dimensional background gauge fields. We also formulate and implement an additional gauge field flow proposal, where the gauge fields satisfy $2n+1$ dimensional equation of motion away from the domain wall, known as the EOM (equation of motion) flow. In both cases, we couple the gauge fields to fermions and demonstrate how current conservation and anomaly inflow work on the lattice.
A recent non-perturbative formulation of $2n$ dimensional chiral gauge theories relies on realizing chiral fermions on the $2n$ dimensional boundary of a $2n+1$ dimensional disk manifold. It also requires extending boundary gauge configurations into the interior of the disk using some flow prescription that preserves 2n dimensional gauge invariance. In this paper we propose a concrete realization of the equation of motion flow with the disk embedded on a square lattice. In addition, we couple the flow gauge field to fermions and demonstrate the mechanism of anomaly inflow and anmaly cancellation at work on the lattice.
We present a continuous-variable adaptive variational quantum eigensolver (CV-ADAPT-VQE). As concrete examples, we consider the ground-state preparation for (i) the Bose-Hubbard model and (ii) the bosonic Kitaev chain, including its extension with an on-site Kerr interaction. The former conserves the total boson number, while the latter conserves global parity. We construct symmetry-preserving operator pools tailored to each case and show, using GPU-based classical simulations, that CV-ADAPT-VQE results in significantly shallower circuits compared to Hamiltonian-based VQE approaches. Our results point toward direct applications in quantum simulations of condensed-matter systems, quantum chemistry, and high-energy physics.
We assess the impact of future measurements of dihadron production in semi-inclusive deep-inelastic scattering from the CLAS12 and proposed SoLID experiments at Jefferson Lab, as well as from the ePIC experiment at the future Electron-Ion Collider (EIC), on the transversity parton distribution functions (PDFs) and the corresponding tensor charges of the nucleon. To this end, we generate pseudo-data for these experiments for a proton target (CLAS12 and ePIC) and a $^3$He target (SoLID and ePIC), and we include these pseudo-data in the JAMDiFF global analysis of existing experimental dihadron data. We find that future data from Jefferson Lab will significantly reduce uncertainties in the transversity PDFs in the region of intermediate-to-large quark momentum fractions $x$, while the EIC will provide strong constraints across the entire range of $x$, allowing for the first experimental test of the predicted small-$x$ behavior of the transversity PDFs. In discussing the reduction of uncertainties in the tensor charges, we also compare the results from the data analyses with those from lattice QCD, highlighting scenarios in which compatibility or tension between the two would arise.