NASA/JPL-Caltech http://hershy314.deviantart.com/ M. Kornmesser Comparing asteroids, comets, moons & planets as WD pollution progenitors Dimitri Veras (University of Warwick) ESO / http://becuo.com
Start with observables Number of confirmed WD planets Mullally et al. (2008) Hogan et al. (2009) 0! Debes et al. (2011) Faedi et al. (2011) Steele et al. (2011) Number of WD dust discs Zuckerman & Becklin (1987); Becklin et al. (2005) 30 Kilic et al. (2005) Reach et al. (2005) Farihi et al. (2009) Number of WD gas discs! Gänsicke et al. (2006) Gänsicke et al. (2007) Gänsicke et al. (2008) 7 Debes et al. (2012) Brinkworth et al. (2012) Melis et al. (2012)
Location of WD discs See Debes talk on Wed! SUN (FOR SCALE) DISC WD VERY COMPACT!
Metal pollution in white dwarf atmospheres Widespread Diverse Sometimes Earth-like Sometimes water-rich Zuckerman et al. (2003, 2010), Koester et al. (2014) See Zuckerman talk on Fri! e.g. Gänsicke et al. (2012) See Gänsicke talk on Thurs! e.g. Klein et al. (2010) Farihi et al. (2013), See Raddi talk on Fri!
Goal: Reproduce Observables Accretion occurs early and late in DA WDs Koester, Gänsicke, Farihi (2014)
Goal: Reproduce Observables Accretion of up to 10 25 g over 1 Myr in DB WDs Girven et al. (2012) DBZ DAZ
What causes WD pollution? NASA/JPL-Caltech Asteroids NASA/JPL-Caltech Comets Planets Gemini Observatory Jon Lomberg Compositionally Dynamically
What causes WD pollution? NASA/JPL-Caltech Asteroids Graham et al. (1990) Jura (2003, 2008) Bonsor et al. (2011) Debes et al. (2012) Bear & Soker (2013) Frewen & Hansen (2014) Veras et al. (2014a,b) NASA/JPL-Caltech Comets Planets Gemini Observatory Jon Lomberg
Asteroid Challenges Unknown accompanying planet(s) One planet + One belt Bonsor, Mustill & Wyatt (2011) Kuiper belt scattered inward Debes, Walsh & Stark (2012) Resonant diffusion towards WD Frewen & Hansen (2014) Planet mass and eccentricity dependence
Asteroid Challenges Unknown accompanying debris Veras, Jacobson, Gänsicke (Submitted to MNRAS) 100m - 10km asymmetric asteroids
Asteroid Challenges Unknown accompanying debris Veras, Jacobson, Gänsicke (Submitted to MNRAS)
Asteroid Challenges Need self-consistent size Wyatt, Farihi, Pringle & Bonsor (2014) Mass power law exponent q t 2 sink + t2 disc (yr)
Asteroid Challenges Depleted during main-sequence Debes, Walsh & Stark (2012) Number of Exo-belts Median mass = 820MSS-belt log (M/M SS belt )
Asteroid Challenges Need to numerically resolve very close pericentre passages Veras, Leinhardt, Bonsor, Gänsicke (Submitted to MNRAS) Disruption Spreading
Asteroid Challenges Need to numerically resolve very close pericentre passages Veras, Leinhardt, Bonsor, Gänsicke (Submitted to MNRAS) time / years Formation timescale of eccentric rings Breakup distance / RWD
What causes WD pollution? NASA/JPL-Caltech Asteroids Alcock et al. (1986) Parriott & Alcock (1998) Jura (2011) Stone et al. (2014) Veras et al. (2014c) NASA/JPL-Caltech Comets Planets Gemini Observatory Jon Lomberg
Comet Challenges Need to incorporate Galactic tides, stellar flybys and mass loss Main sequence TMSyr orbital period / yr 10 9 10 8 10 7 10 6 10 5 10 4 10 3 10 2 10 1 10 0 Regions of Motion for 0.02 Veras, Evans, Wyatt, Tout (2014) [ Great Escape III ] Beyond Hill Axes x z Nonadiabatic PostMS 0.1 0.2 0.5 1 2 3 5 8 12 20 30 50 Rkpc TPAGB Adiabatic Nonadiabatic Galactic MS Orbit Crossin g with Clo sest Stellar Flyby 8M 1M M 0 1M M 0 8M Bulge Disc Halo Distance to Galactic Centre / kpc
Comet Challenges Exo-Oort clouds not massive enough Veras, Shannon, Gänsicke (Submitted to MNRAS)
Comet Challenges Exo-Oort clouds not massive enough Veras, Shannon, Gänsicke (Submitted to MNRAS) GB Start
Distance / au Comet Challenges Composition-dependent dynamics Stone, Metzger & Loeb (Submitted to MNRAS) Sublimation distances comet comet ice comet rock WD cooling age / yr See Stone poster for more details! (#56) See Shannon poster for Oort cloud rockets! (#6)
What causes WD pollution? NASA/JPL-Caltech Asteroids Debes & Sigurdsson (2002) Veras et al. (2013) Voyatzis et al. (2013) Mustill et al. (2014) NASA/JPL-Caltech Comets Planets Gemini Observatory Jon Lomberg Next talk!
What about moons? Contains more mass than asteroids Regular delivery to WD dynamically unlikely (?) Hardest to model numerically (and perhaps analytically) See Payne poster! (#52)
Conclusions Origin of WD debris still mystery Dynamically, exo-oort cloud comets cannot be primary Total mass and delivery timescales still issues for asteroids Post-main-sequence planetary architectures largely unexplored
Semimajor axis in astronomical units 5 4 3 2 1 Tidal engulfment distance for 1 planet Mustill & Villaver (2012) 0.0 0.2 0.4 0.6 Time (Myr) Giant Planets Survive Planets Crash Stellar Radius Pulses
0.01 0.1 1 10 100 1000 10 7 10 5 0.001 0.1 10 Semimajor axis in astronomical units Minimum mass in Jupiter masses Surviving the giant branch phases ENGULFED SURVIVED Exoplanets Solar System Planets Asteroid Belt Kuiper Belt Sedna Will planets be swallowed? Giant Mass in Jupiter masses ME MA V E
Planet motion: Mass- loss 2-body problem Giant Gyldén (1884) ; Mestschersky (1893) ; Jeans (1924)
Regimes of motion Veras, Wyatt, Mustill, Bonsor, Eldridge (2011) Giant Orbital Period << Mass Loss Timescale ADIABATIC REGIME Semimajor axis increases Eccentricity is constant Orbital Period >> Mass Loss Timescale RUNAWAY REGIME Semimajor axis increases more Eccentricity no longer constant
Adiabatic regime Veras, Wyatt, Mustill, Bonsor, Eldridge (2011) Giant Planet Planet AU 100 True Anomaly (deg) fdegrees 270 50 Star Star 350 300 250 200 150 100 50 AU Star 180 50 90 0 10 20 30 40 50 Time/1000 yr Time (1000 yr)
Runaway regime Veras, Wyatt, Mustill, Bonsor, Eldridge (2011) Giant Planet Planet AU 200 150 100 True Anomaly (deg) fdegrees 360 270 AU Star Star Star 1000 800 600 400 200 50 50 180 90 100 0 0 3 6 9 12 15 Time/1000 yr Time (1000 yr)
The Solar System s Critical semimajor axis in 10 3 au a 010 3 AU 10.0 5.0 4.0 3.0 2.0 1.0 Stability Boundary The Solar System s Stability Boundary DANGER OF EJECTION BOUNDED 0.2 1 0.3 2 0.4 3 0.5 4 0.6 5 0.7 6 0.8 7 Thermally Pulsing AGB Veras & Wyatt (2012) Model Η Number Early AGB RGB Giant 10.0 5.0 4.0 3.0 2.0 1.0
Anisotropic mass loss Giant Veras, Hadjidemetriou, Tout (2013)
1.15 1.25 1.35 1.45 1.55 2 4 6 8 10 2:1 5:4 4:3 3:2 5:3 Initial Semimajor Axis Ratio Lo g 10 tinstyr Instability Times for M t 0 6M MS WD Hill Stability Boundary Veras, Mustill, Bonsor, Wyatt (2013) Full-lifetime simulations
Observation Motivation Disc shape change Gänsicke, Koester, Marsh, Southworth, Rebassa-Mansergas (2008) 2004 2008 WD
Link with observations WD Veras, Mustill, Bonsor, Wyatt (2013) EjectionHyper Orbit PlanPlan Hit WD Collision 25 Number of of Unstable unstable Systems WD systems 20 15 10 M 0 5M SIMULATED 5 5 6 7 8 9 Log 10 White Dwarf Cooling Ageyr Log 10 [White Dwarf Cooling Age] /yr An Observed Polluted WD Sample OBSERVED SAMPLE Number of WD polluted systems Number of WD Polluted Systems 15 10 5 8 9 10 Log 10 White Dwarf Cooling Ageyr 8 9 10
Hill Stability Marchal & Bozis (1982) Milani & Nobili (1983) Gladman (1993) Georgakarakos (2008) MS Donnison (2011) Veras, Mustill, Bonsor, Wyatt (2013)
Lagrange Stability Barnes & Greenberg (2006, 2007) Raymond et al. (2009) Veras & Mustill (2013) MS log 10 x 5.2 µ M J /M 0.18
Many planets Chambers et al. (1996) Chatterjee et al. (2008) Smith & Lissauer (2009) MS
How big? Giant Veras, Mustill, Bonsor, Wyatt (2013) Sun
How much mass lost? Giant Veras, Mustill, Bonsor, Wyatt (2013) All stars lose over half of their mass
0 2 4 6 8 Mutual Hill Radii Multi-planet instability 8 Debes & Sigurdsson (2002) Giant Log t (orbits) 6 4 2 Without Mass Loss With Constant Mass Loss
Multi-planet instability from strong mass loss Giant Voyatzis, Hadjidemetriou, Veras, Varvoglis (2013) Mass Loss Stopped Planet 1 Planet 2 Time (kyr)
Multi-planet instability MS Mustill, Veras, Villaver (2014) Giant WD
Initial Semimajor Axis Ratio Veras, Mustill, Bonsor, Wyatt (2013) 30 Survivor Orbits for M t 0 5M Semimajor WD 25 Axes DistanceAU 20 15 10 5 0 Pericentres Max AGB Radius 1.25 1.30 1.35 1.40 1.45 1.50 1.55
1 planet, 1 star, 0 discs Veras, Evans, Wyatt, Tout (2014)
Veras, Shannon, Gänsicke (Submitted to MNRAS 2014) Exo-Oort clouds Extends from 10 4 au 10 5 au Radial distributions a 1,a 1.5,a 2 Inclinations sin i Eccentricities random Orbital angles random # of comets per sim: 5000
Veras, Shannon, Gänsicke (Submitted to MNRAS 2014) Comet escape from close flyby
Observation Motivation Bowler et al. (2010) Absence of planets at low semimajor axes Giant Star Stellar Mass Planets
Observation Motivation Gänsicke, Marsh, Southworth, Rebassa-Mansergas (2006) Gaseous discs Points =Ca II emission line profiles
Observation Motivation Koester (2013) Typical WDs DB He lines DA H lines
Atmospheric sinking time / yr Observation Motivation Wyatt, Farihi, Pringle, Bonsor (2014) Sinking timescales 10 6 10 4 10 2 He-dominated DB H-dominated DA 10 0 10 2 10 8 10 9 Time after becoming white dwarf / yr
Observation Motivation Gänsicke et al. (2012) Metal pollution is diverse
Observation Motivation Klein, Jura, Koester, Zuckerman, Melis (2010) Some pollution is Earth-like White dwarf GD 40
Observation Motivation Farihi, Gänsicke, Koester (2013) Some pollution includes water White Dwarf
Age/yr Fate of Earth Schröder & Connon Smith (2008) Giant Earth Orbit Mass of Sun log 10 (d/r ) Red Giant Branch Asymptotic Giant Branch
Asteroid Challenges Circularising eccentric rings Veras, Leinhardt, Gänsicke (In Prep) Sublimation + PR-Drag 10 Semimajor axis decrease from radiation 1 aau 0.1 0.01 0.001 10 4 R 10 4 m R 10 3 m R 10 2 m R 10 1 m R 10 0 m R 10 2 m R 10 4 m R 10 6 m 10 9 10 7 10 5 10 3 10 1 10 1 10 3 tmyr