arxiv: v1 [astro-ph.ep] 7 Aug 2011

Transcription

1 Rocky Extrasolar Planetary Compositions Derived from Externally-Polluted White Dwarfs B. Klein 1, M. Jura 1, D. Koester 2, B. Zuckerman 1 ABSTRACT arxiv: v1 [astro-ph.ep] 7 Aug 2011 We report Keck High Resolution Echelle Spectrometer data and model atmosphere analysis of two helium-dominated white dwarfs, PG and HS , whose heavy pollutions most likely derive from the accretion of terrestrial-type planet(esimal)s. For each system, the minimum accreted mass is g, that of a large asteroid. In PG , Mg, Cr, Mn, Fe and Ni have abundance ratios similar to bulk Earth values, while we measure four refractory elements, Ca, Sc, Ti and V, all at a factor of 2 3 higher abundance than in the bulk Earth. For HS the swallowed material was compositionally similar to bulk Earth in being more than 85% by mass in the major element species, O, Mg, Si, and Fe, and with abundances in the distinctive proportions of mineral oxides compelling evidence for an origin in a rocky parent body. Including previous studies we now know of four heavily polluted white dwarfs where the measured oxygen and hydrogen are consistent with the view that the parents bodies formed with little ice, interior to any snow-line in their nebular environments. The growing handful of polluted white dwarf systems with comprehensive abundance measurements form a baseline for characterizing rocky exoplanet compositions that can be compared with bulk Earth. Subject headings: Planets and satellites: composition Stars: individual: (HS , PG , G 241-6) Stars: abundances white dwarfs 1. INTRODUCTION Accretion by white dwarfs of tidally disrupted bodies from their planetary systems (Jura 2003) has become the most viable explanation for the presence of high-z (atomic number Z>2) elements in the vast majority of polluted DAZ, DBZ, and DZ white dwarf atmospheres with T eff < 20,000 K (Jura 2008; Zuckerman et al. 2010; Farihi et al. 2010a, and references therein). The basic idea is that smaller planetary bodies such as asteroids or dwarf planets, can acquire highly eccentric orbits by perturbations from larger planetary bodies, especially considering the dynamical rearrangement of the planetary system during stellar evolution and mass loss (Debes & Sigurdsson 2002). Bonsor et al. (2011) demonstrate that the post-stellar mass loss evolution of a planetesimal belt can supply the 1 Department of Physics and Astronomy, University of California, Los Angeles, CA , USA; kleinb, jura, ben@astro.ucla.edu 2 Institut fur Theoretische Physik und Astrophysik, University of Kiel, Kiel, Germany; koester@astrophysik.uni-kiel.de

2 2 required mass influx to explain observed polluted white dwarf (WD) accretion rates. Objects that journey within the Roche radius of the WD will experience tidal shredding, forming a disk of dust and/or gas (Jura 2003, 2008), and subsequently accrete into the star s, otherwise pure, hydrogen and/or helium atmosphere. This natural mechanism offers the potential to measure the distilled elemental constituents of a planetary parent body or bodies a unique and powerful tool, indeed. High-resolution and high-sensitivity spectroscopy, along with state-of-the-art model atmospheres (Koester 2009, 2010; Tremblay & Bergeron 2009; Dufour et al. 2010; Vennes et al. 2010, 2011), has opened the door to detailed studies of WD atmospheres and the high-z material that pollutes them. To date, a handful of WDs with very rich spectra have been studied at high-resolution and high-sensitivity, beginning with the Keck/HIRES follow-up (Zuckerman et al. 2007) of GD 362, whose extreme high-z-rich atmosphere was pointed out by Gianninas et al. (2004) and Kawka & Vennes (2006). The detection of oxygen (O) along with the other major elements (Mg, Si, and Fe) in GD40 s polluting parent body (Klein et al. 2010) established chemical evidence the stoichiometric balance of mineral oxides for rocky extrasolar planetesimals. Ongoing discoveries of polluted WDs displaying all major terrestrial elements (Vennes et al. 2010; Dufour et al. 2010; Zuckerman et al. 2010; Farihi et al. 2011) signify that externally polluted white dwarfs have entered a stage of being practical measuring devices of the elemental compositions of rocky exoplanets. Bond et al. (2010) predict a wide diversity of terrestrial exoplanets including Earth-like, refractory-rich, and carbon-rich objects. Polluted white dwarfs provide real data at the level of detail needed to make discerning measurements (e.g., Zuckerman et al. 2011). In this paper we report the detailed analysis of HIRES spectra of HS (V mag =16.1, DBAZ, T eff 14,400 K; hereafter HS2253) and PG ( K , V mag =14.8, DZAB, T eff 10,800 K; hereafter PG1225). Both stars were previously known to display features of H, Mg, Si, Ca and Fe mainly from low ( 6Å) resolution IUE and optical data (Koester et al. 1990; Friedrich et al. 1999; Wolff et al. 2002). Calcium in PG1225 was also (re)observed and analyzed from the ESO SN Ia progenitor (SPY) survey (Koester et al. 2005b). The HIRES observations provide a new and important detection of oxygen in HS2253, completing the set of major terrestrial elements O, Mg, Si, and Fe as well as H, Ca, Ti, Cr, Mn and possibly, Ni. For PG1225 the HIRES data display six newly identified elements: Sc, Ti, V, Cr, Mn, and Ni. By analyzing all the detected elements together in a given model, with multiple lines for each element, and accounting for abundance sensitivity to errors of the model parameters, we calculate pollutant-to-pollutant abundance ratios with an estimated uncertainty of 40% or less, a significant improvement over previous analyses in which derived abundances relative to helium are typically uncertain by dex. Infrared studies show that HS2253 does not have an infrared excess (Farihi et al. 2009), and PG1225 has just a subtle infrared excess appearing only at 7.9µm (Farihi et al. 2010b; Kilic et al. 2008a). Notwithstanding the lack of detected circumstellar dust, we demonstrate that these systems have accreted an appreciable amount of material with terrestrial planet(esimal) composition. We also report the detection of Ca II H&K emission cores in PG1225, the first time seen in a white dwarf. Other peculiar line profiles observed with HIRES are narrow absorption cores in helium lines of PG1225, GD362 and GD16, which are displayed in this paper.

3 3 There are a number of benefits of working with helium-dominated WDs. One is that, compared to hydrogendominated WD atmospheres, the lower opacity helium environments more readily display tiny concentrations of numerous contaminant elements, e.g. the trace elements vanadium and scandium are detected in PG1225 at orders of magnitude lower abundance than the helium background they are diluted in. Secondly, detailed model atmosphere analysis of helium-atmosphere WDs indicates the hydrogen in these stars must be accreted (Voss et al. 2007; Dufour et al. 2007); thus when helium is the background constituent, one obtains a measurement of the polluting hydrogen abundance, which is important for understanding the potential accretion of water ice (Jura et al. 2009; Jura & Xu 2010). A simplifying aspect of the helium-atmosphere WDs is that in our T eff range of interest (< 20,000 K), all have outer convection zones where the elements are homogeneously mixed. Finally, another benefit of working with helium-rich WDs is that with relatively large convective zones (relatively long settling times), the pollutants are not disappearing so fast from our view. Thus, a higher fraction of polluted WDs can be found (e.g., Zuckerman et al. 2010), and meaningful lower limits on accreted masses can be obtained. Recently, Dufour et al. (2010) and Koester et al. (2011) have established that absolute minimum accreted masses for a number of polluted WDs are g, close to that of Ceres, the largest asteroid in our own Solar system. Gravitational settling times are always less than the age of the WD, but in helium-atmosphere stars of T eff < 20,000K, the settling (diffusion) time scales for most elements are years (Koester 2009). Thus a measurement of atmospheric abundances in a helium-dominated WD is a snapshot of one stage of the evolving interplay of accretion and diffusion in that system. As was found for GD40 (Klein et al. 2010), the patterns of pollutant abundances themselves, can narrow the range of possibilities for the accretion diffusion state of a WD, and we utilize some such constraints in the present work. The article is laid out as follows: observations, data reduction, and measurements of the spectra are described in Section 2. In Section 3 we discuss effective stellar temperature, surface gravity, and the synthetic model atmospheres from which we derive the polluting element abundances presented in Section 4. We describe the effects of accretion and diffusion in Section 5, incorporate that into our results including results of previous studies in Section 6, and give our conclusions in Section OBSERVATIONS & MEASUREMENTS The data were acquired with the High Resolution Echelle Spectrometer (HIRES, Vogt et al. 1994) on the Keck I telescope of Mauna Kea Observatory. Table 1 is the observation log. Typically, a slit width was used, providing a spectral resolution of R = λ/ λ 40,000. To compensate for inferior seeing on 2008 Nov 15, a slit was used, resulting in a resolution of R 28,000 for those spectra. Before combining exposures of the same star with differing setups, the higher resolution spectra were degraded by Gaussian convolution to approximate the lower resolution spectra. Data reduction and absorption line measurements were carried out using PyRAF in a similar manner as described

4 4 forgd40 inklein et al. (2010). For thepresentworkwetook someextrasteps, sinceinadditiontoafullycontinuumnormalized version of each spectrum, we also wanted a version that preserved broad-line (> 1-2 Å) features. For this purpose we used observations of calibration white dwarfs (BD , G191-B2B, or Feige 34) in the fluxcalibration procedure to align the echelle orders (however, absolute flux calibration is not done). As is known for HIRES data (Suzuki et al. 2003), echelle artifacts irregularly appear as bends or waves in the continuum. We visually inspected and compared regions of similar wavelength coverage, and selectively combined only those portions of the spectra that appeared not to suffer significant distortions. Mostly this amounted to omitting some echelle-order edges. Figures 1 and 2 show the high-quality output, and success with broad lines, from these procedures. Residual artifacts are low-level and fairly broad, which for example can be seen as waves in the continuum in the region λ å of Figure 2. For HS2253 the continuum signal-to-noise (SNR) of the data ranges from for wavelengths in Å, and is at wavelengths shorter and longer than that. The spectrum of PG1225 is very high quality with the SNR of the unsmoothed continuum ranging from for wavelengths in Å and at wavelengths shorter and longer than that. The line identifications and equivalent widths are listed in Table 2, and were carried out in a similar manner as described for GD40 in Klein et al. (2010). For PG1225 the many broad lines and high number density of lines was not amenable to our usual method of measurement error estimation. Instead, based on our experience with HIRES data and checks over multiple exposures, for PG1225 we set the EW measurement errors to 20% for most lines, and 30% or more for a few lines where the continuum or interval surrounding the line is especially uncertain. We have also included here the details of the line measurements of G241-6 (SDSS g =15.5 mag, DBZ, T eff =15,300 K, log g=8.0), whose high-z rich nature and pollutant abundances were reported by Zuckerman et al. (2010). Especially in the case of the very rich spectrum of PG1225, there are a number of blended lines, which are indicated with a b superscript in Table 2. In a few regions, such as the Mg I λ3835 triplet, instead of taking the EW from fitting a Lorentzian or Voigt profile, we report the EW from a numerical sum of the flux under the continuum for a given interval surrounding the features. Owing to the excellent SNR of these spectra, the measured line strengths span a remarkably large range a factor of 3000 from Å (V II), to 1 Å (Ti II) to greater than 30 Å (Ca II H&K blend). All observed absorption lines, and those not observed, are consistent with expectations of line strengths for our derived abundances and SNRs of the spectra. Displays of various portions of the HIRES spectra are given in Figures For each transition, the shift of line center from the laboratory wavelength is measured, and converted to a heliocentric velocity. We use the average line velocity for a systemic value, and Doppler shift the models to the measured frame of the WD, which are plotted in red on the figures of spectra. The net systemic heliocentric velocities, including the stars gravitational redshifts, are 4 ± 4 km s 1 for HS2253 and 49 ± 3 km s 1 for PG1225. From the histogram of line velocity measurements shown in Figure 14, for PG1225 we found a range of outliers (filled bars in the figure) that represent a detectable shift of the set of Fe I lines (see also Figure 12). These are not of an interstellar origin, since as indicated in Table 2, most transitions arise from excited atomic levels. An analogous

5 5 shift for Si II lines has been examined and modeled with Stark shifts by Vennes et al. (2011), thus it is conceivable that a similar process is responsible for the effect on Fe I in PG1225; further investigation remains for future work. HS2253 s velocity measurements do not indicate any distinctive sub-populations. With such high-sensitivity observations of PG1225 a new feature has turned up (literally) that we have not seen before in WDs emission cores in the Ca II H&K lines as shown in Figure 15. The emission features are seen in all 3 separate exposures covering the Ca H&K region of PG1225, with no noticeable variations between exposures. Emission cores have been seen in Hα of hydrogen-dominated white dwarfs hotter than 25,000K, and are understood to arise from non-local-thermodynamic-equilibrium (NLTE) effects. The one HIRES exposure covering Hα in PG1225, taken about 2 weeks later, does not display clear Hα core emission with a continuum SNR 50 near 6650Å. Still, NLTE effects could conceivably cause emission in the strong Ca H&K lines since the cores are formed high up in the atmosphere, while the relatively weaker Hα line of PG1225 is likely formed at a deeper location. Alternatively, Ca H&K emission cores, as well as Hα core emission, are well known signatures of chromospheric activity in main sequence stars; further discussion of chromospheric activity in a WD is beyond the scope of this work (see Musielak et al. 2005). One more spectroscopic peculiarity we point out is the presence of a narrow absorption core component at the line center of He I λ5876 in PG1225, as well as in GD362 and GD16 ( HS ) shown in Figure 16. These three WDs have helium-dominated atmospheres, with effective temperatures in the range 10,500 K 11,500 K; they all have detected circumstellar dust (Becklin et al. 2005; Kilic et al. 2005; Farihi et al. 2009), heavy atmospheric pollution (Koester et al. 2005a; Zuckerman et al. 2007), and possess relatively large amounts of hydrogen in their convective envelopes. Two possibly analogous, but weakly detected, features in PG1225 are narrow absorption lines instead of broad ones centered on He I λ3889 and He I λ4472, with the latter shown in Figure 8. The He I λ5876 core feature was first discovered in HIRES data of GD362 by Zuckerman et al. (2007) and noted again by Tremblay et al. (2010); we now see it in HIRES spectra of GD16 and PG1225 as well. Zuckerman et al. (2007) suggested that NLTE effects, similar to those seen in the cores of Hα lines in hydrogen-dominated white dwarfs, could be present. The exact origin of the He core feature is not known at this time and we do not discuss it further in this paper. 3. MODEL ATMOSPHERES The model input physics is described by Koester (2010) with all absorption line data taken from the Vienna Atomic Line Database (VALD). Key inputparameters are the effective temperature, T eff, surface gravity, g, and the atmospheric abundances of elements relative to the dominant constituent, in the current study, helium. In the range of temperatures considered here, helium is mostly neutral, and the high-z elements and hydrogen are significant contributors of electrons in the atmosphere. Thus, we calculated and analyzed our models with hydrogen and high-z elements in a self-consistent way, and included Earth-like abundances for undetected but expected elements (such as O in PG1225). More than 2000 continuous absorption coefficients were incorporated to account for the blanketing

6 6 effect. In tuning the model inputs to reproduce an observed WD spectrum, it is desirable to achieve consistency for T eff, log g and abundance estimates that derive from 1) photometric fitting of the spectral energy distribution (SED), preferably using UV to optical and infrared wavelengths, 2) fitting the hydrogen and helium lines in profile and/or strength, and 3) demanding a balance in the abundances that come from different ionization states of the same element. In principle, the model would meet the requirements of 1) through 3) while also achieving a reasonable fit of the broadening and symmetry of heavy element lines in practice, it can be challenging to reconcile all these points, but we find reasonable consistency for the WDs studied here. Precise parameters are not required for our abundance analysis, since we are interested in the pollutant-to-pollutant ratios, which are relatively insensitive to the particular model (Klein et al. 2010, also Section 4 of this paper). In the analysis of helium dominated atmospheres with effective temperature below 16,000 K, we must deal with the difficulty that the parameters, T eff and gravity, g, are effectively coupled, such that T eff /log g fits to the helium lines are not unique. This has led many authors (present authors included) to set log g to the average value of 8.0, so that one may proceed with a fit to the temperature. In the current analysis we took a different approach, similar to Dufour et al. (2010), who found that a larger than average gravity was required to fit simultaneously the features from two ionization states, Mg I and Mg II, in the extremely polluted helium-rich WD, SDSS J (hereafter, SDSS0738). For each star in the present study we allow for a departure from the average gravity value, and instead work with a set of 3 T eff /log g solutions for T eff and log g which do a reasonably good job of satisfying the spectral energy distribution (SED), He I line strengths, and the ionization balance between neutral and singly ionized species of Fe, Mg and Ca. The model parameters used are listed in Table 3 and span a range of ±300 or 400 K for effective temperature and ±0.3 dex for log g. An additional constraint on the gravity could come from a sufficiently precise parallax measurement e.g. as Kilic et al. (2008b) carried out for GD 362 and help narrow the range of possible log g and T eff values. To our knowledge, PG and HS do not have reported parallaxes at this time HS T eff /log g The SEDs of our current models for HS2253 in the range T eff = K, scaled to 2MASS photometry, are in agreement with the IUE data and GALEX photometry (Figure 17). As Limoges & Bergeron (2010) pointed out, the strengths of some He I lines are especially gravity sensitive, and we make use of He I λ3820, λ4388, λ4713 and λ4922, in that regard. We find the set of models: T eff /log g = 14000/8.1, 14400/8.4 and 14800/8.7, all give reasonable and roughly similar fits for most of the He I lines, shown plotted by the red line in Figure 2. It turns 3 While G241-6 does have a parallax measurement (van Altena et al. 1995), it is not precise enough to improve on the spectroscopically assumed T eff /log g used by Zuckerman et al. (2010).

7 7 out that the over-prediction of He I λ4713 mentioned by Friedrich et al. (1999) is resolved by our use of higher than average gravity, although there remains an unaccounted-for shift of that line s center (top panel Figure 2). Another reason weprefer modelswith ahigher than average gravity, is wenoticed that at thevalue log g=8.0, for a given element, even though we could get a reasonably good fit to the equivalent widths with the same abundance, for almost all lines the model profiles are too deep and too narrow; this appeared to be the case for most elements. Having worked with matching similar models to similar HIRES observations of similar WDs (GD40 and G241-6), we presume that a mismatch in high-z line broadening between the data and instrument-convolved model, is not a problem of accounting for the instrumental broadening, but rather is indication that there is some intrinsic broadening needed in the model, such as a higher gravity. We found that the models listed in Table 3 produced better overall fits to the high-z line profiles than models of average gravity. Similarly, the fit to the Hα profile (Figure 3) is better with higher gravity. Finally, with our choice models, the ionization balance for abundances of Mg and Fe naturally come out of the analysis of HS2253 at our desired consistency of within 0.1 dex between abundances from neutral and singly ionized states (Table 4). Obtaining agreement in the ionization balance has proven challenging in other stars, including GD362, GD40 and PG1225 as described below PG T eff /log g At temperatures <12,000 K, the optical He I lines begin to disappear; He I λ5876 (Figure 8) is the strongest line observed. In this range of T eff, only a very small fraction of the photospheric helium is involved in the formation of optical helium lines, since they arise from levels of high excitation energy, e.g ev above ground. Small changes in the atmosphere may have a large impact on the He line appearance, and we know that there is some unmodeled activity in this stars upper atmosphere as is apparent from the emission cores of the Ca II H&K lines (Figure 15) and the narrow core of He I λ5876 (Figures 8 & 16). For the discussion of atmospheric parameters that follows, we ignore the narrow core of the He I λ5876 feature described in Section 2, and only consider the broad component of He I λ5876. The SEDs of our current models for PG1225, scaled to 2MASS photometry, are in good agreement with the IUE spectrum and GALEX photometry (upper panel Figure 17) for T eff near K, as used in previous work by Wolff et al. (2002) withlog g =8.0. However, looking at theionization balance, a10500/8.0 modelresults ina0.2dex difference between iron abundances derived from Fe I and Fe II. This is in the sense that Fe I lines are over-predicted and the discrepancy increases to 0.4 dex with a higher gravity model at 10500/8.3. Inversely, we find at T eff =10500 K, the two states of Fe are brought into marginal agreement (within 0.03 dex uncertainties Table 5) by lowering the gravity to log g=7.7. At this point in the T eff /logg parameter space, the 10500/7.7 model produces a He I λ5876 line with a predicted EW that matches the measured value at 1.5 Å, but the predicted profile is somewhat narrower and deeper than observed (red curve in Figure 8). Also, the He I λ4472 is not matched. The model helium lines can

8 8 be made shallower by increasing the gravity or decreasing the temperature, but always at the expense of degrading the ionization balance and the EW fit. As was the case for HS2253, there is a degeneracy of model solutions that fit the equivalent width of the helium line for example 10500/7.7, 10800/8.0, and 11100/8.3 all produce a similar strength (and depth) He I 5876 feature. Accepting the imperfections in the helium line profile fitting, we favor these models as our best compromise since they: 1) are in reasonable agreement with the UV/optical/IR spectral energy distribution, 2) fit the He I λ5876 equivalent width, and 3) support Fe I/Fe II abundance agreement to within 0.15 dex. We use this set of three models for establishing pollutant-to-pollutant abundance ratios, which we show in Section 4, are not very sensitive to the particular choice of model. 4. ABUNDANCES Abundances are calculated for each line by linearly scaling the model input abundance by the ratio of EWs measured in the observed spectrum (Table 2) to that of the model. In cases of blended features or otherwise noisy regions, a visual inspection of the profile fit aids the determinations. For each element we use the average of the set of calculated line abundances as the input for the next iteration of the model. It takes a few iterations until the average converges. Individual abundances are given equal weight in the final average element abundance. The EW measurement uncertainties are propagated through to an uncertainty in the final average abundance. We compare that to the mean deviation for the set of lines of an element, and assume the larger error. A challenge to abundance accuracy can come from discrepancies in the atomic line data and broadening parameters (Vennes et al. 2011). However, with multiple lines of each element observed, biases due to errors in the atomic data are minimized. We typically find a low dispersion in derived abundances for elements with multiple line detections (ranging from dex) generally comparable to, or smaller than, the propagated measurement uncertainties suggesting that the errors due to line data (including broadening parameters) are not a big factor in the current results. Tables 4 and 5 list the pollutant abundances relative to helium, derived from the different models considered for each star. This explicitly shows that the choice of model temperature and gravity are a significant source of systematic uncertainty in absolute abundance determinations. The absolute (pollutant-to-helium) abundances vary by dex with the modest range of model changes we apply. Also, estimates of the minimum accreted mass span a factor of ten, due to the fact that the calculated size of the convective zone has a strong dependence on the choice of model parameters as shown in Table 3. Fortunately for our goals, the relative (pollutant-to-pollutant) element abundances are rather homologous from one model to another with ratios varying only 40% or less, a comparable contribution to the uncertainty as from the propagated measurement error. The upshot is that the relative composition of the polluting elements can be measured to a greater precision than that of the absolute mass, or rate, of material accreted. The relative element abundances with respect to the fiducial magnesium are given in Table 6.

9 9 Comparisons of abundances with prior work is possible for H, Mg, Si, Ca and Fe, although previous uncertainties are typically large ( dex), and tend to propagate to a sizable range for the abundance ratios. For PG1225, the HIRES abundances relative to helium are nearly all within the uncertainties of those quoted by Wolff et al. (2002). Within the uncertainty limits, pollutant-to-pollutant ratios agree between previous and present analyses. From the HIRES data of PG1225, we obtain abundances for H, Mg, Ca, Sc, Ti, V, Cr, Mn, Fe, and Ni, and upper limits on O, Na, Al, and Sr, with an upper limit for Si that is consistent with its UV detection and analysis (Wolff et al. 2002); an upper limit for carbon is taken from the UV work (Wolff et al. 2002). For HS2253, the nominal HIRES abundances with respect to helium are all lower than those quoted by Friedrich et al. (1999) by at least 0.6 dex, most likely due to the differences in T eff /logg model inputs, particularly the much larger gravity used in the present work. Nonetheless, the pollutant-to-pollutant abundances derived from HIRES data and Friedrich et al. (1999) are mostly in agreement, with the exception of the Fe/Mg ratio. Previously, HS2253 appeared iron-rich (by number, Fe/Mg = 2, with 1.2 to 3.3 covering the range of uncertainty); the present HIRES analysis yields a ratio that is nominally half as large (Fe/Mg = 0.89 ± 0.18), and similar to Solar and Bulk Earth values. We obtain abundances in HS2253 from the HIRES data for H, O, Mg, Si, Ca, Ti, Cr, Mn, Fe, and a marginal detection of Ni; we derive upper limits on Na, Al, Sc, V, and Sr, while an upper limit for carbon is taken from the UV work (Friedrich et al. 1999). 5. ACCRETION DIFFUSION MODELS As mentioned in Section 1, a helium-dominated polluted WD system has an evolving interplay of accretion and diffusion, since these WDs have relatively long settling times. Calculated settling times specific to our particular models of HS2253 and PG1225 for the various elements studied are given in Tables 4 and 5. As described in Koester (2009), the accretion diffusion situation is roughly categorized by 3 phases: 1) an early phase in which accretion began recently, that is, before differential settling has had time to significantly alter the abundance ratios the observed atmospheric abundances directly represent the accreted ones, 2) a steady state with ongoing accretion, where the observed abundances are related to the accreted ones by ratios of settling times, and 3) a dissipation phase after accretion has slowed, paused or ended, and elements with shorter settling times disappear faster than those with longer settling times. We use these benchmarks to characterize our observations. If accretion has been going on longer than a couple of settling times, then a steady state interpretation may be appropriate. This assumes a semi-continuous flow of material lasting for 1 Myr or more in the systems under study here. To calculate the steady state abundances, for each element we take the average over different models of the ratios of settling times, and use the average ratios to apply a differential settling correction to the abundance ratios of Table 6 as prescribed by Equation 7 of Koester (2009). It is worth noting that analogous to the abundances, while the settling times themselves vary significantly from one model to another, the pollutant-to-pollutant ratios of settling times vary very little (RMS of most < 3%, all < 8%) with the model variations for a given source. For each element, the RMS s of diffusion data among models are propagated in quadrature into the steady state abundance

10 10 uncertainties. The principle result of applying a steady state correction is an increase in the abundance of iron and other heavy elements relative to the lighter elements in the inferred parent body. This effect is apparent in Figure 18 where the early phase ( convective zone = CVZ) abundances are plotted as filled markers, while the the open markers represent parent body abundances for the system in a steady state. For the set of elements plotted in Figure 18, Mg and Si are the lightest and have the longest settling times. Since Mg is the fiducial, applying a settling correction moves all the points of the plot up, relative to Mg, as seen in going from early phase to steady state values. The difference between an early phase and a steady state interpretation is most pronounced in the O/Fe ratio, which can be seen by comparison of the bulk composition bar graphs, Figures 19 and RESULTS AND DISCUSSION From Table 6 we find the accreted abundances are not in Solar proportions for the volatiles H, C, and O (HS2253), verifying that these systems are neither accreting from the interstellar medium nor gas giant planets. Rather for example, the distinctively low carbon abundances in both HS2253 and GD40 are such that C/Mg is more than an order of magnitude below the Solar abundance and therefore more similar to Earth (Jura 2006). It is helpful to compare all measurements to a set of values such as the bulk Earth. In the following we will discuss the abundance ratios of Table 6, mainly through Figure 18, in which the observed abundance ratios are normalized to Mg and the bulk Earth, and Figures 19 and 20 in which we compare the major element compositions of inferred parent bodies with the bulk Earth. HS2253 and PG1225 each have the high-z mass of a large asteroid g currently in their convective zones (Table 3) similar to GD362 (Zuckerman et al. 2007) and GD40 (Klein et al. 2010). HS2253 has one of the highest time-averaged mass accretion rates of all polluted white dwarfs, being comparable to GD362 (Jura et al. 2009; Koester 2009) at g s 1, with SDSS0738 apparently being the current champion at g s 1 based on Dufour et al. (2010) and our settling time calculations parallel to Koester (2009). The time-averaged mass accretion rate of rocky material in PG1225 is a few 10 8 g s 1, about a factor of 5-10 less than GD HS Composition If a WD accretes a rocky planetary body, then the relative element abundances will have the distinctive proportions that arise from mineral oxides (Jura et al. 2009; Klein et al. 2010), made up of MgO, Al 2 O 3, SiO 2, CaO, TiO, Cr 2 O 3, etc. with iron coming from either FeO, Fe 2 O 3, or Fe metal if a differentiated core formed. Another possible vehicle for oxygen in a terrestrial-type planetary body is water ice, which depends on how much hydrogen is present. For HS2253, the HIRES measurement of oxygen plus all major elements, some minor ones, and hydrogen, allows us

11 11 to evaluate the balance of oxides and water in this system according to Equations 2 and 3 of Klein et al. (2010). Looking at the oxide balance in an early phase interpretation for HS2253, we find that the number ratio of oxygen that can be accounted for by oxides with the other elements (using half the iron in FeO, half in Fe 2 O 3, plus aluminum at an Earth-like ratio of Al/Mg = 0.086) is O/Mg = 3.8 ± 0.6. Since the observed amount of oxygen is O/Mg = 5.4 ± 1.1, at the limits of the uncertainties the oxygen can be carried in rocky minerals if iron is predominantly in oxide form; water is not needed to balance the oxides. Generally though, due to the ambiguity of the state of iron and since iron is a major element, there is a sensitive trade-off between metallic iron and water in balancing the oxides. In HS2253 the early phase abundances can be attributed solely to mineral oxides, but could involve up to 15% water by mass (limited by the observed hydrogen abundance). In the steady state for HS2253, the number ratio of oxygen expected to be associated with the observed major elements, Mg, Al, Si, Ca, and Fe, as minerals (half the iron in FeO, half in Fe 2 O 3, plus an Earth-like ratio for Al) is O/Mg = 4.7 ± 0.8, and the ratio of of observed oxygen coming from the parent body is O/Mg = 4.5 ± 0.9. In the steady state the ratios agree very well with a rocky mineral oxide composition. Again, since some Fe may have been metallic and there is a range of uncertainty overlap in the O/Mg ratios, there is still sufficient oxygen to allow for some water; the amount of possible water in the parent body under a steady state interpretation is by mass 7%, limited by the observed hydrogen abundance. While HS2253 does not display an infrared excess, it can be seen from Figures 19 and 20 that this heavily polluted WD has atmospheric abundance ratios that are similar to GD40 an analogous WD that does display an infrared excess. Could the system be in a dissipating phase? Not likely, since the oxides and water are well balanced in both the early phase and steady state interpretations; a dissipative interpretation will eventually run into difficulty with oxides. From the lower panel of Figure 18 in an early phase interpretation, the abundances of Ti, Mn, Cr, Fe, and Ni, are in excellent agreement with bulk Earth values (the Ca overabundance and Si underabundance are similar to GD40). As accretion diffusion proceeds, the ratios move up to the steady state values. If dissipation takes over, the inferred abundance ratios would be driven even further up compared to Mg and Si, while O/Mg decreases, and would imply a parent body that is somehow depleted of the major elements Mg, Si, and O, compared to high abundance ratios of Ca, Ti, Mn, Cr, Fe and Ni. We have no physical explanation for how that may occur. Thus, we prefer the simpler explanation that HS2253 most likely accreted material with abundances predominantly similar to bulk Earth and is somewhere at or between an early phase and a steady state of accretion and diffusion PG Composition In addition to significant high-z element pollution, PG1225 has a considerable amount of hydrogen in its convective zone. Where did it come from? For GD 362, Jura et al. (2009) showed that the presence of circumstellar dust, high-z material and copious hydrogen can simultaneously be explained by the accretion of a large ice-rich planet. Figure 1 of Jura et al. (2009) illustrates that the atmospheric pollution by hydrogen in GD362 and GD16 is

12 12 anomalously high compared to other known helium-dominated externally polluted WDs 4 of their evolutionary epoch, 10,500 K < T eff < 11,500 K, i.e. cooling age Gyr. At a similar cooling age, our set of possible models for PG1225 indicate a hydrogen mass which can range from to grams. Referring to Figure 1 of Jura et al. (2009) this is similar, perhaps on the high-end, but not extraordinary, compared to the average amount of accreted hydrogen in the set of helium-dominated WDs with detected hydrogen at this age. Jura & Xu (2010) suggest that the origin of the hydrogen in these more ordinary WDs may be from the accretion of ice-rich planet(esimal)s which survive post-main sequence evolution. Unfortunately, we cannot evaluate a water and oxide balance for PG1225 since the OI λ7775 lines are not visible in its HIRES spectrum. If the detected elements came from a rocky parent body, then the amount of oxygen (by number) associated with mineral oxides is O/Mg 5, which is far below the sensitivity of our observations that yield an upper limit of O/Mg < 60 (Table 6). With just a small portion of the allowed O/Mg number ratio allocated to rocky oxides, the remaining number ratio that could be associated with water is O/Mg < 55, but that only accounts for less than 10% of the observed hydrogen. In order to explain all of the observed hydrogen as being derived from icy parent bodies, either an ongoing steady state of atmospheric settling from a somewhat large parent body (perhaps a Ceres analog) or previous ice-rich accretion event(s), would be implied. PG1225 s spectrum is very rich with high-z lines and the numerous newly detected elements display an interesting abundance pattern. In the top panel of Figure 18, Mg, Cr, Mn, Fe and Ni, have abundance ratios similar to bulk Earth values, but we find that the refractory elements Ca, Sc, Ti and V are all a factor of 2 3 higher abundance than in the Earth. High Ca and Ti abundance ratios have previously been measured in GD362 (Zuckerman et al. 2007) and GD40 (Klein et al. 2010). As shown in the bottom panel of Figure 18, HS2253 also has a somewhat enhanced Ca/Mg ratio. The enhancement of refractory element abundances with condensation temperatures >1400K is significant and uniform in PG1225, and the upper limits derived for Al and Sr are consistent with the trend. One possible explanation for this observation is that the material accreted onto PG derived from a planetary body that formed in a higher temperature environment than did the Earth. The simulations of Bond et al. (2010) suggest that refractory-rich exoplanets should be a common outcome of planet formation in nebular regions sufficiently close to the host star. PG1225 has just a subtle infrared excess appearing only at 7.9µm (Farihi et al. 2010b; Kilic et al. 2008a). Farihi et al.(2010b)raisethepossibilitythatitcouldbeduetoarelatively narrowdustringwithaninnerholecleared by accretion, suggesting that the system is in a dissipative phase. Referring to Figure 18, since all the refractory elements have shorter settling times than Mg, any correction for settling causes an even greater overabundance of refractory species inferred in the parent body. In other words, contrary to the expectation from a dissipative phase, the heavy elements do not appear sunk. However, the settling times in this WD are quite long Myr, and disk 4 See also Tremblay et al. (2011) for 2 more recently uncovered candidates of this unusual type of WD helium-dominated with large quantities of atmospheric hydrogen.

13 13 lifetimes may be shorter; a large range has been suggested yr (Jura et al. 2009, and references therein). Thus it is conceivable that we could be viewing the system at a point toward the end of an accretion episode but before dissipative differential settling has a significant effect on the atmospheric ratios. Whatever the case regarding circumstellar dust, the atmospheric element ratios imply that PG1225 most likely accreted material with abundances predominantly similar to bulk Earth, but enhanced in refractory species, and the system is somewhere in or near an early phase or steady state of accretion and diffusion Helium-Rich Highly-Polluted White Dwarfs Finally we discuss the set of highly polluted helium-dominated white dwarfs that have been comprehensively studied. To date, 5 helium-dominated WDs have 8 elements reported (GD362, GD40, G241-6, HS , PG ) plus 2 more(sdss0738, GD61) have measurements of all major(o, Mg, Si, and Fe) terrestrial planetary constituents. Though our focus in this paper is on compositions in heavily polluted, helium-dominated WDs, we recognize the importance of the study of heavily polluted, hydrogen-dominated WDs, such as NLTT43806 (Zuckerman et al. 2011) and GALEX J (Vennes et al. 2010). With the detection of circumstellar dust (Debes et al. 2011; Melis et al. 2011) and all the major element species (O, Mg, Si, and Fe), GALEX J (GALEX1931) is almost certainly accreting a planetary body; with settling times on the order of weeks, accretion must be happening now. However we note there is currently a troubling inconsistency between the modeled results of different groups (Vennes et al. 2011;Melis et al. 2011), which is dueinpart to alarge, unexplaineddifference fortheassumedeffective temperature ( T eff > 2,000 K) of GALEX1931, and in part to the treatment of atmospheric settling, which for warm (T eff 20,000 K) hydrogen-dominated WDs, is complicated by the lack of a convective envelope. Figures 19 and 20 represent the compositions by mass of major elements in the helium-dominated polluted WDs in which oxygen has been detected, or a meaningful upper limit set (GD362). Lacking an oxygen detection, PG1225 is not on the bar graphs. There are two different bar graphs showing two different interpretations of the parent body abundances, in the early phase and in the steady state. The graphs are related by ratios of settling times as described in Section 5. We are interested in comparing the bulk accreted compositions, and since hydrogen never settles out of the convective zone, we cannot be certain how much of it is associated with the current heavy element pollution. We do not include hydrogen in Figures 19 and 20, which has a negligible effect on the graph appearance for bulk Earth, HS2253, GD40, SDSS0738 and G241-6, being less than 1% of the accreted mass in those sources. On the other hand, in GD362 hydrogen is over 99% of the accreted mass and in GD61 it is over 75%. As discussed in the preceding section, scenarios have been proposed to explain the presence of hydrogen from the accretion of ice-rich planetary bodies. Now we consider the compositional possibilities from the perspective of oxygen. In each of the WD systems with detected oxygen, there is always enough O to account for the delivery of the

14 14 major 5 elements Mg, Si, Ca, and Fe, in the form of a terrestrial-type planetary body. The upper limit of oxygen in GD362 is just enough to have come from Mg, Si, Ca and Fe in oxides, with the possibility of metallic iron. Since there is very little spare oxygen in GD362, in order to account for the hydrogen, high-z elements, infrared excess, and lack of x-ray luminosity in a unified picture, Jura et al. (2009) calculate a scenario in which settling and accretion have been ongoing for a long time, implying the consumption of a large ice-rich planetary body at least the size of Callisto. On the right side of the bar graphs, the abundance ratios of GD61 imply the accretion of more oxygen than can be carried in rock. An O excess can be the result of dissipative gravitational settling, which would make the abundances seem oxygen-rich and iron-poor. Another way to get an O excess is by the accretion of an ice-rich parent body, as considered in detail by Jura & Xu (2010). The especially high O/Fe and Mg/Fe abundance ratios in the atmosphere of GD61 look suspiciously similar to what is expected in a system where accretion has stopped or paused the heavier elements sink faster than the lighter ones. However, along with measurements of all the major element abundances, Farihi et al. (2011) report an infrared excess for GD61. Since the WD has plenty of hydrogen in its atmosphere, an O excess, and detected circumstellar material, Farihi et al. (2011) raise again the possibility of the accretion of an ice-rich planetary body. Due to abundance measurement uncertainties and the uncertainty of the accretion-diffusion phase, the results are yet inconclusive regarding the composition of the parent body accreted onto GD61. G241-6 also appears to be a case with an O excess, even at the limits of the abundance uncertainties. The upper limit on hydrogen in the atmosphere amounts to less than 0.6% of the mass of detected contaminants, and therefore cannot sufficiently account for the excess oxygen by an origin in ice; the parent body must have been relatively dry. The abundance pattern of G241-6 is similar to GD40 but with Mg, Si, Ca and Fe appearing sunk (compare the GD40 early phase bar of Figure 19 with G241-6 s steady state bar of Figure 20). In other words, starting from GD40 abundances, if we could turn accretion off, then eventually gravitational settling will evolve the abundances into ratios similar to what we observe for G Since we now know that G241-6 does not have an infrared excess (Xu & Jura, Spitzer data, in preparation), it is likely that G241-6 is in a dissipative phase, having originally accreted an ice-free planetary body with composition similar to GD40 and the bulk Earth. The abundances in the remaining three: HS2253, GD40, SDSS0738, are all consistent with the major elements being derived from rocky mineral oxides, and a composition similar to bulk Earth being composed more than 85% by mass in O, Mg, Si, and Fe, for either a steady state and/or early phase interpretation. Overall their parent body compositions are depleted of volatiles; in addition to the dearth of carbon in GD40 and HS2253, all of them have distinctively low abundances of hydrogen (compared to high-z elements) and little, if any, excess oxygen again indicating a dearth of ice in the accreted material. 5 Aluminum may also be expected in comparable mass as Ca, but it is more difficult to detect. Within uncertainties, the undetected presence of Al at an Earth-like abundance will not substantially change the oxide balance.

15 15 While GD362 and GD61 are candidates for exoplanetary systems having accreted a substantial amount of water, the other four stars with oxygen measurements must have accreted relatively dry planetary bodies ( 10% H 2 O by mass). Furthermore, in all four of those cases, the pollution is most plausibly dominated by planet(esimal)s that formed with little or no ice, for reasons as follows. GD40 and SDSS0738 both have an infrared excess, and as discussed by Jura (2008), the existence of a dust disk is best understood as resulting from the disruption of a single parent body, rather than multiple parent bodies which, arriving on trajectories of varied inclinations, tend to result in gaseous disks due to dust grain destruction by sputtering in high-velocity impacts. With a minimum accreted mass of g, and a rock density of 3 g cm 3, the radius of a single polluting parent body is greater than 100 km. Thermal erosion calculations by Jura & Xu (2010) demonstrate that asteroids of radii > 100 km would have retained at least 50% of their internal ice, if present, through red giant evolution. Thus the parent bodies polluting GD40 and SDSS0738 most likely formed without much ice. Could HS2253 and G241-6 be polluted by many small (radius < km), initially ice-rich asteroids whose internal water was lost during the host star s Asymptotic Giant Branch (AGB) phase? Consider a model for HS2253 involving 25 km asteroids that were 25% water by mass before the star was on the AGB. According to Figure 4 of Jura & Xu (2010), unless they originally orbited at less than 5 AU, 25 km asteroids would retain at least 10% of their water during AGB evolution of a 3 M star. At a density of 2.1 g cm 3 (analogous to Ceres), the pre-agb asteroid might have had g of high-z material, and g of hydrogen (H) from ice. After the AGB, the asteroid has g of high-z material, and g of H. With the middle model from Table 4, the convective zone of HS2253 has g of high-z matter. Therefore, it needs to have accreted about 200 asteroids during the most recent 10 5 years, or about 1 asteroid every 500 years. During this time, the white dwarf would acquire g of H, about a quarter of the total H ( g with the middle model from Table 4), which the WD has accumulated over its entire cooling age of approximately 400 Myr. The relatively small amount of H in the convective zone implies that no more than 600 similar asteroids could have arrived in the previous 400 Myr. Why would the asteroids arrive at a rate of yr 1 recently and yr 1 previously? The result of this model exercise is analogous, but more extreme, for G Another difficulty with the above hypothesis is that, even if there is a burstof impacts, one imagines that there would beasize distribution of impactors and, at least in the usual kind of model, most of the mass would be carried in the larger objects. All that is required to produce the observed hydrogen is the accretion of single ice-rich asteroid of radius 100 km. Unless the size distribution of asteroids is very steep, there could easily be one big asteroid for every 200 small ones. Therefore, to explain the pollutions of HS2253 and G241-6 by swarms of initially ice-rich small asteroids, would require both a spike in time of number frequency and and a spike in the size distribution of the asteroids. This double constraint seems implausible. It is most likely that the planet(esimal)s which pollute these four now comprehensively studied WDs HS2253, GD40, SDSS0738 and G241-6 were formed with little ice, and we are measuring the distilled elemental compositions of bodies which condensed out of protoplanetary disks at radii interior to any snow line.

16 16 7. CONCLUSIONS We have obtained and analyzed high-resolution optical spectra of HS and PG with HIRES at Keck Observatory. This work adds two more systems to the small but fast-growing sample of polluted white dwarfs with detailed and comprehensive analyses. Four out of six objects with oxygen measurements most likely accreted planet(esimal)s of mineral oxide constitution, which were nearly ice-free, consistent with formation interior to a snow line. It is clear that extrasolar planetary systems produce rocky bodies that are compositionally similar to terrestrial planets in our own solar system; Earth-like planet(esimal)s apparently do form elsewhere in the galaxy. 8. ACKNOWLEDGEMENTS The authors would like to thank S. Xu and C. Melis for their contributions, and an anonymous referee for suggestions that improved the resulting manuscript. This work has been supported by grants from the NSF and NASA to UCLA. We acknowledge extensive use of VALD, the Vienna Atomic Line Database (Piskunov et al. 1995; Ryabchikova et al. 1997; Kupka et al. 1999, 2000). This publication has made use of NASAs Astrophysics Data System. We thank the Keck Observatory staff for their support. The data presented herein were 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. We recognize and acknowledge the very significant cultural role and reverence that the summit of Mauna Kea has always had within the indigenous Hawaiian community. We are most fortunate to have the opportunity to conduct observations from this mountain. REFERENCES Allègre, C., Manhès, G., & Lewin, E. 2001, Earth and Planetary Science Letters, 185, 49 Becklin, E. E., Farihi, J., Jura, M., Song, I., Weinberger, A. J., & Zuckerman, B. 2005, ApJ, 632, L119 Bond, J. C., O Brien, D. P., & Lauretta, D. S. 2010, ApJ, 715, 1050 Bonsor, A., Mustill, A., & Wyatt, M. 2011, ArXiv e-prints Debes, J. H., Hoard, D. W., Kilic, M., Wachter, S., Leisawitz, D. T., Cohen, M., Kirkpatrick, J. D., & Griffith, R. L. 2011, ApJ, 729, 4 Debes, J. H. & Sigurdsson, S. 2002, ApJ, 572, 556 Dufour, P., Bergeron, P., Liebert, J., Harris, H. C., Knapp, G. R., Anderson, S. F., Hall, P. B., Strauss, M. A., Collinge, M. J., & Edwards, M. C. 2007, ApJ, 663, 1291

17 17 Dufour, P., Kilic, M., Fontaine, G., Bergeron, P., Lachapelle, F., Kleinman, S. J., & Leggett, S. K. 2010, ApJ, 719, 803 Farihi, J., Barstow, M. A., Redfield, S., Dufour, P., & Hambly, N. C. 2010a, MNRAS, 404, 2123 Farihi, J., Brinkworth, C. S., Gänsicke, B. T., Marsh, T. R., Girven, J., Hoard, D. W., Klein, B., & Koester, D. 2011, ApJ, 728, L8+ Farihi, J., Jura, M., Lee, J., & Zuckerman, B. 2010b, ApJ, 714, 1386 Farihi, J., Jura, M., & Zuckerman, B. 2009, ApJ, 694, 805 Friedrich, S., Koester, D., Heber, U., Jeffery, C. S., & Reimers, D. 1999, A&A, 350, 865 Gianninas, A., Dufour, P., & Bergeron, P. 2004, ApJ, 617, L57 Jura, M. 2003, ApJ, 584, L , ApJ, 653, , AJ, 135, 1785 Jura, M., Muno, M. P., Farihi, J., & Zuckerman, B. 2009, ApJ, 699, 1473 Jura, M. & Xu, S. 2010, AJ, 140, 1129 Kawka, A. & Vennes, S. 2006, ApJ, 643, 402 Kilic, M., Farihi, J., Nitta, A., & Leggett, S. K. 2008a, AJ, 136, 111 Kilic, M., Thorstensen, J. R., & Koester, D. 2008b, ApJ, 689, L45 Kilic, M., von Hippel, T., Leggett, S. K., & Winget, D. E. 2005, ApJ, 632, L115 Klein, B., Jura, M., Koester, D., Zuckerman, B., & Melis, C. 2010, ApJ, 709, 950 Koester, D. 2009, A&A, 498, , Mem. Soc. Astron. Italiana, 81, 921 Koester, D., Girven, J., Gänsicke, B. T., & Dufour, P. 2011, ArXiv e-prints Koester, D., Napiwotzki, R., Voss, B., Homeier, D., & Reimers, D. 2005a, A&A, 439, 317 Koester, D., Rollenhagen, K., Napiwotzki, R., Voss, B., Christlieb, N., Homeier, D., & Reimers, D. 2005b, A&A, 432, 1025

18 18 Koester, D., Wegner, G., & Kilkenny, D. 1990, ApJ, 350, 329 Kupka, F., Piskunov, N., Ryabchikova, T. A., Stempels, H. C., & Weiss, W. W. 1999, A&AS, 138, 119 Kupka, F. G., Ryabchikova, T. A., Piskunov, N. E., Stempels, H. C., & Weiss, W. W. 2000, Baltic Astronomy, 9, 590 Limoges, M. & Bergeron, P. 2010, ApJ, 714, 1037 Lodders, K. 2003, ApJ, 591, 1220 Melis, C., Farihi, J., Dufour, P., Zuckerman, B., Burgasser, A. J., Bergeron, P., Bochanski, J., & Simcoe, R. 2011, ApJ, 732, 90 Musielak, Z. E., Winget, D. E., & Montgomery, M. H. 2005, ApJ, 630, 506 Piskunov, N. E., Kupka, F., Ryabchikova, T. A., Weiss, W. W., & Jeffery, C. S. 1995, A&AS, 112, 525 Ryabchikova, T. A., Piskunov, N. E., Kupka, F., & Weiss, W. W. 1997, Baltic Astronomy, 6, 244 Suzuki, N., Tytler, D., Kirkman, D., O Meara, J. M., & Lubin, D. 2003, PASP, 115, 1050 Tremblay, P., Bergeron, P., & Gianninas, A. 2011, ApJ, 730, 128 Tremblay, P., Bergeron, P., Kalirai, J. S., & Gianninas, A. 2010, ApJ, 712, 1345 Tremblay, P.-E. & Bergeron, P. 2009, ApJ, 696, 1755 van Altena, W. F., Lee, J. T., & Hoffleit, E. D. 1995, The general catalogue of trigonometric [stellar] paralaxes (New Haven, CT: Yale University Observatory) Vennes, S., Kawka, A., & Németh, P. 2010, MNRAS, 404, L , MNRAS, 300 Vogt, S. S., Allen, S. L., Bigelow, B. C., Bresee, L., Brown, B., Cantrall, T., Conrad, A., Couture, M., Delaney, C., Epps, H. W., Hilyard, D., Hilyard, D. F., Horn, E., Jern, N., Kanto, D., Keane, M. J., Kibrick, R. I., Lewis, J. W., Osborne, J., Pardeilhan, G. H., Pfister, T., Ricketts, T., Robinson, L. B., Stover, R. J., Tucker, D., Ward, J., & Wei, M. Z. 1994, in Society of Photo-Optical Instrumentation Engineers(SPIE) Conference Series, Vol. 2198, Society of Photo-Optical Instrumentation Engineers (SPIE) Conference Series, ed. D. L. Crawford & E. R. Craine, Voss, B., Koester, D., Napiwotzki, R., Christlieb, N., & Reimers, D. 2007, A&A, 470, 1079 Wolff, B., Koester, D., & Liebert, J. 2002, A&A, 385, 995

19 19 Zuckerman, B., Koester, D., Dufour, P., Melis, C., Klein, B., & Jura, M. 2011, ArXiv e-prints Zuckerman, B., Koester, D., Melis, C., Hansen, B. M., & Jura, M. 2007, ApJ, 671, 872 Zuckerman, B., Melis, C., Klein, B., Koester, D., & Jura, M. 2010, ApJ, 722, 725 This preprint was prepared with the AAS L A TEX macros v5.2.

20 20 Table 1: Keck/HIRES Observations WD UT Date Approx. Seeing Collimator λ range Resolution Exposure λ (visual) (Å) λ (seconds) G Nov Blue Aug Blue x Nov Red HS Aug Blue Nov Red x Nov Red x 1800 PG Feb Blue x Feb 14 not avail Blue Feb Red On 2008 Feb 14 the sky was partly cloudy with an estimated visual extinction of 0.5 mag. The rest of the observations occurred during mostly clear skies. Gaps in blue coverage with our setups occur between Å, and Å. There are many small gaps in red coverage. Table 2:: Absorption Line Measurements of HIRES Spectra Ion λ a χ log gf Equivalent Width (må) (Å) (ev) PG HS G241-6 c H I ,500 (2,000) 850 (120) < 250 H I ,500 (1,500) 120 (30) O I < (90) 250 (50) O I (50) 170 (50) O I (110) 90 (40) O I (190) 210 (50) Na I < 50 < 120 Mg I in 3800 interval 50 (6) 15 (5) Mg I , 0.12 in 3800 interval 220 (40) 75 (14) Mg I , in 3800 interval 470 (40) 106 (14) Mg I (15) 70 (40) 31 (15) Mg I (45) 170 (40) 52 (11) Mg II ,-0.56, (20) 1040 (110) 470 (50) Mg II (110) 110 (25) Continued on Next Page...

21 21 Table 2 Continued Ion λ a χ log gf Equivalent Width (må) (Å) (ev) PG HS G241-6 c Mg II in MgII in MgII Mg II (130) b MgII 250 (40) b MgII Al I < 40 < 60 < 30 Si I < 35 Si II (20) 74 (19) Si II (13) 38 (8) Si II (25) 27 (4) Si II (17) 35 (8) Si II (27) 53 (15) Si II (70) 113 (22) Si II (20) 43 (16) Ca I (60) Ca II (520) 890 (70) 200 (20) Ca II (720) b CaII+CrII 1180 (150) b CaII+CrII 250 (30) Ca II in CaII in CaII (18) Ca II (220) 240 (40) 45 (5) Ca II (340) b FeI 590 (50) 113 (9) Ca II ,000 d 4700 (500) 1300 (150) Ca II in CaII (360) 790 (80) Ca II (400) 570 (80) Ca II (1500) 2340 (320) 780 (90) Ca II (1300) 1720 (160) 290 (90) Sc II (5) Sc II (8) < 80 Sc II in FeI Ti II (8) Ti II in FeII Ti II (4) Ti II (8) Ti II (11) Ti II (30) b FeII 25 (10) Ti II (27) 37 (18) Ti II (15) Continued on Next Page...

22 22 Table 2 Continued Ion λ a χ log gf Equivalent Width (må) (Å) (ev) PG HS G241-6 c Ti II (24) 40 (20) Ti II (9) Ti II (26) 43 (13) Ti II (5) Ti II in 3226 interval Ti II (50) 36 (11) Ti II in 3226 interval Ti II (4) Ti II (120) 140 (30) 34 (7) Ti II (100) 87 (20) 34 (6) Ti II (90) 74 (22) 21 (7) Ti II (15) Ti II (60) 79 (12) Ti II (27) 65 (35) Ti II (18) Ti II (30) b TiII 52 (10) b TiII Ti II in TiII in TiII Ti II (19) 40 (16) Ti II (40) b TiII 65 (29) b TiII Ti II in TiII in TiII Ti II in TiII Ti II (20) b TiII+VII Ti II (9) Ti II (12) Ti II (7) Ti II (15) Ti II (4) Ti II (8) Ti II (50) 72 (15) b FeII Ti II (4) Ti II (30) 38 (10) Ti II (8) Ti II (26) 48 (11) Continued on Next Page...

23 23 Table 2 Continued Ion λ a χ log gf Equivalent Width (må) (Å) (ev) PG HS G241-6 c Ti II (24) 23 (10) Ti II (70) 87 (14) 18 (3) Ti II in TiII (14) 13 (3) Ti II (200) b TiII 180 (20) 60 (5) Ti II (150) 148 (27) 58 (5) Ti II (110) 109 (12) 32 (4) Ti II (24) 42 (13) Ti II (90) 107 (21) 31 (3) Ti II (25) 38 (13) Ti II (25) 27 (8) Ti II (11) Ti II (9) Ti II (6) Ti II (10) b FeI Ti II (16) Ti II (10) Ti II (6) Ti II in TiII in TiII in TiII Ti II (80) b TiII 84 (6) b TiII 20 (2) b TiII Ti II (4) Ti II (100) b FeI 72 (15) 15 (2) Ti II (50) b FeI 67 (23) 24 (3) Ti II (6) Ti II (5) Ti II (7) V II in CrII no data V II (6) < 80 V II in TiII V II (3) Cr II (20) b V II no data 35 (15) Cr II (15) no data 10 (5) Cr II (30) no data 65 (10) Cr II (15) no data Continued on Next Page...

24 24 Table 2 Continued Ion λ a χ log gf Equivalent Width (må) (Å) (ev) PG HS G241-6 c Cr II (10) 43 (15) Cr II in CaII in CaII Cr II (8) 53 (20) Cr II (12) Cr II (11) Cr II (4) Cr II (14) Cr II (3) 36 (13) Cr II (16) Cr II (12) 82 (15) 17 (3) Cr II (6) 42 (6) Cr II (7) 80 (9) 22 (6) Cr II (4) 50 (8) Cr II (6) 72 (8) 23 (7) Cr II (5) 46 (15) Mn II (8) 68 (9) 20 (5) Mn II (9) Mn II in MnII in MnII Mn II (4) b MnII 58 (11) b MnII Fe I (8) b FeI Fe I in FeI Fe I in Ti II Fe I (12) Fe I (25) 45 (15) Fe I (60) 107 (22) Fe I (8) Fe I (20) b NiI Fe I (20) b ScII 47 (19) Fe I (22) 46 (12) Fe I in CaII (15) Fe I in CaII 3737 Fe I (6) Fe I (27) 55 (13) Continued on Next Page...

25 25 Table 2 Continued Ion λ a χ log gf Equivalent Width (må) (Å) (ev) PG HS G241-6 c Fe I in Fe I Fe I in Fe I Fe I (9) Fe I (26) 32 (10) Fe I (10) 46 (15) Fe I in 3827 interval Fe I in 3827 interval Fe I in 3827 interval Fe I (6) Fe I (12) Fe I (8) Fe II (10) b TiII 200 (40) 25 (5) Fe II in Ti II (30) 17 (7) Fe II (16) Fe II (46) Fe II (11) 110 (30) 10 (3) Fe II (6) 80 (12) 14 (5) Fe II (12) b FeII 120 (30) b FeII 25 (10) b FeII Fe II in FeII in FeII in FeII Fe II (15) 103 (15) 16 (3) Fe II (9) 124 (22) 33 (6) Fe II (30) 200 (130) 58 (6) Fe II (50) 410 (60) 90 (10) Fe II (6) Fe II (17) Fe II (5) 75 (53) Fe II (6) 64 (17) Fe II (3) 48 (13) Fe II (5) 68 (24) Fe II (6) 66 (14) Fe II (4) 30 (8) Fe II (10) Fe II in TiII Continued on Next Page...

26 26 Table 2 Continued Ion λ a χ log gf Equivalent Width (må) (Å) (ev) PG HS G241-6 c Fe II (9) Fe II (18) Fe II (7) Fe II (10) Fe II (17) Fe II (13) Fe II (14) 77 (15) Fe II (18) 142 (13) 21 (7) Fe II (45) 210 (20) 42 (4) Fe II (5) Fe II (29) Ni I (12) Ni I (16) Ni I in Fe I Ni II : Sr II < 60 < 90 interval (50) b FeII+TiII see FeII, TiII see FeII interval ,100 (400) b MgI+FeI see MgI, FeI see MgI a Laboratory wavelength in air. Wavelength, χ (excitation energy of lower state) and log gf data from the Vienna Atomic Line Database. b Blended line - contributing element noted in italic with line data given in the table (Table 2). The equivalent width is from the total summed flux of the blend. c G241-6 abundances from Zuckerman et al. (2010) are associated with the EWs of this table. Spectra of G241-6 are shown in Figure 13 of this paper and Figure 4 of Zuckerman et al. (2010). d Combined Ca H&K EW measurement from Koester et al. (2005b) is used since the HIRES data do not fully cover the red wing of this large feature.

27 27 Table 3: White Dwarf Models, Derived Masses and Accretion Rates T eff log g M log M He cvz/m M high Z cvz < Ṁacc > (K) (cgs) (M ) (10 22 g) (g s 1 ) PG1225 model PG1225 model PG1225 model HS2253 model HS2253 model HS2253 model Convective zone (CVZ) helium masses, contaminant (high-z) masses, and time-averaged mass accretion rates (excluding hydrogen) for WD models of varying effective temperature (T eff ) and gravity (g). The (steady state) accretion rates are the sum of individual element rates calculated from the mass of the element in the CVZ divided by its settling time, with abundances and settling times from Tables 4 and 5. For PG the Si upper limit is used, and an abundance of oxygen is assumed as would derive from the dominant mineral oxides: MgO, SiO 2, FeO, CaO. Including carbon at its upper limit (Wolff et al. 2002) increases the high-z masses and accretion rates by 30% for PG , and is negligible for HS Including aluminum at its upper limit increases the high-z masses and accretion rates by 5% for both stars.

28 28 Table 4: HS Atmospheric Abundances by Number and Settling Times for 3 Different T eff /log g Models 14000/ / / / / /8.7 Z log(z/he) log(z/he) log(z/he) log(τ/yrs) log(τ/yrs) log(τ/yrs) H (0.06) (0.06) (0.06) O (0.07) (0.07) (0.07) Mg (0.06) (0.04) (0.03) from Mg I (0.05) (0.06) (0.04) from Mg II (0.05) (0.05) (0.05) Si (0.03) (0.03) (0.04) Ca (0.03) (0.03) (0.05) Ti (0.02) (0.02) (0.02) Cr (0.04) (0.03) (0.03) Mn (0.09) (0.10) (0.09) Fe (0.04) (0.03) (0.02) from Fe I (0.05) (0.04) (0.05) from Fe II (0.03) (0.03) (0.03) Ni : : : C Na < Al < Sc < V < Sr < Element abundances and settling times derived from the equivalent width measurements of Table 2 and different models, as described in Sections 3, 4, and 5. Uncertainties in parentheses are either the standard deviation of the mean for the set of individual line abundances, or the propagated measurement error, whichever is greater. When a significant difference is measured between ionization states of an element (only true for iron), this additional uncertainty is added in. While the absolute abundances and diffusion times vary considerably between models, the relative abundances and relative diffusion times among the listed elements are far less sensitive to changes in the model.

29 29 Table 5: PG Atmospheric Abundances by Number and Settling Times for 3 Different T eff /log g Models 10500/ / / / / /8.3 Z log(z/he) log(z/he) log(z/he) log(τ/yrs) log(τ/yrs) log(τ/yrs) H (0.10) (0.10) (0.10) Mg (0.05) (0.05) (0.05) from Mg I (0.06) (0.06) (0.06) from Mg II (0.10) (0.10) (0.10) Ca (0.03) (0.03) (0.04) from Ca I (0.15) (0.15) (0.15) from Ca II (0.03) (0.03) (0.04) Sc (0.07) (0.07) (0.07) Ti (0.02) (0.02) (0.02) V (0.10) (0.10) (0.10) Cr (0.06) (0.06) (0.06) Mn (0.14) (0.14) (0.14) Fe (0.03) (0.07) (0.08) from Fe I (0.03) (0.03) (0.03) from Fe II (0.03) (0.03) (0.03) Ni (0.14) (0.14) (0.14) C O < < < Na < < < Al < < < Si < < < Sr < < < Similar to Table 4, but for PG

30 30 Table 6: Abundance Ratios By Number a C/Mg from Wolff et al. (2002). b C/Mg from Friedrich et al. (1999). PG HS Solar Convective Zone Convective Zone (Lodders 2003) H/Mg 1700 (500) 3.1 (0.6) C/Mg < 5 a < 0.01 b 6.9 O/Mg < (1.1) 14 Na/Mg < 0.12 < Al/Mg < 0.31 < Si/Mg < (0.10) 0.98 Ca/Mg 0.16 (0.02) 0.13 (0.02) Sc/Mg 9.4± < Ti/Mg 6.7± ± V/Mg 7.1± < Cr/Mg (0.002) (0.003) Mn/Mg 3.1± ± Fe/Mg 0.71 (0.12) 0.89 (0.18) 0.83 Ni/Mg (0.011) 0.06 : Sr/Mg < < Measured atmospheric abundance ratios calculated as the average of 3 models with ratios from Tables 4 and 5. Uncertainties in parentheses are a combination of the propagated abundance ratio measurement uncertainties, and the root-mean-square of ratio values from different models from Tables 4 and 5, added in quadrature.

31 31 #"!!!"+!!"*!,-./ !!")!!"(!!"'!!"&!!"%!!"$!!"#!!"!! %#$!! %$&!! %%(!! %&*!! %(!!! %)$!! %*&!! %+(!! Fig. 1. Combined, echelle-order-merged, shortest wavelength portion (CCD1 with blue collimator) of HIRES spectrum of PG The number density of absorption lines increases dramatically toward UV wavelengths. The absolute flux level and the overall continuum slope are uncalibrated, which is not a problem for the line measurements which are always made relative to their local continuum.

32 32 *+,-./ /.5! *+,-./ /.5! &"!!! B+12!!"%!! A+12! A+122!!"$!! B+12! B+12! 67122! 67122!!"#!! #'#!! ###!! #(#!! #$#!! #)#!! &"!!! B+12! B+12!!"8!! B1!"!"%!! C/122! A+122!!")!! A+122! B+12!!"$!! #%#!! #%8!! #8#!! #88!! (!#!! Fig. 2. Combined, echelle-order-merged HIRES spectrum of HS (black), highly smoothed for clarity in displaying the broad helium features over the large wavelength range plotted. A T eff =14,400 K, log g = 8.4 model, with element abundances from Table 4, is plotted in red.

33 33 '"'! '"!! *+,-./ /.5!!"&!!"%!!"$!?1!"!"#! #(()! #((%! #(#)! #(#%! #($)! #($%! #(%)! 6-0+, :374.;<=>! Fig. 3. Portion of the HIRES data (black) of HS , smoothed by a 5-point boxcar average and displaying Hα. Wavelengths are in air and a heliocentric rest frame. A T eff =14,400 K, log g = 8.4 model, convolved with a Gaussian to approximate the instrumental resolution and with the hydrogen abundance from Table 4, is plotted in red. The echelle order end occurs in the data near 6553Å, thus the continuum blue-ward of the Hα feature is minimally sampled, but is adequate for an EW measurement, which is matched with the model.

34 34 Relative Intensity Cr II He I Cr II Fe II Fe II Fe II Ca II Wavelength (Angstrom) Fig. 4. Portion of the HIRES data (black) of HS , smoothed with a 5 point boxcar average and displaying lines of He I, Ca II, Cr II, and Fe II. Wavelengths are in air and a heliocentric rest frame. A T eff =14,400 K, log g = 8.4 model, convolved with a Gaussian to approximate the instrumental resolution and with element abundances from Table 4, is plotted in red.

35 35 ("(! ("!!,-./ !!"'!!"&!!"%! A-34!!"$! C934! B1344!!"#! )&**! )&)*! )&+*! )&#*! )&$*! Fig. 5. HS Similar to Figure 4, but for Mg I, Si II, and Fe I. This portion of the data was smoothed with a 3 point boxcar average.

36 36 '"'! '"!! )*+,-./*012-*23.-4!!"&!!"%!!"$! >01!!"#! $$((! $$#(! $$$(! $$%(! 5,/*+* :;<=! Fig. 6. HS Similar to Figure 4, but for O I. This portion of the data was smoothed with a 5 point boxcar average.

37 37 &"!!!"#! C=344!,-./ ! D1344!!"!! ''(+! ''#&! ''#+! ''+&! &"&! &"&! &"!!!"%!!"%!!"$! A-344!!")!!"#! B5344!!"$! D1344! E13443F! '(')! '((*! '((+! '#!%! '#&'! '#&$ Fig. 7. HS Similar to Figure 4, but displaying lines of Ti II, Cr II, Mn II, and a marginal detection of Ni II.

38 38 &"&! &"!! &!'!!C$"!! &!'!!C#"#! D,23!!"#! '$(!! '$)!! '$*!! '$$!! '%!!! '%(!! &"!! D,23! F0233! E8233!!"%! ))*(! ))#(! ))$(! ))%( 7.1,-,48/92:;485/<=>?! Fig. 8. Portions of the HIRES spectrum of PG are plotted in black, displaying a helium line at He I λ5876 (top panel, data unsmoothed) and the region near He I λ4472 including the detection of Mg II λ4481 (lower panel, data smoothed by a 3 point boxcar average). Wavelengths are in air and a heliocentric rest frame. The upper panel spans a region across 2 echelle orders causing the gap near λ5896, which is simply filled in with a horizontal line. As described in Section 3.2, it is challenging to simultaneously fit the He I lines and achieve an ionization balance for iron; the red curve is our best compromise. The downward spike at the He line center is not noise, rather our current models do not predict the deep narrow core observed in this star s He I λ5876 line, see also Figure 16 and the last paragraph of Section 2. The sodium upper limit reported in Table 5 is 0.15 dex higher than that of the synthetic Na I D line features near λ5890å shown here. Lower Panel: Instead of a broad shallow helium line, a weak narrow feature (EW 10 må) appears, and is centered at He I λ4472 where no high-z lines are expected. Perhaps this is analogous to the narrow core of the He I 5876 line. A weak narrow line also appears in the data at HeI λ3889. The 10500/8.0 model (green curve) shown for comparison was not used for abundance determinations so the Mg abundance of this model is not matched to the data.

39 39 ("!$! B-34!,-./ !!"'$!!"&$!!"%$! A934!!"#$! )&(!! )&*!! )&)!! )&+!! )&$!! )&#!! )&%! Fig. 9. Unsmoothed portion of the combined HIRES spectrum of PG is plotted in black, displaying a Mg I triplet and many Fe I lines. Wavelengths are in air and a heliocentric rest frame. A T eff =10,800 K log g = 8.0 model, convolved with a Gaussian to approximate the instrumental resolution and with element abundances from Table 5, is plotted in red.

40 40,-./ ! #"$! B=344! #"#! #"!!!"+!!"*!!")!!"(!!"'!!"&!!"%!!"$!!"#! A1344!!"!! %%%)! %%&)! %%')! %%()! %%))! %%*)! %%+)! Fig. 10. PG Similar to Figure 9, but for Ti and Cr.

41 41 +,-./01,234/,450/6! &"!$! 80233! D,233! &"!!!!"%$!!"%!! 7233! 7233! D,233! 80233!!"#$! '()$! '(*!! '(*$! '(#!! '(#$ +,-./01,234/,450/6! &"!!!!"%9! 80233! C023! 80233! 80233!!"##! '$&!! '$&$! '$(!! '$($! '$'!! '$'$! '$9! Fig. 11. PG Similar to Figure 9, but for V and Ni. The data in the lower panel are smoothed by a 3 point boxcar average.

42 42 *+,-./ /.5! &"!!!!"$%! 78122! 6+12! 6+12! 6+12!!"#$! '()!! '()(! '(#!! '(#(! '($!! '($(! '(%! *+,-./ /.5! &"!!!!"%9! D/12! 6+12! 6+12! 78122!!"$$! 6+12! 78122! 6+12! ')!(! ')&!! ')&(! '):!! '):(! ')'!! ')'( Fig. 12. PG Similar to Figure 9, but smoothed by a 3-point (upper panel) and 5-point (lower panel) boxcar average, and displaying Sc II and Fe I lines; a probable blend of a Ni I feature is also labeled. The apparent mis-fit of some Fe I lines is due to a combination of the 0.1 dex discrepancy of Fe I and Fe II abundances, uncertainties in the modeled line data, and an un-modeled shift of Fe I features (see also Figure 14).

43 43 )*+,-./*012-*23.-4! )*+,-./*012-*23.-4! &"!!!"%!!"$! D*011!!"#! '(!%! '(&'! '(&%! '(('! '((%! &"&! &"!!!"6!!"%!!"5! A*01! A*01! C901! B.011!!"$! '%!7! '%&7! '%(7! '%'7! '%#7! '%77! '%$7! '%57! Fig. 13. Unsmoothed portions of the combined HIRES spectrum of G241-6 is plotted in black, displaying Fe II, Mg I, Si II and He I lines. Wavelengths are in air and a heliocentric rest frame. Plotted in red is a T eff =15,300 K log g = 8.0 model, convolved with a Gaussian to approximate the instrumental resolution, and with element abundances from Zuckerman et al. (2010); the model is blue shifted by km s 1, the mean velocity from the full set of measured absorption lines. See also Figure 4 of Zuckerman et al. (2010) for a display of O I lines in G241-6.

44 44 +,-./ :89/5! &%! &$! &#! &"! &!! %! $! #! "!!! &%! &$! &#! &"! &!! %! $! #! "!!! 4<<1<89/5! E/F1<89/5! #!!#&!#"!#'!##!#(!#$!#)!#%!#*!(!!(&!("!('! CD1&""(A!)*! A%! A$! A#! A"!!! "! #! $! %! &!! &"! &# GH1""('I%!"'! Fig. 14. Histograms of (heliocentric) velocity measurements from HIRES data. The velocity bin size is 0.5 km s 1 for PG1225 and 1 km s 1 for HS2253; the larger velocity dispersion for the latter star is likely due to the lower SNR. Fe I lines are a distinguished set in PG1225, whereas the eight Fe I lines in HS span between -3 km s 1 and +12 km s 1.

45 Relative Intensity Wavelength (Angstrom) Relative Intensity Wavelength (Angstrom) Fig. 15. Upper Panel: HIRES observations of Ca II H & K emission cores in the white dwarf PG ; zoomedin display of Figure 1. The local level of the stellar continuum is near 0.55 on this ordinate axis of arbitrary scale. The equivalent widths of the core emission features, relative to the lowest depth of the associated calcium lines, are 50 må centered on the peak at λ Å and 40 må centered on the peak at λ Å. Lower Panel: Core portion of the Hα profile of PG All wavelengths are in air and a heliocentric rest frame.

46 &"!!!!"%$!!"%!!!"#$! ;<3&==$>!(%!!"#!! $#'(! $#(&! $#($! $#(%! $##)! &"!!!!"%$!!"%!!!"#$! <C3&'!!"#!! $#'(! $#(&! $#($! $#(%! $##)! &"!!!!"%$!!"%!!!"#$! <C3)'=!!"#!! $#'(! $#(&! $#($! $#(%! $##)! *+,-.-/012345/061789:! Fig. 16. Narrow cores on the broad He I λ5876 line are detected in HIRES spectra of three helium-dominated, hydrogen-lined, externally polluted WDs in the T eff range of 10,000 K to 12,000 K. PG s observations are from Table 1, GD16 was observed in both February and November of 2008 with no noted variations of the He line, and GD362 s observations are from Zuckerman et al. (2007). The equivalent widths of the core features are: 50 må, 45 må, and 25 må in PG , GD16, and GD362, respectively. For this display each WD spectrum was Doppler shifted to the laboratory rest frame, by its measured redshift, to show the alignment of the core feature with the central wavelength in air λ Å of the helium line.

47 47 "!! &$! &"! %! (! ))*&! +&, )-./)01 +" )2 +& )3 +& 4! #!!! +#! &"#!! &'#!! ""#!! "'#!! &"! &!! %! AB)&"",+!'C! $! #! "!!! &"#!! &'#!! ""#!! "'#!! /:;)*39/2:.<14! Fig. 17. Portions of the spectral energy distributions (SED) of PG (top) and HS (bottom) displaying archival IUE spectra(black, solid line) and GALEX photometry (yellow boxes). The models are convolved with 6Å Gaussians and scaled to match visual photometry. For PG the 10,500 K model is a red, dashed line, and the 11,100 K model is a green, dotted line. For HS the 14,000 K model is a red, dashed line, and the 14,800 K model is a green, dotted line. An extensive list of UV absorption lines is included in the models, with element abundances from the current HIRES measurements (Tables 4 and 5) and carbon from previous UV analyses (Friedrich et al. 1999; Wolff et al. 2002). In both cases the lower temperature model fits the SED well, but if there is some unaccounted-for extinction, then the higher temperature model may be appropriate.

48 48 &" HI#$$'"<2+4G"H=23," HI#$$'"F-,2.G"F-2-," )4" %" F+" )**+,-,.")/01.21*,3"" +,42-56,"-7"89"21.":04;"<2+-=" $" #"!" ##!!" #%!!" #'!!" #(!!" &" %" $" >2"?5" F5" L," F*" J" 81" >+" K5" EF$$'%"<2+4G"H=23," EF$$'%"F-,2.G"F-2-," L," >+" #"!" 81" F5" K5M" >2"?5" )4" ##!!" #%!!" #'!!" #(!! 7 CD" Fig. 18. Parent body abundance ratios from Table 6, with values for Earth from Allègre et al. (2001) and settling times from Tables 4 and 5; 50% condensation temperatures are from Lodders (2003). The fiducial element, Mg, has a 50% T c = 1336 K, nearly the same as that of Fe. Oxygen is off the plot at 50% T c = 180 K. The dotted line at ordinate value 1 is to guide the eye to Earth-like values.

49 49 #!!" <-",-./0"12-34"5-3343"67"5-86."96:;0" +!" *!" )!" (!" '!" &!" %!" $!" #!"!" F?" 5E" D4" B" <-" F?" 5E D4" B" Fig. 19. Accreted parent body mass compositions of major elements not including hydrogen, for an early phase interpretation. For bulk Earth, HS 2253, GD 40, G and SDSS 0738, the omission of hydrogen has a negligible effect on this display, since it is less than 1% of the bulk mass in each of these systems. Conversely, GD 362 and GD 61 have large amounts of hydrogen in their atmospheres, which may or may not be associated with currently observed high-z elements from polluting parent bodies; see Section 6.3. Abundance data for the bulk Earth are from Allègre et al. (2001); GD362 from Zuckerman et al. (2007); GD40 from Klein et al. (2010); SDSS0738 from Dufour et al. (2010); G241-6 from Zuckerman et al. (2010); GD61 from Farihi et al. (2011).