Dark rings with bright rims are the indirect signposts of planets embedded in protoplanetary discs. In a recent first, an azimuthally elongated AU-scale blob, possibly a planet, was resolved with ALMA in TW Hya. The blob is at the edge of a cliff-like rollover in the dust disc rather than inside a dark ring. Here we build time-dependent models of TW Hya disc. We find that the classical paradigm cannot account for the morphology of the disc and the blob. We propose that ALMA-discovered blob hides a Neptune mass planet losing gas and dust. We show that radial drift of mm-sized dust particles naturally explains why the blob is located on the edge of the dust disc. Dust particles leaving the planet perform a characteristic U-turn relative to it, producing an azimuthally elongated blob-like emission feature. This scenario also explains why a 10 Myr old disc is so bright in dust continuum. Two scenarios for the dust-losing planet are presented. In the first, a dusty pre-runaway gas envelope of about 40 Earth mass Core Accretion planet is disrupted, e.g., as a result of a catastrophic encounter. In the second, a massive dusty pre-collapse gas giant planet formed by Gravitational Instability is disrupted by the energy released in its massive core. Future modelling may discriminate between these scenarios and allow us to study planet formation in an entirely new way -- by analysing the flows of dust and gas recently belonging to planets, informing us about the structure of pre-disruption planetary envelopes.

Recent ALMA observations may indicate a surprising abundance of sub-Jovian planets on very wide orbits in protoplanetary discs that are only a few million years old. These planets are too young and distant to have been formed via the Core Accretion (CA) scenario, and are much less massive than the gas clumps born in the classical Gravitational Instability (GI) theory. It was recently suggested that such planets may form by the partial destruction of GI protoplanets: energy output due to the growth of a massive core may unbind all or most of the surrounding pre-collapse protoplanet. Here we present the first 3D global disc simulations that simultaneously resolve grain dynamics in the disc and within the protoplanet. We confirm that massive GI protoplanets may self-destruct at arbitrarily large separations from the host star provided that solid cores of mass around 10-20 Earth masses are able to grow inside them during their pre-collapse phase. In addition, we find that the heating force recently analysed by Masset and Velasco Romero (2017) perturbs these cores away from the centre of their gaseous protoplanets. This leads to very complicated dust dynamics in the protoplanet centre, potentially resulting in the formation of multiple cores, planetary satellites, and other debris such as planetesimals within the same protoplanet. A unique prediction of this planet formation scenario is the presence of sub-Jovian planets at wide orbits in Class 0/I protoplanetary discs.

Direct imaging (DI) surveys suggest that gas giants beyond 20 AU are rare around FGK stars. However, it is not clear what this means for the formation frequency of Gravitational Instability (GI) protoplanets due to uncertainties in gap opening and migration efficiency. Here we combine state-of-the-art calculations of homogeneous planet contraction with a population synthesis code. We find DI constraints to be satisfied if protoplanet formation by GI occurs in tens of percent of systems if protoplanets `super migrate' to small separations. In contrast, GI may occur in only a few percent of systems if protoplanets remain stranded at wide orbits because their migration is `quenched' by efficient gap opening. We then use the frequency of massive giants in radial velocity surveys inside 5 AU to break this degeneracy - observations recently showed that this population does not correlate with the host star metallicity and is therefore suspected to have formed via GI followed by inward migration. We find that only the super-migration scenario can sufficiently explain this population whilst simultaneously satisfying the DI constraints and producing the right mass spectrum of planets inside 5 AU. If massive gas-giants inside 5 AU formed via GI, then our models imply that migration must be efficient and that the formation of GI protoplanets occurs in at least a tens of percent of systems.

We run numerical simulations to study the accretion of gas and dust grains onto gas giant planets embedded into massive protoplanetary discs. The outcome is found to depend on the disc cooling rate, planet mass, grain size and irradiative feedback from the planet. If radiative cooling is efficient, planets accrete both gas and pebbles rapidly, open a gap and usually become massive brown dwarfs. In the inefficient cooling case, gas is too hot to accrete onto the planet but pebble accretion continues and the planets migrate inward rapidly. Radiative feedback from the planet tends to suppress gas accretion. Our simulations predict that metal enrichment of planets by dust grain accretion inversely correlates with the final planet mass, in accordance with the observed trend in the inferred bulk composition of Solar System and exosolar giant planets. To account for observations, however, as much as ~30-50% of the dust mass should be in the form of large grains.

Emission line galaxies (ELGs) are used in several ongoing and upcoming surveys (SDSS-IV/eBOSS, DESI) as tracers of the dark matter distribution. Using a new galaxy formation model, we explore the characteristics of [OII] emitters, which dominate optical ELG selections at $z\simeq 1$. Model [OII] emitters at $0.5<z<1.5$ are selected to mimic the DEEP2, VVDS, eBOSS and DESI surveys. The luminosity functions of model [OII] emitters are in reasonable agreement with observations. The selected [OII] emitters are hosted by haloes with $M_{\rm halo}\geq 10^{10.3}h^{-1}{\rm M}_{\odot}$, with ~90% of them being central star-forming galaxies. The predicted mean halo occupation distributions of [OII] emitters has a shape typical of that inferred for star-forming galaxies, with the contribution from central galaxies, $\langle N \rangle_{\left[OII\right]\, cen}$, being far from the canonical step function. The $\langle N \rangle_{\left[OII\right]\, cen}$ can be described as the sum of an asymmetric Gaussian for disks and a step function for spheroids, which plateaus below unity. The model [OII] emitters have a clustering bias close to unity, which is below the expectations for eBOSS and DESI ELGs. At $z\sim 1$, a comparison with observed g-band selected galaxy, which are expected to be dominated by [OII] emitters, indicates that our model produces too few [OII] emitters that are satellite galaxies. This suggests the need to revise our modelling of hot gas stripping in satellite galaxies.