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.

The $\nu 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 $\nu 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 $\nu 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$\alpha$ and $^{56}$Ni(n,p)$^{56}$Co rates at temperatures $T>3$ GK.

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.

It has been suggested that a $\nu$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{\nu}_\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 $\nu$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 $\nu$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.

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%.

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.

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.

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.

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 $\beta$-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,$\gamma$), $^{64}$Ni(n,$\gamma$), and $^{138}$Ba(n,$\gamma$). We discuss the prospect of reducing uncertainties in the key reactions identified in this study with future experiments.

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 $\beta$-decay.

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 $\beta$-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},\gamma$) rates, tens of per cent on average, effects the production of s-process nuclei along the line of $\beta$-stability, and that the uncertainties in $\beta$-decay from excited state contributions, has the strongest impact on branching points.

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.

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 $\gamma$-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.

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.

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 $\beta$-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 $\beta$-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.

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.

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,$\gamma$) rates, tens of per cent on average, affect the production of s-process nuclei along the $\beta$-stability line, while for $\beta$-decay, for which contributions from excited states enhances the uncertainty, has the strongest impact on branching points.

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.

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.

We investigated the effects of first-forbidden transitions in $\beta$ decays on the production of the r-process $A \sim 195$ peak. The theoretical calculated $\beta$-decay rates with $\beta$-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 $\beta$-decay lifetimes and shifts the abundance peak towards higher masses. Additionally, the inclusion of the $\beta$-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.

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.

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.

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.

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.

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.

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.

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.

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

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.

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.

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.

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.

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 $\beta$-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.