The SPHERE view of the planet-forming disk around HD100546 ††thanks: Based on data collected at the European Southern Observatory, Chile (ESO Programs 095.C-0273(A) and 095.C-0298(A)).
Key Words.:stars: pre-main sequence – planetary systems: protoplanetary disks – planetary systems: planet-disk interactions – ISM: individual object: HD100546 – Techniques: polarimetric
11email: email@example.com 22institutetext: Lunar and Planetary Laboratory, The University of Arizona, Tucson, AZ 85721, USA 33institutetext: Earths in Other Solar Systems Team, NASA Nexus for Exoplanet System Science 44institutetext: Departamento de Astronomía, Universidad de Chile, Casilla 36-D, Santiago, Chile 55institutetext: LESIA, Observatoire de Paris-Meudon, CNRS, Université Pierre et Marie Curie, Université Paris Didierot, 5 Place Jules Janssen, F-92195 Meudon, France 66institutetext: Sterrewacht Leiden, P.O. Box 9513, Niels Bohrweg 2, 2300RA Leiden, The Netherlands 77institutetext: CNRS/CRAL/Observatoire de Lyon/Université de Lyon 1/Ecole Normale Supérieure de Lyon, Lyon, France 88institutetext: Astronomical Institute Anton Pannekoek, University of Amsterdam, PO Box 94249, 1090 GE Amsterdam, The Netherlands 99institutetext: Univ. Grenoble Alpes, Institut de Planétologie et d’Astrophysique de Grenoble (IPAG, UMR 5274), F-38000 Grenoble, France 1010institutetext: CNRS, Institut de Planétologie et d’Astrophysique de Grenoble (IPAG, UMR 5274), F-38000 Grenoble, France 1111institutetext: Laboratoire Lagrange, Université Côte d’Azur, Observatoire de la Côte d’Azur, CNRS, Nice Cedex 4, France 1212institutetext: INAF - Osservatorio Astronomico di Padova, Vicolo dell’Osservatorio 5, 35122 Padova, Italy 1313institutetext: Max Planck Institute for Astronomy, Königstuhl 17, 69117 Heidelberg, Germany 1414institutetext: Department of Astronomy, Stockholm University, AlbaNova University Center, 106 91 Stockholm, Sweden 1515institutetext: UMI-FCA, CNRS/INSU France, and Departamento de Astronomía, Universidad de Chile, Casilla 36-D Santiago, Chile 1616institutetext: Università degli Studi di Padova, dipartimento di Fisica e Astronomia, vicolo dell’osservatorio 3, 35122 Padova 1717institutetext: Aix Marseille Université, CNRS, LAM - Laboratoire d’Astrophysique de Marseille, UMR 7326, 13388, Marseille, France 1818institutetext: European Southern Observatory, Alonso de Cordova 3107, Casilla 19001 Vitacura, Santiago 19, Chile 1919institutetext: Núcleo de Astronomía, Facultad de Ingeniería, Universidad Diego Portales, Av. Ejercito 441, Santiago, Chile 2020institutetext: Geneva Observatory, University of Geneva, Ch. des Maillettes 51, 1290, Versoix, Switzerland 2121institutetext: European Southern Observatory, Karl-Schwarzschild-Strasse 2, D-85748 Garching, Germany
Context:The mechanisms governing planet formation are not fully understood. A new era of high-resolution imaging of protoplanetary disks has recently started, thanks to new instruments such as SPHERE, GPI and ALMA. The planet formation process can now be directly studied by imaging both planetary companions embedded in disks and their effect on disk morphology.
Aims:We image with unprecedented spatial resolution and sensitivity disk features that could be potential signs of planet-disk interaction. Two companion candidates have been claimed in the disk around the young Herbig Ae/Be star HD100546. Thus, this object serves as an excellent target for our investigation of the natal environment of giant planets.
Methods:We exploit the power of extreme adaptive optics operating in conjunction with the new high-contrast imager SPHERE to image HD100546 in scattered light. We obtain the first polarized light observations of this source in the visible (with resolution as fine as 2 AU) and new H and K band total intensity images that we analyze with the pynpoint package.
Results:The disk shows a complex azimuthal morphology, where multiple scattering of photons most likely plays an important role. High brightness contrasts and arm-like structures are ubiquitous in the disk. A double-wing structure (partly due to ADI processing) resembles a morphology newly observed in inclined disks. Given the cavity size in the visible (11 AU), the CO emission associated to the planet candidate c might arise from within the circumstellar disk. We find an extended emission in the K band at the expected location of b. The surrounding large-scale region is the brightest in scattered light. There is no sign of any disk gap associated to b.
Our knowledge of the processes governing planet formation will be greatly enhanced by new generation instruments that recently started operations, like the SPHERE (Spectro-Polarimeter High contrast Exoplanet REsearch, Beuzit et al. 2008) instrument at the Very Large Telescope (VLT) and GPI (Gemini Planet Imager, Macintosh et al. 2014). The high-contrast and -resolution images enabled by these facilities will continue the recent plethora of observations of planet formation caught in the act. These consist of both the detection of planetary companions still embedded in disks (e.g., Kraus & Ireland 2012; Biller et al. 2014; Reggiani et al. 2014; Sallum et al. 2015) and of their imprints on the natal environment (e.g., Hashimoto et al. 2011, 2012; Quanz et al. 2013b; Garufi et al. 2013). Imaging disk features (such as disk cavities, spirals, rings etc.) that can be due to the interaction with forming planets is fundamental to determine the framework (and the timeframe) of the planet formation. Excellent first examples of SPHERE’s capability to image protoplanetary disks were provided by Thalmann et al. (2015) and Benisty et al. (2015).
The disk around HD100546 (B9V star at 97 4 pc, Levenhagen & Leister 2006; van Leeuwen et al. 2007) is one of the best laboratories to study the interaction with forming planets since it hosts two embedded companion candidates. The existence of HD100546b (hereafter planet ) at has been proposed by Quanz et al. (2013a), who detected a bright point source sitting on top of an extended emission in the L band. This detection was later confirmed by Currie et al. (2014). Quanz et al. (2015) obtained estimates on its temperature, emitting radius and luminosity, proposing the existence of a warm circumplanetary disk. More recently, Currie et al. (2015) reported the detection of emission in the H band at the location of b. On the other hand, the presence of HD100546c (planet ) has been claimed by Brittain et al. (2013) via spectroastrometric studies of the CO and OH ro-vibrational line emission. These authors ascribed the asymmetry of the OH line profiles to gas emission in an eccentric orbit at the disk wall and the annual CO line variability to a concentrated source of emission orbiting the star in proximity of the circumstellar disk wall at AU. The latest fit to these variations gives an emitting area of AU at (when at orbital phase from the major axis, Brittain et al. 2014). They noted that this emitting area is only slightly larger than the expected size of a circumplanetary disk around a 5 M planet at this orbital separation. Fedele et al. (2015) showed that the asymmetric profile of the OH lines can be explained by a misalignment of the spectrograph slit with no need to invoke an eccentric gas disk.
An extensive literature on this remarkable protoplanetary disk exists. The large, roughly -inclined disk was firstly resolved in scattered light by Pantin et al. (2000). Peculiar disk structures (as the inner cavity, dark lanes and spiral arms) have been resolved in scattered light by Augereau et al. (2001), Grady et al. (2005), Ardila et al. (2007), Quanz et al. (2011), Boccaletti et al. (2013), and Avenhaus et al. (2014a). A quasi-coplanar inner disk at (sub-)AU scale has been studied by Benisty et al. (2010) and Tatulli et al. (2011) with near-IR interferometry. A rounded disk wall at 11 AU has been claimed by Panić et al. (2014) through mid-IR interferometry. The disk was also resolved with mid-IR imaging (Liu et al. 2003) and at millimeter wavelengths by both ALMA and ATCA. The ALMA images (at moderate angular resolution) suggest that the large dust grains are mainly located in form of a ring between the radial locations of planet and (Walsh et al. 2014; Pineda et al. 2014). Pinilla et al. (2015) modeled these observations and found that if two planets at 10 and 70 AU are responsible for confining the dust in the observed ring, then the outer planet should be at least 2.5 Myr younger than the inner planet. Finally, ATCA observations at 7 mm (with resolution as good as 0.15″) constrained the disk cavity size to be 25 AU (Wright et al. 2015). Those observations also reveal a horseshoe-shaped concentration of grains at the disk inner edge.
We present the first SPHERE images of HD100546 in scattered light. The observations reported in this paper consist of high-contrast polarized light images in the visible as well as total intensity images in the near-IR obtained with the SPHERE sub-instruments ZIMPOL and IRDIS respectively. The ZIMPOL observations are the first polarized light images of HD100546 in the visible and are among the highest-resolution direct imaging data () of a protoplanetary disk ever obtained. This paper is organized as follows. In Sect. 2 we describe observing conditions and data reduction, in Sect. 3 we present the results from both the ZIMPOL and IRDIS images, and in Sect. 4 and 5 we discuss our findings and provide our main conclusions.
2 Observations and data reduction
HD100546 was observed in the context of the Guaranteed Time Observations of the new high-contrast imager SPHERE (Beuzit et al. 2008), operating at the VLT in conjunction with the extreme adaptive optics (AO) system SAXO (Fusco et al. 2006). This paper presents the observations performed with the sub-instruments ZIMPOL (Zurich IMaging POLarimeter, Thalmann et al. 2008) and IRDIS (Infra-Red Dual-beam Imager and Spectrograph, Dohlen et al. 2008). Further IFS (Integral Field Spectrograph, Claudi et al. 2008) observations were taken along with the IRDIS observations (and will be presented by Sissa et al. in prep.). The observing conditions and the reduction of the ZIMPOL and IRDIS data are described in Sect. 2.1 and Sect. 2.2, respectively. We also retrieved archival VLT/NACO and Gemini/NICI datasets of HD100546, which are described in Sect. 2.3.
The ZIMPOL observations were performed on 2015, April 23 (night 1) and May 7 (night 3) in the R band ( nm) in Differential Polarization Imaging (DPI). In this mode, ZIMPOL allows polarimetric diffraction-limited observations with very high polarimetric sensitivity. A fast polarization modulator (kHz) provides quasi-simultaneous observations of opposite polarization states, which are registered on the same ”even” detector rows during two consecutive modulation cycles. The charges from the first polarization cycle are shifted to the ”odd”, covered rows and the information about the second state are then collected in the even rows. Furthermore, a half-wave plate (HWP) controls the polarization orientation and permits the observations of a full polarization cycle, consisting of the Stokes parameters , , , (see e.g., Tinbergen 2005).
All our observations were performed in DPI field stabilized mode. We alternated the object orientation by and dithered its position by 14 pixels on the detector to evaluate any possible instrument artifacts. During night 1, HD100546 was mostly observed in FastPolarimetry mode, which provides the highest polarimetric precision with short exposure time (DIT=1.2 sec). Atmospheric conditions were initially challenging due to the transient presence of thin clouds but converged to good conditions later on (with optical seeing varying from 0.9 to 2.4). During the last observing block of night 1 and during all observations of night 3, we operated in SlowPolarimetry (DIT=20 sec) and employed a 155 mas-diameter classical Lyot coronagraph. The weather conditions during night 3 were significantly affected by a strong wind, which resulted in a largely variable seeing (from 0.9 to 2.2). A summary of the ZIMPOL observations is given in Table 1.
The DPI data reduction follows the technique illustrated by Avenhaus et al. (2014b). This includes the equalization of each pair of quasi-simultaneous polarization states to compensate for possible instrumental polarization. This technique assumes the star to be unpolarized. The astrometric calibration was performed following Ginski et al. (in prep.), i.e. by adopting a plate scale of mas/pixel and a detector orientation of . The final product of our reduction is the pair of polar Stokes parameters (, )111which was sometimes referred to as (, ), (, ), or (, ).. contains the polarized light component tangential to the star on the image plane. provides polarization vectors 45 inclined to the tangential component. For face-on systems only, is equivalent to but it is unbiased (since the signal is not artificially increased by the squares), whereas is not expected to contain any signal. In Sect. 4.2 we discuss the validity of this statement for inclined disks. Since the PSF in many exposures was degraded by the high atmospheric turbulence, the final images were produced by stacking a visual selection of best frames (see Table 1). This operation significantly decreased the effective exposure time but also clearly improved the contrast of the final images.
The IRDIS observations were taken in parallel with IFS. This possibility is enabled in IRDIFS mode by dichroic beam splitters and represents the SPHERE nominal infrared mode (with IFS operating at Y-J wavelengths and IRDIS in the H band). Observations in this mode, with IRDIS working in DBI mode (Dual-band imaging mode, Vigan et al. 2010) with the H2 and H3 filters ( nm and 1667 nm), were taken on 2015, May 28 (night 4). Furthermore, observations in IRDIFSEXT mode (providing IFS coverage up to the H band and IRDIS working in the K1 and K2, nm and 2255 nm) were taken on 2015, May 3 (night 2). Both runs were carried-out in pupil stabilized mode, allowing us to perform Angular Differential Imaging (ADI, Marois et al. 2006). An apodized pupil Lyot coronagraph with 185 mas-diameter was employed during both runs. During both night 2 and 4, the atmospheric conditions were good, with an average seeing of 0.8″. The individual frame integration times were 16 seconds and in total 256 frames were obtained in both H and K bands.
The basic data reduction (e.g., bad pixel cleaning, flat fielding, image alignment) was performed using the SPHERE Data Reduction and Handling (DRH) pipeline (Pavlov et al. 2008). From astrometric reference measurements during the same observing run, we derived a pixel scale of mas/pixel and a true North orientation of (Maire et al. 2015). Visual frame selection was carried-out for both the H and K band data to sort out frames with poor AO correction leading to 121 frames in the H band and 21 frames in the K band that we disregarded from further analysis. The final selection of frames provides a field rotation of in the H band and of in the K band. A summary of the IRDIS observations is given in Table 1.
To subtract the stellar contribution and reveal the circumstellar material, we used the pynpoint package (Amara & Quanz 2012; Amara et al. 2015). pynpoint uses Principal Component Analysis (PCA) to model the stellar contribution in each frame using a basis set of principal components (PCs) that was created from all existing frames. This was done for each filter independently. After stellar subtraction, all frames were rotated to a common sky orientation and mean combined. We varied the number of PCs and the inner/outer radius over which the fit is performed to estimate the impact of these parameters. For the final analysis we eventually settled on simply mean subtracted images (the main PC) as fitting a higher number of PCs led to more flux subtraction and did not reveal any significant additional structure. On the other hand, changing the inner and outer radius of the input images helped to identify features at different spatial scales. In this paper, we present the results obtained with inner/outer diameters equivalent to 0.14″ / 1.0″ (small scale), 0.5″ / 1.5″ (medium scale), 1.25″ / 5.0″ (large scale).
2.3 Literature data
We also make extensive use of archival data available for HD100546. We retrieved H ( nm) and K ( nm) band polarized light images of the source taken with VLT/NACO in April 2006 and in March 2013 (presented by Quanz et al. 2011; Avenhaus et al. 2014a). We also recover NACO ADI L ( nm) and M ( nm) band data taken in April 2013 and published by Quanz et al. (2015). Finally, we make use of Gemini/NICI K band images of HD100546 taken in March 2010 and reduced with the algorithm LOCI (Lafrenière et al. 2007) by Boccaletti et al. (2013).
The ZIMPOL and images resulting from our data reduction are described in Sect. 3.1 and 3.2, respectively. The IRDIS images are illustrated in Sect. 3.3. A comparison between this dataset and the VLT/NACO as well as the Gemini/NICI images is shown in Sect. 3.4.
3.1 Zimpol images
Figure 1a and 1b shows the clear detection of the disk around HD100546 in the images obtained on night 3 and 1 respectively. The images show a globally complex structure. A signal is detected down to radii as small as 0.08″-0.11″ (depending on the azimuthal angle), which is a factor larger than the inner working angle of the non-coronagraphic data. Along the major axis, two bright quasi-symmetric lobes lie outward of the detection inner edge. Farther out, polarized light with an azimuthally complex distribution is detected at radii as large as 0.6″-1.6″. Radial and azimuthal brightness profiles from the images are shown in Fig. 2. These profiles are obtained by averaging the counts contained in a resolution element (i.e. pixels). Errors are extracted from the images as the standard deviation of the same box divided by the square root of the number of pixels therein. In this section, we analyze the brightness distribution of the images. As from previous interpretations (e.g. Quanz et al. 2011), we will refer to the region close to the star with no detectable signal as the disk cavity (Sect. 3.1.1) and to the signal at larger radii as the outer disk (Sect. 3.1.2).
3.1.1 Disk cavity
The cavity is elongated along the NW-SE direction, with the major axis lying at East of North. Along this direction, the inner rim reaches its maximum emission on both sides at 0.13″ but the brightness of the NW side is lower (see Fig. 2a). Along the SW side (i.e. the disk near side, Avenhaus et al. 2014a), the signal is on average lower and the intensity peaks at 0.10″. This spatial difference is consistent with a circular cavity which is inclined. Along the NE side (the far side), the emission is significantly weaker (a median factor 2.3 lower) but the signal also peaks at 0.10″. However, from Fig. 2a it is clear that the radial profiles along the near and the far sides are different. Along the near side, the slope is similar to the major axis, except for an abrupt decrease at 0.19″. On the other hand, the slope along the far side is much more shallow. Furthermore, a brightness bump is visible from 0.17″ to 0.25″. This bump can be also appreciated in Fig. 1e as a spiral arm (see Sect. 3.1.3).
In Fig. 2b, we show the ZIMPOL radial profile of the inner 0.3″ along the major axis and compare it with the NACO H and K band profiles. A small radial offset seems to exist between the two datasets. To evaluate this possible discrepancy, we adopted the technique described by Thalmann et al. (2015). Shortly, we integrated the flux contained in 3-pixel large concentric elliptic annuli, which are obtained from the projection of circular rings ( and , Pineda et al. 2014). We searched for the largest increase in the integrated flux over two contiguous annuli and considered it as the location of the disk inner edge. In parallel, we also generated a disk toy model with surface brightness proportional to outside of a defined radius and only Poisson noise (quantitatively mimicking the image noise) inward. Then, we convolved the synthetic image with a two-dimensional Gaussian with the Full Width Half Maximum (FWHM) of our observations and changed the disk inner edge over a discrete range of values to qualitatively reproduce the signal distribution observed along the major axis.222The elliptic annuli technique may yield different results by changing the adopted inclination and position angle. However, we did not find any relevant difference over a reasonable range of values. The results from the toy model are valid under the assumption that the disk inner edge is sharp, which is not necessarily true (see Mulders et al. 2013b; Panić et al. 2014). However, for the purpose of comparison of the two datasets this uncertainty is not an issue. Both techniques led to a inner radius of . For consistency, we applied the same techniques to the NACO datasets and found for both H and K bands. Even though the error bars from these estimates overlap, the small offset from the radial profiles derived from the two datasets might be real rather than due to the different angular resolution of the observations.
We also extracted the brightness distribution along the inner rim (Fig. 2c) by averaging the contribution from a 5 pixel-wide ellipse at 0.15″, obtained by projecting a circular ring by . We centered the ellipse slightly outside of the intensity peak (0.13″) to include only regions with detectable signal. The location of the main peak well matches the major axis. On the contrary, the peak to the NW is offset toward North (consistently with the NACO dataset, Avenhaus et al. 2014a). The same amount of offset is found between the intensity peak on the near side and the location of the minor axis.
3.1.2 Outer disk
The azimuthal brightness distribution at large radii is highly asymmetric. The NE half of the image (the far side) is significantly brighter than the SW half (the near side). The distribution on the far side seems to be also discontinuous. As highlighted in Fig. 1c, a dark wedge to East separates two wide bright regions (centered at and ). On the near side, a bright wedge (spanning -) stands out. A sharp bright blob is also visible to West (marked as AO artifact in Fig. 1c). This is most probably an instrument artifact as it does not appear in differently sky-oriented frames while it is seen in the total intensity images at the radial location of the main AO ring (corresponding to the AO correction radius at /D, inside of which the AO system provides almost perfect corrections).
In Fig. 2d we show the normalized azimuthal profile from both NACO and SPHERE images, obtained at (i.e. the predicted radial location of b, Quanz et al. 2015) similar to Fig. 2c. Interestingly, the intensity peak from the three wavebands lie at the same azimuthal position, which is farther East of North of HD100546b. The normalized distributions from the three wavebands are consistent within the error bars almost everywhere. The only exception is the SW bright wedge from the SPHERE image, which is roughly 40% brighter than the same structures from the NACO datasets. The azimuthal width of this wedge is similar to what we found along the rim (Fig. 2c) but it is displaced.
3.1.3 Disk structures
No significant sub-structures can easily be spotted from the outer disk. Only a spiral arm to the NW at can be seen from Fig. 1a and Fig. 2a. To reveal any additional elusive features in the disk, we applied an unsharp masking technique to our images. This technique consists in adding a blurred, negative version of the original image to sharpen its details. It has been applied to protoplanetary disks images by e.g., Ardila et al. (2007) and Quanz et al. (2011). The resulting image must be interpreted with caution as, although sharper, it is most certainly not a more accurate representation of the real disk structure.
An unsharp masked version of the SPHERE image was obtained by subtracting the smoothed version (by ) from the original image. Since the noise in the science image is high, we firstly subtracted the variance over a FWHM from the original image. This operation had the effect of decreasing the signal of the inner regions so as to avoid over-subtraction from the smoothed image, and resulted in a much sharper image of the inner 0.5″. Finally, the resulting image has been smoothed by one FWHM. The morphology of the final image strongly depends on the choice of the mentioned parameters. However, some features are visible for a broad range of parameters. As shown in Fig. 1e, these are the inner rim, the spiral to NE, and an arm-like structure to North. The reality of this last feature, although not necessarily compelling by itself, will lend credence from the comparison with other datasets (see Sect. 3.3 and 3.4).
3.2 Zimpol images
The image resulting from our data reduction is shown in Fig. 1d. This image seems to contain a significant signal pattern. The intensity and the morphology of this signal is not likely to be due to an imperfect data reduction as the equalization of polarization states (minimizing instrumental polarization, see Sect. 2.1) cannot eliminate such a signal and the crosstalk effect between Stokes parameters (see Avenhaus et al. 2014b) is a minor effect in SPHERE. Furthermore, the spatial distribution of this signal is consistent through instruments and wavebands (see the NACO images by Avenhaus et al. 2014a), with a positive signal being diffusely detected to the East and to the West, and with a negative signal to the NW and SW. Both the bright and the dark wedge highlighted in Fig. 1c seem to have a counterpart in the SPHERE image. However, a offset exists in the azimuthal location of the bright wedge from the and images (it can also be seen in Fig. 2d).
More quantitatively, the signal in the SPHERE image varies from of the image at 0.2″ to at 0.8″. These estimates are obtained by summing up the absolute values of all counts contained in the respective circular annulus. Similar values are found by comparing the peak-to-peak signal at those radii. We will discuss these findings in Sect. 4.2.
3.3 IRDIS images
The IRDIS images resulting from the reduction described in Sect. 2.2 are shown in Fig. 3. We found no significant differences between the H2 and the H3 bands as well as between the K1 and the K2 bands. Thus, we only show their combinations (i.e. H2H3 and K1K2 - or H and K for simplicity). The overall brightness morphology in the H and K images looks very similar but some minor, though important, differences exist. Different structures are revealed from the reductions with small, medium, and large scales (see Sect. 2.2).
Small scale. The main features in the small scale images (first column in Fig. 3) are two bright arms with wide pitch angles (the wings), which are roughly symmetric around the disk minor axis. Signal is detected almost down to the software mask radius at . A bright knot is visible along the SE arm at (P.A.=) from both wavebands (but it is brighter in H). There is evidence of at least two additional arms to SE with comparable pitch angle. One of these (spanning P.A. ) seems to originate from the bright knot of the SE arm. Interestingly, the spiral arm at 0.2″ from the DPI data (see Fig. 1e) is not detected, whereas the faint northern arm matches the location of the Northern wing.
Medium scale. All arms revealed in the small scale images are also visible at medium scale (second column in Fig. 3). A diffuse enhanced brightness is detected in the K band at P.A. . On top of this extended emission, a more localized emission spans radii from 0.43″ to 0.47″. Interestingly, the region of extended emission roughly lies at the end of the northern arm. We do not detect such an enhanced emission in the H band.
Large scale. Many bright arms are also detected in the large scale images (third column in Fig. 3). The two clearest detections lie to South in both images (the Southern spirals spanning angles ). Two additional smaller arms can be seen at P.A. , close to the easternmost part of the main arms. Moreover, we detect in both bands two similar features to North (the Northern spirals which run in parallel from to ) and also tentatively two arms which are specular to the Southern spirals around the minor axis (the Western spirals). Any other structure close to the software mask should not be considered significant, as that region is dominated by subtraction residuals and none of those features is persistent across different frames and reductions.
3.4 Comparison with near-IR images
Having collected throughout the years a rich set of data from 0.6 m to 4.8 m, both in polarized light and total intensity, we can obtain insight into both the disk geometry and the properties of scattered light by comparing this image collection. Figure 4 is a selection of these comparisons that we discuss in this section.
Figure 4a shows the brightness distribution in the H band, highlighting the differences between the polarized light (from NACO, in color) and the total intensity distribution (from IRDIS, contours). It is clear that our IRDIS images can trace the signal as far down as the inner rim in polarized light, but that this is slightly azimuthally offset toward the near side. It is also evident how different the two brightness distributions are at large radii, with the polarized light being more distributed over the angles. Nonetheless, neither of the datasets reveal any significant emission from the SW side, and the locations of the bright/dark boundaries are very well matched.
In Fig. 4b, the shortest (R) and the longest (M) available wavebands are compared. The two bright large-scale wedges observed in the visible (colors) are also evident from the IR NACO data (contours). Particularly good is the match between the bright/dark edge to East. Also worth mentioning is the tentative detection of the spiral at 0.2″ (Fig. 1e) from the M dataset.
In Fig. 4c, we show the spatial connection between the northern arm spotted from the ZIMPOL image after unsharp masking (colors) and the diffuse emission from the M band (cyan contours) surrounding planet b (Quanz et al. 2015). The same emission from the IRDIS K band is shown (yellow contours). The arm seems to head to the location of b and a blurred counterpart of this arm may be seen in the M band. A minor radial offset between the emission in the K and M band seems to exist but requires further investigation. The inset image to the top right is the K band NACO from 2006 (Quanz et al. 2011) after applying unsharp masking similarly to the ZIMPOL image. It is clear that the location of the spiral is, to first order, the same.
The interpretation of scattered light images is not straightforward. In fact, many physical and geometrical factors may equally contribute to the final appearance of such images. First of all, dust particles scatter photons anisotropically. This scattering angle distribution (referred to as the phase function) depends on the incident radiation as well as on the dust size and shape (see e.g., Min et al. 2012). In addition, scattered photons are polarized also depending on the scattering angle (hereafter the polarizing efficiency). In protoplanetary disks, the scattering angles experienced by the observed photons cannot be precisely determined because of our limited knowledge of the exact disk geometry. Secondly, the amount of radiation incident on the disk may be strongly affected by the geometry of the system at smaller radii (like e.g., warps or belts in proximity of the star). Moreover, multiple scattering on the disk surface may also have an impact on the scattered light distribution. Finally, instrumental contributions or data processing can also significantly alter the final images.
All this might suggest caution in interpreting the observed features, as very often these may have twofold (or multiple) explanations. In this section we discuss our findings in the context of different possible scenarios, focussing on the inner disk (Sect. 4.1), on the global disk morphology (Sect. 4.2), on the disk features (Sect. 4.3), and on the environment around planet b (Sect. 4.4).
4.1 Cavity and inner rim
Our results on the disk cavity (from Sect. 3.1.1) support the previous findings in DPI (Quanz et al. 2011; Avenhaus et al. 2014a), that are that the disk around HD100546 is truncated at roughly 15 AU and that this cavity is consistent with being intrinsically circular. Our data show no sign of offset along the major axis between the center of the cavity and the star, contrarily to what was found by Grady et al. (2005) with spectral observations (). Along the minor axis, it is harder to claim or rule out any offset because the scattered light distribution from the back side of inclined disks inherently differs from that of the near side.
Mulders et al. (2013b) showed with mid-IR interferometry that the disk wall is rounded off over a large radial range and that this creates a broad surface brightness profile (10 to 25 AU) peaking at 12 AU. Our radial profile of Fig. 2a may support this scenario. In fact, the brightness distribution along the major axis is constantly decreasing between the peak at 13 AU and 30 AU. This trend differs from those of e.g. SAO206462 (Garufi et al. 2013) and HD169142 (Quanz et al. 2013b), where an abrupt decline of the brightness profile is seen outward of the intensity peak. However, the exact slope of the brightness distribution at such small radii should not be trusted, because of the PSF smearing effects described by Avenhaus et al. (2014a).
The radial brightness distribution on the near side is very similar to what is seen along the major axis, except for an abrupt discontinuity (almost 50%) at a de-projected radius of 25 AU. This may indicate both a shadowed region or a change in the dust distribution at the disk surface. On the other hand, the profile on the far side is completely different. The shallow slope (clear from Fig. 2a) may suggest that from 0.1″ to 0.2″ our observations (partly) trace the disk wall, which is presumably seen at a angle from face-on. Such an exquisite data quality may inspire a later, robust modeling of the inner disk geometry to investigate this scenario.
Another intriguing result is the possible discrepancy between the location of the disk inner edge in the visible (11 AU) and in the near-IR (13 AU). This incongruity, though marginal, may be the result of a gradual increase of optical depth at the inner rim with disk opacity varying with the wavelength. To explore this scenario, we calculated the spectral index of necessary to explain this discrepancy. We imposed the observed 11 AU and 13 AU edges as the yielding in the visible and near-IR respectively, for a range of realistic values of . This exercise resulted in values spanning from -0.5 to -1.3, which is consistent with the expectations from (sub-)m sized particles (e.g., Backman & Paresce 1993). We note that the disk inner edge inferred from the visible dataset matches the location of the bright ring of emission ( AU) claimed through mid-IR interferometry by Panić et al. (2014).
Interferometric ATCA images of HD100546 (Wright et al. 2015) suggest that the cavity size at millimeter wavelengths is roughly 25 AU, which is more than twice as large as our estimate in the visible. Such a large discrepancy is a probable indication of the different behavior of m- and mm-sized dust grains in gapped disks. In fact, many authors have shown that a pressure bump at the disk inner edge may act as a filter, allowing grains smaller than a certain size to drift inward but efficiently trapping larger grains (e.g., Rice et al. 2006; Brauer et al. 2008; Dong et al. 2012). Pinilla et al. (2012) have shown that two distinct cavity sizes are expected from the interaction of a disk with a planet, one for the gas (closer to the companion) and one for the millimeter particles farther out. This morphology has been observed in a number of disks by comparing scattered light images (tracing m-sized grains which are well coupled to the gas) and millimeter interferometric images (e.g., Garufi et al. 2013; Tsukagoshi et al. 2014; Follette et al. 2015). More specifically, de Juan Ovelar et al. (2013) simulated visible and millimeter images of disks with a cavity carved by planets with different masses and orbital radii. Their results provide a fitting function to relate at different radii the visible/millimeter cavity size ratio with the planet mass. In the case of HD100546, this ratio is observed to be 0.5. Extending their results to AU, one finds that a companion responsible for this dust differentiation should be more massive than 15 . Adopting an analytic solution to this difference as in Garufi et al. (2013) (after the model by Pinilla et al. 2012) leads to an even higher lower limit (). Similar values were estimated by Mulders et al. (2013b) from the aforementioned shape of the disk wall and by Pinilla et al. (2015) by modeling ALMA observations. The scenario where multiple planets contribute to sculpt the disk and to differentiate the dust grains, though a valid explanation, has not been extensively studied because of the degeneracies introduced.
Given these considerations, the potential companion c associated with the compact CO emission by Brittain et al. (2013) is unlikely to be the (unique) cause of the observed cavity. From Fig. 2b, it is clear that such an object would lie within the m-sized particles halo, disfavoring a scenario where these particles are retained at 11 AU because of the interaction with this object. Furthermore, the size of this CO emission would be associated with a planet with a relatively low mass (, Brittain et al. 2014), which is not consistent with the large cavity size nor with the observed m-mm dust differentiation in the framework of the current models. An intriguing possibility is that this object is actually a consequence of the observed cavity, since the accumulation of material is favored by the pressure bump at the disk inner edge (e.g., Pinilla et al. 2012).
One would be tempted to relate the location of c in May 2015 (see Fig. 1e) with both the increased brightness of the SE side of the rim and the possible starting point of the main spiral. However, these associations may be fortuitous as both features were morphologically similar in the NACO dataset from April 2006 (Avenhaus et al. 2014a, see also Sect. 4.3), when the companion candidate was orbiting at angles close to North (Brittain et al. 2013). Recent H band images of HD100546 suggest the putative detection of c at angles of (Currie et al. 2015). The location of this detection roughly matches the bright knot of the SE arm from the IRDIS images at the convergence of another arm (see bottom right panel of Fig. 3). However, the morphological similarity with the NW arm (both in brightness and symmetry around the minor axis) and the absence of such a knot in the K band images may suggest that the signal detected at those locations is pure disk emission. This discourages us from further analyzing and interpreting this feature. All in all, our dataset does not firmly reveal the presence of planet c nor of any disk feature which may be connected to the interaction with the planet candidate.
4.2 Disk morphology and scattering properties
The complex azimuthal distribution of the scattered light from HD100546, from both ZIMPOL and IRDIS, can only partially be explained by the scattering phase function. In this section, we discuss what other factors may contribute to the observations.
vs images. The presence of some diffuse features of the image (like the bright and dark wedges) also in the images (see Fig. 1d) may cast doubt on their intrinsic existence. In fact, the parameter has largely been used to determine the noise level of polarimetric observations (e.g., Garufi et al. 2013; Avenhaus et al. 2014a), since by construction it was not thought to contain any scattering information. However, the amount ( to of the images) and the consistency (between different wavebands and instruments) of the image of HD100546 brought this view into question. Bastien & Menard (1988) showed that in optically thick, inclined disks the effect of multiple scattering may lead to significant deviations from purely azimuthal linear polarization. Recently, Canovas et al. (2015) have produced synthetic images of disks with different inclinations and masses, showing that these images contain a strong signal in disks with . The morphology of the synthetic images resemble that of HD100546, with the near side of the disk showing an alternation of strongly negative and positive signal. All this may suggest that multiple scattering in the disk of HD100546 may act to transfer a fraction of the polarized signal from the to the image and therefore that the presence of features in both images does not necessarily discredit their intrinsic existence.
Inner rim. From Fig. 1a and Fig. 2c, it is clear that the polarized emission at the inner rim is maximized along the major axis. An analogous morphology is observed in other inclined disks (e.g., Hashimoto et al. 2012; Garufi et al. 2014). This is not surprising because the polarizing efficiency of scatterers is typically maximized for scattering angles around (e.g., Murakawa 2010). It is also clear from Fig. 4a that the (non-polarized) emission from IRDIS is maximized at angles somewhat closer to the minor axis. This indicates that dust grains at the inner rim preferentially scatter in a forward direction. However, the net amount of this trend must be small compared to the polarizing efficiency, so as to maximize the polarized light (which is a combination of phase function and polarizing efficiency) along the major axis.
Polarized light from the (bright) far side. At larger radii (), the polarized emission is predominantly distributed on the disk far side. This morphology can be explained by the disk geometry. In fact, for inclined and flared disks (like HD100546) the scattering angles experienced by photons from the back side are closer to . However, the high brightness contrast between the near and the far side and the darkness of the near side also in the IRDIS (non-polarized) images may both suggest backward-scattering particles at the disk surface (see also work by Avenhaus et al. 2014a). Specifically, two brighter regions are clear (at angles around and ) and these are also visible in the M band (see Fig. 4b). These regions are offset from the major axis by a different angles respectively ( to North and to SE). This may indicate a larger disk flaring angle to North, since higher disk scale heights move the -scatters toward the far minor axis333One can relate the scattering angle to the disk geometry through with disk inclination, P.A. position angle from the major axis, and disk opening angle.. This scenario is consistent with the different brightness of the two regions (with the North being twice as bright) and would also explain why the offset with the major axis increases with the radius ( at the rim and up to for higher disk scale heights). On the other hand, the dark wedge to East (see Fig. 1c) is too localized to be described by a change in the disk scale height. A possible explanation for this deficit is a shadow. Shadows in protoplanetary disks have been predicted (e.g. Dullemond et al. 2001) and observed (e.g. Avenhaus et al. 2014b). The compact inner dust belt (e.g., Benisty et al. 2010; Mulders et al. 2013b) may be responsible for such a shadow, similar to HD142527 (Marino et al. 2015). An alternative cause could be the spiral at 20 AU (see Fig. 1e), whose inner (brighter) portion is azimuthally consistent with the dark wedge. In this scenario, the dark wedge should rotate with the spiral (see Sect. 4.3).
Polarized light from the (dark) near side. The dark region to the SW is also challenging to interpret. This discussion must be related to the interpretation of the bright wedge therein (see Fig. 1c) which is brighter than in the NACO dataset (Fig. 2d). Avenhaus et al. (2014a) suggested that the dark region is the result of a particularly backward-peaked scattering phase function (where photons from this region are scattered by angles as small as ). In this scenario, the bright wedge could represent the low tail of a forward-scattering peak, which is typically estimated to be in the range (e.g., Min et al. 2016). These scattering angles are typically not observed in moderately inclined disks. However, a combination of high flaring angle and broad peak of the phase function may enable this in HD100546. Alternatively, this darkness may also be due to a large-scale shadow. Wright et al. (2015) revealed from millimeter observations a horseshoe structure lying at the SW disk inner edge and claimed that such an accumulation of particles may shadow the outer disk. However, our data trace the dark region down to 20 AU whereas the inner edge of the millimeter emission by Wright et al. (2015) is measured to be 25 AU. Furthermore, the presence of the bright wedge (highlighted in Fig. 1c) is difficult to reconcile with this scenario. It could be a region of penumbra due to asymmetries in the millimeter horseshoe, but the presence of this bright wedge stretching down to the disk inner rim (see Fig. 2c) disfavors this view. Nonetheless, these arguments do not rule-out a shadow cast by the aforementioned inner dust belt. A good vehicle to solve this ambiguity is inferring whether the different intensity between the SPHERE (visible) and the NACO (near-IR) datasets is due to the different epochs (which would favor the shadow scenario) or to the different wavebands (which would suggest a small dependency of the phase function on wavelength). Therefore, new near-IR observations (with e.g., SPHERE/IRDIS in DPI mode) will probably clarify the nature of the dark region.
Structures in the ADI images. The brightness distribution from the IRDIS images significantly differs from that of polarized light. From the small and medium scale images of Fig. 3, it appears as a double-wing structure rising from the major axis and arching toward the far side (similarly to the GPI images by Currie et al. 2015). Contrarily to the polarized light, the region between the two wings (to NE) is dark. To explore to what extent this morphology can be generated by the ADI processing, we produced synthetic images of HD100546 in the H band from the model by Mulders et al. (2013a) with Henyey-Greenstein asymmetry parameter (Henyey & Greenstein 1941) spanning the range from 0.5 to -0.2. Then, we convolved these images with the angular resolution of the IRDIS observations, replicated and rotated them to reproduce the field rotation of our observations (see Sect. 2.2) and gave them as input to pynpoint, which was run with the same setup described in Sect. 2.2 (with the small scale boundaries). Some illustrative outputs of this exercise are shown in Fig. 5. It turned out that the ADI processing transforms a continuous brightness distribution into a double-wing structure, with the wings always rising from the major axis and the azimuthal distance between them being strongly dependent on . None of the test models could reproduce the observed morphology. Speculatively, this supports the idea of a strong backward-peaking phase function. However, a deeper discussion is beyond the scope of the current paper. Broadly speaking, this double-wing structure is a new form of feature (see also AK Sco, Janson et al. 2016) that may be recurrent among the observations of inclined disks carried-out with the new generation AO systems. Our exploration has shown that such a peculiar structure may still be interpreted in the context of a disk with an azimuthally continuous brightness distribution.
4.3 Disk structures
Almost all protoplanetary disks imaged so far with high resolution show peculiar features. Spirals (e.g., Muto et al. 2012; Garufi et al. 2013; Benisty et al. 2015) and annular gaps (e.g. Debes et al. 2013; Quanz et al. 2013b; Rapson et al. 2015) are, among these structures, the most intriguing. However, none of these disks seem to show both types of structures. The disk around HD100546 only shows spiral arms. Nonetheless, the absence of annular gaps is not less important in the context of planet-disk interaction and is discussed in Sect. 4.4.
The inner spiral arm at (see Fig. 1e and Avenhaus et al. 2014a) is evident in polarized light but is not in the ADI images. This non-detection is most likely due to the ADI processing, which tends to cancel out any disk feature with high azimuthal symmetry and to over-subtract flux along the disk minor axis (see Fig. 5). The aperture of this spiral significantly differs from that of all others. In fact, if we define the aperture as from (with distance from the star in arcseconds and azimuthal angle in radians444To obtain , we de-projected the disk following Pineda et al. (2014, and ) and considered a disk opening angle of 10.), this spiral is fitted by whereas the other features farther out show as high as . The most interesting finding about this spiral is the absence of rotation on a nine-years timescale (see comparison with NACO data in Fig. 4c). This is surely not consistent with a Keplerian motion around the star. In fact, a spiral at a de-projected separation of 30 AU from a 2.4 star (van den Ancker et al. 1997), has an orbital period of roughly 100 years, which yields an expected rotation of about in nine years.
The outer Southern spirals from IRDIS also show no apparent motion with respect to the detection by Boccaletti et al. (2013) (see Fig. 4d). Given the 5.2 years between the two epochs and a de-projected radius of 120 AU, the Keplerian motion of these spirals projected onto the disk plane is , which translates into 0.04″ at those radii. It is difficult to determine the exact precision of our comparison because of the ADI processing. However, a visual inspection of the relative position of the spirals yields no sign of such a shift, as the edges of these features (whose width is on average 0.06″) match very well between the two datasets.
The explanation for the absence of rotation is not obvious. If we focus on the inner spiral only (where the deviation from Keplerian motion is safe), one could try to relate its motion to that of a perturber object at larger radii. In fact, if this feature is excited by a massive object orbiting farther out (see discussion below), the spiral will be locked to the Keplerian motion of such an object. By conservatively assuming that the spiral has moved by less than in nine years, this scenario would point toward a perturber object at more than 70 AU. Alternatively, the observed feature may only resemble a spiral arm and instead be a geometrical property of the disk (such as the top of the inner wall) whose appearance in scattered light is, therefore, not supposed to change with time. Analogously, one may be tempted to explain the Southern spirals with the geometrical properties of the disk. If one gives credence to the Western spirals of Fig. 3, the global arm-pattern visible from the near side resembles the wing-structure observed at small radii, with multiple arms symmetric around the minor axis which arch toward the far side. In this scenario, the Northern spirals would be morphologically different from the Southern/Western spirals and possibly be in Keplerian rotation. Such an idea can be tested after 2020.
The interpretation of the multiple arms from the IRDIS small/medium scale (left/middle panels in Fig. 3) is also not straightforward. The perfect symmetry around the minor axis suggests that this pattern is connected to the dust scattering properties and to the ADI processing from an inclined disk. As discussed in Sect. 4.2, the morphology of the wing structure can be explained by the ADI processing. It is also nonetheless clear from Fig. 5 that sharp features as the multiple arms observed from IRDIS cannot be generated from the pipeline. Moreover, the spatial consistency between different pynpoint runs, different instruments and observing modes reinforce that these features are real, even though their morphology differs from the nominal spiral arms frequently observed from face-on disks (e.g., Garufi et al. 2013; Wagner et al. 2015). Second epoch observations of these inner features (as early as in 2018) and additional high-contrast ADI images as well as complementary RDI (Reference star Differential Imaging, Mawet et al. 2012) images of inclined disks will help unravel the nature of these arms.
A lively debate on the nature of the increasingly observed spirals in protoplanetary disks is ongoing. First of all, it has not been solved yet whether the detection of these features in scattered light reflects an intrinsic change in the dust distribution/properties down to the disk mid-plane (as proposed for SAO206462 by Pérez et al. 2014; Quanz 2015). Many authors instead favor a scenario where these are due to a change in the pressure scale height (e.g., Juhász et al. 2015; Pohl et al. 2015). Disentangling the spiral morphology is also fundamental to determine their causes. The effects of the gravitational instability (Durisen et al. 2007) on the disk morphology has been studied with hydrodynamical simulations by e.g., Dong et al. (2015a), which successfully reproduce the appearance of SAO206462 and MWC758. However, gravitational instability has been disfavored in many specific cases because of the insufficiently massive nature of these disks (e.g., Boccaletti et al. 2013; Garufi et al. 2013), even though the gas mass in these disks is not firmly constrained. Recently, planet-disk interactions have been investigated quantitatively (e.g. Juhász et al. 2015; Pohl et al. 2015; Dong et al. 2015b). Despite the good agreement with the observed brightness contrast of spirals, this hypothesis has difficulties in reproducing the observed aperture of spirals (unless planets at very large radii or very high disk scale heights are invoked). Thus, the current framework does not allow us to firmly ascribe the spiral structure of HD100546 to any scenario.
4.4 Disk interaction with planet b
Independent multi-epoch and -filter observations of HD100546 in thermal IR have revealed the existence of a point source at , sitting on top of an extended component (Quanz et al. 2013a, 2015; Currie et al. 2014, 2015). The current interpretation of the point source is thermal emission from a young planet (possibly surrounded by a circumplanetary disk) which is released from an effective area with and effective temperature K (Quanz et al. 2015). This claim implied the rejection of a scattered light nature for this emission. The absence of a localized brightness enhancement from our ZIMPOL data in correspondence of the planet supports the validity of this rejection. Nonetheless, the large-scale emission around the point source is spatially consistent with our diffuse polarized emission to North (see Fig. 4b), suggesting that our observed diffuse emission is (mostly) scattered light. A similar feature is detected by Quanz et al. (2015) and Currie et al. (2014) to the SE and it also matches the diffuse enhancement from ZIMPOL.
Our ADI K band image reveals an extended emission at the location of the b planet (see mid-bottom panel of Fig. 3). This emission has no equivalent in the H band, and such a discrepancy is the only significant difference between the two wavebands. This dissimilarity can be either explained by the different optical depth of incident light or by the different extinction experienced through the disk by a local emission. The former explanation implies that the emission is purely scattered light and that a change in the disk morphology occurs between the depths where photons with H and K wavelengths are scattered. However, these two surfaces are probably too close to allow such a change. Moreover, the lack of a similar discrepancy between the H and K band NACO images by Avenhaus et al. (2014a) casts further doubt on the hypothesis. The latter explanation lies on the assumption that a local source of heat is present at a certain disk depth. This thermal emission would be subject to a higher extinction at shorter wavelength and would thus explain the observed difference. The recovery of this emission might still be possible in the H band, as shown by Currie et al. (2015). In any case, it remains to be understood whether the K band detection is actually related to any process of planet formation. A more in depth analysis of this emission, complemented by the simultaneous IFS data, will be presented by Sissa et al. (in prep.).
From the entire dataset analyzed in this paper, any disk interaction with planet b remains fairly elusive. First of all, there is no hint of any disk gap at the planet location. Based on the local noise of our images, we can rule out with 3 confidence any brightness discontinuity more pronounced than . We also integrated the flux from the back side contained in an ellipse passing through b (obtained consistently with the disk geometry) but we could not infer any discontinuity. The absence of any brightness discontinuity does not necessarily imply the absence of gaps in the distribution of m-sized particles, since the high inclination of the source may not enable the observer to see through a narrow hole. From the angular resolution of the ZIMPOL data ( AU), some naive geometrical considerations, and the assumption of a 5 AU disk scale height at 50 AU (Montesinos et al. 2015), one obtains that our images could still detect a gap in the small dust grain distribution larger than AU. Furthermore, the observation of a significant portion of the disk wall at the outer gap edge may result in a significantly different flux distribution (as shown for the inner rim in Fig. 2a). The absence of an observed discontinuity is difficult to reconcile with the large gaseous gap expected to be sculpted by a giant planet in a few tens orbits timeframe (e.g., Crida et al. 2006; Masset 2008), unless the planet is very young ( years, given its orbital period). de Juan Ovelar et al. (2013) showed that small polarized light brightness discontinuities in the R band are expected from interactions with mass planets. An additional possibility to keep in mind is that the polarized signal deriving from the disk surface is blended with the contribution from halo particles. Speculatively, this contribution may be significant in highly inclined systems as suggested by the notion that sharp, prominent features like spirals and annular gaps are hardly observed in polarized light from inclined disks (see e.g., Hashimoto et al. 2012; Takami et al. 2013; Follette et al. 2015; Thalmann et al. 2015).
The only disk properties that can be ascribed to the interaction with b are the polarized light peak at comparable azimuthal angles (see Fig. 2d) and the Northern arm highlighted in Fig. 1e and Fig. 3. As discussed in Sect. 4.2, the large-scale enhancement of scattered light to the North is a possible imprint of a diffusely higher disk scale height. This may be in turn due to the hot environment produced by an enhanced disk accretion rate in correspondence of a luminous giant planet (as shown by the hydrodynamical simulations customized to HD100546 by Montesinos et al. 2015). Finally, the spatial connection between the northern arm detected in both ZIMPOL and IRDIS and the planet b is tantalizing. As discussed in Sect. 4.3, spiral arms are a natural outcome of the planet-disk interaction. However, this arm appears wrapped in the opposite direction of the other spirals and, as commented in Sect. 4.3, it is not clear yet to what extent this and other arm-like features in our datasets can be due to a combination of dust scattering properties and disk geometry.
5 Summary and conclusions
We present the first SPHERE observations of HD100546 in scattered light, obtained with the scientific sub-systems ZIMPOL (polarized R band) and IRDIS (H and K band). Complementary data from previous works enabled us to draw a comprehensive picture of the disk emission from 0.6 m to 4.8 m.
The ZIMPOL images in the visible resemble previous images in the near-IR, with the presence of an ellipsoidal cavity, of two bright lobes at the disk inner rim and of a tremendous brightness contrast between the disk near (dark) and far side (bright). The main findings from the analysis of these images are:
The cavity size in the visible is AU. This is consistent with the estimate from mid-IR interferometry by Panić et al. (2014). This finding implies that the CO emission associated to the planet candidate c (Brittain et al. 2013) is located within the disk. The marginal difference with the cavity size in the near-IR (13 AU) can be ascribed to the different disk opacity at those wavelengths. The large discrepancy with the cavity size at millimeter wavelengths (25 AU, Wright et al. 2015) is qualitatively consistent with the dust differentiation expected from the interaction with a (yet unseen) giant planet in the disk cavity. The amount of such a discrepancy requires a mass companion.
A luminous wedge stands out from the dark side and this is much brighter than seen in previous near-IR images. The cause of this difference may also clarify the nature of the global darkness. If it is a time-related difference, then the wedge might be a penumbra in a globally shadowed region. If it is a wavelength-related difference, then we may be tracing the forward peak of an otherwise backward peaking scattering phase function. We favor the latter scenario and entrust the answer to upcoming near-IR images and custom models.
The spiral arm at AU revealed by Avenhaus et al. (2014a) does not show any proper motion in a nine-years timescale, inconsistent with being in Keplerian motion. This could either mean that the spiral is locked with the orbital motion of a companion at AU or that it is an azimuthally symmetric feature misinterpreted as a spiral.
The effect of multiple scattering might be important. This is suggested by the morphology of the U images, quantitatively resembling the expectations for an optically thick, significantly inclined disk when this phenomenon is taking into account (Canovas et al. 2015).
The near-IR IRDIS images, analyzed in ADI with pynpoint, show a complex multiple-arm structure at all spatial scales. In the H band, a bright knot along one of these arms roughly matches the location of the putative detection of planet c by Currie et al. (2015). This knot is not present in the K band. Moreover, an extended emission in the K band is detected at the location of the b planet (Quanz et al. 2013a). Our conclusions on the analysis of these images are:
The double-wing structure symmetric around the minor axis can be generated from a continuous disk by the ADI processing. However, the multiple-arm morphology is consistent through different reductions and instruments and is therefore real. It is yet unclear whether these features are all spiral arms at the disk surface or marginal anisotropies accentuated by the scattering properties from an inclined disk. Their temporal evolution will provide fruitful insight into their nature.
The extended brightness associated to HD100546b is most likely a thermal emission originating deep in the disk. The non-detection in the H band (and from all polarized light datasets) is in fact difficult to reconcile with the scattered light scenario. An in-depth analysis of this detection will be described by Sissa et al. (in prep.).
The direct comparison between the ZIMPOL and the IRDIS images with previous works also provide some insight into the disk geometry and the scattering properties. The difference between the total and the polarized light in the H band can be mainly ascribed to the ADI processing. However, the different brightness distribution at the disk inner edge may indicate that particles at the inner rim are more prone to forward-scatter photons, contrarily to what is seen at larger radii. Furthermore, the overall similarity of the polarized light at 0.6 m and the total intensity at 4.8 m supports the idea that the latter emission is also (mainly) scattered light.
All things considered, the imprints of the (probable) giant planets around HD100546 on the disk morphology remain elusive. We found no strong evidence of any disk feature which might be due to the interaction with HD100546c. Some of them (i.e. an arm apparently originating along the major axis, the enhanced polarization of a lobe at the inner rim, and a spiral arm at 30 AU) are, most likely, fortuitous. On the other hand, the maximized polarized brightness on a large scale around HD100546b may be speculatively connected to the expected enhanced disk accretion rate in proximity of the planet. Furthermore, the spatial connection between an arm-like structure and the localized thermal emission in the K band may also reveal some yet unclear planet-disk interplay. The absence of a detectable gap in correspondence of planet b also raises unanswered questions on the processes of planet formation. A narrow ( AU) or a shallow (providing a flux decrease) gap, as well as a very young nature for b may reconcile with our observations.
These data are a showcase of the capability of the new generation high-contrast imager SPHERE. However, a lot remains to be understood about the enigmatic planet-forming disk of HD100546. The future confirmation (or rejection) of planet candidates and the recovery of disk features similar to those shown here from other suspected planet-forming disks (with SPHERE, ALMA, or GPI) will also help unravel the mysterious geometry of the disk around HD100546 and the mechanisms governing planet formation.
Acknowledgements.We thank Paola Pinilla and Sean Brittain for the fruitful scientific discussion, and the referee for a thorough review that helped improve the text. We are grateful to the ESO technical operators at the Paranal observatory for their valuable help during the observations. Part of this work has been carried out within the frame of the National Centre for Competence in Research PlanetS supported by the Swiss National Science Foundation. SPQ, HMS, and MRM acknowledge the financial support of the SNSF. HA acknowledges support from the Millennium Science Initiative (Chilean Ministry of Economy), through grant ”Nucleus P10-022-F” and financial support from FONDECYT grant 3150643. This research has made use of the SIMBAD database, operated at the CDS, Strasbourg, France. SPHERE is an instrument designed and built by a consortium consisting of IPAG (Grenoble, France), MPIA (Heidelberg, Germany), LAM (Marseille, France), LESIA (Paris, France), Laboratoire Lagrange (Nice, France), INAF - Osservatorio di Padova (Italy), Observatoire de Geneve (Switzerland), ETH Zurich (Switzerland), NOVA (Netherlands), ONERA (France) and ASTRON (Netherlands) in collaboration with ESO. SPHERE was funded by ESO, with additional contributions from CNRS (France), MPIA (Germany), INAF (Italy), FINES (Switzerland) and NOVA (Netherlands). SPHERE also received funding from the European Commission Sixth and Seventh Framework Programmes as part of the Optical Infrared Coordination Network for Astronomy (OPTICON) under grant number RII3-Ct-2004-001566 for FP6 (2004-2008), grant number 226604 for FP7 (2009-2012) and grant number 312430 for FP7 (2013-2016).
- Amara & Quanz (2012) Amara, A. & Quanz, S. P. 2012, MNRAS, 427, 948
- Amara et al. (2015) Amara, A., Quanz, S. P., & Akeret, J. 2015, Astronomy and Computing, 10, 107
- Ardila et al. (2007) Ardila, D. R., Golimowski, D. A., Krist, J. E., et al. 2007, ApJ, 665, 512
- Augereau et al. (2001) Augereau, J. C., Lagrange, A. M., Mouillet, D., & Ménard, F. 2001, A&A, 365, 78
- Avenhaus et al. (2014a) Avenhaus, H., Quanz, S. P., Meyer, M. R., et al. 2014a, ApJ, 790, 56
- Avenhaus et al. (2014b) Avenhaus, H., Quanz, S. P., Schmid, H. M., et al. 2014b, ApJ, 781, 87
- Backman & Paresce (1993) Backman, D. E. & Paresce, F. 1993, in Protostars and Planets III, ed. E. H. Levy & J. I. Lunine, 1253–1304
- Bastien & Menard (1988) Bastien, P. & Menard, F. 1988, ApJ, 326, 334
- Benisty et al. (2015) Benisty, M., Juhasz, A., Boccaletti, A., et al. 2015, A&A, 578, L6
- Benisty et al. (2010) Benisty, M., Tatulli, E., Ménard, F., & Swain, M. R. 2010, A&A, 511, A75
- Beuzit et al. (2008) Beuzit, J.-L., Feldt, M., Dohlen, K., et al. 2008, in Society of Photo-Optical Instrumentation Engineers (SPIE) Conference Series, Vol. 7014, Society of Photo-Optical Instrumentation Engineers (SPIE) Conference Series, 18
- Biller et al. (2014) Biller, B. A., Males, J., Rodigas, T., et al. 2014, ApJ, 792, L22
- Boccaletti et al. (2013) Boccaletti, A., Pantin, E., Lagrange, A.-M., et al. 2013, A&A, 560, A20
- Brauer et al. (2008) Brauer, F., Henning, T., & Dullemond, C. P. 2008, A&A, 487, L1
- Brittain et al. (2014) Brittain, S. D., Carr, J. S., Najita, J. R., Quanz, S. P., & Meyer, M. R. 2014, ApJ, 791, 136
- Brittain et al. (2013) Brittain, S. D., Najita, J. R., Carr, J. S., et al. 2013, ApJ, 767, 159
- Canovas et al. (2015) Canovas, H., Ménard, F., de Boer, J., et al. 2015, A&A, 582, L7
- Claudi et al. (2008) Claudi, R. U., Turatto, M., Gratton, R. G., et al. 2008, in Society of Photo-Optical Instrumentation Engineers (SPIE) Conference Series, Vol. 7014, Society of Photo-Optical Instrumentation Engineers (SPIE) Conference Series, 3
- Crida et al. (2006) Crida, A., Morbidelli, A., & Masset, F. 2006, Icarus, 181, 587
- Currie et al. (2015) Currie, T., Cloutier, R., Brittain, S., et al. 2015, ApJ, 814, L27
- Currie et al. (2014) Currie, T., Muto, T., Kudo, T., et al. 2014, ApJ, 796, L30
- de Juan Ovelar et al. (2013) de Juan Ovelar, M., Min, M., Dominik, C., et al. 2013, A&A, 560, A111
- Debes et al. (2013) Debes, J. H., Jang-Condell, H., Weinberger, A. J., Roberge, A., & Schneider, G. 2013, ApJ, 771, 45
- Dohlen et al. (2008) Dohlen, K., Langlois, M., Saisse, M., et al. 2008, in Society of Photo-Optical Instrumentation Engineers (SPIE) Conference Series, Vol. 7014, Society of Photo-Optical Instrumentation Engineers (SPIE) Conference Series, 3
- Dong et al. (2015a) Dong, R., Hall, C., Rice, K., & Chiang, E. 2015a, ApJ, 812, L32
- Dong et al. (2012) Dong, R., Rafikov, R., Zhu, Z., et al. 2012, ApJ, 750, 161
- Dong et al. (2015b) Dong, R., Zhu, Z., Rafikov, R. R., & Stone, J. M. 2015b, ApJ, 809, L5
- Dullemond et al. (2001) Dullemond, C. P., Dominik, C., & Natta, A. 2001, ApJ, 560, 957
- Durisen et al. (2007) Durisen, R. H., Boss, A. P., Mayer, L., et al. 2007, Protostars and Planets V, 607
- Fedele et al. (2015) Fedele, D., Bruderer, S., van den Ancker, M. E., & Pascucci, I. 2015, ApJ, 800, 23
- Follette et al. (2015) Follette, K. B., Grady, C. A., Swearingen, J. R., et al. 2015, ApJ, 798, 132
- Fusco et al. (2006) Fusco, T., Petit, C., Rousset, G., et al. 2006, in Society of Photo-Optical Instrumentation Engineers (SPIE) Conference Series, Vol. 6272, Society of Photo-Optical Instrumentation Engineers (SPIE) Conference Series, 0
- Garufi et al. (2013) Garufi, A., Quanz, S. P., Avenhaus, H., et al. 2013, A&A, 560, A105
- Garufi et al. (2014) Garufi, A., Quanz, S. P., Schmid, H. M., et al. 2014, A&A, 568, A40
- Grady et al. (2005) Grady, C. A., Woodgate, B., Heap, S. R., et al. 2005, ApJ, 620, 470
- Hashimoto et al. (2012) Hashimoto, J., Dong, R., Kudo, T., et al. 2012, ApJ, 758, L19
- Hashimoto et al. (2011) Hashimoto, J., Tamura, M., Muto, T., et al. 2011, ApJ, 729, L17
- Henyey & Greenstein (1941) Henyey, L. G. & Greenstein, J. L. 1941, ApJ, 93, 70
- Janson et al. (2016) Janson, M., Thalmann, C., Boccaletti, A., et al. 2016, ApJ, 816, L1
- Juhász et al. (2015) Juhász, A., Benisty, M., Pohl, A., et al. 2015, MNRAS, 451, 1147
- Kraus & Ireland (2012) Kraus, A. L. & Ireland, M. J. 2012, ApJ, 745, 5
- Lafrenière et al. (2007) Lafrenière, D., Marois, C., Doyon, R., Nadeau, D., & Artigau, É. 2007, ApJ, 660, 770
- Levenhagen & Leister (2006) Levenhagen, R. S. & Leister, N. V. 2006, MNRAS, 371, 252
- Liu et al. (2003) Liu, W. M., Hinz, P. M., Meyer, M. R., et al. 2003, ApJ, 598, L111
- Macintosh et al. (2014) Macintosh, B., Graham, J. R., Ingraham, P., et al. 2014, Proceedings of the National Academy of Science, 111, 12661
- Maire et al. (2015) Maire, A.-L., Bonnefoy, M., Ginski, C., et al. 2015, ArXiv e-prints [\eprint[arXiv]1511.04072]
- Marino et al. (2015) Marino, S., Perez, S., & Casassus, S. 2015, ApJ, 798, L44
- Marois et al. (2006) Marois, C., Lafrenière, D., Doyon, R., Macintosh, B., & Nadeau, D. 2006, ApJ, 641, 556
- Masset (2008) Masset, F. S. 2008, in EAS Publications Series, Vol. 29, EAS Publications Series, ed. M.-J. Goupil & J.-P. Zahn, 165–244
- Mawet et al. (2012) Mawet, D., Pueyo, L., Lawson, P., et al. 2012, in Society of Photo-Optical Instrumentation Engineers (SPIE) Conference Series, Vol. 8442, Society of Photo-Optical Instrumentation Engineers (SPIE) Conference Series, 4
- Min et al. (2012) Min, M., Canovas, H., Mulders, G. D., & Keller, C. U. 2012, A&A, 537, A75
- Min et al. (2016) Min, M., Rab, C., Woitke, P., Dominik, C., & Ménard, F. 2016, A&A, 585, A13
- Montesinos et al. (2015) Montesinos, M., Cuadra, J., Perez, S., Baruteau, C., & Casassus, S. 2015, ApJ, 806, 253
- Mulders et al. (2013a) Mulders, G. D., Min, M., Dominik, C., Debes, J. H., & Schneider, G. 2013a, A&A, 549, A112
- Mulders et al. (2013b) Mulders, G. D., Paardekooper, S.-J., Panić, O., et al. 2013b, A&A, 557, A68
- Murakawa (2010) Murakawa, K. 2010, A&A, 518, A63
- Muto et al. (2012) Muto, T., Grady, C. A., Hashimoto, J., et al. 2012, ApJ, 748, L22
- Panić et al. (2014) Panić, O., Ratzka, T., Mulders, G. D., et al. 2014, A&A, 562, A101
- Pantin et al. (2000) Pantin, E., Waelkens, C., & Lagage, P. O. 2000, A&A, 361, L9
- Pavlov et al. (2008) Pavlov, A., Feldt, M., & Henning, T. 2008, in Astronomical Society of the Pacific Conference Series, Vol. 394, Astronomical Data Analysis Software and Systems XVII, ed. R. W. Argyle, P. S. Bunclark, & J. R. Lewis, 581
- Pérez et al. (2014) Pérez, L. M., Isella, A., Carpenter, J. M., & Chandler, C. J. 2014, ApJ, 783, L13
- Pineda et al. (2014) Pineda, J. E., Quanz, S. P., Meru, F., et al. 2014, ApJ, 788, L34
- Pinilla et al. (2012) Pinilla, P., Benisty, M., & Birnstiel, T. 2012, A&A, 545, A81
- Pinilla et al. (2015) Pinilla, P., Birnstiel, T., & Walsh, C. 2015, A&A, 580, A105
- Pohl et al. (2015) Pohl, A., Pinilla, P., Benisty, M., et al. 2015, MNRAS, 453, 1768
- Quanz (2015) Quanz, S. P. 2015, Ap&SS, 357, 148
- Quanz et al. (2015) Quanz, S. P., Amara, A., Meyer, M. R., et al. 2015, ApJ, 807, 64
- Quanz et al. (2013a) Quanz, S. P., Amara, A., Meyer, M. R., et al. 2013a, ApJ, 766, L1
- Quanz et al. (2013b) Quanz, S. P., Avenhaus, H., Buenzli, E., et al. 2013b, ApJ, 766, L2
- Quanz et al. (2011) Quanz, S. P., Schmid, H. M., Geissler, K., et al. 2011, ApJ, 738, 23
- Rapson et al. (2015) Rapson, V. A., Kastner, J. H., Andrews, S. M., et al. 2015, ApJ, 803, L10
- Reggiani et al. (2014) Reggiani, M., Quanz, S. P., Meyer, M. R., et al. 2014, ApJ, 792, L23
- Rice et al. (2006) Rice, W. K. M., Armitage, P. J., Wood, K., & Lodato, G. 2006, MNRAS, 373, 1619
- Sallum et al. (2015) Sallum, S., Follette, K. B., Eisner, J. A., et al. 2015, Nature, 527, 342
- Takami et al. (2013) Takami, M., Karr, J. L., Hashimoto, J., et al. 2013, ApJ, 772, 145
- Tatulli et al. (2011) Tatulli, E., Benisty, M., Ménard, F., et al. 2011, A&A, 531, A1
- Thalmann et al. (2015) Thalmann, C., Mulders, G. D., Janson, M., et al. 2015, ApJ, 808, L41
- Thalmann et al. (2008) Thalmann, C., Schmid, H. M., Boccaletti, A., et al. 2008, in Society of Photo-Optical Instrumentation Engineers (SPIE) Conference Series, Vol. 7014, Society of Photo-Optical Instrumentation Engineers (SPIE) Conference Series, 3
- Tinbergen (2005) Tinbergen, J. 2005, Astronomical Polarimetry
- Tsukagoshi et al. (2014) Tsukagoshi, T., Momose, M., Hashimoto, J., et al. 2014, ApJ, 783, 90
- van den Ancker et al. (1997) van den Ancker, M. E., The, P. S., Tjin A Djie, H. R. E., et al. 1997, A&A, 324, L33
- van Leeuwen et al. (2007) van Leeuwen, F., Feast, M. W., Whitelock, P. A., & Laney, C. D. 2007, MNRAS, 379, 723
- Vigan et al. (2010) Vigan, A., Moutou, C., Langlois, M., et al. 2010, MNRAS, 407, 71
- Wagner et al. (2015) Wagner, K., Apai, D., Kasper, M., & Robberto, M. 2015, ApJ, 813, L2
- Walsh et al. (2014) Walsh, C., Juhász, A., Pinilla, P., et al. 2014, ApJ, 791, L6
- Wright et al. (2015) Wright, C. M., Maddison, S. T., Wilner, D. J., et al. 2015, MNRAS, 453, 414