Following its flyby and first imaging the Pluto-Charon binary, the New Horizons spacecraft visited the Kuiper-Belt-Object (KBO) (486958) 2014 MU69 (Arrokoth). Imaging showed MU69 to be a contact-binary, made of two individual lobes connected by a narrow neck, rotating at low spin period (15.92 h), and having high obliquity (~98 deg), similar to other KBO contact-binaries inferred through photometric observations. The origin of such peculiar configurations is puzzling, and all scenarios suggested for the origins of contact-binaries fail to reproduce such properties and their likely high frequency. Here we show that semi-secular perturbations operating only on ultra-wide (~0.1-0.4 Hill-radius) KBO-binaries can robustly lead to gentle, slow-speed binary mergers at arbitrarily high obliquities, but low rotational velocities, that can reproduce MU69's (and similar oblique contact binaries) characteristics. Using N-body simulations, we find that ~15% of all ultra-wide binaries with cosine-uniform inclination distribution are likely to merge through this process. Moreover, we find that such mergers are sufficiently gentle as to only slightly deform the KBO shape, and can produce the measured rotation speed of MU69. The semi-secular contact-binary formation channel not only explains the observed properties of MU69, but could also apply for other Kuiper/asteroid belt binaries, and for Solar/extra-solar moon systems.

Exomoons orbiting terrestrial or super-terrestrial exoplanets have not yet been discovered; their possible existence and properties are therefore still an unresolved question. Here we explore the collisional formation of exomoons through giant planetary impacts. We make use of smooth particle hydrodynamical (SPH) collision simulations and survey a large phase-space of terrestrial/super-terrestrial planetary collisions. We characterize the properties of such collisions, finding one rare case in which an exomoon forms through a graze&capture scenario, in addition to a few graze&merge or hit&run scenarios. Typically however, our collisions form massive circumplanetary discs, for which we use follow-up N-body simulations in order to derive lower-limit mass estimates for the ensuing exomoons. We investigate the mass, long-term tidal-stability, composition and origin of material in both the discs and the exomoons. Our giant-impact models often generate relatively iron-rich moons, that form beyond the synchronous radius of the planet, and would thus tidally evolve outward with stable orbits, rather than be destroyed. Our results suggest that it is extremely difficult to collisionally form currently-detectable exomoons orbiting super-terrestrial planets, through single giant impacts. It might be possible to form massive, detectable exomoons through several mergers of smaller exomoons, formed by multiple impacts, however more studies are required in order to reach a conclusion. Given the current observational initiatives, the search should focus primarily on more massive planet categories. However, about a quarter of the exomoons predicted by our models are approximately Mercury-mass or more, and are much more likely to be detectable given a factor 2 improvement in the detection capability of future instruments, providing further motivation for their development.

It has been suggested that the comet-like activity of Main Belt Comets is due to the sublimation of sub-surface water-ice that is exposed when these objects are impacted by meter-sized bodies. We recently examined this scenario and showed that such impacts can in fact excavate ice and present a plausible mechanism for triggering the activation of MBCs (Haghighipour et al. 2016). However, because the purpose of that study was to prove the concept and identify the most viable ice-longevity model, the porosity of the object and the loss of ice due to the heat of impact were ignored. In this paper, we extend our impact simulations to porous materials and account for the loss of ice due to an impact. We show that for a porous MBC, impact craters are deeper, reaching to approximately 15 m implying that if the activation of MBCs is due to the sublimation of sub-surface ice, this ice has to be within the top 15 m of the object. Results also indicate that the loss of ice due to the heat of impact is negligible, and the re-accretion of ejected ice is small. The latter suggests that the activities of current MBCs are most probably from multiple impact sites. Our study also indicates that in order for sublimation from multiple sites to account for the observed activity of the currently known MBCs, the water content of MBCs (and their parent asteroids) needs to be larger than the values traditionally considered in models of terrestrial planet formation.

Future large-scale surveys with high resolution imaging will provide us with a few $10^5$ new strong galaxy-scale lenses. These strong lensing systems however will be contained in large data amounts which are beyond the capacity of human experts to visually classify in a unbiased way. We present a new strong gravitational lens finder based on convolutional neural networks (CNNs). The method was applied to the Strong Lensing challenge organised by the Bologna Lens Factory. It achieved first and third place respectively on the space-based data-set and the ground-based data-set. The goal was to find a fully automated lens finder for ground-based and space-based surveys which minimizes human inspect. We compare the results of our CNN architecture and three new variations ("invariant" "views" and "residual") on the simulated data of the challenge. Each method has been trained separately 5 times on 17 000 simulated images, cross-validated using 3 000 images and then applied to a 100 000 image test set. We used two different metrics for evaluation, the area under the receiver operating characteristic curve (AUC) score and the recall with no false positive ($\mathrm{Recall}_{\mathrm{0FP}}$). For ground based data our best method achieved an AUC score of $0.977$ and a $\mathrm{Recall}_{\mathrm{0FP}}$ of $0.50$. For space-based data our best method achieved an AUC score of $0.940$ and a $\mathrm{Recall}_{\mathrm{0FP}}$ of $0.32$. On space-based data adding dihedral invariance to the CNN architecture diminished the overall score but achieved a higher no contamination recall. We found that using committees of 5 CNNs produce the best recall at zero contamination and consistenly score better AUC than a single CNN. We found that for every variation of our CNN lensfinder, we achieve AUC scores close to $1$ within $6\%$.

It has been suggested that the comet-like activity of Main Belt Comets are due to the sublimation of sub-surface water-ice that has been exposed as a result of their surfaces being impacted by m-sized bodies. We have examined the viability of this scenario by simulating impacts between m-sized and km-sized objects using a smooth particle hydrodynamics approach. Simulations have been carried out for different values of the impact velocity and impact angle as well as different target material and water-mass fraction. Results indicate that for the range of impact velocities corresponding to those in the asteroid belt, the depth of an impact crater is slightly larger than 10 m suggesting that if the activation of MBCs is due to the sublimation of sub-surface water-ice, this ice has to exist no deeper than a few meters from the surface. Results also show that ice-exposure occurs in the bottom and on the interior surface of impact craters as well as the surface of the target where some of the ejected icy inclusions are re-accreted. While our results demonstrate that the impact scenario is indeed a viable mechanism to expose ice and trigger the activity of MBCs, they also indicate that the activity of the current MBCs is likely due to ice sublimation from multiple impact sites and/or the water contents of these objects (and other asteroids in the outer asteroid belt) is larger than the 5% that is traditionally considered in models of terrestrial planet formation providing more ice for sublimation. We present details of our simulations and discuss their results and implications.

We carry out three-dimensional Smoothed Particle Hydrodynamics simulations of spherical homogeneous SiO2 dust aggregates to investigate how the mass and the porosity of the aggregates affects their ability to survive an impact at various different collision velocities (between 1 - 27.5m/s). We explore how the threshold velocities for fragmentation vary with these parameters. Crucially, we find that the porosity plays a part of utmost importance in determining the outcome of collisions. In particular, we find that aggregates with filling factors >37% are significantly weakened and that the velocity regime in which the aggregates grow is reduced or even non-existent (instead, the aggregates either rebound off each other or break apart). At filling factors less than ~37% we find that more porous objects are weaker but not as weak as highly compact objects with filling factors >37%. In addition we find that (for a given aggregate density) collisions between very different mass objects have higher threshold velocities than those between very similar mass objects. We find that fragmentation velocities are higher than the typical values of 1m/s and that growth can even occur for velocities as high as 27.5m/s. Therefore, while the growth of aggregates is more likely if collisions between different sized objects occurs or if the aggregates are porous with filling factor <37%, it may also be hindered if the aggregates become too compact.

We perform numerical simulations of a disc-planet system using various grid-based and smoothed particle hydrodynamics (SPH) codes. The tests are run for a simple setup where Jupiter and Neptune mass planets on a circular orbit open a gap in a protoplanetary disc during a few hundred orbital periods. We compare the surface density contours, potential vorticity and smoothed radial profiles at several times. The disc mass and gravitational torque time evolution are analyzed with high temporal resolution. There is overall consistency between the codes. The density profiles agree within about 5% for the Eulerian simulations while the SPH results predict the correct shape of the gap although have less resolution in the low density regions and weaker planetary wakes. The disc masses after 200 orbital periods agree within 10%. The spread is larger in the tidal torques acting on the planet which agree within a factor 2 at the end of the simulation. In the Neptune case the dispersion in the torques is greater than for Jupiter, possibly owing to the contribution from the not completely cleared region close to the planet.

We study the evolution and emission of circumbinary disks around close classical T Tauri binary systems. High resolution numerical hydrodynamical simulations are employed to model a system consisting of a central eccentric binary star within an irradiated accretion disk. A detailed energy balance including viscous heating, radiative cooling and irradiation from the central star is applied to calculate accurately the emitted spectral energy distribution. Numerical simulations using two different methods, the previously developed Dual-Grid technique with a finite difference discretization, and the Smoothed Particle Hydrodynamics method are employed to compare the hydrodynamical features and strengthen our conclusions. Physical parameters of the setup are chosen to model the close systems of DQ Tau and AK Sco. Using the self-consistent models, we are able to fit the observed spectral energy distributions by constraining parameters such as disk mass, density profile and radial extension for those systems. We find that the incorporation of irradiation effects is necessary to obtain correct disk temperatures.

We have performed Smoothed Particle Hydrodynamics (SPH) simulations to study the time evolution of one and two protoplanets embedded in a protoplanetary accretion disc. We investigate accretion and migration rates of a single protoplanet depending on several parameters of the protoplanetary disc, mainly viscosity and scale height. Additionally, we consider the influence of a second protoplanet in a long time simulation and examine the migration of the two planets in the disc, especially the growth of eccentricity and chaotic behaviour. One aim of this work is to establish the feasibility of SPH for such calculations considering that usually only grid-based methods are adopted. To resolve shocks and to prevent particle penetration, we introduce a new approach for an artificial viscosity, which consists of an additional artificial bulk viscosity term in the SPH-representation of the Navier-Stokes equation. This allows for an accurate treatment of the physical kinematic viscosity to describe the shear, without the use of an artificial shear viscosity.