The Nuclear Physics European Collaboration Committee ( NuPECC, this http URL ) hosted by the European Science Foundation represents today a large nuclear physics community from 23 countries, 3 ESFRI (European Strategy Forum for Research Infrastructures) nuclear physics infrastructures and ECT* (European Centre for Theoretical Studies in Nuclear Physics and Related Areas), as well as from 4 associated members and 10 observers. As stated in the NuPECC Terms of Reference one of the major objectives of the Committee is: "on a regular basis, the Committee shall organise a consultation of the community leading to the definition and publication of a Long Range Plan (LRP) of European nuclear physics". To this end, NuPECC has in the past produced five LRPs: in November 1991, December 1997, April 2004, December 2010, and November 2017. The LRP, being the unique document covering the whole nuclear physics landscape in Europe, identifies opportunities and priorities for nuclear science in Europe and provides national funding agencies, ESFRI, and the European Commission with a framework for coordinated advances in nuclear science. It serves also as a reference document for the strategic plans for nuclear physics in the European countries.

Silicon photomultipliers (SiPMs) will be used to read out all calorimeters in the ePIC experiment at the Electron-Ion Collider (EIC). A thorough characterization of the radiation damage expected for SiPMs under anticipated EIC fluences is essential for accurate simulations, detector design, and effective operational strategies. In this study, we evaluate radiation damage for the specific SiPM models chosen for ePIC across the complete fluence range anticipated at the EIC, $10^8$ to $10^{12}$ 1-MeV $n_{\mathrm{eq}}$/cm$^2$ per year, depending on the calorimeter location. The SiPMs were irradiated using a 64 MeV proton beam provided by the University of California, Davis 76" Cyclotron. We measured the SiPM dark-current as a function of fluence and bias voltage and investigated the effectiveness of high-temperature annealing to recover radiation damage. These results provide a comprehensive reference for the design, simulation, and operational planning of all ePIC calorimeter systems.

The experimental observations that led to the quark structure of matter and the development of hadron physics are reviewed with emphasis on the discoveries of mesons and baryons, starting in the 1940s with the pion and kaon which mediate the strong hadronic force. The evidence for an internal structure of the hadrons consisting of two or three elementary spin 1/2 particles is reviewed. The discoveries of hadrons made of the heavier charm and bottom quarks are described. In 2003 more complex multi-quark hadrons began to emerge. The subsequent developments beyond the early 2000s are covered in the Review of Particle Physics (Phys. Rev. D 110 (2024) 030001). Given the very large number of observed hadrons, the choice of key experiments is somewhat subjective.

Measurement of the Chiral Magnetic Effect (CME) has been a popular topic of high-energy nuclear physics in the last decade. The flow correlation $\gamma$ between charged hadron pairs of the same and opposite charges and their difference $\Delta \gamma$ were measured to separate the CME-driven signal from the collective flow background especially second-order elliptic $v_{2}$. The STAR experiment have stepped further to the isobar experiment to compare $\gamma$ and $\Delta \gamma$ between Ru+Ru and Zr+Zr ~\cite{PhysRevC.105.014901}, which were theoretically expected to produce the same elliptic flow background but different CME signals. However, the measured flow backgrounds also differ between Ru+Ru and Zr+Zr, indicating more fine-tuning of RP and centrality definition necessary. This analysis applied the AMPT model~\cite{PhysRevC.72.064901} to simulate the same collision system and energy as the STAR isobar experiment. Since the AMPT model does not include magnetic field effects, we expect comparing its output between Ru+Ru and Zr+Zr collision systems can provide an insight of the possible bias of flow background definition, and help improve the measurement of CME signal in real experiments. Multiple combinations of centrality and flow definition were chosen to study how the $v_2$ and their difference would be affected, especially by varying the particles selection of charge versus neutral properties and broadening (pseudo-)rapidity regions, while STAR CME work relied on charged-only particles at central rapidity.

Despite numerous achievements and recent progress, nuclear physics is often (wrongly) considered an old field of research nowadays. However, developments in theoretical frameworks and reliable experimental techniques have made the field mature enough to explore many new frontiers. In this regard, extending existing knowledge to an emerging field of physics -- where particles interact with a relatively low-energy but high intensity field (intense enough so that multi-particle processes become comparable or more important than one-to-one processes) -- can lead to exciting discoveries. Investigations can be realized under a highly time-compressed beam source (e.g., particle sources generated by laser-matter interaction using high-power laser systems). Here we focus on a new scheme, where high-power laser systems are exploited as a driver to generate energetic ($\gamma$-ray) photons. Together with additional low-energy photons provided by a second, less intense laser, a multi-photon absorption scheme enables a very attainable manipulation of nuclear transitions including isomer pumping and depletion.

The precise determination of the free neutron lifetime is of great significance in modern precision physics. This key observable is linked to the mixing of up and down quarks via the Cabibbo-Kobayashi-Maskawa matrix element $V_{ud}$, and the abundance of primordial elements after the Big-Bang Nucleosynthesis. However, the two leading measurement techniques for the neutron lifetime currently yield incompatible results, a discrepancy referred to as the neutron lifetime puzzle. To address the systematic uncertainties arising from neutron interactions with material walls, the $\tau$SPECT experiment employs a fully magnetic trap for ultra-cold neutrons (UCNs). UCNs are extremely low-energy neutrons with typical velocities below $8\,\textrm{m/s}$, which can be manipulated using magnetic fields, gravity, and suitable material guides, whose surface can reflect them at any angle of incidence. To precisely study and characterize UCN behavior during production, guidance, storage, and detection in $\tau$SPECT, we have developed a dedicated simulation framework. This framework is built upon the externally developed UCN Monte Carlo software package PENTrack and is enhanced with two companion tools: one for flexible and parametrizable upstream configuration of PENTrack such that the simulation's input settings can be adjusted to reproduce the experimental observations. The second package is used for analyzing, visualizing, and animating simulation data. The simulation results align well with experimental data obtained with $\tau$SPECT at the Paul Scherrer Institute and serve as a powerful resource for identifying systematic uncertainties and guiding future improvements to the current experimental setup.

Calculations of neutrinoless $\beta\beta$ ($0\nu\beta\beta$) decay nuclear matrix elements (NMEs) are carried out within the interacting boson model (IBM) that is based on the nuclear energy density functional (EDF) theory. The Hamiltonian of the IBM that gives rise to the energies and wave functions of the ground and excited states of $0\nu\beta\beta$ decay emitter isotopes and corresponding final nuclei is determined by mapping the self-consistent mean-field deformation energy surface obtained with a given EDF onto the corresponding bosonic energy surface. The transition operators are formulated using the generalized seniority scheme, and the pair structure constants are determined by the inputs provided by the self-consistent calculations. The predicted values of the $0\nu\beta\beta$-decay NMEs with the nonrelativistic and relativistic EDFs are compared to those resulting from different many-body methods. Sensitivities of the predicted NMEs to the model parameters and assumptions, including those arising from the choice of the effective interactions, are discussed.

Cosmic-ray physics in the GeV-to-TeV energy range has entered a precision era thanks to recent data from space-based experiments. However, the poor knowledge of nuclear reactions, in particular for the production of antimatter and secondary nuclei, limits the information that can be extracted from these data, such as source properties, transport in the Galaxy and indirect searches for particle dark matter. The Cross-Section for Cosmic Rays at CERN workshop series has addressed the challenges encountered in the interpretation of high-precision cosmic-ray data, with the goal of strengthening emergent synergies and taking advantage of the complementarity and know-how in different communities, from theoretical and experimental astroparticle physics to high-energy and nuclear physics. In this paper, we present the outcomes of the third edition of the workshop that took place in 2024. We present the current state of cosmic-ray experiments and their perspectives, and provide a detailed road map to close the most urgent gaps in cross-section data, in order to efficiently progress on many open physics cases, which are motivated in the paper. Finally, with the aim of being as exhaustive as possible, this report touches several other fields -- such as cosmogenic studies, space radiation protection and hadrontherapy -- where overlapping and specific new cross-section measurements, as well as nuclear code improvement and benchmarking efforts, are also needed. We also briefly highlight further synergies between astroparticle and high-energy physics on the question of cross-sections.

We analyze joint factorial cumulants of protons and antiprotons in relativistic heavy-ion collisions and point out that they obey the scaling $\hat{C}_{nm}^{p,\bar{p}} \propto \langle N_p\rangle^n \langle N_{\bar{p}} \rangle^m$ as a function of acceptance when only long-range correlations are present in the system, such as global baryon conservation and volume fluctuations. This hypothesis can be directly tested experimentally without the need for corrections for volume fluctuations. We show that if correlations among protons and antiprotons are driven by global baryon conservation and volume fluctuations only, the equality $\hat{C}_{2}^{p} / \langle N_p \rangle^2 = \hat{C}_{2}^{\bar{p}} / \langle N_{\bar{p}} \rangle^2$ holds for large systems created in central collisions. We point out that the experimental data of the STAR Collaboration from phase I of RHIC beam energy scan are approximately consistent with the scaling $\hat{C}_{nm}^{p,\bar{p}} \propto \langle N_p \rangle^n \langle N_{\bar{p}} \rangle^m$, but the normalized antiproton correlations are stronger than that of protons, $-\hat{C}_{2}^{\bar{p}} / \langle N_{\bar{p}} \rangle^2 > -\hat{C}_{2}^{p} / \langle N_p \rangle^2$, indicating that global baryon conservation and volume fluctuations alone cannot explain the data. We also discuss high-order factorial cumulants which can be measured with sufficient precision within phase II of RHIC-BES.