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

The interiors of giant planets remain poorly understood. Even for the planets in the Solar System, difficulties in observation lead to major uncertainties in the properties of planetary cores. Exoplanets that have undergone rare evolutionary pathways provide a new route to understanding planetary interiors. We present the discovery of TOI-849b, the remnant core of a giant planet, with a radius smaller than Neptune but an anomalously high mass $M_p=40.8^{+2.4}_{-2.5}M_{\oplus}$ and density of $5.5\pm0.8$gcm$^{-3}$, similar to the Earth. Interior structure models suggest that any gaseous envelope of pure hydrogen and helium consists of no more than $3.9^{+0.8}_{-0.9}$% of the total mass of the planet. TOI-849b transits a late G type star (T$_{\rm mag}=11.5$) with an orbital period of 18.4 hours, leading to an equilibrium temperature of 1800K. The planet's mass is larger than the theoretical threshold mass for runaway gas accretion. As such, the planet could have been a gas giant before undergoing extreme mass loss via thermal self-disruption or giant planet collisions, or it avoided substantial gas accretion, perhaps through gap opening or late formation. Photoevaporation rates cannot provide the mass loss required to reduce a Jupiter-like gas giant, but can remove a few $M_\oplus$ hydrogen and helium envelope on timescales of several Gyr, implying that any remaining atmosphere is likely to be enriched by water or other volatiles from the planetary interior. TOI-849b represents a unique case where material from the primordial core is left over from formation and available to study.

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

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

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

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

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

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.

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

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

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

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

The Transiting Exoplanet Survey Satellite (TESS) is performing a near all-sky survey for planets that transit bright stars. In addition, its excellent photometric precision enables asteroseismology of solar-type and red-giant stars, which exhibit convection-driven, solar-like oscillations. Simulations predict that TESS will detect solar-like oscillations in nearly 100 stars already known to host planets. In this paper, we present an asteroseismic analysis of the known red-giant host stars HD 212771 and HD 203949, both systems having a long-period planet detected through radial velocities. These are the first detections of oscillations in previously known exoplanet-host stars by TESS, further showcasing the mission's potential to conduct asteroseismology of red-giant stars. We estimate the fundamental properties of both stars through a grid-based modeling approach that uses global asteroseismic parameters as input. We discuss the evolutionary state of HD 203949 in depth and note the large discrepancy between its asteroseismic mass ($M_\ast = 1.23 \pm 0.15\,{\rm M}_\odot$ if on the red-giant branch or $M_\ast = 1.00 \pm 0.16\,{\rm M}_\odot$ if in the clump) and the mass quoted in the discovery paper ($M_\ast = 2.1 \pm 0.1\,{\rm M}_\odot$), implying a change $>30\,\%$ in the planet's mass. Assuming HD 203949 to be in the clump, we investigate the planet's past orbital evolution and discuss how it could have avoided engulfment at the tip of the red-giant branch. Finally, HD 212771 was observed by K2 during its Campaign 3, thus allowing for a preliminary comparison of the asteroseismic performances of TESS and K2. We estimate the ratio of the observed oscillation amplitudes for this star to be $A_{\rm max}^{\rm TESS}/A_{\rm max}^{\rm K2} = 0.75 \pm 0.14$, consistent with the expected ratio of $\sim0.85$ due to the redder bandpass of TESS.

A gas giant planet which survives the giant branch stages of evolution at a distance of many au and then is subsequently perturbed sufficiently close to a white dwarf will experience orbital shrinkage and circularization due to star-planet tides. The circularization timescale, when combined with a known white dwarf cooling age, can place coupled constraints on the scattering epoch as well as the active tidal mechanisms. Here, we explore this coupling across the entire plausible parameter phase space by computing orbit shrinkage and potential self-disruption due to chaotic f-mode excitation and heating in planets on orbits with eccentricities near unity, followed by weakly dissipative equilibrium tides. We find that chaotic f-mode evolution activates only for orbital pericentres which are within twice the white dwarf Roche radius, and easily restructures or destroys ice giants but not gas giants. This type of internal thermal destruction provides an additional potential source of white dwarf metal pollution. Subsequent tidal evolution for the surviving planets is dominated by non-chaotic equilibrium and dynamical tides which may be well-constrained by observations of giant planets around white dwarfs at early cooling ages.

Planetary debris is observed in the atmospheres of over 1,000 white dwarfs, and two white dwarfs are now observed to contain orbiting minor planets. Exoasteroids and planetary core fragments achieve orbits close to the white dwarf through scattering with major planets. However, the architectures that allow for this scattering to take place are time-consuming to explore with $N$-body simulations lasting $\sim 10^{10}$ yr; these long-running simulations restrict the amount of phase space that can be investigated. Here we use planar and three-dimensional (spatial) elliptic periodic orbits, as well as chaotic indicators through dynamical stability maps, as quick scale-free analytic alternatives to $N$-body simulations in order to locate and predict instability in white dwarf planetary systems that consist of one major and one minor planet on very long timescales. We then classify the instability according to ejection versus collisional events. We generalized our previous work by allowing eccentricity and inclination of the periodic orbits to increase, thereby adding more realism but also significantly more degrees of freedom to our architectures. We also carried out a suite of computationally expensive 10 Gyr $N$-body simulations to provide comparisons with chaotic indicators in a limited region of phase space. We compute dynamical stability maps that are specific to white dwarf planetary systems and that can be used as tools in future studies to quickly estimate pollution prospects and timescales for one-planet architectures. We find that these maps also agree well with the outcomes of our $N$-body simulations. As observations of metal-polluted white dwarfs mount exponentially, particularly in the era of Gaia, tools such as periodic orbits can help infer dynamical histories for ensembles of systems.

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

The discovery of the intact minor planet embedded in the debris disc orbiting SDSS J1228+1040 raises questions about the dynamical history of the system. Further, the recent passage of the potentially interstellar object 1I/'Oumuamua within the solar system has re-ignited interest in minor body flux through exoplanetary systems. Here, we utilize the new analytical formalism from Grishin, Perets and Avni (2019) to estimate the rate at which the gaseous components of typical white dwarf discs trap an exo-planetesimal. We compare the types of captured orbits which arise from planetesimals originating from the interstellar medium, exo-Kuiper belts and exo-Oort clouds. We find that the rate of interstellar medium injection is negligible, whereas capture of both exo-Kuiper and exo-Oort cloud planetesimals is viable, but strongly size-dependent. For a gaseous disc which extends much beyond its Roche limit, capture is more probable than disruption at the Roche limit. We find that the capture probability linearly increases with the radial extent of the disc. Even in systems without minor planets, capture of smaller bodies will change the disc size distribution and potentially its temporal variability. Our formalism is general enough to be applied to future discoveries of embedded planetesimals in white dwarf debris discs.

Structure in the planet distribution provides an insight into the processes that shape the formation and evolution of planets. The Kepler mission has led to an abundance of statistical discoveries in regards to planetary radius, but the number of observed planets with measured masses is much smaller. By incorporating results from recent mass determination programs, we have discovered a new gap emerging in the planet population for sub-Neptune mass planets with orbital periods less than 20 days. The gap follows a slope of decreasing mass with increasing orbital period, has a width of a few $M_\oplus$, and is potentially completely devoid of planets. Fitting gaussian mixture models to the planet population in this region favours a bimodel distribution over a unimodel one with a reduction in Bayesian Information Criterion (BIC) of 19.9, highlighting the gap significance. We discuss several processes which could generate such a feature in the planet distribution, including a pileup of planets above the gap region, tidal interactions with the host star, dynamical interactions with the disk, with other planets, or with accreting material during the formation process.

The discovery of the intact metallic planetary core fragment orbiting the white dwarf SDSS J1228+1040 within one Solar radius highlights the possibility of detecting larger, unfragmented conducting cores around magnetic white dwarfs through radio emission. Previous models of this decades-old idea focussed on determining survivability of the cores based on their inward Lorentz drift towards the star. However, gravitational tides may represent an equal or dominant force. Here, we couple both effects by assuming a Maxwell rheological model and performing simulations over the entire range of observable white dwarf magnetic field strengths (10^3 -- 10^9 G) and their potential atmospheric electrical conductivities (10^{-1} -- 10^4 S/m) in order to more accurately constrain survivability lifetimes. This force coupling allows us to better pinpoint the physical and orbital parameters which allow planetary cores to survive for over a Gyr, maximizing the possibility that they can be detected. The most robust survivors showcase high dynamic viscosities (>~ 10^{24} Pa*s) and orbit within kG-level magnetic fields.

Practically all known planet hosts will evolve into white dwarfs, and large parts of their planetary systems will survive this transition - the same is true for the solar system beyond the orbit of Mars. Spectroscopy of white dwarfs accreting planetary debris provides the most accurate insight into the bulk composition of exo-planets. Ground-based spectroscopic surveys of ~260, 000 white dwarfs detected with Gaia will identify >1000 evolved planetary systems, and high-throughput high-resolution space-based ultraviolet spectroscopy is essential to measure in detail their abundances. So far, evidence for two planetesimals orbiting closely around white dwarfs has been obtained, and their study provides important constraints on the composition and internal structure of these bodies. Major photometric and spectroscopic efforts will be necessary to assemble a sample of such close-in planetesimals that is sufficiently large to establish their properties as a population, and to deduce the architectures of the outer planetary systems from where they originated. Mid-infrared spectroscopy of the dusty disks will provide detailed mineralogical information of the debris, which, in combination with the elemental abundances measured from the white dwarf spectroscopy, will enable detailed physical modelling of the chemical, thermodynamic, and physical history of the accreted material. Flexible multi-epoch infrared observations are essential to determine the physical nature, and origin of the variability observed in many of the dusty disks. Finally, the direct detection of the outer reservoirs feeding material to the white dwarfs will require sensitive mid- and far-infrared capabilities.

Observational evidence of white dwarf planetary systems is dominated by the remains of exo-asteroids through accreted metals, debris discs, and orbiting planetesimals. However, exo-planets in these systems play crucial roles as perturbing agents, and can themselves be perturbed close to the white dwarf Roche radius. Here, we illustrate a procedure for computing the tidal interaction between a white dwarf and a near-spherical solid planet. This method determines the planet's inward and/or outward drift, and whether the planet will reach the Roche radius and be destroyed. We avoid constant tidal lag formulations and instead employ the self-consistent secular Darwin-Kaula expansions from Bou\'{e} & Efroimsky (2019), which feature an arbitrary frequency dependence on the quality functions. We adopt wide ranges of dynamic viscosities and spin rates for the planet in order to straddle many possible outcomes, and provide a foundation for the future study of individual systems with known or assumed rheologies. We find that: (i) massive Super-Earths are destroyed more readily than minor planets (such as the ones orbiting WD 1145+017 and SDSS J1228+1040), (ii) low-viscosity planets are destroyed more easily than high-viscosity planets, and (iii) the boundary between survival and destruction is likely to be fractal and chaotic.

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.

Understanding the fate of planetary systems through white dwarfs which accrete debris crucially relies on tracing the orbital and physical properties of exo-asteroids during the giant branch phase of stellar evolution. Giant branch luminosities exceed the Sun's by over three orders of magnitude, leading to significantly enhanced Yarkovsky and YORP effects on minor planets. Here, we place bounds on Yarkovsky-induced differential migration between asteroids and planets during giant branch mass loss by modelling one exo-Neptune with inner and outer exo-Kuiper belts. In our bounding models, the asteroids move too quickly past the planet to be diverted from their eventual fate, which can range from: (i) populating the outer regions of systems out to 10^4-10^5 au, (ii) being engulfed within the host star, or (iii) experiencing Yarkovsky-induced orbital inclination flipping without any Yarkovsky-induced semimajor axis drift. In these violent limiting cases, temporary resonant trapping of asteroids with radii of under about 10 km by the planet is insignificant, and capture within the planet's Hill sphere requires fine-tuned dissipation. The wide variety of outcomes presented here demonstrates the need to employ sophisticated structure and radiative exo-asteroid models in future studies. Determining where metal-polluting asteroids reside around a white dwarf depends on understanding extreme Yarkovsky physics.

Solar system, exoplanet and stellar science rely on transits, eclipses and occultations for dynamical and physical insight. Often, the geometry of these configurations are modelled by assuming a particular viewpoint. Here, instead, I derive user-friendly formulae from first principles independent of viewpoint and in three dimensions. I generalise the results of Veras & Breedt (2017) by (i) characterising three-body systems which are in transit but are not necessarily perfectly aligned, and by (ii) incorporating motion. For a given snapshot in time, I derive explicit criteria to determine whether a system is in or out of transit, if an eclipse is total or annular, and expressions for the size of the shadow, including their extreme values and a condition for engulfment. These results are exact. For orbital motion, I instead obtain approximate results. By assuming fixed orbits, I derive a single implicit algebraic relation which can be solved to obtain the frequency and duration of transit events -- including ingresses and egresses -- for combinations of moons, planets and stars on arbitrarily inclined circular orbits; the eccentric case requires the solution of Kepler's equation but remains algebraic. I prove that a transit shadow -- whether umbral, antumbral or penumbral -- takes the shape of a parabolic cylinder, and finally present geometric constraints on Earth-based observers hoping to detect a three-body syzygy (or perfect alignment) -- either in extrasolar systems or within the solar system -- potentially as a double annular eclipse.

As discoveries of multiple planets in the habitable zone of their parent star mount, developing analytical techniques to quantify extrasolar intra-system panspermia will become increasingly important. Here, we provide user-friendly prescriptions that describe the asteroid impact characteristics which would be necessary to transport life both inwards and outwards within these systems within a single framework. Our focus is on projectile generation and delivery and our expressions are algebraic, eliminating the need for the solution of differential equations. We derive a probability distribution function for life-bearing debris to reach a planetary orbit, and describe the survival of micro-organisms during planetary ejection, their journey through interplanetary space, and atmospheric entry.

Previous full-lifetime simulations of single-star multi-planet systems across all phases of stellar evolution have predominately assumed coplanar or nearly-coplanar orbits. Here we assess the consequences of this assumption by removing it and exploring the effect of giant branch mass loss on the stability of two-planet systems with small to moderate non-Kozai (<40 degrees) relative inclinations. We run nearly 10^4 simulations over 14 Gyr for F-star, A-star and B-star planet hosts, incorporating main sequence stellar masses of 1.5, 2.0, 2.5, 3.0 and 5.0 solar masses, and initial planetary semimajor axis ratios that straddle their three-dimensional Hill stability limits. We find that the near-coplanar assumption can approximate well the stability frequencies and critical separations found for higher inclinations, except around strong mean-motion commensurabilities. Late instabilities -- after the star has become a white dwarf -- occur throughout the explored mutual inclination range. Consequently, non-Kozai mutual inclination should not be used as a predictive orbital proxy for determining which white dwarf multi-planet systems discovered by Gaia should represent high-priority follow-up targets for the detection of metal pollution and planetary debris discs.

The most heavily polluted white dwarfs often show excess infrared radiation from circumstellar dust disks, which are modeled as a result of tidal disruption of extrasolar minor planets. Interaction of dust, gas, and disintegrating objects can all contribute to the dynamical evolution of these dust disks. Here, we report on two infrared variable dusty white dwarfs, SDSS J1228+1040 and G29-38. For SDSS J1228+1040, compared to the first measurements in 2007, the IRAC [3.6] and [4.5] fluxes decreased by 20% by 2014 to a level also seen in the recent 2018 observations. For G29-38, the infrared flux of the 10 $\mu$m silicate emission feature became 10% stronger between 2004 and 2007, We explore several scenarios that could account for these changes, including tidal disruption events, perturbation from a companion, and runaway accretion. No satisfactory causes are found for the flux drop in SDSS J1228+1040 due to the limited time coverage. Continuous tidal disruption of small planetesimals could increase the mass of small grains and concurrently change the strength of the 10 $\mu$m feature of G29-38. Dust disks around white dwarfs are actively evolving and we speculate that there could be different mechanisms responsible for the temporal changes of these disks.

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.

'Oumuamua was discovered passing through our Solar System on a hyperbolic orbit. It presents an apparent contradiction, with colors similar to those of volatile-rich Solar System bodies but with no visible outgassing or activity during its close approach to the Sun. Here we show that this contradiction can be explained by the dynamics of planetesimal ejection by giant planets. We propose that 'Oumuamua is an extinct fragment of a comet-like planetesimal born in a planet-forming disk that also formed Neptune- to Jupiter-mass giant planets. On its pathway to ejection 'Oumuamua's parent body underwent a close encounter with a giant planet and was tidally disrupted into small pieces, similar to comet Shoemaker-Levy 9's disruption after passing close to Jupiter. We use dynamical simulations to show that 0.1-1% of cometary planetesimals undergo disruptive encounters prior to ejection. Rocky asteroidal planetesimals are unlikely to disrupt due to their higher densities. After disruption, the bulk of fragments undergo enough close passages to their host stars to lose their surface volatiles and become extinct. Planetesimal fragments such as 'Oumuamua contain little of the mass in the population of interstellar objects but dominate by number. Our model makes predictions that will be tested in the coming decade by LSST.

'Oumuamua, the first bona-fide interstellar planetesimal, was discovered passing through our Solar System on a hyperbolic orbit. This object was likely dynamically ejected from an extrasolar planetary system after a series of close encounters with gas giant planets. To account for 'Oumuamua's detection, simple arguments suggest that ~1 Earth mass of planetesimals are ejected per Solar mass of Galactic stars. However, that value assumes mono-sized planetesimals. If the planetesimal mass distribution is instead top-heavy the inferred mass in interstellar planetesimals increases to an implausibly high value. The tension between theoretical expectations for the planetesimal mass function and the observation of 'Oumuamua can be relieved if a small fraction (~0.1-1%) of planetesimals are tidally disrupted on the pathway to ejection into 'Oumuamua-sized fragments. Using a large suite of simulations of giant planet dynamics including planetesimals, we confirm that 0.1-1% of planetesimals pass within the tidal disruption radius of a gas giant on their pathway to ejection. 'Oumuamua may thus represent a surviving fragment of a disrupted planetesimal. Finally, we argue that an asteroidal composition is dynamically disfavoured for 'Oumuamua, as asteroidal planetesimals are both less abundant and ejected at a lower efficiency than cometary planetesimals.

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.

The atmospheres of between one quarter and one half of observed single white dwarfs in the Milky Way contain heavy element pollution from planetary debris. The pollution observed in white dwarfs in binary star systems is, however, less clear, because companion star winds can generate a stream of matter which is accreted by the white dwarf. Here we (i) discuss the necessity or lack thereof of a major planet in order to pollute a white dwarf with orbiting minor planets in both single and binary systems, and (ii) determine the critical binary separation beyond which the accretion source is from a planetary system. We hence obtain user-friendly functions relating this distance to the masses and radii of both stars, the companion wind, and the accretion rate onto the white dwarf, for a wide variety of published accretion prescriptions. We find that for the majority of white dwarfs in known binaries, if pollution is detected, then that pollution should originate from planetary material.

Doppler-based planet surveys point to an increasing occurrence rate of giant planets with stellar mass. Such surveys rely on evolved stars for a sample of intermediate-mass stars (so-called retired A stars), which are more amenable to Doppler observations than their main-sequence progenitors. However, it has been hypothesised that the masses of subgiant and low-luminosity red-giant stars targeted by these surveys --- typically derived from a combination of spectroscopy and isochrone fitting --- may be systematically overestimated. Here, we test this hypothesis for the particular case of the exoplanet-host star HD 212771 using K2 asteroseismology. The benchmark asteroseismic mass ($1.45^{+0.10}_{-0.09}\:\text{M}_{\odot}$) is significantly higher than the value reported in the discovery paper ($1.15\pm0.08\:\text{M}_{\odot}$), which has been used to inform the stellar mass-planet occurrence relation. This result, therefore, does not lend support to the above hypothesis. Implications for the fates of planetary systems are sensitively dependent on stellar mass. Based on the derived asteroseismic mass, we predict the post-main-sequence evolution of the Jovian planet orbiting HD 212771 under the effects of tidal forces and stellar mass loss.

Although conjunctions and oppositions frequently occur in planetary systems, eclipse-related phenomena are usually described from an Earth-centric perspective. Space missions to different parts of the Solar system, as well as the mounting number of known exo-planets in habitable zones and the possibility of sending featherweight robot spacecraft to them, prompt broader considerations. Here, we derive the geometry of eclipses, transits and occultations from a primarily exo-Earth viewpoint, and apply the formulation to the Solar system and three types of three-body extrasolar planetary systems: with 1 star and 2 planets (Case I), with 2 stars and 1 planet (Case II), and with 1 planet, 1 star and 1 moon (Case III). We derive the general conditions for total, partial and annular eclipses to occur at exo-syzygy, and implement them in each case in concert with stability criteria. We then apply the formalism to the TRAPPIST-1, Kepler-444 and Kepler-77 systems -- the first of which contains multiple potentially habitable planets -- and provide reference tables of both Solar system and TRAPPIST-1 syzygy properties. We conclude by detailing a basic algebraic algorithm which can be used to quickly characterize eclipse properties in any three-body system.

High-metallicity pollution is common in white dwarf (WD) stars hosting remnant planetary systems. However, they rarely have detectable debris accretion discs, possibly because much of the influx is fast steeply-infalling debris in star-grazing orbits, producing a more tenuous signature than a slowly accreting disk. Processes governing such deposition between the Roche radius and photosphere have so far received little attention and we model them here analytically by extending recent work on sun-grazing comets to WD systems. We find that the evolution of cm-to-km size (a_0) infallers most strongly depends on two combinations of parameters, which effectively measure sublimation rate and binding strength. We then provide an algorithm to determine the fate of infallers for any WD, and apply the algorithm to four limiting combinations of hot versus cool (young/old) WDs with snowy (weak, volatile) versus rocky (strong, refractory) infallers. We find: (i) Total sublimation above the photosphere befalls all small infallers across the entire WD temperature (T_WD) range, the threshold size rising with T_WD and 100X larger for rock than snow. (ii) All very large objects fragment tidally regardless of T_WD: for rock, a_0 >= 10^5 cm; for snow, a_0 >= 10^3 -- 3x10^4 cm across all WD cooling ages. (iii) A considerable range of a_0 avoids fragmentation and total sublimation, yielding impacts or grazes with cold WDs. This range narrows rapidly with increasing T_WD, especially for snowy bodies. Finally, we discuss briefly how the various forms of deposited debris may finally reach the photosphere surface itself.

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.

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

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

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

Mounting discoveries of debris discs orbiting newly-formed stars and white dwarfs (WDs) showcase the importance of modeling the long-term evolution of small bodies in exosystems. WD debris discs are in particular thought to form from very long-term (0.1-5.0 Gyr) instability between planets and asteroids. However, the time-consuming nature of $N$-body integrators which accurately simulate motion over Gyrs necessitates a judicious choice of initial conditions. The analytical tools known as \textit{periodic orbits} can circumvent the guesswork. Here, we begin a comprehensive analysis directly linking periodic orbits with $N$-body integration outcomes with an extensive exploration of the planar circular restricted three-body problem (CRTBP) with an outer planet and inner asteroid near or inside of the $2$:$1$ mean motion resonance. We run nearly 1000 focused simulations for the entire age of the Universe (14 Gyr) with initial conditions mapped to the phase space locations surrounding the unstable and stable periodic orbits for that commensurability. In none of our simulations did the planar CRTBP architecture yield a long-timescale ($\gtrsim 0.25$% of the age of the Universe) asteroid-star collision. The pericentre distance of asteroids which survived beyond this timescale ($\approx 35$ Myr) varied by at most about 60%. These results help affirm that collisions occur too quickly to explain WD pollution in the planar CRTBP $2$:$1$ regime, and highlight the need for further periodic orbit studies with the eccentric and inclined TBP architectures and other significant orbital period commensurabilities.

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

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

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

We report the discovery of two transiting planetary systems: a super dense, sub-Jupiter mass planet WASP-86b (\mpl\ = 0.82 $\pm$ 0.06 \mj, \rpl\ = 0.63 $\pm$ 0.01 \rj), and a bloated, Saturn-like planet WASP-102b (\mpl\ = 0.62 $\pm$ 0.04 \mj, \rpl\=1.27 $\pm$ 0.03 \rj). They orbit their host star every $\sim$5.03, and $\sim$2.71 days, respectively. The planet hosting WASP-86 is a F7 star (\teff\ = 6330$\pm$110 K, \feh\ = $+$0.23 $\pm$ 0.14 dex, and age $\sim$0.8--1~Gyr), WASP-102 is a G0 star (\teff\ = 5940$\pm$140 K, \feh\ = $-$0.09$\pm$ 0.19 dex, and age $\sim$1~Gyr). These two systems highlight the diversity of planetary radii over similar masses for giant planets with masses between Saturn and Jupiter. WASP-102b shows a larger than model-predicted radius, indicating that the planet is receiving a strong incident flux which contributes to the inflation of its radius. On the other hand, with a density of $\rho_{pl}$ = 3.24$\pm$~0.3~$\rho_{jup}$, WASP-86b is the densest gas giant planet among planets with masses in the range 0.05 $<M$_{pl}$<$ 2.0 \mj. With a stellar mass of 1.34 M$_{\odot}$ and \feh = $+$0.23 dex, WASP-86 could host additional massive and dense planets given that its protoplanetary disc is expected to also have been enriched with heavy elements. In order to match WASP-86b's density, an extrapolation of theoretical models predicts a planet composition of more than 80\% in heavy elements (whether confined in a core or mixed in the envelope). This fraction corresponds to a core mass of approximately 210\me\ for WASP-86b's mass of \mpl$\sim$260\,\me. Only planets with masses larger than about 2\mj\ have larger densities than that of WASP-86b, making it exceptional in its mass range.

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

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

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

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

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

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

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

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

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

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.

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

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

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

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

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

The solar system started to form about 4.56 Gyr ago and despite the long intervening time span, there still exist several clues about its formation. The three major sources for this information are meteorites, the present solar system structure and the planet-forming systems around young stars. In this introduction we give an overview of the current understanding of the solar system formation from all these different research fields. This includes the question of the lifetime of the solar protoplanetary disc, the different stages of planet formation, their duration, and their relative importance. We consider whether meteorite evidence and observations of protoplanetary discs point in the same direction. This will tell us whether our solar system had a typical formation history or an exceptional one. There are also many indications that the solar system formed as part of a star cluster. Here we examine the types of cluster the Sun could have formed in, especially whether its stellar density was at any stage high enough to influence the properties of today's solar system. The likelihood of identifying siblings of the Sun is discussed. Finally, the possible dynamical evolution of the solar system since its formation and its future are considered.

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

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

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

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

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

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

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

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

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

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.

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

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.

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

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.

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

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

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

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

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

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

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

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

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

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

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

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.

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

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

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

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

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

We present new precision radial velocities and a three-planet Keplerian orbit fit for the V = 8.5, G5 V star HIP 14810. We began observing this star at Keck Observatory as part of the N2K Planet Search Project. Wright et al. (2007) announced the inner two planets to this system, and subsequent observations have revealed the outer planet planet and the proper orbital solution for the middle planet. The planets have minimum masses of 3.9, 1.3, and 0.6 M_Jup and orbital periods of 6.67, 147.7, and 952 d, respectively. We have numerically integrated the family of orbital solutions consistent with the data and find that they are stable for at least 10^6 yr. Our photometric search shows that the inner planet does not transit.

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

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

The eccentricity evolution of multiple planet systems can provide valuable constraints on planet formation models. Unfortunately, the inevitable uncertainties in the current orbital elements can lead to significant ambiguities in the nature of the secular evolution. Integrating any single set of orbital elements inadequately describes the full range of secular evolutions consistent with current observations. Thus, we combine radial velocity observations of HD 12661 with Markov Chain Monte Carlo sampling to generate ensembles of initial conditions for direct n-body integrations. We find that any mean motion resonances are quite weak and do not significantly impact the secular evolution, and that current observations indicate circulation or large amplitude libration of the periapses. The eccentricity of the outer planet undergoes large oscillations for nearly all of the allowed two-planet orbital solutions. This type of secular evolution would arise if planet c had been impulsively perturbed, perhaps due to strong scattering of an additional planet that was subsequently accreted onto the star. Finally, we note that the secular evolution implied by the current orbital configuration implies that planet c spends ~96% of the time following an orbit more eccentric than that presently observed. Either this system is being observed during a relatively rare state, or additional planets are affecting the observed radial velocities and/or the system's secular eccentricity evolution.

The discovery of over 200 extrasolar planets with the radial velocity (RV) technique has revealed that many giant planets have large eccentricities, in striking contrast with most of the planets in the solar system and prior theories of planet formation. The realization that many giant planets have large eccentricities raises a fundamental question: ``Do terrestrial-size planets of other stars typically have significantly eccentric orbits or nearly circular orbits like the Earth?'' Here, we demonstrate that photometric observations of transiting planets could be used to characterize the orbital eccentricities for individual transiting planets, as well the eccentricity distribution for various populations of transiting planets (e.g., those with a certain range of orbital periods or physical sizes). Such characterizations can provide valuable constraints on theories for the excitation of eccentricities and tidal dissipation. We outline the future prospects of the technique given the exciting prospects for future transit searches, such as those to be carried out by the CoRoT and Kepler missions.

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

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

The large eccentricities of many giant extrasolar planets may represent the endpoint of gravitational scattering in initially more crowded systems. If so, the early evolution of the giant planets is likely to be more restrictive of terrestrial planet formation than would be inferred from the current, dynamically quiescent, configurations. Here, we study statistically the extent of the anti-correlation between giant planets and terrestrial planets expected in a scattering model. We use marginally stable systems of three giant planets, with a realistic range of planetary masses, as a simple model for the initial conditions prior to scattering, and show that after scattering the surviving planets reproduce well the known extrasolar planet eccentricities beyond a > 0.5 AU. By tracking the minimum periastron values of all planets during the evolution, we derive the distribution of orbital radii across which strong perturbations (from crossing orbits) are likely to affect low mass planet formation. We find that scattering affects inner planet formation at orbital separations less than 50% of the final periastron distance of the innermost massive planet in approximately 30% of the realizations, and can occasionally influence planet formation at orbital separations less than 20% of the final periastron distance of the innermost massive planet. The domain of influence of the scattering massive planets increases as the mass differential between the massive planets decreases. Observational study of the correlation between massive and terrestrial extrasolar planets in the same system has the potential to constrain the origin of planetary eccentricity.

In most extrasolar planetary systems, the present orbits of known giant planets admit the existence of stable terrestrial planets. Those same giant planets, however, have typically eccentric orbits that hint at violent early dynamics less benign for low mass planet formation. Under the assumption that massive planet eccentricities are the end point of gravitational scattering in multiple planet systems, we study the evolution of the building blocks of terrestrial planets during the scattering process. We find that typically, evolutionary sequences that result in a moderately eccentric giant planet orbiting at a ~ 2.5 AU eject over 95% of the material initially present within the habitable zone. Crossing orbits largely trigger the ejection, and leave the surviving material with a wide dispersion in semi-major axis, eccentricity and inclination. Based on these results, we predict that radial velocity follow-up of terrestrial planet systems found by Kepler will find that these are anti-correlated with the presence of eccentric giant planets orbiting at a few AU.

The significant orbital eccentricities of most giant extrasolar planets may have their origin in the gravitational dynamics of initially unstable multiple planet systems. In this work, we explore the dynamics of two close planets on inclined orbits through both analytical techniques and extensive numerical scattering experiments. We derive a criterion for two equal mass planets on circular inclined orbits to achieve Hill stability, and conclude that significant radial migration and eccentricity pumping of both planets occurs predominantly by 2:1 and 5:3 mean motion resonant interactions. Using Laplace-Lagrange secular theory, we obtain analytical secular solutions for the orbital inclinations and longitudes of ascending nodes, and use those solutions to distinguish between the secular and resonant dynamics which arise in numerical simulations. We also illustrate how encounter maps, typically used to trace the motion of massless particles, may be modified to reproduce the gross instability seen by the numerical integrations. Such a correlation suggests promising future use of such maps to model the dynamics of more coplanar massive planet systems.

We have undertaken a thorough dynamical investigation of five extrasolar planetary systems using extensive numerical experiments. The systems Gl 777 A, HD 72659, Gl 614, 47 Uma and HD 4208 were examined concerning the question of whether they could host terrestrial like planets in their habitable zones (=HZ). First we investigated the mean motion resonances between fictitious terrestrial planets and the existing gas giants in these five extrasolar systems. Then a fine grid of initial conditions for a potential terrestrial planet within the HZ was chosen for each system, from which the stability of orbits was then assessed by direct integrations over a time interval of 1 million years. The computations were carried out using a Lie-series integration method with an adaptive step size control. This integration method achieves machine precision accuracy in a highly efficient and robust way, requiring no special adjustments when the orbits have large eccentricities. The stability of orbits was examined with a determination of the Renyi entropy, estimated from recurrence plots, and with a more straight forward method based on the maximum eccentricity achieved by the planet over the 1 million year integration. Additionally, the eccentricity is an indication of the habitability of a terrestrial planet in the HZ; any value of e>0.2 produces a significant temperature difference on a planet's surface between apoapse and periapse. The results for possible stable orbits for terrestrial planets in habitable zones for the five systems are summarized as follows: for Gl 777 A nearly the entire HZ is stable, for 47 Uma, HD 72659 and HD 4208 terrestrial planets can survive for a sufficiently long time, while for Gl 614 our results exclude terrestrial planets moving in stable orbits within the HZ.

Observations of structure in circumstellar debris discs provide circumstantial evidence for the presence of massive planets at large (several tens of au) orbital radii, where the timescale for planet formation via core accretion is prohibitively long. Here, we investigate whether a population of distant planets can be produced via outward migration subsequent to formation in the inner disc. Two possibilities for significant outward migration are identified. First, cores that form early at radii of around 10 au can be carried to larger radii via gravitational interaction with the gaseous disc. This process is efficient if there is strong mass loss from the disc - either within a cluster or due to photoevaporation from a star more massive than the Sun - but does not require the extremely destructive environment found, for example, in the core of the Orion Nebula. We find that, depending upon the disc model, gas disc migration can yield massive planets (several Jupiter masses) at radii of around 20-50 au. Second, interactions within multiple planet systems can drive the outer planet into a large, normally highly eccentric orbit. A series of scattering experiments suggests that this process is most efficient for lower mass planets within systems of unequal mass ratio. This mechanism is a good candidate for explaining the origin of relatively low mass giant planets in eccentric orbits at large radii.