Showing 1–50 of 178 results for author: Goriely, S

Photoneutron cross section measurements on $^{208}$Pb in the Giant Dipole Resonance region

Authors: I. Gheorghe, S. Goriely, N. Wagner, T. Aumann, M. Baumann, P. van Beek, P. Kuchenbrod, H. Scheit, D. Symochko, T. Ari-izumi, F. L. Bello Garrote, T. Eriksen, W. Paulsen, L. G. Pedersen, F. Reaz, V. W. Ingeberg, S. Belyshev, S. Miyamoto, H. Utsunomiya

Abstract: Photoneutron reactions on $^{208}$Pb in the Giant Dipole Resonance energy region have been investigated at the $γ$-ray beam line of the NewSUBARU facility in Japan. The measurements made use of quasi-monochromatic laser Compton backscattering $γ$-ray beams in a broad energy range, from the neutron threshold up to 38 MeV, and of a flat-efficiency moderated $^3$He neutron detection system along with… ▽ More Photoneutron reactions on $^{208}$Pb in the Giant Dipole Resonance energy region have been investigated at the $γ$-ray beam line of the NewSUBARU facility in Japan. The measurements made use of quasi-monochromatic laser Compton backscattering $γ$-ray beams in a broad energy range, from the neutron threshold up to 38 MeV, and of a flat-efficiency moderated $^3$He neutron detection system along with associated neutron-multiplicity sorting methods. We report absolute cross sections and mean photoneutron energies for the $^{208}$Pb$(γ,\,inX)$ reactions with $i$~=~1 to 4. The fine structure present in the $^{208}$Pb$(γ,\,n)$ cross sections at incident energies lower than 13~MeV has been observed. The photoabsorption cross section has been obtained as the sum of the $(γ,\,inX)$ reaction cross sections. By reproducing the measured ring-ratio values at excitation energies below the two neutron separation energy, we were able to extract estimations on the $^{208}$Pb$(γ,\,n)$ photoneutron energy spectra and on the partial photoneutron cross sections for leaving the residual $^{207}$Pb in its ground and first two excited states. The present results are compared with data from the literature and statistical model calculations. △ Less

Submitted 18 March, 2024; originally announced March 2024.

Giant Dipole Resonance photofission and photoneutron reactions in $^{238}$U and $^{232}$Th

Authors: D. Filipescu, I. Gheorghe, S. Goriely, A. Tudora, K. Nishio, T. Ohtsuki, H. Wang, G. Fan, K. Stopani, F. Suzaki, K. Hirose, M. Inagaki, Y. -W. Lui, T. Ari-izumi, S. Miyamoto, T. Otsuka, H. Utsunomiya

Abstract: New measurements of photofission and photoneutron reactions on $^{238}$U and $^{232}$Th in the Giant Dipole Resonance (GDR) energy region have been performed at the laser Compton-scattering $γ$-ray source of the NewSUBARU synchrotron radiation facility using a high-and-flat efficiency moderated $^3$He detection array. The neutron-multiplicity sorting of high-multiplicity fission neutron coincidenc… ▽ More New measurements of photofission and photoneutron reactions on $^{238}$U and $^{232}$Th in the Giant Dipole Resonance (GDR) energy region have been performed at the laser Compton-scattering $γ$-ray source of the NewSUBARU synchrotron radiation facility using a high-and-flat efficiency moderated $^3$He detection array. The neutron-multiplicity sorting of high-multiplicity fission neutron coincidence events has been performed using a dedicated energy dependent, multiple firing statistical treatment. The photoneutron $(γ,\,in)$ with $i$~=~1~--~3 and photofission $(γ,\,F)$ reactions have been discriminated by considering a Gaussian distribution of prompt-fission-neutron (PFN) multiplicities predicted by the theory of evaporation in sequential neutron emission from excited fission fragments. We report experimental $(γ,\,n)$, $(γ,\,2n)$, $(γ,\,3n)$ and $(γ,\,F)$ cross sections, average energies of PFNs and of $(γ,\,in)$ photoneutrons, as well as the mean number of PFNs per fission act and the width of the PFNs multiplicity distribution. Based on these primary experimental results and combined with reasonable assumptions, we extract also the first- and second-chance fission contributions. The new experimental results are compared with statistical-model calculations performed with the EMPIRE-3.2 Malta and TALYS-1.964 codes on the present data and with prompt fission emission calculations obtained with the Los Alamos model in the frame of the most probable fragmentation approach with and without sequential emission. △ Less

Submitted 20 February, 2024; originally announced February 2024.

The intermediate neutron capture process. V. The i-process in AGB stars with overshoot

Authors: A. Choplin, L. Siess, S. Goriely, S. Martinet

Abstract: The intermediate neutron capture process (i-process) can develop during proton ingestion events (PIE), potentially during the early stages of low-mass low-metallicity asymptotic giant branch (AGB) stars. We examine the impact of overshoot mixing on the triggering and development of i-process nucleosynthesis in AGB stars of various initial masses and metallicities. We computed AGB stellar models, w… ▽ More The intermediate neutron capture process (i-process) can develop during proton ingestion events (PIE), potentially during the early stages of low-mass low-metallicity asymptotic giant branch (AGB) stars. We examine the impact of overshoot mixing on the triggering and development of i-process nucleosynthesis in AGB stars of various initial masses and metallicities. We computed AGB stellar models, with initial masses of 1, 2, 3, and 4 M$_{\odot}$ and metallicities in the $-2.5 \le $ [Fe/H] $\le 0$ range, using the stellar evolution code STAREVOL with a network of 1160 nuclei coupled to the transport equations. We considered different overshooting profiles below and above the thermal pulses, and below the convective envelope. The occurrence of PIEs is found to be primarily governed by the amount of overshooting at the top of pulse ($f_{\rm top}$) and to increase with rising $f_{\rm top}$. For $f_{\rm top} =$ 0, 0.02, 0.04, and 0.1, we find that 0 %, 6 %, 24 %, and 86 % of our 21 AGB models with $-2<$ [Fe/H] $<0$ experience a PIE, respectively. We also find that PIEs leave a $^{13}$C-pocket at the bottom of the pulse that can give rise to an additional radiative s-process nucleosynthesis, and ultimately produce a noticeable mixed i+s chemical signature at the surface. Finally, the chemical abundance patterns of 22 observed r/s-stars candidates with $-2<$ [Fe/H] $<-1$ are found to be in reasonable agreement with our AGB model predictions. The binary status of the dwarfs/giants being unclear, we suggest that these stars have acquired their chemical pattern either from the mass transfer of a now-extinct AGB companion or from an early generation AGB star that polluted the natal cloud. Stricter constraints from multi-dimensional hydrodynamical models on overshoot coefficients could deliver new insights into the contribution of AGB stars to heavy elements in the Universe. △ Less

Submitted 21 February, 2024; v1 submitted 15 February, 2024; originally announced February 2024.

Comments: 21 pages, 15 figures, accepted in A&A

Systematic study of the low-lying electric dipole strength in Sn isotopes and its astrophysical implications

Authors: M. Markova, A. C. Larsen, P. von Neumann-Cosel, E. Litvinova, A. Choplin, S. Goriely, S. Martinet, L. Siess, M. Guttormsen, F. Pogliano, S. Siem

Abstract: The $γ$-ray strength functions (GSF) and nuclear level densities (NLD) below the neutron threshold have been extracted for $^{111-113,116-122,124}$Sn from particle-$γ$ coincidence data with the Oslo method. The evolution of bulk properties of the low-lying electric dipole response has been investigated on the basis of the Oslo GSF data and results of a recent systematic study of electric and magne… ▽ More The $γ$-ray strength functions (GSF) and nuclear level densities (NLD) below the neutron threshold have been extracted for $^{111-113,116-122,124}$Sn from particle-$γ$ coincidence data with the Oslo method. The evolution of bulk properties of the low-lying electric dipole response has been investigated on the basis of the Oslo GSF data and results of a recent systematic study of electric and magnetic dipole strengths in even-even Sn isotopes with relativistic Coulomb excitation. The obtained GSFs reveal a resonance-like peak on top of the tail of the isovector giant dipole resonance, centered at $\approx$8 MeV and exhausting $\approx$2\% of the classical Thomas-Reiche-Kuhn (TRK) sum. In contrast to predictions of the relativistic quasiparticle random-phase and time-blocking approximation calculations (RQRPA and RQTBA), no monotonous increase in the total low-lying $E1$ strength was observed in the experimental data from $^{111}$Sn to $^{124}$Sn, demonstrating rather similar strength distributions in these nuclei. The Oslo GSFs and NLDs were further used as inputs to constrain the cross sections and Maxwellian-averaged cross sections of $(n,γ)$ reactions in the Sn isotopic chain using TALYS. The obtained results agree well with other available experimental data and the recommended values from the JINA REACLIB, BRUSLIB, and KADoNiS libraries. Despite relatively small exhausted fractions of the TRK sum rule, the low-lying electric dipole strength makes a noticeable impact on the radiative neutron-capture cross sections in stable Sn isotopes. Moreover, the experimental Oslo inputs for the $^{121,123}$Sn$(n,γ)$$^{122,124}$Sn reactions were found to affect the production of Sb in the astrophysical $i$-process, providing new constraints on the uncertainties of the resulting chemical abundances from multi-zone low-metallicity Asymptotic Giant Branch stellar models. △ Less

Submitted 15 November, 2023; originally announced November 2023.

Comments: 27 pages, 14 pages. Submitted to Physical Review C journal on 13 November 2023

MELCHIORS: The Mercator Library of High Resolution Stellar Spectroscopy

Authors: P. Royer, T. Merle, K. Dsilva, S. Sekaran, H. Van Winckel, Y. Frémat, M. Van der Swaelmen, S. Gebruers, A. Tkachenko, M. Laverick, M. Dirickx, G. Raskin, H. Hensberge, M. Abdul-Masih, B. Acke, M. L. Alonso, S. Bandhu Mahato, P. G. Beck, N. Behara, S. Bloemen, B. Buysschaert, N. Cox, J. Debosscher, P. De Cat, P. Degroote , et al. (49 additional authors not shown)

Abstract: Over the past decades, libraries of stellar spectra have been used in a large variety of science cases, including as sources of reference spectra for a given object or a given spectral type. Despite the existence of large libraries and the increasing number of projects of large-scale spectral surveys, there is to date only one very high-resolution spectral library offering spectra from a few hundr… ▽ More Over the past decades, libraries of stellar spectra have been used in a large variety of science cases, including as sources of reference spectra for a given object or a given spectral type. Despite the existence of large libraries and the increasing number of projects of large-scale spectral surveys, there is to date only one very high-resolution spectral library offering spectra from a few hundred objects from the southern hemisphere (UVES-POP) . We aim to extend the sample, offering a finer coverage of effective temperatures and surface gravity with a uniform collection of spectra obtained in the northern hemisphere. Between 2010 and 2020, we acquired several thousand echelle spectra of bright stars with the Mercator-HERMES spectrograph located in the Roque de Los Muchachos Observatory in La Palma, whose pipeline offers high-quality data reduction products. We have also developed methods to correct for the instrumental response in order to approach the true shape of the spectral continuum. Additionally, we have devised a normalisation process to provide a homogeneous normalisation of the full spectral range for most of the objects. We present a new spectral library consisting of 3256 spectra covering 2043 stars. It combines high signal-to-noise and high spectral resolution over the entire range of effective temperatures and luminosity classes. The spectra are presented in four versions: raw, corrected from the instrumental response, with and without correction from the atmospheric molecular absorption, and normalised (including the telluric correction). △ Less

Submitted 5 November, 2023; originally announced November 2023.

Comments: 17 pages, 18 figures Preview and access to the library: https://www.royer.se/melchiors.html

Journal ref: A&A 681, A107 (2024)

The intermediate neutron capture process: IV. Impact of nuclear model and parameter uncertainties

Authors: S. Martinet, A. Choplin, S. Goriely, L. Siess

Abstract: We investigate both the systematic and statistical uncertainties associated with theoretical nuclear reaction rates of relevance during the i-process and explore their impact on the i-process elemental production, and subsequently on the surface enrichment, for a low-mass low-metallicity star during the early AGB phase. We use the TALYS reaction code (Koning et al. 2023) to estimate both the model… ▽ More We investigate both the systematic and statistical uncertainties associated with theoretical nuclear reaction rates of relevance during the i-process and explore their impact on the i-process elemental production, and subsequently on the surface enrichment, for a low-mass low-metallicity star during the early AGB phase. We use the TALYS reaction code (Koning et al. 2023) to estimate both the model and parameter uncertainties affecting the photon strength function and the nuclear level densities, hence the radiative neutron capture rates. The STAREVOL code (Siess et al. 2006) is used to determine the impact of nuclear uncertainties on the i-process nucleosynthesis in a 1 $M_{\odot}$ [Fe/H] = - 2.5 model star during the proton ingestion event in the early AGB phase. A large nuclear network of 1160 species coherently coupled to the transport processes is solved to follow the i-process nucleosynthesis. We find that the non-correlated parameter uncertainties lead the surface abundances uncertainties of element with $Z\geq 40$ to range between 0.5 and 1.0 dex, with odd-$Z$ elements displaying higher uncertainties. The correlated model uncertainties are of the same order of magnitude, and both model and parameter uncertainties have an important impact on potential observable tracers such as Eu and La. Both the correlated model and uncorrelated parameter uncertainties need to be estimated coherently before being propagated to astrophysical observables through multi-zone stellar evolution models. Many reactions are found to affect the i-process predictions and will require improved nuclear models guided by experimental constraints. Priority should be given to the reactions influencing the observable tracers. △ Less

Submitted 12 October, 2023; originally announced October 2023.

Comments: Accepted: October 11, 2023 \\ 14 Pages, 14 Figures, 2 Tables

Skyrme-Hartree-Fock-Bogoliubov mass models on a 3D mesh: III. From atomic nuclei to neutron stars

Authors: Guilherme Grams, Wouter Ryssens, Guillaume Scamps, Stephane Goriely, Nicolas Chamel

Abstract: We present BSkG3, the latest entry in the Brussels-Skyrme-on-a-grid series of large-scale models of nuclear structure based on an energy density functional. Compared to its predecessors, the new model offers a more realistic description of nucleonic matter at the extreme densities relevant to neutron stars. This achievement is made possible by incorporating a constraint on the infinite nuclear mat… ▽ More We present BSkG3, the latest entry in the Brussels-Skyrme-on-a-grid series of large-scale models of nuclear structure based on an energy density functional. Compared to its predecessors, the new model offers a more realistic description of nucleonic matter at the extreme densities relevant to neutron stars. This achievement is made possible by incorporating a constraint on the infinite nuclear matter properties at high densities in the parameter adjustment, ensuring in this way that the predictions of BSkG3 for the nuclear Equation of State are compatible with the observational evidence for heavy pulsars with $M > 2 M_{\odot}$. Instead of the usual phenomenological pairing terms, we also employ a more microscopically founded treatment of nucleon pairing, resulting in extrapolations to high densities that are in line with the predictions of advanced many-body methods and are hence more suited to the study of superfluidity in neutron stars. By adopting an extended form of the Skyrme functional, we are able to reconcile the description of matter at high densities and at saturation density: the new model further refines the description of atomic nuclei offered by its predecessors. A qualitative improvement is our inclusion of ground state reflection asymmetry, in addition to the spontaneous breaking of rotational, axial, and time-reversal symmetry. Quantitatively, the model offers lowered root-mean-square deviations on 2457 masses (0.631 MeV), 810 charge radii (0.0237 fm) and an unmatched accuracy with respect to 45 primary fission barriers of actinide nuclei (0.33 MeV). Reconciling the complexity of neutron stars with those of atomic nuclei establishes BSkG3 as a tool of choice for applications to nuclear structure, the nuclear equation of state and nuclear astrophysics in general. △ Less

Submitted 6 October, 2023; v1 submitted 26 July, 2023; originally announced July 2023.

Report number: NT@UW-23-11

Nuclear Level Density and $γ$-ray Strength Function of $^{67}\mathrm{Ni}$ and the impact on the i-process

Authors: V. W. Ingeberg, S. Siem, M. Wiedeking, A. Choplin, S. Goriely, L. Siess, K. J. Abrahams, K. Arnswald, F. Bello Garrote, D. L. Bleuel, J. Cederkäll, T. L. Christoffersen, D. M. Cox, H. De Witte, L. P. Gaffney, A. Görgen, C. Henrich, A. Illana, P. Jones, B. V. Kheswa, T. Kröll, S. N. T. Majola, K. L. Malatji, J. Ojala, J. Pakarinen , et al. (7 additional authors not shown)

Abstract: Proton-$γ$ coincidences from $(\mathrm{d},\mathrm{p})$ reactions between a $^{66}\mathrm{Ni}$ beam and a deuterated polyethylene target have been analyzed with the inverse-Oslo method to find the nuclear level density (NLD) and $γ$-ray strength function ($γ$SF) of $^{67}\mathrm{Ni}$. The $^{66}\mathrm{Ni}(\mathrm{n},γ)$ capture cross section has been calculated using the Hauser-Feshbach model in T… ▽ More Proton-$γ$ coincidences from $(\mathrm{d},\mathrm{p})$ reactions between a $^{66}\mathrm{Ni}$ beam and a deuterated polyethylene target have been analyzed with the inverse-Oslo method to find the nuclear level density (NLD) and $γ$-ray strength function ($γ$SF) of $^{67}\mathrm{Ni}$. The $^{66}\mathrm{Ni}(\mathrm{n},γ)$ capture cross section has been calculated using the Hauser-Feshbach model in TALYS using the measured NLD and $γ$SF as constraints. The results confirm that the $^{66}\mathrm{Ni}(\mathrm{n},γ)$ reaction acts as a bottleneck when relying on one-zone nucleosynthesis calculations. However, the impact of this reaction is strongly damped in multi-zone low-metallicity AGB stellar models experiencing i-process nucleosynthesis. △ Less

Submitted 10 June, 2024; v1 submitted 14 July, 2023; originally announced July 2023.

Comments: Submitted to Phys. Rev. C

Modelling the spectra of the kilonova AT2017gfo -- II: Beyond the photospheric epochs

Authors: J. H. Gillanders, S. A. Sim, S. J. Smartt, S. Goriely, A. Bauswein

Abstract: Binary neutron star mergers are the first confirmed site of element nucleosynthesis by the rapid neutron-capture process (r-process). The kilonova AT2017gfo is the only electromagnetic counterpart of a neutron star merger spectroscopically observed. We analyse the entire spectral sequence of AT2017gfo (from merger to +10.4 days) and identify seven emission-like features. We confirm that the promin… ▽ More Binary neutron star mergers are the first confirmed site of element nucleosynthesis by the rapid neutron-capture process (r-process). The kilonova AT2017gfo is the only electromagnetic counterpart of a neutron star merger spectroscopically observed. We analyse the entire spectral sequence of AT2017gfo (from merger to +10.4 days) and identify seven emission-like features. We confirm that the prominent 1.08 um feature can be explained by the Sr II near-infrared triplet evolving from a P-Cygni profile through to pure emission. We calculate the expected strength of the [Sr II] doublet and show that its absence requires highly clumped ejecta. Near-infrared features at 1.58 and 2.07 um emerge after three days and become more prominent as the spectra evolve. We model these as optically thick P-Cygni profiles and alternatively as pure emission features (with FWHM = 35600 +/- 6600 km/s), and favour the latter interpretation. The profile of the strong 2.07 um emission feature is best reproduced with two lines, centred at 2.059 and 2.135 um. We search for candidate ions for all prominent features in the spectra. Strong, permitted transitions of La III, Ce III, Gd III, Ra II and Ac I are plausible candidates for the emission features. If any of these features are produced by intrinsically weak, forbidden transitions, we highlight candidate ions spanning the three r-process peaks. The second r-process peak elements Te and I have plausible matches to multiple features. We highlight the need for more detailed and quantitative atomic line transition data. △ Less

Submitted 27 November, 2023; v1 submitted 26 June, 2023; originally announced June 2023.

Comments: Updating to the final accepted version in MNRAS. Main text: 20 pages, with 7 figures and 6 tables. Appendices: 9 pages with 3 tables

Does the i-process operate at nearly solar metallicity?

Authors: D. Karinkuzhi, S. Van Eck, S. Goriely, L. Siess, A. Jorissen, A. Choplin, A. Escorza, S. Shetye, H. Van Winckel, .

Abstract: A sample of 895 s-process-rich candidates has been found among the 454180 giant stars surveyed by LAMOST at low spectral resolution (R~1800). In a previous study, taking advantage of the higher resolution (R~86 000) offered by the the HERMES-Mercator spectrograph, we performed the re-analysis of 15 among the brightest stars of this sample. Among these 15 program stars, having close-to-solar metall… ▽ More A sample of 895 s-process-rich candidates has been found among the 454180 giant stars surveyed by LAMOST at low spectral resolution (R~1800). In a previous study, taking advantage of the higher resolution (R~86 000) offered by the the HERMES-Mercator spectrograph, we performed the re-analysis of 15 among the brightest stars of this sample. Among these 15 program stars, having close-to-solar metallicities, 11 showed mild to strong heavy element overabundances. The nucleosynthesis process(es) at the origin of these overabundances were however not questioned in our former study. We derive the abundances in s- and r-process elements of the 15 targets in order to investigate whether some stars also show an i-process signature, as sometimes found in their lower metallicity counterparts (namely, the Carbon-Enhanced Metal-Poor (CEMP)-rs stars). Abundances are derived from the high-resolution HERMES spectra for Pr, Nd, Sm, and Eu, using the TURBOSPECTRUM radiative transfer LTE code with MARCS model atmospheres. Using the new classification scheme proposed in our recent study we find that two stars show overabundances in both s- and r-process elements well above the level expected from the Galactic chemical evolution, an analogous situation to the one of CEMP-rs stars at lower metallicities. We compare the abundances of the most enriched stars with the nucleosynthetic predictions from the STAREVOL stellar evolutionary code and find abundances compatible with an i-process occurring in AGB stars. Despite a larger number of heavy elements to characterize the enrichment pattern, the limit between CEMP-s and CEMP-rs stars remains fuzzy. It is however interesting to note that an increasing number of extrinsic stars are found to have abundances better reproduced by an i-process pattern even at close-to-solar metallicities. △ Less

Submitted 7 May, 2023; originally announced May 2023.

Comments: Accepted for publication in A&A, 9 pages, 9 figures including the two in appendix

Journal ref: A&A 677, A47 (2023)

Experimentally constrained $^{165,166}\text{Ho}(n,γ)$ rates and implications for the $s$ process

Authors: Francesco Pogliano, Ann-Cecilie Larsen, Stephane Goriely, Lionel Siess, Maria Markova, Andreas Görgen, Johannes Heines, Vetle Werner Ingeberg, Robin Grongstad Kjus, Johan Emil Linnestad Larsson, Kevin Ching Wei Li, Elise Malmer Martinsen, Gerard Jordan Owens-Fryar, Line Gaard Pedersen, Gulla Serville Torvund, Artemis Tsantiri

Abstract: The $γ$-ray strength function and the nuclear level density of $^{167}$Ho have been extracted using the Oslo method from a $^{164}\text{Dy}(α,pγ)^{167}$Ho experiment carried out at the Oslo Cyclotron Laboratory. The level density displays a shape that is compatible with %can be approximated with the constant temperature model in the quasicontinuum, while the strength function shows structures indi… ▽ More The $γ$-ray strength function and the nuclear level density of $^{167}$Ho have been extracted using the Oslo method from a $^{164}\text{Dy}(α,pγ)^{167}$Ho experiment carried out at the Oslo Cyclotron Laboratory. The level density displays a shape that is compatible with %can be approximated with the constant temperature model in the quasicontinuum, while the strength function shows structures indicating the presence of both a scissors and a pygmy dipole resonance. Using our present results as well as data from a previous $^{163}\text{Dy}(α,pγ)^{166}$Ho experiment, the $^{165}\text{Ho}(n,γ)$ and $^{166}\text{Ho}(n,γ)$ MACS uncertainties have been constrained. The possible influence of the low-lying, long-lived 6~keV isomer $^{166}$Ho in the $s$ process is investigated in the context of a 2~$M_\odot$, [Fe/H]=-0.5 AGB star. We show that the newly obtained $^{165}\text{Ho}(n,γ)$ MACS affects the final $^{165}$Ho abundance, while the $^{166}\text{Ho}(n,γ)$ MACS only impacts the enrichment of $^{166,167}$Er to a limited degree due to the relatively rapid $β$ decay of the thermalized $^{166}$Ho at typical $s$-process temperatures. △ Less

Submitted 9 June, 2023; v1 submitted 27 April, 2023; originally announced April 2023.

Comments: 11 pages, submitted to Physical Reviews C

Photoneutron emission cross sections for $^{13}$C

Authors: H. Utsunomiya, S. Goriely, M. Kimura, N. Shimizu, Y. Utsuno, G. M. Tveten, T. Renstrøm, T. Ari-izumi, S. Miyamoto

Abstract: Photoneutron emission cross sections were measured for $^{13}$C below $2n$ threshold using quasi-monochromatic $γ$-ray beams produced in laser Compton-scattering at the NewSUBARU synchrotron radiation facility. The data show fine structures in the low-energy tail of the giant-dipole resonance; the integrated strength of the fine structure below 18~MeV is intermediate among the past measurements wi… ▽ More Photoneutron emission cross sections were measured for $^{13}$C below $2n$ threshold using quasi-monochromatic $γ$-ray beams produced in laser Compton-scattering at the NewSUBARU synchrotron radiation facility. The data show fine structures in the low-energy tail of the giant-dipole resonance; the integrated strength of the fine structure below 18~MeV is intermediate among the past measurements with bremsstrahlung and the positron annihilation $γ$ rays. We compare the photoneutron emission data with the {\sf TALYS} statistical model calculation implemented with the simple modified Lorentzian model of $E1$ and $M1$ strengths. We also compare the total photoabsorption cross sections for $^{13}$C with the shell model and antisymmetrized molecular dynamics calculations as well as the statistical model calculation. We further investigate the consistency between the present photoneutron emission and the reverse $^{12}$C(n,$γ$) cross sections through their corresponding astrophysical rate. △ Less

Submitted 7 March, 2023; originally announced March 2023.

End-to-end kilonova models of neutron-star mergers with delayed black-hole formation

Authors: Oliver Just, Vimal Vijayan, Zewei Xiong, Stephane Goriely, Theodoros Soultanis, Andreas Bauswein, Jérôme Guilet, Hans-Thomas Janka, Gabriel Martínez-Pinedo

Abstract: We investigate the nucleosynthesis and kilonova properties of binary neutron-star (NS) merger models which lead to intermediate remnant lifetimes of ~0.1-1seconds until black-hole (BH) formation and describe all components of material ejected during the dynamical merger phase, NS-remnant evolution, and final viscous disintegration of the BH torus after gravitational collapse. To this end we employ… ▽ More We investigate the nucleosynthesis and kilonova properties of binary neutron-star (NS) merger models which lead to intermediate remnant lifetimes of ~0.1-1seconds until black-hole (BH) formation and describe all components of material ejected during the dynamical merger phase, NS-remnant evolution, and final viscous disintegration of the BH torus after gravitational collapse. To this end we employ a combination of hydrodynamics, nucleosynthesis, and radiative-transfer tools to achieve a consistent end-to-end modeling of the system and its observables. We adopt a novel version of the Shakura-Sunyaev scheme allowing to vary the approximate turbulent viscosity inside the NS remnant independently of the surrounding disk. We find that asymmetric progenitors lead to shorter remnant lifetimes and enhanced ejecta masses, although the viscosity affects the absolute values of these characteristics. The integrated production of lanthanides and heavier elements in such binary systems is sub-solar, suggesting that the considered scenarios contribute in a sub-dominant fashion to r-process enrichment. One reason is that BH-tori formed after delayed collapse exhibit less neutron-rich conditions than typically found, and often assumed in previous BH-torus models, for early BH formation. The outflows in our models feature strong anisotropy as a result of the lanthanide-poor polar neutrino-driven wind pushing aside lanthanide-rich dynamical ejecta. Considering the complexity of the models, the estimated kilonova light curves show promising agreement with AT2017gfo after times of several days, while the remaining inconsistencies at early times could possibly be overcome in binary configurations with a more dominant neutrino-driven wind relative to the dynamical ejecta. △ Less

Submitted 2 June, 2023; v1 submitted 21 February, 2023; originally announced February 2023.

Comments: 16 pages, 9 figures, 1 table, accepted to ApJL

Skyrme-Hartree-Fock-Bogoliubov mass models on a 3D mesh. IIb. Fission properties of BSkG2

Authors: Wouter Ryssens, Guillaume Scamps, Stephane Goriely, Michael Bender

Abstract: 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 coo… ▽ More 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. △ Less

Submitted 8 May, 2023; v1 submitted 6 February, 2023; originally announced February 2023.

Comments: Version accepted for publication; 25 pages, 11 figures, 3 tables

Report number: NT@UW-23-01

Journal ref: The European Physical Journal A 59, 96 (2023)

Opacities of Singly and Doubly Ionised Neodymium and Uranium for Kilonova Emission Modeling

Authors: A. Flörs, R. F. Silva, J. Deprince, H. Carvajal Gallego, G. Leck, G. Martínez-Pinedo, J. M. Sampaio, P. Amaro, J. P. Marques, S. Goriely, P. Quinet, P. Palmeri, M. Godefroid

Abstract: Even though the electromagnetic counterpart AT2017gfo to the binary neutron star merger GW170817 is powered by the radioactive decay of r-process nuclei, only few tentative identifications of light r-process elements have been made so far. One of the major limitations for the identification of heavy nuclei is incomplete or missing atomic data. While substantial progress has been made on lanthanide… ▽ More Even though the electromagnetic counterpart AT2017gfo to the binary neutron star merger GW170817 is powered by the radioactive decay of r-process nuclei, only few tentative identifications of light r-process elements have been made so far. One of the major limitations for the identification of heavy nuclei is incomplete or missing atomic data. While substantial progress has been made on lanthanide atomic data over the last few years, for actinides there has been less emphasis, with the first complete set of opacity data only recently published. We perform atomic structure calculations of neodymium $(Z=60)$ as well as the corresponding actinide uranium $(Z=92)$. Using two different codes (FAC and HFR) for the calculation of the atomic data, we investigate the accuracy of the calculated data (energy levels and electric dipole transitions) and their effect on kilonova opacities. For the FAC calculations, we optimise the local central potential and the number of included configurations and use a dedicated calibration technique to improve the agreement between theoretical and available experimental atomic energy levels (AELs). For ions with vast amounts of experimental data available, the presented opacities agree quite well with previous estimations. On the other hand, the optimisation and calibration method cannot be used for ions with only few available AELs. For these cases, where no experimental nor benchmarked calculations are available, a large spread in the opacities estimated from the atomic data obtained with the various atomic structure codes is observed.We find that the opacity of uranium is almost double the neodymium opacity. △ Less

Submitted 21 July, 2023; v1 submitted 3 February, 2023; originally announced February 2023.

Comments: 20 pages, 13 figures. Accepted by MNRAS

New experimental constraint on the $^{185}$W($n,γ$)$^{186}$W cross section

Authors: A. C. Larsen, G. M. Tveten, T. Renstrøm, H. Utsunomiya, E. Algin, T. Ari-izumi, K. O. Ay, F. L. Bello Garrote, L. Crespo Campo, F. Furmyr, S. Goriely, A. Görgen, M. Guttormsen, V. W. Ingeberg, B. V. Kheswa, I. K. B. Kullmann, T. Laplace, E. Lima, M. Markova, J. E. Midtbø, S. Miyamoto, A. H. Mjøs, V. Modamio, M. Ozgur, F. Pogliano , et al. (6 additional authors not shown)

Abstract: In this work, we present new data on the $^{182,183,184}$W($γ,n$) cross sections, utilizing a quasi-monochromatic photon beam produced at the NewSUBARU synchrotron radiation facility. Further, we have extracted the nuclear level density and $γ$-ray strength function of $^{186}$W from data on the $^{186}$W($α,α^\primeγ$)$^{186}$W reaction measured at the Oslo Cyclotron Laboratory. Combining previou… ▽ More In this work, we present new data on the $^{182,183,184}$W($γ,n$) cross sections, utilizing a quasi-monochromatic photon beam produced at the NewSUBARU synchrotron radiation facility. Further, we have extracted the nuclear level density and $γ$-ray strength function of $^{186}$W from data on the $^{186}$W($α,α^\primeγ$)$^{186}$W reaction measured at the Oslo Cyclotron Laboratory. Combining previous measurements on the $^{186}$W($γ,n$) cross section with our new $^{182,183,184}$W($γ,n$) and ($α,α^\primeγ$)$^{186}$W data sets, we have deduced the $^{186}$W $γ$-ray strength function in the range of $1 < E_γ< 6$ MeV and $7 < E_γ< 14$ MeV. Our data are used to extract the level density and $γ$-ray strength functions needed as input to the nuclear-reaction code \textsf{TALYS}, providing an indirect, experimental constraint for the $^{185}$W($n,γ$)$^{186}$W cross section and reaction rate. Compared to the recommended Maxwellian-averaged cross section (MACS) in the KADoNiS-1.0 data base, our results are on average lower for the relevant energy range $k_B T \in [5,100]$ keV, and we provide a smaller uncertainty for the MACS. The theoretical values of Bao \textit{et al.} and the cross section experimentally constrained on photoneutron data of Sonnabend \textit{et al.} are significantly higher than our result. The lower value by Mohr \textit{et al.} is in very good agreement with our deduced MACS. Our new results could have implications for the $s$-process and in particular the predicted $s$-process production of $^{186,187}$Os nuclei. △ Less

Submitted 30 January, 2023; originally announced January 2023.

Comments: 17 pages, 15 figures; to be submitted to Phys. Rev. C

Journal ref: Phys. Rev. C 108, 025804 (2023)

Comprehensive test of nuclear level density models

Authors: Stephane Goriely, Ann-Cecilie Larsen, Dennis Mücher

Abstract: For the last two decades, experimental information on nuclear level densities for about 60 different nuclei has been obtained on the basis of the Oslo method. While each of these measurements has been typically compared to one or a few level density models, a global study including all the measurements has been missing. The present study provides a systematic comparison between Oslo data and six g… ▽ More For the last two decades, experimental information on nuclear level densities for about 60 different nuclei has been obtained on the basis of the Oslo method. While each of these measurements has been typically compared to one or a few level density models, a global study including all the measurements has been missing. The present study provides a systematic comparison between Oslo data and six global level density models for 42 nuclei for which $s$-wave resonance spacings are also available. We apply a coherent normalization procedure to the Oslo data for each of the six different models, all being treated on the same footing. Our quantitative analysis shows that the constant-temperature model presents the best global description of the Oslo data, closely followed by the mean-field plus combinatorial model and Hartree-Fock plus statistical model. Their accuracies are quite similar, so that it remains difficult to clearly favour one of these models. When considering energies above the threshold where the experimental level scheme is complete, all the six models are shown to lead to rather similar accuracies with respect to Oslo data. The recently proposed shape method can, in principle, improve the situation since it provides an absolute estimate of the excitation-energy dependence of the measured level densities. We show for the specific case of $^{112}$Cd that the shape method could exclude the Hartree-Fock plus statistical model. Such an analysis remains to be performed for the bulk of data for which the shape method can be applied to the Oslo measurements before drawing conclusions on the general quality of a given nuclear level density model. △ Less

Submitted 26 January, 2023; originally announced January 2023.

Comments: 11 pages, 7 figures, published in PRC

Journal ref: Physical Review C 106, 044315 (2022)

Progress in Nuclear Astrophysics: a multi-disciplinary field with still many open questions

Authors: S. Goriely, A. Choplin, W. Ryssens, I. Kullmann

Abstract: 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 progres… ▽ More 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 s-and 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. △ Less

Submitted 5 December, 2022; originally announced December 2022.

Comments: 8 pages, 4 figures; Contribution to the proceedings of INPC 2022, Cape Town, South Africa

Can neutron star mergers alone explain the r-process enrichment of the Milky Way?

Authors: Chiaki Kobayashi, Ilya Mandel, Krzysztof Belczynski, Stephane Goriely, Thomas H. Janka, Oliver Just, Ashley J. Ruiter, Dany Van Beveren, Matthias U. Kruckow, Max M. Briel, Jan J. Eldridge, Elizabeth Stanway

Abstract: Comparing Galactic chemical evolution models to the observed elemental abundances in the Milky Way, we show that neutron star mergers can be a leading r-process site only if at low metallicities such mergers have very short delay times and significant ejecta masses that are facilitated by the masses of the compact objects. Namely, black hole-neutron star mergers, depending on the black-hole spins,… ▽ More Comparing Galactic chemical evolution models to the observed elemental abundances in the Milky Way, we show that neutron star mergers can be a leading r-process site only if at low metallicities such mergers have very short delay times and significant ejecta masses that are facilitated by the masses of the compact objects. Namely, black hole-neutron star mergers, depending on the black-hole spins, can play an important role in the early chemical enrichment of the Milky Way. We also show that none of the binary population synthesis models used in this paper, i.e., COMPAS, StarTrack, Brussels, ComBinE, and BPASS, can currently reproduce the elemental abundance observations. The predictions are problematic not only for neutron star mergers, but also for Type Ia supernovae, which may point to shortcomings in binary evolution models. △ Less

Submitted 18 December, 2022; v1 submitted 8 November, 2022; originally announced November 2022.

Comments: Accepted for publication in the Astrophysical Journal Letters, 15 pages, 4 figures, 1 table

PANDORA project: photo-nuclear reactions below $A=60$

Authors: A. Tamii, L. Pellegri, P. -A. Söderström, D. Allard, S. Goriely, T. Inakura, E. Khan, E. Kido, M. Kimura, E. Litvinova, S. Nagataki, P. von Neumann-Cosel, N. Pietralla, N. Shimizu, N. Tsoneva, Y. Utsuno, S. Adachi, P. Adsley, A. Bahini, D. Balabanski, B. Baret, J. A. C. Bekker, S. D. Binda, E. Boicu, A. Bracco , et al. (56 additional authors not shown)

Abstract: Photo-nuclear reactions of light nuclei below a mass of $A=60$ are studied experimentally and theoretically by the PANDORA (Photo-Absorption of Nuclei and Decay Observation for Reactions in Astrophysics) project. Two experimental methods, virtual-photon excitation by proton scattering and real-photo absorption by a high-brilliance gamma-ray beam produced by laser Compton scattering, will be applie… ▽ More Photo-nuclear reactions of light nuclei below a mass of $A=60$ are studied experimentally and theoretically by the PANDORA (Photo-Absorption of Nuclei and Decay Observation for Reactions in Astrophysics) project. Two experimental methods, virtual-photon excitation by proton scattering and real-photo absorption by a high-brilliance gamma-ray beam produced by laser Compton scattering, will be applied to measure the photo-absorption cross sections and the decay branching ratio of each decay channel as a function of the photon energy. Several nuclear models, e.g. anti-symmetrized molecular dynamics, mean-field type models, a large-scale shell model, and ab initio models, will be employed to predict the photo-nuclear reactions. The uncertainty in the model predictions will be evaluated from the discrepancies between the model predictions and the experimental data. The data and the predictions will be implemented in a general reaction calculation code TALYS . The results will be applied to the simulation of the photo-disintegration process of ultra-high-energy cosmic rays in inter-galactic propagation. △ Less

Submitted 18 November, 2022; v1 submitted 7 November, 2022; originally announced November 2022.

Synthesis of thorium and uranium in asymptotic giant branch stars

Authors: A. Choplin, S. Goriely, L. Siess

Abstract: The intermediate neutron capture process (i-process) operates at neutron densities between those of the slow and rapid neutron-capture processes. It can be triggered by the ingestion of protons in a convective helium-burning region. One possible astrophysical site is low-mass low-metallicity asymptotic giant branch (AGB) stars. We study here the possibility that actinides (particularly Th and U) m… ▽ More The intermediate neutron capture process (i-process) operates at neutron densities between those of the slow and rapid neutron-capture processes. It can be triggered by the ingestion of protons in a convective helium-burning region. One possible astrophysical site is low-mass low-metallicity asymptotic giant branch (AGB) stars. We study here the possibility that actinides (particularly Th and U) may be significantly synthesized through i-process nucleosynthesis in AGB stars. We computed a 1 $M_{\odot}$ model at [Fe/H] $= -2.5$ with the stellar evolution code STAREVOL. We used a nuclear network of 1160 species from H to Cf coupled to the transport processes. During the proton ingestion event, the neutron density goes up to $\sim 10^{15}$ cm$^{-3}$. While most of the nuclear flow cycles in the neutron-rich Pb-Bi-Po region, a non-negligible fraction leaks towards heavier elements and eventually synthesizes actinides. The surface enrichment in Th and U is subject to nuclear and astrophysical model uncertainties that could be lowered in the future, in particular by a detailed analysis of the nuclear inputs that affect the neutron capture rates of neutron-rich isotopes between Pb and Pa. One stellar candidate that may confirm the production of actinides by the i-process is the carbon-enhanced metal-poor r/s star J0949-1617, which shows Th lines in its spectrum. Its surface abundance is shown to be reasonably well reproduced by our AGB model. Combined with cosmochronometry, this finding opens the way to dating the i-process event and thus obtaining a lower limit for the age of CEMP-r/s stars. Such a dating is expected to be accurate only if surface abundances of Th and U can be extracted simultaneously. This work shows that actinides can be synthesized in AGB stars through the i-process. As a consequence, the r-process may not be the sole mechanism for the production of U and Th. △ Less

Submitted 7 November, 2022; originally announced November 2022.

Comments: 6 pages, 3 figures, accepted in A&A

Journal ref: A&A 667, L13 (2022)

The mass of odd-odd nuclei in microscopic mass models

Authors: W. Ryssens, G. Scamps, G. Grams, I. Kullmann, M. Bender, S. Goriely

Abstract: 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. Ho… ▽ More 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. △ Less

Submitted 7 November, 2022; originally announced November 2022.

Comments: 6 pages, 2 figures; Contribution to the proceedings of INPC 2022, Cape Town, South Africa

Progress on Brussels-Skyrme atomic mass models on a grid: stiff neutron matter equation of state

Authors: Guilherme Grams, Wouter Ryssens, Guillaume Scamps, Stephane Goriely, Nicolas Chamel

Abstract: 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… ▽ More 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. △ Less

Submitted 7 November, 2022; originally announced November 2022.

Comments: 3 pages, 2 figures

$^{178}$Hg and asymmetric fission of neutron-deficient pre-actinides

Authors: A. Jhingan, C. Schmitt, A. Lemasson, S. Biswas, Y. H. Kim, D. Ramos, A. N. Andreyev, D. Curien, M. Ciemala, E. Clément, O. Dorvaux, B. De Canditiis, F. Didierjean, G. Duchêne, J. Dudouet, J. Frankland, G. Frémont, J. Goupil, B. Jacquot, C. Raison, D. Ralet, B. -M. Retailleau, L. Stuttgé, I. Tsekhanovich, A. V. Andreev , et al. (5 additional authors not shown)

Abstract: Fission at low excitation energy is an ideal playground to probe the impact of nuclear structure on nuclear dynamics. While the importance of structural effects in the nascent fragments is well-established in the (trans-)actinide region, the observation of asymmetric fission in several neutron-deficient pre-actinides can be explained by various mechanisms. To deepen our insight into that puzzle, a… ▽ More Fission at low excitation energy is an ideal playground to probe the impact of nuclear structure on nuclear dynamics. While the importance of structural effects in the nascent fragments is well-established in the (trans-)actinide region, the observation of asymmetric fission in several neutron-deficient pre-actinides can be explained by various mechanisms. To deepen our insight into that puzzle, an innovative approach based on inverse kinematics and an enhanced version of the VAMOS++ heavy-ion spectrometer was implemented at the GANIL facility, Caen. Fission of $^{178}$Hg was induced by fusion of $^{124}$Xe and $^{54}$Fe. The two fragments were detected in coincidence using VAMOS++ supplemented with a new SEcond Detection arm. For the first time in the pre-actinide region, access to the pre-neutron mass and total kinetic energy distributions, and the simultaneous isotopic identification of one the fission fragment, was achieved. The present work describes the experimental approach, and discusses the pre-neutron observables in the context of an extended asymmetric-fission island located south-west of $^{208}Pb. A comparison with different models is performed, demonstrating the importance of this "new" asymmetric-fission island for elaborating on driving effects in fission. △ Less

Submitted 3 November, 2022; originally announced November 2022.

Journal ref: Phys. Rev. C 106, 044607 (2022)

Accreting neutron stars from the nuclear energy-density functional theory. II. Equation of state and global properties

Authors: A. F. Fantina, J. L. Zdunik, N. Chamel, J. M. Pearson, L. Suleiman, S. Goriely

Abstract: The accretion of matter onto the surface of a neutron star in a low-mass X-ray binary triggers X-ray bursts, whose ashes are buried and further processed thus altering the composition and the properties of the stellar crust. In this second paper of a series, the impact of accretion on the equation of state and on the global properties of neutron stars is studied in the framework of the nuclear ene… ▽ More The accretion of matter onto the surface of a neutron star in a low-mass X-ray binary triggers X-ray bursts, whose ashes are buried and further processed thus altering the composition and the properties of the stellar crust. In this second paper of a series, the impact of accretion on the equation of state and on the global properties of neutron stars is studied in the framework of the nuclear energy-density functional theory. Considering ashes made of $^{56}$Fe, we calculated the equations of state using the same Brussels-Montreal nuclear energy-density functionals BSk19, BSk20, and BSk21, as those already employed for determining the crustal heating in our previous study for the same ashes. All regions of accreting neutron stars were treated in a unified and thermodynamically consistent way. With these equations of state, we determined the mass, radius, moment of inertia, and tidal deformability of accreted neutron stars and compared with catalyzed neutron stars for which unified equations of state based on the same functionals are available. The equation of state of accreted neutron stars is found to be significantly stiffer than that of catalyzed matter, with an adiabatic index $Γ\approx 4/3$ throughout the crust. For this reason, accreting neutron stars have larger radii. However, their crustal moment of inertia and their tidal deformability are hardly changed provided density discontinuities at the interface between adjacent crustal layers are properly taken into account. The enhancement of the stiffness of the equation of state of accreting neutron stars is mainly a consequence of nuclear shell effects, thus confirming the importance of a quantum treatment as stressed in our first study. With our previous calculations of crustal heating using the same functionals, we have thus obtained consistent microscopic inputs for simulations of accreting neutron stars. △ Less

Submitted 23 September, 2022; originally announced September 2022.

Comments: 10 pages, 7 figures. The tables of the equations of state are available at the CDS

Journal ref: A&A 665, A74 (2022)

The intermediate neutron capture process. III. The i-process in AGB stars of different masses and metallicities without overshoot

Authors: A. Choplin, L. Siess, S. Goriely

Abstract: Alongside the slow (s) and rapid (r) neutron capture processes, an intermediate neutron capture process (i-process) is thought to exist. It happens when protons are mixed in a convective helium-burning zone, and is referred to as proton ingestion event (PIE). A possible astrophysical site is the asymptotic giant branch (AGB) phase of low-mass low-metallicity stars. We provide i-process yields of a… ▽ More Alongside the slow (s) and rapid (r) neutron capture processes, an intermediate neutron capture process (i-process) is thought to exist. It happens when protons are mixed in a convective helium-burning zone, and is referred to as proton ingestion event (PIE). A possible astrophysical site is the asymptotic giant branch (AGB) phase of low-mass low-metallicity stars. We provide i-process yields of a grid of AGB stars experiencing PIEs. We computed 12 models with initial masses of 1, 2, and 3 $M_{\odot}$ and metallicities of [Fe/H] $=-3.0$, $-2.5$ $-2.3,$ and $-2.0, $ with the stellar evolution code STAREVOL. We used a nuclear network of 1160 species at maximum, coupled to the chemical transport equations. These simulations do not include any extra mixing process. Proton ingestion takes place in six out of our 12 AGB models. These models experience i-process nucleosynthesis characterized by neutron densities of $\simeq 10^{14} -10^{15}$ cm$^{-3}$. Depending on the PIE properties two different evolution paths follow: either the stellar envelope is quickly lost and no more thermal pulses develop or the AGB phase resumes with additional thermal pulses. This behaviour critically depends on the pulse number when the PIE occurs, the mass of the ingested protons, and the extent to which the pulse material is diluted in the convective envelope. The surface enrichment after a PIE is a robust feature of our models and it persists under various convective assumptions. Our models can synthesise heavy elements up to Pb without any parametrized extra mixing process such as overshoot or inclusion of a $^{13}$C-pocket. Nevertheless, it remains to be explored how the i-process depends on mixing processes, such as overshoot, thermohaline, or rotation. △ Less

Submitted 21 September, 2022; originally announced September 2022.

Comments: 18 pages, 16 figures, accepted in A&A

Journal ref: A&A 667, A155 (2022)

Skyrme-Hartree-Fock-Bogoliubov mass models on a 3D mesh: II. Time-reversal symmetry breaking

Authors: Wouter Ryssens, Guillaume Scamps, S. Goriely, Michael Bender

Abstract: Models based on nuclear energy density functionals can provide access to a multitude of observables for thousands of nuclei in a single framework with microscopic foundations. Such models can rival the accuracy of more phenomenological approaches, but doing so requires adjusting parameters to thousands of nuclear masses. To keep such large-scale fits feasible, several symmetry restrictions are gen… ▽ More Models based on nuclear energy density functionals can provide access to a multitude of observables for thousands of nuclei in a single framework with microscopic foundations. Such models can rival the accuracy of more phenomenological approaches, but doing so requires adjusting parameters to thousands of nuclear masses. To keep such large-scale fits feasible, several symmetry restrictions are generally imposed on the nuclear configurations. One such example is time-reversal invariance, which is generally enforced via the Equal Filling Approximation (EFA). Here we lift this assumption, enabling us to access the spin and current densities in the ground states of odd-mass and odd-odd nuclei and which contribute to the total energy of such nuclei through so-called "time-odd" terms. We present here the Skyrme-based BSkG2 model whose parameters were adjusted to essentially all known nuclear masses without relying on the EFA, refining our earlier work [G. Scamps et al., EPJA 57, 333 (2021), arXiv:2011.07904]. Moving beyond ground state properties, we also incorporated information on the fission barriers of actinide nuclei in the parameter adjustment. The resulting model achieves a root-mean-square (rms) deviation of (i) 0.678 MeV on 2457 known masses, (ii) 0.027 fm on 884 measured charge radii, (iii) 0.44 MeV and 0.47 MeV, respectively, on 45 reference values for primary and secondary fission barriers of actinide nuclei, and (iv) 0.49 MeV on 28 fission isomer excitation energies. We limit ourselves here to a description of the model and the study the impact of lifting the EFA on ground state properties such as binding energies, deformation and pairing, deferring a detailed discussion of fission to a forthcoming paper. △ Less

Submitted 12 August, 2022; originally announced August 2022.

Comments: 30 pages, 15 figures

Report number: NT@UW-22-08

Cannibals in the thick disk II -- Radial-velocity monitoring of the young alpha-rich stars

Authors: P. Jofre, A. Jorissen, C. Aguilera-Gomez, S. Van Eck, J. Tayar, M. Pinsonneault, J. Zinn, S. Goriely, H. Van Winckel

Abstract: Determining ages of stars for reconstructing the history of the Milky Way remains one of the most difficult tasks in astrophysics. This involves knowing when it is possible to relate the stellar mass with its age and when it is not. The young $α-$rich (YAR) stars present such a case in which we are still not sure about their ages because they are relatively massive, implying young ages, but their… ▽ More Determining ages of stars for reconstructing the history of the Milky Way remains one of the most difficult tasks in astrophysics. This involves knowing when it is possible to relate the stellar mass with its age and when it is not. The young $α-$rich (YAR) stars present such a case in which we are still not sure about their ages because they are relatively massive, implying young ages, but their abundances are $α-$enhanced, which implies old ages. We report the results from new observations from a long-term radial-velocity-monitoring campaign complemented with high-resolution spectroscopy, as well as new astrometry and seismology of a sample of 41 red giants from the third version of APOKASC, which includes YAR stars. The aim is to better characterize the YAR stars in terms of binarity, mass, abundance trends, and kinematic properties.The radial velocities of HERMES, APOGEE, and Gaia were combined to determine the binary fraction among YAR stars. In combination with their mass estimate, evolutionary status, chemical composition, and kinematic properties, it allowed us to better constrain the nature of these objects. We found that stars with $\mathrm{M} < 1 \mathrm{M}_\odot$ were all single, whereas stars with $\mathrm{M} > 1 \mathrm{M}_\odot$ could be either single or binary. This is in agreement with theoretical predictions of population synthesis models. Studying their [C/N], [C/Fe], and [N/Fe], trends with mass, it became clear that many YAR stars do not follow the APOKASC stars, favoring the scenario that most of them are the product of mass transfer. Abr. △ Less

Submitted 28 November, 2022; v1 submitted 22 July, 2022; originally announced July 2022.

Comments: Accepted in A&A

Journal ref: A&A 671, A21 (2023)

Impact of systematic nuclear uncertainties on composition and decay heat of dynamical and disk ejecta in compact binary mergers

Authors: I. Kullmann, S. Goriely, O. Just, A. Bauswein, H. -T. Janka

Abstract: Theoretically predicted yields of elements created by the rapid neutron capture (r-) process carry potentially large uncertainties associated with incomplete knowledge of nuclear properties and approximative hydrodynamical modelling of the matter ejection processes. We present an in-depth study of the nuclear uncertainties by varying theoretical nuclear input models that describe the experimentall… ▽ More Theoretically predicted yields of elements created by the rapid neutron capture (r-) process carry potentially large uncertainties associated with incomplete knowledge of nuclear properties and approximative hydrodynamical modelling of the matter ejection processes. We present an in-depth study of the nuclear uncertainties by varying theoretical nuclear input models that describe the experimentally unknown neutron-rich nuclei. This includes two frameworks for calculating the radiative neutron capture rates and 14 different models for nuclear masses, $β$-decay rates and fission properties. Our r-process nuclear network calculations are based on detailed hydrodynamical simulations of dynamically ejected material from NS-NS or NS-BH binary mergers plus the secular ejecta from BH-torus systems. The impact of nuclear uncertainties on the r-process abundance distribution and the early radioactive heating rate is found to be modest (within a factor of $\sim20$ for individual $A>90$ abundances and a factor of 2 for the heating rate). However, the impact on the late-time heating rate is more significant and depends strongly on the contribution from fission. We witness significantly larger sensitivity to the nuclear physics input if only a single trajectory is used compared to considering ensembles of $\sim$200-300 trajectories, and the quantitative effects of the nuclear uncertainties strongly depend on the adopted conditions for the individual trajectory. We use the predicted Th/U ratio to estimate the cosmochronometric age of six metal-poor stars and find the impact of the nuclear uncertainties to be up to 2 Gyr. △ Less

Submitted 11 May, 2023; v1 submitted 15 July, 2022; originally announced July 2022.

Comments: 28 pages, 24 figures, accepted to MNRAS

Neutron star crust properties: comparison between the compressible liquid-drop model and the extended Thomas-Fermi approach

Authors: Guilherme Grams, Jerome Margueron, Rahul Somasundaram, Nicolas Chamel, Stephane Goriely

Abstract: We present a detailed analysis of three models predicting the properties of non-uniform matter in the crust of neutron stars: the compressible liquid-drop model, the fourth order Extended Thomas Fermi (ETF) method, and ETF plus Strutinsky integral (ETFSI) correction. The former treats the nuclear clusters as uniform hard spheres, the second takes into account the density distribution which can be… ▽ More We present a detailed analysis of three models predicting the properties of non-uniform matter in the crust of neutron stars: the compressible liquid-drop model, the fourth order Extended Thomas Fermi (ETF) method, and ETF plus Strutinsky integral (ETFSI) correction. The former treats the nuclear clusters as uniform hard spheres, the second takes into account the density distribution which can be different for neutrons and protons, and the last one includes the proton shell effects within the Strutinsky approach. The purpose of this work is to understand the importance of the improvements in the nuclear modeling and to analyze the quantities which are the most sensitive to them. We find that thermodynamic quantities such as pressure, energy and chemical potential, as well as the electron fraction, are in very good agreement among the three models. This confirms previous results where we have shown that the improvement in the finite-size description of the nuclear clusters has a small impact on these quantities, since they are mainly constrained by the bulk properties. The refinements in the finite-size modeling are shown to impact mostly the composition of the nuclear clusters ($Z_{cl}$, $N_{cl}$) in an ordering which ranks according to the leptodermous expansion. This analysis is performed considering both the r-cluster and the e-cluster representations. The proton shell effects are shown to stabilize $Z_{cl}$, which consequently impacts the neutron number $N_{cl}$ as well. △ Less

Submitted 30 May, 2022; originally announced May 2022.

Comments: 8 pages, 4 figures

The p-process in exploding rotating massive stars

Authors: A. Choplin, S. Goriely, R. Hirschi, N. Tominaga, G. Meynet

Abstract: 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 du… ▽ More 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. △ Less

Submitted 30 March, 2022; originally announced March 2022.

Comments: 12 pages, 13 figures, 2 tables, accepted for publication in A&A

Journal ref: A&A 661, A86 (2022)

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

Authors: H. Utsunomiya, Z. R. Hao, S. Goriely, X. G. Cao, G. T. Fan, H. W. Wang

Abstract: 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… ▽ More 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. △ Less

Submitted 24 March, 2022; originally announced March 2022.

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

Authors: J. H. Gillanders, S. J. Smartt, S. A. Sim, A. Bauswein, S. Goriely

Abstract: 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… ▽ More 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. △ Less

Submitted 3 May, 2022; v1 submitted 3 February, 2022; originally announced February 2022.

Comments: 20 pages, 14 figures and 3 tables. Accepted for publication in MNRAS

Tc-rich M stars: platypuses of low-mass star evolution

Authors: Shreeya Shetye, Sophie Van Eck, Alain Jorissen, Lionel Siess, Stephane Goriely

Abstract: The technetium-rich (Tc-rich) M stars reported in the literature (Little-Marenin & Little (1979); Uttenthaler et al. (2013)) are puzzling objects since no isotope of technetium has a half-life longer than a few million years, and 99Tc, the longest-lived isotope along the s-process path, is expected to be detected only in thermally-pulsing stars enriched with other s-process elements (like zirconiu… ▽ More The technetium-rich (Tc-rich) M stars reported in the literature (Little-Marenin & Little (1979); Uttenthaler et al. (2013)) are puzzling objects since no isotope of technetium has a half-life longer than a few million years, and 99Tc, the longest-lived isotope along the s-process path, is expected to be detected only in thermally-pulsing stars enriched with other s-process elements (like zirconium). Carbon should also be enriched, since it is dredged up at the same time, after each thermal pulse on the asymptotic giant branch (AGB). However, these Tc-enriched objects are classified as M stars, meaning that they neither have any significant zirconium enhancement (otherwise they would be tagged as S-type stars) nor any large carbon overabundance (in which case they would be carbon stars). Here we present the first detailed chemical analysis of a Tc-rich M-type star, namely S Her. We first confirm the detection of the Tc lines, and then analyze its carbon and s-process abundances, and draw conclusions on its evolutionary status. Understanding these Tc-rich M stars is an important step to constrain the threshold luminosity for the first occurrence of the third dredge-up and the composition of s-process ejecta during the very first thermal pulses on the AGB. △ Less

Submitted 8 January, 2022; originally announced January 2022.

Comments: 6 pages, 3 figures, Under review for 'The origin of outflows in evolved stars Proceedings IAU Symposium No. 366, 2021'

Multidimensional King-plot analysis for accurate extraction of Cd nuclear charge radii: a challenge for nuclear structure theory

Authors: J. Z. Han, C. Pan, K. Y. Zhang, X. F. Yang, S. Q. Zhang, J. C. Berengut, S. Goriely, H. Wang, Y. M. Yu, J. Meng, J. W. Zhang, L. J. Wang

Abstract: High-accuracy determination of isotope shift factors, which plays a crucial role in accurate extraction of nuclear charge radius, is well-known to be challenging experimentally and theoretically. Nonetheless, based on an accurate measurement of the isotope shifts for the Cd$^+$ $5s~^2S_{1/2}-5p~^2P_{3/2}$ and the Cd $5s^2~^1S_0-5s5p~^1P_1$ transition, a multidimensional King-plot analysis is perfo… ▽ More High-accuracy determination of isotope shift factors, which plays a crucial role in accurate extraction of nuclear charge radius, is well-known to be challenging experimentally and theoretically. Nonetheless, based on an accurate measurement of the isotope shifts for the Cd$^+$ $5s~^2S_{1/2}-5p~^2P_{3/2}$ and the Cd $5s^2~^1S_0-5s5p~^1P_1$ transition, a multidimensional King-plot analysis is performed to extract the atomic field shift and mass shift factors accurately. The results are further cross-checked against results from atomic structure calculations using a high-accuracy configuration interaction and many-body perturbation theory. Combined with previous isotope shifts, nuclear charge radii for $^{100-130}$Cd of the highest accuracy are reported. For the neutron-rich region, accuracies for the charge radii are improved by nearly one order of magnitude. This work provides a coherent and systematic extraction of Cd nuclear charge radii from isotope shifts. The results offer stringent testing and possible challenges to the latest advances in nuclear theory and excellent benchmarking to the atomic structure calculations. △ Less

Submitted 8 November, 2021; originally announced November 2021.

Dynamical ejecta of neutron star mergers with nucleonic weak processes II: Kilonova emission

Authors: Oliver Just, Ina Kullmann, Stephane Goriely, Andreas Bauswein, Hans-Thomas Janka, Christine E. Collins

Abstract: The majority of existing results for the kilonova (or macronova) emission from material ejected during a neutron-star (NS) merger is based on (quasi-)one-zone models or manually constructed toy-model ejecta configurations. In this study we present a kilonova analysis of the material ejected during the first ~10ms of a NS merger, called dynamical ejecta, using directly the outflow trajectories from… ▽ More The majority of existing results for the kilonova (or macronova) emission from material ejected during a neutron-star (NS) merger is based on (quasi-)one-zone models or manually constructed toy-model ejecta configurations. In this study we present a kilonova analysis of the material ejected during the first ~10ms of a NS merger, called dynamical ejecta, using directly the outflow trajectories from general relativistic smoothed-particle hydrodynamics simulations including a sophisticated neutrino treatment and the corresponding nucleosynthesis results, which have been presented in Part I of this study. We employ a multi-dimensional two-moment radiation transport scheme with approximate M1 closure to evolve the photon field and use a heuristic prescription for the opacities found by calibration with atomic-physics based reference results. We find that the photosphere is generically ellipsoidal but augmented with small-scale structure and produces emission that is about 1.5-3 times stronger towards the pole than the equator. The kilonova typically peaks after 0.7-1.5days in the near-infrared frequency regime with luminosities between 3-7x10^40erg/s and at photospheric temperatures of 2.2-2.8x10^3K. A softer equation of state or higher binary-mass asymmetry leads to a longer and brighter signal. Significant variations of the light curve are also obtained for models with artificially modified electron fractions, emphasizing the importance of a reliable neutrino-transport modeling. None of the models investigated here, which only consider dynamical ejecta, produces a transient as bright as AT2017gfo. The near-infrared peak of our models is incompatible with the early blue component of AT2017gfo. △ Less

Submitted 5 May, 2022; v1 submitted 29 September, 2021; originally announced September 2021.

Comments: 23 pages, 15 figures, 1 table, accepted to MNRAS; updated for correct version of Fig. A1, left panel, and corrected Eqs. 7 and 9

Systematical studies of the E1 photon strength functions combining Skyrme-HFB+QRPA model and experimental giant dipole resonance properties

Authors: Y. Xu, S. Goriely, E. Khan

Abstract: Valuable theoretical predictions of nuclear dipole excitations in the whole nuclear chart are of great interest for different applications, including in particular nuclear astrophysics. We present here the systematic study of the electric dipole (E1) photon strength functions (PSFs) combining the microscopic Hartree-Fock-Bogoliubov plus Quasiparticle Random Phase Approximation (HFB+QRPA) model and… ▽ More Valuable theoretical predictions of nuclear dipole excitations in the whole nuclear chart are of great interest for different applications, including in particular nuclear astrophysics. We present here the systematic study of the electric dipole (E1) photon strength functions (PSFs) combining the microscopic Hartree-Fock-Bogoliubov plus Quasiparticle Random Phase Approximation (HFB+QRPA) model and the parametrizations constrained by the available experimental giant dipole resonance (GDR) data. For about 10000 nuclei with 8<Z<124 lying between the proton and the neutron drip-lines on nuclear chart, the particle-hole strength distributions are computed using the HFB+QRPA model under the assumption of spherical symmetry and making use of the BSk27 Skyrme effective interaction derived from the most accurate HFB mass model (HFB-27) so far achieved. Large-scale calculations of the BSk27+QRPA E1 PSFs are performed in the framework of a specific folding procedure, in which three phenomenological improvements are considered. First, two interference factors are introduced and adjusted to reproduce at best the available experimental GDR data. Second, an empirical expression accounting for the deformation effect is applied to describe the peak splitting of the strength function. Third, the width of the strength function is corrected by a temperature-dependent term, which effectively increases the de-excitation photon strength function at low-energy. The E1 PSFs as well as the extracted GDR peaks and widths are compared with available experimental data. A relatively good agreement with data indicates the reliability of the calculations. Eventually, the astrophysical (n,g) rates for all the 10000 nuclei with 8<Z<124 are estimated using the present E1 PSFs. The resulting reaction rates are compared with previous BSk7+QRPA results and Gogny-HFB+QRPA predictions based on the D1M interaction. △ Less

Submitted 28 September, 2021; originally announced September 2021.

Comments: 16 pages, 14 figures, to be published on Phys. Rev. C

Dynamical ejecta of neutron star mergers with nucleonic weak processes I: Nucleosynthesis

Authors: I. Kullmann, S. Goriely, O. Just, R. Ardevol-Pulpillo, A. Bauswein, H. -T. Janka

Abstract: We present a coherent study of the impact of neutrino interactions on the r-process element nucleosynthesis and the heating rate produced by the radioactive elements synthesised in the dynamical ejecta of neutron star-neutron star (NS-NS) mergers. We have studied the material ejected from four NS-NS merger systems based on hydrodynamical simulations which handle neutrino effects in an elaborate wa… ▽ More We present a coherent study of the impact of neutrino interactions on the r-process element nucleosynthesis and the heating rate produced by the radioactive elements synthesised in the dynamical ejecta of neutron star-neutron star (NS-NS) mergers. We have studied the material ejected from four NS-NS merger systems based on hydrodynamical simulations which handle neutrino effects in an elaborate way by including neutrino equilibration with matter in optically thick regions and re-absorption in optically thin regions. We find that the neutron richness of the dynamical ejecta is significantly affected by the neutrinos emitted by the post-merger remnant, in particular when compared to a case neglecting all neutrino interactions. Our nucleosynthesis results show that a solar-like distribution of r- process elements with mass numbers $A \gtrsim 90$ is produced, including a significant enrichment in Sr and a reduced production of actinides compared to simulations without inclusion of the nucleonic weak processes. The composition of the dynamically ejected matter as well as the corresponding rate of radioactive decay heating are found to be rather independent of the system mass asymmetry and the adopted equation of state. This approximate degeneracy in abundance pattern and heating rates can be favourable for extracting the ejecta properties from kilonova observations, at least if the dynamical component dominates the overall ejecta. Part II of this work will study the light curve produced by the dynamical ejecta of our four NS merger models. △ Less

Submitted 20 November, 2021; v1 submitted 6 September, 2021; originally announced September 2021.

Comments: 18 pages, 21 figures, 2 Tables, accepted by MNRAS 2021-11-19

The intermediate neutron capture process II.Nuclear uncertainties

Authors: S. Goriely, L. Siess, A. Choplin

Abstract: Carbon-enhanced metal-poor (CEMP) r/s-stars show surface-abundance distributions characteristic of the so-called intermediate neutron capture process (i-process) of nucleosynthesis. We previously showed that the ingestion of protons in the convective helium-burning region of a low-mass low-metallicity star can explain the surface abundance distribution observed in CEMP r/s stars relatively well. S… ▽ More Carbon-enhanced metal-poor (CEMP) r/s-stars show surface-abundance distributions characteristic of the so-called intermediate neutron capture process (i-process) of nucleosynthesis. We previously showed that the ingestion of protons in the convective helium-burning region of a low-mass low-metallicity star can explain the surface abundance distribution observed in CEMP r/s stars relatively well. Such an i-process requires detailed reaction network calculations involving hundreds of nuclei for which reaction rates have not yet been determined experimentally. We investigate the nuclear physics uncertainties affecting the i-process during the AGB phase of low-metallicity low-mass stars by propagating the theoretical uncertainties in the radiative neutron capture cross sections, as well as the 13C(a,n)16O reaction rate, and estimating their impact on the surface-abundance distribution. It is found that considering systematic uncertainties on the various nuclear ingredients affecting the radiative neutron capture rates, surface elemental abundances are typically predicted within +/-0.4 dex. The 56 < Z < 59 region of the spectroscopically relevant heavy-s elements of Ba-La-Ce-Pr as well as the r-dominated Eu element remain relatively unaffected by nuclear uncertainties. In contrast, the inclusion of the direct capture contribution impacts the rates in the neutron-rich A~45, 100, 160, and 200 regions, and the i-process production of the Z~45 and 65-70 elements. Uncertainties in the photon strength function also impact the overabundance factors by typically 0.2-0.4 dex. Nuclear level densities tend to affect abundance predictions mainly in the Z=74-79 regions. The uncertainties associated with the neutron-producing reaction 13C(a,n)16O and the unknown beta-decay rates are found to have a low impact on the overall surface enrichment △ Less

Submitted 1 September, 2021; originally announced September 2021.

Comments: 9 pages, 8 Figures, A&A to be published

Journal ref: A&A 654, A129 (2021)

Sr and Ba abundance determinations: comparing machine-learning with star-by-star analyses -- High-resolution re-analysis of suspected LAMOST barium stars

Authors: D. Karinkuzhi, S. Van Eck, A. Jorissen, A. Escorza, S. Shetye, T. Merle, L. Siess, S. Goriely, H. Van Winckel

Abstract: A new large sample of 895 s-process-rich candidates out of 454180 giant stars surveyed by LAMOST at low spectral resolution (R ~ 1800) has been reported by Norfolk et al. (2019; hereafter N19). We aim at confirming the s-process enrichment at the higher resolution (R ~ 86000) offered by the HERMES-Mercator spectrograph, for the 15 brightest targets of the previous study sample which consists in 13… ▽ More A new large sample of 895 s-process-rich candidates out of 454180 giant stars surveyed by LAMOST at low spectral resolution (R ~ 1800) has been reported by Norfolk et al. (2019; hereafter N19). We aim at confirming the s-process enrichment at the higher resolution (R ~ 86000) offered by the HERMES-Mercator spectrograph, for the 15 brightest targets of the previous study sample which consists in 13 Sr-only stars and two Ba-only stars. Abundances were derived for elements Li, C, N, O, Na, Mg, Fe, Rb, Sr, Y, Zr, Nb, Ba, La, and Ce. Binarity has been tested by comparing the Gaia DR2 radial velocity with the HERMES velocity obtained 1600 - 1800 days later. Among the 15 programme stars, four show no s-process overabundances ([X/Fe] < 0.2 dex), eight show mild s-process overabundances (at least three heavy elements with 0.2 < [X/Fe] < 0.8), and three have strong overabundances (at least three heavy elements with [X/Fe] > 0.8). Among the 13 stars classified as Sr-only by the previous investigation, four have no s-process overabundances, eight are mild barium stars, and one is a strong barium star. The two Ba-only stars turn out to be both strong barium stars and are actually dwarf barium stars. They also show clear evidence for being binaries. Among the no-s stars, there are two binaries out of four, whereas only one out of the eight mild barium stars show a clear signature of radial-velocity variations. Blending effects and saturated lines have to be considered very carefully when using machine-learning techniques, especially on low-resolution spectra. Among the Sr-only stars from the previous study sample, one may expect about 60% (8/13) of them to be true mild barium stars and about 8% to be strong barium stars, and this fraction is likely close to 100% for the previous study Ba-only stars (2/2). △ Less

Submitted 18 July, 2021; originally announced July 2021.

Comments: Accepted for publication in A&A, 15 pages, 12 figures

Journal ref: A&A 654, A140 (2021)

S stars and s-process in the Gaia era II. Constraining the luminosity of the third dredge-up with Tc-rich S stars

Authors: Shreeya Shetye, Sophie Van Eck, Alain Jorissen, Stephane Goriely, Lionel Siess, Hans Van Winckel, Bertrand Plez, Michel Godefroid, George Wallerstein

Abstract: S stars are late-type giants that are transition objects between M-type stars and carbon stars on the asymptotic giant branch (AGB). They are classified into two types: intrinsic or extrinsic, based on the presence or absence of technetium (Tc). The Tc-rich or intrinsic S stars are thermally-pulsing (TP-)AGB stars internally producing s-process elements (including Tc) which are brought to their su… ▽ More S stars are late-type giants that are transition objects between M-type stars and carbon stars on the asymptotic giant branch (AGB). They are classified into two types: intrinsic or extrinsic, based on the presence or absence of technetium (Tc). The Tc-rich or intrinsic S stars are thermally-pulsing (TP-)AGB stars internally producing s-process elements (including Tc) which are brought to their surface via the third dredge-up (TDU). Tc-poor or extrinsic S stars gained their s-process overabundances via accretion of s-process-rich material from an AGB companion which has since turned into a dim white dwarf. Our goal is to investigate the evolutionary status of Tc-rich S stars by locating them in a Hertzsprung-Russell (HR) diagram using the results of Gaia early Data Release 3 (EDR3). We combine the current sample of 13 Tc-rich stars with our previous studies of 10 Tc-rich stars to determine the observational onset of the TDU in the metallicity range [-0.7; 0]. We also compare our abundance determinations with dedicated AGB nucleosynthesis predictions. The stellar parameters are derived using an iterative tool which combines HERMES high-resolution spectra, accurate Gaia EDR3 parallaxes, stellar evolution models and tailored MARCS model atmospheres for S-type stars. Using these stellar parameters we determine the heavy-element abundances by line synthesis. In the HR diagram, the intrinsic S stars are located at higher luminosities than the predicted onset of the TDU. These findings are consistent with Tc-rich S stars being genuinely TP-AGB stars. The comparison of the derived s-process abundance profiles of our intrinsic S stars with the nucleosynthesis predictions provide an overall good agreement. Stars with highest [s/Fe] tend to have the highest C/O ratios. △ Less

Submitted 7 April, 2021; originally announced April 2021.

Comments: 22 pages, 10 figures, accepted for publication in Section 7. Stellar structure and evolution of Astronomy and Astrophysics (A&A)

Journal ref: A&A 650, A118 (2021)

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

Authors: A. Choplin, L. Siess, S. Goriely

Abstract: 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 investigat… ▽ More 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. △ Less

Submitted 7 June, 2022; v1 submitted 17 February, 2021; originally announced February 2021.

Comments: Published in A&A. This version of the preprint has been corrected (especially Table 2 and Section 6.5) with the modifications appearing in the A&A corrigendum (which is included at the end of the paper)

Journal ref: A&A 662, C3 (2022)

Neutrino absorption and other physics dependencies in neutrino-cooled black-hole accretion disks

Authors: Oliver Just, Stephane Goriely, Hans-Thomas Janka, Shigehiro Nagataki, Andreas Bauswein

Abstract: Black-hole (BH) accretion disks formed in compact-object mergers or collapsars may be major sites of the rapid-neutron-capture (r-)process, but the conditions determining the electron fraction (Y_e) remain uncertain given the complexity of neutrino transfer and angular-momentum transport. After discussing relevant weak-interaction regimes, we study the role of neutrino absorption for shaping Y_e u… ▽ More Black-hole (BH) accretion disks formed in compact-object mergers or collapsars may be major sites of the rapid-neutron-capture (r-)process, but the conditions determining the electron fraction (Y_e) remain uncertain given the complexity of neutrino transfer and angular-momentum transport. After discussing relevant weak-interaction regimes, we study the role of neutrino absorption for shaping Y_e using an extensive set of simulations performed with two-moment neutrino transport and again without neutrino absorption. We vary the torus mass, BH mass and spin, and examine the impact of rest-mass and weak-magnetism corrections in the neutrino rates. We also test the dependence on the angular-momentum transport treatment by comparing axisymmetric models using the standard alpha-viscosity with viscous models assuming constant viscous length scales (l_t) and three-dimensional magnetohydrodynamic (MHD) simulations. Finally, we discuss the nucleosynthesis yields and basic kilonova properties. We find that absorption pushes Y_e towards ~0.5 outside the torus, while inside increasing the equilibrium value Y_e^eq by ~0.05-0.2. Correspondingly, a substantial ejecta fraction is pushed above Y_e=0.25, leading to a reduced lanthanide fraction and a brighter, earlier, and bluer kilonova than without absorption. More compact tori with higher neutrino optical depth, tau, tend to have lower Y_e^eq up to tau~1-10, above which absorption becomes strong enough to reverse this trend. Disk ejecta are less (more) neutron-rich when employing an l_t=const. viscosity (MHD treatment). The solar-like abundance pattern found for our MHD model marginally supports collapsar disks as major r-process sites, although a strong r-process may be limited to phases of high mass-infall rates, Mdot>~ 2 x 10^(-2) Msun/s. △ Less

Submitted 25 October, 2021; v1 submitted 16 February, 2021; originally announced February 2021.

Comments: 37 pages, 22 figures, 5 tables, accepted to MNRAS, minor changes compared to previous version

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

Authors: J. -F. Lemaître, S. Goriely, A. Bauswein, H. -T. Janka

Abstract: 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-ric… ▽ More 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. △ Less

Submitted 6 February, 2021; originally announced February 2021.

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

Authors: K. L. Malatji, K. S. Beckmann, M. Wiedeking, S. Siem, S. Goriely, A. C. Larsen, K. O. Ay, F. L. Bello Garrote, L. Crespo Campo, A. Görgen, M. Guttormsen, V. W. Ingeberg, P. Jones, B. V. Kheswa, P. von Neumann-Cosel, M. Ozgur, G. Potel, L. Pellegri, T. Renstrøm, G. M. Tveten, F. Zeiser

Abstract: 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 en… ▽ More 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. △ Less

Submitted 7 January, 2021; originally announced January 2021.

Journal ref: PRC 103 (2021) 014309

Shell-model based study of the direct capture in neutron-rich nuclei

Authors: K. Sieja, S. Goriely

Abstract: The radiative neutron capture rates for isotopes of astrophysical interest are commonly calculated within the statistical Hauser-Feshbach reaction model. Such an approach, assuming a high level density in the compound system, can be questioned in light and neutron-rich nuclei for which only a few or no resonant states are available. Therefore, in this work we focus on the direct neutron-capture pr… ▽ More The radiative neutron capture rates for isotopes of astrophysical interest are commonly calculated within the statistical Hauser-Feshbach reaction model. Such an approach, assuming a high level density in the compound system, can be questioned in light and neutron-rich nuclei for which only a few or no resonant states are available. Therefore, in this work we focus on the direct neutron-capture process. We employ a shell-model approach in several model spaces with well-established effective interactions to calculate spectra and spectroscopic factors in a set of 50 neutron-rich target nuclei in different mass regions, including doubly-, semi-magic and deformed ones. Those theoretical energies and spectroscopic factors are used to evaluate direct neutron capture rates and to test global theoretical models using average spectroscopic factors and level densities based on the Hartree-Fock-Bogoliubov plus combinatorial method. The comparison of shell-model and global model results reveals several discrepancies that can be related to problems in level densities. All the results show however that the direct capture is non-negligible with respect to the by-default Hauser-Feshbach predictions and can be even 100 times more important for the most neutron-rich nuclei close to the neutron drip line. △ Less

Submitted 2 April, 2021; v1 submitted 25 November, 2020; originally announced November 2020.

Comments: 10 pages, 9 figures, submitted to EPJ A

Journal ref: Eur. Phys. J. A 57, 110 (2021)

Skyrme-Hartree-Fock-Bogoliubov mass models on a 3D Mesh: Effect of triaxial shape

Authors: Guillaume Scamps, Stephane Goriely, Erik Olsen, Michael Bender, Wouter Ryssens

Abstract: The modeling of nuclear reactions and radioactive decays in astrophysical or earth-based conditions requires detailed knowledge of the masses of essentially all nuclei. Microscopic mass models based on nuclear energy density functionals (EDFs) can be descriptive and used to provide this information. The concept of intrinsic symmetry breaking is central to the predictive power of EDF approaches, ye… ▽ More The modeling of nuclear reactions and radioactive decays in astrophysical or earth-based conditions requires detailed knowledge of the masses of essentially all nuclei. Microscopic mass models based on nuclear energy density functionals (EDFs) can be descriptive and used to provide this information. The concept of intrinsic symmetry breaking is central to the predictive power of EDF approaches, yet is generally not exploited to the utmost by mass models because of the computational demands of adjusting up to about two dozen parameters to thousands of nuclear masses. We report on a first step to bridge the gap between what is presently feasible for studies of individual nuclei and large-scale models: we present a new Skyrme-EDF-based model that was adjusted using a three-dimensional coordinate-space representation, for the first time allowing for both axial and triaxial deformations during the adjustment process. To compensate for the substantial increase in computational cost brought by the latter, we have employed a committee of multilayer neural networks to model the objective function in parameter space and guide us towards the overall best fit. The resulting mass model BSkG1 is computed with the EDF model independently of the neural network. It yields a root mean square (rms) deviation on the 2457 known masses of 741 keV and an rms deviation on the 884 measured charge radii of 0.024 fm. △ Less

Submitted 12 October, 2021; v1 submitted 16 November, 2020; originally announced November 2020.

Comments: 25 pages, 15 figures

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

Authors: D. Karinkuzhi, S. Van Eck, S. Goriely, L. Siess, A. Jorissen, T. Merle, A. Escorza, T. Masseron

Abstract: 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 abu… ▽ More 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. △ Less

Submitted 18 January, 2021; v1 submitted 26 October, 2020; originally announced October 2020.

Comments: 35 pages, 26 figures, Accepted for publication in A&A, Final version

Journal ref: A&A 645, A61 (2021)

The impact of turbulent mixing on the galactic r-process enrichment by binary neutron star mergers

Authors: Irina Dvorkin, Frédéric Daigne, Stéphane Goriely, Elisabeth Vangioni, Joseph Silk

Abstract: We study the enrichment of the interstellar medium with rapid neutron capture (r-process) elements produced in binary neutron star (BNS) mergers. We use a semi-analytic model to describe galactic evolution, with merger rates and time delay distributions of BNS mergers consistent with the latest population synthesis models. In order to study the dispersion of the relative abundances of r-process el… ▽ More We study the enrichment of the interstellar medium with rapid neutron capture (r-process) elements produced in binary neutron star (BNS) mergers. We use a semi-analytic model to describe galactic evolution, with merger rates and time delay distributions of BNS mergers consistent with the latest population synthesis models. In order to study the dispersion of the relative abundances of r-process elements and iron, we applied a turbulent mixing scheme, where the freshly synthesized elements are gradually dispersed in the interstellar medium. We show that within our model the abundances observed in Milky-Way stars, in particular the scatter at low metallicities, can be entirely explained by BNS mergers. Our results suggest that binary neutron star mergers could be the dominant source of r-process elements in the Galaxy. △ Less

Submitted 20 July, 2021; v1 submitted 1 October, 2020; originally announced October 2020.

Comments: 16 pages, 19 figures. Modified to match the version accepted by MNRAS (some changes to main results and conclusions, new section added)

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

Authors: S. Shetye, S. Van Eck, S. Goriely, L. Siess, A. Jorissen, A. Escorza, H. Van Winckel

Abstract: 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 addit… ▽ More 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. △ Less

Submitted 19 February, 2020; originally announced February 2020.

Comments: A & A Letters, accepted for publication

Journal ref: A&A 635, L6 (2020)