text:
# THE 𝑟-PROCESS IN SUPERNOVA EXPLOSIONS FROM THE COLLAPSE OF O-Ne-Mg CORES

## 1 INTRODUCTION

The astrophysical origin of the rapid neutron-capture ($`r`$-process) species has been a long-standing mystery. Recently, however, a number of important new clues have been provided by spectroscopic studies of extremely-metal-poor stars in the Galaxy. The appearance of neutron-capture elements in these oldest stars in the Galaxy, including the pure-$`r`$-process origin of elements such as thorium and uranium, strongly suggests that the $`r`$-process nuclei have come from core-collapse supernovae (Sneden et al., 1996, 2000, 2003; Cayrel et al., 2001; Cowan et al., 2002; Hill et al., 2002). Ishimaru & Wanajo (1999) have shown that the large star-to-star dispersion of the observed abundances of neutron-capture elements relative to iron in very metal-poor stars is also naturally explained if the $`r`$-process elements originate from a limited mass range of core-collapse supernovae with little iron production ($`810M_{}`$ or $`30M_{}`$). Qian & Wasserburg (2003; see also Wasserburg & Qian 2000; Qian & Wasserburg 2001, 2002) have proposed that the production of the heavy $`r`$-process nuclei ($`A>130`$) is decoupled from the production of iron-peak and $`\alpha `$ nuclei by comparison of observed abundances among the extremely metal-poor stars. In this view, the production site of the heavy $`r`$-process nuclei is associated with the accretion-induced collapse (AIC) of a (carbon-oxygen or oxygen-neon-magnesium) white dwarf in a binary system (Nomoto & Kondo, 1991), or Type II supernovae from $`810M_{}`$ stars (Nomoto, 1984).

So far, the “neutrino wind” scenario, in which the free nucleons accelerated by the intense neutrino flux near the neutrino sphere of a core-collapse supernova assemble to heavier nuclei, has been believed to be the most promising astrophysical site of the $`r`$-process (Woosley et al., 1994). Even this scenario, however, encounters some difficulties (Qian & Woosley, 1996; Hoffman et al., 1997; Cardall & Fuller, 1997; Otsuki et al., 2000; Thompson, Burrows, & Meyer, 2001; Wanajo et al., 2001, 2002). For example, Wanajo et al. (2001) have shown that an extremely compact proto-neutron star, e.g., $`2.0M_{}`$ and 10 km, must be formed in order to account for the solar $`r`$-process pattern, at least within the framework of a spherically symmetric explosion. Although its possibility cannot be excluded, such a compact remnant is allowed by only a few of the many existing equations of state of nuclear matter.

In addition, recent spectroscopic studies of extremely metal-poor stars in the Galactic halo indicate that the observed abundance patterns of the lighter ($`Z<56`$) and heavier ($`Z>56`$) neutron-capture elements cannot be explained by a single astrophysical site (e.g., neutrino winds); there must exist at least two different $`r`$-process sites (Ishimaru & Wanajo, 2000; Qian & Wasserburg, 2001; Johnson & Bolte, 2002; Sneden & Cowan, 2003). Hence, it is of special importance to consider alternative possibilities for the occurrence of the $`r`$-process in core-collapse supernovae.

Nomoto (1984, 1987) has shown that $`810M_{}`$ stars form an electron-degenerate oxygen-neon-magnesium (O-Ne-Mg) core that does not undergo further nuclear burning; rather, it directly collapses due to electron capture on 24Mg and 20Ne (Miyaji et al. 1980). Hillebrandt, Nomoto, & Wolff (1984) have demonstrated that the collapsing O-Ne-Mg core explodes in a prompt manner, and Wheeler, Cowan, & Hillebrandt (1998) have suggested that the exploding O-Ne-Mg core could be a viable site for the $`r`$-process. It has been pointed out that, if the core exploded hydrodynamically prior to the onset of delayed neutrino heating (i.e., it underwent a prompt explosion), the electron fraction (electron number per baryon), $`Y_e`$, in the innermost layer of the ejecta would approach $`0.2`$ (Hillebrandt, Nomoto, & Wolff, 1984). Earlier works have in fact shown that a robust $`r`$-process proceeds in such conditions (Schramm, 1973; Sato, 1974; Hillebrandt, Takahashi, & Kodama, 1976).

Recently, Sumiyoshi et al. (2001) have demonstrated that the prompt explosion of an $`11M_{}`$ star with an iron core might also be a promising site of the $`r`$-process. They obtained a prompt explosion by an adiabatic core-collapse calculation without inclusion of electron capture and neutrino transport. Thompson, Burrows, & Pinto (2003) have shown, however, that no explosion is obtained with such a progenitor (Woosley & Weaver, 1995), when including electron capture along with a detailed treatment of neutrino transport. Many previous works have suggested that even the lowest mass of core-collapse supernovae ($`10M_{}`$), in which a relatively smaller iron core is formed, may have difficulties in achieving a hydrodynamic explosion (Bowers & Wilson, 1982; Burrows & Lattimer, 1983, 1985; Bruenn, 1989a, b; Baron & Cooperstein, 1990).

On the other hand, the question of whether $`810M_{}`$ stars that form O-Ne-Mg cores can explode hydrodynamically is still open. The possibility that these stars explode promptly remains because of the smaller iron core present at the onset of the core bounce, as well as the smaller gravitational potential of their collapsing cores. Hillebrandt, Nomoto, & Wolff (1984) have obtained a prompt explosion of a $`9M_{}`$ star with a $`1.38M_{}`$ O-Ne-Mg core, while others, using the same progenitor, have not (Burrows & Lattimer, 1985; Baron, Cooperstein, & Kahana, 1987). Mayle & Wilson (1988) obtained an explosion, not by a prompt shock, but by late-time neutrino heating. Similar results can be seen in the studies of AICs. Note that an AIC is an analogous phenomenon to a collapsing O-Ne-Mg core resulted from a single $`810M_{}`$ star, since both consist of electron degenerate cores (Nomoto & Kondo, 1991). Fryer et al. (1999) have obtained an explosion by neutrino heating, while others have not (Baron et al., 1987; Woosley & Baron, 1992). The reason for these different outcomes is due, perhaps, to the application of different equations of state for dense matter, although other physical inputs may also have some influence (Fryer et al., 1999). Thus, even if a star of $`810M_{}`$ exploded, it would be difficult to derive, with confidence, the physical properties as well as the mass of the ejected matter. Given this highly uncertain situation it is necessary to examine the resulting $`r`$-process nucleosynthesis in explosions obtained with different sets of input physics.

The purpose of this study is to investigate conditions necessary for the production of $`r`$-process nuclei obtained in purely hydrodynamical models of prompt explosions of collapsing O-Ne-Mg cores, and to explore some of the consequences if those conditions are met. The core collapse and the subsequent core bounce are simulated by a one-dimensional hydrodynamic code with Newtonian gravity (§ 2). For simplicity, neutrino transport is not taken into account. Hence, we focus on only hydrodynamical explosions just after core bounce, not on the delayed explosions obtained via late-time neutrino heating. As seen in § 2, the explosion is marginal, and no $`r`$-processing is expected. In order to obtain $`r`$-processed material, we find it necessary to force the occurrence of more energetic explosions. The energetic explosions are simulated by artificial enhancements of the shock-heating energy, rather than by application of different sets of input physics, for simplicity. The $`r`$-process nucleosynthesis in these explosions is then calculated with the use of a nuclear reaction network code (§ 3). The resulting contribution of the $`r`$-process material created in these simulations to the early chemical evolution of the Galaxy is discussed in § 4. The results of chronometric age dating, using the U-Th chronometer pair based on our nucleosynthesis results, is discussed in § 5. A summary and conclusions follow in § 6.

## 2 PROMPT EXPLOSION

A pre-supernova model of a $`9M_{}`$ star is taken from Nomoto (1984), which forms a 1.38 $`M_{}`$ O-Ne-Mg core near the end of its evolution (see also Miyaji et al., 1980; Miyaji & Nomoto, 1987; Nomoto, 1987). We link this core to a one-dimensional implicit Lagrangian hydrodynamic code with Newtonian gravity (Bowers & Wilson, 1991). This core is modeled with a finely zoned mesh of 200 mass shells ($`2\times 10^2M_{}`$ to $`0.8M_{}`$, $`5\times 10^3M_{}`$ to $`1.3M_{}`$, and $`5\times 10^31\times 10^7M_{}`$ to the edge of the core).

The equation of state of nuclear matter (EOS) is taken from Shen et al. (1998), which is based on relativistic mean field theory. The equation of state for the electron and positron gas includes arbitrary relativistic pairs as well as arbitrary degeneracy. Electron and positron capture on nuclei, as well as on free nucleons, are included, along with the use of the up-to-date rates from Langanke & Martinez-Pinedo (2000). The capture is suppressed above the neutrino trapping density, taken to be $`3\times 10^{11}`$ g cm-3, since the neutrino transport process is not taken into account in this study. This simplification may be appropriate, since neutrinos at the early epoch of the core bounce do not appear to make a significant correction to the entropy compared to that obtained from shock heating as shown by Hillebrandt, Nomoto, & Wolff (1984). The delayed neutrino heating may not significantly modify the mass trajectories of the outgoing matter either, since the bulk of the ejecta are lifted to $`1000`$ km during the first few 100 ms from the onset of the core bounce (as shown below), which is far from the location of the neutrino sphere (a few tens of kilometers). It is evident, however, that an accurate treatment of neutrino transport will be needed to obtain accurate mass trajectories in future work.

Nuclear burning is implemented in a simplified manner. The composition of the O-Ne-Mg core is held fixed until the temperature in each zone reaches the onset of oxygen-burning, taken to be $`2\times 10^9`$ K (Nomoto, 1987), at which point the matter is assumed to be instantaneously in nuclear statistical equilibrium (NSE). The temperature is then calculated by including its nuclear energy release. It should be noted that we find a weak $`\alpha `$-rich freezeout in the subsequent post-processing nucleosynthesis calculations (§ 3), owing to the entropy, $`10N_Ak`$, in the ejecta. This means that the outgoing ejecta are not in perfect NSE, which is assumed in our hydrodynamic calculations. An improvement that takes the non-NSE matter properly into account will be needed to obtain more accurate trajectories.

We begin the hydrodynamical computations with this pre-supernova model, which has a density of $`4.4\times 10^{10}`$ g cm-3 and temperature of $`1.3\times 10^{10}`$ K at its center. The inner $`0.1M_{}`$ has already burned to NSE. As a result, the central $`Y_e`$ is rather low, 0.37, owing to electron capture. The core bounce is initiated when $`90`$ ms has passed from the start of the calculation. At this time the NSE core contains only $`1.0M_{}`$, which is significantly smaller than the cases of collapsing iron cores ($`1.3M_{}`$). The central density is $`2.2\times 10^{14}`$ g cm-3, significantly lower than that of Hillebrandt, Nomoto, & Wolff (1984, slightly above $`3\times 10^{14}`$ g cm-3), although the temperature ($`=2.1\times 10^{10}`$ K) and $`Y_e`$ (= 0.34), are similar. This difference is perhaps due to the use of a relatively stiff EOS in this study.

We find that a very weak explosion results, with an ejected mass of $`0.008M_{}`$ and an explosion energy of $`2\times 10^{49}`$ ergs (model Q0 in Table 1). The time variations of the radius, temperature, and density of each zone are displayed in Figure 1. The lowest $`Y_e`$ in the outgoing ejecta is 0.45, where no $`r`$-processing is expected given the entropy of $`10N_Ak`$. This is in contrast to the very energetic explosions, with ejected masses of $`0.2M_{}`$, explosion energies of $`2\times 10^{51}`$ ergs, and low $`Y_e`$ of $`0.2`$ obtained by Hillebrandt, Nomoto, & Wolff (1984). This might be a consequence of the lower gravitational energy release owing to the EOS applied in this study.

In order to examine the possible operation of the $`r`$-process in the explosion of this model, we artificially obtain explosions with typical energies of $`10^{51}`$ ergs by application of a multiplicative factor ($`f_{\mathrm{shock}}`$) to the shock-heating term in the energy equation (models Q3, Q5, and Q6 in Table 1). We take this simplified approach in this study, since the main difference between our result and that by Hillebrandt, Nomoto, & Wolff (1984) appear to be the lower central density in ours. If the inner core reached a higher density at the time of core bounce by applying, for example, a softer EOS, the matter would obtain higher shock-heating energy. This is clearly not a self-consistent approach, and a further study is needed to conclude whether such a progenitor star explodes or not, taking into account a more accurate treatment of neutrino transport, as well as with various sets of input physics (like EOSs). It should be emphasized, however, that our purpose in this paper is not to justify the prompt explosions of collapsing O-Ne-Mg cores, but to investigate the conditions necessary for the production of $`r`$-process nuclei from such an event if it occurs.

Table 1 lists the multiplicative factor applied to the shock-heating term ($`f_{\mathrm{shock}}`$), explosion energy ($`E_{\mathrm{exp}}`$), ejected mass ($`M_{\mathrm{ej}}`$), and minimum $`Y_e`$ in the ejecta obtained for each model. Energetic explosions with $`E_{\mathrm{exp}}>10^{51}`$ ergs are obtained for $`f_{\mathrm{shock}}1.5`$ (models Q5 and Q6), in which deeper neutronized zones are ejected by the prompt shock, as can be seen in Figure 2 (model Q6). This is in contrast to the weak explosions with $`E_{\mathrm{exp}}10^{50}`$ ergs (models Q0 and Q3), in which only the surface of the core blows off (Figure 1). Note that the remnant masses for models Q5 and Q6 are $`1.19M_{}`$ and $`0.94M_{}`$, respectively, which are significantly smaller than the typical neutron star mass of $`1.4M_{}`$. We consider it likely that a mass of $`1.4M_{}`$ is recovered by fallback of the once-ejected matter, as discussed in § 4.

In Figure 3 the electron fraction in the ejecta of each model is shown as a function of the ejected mass point, $`M_{\mathrm{ej}}`$. For models Q0 and Q3, $`Y_e`$ decreases steeply with $`M_{\mathrm{ej}}`$, since the duration of electron capturing is long, owing to the slowly expanding ejecta (Figure 1). For models Q5 and Q6, on the other hand, $`Y_e`$ decreases gradually with $`M_{\mathrm{ej}}`$, owing to the fast expansion of the outgoing ejecta. Nevertheless, the inner regions approach very low $`Y_e`$, 0.30 and 0.14 for models Q5 and Q6, respectively, owing to their rather high density ($`10^{11}`$ g cm-3) at the time of core bounce (Figure 2). Note that, for model Q6, $`Y_e`$ increases again for $`M_{\mathrm{ej}}>0.3M_{}`$. This is due to the fact that the positron capture on free neutrons overcomes the electron capture on free protons when the electron degeneracy becomes less effective in the high temperature matter. The innermost, slowly outgoing region suffers from this effect. This can be also seen in the results of Mayle & Wilson (1988), who obtained a neutrino-powered (not prompt) explosion with the same O-Ne-Mg core.

The minimum $`Y_e`$ ($`0.40`$) in Mayle & Wilson (1988) is, however, significantly higher than ours. This is mainly due to neutrino capture on free nucleons, since the matter is driven by intense neutrino flux in their simulation. $`Y_e`$ might increase further once nucleons begin to assemble into $`\alpha `$ particles and heavy nuclei (the “$`\alpha `$-effect”; Meyer, McLaughlin, & Fuller, 1998). In the case of the prompt explosions considered here, however, these neutrino effects may not alter $`Y_e`$ significantly. The reason is that the bulk of the ejecta are lifted to $`1000`$ km at the arrival of the delayed neutrinos (a few 100 ms from the onset of the core bounce), at which the capture timescale of neutrinos on free nucleons is no less than a few seconds (see, e.g., eq. (1) in Qian et al., 1997). The dynamical timescales of the outgoing mass shells in model Q6 at $`T_9=1,3,5`$, and 7 (after the core bounce), defined by $`\tau _{\mathrm{dyn}}=|\rho /(d\rho /dt)|`$, are shown in Figure 4. As can be seen, the dynamical timescale prior to the $`r`$-process phase ($`T_93`$) is significantly smaller than that of neutrino interaction.

The trend of the $`Y_eM_{\mathrm{ej}}`$ relation up to $`M_{\mathrm{ej}}0.2M_{}`$ is similar in models Q5 and Q6, although it is inverted at $`M_{\mathrm{ej}}0.14M_{}`$, owing to the slightly different contribution of the positron and electron capture on free nucleons (Figure 3). Hence, the $`Y_eM_{\mathrm{ej}}`$ relation between the surface and the innermost layer of the ejecta is expected to be similar to that of model Q6, as long as the explosion is sufficiently energetic ($`10^{51}`$ ergs). In the subsequent sections, therefore, we focus only on model Q6, which is taken to be representative of cases where $`r`$-process nucleosynthesis occurs. The ejected mass, $`M_{\mathrm{ej}}`$, is thus taken to be a free parameter, instead of simulating many other models by changing $`f_{\mathrm{shock}}`$. Note that the results by Hillebrandt, Nomoto, & Wolff (1984) and Sumiyoshi et al. (2001) are very similar to the cases with $`M_{\mathrm{ej}}0.20.3M_{}`$ in model Q6.

## 3 THE $`r`$-PROCESS

The yields of $`r`$-process nucleosynthesis species, adopting the model described in § 2 for the physical conditions, are obtained by application of an extensive nuclear reaction network code. The network consists of $`3600`$ species, all the way from single neutrons and protons up to the fermium isotopes ($`Z=100`$). We include all relevant reactions, i.e., $`(n,\gamma )`$, $`(p,\gamma )`$, $`(\alpha ,\gamma )`$, $`(p,n)`$, $`(\alpha ,p)`$, $`(\alpha ,n)`$, and their inverses. Reaction rates are taken from Thielemann (1995, private communication) for nuclei with $`Z46`$ and from Cowan, Thielemann, & Truran (1991) for those with $`Z47`$. The latter used the mass formula by Hilf et al. (1976). The three-body reaction $`\alpha (\alpha n,\gamma )^9`$Be, which is of special importance as the bottleneck reaction to heavier nuclei, is taken from the recent experimental data of Utsunomiya et al. (2001). The weak interactions, such as $`\beta `$-decay, $`\beta `$-delayed neutron emission (up to three neutrons), and electron capture are also included, although the latter is found to be unimportant.

The $`\alpha `$-decay chains and spontaneous fission processes are taken into account only after the freezeout of all other reactions, as in Wanajo et al. (2002). For the latter, all nuclei with $`A256`$ are assumed to decay by spontaneous fission only. The few known nuclei undergoing spontaneous fission for $`A<256`$ are also included, along with their branching ratios. Neutron-induced and $`\beta `$-delayed fissions, as well as the contribution of fission fragments to the lighter nuclei, are neglected. Obviously, these treatments of the fission reactions are oversimplified. Nevertheless, this may be acceptable, at least to first order. We leave more accurate treatment of these matters to future work.

Each calculation is started at $`T_9=9`$ (where $`T_9T/10^9`$ K). The initial composition is taken to be that of NSE with the density and electron fraction at $`T_9=9`$, and consists mostly of free nucleons and alpha particles.

For reasons outlined in § 2, we examine the nucleosynthesis in model Q6 only, in which robust $`r`$-processing is expected. The resulted abundances in several representative Lagrangian mass shells are depicted in Figure 5. The initial electron fractions are 0.29, 0.20, 0.18, 0.16, 0.14, and 0.20 for the corresponding zone numbers 83, 92, 95, 98, 105, and 132, respectively. As can be seen, a robust $`r`$-processing is possible only for $`Y_e0.20`$. In particular, a substantial amount of thorium and uranium are produced only when $`Y_e`$ is less than 0.18.

The mass-integrated abundances from the surface (zone 1) to the zones 83, 92, 95, 98, 105, and 132 are compared with the solar $`r`$-process abundances (Käppeler et al., 1989) in Figure 6 (models Q6a-f in Table 2). The latter is scaled to match the height of the first ($`A=80`$) and third ($`A=195`$) peaks of the abundances in models Q6a-b and Q6c-f, respectively. The ejecta masses of these models are listed in Table 2. The nucleosynthesis result in model Q5 (not presented here) is expected to be similar to that of model Q6a, because of the resemblance of the $`Y_eM_{\mathrm{ej}}`$ profiles between these models (Figure 3). As can be seen in Figure 6, a solar $`r`$-process pattern for $`A130`$ is naturally reproduced in models Q6c-f, while models Q6a-b fail to reproduce the third abundance peak. This implies that the region with $`Y_e<0.20`$ must be ejected to account for production of the third $`r`$-process peak. Furthermore, to account for the solar level of thorium ($`A=232`$) and uranium ($`A=235,238`$) production, the region with rather low $`Y_e`$ ($`<0.18`$) must be ejected. Note that large deficiencies of nuclei at $`A115`$ may be supplied if the slower rates of neutron capture in this region are adopted, as demonstrated by Goriely (1997).

We find that, for models Q6c-e, the lighter $`r`$-process nuclei with $`A<130`$ are somewhat deficient compared to the solar $`r`$-process pattern (Figure 6c-e). This trend can be also seen in the observational abundance patterns of the highly $`r`$-process-enhanced, extremely metal- poor stars CS 22892-052 (Sneden et al., 2003) and CS 31082-001 (Hill et al., 2002). In model Q6f, the deficiency is outstanding because of large ejection of the low $`Y_e`$ matter (Figure 3). This is in contrast to the previous results obtained for the neutrino wind scenario, which significantly overproduce the nuclei with $`A90`$ (Woosley et al., 1994; Wanajo et al., 2001). The nuclei with $`A<130`$ can be supplied by slightly less energetic explosions, like models Q6a-b (Figures 5a-b). It is also possible to consider that these lighter $`r`$-process nuclei originate from “neutrino winds” in more massive supernovae ($`>10M_{}`$). The nuclei with $`A<130`$ can be produced naturally in neutrino winds with a reasonable compactness of the proto-neutron star, e.g., $`1.4M_{}`$ and 10 km (Wanajo et al., 2001).

Figure 6 implies that the production of thorium and uranium differs from model to model, even though the abundance pattern seems to be universal between the second and third $`r`$-process peaks, as seen in models Q6c-f. This is in agreement with recent observational results suggesting that the ratio Th/Eu may exhibit a star-to-star scatter, while the abundance pattern between the second and third peaks is in good agreement with the solar $`r`$-process pattern (Honda et al., in preparation). Thus, the use of Th/Eu as a cosmochronometer should be regarded with caution, at least until the possible variations can be better quantified; U/Th might be a far more reliable chronometer, as discussed further in § 5.

It is interesting to note that the nucleosynthesis results obtained for models Q6c-e are in good agreement with that of the prompt explosion of a $`11M_{}`$ star with an iron core by Sumiyoshi et al. (2001). Our results seem, however, in better agreement with the solar $`r`$-process pattern, in particular near the rare-earth peak ($`A160`$) and the third peak ($`A195`$). This is a consequence of the slower ejection of the innermost region in our results owing to the reduction of pressure by electron capture, which is not taken into account in Sumiyoshi et al. (2001). In particular, the dynamical timescale of outgoing matter becomes significantly long during the epoch of the $`r`$-process ($`T_93`$ to 1), as can be seen in Figure 4. The extremely low $`Y_e`$ ($`<0.2`$) drives matter near the neutron-drip line, where the abundance pattern deviates from the solar one. In our models, however, the solar $`r`$-process pattern is recovered by the “freezeout effect” after the epoch of $`r`$-processing, as discussed in detail by Wanajo et al. (2002). This arises because the matter stays at higher temperatures for a longer time, hence the quasi-equilibrium between neutron emission by photodisintegration and subsequent neutron-capture processes continues to operate (Surman & Engel, 2001).

It is also interesting to note that the mass of the $`r`$-processed ejecta in our results ($`0.05M_{}`$) is more than one order of magnitude larger than that obtained by Sumiyoshi et al. (2001, $`0.003M_{}`$). This is also a consequence of the longer dynamical timescale of outgoing ejecta in our simulation. The $`\alpha `$ particles are mostly consumed prior to the $`r`$-process phase in our calculations, whereas they dominate the final nucleosynthesis yields in Sumiyoshi et al. (2001). These results strongly suggest that an accurate treatment of electron capture, as well as of neutrino transport, will be crucial to derive an accurate prediction of the $`r`$-process abundances from prompt explosions of O-Ne-Mg cores.

## 4 CONTRIBUTION TO CHEMICAL EVOLUTION OF THE GALAXY

We have demonstrated in the previous section that, if stars of $`810M_{}`$ do explode energetically, the solar $`r`$-process pattern can be reproduced quite naturally. It is of importance, however, to see if the nucleosynthetic material contributed to the Galaxy from such stars is consistent with the currently observed elemental abundances in the solar neighborhood.

Figure 7 shows the “production factors” per supernova for model Q6e, defined for each nuclide as the final mass fraction, $`X_{\mathrm{ej}}`$, diluted by the total ejected mass ($`7.9M_{}`$) from the $`9M_{}`$ star, divided by its solar abundance $`X_{}`$ (Anders & Grevesse, 1989). The solid lines connect isotopes of a given element (after nuclear decay). The dotted horizontal lines indicate a “normalization band” between the largest production factor (129Xe) and that by a factor of ten less than that, along with a median value (dashed line). This band is taken to be representative of the uncertainty in the nuclear data for very neutron-rich nuclei (Woosley et al., 1994). As can be seen in Figure 7, a majority of the nuclei with $`A>100`$ fall within the normalization band (except for the large deficiencies near $`A=120`$), which is regarded to be the dominant species produced. The deficiencies near $`A=88,138`$, and 208 may be supplied by $`s`$-process contributions from other sources. Although such a supernova might be considered a likely production site of $`r`$-process nuclei, its contribution to nuclei with $`A<80`$, including the $`\alpha `$ elements and iron-peak species, is negligible.

One of the essential questions raised by previous works is that, if prompt supernova explosions are one of the major sites of $`r`$-process nuclei, would in fact the $`r`$-process nuclei be significantly overproduced (Hillebrandt, Takahashi, & Kodama, 1976). As far as the explosion is purely hydrodynamical, a highly neutronized deeper region must be ejected in order for a successful $`r`$-process to result. It seems inevitable, therefore, that one must avoid an ejection of large amounts of $`r`$-process matter, at least when assuming spherical symmetry. Our result shows that more than $`0.05M_{}`$ of the $`r`$-process matter ($`A120`$) is ejected per event, which reproduces the solar $`r`$-process pattern (models Q6c-f in Table 2). This is about three orders of magnitude larger than the $`5.8\times 10^5M_{}`$ in the neutrino-heated supernova ejecta from a $`20M_{}`$ star obtained by Woosley et al. (1994). This value might would be reduced to some extent by $`\alpha `$-rich freezeout in faster outgoing mass shells than ours as observed in Sumiyoshi et al. (2001). However, good agreement with the solar $`r`$-process pattern would not be achieved in the ejecta expanding too fast, as discussed in § 3.

As discussed by Woosley et al. (1994), production factors must be on the order of $`10`$ for the case that all supernovae contribute equally to $`r`$-process production in the Galaxy. Stars of $`810M_{}`$ would account for $`30\%`$ of all supernovae (if they explode). The ejected masses ($`79M_{}`$) are smaller than a factor of two to three than those from more massive supernovae. Thus, production factors of $`100`$ are allowed in the case that all $`810M_{}`$ stars contribute equally to the Galactic $`r`$-process material. The production factors in our results are, however, about three orders of magnitude higher than this (Figure 7).

It might be argued that this type of event is extremely rare, accounting for only $`0.010.1\%`$ of all core-collapse supernovae. However, observations of extremely metal-poor stars ($`[\mathrm{Fe}/\mathrm{H}]3`$) in the Galactic halo show that at least two, CS 22892-052 and CS 31082-001, out of about a hundred studied at high resolution, imply contributions from highly $`r`$-process-enhanced supernova ejecta (Sneden et al., 2000; Hill et al., 2002). Moreover, such an extremely rare event would result in a much larger dispersion of $`r`$-process elements relative to iron than is currently observed amongst extremely metal-poor stars. Ishimaru & Wanajo (1999) demonstrated that the observed star-to-star dispersion of $`[\mathrm{Eu}/\mathrm{Fe}]`$ over a range $`1`$ to 2 dex, was reproduced by their chemical evolution model if Eu originated from stars of $`810M_{}`$. Recent abundance measurements of Eu in a few extremely metal-poor stars with $`[\mathrm{Fe}/\mathrm{H}]3`$ by SUBARU/HDS further supports their result (Ishimaru et al., in preparation). The requisite mass of Eu in their model is $`10^6M_{}`$ per event. The ejected mass of Eu in our result is more than two orders of magnitude larger (Table 2).

In order to resolve this conflict, we propose that the “mixing-fallback” mechanism operates in this kind of supernova. The peculiar abundance patterns of some extremely iron-deficient stars, including HE 0107-5240 with $`[\mathrm{Fe}/\mathrm{H}]=5.3`$ (Christlieb et al., 2002), is explained successfully with this mechanism, as proposed by Umeda & Nomoto (2002, 2003). If a substantial amount of the hydrogen and helium envelope remains at the onset of the explosion, the outgoing ejecta may undergo large-scale mixing by Rayleigh-Taylor instabilities, which is believed to have happened in SN1987A (Hachisu et al., 1990; Herant & Woosley, 1994; Kifonidis et al., 2003). Thus a tiny amount, say, $`1\%`$, of the $`r`$-process material is mixed into the outer layers and then ejected, but most of the core material may fall back onto the proto-neutron star via the reverse shock arising from the hydrogen-helium layer interface. In this case, the typical mass of the proto-neutron star ($`1.4M_{}`$) is recovered. An asymmetric explosion mechanism, such as that which might arise from rotating cores, may have a similar effect as the ejection of deep-interior material in a small amount (Yamada & Sato, 1994; Fryer & Heger, 2000). This “mixing-fallback” scenario must be further tested by detailed multidimensional-hydrodynamic studies. However, it may provide us with a new paradigm for the nature of supernova nucleosynthesis.

This type of event may be characterized as a “faint” supernova, owing to the weakened explosion energy as well as the reduced amount of 56Ni by “mixing-fallback”. In addition, its ejecta consists mostly of hydrogen and helium with a higher ratio of He/H than that in the solar system. This event can be easily distinguished from the core-collapse supernovae of iron cores resulted from stars of $`>10M_{}`$, which are characterized by abundant $`\alpha `$-elements. It is interesting to note that such an explosion bears a close resemblance to the Crab supernova, SN1054 (Nomoto et al., 1982). There is the possibility, therefore, that direct evidence will be obtained from the Crab nebula (or other similar supernova remnants) through the detection of $`r`$-process elements by future observations (Wallerstein et al., 1995). Detection of the $`\gamma `$-ray lines from decays of radioactive isotopes produced by the $`r`$-process could provide additional direct evidence for this scenario (Qian, Vogel, & Wasserburg, 1998).

It is of special importance to confirm, from spectroscopic studies of extremely metal-poor stars in the Galactic halo, that the $`r`$-process elements are not associated with the production of $`\alpha `$\- and iron-peak elements (Wasserburg & Qian, 2000; Qian & Wasserburg, 2001, 2002). Qian & Wasserburg (2003) have suggested that the $`r`$-process enrichment in extremely metal-poor stars HD 115444, HD 122563, and CS 31082-001 is independent of the production of the elements from Na to Zn (including $`\alpha `$\- and iron-peak elements). The abundances of Na-Zn among these stars are mostly the same, while the level of $`r`$-process enhancement differs from star to star. These authors take this as evidence that the heavy $`r`$-process nuclei originated from AICs or Type II supernovae from $`810M_{}`$ stars, although they prefer AICs for the $`r`$-process site. Note that the nucleosynthetic outcome of an AIC event may be similar to the collapsing O-Ne-Mg core resulting from a single $`810M_{}`$ star. Hence, it is possible that AICs also undergo prompt explosions. The frequency of the AIC events in the Galaxy is expected to be very small, $`10^5`$ yr-1 (Bailyn & Grindlay, 1990). Assuming that the frequency of the core-collapse supernovae in the Galaxy is $`10^2`$ yr-1 (Cappellaro et al., 1997), the production factor per AIC event is $`10^4`$. This is in reasonable agreement with our result of $`10^5`$, as can be seen in Figure 7. In the case of AICs, however, there is no chance to undergo “mixing-fallback”, owing to the absence of the outer envelope. Thus, a much larger dispersion of \[Eu/Fe\] than that observed in the extremely metal-poor stars seems difficult to avoid.

We consider, therefore, that the heavy $`r`$-process nuclei originate from $`810M_{}`$ stars discussed in this study, being independent of whether they are single stars or are in binary systems. The lighter $`r`$-process elements ($`Z<56`$) observed in extremely metal-poor stars may represent simply the interstellar medium from which these stars were formed, which originated from, perhaps, “neutrino winds” in supernovae from stars of $`>10M_{}`$. The observed abundances of Eu relative to iron in extremely metal-poor stars are in fact well reproduced with this assumption in the chemical evolution model by Ishimaru & Wanajo (1999). Future spectroscopic studies of extremely metal-poor stars may be able to distinguish the scenario suggested in this study from AICs, for example, through observations of $`s`$-process abundances originating from mass transfer in binary systems.

## 5 U-Th COSMOCHRONOLOGY

The recently discovered $`r`$-process-enhanced, extremely metal-poor star, CS 31082-001 ($`[\mathrm{Fe}/\mathrm{H}]=2.9`$) has provided a new, potentially quite powerful cosmochronometer, uranium (Cayrel et al., 2001; Hill et al., 2002). Wanajo et al. (2002) have determined the age of this star (more precisely, the time that has passed since the production event of the $`r`$-process species incorporated into this star) using the U-Th chronometer pair to be $`14.1\pm 2.5`$ Gyr, based on the neutrino wind scenario. This is regarded as a hard lower limit of the age of the universe, and is in good agreement with that derived by an site-independent approach (Goriely & Arnould, 2001; Schatz et al., 2002). It is of interest to examine if the same age of this star is also obtained, based on the prompt explosion model of a $`9M_{}`$ star described in the previous sections. Note that the same mass formula (Hilf et al., 1976) as in Wanajo et al. (2002) is used in this study for comparison.

In Figure 8 the available spectroscopic abundance data for CS 31082-001 (dots) are compared with the nucleosynthesis result of model Q6e discussed in § 3 (thick line) and with the solar $`r`$-process pattern (thin line), scaled at Eu ($`Z=63`$). The data for the neutron-capture elements in this star are taken from Hill et al. (2002). An overall agreement of our result with the spectroscopic data up to lead ($`Z=82`$) can be seen, although our result appears somewhat deficient for the lighter elements. This is in contrast to the results in neutrino winds (Wanajo et al., 2002), in which the lighter elements are significantly overproduced. Thus, this model might be a reasonable one for the $`r`$-process events that produce large amounts of thorium and uranium, and whose products now appear in the atmosphere of CS 31082-001.

Figure 9 shows the mass-integrated abundance ratios of Th/Eu (open squares) and U/Th (open circles) from the surface to the mass point $`M_{\mathrm{ej}}`$, as well as the inferred ages of CS 31082-001, $`t^{}(\mathrm{Th}/\mathrm{Eu})`$ (filled squares) and $`t^{}(\mathrm{U}/\mathrm{Th})`$ (filled circles). As discussed by Hill et al. (2002), the age of CS 31082-001 can be inferred by application of the following relations:

$`t^{}(\mathrm{Th}/\mathrm{Eu})`$ $`=`$ $`46.67[\mathrm{log}(\mathrm{Th}/\mathrm{Eu})_0\mathrm{log}(\mathrm{Th}/\mathrm{Eu})_{\mathrm{now}}]\mathrm{Gyr}`$ (1)
$`t^{}(\mathrm{U}/\mathrm{Th})`$ $`=`$ $`21.76[\mathrm{log}(\mathrm{U}/\mathrm{Th})_0\mathrm{log}(\mathrm{U}/\mathrm{Th})_{\mathrm{now}}]\mathrm{Gyr}`$ (2)
where the half lives of 232Th (14.05 Gyr) and 238U (4.468 Gyr), and the subscripts “0” and “now” denote the initial and current values derived by theory and observation, respectively. In principle, $`t^{}(\mathrm{U}/\mathrm{Th})`$ may serve as a more precise chronometer than $`t^{}(\mathrm{Th}/\mathrm{Eu})`$, owing to the smaller coefficients in front of equations (2). Moreover, the ratio U/Th is less dependent on the model parameter $`M_{\mathrm{ej}}`$, since these species are separated by only two units in atomic number. Note that 235U is assumed to have $`\alpha `$-decayed away because of its relatively short half life (0.704 Gyr).

As can be seen in Figure 9, U/Th approaches a constant value ($`=0.51`$) for $`M_{\mathrm{ej}}0.26M_{}`$ ($`Y_e0.17`$), while Th/Eu varies widely. As a result, the age of CS 31082-001 determined by U/Th results in a constant value, $`t^{}(\mathrm{U}/\mathrm{Th})=14.1`$ Gyr, for $`M_{\mathrm{ej}}0.26M_{}`$. The age $`t^{}(\mathrm{Th}/\mathrm{Eu})`$ is sensitive to the parameter $`M_{\mathrm{ej}}`$, ranging from a negative age to 23.8 Gyr, which illustrates that caution must be used in the application of this chronometer pair. It is useful, however, to take $`t^{}(\mathrm{Th}/\mathrm{Eu})`$ as a constraint on the model parameter $`M_{\mathrm{ej}}`$, although the result might be changed if one includes an accurate neutrino transport and other input physics. It is found that the models with $`M_{\mathrm{ej}}=0.30M_{}`$ and $`0.37M_{}`$ (zone numbers 105 and 118, respectively) give the same ages between $`t^{}(\mathrm{Th}/\mathrm{Eu})`$ and $`t^{}(\mathrm{U}/\mathrm{Th})`$. The former corresponds to model Q6e, which would provide a consistent scenario for the origin of the $`r`$-process elements in CS 31082-001, as can be seen in Figure 8. It should be noted that the fission fragments in model Q6e account for 23% of the mass contained in $`A120`$ nuclei (the last column in Table 2), whose contribution to the abundance pattern is neglected in this study. Obviously, more accurate treatment of fission reactions is needed. Nevertheless, the age $`t^{}(\mathrm{U}/\mathrm{Th})`$ would not be altered significantly by this improvement, since the ratio U/Th is at the saturated value ($`=0.51`$) at $`M_{\mathrm{ej}}=0.26M_{}`$, where the mass fraction of fission fragments is only 2%.

It is noteworthy that the inferred age by U/Th in this study is the same as the result in Wanajo et al. (2002), which is based on a different astrophysical scenario (neutrino winds), but uses the same nuclear mass formula (Hilf et al., 1976). Our result confirms the robustness of the age determination using the U-Th chronometer pair, which is mostly independent of the astrophysical conditions considered. Rather than the $`r`$-process site, the nuclear mass formulae adopted, as well as the treatment of fission are crucial as far as the U-Th pair is concerned (Seeger & Schramm, 1970; Goriely & Arnould, 2001; Schatz et al., 2002).

## 6 SUMMARY AND CONCLUSIONS

We have examined the $`r`$-process nucleosynthesis obtained in the prompt explosion arising from the collapse of a $`9M_{}`$ star with an O-Ne-Mg core. The core collapse and subsequent core bounce were simulated with a one-dimensional, implicit, Lagrangian hydrodynamic code with Newtonian gravity. Neutrino transport was neglected for simplicity. We obtained a very weak explosion (model Q0) with an explosion energy of $`2\times 10^{49}`$ ergs, and an ejected mass of $`0.008M_{}`$. No $`r`$-processing occurred in this model, because of the high electron fraction ($`0.45`$) with low entropy ($`10N_Ak`$).

We further simulated energetic explosions by an artificial enhancement of the shock-heating energy, which might be expected from calculations with other sets of input physics, as well as with other pre-supernova models. This resulted in an explosion energy of $`10^{51}`$ ergs and an ejected mass of $`0.2M_{}`$. Highly neutronized matter ($`Y_e0.14`$) was ejected, which led to strong $`r`$-processing (model Q6). Material arising from $`r`$-process nucleosynthesis was calculated with a nuclear reaction network code containing $`3600`$ isotopes with all relevant reactions. The result was in good agreement with the solar $`r`$-process pattern, in particular for nuclei with $`A>130`$. Some of lighter $`r`$-process nuclei ($`A<130`$) were deficient, which is consistent with the abundance patterns of the highly $`r`$-process enhanced, extremely metal-poor stars CS 22892-052 and CS 31082-001. This implies that the lighter $`r`$-process nuclei may originate from another site, which we suggest might be associated with the “neutrino wind” in core-collapsing supernovae of iron cores resulted from more massive stars ($`>10M_{}`$).

The large ejection of $`r`$-process material ($`0.05M_{}`$ per event) conflicts with the level of dispersion of $`r`$-process elements relative to iron observed in extremely metal-poor stars. We suggest, therefore, that only a small fraction ($`1\%`$) of the $`r`$-processed material is ejected, while the bulk of such material falls back onto the proto-neutron star by the “mixing-fallback” mechanism.

The age of the highly $`r`$-process enhanced, extremely metal-poor star, CS 31082-001 was derived by application of the U-Th chronometer pair, and can be regarded as a hard lower limit on the age of the universe. The age obtained is $`14.1\pm 2.4`$ Gyr (the quoted error only includes that arising from the observations), the same as that based on a different astrophysical site, the neutrino-wind scenario (Wanajo et al., 2002), using the same nuclear mass formula. This confirms that the age determined by the U-Th pair is mostly independent of the astrophysical environment considered. The dependence of age dating on different nuclear mass formulae, based on the prompt explosion scenario presented in this paper, will be reported in future work.

It is obvious that more studies, including an accurate treatment of neutrino transport and other input physics, as well as multi-dimensional simulations of prompt supernova explosions, are needed to derive the final conclusion on the $`r`$-process scenario presented in this study. Nevertheless, this scenario is attractive as a promising site of the $`r`$-process, since the solar $`r`$-process pattern can be naturally reproduced without the problematic overproduction of $`A90`$ that appeared in the neutrino wind scenario. This type of event is characterized with the absence of $`\alpha `$\- and iron-peak elements, which can be easily distinguished from that of the core-collapsing iron core resulted from a more massive star ($`>10M_{}`$). Future spectroscopic studies of extremely metal-poor stars in the Galactic halo will reveal if the collapsing O-Ne-Mg cores of $`810M_{}`$ stars are a viable site for the production of $`r`$-process nuclei. Improved observational determination of the U/Th ratio in CS 31082-001, as is presently being pursued, as well as a measured abundance of Pb in this star (as is being obtained with the Hubble Space Telescope), and the identification of a greater number of highly r-process-enhanced, metal-poor stars, also underway, will surely deepen our understanding of the relevant processes involved.

We would like to acknowledge K. Sumiyoshi for providing an EOS table of nuclear matter applied in this study and for helpful discussions. We also thank G. Martinez-Pinedo for providing a data table of weak rates. We also acknowledge the contributions of an anonymous referee, which led to clarification of a number of points in our original presentation. This work was supported by a Grant-in-Aid for Scientific Research (13640245, 13740129, 14047206, 14540223) from the Ministry of Education, Culture, Sports, Science, and Technology of Japan. T.C.B. acknowledges partial support from grants AST 00-98508 and AST 00-98549 awarded by the U.S. National Science Foundation.