Modeling fission properties, such as barriers and rates, is highly challenging. The most microscopic methods available are based on energy density functionals (EDFs) and rely on a limited set of collective coordinates to describe the evolution of a fissioning nucleus from its ground state to scission. Leveraging the efficiency of the MOCCa nuclear structure code and the predictive power of the BSkG3 EDF, we systematically study fission properties of the heaviest nuclei (roughly 3,300) accounting for (1) axial, triaxial and octupole moment; (2) all nuclei, including odd and odd-odd systems; and (3) fission paths determined via the least-action principle. We present the set of primary fission barriers and spontaneous fission half-lives we obtain and discuss their implications for r-process nucleosynthesis.

We demonstrate that energy density functional (EDF) models for atomic nuclei can achieve a level of accuracy in describing fission properties comparable to that of more phenomenological approaches, while maintaining predictive power for nuclear ground-state observables. Using multiple BSkG parameterizations, we conducted a comprehensive benchmark of spontaneous fission (SF) half-lives and compare our predictions with all available experimental data - 136 values - from both ground states and fission isomers. Leveraging the quality of the BSkG functionals and explicitly accounting for more than two collective degrees of freedom in our calculations, we attained one of the closest agreements with experiment reported by an EDF model so far—within less than four orders of magnitude. This level of accuracy, together with the computational efficiency of our method, opens the way to systematic studies of the thousands of neutron-rich nuclei crucial for modeling r-process nucleosynthesis and the inclusion of SF half-lives into the fitting protocol of future EDF models.

The quasiparticle random phase approximation or QRPA has been a foundational many-body technique for decades across quantum chemistry, condensed matter and nuclear physics. Although computing power has increased and the matrix-free Finite Amplitude Method (FAM) exists, the computational complexity of FAM-QRPA calculations remains a limiting factor for the generation of linear response data on atomic nuclei that are crucial for several research fields. In this work, we establish that the FAM-QRPA equations are inherently suited to a reduced order modelling framework and can be emulated efficiently. Moreover, we present a greedy snapshot selection strategy that leverages the reduced cost of FAM-QRPA calculations when the imaginary part of the excitation frequency is large. Even when accounting for its construction, the resulting emulator accelerates strength function calculations by significantly more than an order of magnitude. We demonstrate that this framework and its speed-up generalize to light and heavy nuclei, different numerical representations, and diverse nuclear models including chiral EFT and configuration-interaction shell model approaches, as well as Skyrme, Gogny, and relativistic energy density functionals.

We report on mass measurements of three long-lived states in 114Rh performed with the JYFLTRAP Penning-trap mass spectrometer: the ground state and two isomers with estimated half-lives of about one second. The used Phase-Imaging Ion-Cyclotron-Resonance technique allowed for the discovery of a so far unknown second longlived isomer. All three states were produced directly in proton-induced fission on a uranium target, whereas only the isomeric states were populated in the beta decay of the 114Ru ground state with spin-parity 0+. We propose spin-parity assignments of (6-) for the ground state, and (3+) and (0-) for the isomers. They resolve the puzzle of anomalous fission yields of this isotope despite the existing literature assigning a low angular momentum to the ground state. The experimental evidence is further supported by a detailed analysis based on mean-field calculations with the BSkG3 model. As for many other nuclei in this mass region, considering triaxial shapes is decisive for the interpretation of low-lying states of this nucleus. The discovery of a new isomer in 114Rh and our theoretical work challenge the currently adopted spin-parity assignments in this and several other odd-odd neutron-rich Rh isotopes.

Nuclear physics influences are present throughout the Universe at every scale. Over the past decades, significant efforts have been dedicated to various fields related to nucleosynthesis and stellar evolution. These include both experimental and theoretical nuclear physics, as well as ground-and space-based astronomical observations and astrophysical modeling. Despite numerous successes, major challenges and unresolved questions persist. Concerning nuclear physics, despite the remarkable efforts of experimentalists in studying unstable targets, it remains unlikely that we will be able to measure the structure and interaction properties of all astrophysically relevant nuclei in the near future. Therefore, further progress must rely on theoretical approaches. The necessary ingredients (properties of cold and hot nuclei, nuclear level densities, optical potentials, photon strength functions, fission properties, β-strength functions) should ideally be derived from global, universal and microscopic models. New progress based on mean-field models are described and their impact on nuclear reaction rates of astrophysical relevance discussed and on nucleosynthesis illustrated.

Providing reliable data on the properties of atomic nuclei and infinite nuclear matter to astrophysical applications remains extremely challenging, especially when treating both properties coherently within the same framework. Methods based on energy density functionals (EDFs) enable manageable calculations of nuclear structure throughout the entire nuclear chart and of the properties of infinite nuclear matter across a wide range of densities and asymmetries. To address these challenges, we present BSkG4, the latest Brussels-Skyrme-on-a-Grid model. It is based on an EDF of the extended Skyrme type with terms that are both momentum and density-dependent, and refines the treatment of ^1S_0 nucleon pairing gaps in asymmetric nuclear matter as inspired by more advanced many-body calculations. The newest model maintains the accuracy of earlier BSkGs for known atomic masses, radii and fission barriers with rms deviations of 0.633 MeV w.r.t. 2457 atomic masses, 0.0246 fm w.r.t. 810 charge radii, and 0.36 MeV w.r.t 45 primary fission barriers of actinides. It also improves some specific pairing-related properties, such as the ^1S_0 pairing gaps in asymmetric nuclear matter, neutron separation energies, Q_β values, and moments of inertia of finite nuclei. This improvement is particularly relevant for describing the r-process nucleosynthesis as well as various astrophysical phenomena related to the rotational evolution of neutron stars, their oscillations, and their cooling.

High-energy collisions involving the A=96 isobars ^96Zr and ^96Ru have been performed in 2018 at Brookhaven National Laboratory's Relativistic Heavy Ion Collider (RHIC) as a means to search for the chiral magnetic effect in QCD. This would manifest itself as specific deviations from unity in the ratio of observables taken between ^96Zr+^96Zr and ^96Ru+^96Ru collisions. Measurements of such ratios (released at the end of 2021) indeed reveal deviations from unity, but these are primarily caused by the two collided isobars having different radial profiles and intrinsic deformations. To make progress in understanding RHIC data, nuclear physicists across the energy spectrum gathered in Heidelberg in 2022 as part of an EMMI Rapid Reaction Task Force (RRTF) to address the following question. Does the combined effort of low-energy nuclear structure physics and high-energy heavy-ion physics enable us to understand the observations made in isobar collisions at RHIC?

We present the first measurements with a new collinear laser spectroscopy setup at the Argonne Tandem Linac Accelerator System, utilizing its unique capability to deliver neutron-rich refractory metal isotopes produced by the spontaneous fission of ^{252}Cf. We measured isotope shifts from optical spectra for nine radioactive ruthenium isotopes ^{106-114}Ru, reaching deep into the mid-shell region. The extracted charge radii are in excellent agreement with predictions from the Brussels-Skyrme-on-a-Grid models that account for the triaxial deformation of nuclear ground states. We show that triaxial deformation impacts charge radii in models that feature shell effects, in contrast to what could be concluded from a liquid drop analysis. This indicates that this exotic type of deformation should not be neglected in regions where it is known to occur, even if its presence cannot be unambiguously inferred through laser spectroscopy.

The quantum-mechanical nuclear-shell structure determines the stability and limits of the existence of the heaviest nuclides with large proton numbers Z ≳ 100 (refs. 1-3). Shell effects also affect the sizes and shapes of atomic nuclei, as shown by laser spectroscopy studies in lighter nuclides4. However, experimental information on the charge radii and the nuclear moments of the heavy actinide elements, which link the heaviest naturally abundant nuclides with artificially produced superheavy elements, is sparse5. Here we present laser spectroscopy measurements along the fermium (Z = 100) isotopic chain and an extension of data in the nobelium isotopic chain (Z = 102) across a key region. Multiple production schemes and different advanced techniques were applied to determine the isotope shifts in atomic transitions, from which changes in the nuclear mean-square charge radii were extracted. A range of nuclear models based on energy density functionals reproduce well the observed smooth evolution of the nuclear size. Both the remarkable consistency of model prediction and the similarity of predictions for different isotopes suggest a transition to a regime in which shell effects have a diminished effect on the size compared with lighter nuclei.

Precision mass measurements of 104Y, 106Zr, 104,104m,109Nb, and 111,112Mo have been performed with the JYFLTRAP double Penning trap mass spectrometer at the Ion Guide Isotope Separator On-Line facility. The order of the long-lived states in 104Nb was unambiguously established. The trend in two-neutron separation energies around the N=66 neutron midshell appeared to be steeper with respect to the Atomic Mass Evaluation 2020 extrapolations for the 39Y and 40Zr isotopic chains and less steep for the 41Nb chain, indicating a possible gap opening around Z=40. The experimental results were compared to the BSkG2 model calculations performed with and without vibrational and rotational corrections. All of them predict two low-lying minima for 106Zr. While the unaltered BSkG2 model fails to predict the trend in two-neutron separation energies, selecting the more deformed minima in calculations and removing the vibrational correction, the calculations are more in line with experimental data. The same is also true for the 21+ excitation energies and differences in charge radii in the Zr isotopes. The results stress the importance of improved treatment of collective corrections in large-scale models and further development of beyond-mean-field techniques.

Background: The low-energy enhancement observed recently in the deexcitation gamma -ray strength functions, suggested to arise due to the magnetic dipole (M1) radiation, motivates theoretical efforts to improve the description of M1 strength in available nuclear structure models. Reliable theoretical predictions of nuclear dipole excitations are of interest for different nuclear applications and in particular for nuclear astrophysics, where the calculations of radiative capture cross sections often resort to theoretical gamma strength functions. The quasiparticle random-phase approximation (QRPA) approach is arguably the most widely spread microscopic tool in this context since it can be applied to heavy nuclei and at the scale of the entire chart. Purpose: We aim to benchmark the performance of QRPA calculations with respect to M1 gamma strength functions, with a special emphasis on the description of the low-energy effects observed in the deexcitation strength. Methods: We investigate the zero-temperature and finite-temperature (FT) magnetic dipole strength functions computed within the QRPA and compare them to those obtained from exact diagonalizations of the same Hamiltonian in restricted orbital spaces. Our sample consists of 25 spherical and deformed nuclei, with masses ranging from A = 26 to A = 136, for which the exact diagonalization of the respective effective Hamiltonian in three different valence spaces remains feasible. Results: We find a reasonable agreement for the total photoabsorption strengths between both many-body methods but show that the QRPA distributions are shown to be systematically shifted down in energy with respect to exact results. Photoemission strengths obtained within the FT-QRPA formalism appear insufficient to explain the low-energy enhancement of the M1 strength functions evidenced by the exact diagonalization approach. Conclusions: We ascribe the problems encountered in the zero- and finite-temperature QRPA calculations to the lack of correlations in the nuclear ground state and to the truncation of the many-body space. In particular, the latter prevents obtaining the sufficiently high level density to produce the low-energy enhancement of the M1 strength function, making the (FT-)QRPA approach unsuitable for predictions of such effects across the nuclear chart.

In the framework of nuclear energy density functional (EDF) methods, many nuclear phenomena can be related to the deformation of intrinsic states. Their accurate modeling relies on the correct description of the change of nuclear binding energy with deformation. The two most important contributions to the deformation energy have their origin in shell effects and the surface energy coefficient of nuclear matter. In a first step, we build nine series of parametrizations with a systematically varied surface-energy coefficient a_surf for three frequently-used options for the CM correction (none, one-body term only, full one-body and two-body contributions) combined with three values for the isoscalar effective mass m^*_0/m (0.7, 0.8, 0.85) and analyse how well each of these parametrizations can be adjusted to the properties of spherical nuclei and infinite nuclear matter. In a second step, we performed additional fits without the constraint on surface energy, adding one ``best-fit" parametrization to each of the nine series. We then benchmark these parametrizations to the deformation properties of heavy nuclei by means of three-dimensional Hartree-Fock-Bogoliubov calculations that allow for non-axial and/or non-reflection symmetric configurations. We perform a detailed correlation analysis between surface and volume properties of nuclear matter using the nine series of parametrizations. The best fits out of each series are then benchmarked on the fission barriers of Pu240 and Hg180, as well as on the properties of deformed states at normal and superdeformation for actinides and nuclei in the neutron-deficient Hg region. (see paper for full abstract)

We calculate the magnetic dipole gamma -ray strength functions in a chain of even-mass neodymium isotopes 144-152Nd in the framework of the configuration-interaction (CI) shell model. We infer the strength function by applying the maximum entropy method (MEM) to the exact imaginary-time response function calculated with the shell-model Monte Carlo (SMMC) method. The success of the MEM depends on the choice of a good strength function as a prior distribution. We investigate two choices for the prior strength function: the static path approximation (SPA) and the quasiparticle random-phase approximation (QRPA). We find that the QRPA is a better approximation at low temperatures (i.e., near the ground state), while the SPA is a better choice at finite temperatures. We identify a low-energy enhancement (LEE) in the MEM deexcitation M1 strength functions of the even-mass neodymium isotopes and compare with recent experimental results for the total deexcitation gamma -ray strength functions. The LEE is already seen in the SPA strength function but not in the QRPA strength function, indicating the importance of large-amplitude static fluctuations around the mean field in reproducing the LEE. Our method is currently the only one that can reproduce LEE in heavy open-shell nuclei where conventional CI shell-model calculations are prohibited. With the onset of deformation as number of neutrons increases along the chain of neodymium isotopes, we observe that some of the LEE strength transfers to a low-energy excitation, which we interpret as a finite-temperature scissors mode. We also observe a finite-temperature spin-flip mode.

Large-scale models of nuclear structure are currently the only way to provide consistent datasets for the many properties of thousands of exotic nuclei that are required by nucleosynthesis simulations. In [W. Ryssens et al., Eur. Phys. J. A 58, 246 (2022)], we recently presented the new BSkG2 model based on an energy density functional of the Skyrme type. Relying on a flexible three-dimensional coordinate representation of the nucleus, the model takes into account both triaxial deformation and time-reversal symmetry breaking. BSkG2 achieves a state-of-the-art global description of nuclear ground state (g.s.) properties and reproduces in particular the known masses with a root-mean-square (rms) deviation of 678 keV. Moving beyond g.s. properties, the model also reproduces all empirical values for the primary and secondary barriers as well as isomer excitation energies of actinide nuclei with rms deviations below 500 keV, i.e. with unprecedented accuracy. Here we discuss in detail the extension of our framework to the calculation of the fission barriers of 45 actinide nuclei, including odd-mass and odd-odd systems. We focus in particular on the impact of symmetry breaking which is key to the accuracy of the model: we allow systematically for axial, reflection and time-reversal symmetry breaking. The effect of the latter on the fission properties of odd-mass and odd-odd nuclei is small, but we find that allowing for shapes with triaxial or octupole deformation, as well as shapes with both, is crucial to achieving this accuracy. The numerical accuracy of our coordinate space approach, the variety of nuclear configurations explored and the simultaneous successful description of fission properties and known masses makes BSkG2 the tool of choice for the large-scale study of nuclear structure.

We report on precision mass measurements of $^{113,115,117}$Ru performed with the JYFLTRAP double Penning trap mass spectrometer at the Accelerator Laboratory of University of Jyv\"askyl\"a. The phase-imaging ion-cyclotron-resonance technique was used to resolve the ground and isomeric states in $^{113,115}$Ru and enabled for the first time a measurement of the isomer excitation energies, $E_x(^{113}$Ru$^{m})=100.5(8)$ keV and $E_x(^{115}$Ru$^{m})=129(5)$ keV. The ground state of $^{117}$Ru was measured using the time-of-flight ion-cyclotron-resonance technique. The new mass-excess value for $^{117}$Ru is around 36 keV lower and 7 times more precise than the previous literature value. With the more precise ground-state mass values, the evolution of the two-neutron separation energies is further constrained and a similar trend as predicted by the BSkG1 model is obtained up to the neutron number $N=71$.

We report here the current developments on the Brussels-Skyrme-on-a-Grid (BSkG) atomic mass models. In comparison with our previous models, BSkG3 improves the infinite nuclear matter (INM) properties which opens its applications to neutron stars. The results presented here show that BSkG3 preserve the excellent agreement with experimental nuclear masses and radii, together with fission barriers of actinides obtained by BSkG1 and BSkG2, while the nuclear matter properties are considerably improved.

Accurate estimates of the binding energy of nuclei far from stability that cannot be produced in the laboratory are crucial to our understanding of nuclear processes in astrophysical scenarios. Models based on energy density functionals have shown that they are capable of reproducing all known masses with root-mean-square error better than 800 keV, while retaining a firm microscopic foundation. However, it was recently pointed out in [M. Hukkanen et al., arXiv:2210.10674] that the recent BSkG1 model fails to account for a contribution to the binding energy that is specific to odd-odd nuclei, and which can be studied by using appropriate mass difference formulas. We analyse here the (lacking) performance of three recent microscopic mass models with respect to such formulas and examine possibilities to remedy this deficiency in the future.

Nuclear astrophysics is a multi-disciplinary field with a huge demand for nuclear data. Among its various fields, stellar evolution and nucleosynthesis are clearly the most closely related to nuclear physics. The need for nuclear data for astrophysics applications challenges experimental techniques as well as the robustness and predictive power of present nuclear models. Despite impressive progress for the last years, major problems and puzzles remain. In the present contribution, only a few nuclear astrophysics specific aspects are discussed. These concern some experimental progress related to the measurement of key reactions of relevance for the so-called sand p-processes of nucleosynthesis, the theoretical effort in predicting nuclear properties of exotic neutron-rich nuclei of interest for the r-process nucleosynthesis, and the recent introduction of machine learning techniques in nuclear astrophysics applications.

The low-energy enhancement observed recently in the deexcitation gamma-ray strength functions, suggested to arise due to the magnetic dipole radiation, motivates theoretical efforts to improve the description of M1 strength in available nuclear structure models. Reliable theoretical predictions of nuclear dipole excitations are of interest for different nuclear applications and in particular for nuclear astrophysics, where the calculations of radiative capture cross sections often resort to theoretical strength functions. We aim to benchmark many-body methods in their description of the M1 strength functions, with a special emphasis on the low-energy effects observed in the deexcitation strength. We investigate the zero-temperature and finite-temperature magnetic dipole strength functions computed within the quasiparticle random-phase approximation and compare them to those from exact diagonalizations of the same Hamiltonian in restricted orbital spaces. The study is carried out for a sample of 25 spherical and deformed nuclei which can be described by diagonalization of the respective effective Hamiltonian in three different valence spaces. A reasonable agreement is found for the total photoabsorption strengths while the QRPA distributions are shown to be systematically shifted down in energy with respect to exact results. Photoemission strengths obtained within the FT-QRPA appear insufficient to explain the low-energy enhancement of the M1 strength functions. The problems encountered in QRPA calculations are ascribed to the lack of correlations in the nuclear ground state and to the truncation of the many-body space. In particular, the latter prevents obtaining the sufficiently high level density to produce the low-energy enhancement of the strength function, making the (FT-)QRPA approach unsuitable for predictions of such effects across the nuclear chart.

State-of-the-art hydrodynamic simulations of the quark-gluon plasma are unable to reproduce the elliptic flow of particles observed at the BNL Relativistic Heavy Ion Collider (RHIC) in relativistic ^{238}U+^{238}U collisions when they rely on information obtained from low-energy experiments for the implementation of deformation in the colliding ^{238}U ions. We show that this is due to an inappropriate treatment of well-deformed nuclei in the modeling of the initial conditions of the quark-gluon plasma. Past studies have identified the deformation of the nuclear surface with that of the nuclear volume, though these are different concepts. In particular, a volume quadrupole moment can be generated by both a surface hexadecapole and a surface quadrupole moment. This feature was so far neglected in the modeling of heavy-ion collisions, and is particularly relevant for nuclei like ^{238}U, which is both quadrupole deformed and hexadecapole deformed. With rigorous input from Skyrme density functional calculations, we show that correcting for such effects in the implementation of nuclear deformations in hydrodynamic simulations restores agreement with BNL RHIC data. This brings consistency to the results of nuclear experiments across energy scales, and demonstrates the impact of the hexadecapole deformation of ^{238}U on high-energy collisions.