Host-star and exoplanet compositions: a pilot study usinga wide binary with a polluted white dwarf
Amy Bonsor, Paula Jofre, Oliver Shorttle, Laura K Rogers, Siyi Xu, Carl Melis
MMNRAS , 1–8 (2020) Preprint 16 February 2021 Compiled using MNRAS L A TEX style file v3.0
Host-star and exoplanet compositions: a pilot study using a widebinary with a polluted white dwarf
Amy Bonsor (cid:63) Paula Jofr´e , Oliver Shorttle , , Laura K Rogers ,Siyi Xu( 许 偲 艺 ) and Carl Melis Institute of Astronomy, University of Cambridge, Madingley Road, Cambridge, CB3 0HA, UK N´ucleo de Astronom´ıa, Universidad Diego Portales, Ej´ercito 441, Santiago de Chile Department of Earth Sciences, University of Cambridge, Cambridge CB2 3EQ, UK Gemini Observatory/NSF’s NOIRLab, 670 N. A ’ ohoku Place, Hilo, HI 96720, USA Center for Astrophysics and Space Sciences, University of California, San Diego, CA 92093-0424, USA.
Accepted XXX. Received YYY; in original form ZZZ
ABSTRACT
Planets and stars ultimately form out of the collapse of the same cloud of gas. Whilstplanets, and planetary bodies, readily loose volatiles, a common hypothesis is that theyretain the same refractory composition as their host star. This is true within the SolarSystem. The refractory composition of chondritic meteorites, Earth and other rockyplanetary bodies are consistent with solar, within the observational errors. This workaims to investigate whether this hypothesis holds for exoplanetary systems. If true, theinternal structure of observed rocky exoplanets can be better constrained using theirhost star abundances. In this paper, we analyse the abundances of the K-dwarf, G200-40, and compare them to its polluted white dwarf companion, WD 1425+540. Thewhite dwarf has accreted planetary material, most probably a Kuiper belt-like object,from an outer planetary system surviving the star’s evolution to the white dwarf phase.Given that binary pairs are chemically homogeneous, we use the binary companion,G200-40, as a proxy for the composition of the progenitor to WD 1425+540. We showthat the elemental abundances of the companion star and the planetary materialaccreted by WD 1425+540 are consistent with the hypothesis that planet and host-stars have the same true abundances, taking into account the observational errors.
Key words: planets and satellites: composition < Planetary Systems, planets andsatellites: formation < Planetary Systems, planets and satellites: interiors < PlanetarySystems, stars: abundances < Stars, (stars:) planetary systems < Stars, techniques:spectroscopic < Astronomical instrumentation, methods, and techniques stars: late-type
The processes occuring during planet formation have beenkey to our own planet’s history, structure and its develop-ment of life, just as they have been key in determining thenature of planets around other stars. We lack a full under-standing of how planet formation determines the compo-sition of a planet. Whilst planets clearly form out of thesame material as their host stars, a range of processes mayoccur during planet formation and evolution that alter com-positions. Volatiles are easily lost, or are less readily incor-porated, into planetary bodies at high temperatures. How-ever, not all planetary bodies have experienced high tem-peratures. In fact, many comets in the Solar System retain (cid:63)
E-mail: [email protected] clear signatures of their origin in interstellar gas (Mumma& Charnley 2011).The refractory content of planetary bodies will by def-inition not be influenced by the high temperatures expe-rienced, even in the inner regions of planet-forming discs.In our Solar System, many planetary bodies match the re-fractory composition of the Sun. Modelling of the interiorstructure of Pluto and Charon from New Horizons providedno evidence that the interior had a composition differentfrom solar (McKinnon et al. 2017). The refractory compo-sition of chondritic meteorites and the Sun agree to withinerrors on the solar abundances (Anders & Ebihara 1982).Indeed a chondritic reference model has proved a powerfultool for studying the bulk composition of Earth and theterrestrial planets ( e.g.,
Palme & O’Neill 2003), althoughthere are many subtle variations in elemental and isotopic © a r X i v : . [ a s t r o - ph . E P ] F e b ratios which point to the complex nature of planet forma-tion (Boyet & Carlson 2005; Burkhardt et al. 2011; Klaveret al. 2020). For rocky exoplanets, a match in the refrac-tory content of the planet and host star, can significantlyimprove our ability to determine the interior structure ofan observed planet based on its density alone (Dorn et al.2015, 2017b). Given that for most rocky exoplanets, massand radii measurements will remain the best probe of theirinterior structure, this assumption can make a key differencein our ability to characterise detected planets.Like stars and planets, pairs of binary stars should formout of the collapse of the same cloud of interstellar gas. Thus,even if early in their lifetimes, in the centre of dense stel-lar clusters, scattering and exchange of binary pairs is com-mon (Malmberg & Davies 2009), binary pairs should havevery similar compositions, especially true for wide binariesthat have not interacted. Recent works confirm this trend, e.g., Andrews et al. (2019); Hawkins et al. (2020) show thatwide binary stars are more chemically homogeneous thanrandom pairs of stars in the Galactic disk. These findingsare consistent with recent complementary work on chemicalabundances of open clusters (Bovy 2016; Casamiquela et al.2020), showing that member stars are more chemically ho-mogeneous than field stars.Here, we use this principle of chemical homogeneity instars born together to present a novel way to test whetherthe refractory abundances of host stars match their planets.We compare the composition of exoplanetary bodies andtheir wide binary companion. If the refractory compositionof the planetary bodies matches those of their wide binarycompanion, this indicates that the refractory composition ofthe planetary bodies and their host star are likely to be thesame, given the observed chemical homogeneity of wide bi-nary pairs. We present observations and abundance analysisfor the K-dwarf, G200-40, and compare with abundances de-rived for the planetary material in the atmosphere of its widebinary companion, the polluted white dwarf, WD 1425+540.More than a thousand white dwarfs, known as polluted ,have metals from planetary material in their atmospheres(Coutu et al. 2019). Given that white dwarfs should haveclean hydrogen or helium atmospheres, with metals heavierthan helium sinking out of sight on timescales (days to mil-lions of years) much shorter than the white dwarf coolingtime (Koester 2009), the presence of metals in the atmo-spheres of 30-50% of white dwarfs (Zuckerman et al. 2010;Koester et al. 2014) suggests the recent accretion of plane-tary bodies. Outer planetary systems should survive intactto the white dwarf phase and planetary bodies scattered in-wards can be tidally disrupted and accreted by the whitedwarf (Veras 2016). Observations of polluted white dwarfscan tell us about the range of compositions present in ex-oplanetary bodies, as well as the geology of rocky exoplan-ets and processes that alter planetary compositions (Jura &Young 2014). Evidence exists for white dwarfs that have ac-creted icy planetary bodies, rocky planetary bodies and evenplanetary bodies very rich in highly refractory material sim-ilar to calcium aluminium rich inclusions ( e.g.,
Farihi et al.2011; Zuckerman et al. 2011; Harrison et al. 2018).WD 1425+540 is thought to have accreted an icy,volatile-rich body, whose best analog in the Solar System iscomet Halley (Xu et al. 2017). Abundances for the elementsC, N, O, Mg, Si, S, Ca, Fe, and Ni are derived from its spec- tra (Xu et al. 2017). WD 1425+540 has a helium dominatedatmosphere and an effective temperature of 14,490K (Berg-eron et al. 2011). In this work, we present observations of theK-dwarf, G200-40, whose abundances are to be comparedto its wide binary companion, the highly polluted whitedwarf, WD 1425+540. The two stars are separated by ∼ §
2. Thisis followed by a comparison with the abundances derivedpreviously for WD 1425+540 in § § Optical spectroscopic observations were conducted at LickObservatory with the 2.4 m Automated Planet Finder tele-scope (APF; Vogt et al. 2014). APF feeds the Levy Spectro-graph, a stabilized instrument with a minimum of moveableparts that covers most of the optical wavelength range. A1.0 (cid:48)(cid:48) slit was used resulting in resolving powers of (cid:38) ≈
60 ata wavelength of 6000 ˚A. At the end of the night, exposuresof similar K-dwarf, 61 Cyg B were obtained with the inte-gration time limited by the photon integrator (set to obtainsignal-to-noise >
100 in most orders).Data reduction is performed with standard
IRAF echelletasks. Two-dimensional spectral image frames are bias sub-tracted, flat-fielded, extracted, and finally wavelength cali-brated with ThAr arclamp spectra.The spectral analysis was done using the iSpec frame-work (Blanco-Cuaresma et al. 2014), which performs syn-thesis on-the-fly on pre-defined spectral regions until a goodfit to the observed spectrum is found. iSpec is a wrappercode that has the option to call several state-of-the-art ra-diative transfer codes and model atmospheres to performthe syntheses. Here, we employed the code TURBOSPEC-TRUM (Plez 2012) and the MARCS models (Gustafssonet al. 2008), which consider local thermodynamical equilib-rium and one-dimensional atmospheric layers. We consid-ered the atomic data from VALD as well as the line listdeveloped for the Gaia-ESO survey (Heiter et al. 2020).A crucial part of the analysis is to find the set of suit-able regions in the spectrum for the determination of atmo-spheric parameters and then for the abundances ( e.g.,
Jofr´eet al. 2019). These regions need to show dependency on tem-perature, surface gravity, metallicity as well as the elementsmeasured in the white dwarf. Therefore, we used the spectra
MNRAS000
MNRAS000 , 1–8 (2020)
Name WD 1425+540 G200-40 61 Cyg BSource ID 1608497864040134016 1608497623521965184 1872046574983497216 α (deg) 216.89821430610 216.90195396444. 316.75293147546 δ (deg) +53.80838190929 +53.79756047297. +38.75563452356 Table 1.
Gaia
DR2 astrometry for the target K-dwarf, G200-40, its white dwarf companion, WD 1425+540 and the comparison K-dwarf,61 Cyg B. α and δ are barycentric and in ICRS at Gaia
Epoch 2015.5. of the K-dwarf Gaia benchmark star 61 Cyg B (Jofr´e et al.2018), as well as the Sun, to help us to identify the spectralregions for our analysis. These stars are widely used as refer-ence pillars for spectroscopic analysis of automatic pipelinesin stellar surveys. The spectral regions included the wingsof strong lines, in addition to several lines of iron-peak andalpha elements.For G200-40, we determined the following atmosphericparameters. We obtained an effective temperature T eff =4036 ± g = 4 . ± .
48 dex, atotal metal content [M / H] = − . ± .
18, and microtur-bulence parameter v t = 0 . ± .
48 km s − . We further de-termined a stellar mass of 0 . ± . M (cid:12) and radius of0 . ± . R (cid:12) . The uncertainties are internal only, andtheir large values reflect the difficulties in obtaining a goodfit for this spectrum. The cool temperatures imply the spec-trum is crowded with molecules, making it difficult to nor-malise the continuum and find regions free of molecules (see e.g., Lebzelter et al. 2012; Jofr´e et al. 2015). Although weinclude molecules in our line list, the list is still incompletefor such cool stars (Masseron et al. 2014).To validate our results, we used the same procedure toanalyse 61 Cyg B finding parameters of T eff = 4064 ± g = 4 . ± . M/H ] = − . ± . v t = 0 . ± . − , where uncertainties are estimated, but the finalvalues agree within the uncertainties with those reported byJofr´e et al. (2018). We note, however, that our metal contentis lower. This can be explained because the value determinedby Jofr´e et al. (2018) considered a line list without moleculesfor the iron abundance determination. That work reporteda final result which was a combination of several differentmethods. One such methods included molecules, and thatresult yields a metallicity with excellent agreement with ourresult (see Jofr´e et al. 2014, for details about this methodand the procedure in). We, therefore, conclude that althoughour parameters have high uncertainties, our spectral analysisyields accurate results.To derive the abundances, we searched by eye all pos-sible atomic spectral lines in the visual Atlas of Acturus(Hinkle et al. 2000) that are visible in our spectrum. Then,we selected only the lines that had a fit with a conserva-tive uncertainty below 1.5 dex and a root-mean-square ofthe difference between synthesis and observation for a givenline below 0.1. We could not derive abundances of S, C andO from atomic lines, their uncertainties were too large or For solar-scaled abundances, the metal content [M/H] andthe iron content [Fe/H] is the same, but for metal-poor alpha-enhanced stars, [M/H] is slightly higher than [Fe/H]. When de-riving metallicities from alpha-capture elements, it is importantto be aware of this potential difference if the abundances are notsolar-scaled. Element Atoms No. of Molecules No. oflines linesG200-40[C / H] − . ± .
14 1 − . ∗ ± .
07 6[O / H] − . ± .
07 13[Ca / H] − . ± .
04 22[Fe / H] − . ± .
03 101[Mg / H] − . ∗ ± .
14 2 − . ± .
05 14[Ni / H] − . ± .
09 13[Si / H] − . ± .
18 261 Cyg B[C / H] − . ± .
12 1 − . ∗ ± .
05 5[O / H] − . ± .
09 18[Ca / H] − . ± .
03 22[Fe / H] − . ± .
03 101[Mg / H] − . ∗ ± .
24 2 − . ± .
05 15[Ni / H] − . ± .
20 13[Si / H] − . ± .
32 2
Table 2.
Abundances for G200-40 and 61 Cyg B. * indicates theadopted abundances. the lines were too blended. Therefore, we performed synthe-ses for a star with stellar parameters obtained before butvarying O, C, N, Mg and Si to identify possible molecularregions that could serve as alternative. We were still unableto derive abundances for N and S. Our final values are shownin Table 2 for G200-40 and 61 Cyg B, indicating the num-ber of lines or regions for the respective abundances. Thefinal abundances are considered to be the median with thestandard deviation of the distribution as the uncertainty.For C and Mg we can compare our results obtainedfrom atoms and molecules. The value with an asterisk inTable 2 is the adopted one. For carbon, the only atomicline (5052.19 ˚A) was quite blended and yielded an uncertainresult. Furthermore, the resulting value is 0.1 dex higherthan the result obtained with molecules and too high for ourexpectations of 61 Cyg B, given its metallicity. Therefore, weconcluded that the C abundance from that atomic line wasless accurate than the result from molecular features.For the case of Mg, however, we adopted the atomicvalue, even if its uncertainties might seem larger. The Mglines adopted here (5528.43 ˚A and 5711.09 ˚A) are widelyused for Mg determination for galactic stellar populations.On the other hand, very little has been reported in terms ofabundance determination from Mg molecular features. Sincethe values do not agree between the atomic and molecularMg obtained for G200-40, and the atomic Mg value obtainedfor 61 Cyg B agrees well with our reference taken from Jofr´eet al. (2015), we adopted the conservative value, namely theatomic one.
MNRAS , 1–8 (2020) [Fe/H] [ N i / F e ] [Fe/H] [ C a / F e ] [Fe/H] [ M g / F e ] [Fe/H] [ S i / F e ] [Fe/H] [ C / F e ] [Fe/H] [ O / F e ] WD 1425B&FG200-40Cyg B
Figure 1.
The abundances relative to solar determined for G200-40, 61 Cyg B as a function of their metallicity (Table 2), with 1 − σ errorbars. The abundances for WD 1425 +540 are shown as a black line and blue shaded region (1 − σ error bars), as the [Fe/H] cannot becompared. The abundances have not been adjusted for relative sinking of different species, as we assume that the planetary material isaccreting in build-up phase. Stars from Brewer et al. (2016) are included in orange for comparison, noting that typical errors are quotedas being between 0.01 and 0.04 dex. G200-40 has similar abundances to other metal poor stars. The final abundances for 61 Cyg B and G200-40 are pre-sented in Table. 2. As both 61 Cyg B and G200-40 havevery similar stellar parameters, and the same techniqueswere used to derive their abundances, we can compare theirrelative abundances to higher accuracy than to stars withvery different properties (see e.g.,
Jofr´e et al. 2015, 2019;Casamiquela et al. 2020). Fig. 1 shows the abundances rela-tive to solar plotted against their stellar metallicity [Fe/H].Both G200-40 and 61 Cyg B are relatively metal poor([Fe / H] = − .
49 or − .
66 or [M / H] = − . , − .
63) andtheir abundances are fairly consistent with those of othersimilar metal-poor stars in Brewer et al. (2016), shown in or-ange on Fig. 1 for comparison. We note here that the Mg forG200-40 and O for 61 Cyg B do not sit as comfortably withthe Brewer et al. (2016) sample as other elements, poten-tially due to uncertainties in our abundance determination. Jofr´e et al. (2015) discuss in detail the challenges involvedin deriving abundances for the cool dwarfs analysed here incomparison with warmer solar-type stars such as those anal-ysed by Brewer et al. (2016). Issues of note include that thetwo Mg lines detected here were very blended and that ourspectra do not reach the 777 nm range required to determineO abundance from the O triplet and as a result molecularfeatures must be used instead.The abundances for WD 1425+540 are listed in Table 3.We note here that we use the abundances labelled Model Iin Xu et al. (2017), derived using log n(H)/n(He) = -4.2,as recent studies show that the hydrogen abundance derivedusing Balmer lines is more reliable (G¨ansicke et al. 2018;Allard et al. 2020). An additional complication regards thesinking of elemental species out of the white dwarf atmo-sphere. If the accretion started recently (compared to typ-ical sinking timescales on the order of 1 Myr), the systemwill be in build-up phase and the abundances will reflect
MNRAS000
MNRAS000 , 1–8 (2020)
Element log n (Z) /n (He)H − . ± . − . ± . − . ± . − . ± . − . ± . − . ± . − . ± . − . ± . − . ± . − . ± . Table 3.
The log number abundances relative to He derived forthe planetary material accreted by WD 1425+540 taken from Xuet al. (2017), assuming that log n (H) /n (He) = − .
2, as derivedfrom Balmer lines. C a / F e M g / F e S i / F e N i / F e C / F e O / F e Log R e l a t i v e N u m be r A bundan c e s G200-40WD 1425
Figure 2.
The log number abundances relative to iron for WD1425+540 and G200-40, alongside 1 − σ error bars. The whitedwarf abundances have not been adjusted for the relative sinking,as we assume that the planetary material is accreting in build-upphase. All abundances, except C/Fe and Ni/Fe are consistent towithin 1 − σ . those of the accreted body. On the other hand, if accretionhas reached a steady-state between diffusion and accretion,the abundances of the accreted planetary material will havebeen adjusted by its sinking timescale.Harrison et al. (2018) model in detail the compositionof the accreted material and find that the abundances arebest explained if accretion is in the build-up phase. If WD1425+540 were accreting in steady-state, Ca/Fe in the ac-creted planetary body would be significantly higher thanCa/Fe in nearby stars. Ca/Fe can be enhanced in plane-tary bodies due to processing at high temperatures and thedepletion of iron-rich minerals. However, this would be atodds with the high volatile abundance (C, N, O and S) ofthe accreted planetary body. Ca/Fe is also enriched in themantles of differentiated planetary bodies. However, in thiscase, Mg/Fe and other lithophile/sidereophile abundance ra-tios would be enhanced, which is not the case. For the restof this paper, therefore, we assume that WD 1425+540 isaccreting in the build-up phase.In order to compare the abundances observed for theplanetary material in the atmosphere of the white dwarf WD1425+540 with those of its companion, we consider metalsonly (Z > He). The abundances for WD 1425+540 are shownas black horizontal lines, with 1 σ error regions shaded in blue, on Fig. 1. Fig. 2 shows the log number abundances rel-ative to iron for all 6 species observed in both WD 1425+540and G200-40. All number abundances are consistent within1 σ , apart from C/Fe and Ni/Fe, which are consistent within2 σ . Fig. 1 highlights how the comparison is limited by theuncertainties on both the white dwarf and K-dwarf abun-dances. Often the 1 σ uncertainties on the white dwarf abun-dances span many of the typical stellar abundances fromBrewer et al. (2016). Notably the white dwarf abundances sitcentrally compared to the stars in Brewer et al. (2016), apartfrom Ni/Fe, a ratio that can be altered due to planetary pro-cessing, for example segregation into iron melt. Fig. 3 showsthe abundance ratios, alongside their correlated error struc-tures. There is significant overlap in the abundances for WD1425+540 and G200-40, whilst between WD 1425+540 and61 CyB the overlap is more limited, particularly in such ra-tios as C/Fe, O/Fe, C/O and Si/O.We next quantify how likely it is that the true abun-dances of the star G200-40 match those of the material pol-luting its white dwarf companion, given the observationalerrors, and across multiple elemental ratios. For this analysiswe create a set of synthetic abundances, based on the obser-vational errors and an assumption that the true abundancesmatch and ask how often these synthetic observed abun-dances lie closer together than the actual observed abun-dances.In order to make this comparison, we must first trans-form the data via isometric log ratio, such that we take intoaccount all elemental ratios equally, and we consider the fullcovariance error matrix, similar to those shown in Fig. 3.We find that 70% of 10 randomly sampled synthetic abun-dances are closer together (in terms of their Mahalanobisdistance) than the observed abundances. Thus, if the ma-terial accreted by the white dwarf and its companion trulyhad the same elemental abundances, then given the error onthe observations, 70% of the time that they were observed,we would find abundances that are less similar than thosereported here.In other words, this supports G200-40 having the samerefractory element composition as the planetary material ac-creted by its white dwarf companion. The same exercise ap-plied to 61 Cyg B finds that <
25% of the samples are asfar, or further, apart than the real data. Thus, we can rejectthe hypothesis that the abundances observed for 61 Cyg Bmatch WD 1425 +540 at a 75% confidence level.
In this work, we present a pilot study to investigate whetherthe composition of rocky exoplanets matches that of theirhost stars. We present abundance analysis for the K-dwarf,G200-40, and compare with the abundances of the plane-tary material found in the atmosphere of its white dwarfcompanion, WD 1425+540. Assuming that the abundancesof G200-40 are a good proxy for those of the progenitor toWD 1425+540, this technique can be used to compare theabundances of rocky exoplanets and their host-stars.
MNRAS , 1–8 (2020)
Ca/Fe M g / F e Si/Fe N i / F e C/Fe O / F e C/O S i / O WD 1425
61 Cyg B G200-40 Figure 3.
The relative log number abundance ratios for G200-40 (blue), its white dwarf companion, WD 1145+540 (red), and a similarstar, 61 Cyg B (green), listed in Table 2. The pink, blue and green shaded regions indicate the correlated error structure. Seven observedelements leads to six independent ratios to Fe, which are shown on the three left-most panels. In addition to which, the right-hand panelshows, Si/O against C/O.
This work is based on the premise that the abundances mea-sured in the atmospheres of white dwarfs provide good trac-ers of the composition of planetary bodies from an outerplanetary system. Whilst our understanding of the accre-tion of planetary bodies by white dwarfs is good, many de-tails are missing (Farihi 2016). For the purpose of this work,we assume that abundances are not altered during accre-tion, although current models do not yet rule this possibilityout. Whether multiple, or a single planetary body is seen inthe white dwarf atmosphere simultaneously does not changeour conclusions, as the sum of the refractory abundances ofmany bodies should still match that of the host-star. Wenoted earlier in the text, however, that the phase of accre-tion is crucial. We hypothesise that WD 1425+540 startedaccreting planetary material recently (in the last 1 Myr) andis accreting in build-up phase, such that the accreted abun-dances match those of known planetary bodies (Harrisonet al. 2018). However, if the system could be shown to beaccreting in steady-state, the high Ca/Fe and Ni/Fe ratios ofthe accreted material (adjusted for their relative differentialsinking) would be very different to Ca/Fe Ni/Fe observedfor G200-40, or in fact any other stars. Planetary processesthat alter these ratios would then need to be invoked, andit is not clear what these would be.
The analysis presented here points towards the many dif-ficulties in comparing abundances derived for cool main-sequence stars to those derived from atmospheric modellingof white dwarfs. Studies of wide binaries that compare abun-dances between companions can derive very precise differ-ential abundances, particularly if the components of the bi-naries are very similar to each other (Ram´ırez et al. 2019;Tucci-Maia et al. 2019; Hawkins et al. 2020). Here, our starsare so different from each other that we are limited by var-ious large systematic uncertainties, notably on the whitedwarf abundances (Xu et al. 2017) and the K-dwarf abun-dances (see extensive discussions in e.g.,
Jofr´e et al. 2015). The uncertainties on the white dwarf abundances result fromdifferences in the derived abundances depending on the stel-lar parameters (as shown in Table 4 and 5 of (Klein et al.2011)) and abundances that differ depending on the linesused, potentially probing the atmospheric structure of thewhite dwarf (see e.g.,
Table 3 and 4 of (Jura et al. 2012)).We note here that the statistical analysis in the composi-tion comparison was performed assuming random and notsystematic errors.The difference in elemental abundances (Z/Fe) betweenG200-40 and WD 1425+540 is within 1 σ for all elements,except C/Fe. This fits with a model in which volatilesare depleted in planetary bodies, and indicates a strongmatch. However, we note here that the abundances for WD1425+540 would be consistent with many stars selected atrandom. Our analysis finds that the hypothesis that G200-40and its companion WD 1425+540 share the same true abun-dances is more likely than the hypothesis that WD 1425+540shares the true abundances of a random star analysed in asimilar manner in this work, 61 Cyg B. For the former > <
25% for 61 Cyg B. Thisconclusion holds if we exclude O (or if we exclude Mg), wherethe measurement for 61 Cyg B (G200-40) find abundanceswhich are less consistent with other stars of similar metallic-ity. However, considering those elements where we are mostconfident of abundance determinations (Ca, Fe, Si and Nionly), 55% (70%) of the random draws are further apartthan the observed abundances, given the hypothesis thatthe true abundances of G200-40 (61 Cyg B) match those ofWD 1425+540. In other words, a firm conclusion cannot bereached given the large uncertainties on the observed abun-dances and further data are required.
If the refractory abundances of an observed exoplanet areknown to match those of its host-star, the stellar abundancescan be used to improve our knowledge of the planet’s com-position (Dorn et al. 2015), often constraining in particular
MNRAS000
MNRAS000 , 1–8 (2020) the planet’s core size (Dorn et al. 2017a). The analysis pre-sented here provides observational evidence in support ofthis assumption and highlights how stellar abundances caninform our knowledge and understanding of exoplanets.
In this paper, we present a novel means to investigatewhether planetary bodies form with the same composition astheir host-stars. Wide binary pairs are found, in general, tobe chemically homogeneous (Andrews et al. 2019; Hawkinset al. 2020). This means that wide binary companions canbe used as a proxy for the composition of planet-host stars.Analysis of polluted white dwarfs provide the bulk compo-sition of the exoplanetary material that they have accreted.We compare the abundances of the K-dwarf, G200-40, to theplanetary material accreted by its white dwarf companion,WD 1425+540, as a means to compare the composition ofrocky exoplanetary material and their host-stars. Elementalabundances (Ca, Ni, Fe, Mg, Si, C and O) of G200-40 areconsistent with those of the exoplanetary bodies accreted byits white dwarf companion, WD 1425+540, within the ob-servational errors (1 − σ ). In fact, given the observationalerrors, our analysis finds that in 7/10 observations, the de-rived abundances for the two objects would on average befurther apart than in the observations presented here, evenin the case that both the K-dwarf and planetary bodies ac-creted by its white dwarf companion have the same trueabundances. This is in stark contrast to the same analysisfor 61 Cyg B, a similar K-dwarf also analysed in this work,where the null hypothesis that the abundances of 61 Cyg Band WD 1425+540 match can be ruled out at a 75% confi-dence level. We consider this to be evidence in favour of thehypothesis that exoplanetary bodies have the same refrac-tory composition as their host stars, although noting thatour conclusions are limited by the large uncertainties on thedata. Our work supports the idea that host-star abundancescan be used to improve the determination of the interior ofobserved rocky exoplanets. AB acknowledges the support of a Royal Society DorothyHodgkin Fellowship. LR is grateful to STFC and the In-stitute of Astronomy, University of Cambridge for fundingher PhD studentship. We thank Marcelo Tucci Maia forpreparing the spectrum of G200 for analysis, as well as LaiaCasamiquela and Francisca Rojas Espinoza for fruitful dis-cussions. PJ acknowledges partial support of FONDECYTIniciaci´on grant Number 11170174 and ECOS-ANID collab-oration grant Number 180049. SX acknowledges the sup-ported from the international Gemini Observatory, a pro-gram of NSF ’ s NOIRLab, which is managed by the Associ-ation of Universities for Research in Astronomy (AURA) un-der a cooperative agreement with the National Science Foun-dation, on behalf of the Gemini partnership of Argentina,Brazil, Canada, Chile, the Republic of Korea, and the UnitedStates of America. C.M. acknowledges support from theUS National Science Foundation grant SPG-1826583. Theauthors wish to thank Bradford Holden for assistance in scheduling the APF observations presented in this paper.Research at Lick Observatory is partially supported by agenerous gift from Google. The data underlying this article will be shared on reasonablerequest to the corresponding author.
REFERENCES
Allard N. F., Kielkopf J. F., Xu S., Guillon G., Mehnen B.,Linguerri R., Al Mogren M. M., Hochlaf M., Hubeny I., 2020,MNRAS, 494, 868Anders E., Ebihara M., 1982, Geochimica Cosmochimica Acta,46, 2363Andrews J. J., Anguiano B., Chanam´e J., Ag¨ueros M. A., LewisH. M., Hayes C. R., Majewski S. R., 2019, ApJ, 871, 42Bergeron P., Wesemael F., Dufour P., Beauchamp A., Hunter C.,Saffer R. A., Gianninas A., Ruiz M. T., Limoges M.-M., Du-four P., Fontaine G., Liebert J., 2011, ApJ, 737, 28Blanco-Cuaresma S., Soubiran C., Heiter U., Jofr´e P., 2014, A&A,569, A111Bovy J., 2016, ApJ, 817, 49Boyet M., Carlson R. W., 2005, Science, 309, 576Brewer J. M., Fischer D. A., Valenti J. A., Piskunov N., 2016,ApJS, 225, 32Burkhardt C., Kleine T., Oberli F., Pack A., Bourdon B., WielerR., 2011, Earth and Planetary Science Letters, 312, 390Casamiquela L., Tarricq Y., Soubiran C., Blanco-Cuaresma S.,Jofr´e P., Heiter U., Tucci Maia M., 2020, A&A, 635, A8Coutu S., Dufour P., Bergeron P., Blouin S., Loranger E., AllardN. F., Dunlap B. H., 2019, ApJ, 885, 74Dorn C., Hinkel N. R., Venturini J., 2017a, A&A, 597, A38Dorn C., Khan A., Heng K., Connolly J. A. D., Alibert Y., BenzW., Tackley P., 2015, A&A, 577, A83Dorn C., Venturini J., Khan A., Heng K., Alibert Y., Helled R.,Rivoldini A., Benz W., 2017b, A&A, 597, A37Farihi J., 2016, New Astron. Rev., 71, 9Farihi J., Brinkworth C. S., G¨ansicke B. T., Marsh T. R., GirvenJ., Hoard D. W., Klein B., Koester D., 2011, ApJ, 728, L8G¨ansicke B. T., Koester D., Farihi J., Toloza O., 2018, MNRAS,481, 4323Gustafsson B., Edvardsson B., Eriksson K., Jørgensen U. G.,Nordlund ˚A., Plez B., 2008, A&A, 486, 951Harrison J. H. D., Bonsor A., Madhusudhan N., 2018, MNRAS,479, 3814Hawkins K., Lucey M., Ting Y.-S., Ji A., Katzberg D., Thomp-son M., El-Badry K., Teske J., Nelson T., Carrillo A., 2020,MNRAS, 492, 1164Heiter U., Lind K., Bergemann M., Asplund M., Mikolaitis ˇS.,Barklem P. S., Masseron T., de Laverny P., Magrini L., Ed-vardsson B., J¨onsson H., Pickering J. C., Ryde N., Bayo Ar´anA., Bensby T., Casey A. R., Feltzing S., Jofr´e P., Korn A. J.,Pancino E., Damiani F., Lanzafame A., Lardo C., Monaco L.,Morbidelli L., Smiljanic R., Worley C., Zaggia S., Randich S.,Gilmore G. F., 2020, arXiv e-prints, arXiv:2011.02049Hinkle K., Wallace L., Valenti J., Harmer D., 2000, Visible andNear Infrared Atlas of the Arcturus Spectrum 3727-9300 AJofr´e P., Heiter U., Soubiran C., 2019, ARA&A, 57, 571Jofr´e P., Heiter U., Soubiran C., Blanco-Cuaresma S., MasseronT., Nordlander T., Chemin L., Worley C. C., Van Eck S.,Hourihane A., Gilmore G., Adibekyan V., Bergemann M.,Cantat-Gaudin T., Delgado-Mena E., Gonz´alez Hern´andezJ. I., Guiglion G., Lardo C., de Laverny P., Lind K., Ma-grini L., Mikolaitis S., Montes D., Pancino E., Recio-BlancoMNRAS , 1–8 (2020)
A., Sordo R., Sousa S., Tabernero H. M., Vallenari A., 2015,A&A, 582, A81Jofr´e P., Heiter U., Soubiran C., Blanco-Cuaresma S., WorleyC. C., Pancino E., Cantat-Gaudin T., Magrini L., BergemannM., Gonz´alez Hern´andez J. I., Hill V., Lardo C., de LavernyP., Lind K., Masseron T., Montes D., Mucciarelli A., Nord-lander T., Recio Blanco A., Sobeck J., Sordo R., Sousa S. G.,Tabernero H., Vallenari A., Van Eck S., 2014, A&A, 564, A133Jofr´e P., Heiter U., Tucci Maia M., Soubiran C., Worley C. C.,Hawkins K., Blanco-Cuaresma S., Rodrigo C., 2018, ResearchNotes of the American Astronomical Society, 2, 152Jura M., Xu S., Klein B., Koester D., Zuckerman B., 2012, ApJ,750, 69Jura M., Young E. D., 2014, Annual Review of Earth and Plan-etary Sciences, 42, 45Klaver M., Ionov D. A., Takazawa E., Elliott T., 2020, GeochimicaCosmochimica Acta, 268, 405Klein B., Jura M., Koester D., Zuckerman B., 2011, ApJ, 741, 64Koester D., 2009, A&A, 498, 517Koester D., G¨ansicke B. T., Farihi J., 2014, A&A, 566, A34Lebzelter T., Heiter U., Abia C., Eriksson K., Ireland M., NeilsonH., Nowotny W., Maldonado J., Merle T., Peterson R., PlezB., Short C. I., Wahlgren G. M., Worley C., Aringer B., BladhS., de Laverny P., Goswami A., Mora A., Norris R. P., Recio-Blanco A., Scholz M., Th´evenin F., Tsuji T., Kordopatis G.,Montesinos B., Wing R. F., 2012, A&A, 547, A108Malmberg D., Davies M. B., 2009, MNRAS, 394, L26Masseron T., Plez B., Van Eck S., Colin R., Daoutidis I., Gode-froid M., Coheur P. F., Bernath P., Jorissen A., Christlieb N.,2014, A&A, 571, A47McKinnon W. B., Stern S. A., Weaver H. A., Nimmo F., BiersonC. J., Grundy W. M., Cook J. C., Cruikshank D. P., ParkerA. H., Moore J. M., Spencer J. R., Young L. A., Olkin C. B.,Ennico Smith K., New Horizons Geology G., Imaging, Com-position Theme Teams, 2017, Icarus, 287, 2Mumma M. J., Charnley S. B., 2011, ARA&A, 49, 471Palme H., O’Neill H. S. C., 2003, Treatise on Geochemistry, Hol-land H. D., Turekian K. K., eds. ElsevierPlez B., 2012, Turbospectrum: Code for spectral synthesisRam´ırez I., Khanal S., Lichon S. J., Chanam´e J., Endl M., Mel´en-dez J., Lambert D. L., 2019, MNRAS, 490, 2448Tucci-Maia M. T., Mel´endez J., Lorenzo-Oliveira D., Spina L.,Jofr´e P., 2019, A&A, 628, A126Veras D., 2016, Royal Society Open Science, 3, 150571Vogt S. S., Radovan M., Kibrick R., Butler R. P., Alcott B.,Allen S., Arriagada P., Bolte M., Burt J., Cabak J., ChlorosK., Cowley D., Deich W., Dupraw B., Earthman W., EppsH., Faber S., Fischer D., Gates E., Hilyard D., Holden B.,Johnston K., Keiser S., Kanto D., Katsuki M., Laiterman L.,Lanclos K., Laughlin G., Lewis J., Lockwood C., Lynam P.,Marcy G., McLean M., Miller J., Misch T., Peck M., PfisterT., Phillips A., Rivera E., Sand ford D., Saylor M., StoverR., Thompson M., Walp B., Ward J., Wareham J., Wei M.,Wright C., 2014, PASP, 126, 359Wegner G., 1981, AJ, 86, 264Xu S., Zuckerman B., Dufour P., Young E. D., Klein B., Jura M.,2017, ApJ, 836, L7Zuckerman B., Koester D., Dufour P., Melis C., Klein B., JuraM., 2011, ApJ, 739, 101Zuckerman B., Melis C., Klein B., Koester D., Jura M., 2010,ApJ, 722, 725This paper has been typeset from a TEX/L A TEX file prepared bythe author. MNRAS000