Articles

CHEMICAL ABUNDANCES IN FIELD RED GIANTS FROM HIGH-RESOLUTION H-BAND SPECTRA USING THE APOGEE SPECTRAL LINELIST

, , , , , , , , , , and

Published 2013 February 13 © 2013. The American Astronomical Society. All rights reserved.
, , Citation Verne V. Smith et al 2013 ApJ 765 16 DOI 10.1088/0004-637X/765/1/16

0004-637X/765/1/16

ABSTRACT

High-resolution H-band spectra of five bright field K, M, and MS giants, obtained from the archives of the Kitt Peak National Observatory Fourier transform spectrometer, are analyzed to determine chemical abundances of 16 elements. The abundances were derived via spectrum synthesis using the detailed linelist prepared for the Sloan Digital Sky Survey III Apache Point Galactic Evolution Experiment (APOGEE), which is a high-resolution near-infrared spectroscopic survey to derive detailed chemical abundance distributions and precise radial velocities for 100,000 red giants sampling all Galactic stellar populations. The red giant sample studied here was chosen to probe which chemical elements can be derived reliably from the H-band APOGEE spectral region. These red giants consist of two K-giants (α Boo and μ Leo), two M-giants (β And and δ Oph), and one thermally pulsing asymptotic giant branch (TP-AGB) star of spectral type MS (HD 199799). Measured chemical abundances include the cosmochemically important isotopes 12C, 13C, 14N, and 16O, along with Mg, Al, Si, K, Ca, Ti, V, Cr, Mn, Fe, Co, Ni, and Cu. The K and M giants exhibit the abundance signature of the first dredge-up of CN-cycle material, while the TP-AGB star shows clear evidence of the addition of 12C synthesized during 4He-burning thermal pulses and subsequent third dredge-up. A comparison of the abundances derived here with published values for these stars reveals consistent results to ∼0.1 dex. The APOGEE spectral region and linelist is thus well suited for probing both Galactic chemical evolution, as well as internal nucleosynthesis and mixing in populations of red giants via high-resolution spectroscopy.

Export citation and abstract BibTeX RIS

1. INTRODUCTION

The Apache Point Observatory Galactic Evolution Experiment (APOGEE) is one of four experiments that are part of the Sloan Digital Sky Survey III (SDSS-III; Eisenstein et al. 2011). APOGEE is obtaining high-resolution (R ∼ 22, 300), high signal-to-noise ratio (S/N ⩾ 100 per pixel), H-band (λ1.51–1.69 μm) spectra of evolved, late-type stars, with the goal being to measure chemical abundances of 15 elements per star. Coupled with radial velocities that are accurate to ∼100 m s−1, APOGEE is creating the first high-precision spectroscopic and radial-velocity survey of all Galactic stellar populations (bulge, bar, disks, halo) using a uniform set of stellar tracers and spectral diagnostics, with a plan to observe and derive these parameters for 100,000 red giants by the end of 2014.

The physical and chemical parameters of the APOGEE stars are derived from a suite of software packages that, together, are called the APOGEE Stellar Parameters and Chemical Abundances Pipeline (ASPCAP); the details of this software and its analysis techniques will be discussed by A. Garcìa Pèrez et al. (in preparation). The first-generation output of ASPCAP, which will be based on one-dimensional (1-D) stellar atmosphere analyses in LTE, will consist of the stellar parameters of effective temperature, Teff, surface gravity (noted as log g), and the 1-D microturbulence parameter (ξ), along with the abundances of up to 16 elements (Fe, 12C, 14N, 16O, Mg, Al, Si, K, Ca, Ti, V, Cr, Mn, Co, Ni, or Cu) and, for many stars, the 13C isotopic abundance.

One of the central components of the ASPCAP machinery is a spectral linelist constructed to produce the synthetic spectra that are used to compare to the observed ones; the APOGEE linelist will be presented and discussed in detail by M. Shetrone et al. (in preparation).

The goal of this study is to analyze a small number nearby field red giants with well-defined stellar parameters and use the APOGEE linelist to derive detailed chemical abundances via standard "manual" abundance analysis techniques (Ramirez & Allende Prieto 2011 is a recent example using such techniques). The derived stellar parameters and chemical abundances can be compared to the same quantities derived via ASPCAP using the same spectra, linelist, and model atmosphere grid. This comparison will be an important test of ASPCAP. This work will also present new abundances in these stars based on an extensive and up-to-date spectral linelist in the near-IR H band.

Abundances will be presented in this paper using two different types of nomenclature, with one abundance scale defined as A(X) = log(N(X)/N(H)) + 12.0. The other scale compares the element ratio of N(X)/N(Q) in a target star to the Sun via [X/Q] = (A(X)TargetStarA(X))–(A(Q)TargetStarA(Q)).

2. OBSERVATIONAL DATA: THE FOURIER TRANSFORM SPECTRA OF BRIGHT FIELD RED GIANTS

The APOGEE high-resolution spectra cover the wavelength range from ∼15100 to 16900 Å and to analyze bright red giants over this spectral interval, data were obtained from the archives of the Fourier transform spectrometer (FTS; Hall et al. 1979) that operated at the coude focus of the Kitt Peak National Observatory Mayall 4 m reflector. A more complete description of observing techniques and the original data is discussed in Hinkle et al. (1995), in particular for the spectral atlas of α Boo.

The originally reduced spectra of selected stars were recorded as flux versus wavenumber (cm−1) over different dates and runs and covered a broad region in the H band. The spectra analyzed here were resampled in wavelength space and restricted to a wavelength range of 15000–17000 Å. The resolution of the spectra for different stars varies from R = 45,000 to 100,000, but were all higher than what is being observed for APOGEE. The higher resolution allows for a detailed investigation of the limits of what can be derived high-quality (i.e., high-resolution and high signal-to-noise) red giant spectra when using the APOGEE linelist.

3. THE "STANDARD STAR" PARAMETERS

3.1. The Choices

There are five stars analyzed here: the well-studied α Boo, a mildly metal-poor K1.5 giant, μ Leo, the prototype metal-rich star (Spinrad & Taylor 1969), which is a K1 giant, two near-solar-metallicity M giants, β And (M0III) and δ Oph (M0III), and HD199799, an asymptotic giant branch (AGB) star of spectral type MS with a mild enhancement of 12C and a slightly elevated ratio of C/O (relative to solar) due to third dredge-up (Smith & Lambert 1990). Figure 1 illustrates a small portion of the APOGEE wavelength coverage for four of the red giants here. This region from 15560 to 15590 Å includes representative species and how they vary with spectral type (or effective temperature) and metallicity. Lines from the main CNO-containing molecules, CO, OH, and CN fall within this region, as noted in the figure. Of particular interest is the metal-rich K-giant μ Leo, which exhibits quite strong CN lines; much of the variation in the spectrum of μ Leo compared to the other giants result from the high line density of CN. This figure also points to the desirability of analyzing large numbers of red giants via spectrum synthesis techniques using a detailed linelist.

Figure 1.

Figure 1. Sample spectra of red giants in the region near the 12C16O (3–0) bandhead. These stars span a range in temperature, from Teff = 4550 K (μ Leo) down to Teff = 3825 K (β And) and in metallicity, with α Boo having [Fe/H] = −0.52, while μ Leo has [Fe/H] = +0.31 dex. These panels illustrate both the quality of the spectra and the variation of molecular and atomic line absorption with primarily Teff.

Standard image High-resolution image

With a significant fraction of its red giant targets in the thin and thick disks, plus the inner bulge and bar, the five types of red giants studied here span a range of temperature and metallicity that provides a good test for ASPCAP. The APOGEE survey targets red giants with effective temperatures from Teff = 3400 K–5000 K, having surface gravities from log g = 3.0 to −0.5, and metallicities from [Fe/H] = −5.0 to +1.0. These nearby giants have well-measured photometry and parallaxes, so temperatures, luminosities, and masses are all rather well constrained. In addition, four of the five stars have previous abundance analyses in the literature and thus provide comparison tests for the abundances derived from the APOGEE linelist.

3.2. Basic Data

Deriving stellar chemical abundances from high-resolution spectra requires as basic input the effective temperature (Teff), surface gravity (characterized as log g with g as cm s−2), and overall metallicity as the fundamental parameters of the model atmosphere. In addition, high-resolution abundance analyses of cool stars that use static 1-D stellar atmospheres also require the derivation of a non-thermal Doppler-like broadening term called microturbulence (ξ), which is determined as part of the spectral analysis. The first step is the determination of Teff and, for these relatively nearby red giants, is set by the near-infrared magnitudes of J and K. This particular choice is used because the APOGEE targets are selected from the Two Micron All Sky Survey (2MASS) catalog (Skrutskie et al. 2006) and therefore all have (JK) available. The aim is to focus on utilizing near-IR data and then compare to published literature results.

All of the red giants here have parallaxes measured by Hipparcos and thus have distances known to various levels of accuracy. The parallax is used to set the distance and thus MK. Bolometric corrections from Bessell et al. (1998) are used in conjunction with (JK) to determine Mbol and consequently stellar luminosity. Basic data for the red giants are presented in Table 1.

Table 1. Red Giant Standard Stars: Observed Properties

Star HR SpT π d JKSb Mbolc
(mas)a (pc)
β And 337 M0III 16.4 ± 0.8 61 ± 3 0.98 −3.06 ± 0.11
μ Leo 3905 K2III 26.3 ± 0.2 38 ± 1 0.68 0.41 ± 0.05
α Boo 5340 K2III 88.8 ± 0.5 11.3 ± 0.1 0.79 −0.92 ± 0.03
δ Oph 6056 M0III 19.1 ± 0.2 52.5 ± 0.5 0.97 −2.16 ± 0.05
HD 199799 ... M2S 2.6 ± 0.6 386 ± 92 1.25 −3.64 ± 0.54

Notes. aParallax from van Leeuwen (2007). bJohnson (1966) transformed to the 2MASS system using Carpenter (2001). cCalculated from MK with bolometric corrections from Bessell et al. (1998).

Download table as:  ASCIITypeset image

The most important derived quantity affecting the overall shape of the spectrum is Teff, which here is based on an average of two Teff–(JK) calibrations; one from Gonzalez Hernandez & Bonifacio (2009) and the other from Bessell et al. (1998). The Gonzalez Hernandez & Bonifacio calibration is defined in the 2MASS photometric system of (JKS), so the original Johnson (1966) (JK) colors and K magnitudes for these bright red giants (which have saturated magnitudes in the 2MASS catalog) were transformed to 2MASS values using the prescription from Carpenter (2001). The differences between the calibrations of Bessell et al. (1998) and Gonzalez Hernandez & Bonifacio (2009), for a given metallicity, are less than 30–50 K for Teff ⩽4500 K, but above this temperature the Bessell et al. (1998) becomes cooler by about 100 K for a given (JKS) (where we have transformed the Bessell et al. scale onto the 2MASS system). There is also a small metallicity effect on the (JKS) calibration which, for a given (JKS), can be ∼100 K over the broad metallicity range of [m/H] = −2.0 to +0.0 (where m is taken to represent the generic abundance of metals, i.e., those elements other than H and He). Given the differences between the two independent Teff calibrations of less than 30–50 K over the temperature range of the target stars, as well as their modest range in metallicity, and well-determined magnitudes with ∼0.02 uncertainty in (JKS) (which translates to ΔTeff = 30 K), the derived values of Teff have uncertainties of ±50 K or less.

Given the stellar luminosity and effective temperature, these fundamental properties are then compared to stellar evolution models presented in Bertelli et al. (2008, 2009; with data from these papers taken from the Web site http://stev.oapd.inaf.it/YZVAR) to estimate stellar masses and, ultimately, surface gravities through

Determining the surface gravity from evolutionary tracks is an iterative process because it requires knowledge of the stellar metallicity. This is illustrated in Figure 2 where the stellar luminosity and Teff are plotted along with evolutionary models from a range of masses, with each panel having a different metallicity: Z = 0.017 (top panel), 0.008 (middle panel), and 0.004 (bottom panel), corresponding to overall metallicities of [m/H] = +0.07, −0.26, and −0.56, respectively, when taking Z = 0.0145 for the Sun, as suggested by Lodders (2010). The illustrated evolutionary tracks follow stellar evolution up to the tip of the first-ascent red giant branch (which we label RGB). Initial model metallicities for the stars studied here are based on previously published results, with β And, δ Oph, and HD199799 started with solar metallicity, α Boo with [Fe/H] = −0.5, and μ Leo with [Fe/H] = +0.3. As there are numerous Fe i lines in this spectral region, a sample of Fe i lines are used to determine the iron abundance. Given the limited range of metallicity spanned by the nearby red giants studied here, the Fe abundance is a good proxy for the overall stellar metallicity. If the derived Fe abundance differed by more than 0.1 dex from the assumed value, a new surface gravity was derived and a new model generated in order to re-analyze the Fe i lines until convergence between model abundance and derived abundance. For these relatively well-studied red giants, no more than two iterations were required until convergence. Table 2 presents, for these red giants, the derived parameters Teff, log (L/L), mass, log g, microturbulent velocity (ξ; see Section 4.2), and overall metallicity, which is represented by [Fe/H] as discussed in Section 4.2.

Figure 2.

Figure 2. Positions of the field red giants in the log (L/L)–Teff plane, which is a modified H-R diagram. The continuous curves represent stellar evolution tracks for different mass stars with each panel representing a different heavy-element mass fraction (Z); this is equivalent to overall stellar metallicity with the Sun having Z = 0.0145 (Lodders 2010). The model tracks are used to set log g for each red giant studied here.

Standard image High-resolution image

Table 2. Red Giant Standard Stars: Derived Parameters

Star log M/M Teff log g ξ [Fe/H]a
(L/L) (K) (cm s−2) (km s−1)
β And 3.11 ± 0.05 1.9 ± 0.2 3825 ± 75 0.9 ± 0.1 2.20 ± 0.10 −0.2
μ Leo 2.68 ± 0.03 2.0 ± 0.3 4550 ± 50 2.1 ± 0.1 1.80 ± 0.15 +0.3
α Boo 2.24 ± 0.02 1.1 ± 0.2 4275 ± 50 1.7 ± 0.1 1.85 ± 0.05 −0.4
δ Oph 2.75 ± 0.03 1.5 ± 0.2 3850 ± 75 1.2 ± 0.1 1.90 ± 0.10 +0.0
HD 199799 3.34 ± 0.22 1.5 ± 0.2 3400 ± 75 0.5 ± 0.2 2.40 ± 0.20 −0.2

Note. a[Fe/H] represents the overall metallicity used in the final model atmosphere. Final Fe abundances are given in Table 6.

Download table as:  ASCIITypeset image

3.3. Stellar Parameters and Evolutionary State

Given the position of the stars in the log L versus Teff plane shown in Figure 2, it can be ascertained that four stars (α Boo, μ Leo, δ Oph, and β And) are most likely first-ascent red giants of low to moderate mass. Based on its warmer temperature and somewhat lower luminosity, it may be more likely that μ Leo is in or very near the core-He burning phase (clump giant) and has thus already ascended the RGB and experienced the core-He flash. It is also possible that the other three stars already noted above may be post-He core-burning stars ascending the AGB; however, lifetimes on the two separate evolutionary sequences would favor them as being on the RGB. All fall above the "luminosity bump" (Fusi Pecci et al. 1990) and would have thus experienced any extra mixing that may occur during that phase of stellar evolution.

The coolest and most luminous red giant in this sample, HD199799, is an AGB star that is experiencing the third dredge-up, as evidenced by the observation of Tc I in its spectrum and its elevated s-process abundances, such as [Y/Fe] ∼ +0.5 or [Nd/Fe] ∼ +0.4 (Smith & Lambert 1988, 1990).

4. THE ABUNDANCE ANALYSIS

As with the initial abundances to be derived from APOGEE, the analysis here will rely on 1-D model atmospheres computed in LTE. Within APOGEE, the choice was made to base the initial abundance determinations on an extensive model atmosphere database generated using the ATLAS9 code (Kurucz 1993), thus the model atmospheres used here are based on these APOGEE models. ATLAS9 is widely used as a universal LTE 1D plane-parallel atmosphere modeling code. The details of constructing the atmosphere database for APOGEE using ATLAS9 are presented in Meszaros et al. (2012).

The individual elemental abundances themselves will also be derived under the assumption of LTE with the spectrum synthesis code MOOG as the analysis tool (Sneden 1973 is the original reference, but updated information, data, and the most recent versions of the code can be found at http://www.as.utexas.edu/~chris/moog.html).

An extensive effort was made to produce as complete a spectral linelist as possible, so the analysis technique used here, as in APOGEE in general, is to utilize spectrum synthesis in a quantitative comparison between synthetic and observed spectra. All abundances are determined for each spectral line, or group of closely spaced molecular lines, via minimizing the residuals between observed and synthetic spectra as a function of the abundance in question.

4.1. Choice of Fe i Lines and the Microturbulence

Due to both its relatively large cosmic abundance and large number of energy levels, iron has historically been used as a diagnostic species for certain stellar parameters, as well as an overall metallicity indicator, so the results for Fe i lines selected from the APOGEE linelist (M. Shetrone et al., in preparation) are discussed here. In particular, Fe i lines are well suited for setting the microturbulent velocity, ξ, which is required for 1-D model atmosphere abundance analyses.

A large number of Fe i lines fall within the APOGEE spectral window (15100–16900 Å), although there are no Fe ii lines strong enough to be detected in these red giants. Since all abundance determinations are done via spectrum synthesis, the Fe i lines were culled to include those lines that contain only contributions from Fe i (this was determined by preliminary synthesis of the red giant spectra, then setting the Fe abundance to zero and ensuring that the spectral feature in question vanished); since the Fe i lines are so numerous, this stringent selection still results in an adequate list of lines. Table 3 presents the Fe i lines used for the iron abundance determinations. As is true of most of the atomic lines in the H band, the excitation energies are dominated by rather high-excitation lines, although for Fe i there are two lower-excitation lines.

Table 3. Fe i Lines used in the Abundance Determinations

λ χ log gf α Boo β And δ Oph μ Leo HD 199799
(Å) (eV)
15194.492 2.223 −4.779 7.04 7.22 7.37 7.83 7.05
15207.526 5.385 +0.080 7.04 7.29 7.40 bla 7.20
15395.718 5.620 −0.341 6.96 7.27 bla 7.88 7.42
15490.339 2.198 −4.807 6.99 7.23 7.55 7.76 7.50
15648.510 5.426 −0.701 6.97 7.20 7.35 7.75 7.03
15964.867 5.921 −0.128 6.93 7.22 7.35 7.65 7.33
16040.657 5.874 +0.066 6.94 7.20 bla 7.69 7.14
16153.247 5.351 −0.743 6.94 7.22 7.53 7.75 7.32
16165.032 6.319 +0.723 6.98 7.24 7.56 7.73 7.30

Note. abl: feature too blended to use in cooler or very metal-rich red giants.

Download table as:  ASCIITypeset image

The Fe i lines span a large enough range in line strength that they can be used to set the microturbulent velocity. This parameter is defined by the value of ξ that produces no trend in the Fe abundance as a function of line strength. For a range of microturbulent velocities Fe abundances were determined from each Fe i line and the adopted value of ξ was one in which a linear regression between A(Fe) as a function of log(Wλ/λ) yielded a slope of zero. This procedure is illustrated graphically in Figure 3, where the abundances derived for each individual line are plotted versus ξ (the continuous lines); the circle with error bars shows where the slope of A(Fe) with line strength goes to zero.

Figure 3.

Figure 3. Individual Fe i abundances as a function of microturbulent velocity for μ Leo (top panel) and β And (bottom panel). Each curve represents a single Fe i line. Stronger lines are more sensitive to changes in ξ than weaker lines. The value of ξ that yields no significant abundance differences as a function of reduced line strength (log Wλ/λ) is taken as the characteristic microturbulent velocity for the star. The derived microturbulent velocity as displayed here also corresponds to the smallest scatter in A(Fe).

Standard image High-resolution image

The microturbulent velocities set by the Fe i lines, as well as the Fe abundances are listed in Table 2 as part of the derived set of stellar parameters. The individual line-by-line abundances as determined for each Fe i line are also listed in Table 3. The Fe abundances were used to set the final overall metallicities of the ATLAS9 model atmospheres for the analysis of additional elements via atomic or molecular lines.

4.2. The C, N, and O Abundances

Abundances for 12C, 14N, 16O, and the minor carbon isotope 13C are derived from combinations of vibration–rotation (V-R) lines of CO (X1Σ+) and OH (X2Π), along with electronic transitions of CN (A2Π–X2Σ). Although the details of the ingredients that went into these molecular lines in the APOGEE linelist will be found in M. Shetrone et al. (in preparation), a few highlights are noted here. The adopted dissociation energies (D0) are 11.092 eV for CO, 4.395 eV for OH, and 7.76 eV for CN. The gf-values are from Goorvitch (1994) for CO, Goldman et al. (1998) for OH, and the CN gf-values are taken from the baseline linelist of Kurucz (1993), with updates to these values from the prescription given in Melendez & Barbuy (1999).

The procedure used for the CNO analysis, since all red giants studied here have C/O ⩽ 1, is to use CO to set the carbon abundance. With that carbon abundance, OH provides an O abundance, which, if different from the initial value (taken to be the solar value scaled by the stellar value of [Fe/H]), is used to re-analyze the CO lines. This process is repeated until the C and O abundances yield consistent abundances from CO and OH. With these values of C and O, the CN lines are used to derive the abundance of nitrogen. In general, the abundance of N has little to no effect on the CO and OH lines, however, the respective final C, N, and O abundances all provide self-consistent results from CO, OH, and CN.

All spectral features were synthesized and the molecular regions that were used to extract the abundances are listed in Table 4. In this case, if a range of wavelengths is given it indicates that a spectral window containing numerous molecular lines was used (for CO and OH), while for CN individual transitions were identified and synthesized. As with Table 3, the abundances derived for each feature are also listed in Table 4. In the case of 13C, as has been customary for minor isotopes, the abundance is listed as the ratio of 12C/13C, with A(12C) being set by the mean abundance of carbon-12.

Table 4. Molecular Lines and Features used for 12C, 13C, 14N, and 16O

Molecular Lines λ–interval α Boo β And δ Oph μ Leo HD 199799
(Å)
12C from 12C16O lines
(3–0) V-R 15578–15586 7.86 8.05 8.27 wa 8.51
(4–1) V-R 15774–15787 7.97 8.11 8.24 8.46 8.51
(5–2) V-R 15976–16000 7.90 8.07 8.23 8.56 8.38
(6–3) V-R 16183–16196 7.90 8.03 8.22 8.54 8.43
12C/13C ratios from 13C16O and 13C14N lines
13C16O(3–0) V-R 15922–15926 7 17 14 wa 30
13C16O(4–1) V-R 16120–16125 6 15 11 wa 32
13C16O(6–3) V-R 16740–16747 6 13 11 wa 20
13C14N 15314–15315 wa wa wa 22 wa
13C14N 15354–15356 wa wa wa 17 wa
16O from 16OH lines
(2–0) P1 9.5 15277–15282 8.55 8.83 8.78 9.10 8.80
(3–1) P2 3.5 15390–15392 8.54 8.78 8.76 8.93 8.60
(3–1) P2 4.5 15504–15507 8.48 8.71 8.73 9.15 8.58
(2–0) P2 11.5 15568–15573 8.51 8.80 8.78 8.98 8.52
(3–1) P2 9.5 16189–16193 8.47 8.78 8.81 9.07 8.63
14N from 12C14N lines
(1–2) Q2 41.5 15260. 7.61 8.10 8.23 8.83 8.31
(1–2) P2 34.5 15322. 7.68 8.07 8.21 8.88 8.35
(1–2) R2 56.5 15397. 7.65 8.06 8.26 8.63 8.28
(0–1) R1 68.5 15332. 7.63 7.95 blb 8.68 8.35
(0–1) P2 49.5 15410. 7.66 8.03 8.13 8.80 blb
(0–1) Q2 59.5 15447. 7.64 7.98 8.28 8.73 blb
(0–1) Q1 60.5 15466. 7.65 8.06 8.23 8.51 blb
(1–2) P2 38.5 15472. 7.71 8.13 8.11 8.68 8.25
(0–1) P1 51.5 15482. 7.56 8.03 8.15 8.68 8.25

Notes. aw: feature too weak to use. bbl: feature too blended by nearby lines.

Download table as:  ASCIITypeset image

4.3. The Other Atomic Lines and Elemental Abundances

With the Fe and C, N, and O abundances in hand, additional elements are derived via atomic spectral lines in the wavelength interval 15100–16900 Å. The APOGEE linelist was searched for suitable species, all neutral over the Teff range of the red giants here, and some 12 other elements were deemed detectable and able to be analyzed quantitatively: Mg i, Al i, Si i, K i, Ca i, Ti i, V i, Cr i, Mn i, Co i, Ni i, and Cu i. The additional atomic lines are listed in Table 5. As with the molecular lines, if a precise wavelength is shown, the transition in question is represented by a single line, however, if the wavelength is approximated by an integer value it indicates multiple transitions due to hyperfine (hfs) splitting or isotopic splitting. In the case of isotopic components, solar-system isotopic fractions are assumed and, in general, this assumption has no significant effect on the derived total atomic abundances. Figure 4 illustrates the fitting procedure for an atomic line from Mn I in δ Oph, with the underlying hfs components illustrated, as well as nearby blending features, with several synthetic spectra plotted which have differing Mn abundances. As with the previous Tables 3 and 4, individual abundances are provided for each line or feature. As noted previously, the details of the generation of the APOGEE linelist will be found in M. Shetrone et al. (in preparation).

Figure 4.

Figure 4. Comparison of the observed spectrum near the Mn i line at 15262 Å in δ Oph (blue open circles) with synthetic spectra which were computed for Mn abundances separated in intervals of 0.1 dex. This Mn i transition consists of a number of hfs components of differing wavelengths, as shown, with the vertical length of each component proportional to the gf-value of that component transition. Note that Mn I is blended with a CN feature containing both 12C14N and 13C14N lines, thus the value of 12C/13C ratio must be included. Both the hfs structure and CN blending illustrate the need for the APOGEE analysis to be based on spectrum synthesis with an accurate linelist.

Standard image High-resolution image

Table 5. Atomic Lines used and Derived Abundance

Element λ χ log gf α Boo β And δ Oph μ Leo HD 199799
(Å) (eV)
Mg i 15740.716 5.931 −0.262 7.15 7.25 7.53 7.80 7.29
  15748.9 5.932 +0.276 7.09 sa 7.57 sa 7.23
  15765.8 5.933 +0.504 7.04 sa sa 7.87 blb
  15879.5 5.946 −1.248 7.13 7.20 7.45 7.82 blb
  15886.2 5.946 −1.555 7.22 7.21 7.53 blb 7.18
  15889.485 5.946 −2.013 blb 7.30 7.55 7.90 wc
  15954.477 6.588 −0.807 7.25 7.36 7.61 blb 7.45
Al i 16718.957 4.085 +0.290 6.15 6.15 sa 6.85 blb
  16763.359 4.087 −0.524 6.16 6.09 6.45 6.94 6.10
Si i 15361.161 5.954 −1.925 6.99 blb blb blb wc
  15376.831 6.222 −0.649 7.18 7.17 7.57 7.66 wc
  15833.602 6.222 −0.168 6.93 6.95 7.28 7.72 7.20
  15960.063 5.984 +0.107 7.07 7.14 7.43 7.88 wc
  16060.009 5.954 −0.566 7.21 7.25 7.59 7.74 7.33
  16094.787 5.964 −0.168 7.11 7.06 7.46 7.66 7.32
  16215.670 5.964 −0.665 7.15 7.23 7.62 blb blb
  16680.770 5.984 −0.140 7.11 7.25 7.62 7.83 7.35
  16828.159 5.984 −1.102 7.13 7.17 7.45 7.85 blb
K i 15163.067 2.670 +0.524 4.81 4.82 5.18 5.63 5.20
  15168.376 2.670 +0.347 4.77 4.89 blb blb blb
Ca i 16136.823 4.531 −0.552 5.75 6.09 6.23 6.69 6.06
  16150.763 4.532 −0.229 5.88 5.95 6.21 6.53 6.07
  16155.236 4.532 −0.619 5.82 6.02 6.30 6.61 6.16
  16157.364 4.554 −0.208 5.90 6.02 6.23 6.63 6.16
Ti i 15543.756 1.879 −1.120 4.53 4.74 5.05 5.58 4.74
  15602.842 2.267 −1.643 4.66 4.74 5.03 5.43 blb
  15698.979 1.887 −2.060 4.52 4.69 5.03 5.13 blb
  15715.573 1.873 −1.250 4.57 4.71 5.14 blb 4.70
  16635.161 2.345 −1.807 4.65 blb 5.10 5.46 4.90
V i 15924. 2.138 −1.175 3.61 3.66 3.86 4.18 blb
Cr i 15680.063 4.697 +0.179 5.08 5.38 5.62 6.11 5.20
  15860.214 4.697 −0.012 5.14 5.39 5.57 6.16 5.54
Mn i 15159. 4.889 +0.544 4.88 5.19 blb 5.82 blb
  15217. 4.889 +0.414 4.91 5.26 5.47 5.83 5.32
  15262. 4.889 +0.345 4.78 5.08 5.21 5.71 5.30
Co i 16757.7 3.409 −1.230 4.50 4.73 4.85 5.23 4.50
Ni i 15605.680 5.299 −0.376 5.76 6.08 6.18 6.69 6.05
  15632.654 5.305 −0.106 5.82 5.99 6.20 6.65 6.02
  16584.439 5.299 −0.528 5.71 6.00 6.21 6.58 6.12
  16589.295 5.469 −0.600 5.81 6.03 6.23 6.72 5.98
  16673.711 6.029 +0.317 5.72 blb blb 6.39 6.12
  16815.471 5.305 −0.606 5.79 5.94 6.06 6.58 5.95
  16818.760 6.039 +0.311 5.77 blb 6.22 blb 5.95
Cu i 16005.7 5.348 −0.131 3.55 4.01 4.18 4.41 4.20

Notes. as: feature too strong to use. bbl: feature too blended by nearby lines. cw: feature too weak.

Download table as:  ASCIITypeset image

Table 6 provides a summary table of all abundances listed as their mean values and standard deviations from the lines or spectral intervals shown in Tables 3– 5. In cases where only one line or region was observed, no standard deviation is given. The abundance of 13C is given as the isotopic ratio of 12C/13C.

Table 6. Chemical Abundances

Element α Boo β And δ Oph μ Leo HD 199799
Fe 6.98 ± 0.04 7.23 ± 0.03 7.44 ± 0.09 7.76 ± 0.08 7.25 ± 0.12
12C 7.96 ± 0.05 8.06 ± 0.03 8.24 ± 0.02 8.52 ± 0.05 8.46 ± 0.06
12C/13C 6.3 ± 0.6 15 ± 2 12 ± 2 20 ± 3 27 ± 6
14N 7.64 ± 0.04 8.05 ± 0.06 8.20 ± 0.06 8.71 ± 0.11 8.31 ± 0.05
16O 8.64 ± 0.04 8.78 ± 0.04 8.77 ± 0.03 9.05 ± 0.09 8.63 ± 0.09
Mg 7.15 ± 0.08 7.26 ± 0.07 7.54 ± 0.05 7.85 ± 0.05 7.28 ± 0.11
Al 6.16 ± 0.01 6.12 ± 0.03 6.45 6.90 ± 0.05 6.10
Si 7.12 ± 0.07 7.18 ± 0.07 7.53 ± 0.08 7.76 ± 0.09 7.30 ± 0.07
K 4.79 ± 0.02 4.86 ± 0.04 5.18 5.63 5.20
Ca 5.84 ± 0.07 6.02 ± 0.06 6.24 ± 0.04 6.62 ± 0.07 6.11 ± 0.05
Ti 4.59 ± 0.07 4.72 ± 0.02 5.07 ± 0.05 5.40 ± 0.19 4.78 ± 0.10
V 3.61 3.66 3.86 4.18 ...
Cr 5.11 ± 0.03 5.38 ± 0.01 5.60 ± 0.03 6.14 ± 0.03 5.37 ± 0.17
Mn 4.86 ± 0.07 5.18 ± 0.09 5.34 ± 0.13 5.79 ± 0.07 5.31 ± 0.01
Co 4.50 4.73 4.85 5.23 4.50
Ni 5.77 ± 0.04 6.01 ± 0.05 6.18 ± 0.06 6.60 ± 0.12 6.03 ± 0.07
Cu 3.55 4.01 4.18 4.41 4.20

Download table as:  ASCIITypeset image

4.4. Abundance Sensitivity to Stellar Parameters

Most of the abundances presented here are derived from more than one atomic or molecular line or features of lines, thus internal consistency in the average abundances is defined by the scatter from the individual features. This scatter is characterized by the standard deviations of Table 6. Because multiple features arise from differing line strengths (weak versus strong), excitation potentials, or ionization fractions, these respective standard deviations provide some estimate as to how well the 1-D LTE analysis can recover the details of the underlying stellar spectrum from the basic model atmosphere parameters Teff, log g, metallicity, and derived value of ξ. The typical standard deviations from Table 6 are ∼ ±0.05 dex. Of note is the metal-rich K-giant μ Leo, which exhibits exceptionally strong CN lines, that are numerous throughout this H-band window. The background CN absorption creates blending which lowers the overall signature of other abundances and increases the uncertainty in these other abundances; this is reflected in the generally larger standard deviations found in the μ Leo abundances.

Besides internal consistency from the various spectral lines of different species, it is useful to explore their sensitivities to stellar parameters in order to characterize the magnitudes of possible systematic effects due to uncertainties in the four primary parameters of effective temperature, surface gravity, metallicity, and 1-D derived microturbulent velocity. For this discussion Teff is written as T, log g as G, model atmosphere metallicity ([m/H]) as m, and microturbulent velocity as ξ. Because the abundance for a given element, A = log [N(A)/N(H)] + 12. = A(T, G, m, ξ), the incremental change in this abundance due to small variations in the parameters is given by

The scatter in the abundance caused by small changes in the stellar parameters is then defined as

and dA2 can be written as

The partial derivative terms are calculated by changing each model parameter and determining the resulting change in logarithmic abundance, A. Since the perturbation in stellar parameters is relatively small, the changes are approximated by linear trends. As this paper is a limited study of a few stars as a test of the APOGEE linelist, a test of abundance sensitivity to stellar parameters is carried out for a representative model taken with Teff = 4000 K, log g = 1.3, [m/H] = +0.0, with ξ = 2.0 km s−1. The predicted equivalent widths were generated for the lines used in the abundance analysis and the stellar parameters were then perturbed by dT = +50 K, dG = +0.2 dex, dm = +0.1 dex, and dξ = +0.2 km s−1, with new abundances derived for each separate perturbation, giving the change in abundance for each parameter change. These coefficients, which reflect the change in abundance, are listed in Table 7 and provide an overall view of how sensitive the derived abundances are to changes in stellar parameters that are representative of realistic errors for these particular red giants.

Table 7. Abundance Sensitivity: Teff = 4000 K, log g = 1.3, ξ = 2.0 km s−1, [m/H] = 0.0

Species A/∂T A/∂G A/∂ξ A/∂m dA' dA
(+50 K) (+0.2 dex) (+0.2 km s−1) (+0.1 dex)
Fe i +0.015 +0.003 −0.056 +0.032 0.066 0.087
CO +0.032 +0.072 −0.011 +0.041 0.089 0.167
CN −0.018 +0.064 −0.006 +0.038 0.077 0.127
OH +0.124 −0.022 −0.079 +0.111 0.186 0.363
Mg i +0.017 −0.031 −0.042 +0.029 0.062 0.098
Al i +0.055 −0.077 −0.052 +0.027 0.111 0.207
Si i +0.003 −0.013 −0.039 +0.038 0.056 0.065
K i +0.042 −0.046 −0.048 +0.014 0.081 0.143
Ca i +0.045 −0.046 −0.032 +0.017 0.074 0.143
Ti i +0.082 −0.004 −0.111 +0.015 0.139 0.168
V i +0.076 +0.009 −0.054 +0.011 0.094 0.132
Cr i +0.038 −0.018 −0.031 +0.016 0.055 0.100
Mn i +0.032 −0.057 −0.056 +0.041 0.095 0.163
Co i +0.026 +0.033 −0.075 +0.034 0.092 0.136
Ni i −0.005 +0.017 −0.043 +0.031 0.056 0.068
Cu i +0.001 +0.032 −0.013 +0.019 0.039 0.044

Download table as:  ASCIITypeset image

Columns 6 and 7 in Table 7 provide the values of dA two cases: one using only the first four terms which contains the change in abundance due to changes in each primary parameter (dA' in Column 6, which is valid if all parameters are independent of each other), while Column 7 includes the covariant terms. Within the analysis technique applied here, it was found that the microturbulent velocity did not measurably depend on the derived effective temperature or gravity, and the dTdξ and dGdξ cross terms are not included. As can be seen by the comparison of dA' and dA, the covariant terms contribute significantly, in many species, to the uncertainty. The largest covariant terms are typically those between Teff and log g, caused by the slope of the RGB in the Teff–log(L/L) plane (since LT4eff/g).

4.5. Comparison to Previously Published Abundances

The abundances derived here rely on the APOGEE H-band spectral window from 15100 to 16900 Å, along with parallaxes and photometry. Historically, most detailed abundance patterns obtained from high-resolution spectroscopy have tended to be based on wavelengths below 10000 Å, so it is useful to compare the abundances based on the APOGEE linelist with those from other studies for stars in common. Within the sample of red giants discussed here, previous high-resolution abundance analyses have been conducted on α Boo, μ Leo, β And, and HD 199799.

4.5.1. α Boo

Being the nearest red giant, α Boo is often used as the standard star of choice for abundance comparisons. Ramirez & Allende Prieto (2011) have produced a recent and detailed high-resolution optical spectral analysis, reviewing both fundamental stellar parameters for this star, as well as deriving an extensive set of elemental abundances. They used the observed spectral energy distribution and compared to theoretical ones from the Kurucz grid of no-overshoot model atmospheres with α-element enhanced compositions ([α/Fe] = +0.4; Castelli & Kurucz 2003) in order to set the Teff for α Boo. The flux from λ ∼0.3–10 μm was used with the resultant best-fit Teff = 4286 ± 30 K. This value is entirely consistent with that obtained here (Teff = 4275 ± 50 K) based on the (JK) color and calibration used here from Bessell et al. (1998) and Gonzalez Hernandez & Bonifacio (2009).

Ramirez & Allende Prieto (2011) used stellar evolutionary model tracks to estimate the mass (as well as age) of α Boo, as the effective temperature and parallax are known to high precision. The mass was set by comparison to Yonsei-Yale isochrones (Yi et al. 2001; Kim et al. 2002) and was found to be M = 1.08 ± 0.06 M, which yields a surface gravity of log g = 1.66 ± 0.05. In this study the Teff = 4275 K value from the (JK) color is adopted, while the Ramirez & Allende Prieto log g is used. These stellar parameters were used to generate a model atmosphere as described in Section 4 with the ATLAS9 code in the same manner as is being done for the APOGEE targets (Meszaros et al. 2012).

The microturbulence, ξ, derived by Ramirez & Allende Prieto was set by Fe i lines and was found to be ξ = 1.74 km s−1. The Fe i lines used from the APOGEE linelist yielded a value of ξ = 1.85 ± 0.05 km s−1: very close to that derived from optical Fe i lines.

A direct comparison of the derived abundances here with those from Ramirez & Allende Prieto (2011) using completely different sets of lines leads to the following mean difference and standard deviation for 13 elements in common of ΔA(x)(this paper; Ramirez & Allende Prieto) = +0.03 ± 0.10 dex. There is no significant offset in the overall abundance patterns and the scatter is that expected given the uncertainties in the log gf scales and internal line-to-line scatter within each study. This comparison suggests that abundances derived from the APOGEE linelist can be used for comparisons to those found from analyses at visual wavelengths when high-quality spectra are analyzed.

4.5.2. μ Leo

This K-giant is often cited as the "prototype" for metal-rich stars, having been noted as such by Spinrad & Taylor (1969), and some of its spectral regions contain numerous and strong CN lines. Two older abundance studies in the literature (Gratton & Sneden 1990; Smith & Ruck 2000) are used to compare to the abundances derived here, as well as a more recent analysis by Fulbright et al. (2007; FMR). As with the α Boo comparison in Section 4.5.1, both of these published analyses are based on high-resolution spectra in the visual. Both studies used Teff = 4540 K for μ Leo, which is similar to the value 4550 K adopted here. The surface gravities were slightly different, with log g = 2.3 ± 0.3 in Gratton & Sneden and 2.2 ± 0.1 in Smith & Ruck; both values agree with the value here, within the uncertainties, of log g = 2.1 derived from the luminosity and estimate of stellar mass.

The value of ξ = 1.8 km s−1 found from the H-band Fe i lines is slightly higher than 1.2 km s1 found by both Gratton & Sneden (1990) and Smith & Ruck (2000), but with all studies having uncertainties of 0.3–0.5 km s−1, this offset cannot be claimed to be significant. A comparison of abundances between here and Gratton & Sneden (1990) for 12 elements finds a mean ΔA(x)(this paper; GS) = −0.05 ± 0.14 dex, while the same comparison with Smith & Ruck (2000) yields ΔA(x) = −0.05 ± 0.02 dex for the three elements analyzed by them: Fe, Mg, and Ca. The abundance offsets are not significant and the scatter is within what is expected.

The more recent abundance analysis of μ Leo by FMR was part of their study of bulge red giants. Their derived stellar parameters were quite similar to those found here, with Teff(FMR) = 4520 K, (ΔT(this study; FMR) = +30 K), log g = 2.33 (Δlog g(this study; FMR) = −0.23), and ξ = 1.50 km s−1 (Δξ(this study; FMR) = +0.3 km s−1). There are seven elements in common (Fe, O, Mg, Al, Si, Ca, and Ti) and the mean difference is ΔA(x)(this study; FMR) = −0.10 ± 0.13 dex. A comparison among the three μ Leo studies points to rather small offsets in stellar parameters (∼30 K in Teff, 0.2 dex in log g, 0.3–0.5 km s−1 in ξ, and 0.1 dex in [m/H]) with the conclusion that H-band spectroscopy can be compared reliably to optical analyses with reasonable accuracies, even at this higher metallicity.

4.5.3. β And

This M-giant was analyzed previously by Smith & Lambert (1985) who derived Teff = 3800 K (from its (VK) color), log g = 1.6, and ξ = 2.1 km s−1. The derived effective temperature here of 3825 K is very close to Smith & Lambert, as is the microturbulence of 2.2 km s−1. The surface gravity from Smith & Lambert is larger by 0.7 dex, but their value relied on the distance (as it does in this paper), for which they used a luminosity calibration for the Ca ii K-line reversal that was tied to a pre-Hipparcos Hyades distance that was too small. Making the estimated log g corrections to the Smith & Lambert (1985) abundances (based on their published sensitivities of abundance with log g) for 12C, 14N, 16O, Ti, Fe, and Ni we find a mean difference and standard deviation of ΔA(x) (this paper; Smith & Lambert) = +0.02 ± 0.14 dex.

Beta And was also used as a comparison star by Chou et al. (2010) in their study of Sgr dwarf-galaxy tidal streams in the Galactic halo. Chou et al. (2010) applied a similar technique as that employed here to estimate temperatures and gravities and found Teff = 3850 ± 75 K, log g = 0.9, and ξ = 1.96 km s−1, giving differences of only 25 K in Teff, 0.0 in log g, and +0.24 km s−1 in ξ. Their analysis used the spectral region around 7500 Å containing a sample of Fe i and Ti i lines and determined A(Fe)=7.12 ± 0.06 and A(Ti) = 4.55 ± 0.12, which can be compared to values of 7.23 ± 0.03 and 4.72 ± 0.02, respectively, here, yielding differences of +0.11 dex in Fe and +0.17 dex in Ti.

4.5.4. HD 199799

HD 199799 is classified as spectral type M2S, which means that it displayed enhanced bands of ZrO in classification spectra. Smith & Lambert (1988) detected the radioactive s-process element Tc I in its spectrum, identifying this star as an intrinsic thermally pulsing asymptotic giant branch (TP-AGB) star undergoing third dredge-up. A quantitative abundance analysis was then carried out by Smith & Lambert (1990) and they derived similar stellar parameters to those derived here: the same Teff = 3400 K, with a slightly lower value of log g = +0.3, compared to +0.5 here. The microturbulent velocity determined here from the H-band spectra is ξ = 2.4 km s−1, which is 0.7 km s−1 larger than Smith & Lambert (1990), who used Fe i lines near 7500 Å. Taking into account the slightly lower gravity and smaller ξ, the abundances from Smith & Lambert (1990) can be compared to those here for Fe, 12C, 14N, 16O, Ti, Cr, Co, and Ni, with the mean difference being ΔA(x)(Us—SL90) = −0.08 ± 0.13. The respective carbon isotopic ratios compare well, with Smith & Lambert (1990) finding 12C/13C = 28, while a value of 27 is obtained here. The derived C/O ratios are the same in the two studies, with 12C/16O = 0.68.

The comparisons presented here are between studies that have possibly used somewhat different families of model atmospheres (e.g., MARCS or ATLAS), different spectral regions, and linelists derived from heterogeneous sources. Given the variety of data and analysis techniques, the derived abundances, in all cases, compare well at approximately the 0.10–0.15 dex level. A glance at Table 7 reveals that such differences fall within this range of uncertainty for most elements given modest changes in fundamental model atmosphere parameters. Since the derived stellar parameters themselves rely to some degree on the type of analysis, the generally good agreement between the derived abundances from the various published sources, along with the expected uncertainties in this study, as presented in Table 7, is encouraging. The gist of the comparisons, since many of the studies are based on optical spectra, is that H-band spectroscopy can be used to define a precise internal abundance scale that can also be extended to abundances derived from other, more common spectral windows, such as in the optical.

5. DISCUSSION

The abundances from Table 6 and an estimation of their respective uncertainties from Table 7, along with the stellar parameters from Table 2 represent, to some extent, what ASPCAP will produce from the APOGEE spectra. The following discussion highlights these abundances in light of internal red giant stellar evolution, as well as Galactic chemical evolution and stellar populations. These discussion topics anticipate what sorts of information will be gleaned from the APOGEE targets.

5.1. CNO Abundances and Red Giant Evolution

Based on their luminosities and effective temperatures (Figure 2), four of the red giants studied here are either first-ascent giants or, perhaps for μ Leo, a clump red giant. HD 199799 is in a more advanced stage of stellar evolution where it is undergoing thermal pulses and third dredge-up on the AGB. The combination of nucleosynthesis and convective mixing in evolved red giants should result in surface abundances of 12C, 13C, and 14N that have been altered from their initial (main sequence) values due to CN-cycle H-burning and first dredge-up.

The result of CN-cycle mixing is seen mostly easily in the 12C/13C ratios which are much higher in main-sequence stars, with the Sun having a value of 89. All four of the RGB or clump giants have quite low values, ranging from 12C/13 C ∼ 6–20. Gilroy (1989) and many subsequent studies of both open and globular clusters (e.g., Mikolaitis et al. 2012 or Briley et al. 1995) have noted that the carbon isotopic ratios correlate with stellar mass in stars above the luminosity bump on the RGB or in the clump; red giants more massive than M ∼ 2.5 M have values of 12C/13C ∼ 22–30, while lower-mass red giants exhibit decreasing isotopic ratios with declining mass. The mass estimates derived here for field stars from distance–luminosity and stellar evolutionary tracks are modestly well defined and the top panel of Figure 5 shows the 12C/13C values versus the estimated stellar mass. Also plotted are representative results for open clusters from Gilroy (1989) and Mikolaitis et al. (2012), and the globular cluster M71 from Briley et al. (1995). The cluster red giants have recently evolved from the cluster turn-offs, which have a well-determined mass. For the lower-mass stars included in these cluster studies, evolution along the RGB is rapid (∼107 yr) when compared to main-sequence lifetimes (∼108–109 yr) so the red giant mass is similar to the turn-off mass. This can be demonstrated using the evolutionary tracks from Bertelli et al. (2009), where an isochrone with a turn-off mass of 4.0 M has a red giant tip mass of 4.2 M, or a ratio of RGB to turn-off (TO) masses of 1.05. This ratio of masses becomes even closer to 1.0 for lower masses, which include all of the points in Figure 5. The points from Gilroy represent the mean values of 12C/13C for 19 open clusters with a total of 55 red giants in her sample. The mean carbon isotopic ratios for two other open clusters (with a total of 10 red giants from Cr261 or NGC 6253) from Mikolaitis et al. (2012) are also shown (open blue squares). A representative globular cluster, M71, which has a lower main-sequence mass (M ∼ 0.9 M) from Briley et al. (1995) is plotted as the open blue triangle and represents the mean of two of the CN-weak stars in this cluster. The CN-strong giants that Briley et al. also studied are not plotted, as they exhibit slightly lower isotopic ratios but almost certainly represent second-generation stars in M71 that were formed from large fractions of the ejecta of first-generation intermediate-mass red giants and thus are not representative of single-star red giant evolution.

Figure 5.

Figure 5. Behavior of the nuclei 12C, 13C, and 14N with stellar mass on the RGB as characterized by 12C/13C and C/N (where C = 12C+13C and N = 14N). The relative abundances of these nuclei are altered by H-burning via the CN-cycle and convective mixing on the RGB dredges up the processed abundances. The TP-AGB star HD 199799 is identified, as it has significantly increased its 12C abundance from He-burning and third dredge-up and has thus increased its 12C/13C relative to when it was on the RGB.

Standard image High-resolution image

The curves in the top panel of Figure 5 represent stellar evolutionary model results from Charbonnel & Lagarde (2010); the solid black line illustrates their predictions for standard RGB dredge-up (no "extra mixing" or rotation included). Clearly, standard RGB first dredge-up does not fit the observed trend of 12C/13C for red giants with M ⩽ 2.2 M. This poor comparison between theory and observation has led to work in trying to identify other types of "extra mixing" mechanisms that would lead to lower carbon isotopic ratios in the lower-mass red giants. Two such hypothesized mechanisms are thermohaline mixing or thermohaline mixing coupled with rotational mixing, as discussed recently by Charbonnel & Lagarde (2010). The dashed curves present results from Charbonnel & Lagarde also, with the short dashes for thermohaline mixing only and the long dashes for thermohaline mixing plus rotational mixing for models with initial equatorial rotational velocities of 110 km s−1. The stellar models that incorporate these particular types of extra mixing provide fair agreement with the observations, although these are not the only types of processes that have been studied and are used here to merely illustrate the effect of extra mixing. The trends found in the clusters also overlap well with the field stars analyzed here, except for HD 199799, which is in a different evolutionary phase than the other red giants. This star has enriched its surface in 12C that was synthesized by He-burning as a TP-AGB star. In the simplest case of pure 12C dredge-up, the surface 12C/13C ratio would increase, which is where HD 199799 lies in the top panel of Figure 5 relative to the other red giants. Such an effect has been noted by both Lambert et al. (1986) and Smith & Lambert (1990), where these studies find increasing values of 12C/13C as TP-AGB stars evolve from M to MS to S to C stars (with increasing values of 12C/16O).

The bottom panel of Figure 5 plots the abundance by number ratio of N(C)/N(N) versus stellar mass, where the carbon abundance is now the total of 12C + 13C: for reference, the solar ratio is 3.80. Nitrogen is also the total abundance, however, in the case of the observations, 15N is ignored as its abundance is very low (and not detected in these red giants) when compared to 14N. There is an observed decrease in C/N with increasing mass as predicted by both standard RGB dredge-up and dredge-up plus thermohaline mixing (Charbonnel & Lagarde 2010), which are shown as the solid and dashed curves, respectively, and are from the same models as in the top panel. In the case of the C/N ratios, the total carbon and nitrogen abundances are not as strongly influenced by non-standard mixing as the carbon isotopic ratios.

Figure 6 combines the C/N ratios versus 12C/13C ratios for both the observed red giants and the models from Charbonnel & Lagarde (2010). When displayed this way, the different behavior between C/N versus 12C/13C is striking and these two values are purely spectroscopically derived quantities. The trend of C/N with 12C/13 seems to provide a sensitive indicator of red giant mixing processes, as well as yields a mass estimate for stars on the RGB. Again, HD 199799 stands apart from the other red giants as it has internally enhanced its 12C abundance by about a factor of two since it was an RGB star, based on its 12C/13C and mass, relative to other RGB stars in the top panel of Figure 5. The dash-dotted line connects to the region in this diagram where HD 199799 would fall if its surface 12C abundance were reduced by a factor of two, which is a simplified sketch of how TP-AGB stars might evolve across such a diagram from the RGB.

Figure 6.

Figure 6. Observed and model predicted behavior of N(C)/N(N) vs. 12C/13C; the model curves show the so-called standard dredge-up (solid line curve) and a thermohaline mixing model as a function of stellar mass from Charbonnel & Lagarde (2010; dashed line curve). Standard dredge-up does not predict the observed behavior of N(C)/N(N) as a function of 12C/13C and some type of additional mixing process is needed. The TP-AGB star HD 199799 is identified here as it has dredged-up enough 12C to have approximately doubled its abundance relative to its abundance on the RGB; the dash-dotted line moves HD 199799 to its approximate position when it was evolving along the RGB.

Standard image High-resolution image

The assumption of only CN-cycle mixing taking place in the red giants is tested in Figure 7, where the total carbon abundance (12C + 13C) is plotted versus the nitrogen abundance (here 14N). The solar values of C and N are also shown, along with the dashed line illustrating scaled solar abundances for both C and N. The continuous blue curves represent decreasing carbon and increasing nitrogen abundances from the scaled solar lines, with the constraint that the number of C + N nuclei is conserved, as in the CN-cycle (with no leakage into the ON-cycle). Each curve is scaled in its initial C and N abundances by the stellar value of [Fe/H] for α Boo, β And, δ Oph, and μ Leo. CN-cycle mixing alone is sufficient to represent the derived abundances of C and N with initial abundances having [C + N/Fe] ∼ 0.0 for each RGB star. Due to having undergone a third dredge-up and 12C enrichment, HD 199799 will be offset from its CN curve, so no curve is plotted for it. The dash-dotted line simply shows the shift in A(C) if HD 199799 has doubled its RGB 12C abundance. Having nearly the same metallicity as β And, HD 199799 might be expected to have been near the CN curve for β And before it began its evolution as a TP-AGB star.

Figure 7.

Figure 7. Testing if total N(C) + N(N) is conserved, as would be expected for pure CN-cycle dredge-up. The dashed line represents metallicity scaled values of A(C) and A(N) with a solar C/N ratio. The blue curves map constant values of C+N scaled to the metallicity (taken to be [Fe/H]) of each star. Clearly, initial (main-sequence) values of [C+N/Fe] = 0, followed by simple CN-cycle dredge-up, are adequate representations of what is observed in these field RGB stars. The TP-AGB star HD 199799 has increased it 12C abundance by approximately 0.35 dex and this shift in A(C) is shown by the dash-dotted line.

Standard image High-resolution image

The APOGEE H-band wavelength region presents good opportunities for probing many nuclear and mixing processes that occur over the various phases of red giant evolution.

5.2. Element Ratios and Galactic Chemical Evolution

Other than the C and N isotopes, the vast majority of APOGEE red giant targets will not have altered their surface abundances of the other studied elements, so the remaining elements are dominated by chemical evolution within the various stellar populations of the Milky Way (or its accreted systems). The derived abundances of each of the other elements are discussed in the light of observed trends of these elements with A(Fe) from other studies in the literature. The number of published abundances is now enormous; however, in this first exploratory study using the APOGEE linelist, a comparison with a small number of selected literature sources is sufficient. For comparison purposes, in this paper a stellar sample representative of the Milky Way stellar populations was constructed using Reddy et al. (2003; thin disk with some thick disk and having 58 stars), Reddy et al. (2006; thick disk with 176 stars), Johnson (2002; halo with 23 stars), and Fulbright (2002; halo with 178 stars).

The investigation of elements derived from the APOGEE linelist is via a number of figures, beginning with Figure 8, which presents values of O/Fe, Mg/Fe, and Al/Fe versus A(Fe); the shorthand notation here of x/Fe denotes the logarithmic number ratio of element x to iron, log[N(x)/N(Fe)], or A(x)–A(Fe). This nomenclature is chosen over the [x/Fe] in the plots as it presents a direct logarithmic number ratio with the solar abundance ratio being noted in each panel, thus no reference to the underlying solar abundance is needed, but the relative scales are equivalent to [x/Fe]. The top panel of Figure 8 shows the behavior of O/Fe and, as has been well known for decades, values of O/Fe are elevated in the thick disk and halo stars. Among the field red giants here, the slightly metal-poor objects α Boo and β And both have O/Fe enhanced by ∼ +0.2–+0.3. The near-solar metallicity M-giant δ Oph has a near-solar ratio of O/Fe, while the metal-rich giant μ Leo also exhibits a near-solar value O/Fe.

Figure 8.

Figure 8. Abundances of oxygen, magnesium, and aluminum, relative to iron, shown as log(N(X)/N(H)) vs. A(Fe). The five stars studied here are represented by the filled red circles. The small blue symbols are results from Reddy et al. (2003, 2006) and the lower metallicity small magenta symbols are from Fulbright (2002) and Johnson (2002). The solar position is indicated with the dashed line showing the solar abundance ratio. The abundances derived here using the APOGEE linelist are well mapped into the general trends observed for Galactic field populations. The error bars are taken from the standard deviations in Table 6. Errors for the ratios are quadrature sums of the element in question and Fe i values of σ.

Standard image High-resolution image

The behavior of Mg/Fe is shown in the middle panel of Figure 8 and, like oxygen, Mg/Fe ratios are typically elevated in the thick disk and halo populations. The four red giants closest to solar in A(Fe) all exhibit basically solar Mg/Fe ratios, with α Boo having a slight enhancement of ∼ +0.1 dex, which overlaps the thick disk populations.

Aluminum abundances are plotted in the bottom panel of Figure 8 and the behavior of Al/Fe versus A(Fe) is not so well defined, with relatively large scatter. Mu Leo is found to be modestly enhanced in Al/Fe by about +0.10–+0.15 dex, as is α Boo, with an enhancement of ∼+0.2.

The abundances for the rest of the elements studied here are illustrated in Figure 9 (for Si, K, and Ca), Figure 10 (for Ti, V, and Cr), Figure 11 (for Mn, Co, and Ni), and Figure 12 (for Cu). In all elements, the general trends found in the general Galactic populations are recovered in the small sample of field red giants included here. The one exception is Ti in δ Oph and μ Leo, both of which are found to exhibit enhanced values of Ti/Fe of +0.2 dex. There is no obvious explanation for this, however, the continuing analysis of APOGEE spectra with this linelist will reveal if this possible effect appears in significant numbers of other red giants. The results presented and discussed here indicate that H-band spectra analyzed in LTE with the APOGEE linelist and 1-D model atmospheres can provide accurate results for chemical abundances.

Figure 9.

Figure 9. Abundances of silicon, potassium, and calcium, relative to iron, vs. A(Fe). The plotting format and symbols are the same as in Figure 8.

Standard image High-resolution image
Figure 10.

Figure 10. Abundances of titanium, vanadium, and chromium, relative to iron, vs. A(Fe). The plotting format and symbols are the same as in Figure 8.

Standard image High-resolution image
Figure 11.

Figure 11. Abundances of manganese, cobalt, and nickel, relative to iron, vs. A(Fe). The plotting format and symbols are the same as in Figure 8.

Standard image High-resolution image
Figure 12.

Figure 12. Abundance of copper relative to iron vs. A(Fe). The plotting format and symbols are the same as in Figure 8.

Standard image High-resolution image

6. CONCLUSIONS

Archival high-resolution FTS spectra from 15100 to 16900 Å are analyzed to derive detailed chemical abundance distributions for five field red giants using the APOGEE linelist and model atmosphere grid. The red giants span a range in Teff from 3400 K to 4540 K and log g from 0.5 to 2.1, covering a significant part of the effective temperature–surface gravity range being observed in the APOGEE survey.

Chemical abundances are determined for 16 elements, along with the minor carbon isotope, 13C. Comparisons of the abundances derived here with those from previously published studies of four of the target stars here, mostly from visual-wavelength high-resolution spectra, find abundance differences of ∼0.1 dex. The APOGEE wavelength window with its derived linelist and model atmosphere grid, coupled to spectrum synthesis modeling, can be used to analyze high-resolution H-band spectra of red giants to probe both Galactic chemical evolution in stellar populations, as well as internal red giant nucleosynthesis and mixing.

We thank Ken Hinkle for helping to make the archival FTS spectral available and to Luan Ghezzi for useful comments on an initial draft of this paper. V.S. acknowledges partial support for this research from the National Science Foundation (AST1109888) and CAPES, Brazil. This publication makes use of data products from the Two Micron All Sky Survey, which is a joint project of the University of Massachusetts and the Infrared Processing and Analysis Center/California Institute of Technology, funded by the National Aeronautics and Space Administration and the National Science Foundation. This research has made use of the SIMBAD database, operated at CDS, Strasbourg, France.

Funding for SDSS-III has been provided by the Alfred P. Sloan Foundation, the Participating Institutions, the National Science Foundation, and the U.S. Department of Energy Office of Science. The SDSS-III Web site is http://www.sdss3.org/.

SDSS-III is managed by the Astrophysical Research Consortium for the Participating Institutions of the SDSS-III Collaboration including the University of Arizona, the Brazilian Participation Group, BrookhavenNational Laboratory, University of Cambridge, Carnegie Mellon University, University of Florida, the French Participation Group, the German Participation Group, Harvard University, the Instituto de Astrofisica de Canarias, the Michigan State/Notre Dame/JINA Participation Group, Johns Hopkins University, Lawrence Berkeley National Laboratory, Max Planck Institute for Astrophysics, New Mexico State University, New York University, Ohio State University, Pennsylvania State University, University of Portsmouth, Princeton University, the Spanish Participation Group, University of Tokyo, University of Utah, Vanderbilt University, University of Virginia, University of Washington, and Yale University.

Facility: Mayall (FTS) - Kitt Peak National Observatory's 4 meter Mayall Telescope

Please wait… references are loading.
10.1088/0004-637X/765/1/16