-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
Astronuclear Physics: a Tale of the Atomic Nuclei in the Skies
Authors:
M. Arnould,
S. Goriely
Abstract:
A century ago, nuclear physics entered astrophysics, giving birth to a new field of science referred to as "Nuclear Astrophysics". With time, it developed at an impressive pace into a vastly inter- and multidisciplinary discipline bringing into its wake not only astronomy and cosmology, but also many other sub-fields of physics, especially particle, solid-state and computational physics, as well a…
▽ More
A century ago, nuclear physics entered astrophysics, giving birth to a new field of science referred to as "Nuclear Astrophysics". With time, it developed at an impressive pace into a vastly inter- and multidisciplinary discipline bringing into its wake not only astronomy and cosmology, but also many other sub-fields of physics, especially particle, solid-state and computational physics, as well as chemistry, geology and even biology. The present Astronuclear Physics review focusses primarily on the facets of nuclear physics that are of relevance to astronomy and astrophysics, the theoretical aspects being of special concern here.
△ Less
Submitted 30 January, 2020;
originally announced January 2020.
-
Role of the Symmetry Energy on the Structure of Neutron Stars with Unified Equations of State
Authors:
Nicolas Chamel,
John Michael Pearson,
Alexander Y. Potekhin,
Anthea F. Fantina,
Camille Ducoin,
Anup Dutta,
Stéphane Goriely,
Loïc Perot
Abstract:
The role of the symmetry energy on the internal constitution and the global structure of a cold nonaccreted neutron star is studied using a set of unified equations of state. Based on the nuclear energy-density functional theory, these equations of state provide a thermodynamically consistent treatment of all regions of the star and were calculated using the four different Brussels-Montreal functi…
▽ More
The role of the symmetry energy on the internal constitution and the global structure of a cold nonaccreted neutron star is studied using a set of unified equations of state. Based on the nuclear energy-density functional theory, these equations of state provide a thermodynamically consistent treatment of all regions of the star and were calculated using the four different Brussels-Montreal functionals BSk22, BSk24, BSk25 and BSk26. Our predictions are compared to various constraints inferred from astrophysical observations including the recent detection of the gravitational wave signal GW170817 from a binary neutron-star merger.
△ Less
Submitted 29 April, 2019;
originally announced April 2019.
-
Observational evidence of third dredge-up occurrence in S-type stars with initial masses around 1 Msun
Authors:
Shreeya Shetye,
Stephane Goriely,
Lionel Siess,
Sophie Van Eck,
Alain Jorissen,
Hans Van Winckel
Abstract:
Context- S stars are late-type giants with spectra showing characteristic molecular bands of ZrO in addition to the TiO bands typical of M stars. Their overabundance pattern shows the signature of s-process nucleosynthesis. Intrinsic, technetium (Tc)-rich S stars are the first objects, on the Asymptotic Giant Branch (AGB), to undergo third dredge-up (TDU) events. Gaia exquisite parallaxes now allo…
▽ More
Context- S stars are late-type giants with spectra showing characteristic molecular bands of ZrO in addition to the TiO bands typical of M stars. Their overabundance pattern shows the signature of s-process nucleosynthesis. Intrinsic, technetium (Tc)-rich S stars are the first objects, on the Asymptotic Giant Branch (AGB), to undergo third dredge-up (TDU) events. Gaia exquisite parallaxes now allow to precisely locate these stars in the Hertzsprung-Russell (HR) diagram. Here we report on a population of low-mass, Tc-rich S stars, previously unaccounted for by stellar evolution models. Aims- Our aim is to derive parameters of a sample of low-mass Tc-rich S stars and then, by comparing their location in the HR diagram with stellar evolution tracks, to derive their masses and to compare their measured s-process abundance profiles with recently derived STAREVOL nucleosynthetic predictions for low-mass AGB stars. Methods- The stellar parameters were obtained using a combination of HERMES high-resolution spectra, accurate Gaia Data Release 2 (Gaia-DR2) parallaxes, stellar-evolution models and newly-designed MARCS model atmospheres for S-type stars. Results- We report on 6 Tc-rich S stars lying close to the 1 Msun (initial mass) tracks of AGB stars of the corresponding metallicity and above the predicted onset of TDU, as expected. This provides direct evidence for TDUs occurring in AGB stars with initial masses as low as ~ 1 Msun and at low luminosity, i.e. at the start of the thermally-pulsing AGB. We present AGB models producing TDU in those stars with [Fe/H] in the range -0.25 to -0.5. There is a reasonable agreement between the measured and predicted s-process abundance profiles. For 2 objects however (CD -29 5912 and BD +34 1698), the predicted C/O ratio and s-process enhancements do not match simultaneously the measured ones.
△ Less
Submitted 8 April, 2019;
originally announced April 2019.
-
Unified equations of state for cold non-accreting neutron stars with Brussels-Montreal functionals. I. Role of symmetry energy
Authors:
J. M. Pearson,
N. Chamel,
A. Y. Potekhin,
A. F. Fantina,
C. Ducoin,
A. K. Dutta,
S. Goriely
Abstract:
The theory of the nuclear energy-density functional is used to provide a unified and thermodynamically consistent treatment of all regions of cold non-accreting neutron stars. In order to assess the impact of our lack of complete knowledge of the density dependence of the symmetry energy on the constitution and the global structure of neutron stars, we employ four different functionals. All of the…
▽ More
The theory of the nuclear energy-density functional is used to provide a unified and thermodynamically consistent treatment of all regions of cold non-accreting neutron stars. In order to assess the impact of our lack of complete knowledge of the density dependence of the symmetry energy on the constitution and the global structure of neutron stars, we employ four different functionals. All of them were precision fitted to essentially all the nuclear-mass data with the Hartree-Fock-Bogoliubov method and two different neutron-matter equations of state based on realistic nuclear forces. For each functional, we calculate the composition, the pressure-density relation, and the chemical potentials throughout the star. We show that uncertainties in the symmetry energy can significantly affect the theoretical results for the composition and global structure of neutron stars. To facilitate astrophysical applications, we construct analytic fits to our numerical results.
△ Less
Submitted 12 March, 2019;
originally announced March 2019.
-
Catching Element Formation In The Act
Authors:
Chris L. Fryer,
Frank Timmes,
Aimee L. Hungerford,
Aaron Couture,
Fred Adams,
Wako Aoki,
Almudena Arcones,
David Arnett,
Katie Auchettl,
Melina Avila,
Carles Badenes,
Eddie Baron,
Andreas Bauswein,
John Beacom,
Jeff Blackmon,
Stephane Blondin,
Peter Bloser,
Steve Boggs,
Alan Boss,
Terri Brandt,
Eduardo Bravo,
Ed Brown,
Peter Brown,
Steve Bruenn. Carl Budtz-Jorgensen,
Eric Burns
, et al. (194 additional authors not shown)
Abstract:
Gamma-ray astronomy explores the most energetic photons in nature to address some of the most pressing puzzles in contemporary astrophysics. It encompasses a wide range of objects and phenomena: stars, supernovae, novae, neutron stars, stellar-mass black holes, nucleosynthesis, the interstellar medium, cosmic rays and relativistic-particle acceleration, and the evolution of galaxies. MeV gamma-ray…
▽ More
Gamma-ray astronomy explores the most energetic photons in nature to address some of the most pressing puzzles in contemporary astrophysics. It encompasses a wide range of objects and phenomena: stars, supernovae, novae, neutron stars, stellar-mass black holes, nucleosynthesis, the interstellar medium, cosmic rays and relativistic-particle acceleration, and the evolution of galaxies. MeV gamma-rays provide a unique probe of nuclear processes in astronomy, directly measuring radioactive decay, nuclear de-excitation, and positron annihilation. The substantial information carried by gamma-ray photons allows us to see deeper into these objects, the bulk of the power is often emitted at gamma-ray energies, and radioactivity provides a natural physical clock that adds unique information. New science will be driven by time-domain population studies at gamma-ray energies. This science is enabled by next-generation gamma-ray instruments with one to two orders of magnitude better sensitivity, larger sky coverage, and faster cadence than all previous gamma-ray instruments. This transformative capability permits: (a) the accurate identification of the gamma-ray emitting objects and correlations with observations taken at other wavelengths and with other messengers; (b) construction of new gamma-ray maps of the Milky Way and other nearby galaxies where extended regions are distinguished from point sources; and (c) considerable serendipitous science of scarce events -- nearby neutron star mergers, for example. Advances in technology push the performance of new gamma-ray instruments to address a wide set of astrophysical questions.
△ Less
Submitted 7 February, 2019;
originally announced February 2019.
-
Probing stellar evolution with S stars and Gaia
Authors:
Shreeya Shetye,
Sophie Van Eck,
Alain Jorissen,
Hans Van Winckel,
Lionel Siess,
Stephane Goriely
Abstract:
S-type stars are late-type giants enhanced with s-process elements originating either from nucleosynthesis during the Asymptotic Giant Branch (AGB) or from a pollution by a binary companion. The former are called intrinsic S stars, and the latter extrinsic S stars. The intrinsic S stars are on the AGB and have undergone third dredge-up events. The atmospheric parameters of S stars are more numerou…
▽ More
S-type stars are late-type giants enhanced with s-process elements originating either from nucleosynthesis during the Asymptotic Giant Branch (AGB) or from a pollution by a binary companion. The former are called intrinsic S stars, and the latter extrinsic S stars. The intrinsic S stars are on the AGB and have undergone third dredge-up events. The atmospheric parameters of S stars are more numerous than those of M-type giants (C/O ratio and s-process abundances affect the thermal structure and spectral synthesis), and hence they are more difficult to derive. These atmospheric parameters are also entangled within each other. Nevertheless, high-resolution spectroscopic data of S stars combined with the Gaia Data Release 2 (GDR2) parallaxes and with the MARCS model atmospheres for S-type stars were used to derive effective temperatures, surface gravities, and luminosities. These parameters not only allow to locate the intrinsic and extrinsic S stars in the Hertzsprung-Russell (HR) diagram but also allow the accurate abundance analysis of the s-process elements.
△ Less
Submitted 17 October, 2018;
originally announced October 2018.
-
S stars and s-process in the Gaia era I. Stellar parameters and chemical abundances in a sub-sample of S stars with new MARCS model atmospheres
Authors:
Shreeya Shetye,
Sophie Van Eck,
Alain Jorissen,
Hans Van Winckel,
Lionel Siess,
Stephane Goriely,
Ana Escorza,
Drisya Karinkuzhi,
Bertrand Plez
Abstract:
S stars are transition objects between M-type giants and carbon stars on the asymptotic giant branch (AGB). They are characterized by overabundances of s-process elements. Roughly half of them are enhanced in technetium (Tc), an s-process element with no stable isotope, while the other half lack technetium. This dichotomy arises from the fact that Tc-rich S stars are intrinsically producing s-proc…
▽ More
S stars are transition objects between M-type giants and carbon stars on the asymptotic giant branch (AGB). They are characterized by overabundances of s-process elements. Roughly half of them are enhanced in technetium (Tc), an s-process element with no stable isotope, while the other half lack technetium. This dichotomy arises from the fact that Tc-rich S stars are intrinsically producing s-process elements and have undergone third dredge-up (TDU) events, while Tc-poor S stars owe their s-process overabundances to a past pollution by a former AGB companion. Our aim is to analyse the abundances of S stars and gain insights into their evolutionary status and on the nucleosynthesis of heavy s-process elements taking place in their interior. In particular, the location of extrinsic and intrinsic S stars in the HR diagram will be compared with the theoretical onset of the TDU on the thermally-pulsing AGB. A sample of 19 S-type stars was analysed by combining HERMES high-resolution spectra, accurate Gaia Data Release 2 (GDR2) parallaxes, stellar-evolution models, and newly-designed MARCS model atmospheres for S-type stars. Combining the derived parameters with GDR2 parallaxes allows a joint analysis of the location of the stars in the Hertzsprung-Russell diagram and of their surface abundances. For all 19 stars, Zr and Nb abundances are derived, complemented by abundances of other s-process elements for the three Tc-rich S stars. These abundances agree within the uncertainties with nucleosynthesis predictions for stars of corresponding mass, metallicity and evolutionary stage. Most extrinsic S stars lie close to the tip of the red giant branch (RGB), and a few are located along the early AGB. The location of intrinsic S stars in the HR diagram is compatible with them being thermally-pulsing AGB stars.
△ Less
Submitted 16 October, 2018;
originally announced October 2018.
-
Determination of the photodisintegration reaction rates involving charged particles: systematical calculations and proposed measurements based on Extreme Light Infrastructure - Nuclear Physics (ELI-NP)
Authors:
H. Y. Lan,
Y. Xu,
W. Luo,
D. L. Balabanski,
S. Goriely,
M. La Cognata,
C. Matei,
A. Anzalone,
S. Chesnevskaya,
G. L. Guardo,
D. Lattuada,
R. G. Pizzone,
S. Romano,
C. Spitaleri,
A. Taffara,
A. Tumino,
Z. C. Zhu
Abstract:
Photodisintegration reaction rates involving charged particles are of relevance to the p-process nucleosynthesis that aims at explaining the production of the stable neutron-deficient nuclides heavier than iron. In this study, the cross sections and astrophysical rates of (g,p) and (g,a) reactions for about 3000 target nuclei with 10<Z<100 ranging from stable to proton dripline nuclei are computed…
▽ More
Photodisintegration reaction rates involving charged particles are of relevance to the p-process nucleosynthesis that aims at explaining the production of the stable neutron-deficient nuclides heavier than iron. In this study, the cross sections and astrophysical rates of (g,p) and (g,a) reactions for about 3000 target nuclei with 10<Z<100 ranging from stable to proton dripline nuclei are computed. To study the sensitivity of the calculations to the optical model potentials (OMPs), both the phenomenological Woods-Saxon and the microscopic folding OMPs are taken into account. The systematic comparisons show that the reaction rates, especially for the (g,a) reaction, are dramatically influenced by the OMPs. Thus the better determination of the OMP is crucial to reduce the uncertainties of the photodisintegration reaction rates involving charged particles. Meanwhile, a gamma-beam facility at ELI-NP is being developed, which will open new opportunities to experimentally study the photodisintegration reactions of astrophysics interest. Considering both the important reactions identified by the nucleosynthesis studies and the purpose of complementing the experimental results for the reactions involving p-nuclei, the measurements of six (g,p) and eight (g,a) reactions based on the gamma-beam facility at ELI-NP and the ELISSA detector for the charged particles detection are proposed, and the GEANT4 simulations are correspondingly performed. The minimum required energies of the gamma-beam to measure these reactions are estimated. It is shown that the direct measurements of these photonuclear reactions within the Gamow windows at T_9=2.5 for p-process are fairly feasible and promising at ELI-NP. The expected experimental results will be used to constrain the OMPs of the charged particles, which can eventually reduce the uncertainties of the reaction rates for the p-process nucleosynthesis.
△ Less
Submitted 11 October, 2018;
originally announced October 2018.
-
When binaries keep track of recent nucleosynthesis: The Zr - Nb pair in extrinsic stars as an s-process diagnostic
Authors:
D. Karinkuzhi,
S. Van Eck,
A. Jorissen,
S. Goriely,
L. Siess,
T. Merle,
A. Escorza,
M. Van der Swaelmen,
H. M. J. Boffin,
T. Masseron,
S. Shetye,
B. Plez
Abstract:
Barium stars are s-process enriched giants. They owe their chemical peculiarities to a past mass transfer phase. During which they were polluted by their binary companion, which at the time was an AGB star, but is now an extinct white dwarf. Barium stars are thus ideal targets for understanding and constraining the s-process in low and intermediate-mass AGB stars. We derive the abundances of a lar…
▽ More
Barium stars are s-process enriched giants. They owe their chemical peculiarities to a past mass transfer phase. During which they were polluted by their binary companion, which at the time was an AGB star, but is now an extinct white dwarf. Barium stars are thus ideal targets for understanding and constraining the s-process in low and intermediate-mass AGB stars. We derive the abundances of a large number of heavy elements in order to shed light on the conditions of operation of the neutron source responsible for the production of s-elements in the former companions of the barium stars. Adopting a recently used methodology(Neyskens et al. 2015), we analyse a sample of 18 highly enriched barium stars observed with the high-resolution HERMES spectrograph mounted on the MERCATOR telescope (La Palma). We determine the stellar parameters and abundances using MARCS model atmospheres. In particular, we derive the Nb/Zr ratio which was previously shown to be a sensitive thermometer for the s-process nucleosynthesis. Indeed, in barium stars, $^{93}Zr$ has fully decayed into mono-isotopic $^{93}Nb$ , so Nb/Zr is a measure of the temperature-sensitive $^{93}Zr/Zr$ isotopic ratio. HD 28159, previously classified as K5III and initially selected to serve as a reference cool K star for our abundance analysis, turns out to be enriched in s-process elements, and as such is a new barium star. Four stars, characterised by high nitrogen abundances, also tend to have high [Nb/Zr] and [hs/ls] ratios. The derived Zr and Nb abundances provide more accurate constraints on the s-process neutron source, identified to be $^{13}C(alpha,n)^{16}O$ for barium stars. The comparison with stellar evolution and nucleosynthesis models shows that the investigated barium stars were polluted by a low-mass (2-3 Solar mass) AGB star. HD 100503 is potentially identified as the highest metallicity CEMP-rs star yet discovered.
△ Less
Submitted 17 July, 2018;
originally announced July 2018.
-
Crustal heating in accreting neutron stars from the nuclear energy-density functional theory. I. Proton shell effects and neutron-matter constraint
Authors:
A. F. Fantina,
J. L. Zdunik,
N. Chamel,
J. M. Pearson,
P. Haensel,
S. Goriely
Abstract:
Observations of soft X-ray transients in quiescence suggest the existence of heat sources in the crust of accreting neutron stars. Heat is thought to be released by electroweak and nuclear processes triggered by the burying of ashes of X-ray bursts. The heating is studied using a fully quantum approach taking consistently into account nuclear shell effects. We have followed the evolution of ashes…
▽ More
Observations of soft X-ray transients in quiescence suggest the existence of heat sources in the crust of accreting neutron stars. Heat is thought to be released by electroweak and nuclear processes triggered by the burying of ashes of X-ray bursts. The heating is studied using a fully quantum approach taking consistently into account nuclear shell effects. We have followed the evolution of ashes made of $^{56}$Fe employing the nuclear energy-density functional theory. Both the outer and inner crusts are described using the same functional, thus ensuring a unified and thermodynamically consistent treatment. To assess the role of the neutron-matter constraint, we have employed the set of accurately calibrated Brussels-Montreal functionals BSk19, BSk20, and BSk21 and for comparison the SLy4 functional. Due to nuclear shell effects, the fully accreted crust is found to be much less stratified than in previous studies. In particular, large regions of the inner crust contain clusters with the magic number $Z=14$. The heat deposited in the outer crust is tightly constrained by experimental atomic mass data. The shallow heating we obtain does not exceed $0.2$~MeV and is therefore not enough to explain the cooling of some soft X-ray transients. The total heat released in the crust is very sensitive to details of the nuclear structure and is predicted to lie in the range from $1.5$~MeV to $1.7$~MeV. The evolution of an accreted matter element and therefore the location of heat sources are governed to a large extent by the existence of nuclear shell closures. Ignoring these effects in the inner crust, the total heat falls to $\sim 0.6$~MeV. The neutron-matter constraint is also found to play a key role. The large amount of heat obtained by Steiner et al. (2012) could thus be traced back to unrealistic neutron-matter equations of state.
△ Less
Submitted 11 June, 2018;
originally announced June 2018.
-
Large-scale deformed quasiparticle random-phase approximation calculations of the $γ$-ray strength function using the Gogny force
Authors:
M. Martini,
S. Péru,
S. Hilaire,
S. Goriely,
F. Lechaftois
Abstract:
Valuable theoretical predictions of nuclear dipole excitations in the whole chart are of great interest for different nuclear applications, including in particular nuclear astrophysics. Here we present large-scale calculations of the $E1$ $γ$-ray strength function obtained in the framework of the axially-symmetric deformed QRPA based on the finite-range Gogny force. This approach is applied to eve…
▽ More
Valuable theoretical predictions of nuclear dipole excitations in the whole chart are of great interest for different nuclear applications, including in particular nuclear astrophysics. Here we present large-scale calculations of the $E1$ $γ$-ray strength function obtained in the framework of the axially-symmetric deformed QRPA based on the finite-range Gogny force. This approach is applied to even-even nuclei, the strength function for odd nuclei being derived by interpolation. The convergence with respect to the adopted number of harmonic oscillator shells and the cut-off energy introduced in the 2-quasiparticle (2-$qp$) excitation space is analyzed. The calculations performed with two different Gogny interactions, namely D1S and D1M, are compared. A systematic energy shift of the $E1$ strength is found for D1M relative to D1S, leading to a lower energy centroid and a smaller energy-weighted sum rule for D1M. When comparing with experimental photoabsorption data, the Gogny-QRPA predictions are found to overestimate the giant dipole energy by typically $\sim$2 MeV. Despite the microscopic nature of our self-consistent Hartree-Fock-Bogoliubov plus QRPA calculation, some phenomenological corrections need to be included to take into account the effects beyond the standard 2-$qp$ QRPA excitations and the coupling between the single-particle and low-lying collective phonon degrees of freedom. For this purpose, three prescriptions of folding procedure are considered and adjusted to reproduce experimental photoabsorption data at best. All of them are shown to lead to rather similar predictions of the $E1$ strength, both at low energies and for exotic neutron-rich nuclei. Predictions of $γ$-ray strength functions and Maxwellian-averaged neutron capture rates for the whole Sn isotopic chain are also discussed and compared with previous theoretical calculations.
△ Less
Submitted 28 July, 2016;
originally announced July 2016.
-
Solar r-process-constrained actinide production in neutrino-driven winds of supernovae
Authors:
S. Goriely,
H. -Th. Janka
Abstract:
Long-lived radioactive nuclei play an important role as nucleo-cosmochronometers and as cosmic tracers of nucleosynthetic source activity. In particular nuclei in the actinide region like thorium, uranium, and plutonium can testify to the enrichment of an environment by the still enigmatic astrophysical sources that are responsible for the production of neutron-rich nuclei by the rapid neutron-cap…
▽ More
Long-lived radioactive nuclei play an important role as nucleo-cosmochronometers and as cosmic tracers of nucleosynthetic source activity. In particular nuclei in the actinide region like thorium, uranium, and plutonium can testify to the enrichment of an environment by the still enigmatic astrophysical sources that are responsible for the production of neutron-rich nuclei by the rapid neutron-capture process (r-process). Supernovae and merging neutron-star (NS) or NS-black hole binaries are considered as most likely sources of the r-nuclei. But arguments in favour of one or the other or both are indirect and make use of assumptions; they are based on theoretical models with remaining simplifications and shortcomings. An unambiguous observational determination of a production event is still missing. In order to facilitate searches in this direction, e.g.\ by looking for radioactive tracers in stellar envelopes, the interstellar medium or terrestrial reservoirs, we provide improved theoretical estimates and corresponding uncertainty ranges for the actinide production ($^{232}$Th, $^{235,236,238}$U, $^{237}$Np, $^{244}$Pu, and $^{247}$Cm) in neutrino-driven winds of core-collapse supernovae. Since state-of-the-art supernova models do not yield r-process viable conditions ---but still lack, for example, the effects of strong magnetic fields--- we base our investigation on a simple analytical, Newtonian, adiabatic and steady-state wind model and consider the superposition of a large number of contributing components, whose nucleosynthesis-relevant parameters (mass weight, entropy, expansion time scale, and neutron excess) are constrained by the assumption that the integrated wind nucleosynthesis closely reproduces the solar system distribution of r-process elements. We also test the influence of uncertain nuclear physics.
△ Less
Submitted 19 April, 2016; v1 submitted 14 March, 2016;
originally announced March 2016.
-
The temperature and chronology of heavy-element synthesis in low-mass stars
Authors:
Pieter Neyskens,
Sophie Van Eck,
Alain Jorissen,
Stephane Goriely,
Lionel Siess,
Bertrand Plez
Abstract:
Roughly half of the heavy elements (atomic mass greater than that of iron) are believed to be synthesized in the late evolutionary stages of stars with masses between 0.8 and 8 solar masses. Deep inside the star, nuclei (mainly iron) capture neutrons and progressively build up (through the slow-neutron-capture process, or s-process) heavier elements that are subsequently brought to the stellar sur…
▽ More
Roughly half of the heavy elements (atomic mass greater than that of iron) are believed to be synthesized in the late evolutionary stages of stars with masses between 0.8 and 8 solar masses. Deep inside the star, nuclei (mainly iron) capture neutrons and progressively build up (through the slow-neutron-capture process, or s-process) heavier elements that are subsequently brought to the stellar surface by convection. Two neutron sources, activated at distinct temperatures, have been proposed: 13C and 22Ne, each releasing one neutron per alpha-particle (4He) captured. To explain the measured stellar abundances, stellar evolution models invoking the 13C neutron source (which operates at temperatures of about one hundred million kelvin) are favoured. Isotopic ratios in primitive meteorites, however, reflecting nucleosynthesis in the previous generations of stars that contributed material to the Solar System, point to higher temperatures (more than three hundred million kelvin), requiring at least a late activation of 22Ne. Here we report a determination of the s-process temperature directly in evolved low-mass giant stars, using zirconium and niobium abundances, independently of stellar evolution models. The derived temperature supports 13C as the s-process neutron source. The radioactive pair 93Zr-93Nb used to estimate the s-process temperature also provides, together with the pair 99Tc-99Ru, chronometric information on the time elapsed since the start of the s-process, which we determine to be one million to three million years.
△ Less
Submitted 21 January, 2016;
originally announced January 2016.
-
Detailed homogeneous abundance studies of 14 Galactic s-process enriched post-AGB stars: In search of lead (Pb)
Authors:
Kenneth De Smedt,
Hans Van Winckel,
Devika Kamath,
Lionel Siess,
Stephane Goriely,
Amanda I. Karakas,
Rajeev Manick
Abstract:
This paper is part of a larger project in which we systematically study the chemical abundances of Galactic and extragalactic post-asymptotic giant branch (post-AGB) stars. Lead (Pb) is the final product of the s-process nucleosynthesis and is predicted to have large overabundances with respect to other s-process elements in AGB stars of low metallicities. However, Pb abundance studies of enriched…
▽ More
This paper is part of a larger project in which we systematically study the chemical abundances of Galactic and extragalactic post-asymptotic giant branch (post-AGB) stars. Lead (Pb) is the final product of the s-process nucleosynthesis and is predicted to have large overabundances with respect to other s-process elements in AGB stars of low metallicities. However, Pb abundance studies of enriched post-AGB stars in the Magellanic Clouds show the Pb overabundance is not observed. We used high-resolution UVES and HERMES spectra for detailed spectral abundance studies of our sample of 14 Galactic post-AGB stars. We do not find any clear evidence of Pb overabundances in our sample. Stars with T(eff) > 7500 K do not provide strong constraints on the Pb abundance. We conclude that the discrepancy between theory and observation increases towards lower metallicities. All objects, except IRAS 17279-1119, confirm the relation between neutron exposure [hs/ls] and third dredge-up efficiency [s/Fe], whereas no relation between metallicity and neutron exposure is detected within the metallicity range of our total sample (-1.4 < [Fe/H] < -0.2). The mild enrichment of IRAS 17279-1119 can probably be attributed to a cut-off of the AGB evolution due to binary interactions. To our knowledge, IRAS 17279-1119 is the first s-process enhanced Galactic post-AGB star known in a binary system and is a possible precursor of the extrinsic Ba dwarf stars. Lead-rich stars have yet to be found in post-AGB stars. extrinsic Ba dwarf stars. Lead-rich stars are yet to be found in post-AGB stars.
△ Less
Submitted 16 December, 2015;
originally announced December 2015.
-
Nucleosynthesis in dynamical and torus ejecta of compact binary mergers
Authors:
Oliver Just,
Andreas Bauswein,
Ricard Ardevol Pulpillo,
Stephane Goriely,
H. -Thomas Janka
Abstract:
We present a comprehensive study of r-process element nucleosynthesis in the ejecta of compact binary mergers (CBMs) and their relic black-hole (BH)-torus systems. The evolution of the BH-accretion tori is simulated for seconds with a Newtonian hydrodynamics code including viscosity effects, pseudo-Newtonian gravity for rotating BHs, and an energy-dependent two-moment closure scheme for the transp…
▽ More
We present a comprehensive study of r-process element nucleosynthesis in the ejecta of compact binary mergers (CBMs) and their relic black-hole (BH)-torus systems. The evolution of the BH-accretion tori is simulated for seconds with a Newtonian hydrodynamics code including viscosity effects, pseudo-Newtonian gravity for rotating BHs, and an energy-dependent two-moment closure scheme for the transport of electron neutrinos and antineutrinos. The investigated cases are guided by relativistic double neutron star (NS-NS) and NS-BH merger models, producing ~3-6 Msun BHs with rotation parameters of A~0.8 and tori of 0.03-0.3 Msun. Our nucleosynthesis analysis includes the dynamical (prompt) ejecta expelled during the CBM phase and the neutrino and viscously driven outflows of the relic BH-torus systems. While typically ~20-25% of the initial accretion-torus mass are lost by viscously driven outflows, neutrino-powered winds contribute at most another ~1%, but neutrino heating enhances the viscous ejecta significantly. Since BH-torus ejecta possess a wide distribution of electron fractions (0.1-0.6) and entropies, they produce heavy elements from A~80 up to the actinides, with relative contributions of A>130 nuclei being subdominant and sensitively dependent on BH and torus masses and the exact treatment of shear viscosity. The combined ejecta of CBM and BH-torus phases can reproduce the solar abundances amazingly well for A>90. Varying contributions of the torus ejecta might account for observed variations of lighter elements with 40<Z<56 relative to heavier ones, and a considerable reduction of the prompt ejecta compared to the torus ejecta, e.g. in highly asymmetric NS-BH mergers, might explain the composition of heavy-element deficient stars.
△ Less
Submitted 12 October, 2015; v1 submitted 21 April, 2015;
originally announced April 2015.
-
Impact of weak interactions of free nucleons on the r-process in dynamical ejecta from neutron-star mergers
Authors:
Stephane Goriely,
Andreas Bauswein,
Oliver Just,
Else Pllumbi,
Hans-Thomas Janka
Abstract:
We investigate beta-interactions of free nucleons and their impact on the electron fraction (Y_e) and r-process nucleosynthesis in ejecta characteristic of binary neutron star mergers (BNSMs). For that we employ trajectories from a relativistic BNSM model to represent the density-temperature evolutions in our parametric study. In the high-density environment, positron captures decrease the neutron…
▽ More
We investigate beta-interactions of free nucleons and their impact on the electron fraction (Y_e) and r-process nucleosynthesis in ejecta characteristic of binary neutron star mergers (BNSMs). For that we employ trajectories from a relativistic BNSM model to represent the density-temperature evolutions in our parametric study. In the high-density environment, positron captures decrease the neutron richness at the high temperatures predicted by the hydrodynamic simulation. Circumventing the complexities of modelling three-dimensional neutrino transport, (anti)neutrino captures are parameterized in terms of prescribed neutrino luminosities and mean energies, guided by published results and assumed as constant in time. Depending sensitively on the adopted neutrino-antineutrino luminosity ratio, neutrino processes increase Y_e to values between 0.25 and 0.40, still allowing for a successful r-process compatible with the observed solar abundance distribution and a significant fraction of the ejecta consisting of r-process nuclei. If the electron neutrino luminosities and mean energies are relatively large compared to the antineutrino properties, the mean Y_e might reach values >0.40 so that neutrino captures seriously compromise the success of the r-process. In this case, the r-abundances remain compatible with the solar distribution, but the total amount of ejected r-material is reduced to a few percent, because the production of iron-peak elements is favored. Proper neutrino physics, in particular also neutrino absorption, have to be included in BNSM simulations before final conclusions can be drawn concerning r-processing in this environment and concerning observational consequences like kilonovae, whose peak brightness and color temperature are sensitive to the composition-dependent opacity of the ejecta.
△ Less
Submitted 16 April, 2015;
originally announced April 2015.
-
Cosmic Neutron Star Merger Rate and Gravitational Waves constrained by the R Process Nucleosynthesis
Authors:
E. Vangioni,
S. Goriely,
F. Daigne,
P. Francois,
K. Belczynski
Abstract:
The cosmic evolution of the neutron star merger (NSM) rate can be deduced from the observed cosmic star formation rate. This allows to estimate the rate expected in the horizon of the gravitational wave detectors advanced Virgo and ad LIGO and to compare those rates with independent predictions. In this context, the rapid neutron capture process, or r process, can be used as a constraint assuming…
▽ More
The cosmic evolution of the neutron star merger (NSM) rate can be deduced from the observed cosmic star formation rate. This allows to estimate the rate expected in the horizon of the gravitational wave detectors advanced Virgo and ad LIGO and to compare those rates with independent predictions. In this context, the rapid neutron capture process, or r process, can be used as a constraint assuming NSM is the main astrophysical site for this nucleosynthetic process. We compute the early cosmic evolution of a typical r process element, Europium. Eu yields from NSM are taken from recent nucleosynthesis calculations. The same approach allows to compute the cosmic rate of Core Collapse SuperNovae (CCSN) and the associated evolution of Eu. We find that the bulk of Eu observations at high iron abundance can be rather well fitted by either CCSN or NSM scenarios. However, at lower metallicity, the early Eu cosmic evolution favors NSM as the main astrophysical site for the r process. A comparison between our calculations and spectroscopic observations at very low metallicities allows to constrain the coalescence timescale in the NSM scenario to about 0.1 to 0.2 Gyr. These values are in agreement with the coalescence timescales of some observed binary pulsars. Finally, the cosmic evolution of Eu is used to put constraints on the NSM rate, the merger rate in the horizon of the gravitational wave detectors advanced Virgo/ad LIGO, as well as the expected rate of electromagnetic counterparts to mergers (kilonovae) in large near-infrared surveys.
△ Less
Submitted 2 October, 2015; v1 submitted 6 January, 2015;
originally announced January 2015.
-
Neutron-powered precursors of kilonovae
Authors:
Brian D. Metzger,
Andreas Bauswein,
Stephane Goriely,
Daniel Kasen
Abstract:
The merger of binary neutron stars (NSs) ejects a small quantity of neutron rich matter, the radioactive decay of which powers a day to week long thermal transient known as a kilonova. Most of the ejecta remains sufficiently dense during its expansion that all neutrons are captured into nuclei during the r-process. However, recent general relativistic merger simulations by Bauswein and collaborato…
▽ More
The merger of binary neutron stars (NSs) ejects a small quantity of neutron rich matter, the radioactive decay of which powers a day to week long thermal transient known as a kilonova. Most of the ejecta remains sufficiently dense during its expansion that all neutrons are captured into nuclei during the r-process. However, recent general relativistic merger simulations by Bauswein and collaborators show that a small fraction of the ejected mass (a few per cent, or ~1e-4 Msun) expands sufficiently rapidly for most neutrons to avoid capture. This matter originates from the shocked-heated interface between the merging NSs. Here we show that the beta-decay of these free neutrons in the outermost ejecta powers a `precursor' to the main kilonova emission, which peaks on a timescale of a few hours following merger at U-band magnitude ~22 (for an assumed distance of 200 Mpc). The high luminosity and blue colors of the neutron precursor render it a potentially important counterpart to the gravitational wave source, that may encode valuable information on the properties of the merging binary (e.g. NS-NS versus NS-black hole) and the NS equation of state. Future work is necessary to assess the robustness of the fast moving ejecta and the survival of free neutrons in the face of neutrino absorptions, although the precursor properties are robust to a moderate amount of leptonization. Our results provide additional motivation for short latency gravitational wave triggers and rapid follow-up searches with sensitive ground based telescopes.
△ Less
Submitted 8 November, 2014; v1 submitted 1 September, 2014;
originally announced September 2014.
-
Stellar origin of the 182Hf cosmochronometer and the presolar history of solar system matter
Authors:
Maria Lugaro,
Alexander Heger,
Dean Osrin,
Stephane Goriely,
Kai Zuber,
Amanda I. Karakas,
Brad K. Gibson,
Carolyn L. Doherty,
John C. Lattanzio,
Ulrich Ott
Abstract:
Among the short-lived radioactive nuclei inferred to be present in the early solar system via meteoritic analyses, there are several heavier than iron whose stellar origin has been poorly understood. In particular, the abundances inferred for 182Hf (half-life = 8.9 million years) and 129I (half-life = 15.7 million years) are in disagreement with each other if both nuclei are produced by the rapid…
▽ More
Among the short-lived radioactive nuclei inferred to be present in the early solar system via meteoritic analyses, there are several heavier than iron whose stellar origin has been poorly understood. In particular, the abundances inferred for 182Hf (half-life = 8.9 million years) and 129I (half-life = 15.7 million years) are in disagreement with each other if both nuclei are produced by the rapid neutron-capture process. Here, we demonstrate that contrary to previous assumption, the slow neutron-capture process in asymptotic giant branch stars produces 182Hf. This has allowed us to date the last rapid and slow neutron-capture events that contaminated the solar system material at roughly 100 million years and 30 million years, respectively, before the formation of the Sun.
△ Less
Submitted 9 August, 2014;
originally announced August 2014.
-
Nucleosynthesis constraints on the neutron star-black hole merger rate
Authors:
Andreas Bauswein,
Ricard Ardevol Pulpillo,
H. -Thomas Janka,
Stephane Goriely
Abstract:
We derive constraints on the time-averaged event rate of neutron star-black hole (NS-BH) mergers by using estimates of the population-integrated production of heavy rapid neutron-capture (r-process) elements with nuclear mass numbers A > 140 by such events in comparison to the Galactic repository of these chemical species. Our estimates are based on relativistic hydrodynamical simulations convolve…
▽ More
We derive constraints on the time-averaged event rate of neutron star-black hole (NS-BH) mergers by using estimates of the population-integrated production of heavy rapid neutron-capture (r-process) elements with nuclear mass numbers A > 140 by such events in comparison to the Galactic repository of these chemical species. Our estimates are based on relativistic hydrodynamical simulations convolved with theoretical predictions of the binary population. This allows us to determine a strict upper limit of the average NS-BH merger rate of ~6*10^-5 per year. We quantify the uncertainties of this estimate to be within factors of a few mostly because of the unknown BH spin distribution of such systems, the uncertain equation of state of NS matter, and possible errors in the Galactic content of r-process material. Our approach implies a correlation between the merger rates of NS-BH binaries and of double NS systems. Predictions of the detection rate of gravitational-wave signals from such compact-object binaries by Advanced LIGO and Advanced Virgo on the optimistic side are incompatible with the constraints set by our analysis.
△ Less
Submitted 25 September, 2014; v1 submitted 8 August, 2014;
originally announced August 2014.