We review model calculations of exclusive vector meson production in ultraperipheral heavy ion collisions. We highlight differences and similarities between different dipole models and leading twist shadowing calculations. Recent color glass condensate calculations are presented with focus on effects from nuclear structure and azimuthal anisotropies driven by interference effects.

Evidences are updated and strengthened for the two-scales picture of low-energy nucleon structure as a compact `hard' valence quark core surrounded by a `soft' cloud of quark-antiquark pairs (the meson cloud). These considerations are quantified by a spectral analysis of the mean-squared radii associated with the isoscalar and isovector electric form factors of the nucleon. Further supporting arguments come from corresponding studies of the axial and mass form factors and their inferred radii. Separating low-mass (mesonic) and high-mass (short-range) contributions in the spectral representations of each of these form factors, we conclude that a central core with an r.m.s. radius of about 1/2 fm results consistently as the common feature in all cases. Implications are discussed for baryonic matter at densities beyond that of equilibrium nuclear matter.

Starting from a complete set of relativistic nucleon-nucleon contact operators up to order $O(p^4)$ of the expansion in the soft (relative or nucleon) momentum $p$, we show that non-relativistic expansions of relativistic operators involve twenty-six independent combinations, two starting at $O(p^0)$, seven at order $O(p^2)$ and seventeen at order $O(p^4)$. This demonstrates the existence of two low-energy free constants that parameterize interactions dependent on the total momentum of the pair of nucleons $P$. The latter, through the use of a unitary transformation, can be removed in the two-nucleon fourth-order contact interaction of the Chiral Effective Field Theory, generating a three-nucleon interaction at the same order. Within a hybrid approach in which this interaction is considered together with the phenomenological potential AV18, we show that the LECs involved can be used to fit very accurate data on the polarization observables of the low-energy $p-d$ scattering, in particular the $A_y$ asymmetry.

Numerous phenomenological nuclear models have been proposed to describe specific observables within different regions of the nuclear chart. However, developing a unified model that describes the complex behavior of all nuclei remains an open challenge. Here, we explore whether novel symbolic Machine Learning (ML) can rediscover traditional nuclear physics models or identify alternatives with improved simplicity, fidelity, and predictive power. To address this challenge, we developed a Multi-objective Iterated Symbolic Regression approach that handles symbolic regressions over multiple target observables, accounts for experimental uncertainties and is robust against high-dimensional problems. As a proof of principle, we applied this method to describe the nuclear binding energies and charge radii of light and medium mass nuclei. Our approach identified simple analytical relationships based on the number of protons and neutrons, providing interpretable models with precision comparable to state-of-the-art nuclear models. Additionally, we integrated this ML-discovered model with an existing complementary model to estimate the limits of nuclear stability. These results highlight the potential of symbolic ML to develop accurate nuclear models and guide our description of complex many-body problems.

Chiral effective field theory ($\chi$EFT) is a powerful tool for studying electroweak processes in nuclei. I discuss $\chi$EFT calculations of three key nuclear electroweak processes: primordial deuterium production, proton-proton fusion, and magnetic dipole excitations of $^{48}\mathrm{Ca}$. This article showcases $\chi$EFT's ability to quantify theory uncertainties at the appropriate level of rigor for addressing the different precision demands of these three processes.

The validation, verification, and uncertainty quantification of computationally expensive theoretical models of quantum many-body systems require the construction of fast and accurate emulators. In this work, we develop emulators for auxiliary field diffusion Monte Carlo (AFDMC), a powerful many-body method for nuclear systems. We introduce a reduced-basis method (RBM) emulator for AFDMC and study it in the simple case of the deuteron. Furthermore, we compare our RBM emulator with the recently proposed parametric matrix model (PMM) that combines elements of RBMs with machine learning. We contrast these two approaches with a traditional Gaussian Process emulator. All three emulators constructed here are based on a very limited set of 5 training points, as expected for realistic AFDMC calculations, but validated against $\mathcal{O}(10^3)$ exact solutions. We find that the PMM, with emulator errors of only $\approx 0.1 \%$ and speed-up factors of $\approx 10^7$, outperforms the other two emulators when applied to AFDMC.

The strongly interacting system created in ultrarelativistic nuclear collisions behaves almost as an ideal fluid with rich patterns of the velocity field exhibiting strong vortical structure. Vorticity of the fluid, via spin-orbit coupling, leads to particle spin polarization. Due to the finite orbital momentum of the system, the polarization on average is not zero; it depends on the particle momenta reflecting the spatial variation of the local vorticity. In the last few years, this field experienced a rapid growth due to experimental discoveries of the global and local polarizations. Recent measurements triggered further development of the theoretical description of the spin dynamics and suggestions of several new mechanisms for particle polarization. In this review, we focus mostly on the experimental results. We compare the measurements with the existing theoretical calculations but try to keep the discussion of possible underlying physics at the qualitative level. Future measurements and how they can help to answer open theoretical questions are also discussed. We pay a special attention to the employed experimental methods, as well as to the detector effects and associated corrections to the measurements.

In this work we perform an exploratory study of the photoproduction of singlet QED bound states $(l^-l^+)_S$ in electron - ion collisions at the EicC, EIC and LHeC energies. The total cross - sections, event rates per year and rapidity distributions associated with the parapositronium, paramuonium and paratauonium production are estimated. Moreover, we consider the decay of these states in a two - photon system and implement kinematical cuts on the rapidities and energies of the photons in the final state. We demonstrate the paramuonium can, in principle, be observed for the first time in all these colliders and that the EIC is a potential collider to discover the paratauonium state.