Papers

We investigate the specific heat and isentropic bulk modulus of finite-temperature pure SU(3) gauge matter within a lattice-data-driven phenomenological framework. The equation of state is formulated in terms of a temperature-dependent effective gluon mass constrained { by lattice QCD pressure data as input, allowing the pressure}, trace anomaly, gluon number density, energy per thermally active gluonic mode, and derivative-sensitive response functions to be derived in a thermodynamically consistent manner. The resulting pressure and trace anomaly reproduce the characteristic lattice behavior across the deconfinement region, while the effective gluonic degrees of freedom increase rapidly above $T_c$. The normalized specific heat $C_V/T^3$ develops a pronounced enhancement in the vicinity of $T_c$, reflecting the rapid temperature variation of the energy density across the deconfinement region. The isentropic bulk modulus $K_S/T^4$ also rises sharply across the transition region, indicating a substantial stiffening of the equation of state. At high temperatures, both response functions gradually approach values close to their massless conformal Stefan–Boltzmann reference values, with $\left(C_V/T^3\right)_{\rm SB}=32\pi^2/15\simeq 21.06$ and $\left(K_S/T^4\right)_{\rm SB}=32\pi^2/135\simeq 2.34$. These findings indicate that the specific heat and isentropic bulk modulus provide complementary constraints on the temperature evolution of nonconformal dynamics in pure SU(3) gauge matter.

Though quantum computing naturally offers an advantage for simulations of real-time quantum systems, implementing Imaginary-Time Evolution (ITE) is comparatively more difficult. Nevertheless, quantum implementations of ITE are useful both in cases where classical implementations have associated sign problems, and also as an exact method of preparing eigenstates on quantum computers. In this work we present an algorithm to obtain ITE of a generic quantum Hamiltonian by performing analytic continuation of measured real-time correlation functions. We present simulations and demonstrations on IBM quantum hardware of 1D Fokker-Planck equations for classical diffusion process and imaginary-time evolution of integrated correlation functions in 1D quantum mechanical scattering as examples to demonstrate the effectiveness of the method.

We investigate the skewness dependence of the pion generalized parton distribution (GPD) at different momentum transfers $-t$ within the covariant Nambu-Jona-Lasinio (NJL) model, employing the Schwinger proper-time regularization scheme to regulate ultraviolet divergences and implement an effective description of quark confinement. We evolve the pion GPDs for various values of $-t$ and $\xi$ to the scales $\mu^2=4$ and $27~\mathrm{GeV}^2$. We find that the valence-quark distributions are suppressed with increasing $-t$ and $\xi$, with the suppression becoming more pronounced at the lower scale, $\mu^2=4~\mathrm{GeV}^2$. In the forward limit, $\xi=0$ and $-t=0$, the predicted pion valence-quark distribution at $\mu^2=27~\mathrm{GeV}^2$ is in good agreement with the available experimental data and the global JAM analysis. The corresponding pion gluon distribution at $\mu^2=4~\mathrm{GeV}^2$, obtained through next-to-leading-order (NLO) DGLAP evolution, is also consistent with the JAM analysis. We further find that the pion and kaon valence-quark distributions exhibit a strong dependence on $-t$ but only a weak dependence on $\xi$, whereas the corresponding gluon distributions show weak dependence on both $-t$ and $\xi$. The generalized form factors of the pion and kaon at $\mu^2=4~\mathrm{GeV}^2$ are in good agreement with the corresponding lattice-QCD results.

This memorial tribute to Mannque Rho follows a line of thoughts and ideas that he continuously inspired and shaped over many decades: from the two-scales picture of low-energy nucleon structure to dense and cold baryonic matter as it is realized in the cores of neutron stars. Early groundbreaking concepts are recalled and updated by recent advanced analyses of the `core' and `cloud' sizes of the nucleon. Implications for dense nuclear matter are then discussed and confronted with empirical information from Bayesian inference analyses of neutron star observables.

We show that color triality constrains the nonperturbative states that screen a Wilson-line endpoint, allowing the joint flows of net electric charge, baryon number, and strangeness to probe QCD vacuum hadronization. Whether evaluated from resolved hadrons in a jet initiated by a quark of flavor $f$ or from the corresponding charge correlators, these flows satisfy the Gell-Mann–Nishijima relation $\langle Q_{\rm flow}\rangle_f \simeq I_{3,f} +\frac{1}{2}\left( \langle S_{\rm flow}\rangle_f +\langle B_{\rm flow}\rangle_f \right)$, where $I_{3,f}$ is the third component of the initiating-quark isospin. The net jet baryon number, $\langle B_{\rm flow}\rangle_f\simeq \frac{r_{qq}}{1+r_{qq}}$, directly probes the relative probability of a diquark–antidiquark vacuum excitation, with $r_{qq}$ the diquark-to-quark production ratio. Thus, for $r_{qq}\ll1$, the baryon number carried by the jet is substantially suppressed relative to the initiating-quark value $B_f=1/3$. Likewise, the strange-to-light pair-production ratio, defined by $u\bar u:d\bar d:s\bar s=1:1:r_s$, is encoded in the measured jet strangeness $\langle S_{\rm flow}\rangle_f \simeq S_f+\frac{3r_s}{2+r_s} \left(\frac{1}{3}- \langle B_{\rm flow}\rangle_f\right)$, predicting a nonzero mean net strangeness even in $u$- and $d$-initiated jets. The conserved-current moments appearing in these relations are independent of the renormalization scale. Their simultaneous measurement therefore provides a direct, flavor-resolved probe of vacuum pair production and quantum-number transport during hadronization.

One-proton emission is a direct probe of nuclear structure near the proton drip line and plays a critical role in understanding exotic decay modes and nucleosynthesis processes. In this study, we investigate the half-lives of one-proton emitters for $71 \leq Z \leq 83$ odd-$Z$ nuclei by employing the WKB approximation with nuclear potentials obtained from the deformed relativistic Hartree-Bogoliubov theory in continuum (DRHBc) and, for comparison, the relativistic continuum Hartree-Bogoliubov theory (RCHB). We first compare the calculated half-lives with available experimental data. The inclusion of quadrupole deformation via the DRHBc hardly contributes to improving the predictions of half-lives for the deformed nuclei. We find that all the studied nuclei exhibit ground states with $|\beta_{2,{\rm DRHBc}}| < 0.15$, and within this limited deformation range the spectroscopic factor provides the dominant contribution to the half-life, compared to the decay width. In particular, for nuclei exhibiting shape coexistence in DRHBc, such as $^{170}$Au, where the half-life varies significantly with the quadrupole deformation through its effect on the spectroscopic factor, we expect shape coexistence to exert a substantial influence on the variation of half-lives. Finally, we discuss the half-lives in consideration of shape coexistence. Our results indicate that the calculated half-life is governed not by the total-energy difference between coexisting minima but rather by the spectroscopic factor influenced by the deformation.

We present the discovery and analysis of three rotational glitches in the young pulsar PSR J1637$-$4642. The timing observations span from 19 February 2009 to 6 October 2024 (MJD 54881$-$60589) from the Murriyang radio telescope of the Parkes Observatory. The first and strongest glitch occurred around MJD 58352 with a fractional frequency change of $\Delta\nu/\nu \sim 2.7 \times 10^{-6}$, while two additional smaller glitches were detected at MJD 59443 and MJD 60445 with fractional changes of $2.2 \times 10^{-9}$ and $2.8 \times 10^{-8}$, respectively. Prior to this, the pulsar had shown no glitch activity since its discovery in the Parkes Multibeam survey. Only the first glitch exhibits detectable exponential recovery, with a decay timescale of $\sim$100 days and a small recovery fraction $\approx 0.015$, accompanied by a permanent increase in the magnitude of the spin-down rate. Modeling the post-glitch evolution of $\dot{\nu}$ within the vortex-creep framework using Bayesian inference gives a superfluid moment-of-inertia fraction $\approx 0.0187$, consistent with the inner-crust superfluid. These results reinforce the standard superfluid glitch paradigm and demonstrate that even “quiet” pulsars can still host substantial glitch activity.

Hydrodynamic universality refers to the property that different formulations of relativistic dissipative hydrodynamics yield identical predictions in the asymptotic long-wavelength regime. We show that neutron-star oscillations need not reach this regime. Because the finite stellar radius limits the accessible wavelengths and causality imposes a lower bound on microscopic relaxation times, the separation between microscopic and macroscopic scales can become insufficient for hydrodynamic universality to emerge. Comparing linear oscillations in relativistic Navier–Stokes and Israel–Stewart theories, we find that realistic bulk viscosities can produce sizeable modifications of the oscillation spectrum, even at the longest wavelengths. These results identify neutron-star oscillations as a direct probe of microscopic nonequilibrium dynamics beyond the leading hydrodynamic description.

We present the first application of the unbiased exponential resummation method to the determination of the QCD pseudocritical temperature $T_{\rm pc}$ at vanishing baryon chemical potential $\mu_B$. By reconstructing the finite-density baryon number susceptibility $\chi_2^B$, we define two thermal-derivative observables whose peak positions yield $T_{\rm pc}^T=160.2(4)(3)$ MeV and $T_{\rm pc}^C=158.1(5)(3)$ MeV. We also show that the pseudocritical temperature, $T_{\rm pc}^\chi$ obtained from the peak of the fourth-order baryon number susceptibility $\chi_4^B$ is consistent with these determinations and with previous lattice-QCD results based on Taylor expansions \cite{Bollweg:2022Pade}. Further, we demonstrate that scaling relations \cite{Bollweg:2022fqq} between $\chi_4^B$ and the thermal derivatives of $\chi_2^B$ yield mutually consistent estimates of the leading second order curvature coefficient $\kappa_2^B$ at the corresponding pseudocritical temperatures. A direct analysis of the $\mu_B$ dependence of the pseudocritical temperatures provides an independent but substantially less constrained determination of the curvature, which remains statistically consistent with the scaling-based estimates. These results demonstrate the consistency of unbiased exponential resummation with the expected scaling behavior and establish it as a complementary approach for probing the QCD crossover at small finite baryon density.

The dynamics of heavy quarks in quark-gluon plasma (QGP) formed in heavy ion collisions provide a unique window to characterize its properties. Existing approaches to describe heavy quarks in medium rely either on quasiparticle-based models of QGP, or on assuming that the momentum transfer from the medium follows Gaussian statistics. However, neither of these assumptions can be taken for granted in QCD. In fact, in the prototypical theory that is strongly coupled $\mathcal{N}=4$ SYM, it has been long known there are no momentum-carrying quasiparticles, and furthermore, we have recently shown that the momentum transfer from the medium is far from being Gaussian [arXiv:2501.06289]. Since then, we showed that the asymmetry of said momentum transfer between energy loss and energy gain – which affects all moments of the distribution, not only its Gaussian characteristics – is, in fact, universal [arXiv:2504.21139]. Therefore, in order to connect the initial heavy quark production cross section with the final hadron spectrum in a way that is consistent with QFT principles, new methods are needed. In this talk, we present a new transport description of heavy quarks that encodes all of the non-Gaussian features of the momentum transfer from the medium, all of which can be defined and in principle calculated in QCD, without relying on any assumptions regarding the strength of the coupling. As a demonstrative example, we discuss heavy quark equilibration in $\mathcal{N}=4$ SYM. This paves the way towards extracting novel information about fundamental properties of QGP.

Quantum magic, or nonstabilizerness, is a quantum resource associated with computational advantage in quantum systems. In high energy collisions, Quantum Electrodynamics (QED) is inefficient at generating magic while the weak mixing angle, a fundamental constant of nature, sits near a value that minimizes magic production in charged-lepton scattering. These observations were made in the laboratory (lab) basis, in which spin is projected along the incoming beam axis. An alternative choice is the helicity basis, in which spin is projected along the direction of motion of each particle. The transformation between these two bases is, in general, not a Clifford operation and therefore can change the amount of magic. We present a detailed study of magic production in both bases for QED and electroweak processes, and compare these results with the basis-invariant non-local magic. In the ultra-relativistic limit, magic production is generally smaller in the helicity basis due to helicity selection rules, while the lab basis generally yields less magic in the non-relativistic regime. We provide circuit realizations of the ultra-relativistic Bhabha amplitudes using linear combinations of unitaries and show that the lab basis construction contains a larger $T$-gate count at generic scattering angles. Interestingly, in both bases the physical weak mixing angle lies close to the value that minimizes magic production.

We investigate uncertainty quantification for the neutron-star equation of state (EOS) by comparing point-estimation and distributional inference approaches using the same Chebyshev and piecewise-linear parameterizations. We combine neutron-star mass–radius and gravitational-wave tidal-deformability information within Bayesian, multilayer-perceptron (MLP), and normalizing-flow frameworks. Although the methods yield similar mean EOS behavior, the deterministic MLP produces substantially narrower uncertainty bands at high densities. We show that this behavior is associated with the point-estimation objective, which maps degenerate solutions toward the conditional mean rather than representing the full parameter posterior. By contrast, the normalizing flow yields distributions more consistent with the Bayesian inference. Our results demonstrate that reliable uncertainty quantification of the high-density EOS requires methods that represent conditional probability distributions rather than only point estimates.

QCD, the theory of the strong interactions, is formulated in terms of quarks and gluons, while low-energy nuclear physics deals with hadrons such as protons, neutrons, and pions. Symmetries establish a systematic connection between these two descriptions of strongly-interacting systems. The objective of this article is to review the symmetries of QCD and to explain how they constrain hadronic interactions. Chiral symmetry, which emerges in QCD in the limit of massless quarks, is of particular importance for low-energy nuclear physics. Together with its explicit and spontaneous breaking, chiral symmetry provides the basis for chiral perturbation theory, the effective field theory describing pions and nucleons at low energies.

First-principles calculations of atomic nuclei are necessarily incomplete as the Schrödinger equation is solved using approximate methods and due to the finite dimension of the employed Hilbert space. By balancing many-body truncation and basis-size uncertainties, we formalize a criterion for the optimal one-body basis dimension in a given ab initio nuclear structure computation. Next, it is demonstrated that higher-order many-body contributions can be computed using smaller basis sizes than used for the lower orders when a consistent accuracy in the calculation is targeted. Our findings are empirically validated using many-body perturbation theory and coupled-cluster calculations of nuclei spanning a large portion of the nuclear chart using two sets of chiral two- and three-nucleon interactions. The results suggest that considerable computational savings can be obtained using many-body-order-dependent one-body basis sizes.

We review the connection between constituent-quark Hamiltonians and QCD and investigate the long-standing difficulty of describing meson and baryon spectra with one common two-body interaction. A Hamiltonian calibrated to ground-state mesons leaves systematic baryon mass residuals, largest in the light-quark sector and decreasing toward heavier flavors. We show that a short-range, color-spin-dependent connected three-quark interaction substantially reduces this incompatibility. Mass-scaled finite-range profiles yield high-accuracy baryon spectra, whereas flavor-independent common-range profiles do not remove the residual flavor pattern. The result is tested on additional ground-state baryons outside the calibration set and through meson–baryon compatibility analyses across several alternative quark-model Hamiltonians. We also benchmark radial and orbital excitations to identify the regime in which a static compact valence Hamiltonian remains reliable, and provide explicit color-spin matrix elements for two- and three-body operators in baryons and multiquark configurations. Within the tested valence-space representations, the results indicate that a mass-dependent short-range connected three-quark interaction provides the missing contribution required for a consistent simultaneous description of meson and baryon ground-state spectra.

Heavy quarkonium production inside jets offers a sensitive probe of QCD dynamics and bound-state formation mechanisms. While recent studies demonstrate that charmonium-in-jet observables effectively discriminate among competing nonrelativistic QCD (NRQCD) long-distance matrix element (LDME) sets, whether this discriminating power persists in the bottomonium sector remains an open question. Here, we present the first phenomenological study of $\Upsilon(1S)$, $\Upsilon(2S)$, and $\Upsilon(3S)$ production inside jets using the fragmenting jet function (FJF) framework at next-to-leading order (NLO), incorporating DGLAP evolution, threshold resummation, and feeddown contributions from higher bottomonium states. In sharp contrast to charmonium, we find that bottomonium-in-jet momentum-fraction ($z_H$) distributions exhibit a universal shape that is remarkably insensitive to the choice of LDME sets. We show that this universality stems from the strong dominance of the S-wave spin-triplet color-octet ($^3S_1^{[8]}$) production mechanism reinforced by $\chi_b$ feeddown transitions. Our predictions capture both the characteristic large-$z_H$ peak and the spectral broadening with increasing jet transverse momentum observed in recent CMS measurements. These results establish a clear physical distinction between charmonium and bottomonium fragmentation inside jets, providing a theoretical benchmark for future high-precision measurements at the LHC.

The antineutron$\unicode{x2013}$nucleus scattering length is currently known only indirectly, via antiproton$\unicode{x2013}$nucleus optical potentials fitted to level shifts and widths of antiprotonic atoms. The antineutron$\unicode{x2013}$nucleus and antiproton$\unicode{x2013}$nucleus potentials are related to each other through charge symmetry. We calculate the scattering length from optical potentials proposed for antiprotonic atoms using nucleon density distributions as input. We find that the scattering length for an $N>Z$ nuclide is largely affected by the poorly constrained neutron density distribution and by a possible isovector interaction, one of the mechanisms introduced to reproduce the isotope dependence of antiprotonic $^{40,48}\mathrm{Ca}$ data. For $^{48}\mathrm{Ca}$, the isovector term modifies the scattering length by $0.3\unicode{x2013}0.4\,\mathrm{fm}$ ($0.4\unicode{x2013}0.5\,\mathrm{fm}$) for the real (imaginary) part, an order of magnitude beyond the uncertainty propagated from the isoscalar potential. As no antineutron$\unicode{x2013}$nucleus scattering data are available below $76\,\mathrm{MeV}/c$, a direct measurement with recently proposed low-energy antineutron beams would provide the first access to the antinucleon$\unicode{x2013}$nucleus interaction in the $s$-wave regime.

Given a certain specific, $M$ by $M$ matrix form of the Feshbach's effective (i.e., energy-dependent) Hamiltonian $H^{(M)}_{e\!f\!f}(E)$, the inverse-problem reconstruction of the full-space, $N$ by $N$ matrix Hamiltonian $H^{(N)}$ is considered and reduced to the solution of a coupled set of polynomial algebraic equations. Using computer-assisted symbolic manipulations, an explicit algebraic reconstruction of $H^{(N)}$ is found feasible at not too large $K = N-M$.

Calculating cross sections in the Color Glass Condensate effective theory requires Fourier transforms of multi-Wilson-line correlators from transverse coordinate space to transverse momentum space. Under the common assumption of impact-parameter independence, each transform reduces to a set of Hankel transforms whose Bessel-function kernels oscillate rapidly at phenomenologically relevant momenta, making direct quadrature prohibitively expensive. We present a Filon-type quadrature, applicable to any integrand, that integrates these oscillatory factors in closed form on the stored coordinate grid, reducing each Hankel transform to a precomputed weight vector and the full nested transform chain to a sequence of matrix products. We develop and validate the method on the deep inelastic scattering dijet cross section beyond the correlation-limit approximation, where an exprel-based reformulation of the quadrupole Wilson-line correlator removes a numerical $0/0$ instability inherent to its standard parametrization. Porting the calculation to the Graphics Processing Unit (GPU), with custom CUDA kernels that fuse the momentum-space contraction directly into the correlator evaluation, brings the runtime for one dipole input down to about two minutes on a single NVIDIA A800, from several hours on a multi-core Central Processing Unit (CPU). We further generalize the algorithm to three sequential Hankel transforms and validate the resulting six-dimensional transform against an analytic Gaussian integrand family with closed-form results at every stage. This general, process-independent algorithm is directly applicable to next-to-leading-order proton-nucleus and electron-ion scattering cross-section calculations performed without the correlation-limit approximation. The code is publicly available at https://github.com/CCNU-CGC-py/FFT_filon.

Reaction cross sections $\sigma_{\rm R}$ are widely used to infer matter radii, yet their sensitivity to nuclear structure beyond radial one-body distributions is less well understood. We combine complete $A$-body nucleon configurations sampled from ab initio nuclear lattice effective field theory (NLEFT) with event-by-event Monte Carlo Glauber calculations, thereby retaining the many-body correlations encoded in NLEFT. Using a fixed binary-collision prescription determined by the measured energy- and isospin-dependent total nucleon-nucleon cross sections, the calculations capture the overall magnitudes and energy dependence simultaneously for the available data on $^{12}$C and $^{9}$Be projectiles on carbon and hydrogen. Controlled randomization of angular correlations at fixed matter root-mean-square radius and spherically averaged one-body radial density produces only a weak change in $\sigma_{\rm R}$ for $^{12}$C but approximately a $10\%$ increase for $^{9}\mathrm{Be}+{}^{1}\mathrm{H}$. The calculations also capture the measured rise–plateau–sharp-rise–reduction trend across $^{7,9\text{--}12}$Be, a distinctive pattern reflecting the evolution of cluster and halo structures along the isotopic chain. These results show that $\sigma_{\rm R}$ retains sensitivity to intrinsic many-body geometry beyond a single inferred matter radius, opening a route to studies of exotic $\alpha$-cluster geometries and spatial nucleon correlations through reaction cross sections.

The performance of five model averaging methods, including the arithmetic mean (AM), weighted mean (WM), naive Bayesian model averaging (NBMA), principal component analysis (PCA), and power-moderated mean (PMM) methods, in nuclear charge radius predictions is investigated. Five commonly used nuclear charge radius models are adopted as inputs for the averaging procedures. The charge radius differences between the experimental data and the original nuclear models are analyzed and the results after considering the model averaging methods are also discussed. The calculations show that the NBMA method can provide the best root-mean-square (rms) deviation among these five model averaging methods. The PCA method can extract useful physical information and not only helps to interpret the model differences but also offers a feasible way to construct improved empirical models by recombining the principal components. In contrast to the other methods, whose results worsen upon including a new model with a larger rms deviation, the rms deviation of the PCA method remains almost unaffected. The PMM method is capable of integrating the strengths of various nuclear models and delivering reasonable uncertainty estimates not only in known regions but also in unknown ones. This method can automatically adjust data uncertainties to achieve consistency, and it can provide a tool for a smooth transition of the nuclear charge radius prediction from the WM to the AM. The extrapolation ability of these model averaging methods is checked by 66 newly observed data after year 2021. The calculations show that model averaging offers a reliable strategy for nuclear charge radius predictions, combining high accuracy on known data with robust extrapolation to new measurements. The charge radii and the odd-even staggering in calcium isotopes are also discussed.

Precision experiments of atoms and molecules have become a powerful probe of the electroweak structure of atomic nuclei and of physics beyond the Standard Model. We review how the interaction between a nucleus and its surrounding bound electrons can be exploited to precisely measure the electromagnetic, parity-violating, and CP-violating properties of nuclei and their fundamental constituents. We focus on rare, unstable isotopes, surveying the experimental techniques and facilities developed in recent years that have extended these measurements to the most exotic regions of the nuclear chart. Recent advances in the precision control and interrogation of single molecules, together with direct laser excitation of nuclear transitions, are opening new frontiers in nuclear and particle physics. At the same time, progress in nuclear theory, machine learning, and high-performance computing is strengthening the connection between our microscopic description of nature and laboratory observables. In many cases, the precision with which nuclear and particle physics properties can be extracted is now limited not by experiment, but by the molecular, atomic, or nuclear theory required to interpret the measurements. This challenge presents a major opportunity for combined theoretical and experimental advances that will enable future discoveries.

We define detector operators in the nonrelativistic conformal field theory describing fermions at unitarity. We reduce the problem of computing the momentum distribution and correlation between final particles produced by a local source ("point-produced") to the computation of correlation functions involving the detector operators. The general formalism is applied to the point production of three unitary fermions, where we find the momentum and angular distribution of final particles. We discuss a nonrelativistic version of celestial holography, which maps the asymptotic out-state to a quantum wave function in the so-called "oscillator frame."

We present a density functional-based equation of state for warm, dense nuclear matter with a transition to deconfined quark matter for applications to simulations of supernova explosions and neutron star mergers, but also for the cosmological evolution of Q-balls. For the quark matter equation of state, we employ a recently developed confining density functional approach while nuclear matter is described within a relativistic density functional model of the DD2 class. The phase transition is obtained by a Maxwell construction at constant entropy per baryon. We discuss the solutions of TOV equations for isentropic hybrid stars for two hybrid equation of state models: DDf-SFM without and DD2-$\chi$CDF with color superconductivity. We find that at finite temperatures above a critical value of entropy per baryon sequences of disconnected third family branches ("thermal twin stars") may appear for the DDf-SFM model, while they are absent for the color superconducting model and at $T=0$. We discuss the relation of this critical entropy per baryon to the Seidov criterion of gravitational instability for $T=0$ and find that it is a good guide. We suggest that the presence of thermal twin stars may be regarded as an indicator for the core-collapse supernova explodability of massive blue supergiant stars with solar metallicity and thus serve as a new criterion for the reliability of hybrid equation of state models. By this argument, strong color superconductivity shall be excluded and it remains to be shown whether models with moderate diquark pairing could fulfill the thermal twin constraint. For the case of symmetric matter, we compare the resulting hybrid EOS with the flow constraint by Danielewicz et al. and find a a sensitivity of the onset density for deconfinement to the presence or absence of color superconductivity.

The familiar Bethe-Heitler process on the proton target $e+p\to e+p+\gamma$ is transformed into a laboratory for studying multiparticle entanglement. We discuss how bipartite and genuine tripartite entanglement between the final state electron, proton and photon are built up by successive $1\to 2$ and $2\to 2$ elementary interactions. We validate our argument by simulating events. Below 5 GeV center-of-mass energy, we identify more than 900 Greenberger-Horne-Zeilinger (GHZ) states and 1200 W states, each with a fidelty exceeding 99%.

Energy loss by fast partons induces a Mach-cone-like medium response as they propagate inside the hot quark-gluon plasma (QGP) in high-energy heavy-ion collisions. Because the QGP is nonuniform and its initial gradients generate collective flow, jet-induced medium response in this evolving system is also distorted by the flow and density gradient. This distortion leads to a broadened jet wake whose transverse width depends on the azimuthal angle of the jet propagation due to the elliptic anisotropy of the density gradient and the flow velocity in noncentral heavy-ion collisions. We propose and calculate the difference between the azimuth-dependent jet-hadron correlations for soft charged hadrons in in-plane and out-plane $\gamma$-jets as a measure of the elliptic broadening of the wake front and the deepening of the diffusion wake. We also study the sensitivity of this observable to the shear viscosity of the QGP. Experimental measurements of the azimuthal modulation of jet wakes induced by the wind of the elliptic flow at RHIC and LHC can provide additional constraints on the transport properties of the QGP.

Symbolic regression provides an interpretable machine-learning approach for constructing explicit analytic relations between physical inputs and observables. In this work, we develop symbolic-regression emulators for the isospin-dependent Boltzmann-Uehling-Uhlenbeck (IBUU) transport model and compare their performance with deep neural network (DNN) emulators. Using the same transport-model data employed in our previous emulator studies, we show that symbolic regression can reproduce the proton mid-rapidity slope $F_1$ of transverse flow $v_1$ and elliptic flow $v_2$ with accuracy comparable to that of DNNs, while providing explicit analytic expressions and substantially faster prediction once trained. We further demonstrate the use of symbolic regression in the reverse direction by constructing analytic relations that predict the in-medium nucleon-nucleon cross-section modification factor $X$ from the flow observables. Although the symbolic-regression models require substantially longer training times and exhibit greater run-to-run variation than DNNs, their analytic form and rapid evaluation make them promising tools for future transport-model sensitivity and uncertainty analyses.

Recent advances in ab initio nuclear theory, machine learning, and Bayesian inference, coupled with increasingly precise nuclear experiments and astrophysical observations, have enabled more robust constraints on fundamental descriptions of the nuclear interaction. Although well-established nonlinear regression methods can identify best-fit sets of low-energy constants at fixed resolution scale, they provide limited insight into the full underlying probability distributions of those constants. A central remaining challenge in nuclear theory is therefore to characterize full probability distributions of nuclear forces across resolution scales. In this work, we employ normalizing flows, a class of expressive generative machine learning models, to infer the joint probability distribution of two-body low-energy constants (LECs) in chiral effective field theory over a wide range of resolution scales. The resulting LEC distributions are shown to accurately reproduce experimental neutron-proton scattering phase-shift distributions. Furthermore, strong non-Gaussian correlations among LECs are revealed, indicating a nontrivial interplay among distinct short-range nuclear dynamics. This work establishes a general framework for constructing statistical ensembles of nuclear interactions that can be systematically constrained by future nuclear experiments and astrophysical observations.

The thermonuclear $^{45}$V($p,\gamma$)$^{46}$Cr reaction is the primary leakage pathway from the $^{44}$Ti–$^{45}$V quasi-equilibrium cluster during $\alpha$-rich freeze-out in core-collapse supernovae (CCSN), governing the final abundance of the $\gamma$-ray-emitting isotope $^{44}$Ti. A recent high-resolution $\gamma$-ray study [C. Cousins et al., Phys. Rev. Lett. 136, 252701 (2026)] identified ten previously unknown low-spin proton-unbound states in $^{46}$Cr, enabling the first experimentally constrained $^{45}$V($p,\gamma$)$^{46}$Cr reaction rate using the AME2020 mass excess, $\text{ME}(^{46}\text{Cr}) = -29472(11)$~keV. Here, we adopt the four-fold more precise CSRe mass excess $\text{ME}(^{46}\text{Cr}) = -29477.2(2.6)$~keV [M.~Wang et al., Phys. Rev. C 106, L051301 (2022)] to recalculate the reaction rate. Including proton capture on the ground and first two excited states of $^{45}$V alongside new shell-model proton spectroscopic factors, we reduce mass-related rate uncertainties to a subdominant level. The revised rate is up to 69% higher than that of Cousins et al. at $\alpha$-rich freeze-out temperatures ($T \simeq 1.5$–$2$~GK). CCSN nucleosynthesis calculations show this revised rate increases the ejected $^{44}$Ti yield by $\sim$26% in a $20\,M_\odot$ model compared to The et al. [ApJ 504, 500 (1998)], while causing negligible changes for the SN~1987A trajectory. We demonstrate that $^{44}$Ti production sensitivity is dictated by the ejecta electron fraction ($Y_e$): the reaction significantly affects proton-rich ejecta ($Y_e \approx 0.50$) but has little impact on neutron-rich ejecta ($Y_e \approx 0.496$), where lower free-proton abundances suppress reaction flow. This reconciles conflicting results from past sensitivity studies.

In this work, we investigate possible $\chi_{cJ}$–nuclear bound states with $J=0,1,2$ using in-medium mass shifts generated by virtual $D^{(*)}\bar{D}^{(*)}$ loops within an unquenched framework. The resulting $\chi_{cJ}$–nucleus potentials are constructed in the local density approximation, and the bound state spectra are calculated for $^{12}{\rm C}$, $^{16}{\rm O}$, $^{40}{\rm Ca}$, $^{90}{\rm Zr}$, $^{197}{\rm Au}$, and $^{208}{\rm Pb}$. Bound states are obtained for all systems considered. The $\chi_{c0}(1P)$ and $\chi_{c1}(1P)$ spectra are nearly degenerate, whereas the larger in-medium mass shift of $\chi_{c2}(1P)$ leads to deeper binding. Although the absolute bound state energies depend appreciably on the cutoff parameter, the energy differences relative to the $1s$ level are considerably less sensitive to it and exhibit a regular pattern that decreases approximately as $A^{-2/3}$ with increasing nuclear mass number. A cosh-type potential with a common nuclear geometry provides a compact description of these spectra. The predicted bound-state structures and level-spacing systematics could be investigated in future high-statistics near-threshold photoproduction experiments at the upgraded JLab facility.