Spiral arms in protoplanetary discs are thought to be linked to the presence of companions. We test the hypothesis that the double spiral arm morphology observed in the transition disc MWC 758 can be generated by an $\approx 10$ M$_{\rm Jup}$ companion on an eccentric orbit internal to the spiral arms. Previous studies on MWC 758 have assumed an external companion. We compare simulated observations from three dimensional hydrodynamics simulations of disc-companion interaction to scattered light, infrared and CO molecular line observations, taking into account observational biases. The inner companion hypothesis is found to explain the double spiral arms, as well as several additional features seen in MWC 758 -- the arc in the northwest, substructures inside the spiral arms, the cavity in CO isotopologues, and the twist in the kinematics. Testable predictions include detection of fainter spiral structure, detection of a point source south-southeast of the primary, and proper motion of the spiral arms.

We present Adaptive Optics (AO) near infrared (NIR) observations using VLT/NACO and Keck/NIRC2 of ODISEA targets. ODISEA is an ALMA survey of the entire population of circumstellar discs in the Ophiuchus molecular cloud. From the whole sample of ODISEA we select all the discs that are not already observed in the NIR with AO and that are observable with NACO or NIRC2. The NIR-ODISEA survey consists of 147 stars observed in NIR AO imaging for the first time, as well as revisiting almost all the binary systems of Ophiuchus present in the literature (20 out of 21). In total, we detect 20 new binary systems and one triple system. For each of them we calculate the projected separation and position angle of the companion, as well as their NIR and millimeter flux ratios. From the NIR contrast we derived the masses of the secondaries, finding that 9 of them are in the sub-stellar regime (30-50 \MJup). Discs in multiple systems reach a maximum total dust mass of $\sim$ 50 M$_{\oplus}$, while discs in single stars can reach a dust mass of 200 M$_{\oplus}$. Discs with masses above 10 M$_{\oplus}$ are found only around binaries with projected separations larger than $\sim$ 110 au. The maximum disc size is also larger around single star than binaries. However, since most discs in Ophiuchus are very small and low-mass, the effect of visual binaries is relatively weak in the general disc population.

Massive stars like company. However, low-mass companions have remained extremely difficult to detect at angular separations ($\rho$) smaller than 1" (~1000-3000 au considering typical distance to nearby massive stars) given the large brightness contrast between the companion and the central star. Constraints on the low-mass end of the companions mass-function for massive stars are however needed, for example to help distinguishing between various scenarios for the formation of massive stars. To obtain statistically significant constraint on the presence of low-mass companions beyond the typical detection limit of current surveys ($\Delta \mathrm{mag} \lesssim 5$ at $\rho \lesssim 1$"), we initiated a survey of O and Wolf-Rayet stars in the Carina region using the SPHERE coronagraphic instrument on the VLT. In this first paper, we aim to introduce the survey, to present the methodology and to demonstrate the capability of SPHERE for massive stars using the multiple system QZ~Car. High-contrast imaging techniques, such as angular- and spectral-differential imaging techniques as well as PSF-fitting, were applied to detect and measure the relative flux of companions in each spectral channel of the instrument. We detected 19 sources around the QZ~Car system with detection limits of 9~mag at $\rho > 200$~mas for IFS and as faint as 13~mag at $\rho$ > 1".8 for IRDIS (corresponding to sub-solar masses for potential companions). All but two are reported here for the first time. Based on this proof of concept, we showed that VLT/SPHERE allows us to reach the sub-solar mass regime of the companion mass function. This paves the way for this type of observation with a large sample of massive stars to provide novel constraints on the multiplicity of massive stars in a region of the parameter space that has remained inaccessible so far.

Recent mm-wavelength surveys performed with the Atacama Large Millimeter Array (ALMA) have revealed protoplanetary discs characterized by rings and gaps. A possible explanation for the origin of such rings is the tidal interaction with an unseen planetary companion. The protoplanetary disc around DS Tau shows a wide gap in the ALMA observation at 1.3 mm. We construct a hydrodynamical model for the dust continuum observed by ALMA assuming the observed gap is carved by a planet between one and five Jupiter masses. We fit the shape of the radial intensity profile along the disc major axis varying the planet mass, the dust disc mass, and the evolution time of the system. The best fitting model is obtained for a planet with $M_{\rm p}=3.5\,M_{\rm Jup}$ and a disc with $M_{\rm dust}= 9.6\cdot10^{-5}\,M_{\odot}$. Starting from this result, we also compute the expected signature of the planet in the gas kinematics, as traced by CO emission. We find that such a signature (in the form of a "kink" in the channel maps) could be observed by ALMA with a velocity resolution between $0.2-0.5\,\rm{kms}^{-1}$ and a beam size between 30 and 50 mas.

We report the discovery of a low-mass stellar companion around the young Herbig Be star MWC 297. We performed multi-epoch high-contrast imaging in the near infrared (NIR) with the Very Large Telescope (VLT)/Spectro-Polarimetric High-contrast Exoplanet REsearch (SPHERE) instrument. The companion is found at projected separation of 244.7$\pm$13.2 au and a position angle of 176.4$\pm$0.1 deg. The large separation supports formation via gravitational instability. From the spectrum, we estimate a mass of 0.1-0.5 M$_{\odot}$, the range conveying uncertainties in the extinction of the companion and in evolutionary models at young ages. The orbit coincides with a gap in the dust disk inferred from the Spectral Energy Distribution (SED). The young age ($\lesssim$ 1 Myr) and mass ratio with the central star ($\sim 0.01$) makes the companion comparable to PDS 70~b, suggesting a relation between formation scenarios and disk dynamics.

We present $L^\prime$-band Keck/NIRC2 imaging and $H$-band Subaru/AO188+HiCIAO polarimetric observations of CQ Tau disk with a new spiral arm. Apart from the spiral feature our observations could not detect any companion candidates. We traced the spiral feature from the $r^2$-scaled HiCIAO polarimetric intensity image and the fitted result is used for forward modeling to reproduce the ADI-reduced NIRC2 image. We estimated the original surface brightness after throughput correction in $L^\prime$-band to be $\sim126$ mJy/arcsec$^2$ at most. We suggest that the grain temperature of the spiral may be heated up to $\sim$200 K in order to explain both of the $H$- and $L^{\prime}$-bands results. The $H$-band emission at the location of the spiral originates from the scattering from the disk surface while both scattering and thermal emission may contribute to the $L^{\prime}$-band emission. If the central star is only the light source of scattered light, the spiral emission at $L^\prime$-band should be thermal emission. If an inner disk also acts as the light source, the scattered light and the thermal emission may equally contribute to the $L^\prime$-band spiral structure.

We present evidence for localised deviations from Keplerian rotation, i.e., velocity "kinks", in 8 of 18 circumstellar disks observed by the DSHARP program: DoAr 25, Elias 2-27, GW Lup, HD 143006, HD 163296, IM Lup, Sz 129 and WaOph 6. Most of the kinks are detected over a small range in both radial extent and velocity, suggesting a planetary origin, but for some of them foreground contamination prevents us from measuring their spatial and velocity extent. Because of the DSHARP limited spectral resolution and signal-to-noise in the 12CO J=2-1 line, as well as cloud contamination, the kinks are usually detected in only one spectral channel, and will require confirmation. The strongest circumstantial evidence for protoplanets in the absence of higher spectral resolution data and additional tracers is that, upon deprojection, we find that all of the candidate planets lie within a gap and/or at the end of a spiral detected in dust continuum emission. This suggests that a significant fraction of the dust gaps and spirals observed by ALMA in disks are caused by embedded protoplanets.

The mechanisms of planet formation are still under debate. We know little about how planets form, even if more than 4000 exoplanets have been detected to date. Recent investigations target the cot of newly born planets: the protoplanetary disk. At the first stages of their life, exoplanets still accrete material from the gas-rich disk in which they are embedded. Transitional disks are indeed disks that show peculiarities, such as gaps, spiral arms, and rings, which can be connected to the presence of substellar companions. To investigate what is responsible for these features, we selected all the known transitional disks in the solar neighborhood (<200 pc) that are visible from the southern hemisphere. We conducted a survey of 11 transitional disks (TDs) with the SPHERE instrument at the VLT. This is the largest Halpha survey that has been conducted so far to look for protoplanets. The observations were performed with the Halpha filter of ZIMPOL in order to target protoplanets that are still in the accretion stage. All the selected targets are very young stars, less than 20 Myr, and show low extinction in the visible. We reduced the ZIMPOL pupil stabilized data by applying the method of the angular spectral differential imaging (ASDI), which combines both techniques. The datacubes are composed of the CntHalpha and the narrow band filter Halpha, which are taken simultaneously to permit the suppression of the speckle pattern. The principal component analysis (PCA) method was employed for the reduction of the data. For each dataset, we derived the 5sigma contrast limit and converted it in upper limits on the accretion luminosity. We do not detect any new accreting substellar companions around the targeted transition disks down to an average contrast of 12 magnitudes at 0.2 arcsec from the central star (continues in the manuscript).

We present polarized light observations of the transitional disk around Sz 91 acquired with VLT/NaCo at $H$ (1.7$\mu$m) and $K_s$ (2.2$\mu$m) bands. We resolve the disk and detect polarized emission up to $\sim$0."5 ($\sim$80 au) along with a central cavity at both bands. We computed a radiative transfer model that accounts for the main characteristics of the polarized observations. We found that the emission is best explained by small, porous grains distributed in a disk with a $\sim$45 au cavity. Previous ALMA observations have revealed a large sub-mm cavity ($\sim$83 au) and extended gas emission from the innermost (<16 au) regions up to almost 400 au from the star. Dynamical clearing by multiple low-mass planets arises as the most probable mechanism for the origin of Sz 91's peculiar structure. Using new $L'$ band ADI observations we can rule out companions more massive than $M_p$ $\geq$ 8 $M_\mathrm{Jup}$ beyond 45 au assuming hot-start models. The disk is clearly asymmetric in polarized light along the minor axis, with the north side brighter than the south side. Differences in position angle between the disk observed at sub-mm wavelengths with ALMA and our NaCo observations were found. This suggests that the disk around Sz 91 could be highly structured. Higher signal-to-noise near-IR and sub-mm observations are needed to confirm the existence of such structures and to improve the current understanding in the origin of transitional disks.

Tidal encounters in star clusters perturb discs around young protostars. In Cuello et al. (2019a, Paper I) we detailed the dynamical signatures of a stellar flyby in both gas and dust. Flybys produce warped discs, spirals with evolving pitch angles, increasing accretion rates, and disc truncation. Here we present the corresponding observational signatures of these features in optical/near-infrared scattered light and (sub-) millimeter continuum and CO line emission. Using representative prograde and retrograde encounters for direct comparison, we post-process hydrodynamical simulations with radiative transfer methods to generate a catalogue of multi-wavelength observations. This provides a reference to identify flybys in recent near-infrared and sub-millimetre observations (e.g., RW Aur, AS 205, HV Tau & DO Tau, FU Ori, V2775 Ori, and Z CMa).

We test the hypothesis that the disc cavity in the `transition disc' Oph IRS 48 is carved by an unseen binary companion. We use 3D dust-gas smoothed-particle hydrodynamics simulations to demonstrate that marginally coupled dust grains concentrate in the gas over-density that forms in in the cavity around a low binary mass ratio binary. This produces high contrast ratio dust asymmetries at the cavity edge similar to those observed in the disc around IRS 48 and other transition discs. This structure was previously assumed to be a vortex. However, we show that the observed velocity map of IRS 48 displays a peculiar asymmetry that is not predicted by the vortex hypothesis. We show the unusual kinematics are naturally explained by the non-Keplerian flow of gas in an eccentric circumbinary cavity. We further show that perturbations observed in the isovelocity curves of IRS 48 may be explained as the product of the dynamical interaction between the companion and the disc. The presence of a $\sim$0.4 M$_{\odot}$ companion at a $\sim$10 au separation can qualitatively explain these observations. High spatial resolution line and continuum imaging should be able to confirm this hypothesis.

We still do not understand how planets form, or why extra-solar planetary systems are so different from our own solar system. But the last few years have dramatically changed our view of the discs of gas and dust around young stars. Observations with the Atacama Large Millimeter/submillimeter Array (ALMA) and extreme adaptive-optics systems have revealed that most --- if not all --- discs contain substructure, including rings and gaps, spirals, azimuthal dust concentrations, and shadows cast by misaligned inner discs. These features have been interpreted as signatures of newborn protoplanets, but the exact origin is unknown. Here we report the kinematic detection of a few Jupiter-mass planet located in a gas and dust gap at 130 au in the disc surrounding the young star HD 97048. An embedded planet can explain both the disturbed Keplerian flow of the gas, detected in CO lines, and the gap detected in the dust disc at the same radius. While gaps appear to be a common feature in protoplanetary discs, we present a direct correspondence between a planet and a dust gap, indicating that at least some gaps are the result of planet-disc interactions.

We present the first observational evidence for a circumplanetary disk around the protoplanet PDS~70~b, based on a new spectrum in the $K$ band acquired with VLT/SINFONI. We tested three hypotheses to explain the spectrum: Atmospheric emission from the planet with either (1) a single value of extinction or (2) variable extinction, and (3) a combined atmospheric and circumplanetary disk model. Goodness-of-fit indicators favour the third option, suggesting circumplanetary material contributing excess thermal emission --- most prominent at $\lambda \gtrsim 2.3 \mu$m. Inferred accretion rates ($\sim 10^{-7.8}$--$10^{-7.3} M_J$ yr$^{-1}$) are compatible with observational constraints based on the H$\alpha$ and Br$\gamma$ lines. For the planet, we derive an effective temperature of 1500--1600 K, surface gravity $\log(g)\sim 4.0$, radius $\sim 1.6 R_J$, mass $\sim 10 M_J$, and possible thick clouds. Models with variable extinction lead to slightly worse fits. However, the amplitude ($\Delta A_V \gtrsim 3$mag) and timescale of variation ($\lesssim$~years) required for the extinction would also suggest circumplanetary material.

Transition discs are prime targets to look for protoplanets and study planet-disc interactions. We present VLT/SINFONI observations of PDS~70, a transition disc with a recently claimed embedded protoplanet. We take advantage of the angular and spectral diversity present in our data for an optimal PSF modeling and subtraction using principal component analysis (PCA). We report the redetection of PDS 70 b, both the front and far side of the outer disc edge, and the detection of several extended features in the annular gap. We compare spectral differential imaging applied before (PCA-SADI), and after (PCA-ASDI) angular differential imaging. Our tests suggest that PCA-SADI better recovers extended features, while PCA-ASDI is more sensitive to point sources. We adapted the negative fake companion (NEGFC) technique to infer the astrometry of the companion, and derived $r = 193.5 \pm 4.9 \mathrm{mas}$ and PA = 158.7deg $\pm$ 3.0deg. We used both NEGFC and ANDROMEDA to infer the $K$-band spectro-photometry of the protoplanet, and found results consistent with recent VLT/SPHERE observations, except for their 2018/02 epoch measurement in the $K2$ filter. Finally, we derived an upper limit of $\dot{M_b} < 1.26 \times 10^{-7} \big[ \frac{5 M_{\rm Jup}}{M_b} \big] \big[ \frac{R_b}{R_{\rm Jup}}\big] M_{\rm Jup} $ yr$^{-1}$ for the accretion rate of the companion based on an adaptation of PCA-SADI/PCA-ASDI around the Br$\gamma$ line.

Resolved ALMA and VLA observations indicate the existence of two dust traps in the protoplanetary disc MWC 758. By means of 2D gas+dust hydrodynamical simulations post-processed with 3D dust radiative transfer calculations, we show that the spirals in scattered light, the eccentric, asymmetric ring and the crescent-shaped structure in the (sub)millimetre can all be caused by two giant planets: a 1.5-Jupiter mass planet at 35 au (inside the spirals) and a 5-Jupiter mass planet at 140 au (outside the spirals). The outer planet forms a dust-trapping vortex at the inner edge of its gap (at ~85 au), and the continuum emission of this dust trap reproduces the ALMA and VLA observations well. The outer planet triggers several spiral arms which are similar to those observed in polarised scattered light. The inner planet also forms a vortex at the outer edge of its gap (at ~50 au), but it decays faster than the vortex induced by the outer planet, as a result of the disc's turbulent viscosity. The vortex decay can explain the eccentric inner ring seen with ALMA as well as the low signal and larger azimuthal spread of this dust trap in VLA observations. Finding the thermal and kinematic signatures of both giant planets could verify the proposed scenario.

The large crescents imaged by ALMA in transition disks suggest that azimuthal dust trapping concentrates the larger grains, but centimetre-wavelengths continuum observations are required to map the distribution of the largest observable grains. A previous detection at ~1cm of an unresolved clump along the outer ring of MWC758 (Clump1), and buried inside more extended sub-mm continuum, motivates followup VLA observations. Deep multiconfiguration integrations reveal the morphology of Clump 1 and additional cm-wave components which we characterize via comparison with a deconvolution of recent 342GHz data (~1mm). Clump1, which concentrates ~1/3 of the whole disk flux density at ~1cm, is resolved as a narrow arc with a deprojected aspect ratio Chi>5.6, and with half the azimuthal width than at 342 GHz. The spectral trends in the morphology of Clump1 are quantitatively consistent with the Lyra-Lin prescriptions for dust trapping in an anticyclonic vortex, provided with porous grains (f~0.2+-0.2) in a very elongated (Chi~14+-3) and cold (T~23+-2K) vortex. The same prescriptions constrain the turbulence parameter alpha and the gas surface density Sigma_g through log10( alpha x Sigma_g /g/cm2)~-2.3+-0.4, thus requiring values for Sigma_g larger than a factor of a few compared to that reported in the literature from the CO isotopologues, if alpha <~ 1E-3. Such physical conditions imply an appreciably optically thick continuum even at cm-wavelengths (tau(33GHz)~0.2). A secondary and shallower peak at 342GHz is about twice fainter relative to Clump1 at 33GHz. Clump2 appears to be less efficient at trapping large grains.

We present the most sensitive direct imaging and radial velocity (RV) exploration of $\epsilon$ Eridani to date. $\epsilon$ Eridani is an adolescent planetary system, reminiscent of the early Solar system. It is surrounded by a prominent and complex debris disk which is likely stirred by one or several gas giant exoplanets. The discovery of the RV signature of a giant exoplanet was announced 15 years ago, but has met with scrutiny due to possible confusion with stellar noise. We confirm the planet with a new compilation and analysis of precise RV data spanning 30 years, and combine it with upper limits from our direct imaging search, the most sensitive ever performed. The deep images were taken in the Ms band (4.7$\mu$m) with the vortex coronagraph recently installed in W.M. Keck Observatory's infrared camera NIRC2, which opens a sensitive window for planet searches around nearby adolescent systems. The RV data and direct imaging upper limit maps were combined in an innovative joint Bayesian analysis, providing new constraints on the mass and orbital parameters of the elusive planet. $\epsilon$ Eridani b has a mass of $0.78^{+0.38}_{-0.12}$ $M_{Jup}$ and is orbiting $\epsilon$ Eridani at about $3.48\pm 0.02$ AU with a period of $7.37 \pm 0.07$ years. The eccentricity of $\epsilon$ Eridani b's orbit is $0.07^{+0.06}_{-0.05}$, an order of magnitude smaller than early estimates and consistent with a circular orbit. We discuss our findings from the standpoint of planet-disk interactions and prospects for future detection and characterization with the James Webb Space Telescope.

We introduce the Ophiuchus DIsc Survey Employing ALMA (ODISEA), a project aiming to study the entire population of Spitzer-selected protoplanetary discs in the Ophiuchus Molecular Cloud (~300 objects) from both millimeter continuum and CO isotopologues data. Here we present 1.3 mm/230 GHz continuum images of 147 targets at 0.2" (28 au) resolution and a typical rms of 0.15 mJy. We detect a total of 133 discs, including the individual components of 11 binary systems and 1 triple system. Fifty-three of these discs are spatially resolved. We find clear substructures (inner cavities, rings, gaps, and/or spiral arms) in 8 of the sources and hints of such structures in another 4 discs. We construct the disc luminosity function for our targets and perform comparisons to other regions. A simple conversion between flux and dust mass (adopting standard assumptions) indicates that all discs detected at 1.3 mm are massive enough to form one or more rocky planets. In contrast, only ~50 discs (~1/3 of the sample) have enough mass in the form of dust to form the canonical 10 M_Earth core needed to trigger runaway gas accretion and the formation of gas giant planets, although the total mass of solids already incorporated into bodies larger than cm scales is mostly unconstrained. The distribution in continuum disc sizes in our sample is heavily weighted towards compact discs: most detected discs have radii < 15 au, while only 23 discs (~15% of the targets) have radii > 30 au.

The circumstellar disk of the Herbig Fe star HD 142527 is host to several remarkable features including a warped inner disk, a 120 au-wide annular gap, a prominent dust trap and several spiral arms. A low-mass companion, HD 142527 B, was also found orbiting the primary star at $\sim$14 au. This study aims to better characterize this companion, which could help explain its impact on the peculiar geometry of the disk. We observed the source with VLT/SINFONI in $H$+$K$ band in pupil-tracking mode. Data were post-processed with several algorithms based on angular differential imaging (ADI). HD 142527 B is conspicuously re-detected in most spectral channels, which enables us to extract the first medium-resolution spectrum of a low-mass companion within 0.1'' from its central star. Fitting our spectrum with both template and synthetic spectra suggests that the companion is a young M2.5$\pm$1.0 star with an effective temperature of $3500\pm100$ K, possibly surrounded with a hot (1700 K) circum-secondary environment. Pre-main sequence evolutionary tracks provide a mass estimate of $0.34\pm0.06 M_{\odot}$, independent of the presence of a hot environment. However, the estimated stellar radius and age do depend on that assumption; we find a radius of $1.37 \pm 0.05 R_{\odot}$ (resp. $1.96 \pm 0.10 R_{\odot}$) and an age of $1.8^{+1.2}_{-0.5}$ Myr (resp. $0.75 \pm 0.25$ Myr) in the case of the presence (resp. absence) of a hot environment contributing in $H$+$K$. Our new values for the mass and radius yield a mass accretion rate of $\sim$5 $\times 10^{-9} M_{\odot}$ yr$^{-1}$ (2-3% that of the primary). Our results illustrate thus the potential for SINFONI+ADI to characterize faint close-in companions. The new spectral type makes HD 142527 B a twin of the well known TW Hya T-Tauri star, and the revision of its mass to higher values further supports its role in shaping the disk.

We present observations of the nearby (D$\sim$100\,pc) Herbig star HD~163296 taken with the vortex coronograph at Keck/NIRC2 in the L' band (3.7~$\mu$m), to search for planetary mass companions in the ringed disc surrounding this pre-main sequence star. The images reveal an arc-like region of scattered light from the disc surface layers that is likely associated with the first bright ring detected with ALMA in the $\lambda$=1.3mm dust continuum at $\sim$65~au. We also detect a point-like source at $\sim$0\farcs5 projected separation in the North-East direction, close to the inner edge of the second gap in the millimetre images. Comparing the point source photometry with the atmospheric emission models of non-accreting giant planets, we obtain a mass of 6--7~M$_J$ for a putative protoplanet, assuming a system age of 5~Myr. Based on the contrast at a 95\% level of completeness calculated on the emission-free regions of our images, we set upper limits for the masses of giant planets of 8--15~M$_J$, 4.5--6.5~M$_J$ and 2.5-4.0~M$_J$ at the locations of the first, second and third gap in the millimetre dust continuum, respectively. Further deep, high resolution thermal IR imaging of the HD~163296 system are warranted, to confirm the presence and nature of the point source and to better understand the structure of the dust disc.

We present 3D hydrodynamical models of the HD142527 protoplanetary disc, a bright and well studied disc that shows spirals and shadows in scattered light around a 100 au gas cavity, a large horseshoe dust structure in mm continuum emission, together with mysterious fast radial flows and streamers seen in gas kinematics. By considering several possible orbits consistent with the observed arc, we show that all of the main observational features can be explained by one mechanism - the interaction between the disc and the observed binary companion. We find that the spirals, shadows and horseshoe are only produced in the correct position angles by a companion on an inclined and eccentric orbit approaching periastron - the 'red' family from Lacour et al. (2016). Dust-gas simulations show radial and azimuthal concentration of dust around the cavity, consistent with the observed horseshoe. The success of this model in the HD142527 disc suggests other mm-bright transition discs showing cavities, spirals and dust asymmetries may also be explained by the interaction with central companions.

Transition disks offer the extraordinary opportunity to look for newly born planets and investigate the early stages of planet formation. In this context we observed the Herbig A5 star MWC 758 with the L band vector vortex coronagraph installed in the near-infrared camera and spectrograph NIRC2 at the Keck II telescope, with the aim of unveiling the nature of the spiral structure by constraining the presence of planetary companions in the system. Our high-contrast imaging observations show a bright (delta L=7.0+/-0.3 mag) point-like emission, south of MWC 758 at a deprojected separation of about 20 au (r=0.111+/- 0. 004 arcsec) from the central star. We also recover the two spiral arms (south-east and north-west), already imaged by previous studies in polarized light, and discover a third one to the south-west of the star. No additional companions were detected in the system down to 5 Jupiter masses beyond 0.6 arcsec from the star. We propose that the bright L band emission could be caused by the presence of an embedded and accreting protoplanet, although the possibility of it being an asymmetric disk feature cannot be excluded. The spiral structure is probably not related to the protoplanet candidate, unless on an inclined and eccentric orbit, and it could be due to one (or more) yet undetected planetary companions at the edge of or outside the spiral pattern. Future observations and additional simulations will be needed to shed light on the true nature of the point-like source and its link with the spiral arms.

We present ALMA 1.3 mm continuum observations at 0.2" (25 au) resolution of Elias 2-24, one of the largest and brightest protoplanetary disks in the Ophiuchus Molecular Cloud, and report the presence of three partially resolved concentric gaps located at ~20, 52, and 87 au from the star. We perform radiative transfer modeling of the disk to constrain its surface density and temperature radial profile and place the disk structure in the context of mechanisms capable of forming narrow gaps such as condensation fronts and dynamical clearing by actively forming planets. In particular, we estimate the disk temperature at the locations of the gaps to be 23, 15, and 12 K (at 20, 52, and 87 au respectively), very close to the expected snow-lines of CO (23-28 K) and N2 (12-15 K). Similarly, by assuming that the widths of the gaps correspond to 4-8 x the Hill radii of forming planets (as suggested by numerical simulations), we estimate planet masses in the range of 0.2-1.5 M_Jup, 1.0-8.0 M_Jup, and 0.02-0.15 M_Jup for the inner, middle, and outer gap, respectively. Given the surface density profile of the disk, the amount of "missing mass" at the location of each one of these gaps (between 4 and 20 M_Jup) is more than sufficient to account for the formation of such planets.

We describe a joint high contrast imaging survey for planets at Keck and VLT of the last large sample of debris disks identified by the Spitzer Space Telescope. No new substellar companions were discovered in our survey of 30 Spitzer-selected targets. We combine our observations with data from four published surveys to place constraints on the frequency of planets around 130 debris disk single stars, the largest sample to date. For a control sample, we assembled contrast curves from several published surveys targeting 277 stars which do not show infrared excesses. We assumed a double power law distribution in mass and semi-major axis of the form f(m,a) = $Cm^{\alpha}a^{\beta}$, where we adopted power law values and logarithmically flat values for the mass and semi-major axis of planets. We find that the frequency of giant planets with masses 5-20 $M_{\rm Jup}$ and separations 10-1000 AU around stars with debris disks is 6.27% (68% confidence interval 3.68 - 9.76%), compared to 0.73% (68% confidence interval 0.20 - 1.80%) for the control sample of stars without disks. These distributions differ at the 88% confidence level, tentatively suggesting distinctness of these samples.

The substellar companion HD206893b has recently been discovered by direct imaging of its disc-bearing host star with the SPHERE instrument. We investigate the atypical properties of the companion, which has the reddest near-infrared colours among all known substellar objects, either orbiting a star or isolated, and we provide a comprehensive characterisation of the host star-disc-companion system. We conducted a follow-up of the companion with adaptive optics imaging and spectro-imaging with SPHERE, and a multiinstrument follow-up of its host star. We obtain a R=30 spectrum from 0.95 to 1.64 micron of the companion and additional photometry at 2.11 and 2.25 micron. We carried out extensive atmosphere model fitting for the companions and the host star in order to derive their age, mass, and metallicity. We found no additional companion in the system in spite of exquisite observing conditions resulting in sensitivity to 6MJup (2MJup) at 0.5" for an age of 300 Myr (50 Myr). We detect orbital motion over more than one year and characterise the possible Keplerian orbits. We constrain the age of the system to a minimum of 50 Myr and a maximum of 700 Myr, and determine that the host-star metallicity is nearly solar. The comparison of the companion spectrum and photometry to model atmospheres indicates that the companion is an extremely dusty late L dwarf, with an intermediate gravity (log g 4.5-5.0) which is compatible with the independent age estimate of the system. Though our best fit corresponds to a brown dwarf of 15-30 MJup aged 100-300 Myr, our analysis is also compatible with a range of masses and ages going from a 50 Myr 12MJup planetary-mass object to a 50 MJup Hyades-age brown dwarf...

We present the Vortex Image Processing (VIP) library, a python package dedicated to astronomical high-contrast imaging. Our package relies on the extensive python stack of scientific libraries and aims to provide a flexible framework for high-contrast data and image processing. In this paper, we describe the capabilities of VIP related to processing image sequences acquired using the angular differential imaging (ADI) observing technique. VIP implements functionalities for building high-contrast data processing pipelines, encompass- ing pre- and post-processing algorithms, potential sources position and flux estimation, and sensitivity curves generation. Among the reference point-spread function subtraction techniques for ADI post-processing, VIP includes several flavors of principal component analysis (PCA) based algorithms, such as annular PCA and incremental PCA algorithm capable of processing big datacubes (of several gigabytes) on a computer with limited memory. Also, we present a novel ADI algorithm based on non-negative matrix factorization (NMF), which comes from the same family of low-rank matrix approximations as PCA and provides fairly similar results. We showcase the ADI capabilities of the VIP library using a deep sequence on HR8799 taken with the LBTI/LMIRCam and its recently commissioned L-band vortex coronagraph. Using VIP we investigated the presence of additional companions around HR8799 and did not find any significant additional point source beyond the four known planets. VIP is available at http://github.com/vortex-exoplanet/VIP and is accompanied with Jupyter notebook tutorials illustrating the main functionalities of the library.

HD 141569 A is a pre-main sequence B9.5 Ve star surrounded by a prominent and complex circumstellar disk, likely still in a transition stage from protoplanetary to debris disk phase. Here, we present a new image of the third inner disk component of HD 141569 A made in the L' band (3.8 micron) during the commissioning of the vector vortex coronagraph recently installed in the near-infrared imager and spectrograph NIRC2 behind the W.M. Keck Observatory Keck II adaptive optics system. We used reference point spread function subtraction, which reveals the innermost disk component from the inner working distance of $\simeq 23$ AU and up to $\simeq 70$ AU. The spatial scale of our detection roughly corresponds to the optical and near-infrared scattered light, thermal Q, N and 8.6 micron PAH emission reported earlier. We also see an outward progression in dust location from the L'-band to the H-band (VLT/SPHERE image) to the visible (HST/STIS image), likely indicative of dust blowout. The warm disk component is nested deep inside the two outer belts imaged by HST NICMOS in 1999 (respectively at 406 and 245 AU). We fit our new L'-band image and spectral energy distribution of HD 141569 A with the radiative transfer code MCFOST. Our best-fit models favor pure olivine grains, and are consistent with the composition of the outer belts. While our image shows a putative very-faint point-like clump or source embedded in the inner disk, we did not detect any true companion within the gap between the inner disk and the first outer ring, at a sensitivity of a few Jupiter masses.

Uncovering the ingredients and the architecture of planetary systems is a very active field of research that has fuelled many new theories on giant planet formation, migration, composition, and interaction with the circumstellar environment. We aim at discovering and studying new such systems, to further expand our knowledge of how low-mass companions form and evolve. We obtained high-contrast H-band images of the circumstellar environment of the F5V star HD206893, known to host a debris disc never detected in scattered light. These observations are part of the SPHERE High Angular Resolution Debris Disc Survey (SHARDDS) using the InfraRed Dual-band Imager and Spectrograph (IRDIS) installed on VLT/SPHERE. We report the detection of a source with a contrast of 3.6x10^{-5} in the H-band, orbiting at a projected separation of 270 milliarcsecond or 10 au, corresponding to a mass in the range 24 to 73 Mjup for an age of the system in the range 0.2 to 2 Gyr. The detection was confirmed ten months later with VLT/NaCo, ruling out a background object with no proper motion. A faint extended emission compatible with the disc scattered light signal is also observed. The detection of a low-mass companion inside a massive debris disc makes this system an analog of other young planetary systems such as beta Pictoris, HR8799 or HD95086 and requires now further characterisation of both components to understand their interactions.

For several years, we have been developing vortex phase masks based on sub-wavelength gratings, known as Annular Groove Phase Masks. Etched onto diamond substrates, these AGPMs are currently designed to be used in the thermal infrared (ranging from 3 to 13 {\mu}m). Our AGPMs were first installed on VLT/NACO and VLT/VISIR in 2012, followed by LBT/LMIRCam in 2013 and Keck/NIRC2 in 2015. In this paper, we review the development, commissioning, on-sky performance, and early scientific results of these new coronagraphic modes and report on the lessons learned. We conclude with perspectives for future developments and applications.

Transition disks correspond to a short stage between the young protoplanetary phase and older debris phase. Along this evolutionary sequence, the gas component disappears leaving room for a dust-dominated environment where already-formed planets signpost their gravitational perturbations. We endeavor to study the very inner region of the well-known and complex debris, but still gas-rich disk, around HD 141569A using the exquisite high-contrast capability of SPHERE at the VLT. Recent near-infrared (IR) images suggest a relatively depleted cavity within ~200 au, while former mid-IR data indicate the presence of dust at separations shorter than ~100 au. We obtained multi-wavelength images in the near-IR in J, H2, H3 and Ks bands with the IRDIS camera and a 0.95-1.35 micrometers spectral data cube with the IFS. Data were acquired in pupil-tracking mode, thus allowing for angular differential imaging. We discovered several new structures inside 1", of which the most prominent is a bright ring with sharp edges (semi-major axis: 0.4") featuring a strong north-south brightness asymmetry. Other faint structures are also detected from 0.4" to 1" in the form of concentric ringlets and at least one spiral arm. Finally, the VISIR data at 8.6 micrometers suggests the presence of an additional dust population closer in. Besides, we do not detect companions more massive than 1-3 mass of Jupiter. The performance of SPHERE allows us to resolve the extended dust component, which was previously detected at thermal and visible wavelengths, into very complex patterns with strong asymmetries ; the nature of these asymmetries remains to be understood. Scenarios involving shepherding by planets or dust-gas interactions will have to be tested against these observations.

Circumstellar asymmetries such as central warps have recently been shown to cast shadows on outer disks. We investigate the hydrodynamical consequences of such variable illumination on the outer regions of a transition disk, and the development of spiral arms. Using 2D simulations, we follow the evolution of a gaseous disk passively heated by the central star, under the periodic forcing of shadows with an opening angle of $\sim$28$^\circ$. With a lower pressure under the shadows, each crossing results in a variable azimuthal acceleration, which in time develops into spiral density waves. Their pitch angles evolve from $\Pi \sim 15^\circ-22^\circ$ at the onset, to $\sim$11$^\circ$-14$^\circ$, over $\sim$65~AU to 150~AU. Self-gravity enhances the density contrast of the spiral waves, as also reported previously for spirals launched by planets. Our control simulations with unshadowed irradiation do not develop structures, except for a different form of spiral waves seen at later times only in the gravitationally unstable control case. Scattered light predictions in the $H$-band show that such illumination spirals should be observable. We suggest that spiral arms in the case-study transition disk HD~142527 could be explained as a result of shadowing from the tilted inner disk.

The finding of residual gas in the large central cavity of the HD142527 disk motivates questions on the origin of its non-Keplerian kinematics, and possible connections with planet formation. We aim to understand the physical structure that underlies the intra-cavity gaseous flows, guided by new molecular-line data in CO(6-5) with unprecedented angular resolutions. Given the warped structure inferred from the identification of scattered-light shadows cast on the outer disk, the kinematics are consistent, to first order, with axisymmetric accretion onto the inner disk occurring at all azimuth. A steady-state accretion profile, fixed at the stellar accretion rate, explains the depth of the cavity as traced in CO isotopologues. The abrupt warp and evidence for near free-fall radial flows in HD 142527 resemble theoretical models for disk tearing, which could be driven by the reported low mass companion, whose orbit may be contained in the plane of the inner disk. The companion's high inclination with respect to the massive outer disk could drive Kozai oscillations over long time-scales; high-eccentricity periods may perhaps account for the large cavity. While shadowing by the tilted disk could imprint an azimuthal modulation in the molecular-line maps, further observations are required to ascertain the significance of azimuthal structure in the density field inside the cavity of HD142527.

A pathway to the formation of planetesimals, and eventually giant planets, may occur in concentrations of dust grains trapped in pressure maxima. Dramatic crescent-shaped dust concentrations have been seen in recent radio images at sub-mm wavelengths. These disk asymmetries could represent the initial phases of planet formation in the dust trap scenario, provided that grain sizes are spatially segregated. A testable prediction of azimuthal dust trapping is that progressively larger grains should be more sharply confined and furthermore the trapped grains should follow a distribution that is markedly different from the gas. However, gas tracers such as CO and the infrared emission from small grains are both very optically thick where the submm continuum originates, so observations have been unable to test the trapping predictions or to identify compact concentrations of larger grains required for planet formation by core-accretion. Here we report multifrequency observations of HD142527, from 34GHz to 700GHz, that reveal a compact concentration of ~cm-sized grains, with a few Earth masses, embedded in a large-scale crescent of ~mm-sized particles. The emission peaks at wavelengths shorter than ~1mm are optically thick and trace the temperature structure resulting from shadows cast by the inner regions. Given this temperature structure, we infer that the largest dust grains are concentrated in the 34 GHz clump. We conclude that dust trapping is efficient for approximately cm-sized grains and leads to enhanced concentrations, while the smaller grains largely reflect the gas distribution.

The study of dynamical processes in protoplanetary disks is essential to understand planet formation. In this context, transition disks are prime targets because they are at an advanced stage of disk clearing and may harbor direct signatures of disk evolution. In this paper, we aim to derive new constraints on the structure of the transition disk MWC 758, to detect non-axisymmetric features and understand their origin. We obtained infrared polarized intensity observations of the protoplanetary disk MWC 758 with SPHERE/VLT at 1.04 microns to resolve scattered light at a smaller inner working angle (0.093") and a higher angular resolution (0.027") than previously achieved. We observe polarized scattered light within 0.53" (148 au) down to the inner working angle (26 au) and detect distinct non-axisymmetric features but no fully depleted cavity. The two small-scale spiral features that were previously detected with HiCIAO are resolved more clearly, and new features are identified, including two that are located at previously inaccessible radii close to the star. We present a model based on the spiral density wave theory with two planetary companions in circular orbits. The best model requires a high disk aspect ratio (H/r~0.20 at the planet locations) to account for the large pitch angles which implies a very warm disk. Our observations reveal the complex morphology of the disk MWC758. To understand the origin of the detected features, the combination of high-resolution observations in the submillimeter with ALMA and detailed modeling is needed.

Inner cavities and annular gaps in circumstellar disks are possible signposts of giant planet formation. The young star HD 142527 hosts a massive protoplanetary disk with a large cavity that extends up to 140 au from the central star, as seen in continuum images at infrared and millimeter wavelengths. Estimates of the survival of gas inside disk cavities are needed to discriminate between clearing scenarios. We present a spatially and spectrally resolved carbon monoxide isotopologue observations of the gas-rich disk HD 142527, in the J=2-1 line of 12CO, 13CO and C18O, obtained with the Atacama Large Millimeter Array (ALMA). We detect emission coming from inside the dust-depleted cavity in all three isotopologues. Based on our analysis of the gas in the dust cavity, the 12CO emission is optically thick, while 13CO and C18O emission are both optically thin. The total mass of residual gas inside the cavity is about 1.5-2 Jupiter masses. We model the gas with an axisymmetric disk model. Our best fit model shows that the cavity radius is much smaller in CO than it is in millimeter continuum and scattered light observations, with a gas cavity that does not extend beyond 105 au (at 3-sigma). The gap wall at its outer edge is diffuse and smooth in the gas distribution, while in dust continuum it is manifestly sharper. The inclination angle, as estimated from the high velocity channel maps, is 28+/-0.5 degrees, higher than in previous estimates, assuming a fix central star mass of 2.2 Solar masses.

(abridged) Vortex coronagraphs are among the most promising solutions to perform high contrast imaging at small angular separations. They feature a very small inner working angle, a clear 360 degree discovery space, have demonstrated very high contrast capabilities, are easy to implement on high-contrast imaging instruments, and have already been extensively tested on the sky. Since 2005, we have been designing, developing and testing an implementation of the charge-2 vector vortex phase mask based on concentric subwavelength gratings, referred to as the Annular Groove Phase Mask (AGPM). Science-grade mid-infrared AGPMs were produced in 2012 for the first time, using plasma etching on synthetic diamond substrates. They have been validated on a coronagraphic test bench, showing broadband peak rejection up to 500:1 in the L band, which translates into a raw contrast of about $6\times 10^{-5}$ at $2 \lambda/D$. Three of them have now been installed on world-leading diffraction-limited infrared cameras (VLT/NACO, VLT/VISIR and LBT/LMIRCam). During the science verification observations with our L-band AGPM on NACO, we observed the beta Pictoris system and obtained unprecedented sensitivity limits to planetary companions down to the diffraction limit ($0.1''$). More recently, we obtained new images of the HR 8799 system at L band during the AGPM first light on LMIRCam. After reviewing these first results obtained with mid-infrared AGPMs, we will discuss the short- and mid-term goals of the on-going VORTEX project, which aims to improve the performance of our vortex phase masks for future applications on second-generation high-contrast imagers and on future extremely large telescopes (ELTs).

The disk around HD 142527 attracts a lot of attention, amongst others because of its resolved (sub) mm dust continuum that is concentrated into a horseshoe-shape towards the north of the star. In this manuscript we present spatially resolved ALMA detections of the HCN J=4-3 and CS J=7-6 emission lines. These lines give us a view deeper into the disk compared to the (optically thicker) CO isotopes. This is the first detection of CS J=7-6 coming from a protoplanetary disk. Both emission lines are azimuthally asymmetric and are suppressed under the horseshoe-shaped continuum emission peak. A possible mechanism to explain the decrease under the horseshoe-shaped continuum is the increased opacity coming from the higher dust concentration at the continuum peak. Lower {\gr dust and/or gas} temperatures and an optically thick radio-continuum reduce line emission by freeze-out and shielding of emission from the far side of the disk.

In view of both the size of its gap and the previously reported asymmetries and near-infrared spiral arms, the transition disk of the Herbig Fe star HD 142527 constitutes a remarkable case study. This paper focuses on the morphology of the outer disk through ALMA observations of $^{12}$CO J=2-1, $^{12}$CO J=3-2 and $^{13}$CO J=2-1. Both $^{12}$CO J=2-1 and $^{12}$CO J=3-2 show spiral features of different sizes. The innermost spiral arm (S1) is a radio counterpart of the first near-infrared spiral observed by Fukagawa et al. (2006), but it is shifted radially outward. However, the most conspicuous CO spiral arm (S2) lies at the outskirts of the disk and had not been detected before. It corresponds to a cold density structure, with both brightness and excitation temperatures of order 13$\pm$2 K and conspicuous in the $^{12}$CO J=2-1 peak-intensity map, but faint in $^{12}$CO J=3-2. There is also a faint counterarm (S3), point-symmetrical of S2 with respect to the star. These three spirals are modelled separately with two different formulae that approximate the loci of density maxima in acoustic waves due to embedded planets. S1 could be fit relatively well with these formulae, compared to S2 and S3. Alternative scenarios such as gravitational instability or external tidal interaction are discussed. The impact of channelization on spectrally and spatially resolved peak intensity maps is also briefly addressed.

Coronagraphy is a powerful technique to achieve high contrast imaging and hence to image faint companions around bright targets. Various concepts have been used in the visible and near-infrared regimes, while coronagraphic applications in the mid-infrared remain nowadays largely unexplored. Vector vortex phase masks based on concentric subwavelength gratings show great promise for such applications. We aim at producing and validating the first high-performance broadband focal plane phase mask coronagraphs for applications in the mid-infrared regime, and in particular the L band with a fractional bandwidth of ~16% (3.5-4.1 \mu m). Based on rigorous coupled wave analysis, we designed an annular groove phase mask (AGPM) producing a vortex effect in the L band, and etched it onto a series of diamond substrates. The grating parameters were measured by means of scanning electron microscopy. The resulting components were then tested on a mid-infrared coronagraphic test bench. A broadband raw null depth of 2 x 10^{-3} was obtained for our best L-band AGPM after only a few iterations between design and manufacturing. This corresponds to a raw contrast of about 6 x 10^{-5} (10.5 mag) at 2\lambda/D. This result is fully in line with our projections based on rigorous coupled wave analysis modeling, using the measured grating parameters. The sensitivity to tilt and focus has also been evaluated. After years of technological developments, mid-infrared vector vortex coronagraphs finally become a reality and live up to our expectations. Based on their measured performance, our L-band AGPMs are now ready to open a new parameter space in exoplanet imaging at major ground-based observatories.

Context. High contrast imaging has thoroughly combed through the limited search space accessible with first-generation ground-based adaptive optics instruments and the Hubble Space Telescope. Only a few objects were discovered, and many non-detections reported and statistically interpreted. The field is now in need of a technological breakthrough. Aim. Our aim is to open a new search space with first-generation systems such as NACO at the Very Large Telescope, by providing ground-breaking inner working angle (IWA) capabilities in the L' band. The L' band is a sweet spot for high contrast coronagraphy since the planet-to-star brightness ratio is favorable, while the Strehl ratio is naturally higher. Methods. An annular groove phase mask (AGPM) vector vortex coronagraph optimized for the L' band, made from diamond subwavelength gratings was manufactured and qualified in the lab. The AGPM enables high contrast imaging at very small IWA, potentially being the key to unexplored discovery space. Results. Here we present the installation and successful on-sky tests of an L'-band AGPM coronagraph on NACO. Using angular differential imaging, which is well suited to the rotational symmetry of the AGPM, we demonstrated a \Delta L' > 7.5 mag contrast from an IWA ~ 0".09 onwards, during average seeing conditions, and for total integration times of a few hundred seconds.

NACO is the famous and versatile diffraction limited NIR imager and spectrograph with which ESO celebrated 10 years of Adaptive Optics at the VLT. Since two years a substantial effort has been put in to understanding and fixing issues that directly affect the image quality and the high contrast performances of the instrument. Experiments to compensate the non-common-path aberrations and recover the highest possible Strehl ratios have been carried out successfully and a plan is hereafter described to perform such measurements regularly. The drift associated to pupil tracking since 2007 was fixed in October 2011. NACO is therefore even better suited for high contrast imaging and can be used with coronagraphic masks in the image plane. Some contrast measurements are shown and discussed. The work accomplished on NACO will serve as reference for the next generation instruments on the VLT, especially those working at the diffraction limit and making use of angular differential imaging (i.e. SPHERE, VISIR, possibly ERIS).