The inner regions of the discs of high-mass young stellar objects (HMYSOs) are still poorly known due to the small angular scales and the high visual extinction involved. We deploy near-infrared (NIR) spectro-interferometry to probe the inner gaseous disc in HMYSOs and investigate the origin and physical characteristics of the CO bandhead emission (2.3-2.4 $\mu$m). We present the first GRAVITY/VLTI observations at high spectral (R=4000) and spatial (mas) resolution of the CO overtone transitions in NGC 2024 IRS2. The continuum emission is resolved in all baselines and is slightly asymmetric, displaying small closure phases ($\leq$8$^{\circ}$). Our best ellipsoid model provides a disc inclination of 34$^{\circ}$$\pm$1$^{\circ}$, a disc major axis position angle of 166$^{\circ}$$\pm$1$^{\circ}$, and a disc diameter of 3.99$\pm$0.09 mas (or 1.69$\pm$0.04 au, at a distance of 423 pc). The small closure phase signals in the continuum are modelled with a skewed rim, originating from a pure inclination effect. For the first time, our observations spatially and spectrally resolve the first four CO bandheads. Changes in visibility, as well as differential and closure phases across the bandheads are detected. Both the size and geometry of the CO-emitting region are determined by fitting a bidimensional Gaussian to the continuum-compensated CO bandhead visibilities. The CO-emitting region has a diameter of 2.74$\pm^{0.08}_{0.07}$ mas (1.16$\pm$0.03 au), and is located in the inner gaseous disc, well within the dusty rim, with inclination and $PA$ matching the dusty disc geometry, which indicates that both dusty and gaseous discs are coplanar. Physical and dynamical gas conditions are inferred by modelling the CO spectrum. Finally, we derive a direct measurement of the stellar mass of $M_*\sim$14.7$^{+2}_{-3.6}$ M$_{\odot}$ by combining our interferometric and spectral modelling results.

The formation and the evolution of protoplanetary disks are important stages in the lifetime of stars. The processes of disk evolution and planet formation are intrinsically linked. We spatially resolve with GRAVITY/VLTI in the K-band the sub au-scale region of 27 stars to gain statistical understanding of their properties. We look for correlations with stellar parameters, such as luminosity, mass, temperature and age. Our sample also cover a range of various properties in terms of reprocessed flux, flared or flat morphology, and gaps. We developed semi-physical geometrical models to fit our interferometric data. Our best models correspond to smooth and wide rings, implying that wedge-shaped rims at the dust sublimation edge are favored, as found in the H-band. The closure phases are generally non-null with a median value of ~10 deg, indicating spatial asymmetries of the intensity distributions. Multi-size grain populations could explain the closure phase ranges below 20-25 deg but other scenarios should be invoked to explain the largest ones. Our measurements extend the Radius-Luminosity relation to ~1e4 Lsun and confirm the significant spread around the mean relation observed in the H-band. Gapped sources exhibit a large N-to-K band size ratio and large values of this ratio are only observed for the members of our sample that would be older than 1 Ma, less massive, and with lower luminosity. In the 2 Ms mass range, we observe a correlation in the increase of the relative age with the transition from group II to group I, and an increase of the N-to-K size ratio. However, the size of the current sample does not yet permit us to invoke a clear universal evolution mechanism across the HAeBe mass range. The measured locations of the K-band emission suggest that these disks might be structured by forming young planets, rather than by depletion due to EUV, FUV, and X-ray photo-evaporation.

The Earth is strongly depleted in carbon compared to the dust in the ISM, implying efficient removal of refractory carbon before parent body formation. It has been argued that grains get rid of their carbon through oxidation and photolysis in the exposed upper disk layers. We assess the efficacy of these C-removal mechanisms accounting for the vertical and radial transport of grains. We obtain the carbon and carbon free mass budget of solids by solving two 1D advection-diffusion equations, accounting for the dust grain size distribution and radial transport. The carbon removal acts on the fraction of the grains that are in the exposed layer and requires efficient vertical transport. In models without radial transport, oxidation and photolysis can destroy most of the refractory carbon in terrestrial planet formation region. But it only reaches the observed depletion levels for extreme parameter combinations and requires that parent body formation was delayed by 1 Myr. Adding radial transport of solids prevents the depletion entirely, leaving refractory carbon equally distributed throughout the disk. It is unlikely that the observed carbon depletion can ultimately be attributed to mechanisms operating on small grains in the disk surface layers. Other mechanisms need to be studied, for example flash heating events or FU Ori outbursts in order to remove carbon quickly and deeply. However, a sustained drift barrier or strongly reduced radial grain mobility are necessary to prevent replenishment of carbon from the outer disk.

This paper investigates how the far-IR water ice features can be used to infer properties of disks around T Tauri stars and the water ice thermal history. We explore the power of future observations with SOFIA/HIRMES and SPICA's proposed far-IR instrument SAFARI. A series of detailed radiative transfer disk models around a representative T Tauri star are used to investigate how the far-IR water ice features at 45 and 63 micron change with key disk properties: disk size, grain sizes, disk dust mass, dust settling, and ice thickness. In addition, a series of models is devised to calculate the water ice emission features from warmup, direct deposit and cooldown scenarios of the water ice in disks. Photodesorption from icy grains in disk surfaces weakens the mid-IR water ice features by factors 4-5. The far-IR water ice emission features originate from small grains at the surface snow line in disks at distance of 10-100 au. Unless this reservoir is missing in disks (e.g. transitional disks with large cavities), the feature strength is not changing. Grains larger than 10 micron do not contribute to the features. Grain settling (using turbulent description) is affecting the strength of the ice features by at most 15%. The strength of the ice feature scales with the disk dust mass and water ice fraction on the grains, but saturates for dust masses larger than 1.e-4 Msun and for ice mantles that increase the dust mass by more than 50%. The various thermal histories of water ice leave an imprint on the shape of the features (crystalline/amorphous) as well as on the peak strength and position of the 45 micron feature. SOFIA/HIRMES can only detect crystalline ice features much stronger than simulated in our standard T Tauri disk model in deep exposures (1 hr). SPICA/SAFARI can detect the typical ice features in our standard T Tauri disk model in short exposures (10 min). (abbreviated)

Context. Multi-wavelength observations are indispensable in studying disk geometry and dust evolution processes in protoplanetary disks. Aims. We aimed to construct a 3-dimensional model of HD 163296 capable of reproducing simultaneously new observations of the disk surface in scattered light with the SPHERE instrument and thermal emission continuum observations of the disk midplane with ALMA. We want to determine why the SED of HD 163296 is intermediary between the otherwise well-separated group I and group II Herbig stars. Methods. The disk was modelled using the Monte Carlo radiative transfer code MCMax3D. The radial dust surface density profile was modelled after the ALMA observations, while the polarized scattered light observations were used to constrain the inclination of the inner disk component and turbulence and grain growth in the outer disk. Results. While three rings are observed in the disk midplane in millimeter thermal emission at $\sim$80, 124 and 200 AU, only the innermost of these is observed in polarized scattered light, indicating a lack of small dust grains on the surface of the outer disk. We provide two models capable of explaining this difference. The first model uses increased settling in the outer disk as a mechanism to bring the small dust grains on the surface of the disk closer to the midplane, and into the shadow cast by the first ring. The second model uses depletion of the smallest dust grains in the outer disk as a mechanism for decreasing the optical depth at optical and NIR wavelengths. In the region outside the fragmentation-dominated regime, such depletion is expected from state-of-the-art dust evolution models. We studied the effect of creating an artificial inner cavity in our models, and conclude that HD 163296 might be a precursor to typical group I sources.

To understand the chemical composition of planets, it is important to know the chemical composition of the region where they form in protoplanetary disks. Due to its fundamental role in chemical and biological processes, carbon is a key element to trace. We aim to identify the carriers and processes behind the extended NIR flux observed around several Herbig stars. We compare the extended NIR flux from objects in the PIONIER Herbig Ae/Be survey with their flux in the PAH features. HD 100453 is used as a benchmark case to investigate the influence of quantum heated particles, like PAHs or very small carbonaceous grains, in more detail. We use the Monte Carlo radiative transfer code MCMax to do a parameter study of the QHP size and scale- height and examine the influence of quantum heating on the amount of extended flux in the NIR visibilities. There is a correlation between the PAH feature flux of a disk and the amount of its extended NIR flux. We find that very small carbonaceous grains create the observed extended NIR flux around HD 100453 and still lead to a realistic SED. These results can not be achieved without using quantum heating effects, e.g. only with scattered light and grains in thermal equilibrium. It is possible to explain the extended NIR emission around Herbig stars with the presence of carbonaceous, quantum heated particles. Interferometric observations can be used to constrain the spatial distribution and typical size of carbonaceous material in the terrestrial planet forming region.

We combined hydrodynamical simulations of planet-disk interactions with dust evolution models that include coagulation and fragmentation of dust grains over a large range of radii and derived observational properties using radiative transfer calculations. We studied the role of the snow line in the survival of the inner disk of transition disks. Inside the snow line, the lack of ice mantles in dust particles decreases the sticking efficiency between grains. As a consequence, particles fragment at lower collision velocities than in regions beyond the snow line. This effect allows small particles to be maintained for up to a few Myrs within the first astronomical unit. These particles are closely coupled to the gas and do not drift significantly with respect to the gas. For lower mass planets (1$M_{\rm{Jup}}$), the pre-transition appearance can be maintained even longer because dust still trickles through the gap created by the planet, moves invisibly and quickly in the form of relatively large grains through the gap, and becomes visible again as it fragments and gets slowed down inside of the snow line. The global study of dust evolution of a disk with an embedded planet, including the changes of the dust aerodynamics near the snow line, can explain the concentration of millimetre-sized particles in the outer disk and the survival of the dust in the inner disk if a large dust trap is present in the outer disk. This behaviour solves the conundrum of the combination of both near-infrared excess and ring-like millimetre emission observed in several transition disks.