To understand the formation of planetary systems, one needs to understand the initial conditions of planet formation, i.e. the young gas-rich planet forming disks. Spatially resolved high-contrast observations are of particular interest, since substructures in disks, linked to planet formation, can be detected and close companions or even planets in formation embedded in the disk can be revealed. In this study we present the first result of the DESTINYS survey (Disk Evolution Study Through Imaging of Nearby Young Stars). DESTINYS is an ESO/SPHERE large program that aims at studying disk evolution in scattered light, mainly focusing on a sample of low-mass stars (<1$M_\odot$) in nearby (~200 pc) star-forming regions. In this particular study we present the observations of the ET Cha (RECX 15) system, a nearby 'old' classical T Tauri star (5-8 Myr, ~100 pc), which is still strongly accreting. We use SPHERE/IRDIS in H-band polarimetric imaging mode to obtain high contrast images of the ET Cha system to search for scattered light from the circumstellar disk as well as thermal emission from close companions. We additionally employ VLT/NACO total intensity archival data taken in 2003. We report here the discovery of a low-mass (sub)stellar companion with SPHERE/IRDIS to ET Cha. We are estimating the mass of this new companion based on photometry. Depending on the system age it is a 5 Myr, 50 $M_{Jup}$ brown dwarf or an 8 Myr, 0.10 $M_\odot$ M-type pre-main-sequence star. We explore possible orbital solutions and discuss the recent dynamic history of the system. Independent of the precise companion mass we find that the presence of the companion likely explains the small size of the disk around ET Cha. The small separation of the binary pair indicates that the disk around the primary component is likely clearing from the outside in, explaining the high accretion rate of the system.

Disc-driven migration is a key evolutionary stage of supermassive black-hole binaries hosted in gas-rich galaxies. Besides promoting the inspiral, viscous interactions tend to align the spins of the black holes with the orbital angular momentum of the disc. We present a critical and systematic investigation of this problem, also known as the Bardeen-Petterson effect. We design a new iterative scheme to solve the non-linear dynamics of warped accretion discs under the influence of both relativistic frame dragging and binary companion. We characterize the impact of the disc "critical obliquity", which marks regions of the parameter space where stationary solutions do not exist. We find that black-hole spins reach either complete alignment or a critical configuration. Reaching the critical obliquity might imply that the disc breaks as observed in hydrodynamical simulations. Our findings are important to predict the spin configurations with which supermassive black-hole binaries enter their gravitational-wave driven regime and become detectable by LISA.

We report on the discovery of periodic dips in the multiband lightcurve of ISO-ChaI 52, a young stellar object in the Chamaeleon I dark cloud. This is one among the peculiar objects that display very low or negligible accretion both in their UV continuum and spectral lines, although they present a remarkable infrared excess emission characteristic of optically-thick circumstellar disks. We have analyzed a VLT/X-Shooter spectrum with the tool ROTFIT to determine the stellar parameters. The latter, along with photometry from our campaign with the REM telescope and from the literature, have allowed us to model the spectral energy distribution and to estimate the size and temperature of the inner and outer disk. From the rotational period of the star/disk system of 3.45 days we estimate a disk inclination of 36$^\circ$. The depth of the dips in different bands has been used to gain information about the occulting material. A single extinction law is not able to fit the observed behavior, while a two-component model of a disk warp composed of a dense region with a gray extinction and an upper layer with an ISM-type extinction provides a better fit of the data.

In recent years evidence has been building that planet formation starts early, in the first $\sim$ 0.5 Myr. Studying the dust masses available in young disks enables understanding the origin of planetary systems since mature disks are lacking the solid material necessary to reproduce the observed exoplanetary systems, especially the massive ones. We aim to determine if disks in the embedded stage of star formation contain enough dust to explain the solid content of the most massive exoplanets. We use Atacama Large Millimeter/submillimeter Array (ALMA) Band 6 observations of embedded disks in the Perseus star-forming region together with Very Large Array (VLA) Ka-band (9 mm) data to provide a robust estimate of dust disk masses from the flux densities. Using the DIANA opacity model including large grains, with a dust opacity value of $\kappa_{\rm 9\ mm}$ = 0.28 cm$^{2}$ g$^{-1}$, the median dust masses of the embedded disks in Perseus are 158 M$_\oplus$ for Class 0 and 52 M$_\oplus$ for Class I from the VLA fluxes. The lower limits on the median masses from ALMA fluxes are 47 M$_\oplus$ and 12 M$_\oplus$ for Class 0 and Class I, respectively, obtained using the maximum dust opacity value $\kappa_{\rm 1.3mm}$ = 2.3 cm$^{2}$ g$^{-1}$. The dust masses of young Class 0 and I disks are larger by at least a factor of 10 and 3, respectively, compared with dust masses inferred for Class II disks in Lupus and other regions. The dust masses of Class 0 and I disks in Perseus derived from the VLA data are high enough to produce the observed exoplanet systems with efficiencies acceptable by planet formation models: the solid content in observed giant exoplanets can be explained if planet formation starts in Class 0 phase with an efficiency of $\sim$ 15%. Higher efficiency of $\sim$ 30% is necessary if the planet formation is set to start in Class I disks.

Interferometric observations of the mm dust distribution in protoplanetary discs are now showing a ubiquity of annular gap and ring substructures. Their identification and accurate characterization is critical to probing the physical processes responsible. We present Frankenstein (frank), an open source code that recovers axisymmetric disc structures at sub-beam resolution. By fitting the visibilities directly, the model reconstructs a disc's 1D radial brightness profile nonparametrically using a fast (<~1 min) Gaussian process. The code avoids limitations of current methods that obtain the radial brightness profile by either extracting it from the disc image via nonlinear deconvolution at the cost of reduced fit resolution, or by assumptions placed on the functional forms of disc structures to fit the visibilities parametrically. We use mock ALMA observations to quantify the method's intrinsic capability and its performance as a function of baseline-dependent signal-to-noise. Comparing the technique to profile extraction from a CLEAN image, we motivate how our fits accurately recover disc structures at a sub-beam resolution. Demonstrating the model's utility in fitting real high and moderate resolution observations, we conclude by proposing applications to address open questions on protoplanetary disc structure and processes.

The evolution of protoplanetary disks is dominated by the conservation of angular momentum, where the accretion of material onto the central star is driven by viscous expansion of the outer disk or by disk winds extracting angular momentum without changing the disk size. Studying the time evolution of disk sizes allows us therefore to distinguish between viscous stresses or disk winds as the main mechanism of disk evolution. Observationally, estimates of the disk gaseous outer radius are based on the extent of the CO rotational emission, which, during the evolution, is also affected by the changing physical and chemical conditions in the disk. We use physical-chemical DALI models to study how the extent of the CO emission changes with time in a viscously expanding disk and investigate to what degree this observable gas outer radius is a suitable tracer of viscous spreading and whether current observations are consistent with viscous evolution. We find that the gas outer radius (R_co) measured from our models matches the expectations of a viscously spreading disk: R_co increases with time and for a given time R_co is larger for a disk with a higher viscosity alpha_visc. However, in the extreme case where the disk mass is low (less than 10^-4 Msun) and alpha_visc is high (larger than 10^-2), R_co will instead decrease with time as a result of CO photodissociation in the outer disk. For most disk ages R_co is up to 12x larger than the characteristic size R_c of the disk, and R_co/R_c is largest for the most massive disk. As a result of this difference, a simple conversion of R_co to alpha_visc will overestimate the true alpha_visc of the disk by up to an order of magnitude. We find that most observed gas outer radii in Lupus can be explained using a viscously evolving disk that starts out small (R_c = 10 AU) and has a low viscosity (alpha_visc = 10^-4 - 10^-3).

Determining the mechanisms that drive the evolution of protoplanetary disks is a necessary step to understand how planets form. Here we measured the mass accretion rate for young stellar objects at age >5 Myr, a critical test for the current models of disk evolution. We present the analysis of the spectra of 36 targets in the ~5-10 Myr old Upper Scorpius region for which disk masses were measured with ALMA. We find that the mass accretion rates in this sample of old but still survived disks are similarly high as those of the younger (<3 Myr old) star-forming regions of Lupus and Cha I, when considering the dependence on stellar and disk mass. In particular, several disks show high mass accretion rates >10^-9 Msun/yr while having low disk masses. Furthermore, the median values of the measured mass accretion rates in the disk mass ranges where our sample is complete at a level ~60-80% are compatible in these three regions. At the same time, the spread of mass accretion rates at any given disk mass is still >0.9 dex even at age>5 Myr. These results are in contrast with simple models of viscous evolution, which would predict that the values of the mass accretion rate diminish with time, and a tighter correlation with disk mass at age>5 Myr. Similarly, simple models of internal photoevaporation cannot reproduce the observed mass accretion rates, while external photoevaporation might explain the low disk masses and high accretion rates. A partial possible solution to the discrepancy with the viscous models is that the gas-to-dust ratio of the disks at >5 Myr is significantly different and higher than the canonical 100, as suggested by some dust and gas disk evolution models. The results shown here require the presence of several inter-playing processes, such as detailed dust evolution, external photoevaporation and possibly MHD winds, to explain the secular evolution of protoplanetary disks.

When imaged at high-resolution, many proto-planetary discs show gaps and rings in their dust sub-mm continuum emission profile. These structures are widely considered to originate from local maxima in the gas pressure profile. The properties of the underlying gas structures are however unknown. In this paper we present a method to measure the dust-gas coupling $\alpha/St$ and the width of the gas pressure bumps affecting the dust distribution, applying high-precision techniques to extract the gas rotation curve from emission lines data-cubes. As a proof-of-concept, we then apply the method to two discs with prominent sub-structure, HD163296 and AS 209. We find that in all cases the gas structures are larger than in the dust, confirming that the rings are pressure traps. Although the grains are sufficiently decoupled from the gas to be radially concentrated, we find that the degree of coupling of the dust is relatively good ($\alpha/St \sim 0.1$). We can therefore reject scenarios in which the disc turbulence is very low and the dust has grown significantly. If we further assume that the dust grain sizes are set by turbulent fragmentation, we find high values of the $\alpha$ turbulent parameter ($\alpha \sim 10^{-2}$). Alternatively, solutions with smaller turbulence are still compatible with our analysis if another process is limiting grain growth. For HD163296, recent measurements of the disc mass suggest that this is the case if the grain size is 1mm. Future constraints on the dust spectral indices will help to discriminate between the two alternatives.

Annular structures in proto-planetary discs, such as gaps and rings, are now ubiquitously found by high-resolution ALMA observations. Under the hypothesis that they are opened by planets, in this paper we investigate how the minimum planet mass needed to open a gap varies across different stellar host masses and distances from the star. The dependence on the stellar host mass is particularly interesting because, at least in principle, gap opening around low mass stars should be possible for lower mass planets, giving us a look into the young, low mass planet population. Using dusty hydrodynamical simulations, we find however the opposite behaviour, as a result of the fact that discs around low mass stars are geometrically thicker: gap opening around low mass stars can require more massive planets. Depending on the theoretical isochrone employed to predict the relationship between stellar mass and luminosity, the gap opening planet mass could also be independent of stellar mass, but in no case we find that gap opening becomes easier around low mass stars. This would lead to the expectation of a {\it lower} incidence of such structures in lower mass stars, since exoplanet surveys show that low mass stars have a lower fraction of giant planets. More generally, our study enables future imaging observations as a function of stellar mass to be interpreted using information on the mass vs. luminosity relations of the observed samples.

We predict magnitudes for young planets embedded in transition discs, still affected by extinction due to material in the disc. We focus on Jupiter-size planets at a late stage of their formation, when the planet has carved a deep gap in the gas and dust distributions and the disc starts being transparent to the planet flux in the infrared (IR). Column densities are estimated by means of three-dimensional hydrodynamical models, performed for several planet masses. Expected magnitudes are obtained by using typical extinction properties of the disc material and evolutionary models of giant planets. For the simulated cases located at $5.2$ AU in a disc with local unperturbed surface density of $127$ $\mathrm{g} \cdot \mathrm{cm}^{-2}$, a $1$ $M_J$ planet is highly extincted in J-, H- and K-bands, with predicted absolute magnitudes $\ge 50$ mag. In L- and M-bands extinction decreases, with planet magnitudes between $25$ and $35$ mag. In the N-band, due to the silicate feature on the dust opacities, the expected magnitude increases to $40$ mag. For a $2$ $M_J$ planet, the magnitudes in J-, H- and K-bands are above $22$ mag, while for L-, M- and N-bands the planet magnitudes are between $15$ and $20$ mag. For the $5$ $M_J$ planet, extinction does not play a role in any IR band, due to its ability to open deep gaps. Contrast curves are derived for the transition discs in CQ Tau, PDS70, HL Tau, TW Hya and HD163296. Planet mass upper-limits are estimated for the known gaps in the last two systems.

In this paper we analyse giant gap-opening planet migration in protoplanetary discs, focusing on the type II migration regime. According to standard type II theory, planets migrate at the same rate as the gas in the disc, as they are coupled to the disc viscous evolution; however, recent studies questioned this paradigm, suggesting that planets migrate faster than the disc material. We study the problem through 2D long-time simulations of systems consistent with type II regime, using the hydrodynamical grid code FARGO3D. Even though our simulations confirm the presence of an initial phase characterised by fast migration, they also reveal that the migration velocity slows down and eventually reaches the theoretical prediction if we allow the system to evolve for enough time. We find the same tendency to evolve towards the theoretical predictions at later times when we analyse the mass flow through the gap and the torques acting on the planet. This transient is related to the initial conditions of our (and previous) simulations, and is due to the fact that the shape of the gap has to adjust to a new profile, once the planet is set into motion. Secondly, we test whether the type II theory expectation that giant planet migration is driven by viscosity is consistent with our simulation by comparing simulations with the same viscosity and different disc mass (or viceversa). We find a good agreement with the theory, since when the discs are characterised by the same viscosity, the migration properties are the same.

Protoplanetary discs are the site of star and planet formation, and their evolution and consequent dispersal deeply affect the formation of planetary systems. In the standard scenario they evolve on timescales ~Myr due to the viscous transport of angular momentum. The analytical self-similar solution for their evolution predicts also specific disc isochrones in the accretion rate - disc mass plane. However, photoevaporation by radiation emitted by the central star is likely to dominate the gas disc dispersal of the innermost region, introducing another (shorter) timescale for this process. In this paper, we include the effect of internal (X and EUV) photoevaporation on the disc evolution, finding numerical solutions for a population of protoplanetary discs. Our models naturally reproduce the expected quick dispersal of the inner region of discs when their accretion rates match the rate of photoevaporative mass loss, in line with previous studies. We find that photoevaporation preferentially removes the lightest discs in the sample. The net result is that, counter-intuitively, photoevaporation increases the average disc mass in the sample, by dispersing the lightest discs. At the same time, photoevaporation also reduces the mass accretion rate by cutting the supply of material from the outer to the inner disc. In a purely viscous framework, this would be interpreted as the result of a longer viscous evolution, leading to an overestimate of the disc age. Our results thus show that photoevaporation is a necessary ingredient to include when interpreting observations of large disc samples with measured mass accretion rates and disc masses. Photoevaporation leaves a characteristic imprint on the shape of the isochrone. Accurate data in the accretion rate - disc mass plane in the low disc mass region therefore give clues on the typical photoevaporation rate.

Scattered light high-resolution imaging of the proto-planetary disc orbiting HD100453 shows two symmetric spiral arms, possibly launched by an external stellar companion. In this paper we present new, sensitive high-resolution ($\sim$30 mas) Band 7 ALMA observations of this source. This is the first source where we find counterparts in the sub-mm continuum to both scattered light spirals. The CO J=3-2 emission line also shows two spiral arms; in this case they can be traced over a more extended radial range, indicating that the southern spiral arm connects to the companion position. This is clear evidence that the companion is responsible for launching the spirals. The pitch angle of the sub-millimeter continuum spirals ($\sim 6 ^{\circ}$) is lower than the one in scattered light ($\sim 16 ^{\circ}$). We show that hydrodynamical simulations of binary-disc interaction can account for the difference in pitch angle only if one takes into account that the midplane is colder than the upper layers of the disc, as expected for the case of externally irradiated discs.

Binary stars often move through an ambient medium from which they accrete material and angular momentum, as in triple-star systems, star-forming clouds, young globular clusters and in the centres of galaxies. A binary form of Bondi-Hoyle-Lyttleton accretion results whereby the accretion rate depends on the binary properties: the stellar masses and separation, and the relative wind speed. We present the results of simulations performed with the hydrodynamic code GANDALF, to determine the mass accretion rates over a range of binary separations, inclinations and mass ratios. When the binary separation is short, the binary system accretes like a single star, while accretion onto stars in wide binaries is barely affected by their companion. We investigate intermediate-separation systems in some detail, finding that as the binary separation is increased, accretion rates smoothly decrease from the rate equal to that of a single star to the rate expected from two isolated stars. The form of this decrease depends on the relative centre-of-mass velocity of the binary and the gas, with faster-moving binaries showing a shallower decrease. Accretion rates vary little with orbital inclination, except when the orbit is side-on and the stars pass through each others' wakes. The specific angular momentum accretion rate also depends on the inclination but is never sufficient to prevent the binary orbit from contracting. Our results may be applied to accretion onto protostars, pollution of stars in globular and nuclear clusters, and wind mass-transfer in multiple stellar systems.

Protoplanetary discs (PPDs) in the Orion Nebula Cluster (ONC) are irradiated by UV fields from the massive star $\theta^1$C. This drives thermal winds, inducing mass loss rates of up to $\dot{M}_\mathrm{wind}\sim 10^{-7}\,M_\odot$/yr in the `proplyds' (ionised PPDs) close to the centre. For the mean age of the ONC and reasonable initial PPD masses, such mass loss rates imply that discs should have been dispersed. However, ~80% of stars still exhibit a NIR excess, suggesting that significant circumstellar mass remains. This `proplyd lifetime problem' has persisted since the discovery of photoevaporating discs in the core of the ONC by O'Dell & Wen (1994). In this work, we demonstrate how an extended period of star formation can solve this problem. Coupling N-body calculations and a viscous disc evolution model, we obtain high disc fractions at the present day. This is partly due to the migration of older stars outwards, and younger stars inwards such that the most strongly irradiated PPDs are also the youngest. We show how the disc mass distribution can be used to test the recent claims in the literature for multiple stellar populations in the ONC. Our model also explains the recent finding that host mass and PPD mass are only weakly correlated, in contrast with other regions of similar age. We conclude that the status of the ONC as the archetype for understanding the influence of environment on planet formation is undeserved; the complex star formation history (involving star formation episodes within ~0.8 Myr of the present day) results in confusing signatures in the PPD population.

The large majority of protoplanetary disks have very compact ($\lesssim15\,$AU) continuum emission at mm wavelengths. However, high angular resolution observations that resolve these small disks are still lacking, due to their intrinsically fainter emission compared with large bright disks. In this letter, we present $1.3\,$mm ALMA data of the faint ($\sim10\,$mJy) disk orbiting the TTauri star CX Tau at a resolution of $\sim40\,$mas, $\sim5\,$AU in diameter. The mm-dust disk is compact, with a 68$\%$ enclosing flux radius of 14$\,$AU, and the intensity profile exhibits a sharp drop between 10-20$\,$AU, and a shallow tail between 20-40$\,$AU. No clear signatures of substructure in the dust continuum are observed, down to the same sensitivity level of the DSHARP large program. However, the angular resolution does not allow to detect substructures at a scale of the disk aspect ratio in the inner regions. The radial intensity profile resembles well the inner regions of more extended disks imaged at the same resolution in DSHARP, but with no rings present in the outer disk. No inner cavity is detected, even though the disk has been classified as a transition disk from the spectral energy distribution in the near infrared. The emission of $^{12}$CO is much more extended, with a 68$\%$ enclosing flux radius of 75$\,$AU. The large difference of the mm dust and gas extents ($>5$) strongly points to the occurrence of radial drift, and matches well the predictions of theoretical models.

The combination of high resolution and sensitivity offered by ALMA is revolutionizing our understanding of protoplanetary discs, as their bulk gas and dust distributions can be studied independently. In this paper we present resolved ALMA observations of the continuum emission ($\lambda=1.3$ mm) and CO isotopologues ($^{12}$CO, $^{13}$CO, C$^{18}$O $J=2-1$) integrated intensity from the disc around the nearby ($d = 162$ pc), intermediate mass ($M_{\star}=1.67\,M_{\odot}$) pre-main-sequence star CQ Tau. The data show an inner depression in continuum, and in both $^{13}$CO and C$^{18}$O emission. We employ a thermo-chemical model of the disc reproducing both continuum and gas radial intensity profiles, together with the disc SED. The models show that a gas inner cavity with size between 15 and 25 au is needed to reproduce the data with a density depletion factor between $\sim 10^{-1}$ and $\sim 10^{-3}$. The radial profile of the distinct cavity in the dust continuum is described by a Gaussian ring centered at $R_{\rm dust}=53\,$au and with a width of $\sigma=13\,$au. Three dimensional gas and dust numerical simulations of a disc with an embedded planet at a separation from the central star of $\sim20\,$au and with a mass of $\sim 6\textrm{-} 9\,M_{\rm Jup}$ reproduce qualitatively the gas and dust profiles of the CQ Tau disc. However, a one planet model appears not to be able to reproduce the dust Gaussian density profile predicted using the thermo-chemical modeling.

Proto-planetary disc surveys conducted with ALMA are measuring disc radii in multiple star forming regions. The disc radius is a fundamental quantity to diagnose whether discs undergo viscous spreading, discriminating between viscosity or angular momentum removal by winds as drivers of disc evolution. Observationally, however, the sub-mm continuum emission is dominated by the dust, which also drifts inwards, complicating the picture. In this paper we investigate, using theoretical models of dust grain growth and radial drift, how the radii of dusty viscous proto-planetary discs evolve with time. Despite the existence of a sharp outer edge in the dust distribution, we find that the radius enclosing most of the dust $\textit{mass}$ increases with time, closely following the evolution of the gas radius. This behaviour arises because, although dust initially grows and drifts rapidly onto the star, the residual dust retained on Myr timescales is relatively well coupled to the gas. Observing the expansion of the dust disc requires using definitions based on high fractions of the disc $\textit{flux}$ (e.g. 95 per cent) and very long integrations with ALMA, because the dust grains in the outer part of the disc are small and have a low sub-mm opacity. We show that existing surveys lack the sensitivity to detect viscous spreading. The disc radii they measure do not trace the mass radius or the sharp outer edge in the dust distribution, but the outer limit of where the grains have significant sub-mm opacity. We predict that these observed radii should shrink with time.

A correlation between proto-planetary disc radii and sub-mm fluxes has been recently reported. In this Letter we show that the correlation is a sensitive probe of grain growth processes. Using models of grain growth and drift, we have shown in a companion paper that the observed disc radii trace where the dust grains are large enough to have a significant sub-mm opacity. We show that the observed correlation emerges naturally if the maximum grain size is set by radial drift, implying relatively low values of the viscous $\alpha$ parameter $ \lesssim 0.001$. In this case the relation has an almost universal normalisation, while if the grain size is set by fragmentation the flux at a given radius depends on the dust-to-gas ratio. We highlight two observational consequences of the fact that radial drift limits the grain size. The first is that the dust masses measured from the sub-mm could be overestimated by a factor of a few. The second is that the correlation should be present also at longer wavelengths (e.g. 3mm), with a normalisation factor that scales as the square of the observing frequency as in the optically thick case.

PDS 70b is the most robustly detected young planet imaged in the gap of a transition disk so far, found at a separation of ~195 mas (~22 au) from its host star and at a position angle of about 155 deg. This system is therefore a unique laboratory to characterize the properties of young planetary systems at the stage of their formation. We aim to trace direct and indirect imprints of PDS 70b on the gas and dust emission of the circumstellar disk in order to study the properties of this ~5 Myr young planetary system. We obtained ALMA band 7 observations of PDS 70 in dust continuum and $^{12}CO$ (3-2) and combined them with archival data resulting in an unprecedented angular resolution of about 70 mas (~8 au). We derive an upper limit on circumplanetary material at the location of PDS 70b of ~0.01 $M_{\oplus}$ and find a highly structured circumstellar disk in both dust and gas. The outer dust ring peaks at ~0.65" (74 au) and reveals a possible second unresolved peak at ~0.53" (60 au). The integrated intensity of CO also shows evidence of a depletion of emission at ~0.2" (23 au) with a width of ~0.1" (11 au). The gas kinematics show evidence of a deviation from Keplerian rotation inside $\lesssim$0.8" (91 au). This implies a pressure gradient that can account for the location of the dust ring well beyond the location of PDS 70b. Further in, we detect an inner disk which appears to be connected to the outer disk by a possible bridge feature in the North-West region in both gas and dust. We compare the observations to hydrodynamical simulations including a planet with different masses covering the estimated mass range previously derived from NIR photometry (~5-9 M$_{Jup}$). We find that even a planet with a mass of 10 M$_{Jup}$ may not be sufficient to explain the extent of the wide gap and an additional low-mass companion may be needed to account for the observed disk morphology.

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.

Many stars form in regions of enhanced stellar density, wherein the influence of stellar neighbours can have a strong influence on a protoplanetary disc (PPD) population. In particular, far ultraviolet (FUV) flux from massive stars drives thermal winds from the outer edge of PPDs, accelerating disc destruction. In this work, we present a novel technique for constraining the dynamical history of a star forming environment using PPD properties in a strongly FUV irradiated environment. Applying recent models for FUV induced mass loss rates to the PPD population of Cygnus OB2, we constrain how long ago primordial gas was expelled from the region; $ 0.5$ Myr ago if the Shakura & Sunyaev $\alpha$-viscosity parameter is $\alpha = 10^{-2}$ (corresponding to a viscous timescale of $\tau_\mathrm{visc} \approx 0.5$ Myr for a disc of scale radius $40$ au around a $1\, M_\odot$ star). This value of $\alpha$ is effectively an upper limit, since it assumes efficient extinction of FUV photons throughout the embedded phase. With this gas expulsion timescale we are able to produce a full dynamical model that fits kinematic and morphological data as well as disc fractions. We suggest Cygnus OB2 was originally composed of distinct massive clumps or filaments, each with a stellar mass $\sim 10^4 \, M_\odot$. Finally we predict that in regions of efficient FUV induced mass loss, disc mass $M_\mathrm{disc}$ as a function of stellar host mass $m_\mathrm{star}$ follows a power law with $M_\mathrm{disc} \propto m_\mathrm{star}^\beta$, where $\beta \gtrsim 2.7$ (depending on disc initial conditions and FUV exposure). This is steeper than observed correlations in regions of moderate FUV flux ($1 < \beta <1.9$), and offers a promising diagnostic to establish the influence of external photoevaporation in a given region.

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.

Extensive surveys of star-forming regions with Spitzer have revealed populations of disk-bearing young stellar objects. These have provided crucial constraints, such as the timescale of dispersal of protoplanetary disks, obtained by carefully combining infrared data with spectroscopic or X-ray data. While observations in various regions agree with the general trend of decreasing disk fraction with age, the Lupus V and VI regions appeared to have been at odds, having an extremely low disk fraction. Here we show, using the recent Gaia data release 2 (DR2), that these extremely low disk fractions are actually due to a very high contamination by background giants. Out of the 83 candidate young stellar objects (YSOs) in these clouds observed by Gaia, only five have distances of 150 pc, similar to YSOs in the other Lupus clouds, and have similar proper motions to other members in this star-forming complex. Of these five targets, four have optically thick (Class II) disks. On the one hand, this result resolves the conundrum of the puzzling low disk fraction in these clouds, while, on the other hand, it further clarifies the need to confirm the Spitzer selected diskless population with other tracers, especially in regions at low galactic latitude like Lupus V and VI. The use of Gaia astrometry is now an independent and reliable way to further assess the membership of candidate YSOs in these, and potentially other, star-forming regions.

Gas giants' early ($\lesssim 5$ Myr) orbital evolution occurs in a disc losing mass in part to photoevaporation driven by high energy irradiance from the host star. This process may ultimately overcome viscous accretion to disperse the disc and halt migrating giants by starving their orbits of gas, imprinting on giant planet separations in evolved systems. Inversion of this distribution could then give insight into whether stellar FUV, EUV or X-ray flux dominates photoevaporation, constraining planet formation and disc evolution models. We use a 1D hydrodynamic code in population syntheses for gas giants undergoing Type II migration in a viscously evolving disc subject to either a primarily FUV, EUV or X-ray flux from a pre-solar T Tauri star. The photoevaporative mass loss profile's unique peak location and width in each energetic regime produces characteristic features in the distribution of giant separations: a severe dearth of $\lesssim$ 2 M$_{\rm J}$ planets interior to 5 AU in the FUV scenario, a sharp concentration of $\lesssim$ 3 M$_{\rm J}$ planets between $\approx 1.5 - 2$ AU in the EUV case, and a relative abundance of $\approx 2 - 3.5$ M$_{\rm J}$ giants interior to 0.5 AU in the X-ray model. These features do not resemble the observational sample of gas giants with mass constraints, though our results do show some weaker qualitative similarities. We thus assess how the differing photoevaporative profiles interact with migrating giants and address the effects of large model uncertainties as a step to better connect disc models with trends in the exoplanet population.

Most stars form and spend their early life in regions of enhanced stellar density. Therefore the evolution of protoplanetary discs (PPDs) hosted by such stars are subject to the influence of other members of the cluster. Physically, PPDs might be truncated either by photoevaporation due to ultraviolet flux from massive stars, or tidal truncation due to close stellar encounters. Here we aim to compare the two effects in real cluster environments. In this vein we first review the properties of well studied stellar clusters with a focus on stellar number density, which largely dictates the degree of tidal truncation, and far ultraviolet (FUV) flux, which is indicative of the rate of external photoevaporation. We then review the theoretical PPD truncation radius due to an arbitrary encounter, additionally taking into account the role of eccentric encounters that play a role in hot clusters with a 1D velocity dispersion $\sigma_v > 2$ km/s. Our treatment is then applied statistically to varying local environments to establish a canonical threshold for the local stellar density ($n_{c} > 10^4$ pc$^{-3}$) for which encounters can play a significant role in shaping the distribution of PPD radii over a timescale $\sim 3$ Myr. By combining theoretical mass loss rates due to FUV flux with viscous spreading in a PPD we establish a similar threshold for which a massive disc is completely destroyed by external photoevaporation. Comparing these thresholds in local clusters we find that if either mechanism has a significant impact on the PPD population then photoevaporation is always the dominating influence.

The majority of stars form in a clustered environment. This has an impact on the evolution of surrounding protoplanetary discs (PPDs) due to either photoevaporation or tidal truncation. Consequently, the development of planets depends on formation environment. Here we present the first thorough investigation of tidally induced angular momentum loss in PPDs in the distant regime, partly motivated by claims in the literature for the importance of distant encounters in disc evolution. We employ both theoretical predictions and dynamical/hydrodynamical simulations in 2D and 3D. Our theoretical analysis is based on that of Ostriker (1994) and leads us to conclude that in the limit that the closest approach distance $x_{min} \gg r$, the radius of a particle ring, the fractional change in angular momentum scales as $(x_{min}/r)^{-5}$. This asymptotic limit ensures that the cumulative effect of distant encounters is minor in terms of its influence on disc evolution. The angular momentum transfer is dominated by the $m=2$ Lindblad resonance for closer encounters and by the $m=1$, $\omega = 0$ Lindblad resonance at large $x_{min}/r$. We contextualise these results by comparing expected angular momentum loss for the outer edge of a PPD due to distant and close encounters. Contrary to the suggestions of previous works we do not find that distant encounters contribute significantly to angular momentum loss in PPDs. We define an upper limit for closest approach distance where interactions are significant as a function of arbitrary host to perturber mass ratio $M_2/M_1$.

During the process of planet formation, the planet-discs interactions might excite (or damp) the orbital eccentricity of the planet. In this paper, we present two long ($t\sim 3\times 10^5$ orbits) numerical simulations: (a) one (with a relatively light disc, $M_{\rm d}/M_{\rm p}=0.2$) where the eccentricity initially stalls before growing at later times and (b) one (with a more massive disc, $M_{\rm d}/M_{\rm p}=0.65$) with fast growth and a late decrease of the eccentricity. We recover the well-known result that a more massive disc promotes a faster initial growth of the planet eccentricity. However, at late times the planet eccentricity decreases in the massive disc case, but increases in the light disc case. Both simulations show periodic eccentricity oscillations superimposed on a growing/decreasing trend and a rapid transition between fast and slow pericentre precession. The peculiar and contrasting evolution of the eccentricity of both planet and disc in the two simulations can be understood by invoking a simple toy model where the disc is treated as a second point-like gravitating body, subject to secular planet-planet interaction and eccentricity pumping/damping provided by the disc. We show how the counterintuitive result that the more massive simulation produces a lower planet eccentricity at late times can be understood in terms of the different ratios of the disc-to-planet angular momentum in the two simulations. In our interpretation, at late times the planet eccentricity can increase more in low-mass discs rather than in high-mass discs, contrary to previous claims in the literature.

Spiral arms have been observed in nearly a dozen protoplanetary discs in near-infrared scattered light and recently also in the sub-millimetre continuum. While one of the most compelling explanations is that they are driven by planetary or stellar companions, in all but one cases such companions have not yet been detected and there is even ambiguity on whether the planet should be located inside or outside the spirals. Here we use 3D hydrodynamic simulations to study the morphology of spiral density waves launched by embedded planets taking into account the vertical temperature gradient, a natural consequence of stellar irradiation. Our simulations show that the pitch angle of the spirals in thermally stratified discs is the lowest in the disc mid-plane and increases towards the disc surface. We combine the hydrodynamic simulations with 3D radiative transfer calculations to predict that the pitch-angle of planetary spirals observed in the near-infrared is higher than in the sub-millimetre. We also find that in both cases the spirals converge towards the planet. This provides a new powerful observational method to determine if the perturbing planet is inside our outside the spirals, as well as map the thermal stratification of the disc.

A large fraction of stars are in binary systems, yet the evolution of proto-planetary discs in binaries has been little explored from the theoretical side. In this paper we investigate the evolution of the discs surrounding the primary and secondary components of binary systems on the assumption that this is driven by photoevaporation induced by X-rays from the respective star. We show how for close enough separations (20-30 AU for average X-ray luminosities) the tidal torque of the companion changes the qualitative behaviour of disc dispersal from inside out to outside in. Fewer transition discs created by photoevaporation are thus expected in binaries. We also demonstrate that in close binaries the reduction in viscous time leads to accelerated disc clearing around both components, consistent with $\textit{unresolved}$ observations. When looking at the $\textit{differential}$ disc evolution around the two components, in close binaries discs around the secondary clear first due to the shorter viscous timescale associated with the smaller outer radius. In wide binaries instead the difference in photo-evaporation rate makes the secondaries longer lived, though this is somewhat dependent on the assumed scaling of viscosity with stellar mass. We find that our models are broadly compatible with the growing sample of $\textit{resolved}$ observations of discs in binaries. We also predict that binaries have higher accretion rates than single stars for the same disc mass. Thus binaries probably contribute to the observed scatter in the relationship between disc mass and accretion rate in young stars.

GANDALF is a new hydrodynamics and N-body dynamics code designed for investigating planet formation, star formation and star cluster problems. GANDALF is written in C++, parallelised with both OpenMP and MPI and contains a python library for analysis and visualisation. The code has been written with a fully object-oriented approach to easily allow user-defined implementations of physics modules or other algorithms. The code currently contains implementations of Smoothed Particle Hydrodynamics, Meshless Finite-Volume and collisional N-body schemes, but can easily be adapted to include additional particle schemes. We present in this paper the details of its implementation, results from the test suite, serial and parallel performance results and discuss the planned future development. The code is freely available as an open source project on the code-hosting website github at https://github.com/gandalfcode/gandalf and is available under the GPLv2 license.

The streaming instability is often invoked as solution to the fragmentation and drift barriers in planetesimal formation, catalyzing the aggregation of dust on kyr timescales to grow km-sized cores. However there remains a lack of consensus on the physical mechanism(s) responsible for initiating it. One potential avenue is disc photoevaporation, wherein the preferential removal of relatively dust-free gas increases the disc metallicity. Late in the disc lifetime, photoevaporation dominates viscous accretion, creating a gradient in the depleted gas surface density near the location of the gap. This induces a local pressure maximum that collects drifting dust particles, which may then become susceptible to the streaming instability. Using a one-dimensional viscous evolution model of a disc subject to internal X-ray photoevaporation, we explore the efficacy of this process to build planetestimals. Over a range of parameters we find that the amount of dust mass converted into planetesimals is often < 1 M_Earth and at most a few M_Earth spread across tens of AU. We conclude that photoevaporation may at best be relevant for the formation of debris discs, rather than a common mechanism for the formation of planetary cores. Our results are in contrast to a recent, similar investigation that considered an FUV-driven photoevaporation model and reported the formation of tens of M_Earth at large (> 100 AU) disc radii. The discrepancies are primarily a consequence of the different photoevaporation profiles assumed. Until observations more tightly constrain photoevaporation models, the relevance of this process to the formation of planets remains uncertain.

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.

We explore how measurements of protoplanetary disc masses and accretion rates provided by surveys of star forming regions can be analysed via the dimensionless accretion parameter, which we define as the product of the accretion rate and stellar age divided by the disc mass. By extending and generalising the study of Jones et al (2012), we demonstrate that this parameter should be less than or of order unity for a wide range of evolutionary scenarios, rising above unity only during the final stages of outside in clearing by external photoevaporation. We use this result to assess the reliability of disc mass estimates derived from CO isotopologues and submm continuum emission by examining the distribution of accretion efficiencies in regions which are not subject to external photoevaporation. We find that while dust based mass estimates produce results compatible with theoretical expectations assuming canonical dust/gas ratio, the systematically lower CO based estimates yield accretion efficiencies significantly above unity in contrast with the theory. This finding provides additional evidence that CO based disc masses are an under-estimate, in line with arguments that have been made on the basis of chemical modelling of relatively small samples. On the other hand, we demonstrate that dust based mass estimates are sufficiently accurate to reveal distinctly higher accretion efficiencies in the Trapezium cluster, where this result is expected given the evident importance of external photoevaporation. We therefore propose the dimensionless accretion parameter as a new diagnostic of external photoevaporation in other star forming regions.

RW Aur is a young binary system showing strong signatures of a recent tidal encounter between the circumprimary disk and the secondary star. The primary star has recently undergone two major dimming events ($\Delta$mag $\approx$ 2 in V-band), whose origin is still under debate. To shed light on the mechanism leading to the dimming events, we study the extinction properties, accretion variability, and gas kinematics using absorption lines from the material obscuring star RW Aur A. We compare our moderate resolution X-Shooter spectra of the dim state of RW Aur A with other spectral observations. In particular, we analyse archival high resolution UVES spectra obtained during the bright state of the system, in order to track the evolution of the spectral properties across the second dimming event. The spectrum obtained during the dim state shows narrow absorption lines in the Na and K optical doublets, where the former is saturated. With a velocity of -60 km/s these lines indicate that during the dim state the disk wind is either enhanced, or significantly displaced into the line of sight. The photometric evolution across the dimming event shows a gray extinction, and is correlated with a significant reduction of the EW of all photospheric lines. Emission lines tracing accretion do not vary significantly across the dimming. We conclude that the dimming event is related to a major perturbation on the inner disk. We suggest that the inner disk is occulting (most of) the star, and thus its photosphere, but is not occulting the accretion regions within a few stellar radii. Since observations of the outer disk indicate that the disk is modestly inclined (45 - 60 deg), we propose that the inner disk might be warped by a yet unseen (sub-)stellar companion, which may also explain the 2.77 day periodic variability of the spectral lines.

The dispersal of the circumstellar discs of dust and gas surrounding young low- mass stars has important implications for the formation of planetary systems. Photo- evaporation from energetic radiation from the central object is thought to drive the dispersal in the majority of discs, by creating a gap which disconnects the outer from the inner regions of the disc and then disperses the outer disc from the inside-out, while the inner disc keeps draining viscously onto the star. In this Letter we show that the disc around TW Hya, the closest protoplanetary disc to Earth, may be the first object where a photoevaporative gap has been imaged around the time at which it is being created. Indeed the detected gap in the ALMA images is consistent with the expectations of X-ray photoevaporation models, thus not requiring the presence of a planet. The photoevaporation model is also consistent with a broad range of properties of the TW Hya system, e.g. accretion rate and the location of the gap at the onset of dispersal. We show that the central, unresolved 870 {\mu}m continuum source might be produced by free free emission from the gas and/or residual dust inside the gap.

Following the recent discovery of the first radial velocity planet in a star still possessing a protoplanetary disc (CI Tau), we examine the origin of the planet's eccentricity (e $\sim 0.3$). We show through long timescale ($10^5$ orbits) simulations that the planetary eccentricity can be pumped by the disc, even when its local surface density is well below the threshold previously derived from short timescale integrations. We show that the disc may be able to excite the planet's orbital eccentricity in $<$ a Myr for the system parameters of CI Tau. We also perform two planet scattering experiments and show that alternatively the observed planet may plausibly have acquired its eccentricity through dynamical scattering of a migrating lower mass planet, which has either been ejected from the system or swallowed by the central star. In the latter case the present location and eccentricity of the observed planet can be recovered if it was previously stalled within the disc's magnetospheric cavity.

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.

A relation between the mass accretion rate onto the central young star and the mass of the surrounding protoplanetary disk has long been theoretically predicted and observationally sought. For the first time, we have accurately and homogeneously determined the photospheric parameters, mass accretion rate, and disk mass for an essentially complete sample of young stars with disks in the Lupus clouds. Our work combines the results of surveys conducted with VLT/X-Shooter and ALMA. With this dataset we are able to test a basic prediction of viscous accretion theory, the existence of a linear relation between the mass accretion rate onto the central star and the total disk mass. We find a correlation between the mass accretion rate and the disk dust mass, with a ratio that is roughly consistent with the expected viscous timescale when assuming an interstellar medium (ISM) gas-to-dust ratio. This confirms that mass accretion rates are related to the properties of the outer disk. We find no correlation between mass accretion rates and the disk mass measured by CO isotopologues emission lines, possibly owing to the small number of measured disk gas masses. This suggests that the mm-sized dust mass better traces the total disk mass and that masses derived from CO may be underestimated, at least in some cases.

We investigate the minimum planet mass that produces observable signatures in infrared scattered light and submm continuum images and demonstrate how these images can be used to measure planet masses to within a factor of about two. To this end we perform multi-fluid gas and dust simulations of discs containing low mass planets, generating simulated observations at $1.65 \mu$m, $10 \mu$m and $850 \mu$m. We show that the minimum planet mass that produces a detectable signature is $\sim 15 M_\oplus$: this value is strongly dependent on disc temperature and changes slightly with wavelength (favouring the submm). We also confirm previous results that there is a minimum planet mass of $\sim 20 M_\oplus$ that produces a pressure maximum in the disc: only planets above this threshold mass generate a dust trap that can eventually create a hole in the submm dust. Below this mass, planets produce annular enhancements in dust outward of the planet and a reduction in the vicinity of the planet. These features are in steady state and can be understood in terms of variations in the dust radial velocity, imposed by the perturbed gas pressure radial profile, analogous to a traffic jam. We also show how planet masses can be derived from structure in scattered light and sub-mm images. We emphasise that simulations with dust need to be run over thousands of planetary orbits so as to allow the gas profile to achieve a steady state and caution against the estimation of planet masses using gas only simulations.

Photo-evaporation and planet formation have both been proposed as mechanisms responsible for the creation of a transition disc. We have studied their combined effect through a suite of 2d simulations of protoplanetary discs undergoing X-ray photoevaporation with an embedded giant planet. In a previous work we explored how the formation of a giant planet triggers the dispersal of the inner disc by photo-evaporation at earlier times than what would have happened otherwise. This is particularly relevant for the observed transition discs with large holes and high mass accretion rates that cannot be explained by photo-evaporation alone. In this work we significantly expand the parameter space investigated by previous simulations. In addition, the updated model includes thermal sweeping, needed for studying the complete dispersal of the disc. After the removal of the inner disc the disc is a non accreting transition disc, an object that is rarely seen in observations. We assess the relative length of this phase, to understand if it is long lived enough to be found observationally. Depending on the parameters, especially on the X-ray luminosity of the star, we find that the fraction of time spent as a non-accretor greatly varies. We build a population synthesis model to compare with observations and find that in general thermal sweeping is not effective enough to destroy the outer disc, leaving many transition discs in a relatively long lived phase with a gas free hole, at odds with observations. We discuss the implications for transition disc evolution. In particular, we highlight the current lack of explanation for the missing non-accreting transition discs with large holes, which is a serious issue in the planet hypothesis.

Interactions between a supermassive black hole binary and the surrounding accretion disc can both assist the binary inspiral and align the black hole spins to the disc angular momentum. While binary migration is due to angular-momentum transfer within the circumbinary disc, the spin-alignment process is driven by the mass accreting on to each black hole. Mass transfer between different disc components thus couples the inspiral and the alignment process together. Mass is expected to leak through the cavity cleared by the binary, and preferentially accretes on to the lighter (secondary) black hole which orbits closer to the disc edge. Low accretion rate on to the heavier (primary) black hole slows the alignment process down. We revisit the problem and develop a semi-analytical model to describe the coupling between gas-driven inspiral and spin alignment, finding that binaries with mass ratio q<~0.2 approach the gravitational-wave driven inspiral in differential misalignment: light secondaries prevent primaries from aligning. Binary black holes with misaligned primaries are ideal candidates for precession effects in the strong-gravity regime and may suffer from moderately large (~1500 km/s) recoil velocities.

We investigate the influence of photoevaporation of protoplanetary discs on the final distribution of exoplanets semi-major axis distances. We model giant planet migration in viscous discs affected by photoevaporation driven by either pure EUV or soft X-ray radiation (XEUV). We show that the final exoplanet distributions are strongly dependant on the choice of the photoevaporation model. In particular, we find that XEUV is more efficient than pure EUV radiation at parking planets at approximately 1-2 AU distance from their central star, hence roughly reproducing the observed peak in the exoplanets semi-major axis distributions. We note however that a more quantitative comparison with the observations is hindered by the oversimplified treatment of planetary accretion, which severely affects migration rates. For this reason, caution should be used when using these models to constrain details of disc clearing and/or migration from the observations. Nevertheless our results indicate that disc dispersal by photoevaporation may be the main driver of the features in the exoplanets semi-major axis distribution observed by recent surveys.

Transitional disks (TDs) are thought to be a late evolutionary stage of protoplanetary disks with dust depleted inner regions. The mechanism responsible for this depletion is still under debate. To constrain the models it is mandatory to have a good understanding of the properties of the gas content of the inner disk. Using X-Shooter broad band -UV to NIR- medium resolution spectroscopy we derive the stellar, accretion, and wind properties of a sample of 22 TDs. The analysis of these properties allows us to put strong constraints on the gas content in a region very close to the star (<0.2 AU) which is not accessible with any other observational technique. We fit the spectra with a self-consistent procedure to derive simultaneously SpT,Av,and mass accretion rates (Macc) of the targets. From forbidden emission lines we derive the wind properties of the targets. Comparing our findings to values for cTTs, we find that Macc and wind properties of 80% of the TDs in our sample, which is strongly biased towards strongly accreting objects, are comparable to those of cTTs. Thus, there are (at least) some TDs with Macc compatible with those of cTTs, irrespective of the size of the dust inner hole.Only in 2 cases Macc are much lower, while the wind properties are similar. We do not see any strong trend of Macc with the size of the dust depleted cavity, nor with the presence of a dusty optically thick disk close to the star. In the TDs in our sample there is a gas rich inner disk with density similar to that of cTTs disks. At least for some TDs, the process responsible of the inner disk clearing should allow for a transfer of gas from the outer disk to the inner region. This should proceed at a rate that does not depend on the physical mechanism producing the gap seen in the dust emission and results in a gas density in the inner disk similar to that of unperturbed disks around stars of similar mass.

Young stars show evidence of accretion discs which evolve quickly and disperse with an e-folding time of $\sim$ 3Myr. This is in striking contrast with recent observations that suggest evidence for numerous $>30$ Myr old stars with an accretion disc in large star-forming complexes. We consider whether these observations of apparently old accretors could be explained by invoking Bondi-Hoyle accretion to rebuild a new disc around these stars during passage through a clumpy molecular cloud. We combine a simple Monte Carlo model to explore the capture of mass by such systems with a viscous evolution model to infer the levels of accretion that would be observed. We find that a significant fraction of stars may capture enough material via the Bondi-Hoyle mechanism to rebuild a disc of mass $\gtrsim$ 1 minimum-mass solar nebula, and $\lesssim 10\%$ accrete at observable levels at any given time. A significant fraction of the observed old accretors may be explained with our proposed mechanism. Such accretion may provide a chance for a second epoch of planet formation, and have unpredictable consequences for planetary evolution.

We present new 3 mm ATCA data of two Class I Young Stellar Objects in the Ophiucus star forming region: Elias29 and WL12. For our analysis we compare them with archival 1.1 mm SMA data. In the (u,v) plane the two sources present a similar behavior: a nearly constant non-zero emission at long baselines, which suggests the presence of an unresolved component and an increase of the fluxes at short baselines, related to the presence of an extended envelope. Our data analysis leads to unusually low values of the spectral index $\alpha_{\rm 1.1-3mm}$, which may indicate that mm-sized dust grains have already formed both in the envelopes and in the disk-like structures at such early stages. To explore the possible scenarios for the interpretation of the sources we perform a radiative transfer modeling using a Monte Carlo code, in order to take into account possible deviations from the Rayleigh-Jeans and optically thin regimes. Comparison between the model outputs and the observations indicates that dust grains may form aggregates up to millimeter size already in the inner regions of the envelopes of Class I YSOs. Moreover, we conclude that the embedded disk-like structures in our two Class Is are probably very compact, in particular in the case of WL12, with outer radii down to tens of AU.

Most stars form in a clustered environment. Therefore, it is important to assess how this environment influences the evolution of protoplanetary discs around young stars. In turn, this affects their ability to produce planets and ultimately life. We present here for the first time 3D SPH/N-body simulations that include both the hydrodynamical evolution of the discs around their natal stars, as well as the dynamics of the stars themselves. The discs are viscously evolving, accreting mass onto the central star and spreading. We find penetrating encounters to be very destructive for the discs as in previous studies, although the frequency of such encounters is low. We also find, however, that encounter influence the disc radii more strongly than other disc properties such as the disc mass. The disc sizes are set by the competition between viscous spreading and the disruptive effect of encounters. As discs spread, encounters become more and more important. In the regime of rapid spreading encounters simply truncate the discs, stripping the outer portions. In the opposite regime, we find that the effect of many distant encounters is able to limit the disc size. Finally, we predict from our simulations that disc sizes are limited by encounters at stellar densities exceeding $\sim 2-3 \times 10^3 \ \mathrm{pc}^{-2}$.

We analyse current measurements of accretion rates onto pre-main sequence stars as a function of stellar mass, and conclude that the steep dependance of accretion rates on stellar mass is real and not driven by selection/detection threshold, as has been previously feared. These conclusions are reached by means of statistical tests including a survival analysis which can account for upper limits. The power-law slope of the Mdot-M* relation is found to be in the range of 1.6-1.9 for young stars with masses lower than 1 Msun. The measured slopes and distributions can be easily reproduced by means of a simple disc model which includes viscous accretion and X-ray photoevaporation. We conclude that the Mdot-M* relation in pre-main sequence stars bears the signature of disc dispersal by X-ray photoevaporation, suggesting that the relation is a straight- forward consequence of disc physics rather than an imprint of initial conditions.

We assess the potential of planet formation instigating the early formation of a photoevaporation driven gap, up to radii larger than typical for photoevaporation alone. For our investigation we make use of hydrodynamics models of photoevaporating discs with a giant planet embedded. We find that, by reducing the mass accretion flow onto the star, discs that form giant planets will be dispersed at earlier times than discs without planets by X-ray photoevaporation. By clearing the portion of the disc inner of the planet orbital radius, planet formation induced photoevaporation (PIPE) is able to produce transition disc that for a given mass accretion rate have larger holes when compared to standard X-ray photoevaporation. This constitutes a possible route for the formation of the observed class of accreting transition discs with large holes, which are otherwise difficult to explain by planet formation or photoevaporation alone. Moreover, assuming that a planet is able to filter dust completely, PIPE produces a transition disc with a large hole and may provide a mechanism to quickly shut down accretion. This process appears to be too slow however to explain the observed desert in the population of transition disc with large holes and low mass accretion rates.

We investigate the evolution of the surface density of a circumbinary accretion disc after the mass loss induced by the merger of two supermassive black holes. We first introduce an analytical model, under the assumption of a disc composed of test particles, to derive the surface density evolution of the disc following the mass loss. The model predicts the formation of sharp density peaks in the disc; the model also allows us to compute the typical timescale for the formation of these peaks. To test and validate the model, we run numerical simulations of the process using the Smoothed Particle Hydrodynamics (SPH) code PHANTOM, taking fluid effects into account. We find good agreement in the shape and position of the peaks between the model and the simulations. In a fluid disc, however, the epicyclic oscillations induced by the mass loss can dissipate, and only some of the predicted peaks form in the simulation. To quantify how fast this dissipation proceeds, we introduce an appropriate parameter, and we show that it is effective in explaining the differences between the analytical, collisionless model and a real fluid disc.