BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

S. Goriely

S. Goriely appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
nucl-thastro-ph.HEastro-ph.SRnucl-exastro-phcond-mat.supr-conphysics.ins-det

Catalog footprint

What is connected

58works
7topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

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

Published work

58 published item(s)

preprint2026arXiv

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

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

0 citations0 reviewsNo code linkedOpen access
preprint2026arXiv

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

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

0 citations0 reviewsNo code linkedOpen access
preprint2026arXiv

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

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

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

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

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

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

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

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

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

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

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

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

The p-process in exploding rotating massive stars

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

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

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

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

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

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

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

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

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

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

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

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

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

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

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

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

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

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

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

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Experimentally constrained ($p,γ$)$^{89}$Y and ($n,γ$)$^{89}$Y reaction rates relevant to the $p$-process nucleosynthesis

The nuclear level density and the $γ$-ray strength function have been extracted for $^{89}$Y, using the Oslo Method on $^{89}$Y($p,p' γ$)$^{89}$Y coincidence data. The $γ$-ray strength function displays a low-energy enhancement consistent with previous observations in this mass region ($^{93-98}$Mo). Shell-model calculations give support that the observed enhancement is due to strong, low-energy $M1$ transitions at high excitation energies. The data were further used as input for calculations of the $^{88}$Sr($p,γ$)$^{89}$Y and $^{88}$Y($n,γ$)$^{89}$Y cross sections with the TALYS reaction code. Comparison with cross-section data, where available, as well as with values from the BRUSLIB library, shows a satisfying agreement.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Large-scale deformed quasiparticle random-phase approximation calculations of the $γ$-ray strength function using the Gogny force

Valuable theoretical predictions of nuclear dipole excitations in the whole chart are of great interest for different nuclear applications, including in particular nuclear astrophysics. Here we present large-scale calculations of the $E1$ $γ$-ray strength function obtained in the framework of the axially-symmetric deformed QRPA based on the finite-range Gogny force. This approach is applied to even-even nuclei, the strength function for odd nuclei being derived by interpolation. The convergence with respect to the adopted number of harmonic oscillator shells and the cut-off energy introduced in the 2-quasiparticle (2-$qp$) excitation space is analyzed. The calculations performed with two different Gogny interactions, namely D1S and D1M, are compared. A systematic energy shift of the $E1$ strength is found for D1M relative to D1S, leading to a lower energy centroid and a smaller energy-weighted sum rule for D1M. When comparing with experimental photoabsorption data, the Gogny-QRPA predictions are found to overestimate the giant dipole energy by typically $\sim$2 MeV. Despite the microscopic nature of our self-consistent Hartree-Fock-Bogoliubov plus QRPA calculation, some phenomenological corrections need to be included to take into account the effects beyond the standard 2-$qp$ QRPA excitations and the coupling between the single-particle and low-lying collective phonon degrees of freedom. For this purpose, three prescriptions of folding procedure are considered and adjusted to reproduce experimental photoabsorption data at best. All of them are shown to lead to rather similar predictions of the $E1$ strength, both at low energies and for exotic neutron-rich nuclei. Predictions of $γ$-ray strength functions and Maxwellian-averaged neutron capture rates for the whole Sn isotopic chain are also discussed and compared with previous theoretical calculations.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Solar r-process-constrained actinide production in neutrino-driven winds of supernovae

Long-lived radioactive nuclei play an important role as nucleo-cosmochronometers and as cosmic tracers of nucleosynthetic source activity. In particular nuclei in the actinide region like thorium, uranium, and plutonium can testify to the enrichment of an environment by the still enigmatic astrophysical sources that are responsible for the production of neutron-rich nuclei by the rapid neutron-capture process (r-process). Supernovae and merging neutron-star (NS) or NS-black hole binaries are considered as most likely sources of the r-nuclei. But arguments in favour of one or the other or both are indirect and make use of assumptions; they are based on theoretical models with remaining simplifications and shortcomings. An unambiguous observational determination of a production event is still missing. In order to facilitate searches in this direction, e.g.\ by looking for radioactive tracers in stellar envelopes, the interstellar medium or terrestrial reservoirs, we provide improved theoretical estimates and corresponding uncertainty ranges for the actinide production ($^{232}$Th, $^{235,236,238}$U, $^{237}$Np, $^{244}$Pu, and $^{247}$Cm) in neutrino-driven winds of core-collapse supernovae. Since state-of-the-art supernova models do not yield r-process viable conditions ---but still lack, for example, the effects of strong magnetic fields--- we base our investigation on a simple analytical, Newtonian, adiabatic and steady-state wind model and consider the superposition of a large number of contributing components, whose nucleosynthesis-relevant parameters (mass weight, entropy, expansion time scale, and neutron excess) are constrained by the assumption that the integrated wind nucleosynthesis closely reproduces the solar system distribution of r-process elements. We also test the influence of uncertain nuclear physics.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Analytical representations of unified equations of state for neutron-star matter

Context. A unified equation of state (EoS) should describe the crust and the core of a neutron star using the same physical model. The Brussels-Montreal group has recently derived a family of such EoSs based on the nuclear energy-density functional theory with generalized Skyrme effective forces, fitted to the available mass data. At the same time, these forces were constrained to reproduce microscopic calculations of homogeneous neutron matter based on realistic two- and three-nucleon forces. Aims. We represent basic physical characteristics of the latest Brussels-Montreal EoS models by analytical expressions to facilitate their inclusion in astrophysical simulations. Methods. We consider three EoS models, which significantly differ by stiffness: BSk19, BSk20, and BSk21. For each of them we constructed two versions of the EoS parametrization. In the first version, pressure P and gravitational mass density ρare given as functions of the baryon number density n_b. In the second version, P, ρ, and n_b are given as functions of pseudo-enthalpy, which is useful for two-dimensional calculations of stationary rotating configurations of neutron stars. In addition to the EoS, we derived analytical expressions for several related quantities that are required in neutron-star simulations: number fractions of electrons and muons in the stellar core, nucleon numbers per nucleus in the inner crust, and equivalent radii and shape parameters of the nuclei in the inner crust. Results. We obtain analytical representations for the basic characteristics of the models of cold dense matter, which are most important for studies of neutron stars. We demonstrate the usability of our results by applying them to calculations of neutron-star mass-radius relations, maximum and minimum masses, thresholds of direct Urca processes, and the electron conductivity in the neutron-star crust.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Cosmic Neutron Star Merger Rate and Gravitational Waves constrained by the R Process Nucleosynthesis

The cosmic evolution of the neutron star merger (NSM) rate can be deduced from the observed cosmic star formation rate. This allows to estimate the rate expected in the horizon of the gravitational wave detectors advanced Virgo and ad LIGO and to compare those rates with independent predictions. In this context, the rapid neutron capture process, or r process, can be used as a constraint assuming NSM is the main astrophysical site for this nucleosynthetic process. We compute the early cosmic evolution of a typical r process element, Europium. Eu yields from NSM are taken from recent nucleosynthesis calculations. The same approach allows to compute the cosmic rate of Core Collapse SuperNovae (CCSN) and the associated evolution of Eu. We find that the bulk of Eu observations at high iron abundance can be rather well fitted by either CCSN or NSM scenarios. However, at lower metallicity, the early Eu cosmic evolution favors NSM as the main astrophysical site for the r process. A comparison between our calculations and spectroscopic observations at very low metallicities allows to constrain the coalescence timescale in the NSM scenario to about 0.1 to 0.2 Gyr. These values are in agreement with the coalescence timescales of some observed binary pulsars. Finally, the cosmic evolution of Eu is used to put constraints on the NSM rate, the merger rate in the horizon of the gravitational wave detectors advanced Virgo/ad LIGO, as well as the expected rate of electromagnetic counterparts to mergers (kilonovae) in large near-infrared surveys.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Experimental constraints on the $γ$-ray strength function in $^{90}$Zr using partial cross sections of the $^{89}$Y(p,$γ$)$^{90}$Zr reaction

Partial cross sections of the $^{89}$Y(p,$γ$)$^{90}$Zr reaction have been measured to investigate the $γ$-ray strength function in the neutron-magic nucleus $^{90}$Zr. For five proton energies between $E_p=3.65$ MeV and $E_p=4.70$ MeV, partial cross sections for the population of seven discrete states in $^{90}$Zr have been determined by means of in-beam $γ$-ray spectroscopy. Since these $γ$-ray transitions are dominantly of $E1$ character, the present measurement allows an access to the low-lying dipole strength in $^{90}$Zr. A $γ$-ray strength function based on the experimental data could be extracted, which is used to describe the total and partial cross sections of this reaction by Hauser-Feshbach calculations successfully. Significant differences with respect to previously measured strength functions from photoabsorption data point towards deviations from the Brink-Axel hypothesis relating the photo-excitation and de-excitation strength functions.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Low-energy enhancement in the γ-ray strength functions of $^{73,74}$Ge

The $γ$-ray strength functions and level densities of $^{73,74}$Ge have been extracted up to the neutron separation energy S$_n$ from particle-$γ$ coincidence data using the Oslo method. Moreover, the $γ$-ray strength function of $^{74}$Ge above S$_n$ has been determined from photo-neutron measurements, hence these two experiments cover the range of E$_γ\approx$ 1-13 MeV for $^{74}$Ge. The obtained data show that both $^{73,74}$Ge display an increase in strength at low $γ$ energies. The experimental $γ$-ray strength functions are compared with $M1$ strength functions deduced from average $B(M1)$ values calculated within the shell model for a large number of transitions. The observed low-energy enhancements in $^{73,74}$Ge are adopted in the calculations of the $^{72,73}$Ge(n,$γ$) cross sections, where there are no direct experimental data. Calculated reaction rates for more neutron-rich germanium isotopes are shown to be strongly dependent on the presence of the low-energy enhancement.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Gamow-Teller strength in deformed nuclei within the self-consistent charge-exchange quasi-particle random-phase approximation with the Gogny force

The charge-exchange excitations in nuclei are studied within the fully self-consistent proton-neutron quasiparticle random-phase approximation using the finite-range Gogny interaction. No additional parameters beyond those included in the effective nuclear force are included. Axially symmetric deformations are consistently taken into account, both in the description of the ground states and spin-isospin excitations. We focus on the isobaric analog and Gamow-Teller resonances. A comparison of the predicted strength distributions to the existing experimental data is presented and the role of nuclear deformation analyzed. The Gamow-Teller strength is used to estimate the beta-decay half-life of nuclei for which experimental data exist. A satisfactory agreement with experimental half-lives is found and justifies the additional study of the exotic neutron-rich N=82, 126 and 184 isotonic chains of relevance for the r-process nucleosynthesis.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Linear response theory and neutrino mean free path using Brussels-Montreal Skyrme functionals

The Brussels-Montreal Skyrme functionals have been successful to describe properties of both finite nuclei and infinite homogeneous nuclear matter. In their latest version, these functionals have been equipped with two extra density-dependent terms in order to reproduce simultaneously ground state properties of nuclei and infinite nuclear matter properties while avoiding at the same time the arising of ferromagnetic instabilities. In the present article, we extend our previous results of the linear response theory to include such extra terms at both zero and finite temperature in pure neutron matter. The resulting formalism is then applied to derive the neutrino mean free path. The predictions from the Brussels-Montreal Skyrme functionals are compared with ab-initio methods.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

The lead discrepancy in intrinsically s-process enriched post-AGB stars in the Magellanic Clouds

Context: Our understanding of the s-process nucleosynthesis in asymptotic giant branch (AGB) stars is incomplete. AGB models predict, for example, large overabundances of lead (Pb) compared to other s-process elements in metal-poor low-mass AGB stars. This is indeed observed in some extrinsically enhanced metal-poor stars, but not in all. An extensive study of intrinsically s-process enriched objects is essential for improving our knowledge of the AGB third dredge-up and associated s-process nucleosynthesis. Aims: We compare the spectral abundance analysis of the SMC post-AGB star J004441.04-732136.4 with state-of-the-art AGB model predictions with a main focus on Pb. The low S/N in the Pb line region made the result of our previous study inconclusive. We acquired additional data covering the region of the strongest Pb line. Methods: By carefully complementing re-reduced previous data, with newly acquired UVES optical spectra, we improve the S/N of the spectrum around the strongest Pb line. Therefore, an upper limit for the Pb abundance is estimated from a merged weighted mean spectrum using synthetic spectral modeling. We then compare the abundance results from the combined spectra to predictions of tailored AGB evolutionary models from two independent evolution codes. In addition, we determine upper limits for Pb abundances for three previously studied LMC post-AGB objects. Results: Although theoretical predictions for J004441.04-732136.4 match the s-process distribution up to tungsten (W), the predicted very high Pb abundance is clearly not detected. The three additional LMC post-AGB stars show a similar lack of a very high Pb abundance. Conclusion: From our study, we conclude that none of these low-mass, low-metallicity post-AGB stars of the LMC and SMC are strong Pb producers. This conflicts with current theoretical predictions.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Cross-Section Measurements of the 86Kr(g,n) Reaction to Probe the s-Process Branching at 85Kr

We have carried out photodisintegration cross-section measurements on 86Kr using monoenergetic photon beams ranging from the neutron separation energy, S_n = 9.86 MeV, to 13 MeV. We combine our experimental 86Kr(g,n)85Kr cross section with results from our recent 86Kr(g,g') measurement below the neutron separation energy to obtain the complete nuclear dipole response of 86Kr. The new experimental information is used to predict the neutron capture cross section of 85Kr, an important branching point nucleus on the abundance flow path during s-process nucleosynthesis. Our new and more precise 85Kr(n,g)86Kr cross section allows to produce more precise predictions of the 86Kr abundance from s-process models. In particular, we find that the models of the s-process in asymptotic giant branch stars of mass < 1.5 Msun, where the 13C neutron source burns convectively rather than radiatively, represent a possible solution for the highest 86Kr/82Kr ratios observed in meteoritic stardust SiC grains.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Does a hadron-quark phase transition in dense matter preclude the existence of massive neutron stars ?

We study the impact of a hadron-quark phase transition on the maximum neutron-star mass. The hadronic part of the equation of state relies on the most up-to-date Skyrme nuclear energy density functionals, fitted to essentially all experimental nuclear mass data and constrained to reproduce the properties of infinite nuclear matter as obtained from microscopic calculations using realistic forces. We show that the softening of the dense matter equation of state due to the phase transition is not necessarily incompatible with the existence of massive neutron stars like PSR J1614-2230.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

New fission fragment distributions and r-process origin of the rare-earth elements

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, but calculations sensitively depend on nuclear physics. In particular the fission fragment yields determine the creation of 110 < A < 170 nuclei. Here we apply a new scission-point model, called SPY, to derive the fission fragment distribution (FFD) of all relevant neutron-rich, fissioning nuclei. The model predicts a doubly asymmetric FFD in the abundant A ~ 278 mass region that is responsible for the final recycling of the fissioning material. Using ejecta conditions based on relativistic NS merger calculations we show that this specific FFD leads to a production of the A ~ 165 rare-earth peak that is nicely compatible with the abundance patterns in the Sun and metal-poor stars. This new finding further strengthens the case of NS mergers as possible dominant origin of r-nuclei with A > 140.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Nuclear charge radii of potassium isotopes beyond N=28

We report on the measurement of optical isotope shifts for $^{38,39,42,44,46\text{-}51}$K relative to $^{47}$K from which changes in the nuclear mean square charge radii across the N=28 shell closure are deduced. The investigation was carried out by bunched-beam collinear laser spectroscopy at the CERN-ISOLDE radioactive ion-beam facility. Mean square charge radii are now known from $^{37}$K to $^{51}$K, covering all $νf_{7/2}$-shell as well as all $νp_{3/2}$-shell nuclei. These measurements, in conjunction with those of Ca, Cr, Mn and Fe, provide a first insight into the $Z$ dependence of the evolution of nuclear size above the shell closure at N=28.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Phase transitions in dense matter and the maximum mass of neutron stars

The recent precise measurement of the mass of pulsar PSR J1614$-$2230, as well as observational indications of even more massive neutron stars, has revived the question of the composition of matter at the high densities prevailing inside neutron-star cores. We study the impact on the maximum possible neutron-star mass of an "exotic" core consisting of non-nucleonic matter. For this purpose, we study the occurrence of a first-order phase transition in nucleonic matter. Given the current lack of knowledge of non-nucleonic matter, we consider the stiffest possible equation of state subject only to the constraints of causality and thermodynamic stability. The case of a hadron-quark phase transition is discussed separately. The purely nucleonic matter is described using a set of unified equations of state that have been recently developed to permit a consistent treatment of both homogeneous and inhomogeneous phases. We then compute the mass-radius relation of cold nonaccreting neutron stars with and without exotic cores from the Tolman-Oppenheimer-Volkoff equations. We find that even if there is a significant softening of the equation of state associated with the actual transition to an exotic phase, there can still be a stiffening at higher densities closer to the center of the star that is sufficient to increase the maximum possible mass. However with quarks the maximum neutron-star mass is always reduced by assuming that the sound speed is limited by $c/\sqrt{3}$ as suggested by QCD calculations. In particular, by invoking such a phase transition, it becomes possible to support PSR J1614$-$2230 with a nucleonic equation of state that is soft enough to be compatible with the kaon and pion production in heavy-ion collisions.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Properties of the outer crust of strongly magnetized neutron stars from Hartree-Fock-Bogoliubov atomic mass models

The equilibrium properties of the outer crust of cold nonaccreting magnetars (i.e. neutron stars endowed with very strong magnetic fields) are studied using the latest experimental atomic mass data complemented with a microscopic atomic mass model based on the Hartree-Fock-Bogoliubov method. The Landau quantization of electron motion caused by the strong magnetic field is found to have a significant impact on the composition and the equation of state of crustal matter. It is also shown that the outer crust of magnetars could be much more massive than that of ordinary neutron stars.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

STARLIB: A Next-Generation Reaction-Rate Library for Nuclear Astrophysics

STARLIB is a next-generation, all-purpose nuclear reaction-rate library. For the first time, this library provides the rate probability density at all temperature grid points for convenient implementation in models of stellar phenomena. The recommended rate and its associated uncertainties are also included. Currently, uncertainties are absent from all other rate libraries, and, although estimates have been attempted in previous evaluations and compilations, these are generally not based on rigorous statistical definitions. A common standard for deriving uncertainties is clearly warranted. STARLIB represents a first step in addressing this deficiency by providing a tabular, up-to-date database that supplies not only the rate and its uncertainty but also its distribution. Because a majority of rates are lognormally distributed, this allows the construction of rate probability densities from the columns of STARLIB. This structure is based on a recently suggested Monte Carlo method to calculate reaction rates, where uncertainties are rigorously defined. In STARLIB, experimental rates are supplemented with: (i) theoretical TALYS rates for reactions for which no experimental input is available, and (ii) laboratory and theoretical weak rates. STARLIB includes all types of reactions of astrophysical interest to Z = 83, such as (p,g), (p,a), (a,n), and corresponding reverse rates. Strong rates account for thermal target excitations. Here, we summarize our Monte Carlo formalism, introduce the library, compare methods of correcting rates for stellar environments, and discuss how to implement our library in Monte Carlo nucleosynthesis studies. We also present a method for accessing STARLIB on the Internet and outline updated Monte Carlo-based rates.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Symmetry energy: nuclear masses and neutron stars

We describe the main features of our most recent Hartree-Fock-Bogoliubov nuclear mass models, based on 16-parameter generalized Skyrme forces. They have been fitted to the data of the 2012 Atomic Mass Evaluation, and favour a value of 30 MeV for the symmetry coefficient J, the corresponding root-mean square deviation being 0.549 MeV. We find that this conclusion is compatible with measurements of neutron-skin thickness. By constraining the underlying interactions to fit various equations of state of neutron matter calculated {\it ab initio} our models are well adapted to a realistic and unified treatment of all regions of neutron stars. We use our models to calculate the composition, the equation of state, the mass-radius relation and the maximum mass. Comparison with observations of neutron stars again favours a value of J = 30 MeV.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Systematics of dynamical mass ejection, nucleosynthesis, and radioactively powered electromagnetic signals from neutron-star mergers

We investigate systematically the dynamical mass ejection, r-process nucleosynthesis, and properties of electromagnetic counterparts of neutron-star (NS) mergers in dependence on the uncertain properties of the nuclear equation of state (EoS) by employing 40 representative, microphysical high-density EoSs in relativistic, hydrodynamical simulations. The crucial parameter determining the ejecta mass is the radius R_1.35 of a 1.35 M_sun NS. NSs with smaller R_1.35 ("soft" EoS) eject systematically higher masses. These range from ~10^-3 M_sun to ~10^-2 M_sun for 1.35-1.35 M_sun binaries and from ~5*10^-3 M_sun to ~2*10^-2 M_sun for 1.2-1.5 M_sun systems (with kinetic energies between ~5*10^49 erg and 10^51 erg). Correspondingly, the bolometric peak luminosities of the optical transients of symmetric (asymmetric) mergers vary between 3*10^41 erg/s and 14*10^41 erg/s (9*10^41 erg/s and 14.5*10^41 erg/s) on timescales between ~2 h and ~12 h. If these signals with absolute bolometric magnitudes from -15.0 to -16.7 are measured, the tight correlation of their properties with those of the merging NSs might provide valuable constraints on the high-density EoS. The r-process nucleosynthesis exhibits a remarkable robustness independent of the EoS, producing a nearly solar abundance pattern above mass number 130. By the r-process content of the Galaxy and the average production per event the Galactic merger rate is limited to 4*10^-5/yr (4*10^-4/yr) for a soft (stiff) NS EoS, if NS mergers are the main source of heavy r-nuclei. The production ratio of radioactive 232Th to 238U attains a stable value of 1.64-1.67, which does not exclude NS mergers as potential sources of heavy r-material in the most metal-poor stars.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Transitional $γ$ strength in Cd isotopes

The level densities and $γ$-ray strength functions of $^{105,106,111,112}$Cd have been extracted from particle-$γ$ coincidence data using the Oslo method. The level densities are in very good agreement with known levels at low excitation energy. The $γ$-ray strength functions display no strong enhancement for low $γ$ energies. However, more low-energy strength is apparent for $^{105,106}$Cd than for $^{111,112}$Cd. For $γ$ energies above $\approx$ 4 MeV, there is evidence for some extra strength, similar to what has been previously observed for the Sn isotopes. The origin of this extra strength is unclear; it might be due to $E1$ and $M1$ transitions originating from neutron skin oscillations or the spin-flip resonance, respectively.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Unified description of dense matter in neutron stars and magnetars

We have recently developed a set of equations of state based on the nuclear energy density functional theory providing a unified description of the different regions constituting the interior of neutron stars and magnetars. The nuclear functionals, which were constructed from generalized Skyrme effective nucleon-nucleon interactions, yield not only an excellent fit to essentially all experimental atomic mass data but were also constrained to reproduce the neutron-matter equation of state as obtained from realistic many-body calculations.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Equation of state of magnetar crusts from Hartree-Fock-Bogoliubov atomic mass models

The equation of state (EoS) of the outer crust of a cold non-accreting magnetar has been determined using the model of Lai and Shapiro (1991). For this purpose, we have made use of the latest experimental atomic mass data complemented with a Hartree-Fock-Bogoliubov (HFB) mass model. Magnetar crusts are found to be significantly different from the crusts of ordinary neutron stars.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Impact of a low-energy enhancement in the gamma-ray strength function on the radiative neutron-capture

A low-energy enhancement of the gamma-ray strength function in several light and medium-mass nuclei has been observed recently in 3He-induced reactions. The effect of this enhancement on (n,gamma) cross-sections is investigated for stable and unstable neutron-rich Fe, Mo and Cd isotopes. Our results indicate that the radiative neutron capture cross sections may increase considerably due to the low-energy enhancement when approaching the neutron drip line. This could have non-negligible consequences on r-process nucleosynthesis calculations.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Inner crust of neutron stars with mass-fitted Skyrme functionals

The equation of state and composition of the inner crust of neutron stars at zero temperature are calculated, using the T = 0 version of the TETFSI (temperature-dependent extended Thomas-Fermi plus Strutinsky integral) method, for each of a family of three functionals based on Skyrme-type forces BSk19, BSk20 and BSk21, which are characterized by different degrees of symmetry-energy stiffness, and also for the SLy4 functional. We also solve the Tolman-Oppenheimer-Volkoff equations to calculate the distribution of mass within the inner crust. Qualitatively similar results are found for all four functionals, and in particular the number of protons per Wigner-Seitz cell is in all cases equal to 40 throughout the inner crust.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Large-scale Continuum Random Phase Approximation predictions of dipole strength for astrophysical applications

Large-scale calculations of the E1 strength are performed within the random phase approximation (RPA) based on the relativistic point-coupling mean field approach in order to derive the radiative neutron capture cross sections for all nuclei of astrophysical interest. While the coupling to the single-particle continuum is taken into account in an explicit and self-consistent way, additional corrections like the coupling to complex configurations and the temperature and deformation effects are included in a phenomenological way to account for a complete description of the nuclear dynamical problem. It is shown that the resulting E1-strength function based on the PCF1 force is in close agreement with photoabsorption data as well as the available experimental E1 strength data at low energies. For neutron-rich nuclei, as well as light neutron-deficient nuclei, a low-lying so-called pygmy resonance is found systematically in the 5-10 MeV region. The corresponding strength can reach 10% of the giant dipole strength in the neutron-rich region and about 5% in the neutron-deficient region, and is found to be reduced in the vicinity of the shell closures. Finally, the neutron capture reaction rates of neutron-rich nuclei is found to be about 2-5 times larger than those predicted on the basis of the nonrelativistic RPA calculation and about a factor 50 larger than obtained with traditional Lorentzian-type approaches.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Primary gamma-ray spectra in 44Ti of astrophysical interest

Primary gamma-ray spectra for a wide excitation-energy range have been extracted for 44Ti from particle-gamma coincidence data of the 46Ti(p,t gamma)44Ti reaction. These spectra reveal information on the gamma-decay pattern of the nucleus, and may be used to extract the level density and radiative strength function applying the Oslo method. Models of the level density and radiative strength function are used as input for cross-section calculations of the 40Ca(alpha,gamma)44Ti reaction. Acceptable models should reproduce data on the 40Ca(alpha,gamma)44Ti reaction cross section as well as the measured primary gamma-ray spectra. This is only achieved when a coherent normalization of the slope of the level density and radiative strength function is performed. Thus, the overall shape of the experimental primary gamma-ray spectra puts a constraint on the input models for the rate calculations.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Structure of neutron stars with unified equations of state

We present a set of three unified equations of states (EoSs) based on the nuclear energy-density functional (EDF) theory.These EoSs are based on generalized Skyrme forces fitted to essentially all experimental atomic mass data and constrained to reproduce various properties of infinite nuclear matter as obtained from many-body calculations using realistic two- and three-body interactions. The structure of cold isolated neutron stars is discussed in connection with some astrophysical observations.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Superfluidity and entrainment in neutron-star crusts

Despite the absence of viscous drag, the neutron superfluid permeating the inner crust of a neutron star can still be strongly coupled to nuclei due to non-dissipative entrainment effects. Neutron superfluidity and entrainment have been systematically studied in all regions of the inner crust of a cold non-accreting neutron star in the framework of the band theory of solids. It is shown that in the intermediate layers of the inner crust a large fraction of "free" neutrons are actually entrained by the crust. The results suggest that a revision of the interpretation of many observable astrophysical phenomena might be necessary.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Masses of neutron stars and nuclei

We calculate the maximum mass of neutron stars for three different equations of state (EOS) based on generalized Skyrme functionals that are simultaneously fitted to essentially all the 2003 nuclear mass data (the rms deviation is 0.58 MeV in all three cases) and to one or other of three different equations of state of pure neutron matter, each determined by a different many-body calculation using realistic two- and three-body interactions but leading to significantly different degrees of stiffness at the high densities prevailing in neutron-star interiors. The observation of a neutron star with mass 1.97 $\pm$ 0.04 $\mathcal{M}_{\odot}$ eliminates the softest of our models (BSk19), but does not discriminate between BSk20 and BSk21. However, nuclear-mass measurements that have been made since our models were constructed strongly favor BSk21, our stiffest functional.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Self-consistent calculations of the strength function and radiative neutron capture cross section for stable and unstable tin isotopes

The E1 strength function for 15 stable and unstable Sn even-even isotopes from A=100 till A=176 are calculated using the self-consistent microscopic theory which, in addition to the standard (Q)RPA approach, takes into account the single-particle continuum and the phonon coupling. Our analysis shows two distinct regions for which the integral characteristics of both the giant and pygmy resonances behave rather differently. For neutron-rich nuclei, starting from $^{132}$Sn, we obtain a giant E1 resonance which significantly deviates from the widely-used systematics extrapolated from experimental data in the $β$-stability valley. We show that the inclusion of the phonon coupling is necessary for a proper description of the low-energy pygmy resonances and the corresponding transition densities for $A<$132 nuclei, while in the $A>132$ region the influence of phonon coupling is significantly smaller. The radiative neutron capture cross sections leading to the stable $^{124}$Sn and unstable $^{132}$Sn and $^{150}$Sn nuclei are calculated with both the (Q)RPA and the beyond-(Q)RPA strength functions and shown to be sensitive to both the predicted low-lying strength and the phonon coupling contribution. The comparison with the widely-used phenomenological Generalized Lorentzian approach shows considerable differences both for the strength function and the radiative neutron capture cross section. In particular, for the neutron-rich $^{150}$Sn, the reaction cross section is found to be increased by a factor greater than 20. We conclude that the present approach may provide a complete and coherent description of the $γ$-ray strength function for astrophysics applications. In particular, such calculations are highly recommended for a reliable estimate of the electromagnetic properties of exotic nuclei.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

The decompression of the outer neutron star crust and r-process nucleosynthesis

The rapid neutron-capture process, or r-process, is known to be fundamental for explaining the origin of approximately half of the A>60 stable nuclei observed in nature. In recent years nuclear astrophysicists have developed more and more sophisticated r-process models, by adding new astrophysical or nuclear physics ingredients to explain the solar system composition in a satisfactory way. Despite these efforts, the astrophysical site of the r-process remains unidentified. The composition of the neutron star outer crust material is investigated after the decompression that follows its possible ejection. The composition of the outer crust of a neutron star is estimated before and after decompression. Two different possible initial conditions are considered, namely an idealized crust composed of cold catalyzed matter and a crust initially in nuclear statistical equilibrium at temperatures around 10 GK. We show that in this second case before decompression and at temperatures typically corresponding to 8 GK, the Coulomb effect due to the high densities in the crust leads to an overall composition of the outer crust in neutron-rich nuclei with a mass distribution close to the solar system r-abundance distribution. Such distributions differ, however, from the solar one due to a systematic shift in the second peak to lower values. After decompression, the capture of the few neutrons per seed nucleus available in the hot outer crust leads to a final distribution of stable neutron-rich nuclei with a mass distribution of 80 < A < 140 nuclei in excellent agreement with the solar distribution, provided the outer crust is initially at temperatures around 8 GK and all layers of the outer crust are ejected. The decompression of the neutron star matter from the outer crust provides suitable conditions for a robust r-processing of the light species, i.e., r-nuclei with A < 140.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Breathing-mode measurements in Sn isotopes and isospin dependence of nuclear incompressibility

T. Li {\it et al.}[Phys. Rev. C {\bf 81}, 034309 (2010)] have analyzed their measured breathing-mode energies of some tin isotopes in terms of a first-order leptodermous expansion, and find for the symmetry-incompressibility coefficient $K_τ$ the value of -550 $\pm$ 100 MeV. Removing an approximation that they made, we find that the first-order estimate of $K_τ$ shifts to -661 $\pm$ 144 MeV. However, taking into account higher-order terms in the leptodermous expansion shows that the data are compatible with the significantly lower magnitudes indicated by both another experiment and some theoretical estimates.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Further explorations of Skyrme-Hartree-Fock-Bogoliubov mass formulas. XII: Stiffness and stability of neutron-star matter

We construct three new Hartree-Fock-Bogoliubov (HFB) mass models, labeled HFB-19, HFB-20, and HFB-21, with unconventional Skyrme forces containing $t_4$ and $t_5$ terms, i.e., density-dependent generalizations of the usual $t_1$ and $t_2$ terms, respectively. The new forces underlying these models are fitted respectively to three different realistic equations of state of neutron matter for which the density dependence of the symmetry energy ranges from the very soft to the very stiff, reflecting thereby our present lack of complete knowledge of the high-density behavior of nuclear matter. All unphysical instabilities of nuclear matter, including the transition to a polarized state in neutron-star matter, are eliminated with the new forces. At the same time the new models fit essentially all the available mass data with rms deviations of 0.58 MeV and give the same high quality fits to measured charge radii that we obtained in earlier models with conventional Skyrme forces. Being constrained by neutron matter, these new mass models, which all give similar extrapolations out to the neutron drip line, are highly appropriate for studies of the $r$-process and the outer crust of neutron stars. Moreover, the underlying forces, labeled BSk19, BSk20 and BSk21, respectively, are well adapted to the study of the inner crust and core of neutron stars. The new family of Skyrme forces thus opens the way to a unified description of all regions of neutron stars.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Spin and spin-isospin instabilities in asymmetric nuclear matter at zero and finite temperatures using Skyrme functionals

Self-consistent mean field methods based on phenomenological Skyrme effective interactions are known to exhibit spurious spin and spin-isospin instabilities both at zero and finite temperatures when applied to homogeneous nuclear matter at the densities encountered in neutron stars and in supernova cores. The origin of these instabilities is revisited in the framework of the nuclear energy density functional theory and a simple prescription is proposed to remove them. The stability of several Skyrme parametrizations is reexamined.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Unified description of neutron superfluidity in the neutron-star crust with analogy to anisotropic multi-band BCS superconductors

The neutron superfluidity in the inner crust of a neutron star has been traditionally studied considering either homogeneous neutron matter or only a small number of nucleons confined inside the spherical Wigner-Seitz cell. Drawing analogies with the recently discovered multi-band superconductors, we have solved the anisotropic multi-band BCS gap equations with Bloch boundary conditions, thus providing a unified description taking consistently into account both the free neutrons and the nuclear clusters. Calculations have been carried out using the effective interaction underlying our recent Hartree-Fock-Bogoliubov nuclear mass model HFB-16. We have found that even though the presence of inhomogeneities lowers the neutron pairing gaps, the reduction is much less than that predicted by previous calculations using the Wigner-Seitz approximation. We have studied the disappearance of superfluidity with increasing temperature. As an application we have calculated the neutron specific heat, which is an important ingredient for modeling the thermal evolution of newly-born neutron stars. This work provides a new scheme for realistic calculations of superfluidity in neutron-star crusts.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

A holistic approach to carbon-enhanced metal-poor stars

By considering the various CEMP subclasses separately, we try to derive, from the specific signatures imprinted on the abundances, parameters (such as metallicity, mass, temperature, and neutron source) characterizing AGB nucleosynthesis from the specific signatures imprinted on the abundances, and separate them from the impact of thermohaline mixing, first dredge-up, and dilution associated with the mass transfer from the companion.To put CEMP stars in a broad context, we collect abundances for about 180 stars of various metallicities, luminosity classes, and abundance patterns, from our own sample and from literature. First, we show that there are CEMP stars which share the properties of CEMP-s stars and CEMP-no stars (which we call CEMP-low-s stars). We also show that there is a strong correlation between Ba and C abundances in the s-only CEMP stars. This strongly points at the operation of the 13C neutron source in low-mass AGB stars. For the CEMP-rs stars (seemingly enriched with elements from both the s- and r-processes), the correlation of the N abundances with abundances of heavy elements from the 2nd and 3rd s-process peaks bears instead the signature of the 22Ne neutron source. Adding the fact that CEMP-rs stars exhibit O and Mg enhancements, we conclude that extremely hot conditions prevailed during the thermal pulses of the contaminating AGB stars. Finally, we argue that most CEMP-no stars (with no overabundances for the neutron-capture elements) are likely the extremely metal-poor counterparts of CEMP neutron-capture-rich stars. We also show that the C enhancement in CEMP-no stars declines with metallicity at extremely low metallicity ([Fe/H]~< -3.2). This trend is not predicted by any of the current AGB models.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Recent breakthroughs in Skyrme-Hartree-Fock-Bogoliubov mass formulas

We review our recent achievements in the construction of microscopic mass tables based on the Hartree-Fock-Bogoliubov method with Skyrme effective interactions. In the latest of our series of HFB-mass models, we have obtained our best fit ever to essentially all the available mass data, by treating the pairing more realistically than in any of our earlier models. The rms deviation on the 2149 measured masses of nuclei with N and Z>8 has been reduced for the first time in a mean field approach to 0.581 MeV. With the additional constraint on the neutron-matter equation of state, this new force is thus very well-suited for the study of neutron-rich nuclei and for the description of astrophysical environments like supernova cores and neutron-star crusts.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Skyrme-Hartree-Fock-Bogoliubov nuclear mass formulas: Crossing the 0.6 MeV threshold with microscopically deduced pairing

We present a new Skyrme-Hartree-Fock-Bogoliubov nuclear-mass model in which the contact-pairing force is constructed from microscopic pairing gaps of symmetric nuclear matter and neutron matter calculated from realistic two- and three-body forces, with medium-polarization effects included. With the pairing being treated more realistically than in any of our earlier models, the rms deviation with respect to essentially all the available mass data falls to 0.581 MeV, the best value ever found within the mean-field framework. Since our Skyrme force is also constrained by the properties of pure neutron matter, this new model is particularly well suited for application to astrophysical problems involving a neutron-rich environment, such as the elucidation of the r process of nucleosynthesis, and the description of supernova cores and neutron-star crusts.

0 citations0 reviewsNo code linkedOpen access
preprint2008arXiv

Improved predictions of nuclear reaction rates with the TALYS reaction code for astrophysical applications

Nuclear reaction rates of astrophysical applications are traditionally determined on the basis of Hauser-Feshbach reaction codes. These codes adopt a number of approximations that have never been tested, such as a simplified width fluctuation correction, the neglect of delayed or multiple-particle emission during the electromagnetic decay cascade, or the absence of the pre-equilibrium contribution at increasing incident energies. The reaction code TALYS has been recently updated to estimate the Maxwellian-averaged reaction rates that are of astrophysical relevance. These new developments enable the reaction rates to be calculated with increased accuracy and reliability and the approximations of previous codes to be investigated. The TALYS predictions for the thermonuclear rates of relevance to astrophysics are detailed and compared with those derived by widely-used codes for the same nuclear ingredients. It is shown that TALYS predictions may differ significantly from those of previous codes, in particular for nuclei for which no or little nuclear data is available. The pre-equilibrium process is shown to influence the astrophysics rates of exotic neutron-rich nuclei significantly. For the first time, the Maxwellian-averaged (n,2n) reaction rate is calculated for all nuclei and its competition with the radiative capture rate is discussed. The TALYS code provides a new tool to estimate all nuclear reaction rates of relevance to astrophysics with improved accuracy and reliability.

0 citations0 reviewsNo code linkedOpen access
preprint2005arXiv

BRUSLIB and NETGEN: the Brussels nuclear reaction rate library and nuclear network generator for astrophysics

Nuclear reaction rates are quantities of fundamental importance in astrophysics. Substantial efforts have been devoted in the last decades to measure or calculate them. The present paper presents for the first time a detailed description of the Brussels nuclear reaction rate library BRUSLIB and of the nuclear network generator NETGEN so as to make these nuclear data packages easily accessible to astrophysicists for a large variety of applications. BRUSLIB is made of two parts. The first one contains the 1999 NACRE compilation based on experimental data for 86 reactions with (mainly) stable targets up to Si. The second part of BRUSLIB concerns nuclear reaction rate predictions calculated within a statistical Hauser-Feshbach approximation, which limits the reliability of the rates to reactions producing compound nuclei with a high enough level density. These calculations make use of global and coherent microscopic nuclear models for the quantities entering the rate calculations. The use of such models is utterly important, and makes the BRUSLIB rate library unique. A description of the Nuclear Network Generator NETGEN that complements the BRUSLIB package is also presented. NETGEN is a tool to generate nuclear reaction rates for temperature grids specified by the user. The information it provides can be used for a large variety of applications, including Big Bang nucleosynthesis, the energy generation and nucleosynthesis associated with the non-explosive and explosive hydrogen to silicon burning stages, or the synthesis of the heavy nuclides through the s-, alpha- and r-, rp- or p-processes.

0 citations0 reviewsNo code linkedOpen access
preprint2005arXiv

UHE nuclei propagation and the interpretation of the ankle in the cosmic-ray spectrum

We consider the stochastic propagation of high-energy protons and nuclei in the cosmological microwave and infrared backgrounds, using revised photonuclear cross-sections and following primary and secondary nuclei in the full 2D nuclear chart. We confirm earlier results showing that the high-energy data can be fit with a pure proton extragalactic cosmic ray (EGCR) component if the source spectrum is \propto E^{-2.6}. In this case the ankle in the CR spectrum may be interpreted as a pair-production dip associated with the propagation. We show that when heavier nuclei are included in the source with a composition similar to that of Galactic cosmic-rays (GCRs), the pair-production dip is not present unless the proton fraction is higher than 85%. In the mixed composition case, the ankle recovers the past interpretation as the transition from GCRs to EGCRs and the highest energy data can be explained by a harder source spectrum \propto E^{-2.2} - E^{-2.3}, reminiscent of relativistic shock acceleration predictions, and in good agreement with the GCR data at low-energy and holistic scenarios.

0 citations0 reviewsNo code linkedOpen access
preprint2003arXiv

More lead stars

The standard model for the operation of the s-process in asymptotic giant branch (AGB) stars predicts that low-metallicity ([Fe/H] <= -1) AGB stars should exhibit large overabundances of Pb and Bi as compared to other s-elements. Eight CH stars (with [Fe/H] ranging from -1.5 to -2.5) are studied in the present paper using high-resolution spectra (R = lambda/Delta lambda = 135 000), and five of them appear to be enriched in lead (with [Pb/Ce] ~ 0.7). Abundances of heavy elements in these `lead stars' (defined as stars enriched in s-elements and with [Pb/hs] >= 1, hs being any of Ba, La or Ce) are consistent with the predictions for the s-process operating in low-metallicity AGB stars as a consequence of the `partial mixing' of protons below the convective hydrogen envelope. Two among the studied stars (HD 189711 and V Ari) add to a growing number of low-metallicity stars which do not conform to these predictions, however. Variations on the canonical proton-mixing scenario for the operation of the s-process in low-metallicity stars, that could account for these discrepant stars, are briefly discussed.

0 citations0 reviewsNo code linkedOpen access
preprint2001arXiv

The p- and r-processes: reviews and other views

A review is presented of the p-process in Type II supernovae, one of its goals being to enlighten the changes in views on this nucleosynthesis mechanism since the work of Jean and Jim on the subject in 1975. Specific discussions are also devoted to cases of particular interest, like the light Mo and Ru stable isotopes, the rare nuclide 138La or the radionuclide 146Sm. Some comments of diverse natures are also made on the r-process. These considerations do not aim at really providing an exhaustive review of the many nuclear physics and astrophysics intricacies of this process. In contrast, they are hoped to complement or to put in perspective other views that are often expressed in relation with this nucleosynthesis mechanism

0 citations0 reviewsNo code linkedOpen access