Arxiv:2102.01835V1 [Astro-Ph.EP] 3 Feb 2021 Than Helium Sink Rapidly Below Their Surfaces
Total Page:16
File Type:pdf, Size:1020Kb
Draft version February 4, 2021 Typeset using LATEX twocolumn style in AASTeX63 ICY EXOMOONS EVIDENCED BY SPALLOGENIC NUCLIDES IN POLLUTED WHITE DWARFS Alexandra E. Doyle,1 Steven J. Desch,2 and Edward D. Young1 1Earth, Planetary, and Space Sciences University of California, Los Angeles Los Angeles, CA 90095, USA 2School of Earth and Space Exploration Arizona State University Tempe, AZ 85287, USA (Accepted to ApJL, January 6, 2021) ABSTRACT We present evidence that excesses in Be in polluted white dwarfs (WDs) are the result of accretion of icy exomoons that formed in the radiation belts of giant exoplanets. Here we use excess Be in the white dwarf GALEX J2339-0424 as an example. We constrain the parent body abundances of rock-forming elements in GALEX J2339-0424 and show that the over abundance of beryllium in this WD cannot be accounted for by differences in diffusive fluxes through the WD outer envelope nor by chemical fractionations during typical rock-forming processes. We argue instead that the Be was produced by energetic proton irradiation of ice mixed with rock. We demonstrate that the MeV proton fluence required to form the high Be/O ratio in the accreted parent body is consistent with irradiation of ice in the rings of a giant planet within its radiation belt, followed by accretion of the ices to form a moon that is later accreted by the WD. The icy moons of Saturn serve as useful analogs. Our results provide an estimate of spallogenic nuclide excesses in icy moons formed by rings around giant planets in general, including those in the solar system. While excesses in Be have been detected in two polluted WDs to date, including the WD described here, we predict that excesses in the other spallogenic elements Li and B, although more difficult to detect, should also be observed, and that such detections would also indicate pollution by icy exomoons formed in the ring systems of giant planets. Keywords: white dwarfs: exoplanets 1. INTRODUCTION by heavy elements resulting from accretion of asteroid- White dwarfs represent the last stage of stellar evo- like or comet-like bodies (Jura 2003). Many WDs have lution. These stellar remnants are extremely dense and observable debris disks from shredded rocky remnants, have extraordinary gravity such that elements heavier and a few WDs possess evidence of transiting rocky bod- ies (Vanderburg et al. 2015; Manser et al. 2019; Vander- arXiv:2102.01835v1 [astro-ph.EP] 3 Feb 2021 than helium sink rapidly below their surfaces. One would expect to observe only H and He at the sur- bosch et al. 2019). The geochemical compositions of ex- faces of WDs. However, 25 - 50% of WDs exhibit ele- trasolar rocky bodies accreting onto WDs is a burgeon- ments heavier than helium (Zuckerman et al. 2003, 2010; ing field unto itself (e.g., Klein et al. 2010; Zuckerman Koester et al. 2014). These white dwarfs are `polluted' et al. 2010; Jura et al. 2012; Dufour et al. 2012; Vennes et al. 2010; Melis et al. 2011; Farihi et al. 2011; Gaen- sicke et al. 2012; Jura & Young 2014; Xu et al. 2017; Corresponding author: Alexandra E. Doyle Harrison et al. 2018; Hollands et al. 2018; Doyle et al. [email protected] 2019; Swan et al. 2019; Bonsor et al. 2020). Corresponding author: Edward D. Young Recently, exceptionally high and robust over- [email protected] abundances of Be relative to other rock-forming ele- 2 Doyle et al. Table 1. Abundances by number for GALEX J2339-0424 (GALEX 2667197548689621056) from Klein et al. (2021). GALEX J2339-0424 z n(z)/n(He) σspread n(z)/n(Fe) σspread (10−8) (10−8) Be 0.0041 0.0010 3:98 × 10−4 1:40 × 10−4 O 298.9 25.7 29.02 7.74 Mg 26.42 4.0 2.56 0.75 Si 25.6 4.3 2.49 0.75 Ca 0.94 0.28 9:13 × 10−2 3:56 × 10−2 Ti 0.027 0.007 2:62 × 10−3 9:49 × 10−4 V <0.0046 < 4:47 × 10−4 Cr 0.19 0.03 1:85 × 10−2 5:49 × 10−3 Mn 0.094 0.006 9:13 × 10−3 2:38 × 10−3 Fe 10.3 2.6 1.00 0.36 ments (e.g., Fe, Mg, and O) were discovered in two pol- face, and age (e.g., Somayajulu 1977), affording an age luted WDs, GALEX J2339-0424 and GD 378 (Klein et dating technique. On Earth, the formation of primor- al., 2021), and by inference in the planetary materials dial and cosmogenic nuclides by cosmic ray spallation polluting them. Lithium (Li), boron (B) and Be share occurs in the upper atmosphere, and the radio-isotope the characteristic of being the products of spallation re- products, including 10Be, that precipitate onto Earth's actions; all three elements are under-abundant in terms surface are often used to date deep-sea sediments (e.g of cosmic abundances but enriched by cosmic rays (CRs) Arnold 1956; Lal & Peters 1967) and the recycling of sed- as the result of spallation reactions involving collisions iments through volcanoes (Morris et al. 1990). An anal- between protons and carbon and oxygen atoms in the ogous means of using radio-isotope spallation products interstellar medium (ISM). Rather than being produced (collectively referred to as cosmogenic nuclides) has been by stellar nucleosynthesis, Li is destroyed, or astrated, contemplated for dating the surfaces of icy moons in the in stellar interiors at temperatures > 2:5×106 K, and Be solar system (Nordheim et al. 2019; Hedman 2019). and B are destroyed at temperatures > 3:5 × 106 K and In this paper we consider the possible mechanisms for > 5:3 × 106 K(Vangioni-Flam et al. 2000). Lithium is enriching a rocky or icy body in spallogenic nuclides. We easily ionized and thus is more difficult to observe in the consider various sources of MeV protons and evaluate optical regime for stars with higher Teff , such as GALEX the likelihood that these sources could have produced J2339-0424 and GD 378. However, the first detection of the high Be/O observed in bodies accreted by polluted Li in polluted WDs was just recently reported for two WDs. ultra-cool white dwarfs with Teff < 4500K (Kaiser et al. Using Saturn as a model, we find that rings composed 2020). The strongest line of boron can be found in the mainly of water ice within the magnetosphere of a gi- UV (1362.461 A)˚ but sufficient data in the UV have not ant planet satisfy the constraints imposed by the ex- yet been acquired for the polluted WDs discussed here. cess Be concentrations exhibited by the polluted WDs. Nonetheless, the recent detections of Be and Li suggest Mid-sized icy moons of Saturn evidently formed from that a detection of B is highly likely, given sufficient data rings (Charnoz et al. 2009), and indeed accretion of ex- quality. omoons by WDs was anticipated. Moons stripped from Of the isotopes of these elements, only 7Li is produced their host planets were predicted to be a likely source in significant amounts by stellar nucleosynthesis and as of rocky/icy material for pollution of WDs based on an a Big-Bang relic (e.g., Clayton 2003). The stable iso- analysis of post-main-sequence scattering in WD plane- tope of Be is 9Be and it is produced by the reaction tary systems (Payne et al. 2016, 2017). It is worth re- 16O(p,X)9Be where X, the ejected particles, in this case marking that the study of Li and Be as markers of stellar refers to 3pαn. We emphasize again that no 9Be is pro- pollution by planets has a rich history (e.g., Tucci Maia duced in stars. Beryllium is a rare element in the Earth's et al. 2019; Deliyannis et al. 1997). We suggest that crust, as well as in the universe, but concentrations of Be, and perhaps Li and B, excesses in WDs polluted by 10 the rare radionuclide Be (t1=2 = 1:4 Myr) in rocks rocky and icy bodies are signatures of accretion of icy vary with cosmic ray intensities, depth below the sur- exomoons. Icy exomoons evidenced by spallogenic nuclides in polluted white dwarfs 3 We focus our study on GALEX J2339-0424 as an ex- ample of a polluted WD with evidence for excess spallo- o M τz h i genic nuclides, but our general conclusions also apply to PB;z −t/τdisk −t/τz MCV;z (t) = e − e ; (1) GD 378. An evaluation of settling effects on estimates of τdisk − τz the composition of the polluting parent body and a dis- where M o is the initial mass of z in the parent body cussion of the duration of accretion for GALEX J2339- PB;z that forms the circumstellar disk, τ is the e-folding 0424 is detailed in Section2. In Section3 we outline the disk time for the depleting disk mass of parent body mate- proposed scenario for acquiring excess Be in view of the rial, and τ is the e-folding time for diffusive settling of various alternatives. Section4 provides a brief summary z element z through the WD's convective zone. We use of our conclusions. this model to explore the effects of elemental settling through the WD envelope on calculated element ratios 2. THE PARENT BODY ACCRETED BY WD as a function of time. GALEX J2339-0424 Differential settling through the WD envelope may 2.1. The Effects of Settling on Element Ratios cause lighter elements to appear in excess, altering the As reported by Klein et al.