The study of nuclear shape evolution and coexistence provides key insight into the nuclear interaction, particularly in semi-magic systems expected to be spherical. A high-precision Coulomb-excitation measurement of $^{116}$Sn yields a comprehensive set of electromagnetic matrix elements, enabling the determination of quadrupole moments of the $2_{1,2,3}^+$ states and intrinsic deformations of the $0_{1,2,3}^+$ states. The results provide an unambiguous and direct proof of the rare phenomenon of multiple-shape coexistence, revealing a weakly deformed ground state incompatible with the spherical shape.

By integrating the quantum Monte Carlo technique into the similarity renormalization group (SRG), we have developed a stochastic SRG framework (SRGQMC) capable of both free-space two-body and in-medium many-body evolutions. This approach circumvents the combinatorial tensor-space explosion of many-body flow equations by mapping continuous unitary transformations onto an ensemble of signed random walkers. We benchmark the SRGQMC against deterministic free-space SRG evolutions of realistic nucleon-nucleon (NN) interactions, as well as against in-medium SRG (IMSRG) many-body calculations with the Richardson pairing model at two- and three-body levels [IMSRG(2)/(3)]. While a deterministic extension to the four-body level [IMSRG(4)] remains unfeasible due to prohibitive computational costs, we have achieved the first IMSRG(4) calculation by using the stochastic technique, demonstrating a substantial improvement toward the full configuration-interaction limit. This stochastic framework provides a practical pathway to higher-order IMSRG calculations.

The p-process is responsible for the production of the stable neutron-deficient nuclei heavier than iron observed in the solar system. However, important nuclear uncertainties still limit our understanding of this nucleosynthesis process. Among the most significant are the nuclear level densities (NLDs) and photon strength functions (PSFs) entering the calculation of photodisintegration rates under supernova conditions. We investigate both model (systematic) and parameter (statistical) uncertainties affecting NLDs and PSFs and quantify their impact on p-process nucleosynthesis in type-Ia and type-II supernovae. Correlated model uncertainties are estimated using several NLD and PSF models that reproduce available experimental observables. Uncorrelated parameter uncertainties are evaluated with a backward-forward Monte Carlo approach, in which parameter variations are constrained by measured reaction rates before being propagated to unknown cross sections of neutron-deficient nuclei. The resulting uncertainties are propagated through p-process calculations while preserving model correlations. To identify the reactions driving abundance uncertainties, we combine regularized linear-response modeling, stability analysis, and contribution and interaction decompositions. We find that photoneutron-emission uncertainties dominate the overall uncertainty budget. The leading source of uncertainty arises from local parameter variations still compatible with current experimental constraints, highlighting the lack of constraining nuclear data in the neutron-deficient region. For many p-nuclei, the dominant contribution originates either from the photoneutron emission of the p-nucleus itself or from a nearby $(\gamma,n)$ reaction along the same isotopic chain. While improved nuclear models remain important, many key reactions involve stable or near-stable nuclei and should be experimentally accessible.

We report on the performance of six custom-made tapered LYSO crystals of unprecedented volume, which constitute a sector of the 19 radiation length electromagnetic calorimeter planned for the PIONEER experiment. The longitudinal response uniformity of each crystal was measured using radioactive sources before characterizing the energy and time resolution of the crystals in an array using a 20 to 80 MeV positron beam at the Paul Scherrer Institute. The array demonstrated an energy resolution better than 2 % for energies above 40 MeV and a time resolution better than 130 ps for energies above 30 MeV. The spatial resolution was measured in the central region of the array to be 4.9 mm at 70 MeV, and was extrapolated to 5.4 mm across a 30 mm radius region using Geant4 simulation. The measured properties satisfy the key design parameters of the PIONEER calorimeter for the measurement of rare pion decays.

Tensor-polarized parton distribution functions (PDFs) are calculated for the deuteron by using a convolution formalism, where the tensor-polarized PDFs are given by the corresponding nucleon's unpolarized PDFs convoluted with the tensor-polarized nucleon momentum distribution in the deuteron. These distributions are obtained at $Q^2=2.5$ GeV$^2$ in order to compare with the tensor-polarized PDFs which were determined by HERMES $b_1$ data. The obtained distributions are very different from the ones determined from the HERMES data, which indicates further studies are needed to clarify the difference, possibly by considering a new mechanism beyond the simple bound system of a proton and a neutron. The obtained PDFs $\delta_T q $ and $\delta_T \bar q$ are converted to the PDFs of the Trento convention $f_{1LL}^{\, q}$ and $f_{1LL}^{\, \bar q}$, and they are used for estimating the twist-3 PDFs $f_{LT}^{\, q}$ and $f_{LT}^{\, \bar q}$ by using a Wandzura-Wilczek-like relation. Because deep-inelastic-scattering experiments are under preparation for structure functions with a tensor-polarized deuteron target at the Thomas Jefferson National Accelerator Facility, and a Drell-Yan experiment will be possible at hadron accelerator facilities, such as the Fermi National Accelerator Laboratory, the obtained tensor-polarized PDFs will be tested experimentally.

We perform a Diffusion Monte Carlo study of fully light hexaquark containing three $u$ quarks and three $d$ quarks within a constituent-quark model. Both compact and baryon--baryon-like arrangements were considered separately. All compact hexaquark configurations are found well above their corresponding baryon--baryon thresholds, suggesting that deeply bound compact six-quark states are not favored in the light-quark sector within this model. By contrast, several dibaryon-like configurations lie close, but above, to the $NN$, $N\Delta$, and $\Delta\Delta$ thresholds and show spatial structures compatible with molecular states. One configuration exhibits two well-defined nucleon-like subclusters separated by several femtometers, closely resembling the spatial structure of the deuteron, although its calculated energy remains slightly above the corresponding threshold.

We study the effects of the running of the QCD coupling on the anti-collinear resummation in JIMWLK evolution in the linear (BFKL) regime. We determine the appropriate scale choice for the coupling entering the JIMWLK kernel, and derive the anti-collinearly resummed BFKL kernel, which includes running-coupling effects both in the resummation equation (i.e. DGLAP) and in the JIMWLK kernel proper. We find that the running of the coupling generally further slows down BFKL evolution, as expected. Surprisingly however the value of the generalized characteristic function at $\gamma=1$ is unaffected by the running coupling owing to subtle cancellations. We develop an approximation that allows us to use the generalized characteristic function to study the BFKL Green's function. Within this approximation we find that the Pomeron intercept in the anti-collinear regime is significantly reduced by both, the resummation and the running of the coupling.

We present a closed analytical description of the spin distributions of the fragments produced in low-energy induced and spontaneous fission. In our model the high fragment spins and the relative orbital angular momentum arise from the zero-point transverse wriggling and bending oscillations of the two pre-fragments, under the postulate that the fissioning system remains ``cold'' up to scission -- its available energy being stored as non-equilibrium deformation rather than as heat. From the probability distributions of the two modes we derive a closed expression for the spin distribution of each fragment and for its mean value. The decisive quantities are the fragment moments of inertia, which we evaluate in the hydrodynamic model from the non-equilibrium scission deformations reconstructed from the measured prompt-neutron multiplicities. Confronted with the recent data on $\rm ^{232}Th(n, f)$, $\rm ^{238}U(n, f)$, and $\rm ^{252}Cf(sf)$, the model reproduces both the magnitude of the mean spins and their characteristic sawtooth dependence on the fragment mass. Comparison with the statistical and microscopic approaches indicates that the differences for individual fragments can be traced largely to the deformation dependence of the moments of inertia.

Some atomic nuclei exhibit enhanced octupole collectivity, reflected in finite reflection-asymmetric multipole correlations rather than necessarily in a rigid static pear-shaped ground state. Low-energy studies indicate finite octupole strength in uranium-238, commonly interpreted as soft or vibrational in nature, in addition to its large prolate quadrupole collectivity~\cite{MCGOWAN1994569,KIBEDI:2002wxc}, in addition to its large prolate quadrupole collectivity. Here we investigate how such octupole correlations can be encoded in the initial geometry of relativistic heavy-ion collisions and mapped to final-state flow observables. Using state-of-the-art hydrodynamic calculations, we demonstrate quantitative sensitivity to octupole-induced features encoded in the initial-state geometry and suggest a modest octupole collectivity in uranium-238, confirmed by the latest high-energy experimental measurements~\cite{STAR:2025elk}. These findings provide as a complementary probe of odd-order nuclear collectivity and help constrain quark-gluon plasma initial conditions.