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

Alexander J. Mustill

Alexander J. Mustill 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.GAphysics.class-ph

Catalog footprint

What is connected

23works
4topics
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

23 published item(s)

preprint2022arXiv

A low-eccentricity migration pathway for a 13-h-period Earth analogue in a four-planet system

It is commonly accepted that exoplanets with orbital periods shorter than 1 day, also known as ultra-short period (USP) planets, formed further out within their natal protoplanetary disk, before migrating to their current-day orbits via dynamical interactions. One of the most accepted theories suggests a violent scenario involving high-eccentricity migration followed by tidal circularization. Here, we present the discovery of a four planet system orbiting the bright (V=10.5) K6 dwarf star TOI-500. The innermost planet is a transiting, Earth-sized USP planet with an orbital period of $\sim$ 13 hours, a mass of 1.42 $\pm$ 0.18 M$_{\oplus}$, a radius of $1.166^{0.061}_{-0.058}$ R$_{\oplus}$, and a mean density of 4.89$^{+1.03}_{-0.88}$ gcm$^{-3}$. Via Doppler spectroscopy, we discovered that the system hosts three outer planets on nearly circular orbits with periods of 6.6, 26.2, and 61.3d and minimum masses of 5.03 $\pm$ 0.41 M$_{\oplus}$, 33.12 $\pm$ 0.88 M$_{\oplus}$ and 15.05$^{+1.12}_{-1.11}$ M$_{\oplus}$, respectively. The presence of both a USP planet and a low-mass object on a 6.6-day orbit indicates that the architecture of this system can be explained via a scenario in which the planets started on low-eccentricity orbits, then moved inwards through a quasi-static secular migration. Our numerical simulations show that this migration channel can bring TOI-500 b to its current location in 2 Gyrs, starting from an initial orbit of 0.02au. TOI-500 is the first four planet system known to host a USP Earth analog whose current architecture can be explained via a non-violent migration scenario.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Hot Jupiters, cold kinematics: High phase space densities of host stars reflect an age bias

Context. The birth environments of planetary systems are thought to influence planet formation and orbital evolution, through external photoevaporation and stellar flybys. Recent work has claimed observational support for this, in the form of a correlation between the properties of planetary systems and the local Galactic phase space density of the host star. In particular, Hot Jupiters are found overwhelmingly around stars in regions of high phase space density, which may reflect a formation environment with high stellar density. Aims. We instead investigate whether the high phase space density may have a galactic kinematic origin: Hot Jupiter hosts may be biased towards being young and therefore kinematically cold, because tidal inspiral leads to the destruction of the planets on Gyr timescales, and the velocity dispersion of stars in the Galaxy increases on similar timescales. Methods. We use 6D positions and kinematics from Gaia for the Hot Jupiter hosts and their neighbours, and construct distributions of the phase space density. We investigate correlations between the stars' local phase space density and peculiar velocity. Results. We find a strong anticorrelation between the phase space density and the host star's peculiar velocity with respect to the Local Standard of Rest. Therefore, most stars in "high-density" regions are kinematically cold, which may be caused by the aforementioned bias towards detecting Hot Jupiters around young stars before the planets' tidal destruction. Conclusions. We do not find evidence in the data for Hot Jupiter hosts preferentially being in phase space overdensities compared to other stars of similar kinematics, nor therefore for their originating in birth environments of high stellar density.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Metal pollution of the solar white dwarf by solar system small bodies

White dwarfs (WDs) often show metal lines in their spectra, indicating accretion of asteroidal material. Our Sun is to become a WD in several Gyr. Here, we examine how the solar WD accretes from the three major small body populations: the main belt asteroids (MBAs), Jovian trojan asteroids (JTAs), and trans-Neptunian objects (TNOs). Owing to the solar mass loss during the giant branch, 40\% of the JTAs are lost but the vast majority of MBAs and TNOs survive. During the WD phase, objects from all three populations are sporadically scattered onto the WD, implying ongoing accretion. For young cooling ages $\lesssim 100$ Myr, accretion of MBAs predominates; our predicted accretion rate $\sim10^6$ g/s falls short of observations by two orders of magnitude. On Gyr timescales, thanks to the consumption of the TNOs that kicks in $\gtrsim 100$ Myr, the rate oscillates around $10^6-10^7$ g/s until several Gyr and drops to $\sim10^5$ g/s at 10 Gyr. Our solar WD accretion rate from 1 Gyr and beyond agrees well with those of the extrasolar WDs. We show that for the solar WD, the accretion source region evolves in an inside-out pattern. Moreover, in a realistic small body population with individual sizes covering a wide range as WD pollutants, the accretion is dictated by the largest objects. As a consequence, the accretion rate is lower by an order of magnitude than that from a population of bodies of a uniform size and the same total mass and shows greater scatter.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Capture of satellites during planetary encounters A case study of the Neptunian moons Triton and Nereid

Single-binary scattering may lead to an exchange where the single object captures a component of the binary, forming a new binary. This has been well studied in encounters between a star--planet pair and a single star. Here we explore the application of the exchange mechanism to a planet--satellite pair and another planet in the gravitational potential of a central star. As a case study, we focus on encounters between a satellite-bearing object and Neptune. We investigate whether Neptune can capture satellites from that object and if the captured satellites have orbits analogous to the Neptunian moons Triton and Nereid. Using $N$-body simulations, we study the capture probability at different encounter distances. Post-capture, we use a simple analytical argument to estimate how the captured orbits evolve under collisional and tidal effects. We find that the average capture probability reaches $\sim$$10\%$ if Neptune penetrates the donor planet's satellite system. Most moons grabbed by Neptune acquire highly eccentric orbits. Post-capture, around half of those captured, especially those on tight orbits, can be circularised, either by tides only or by collisions+tides, turning into Triton-like objects. Captures further out, on the other hand, stay on wide and eccentric orbits like that of Nereid. Both moon types can be captured in the same encounter and they have wide distributions in orbital inclination. Therefore, Triton naturally has a $\sim$50\% chance of being retrograde. A similar process potentially applies to an exoplanetary system, and our model predicts that exomoons can jump from one planet to another during planetary scattering. Specifically, there should be two distinct populations of captured moons: one on close-in circular orbits and the other on far-out eccentric orbits. The two populations may have highly inclined prograde or retrograde orbits.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Effects of capturing a wide-orbit planet on planetary systems: system stability and Habitable Zone bombardment rates

A large fraction of stars are formed in dense clusters. In the cluster, close encounters between stars at distances less than 100 au are common. It has been shown that during close encounters planets can transfer between stars. Such captured planets will be on different orbits compared to planets formed in the system, often on very wide, eccentric and inclined orbits. We examine how these captured planets affect Kuiper-belt like asteroid belts in their new systems, and how this affects habitable planets in the system. We show that these captured planets can destabilize the asteroid belt, and we show that the fraction of the asteroid that make it past the giant planets into the system to impact the habitable planet is independent of the captured planets orbital plane, whereas the fraction of the asteroids that are removed and the rate at which they are removed depend strongly on the captured planets pericentre and inclination. We then examine all possible outcomes of planet capture and find that when a Jupiter-mass planet is captured it will in 40\% of cases destabilize the planets in the system, in 40\% of cases deplete the asteroid belt in a few Myr, i.e. not posing much risk to life on terrestrial planets which would be expected to develop later. In the final 20\% of cases the result will be a flux of impactors 5-10 times greater than that on Earth that can persist for several Gyr, quite detrimental to the development of life on the planet.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Fly-by encounters between two planetary systems II: Exploring the interactions of diverse planetary system architectures

Planetary systems formed in clusters may be subject to stellar encounter flybys. Here we create a diverse range of representative planetary systems with different orbital scales and planets' masses and examine encounters between them in a typical open cluster. We first explore the close-in multi-super earth systems $\lesssim0.1$ au. They are resistant to flybys in that only ones inside a few au can destabilise a planet or break the resonance between such planets. But these systems may capture giant planets onto wide orbits from the intruding star during distant flybys. If so, the original close-in small planets' orbits may be tilted together through Kozai--Lidov mechanism, forming a ``cold'' system that is significantly inclined against the equator of the central host. Moving to the intermediately-placed planets around solar-like stars, we find that the planets' mass gradient governs the systems' long-term evolution post-encounter: more massive planets have better chances to survive. Also, a system's angular momentum deficit, a quantity describing how eccentric/inclined the orbits are, measured immediately after the encounter, closely relates to the longevity of the systems -- whether or not and when the systems turn unstable in the ensuing evolution millions of years post-encounter. We compare the orbits of the surviving planets in the unstable systems through (1) the immediate consequence of the stellar fly or (2) internal interplanetary scattering long post-encounter and find that those for the former are systematically colder. Finally, we show that massive wide-orbit multi-planet systems like that of HR 8799 can be easily disrupted and encounters at a few hundreds of au suffice.

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

Two intermediate-mass transiting brown dwarfs from the TESS mission

We report the discovery of two intermediate-mass brown dwarfs (BDs), TOI-569b and TOI-1406b, from NASA's Transiting Exoplanet Survey Satellite mission. TOI-569b has an orbital period of $P = 6.55604 \pm 0.00016$ days, a mass of $M_b = 64.1 \pm 1.9 M_J$, and a radius of $R_b = 0.75 \pm 0.02 R_J$. Its host star, TOI-569, has a mass of $M_\star = 1.21 \pm 0.03 M_\odot$, a radius of $R_\star = 1.47 \pm 0.03 R_\odot$, $\rm [Fe/H] = +0.29 \pm 0.09$ dex, and an effective temperature of $T_{\rm eff} = 5768 \pm 110K$. TOI-1406b has an orbital period of $P = 10.57415 \pm 0.00063$ days, a mass of $M_b =46.0 \pm 2.7 M_J$, and a radius of $R_b = 0.86 \pm 0.03 R_J$. The host star for this BD has a mass of $M_\star =1 .18 \pm 0.09 M_\odot$, a radius of $R_\star = 1.35 \pm 0.03 R_\odot$, $ \rm [Fe/H] = -0.08 \pm 0.09$ dex and an effective temperature of $T_{\rm eff} = 6290 \pm 100K$. Both BDs are in circular orbits around their host stars and are older than 3 Gyr based on stellar isochrone models of the stars. TOI-569 is one of two slightly evolved stars known to host a transiting BD (the other being KOI-415). TOI-1406b is one of three known transiting BDs to occupy the mass range of $40-50 M_J$ and one of two to have a circular orbit at a period near 10 days (with the first being KOI-205b).Both BDs have reliable ages from stellar isochrones in addition to their well-constrained masses and radii, making them particularly valuable as tests for substellar isochrones in the BD mass-radius diagram.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

On the inclinations of the Jupiter Trojans

Jupiter Trojans are are characterized by dark photometric colors, high inclinations and an asymmetry in number of bodies between the two swarms. Different models have been proposed to explain the high inclination of the Trojans and to interpret their relation with the TNOs, but none of them can also satisfactorily explain the asymmetry ratio. Recently it has been found that the asymmetry can arise if Jupiter has migrated inwards by at least a few au during its growth. The asymmetry and the dark colors of the Trojans are natural outcomes of this model, but simulations with massless unperturbed disc particles led to a flat distribution of the Trojan inclinations and a final total mass that was 3-4 orders of magnitude larger than the current one. In our work, we investigate the possible origin of the peculiar inclination distribution of the Trojans in the scenario where Jupiter migrates inwards. We analyze different possibilities: (a) the secular evolution of an initially flat Trojan population, (b) the presence of planetary embryos among the Trojans and (c) capture of the Trojans from a pre-stirred planetesimal population. We find that the secular evolution of the Trojans and Saturn do not affect the Trojan inclinations appreciably, nor is there any significant mass depletion. Embryos embedded in the swarms, in contrast, can stir the Trojan inclinations and can also deplete the swarms efficiently, but it turns out that it is very difficult to get rid of all of the massive bodies. We propose that the disc where Jupiter's core was forming was already excited by the presence of other embryos competing in Jupiter's core's feeding zone. We show that the trapped Trojans preserve their high inclination through the gas phase of the disc and that Saturn's perturbations are more effective on highly inclined Trojans, leading to a lower capture efficiency and to a substantial depletion of the swarms.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Resilient habitability of nearby exoplanet systems

We investigate the possibility of finding Earth-like planets in the habitable zone of 34 nearby FGK-dwarfs, each known to host one giant planet exterior to their habitable zone detected by RV. First we simulate the dynamics of the planetary systems in their present day configurations and determine the fraction of stable planetary orbits within their habitable zones. Then, we postulate that the eccentricity of the giant planet is a result of an instability in their past during which one or more other planets were ejected from the system. We simulate these scenarios and investigate whether planets orbiting in the habitable zone survive the instability. Explicitly we determine the fraction of test particles, originally found in the habitable zone, which remain in the habitable zone today. We label this fraction the resilient habitability of a system. We find that for most systems the probability of planets existing [or surviving] on stable orbits in the habitable zone becomes significantly smaller when we include a phase of instability in their history. We present a list of candidate systems with high resilient habitability for future observations. These are: HD 95872, HD 154345, HD 102843, HD 25015, GJ 328, HD 6718 and HD 150706. The known planets in the last two systems have large observational uncertainties on their eccentricities, which propagate into large uncertainties on their resilient habitability. Further observational constraints of these two eccentriciti

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Super-Earth ingestion can explain the anomalously high metal abundances of M67 Y2235

We investigate the hypothesis that ingestion of a terrestrial or super-Earth planet could cause the anomalously high metal abundances seen in a turn-off star in the open cluster M67, when compared to other turn-off stars in the same cluster. We show that the mass in convective envelope of the star is likely only $3.45\,\times 10^{-3}\,{\rm M}_\odot$, and hence $5.2\,{\rm M}_\oplus$ of rock is required to obtain the observed 0.128 dex metal enhancement. Rocky planets dissolve entirely in the convective envelope if they enter it with sufficiently tangential orbits: we find that the critical condition for dissolution is that the planet's radial speed must be less than 40% of its total velocity at the stellar surface; or, equivalently, the impact parameter must be greater than about 0.9. We model the delivery of rocky planets to the stellar surface both by planet-planet scattering in a realistic multi-planet system, and by Lidov-Kozai cycles driven by a more massive planetary or stellar companion. In both cases almost all planets that are ingested arrive at the star on grazing orbits and hence will dissolve in the surface convection zone. We conclude that super-Earth ingestion is a good explanation for the metal enhancement in M67 Y2235, and that a high-resolution spectroscopic survey of stellar abundances around the turn-off and main sequence of M67 has the potential to constrain the frequency of late-time dynamical instability in planetary systems.

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

Is there an exoplanet in the Solar System?

We investigate the prospects for the capture of the proposed Planet 9 from other stars in the Sun's birth cluster. Any capture scenario must satisfy three conditions: the encounter must be more distant than ~150 au to avoid perturbing the Kuiper belt; the other star must have a wide-orbit planet (a>~100au); the planet must be captured onto an appropriate orbit to sculpt the orbital distribution of wide-orbit Solar System bodies. Here we use N-body simulations to show that these criteria may be simultaneously satisfied. In a few percent of slow close encounters in a cluster, bodies are captured onto heliocentric, Planet 9-like orbits. During the ~100 Myr cluster phase, many stars are likely to host planets on highly-eccentric orbits with apastron distances beyond 100 au if Neptune-sized planets are common and susceptible to planet--planet scattering. While the existence of Planet 9 remains unproven, we consider capture from one of the Sun's young brethren a plausible route to explain such an object's orbit. Capture appears to predict a large population of Trans-Neptunian Objects (TNOs) whose orbits are aligned with the captured planet, and we propose that different formation mechanisms will be distinguishable based on their imprint on the distribution of TNOs.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

The origin of the eccentricity of the hot Jupiter in CI Tau

Following the recent discovery of the first radial velocity planet in a star still possessing a protoplanetary disc (CI Tau), we examine the origin of the planet's eccentricity (e $\sim 0.3$). We show through long timescale ($10^5$ orbits) simulations that the planetary eccentricity can be pumped by the disc, even when its local surface density is well below the threshold previously derived from short timescale integrations. We show that the disc may be able to excite the planet's orbital eccentricity in $<$ a Myr for the system parameters of CI Tau. We also perform two planet scattering experiments and show that alternatively the observed planet may plausibly have acquired its eccentricity through dynamical scattering of a migrating lower mass planet, which has either been ejected from the system or swallowed by the central star. In the latter case the present location and eccentricity of the observed planet can be recovered if it was previously stalled within the disc's magnetospheric cavity.

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

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

The destruction of inner planetary systems during high-eccentricity migration of gas giants

Hot Jupiters are giant planets on orbits a few hundredths of an AU. They do not share their system with low-mass close-in planets, despite these latter being exceedingly common. Two migration channels for hot Jupiters have been proposed: through a protoplanetary gas disc or by tidal circularisation of highly-eccentric planets. We show that highly-eccentric giant planets that will become hot Jupiters clear out any low-mass inner planets in the system, explaining the observed lack of such companions to hot Jupiters. A less common outcome of the interaction is that the giant planet is ejected by the inner planets. Furthermore, the interaction can implant giant planets on moderately-high eccentricities at semimajor axes $<1$ AU, a region otherwise hard to populate. Our work supports the hypothesis that most hot Jupiters reached their current orbits following a phase of high eccentricity, possibly excited by other planetary or stellar companions.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Hot Jupiters and Cool Stars

Close-in planets are in jeopardy as their host stars evolve off the main sequence to the subgiant and red giant phases. In this paper, we explore the influences of the stellar mass (in the range 1.5--2\Mso ), mass-loss prescription, planet mass (from Neptune up to 10 Jupiter masses), and eccentricity, on the orbital evolution of planets as their parent stars evolve to become subgiants and Red Giants. We find that planet engulfment during the Red Giant Branch is not very sensitive to the stellar mass or mass-loss rates adopted in the calculations, but quite sensitive to the planetary mass. The range of initial separations for planet engulfment increases with decreasing mass-loss rates or stellar mass and increasing planetary masses. Regarding the planet's orbital eccentricity, we find that as the star evolves into the red giant phase, stellar tides start to dominate over planetary tides. As a consequence, a transient population of moderately eccentric close-in Jovian planets is created, that otherwise would have been expected to be absent from main sequence stars. We find that very eccentric and distant planets do not experience much eccentricity decay, and that planet engulfment is primarily determined by the pericenter distance and the maximum stellar radius.

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

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
preprint2011arXiv

Dependence of a planet's chaotic zone on particle eccentricity: the shape of debris disc inner edges

The orbit of a planet is surrounded by a chaotic zone wherein nearby particles' orbits are chaotic and unstable. Wisdom (1980) showed that the chaos is driven by the overlap of mean motion resonances which occurs within a distance (da/a)_chaos = 1.3 mu^2/7 of the planet's orbit. However, the width of mean motion resonances grows with the particles' eccentricity, which will increase the width of the chaotic zone at higher eccentricities. Here we investigate the width of the chaotic zone using the iterated encounter map and N-body integrations. We find that the classical prescription works well for particles on low-eccentricity orbits. However, above a critical eccentricity, dependent upon the mass of the planet, the width of the chaotic zone increases with eccentricity. An extension of Wisdom's analytical arguments then shows that, above the critical eccentricity, the chaotic zone width is given by (da/a)_chaos = 1.8 e^1/5 mu^1/5, which agrees well with the encounter map results. The critical eccentricity is given by e_cr = 0.21 mu^3/7. This extended chaotic zone results in a larger cleared region when a planet sculpts the inner edge of a debris disc comprised of eccentric planetesimals. Hence, the planet mass estimated from the classical chaotic zone may be erroneous. We apply this result to the HR8799 system, showing that the masses of HR8799b inferred from the truncation of the disc may vary by up to 50% depending on the disc particles' eccentricities. With a disc edge at 90AU, the necessary mass of planet b to cause the truncation would be 8-10 Jovian masses if the disc particles have low eccentricities (<0.02), but only 4-8 Jovian masses if the disc particles have higher eccentricities. Our result also has implications for the ability of a planet to feed material into an inner system, a process which may explain metal pollution in White Dwarf atmospheres.

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

A general model of resonance capture in planetary systems: First and second order resonances

Mean motion resonances are a common feature of both our own Solar System and of extrasolar planetary systems. Bodies can be trapped in resonance when their orbital semi-major axes change, for instance when they migrate through a protoplanetary disc. We use a Hamiltonian model to thoroughly investigate the capture behaviour for first and second order resonances. Using this method, all resonances of the same order can be described by one equation, with applications to specific resonances by appropriate scaling. We focus on the limit where one body is a massless test particle and the other a massive planet. We quantify how the the probability of capture into a resonance depends on the relative migration rate of the planet and particle, and the particle's eccentricity. Resonant capture fails for high migration rates, and has decreasing probability for higher eccentricities. More massive planets can capture particles at higher eccentricities and migration rates. We also calculate libration amplitudes and the offset of the libration centres for captured particles, and the change in eccentricity if capture does not occur. Libration amplitudes are higher for larger initial eccentricity. The model allows for a complete description of a particle's behaviour as it successively encounters several resonances. We discuss implications for several scenarios: (i) Planet migration through gas discs trapping other planets or planetesimals in resonances. (ii) Planet migration through a debris disc. (iii) Dust migration through PR drag. The Hamiltonian model will allow quick interpretation of the resonant properties of extrasolar planets and Kuiper Belt Objects, and will allow synthetic images of debris disc structures to be quickly generated, which will be useful for predicting and interpreting disc images made with ALMA, Darwin/TPF or similar missions. [Abridged]

0 citations0 reviewsNo code linkedOpen access