We simulate the passage through the Sun-Jupiter system of interstellar objects (ISOs) similar to 1I/`Oumuamua or 2I/Borisov. Capture of such objects is rare and overwhelmingly from low incoming speeds onto orbits akin to those of known long-period comets. This suggests that some of these comets could be of extra-solar origin, in particular inactive ones. Assuming ISOs follow the local stellar velocity distribution, we infer a volume capture rate of $0.051\,\mathrm{au}^3 \mathrm{yr}^{-1}$. Current estimates for orbital lifetimes and space densities then imply steady-state captured populations of $\sim10^2$ comets and $\sim10^5$ `Oumuamua-like rocks, of which 0.033% are within 6au at any time.

We perform $N$-body simulations of the early phases of open cluster evolution including a large population of planetesimals, initially arranged in Kuiper-belt like discs around each star. Using a new, 4th-order and time-reversible $N$-body code on Graphics Processing Units (GPUs), we evolve the whole system under the stellar gravity, i.e. treating planetesimals as test particles, and consider two types of initial cluster models, similar to IC348 and the Hyades, respectively. In both cases, planetesimals can be dynamically excited, transferred between stars or liberated to become free-floating (such as A/2017 U1 or 'Oumuamua) during the early cluster evolution. We find that planetesimals captured from another star are not necessarily dynamically distinct from those native to a star. After an encounter both native and captured planetesimals can exhibit aligned periastrons, qualitatively similar to that seen in the Solar system and commonly thought to be the signature of Planet 9. We discuss the implications of our results for both our Solar system and exoplanetary systems.

We hypothesise that differences in the temperatures at which the rocky material condensed out of the nebula gas can lead to differences in the composition of key rocky species (e.g., Fe, Mg, Si, Ca, Al, Na) and thus planet bulk density. Such differences in the observed bulk density of planets may occur as a function of radial location and time of planet formation. In this work we show that the predicted differences are on the cusp of being detectable with current instrumentation. In fact, for HD 219134, the 10 % lower bulk density of planet b compared to planet c could be explained by enhancements in Ca, Al rich minerals. However, we also show that the 11 % uncertainties on the individual bulk densities are not sufficiently accurate to exclude the absence of a density difference as well as differences in volatile layers. Besides HD 219134 b, we demonstrate that 55 Cnc e and WASP-47 e are similar candidates of a new Super-Earth class that have no core and are rich in Ca and Al minerals which are among the first solids that condense from a cooling proto-planetary disc. Planets of this class have densities 10-20% lower than Earth-like compositions and may have very different interior dynamics, outgassing histories and magnetic fields compared to the majority of Super-Earths.

The large-scale magnetic fields observed in spiral disc galaxies are often thought to result from dynamo action in the disc plane. However, the increasing importance of Faraday depolarization along any line of sight towards the galactic plane suggests that the strongest polarization signal may come from well above (~0.3-1 kpc) this plane, from the vicinity of the warm interstellar medium (WIM)/halo interface. We propose (see also Henriksen & Irwin 2016) that the observed spiral fields (polarization patterns) result from the action of vertical shear on an initially poloidal field. We show that this simple model accounts for the main observed properties of large-scale fields. We speculate as to how current models of optical spiral structure may generate the observed arm/interarm spiral polarization patterns.

We present two-dimensional hydrodynamical simulations of pairs of planets migrating simultaneously in the Type I regime in a protoplanetary disc. Convergent migration naturally leads to the trapping of these planets in mean-motion resonances. Once in resonance the planets' eccentricity grows rapidly, and disc-planet torques cause the planets to escape resonance on a time-scale of a few hundred orbits. The effect is more pronounced in highly viscous discs, but operates efficiently even in inviscid discs. We attribute this resonance-breaking to overstable librations driven by moderate eccentricity damping, but find that this mechanism operates differently in hydrodynamic simulations than in previous analytic calculations. Planets escaping resonance in this manner can potentially explain the observed paucity of resonances in Kepler multi-transiting systems, and we suggest that simultaneous disc-driven migration remains the most plausible means of assembling tightly-packed planetary systems.

The limited completeness of the Kepler sample for planets with orbital periods $\gtrsim$ 1 yr leaves open the possibility that exoplanetary systems may host undetected giant planets. Should such planets exist, their dynamical interactions with the inner planets may prove vital in sculpting the final orbital configurations of these systems. Using an $N$-body code with additional forces to emulate the effects of a protoplanetary disc, we perform simulations of the assembly of compact systems of super-Earth-mass planets with unseen giant companions. The simulated systems are analogous to Kepler-11 or Kepler-32 in that they contain 4 or 5 inner super-Earths, but our systems also contain longer-period giant companions which are unlikely to have been detected by Kepler. We find that giant companions tend to break widely-spaced, first-order mean-motion resonances, allowing the inner planets to migrate into tighter resonances. This leads to more compact architectures and increases the occurrence rate of Laplace resonant chains.

The Kepler mission has recently discovered a number of exoplanetary systems, such as Kepler-11 and Kepler-32, in which ensembles of several planets are found in very closely packed orbits (often within a few percent of an AU of one another). These compact configurations present a challenge for traditional planet formation and migration scenarios. We present a dynamical study of the assembly of these systems, using an N-body method which incorporates a parametrized model of planet migration in a turbulent protoplanetary disc. We explore a wide parameter space, and find that under suitable conditions it is possible to form compact, close-packed planetary systems via traditional disc-driven migration. We find that simultaneous migration of multiple planets is a viable mechanism for the assembly of tightly-packed planetary systems, as long as the disc provides significant eccentricity damping and the level of turbulence in the disc is modest. We discuss the implications of our preferred parameters for the protoplanetary discs in which these systems formed, and comment on the occurrence and significance of mean-motion resonances in our simulations.