A&A 467, 665-677 (2007)
DOI: 10.1051/0004-6361:20065797
S. Feltzing1 - M. Fohlman1 - T. Bensby2
1 - Lund Observatory, Box 43, 221 00 Lund, Sweden
2 -
Department of Astronomy,
University of Michigan,
Ann Arbor, MI,
USA
Received 10 June 2006 / Accepted 18 December 2006
Abstract
Context. Manganese is an iron-peak element and although the nucleosynthesis path that leads to its formation is fairly well understood, it remains unclear which objects, SN II and/or SN Ia, that contribute the majority of Mn to the interstellar medium. It also remains unclear to which extent the supernovae Mn yields depend on the metallicity of the progenitor star or not.
Aims. By using a well studied and well defined sample of 95 dwarf stars we aim at further constraining the formation site(s) of Mn.
Methods. We derive Mn abundances through spectral synthesis of four Mn I lines at 539.4, 549.2, 601.3, and 601.6 nm. Stellar parameters and data for oxygen are taken from Bensby et al. (2003, 2004, 2005).
Results. When comparing our Mn abundances with O abundances for the same stars we find that the abundance trends in the stars with kinematics typical of the thick disk can be explained by metallicity dependent yields from SN II. We go on and combine our data for dwarf stars in the disks with data for dwarf and giant stars in the metal-poor thick disk and halo from the literature. We find that dwarf and giant stars show the same trends, which indicates that neither non-LTE nor evolutionary effects are a major concern for Mn. Furthermore, the [Mn/O] vs. [O/H] trend in the halo is flat.
Conclusions. We conclude that the simplest interpretation of our data is that Mn is most likely produced in SN II and that the Mn yields for such SNae must be metallicity dependent. Contribution from SN Ia in the metal-rich thin disk can not, however, be excluded.
Key words: stars: abundances - Galaxy: abundances - stars: late-type
Manganese is an element in the lower iron group. It has only one stable isotope, 55Mn, which is thought to be mainly produced in the processes of explosive silicon burning and nuclear statistical burning (Woosley & Weaver 1995). Even though the nuclear path that leads to 55Mn is fairly well understood it remains unclear which objects that are the main contributors to the chemical enrichment for this element and if the Mn yields from either supernovae type Ia (SN Ia) or supernovae type II (SN II) are metallicity dependent (e.g. McWilliam et al. 2003; Carretta et al. 2004; Woosley & Weaver 1995).
Gratton (1989) investigated metal-poor stars and showed that the run
of [Mn/Fe] vs. [Fe/H] mimics that of the
-elements but
"up-side-down'', i.e. that [Mn/Fe] is under-abundant in
low-metallicity stars but that the trend subsequently starts to
increase towards the solar value. In the study by Gratton (1989) this
increase starts at [Fe/H] = -1. Based on these observations he
suggested that Mn is made in SN Ia to a much larger degree than in
SN II, hence this leads to an increase at higher metallicities when
the SN Ia starts to contribute to the chemical enrichment of the
interstellar medium from which the subsequent generations of stars are
formed.
However, nucleosynthesis calculations (e.g. Arnett 1971; Woosley & Weaver 1995; Chieffi & Limongi 2004; Limongi & Chieffi 2005) show that the yields from SN II likely are metallicity dependent. In such a case the trends observed by e.g. Gratton (1989) could be explained by higher yields from SN II with metal-rich progenitor stars.
McWilliam et al. (2003) considered Mn abundances in three different stellar populations: a sample of Galactic bulge K giants, giant stars belonging to the Sagittarius stream, and a sample of solar neighbourhood stars. The latter was taken from Nissen et al. (2000) but with improved treatment of hyperfine structure (Prochaska & McWilliam 2000). From their comparison of these three stellar populations McWilliam et al. (2003) concluded that Mn is produced both in SN Ia as well as in SN II and the Mn yields from both types of SNe are metallicity dependent.
Carretta et al. (2004) revisited the arguments of McWilliam et al. (2003) and argued that by including two red giant stars in the Sagittarius (Sgr) dwarf spheroidal galaxy, studied by Bonifacio et al. (2000), the case for metallicity dependent SN Ia yields was considerably weakened.
In this paper we study two local stellar samples, representative of the thin and the thick disk, respectively, that have been shown to have different elemental abundance trends for a number of elements (Bensby et al. 2003, 2005; Bensby & Feltzing 2006).
Of particular interest for our study is that the trends of oxygen abundances for these stars are well separated and tight (Bensby et al. 2004). This means, since oxygen is only produced in SN II, that we have access to a "clock'' that is independent of the SN Ia time-scale. Combining this information with the new Mn abundances and data for halo stars and stars in dwarf spheroidal galaxies (dSph) we investigate the origin of Mn.
The paper is organized as follows. In Sect. 2 the two stellar samples are defined and the observations described. Section 3 describes the abundance analysis in detail including an extended discussion of error sources. After that Sect. 4 describes the resulting abundance trends and these are discussed and interpreted in Sect. 5. Finally, Sect. 6 summarizes our findings.
Table 1: Stellar parameters for the stars.
The stars were selected purely based on their kinematic properties. Briefly, we assume that the velocity components for the halo, thick disk, and thin disk all have Gaussian distributions. Allowing for the different asymmetric drifts of the three populations we calculate the probability for each star that it belongs to the halo, the thin disk, and the thick disk. For the thick disk sample we then selected those stars that were much more likely to be thick than thin disk stars and vice versa. In this way we selected two fairly extreme kinematic samples of disk stars. The full selection procedure is detailed in Bensby et al. (2003, 2005).
The stars were observed during four observing runs. Two using the
FEROS spectrograph on the ESO 1.5 m telescope (in Sept. 2000 and
Sept. 2001) and two runs using the SOFIN spectrograph on the NOT
telescope (Aug. 2002 and Oct./Nov. 2002). The spectra from FEROS have
a resolution of
and those from SOFIN have
.
For details about data reduction see Bensby et al. (2003, 2005).
At each observing run we also obtained spectra of scattered solar light. At the FEROS runs we obtained an integrated solar spectrum by observing the sky in the afternoon. In the SOFIN runs we obtained a solar spectrum by observing the Moon during the night.
Table 1 lists the stars used in this study, their fundamental parameters, and which spectrograph that was used to obtain the stellar spectra.
We have performed a standard Local Thermodynamic Equilibrium (LTE) analysis to derive chemical abundances with the help of synthetic spectra and one-dimensional, plane-parallel model atmospheres. We have used the MARCS model atmosphere code (Gustafsson et al. 1974; Edvardsson et al. 1993; Asplund et al. 1997) to generate the stellar atmospheres used in the abundance analyses.
Mn is an odd-Z element and hence its lines are subject to hyperfine splitting (hfs) and at least for strong lines this effect must be taken into account. As our stars span a large range in [Fe/H] and as we know from previous studies (e.g. Gratton 1989) that Mn abundances increases as [Fe/H] increases the Mn I lines that are weak at low [Fe/H] will be rather strong at high [Fe/H]. Hence, in order to derive reliable abundance trends for our two samples of stars we must include the hfs in the analysis of our spectra.
Stellar parameters are taken from Bensby et al. (2003, 2005). The
effective temperatures (
)
are based on excitation
equilibrium whilst the
values are derived using Hipparcos
parallaxes. Full details of the iterative procedure to tune the
stellar parameters can be found in Bensby et al. (2003).
We selected Mn I lines for our analysis based on two criteria: the lines should be easily identified and analyzed in our FEROS spectra and they should have data for the hfs readily available in the literature. Based on this we selected 5 lines from the list of Prochaska et al. (2000). Eventually one line (at 602.1 nm) was discarded (see Sect. 3.4 in paragraph "Erroneous hyper fine structure (hfs)''). The remaining four lines are listed in Table 2.
Atomic data for other lines in the vicinity of the Mn I lines that are needed for the line-synthesis were extracted from the Vienna Atomic Line Database (VALD) (Kupka et al. 1999; Ryabchikova et al. 1999; Piskunov et al. 1995).
The
-values used in our study are the same as in
Prochaska et al. (2000), i.e. from Booth et al. (1984). These
-values
have been experimentally determined and have a relative accuracy (as
quoted by Booth et al. 1984) of better than 2%. Furthermore, Booth
et al. (1984) estimate that the absolute scale for their values has
an accuracy of about 7% for the ground state lines and 3% for
the excited state lines. Although more recent measurements have been
done for some Mn I lines that are observable in stellar spectra
no laboratory work has been done on the lines used in this study. We
thus have no possibility for an independent, empirical assessment of
the accuracy of the
-values.
We have, however, made comparisons of the solar abundances derived
from the lines from multiplet 1 and 27 separately. The results are
given in Eqs. (1)-(3). For the lines in multiplet 1
the differences are small, but for the lines in multiplet 27 we see a
real offset (Eq. (2)). To first order we attribute this to an
erroneous
-value (see discussion in Sect. 3.4.1,
paragraph "Absolute abundances'').
![]() |
Figure 1:
Upper panel shows the observed spectrum, |
| Open with DEXTER | |
Our calculations of the stellar spectra take into account the broadening of the stellar lines caused by collisions with hydrogen atoms (van der Waals damping). Whenever possible, we use the data from Barklem et al. (2000) for these calculations. When data are not available from Barklem et al. (2000), we apply an enhancement factor to the classical Unsöld approximation of the van der Waals damping, which for most elements were set to 2.5, following Mäckle et al. (1975). For Fe I we take the correction terms from Simmons & Blackwell (1982), but for Fe I lines with a lower excitation potential greater than 2.6 eV we follow Chen et al. (2000) and adopt a value of 1.4. For Fe II we adopt a constant value of 2.5 (Holweger et al. 1990).
Two of our Mn I lines have data in Barklem et al. (2000); the lines at 543.2 and 539.4 nm. Experimentation shows, however, that it does not matter for the abundance analysis which values (i.e. the classical Unsöld approximation or Barklem et al. data) we use for the synthesis as these two lines are, in the Sun, completely insensitive to any reasonable (and unreasonable) changes of the enhancement factor or using the Barklem et al. (2000) data.
![]() |
Figure 2:
Demonstration of the determination of the Mn abundance using
line synthesis. Here for the Mn I line at 539.4 nm in
HIP 91438. The star has [Fe/H] =-0.24 and the spectrum a
|
| Open with DEXTER | |
Table 2:
Selected Mn I lines used in the abundance analysis
(see Sects. 3.2 and 3.3). In the
first column we give the wavelength, in the second the excitation
energy of the lower level in the atom involved in the line
formation, in column three the
-values from Prochaska et al. (2000), column four lists the multiplet the line belongs to,
and column five lists the Fe I lines used to infer the
combined macro turbulence and rotation,
(Sect. 3.3 and Fig. 1).
For the two lines in the red the situation is different. They are both
very sensitive to changes in the enhancement factor. However, these
lines are not included in the calculations by Barklem et al. (2000). We can therefore only conclude that the lines are
sensitive to changes in the damping constant. However, we feel that
the value we use does give reasonable results in that the
(i.e. the line-to-line scatter) does for example not increase as
[Fe/H] increases (see Sect. 3.4 paragraph Line-to-line
scatter).
In hot stars the rotational broadening often dominates over the macro turbulence but in cool stars the two broadening effects are comparable in size (Gray 1992). Therefore the two effects can be modeled simultaneously. To determine the rotational velocity and the macro turbulence for the stars un-blended Fe-lines in the vicinity of our Mn I lines were synthesized. Data for the Fe I lines and adjacent lines were extracted from the VALD database (compare Sect. 3.2). Wavelengths for these Fe I lines are listed in Table 2.
The synthetic spectra were first convolved with a Gaussian profile in
order to imitate the instrumental broadening and secondly by a
radial-tangential (rad-tan) profile to simulate the combined effects
from macro turbulence and rotation (Gray 1992). Before the combined
rotational and macro turbulence velocity (
)
could be
determined the
-values for the Fe I lines had to be
adjusted in order to achieve a good fit to the observed spectrum.
The fit was deemed good when the intersection points between the
observed spectra and the synthetic spectra lied at zero and the
difference between synthetic spectra and observed spectra, where the
wings reach the continuum, also were at zero (see
Fig. 1 and compare e.g. Bensby et al. 2004).
To derive the
we calculated seven synthetic spectra
with different
(in steps of
)
for each Fe I line. The synthetic spectra were
then compared with the observed spectrum and an un-normalised (![]()
)
was calculated for each comparison (i.e. observed-synthetic) and the
was read off as the value that
minimizes ![]()
.
The procedure is illustrated in
Fig. 1.
![]() |
Figure 3:
Mn I lines used in the final abundance analysis. Here
we show the solar spectrum with |
| Open with DEXTER | |
The synthetic spectra for the Mn I lines were first convolved
with a Gaussian line profile to compensate for the instrumental
broadening and after that they were convolved with a rad-tan profile
to account for the broadening due to the combined effects of rotation
and macro turbulence (see Sect. 3.3.1). After the
continuum was set a ![]()
was calculated in the same way as
when deriving the
and the ![]()
was minimized to
estimate the abundance. An example for how the Mn abundance was
derived for the Mn I line at 539.4 nm is shown in
Fig. 2.
Finally, all the derived abundances were normalized, line by line, to the solar values that we obtained from our own spectra of scattered solar light (Sect. 2). This means that the Sun will always have [Mn/H] =0. It is the mean value of these normalized values that we are using in the Sects. 4 and 5 and which are the basis for the [Mn/H] values quoted in Table 1.
Table 3:
Effects on derived Mn abundances due to errors in the stellar
parameters. For four stars we show the effect on the resulting abundance,
line by line (columns two to five),
when
is increased by 70 K, when
is increased
by 0.1 dex, and when [Fe/H] is increased by 0.1 dex.
In this work we are concerned with a differential study. This means
that the Mn abundance derived from each Mn I line is normalized
to the abundance we derived from our own solar spectra for that
particular line (see Sect. 3.3.2). Furthermore, we are
studying stars that span a fairly narrow range in
as
well as in
(see e.g. Table 1). This means
that when we compare two stars at a given [Fe/H] then any differences
in the derived Mn are, to first order, independent of the assumptions
made for the stellar parameters and in the modeling of the line.
There are two types of errors that we in general are concerned with - statistical and systematic. We will here quantify both of them, however, bearing in mind that for the interpretation of our results in Sect. 5 it is mainly the statistical errors that concerns us.
In Fig. 4 we show the
([Mn/H]) arising from the
line-to-line scatter for each star using the values normalized to the
solar values (see Sect. 3.2). As can be seen the
resulting
are small, generally below 0.06 dex with a mean
error for the FEROS observations of 0.047. This gives a formal error
in the mean
of 0.024 dex. For
stars observed with SOFIN only two lines are used and hence the
are of less value to estimate the internal
consistency. However, for most stars observed with SOFIN the two lines
are in good agreement. Furthermore, there is no sign
that the scatter increases with increasing Mn abundance (and hence
increasing Fe abundance). In summary we find that the internal
statistical errors due to line-to-line scatter are small and an easily
quantifiable source of scatter in any abundance trend.
We have used an un-normalised
to estimate the final
abundances rather than a
-by-eye procedure. Our procedure would
give a best value also when e.g. the hfs in the lines are not well
reproduced in the line-synthesis (i.e. due to incomplete or erroneous
atomic data). However, since we always do additional visual
inspections of the fits we believe that for the data-set presented
in this work this type of systematic error is rather small which is
evidenced in e.g. Figs. 2 and 4. We have
also extensively checked our derived abundances, line-by-line, as a
function of the stellar parameters and found no obvious trends (see
Sect. 3.4.1). This together with the fact that the
internal scatter between Mn abundances derived from different lines is
small (see paragraph above and Fig. 4) implies that the
hfs, as well as other aspects of the shape of the lines have been
taken care of.
![]() |
Figure 4:
|
| Open with DEXTER | |
The exact placement of the continuum is sometimes difficult. Especially so when there are many lines close to the line that is being analyzed which results in an overall depression of the apparent continuum. Figure 3 shows the spectrum in the vicinity of each of the four Mn I lines used in this study. As can be seen these, small, regions of the spectra would (in most cases) not be enough to judge the continuum level wherefore we relied on a larger wavelength range to set the continua. In general the continua are harder to set in the most metal-rich stars as the lines are getting stronger in such stars. However, this appears to be of no great concern in our study as the line-to-line scatter does in fact not increase when [Mn/H], and thus also [Fe/H], increases. We note that from the SOFIN spectra of the Sun, for lines at 539.4 and 543.2 nm, we derive somewhat lower Mn abundances than from the FEROS spectra (see Table 3). This could be due to that we placed the continua too low in the SOFIN spectra. However, we have not been able to quantify this. As we normalize our final abundances to that of our own solar abundance, line-by-line, this will not be a concern for our final conclusions.
Prochaska & McWilliam (2000) showed that the hfs data used by Nissen
et al. (2000) gave erroneous abundances. Figure 2 in Prochaska &
McWilliam (2000) shows the effect of applying the correct hfs to the
data by Nissen et al. (2000). In particular they found that [Mn/Fe]
increased by
0.05 dex and that the effect was larger below
.
Recently, del Peloso et al. (2005) have compared
different sets of hfs. In particular they compare the hfs from
R. Kurucz's electronically available linelists
and
Steffen (1985). Nissen et al. (2000) used the data by
Steffen (1985). In good agreement with Prochaska &
McWilliam (2000), del Peloso et al. (2005) find that the two sets
of hfs result in Mn abundances that show increasingly different
[Mn/H] as [Fe/H] decreases and that the differences increases
progressively from
in the sense that the data from
Steffen (1985) give the smaller [Mn/H]. As discussed in
Sect. 3.2 we use the hfs (both as concerns wavelength
as well as the strength of the different hfs components) listed by
Prochaska et al. (2000) in their Table 19
. In our analysis we have payed special attention to
how well the structure of each Mn I line is fitted. Through
inspections of the fitting of the line at 602.1 nm in
15 stars
with a large span in [Fe/H] we came to the conclusion that the hfs for
this particular line as listed in Prochaska et al. (2000) (which is
based on R. Kurucz's electronically available hfs lists) is not a
good representation of the actual shape of the line. The miss-match
was particularly obvious in stars with small values for the combined
rotation and macro turbulence broadening, i.e. where lines are
sharper. This line was hence discarded from further analysis. For the
other lines we found no such discrepancies. We thus believe that, at
the accuracy of our observed spectra, the hfs we use is a good
representation of the actual line profiles.
![]() |
Figure 5:
Panels a)-c) show the differences between the Mn
abundance derived from the Mn I lines at 539.4 and 543.2 nm
as a function of
|
| Open with DEXTER | |
Our current analysis does not include a full non-LTE analysis of the line formation for the Mn I lines. However, in order to exclude that our results should be compromised by this omission we have done extensive checks to exclude that departures from LTE is a major concern for the stars that we study.
The lines we kept for our final analysis divide into two pairs. One set of lines with a lower excitation energy of 0.0 eV and one with a lower excitation energy of 3.07 eV (see Table 2). We started by investigating these two line pairs separately.
The differences in derived Mn abundance from the lines at 539.4 nm
and 549.2 nm show no trends with
,
,
or
[Fe/H]. Furthermore, the differences between the two lines are small
with a low scatter (see Eqs. (1) and (3)) (Fig. 5a-c).
For the line pair at
601.3 and 601.6 nm the comparison,
again, showed no obvious trends with either
,
,
or [Fe/H] (Fig. 5 d-e, note that there are no
SOFIN observations for these lines). However, there is a significant
offset between the Mn abundance derived from the two lines
(Eq. (2)). Both of these two lines also show offset with
respect to the abundances derived from the lines at 539.4 and
543.2 nm. The abundance derived from the line at 601.3 nm is offset
by -0.074 dex from the mean Mn abundance derived from the two lines
at 539.4 and 543.2 nm, while the abundance derived from the line at
601.6 nm is offset by +0.059 dex. The standard deviation is in both
cases 0.06 dex.
We find no trends with either [Fe/H] or
for the difference
between the mean of the 539.4 nm and 543.2 nm and the mean of the
601.3 nm and 601.6 nm lines (see Figs. 6b. and c).
When we do the same comparison but now as a function of
we do find what might be a pattern, Fig. 6a., such
that there appears to be a "dip'' in the difference at
5400 K,
spanning perhaps
250 K. For stars hotter than 5750 K there
appears to be hardly any trend. Also for stars with temperatures below
5100 K the trend appears flat as well, although here we are
dealing with only a few stars so the conclusion will be less robust.
Note, that also when we compare the two lines in multiplet 27 one by one
to the mean value derived from the two lines of multiplet 1 the same
pattern shows up.
These comparisons show that the excitation balance, most probably, is not influenced by departure from LTE. However, departures from ionizational balance may still be present.
For our observation of the Sun with SOFIN we derive
and
for the lines at 539.4 and 543.2 nm,
respectively. The internal consistency is also here good, in fact
better than for the FEROS spectra. However, the absolute values
deviate more from the standard value than the abundances based on
the FEROS spectra.
For the lines at 601.3 nm and 601.6 nm (multiplet 27)
we derive
(Mn) =
and
dex, respectively. Again, the value derived from the line
at 601.6 nm agrees very well with the standard value. Could the line
at 601.3 nm suffer from departures from LTE? If this was the case it
would in general be likely that the line that suffered from non-LTE
effects should show changes as a function of e.g.
or
another atmosphere parameter. In this case we see no such thing, only
an offset from the other line in the same multiplet and, to first
order, one would attribute the offset to an erroneous loggf-value
for the line at 601.3 nm.
In summary, it appears that the three lines at 539.4 nm, 543.2 nm,
and 601.6 nm all yield Mn abundances that are, within the errors, in
accordance with the standard value for the solar photosphere, whilst
the line at 601.3 nm gives
(Mn) about 0.15 dex lower
than the other lines. A first interpretation of this deviation is that
the offset is due to the errors in the derivation of the
loggf-value.
However, the size of the changes in equivalent width are small. Moreover, although we find a difference between the abundances derived for the Sun from the spectra recorded by the FEROS and SOFIN spectrographs, which are taken about two years apart, once the data for all stars have been normalized to the solar value (i.e. putting the Sun at 0 dex) there is no difference between the data recorded with the two different spectrographs. This indicates that the size of the variations induced by the chromospheric activity, if real, are sufficiently small that it does not effect elemental abundance analysis at the level of resolution and S/N that our data have.
![]() |
Figure 6:
Differences between the mean Mn abundances derived from the
Mn I lines at 539.4 and 543.2 nm and the two lines at 601.3
and 601.6 nm as a function of
|
| Open with DEXTER | |
Table 4: Mn abundances line by line for stars observed with FEROS. The wavelengths of the lines used are indicated in the header row. The last column indicates if the star is classified as thin or thick disk star (see Bensby et al. 2003, 2005).
In summary we find that the statistical errors in our derived abundances are small. The line-to-line scatter is less than 0.06 dex and the mean error in the mean is 0.047 dex. Statistical errors in the stellar parameters also give low errors, on the same order as the line-to-line scatter.
Systematic errors are more difficult to assess. We find one
-value that might be erroneous. Further measurements of the
-values for these lines are desirable.
Errors arising from the limitations in the modeling include e.g. non-LTE effects and neglect of chromospheric activity. Both are possibilities for the lines we use. From an empirical inspection of our results we find no basis for claiming that such effects are present in the spectra we have analysed.
Table 5: Mn abundances line by line for stars observed with SOFIN. The wavelengths of the lines used are indicated in the header row. The last column indicates if the star is classified as thin or thick disk star (see Bensby et al. 2003, 2005).
![]() |
Figure 7:
a) [Mn/Fe] vs. [Fe/H] trends for the sample with
kinematics typical of the thin disk. |
| Open with DEXTER | |
![]() |
Figure 8: [Mn/O] as a function of [O/H]. The oxygen data are taken from Bensby et al. (2004) and Bensby et al. (2005). Stars with kinematics typical of the thick disk are shown as filled circles and the stars with kinematics typical of the thin disk are shown as open symbols. Typical error-bars are shown in the lower right hand corner. |
| Open with DEXTER | |
The abundance results for all our stars are summarized in Tables 1, 4, and 5 and shown in Fig. 7. For stars with kinematics typical of the thick disk we see a steady increase in the [Mn/Fe] as a function of [Fe/H], while for the stars with kinematics typical of the thin disk we see an essentially flat trend below [Fe/H] = 0 and a shallow increase at higher metallicities.
We also combine our new Mn abundances with oxygen abundances from Bensby et al. (2004, 2005). In Fig. 8 the trend of [Mn/O] vs. [O/H] is shown. For the stars with kinematics typical of the thin disk we simply see a steady increase in [Mn/O] as [O/H] increases.
Stars with kinematics typical of the thick disk show a similar trend, albeit offset, for [O/H] < 0. For higher [O/H] there is a hint of a faster increase in [Mn/O]. However, based on the current data set this must remain a tentative observation.
In summary we find that when plotting [Mn/O] vs. [O/H] the two samples show trends with similar slopes but an offset of about 0.3 dex in [Mn/O].
![]() |
Figure 9:
Mn and O abundances for several samples collected from the
literature. a) Our own samples (compare
Fig. 8). |
| Open with DEXTER | |
That Mn is made in explosive nucleosynthesis appears well established through theoretical modeling, e.g. Arnett (1996). However, it is still uncertain how much Mn is expelled from the SNe, i.e. the yield, and it also remains uncertain which type of SN are the main contributor to the Mn enrichment, i.e. core collapse or degenerate systems. It is also controversial how metallicity dependent the yields may or may not be (see e.g. McWilliam et al. 2003 and Carretta et al. 2004 for a discussion of the observational evidence).
McWilliam et al. (2003) extensively discussed the available data for Mn and concluded that it is most likely made in SN II and that the yields are metallicity dependent. Furthermore, from the data for the Sgr dSph they found preliminary evidence that Mn may also be made in SN Ia. Carretta et al. (2004) used Mn as an example to explore the origins of elements with odd atomic numbers. By including two Sgr dSph stars analysed by Bonifacio et al. (2000) they concluded that the differences between the galactic field stars discussed in McWilliam et al. (2003) and their Sgr dSph RGB stars were less compelling and thus that the case for significant production of Mn in SN Ia was weakened. Our new data for two kinematically selected samples significantly increases the information available for Mn at disk-like metallicities. We will therefore revisit the data discussed in McWilliam et al. (2003) and Carretta et al. (2004) and combine that with our new data as well as data for halo stars in the Milky Way to infer the origin of Mn.
McWilliam et al. (2003) revisited the Mn data from Nissen et al. (2000). In particular they applied a correction to the Nissen et al. (2000) Mn values to correct for an erroneous hyper fine structure
that had been used in that study. In McWilliam et al. (2003) they
noted that for the disk sample they saw a jump in [Mn/Fe]. This jump
appear at [Fe/H]
.
McWilliam et al. (2003) question if this
apparent jump is real. Our data shows that the jump can be explained
as a result of the selection process in Nissen et al. (2000). Our
study consists of two stellar samples: one with kinematics that are
typical for the thick disk and one which has kinematics typical of the
thin disk in our Galaxy. Our stars are confined to the metallicity
range which is typical of the Galactic disks, i.e. with an [Fe/H] of
-1 to +0.5 dex. The sample of stars in Nissen et al. (2000) are
from the list by Chen et al. (2000). In Chen et al. (2000) the stars
were selected to represent the thin and the thick disks, respectively.
However, the criteria they used are such that they selected
exclusively thin disk stars for
and only thick disk
stars for
(see also discussion in Bensby et al. 2003
and their Fig. 14). If we from our sample,
Fig. 7a, take only thick disk stars for
and only thin disk stars for
then the
resulting abundance trend would be very similar to the one in Nissen
et al. (2000) as seen in the corrected abundances from Prochaska &
McWilliam (2000). We thus conclude that the jump in [Mn/Fe] vs. [Fe/H]
found in previous studies is due to a sample selection effect and not
an intrinsic property of the disks. It is, however, worth noting that
the Nissen et al. (2000) indeed is the first study that indicated
that the Mn trends for the thin and the thick disks do differ. A
finding that we confirm and extend to a much larger [Fe/H] range, in
particular for the thick disk.
![]() |
Figure 10:
Toomre diagram showing the kinematics for the stars with Mn
abundances from Sobeck et al. (2006). Stars with kinematics typical
of the thin disk are shown as |
| Open with DEXTER | |
An
-element, like oxygen, is a better reference element
than Fe if we want to investigate which supernovae an element (mainly)
comes from. The main reason for this is that oxygen is only made in
SN II whilst Fe is made in both SN II and in SN Ia (e.g. Timmes et al. 1995). With O as the reference element the interpretation of any
abundance trend hence becomes simpler. Oxygen is preferred over
-elements such as Ca and Mg since those elements have a small
contribution from SN Ia (Timmes et al. 1995).
The main reason why previous studies have not used O extensively as the reference element is because it is difficult to derive O abundances from stellar spectra for three reasons: the O lines are few and hence not always present in the available stellar spectra; the stronger O lines (the triplet at 777 nm) have significant non-LTE effects (see e.g. Kiselman 1993); and the lines that do not suffer from non-LTE effects are very weak, in fact so weak that they may not be observable at all in more metal-poor stars.
Over the last years an extensive discussion has taken place in the
literature about the exact [O/Fe] vs [Fe/H] trend at metallicities
below [Fe/H]
(see for example volume 45 issue 8 of New
Astronomy and discussions in Nissen et al. 2002 and Asplund et al. 2004). However, with the advent of larger telescopes and a
better understanding of how to compensate for the non-LTE effects that
are present in the triplet lines the number of stars with reliable O
abundances have increased enormously.
In order to study the origin of Mn, i.e. which type of supernovae are
the main contributors and if the yields are metallicity dependent, we
have collected a representative set of studies of stars with Mn as
well as O data available in the literature. In particular we have
selected 42 dwarf and giant stars from the new study of Mn by Sobeck
et al. (2006) [Mn data kindly provided by the authour prior to
publication] which have O abundances in Fulbright & Johnson (2003);
14 red giant stars in the Sgr dSph from the study by McWilliam &
Smecker-Hane (2005) [data kindly provided by the authors]. In
addition to these we have included additional data for stars in dSph
galaxies as well as for halo dwarf stars. Below we briefly comment on
the different studies and how many stars we have selected from each.
![]() |
Figure 11:
|
| Open with DEXTER | |
![]() |
Figure 12:
a) Shows the [Mn/O] vs [O/H] trend for dwarf stars in
the two disks and the halo. |
| Open with DEXTER | |
![]() |
Figure 13:
[Mn/O] vs. [O/H] for the Milky Way. The data are the same as
in Figs. 9a,b,c and Figs. 12a,b but here we distinguish the different
kinematic components instead of distinguishing the different
studies. |
| Open with DEXTER | |
In Fig. 11 we show the stellar parameters (
and
)
for all of the stars described above. As can be seen the
data span a large range in
.
In Fig. 12 we
have divided the stars according to their evolutionary status to see
if any, obvious, effects/trends in the abundance ratios are due to
stellar evolution. We used
to divide the stars into dwarfs
and giants. An arbitrary cut was imposed at 3.75 dex (compare
Fig. 11). From Fig. 12 we conclude that
although the dwarf star samples show much tighter and better defined
trends the more evolved stars overall follow the same trends.
Figure 12c shows that the RGB stars in
the dSph galaxies are enhanced in [Mn/O] relative to the Milky Way
giants in Fig. 12b. For example the dSph RGB stars
from all galaxies studied by Shetrone et al. (2003) have
and the RGB stars in the Sgr dSph have
(McWilliam et al. 2005) while the Milky Way halo giants from the compilation of data from Sobeck et al. (2006) (Mn
data) and Fulbright & Johnson (2003) (oxygen data) have
.
We note that the giant stars studied by Sobeck et al. (2006)
that
have
do not show enhanced [Mn/O] relative to the less
evolved stars in their sample. In fact, the stars with
are some of the stars with the lowest [Mn/O] in that sample. Hence,
it does appear that the dSph stars show an abundance pattern that is
enhanced in relation to what we see for similar stars in the Milky Way
and that the enhancement, probably, should be sought in environmental
effects (e.g. a slower star formation history) rather than being a
result of stellar evolution. However, to establish if the effect is
real further studies of halo and dSph RGB stars should be
undertaken. These studies should be done in a differential manner to,
to first order, be able to exclude modeling errors.
We note that one star with thin disk kinematics (Mn data from Sobeck et al. 2006) have an unusually high [Mn/O] (Fig. 9c). We have no explanation for this.
We use the final compilation of data for Milky Way stars in
Fig. 13 to investigate the origin of Mn in the Milky
Way. For the halo and metal-poor thick disk,
,
the
[Mn/O] trend is flat. This indicates that the production of Mn and O
are well balanced. Moreover, we know from the study of Bensby et al. (2004) that the archetypal signature of SN Ias in the thick disk
do not occur until [O/H] = 0. Hence the up-going trend we see after
must be interpreted as being due to metallicity
dependent Mn yields in SN II. The rising trend seen for the thin disk
sample could also be interpreted in this fashion. Although here we do
know that SN Ia contribute to the chemical enrichment and hence the
increase might also be due to these objects.
The dSph galaxies provide an interesting comparison. As discussed in Sect. 5.2, the RGB stars in the dSph appear to be genuinely more enhanced in [Mn/O] at a give [O/H] than the Milky Way halo and thick disk giants. The higher [Mn/O] values for the RGB stars in the dSph could be interpreted as being evidence for a slower star formation in those galaxies as compared to the star formation in our own halo and thick disk. For the Sgr dSph we note that the [Mn/O] vs. [O/H] trend does appear to resemble the trend for the thick disk (our data plus data from the literature and as shown in Fig. 12). Perhaps indicating that the Sgr dSph has a more Milky Way like star formation history than the other dSph. But, these are tentative conclusions and need further confirmation based on larger, and differential, abundance studies.
We have analysed four Mn I lines in 95 dwarf stars previously studied by Bensby et al. (2003, 2005). The stars were selected to have kinematics typical of the thick or the thin disk. Using these two well defined and well studied stellar samples we find that the abundance trends in the two samples differ such that the stars with thin disk kinematics are enhanced in Mn relative to stars with thick disk kinematics.
We also find that the previously reported "step'' in the [Mn/Fe] vs. [Fe/H] trend for disk stars in the Milky Way (see e.g. McWilliam et al. 2003) is in fact an artifact due to incomplete sampling of the two disk populations. Thus there is no need to invoke a large spread in the age-metallicity relation for the thin disk needed to explain such a "step'' (compare discussion in McWilliam et al. 2003).
Furthermore, when comparing the Mn abundances with Fe abundances the thick disk stars show a steadily increasing trend of [Mn/Fe] vs. [Fe/H] whilst the stars with kinematics typical for the thin disk show a flat trend up and until [Fe/H] = 0 and after that an increasing trend.
In order to further study the origin of Mn we have combined our new Mn abundances with O abundances. Fe is made both in SN II and in SN Ia. By using O, which is only made in SN II, as the reference element we simplify the interpretation of the abundance data. For our stars we took the O abundances from Bensby et al. (2004) and added data from a number of other studies of (mainly) giant stars in the disks and halo of the Milky Way as well as giant stars in dSph galaxies.
Our interpretation is that these data, to first order, can be explained by metallicity dependent yields in SN II. This is, essentially, in agreement with the conclusions in McWilliam et al. (2003).
However, it is not possible from these data to exclude that SN Ia contribute to the Mn production at least in the thin disk. Observations that would be helpful to further clarify the nature of Mn includes:
Acknowledgements
We would like to thank Bengt Gustafsson, Bengt Edvardsson, Kjell Eriksson, and Martin Asplund for usage of the MARCS model atmosphere program and their suite of stellar abundance (EQWIDTH) and synthetic spectrum generating programs. Jennifer Sobeck and co-authours are thanked for providing their Mn data prior to publication. Tammy Smecker-Hane and Andy McWilliam are thanked for sharing their Mn and O data for the Sgr dSph. SF is a Royal Swedish Academy of Sciences Research Fellow supported by a grant from the Knut and Alice Wallenberg Foundation. TB acknowledges support from the National Science Foundation, grant AST-0448900. This work has made use of the SIMBAD database operated at CDS, Strasbourg, France, and the VALD database at http://www.astro.uu.se/~vald.