Researcher profile

S. Goriely

S. Goriely contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

astro-ph.SR nucl-th astro-ph.HE nucl-ex physics.ins-det

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author

13works

0followers

5topics

4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

BZPEER is loading the nearby papers, people, topics and institutions for this page.

preprint2026arXiv

Improved microscopic nuclear level densities within the triaxial Hartree-Fock-Bogoliubov plus combinatorial method

New developments have been brought to our energy-, spin- and parity-dependent nuclear level densities based on the microscopic combinatorial method. Our new calculation is based on the BSkG3 mean-field model which relies on a three-dimensional coordinate-space representation of the nucleus, allowing for the spontaneous breaking of ground state rotational, axial and reflection symmetry. In particular, we now account for the impact of possible triaxial deformation of nuclear ground states on the level density. This has two effects on our calculations: the additional freedom of the single-particle levels affects the intrinsic level density while the absence of a rotational symmetry axis results in a larger collective correction. The present model reproduces the experimental s- and p-wave neutron resonance spacings with a degree of accuracy comparable to that of the best global models available. It is also shown that the model gives a reliable extrapolation at low energies where experimental data on the cumulative number of levels can be extracted. The predictions are also in good agreement with the experimental data extracted from the Oslo method. Total level densities for more than 8500 nuclei are made available in a table format for practical applications. For the nuclei for which experimental s-wave spacings and enough low-lying states exist, renormalization factors are provided to reproduce simultaneously both observables. The same combinatorial method is used to estimate the nuclear level densities at the fission saddle points of actinides and at the shape isomer deformation. Finally, the new nuclear level densities are applied to the calculation of radiative neutron capture cross sections and compared with those obtained with our previous combinatorial model.

preprint2026arXiv

Skyrme-Hartree-Fock-Bogoliubov mass models on a 3D mesh: V. The N2LO extension of the Skyrme EDF

We present BSkG5, the latest entry in the Brussels-Skyrme-on-a-Grid (BSkG) series and the first large-scale nuclear structure model based on next-to-next-to-leading order (N2LO) Skyrme energy density functional (EDF). By extending the traditional Skyrme EDF ansatz with central terms containing up to four gradients, we are able to combine an excellent global description of nuclear ground state properties with a stiff equation of state for pure neutron matter that is consistent with all astronomical observations of neutron stars. More precisely, the new model matches the accuracy of earlier BSkG models but with two parameters less: we achieve root-mean-square deviations of 0.649 MeV for 2457 atomic masses, 0.0267 fm for 810 charge radii, and 0.43 MeV for 45 primary fission barriers of actinide nuclei. We demonstrate that the complexities of N2LO EDFs are not insurmountable, even for demanding many-body calculations.

preprint2026arXiv

The intermediate neutron capture process. VI. Proton ingestion and i-process in rotating magnetic asymptotic giant branch stars

The intermediate neutron-capture process (i-process) can occur during proton ingestion events (PIEs), which may take place in the early evolutionary phases of asymptotic giant branch (AGB) stars. We investigate the impact of rotational and magnetic mixing on i-process nucleosynthesis in low-metallicity, low-mass AGB stars. We computed AGB models with [Fe/H] = $-2.5$ and $-1.7$ and initial masses of 1 and 1.5 $M_{\odot}$ using the STAREVOL code, including a network of 1160 nuclei coupled to transport equations. Rotating models incorporate a calibrated Tayler-Spruit (TS) dynamo to account for core rotation rates inferred from asteroseismic observations of solar-metallicity sub-giants and giants. Initial rotation velocities of 0, 30, and 90 km s$^{-1}$ were considered, along with varying assumptions for magnetic mixing. We find that rotation without magnetic fields strongly suppresses the i-process due to the production of primary $^{14}$N, which is subsequently converted into $^{22}$Ne $-$ a potent neutron poison during the PIE. Including magnetic fields via the TS dynamo restores the models close to their non-rotating counterparts: strong core-envelope coupling suppresses shear mixing and prevents primary $^{14}$N synthesis, yielding i-process nucleosynthesis similar to non-rotating models. We also find that rotational mixing during the AGB phase is insufficient to affect the occurrence of PIEs. Proton ingestion event-driven nucleosynthesis proceeds similarly in asteroseismic-calibrated magnetic rotating AGB stars and non-rotating stars, producing identical abundance patterns.

preprint2022arXiv

Modelling the spectra of the kilonova AT2017gfo -- I: The photospheric epochs

The kilonova (KN) associated with the binary neutron star (BNS) merger GW170817 is the only known electromagnetic counterpart to a gravitational wave source. Here we produce a sequence of radiative transfer models (using $\textsc{tardis}$) with updated atomic data, and compare them to accurately calibrated spectra. We use element compositions from nuclear network calculations based on a realistic hydrodynamical simulation of a BNS merger. We show that the blue spectrum at +1.4 days after merger requires a nucleosynthetic trajectory with a high electron fraction. Our best-fitting model is composed entirely of first $r$-process peak elements (Sr & Zr) and the strong absorption feature is reproduced well by Sr$\,\textsc{ii}$ absorption. At this epoch, we set an upper limit on the lanthanide mass fraction of $X_{\textsc{ln}} \lesssim 5 \times 10^{-3}$. In contrast, all subsequent spectra from $+2.4 - 6.4$ days require the presence of a modest amount of lanthanide material ($X_{\textsc{ln}} \simeq 0.05^{+0.05}_{-0.02}$), produced by a trajectory with $Y_{\rm e} = 0.29$. This produces lanthanide-induced line blanketing below 6000$\,$Å, and sufficient light $r$-process elements to explain the persistent strong feature at $\sim 0.7 - 1.0 \, \micron$ (Sr$\,\textsc{ii}$). The composition gives good matches to the observed data, indicating that the strong blue flux deficit results in the near-infrared (NIR) excess. The disjoint in composition between the first epoch and all others indicates either ejecta stratification, or the presence of two distinct components of material. This further supports the `two-component' kilonova model, and constrains the element composition from nucleosynthetic trajectories. The major uncertainties lie in availability of atomic data and the ionisation state of the expanding material.

preprint2022arXiv

p-process chaser detector in $n$-$γ$ coincidences

We propose two types of neutron-$γ_1$-$γ_2$ triple coincidence detectors (not constructed) to chase gamma transitions to produce p-nuclei following the neutron emission in the $(γ, n)$ reaction. Neutrons are detected with 24 $^3$He counters embedded in a polyethylene moderator in Type I detector and with 6 liquid scintillation detectors in Type II detector, respectively. $γ$ rays are detected with two high-purity germanium detectors and four LaBr$_3$(Ce) detectors. The detector which is referred to as p-process chaser detector is used to search for mediating states in $^{180}$Ta through which the isomeric and ground states in $^{180}$Ta are thermalized in the p-process. A search is made for both resonant states and unresolved states in high nuclear-level-density domain.

preprint2022arXiv

The intermediate neutron capture process. I. Development of the i-process in low-metallicity low-mass AGB stars

Results from observations report a growing number of metal-poor stars showing an abundance pattern midway between the s- and r-processes. These so-called r/s-stars raise the need for an intermediate neutron capture process (i-process), which is thought to result from the ingestion of protons in a convective helium-burning region, but whose astrophysical site is still largely debated. We investigate whether an i-process during the asymptotic giant branch (AGB) phase of low-metallicity low-mass stars can develop and whether it can explain the abundances of observed r/s-stars. At the beginning of the AGB phase, during the third thermal pulse, the helium driven convection zone is able to penetrate the hydrogen-rich layers. The subsequent proton ingestion leads to a strong neutron burst with neutron densities of $\approx 4.3 \times 10^{14}$ cm$^{-3}$ at the origin of the synthesis of i-process elements. The nuclear energy released by proton burning in the helium-burning convective shell strongly affects the internal structure: the thermal pulse splits and after approximately ten years the upper part of the convection zone merges with the convective envelope. The surface carbon abundance is enhanced by more than 3 dex. This leads to an increase in the opacity, which triggers a strong mass loss and prevents any further thermal pulse. We show that specific isotopic ratios of Ba, Nd, Sm, and Eu can represent good tracers of i-process nucleosynthesis. Finally, an extended comparison with 14 selected r/s-stars show that the observed composition patterns can be well reproduced by our i-process AGB model.

preprint2022arXiv

The p-process in exploding rotating massive stars

The p-process nucleosynthesis can explain proton-rich isotopes that are heavier than iron, which are observed in the Solar System, but discrepancies still persist and important questions concerning the astrophysical site(s) of the p-process remain unanswered. We investigate how the p-process operates in exploding rotating massive stars that have experienced an enhanced s-process nucleosynthesis during their life through rotational mixing. We computed 25 $M_{\odot}$ stellar models at a metallicity of $Z=10^{-3}$ with different initial rotation velocities and rates for the uncertain $^{17}$O($α$,$γ$)$^{21}$Ne reaction. The nucleosynthesis calculation, followed with a network of 737 isotopes, was coupled to stellar evolution, and the p-process nucleosynthesis was calculated in post-processing during both the final evolutionary stages and spherical explosions of various energies. In our models, the p-nuclides are mainly synthesized during the explosion, but not much during the ultimate hydrostatic burning stages. The p-process yields mostly depend on the initial number of trans-iron seeds, which in turn depend on the initial rotation. We found that the impact of rotation on the p-process is comparable to the impact of rotation on the s-process. From no to fast rotation, the s-process yields of nuclides with mass number $A<140$ increase by $3-4$ dex, and so do the p-process yields. Fast rotation with a lower $^{17}$O($α,γ$) rate significantly produces s- and p-nuclides with $A\geq140$. Our results suggest that the contribution of core-collapse supernovae from massive stars to the solar (and Galactic) p-nuclei has been underestimated in the past, and more specifically, that the contribution from massive stars with sub-solar metallicities may even dominate. A more detailed study including stellar models with a wide range of masses and metallicities remains to be performed.

preprint2021arXiv

Fission fragment distributions and their impact on the r-process nucleosynthesis in neutron star mergers

Neutron star (NS) merger ejecta offer a viable site for the production of heavy r-process elements with nuclear mass numbers A >140. The crucial role of fission recycling is responsible for the robustness of this site against many astrophysical uncertainties. Here, we introduce new improvements to our scission-point model, called SPY, to derive the fission fragment distribution for all neutron-rich fissioning nuclei of relevance in r-process calculations. These improvements include a phenomenological modification of the scission distance and a smoothing procedure of the distribution. Such corrections lead to a much better agreement with experimental fission yields. Those yields are also used to estimate the number of neutrons emitted by the excited fragments on the basis of different neutron evaporation models. Our new fission yields are extensively compared to those predicted by the so-called GEF model. The impact of fission on the r-process nucleosynthesis in binary neutron mergers is also reanalyzed. Two scenarios are considered, the first one with low initial electron fraction is subject to intense fission recycling, in contrast to the second one which includes weak interactions on nucleons. The various regions of the nuclear chart responsible for fission recycling during the neutron irradiation as well as after freeze-out are discussed. The contribution fission processes may have to the final abundance distribution is also studied in detail in the light of newly defined quantitative indicators describing the fission recycling, the fission seeds and the fission progenitors. In particular, those allow us to estimate the contribution of fission to the final abundance distribution stemming from specific heavy nuclei. Calculations obtained with SPY and GEF fission fragment distributions are compared for both r-process scenarios.

preprint2021arXiv

Low-mass low-metallicity AGB stars as an efficient i-process site explaining CEMP-rs stars

Among Carbon-Enhanced Metal-Poor (CEMP) stars, some are found to be enriched in s-process elements (CEMP-s), in r-process elements (CEMP-r) or in both s- and r-process elements (CEMP-rs). The origin of the abundance differences between CEMP-s and CEMP-rs stars is presently unknown. It has been claimed that the i-process, whose site still remains to be identified, could better reproduce CEMP-rs abundances than the s-process. We analyze high-resolution spectra of 25 metal-poor stars, observed with the high-resolution HERMES spectrograph mounted on the Mercator telescope, La Palma, or with the UVES/VLT and HIRES/KECK spectrographs. We propose a new, robust classification method for CEMP-s and CEMP-rs stars using eight heavy element abundances. The abundance profiles of CEMP-s and CEMP-rs stars are derived and there appears to be an abundance continuum between the two stellar classes. CEMP-rs stars present most of the characteristics of extrinsic stars such as CEMP-s, CH, Barium and extrinsic S stars, with an even larger binarity rate among CEMP-rs stars than among CEMP-s stars. Stellar evolutionary tracks of an enhanced carbon composition (consistent with our abundance determinations) are necessary to explain the position of CEMP-s and CEMP-rs stars in the HR diagram using Gaia DR2 parallaxes; they are found to lie mostly on the RGB. CEMP-rs stars can be explained as being polluted by a low-mass, low-metallicity TP-AGB companion experiencing i-process nucleosynthesis after proton ingestion during its first convective thermal pulses. The global fitting of our i-process models to CEMP-rs stars is as good as the one of our s-process models to CEMP-s stars. As such, CEMP-rs stars could be renamed as CEMP-sr stars, since they represent a particular manifestation of the s-process at low-metallicities. For these objects a call for an exotic i-process site may not necessarily be required anymore.

preprint2021arXiv

Statistical properties of the well deformed $^{153,155}$Sm nuclei and the scissors resonance

The Nuclear Level Densities (NLDs) and the $γ$-ray Strength Functions ($γ$SFs) of $^{153,155}$Sm have been extracted from (d,p$γ$) coincidences using the Oslo method. The experimental NLD of $^{153}$Sm is higher than the NLD of $^{155}$Sm, in accordance with microscopic calculations. The $γ$SFs of $^{153,155}$Sm are in fair agreement with QRPA calculations based on the D1M Gogny interaction. An enhancement is observed in the $γ$SF for both $^{153,155}$Sm nuclei around 3 MeV in excitation energy and is attributed to the M1 Scissors Resonance (SR). Their integrated strengths were found to be in the range 1.3 - 2.1 and 4.4 - 6.4 $μ^{2}_{N}$ for $^{153}$Sm and $^{155}$Sm, respectively. The strength of the SR for $^{155}$Sm is comparable to those for deformed even-even Sm isotopes from nuclear resonance fluorescence measurements, while that of $^{153}$Sm is lower than expected.

preprint2020arXiv

Astronuclear Physics: a Tale of the Atomic Nuclei in the Skies

A century ago, nuclear physics entered astrophysics, giving birth to a new field of science referred to as "Nuclear Astrophysics". With time, it developed at an impressive pace into a vastly inter- and multidisciplinary discipline bringing into its wake not only astronomy and cosmology, but also many other sub-fields of physics, especially particle, solid-state and computational physics, as well as chemistry, geology and even biology. The present Astronuclear Physics review focusses primarily on the facets of nuclear physics that are of relevance to astronomy and astrophysics, the theoretical aspects being of special concern here.

preprint2020arXiv

Discovery of technetium- and niobium-rich S stars: the case for bitrinsic stars

S stars are late-type giants with overabundances of s-process elements. They come in two flavours depending on the presence or not of technetium (Tc), an element without stable isotopes. Intrinsic S stars are Tc-rich and genuine asymptotic giant branch (AGB)stars while extrinsic S stars owe their s-process overabundances to the pollution from a former AGB companion, now a white dwarf(WD). In addition to Tc, another distinctive feature between intrinsic and extrinsic S stars is the overabundance of niobium (Nb) in the latter class. We discuss the case of the S stars BD+79 156 and o1Ori whose specificity is to share the distinctive features of both intrinsic and extrinsic S stars, namely the presence of Tc along with a Nb overabundance. BD+79 156 is the first clear case of a bitrinsic star, i.e. a doubly s-process-enriched object, first through mass transfer in a binary system, and then through internal nucleosynthesis (responsible for the Tc-enrichment in BD+79 156 which must therefore have reached the AGB phase of its evolution). This hybrid nature of the s-process pattern in BD+79 156 is supported by its binary nature and its location in the HR diagram just beyond the onset of the third dredge-up on the AGB. The Tc-rich, binary S-star o1Ori with a WD companion was another long-standing candidate for a similar hybrid s-process enrichment. However the marginal overabundance of Nb derived ino1Ori does not allow to trace unambiguously the evidence of a large pollution coming from the AGB progenitor of its current WD companion. As a side product, the current study offers a new way of detecting binary AGB stars with WD companions by identifying their Tc-rich nature along with a Nb overabundance.

preprint2019arXiv

Unified equations of state for cold non-accreting neutron stars with Brussels-Montreal functionals. I. Role of symmetry energy

The theory of the nuclear energy-density functional is used to provide a unified and thermodynamically consistent treatment of all regions of cold non-accreting neutron stars. In order to assess the impact of our lack of complete knowledge of the density dependence of the symmetry energy on the constitution and the global structure of neutron stars, we employ four different functionals. All of them were precision fitted to essentially all the nuclear-mass data with the Hartree-Fock-Bogoliubov method and two different neutron-matter equations of state based on realistic nuclear forces. For each functional, we calculate the composition, the pressure-density relation, and the chemical potentials throughout the star. We show that uncertainties in the symmetry energy can significantly affect the theoretical results for the composition and global structure of neutron stars. To facilitate astrophysical applications, we construct analytic fits to our numerical results.

S. Goriely

Quick read

Decide how to stay connected

How to connect with this researcher

Open a focused conversation when the fit is right

See the researcher in context

Building this graph slice

13 published item(s)

Improved microscopic nuclear level densities within the triaxial Hartree-Fock-Bogoliubov plus combinatorial method

Skyrme-Hartree-Fock-Bogoliubov mass models on a 3D mesh: V. The N2LO extension of the Skyrme EDF

The intermediate neutron capture process. VI. Proton ingestion and i-process in rotating magnetic asymptotic giant branch stars

Modelling the spectra of the kilonova AT2017gfo -- I: The photospheric epochs

p-process chaser detector in $n$-$γ$ coincidences

The intermediate neutron capture process. I. Development of the i-process in low-metallicity low-mass AGB stars

The p-process in exploding rotating massive stars

Fission fragment distributions and their impact on the r-process nucleosynthesis in neutron star mergers

Low-mass low-metallicity AGB stars as an efficient i-process site explaining CEMP-rs stars

Statistical properties of the well deformed $^{153,155}$Sm nuclei and the scissors resonance

Astronuclear Physics: a Tale of the Atomic Nuclei in the Skies

Discovery of technetium- and niobium-rich S stars: the case for bitrinsic stars

Unified equations of state for cold non-accreting neutron stars with Brussels-Montreal functionals. I. Role of symmetry energy