BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Dimitri Veras

Dimitri Veras appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
astro-ph.EPastro-ph.SRastro-ph.GAastro-ph.IMastro-phmath-phmath.MPphysics.pop-phApplicationsComputationgr-qcMethodologyphysics.class-phphysics.ed-ph

Catalog footprint

What is connected

83works
14topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

83 published item(s)

preprint2022arXiv

A white dwarf accreting planetary material determined from X-ray observations

The atmospheres of a large proportion of white dwarf stars are polluted by heavy elements that are expected to sink out of visible layers on short timescales. This has been interpreted as a signature of ongoing accretion of debris from asteroids, comets, and giant planets. This scenario is supported by the detection of debris discs and transits of planetary fragments around some white dwarfs. However, photospheric metals are only indirect evidence for ongoing accretion, and the inferred accretion rates and parent body compositions heavily depend on models of diffusion and mixing processes within the white dwarf atmosphere. Here we report a 4.4$σ$ detection of X-rays from a polluted white dwarf, G29$-$38, using a 106 ks exposure with the Chandra X-ray Observatory, demonstrating directly that the star is currently accreting. From the measured X-ray luminosity, we find an instantaneous accretion rate of $\dot{M_{\rm X}}=1.63^{+1.29}_{-0.40}\times 10^{9}\mathrm{\,g\,s^{-1}}$. This is the first direct measurement of the accretion rate onto the white dwarf, which is independent of stellar atmosphere models. This rate exceeds estimates based on past studies of the photospheric abundances by more than a factor two, and implies that convective overshoot has to be accounted for in modelling the spectra of debris-accreting white dwarfs. We measure a low plasma temperature of $kT=0.5\pm0.2\,\mathrm{keV}$, corroborating the predicted bombardment solution for white dwarfs accreting at low accretion rates. Offering a new method for studying evolved planetary systems, these observations provide the opportunity to independently measure the instantaneous accretion rate of planetary material, and therefore investigate the timescale of accretion onto white dwarfs, and the evolution and replenishment of debris disks.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Birth cluster simulations of planetary systems with multiple super-Earths: initial conditions for white dwarf pollution drivers

Previous investigations have revealed that eccentric super-Earths represent a class of planets which are particularly effective at transporting minor bodies towards white dwarfs and subsequently polluting their atmospheres with observable chemical signatures. However, the lack of discoveries of these planets beyond a few astronomical units from their host stars prompts a better understanding of their orbital architectures from their nascent birth cluster. Here, we perform stellar cluster simulations of 3-planet and 7-planet systems containing super-Earths on initially circular, coplanar orbits. We adopt the typical stellar masses of main-sequence progenitors of white dwarfs ($1.5\,\mathrm{M}_{\odot}$-$2.5\,\mathrm{M}_{\odot}$) as host stars and include 8,000 main-sequence stars following a Kroupa initial mass function in our clusters. Our results reveal that about 30 per cent of the simulated planets generate eccentricities of at least 0.1 by the time of cluster dissolution, which would aid white dwarf pollution. We provide our output parameters to the community for potential use as initial conditions for subsequent evolution simulations.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Chaos over Order: Mapping 3D Rotation of Triaxial Asteroids and Minor Planets

Celestial bodies approximated with rigid triaxial ellipsoids in a two-body system can rotate chaotically due to the time-varying gravitational torque from the central mass. At small orbital eccentricity values, rotation is short-term orderly and predictable within the commensurate spin-orbit resonances, while at eccentricity approaching unity, chaos completely takes over. Here, we present the full 3D rotational equations of motion around all three principle axes for triaxial minor planets and two independent methods of numerical solution based on Euler rotations and quaternion algebra. The domains of chaotic rotation are numerically investigated over the entire range of eccentricity with a combination of trial integrations of Euler's equations of motion and the GALI($k$) method. We quantify the dependence of the order--chaos boundaries on shape by changing a prolateness parameter, and find that the main 1:1 spin-orbit resonance disappears for specific moderately prolate shapes already at eccentricities as low as 0.3. The island of short-term stability around the main 1:1 resonance shrinks with increasing eccentricity at a fixed low degree of prolateness and completely vanishes at approximately 0.8. This island is also encroached by chaos on longer time scales indicating longer Lyapunov exponents. Trajectories in the close vicinity of the 3:2 spin-orbit resonance become chaotic at smaller eccentricities, but separated enclaves of orderly rotation emerge at eccentricities as high as 0.8. Initial perturbations of rotational velocity in latitude away from the exact equilibrium result in a spectrum of free libration, nutation, and polar wander, which is not well matched by the linearized analysis omitting the inertial terms.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Mathematical encoding within multi-resonant planetary systems as SETI beacons

How might an advanced alien civilization manipulate the orbits within a planetary system to create a durable signpost that communicates its existence? While it is still debated whether such a purposeful advertisement would be prudent and wise, we propose that mean-motion resonances between neighboring planets -- with orbital periods that form integer ratios -- could in principle be used to encode simple sequences that one would not expect to form in nature. In this Letter we build four multi-resonant planetary systems and test their long-term orbital stability. The four systems each contain 6 or 7 planets and consist of: (i) consecutive integers from 1 to 6; (ii) prime numbers from 2 to 11; (iii) the Fibonacci sequence from 1 to 13; and (iv) the Lazy Caterer sequence from 1 to 16. We built each system using N-body simulations with artificial migration forces. We evaluated the stability of each system over the full 10 Gyr integration of the Sun's main sequence phase. We then tested the stability of these systems for an additional 10 Gyr, during and after post-main sequence evolution of the central stars (assumed to be Sun-like) to their final, white dwarf phase. The only system that was destabilized was the consecutive integer sequence (system i). The other three sequences therefore represent potential SETI beacons.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

TESS Giants Transiting Giants III: An eccentric warm Jupiter supports a period-eccentricity relation for giant planets transiting evolved stars

The fate of planets around rapidly evolving stars is not well understood. Previous studies have suggested that relative to the main sequence population, planets transiting evolved stars ($P$ $<$ 100 d) tend to have more eccentric orbits. Here we present the discovery of TOI-4582 b, a 0.94 $\pm$ 0.12 R$_\mathrm{J}$, 0.53 $\pm$ 0.05 M$_\mathrm{J}$ planet orbiting an intermediate-mass subgiant star every 31.034 days. We find that this planet is also on a significantly eccentric orbit ($e$ = 0.51 $\pm$ 0.05). We then compare the population of planets found transiting evolved (log$g$ $<$ 3.8) stars to the population of planets transiting main sequence stars. We find that the rate at which median orbital eccentricity grows with period is significantly higher for evolved star systems than for otherwise similar main sequence systems, particularly for systems with only one planet detected. In general, we observe that mean planet eccentricity $<e>$ = $a$ + $b$log$_{10}$($P$) for the evolved population with a single transiting planet where $a$ = (-0.18 $\pm$ 0.08) and $b$ = (0.38 $\pm$ 0.06), significantly distinct from the main sequence planetary system population. This trend is seen even after controlling for stellar mass and metallicity. These systems do not appear to represent a steady evolution pathway from eccentric, long-period planetary orbits to circular, short period orbits, as orbital model comparisons suggest inspiral timescales are uncorrelated with orbital separation or eccentricity. Characterization of additional evolved planetary systems will distinguish effects of stellar evolution from those of stellar mass and composition.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Rapid destruction of planetary debris around white dwarfs through aeolian erosion

The discovery of numerous debris disks around white dwarfs (WDs), gave rise to extensive study of such disks and their role in polluting WDs, but the formation and evolution of these disks is not yet well understood. Here we study the role of aeolian (wind) erosion in the evolution of solids in WD debris disks. Aeolian erosion is a destructive process that plays a key role in shaping the properties and size-distribution of planetesimals, boulders and pebbles in gaseous protoplanetary disks. Our analysis of aeolian erosion in WD debris disks shows it can also play an important role in these environments. We study the effects of aeolian erosion under different conditions of the disk, and its erosive effect on planetesimals and boulders of different sizes. We find that solid bodies smaller than $\sim 5 \rm{km}$ will be eroded within the short disk lifetime. We compare the role of aeolian erosion in respect to other destructive processes such as collisional fragmentation and thermal ablation. We find that aeolian erosion is the dominant destructive process for objects with radius $\lesssim 10^3 \rm{cm}$ and at distances $\lesssim 0.6 R_\odot$ from the WD. Thereby, aeolian erosion constitutes the main destructive pathway linking fragmentational collisions operating on large objects with sublimation of the smallest objects and Poynting-Robertson drag, which leads to the accretion of the smallest particles onto the photosphere of WDs, and the production of polluted WDs.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

The Pedagogical Representation of Mass Functions with LEGO and their Origin

We promote the teaching of mass functions as an integral part of an interdisciplinary science education. Mass functions characterize the frequency distributions of objects with different masses on all cosmic scales. We intend to enhance experiential learning of this concept with a creative LEGO brick experiment for a diverse student audience. To our surprise, the LEGO mass function is not only qualitatively but also quantitatively comparable to mass functions found across the Universe. We also discuss the relation between gravitation and mass distributions as a possible explanation for the continuity of the universal mass function.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

A remnant planetary core in the hot-Neptune desert

The interiors of giant planets remain poorly understood. Even for the planets in the Solar System, difficulties in observation lead to large uncertainties in the properties of planetary cores. Exoplanets that have undergone rare evolutionary processes provide a route to understanding planetary interiors. Planets found in and near the typically barren hot-Neptune 'desert' (a region in mass-radius space that contains few planets) have proved to be particularly valuable in this regard. These planets include HD149026b, which is thought to have an unusually massive core, and recent discoveries such as LTT9779b and NGTS-4b, on which photoevaporation has removed a substantial part of their outer atmospheres. Here we report observations of the planet TOI-849b, which has a radius smaller than Neptune's but an anomalously large mass of $39.1^{+2.7}_{-2.6}$ Earth masses and a density of $5.2^{+0.7}_{-0.8}$ grams per cubic centimetre, similar to Earth's. Interior structure models suggest that any gaseous envelope of pure hydrogen and helium consists of no more than $3.9^{+0.8}_{-0.9}$ per cent of the total planetary mass. The planet could have been a gas giant before undergoing extreme mass loss via thermal self-disruption or giant planet collisions, or it could have avoided substantial gas accretion, perhaps through gap opening or late formation. Although photoevaporation rates cannot account for the mass loss required to reduce a Jupiter-like gas giant, they can remove a small (a few Earth masses) hydrogen and helium envelope on timescales of several billion years, implying that any remaining atmosphere on TOI-849b is likely to be enriched by water or other volatiles from the planetary interior. We conclude that TOI-849b is the remnant core of a giant planet.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Constraining planet formation around 6$M_{\odot}$-8$M_{\odot}$ stars

Identifying planets around O-type and B-type stars is inherently difficult; the most massive known planet host has a mass of only about $3M_{\odot}$. However, planetary systems which survive the transformation of their host stars into white dwarfs can be detected via photospheric trace metals, circumstellar dusty and gaseous discs, and transits of planetary debris crossing our line-of-sight. These signatures offer the potential to explore the efficiency of planet formation for host stars with masses up to the core-collapse boundary at $\approx 8M_{\odot}$, a mass regime rarely investigated in planet formation theory. Here, we establish limits on where both major and minor planets must reside around $\approx 6M_{\odot}-8M_{\odot}$ stars in order to survive into the white dwarf phase. For this mass range, we find that intact terrestrial or giant planets need to leave the main sequence beyond approximate minimum star-planet separations of respectively about 3 and 6 au. In these systems, rubble pile minor planets of radii 10, 1.0, and 0.1 km would have been shorn apart by giant branch radiative YORP spin-up if they formed and remained within, respectively, tens, hundreds and thousands of au. These boundary values would help distinguish the nature of the progenitor of metal-pollution in white dwarf atmospheres. We find that planet formation around the highest mass white dwarf progenitors may be feasible, and hence encourage both dedicated planet formation investigations for these systems and spectroscopic analyses of the highest mass white dwarfs.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Constraining the origin of the planetary debris surrounding ZTF J0139+5245 through rotational fission of a triaxial asteroid

White dwarfs containing orbiting planetesimals or their debris represent crucial benchmarks by which theoretical investigations of post-main-sequence planetary systems may be calibrated. The photometric transit signatures of likely planetary debris in the ZTF J0139+5245 white dwarf system has an orbital period of about 110 days. An asteroid which breaks up to produce this debris may spin itself to destruction through repeated close encounters with the star without entering its Roche radius and without influence from the white dwarf's luminosity. Here, we place coupled constraints on the orbital pericentre ($q$) and the ratio ($β$) of the middle to longest semiaxes of a triaxial asteroid which disrupts outside of this white dwarf's Roche radius ($r_{\rm Roche}$) soon after attaining its 110-day orbit. We find that disruption within tens of years is likely when $β\lesssim 0.6$ and $q\approx 1.0-2.0r_{\rm Roche}$, and when $β\lesssim 0.2$ out to $q\approx 2.5r_{\rm Roche}$. Analysing the longer-timescale disruption of triaxial asteroids around ZTF J0139+5245 is desirable but may require either an analytical approach relying on ergodic theory or novel numerical techniques.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

DebrisWatch I: A survey of faint geosynchronous debris

Recent anomalies exhibited by satellites and rocket bodies have highlighted that a population of faint debris exists at geosynchronous (GEO) altitudes, where there are no natural removal mechanisms. Despite previous optical surveys probing to around 10-20 cm in size, regular monitoring of faint sources at GEO is challenging, thus our knowledge remains sparse. It is essential that we continue to explore the faint debris population using large telescopes to better understand the risk posed to active GEO satellites. To this end, we present photometric results from a survey of the GEO region carried out with the 2.54 m Isaac Newton Telescope in La Palma, Canary Islands. We probe to 21st visual magnitude (around 10 cm, assuming Lambertian spheres with an albedo of 0.1), uncovering 129 orbital tracks with GEO-like motion across the eight nights of dark-grey time comprising the survey. The faint end of our brightness distribution continues to rise until the sensitivity limit of the sensor is reached, suggesting that the modal brightness could be even fainter. We uncover a number of faint, uncatalogued objects that show photometric signatures of rapid tumbling, many of which straddle the limiting magnitude of our survey over the course of a single exposure, posing a complex issue when estimating object size. This work presents the first instalment of DebrisWatch, an ongoing collaboration between the University of Warwick and the Defence Science and Technology Laboratory (UK) investigating the faint population of GEO debris.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Generating metal-polluting debris in white dwarf planetary systems from small-impact crater ejecta

Metal pollution in white dwarf photospheres originates from the accretion of some combination of planets, moons, asteroids, comets, boulders, pebbles and dust. When large bodies reside in dynamically stagnant locations -- unable themselves to pollute nor even closely approach the white dwarf -- then smaller reservoirs of impact debris may become a complementary or the primary source of metal pollutants. Here, we take a first step towards exploring this possibility by computing limits on the recoil mass that escapes the gravitational pull of the target object following a single impact onto an atmosphere-less surface. By considering vertical impacts only with the full-chain analytical prescription from Kurosawa & Takada (2019), we provide lower bounds for the ejected mass for basalt, granite, iron and water-rich target objects across the radii range 10^{0-3} km. Our use of the full-chain prescription as opposed to physical experiments or hydrocode simulations allows us to quickly sample a wide range of parameter space appropriate to white dwarf planetary systems. Our numerical results could be used in future studies to constrain freshly-generated small debris reservoirs around white dwarfs given a particular planetary system architecture, bombardment history, and impact geometries.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Gyr-timescale destruction of high-eccentricity asteroids by spin and why 2006 HY51 has been spared

Asteroids and other small celestial bodies have markedly prolate shapes, and the perturbative triaxial torques which are applied during pericenter passages in highly eccentric orbits trigger and sustain a state of chaotic rotation. Because the prograde spin rate around the principal axis of inertia is not bounded from above, it can accidentally reach the threshold value corresponding to rotational break-up. Previous investigations of this process were limited to integrations of $\sim 10^3$ orbits because of the stiff equation of motion. We present here a fast 1D simulation method to compute the evolution of this spin rate over $\sim 10^9$ orbits. We apply the method to the most eccentric solar system asteroid known, 2006 HY51 (with $e = 0.9684$), and find that for any reasonably expected shape parameters, it can never be accelerated to break-up speed. However, primordial solar system asteroids on more eccentric orbits may have already broken up from this type of rotational fission. The method also represents a promising opportunity to investigate the long-term evolution of extremely eccentric triaxial exo-asteroids ($e > 0.99$), which are thought to be common in white dwarf planetary systems

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

How Jupiters save or destroy inner Neptunes around evolved stars

In about 6 Giga years our Sun will evolve into a red giant and finally end its life as a white dwarf. This stellar metamorphosis will occur to virtually all known host stars of exo-planetary systems and is therefore crucial for their final fate. It is clear that the innermost planets will be engulfed and evaporated during the giant phase and that planets located farther out will survive. However, the destiny of planets in-between, at ~ 1 - 10 au, has not yet been investigated with a multi-planet tidal treatment. We here combine for the first time multi-planet interactions, stellar evolution, and tidal effects in an N-body code to study the evolution of a Neptune-Jupiter planetary system. We report that the fate of the Neptune-mass planet, located closer to the star than the Jupiter-mass planet, can be very different from the fate of a single Neptune. The simultaneous effects of gravitational interactions, mass loss and tides can drive the planetary system towards mean motion resonances. Crossing these resonances affects particularly the eccentricity of the Neptune and thereby also its fate, which can be engulfment, collision with the Jupiter-mass planet, ejection from the system, or survival at a larger separation.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Linking the formation and fate of exo-Kuiper belts within solar system analogues

Escalating observations of exo-minor planets and their destroyed remnants both passing through the solar system and within white dwarf planetary systems motivate an understanding of the orbital history and fate of exo-Kuiper belts and planetesimal discs. Here we explore how the structure of a 40-1000 au annulus of planetesimals orbiting inside of a solar system analogue that is itself initially embedded within a stellar cluster environment varies as the star evolves through all of its stellar phases. We attempt this computationally challenging link in four parts: (1) by performing stellar cluster simulations lasting 100 Myr, (2) by making assumptions about the subsequent quiescent 11 Gyr main-sequence evolution, (3) by performing simulations throughout the giant branch phases of evolution, and (4) by making assumptions about the belt's evolution during the white dwarf phase. Throughout these stages, we estimate the planetesimals' gravitational responses to analogues of the four solar system giant planets, as well as to collisional grinding, Galactic tides, stellar flybys, and stellar radiation. We find that the imprint of stellar cluster dynamics on the architecture of $\gtrsim 100$ km-sized exo-Kuiper belt planetesimals is retained throughout all phases of stellar evolution unless violent gravitational instabilities are triggered either (1) amongst the giant planets, or (2) due to a close ($\ll 10^3$ au) stellar flyby. In the absence of these instabilities, these minor planets simply double their semimajor axis while retaining their primordial post-cluster eccentricity and inclination distributions, with implications for the free-floating planetesimal population and metal-polluted white dwarfs.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Necroplanetology: Simulating the Tidal Disruption of Differentiated Planetary Material Orbiting WD 1145+017

The WD 1145+017 system shows irregular transit features that are consistent with the tidal disruption of differentiated asteroids with bulk densities $< 4$ g cm$^{-3}$ and bulk masses $\lesssim 10^{21}$ kg (Veras et al., 2017; arXiv:1610.06926). We use the open-source N-body code REBOUND (Rein and Liu, 2011; arXiv:1110.4876) to simulate this disruption with different internal structures: varying the core volume fraction, mantle/core density ratio, and the presence/absence of a thin low-density crust. We allow the rubble pile to partially disrupt and capture lightcurves at a specific point during the disruption at cadences comparable to those from ground-based photometry. As a proof-of-concept we show that varying these structural parameters have observationally distinguishable effects on the transit light curve as the asteroid is disrupted and compare the simulation-generated lightcurves to data from Gary et al., 2017 (arXiv:1608.00026). With the caveat that our simulations do not model the sublimation in detail or account for its effects on orbital evolution, we find that a low core fraction and low mantle/core density ratio asteroid is most consistent with the stable transit feature present for multiple weeks circa April 2016 (referred to as G6121 in Gary et al., 2017 (arXiv:1608.00026) and A1 in Hallakoun et al., 2017 (arXiv:1702.05483). Connecting tidal disruption simulations to photometry suggests characteristics for the interior structure and composition of an exoplanetary body, information that is only possible because we are observing the death of the planetary system in action. All-sky survey missions such as TESS and LSST will be able to detect other systems like WD 1145+017, creating a sample of subjects for a new subfield of planetary science: necroplanetology.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Post-main-sequence debris from rotation-induced YORP break-up of small bodies II: multiple fissions, internal strengths and binary production

Over one quarter of white dwarfs contain observable metallic debris from the breakup of exo-asteroids. Understanding the physical and orbital history of this debris would enable us to self-consistently link planetary system formation and fate. One major debris reservoir is generated by YORP-induced rotational fission during the giant branch phases of stellar evolution, where the stellar luminosity can exceed the Sun's by four orders of magnitude. Here, we determine the efficacy of the giant branch YORP effect for asteroids with nonzero internal strength, and model post-fission evolution by imposing simple analytic fragmentation prescriptions. We find that even the highest realistic internal strengths cannot prevent the widespread fragmentation of asteroids and the production of a debris field over 100 au in size. We compute the number of successive fission events as they occur in progressively smaller time intervals as the star ascends the giant branches, providing a way to generate size distributions of asteroid fragments. The results are highly insensitive to progenitor stellar mass. We also conclude that the ease with which giant branch YORP breakup can generate binary asteroid subsystems is strongly dependent on internal strength. Formed binary subsystems in turn could be short-lived due to the resulting luminosity-enhanced BYORP effect.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Short-term stability of particles in the WD J0914+1914 white dwarf planetary system

Nearly all known white dwarf planetary systems contain detectable rocky debris in the stellar photosphere. A glaring exception is the young and still evolving white dwarf WD J0914+1914, which instead harbours a giant planet and a disc of pure gas. The stability boundaries of this disc and the future prospects for this white dwarf to be polluted with rocks depend upon the mass and orbit of the planet, which are only weakly constrained. Here we combine an ensemble of plausible planet orbits and masses to determine where observers should currently expect to find the outer boundary of the gas disc. We do so by performing a sweep of the entire plausible phase space with short-term numerical integrations. We also demonstrate that particle-star collisional trajectories, which would lead to the (unseen) signature of rocky metal pollution, occupy only a small fraction of the phase space, mostly limited to particle eccentricities above 0.75. Our analysis reveals that a highly inflated planet on a near-circular orbit is the type of planet which is most consistent with the current observations.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

TESS Asteroseismic Analysis of the Known Exoplanet Host Star HD 222076

The Transiting Exoplanet Survey Satellite (TESS) is an all-sky survey mission aiming to search for exoplanets that transit bright stars. The high-quality photometric data of TESS are excellent for the asteroseismic study of solar-like stars. In this work, we present an asteroseismic analysis of the red-giant star HD~222076 hosting a long-period (2.4 yr) giant planet discovered through radial velocities. Solar-like oscillations of HD~222076 are detected around $203 \, μ$Hz by TESS for the first time. Asteroseismic modeling, using global asteroseismic parameters as input, yields a determination of the stellar mass ($M_\star = 1.12 \pm 0.12\, M_\odot$), radius ($R_\star = 4.34 \pm 0.21\,R_\odot$), and age ($7.4 \pm 2.7\,$Gyr), with precisions greatly improved from previous studies. The period spacing of the dipolar mixed modes extracted from the observed power spectrum reveals that the star is on the red-giant branch burning hydrogen in a shell surrounding the core. We find that the planet will not escape the tidal pull of the star and be engulfed into it within about $800\,$Myr, before the tip of the red-giant branch is reached.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

The dynamical history of the evaporating or disrupted ice giant planet around white dwarf WD J0914+1914

Robust evidence of an ice giant planet shedding its atmosphere around the white dwarf WD J0914+1914 represents a milestone in exoplanetary science, allowing us to finally supplement our knowledge of white dwarf metal pollution, debris discs and minor planets with the presence of a major planet. Here, we discuss the possible dynamical origins of this planet, WD J0914+1914 b. The very young cooling age of the host white dwarf (13 Myr) combined with the currently estimated planet-star separation of about 0.07 au imposes particularly intriguing and restrictive coupled constraints on its current orbit and its tidal dissipation characteristics. The planet must have been scattered from a distance of at least a few au to its current location, requiring the current or former presence of at least one more major planet in the system in the absence of a hidden binary companion. We show that WD J0914+1914 b could not have subsequently shrunk its orbit through chaotic f-mode tidal excitation (characteristic of such highly eccentric orbits) unless the planet was or is highly inflated and possibly had partially thermally self-disrupted from mode-based energy release. We also demonstrate that if the planet is currently assumed to reside on a near-circular orbit at 0.07 au, then non-chaotic equilibrium tides impose unrealistic values for the planet's tidal quality factor. We conclude that WD~J0914+1914~b either (i) actually resides interior to 0.07 au, (ii) resembles a disrupted "Super-Puff" whose remains reside on a circular orbit, or (iii) resembles a larger or denser ice giant on a currently eccentric orbit. Distinguishing these three possibilities strongly motivates follow-up observations.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

The grain size survival threshold in one-planet post-main-sequence exoplanetary systems

The size distribution and orbital architecture of dust, grains, boulders, asteroids, and major planets during the giant branch phases of evolution dictate the preponderance and observability of the eventual debris, which have been found to surround white dwarfs and pollute their atmospheres with metals. Here, we utilize the photogravitational planar restricted three-body problem in one-planet giant branch systems in order to characterize the orbits of grains as the parent star luminosity and mass undergo drastic changes. We perform a detailed dynamical analysis of the character of grain orbits (collisional, escape, or bounded) as a function of location and energy throughout giant branch evolution. We find that for stars with main-sequence masses of $2.0M_{\odot}$, giant branch evolution, combined with the presence of a planet, ubiquitously triggers escape in grains smaller than about 1 mm, while leaving grains larger than about 5 cm bound to the star. This result is applicable for systems with either a terrestrial or giant planet, is largely independent of the location of the planet, and helps establish a radiative size threshold for escape of small particles in giant branch planetary systems.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

The lifetimes of planetary debris discs around white dwarfs

The lifetime of a planetary disc which orbits a white dwarf represents a crucial input parameter into evolutionary models of that system. Here we apply a purely analytical formalism to estimate lifetimes of the debris phase of these discs, before they are ground down into dust or are subject to sublimation from the white dwarf. We compute maximum lifetimes for three different types of white dwarf discs, formed from (i) radiative YORP breakup of exo-asteroids along the giant branch phases at 2-100 au, (ii) radiation-less spin-up disruption of these minor planets at $\sim 1.5-4.5R_{\odot}$, and (iii) tidal disruption of minor or major planets within about $1.3R_{\odot}$. We display these maximum lifetimes as a function of disc mass and extent, constituent planetesimal properties, and representative orbital excitations of eccentricity and inclination. We find that YORP discs with masses up to $10^{24}$ kg live long enough to provide a reservoir of surviving cm-sized pebbles and m- to km-sized boulders that can be perturbed intact to white dwarfs with cooling ages of up to 10 Gyr. Debris discs formed from the spin or tidal disruption of these minor planets or major planets can survive in a steady state for up to respectively 1 Myr or 0.01 Myr, although most tidal discs would leave a steady state within about 1 yr. Our results illustrate that dust-less planetesimal transit detections are plausible, and would provide particularly robust evolutionary constraints. Our formalism can easily be adapted to individual systems and future discoveries.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

The white dwarf planet WD J0914+1914 b: barricading potential rocky pollutants?

An ice giant planet was recently reported orbiting white dwarf WD J0914+1914 at an approximate distance of 0.07 au. The striking non-detection of rocky pollutants in this white dwarf's photosphere contrasts with the observations of nearly every other known white dwarf planetary system. Here, I analyze the prospects for exterior extant rocky asteroids, boulders, cobbles and pebbles to radiatively drift inward past the planet due to the relatively high luminosity (0.1 L_Sun) of this particularly young (13 Myr) white dwarf. Pebbles and cobbles drift too slowly from Poynting-Robertson drag to bypass the planet, but boulders and asteroids are subject to the much stronger Yarkovsky effect. In this paper, I (i) place lower limits on the timescales for these objects to reach the planet's orbit, (ii) establish 3 m as the approximate limiting radius above which a boulder drifts too slowly to avoid colliding with the planet, and (iii) compute bounds on the fraction of boulders which succeed in traversing mean motion resonances and the planet's Hill sphere to eventually pollute the star. Overall, I find that the planet acts as a barrier against rather than a facilitator for radiatively-driven rocky pollution, suggesting that future rocky pollutants would most likely originate from distant scattering events.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Chaotic rotation and evolution of asteroids and small planets in high-eccentricity orbits around white dwarfs

Observed planetary debris in white dwarf atmospheres predominately originate from the destruction of small bodies on highly eccentric ($>0.99$) orbits. Despite their importance, these minor planets have coupled physical and orbital evolution which has remained largely unexplored. Here, we present a novel approach for estimating the influence of fast chaotic rotation on the orbital evolution of high-eccentricity triaxial asteroids, and formally characterize the propagation of their angular rotation velocities and orbital elements as random time processes. By employing the impulse approximation, we demonstrate that the violent gravitational interactions during periastron passages transfer energy between the orbit and asteroid's rotation. If the distribution of spin impulses were symmetric around zero, then the net result would be a secular decrease of the semimajor axis and a further increase of the eccentricity. We find evidence, however, that the chaotic rotation may be self-regulated in such a manner that these effects are reduced or nullified. We discover that asteroids on highly eccentric orbits can break themselves apart --- in a type of YORP-less rotational fission --- without actually entering the Roche radius, with potentially significant consequences for the distribution of debris and energy requirements for gravitational scattering in metal-polluted white dwarf planetary systems. This mechanism provides a steady stream of material impacting a white dwarf without rapidly depleting the number of small bodies in the stellar system.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

A simple model to describe intrinsic stellar noise for exoplanet detection around red giants

In spite of the huge advances in exoplanet research provided by the NASA Kepler Mission, there remain only a small number of transit detections around evolved stars. Here we present a reformulation of the noise properties of red-giant stars, where the intrinsic stellar granulation, and the stellar oscillations described by asteroseismology play a key role. The new noise model is a significant improvement on the current Kepler results for evolved stars. Our noise model may be used to help understand planet detection thresholds for the ongoing K2 and upcoming TESS missions, and serve as a predictor of stellar noise for these missions. As an application of our noise model, we explore the minimum detectable planet radii for red giant stars, and find that Neptune sized planets should be detectable around low luminosity red giant branch stars.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Binary star influence on post-main-sequence multi-planet stability

Nearly every star known to host planets will become a white dwarf, and nearly 100 planet-hosts are now known to be accompanied by binary stellar companions. Here, we determine how a binary companion triggers instability in otherwise unconditionally stable single-star two-planet systems during the giant branch and white dwarf phases of the planet host. We perform about 700 full-lifetime (14 Gyr) simulations with A0 and F0 primary stars and secondary K2 companions, and identify the critical binary distance within which instability is triggered at any point during stellar evolution. We estimate this distance to be about seven times the outer planet separation, for circular binaries. Our results help characterize the fates of planetary systems, and in particular which ones might yield architectures that are conducive to generating observable heavy metal pollution in white dwarf atmospheres.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Dynamical mass and multiplicity constraints on co-orbital bodies around stars

Objects transiting near or within the disruption radius of both main sequence (e.g. KOI 1843) and white dwarf (WD 1145+017) stars are now known. Upon fragmentation or disintegration, these planets or asteroids may produce co-orbital configurations of nearly equal-mass objects. However, as evidenced by the co-orbital objects detected by transit photometry in the WD 1145+017 system, these bodies are largely unconstrained in size, mass, and total number (multiplicity). Motivated by potential future similar discoveries, we perform N-body simulations to demonstrate if and how debris masses and multiplicity may be bounded due to second-to-minute deviations and the resulting accumulated phase shifts in the osculating orbital period amongst multiple co-orbital equal point masses. We establish robust lower and upper mass bounds as a function of orbital period deviation, but find the constraints on multiplicity to be weak. We also quantify the fuzzy instability boundary, and show that mutual collisions occur in less than 5%, 10% and 20% of our simulations for masses of 10^{21}, 10^{22} and 10^{23} kg. Our results may provide useful initial rough constraints on other stellar systems with multiple co-orbital bodies.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Explaining the variability of WD 1145+017 with simulations of asteroid tidal disruption

Post-main-sequence planetary science has been galvanised by the striking variability, depth and shape of the photometric transit curves due to objects orbiting white dwarf WD 1145+017, a star which also hosts a dusty debris disc and circumstellar gas, and displays strong metal atmospheric pollution. However, the physical properties of the likely asteroid which is discharging disintegrating fragments remain largely unconstrained from the observations. This process has not yet been modelled numerically. Here, we use the N-body code PKDGRAV to compute dissipation properties for asteroids of different spins, densities, masses, and eccentricities. We simulate both homogeneous and differentiated asteroids, for up to two years, and find that the disruption timescale is strongly dependent on density and eccentricity, but weakly dependent on mass and spin. We find that primarily rocky differentiated bodies with moderate (~3-4 g/cm^3) bulk densities on near-circular (e <~ 0.1) orbits can remain intact while occasionally shedding mass from their mantles. These results suggest that the asteroid orbiting WD 1145+017 is differentiated, resides just outside of the Roche radius for bulk density but just inside the Roche radius for mantle density, and is more akin physically to an asteroid like Vesta instead of one like Itokawa.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Full-lifetime simulations of multiple unequal-mass planets across all phases of stellar evolution

We know that planetary systems are just as common around white dwarfs as around main sequence stars. However, self-consistently linking a planetary system across these two phases of stellar evolution through the violent giant branch poses computational challenges, and previous studies restricted architectures to equal-mass planets. Here, we remove this constraint and perform over 450 numerical integrations over a Hubble time (14 Gyr) of packed planetary systems with unequal-mass planets. We characterize the resulting trends as a function of planet order and mass. We find that intrusive radial incursions in the vicinity of the white dwarf become less likely as the dispersion amongst planet masses increases. The orbital meandering which may sustain a sufficiently dynamic environment around a white dwarf to explain observations is more dependent on the presence of terrestrial-mass planets than any variation in planetary mass. Triggering unpacking or instability during the white dwarf phase is comparably easy for systems of unequal-mass planets and systems of equal-mass planets; instabilities during the giant branch phase remain rare and require fine-tuning of initial conditions. We list the key dynamical features of each simulation individually as a potential guide for upcoming discoveries.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

High-speed photometry of the disintegrating planetesimals at WD1145+017: evidence for rapid dynamical evolution

We obtained high-speed photometry of the disintegrating planetesimals orbiting the white dwarf WD1145+017, spanning a period of four weeks. The light curves show a dramatic evolution of the system since the first observations obtained about seven months ago. Multiple transit events are detected in every light curve, which have varying durations(~3-12min) and depths (~10-60%). The time-averaged extinction is ~11%, much higher than at the time of the Kepler observations. The shortest-duration transits require that the occulting cloud of debris has a few times the size of the white dwarf, longer events are often resolved into the superposition of several individual transits. The transits evolve on time scales of days, both in shape and in depth, with most of them gradually appearing and disappearing over the course of the observing campaign. Several transits can be tracked across multiple nights, all of them recur on periods of ~4.49h, indicating multiple planetary debris fragments on nearly identical orbits. Identifying the specific origin of these bodies within this planetary system, and the evolution leading to their current orbits remains a challenging problem.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Liberating exomoons in white dwarf planetary systems

Previous studies indicate that more than a quarter of all white dwarf (WD) atmospheres are polluted by remnant planetary material, with some WDs being observed to accrete the mass of Pluto in 10^6 years. The short sinking timescale for the pollutants indicate that the material must be frequently replenished. Moons may contribute decisively to this pollution process if they are liberated from their parent planets during the post-main-sequence evolution of the planetary systems. Here, we demonstrate that gravitational scattering events among planets in WD systems easily triggers moon ejection. Repeated close encounters within tenths of a planetary Hill radii are highly destructive to even the most massive, close-in moons. Consequently, scattering increases both the frequency of perturbing agents in WD systems, as well as the available mass of polluting material in those systems, thereby enhancing opportunities for collision and fragmentation and providing more dynamical pathways for smaller bodies to reach the WD. Moreover, during intense scattering, planets themselves have pericenters with respect to the WD of only a fraction of an AU, causing extreme Hill-sphere contraction, and the liberation of moons into WD-grazing orbits. Many of our results are directly applicable to exomoons orbiting planets around main sequence stars.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Mass and eccentricity constraints on the planetary debris orbiting the white dwarf WD 1145+017

Being the first of its kind, the white dwarf WD 1145+017 exhibits a complex system of disintegrating debris which offers a unique opportunity to study its disruption process in real time. Even with plenty of transit observations there are no clear constraints on the masses or eccentricities of such debris. Using $N$-body simulations we show that masses greater than approximately $10^{20}$ kg (a tenth of the mass of Ceres) or orbits that are not nearly circular ($\mathrm{eccentricity}>10^{-3}$) dramatically increase the chances of the system becoming unstable within two years, which would contrast with the observational data over this timespan. We also provide a direct comparison between transit phase shifts detected in the observations and by our numerical simulations.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Post-main-sequence planetary system evolution

The fates of planetary systems provide unassailable insights into their formation and represent rich cross-disciplinary dynamical laboratories. Mounting observations of post-main-sequence planetary systems necessitate a complementary level of theoretical scrutiny. Here, I review the diverse dynamical processes which affect planets, asteroids, comets and pebbles as their parent stars evolve into giant branch, white dwarf and neutron stars. This reference provides a foundation for the interpretation and modelling of currently known systems and upcoming discoveries.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Relating binary-star planetary systems to central configurations

Binary-star exoplanetary systems are now known to be common, for both wide and close binaries. However, their orbital evolution is generally unsolvable. Special cases of the N-body problem which are in fact completely solvable include dynamical architectures known as central configurations. Here, I utilize recent advances in our knowledge of central configurations to assess the plausibility of linking them to coplanar exoplanetary binary systems. By simply restricting constituent masses to be within stellar or substellar ranges characteristic of planetary systems, I find that (i) this constraint reduces by over 90 per cent the phase space in which central configurations may occur, (ii) both equal-mass and unequal-mass binary stars admit central configurations, (iii) these configurations effectively represent different geometrical extensions of the Sun-Jupiter-Trojan-like architecture, (iv) deviations from these geometries are no greater than ten degrees, and (v) the deviation increases as the substellar masses increase. This study may help restrict future stability analyses to architectures which resemble exoplanetary systems, and might hint at where observers may discover dust, asteroids and/or planets in binary star systems.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

The Detection of Dust around NN Ser

Eclipse timing variations observed from the post common-envelope binary (PCEB) NN Ser offer strong evidence in favour of circumbinary planets existing around PCEBs. If real, these planets may be accompanied by a disc of dust. We here present the ALMA detection of flux at 1.3 mm from NN Ser, which is likely due to thermal emission from a dust disc of mass $\sim0.8\pm0.2\,M_{\oplus}$. We performed simulations of the history of NN Ser to determine possible origins of this dust, and conclude that the most likely origin is, in fact, common-envelope material which was not expelled from the system and instead formed a circumbinary disc. These discs have been predicted by theory but previously remained undetected. While the presence of this dust does not prove the existence of planets around NN Ser, it adds credibility to the possibility of planets forming from common-envelope material in a 'second-generation' scenario.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

The Fate of Exomoons in White Dwarf Planetary Systems

Roughly 1000 white dwarfs are known to be polluted with planetary material, and the progenitors of this material are typically assumed to be asteroids. The dynamical architectures which perturb asteroids into white dwarfs are still unknown, but may be crucially dependent on moons liberated from parent planets during post-main-sequence gravitational scattering. Here, we trace the fate of these exomoons, and show that they more easily achieve deep radial incursions towards the white dwarf than do scattered planets. Consequently, moons are likely to play a significant role in white dwarf pollution, and in some cases may be the progenitors of the pollution itself.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

The fates of Solar system analogues with one additional distant planet

The potential existence of a distant planet ("Planet Nine") in the Solar system has prompted a re-think about the evolution of planetary systems. As the Sun transitions from a main sequence star into a white dwarf, Jupiter, Saturn, Uranus and Neptune are currently assumed to survive in expanded but otherwise unchanged orbits. However, a sufficiently-distant and sufficiently-massive extra planet would alter this quiescent end scenario through the combined effects of Solar giant branch mass loss and Galactic tides. Here, I estimate bounds for the mass and orbit of a distant extra planet that would incite future instability in systems with a Sun-like star and giant planets with masses and orbits equivalent to those of Jupiter, Saturn, Uranus and Neptune. I find that this boundary is diffuse and strongly dependent on each of the distant planet's orbital parameters. Nevertheless, I claim that instability occurs more often than not when the planet is as massive as Jupiter and harbours a semimajor axis exceeding about 300 au, or has a mass of a super-Earth and a semimajor axis exceeding about 3000 au. These results hold for orbital pericentres ranging from 100 to at least 400 au. This instability scenario might represent a common occurrence, as potentially evidenced by the ubiquity of metal pollution in white dwarf atmospheres throughout the Galaxy.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

The unstable fate of the planet orbiting the A-star in the HD 131399 triple stellar system

Validated planet candidates need not lie on long-term stable orbits, and instability triggered by post-main-sequence stellar evolution can generate architectures which transport rocky material to white dwarfs, polluting them. The giant planet HD 131399Ab orbits its parent A star at a projected separation of about 50-100 au. The host star, HD 131399A, is part of a hierarchical triple with HD 131399BC being a close binary separated by a few hundred au from the A star. Here, we determine the fate of this system, and find that (i) stability along the main sequence is achieved only for a favourable choice of parameters within the errors, and (ii) even for this choice, in almost every instance the planet is ejected during the transition between the giant branch and white dwarf phases of HD 131399A. This result provides an example of both how the free-floating planet population may be enhanced by similar systems, and how instability can manifest in the polluted white dwarf progenitor population.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

A wide binary trigger for white dwarf pollution

Metal pollution in white dwarf atmospheres is likely to be a signature of remnant planetary systems. Most explanations for this pollution predict a sharp decrease in the number of polluted systems with white dwarf cooling age. Observations do not confirm this trend, and metal pollution in old (1-5 Gyr) white dwarfs is difficult to explain. We propose an alternative, time-independent mechanism to produce the white dwarf pollution. The orbit of a wide binary companion can be perturbed by Galactic tides, approaching close to the primary star for the first time after billions of years of evolution on the white dwarf branch. We show that such a close approach perturbs a planetary system orbiting the white dwarf, scattering planetesimals onto star-grazing orbits, in a manner that could pollute the white dwarf's atmosphere. Our estimates find that this mechanism is likely to contribute to metal pollution, alongside other mechanisms, in up to a few percent of an observed sample of white dwarfs with wide binary companions, independent of white dwarf age. This age independence is the key difference between this wide binary mechanism and others mechanisms suggested in the literature to explain white dwarf pollution. Current observational samples are not large enough to assess whether this mechanism makes a significant contribution to the population of polluted white dwarfs, for which better constraints on the wide binary population are required, such as those that will be obtained in the near future with Gaia.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Doppler-imaging of the planetary debris disc at the white dwarf SDSS J122859.93+104032.9

Debris discs which orbit white dwarfs are signatures of remnant planetary systems. We present twelve years of optical spectroscopy of the metal-polluted white dwarf SDSS J1228+1040, which shows a steady variation in the morphology of the 8600 Å Ca II triplet line profiles from the gaseous component of its debris disc. We identify additional emission lines of O I, Mg I, Mg II, Fe II and Ca II in the deep co-added spectra. These emission features (including Ca H & K) exhibit a wide range in strength and morphology with respect to each other and to the Ca II triplet, indicating different intensity distributions of these ionic species within the disc. Using Doppler tomography we show that the evolution of the Ca II triplet profile can be interpreted as the precession of a fixed emission pattern with a period in the range 24-30 years. The Ca II line profiles vary on time-scales that are broadly consistent with general relativistic precession of the debris disc.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Formation of planetary debris discs around white dwarfs II: Shrinking extremely eccentric collisionless rings

The formation channel of the tens of compact debris discs which orbit white dwarfs (WDs) at a distance of one Solar radius remains unknown. Asteroids that survive the giant branch stellar phases beyond a few au are assumed to be dynamically thrust towards the WD and tidally disrupted within its Roche radius, generating extremely eccentric (e>0.98) rings. Here, we establish that WD radiation compresses and circularizes the orbits of super-micron to cm-sized ring constituents to entirely within the WD's Roche radius. We derive a closed algebraic formula which well-approximates the shrinking time as a function of WD cooling age, the physical properties of the star and the physical and orbital properties of the ring particles. The shrinking timescale increases with both particle size and cooling age, yielding age-dependent WD debris disc size distributions.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

One of the closest exoplanet pairs to the 3:2 Mean Motion Resonance: K2-19b \& c

The K2 mission has recently begun to discover new and diverse planetary systems. In December 2014 Campaign 1 data from the mission was released, providing high-precision photometry for ~22000 objects over an 80 day timespan. We searched these data with the aim of detecting further important new objects. Our search through two separate pipelines led to the independent discovery of K2-19b \& c, a two-planet system of Neptune sized objects (4.2 and 7.2 $R_\oplus$), orbiting a K dwarf extremely close to the 3:2 mean motion resonance. The two planets each show transits, sometimes simultaneously due to their proximity to resonance and alignment of conjunctions. We obtain further ground based photometry of the larger planet with the NITES telescope, demonstrating the presence of large transit timing variations (TTVs), and use the observed TTVs to place mass constraints on the transiting objects under the hypothesis that the objects are near but not in resonance. We then statistically validate the planets through the \texttt{PASTIS} tool, independently of the TTV analysis.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Prospects for detecting decreasing exoplanet frequency with main sequence age using PLATO

The space mission PLATO will usher in a new era of exoplanetary science by expanding our current inventory of transiting systems and constraining host star ages, which are currently highly uncertain. This capability might allow PLATO to detect changes in planetary system architecture with time, particularly because planetary scattering due to Lagrange instability may be triggered long after the system was formed. Here, we utilize previously published instability timescale prescriptions to determine PLATO's capability to detect a trend of decreasing planet frequency with age for systems with equal-mass planets. For two-planet systems, our results demonstrate that PLATO may detect a trend for planet masses which are at least as massive as super-Earths. For systems with three or more planets, we link their initial compactness to potentially detectable frequency trends in order to aid future investigations when these populations will be better characterized.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Stellar archaeology with Gaia: the Galactic white dwarf population

Gaia will identify several 1e5 white dwarfs, most of which will be in the solar neighborhood at distances of a few hundred parsecs. Ground-based optical follow-up spectroscopy of this sample of stellar remnants is essential to unlock the enormous scientific potential it holds for our understanding of stellar evolution, and the Galactic formation history of both stars and planets.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Sublimation-induced orbital perturbations of extrasolar active asteroids and comets: application to white dwarf systems

The metal budgets in some white dwarf (WD) atmospheres reveal that volatile-rich circumstellar bodies must both exist in extrasolar systems and survive the giant branch phases of stellar evolution. The resulting behaviour of these active asteroids or comets which orbit WDs is not well-understood, but may be be strongly influenced by sublimation due to stellar radiation. Here we develop a model, generally applicable to any extrasolar system with a main sequence or WD star, that traces sublimation-induced orbital element changes in approximately km-sized extrasolar minor planets and comets traveling within hundreds of au. We derive evolution equations on orbital timescales and for arbitrarily steep power-law sublimation dependencies on distance, and place our model in a Solar system context. We also demonstrate the importance of coupling sublimation and general relativity, and the orbital consequences of outgassing in arbitrary directions. We prove that nongravitational accelerations alone cannot result in orbit crossing with the WD disruption radius, but may shrink or expand the orbit by up to several au after a single pericentre passage, potentially affecting subsequent interactions with remnant debris and planets. Our analysis suggests that extant planets must exist in polluted WD systems.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

The bulk composition of exo-planets

Priorities in exo-planet research are rapidly moving from finding planets to characterizing their physical properties. Of key importance is their chemical composition, which feeds back into our understanding of planet formation. For the foreseeable future, far-ultraviolet spectroscopy of white dwarfs accreting planetary debris remains the only way to directly and accurately measure the bulk abundances of exo-planetary bodies. The exploitation of this method is limited by the sensitivity of HST, and significant progress will require a large-aperture space telescope with a high-throughput ultraviolet spectrograph.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

The orbital evolution of asteroids, pebbles and planets from giant branch stellar radiation and winds

The discovery of over 50 planets around evolved stars and more than 35 debris discs orbiting white dwarfs highlight the increasing need to understand small body evolution around both early and asymptotic giant branch (GB) stars. Pebbles and asteroids are susceptible to strong accelerations from the intense luminosity and winds of GB stars. Here, we establish equations that can model time-varying GB stellar radiation, wind drag and mass loss. We derive the complete three-dimensional equations of motion in orbital elements due to (1) the Epstein and Stokes regimes of stellar wind drag, (2) Poynting-Robertson drag, and (3) the Yarkovsky drift with seasonal and diurnal components. We prove through averaging that the potential secular eccentricity and inclination excitation due to Yarkovsky drift can exceed that from Poynting-Robertson drag and radiation pressure by at least three orders of magnitude, possibly flinging asteroids which survive YORP spin-up into a widely dispersed cloud around the resulting white dwarf. The GB Yarkovsky effect alone may change an asteroid's orbital eccentricity by ten per cent in just one Myr. Damping perturbations from stellar wind drag can be just as extreme, but are strongly dependent on the highly uncertain local gas density and mean free path length. We conclude that GB radiative and wind effects must be considered when modelling the post-main-sequence evolution of bodies smaller than about 1000 km.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

A simple bound for the variation at closest approach of a small body and star due to general relativity

As a comet, asteroid or planet approaches its parent star, the orbit changes shape due to the curvature of spacetime. For comets in particular, the deviation at the pericentre may noticeably change their ephemerides and affect the dynamics of outgassing, tidal disruption or other processes which act on orbital timescales and are assumed to follow Newtonian gravity. By obtaining and analysing the unaveraged equations of motion in orbital elements due to the dominant post-Newtonian contribution (1PN), I derive a simple analytic expression for the maximum deviation in terms of only the stellar mass and eccentricity of the orbit. This relation can be used to assess the potential importance of including short-period relativistic terms in models containing comets, asteroids or planets, and help determine the level of precision needed in numerical integrations. The magnitude of the deviation in systems with Solar-like stars is typically comparable to the size of comet nuclei, and the direction of the deviation is determined by the eccentricity. I show that for eccentricities above a critical value of approximately 0.359, the direction is away from the star.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Debris from giant impacts between planetary embryos at large orbital radii

We consider the observational signatures of giant impacts between planetary embryos. While the debris released in the impact remains in a clump for only a single orbit, there is a much longer lasting asymmetry caused by the fact that all debris must pass through the collision-point. The resulting asymmetry is stationary, it does not orbit the star. The debris is concentrated in a clump at the collision-point, with a more diffuse structure on the opposite side. The asymmetry lasts for typically around 1000 orbital periods of the progenitor, which can be several Myr at distances of ~50 AU. We describe how the appearance of the asymmetric disc depends on the mass and eccentricity of the progenitor, as well as viewing orientation. The wavelength of observation, which determines the grain sizes probed, is also important. Notably, the increased collision rate of the debris at the collision-point makes this the dominant production site for any secondary dust and gas created. For dust small enough to be removed by radiation pressure, and gas with a short lifetime, this causes their distribution to resemble a jet emanating from the (stationary) collision-point. We suggest that the asymmetries seen at large separations in some debris discs, like Beta Pictoris, could be the result of giant impacts. If so this would indicate that planetary embryos are present and continuing to grow at several tens of AU at ages of up to tens of Myr.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Detectable close-in planets around white dwarfs through late unpacking

Although 25%-50% of white dwarfs (WDs) display evidence for remnant planetary systems, their orbital architectures and overall sizes remain unknown. Vibrant close-in (~1 Solar radius) circumstellar activity is detected at WDs spanning many Gyrs in age, suggestive of planets further away. Here we demonstrate how systems with 4 and 10 closely-packed planets that remain stable and ordered on the main sequence can become unpacked when the star evolves into a WD and experience pervasive inward planetary incursions throughout WD cooling. Our full-lifetime simulations run for the age of the Universe and adopt main sequence stellar masses of 1.5, 2.0 and 2.5 Solar masses, which correspond to the mass range occupied by the progenitors of typical present-day WDs. These results provide (i) a natural way to generate an ever-changing dynamical architecture in post-main-sequence planetary systems, (ii) an avenue for planets to achieve temporary close-in orbits that are potentially detectable by transit photometry, and (iii) a dynamical explanation for how residual asteroids might pollute particularly old WDs.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Eight billion asteroids in the Oort cloud

The Oort cloud is usually thought of as a collection of icy comets inhabiting the outer reaches of the Solar system, but this picture is incomplete. We use simulations of the formation of the Oort cloud to show that ~4% of the small bodies in the Oort cloud should have formed within 2.5 au of the Sun, and hence be ice-free rock-iron bodies. If we assume these Oort cloud asteroids have the same size distribution as their cometary counterparts, the Large Synoptic Survey Telescope should find roughly a dozen Oort cloud asteroids during ten years of operations. Measurement of the asteroid fraction within the Oort cloud can serve as an excellent test of the Solar system's formation and dynamical history. Oort cloud asteroids could be of particular concern as impact hazards as their high mass density, high impact velocity, and low visibility make them both hard to detect and hard to divert or destroy. However, they should be a rare class of object, and we estimate globally catastrophic collisions should only occur about once per billion years.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Explicit evolution relations with orbital elements for eccentric, inclined, elliptic and hyperbolic restricted few-body problems

Planetary, stellar and galactic physics often rely on the general restricted gravitational N-body problem to model the motion of a small-mass object under the influence of much more massive objects. Here, I formulate the general restricted problem entirely and specifically in terms of the commonly-used orbital elements of semimajor axis, eccentricity, inclination, longitude of ascending node, argument of pericentre, and true anomaly, without any assumptions about their magnitudes. I derive the equations of motion in the general, unaveraged case, as well as specific cases, with respect to both a bodycentric and barycentric origin. I then reduce the equations to three-body systems, and present compact singly- and doubly-averaged expressions which can be readily applied to systems of interest. This method recovers classic Lidov-Kozai and Laplace-Lagrange theory in the test particle limit to any order, but with fewer assumptions, and reveals a complete analytic solution for the averaged planetary pericentre precession in coplanar circular circumbinary systems to at least the first three nonzero orders in semimajor axis ratio. Finally, I show how the unaveraged equations may be used to express resonant angle evolution in an explicit manner that is not subject to expansions of eccentricity and inclination about small nor any other values.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Formation of planetary debris discs around white dwarfs I: Tidal disruption of an extremely eccentric asteroid

25%-50% of all white dwarfs (WDs) host observable and dynamically active remnant planetary systems based on the presence of close-in circumstellar dust and gas and photospheric metal pollution. Currently-accepted theoretical explanations for the origin of this matter include asteroids that survive the star's giant branch evolution at au-scale distances and are subsequently perturbed onto WD-grazing orbits following stellar mass loss. In this work we investigate the tidal disruption of these highly-eccentric (e > 0.98) asteroids as they approach and tidally disrupt around the WD. We analytically compute the disruption timescale and compare the result with fully self-consistent numerical simulations of rubble piles by using the N-body code PKDGRAV. We find that this timescale is highly dependent on the orbit's pericentre and largely independent of its semimajor axis. We establish that spherical asteroids readily break up and form highly eccentric collisionless rings, which do not accrete onto the WD without additional forces such as radiation or sublimation. This finding highlights the critical importance of such forces in the physics of WD planetary systems.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Hydrogen delivery onto white dwarfs from remnant exo-Oort cloud comets

The origin of trace hydrogen in white dwarfs (WDs) with He-dominated atmospheres is a long-standing problem, one that cannot satisfactorily be explained by the historically-favoured hypothesis of accretion from the interstellar medium. Here we explore the possibility that the gradual accretion of exo-Oort cloud comets, which are a rich source of H, contributes to the apparent increase of trace H with WD cooling age. We determine how often remnant exo-Oort clouds, freshly excited from post-main-sequence stellar mass loss, dynamically inject comets inside the WD's Roche radius. We improve upon previous studies by considering a representative range of single WD masses (0.52-1.00 Solar masses) and incorporating different cloud architectures, giant branch stellar mass loss, stellar flybys, Galactic tides and a realistic escape ellipsoid in self-consistent numerical simulations that integrate beyond 8 Gyr ages of WD cooling. We find that about 10^{-5} of the material in an exo-Oort cloud is typically amassed onto the WD, and that the H deposits accumulate even as the cloud dissipates. This accumulation may account for the relatively large amount of trace H, 10^{22}-10^{25} g, that is determined frequently among WDs with cooling ages >= 1 Gyr. Our results also reaffirm the notion that exo-Oort cloud comets are not the primary agents of the metal budgets observed in polluted WD atmospheres.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Post-main-sequence debris from rotation-induced YORP break-up of small bodies

Although discs of dust and gas have been observed orbiting white dwarfs, the origin of this circumstellar matter is uncertain. We hypothesize that the in-situ breakup of small bodies such as asteroids spun to fission during the giant branch phases of stellar evolution provides an important contribution to this debris. The YORP effect, which arises from radiation pressure, accelerates the spin rate of asymmetric asteroids, which can eventually shear themselves apart. This pressure is maintained and enhanced around dying stars because the outward push of an asteroid due to stellar mass loss is insignificant compared to the resulting stellar luminosity increase. Consequently, giant star radiation will destroy nearly all bodies with radii in the range 100 m - 10 km that survive their parent star's main sequence lifetime within a distance of about 7 au; smaller bodies are spun apart to their strongest, competent components. This estimate is conservative, and would increase for highly asymmetric shapes or incorporation of the inward drag due to giant star stellar wind. The resulting debris field, which could extend to thousands of au, may be perturbed by remnant planetary systems to reproduce the observed dusty and gaseous discs which accompany polluted white dwarfs.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

A simple scaling for the minimum instability time-scale of two widely spaced planets

Long-term instability in multi-planet exosystems is a crucial consideration when confirming putative candidates, analyzing exoplanet populations, constraining the age of exosystems, and identifying the sources of white dwarf pollution. Two planets which are Hill stable are separated by a wide-enough distance to ensure that they will never collide. However, Hill stable planetary systems may eventually manifest Lagrange instability when the outer planet escapes or the inner planet collides with the star. We show empirically that for two nearly coplanar Hill stable planets with eccentricities less than about 0.3, instability can manifest itself only after a time corresponding to X initial orbits of the inner planet, where log_{10}(X) is of the order of 5.2 mu^{-0.18} and mu is the planet-star mass ratio measured in (Jupiter mass/Solar mass). This relation applies to any type of equal-mass secondaries, and suggests that two low-eccentricity Hill stable terrestrial-mass or smaller-mass planets should be Lagrange stable throughout the main sequence lifetime of any white dwarf progenitor. However, Hill stable giant planets are not guaranteed to be Lagrange stable, particularly within a few tens of percent beyond the critical Hill separation. Our scaling represents a useful "rule of thumb" for planetary population syntheses or individual systems for which performing detailed long-term integrations is unfeasible.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

An Exoplanet's Response to Anisotropic Stellar Mass-Loss During Birth and Death

The birth and death of planets may be affected by mass outflows from their parent stars during the T-Tauri or post-main-sequence phases of stellar evolution. These outflows are often modelled to be isotropic, but this assumption is not realistic for fast rotators, bipolar jets and supernovae. Here we derive the general equations of motion for the time evolution of a single planet, brown dwarf, comet or asteroid perturbed by anisotropic mass loss in terms of a complete set of planetary orbital elements, the ejecta velocity, and the parent star's co-latitude and longitude. We restrict our application of these equations to 1) rapidly rotating giant stars, and 2) arbitrarily-directed jet outflows. We conclude that the isotropic mass-loss assumption can safely be used to model planetary motion during giant branch phases of stellar evolution within distances of hundreds of au. In fact, latitudinal mass loss variations anisotropically affect planetary motion only if the mass loss is asymmetric about the stellar equator. Also, we demonstrate how constant-velocity, asymmetric bipolar outflows in young systems incite orbital inclination changes. Consequently, this phenomenon readily tilts exoplanetary orbits external to a nascent disc on the order of degrees.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Long-term evolution of three-planet systems to the post-Main Sequence and beyond

We study the stability of systems of three giant planets orbiting 3 - 8 M_Sol stars at orbital distances of >10 au as the host star ages through the Main Sequence (MS) and well into the White Dwarf (WD) stage. Systems are stable on the MS if the planets are separated by more than ~9 Hill radii. Most systems surviving the MS will remain stable until the WD phase, although planets scattered onto small pericentres in unstable systems can be swallowed by the expanding stellar envelope when the star ascends the giant branches. Mass loss at the end of the asymptotic giant branch triggers delayed instability in many systems, leading to instabilities typically occurring at WD cooling ages of a few 100 Myr. This instability occurs both in systems that survived the star's previous evolution unscathed, and in systems that previously underwent scattering instabilities. The outcome of such instability around WDs is overwhelmingly the ejection of one of the planets from the system, with several times more ejections occurring during the WD phase than during the MS. Furthermore, few planets are scattered close to the WD, just outside the Roche limit, where they can be tidally circularised. Hence, we predict that planets in WD systems rarely dynamically evolve to become "hot Jupiters". Nor does it appear that the observed frequency of metal pollution in WD atmospheres can be entirely explained by planetesimals being destabilised following instability in systems of multiple giant planets, although further work incorporating low-mass planets and planetesimals is needed.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Main-sequence progenitor configurations of the NN Ser candidate circumbinary planetary system are dynamically unstable

Recent observations of the NN Serpentis post-common envelope binary system have revealed eclipse timing variations that have been attributed to the presence of two Jovian-mass exoplanets. Under the assumption that these planets are real and survived from the binary's Main Sequence state, we reconstruct initial binaries that give rise to the present NN Ser configuration and test the dynamical stability of the original system. Under standard assumptions about binary evolution, we find that survival of the planets through the entire Main Sequence life-time is very unlikely. Hence, we conclude that the planets are not survivors from before the Common Envelope phase, implying that either they formed recently out of material ejected from the primary, or that the observed signals are of non-planetary origin.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Multi-Planet Destabilisation and Escape in Post-Main Sequence Systems

Discoveries of exoplanets orbiting evolved stars motivate critical examinations of the dynamics of $N$-body systems with mass loss. Multi-planet evolved systems are particularly complex because of the mutual interactions between the planets. Here, we study the underlying dynamical mechanisms which can incite planetary escape in two-planet post-main sequence systems. Stellar mass loss alone is unlikely to be rapid and high enough to eject planets at typically-observed separations. However, the combination of mass loss and planet-planet interactions can prompt a shift from stable to chaotic regions of phase space. Consequently, when mass loss ceases, the unstable configuration may cause escape. By assuming a constant stellar mass loss rate, we utilize maps of dynamical stability to illustrate the distribution of regular and chaotic trajectories in phase space. We show that chaos can drive the planets to undergo close encounters, leading to the ejection of one planet. Stellar mass loss can trigger the transition of a planetary system from a stable to chaotic configuration, subsequently causing escape. We find that mass loss non-adiabatically affects planet-planet interaction for the most massive progenitor stars which avoid the supernova stage. For these cases, we present specific examples of planetary escape.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Simulations of two-planet systems through all phases of stellar evolution: implications for the instability boundary and white dwarf pollution

Exoplanets have been observed at many stages of their host star's life, including the main sequence (MS), subgiant and red giant branch stages. Also, polluted white dwarfs (WDs) likely represent dynamically active systems at late times. Here, we perform 3-body simulations which include realistic post-MS stellar mass loss and span the entire lifetime of exosystems with two massive planets, from the endpoint of formation to several Gyr into the WD phase of the host star. We find that both MS and WD systems experience ejections and star-planet collisions (Lagrange instability) even if the planet-planet separation well-exceeds the analytical orbit-crossing (Hill instability) boundary. Consequently, MS-stable planets do not need to be closely-packed to experience instability during the WD phase. This instability may pollute the WD directly through collisions, or, more likely, indirectly through increased scattering of smaller bodies such as asteroids or comets. Our simulations show that this instability occurs predominately between tens of Myr to a few Gyrs of WD cooling.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

The Great Escape III: Placing post-main-sequence evolution of planetary and binary systems in a Galactic context

Our improving understanding of the life cycle of planetary systems prompts investigations of the role of the Galactic environment before, during and after Asymptotic Giant Branch (AGB) stellar evolution. Here, we investigate the interplay between stellar mass loss, Galactic tidal perturbations, and stellar flybys for evolving stars which host one planet, smaller body or stellar binary companion and reside in the Milky Way's bulge or disc. We find that the potential evolutionary pathways from a main sequence (MS) to a white dwarf (WD) planetary system are a strong function of Galactocentric distance only with respect to the prevalence of stellar flybys. Planetary ejection and collision with the parent star should be more common towards the bulge. At a given location anywhere in the Galaxy, if the mass loss is adiabatic, then the secondary is likely to avoid close flybys during AGB evolution, and cannot eventually escape the resulting WD because of Galactic tides alone. Partly because AGB mass loss will shrink a planetary system's Hill ellipsoid axes by about 20 to 40 per cent, Oort clouds orbiting WDs are likely to be more depleted and dynamically excited than on the MS.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

The Long-Term Dynamical Evolution of Planetary Systems

This chapter concerns the long-term dynamical evolution of planetary systems from both theoretical and observational perspectives. We begin by discussing the planet-planet interactions that take place within our own Solar System. We then describe such interactions in more tightly-packed planetary systems. As planet-planet interactions build up, some systems become dynamically unstable, leading to strong encounters and ultimately either ejections or collisions of planets. After discussing the basic physical processes involved, we consider how these interactions apply to extrasolar planetary systems and explore the constraints provided by observed systems. The presence of a residual planetesimal disc can lead to planetary migration and hence cause instabilities induced by resonance crossing; however, such discs can also stabilise planetary systems. The crowded birth environment of a planetary system can have a significant impact: close encounters and binary companions can act to destabilise systems, or sculpt their properties. In the case of binaries, the Kozai mechanism can place planets on extremely eccentric orbits which may later circularise to produce hot Jupiters.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Disrupting Primordial Planet Signatures: The Close Encounter of Two Single-Planet Exosystems in the Galactic Disc

During their main sequence lifetimes, the majority of all Galactic Disc field stars must endure at least one stellar intruder passing within a few hundred AU. Mounting observations of planet-star separations near or beyond this distance suggest that these close encounters may fundamentally shape currently-observed orbital architectures and hence obscure primordial orbital features. We consider the commonly-occurring fast close encounters of two single-planet systems in the Galactic Disc, and investigate the resulting change in the planetary eccentricity and semimajor axis. We derive explicit 4-body analytical limits for these variations and present numerical cross-sections which can be applied to localized regions of the Galaxy. We find that each wide-orbit planet has a few percent chance of escape and an eccentricity that will typically change by at least 0.1 due to these encounters. The orbital properties established at formation of millions of tight-orbit Milky Way exoplanets are likely to be disrupted.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Exoplanets Beyond the Solar Neighbourhood: Galactic Tidal Perturbations

The majority of Milky Way extrasolar planets likely reside within a few kpc of the Galactic centre. The Galactic tidal forces acting on planets scale inversely with radius in the Galaxy and so are much greater in the inner Galaxy than in the Solar neighbourhood. Within a range of 3.5 to 10 kpc, the vertical tide from the Galactic disc is predominant. Interior to 3.5 kpc, the effects of the Galactic bulge cannot be neglected and the in-plane tidal components are as important as the vertical ones. Here, we quantify the orbital changes induced by these tides. We find that the greatest perturbations occur when the planetary orbit is severely misaligned to the parent star's orbit. When both planes are perpendicular, the eccentricity of the planet is driven to unity, although the semimajor axis is secularly unaffected. When both planes are coincident, the effect from Galactic tides is minimized, but remains non-zero. In these cases, we provide estimates for the survival times, as well as the minimum baseline eccentricity variation for all Milky Way exoplanets as a function of Galactic parameters. Inclinations similar to the Solar System's (about 60 degrees) can easily cause eccentric Neptunes (at about 30 AU) around host stars deep within the Galactic bulge (within 50 pc) to experience eccentricity variations of several tenths, and cause the exoplanets with the widest-known separations (at about 1000 AU) to experience similar variations in the Galactic disc. These variations occur on timescales of a few Gyr, a fraction of a typical main sequence lifetime.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Exoplanets Bouncing Between Binary Stars

Exoplanetary systems are found not only among single stars, but also binaries of widely varying parameters. Binaries with separations of 100--1000 au are prevalent in the Solar neighborhood; at these separations planet formation around a binary member may largely proceed as if around a single star. During the early dynamical evolution of a planetary system, planet--planet scattering can eject planets from a star's grasp. In a binary, the motion of a planet ejected from one star has effectively entered a restricted three-body system consisting of itself and the two stars, and the equations of motion of the three body problem will apply as long as the ejected planet remains far from the remaining planets. Depending on its energy, escape from the binary as a whole may be impossible or delayed until the three-body approximation breaks down, and further close interactions with its planetary siblings boost its energy when it passes close to its parent star. Until then this planet may be able to transition from the space around one star to the other, and chaotically `bounce' back and forth. In this paper we directly simulate scattering planetary systems that are around one member of a circular binary, and quantify the frequency of bouncing in scattered planets. We find that a great majority (70 to 85 per cent) of ejected planets will pass at least once through the space of it's host's binary companion, and depending on the binary parameters about 45 to 75 per cent will begin bouncing. The time spent bouncing is roughly log-normally distributed with a peak at about $10^4$ years, with only a small percentage bouncing for more than a Myr. This process may perturb and possibly incite instability among existing planets around the companion star. In rare cases, the presence of multiple planets orbiting both stars may cause post-bouncing capture or planetary swapping.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Planet-Planet Scattering Alone Cannot Explain the Free-Floating Planet Population

Recent gravitational microlensing observations predict a vast population of free-floating giant planets that outnumbers main sequence stars almost twofold. A frequently-invoked mechanism for generating this population is a dynamical instability that incites planet-planet scattering and the ejection of one or more planets in isolated main sequence planetary systems. Here, we demonstrate that this process alone probably cannot represent the sole source of these galactic wanderers. By using straightforward quantitative arguments and N-body simulations, we argue that the observed number of exoplanets exceeds the plausible number of ejected planets per system from scattering. Thus, other potential sources of free-floaters, such as planetary stripping in stellar clusters and post-main-sequence ejection, must be considered.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Planetary Orbital Equations in Externally-Perturbed Systems: Position and Velocity-Dependent Forces

The increasing number and variety of extrasolar planets illustrates the importance of characterizing planetary perturbations. Planetary orbits are typically described by physically intuitive orbital elements. Here, we explicitly express the equations of motion of the unaveraged perturbed two-body problem in terms of planetary orbital elements by using a generalized form of Gauss' equations. We consider a varied set of position and velocity-dependent perturbations, and also derive relevant specific cases of the equations: when they are averaged over fast variables (the "adiabatic" approximation), and in the prograde and retrograde planar cases. In each instance, we delineate the properties of the equations. As brief demonstrations of potential applications, we consider the effect of Galactic tides. We measure the effect on the widest-known exoplanet orbit, Sedna-like objects, and distant scattered disk objects, particularly with regard to where the adiabatic approximation breaks down. The Mathematica code which can help derive the equations is freely available upon request.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

The Great Escape II: Exoplanet Ejection from Dying Multiple Star Systems

Extrasolar planets and belts of debris orbiting post-main-sequence single stars may become unbound as the evolving star loses mass. In multiple star systems, the presence or co-evolution of the additional stars can significantly complicate the prospects for orbital excitation and escape. Here, we investigate the dynamical consequences of multi-phasic, nonlinear mass loss and establish a criterion for a system of any stellar multiplicity to retain a planet whose orbit surrounds all of the parent stars. For single stars which become white dwarfs, this criterion can be combined with the Chandrasekhar Limit to establish the maximum allowable mass loss rate for planet retention. We then apply the criterion to circumbinary planets in evolving binary systems over the entire stellar mass phase space. Through about 10^5 stellar evolutionary track realizations, we characterize planetary ejection prospects as a function of binary separation, stellar mass and metallicity. This investigation reveals that planets residing at just a few tens of AU from a central concentration of stars are susceptible to escape in a wide variety of multiple systems. Further, planets are significantly more susceptible to ejection from multiple star systems than from single star systems for a given system mass. For system masses greater than about 2 Solar masses, multiple star systems represent the greater source of free-floating planets.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

The Solar System's Post-Main Sequence Escape Boundary

The Sun will eventually lose about half of its current mass nonlinearly over several phases of post-main sequence evolution. This mass loss will cause any surviving orbiting body to increase its semimajor axis and perhaps vary its eccentricity. Here, we use a range of Solar models spanning plausible evolutionary sequences and assume isotropic mass loss to assess the possibility of escape from the Solar System. We find that the critical semimajor axis in the Solar System within which an orbiting body is guaranteed to remain bound to the dying Sun due to perturbations from stellar mass loss alone is approximately 1,000 AU - 10,000 AU. The fate of objects near or beyond this critical semimajor axis, such as the Oort Cloud, outer scattered disc and specific bodies such as Sedna, will significantly depend on their locations along their orbits when the Sun turns off of the main sequence. These results are applicable to any exoplanetary system containing a single star with a mass, metallicity and age which are approximately equal to the Sun's, and suggest that few extrasolar Oort Clouds could survive post-main sequence evolution intact.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Traditional formation scenarios fail to explain 4:3 mean motion resonances

At least two multi-planetary systems in a 4:3 mean motion resonance have been found by radial velocity surveys. These planets are gas giants and the systems are only stable when protected by a resonance. Additionally the Kepler mission has detected at least 4 strong candidate planetary systems with a period ratio close to 4:3. This paper investigates traditional dynamical scenarios for the formation of these systems. We systematically study migration scenarios with both N-body and hydro-dynamic simulations. We investigate scenarios involving the in-situ formation of two planets in resonance. We look at the results from finely tuned planet-planet scattering simulations with gas disk damping. Finally, we investigate a formation scenario involving isolation-mass embryos. Although the combined planet-planet scattering and damping scenario seems promising, none of the above scenarios is successful in forming enough systems in 4:3 resonance with planetary masses similar to the observed ones. This is a negative result but it has important implications for planet formation. Previous studies were successful in forming 2:1 and 3:2 resonances. This is generally believed to be evidence of planet migration. We highlight the main differences between those studies and our failure in forming a 4:3 resonance. We also speculate on more exotic and complicated ideas. These results will guide future investigators toward exploring the above scenarios and alternative mechanisms in a more general framework.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

A Bayesian Surrogate Model for Rapid Time Series Analysis and Application to Exoplanet Observations

We present a Bayesian surrogate model for the analysis of periodic or quasi-periodic time series data. We describe a computationally efficient implementation that enables Bayesian model comparison. We apply this model to simulated and real exoplanet observations. We discuss the results and demonstrate some of the challenges for applying our surrogate model to realistic exoplanet data sets. In particular, we find that analyses of real world data should pay careful attention to the effects of uneven spacing of observations and the choice of prior for the "jitter" parameter.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Identifying Non-Resonant Kepler Planetary Systems

The Kepler mission has discovered a plethora of multiple transiting planet candidate exosystems, many of which feature putative pairs of planets near mean motion resonance commensurabilities. Identifying potentially resonant systems could help guide future observations and enhance our understanding of planetary formation scenarios. We develop and apply an algebraic method to determine which Kepler 2-planet systems cannot be in a 1st-4th order resonance, given the current, publicly available data. This method identifies when any potentially resonant angle of a system must circulate. We identify and list 70 near-resonant systems which cannot actually reside in resonance, assuming a widely-used formulation for deriving planetary masses from their observed radii and that these systems do not contain unseen bodies that affect the interactions of the observed planets. This work strengthens the argument that a high fraction of exoplanetary systems may be near resonance but not actually in resonance.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

The California Planet Survey III. A Possible 2:1 Resonance in the Exoplanetary Triple System HD 37124

We present new radial velocities from Keck Observatory and both Newtonian and Keplerian solutions for the triple-planet system orbiting HD 37124. The orbital solution for this system has improved dramatically since the third planet was first reported in Vogt et al. 2005 with an ambiguous orbital period. We have resolved this ambiguity, and the outer two planets have an apparent period commensurability of 2:1. A dynamical analysis finds both resonant and non-resonant configurations consistent with the radial velocity data, and constrains the mutual inclinations of the planets to be less than about 30 degrees. We discuss HD 37124 in the context of the other 19 exoplanetary systems with apparent period commenserabilities, which we summarize in a table. We show that roughly one in three well-characterized multiplanet systems has a apparent low-order period commensuribility, which is more than would naively be expected if the periods of exoplanets in known multiplanet systems were drawn randomly from the observed distribution of planetary orbital periods.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

The Great Escape: How Exoplanets and Smaller Bodies Desert Dying Stars

Mounting discoveries of extrasolar planets orbiting post-main sequence stars motivate studies aimed at understanding the fate of these planets. In the traditional "adiabatic" approximation, a secondary's eccentricity remains constant during stellar mass loss. Here, we remove this approximation, investigate the full two-body point-mass problem with isotropic mass loss, and illustrate the resulting dynamical evolution. The magnitude and duration of a star's mass loss combined with a secondary's initial orbital characteristics might provoke ejection, modest eccentricity pumping, or even circularisation of the orbit. We conclude that Oort clouds and wide-separation planets may be dynamically ejected from 1-7 Solar-mass parent stars during AGB evolution. The vast majority of planetary material which survives a supernova from a 7-20 Solar-mass progenitor will be dynamically ejected from the system, placing limits on the existence of first-generation pulsar planets. Planets around >20 Solar-mass black hole progenitors may easily survive or readily be ejected depending on the core collapse and superwind models applied. Material ejected during stellar evolution might contribute significantly to the free-floating planetary population.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Quantifying the challenges of detecting unseen planetary companions with transit timing variations

Both ground and space-based transit observatories are poised to significantly increase the number of known transiting planets and the number of precisely measured transit times. The variation in a planet's transit times may be used to infer the presence of additional planets. Deducing the masses and orbital parameters of such planets from transit time variations (TTVs) alone is a rich and increasingly relevant dynamical problem. In this work, we evaluate the extent of the degeneracies in this process, systematically explore the dependence of TTV signals on several parameters and provide phase space plots that could aid observers in planning future observations. Our explorations are focused on a likely-to-be prevalent situation: a known transiting short-period Neptune or Jupiter-sized planet and a suspected external low-mass perturber on a nearly-coplanar orbit. Through approximately 10^7 N-body simulations, we demonstrate how TTV signal amplitudes may vary by orders of magnitude due to slight variations in any one orbital parameter (0.001 AU in semimajor axis, 0.005 in eccentricity, or a few degrees in orbital angles), and quantify the number of consecutive transit observations necessary in order to obtain a reasonable opportunity to characterize the unseen planet (approximately greater or equal to 50 observations). Planets in or near period commensurabilities of the form p:q, where p < 21 and q < 4, produce distinct TTV signatures, regardless of whether the planets are actually locked in a mean motion resonance. We distinguish these systems from the secular systems in our explorations. Additionally, we find that computing the autocorrelation function of a TTV signal can provide a useful diagnostic for identifying possible orbits for additional planets and suggest that this method could aid integration of TTV signals in future studies of particular exosystems.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Secular Orbital Dynamics of Hierarchical Two Planet Systems

The discovery of multi-planet extrasolar systems has kindled interest in using their orbital evolution as a probe of planet formation. Accurate descriptions of planetary orbits identify systems which could hide additional planets or be in a special dynamical state, and inform targeted follow-up observations. We combine published radial velocity data with Markov Chain Monte Carlo analyses in order to obtain an ensemble of masses, semimajor axes, eccentricities and orbital angles for each of 5 dynamically active multi-planet systems: HD 11964, HD 38529, HD 108874, HD 168443, and HD 190360. We dynamically evolve these systems using 52,000 long-term N-body integrations that sample the full range of possible line-of-sight and relative inclinations, and we report on the system stability, secular evolution and the extent of the resonant interactions. We find that planetary orbits in hierarchical systems exhibit complex dynamics and can become highly eccentric and maybe significantly inclined. Additionally we incorporate the effects of general relativity in the long-term simulations and demonstrate that can qualitatively affect the dynamics of some systems with high relative inclinations. The simulations quantify the likelihood of different dynamical regimes for each system and highlight the dangers of restricting simulation phase space to a single set of initial conditions or coplanar orbits.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Transit Timing Variations for Inclined and Retrograde Exoplanetary Systems

We perform numerical calculations of the expected transit timing variations (TTVs) induced on a Hot-Jupiter by an Earth-mass perturber. Motivated by the recent discoveries of retrograde transiting planets, we concentrate on an investigation of the effect of varying relative planetary inclinations, up to and including completely retrograde systems. We find that planets in low order (E.g. 2:1) mean-motion resonances (MMRs) retain approximately constant TTV amplitudes for $0<\,^{\circ}i<170\,^{\circ}$, only reducing in amplitude for $i > 170\,^{\circ}$. Systems in higher order MMRs (E.g. 5:1) increase in TTV amplitude as inclinations increase towards $45\,^{\circ}$, becoming approximately constant for $45 < i < 135$, and then declining for $i > 135\,^{\circ}$. Planets away from resonance slowly decrease in TTV amplitude as inclinations increase from 0 to 180, where-as planets adjacent to resonances can exhibit a huge range of variability in TTV amplitude as a function of both eccentricity and inclination. For highly retrograde systems ($135\,^{\circ} < i \leq 180\,^{\circ}$), TTV signals will be undetectable across almost the entirety of parameter space, with the exceptions occurring when the perturber has high eccentricity or is very close to a MMR. This high inclination decrease in TTV amplitude (on and away from resonance) is important for the analysis of the known retrograde and multi-planet transiting systems, as inclination effects need to be considered if TTVs are to be used to exclude the presence of any putative planetary companions: absence of evidence is not evidence of absence.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Formation, Survival, and Detectability of Planets Beyond 100 AU

Direct imaging searches have begun to detect planetary and brown dwarf companions and to place constraints on the presence of giant planets at large separations from their host star. This work helps to motivate such planet searches by predicting a population of young giant planets that could be detectable by direct imaging campaigns. Both the classical core accretion and the gravitational instability model for planet formation are hard-pressed to form long-period planets in situ. Here, we show that dynamical instabilities among planetary systems that originally formed multiple giant planets much closer to the host star could produce a population of giant planets at large (~100 AU - 100000 AU) separations. We estimate the limits within which these planets may survive, quantify the efficiency of gravitational scattering into both stable and unstable wide orbits, and demonstrate that population analyses must take into account the age of the system. We predict that planet scattering creates a population of detectable giant planets on wide orbits that decreases in number on timescales of ~10 Myr. We demonstrate that several members of such populations should be detectable with current technology, quantify the prospects for future instruments, and suggest how they could place interesting constraints on planet formation models.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Planet-planet scattering leads to tightly packed planetary systems

The known extrasolar multiple-planet systems share a surprising dynamical attribute: they cluster just beyond the Hill stability boundary. Here we show that the planet-planet scattering model, which naturally explains the observed exoplanet eccentricity distribution, can reproduce the observed distribution of dynamical configurations. We calculated how each of our scattered systems would appear over an appropriate range of viewing geometries; as Hill stability is weakly dependent on the masses, the mass-inclination degeneracy does not significantly affect our results. We consider a wide range of initial planetary mass distributions and find that some are poor fits to the observed systems. In fact, many of our scattering experiments overproduce systems very close to the stability boundary. The distribution of dynamical configurations of two-planet systems actually may provide better discrimination between scattering models than the distribution of eccentricity. Our results imply that, at least in their inner regions which are weakly affected by gas or planetesimal disks, planetary systems should be "packed", with no large gaps between planets.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

The Formation Mechanism of Gas Giants on Wide Orbits

The recent discoveries of massive planets on ultra-wide orbits of HR 8799 (Marois et al. 2008) and Fomalhaut (Kalas et al. 2008) present a new challenge for planet formation theorists. Our goal is to figure out which of three giant planet formation mechanisms--core accretion (with or without migration), scattering from the inner disk, or gravitational instability--could be responsible for Fomalhaut b, HR 8799 b, c and d, and similar planets discovered in the future. This paper presents the results of numerical experiments comparing the long-period planet formation efficiency of each possible mechanism in model A star, G star and M star disks. First, a simple core accretion simulation shows that planet cores forming beyond 35 AU cannot reach critical mass, even under the most favorable conditions one can construct. Second, a set of N-body simulations demonstrates that planet-planet scattering does not create stable, wide-orbit systems such as HR 8799. Finally, a linear stability analysis verifies previous work showing that global spiral instabilities naturally arise in high-mass disks. We conclude that massive gas giants on stable orbits with semimajor axes greater than 35 AU form by gravitational instability in the disk. We recommend that observers examine the planet detection rate as a function of stellar age, controlling for planet dimming with time. If planet detection rate is found to be independent of stellar age, it would confirm our prediction that gravitational instability is the dominant mode of producing detectable planets on wide orbits. We also predict that the occurrence ratio of long-period to short-period gas giants should be highest for M dwarfs due to the inefficiency of core accretion and the expected small fragment mass in their disks.

0 citations0 reviewsNo code linkedOpen access
preprint2007arXiv

A resonant-term-based model including a nascent disk, precession, and oblateness: application to GJ 876

Investigations of two resonant planets orbiting a star or two resonant satellites orbiting a planet often rely on a few resonant and secular terms in order to obtain a representative quantitative description of the system's dynamical evolution. We present a semianalytic model which traces the orbital evolution of any two resonant bodies in a first- through fourth-order eccentricity or inclination-based resonance dominated by the resonant and secular arguments of the user's choosing. By considering the variation of libration width with different orbital parameters, we identify regions of phase space which give rise to different resonant ''depths,'' and propose methods to model libration profiles. We apply the model to the GJ 876 extrasolar planetary system, quantify the relative importance of the relevant resonant and secular contributions, and thereby assess the goodness of the common approximation of representing the system by just the presumably dominant terms. We highlight the danger in using ''order'' as the metric for accuracy in the orbital solution by revealing the unnatural libration centers produced by the second-order, but not first-order, solution, and by demonstrating that the true orbital solution lies somewhere ''in-between'' the third- and fourth-order solutions. We also present formulas used to incorporate perturbations from central-body oblateness and precession, and a protoplanetary or protosatellite thin disk with gaps, into a resonant system. We quantify these contributions to the GJ 876 system, and thereby highlight the conditions which must exist for multi-planet exosystems to be significantly influenced by such factors. We find that massive enough disks may convert resonant libration into circulation; such disk-induced signatures may provide constraints for future studies of exoplanet systems.

0 citations0 reviewsNo code linkedOpen access
preprint2007arXiv

Extrasolar planetary dynamics with a generalized planar Laplace-Lagrange secular theory

The dynamical evolution of nearly half of the known extrasolar planets in multiple-planet systems may be dominated by secular perturbations. The commonly high eccentricities of the planetary orbits calls into question the utility of the traditional Laplace-Lagrange (LL) secular theory in analyses of the motion. We analytically generalize this theory to fourth-order in the eccentricities, compare the result with the second-order theory and octupole-level theory, and apply these theories to the likely secularly-dominated HD 12661, HD 168443, HD 38529 and Ups And multi-planet systems. The fourth-order scheme yields a multiply-branched criterion for maintaining apsidal libration, and implies that the apsidal rate of a small body is a function of its initial eccentricity, dependencies which are absent from the traditional theory. Numerical results indicate that the primary difference the second and fourth-order theories reveal is an alteration in secular periodicities, and to a smaller extent amplitudes of the planetary eccentricity variation. Comparison with numerical integrations indicates that the improvement afforded by the fourth-order theory over the second-order theory sometimes dwarfs the improvement needed to reproduce the actual dynamical evolution. We conclude that LL secular theory, to any order, generally represents a poor barometer for predicting secular dynamics in extrasolar planetary systems, but does embody a useful tool for extracting an accurate long-term dynamical description of systems with small bodies and/or near-circular orbits.

0 citations0 reviewsNo code linkedOpen access