The positive-parity structure of 152Sm has been investigated through high-statistics {\gamma}-ray spectroscopy following the (150Nd({\alpha},2n)152Sm reaction at Elab = 26 MeV. Several collective structures built on excited 0+ states have been extended through the observation of new levels and {\gamma}-ray transitions, and spin-parity assignments have been established using directional-correlation and linear-polarization measurements. Electromagnetic transition strengths (B(E2)), deduced from measured branching ratios and known level lifetimes, reveal pronounced collectivity among the excited configurations. The resulting level scheme provides evidence for a sequence of excited collective bands extending beyond the well-known ground-state and first excited 0+ structures. The excitation energies and transition strengths are examined within the framework of the X(5) critical-point description of the first-order U(5)-SU(3) shape-phase transition. In addition to the established X(5)-like features of the low-lying spectrum, the observed systematics of the higher-lying bands are found to be consistent with excited collective structures exhibiting X(5)-like characteristics. The results provide new constraints on the realization of critical-point behavior in finite nuclei and on the evolution of collectivity in the N=90 region.

The STAR experiment reports the first measurement of the $^4_{\Lambda}\hbox{He}$ hyper-nuclei yield as a function of rapidity and transverse momentum in 0-50% central Au+Au collisions at $\sqrt{s_{\rm{NN}}} =$ 3 GeV. The $^4_{\Lambda}\hbox{He}$ is reconstructed through its three-body decay channel, $^4_{\Lambda}\rm{He} \rightarrow {}^{3}\rm{He} + \rm{p} + \pi^-$, with a statistical significance of about 9.5 standard deviations. We find that the yield of $^4_{\Lambda}\hbox{He}$ as a function of rapidity is consistent with that of $^4_{\Lambda}\hbox{H}$, and the rapidity-dependent yield ratio of $^4_{\Lambda}\hbox{He}$/$^4_{\Lambda}\hbox{H}$ is consistent with that of $^3$He/t. All the measurements, as well as the transverse-momentum spectra, can be reasonably described by the JAM with a coalescence afterburner, suggesting a coalescence-based formation scenario for hyper-nuclei at this energy. The canonical thermal model reproduces the observed yield ratios but overpredicts the absolute hyper-nuclei yields.

A clear spectral enhancement, known as the ``$\Sigma N$ cusp'', has been observed near the $\Sigma N$ threshold in the $d(K^-, \pi^-)$ reaction. To understand the dynamical origin of this enhancement, the J-PARC E90 experiment aims to investigate the missing-mass spectrum with an unprecedented resolution of 0.4 MeV ($\sigma$). In this experiment, a Hyperon Time Projection Chamber (HypTPC) is utilized to detect charged decay products and suppress severe contamination from quasi-free (QF) background processes. While a conventional track multiplicity condition of three (Mt=3) effectively suppresses these QF events, it restricts the signal statistics to approximately 17\% and introduces a mass-dependent acceptance bias that distorts the spectrum. In contrast, events with a track multiplicity of two (Mt=2) offer roughly double the statistical power ($\sim$39\%) with minimal mass dependence, but they suffer from heavy background contamination. To fully exploit the Mt=2 events, we developed an innovative background suppression technique based on a neural network. By constructing a binary classification model using the HypTPC track topology and energy loss ($dE/dx$) as input features, we successfully discriminated the signal from QF backgrounds. This machine learning approach achieves a signal-to-noise ratio comparable to the strict Mt=3 condition while preserving the integrity of the spectral shape. By combining this independent ML-selected Mt=2 sample with the conventional Mt=3 sample, the total usable statistics are effectively doubled compared to traditional methods, significantly enhancing the sensitivity for determining the $\Sigma N$ cusp parameters.

The excitation scheme of the $^{65}$Zn ($Z = 30, N = 35$) nucleus has been probed following its population in the $^{63}$Cu($\alpha$,pn) reaction at E$_{beam}$ = 30 MeV and using an array of Compton suppressed HPGe clovers as the detection system. This work has identified several new transitions of the nucleus and have modified the placements of some of the previously known ones. The multipolarities and the electric/ magnetic nature of the observed $\gamma$-ray rays have been measured, using the conventional methodologies. The spin-parity assignments for the levels have consequently been made; some of the spin-parities are new while others are either validation of the existing values or are modified results based on the present analysis. The experimental level scheme exhibits collective as well as single particle structures. The measured level energies have been compared with those calculated in the framework of the large basis shell model using a model space of $p_{3/2}, f_{5/2}, p_{1/2}, g_{9/2}$ orbitals and two different interactions. The collective excitations of the nucleus were probed through the properties of its band structures and through the calculations of the Total Routhian Surface (TRS) for the associated deformations/ shapes. The results of this study brings out the essential features of evolving structural characteristics and developing collectivity with increasing number of nucleons outside a doubly-magic core and with their occupancy of deformation driving high-$j$ orbitals.

The era of precision ab initio nuclear theory has arrived, enabling uncertainty-quantified predictions for nuclear structure and for interactions with external probes directly from the underlying nuclear force and electroweak currents. This review highlights recent breakthroughs that extend ab initio calculations to the heavy nucleus $^{208}$Pb, to medium-mass systems with complex deformation, and to weakly-bound nuclei near the driplines. We also summarize ab initio calculations of nuclear responses for dark matter direct detection. Together, these advances demonstrate how ab initio methods can substantially reduce nuclear-physics uncertainties in searches for physics beyond the Standard Model, providing a more robust interpretation of current and forthcoming precision experiments.

We present the characterization of two fast, crystalline inorganic scintillators, silicon-doped gallium nitride (GaN:Si) and gallium-doped zinc oxide (ZnO:Ga), and compare their performance with cerium-doped yttrium aluminium perovskite (YAP:Ce) for in-vacuum alpha-detection applications that require high-performance timing, position, and energy resolution, such as 3D elemental mapping, medical imaging, and homeland security applications. In this paper, we propose ZnO:Ga and GaN:Si as high-performance drop-in replacements for the alpha detector in Associated Particle Imaging (API) systems. However, the results reported here also have wide applicability. Prior work has reported on polycrystalline forms of ZnO:Ga, which suffer from self-absorption. To our knowledge, GaN:Si has not been proposed to be used in API systems. We present room-temperature scintillation time constants obtained via X-ray-induced time-correlated single-photon counting for both proposed materials. They both exhibit exceedingly fast rise times of <15ps, and high brightness >1000ph/MeV with resolved alpha-peaks. Single-crystal ZnO:Ga and single-crystal GaN:Si yield single-component decays of 805ps and 32ps, respectively. Using a plastic scintillator reference setup, coincidence timing resolution (CTR) and detector timing resolution (DTR) measurements demonstrate a >3x improvement in timing resolution compared to traditional YAP:Ce. GaN:Si and ZnO:Ga exhibit (35(9))ps and (49(5))ps DTR, respectively, compared to(144(2))ps for conventional, single-crystal YAP:Ce. Finally, we evaluate their position resolution in an experimental setup designed for API and measure better than 0.2mm for YAP:Ce and approximately 1mm for GaN:Si. We obtain a position resolution of 0.3mm for ZnO:Ga from simulations. We also present alpha-induced ionoluminescence emission spectra that reveal direct, red-shifted near-bandgap emission.

The dimensionless tidal deformability, $\Lambda$, of neutron stars (NSs), inferred from gravitational-wave (GW) observations, has thus far been used primarily to constrain the pressure of dense matter near twice nuclear saturation density, leaving the core equation of state (EOS) largely inaccessible to inspiral-phase GW observations. We show that the core EOS can be probed directly through $\Lambda$ using a perturbative analysis of the dimensionless stellar-structure and tidal-response equations formulated in terms of scaled intrinsic variables, without invoking any specific EOS model. We uncover a remarkable EOS-insensitive scaling relation between $\Lambda$ and the central EOS parameter $\mathrm{X}\equiv P_{\rm c}/\varepsilon_{\rm c}$, where $P_{\rm c}$ and $\varepsilon_{\rm c}$ denote the central pressure and energy density, respectively. The relation is validated against a broad ensemble of physically viable EOSs. Applying it to tidal deformabilities inferred from events such as GW170817 enables a direct determination of $\mathrm{X}$. We further derive a tight lower bound, $\Lambda_{\rm{TOV}}\gtrsim 9.2^{+1.2}_{-1.2}$, for maximum-mass NSs along stable mass-radius sequences, quantitatively demonstrating that even the most compact stable NSs remain distinctly separated from black holes, for which $\Lambda_{\rm{BH}}=0$. These findings reveal a previously unrecognized connection between inspiral-phase tidal deformability and the core EOS, establishing a direct link between GW observables and the microphysics of ultradense matter in the strong-gravity regime. The resulting scaling establishes inspiral-phase tidal deformability as a direct and largely model-insensitive probe of the EOS of NS cores.

The Fayans energy density functional (EDF), and in particular its pairing sector, have been claimed to be able to reproduce the experimental data of charge radii in many instances. A particularly intriguing case is that of the $ \mathrm{Ca} $ isotopes between $ A = 40 $ and $ 48 $, where charge radii exhibit a "bell shape". In our work, we examine the microscopic origin of this behaviour. We prepare in total $ 25 $ paramerizations of the Fayans-like pairing interaction, that are equivalent in fulfilling the same criteria for the reproduction of empirical pairing gaps. We find that both the density and the density-gradient dependence of the pairing interaction are important to reproduce the well-known enhancement of charge radii in the open-shell nuclei, leading to the "bell shape" behaviour of $ \mathrm{Ca} $ isotopes. In particular, this originates from the repulsive nature of the rearrangement potential, and cannot simply be mocked up by a refit of the pairing strength. At the same time, we notice some drawbacks of the Fayans standard EDFs, that may call for investigating a more general form of it.

We propose a noninvasive SQUID-based polarimeter for the polarized proton beam in the Electron-Ion Collider (EIC) Hadron Storage Ring (HSR), exploiting the collective magnetic dipole moment of the bunches rather than scattering. The six-snake HSR lattice has synchronous-particle spin tune $\nu_s = 1/2$, placing the in-plane spin-precession signal at half the revolution frequency ($\sim$39 kHz), in the DC SQUID band. Three pickup channels (cosine-$\theta$ and sine-$\theta$ saddle loops for the transverse components, a coaxial axial gradiometer for the longitudinal one) reconstruct the full polarization vector $(P_x, P_y, P_z)$ in two complementary modes. Static mode, the default for continuous noninvasive monitoring, reads all three components: $P_y$ at the revolution frequency and the residual in-plane components at $\nu_s f_\mathrm{rev}$, bunch by bunch over an hours-long fill, including $P_z$, inaccessible to single-spin scattering polarimetry by parity conservation. Dynamic mode gives a precise polarization-magnitude measurement: a longitudinal kicker tips a small fraction of the polarization into the horizontal (ring) plane to produce a free-induction-decay (FID) signal, and many phase-locked tip-$\pi$-echo-restore cycles are summed coherently via a matched filter across all bunches, with $\mathcal{O}(\alpha^2/\pi^2) \sim 10^{-4}$ loss per cycle, negligible over a full $\delta P/P = 1\%$ measurement. For tipping angle $\alpha = 30$ mrad, polarization $P = 0.7$, and effective rms spin-tune spread $\sigma_{\nu_s}^\mathrm{eff} = 10^{-3}$ (coherence time $\sim$2 ms), the integration time to reach $\delta P/P = 1\%$ is about 18 s at injection and 5 min at flattop. The architecture extends to deuteron and $^3$He beams via species-specific spin-magnetic factors, with applications to storage-ring EDM searches.

We introduce novel full configuration interaction quantum Monte Carlo (FCIQMC) as an accurate many-body solver for $\textit{ab initio}$ nuclear structure calculations. This stochastic approach directly samples the exact wave function in the full configuration space, enabling high-fidelity treatment of high-order many-body correlations in strongly interacting nuclear systems. Using interactions from chiral effective field theory, we have computed ground-state energies and charge radii of $^4$He, $^8$Be, $^{12}$C and $^{16}$O with sub-percent-level many-body uncertainties. These results establish FCIQMC as a stochastic full-configuration-space solver capable of treating systems beyond the reach of the conventional no-core shell model, and as an accurate benchmark for truncated many-body expansion methods.

Pion exchange plays a fundamental role in nuclear structure and is responsible for tensor correlations and high-momentum components in nuclei. The $(p,p'\pi)$ reaction provides a unique opportunity to investigate pion dynamics under large-momentum-transfer conditions. Its three-body kinematics allows large momentum transfer to be achieved while keeping the excitation energy of the residual nucleus low. We investigate the kinematical properties of the $^{12}\mathrm{C}(p,p'\pi^+)^{12}\mathrm{B}$ reaction using Lorentz-invariant three-body phase-space calculations. The calculations were performed for a 392-MeV proton beam assuming a constant transition amplitude. The resulting momentum-transfer map and phase-space distribution identify experimentally accessible regions of large momentum transfer and provide guidance for optimizing a double-arm spectrometer experiment at RCNP. The present study establishes a model-independent kinematical foundation for future investigations of pion-induced correlations, high-momentum components, and pion dynamics in nuclei.

Precise knowledge of nuclear structure is essential across fundamental physics, yet probing these structures is notoriously difficult. To address this challenge, ultra-peripheral collisions (UPCs) provide a femtoscopic tomography for imaging the atomic nucleus. UPCs offer a pristine electromagnetic pathway: coherent vector-meson photoproduction generates patterns of diffraction and two-source interference that directly encode the nuclear spatial density. Turning these patterns into quantitative constraints is, however, a challenging inverse problem, complicated by correlated sensitivities to deformation and neutron skin, phase smearing, and experimental backgrounds. Here we introduce an interpretable Multitask deep-learning framework that maps transverse momentum distributions to multiple nuclear-structure indicators simultaneously and identifies the kinematic regions driving each inference. We demonstrate the approach with coherent $J/\psi$ photoproduction in $^{96}_{40}\text{Zr} + ^{96}_{40}\text{Zr}$ collisions, showing that the learned features separate diffraction-dominated and interference-dominated information and provide analysis-ready observables for future high-luminosity data.

The flavor diagonal CP violation of elementary particle physics contributes to the atomic, nuclear, and nucleon electric dipole moments (EDMs), T-violating neutron optics, and to the angular correlations of beta decay. In this contribution, we review the basics and the importance of CP violation in the search for new physics beyond the standard model, the recent progress in the quantification of the hadron level CP violation contributing to the aforementioned observables, and finally the current attempt to solve the strong CP problem without additional interactions and fields.

The scaling properties of the $\Delta\gamma$ correlator, guided by calculations from the Anomalous Viscous Fluid Dynamics (AVFD) model, are used to investigate charge separation in $p$+Au, $d$+Au, Ru+Ru, Zr+Zr, and Au+Au collisions at $\sqrt{s_{\mathrm{NN}}}=200$~GeV, and in $p$+Pb and Pb+Pb collisions at $\sqrt{s_{\mathrm{NN}}}=5.02$ and $2.76$~TeV. The ratio $\Delta\gamma/v_2$, where $v_2$ is the elliptic-flow coefficient, exhibits a common approximate $1/N_{\rm ch}$ scaling behavior for $p$+Au, $d$+Au, $p$+Pb, and Pb+Pb collisions, establishing a common multiplicity-dilution baseline for background-driven charge correlations. In contrast, Ru+Ru, Zr+Zr, and Au+Au collisions show significant deviations from this scaling trend. These violations are qualitatively similar to those obtained in AVFD calculations with an input chiral magnetic effect (CME) signal and point to an additional charge-separation component beyond that expected from the observed background scaling. Quantitative estimates indicate that, in mid-central collisions, the corresponding CME-sensitive fraction of the measured $\Delta\gamma/v_2$ is approximately 27\% for Au+Au collisions and roughly a factor of two smaller for Ru+Ru and Zr+Zr collisions, which exhibit similar magnitudes within uncertainties. The extracted magnitudes imply an expected Ru+Ru--Zr+Zr signal difference of only $\sim1.3\%$, suggesting limited sensitivity of the $\Delta\gamma$ correlator to the small difference expected between the isobar signals.

The internal charge and magnetization distributions of the proton are characterized by electromagnetic form factors GE and GM. They are experimentally extracted via Rosenbluth separation, which measures the elastic scattering of electrons and protons at multiple beam energies and angles at fixed momentum transfer Q2. Conventionally, form factor values are obtained by plotting reduced cross sections against the virtual photon polarization parameter epsilon and then extracting the slope and intercept of the best fit lines. An alternative visualization method, proposed by Hofstadter and Herman in 1960, plots GM2 vs. GE2 curves instead. The best fit values of GE2 and GM2 are immediately visible from the intersection region of the curves and their uncertainty bands. In this work, we apply both conventional and Hofstadter-Herman visualizations to classic 1994 SLAC elastic scattering data. We demonstrate that the Hofstadter-Herman method reveals previously obscured regions of form factor parameter space and highlights subtle experimental discrepancies among data sets. Our results motivate adopting this visualization method as a routine diagnostic cross-check at the Electron-Ion Collider and elsewhere to flag normalization shifts and related adjustments before they enter global fits.

Neutrinoless double $\beta$-decay (0$\nu\beta\beta$), if observed, would provide unequivocal evidence of physics beyond the Standard Model. $^{48}$Ca is an interesting candidate system to study because it has the largest Q value among all 2$\beta$ transitions and is also unstable against single $\beta$-decay. The observation of both $\beta$ and 2$\beta$-decay in the same isotope would provide a unique opportunity to benchmark theoretical calculations of $\beta$ and 2$\beta$-decay matrix elements and could provide insight on the quenching of the axial vector coupling constant, g$_A$. We performed a precise measurement of the $^{48}$Ca $\beta$-decay Q value using the TITAN Penning trap mass spectrometer at the TRIUMF facility. This was achieved through cyclotron frequency ratio measurements of $^{48}$Ca$^{+}$/$^{48}$Sc$^{+}$ and $^{48}$Sc$^{+}$/$^{48}$Ti$^{+}$ using the Time-of-Flight Ion Cyclotron Resonance technique. The $^{48}$Ca $\beta$-decay Q value was determined to be 279.14(50) keV, a factor of 10 more precise than the previous value given in the 2020 Atomic Mass Evaluation [Chin. Phys. C 45, 030003 (2021)]. This Q value was used to determine the $^{48}$Ca $\beta$-decay partial half-life, with the result $T_{1/2}^{\beta}$ = 5.09(5) x 10$^{20}$ ($g_{A}^{-2}$) y. Our $^{48}$Ca $\beta$-decay half-life was determined to a precision of 1%, a factor of 30 improvement compared to calculations with the previous Q value. Our result is marginally closer to the experimental lower limit $T_{1/2}^{\beta}$ > 1.1 x 10$^{20}$ y, but still a factor 5 longer. It is also a factor of 10 longer than the observed 2$\nu\beta\beta$ decay mode with $T_{1/2}^{2\nu\beta\beta} = 5.96^{+1.39}_{-1.08}$ x 10$^{19}$ y. Hence, it could be possible to observe $^{48}$Ca $\beta$-decay in future experiments, strengthening the potential importance of $^{48}$Ca to benchmark nuclear structure and 2$\beta$-decay studies.

Ultra-relativistic heavy-ion collisions create a hot and dense medium of deconfined quarks and gluons, the quark$-$gluon plasma (QGP), in which parton energy loss ("jet quenching") is a key probe of hot medium properties. While parton energy loss has been firmly established in large systems such as Pb$-$Pb and Au$-$Au collisions, no unambiguous direct evidence exists in smaller systems such as high-multiplicity p$-$Pb and pp collisions. To probe the onset of parton energy loss at intermediate system size, measurements of neutral-pion production are presented in this Letter for oxygen$-$oxygen (OO) and proton$-$oxygen (pO) collisions recorded with the ALICE detector in July 2025, relative to a pp baseline. The nuclear modification factor $R_{\rm OO}$ is suppressed relative to unity with a transverse-momentum dependence similar to that observed in Pb$-$Pb collisions, consistent with a previous CMS measurement in OO collisions with charged particles. As $R_{\rm OO}$ contains contributions from both cold and hot nuclear matter effects, $R_{\rm pO}$ is also presented in order to constrain cold nuclear matter (CNM) contributions. $R_{\rm pO}$ is found to be compatible with unity, indicating that CNM effects alone cannot account for the suppression observed in $R_{\rm OO}$. Final-state effects are isolated using the measured double ratio $R_{\rm OO} \left/ R_{\rm pO}^2 \right.$, which largely cancels CNM contributions and exhibits a significant suppression relative to expectations without energy loss at a 4.9$\sigma$ level. Theoretical models incorporating parton energy loss via different mechanisms predict a significant suppression of the $R_{\rm OO} \left/ R_{\rm pO}^2 \right.$ relative to unity, consistent with the data. These findings establish parton energy loss in OO collisions, extending experimental evidence for jet quenching to the smallest nuclear system studied to date.

Low-energy vibrational excitations associated with the fluctuation of quadrupole deformed shapes are discussed within the frame of state-of-the-art Configuration Interaction calculations, actually performed via the Quasi-particle Vacua Shell Model version of the Monte Carlo Shell Model. Recently, low-lying $\gamma$ bands in heavy strongly deformed nuclei were shown to be rotational $K^P$ = 2$^+$ excitations of triaxially deformed states (see T. Otsuka \etal, Eur. Phys. J. A 61, 126 (2025)) rather than vibrational excitations as traditionally interpreted. In this context, it is important to identify possible low-lying vibrational excitations and to characterize the excitation energy at which they emerge. Focusing on two typical examples, $^{166}$Er and $^{162}$Dy, vibrational states are indeed identified above the $\gamma$ band using an extended version of the so-called T-plot. The phenomenon of shape coexistence is also shown to produce low-lying states below such vibrational band heads. These results suggest novel and rich structures in heavy deformed nuclei. While experimental counterparts are seen for some of such states, others are predictions opening doors to future dedicated experiments.

Measurements of SIDIS multiplicities for $\pi^+$ and $\pi^-$ from proton and deuteron targets are reported on a grid of hadron kinematic variables $z$, $P_{T}$, and $\phi^{*}$ for leptonic kinematic variables in the range $0.3

Dark matter remains one of the most compelling open problems in modern physics, motivating experimental searches for new light, weakly coupled particles beyond the Standard Model. Despite extensive efforts employing diverse detection strategies, large regions of parameter space remain unexplored. We report a high-statistics laboratory search for invisible decay modes in nuclear $\gamma$-ray cascades using approximately $100~\mathrm{kg}$ of CsI(Tl) scintillators operated at Texas A\&M University. The experiment employs a high-activity ${}^{46}$Sc radioactive source and a ``missing-$\gamma$'' technique, in which the absence of a photon from a well-identified cascade serves as a signature of new physics. Unlike appearance-disappearance experiments, this approach requires only a single photon conversion into a dark-sector particle, enabling sensitivity to significantly weaker couplings. The setup provides simultaneous sensitivity to a broad class of light dark-sector candidates, including axions and axion-like particles, dark scalars, and dark photons in the $0.1 - 1 \text{ MeV}$ mass region. Through careful control of detector containment, energy resolution, and environmental backgrounds, we exclude certain regions on the previously explored parameter space. With foreseeable improvements in detector volume and systematic uncertainty control, this technique has the potential to probe currently unexplored parameter space for axion-like particles and light dark scalars.

Three-nucleon force (3NF) plays an important role in understanding the structure of finite nuclei and the saturation properties of infinite nuclear matter. More specifically, 3NF should be necessary for each two-nucleon force (2NF) to obtain more accurate description of nuclear systems. 3NF derived from the chiral effective field theory has been successful in $\textit{ab initio}$ calculations of atomic nuclei. Most of established chiral nuclear forces have a nonlocal form in the momentum space. In this work, we construct a companion chiral 3NF specifically tailored to the new Idaho local position-space 2NF, and calculate binding energies and radii of nuclei up to $^{132}$Sn. We find that a chiral 3NF with hybrid local and nonlocal regulators has advantages in improving the nuclear structure calculations of both binding energies and radii with the new Idaho 2NF. The two low-energy constants of 3NF are constrained by the ground-state energies of $^3$H and $^{16}$O as suggested in a recent work.

Adiabatic fast passage (AFP) is a rapid method for reversing nuclear polarization and manipulating spin populations in polarized solid targets, avoiding the long repolarization times associated with dynamic nuclear polarization (DNP). We report AFP measurements in a 5~T, 1~K polarized-target system for irradiated $^{15}$NH$_3$, irradiated $^{14}$ND$_3$, and butanol-based materials prepared either with TEMPO doping or by irradiation. We also present a joint manipulated-lineshape analysis for spin-1 targets and demonstrate that vector and tensor polarizations can be extracted from AFP-manipulated deuteron NMR spectra even when the populations are not described by a single Boltzmann spin temperature. Finally, we report a reproducible polarization- and direction-dependent AFP response in a large irradiated $^{15}$NH$_3$ sample. These ammonia results are presented as empirical observations under the specific sample-coil conditions of the experiment, with possible circuit-mediated mechanisms such as radiation damping or superradiant behavior discussed but not assigned as a definitive cause.

We present the TQ4Q2.0 fragmentation functions for the production of all-heavy (fully heavy) $S$-wave tetraquarks ($T_{4Q}$) with scalar ($0^{++}$), axial-vector ($1^{+-}$), and tensor ($2^{++}$) quantum numbers in high-energy hadronic collisions. This work extends the previous TQ4Q1.1 framework by incorporating nonconstituent heavy-quark contributions and introducing a replica-based uncertainty-quantification strategy derived from multi-scale variations (MHOUs). The construction follows a nonrelativistic QCD factorization approach, combining gluon- and heavy-quark-initiated fragmentation channels at leading power. Initial-scale inputs are modeled through updated potential-inspired wave functions, while the subsequent DGLAP evolution is performed via the threshold-aware HF-NRevo scheme. A comprehensive systematic analysis of uncertainties is carried out, with contributions from color-composite long-distance matrix elements (LDMEs) and perturbative multiscale inputs. The resulting TQ4Q2.0 grids, publicly released in LHAPDF6 format, provide the first complete phenomenological set for all-heavy exotics, enabling precise studies of all-charm tetraquark production and jet-associated observables within the JETHAD environment. This article completes the high-energy resummation-driven generation of the TQ4Q program and establishes a definitive baseline for future collider-oriented analyses of all-heavy multiquark dynamics.