Press Release
Open Access
Issue
A&A
Volume 637, May 2020
Article Number L5
Number of page(s) 6
Section Letters to the Editor
DOI https://doi.org/10.1051/0004-6361/202038008
Published online 20 May 2020

© A. Boccaletti et al. 2020

Licence Creative CommonsOpen Access article, published by EDP Sciences, under the terms of the Creative Commons Attribution License (https://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

1. Introduction

Planets are believed to form in protoplanetary disks over a few million years. Images at all wavelengths exhibit a variety of structures, such as asymmetries, clumps, rings, or spirals. Most of these features, which hypothetically betray the presence of planets via gravitational or hydrodynamical effects, have been observed in disks orbiting Herbig Ae/Be stars because of their suitable brightness, particularly in the near-infrared (NIR) (see for instance Grady et al. 2001; Clampin et al. 2003; Fukagawa et al. 2006; Isella et al. 2010; Christiaens et al. 2014). Planet formation is generally favored (Zhu 2015) over instabilities in a massive self-gravitating disk (Rice et al. 2005). In the early stage of planet formation, hydrodynamical simulations indicate that the accretion process generates at the planet location an inner and outer spiral pattern due to Lindblad resonances induced by disk-planet interactions (Gressel et al. 2013). While this crucial step is well-documented by theoretical works (Dong et al. 2015; Bae & Zhu 2018), observational evidence is rare and not fully conclusive. Recently, Muro-Arena et al. (2020) identified a kinked spiral in the protoplanetary disk SR 21, suggesting that it may be a promising system for potentially witnessing ongoing planet formation.

With a spectral type A0 (2.4 ± 0.2 M, DeWarf et al. 2003), AB Aurigae is one of the closest (d = 162.9 ± 1.5 pc, Gaia Collaboration 2018), and one of the most intensively studied Herbig Ae star. Optical scattered light images of its protoplanetary disk have been obtained in the visible (Grady et al. 1999; Fukagawa et al. 2004), and in the NIR (Perrin et al. 2009; Hashimoto et al. 2011), revealing a moderately inclined (i ∼ 30°), flared disk and prominent spiral patterns extending from 200 to 450 au. In fact, AB Aurigae exhibits among the most spectacular spirals imaged so far in scattered light and high-contrast polarimetric imaging.

Millimetric observations have also been key in constraining the distribution of gas and dust in the AB Aur system. Using the IRAM Plateau de Bure Interferometer, Piétu et al. (2005) reported the presence of a large CO and dust rotating disk with a central cavity of an inner radius of ∼70 au. Surprisingly, Tang et al. (2012) found counter-rotating CO spirals in the outer disk, which has been explained by projection effects of accretion flows arising from above the disk’s midplane. The high accretion rate (∼10−7M yr−1, Salyk et al. 2013) measured for this 1 Myr object reinforces this interpretation.

More recently, Tang et al. (2017) studied the CO 2-1 and dust continuum distribution at very high angular scale (0.08″ or 13 au) using ALMA and revealed two spectacular CO spirals residing inside the large disk cavity. The authors speculated that these spirals are driven by planet formation and could be generated by two planets respectively located at ∼30 and 60–80 au from the star. Motivated by these results, we conducted Very Large Telescope (VLT)/SPHERE observations of the AB Aurigae system and present them in this Letter. Here, we describe the observations, analyze the data and discuss the possible origin of the inner spiral pattern.

2. Observations

AB Aur (V = 7.05, H = 5.06, K = 4.23) has been observed with SPHERE (Beuzit et al. 2019) in polarimetry (Dec. 2019) and in spectro-photometry (Jan. 2020). The first data set was obtained with the infrared camera IRDIS (Dohlen et al. 2008; de Boer et al. 2020) using Pupil Tracking (PT) in the H band, which allows us to retrieve, at the same time, the polarized intensity of the disk and the total intensity by making use of differential polarimetric imaging (DPI) and of angular differential imaging (ADI). The second epoch used the IRDIFS-EXT mode of SPHERE, which combines IRDIS in the K1K2 filters, and the spectrograph (IFS, Claudi et al. 2008) in low-resolution mode YJH (R = 30), also processed with ADI. The observing log is provided in Table A.1.

Both DPI and ADI increase the contrast by rejecting the starlight, exploiting either the fact that the stellar light is unpolarized or that off-axis objects rotate with the field in PT mode. While DPI provides unbiased polarimetric images of disks, ADI induces very strong photometric and astrometric artifacts owing to the self-subtraction effect (Milli et al. 2012). The data were processed at the SPHERE Data Center following Pavlov et al. (2008), Delorme et al. (2017), Maire et al. (2016). The calibrated data cubes were then processed with SpeCal for total intensity data (Galicher et al. 2018) and with a custom pipeline for polarimetric data (as in Bhowmik et al. 2019). The total intensity data were reduced with Principal Component Analysis (Soummer et al. 2012), whereas for polarimetry, we calculated the azimuthal Stokes parameters Qϕ and Uϕ maps following Schmid et al. (2006). For a sanity check, we compared our polarimetric reductions with those of the IRDAP pipeline (van Holstein et al. 2020), which implements a careful treatment of the instrumental polarization. We found that since the instrumental effects are much fainter than the disk structures we further discuss below, the Qϕ images are identical in both pipelines. The first epoch is significantly better in terms of starlight rejection, therefore the following sections are mainly based on these observations.

3. Disk morphology

As the disk is moderately inclined, the ADI process strongly and irremediably attenuates the disk signal (Fig. 1c). Therefore, the following description is based on the polarimetric H-band Qϕ image (Figs. 1a, b). The disk is visible all across the IRDIS 11″ × 11″ field of view. We identify three main regions: the dust ring, the outer spirals, and the inner spirals.

thumbnail Fig. 1.

Images of the AB Aur system obtained with SPHERE in polarized light (a, b) and unpolarized light (c). A large field of view (10″) is shown in panel a, where the polarized intensity has been multiplied with the square of the stellocentric distance (Qϕ × r2) to visually enhance the outer part of the disk. A narrower field of view (2″) of the Qϕ map is displayed in panel b, in which the inner spirals are labelled S1 and S2. Comparable ADI-processed image is shown in panel c for total intensity in the H band filter. Features f1 and f2 are discussed in Sect. 4. North is up, east is left.

The position of the ring detected with ALMA in the continuum (Tang et al. 2017) is shown in Fig. 2a. In the sub-mm, it encompasses an elliptical region oriented at about 60° north to east, from a minimal radius of 0.77″ to a maximal separation of 1.19″, corresponding to a ring of ∼140 au on average (yellow ellipses in Fig. 2a). The aspect ratio is consistent with an inclination of 30°. In scattered light, the ring is not homogeneous with brighter parts in the west and the northeast, as opposed to the smooth 1.3 mm image, which instead peaks to the southwest. SPHERE resolves some structures in the ring that were already reported in scattered light by Oppenheimer et al. (2008) and Hashimoto et al. (2011). In particular, a depleted region in the polarized flux, formerly reported by Oppenheimer et al. (2008), is also observed at PA = 343° although the contrast is low and its visibility is mostly a matter of intensity cut. We do not detect any point-source at this location, confirming the findings of Perrin et al. (2009).

thumbnail Fig. 2.

Full-field polarimetric image (same as in Fig. 1, Qϕ × r2) displayed in panel a, on top of which several structures previously detected are overlaid: Hashimoto et al. (2011) spirals (orange lines), out of plane accreting spirals from Tang et al. (2012; blue lines), and the dust ring (yellow ellipses) from the 0.9 mm ALMA image (Tang et al. 2017). The inner spiral arms (S1_CO, S2_CO) detected by Tang et al. (2017) are superimposed to the SPHERE polarimetric image Qϕ, further processed with a gradient algorithm to enhance the inner spirals (b), and the ALMA CO map (c), in a 2″ field of view. Features f1 and f2 are discussed in Sect. 4. North is up, east is left.

Outward of the ring, there are several spiral-like structures which, given the sensitivity achieved with SPHERE DPI, could in fact be the trace of one single grand spiral with several ramifications. This pattern likely starts from the south and rolls clockwise while it gets broader with manifold branches (Fig. 1a). It is particularly pronounced at the east and matches well to all the spirals identified in 2009 H-band data by Hashimoto et al. (2011) at Subaru (orange lines in Fig. 2a), indicating no perceptible motion in ten years, as already concluded by Lomax et al. (2016) on a shorter timescale. Therefore, this pattern is unlikely to be driven by planets inside the ring. Out-of-plane CO spirals (Tang et al. 2012), which have no obvious counterpart in scattered light, are indicated in cyan (Fig. 2a).

Inward of the ring (< 0.6″), we detect two main spirals (S1 and S2 in Fig. 1b), starting as close as the edge of the coronagraphic mask (about 95 mas in radius) and in good agreement with the CO spirals identified by Tang et al. (2017, and shown as green lines labeled S1_CO and S2_CO in Figs. 2b, c). The angular resolution brought by SPHERE evidences that the two inner spirals seem to intersect at the south of the star exactly at the location of an intensity peak (ρ = 0.21″, PA = 177.7°), which itself is found at a very similar position than the CO peak. S2 matches well the western CO spiral, while S1 seems to deviate from the eastern CO spiral. We suspect that this mismatch is due to the difference of angular resolution, which for ALMA is about 0.11″, while SPHERE provides in H band a resolution of 0.04″. As a consequence, it is difficult to disentangle S1 from S2 at the north in the CO map. S2 is more diffuse than S1 and appears larger in the west. High-pass filtering of the polarimetric image (Fig. 2b) shows higher complexity with clumps and several arms.

While the CO map contains a cavity inward of the CO spirals the SPHERE observations reveal a third component (PA ∼ 270° in Fig. 1b) connecting the center of the image to S2, as a kind of “bridge” of 25 − 30 au long. This is not a post-processing artifact since this feature is seen both in the polarimetry image and in the total intensity ADI image, and was in fact already distinguishable in Hashimoto et al. (2011), but also in the gas line HCO+ (Rivière-Marichalar et al. 2019). Whether it is a true connection or a projection effect cannot be determined with these observations alone. On each side of the “bridge”, the polarized intensity is rather low, even lower than in between the spiral arms. Interestingly, in the direction of this feature the intensity of the spiral S2 and of the ring is fainter, possibly tracing a shadow pattern. If so, it will favor a feature near the midplane of the disk. Finally, S1 shows an increase of its width at PA ∼ 115° as a compact clump that is ∼0.12″ long (Figs. 1b and 2b), complementary images are displayed in Appendix A.

4. Point-like sources

Many structures are observed in the disk but two features have a strong interest. The first one (f1) is the strongest signal in Figs. 1b and c, south of the coronagraphic mask edge, and a second one (f2) is a nearly point source further out to the north. Interestingly, both show up in DPI and ADI images (Fig. 3). They are also detected at the second epoch in total intensity (ADI), but with a much lower signal to noise ratio (Fig. A.1).

thumbnail Fig. 3.

Positions of features f1 and f2 overlaid on the polarimetric image shown for two intensity thresholds. The dusty ring is materialized with the dashed ellipses. For visualization purposes, the polarized intensity has been high-pass filtered and multiplied with the stellocentric distance (Qϕ × r).

We extracted the photometry and astrometry with SpeCal. The positions of the sources are based on the ADI reduction using the PCA 5 modes algorithm (Fig. 2b) and the negative fake planet injection method (for details see Galicher et al. 2018), but the values are consistent within error bars for all algorithms. f1 is located at ρ = 0.160 ± 0.006″, PA = 203.9 ± 2.9°. It coincides with the root of the eastern spiral identified as S1 in the ALMA image, where the two spirals S1 and S2 observed with SPHERE appear to intersect. The shape of this feature in the Qϕ image (Fig. 3a) is remarkable and we discuss it further in the next section. The second source, f2, lies at the very edge of the dust ring inside the cavity, in the prolongation of the eastern spiral. Its position is ρ = 0.681 ± 0.006″, PA = 7.6 ± 1.8°.

Feature f1 is elongated while f2 is nearly point-like. This characteristic is confirmed by the photometric analysis. The contrast relative to the star derived from IRDIS data in the H band and measured at f1 is quite variable from one ADI reduction to another, ranging from 6.5 ± 0.8 × 10−4 to 6.8 ± 0.6 × 10−5, depending on the effectiveness (ability to suppress the starlight) of the ADI algorithm.

Similarly, the photometry in the K1K2 IRDIS image obtained at the second epoch varies from one ADI reduction to the other, but with a small difference of contrast between the two narrow band filters: 3.4 ± 0.5 × 10−4 and 4.3 ± 0.5 × 10−4, respectively. Whether it is indicative of a red color remains to be confirmed with higher quality data in the K band. Overall, we cannot consider the total intensity contrast measurements of f1 reliable as it cannot be assimilated to a point source. On the other hand, f2 has a stable contrast in the H band of 4.9 ± 0.5 × 10−6, and 4.8 ± 0.8 × 10−6 in K1K2, indicative of a bona-fide point source, but with a grey color. Photometric measurements with the IFS (not shown here) in the YJH bands are not much more conclusive regarding f1, while f2 is actually not detected.

The fact that f1 and f2 are detectable in polarimetry and total intensity strongly suggests that these signals cannot result from the photosphere of some protoplanets. Instead, they are likely to be the telltale of dust particles detected from scattering or from other thermally induced mechanisms.

5. Spirals

The presence of spiral arms in protoplanetary disks has been usually associated to the indirect clues of planets provided that the disks are not massive enough to undergo self-gravity (Muto et al. 2012; Boccaletti et al. 2013; Benisty et al. 2015, for instance). This link with a gravitational perturber has never been definitely established, except for the prominent spiral arms in HD 100453, very likely launched by a stellar M-dwarf companion (Dong et al. 2016). In that respect, the ALMA observations of AB Aur (Tang et al. 2017) are likely the only case, together with GG Tau A (Phuong et al. 2020), to provide a strong suspicion for the presence of a protoplanet associated to spiral arms.

To model the morphology of the inner spirals we rely on the formulation of Muto et al. (2012) described in Appendix B. A direct fit of this analytical relation (Eq. (B.1)) onto the spiral positions was found too unstable given the difficulty to isolate the spine of the spiral arms against the other disk structures. Instead, we proceed with a manual adjustment of the theoretical spiral location taking into account the disk inclination (i = 30°), the disk position angle (PA = 60°), the gas disk vertical scale-height (H0 = 25 au at 100 au) and the flaring index (r1.2), following di Folco et al. (2009), Hashimoto et al. (2011), Tang et al. (2012), Li et al. (2016). Because the spirals are pressure bumps and the NIR emission is optically thick, the net effect is a local variation of the disk thickness, hence these spirals materialize at the disk surface and not in the mid-plane (Juhasz et al. 2015). Therefore, we assume that the altitude of the disk layer at which the optical depth (τν) is equal to one is H(τ = 1) = 1.5 × H0. These values are used to project the analytical relation θ(r) onto the plane of the sky and are known to some level of accuracy. Therefore, the solutions presented here provide qualitative results, which may not be unique. We note, in particular, a degeneracy between the inclination, the disk scale-height, and the radial distance of the putative planet, whereas the flaring index has only a weak impact on the shape of the projected spiral.

To account for S1 we considered two cases for which we forced a planet to be located close to the positions of the features discussed in Sect. 4. A planet located at rc = 0.184″, θ0 = 196.8° (in the orbital plane of the disk) reproduces the shape of the spiral and is also in accordance with the location of f1 as seen in Figs. 4a, and b (green line). The analytical model of Eq. (B.1), while simple, matches the “S” shape expected for a spiral density wave surprisingly well as it goes on to trigger the accretion of gas onto a protoplanet. On the contrary, imposing the planet location at f2 matches only the external part of the spiral S1 (blue line in Fig. 4a). For comparison, the contours of the CO peak in the ALMA image is shown in Fig. 4b. Being localized at the surface, the spiral twists appear slightly offset from the perturber positions (circles), which are in the midplane.

thumbnail Fig. 4.

Same data as Fig. 3 with three different models of spirals matching S1 (green and blue lines) and S2 (purple line). The positions of the perturbers are indicated with a circle (panel a). Panel b: zoom in version of panel a, where only the best model reproducing the spiral twist in S1 is shown. Panel c: the contours of the CO map at the location of the candidate protoplanet are overlaid, together with the spiral model offset by 14.1° (dashed green line), corresponding to the 4 yr timelapse between ALMA and SPHERE observations.

For the western spiral S2, since there is no obvious features nor strong signal which we can consider as a source, we blindly searched for the best matching model. The purple line in Fig. 4a corresponds to rc = 0.172″, θ0 = 6.0°. This model reproduces well the external part of the western spiral, which presents a more diffuse aspect or broader shape than the eastern one, but it fails to match the innermost, eastern arm which is less tightly wound as it approaches the occulting spot (0.1″ <  r <  0.16″).

6. Discussions

We provide a further exploratory analysis of the two identified objects, noting that they have both been detected in polarimetry and total intensity.

The object labeled f1 is fitting the source of the S1 spiral which at this location features an “S” shape, suggestive of the stationary pattern composed of inward and outward spiral density waves excited by the gravitational potential of an accreting protoplanet (see Kley & Nelson 2012, for a review). The spirals are trailing, as predicted by theory, since the near side is the southeast part of the disk (Fukagawa et al. 2004) and the resolved velocity field (Piétu et al. 2005) implies that the disk rotates counterclockwise. However, it is not a straightforward exercise to determine whether the signal is from pure scattering (polarized) or a mix of scattering and emission (usually unpolarized). Zhu (2015) has calculated the spectral energy distribution emitted from an accreting circumplanetary disk, which can be compared to our photometric measurements, assuming the signal of f1 is dominated by emission. The luminosity varies as the product of the planet mass with the accretion rate, and also depends on the inner radius of the circumplanetary disk. These three parameters necessarily lead to a high degree of degeneracy. Although it is difficult to overcome the ADI bias in total intensity measured for the particular case of f1, as it is not point-like, we measured a minimum contrast of 6.5 × 10−4 in the H band. If we assume a moderate accretion rate of 10−8M yr−1 as in Zhu (2015), it translates into a planet mass of 4 − 13MJ. We note that the accretion rates estimated for PDS 70b are orders of magnitude lower (Haffert et al. 2019; Hashimoto et al. 2020), which then would imply stellar mass for f1. As a consequence, the flux of f1 is likely not to be dominated by thermal emission.

On the contrary, the object f2 is not directly associated to a spiral, so we assume its flux can be attributed to pure scattering. Mass constraints can be derived from dynamical argument. Indeed, the inner edge of the disk cavity at ∼140 au can be sculpted by the chaotic zone of a planetary mass object. Assuming a circular orbit for the sake of simplicity, Wisdom (1980) provides a relation between the planet/star mass ratio, and the inner edge/planet distance ratio, from which we derived a mass of ∼3 MJ. The AB Aur cavity is very large, such that, if resulting from dynamical interplay with planets, several bodies including f1 should be involved.

Finally, the object f1 appears rotated counterclockwise from the CO peak in the ALMA image which, considering the stellar mass, the 4 yr interval, and the physical stellocentric distance of f1, is in excellent agreement with the expected rotation (14.1 ± 1.1° in the orbital plane) of the spiral due to the Keplerian motion of the protoplanet candidate (Fig. 4c). However, the peak of the CO emission is at a radial distance slightly larger than expected. There are a few technical reasons accounting for this difference, which remain to be investigated thoroughly. Nevertheless, to our knowledge, this is the first time that we can confidently associate the rotation of a spiral pattern with the orbital motion of a planet candidate, which reinforces our confidence in the nature of the detected spiral twist.

In conclusion, the SPHERE observations of AB Aur in scattered light combined to the ALMA data in the thermal regime provide strong evidence that we are actually witnessing ongoing planet formation revealed by its associated spiral arms. Further observations would be required to confirm this result and to derive better mass estimates for potential planets in this location.

Acknowledgments

French co-authors acknowledge financial support from the Programme National de Planétologie (PNP) and the Programme National de Physique Stellaire (PNPS) of CNRS-INSU in France. This work has made use of the SPHERE Data Centre, jointly operated by OSUG/IPAG (Grenoble), PYTHEAS/LAM/CESAM (Marseille), OCA/Lagrange (Nice) and Observatoire de Paris/LESIA (Paris). We thank P. Delorme (SPHERE Data Centre) for his efficient help during the data reduction process. YWT acknowledges support through MoST grant 108-2112-M-001-004-MY2. ALM acknowledges the financial support of the F.R.S-FNRS through a postdoctoral researcher grant.

References

  1. Bae, J., & Zhu, Z. 2018, ApJ, 859, 118 [NASA ADS] [CrossRef] [Google Scholar]
  2. Benisty, M., Juhasz, A., Boccaletti, A., et al. 2015, A&A, 578, L6 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  3. Beuzit, J. L., Vigan, A., Mouillet, D., et al. 2019, A&A, 631, A155 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  4. Bhowmik, T., Boccaletti, A., Thébault, P., et al. 2019, A&A, 630, A85 [Google Scholar]
  5. Boccaletti, A., Pantin, E., Lagrange, A. M., et al. 2013, A&A, 560, A20 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  6. Christiaens, V., Casassus, S., Perez, S., van der Plas, G., & Ménard, F. 2014, ApJ, 785, L12 [Google Scholar]
  7. Clampin, M., Krist, J. E., Ardila, D. R., et al. 2003, AJ, 126, 385 [NASA ADS] [CrossRef] [Google Scholar]
  8. Claudi, R., Beuzit, J. L., Feldt, M., et al. 2008, Eur. Planet. Sci. Congr., 2008, 875 [Google Scholar]
  9. de Boer, J., Langlois, M., van Holstein, R. G., et al. 2020, A&A, 633, A63 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  10. Delorme, P., Meunier, N., Albert, D., et al. 2017, Proc. SF2A, 347 [Google Scholar]
  11. DeWarf, L. E., Sepinsky, J. F., Guinan, E. F., Ribas, I., & Nadalin, I. 2003, ApJ, 590, 357 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  12. di Folco, E., Dutrey, A., Chesneau, O., et al. 2009, A&A, 500, 1065 [Google Scholar]
  13. Dohlen, K., Langlois, M., Saisse, M., et al. 2008, in Proc. SPIE, eds. I. S. McLean, & M. M. Casali, 118 [Google Scholar]
  14. Dong, R., Zhu, Z., Rafikov, R. R., & Stone, J. M. 2015, ApJ, 809, L5 [NASA ADS] [CrossRef] [Google Scholar]
  15. Dong, R., Zhu, Z., Fung, J., et al. 2016, ApJ, 816, L12 [NASA ADS] [CrossRef] [Google Scholar]
  16. Fukagawa, M., Hayashi, M., Tamura, M., et al. 2004, ApJ, 605, L53 [NASA ADS] [CrossRef] [Google Scholar]
  17. Fukagawa, M., Tamura, M., Itoh, Y., et al. 2006, ApJ, 636, L153 [NASA ADS] [CrossRef] [Google Scholar]
  18. Gaia Collaboration (Brown, A. G. A., et al.) 2018, A&A, 616, A1 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  19. Galicher, R., Boccaletti, A., Mesa, D., et al. 2018, A&A, 615, A92 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  20. Grady, C. A., Woodgate, B., Bruhweiler, F. C., et al. 1999, ApJ, 523, L151 [NASA ADS] [CrossRef] [Google Scholar]
  21. Grady, C. A., Polomski, E. F., Henning, T., et al. 2001, AJ, 122, 3396 [NASA ADS] [CrossRef] [Google Scholar]
  22. Gressel, O., Nelson, R. P., Turner, N. J., & Ziegler, U. 2013, ApJ, 779, 59 [NASA ADS] [CrossRef] [Google Scholar]
  23. Haffert, S. Y., Bohn, A. J., de Boer, J., et al. 2019, Nat. Astron., 3, 749 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  24. Hashimoto, J., Tamura, M., Muto, T., et al. 2011, ApJ, 729, L17 [NASA ADS] [CrossRef] [Google Scholar]
  25. Hashimoto, J., Aoyama, Y., Konishi, M., et al. 2020, AJ, 159, 222 [NASA ADS] [CrossRef] [Google Scholar]
  26. Isella, A., Natta, A., Wilner, D., Carpenter, J. M., & Testi, L. 2010, ApJ, 725, 1735 [NASA ADS] [CrossRef] [Google Scholar]
  27. Juhasz, A., benisty, M., Pohl, A., et al. 2015, MNRAS, 451, 1147 [NASA ADS] [CrossRef] [Google Scholar]
  28. Kadam, K., Vorobyov, E., Regály, Z., Kóspál, Á., & Ábrahám, P. 2019, ApJ, 882, 96 [NASA ADS] [CrossRef] [Google Scholar]
  29. Kley, W., & Nelson, R. P. 2012, ARA&A, 50, 211 [NASA ADS] [CrossRef] [Google Scholar]
  30. Li, D., Pantin, E., Telesco, C. M., et al. 2016, ApJ, 832, 18 [NASA ADS] [CrossRef] [Google Scholar]
  31. Lomax, J. R., Wisniewski, J. P., Grady, C. A., et al. 2016, ApJ, 828, 2 [NASA ADS] [CrossRef] [Google Scholar]
  32. Maire, A.-L., Langlois, M., Dohlen, K., et al. 2016, Proc. SPIE, 9908, 990834 [CrossRef] [Google Scholar]
  33. Milli, J., Mouillet, D., Lagrange, A. M., et al. 2012, A&A, 545, A111 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  34. Muro-Arena, G. A., Ginski, C., Dominik, C., et al. 2020, A&A, 636, L4 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  35. Muto, T., Grady, C. A., Hashimoto, J., et al. 2012, ApJ, 748, L22 [NASA ADS] [CrossRef] [Google Scholar]
  36. Oppenheimer, B. R., Brenner, D., Hinkley, S., et al. 2008, ApJ, 679, 1574 [Google Scholar]
  37. Pavlov, A., Möller-Nilsson, O., Feldt, M., et al. 2008, Proc. SPIE, 7019, 39 [Google Scholar]
  38. Perrin, M. D., Schneider, G., Duchêne, G., et al. 2009, ApJ, 707, L132 [NASA ADS] [CrossRef] [Google Scholar]
  39. Phuong, N. T., Dutrey, A., di Folco, E., et al. 2020, A&A, 635, L9 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  40. Piétu, V., Guilloteau, S., & Dutrey, A. 2005, A&A, 443, 945 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  41. Rafikov, R. R. 2002, ApJ, 569, 997 [NASA ADS] [CrossRef] [Google Scholar]
  42. Rice, W. K. M., Lodato, G., & Armitage, P. J. 2005, MNRAS, 364, L56 [NASA ADS] [CrossRef] [Google Scholar]
  43. Rivière-Marichalar, P., Fuente, A., Baruteau, C., et al. 2019, ApJ, 879, L14 [NASA ADS] [CrossRef] [Google Scholar]
  44. Salyk, C., Herczeg, G. J., Brown, J. M., et al. 2013, ApJ, 769, 21 [NASA ADS] [CrossRef] [Google Scholar]
  45. Schmid, H. M., Joos, F., & Tschan, D. 2006, A&A, 452, 657 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  46. Soummer, R., Pueyo, L., & Larkin, J. 2012, ApJ, 755, L28 [NASA ADS] [CrossRef] [Google Scholar]
  47. Tang, Y. W., Guilloteau, S., Piétu, V., et al. 2012, A&A, 547, A84 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  48. Tang, Y.-W., Guilloteau, S., Dutrey, A., et al. 2017, ApJ, 840, 32 [NASA ADS] [CrossRef] [Google Scholar]
  49. van Holstein, R. G., Girard, J. H., de Boer, J., et al. 2020, A&A, 633, A64 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  50. Wisdom, J. 1980, AJ, 85, 1122 [NASA ADS] [CrossRef] [Google Scholar]
  51. Zhu, Z. 2015, ApJ, 799, 16 [NASA ADS] [CrossRef] [Google Scholar]

Appendix A: Complementary data

thumbnail Fig. A.1.

Complementary images of the AB Aur system. (a) ADI image (H band) in a large field of view (10″) multiplied with the square of the stellocentric distance. (b) High-pass filtering (gradient) of the DPI image to enhance the spirals S1 and S2. (c) Second epoch ADI observations in the K1K2 filters.

Table A.1.

Log of SPHERE observations.

Appendix B: Expression of a spiral arm

According to the theory of density waves generated by a gravitational perturber (Rafikov 2002) in protoplanetary disks, Muto et al. (2012) proposed an approximate linear relation to describe the shape of a spiral θ(r), as a function of the planet location (rc, θ0) and disk properties, as follows:

θ ( r ) = θ 0 sgn ( r r c ) h c × [ ( r r c ) 1 + β { 1 1 + β 1 1 α + β ( r r c ) α } ( 1 1 + β 1 1 α + β ) ] , $$ \begin{aligned} \theta (r)&= \theta _0-\frac{\text{ sgn}(r-r_c)}{h_c} \nonumber \\&\quad \times \left[ \left(\frac{r}{r_c}\right)^{1+\beta } \left\{ \frac{1}{1+\beta }-\frac{1}{1-\alpha +\beta } \left(\frac{r}{r_c}\right)^{-\alpha } \right\} \right. \nonumber \\&\quad \left. -\left( \frac{1}{1+\beta } - \frac{1}{1-\alpha +\beta } \right)\right] , \end{aligned} $$(B.1)

where hc is the disk aspect ratio at the planet location, α and β are the power-law exponents of the angular frequency and temperature profile dependence with r. In this paper, we assumed α = 1.5 for the Keplerian velocity and β = 0.25 to account for the temperature profile which varies as r−2β.

This formulation is now commonly used to interpret high-contrast images in which spiral patterns are observed, but it involves a number of assumptions, in particular, the idea that small grains observed in scattered light are coupled to the gas and that a single planet induces a single spiral arm. Clearly, this is not always supported by hydro-dynamical simulations, which result in much more complex morphologies. In a more realistic situation, a planet could produce several spiral arms (Dong et al. 2015) or planet formation could inevitably lead to multiple cores and, hence, to manifold spirals (Kadam et al. 2019).

All Tables

Table A.1.

Log of SPHERE observations.

All Figures

thumbnail Fig. 1.

Images of the AB Aur system obtained with SPHERE in polarized light (a, b) and unpolarized light (c). A large field of view (10″) is shown in panel a, where the polarized intensity has been multiplied with the square of the stellocentric distance (Qϕ × r2) to visually enhance the outer part of the disk. A narrower field of view (2″) of the Qϕ map is displayed in panel b, in which the inner spirals are labelled S1 and S2. Comparable ADI-processed image is shown in panel c for total intensity in the H band filter. Features f1 and f2 are discussed in Sect. 4. North is up, east is left.

In the text
thumbnail Fig. 2.

Full-field polarimetric image (same as in Fig. 1, Qϕ × r2) displayed in panel a, on top of which several structures previously detected are overlaid: Hashimoto et al. (2011) spirals (orange lines), out of plane accreting spirals from Tang et al. (2012; blue lines), and the dust ring (yellow ellipses) from the 0.9 mm ALMA image (Tang et al. 2017). The inner spiral arms (S1_CO, S2_CO) detected by Tang et al. (2017) are superimposed to the SPHERE polarimetric image Qϕ, further processed with a gradient algorithm to enhance the inner spirals (b), and the ALMA CO map (c), in a 2″ field of view. Features f1 and f2 are discussed in Sect. 4. North is up, east is left.

In the text
thumbnail Fig. 3.

Positions of features f1 and f2 overlaid on the polarimetric image shown for two intensity thresholds. The dusty ring is materialized with the dashed ellipses. For visualization purposes, the polarized intensity has been high-pass filtered and multiplied with the stellocentric distance (Qϕ × r).

In the text
thumbnail Fig. 4.

Same data as Fig. 3 with three different models of spirals matching S1 (green and blue lines) and S2 (purple line). The positions of the perturbers are indicated with a circle (panel a). Panel b: zoom in version of panel a, where only the best model reproducing the spiral twist in S1 is shown. Panel c: the contours of the CO map at the location of the candidate protoplanet are overlaid, together with the spiral model offset by 14.1° (dashed green line), corresponding to the 4 yr timelapse between ALMA and SPHERE observations.

In the text
thumbnail Fig. A.1.

Complementary images of the AB Aur system. (a) ADI image (H band) in a large field of view (10″) multiplied with the square of the stellocentric distance. (b) High-pass filtering (gradient) of the DPI image to enhance the spirals S1 and S2. (c) Second epoch ADI observations in the K1K2 filters.

In the text

Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.

Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.

Initial download of the metrics may take a while.