D. Argast 1,2 - M. Samland 2 - F.-K. Thielemann 1 - Y.-Z. Qian 3
1 - Institut für Physik der Universität Basel,
Klingelbergstrasse 82, 4056 Basel, Switzerland
2 - Astronomisches Institut der Universität Basel,
Venusstrasse 7, 4102 Binningen, Switzerland
3 - School of Physics and Astronomy, University of Minnesota,
Minneapolis, Minnesota 55455, USA
Received 2 September 2003 / Accepted 6 November 2003
Abstract
The astrophysical nature of r-process sites is a long-standing
mystery and many probable sources have been suggested, among them lower-mass
core-collapse supernovae (in the range 8-
),
higher-mass core-collapse supernovae (with masses
20
)
and neutron star mergers. In this work, we present a
detailed inhomogeneous chemical evolution study that considers for the first
time neutron star mergers as major r-process sources, and compare this
scenario to the ones in which core-collapse supernovae act as dominant
r-process sites. We conclude that, due to the lack of reliable iron and
r-process yields as a function of progenitor mass, it is not possible at
present to distinguish between the lower-mass and higher-mass supernovae
scenarios within the framework of inhomogeneous chemical evolution. However,
neutron-star mergers seem to be ruled out as the dominant r-process
source, since their low rates of occurrence would lead to r-process enrichment
that is not consistent with observations at very low metallicities.
Additionally, the considerable injection of r-process material by a single
neutron-star merger leads to a scatter in [r-process/Fe] ratios at later times
which is much too large compared to observations.
Key words: nuclear reactions, nucleosynthesis, abundances - stars: abundances - ISM: abundances - Galaxy: abundances - Galaxy: evolution - Galaxy: halo
The principal production mechanisms of heavy elements beyond the iron peak
have been known since the classical paper by Burbidge et al. (1957,
B2FH): a major fraction of these elements is formed by slow or rapid
capture of neutrons on seed nuclei. These nucleosynthesis channels are denoted
as the s-process, if the capture of neutrons occurs on timescales slower than
the -decay timescales of the newly formed nuclei, or as the r-process,
if neutron capture operates on timescales much shorter than the relevant
-decay timescales. Although the physical requirements for the
occurrence of r-process nucleosynthesis are well understood (see e.g. Pfeiffer
et al. 2001), the astrophysical nature of the dominant r-process site
is still unknown.
Ultra metal-poor halo stars were formed when the s-process had little time to
make any significant contribution to the interstellar medium (ISM). The
abundance pattern of neutron capture elements heavier than Ba ()
in
these stars matches the scaled solar system r-process abundances remarkably
well. This suggests that, contrary to s-process nucleosynthesis, the synthesis
of r-process elements started early in galactic evolution and that the
r-process is robust for elements heavier than Ba, i.e. originates from a
single astrophysical site or at least occurs under well-defined physical
conditions (e.g. Sneden et al. 2000a; Westin et al. 2000; Cowan
et al. 2002; but see also Hill et al. 2002). On the other hand,
lighter neutron capture elements (
30 < Z < 56) in ultra metal-poor halo
stars show significant deviations from the scaled solar system r-process
abundance curve, indicating the possible existence of a second r-process
source. The diversity of r-process sources was first proposed based on
meteorite data, which seem to require two r-process sources operating on
different timescales, namely
107 years for the production of
nuclei and
108 years for Z < 56 nuclei (Wasserburg et al. 1996). In principle, these two r-process sites could either be two
types (frequent and less common) of core-collapse supernova (SN II) events
(Wasserburg & Qian 2000; Qian 2000, 2001), or a mixture
of SNe and neutron star mergers (Rosswog et al. 1999, 2000). On
the other hand, Cameron (2001) showed that it is possible to reproduce
the abundance distribution of both light and heavy neutron capture elements in
different regions of a single source, namely the accretion-disk and jet
forming near the proto-neutron star in SN II events. However, the evidence
that the r-process beyond Ba (at least in the range
)
is
robust is very convincing and we will use this property in the following to
gain some insight into the enrichment of the Galaxy with r-process elements
beyond Ba.
A number of possible astrophysical sites responsible for the robust r-process
were put forth in the past, most of them linked to the violent death of
massive stars in supernova events. Wheeler et al. (1998) suggested
that, during the prompt explosion of a massive star in the range 8-
,
physical conditions in the innermost mass layers are
sufficiently extreme for the r-process to work and that the amount of ejected
r-process matter may be consistent with observed Galactic r-process
abundances. Detailed studies of artificially induced prompt explosions by
Sumiyoshi et al. (2001) and Wanajo et al. (2003) indeed show
that r-process nucleosynthesis in lower-mass SNe II may match the solar
r-process signature - if such prompt explosions really occur. On the other
hand, the mass of the ejected r-process material is too large to be consistent
with chemical evolution studies. A possible solution to this problem may be
the fall-back of r-process material onto the stellar remnant (Wanajo et al. 2003). However, there are major objections to the prompt explosion
mechanism from detailed studies by Bruenn (1989a,b). Optimistically speaking, if prompt explosions of
massive stars may occur in reality, lower-mass progenitors are the only ones
that show any hope of doing so (e.g. Wheeler et al. 1998; Sumiyoshi et al. 2001; but see also Liebendörfer et al. 2001).
On the other hand, neutrino-driven winds from nascent neutron-stars were also
proposed as a promising site for r-process nucleosynthesis (e.g. Woosley &
Hoffman 1992; Takahashi et al. 1994; Woosley et al. 1994;
Qian & Woosley 1996; Thompson et al. 2001; Wanajo et al. 2001; Terasawa et al. 2002). Following the delayed
explosion of a very massive star (
), neutrinos
diffusing out of the contracting proto-neutron star heat and ablate material
from its surface and a neutrino-driven wind develops. Such winds may exhibit
sufficiently high entropies and/or sufficiently short dynamical timescales for
the r-process to occur. However, r-process yields consistent with observed
r-process abundances in stars may be obtained in this scenario only for
extreme assumptions such as massive neutron stars of
or more, which makes this not a very likely scenario. Furthermore, it seems
questionable whether the high entropies required to reproduce the solar system
r-process signature (e.g. Meyer & Brown 1997; Freiburghaus et al. 1999a) can be provided by Type II supernovae (Witti et al. 1994; Qian & Woosley 1996; Thompson et al. 2001).
The r-process scenarios associated with core-collapse supernovae discussed
above are affected by considerable theoretical uncertainties as the physics of
these events is not really understood (e.g. Rampp & Janka 2000;
Liebendörfer et al. 2001). Recently, Freiburghaus et al.
(1999b) presented for the first time detailed r-process calculations
for neutron star mergers (NSM). Coalescing neutron stars potentially can
provide in a natural way the large neutron fluxes required for the build-up of
heavy elements through rapid neutron capture. This scenario was subsequently
explored and refined by Rosswog et al. (1999; 2000). Their
calculations showed that a few times 10-3-
of r-process matter might be ejected in a merger event. This amount is
significantly larger than the typical 10-5-
of r-process material thought to be ejected in each core-collapse SN event
(Cowan et al. 1991; Woosley et al. 1994; Qian 2000;
Wanajo et al. 2001). Since the rate of NSM in the Galaxy is
significantly lower than that of Type II SNe (e.g. Tamman et al. 1994;
Belczynski et al. 2002), either of these two sources may account for
the total amount of r-process matter in the Galaxy (Qian & Woosley
1996; Rosswog et al. 1999; Rosswog & Davies 2002;
Thielemann et al. 2002).
However, Qian (2000) argued that if NSM were a major r-process source,
the low coalescence rate of binary compact objects would prevent any
correlation between abundances of r-process elements and iron, which is
clearly in disagreement with observed r-process abundances in stars more
metal-rich than [Fe/H
.
In addition, the large amount of r-process
ejecta from NSM would lead to r-process abundances in metal-poor halo stars
that would be too large compared to observations. These arguments deserve
closer examination, especially since few detailed chemical evolution
calculations exist to date that assume NSM to be the major r-process source.
Chemical evolution studies that treat SN with either lower-mass (8-
)
or higher-mass (
)
progenitors
as the major r-process source were carried out by Ishimaru & Wanajo
(1999), Travaglio et al. (1999) and Tsujimoto et al.
(2000). Unfortunately, each of these studies favoured different
r-process sites: lower-mass SNe (Travaglio et al. 1999), higher-mass
SNe (Tsujimoto et al. 2000), or both (Ishimaru & Wanajo 1999)
are acceptable from the viewpoint of chemical evolution. Recently, De Donder
& Vanbeveren (2003) studied the impact of NSM as major r-process
sources on Galactic chemical evolution under the assumption that the ISM is
well mixed at all times. In this scenario, enrichment of r-process elements by
NSM seem to be consistent with observations. The basis of the present paper is
that we do not assume such an instantaneous mixing.
The purpose of this work is twofold. First, neutron star mergers are included as major r-process sites in detailed inhomogeneous chemical evolution calculations for the first time. Second, we compare the r-process enrichment of the ISM under the assumption that neutron star mergers, lower-mass SNe II or higher-mass SNe II are the dominant r-process sites in the framework of inhomogeneous chemical evolution. In Sect. 2, the inhomogeneous chemical evolution model is presented. The impact of SNe II and NSM as r-process sites on Galactic chemical evolution is studied in Sect. 3 and conclusions are presented in Sect. 4.
The large scatter seen in the abundances of neutron capture elements in metal-poor halo stars is generally attributed to local chemical inhomogeneities of the interstellar medium during the earliest stages of Galaxy formation and enrichment (e.g. McWilliam et al. 1995a,b; Ryan et al. 1996; Burris et al. 2000; Mishenina & Kovtyukh 2001). Several authors developed chemical evolution models (Ishimaru & Wanajo 1999; Raiteri et al. 1999; Tsujimoto et al. 1999, 2000; Argast et al. 2000; Travaglio et al. 2001; Fields et al. 2002) that investigate the impact of local chemical inhomogeneities on the enrichment of the ISM and the scatter in element abundances. Although the methods employed by these authors differ significantly, the model results are in qualitative agreement. The notion of inhomogeneous chemical evolution therefore seems to be well established.
In the present work we use the stochastic chemical evolution model presented in Argast et al. (2000, hereafter Paper I) to investigate the enrichment of the ISM with r-process elements. In the following, we give a short summary of the model and discuss some major changes (mass infall, star formation rate, SN Ia events and r-process sites) that were implemented for this work.
We model the chemical enrichment of the halo ISM within a volume of (2.0
kpc)3, down to a resolution of (50 pc)3. Primordial matter is assumed to
be falling into this volume, following an infall law of the form
![]() |
(1) |
![]() |
(2) |
![]() |
(3) |
A crucial ingredient of chemical evolution models is the treatment of star
formation (SF). Unfortunately, it is also one of the least understood. In this
work we adopt a simple SF description based on Schmidt's law (Schmidt
1959). The number n(t) of stars that are formed per
time-step (of duration 106 yr) is determined by the product
The mass of a newly formed star is chosen randomly, subject to the condition
that the mass distribution of all stars follows a Salpeter initial mass
function (
)
with lower and upper mass limits of
and
,
respectively. The average mass
of a star is given by the
initial mass function (IMF) as
![]() |
(5) |
Low mass stars (
)
do not evolve significantly
during the considered time but serve to lock up part of the total mass, thus
affecting the abundances of elements with respect to hydrogen. Stars of
intermediate mass (
1-
)
return most of their
mass at the end of their stellar lifetime, leaving a white dwarf as stellar
remnant. Stars in the mass range of 10-
are assumed
to explode as SNe II, polluting the neighbouring ISM with their highly
enriched ejecta. Stellar lifetimes are taken from the Geneva Stellar
Evolution and Nucleosynthesis Group (cf. Schaller et al. 1992;
Schaerer et al. 1993a,b; Charbonnel et al. 1993). Stellar yields of O, Mg, Si, Ca and Fe are taken from
Thielemann et al. (1996) and Nomoto et al. (1997) and are scaled
according to Samland (1998) to account for the global chemical
enrichment of the Galaxy. In particular, Fe yields are reduced by a factor of
two.
SN events pollute the neighbouring ISM with their nucleosynthesis products and
sweep up the material in a spherical, chemically well mixed shell. Here, it is
assumed that each SN pollutes
of
ISM (Ryan et al. 1996; Shigeyama & Tsujimoto 1998). Stars which
form out of material enriched by a single SN II show an element
abundance pattern that is characteristic of the yields for the particular
progenitor of this SN II. This will lead to a large scatter in element
abundance ratios ([el/Fe]) as long as local inhomogeneities caused by SN II
events dominate the halo ISM. As time progresses, supernova remnants overlap
and the abundance pattern in each cell approaches the average defined by SN II
yield patterns for different progenitor masses and the IMF. This leads to a
decrease of the scatter in element abundance ratios at later times.
Table 1: Parameter values of the standard model.
To determine which intermediate mass stars form Type Ia SN events, we adopt
the following simple procedure: with probability
a newly
formed intermediate mass star has a companion in the same mass range and one
of these two stars will end its life as SN Ia. The mass of the companion
(again in the range 1-
and following a Salpeter IMF)
is determined randomly and the SN Ia event occurs after the less massive of
the two stars enters the red giant phase. Although this procedure admittedly
is rather simple, it has the advantage that the SN Ia frequency is determined
by only one free parameter, namely the probability
.
This
parameter is chosen in such a way that the slope observed in [
/Fe]
abundance ratios at
is reproduced. We chose
as upper mass limit for SN Ia progenitors to be consistent
with the lower mass limit of SN II. Since we use such a simplified description
to determine the occurrence of SN Ia, the upper mass limit of SN Ia
progenitors is not very relevant: assuming a Salpeter IMF, more than 99% of
the binary systems will consist of stars with masses
each. Finally, the yields of Type Ia SNe are taken from
Iwamoto et al. (1999, Model CDD2). Note, that if revised,
i.e. reduced electron capture rates are applied (see Brachwitz et al. 2000), the yields of model CDD2 correspond to those of model
WDD2, which are more accurate in terms of central conditions.
In Table 1 we list the parameter values adopted for our standard model.
An ideal element to trace the r-process enrichment of the ISM is the pure
r-process element Europium; approximately 97% of the solar Eu abundance was
produced in r-process events (Burris et al. 2000). Unfortunately, only
a small sample of Eu abundances at very low metallicities (i.e. [Fe/H] )
are available to date. In order to trace the r-process enrichment at
lower metallicities, the well studied element Barium is also used in this
investigation. Ba abundances in stars are dominated by the s-process, and only
15% of the solar Ba abundance was formed by rapid neutron capture
(Burris et al. 2000). However, the r-process fraction [Ba
/Fe] of Ba
abundances in stars can be easily calculated by subtracting the s-process
component from its total Ba abundance (cf. Burris et al. 2000). Since
we did not include s-process sources in our model, displayed Ba
abundances
in Sect. 3 of model and halo stars only show the r-process
contribution to the total Ba abundances.
In the following, r-process yields of Eu and Ba are estimated under the
assumption of a robust r-process for nuclei more massive than Ba, i.e. only
one source is responsible for the enrichment of the ISM with r-process
elements beyond Z=56. In this work we focus only on these heavier neutron
capture elements, since the enrichment of lighter neutron capture elements may
require a second, independent source (e.g. Wasserburg et al. 1996).
Thus, r-process matter ejected in an r-process event is considered to consist
of pure r-process nuclei more massive than Ba. First, the mass fractions of Ba
and Eu are calculated with the help of Table 5 in Burris et al. (2000)
and the standard solar system element abundances of Anders & Grevesse
(1989):
![]() |
(7) |
![]() |
(8) |
We now shortly discuss the implementation of lower-mass SNe (8-
), higher-mass SNe (
)
and NSM
as r-process sites in our model. We did not consider a hybrid model, where
different r-process sites contribute to the enrichment of the ISM with neutron
capture elements. Instead, it is assumed that for each of the three cases the
corresponding site is responsible for the whole r-process element inventory
(with
)
of the Galaxy, i.e. is the major source of r-process
elements beyond Ba.
Lower mass core-collapse SNe (8-10
)
are suspected to be
major r-process sources without contributing much iron-peak elements to the
enrichment of the ISM (Wheeler et al. 1998). Recently, Sumiyoshi et al. (2001) presented results of r-process nucleosynthesis occurring
during the artificially induced prompt explosion of an
star. The resulting distribution of r-process yields is in
reasonable agreement with solar r-process abundances for the heavier elements,
whereas lighter r-process elements (A <100) are underproduced. Since
elements in this mass-regime are often overproduced in neutrino-driven wind
models, the authors speculate that lower-mass SNe II might be responsible for
the production of massive r-process nuclei, whereas higher-mass SNe II account
for the less massive r-process nuclei. Wanajo et al. (2003) also find
r-process nucleosynthesis during the prompt explosion of progenitors with
masses in the range 8-10
which nicely reproduces the
second and third peak of the solar r-process abundances. Note however, that
it is still doubtful if the prompt explosion of a massive star occurs in
reality (e.g. Bruenn 1989a,b; Liebendörfer et al. 2001). Chemical evolution models including lower mass SNe as r-process source were presented by Ishimaru & Wanajo (1999) and Travaglio et al. (1999).
In the treatment of such lower-mass SNe, we proceed similarly to Ishimaru &
Wanajo (1999): in model SN810, r-process nucleosynthesis is assumed to
occur in the mass range 8-10
with constant Ba and Eu
yields over the whole mass range and assuming a Salpeter IMF. The yields then
are deduced from the average [Eu/Fe] and [Ba
/Fe] ratios (both
)
of metal-poor halo stars:
and
,
which equals a total of
of ejected r-process matter per event (Eqs. (9) and
(10)). Furthermore, we assume that the amount of
- and iron peak
elements synthesized in these SN II events are negligible (Hillebrandt et al. 1984).
Nucleosynthesis of r-process elements might also occur in neutrino-driven
winds or jets from nascent neutron stars during the delayed explosion of
high-mass stars with masses
(e.g. Woosley &
Hoffman 1992; Takahashi et al. 1994; Woosley et al. 1994;
Qian & Woosley 1996; Thompson et al. 2001; Nagataki
2001; Wanajo et al. 2001, 2002; Terasawa et al.
2002). The physical conditions needed for a robust r-process in
neutrino driven winds (such as high entropies, low electron fraction and/or
short dynamical timescales) are hard to achieve in present models, and lighter
r-process nuclei are often overproduced (e.g. Thompson et al. 2001). As
an example, Wanajo et al. (2002) require a very massive and compact
proto-neutron star of
to reproduce the solar system
r-process abundance pattern of heavier neutron capture elements. By allowing
for strong magnetic fields, Thompson (2003) finds that the required
physical conditions are much more easily achieved than in previous
models. Yet, there exist still some considerable theoretical uncertainties and
only very high-mass stars seem to be capable of providing the environment for
a robust r-process which reproduces the abundance pattern of neutron-capture
elements beyond Ba.
Williams (1987) reported the detection of Ba absorption lines in the
spectra of SN 1987A (progenitor mass
), and
Mazzali et al. (1992) noted that Ba was lacking at the surface of the
ejecta. Because of this remarkable feature, Tsujimoto & Shigeyama
(2001) identify SN 1987A as r-process site and estimate that
of Ba were synthesized during the SN
event. However, the poor knowledge of the far UV radiation field in the
envelope of SN 1987A makes it difficult to determine a truly reliable Ba
abundance in its ejecta (Utrobin & Chugai 2002) and it is not yet
established without doubt that the detected Ba was formed in the SN
event. Complementary to the observation of Ba in SN 1987A, Tsujimoto et
al. (2000) deduce Ba and Eu yields from observations of metal-poor halo
stars and inhomogeneous chemical evolution models, suggesting that
core-collapse SNe in the mass range 20-25
dominate the
production of r-process elements. They propose a Ba yield of
for a
and
for a
progenitor and a
Eu yield of
and
for a 20 and
progenitor,
respectively.
In a first model (SN2025), we adopt yields similar to the ones given by
Tsujimoto et al. (2000) and assume that SNe outside the indicated mass
range (i.e. 20-
)
do not contribute significantly to
r-process nucleosynthesis. In this model, a
star
produces
of Ba and
of Eu, whereas a
star yields
of Ba and
of Eu. For stars between 20 and
,
the yields are interpolated linearly. Thus, a SN II in
this mass range ejects on average
of r-process matter. Compared to the r-process yields
proposed by Tsujimoto et al. (2000) our yields are reduced by a factor
of two. This is the case, since we had to scale the Fe yields of Thielemann et al. (1996) by the same factor to account for the global chemical enrichment of the Galaxy.
As an alternative to model SN2025, we also calculated the chemical evolution
of the ISM with r-process yields from SNe II in the mass range 20-
(model SN2050). The yields were chosen in such a way that
the whole range of r-process abundances in metal-poor halo stars are
reproduced in the model. The r-process yield of a
star was set to
,
which is about three times the r-process
yield of a
star in model SN2025. This large amount
of ejected r-process matter requires a rapid decline in the r-process yields
of progenitors with masses
.
Otherwise, the average [r/Fe] ratios of model stars do not
reproduce the observed average abundances of metal-poor halo stars ([r/Fe
-0.5). Yields of progenitors in the mass range 28-
then are assumed to be constant. On average,
of r-process matter are ejected in each event.
Table 2 lists our Ba yields
and Eu
yields
and total ejected r-process matter
as function of progenitor mass m of models SN810, SN2025 and SN2050.
We point out that the yields adopted in this work are chosen in such a way
that observational constraints from metal-poor halo stars are respected.
Whether these yields are feasible would have to be tested with the help of
self-consistent core-collapse SN models, which, unfortunately, are not
available to date.
Table 2:
Adopted Ba (
)
and Eu
(
)
yields and ejected r-process
matter (
)
as function of progenitor mass m of models SN810,
SN2025 and SN2050. Yields in the mass range 8-
are assumed to be constant (SN810). Yields in the
mass range 20-
are linearly interpolated
(SN2025). Yields of model SN2050 exhibit a more complex behaviour: a
rapid decline in the mass range 20-
is
assumed. Progenitors with masses 28-
all
have the same, constant yield.
Another major source of r-process elements might be neutron star mergers
(Freiburghaus et al. 1999b; Rosswog et al. 1999, 2000). A
scheme similar to the determination of SN Ia events was applied for the
determination of their occurrence (cf. Sect. 2.1): with probability
a newly formed high mass star (in the range 10-
)
has a companion in the same mass range. It is assumed
that, subsequent to both SN II events, the remaining neutron stars will
eventually merge. The time
that is needed for the coalescence of the two
neutron stars and the probability
for the occurrence of NSM
events are treated as free parameters which can be adjusted within given
observational estimates on the coalescence rate of binary compact objects and
merger timescales (van den Heuvel & Lorimer 1996; Kalogera & Lorimer
2000; Belczynski et al. 2002). As was the case for SN Ia events,
this treatment is very simple, but has the advantage that all difficulties
associated with the determination of formation channels of neutron star
binaries (or neutron star / black hole binaries) and the corresponding
formation rates are combined in only one free parameter. Furthermore, the
lower mass boundary for neutron star formation (here
)
is of little importance: the decisive term for the
enrichment of the ISM with r-process material in our model is the NSM
rate. Since a decrease of the lower mass boundary results in an increase of
the NSM rate, a simple adjustment of the parameter
is
sufficient to retrieve a NSM rate consistent with constraints on the galactic
rate. In the model the ejected r-process matter is also mixed with
of ISM and is assumed to be distributed in a
spherical, chemically mixed shell. This might not be true, however, since it
is well conceivable that r-process matter is ejected only in the orbital plane
of the merging neutron stars. In this case, the ejecta would be distributed
over a smaller volume and would consequently mix with a smaller amount of ISM
than the one assumed in our model. As will be seen below (Sect. 3.2),
this would lead to even higher [r/Fe] ratios in the computed model
stars. Thus, our assumptions are in favour of the NSM scenario and changing
the mixing mechanism in our model would only aggravate the problems associated
with their large r-process yields. NSM as major r-process sources have also
been studied by De Donder & Vanbeveren (2003). Contrary to their work,
we investigate the ISM enrichment by NSM in the framework of inhomogeneous
chemical evolution.
Coalescence timescales for neutron star mergers are typically estimated to be
of the order of 100-1000 Myr (e.g. Portegies Zwart & Yungel'son
1998; Fryer et al. 1999). Recently, Belczynski et al. (2002) suggested a dominating population of short-lived
neutron star binaries with merger times less than 1 Myr. This population of
neutron star mergers might be formed through channels involving mass-transfer
episodes from helium stars, leading to tightly bound binary systems with very
short orbit decay timescales. However, this result depends on the occurrence
of a common envelope (CE) phase of the progenitor He stars, which was treated
in a simplified manner by Belczynski et al. (2002). The authors caution
that detailed hydrodynamical calculations of the CE phase still have to
confirm these results. In view of these uncertainties, we adopt two different
coalescence timescales ,
namely 1 and 100 Myr. In each case, it is
assumed that all neutron star binaries merge within this
time. Coalescence timescales of the order 1 Gyr have not been considered since
they are not consistent with observations of neutron capture elements in ultra
metal-poor stars (cf. Sect 3.2).
Table 3:
Constraints on the Galactic NSM rate
and ejected
r-process matter
.
Constraints on the Galactic NSM rate are also controversial. Van den Heuvel &
Lorimer (1996) estimate a rate of approximately
,
Kalogera & Lorimer (2000) give an upper limit of
whereas Belczynski et al. (2002) get rates in the range
.
In Newtonian calculations, the
amount
of r-process matter ejected in a NSM event is of the order of a
few times 10-3-
,
depending on the initial
configuration of the binary system (Rosswog et al. 1999). Taking
general relativistic effects into account, Oechslin et al. (2002) find
-
as lower limit for
.
All constraints are again listed in Table 3.
Table 4:
Parameter values adopted for NSM events (see text for details). For
each value of
,
two models with the indicated
coalescence timescales
have been calculated.
![]() |
Figure 1:
Evolution of [Eu/Fe] and [Ba![]() ![]() |
Open with DEXTER |
Table 4 lists the parameter values adopted for the NSM
probability
,
the resulting average NSM rate
,
the amount of r-process matter
ejected in each event
and coalescence timescales
.
After a value for
has
been chosen, the coalescence timescale is the only free parameter remaining in
the model since the actual NSM rate is determined by the NSM probability and
the SF rate (resulting from the parameter values given in Table 1).
In addition, the NSM rate and the amount of ejected r-process matter are
tightly correlated, since the total amount of r-process matter in the Galaxy
(
,
Wallerstein et al. 1997) has to
be reproduced. Thus, higher NSM rates require that less r-process matter is
ejected in each event, and vice versa.
Comparing the values in Table 4 with the constraints listed in
Table 3 reveals that the NSM rate of the first row in
Table 4 is too high by a factor of 10-100 and that the
corresponding low value of
of ejected r-process
matter is close to the lower limit given by Oechslin et al. (2002). NSM
rates listed in the three lower rows seem to be consistent with the galactic
NSM rate. However, the ejected r-process matter in the last row is clearly
beyond the upper limit of allowed values.
![]() |
Figure 2:
Evolution of [Eu/Fe] and [Ba![]() ![]() |
Open with DEXTER |
![]() |
Figure 3:
Evolution of [Eu/Fe] and [Ba![]() ![]() |
Open with DEXTER |
In this section, the enrichment of the ISM with neutron capture elements is
discussed under the assumption, that the dominating r-process sources are
either lower-mass SNe II (8-
)
or higher-mass SNe II
(>
). The results of models SN810, SN2025 and SN2050
are shown in Figs. 1-3, respectively. The
figures show the evolution of [Eu/Fe] and [Ba
/Fe] as a function of
metallicity [Fe/H]. Model stars are shown as black dots, whereas observations
are indicated by filled red squares and diamonds. Observations are taken from
Peterson et al. (1990); Gratton & Sneden (1991a,b);
Ryan et al. (1991); Edvardsson et al. (1993); François et al. (1993); Beveridge & Sneden (1994); McWilliam et al. (1995a); Ryan et al. (1996); Jehin et al. (1999);
Aoki et al. (2000); Burris et al. (2000); Mashonkina & Gehren
(2000, 2001); Sneden et al. (2000a); Mishenina & Kovtyukh (2001); Koch & Edvardsson
(2002) and Stephens & Boesgaard (2002). Single observations of
stars are marked by a square. In the case where multiple observations of a
star are present, we plot the most recent one if all observations were
published before the year 2000 (also marked by squares). If several more
recent observations are available, the given element abundances are averaged
(marked by diamonds). The yellow line shows the average element abundances in
the model ISM and can directly be compared to classical chemical evolution
models, which assume that the ISM is well mixed at all times. Purple filled
circles with error bars denote the average [r/Fe] ratios of model stars and
their standard deviation in [Fe/H] bins with binsize 0.1 dex. Note, that the
[Ba
/Fe] plots only show the r-process contribution to the total Ba
abundances of halo stars. According to Burris et al. (2000), r-process
Ba abundances can be computed by removing the s-process contribution to Ba in
stars with [Fe/H] > -2.5, if Eu abundances have also been determined. For
stars with metallicities [Fe/H] < -2.5, it can be assumed that the whole Ba
inventory is of pure r-process origin. Thus, published Ba abundances of such
metal-poor halo stars simply have to be renormalized to the level of the solar
r-process fraction (i.e.
15% of the total solar Ba abundance). Care
has been taken to remove known carbon stars from our sample. Such stars mostly
show unusually large Ba abundances, which are thought to originate from mass
transfer of s-process enriched matter in binary systems.
The evolution of r-process elements shown in
Figs. 1-3 are qualitatively very similar. At
very low metallicities ([Fe/H
), a large scatter in abundances of
model stars is visible. This scatter is due to chemical inhomogeneities in the
early ISM (cf. Paper I). The scatter decreases as the mixing of the ISM
improves and finally reaches the IMF averaged mean. At this stage, the ISM can
be considered well mixed and the further evolution is comparable to the one of
classic chemical evolution models. In the following, we point out some
important features of the ISM enrichment resulting from models SN810, SN2025
and SN2050:
On the other hand, Qian & Wasserburg (2001, 2002) argue
that the highest Eu and Ba abundances in metal-poor halo stars may originate
from the dumping of r-process enriched matter from a high-mass companion onto
the surface of a low mass star. This explanation requires that during the
SN II event of the high-mass companion, r-process elements but no iron-peak
nuclei are produced. Possible candidates for this kind of SN II events are
progenitors in the mass range 8-10
;
see e.g. Qian &
Wasserburg (2002) for a discussion of this mechanism. To date, our
model is not able to handle the occurrence of surface contamination in binary
systems, so that the results of model SN810 do not directly apply to this
scenario. We speculate however, that if this mechanism would be at work during
the evolution of the Galaxy, we might also expect stars with high [r/Fe]
abundances at metallicities [Fe/H] > -2. No such stars have been observed to
date. Yet, since many abundance studies primarily target metal-poor stars, it
is possible that r-process enriched, metal-rich stars simply have been
missed. Note, that if iron is produced in the SN II event of the companion,
the resulting surface abundance will be a mixture of the intrinsic r-process
abundances and the yields of the SN II. In this case, large r-process yields
are still required to enrich the companion to the high level seen in
observations;
![]() |
Figure 4:
Evolution of [Eu/Fe] and [Ba![]() ![]() ![]() |
Open with DEXTER |
![]() |
Figure 5:
[Ba![]() ![]() ![]() ![]() ![]() ![]() |
Open with DEXTER |
We briefly note here, that the fit of model SN810 to observations of
metal-poor halo stars can be improved, if infall of pre-enriched ISM with
[Fe/H] =-3 is assumed (e.g. Qian & Wasserburg 2002). In this case,
model stars with high [r/Fe] abundances at [Fe/H
are shifted to
[Fe/H] =-3, where such r-process enriched stars are observed. This is the
case only for model SN810. The fit to observations of models SN2025, SN2050
and the NSM models (see Sect. 3.2) are not changed significantly by
the assumption of metal-rich infall;
This is demonstrated in Fig. 1 (model SN810). Although the
calculated standard deviation of model stars is large (due to the low number
of model stars with metallicities below [Fe/H] <-3), the average model star
abundances rise with decreasing metallicity. Thus, model SN810 would predict a
rising slope in [Ba/Fe], contrary to the one present in observations of
metal-poor halo stars. In model SN2025, the average model abundances stay more
or less constant, whereas in model SN2050, average abundances give a nice fit
to the observational data. Note, that the empty bin at [Fe/H] = -4 in
Fig. 2 only contains one model star.
![]() |
Figure 6:
Same as Fig. 5 but with
![]() |
Open with DEXTER |
The enrichment of the ISM with r-process elements in the case of NSM acting as
major r-process sources is discussed in the following. As already mentioned in
Sect. 2.2.3, eight model-runs were carried out in total, assuming
NSM rates ranging from
to
and coalescence time-scales of 1 and 100 Myr (see
Table 4). A representative example of the evolution of [Eu/Fe]
and [Ba
/Fe] as a function of [Fe/H] in the NSM case is shown in
Fig. 4. A condensed overview of all models displaying only
the evolution of [Ba
/Fe] is shown in Figs. 5 and
6. The evolution of r-process abundances is strikingly
different from the case in which r-process nucleosynthesis occurs in SNe II
(see Figs. 1-3). All models with NSM as
r-process sources fail to reproduce observations.
In Fig. 4, a NSM rate of
,
a total of
of ejected
r-process matter and a coalescence time-scale
of 1 Myr was adopted.
Note, that the case discussed here assumes a NSM rate which is at the upper
limit given by present estimates of the galactic NSM rate: estimates range
from
(van den Heuvel & Lorimer
1996) to
(Belczynski et al. 2002) and Kalogera & Lorimer (2000) even
give an upper limit of
.
Furthermore, the coalescence timescale
adopted in this model is only 1 Myr, which is two to three orders of magnitude
lower than the classical estimate of 100-1000 Myr (e.g. Portegies Zwart &
Yungel'son 1998; Fryer et al. 1999). However, Belczynski et al.
(2002) suggest the existence of a dominating population of short
lived neutron star binaries with merger times less than 1 Myr.
![]() |
Figure 7:
![]() ![]() ![]() ![]() ![]() ![]() |
Open with DEXTER |
The following qualitative differences are immediately visible, if Figs. 1-3 are compared to Fig. 4:
In the upper left panels of both figures (NSM rate
), a few model stars with r-process abundances first appear
around [Fe/H
,
in contrast to [Fe/H
in the panels
at the lower right (NSM rate
).
Simultaneously, the scatter in [Ba
/Fe] and [Eu/Fe] at solar metallicity,
which is of the same order as the one observed in the upper left panels,
increases to almost 3 dex in the lower right panels, clearly not consistent
with observations. Furthermore, the tail of model stars with low [r/Fe]
abundances gets more pronounced and concentrated and is shifted to higher
metallicities.
Note, that the upper left panel of Fig. 5 gives the best fit to
observations in the NSM case, although model stars with r-process
contributions at [Fe/H] < -3 are still somewhat scarce. This model was
calculated with a NSM rates of
,
of ejected r-process matter and a coalescence timescale
of 1 Myr. It is not very surprising that this model gives the best fit to
observations, since its parameters (i.e. rates, ejected r-process matter) are
comparable to the SN II models. A further increase of the NSM rate and,
consequently, a decrease of the amount of ejected r-process matter, would even
improve the fit. However, a comparison with Table 3 shows, that
the adopted NSM rate in this case is one order of magnitude larger than the
estimate by Belczynski et al. (2002), whereas the amount of ejected
r-process matter is of the same order as the lower limit given by Oechslin et al. (2002).
Contrary to expectations, the impact of the coalescence timescale
on the
distribution of r-process nuclei is rather small. This is surprising, since a
large value of
results in a delayed injection of r-process matter into
the ISM. However, this can be understood by noting that in the NSM models, the
enrichment of the ISM was adjusted to proceed very slowly at the beginning of
Galaxy formation: it takes
200 Myr to enrich the ISM from [Fe/H] =
-4 to [Fe/H] = -3. Thus, a coalescence timescale of 100 Myr is of the same
order as the halo enrichment timescale, resulting in only a slight shift of
the first r-process matter injection to higher metallicities. The slow
enrichment was assumed in order to enable r-process nuclei injection from NSM
even at low metallicities. However, as Figs. 5 and
6 show, this is not the case even with the help of a slow
enrichment history. Thus, a fast halo enrichment would shift the moment of
first r-process matter injection to even higher metallicities. Note that, on
the basis of these considerations, coalescence timescales of the order 1 Gyr
clearly are not consistent with the occurrence of r-process nuclei early on
during Galaxy formation.
Figure 7 demonstrates another qualitative difference between the
SN II and NSM scenarios. In this figure,
vs. [Fe/H] of model SN2050 (left panel) and NSM occurring at a
rate of
yr-1 (right panel) are shown. In the SN II
scenario, the scatter in
at low
metallicities ([Fe/H] < -2) is primarily determined by the progenitor mass
dependent r-process yields of SNe II (indicated by numbered circles). Model
stars at [Fe/H] < -2 with Ba
abundances above
are formed out of material enriched by a single
SN II whereas some model stars with
are formed out of material which underwent subsequent mixing
with primordial ISM. Stars with such low r-process abundances disappear at
[Fe/H] > -2 due to the increased enrichment of the ISM with r-process
material.
In the NSM case, the scatter in
is solely determined by subsequent mixing of r-process enriched material with
r-process poor ISM. Under the assumption that a NSM event is able to enrich
of ISM with
of
r-process matter, model stars with typically
or more can be formed (indicated by the labelled
circle in the right-hand panel of Fig. 7). Thus, the distinct
peaks in the distribution of Ba
abundances of model stars are telltale
signs of recent NSM events. However, since NSM events are rare, there is time
enough to thoroughly mix the r-process enriched material in the vicinity of a
recent NSM with r-process poor ISM. In consequence of this mixing and the low
NSM rate, the bulk of model stars with
are formed.
We conclude that NSM have to respect stringent requirements if they are to
operate as dominating r-process contributors: High NSM rates of the order
are required to account for the early appearance
of r-process nuclei in the ISM. Consequently, the amount of r-process matter
ejected in a NSM event has to be of the order
,
lest the total inventory of r-process nuclei in the Galaxy
is exceeded. In addition, the coalescence timescale of a large fraction of
neutron star binaries has to be of the order 1-10 Myr or at the most 100
Myr. If these constraints are indeed fulfilled in reality, then the properties
of NSM (i.e. rates and ejected r-process matter) are comparable to the ones of
SNe II that (may) give rise to r-process nucleosynthesis from the point of
view of chemical evolution. Consequently, the enrichment of the ISM with
r-process matter from NSM would be similar to the one where core-collapse SNe
are the source of r-process elements (compare
Figs. 1-2 with the upper left panels of
Figs. 5 and 6).
At present, it seems unlikely that NSM fulfil these conditions, which suggests
that NSM are ruled out as the major r-process source. It should be noted,
however, that this conclusion is strictly based upon the current knowledge of
NSM rates, coalescence timescales, the amount of ejected r-process matter,
etc. If these estimates should undergo mayor revisions in the future, a
reanalysis of NSM as major r-process sources may also be required. Yet, even
with the present estimates, NSM occurring at low rates and with low ejecta
masses of the order 10-4-
may still
contribute to the enrichment of the ISM with r-process nuclei if their
r-process signature is similar to the one generated by SNe II. In this case,
the impact of NSM on r-process nuclei enrichment would be negligible compared
to the fast injection of r-process matter by SNe II.
In this work, we study the enrichment of the interstellar medium with
r-process elements in the framework of inhomogeneous chemical evolution. We
present a detailed comparison of the impact of lower-mass SNe II (8-
,
model SN810), higher-mass SNe II (
,
models SN2025 and SN2050) and NSM as major r-process sites on the enrichment history of the early Galaxy.
Among the SNe II scenarios, model SN2050 gives the best fit to observations,
since it reproduces the trend of r-process abundances at ultra low
metallicities and at the same time fits observations with the highest [r/Fe]
abundance ratios at [Fe/H
.
Neither model SN2025 nor SN810 are
able to reproduce observations of the most r-process enriched stars. In
addition, model SN810 gives the worst fit to the trend of r-process abundances
at ultra low metallicities. However, we conclude that, due to the large
uncertainties inherent in the progenitor mass dependence of iron yields of
SNe II, it is not possible to clearly rule out either lower-mass SNe II or
higher-mass SNe II as dominant r-process sites from the point of view of
inhomogeneous chemical evolution. Additional uncertainties are introduced by
the fact that reliable r-process yields from SNe II are unavailable as
yet. Here, they were deduced in such a way that the average [r/Fe] abundances
in metal-poor halo stars are reproduced. Clearly, iron and r-process yields
and their dependence on progenitor mass from self-consistent core-collapse SN
calculations are required before any decisive conclusion can be reached.
On the other hand, NSM seem to be ruled out as major r-process sources for the following reasons:
Thus, we conclude, that the exact astrophysical nature of r-process
sites still remains a mystery, since it is not possible to clearly distinguish
between neutron capture element abundance patterns resulting from lower-mass
SNe II (8-10
)
and the ones from higher-mass SNe II
(
)
in the framework of inhomogeneous chemical
evolution. However, the present investigation suggests that core-collapse SNe
are much more likely to be the dominant r-process sites than coalescing
neutron star binaries, which at least reduces the list of possible major
contributors of r-process nuclei to the enrichment of the interstellar medium.
Yet, it remains to be seen how SNe II can actually produce the required
r-process yields.
Acknowledgements
The authors thank the referee John J. Cowan for his constructive remarks which helped to improve this paper. We would also like to thank D. Vanbeveren and S. Wanajo, who have provided helpful discussions. This work was supported in part by the Swiss National Science Foundation and by US DOE grants DE-FG02-87ER40328 and DE-FG02-00ER41149 at Minnesota.