Search

Asymptotic giant branch (AGB) stars can experience proton ingestion events (PIEs), leading to a rich nucleosynthesis. During a PIE, the intermediate neutron capture process (i-process) develops, leading to the production of trans-iron elements. It is also suggested that lithium is produced during these events. We investigate the production of lithium and trans-iron elements in AGB stars experiencing a PIE with $1

J-type stars are a subclass of carbon stars that are generally Li-rich, not enriched in s-elements, and have low $^{12}$C/$^{13}$C ratios. They were suggested to be the manufacturers of the pre-solar grains of type AB2 (having low $^{12}$C/$^{13}$C and supersolar $^{14}$N/$^{15}$N). In this Letter, we investigate the possibility that J-type stars are early asymptotic giant branch (AGB) stars that experienced a proton ingestion event (PIE). We used the stellar evolution code STAREVOL to compute AGB stellar models with initial masses of 1, 2, and 3 $M_{\odot}$ and metallicities [Fe/H] $= -0.5$ and 0.0. We included overshooting above the thermal pulse and used a network of 1160 nuclei coupled to the transport equations. In solar-metallicity AGB stars, PIEs can be triggered if a sufficiently high overshoot is considered. These events lead to low $^{12}$C/$^{13}$C ratios, high Li abundances, and no enrichment in s-elements. We find that the $2-3$ $M_{\odot}$ AGB models experiencing a PIE can account for most of the observational features of J-type stars and AB2 grains. The remaining tensions between models and observations are (1) the low $^{14}$N/$^{15}$N ratio of some AB2 grains and of 2 out of 13 J-type stars, (2) the high $^{26}$Al/$^{27}$Al of some AB2 grains, and (3) the J-type stars with A(Li) $<2$. Extra mixing mechanisms can alleviate some of these tensions, such as thermohaline or rotation. This work highlights a possible match between AGB stellar models that undergo a PIE and J-type stars and AB2 grains. To account for other types of carbon stars, such as N-type stars, PIEs should only develop in a fraction of solar-metallicity AGB stars. Additional work is needed to assess how the occurrence of PIEs depends on mixing parameters and initial conditions, and therefore to further confirm or exclude the proposed scenario., Comment: Accepted in A&A

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

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

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., Comment: 21 pages, 15 figures, accepted in A&A

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

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)., Comment: 17 pages, 18 figures Preview and access to the library: https://www.royer.se/melchiors.html

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., Comment: Accepted: October 11, 2023 \\ 14 Pages, 14 Figures, 2 Tables

Proton-$\gamma$ 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 $\gamma$-ray strength function ($\gamma$SF) of $^{67}\mathrm{Ni}$. The $^{66}\mathrm{Ni}(n,\gamma)$ capture cross section has been calculated using the Hauser-Feshbach model in TALYS using the measured NLD and $\gamma$SF as constraints. The results confirm that the $^{66}\mathrm{Ni}(n,\gamma)$ reaction acts as a bottleneck when relying on one-zone nucleosynthesis calculations. However, the impact of this reaction is strongly dampened in multi-zone models of low-metallicity AGB stars experiencing i-process nucleosynthesis., Comment: Submitted to Phys. Rev. C

Binary neutron star mergers are the first confirmed site of rapid neutron capture (r-process) element nucleosynthesis. 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.08um 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.07um 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.07um emission feature is best reproduced with two lines, centred at 2.059 and 2.135um. 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., Comment: Main text: 19 pages, 7 figures, 6 tables. Appendices: 10 pages, 3 tables. Submitted to MNRAS. Comments welcome

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., Comment: Accepted for publication in A&A, 9 pages, 9 figures including the two in appendix

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

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., Comment: 20 pages, 13 figures. Accepted by MNRAS

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

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., Comment: 8 pages, 4 figures; Contribution to the proceedings of INPC 2022, Cape Town, South Africa

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.

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., Comment: 6 pages, 3 figures, accepted in A&A

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

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

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 $\Gamma \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., Comment: 10 pages, 7 figures. The tables of the equations of state are available at the CDS

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., Comment: 18 pages, 16 figures, accepted in A&A

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

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 $\alpha-$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 $\alpha-$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., Comment: Accepted in A&A

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, $\beta$-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., Comment: 28 pages, 24 figures, accepted to MNRAS

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($\alpha$,$\gamma$)$^{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($\alpha,\gamma$) 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., Comment: 12 pages, 13 figures, 2 tables, accepted for publication in A&A

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

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$\,$\AA, 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., Comment: 20 pages, 14 figures and 3 tables. Accepted for publication in MNRAS

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

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

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., Comment: 18 pages, 21 figures, 2 Tables, accepted by MNRAS 2021-11-19

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, Comment: 9 pages, 8 Figures, A&A to be published

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)., Comment: Accepted for publication in A&A, 15 pages, 12 figures

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., Comment: Published in A&A. This version of the preprint has been corrected (especially Table 2 and Section 6.5) with the modifications appearing in the A&A corrigendum (which is included at the end of the paper)

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

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