Unraveling the 3D physical structure, the temperature and density distribution, of protoplanetary discs is an essential step if we are to confront simulations of embedded planets or dynamical instabilities. In this paper we focus on Submillimeter Array observations of the edge-on source, Gomez's Hamburger, believed to host an over-density hypothesised to be a product of gravitational instability in the disc, GoHam b. We demonstrate that, by leveraging the well characterised rotation of a Keplerian disc to deproject observations of molecular lines in position-position-velocity space into disc-centric coordinates, we are able to map out the emission distribution in the (r, z) plane and (x, y, z) space. We show that 12CO traces an elevated layer of $z\,/\,r \sim 0.3$, while 13CO traces deeper in the disc at $z\,/\,r \lesssim 0.2$. We localize emission associated with GoHam b, finding it at deprojected radius of approximately 500 au and at polar angle of +\- 30 degrees from the disc major axis. At the spatial resolution of $\sim 1.5^{\prime\prime}$, GoHam b is spatially unresolved, with an upper limit to its radius of $<190$~au.

We present the first detection of the 13C17O J=3-2 transition toward the HL Tau protoplanetary disc. We find significantly more gas mass (at least a factor of ten higher) than has been previously reported using C18O emission. This brings the observed total disc mass to 0.2 M, which we consider to be a conservative lower limit. Our analysis of the Toomre Q profile suggests that this brings the disc into the regime of gravitational instability. The radial region of instability (50-110 au) coincides with the location of a proposed planet-carved gap in the dust disc and a spiral in the gas. We, therefore, propose that if the origin of the gap is confirmed to be due to a forming giant planet, then it is likely to have formed via the gravitational fragmentation of the protoplanetary disc.

Measurements of the gas mass are necessary to determine the planet formation potential of protoplanetary disks. Observations of rare CO isotopologues are typically used to determine disk gas masses; however, if the line emission is optically thick this will result in an underestimated disk mass. With ALMA we have detected the rarest stable CO isotopologue, 13C17O, in a protoplanetary disk for the first time. We compare our observations with the existing detections of 12CO, 13CO, C18O and C17O in the HD163296 disk. Radiative transfer modelling using a previously benchmarked model, and assuming interstellar isotopic abundances, significantly underestimates the integrated intensity of the 13C17O J=3-2 line. Reconciliation between the observations and the model requires a global increase in CO gas mass by a factor of 3.5. This is a factor of 2-6 larger than previous gas mass estimates using C18O. We find that C18O emission is optically thick within the CO snow line, while the 13C17O emission is optically thin and is thus a robust tracer of the bulk disk CO gas mass.

Observations of different molecular lines in protoplanetary disks provide valuable information on the gas kinematics, as well as constraints on the radial density and temperature structure of the gas. With ALMA we have detected H13CO+ (J=4-3) and HC15N (J=4-3) in the HD97048 protoplanetary disk for the first time. We compare these new detections to the ringed continuum mm-dust emission and the spatially resolved CO (J=3-2) and HCO+ (J=4-3) emission. The radial distributions of the H13CO+ and HC15N emission show hints of ringed sub-structure whereas, the optically thick tracers, CO and HCO+, do not. We calculate the HCO+/H13CO+ intensity ratio across the disk and find that it is radially constant (within our uncertainties). We use a physio-chemical parametric disk structure of the HD97048 disk with an analytical prescription for the HCO+ abundance distribution to generate synthetic observations of the HCO+ and H13CO+ disk emission assuming LTE. The best by-eye fit models require radial variations in the HCO+/H13CO+ abundance ratio and an overall enhancement in H13CO+ relative to HCO+. This highlights the need to consider isotope selective chemistry and in particular low temperature carbon isotope exchange reactions. This also points to the presence of a reservoir of cold molecular gas in the outer disk (T < 10K, R > 200au). Chemical models are required to confirm that isotope-selective chemistry alone can explain the observations presented here. With these data, we cannot rule out that the known dust substructure in the HD97048 disk is responsible for the observed trends in molecular line emission, and higher spatial resolution observations are required to fully explore the potential of optically thin tracers to probe planet-carved dust gaps. We also report non-detections of H13CO+ and HC15N in the HD100546 protoplanetary disk.

The composition of gas and solids in protoplanetary discs sets the composition of planets that form out of them. Recent chemical models have shown that the composition of gas and dust in discs evolves on Myr time-scales, with volatile species disappearing from the gas phase. However, discs evolve due to gas accretion and radial drift of dust on time-scales similar to these chemical time-scales. Here we present the first model coupling the chemical evolution in the disc mid-planes with the evolution of discs due to accretion and radial drift of dust. Our models show that transport will always overcome the depletion of CO$_2$ from the gas phase, and can also overcome the depletion of CO and CH$_4$ unless both transport is slow (viscous $\alpha \lesssim 10^{-3}$) and the ionization rate is high ($\zeta \approx 10^{-17}$). Including radial drift further enhances the abundances of volatile species because they are carried in on the surface of grains before evaporating left at their ice lines. Due to large differences in the abundances within 10 au for models with and without efficient radial drift, we argue that composition can be used to constrain models of planet formation via pebble accretion.

Modern astronomy has finally been able to observe protoplanetary disks in reasonable resolution and detail, unveiling the processes happening during planet formation. These observed processes are understood under the framework of disk-planet interaction, a process studied analytically and modeled numerically for over 40 years. Long a theoreticians' game, the wealth of observational data has been allowing for increasingly stringent tests of the theoretical models. Modeling efforts are crucial to support the interpretation of direct imaging analyses, not just for potential detections but also to put meaningful upper limits on mass accretion rates and other physical quantities in current and future large-scale surveys. This white paper addresses the questions of what efforts on the computational side are required in the next decade to advance our theoretical understanding, explain the observational data, and guide new observations. We identified the nature of accretion, ab initio planet formation, early evolution, and circumplanetary disks as major fields of interest in computational planet formation. We recommend that modelers relax the approximations of alpha-viscosity and isothermal equations of state, on the grounds that these models use flawed assumptions, even if they give good visual qualitative agreement with observations. We similarly recommend that population synthesis move away from 1D hydrodynamics. The computational resources to reach these goals should be developed during the next decade, through improvements in algorithms and the hardware for hybrid CPU/GPU clusters. Coupled with high angular resolution and great line sensitivity in ground based interferometers, ELTs and JWST, these advances in computational efforts should allow for large strides in the field in the next decade.

Consistent modeling of protoplanetary disks requires the simultaneous solution of both continuum and line radiative transfer, heating/cooling balance between dust and gas and, of course, chemistry. Such models depend on panchromatic observations that can provide a complete description of the physical and chemical properties and energy balance of protoplanetary systems. Along these lines we present a homogeneous, panchromatic collection of data on a sample of 85 T Tauri and Herbig Ae objects for which data cover a range from X-rays to centimeter wavelengths. Datasets consist of photometric measurements, spectra, along with results from the data analysis such as line fluxes from atomic and molecular transitions. Additional properties resulting from modeling of the sources such as disc mass and shape parameters. dust size and PAH properties are also provided for completeness. Targets were selected based on their properties data availability. Data from more than 50 different telescopes and facilities were retrieved and combined in homogeneous datasets directly from public data archives or after being extracted from more than 100 published articles. X-ray data for a subset of 56 sources represent an exception as they were reduced from scratch and are presented here for the first time. Compiled datasets along with a subset of continuum and emission-line models are stored in a dedicated database and distributed through a publicly accessible online system. All datasets contain metadata descriptors that allow to backtrack them to their original resources. The graphical user interface of the online system allows the user to visually inspect individual objects but also compare between datasets and models. It also offers to the user the possibility to download any of the stored data and metadata for further processing.

The European FP7 project DIANA has performed a coherent analysis of a large set of observations from protoplanetary disks by means of thermo-chemical disk models. The collected data include extinction-corrected stellar UV and X-ray input spectra (as seen by the disk), photometric fluxes, low and high resolution spectra, interferometric data, emission line fluxes, line velocity profiles and line maps. We define and apply a standardized modelling procedure to simultaneously fit all these data by state-of-the-art modelling codes (ProDiMo, MCFOST, MCMax) which solve the continuum and line radiative transfer, disk chemistry, and the heating & cooling balance for both the gas and the dust. We allow for up to two radial disk zones to obtain our best-fitting models that have about 20 free parameters. This approach is novel and unique in its completeness and level of consistency. In this paper, we present the results from pure SED fitting for 27 objects and from the all inclusive DIANA-standard models for 14 objects. We fit most infrared to millimeter emission line fluxes within a factor better than 3, simultaneously with SED, PAH features and radial brightness profiles extracted from images at various wavelengths. Our analysis shows a number of Herbig Ae and T Tauri stars with very cold and massive outer disks which are situated at least partly in the shadow of a tall and gas-rich inner disk. The disk masses derived are often in excess to previously published values, since these disks are partially optically thick even at millimeter wavelength and so cold that they emit less than in the Rayleigh-Jeans limit. Some line observations cannot be reproduced by the models, probably caused by foreground cloud absorption or object variability. Our data collection, the fitted physical disk parameters as well as the full model output are available at an online database (http://www.univie.ac.at/diana).

We present high resolution ($\sim$300 au) Atacama Large Millimeter/submillimeter Array (ALMA) observations of the massive young stellar object G11.92-0.61 MM 1. We resolve the immediate circumstellar environment of MM 1 in 1.3 mm continuum emission and CH$_{3}$CN emission for the first time. The object divides into two main sources - MM 1a, which is the source of a bipolar molecular outflow, and MM 1b, located 0.57'' (1920 au) to the South-East. The main component of MM 1a is an elongated continuum structure, perpendicular to the bipolar outflow, with a size of $0.141'' \times 0.050''$ ($480\times170$ au). The gas kinematics toward MM 1a probed via CH$_{3}$CN trace a variety of scales. The lower energy $J=12-11$ $K=3$ line traces extended, rotating gas within the outflow cavity, while the $v$8=1 line shows a clearly-resolved Keplerian rotation signature. Analysis of the gas kinematics and dust emission shows that the total enclosed mass in MM 1a is $40\pm5$ M$_{\odot}$ (where between 2.2-5.8 M$_{\odot}$ is attributed to the disk), while MM 1b is $<0.6$ M$_{\odot}$. The extreme mass ratio and orbital properties of MM 1a and MM 1b suggest that MM 1b is one of the first observed examples of the formation of a binary star via disk fragmentation around a massive young (proto)star.

Simulations of massive star formation predict the formation of discs with significant substructure, such as spiral arms and clumps due to fragmentation. Here we present a semi-analytic framework for producing synthetic observations of discs with substructure, in order to determine their observability in interferometric observations. Unlike post-processing of hydrodynamical models, the speed inherent to our approach permits us to explore a large parameter space of star and disc parameters, and thus constrain properties for real observations. We compute synthetic dust continuum and molecular line observations probing different disc masses, distances, inclinations, thermal structures, dust distributions, and number and orientation of spirals and fragments. With appropriate spatial and kinematic filtering applied, our models predict that ALMA observations of massive YSOs at <5 kpc distances should detect spirals in both gas and dust in strongly self-gravitating discs (i.e. discs with up to two spiral arms and strong kinematic perturbations). Detecting spirals will be possible in discs of arbitrary inclination, either by directly spatially resolving them for more face-on discs (inclinations up to ~50 degrees), or through a kinematic signature otherwise. Clumps resulting from disc fragmentation should be detectable in the continuum, if the clump is sufficiently hotter than the surrounding disc material.

At the time of formation, protoplanetary discs likely contain a comparable mass to their host protostars. As a result, gravitational instabilities (GIs) are expected to play a pivotal role in the early phases of disc evolution. However, as these young objects are heavily embedded, confirmation of GIs has remained elusive. Therefore, we use the radiative transfer code LIME to produce line images of a $0.17\,\mathrm{M}_{\odot}$ self-gravitating protosolar-like disc. We note the limitations of using LIME and explore methods to improve upon the default gridding. We synthesise noiseless observations to determine the sensitivities required to detect the spiral flux, and find that the line flux distribution does not necessarily correlate to the abundance density distribution; hence performing radiative transfer calculations is imperative. Moreover, the spiral features are seen in absorption, due to the GI-heated midplane and high extinction, which could be indicative of GI activity. If a small beamsize and appropriate molecular line are used then spatially resolving spirals in a protosolar-like disc should be possible with ALMA for an on-source time of 30 hr. Spectrally resolving non-axisymmetric structure takes only a tenth as long for a reasonable noise level, but attributing this structure to GI-induced activity would be tentative. Finally, we find that identifying finger-like features in PV diagrams of nearly edge-on discs, which are a direct indicator of spirals, is feasible with an on-source time of 19 hr, and hence likely offers the most promising means of confirming GI-driven spiral structure in young, embedded protoplanetary discs.

Recent high-sensitivity observations carried out with ALMA have revealed the presence of complex organic molecules (COMs) such as methyl cyanide (CH$_{\rm 3}$CN) and methanol (CH$_{\rm 3}$OH) in relatively evolved protoplanetary discs. The behaviour and abundance of COMs in earlier phases of disc evolution remains unclear. Here we combine a smoothed particle hydrodynamics simulation of a fragmenting, gravitationally unstable disc with a gas-grain chemical code. We use this to investigate the evolution of formamide (NH$_{\rm 2}$CHO), a pre-biotic species, in both the disc and in the fragments that form within it. Our results show that formamide remains frozen onto grains in the majority of the disc where the temperatures are $<$100 K, with a predicted solid-phase abundance that matches those observed in comets. Formamide is present in the gas-phase in three fragments as a result of the high temperatures ($\geq$200\,K), but remains in the solid-phase in one colder ($\leq$150 K) fragment. The timescale over which this occurs is comparable to the dust sedimentation timescales, suggesting that any rocky core which is formed would inherit their formamide content directly from the protosolar nebula.

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.

To date, there is no explanation as to why disc-tracing CO first overtone (or `bandhead') emission is not a ubiquitous feature in low- to medium-resolution spectra of massive young stellar objects, but instead is only detected toward approximately 25 per cent of their spectra. In this paper, we investigate the hypothesis that only certain mass accretion rates result in detectable bandhead emission in the near infrared spectra of MYSOs. Using an analytic disc model combined with an LTE model of the CO emission, we find that high accretion rates ($\gtrsim 10^{-4}\,{\rm M}_{\odot}{\mathrm{yr}}^{-1}$) result in large dust sublimation radii, a larger contribution to the $K$-band continuum from hot dust at the dust sublimation radius, and therefore correspondingly lower CO emission with respect to the continuum. On the other hand, low accretion rates ($\lesssim10^{-6}\,{\rm M}_{\odot}{\mathrm{yr}}^{-1}$) result in smaller dust sublimation radii, a correspondingly smaller emitting area of CO, and thus also lower CO emission with respect to the continuum. In general, moderate accretion rates produce the most prominent, and therefore detectable, CO first overtone emission. We compare our findings to a recent near-infrared spectroscopic survey of MYSOs, finding results consistent with our hypothesis. We conclude that the detection rate of CO bandhead emission in the spectra of MYSOs could be the result of MYSOs exhibiting a range of mass accretion rates, perhaps due to the variable accretion suggested by recent multi-epoch observations of these objects.

Protoplanetary discs around brown dwarfs and very low mass stars offer some of the best prospects for forming Earth-sized planets in their habitable zones. To this end, we study the nature of the disc around the very low mass star V410 X-ray 1, whose SED is indicative of an optically thick and very truncated dust disc, with our modelling suggesting an outer radius of only 0.6 au. We investigate two scenarios that could lead to such a truncation, and find that the observed SED is compatible with both. The first scenario involves the truncation of both the dust and gas in the disc, perhaps due to a previous dynamical interaction or the presence of an undetected companion. The second scenario involves the fact that a radial location of 0.6 au is close to the expected location of the H$_2$O snowline in the disc. As such, a combination of efficient dust growth, radial migration, and subsequent fragmentation within the snowline leads to an optically thick inner dust disc and larger, optically thin outer dust disc. We find that a firm measurement of the CO $J=2$--1 line flux would enable us to distinguish between these two scenarios, by enabling a measurement of the radial extent of gas in the disc. Many models we consider contain at least several Earth-masses of dust interior to 0.6 au, suggesting that V410 X-ray 1 could be a precursor to a system with tightly-packed inner planets, such as TRAPPIST-1.

Our understanding of protoplanetary disks has greatly improved over the last decade due to a wealth of data from new facilities. Unbiased dust surveys with Spitzer leave us with good constraints on the dust dispersal timescale of small grains in the terrestrial planet forming region. In the ALMA era, this can be confronted for the first time also with evolutionary timescales of mm grains in the outer disk. Gas surveys in the context of the existing multi-wavelength dust surveys will be a key in large statistical studies of disk gas evolution. Unbiased gas surveys are limited to ALMA CO submm surveys, where the quantitative interpretation is still debated. Herschel gas surveys have been largely biased, but [OI] 63 mic surveys and also accretion tracers agree qualitatively with the evolutionary timescale of small grains in the inner disk. Recent advances achieved by means of consistent multi-wavelength studies of gas AND dust in planet forming disks reveal the subtleties of the quantitative interpretation of gas surveys. Observational methods to determine disk masses e.g. from CO submm lines require the knowledge of the dust properties in the disk. Understanding not only the gas evolution, but also its chemical composition will provide crucial input for planet formation models. Kinetic chemical results give profoundly different answers than thermodynamic equilibrium in terms of the C/O ratios as well as the water ice/rock ratios. Again, dust has a key impact on the chemical evolution and composition of the gas. Grain growth for example affects freeze-out processes and strongly increases the cosmic ray induced UV field.

In this paper, we model the chemical evolution of a 0.25 M$_{\odot}$ protoplanetary disc surrounding a 1 M$_{\odot}$ star that undergoes fragmentation due to self-gravity. We use Smoothed Particle Hydrodynamics including a radiative transfer scheme, along with time-dependent chemical evolution code to follow the composition of the disc and resulting fragments over approximately 4000 yrs. Initially, four quasi-stable fragments are formed, of which two are eventually disrupted by tidal torques in the disc. From the results of our chemical modelling, we identify species that are abundant in the fragments (e.g. H$_{\rm 2}$O, H$_{\rm 2}$S, HNO, N$_{\rm 2}$, NH$_{\rm 3}$, OCS, SO), species that are abundant in the spiral shocks within the disc (e.g. CO, CH$_{\rm 4}$, CN, CS, H$_{\rm 2}$CO), and species which are abundant in the circumfragmentary material (e.g. HCO$^{\rm +}$). Our models suggest that in some fragments it is plausible for grains to sediment to the core before releasing their volatiles into the planetary envelope, leading to changes in, e.g., the C/O ratio of the gas and ice components. We would therefore predict that the atmospheric composition of planets generated by gravitational instability should not necessarily follow the bulk chemical composition of the local disc material.

Gravitational instabilities (GIs) are most likely a fundamental process during the early stages of protoplanetary disc formation. Recently, there have been detections of spiral features in young, embedded objects that appear consistent with GI-driven structure. It is crucial to perform hydrodynamic and radiative transfer simulations of gravitationally unstable discs in order to assess the validity of GIs in such objects, and constrain optimal targets for future observations. We utilise the radiative transfer code LIME to produce continuum emission maps of a $0.17\,\mathrm{M}_{\odot}$ self-gravitating protosolar-like disc. We note the limitations of using LIME as is and explore methods to improve upon the default gridding. We use CASA to produce synthetic observations of 270 continuum emission maps generated across different frequencies, inclinations and dust opacities. We find that the spiral structure of our protosolar-like disc model is distinguishable across the majority of our parameter space after 1 hour of observation, and is especially prominent at 230$\,$GHz due to the favourable combination of angular resolution and sensitivity. Disc mass derived from the observations is sensitive to the assumed dust opacities and temperatures, and therefore can be underestimated by a factor of at least 30 at 850$\,$GHz and 2.5 at 90$\,$GHz. As a result, this effect could retrospectively validate GIs in discs previously thought not massive enough to be gravitationally unstable, which could have a significant impact on the understanding of the formation and evolution of protoplanetary discs.

Chemical compositions of giant planets provide a means to constrain how and where they form. Traditionally, super-stellar elemental abundances in giant planets were thought to be possible due to accretion of metal-rich solids. Such enrichments are accompanied by oxygen-rich compositions (i.e. C/O below the disc's value, assumed to be solar, C/O=0.54). Without solid accretion the planets are expected to have sub-solar metallicity, but high C/O ratios. This arises because the solids are dominated by oxygen-rich species, e.g. H$_2$O and CO$_2$, which freeze out in the disk earlier than CO, leaving the gas metal poor but carbon-rich. Here we demonstrate that super-solar metallicities can be achieved by gas accretion alone when growth and radial drift of pebbles are considered in protoplanetary discs. Through this mechanism planets may simultaneously acquire super-solar metallicities and super-solar C/O ratios. This happens because the pebbles transport volatile species inward as they migrate through the disc, enriching the gas at snow lines where the volatiles sublimate. Furthermore, the planet's composition can be used to constrain where it formed. Since high C/H and C/O ratios cannot be created by accreting solids, it may be possible to distinguish between formation via pebble accretion and planetesimal accretion by the level of solid enrichment. Finally, we expect that Jupiter's C/O ratio should be near or above solar if its enhanced carbon abundance came through accreting metal rich gas. Thus Juno's measurement of Jupiter's C/O ratio should determine whether Jupiter accreted its metals from carbon rich gas or oxygen rich solids.

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 Herbig Ae/Be stars (HAeBes) allow an exploration of the properties of Pre-Main Sequence(PMS) stars above the low-mass range ($<2{\rm\thinspace M_{\odot}}$) and those bordering the high-mass range ($>8{\rm\thinspace M_{\odot}}$). This paper is the second in a series exploring accretion in 91 HAeBes with Very Large Telescope/X-shooter spectra. Equivalent width measurements are carried out on 32 different lines, spanning the UV to NIR, in order to obtain their line luminosities. The line luminosities were compared to accretion luminosities, which were determined directly from measurements of an UV-excess. When detected, emission lines always demonstrate a correlation with the accretion luminosity, regardless of detection frequency. The average relationship between accretion luminosity and line luminosity is found to be ${\thinspace L_{\rm acc}}\propto{\thinspace L_{\rm line}}^{1.16 \pm 0.15}$. This is in agreement with the findings in Classical T Tauri stars, although the HAeBe relationship is generally steeper, particularly towards the Herbig Be mass range. Since all observed lines display a correlation with the accretion luminosity, all of them can be used as accretion tracers. This has increased the number of accretion diagnostic lines in HAeBes tenfold. However, questions still remain on the physical origin of each line, which may not be due to accretion.

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 formation process of massive stars is not well understood, and advancement in our understanding benefits from high resolution observations and modelling of the gas and dust surrounding individual high-mass (proto)stars. Here we report sub-arcsecond (<1550 au) resolution observations of the young massive star G11.92-0.61 MM1 with the SMA and VLA. Our 1.3 mm SMA observations reveal consistent velocity gradients in compact molecular line emission from species such as CH$_3$CN, CH$_3$OH, OCS, HNCO, H$_2$CO, DCN and CH$_3$CH$_2$CN, oriented perpendicular to the previously reported bipolar molecular outflow from MM1. Modelling of the compact gas kinematics suggests a structure undergoing rotation around the peak of the dust continuum emission. The rotational profile can be well fit by a model of a Keplerian disc, including infall, surrounding an enclosed mass of 30-60M$_{\odot}$, of which 2-3M$_{\odot}$ is attributed to the disc. From modelling the CH$_3$CN emission, we determine that two temperature components, of 150 K and 230 K, are required to adequately reproduce the spectra. Our 0.9 and 3.0cm VLA continuum data exhibit an excess above the level expected from dust emission; the full centimetre-submillimetre wavelength spectral energy distribution of MM1 is well reproduced by a model including dust emission, an unresolved hypercompact H{\i}{\i} region, and a compact ionised jet. In combination, our results suggest that MM1 is an example of a massive proto-O star forming via disc accretion, in a similar way to that of lower mass stars.

There have been several recent detections of candidate Keplerian discs around massive young protostars. Given the relatively large disc-to-star mass ratios in these systems, and their young ages, it is worth investigating their propensity to becoming self-gravitating. To this end, we compute self-consistent, semi-analytic models of putative self-gravitating discs for five candidate disc systems. Our aim is not to fit exactly the observations, but to demonstrate that the expected dust continuum emission from marginally unstable self-gravitating discs can be quite weak, due to high optical depth at the midplane even at millimetre wavelengths. In the best cases, the models produce "observable" disc masses within a factor of <1.5 of those observed, with midplane dust temperatures comparable to measured temperatures from molecular line emission. We find in two cases that a self-gravitating disc model compares well with observations. If these discs are self-gravitating, they satisfy the conditions for disc fragmentation in their outer regions. These systems may hence have as-yet-unresolved low mass stellar companions, and are thus promising targets for future high angular resolution observations.

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 propose a set of standard assumptions for the modelling of Class II and III protoplanetary disks, which includes detailed continuum radiative transfer, thermo-chemical modelling of gas and ice, and line radiative transfer from optical to cm wavelengths. We propose new standard dust opacities for disk models, we present a simplified treatment of PAHs sufficient to reproduce the PAH emission features, and we suggest using a simple treatment of dust settling. We roughly adjust parameters to obtain a model that predicts typical Class II T Tauri star continuum and line observations. We systematically study the impact of each model parameter (disk mass, disk extension and shape, dust settling, dust size and opacity, gas/dust ratio, etc.) on all continuum and line observables, in particular on the SED, mm-slope, continuum visibilities, and emission lines including [OI] 63um, high-J CO lines, (sub-)mm CO isotopologue lines, and CO fundamental ro-vibrational lines. We find that evolved dust properties (large grains) often needed to fit the SED, have important consequences for disk chemistry and heating/cooling balance, leading to stronger emission lines in general. Strong dust settling and missing disk flaring have similar effects on continuum observations, but opposite effects on far-IR gas emission lines. PAH molecules can shield the gas from stellar UV radiation because of their strong absorption and negligible scattering opacities. The observable millimetre-slope of the SED can become significantly more gentle in the case of cold disk midplanes, which we find regularly in our T Tauri models. We propose to use line observations of robust chemical tracers of the gas, such as O, CO, and H2, as additional constraints to determine some key properties of the disks, such as disk shape and mass, opacities, and the dust/gas ratio, by simultaneously fitting continuum and line observations.

We present spatially and spectrally resolved Br-gamma emission around the planet-hosting, transitional Herbig Ae/Be star HD 100546. Aiming to gain insight into the physical origin of the line in possible relation to accretion processes, we carried out Br-gamma spectro-interferometry using AMBER/VLTI from three different baselines achieving spatial and spectral resolutions of 2-4 mas and 12000. The Br-gamma visibility is larger than that of the continuum for all baselines. Differential phases reveal a shift between the photocentre of the Br-gamma line -displaced 0.6 mas (0.06 au at 100 pc) NE from the star- and that of the K-band continuum emission -displaced 0.3 mas NE from the star. The photocentres of the redshifted and blueshifted components of the Br-gamma line are located NW and SE from the photocentre of the peak line emission, respectively. Moreover, the photocentre of the fastest velocity bins within the spectral line tends to be closer to that of the peak emission than the photocentre of the slowest velocity bins. Our results are consistent with a Br-gamma emitting region inside the dust inner rim (<0.25 au) and extending very close to the central star, with a Keplerian, disc-like structure rotating counter-clockwise, and most probably flared (25 deg). Even though the main contribution to the Br-gamma line does not come from gas magnetically channelled on to the star, accretion on to HD 100546 could be magnetospheric, implying a mass accretion rate of a few 10^-7 Msun/yr. This value indicates that the observed gas has to be replenished on time-scales of a few months to years, perhaps by planet-induced flows from the outer to the inner disc as has been reported for similar systems.

To date, most simulations of the chemistry in protoplanetary discs have used 1+1D or 2D axisymmetric $\alpha$-disc models to determine chemical compositions within young systems. This assumption is inappropriate for non-axisymmetric, gravitationally unstable discs, which may be a significant stage in early protoplanetary disc evolution. Using 3D radiative hydrodynamics, we have modelled the physical and chemical evolution of a 0.17 M$_{\odot}$ self-gravitating disc over a period of 2000 yr. The 0.8 M$_{\odot}$ central protostar is likely to evolve into a solar-like star, and hence this Class 0 or early Class I young stellar object may be analogous to our early Solar System. Shocks driven by gravitational instabilities enhance the desorption rates, which dominate the changes in gas-phase fractional abundances for most species. We find that at the end of the simulation, a number of species distinctly trace the spiral structure of our relatively low-mass disc, particularly CN. We compare our simulation to that of a more massive disc, and conclude that mass differences between gravitationally unstable discs may not have a strong impact on the chemical composition. We find that over the duration of our simulation, successive shock heating has a permanent effect on the abundances of HNO, CN and NH$_3$, which may have significant implications for both simulations and observations. We also find that HCO$^+$ may be a useful tracer of disc mass. We conclude that gravitational instabilities induced in lower mass discs can significantly, and permanently, affect the chemical evolution, and that observations with high-resolution instruments such as ALMA offer a promising means of characterising gravitational instabilities in protosolar discs.

Herbig Ae/Be stars span a key mass range that links low and high mass stars, and thus provide an ideal window from which to explore their formation. This paper presents VLT/X-Shooter spectra of 91 Herbig Ae/Be stars, HAeBes; the largest spectroscopic study of HAeBe accretion to date. A homogeneous approach to determining stellar parameters is undertaken for the majority of the sample. Measurements of the ultra-violet (UV) are modelled within the context of magnetospheric accretion, allowing a direct determination of mass accretion rates. Multiple correlations are observed across the sample between accretion and stellar properties: the youngest and often most massive stars are the strongest accretors, and there is an almost 1:1 relationship between the accretion luminosity and stellar luminosity. Despite these overall trends of increased accretion rates in HAeBes when compared to classical T Tauri stars, we also find noticeable differences in correlations when considering the Herbig Ae and Herbig Be subsets. This, combined with the difficulty in applying a magnetospheric accretion model to some of the Herbig Be stars, could suggest that another form of accretion may be occurring within the Herbig Be mass range.

Herbig Ae/Be stars lie in the mass range between low and high mass young stars, and therefore offer a unique opportunity to observe any changes in the formation processes that may occur across this boundary. This paper presents medium resolution VLT/X-Shooter spectra of six Herbig Ae/Be stars, drawn from a sample of 91 targets, and high resolution VLT/CRIRES spectra of five Herbig Ae/Be stars, chosen based on the presence of CO first overtone bandhead emission in their spectra. The X-Shooter survey reveals a low detection rate of CO first overtone emission (7 per cent), consisting of objects mainly of spectral type B. A positive correlation is found between the strength of the CO v=2-0 and Br {\gamma} emission lines, despite their intrinsic linewidths suggesting a separate kinematic origin. The high resolution CRIRES spectra are modelled, and are well fitted under the assumption that the emission originates from small scale Keplerian discs, interior to the dust sublimation radius, but outside the co-rotation radius of the central stars. In addition, our findings are in very good agreement for the one object where spatially resolved near-infrared interferometric studies have also been performed. These results suggest that the Herbig Ae/Be stars in question are in the process of gaining mass via disc accretion, and that modelling of high spectral resolution spectra is able to provide a reliable probe into the process of stellar accretion in young stars of intermediate to high masses.

We present medium resolution ($R\approx5300$) $K$-band integral field spectroscopy of six MYSOs. The targets are selected from the RMS survey, and we used the NIFS on the Gemini North telescope. The data show various spectral line features including Br$\gamma$, CO, H$_2$, and \mbox{He\,{\sc i}}. The Br$\gamma$ line is detected in emission in all objects with $v_\mathrm{FWHM}\sim100$ -- 200 kms$^{-1}$. V645 Cyg shows a high-velocity P-Cygni profile between -800 kms$^{-1}$ and -300 kms$^{-1}$. We performed three-dimensional spectroastrometry to diagnose the circumstellar environment in the vicinity of the central stars using the Br$\gamma$ line. We measured the centroids of the velocity components with sub-mas precision. The centroids allow us to discriminate the blueshifted and redshifted components in a roughly east--west direction in both IRAS 18151--1208 and S106 in Br$\gamma$. This lies almost perpendicular to observed larger scale outflows. We conclude, given the widths of the lines and the orientation of the spectroastrometric signature, that our results trace a disc wind in both IRAS 18151--1208 and S106. The CO $\nu=2-0$ absorption lines at low $J$ transitions are detected in IRAS 18151--1208 and AFGL 2136. We analysed the velocity structure of the neutral gas discs. In IRAS 18151--1208, the absorption centroids of the blueshifted and redshifted components are separated in a direction of north-east to south-west, nearly perpendicular to that of the larger scale $H_2$ jet. The position-velocity relations of these objects can be reproduced with central masses of 30 $M_{\sun}$ for IRAS 18151--1208 and 20 $M_{\sun}$ for AFGL 2136. We also detect CO $\nu=2-0$ bandhead emission in IRAS 18151--1208, S106 and V645 Cyg. The results can be fitted reasonably with a Keplerian rotation model, with masses of 15, 20 and 20 $M_{\sun}$ respectively.

The formation and earliest stages of protoplanetary discs remain poorly constrained by observations. ALMA will soon revolutionise this field. Therefore, it is important to provide predictions which will be valuable for the interpretation of future high sensitivity and high angular resolution observations. Here we present simulated ALMA observations based on radiative transfer modelling of a relatively massive (0.39 M_solar) self-gravitating disc embedded in a 10 M_solar dense core, with structure similar to the pre-stellar core L1544. We focus on simple species and conclude that C17O 3-2, HCO+ 3-2, OCS 26-25 and H2CO 404-303 lines can be used to probe the disc structure and kinematics at all scales.

We present a detailed spectroscopic investigation of the young eruptive star PV Cep, to improve our understanding of its nature and characterise its circumstellar environment after its last outburst in 2004. The analysis of our medium-resolution spectroscopy in the near-IR (0.9-2.35 um), collected in 2012 at the Large Binocular Telescope with the IR spectrograph LUCIFER, allows us to infer the main stellar parameters (visual extinction, accretion luminosity, mass accretion and ejection rates), and model the inner disc, jet, and wind. The NIR spectrum displays several strong emission lines associated with accretion/ejection activity and circumstellar environment. Our analysis shows that the brightness of PV Cep is fading, as well as the mass accretion rate (2x10^-7 Msun/yr^-1 in 2012 vs ~5x10^-6 Msun/yr^-1 in 2004), which is more than one order of magnitude lower than in the outburst phase. Among the several emission lines, only the [FeII] intensity increased after the outburst. The observed [FeII] emission delineates blue- and red-shifted lobes, both with high- and low-velocity components, which trace an asymmetric jet and wind, respectively. The observed emission in the jet has a dynamical age of ~8 years, indicating that it was produced during the last outburst. The mass ejection rate in both lobes is 1.5x10^-7 Msun/yr^-1, approximately matching the high accretion rate observed during and immediately after the outburst . The observed jet/outflow asymmetries are consistent with an inhomogeneous medium. Our modelling of the CO emission hints at a small-scale gaseous disc ring, extending from ~0.2-0.4 AU to ~3 AU from the source, with an inner temperature of ~3000 K. Our HI lines modelling indicates that most of the observed emission comes from an expanding disc wind at Te=10000 K. The line profiles are strongly affected by scattering, disc screening, and outflow self-absorption.

Massive stars play an important role in many areas of astrophysics, but numerous details regarding their formation remain unclear. In this paper we present and analyse high resolution (R ~ 30,000) near-infrared 2.3 micron spectra of 20 massive young stellar objects from the RMS database, in the largest such study of CO first overtone bandhead emission to date. We fit the emission under the assumption it originates from a circumstellar disc in Keplerian rotation. We explore three approaches to modelling the physical conditions within the disc - a disc heated mainly via irradiation from the central star, a disc heated mainly via viscosity, and a disc in which the temperature and density are described analytically. We find that the models described by heating mechanisms are inappropriate because they do not provide good fits to the CO emission spectra. We therefore restrict our analysis to the analytic model, and obtain good fits to all objects that possess sufficiently strong CO emission, suggesting circumstellar discs are the source of this emission. On average, the temperature and density structure of the discs correspond to geometrically thin discs, spread across a wide range of inclinations. Essentially all the discs are located within the dust sublimation radius, providing strong evidence that the CO emission originates close to the central protostar, on astronomical unit scales. In addition, we show that the objects in our sample appear no different to the general population of MYSOs in the RMS database, based on their near- and mid-infrared colours. The combination of observations of a large sample of MYSOs with CO bandhead emission and our detailed modelling provide compelling evidence of the presence of small scale gaseous discs around such objects, supporting the scenario in which massive stars form via disc accretion.

HD 327083 is a sgB[e] star that forms a binary system with an orbital semi-major axis of ~1.7 AU. Our previous observations using the VLTI and AMBER in the medium resolution K-band mode spatially resolved the environment of HD 327083. The continuum visibilities obtained indicate the presence of a circumbinary disc. CO bandhead emission was also observed. However, due to the limited spectral resolution of the previous observations, the kinematic structure of the emitting material was not constrained. In this paper, we address this and probe the source of the CO emission with high spectral resolution and spatial precision. We have observed HD 327083 with high spectral resolution (25 & 6 km/s) using AMBER and CRIRES. The observations are compared to kinematical models to constrain the source of the emission. It is shown that the CO bandhead emission can be reproduced using a model of a Keplerian disc with an inclination and size consistent with our previous VLTI observations. The model is compared to AMBER differential phase measurements, which have a precision as high as 30-micro-arcseconds. A differential phase signal corresponding to 0.15 milli-arcseconds (~5 sigma) is seen over the bandhead emission, which is in excellent agreement with the model that fits the CRIRES observations. In comparison, a model of an equatorial outflow, as envisaged in the standard sgB[e] scenario, does not reproduce the observations well. The excellent agreement between the disc model and observations in the spatial and spectral domains is compelling evidence that the CO bandhead emission of HD 327083 originates in a circumbinary Keplerian disc. In contrast, the model of an equatorial outflow cannot reproduce the observations well. This suggests that the standard sgB[e] scenario is not applicable to HD 327083, which supports the hypothesis that the B[e] behaviour of HD 327083 is due to binarity (ABRIDGED).

Until now, axisymmetric, alpha-disc models have been adopted for calculations of the chemical composition of protoplanetary discs. While this approach is reasonable for many discs, it is not appropriate when self-gravity is important. In this case, spiral waves and shocks cause temperature and density variations that affect the chemistry. We have adopted a dynamical model of a solar-mass star surrounded by a massive (0.39 Msun), self-gravitating disc, similar to those that may be found around Class 0 and early Class I protostars, in a study of disc chemistry. We find that for each of a number of species, e.g. H2O, adsorption and desorption dominate the changes in the gas-phase fractional abundance; because the desorption rates are very sensitive to temperature, maps of the emissions from such species should reveal the locations of shocks of varying strengths. The gas-phase fractional abundances of some other species, e.g. CS, are also affected by gas-phase reactions, particularly in warm shocked regions. We conclude that the dynamics of massive discs have a strong impact on how they appear when imaged in the emission lines of various molecular species.