Metal pollution onto white dwarfs is a widespread phenomenon that remains puzzling. Some of these white dwarfs also harbor gaseous debris disks. Emission lines from these disks open a unique window to the physical properties of the polluting material, lending insights into their origin. We model the emission line kinematics for the gas disk around SDSS J1228+1040, a system that has been monitored for over two decades. We show that the disk mass is strongly peaked at 1 R⊙ (modulo the unknown inclination), and the disk eccentricity decreases from a value of 0.44 at the inner edge, to nearly zero at the outer edge. This eccentricity profile is exactly what one expects if the disk is in a global eccentric mode, precessing rigidly under general relativity and gas pressure. The precession period is about two decades. We infer that the mass of the gas disk is roughly equivalent to that of a 50 km rocky body, while the mass of the accompanying dust disk is likely insignificant. The disk eccentricity confirms an origin in tidal disruption, while the short disk diffusion time suggests that the disruption event happened a few centuries ago. Moreover, we argue that the initial orbit for the disrupted body, and that of its putative planetary perturber, fall within an astronomical unit around the white dwarf. The total mass of the source population is likely orders of magnitude more massive than our own Asteroid belt and does not seem to exist around main-sequence stars.
1.Introduction
About a third of all white dwarfs show signs of ongoing or recent accretion of heavy metals (e.g., Zuckerman et al. 2003). Starting from the first example of white dwarf G29-38 (Koester 1987; Graham et al. 1990; Jura 2003), many are now known to also exhibit infrared excesses, signs of circ*mstellar dust disks (e.g., Kilic et al. 2005, 2006; Farihi et al. 2009; Debes et al. 2011, 2012). These disks are likely metallic in composition and are responsible for the pollution. It is now commonly believed that large asteroids (and/or comets) around these stars are, for some reason, excited to high eccentricities and are tidally disrupted when they approach the white dwarfs within the Roche radius. The resulting debris forms the dust disk. However, many key elements in this story, including the source and the orbital excitation for these bodies, remain mysterious (for reviews, see Farihi 2016; Veras2016, 2021).
Interestingly, a few percent of these dusty white dwarfs are also known to possess gaseous debris disks (Manser et al. 2020). First discovered from SDSS data by Gänsicke et al. (2006) around the white dwarf SDSS J122859.93+104032.9 (short-named as J1228 below), about two dozen such disks are now known (Gänsicke et al. 2006, 2007, 2008; Farihi et al. 2012; Gentile Fusillo et al. 2021). These disks manifest as double-peaked emission lines in the spectra, most conspicuously in Ca ii infrared triplets. They are found exclusively around white dwarfs hotter than ∼13,000 K, likely because only such stars are luminous enough to sublimate dust at a distance of ∼1 R⊙. The rarity of gaseous disks is likely explained by the rarity of these hot stars (constituting only a few percent of all white dwarfs).
Unlike dust disks which reveal little about their kinematics, compositions, or density distributions (e.g., see a discussion in Nixon et al. 2020), gaseous disks open a lucky window. The characteristic double-peaked emission lines from these disks contain information about their radial extent, eccentricity, surface density, and temperature profiles. Interestingly, some of them exhibit asymmetric lines that are most easily interpreted as the signature of an eccentric disk. Moreover, these disks also appear to be time variable (e.g., Wilson et al. 2014, 2015; Dennihy et al. 2018; Manser et al. 2021). For instance, J1228, the best monitored system (Manser et al. 2016, 2019), shows gradual variations of its line profile over a timescale of decades. This is much longer than the local Keplerian timescale, which is of order hours.
This kinematic information offers the hope of understanding the origin of these debris disks. Here, we undertake a study of the J1228 gaseous disk with the aim of answering the following specific questions.
First, how does the disk manage to retain its eccentric shape? This disk is known to extend radially(with a width Δr ∼ r). If so, general relativistic effects precess the gas differentially and would have led to its total circularization within a few decades. In fact, such a consideration motivates a number of previous works (Hartmann et al. 2011; Metzger et al. 2012; Cauley et al. 2018; Fortin-Archambault et al. 2020) to propose that the observed line profiles are not due to an eccentric disk but are instead due to nonaxisymmetric brightness patterns (a vortex, or a spiral wave) on a circular disk. An alternative possibility is laid out by Miranda & Rafikov (2018) wherein they argue that gas pressure in the disk can resist general relativity (GR) and maintain the disk in a state of rigid precession. We resolve this issue by first fitting a physically motivated disk model (Section 2) against detailed observations (Section 3) and then show that the disk is indeed in a state of pressure-maintained rigid precession (Section 4), as proposed by Miranda & Rafikov (2018).
Second, the kinematic information we extract from the emission lines allows us to set unique constraints on the origin of the gas disk and, more generally, on the origin of white dwarf pollution (Section 5).
2.Physical Model
We will construct a Keplerian disk model to reproduce the observed double-peaked line profiles in the Ca ii infrared triplet. We assume that the gaseous disk region is free of solid grains. As we argue below (Section 5.3), one should not expect dust to coexist with an eccentric gas disk.
In such an exercise, as one only measures the line-of-sight velocity and is ignorant of the orbital period, one can only determine the length combination, , where i is the orbital inclination relative to the line of sight (with i = 0 being face-on). This differs from the usual radial velocity literature, where one also knows the orbital period and can determine
. We suppress the factor
in this section but reintroduce it in later discussions.
Furthermore, in order to translate line flux into gas density, we need to consider the physics of emission. The Ca ii triplets are most likely recombination lines, namely, spontaneous emission from Ca ii ions at excited levels, arriving after photoionization and recombination. The line fluxes should, therefore, scale linearly with the rate of recombination, which, at equilibrium, equals the rate of photoionization. So the emissivity should scale as ∝nγ nCa II, where nγ is the number density of ionizing photons from the white dwarf, and nCa II is that of Ca ii. The above scaling remains valid even when the disk is very optically thick to the recombination lines (as is the case for our disk).
This consideration allows us to determine the local gas density under some simplifying assumptions. In particular, we will assume that nγ depends on the radial distance from the white dwarf as nγ ∝ r−q , with q being a free parameter. We expect q = 2 when the disk is optically thin to the ionizing photons, while q > 2 when the radial optical depth to these photons is larger than unity. We will also assume that the disk has a negligible vertical extent (H) and is not being viewed nearly edge-on (i.e., ) These allow us to determine the local column density from the local height-integrated emission, without worrying about the disk vertical structure. Lastly, we assume that Ca ii perfectly traces the local gas density.
We model the gas disk around J1228 as an assembly of 20 tightly packed, cofocal elliptical Keplerian rings. Their semimajor axes are evenly distributed between ain and aout. The gas surface density profile is assumed to be a broken power law with a transition radius abreak:
Such a broken profile is motivated by the observed line profiles, which also suggest that p1 > 0 and p2 < 0.
As the simplest approximation, we assume that the ring eccentricities vary linearly with the semimajor axis as
We do not specify the sign of the eccentricity gradient (de/da). The rings are assumed to remain apse aligned at all times. This requires the disk to precess rigidly, a working assumption we justify in Section 4.
The line emissivity , from a ring segment of length d ℓ, is therefore
with v being the Keplerian velocity of the segment. The first factor (r−q ) describes the radial dependence for the ionizing flux, and the last factor describes the fractional mass within the line segment. This scales inversely with the local Keplerian velocity as mass is conserved along a Keplerian streamline. This behavior is behind the so-called "apocenter glow" in debris disks (Pan et al. 2016; MacGregor et al. 2017). The overall normalization constant is discussed in Section 5, where we show that our model disk can also account for the total observed flux in addition to the line profile.
To assign a Doppler velocity to the above line segment, we introduce a phase angle ϕ, where ϕ = 0 corresponds to the case where the orbital long axis lies on the plane of the sky. The resultant line profile is symmetric at this phase, while the line is at its most asymmetric when ϕ = π/2.
All our model parameters are illustrated in Figure 1. In the following, we proceed to search for the best model parameters for the observed Ca ii line data from Manser et al. (2016).
3.Markov Chain Monte Carlo and Results
3.1.Markov Chain Monte Carlo
Manser et al. (2016) have gathered spectra of J1228 from 2003 March to 2015 May in a total of 18 epochs using the Very Large Telescope, the SDSS telescope, and the William Herschel Telescope. C. Manser has kindly provided us with the data. Some subsequent observations are presented in Manser et al. (2019) but are not used for fitting.
To prepare the emission profiles for analysis, we convert the Ca ii triplet data from wavelength (Å) to velocity (km s−1) using the atomic rest wavelengths and a systemic velocity of +22 km s−1. This value is within the range reported by Manser et al. (2016): +19 ± 4 km s−1. 4 It is chosen so that the emission lines at the June 2007 epoch, which have very similar amplitudes in the blue- and redshifted peaks, are also symmetric in the velocity space. We coadd the three lines to produce a joint line profile.
The precession of the disk means we have the good fortune to observe it from different vantage points, each giving some unique constraints on the disk model. We decided, initially, to focus on data from three (equally spaced) epochs: 2007 June, 2011 June, and 2015 May. We assign a phase of ϕ = π/2 to the first epoch (most symmetric) and phases of 0.31 and −0.92 to the remaining two. This choice is motivated by our inferred precession period of ∼20 yr (see below). Our initial attempt did not produce a satisfactory fit to the June 2007 data, so we proceeded to include two more epochs, 2007 April and 2007 July, with their corresponding ϕ values, into the procedure. This serves to strengthen the model constraint around ϕ = π/2.
We now determine the best-fit model parameters using the emcee (Foreman-Mackey et al. 2013) implementation of the Markov Chain Monte Carlo (mcmc) method.
This procedure requires appropriate priors. For each of our 8 parameters, we choose a flat prior over a wide range (Table 1). Our prior on the eccentricity profile warrants some comments. First, we posit that e ∈ [0, 1) everywhere. Second, we insist that streamlines cannot cross within the disk. As streamline crossing likely occurs at highly supersonic speed, the resultant shock will remove much of the orbital energy from the gas and cause rapid in-spiral, effectively truncating the disk (i.e., Artymowicz & Lubow 1994). This leads us to impose the condition (Goldreich & Tremaine 1979)
We run emcee with 100 walkers and iterate each for 4000 steps. This ensures that the autocorrelation time is a sufficiently small fraction of the total run. We then trim the first 2000 steps to minimize the effects of the initial conditions. The full results are presented in Figure 4. There is a good convergence for all parameters, though some show slight (but unimportant) bimodal distributions in their posteriors. The maximum-likelihood parameters, and their corresponding uncertainties, are presented in Table 1. Furthermore, the resultant line profiles for the three chosen epochs are presented in Figure 2; while Figure 3 illustrates them for a continuum of phases.
Table 1.Most Likely Model Parameters and their 1σ Uncertainties
Parameter | Prior | Solution | |
---|---|---|---|
Mean (μ) | Uncertainty (σ) | ||
ain [R⊙] | ∈[0.3, 4] | 0.57 | 5.5% |
abreak [R⊙] | ∈[ain, aout] | 1.0 | 2.8% |
aout [R⊙] | ∈[ain, 7] | 1.7 | 4.1% |
ein | ∈[0, 1) | 0.44 | 3.9% |
∇e | Equation (4) | −0.42 | 14% |
q | >0 | 2.4 | 6.8% |
p1 | ∈(0, 5) | 1.8 | 17% |
p2 | ∈( − 5, 0) | −1.9 | 10% |
Download table as: ASCIITypeset image
While the overall comparison is satisfactory, we note that the data show a conspicuous shortage of emission in the blue wing at ϕ = π/2 (bottom panel in Figure 2, yellow box in Figure 3). This is the phase where we expect symmetric emission, and indeed the blue and the red peaks do look symmetric. The deficit is only in the wing, and it persists in observations from nearby epochs (2007 April and July). The latter rules out the possibility that the choice of our ϕ = π/2 epoch is the cause of such an asymmetry. We have no explanation for this deficit.
3.2.Disk Properties
We now review the properties of our best-fit model. Gänsicke et al. (2006) and Hartmann et al. (2016) have previously determined the inner and outer radii for the gas disk: ain ∼ 0.6 R⊙ and aout ∼ 1.2 R⊙. Our solutions are broadly consistent with their values, with ain = 0.57 R⊙ and aout = 1.7 R⊙. Bear in mind that the physical lengths are smaller than these by . Our value for the inner eccentricity (ein = 0.44) is also consistent with that inferred by Manser et al. (2016, 2019) This value is higher than the e = 0.021 value in the original discovery paper (Gänsicke et al. 2006) because they happened to catch the lines when they were more symmetric.
Our most interesting result from the MCMC fit is the eccentricity gradient. We find a significant and negative eccentricity gradient: de/da = −0.42 ± 0.059. Compared to the inner disk, the outer disk is substantially more circular, and is in fact consistent with being circular. This result can be intuitively understood by looking at the top panel of Figure 2: the sharp spike in the blue wing comes about because the rings are compressed together at apoapse (both in physical space and in velocity space). To do so, the inner part of the disk has to be more eccentric. A negative eccentricity gradient is a characteristic feature of a rigidly precessing gas disk (Miranda & Rafikov 2018). We return to this gradient in Section 4.
With our inferred surface density profile (power-law indexes p1 ∼ 1.8, p2 ∼ −1.9), the mass of the disk is strongly concentrated around abreak ∼ 1 R⊙. We also find that q ∼ 2.4 is slightly steeper than our expectation for an optically thin disk (q = 2) and represents a disk that is radially optically thick to the ionizing photons. It is worth commenting that, while one naively expects a degeneracy between q and (p1, p2), since all of them describe radial dependencies (q on the radius, while p on the semimajor axis), Figure 4 convincingly shows that the degeneracy is broken, likely due to the fact that the rings are substantially eccentric.
4.Rigidly Precessing Disk
Here, we first argue that the disk cannot be differentially precessing, as general relativity would have it. We then show, using the tool of linear eigenmode calculations, that it has the appropriate gas pressure to resist differential GR precession. In fact, both the observed precession period and the eccentricity profile agree with theoretical expectations. We thus firmly establish a long-suspected behavior, that the J1228 disk is rigidly precessing.
We first establish a new estimate for the observed precession period. Previously, Manser et al. (2016) reported a period of 24–30 yr based on data up to 2015 May. More recent monitoring extends the data to 2018 May (Manser et al. 2019). From these, one infers that the triplet evolves through a symmetric profile around Oct. 2017. The last time they did so was around 2007 July. This leads us to refine the precession period to a value of 20.5 yr.
We now also insert the factor of where necessary. Previously, Gänsicke et al. (2006) have assigned an inclination of i ∼ 70° to the disk based on the crude arguments that the disk is far from being face-on (double-peaked emission) and is also not edge-on (no self-absorption). This value remains uncertain, so we keep it as a variable.
4.1.General Relativity Makes It Differential
An eccentric ring that is in proximity to the white dwarf experiences GR precession. Let the complex eccentricity be
where ϖ is the longitude of pericentre measured relative to a fixed direction in space. GR acts to advance ϖ at a rate
where is the Keplerian frequency.
Since our gas disk extends radially(with a width Δr ∼ r), the above equation suggests that the eccentric disk should have been markedly twisted after only 10 yr, the precession period for the innermost orbit. If so, streamlines from different orbits could have crossed, and the resulting dissipation should have circularized and shrunk the disk. In contrast, the sharp spikes seen in the line profiles suggest a significant eccentricity. In fact, the disk has been observed to remain in the same eccentric shape for over 20 yr.
For this reason, previous studies have argued that J1228 and other similar white dwarfs do not harbor eccentric disks but instead host circular disks with nonaxisymmetric brightness patterns (e.g., spiral wave, vortex, Hartmann et al. 2011; Metzger et al. 2012). There are also suggestions of eccentric disks but with misaligned apses (Cauley et al. 2018; Fortin-Archambault et al. 2020). These proposals, while being able to provide reasonable fits to the data, are not physically motivated. An asymmetric pattern on a circular disk can be rapidly sheared out on the Keplerian timescale (even faster than GR), and an eccentric disk with misaligned apses is not known to be self-sustaining.
4.2.Pressure Keeps It Rigid
In our work, we opt to model the disk as a series of apse-aligned rings that rigidly precess. We now confirm that this is physically motivated.
First, the radial pressure gradient also causes precession, at a rate that is, to order of magnitude,
for a ring with width Δr ∼ r (Goodchild & Ogilvie 2006; Teyssandier & Ogilvie 2016). To be competitive against GR (Equation (6)), we only need cs /vkep ≥ 0.0024, or a gas temperature 5
where we have evaluated at a = 1 R⊙ and have scaled the mean-molecular weight against that for singly ionized metallic gas (see below). The temperature of the gas disk is likely controlled by photoionization and ranges from 5000 to 9000 K (Melis et al. 2010). We have also confirmed this independently using the spectral synthesis code CLOUDY (Ferland et al. 2017). So the whole disk can easily communicate via pressure, and can smooth out any precessional misdemeanor.
The negative eccentricity gradient we report here supports the hypothesis that the disk is rigidly precessing. Such a configuration means that the rings are more compressed together at their apocenters. The radial pressure gradient there tends to precess the inner streamline backward, while the outer streamlines precess forward. This equilibrates their differential GR rates (Equation (6)).
We now make the above arguments more quantitative. We follow Miranda & Rafikov (2018) to compute the eccentricity eigenmode, the global coherent response of the disk to an eccentricity perturbation. Teyssandier & Ogilvie (2016) have studied the linear response of a locally isothermal, 3D disk. For the case of a power-law disk, where the surface density scales as Σ ∝ rp , and where the temperature also obeys a power law,
their result is simplified by Miranda & Rafikov (2018) into a 1D ordinary differential equation,
where E = E(r) is the (apse-aligned) eccentricity eigenfunction, and wprec the frequency of global precession. The isothermal sound speed is and we adopt μ = 9 (see Section 5).
For our problem, since we only measure the length combination , we transform the above equation to
where and
.
One seems to have some flexibility in choosing the boundary condition (see, e.g., Miranda & Rafikov 2018). We adopt the following set
These differ from that in Miranda & Rafikov (2018), where they also took ∂E/∂r = 0 at the outer boundary. Our adopted one is more descriptive of our best-fit solution. In any case, this does not much affect the precession rate.
For our broken power-law disk, we integrate Equation (11) from the two boundaries toward and insist that E and ∂E/∂r remain continuous across
. We then look for eigenmodes of the lowest radial order. These have the smoothest eccentricity profiles and hence the lowest dissipation rates. They are most likely to persist.
We present in Figure 5 the results of these calculations, adopting the best-fit model parameters from Table 1. As we do not have information on the values of Tin and γ, we experiment within some sensible ranges. These calculations are performed for i = 90°, but the conclusion remains largely the same, given our uncertainties in Tin and γ, for inclinations as low as ∼60o .
We conclude with two findings. First, within the relevant range of Tin, and sensible temperature profiles (γ), we find that the theoretical modes have periods comparable to the observed value (20.5 yr). Second, even though we have only modeled the disk using a linear profile (the simplest choice), this profile agrees very well with the shape of the linear eigenfunctions. These two quantitative agreements strongly support the hypothesis that J1228 hosts a rigidly precessing gas disk under the combined effects of GR and gas pressure.
5.Insights on Origin
We now discuss what our results imply for the origin of the gas disk around J1228, as well as for the pollution of white dwarfs in general. We base our discussions on the current favorite model for white dwarf pollution, tidal disruption of highly eccentric planetesimals (see reviews by Farihi 2016; Brouwers et al. 2022).
5.1.Disk Mass and the Progenitor Mass
We estimate a mass for the gas disk, Mgas, by assuming that the viscous spreading of the gas disk supplies the observed accretion onto J1228.
The gas scale height is
If the disk is accreting under a constant viscosity parameter α (Shakura & Sunyaev 1973), the viscous diffusion time is
We adopt the accretion rate as determined by Dwomoh & Bauer (2023) 6 and estimate the current disk mass by assuming , to obtain
We argue that this estimate likely also reflects the total mass of the original disrupted planetesimal. The current disk is significantly eccentric — it avoids rapid streamline crossing and circularization by organizing itself into a coherent eccentric mode. But over the viscous timescale, the eccentricity should be gradually damped as the disk spreads radially. So the current disk has likely weathered no more than a few viscous times. Or its current mass is close to its original mass. If true, this means the original disrupted planetesimals have a radius Rp ∼ 50 km (assuming a bulk density of ρp ∼ 5 g cm−3).
To substantiate the above estimate for the disk mass, we also consider whether it is consistent with the observed Ca ii emissions. Let us adopt the same chemical composition as that for CI chondrite from Palme et al. (2014), and assume that all metals are singly ionized, the mean nuclei weight of the gas is 15 and the mean-molecular weight is 9 (in units of hydrogen mass). So the above mass estimate corresponds to a surface density Σ ∼ Mgas/r2 ∼ 0.13 g cm−2, and a midplane density of ρ ∼ 3 × 10−10 g cm−3. In the meantime, Ca nuclei are only 1/280 of the total nucleus number. We therefore arrive at a radial column density for Ca of ∼3 × 1021 cm−2, and a vertical one at ∼2 × 1019 cm−2. We find:
1.
if all Ca is in Ca ii, 7 , with a photoionization cross section of 5 × 10−19 cm2, the optical depth for a Ca ii ionizing photon is τ ∼ 1600. So the gas midplane is very opaque to the Ca ii ionizing photons. This may explain why we obtain a steeper fall-off of the ionizing flux (nγ ∝ r−2.4) than is expected for an optically thin disk (nγ ∝ r−2).
2.
the vertical optical depth for the Ca ii infrared triplet is ∼106 (CLOUDY result). This explains why the line ratios within the Ca ii triplet do not reflect their individual strengths but approach those of a blackbody (Melis et al. 2010).
3.
the above mass estimate allows us to explain the total flux observed in Ca ii triplet. For J1228, about 5.6% of its energy is in photons that can ionize Ca ii (ionization energy 11.87 eV). Let the disk be optically thick to these photons from the midplane up to n scale heights. The total energy intercepted by Ca ii is then 2n × (H/r) × 5.6% × Lwd. As Ca ii is ionized and recombined, a fraction of the ionization energy is emitted in Ca ii triplets (photon energy ∼1.5 eV). So we expect a total line flux of ∼2n(H/r) × 5.6% × (1.5/11.87) × Lwd. The observed total line flux is ∼3 × 10−4 Lwd (Manser et al. 2016), or we require n ∼ 3.5. In comparison, using the above midplane radial optical depth (τ ∼ 1600) and assuming that the disk is in vertical hydrostatic equilibrium, we find that the disk can capture these photons up to a height multiple of n ∼ 3.8. In other words, the observed Ca ii line fluxes can be explained by our inferred disk density.
5.2.The Three Radii
We aim to draw clues on the progenitor by considering the three radii we inferred for the disk, ain, abreak, and aout. With our values of p1 and p2, most of the disk mass lies closely around . This suggests that abreak may be a special radius for the progenitor body. Gas deposited here may then viscously spread both outward and inward, forming the extended disk. In particular, Metzger et al. (2012) and Rafikov (2016) showed that the surface density of an isothermal accretion disk should scale as r−2, similar to our value of p2 = −1.9.
We then consider two physical radii, one for dust sublimation and the other for tidal disruption.
Using the most updated stellar parameters from Koester et al. (2014), M* = 0.705M⊙, R* = 0.0117 R⊙, Teff = 20, 713K, a blackbody at distance r is heated to a temperature
where 1300 K is roughly the sublimation temperature of silicate grains under our disk midplane density, ρ ∼ 3 × 10−10 g cm−3 (Pollack et al. 1994). 8 This will place the sublimation radius at around 1.5 R⊙, much beyond our inferred break radius (), but close to our inferred outer edge for the gas disk,
. It is therefore likely that the gas disk is truncated near the sublimation boundary. There may well be a dust component lying beyond it (see discussions below).
On the other hand, the tidal disruption radius typically lies closer than the sublimation radius. For a body with negligible internal strength 9 and on a parabolic orbit (e ≈ 1), this is located at a pericenter distance of (Sridhar & Tremaine 1992; Watanabe & Miyama 1992)
where ρp is the body's bulk density. If the orbit is less eccentric, this radius moves outward, reaching a maximum of 1.5 R⊙ at e = 0 (Chandrasekhar 1961).
Given our observed break radius (), it seems reasonable to suggest that the progenitor enters the zone of tidal disruption, is threaded apart, and its debris then sublimates to form the observed disk (also see McDonald & Veras 2021). There is much energy dissipation during this process, so the resultant disk is expected to be more compact and less eccentric. This so-called "circularization" process has been discussed by Veras et al. (2015), O'Connor & Lai (2020), Malamud et al. (2021) and Trevascus et al. (2021). Of particular note is the Trevascus et al. (2021) work, where they investigated the formation of a gaseous disk under an eccentric, outgassing parent body. Unfortunately, current simulations can only integrate for hundreds of orbits, or in physical time, ∼days. This is too short to predict the final disk shape (eccentricity and mass), which evolves on precession timescales (decades).
Lastly, the observed line profiles clearly indicate that the gas disk has an abrupt inner cutoff at . What is the origin of this cutoff? Metzger et al. (2012) suggested that a stellar magnetic field may be able to truncate the disk, much like that around in T Tauri stars. However, we suspect a different explanation may be at work here. In the inner part of the disk, rigid precession under the stronger GR precession demands a steeper eccentricity gradient. If the disk extends closer to the white dwarf than is observed, the implied high eccentricity will be challenging for its survival—nonlinear effects may disrupt the pattern of rigid precession and cause streamlined crossing. We hypothesize that disk eccentricity, rather than stellar magnetic field, truncates our disk at the observed ain. A more detailed study is required.
5.3.A Cospatial Dust Disk?
The spectral energy distribution (SED) of J1228 shows the presence of a dusty component, with a total luminosity of Ldust ∼ 6 × 10−3 Lwd, and blackbody temperatures that range from 450 to 1700 K (Brinkworth et al. 2009). If the dust lies in a geometrically flat disk and sees the central star unobstructed, it should be illuminated to a temperature (Chiang & Goldreich1997; Jura 2003)
So the above dust temperatures translate to a range of 0.2 R⊙ to 1.2 R⊙.
But such a model would place the dust component in the same annulus as the gas disk. The eccentricity of the gas disk makes this problematic. The gas and dust components, if orbiting at different eccentricities and precessing independently (i.e., dust experiences GR but not gas pressure), will encounter each other at enormous speeds, of order a few hundred km/s. This would lead to the circularization of the gas disk if the dust mass is high enough. In fact, this unwelcomed prospect led Metzger et al. (2012) to suggest that the gas disk cannot be eccentric, a proposition now amply refuted by our analysis.
One way to resolve this is if the dust component does not lie in a flat, opaque disk, a proposal that is also supported by dust observations of GD 56 (Jura et al. 2007), WD J0846+5703, (Gentile Fusillo et al. 2021) and G29-38 (Ballering et al. 2022). If the dust grains are not in a flat disk and do not block each other's view of the star, their temperatures will then be described by Equation (16), and they can remain hot out to larger distances. The observed blackbody may then arise from a region beyond 0.9 R⊙, largely avoiding the most eccentric part of the gas disk. These grains can lie even further away if they are smaller than the wavelengths of their own thermal radiation and are thus superheated.
Such a situation (free-floating grains) can arise if the grains are short lived and have not yet undergone collisional flattening (into a thin disk). Because of the proximity of the gas disk to the sublimation radius, these grains may be in condensation/sublimation equilibrium with the gas disk and are transiently formed and destroyed. In this case, one can obtain a lower limit to the dust mass by assuming that the observed dust luminosity is produced by grains of size s, bulk density ρbulk and temperature T,
In other words, only a minute amount of dust is needed to reproduce the observed SED. Most of the progenitor mass may instead lie in the gas disk.
In this scenario, the dust component (and its Poynting–Robertson drag) will be irrelevant to the evolution and accretion of the gas disk, differing from the proposal by Rafikov (2011) and Metzger et al. (2012).
5.4.The Progenitor
Here, we remark on how the J1228 disk informs on the origin of white dwarf pollution.
The observed gas disk is markedly eccentric. Barring the possibility that the eccentricity is excited after formation (see, e.g., the proposal from Miranda & Rafikov 2018), this points to a very eccentric orbit for the progenitor. This is expected in the hypothesis of tidal disruption.
We can infer the original pericenter approach by assuming that the orbital angular momentum () is largely conserved when the tidal debris is circularized into the observed disk. This yields a pericenter distance, rp = aorig(1 − eorig)
, for eorig ≈ 1 and ebreak ∼ 0.26 at abreak. This is a couple of times smaller than the Roche radius (Equation (17)).
Who can place the progenitor on such an odd orbit, with an implausibly small pericenter approach? The most likely scenario left is planetary perturbations (see review by Veras 2021). One often invoked mechanism is the so-called Kozai–Lidov oscillations (Kozai 1962; Lidov 1962), wherein secular torque from a strongly misaligned planetary perturber kneads the progenitor orbit into an almost needle-like radial orbit. 10 However, such a process could be easily suppressed by GR precession (e.g., Eggleton & Kiseleva-Eggleton 2001). To avoid suppression, the secular perturbations may need to act in a cohort with mean-motion interactions and/or close encounters (the so-called planet scattering; see, e.g., Nagasawa et al. 2008). This then requires that the orbital separation of the progenitor be within a factor of a few from that of the planet.
Where could the putative planet lie? We invoke the following estimate to argue that it must not be much further than an astronomical unit. Such a new constraint, unseen before in literature, is enabled by the very short diffusion time of the gas disk. Things have to happen in a hurry.
After the progenitor is shredded apart, its debris flies on different orbits, but ones that have a long axis well aligned with the original one. They then start the process of circularization. In principle, differential GR precession can scramble the orbital orientations and produce crossing orbits, which then circularize as mutual collisions dissipate the kinetic energy. However, if the progenitor's semimajor axis is a, and its pericentre distance rp ∼ 0.47 R⊙, the GR precession time is of order ∼2.5 × 105(a/au)3/2 yr (Equation (6)). This is so long that it seems hard to circularize the debris disk before the gaseous component diffuses away (Equation (14), ∼200 yr). In other words, GR alone will not circularize the debris disk fast enough to allow much gas accumulation.
On the other hand, if a planetary perturber lies within an astronomical unit, it is able to quickly scramble the debris orbits. We come to this conclusion as follows. To reduce the orbital energy by order unity (to achieve an order unity shrinkage of the semimajor axis), the debris orbits need to be misaligned by of order unity. This can be induced either by secular perturbations from the planet or by successive scatterings. Both take about the same time for a particle that crosses the planet's orbit. So we will simplify the argument by setting the scrambling time to be the secular time,
This is consistent with results from Li et al. (2021), where a Neptune-massed planet at 10 au takes ∼106 yr to scramble the orbits.
To replenish the gas disk within its short diffusion time (200 yr), we therefore require the scrambling time to be shorter than the diffusion time, or
where MJup is the mass of Jupiter.
Such a close-in planet is at tension with planet survival around an evolved star. Previous calculations (Villaver & Livio 2009; Kunitomo et al. 2011; Mustill & Villaver 2012; Adams & Bloch 2013; Nordhaus & Spiegel 2013; Villaver et al. 2014; Madappatt et al. 2016; Ronco et al. 2020) have shown that as a star evolves through the giant phases and then loses its envelope, the closest distance one expects a planet (or planetesimals) to survive should lie outside a few astronomical unit—reaching ∼7 au for an initial stellar mass of 3M⊙ (Figure 10 of Mustill & Villaver 2012).
It is unclear yet how to resolve this tension. If the perturber lies much further away, one has to come up with an alternative way to rapidly circularize the debris disk. Past proposals have included radiative drag (Veras et al. 2015) or collisional drag with a pre-existing disk (Malamud et al. 2021). However, these are not effective unless the debris has been ground down to microscopic sizes, and even then, it is hard to satisfy the 200 yr constraint.
We now turn to the source population that gives rise to the progenitor.
If the gas disk disappears in a few diffusion times (∼200 yr), that means the disruption event must have occurred only recently, within the past few centuries. We then estimate how common such an event occurs during the lifetime of J1228. White dwarfs that are of a similar temperature as J1228 (so that they can sublimate grains and maintain a gaseous accretion disk) have a space density of ∼4 × 10−5 pc−3 mag−1 (Leggett et al. 1998). So within the distance of J1228 (∼127 pc), there are about 300 such white dwarfs. Among these, about two dozen gas disks have been reported (Gänsicke et al. 2006, 2007, 2008; Farihi et al. 2012; Gentile Fusillo et al. 2021). So the occurrence rate of gas disks among hot white dwarfs is ∼10%, not dissimilar to the occurrence rate of dust disks among all white dwarfs. 11 So the total mass of (disrupted) progenitors (all assumed to be Rp = 50 km) over the cooling age of J1228 (∼107 yr) 12 is,
This mass is similar to the total mass of our own asteroid belt (∼2 × 1024 g). However, the disrupted bodies are those that are able to reach the Roche radius and so likely only constitute a small fraction of the total source disk. Moreover, evidence on older white dwarfs suggests that metal pollution continues well past the cooling age of J1228. Both these arguments imply that, for us to witness a rapidly draining gas disk around J1228, there must exist planetesimal disks around these white dwarfs that are much more massive than what one expects from their main-sequence counterparts. This, a massive belt at au distance, is again in tension with our current knowledge. While a planetesimal belt at the Kuiper Belt distance can remain massive for longer due to a slower collision timescale, a belt at the Asteroid belt distance should grind down quickly during the main-sequence lifetime of a star, even if it started massive (see, e.g., Wyatt et al. 2007).
Lastly, we comment on the report of a planetesimal embedded in the gas disk. Manser et al. (2019) discovered periodic variations in the line profiles of J1228, with a period of 123.4 min (corresponding to an orbital semimajor axis of 0.73 R⊙). They interpreted these disturbances as due to an outgassing planetesimal on an eccentric orbit that overlaps with the gas disk. If such an object does exist and is massive (radius R ≥ 50 km), our above arguments on the disk origin can fail because the disk can now be continuously fed and can survive much longer than the diffusion time. 13 However, there are multiple reasons to suspect the planetesimal interpretation. First, as remarked in Manser et al. (2019), the survival of such a body well inside the nominal Roche radius is problematic: the body has to have an unphysically high bulk density (ρp ≥ 40 g cm−3 for circular orbit; higher value needed for eccentric orbit; see Equation (17)); or the body has to have a high internal strength (see Equation 4.78 of Murray & Dermott 1999), comparable to that of a monolithic rock. Moreover, if the planetesimal precesses at a different rate than the gas disk (which is physically likely), there will be continuous streamline crossings between the existing gas disk and the newly sublimated gas, resulting in a disk that is nearly circular and much more compact than the parent body. Lastly, it is hard to imagine how a massive planetesimal can evolve from a nearly parabolic orbit to its current close orbit, short of any sensible dissipation mechanism.
6.Conclusion
We undertake a detailed modeling of the gaseous disk around the white dwarf SDSS J1228+1040. We find that the disk has a surface density profile that peaks around , and an eccentricity profile that decreases outward. The latter, we show, uncannily reproduces the theoretical profile of a disk that precesses rigidly under the combined forces of general relativity and gas pressure. In other words, the observed disk is in an eccentric eigenstate. This explains why the disk can be eccentric yet long lived. As we expect the eccentricity to be dissipated in the viscous timescale (∼200 yr), the current disk should have formed fairly recently.
Based on the high accretion rate of the white dwarf (current estimate ), we infer a mass of ∼1021 g for the gaseous disk. Such a mass estimate is also consistent with the emission measures in the Ca ii triplets. The young age of the current disk then implies that it still contains most of the source mass, pegging the progenitor at a size of R ∼ 50 km.
Given the eccentricity of the gas disk, there is unlikely to be a massive dusty debris disk that is cospatial. Rather, we suggest that the observed dust emission may arise from small amounts of grains that are in condensation/sublimation equilibrium with the gas disk.
The eccentricity of the remnant disk supports the hypothesis of tidal disruption. The progenitor body, after being excited to a nearly parabolic orbit and reaching a pericenter distance of , was tidally disrupted by the white dwarf. The debris then circularizes into the current moderately eccentric and compact disk.
The short diffusion time of the gas disk allows us to set unique and stringent constraints on the origin of the system. By insisting that the tidal debris circularizes into the current size in a time shorter than the gas disk diffusion time, we find that both the planetesimal and its planetary perturber likely orbit around the white dwarf within an astronomical unit. This is in tension with our current knowledge of planet survival around white dwarf progenitors.
Moreover, to account for the observed rate of gaseous disks among similarly hot white dwarfs, we estimate that the source disk needs to contain at least the mass of our asteroid belt and likely many orders of magnitude more in reality. Such a massive disk does not seem to exist around the main-sequence counterparts of these white dwarfs.
Our study in this work is preliminary in nature. We have not investigated in detail the temperature structure and the emission mechanism of the metallic disk. Emission line diagnostics may be used to constrain the disk inclination, which may lead to further insights. We also fall short of analyzing the "circularization" process after tidal disruption. This latter seems a promising route to infer the nature of the progenitor and the architecture of the planetary system around J1228.
Such careful studies are clearly warranted. J1228 is likely not unusual. First, many white dwarfs with gaseous disks show variable emissions, indicating eccentric, precessing disks (Gänsicke et al. 2008; Melis et al. 2010; Wilson et al. 2014; Cauley et al. 2018; Manser et al. 2021). Similar dynamics as we reveal here for J1228 may be in play in all of these disks. Second, while J1228 is hot and can sublimate rocks at around 1 R⊙, cooler white dwarfs will only harbor fully dusty disks. Such a disk reveals no information on its kinematics and surface density. J1228 offers us a lucky window into these otherwise obscure disks and may well be the Rosetta stone to decipher the mystery of white dwarf pollution.
Acknowledgments
We acknowledge NSERC for funding. We also thank C. Manser for providing the line data, and Renu Malhotra and Elliot Lynch for discussions. Lastly, we thank the anonymous referee for a careful critique of the paper.