-
Constraints on the Neutron-Star Structure from the Clocked X-Ray Burster 1RXS J180408.9$-$342058
Authors:
Akira Dohi,
Wataru Iwakiri,
Nobuya Nishimura,
Tsuneo Noda,
Shigehiro Nagataki,
Masa-aki Hashimoto
Abstract:
Type-I X-ray bursts are rapid-brightening transient phenomena on the surfaces of accreting neutron stars (NSs). Some X-ray bursts, called {\it clocked bursters}, exhibit regular behavior with similar light curve profiles in their burst sequences. The periodic nature of clocked bursters has the advantage of constraining X-ray binary parameters and physics inside the NS. In the present study, we com…
▽ More
Type-I X-ray bursts are rapid-brightening transient phenomena on the surfaces of accreting neutron stars (NSs). Some X-ray bursts, called {\it clocked bursters}, exhibit regular behavior with similar light curve profiles in their burst sequences. The periodic nature of clocked bursters has the advantage of constraining X-ray binary parameters and physics inside the NS. In the present study, we compute numerical models, based on different equations of state and NS masses, which are compared with the observation of a recently identified clocked burster, 1RXS J180408.9$-$342058. We find that the relation between accretion rate and recurrence time is highly sensitive to the NS mass and radius. We determine, in particular, that 1RXS J180408.9$-$342058 appears to possess a mass less than $1.7M_{\odot}$ and favors a stiffer nuclear equation of state (with an NS radius $\gtrsim12.7{\rm km}$). Consequently, the observations of this new clocked burster may provide additional constraints for probing the structure of NSs.
△ Less
Submitted 5 November, 2023;
originally announced November 2023.
-
Postfission properties of uranium isotopes: A hybrid method with Langevin dynamics and the Hauser-Feshbach statistical model
Authors:
Shyoya Tanaka,
Nobuya Nishimura,
Futoshi Minato,
Yoshihiro Aritomo
Abstract:
Background: Precise understanding of nuclear fission is crucial for experimental and theoretical nuclear physics, astrophysics, and industrial applications; however, the complete physical mechanics is unresolved due to the complexities.
Purpose: In this study, we present a new method to describe the dynamical-fission process and following prompt-neutron emission, where we combine the dynamical f…
▽ More
Background: Precise understanding of nuclear fission is crucial for experimental and theoretical nuclear physics, astrophysics, and industrial applications; however, the complete physical mechanics is unresolved due to the complexities.
Purpose: In this study, we present a new method to describe the dynamical-fission process and following prompt-neutron emission, where we combine the dynamical fission calculation based on the Langevin method and the Hauser-Feshbach statistical model.
Methods: Two methods are connected smoothly within the universal charge distribution and the energy conservation, allowing us to calculate a sequence of fission dynamics and post-fission phase, including prompt neutron emission.
Results: Using a certain set of model parameters, we successfully reproduce the experimental primary-fission yields, total kinetic energy, independent-fission yields, and prompt neutron emissions for the neutron induced fission of ${}^{236}$U, a compound nucleus of ${\rm n} + {}^{235}{\rm U}$. We elucidate the physical mechanism of the characteristic features observed in previous experiments, such as shell properties. Additionally, we apply our calculation to two very neutron-rich uranium isotopes, i.e., ${}^{250}$U and ${}^{255}$U, which are not experimentally confirmed but are important for r-process nucleosynthesis. Theoretical results indicate that ${}^{250}$U exhibits an asymmetric multiple-peak fission yield distribution, while the neutron-rich ${}^{255}$U has a single peak due to symmetric fission. Our method predicts post-neutron emission fragments, where ${}^{250}$U shows a stronger neutron emissivity than ${}^{255}$U.
Conclusions: Our framework is highly reproducible in the experiments and shows that the number of emitted neutrons after fission differs significantly in neutron-rich uranium fission depending on distributions of fission variables.
△ Less
Submitted 23 November, 2023; v1 submitted 18 July, 2023;
originally announced July 2023.
-
The Impacts of Neutron-Star Structure and Base Heating on Type I X-Ray Bursts and Code Comparison
Authors:
Guoqing Zhen,
Guoliang Lv,
Helei Liu,
Akira Dohi,
Nobuya Nishimura,
Chunhua Zhu,
Liyu Song,
Weiyang Wang,
Renxin Xu
Abstract:
Type I X-ray bursts are rapidly brightening phenomena triggered by thermonuclear burning on accreting layer of a neutron star (NS). The light curves represent the physical properties of NSs and the nuclear reactions on the proton-rich nuclei. The numerical treatments of the accreting NS and physics of the NS interior are not established, which shows uncertainty in modelling for observed X-ray ligh…
▽ More
Type I X-ray bursts are rapidly brightening phenomena triggered by thermonuclear burning on accreting layer of a neutron star (NS). The light curves represent the physical properties of NSs and the nuclear reactions on the proton-rich nuclei. The numerical treatments of the accreting NS and physics of the NS interior are not established, which shows uncertainty in modelling for observed X-ray light curves. In this study, we investigate theoretical X-ray-burst models, compared with burst light curves with GS~1826-24 observations. We focus on the impacts of the NS mass, the NS radius, and base-heating on the NS surface using the MESA code. We find a monotonic correlation between the NS mass and the parameters of the light curve. The higher the mass, the longer the recurrence time and the greater the peak luminosity. While the larger the radius, the longer the recurrence time, the peak luminosity remains nearly constant. In the case of increasing base heating, both the recurrence time and peak luminosity decrease. We also examine the above results using with a different numerical code, HERES, based on general relativity and consider the central NS. We find that the burst rate, burst energy and burst strength are almost same in two X-ray burst codes by adjusting the base-heat parameter in MESA (the relative errors $\lesssim5\%$), while the duration time and the rise time are significantly different between (the relative error is possibly $\sim50\%$). The peak luminosity and the e-folding time are ragged between two codes for different accretion rates.
△ Less
Submitted 14 April, 2023;
originally announced April 2023.
-
Cross Sections of the $^{83}$Rb(p,$γ)^{84}$Sr and $^{84}$Kr(p,$γ)^{85}$Rb Reactions at Energies Characteristic of the Astrophysical $γ$ Process
Authors:
M. Williams,
B. Davids,
G. Lotay,
N. Nishimura,
T. Rauscher,
S. A. Gillespie,
M. Alcorta,
A. M. Amthor,
G. C. Ball,
S. S. Bhattacharjee,
V. Bildstein,
W. N. Catford,
D. T. Doherty,
N. E. Esker,
A. B. Garnsworthy,
G. Hackman,
K. Hudson,
A. Lennarz,
C. Natzke,
B. Olaizola,
A. Psaltis,
C. E. Svensson,
J. Williams,
D. Walter,
D. Yates
Abstract:
We have measured the cross section of the $^{83}$Rb(p,$γ)^{84}$Sr radiative capture reaction in inverse kinematics using a radioactive beam of $^{83}$Rb at incident energies of 2.4 and $2.7 A$ MeV. Prior to the radioactive beam measurement, the $^{84}$Kr(p,$γ)^{85}$Rb radiative capture reaction was measured in inverse kinematics using a stable beam of $^{84}$Kr at an incident energy of $2.7 A$ MeV…
▽ More
We have measured the cross section of the $^{83}$Rb(p,$γ)^{84}$Sr radiative capture reaction in inverse kinematics using a radioactive beam of $^{83}$Rb at incident energies of 2.4 and $2.7 A$ MeV. Prior to the radioactive beam measurement, the $^{84}$Kr(p,$γ)^{85}$Rb radiative capture reaction was measured in inverse kinematics using a stable beam of $^{84}$Kr at an incident energy of $2.7 A$ MeV. The effective relative kinetic energies of these measurements lie within the relevant energy window for the $γ$ process in supernovae. The central values of the measured partial cross sections of both reactions were found to be $0.17-0.42$ times the predictions of statistical model calculations. Assuming the predicted cross section at other energies is reduced by the same factor leads to a slightly higher calculated abundance of the $p$ nucleus $^{84}$Sr, caused by the reduced rate of the $^{84}$Sr($γ$,p)$^{83}$Rb reaction derived from the present measurement.
△ Less
Submitted 17 March, 2023;
originally announced March 2023.
-
Impacts of the direct URCA and Superfluidity inside a Neutron Star on Type-I X-Ray Bursts and X-Ray Superbursts
Authors:
Akira Dohi,
Nobuya Nishimura,
Hajime Sotani,
Tsuneo Noda,
Helei Liu,
Shigehiro Nagataki,
Masa-aki Hashimoto
Abstract:
We investigate the impacts of neutrino cooling mechanism inside the neutron star (NS) core on the light curves of type-I X-ray bursts and X-ray superbursts. From several observations of NS thermal evolution, physical processes of fast neutrino cooling, such as the direct Urca (DU) process, are indicated. They significantly decrease the surface temperature of NSs, though the cooling effect could be…
▽ More
We investigate the impacts of neutrino cooling mechanism inside the neutron star (NS) core on the light curves of type-I X-ray bursts and X-ray superbursts. From several observations of NS thermal evolution, physical processes of fast neutrino cooling, such as the direct Urca (DU) process, are indicated. They significantly decrease the surface temperature of NSs, though the cooling effect could be suppressed by nucleon superfluidity. In the present study, focusing on the DU process and nucleon superfluidity, we investigate the effects of NS cooling on the X-ray bursts using a general-relativistic stellar-evolution code. We find that the DU process leads find the longer recurrence time and the higher peak luminosity, which could be obstructed by the neutrons superfluidity. We also apply our burst models to the comparison with {\it Clocked burster} GS 1826$-$24, and to the recurrence time of superburst triggered by carbon ignition. These effects are significant within a certain range of binary parameters and uncertainty of the NS equation of state.
△ Less
Submitted 20 March, 2023; v1 submitted 31 August, 2022;
originally announced August 2022.
-
Horizons: Nuclear Astrophysics in the 2020s and Beyond
Authors:
H. Schatz,
A. D. Becerril Reyes,
A. Best,
E. F. Brown,
K. Chatziioannou,
K. A. Chipps,
C. M. Deibel,
R. Ezzeddine,
D. K. Galloway,
C. J. Hansen,
F. Herwig,
A. P. Ji,
M. Lugaro,
Z. Meisel,
D. Norman,
J. S. Read,
L. F. Roberts,
A. Spyrou,
I. Tews,
F. X. Timmes,
C. Travaglio,
N. Vassh,
C. Abia,
P. Adsley,
S. Agarwal
, et al. (140 additional authors not shown)
Abstract:
Nuclear Astrophysics is a field at the intersection of nuclear physics and astrophysics, which seeks to understand the nuclear engines of astronomical objects and the origin of the chemical elements. This white paper summarizes progress and status of the field, the new open questions that have emerged, and the tremendous scientific opportunities that have opened up with major advances in capabilit…
▽ More
Nuclear Astrophysics is a field at the intersection of nuclear physics and astrophysics, which seeks to understand the nuclear engines of astronomical objects and the origin of the chemical elements. This white paper summarizes progress and status of the field, the new open questions that have emerged, and the tremendous scientific opportunities that have opened up with major advances in capabilities across an ever growing number of disciplines and subfields that need to be integrated. We take a holistic view of the field discussing the unique challenges and opportunities in nuclear astrophysics in regards to science, diversity, education, and the interdisciplinarity and breadth of the field. Clearly nuclear astrophysics is a dynamic field with a bright future that is entering a new era of discovery opportunities.
△ Less
Submitted 16 May, 2022;
originally announced May 2022.
-
New constraints on the neutron-star mass and radius relation from terrestrial nuclear experiments
Authors:
Hajime Sotani,
Nobuya Nishimura,
Tomoya Naito
Abstract:
The determination of the equation of state (EOS) for nuclear matter has been one of the biggest problems in nuclear astrophysics, because the EOS is essential for determining the properties of neutron stars. To constrain the density-dependence of the nuclear symmetry energy, several nuclear experiments, e.g., reported by the S$π$RIT and PREX-II collaborations, have recently been performed. However…
▽ More
The determination of the equation of state (EOS) for nuclear matter has been one of the biggest problems in nuclear astrophysics, because the EOS is essential for determining the properties of neutron stars. To constrain the density-dependence of the nuclear symmetry energy, several nuclear experiments, e.g., reported by the S$π$RIT and PREX-II collaborations, have recently been performed. However, since their uncertainties are still large, additional constraints such as astronomical observations must be crucial. In addition, it is interesting to see the effect of their reported value on neutron star properties. In this study, focusing on the relatively lower density region, we investigate the allowed area of the neutron-star mass and radius relation by assuming the constraints from S$π$RIT and PREX-II. Each region predicted by these experiments is still consistent with the allowed area constrained by the various astronomical observations. Our results show that terrestrial nuclear experiments must provide further constraints on the EOS for neutron stars, complementing astronomical observations.
△ Less
Submitted 16 March, 2022; v1 submitted 20 February, 2022;
originally announced March 2022.
-
Effects of Nuclear Equation of State on Type-I X-ray Bursts: Interpretation of the X-ray Bursts from GS 1826-24
Authors:
A. Dohi,
N. Nishimura,
M. Hashimoto,
Y. Matsuo,
T. Noda,
S. Nagataki
Abstract:
Type I X-ray bursts are thermonuclear explosions on the neutron star (NS) surface caused by mass accretion from a companion star. Observations of X-ray bursts provide valuable information on X-ray binary systems, e.g., binary parameters, the chemical composition of accreted matter, and the nuclear equation of state (EOS) of NSs. There have been several theoretical studies to constrain the physics…
▽ More
Type I X-ray bursts are thermonuclear explosions on the neutron star (NS) surface caused by mass accretion from a companion star. Observations of X-ray bursts provide valuable information on X-ray binary systems, e.g., binary parameters, the chemical composition of accreted matter, and the nuclear equation of state (EOS) of NSs. There have been several theoretical studies to constrain the physics of X-ray bursters. However, they have mainly focused on the burning layers above the solid crust of the NS, which brings up issues of the treatment of NS gravitation and internal energy. In this study, focusing on the microphysics inside NSs, we calculate a series of X-ray bursts using a general-relativistic stellar-evolution code with several NS EOSs. We compare the X-ray-burst models with the burst parameters of a clocked burster associated with GS 1826-24. We find a monotonic correlation between the NS radius and the light-curve profile. A larger radius shows a higher recurrence time and a large peak luminosity. In contrast, the dependence of light curves on the NS mass becomes more complicated, where the neutrino cooling suppresses the efficiency of nuclear ignition. We also constrain the EOS and mass of GS 1826-24, i.e., stiffer EOSs, corresponding to larger NS radii, are not preferred due to too-high peak luminosity. The EOS and the cooling and heating of NSs are important to discuss the theoretical and observational properties of X-ray bursts.
△ Less
Submitted 5 January, 2022; v1 submitted 28 May, 2021;
originally announced May 2021.
-
Fission fragment distributions of neutron-rich nuclei based on Langevin calculations: toward r-process simulations
Authors:
Mizuki Okubayashi,
Shoya Tanaka,
Yoshihiro Aritomo,
Shoma Ishizaki,
Shota Amano,
Nobuya Nishimura
Abstract:
The nuclear fission of very neuron-rich nuclei related to the r-process is essential for the termination of nucleosynthesis flows on the nuclear chart and the final abundances. Nevertheless, most of the available fission data for the r-process calculations are based on theory predictions, including phenomenological treatments. In this study, we calculated a series of nuclear fission distribution f…
▽ More
The nuclear fission of very neuron-rich nuclei related to the r-process is essential for the termination of nucleosynthesis flows on the nuclear chart and the final abundances. Nevertheless, most of the available fission data for the r-process calculations are based on theory predictions, including phenomenological treatments. In this study, we calculated a series of nuclear fission distribution for neutron-rich nuclei away from the beta-stability line. As most of these nuclei are experimentally unknown, we are based on theoretical calculations based on the dynamical fission model with the Langevin method. We performed fission distribution calculations for neutron-rich actinoid nuclei, applicable to the r-process nucleosynthesis simulations. In the present paper, we compared the obtained mass and charge distributions with experimental data. We also show the results of the systematic behaviour of mass distribution for neutron-rich U and Fm isotopes.
△ Less
Submitted 19 May, 2021;
originally announced May 2021.
-
Statistical Analysis of the Relation between Coronal Mass Ejections and Solar Energetic Particles
Authors:
Kosuke Kihara,
Yuwei Huang,
Nobuhiko Nishimura,
Nariaki V. Nitta,
Seiji Yashiro,
Kiyoshi Ichimoto,
Ayumi Asai
Abstract:
To improve the forecasting capability of impactful solar energetic particle (SEP) events, the relation between coronal mass ejections (CMEs) and SEP events needs to be better understood. Here we present a statistical study of SEP occurrences and timescales with respect to the CME source locations and speeds, considering all 257 fast ($v_{CME}$ $\ge$ 900 km/s) and wide (angular width $\ge$ 60…
▽ More
To improve the forecasting capability of impactful solar energetic particle (SEP) events, the relation between coronal mass ejections (CMEs) and SEP events needs to be better understood. Here we present a statistical study of SEP occurrences and timescales with respect to the CME source locations and speeds, considering all 257 fast ($v_{CME}$ $\ge$ 900 km/s) and wide (angular width $\ge$ 60$^{\circ}$) CMEs that occurred between December 2006 and October 2017. We associate them with SEP events at energies above 10 MeV. Examination of the source region of each CME reveals that CMEs more often accompany a SEP event if they originate from the longitude of E20-W100 relative to the observer. However, a SEP event could still be absent if the CME is $<$ 2000 km/s. For the associated CME-SEP pairs, we compute three timescales for each of the SEP events, following Kahler (2005, 2013); namely the timescale of the onset (TO), the rise time (TR), and the duration (TD). They are correlated with the longitude of the CME source region relative to the footpoint of the Parker spiral ($ΔΦ$) and $v_{CME}$. The TO tends to be short for $|ΔΦ|$ $<$ 60$^{\circ}$ . This trend is weaker for TR and TD. The SEP timescales are only weakly correlated with $v_{CME}$. Positive correlations of both TR and TD with $v_{CME}$ are seen in poorly connected (large $|ΔΦ|$) events. Additionally, TO appears to be negatively correlated with $v_{CME}$ for events with small $|ΔΦ|$.
△ Less
Submitted 15 July, 2020;
originally announced July 2020.
-
r-Process enrichment in the Galactic halo characterized by nucleosynthesis variation in the ejecta of coalescing neutron star binaries
Authors:
Takuji Tsujimoto,
Nobuya Nishimura,
Koutarou Kyutoku
Abstract:
A large star-to-star variation in the abundances of r-process elements, as seen in the [Eu/Fe] ratio for Galactic halo stars, is a prominent feature that is distinguishable from other heavy elements. It is, in part, caused by the presence of highly r-process enriched stars, classified as r-II stars ([Eu/Fe]>+1). In parallel, halo stars show that the ratio of a light r-process element (Y) to Eu is…
▽ More
A large star-to-star variation in the abundances of r-process elements, as seen in the [Eu/Fe] ratio for Galactic halo stars, is a prominent feature that is distinguishable from other heavy elements. It is, in part, caused by the presence of highly r-process enriched stars, classified as r-II stars ([Eu/Fe]>+1). In parallel, halo stars show that the ratio of a light r-process element (Y) to Eu is tightly correlated with [Eu/Fe], giving the lowest [Y/Eu] ratio that levels off at r-II stars. On the other hand, recent hydrodynamical simulations of coalescing double neutron stars (cNSNSs) have suggested that r-process sites may be separated into two classes providing different electron-fraction distributions: tidally-driven dynamical ejecta and (dynamical or postmerger) non-tidal ejecta. Here, we show that a widely spanning feature of [Eu/Fe] can be reproduced by models that consider the different masses of tidally-driven dynamical ejecta from both cNSNSs and coalescing black hole/neutron star binaries (cBHNSs). In addition, the observed [Y/Eu] trend is explained by the combined nucleosynthesis in two kinds of ejecta with varying mass asymmetry in double NS systems. Our scenario suggests that massive tidally-driven dynamical ejecta accompanied by massive non-tidal part from cNSNSs or cBHNSs could alone accommodate r-II abundances, including an actinide boost in some cases. The event rate for cNSNSs estimated from our study agrees with the latest result of ~1000 (90% confidence interval of 110-3840) Gpc$^{-3}$yr$^{-1}$ by gravitational-wave detection, and a few events per Gpc$^3$ per year of cBHNSs associated with r-process production are predicted to emerge.
△ Less
Submitted 23 December, 2019;
originally announced December 2019.
-
Impact of Uncertainties in Astrophysical Reaction Rates on Nucleosynthesis in the $νp$ Process
Authors:
T. Rauscher,
N. Nishimura,
G. Cescutti,
R. Hirschi,
A. St. J. Murphy,
C. Fröhlich
Abstract:
The $νp$ process appears in proton-rich, hot matter which is expanding in a neutrino wind and may be realised in explosive environments such as core-collapse supernovae or in outflows from accretion disks. The impact of uncertainties in nuclear reaction cross sections on the finally produced abundances has been studied by applying Monte Carlo variation of all astrophysical reaction rates in a larg…
▽ More
The $νp$ process appears in proton-rich, hot matter which is expanding in a neutrino wind and may be realised in explosive environments such as core-collapse supernovae or in outflows from accretion disks. The impact of uncertainties in nuclear reaction cross sections on the finally produced abundances has been studied by applying Monte Carlo variation of all astrophysical reaction rates in a large reaction network. As the detailed astrophysical conditions of the $νp$ process still are unknown, a parameter study was performed, with 23 trajectories covering a large range of entropies and $Y_\mathrm{e}$. The resulting abundance uncertainties are given for each trajectory. The $νp$ process has been speculated to contribute to the light $p$ nuclides but it was not possible so far to reproduce the solar isotope ratios. It is found that it is possible to reproduce the solar $^{92}$Mo/$^{94}$Mo abundance ratio within nuclear uncertainties, even within a single trajectory. The solar values of the abundances in the Kr-Sr region relative to the Mo region, however, cannot be achieved within a single trajectory. They may still be obtained from a weighted superposition of different trajectories, though, depending on the actual conditions in the production site. For a stronger constraint of the required conditions, it would be necessary to reduce the uncertainties in the 3$α$ and $^{56}$Ni(n,p)$^{56}$Co rates at temperatures $T>3$ GK.
△ Less
Submitted 12 September, 2019; v1 submitted 7 September, 2019;
originally announced September 2019.
-
Uncertainties in $ν$p-process nucleosynthesis from Monte Carlo variation of reaction rates
Authors:
N. Nishimura,
T. Rauscher,
R. Hirschi,
G. Cescutti,
A. St. J. Murphy,
C. Fröhlich
Abstract:
It has been suggested that a $ν$p process can occur when hot, dense, and proton-rich matter is expanding within a strong flux of anti-neutrinos. In such an environment, proton-rich nuclides can be produced in sequences of proton captures and (n,p) reactions, where the free neutrons are created in situ by $\overlineν_\mathrm{e}+\mathrm{p} \rightarrow \mathrm{n}+\mathrm{e}^+$ reactions. The detailed…
▽ More
It has been suggested that a $ν$p process can occur when hot, dense, and proton-rich matter is expanding within a strong flux of anti-neutrinos. In such an environment, proton-rich nuclides can be produced in sequences of proton captures and (n,p) reactions, where the free neutrons are created in situ by $\overlineν_\mathrm{e}+\mathrm{p} \rightarrow \mathrm{n}+\mathrm{e}^+$ reactions. The detailed hydrodynamic evolution determines where the nucleosynthesis path turns off from N = Z line and how far up the nuclear chart it runs. In this work, the uncertainties on the final isotopic abundances stemming from uncertainties in the nuclear reaction rates were investigated in a large-scale Monte Carlo approach, simultaneously varying ten thousand reactions. A large range of model conditions was investigated because a definitive astrophysical site for the $ν$p process has not yet been identified. The present parameter study provides, for each model, identification of the key nuclear reactions dominating the uncertainty for a given nuclide abundance. As all rates appearing in the $ν$p process involve unstable nuclei, and thus only theoretical rates are available, the final abundance uncertainties are larger than those for nucleosynthesis processes closer to stability. Nevertheless, most uncertainties remain below a factor of three in trajectories with robust nucleosynthesis. More extreme conditions allow production of heavier nuclides but show larger uncertainties because of the accumulation of the uncertainties in many rates and because the termination of nucleosynthesis is not at equilibrium conditions. It is also found that the solar ratio of the abundances of ${}^{92}$Mo and ${}^{94}$Mo could be reproduced within uncertainties.
△ Less
Submitted 30 July, 2019;
originally announced July 2019.
-
Impact of Uncertainties in Nuclear Reaction Cross Sections on p-Nucleosynthesis in Thermonuclear Supernovae
Authors:
T. Rauscher,
N. Nishimura,
G. Cescutti,
R. Hirschi,
A. St. J. Murphy,
C. Travaglio
Abstract:
The propagation of uncertainties in reaction cross sections and rates of neutron-, proton-, and alpha-induced reactions into the final isotopic abundances obtained in nucleosynthesis models is an important issue in studies of nucleosynthesis and Galactic Chemical Evolution. We developed a Monte Carlo method to allow large-scale postprocessing studies of the impact of nuclear uncertainties on nucle…
▽ More
The propagation of uncertainties in reaction cross sections and rates of neutron-, proton-, and alpha-induced reactions into the final isotopic abundances obtained in nucleosynthesis models is an important issue in studies of nucleosynthesis and Galactic Chemical Evolution. We developed a Monte Carlo method to allow large-scale postprocessing studies of the impact of nuclear uncertainties on nucleosynthesis. Temperature-dependent rate uncertainties combining realistic experimental and theoretical uncertainties are used. From detailed statistical analyses uncertainties in the final abundances are derived as probability density distributions. Furthermore, based on rate and abundance correlations an automated procedure identifies the most important reactions in complex flow patterns from superposition of many zones or tracers. The method so far was already applied to a number of nucleosynthesis processes. Here we focus on the production of p-nuclei in white dwarfs exploding as thermonuclear (type Ia) supernovae. We find generally small uncertainties in the final abundances despite of the dominance of theoretical nuclear uncertainties. A separate analysis of low- and high-density regions indicates that the total uncertainties are dominated by the high-density regions.
△ Less
Submitted 16 August, 2019; v1 submitted 22 July, 2019;
originally announced July 2019.
-
The s-process nucleosynthesis in low mass stars: impact of the uncertainties in the nuclear physics determined by Monte Carlo variations
Authors:
G. Cescutti,
R. Hirschi,
N. Nishimura,
T. Rauscher,
J. den Hartogh,
A. St. J. Murphy,
S. Cristallo
Abstract:
We investigated the impact of uncertainties in neutron-capture and weak reactions (on heavy elements) on the s-process nucleosynthesis in low-mass stars using a Monte-Carlo based approach. We performed extensive nuclear reaction network calculations that include newly evaluated temperature-dependent upper and lower limits for the individual reaction rates. Our sophisticated approach is able to eva…
▽ More
We investigated the impact of uncertainties in neutron-capture and weak reactions (on heavy elements) on the s-process nucleosynthesis in low-mass stars using a Monte-Carlo based approach. We performed extensive nuclear reaction network calculations that include newly evaluated temperature-dependent upper and lower limits for the individual reaction rates. Our sophisticated approach is able to evaluate the reactions that impact more significantly the final abundances. We found that beta-decay rate uncertainties affect typically nuclides near s-process branchings, whereas most of the uncertainty in the final abundances is caused by uncertainties in neutron capture rates, either directly producing or destroying the nuclide of interest. Combined total nuclear uncertainties due to reactions on heavy elements are approximately 50%.
△ Less
Submitted 10 October, 2018;
originally announced October 2018.
-
Early chemical evolution of Zn driven by magnetorotational supernovae and the pathway to the solar Zn composition
Authors:
Takuji Tsujimoto,
Nobuya Nishimura
Abstract:
The site of Zn production remains an elusive and challenging problem in astrophysics. A large enhancement of the [Zn/Fe] ratios of very metal-poor stars in the Galactic halo suggests the death of short-lived massive stars, i.e., core-collapse supernovae (CCSNe), as one major site for Zn production. Previous studies have claimed that some specific CCSNe can produce Zn in sufficient quantities. Howe…
▽ More
The site of Zn production remains an elusive and challenging problem in astrophysics. A large enhancement of the [Zn/Fe] ratios of very metal-poor stars in the Galactic halo suggests the death of short-lived massive stars, i.e., core-collapse supernovae (CCSNe), as one major site for Zn production. Previous studies have claimed that some specific CCSNe can produce Zn in sufficient quantities. However, it remains unclear which models can withstand the critical test of observations. Using a Zn abundance feature similar to that of r-process elements in faint satellite galaxies, we find evidence that Zn production took place through much rarer events than canonical CCSNe. This finding can be unified with the implied decrease in the rate of Zn production with an increasing metallicity for Galactic halo stars, which narrows down the major site of Zn production in the early galaxy to magneto-rotational SNe (MR-SNe). On the other hand, in the later phase of galactic evolution, we predict that the major Zn-production site switched from MR-SNe to thermonuclear SNe (SNe Ia). According to this scenario, an accumulation of the contributions from two types of SNe eventually led to the solar isotope composition of Zn which mainly owes 66,68Zn to MR-SNe and 64Zn to SNe Ia triggered by He-detonation. The requirement of Zn production in SNe Ia sheds a new light on the hot debate on the scenario for SN Ia progenitors, suggesting that a He-detonation model might be one major channel for SNe Ia.
△ Less
Submitted 7 August, 2018;
originally announced August 2018.
-
Uncertainties in the production of p nuclides in SN Ia determined by Monte Carlo variations
Authors:
T. Rauscher,
N. Nishimura,
R. Hirschi,
A. St. J. Murphy,
G. Cescutti,
C. Travaglio
Abstract:
Several thousand tracers from a 2D model of a thermonuclear supernova were used in a Monte Carlo post-processing approach to determine p-nuclide abundance uncertainties originating from nuclear physics uncertainties in the reaction rates.
Several thousand tracers from a 2D model of a thermonuclear supernova were used in a Monte Carlo post-processing approach to determine p-nuclide abundance uncertainties originating from nuclear physics uncertainties in the reaction rates.
△ Less
Submitted 28 July, 2018;
originally announced July 2018.
-
Uncertainties in s-process nucleosynthesis in low mass stars determined from Monte Carlo variations
Authors:
G. Cescutti,
R. Hirschi,
N. Nishimura,
J. W. den Hartogh,
T. Rauscher,
A. St. J. Murphy,
S. Cristallo
Abstract:
The main s-process taking place in low mass stars produces about half of the elements heavier than iron. It is therefore very important to determine the importance and impact of nuclear physics uncertainties on this process. We have performed extensive nuclear reaction network calculations using individual and temperature-dependent uncertainties for reactions involving elements heavier than iron,…
▽ More
The main s-process taking place in low mass stars produces about half of the elements heavier than iron. It is therefore very important to determine the importance and impact of nuclear physics uncertainties on this process. We have performed extensive nuclear reaction network calculations using individual and temperature-dependent uncertainties for reactions involving elements heavier than iron, within a Monte Carlo framework. Using this technique, we determined the uncertainty in the main s-process abundance predictions due to nuclear uncertainties link to weak interactions and neutron captures on elements heavier than iron. We also identified the key nuclear reactions dominating these uncertainties. We found that $β$-decay rate uncertainties affect only a few nuclides near s-process branchings, whereas most of the uncertainty in the final abundances is caused by uncertainties in neutron capture rates, either directly producing or destroying the nuclide of interest. Combined total nuclear uncertainties due to reactions on heavy elements are in general small (less than 50%). Three key reactions, nevertheless, stand out because they significantly affect the uncertainties of a large number of nuclides. These are $^{56}$Fe(n,$γ$), $^{64}$Ni(n,$γ$), and $^{138}$Ba(n,$γ$). We discuss the prospect of reducing uncertainties in the key reactions identified in this study with future experiments.
△ Less
Submitted 3 May, 2018;
originally announced May 2018.
-
Sensitivity to neutron captures and beta-decays of the enhanced s-process in rotating massive stars at low metallicities
Authors:
N. Nishimura,
R. Hirschi,
T. Rauscher
Abstract:
The s-process in massive stars, producing nuclei up to $A\approx 90$, has a different behaviour at low metallicity if stellar rotation is significant. This enhanced s-process is distinct from the s-process in massive stars around solar metallicity, and details of the nucleosynthesis are poorly known. We investigated nuclear physics uncertainties in the enhanced s-process in metal-poor stars within…
▽ More
The s-process in massive stars, producing nuclei up to $A\approx 90$, has a different behaviour at low metallicity if stellar rotation is significant. This enhanced s-process is distinct from the s-process in massive stars around solar metallicity, and details of the nucleosynthesis are poorly known. We investigated nuclear physics uncertainties in the enhanced s-process in metal-poor stars within a Monte-Carlo framework. We applied temperature-dependent uncertainties of reaction rates, distinguishing contributions from the ground state and from excited states. We found that the final abundance of several isotopes shows uncertainties larger than a factor of 2, mostly due to the neutron capture uncertainties. A few nuclei around branching points are affected by uncertainties in the $β$-decay.
△ Less
Submitted 15 February, 2018;
originally announced February 2018.
-
Impacts of nuclear-physics uncertainties in the s-process determined by Monte-Carlo variations
Authors:
N. Nishimura,
G. Cescutti,
R. Hirschi,
T. Rauscher,
J. den Hartogh,
A. St. J. Murphy
Abstract:
The s-process, a production mechanism based on slow-neutron capture during stellar evolution, is the origin of about half the elements heavier than iron. Abundance predictions for s-process nucleosynthesis depend strongly on the relevant neutron-capture and $β$-decay rates, as well as on the details of the stellar model being considered. Here, we have used a Monte-Carlo approach to evaluate the nu…
▽ More
The s-process, a production mechanism based on slow-neutron capture during stellar evolution, is the origin of about half the elements heavier than iron. Abundance predictions for s-process nucleosynthesis depend strongly on the relevant neutron-capture and $β$-decay rates, as well as on the details of the stellar model being considered. Here, we have used a Monte-Carlo approach to evaluate the nuclear uncertainty in s-process nucleosynthesis. We considered the helium burning of massive stars for the weak s-process and low-mass asymptotic-giant-branch stars for the main s-process. Our calculations include a realistic and general prescription for the temperature dependent uncertainty for the reaction cross sections. We find that the adopted uncertainty for (${\rm n},γ$) rates, tens of per cent on average, effects the production of s-process nuclei along the line of $β$-stability, and that the uncertainties in $β$-decay from excited state contributions, has the strongest impact on branching points.
△ Less
Submitted 15 February, 2018;
originally announced February 2018.
-
Uncertainties in the production of p nuclides in thermonuclear supernovae determined by Monte Carlo variations
Authors:
N. Nishimura,
T. Rauscher,
R. Hirschi,
A. St. J. Murphy,
G. Cescutti,
C. Travaglio
Abstract:
Thermonuclear supernovae originating from the explosion of a white dwarf accreting mass from a companion star have been suggested as a site for the production of $p$ nuclides. Such nuclei are produced during the explosion, in layers enriched with seed nuclei coming from prior strong $s$ processing. These seeds are transformed to proton-richer isotopes mainly by photodisintegration reactions. Sever…
▽ More
Thermonuclear supernovae originating from the explosion of a white dwarf accreting mass from a companion star have been suggested as a site for the production of $p$ nuclides. Such nuclei are produced during the explosion, in layers enriched with seed nuclei coming from prior strong $s$ processing. These seeds are transformed to proton-richer isotopes mainly by photodisintegration reactions. Several thousand trajectories from a 2D explosion model were used in a Monte Carlo approach. Temperature-dependent uncertainties were assigned individually to thousands of rates varied simultaneously in post-processing in an extended nuclear reaction network. The uncertainties in the final nuclear abundances originating from uncertainties in the astrophysical reaction rates were determined. In addition to the 35 classical $p$ nuclides, abundance uncertainties were also determined for the radioactive nuclides $^{92}$Nb, $^{97,98}$Tc, $^{146}$Sm, and for the abundance ratios $Y$(${}^{92}$Mo)/$Y$(${}^{94}$Mo), $Y$(${}^{92}$Nb)/$Y$(${}^{92}$Mo), $Y$(${}^{97}$Tc)/$Y$(${}^{98}$Ru), $Y$(${}^{98}$Tc)/$Y$(${}^{98}$Ru), and $Y$(${}^{146}$Sm)/$Y$(${}^{144}$Sm), important for Galactic Chemical Evolution studies. Uncertainties found were generally lower than a factor of two, although most nucleosynthesis flows mainly involve predicted rates with larger uncertainties. The main contribution to the total uncertainties comes from a group of trajectories with high peak density originating from the interior of the exploding white dwarf. The distinction between low-density and high-density trajectories allows more general conclusions to be drawn, also applicable to other simulations of white dwarf explosions.
△ Less
Submitted 11 December, 2017; v1 submitted 24 November, 2017;
originally announced November 2017.
-
Mass Ejection from the Remnant of a Binary Neutron Star Merger: Viscous-Radiation Hydrodynamics Study
Authors:
Sho Fujibayashi,
Kenta Kiuchi,
Nobuya Nishimura,
Yuichiro Sekiguchi,
Masaru Shibata
Abstract:
We perform long-term general relativistic neutrino radiation hydrodynamics simulations (in axisymmetry) for a massive neutron star (MNS) surrounded by a torus, which is a canonical remnant formed after the binary neutron star merger. We take into account the effects of viscosity, which is likely to arise in the merger remnant due to magnetohydrodynamical turbulence. As the initial condition, we em…
▽ More
We perform long-term general relativistic neutrino radiation hydrodynamics simulations (in axisymmetry) for a massive neutron star (MNS) surrounded by a torus, which is a canonical remnant formed after the binary neutron star merger. We take into account the effects of viscosity, which is likely to arise in the merger remnant due to magnetohydrodynamical turbulence. As the initial condition, we employ the azimuthally averaged data of the MNS-torus system derived in a three-dimensional, numerical-relativity simulation for the binary neutron star merger. The viscous effect plays key roles for the remnant evolution and mass ejection from it in two phases of the evolution. In the first $t\lesssim10$ ms, a differential rotation state of the MNS is changed to a rigidly rotating state, and as a result, a sound wave, which subsequently becomes a shock wave, is formed in the vicinity of the MNS due to the variation of the quasi-equilibrium state of the MNS. The shock wave induces significant mass ejection of mass $\sim(0.5-2.0)\times 10^{-2}M_\odot$ for the alpha viscosity parameter of $0.01-0.04$. For the longer-term evolution with $\sim 0.1-10$ s, a significant fraction of the torus material is ejected. The ejecta mass is likely to be of order $10^{-2}M_\odot$, so that the total mass of the viscosity-driven ejecta could dominate that of the dynamical ejecta of mass $\lesssim 10^{-2}M_\odot$. The electron fraction, $Y_e$, of the ejecta is always high enough ($Y_e\gtrsim0.25$) that this post-merger ejecta is lanthanide-poor; hence, the opacity of the ejecta is likely to be $\sim 10-100$ times lower than that of the dynamical ejecta. This indicates that the electromagnetic signal from the ejecta would be rapidly evolving, bright, and blue if it is observed from a small viewing angle ($\lesssim 45^\circ$) for which the effect of the dynamical ejecta is minor.
△ Less
Submitted 5 June, 2018; v1 submitted 6 November, 2017;
originally announced November 2017.
-
Abundance Uncertainties Obtained With the PizBuin Framework For Monte Carlo Reaction Rate Variations
Authors:
T. Rauscher,
N. Nishimura,
G. Cescutti,
R. Hirschi,
A. St. J. Murphy
Abstract:
Uncertainties in nucleosynthesis models originating from uncertainties in astrophysical reaction rates were estimated in a Monte Carlo variation procedure. Thousands of rates were simultaneously varied within individual, temperature-dependent errors to calculate their combined effect on final abundances. After a presentation of the method, results from application to three different nucleosynthesi…
▽ More
Uncertainties in nucleosynthesis models originating from uncertainties in astrophysical reaction rates were estimated in a Monte Carlo variation procedure. Thousands of rates were simultaneously varied within individual, temperature-dependent errors to calculate their combined effect on final abundances. After a presentation of the method, results from application to three different nucleosynthesis processes are shown: the $γ$-process and the s-process in massive stars, and the main s-process in AGB stars (preliminary results). Thermal excitation of nuclei in the stellar plasma and the combined action of several reactions increase the final uncertainties above the level of the experimental errors. The total uncertainty, on the other hand, remains within a factor of two even in processes involving a large number of unmeasured rates, with some notable exceptions for nuclides whose production is spread over several stellar layers and for s-process branchings.
△ Less
Submitted 3 September, 2017;
originally announced September 2017.
-
The s-process nucleosynthesis: impact of the uncertainties in the nuclear physics determined by Monte Carlo variations
Authors:
G. Cescutti,
N. Nishimura,
R. Hirschi,
T. Rauscher,
J. W. den Hartogh,
A. St. J. Murphy
Abstract:
We investigated the impact of uncertainties in neutron-capture and weak reactions (on heavy elements) on the s-process nucleosynthesis in low-mass stars and massive stars using a Monte-Carlo based approach. We performed extensive nuclear reaction network calculations that include newly evaluated temperature-dependent upper and lower limits for the individual reaction rates. We found beta-decay rat…
▽ More
We investigated the impact of uncertainties in neutron-capture and weak reactions (on heavy elements) on the s-process nucleosynthesis in low-mass stars and massive stars using a Monte-Carlo based approach. We performed extensive nuclear reaction network calculations that include newly evaluated temperature-dependent upper and lower limits for the individual reaction rates. We found beta-decay rate uncertainties affect only a few nuclei near s-process branchings, whereas most of the uncertainty in the final abundances is caused by uncertainties in the neutron capture rates. We suggest a list of uncertain rates as candidates for improved measurement by future experiments.
△ Less
Submitted 2 August, 2017;
originally announced August 2017.
-
Impacts of nuclear-physics uncertainty in stellar temperatures on the s-process nucleosynthesis
Authors:
N. Nishimura,
G. Cescutti,
R. Hirschi,
T. Rauscher,
J. Den Hartogh,
A. St. J. Murphy
Abstract:
We evaluated the uncertainty relevant to s-process nucleosynthesis using a Monte-Carlo centred approach. We are based on a realistic and general prescription of temperature dependent uncertainty for the reactions. We considered massive stars for the weak s-process and AGB stars for the main s-process. We found that the adopted uncertainty for (n,$γ$) rates, tens of per cent on average, affect the…
▽ More
We evaluated the uncertainty relevant to s-process nucleosynthesis using a Monte-Carlo centred approach. We are based on a realistic and general prescription of temperature dependent uncertainty for the reactions. We considered massive stars for the weak s-process and AGB stars for the main s-process. We found that the adopted uncertainty for (n,$γ$) rates, tens of per cent on average, affect the production of s-process nuclei along the $β$-stability line, while for $β$-decay, for which contributions from excited states enhances the uncertainty, has the strongest impact on branching points.
△ Less
Submitted 24 January, 2017;
originally announced January 2017.
-
Uncertainties in s-process nucleosynthesis in massive stars determined by Monte Carlo variations
Authors:
Nobuya Nishimura,
Raphael Hirschi,
Thomas Rauscher,
Alexander St. J. Murphy,
Gabriele Cescutti
Abstract:
The $s$-process in massive stars produces the weak component of the $s$-process (nuclei up to $A \sim 90$), in amounts that match solar abundances. For heavier isotopes, such as barium, production through neutron capture is significantly enhanced in very metal-poor stars with fast rotation. However, detailed theoretical predictions for the resulting final $s$-process abundances have important unce…
▽ More
The $s$-process in massive stars produces the weak component of the $s$-process (nuclei up to $A \sim 90$), in amounts that match solar abundances. For heavier isotopes, such as barium, production through neutron capture is significantly enhanced in very metal-poor stars with fast rotation. However, detailed theoretical predictions for the resulting final $s$-process abundances have important uncertainties caused both by the underlying uncertainties in the nuclear physics (principally neutron capture reaction and $β$-decay rates) as well as by the stellar evolution modeling. In this work, we investigated the impact of nuclear-physics uncertainties relevant to the $s$-process in massive stars. Using a Monte-Carlo based approach, we performed extensive nuclear reaction network calculations that include newly evaluated upper and lower limits for the individual temperature dependent reaction rates. We found that most of the uncertainty in the final abundances is caused by uncertainties in the neutron capture rates, while $β$-decay rate uncertainties affect only a few nuclei near $s$-process branchings. The $s$-process in rotating metal-poor stars shows quantitatively different uncertainties and key reactions, although the qualitative characteristics are similar. We confirmed that our results do not significantly change at different metallicities for fast rotating massive stars in the very low metallicity regime. We highlight which of the identified key reactions are realistic candidates for improved measurement by future experiments.
△ Less
Submitted 24 May, 2017; v1 submitted 2 January, 2017;
originally announced January 2017.
-
The intermediate r-process in core-collapse supernovae driven by the magneto-rotational instability
Authors:
Nobuya Nishimura,
Hidetomo Sawai,
Tomoya Takiwaki,
Shoichi Yamada,
Friedrich-Karl Thielemann
Abstract:
We investigated r-process nucleosynthesis in magneto-rotational supernovae, based on a new explosion mechanism induced by the magneto-rotational instability. A series of axisymmetric magneto-hydrodynamical simulations with detailed microphysics including neutrino heating is performed, numerically resolving the magneto-rotational instability. Neutrino-heating dominated explosions, enhanced by magne…
▽ More
We investigated r-process nucleosynthesis in magneto-rotational supernovae, based on a new explosion mechanism induced by the magneto-rotational instability. A series of axisymmetric magneto-hydrodynamical simulations with detailed microphysics including neutrino heating is performed, numerically resolving the magneto-rotational instability. Neutrino-heating dominated explosions, enhanced by magnetic fields, showed mildly neutron-rich ejecta producing nuclei up to $A \sim 130$ (i.e. the weak r-process), while explosion models with stronger magnetic fields reproduce a solar-like r-process pattern. More commonly seen abundance patterns in our models are in between the weak and regular r-process, producing lighter and intermediate mass nuclei. These {\it intermediate r-processes} exhibit a variety of abundance distributions, compatible with several abundance patterns in r-process-enhanced metal-poor stars. The amount of Eu ejecta $\sim 10^{-5} M_\odot$ in magnetically-driven jets agrees with predicted values in the chemical evolution of early galaxies. In contrast, neutrino-heating dominated explosions have a significant amount of Fe ($^{56}{\rm Ni}$) and Zn, comparable to regular supernovae and hypernovae, respectively. These results indicate magneto-rotational supernovae can produce a wide range of heavy nuclei from iron-group to r-process elements, depending on the explosion dynamics.
△ Less
Submitted 2 February, 2017; v1 submitted 7 November, 2016;
originally announced November 2016.
-
Production Uncertainties of p-Nuclei in the $γ$-Process in Massive Stars Using a Monte Carlo Approach
Authors:
T. Rauscher,
N. Nishimura,
R. Hirschi,
G. Cescutti,
A. St. J. Murphy,
A. Heger
Abstract:
Proton-rich nuclei, the so-called p-nuclei, are made in photodisintegration processes in outer shells of massive stars in the course of the final supernova explosion. Nuclear uncertainties in the production of these nuclei have been quantified in a Monte Carlo procedure. Bespoke temperature-dependent uncertainties were assigned to different types of reactions involving nuclei from Fe to Bi and all…
▽ More
Proton-rich nuclei, the so-called p-nuclei, are made in photodisintegration processes in outer shells of massive stars in the course of the final supernova explosion. Nuclear uncertainties in the production of these nuclei have been quantified in a Monte Carlo procedure. Bespoke temperature-dependent uncertainties were assigned to different types of reactions involving nuclei from Fe to Bi and all rates were varied randomly within the uncertainties. The resulting total production uncertainties of the p-nuclei are below a factor of two, with few exceptions. Key reactions dominating the final uncertainties have been identified in an automated procedure using correlations between rate and abundance uncertainties. Our results are compared to those of a previous study manually varying reaction rates.
△ Less
Submitted 14 October, 2016;
originally announced October 2016.
-
Uncertainties in the production of $p$ nuclei in massive stars obtained from Monte Carlo variations
Authors:
T. Rauscher,
N. Nishimura,
R. Hirschi,
G. Cescutti,
A. St. J. Murphy,
A. Heger
Abstract:
Nuclear uncertainties in the production of $p$ nuclei in massive stars have been quantified in a Monte Carlo procedure. Bespoke temperature-dependent uncertainties were assigned to different types of reactions involving nuclei from Fe to Bi. Their simultaneous impact was studied in postprocessing explosive trajectories for three different stellar models. It was found that the grid of mass zones in…
▽ More
Nuclear uncertainties in the production of $p$ nuclei in massive stars have been quantified in a Monte Carlo procedure. Bespoke temperature-dependent uncertainties were assigned to different types of reactions involving nuclei from Fe to Bi. Their simultaneous impact was studied in postprocessing explosive trajectories for three different stellar models. It was found that the grid of mass zones in the model of a 25 $M_\odot$ star, which is widely used for investigations of $p$ nucleosynthesis, is too crude to properly resolve the detailed temperature changes required for describing the production of $p$ nuclei. Using models with finer grids for 15 $M_\odot$ and 25 $M_\odot$ stars with initial solar metallicity, it was found that most of the production uncertainties introduced by nuclear reaction uncertainties are smaller than a factor of two. Since a large number of rates were varied at the same time in the Monte Carlo procedure, possible cancellation effects of several uncertainties could be taken into account. Key rates were identified for each $p$ nucleus, which provide the dominant contribution to the production uncertainty. These key rates were found by examining correlations between rate variations and resulting abundance changes. This method is superior to studying flow patterns, especially when the flows are complex, and to individual, sequential variation of a few rates.
△ Less
Submitted 6 September, 2016; v1 submitted 17 June, 2016;
originally announced June 2016.
-
Impact of the first-forbidden $β$ decay on the production of $A \sim 195$ r-process peak
Authors:
Nobuya Nishimura,
Zsolt Podolyák,
Dong-Liang Fang,
Toshio Suzuki
Abstract:
We investigated the effects of first-forbidden transitions in $β$ decays on the production of the r-process $A \sim 195$ peak. The theoretical calculated $β$-decay rates with $β$-delayed neutron emission were examined using several astrophysical conditions. As the first-borbidden decay is dominant in $N \sim 126$ neutron-rich nuclei, their inclusion shortens $β$-decay lifetimes and shifts the abun…
▽ More
We investigated the effects of first-forbidden transitions in $β$ decays on the production of the r-process $A \sim 195$ peak. The theoretical calculated $β$-decay rates with $β$-delayed neutron emission were examined using several astrophysical conditions. As the first-borbidden decay is dominant in $N \sim 126$ neutron-rich nuclei, their inclusion shortens $β$-decay lifetimes and shifts the abundance peak towards higher masses. Additionally, the inclusion of the $β$-delayed neutron emission results in a wider abundance peak, and smoothens the mass distribution by removing the odd-even mass staggering. The effects are commonly seen in the results of all adopted astrophysical models. Nevertheless there are quantitative differences, indicating that remaining uncertainty in the determination of half-lives for $N=126$ nuclei is still significant in order to determine the production of the r-process peak.
△ Less
Submitted 16 March, 2016;
originally announced March 2016.
-
Theory Considerations For Nucleosynthesis Beyond Fe With Special Emphasis On p-Nuclei In Massive Stars
Authors:
T. Rauscher,
N. Nishimura,
R. Hirschi
Abstract:
Nucleosynthesis of heavy elements requires the use of different experimental and theoretical methods to determine astrophysical reaction rates than light element nucleosynthesis. Additionally, there are also larger uncertainties involved in the astrophysical models, both because the sites are not well known and because of differing numerical treatments in different models. As an example for the la…
▽ More
Nucleosynthesis of heavy elements requires the use of different experimental and theoretical methods to determine astrophysical reaction rates than light element nucleosynthesis. Additionally, there are also larger uncertainties involved in the astrophysical models, both because the sites are not well known and because of differing numerical treatments in different models. As an example for the latter, the production of p-nuclei is compared in two different stellar models, demonstrating that a model widely used for postproduction calculations may have a zone grid too coarse to follow the synthesis of p-nuclei in detail.
△ Less
Submitted 24 October, 2015; v1 submitted 16 October, 2015;
originally announced October 2015.
-
The r-process in Magnetorotational Supernovae
Authors:
Takuji Tsujimoto,
Nobuya Nishimura
Abstract:
One of the hottest open issues involving the evolution of r-process elements is fast enrichment in the early Universe. Clear evidence for the chemical enrichment of r-process elements is seen in the stellar abundances of extremely metal poor stars in the Galactic halo. However, small-mass galaxies are the ideal testbed for studying the evolutionary features of r-process enrichment given the potent…
▽ More
One of the hottest open issues involving the evolution of r-process elements is fast enrichment in the early Universe. Clear evidence for the chemical enrichment of r-process elements is seen in the stellar abundances of extremely metal poor stars in the Galactic halo. However, small-mass galaxies are the ideal testbed for studying the evolutionary features of r-process enrichment given the potential rarity of production events yielding heavy r-process elements. Their occurrences become countable and thus an enrichment path due to each event can be found in the stellar abundances. We examine the chemical feature of Eu abundance at an early stage of ${\rm[Fe/H]} \lesssim -2$ in the Draco and Sculptor dwarf spheroidal (dSph) galaxies. Accordingly, we constrain the properties of the Eu production in the early dSphs. We find that the Draco dSph experienced a few Eu production events, whereas Eu enrichment took place more continuously in the Sculptor dSph due to its larger stellar mass. The event rate of Eu production is estimated to be about one per $100$-$200$ core-collapse supernovae, and a Eu mass of $\sim (1-2) \times 10^{-5}$\ms per single event is deduced by associating this frequency with the observed plateau value of ${\rm [Eu/H]} \sim -1.3$ for ${\rm [Fe/H]} \gtrsim -2$. The observed plateau implies that early Eu enrichment ceases at ${\rm [Fe/H]} \approx -2$. Such a selective operation only in low-metallicity stars supports magnetorotational supernovae, which require very fast rotation, as the site of early Eu production. We show that the Eu yields deduced from chemical evolution agree well with the nucleosynthesis results from corresponding supernovae models.
△ Less
Submitted 24 September, 2015; v1 submitted 31 August, 2015;
originally announced September 2015.
-
The r-process nucleosynthesis in the various jet-like explosions of magnetorotational core-collapse supernovae
Authors:
Nobuya Nishimura,
Tomoya Takiwaki,
Friedrich-Karl Thielemann
Abstract:
The r-process nucleosynthesis in core-collapse supernovae (CC-SNe) is studied, with a focus on the explosion scenario induced by rotation and strong magnetic fields. Nucleosynthesis calculations are conducted based on magneto-hydrodynamical explosion models with a wide range of parameters for initial rotation and magnetic fields. The explosion models are classified in two different types: i.e., pr…
▽ More
The r-process nucleosynthesis in core-collapse supernovae (CC-SNe) is studied, with a focus on the explosion scenario induced by rotation and strong magnetic fields. Nucleosynthesis calculations are conducted based on magneto-hydrodynamical explosion models with a wide range of parameters for initial rotation and magnetic fields. The explosion models are classified in two different types: i.e., prompt-magnetic-jet and delayed-magnetic-jet, for which the magnetic fields of proto-neutron stars (PNSs) during collapse and the core-bounce are strong and comparatively moderate, respectively. Following the hydrodynamical trajectories of each explosion model, we confirmed that r-processes successfully occur in the prompt-magnetic-jets, which produce heavy nuclei including actinides. On the other hand, the r-process in the delayed-magnetic-jet is suppressed, which synthesizes only nuclei up to the second peak ($A \sim 130$). Thus, the r-process in the delayed-magnetic-jets could explain only "weak r-process" patterns observed in metal-poor stars rather than the "main r-process", represented by the solar abundances. Our results imply that core-collapse supernovae are possible astronomical sources of heavy r-process elements if their magnetic fields are strong enough, while weaker magnetic explosions may produce "weak r-process" patterns ($A \lesssim 130$). We show the potential importance and necessity of magneto-rotational supernovae for explaining the galactic chemical evolution, as well as abundances of r-process enhanced metal-poor stars. We also examine the effects of the remaining uncertainties in the nature of PNSs due to weak interactions that determine the final neutron-richness of ejecta. Additionally, we briefly discuss radioactive isotope yields in primary jets (e.g., $^{56}$Ni), with relation to several optical observation of SNe and relevant high-energy astronomical phenomena.
△ Less
Submitted 31 August, 2015; v1 submitted 26 January, 2015;
originally announced January 2015.
-
Production of 92Nb, 92Mo, and 146Sm in the gamma-process in SNIa
Authors:
T. Rauscher,
C. Travaglio,
R. Gallino,
N. Nishimura,
R. Hirschi
Abstract:
The knowledge of the production of extinct radioactivities like 92Nb and 146Sm by photodisintegration processes in ccSN and SNIa models is essential for interpreting abundances in meteoritic material and for Galactic Chemical Evolution (GCE). The 92Mo/92Nb and 146Sm/144Sm ratios provide constraints for GCE and production sites. We present results for SNIa with emphasis on nuclear uncertainties.
The knowledge of the production of extinct radioactivities like 92Nb and 146Sm by photodisintegration processes in ccSN and SNIa models is essential for interpreting abundances in meteoritic material and for Galactic Chemical Evolution (GCE). The 92Mo/92Nb and 146Sm/144Sm ratios provide constraints for GCE and production sites. We present results for SNIa with emphasis on nuclear uncertainties.
△ Less
Submitted 22 December, 2014;
originally announced December 2014.
-
Code dependencies of pre-supernova evolution and nucleosynthesis in massive stars: Evolution to the end of core helium burning
Authors:
Samuel Jones,
Raphael Hirschi,
Marco Pignatari,
Alexander Heger,
Cyril Georgy,
Nobuya Nishimura,
Chris Fryer,
Falk Herwig
Abstract:
Massive stars are key sources of radiative, kinetic, and chemical feedback in the universe. Grids of massive star models computed by different groups each using their own codes, input physics choices and numerical approximations, however, lead to inconsistent results for the same stars. We use three of these 1D codes---GENEC, KEPLER and MESA---to compute non-rotating stellar models of…
▽ More
Massive stars are key sources of radiative, kinetic, and chemical feedback in the universe. Grids of massive star models computed by different groups each using their own codes, input physics choices and numerical approximations, however, lead to inconsistent results for the same stars. We use three of these 1D codes---GENEC, KEPLER and MESA---to compute non-rotating stellar models of $15~\mathrm{M}_\odot$, $20~\mathrm{M}_\odot$, and $25~\mathrm{M}_\odot$ and compare their nucleosynthesis. We follow the evolution from the main sequence until the end of core helium burning. The GENEC and KEPLER models hold physics assumptions used in large grids of published models. The MESA code was set up to use convective core overshooting such that the CO core masses are consistent with those obtained by GENEC. For all models, full nucleosynthesis is computed using the NuGrid post-processing tool MPPNP.
We find that the surface abundances predicted by the models are in reasonable agreement. In the helium core, the standard deviation of the elemental overproduction factors for Fe to Mo is less than $30\,\%$---smaller than the impact of the present nuclear physics uncertainties. For our three initial masses, the three stellar evolution codes yield consistent results. Differences in key properties of the models, e.g., helium and CO core masses and the time spent as a red supergiant, are traced back to the treatment of convection and, to a lesser extent, mass loss. The mixing processes in stars remain the key uncertainty in stellar modelling. Better constrained prescriptions are thus necessary to improve the predictive power of stellar evolution models.
△ Less
Submitted 19 December, 2014;
originally announced December 2014.
-
Production of all the r-process nuclides in the dynamical ejecta of neutron star mergers
Authors:
Shinya Wanajo,
Yuichiro Sekiguchi,
Nobuya Nishimura,
Kenta Kiuchi,
Koutarou Kyutoku,
Masaru Shibata
Abstract:
Recent studies suggest that binary neutron star (NS-NS) mergers robustly produce the heavy r-process nuclei above the atomic mass number A ~ 130 because of their ejecta consisting of almost pure neutrons (electron fraction of Y_e < 0.1). However, little production of the lighter r-process nuclei (A = 90-120) conflicts with the spectroscopic results of r-process-enhanced Galactic halo stars. We pre…
▽ More
Recent studies suggest that binary neutron star (NS-NS) mergers robustly produce the heavy r-process nuclei above the atomic mass number A ~ 130 because of their ejecta consisting of almost pure neutrons (electron fraction of Y_e < 0.1). However, little production of the lighter r-process nuclei (A = 90-120) conflicts with the spectroscopic results of r-process-enhanced Galactic halo stars. We present, for the first time, the result of nucleosynthesis calculations based on the fully general-relativistic simulation of a NS-NS merger with approximate neutrino transport. It is found that the bulk of the dynamical ejecta are appreciably shock-heated and neutrino-processed, resulting in a wide range of Y_e (= 0.09-0.45). The mass-averaged abundance distribution of calculated nucleosynthesis yields is in reasonable agreement with the full-mass range (A = 90-240) of the solar r-process curve. This implies, if our model is representative of such events, that the dynamical ejecta of NS-NS mergers can be the origin of the Galactic r-process nuclei. Our result also shows that the radioactive heating after ~ 1 day from the merging, giving rise to r-process-powered transient emission, is dominated by the beta-decays of several species close to stability with precisely measured half-lives. This implies that the total radioactive heating rate for such an event can be well constrained within about a factor of two if the ejected material has a solar-like r-process pattern.
△ Less
Submitted 10 June, 2014; v1 submitted 28 February, 2014;
originally announced February 2014.
-
Revisiting Impacts of Nuclear Burning for Reviving Weak Shocks in Neutrino-Driven Supernovae
Authors:
Ko Nakamura,
Tomoya Takiwaki,
Kei Kotake,
Nobuya Nishimura
Abstract:
We revisit potential impacts of nuclear burning on the onset of the neutrino-driven explosions of core-collapse supernovae. By changing the neutrino luminosity and its decay time to obtain parametric explosions in one-(1D) and two-dimensional (2D) models with or without a 13-isotope alpha network, we study how the inclusion of nuclear burning could affect the postbounce dynamics for four progenito…
▽ More
We revisit potential impacts of nuclear burning on the onset of the neutrino-driven explosions of core-collapse supernovae. By changing the neutrino luminosity and its decay time to obtain parametric explosions in one-(1D) and two-dimensional (2D) models with or without a 13-isotope alpha network, we study how the inclusion of nuclear burning could affect the postbounce dynamics for four progenitor models; three for 15.0 Msun stars, one for an 11.2 Msun star. We find that the energy supply due to nuclear burning of infalling material behind the shock can energize the shock expansion especially for models that produce only marginal explosions in the absence of nuclear burning. These models are energized by nuclear energy deposition when the shock front passes through the silicon-rich layer and/or later it touches the oxygen-rich layer. Depending on the neutrino luminosity and its decay time, a diagnostic energy of explosion increases up to a few times 10^50 erg for models with nuclear burning compared to the corresponding models without. We point out that these features are most remarkable for the Limongi-Chieffi progenitor in both 1D and 2D, because the progenitor model possesses a massive oxygen layer with its inner-edge radius being smallest among the employed progenitors, so that the shock can touch the rich fuel on a shorter timescale after bounce. The energy difference is generally smaller (~0.1-0.2 times 10^51 erg) in 2D than in 1D (at most ~0.6 times 10^51 erg). This is because neutrino-driven convection and the shock instability in 2D models enhance the neutrino heating efficiency, which makes the contribution of nuclear burning relatively smaller compared to 1D models. Considering uncertainties in progenitor models, our results indicate that nuclear burning should remain as one of the important ingredients to foster the onset of neutrino-driven explosions.
△ Less
Submitted 8 January, 2014; v1 submitted 25 July, 2012;
originally announced July 2012.
-
Impact of New beta-decay Half-lives on r-process Nucleosynthesis
Authors:
Nobuya Nishimura,
Toshitaka Kajino,
Grant J. Mathews,
Shunji Nishimura,
Toshio Suzuki
Abstract:
We investigate the effects of newly measured beta-decay half-lives on r-process nucleosynthesis. These new rates were determined by recent experiments at the radioactive isotope beam factory facility in the RIKEN Nishina Center. We adopt an r-process nucleosynthesis environment based on a magnetohydrodynamic supernova explosion model that includes strong magnetic fields and rapid rotation of the p…
▽ More
We investigate the effects of newly measured beta-decay half-lives on r-process nucleosynthesis. These new rates were determined by recent experiments at the radioactive isotope beam factory facility in the RIKEN Nishina Center. We adopt an r-process nucleosynthesis environment based on a magnetohydrodynamic supernova explosion model that includes strong magnetic fields and rapid rotation of the progenitor. A number of the new beta-decay rates are for nuclei on or near the r-process path, and hence they affect the nucleosynthesis yields and time scale of the r-process. The main effect of the newly measured beta-decay half-lives is an enhancement in the calculated abundance of isotopes with mass number A = 110 -- 120 relative to calculated abundances based upon beta-decay rates estimated with the finite-range droplet mass model. This effect slightly alleviates, but does not fully explain, the tendency of r-process models to underproduce isotopes with A = 110 -- 120 compared to the solar-system r-process abundances.
△ Less
Submitted 4 April, 2012; v1 submitted 23 March, 2012;
originally announced March 2012.
-
Magnetorotationally driven Supernovae as the origin of early galaxy $r$-process elements?
Authors:
Christian Winteler,
Roger Kaeppeli,
Albino Perego,
Almudena Arcones,
Nicolas Vasset,
Nobuya Nishimura,
Matthias Liebendoerfer,
Friedrich-Karl Thielemann
Abstract:
We examine magnetorotationally driven supernovae as sources of $r$-process elements in the early Galaxy. On the basis of thermodynamic histories of tracer particles from a three-dimensional magnetohydrodynamical core-collapse supernova model with approximated neutrino transport, we perform nucleosynthesis calculations with and without considering the effects of neutrino absorption reactions on the…
▽ More
We examine magnetorotationally driven supernovae as sources of $r$-process elements in the early Galaxy. On the basis of thermodynamic histories of tracer particles from a three-dimensional magnetohydrodynamical core-collapse supernova model with approximated neutrino transport, we perform nucleosynthesis calculations with and without considering the effects of neutrino absorption reactions on the electron fraction ($Y_{e}$) during post-processing. We find that the peak distribution of $Y_{e}$ in the ejecta is shifted from $\sim0.15$ to $\sim0.17$ and broadened toward higher $Y_{e}$ due to neutrino absorption. Nevertheless, in both cases the second and third peaks of the solar $r$-process element distribution can be well reproduced. The rare progenitor configuration that was used here, characterized by a high rotation rate and a large magnetic field necessary for the formation of bipolar jets, could naturally provide a site for the strong $r$-process in agreement with observations of the early galactic chemical evolution.
△ Less
Submitted 14 April, 2012; v1 submitted 3 March, 2012;
originally announced March 2012.
-
Nucleosynthesis in core-collapse supernova explosions triggered by a quark-hadron phase transition
Authors:
Nobuya Nishimura,
Tobias Fischer,
Friedrich-Karl Thielemann,
Carla Fröhlich,
Matthias Hempel,
Roger Käppeli,
Gabriel Martínez-Pinedo,
Thomas Rauscher,
Irina Sagert,
Christian Winteler
Abstract:
We explore heavy element nucleosynthesis in the explosion of massive stars which are triggered by a quark-hadron phase transition during the early post bounce phase of core-collapse supernovae. The present study is based on general relativistic radiation hydrodynamics simulations with three-flavor Boltzmann neutrino transport in spherical symmetry, which utilize a quark-hadron hybrid equation of s…
▽ More
We explore heavy element nucleosynthesis in the explosion of massive stars which are triggered by a quark-hadron phase transition during the early post bounce phase of core-collapse supernovae. The present study is based on general relativistic radiation hydrodynamics simulations with three-flavor Boltzmann neutrino transport in spherical symmetry, which utilize a quark-hadron hybrid equation of state based on the MIT bag model for strange quark matter. The quark-hadron phase transition inside the stellar core forms a shock wave propagating towards the surface of the proto-neutron star. The shock wave results in an explosion and ejects neutron-rich matter which is piled up or accreting on the proto-neutron star. Later, during the cooling phase, the proto-neutron star develops a proton-rich neutrino-driven wind. We present a detailed analysis of the nucleosynthesis outcome in both neutron-rich and proton-rich ejecta and compare our integrated nucleosynthesis with observations of metal poor stars.
△ Less
Submitted 24 August, 2012; v1 submitted 23 December, 2011;
originally announced December 2011.
-
Observational Constraint on Heavy Element Production in Inhomogeneous Big Bang Nucleosynthesis
Authors:
Riou Nakamura,
Masa-aki Hashimoto,
Shin-ichiro Fujimoto,
Nobuya Nishimura,
Katsuhiko Sato
Abstract:
Based on a scenario of the inhomogeneous big-bang nucleosynthesis (IBBN), we investigate the detailed nucleosynthesis that includes the production of heavy elements beyond Li-7. From the observational constraints on light elements of He4 and D for the baryon-to-photon ratio given by WMAP, possible regions found on the plane of the volume fraction of the high density region against the ratio betwee…
▽ More
Based on a scenario of the inhomogeneous big-bang nucleosynthesis (IBBN), we investigate the detailed nucleosynthesis that includes the production of heavy elements beyond Li-7. From the observational constraints on light elements of He4 and D for the baryon-to-photon ratio given by WMAP, possible regions found on the plane of the volume fraction of the high density region against the ratio between high- and low-density regions.
In these allowed regions, we have confirmed that the heavy elements beyond Fe can be produced appreciably, where p- and/or r-process elements are produced well simultaneously compared to the solar system abundances. We suggest that recent observational signals such as He4 overabundance in globular clusters and high metallicity abundances in quasars could be partly due to the results of IBBN. Possible implications are given for the formation of the first generation stars
△ Less
Submitted 4 July, 2013; v1 submitted 2 July, 2010;
originally announced July 2010.
-
Nucleosynthesis in Magnetically Driven Jets from Collapsars
Authors:
Shin-ichiro Fujimoto,
Nobuya Nishimura,
Masa-aki Hashimoto
Abstract:
We have made detailed calculations of the composition of magnetically driven jets ejected from collapsars, or rapidly rotating massive stars, based on long-term magnetohydrodynamic simulations of their core collapse with various distributions of magnetic field and angular momentum before collapse. We follow the evolution of the abundances of about 4000 nuclides from the collapse phase to the eje…
▽ More
We have made detailed calculations of the composition of magnetically driven jets ejected from collapsars, or rapidly rotating massive stars, based on long-term magnetohydrodynamic simulations of their core collapse with various distributions of magnetic field and angular momentum before collapse. We follow the evolution of the abundances of about 4000 nuclides from the collapse phase to the ejection phase and through the jet generation phase using a large nuclear reaction network. We find that the r-process successfully operates only in energetic jets (> 1e51 ergs), such that U and Th are synthesized abundantly, even when the collapsar has a relatively weak magnetic field (1e10 G) and a moderately rotating core before the collapse. The abundance patterns inside the jets are similar to those of the r-elements in the solar system. About 0.01-0.06 Msun neutron-rich, heavy nuclei are ejected from a collapsar with energetic jets. The higher energy jets have larger amounts of Ni56, varying from 0.00037 to 0.06Msun. Less energetic jets, which eject small amounts of Ni56, could induce a gamma-ray burst (GRB) a supernova, such as GRB 060505 or GRB 060614. Considerable amounts of r-elements are likely to be ejected from GRBs with hypernovae, if both the GRB and hypernova are induced by jets that are driven near the black hole.
△ Less
Submitted 7 April, 2008;
originally announced April 2008.
-
R-Process Nucleosynthesis in MHD Jet Explosions of Core-Collapse Supernovae
Authors:
Sunao Nishimura,
Kei Kotake,
Masa-aki Hashimoto,
Shoichi Yamada,
Nobuya Nishimura,
Shinichiro Fujimoto,
Katsuhiko Sato
Abstract:
We investigate $r$-process nucleosynthesis during the magnetohydrodynamical (MHD) explosion of supernova in a massive star of 13 $M_{\odot}$. Contrary to the case of the spherical explosion, jet-like explosion due to the combined effects of the rotation and magnetic field lowers the electron fraction significantly inside the layers above the iron core. We find that the ejected material of low el…
▽ More
We investigate $r$-process nucleosynthesis during the magnetohydrodynamical (MHD) explosion of supernova in a massive star of 13 $M_{\odot}$. Contrary to the case of the spherical explosion, jet-like explosion due to the combined effects of the rotation and magnetic field lowers the electron fraction significantly inside the layers above the iron core. We find that the ejected material of low electron fraction responsible for the $r$-process comes out from the silicon rich layer of the presupernova model. This leads to the production up to the third peak in the solar $r$-process elements. We examine whether the fission affects the $r$-process paths by using the full nuclear reaction network with both the spontaneous and $β$-delayed fission included. Moreover, we pay particular attention how the mass formula affects the $r$-process peaks with use of two mass formulae. It is found that both formulae can reproduce the global abundance pattern up to the third peak though detailed distributions are rather different. We point out that there are variations in the $r$-process nucleosynthesis if the MHD effects play an important role in the supernova explosion.
△ Less
Submitted 5 April, 2005;
originally announced April 2005.