Arxiv:1911.08883V2 [Hep-Ph] 17 Jun 2020 2

White Dwarf Bounds on CHAMPs

1, 1, 2, Michael A. Fedderke, ∗ Peter W. Graham, † and Surjeet Rajendran ‡ 1Stanford Institute for Theoretical Physics, Department of Physics, Stanford University, Stanford, CA 94305, USA 2Department of Physics & Astronomy, Johns Hopkins University, Baltimore, MD 21218, USA (Dated: June 19, 2020) White dwarfs (WD) effectively act as high-gain amplifiers for relatively small energy deposits within their volume via their supernova instability. In this paper, we consider the ways a galactic abundance of (1)-charged massive relics (i.e., CHAMPs) could trigger this instability, thereby destroying oldO WD. The dense central core structure formed inside the WD when heavy CHAMPs sink to its center can trigger a supernova via injection of energy during collapse phases, via direct density-enhanced (pycnonuclear) fusion processes of carbon nuclei dragged into the core by the CHAMPs, or via the formation of a black hole (BH) at the center of the WD. In the latter scenario, Hawking radiation from the BH can ignite the star if the BH forms with a sufficiently small mass; if the BH instead forms at large enough mass, heating of carbon nuclei that accrete onto the BH as it grows in size may be able to achieve the same outcome (with the conservative alternative being simply that the WD is devoured by the BH). The known existence of old WD that have not been destroyed by these mechanisms allows us to improve by many orders of magnitude on the 11 18 existing CHAMP abundance constraints in the regime of large CHAMP mass, mX 10 –10 GeV. Additionally, in certain regions of parameter space, we speculate that this setup∼ could provide a trigger mechanism for the calcium-rich gap transients: a class of anomalous, sub-luminous supernova events that occur far outside of a host galaxy.

CONTENTS 1. Positively charged CHAMPs: X+ 22 2. Negatively charged CHAMPs: X− 23 I. Introduction2 3. Comment on applicability 26 + 4. Combined X and X− case 26 II. Executive summary3 F. Super-Chandrasekhar mass of CHAMPs: BH formation, and timescales 27 III. White dwarf ignition4 A. General discussion5 VII. Black holes in white dwarfs 28 B. Trigger length5 A. Dynamical processes governing BH evolution in a WD 28 IV. CHAMP chemistry8 1. Hawking radiation from the BH 28 2. Accretion of WD material onto the BH 29 V. CHAMPs in the Galaxy9 3. CHAMP accretion onto the BH. Case I: VI. CHAMPs in white dwarfs 11 Ms.g., X < MChand., X . 30 A. Primoridal CHAMP contamination 11 4. CHAMP accretion onto the BH. Case II: 1. Population estimate 11 Ms.g., X & MChand., X . 31 2. Behavior of CHAMPs after protostellar B. Implications for BH evolution 32 cloud contamination 12 C. Implications of BH evolution for old WD 34 B. Accumulated CHAMP contamination 14 1. Evaporating case 34 1. Population estimate 15 2. Accreting case 35

arXiv:1911.08883v2 [hep-ph] 17 Jun 2020 2. Eﬃciency of accretion 17 3. Behavior of accumulated CHAMPs 18 VIII. Limits from the existence of old white dwarfs 36 C. Sinking, stratiﬁcation, and timescales 18 A. Limits on X+ 37 D. Thermally supported B. Limits on X− 39 + CHAMP-contaminated WD structure: C. Limits when both X and X− are present 39 self-gravitating collapse and timescales 20 D. Comment on assumed uniform CHAMP E. Degeneracy supported distribution 41 CHAMP-contaminated WD structure: sub-Chandrasekhar case 22 IX. Speculations: Ca-rich gap transients 41

X. Conclusion 44

∗ [email protected] Acknowledgments 45 † [email protected] ‡ [email protected] A. (NX) binding energies 45 2

B. Pycnonuclear fusion rate estimate 47 accumulated inside a NS will collapse to form a black hole (BH) that can then accrete matter from inside the C. Full expressions for accreting CHAMPs 50 NS, destroying it. Specifically, Ref. [50] considered that halo CHAMPs will pass through and be captured by a D. White dwarf structure 51 protostellar cloud as it collapses to form the M 10– 1. General problem statement 52 30M NS progenitor star. This population of CHAMPs∼ 2. Non-dimensional governing equations 52 remains in the star throughout its lifetime, eventually 3. Chandrasekhar mass 53 contaminating the NS as it is formed. 4. CHAMP-modified WD structure 54 Being extremely heavy, CHAMPs rapidly sink to the center of the NS after its formation,2 and form an ap- E. Electron heat conduction 55 proximately isothermal, thermal-pressure-supported central structure3 that is stable so long as the CHAMPs F. Free-free opacity 56 do not locally come to dominate over the density of NS References 57 material. If however the CHAMPs do come to dominate the NS matter, the CHAMP structure becomes self-gravitating and unstable to collapse (the so-called I. INTRODUCTION ‘gravothermal catastrophe’; see, e.g., Ref. [54]). For a small enough CHAMP population, quantum degeneracy pressure would eventually re-stabilize the structure at a Charged Massive Particles (CHAMPs),1 defined here much smaller radius; however, if a super-Chandrasekhar as massive early-universe relics with (1) electrical X O mass of CHAMPs (MChand. M since mX GeV) is charge, appear in many theories of physics beyond the present, re-stabilization becomes impossible [55], and the Standard Model; for instance, the (N)LSP in R-parity collapse proceeds to BH formation on a short timescale. conserving supersymmetric extensions of the SM can be NS If the BH thus formed has a mass Mbh & Mbh, crit., electrically charged and (meta)stable (e.g., Refs. [1,2]); 20 ∼ 4 10− M , Ref. [50] found that accretion of mate- theories of universal extra dimensions may have as their × lightest KK-odd state a charged state (e.g., Ref. [3,4]); rial from the NS onto the BH, even at Eddington-limited and exotic stable composite bound states [5–8] may have rates, exceeds mass loss from the BH due to the Hawking an (1) net electrical charge without their stability be- process [56], and the BH rapidly grows in size, consuming ingO dramatically impacted. Early interest in superheavy the entirety of the NS on an extremely short timescale. charged particles (e.g., Ref. [9]) long predated their pro- The existence of observed old NSs can thus strongly con- strain the abundance of CHAMPs. On the other hand, posal as an early dark matter (DM) candidate [10, 11] NS (since realized to be ruled out), and they have an ex- if the BH has an initial mass Mbh . Mbh, crit.,, the mass tremely rich phenomenology (see, e.g., Refs. [12–15] for loss by the BH due to Hawking radiation dominates the some recent reviews). mass gain from accretion of NS matter, causing the BH Across large regions of their mass range, a number of to shrink in size and eventually evaporate. This evapora- strong observational constraints limit the abundance of tion however inflicts no externally observable structural CHAMPs to at most a small fraction of the dark matter damage to the NS. As a result, Ref. [50] found that NS density: for instance, limits have been considered from destruction bounds on the CHAMP abundance weaken considerably as the CHAMP mass increases above mX & the absence of anomalously heavy nuclei in bulk terres- 11 trial samples (e.g., Refs. [16–33], and see Ref. [13] for 10 GeV. further references), from their impact on Big Bang Nucle- Recently, there has been much interest in the literature osynthesis [BBN] (e.g., Refs. [9, 34–48]), from the absence in white dwarfs (WD) as unconventional, large space- of accelerated ‘CHAMP cosmic rays’ (see, e.g., Ref. [49] time volume, bolometric particle detectors [57–61]. It for recent work), and from a variety of astrophysical observations (e.g., Refs. [12–14, 50, 51]). For at least some ranges of CHAMP masses, terrestrial bounds have been 2 It is of course conservative to assume a uniform contamination of subject to question owing to the complicated dynamics of the CHAMPs exists throughout the NS precursor stellar mate- charged particles in galactic magnetic fields [11, 49, 52]. rial just prior to NS formation [50]; sinking of sufficiently heavy Some of the more interesting and powerful constraints CHAMPs during the stellar evolution will almost certainly occur. 3 Although the properties of NS material near its center are not on the abundance of CHAMPs, X, at very large masses known with any degree of certainty, it is not entirely comprised arise from the consideration of their impact on the sur- of neutrons: at a minimum, a sizable population of protons vival of old neutron stars (NS) [50]. The crucial idea ad- and electrons co-exists in equilibrium with the neutrons [53]. vanced in Ref. [50] is that a sufficient mass of CHAMPs Although central NS temperatures and densities are much too high to allow the formation of electrostatic bound states of the CHAMP with any particle, the existence of these superposed charged plasmas permits screening of the CHAMP charges, allowing the formation of the thermal-pressure-supported CHAMP 1 Also known more recently as charged Stable Massive Particles structure despite the mutually repulsive electrostatic forces ex- (SMPs). isting between CHAMPs. 3 has of course long been known that white dwarfs—“the cretion, the BH will shrink, but here the outcome is dra- biggest powder keg[s] in the Universe” [62]—can be pro- matically different from the BH-in-a-NS case of Ref. [50]. voked into a thermal runaway leading to a type-Ia–like As the BH evaporates, the Hawking temperatures rises (SNIa) supernova explosion (visible at cosmological dis- and the energy loss rate increases. Provided that certain tances) by the concentrated local deposition of a suffi- minimal assumptions are met, the energy deposition rate cient large amount of energy inside the WD volume [63]. from the Hawking process will eventually raise a suffi- For a carbon–oxygen (CO) WD [64] in the mass range ciently large volume of WD material above the critical 24 Mwd 0.8–1.35M , an energy deposition of (5 10 )– temperature to trigger thermal runaway for WD in this (8 10∼16) GeV is sufficient to raise a region near× the cen- mass range [60, 61]. The authors of Ref. [61] verify that × ter of the WD with a physical size (‘trigger length’) of this latter process completes within . Gyr of BH forma- 2 5 λT (6 10− )–(7 10− ) cm to a critical temperature tion in most regions of parameter space, with only some of around∼ × 0.5 MeV. Such× conditions result [57, 63] in the minor exceptions possible. birth of a stable propagating flame front that traverses the entire WD, burning approximately a solar mass of CO mixture to nuclear statistical equilibrium (NSE), and II. EXECUTIVE SUMMARY releasing 1054 GeV of energy. ∼ WDs thus behave as extremely high-gain natural am- In this paper, we consider the scenario in which plifiers for sufficiently large local energy depositions oc- CHAMPs have accumulated in a WD at first via a pri- curring anywhere within a large fraction of their volume. mordial contamination of the protostellar gas cloud of the This observation has been used to place limits on the WD progenitor star àla Ref. [50], and later in their evo- abundance of primordial black holes which transit thor- lution via direct, gravitationally boosted accretion from ough WD, inducing local heating by dynamical friction a population of halo CHAMPs over the first Gyr of the [57] (but see Ref. [65] for a recent reappraisal); on (asym- WD lifetime. These CHAMPs will settle to the∼ center of metric) dark matter which is captured in a WD, forms a the WD, initially forming a thermal-pressure-supported core, and deposits energy in the WD material once the structure. core becomes self-gravitating and collapses [58, 60]; on If a sufficient total mass of CHAMPs, is present in dark matter (either transiting or captured) which trans- this structure, it will become self-gravitating and undergo fers sufficient energy locally to the WD by scattering, de- gravothermal collapse. While it is possible that the en- cay, or annihilation [59]; and on dark matter that, once ergy released by this collapse alone could trigger thermal captured in a WD, forms a core that collapses to a black runaway, it is conservative to assume that the WD sur- hole, which eventually triggers a supernova [60, 61]. It vives this initial collapse phase. is this latter work on black-hole-induced SNIa-like super- Under that assumption, if a sub-Chandrasekhar mass novae that is of particular relevance to the present work. of CHAMPs participates in the gravothermal collapse, The results of Refs. [60, 61] indicate that, with some the CHAMPs merely re-establish a hydrostatic equilib- minor exceptions, the formation of a BH near the center rium configuration at much smaller radius, supported by of a sufficiently massive WD will either satisfy the heating the degeneracy pressure of relativistic electrons in the requirements to trigger thermal runaway, or will simply WD core, communicated to the CHAMPs by electrostatic consume the WD. The former outcome can occur in one forces. of two ways, depending on which of accretion of mat- Later addition of further CHAMPs to this core, ei- ter or Hawking radiation dominates the BH evolution. ther because the entire primordial CHAMP abundance If the accretion dominates, the BH formed will grow in in the star did not form part of the thermal structure by mass and size. Conservatively, this process would merely the time the gravothermal collapse was initiated, or be- consume the entire WD. However, an second possibility cause further CHAMPs are accreted onto the star from exists [61]: the increasing mass of the BH of course leads the halo, will result in its central density increasing over to ever higher surface gravity, which will eventually lead time, until one of two possible instabilities (see below) to the sufficient acceleration of the carbon ions in the occurs once the core mass reaches the trans- or super- WD material in the immediate vicinity of the BH (or, Chandrasekhar mass regime. more specifically, in the vicinity of the capture radius for Similarly, if a trans- or super-Chandrasekhar mass of accretion onto the BH) to enable local heating the WD CHAMPs participates in the initial gravothermal col- above the critical threshold temperature, leading to ther- lapse, that collapse may initially be (momentarily) stalled mal runway. Alternatively, the increase in the local WD by electron degeneracy pressure causing the reestablish- material density around the BH could allow pycnonuclear ment of a quasi-static core structure, but as the collapsing fusion to proceed even absent significant heating [61, 66].4 CHAMP cloud continues to add more CHAMPs to this On the other hand, if Hawking losses dominate the ac-

to collapse, although the rate estimates for this process are suf- 4 Note that this might also occur in the dense CHAMP-dominated ficiently uncertain that it is not clear that this can occur. See BH-progenitor object that forms at the core of the WD just prior discussion in Sec. VI E 2 b. 4 core structure, its central density will also increase over viewing the ‘chemistry’ of CHAMPs generally in Sec.IV, time, again until one of the two aforementioned instabil- before turning to a discussion of their behavior in galaxies ities occurs. in Sec.V. Following this, in Sec.VI we examine in de- For negatively charged CHAMPs, the first possible tail the accumulation of CHAMPs in white dwarfs, and instability is that carbon ions drawn into the core their behavior in such dense objects generally, with a along with the CHAMPs may begin to fuse immediately particular focus on the dense central core structures they via density-enhanced pycnonuclear fusion mechanisms of form in the WD, as well as the maximum supportable the carbon ions once trans-Chandrasekhar masses are masses allowed before such cores undergo gravothermal reached; this would would directly trigger thermal run- collapse, or collapse to a BH. We give estimates for rele- away. vant timescales in this discussion. This discussion of the Alternatively, if this does not occur, or if the CHAMPs central core structure also includes an examination of a are positively charged, then the growing degenerate core possible earlier trigger mechanism of the SNIa runaway simply reaches and then eventually exceeds the maximum of the WD by the compression of carbon nuclei in the density that is supportable by degeneracy pressure as its core of the WD contaminated by X−. Assuming BH for- mass reaches the trans-Chandrasekhar regime, leading mation, and that WD survives until BH formation, in inevitable to collapse to a black hole. If this second al- Sec.VII we review the subsequent dynamical evolution ternative occurs, the BH-induced WD destruction mech- of the BH showing that, even in this case, the WD is still anisms proposed in Refs. [60, 61] become operational destroyed in large regions of parameter space. We turn and can eventually trigger thermal runaway, unless the these considerations into limits on the galactic abundance evolutionary timescales involved for this process are too of CHAMPs in Sec.VIII by virtue of the known existence long. Except for the caveat about timescales, sufficiently of old, high mass, local WD with low surface magnetic CHAMP-contaminated, sufficiently old WDs are thus de- field. More speculatively, in Sec.IX, we offer some ob- stroyed. servations that could allow this WD ignition mechanism Using the known existence of very old, high mass WD to explain the calcium-rich gap transients. Finally, we of low surface magnetic field [67], we extend and improve conclude in Sec.X. There are a number of appendices upon the upper bounds placed on the CHAMP popula- that offer extra information: AppendixA is a discussion tion in our Galaxy set by Ref. [50], finding that much of the binding energies of X− with positively charged stronger constraints can be set at large CHAMP masses nuclei; AppendixB contains an estimate of the pycnonu- 11 18 clear fusion rate of carbons bound in (CX) bound states (mX 10 –10 GeV); see Figs.4 and5 for our major results.∼ in the central dense CHAMP-contaminated core of the WD; AppendixC gives the full expressions we use to Additionally, making use of this supernova trigger compute the CHAMP abundance accreting onto a WD mechanism at lower CHAMP masses (m 107– X over its lifetime; AppendixD contains a detailed consid- 1011 GeV; see Fig.6), where the considerations advanced∼ eration both of canonical WD structure and CHAMP- in, e.g., Ref. [52] indicate that the majority of CHAMPs contaminated WD structure; AppendicesE andF con- may possibly have been evacuated the central regions tain the expressions we have used for the electron heat of galaxies and magnetically inhibited from re-entry, we conduction and free-free opacity, respectively, that are advance a speculative explanation for a class of sub- used in computing the WD trigger conditions. luminous supernova events that occur preferentially displaced from the center of their host systems, the so-called calcium-rich gap transients (CaRGT). CHAMPs, preferentially present in the outskirts of galaxies after having III. WHITE DWARF IGNITION been ejected from the baryon-rich central regions by supernova shockwaves, can accrete onto WD found in such In order to understand the physical conditions required outer regions, thereby supplying the correct spatial mor- to initiate WD destruction, we begin our discussion by phology for these events, as supernovae of this type are reviewing in detail the arguments advanced in a series of then preferentially triggered in the outskirts of galaxies. papers [57–61, 63, 65] which have examined the condi- Moreover, the CHAMP-induced thermal runaway trigger tions required to initiate thermal runaway in a CO WD, mechanisms we discuss would be able to trigger those su- and which most recently culminated in the proposal of pernovae in the fairly abundant sub-Chandrasekhar pop- Refs. [60, 61] for a trigger mechanism for SNIa-like super- ulation [68, 69] of WD with masses Mwd 0.8M ; recent nova events in WD due to the formation of a BH inside ∼ modeling [70] shows that such supernovae events may the WD (or, more conservatively, for destruction of the give roughly the correct luminosity and spectral charac- WD by accretion onto the BH). Sec.IIIA contains a gen- teristics for the CaGRT. eral discussion of ignition in a WD; we discuss the trigger The remainder of this paper is structured as follows: criteria in detail in Sec.IIIB. in Sec.III we review generally the conditions for WD While we independently re-estimate here various nu- ignition by local energy deposition in the WD material merical values for important physical properties, the bulk in order to establish the physical conditions we will be of this section is based on the discussions in Refs. [57– seeking to attain in this work. We then change tack, re- 61, 63, 65]. 5

A. General discussion to finally violently expand as a result of the significant heating [66]. This catastrophic energy release, being also In ordinary stars, the pressure support necessary to in excess of the total gravitational binding energy of the star, few 1050 erg, leads to the total disruption of the maintain the local hydrostatic equilibrium conditions re- ∼ × quired for a stable stellar configuration is supplied by WD, manifesting itself as a SNIa supernova [66], visible 7 thermal motion (see generally Ref. [66]). An abrupt local at cosmological distances (e.g., Ref. [73]). energy deposition occurring in the stellar material as a result of, e.g., a high-energy particle interaction within the stellar volume will result in a local temperature increase, B. Trigger length causing an increase in the highly temperature-dependent rate of local nuclear burning, leading to further energy The energy deposition necessary for the birth of a sta- deposition and heating. However, this process is self- ble flame front in WD material is discussed in detail in regulating: the temperature increase causes an initial lo- Refs. [57–61, 65], following from the detailed numerical cal over-pressure condition in the heated volume, which work of Ref. [63]. relaxes by causing the heated stellar material to expand As already stated, the basic criterion is that a volume quasi-adiabatically to a larger volume, lowering its tem- of material of a certain minimum characteristic size λT , perature, and with it the rate of nuclear burning in the the trigger length, should be heated to at least a minimum perturbed material. Ordinary stars are thus highly sta- temperature Tcrit, the trigger temperature, such that the ble against having their structure significantly disrupted rate of energy loss by thermal diffusion from this volume by small local energy depositions. to the surrounding medium, E˙ diff.(Tcrit, λT ), is smaller In a WD however, the hydrostatic pressure support than the rate of energy deposition− into that volume by required for stability is supplied almost entirely by the nuclear fusion, E˙ nucl.(Tcrit)[63]. If this criterion is not quantum degeneracy pressure of the (possibly quite rela- satisfied, any localized initial temperature perturbation tivistic) electrons [64], communicated to the heavier ion- impressed into the material by an external source (e.g., ized nuclei by electrostatic forces. This pressure is sen- a high-energy particle interaction occurring in the WD) sitive to the electron density, but is highly insensitive to may cause local nuclear burning to proceed all the way 5 local temperature perturbations [71]. Local energy de- to NSE locally, but will ultimately simply diffuse away positions in WD material resulting in a local heating of harmlessly into the bulk WD medium without triggering the ions are thus not subject to the self-regulating adia- the formation of a self-sustaining, propagating nuclear batic expansion mechanism as ordinary stellar material. flame front. This makes WD highly susceptible to a thermal run- The criterion may also be alternatively phrased as re- away condition [62, 64]: some local energy deposit leads quiring the deposition of a certain total minimum energy to local heating, which leads to a (highly) increased local ET , the trigger energy, into the trigger volume within nuclear burning and further local heating, and so on. The the characteristic time for heat diffusion from the trig- only temperature regulation mechanism available is ther- ger volume τdiff.; energy deposition into a volume larger mal diffusion, either by conduction by degenerate elec- than the trigger volume must be proportionally larger on trons, or by radiative transport. If the rate of thermal a volumetric basis. energy flow by diffusion out of some (sufficiently large) We will be mostly concerned with cases where the tem- perturbed volume is lower than the rate of energy in- perature of the trigger volume, which we assume to be lo- jection from nuclear burning in that volume, the local cated near the center of a CO WD lying in the mass range temperature in the perturbed volume simply continues Mwd 0.8–1.35M , is raised at least as high as Tcrit. to rise as the nuclear burning rate increases with increas- 0.5 MeV∼ [63]. The corresponding densities of the WD ma-∼ ing temperature. Moreover, the region of increased nu- 6 terial are obtained by solving the Tolman–Oppenheimer– clear burning will propagate outward as a flame front, Volkhoff (TOV) equation [75, 76] to find the WD-mass– which eventually traverses and consumes the whole star central-density relationship, ρ (M ), which for WDs [62, 66]. The total energy released by the nuclear burn- C wd with masses Mwd [0.1, 1.35]M can be given approxi- ing of an (1) fraction of the material in the WD to ∈ NSE, 10O51 erg, easily exceeds the energy required to lift the∼ electron degeneracy, allowing the stellar material

7 The ‘canonical’ trigger scenario for a SNIa is by mass accretion from a non-degenerate binary companion star, resulting in the WD mass increasing toward the Chandrasekhar limit [55], be- 5 So long as means ∆T E where E is the Fermi energy of coming ever more unstable to perturbation in the process, with | | F F the degenerate electrons; typically, EF (few MeV). A tem- ignition occurring shortly before the limit is reached owing to perature increase ∆T MeV would lift∼ degeneracy O and return a combination of core heating and compression as a result of the stellar material to a normal gas phase. the mass accretion [66, 72, 74]. Alternative scenarios such as 6 Whether this occurs subsonically (a deﬂagration), supersonically double-degenerate collisions [72], in which two degenerate, sub- (a detonation), or some combination (deﬂagration-to-detonation Chandrasekhar WD collide and trigger the necessary runaway, transition) remains an open problem [72]. have also been proposed. 6

mately by8 E˙ (T , λ ) = +4πkλ2 ( ∂ u). We approximate − diff. crit T T − r ∂ru ∆T/λT for a region of radius heated λT heated 6 3 2 3/2 − ∼ ρ 1.95 10 g/cm α(M )− 1 (1) to a temperature Tcrit = T0 +∆T above the ambient tem- C ∼ × wd − 2 3 perature T0, so that E˙ diff.(Tcrit, λT ) 4πkλT ∆T . For α(Mwd) 1.0033 0.3087x 1.1652x + 2.0211x − ∼ ≈ − − a region so heated under isobaric conditions, the excess 2.0604x4 + 1.1687x5 0.2810x6 energy is − − x M /M [0.1, 1.35]; (2) wd Z T0+∆T ≡ ∈ 3 ∆E = M dT 0cp(T 0) = (4πλT /3)ρc¯p(Tcrit)∆T, the central densities of WDs in the mass range Mwd T0 0.8 1.35M thus vary from 107 g/cm3 to 109 g/cm∼3. (4) To− evaluate the trigger criterion,∼ we must∼ know the rate of diffusive heat transport. Diffusive heat flow where cp(T ) is the constant-pressure specific (per-mass) in presence of a temperature gradient is governed by heat capacity, and we have implicitly defined its average Fourier’s Law: Q = k∇u, where Q is the (vector) heat value in the relevant temperature range. The diffusive flow, u = u(t, x) is− the temperature field, and k is the cooling timescale is thus thermal conductivity [78]. The heat flow mechanisms in τdiff. ∆E/ E˙ diff. (5) dense WD matter are dominantly radiative transport and ≡ |" | # electron conduction [57, 63, 66],9 so the thermal conduc- ρλ2 1 Z T0+∆T = T dT c (T ) (6) tivity is given by [66] 0 p 0 3k ∆T T0 2 3 2 4π T 1 1 ρλT c¯p(Tcrit)/(3k). (7) k = krad + kcd = κrad− + κcd− , (3) ∼ 45 ρ The isobaric specific heat capacity can be approximated by as a sum of three terms, independently account- where κrad, cd are, respectively, the radiative and conductive opacities; these in general depend on the density, ing for the ionic, electronic, and radiative contributions. temperature, and chemical composition of the material. We treat the ions as a free ideal gas, and keep only the Electron thermal conductivities are given in Ref. [79] and leading term in the Sommerfeld expansion [71] for the electronic contribution, valid in the limit T E . That do not have a simple scaling at all densities, although for F 8 3 is, we take ρ & 8 10 g/cm and at T 0.5 MeV, the scaling is × 1.4 ∼ approximately κcd ρ− . The radiative opacity can be ions elec. rad. cp(T, ρ) c (T, ρ) + c (T, ρ) + c (T, ρ) (8) taken to be given by∝ the free-free Kramers opacity which ≡ p p p 5 X Xi scales as κrad ρ, corrected by a suitable Gaunt fac- cions(T, ρ) (9) ∝ p ≡ 2µ A tor [66, 80–82] (although we ignore this correction here); a i i see AppendicesE andF for details. For T 0.5 MeV, " # 1 ∼ 2 2 − the electron conduction contribution to the thermal con- elec. π T me 8 3 cp (T, ρ) 1 (10) ductivity dominates at densities above ρ 8 10 g/cm , ≡ µaµe EF − EF ∼ × while the radiative contribution dominates at lower den- 3/2 9 3 4 3 " 2#− sities; outside the range ρ 0.3–3 10 /cm it is a good rad. 4π T me ∼ × cp (T, ρ) 1 (11) approximation that k scales with density as implied by ≡ 5µaµe EF − EF the dominant of the two individual contributions. The total rate of thermal energy loss from a spheri- where µa is the atomic mass unit, EF = [1 + 10 2 3 2/3 1/2 cal region of radius λT by heat flow through its surface (3π ne/me) ] is the electron Fermi energy, ne = ˙ 2 R is Ediff.(Tcrit, λT ) = λT dΩ rˆ Q(Tcrit). Assuming a ρ/(µaµe) is the electron number density, µe − · 1 ≡ radial temperature profile such that u = ∂rurˆ we have P ∇ ( i XiZi/Ai)− is the mean molecular mass per electron, and Xi, Zi, and Ai are, respectively, the mass fraction, charge, and atomic mass number of ion species i. It follows that 8 We obtained this result assuming a fully degenerate electron equation of state (EoS) with electrons assumed to have their 1 Z T0+∆T c¯ (T ) dT 0c (T 0) (12) full dispersion relation (i.e., we did not use the non- or ultra- p crit ≡ ∆T p relativistic approximations) and exact local charge neutrality of T0 the WD plasma; we did not however include Coulomb corrections 1 cions(T , ρ) + celec.(T , ρ) to the EoS (see, e.g., Ref. [77]). We then fitted a polynomial func- ≈ p crit 2 p crit tion to the numerical parameter α which completely dictates the 1 WD structure. See AppendixD. rad. + cp (Tcrit, ρ), (13) 9 Losses by neutrino emission are subdominant [63, 65], except 4 during the final stages of thermal runaway [66]. assuming T ∆T T . 10 This argument applies parametrically to any convex shape, with crit ∼ 0 the radius of the sphere replaced by the smallest characteristic Assuming a specific (again, per-mass) nuclear energy scale over with the heated region has a significant temperature generation rate of S˙nucl.(Tcrit), the total nuclear en- gradient. ergy generation rate in the same volume of material 7

˙ ˙ 3 ˙ is Enucl.(Tcrit) = MSnucl. = (4π/3)λT ρSnucl., so that the passage of a flame front necessarily involves rapid timescale for nuclear burning to double the excess ther- changes in both temperature and density, Eq. (19) does mal energy in a region of mass M is not necessarily give a very good numerical approximation to the actual trigger length [71, 85]. There is more- τnucl. ∆E/E˙ nucl. (14) over an ambiguity as to the exact temperature at which ≡ the nuclear energy generation rate should be evaluated 1 Z T0+∆T = dT c (T ) (15) [63, 71, 85]: although we have indicated this temperature ˙ 0 p 0 Snucl. T0 as Tcrit., we remind the reader that screened carbon fu- c¯p(Tcrit.)∆T/S˙nucl.(Tcrit.). (16) sion proceeds as a tunneling process for T . MeV, and ∼ the specific nuclear energy generation rates scale as fast as d ln S˙ /d ln T 24 in the temperature and density A parametric, order-of-magnitude estimate for the nucl. range relevant for WD∼ flame propagation [83, 84], making trigger condition can then be phrased as even (1) errors in temperature highly relevant. TheO correct procedure to determine the trigger length E˙ (T ) E˙ (T , λ ) τ = τ (17) nucl. crit ∼ − diff. crit T ⇔ nucl. diff. and trigger mass is to perform a numerical simulation of S˙ (T ) 3k∆T/(ρλ2 ) 3kT /(ρλ2 ), (18) flame propagation in the style of Ref. [63], using a large ⇒ nucl. crit ∼ T ∼ crit T network of nuclear reactions to accurately evaluate S˙nucl.. where the final expression holds in the physically relevant The results of Ref. [63] for the trigger mass cover only a limit T T . Therefore, the trigger length estimate small range of densities ρ [0.2, 10] 109 g/cm3 assuming crit 0 ∈ × is given by XC = XO = 1/2, appropriate to the central densities of WD in the mass range Mwd [1.25, 1.41]M , which does q not cover the full (i.e., lower)∈ WD mass range of inter- λ 3kT /ρS˙ (T ); (19) T ∼ crit nucl. crit est to us. We therefore follow the procedure of Ref. [57] and analytically scale the results of Ref. [63] to densi- see also Refs. [63, 71]. ties lower than those numerically sampled in that refer- However, Eq. (19) should be employed with great care: ence by making use of the parametric scalings implied 11 the specific nuclear energy generation rates in the WD by Eq. (19) [and assuming that S˙nucl.(Tcrit.) ρ for 9 3 2 8 3∝ 12 temperature range are extremely fast functions of tem- ρ . 10 g/cm ]: λT ρ− for ρ 8 10 g/cm [57]. perature, and become also exponentially fast functions of Explicitly, and although∝ the density scaling× that results density not too far above the unperturbed WD central has a discontinuous slope at the boundary of the follow- densities for the WD we consider, owing to strong screening piecewise definition, we scale the results of Ref. [63] ing effects and pycnonuclear processes [66, 83, 84]. Since for T 0.5 MeV as follows: ∼

 4 2 1.3 10− cm (ρ/ρ1)− ρ ρ1  × × ≤  λ0.5 MeV λ (ρ/ρ )ln(λ2/λ1)/ ln(ρ2/ρ1) ρ < ρ < ρ , (20) T ≈ 1 1 1 2   5 2.5 10− cm ρ = ρ × 2

where in the range of numerical values sampled by Ref. [63] we have made a simple log-log linear interpo- 13 11 This assumption is already mildly violated at the upper end lation, owing to the graphical results of Ref. [63] being of this density range (for T 0.5 MeV): owing to strong ∼ difficult to reliably read off in the intermediate regime. screening effects that increase the parametric scaling of the 8 3 ˙ Here, λ1,2 λT (ρ1,2), with ρ1 = 2 10 g/cm and rate with increasing density, we have d ln Snucl./d ln ρ 1.1 by ≡10 3 × 8 3 ˙ ∼9 3 ρ = 1 10 g/cm being the endpoints of the numeri- ρ 10 g/cm , and d ln Snucl./d ln ρ 1.3 for ρ 10 g/cm ; 2 × this∼ is actually indicative of the the onset∼ of the∼ pycnonuclear cally sampled values in Ref. [63]. regime, in which the reaction rates become exponentially sensitive to density [66, 83, 84]. 12 Note that if we had to extrapolate above the numerically sam- 0.8 pled range of Ref. [63], we would find λT ρ− for ρ 8 3 ∝ 0.5 8 10 g/cm , in disagreement with the scaling λT ρ− higher densities is unnecessary given that the numerical results used× in Ref. [57]. This difference arises directly from∝ the dif- of Ref. [63] range up almost to the extremal central density of ference in the scaling we find here for κcd from Ref. [79]; we do the Chandrasekhar-mass WD. 13 not however need to resolve this discrepancy as extrapolation to That is, a linear interpolation of log λT as a function of log ρ. 8

by CHAMPs contaminating a WD, we will now change 1024 tack and consider the behavior of CHAMPs in various relevant settings. We begin that discussion in this sec- ] 1022 tion by recalling come basic properties of CHAMPs.

GeV In this paper we consider CHAMPs to be either singly [

T 20 + 10 positively (X ) or singly negatively (X−) charged; many E of our constraints in principle operate for CHAMPs with 1018 other charges, but such particles lie beyond our explicit consideration here. The relative abundance of the population of the X± 11 10− CHAMP states depends on details of the CHAMP production model, and whether CP is violated in the pro- ] s

[ duction; as we are entirely agnostic to the production

diﬀ. model, in this paper we consider any generic scenario in

τ 12 10− which the species are produced either symmetrically or asymmetrically (but always non-thermally, to avoid over- closure constraints [86]) in the early Universe, and survive until today. Having (1) electric charge, CHAMPs 0.8 0.9 1.0 1.1 1.2 1.3 are of course subject to richO and complicated dynamics Mwd[M ] and ‘chemistry’ after they are produced; see generally Refs. [9, 10, 49, 50, 52]. FIG. 1. Upper Panel: Energy deposition required with a Positively charged CHAMPs X+ have quite straight- trigger volume within a trigger time in order to initiate ther- forward chemistry, as the only stable negatively charged mal runaway, as a function of the WD mass. We assume particles available for them to bind with are electrons; Tcrit. 0.5 MeV and that the BH forms near the center of they thus have atomic-sized cross sections. the WD.∼ Lower Panel: Diffusion timescale for the trigger Negatively charged CHAMPs X on the other hand volume as a function of the WD mass, under the same as- − sumptions. have much more complicated behavior: they are able to form bound states [denoted (NX)] directly with nuclei [generically denoted N]14 [9, 10, 34]. The (NX) for N Armed thus with parameterizations for the ther- heavier than He are in particular are very deeply bound, mal conductivity [Eq. (3)], isobaric specific heat ca- with typical ground-state binding energies of a few MeV, pacity [Eqs. (8)–(13)], trigger length assuming that and average ground-state radii of order the nuclear size or smaller; see Tab.II and AppendixA. These bound Tcrit. = 0.5 MeV (which is likely a conservatively high value) [Eq. (20)], WD-mass–central-density relationship states, being net positively charged in all cases except [Eqs. (1) and (2)], trigger energy [Eq. (4)], and diffusion (pX), are thereafter free to bind with electrons, and in timescale [Eq. (7)], it is now possible to compute the general have atomic-sized interaction cross sections with trigger energy required to be injected to initiate thermal electrically charged Standard Model particles. runaway. We show our results for the trigger energy and The case of (pX) is unique, as the state is neutral, has a 25 keV binding energy, and a ground-state Bohr radius diffusion timescale in Fig.1, for WD in the mass range 2 of rB r 30 fm. As such, this state behaves almost Mwd [0.8, 1.35]M ; the trigger energy results we find ∼ 3 h i ∼ are more∈ conservative than those in Fig. 1 of Ref. [59] by as a heavy neutron, with a highly suppressed electro- up to 1.5 orders of magnitude. magnetic interaction cross section with matter, of order We will be primarily interested in massive WD with its geometrical size [10, 11]. However, in dense environments, the (pX) state is susceptible to disruption via the Mwd 0.85M or Mwd 1.1M . We show representative∼ physical characteristics∼ for such WD (and one energetically favorable and classically allowed ‘exchange’ reactions (pX) + N p + (NX); combined with the other) in Tab.I. For the remainder of this paper, we will → sometimes quote results for critical masses and timescales fact that the (pX) are the last states to recombine as the early Universe cools, this results in only a small frac- both generally, and also specifically for WD of these two 4 tion 10− of the total X− abundance being bound in masses; in the latter cases, we will denote these specific ∼ values as ‘WD 0.85’ and ‘WD 1.1’, respectively. (pX) states primordially; see, e.g., Refs. [42, 46–48]. An (1) fraction of the CHAMPs become bound in (HeX) primordiallyO [42, 46–48].15 IV. CHAMP CHEMISTRY In all cases however, given their extremely large mass,

In the previous section we reviewed the physical conditions required to bring about thermal runaway in a WD 14 We distinguish a generic nucleus N from a nitrogen nucleus N by localized energy deposition in the WD volume. As we by the use of italic and Roman fonts. will be arguing that such conditions can be brought about 15 Note that this is in contrast with claims in the earlier literature 9

TABLE I. Physical characteristics for two representative classes of WD. Quoted are the mass in solar masses, radius in solar 3 radii, central density in g/cm , sound speed cs at center of the WD as a fraction of the speed of light c, trigger energy in GeV, and diﬀusion time for a trigger volume in seconds. Rwd, ρc, and cs are obtained using information from the solution of the TOV equation, assuming a fully degenerate electron EoS, without Coulomb corrections (see AppendixD).

3 Mwd [M ] Rwd [R ] ρc [g/cm ] cs/c ET [GeV] τdiﬀ. [s] 2 7 2 24 11 0.85 1 10− 1.2 10 2 10− 1 10 1.5 10− × 3 × 7 × 2 × 20 × 12 1.1 7 10− 5.5 10 2.8 10− 2 10 3 10− × 3 × 8 × 2 × 18 × 13 1.2 6 10− 1.3 10 3 10− 2 10 3 10− × × × × ×

CHAMPs have highly suppressed cross section-to-mass On the other hand, CHAMPs that have been previ- ratios, so their dynamics are vastly simpler than those of ously processed through a stellar environment and then ordinary baryonic matter; they are more akin to those of re-released to the galaxy either during the late stages dark matter, with some notable exceptions. For instance, of stellar evolution, or in a previous galactic supernova heavy CHAMPs are prevented from collapsing into dif- event, will of course in general be bound up to heav- fuse gas clouds in galaxies [49, 50, 52] and can be present ier nuclei as (NX): because nuclear processes with rates in galactic haloes, distributed much like dark matter (at comparable to reactions such as (HeX)+N (NX)+He 11 → least for mX & 10 GeV). clearly do occur at appreciable rates in stars, so too do Because the vast majority of X− are bound primor- such exchange reactions; and the supernova event itself dially to He, and such CHAMPs (if sufficiently heavy) may involve many highly energetic processes capable of do not become bound to diffuse collapsed structures in overcoming any Coulomb barriers and altering the chem- galaxies, the majority of X− that ever manage to enter a ical state of the CHAMPs. Nevertheless, the fraction of stellar environment (see discussion below) will do so for the total number of CHAMPs in any galaxy so processed the first time as a (HeX) bound state from a halo distri- will be small, for the simple reason that the CHAMPs 11 bution of CHAMPs. For heavy X− (mX & 10 GeV), should be roughly homogeneously mixed with the hydro- we will thus make the assumption that the entire X− gen and helium gas primordially in the early Universe, abundance that contaminates or accretes onto a star will and the total fraction of hydrogen and helium in a galaxy initially enter that star as a (HeX) bound state. that gets processed through stars into heavier elements 5 11 For lower masses, mX 10 GeV–10 GeV, how- is small. In any Population I or Population II star then, ever, baryonic dynamics such∼ as supernova shockwaves the total abundance of CHAMPs present that may be in galaxies can dramatically impact the spatial morphol- bound to heavy nuclei as a result of processes occurring ogy and momentum distribution of CHAMPs (both X+ before the birth of that particular star will be no larger and X−)[49, 52], as we review in Sec.V. Even in this than a fraction approximately at the level of the metallic- case, it is highly improbable that any processes occurring ity of such a star (i.e., roughly at the percent level), and outside of a stellar environment could alter the ‘chemi- therefore negligible. Of course, CHAMPs that do enter cal’ nature of any X− that exist in the form (HeX): pro- a star are almost all subsequently processed in that star cesses than can disturb the X− chemical nature, e.g., itself over its lifetime, but the above considerations imply (HeX) + N (NX) + He will necessarily involve a that, as the starting point for considering the evolution Coulomb barrier→ similar in size to that for fusion of a of CHAMPs in a star, it is reasonable to assume that proton onto a heavy nucleus N, which is not a process all CHAMPs initially enter it either as X+ or bound as that happens spontaneously at any appreciable rate out- (HeX), with some exceptions that we discuss when they 16 side dense stellar environments. Even these lighter X− may be important. thus likely remain in the form (HeX) until then enter a dense stellar environment for the first time. For lower masses still, the CHAMPs in the form (HeX) may sim- V. CHAMPS IN THE GALAXY ply collapse into the baryonic disk [49, 50, 52]; however, unless these CHAMPs again enter a dense stellar envi- As we have already mentioned, for sufficiently mas- ronment, they too likely survive as (HeX) for similar sive CHAMPs (m 1011 GeV [49, 52]) the spatial reasons. X & morphology and momentum distribution of CHAMPs in the Galaxy is not significantly impacted by baryonic dynamics; such massive CHAMPs should thus have a halo that a much larger fractional abundance of (pX) exists primor- distribution akin to that of the dark matter. Lighter dially [10, 50]. CHAMPs however can be significantly impacted by such 16 Moreover, the reduced mass of the (HeX)–N system is larger baryonic dynamics [11, 49, 52]; see also discussion in than for p fusion onto N, being µ mN rather than µ mp; this further suppresses the tunneling≈ rate through the Coulomb≈ Refs. [50, 51, 88, 89]. Such considerations are themselves barrier, which is what controls the rate of escape of the X− to subject to significant uncertainty and there is some vari- the heavier nucleus in such an exchange process. ation in the literature as to the fate of such CHAMPs. 10

simplified picture of the dynamics of the diffuse CHAMP TABLE II. Numerically computed ground-state binding en- plasma under the combined action of complex shock- ergies EB [MeV] and average radii [fm] of the (NX) system, where the nucleus N is modeled as a uniformly charged wave dynamics, and the (turbulent) magnetic fields in 1/3 galaxies. More recent results, e.g., Ref. [49], find that al- ball of radius R = A r0 with r0 = 1.22 fm (see, e.g., though (1)-charged CHAMPs in the mass range m Refs. [9, 34, 87]); see AppendixA. Also shown are the na¨ıve O X ∼ estimated hydrogenic Coulomb binding energies that would 105 GeV–1010 GeV are sufficiently shock accelerated such be obtained ignoring the finite charge radius of the nucleus; that an (1) fraction of those CHAMPs in the mass these significantly overestimate the true binding energies for range indicatedO above may be evacuated, their diffusive N heavier than helium, owing to the bound state being local- re-entry into the disk is not completely inhibited by the ized within the nuclear volume. magnetic fields; instead, an equilibrium between accel- na¨ıve NEB [MeV] E [MeV] r [fm] R [fm] eration and diffusion is reached, and the CHAMP mo- B h i p 0.025 0.025 43 — mentum spectrum is significantly altered. Their results 10 4He 0.35 0.40 6.1 1.9 also indicate that for mX & 10 QX GeV, such CHAMPs 8Be 1.6 3.2 2.6 2.4 form a DM-like halo that is not disturbed by baryonic dy- 12 5 2 C 2.9 11 2.1 2.8 namics, and that for mX . 10 QX GeV, such CHAMPs 16O 4.1 25 1.8 3.1 will collapse into the disk of the MW as it forms. 24Mg 6.1 86 1.7 3.5 56Fe 10.0 940 4.1 4.7 Refs. [49] and [52] thus agree that sufficiently massive (1)-charged CHAMPs will form a DM-like halo, but theyO disagree on the details of how the CHAMP abundance, momentum distribution, and spatial morphology Ref. [52] argued that efficient evacuation from, e.g., for lower CHAMP masses is impacted by baryonic dy- the Milky Way (MW) disk, of CHAMPs that are not in 4 namics. Taken together, the conservative conclusion to the form (pX) (which are in any event a . 10− fraction draw from these results is that any CHAMP bounds that of CHAMPs [42, 46–48]) is possible by the combined ac- rely on knowing the galactic abundance or momentum tion of supernova shockwaves, and the confining effects of 5 11 distribution of CHAMPs for mX 10 GeV–10 GeV galactic magnetic fields. Their argument is that a (1)- are subject to significant uncertainty.∼ charged CHAMP that is heavier than m 105Q2 OGeV, X ∼ X where QX is either the X charge, or the charge of the We do however stress that those CHAMPs bound in state in which X finds itself bound, will be accelerated the form (pX) and which have not undergone some re- by the periodic passage of a sufficiently intense supernova action that converts them to (NX)[N = p] would not shockwave through the MW disk, and be unable to radi- get significantly impacted by supernova6 shockwaves or ately dissipate the kinetic energy thus gained sufficiently magnetic field confinement, as they are not charged and quickly to avoid being accelerated above the disk escape hence can be neither efficiently entrained in a SN shock- speed on a timescale of (100 Myr). Moreover, the galac- wave, nor deflected by a magnetic field. Hence, for a tic disk magnetic fields,O which are mostly in the plane of CHAMP population which has any significant X− com- the disk, magnetically confine virialized CHAMPs in and ponent, there is always a residual fraction of X− in the around the disk to within approximately their gyrora- galaxy; this is however estimated to be no larger than dius. Given the MW CHAMP virial speed and typical 4 10− of the X− abundance [42, 46–48]. 1–10 µG magnetic fields, virialized CHAMPs lighter ∼ 11 + than mX 10 QX GeV have gyroradii smaller than the Finally, we note that annihilation of X and X in ∼ − MW disk thickness. Ref. [52] thus argued that CHAMPs diffuse structures is inhibited: the X+ and (HeX) states not in the form (pX) which thus find themselves ini- which are assumed to be the dominant forms of CHAMPs tially in the MW disk would be rapidly accelerated out in diffuse structures at all times in the late Universe are of the disk by supernova shockwaves, and that virialized both positively charged. Therefore, the rate for bringing CHAMPs in the MW halo would be prevented from en- an X+ and a (HeX) to within fm of each other, in or- tering the disk by the confining action of the galactic der to allow the X+ wavefunction∼ to significantly overlap disk magnetic fields, provided their mass is in the range with the highly localized X− wavefunction, is exponen- 5 2 11 10 QX GeV . mX . 10 QX GeV. Similar arguments tially suppressed by a Coulomb barrier. The suppression presumably hold for other disk galaxies, and some cog- is moreover exponentially more severe than the tunnel- nate argument is also likely to hold for the central regions ing suppression factor for pp fusion, since the tunneling of elliptical galaxies. exponential is suppressed by the reduced mass (see, e.g., We note again that since the majority of the CHAMPs Ref. [90]), which is here µ mX /2 mp/2. Combined present in non-collapsed structures in a galaxy will still with intrinsically low CHAMP∼ densities in diffuse struc- be in the form X+ or (HeX), with or without an electron tures, any annihilation process that would tend to remove bound to them, it is appropriate to take QX = +1 for a symmetric component of the CHAMP abundance is + both cases X±, when evaluating the boundaries of the simply much too slow to be relevant, and the X and region in which Ref. [52] claims an effect. X− bound as (HeX) survive independently throughout The results of Ref. [52] are however based on a greatly the era when they are distributed diffusely. 11

VI. CHAMPS IN WHITE DWARFS A. Primoridal CHAMP contamination

We consider ﬁrst the case where CHAMPs accumulate in the protostellar cloud from which the WD progenitor Having discussed the chemistry, history, and behavior forms. We begin with a review of the arguments ad- of CHAMPs outside stellar environments, we can now vanced in Ref. [50] leading to the population estimate, turn to the question of how CHAMPs come to contami- and then turn to the question of how the CHAMPs be- nate WD, and their impact on WD dynamics. Our dis- have as the star evolves. cussion is guided in outline initially by that of Ref. [50].

1. Population estimate In this section, we discuss first how CHAMPs come to be present in WDs and their evolution over time, con- The X, in whatever state they find themselves after sidering in turn two mechanisms: (1) in Sec.VIA, we recombination, will collapse into galaxies during early- consider the same case as Ref. [50], in which CHAMPs ac- Universe structure formation. For both the X+ and cumulate in the protostellar cloud that collapses to form (NX) [i.e., (HeX)] forms of CHAMPs, σsm-X /mX will the stellar progenitor of the WD, leading to a contami- diffuse capture be small enough, if m m , that even the nation of the WD material as the WD is born. And (2), X & X multiple orbits executed through the MW galactic disk in Sec.VIB, we consider the case where halo CHAMPs since galaxy formation would not be efficient in capturing additionally accumulate onto the WD as during the first the CHAMPs into diffuse disk gas clouds [50]. Estimate 1 Gyr of the WD existence (timescale to be discussed diffuse capture diffuse capture 5 below).∼ for mX vary from mX 10 GeV diffuse capture 7 ∼ [49] to mX 10 GeV [50]. For lighter X, such trapping would be∼ efficient, but magnetic heating Having considered the population of CHAMPs that effects could conceivably eject from the disk any such X can be present in a WD, we then estimate the timescale which are trapped in diffuse clouds [50]. It is therefore for that population to sink to the center of the WD unclear if a significant fraction of the X in, e.g., the MW [Sec.VIC], and the evolution of the initial structure end up trapped in parts of the disk containing only dif- diffuse capture formed by the CHAMPs at the center of the WD fuse gas if mX . mX . On the other hand, X [Sec.VID]. We then consider in turn the fate of the can be trapped by denser collapsing protostellar clouds; CHAMPs in the cases of sub- or trans-Chandrasekhar Ref. [50] estimates that a protostellar cloud of molecular

[Sec.VIE], and super-Chandrasekhar [Sec.VIF] total gas of mass Mcloud would capture a total mass MX± of CHAMP masses. For the latter case, this discussion will the respective charge species of CHAMP, X±, giving rise naturally evolve into a discussion of black hole dynamics to star fractionally polluted at birth by each of X± by in a WD, which is the topic of the immediately following an amount η MX± /Mcloud given by (see below for section [Sec.VII]. further discussion± ≡ of the net contamination) [50]

 √ 3/2 1/2 1/4 192 2π α MPl. f ρX Mcloud  ± g(ycrit.) ycrit. ymin.  3 2 3/2 7/4  7e v m me m ≥ η rot. p X , (21) ± ≈ 13/6 2/3 1/3 4 2 3 2 π α MPl.me f ρX 7 4 4 7  ± 1 ymin. + ymin. ycrit. < ymin.  7e 3 19/6 min 7/6 1/3 − 3 3 vrot.mp (nH ) Mcloud where

( Z 1 4 ) 1 7/4 7 4 4 7 x h 4q4 i g(x) ln x− − 1 x + x + 7 dq x− 1 ; (22) ≡ − 3 3 ymin/x q −

ycrit is deﬁned by the equation whose solution is given in terms of the Lambert W - 17 1 1/4 4 function: ycrit. = (4β)− W (4β) ; ymin. dictates the ycrit. exp βycrit. = 1 (23) 2 4 2 β = 4α MPl./(mpmX Mcloud), (24) 17 1 The asymptotic expansions of which are given by x− W (x) 1 x + 3 x2 for x 1 and W (x) ln(x) ln(ln(x)) for x 1.≈ − 2 ∼ − 12

min onset of X-capture when the cloud reaches nH = nH and can thus be ignored. 2 3 ∼ 10 cm− , becoming sufficiently UV-shielded to be molecular: The transition between the low-mX regime (ycrit. < ymin.) where η is independent of mX , and the high-mX 1/6 2/3 min 1/6 1/3 7/4 4π m (n ) M regime (ycrit. ymin.) where η m− occurs because y p H cloud ; (25) ≥ ∝ X min. 1/2 the phase-space available for capture of more massive X ≡ 3 M m α Pl. e is suppressed relative to that for lighter X until the cloud min v 220 km/s is the local circular speed in the MW; grows somewhat denser than nH and is hence physically rot. ∼ ρX is the total CHAMP mass density at the location smaller at the onset of efficient capture [50]. of cloud collapse; and f are the fractions of the total ± Note that since the majority of the CHAMPs present CHAMP abundance in the forms X±, respectively. In in non-collapsed structures in a galaxy will still be in the computing f , we include all X− bound in (NX) states, form X+ or (HeX), it is appropriate to take q = +1 for except for (pX− )[50], for which this specific estimate does X both cases X±, when evaluating η . not apply because it is based on an electromagnetic in- ± teraction between the CHAMP and the diffuse gas which For an Mcloud = 4M protostellar cloud (appropriate is absent for the neutral (pX) state; however, they only for an Mwd 0.8M WD [91]), the numerical estimate 4 ∼ constitute . 10− of the total X− abundance [42, 46–48] for η , 4 η (Mcloud = 4M ) is ± ≡ ±

 2 f ρX 4 2.2 10− ± mX . 1.7 10 GeV  × ρhalo × η , 4 , (26) ± ≈ 7/4  2 f ρX mX − 4 2.2 10− ± gˆ(m ) m 1.7 10 GeV  × ρ 1.7 104 GeV X X & × halo × where ρ 0.3 GeV/cm3 is the local MW DM this assumption, the total abundance of CHAMPs that halo ∼ halo density, andg ˆ(mX ) g[ycrit.(mX )]/g[ycrit.(1.7 are expected to be present in the WD from this primor- 4 ≡ × 10 GeV)]. Given the asymptotic scalings of ycrit., it fol- dial contamination can be estimated as lows that—for this size cloud—ˆg(mX ) supplies only a log prim. 3/4 M η Mwd. (27) 12 X± correction for mX . 10 GeV, but thatg ˆ(mX ) mX ≈ ± 14 ∝ for mX & 10 GeV, which changes the power-law scaling We assume further [50] that the X contamination is ap- of η for large mX . proximately uniform throughout the WD core at the time CHAMPs that accumulate in protostellar clouds will of its formation. inevitably be incorporated fairly uniformly into the stars In the next subsection, we examine further aspects of formed by such clouds, as they are fairly well mixed; as the behavior of the CHAMPs between the time of proto- a result, those stars are born with a baked-in CHAMP stellar capture and WD birth that can impact the validity contamination fraction η as estimated by Eq. (21)[50]. of the approximation Eq. (27). To turn this estimate into an estimate for the total abundance of CHAMPs that are present in the WD at the end of the stellar evolution, we however need to make further 2. Behavior of CHAMPs after protostellar cloud contamination assumptions. Since we assume that mX mp, some gravitational sinking of the population of X+ and (NX) toward the center of the star will undoubtedly occur dur- X+ will merely reside in the star without much change, ing the lifetime of the WD-progenitor star leading up to whereas (NX) [mostly beginning as (HeX)] bound states the formation of a WD-progenitor CO core, which would can undergo significant processing in the core of the potentially boost the fractional mass contamination of star, as the nuclear reaction Q-values for a large num- that core relative to that of the initial protostellar cloud. ber of steps in a large number of different stellar burn- However, large-scale convective processes that can oper- ing reactions exceed the (NX) binding energies for the ate in certain regions, and during certain phases of the nuclei involved (e.g., the modified CNO cycle reaction evolution, of stars that actively burn nuclear fuel make p + (14NX) 15O + X is possible with Q 3.8 MeV), it challenging to give a quantitative estimate of this ef- which would→ enable ejection and subsequent≈ recapture of fect, and we follow Ref. [50] in making the conservative the X− by a different nuclear species. Catalyzed nuclear assumption that the fractional mass pollution of the WD- such as (4HeX) + 14N 16O + X, or exchange reac- progenitor CO core, and hence the WD itself upon forma- tions such as (4HeX) +→14N 14(NX) + 4He are also tion, is just that of the initial protostellar cloud. Under in principle energetically allowed→ (although they proceed 13 as tunneling processes through a Coulomb barrier), and for, e.g., (CX) + O C + (OX) exchange reactions (or could result in the transfer the X out of (4HeX) and exchange reactions→ with heavier nuclei) to occur after into heavier bound states. However, note that even if (CX) formation, such reactions are highly tunneling sup- the X− bound to helium do not exchange onto other pressed: the C wavefunction of the (CX) bound state heavier nuclei, the limits we will ultimately set will still (treating the X is stationary as it is so heavy) is local- conservative; see discussion in Sec. VI E 2. ized within fm of the X position, so the incoming O must also approach∼ to within fm of the X to be cap- As it is most energetically favorable for an X− to bind ∼ with nuclei of higher charge (and hence greater mass) the tured. But since the O nucleus simply sees a Q = +5 natural end result of these processes would be the even- object until it is within fm of the (CX) bound state, there is a Coulomb barrier∼ of few MeV for this to oc- tual accumulation in old, highly evolved stellar cores of ∼ (N 0X) bound states where N 0 is the heaviest nucleus cur. Moreover, since the reduced mass of the (CX) + O present in significant quantity. The approach to the fi- system is about twice that of a C+O system, the tunnel- nal distribution of X in various heavy nuclei is in gen- ing suppression is much more pronounced that even for eral complicated, and detailed study would require ex- C + O fusion (even though the charge is slightly lower, tensive numerical modeling utilizing a large reaction net- the reduced mass change more than compensates; see the work populated with rates for all the reactions including discussion of Gamow energies in, e.g., Ref. [90]), which is the X in addition to the usual nuclear reaction rates, as itself suppressed compared to the C + C fusion probabil- is done in ordinary solar modeling (e.g., Refs. [92–94]); ity. Since the majority of C-ions in a non-extremal WD such an effort is far beyond the intended scope of this have not undergone a fusion reaction since CO core for- work. In particular, in the highly evolved M M WD- mation (else the core would have already burned), this progenitor CO stellar core of a red giant with∼ an initial implies that reactions like (CX) + O C + (OX) are exceedingly unlikely to have occurred. → mass in the appropriate range (i.e., Mstar . 9M [95]), 12 16 the X− will be present in the form ( CX), ( OX), or We do however mention that a CNO-like reaction such possibly bound to some heavier trace species (e.g., Ne, as p+(12CX) 13N+X could still lead to the ionization → Mg, etc.) that was able to form in smaller quantity dur- of the X; however, in highly evolved stellar cores (either ing the stellar evolution; all these (NX) bound states in the He-rich environment where the ersatz triple-α pro- have binding energies of at least 3 MeV. cesses occurs, or already in the degenerate CO core) the p Moreover, even if we were to assume that some signif- abundance has already been almost completely depleted 4 by ordinary stellar burning [66], making this possibility icant population of ( HeX) were to survive unscathed 18 during earlier burning phases in the stellar evolution, rare. Thus, even if such reactions were to occur in once such states become part of the dense He core of the He-rich core prior to the He flash, the X would still an old, massive red giant shortly before the helium flash that leads to the formation of the WD-progenitor core of CO material, the X present in (HeX) can catalyze fusion reactions such as (HeX) + He (BeX) + γ, which 18 Incidentally, we note that catalyzed fusion reactions such as 8 → (pX) + 14N 15C + X could ostensibly dramatically speed up is allowed since (4BeX) is stable (EB 1.6 MeV; see → 8 ≈ certain steps in the CNO solar energy generation cycle in our Tab.II), in contrast to 4Be, which is famously unsta- own Sun. This is because the Coulomb barrier for the p to fuse ble (EB 92 keV) [96]. Such stable (BeX) can then onto the 14N nucleus only becomes apparent at distances on the fuse with≈ a further− He nucleus giving a pathway to (CX) order of tens of fermi [the (pX) bound state size], rather than that is unsuppressed by the three-body nature of the at a couple of thousand fermi [the normal classical turning point triple-α process; it is thus highly likely that even X still under the thermodynamic conditions prevalent in the Sun]. In principle this change could be expected to dramatically alter the bound as (HeX) after the earlier stellar burning phases relative neutrino yields from CNO vs. pp cycles to be expected get into (CX) bound states quite early via this ersatz from the Sun. However, this fails for subtle reasons. Many of the triple-α route, even before the full degenerate CO core steps in the CNO cycle have sufficient energy to directly eject the itself forms as a result of the helium flash. Such X will X in the final state; such an ejected X will find and bind with another nucleus, which is most likely a p (on number-abundance be unlikely to escape to other nuclei until the WD ther- grounds alone), which would allow a cycle of rapid fusions to mal runaway since the binding energy of the (CX) state, proceed. However, there is a regulator: on number-abundance 3 MeV, is much larger than the temperature in the CO arguments alone, once every 104 times the X is ejected, it will ∼ find and bind with an 16O nucleus∼ (present in 10 4 number- core (T . 10 keV; reaching T . 60 keV only in a near ex- ∼ − tremal CO WD just before thermal runaway [97]). More- abundance in the Sun), where it will remain until a proton fuses with the (OX) nucleus in a side-branch of the CNO cycle, possi- over, even though the degenerate electrons at the center bly ejecting the X in the final state. Unfortunately, we estimate of the near-extremal core WD do become relativistic as that the rate for that latter process is roughly as slow as the the Chandrasekhar mass is approached and thus have rate-limiting p + 14N step in the main uncatalyzed CNO cycle. Fermi energies that can exceed the (NX) binding ener- Therefore, given the small number of X in the Sun and the fact gies, they are of course still Pauli blocked from transfer- that any one X can only catalyze a limited number of reactions before being captured by a contaminant that keeps it bound up ring energy much larger than the temperature to (and and unable to catalyze further fusions for about as long as the thereby disrupting) the (CX) bound states. slowest CNO cycle step, we estimate that the overall impact on Although it would of course be energetically favorable the CNO cycle neutrino output is negligible. 14 likely be recaptured by another He nucleus, and the cy- occur during the active nuclear burning phases of the cle would repeat until the X got stuck in another heavy stellar lifetime do have reaction Q-values large enough nucleus such as C or O. If, on the other hand, such re- to potentially disrupt all the likely (NX) bound states actions occurred in the CO core after the He flash, the involved in various reactions, the X− will spend some X− would again simply be recaptured by another heavy portion of their lifetime during the active burning phases nucleus, likely C or O. All of which is by way of saying outside the protection of the positively charged nucleus, that once a X− is bound to C or O, it is highly unlikely and would during such periods be susceptible to annihila- + to be ionized again; or, if it is, it is likely to end up bound tion with X . However, the dilute distribution of the X± to another C or O nucleus. in the ordinary stellar matter (assuming fairly uniform In this paper, we will for simplicity therefore assume distribution in the WD-progenitor star) makes annihila- that the entirety of the X− population in the WD is tion a much less likely scenario than, e.g., the X− just be- bound to the heaviest nuclei that make up an (1) frac- ing immediately re-captured by another nearby nucleus. tion of the CO core of the WD-progenitor objectO that To judge the complete implausibility of any alternative forms at the center of the red giant: that is, they will outcome, consider a simple Γ nσv argument: if the ∼ be bound to either 12C or 16O, forming (CX) or (OX), CHAMPs were to constitute 100% of the local dark mat- respectively. ter abundance (a deeply excluded possibility), the maxi- We also note that Ref. [50] explicitly considers only mum mass of CHAMPs primordially present in the WD + for a m 1011 GeV CHAMP (which maximizes η in the net residual contamination of X and X− present X ∼ in the star after an assumed annihilation of the accumu- the regime where we know that the CHAMP halo abun- lated X+ and X ; that is, their result for the CHAMP dance is undisturbed by galactic dynamics) is of order − 13 contamination, which we claimed to have stated above MX 10− M , which yields an average number abun- ∼ 23 as Eq. (21), is actually given only in terms of the net dance in a Mwd 1.2M CO WD of nX 10− nion, where n is the∼ number density of the C∼ and O ions residual contamination η η+ η (i.e., per Ref. [50], ion η M /M is given≡ |by the− same−| expression as at in the CO WD. Even charitably allowing vX as large as X, net cloud 2 therm. ≡ the WD escape speed 10− 10v , it would still Eq. (21), but with f f+ f ). Since Ref. [50] ∼ ∼ C is ultimately interested± → in | CHAMPs− −| accumulating in require an enormous ratio of the CHAMP annihilation stars whose evolved cores later collapse to form neutron cross section to the radiative capture cross section of the 22 stars, in the cores of which the ambient temperatures are carbon ion on the CHAMP on the order of 10 for the T MeV [98] and the ambient densities are nuclear, this rate for an annihilation to be equal to the rate for ra- ∼ diative capture in any one instance when the X is thus is a justified assumption in their case as the X− are then − unable to bind to nuclei, and are free to annihilate with liberated in a nuclear reaction in the star. Of course, X+ when the CHAMPs collapse to an extremely dense the other way around this is if the X− participates in a state at the center of the star. truly enormous number of nuclear interactions in the life- However, we will be concerned with CHAMP contam- time of the star and therefore spends a large amount of ination of WDs, the interiors of which can only reach time outside the protection of the positively charged ion maximum temperatures (just prior to triggering thermal nuclear charge cloud; this appears exceedingly unlikely runaway in a near-extremal WD resulting in a full ‘ordi- on the grounds of the sheer number of nuclear reactions nary’ type-Ia SN) of T 7 108 K 60 keV (see Ref. [97] required for this to be an issue. Once again, detailed and references therein)∼ and× much lower≈ densities. This is evolutionary stellar modeling beyond the intended scope much smaller than the 3 MeV binding energies of typ- of this work would be required to fully resolve this is- & sue; but on the arguments advanced here, we will simply ical [e.g., (CX), (OX)] states into which X− are bound. + As a result, the (NX) bound states which are formed in state our results assuming that any X –X− annihilation the WD core will not be disrupted by the ambient con- which occurs during the active burning phases will be a ditions. Therefore, similar to the argument already ad- small effect. vanced for the diffuse galactic CHAMPs, the X− will be shielded (during the periods when they most densely accumulate) from annihilating against the X+ by the presence of the large Coulomb barrier between the Q = +1 + state X , and the Q = +(QN 1) 5–7 state (NX). B. Accumulated CHAMP contamination As such, we will not assume annihilation− ∼ of the symmetric CHAMP component, and continue to track both the + X and X− contaminations independently; see further An additional mechanism exists to populate WD with discussion in Secs. VI E 4 andVIIIC on this point and CHAMPs: accumulation of CHAMPs from the halo how it impacts the final limits we are able to set, which which pass within the gravitational capture radius of are in principle stronger than those of Ref. [50] in the the WD after its formation. In this subsection, we dis- CP -symmetric case. cuss the population estimate, estimate when this process An important caveat to the discussion in the preceding is efficient, and discuss the subsequent behavior of the paragraph is that since individual nuclear processes that CHAMPs. 15

1. Population estimate frame [99]: 2 1 v + vwd f(v) = − exp | | Θ[v v + v ] N − v2 esc, mw − | wd| Owing to their unit charge, even very heavy CHAMPs 0 can be stopped eﬃciently by stellar plasma, particularly (29) the extremely dense CO plasma that exists in WD in " 3/2 3 vesc, mw the appropriate mass range. Given the total WD mass = π v0 erf N v0 Mwd, radius Rwd, and fractional eﬃciency (v) in captur- " 2 ## ing a CHAMP with speed v in the WD rest frame pass- 2vesc, mw vesc, mw ing within the gravitational capture radius Rcapture(v) exp 2 , (30) − √πv0 − v 2 2 1/2 p ∼ 0 Rstar[1 + vesc, wd/v ] where vesc. = 2GMstar/Rstar is the WD escape velocity, the total accumulated X where v0 220 km/s is the local circular speed in the ± ≈ CHAMP masses obtained within an accumulation time MW; vesc, mw 540 km/s is the MW escape speed; and ≈ τaccum. are given by vwd is the WD velocity in the galactic rest frame, which we will take to have a magnitude vwd v0. We take ∼

(v) = Θ(vmax v), (31) M accum.(τ ) − X± accum. Z τaccum. Z where we estimate the maximum speed (in the WD rest 3 = dt d vf(v) (v) f ρX (t) frame), vmax, of a CHAMP far from the WD that will be- 0 ± come gravitationally bound to the WD after one passage " 2 # 2 vesc, wd through the WD as follows: a CHAMP hitting the sur- πRwd 1 + v, (28) 1 2 2 × v2 face of the WD carries energy Ei 2 mX vesc, wd + vX,0 3 ∼ where vX,0 . 10− is the CHAMP speed far from the WD in the halo (in the WD rest frame). This energy is lost to ions in the WD via elastic scattering events, each of which carries away a momentum of order ∆p m v where f(v) is the CHAMP velocity distribution far from ∼ ion rel R where vrel is the relative CHAMP-ion speed, and which the star (see below), normalized to ∞ 4πv2f(v)dv 1, 0 can be approximated by v max [v , v ]. In and the [ ]-factor captures the gravitational focusing≡ rel X ion, therm. order to become gravitationally∼ bound to the star, the enhancement··· to the geometrical capture cross section. CHAMP must lose enough energy that its energy drops below E 1 m v2 after traversing an average dis- We again assume that the symmetric CHAMP abun- f ∼ 2 X esc, wd tance of order Rwd within the WD. Approximating the dance is prevented from annihilating away by the large 3 WD as a uniform sphere of densityρ ¯ Mwd/(4πR /3), Coulomb barrier between (NX) and X+, so we will still p ∼ wd and noting that v = 2M /M 2 R 2 10 2 track the individual contamination abundances. Note esc, wd wd Pl. WD − while v v 10 3 (assuming T∼ ×keV), that f here should count the (pX) abundance, as the halo ion, therm. − we are in the∼ regime where∼ v v , and so∼ the the transfer− mechanisms of X from (pX) to (NX) are rel X − CHAMP energy decreases exponentially∼ with distance likely quite efficient given WD densities; this will be rel- travelled through the WD: E(x) E exp ( 2¯ρσx/m ). evant if we discuss accumulation of CHAMPs for masses i X The largest initial speed the CHAMP∼ can have− far from m 1011 GeV where it is possible that CHAMPs other X . the WD and still be captured after traversing an average than (pX) have been evacuated from the galaxy [52] (al- distance R 19 inside the WD is thus though, as we have already noted, this is subject to some ∼ wd controversy [49]). The contamination fraction of the WD 2¯ρσ 1/2 from this accumulation of halo CHAMPs is given by vmax vesc, wd exp Rwd 1 , (32) ∼ mX − η = MX /Mwd. ± ± We will assume that the star experiences a roughly 19 constant CHAMP density throughout its lifetime so that This is of course approximate; a more refined computation would solve for the trajectory of the incoming CHAMPs (specified com- ρ (t) const. Moreover, we will assume that CHAMPs X ≈ pletely by their initial speed far from the WD, and their initial are distributed in momentum space in the same way as direction of motion relative to the WD velocity) from infinite dis- fully virialized dark matter, which applies for mX & tance up to the surface of the WD, and then from that surface 1011 GeV [49, 52]; it turns out (see Sec.VIII) that this onward as they pass through the WD losing energy, following a is the only mass range in which these additional ac- non-conservative trajectory which would likely need to be solved for numerically; such a computation could also make use of the cumulated CHAMPs are relevant for setting limits, so true density profile of the WD to refine the estimate. We expect this is not a very strong additional assumption. More however that such a computation would yield results that are specifically, we will take f(v) to be given by a truncated not dramatically different from the result we give here, and we Maxwellian speed distribution in the MW galactic rest have not performed such a refined computation. 16

where σ is the total cross section for momentum trans- are required to take account of the maximum stopping fers of order ∆p mionvrel, which we very conserva- speed, which imposes an upper bound v vmax. of the tively20 estimate to∼ be approximately the geometrical nu- dv integral in Eq. (28), as implied by Eq.≤ (31). 2 clear cross section σ π 1.22 fm A1/3 170 mb The final result of a careful treatment of the integrals in 200 mb, with A 7∼ consistent with× the mean-ion∼ ap-∼ Eq. (28) is a set of complicated functions for the accreted proach we utilize∼ throughout; see Ref. [59] for a more CHAMP abundance, which we show in AppendixC; we detailed discussion of cross sections for stopping in WD use the full expressions Eqs. (C1)–(C3) for all quantitative results in this paper. It is however useful to develop a material. Demanding that vmax & vesc, mw + v0 conservatively guarantees that the CHAMP will become bound simple approximate estimate for the accreted mass; this after one passage through a distance Rwd of WD mate- can easily be done directly from Eq. (28) in the case of m mmax.. Suppose we ignore the fact that the rial; a conservative estimate for the maximum CHAMP X X mass that is thus guaranteed to become bound to the WD CHAMPs have a velocity distribution, and approximate is given by f(v) = δ(v v0)/4π, and also take = 1; then we have that − max 2¯ρσRwd accum. approx, III mX M (τaccum.) M 2 X± X± ∼ vesc, mw+v ≈ ◦ 2 ln 1 + v esc, wd 2 vesc, wd τaccum.f ρX πRwd (34) σ ≡ ± v 1.5 1016 GeV , (33) 0 ∼ × × 200 mb 2 2 where we taken vesc, mw v0. Numerically, this estimate is a factor of roughly 1/erf(1) 1.2 larger than the full where we again remind the reader that the cross section ∼ assumed is conservatively small given that the CHAMPs result given at Eq. (C3) [which is well within the uncer- are electrically charged. While Eq. (33) is a good es- tainty on the estimate itself]; normalized to this approx- timate for where these effects will begin to become rele- imate result M approx, III, we plot the full result Eq. (C3) X± vant, an (1) fraction of CHAMPs that traverse the WD in Fig.2. The fractional suppression of the accumulat- will stillO be efficiently captured even for CHAMP masses ing CHAMP abundance at large mass scales21 roughly as max 1 somewhat above mX ; see Fig.2. mX− . For CHAMPs lighter than mmax, it is safe to assume Suppose we consider either a M 0.85M WD or X ∼ = 1 in Eq. (28) and to perform the integral without a M 1.1M WD, take the accumulation duration to ∼ regard to the maximum stopping speed; this actually suf- be the timescale for cooling of a WD to the point of fices for the majority of the mass range where we will be crystallization (see discussion below) τ 109 yr accum. ∼ interested in the accumulation mechanism (see discus- [100], and normalize to the local MW DM density ρ halo ∼ sion below and in Sec.VIII). For heavier CHAMPs, we 0.3 GeV/cm3, then for m 1.5 1016 GeV, we have X . × 13 f ρX η , accum.(0.85M WD) 1.9 10− ± (35) ± ∼ × ρhalo 20 While energy loss by charged particles in dense, degenerate 13 f ρX WD material is in general a fairly complicated problem [59], an η , accum.(1.1M WD) 1.3 10− ± , (36) ± alternative estimate for the linear stopping power of the WD ∼ × ρhalo material can be obtained from considering elastic electromag- which exceeds the CHAMP contamination fraction ob- netic scattering of the CHAMPs off the CO ions, represented tained via the protostellar channel (for a Mcloud = 7M by a Thomas–Fermi screened potential [59]. For particles which are on the borderline between becoming bound to the protostellar cloud, as appropriate for a Mwd 1.1M 22 11 ∼ WD and not becoming bound, we would have vX vesc, mw WD [91]) for mX 1.6 10 GeV, assuming for the fairly constant within (1) factors during the traverse∼ of the ∼ × O3 2 WD since vX,0 10− is smaller than vesc, mw 2 10− ; as such, the Coulomb∼ stopping power can be estimated∼ × as 2 2 2 21 R vmax 2 2 2 1 dE/dx 2πn Z α / m v Λ[59] where Λ The integral scales as dv v v v− vmax m− , |Coulomb ∼ − ion ion ion X × 0 × × ∝ ∝ X is a Coulomb log of order Λ 10–15 for the relevant parameters. where the factor of v2 comes from the integration measure, the ∼ By comparison, the linear stopping power we have estimated with factor of v comes converting the CHAMP abundance to a flux, the geometric cross section is dE/dx n m σv2 . and the factor of v 2 comes from the enhanced gravitational |geometric ∼ − ion ion X − The ratio is (dE/dx ) / dE/dx capture radius for slow-moving CHAMPs. The final scaling with |Coulomb |geometric ∼ 2 2 2 4 mX is obtained by expanding the exponential in Eq. (32) for 2πZionα / mionvX σ Λ. Consistent with the mean- × 2 16 ion approach, if we take Zion 7, mion 7µa, vX 2 10− , mX 1.5 10 GeV [see Eq. (33)]. ∼ ∼ ∼ × 22 × and Λ 15, we find (dE/dx Coulomb) / dE/dx geometric Note that while the mass of the WD progenitor star will be ∼ 1 | | ≈ 70 (σ/200mb)− . Our geometrical cross section estimate is comparable to the collapsing protostellar cloud mass we input, thus× perhaps too conservative by about a factor of 70. Taking significant shedding of the outer envelope of a highly evolved red the more aggressive cross section estimate would increase giant occurs after the helium flash reduces the associated WD the maximum stopping mass linearly by the same factor; see progenitor CO core mass to something M [95]; see Ref. [91] ∼ Eq. (33). The only change to the limits would be that the for a discussion of this initial-mass–final-mass relationship. Our correction factor estimated in Fig.2 would remain 1 until results are actually relatively insensitive to the exact assumed ∼ 1/3 m 1 1018 GeV, so the limits we present in Figs.4 and5 value of M , as η M − for large m [see Eq. (21) for X ∼ × cloud ∝ cloud X would be slightly stronger at the largest masses we consider. ycrit. < ymin]. 17

2. Eﬃciency of accretion 100 1

± We note that in order for CHAMPs accumulated by mX− approx, III X ∝ this mechanism to sink to the center of the WD and pos-

/M sibly feed a dense central core object, we should assume ) 10 1 that the accretion is occurring while the interior of the − WD is still in a liquid/non-solid phase, otherwise the accum.

τ sinking of the X± is likely significantly inhibited (see, ( e.g., Ref. [100], in which the diffusion co-efficient for 22Ne

± 2 contamination in a CO WD is set to zero after crystalliza- accum.

X 10 − 1016 1017 1018 tion). Since WD cooling models (e.g., Ref. [101]) indicate M 9 mX [GeV] that that crystallization only occurs at ages τ & 10 yr for CO WDs with masses around MWD 0.85M , we restrict our attention to CHAMP accretion∼ occurring FIG. 2. The CHAMP abundance that becomes gravita- within at most the first τaccum. Gyr after formation. accum. tionally bound to a WD, M (τaccum.) [Eqs. (C1)–(C3)], ∼ X± Note that it is not necessarily the case that sufficiently compared to the na¨ıve estimate M approx, III given at Eq. (34), X± heavy CHAMPs would be prevented from sinking even when the CHAMP speed distribution and maximum stop- after crystallization,23 so this is conservative. ping speed are taken into account. This plot assumes that Note that a short initial interval of CHAMP accumu- the CHAMP passes through a distance R of WD ma- wd lation may be complicated by the persisting existence terial with a density equal to the average∼ WD density for a of an ionized planetary nebula blown off by the late- Mwd 1M WD with Rwd 0.01R , and that the CHAMP velocity∼ distribution far from∼ the WD is given by a truncated stage evolution of the WD progenitor star. Given that Maxwellian distribution, Eq. (29). We assume that the WD planetary nebulae do not remain ionized for much more moves with respect to the CHAMP distribution at the local than few 104 years (see, e.g., Ref. [102])—orders of ∼ × circular speed in the MW, vwd v0 220 km/s, and that the magnitude shorter than the crystallization time for the ∼ ∼ MW escape speed is vesc, MW 540km/s. We further assume WD—this will likely be a negligible effect; we have not ∼ that σ = 200 mb; see text. Note that the accumulated mass attempted to account for it. only begins to significantly deviate from its maximum value 16 max. Relatedly, the results of Ref. [49] indicate accumulation once mX & 3 10 GeV 2mX ; because the estimate of CHAMPs (here, onto the WD) in the presence of stel- × ∼ all at Eq. (33) was conservative in requiring CHAMPs to ac- lar winds and magnetic fields can be complicated by the crete, a significant fraction still accrete even somewhat about this mass. entrainment of inwardly diffusing CHAMPs in magnetic field lines moving outward with the charged wind. How- ever, an absence of solar winds is expected for cooling 7 2 WD with sufficiently high surface gravity g & 10 cm/s 4 sake of argument the same ρX (and ignoring any differ- and Teff . 5 10 K surface temperatures [103]; the sur- × 8 2 ences in (pX) accumulation). face gravity of a Mwd & 0.8M WD is g & 2 10 cm/s × Note that this estimate will be subject to large Pois- (and higher for more massive WD, which have smaller son uncertainties once the CHAMP mass becomes large radii), so as long as the WD we are interested in have enough that, given the CHAMP density, the number of low enough surface temperatures, no winds should be CHAMPs that pass through the capture area in the stel- expected. accum. accum. On the other hand, certain WD are known to possess lar lifetime, N = M /mX , becomes small. X± X± ± high magnetic fields [104] (see also Ref. [67] and refer- √ Demanding conservatively that the fractional 1/ NX± uncertainty in the estimate is less than 10% demands 2 that NX± & 10 . In order not to be subject this Poisson uncertainty, we must require that 23 If the differential gravitational force acting on the CHAMP and a neighboring ion exceeds the electrostatic repulsive force exerted by the ion on the CHAMP at approximately the lattice spacing  38 mX 17 then, in a lattice configuration, the CHAMP would likely sink 6 10− 105 GeV mX . 10 GeV f ρX  × through the lattice. To get an idea of the mass scale involved, ± & 2 , + m consider an X located in a cold Mwd 1.2M WD, approxi- ρhalo Poisson 22 X 18 1.2 10− 2 mX & 10 GeV ∼ 1 /3 mPl. mate the lattice spacing by a (3mion/(4πρ)) , and approxi- × 3 ∼ (37) mate ρ 3Mwd/(4πRwd); taking mion 12µa and Qion 6, we ∼ ∼ 17 ∼ find that this criterion is reached for mX & 10 GeV, so while the very heaviest X+ CHAMPs could potentially sink through a with a small transition region in the scaling in the solid lattice in a WD, lighter CHAMPs may be prevented from doing so. The cognate mass estimate for X− bound as (NX) is intermediate mass range. See further discussion in higher by the charge of the (NX) state. Note however that this Sec.VIII where we consider the limits the accumula- estimate is very rough, which is why we conservatively truncated tion of CHAMPs imposes on the galactic abundance of our accumulation at a timescale such that the WD has not yet CHAMPs. crystallized. 18

ences therein). Moreover, even if the X+ or (NX) bound Finally, note that if all the above caveats about this state are neutral by virtue of having captured electrons, process being efficient are satisfied then, until such time we estimate that the UV luminosity of even an old WD as the WD crystalizes, this accumulation of CHAMPs is sufficient to ionize at least one or more electrons from is in addition to any primordial abundance of CHAMPs the X+ or (NX) state;24 as a result, the WD magnetic that may already have been present in the WD. fields could deflect incoming CHAMPs, making the accumulation estimate too aggressive. However, this depends on the rigidity of the CHAMPs: we estimate that if a 3. Behavior of accumulated CHAMPs + + 11 X or singly-ionized (NX) of mass mX 10 GeV p ∼ + moving at a speed v 2GMwd/(2Rwd) were to expe- An X incident on the WD is simply stopped, then ∼ rience a 7.5 MG magnetic field roughly one WD radius diffusively sinks to the core of the WD (see the follow- ∼ from the surface of a Mwd 1.1M WD, its Larmor ing sections). The fate of an X− may be different: if ∼ radius rL p/qB would be of order the WD radius, it is incident as (pX), it can first be stopped and then ∼ which would likely result in the particle being deflected rapidly undergo an exchange reaction, becoming bound enough to miss the WD surface. Much weaker magnetic to a heaver nucleus (NX), where N is either He (in the fields at this radius would not deflect the particle suffi- WD atmosphere, if any) or is most likely C or O, which ciently to cause it to miss the surface; moreover, since any form the bulk of the WD. Alternatively, the exchange re- magnetic field drops off at large distance from the WD action can occur during the stopping process. In either 3 at least as rapidly as a dipole field B r− , the Lar- case, after exchange, the heavier nucleus bound state will, ∝ mor radius at large distance grows much more rapidly after perhaps having undergone some further stopping, than the distance from the WD, so it would likely not merely diffusively sink to the core of the WD. If the in- cause the particle to miss the WD owing to its effect coming X are in bound states (NX) where N is heavier further out from the WD. Although it is not necessar- than He, it is possible that the X remains bound to that ily the case that fields much stronger than the one es- nucleus, and merely gets stopped and sinks. timated here would cause the particle to miss the WD either (in a very strong field, the CHAMP would likely simply spiral down the magnetic field lines, still likely C. Sinking, stratification, and timescales hitting the WD surface) this is a much more challenging situation to analyze. We will thus wish to consider + CHAMPs, whether X or X− bound as (NX), have a only WD with low magnetic field 7.5 MG when con- much smaller charge-to-mass ratio Q /m than the CO 11 X X sidering mX & 10 GeV CHAMPs accumulating onto a material in which they are interspersed when the WD is Mwd 1.1M WD. On the other hand, if we were to born. As such, they would be expected to sink (diffu- consider∼ a m 3 107 GeV CHAMP (about as light X ∼ × sively) in the WD toward the center of the star. More- a CHAMP as we will concerned with in the accumulat- over, once formed, a WD experiences no internal large- ing case) accumulating onto a Mwd 0.85M WD, a B ∼ scale nuclear burning processes that would trigger con- field of 1.1kG at one WD radius above the surface would vection until/unless thermal runaway is triggered, so con- be problematic; very low magnetic field WD are thus re- vective mixing of the WD contents is unlikely to disrupt quired for bounds based on such accumulating WD to be this sinking. In this subsection, we discuss the timescale robust. for the sinking process. A very rough argument for why q/m is the relevant quantity [100, 105–107] when considering sinking is that the (non-degenerate) ions present at some radius in the 24 3 Consider a WD with a luminosity of Lwd 10− L , which ∼ WD would sink compared to the electrons absent a very is typical for a Mwd M WD less than Gyr old [66], and conservative for the∼ luminosity at an earlier∼ age in the small net positive charge overdensity interior to that ra- WD existence since WD cool over time. The surface temper- dius which, in hydrostatic equilibrium, supplies exactly ature is then of order T eV, and so the fraction of the the correct electric field to balance the gravitational force number of photons emitted∼ from the WD that have an energy on the ions at that radius. However, since FE q and sufficient to ionise the outer electron from the CHAMP state ∼ Fg m, such a bulk electrostatic force balance can only (which we take to be 10 eV) is, conservatively, f>10 eV 3 ∼ ∼ ∼ 10− . Approximating the hydrogenic photoionization cross sec- work for one value of q/m: particles with smaller than 2 11 tion as σ 4παrB where rB 5 10− m is the Bohr average q/m at a given radius must sink, while those with radius, we∼ find that at a distance∼ of×R = 10R from the 12 16 wd a larger q/m must rise. Indeed, the 6 C– 8 O(q/m = 0.5) WD surface, the timescale for an ionization interaction to oc- 1 2 2 mixture in a WD does not stratify (while it remains liq- cur is τionize− nγ σ (f>10 eVLwd/(4πR )/T )(4παrB ) 22 10 ∼ ∼2 ∼ uid), but contaminants, e.g., Ne (q/m = 0.45 < 0.5) 10− (L /eV)(rB /Rwd) 1/(0.05s). However, even moving 10 ∼ at escape speed, a CHAMP would take a time t R/vesc sink on cosmologically long timescales [100, 106, 107]. q ∼ ∼ 2 Here we follow the discussion of Ref. [100, 107] (see R RwdM /2Mwd 8 s τ to traverse the remaining Pl. ∼ ionize distance to the WD. Therefore, we estimate that the UV luminos- also Ref. [105, 106, 108, 109]) to provide an estimate of ity of the WD is sufficient to ionize electrons from the CHAMP the timescale for the sinking. A first very approximate (1) state at least a distance R & 10Rwd from the surface of the WD. estimate for the timescale τ for sinking an (1) frac- sink O 19

(1) 21 27 tion of the WD radius is τsink Rwd/wX (Rwd), where we then estimate is [clearly for mX 10 GeV] w (R ) is the CHAMP diffusion∼ velocity at the sur- X wd m face of the WD. Assuming that the CHAMPs are only a 7 X wX (Rwd) 1 10− m/s 5 (44) trace constituent in the WD background (easily satisfied ∼ × × 10 GeV m during the initial sinking phase), and that the charge-to- 16 X 3 10− c 5 (45) mass ratio of the X is very small, the (terminal) diffusion ≈ × × 10 GeV 25 m 1 velocity can be estimated as τ (1) 1.6 106 years X − . (46) sink ∼ × × 105 GeV

3 For mX & 10 GeV, this is less than the old WD lifetimes w (r) gDm /T, (38) | X | ≈ X and/or crystallization times of Gyr we consider, but the estimate here is crude because∼ we have taken g and hence wX to be radially independent. 2 2 where g = M(r)/(r MPl.) is the local acceleration A better but still highly approximate analytically due to gravity, and D is a diffusion coefficient, which tractable estimate is obtained by taking into account Refs. [100, 107] indicate can be estimated as the self- the radial dependence of wX , making the assump- diffusion coefficient for the CO plasma: tion that the WD is a sphere of uniform density (see, e.g., Ref. [50]). Then, we have that g(r) = 2 2 3 2 2 Menc.(r)/(r MPl.) = Mwd(r/Rwd) /(r MPl.) = 2 4/3 D 3ωpa Γ− (39) g(R ) (r/R ), implying that ≈ wd × wd r wX (r) = wX (Rwd) , (47) where × Rwd so that the timescale estimate to move from radius r = R R to r = R becomes 2 2 2 i ∼ wd f ωp = 4πZionρα/mion (40) 1/3 Z Rwd dr R a (3mion/4πρ) (41) (2) (1) wd ≡ τsink = τsink ln , (48) 2 ∼ Rf wX (r) Rf Γ = αZion/(aT ), (42) 1/9 1/9 4/3 which is longer than the previous estimate by a logarith- 9 3 mion T D , (43) mic factor that depends on the final radius. In the next ⇒ ≈ 2π11/18 4 Z5/3α5/6ρ11/18 ion subsection, we estimate the initial radius of the CHAMP structure that forms at the core of the WD as a result of the sinking, which will show that this logarithmic factor where ‘ion’ here refers to an ion in the mixture through is never large enough to cause the sinking timescale to which the CHAMPs are sinking. We will adopt the become unacceptably long. widely used ‘mean ion’ approach and set m 14µ ion ∼ a Moreover, as we are also interested in the case of and Zion 7, assuming roughly equal abundances of C 26 ∼ sinking where the star also contains an extremely com- and O. pact core object (i.e., a BH) at the center of the For the purposes of an initial rough estimate, we WD, it is worth considering how the above estimate 3 will take T 107 K 1 keV (WD are approximately is modified in the case where Menc = Mwd(r/Rwd) + isothermal, and∼ this is∼ a typical WD core temperature), 3 Mwd 1.1M , Rwd 7 10− R , and we will approximate∼ the density with∼ the× average WD density: ρ ρ¯ 3 6 3 ≈ ≡ 27 We note that the terminal velocity estimate, and hence sink- 3Mwd/(4πRwd) 4.5 10 g/cm . The sinking timescale ≈ × ing time, are dramatically different from the na¨ıve estimates obtained from setting the viscous drag force F ρ v σw drag ∼ ion ion X equal to the gravitational force [50], which would yield wX p ∼ gm /(ρ v σ). If we take v T/m and estimate X ion ion ion ∼ ion σ 200 mb as we have throughout, we would obtain wX ∼ 5 ∼ 15 cm/s (mX /10 GeV), assuming the same average WD density as in× the main text. The difference is likely ascribable to the 25 Note that there is no QX dependence in wX . For X a trace WD interior being a strongly coupled plasma which acts like a element of mass mX sinking through a background plasma with liquid, instead of a rarefied gas, so that the low nuclear cross q/m = 2, the full expressions of Ref. [100] give w = 2Q section estimate we have used up to now (always∼ in a way thus | X | | X − m /µa (mpgD/T ), which yields Eq. (38) in the limit m µa far that was conservative) potentially supplies a dramatic under- X | X regardless of QX (if we approximate mp µa, which is allowed estimate of the viscous drag force on the CHAMPs, and yields at the level of accuracy of these computations).≈ a much too aggressively short estimate of the sinking time. The 26 At this level of approximation, it is irrelevant whether we as- extreme discrepancy between these estimates does however give sume the number or mass abundances are equal. Technically, some pause, and we note that the timescales we give here are this assumes the number abundances are equal. possibly conservatively long. 20

3 2 Mcore 1 (r/Rwd) for Mcore Mwd. This implies for Rf . 10− Rwd, that that − Z Rwd (4),1.1M dr τ (54) sink ∼ w (r) " Rf X r M w (r) = w (R ) 1 core (1) ρ¯ X X wd 2τsink Rwd − Mwd ∼ × ρC 2 # ( " 1/3#) M R Rwd ρC Mwd + core wd (49) min ln , ln Mwd r × Rf ρ¯ Mcore (55) m 1 4 106 years X − If Mcore Mwd, then for r & rcross ∼ × × 105 GeV 1/3 ≡ ( " 20 1/3#) Rwd (Mcore/Mwd) , the dynamics are still dominated 1 ρC /ρ¯ 10− M by the WD material enclosed at radius r, and the esti- 1 + ln × 13 13 Mcore mate at Eq. (48) holds; however, for r . rcross, the dy- (56) namics are dominated by the core, implying that wX 2 ∼ r− , which regulates the logarithmic divergence in the assuming always that Rf Rwd, and where we took total sinking time estimate at Eq. (48): 3 ρC /ρ¯ 13 for this WD. Again, for mX & 10 GeV this is less≈ than WD lifetimes and/or crystallization times of Gyr. Similar estimates hold for WD throughout the Z Rwd (3) dr range∼ of masses in which we are interested in this paper. τsink (50) ∼ Rf wX (r) Note that the logarithmic factor shown at Eq. (56) dif- ( " 1/3#) fers from that at Eq. (51) because the increased central (1) Rwd Mwd density of WD material implies a smaller cross-over ra- τsink min ln , ln . ∼ × Rf Mcore dius if a massive core is present. Finally, note that the (51) core mass assumed here is, within an order of magnitude 10 or so, the Chandrasekhar mass for a mX 10 GeV CHAMP, and that the logarithmic factor is∼ clearly not The maximum total time to sink to the center of the star very sensitive to the exact assumed core mass. in this case is thus We must also consider whether it is actually possible to achieve the terminal diffusion velocity under the gravitational acceleration prevailing in the WD; if not, " # then the above estimates could be incorrect. To see that M 1/3 (3) (1) wd this is easily possible, consider that the maximum ve- max τsink τsink ln (52) ≈ Mcore locity achievable by a particle free-falling through a uni- 1 form density sphere of mass M and radius R is given 7 mX − p 2 2.5 10 years 5 up to (1) factors by v M/(RM ). For ∼ × × 10 GeV O max. ∼ Pl. ( " 20 1/3#) the Mwd M WD discussed above, this estimate is 1 10− M ∼ 2 1 + ln , (53) vmax. 10− c, which is much greater than the termi- 15 M ∼ × core nal diffusion velocity for all mX . MPl.; the terminal velocity is thus always reached. Note that we will conservatively elect not to consider 3 3 CHAMPs lighter than mX 10 GeV, where the sinking which for mX & 10 GeV is less than WD lifetimes . and/or crystallization times of Gyr. time estimates given here approach WD lifetimes and/or ∼ crystallization times of Gyr; see Figs.4 and5. Moreover, Eqs. (48) and (53) are still over-estimates of ∼ the sinking timescale because the density of the WD increases significantly above the average density as the core D. Thermally supported CHAMP-contaminated is approached, which makes Menc(r) and hence wX (r) WD structure: self-gravitating collapse and timescales (in the region r > rcross) larger than that estimates assuming the uniform sphere. Indeed, for a WD of mass + Mwd = 0.8–1.2M , the density exceeds the average den- The X and (NX−) present will sink diffusively to the sity for r/Rwd . 0.6–0.7, with the central density even- center of the newly formed WD in a characteristic time tually reaching a value roughly a factor of 9–16 larger τsink (see Sec.VIC), until (or if) they encounter suffi- than the average density. Recomputing the timescale es- cient pressure to halt this collapse/sinking process and timate with a realistic WD density profile, we find that establish a stable hydrostatic equilibrium (should such for a Mwd = 1.1M WD we have, within a factor of 2 an equilibrium exist). In this subsection, we consider the ∼ 21

initial formation of a CHAMP structure at the center of their Fig. 1, ‘with neutrinos’), corresponding to a core the WD; see also the discussions in Refs. [50, 60, 61]. temperature in the range T 107–106.3 K, or T 0.8– Initially, so long as the WD material still dominates 0.2 keV. ∼ ∼ the central mass density of the WD, the CHAMPs Note that r also supplies an alternative natural cut- will form an approximately isothermal thermal-pressure- off to the logarithmic∗ divergence in Eq. (48): even for 28 18 supported structure at the center of the WD mX 10 GeV (as massive a CHAMP as we con- 2 ∼ 4 11 ρX (r) ρX (0) exp[ (r/r ) ], with a characteristic scale sider; see Sec.VIII), r 10− m 10− Rwd (for ∼ − ∗ 3 ∗ ∼ ∼ height [50] Rwd 7 10− R ), yielding a logarithmic enhance- ∼ × ment of the sinking time by a factor of log(r /Rwd) 25; s ∗ ∼ 3TM 2 the sinking time estimate Eq. (48) for a CHAMP of this r = Pl. (57) mass thus becomes 102 s. On the other hand, for ∗ 2πmX ρwd 5 ∼ mX 10 GeV, we have log(r /Rwd) 10, so the sink- 1/2 ∼ ∗ ∼ T m 1/2 ing timescale Eq. (48) is 107 years. In both cases, these 350 m X − ∼ ∼ × 1 keV × 105 GeV timescales are sufficiently short, and are within a factor of 5–10 of the shorter estimate Eq. (56). ! 1/2 ρ − This∼ thermal-pressure-supported structure is stable as wd , (58) × 5.5 1010 g/cm3 long as ρX (0) . ρwd; at densities above this, the CHAMP × structure will begin to dominate the mass density at and total mass the center of the WD, and the CHAMP structure becomes unstable to a collapse mode in which the CHAMP X 3/2 3 M = π ρX (0)r (59) configuration loses total energy, heats up, and contracts ∗ ∗ 3/2 [50, 61], also known as the ‘gravothermal catastrophe’ 12 ρX (0) T 29 6 10− M [54]. This implies a maximum stable mass for the core, ∼ × × ρwd × 1 keV 3/2 ! 1/2 2 ! 3/2 − mX ρwd 3/2 3 3TMPl. − Ms.g., X = π ρwdr = (61) 5 3 ∗ 1/3 × 10 GeV × 5.5 1010 g/cm 2mX ρwd × (60) 3/2 12 T 6 10− M ∼ × × 1 keV where ρwd is the central WD density, for which we have 3/2 taken the fiducial value for a Mwd 1.1M WD (see mX − Tab.I). T is the central WD temperature,∼ which we con- × 105 GeV servatively estimate as T 1 keV (a higher temperature ! 1/2 ρ − increases the self-gravitating∼ mass, and would weaken our wd 3 (62) constraints; see Sec.VIII): Ref. [110] indicates that a × 5.5 1010 g/cm 8 × Mwd 1M WD with an age in the range τwd 10 – 9 ∼ 2.25 3∼.5 10 yrs has a luminosity Lwd/L 10− –10− (see [i.e., Eq. (60) with ρX (0) = ρwd]; core masses above this ∼ value will collapse spontaneously (see discussion in following sections). The timescale for this collapse could be limited ei- 28 There is a significant caveat to this discussion. For the picture of ther by the timescale for the CHAMPs to lose en- the thermal-pressure-supported structure that we advance in this ergy to the WD material to allow the collapse to oc- section to be correct, the X must contribute an ideal gas term cur, or by the drift time for the CHAMPs to sink in- to the pressure P n (r)T and this must be the only term ⊃ X wards from r . The energy transfer timescale is ex- in the pressure acting on the CHAMPs that varies significantly ∗ over the length scale r (only pressure gradients hydrostatically tremely short: assuming that vX vion, our discus- ∗ support structures against gravitational collapse). The X± are sion of energy loss in Sec. VI B 1 is applicable, and the however electrically charged and immersed in strongly coupled, energy loss rate for a single CHAMP can be estimated charged non-degenerate CO and degenerate electron plasmas, so this assumption is likely a gross approximation. It is therefore possible that the X are not stalled in this thermal structure, but instead simply continue to diffusively sink toward the center of the star, directly forming the denser core structure at the center 29 A virial population of N self-gravitating particles with average of the WD that we discuss in Sec.VIE, which has a maximum total kinetic energy EK , and temperature T EK /N obeys h i ∼ h i 2 Chandrasekhar mass before it too must collapse. However, be- the virial theorem EK (1/2) U where U GM /R cause we will set limits (see Sec.VIII) requiring the presence is the average totalh potentiali ∼ − energyh ofi the particlesh i ∼ if − their to- of the larger of the self-gravitating mass or the Chandrasekhar tal mass is M and the characteristic radius of the configuration mass of CHAMPs in the core structure, and because the self- is R. The average total energy of the system is thus E 2 h i ∼ gravitating mass exceeds the Chandrasekhar mass for large mX EK + U (1/2) U EK NT GM /(2R). The [Eq. (68)], it is conservative to assume that this thermal struc- hconfigurationi h i ∼ thus hash i negative ∼ −h heati ∼ − capacity∼ −dE/dT N, so ture must form and become self-gravitating before the CHAMPs an energy decrease causes a temperature increase and a∼ decrease − can sink further. in R for fixed M, implying a contraction. 22

2 as dE/dt ρionσvX vion. The characteristic timescale We begin in this subsection by considering the case for the CHAMP∼ − to lose (1) of its kinetic energy is where degeneracy pressure can re-stabilize the collapse 2 O thus τenergy mX vX /(dE/dt) mX /(ρionσvion) of the thermal structure, examining the impact on the 16 ∼ 5 ∼ 1 ∼ 1.5 10− s (mX /10 GeV) (σ/200 mb)− , where we mechanical structure of a WD and explaining important used× the same× cross section estimate× as in Sec. VI B 1; evolutionary changes in the CHAMP chemistry that oc- see also Ref. [50, 61]. cur. On the other hand, if we track a particle collapsing inward with the collapsing distribution of CHAMPs of + fixed total mass Ms.g., X from r to smaller distances, we 1. Positively charged CHAMPs: X find that the timescale for infall∗ assuming the estimate 30 of the diffusive sinking speed given by Eq. (38) with We discuss first the X+ case. At typical (non- 2 2 7 g(r) = Ms.g., X /(r MPl.) is given by extremal) CO WD core temperatures of T 10 K 1 keV, the X+ have thermal de Broglie wavelengths∼ ∼ Z r 2 " 3# p 15 p 5 ∗ dr TM Rf λ 2π/(m T ) 2 10− m 10 GeV/m τ core = Pl. 1 th. dB ≡ X ∼ × X coll. 3/2 which are much smaller than the typical inter-ion spac- ∼ Rf wX (r) 3π ρwdmX D − r ∗ ing δ (3/(4πn))1/3 2 10 13 m. Here, we have used (63) − the particle≡ number density∼ × at the core of an M M 2 31 3 ∼ + TMPl. WD, nion 3 10 cm− , neglecting that the X 3/2 [Rf r ] (64) ∼ × ≈ 3π ρwdmX D ∗ themselves will impact the mechanical structure of the m 1 WD, increasing the central number density and decreas- 1.5 105 yr X − , ∼ × × 105 GeV ing δ. Nevertheless, these estimates indicate that the (65) central density would have to be increased by a fac- 6 5 3/2 tor of 10 (mX /10 GeV) over that for a nor- indicating that the diffusive sinking time is by some or- mal WD∼ before× X+ degeneracy begins to be a concern. ders of magnitude the limiting timescale for all masses Such density increases are however not observed to oc- 18 of interest to us (mX . 10 GeV; see Sec.VIII). cur in numerical solutions of WD structure (see Ap- Note also that this is much shorter (by a factor of pendixD): assuming that the CHAMPs form a strati- 25) than the total diffusive sinking time Eq. (56) for fied core structure inside the WD, the maximum sup- CHAMPs∼ distributed throughout the WD volume, so if portable number density in that core structure scales prim. MX Ms.g., X , the core shrinking process completes up from the ordinary WD density as n mX as the more rapidly than any remaining CHAMPs distributed core becomes more and more extremal, causing∝ the inter- throughout the WD in excess of M sink to the core. + 1/3 s.g., X X spacing to decrease as δ mX− . Assuming 31 3 5 ∼ 35 3 We note in passing that if the energy released in the n 3 10 cm− (10 GeV/mC ) 3 10 cm− ∼ × ×14 ∼ × self-gravitating collapse is sufficient to heat the star above implies that δ 10− m, so the de Broglie wavelength the trigger criteria, this is an opportunity for early de- computed above∼ is still smaller than the inter-ion spacing 5 struction of the WD. at the center of the denser core for mX = 10 GeV. More- The endpoint of the self-gravitating collapse phase de- 1/2 over, because λth. dB mX− , which is faster than the pends on the total mass of CHAMPs and the CHAMP inter-ion spacing decreases∝ under these circumstances, mass; we survey the possibilities in the following subsec- this hierarchy persists to higher mX , and we can thus tions. safely take the X+ themselves to be non-degenerate in the WD at all times. Being charged, the X+ are however electromagneti- E. Degeneracy supported CHAMP-contaminated cally tightly coupled to the highly degenerate electron WD structure: sub-Chandrasekhar case e plasma (for a typical M M CO WD, λth. dB 40δe), which thus supplies electron∼ degeneracy pressure∼ support Assuming that the self-gravitating collapse discussed to the X+, in exactly in the same fashion as it does for the in the previous subsection proceeds, the endpoint of that ordinary positive C and O ions in a WD. The electron de- 4/3 collapse depends on the total mass of CHAMPs present generacy pressure scales as Pdeg. ne in the relativistic in the star, and also the CHAMP mass (see Sec. VI E 3 limit,31 but is independent of the∝ mass of the heavy ion for a detailed discussion of when the discussion in this to which this pressure support is communicated. At the section is applicable).

31 5/3 The degeneracy pressure scales as ne for non-relativistic elec- 30 This estimate is still appropriate to use in this context, provided trons, which increases parametrically even faster with increas- that mX GeV, since the CHAMPs can still dominate the mass ing ne; the electrons at the center of even a M M WD al- ∼ density in the core while being a trace element at the onset of ready have EF = 2.5me, which is already into a fairly relativistic the sinking. regime, where E scales approximately linearly with increasing p. 23

onset of self-gravitating collapse, when the central mass which must then be displaced upwards to maintain over- fraction of the CHAMPs and CO mixture are approx- all charge neutrality of the plasma (this is again in con- + imately equal, XX XC/O, the X constitute only a trast to the neutron star case, where both the neutron small number fraction∼ of the central ions: x x star matter and the X can co-exist without one displac- X C/O because mX mC/O. As such, the mean molecular mass ing the other [50]). Since the mean molecular weight per + per free electron [Eq. (D10)] is almost unperturbed by the free electron µe [Eq. (D10)] of the X /e− fluid mixture presence of the CHAMPs, and the electron density and is larger than for the CO mixture, this core will have a spatial distribution in the WD is similarly undisturbed. central ion number density higher than the central num- However, as the self-gravitating collapse proceeds, the ber density of the pure canonical CO WD (indeed, the number fractions of the CHAMPs and C/O ions eventu- extremal core [see below] central number density scales ally become comparable at the center of the WD. At this as nX mX ). See Fig.3 for an example number density point, the electron number density and radial distribu- profile∝ for such a cored WD, obtained by numerical solution begin to respond, increasing the pressure support to tion of the Tolman–Oppenheimer–Volkoff (TOV) equa- the X+ and, provided the total X+ CHAMP contamination [75, 76] assuming as input to the computation that tion MX+ is sufficiently small (to be quantified below), complete stratification has occurred (see AppendixD). re-stabilizing the X+ configuration at a much smaller ra- Note that the unphysical discontinuity in the ion num- dius. We expect that the endpoint of the self-gravitating ber density at the core boundary shown in Fig.3 arises collapse in this scenario is the formation of a core at the due to the continuity of the electron number density ne + center of the WD comprised of the two-fluid X /e− mix- (and hence degeneracy pressure) at the interface, com- 12 16 ture, with an overburden of the multi-fluid C/ O/e− bined with the abrupt change of charge per ion from +6 mixture of the canonical CO WD, with a transition re- for r > RX to +1 for r < RX , which allows 6 times more gion between these two stratified layers whose thickness ions to be present per unit volume just below the inter- is dictated by thermal effects.32 face, compared to just above; this discontinuity would We can roughly estimate the thickness h of the tran- be smoothed by the thermal thickness of the transition sition region by setting the thermal kinetic energy of a region, which we have have approximated to be zero con- CHAMP EK (3/2)T equal to the gravitational poten- sistent with our estimates above. ∼ + tial energy gained upon rising a distance h Rcore above While less extremal X /e− cores are ‘squashed’ to 2 2 the core boundary ∆E hmX Mcore/(MPlRcore), leading some extent by the overburden of CO WD material ∼ 2 to the estimate h/Rcore (3/2)TRcoreMPl./ (mX Mcore). (see, e.g., Figs.3 and9), as the central core becomes ∼ 5 For instance, with mX 10 GeV, an Mcore 4 more massive, it more and more closely approximates 10 ∼ ∼ × 10− M core (about 70% of the Chandrasekhar mass the physical structure of an isolated WD comprised of 8 + for this CHAMP mass) with Rcore 4.7 10− R a pure X /e− mixture; see AppendixD for a detailed ∼ × ∼ 30 m (see, e.g., Fig.3) inside a Mwd M WD gives discussion. In particular, the stratified core has a max- ∼ h 30 cm Rcore for T 10 keV. If the core is near- imum stable Chandrasekhar mass [55], which is equal ∼ 2 ∼ + extremal, Mcore mX− , while Rcore decreases roughly as to the Chandrasekhar mass for a isolated, pure X /e− 4/3 ∝ X+ 2 m− (see AppendixD), so it will always be the case WD: MChand.,X+ 5.7(µe )− M M , where X X+ ≈ that h Rcore for near-extremal cores. To give just one µ = m + /µ 1; see again AppendixD for a de- e X a more explicitly computed set of values to verify this scal- tailed technical argument as to why the Chandrasekhar 13 26 ing, for mX 10 GeV, an Mcore 3.7 10− M core mass scales in this fashion for µ 2. Therefore, if ∼ ∼ × e (about 70% of the Chandrasekhar mass for this CHAMP MX+ > MChand.,X+ , the central core cannot be suffi- 16 7 mass) with Rcore 4.6 10− R 3 10− m gives ciently pressure supported and, in the most conservative 9 ∼ × ∼ × h 3 10− m Rcore for T 10 keV. Since the tran- possible picture, will inevitably collapse to form a BH at ∼ × ∼ sition layer is thin, we will ignore it and approximate the the core of the WD. transition as immediate. Note that the picture of near-complete stratification depends on (a) the very different charge-to-mass ratios 2. Negatively charged CHAMPs: X− [105] of the CO mixture and the X+ (see Sec.VIC), and (b) the fact that the X+ will ‘use up’ (for lack of better For sufficiently small total X− CHAMP contamina- terminology) the ability of the electrons in the center of tion, MX− (again to be quantified below), the X− case the WD to supply pressure support, making them un- is na¨ıvely broadly similar to the X+ case, but a detailed available to supply pressure support to the CO mixture, consideration of the putative stable structure that would form at the center of the WD in this case leads to another possible avenue for triggering WD runaway upon BH formation. 32 Note in the alternative that if the thermal-pressure-supported structure does not form for any reason (see, e.g., footnote 28), As we have already discussed, X− in the Universe take then the CHAMPs sinking out of the full WD will simply ac- the form of (NX), where inside a WD we can assume cumulate in the core, and give rise to this stratified structure that N is either 12C or 16O, in roughly equal fractions; directly, once sufficiently many CHAMPs are present. see discussion in Sec. VI A 2. For the first part of this sub- 24

100 7.162 10 3 × − 1 10−

3 (0) 7.161 10− 5 4 2 × 10− 10− ion 10− /n )

r 0 ( 3 MWD = 1.0 10 M 10− × ion M = 4.0 10 10M

n X − × 3 4 RWD = 8.2 10− R 10− × 32 3 n (0) = 2.0 10 cm− ion × 5 10− 8 7 6 5 4 3 2 1 0 10− 10− 10− 10− 10− 10− 10− 10− 10 r/RWD

32 3 FIG. 3. The number density of ions as a function of radius, normalized to the core number density (nion 2.0 10 cm− ), + 10 5 ≈ × of an M = M pure-C WD with fractional X contamination η = 4 10− , for an mX = 10 GeV CHAMP (red, solid curve). We assume as input that the X+ have all settled to the core of× the star (complete stratification), forming an inner + 6 core of pure-X with mass MX+ 0.8MChand.,X , which is found to have radius RX 5.4 10− RWD (indicated by the ≈ 3 ≈ × vertical dotted black line); the whole star is found to have radius RWD = 8.2 10− R . We ignore thermal effects, which × would presumably mildly smooth out the discontinuity in ion number density at r = RX (see text). Also shown are (a) the ion number density profile for an equal-mass pure-C WD without X+ contamination (green, short-dashed line), which is found to + have an almost identical radius to the contaminated star, and (b) the ion number density for a pure X /e− WD of mass equal to the core mass (blue, long-dashed line). Note that, for ease of comparison, the latter two profiles are shown here normalized to the central density and radius of the X+-contaminated star, and not to their own central number densities. A comparison + of these results indicates that, because MX+ Mwd, the X have very little impact on the mechanical structure of the star outside the core (see inset). Moreover, because this particular core is not quite extremal, it does not well approximate the + profile of the isolated, equal-mass, pure X /e− WD, as we expect a more extremal core would; see discussion in AppendixD.

section, we make the assumption that the X− remain im- mass for the homogeneous (NX) core at the center of mutably bound in such structures even when a CHAMP the WD is 62 = 36 times larger than that for an X+ core forms inside a WD. In the second part of this subsec- CHAMP of≈ the same mass; essentially this can be un- tion, we will comment on the ions densities reached by derstood by virtue of that fact that, on average, each such a configuration and whether modifications to this electron in the core contributes pressure support for only picture are required. 1/6 of the mass of the CHAMP, as compared to the X+ a. X− bound as (NX) Since the (CX) and (OX) case, where the average electron supplies pressure sup- bound states have charge-to-mass ratios that differ by port for the full CHAMP mass. only 40%, it is not guaranteed that they will stratify The mechanical structure of the (NX) core is how- within the lifetime of the WD in the same way as the ever largely similar to that shown in Fig.3, except that CHAMP core stratifies from the CO overburden. As the discontinuity in ion number density across the core such, we will assume that the (CX) and (OX) bound boundary is reduced, since the average charge-per-ion in states form an approximately homogeneous mixture in the CO mixture if 7, whereas that in the (NX) mixture the core of the WD. This means that for the (NX) is 6. case, we have the mean molecular mass per free electron We conclude that if M , the total mass of X− in [Eq. (D10)] X− the WD core exceeds the Chandrasekhar mass for the X− 2 1 5X 7X 6µ 6 (NX)/e− fluid: MChand.,X− 5.7(µe )− M M , C O a ≈ = + + , (66) X− X− X where µ m /(Q µ ) 1 is the mean molec- µe mX /µa mX /µa ≈ mX ≈ µe e ≈ X− (NX) a ular weight per free electron for the (NX)/e− fluid and where we assumed that the WD has composition XC Q(NX) = 6 is the ‘mean-ion charge’, a central core of X 0.5. This is approximately equivalent to assum-≈ (NX) cannot be sufficiently pressure supported and, in O ≈ ing that the X− is bound to the ‘mean ion’ with charge the most conservative possible picture, will inevitably Q¯ 7, implying Q = 6. As such, the Chandrasekhar collapse to form a black hole at the core of the WD; ∼ (NX) 25

see also the discussion in AppendixD. In our modified case, we have a similar concern: once We are also finally able to see why any limits we the central ion density becomes too large, other processes ultimately set would be conservative (see comment in may begin to occur that would disrupt the structure we Sec. VI A 2) even if we were to assume that the X− that have computed here assuming hydrostatic equilibrium enter the star bound to He were not to exchange onto and a fixed equation of state. Similar to the ordinary a heavier nucleus: the Chandrasekhar mass for a (HeX) CO WD, pycnonuclear fusion of carbon ions brought into 2 2 2 core is a factor of Q(HeX)/Q(NX) = 6 = 36 smaller close proximity by the increased density caused by the than the Chandrasekhar mass for an (NX) core (assum- CHAMPs is the process most likely to disrupt the pic- ing Q(NX) = 6 in the mean-ion approach), because (by ture we have already outlined. Obtaining an accurate the same argument as we have just made above for the rate for this process is challenging, as it occurs as a tun- X+ state) the (HeX) state receives less pressure sup- neling process between carbon nuclei, each individually port than the (NX) state while retaining approximately bound to a CHAMP, within a larger ion lattice in which the same mass (in the limit mX µa), all other things other nuclear species (e.g., O) are also present; more- being equal. As such, if we take as our limit criterion over, electron screening enhancements of the fusion rate that a Chandrasekhar mass of CHAMPs are present (see are also present. In AppendixB we develop a very sim- Sec.VIII for further discussion), then we are clearly con- ple estimate of when this process may become impor- servative in assuming that the CHAMPs are all on the tant, finding [c.f., Eq. (B25)] that the density at which q high-charge nuclei. On the other hand, if33 the self- it becomes relevant is roughly nboom 0.2 n nboom, [CX] ∼ nucl [C] gravitating collapse is required to access this dense core 37 3 where nnucl 2 10 cm− is a nuclear number den- structure per the discussion in Sec.VID, and the self- boom ∼ × 32 3 sity, n 5 10 cm− is the number density at gravitating mass is larger than the Chandrasekhar mass [C] ∼ × which pycnonuclear processes become important in ordi- [see discussion around Eq. (68) below], then the charge boom 34 3 nary CO WD material [74], and n 10 cm− is of the ion in which the X− is bound is irrelevant anyway. [CX] ∼ b. High-sub-Chandrasekhar/trans-Chandrasekhar our estimate of when they become relevant in the (NX) mass of X−: the density of ions and pycnonuclear core of the WD. processes The results of AppendixD make clear that We thus see, in the trans-Chandrasekhar mass regime, the central ion density of the fully stratified core of the increasingly dense core may first undergo fusion pro- (NX) material at the central of the N star always cesses that would trigger thermal runaway, destroying the exceeds both the central ion density of an equal-mass WD before the core could even collapse to a BH. How- isolated pure (NX) WD, and the central density of the ever, as there is significant uncertainty in the estimate of uncontaminated pure N star (see Fig.8); in particular, the density at which these fusion processes first become the central density of the core can increase to be some relevant, we will in the following sections also proceed orders of magnitude above the central density of the conservatively to discuss in detail the case in which the CO star as the core becomes more and more extremal: core will proceed to collapse to a BH before it can be max nX mX . This however raises the possibility that, disrupted by fusion events (we must in any event discuss as the∝ core mass becomes trans-Chandrasekhar, it is that case, as these pycnonuclear processes are inapplica- no longer self-consistent to compute the structure of ble to the case of a X+ core). In this way, we will show the contaminated WD assuming the core is comprised that the inevitable outcome (with some minor exceptions immutably of (NX) material. in some regions of parameter space) is destruction of the Already in an ordinary CO WD, it is the case that WD, whether or not the pycnonuclear fusion process dis- three distinct outcomes can occur as the central den- cussed here can trigger the runaway. sity of the WD material increases [74]: (1) carbon- Note also that other processes would of course occur at carbon fusion reactions can become so efficient as the higher central core densities still but, again, these would core approaches the Chandrasekhar mass that the super- still lead to early destruction of the WD: for instance, by nova instability is triggered before the Chandrasekhar the time the ion density reaches nuclear density, fusion mass is actually reached, (2) the core can exceed the would be almost entirely unsuppressed, but there is ad- Chandrasekhar mass and collapse (likely triggering the ditionally the consideration that, at such large densities, supernova instability as it collapses), or (3) the core the X− would essentially no longer be individually asso- can begin to neutronize, which modifies the equation of ciated with single positively charged nuclei, as the radii state. In an ordinary CO WD, outcome (1) is believed of the (NX) bound states are of order the nuclear radius. 32 3 to prevail at densities around n 5 10 cm− , al- C ∼ × Instead, the nuclei N may become delocalized in an X− though the central densities at which (2) and (3) occur— lattice, akin to a metal structure. In this case, if any X+ 33 3 33 3 n 1.3 10 cm− and n 2 10 cm− [74], C ∼ × C ∼ × are present, they would no longer be automatically pre- respectively—are both within a factor of 3–4 of the cen- vented by a large Coulomb barrier from annihilating with tral density at which (1) occurs [74]. an X−, which could lead to an additional large energy deposition in the core, again increasing the likelihood of triggering the supernova runaway. 33 See footnote 28. Finally, we comment that we have not included any 26

prim. Coulomb corrections to the EoS of the material in the M > MChand., X [implying t < τsink]; and (2) X± ± WD, or in the CHAMP-contaminated core; at the ex- s.g. prim. ∗ prim. mX < m , and M < MChand., X , and M + X X± ± X± tremely high charge densities implied by nX mX , this accum. M (τaccum.) > MChand., X± [implying that τsink . may be a poor approximation, which adds additional∝ un- X± t . τaccum.]. certainty to any rate estimate we could give for any of ∗ It additionally holds for ts.g. . τstrat. < t < t where these processes (it could not however modify the ultimate t is the time for the onset of self-gravitating∗ col- existence of a Chandrasekhar limit). s.g. lapse, τstrat. is the time for the formation of a stratified core during the self-gravitating collapse, and t is ∗ as above, under the following circumstances: (3) mX > 3. Comment on applicability s.g. prim. mX , and MX > Ms.g., X > MChand., X± [imply- ± core s.g. ing t < ts.g. + τcoll. . τsink]; and (4) mX > mX , We can now finally quantify our earlier caveats in ∗ prim. and M < Ms.g., X [implying ts.g. > τsink], and Secs. VI E 1 and VI E 2 that the discussion of the equi- X± prim. accum. librium structure in those preceding sections is applica- MX +MX (τaccum.) > Ms.g., X > MChand., X± [im- ± ± core ble ‘[f]or sufficiently small total X± CHAMP contamina- plying t < ts.g. + τcoll. . τaccum.]. ∗ tion’. We have now discovered that this condition is M X± ≤ MChand., X± , where the Chandrasekhar mass here is the + 4. Combined X and X− case one applicable for the X±, and MX± is the total mass of X± in the stratified core at the time t at which the structure is considered. In particular, this discussion definitely prim. accum. Absent a near-complete charge asymmetry between holds if M + M (τaccum.) < MChand., X . + X± X± ± X and X (which could occur in some production To understand where else the discussion might hold, − s.g. mechanisms, or if all but a small residual charge asym- consider that there is a CHAMP mass mX such that for metric population of CHAMPs annihilates away), it is ρwd highly likely that the physical case will correspond to a m ms.g. 0.85π Q4 (67) X & X ≈ T 3 X WD that is contaminated by some abundance of both 3 + T − X , and X− bound as (NX). As we have discussed in 6.4 108 GeV Q4 Sec.IV, the large Coulomb barrier between X+ and the ∼ × × X × 1 keV positively charged (NX) state will prevent the X+ and ρwd (68) X− annihilating, so both species can co-exist in the core. × 5.5 107 g/cm3 × This leads to two possible outcomes: (1) the core is a roughly homogenous mixture of X+ and (NX), or (2) the (with the numerical coefficient found using the re- + + sults of AppendixD), the mass required for the ini- X and (NX) stratify, with the X sinking to the core tial thermal-pressure-supported structure34 to undergo (smaller charge to mass ratio; see Sec.VIC). To estimate self-gravitating collapse exceeds the Chandrasekhar mass which occurs, let us return to a timescale estimate for the s.g. sinking of X+ through a (NX) fluid similar to that we for the same mX : i.e., Ms.g.,X (mX & mX ) M (m ms.g.). However, the discussion of≥ made in Sec.VIC. Assuming nX+ n(NX) 0.5 as an Chand.,X± X & X ¯ ∼ ∼ the core structure in this subsection can still hold in- approximation, Q(NX) = QN QX with QX = 1 and ¯ − | | | | stantaneously for this case, as the timescale for the QN 7 in the mean-ion approach, and mX GeV, the ∼ thermal-pressure-supported structure to collapse is rea- sinking terminal velocity for the X+ in the core is given sonably long [Eq. (65)] compared to the sound-crossing by [100, 107] time of the core (which sets the dynamical timescale over which significant structural alterations can be compen- core g(r)DmX QX 5 g(r)DmX sated for), and a sub-Chandrasekhar core will still form wX+ (r) 2| | 1 , ∼ T Q¯N − ∼ −7 T at some point during this collapse, before the accretion of (69) the additional mass from the thermal-pressure-supported structure onto the core drives the core mass over the 2 Chandrasekhar limit; see Sec.VII for further discussion where g(r) = Menc(r)/(rMPl.) with Menc(r) 3 ≈ of the subsequent dynamics. Mcore(r/Rcore) approximated as for a uniform density Therefore, the discussion in Secs. VI E 1 and VI E 2 also spherical core, and D approximated by the self-diffusion holds instantaneously for times t < t where t is the time coefficient Eq. (43). Note however that in evaluating D ∗ ∗ we must now take m m and Z 0.5Q + at which the mass of X± in the stratified core equals ion ∼ X ion ∼ (NX) deg. 0.5 QX 3.5 as appropriate for an initially homoge- the Chandrasekhar mass, M (t ) = MChand., X , un- X± ± | | ∼ + der the following circumstances:∗ (1) m < ms.g., and neous mixture in the core of X and (NX) with equal X X number densities through which the X+ sink. Assuming Mcore and Rcore to be the mass and radius of a near- extremal ‘mini-WD’ degenerate core [see Eq. (D23) and 34 If it forms; see footnote 28. (D24)], we find that for T keV we have a sinking ∼ 27

sink s.g. timescale of [rX (t) = rX (0) exp( t/τX+ )] star if mX & mX , so that self-gravitating collapse of − the CHAMPs can occur and give rise to the stratified 17/9 1/3 sink mX − T − structure (see Secs.VID and VI E 3). In either case, this τ + 3 yr . (70) X ∼ × 105 GeV × keV abundance of CHAMPs could either be present primor- This estimate of the stratification timescale indicates dially, or could accrete onto the WD before the crystal- that an initially mixed core should stratify almost im- lization time. mediately over most of our parameter space; in particu- The timescale for formation of the BH during the col- lar, this estimate is faster than the thermal core collapse lapse of a trans-Chandrasekhar extremal core initially supported by degeneracy pressure is extremely short, of time Eq. (65) [relevant if Ms.g.,X > MChand.] or the sinking time Eq. (56). The equilibrium structure in this case the order of the gravitational free-fall timescale (see, e.g., + the discussion in 36.3.1 of Ref. [66]): is of course more complicated than either of the X or § (NX) cores alone, consisting now of three nested spher- s M 2 ical shells. However, since the Chandrasekhar mass for τ Pl. (71) X+ is a factor of 36 lower than for (NX), it will a rea- collapse. ∼ M /R3 ∼ Chand. Chand. sonable approximation to estimate the maximum mass of 5 + + 10 GeV the X core to be that of the isolated X core in most 1.5 µs QX (72) cases unless the total (NX) mass exceeds 36 times the ∼ × × mX X+ mass. Likewise, when the total X+ ∼mass is much where QX is the charge of the object in which X appears below 1/36 times the total (NX) mass, it is likely a + in the WD [Q = Q + for X ; Q = Q for X−]. ∼ + X X X (NX) good approximation to ignore the central X core and Note that this collapse provides another possible av- estimate the maximum allowed (NX) mass as one would enue to trigger thermal runaway: the gravitational bind- for the isolated (NX) core. In the intermediate regime ing energy released by the collapse can heat the WD. where the total mass of (NX) is roughly equal to 36 Suppose we only track the collapse as far as when the + times the mass of the X , a more complicated situation core has collapsed to half its initial radius, at which point will arise, but this is a tuned region of parameter space, energy of order of its initial binding energy has been re- and the allowed total mass will be at most an (1) factor O leased by the collapsing CHAMPs. The energy can be different from the maximum allowed (NX) mass in this estimated to be case. 2 init. MChand. EB 2 (73) ∼ MPl.RChand. F. Super-Chandrasekhar mass of CHAMPs: BH 5 8/3 formation, and timescales 41 8/3 10 GeV 8 10 GeV QX (74) ∼ × × × mX As noted in Sec. VI E 3, our discussion thus far has If a fraction of this energy could somehow converted to been largely limited to the impact of a total mass heating of the WD, thermal runaway may possibly be of CHAMPs in the stratified core of the WD, M core, X± triggered. Assuming that any such energy released would that is sub-Chandrasekhar at all times (with the excep- be deposited within a volume parametrically of size the tion of our comments in Sec. VI E 2 b, which considered trigger volume, a na¨ıve criterion for the runaway to be the case of M core in the high-sub-Chandrasekhar/trans- X± initiated is that Chandrasekhar mass range; and our comments in init. Sec. VI E 3). However, if a super-Chandrasekhar mass EB of CHAMPs is able to collect in the stratified core, the Edeposited ζ min [τdiff., τcollapse.] & ET , (75) ∼ τcollapse. quiescent structure discussed in Sec.VIE will not be the equilibrium structure (or even the instantaneous quasi- where ζ 1 controls the fraction of available energy that ≤ static structure). Instead, provided either that we con- is deposited in the WD as heating. For a Mwd 1.1M + 10 ∼ sider X , or the considerations of Sec. VI E 2 b are inap- WD, we find that for mX & 5 10 GeV QX , we × 7 × 8/3 plicable for X− in the high-sub-Chandrasekhar regime, have τ τ , so that ζ 4 10− Q− diff. & collapse & × × X × the general relativistic instability to collapse of a suffi- 10 8/3 mX /5 10 GeV is required, which is possible with ciently dense central core structure implies that once the × 13 sinking CHAMPs accumulate a Chandrasekhar mass ζ . 1 for mX . 10 GeV QX . For mX . 5 ∼ 1010 GeV Q , we have τ ×τ , and we find that× in the stratified core, that core will collapse to a BH × X diff. . collapse 7 5/3 10 5/3 with an initial mass equal to the relevant Chandrasekhar we need ζ 4 10− Q− m /5 10 GeV < & × × X × X × mass. Note however that in order to obtain such a mas- 7 5/3 4 10− QX− . We therefore see that, for mX . sive stratified core, a larger total mass of CHAMPs than 10×13 GeV × Q , there is a possibility of that some frac- 35 X the Chandrasekhar mass may be required to be in the tion ζ ×1 (and possibly ζ 1) of this initial gravitational≤ binding energy release being deposited in the WD and leading to heating would be sufficient to trig- 35 See again footnote 28. ger thermal runaway. Of course, our estimates here are 28

schematic and approximate, although they are conserva- mechanism for SNIa-like supernova events in WD due to tive in the sense that further energy is released as the the formation of a BH inside the WD (or, more conser- core continues to collapse. All we intend to argue here vatively, for destruction of the WD by accretion onto the is that there is yet another plausible alternative WD de- BH), while our discussion of WD outcomes draws also on struction mechanism that could be operative for some a longer series of prior studies in this space [57–61, 63, 65]. range of CHAMP masses, even before the BH is reached; We note however that all numerical or quantitative esti- see also Refs. [60, 61] for similar arguments in a diﬀer- mates here are performed independently of past work to ent context. Of course, if no such mechanism exists to the extent this is possible. cause a concomitant heating of the WD, the core simply Finally, we note that some of the accretion rate esti- collapses to a BH without possibly triggering thermal mates we develop in this section are quite approximate runaway. models for the complex accretion dynamics; where rele- Once the BH forms at the Chandrasekhar mass, its vant, we note where our physical conclusions depend on subsequent dynamics and ultimate fate depend on a num- details of uncertain estimates. Resolution of these uncer- ber of diﬀerent physical processes. Moreover, the ulti- tainties would likely require numerical modeling. mate fate of the WD in which the BH forms also depends on which of the various BH dynamical mechanisms dominate at various points in the BH evolution. In Sec.VII, A. Dynamical processes governing BH evolution in we undertake a detailed consideration of these points. a WD

There are three possible dynamical processes that im- VII. BLACK HOLES IN WHITE DWARFS pact the evolution of the BH mass: (1) Hawking radiation by the BH, assuming that this process occurs in the usual In Sec.VIF, we noted that if a sufficient mass of fashion it is believed to occur (see, e.g., Ref. [111] for one CHAMPs is present primordially or can accrete onto the alternative viewpoint); (2) Bondi (or Eddington) accre- WD within the WD crystallization time τaccum., a BH is tion of CO WD material from the core of the WD; and born at the center of the WD. Here, a ‘sufficient mass’ (3) accretion of CHAMPs from the WD, which occurs at means the larger of the Chandrasekhar mass and the self- various rates throughout the BH evolution, depending on gravitating mass: gravothermal collapse is necessary to a variety of criteria. In this subsection, we consider each form a degenerate core that can later be pushed over the contribution in turn. Chandrasekhar limit in order to trigger the gravitational collapse instability required to give birth to a BH.36 This holds true even in the case where the self-gravitating 1. Hawking radiation from the BH mass exceeds the Chandrasekhar mass, because the thermal structure discussed in Sec.VID takes some time to Once the BH forms, it is widely believed that there fully collapse, and a degenerate core must be born dur- will be a dynamical contribution to its evolution that ing that collapse, even if it is short lived. We will assume will cause it to lose mass: the Hawking process [56]. throughout that the BH evolution begins with a BH born This process causes the BH event horizon (which for a at the Chandrasekhar mass at time tbh, but that various non-rotating BH lies at the Schwarzschild radius, RS = 2 phases of evolution of the BH dynamics can follow. In 2Mbh/MPl.) to emit a thermal spectrum of particles at a 2 this section, we turn to an examination of the dynamics temperature given by TH = MPl./(8πMbh). The Hawk- of the BH, and the implications for old WD. ing mass loss rate is given approximately by We begin in Sec.VIIA by examining in detail the var- 2 ˙ π geff.(TH ) 4 2 ious dynamical processes that govern the BH mass evo- Mbh, H TH (4πRS) (76) lution after BH formation, under various assumptions. ∼ − 60 2 4 Thus armed with an understanding of the relevant dy- (geff.(TH )/2)MPl. = 2 , (77) namical processes, we turn in Sec.VIIB to a consid- − 15360πMbh eration of the temporal evolution of the BH mass under the action of these dynamical contributions, mapping out where geff.(T ) is the effective number of relativistic de- how the BH mass evolves qualitatively (and, where nec- grees of freedom at temperature T . essary, quantitatively) in various regions of parameter If the Hawking radiation completely dominates the BH space. Finally, in Sec.VIIC, we turn to a consideration evolution, the BH mass will decrease, and the timescale of the impact of the BH evolution on the fate of old WD. to radiate away the entirety of the BH mass is finite; it is Our discussion of the BH dynamics and WD outcomes dominated by the time spent at the largest masses [61]. Making the conservative assumption that geff(TH ) re- is guided by Refs. [60, 61], which considered the trigger 37 mains constant at some initial value geff, 0, the timescale

36 Although see again the caveat at footnote 28. 37 This is conservative because as the BH mass drops, the Hawking 29 can be estimated by It is also worth noting the mass at which the Bondi and Hawking rates balance (ignoring any CHAMP con- 15360π 3 τbh, H 4 Mbh, 0, (78) tribution; see below): ∼ 3(geff, 0/2)MPl. 8 3 1/4 where Mbh, 0 is the initial mass. Numerically, we find B H (geff., 0/2)MPl.cs Mbh∼ 2 (83) 3 ∼ 61440π λρwd Mbh, 0 τ 3.4 Gyr , (79) 3/4 bh, H 19 c ∼ × 10− M 38 s = 1.5 10 GeV 2 × × 2.8 10− where we took geff, 0 51/4, as appropriate for tempera- × ≈ 1/4 1 tures TH . MeV. geff., 0 λ − × 51/4 × 1 ! 2 2. Accretion of WD material onto the BH ρ − wd (84) × 5.5 107 g/cm3 The BH will also accrete CO matter from the surround- × ing WD. We assume that the accretion rate of WD matter where we assumed fiducial values for a Mwd 1.1M on the BH is given (at least initially; see discussion be- WD (see Tab.I). Assuming those values, we∼ have B H 19 low) approximately by the Bondi accretion rate [74, 112]: Mbh∼ 1.3 10− M , which is interestingly and coin- cidentally∼ × 2 a region where the intrinsic growth and evapo- ˙ Mbh ration timescales for the BH happen to be about as long Mbh, B +4πλ 2 2 ρwdcs, (80) ∼ MPl.cs as current old WD ages, few Gyr. a. Eddington-limited∼ accretion× However, accretion where λ (1), c is the WD sound speed, and ρ is s wd can become (Eddington) limited if radiative backreac- the WD∼ density; O the latter two quantities are as evalu- tion on the in-falling matter in the region near the ated at the center of the unperturbed WD (i.e., far from sound horizon supplies sufficient pressure support to the BH event horizon, but still near the center of the re-establish hydrostatic equilibrium conditions near the WD). Parametrically, this estimate follows, up to (1) sound horizon [50, 61, 66, 74, 113]. The Eddington lu- numerical factors, by assuming that particles crossingO the minosity is given by L = 4πM /(κ M 2 ), where sonic radius R R /c2 GM /c2 at which the es- edd bh rad Pl. sonic s s bh s κ = σ/m is the radiative opacity with σ the cape speed from the∼ BH is∼ equal to the sound speed in rad ion photon scattering cross section with the stellar mat- the (unperturbed) WD will continue to accrete onto the ter. Following Ref. [61], we will assume that σ = BH. max 1 mb, 100 mb (ω/MeV) 1 σ 200 mb for If Bondi accretion completely dominates the BH evolu- − & ω 0.5 MeV ( T × ). ∗ ≡ tion, the BH grows, and the timescale for it to accrete up . crit. Assuming that∼ a fraction of the in-falling matter is to the full WD mass (assuming nothing cuts this evolu- converted into outgoing radiation, the Eddington-limited tion off) is again finite; it is dominated by the time spent mass accretion rate is at the smallest masses [61]. The timescale for the BH to accrete from an initial mass Mbh, 0 to the full mass of the M˙ 4πM m /(σM 2 ). (85) WD, assuming the Bondi accretion rate holds through- bh, E ∼ bh ion Pl. out, is given approximately by Equating this to the Bondi accretion rate, we estimate 3 4 that accretion becomes Eddington limited for csMPl. 1 1 τbh, B Mbh− , 0 Mwd− , (81) ∼ 4πλρwd − 3 2 B E csmionMPl. M ∼ (86) where λ is the same (1) constant as in Eq. (80), bh λσρ O ∼ wd Mbh, 0 M is the initial BH mass, and Mwd M 3 1 ∼ cs is the WD mass. Taking parameters appropriate for 3.5 1042 GeV − ∼ × × 2.8 10 2 × 0.1 the center of a Mwd 1.1M WD (and treating ρwd − ∼ × as constant, which is a reasonable approximation as the 1 1 λ − σ − timescale is dominated by the time spent at the smallest × 1 × 200 mb masses for the BH, when the WD is mostly unperturbed), ! 2 we find numerically that ρ − wd , (87) 1 × 5.5 107 g/cm3 Mbh, 0 − × τbh, B 3.1 Gyr 18 . (82) ∼ × 10− M where we assumed fiducial values for a Mwd 1.1M WD (see Tab.I). Note that, assuming the fiducial∼ values, B E 15 Mbh∼ 3 10− M . ∼ × B E temperature increases, enabling the relativistic participation of For Mbh . Mbh∼ , the accretion of WD material is more SM species, increasing geff. and the rate of emission. B E at the Bondi rate; for Mbh & Mbh∼ the accretion is 30

Eddington limited. We take assume that this pre-BH evolutionary phase completes h i before the BH could form, so there is no diffuse over- M˙ bh, WD material = min M˙ bh, B, M˙ bh, E (88) dense CHAMP structure at the center of the WD to feed the BH accretion at the time of BH formation. Note however that the Eddington-limited rates will be Consider first the CHAMPs present in the WD pri- essentially irrelevant for most of our considerations, as- mordially. The rate at which these CHAMPs accrete suming fiducial values for a Mwd 1.1M WD. Even onto the newly formed BH is somewhat uncertain, and ignoring any (necessarily positive) CHAMP∼ contribution we give two distinct estimates for it. The first estimate is na¨ıve: we simply assume that for to M˙ bh, the mass accretion rate onto the BH assuming Bondi accretion of WD material already exceeds the t . τsink, the primordial CHAMPs accrete onto the BH Hawking rate in magnitude at a BH mass three orders of at a constant rate equal to the average accretion rate of magnitude smaller than the mass at which the Bondi rate all the primordial CHAMPs: becomes Eddington limited. The impact of the Edding- h i(1) ˙ 1 prim. ton rate limitation is thus only felt deep into the regime Mbh(t) τsink− MX Θ(τsink t), (92) prim., X ≈ ± − where the BH is rapidly accreting matter from the WD: the timescale to accrete to the full WD mass from M B E prim. bh∼ where τsink is given approximately by Eq. (91), M is X± is 3 Myr assuming Bondi accretion, while the cognate given by Eq. (27), and estimate∼ for the Eddington-limited accretion rate is ( 2 1 x > 0 σMPl. Mwd Θ(x) = (93) τbh, E ln B Edd. (89) 0 x 0 ∼ 4πmion Mbh∼ ≤ σ prim. 33 Myr is the Heaviside theta function. Note that if M < ∼ × 0.1 × 200 mb X± ! MChand., X , then this rate is zero since τsink < tbh, where ρ t is the BH formation time. wd . (90) bh × 5.5 107 g/cm3 The second estimate is more refined: assuming an ini- × tial uniform distribution of the CHAMPs, in the time While the Eddington accretion timescale here is (10) interval t [t?, t? + dt?], those CHAMPs present in the O radial slice∈r [r , r + dr ] will have sunk to the core of times longer than the Bondi timescale, the absolute ∈ ? ? ? timescale is still extremely short, τbh, E Gyr τwd, the WD, and will be accreted onto BH. The values of t? ∼ 38 so this small additional amount of evolution time can be and r? = r?(t?) are given by neglected completely in what follows. These conclusions Z r? dr0 agree with those of Ref. [61]. t? = , (94) Rf wX (r0) where, to make the estimate tractable, we take w to be 3. CHAMP accretion onto the BH. Case I: X given by the analytical estimate for the sinking speed of Ms.g., < MChand., . X X the CHAMPs at Eq. (49). It follows from the Leibniz rule that It is our assumption throughout this paper that the initial CHAMP distribution in the WD upon formation is dr? dt? = . (95) homogeneous (see later comments in Sec.VIIID). These wX (r?) CHAMPs will however begin to sink toward the center of The mass accretion rate is thus given by the WD after WD formation, as discussed in Sec.VIC. 2 Our best estimate for the timescale for the sinking of dM(t?) 4πr?ρX dr? 2 CHAMPs from the outskirts of the WD to the core is = = 4πr?ρX wX (r?), (96) dt dr?/wX (r?) given by Eq. (56), which we will approximate here as where ρX is the average (assumed uniform) mass den- 5 6 10 GeV sity of the primordial CHAMPs in the star: ρX = τsink 4 10 yr , (91) 3 ∼ × × mX 3MX /(4πRwd). Since the central object at time t? con- sists of all CHAMPs accreted up to radius r?, we have ignoring the small logarithmic correction. 3 3 In this subsection we will assume that we are in the 4π 3 r? r? Mcore M(t?) ρX r? MX 3 ηMwd 3 , regime where the self-gravitating mass of the initial ther- ∼ ∼ 3 ∼ R ∼ Rwd mal structure (if any) of CHAMPs formed at the core (97) of the WD Ms.g., X has a mass smaller than the Chan- drasekhar mass. If this is the case, the mass contained in the self-gravitating core collapses to a stratified core ob- 38 ject below the Chandrasekhar mass in the time Eq. (65), Note that r? as defined here should not be confused with r refined by Eq. (58). ∗ which is somewhat shorter than τsink or τaccum.. We thus 31 so that in Eq. (49) we can take of τsink for the CHAMPs first accreted onto the WD to " # sink to the core and a steady-state flow to be established; M M R 3 none of our results depend sensitively on this assumption. 1 core + core wd − Mwd Mwd r? M 1 core + η 1, (98) 4. CHAMP accretion onto the BH. Case II: ∼ − Mwd ∼ Ms.g., X & MChand., X . leading to If the self-gravitating mass Ms.g., X exceeds the Chan- 3 dM(t?) 4πr?ρX wX (Rwd) wX (Rwd) drasekhar mass MChand., X there is a modification to = 3 M(t?). (99) the picture advanced in Sec. VII A 3; although see foot- dt ∼ R Rwd note 28. In this case, when the BH is initially born Given that the central density of a WD is higher than at the Chandrasekhar mass, it is surrounded by a sig- the average density, our discussion in Sec.VIC leads nificant overdensity of CHAMPs in the collapsing self- us to conclude that this estimate is too small by, con- gravitating cloud; this can temporarily boost the accre- servatively, a factor of 6.5 [see, e.g., Eq. (55)] for tion rate of CHAMPs just as the BH is born, and until ∼ the central density and core masses considered in arriv- the self-gravitating cloud has been accreted in a time ing at Eq. (56); in the same approximation, we have core τcoll. [Eq. (65)]. We must therefore modify our estimates wX (Rwd)/Rwd 2/τsink where τsink is as given by for the CHAMP mass accretion rates Eqs. (92) or (100), ∼ Eq. (91). Therefore, we can take the second estimate or Eq. (101) [depending on which process is the prevail- of the accretion rate to be ing rate at the time of self-gravitating collapse, ts.g., and (2) 1 which estimate we use for it] during this initial time. M˙ (t) 40τ − M(t), (100) prim., X ∼ sink We provide a simple estimate for this initial accretion rate, assuming just that the self-gravitating core is ac- where τsink is as given by Eq. (91) and where the numerical coefficient is a reasonable approximation when creted at an average rate until it is completely depleted: 10 M(t) MChand.,X (mX 10 GeV) Mwd, with only ˙ core 1 core ∼ ∼ Ms.g., X (t) (τ )− Ms.g., X Θ(ts.g. + τ t), logarithmic dependence on this assumption so long as the coll. ± coll. ≈ − (102) approximations used in deriving this estimate remain satisfied. This accretion rate is of course cut off when the which holds during the time ts.g. . tbh < t < ts.g. + entire primordial abundance of CHAMPs has accreted core τcoll. . This is a reasonable estimate because the mass onto the central object. in the self-gravitating cloud is initially concentrated in Note that the second estimate [Eq. (100)] is much a Gaussian profile near the WD core, and not spread smaller than the first [Eq. (92)] when M(t) M prim. X diffusely throughout the star. [e.g., when M(t) MChand., X ], and also that it depends core sink ∼ Because 25τcoll. τcoll. [cf., Eqs. (65) and (91)], the explicitly on the mass of primordial CHAMPs already in accretion rate estimate∼ Eq. (102) is comparable to the the central object. We comment below on the differing estimate Eq. (100) if M(t) Ms.g., X± is assumed in implications of the two rate estimates. the latter rate. Therefore, assuming∼ that the full self- Provided that the WD is located where a significant gravitating mass were to collapse into the BH (see dis- galactic halo CHAMP abundance is present, there is a cussion in Sec.VIIB below), the matching between the second CHAMP accretion mechanism that operates for accretion rates Eqs. (102) and (100) would be roughly a much longer time period. As already discussed in smooth [up to (1) factors] at the end point of the self- Sec.VIB, halo CHAMPs accrete onto the WD (provided gravitating collapse.O Of course, for earlier times [e.g., that the WD has a sufficiently small magnetic field; see supposing that M(t) MChand., X± Ms.g., X± ], the discussion in Sec. VI B 2). If we assume that there is a rate at Eq. (102) exceeds∼ that at Eq. (100 ) [as expected], quasi-steady-state period of WD evolution during which by roughly the ratio Ms.g., X /MChand., X . Therefore, these CHAMPs simply sink through the star and are cap- ± ± if ts.g. < τsink, our refined approach will be modified by tured by the BH, then the mass accretion rate onto the assuming only the accretion rate Eq. (102) for the initial BH during this quasi-steady-state period is simply equal time period of the self-gravitating cloud collapse, then to the CHAMP mass accretion rate onto the WD: reverting to the rate Eq. (100) until t τ .39 ∼ sink 1 accum. M˙ accum., X (t) τ − M (τaccum.) ≈ accum. X± Θ(τ t)Θ(t τ ), (101) × accum. − − sink accum. 39 Roughly, the derivation of Eq. (100) still holds after the initial where M (τaccum.) is given by the expressions in X± phase where the CHAMPs are stalled in the thermal structure, Sec. VI B 1 and AppendixC, and τaccum. is the timescale provided we make the reasonable assumption that the dynamics for WD crystallization (see the discussion in Sec. VI B 1). of the CHAMPs sinking from radii large enough that they arrive Note that we have set this rate to zero for times t . in the vicinity of the core after the collapse of the self-gravitating cloud are not impacted by the earlier presence of that cloud. τsink because it will generally take a time on the order 32

However, Eq. (102) is smaller than the na¨ıve accretion on the WD sufficient to trigger the supernova instability prim. estimate Eq. (92) when M 25Ms.g., X . If the [60, 61], thereby terminating the BH evolution and giv- X± & ± na¨ıve approach is used, and ts.g. < τsink, we will instead ing a clear observable signal. However, in case (3), there assume that the accretion rate is the larger of Eq. (102) or will (except for highly tuned regions of parameter space) Eq. (92) during the initial self-gravitating cloud collapse, be no observable impact on the WD in which the BH and if necessary then revert to the rate Eq. (92) until resides. t τ . Because case (3) is an important exception, we discuss ∼ sink In both cases, if τsink < ts.g. < τaccum., we assume only the requirements to be in this region. For concreteness, the accretion rate Eq. (102) for the initial time period the qualitative discussion that follows is correct in de- of the self-gravitating cloud collapse, and then revert in- tail (and has been checked quantitatively) for CHAMPs + stead to the rate estimate Eq. (101) until t τaccum.. X present in a WD with mass Mwd 1.1M that has ∼ This piecewise approach to the accretion∼ rate esti- characteristics as in Tab.I and which has an assumed mates is manifestly approximate, but it roughly cap- crystallization time of τaccum. Gyr and an assumed ∼ tures the prevailing dynamics of the major epochs of age/lifetime of τwd 2 Gyr (see Sec.VIII and Tab.III). ∼ CHAMP accretion under the various assumptions about Similar considerations, albeit with some modifications as the accretion dynamics. The exact accretion behavior in to whether certain regions of behavior exist, are appli- the regimes where various CHAMP reservoirs are nearly cable for CHAMPs X− bound as (NX), or for different exhausted (e.g., when the self-gravitating structure is WD parameters. nearly fully accreted) will of course be more complicated It turns out that, for X+, the conditions under which and our model for the accretion rate would not fully cap- outcome (3) can occur are such that the total mass of ture all the subtleties of the dynamics during these times. CHAMPs present primordially in the WD exceeds the We expect that the gross picture of the BH dynamics self-gravitating mass, which in turn exceeds the Chan- is captured by the approach we have outlined here, al- drasekhar mass (as opposed to a situation in which the though details will be lost. primordial mass of CHAMPs is too small to exceed the larger of the Chandrasekhar or self-gravitating masses, but the slow accretion of CHAMPs onto the WD triggers B. Implications for BH evolution BH formation at a later time, but before the WD crystallization time). To maintain concreteness, we discuss only Of the dynamical contributions to the BH mass evolu- this ordering of the self-gravitating and Chandrasekhar tion discussed in the Sec.VIIA, the WD and CHAMP masses in what follows, assuming the existence of the accretion processes tend to increase the BH mass, while thermal-pressure-supported phase (see footnote 28). the Hawking process tends to decrease it. Moreover, the As the CHAMPs in the star collect in the central various accretion and emission processes have rates that thermal-pressure-supported structure, gravothermal col- are explicit functions of time, as well as functions of the lapse will eventually be triggered at some time ts.g. < BH mass in many cases. In general, this leads to a rich τsink. As the thermal cloud of CHAMPs collapses in a core and complicated dynamics of the BH, leading to a sur- time τcoll. thereafter, a stratified core will form, which prising complexity of possible BH evolutionary trajecto- will accrete CHAMPs from this initial thermal cloud ries. structure. Once sufficiently many primordial CHAMPs There are however only three distinct physical out- have accreted onto this core, a BH will form at time tbh core comes: (1) the BH will reach an evolutionary phase dur- where ts.g. . tbh . ts.g. + τcoll. , with a mass at the Chan- ing which it undergoes runaway accretion with a rate suf- drasekhar mass. ficient in principle to accrete the entire WD mass within At this point, the remainder of the thermal cloud is still the WD lifetime, the countervailing effect of Hawking ra- collapsing onto the newly formed BH, supplying a large diation notwithstanding; (2) the BH will reach a phase positive contribution to the BH mass accretion rate. Nev- during which it undergoes Hawking evaporation at a rate ertheless, it is possible that the Chandrasekhar mass is sufficient in principle to radiate away its entire mass so small that, including the effects of Hawking radiation ˙ within the WD lifetime, the countervailing accretion of from the newly formed BH, Mbh(Mbh = MChand.; tbh . core WD material and/or CHAMPs notwithstanding; or (3) t . ts.g. + τcoll. ) < 0, and the BH would immediately the BH will follow an evolutionary trajectory along which begin to lose mass upon formation; outside of a highly it can eventually be either dominantly accreting or dom- tuned region, the timescale for this process is so fast that inantly evaporating, but either owing to a balancing of the BH would evaporate well within the WD lifetime. rates of accretion and evaporation, or simply the intrin- This avoids outcome (3). ˙ core sic timescale for the dominant process being too long, the If Mbh(Mbh = MChand.; tbh . t . ts.g. +τcoll. ) 0 (in a timescale for completion of the accretion or evaporation highly tuned region), the BH will neither grow nor≈ shrink processes is longer than the observed age of the WD. In until such time as all the CHAMPs in the thermal cloud scenarios (1) or (2), we will find that the accretion or have been depleted, at which point the BH is still ap- evaporation likely does not actually need to proceed all proximately at the Chandrasekhar mass. However, once the way to completion in order to trigger a back-reaction the CHAMPs initially in the thermal cloud are all ac- 33

˙ core creted, Mbh(Mbh = MChand.; ts.g + τcoll. < t < τsink) τaccum.) > 0, two outcomes are again possible: in the first may (depending on parameters and on which of the case, the intrinsic timescale for accretion of CHAMPs two estimates discussed in Sec. VII A 4 we give for the and WD material is so fast that the BH will in princi- accretion rates) lose a positive contribution present in ple be able to accrete up to the full WD mass within ˙ core Mbh(Mbh = MChand.; tbh . t . ts.g. + τcoll. ) and turn the WD crystallization time; this avoids outcome (3). In ˙ prim. negative; this is particularly true if we use the refined the second case [when Mbh(Mbh = MX ; τsink < t < approach of Sec. VII A 4 because Ms.g., X > MChand.,X . τaccum.) 0, but slightly positive], the timescale for ac- The intrinsic timescales involved for the Hawking pro- cretion of∼ CHAMPs and WD material is too long to al- cess are sufficiently rapid that, outside a highly tuned low the BH to accrete up to the full WD mass before region, the BH will then radiate away in a fraction of the WD crystallization time, but the BH can still accrete the remaining lifetime of the WD, and outcome (3) is CHAMPs up to the full abundance that accumulate onto avoided. Alternatively, if we use the na¨ıve approach the WD before the crystallization time; such a BH can detailed Sec. VII A 4 and Eq. (92) is the accretion es- still end up with outcome (3). timate throughout the thermal cloud collapse and sink- Outside of tuned regions of parameter space, we thus ing epochs, then no such positive contribution to M˙ bh have two possible BH that have survived to the crystal- is lost after t t + τ core, and the BH will continue lization time: those that had M˙ (M = M prim.; τ < ∼ s.g coll. bh bh X sink to sit at the Chandrasekhar mass until t = τsink, at t < τ ) 0 but slightly negative, and arrive at the accum. ∼ which time it definitely loses a large positive contri- crystallization time with a mass still M M prim.; and ˙ bh X bution to Mbh, and will Hawking radiate away within ˙ prim. ∼ those that had Mbh(Mbh = MX ; τsink < t < τaccum.) the WD lifetime. Outcome (3) is avoided in either al- 0 but slightly positive, and arrive at the crystallization∼ ternative. Note however that the boundary at which time with a mass M M prim. + M accum. (ignoring ˙ core bh X X Mbh(Mbh = MChand.; tbh . t . ts.g. + τcoll. ) = 0 de- any WD material accumulated,∼ because there is almost pends on which of the two sets of rate estimates from no parameter space in which that contribution would be Sec. VII A 4 we use, and this has some impact on the re- significant, but the WD is not also fully accreted by this gion in which outcome (3) is ultimately avoided overall; point). we will discuss this in Sec.VIII below. Once again, a positive contribution to M˙ from the Finally, if M˙ (M = M ; t t t + bh bh bh Chand. bh . . s.g. slow CHAMP accumulation turns off at the crystal- τ core) > 0 the initial BH will begin to increase in mass, coll. lization time. For those BH that had M˙ (M = eventually (outside of highly tuned regions where the bh bh M prim.; τ < t < τ ) 0 but slightly neg- Hawking process converts a significant portion of the X sink accum. ∼ ˙ sinking CHAMPs to radiation) reaching the full mass of ative before the crystallization time, Mbh(Mbh = core M prim.; τ < t < τ ) < 0 and there are two the thermal cloud at a time t = ts.g. + τcoll. . It turns out X accum. wd + ˙ outcomes: firstly, the Hawking evaporation can be fast that, for X , whenever Mbh(Mbh = MChand.; tbh . t . core ˙ enough that the BH evaporates within the remaining ts.g. + τcoll. ) > 0, we also have Mbh(Mbh = Ms.g.; ts.g. + core WD lifetime; this avoids outcome (3). Alternatively, τcoll. < t < τsink) > 0, and the BH will always continue to accrete up in mass until it reaches the full primordial the BH can take too long to evaporate in the remain- mass of CHAMPs present in the star around t τ . If ing WD lifetime, resulting in outcome (3). On the other ∼ sink ˙ prim. we are in this region, we have thus far avoided BH evapo- hand, for BHs that had Mbh(Mbh = MX ; τsink < t < τ ) 0 but slightly positive before the crystalliza- ration on too fast a timescale, and we can still ultimately accum. ∼ reach outcome (3). tion time and thus arrive at this point with mass Mbh prim. accum. ∼ To make further progress, we must consider the value MX + MX , there are three outcomes. First, if ˙ prim. M˙ (M = M prim. + M accum.; τ < t < τ ) < 0, of Mbh(Mbh = MX ; τsink < t < τaccum.), which bh bh X X accum. wd now has no further contribution from the primordial the BH begins to Hawking evaporate. It turns out, given ˙ prim. the assumed parameters, that any such BH takes too CHAMPs. If Mbh(Mbh = MX ; τsink < t < τaccum.) < 0, two outcomes are possible: in the first case, the intrin- long to Hawking evaporate within the remaining WD sic timescale for Hawking evaporation is fast enough that lifetime (this of course depends on the WD lifetime and prim. the BH radiates away from Mbh = M to Mbh = 0 could be violated with the WD lifetime was longer, say X ˙ before the WD crystallization time; this again avoids out- 4 Gyr), resulting in outcome (3). Second, Mbh(Mbh = prim. accum. come (3). Alternatively, in the second case, the timescale MX + MX ; τaccum. < t < τwd) > 0 and sufficiently for the Hawking radiation of the BH can be too long large that the BH accretes up to the full WD mass in the ˙ prim. [when Mbh(Mbh = MX ; τsink < t < τaccum.) 0, but remaining WD lifetime; this avoids outcome (3). Third, slightly negative], such that it survives until at∼ least the ˙ prim. accum. Mbh(Mbh = MX + MX ; τaccum. < t < τwd) > 0 WD crystallization time (likely without radiating much but smaller, and the BH takes too long to accrete up to of its mass away, outside of tuned regions of parameter the full WD mass, resulting in outcome (3). space, owing to the strong BH-mass dependence of the This detailed step-by-step analysis, although it is still Hawking rate); such a BH can still end up with outcome crude and can miss some edge cases, provides a series of (3). conditions that must be satisfied for the BH evolution ˙ prim. On the other hand, if Mbh(Mbh = MX ; τsink < t < to be result in outcome (3): i.e., be too slow to either 34

Hawking evaporate away, or accrete up to the full WD In this subsection, we will discuss the backreaction on mass, in the WD lifetime. When we present limits below, the WD as a result of these processes, and their implica- we will exclude the region in which the following are all tions of the survival of old WD. Our discussion roughly satisfied as, taken together, they approximately cover the follows that of Refs. [60, 61]. region in which outcome (3) obtains40 ˙ core (a) Mbh(Mbh = MChand.; tbh . t . ts.g. + τcoll. ) > 0; 1. Evaporating case this is possibly two different conditions depending on which of the rate estimates discussed in In this section we assume that the BH follows an evo- Sec. VII A 4 is used, and we will show the bound- lutionary trajectory such that the final phase of the BH aries of the regions of parameter space for which evolution is one in which the BH dynamics are dominated both of these conditions obtain in turn in Sec.VIII. by evaporation rapid enough in principle to release the ˙ prim. entire BH mass within the remaining WD lifetime [out- (b) Mbh(Mbh = MX ; τsink < t < τaccum.) < 0, evap. prim. come (2) of Sec.VIIB]. We see from the intrinsic Hawk- and τ (M = M 0) > τ , and H bh X → accum. ing radiation timescale estimate Eq. (79) that a necessary τ evap.(M = M prim. 0) > τ τ . Note H bh X → wd − accum. (although, as the discussion at Sec.VIIB makes clear, that is is slightly conservative to impose both the not sufficient) condition is for the BH to enter this phase timescale conditions independently in this fashion, 38 with Mbh, 0 . 10 GeV, yielding an intrinsic evaporation rather than as a single condition. timescale τH . few Gyr, of order the WD lifetime. ˙ prim. A critical result of Ref. [61] is that the evaporation (c) Mbh(Mbh = MX ; τsink < t < τaccum.) > grow prim. of such a BH inside a WD, assuming additionally only 0, and τB+accum.(Mbh = MX Mwd) > → that Mbh, 0 > ET (Mwd) for its WD host, will always τaccum. (timescale assuming Bondi accretion and prim. eventually deposit sufficient energy sufficiently rapidly to slow CHAMP accretion), and τB(Mbh = MX + satisfy the trigger criteria for WD ignition discussed in M accum. M ) > τ τ . X → wd wd − accum. Sec.III. Indeed, assuming that the remaining BH life- Note that, in principle, there is a region satisfying the fol- time is longer than the diffusion timescale for the WD lowing conditions which would not need to be excluded: trigger radius, then even for the largest trigger energies required (for the lightest WD possibly of interest ˙ prim. 25 (d) Mbh(Mbh = MX ; τsink < t < τaccum.) > 0, to us), ET (Mwd 0.8M ) 10 GeV (see Fig.1), we ˙ prim. accum. ∼ ∼ Mbh(Mbh = M + MX ; τaccum. < t < find conservatively that once the BH mass drops below X wd trig. evap. prim. accum. 30 τwd) < 0, and τ (Mbh = M + M Mbh 10 GeV, the total energy Hawking radi- H X X → ∼ 41,42 Mwd) < τwd τaccum.. ated from the BH in one diffusion time exceeds ET . − Since, as discussed at length in Ref. [59], thermalization However, with the lifetimes we have assumed in this dis- of all particles emitted by the BH, except gravitons and cussion and for these WD parameters, region (d) does possibly neutrinos, occurs within a volume not paramet- not exist. rically larger than the trigger volume, this would satisfy the trigger criterion, leading to thermal runaway. For the two specific WD masses that will be of interest, the trig- C. Implications of BH evolution for old WD wd 0.85 trig. 30 ger condition is satisfied for Mbh 4 10 GeV, wd 1.1 trig. 32 ∼ × and Mbh 1.8 10 GeV, respectively. Having discussed in detail the region in which the BH Given that the∼ minimum× diffusion timescale for a evolution does not have sufficient time to proceed to ei- WD in the mass range Mbh [0.8, 1.35]M is τdiff. 13 ∈ ∼ ther devour the entire WD within the WD lifetime, or 2 10− s (see Fig.1), which is the lifetime for a BH of × 28 wd trig. Hawking radiate away to zero mass within the same time- mass Mbh 2 10 GeV < Mbh , it is appropriate frame, we now switch focus to the case where either out- here to integrate∼ × the energy deposition over the full dif- come could in principle happen, absent backreaction on fusion time. Nevertheless, if for any reason the BH were 28 the WD during the BH evolution. either to drop below Mbh 2 10 GeV without yet triggering thermal runaway, or∼ if× it simply forms below that

40 Evaporation timescales are estimated using only the intrinsic Hawking rate, neglecting any lengthening owing to partial cancel- 41 lation against the positive contributions to M˙ , resulting in some TH is so high at this mass that all SM species are radiated ultra- small error in the boundary of the region. Similarly, accretion relativistically; we took geff. = 101.5 in Eq. (77) to conservatively timescales are estimated using only the intrinsic Bondi (and, if exclude gravitons and neutrinos. 42 Ref. [61] quotes a value for a 1.25M WD of applicable, slow CHAMP accretion rates), neglecting any length- wd trig., [61] 35 ening owing to partial cancellation against the negative Hawking Mbh (Mwd = 1.25M ) 2 10 GeV. We would ˙ wd trig. ∼ × 34 contribution to M, resulting in some small error in the bound- find a value of Mbh (Mwd = 1.25M ) 10 GeV, largely ary of the region. We also neglect the Bondi/Eddington crossover because we have a slightly more conservative ∼ trigger condition. here, per the discussion in Sec. VII A 2 a. There is no qualitative difference here. 35

mass, its (remaining) lifetime to radiate away its entire must of course happen prior to the WD being completely (remaining) mass is less than one diffusion time, and it is devoured, the timescale for this to occur must necessar- possible to draw a robust conclusion that, as long as the ily be faster than the accretion timescale estimates given initial mass of the evaporating BH is above the trigger above. Per Ref. [61], for the direct heating of the mate- energy for the host WD (see Fig.1 and Tab.I), runaway rial near around the sonic radius to have a chance to trig- will always be triggered. Because a number of the evolu- ger thermal runaway, that radius would need to exceed 47 tionary trajectories of BHs formed inside the WD result the trigger length; we find that once Mbh & 10 GeV, in the BH immediately Hawking radiating away upon we have λT & RB for a Mwd 0.85M WD (where X ∼ formation (see Sec.VIIB), so that Mbh, 0 MChand., we have estimated the trigger length using the unper- note that the smallest Chandrasekhar mass∼ of any BH turbed density, which is conservative). Similarly, for a + 18 46 we consider occurs for X for mX 10 GeV, and is Mwd 1.1M WD, Mbh & 10 GeV is required for X+ 18 21∼ λ R∼ . Since c 2–3 10 2 at the sonic radius and MChand.(mX = 10 GeV) 5 10 GeV > ET (Mwd = T & B s − ∼ × 1 2 ∼ × 1.1M ). EK = 2 mioncs 2–5 MeV 10Tcrit., the motion of the ∼ ∼ In summary, we have concluded in part that a suf- carbon ions assuming Bondi accretion would in principle ficiently light BH which forms inside a WD of mass be sufficiently fast that thermal runaway could be trig- Mwd & 1.1M would evaporate in a WD lifetime of gered in this ballpark (although, as discussed in Ref. [61], 2 Gyr in large regions of parameter space (see discus- some significant non-radial flow would be required to trig- sion∼ in Sec.VIIB) and, so long as its initial mass upon ger the SN; this is however entirely plausible as carbon beginning to evaporate is above the trigger energy for the ions outside the sonic radius are not necessarily collaps- host WD (see Fig.1 and Tab.I), the evaporation will de- ing radially to the center of the star). posit sufficient energy to trigger thermal runaway in the WD, destroying it in a SNIa-like supernova event. The existence of WD in this mass range older than 2 Gyr However, the above argument ignores the fact that for can thus place limits on this process having occurred.∼ BH of these masses in WD of the relevant central densi- On the other hand, one could also in principle search ties, the BH mass at which RB & λT occurs is already for the SNIa-like destruction events directly. Apart from high enough that we estimate that the accretion is Ed- the caveats regarding the timescale and BH evolutionary dington limited; see Sec. VII A 2 a. As a result, the flow trajectories, this agrees with the results of Ref. [60, 61]. of material onto the WD is choked by radiative backreaction, and cs is an overestimate for the local speed of the ions. As such, it is not clear that the trigger con- 2. Accreting case dition will be reached; a more detailed consideration of this point is beyond the scope of this work, and would In this section we assume that the BH follows an evo- likely require simulations of the accretion dynamics. lutionary trajectory such that the final phase of the BH evolution is one in which the BH dynamics are dominated by accretion of WD material rapid enough in The other mechanism proposed by Ref. [61] is that the principle to devour the entire WD mass within the re- density increase of the WD material in the vicinity of maining WD lifetime [outcome (1) of Sec.VIIB]. We the BH could trigger pycnonuclear fusion. This is how- see from the intrinsic Bondi accretion timescale estimate ever much more speculative. The density increase near Eq. (82) (see also the discussion of the irrelevance of the the BH event horizon under Bondi accretion conditions 2 3 Eddington limitation in Sec. VII A 2 a) that a necessary is 1/cs 10 , so for WD near Mwd 0.8M with ∼ ∼ ∼ (although, as the discussion at Sec.VIIB makes clear, ρ 107 g/cm3, the density could be boosted to near ∼ not sufficient) condition is for the BH to enter this phase ρ 1010 g/cm3 near the BH event horizon, which could 38 ∼ with Mbh, 0 & 10 GeV, yielding an intrinsic accretion trigger pycnonuclear burning [61, 66, 83, 84]. However, it timescale τH . few Gyr, of order the WD lifetime. is unclear that the heated region is either large enough, If this is the way the WD ends, it is not with a bang but or that the flame front would propagate outward rapidly a whimper:43 the BH simply devours the whole WD, but enough, to ignite the star [61]. there is no immediate observational signature of this process, other than the absence of WD older than few Gyr. ∼ The less conservative outcome proposed by Ref. [61] is We therefore reach conclusions similar to, although that the accretion of WD matter onto the BH will be suf- slightly more conservative than, Ref. [61]: if the BH tra- ficiently violent that, in the vicinity of sonic radius from jectory leads it to a final evolutionary stage of accretion, which Bondi accretion occurs, heating or compression of a BH near formed near the center of a WD can destroy the WD material will trigger thermal runaway. Since this the WD within (Gyr), either by direct accretion into the BH, or possiblyO via the triggering of thermal runaway leading to a supernova explosion; a necessary condition is that the initial mass of the BH upon entering this evo- 43 With apologies to T.S. Eliot. 38 lutionary phase is Mbh, 0 & 10 GeV. 36

TABLE III. Some representative old WD from the database in Ref. [67] (and references therein). We give some relevant physical characteristics: the mass in solar masses, surface magnetic field in MG (if known), cooling age in Gyr, effective surface temperature in Kelvin, surface gravitational acceleration in cm/s2, and distance from Earth in parsec. Where multiple values for a specific physical parameter are shown in Ref. [67], we quote the range of available values and refer the reader to the reference for more details. The WD labelled 1–4 are representative non-magnetic WD with masses similar to those of the fiducial WD we considered in Tab.I and cooling ages >Gyr; the WD labelled 5–7 are high-mass magnetic WD with B fields below those discussed in Sec. VI B 2, and cooling ages on the order of Gyr. With the exception of WD 5, which has one of the smallest measured B fields of any WD in this mass range reported in Ref. [67], we have selected these examples largely at random. 4 2 # Name Mwd [M ] B [MG] tcool [Gyr] Teff. [10 K] log g[cm/s ] D [pc] Ref. 10 1 WDJ062144.86+753011.67 1.18–1.23 — 4.1 0.6 9.0 67 [114] 2 WDJ013839.12-254233.40 1.17–1.22 — 4.2 0.7 9.0 70 [115] 3 LP642-052 0.84 — 5.8 0.6 8.4 62 [116] 4 WDJ054706.58+753103.10 0.76–0.87 — 4.0–5.2 0.6 8.3–8.4 75 [117] 5 WD 2051-208 1.21–1.24 0.20–0.26 0.62 1.9–2.1 9.0–9.1 31 [118] 6 WD 0903+083 1.16 6.0 1.24 1.3–4.9 8.0–9.9 142 [119] 7 WD 2202-000 1.08 1.0 2.19 1.0–2.2 8.0–9.0 152 [120]

TABLE IV. Parameters assumed in setting limits. Parameter Symbol Value WD mass Mwd 1.1M 3 WD radius Rwd 7 10− R × 7 3 WD central density ρc 5.5 10 g/cm × WD crystallization time τaccum. 1 Gyr WD (cooling) age/lifetime τwd 2 Gyr WD central temperature T 1 keV 4 WD surface effective temperature Teff 5 10 K ×7 2 WD surface gravitational acceleration g & 10 cm/s WD progenitor cloud mass Mcloud 7 M (range: 6–8 M )[91] 2 Speed of sound in WD center cs 2.8 10− × Effective degrees of freedom for Hawking radiation [Eq. (77)] geff 51/4 (constant) Bondi accretion parameter [Eq. (80)] λ 1 Location — In the disk, Earth vicinity 3 Local DM density ρ0 0.3 GeV/cm Local circular speed vrot. 220km/s WD speed relative to CHAMP halo vwd 220km/s

VIII. LIMITS FROM THE EXISTENCE OF OLD For concreteness, the exact parameters we assume in WHITE DWARFS setting our limits are shown in Tab.IV. Note that our parameter choices are conservative in a number of ways: (1) the WD mass (and hence central density) is not as In this section, we present limits on the galactic high as that of some other old WD [e.g., #1, 2, 5], (2) CHAMP abundance by considering the regions of param- the age of the WD is not as large as that of other old WD eter space in which the WD destruction mechanisms out- [e.g., #1–4], and (3) the magnetic field is large enough lined thus far would have destroyed a number of known to be measurable: if the absence of magnetic field mea- old WD [see, e.g., Tab.III] had such CHAMPs been surements for the older, more massive WD in Tab.III present. is indicative of their magnetic fields being too small to In setting our bounds, we will conservatively consider measure, then those WD could present even more robust a WD such as that labelled #7 in Tab.III: such a WD is bounds. fairly massive, quite old, has a measured magnetic field which is sufficiently small to avoid the concerns expressed in Sec. VI B 2 about the accretion of halo CHAMPs onto the WD over its lifetime up to the crystallization time (at 10 least for masses mX & 10 GeV), and has low enough The basic criterion for obtaining a bound on the exis- effective surface temperature and high enough surface tence of CHAMPs is that the WD at some point in its gravity to avoid the concerns about stellar winds dis- evolution contains at least the larger of a Chandrasekhar cussed in Sec. VI B 2. mass worth of CHAMPs and a self-gravitating mass of 37

CHAMPs (see Sec.VID), 44 either from the primordial refined estimates discussed there are used instead, this contamination, or by accretion of the CHAMPs onto the region does evolve sufficiently rapidly to destroy the WD WD.45 We first consider in turn the two cases where the in its lifetime; see discussion in Sec.VIIB. + WD is contaminated purely by X , and purely by X− The region below the thin dashed red line above mX & in the form of (NX), then we consider what happens if 1011 GeV is where the self-gravitating mass (which ex- both X+ and (NX) are present. ceeds the Chandrasekhar mass in this region) is only ex- Our limits will be presented graphically on the param- ceeded owing to the slow accretion of CHAMPs onto the + eter space (mX , ρX /ρ0), for either fX = 1 or fX− = 1, WD before the crystallization time; in the remainder of and we make comments about how these limits can be the region, the primordial abundance of CHAMPs suf- + applied to the mixed X /X− case. fices to meet the mass criterion. Note in particular that 11 this means that the bounds for mX . 10 GeV are in principle not open to question on the grounds of the WD + A. Limits on X magnetic field (see discussion in Sec. VI B 2) because the bounds in that region of parameter space do not rely on + For the case of pure X contamination, fX+ = 1, we the slow accumulation of CHAMPs over the WD life- have argued that enough energy can in principle be re- time being efficient. We also show as a faint red dotted 8 leased in the collapse of the Chandrasekhar mass core to line for mX & 6.4 10 GeV the region of parameter the BH that the WD trigger criteria can be reached (pos- space in which the total× accreted and primordial CHAMP sibly also true during the earlier self-gravitating cloud abundance would exceed the Chandrasekhar mass, even collapse); conservatively ignoring any early ignition sce- though it is smaller there than the self-gravitating mass, narios though, and assuming that the BH successfully both as an indication of the degree to which limits could forms from the initial mini-WD-like core structure in the be strengthened if a lower core temperature were as- WD, we must demand that the evolutionary timescales sumed, and also because of the caveat expressed in foot- for the BH to either accrete the WD mass (or until a trig- note 28 about the necessity of considering the thermal- ger criterion is possibly reached), or radiate to effectively pressure-supported phase. 5 11 zero mass (or until the trigger criterion is reached), must However, in the region 10 GeV . mX . 10 GeV be sufficiently short. (delimited by the thin long-dashed vertical black lines in The shaded dark red region in Fig.4 indicates the Fig.4), we have indicated our bounds by dotted lines region of parameter space where the larger of the (and lighter shading) to indicate that, in this region, the Chandrasekhar mass and the self-gravitating mass of results of Ref. [52] (and to some extent Ref. [49]) call CHAMPs can accumulate in the WD, and where the dy- into question whether CHAMPs will be present in the namical timescales for the BH evolution are sufficiently region of the MW in which the WD-progenitor proto- short to guarantee WD destruction within the WD life- stellar cloud is collapsing, and thus whether this bound time [the local DM density ρ0 is used here purely as a is actually applicable in that region of parameter space; convenient normalization for the CHAMP density]. The see discussion in Sec.V. In particular, we do not claim 9 unshaded shark-fin-like region between mX 10 GeV a robust bound in this region; see also the discussion in 12 ∼ and mX 10 GeV indicates the region in which a self- Sec.IX below. ∼ gravitating mass of CHAMPs (which exceeds the Chan- Fig.4 is drawn assuming that the primordial CHAMP drasekhar mass in this region) would be present in the contamination η is obtained using a central value for WD (already primordially, as it turns out), but where the the mass of the progenitor± protostellar cloud for this WD evolutionary timescales are too long to reliably conclude of Mcloud 7M [91] (implicitly assuming that the zero- that the WD would be destroyed within its observed age-main-sequence∼ star has the same mass as the collaps- cooling age. The tip of the shark-fin region is shaded ing protostellar cloud used in the computation of η ); green and indicates a region in which the evolutionary however, owing to a fair degree of scatter in the initial-± timescales are too long if the na¨ıve initial CHAMP ac- mass–final-mass relationship for WDs (see, e.g., Figures cretion estimates of Sec. VII A 4 are used; however, if the 1 and 2 of Ref. [91]), this value could easily vary in the range Mcloud 6–8M . Such changes only move the boundaries of∼ various excluded regions by (1) numeri- 44 But note also the caveat at footnote 28; requiring the self- cal factors. O gravitating mass may not be necessary, and it is conservative Note also that in Fig.4, we have conservatively plotted to demand it. 18 45 We will present somewhat conservative limits assuming that T all our bounds only up to mX 10 GeV, at which point ∼ ∼ 18 1 keV is the WD central temperature; see discussion in Sec.VID. the limit on ρX+ fX+ /ρ0 is (ρX+ fX+ /ρ0)limit 3 10− . In principle, the limits we present could thus be strengthened if This is still abundantly safe from the Poisson∼ uncertainty× very old, cold WD have core temperatures T 1 keV, but they on the accumulated number of CHAMPs that was dis- can only be strengthened up to the limitation imposed by requir- cussed Sec. VI B 1: a 10% Poisson uncertainty required ing at least a Chandrasekhar mass of CHAMPs. The maximum ∼ 24 achievable improvement in the parameter range where we give only (ρX+ fX+ /ρ0)Poisson & 10− at this mass. + 18 limits occurs for X at mX 10 GeV, and is a factor of Also shown in Fig.4 is the the region delineated by 4 ∼ (3 10 ) improvement on the bound on ρ /ρ0; see Fig.4. O × X thick gray lines (and shaded gray where not otherwise al- 38

100

2 10−

4 10−

WD 6 M prim. > τ X = 10− MChand., X τ desctruction 10 8 NS limit 0 − M prim. X = /ρ Ms.g.,

x X ρ 10 10−

12 10 M − prim. X CHAMP evacuation region? + ←− −→ M M accum. prim. 14 X X 10− + = M accum. M X s.g., X 16 + = 10− + M X , fX = 1 Chand.,

X 10 18 − 104 106 108 1010 1012 1014 1016 1018 mx [GeV]

+ FIG. 4. Constraints on the abundance of CHAMPs X , ρX /ρ0, as a function of CHAMP mass mX in GeV. This plot is drawn assuming that fX+ = 1 (see discussion in text). Here, ρX is the local CHAMP mass density, and ρ0 is taken to be the 3 local halo DM abundance ρ0 = 0.3 GeV/cm . The dark red shaded region bounded by the thick dark red lines (annotated as prim. accum. prim. ‘MX + MX = Ms.g., X ’ at large mX , and annotated as ‘MX = MChand., X ’ at small mX ) denotes the region where prim. (1) the sum of the primordial CHAMP abundance, MX+ , and the CHAMP abundance accumulated onto the WD from the accum. + galactic CHAMP halo before the WD crystallization time, MX+ , exceed the larger of the Chandrasekhar mass for X or the self-gravitating mass Eq. (62), leading to BH formation; and (2) the subsequent BH evolutionary trajectory ends either with the BH accreting the entire WD mass (or triggering thermal runaway in the process of doing so) [see Sec. VII C 2], or with the BH Hawking evaporating away, and triggering thermal runaway in the process [see Sec. VII C 1]. This bound has a kink at s.g. 8 m 6.4 10 GeV owing to a cross-over between the Chandrasekhar and self-gravitating masses (with Ms.g., X > MChand., X X ∼ × above this value of mX , and smaller below). Nevertheless, we also continue to show the parameters for which the Chandrasekhar prim. accum. prim. accum. mass is exceeded at larger mX by MX +MX with the faint red dotted line (annotated ‘MX +MX = MChand., X ’ at large mX ), as this is the maximum extent to which the limits could improve if either (a) the conservatively high temperature we 3/2 have assumed for the thermal structure (T keV) is lowered [Ms.g., X T ], or (b) the thermal-pressure-supported structure ∼ ∝ 9 12 does not actually form; see footnote 28. The shark-fin-shaped region around mX 10 GeV–10 GeV which is not shaded red ∼ and is annotated ‘τdestruction > τwd’ is the region in which a sufficient CHAMP abundance is present (primordially) in the WD to cause BH formation, but the evolutionary timescales for the BH are too long to destroy the WD within its assumed lifetime. The green shaded part of this shark-fin-shaped region would however evolve rapidly enough to destroy the BH (by Hawking evaporation) if the initial primordial CHAMP accretion rate is assumed to be given by the refined estimates in Sec. VII A 4, prim. instead of by the na¨ıve estimates discussed there. The thin red dashed line (annotated as ‘MX = Ms.g., X ’ at large mX ) delineates the regions of parameter space where the primordial abundance alone is sufficient to exceed the self-gravitating mass (above the red dashed line), and where the primordial abundance must be augmented by the slow accretion of CHAMPs over the WD lifetime (until the crystallization time) in order to exceed a self-gravitating mass of CHAMPs in the star (below the red dashed line). These bounds all assume the WD parameters shown in Tab.IV. Also shown are the NS-derived bounds reported by Ref. [50], with the region above the thick solid grey line excluded (annotated as ‘NS limit’ and shaded grey except—for 5 11 clarity—in the shaded red region). The region 10 GeV . mX . 10 GeV, where the results of Ref. [52] call into question the presence of CHAMPs in the MW, is indicated by the thin long-dashed vertical black lines; in this region, we present all results as dotted lines (and the shading is lighter) to indicate the uncertainty as to their applicability; we do not claim robust bounds 3 here (see discussion in text). CHAMPs lighter than mX . 10 GeV have a sinking timescale in the WD estimated to approach or exceed the WD lifetime and/or crystallization time of Gyr and we conservatively do not present limits here; see Sec.VIC. ∼ 39

ready shaded) where Ref. [50] claim a bound from the de- sults of Ref. [52] call into question our bounds on the struction of old NS owing to the formation of a BH inside grounds that CHAMPs may get blown out of the galaxy the NS, assuming a neutron star mass of MNS = 2M by supernova shock waves, the (pX) may not be so

with a progenitor protostellar cloud with Mcloud = 10M evacuated because, being neutral with binding energies (note that the same objection to the bounds of Ref. [50 ] 25 keV, they may not become efficiently entrained in 5 11 ∼ for 10 GeV . mX . 10 GeV arises as for our WD supernova shock fronts [52]. Moreover, the capture of bounds based on the arguments of Ref. [52]; we indicate (pX) in the initial protostellar cloud is highly inefficient the bounds there as thick, dotted lines, and shade the re- owing to their small cross sections (related to their charge gion more lightly). As already discussed in Sec.I, these neutrality); as such, X− bound as (pX) are unlikely to 10 bounds weaken sharply around mX 10 GeV owing primordially accumulate in the protostellar cloud (see to the Hawking process causing a BH∼ with initial mass discussion in Ref. [50]). However, we have argued that 20 MBH . 4 10− M to evaporate before it can devour the presence of exchange reactions (pX)+N (NX)+p the whole× NS, which in the NS case does not inflict exter- combined with very high WD material densities→ make the nally observable damage on the NS (absent this consid- stopping of (pX) in a WD by the accumulation mecha- eration, their bound would continue to higher masses). nism discussed in Sec.VIB efficient. Moreover, if they approach the WD still in the form of (pX), they are not deflected by WD magnetic fields to nearly the same ex- B. Limits on X− tent as charged particles, so even fairly light CHAMPs could accrete onto WD with large surface B fields. As For (NX) we have argued that there are at least two such, we argue that a bound where the abundance of points46 in the evolution prior to the BH formation where (pX) accumulated onto the WD before the crystallization the WD thermal runaway could be triggered: (1) by py- time exceeds the Chandrasekhar mass (or self-gravitating cnonuclear processes in the trans-Chandrasekhar mass mass, whichever is larger) should still apply, but it is core of the WD, or (2) directly during collapse of the weaker than the cognate accumulation-only bound as- stratified core structure to a BH, if a fraction of the grav- suming a density ρX of X− accumulate onto the WD itational binding energy is transferred to the stellar mat- by the ratio of abundance of (pX) to other (NX) states. ter. We again conservatively ignore these early ignition This fraction is of course difficult to estimate robustly, as mechanisms, and assume that the BH successfully forms it depends on the evolutionary history of the (pX) states. from the initial stratified mini-WD-like core structure in Nevertheless, to indicate the region where this bound the WD. In order to destroy an old WD, we must again would likely be applicable, we draw an orange shaded re- demand a sufficient mass of CHAMPs (as discussed in gion in Fig.5 indicating where at least a Chandrasekhar Sec.VIIIA), and that the evolutionary timescales for the mass of X− that are galactically present as (pX) would BH to either accrete the WD mass (or until a trigger cri- accumulate onto the WD before the crystallization time, 4 terion is possibly reached), or radiate to effectively zero assuming that the (pX) make up a fraction 10− of all mass (or until the trigger criterion is reached), must be (NX) bound states, which is the early-Universe evolu- sufficiently short. tionary estimate (see, e.g., Refs. [42, 46–48]). Imposing these requirements, we find limits on X−

bound as (NX) that, for fX− = 1, are broadly similar to those for X+, but are weaker by a significant factor at s.g. 11 C. Limits when both + and are present low mX . mX 8.3 10 GeV; see Fig.5. This is be- X X− cause the bound∼ in this× region is controlled by the Chan- drasekhar mass, and the (NX) Chandrasekhar mass is 2 The case of most widely applicable physical relevance a factor of 6 36 larger for the homogeneous (NX) + is some mixture of X and X− [bound as (NX)] in mixture as compared∼ ∼ to the X+ case, assuming the mean s.g. the WD, with both species co-existing without anni- ion approach we have thus far taken. For mX & mX , + hilating by virtue of the large Coulomb barrier be- the bounds for (NX) and X are identical, because the + tween the X and the X− bound within the positively self-gravitating mass controls the bound, and the self- charged N nuclear volume. However, our discussion in gravitating mass does not depend on the charge of the Sec. VI E 4 shows that the X+ and (NX) will rapidly CHAMP (or CHAMP bound state); cf. Eq. (62). Note stratify with the X+ sinking toward the center of the that boundaries of the shark-fin region where the evo- core, and it would be conservative to demand that the lutionary timescales are too long also shift slightly as + + + total X mass alone exceeds the X Chandrasekhar compared to the X case. mass. Given that the X+ bounds are, roughly speak- One important point to bear in mind however is that, 2 5 11 ing, a factor of Q(NX) 6 = 36 times stronger than in the region 10 GeV . mX . 10 GeV where the re- ∼ 2 the (NX) bounds (by virtue of the MChand., X QX 8 ∝ + scaling) at low mX . 6.4 10 GeV [where the X Chandrasekhar and self-gravitating× masses cross], so long 46 2 And possibly a third: during the collapse of the thermal-pressure- as f + (Q + /Q ) f f + 1/37 (since X & X (NX) X− ⇒ X & supported structure to a stratified core. + fX + fX− = 1), the limits on ρX /ρ0 are to a good 40

100 ( pX ) 2 10−

4 10− M prim. X = MChand., X WD > τ 6 10− τdesctruction

10 8 NS limit 0 − M prim. X = /ρ Ms.g.,

x X ρ 10 10−

12 10 M − prim. X CHAMP evacuation region? + ←− −→ M M accum. prim. 14 X X 10− + = M accum. M X s.g., X 16 = 10− M X−, fX− = 1 Chand.,

X 10 18 − 104 106 108 1010 1012 1014 1016 1018 mx [GeV]

FIG. 5. As for Fig.4, but showing the constraints on the abundance of CHAMPs X− immutably bound as (NX), assuming s.g. 11 f = 1 and Q(NX) = +6. The kink in the bound at m 8.3 10 GeV is where the Chandrasekhar and self-gravitating X− X ∼ × masses cross for Q(NX) = +6. The additional orange shaded region [annotated as ‘(pX)’] as compared to Fig.4 indicates 5 11 4 where, in the region 10 GeV . mX . 10 GeV, an abundance of (pX) at the level of 10− of the total number abundance of X− would accumulate to a super-Chandrasekhar mass of CHAMPs (larger than the self-gravitating mass for these mX ) in the WD within the crystallization time, assuming that these neutral particles are not evacuated from the galactic disk by supernova shockwaves (and magnetically inhibited from re-entry) [52], but are efficiently captured by dense WD material owing to exchange reactions such as (pX) + N (NX) + p [see discussion in Sec.VIB]. Of course, the actual limit would not → 11 jump discontinuously from the red line to the orange region at mX 10 GeV; there would be some smooth crossover region. 3 ∼ CHAMPs lighter than mX . 10 GeV have a sinking timescale in the WD estimated to approach or exceed the WD lifetime and/or crystallization time of Gyr and we conservatively do not present limits here; see Sec.VIC. ∼ approximation given by so small compared to the total mass of (NX) that the impact of the small central X+ core structure will likely ρX 1 ρX not significantly impact the structure of the (NX) layer + ρ0 limit ∼ fX+ × ρ0 limit,X+ above it. Furthermore, the dense X core would likely 2 2 only aid to make the (NX) core collapse earlier, so it for QX /Q(NX) . fX+ 1 ≤ is conservative in this case to fix the limit on ρX /ρ0 by 8 and mX . 6.4 10 GeV (103) requiring a Chandrasekhar mass of (NX) in the WD, × which leads to the limit X+ where [ρX /ρ0]limit is the limit assuming fX+ = 1, as shown in Fig.4. ρX 1 ρX ρ0 ∼ f × ρ0 Likewise, once f + 1/37, the X− bounds on ρ /ρ limit X− limit,X− X X 0 obtained by requiring a Chandrasekhar mass of (NX) for 1 Q2 /Q2 f 1 − X (NX) . X− ≤ in the core would be stronger than those obtained from 8 + + and m 6.4 10 GeV. (104) X ; moreover, the total mass of X is in this regime X . × 41

In the tuned cross-over region where fX+ anism, our bounds are approximately 2 8 ∼ + (QX /Q(NX)) fX− (and mX . 6.4 10 GeV), one would need to be more careful, but a limit× correct within ρ ρ X 2 X , (106) a factor that is (1) would be set by requiring that an ρ ∼ ρ + O 0 limit 0 limit,X X− Chandrasekhar mass worth of CHAMPs is present in total in the star. whereas the bounds from Ref. [50] may possibly be sig- We propose the limit in the mixed X+/(NX) case can nificantly weakened (although they would not disappear be given conservatively by entirely except maybe in a highly tuned region of parameter space, because of possible differences in accretion " efficiency, etc. for the X [50]). ρ 1 ρ ± X min X , ρ0 limit ∼ fX+ ρ0 limit,X+ # D. Comment on assumed uniform CHAMP 1 ρX (105) distribution 1 fX+ ρ0 − limit,X− We have assumed throughout this paper that the pri- which will capture the true limit correct to within (1) mordial CHAMP abundance is uniformly distributed in 8 O factors, for mX . 6.4 10 GeV. the WD at the time of WD formation. However, some × 11 In the region mX & 8.3 10 GeV [where the (NX) early sinking is to be expected, possibly even during the Chandrasekhar and self-gravitating× masses cross], the evolution of the WD-progenitor star lifetime. Taken to limits from X+ and (NX) become equal under our the logical extreme, this could even cause the primordial conservative assumptions, because the self-gravitating abundance of CHAMPs to collapse to a BH earlier in mass fixes the bound, and it does not depend on the the evolution of the WD progenitor star, before the CO charge of the CHAMP state; cf. Eq. (62), but see core of this progenitor even becomes degenerate. In this footnote 28. In this case, it is again conservative to case, if the BH thus formed evaporates before the core is fix the bound to be given by Eq. (105), which, since degenerate, it may not trigger a thermal runaway at all.

[ρ /ρ ] = [ρ /ρ ] + in this mass range, is However, in almost the entire region where BH evap- X 0 limit,X− X 0 limit,X at most a factor of 2 weaker than either individual bound: oration is responsible for destroying the WD per our maxf min [1/fX+ , 1/(1 fX+ )] = 2. canonical picture, the limits set by the additional X+ { − } 8 CHAMPs accumulated onto the WD after formation are In the intermediate region, 6.4 10 GeV . mX . 8.3 1011 GeV, Eq. (105) will still× give a conservative actually stronger than those coming from the primordial estimate× for the lower envelope of the possible bounded abundance alone. Those limits would not be impacted region; however, this is also roughly the range of param- by early CHAMP collapse and BH evaporation: even if eter space in which the BH evolutionary timescales are the primordial abundance of CHAMPs collapsed early to too long to destroy the WD within its age, so some care is form a BH which radiated away, a second BH could still required to set a rigorous bound in this region, and ob- form long after the WD is born (but before the crystal- tain the correct shark-ﬁn-shaped non-excludable region lization time) owing to accretion of additional CHAMPs (see Figs.4 and5). Conservatively, one can of course onto the WD. choose to ﬁx the lower edge of the bounded region by Moreover, in the region of parameter space where lim- Eq. (105), but avoid setting a bound in the union of the its from the primordial CHAMP abundance are much regions where the evolutionary timescales are too long stronger than those from the additional accumulated assuming pure X+ contamination, and where they are CHAMPs, it turns out that the BH dynamics dictate that too long assuming (NX) contamination. While not ex- the BH would grow in size after formation to devour the act, this will give an estimate of the actual bound that is whole star, instead of evaporating away. In this case too, likely no worse than the uncertainties in where the bound any early sinking is a non-issue, because the BH would should lie owing to factors such as, e.g., the mass of the still be present when the WD is born (if the progenitor progenitor gas cloud, etc. star even survives that long). Our bounds are thus quite robust to having the CHAMPs sink in the WD progenitor Finally, we note that the NS limits from Ref. [50] star prior to the formation of a degenerate core. bound only the asymmetry f f+ f because, in + ≡ | − −| the NS, the X and X− are not inhibited from annihilating by a Coulomb barrier as they are in our case (see discussion in Sec.IV). We place our overall bounds us- IX. SPECULATIONS: CA-RICH GAP ing the more restrictive of the bounds on the individual TRANSIENTS + X and X− contamination components, without assuming any annihilation. In particular, this means that for In this section, we adopt a much more speculative at-

+ fX fX− galactically, which is a highly motivated case titude to make some comments on the so-called ‘Ca-rich that would∼ occur with a CP conserving production mech- Gap Transients’ (hereinafter ‘CaRGT’). 42

The CaRGTs are a class of approximately ten observed drasekhar mass (or a self-gravitating mass, whichever is anomalous supernova-like events which are found to oc- larger) worth of CHAMPs within a crystallization time 9 cur preferentially displaced from their most-likely host τaccum. 10 yr of initial WD formation, a BH would galaxies; see, e.g., Refs. [121–125] and references therein. form in∼ the WD (unless the WD is destroyed earlier as Events are included in this anomalous class based on spe- we have discussed throughout), and in the interesting re- cific criteria (see e.g., Refs. [122–124]) identifying them gion of parameter, accrete up in mass over time. A SN as rapidly evolving, calcium-rich, and faint (luminosities could then be triggered if the explosive mechanisms upon in the ‘gap’ between novae and supernovae) transients; accretion of WD material [61] discussed in Sec. VII C 2 the statistically significant [123, 126] preferentially large are operative. As this picture envisages CHAMP accre- spatial offset from the most likely host is however not one tion that can only occur on the outskirts of galaxies, it of the defining characteristics for inclusion, but is rather could naturally explain the observed spatial distribution feature of the class to be explained. of the CaRGT, as well as their intermediate luminosities, The progenitors of these transient events are as-yet as the SN is triggered in a sub-Chandrasekhar progenitor. not known, although several ‘conventional astrophysics’ Although detailed modeling of the light curves and nucle- explanations have been advanced (see Ref. [124] and ref- osynthetic abundances of such an event are well beyond erences therein), such as low-mass He WDs which are the intended scope of this paper, it is at least plausible detonated by some external perturbation (e.g., tidal de- that such a sub-Chandrasekhar progenitor could explain formation by a NS or BH in a close binary with the WD), the high Ca yield, per recent modeling [70, 127, 128]. and He shell detonation on a relatively low-mass CO WD This picture may however be challenged by the some re- core. cent high estimates for the inferred rates for these events Here, we advance a possible unconventional explana- [125, 126], as it relies on some fairly rare events to be tion: that these events occur as a result of CHAMPs in successful. 5 11 the mass range 10 GeV < mX < 10 GeV accumulat- In order to roughly estimate the plausible region of ing in sufficient quantity onto a sub-Chandrasekhar mass parameter space in which this picture could operate, we WD that they trigger the thermal runaway instability, consider a MWD = 0.85M WD with an extremely low destroying the WD in a supernova explosion. We sug- magnetic field (. 1 kG; see Sec. VI B 2) [see Tab.I for gest this CHAMP mass range as it is where the results of other WD parameters], assuming that the progenitor of Ref. [52] suggest that the CHAMP population in the cen- this WD formed in a region where the CHAMP abun- ter of galaxies (or in the disk) could be depleted [although dance [except perhaps for (pX)] is zero owing to the 47 see comments about (pX) in Secs.V andVIIIB]. mechanism of Ref. [52], so that the progenitor had zero We imagine a WD allowed to form in the inner regions initial CHAMP contamination [the (pX) are not effi- of a galaxy where the CHAMP density is low. This WD ciently captured by the protostellar cloud [50]]. We as- can thus avoid being destroyed initially by protostellar sume that this WD then gets gravitationally ejected from or accreted CHAMP contamination. If this WD then the initial CHAMP-depleted region into an outer region happens to be gravitationally ejected from the inner re- of its parent galaxy (at galactocentric distancer ˜),where 48 gions of the galaxy, it will enter a region in which the the CHAMP abundance is non-zero. For simplicity, we CHAMP density is not depleted, where it can begin to adopt a somewhat na¨ıve model in which the WD, once accrete CHAMPs from the remaining (unejected) virial- ejected, experiences a constant, virialized CHAMP abun- ized CHAMP density in the outer halo. Alternatively, a 9 dance ρX (˜r) for a period of τ τaccum. 10 yr, with WD could simply be born in the outer halo, but these the total attainable CHAMP mass∼ in the star∼ taken to be regions of a galaxy are star-poor. that given by the approximate expression Eq. (34). This Should this WD come to accrete more than a Chan- is of course very schematic. In particular, the CHAMP abundance will not just rapidly turn on as the WD moves in the galaxy, so by approximating the temporal integral in Eq. (28) by ρ (˜r) τ , we make some error which 47 X accum. Of course, the results of Ref. [52] were obtained specifically us- is difficult to quantify× given that we do not know the spa- ing parameters for the MW, so the boundaries of the CHAMP evacuation region may differ in a different system; as this section tial profile of the evacuated CHAMPs (see below) and we is in any event speculative, we simply adopt the MW results and do not specify the WD trajectory. show the allowed parameter space under that assumption. Any Moreover, the approximation Eq. (34) ignores the ve- follow-up study that further examines our speculations would of course have to take this into account. locity distribution of the CHAMPs [which would in any 48 Note that the long timeframe envisaged here provides ample time event likely not be a Maxwellian distribution in the galac- for the WD to move a significant distance in the galaxy: suppose tic rest frame if the CHAMPs have been significantly im- we impose that the WD may take only 10% of the crystalliza- pacted by the expulsion dynamics (see, e.g., Ref. [49])]; 8 tion time ( 10 yrs) to move from the inner CHAMP-depleted however, the exact assumed velocity distribution only region to the∼ outer CHAMP-rich region, so that it still has 90% of the crystallization time available to accrete CHAMPs in the matters when considering much more massive CHAMPs 16 CHAMP-rich region. Assuming that the WD moves with typical (mX & 10 GeV) where there is a question of whether 3 galactic speeds v 10− , it moves a distance of 30 kpc in this the entire distribution can be efficiently captured onto ∼ ∼ time, which is sufficient (see text). the WD. As our goal here is not to be exact, but rather 43

to advance a plausibility argument that some region of parametrize parameter space could be available, we judge these rough approximations to be fit for purpose. ρˆX ρdm(˜r) ρX (˜r) ρ0, (108) In estimating ρX (˜r) we will make the approximation ≡ ρ0 ρ0 that although the inner region of the galaxy is actually whereρ ˆX is an arbitrary constant density normaliza- evacuated of CHAMPs, the CHAMP density in the outer 3 tion,50 and ρ = 0.3 GeV/cm is the local DM density in region still roughly tracks the spatial distribution of the 0 the MW, which is simply used here as a convenient nor- DM halo at r =r ˜ 30 kpc, a representative median pro- malization. Although none of the CaRGT events have jected host–transient∼ offset for the known CaRGTs (see, been observed to occur in our galaxy, we will take MW- e.g., Figure 11 of Ref. [123]), just with a different nor- like parameters for ρ (˜r), assuming that has approxi- malization to the DM; that is, we crudely approximate dm mately the same radial dependence and normalization as that the MW DM abundance (i.e., ρdm(r ) ρ0, where r is the distance from the Earth to the MW⊕ ∼ galactic center);⊕ ρX (r) ρdm(r)Θ[r revac.], (107) ∝ − this is of course a guess about the DM abundance in the where revac. is the radius within which we assume the CaRGT host galaxy, but the observed CaRGT events oc- CHAMPs would be evacuated by the SN shockwaves. cur in galaxies of comparable sizes to the MW, so it is Here we are implicitly assuming either that (a) if the host a well-motivated guess. Note that under these assump- galaxy is a disk galaxy, the evacuation of CHAMPs [52] tions, the value of ρ (˜r) is insensitive at the 20% level dm ∼ (although see Ref. [49]) above and below the disk is still to the choice of NFWc(γ = 1, rs = 20 kpc), NFWc(γ = efficient even absent a large density of supernova shock- 1.3, rs = 20 kpc), Einasto(α = 0.17, rs = 20 kpc), or waves in these baryon-poor regions, or (b) the host galaxy isothermal (rs = 5 kpc) DM halo profile models (see, e.g., is elliptical. These assumptions are required in order that Ref. [129] and references therein for definitions), and can is plausible that a full 3D volume of CHAMPs have been be approximated as evacuated (see below). However, in the latter case this ρdm(˜r) 2 does of course lead to a mismatch: the mechanism of 9 10− . (109) Ref. [52] on which we would need to rely to create this ρ0 ∼ × evacuated region specifically considered the properties of Finally, in making our estimates, we assume that v a MW-like disk galaxy. Although SN shockwaves will wd v (˜r) pGM /r˜ 10 3, which must be correct∼ certainly accelerate CHAMPs regardless of the galaxy 0 galaxy − within∼ an order of magnitude∼ unless the WD has been type, there is a question as to whether there are appro- expelled from the inner region of the galaxy so fast that priate and sufficiently large magnetic fields in elliptical it is no longer gravitationally bound to the galaxy (any galaxies to magnetically inhibit the entry of CHAMPs slower, and the WD would likely not move a sufficient initially external to the evacuation volume into that vol- distance in the galaxy to make our speculative mecha- ume. We simply assume that this is the case and that nism operable, unless the WD being ignited are formed the same range of m are subject to the mechanism in X in the outskirts of the galaxy). the galaxies in which the CaRGT are observed as in the Under these assumptions, we compute the region of MW; if either of these assumptions is not the case, the parameter space (m ,ρ ˆ /ρ ), for the cases of pure X+ picture presented here may break down, and/or require X X 0 and X [bound as (NX)] contaminations, in which an refinement.49 Note further that it is important for a full − initially CHAMP-free WD as described above would ac- 3D region of CHAMPs to be evacuated, because if only crete at least the larger of a Chandrasekhar mass and a a thin 2D slice (e.g., the MW disk) is evacuated then, self-gravitating mass of CHAMPs before the crystalliza- once projection effects are accounted for, it is likely that tion time; see Fig.6, where we show these regions (blue) the observed SNIa-like events triggered by the accumu- with a fuzzy boundary to emphasize that there are large lation of CHAMPs would not exhibit the correct spatial uncertainties in their exact location. Note also that this morphology: a WD would only need to be ejected above region is drawn without regard to the region in which the or below the disk to experience a large CHAMP abun- evolutionary timescales for WD destruction are sufficient dance, regardless if the radial distance from the center of rapid [the region for Mwd 0.85 M being somewhat the galaxy. ∼ different than that for Mwd 1.1M ]. It is still necessary to settle on a specific parametriza- We see that there is plausibly∼ some region of pa- tion convention for the DM abundance. We will rameter space in which this mechanism could clear the first hurdle and provide a trigger mechanism for a sub- Chandrasekhar WD in the spatial location at which 49 Alternatively, and optimistically, there may even be some dis- criminating power here: if different types of galaxies evacuate CHAMPs to differing degrees of efficiency, or in different spatial volumes, a high statistics sample of CaGRT events could even 50 Although not completely arbitrary: if the CHAMP abundance show distinct spatial morphologies of these events in different thus parametrized where evaluated in the MW, and the CHAMPs galaxy types if this is the trigger mechanism. were not blown out of the MW disk, then ρX (r ) =ρ ˆX . ⊕ 44

100 100 (pX ) 2 2 10− 10− 0 0 /ρ /ρ

x 4 x 4

ρ 10− ρ 10− or or

0 6 0 6 10− 10− /ρ /ρ x x ˆ ˆ ρ ρ 8 8 10− 10− + + X , fX = 1 X−, fX− = 1 10 10 10 10 − 107 109 1011 − 107 109 1011 mx [GeV] mx [GeV]

FIG. 6. These plots indicate the regions (blue shaded) of (mX , ρˆX /ρ0) parameter space where our speculations about a possible + trigger mechanism for the CaRGTs are applicable, both for X (left panel) and X− bound as (NX) (right panel). These shaded regions are shown with very fuzzy edges to emphasize the large uncertainties associated with our estimates, and their highly 11 speculative nature. Also shown for mX & 10 GeV are the bounds on (mX , ρX /ρ0) from Figs.4 and5 [see captions there 11 for explanations of the shading]. The limits from Figs.4 and5 in the region mX . 10 GeV have been removed as, for our CaRGT picture to be applicable, the CHAMP expulsion mechanisms of Ref. [52] must be operative, and this invalidates the limits from Figs.4 and5 in that mass range (but see footnote 47 for an important caveat about the mass range of interest); we have nevertheless included very light shaded lines that follow the outlines of those limits in order to guide the eye in a comparison of the regions of parameter space. Note that the CaRGT regions are shown forρ ˆX /ρ0, while the limits from Figs.4 and5 are on ρX /ρ0; in order to show both sets of results on the same plot and for the comparison to be meaningful, we have had to identifyρ ˆX and ρX [see Sec.IX and footnote 50]. The blue shaded regions assume Mwd = 0.85M and τaccum. Gyr; see Tab.I for other parameters. ∼

CaRGT events occur; of course, this is not dispositive, course of its lifetime. These CHAMPs, being extremely both because of the large uncertainties on our estimates massive, sink in the WD, most probably forming first here, and because a careful estimate for the rate of this a thermally pressure-supported structure at the center, process would be required. Nevertheless, taken at face- then later undergoing self-gravitating collapse to a minia- value, the fact that there appears to be some allowed ture CHAMP-contaminated WD-like object at the cen- region is interesting and allows us to speculate that this ter of the WD, provided a sufficient mass of CHAMPs is may be an explanation for these events that would some- present. If a sufficient mass of CHAMPs is present in the what naturally explain their spatial morphology. WD, either primordially upon WD formation, or accreted up over the WD lifetime, this central dense CHAMP- contaminated WD-like object eventually collapses to a black hole. This black hole in general has complicated dy- X. CONCLUSION namics because it can accrete WD matter and CHAMPs to increase its mass, as well as Hawking radiate to reduce In this paper we have considered the impact of an its mass. abundance of extremely massive, stable (or cosmologically long-lived), early-Universe relics with (1) electri- Our detailed study of the possible BH trajectories leads O cal charge on the late-stage evolution and survival of old us to conclude that there are only three physical out- white dwarfs (WD). comes (see also Refs. [60, 61]): (1) the BH survives for a Such charged massive particles (CHAMPs) can come length of time of order the WD age by virtue of a bal- to contaminate old WD in at least two ways: by con- ancing between the mass accretion and mass loss rates taminating the protostellar cloud which collapses to form and intrinsically long timescales for evolution, and noth- the main-sequence WD-progenitor star, thereby ending ing interesting happens observationally for an external up in the WD at the end of the evolutionary phase of observer; (2) the BH will evaporate within the WD life- that massive progenitor, or—provided the WD has a time if it ever forms, with the increasing high Hawking small enough magnetic field to not deflect the incom- radiation power emitted by the BH serving to trigger the ing CHAMPs—by direct accretion on the WD over the thermal runaway instability of the WD material, lead- 45 ing to a SNIa-like explosion; or (3) the BH will accrete tional astrophysics explanations for these events may suf- enough matter within the WD lifetime to either conser- fice. vatively devour the entire WD (with no directly observ- In summary, our work improves astrophysical bounds able signature) or alternatively trigger the thermal run- on the allowed galactic CHAMP abundance by many or- away instability of the WD by heating of accreting carbon ders of magnitude at the highest CHAMP masses mX & ions (or possibly other mechanisms), again resulting in a 1011 GeV, and advances a speculative explanation for a SNIa-like supernova explosion. class of anomalous supernova events. With the exception of the region of parameter space where the evolutionary timescales are too long, the WD destruction mechanisms outlined above impose severe ACKNOWLEDGMENTS bounds (see Figs.4 and5 for our main results) on the al- + lowed galactic abundance of CHAMPs (either X or X−) We thank David Dunsky, Ryan Janish, Daniel Kasen, 11 that for mX & 10 GeV are many orders of magnitude Vijay Narayan, and Stanley Woosley for productive and stronger than existing astrophysical bounds on such par- informative conversations. M.A.F. would like to thank ticles derived from their destruction of old neutron stars the Berkeley Center for Theoretical Physics at the Uni- via a similar accretion mechanism (although without the versity of California Berkeley and Lawrence Berkeley Na- possibility of a SNIa-like explosion). We are able to place tional Laboratory for their long-term hospitality during + bounds on the abundances of X and X− separately, and which much of this work was completed. M.A.F. and not on the residual asymmetry after their annihilation P.W.G. were supported by DOE Grant DE-SC0012012, in the dense core of the WD, because the X− become by NSF Grant PHY-1720397, the Heising-Simons Foun- deeply bound to nuclei in the WD, resulting in nuclear- dation Grants 2015-037 and 2018-0765, DOE HEP Quan- sized bound states that are net positively charged and tISED award #100495, and the Gordon and Betty Moore + which therefore prevent the X and X− from approach- Foundation Grant GBMF7946. S.R. was supported in ing close enough to capture and annihilate. This makes part by the NSF under grants PHY-1818899 and PHY- our bounds even stronger than existing constraints in the 1638509, the Simons Foundation Award 378243 and the regime where the net charge asymmetry of the galactic Heising-Simons Foundation grant 2015-038. M.A.F. ac- CHAMP abundance is zero or small. knowledges the Galileo Galilei Institute for Theoretical Some variations on the above picture are possible, but Physics in Florence, Italy for their support and hospital- they too lead to WD destruction. For instance, if the en- ity, which enabled part of this work to be completed. ergy injection into the WD material at earlier phases of the evolution of the CHAMP core results in enough en- Appendix A: binding energies ergy deposition into the WD material during any of the (NX) collapse phases (thermal-pressure-supported structure to mini-WD-like core structure, or mini-WD-like core struc- In this Appendix, we outline the computation of the ture to BH), the thermal runaway can be triggered. Or, binding energy of a heavy, negatively charged CHAMP 12 in the case of X−, the ( CX−) bound states which form X− with a heavy nucleus N; see also Refs. [9, 34, 87]. the bulk of the dense pre-BH mini-WD-like core structure We model the nucleus as a uniform charged ball of radius R R A1/3 with R = 1.22 fm [130], and charge inside the WD may become so dense that pycnonuclear ≈ 0 0 (density-enhanced) fusion processes between carbon ions QN . Let r be the relative co-ordinate separating the point-like CHAMP of charge Q and the center of can occur at a sufficient rate to lead to thermal runaway; −| X | this can only happen just prior to the collapse of the core the nucleus. The electrostatic potential energy of this to the BH, if at all. configuration is then Additionally, we have speculated that in certain regions of parameter space, the possible WD trigger mech-  αˆ  r > R anism provided by an accreting BH could give a natu- − r V (r) = , (A1) ral explanation for the so-called calcium-rich gap tran- ⇒ αˆ r2 sient supernova events, because this mechanism can trig-  3 0 r R −2R − R2 ≤ ≤ ger sub-Chandrasekhar WD to go off as a supernova, which could naturally explain the Ca-rich spectra and whereα ˆ α QX QN . The fact that the potential energy sub-luminous nature of these events [70]. Moreover, transitions≡ from| | a 1/r potential for r > R to a shifted if the supernova shocks are efficient at depleting the harmonic oscillator potential for 0 r R implies that CHAMP abundance in baryon-rich regions of galaxies the binding energies will be reduced≤ from≤ the na¨ıve Bohr [with magnetic field inhibiting the (re-)entry of (expelled) atom binding energies. For reference, the na¨ıve Bohr CHAMPs] as has been argued in the literature [52], the atom ground state would have Bohr radius and binding spatial morphology of these events, which are observed energy given respectively by [131] to occur preferentially far from the center of their host galaxies, could be naturally explained. This mechanism 1 point 1 2 a = EB = αˆ µ, (A2) is of course highly speculative, and other more conven- αµˆ 2 46 where µ = mX mN /(mX + mN ) mN is the reduced 12.6 of Ref. [131]): mass of the system. ≈ § 1 αˆ u00 u + E u = 0 (A4) −2µ 1 − r 1 B 1 Setting the wavefunction for the system to be Ψ(t, r) = 1 αˆ r2 iEt u(r) m u200 3 2 u2 + EBu2 = 0. (A5) e− r Yl (θ, φ), the 3D Schrödingerequation (SE) of −2µ − 2R − R course reads Suppose we rescale x r/R, define β R/a where a is the Bohr radius as defined≡ for the point≡ charge setup, 2 point 2 and define = 2µR EB EB = EB (/β ), where point ⇒ × EB is the binding energy of the point charge setup. We then have 1 l(l + 1) u00 + V (r) + u + E u = 0,, (A3) 2β −2µ 2µr2 B (u ) + u = u (A6) 1 xx x 1 1 (u ) βx2u = ( 3β)u (A7) 2 xx − 2 − 2 u ( ) = u (0) = 0 (A8) 1 ∞ 2 u1(1) = u2(1) = 0 (A9) u0 (1) = u0 (1) = 0. (A10) where we set E = E with E > 0 the binding energy 1 2 − B B of the system. The exterior solution which obeys the BC is x√ β u1(x) = C1xe− U 1 , 2, 2x√ , (A11) Consider then the 3D SE for l = 0 with u = u1 for − √ 1 r > R and u = u2 for 0 r R, with C -smoothness ≤ ≤ 1 R ∞ zt a 1 b a 1 between u1 and u2 imposed at r = R; the boundary where U(a, b, z) Γ(a)− 0 e− t − (1 + t) − − dt is conditions (BCs) to be imposed on u are u ( ) = 0 and the confluent hypergeometric≡ function. The interior so- 1 ∞ u2(0) = 0 (for discussion of the latter condition see, e.g., lution which obeys the BC is

  1 3β 1/4  D − 1 , x√2β   2 √β −  u2(x) = C2 r , (A12)  2 π 3β 1 3β 1 3β   cos − 1 Γ − + 1 D − 1 , ix√2β1/4   − π 4 √β − 2 √β 2 − √β −  where D(ν, z) is the parabolic cylinder function. Continuity at x = 1 demands that C1 √ β 1 3β 1/4 e− U 1 , 2, 2√ = D − 1 , √2β C2 − √ 2 √β − r 2 π 3β 1 3β 1 3β cos − 1 Γ − + 1 D − 1 , i√2β1/4 , − π 4 √β − 2 √β 2 − √β − (A13) while continuity of the derivative at x = 1 imposes   U 2 β , 3, 2√ C1 √ β − √ e− U 1 , 2, 2√  √ 1 + 2(√ β)  C2 − √ − − U 1 β , 2, 2√ − √ p 1 3β 1 3β = βD − 1 , √2β1/4 √2β1/4D − + 1 , √2β1/4 2 √β − − 2 √β r 2 π 3β 1 3β p 1 3β + cos − 1 Γ − + 1 βD − 1 , i√2β1/4 π 4 √β − 2 √β 2 − √β − r 2 π 3β 1 3β 1 3β + cos − 1 Γ − + 1 i√2β1/4D − + 1 , i√2β1/4 . (A14) π 4 √β − 2 √β 2 − √β 47

h √ i Using Eq. (A13) to eliminate the C1/C2e− U( ) ··· TABLE V. Ground-state binding energies EB [MeV] com- expression that appears on the LHS of Eq. (A14) in favor puted per Eq. (A15), compared to the ground state en- harmonic of the RHS of Eq. (A13) yields a single transcendental ergy level (EB )0 of the approximate harmonic oscillator eigenvalue equation for as a function of β; this must be treatment deﬁned at Eq. (A16). For A & 8 and Z & 4, it is a harmonic good approximation to take EB (E )0; the approxi- solved numerically to ﬁnd the allowed ground state (β), ≈ B with the binding energy constructed as mation is poor for (A, Z) = (4, 2) as the true ground state is localized mostly outside the nuclear volume.

(β) harmonic NEB [MeV] (E )0 [MeV] EB = 2 . (A15) B 2µR β= Q Q αµR 4 | X | N He 0.35 < 0 8 While the above procedure has the virtue of being an Be 1.6 1.4 12C 2.9 2.9 accurate solution to the problem of finding the binding 16O 4.1 4.1 energies of the system with the potential described by 20Ne 5.2 5.2 Eq. (A1), it is numerically cumbersome. We can de- 24Mg 6.1 6.1 velop a relatively accurate approximation by returning 56Fe 10.0 9.9 to Eq. (A1), and considering only the part of the potential for r < R, which looks like the harmonic oscillator potential for a 3D oscillator with fundamental fre- p 3 quency ω0 = α/ˆ (mN R ), but with its energies offset by We consider a simple one-dimensional tunneling prob- 3ˆα/(2R). If we take this picture literally, the binding lem in the following setup: let CHAMPs X be located − − energies of the system should be the energies of the 3D at x = 0, L, and assume that carbon ions C are bound oscillator to each CHAMP± site. We denote the carbon nuclear radius as R; each C wavefunction is thus localized within 3 αˆ 3 r αˆ harmonic R of the locations x = 0, L. Our approximation will (EB )n = + n + 3 . (A16) −2 R 2 mN R treat∼ the (CX) bound states± at x = L as immutable ob- ± However, the classical turning points for the ground jects of charge Q(CX) = +5, but we will track the X at state motion of this system are given by x = x = 0 and the C initially localized around x = 0 individu- 3 1/4 1/4 ∗ ally. Because it is extremely massive, we will treat the X √3(R /(ˆαmN )) R√3(ˆαmN R)− . Butαm ˆ N R αQ Aµ R A1/3 ∼Q A4/3/24 A7/3/48, where we as-∼ at x = 0 as stationary; additionally, because we actually N a 0 ∼ N ∼ imagine that the massive (CX) bound states at x = L sumed QX = 1 and at the last step we assumed QN ± | | 1/4 ∼ are actually localized in a quasi-periodic structure in A/2; therefore, we have x /R √3(ˆαmN R)− ∗ which (CX) are present at x = kL for k , we imagine 4.6A 7/12 1.4(A/8) 7/12 1.07(∼ A/12) 7/12. While∼ k Z − − − also that the (CX) bound states at x = ∈ L are station- the classical∼ turning points of∼ the ground state motion ary. Within the context of a Wentzel–Kramers–Brillouin± for small A therefore lie outside the region where the (WKB) approximation, we will ask for the probability harmonic oscillator treatment is appropriate, we might that the carbon ion initially localized around x = 0 is still hope that the ground state binding energies are rea- able to tunnel to a location x = L 2R, at which lo- sonable numerical approximations to the actual binding cation it is within a nuclear diameter− of the C in the energies Eq. (A15) already for A 8, with the accu- (CX) bound state∼ at x = L, and can undergo a fusion racy of the approximation improving∼ for larger A. This reaction. This approximation is manifestly crude, but it turns out to be a correct conclusion; see Tab.V. We can should obtain the appropriate parametric scalings of the therefore approximate the binding energies as exponential tunneling suppression factor. harmonic EB (EB )0 Let the x-coordinate of the dynamical carbon ion C ≈ that is initially localized around x = 0 be x; accounting 3α QX QN h 1/2i = | | 1 (α QX QN mN R)− . (A17) for the finite charge radius of the C ion in its binding 2R − | | with the X− located at x = 0, the terms in the electrostatic potential that are of interest in this tunneling computation are (for x [ L + 2R,L 2R]) Appendix B: Pycnonuclear fusion rate estimate ∈ − − α Q Q V V (x) = | X | C β(x) In order to estimate the ion number density in the ⊃ − x (NX) core at which pycnonuclear fusion may become | | αQC Q αQC Q relevant, we develop here an estimate of the tunneling + (CX) + (CX) L x L + x suppression for this process; this estimate is similar to − 2αQC Q that developed in Ref. [74], which correctly captures the (CX) (B1) exponential suppression of the pycnonuclear fusion rate − L in the case of an ordinary CO WD (although it obtains the wrong prefactor). where the first line is the interaction energy of the dy- 48

namical C ion with the X− at x = 0 and point for the motion of the ion, and is the tunneling T  exponential. Therefore, the dynamical C ion initially lo- 1 x > R calized around x = 0 will have a probability per unit time  | | to fuse with the C ion bound to the X at x = L of β(x) 2 ; (B2) − ≡ 1 x x  | | 3 x R 2 R − R2 | | ≤ r 2 2 W ψC (x ) , (B4) the second line is the interaction energy of the dynamical ∝ EK mC | ∗ | T C ion with the (CX) states at x = L, and the third is a constant oﬀset in the potential energy± conveniently where we have dropped the S factor because it varies chosen to zero out the potential energy contribution from slowly with energy, and where x is the classical turning ∗ the neighboring (CX) ions at x = 0. point of the C ion motion in the potential Eq. (B1) as- Per the arguments advanced in Ref. [74], the probabil- suming that the ion has the ground state binding energy ity per unit time per ion pair for a pycnonuclear fusion EB appropriate for the (CX) bound state. The tunnel- reaction to occur is given by ing exponential is (note in connection with the sign under the square-rootT that the energy of the system at 51 S(E) 2 the classical turning point is E = EB < 0) W vinc. ψinc. , (B3) − ∼ E | | T Z L 2R where S(E) is the nuclear reaction S-factor which en- 1 − p ln = √2mC dx0 V (x0) + EB. (B5) codes all the nuclear physics, and is usually a slowly −2 T x ∗ varying function of energy (absent resonances), E is the (kinetic) energy of the ion which must tunnel though the Because the C ion in question is in a bound state with p Coulomb barrier to trigger the fusion, vinc 2E/m is the X− located at x = 0, we expect that x R; see the corresponding speed of that ion (in the usual∼ symmet- AppendixA. To make this more precise, let us∗ ∼ expand ric tunneling case, this is technically the relative speed of Eq. (B7) for x . R L [the convenience of the the 2 | | the pair, so m is replaced by the reduced mass µ), ψinc. choice of the constant potential energy oﬀset in Eq. (B1) is the ion wavefunction evaluated at the classical turning| | is now manifest]:

2 2 α Q Q 1 x 2αQC Q(CX) x V ( x R L) | X | C 3 + (B6) | | . ≈ − R 2 − R2 L L2 3 3 α Q Q 1 α Q Q Q(CX) R = | X | C + | X | C 1 + 4 x2; (B7) −2 R 2 R3 Q L3 | X |

2 3 for R L, the correction term in the [ ] bracket on ψC (x ) (√πR)− as the value of the ground state the second line above can be ignored, and··· the potential wavefunction| ∗ | ∼ of the 3D harmonic oscillator evaluated reduces to one of the same form as that we already con- around the classical turning point (taken for this part of sidered at the end of AppendixA in developing the ap- the computation to be x R per AppendixA); (2) we ∗ ∼ 2 proximate treatment of the ground state binding energy estimate mC EK (3/2)mC ω0 1/R [this follows from of the (NX) state, so x R is valid. the discussion around∼ Eq. (A16∝)]; and (3) we will use There is however a mismatch∗ ∼ with the treatment in the approximate ground state binding energy estimate AppendixA, because we are treating the tunneling part Eq. (A17) developed for the 3D problem to estimate EB of the problem here as a 1D problem, whereas we treated in Eq. (B5); but (4) we will otherwise continue to com- the problem using the 3D SE in AppendixA. Follow- pute the tunneling exponential in the 1D approach. ing the treatment of the pycnonuclear rate estimate in Ref. [74], we will largely gloss over this mismatch and For the tunneling part of the problem, we need V (x)+ use a blend of the 1D and 3D results: (1) we estimate E with V (x) from Eq. (B1) for x [R,L 2R]: B ∈ −

51 Note that since we are treating the other carbon ions as rigidly the mass that appears in the tunneling exponential is the carbon ﬁxed to their respective massive CHAMPs in this approximation, mass, not the reduced mass of the carbon-carbon system. 49

2 α Q Q α Q Q 2αQC Q(CX) x V (x [R,L 2R]) + E = | X | C Λ | X | C + (B8) ∈ − B R − x L L2 x2 − α Q Q R R (x/L)2 = | X | C Λ + 2 (1 q) , (B9) R − x L − 1 (x/L)2 − where we have defined nuclei of interest) and have used that Q = Q (CX) C − QX , and defined q QX /QC . Therefore, if we 3 h 1/2i | | ≡ | | Λ 1 (α QX QC mC R)− 1, (B10) parametrize x γR with γ 1, and define u x/L ≡ 2 − | | ∼ and u R/L∗ <≡1 (indeed, typically,∼ u 1), then≡ 0 ≡ 0 (the approximation being numerically satisfied for the

Z 1 2u0 2 1/2 1 p − u0 u ln = 2α Q Q m L2/R du Λ + 2u (1 q) . (B11) −2 T | X | C C − u 0 − 1 u2 γu0 −

Note that for this piece of the computation, we must set where the exponent is correct up to an (1) factor. γ to satisfy Putting together Eq. (B15) with the theO points (1) and (2) just below Eq. (B7), we expect parametrically that 2 1 γ the rate per unit time per unit volume for this fusion Λ + 2u3 (1 q) = 0; (B12) − γ 0 − 1 γ2u2 process is − 0 2 2 2 since u 1 is assumed, the last term can be neglected, ΓC(CX)/V n WC(CX) n R− 0 ∝ (CX) ∝ (CX) and this gives γ Λ, which is again of (1), approxi- h p i ≈ O exp 2 2α Q Q m L2Λ/R . mately consistent with the estimate from the harmonic × − | X | C C oscillator approximation (although the exact numerical (B16) values for x derived in the two different approximations will differ).∗ Note also that the final term in the square Let us trade out R for the nuclear number density: 3 1/3 root becomes maximum at the upper end of the integra- (4π/3)nnucl(2R) 1 R [3/(32πnnucl)] ; and L ∼ ⇒ ∼ 3 tion range, where we can estimate its size as for the (CX) number density: (4π/3)n(CX)L 1 1/3 ∼ ⇒ L 3/(4πn ) , yielding (we drop numerical pref- u2 (CX) actors)∼ 2u0 (1 q) 2 (B13) − 1 u u=1 2u − − 0 2 Γ /V 1 1 4u0 + 4u0 1 q 1 C(CX) = (1 q) − − < . (B14) 2/3 2 − 1 u0 ≈ 2 2 n2 n − ∝ (CX) nucl q Since (1) the integrand in Eq. (B11) vanishes at the lower 2/3 1/3 1/3 exp 4 α QX QC mC Λn− n (3/(4π)) . integration limit owing to a cancellation of the first and × − | | (CX) nucl second terms in the square-root (with the third term be- (B17) ing negligible); (2) this cancellation persists only until u few u0 1, at which point the integrand be- On the other hand, the cognate estimate for the or- comes∼ dominated× by the first term in the square-root; dinary pycnonuclear C–C fusion process from Ref. [74] and (3) this domination by the first term in the square- yields root persists at the (1) level until the upper limit of the integration owingO to Eqs. (B10) and (B14); we 2 2 L0 2 2 ΓCC/V nC WCC nC 3 exp 2(L0) /(x0) , can estimate the integral in Eq. (B11) parametrically as ∝ ∝ (x0) × − R ··· du 1/2 √Λ, up to an (1) factor. Therefore, (B18) parametrically,··· {···} ∼ O where L0 is the average distance between carbon ions, 1/4 h p i √ − exp 2 2α Q Q m L2Λ/R , (B15) and x0 L0/ 2QC (αmC L0/2) is the classical T ∼ − | X | C C turning point≈ for the C× tunneling motion in the cognate 50 computation, such that We would like to know where the volumetric rates are approximately equal: ΓCC/V ΓC(CX)/V . This will al- 2 2 2 3/4 ∼ Γ /V n (L0)− 2Q αm L0 low us to estimate the number density at which the X− CC ∝ C C C q catalyzed pycnonuclear fusion process could trigger run- exp 2 2Q2 αm L ; (B19) away in the collapsing core based on the number density × − C C 0 at which this occurs in an ordinary CO WD. We therefore 1/3 seek if we similarly trade out L0 [3/(4πnC )] , we have (dropping numerical prefactors)∼

3/4 8/3 2 1/3 1/3 Γ /V n 2Q αm n− (3/(4π)) CC ∝ C C C C q 1/3 1/3 exp 2 2Q2 αm n− (3/(4π)) . × − C C C (B20)

3/4 q 8/3 2 1/3 1/3 1/3 1/3 n 2Q αm n− (3/(4π)) exp 2 2Q2 αm n− (3/(4π)) C C C C − C C C q 2 2/3 2/3 1/3 1/3 n n exp 4 α Q Q m Λn− n (3/(4π)) (B21) ∼ (CX) nucl − | X | C C (CX) nucl 3/4 8/3 2/3 2 2 1/3 1/3 ln n n− n− 2Q αm n− (3/(4π)) ⇒ C nucl (CX) C C C q 1/3 1/3 q 2/3 1/3 1/3 2 2Q2 αm n− (3/(4π)) 1 2( Q /Q )Λn− n n . (B22) ∼ C C C − | X | C (CX) nucl C

To make progress, note that if we take the extreme car- additional Heaviside theta function in the expression as- bon number densities shortly before pycnonuclear fusion sumed for (v) [Eq. (31)], dictate that caution must be 10 3 is relevant [74] of around ρ 10 g/cm nC 5 exercised in the integration over velocity in Eq. (28). 32 3 2 ∼ 1/3 ⇒ 1/3 ∼ ×3 The integral over the azimuthal angle is always with- 10 cm− , the factor 2QC αmC nC− (3/(4π)) 10 , while the largest that the logarithmic factor becomes∼ as- out bound and yields a factor of 2π; however, the po- 37 3 lar angle defined by v v = vv cos θ may have a suming nC . n(CX) . nnucl and that nnucl 10 cm− , wd · wd is (1). Therefore, to good approximation,∼ we can sim- restriction, depending on the value of v. There are O ply set the [ ]-bracket on the RHS of the final line of three cases. Case I: for vmax vesc, mw vwd, there is ··· ≤ R 3 − R vmax Eq. (B22) to zero to find the condition for equal volu- no angular restriction on θ, and d v 4π 0 dv in Eq. (28). Case II: for v v → v metric rates: esc, mw − wd ≤ max ≤ vesc, mw + vwd, the integral over v must be broken into 3/2 n(CX) [2( QX /QC )Λ] √nnuclnC (B23) two domains. For v v v there is no an- ∼ | | ≤ esc, mw − wd 0.2√n n (B24) gular restriction, while for vesc, mw vwd < v vmax ∼ nucl C − ≤ 34 3 there is a maximum allowed value of cos θ cos θ (v) 10 cm− , (B25) 2 2 2 ≤ ∗ ≡ ∼ vesc, mw vwd v /(2vwdv). The integral in Eq. (28) − − R R vesc, mw vwd where we used n 5 1032cm 3 as the usual carbon must thus be performed as d3v 4π − dv + X − → 0 ∼ × 37 3 R vmax R cos θ (v) pycnonuclear fusion density [74], and nnucl 10 cm− . ∗ 2π v v dv 1 d cos θ. Case III: for vmax ∼ esc, mw− wd − ≥ vesc, mw + vwd, the integral over v would in principle need to be split into three domains, but one is identically Appendix C: Full expressions for accreting CHAMPs zero. For v v v there is no angular re- ≤ esc, mw − wd striction, while for vesc, mw vwd < v vesc, mw + vwd We give here the full expressions for the mass of there is a maximum allowed− value of cos≤θ cos θ (v) CHAMPs that accrete over the lifetime of the WD; see 2 2 2 ≤ ∗ ≡ vesc, mw vwd v /(2vwdv). The integral in Eq. (28) Sec. VI B 1. − − R R vesc, mw vwd must thus be performed as d3v 4π − dv + The form of the truncated Maxwellian distribution → 0 R vesc, mw+vwd R cos θ (v) 2π dv ∗ d cos θ. The angular integral Eq. (29), and more specifically the Heaviside theta func- vesc, mw vwd 1 tion appearing therein, as well as the appearance of the − − 51 is identically zero for v + v < v v , and so Taking this into account, the accumulated mass of esc, mw wd ≤ max the final result is independent of vmax. CHAMPs for Cases I, III, and II, are, respectively,

 v v 1 v v 1 v + v 1 v2 v2   0 erf wd + erf max wd erf max wd 1 + 0 + wd   − 2 2   vwd v0 2 v0 − 2 v0 × 2 vesc, wd vesc, wd       1 v2 v2 1 v2 (v v )2 v  M I = M approx, III + 0 exp wd 0 exp max − wd 1 + max X± X± 2 2 2 2 ×  √π vesc, wd − v0 − 2√π vesc, wd − v0 × vwd       1 v2 (v + v )2 v   0 exp max wd 1 max   2 2  −2√π vesc, wd − v0 × − vwd 1 ( " 2 #)− vesc, mw 2vesc, mw vesc, mw erf exp 2 (C1) × v0 − √πv0 − v0

 v v 1 v v 1 v 1 v2 v2   0 erf wd + erf max − wd erf esc, mw 1 + 0 + wd   2 2   vwd v0 2 v0 − 2 v0 × 2 vesc, wd vesc, wd       1 v2 v2 1 v2 (v v )2 v   0 wd 0 max wd max   + 2 exp 2 2 exp −2 1 +   √π v − v − 2√π v − v × vwd   esc, wd 0 esc, wd 0      2 2    vesc, mw v 1 vesc, mw  M II = M approx, III 1 + + 0 1 X± X 2 2 ± ×   vesc, wd vesc, wd × − 2 vwd        " #     2  2 ! 2    1 vesc, mw  v 1 v v 1 v    exp  esc, mw esc, mw max max    2 1 + 2 + 1 + 2   −√π − v0 ×  − vwd 3 vesc, wd vwd 3 vesc, wd                 vesc, mwvwd 1 vwd    2 1  − vesc, wd − 3 vesc, mw 1 ( " 2 #)− vesc, mw 2vesc, mw vesc, mw erf exp 2 (C2) × v0 − √πv0 − v0

 v v 1 v2 v2 1 v2 v2   0 erf wd 1 + 0 + wd + 0 exp wd   v v × 2 v2 v2 √π v2 − v2   wd 0 esc, wd esc, wd esc, wd 0  M III = M approx, III X± X± " # ! ×  2 v2 v2 v2 1 v2   exp esc, mw 1 + esc, mw + 0 + wd   2 2 2 2   −√π − v0 vesc, wd vesc, wd 3 vesc, wd  1 ( " 2 #)− vesc, mw 2vesc, mw vesc, mw erf exp 2 (C3) × v0 − √πv0 − v0

where M approx, III is deﬁned at Eq. (34), and all other Appendix D: White dwarf structure X± variables are as deﬁned in Sec. VI B 1.

In this Appendix, we review the computation of white dwarf mechanical structure, and the derivation of the Chandrasekhar limit. We also discuss modiﬁcations in- troduced by CHAMP contamination. 52

1. General problem statement The electron contribution to the total energy density is given by

The mechanical structure of a white dwarf of mass Mwd and radius R is given by the stable, spherically sym- wd ρe(r) = 3Af(xF (r)) (D6) metric equilibrium solutions to Einstein’s equations (al- 2 p though in practice the GR corrections are small, except f(x) x(2x + 1) 1 + x2 arcsinh x. (D7) ≡ − for a near-extremal, Chandrasekhar-mass WD) for a per- fect fluid—T µν = (ρ + P )uµuν P gµν with u2 = 1— − − Imposing plasma neutrality gives the much larger (for all with the equation of state given to good approximation densities relevant for the present work) ion energy density by the fully degenerate electron equation of state, with contribution: overall plasma neutrality additionally assumed to prevail everywhere throughout the star in order to guarantee the 3 electrostatic communication of the electron degeneracy ρi Bµê(r)[xF (r)] (D8) ≡ pressure to the much heavier ions that form the bulk of 3 R gemeµaµe the mass density of the star; see, e.g., Refs. [55, 66]. B , (D9) ≡ 6π2 Working in a spherical coordinate system (t, r, θ, φ) 2 ν(r) 2 1 2 with the line-element ds = e dt (1 + h(r))− dr and we have defined the mean molecular weight per free r2dΩ2 leads to the two independent− equations: the first− electron µ (r) µRµˆ (r) with µR µ (R) by being for the mass enclosed at radial coordinate r, and e ≡ e e e ≡ e the second, the Tolmann–Oppenheimer–Volkhoff (TOV) equation [75, 76], being the GR-corrected version of the 1 X Xj(r)Zj , (D10) Newtonian hydrostatic equilibrium equation: µ (r) ≡ A e j j dM(r) = 4πr2ρ(r) (D1) dr where Zj,Aj, and Xj(r) are, respectively, the charge, dP (r) 1 (ρ(r) + P (r)) M(r) + 4πr3P (r) mass number, and mass fraction at radius r of ion species = . dr −r2M 2 2M(r) j; this definition assumes full ionization, but could be Pl. 1 rM 2 − Pl. corrected for partial ionization if desired. The total mass (D2) density is ρ(r) = ρe(r)+ρi(r). In this version of the computation, the chemical composition must be assumed as The boundary conditions are M(0) = 0 and P (0) = P0, an input; for a pure CO WD, we will assume a completely where P0 is a free parameter giving the central pres- homogeneous mixture, Xj(r) = const., so thatµ ê(r) 1. sure; this defines a one-parameter class of solutions. The ≡ radial-coordinate extent of the star and the Arnowitt– Deser–Misner (ADM) mass are then fixed in terms of the central pressure via the condition P (Rwd) = 0 and M = M(R ); for r R , the solution matches onto 2. Non-dimensional governing equations wd wd ≥ wd the free-space Schwarzschild metric with mass Mwd. Defining the radial-coordinate-dependent electron In practice it is preferable to re-define variables r Fermi momentum to be p (r) = m x (r), so that the p F e F r ξ, M(r) M Mˆ (ξ), φ(ξ) α 1 + [x (r)]2 with α ≡ p 2 0 0 F Fermi energy is EF (r) me 1 + [xF (r)] and the elec- ≡ ≡ ∈ ≡ 3 2 3 (0, 1] a free parameter, and where tron density is ne(r) = (geme/6π )[xF (r)] with ge = 2, the total pressure is given by the fully degenerate elec- r tron pressure (we ignore the ion thermal contribution, αMPl. 2A and the T = 0 corrections to the electron pressure) r0 , (D11) 6 ≡ B π 3 r P (r) = Ag(x (r)) (D3) 4πr0B 2A A 3 F M0 = = 8 MPl., (D12) p α3 π B2 g(x) 3 arcsinh x + x(2x2 3) 1 + x2 (D4) ≡ − g m4 A e e . (D5) in terms of which the governing equations Eqs. (D1) and ≡ 48π2 (D2) can be written as 53

ˆ 3 3 dφ(ξ) 1 (ˆµe(ξ) + 8Cφ(ξ)/α) M(ξ) + Cα ξ gˆ[φ(ξ)/α] = (D13) dξ −ξ2 1 (16C/α)(Mˆ (ξ)/ξ) − dMˆ (ξ) 3/2 = α3ξ2 µˆ (ξ) (φ(ξ)/α)2 1 + 3Cfˆ[φ(ξ)/α] (D14) dξ e −

where equations Eqs. (D13) and (D14) [but retain A, B = 0 in the scaling factors r ,M relating the dimensional6 and hp i 0 0 fˆ(y) = f y2 1 non-dimensional solutions], and assume that the chem- − p p ical composition is homogeneous. Then the governing = y(2y2 1) y2 1 arcsinh y2 1, (D15) − − − − equations reduce to a single simplified equation: and d dφ(ξ) 3/2 ξ2 = ξ2 [φ(ξ)]2 α2 , (D21) hp i dξ dξ − − gˆ(y) = g y2 1 − p p with the boundary conditions φ(0) = 1 and dφ(0)/dξ = = 3 arcsinh y2 1 + y(2y2 5) y2 1, (D16) − − − 0, and with ξ still defined by φ(ξ ) = α. Moreover, ∗ 2 ∗ R 5 R in this regime, Mˆ (ξ) = ξ dφ(ξ)/dξ. For α 1, and C A/B = me/(8µaµe ) 3.4 10− (2/µe ) − 1; the≡ boundary conditions are≈ fixed× by φ(0)× = 1 and as relevant for the extremal configurations of the star, φ(ξ) 1 α for most ξ < ξ , so the approximate gov- Mˆ (0) = 0. The one-parameter class is solution is now erning∼ equation, Eq. (D21), becomes∗ independent of α for defined by the value of α, with the WD radial-coordinate most ξ (indeed, it reduces to the Lane–Emden equation extent and ADM mass given by with a polytropic index of 3 [66]); it is only in the regime ξ ξ , where φ(ξ) α, that the dependence on α enters. Rwd = r0ξ (D17) ∗ ∗ However,≈ in this regime,∼ the governing equation tells us Mwd = M0Mˆ (ξ ), (D18) ˆ 2 ∗ that dM(ξ)/dξ = (d/dξ) ξ dφ(ξ)/dξ 0; therefore, in the only regime where− the solution does∼ depend on where φ(ξ ) α defines ξ . The central density is given α 1, the value of Mˆ (ξ) is approximately stationary. by [cf., Eq.∗ (≡1)] ∗ Taken together, these observations imply that Mˆ (ξ ) is ∗ h i approximately independent of α when α 1, and takes 2 3/2 ˆ ρC ρ(0) = B µê(0)(α− 1) + 3Cf(1/α) (D19) the value Mˆ (ξ ) 1.96. As such, and because the scaling ≡ − ∗ ≈2 2 3/2 factor M0 µe− also does not depend on α explicitly, Bµê(0)(α− 1) . (D20) ∝ ≈ − the limiting physical mass of a WD with µ = 2 can be e 6 The governing equations Eqs. (D13) and (D14) must given to a very good approximation by still be solved numerically; see Eq. (2) for a numerical fit r 2 that gives α to an accuracy of better than 1% for masses 1 6π MPl. ˆ MChand. = MPl. R M(ξ ) (D22) in the range Mwd [0.1, 1.35]M , assumingµ ê = 1 and 2 ge µaµe ∗ R ∈ 2 µe = 2. 1.4 M (2/µe) , (D23) ≈ ×

even for µe 2. We have checked explicitly with the 3. Chandrasekhar mass numerical solution that this scaling is obeyed well, even 19 for masses as large as mX 10 GeV. It is of course well known that a WD has a maximum It is more diﬃcult to develop∼ an understanding of the stable mass [55], the Chandrasekhar mass. For µe = 2 (as scaling of the radius of the extremal star with mX , as we in a CO WD), the maximum mass is found numerically to need to know both how the value of α required to obtain 2 be MChand. 1.42 M , corresponding to α? 4.4 10− . the critical star scales, since r0 α/µe, and also how ≈ 2 4 R 4≈1/3 × ∝ However, since 8C/α (6π ρC /geµa(µe ) ) 1.0 the value of ξ depends on α. However, on the basis of 3 8 3 ≈1/3 R 4/3 ≈ × ∗ 10− (ρC /10 g cm− ) (2/µe ) 1, the govern- the arguments advanced above about the independence ing equations× Eqs. (D13) and× (D14) depend only weakly of the solution to the Lane–Emden equation to the value on the mean molecular mass per free electron, and the of α for α 1 except in the region where ξ ξ , when ≈ ∗ dependence of the mass of the physical solution on µe Mˆ is approximately stationary, we can argue that the is dominated by the scaling factor M0 that relates the value of ξ is a fairly weak function of parameters when non-dimensionalized solution to the physical one. To α 1. It∗ remains to understand the scaling of α with see this, suppose we send C/α 0 in the governing m . → X 54

CHAMPs co-exist, or are bound) is a complete strat- 3 10− ification of the X material (for QX > 0) or the ho- 52 mogeneously mixed (NjX) material (for QX < 0), 4 10− and the ordinary homogeneous CO mixture; this occurs because of the relative buoyancy (more specifically, ? 5 α 10− the much smaller charge-to-mass ratio) of the CHAMP- 1/3 contaminated material as compared to the standard WD 6 mX− 10− ∝ material; see discussion in Sec.VIC. We thus assume as input that the mean molecular mass per electron,µ ê(ξ), 7 10− undergoes a sharp transition (in reality, this would be 106 108 1010 1012 1014 1016 smoothed by thermal effects; see Sec. VI E 1) at a radius ξ = ξ between its value for ordinary WD mat- mX [GeV] X ter (for ξX < ξ ξ ), and its value for the CHAMP- ≤ ∗ contaminated matter (for 0 ξ ξX ). The value of ξX FIG. 7. The value of α necessary to obtain the critical star ≤ ≤ ? is found consistently such that MX M(ξX ) is the total assuming a pure (NX) configuration with Q(NX) = +6, as CHAMP mass in the stratified core.≡ a function of the CHAMP mass mX (black circles). Over a 1/3 In the heavy CHAMP limit, we obtain the mean molec- wide range of masses, α m− (red line), as discussed in ? X ular mass per electron for the CHAMP-contaminated the text. ∝ mixture using AX mX /µa; ZX = QX (for QX > 0) or Z = Q ≡Q (for Q < 0); and X = 1 (Nj X) j − | X | X X The critical stellar configuration is obtained when the (for QX > 0) or X(Nj X) = yj (for QX < 0), where yj EoS for the electrons at the stellar core comes too close is the fraction of negatively charged CHAMPs bound to 4/3 P to the extreme relativistic EoS, P ρ (see, e.g., ion species j (satisfying j yj = 1): Ref. [74]); this in turn means that the∝ condition is re-  ally one on the central electron number density of the mX 1  QX− QX > 0 star, and thus on the central mass density of the star (for  µa µ (0 ξ < ξ ) fixed charge-to-mass ratio massive constituents). Indeed, e X ≤ ≡ mX ¯ 1 a GR fluid stability analysis (see, e.g., Ref. [74]) shows  (QN QX )− QX < 0, µa − | | that the relevant condition is that the central density of (D25) 2 2 2 the extremal star scales as ρC µe mX /QX . From ∝ ∝ 3 Eq. (D20), we see that for α 1, ρ µ α− ; we thus P C e where Q¯N yjQj is the mean charge of the ions to 1/3 1/3 ∝ j expect that α Q m− is a good approximation ≡ ? ∝ X X which the CHAMPs are bound (weighted by the fraction to the scaling of critical value of α with mX , provided of CHAMPs bound to each ion species). Given that we that α 1. We verify numerically over a broad range of assume a CO mixture, Q¯N [6, 8]; moreover, we assume massesm that the m scaling is obeyed within an (1) ∈ X X O a composition of equal mass-fraction abundances, so we numerical factor; see Fig.7. We also verify the QX scal- will assume throughout this work as an approximation ing numerically. As a result, we expect that the radius that Q¯N 7. This approximation of course ignores the 4/3 4/3 ≈ of the extremal star obeys roughly Rwd QX mX− to small difference between the mass-fraction and number within an (1) numerical factor, or ∝ fraction of the ions, and also ignores the differential affin- O ity of CHAMPs to bind to each ion species given their 3 4/3 RChand. 1.5 10− R (2/µe) , (D24) differing binding energies; however, since Q¯ is in any ≈ × × N event bounded in a small range, choosing the middle of where we used the Chandrasekhar radius for an ex- that range introduces only an (1)-factor error in our tremal CO WD as the fiducial value. Note also results. O that this is self-consistent with the scaling of the cen- For the case where the central pressure in the core tral (technically, average) density ρ M /R3 C wd wd vastly exceeds the ambient pressure in the CO material 2 2 4/3 4/3 3 2 2 ∝ ∝ QX mX− /(QX mX− ) mX /QX . Note that this implies that the central ion∝ number density of the extremal 2 star increases as nC mX /QX . We also that the∝ electron mass does not enter 52 Strictly speaking, given sufficient time to settle, the homoge- Eq. (D23) explicitly. neous (Nj X) material mixture will itself stratify into individual layers of pure (Nj X) for each j due to the small differences in the ion charges (and, with the mass being dominated by the 4. CHAMP-modified WD structure CHAMP, the charge-to-mass ratios) of each species; we will ignore this additional complication in this work, as the timescale for a complete stratification of this type is likely extremely long The equilibrium configuration with a heavy CHAMP (the charge differences are only on the order of 20%), and we will (mX 20 GeV, so that we can ignore the C and O ion assume instead that all the (Nj X) species exist as one entirely masses relative to the X mass whenever the ions and homogeneous mixture at the center of the star. 55

105 100 5 105 GeV 10 GeV 4 6 10 106 GeV 10 GeV ] 7 107 GeV 10 GeV NX

[ 3 10 8 8 isolated 10 GeV

/n 10 GeV ] 10 1 /R

N − 2 + 10 core NX R [ n 101

0 10 2 10− 5 3 1 1 3 5 2 1 0 10− 10− 10− 10 10 10 10− 10− 10 M/MChand. n[NX]/n[N]

FIG. 9. For an X− of varying mass m as indicated in the FIG. 8. For an X− of varying mass as indicated in the X legend, this plot shows a comparison of the central ion num- legend, this plot shows a comparison of the radius of the cen- ber densities in the central core of a fully stratified CO star tral core of a fully stratified CO star containing some CHAMP contamination, with the radius of an equal-mass (denoted M) containing some CHAMP contamination [denoted n[NX+N]], with the central ion number densities in an isolated pure- isolated pure-(NX) star. In both cases, the WD structure is computed in the mean-ion approach assuming the ‘mean ion’ (NX) star of the same mass as the core [denoted n[NX]], and the central ion number densities in the uncontaminated (CO) in the unperturbed CO star has charge +7 [resulting in a core of heavy charge +6 (NX) objects]. These results are obtained star [denoted n[N]]. In all cases, the WD structure is computed in the mean-ion approach assuming the ‘mean ion’ in by fixing the total mass of the WD in the stratified case to be Mwd = 1.2M . It is clear that the central core object in the unperturbed CO star has charge +7 (resulting in a core of heavy charge +6 (NX) objects). These results are ob- the stratified star is always more compact than the isolated tained by fixing the total mass of the WD in the stratified completely contaminated star, as a result of the overburden in the stratified case ‘squashing’ the central core; however, as the and uncontaminated cases to be Mwd = 1.2M . It is clear that the central core object in the stratified star is always mass of the core in the stratified case approaches the Chan- more dense than that of either the uncontaminated star or drasekhar mass of the isolated completely contaminated star, the isolated completely contaminated star, as a result of the its radius and that of the isolated completely contaminated overburden in the stratified case. However, as the central star converge. density in the isolated completely contaminated star becomes much larger [respectively, smaller] than the central ion density in the uncontaminated star, the central number density the center of the WD. On the other hand, somewhat un- in the equal-mass core in the partially contaminated strat- surprisingly, less extremal cores behave quite differently ified case behaves more and more like that of the isolated from an their equal-mass isolated WD brethren com- completely contaminated star (the black dotted line shows prised of pure CHAMP-contaminated material: the CO n[NX+N] = n[NX]) [respectively, like that of the uncontam- overburden significantly reduces the radius of the core, inated star (the black dashed line shows n[NX+N] = n[N])]. and forces the central density of the core to be main- The value of the CHAMP mass m appears to be irrelevant X tained at a much higher value than the isolated object; for this comparison, assuming it is large enough to guarantee complete stratification in the partially contaminated case see again Figs.3,8, and9. case. For CHAMPs that are closer in mass to the ion species’ masses, the picture will be more complicated; we do not treat this case. just outside the core (i.e., when the core itself is near- extremal), the solutions that are obtained via the above procedure have the approximate appearance for 0 Appendix E: Electron heat conduction ξ ξ of a small isolated WD comprised of CHAMP-≤ ≤ X contaminated material (i.e., the CO overburden has little We take the electron heat conductivity kcd from impact on the interior solution for near-extremal cores); Ref. [79]. For completeness, we reproduce the relevant see Figs.3,8, and9. This follows because the boundary expressions from Ref. [79] in full here; the formulae in this condition for an isolated WD, P (r = R) = 0, is more section hold in the natural unit system, ~ = c = kB = 1. closely approximated as Pambient/Pcore 0. Therefore, The thermal conductivity for electron conduction is given to maintain a stable hydrostatic equilibrium,→ the strat- by ified inner CHAMP-contaminated material must obey a 2 separate Chandrasekhar limit that can be given approx- π T ne 2 kcd = , (E1) imately as MX . 1.4 M (2Q0µa/mX ) 1.4 M , 3EF νk × where Q0 = QX for QX > 0 and Q0 = Q¯N QX − | | ≈ 3 2 3 7 QX for QX < 0. Stratified cores of CHAMPs that where ne = (geme/6π )xF is the electron density, T is − | | p 2 violate this bound will thus collapse to a black hole in the temperature, EF = me 1 + xF is the electron Fermi 56 energy, and the dominant contribution to the parameter umn on page 6 of Ref. [79] to handle the MCP case], with ei ee νk = νk + νk is given by electron-ion scattering: Zj and nj (respectively) the charge and number density of ion species j.Λk,j is the Coulomb logarithm for ion 2 4πα EF X species j, defined by νei = Z2n Λ , (E2) k p3 j j k,j F j " 2 # 1 pF 2 where α 1/137 is the fine structure constant, p = Λk,j Λk(sj, wj) + Λk(sj, wj) ≈ F ≡ EF mexF is the electron Fermi momentum, and the sum G D , (E3) runs over all ion species present in the multi-component × k,j j plasma (MCP) [we adopt the approximate prescription detailed in the second-last paragraph in the left-hand col- where

1 1 w swh i 2Λ (s, w) s(1 + s)− 1 e− + (1 + sw)e Ei[ sw] Ei[ (1 + s)w] + ln(1 + 1/s), (E4) k ≡ − − − − 2 1 2 1 w swh i 2Λ (s, w) w− + s (1 + s)− 1 e− 1 + s(2 + sw)e Ei[ sw] Ei[ (1 + s)w] + 2s ln(1 + 1/s), (E5) k ≡ − − − − − ηj 2 1 3 ηj Gk,j r 1 + 0.122 βj + 0.0105(1 Zj− ) 1 + (pF /EF ) βj , (E6) ≡ h i2 − 2 3/2 2 1/6 ηj + 0.0081 ηj + 0.19Zj− 4p2 a2 9.1ηj F j Dj exp α0,ju 1e− /4 , α0,j (E7) ≡ − − ≡ 3Γjηj

R ∞ 1 t with Ei(x) x t− e− dt being the exponential inte- and gral function,≡ − and− where 6 2 J(x, y) 1 + + s (q /2p )2, (E8) ≈ 5x2 5x4 j ≡ s,j F 2 " 3 2 wj u 2 (2pF /qD,j) (1 + βj/3) , (E9) y 2.81 0.81 p ≡ − ln 1 + F × 3(1 + 0.07414y)3 y − y E2 with F # π5 y4 η T/ω , β παZ p /E , (E10) + . (E17) j ≡ p,j j ≡ j F F 6 (13.91 + y)4 5/3 αZ 4πn 1/3 a = (3/4πn )1/3, Γ j e , (E11) j j j ≡ T 3 Appendix F: Free-free opacity

p 4α q = 3Γ /a , k2 = E p , The formulae in this section hold in the natural unit D,j j j TF π F F (E12) system, ~ = c = kB = 1. The Rosseland mean free-free opacity is given by the 2 2 β 2 2 1 q = (q + k )e− j , ω 4παZ n m− s,j i,j TF p,j ≡ j j j Kramers opacity [66], suitably corrected [63, 82] by a (E13) velocity averaged (i.e., thermally averaged) Gaunt factor g¯ (T ) that is also appropriately averaged over fre- 2 2 √Γj ﬀ qi,j = qD,j(1 + 0.06Γj)e− . (E14) quencyh i [80, 82]:

Finally, the constants u 1 = 2.8 and u 2 = 13.0. r − − 2 2 2π α3 We also include the sub-dominant contribution to the κff, rad g¯ff (T ) ≈ 1 + 945ζ(7)/π6 15 3 3/2 2 h i parameter νk from electron-electron scattering [79, 132]: me µa 7/2 2 2 3 (1 + X)(X + Y + B)ρT − (F1) 12α T pF × νee = J(x , y), (E15) 2 1 k 3 F 0.1175 cm g− g¯ff (T ) F (ρ, T ) π EF kTF ≈ × h i × (1 + X)(X + Y + B) × where 7/2 ρ T − , (F2) √ 2 8 3 9 y = 3 ωp,e/T, ωp,e = 4παne/EF , (E16) × 10 g cm− × 10 K 57 where we have assumed that Z/A = 1/2 for all elements this definition ensures that the total free-free radiative in the chemical composition, except hydrogen; X and opacity is the Rosseland mean [66, 82]. The thermally av- Y are, respectively, the helium and hydrogen mass frac- eraged Gaunt factorsg ¯ff(T, ν) are tabulated in Ref. [133], P 2 tions (both zero in our case); B = i Zi Xi/Ai is the which results indicate that this correction is usually (1), heavy-element contribution (with the sum running over or within an order of magnitude or two thereof; weO omit all elements heavier than helium), which takes the value this Gaunt factor correction. 12 16 B = 3.5 for a X( C) = X( O) = 0.5 mixture; me is the electron mass; and µa is the atomic mass unit. The fraction in the ( )-bracket on the first line of Eq. (F1) is numerically equal··· to 1.0044. The appropriate frequency-averaged value g¯ (T ) of h ffi the velocity-averaged free-free Gaunt factorsg ¯ff(T, ν) is given by

R ∞ 7 2x x 3 0 dx x e (e 1)− g¯ﬀ(T, ν = xT/2π) g¯ﬀ (T ) R − ; h i ≡ ∞ dx x7e2x(ex 1) 3 0 − − (F3)

[1] J. R. Ellis, J. S. Hagelin, D. V. Nanopoulos, K. A. Olive [14] F. J. Sánchez-Salcedo, E. Mart´ınez-Gómez and and M. Srednicki, Supersymmetric Relics from the Big J. Magaña, On the fraction of dark matter in charged Bang, Nucl. Phys. B238 (1984) 453–476. massive particles (CHAMPs), JCAP 2010 (2010) 031 [2] M. Byrne, C. F. Kolda and P. Regan, Bounds [arXiv:1002.3145]. on charged, stable superpartners from cosmic [15] M. Fairbairn, A. C. Kraan, D. A. Milstead, T. Sjos- ray production, Phys. Rev. D66 (2002) 075007 trand, P. Z. Skands and T. Sloan, Stable Massive [arXiv:hep-ph/0202252]. Particles at Colliders, Phys. Rept. 438 (2007) 1–63 [3] T. Appelquist, H.-C. Cheng and B. A. Dobrescu, [arXiv:hep-ph/0611040]. Bounds on universal extra dimensions, Phys. Rev. D64 [16] T. Alvägerand R. A. Naumann, Search for stable heavy (2001) 035002[ arXiv:hep-ph/0012100]. massive particles of positive integral charge, Phys. Lett. [4] M. Byrne, Universal extra dimensions and 24B (1967) 647–648. charged LKPs, Phys. Lett. B583 (2004) 309–314 [17] R. Middleton, R. W. Zurmuhle, J. Klein and R. V. Kol- [arXiv:hep-ph/0311160]. larits, Search for an Anomalously Heavy Isotope of Oxy- [5] M. B. Wise and Y. Zhang, Stable Bound States of Asym- gen, Phys. Rev. Lett. 43 (1979) 429–431. metric Dark Matter, Phys. Rev. D90 (2014) 055030 [18] P. F. Smith and J. R. J. Bennett, A Search for Heavy [arXiv:1407.4121]. Stable Particles, Nucl. Phys. B149 (1979) 525–533. [6] M. I. Gresham, H. K. Lou and K. M. Zurek, Nuclear [19] P. F. Smith, J. R. J. Bennett, G. J. Homer, J. D. Lewin, Structure of Bound States of Asymmetric Dark Matter, H. E. Walford and W. A. Smith, A Search for Anoma- Phys. Rev. D96 (2017) 096012[ arXiv:1707.02313]. lous Hydrogen in Enriched D2O, using a Time-Of-Flight [7] D. M. Grabowska, T. Melia and S. Rajendran, De- Spectrometer, Nucl. Phys. B206 (1982) 333–348. tecting Dark Blobs, Phys. Rev. D98 (2018) 115020 [20] A. Turkevich, K. Wielgoz and T. E. Economou, Search- [arXiv:1807.03788]. ing for supermassive Cahn-Glashow particles, Phys. [8] Y. Bai, A. J. Long and S. Lu, Dark Quark Nuggets, Rev. D30 (1984) 1876–1880. Phys. Rev. D99 (2019) 055047[ arXiv:1810.04360]. [21] W. J. Dick, G. W. Greenlees and S. L. Kaufman, Search [9] R. N. Cahn and S. L. Glashow, Chemical Signatures for for Anomalous Isotopes of Sodium and Possible Super- Superheavy Elementary Particles, Science 213 (1981) heavy Stable Particles, Phys. Rev. Lett. 53 (1984) 431. 607–611. [22] D. Nitz, D. Ciampa, T. Hemmick, D. Elmore, P. W. [10] A. De Rujula, S. L. Glashow and U. Sarid, Charged Dark Kubik, S. L. Olsen et al., A search for anomalously heavy Matter, Nucl. Phys. B333 (1990) 173–194. isotopes of low Z nuclei, AIP Conf. Proc. 150 (1986) [11] S. Dimopoulos, D. Eichler, R. Esmailzadeh and G. D. 1143–1146. Starkman, Getting a Charge Out of Dark Matter, Phys. [23] W. J. Dick, G. W. Greenlees and S. L. Kaufman, Search Rev. D41 (1990) 2388. for Possible Superheavy Particles in Sodium Nuclei, [12] M. L. Perl, P. C. Kim, V. Halyo, E. R. Lee, I. T. Lee, Phys. Rev. D33 (1986) 32. D. Loomba et al., The Search for stable, massive, ele- [24] E. B. Norman, S. B. Gazes and D. A. Bennett, Searches mentary particles, Int. J. Mod. Phys. A16 (2001) 2137– for Supermassive X− Particles in Iron, Phys. Rev. Lett. 2164[ arXiv:hep-ex/0102033]. 58 (1987) 1403–1406. [13] S. Burdin, M. Fairbairn, P. Mermod, D. Milstead, [25] B. Pichard, J. Rich, M. Spiro, F. Biraben, G. Grynberg, J. Pinfold, T. Sloan et al., Non-collider searches for P. Verkerk et al., On The Possibility Of Detecting Su- stable massive particles, Phys. Rept. 582 (2015) 1–52 perheavy Hydrogen Through Centrifugation And Atomic [arXiv:1410.1374]. Spectroscopy, Phys. Lett. B193 (1987) 383–388. 58

[26] E. B. Norman, R. B. Chadwick, K. T. Lesko, R. M. thesis, Technische UniversitätMünchen, 2009. Larimer and D. C. Hoffman, Search for Supermassive [45] K. Jedamzik and M. Pospelov, Big Bang Nucleosynthe- Cahn-Glashow Particles in Lead, Phys. Rev. D39 (1989) sis and Particle Dark Matter, New J. Phys. 11 (2009) 2499. 105028[ arXiv:0906.2087]. [27] S. Polikanov, C. S. Sastri, G. Herrmann, [46] M. Kusakabe, T. Kajino, T. Yoshida and G. J. Math- K. Lutzenkirchen, M. Overbeck, N. Trautmann ews, New results on catalyzed BBN with a long-lived et al., Search for supermassive nuclei in nature, Z. negatively-charged massive particle, Phys. Rev. D81 Phys. A338 (1990) 357–361. (2010) 083521[ arXiv:1001.1410]. [28] T. K. Hemmick et al., A Search for Anomalously Heavy [47] M. Pospelov and J. Pradler, Big Bang Nucleosynthesis Isotopes of Low Z Nuclei, Phys. Rev. D41 (1990) 2074– as a Probe of New Physics, Ann. Rev. Nucl. Part. Sci. 2080. 60 (2010) 539–568[ arXiv:1011.1054]. [29] P. Verkerk, G. Grynberg, B. Pichard, M. Spiro, S. Zyl- [48] M. Kusakabe, G. J. Mathews, T. Kajino and M.- berajch, M. E. Goldberg et al., Search for superheavy K. Cheoun, Review on Effects of Long-lived Nega- hydrogen in sea water, Phys. Rev. Lett. 68 (1992) 1116– tively Charged Massive Particles on Big Bang Nucle- 1119. osynthesis, Int. J. Mod. Phys. E26 (2017) 1741004 [30] T. Yamagata, Y. Takamori and H. Utsunomiya, Search [arXiv:1706.03143]. for anomalously heavy hydrogen in deep sea water at [49] D. Dunsky, L. J. Hall and K. Harigaya, CHAMP Cosmic 4000-m, Phys. Rev. D47 (1993) 1231–1234. Rays, JCAP 1907 (2019) 015[ arXiv:1812.11116]. [31] D. Javorsek, D. Elmore, E. Fischbach, T. Miller, [50] A. Gould, B. T. Draine, R. W. Romani and S. Nussi- D. Oliver and V. Teplitz, Experimental limits on the nov, Neutron Stars: Graveyard of Charged Dark Matter, existence of strongly interacting massive particles bound Phys. Lett. B238 (1990) 337–343. to gold nuclei, Phys. Rev. D64 (2001) 012005. [51] F. J. Sanchez-Salcedo and E. Martinez-Gomez, Galactic [32] P. Mueller, L. B. Wang, R. J. Holt, Z. T. Lu, constraints on CHAMPs, arXiv:0812.0797. T. P. O’Connor and J. P. Schiffer, Search for [52] L. Chuzhoy and E. W. Kolb, Reopening the win- anomalously heavy isotopes of helium in the earth’s dow on charged dark matter, JCAP 0907 (2009) 014 atmosphere, Phys. Rev. Lett. 92 (2004) 022501 [arXiv:0809.0436]. [arXiv:nucl-ex/0302025]. [53] N. Takibayev and K. Boshkayev, eds., Neutron Stars: [33] Particle Data Group, M. Tanabashi et al., Review Physics, Properties and Dynamics. Nova Science Pub- of particle physics (other particle searches), Phys. Rev. lishers, New York, 2017. D98 (2018) 030001. [54] D. Lynden-Bell and R. Wood, The gravo-thermal catas- [34] M. Pospelov, Particle physics catalysis of thermal trophe in isothermal spheres and the onset of red-giant Big Bang Nucleosynthesis, Phys. Rev. Lett. 98 (2007) structure for stellar systems, Mon. Not. Roy. Astron. 231301[ arXiv:hep-ph/0605215]. Soc. 138 (1968) 495. [35] K. Kohri and F. Takayama, Big bang nucleosynthesis [55] S. Chandrasekhar, The maximum mass of ideal white with long lived charged massive particles, Phys. Rev. dwarfs, Astrophys. J. 74 (1931) 81. D76 (2007) 063507[ arXiv:hep-ph/0605243]. [56] S. W. Hawking, Black hole explosions?, Nature 248 [36] M. Kaplinghat and A. Rajaraman, Big Bang Nu- (1974) 30–31. cleosynthesis with Bound States of Long-lived [57] P. W. Graham, S. Rajendran and J. Varela, Dark Matter Charged Particles, Phys. Rev. D74 (2006) 103004 Triggers of Supernovae, Phys. Rev. D92 (2015) 063007 [arXiv:astro-ph/0606209]. [arXiv:1505.04444]. [37] C. Bird, K. Koopmans and M. Pospelov, Pri- [58] J. Bramante, Dark matter ignition of type Ia mordial Lithium Abundance in Catalyzed Big Bang supernovae, Phys. Rev. Lett. 115 (2015) 141301 Nucleosynthesis, Phys. Rev. D78 (2008) 083010 [arXiv:1505.07464]. [arXiv:hep-ph/0703096]. [59] P. W. Graham, R. Janish, V. Narayan, S. Rajendran [38] M. Kawasaki, K. Kohri and T. Moroi, Big-Bang Nucle- and P. Riggins, White Dwarfs as Dark Matter Detectors, osynthesis with Long-Lived Charged Slepton, Phys. Lett. Phys. Rev. D98 (2018) 115027[ arXiv:1805.07381]. B649 (2007) 436–439[ arXiv:hep-ph/0703122]. [60] J. F. Acevedo and J. Bramante, Supernovae Sparked By [39] K. Jedamzik, The cosmic Li-6 and Li-7 problems and Dark Matter in White Dwarfs, Phys. Rev. D100 (2019) BBN with long-lived charged massive particles, Phys. 043020[ arXiv:1904.11993]. Rev. D77 (2008) 063524[ arXiv:0707.2070]. [61] R. Janish, V. Narayan and P. Riggins, Type Ia super- [40] K. Jedamzik, Bounds on long-lived charged massive novae from dark matter core collapse, Phys. Rev. D100 particles from Big Bang nucleosynthesis, JCAP 0803 (2019) 035008[ arXiv:1905.00395]. (2008) 008[ arXiv:0710.5153]. [62] S. E. Woosley, Type I Supernovae: Carbon Deflagration [41] M. Pospelov, Bridging the primordial A=8 divide with and Detonation, in Supernovae (A. G. Petschek, ed.). catalyzed big bang nucleosynthesis, arXiv:0712.0647. Springer, 1990. [42] M. Pospelov, J. Pradler and F. D. Steffen, Constraints [63] F. X. Timmes and S. E. Woosley, The conductive prop- on Supersymmetric Models from Catalytic Primordial agation of nuclear flames. I—Degenerate C + O and O Nucleosynthesis of Beryllium, JCAP 0811 (2008) 020 + Ne + Mg white dwarfs, Astrophys. J. 396 (1992) 649. [arXiv:0807.4287]. [64] D. Koester and G. Chanmugam, Physics of white dwarf [43] M. Kamimura, Y. Kino and E. Hiyama, Big-Bang stars, Reports on Progress in Physics 53 (1990) 837. Nucleosynthesis Reactions Catalyzed by a Long-Lived [65] P. Montero-Camacho, X. Fang, G. Vasquez, M. Silva Negatively-Charged Leptonic Particle, Prog. Theor. and C. M. Hirata, Revisiting constraints on asteroid- Phys. 121 (2009) 1059–1098[ arXiv:0809.4772]. mass primordial black holes as dark matter candidates, [44] J. Pradler, The long-lived stau as a thermal relic. PhD JCAP 2019 (2019) 031[ arXiv:1906.05950]. 59

[66] R. Kippenhahn, A. Weigert and A. Weiss, Stellar Struc- the Mass and Radius of Dark Matter Particles, Phys. ture and Evolution. Springer, second ed., 2012. Rev. Lett. 64 (1990) 615. [67] P. Dufour, S. Blouin, S. Coutu, M. Fortin-Archambault, [87] B. C. Tiburzi and B. R. Holstein, Bound states of a C. Thibeault, P. Bergeron et al., The Montreal White uniform spherical charge distributions–revisited!, Am. J. Dwarf Database: A Tool for the Community, in 20th Phys. 68 (2000) 640. European White Dwarf Workshop (P. E. Tremblay, [88] R. Foot, Do magnetic fields prevent mirror particles B. Gaensicke and T. Marsh, eds.), vol. 509 of As- from entering the galactic disk?, Phys. Lett. B699 tronomical Society of the Pacific Conference Series, (2011) 230–232[ arXiv:1011.5078]. p. 3, Sheridan Books (Ann Arbor, Michigan), 2017. [89] S. D. McDermott, H.-B. Yu and K. M. Zurek, Turning [arXiv:1610.00986]. off the Lights: How Dark is Dark Matter?, Phys. Rev. [68] P. E. Tremblay, J. Cummings, J. S. Kalirai, B. T. D83 (2011) 063509[ arXiv:1011.2907]. Gänsicke, N. Gentile-Fusillo and R. Raddi, The field [90] D. Clayton, Principles of Stellar Evolution and Nucle- white dwarf mass distribution, Mon. Not. Roy. Astron. osynthesis. University of Chicago Press, 1983. Soc. 461 (2016) 2100–2114[ arXiv:1606.05292]. [91] S. Catalan, J. Isern, E. Garcia-Berro and I. Ribas, [69] S. O. Kepler, S. J. Kleinman, A. Nitta, D. Koester, The initial-final mass relationship of white dwarfs re- B. G. Castanheira, O. Giovannini et al., White Dwarf visited: effect on the luminosity function and mass dis- Mass Distribution in the SDSS, Mon. Not. Roy. Astron. tribution, Mon. Not. Roy. Astron. Soc. 387 (2008) 1693 Soc. 375 (2007) 1315–1324[ arXiv:astro-ph/0612277]. [arXiv:0804.3034]. [70] A. Polin, P. Nugent and D. Kasen, Nebular models of [92] J. N. Bahcall, M. H. Pinsonneault and S. Basu, So- sub-Chandrasekhar mass type Ia supernovae: Clues to lar models: Current epoch and time dependences, neu- the origin of Ca-rich transients, arXiv:1910.12434. trinos, and helioseismological properties, Astrophys. J. [71] L. D. Landau and E. M. Lifshitz, Course of Theoretical 555 (2001) 990–1012[ arXiv:astro-ph/0010346]. Physics, Volume 5: Statistical Physics, Part 1. Perga- [93] A. M. Serenelli, W. C. Haxton and C. Pena-Garay, mon Press, 3rd ed., 1980. Solar models with accretion. I. Application to the so- [72] W. Hillebrandt and J. C. Niemeyer, Type Ia supernova lar abundance problem, Astrophys. J. 743 (2011) 24 explosion models, Annual Review of Astronomy and As- [arXiv:1104.1639]. trophysics 38 (2000) 191–230. [94] N. Vinyoles, A. M. Serenelli, F. L. Villante, S. Basu, [73] D. O. Jones et al., The Discovery of the Most Distant J. Bergström,M. C. Gonzalez-Garcia et al., A new Gen- Known Type Ia Supernova at Redshift 1.914, Astrophys. eration of Standard Solar Models, Astrophys. J. 835 J. 768 (2013) 166[ arXiv:1304.0768]. (2017) 202[ arXiv:1611.09867]. [74] S. L. Shapiro and S. A. Teukolsky, Black Holes, White [95] A. Heger, C. L. Fryer, S. E. Woosley, N. Langer Dwarfs, and Neutron Stars. Wiley-VCH, Weinheim, and D. H. Hartmann, How massive single stars 2000. end their life, Astrophys. J. 591 (2003) 288–300 [75] R. C. Tolman, Static solutions of Einstein’s field equa- [arXiv:astro-ph/0212469]. tions for spheres of fluid, Phys. Rev. 55 (1939) 364–373. [96] R. B. Firestone, Table of Isotopes. Wiley, 1st CD- [76] J. R. Oppenheimer and G. M. Volkoff, On massive neu- ROM ed., 1996. tron cores, Phys. Rev. 55 (1939) 374–381. [97] S. E. Woosley, S. Wunsch and M. Kuhlen, Carbon ig- [77] E. E. Salpeter, Energy and Pressure of a Zero- nition in type Ia supernovae: an analytic model, Astro- Temperature Plasma, Astrophys. J. 134 (1961) 669–682. phys. J. 607 (2004) 921[ arXiv:astro-ph/0307565]. [78] C. Kittel and H. Kroemer, Thermal Physics. W.H. Free- [98] G. Baym and C. Pethick, Physics of neutron stars, man and Company, United States, 1980. Annu. Rev. Astron. Astrophys. 17 (1979) 415–443. [79] A. Y. Potekhin, D. A. Baiko, P. Haensel and D. G. [99] N. W. Evans, C. A. J. O’Hare and C. McCabe, Re- Yakovlev, Transport properties of degenerate electrons finement of the standard halo model for dark matter in neutron star envelopes and white dwarf cores, Astron. searches in light of the Gaia Sausage, Phys. Rev. D99 Astrophys. 346 (1999) 345 [arXiv:astro-ph/9903127]. (2019) 023012[ arXiv:1810.11468]. [80] J. P. Cox and R. T. Giuli, Principles of Stellar Struc- [100] E. Garc´ıa-Berro,L. G. Althaus, A. H. Córsicoand J. Is- ture. Volume 1: Physical Principles. Gordon and ern, Gravitational settling of Ne-22 and white dwarf evo- Breach, London, 1968. lution, Astrophys. J. 677 (2008) 473[ arXiv:0712.1212]. [81] J. P. Cox and R. T. Giuli, Principles of Stellar Struc- [101] I. Renedo, L. G. Althaus, M. M. Miller Bertolami, A. D. ture. Volume 2: Applications to Stars. Gordon and Romero, A. H. Córsico,R. D. Rohrmann et al., New Breach, London, 1968. cooling sequences for old white dwarfs, Astrophys. J. 717 [82] G. B. Rybicki and A. P. Lightman, Radiative Processes (2010) 183–195[ arXiv:1005.2170]. in Astrophysics. Wiley-VCH, Weinheim, 2004. [102] S. Kwok, Planetary nebulae: A modern view, Publi- [83] L. R. Gasques, A. V. Afanasjev, E. F. Aguilera, cations of the Astronomical Society of the Pacific 106 M. Beard, L. C. Chamon, P. Ring et al., Nu- (1994) 344. clear fusion in dense matter: Reaction rate and [103] K. Unglaub, Mass-loss and diffusion in subdwarf b stars carbon burning, Phys. Rev. C72 (2005) 025806 and hot white dwarfs: do weak winds exist?, A&A 486 [arXiv:astro-ph/0506386]. (2008) 923–940. [84] D. G. Yakovlev, L. R. Gasques, M. Beard, M. Wiescher [104] L. Ferrario, D. de Martino and B. T. Gänsicke, Mag- and A. V. Afanasjev, Fusion reactions in multicom- netic White Dwarfs, Space Sci. Rev. 191 (2015) 111–169 ponent dense matter, Phys. Rev. C74 (2006) 035803 [arXiv:1504.08072]. [arXiv:astro-ph/0608488]. [105] E. Bravo, J. Isern, R. Canal and J. Labay, On the con- [85] S. E. Woosley. Private communication, 2019. tribution of Ne-22 to the synthesis of Fe-54 and Ni-58 [86] K. Griest and M. Kamionkowski, Unitarity Limits on in thermonuclear supernovae, Astron. Astrophys. 257 60

(1992) 534–538. 080715.9.html (accessed November 5, 2019). [106] L. Bildsten and D. M. Hall, Gravitational settling of [120] http://www.montrealwhitedwarfdatabase.org/WDs/ 22ne in liquid white dwarf interiors, Astrophys. J. 549 SDSS%20J220435.05+001242.9/SDSS%20J220435.05+ (2001) L219[ arXiv:astro-ph/0101365]. 001242.9.html (accessed November 5, 2019). [107] C. J. Deloye and L. Bildsten, Gravitational settling of [121] H. B. Perets, A. Gal-Yam, P. A. Mazzali, D. Arnett, 22Ne in liquid white dwarf interiors-cooling and seis- D. Kagan, A. V. Filippenko et al., A faint type of super- mological eﬀects, Astrophys. J. 580 (2002) 1077–1090 nova from a white dwarf with a helium-rich companion, [arXiv:astro-ph/0207623]. Nature 465 (2010) 322–325[ arXiv:0906.2003]. [108] J. Burgers, Flow Equations for Composite Gases. Aca- [122] M. M. Kasliwal, S. R. Kulkarni, A. Gal-Yam, P. E. Nu- demic Press, New York and London, 1969. gent, M. Sullivan, L. Bildsten et al., Calcium-rich gap [109] E. Olson and M. Bailyn, Internal Structure of Multicom- transients in the remote outskirts of galaxies, Astrophys. ponent Static Spherical Gravitating Fluids, Phys. Rev. J. 755 (2012) 161[ arXiv:1111.6109]. D12 (1975) 3030–3036. [123] R. Lunnan et al., Two New Calcium-Rich Gap Tran- [110] H. M. van Horn, Cooling of White Dwarfs, in White sients in Group and Cluster Environments, Astrophys. Dwarfs (W. J. Luyten, ed.), Proceedings of the IAU J. 836 (2017) 60[ arXiv:1612.00454]. Symposium no. 42, p. 97, Springer-Verlag (Dordrecht), [124] K. De et al., iPTF 16hgs: A double-peaked Ca-rich gap 1971. transient in a metal poor, star forming dwarf galaxy, [111] D. E. Kaplan and S. Rajendran, Firewalls in Astrophys. J. 866 (2018) 72[ arXiv:1806.10623]. General Relativity, Phys. Rev. D99 (2019) 044033 [125] K. J. Shen, E. Quataert and R. Pakmor, The Progen- [arXiv:1812.00536]. itors of Calcium-Strong Transients, Astrophys. J. 887 [112] H. Bondi, On spherically symmetrical accretion, (2019) 180[ arXiv:1908.08056]. Monthly Notices of the Royal Astronomical Society 112 [126] C. Frohmaier, M. Sullivan, K. Maguire and P. E. Nu- (1952) 195–204. gent, The volumetric rate of calcium-rich transients [113] A. S. Eddington, The Internal Constitution of the in the local universe, Astrophys. J. 858 (2018) 50 Stars. Cambridge University Press, Cambridge, Eng- [arXiv:1804.03103]. land, 1926. [127] S. E. Woosley and D. Kasen, Sub-Chandrasekhar Mass [114] http://www.montrealwhitedwarfdatabase.org/WDs/ Models for Supernovae, Astrophys. J. 734 (2011) 38 Gaia%20DR2%201116006982252923136/Gaia%20DR2% [arXiv:1010.5292]. 201116006982252923136.html (accessed November 5, [128] D. Kasen. Private communication, 2017, 2019. 2019). [129] M. A. Fedderke, E. W. Kolb, T. Lin and L.-T. Wang, [115] http://www.montrealwhitedwarfdatabase.org/WDs/ Gamma-ray constraints on dark-matter annihilation to Gaia%20DR2%205037670392494744320/Gaia%20DR2% electroweak gauge and Higgs bosons, JCAP 1401 (2014) 205037670392494744320.html (accessed November 5, 001[ arXiv:1310.6047]. 2019). [130] K. S. Krane, Introductory Nuclear Physics. John Wiley [116] http://www.montrealwhitedwarfdatabase.org/WDs/ & Sons, Hoboken, NJ, 1988. LP642-052/LP642-052.html (accessed November 5, [131] R. Shankar, Principles of Quantum Mechanics. 2019). Springer, second ed., 1994. [117] http://www.montrealwhitedwarfdatabase.org/WDs/ [132] A. Y. Potekhin, G. Chabrier and D. G. Yakovlev, In- Gaia%20DR2%20502957442209211136/Gaia%20DR2% ternal temperatures and cooling of neutron stars with 20502957442209211136.html (accessed November 5, accreted envelopes., Astron. Astrophys. 323 (1997) 415– 2019). 428 [arXiv:astro-ph/9706148]. [118] http://www.montrealwhitedwarfdatabase.org/ [133] N. Itoh, K. Kojo and M. Nakagawa, Relativistic free- WDs/WD%202051-208/WD%202051-208.html (accessed free Gaunt factor of the dense high-temperature stellar November 5, 2019). plasma. II - Carbon and oxygen plasmas, Astrophys. J. [119] http://www.montrealwhitedwarfdatabase.org/WDs/ Suppl. Ser. 74 (1990) 291–314. SDSS%20J090632.65+080715.9/SDSS%20J090632.65+