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.

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.

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.

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.

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

The atmospheres of close-in planets are strongly influenced by mass loss driven by the high-energy (X-ray and extreme ultraviolet, EUV) irradiation of the host star, particularly during the early stages of evolution. We recently developed a framework to exploit this connection and enable us to recover the past evolution of the stellar high-energy emission from the present-day properties of its planets, if the latter retains some remnants of their primordial hydrogen-dominated atmospheres. Furthermore, the framework can also provide constraints on planetary initial atmospheric mass fractions. The constraints on the output parameters improve when more planets can be simultaneously analysed. This makes the Kepler-11 system, which hosts six planets with bulk densities between 0.66 and 2.45g cm^{-3}, an ideal target. Our results indicate that the star has likely evolved as a slow rotator (slower than 85\% of the stars with similar masses), corresponding to a high-energy emission at 150 Myr of between 1-10 times that of the current Sun. We also constrain the initial atmospheric mass fractions for the planets, obtaining a lower limit of 4.1% for planet c, a range of 3.7-5.3% for planet d, a range of 11.1-14% for planet e, a range of 1-15.6% for planet f, and a range of 4.7-8.7% for planet g assuming a disc dispersal time of 1 Myr. For planet b, the range remains poorly constrained. Our framework also suggests slightly higher masses for planets b, c, and f than have been suggested based on transit timing variation measurements. We coupled our results with published planet atmosphere accretion models to obtain a temperature (at 0.25 AU, the location of planet f) and dispersal time of the protoplanetary disc of 550 K and 1 Myr, although these results may be affected by inconsistencies in the adopted system parameters.

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.

Statistical analyses from exoplanet surveys around low-mass stars indicate that super-Earth and Neptune-mass planets are more frequent than gas giants around such stars, in agreement with core accretion theory of planet formation. Using precise radial velocities derived from visual and near-infrared spectra, we report the discovery of a giant planet with a minimum mass of 0.46 Jupiter masses in an eccentric 204-day orbit around the very low-mass star GJ 3512. Dynamical models show that the high eccentricity of the orbit is most likely explained from planet-planet interactions. The reported planetary system challenges current formation theories and puts stringent constraints on the accretion and migration rates of planet formation and evolution models, indicating that disc instability may be more efficient in forming planets than previously thought.

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.

Stars formed in clusters can encounter other stars at close distances. In typical open clusters in the Solar neighbourhood containing hundreds or thousands of member stars, ten to twenty per cent of Solar-mass member stars are expected to encounter another star at distances closer than 100 au. These close encounters strongly perturb the planetary systems, directly causing ejection of planets or their capture by the intruding star, as well as exciting the orbits. Using extensive $N$-body simulations, we study such fly-by encounters between two Solar System analogues, each with four giant planets from Jupiter to Neptune. We quantify the rates of loss and capture immediately after the encounter, e.g., the Neptune analogue is lost in one in four encounters within 100 au, and captured by the flying-by star in one in twelve encounters. We then perform long-term (up to 1 Gyr) simulations investigating the ensuing post-encounter evolution. We show that large numbers of planets are removed from systems due to planet--planet interactions and that captured planets further enhance the system instability. While encounters can initially leave a planetary system containing more planets by inserting additional ones, the long-term instability causes a net reduction in planet number. A captured planet ends up on a retrograde orbit in half of the runs in which it survives for 1 Gyr; also, a planet bound to its original host star but flipped during the encounter may survive. Thus, encounters between planetary systems are a channel to create counter-rotating planets, This would happen in around 1% of systems, and such planets are potentially detectable through astrometry or direct imaging.

Our aim is to investigate the nature and formation of brown dwarfs by adding a new well-characterised object to the small sample of less than 20 transiting brown dwarfs. One brown dwarf candidate was found by the KESPRINT consortium when searching for exoplanets in the K2 space mission Campaign 16 field. We combined the K2 photometric data with a series of multi-colour photometric observations, imaging and radial velocity measurements to rule out false positive scenarios and to determine the fundamental properties of the system. We report the discovery and characterisation of a transiting brown dwarf in a 5.17 day eccentric orbit around the slightly evolved F7V star EPIC 212036875. We find a stellar mass of 1.15+/-0.08 M$_\odot$, a stellar radius of 1.41+/-0.05 R$_\odot$, and an age of 5.1+/-0.9 Gyr. The mass and radius of the companion brown dwarf are 51+/-2 MJ and 0.83+/-0.03 RJ, respectively, corresponding to a mean density of 108+15-13 g cm-3. EPIC 212036875 b is a rare object that resides in the brown dwarf desert. In the mass-density diagram for planets, brown dwarfs and stars, we find that all giant planets and brown dwarfs follow the same trend from ~0.3 MJ to the turn-over to hydrogen burning stars at ~73 MJ. EPIC 212036875 b falls close to the theoretical model for mature H/He dominated objects in this diagram as determined by interior structure models, as well as the empirical fit. We argue that EPIC 212036875 b formed via gravitational disc instabilities in the outer part of the disc, followed by a quick migration. Orbital tidal circularisation may have started early in its history for a brief period when the brown dwarf's radius was larger. The lack of spin-orbit synchronisation points to a weak stellar dissipation parameter which implies a circularisation timescale of >23 Gyr, or suggests an interaction between the magnetic and tidal forces of the star and the brown dwarf.

Many white dwarf stars show signs of having accreted smaller bodies, implying that they may host planetary systems. A small number of these systems contain gaseous debris discs, visible through emission lines. We report a stable 123.4min periodic variation in the strength and shape of the CaII emission line profiles originating from the debris disc around the white dwarf SDSSJ122859.93+104032.9. We interpret this short-period signal as the signature of a solid body held together by its internal strength.

Pebble accretion is an efficient mechanism able to build up the core of the giant planets within the lifetime of the protoplanetary disc gas-phase. The core grows via this process until the protoplanet reaches its pebble isolation mass and starts to accrete gas. During the growth, the protoplanet undergoes a rapid, large-scale, inward migration due to the interactions with the gaseous protoplanetary disc. In our work, we investigate how this early migration would have affected the minor body populations in our solar system. In particular, we focus on the Jupiter Trojans and the Hildas asteroids. We found that a massive and eccentric Hilda group is captured during the migration from a region between 5 and 8 au and subsequently depleted during the late instability of the giant planets. Our simulations also show that inward migration of the giant planets always produces a Jupiter Trojans' leading swarm more populated than the trailing one, with a ratio comparable to the current observed Trojan asymmetry ratio. The in situ formation of Jupiter, on the other hand, produces symmetric leading/trailing swarms. The reason for the asymmetry is the relative drift between the migrating planet and the particles in the coorbital resonance. The capture happens during the growth of Jupiter's core and Trojan asteroids are afterwards carried along during the giant planet's migration to their final orbits. The asymmetry and eccentricity of the captured Trojans correspond well to observations, but their inclinations are near zero and their total mass is 3-4 orders of magnitude higher than the current population. Future modelling will be needed to understand whether the dynamical evolution of the Trojans over billions of years will raise the inclinations and deplete the masses to observed values.

WD 1145+017 is currently the only white dwarf known to exhibit periodic transits of planetary debris as well as absorption lines from circumstellar gas. We present the first simultaneous fast optical spectrophotometry and broad-band photometry of the system, obtained with the Gran Telescopio Canarias (GTC) and the Liverpool Telescope (LT), respectively. The observations spanned $5.5$ h, somewhat longer than the $4.5$-h orbital period of the debris. Dividing the GTC spectrophotometry into five wavelength bands reveals no significant colour differences, confirming grey transits in the optical. We argue that absorption by an optically thick structure is a plausible alternative explanation for the achromatic nature of the transits that can allow the presence of small-sized ($\sim\mu$m) particles. The longest ($87$ min) and deepest ($50$ per cent attenuation) transit recorded in our data exhibits a complex structure around minimum light that can be well modelled by multiple overlapping dust clouds. The strongest circumstellar absorption line, Fe II $\lambda$5169, significantly weakens during this transit, with its equivalent width reducing from a mean out-of-transit value of $2$ \AA\ to $1$ \AA\ in-transit, supporting spatial correlation between the circumstellar gas and dust. Finally, we made use of the Gaia Data Release 2 and archival photometry to determine the white dwarf parameters. Adopting a helium-dominated atmosphere containing traces of hydrogen and metals, and a reddening $E(B-V)=0.01$ we find $T_\mathrm{eff}=15\,020 \pm 520$ K, $\log g=8.07\pm0.07$, corresponding to $M_\mathrm{WD}=0.63\pm0.05\ \mbox{$\mathrm{M}_{\odot}$}$ and a cooling age of $224\pm30$ Myr.

Context. Gaia DR2 offers unparalleled precision on stars' parallaxes and proper motions. This allows the prediction of microlensing events for which the lens stars (and any planets they possess) are nearby and may be well studied and characterised. Aims. We identify a number of potential microlensing events that will occur before the year 2035.5, 20 years from the Gaia DR2 reference epoch. Methods. We query Gaia DR2 for potential lenses within 100 pc, extract parallaxes and proper motions of the lenses and background sources, and identify potential lensing events. We estimate the lens masses from Priam effective temperatures, and use these to calculate peak magnifications and the size of the Einstein radii relative to the lens stars' habitable zones. Results. We identify 7 future events with a probability > 10% of an alignment within one Einstein radius. Of particular interest is DR2 5918299904067162240 (WISE J175839.20-583931.6), magnitude G = 14.9, which will lens a G = 13.9 background star in early 2030, with a median 23% net magnification. Other pairs are typically fainter, hampering characterisation of the lens (if the lens is faint) or the ability to accurately measure the magnification (if the source is much fainter than the lens). Of timely interest is DR2 4116504399886241792 (2MASS J17392440-2327071), which will lens a background star in July 2020, albeit with weak net magnification (0.03%). Median magnifications for the other 5 high-probability events range from 0.3% to 5.3%. The Einstein radii for these lenses are 1-10 times the radius of the habitable zone, allowing these lensing events to pick out cold planets around the ice line, and filling a gap between transit and current microlensing detections of planets around very low-mass stars. Conclusions. We provide a catalogue of the predicted events to aid future characterisation efforts... [abridged]

Planet-planet collisions are a common outcome of instability in systems of transiting planets close to the star, as well as occurring during in-situ formation of such planets from embryos. Previous N-body studies of instability amongst transiting planets have assumed that collisions result in perfect merging. Here, we explore the effects of implementing a more realistic collision prescription on the outcomes of instability and in-situ formation at orbital radii of a few tenths of an au. There is a strong effect on the outcome of the growth of planetary embryos, so long as the debris thrown off in collisions is rapidly removed from the system (which happens by collisional processing to dust, and then removal by radiation forces) and embryos are small ($< 0.1\mathrm{\,M}_\oplus$). If this is the case, then systems form fewer detectable ($\gtrsim1\mathrm{\,M}_\oplus$) planets than systems evolved under the assumption of perfect merging in collisions. This provides some contribution to the "Kepler Dichotomy": the observed over-abundance of single-planet systems. The effects of changing the collision prescription on unstable mature systems of super-Earths are less pronounced. Perfect mergers only account for a minority of collision outcomes in such systems, but most collisions resulting in mass loss are grazing impacts in which only a few per cent. of mass is lost. As a result, there is little impact on the final masses and multiplicities of the systems after instability when compared to systems evolved under the assumption that collisions always result in perfect merging.

Differences in the elemental abundances of planet hosting stars in binary systems can give important clues and constraints about planet formation and evolution. In this study we performed a high-precision, differential elemental abundance analysis of a wide binary system, HD 80606/80607, based on high-resolution, high signal-to-noise ratio Keck/HIRES spectra. HD 80606 is known to host a four Jupiter mass giant planet while no planet has yet been detected around HD 80607. We determined stellar parameters as well as abundances for 23 elements for these two stars with extremely high precision. Our main results are: (i) we confirmed that the two components share very similar chemical compositions, but HD 80606 is marginally more metal-rich than HD 80607 with an average difference of +0.013 $\pm$ 0.002 dex ($\sigma$ = 0.009 dex) and (ii) there is no obvious trend between abundance differences and condensation temperature. Assuming this binary formed from material with the same chemical composition, it is difficult to understand how giant planet formation could produce the present-day photospheric abundances of the elements we measure. We can not exclude the possibility that HD 80606 might have accreted about 2.5 to 5 $M_{\rm Earth}$ material onto its surface, possibly from a planet destabilised by the known highly-eccentric giant.

At least 25% of white dwarfs show atmospheric pollution by metals, sometimes accompanied by detectable circumstellar dust/gas discs or (in the case of WD 1145+017) transiting disintegrating asteroids. Delivery of planetesimals to the white dwarf by orbiting planets is a leading candidate to explain these phenomena. Here, we study systems of planets and planetesimals undergoing planet-planet scattering triggered by the star's post-main sequence mass loss, and test whether this can maintain high rates of delivery over the several Gyr that they are observed. We find that low-mass planets (Earth to Neptune mass) are efficient deliverers of material and can maintain the delivery for Gyr. Unstable low-mass planetary systems reproduce the observed delayed onset of significant accretion, as well as the slow decay in accretion rates at late times. Higher-mass planets are less efficient, and the delivery only lasts a relatively brief time before the planetesimal populations are cleared. The orbital inclinations of bodies as they cross the white dwarf's Roche limit are roughly isotropic, implying that significant collisional interactions of asteroids, debris streams and discs can be expected. If planet-planet scattering is indeed responsible for the pollution of white dwarfs, many such objects, and their main-sequence progenitors, can be expected to host (currently undetectable) super-Earth planets on orbits of several au and beyond.

Unexpected clustering in the orbital elements of minor bodies beyond the Kuiper belt has led to speculations that our solar system actually hosts nine planets, the eight established plus a hypothetical "Planet Nine". Several recent studies have shown that a planet with a mass of about 10 Earth masses on a distant eccentric orbit with perihelion far beyond the Kuiper belt could create and maintain this clustering. The evolutionary path resulting in an orbit such as the one suggested for Planet Nine is nevertheless not easily explained. Here we investigate whether a planet scattered away from the giant-planet region could be lifted to an orbit similar to the one suggested for Planet Nine through dynamical friction with a cold, distant planetesimal belt. Recent simulations of planetesimal formation via the streaming instability suggest that planetesimals can readily form beyond 100au. We explore this circularisation by dynamical friction with a set of numerical simulations. We find that a planet that is scattered from the region close to Neptune onto an eccentric orbit has a 20-30% chance of obtaining an orbit similar to that of Planet Nine after 4.6Gyr. Our simulations also result in strong or partial clustering of the planetesimals; however, whether or not this clustering is observable depends on the location of the inner edge of the planetesimal belt. If the inner edge is located at 200au the degree of clustering amongst observable objects is significant.

Radial velocity observations from three instruments reveal the presence of a 4 M_jup planet candidate orbiting the K giant HD 76920. HD 76920b has an orbital eccentricity of 0.856$\pm$0.009, making it the most eccentric planet known to orbit an evolved star. There is no indication that HD 76920 has an unseen binary companion, suggesting a scattering event rather than Kozai oscillations as a probable culprit for the observed eccentricity. The candidate planet currently approaches to about four stellar radii from its host star, and is predicted to be engulfed on a $\sim$100 Myr timescale due to the combined effects of stellar evolution and tidal interactions.

The star EPIC 210894022 has been identified from a light curve acquired through the K2 space mission as possibly orbited by a transiting planet. Our aim is to confirm the planetary nature of the object and derive its fundamental parameters. We combine the K2 photometry with reconnaissance spectroscopy and radial velocity (RV) measurements obtained using three separate telescope and spectrograph combinations. The spectroscopic synthesis package SME has been used to derive the stellar photospheric parameters that were used as input to various stellar evolutionary tracks in order to derive the parameters of the system. The planetary transit was also validated to occur on the assumed host star through adaptive imaging and statistical analysis. The star is found to be located in the background of the Hyades cluster at a distance at least 4 times further away from Earth than the cluster itself. The spectrum and the space velocities of EPIC 210894022 strongly suggest it to be a member of the thick disk population. We find that the star is a metal poor ([Fe/H]=-0.53+/-0.05 dex) and alpha-rich somewhat evolved solar-like object of spectral type G3 with Teff=5730+/-50 K, logg=4.15+/-0.1 (cgs), radius of 1.3+/-0.1 R_Sun, and mass of 0.88+/-0.02 M_Sun. The RV detection together with the imaging confirms with a high level of significance that the transit signature is caused by a super-Earth orbiting the star EPIC 210894022. We measure a mass of 8.6+/-3.9 M_Earth and a radius of 1.9+/-0.2 R_Earth. A second more massive object with a period longer than about 120 days is indicated by a long term linear acceleration. With an age of > 10 Gyrs this system is one of the oldest where planets is hitherto detected. Further studies of this planetary system is important since it contains information about the planetary formation process during a very early epoch of the history of our Galaxy.

We study how close-in systems such as those detected by Kepler are affected by the dynamics of bodies in the outer system. We consider two scenarios: outer systems of giant planets potentially unstable to planet--planet scattering, and wide binaries that may be capable of driving Kozai or other secular variations of outer planets' eccentricities. Dynamical excitation of planets in the outer system reduces the multiplicity of Kepler-detectable planets in the inner system in $\sim20-25\%$ of our systems. Accounting for the occurrence rates of wide-orbit planets and binary stars, $\approx18\%$ of close-in systems could be destabilised by their outer companions in this way. This provides some contribution to the apparent excess of systems with a single transiting planet compared to multiple, however, it only contributes at most $25\%$ of the excess. The effects of the outer dynamics can generate systems similar to Kepler-56 (two coplanar planets significantly misaligned with the host star) and Kepler-108 (two significantly non-coplanar planets in a binary). We also identify three pathways to the formation of eccentric warm Jupiters resulting from the interaction between outer and inner systems: direct inelastic collision between an eccentric outer and an inner planet, secular eccentricity oscillations that may "freeze out" when scattering resolves in the outer system; and scattering in the inner system followed by "uplift", where inner planets are removed by interaction with the outer planets. In these scenarios, the formation of eccentric warm Jupiters is a signature of a past history of violent dynamics among massive planets beyond $\sim1$ au.

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.

We present an intensive monitoring of high-resolution spectra of the Ca {\sc ii} K line in the A7IV shell star $\Phi$ Leo at very short (minutes, hours), short (night to night), and medium (weeks, months) timescales. The spectra show remarkable variable absorptions on timescales of hours, days, and months. The characteristics of these sporadic events are very similar to most that are observed toward the debris disk host star $\beta$ Pic, which are commonly interpreted as signs of the evaporation of solid, comet-like bodies grazing or falling onto the star. Therefore, our results suggest the presence of solid bodies around $\Phi$ Leo. To our knowledge, with the exception of $\beta$ Pic, our monitoring has the best time resolution at the mentioned timescales for a star with events attributed to exocomets. Assuming the cometary scenario and considering the timescales of our monitoring, our results indicate that $\Phi$ Leo presents the richest environment with comet-like events known to date, second only to $\beta$ Pic.

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.

HR 8799 is a star accompanied by four massive planets on wide orbits. The observed planetary configuration has been shown to be unstable on a timescale much shorter than the estimated age of the system (~ 30 Myr) unless the planets are locked into mean motion resonances. This condition is characterised by small-amplitude libration of one or more resonant angles that stabilise the system by preventing close encounters. We simulate planetary systems similar to the HR 8799 planetary system, exploring the parameter space in separation between the orbits, planetary masses and distance from the Sun to the star. We find systems that look like HR 8799 and remain stable for longer than the estimated age of HR 8799. None of our systems are forced into resonances. We find, with nominal masses and in a narrow range of orbit separations, that 5 of 100 systems match the observations and lifetime. Considering a broad range of orbit separations, we find 12 of 900 similar systems. The systems survive significantly longer because of their slightly increased initial orbit separations compared to assuming circular orbits from the observed positions. A small increase in separation leads to a significant increase in survival time. The low eccentricity the orbits develop from gravitational interaction is enough for the planets to match the observations. With lower masses, but still comfortably within the estimated planet mass uncertainty, we find 18 of 100 matching and long-lived systems in a narrow orbital separation range. In the broad separation range, we find 82 of 900 matching systems. Our results imply that the planets in the HR 8799 system do not have to be in strong mean motion resonances.

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.

We explore the role of dynamics in shaping planetary system multiplicities, focussing on two particular problems. (1) We propose that the lack of close-in super-Earths in hot Jupiter systems is a signature of the migration history of the hot Jupiters and helps to discriminate between different mechanisms of migration. We present N-body simulations of dynamical migration scenarios where proto-hot Jupiters are excited to high eccentricities prior to tidal circularisation and orbital decay. We show that in this scenario, the eccentric giant planet typically destroys planets in the inner system, in agreement with the observed lack of close super-Earth companions to hot Jupiters. (2) We explore the role of the dynamics of outer systems in affecting the multiplicities of close-in systems such as those discovered by Kepler. We consider specifically the effects of planet--planet scattering and Kozai perturbations on an exterior giant planet on the architecture of the inner system, and evaluate the ability of such scenarios to reduce the inner system's multiplicity and contribute to the observed excess of single Kepler planets.

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.

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.

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.

Dust grains migrating under Poynting-Robertson drag may be trapped in mean-motion resonances with planets. Such resonantly trapped grains are observed in the solar system. In extrasolar systems, the exozodiacal light produced by dust grains is expected to be a major obstacle to future missions attempting to directly image terrestrial planets. The patterns made by resonantly trapped dust, however, can be used to infer the presence of planets, and the properties of those planets, if the capture and evolution of the grains can be modelled. This has been done with N-body methods, but such methods are computationally expensive, limiting their usefulness when considering large, slowly evolving grains, and for extrasolar systems with unknown planets and parent bodies, where the possible parameter space for investigation is large. In this work, we present a semi-analytic method for calculating the capture and evolution of dust grains in resonance, which can be orders of magnitude faster than N-body methods. We calibrate the model against N-body simulations, finding excellent agreement for Earth to Neptune mass planets, for a variety of grain sizes, initial eccentricities, and initial semimajor axes. We then apply the model to observations of dust resonantly trapped by the Earth. We find that resonantly trapped, asteroidally produced grains naturally produce the `trailing blob' structure in the zodiacal cloud, while to match the intensity of the blob, most of the cloud must be composed of cometary grains, which owing to their high eccentricity are not captured, but produce a smooth disk.

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.

We report the discovery of WTS-2 b, an unusually close-in 1.02-day hot Jupiter (Mp=1.12MJ, Rp=1.363RJ) orbiting a K2V star, which has a possible gravitationally-bound M-dwarf companion at 0.6 arcsec separation contributing ~20 percent of the total flux in the observed J-band light curve. The planet is only 1.5 times the separation from its host star at which it would be destroyed by Roche lobe overflow, and has a predicted remaining lifetime of just ~40 Myr, assuming a tidal dissipation quality factor of Q'*=10^6. Q'* is a key factor in determining how frictional processes within a host star affect the orbital evolution of its companion giant planets, but it is currently poorly constrained by observations. We calculate that the orbital decay of WTS-2 b would correspond to a shift in its transit arrival time of T_shift~17 seconds after 15 years assuming Q'*=10^6. A shift less than this would place a direct observational constraint on the lower limit of Q'* in this system. We also report a correction to the previously published expected T_shift for WASP-18 b, finding that T_shift=356 seconds after 10 years for Q'*=10^6, which is much larger than the estimated 28 seconds quoted in WASP-18 b discovery paper. We attempted to constrain Q'* via a study of the entire population of known transiting hot Jupiters, but our results were inconclusive, requiring a more detailed treatment of transit survey sensitivities at long periods. We conclude that the most informative and straight-forward constraints on Q'* will be obtained by direct observational measurements of the shift in transit arrival times in individual hot Jupiter systems. We show that this is achievable across the mass spectrum of exoplanet host stars within a decade, and will directly probe the effects of stellar interior structure on tidal dissipation.

In the last few years, a number of planets have been proposed to orbit several post main-sequence binary star systems on the basis of observed variations in the timing of eclipses between the binary components. A common feature of these planet candidates is that the best-fit orbits are often highly eccentric, such that the multiple planet systems proposed regularly feature mutually crossing orbits - a scenario that almost always leads to unstable planetary systems. In this work, we present the results of dynamical studies of all multiple-planet systems proposed to orbit these highly evolved binary stars, finding that most do not stand up to dynamical scrutiny. In one of the potentially stable cases (the NN Serpentis 2-planet system), we consider the evolution of the binary star system, and show that it is highly unlikely that planets could survive from the main sequence to obtain their current orbits - again casting doubt on the proposed planets. We conclude by considering alternative explanations for the observed variation in eclipse timings for these systems.

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.

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.

In recent years, a number of planetary systems have been proposed to orbit evolved binary star systems. The presence of planets is invoked to explain observed variations in the timing of mutual eclipses between the primary and secondary components of the binary star system. The planets recently proposed orbiting the cataclysmic variable system QS Virginis are the latest in this on-going series of "extreme planets". The two planets proposed to orbit QS Virginis would move on mutually crossing orbits - a situation that is almost invariably unstable on very short timescales. In this work, we present the results of a detailed dynamical study of the orbital evolution of the two proposed planets, revealing that they are dynamically unstable on timescales of less than one thousand years across the entire range of orbital elements that provide a plausible fit to the observational data, and regardless of their mutual orbital inclination. We conclude that the proposed planets around the cataclysmic variable QS Virginis simply cannot exist.

(abridged) Infrared excesses associated with debris disk host stars detected so far peak at wavelengths around ~100{\mu}m or shorter. However, six out of 31 excess sources in the Herschel OTKP DUNES have been seen to show significant - and in some cases extended - excess emission at 160{\mu}m, which is larger than the 100{\mu}m excess. This excess emission has been suggested to stem from debris disks colder than those known previously. Using several methods, we re-consider whether some or even all of the candidates may be associated with unrelated galactic or extragalactic emission and conclude that it is highly unlikely that none of the candidates represents a true circumstellar disk. For true disks, both the dust temperatures inferred from the SEDs and the disk radii estimated from the images suggest that the dust is nearly as cold as a blackbody. This requires the grains to be larger than ~100{\mu}m, regardless of their material composition. To explain the dearth of small grains, we explore several conceivable scenarios: transport-dominated disks, disks of low dynamical excitation, and disks of unstirred primordial macroscopic grains. Our qualitative analysis and collisional simulations rule out the first two of these scenarios, but show the feasibility of the third one. We show that such disks can survive for gigayears, largely preserving the primordial size distribution. They should be composed of macroscopic solids larger than millimeters, but smaller than kilometers in size. Thus planetesimal formation, at least in the outer regions of the systems, has stopped before "cometary" or "asteroidal" sizes were reached.

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.

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.

The search for planets around White Dwarf stars, and evidence for dynamical instability around them in the form of atmospheric pollution and circumstellar discs, raises questions about the nature of planetary systems that can survive the vicissitudes of the Asymptotic Giant Branch (AGB). We study the competing effects, on planets at several AU from the star, of strong tidal forces arising from the star's large convective envelope, and of the planets' orbital expansion due to stellar mass loss. We, for the first time, study the evolution of planets while following each thermal pulse on the AGB. For Jovian planets, tidal forces are strong, and can pull into the envelope planets initially at ~3 AU for a 1M_Sol star and ~5 AU for a 5M_Sol star. Lower-mass planets feel weaker tidal forces, and Terrestrial planets initially within 1.5-3 AU enter the stellar envelope. Thus, low-mass planets that begin inside the maximum stellar radius can survive, as their orbits expand due to mass loss. The inclusion of a moderate planetary eccentricity slightly strengthens the tidal forces experienced by Jovian planets. Eccentric Terrestrial planets are more at risk, since their eccentricity does not decay and their small pericentre takes them inside the stellar envelope. We also find the closest radii at which planets will be found around White Dwarfs, assuming that any planet entering the stellar envelope is destroyed. Planets are in that case unlikely to be found inside ~1.5 AU of a White Dwarf with a 1M_Sol progenitor and ~10 AU of a White Dwarf with a 5M_Sol progenitor.

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.

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.

A quarter of DA white dwarfs are metal polluted, yet elements heavier than helium sink down through the stellar atmosphere on timescales of days. Hence, these white dwarfs must be currently accreting material containing heavy elements. Here, we consider whether the scattering of comets or asteroids from an outer planetary system, following stellar mass loss on the asymptotic giant branch, can reproduce these observations. We use N-body simulations to investigate the effects of stellar mass loss on a simple system consisting of a planetesimal belt whose inner edge is truncated by a planet. Our simulations find that, starting with a planetesimal belt population fitted to the observed main sequence evolution, sufficient mass is scattered into the inner planetary system to explain the inferred heavy element accretion rates. This assumes that some fraction of the mass scattered into the inner planetary system ends up on star-grazing orbits, is tidally disrupted and accreted onto the white dwarf. The simulations also reproduce the observed decrease in accretion rate with cooling age and predict accretion rates in old (>1Gyr) white dwarfs, in line with observations. The efficiency we assumed for material scattered into the inner planetary system to end up on star-grazing orbits is based on a Solar-like planetary system, since the simulations show that a single planet is not sufficient. Although the correct level of accretion is reproduced, the simulations predict a higher fraction of accreting white dwarfs than observed. This could indicate that evolved planetary systems are less efficient at scattering bodies onto star-grazing orbits or that dynamical instabilities post-stellar mass loss cause rapid planetesimal belt depletion for a significant fraction of systems.

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]

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, although for certain migration rates, capture probability peaks at a finite eccentricity. 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. The model is applicable to many scenarios, including (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. Full details can be found in \cite{2010submitted}. (Abridged)

Detectable debris discs are thought to require dynamical excitation (`stirring'), so that planetesimal collisions release large quantities of dust. We investigate the effects of the secular perturbations of a planet, which may lie at a significant distance from the planetesimal disc, to see if these perturbations can stir the disc, and if so over what time-scale. The secular perturbations cause orbits at different semi-major axes to precess at different rates, and after some time t_cross initially non-intersecting orbits begin to cross. We show that t_cross is proportional to a_disc^(9/2)/(m_pl e_pl a_pl^3), where m_pl, e_pl and a_pl are the mass, eccentricity, and semi-major axis of the planet, and a_disc is the semi-major axis of the disc. This time-scale can be faster than that for the growth of planetesimals to Pluto's size within the outer disc. We also calculate the magnitude of the relative velocities induced amongst planetesimals and infer that a planet's perturbations can typically cause destructive collisions out to 100's of AU. Recently formed planets can thus have a significant impact on planet formation in the outer disc which may be curtailed by the formation of giant planets much closer to the star. The presence of an observed debris disc does not require the presence of Pluto-sized objects within it, since it can also have been stirred by a planet not in the disc. For the star epsilon Eridani, we find that the known RV planet can excite the planetesimal belt at 60 AU sufficiently to cause destructive collisions of bodies up to 100 km in size, on a time-scale of 40 Myr.