Search

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

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

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

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., Comment: 16 pages, 9 figures, 1 table, accepted to ApJL

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

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

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., Comment: Accepted for publication in the Astrophysical Journal Letters, 15 pages, 4 figures, 1 table

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

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

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., Comment: 6 pages, 3 figures, Under review for 'The origin of outflows in evolved stars Proceedings IAU Symposium No. 366, 2021'

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., Comment: 23 pages, 15 figures, 1 table, accepted to MNRAS; updated for correct version of Fig. A1, left panel, and corrected Eqs. 7 and 9

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., Comment: 22 pages, 10 figures, accepted for publication in Section 7. Stellar structure and evolution of Astronomy and Astrophysics (A&A)

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., Comment: 37 pages, 22 figures, 5 tables, accepted to MNRAS, minor changes compared to previous version

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

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., Comment: 16 pages, 19 figures. Modified to match the version accepted by MNRAS (some changes to main results and conclusions, new section added)

We present recent shell-model calculations of the gamma-decay in sd-pf and pf-shell nuclei. We focus on the M1 part of the dipole strength which was shown to exhibit interesting low-energy effects, in particular a low-energy enhancement which can have a considerable impact on the radiative neutron capture. We discuss the persistence of the shell effects in the nuclear quasi-continuum and the relation between the shape of the strength function at low energy and nuclear deformation., Comment: 8 pages, 4 figures

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., Comment: 15 pages, 6 figures. To appear in the AIP Conference Proceedings of the Xiamen-CUSTIPEN Workshop on the EOS of Dense Neutron-Rich Matter in the Era of Gravitational Wave Astronomy (January 3 - 7, 2019, Xiamen, China)

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., Comment: Recommended for publication in A&A letters

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., Comment: 14 pages including 3 figures

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., Comment: Proceedings IAU Symposium No. 343: Why galaxies care about AGB stars?, 2018

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., Comment: 23 pages, 19 figures

The nuclear level densities and $\gamma$-ray strength functions of $^{138,139,140}$La were measured using the $^{139}$La($^{3}$He, $\alpha$), $^{139}$La($^{3}$He, $^{3}$He$^\prime$) and $^{139}$La(d, p) reactions. The particle-$\gamma$ coincidences were recorded with the silicon particle telescope (SiRi) and NaI(Tl) (CACTUS) arrays. In the context of these experimental results, the low-energy enhancement in the A$\sim$140 region is discussed. The $^{137,138,139}$La($n, \gamma)$ cross sections were calculated at $s$- and $p$-process temperatures using the experimentally measured nuclear level densities and $\gamma$-ray strength functions. Good agreement is found between $^{139}$La($n, \gamma)$ calculated cross sections and previous measurements.

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., Comment: 30 pages, 10 figures

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., Comment: 26 pages, 40 figures, Accepted for publication in Astronomy and Astrophysics, main journal

A survey of some results in the modern microscopic theory of properties of nuclear reactions with gamma-rays is given. First of all, we discuss the impact of phonon coupling (PC) on the photon strength function (PSF) because it represents the most natural physical source of additional strength found for Sn isotopes in recent experiments that could not be explained within the stan- dard HFB+QRPA approach. The self-consistent version of the Extended Theory of Finite Fermi Systems in the Quasiparticle Time Blocking Approximation, or simply QTBA, is applied. It uses the HFB mean field and includes both the QRPA and PC effects on the basis of the SLy4 Skyrme force. With our microscopic E1 PSFs, the following properties have been calculated for many stable and unstable even-even semi-magic Sn and Ni isotopes as well as for double-magic 132Sn and 208Pb using the reaction codes EMPIRE and TALYS with several nuclear level density (NLD) models: 1) the neutron capture cross sections, 2) the corresponding neutron capture gamma spectra, 3) the av- erage radiative widths of neutron resonances. In all the properties considered, the PC contribution turned out to be significant, as compared with the standard QRPA one, and necessary to explain the available experimental data. The results with the phenomenological so-called generalized super- fluid NLD model turned out to be worse, on the whole, than those obtained with the microscopic HFB+combinatorial NLD model. Finally, we also discuss the modern microscopic NLD models based on the self-consistent HFB method and show their relevance to explain experimental data as compared with the phenomeno- logical models. The use of these self-consistent microscopic approaches is of particular relevance for nuclear astrophysics, but also for the study of double-magic nuclei., Comment: 13 pages, 14 figures, a survey given as a plenary talk to the Intern. Conference "NUCLEUS 2015" (June 29 - July 3, 2015, Saint-Petersburg, Russia). To be published in Phys. Atom. Nucl

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

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., Comment: 12 pages, 9 figures; submitted to MNRAS

The pygmy dipole resonances and photon strength functions in stable and unstable Ni and Sn isotopes are calculated within the microscopic self-consistent version of the extended theory of finite fermi systems which includes phonon coupling effects.The Skyrme forces SLy4 is used. A pygmy dipole resonance in $^{72}$Ni is predicted at the mean energy of 12.4 MeV exhausting 25.7\% of the total energy-weighted sum rule. The microscopically obtained photon E1 strength functions are compared with available experimental data and used to calculate nuclear reaction properties. Average radiative widths and radiative neutron capture cross sections have been calculated taking the phonon coupling into account as well as the uncertainties caused by various microscopic level density models. In all three quantities considered, the contribution of phonon coupling turned out to be significant and is found necessary to explain available experimental data., Comment: 4 pages, 3 figures Submitted to Phys. Rev. C

The pygmy-dipole resonances and photon strength functions in stable and unstable Ni and Sn isotopes are calculated within the microscopic self-consistent version of the extended theory of finite fermi systems which includes the QRPA and phonon coupling effects and uses the known Skyrme forces SLy4. The pygmy dipole resonance in $^{72}Ni$ is predicted with the mean energy of 12.4 MeV and the energy-weighted sum rule exhausting 25.6\% of the total strength. The microscopically obtained photon E1 strength functions are used to calculate nuclear reaction properties, i.e the radiative neutron capture cross section, gamma-ray spectra, and average radiative widths. Our main conclusion is that in all these quantities it is necessary to take the phonon coupling effects into account., Comment: 4 pages, 5 figures, 2 tables. Talk at 15-th International Symposium on Capture Gamma-Ray Spectroscopy and Related Topics (CGS15), Dresden, August 2o14

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., Comment: 6 pages, 5 figures, accepted to MNRAS main journal

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., Comment: 24 pages, 5 figures, published on Science Vol. 345, pp. 650-653

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., Comment: 5 pages, 3 figures; accepted for publication in ApJL

We present the first 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

This "Comment" submitted to ADNDT on December 13, 2013 concerns a publication entitled "Atomic Mass Compilation 2012", which appeared in the March 2014 issue of the journal Atomic Data and Nuclear Data Tables (available online on September 6, 2013)., Comment: 4 pages

An update of the NACRE compilation [Angulo et al., Nucl. Phys. A 656 (1999) 3] is presented. This new compilation, referred to as NACRE II, reports thermonuclear reaction rates for 34 charged-particle induced, two-body exoergic reactions on nuclides with mass number $A<16$, of which fifteen are particle-transfer reactions and the rest radiative capture reactions. When compared with NACRE, NACRE II features in particular (1) the addition to the experimental data collected in NACRE of those reported later, preferentially in the major journals of the field by early 2013, and (2) the adoption of potential models as the primary tool for extrapolation to very low energies of astrophysical $S$-factors, with a systematic evaluation of uncertainties. As in NACRE, the rates are presented in tabular form for temperatures in the $10^{6}$ $\simeq\leq$ T $\leq$ $10^{10}$ K range. Along with the 'adopted' rates, their low and high limits are provided. The new rates are available in electronic form as part of the Brussels Library (BRUSLIB) of nuclear data. The NACRE II rates also supersede the previous NACRE rates in the Nuclear Network Generator (NETGEN) for astrophysics. [http://www.astro.ulb.ac.be/databases.html.], Comment: 86 figures