A&A 454, 917-931 (2006)
DOI: 10.1051/0004-6361:20064968
Instituto de Astronomia, Geofísica e Ciências Atmosféricas, Universidade de São Paulo, Rua do Matão 1226, 05508-900 São Paulo, Brazil
Received 4 February 2006 / Accepted 16 March 2006
Abstract
Context. Barium stars show enhanced abundances for the slow neutron capture (s-process) heavy elements, so they are suitable objects for studying s-process elements.
Aims. The aim of this work is to quantify the contributions of the s-, r-, and p-processes for the total abundance of heavy elements from abundances derived for a sample of 26 barium stars. The abundance ratios between these processes and neutron exposures were studied.
Methods. The abundances of the sample stars were compared to those of normal stars, thus identifying the fraction relative to the main component of the s-process.
Results. The fittings of the
curves (neutron-capture cross-section times abundance, plotted against atomic mass number) for the sample stars suggest that the material from the companion asymptotic giant branch star had approximately the solar isotopic composition as concerns fractions of abundances relative to the s-process main component. The abundance ratios of heavy elements, hs, ls, and s and the computed neutron exposure are similar to those of post-AGB stars. For some sample stars, an exponential neutron exposure fits the observed data well, whereas a single neutron exposure provides a better fit for others.
Conclusions. The comparison of barium and AGB stars supports the hypothesis of binarity for the barium star formation. Abundances of r-elements that are part of the s-process path in barium stars are usually higher than those in normal stars, so barium stars also seemed to be enriched in r-elements, although to a lower degree than s-elements. No dependence on luminosity classes was found in the abundance-ratio behaviour among the dwarfs and giants of the sample of barium stars.
Key words: nuclear reactions, nucleosynthesis, abundances - stars: abundances
Three of the eight processes defined by B2FH are the focus of the present work: the s-, r-, and p-processes. We present the abundance ratios between these processes and abundance distribution for a number of heavy elements in a sample of barium stars. The abundance determination for the sample stars was described in Allen & Barbuy (2006, Paper I hereafter).
The abundance ratios for the best representatives of s- and r-processes provide clues to the formation and evolution of the Galaxy, since each of them is related to a different formation site, including stars of different characteristics and evolutionary stages. As an example, Mashonkina et al. (2003) estimate the timescale for the thick disk and halo formation based on the abundance ratios of [Eu/Ba], [Mg/Fe], and [Eu/Fe] and on the calculations of chemical evolution by Travaglio et al. (1999). Burris et al. (2000) conclude from abundances of metal-poor giant stars that the contributions of s-processes can be seen at metallicities as low as [Fe/H] = -2.75, and they are present in stars of [Fe/H] > -2.3, indicating that the s-process acts at lower metallicities than predicted by previous works.
Considering the scenario of barium star formation through enriched material transfer from a companion AGB star, it is worth comparing abundances of AGBs and barium stars. Such study can provide clues to the enrichment of barium stars in heavy elements, given that they should conserve the surface characteristics of AGB stars.
This work intends to address the following questions:
1. What are the s-element abundances of normal stars?
2. Considering the hypothesis of the transfer of material enriched in carbon and s-elements from a more evolved companion for barium star formation, is it possible to estimate from observed abundances which proportions these elements are received?
3. What is the isotopic composition of the material received? Is it possible to consider the solar system mix?
4. How is the
curve behaviour for barium stars?
5. What kind of neutron exposures were involved in the nucleosynthesis of the transferred material?
This paper is organised as follows. Section 2 is a brief description of atmospheric parameters and of the abundance determination carried out in Paper I; Sect. 3 briefly explains the s-, r-, and p-processes; Sect. 4 shows the contributions of s-, r-, and p-processes in normal stars compared to barium stars; Sect. 5 evaluates two kinds of neutron exposures and shows ratios involving s, hs, and ls. In Sect. 6 conclusions are drawn.
Photometric data were taken in several runs at the ZEISS 60 cm telescope at the LNA (Laboratório Nacional de Astrofísica) and from the literature. Optical spectra were obtained at the 1.52 m telescope at ESO, La Silla, using the Fiber Fed Extended Range Optical Spectrograph (FEROS) (Kaufer et al. 2000). A set of atmospheric parameters (temperature, surface gravity, metallicities, and microturbulent velocity) was obtained in an iterative way.
Photospheric 1D models were extracted from the NMARCS grid (Plez et al. 1992)
for gravities
,
and from Edvardsson et al. (1993) for less evolved stars,
with
.
The LTE abundance analysis and the spectrum synthesis calculations
were performed using the codes by
Spite (1967, and subsequent improvements in the past thirty years),
described in Cayrel et al. (1991) and Barbuy et al. (2003). Line lists and respective
parameters are given in Paper I.
![]() |
Figure 1:
![]() ![]() ![]() ![]() ![]() ![]() |
Open with DEXTER |
![]() |
Figure 2: Same as Fig. 1 for other elements. |
Open with DEXTER |
Typical s-elements are mainly produced by the s-process, but a smaller fraction of their abundance is due to the r-process. Analogously, the major production of r-elements is due to the r-process with a smaller contribution from s-process. Beyond these main processes, the p-process contributes with a small fraction of heavy elements production, as shown in Fig. 1 of Meyer (1994).
In the s-process, neutrons are captured by seed nuclei, which are the iron-peak
elements, on a long timescale relative to
decay, denominated "s'' (slow)
by B
FH.
This process has been subdivided into three components according to the site and the
nucleosynthetic products: main, weak, and strong.
The s-process main component is responsible for many nuclides in the range of
atomic mass
.
Abundance peaks can be observed near
A = 90, 138, and 208.
In the classical analysis (Käppeler et al. 1989) the element formation through the s-process
occurs in a chain, starting with the seed nucleus
Fe.
Following this analysis, it is possible to obtain an expression for the
,
where
is the neutron capture cross-section and N the abundances due
to the s-process only for each nuclide. This expression involves the mean
distribution of neutron exposure
,
the abundance fraction of
Fe required
as seed to s-process G, and the solar abundance of
Fe N
,
For the solar system, the corresponding curve from classical analysis represents the abundances of s-process nuclei for A > 90 very well. For lighter nuclides, this line appears below the empirical data, suggesting another form of synthesis of s-process nuclides in stars, with the first one called the main component and the second, the weak component. The good agreement between the classical model and the observed data of this curve for the solar system represents an interesting characteristic of the s-process, taking the large number of nuclides between Fe and Bi into account.
The s-process main component is believed to occur during the
thermal pulses - asymptotic giant branch (TP-AGB) phase of intermediate or
low mass stars. In this phase, the star consists of an inert CO core and
the He and H burning shells. The region between the two shells is
the so-called He intershell, where neutrons released by the
13C (, n)16O reaction during the interpulse period and
by the 22Ne (
, n)25Mg reaction during the convective
thermal pulse are captured by iron-peak nuclei.
The s-process weak component is responsible for part of the abundance of nuclides with
atomic mass in the range
(Raiteri et al. 1993; Lamb et al. 1977). The
nucleosynthetic site is probably the core helium burning of stars with masses
10
,
where temperature is high enough for
the main neutron source to be the 22Ne (
, n)25Mg reaction.
Neutron density is low compared to the main component.
The s-process strong component was postulated in order to provide part of the Pb abundance (see Käppeler et al. 1989). However, according to Busso et al. (1999), it is possible to explain the 208Pb galactic abundance without it.
The r-process occurs in an environment that is rich in neutrons, where several
of them are
captured by nuclides on a short timescale compared to
decay, and for this
reason, this process was denominated as "r'' (rapid) by B
FH. In this case,
the neutron density is higher than for the s-process. Sites that favour
such high neutron density are the final stages of massive stars
such as core collapse supernovae (SN II, Ib, Ic)
(Qian 2000; Wasserburg & Qian 2000; Qian 2001) or those involving neutron stars
(Freiburghaus et al. 1999; Woosley & Hoffman 1992; Meyer 1994; Rosswog et al. 2000,1999).
The fact that r-elements are observed in very metal-poor stars suggests that these
elements were produced in supernovae events resulting from the evolution of the
first massive stars in the Galaxy (Cowan et al. 2002; Sneden et al. 2003; Hill et al. 2002; Sneden et al. 1996; Ishimaru et al. 2004; Honda et al. 2004).
Despite the promise of these scenarios, some difficulties have been found; for instance,
Wanajo et al. (2001) show that to reproduce a solar abundance of r-elements,
proto neutron stars must have 2
and 10 km of radius, characteristics
not observed so far.
The p-process forms nuclei rich in protons. Some s- or r-nuclei, where s- and
r-processes were blocked, capture protons with emission (p,
).
The p-nuclei may also be synthesised by photodesintegration (
,
n) of
a pre-existing nucleus rich in neutrons (especially s-nuclei), followed by
possible cascades of (
,
p) and/or (
,
)
reactions.
The p-process site should be rich in hydrogen, with proton density
10
g/cm
,
at temperatures of
K.
This process site is most likely related to SN II, according to Arnould (1976),
Woosley & Howard (1978, and references therein), Arnould et al. (1992), and
Rayet et al. (1993) or the explosion of a moderately massive white dwarf
due to the accretion of He-rich matter
(Goriely et al. 2005, and references therein).
A quantitative study of the p-process is presented in Rayet et al. (1995).
One way to quantify the enrichment of s-elements in barium stars is to compare
their abundances with normal stars.
For the s-, r-, and p-elements, the total abundance for the element can be described
by the sum of the abundances corresponding to the three nucleosynthetic processes,
taking into account all contributing isotopes "i'':
Table 1:
Results obtained from least-square fitting of
vs. [Fe/H] of normal stars:
= A[Fe/H] + B;
"cov'' is
the covariance between A and B; "d.o.f.'' is the number of degrees of freedom.
From
vs. [Fe/H] fittings, it is possible
to determine the total abundance of a certain element in a normal star of a
given metallicity. The first line of Table 2 shows
obtained directly from fittings for normal stars with metallicities corresponding to
the barium star, indicated in the header.
Data for molibdenium are rarely available in the literature in the same metallicity range
as the present sample, and for this reason the least-square fitting was done
using data derived for the globular cluster
Centauri by Smith et al. (2000), with results
shown in Col. 2 of Table 3. In order to verify the reliability
of these results, Mo abundances were calculated through
(Mo) =
(Mo) + [Fe/H] + [Mo/Fe], considering
[Mo/Fe] = 0, and results are shown in Col. 3 of Table 3. The agreement
between results from these two columns means that, according to the fitting,
[Mo/Fe]
0 at metallicities near solar. Data for Gd and Pb are also
rare in the literature for normal stars in the range of metallicities
of the present sample.
For this reason,
(Gd) and
(Pb) were
determined by considering [Gd/Fe] = [Pb/Fe] =
,
near the solar metallicity.
Gadolinium is expected to behave like dysprosium given that both are produced mainly through the r-process in
almost the same proportions (Arlandini et al. 1999). Columns 4 and 5 of Table 3
show
(Dy) values determined by least-square fitting and by
summing the star metallicity to the solar value, and Cols. 6 and 7 show
(Gd) and
(Pb).
Table 2:
Abundances obtained from least-square fittings for normal stars with metallicities
corresponding to our sample of barium stars. Each set of five lines corresponds to the normal star
with the same metallicity as the barium star indicated in parenthesis.
corresponds to
and
corresponds to
from
Eq. (2).
is from Eq. (6). The full table is only available in electronic form at the CDS.
Table 3:
Abundances [
(X)] of Mo, Dy, Gd, and Pb for normal stars.
Symbol "*'' indicates that the
column was obtained from a least-square fitting.
Considering that the total abundance of an element is the sum of the contributions of
s-, r-, and p-processes, one can write:
Table 6:
Results relative to the s-, r-, and p-processes for barium stars.
,
,
and
:
abundance fractions of s-process
main, weak, and strong components;
and
:
abundance fractions of r- and p-processes;
and
(Si): cross-section times abundance fraction corresponding to the
s-process main component taking the overabundance of barium stars into account on the usual
scale and Si scale;
:
abundance fractions of the s-process main component for normal stars;
diff:
;
:
cross-section times abundance fraction corresponding to the
s-process main component without the overabundance of barium stars.
The full table is only available in electronic form at the CDS.
The uncertainty on
is calculated with
![]() |
(4) |
![]() |
(5) |
Finally, the sum of abundances from other processes except the s-process main component
was computed, i.e., the sum of fractions
sw, st, r, and p, indicated by subscript "o'' and its uncertainty are
![]() |
The abundance fraction corresponding to the s-process main component of an element is
Normal stars are expected to have lower abundances of heavy elements than barium stars
(
); however, for some elements in some stars,
abundances in barium stars obtained for the present sample were
much lower than those of normal stars. This is the case of
Mo (HD 27271, HD 116869, HD 123396, HD 210709,
HD 210910, HD 223938),
Eu (HR 107), Gd (HD 210709), Dy (HD 89948), and Pb (HD 22589, HD 210910).
It is not clear why these barium stars show such low abundances for these
elements. For these cases, the Eq. (8)
does not apply, and there is no
for them (Table 5).
After computing the abundance relative to the s-process main component for an element
in barium stars using Eq. (8), the abundance of each isotope
and for each process can be determined with
![]() |
(9) |
![]() |
(10) |
![]() |
(11) |
![]() |
(12) |
![]() |
(13) |
![]() |
Figure 3:
Total abundances of
heavy elements for the barium (
![]() ![]() ![]() |
Open with DEXTER |
![]() |
Figure 4:
Abundance fraction of heavy elements
due to all processes except the s-process main component (
![]() ![]() ![]() |
Open with DEXTER |
Some elements of Table 2 are formed in larger amounts through the r-process. They are
Eu (94.2%), Gd (84.54%), Dy (84.8%), and Sm (67.4%), according to Arlandini et al. (1999).
For these elements, barium stars are expected to have abundance values closer to normal
stars than for s-elements. It can be verified that
this supposition is valid for Eu by comparing data from Tables 16 and 17 of
Paper I to those from Table 2. Figure 3 shows the behaviour
of
of normal stars calculated by least-square fitting with
of
barium stars, and Fig. 4 shows the behaviour of
for barium
stars with
.
It is important to point out that Figs.
3 to 10, the elements were arranged in increasing order of
the contribution by the s-process main component, following Arlandini et al. (1999),
Eu, Gd, Dy, Sm, Pb, Pr, Mo, Nd, La, Ce, Ba, Zr, Sr, and Y. Behaviours tend to be approximately
constant; however, Eu data are very close to a straight line with tangent = 1,
differently from the other elements. The larger the distance of the data from tangent = 1,
the larger the abundance of barium stars compared to normal stars. If the fraction
corresponding to the s-process main component
is withdrawn as in Fig. 4, the behaviours of Eu, Gd, and Dy are almost
unaltered relative to Fig. 3; however, the change in ordinates
due to the missing abundance is remarkable for the other elements.
![]() |
Figure 5:
Abundances corresponding to the s-process
main component for heavy elements of barium stars (
![]() ![]() ![]() |
Open with DEXTER |
![]() |
Figure 6:
![]() ![]() |
Open with DEXTER |
Figure 5 shows the behaviour of
with
,
which are the abundance fractions of the s-process main component of barium stars
(Eq. (8)) and normal stars of the same metallicity, respectively.
Unlike in Fig. 4, the changes in ordinates due to the missing
fraction r in Eu, Gd, and Dy abundances in Fig. 5 are remarkable in comparison
with Fig. 3. Figure 5 shows similar differences from tangent = 1
as in Figs. 3 and 4, making it clear that the fraction s of the
abundance
of barium stars leads their behaviour. The difference between the maximum and minimum values of
the abundance fraction corresponding to the s-process main component of normal stars is
in the range of metallicities of the sample stars.
For barium stars this difference is larger,
.
According to Fig. 6, the abundance relative to the s-process main component of heavy elements is essentially independent of [Fe/H] for the present sample of barium stars.
![]() |
Figure 7:
![]() ![]() ![]() ![]() |
Open with DEXTER |
It is usual to investigate the starting point of the s-process contribution by analyzing
the [Eu/Ba] abundance ratio, since they are, respectively, the best representatives
of the r- and s-processes. Figure 6 from Burris et al. (2000) shows the
behaviour of Ba relative to Eu for normal stars of metallicities
.
From the moment that the s-process starts to produce Ba, an increase
in its abundance relative to Eu can be seen, since the s-process is
responsible for
81% of its production, according to Arlandini et al. (1999).
Figure 7 shows
(X) vs. [Fe/H] for barium stars from the
present sample and post-AGB stars from Reyniers et al. (2004) and van Winckel & Reyniers (2000).
Abundances of post-AGBs are larger or similar to those of barium stars,
except for Sm in one star. This is expected considering that barium stars were
enriched by an AGB companion.
Figure 8 shows
(Eu) vs.
,
where X are
heavy elements other than Eu for the present sample of barium stars and stars from the literature.
From Nd to Y, for which the s-process contribution is larger,
barium stars and post-AGBs are clearly overabundant, located on the high abundance end.
In the Ba panel of Fig. 8, the lowest values of
(Eu) corresponding to normal stars are close to the upper line where
[Ba/Eu] = -0.70, representing Ba production only by r-process (Mashonkina et al. 2003).
As
(Eu) increases, data become closer to the lower line, where data are
compatible with [Ba/Eu] values for the solar ratio. Both barium and post-AGB stars values
are very different from the solar ratio line.
In Fig. 9, the
r-process fraction of Eu abundance [
(Eu))] correlates with
the s-process main component abundance fraction of other elements in barium stars. In
Fig. 10, the r-process fraction of Eu is plotted against fractions
of r-process of other elements, for the data given in Table 2.
In the second, the correlation is linear with no
dispersion, indicating that the s-process main component
contribution causes a scatter in earlier figures. This linear correlation
in Fig. 10 is expected, given that only the r part of the
abundances for all involved elements were used. Strontium was not included
in this figure because it does not have an r-process contribution
(see Table 4).
![]() |
Figure 8:
![]() ![]() |
Open with DEXTER |
Table 7 and Fig. 11 show
[X/Eu] vs. [Fe/H] for barium and post-AGB stars. The star HD 210910 shows lower
values of [X/Eu] for Y, Zr, Ba, La, Ce, Nd, Pr, and Pb, reaching negative values
for [Zr/Eu] = -0.18. The value of [Sr/Eu] is low for the star HD 123396,
[Sr/Eu] = -0.1. This result is expected from
Fig. 13 of Paper I, which shows that the stars HD 123396 and HD 210910
present the lowest values of [SrII/Fe] and [ZrII/Fe], respectively, and their
[Eu/Fe] = 0.50 and 0.54, thus very close. For other stars, [X/Eu] is
in the range of
.
The post-AGBs are also
mainly in this range, the lowest value being [Zr/Eu] = -0.01. For all elements,
[X/Eu] vs. [Fe/H] is approximately constant in the range of metallicities of
the sample stars.
For Mo, which has a contribution of 49.76% of the s-process main component,
26.18% of the r-process and 24.06% of the p-process, [Mo/Eu] are mainly in the range
.
Relative to Ba, the values are lower, in the range
.
Regarding the ratios involving Pr, that has 49% of the contribution from the s-process main
component and 51% from the r-process, data are mainly in the ranges
and
.
There is no r-process contribution in the Pb abundance, according to
Arlandini et al. (1999), while 46% of the contribution comes from
s-process main component and 54% has been attributed in previous work to the
s-process strong component. Data are in the ranges
and
.
Abundances of Pb obtained for barium stars are
usually higher than those of Mo, which are higher than those of Pr, and the same
is true for solar abundances. Only for 5 stars the Mo abundance is
higher than those of Pb and for one star Pr abundance is higher than that of Mo.
However, the [X/Fe] values are such that the values of [Mo/Ba, Eu] are usually lower than
[Pr, Pb/Ba, Eu]. As a consequence, the range for Mo involves lower values in
Fig. 11.
Elements formed mainly by the s-process main component can be divided into two groups:
light s-elements around the magic neutron number 50 and heavy s-elements
around the magic neutron number 82. In this work, Sr, Y, Zr were included in
light s-elements and Ba, La, Ce, Nd in the heavy s-elements groups. For these
elements, the s-process main component contribution is larger than 50% according
to Arlandini et al. (1999). It is
worth emphasising that Sm is rather an r-element because its s-process contribution
is less than 30%, while the r-process contributes with 67.4% of its production. The
presence of light s-elements in very metal-poor stars cannot be
explained entirely by an r-process contribution.
For instance, there is no production of Sr by the r-process, as shown in
Table 4, and it is observed in the most metal-poor stars. Another
nucleosynthetic process related to massive stars is needed in order to explain
such observed abundances. Figure 9, from Burris et al. (2000), shows an increasing trend of
[Sr/Ba] toward lower metallicities
.
It has been suggested
that beyond the observational uncertainties, another
nucleosynthesis source would explain the increasing [Sr/Ba] at lower
metallicities, and this source could be the s-process weak component.
According to Table 4, 14.9% of the solar abundance of Sr comes from
the s-process weak component and 0.56% from the p-process. At metallicities as
low as
,
the Sr production is expected to be low, given that
the s-process weak component is secondary and that in such environment there is
a lack of pre-existing seed nuclei.
Abundance ratios of
light s-elements relative to Ba are shown in Table 7 and
Fig. 11. In the range of metallicities studied, the ratios
are approximately constant in the ranges
,
,
and
.
The AGB stars are included in the same range for [Zr/Ba], while for
[Y/Ba] two of them show higher values, 0.25 and 0.66. Figure 9, from
Burris et al. (2000), shows a dispersion around
with
,
similar to the present sample of
barium stars.
![]() |
Figure 9:
![]() ![]() ![]() ![]() |
Open with DEXTER |
![]() |
Figure 10:
![]() ![]() ![]() ![]() |
Open with DEXTER |
![]() |
Figure 11: Abundance ratios involving s- and r-processes, light and heavy s-elements. Uncertainties indicated are the highest values of Table 7 for each panel. Symbols are the same as in Fig. 7. |
Open with DEXTER |
The s-element nucleosynthesis depends on the neutron exposure to which the seed nuclei were submitted inside the AGB star. Considering the scenario of material transfer for the barium star formation, it is reasonable to expect that the abundance patterns of barium stars show signatures of the neutron exposure during the occurrence of the s-process in the AGB companion.
Models trying to reproduce the resulting abundances of the s-process were presented by Cowley & Downs (1980), who calculated theoretical abundances by using exact solutions from Clayton & Ward (1974) for models of exponential distribution and the approximate solution from Clayton et al. (1961) for a model of single exposure.
Malaney (1987a) presented theoretical predictions of
abundances from the s-process starting on iron seeds, considering single neutron
exposure for several values of
and neutron density
= 10
cm
.
Malaney (1987b) also provides
theoretical predictions, but considering an exponential neutron distribution
for several values of
and two different values of neutron density,
108 cm
and 1012 cm
.
Abundances resulting from theoretical predictions are usually normalised for
(Sr) = 20. Tomkin & Lambert (1983) provide an expression to transform
observational data to this scale
In order to check the fit, Cowley & Downs (1980) used an expression largely
used in the literature (e.g. Pereira & Junqueira 2003; van Winckel & Reyniers 2000). There
![]() |
(17) |
![]() |
(18) |
The best fit between theoretical predictions and observational data corresponds
to the lowest value of ,
where
was calculated for all P values from
Malaney (1987a) and Malaney (1987b) tables. Table 8 shows the best
fittings of
and
for each table in Malaney.
Figures 12 to 15 show the best fittings obtained for each
present sample of barium star, showing the values of
and
neutron density
,
in Cols. 8 and 9 of Table 8.
![]() |
Figure 12: Fittings of observed data to theoretical predictions by Malaney (1987b,a). |
Open with DEXTER |
![]() |
Figure 13: Same as Fig. 12 for 6 other sample stars. |
Open with DEXTER |
![]() |
Figure 14: Same as Fig. 12 for 6 other sample stars. |
Open with DEXTER |
![]() |
Figure 15: Same as Fig. 12 for 8 other sample stars. |
Open with DEXTER |
For some elements of some barium stars, the differences between abundances
of barium and normal stars were too small, although
.
This is the case for
Mo (HD 749), Eu (HD 89948, HD 116869, HD 210709), and Gd (HD 22589).
In these cases, the uncertainty on
is too large, decreasing
the quality of the fit. They were withdrawn before using Eq. (16),
and are easily seen in Figs. 12 to 15 due their discrepant
values with respect to theoretical predictions, and there are no error bars for them.
Except for such cases, there is good agreement between observed and
theoretical data shown in Figs. 12 to 15. For most stars,
the best fit was that for which theoretical predictions consider an exponential
neutron exposure distribution, with neutron density 10
cm
for 5 stars and 10
cm
for 14 stars. For 7 stars, the best fit was
found for a single neutron exposure and
neutron density 10
cm
.
This result is curious since the
exponential neutron exposure distribution is better accepted in the literature.
However, Busso et al. (1999) suggest that s-process is a result of a number of
single exposures instead of an
exponential distribution. Pereira & Junqueira (2003) also obtained a single exposure
and neutron density 10
cm
as the best fit for 2 barium stars.
One of them, HD 87080 for which they found
,
is in common with the
present sample, for which we found a best fit with an exponential distribution of
and neutron density 10
cm
.
The reason for this
difference is unclear. The metallicities adopted were very close (see Tables 9 and 11
from Paper I); however, values of
in this work are usually
higher than those from Pereira & Junqueira (2003). Least-square fittings were used here
to compute abundances of normal stars (Eq. (15)), while
Pereira & Junqueira (2003) use the sum of solar abundance with
star metallicity as a reference.
![]() |
Figure 16:
Example of the fit of parameters G and
![]() ![]() ![]() |
Open with DEXTER |
The uncertainty on
is given by the expression
![]() |
(19) |
![]() |
(20) |
![]() |
(21) |
For comparison purposes, the distribution of abundance in several nucleosynthetic
processes were also done. The total abundance of the element of a
barium star was distributed in solar proportions, with fractions corresponding to
s-, r-, and p-processes, shown in Cols. 3 to 7 in Table 4.
Results for this
curve are shown in Col. 18 of Table 6, and
they are represented by open symbols in Figs. 18 to 21.
These results are lower than others because the overabundance derived from the s-process
is neglected.
A theoretical
curve is calculated with Eq. (1), where
two parameters (
and G) must be determined by fitting observed data.
A robust statistics was used, which, in this work, consists in finding values of G
and
that simultaneously minimise the sum of absolute deviations
represented by
![]() |
Figure 17:
Solar ![]() |
Open with DEXTER |
The method was tested for the Sun by using data from Arlandini et al. (1999). The
resulting line from the fitting is shown in Fig. 17, in which
mb
and G =
%.
Figure 19b from Käppeler et al. (1989) shows their solar
curve,
resulting in
mb
and
%.
Figures 18 to 21 show that the uncertainties,
represented by error bars, are very large
for barium stars, so they were neglected in computing the quality
of the fit. Branching and the too discrepant data were also neglected.
For the present sample, theoretical
curves fit the observed data very well,
as shown in Figs. 18 to 21. This confirms that the solar isotopic
composition is adequate for barium stars and that the transfer of enriched material
keeps approximately the
solar proportions of each nuclide. In order to build the theoretical curve, one
considers that abundances of all elements heavier than iron are null at the
beginning of the neutron capture chain. This supposition becomes less reasonable as the metallicity increases. Furthermore, the fit creates more than one
difficulty: a) there is a lack of elements along the curve due to the
difficulty in finding lines in the spectrum, thereby preventing a
higher quality of the fit between theoretical and observed curves; b) barium star
data show a large dispersion; c) branching values were not considered; d) the numerical
solution of Eq. (1) is difficult, given that two parameters to be
fitted appear in non linear form in the equation.
The good agreement between theoretical
and observed curves is interesting, taking all assumptions
and difficulties found in their building into account.
One way to relate light and heavy s-elements is through hs and ls indices, defined as the mean of the abundances of light and heavy elements, respectively. According to Sect. 4.3, Sr, Y, and Zr were included in ls, and Ba, La, Ce, and Nd were included in hs. For the s index, all these heavy and light s-elements were considered. In the case of missing abundances, such an element was excluded from the index. Figure 22 and Table 9 show [hs/ls], [ls/Fe], [hs/Fe], and [s/Fe] vs. [Fe/H], with [hs/ls] = [hs/Fe] - [ls/Fe]. Uncertainties on [s/Fe], [ls/Fe], and [hs/Fe] have to take into account the contribution of the uncertainties on abundances of each element included in the index.
A behaviour of [hs/ls] as a function of
can be inferred from the slope
of the
curve. The larger
,
the smaller the slope, and,
at the same time, a smaller slope means
larger abundance of heavy s-elements located at the right end of the curve.
The conclusion is that the larger
,
the larger [hs/ls], as
shown in Fig. 17 of Wallerstein et al. (1997).
This is reasonable considering that the chain of formation of s-elements cannot
go far if the neutron flux is low.
![]() |
Figure 22: [hs/ls], [ls/Fe], [hs/Fe], and [s/Fe] vs. [Fe/H]. Filled symbols indicate barium stars, as in Fig. 7; open circles are post-AGBs; crosses are data from Junqueira & Pereira (2001); open squares are data from Luck & Bond (1991); open triangles are data from North et al. (1994). The uncertainties indicated are the higher values shown in Table 9. |
Open with DEXTER |
During the third dredge-up, a given amount of protons is introduced into the
intershell, a region composed of helium located between
helium and hydrogen burning shells. These
protons are captured by 12C and form 13C through the reaction
12C(p, )13N(
)13C or 14N through
13C(p,
)14N, creating the 13C pocket.
If the thermal pulse is independent of metallicity, protons are
introduced into the intershell in the same amounts for higher or lower metallicities.
It means that the neutron source
C(
, n)
O,
probably the main neutron source in AGB stars, is independent of metallicity.
However, the neutron number by seed iron nucleus
will be larger at low metallicities (Clayton 1988). If
is
proportional to [hs/ls] and inversely proportional to [Fe/H], [hs/ls]
will also be inversely proportional to [Fe/H], so it is expected that
the lower the metallicity, the higher the ratio [hs/ls]. Figure 24
confirms a correlation between
and [hs/ls] for barium
and post-AGB stars. However, the anticorrelation between [hs/ls] and [Fe/H]
is weak for barium and post AGBs stars, as shown in Fig. 22.
In the same way, the anticorrelation between
and [Fe/H] is not
confirmed in Fig. 24.
The yield of all s-elements to decrease with decreasing metallicities is
compatible with the secondary characteristic of the s-process, which requires
pre-existing seed nuclei.
At intermediate metallicities (
-0.8), the Ba peak is dominant among
s-process products in AGB models (Busso et al. 1999). For higher metallicities the
Zr peak dominates.
If giant, as well as dwarf, barium stars have the same physical origin for the
accretion of enriched material from a more evolved companion, it is reasonable to
expect that if the neutron exposure is higher for dwarf and more metal-deficient stars,
the same occurs for giants.
In Fig. 24,
derived from the
curve are mainly
within 1.8 <
< 0.6, with less spread as compared to those derived
from theoretical predictions by Malaney (1987b,a), although
both show the same trend. For giant stars,
was found from
theoretical predictions, but not from the
curves.
Solid and dashed lines in Fig. 24 are least-square fittings, with parameters:
Table 9: s-process indices: s, ls, and hs for the present sample of barium stars.
Table 10: s-process indices s, ls, and hs for barium stars collected in the literature. References: N94 - North et al. (1994); LB91 - Luck & Bond (1991).
Figure 10 from van Winckel & Reyniers (2000) shows [hs/ls] vs. [Fe/H], including data from several previous work. Despite the dispersion, [hs/ls] increases toward higher [Fe/H]. The dependence of neutron exposure on [Fe/H] could indicate that other important parameters affect dredge-up events and nucleosynthetic processes along the red giant branch evolution.
![]() |
Figure 23: [hs/ls] vs. [s/Fe]. Symbols are the same as in Fig. 7. Error bars in the bottom right corner show the largest uncertainties. |
Open with DEXTER |
![]() |
Figure 24:
[hs/ls] vs.
![]() ![]() ![]() ![]() ![]() ![]() ![]() |
Open with DEXTER |
For the present sample of barium stars, [s, ls, hs/Fe] and [hs/ls] decrease slightly
toward higher metallicities in the range
0.45 [s/Fe]
1.6,
0.4
[ls/Fe]
1.6,
0.7
[hs/Fe]
1.75,
and -0.2
[hs/ls]
0.7.
According to Wheeler et al. (1989), [s/Fe] = 0 for less evolved stars
in this same range of metallicities; therefore, the results found for
the present sample show the overabundance of s-elements in barium stars.
One has to be careful in comparing the indices of the present work to those of post-AGBs in Figs. 11, 22, and 23, given that Reyniers et al. (2004) and van Winckel & Reyniers (2000) included Sm and not Ce in s and hs, and only Y and Zr in ls. Figure 23 shows that [hs/ls] vs. [s/Fe] has a large dispersion for barium stars, whereas there is a linear correlation for post-AGBs. Figure 8 from Reyniers et al. shows the same correlation but using different combinations of the elements included in the indices. One of them uses Ce, but the linear correlation is still present comparing with using Ba, La, Nd, and Sm in hs. Other configurations show larger dispersion, hence the presence of Ce in hs and s is not the origin of the dispersion for barium stars.
According to Table 7 and Fig. 11, there is no difference as a function of log g, confirming that the overabundances characteristic of a barium star do not depend on luminosity classes.
Abundances of elements with a lower contribution from the s-process main component are closer to the normal stars, whereas elements with a higher s-contribution are more overabundant. However, r-elements such as Sm, Eu, Gd, and Dy are also enriched in barium stars, sometimes at similar magnitudes to s-elements, and this occurs because the s-process main component chain includes the r-elements with A < 209.
The ratios involving light s-elements (Sr, Y, and Zr) and Ba, a heavy s-element, are approximately constant in the range of metallicities of the present sample, with a dispersion similar to normal stars and post-AGBs.
Considering that the abundance fraction due to all processes, except the
s-process main component, was present in the proto-barium star, it was possible
to isolate the fraction
corresponding to the s-process main component for barium stars.
The solar isotopic distribution of this abundance fraction was used to
build the observed
curves. Their fittings to theoretical curves indicate
that it is possible to use a solar isotopic mix to estimate the
contribution of the s-process main component for each isotope of an element.
The derivation of
was obtained by fitting observed data to theoretical
predictions and
curves. Theoretical predictions for
abundances starting with Sr fit the observed data very well for the
present sample, thereby providing an estimation for neutron exposure that occurred in
AGB supplying the s-process. The [hs/ls] vs.
considering these theoretical
predictions increases linearly, whereas no conclusion was reached
for
vs. [Fe/H].
The
values from
fittings also provide linear and increasing
vs. [hs/ls], but data in this case are
in the range 0.2
0.8, showing
lower dispersion than those for theoretical predictions.
Abundances obtained for barium stars are close to those for AGBs stars, but they are usually lower. This is reasonable considering that only part of the surface material of AGBs is transferred to the companion, which becomes a barium star.
Acknowledgements
We acknowledge partial financial support from the Brazilian Agencies CNPq and FAPESP. DMA acknowledges a FAPESP PhD fellowship No. 00/10405-8 and a FAPERJ post-doctoral fellowship No. 152.680/2004. We are grateful to Marcelo Porto Allen for making his robust statistics code available, and to the referee, Nils Ryde, for useful comments.
![]() |
Figure 18:
The ![]() |
Open with DEXTER |
![]() |
Figure 19: Same as Fig. 18 for 6 other sample stars. |
Open with DEXTER |
![]() |
Figure 20: Same as Fig. 18 for 6 other sample stars. |
Open with DEXTER |
![]() |
Figure 21: Same as Fig. 18 for 8 other sample stars. |
Open with DEXTER |
Table 4:
Cross-sections in mb
for 30, 23, and 8 keV and abundance fractions
for s-, r-, and p-processes relative to total abundance and to total abundance of each process.
Table 5:
Abundance fractions corresponding to the s-process main component for
barium stars calculated with Eq. (8) (upper table) and percentage of abundance
(lower table) relative to the s-process main component of a barium star compared to normal stars
of a same metallicity (
).
Table 7: [X/Eu] ( upper table) and [X/Ba] ( lower table) for the barium stars.
Table 8:
Results of neutron exposures. Columns 2 and 3 correspond to results of the
fittings to the theoretical predictions of Malaney (1987a) and Cols. 4 to 7 to those of Malaney (1987b).
(u),
(e8),
(e12) are the fit quality, respectively, for single exposure,
exponential exposure under neutron density 10
cm
,
and exponential exposure under
neutron density 10
cm
.
,
,
and
are the mean
distributions of neutron exposure corresponding to
(u),
(e8), and
(e12),
respectively.
(f) and
are results corresponding to smaller
.
Column 10
shows if the best fit corresponds to exponential (exp) or single (sing) exposure. Columns 11 to 13
correspond to the results of the fit of
curves, with
as the quality of the fit.
Numbers in parenthesis are errors in last decimals and were estimated directly from the fittings.
The uncertainties on
,
,
and
were estimated in
0.1 mb
.
For all cases,
is given in mb
.