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.

We carry out three-dimensional SPH simulations to study whether planets can survive in self-gravitating protoplanetary discs. The discs modelled here use a cooling prescription that mimics a real disc which is only gravitationally unstable in the outer regions. We do this by modelling the cooling using a simplified method such that the cooling time in the outer parts of the disc is shorter than in the inner regions, as expected in real discs. We find that both giant (> M_Sat) and low mass (< M_Nep) planets initially migrate inwards very rapidly, but are able to slow down in the inner gravitationally stable regions of the disc without needing to open up a gap. This is in contrast to previous studies where the cooling was modelled in a more simplified manner where regardless of mass, the planets were unable to slow down their inward migration. This shows the important effect the thermodynamics has on planet migration. In a broader context, these results show that planets that form in the early stages of the discs' evolution, when they are still quite massive and self-gravitating, can survive.

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.

We use a suite of SPH simulations to investigate the susceptibility of protoplanetary discs to the effects of self-gravity as a function of star-disc properties. We also include passive irradiation from the host star using different models for the stellar luminosities. The critical disc-to-star mass ratio for axisymmetry (for which we produce criteria) increases significantly for low-mass stars. This could have important consequences for increasing the potential mass reservoir in a proto Trappist-1 system, since even the efficient Ormel et al. (2017) formation model will be influenced by processes like external photoevaporation, which can rapidly and dramatically deplete the dust reservoir. The aforementioned scaling of the critical $M_d/M_*$ for axisymmetry occurs in part because the Toomre $Q$ parameter has a linear dependence on surface density (which promotes instability) and only an $M_*^{1/2}$ dependence on shear (which reduces instability), but also occurs because, for a given $M_d/M_*$, the thermal evolution depends on the host star mass. The early phase stellar irradiation of the disc (for which the luminosity is much higher than at the zero age main sequence, particularly at low stellar masses) can also play a key role in significantly reducing the role of self-gravity, meaning that even Solar mass stars could support axisymmetric discs a factor two higher in mass than usually considered possible. We apply our criteria to the DSHARP discs with spirals, finding that self-gravity can explain the observed spirals so long as the discs are optically thick to the host star irradiation.

The Transiting Exoplanet Survey Satellite, TESS, is currently carrying out an all-sky search for small planets transiting bright stars. In the first year of the TESS survey, steady progress was made in achieving the mission's primary science goal of establishing bulk densities for 50 planets smaller than Neptune. During that year, TESS's observations were focused on the southern ecliptic hemisphere, resulting in the discovery of three mini-Neptunes orbiting the star TOI-125, a V=11.0 K0 dwarf. We present intensive HARPS radial velocity observations, yielding precise mass measurements for TOI-125b, TOI-125c and TOI-125d. TOI-125b has an orbital period of 4.65 days, a radius of $2.726 \pm 0.075 ~\mathrm{R_{\rm E}}$, a mass of $ 9.50 \pm 0.88 ~\mathrm{M_{\rm E}}$ and is near the 2:1 mean motion resonance with TOI-125c at 9.15 days. TOI-125c has a similar radius of $2.759 \pm 0.10 ~\mathrm{R_{\rm E}}$ and a mass of $ 6.63 \pm 0.99 ~\mathrm{M_{\rm E}}$, being the puffiest of the three planets. TOI-125d, has an orbital period of 19.98 days and a radius of $2.93 \pm 0.17~\mathrm{R_{\rm E}}$ and mass $13.6 \pm 1.2 ~\mathrm{M_{\rm E}}$. For TOI-125b and TOI-125d we find unusual high eccentricities of $0.19\pm 0.04$ and $0.17^{+0.08}_{-0.06}$, respectively. Our analysis also provides upper mass limits for the two low-SNR planet candidates in the system; for TOI-125.04 ($R_P=1.36 ~\mathrm{R_{\rm E}}$, $P=$0.53 days) we find a $2\sigma$ upper mass limit of $1.6~\mathrm{M_{\rm E}}$, whereas TOI-125.05 ( $R_P=4.2^{+2.4}_{-1.4} ~\mathrm{R_{\rm E}}$, $P=$ 13.28 days) is unlikely a viable planet candidate with upper mass limit $2.7~\mathrm{M_{\rm E}}$. We discuss the internal structure of the three confirmed planets, as well as dynamical stability and system architecture for this intriguing exoplanet system.

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.

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.

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.

Recent observations show that rings and gaps are ubiquitous in protoplanetary discs. These features are often interpreted as being due to the presence of planets; however, the effect of planetary migration on the observed morphology has not been investigated hitherto. In this work we investigate whether multiwavelength mm/submm observations can detect signatures of planet migration, using 2D dusty hydrodynamic simulations to model the structures generated by migrating planets and synthesising ALMA continuum observations at 0.85 and 3 mm. We identify three possible morphologies for a migrating planet: a slowly migrating planet is associated with a single ring outside the planet's orbit, a rapidly migrating planet is associated with a single ring inside the planet's orbit while a planet migrating at intermediate speed generates one ring on each side of the planet's orbit. We argue that multiwavelength data can distinguish multiple rings produced by a migrating planet from other scenarios for creating multiple rings, such as multiple planets or discs with low viscosity. The signature of migration is that the outer ring has a lower spectral index, due to larger dust grains being trapped there. Of the recent ALMA observations revealing protoplanetary discs with multiple rings and gaps, we suggest that Elias 24 is the best candidate for a planet migrating in the intermediate speed regime.

We report on the confirmation and mass determination of a transiting planet orbiting the old and inactive G7 dwarf star HD219666 (Mstar = 0.92 +/- 0.03 MSun, Rstar = 1.03 +/- 0.03 RSun, tau_star = 10 +/- 2 Gyr). With a mass of Mb = 16.6 +/- 1.3 MEarth, a radius of Rb = 4.71 +/- 0.17 REarth, and an orbital period of P ~ 6 days, HD219666b is a new member of a rare class of exoplanets: the hot-Neptunes. The Transiting Exoplanet Survey Satellite (TESS) observed HD219666 (also known as TOI-118) in its Sector 1 and the light curve shows four transit-like events, equally spaced in time. We confirmed the planetary nature of the candidate by gathering precise radial-velocity measurements with HARPS@ESO3.6m. We used the co-added HARPS spectrum to derive the host star fundamental parameters (Teff = 5527 +/- 65 K, log g = 4.40 +/- 0.11 (cgs), [Fe/H]= 0.04 +/- 0.04 dex, log R'HK = -5.07 +/- 0.03), as well as the abundances of many volatile and refractory elements. The host star brightness (V = 9.9) makes it suitable for further characterisation by means of in-transit spectroscopy. The determination of the planet orbital obliquity, along with the atmospheric metal-to-hydrogen content and thermal structure could provide us with important clues on the formation mechanisms of this class of objects.

Gas clumps formed within massive gravitationally unstable circumstellar discs are potential seeds of gas giant planets, brown dwarfs and companion stars. Simulations show that competition between three processes -- migration, gas accretion and tidal disruption -- establishes what grows from a given seed. Here we investigate the robustness of numerical modelling of clump migration and accretion with the codes PHANTOM, GADGET, SPHINX, SEREN, GIZMO-MFM, SPHNG and FARGO. The test problem comprises a clump embedded in a massive disc at an initial separation of 120 AU. There is a general qualitative agreement between the codes, but the quantitative agreement in the planet migration rate ranges from $\sim 10$% to $\sim 50$%, depending on the numerical setup. We find that the artificial viscosity treatment and the sink particle prescription may account for much of the differences between the codes. In order to understand the wider implications of our work, we also attempt to reproduce the planet evolution tracks from our hydrodynamical simulations with prescriptions from three previous population synthesis studies. We find that the disagreement amongst the population synthesis models is far greater than that between our hydrodynamical simulations. The results of our code comparison project are therefore encouraging in that uncertainties in the given problem are probably dominated by the physics not yet included in the codes rather than by how hydrodynamics is modelled in them.

We present long baseline Atacama Large Millimeter/submillimeter Array (ALMA) observations of the 870$\,\mu$m dust continuum emission and CO (3-2) from the protoplanetary disk around the Herbig Ae/Be star HD 100546, which is one of the few systems claimed to have two young embedded planets. These observations achieve a resolution of 4 au (3.8 mas), an rms noise of 66$\mu$Jy/beam, and reveal an asymmetric ring between $\sim$20-40 au with largely optically thin dust continuum emission. This ring is well fit by two concentric and overlapping Gaussian rings of different widths and a Vortex. In addition, an unresolved component is detected at a position consistent with the central star, which may trace the central inner disk ($<$2au in radius). We report a lack of compact continuum emission at the positions of both claimed protoplanets. We use this result to constrain the circumplanetary disk (CPD) mass and size of 1.44M$_{\rm Earth}$ and 0.44au in the optically thin and thick regime, respectively, for the case of the previously directly imaged protoplanet candidate at $\sim$55 au (HD100546 b). We compare these empirical CPD constraints to previous numerical simulations. This suggests that HD100546 b is inconsistent with several planet accretion models, while gas-starved models are also still compatible. We estimate the planetary mass as 1.65 M$_J$ by using the relation between planet, circumstellar, and circumplanetary masses derived from numerical simulations. Finally, the CO integrated intensity map shows a possible spiral arm feature that would match the spiral features identified in Near-Infrared scattered light polarized emission, which suggests a real spiral feature in the disk surface that needs to be confirmed with further observations.

Planet migration in protoplanetary discs plays an important role in the longer term evolution of planetary systems, yet we currently have no direct observational test to determine if a planet is migrating in its gaseous disc. We explore the formation and evolution of dust rings - now commonly observed in protoplanetary discs by ALMA - in the presence of relatively low mass (12-60 Earth masses) migrating planets. Through two dimensional hydrodynamical simulations using gas and dust we find that the importance of perturbations in the pressure profile interior and exterior to the planet varies for different particle sizes. For small sizes a dust enhancement occurs interior to the planet, whereas it is exterior to it for large particles. The transition between these two behaviours happens when the dust drift velocity is comparable to the planet migration velocity. We predict that an observational signature of a migrating planet consists of a significant outwards shift of an observed midplane dust ring as the wavelength is increased.

We present high resolution millimeter continuum imaging of the disc surrounding the young star CI Tau, a system hosting the first hot Jupiter candidate in a protoplanetary disc system. The system has extended mm emission on which are superposed three prominent annular gaps at radii ~ 13, 39 and 100 au. We argue that these gaps are most likely to be generated by massive planets so that, including the hot Jupiter, the system contains four gas giant planets at an age of only 2 Myr. Two of the new planets are similarly located to those inferred in the famous HL Tau protoplanetary disc; in CI Tau, additional observational data enables a more complete analysis of the system properties than was possible for HL Tau. Our dust and gas dynamical modeling satisfies every available observational constraint and points to the most massive ensemble of exo-planets ever detected at this age, with its four planets spanning a factor 1000 in orbital radius.Our results show that the association between hot Jupiters and gas giants on wider orbits, observed in older stars, is apparently in place at an early evolutionary stage.

The spiral waves detected in the protostellar disc surrounding Elias 2-27 have been suggested as evidence of the disc being gravitationally unstable. However, previous work has shown that a massive, stable disc undergoing an encounter with a massive companion are also consistent with the observations. We compare the spiral morphology of smoothed particle hydrodynamic simulations modelling both cases. The gravitationally unstable disc produces symmetric, tightly wound spiral arms with constant pitch angle, as predicted by the literature. The companion disc's arms are asymmetric, with pitch angles that increase with radius. However, these arms are not well-fitted by standard analytic expressions, due to the high disc mass and relatively low companion mass. We note that differences (or indeed similarities) in morphology between pairs of spirals is a crucial discriminant between scenarios for Elias 2-27, and hence future studies must fit spiral arms individually. If Elias 2-27 continues to show symmetric tightly wound spiral arms in future observations, then we posit that it is the first observed example of a gravitationally unstable protostellar disc.

For grain growth to proceed effectively and lead to planet formation a number of barriers to growth must be overcome. One such barrier, relevant for compact grains in the inner regions of the disc, is the `bouncing barrier' in which large grains ($\sim$ mm size) tend to bounce off each other rather than sticking. However, by maintaining a population of small grains it has been suggested that cm-size particles may grow rapidly by sweeping up these small grains. We present the first numerically resolved investigation into the conditions under which grains may be lucky enough to grow beyond the bouncing barrier by a series of rare collisions leading to growth (so-called `breakthrough'). Our models support previous results, and show that in simple models breakthrough requires the mass ratio at which high velocity collisions transition to growth instead of causing fragmentation to be low, $\phi \lesssim 50$. However, in models that take into account the dependence of the fragmentation threshold on mass-ratio, we find breakthrough occurs more readily, even if mass transfer is relatively inefficient. This suggests that bouncing may only slow down growth, rather than preventing growth beyond a threshold barrier. However, even when growth beyond the bouncing barrier is possible, radial drift will usually prevent growth to arbitrarily large sizes.

We carry out simulations of gravitationally unstable disks using smoothed particle hydrodynamics(SPH) and the novel Lagrangian meshless finite mass (MFM) scheme in the GIZMO code (Hopkins 2015). Our aim is to understand the cause of the non-convergence of the cooling boundary for fragmentation reported in the literature. We run SPH simulations with two different artificial viscosity implementations, and compare them with MFM, which does not employ any artificial viscosity. With MFM we demonstrate convergence of the critical cooling time scale for fragmentation at \beta_{crit} =3.. Non-convergence persists in SPH codes, although it is significantly mitigated with schemes having reduced artificial viscosity such as inviscid SPH (ISPH) (Cullen & Dehnen 2010). We show how the non-convergence problem is caused by artificial fragmentation triggered by excessive dissipation of angular momentum in domains with large velocity derivatives. With increased resolution such domains become more prominent. Vorticity lags behind density due to numerical viscous dissipation in these regions, promoting collapse with longer cooling times. Such effect is shown to be dominant over the competing tendency of artificial viscosity to diminish with increasing resolution. When the initial conditions are first relaxed for several orbits, the flow is more regular, with lower shear and vorticity in non-axisymmetric regions, aiding convergence. Yet MFM is the only method that converges exactly. Our findings are of general interest as numerical dissipation via artificial viscosity or advection errors can also occur in grid-based codes. Indeed for the FARGO code values of \beta_{crit} significantly higher than our converged estimate have been reported in the literature. Finally, we discuss implications for giant planet formation via disk instability.

It is likely that most protostellar systems undergo a brief phase where the protostellar disc is self-gravitating. If these discs are prone to fragmentation, then they are able to rapidly form objects that are initially of several Jupiter masses and larger. The fate of these disc fragments (and the fate of planetary bodies formed afterwards via core accretion) depends sensitively not only on the fragment's interaction with the disc, but with its neighbouring fragments. We return to and revise our population synthesis model of self-gravitating disc fragmentation and tidal downsizing. Amongst other improvements, the model now directly incorporates fragment-fragment interactions while the disc is still present. We find that fragment-fragment scattering dominates the orbital evolution, even when we enforce rapid migration and inefficient gap formation. Compared to our previous model, we see a small increase in the number of terrestrial-type objects being formed, although their survival under tidal evolution is at best unclear. We also see evidence for disrupted fragments with evolved grain populations - this is circumstantial evidence for the formation of planetesimal belts, a phenomenon not seen in runs where fragment-fragment interactions are ignored. In spite of intense dynamical evolution, our population is dominated by massive giant planets and brown dwarfs at large semimajor axis, which direct imaging surveys should, but only rarely, detect. Finally, disc fragmentation is shown to be an efficient manufacturer of free floating planetary mass objects, and the typical multiplicity of systems formed via gravitational instability will be low.

The young star Elias 2-27 has recently been observed to posses a massive circumstellar disc with two prominent large-scale spiral arms. In this Letter we perform three-dimensional Smoothed Particle Hydrodynamics simulations, radiative transfer modelling, synthetic ALMA imaging and an unsharped masking technique to explore three possibilities for the origin of the observed structures -- an undetected companion either internal or external to the spirals, and a self-gravitating disc. We find that a gravitationally unstable disc and a disc with an external companion can produce morphology that is consistent with the observations. In addition, for the latter, we find that the companion could be a relatively massive planetary mass companion (less than approximately 10 - 13 MJup) and located at large radial distances (between approximately 300 - 700 au). We therefore suggest that Elias 2-27 may be one of the first detections of a disc undergoing gravitational instabilities, or a disc that has recently undergone fragmentation to produce a massive companion.

The Protoplanetary Discussions conference --- held in Edinburgh, UK, from 7th --11th March 2016 --- included several open sessions led by participants. This paper reports on the discussions collectively concerned with the multiphysics modelling of protoplanetary discs, including the self-consistent calculation of gas and dust dynamics, radiative transfer and chemistry. After a short introduction to each of these disciplines in isolation, we identify a series of burning questions and grand challenges associated with their continuing development and integration. We then discuss potential pathways towards solving these challenges, grouped by strategical, technical and collaborative developments. This paper is not intended to be a review, but rather to motivate and direct future research and collaboration across typically distinct fields based on \textit{community driven input}, to encourage further progress in our understanding of circumstellar and protoplanetary discs.

The Planet Formation Imager (PFI) project aims to provide a strong scientific vision for ground-based optical astronomy beyond the upcoming generation of Extremely Large Telescopes. We make the case that a breakthrough in angular resolution imaging capabilities is required in order to unravel the processes involved in planet formation. PFI will be optimised to provide a complete census of the protoplanet population at all stellocentric radii and over the age range from 0.1 to about 100 Myr. Within this age period, planetary systems undergo dramatic changes and the final architecture of planetary systems is determined. Our goal is to study the planetary birth on the natural spatial scale where the material is assembled, which is the "Hill Sphere" of the forming planet, and to characterise the protoplanetary cores by measuring their masses and physical properties. Our science working group has investigated the observational characteristics of these young protoplanets as well as the migration mechanisms that might alter the system architecture. We simulated the imprints that the planets leave in the disk and study how PFI could revolutionise areas ranging from exoplanet to extragalactic science. In this contribution we outline the key science drivers of PFI and discuss the requirements that will guide the technology choices, the site selection, and potential science/technology tradeoffs.

We carry out three dimensional radiation hydrodynamical simulations of gravitationally unstable discs to explore the movement of mass in a disc following its initial fragmentation. We find that the radial velocity of the gas in some parts of the disc increases by up to a factor of approximately 10 after the disc fragments, compared to before. While the movement of mass occurs in both the inward and outward directions, the inwards movement can cause the inner spirals of a self-gravitating disc to become sufficiently dense such that they can potentially fragment. This suggests that the dynamical behaviour of fragmented discs may cause subsequent fragmentation to occur at smaller radii than initially expected, but only after an initial fragment has formed in the outer disc.

We perform two-dimensional hydrodynamical simulations to quantitatively explore the torque balance criterion for gap-opening (as formulated by Crida et al. 2006) in a variety of disks when considering a migrating planet. We find that even when the criterion is satisfied, there are instances when planets still do not open gaps. We stress that gap-opening is not only dependent on whether a planet has the ability to open a gap, but whether it can do so quickly enough. This can be expressed as an additional condition on the gap-opening timescale versus the crossing time, i.e. the time it takes the planet to cross the region which it is carving out. While this point has been briefly made in the previous literature, our results quantify it for a range of protoplanetary disk properties and planetary masses, demonstrating how crucial it is for gap-opening. This additional condition has important implications for the survival of planets formed by core accretion in low mass disks as well as giant planets or brown dwarfs formed by gravitational instability in massive disks. It is particularly important for planets with intermediate masses susceptible to Type III-like migration. For some observed transition disks or disks with gaps, we expect that estimates on the potential planet masses based on the torque balance gap-opening criterion alone may not be sufficient. With consideration of this additional timescale criterion theoretical studies may find a reduced planet survivability or that planets may migrate further inwards before opening a gap.

We propose a mechanism by which dust rings in protoplanetary disks can form and be long-lasting compared to gas rings. This involves the existence of a pressure maximum which traps dust either in between two gap-opening planets or at the outermost gap edge of a single or multiple planet system, combined with the decoupling of large dust particles from the gas. We perform 2D gas hydrodynamical simulations of disks with one and two giant planets which may open deep or partial gaps. A gas ring forms in between two planets such that the surface mass density is higher than on either side of it. This ring is a region of pressure maximum where we expect large grains, which are marginally coupled to the gas and would otherwise be subject to radial drift, to collect. Such a pressure maximum also occurs at the outermost gap edge in a disk with one or more planets. We infer the dust evolution in these regions as the gas disk evolves, to understand the longer term behavior of the resulting dust rings. Over time the gas surface density in the ring(s) decreases, which may cause the larger trapped particles to decouple. Consequently, these particles are expected to stay in ring structure(s) longer than the gas. For a Minimum Mass Solar Nebula model, we expect that millimeter and centimeter-sized grains in the outer O(10) au would most likely undergo this trapping and decoupling process.

We present L' and J-band high-contrast observations of HD169142, obtained with the VLT/NACO AGPM vector vortex coronagraph and the Gemini Planet Imager, respectively. A source located at 0".156+/-0".032 north of the host star (PA=7.4+/-11.3 degrees) appears in the final reduced L' image. At the distance of the star (~145 pc), this angular separation corresponds to a physical separation of 22.7+/-4.7 AU, locating the source within the recently resolved inner cavity of the transition disk. The source has a brightness of L'=12.2+/-0.5 mag, whereas it is not detected in the J band (J>13.8 mag). If its L' brightness arose solely from the photosphere of a companion and given the J-L' color constraints, it would correspond to a 28-32 MJupiter object at the age of the star, according to the COND models. Ongoing accretion activity of the star suggests, however, that gas is left in the inner disk cavity from which the companion could also be accreting. In this case the object could be lower in mass and its luminosity enhanced by the accretion process and by a circumplanetary disk. A lower mass object is more consistent with the observed cavity width. Finally, the observations enable us to place an upper limit on the L'-band flux of a second companion candidate orbiting in the disk annular gap at ~50 AU, as suggested by millimeter observations. If the second companion is also confirmed, HD169142 might be forming a planetary system, with at least two companions opening gaps and possibly interacting with each other.

The disk around the Herbig Ae/Be star HD 100546 has been extensively studied and it is one of the systems for which there are observational indications of ongoing and/or recent planet formation. However, up until now no resolved image of the millimeter dust emission or the gas has been published. We present the first resolved images of the disk around HD 100546 obtained in Band 7 with the ALMA observatory. The CO (3-2) image reveals a gas disk that extends out to 350 au radius at the 3-sigma level. Surprisingly, the 870um dust continuum emission is compact (radius <60 au) and asymmetric. The dust emission is well matched by a truncated disk with outer radius of $\approx$50 au. The lack of millimeter-sized particles outside the 60 au is consistent with radial drift of particles of this size. The protoplanet candidate, identified in previous high-contrast NACO/VLT L' observations, could be related to the sharp outer edge of the millimeter-sized particles. Future higher angular resolution ALMA observations are needed to determine the detailed properties of the millimeter emission and the gas kinematics in the inner region (<2arcsec). Such observations could also reveal the presence of a planet through the detection of circumplanetary disk material.

Transitional disks represent a short stage of the evolution of circumstellar material. Studies of dust grains in these objects can provide pivotal information on the mechanisms of planet formation. Dissimilarities in the spatial distribution of small (micron-size) and large (millimeter-size) dust grains have recently been pointed out. Constraints on the small dust grains can be obtained by imaging the distribution of scattered light at near-infrared wavelengths. We aim at resolving structures in the surface layer of transitional disks (with particular emphasis on the inner 10 - 50 AU), thus increasing the scarce sample of high resolution images of these objects. We obtained VLT/NACO near-IR high-resolution polarimetric differential imaging observations of SAO 206462 (HD135344B). This technique allows one to image the polarized scattered light from the disk without any occulting mask and to reach an inner working angle of 0.1''. A face-on disk is detected in H and Ks bands between 0.1'' and 0.9''. No significant differences are seen between the H and Ks images. In addition to the spiral arms, these new data allow us to resolve for the first time an inner cavity for small dust grains. The cavity size (about 28 AU) is much smaller than what is inferred for large dust grains from (sub)mm observations (39 to 50 AU). The interaction between the disk and potential orbiting companion(s) can explain both the spiral arm structure and the discrepant cavity sizes for small and large dust grains. One planet may be carving out the gas (and, thus, the small grains) at 28 AU, and generating a pressure bump at larger radii (39 AU), which holds back the large grains. We analytically estimate that, in this scenario, a single giant planet (with a mass between 5 and 15 Jupiter masses) at 17 to 20 AU from the star is consistent with the observed cavity sizes.

The large number of detected giant exoplanets offers the opportunity to improve our understanding of the formation mechanism, evolution, and interior structure of gas giant planets. The two main models for giant planet formation are core accretion and disk instability. There are substantial differences between these formation models, including formation timescale, favorable formation location, ideal disk properties for planetary formation, early evolution, planetary composition, etc. First, we summarize the two models including their substantial differences, advantages, and disadvantages, and suggest how theoretical models should be connected to available (and future) data. We next summarize current knowledge of the internal structures of solar- and extrasolar- giant planets. Finally, we suggest the next steps to be taken in giant planet exploration.

We perform coagulation & fragmentation simulations to understand grain growth in T Tauri & brown dwarf discs. We present a physically-motivated approach using a probability distribution function for the collision velocities and separating the deterministic & stochastic velocities. We find growth to larger sizes compared to other models. Furthermore, if brown dwarf discs are scaled-down versions of T Tauri discs (in terms of stellar & disc mass, and disc radius), growth at the same location with respect to the outer edge occurs to similar sizes in both discs.

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 coagulation and fragmentation simulations using the new physically-motivated model by Garaud et al. (2013) to determine growth locally in brown dwarf disks. We show that large grains can grow and that if brown dwarf disks are scaled down versions of T Tauri disks (in terms of stellar mass, disk mass and disk radius) growth at an equivalent location with respect to the disk truncation radius can occur to the same size in both disks. We show that similar growth occurs because the collisional timescales in the two disks are comparable. Our model may therefore potentially explain the recent observations of grain growth to millimetre sizes in brown dwarf disks, as seen in T Tauri disks.

We carry out three dimensional radiation hydrodynamical simulations of gravitationally unstable discs using to explore the movement of mass in a disc following its fragmentation. Compared to a more quiescent state before it fragments, the radial velocity of the gas increases by up to a factor of ~2-3 after fragmentation. While the mass movement occurs both inwards and outwards, the inwards motion can cause the inner spirals to be sufficiently dense that they may become unstable and potentially fragment. Consequently, the dynamical behaviour of fragmented discs may cause subsequent fragmentation at smaller radii after an initial fragment has formed in the outer disc.

Planet formation occurs within the gas and dust rich environments of protoplanetary disks. Observations of these objects show that the growth of primordial sub micron sized particles into larger aggregates occurs at the earliest stages of the disks. However, theoretical models of particle growth that use the Smoluchowski equation to describe collisional coagulation and fragmentation have so far failed to produce large particles while maintaining a significant populations of small grains. This has been generally attributed to the existence of two barriers impeding growth due to bouncing and fragmentation of colliding particles. In this paper, we demonstrate that the importance of these barriers has been artificially inflated through the use of simplified models that do not take into account the stochastic nature of the particle motions within the gas disk. We present a new approach in which the relative velocities between two particles is described by a probability distribution function that models both deterministic motion and stochastic motion. Taking both into account can give quite different results to what has been considered recently in other studies. We demonstrate the vital effect of two "ingredients" for particle growth: the proper implementation of a velocity distribution function that overcomes the bouncing barrier and, in combination with mass transfer in high-mass-ratio collisions, boosts the growth of larger particles beyond the fragmentation barrier. A robust result of our simulations is the emergence of two particle populations (small and large), potentially explaining simultaneously a number of long-standing problems in protoplanetary disks, including planetesimal formation close to the central star, the presence of mm to cm size particles far out in the disk, and the persistence of micron-size grains for millions of years.

We carry out simulations of gravitationally unstable discs using a Smoothed Particle Hydrodynamics (SPH) code and a grid-based hydrodynamics code, FARGO, to understand the previous non-convergent results reported by Meru & Bate (2011a). We obtain evidence that convergence with increasing resolution occurs with both SPH and FARGO and in both cases we find that the critical cooling timescale is larger than previously thought. We show that SPH has a first-order convergence rate while FARGO converges with a second-order rate. We show that the convergence of the critical cooling timescale for fragmentation depends largely on the numerical viscosity employed in both SPH and FARGO. With SPH, particle velocity dispersion may also play a role. We show that reducing the dissipation from the numerical viscosity leads to larger values of the critical cooling time at a given resolution. For SPH, we find that the effect of the dissipation due to the numerical viscosity is somewhat larger than had previously been appreciated. In particular, we show that using a quadratic term in the SPH artificial viscosity (beta_{SPH}) that is too low appears to lead to excess dissipation in gravitationally unstable discs, which may affect any results that sensitively depend on the thermodynamics, such as disc fragmentation. We show that the two codes converge to values of the critical cooling timescale, beta_{crit} > 20 (for a ratio of specific heats of gamma=5/3), and perhaps even as large as beta_{crit} \approx 30. These are approximately 3-5 times larger than has been found by most previous studies. This is equivalent to a maximum gravitational stress that a disc can withstand without fragmenting of alpha_{GI,crit} \approx 0.013-0.02, which is much smaller than the values typically used in the literature. It is therefore easier for self-gravitating discs to fragment than has been concluded from most past studies.

In 2D simulations of thin gaseous disks with embedded planets or self-gravity the gravitational potential needs to be smoothed to avoid singularities in the numerical evaluation of the gravitational potential or force. In order to correctly resemble the realistic case of vertically extended 3D disks the softening prescription used in 2D needs to be adjusted properly. In this paper we analyze the embedded planet and the self-gravity case and provide a method to evaluate the required smoothing in 2D simulations of thin disks. Starting from the averaged hydrodynamic equations and using a vertically isothermal disk model, we calculate the force. We compare our results to the often used Plummer form of the potential which runs as \propto 1/(r^2+\epsilon^2)^{1/2}. For that purpose we compute the required smoothing length \epsilon as a function of distance r to the planet or to a disk element within a self-gravitating disk. We find that for larger distances \epsilon is determined solely by the vertical disk thickness H. For the planet case we find that outside r\simH a value of \epsilon = 0.7H describes the averaged force very well, while in the self-gravitating disk the value needs to be larger, \epsilon = 1.2H. For smaller distances the smoothing needs to be reduced significantly. Comparing torque densities of 3D and 2D simulations we show that the modification to the vertical density stratification as induced by an embedded planet needs to be taken into account to obtain agreeing results. In disk fragmentation simulations the choice of \epsilon can determine whether a disk will fragment or not. Because a wrong smoothing length can change even the direction of migration, it is very important to include the effect of the planet on the local scale height in 2D planet-disk simulations. We provide an approximate and fast method for this purpose which gives very good agreement to full 3D simulations.

We study the numerical convergence of hydrodynamical simulations of self-gravitating accretion discs, in which a simple cooling law is balanced by shock heating. It is well-known that there exists a critical cooling time scale for which shock heating can no longer compensate for the energy losses, at which point the disc fragments. The numerical convergence of previous results of this critical cooling time scale was questioned recently using Smoothed Particle Hydrodynamics (SPH). We employ a two-dimensional grid-based code to study this problem, and find that for smooth initial conditions, fragmentation is possible for slower cooling as the resolution is increased, in agreement with recent SPH results. We show that this non-convergence is at least partly due to the creation of a special location in the disc, the boundary between the turbulent and the laminar region, when cooling towards a gravito-turbulent state. Converged results appear to be obtained in setups where no such sharp edges appear, and we then find a critical cooling time scale of ~ 4 $\Omega^{-1}$, where $\Omega$ is the local angular velocity.

The observation of massive exoplanets at large separation from their host star, like in the HR 8799 system, challenges theories of planet formation. A possible formation mechanism involves the fragmentation of massive self-gravitating discs into clumps. While the conditions for fragmentation have been extensively studied, little is known of the subsequent evolution of these giant planet embryos, in particular their expected orbital migration. Assuming a single planet has formed by fragmentation, we investigate its interaction with the gravitoturbulent disc it is embedded in. Two-dimensional hydrodynamical simulations are used with a simple prescription for the disc cooling. A steady gravitoturbulent disc is first set up, after which simulations are restarted including a planet with a range of masses approximately equal to the clump's initial mass expected in fragmenting discs. Planets rapidly migrate inwards, despite the stochastic kicks due to the turbulent density fluctuations. We show that the migration timescale is essentially that of type I migration, with the planets having no time to open a gap. In discs with aspect ratio ~ 0.1 at their forming location, planets with a mass comparable to, or larger than Jupiter's can migrate in as short as 10000 years, that is, about 10 orbits at 100 AU. Massive planets formed at large separation from their star by gravitational instability are thus unlikely to stay in place, and should rapidly migrate towards the inner parts of protoplanetary discs, regardless of the planet mass.

Within the collision growth scenario for planetesimal formation, the growth step from centimetre sized pre-planetesimals to kilometre sized planetesimals is still unclear. The formation of larger objects from the highly porous pre-planetesimals may be halted by a combination of fragmentation in disruptive collisions and mutual rebound with compaction. However, the right amount of fragmentation is necessary to explain the observed dust features in late T Tauri discs. Therefore, detailed data on the outcome of pre-planetesimal collisions is required and has to be presented in a suitable and precise format. We propose and apply a new classification scheme for pre-planetesimal collisions based on the quantitative aspects of four fragment populations: the largest and second largest fragment, a power-law population, and a sub-resolution population. For the simulations of pre-planetesimal collisions, we adopt the SPH numerical scheme with extensions for the simulation of porous solid bodies. By means of laboratory benchmark experiments, this model was previously calibrated and tested for the correct simulation of the compaction, bouncing, and fragmentation behaviour of macroscopic highly porous silica dust aggregates. It is shown that previous attempts to map collision data were much too oriented on qualitatively categorising into sticking, bouncing, and fragmentation events. We show that the four-populations model encompasses all previous categorisations and in addition allows for transitions. This is because it is based on quantitative characteristic attributes of each population such as the mass, kinetic energy, and filling factor. As a demonstration of the applicability and the power of the four-populations model, we utilise it to present the results of a study on the influence of collision velocity in head-on collisions of intermediate porosity aggregates.

We carry out a resolution study on the fragmentation boundary of self-gravitating discs. We perform three-dimensional Smoothed Particle Hydrodynamics (SPH) simulations of discs to determine whether the critical value of the cooling timescale in units of the orbital timescale, beta_{crit}, converges with increasing resolution. Using particle numbers ranging from 31,250 to 16 million (the highest resolution simulations to date) we do not find convergence. Instead, fragmentation occurs for longer cooling timescales as the resolution is increased. These results certainly suggest that beta_{crit} is larger than previously thought. However, the absence of convergence also questions whether or not a critical value exists. In light of these results, we caution against using cooling timescale or gravitational stress arguments to deduce whether gravitational instability may or may not have been the formation mechanism for observed planetary systems.

We carry out a resolution study on the fragmentation boundary of self-gravitating discs. We perform three-dimensional Smoothed Particle Hydrodynamics simulations of discs to determine whether the critical value of the cooling timescale in units of the orbital timescale, beta_{crit}, converges with increasing resolution. Using particle numbers ranging from 31,250 to 16 million (the highest resolution simulations to date) we do not find convergence. Instead, fragmentation occurs for longer cooling timescales as the resolution is increased. These results suggest that at the very least, the critical value of the cooling timescale is longer than previously thought. However, the absence of convergence also raises the question of whether or not a critical value exists. In light of these results, we caution against using cooling timescale or gravitational stress arguments to deduce whether gravitational instability may or may not have been the formation mechanism for observed planetary systems.

We investigate the fragmentation criterion in massive self-gravitating discs. We present new analysis of the fragmentation conditions which we test by carrying out global three-dimensional numerical simulations. Whilst previous work has placed emphasis on the cooling timescale in units of the orbital timescale, \beta , we find that at a given radius the surface mass density (i.e. disc mass and profile) and star mass also play a crucial role in determining whether a disc fragments or not as well as where in the disc fragments form. We find that for shallow surface mass density profiles (p<2, where \Sigma \propto R^{-p}), fragments form in the outer regions of the disc. However for steep surface mass density profiles (p \gtrsim 2), fragments form in the inner regions of a disc. In addition, we also find that the critical value of the cooling timescale in units of the orbital timescale found in previous simulations is only applicable to certain disc surface mass density profiles and for particular disc radii and is not a general rule for all discs. We find an empirical fragmentation criteria between the cooling timescale in units of the orbital timescale, \beta , the surface mass density, the star mass and the radius.

We carry out global three-dimensional radiation hydrodynamical simulations of self-gravitating accretion discs to determine if, and under what conditions, a disc may fragment to form giant planets. We explore the parameter space (in terms of the disc opacity, temperature and size) and include the effect of stellar irradiation. We find that the disc opacity plays a vital role in determining whether a disc fragments. Specifically, opacities that are smaller than interstellar Rosseland mean values promote fragmentation (even at small radii, R < 25AU) since low opacities allow a disc to cool quickly. This may occur if a disc has a low metallicity or if grain growth has occurred. With specific reference to the HR 8799 planetary system, given its star is metal-poor, our results suggest that the formation of its imaged planetary system could potentially have occurred by gravitational instability. We also find that the presence of stellar irradiation generally acts to inhibit fragmentation (since the discs can only cool to the temperature defined by stellar irradiation). However, fragmentation may occur if the irradiation is sufficiently weak that it allows the disc to attain a low Toomre stability parameter.

We present preliminary results on the ability of self-gravitating discs to cool in response to their internal heating. These discs are modelled using a Smoothed Particle Hydrodynamics (SPH) code with radiative transfer (Whitehouse, Bate & Monaghan 2005) and we investigate the ability of these discs to maintain a state of thermal equilibrium with their boundaries as an indication of their likelihood to fragment.

In this paper, we use high-resolution smoothed particle hydrodynamics (SPH) simulations to investigate the response of a marginally stable self-gravitating protostellar disc to a close parabolic encounter with a companion discless star. Our main aim is to test whether close brown dwarfs or massive planets can form out of the fragmentation of such discs. We follow the thermal evolution of the disc by including the effects of heating due to compression and shocks and a simple prescription for cooling and find results that contrast with previous isothermal simulations. In the present case we find that fragmentation is inhibited by the interaction, due to the strong effect of tidal heating, which results in a strong stabilization of the disc. A similar behaviour was also previously observed in other simulations involving discs in binary systems. As in the case of isolated discs, it appears that the condition for fragmentation ultimately depends on the cooling rate.

We present measurements of the effective temperatures and surface gravities for a sample of hot DA white dwarfs, using the Lyman line data available from the HUT, ORFEUS and FUSE far-UV space missions. Comparing the results with those from the standard Balmer line technique, we find that there is a general good overall agreement between the two methods. However, significant differences are found for a number of stars, but not always of a consistent nature in that sometimes the Balmer temperature exceeds that derived from the Lyman lines and in other instances is lower. We conclude that, with the latest model atmosphere calculations, these discrepancies probably do not arise from an inadequate theoretical treatment of the Lyman lines but rather from systematic effects in the observation and data reduction processes, which dominate the statistical errors in these spectra. If these systematic data reduction effects can be adequately controlled, the Lyman line temperature and gravity measurements are consistent with those obtained from the Balmer lines when allowance is made for reasonable observational uncertainties.