Arxiv:1106.5449V1 [Astro-Ph.SR] 27 Jun 2011 Aepsil Ytegnru Nnilspoto H .M Ekf Keck M
Total Page:16
File Type:pdf, Size:1020Kb
Carbon and Oxygen in Nearby Stars: Keys to Protoplanetary Disk Chemistry0 Erik A. Petigura Astronomy Department, University of California, Berkeley, CA 94720 [email protected] Geoffrey W. Marcy Astronomy Department, University of California, Berkeley, CA 94720 Received ; accepted arXiv:1106.5449v1 [astro-ph.SR] 27 Jun 2011 0Based in part on observations obtained at the W. M. Keck Observatory, which is operated as a scientific partnership among the California Institute of Technology, the University of California, and the National Aeronautics and Space Administration. The Observatory was made possible by the generous financial support of the W. M. Keck Foundation. –2– ABSTRACT We present carbon and oxygen abundances for 941 FGK stars—the largest such catalog to date. We find that planet-bearing systems are enriched in these elements. We self-consistently measure NC /NO, which is thought to play a key role in planet formation. We identify 46 stars with NC /NO 1.00 as poten- ≥ tial hosts of carbon-dominated exoplanets. We measure a downward trend in [O/Fe] versus [Fe/H] and find distinct trends in the thin and thick disk, sup- porting the work of Bensby et al. (2004). Finally, we measure sub-solar NC /NO +0.11 = 0.40 0.07 for WASP-12, a surprising result as this star is host to a transiting − hot Jupiter whose dayside atmosphere was recently reported to have NC /NO ≥ 1 by Madhusudhan et al. (2011). Our measurements are based on 15,000 high signal-to-noise spectra taken with the Keck 1 telescope as part of the California Planet Search. We derive abundances from the [OI] and CI absorption lines at λ = 6300 and 6587 A˚ using the SME spectral synthesizer. Subject headings: catalogs — stars: abundances — stars: fundamental parameters — techniques: spectroscopic — planetary systems –3– 1. Introduction After primordial hydrogen and helium, carbon and oxygen are the most abundant elements in the cosmos. Life on earth is built upon the versatility of carbon’s four valence electrons and is powered by metabolizing nutrients with oxygen. The prevalence of carbon and oxygen gives them a prominent role in stellar interiors, opacities, and energy generation. As a result, studying their abundances helps to reveal the nucleosynthetic chemical evolution of galaxies. The interstellar medium is thought to be enriched with oxygen by Type II supernovae. Taken with iron, which is produced in both Type Ia and Type II supernovae, oxygen provides a record of galactic chemical enrichment and star formation rate (Bensby et al. 2004). It is well known that stars synthesize helium into carbon through the triple alpha reaction. However, it is still unclear which stars dominate carbon production in the galaxy. For a discussion of the possible sites of carbon synthesis see Gustafsson et al. (1999). The ratio of carbon to oxygen (NC /NO) is thought to play a critical role in the bulk properties of terrestrial extrasolar planets. Kuchner & Seager (2005) and Bond et al. (2010) predict that above a threshold ratio of NC /NO near unity, planets transition from silicate- to carbide-dominated compositions. We present the oxygen and carbon abundances derived from the [OI] line at 6300 A˚ and the CI line at 6587 A˚ for 941 stars in the California Planet Search (CPS) catalog. We compute the abundances with the Spectroscopy Made Easy (SME) spectral synthesizer (Valenti & Piskunov 1996). Using SME, we self-consistently account for the NiI contamination in [OI] and report detailed Monte Carlo-based errors. Others have measured stellar carbon and oxygen before. Edvardsson et al. (1993) measured oxygen in 189 F and G dwarfs, and Gustafsson et al. (1999) measured carbon in 80 of these stars. –4– More recent studies include, Bensby et al. (2005), Luck & Heiter (2006), and Ram´ırez et al. (2007). However, the shear number (15,000) of CPS spectra give us a unique opportunity to measure the distributions of these important elements in a large sample. 2. Observations 2.1. Stellar Sample The stellar sample is drawn from the Spectroscopic Properties of Cool Stars (SPOCS) catalog (Valenti & Fischer 2005, hereafter VF05) and from the N2K (“Next 2000”) sample (Fischer et al. 2005). We include 533 N2K stars and 537 VF05 stars for a total of 1070 stars. We adopt stellar atmospheric parameters for each star from VF05 and from the identical analysis for the N2K targets (D. Fischer 2008, private communication). These parameters are: effective temperature, Teff ; gravity, log g; metallicity, [M/H]; rotational broadening, v sin i; macroturbulent broadening, vmac; microturbulent broadening, vmic; and abundances of Na, Si, Ti, Fe, and Ni. Metallically includes all elements heavier than helium. A star’s abundance distribution is the solar abundance pattern from Grevesse & Sauval (1998) scaled by the star’s metallicity. Na, Si, Ti, Fe, and Ni abundances are computed independently from [M/H] and are allowed vary from scaled solar [M/H]. 2.2. Spectra Our spectra were taken with HIRES, the High Resolution Echelle Spectro- graph (Vogt et al. 1994) between August, 2004 and April, 2010 on the Keck 1 Telescope. The spectra were originally obtained by the CPS to detect exoplanets. For a more complete –5– description of the CPS and its goals, see Marcy et al. (2008). The CPS uses the same detector setup each observing run and employs the HIRES exposure meter (Kibrick et al. 2006) to set exposure times, ensuring consistent and high quality spectra across years of data collection. The spectra have resolution R = 50,000 and S/N 200 at 6300 and ∼ 6587 A.˚ This analysis deals with three classes of observations: 1. Iodine cell in. For the majority of its observations, the CPS passes starlight through an iodine cell (Marcy & Butler 1992), which imprints lines between 5000 and 6400 A˚ that serve as a wavelength fiducial. We discuss how we remove these lines and their effect on oxygen measurements in 3.4. § 2. Iodine cell out. Calibration spectra taken without the iodine cell. 3. Iodine reference. At the beginning and end of each observing night, the CPS takes reference spectra of the iodine cell using an incandescent lamp. 3. Spectroscopic Analysis 3.1. Line Synthesis We use the SME suite of routines to fine-tune line lists based on the solar spectrum, determine global stellar parameters, and measure carbon and oxygen. To generate a synthetic spectrum, SME first constructs a model atmosphere by interpolating between the Kurucz (1992) grid of model atmospheres. Then, SME solves the equations of radiative transfer assuming Local Thermodynamic Equilibrium (LTE). Finally, SME applies line-broadening to account for photospheric turbulence, stellar rotation, and instrument profile. For a more complete description of SME, please consult Valenti & Piskunov (1996) and VF05. We emphasize that SME solves molecular and ionization equilibrium for a core –6– group (around 400) of species that includes CO (N. Piskunov 2011, private communication). 3.2. Atomic Parameters Measuring stellar oxygen is notoriously difficult because of the limited number of indicators in visible wavelengths. The general consensus is that the weak, forbidden [OI] transition at 6300 A˚ is the best indicator because it is less sensitive to departures from local thermodynamic equilibrium than other indicators. In dwarf stars, this line suffers from a significant NiI blend, which is an isotopic splitting of 58Ni and 60Ni (Johansson et al. 2003). The NiI feature was first noted by Lambert (1978), but only recently included in abundance studies (Allende Prieto et al. 2001). Carbon is more generous to visual spectroscopists. We select the CI line at 6587 A˚ because it sits relatively far from neighboring lines and is in a wavelength region with weak iodine lines (see 3.4). § Line lists are initially drawn from the Vienna Astrophysics Line Database (Piskunov et al. 1995). We tune line parameters by fitting the disk-integrated National Solar Observatory (NSO) solar spectrum of Kurucz et al. (1984) with the SME model of the solar atmosphere. Table 1 lists the atmospheric parameters adopted when modeling the sun. We fit a broad spectral range from 6295 to 6305 A˚ surrounding the [OI] line and 6584 to 6591 A˚ surrounding the CI line. We adopt the solar abundances of Grevesse & Sauval (1998) except for O and 1 Ni where we adopt log ǫO = 8.70 and log ǫNi =6.17 (Scott et al. 2009) and log ǫC =8.50 (Caffau et al. 2010). We adjust line centers, van der Waals broadening parameters (Γ6), and oscillator strengths (log gf) so our synthetic spectra best match the NSO atlas. Table 2 shows the best fit atomic parameters after fitting the NSO solar atlas. Given the high quality of the solar spectrum, solar abundances and line parameters 1 log ǫX = log10(NX /NH )+12 –7– are often measured using sophisticated three-dimensional, hydrodynamical, non-LTE codes. For this work, however, we are more interested in self-consistently determining line parameters using SME than from a more sophisicated solar model. As a result, the line parameters in Table 2 are not in tight agreement with the best laboratory measurements. For example, Johansson et al. (2003) measured log gf = 2.11 for the NiI blend in contrast − to log gf = 1.98 in this work. The purpose of fitting the atomic parameters in the sun − is to determine the best parameters given our atmospheric code and our adopted solar abundance distribution. We show the fitted NSO spectrum for both wavelength regions in Figures 1 and 2. The shaded regions (6300.0-6300.6 and 6587.4-6587.8 A)˚ represent the fitting region. Only points in the fitting region are used in the χ2 minimization routines (see 3.5). § Figure 3 shows a close up view of the [OI]/NiI blend in the sun.