© ESO 2019

1 Introduction

Planet formation occurs in disks around young stellar objects. Interactions between planets and disks are very complex. Young planets are expected to cause rings, cavities, spirals, and disturbances in the velocity field and other features in the disk, which in turn may be used to infer the presence of these young planets. In the past few years, much evidence about this phase of planet formation has been accumulated because high-resolution images in the millimeter and sub-millimeter wavelength ranges have been provided by the Very Large Array (VLA) and the Atacama Large Millimeter Array (ALMA; see e.g. the case of HL Tau; ALMA Partnership 2015), and by high-contrast imagers such as the Gemini Planet Imager (GPI, Macintosh et al. 2014) and SPHERE (Spectro- Polarimetic High contrast imager for Exoplanets REsearch, Beuzit et al. 2008; see, e.g., Avenhaus et al. 2018). The literature on indirect evidence of the presence of planets is now becoming very rich, and nearby young stars surrounded by gas-rich disks are intensively studied for this purpose. In most cases, available data cannot fully eliminate alternative hypotheses, or the data have ambiguous interpretations (see, e.g., Bae et al. 2018 and Dong et al. 2018), although strong indirect evidence of the presence of planets from local disturbances of the velocity field have recently been considered for the case of HD 163296 (Pinte et al. 2018; Teague et al. 2018). In general, small grains are thought to be more strongly coupled with gas and are thus less sensitive to radial drift and concentration that can strongly affect large grains (see the discussion in Dipierro et al. 2018). For this reason, observations at short wavelengths provide an important complementary view of what can be seen with ALMA. On the other hand, a direct detection of still-forming planets embedded within primordial gas-rich disks, which is expected to be possible with high-contrast imaging in the near infrared (NIR), is still scarce; remarkable cases are LkCa-15 (Kraus & Ireland 2012; Sallum et al. 2015) and PDS-70 (Keppler et al. 2018; Müller et al. 2018; Wagner et al. 2018). In particular, in this second case, a clear detection of an accreting planet in the cavity between the inner and outer ring was obtained, making it an archetype for planet formation and planet-disk interactions. However, many cases remain ambiguous; a classical example is HD 100546 (see, e.g., Quanz et al. 2013a, 2015; Currie et al. 2014, 2015; Rameau et al. 2017; Sissa et al. 2018).

HD 169142 is a very young Herbig Ae-Be star with a mass of 1.65–2 M_⊙ and an age of 5–11 Myr (Blondel & Djie 2006; Manoj et al. 2007) that is surrounded by a gas-rich disk (i = 13°; Raman et al. 2006; PA = 5°; Fedele et al. 2017) that is seen almost face-on. The parallax is 8.77 ± 0.06 mas (Gaia DR2 2018). Disk structures dominate the inner regions around HD 169142 (see, e.g., Ligi et al. 2018). Figure 1 shows the view obtained from polarimetric observations: the left panel shows the Q_Φ image in the J band obtained by Pohl et al. (2017) using SPHERE on a linear scale, and the two rings are clearly visible. The right panel shows a pseudo-ADI image of the inner regions obtained by differentiating the Q_Φ image (see Ligi et al. 2018, for more details). Biller et al. (2014) and Reggiani et al. (2014) discussed the possible presence of a point source candidate at small separation (<0.2 arcsec from the star). However, the analysis by Ligi et al. (2018) based on SPHERE data does not support or refute these claims; in particular, they suggested that the candidate identified by Biller et al. (2014) might be a disk feature rather than a planet. Polarimetricimages with the adaptive optics system NACO at the Very Large Telescope (VLT; Quanz et al. 2013b), SPHERE (Pohl et al. 2017; Bertrang et al. 2018) and GPI (Monnier et al. 2017) show a gap at around 36 au, with an outer ring at a separation >40 au from the star. This agrees very well with the position of the rings obtained from ALMA data (Fedele et al. 2017); similar results were obtained from VLA data (Osorio et al. 2014; Macías et al. 2017). We summarize this information about the disk structure in Table 1 and call the ring at 0.17–0.28 arcsec from the star Ring 1 and the ring at 0.48–0.64 arcsec Ring 2. We remark that in addition to these two rings, both the spectral energy distribution (Wagner et al. 2015) and interferometric observations (Lazareff et al. 2017; Chen et al. 2018) show an inner disk at a separation smaller than 3 au. This inner disk isunresolved from the star in high-contrast images and consistent with ongoing accretion from it onto the young central star.While the cavities between the rings seem devoid of small dust, some gas is present there (Osorio et al. 2014; Macías et al. 2017; Fedele et al. 2017). Fedele et al. (2017) and Bertrang et al. (2018) have suggested the possibility that the gap between Rings 1 and 2 is caused by a planet with a mass slightly higher than that of Jupiter. However, this planet has not yet been observed, possibly because it is at the limit of or beyond current capabilities of high-contrast imagers. On the other hand, Bertrang et al. (2018) found a radial gap in Ring 1 at PA ~ 50° that might correspond to a similar radial gap found by Quanz et al. (2013b) at PA ~ 80°. The authors noted that if this correspondence were real, then this gap might be caused by a planet at about 0.14 arcsec from the star. So far, this planet has not been unambiguously detected either.

In this paper, we pursue a new view on the subject through analyzing high-contrast images. In particular, we underline that while polarimetric observations in the NIR and millimeter observations are best to reveal the overall structure of the disk, pupil-stabilized NIR observations where angular differential imaging can be applied may reveal fainter structures on a smaller scale. The risk of false alarms inherent to the image-processing procedures used in high-contrast imaging can be mitigated by comparing different sets of observations taken at intervals of months or years. In the case of HD 169142, this is exemplified by the study of Ligi et al. (2018), who identified a number of blobs within Ring 1. We have now accumulated a quite consistent series of observationsof this star with SPHERE that extends the set of data considered by Ligi et al. (2018). The observations have a comparable limiting contrast so that we may try to combine this whole data set to improve our knowledge of this system. Thecombination of different data sets acquired over a few years offers several advantages. In addition to verification of previous claims, we might try to detect persistent features around HD 169142 using a coincidence method to obtain a combined image that is deeper than the individual images and allows a quantitative discussion of the false-alarm probability of detected features. The expected orbital motion needs to be taken into account in this.

In Sect. 2 we describe observation and analysis methods. In Sect. 3 we present the main results about the blobs we detect around the star. In Sect. 4 we discuss the spiral arms within the disk and the possible connection to the blobs. Conclusions are given in Sect. 5.

2 Observation and data analysis

Data were acquired with the SPHERE high-contrast imager (Beuzit et al. 2008) at the ESO VLT Unit Telescope 3 within the guaranteed time observations used for the SHINE (SpHere INfrared survey for Exoplanets) survey (Chauvin et al. 2017). Data acquired up to 2017 have been described in Ligi et al. (2018). Here we add new data acquired in 2018 and study the system anew using different ways to combine different images. In these observations, we used SPHERE with both the Integral Field Spectrograph (IFS, Claudi et al. 2008) and the Infra-Red Dual Imaging and Spectrograph (IRDIS, Dohlen et al. 2008; Vigan et al. 2010) simultaneously. IFS was used in two modes: Y–J, that is, with spectra from 0.95 to 1.35 μm and a resolution of R ~ 50; and Y–H, with spectra from 0.95 to 1.65 μm and a resolution of R ~ 30. When IFS was in Y–J mode, IRDIS observed in the H2–H3 narrow bands (1.59 and 1.66 μm, respectively); when IFS was in Y–H mode, IRDIS observed in K1–K2 bands (2.09 and 2.25 μm, respectively). Hereafter we mainly consider data acquired with the IFS; IRDIS data are considered for the photometry in the K1-K2 bands. We considered the six best observations obtained for HD 169142 (see Table 2). Most of the epochs were obtained with an apodized Lyot coronagraph (Boccaletti et al. 2008; the field mask in YJH has a radius of 92 mas, and of 120 mas for the K-band coronagraph). Two of the observations (obtained in better observing conditions) were acquired without the coronagraph in order to study thevery central region around the star. The use of the coronagraph allows a better contrast at separation larger than ~ 0.1 arcsec. For all data sets, the observations were acquired in pupil-stabilized mode. In addition to the science data, we acquired three kind of on-sky calibrations: (i) a flux calibration obtained by offsetting the star position by about 0.5 arcsec, that is, out of the coronagraphic mask (point spread function, PSF, calibration). (ii) An image acquired by imprinting a bidimensional sinusoidal pattern (waffle calibration) on the deformable mirror. The symmetric replicas of the stellar images obtained by this second calibration allow an accurate determination of the star centers even when the coronagraphic field mask is in place. These calibrations were obtained both before and after the science observation, and the results were averaged. (iii) Finally, an empty field was observed at the end of the whole sequence to allow proper sky subtraction. This is relevant in particular for the K2 data sets.

Data were reduced to a 4D datacube (x, y, time, and λ) at the SPHERE Data Center in Grenoble (Delorme et al. 2017) using the standard procedures in the SPHERE pipeline (DRH, Pavlov et al. 2008) and special routines that recenter individual images using the satellite spot calibration, and correct for anamorphism, true north, and filter transmission. Faint structures can be detected in these images using differential imaging. Various differential imaging procedures were run on these data sets. We used here results obtained with a principal component analysis (PCA; see Soummer et al. 2012) applied to the whole 4D datacubes, which combines both angular and spectral differential imaging in a single step (ASDI-PCA, see Mesa et al. 2015). The PCA algorithm we used is the singular-value decomposition that generates the eigenvectors and eigenvalues that are used to reconstruct the original data. A principal components subset was used to generate an image with the quasi-static noise pattern that can then be subtracted from the original image. Clearly, the larger the number of principal components, the better the noise subtraction, but this also means that the signal from possible faint companion objects is cancelled out more strongly. Most of the results were obtained using 50 modes, but we also considered other numbers of modes (10, 25, 100, and 150 modes). To avoid spectrum distortion characteristics of the ASDI-PCA, photometry was obtained using a monochromatic PCA with only two modes for each spectral channel. Photometry was obtained with respect to the maximum of the PSF calibration corrected for the attenuation inherent to the PCA. The final step of the procedure was to obtain signal-to-noise ratio (S/N) maps from the IFS images obtained by making a median over wavelengths.

Finally, we also used the Q_Φ image obtained by Pohl et al. (2017) for astrometry, reduced as described in that paper and in Ligi et al. (2018). We note that this data set was obtained with the YJ field mask, whose radius is 72.5 mas.

3 Results

The following discussion is based on the application of differential imaging algorithms that allow detecting faint structure that is not otherwise easily detectable in the images. The typical contrast of Ring 1 that we were able to measure using simple subtraction of a reference image is about 1.5 × 10⁻³. The structures we consider in this paper are more than an order of magnitude fainter. They represent small fluctuations of the signal that cannot be detected without differential imaging.

Figure 2 shows the S/N maps obtained by applying the PCA ASDI algorithm to the IFS data for the individual epochs. The images have a linear scale from S∕N = 0 (dark) to S∕N = 5 (bright). These figures clearly show a similar pattern of bright spots, as well as a rotation of these spots with time. This suggests that a combination of the images that takes into account a Keplerian motion around the star should improve detection of the real pattern present in the data. The full solution is quite complex, leaving many free parameters, and may be attempted using an approach such as that considered by K-stacker (Nowak et al. 2018). However, a simplified approach that greatly reduces the number of free parameters is to assume that the system is seen face-on and that the orbits are circular: if the distance is known, the only free parameter is the stellar mass. This appears to be a reasonable approximation for disk-related features around HD 169142 because in this case, we only consider a fraction of the orbit. On the other hand, observations spread over a few years enableseparating static features that are due to radiative transfer effects from scattered-light fluctuations that are due to moving clouds or sub-stellar objects.

The upper panels of Fig. 3 show images of HD 169142 obtained by combining the six individual images, assuming the distance given by Gaia DR2, a mass of 1.7 M_⊙ (see below), and circular orbits.

Fig. 1 View of the surroundings of HD 169142 obtained from polarimetric observations. Left panel: QΦ image in the J band acquired with SPHERE (Pohl et al. 2017) on a linear scale. The two rings are clearly visible. Right panel: pseudo-ADI image of the inner regions obtained by differentiating the QΦ image (see Ligi et al. 2018, for more details). The white cross marks the position of the star, and the cyan circle shows the size of the coronagraphic mask. The other labels refer to the blobs we discuss in this paper that are also visible in these images. Thecolor scale of the differential image is five times less extended to show the faint structures better. In both panels, N is up and E to the left; a segment represents 1 and 0.5 arcsec in the left and right panel, respectively.

Open with DEXTER

Fig. 2 Signal-to-noise ratio (S/N) maps for the individual epochs for IFS. The individual S/N maps are obtained from the ASDI PCA algorithm using 50 components and making a median over the wavelength (see Mesa et al. 2015). Upper row, left panel: JD = 57180.17, middle panel: JD = 57201.12, and right panel: JD = 57499.34. Lower row, left panel: JD = 57566.15, middle panel: JD = 57873.30, and right panel: JD = 58288.19. In all panels, the central 0.1 arcsec is masked, the solid white line at the bottom represents 1 arcsec, the white cross represents the position of the star. N is up and E to the left.

Open with DEXTER

3.1 Coincidence images

To improve our ability of discerning faint signals, we combined data from different epochs using a coincidence map (see the lower panel of Fig. 3). The principle of this coincidence map is to start with S/N maps for individual epochs. We used S/N maps after correcting for the small-number-statistics effect using the formula by Mawet et al. (2014). The maps were then multiplied by each other. The S/N maps average to zero, with both positive and negative values for individual pixels. Of course, this may result in a false-positive signal for an even number of negative signals in the individual S/N maps. To avoid this problem, we arbitrarily set to negative the result for a given pixel when the signal for that pixel was negative for at least one epoch. Of course, this is not a realistic flux map: the aim is merely to identify consistent signals throughout all individual images.

To consider the orbital motion around the star over the three years covered by our observations, we divided the field into 65 rings, each one 2 IFS pixels wide (15 mas, i.e., about 1.8 au at the distance of HD 169142). For each ring, we rotated the S/N maps obtained at different epochs with respect to the first reference image according to Kepler’s third law. When the distance to the star was fixed, the only free parameter remaining in this model is the (dynamical) stellar mass. If there is a companion orbiting the star, the signal is maximized for a value of the mass that, if the assumptions made (circular motion seen face on) are correct, is the dynamical mass of the star. If these hypotheses are not correct, the estimate of the mass is incorrect, by a value that depends on the real orbital parameters. We adopted a mass of 1.7 M_⊙, the Gaia parallax, and clockwise rotation, as indicated by the analysis of motion of disk features in Ligi et al. (2018); see also Macías et al. (2017). Figure 3 shows the coincidence map and a mean of the S/N maps for the six epochs for this value of the mass.

Fig. 3 Upperleft panel: median over time of the wavelength-collapsed images of HD 169142 obtained with an ADI PCA algorithm, one mode per wavelength. Upper right panel: image obtained by averaging the S/N maps for the individual epochs for IFS. The individual S/N maps are obtained from the ASDI PCA algorithm using 50 components and making a median of the wavelength (see Mesa et al. 2015), and they were rotated for Keplerian motion to the last image before making the median. Lower panel: coincidence image obtained from the same data set. In all panels, the solid line at the bottom of each panel represents 1 arcsec, and a white cross shows the position of the star. N is up and E to the left.

Open with DEXTER

3.2 Blob detection

A quite large number of blobs can be found around HD 169142. Several of them are found consistently in all individual images and are also visible in the J-band Q_Φ image seen in Fig. 1; some of them have been identified and discussed by Ligi et al. (2018). We fixed our attention on four of them (see Figs. 1 and 4 for their definition). The two brightest blobs (B and C) are within Ring 1 (see Fig. 4) and have been identified by Ligi et al. (2018); they called them blobs A and B, respectively. Our blob A is closer to the star than Ring 1. Blob D is between Rings 1 and 2. All of these blobs appear to be slightly extended. We verified in the individual images that this is not an artifact caused by combining individual images. For instance, when we consider the best set of data (the last set from June 2018), the FWHM of blobs B and C can be measured with reasonable accuracy at about 40 mas, which is significantly larger than expected for a point source at this separation (about 26 mas, after applying differential imaging). To better estimate the physical size of the blobs, we compared the FWHM measured inour differential images with the FWHMs obtained for fake blobs that are the result of convolving Gaussian profiles with the observed PSF inserted into the images at the same separation but at a different position angle, and processed through the same differential imaging procedure. We repeated this procedure for the images obtained considering 50 modes (best image for detection) and with a less aggressive image where only 25 modes were considered, which better conserves the original shape of the blobs. In this way, we found that the FWHM of blobs B and C is the same as that of fake Gaussian blobs with an intrinsic FWHM of 42 and 30 mas for blobs B and C, respectively. However, these are average values for tangential and radial profiles (with respect to the star): both blobs appear elongated in the tangential (rotation) direction with axis ratios of 1.4 (blob B) and 1.9 (blob C)¹ . The uncertainty on this size estimate is of about 7 mas, as obtained by comparing the results obtained in individual images. This size corresponds to ~ 4–5 au, with an uncertainty of about 1 au. This result should be considered with some caution because the light distribution of the blobs might be not well reproduced by Gaussians.

In Fig. 5 we show a zoom of the region around blob D in a two-color image. The blue structures visible in the image are obtained through IFSin the Y , J, and H bands. They could be interpreted as stellar light scattered by a (dusty) spiral structure around a protoplanet that is accreting material funnelled through the spiral arm from the disk. This interpretation agrees with the detection in the Q_Φ image (see Fig. 1). In the same image, the red structures are obtained through IRDIS observing in the K2 band. In particular, the structure in the green circle that appears to be much more similar to a point source might indicate a planetary photosphere. The position of the blob in the K2 image is Sep = 332 mas, PA = 34.9°, which is notthe photocenter at shorter wavelengths. Even if this interpretation is speculative (there are other structures in this image that we consider as noise), various circumstantial arguments discussed below possibly support it. We return to this point in the next section.

There is of course some probability that these detections are spurious. In order to estimate the false-alarm probability (FAP), we proceeded as follows. First, we fixed the stellar mass at the value given by fitting isochrones (1.7 M_⊙). With this assumption, the prediction for the orbital motion is fully independent of the SHINE data set. We derotated the individual images to the same epoch using the same approach as described above (ring by ring). We searched for signals in the final coincidence data set using the FIND procedure in IDL. We recovered the detection of the candidate. We ordered the different epochs according to the value of the S/N at the candidate positions (separately for each candidate). We then used binomial statistics on the remaining epochs (i.e., excluding the reference with the highest S/N), considering as number of trials the number of pixels with an S/N higher than the S/N measured in the candidate position in the image giving the highest S/Nvalue at this position. To estimate the probability in the binomial statistics, we considered the product ∏ _c of the S/N rankings in the pixel corresponding to the candidate position in the remaining images, and compared this product to a similar product ∏ _r obtained from random extractions. We repeatedthe random extraction 10⁷ times, and assumed that the probability of success is given by the fraction of cases where ∏ _r < ∏ _c.

With this approach, we obtained the FAP values listed in the second column of Table 4; values for blobs B, C, and D are highly significant.

	Fig. 4 Same as the lower panel of Fig 3, but showing the edges of the two disk rings (Ring 1 in green, Ring 2 in magenta). The ring edges are drawn according to Fedele et al. (2017). The blobs are labeled. The white solid line represents 1 arcsec, and the white cross shows the position of the star. N is up and E to the left.
Open with DEXTER

Fig. 5 Zoom of a two-color image of the region around blob D. This image was constructed using the K2 observation of JD = 2458288.19 (red) and the weighted sum of all the IFS images (collapsed against wavelength) and rotated for a Keplerian motion assuming that the star has a mass of 1.7 M⊙ (blue). This last image is for the same epoch as the K2 observation. For clarity, the region within 0.28 arcsec from the star (i.e., within the outer edge of Ring 1) was masked in the K2-band image. The green circle is centered on the position of the blob measured in the K2 image. We note the different aspect and small offset between the position of the blob in the K2 image with respect to that at shorter wavelength. The white solid line represents 1 arcsec. N is up and E to the left. The white cross marks the position of the star.

Open with DEXTER

3.3 Blob astrometry

All these blobs rotate around the star. This can be shown by measuring their position in the individual images (see Table 3). We used the IDL FIND algorithm that uses marginal Gaussian distributions to measure the position of the spot centers in the ASDI 50 components images obtained at the various epochs. In addition, we also measured the blob positions in the polarimetric image. We found that the rotational speed decreases with separation (see Table 4 and Fig. 6), as expected for Keplerian motions. Figure 7 shows the run of the angular speed as a function of separation for the four blobs around HD 169142; these values are on the disk plane. Overimposed are predictions for circular Keplerian motion for three different values for the mass of the star (1.5, 2.0, and 2.5 M_⊙ ).

When we interpret the observed angular motion as Keplerian circular orbits in the disk plane, we can determine the mass of HD 169142. When we use the three blobs B, C, and D, the mass of the star is 1.85 ± 0.25M_⊙ (we did not use blob A here because it has too few astrometric points). When we add the uncertainties that are due to parallax (0.05 M_⊙ ) and disk inclination (0.09 M_⊙), the result is 1.85 ± 0.27M_⊙.

The mass estimated by this procedure might be underestimated because the photocenter of the blobs might be closer to the star than their center of mass and the mass determination depends on the cube of the separation. In particular, as noted above, we may interpret blob D as the accretion flows on a planet along a spiral arm (see the next section); in this case, the photocenter is dominated by the leading arm, which is at about 313 mas from the star, while the trailing arm is at 347 mas when it is deprojected on the disk plane. The putative planet would be in the middle of the two arms, that is, at 330 mas from the star, yielding a mass estimate that is ~ 10% higher than listed in Table 4. We note that the separation measured in the K2 band agrees very well with this interpretation.

The mass determined from blob motion is slightly higher than but in agreement within the error bars with the mass that fits photometry. To show this, we determined the stellar mass by minimizing the χ² with respect to the main-sequence values considered by Pecaut & Mamajek (2013). We considered the Gaia DR2 parallax and included an absorption term A_V multiplied for the reddening relation by Cardelli et al. (1989). We also left free the ratio between the stellar and the main-sequence radius. The best match is with an F0V star (T_eff = 7220 K), with A_V = 0.25 mag and a radiusthat is 0.97 times the radius of the main-sequence star. According to Pecaut & Mamajek (2013), the mass of an F0V star is 1.59 M_⊙. This spectral type compares quite well with the most recent determinations (A7V: Dent et al. 2013; A9V: Vieira et al. 2003; F0V: Paunzen et al. 2001; F1V: Murphy et al. 2015) and with the temperature determined by Gaia (T_eff = 7320 ± 150 K), but it is much later than the B9V spectra type proposed by Wright (2003).

For comparison, other determinations of the mass of HD 169142 are 2.0 M_⊙ (Manoj & Bhatt 2005), 2.28 M_⊙ (Maaskant et al. 2013), 1.8 M_⊙ (Salyk et al. 2013), and 2.0 M_⊙ (Vioque et al. 2018). The mass adopted by Ligi et al. (2018) is 1.7 M_⊙. We note that these values were obtained assuming distances different from the distance given by Gaia DR2: for instance, Maaskant et al. (2013) adopted a distance of 145 pc, which is 27% larger than the Gaia DR2 value considered here. On the other hand, the value used by Ligi et al. (2018) was taken from Gaia DR1 and it is only 3% longer than that from Gaia DR2. Hereafter, we adopt a mass of 1.7 M_⊙ for HD 169142, which is the same value as was considered by Ligi et al. (2018).

We also note that the projected rotational velocity of the star V   sini = 50.3 ± 0.8 km s⁻¹ determined from the HARPS spectra (see Appendix A) is high when we consider that the star is likely seen close to the pole. This value agrees quite well with literature values (V   sini = 55 ± 2 km s⁻¹; Dunkin et al. 1997a,b). When we assume that the stellar rotation is aligned with the disk, the equatorial rotational velocity is 224 km s⁻¹, which is at the upper edge of the distribution for F0 stars. For a discussion, see Grady et al. (2007).

	Fig. 6 Variation in PA of the blobs with time. Upper left panel: blob A. Upper right panel: blob B. Lower left panel: blob C. Lower right panel: blob D. The dashed lines are best-fit lines through the points. Dash-dotted lines are predictions for circular orbits assuming a mass of 1.85 M⊙ for the star.
Open with DEXTER

	Fig. 7 Run of the angular speed as a function of separation for the four blobs around HD 169142; these values are on the disk plane. Overimposed we show predictions for circular Keplerian motion for three different values for the stellar mass (1.5, 2.0, and 2.5 M⊙ ).
Open with DEXTER

3.4 Blob photometry

We measured the magnitude of the sources in various bands by weighting the results obtained from the different epochs according to the quality of the images. The magnitudes refer to a 3 × 3 pixel area centered on each object and are obtained by comparison with those of simulated planets inserted into the image at 0.2 and 0.3 arcsec from the star and run through the same differential imaging algorithm. The underlying assumption is that the blobs are point sources, while they are likely slightly extended. These results should then be taken with caution. Using the fake blob procedure described in Sect. 3.2, we estimated that the brightness is underestimated by about a factor of ~ 2.8 because of this effect for blobs B and C, that is, these blobs are likely ~1.1 magnitude brighter than estimated when we assume that they are point sources. The effect is likely slightly smaller for blob D because it isfarther away from the star. We summarize the results in Table 5; error bars are the standard deviation of the mean of the results obtained at different epochs. All the blobs have a rather flat, only slightly reddish contrast with respectto the star (see Fig. 8). Results are consistent with stellar light scattered by grains with a size on the order of a micron or smaller if stellar light is extinguished between the star and blobs or the blobs and us. Under the hypothesis (not demonstrated) that they are optically thick, the albedo required to reproduce observations of blobs A, B, and C is about 0.1.

Blob D is about two magnitudes fainter than expected from this consideration, suggesting that either it receives less light from the star (e.g., because of absorption by Ring 1) or it is not optically thick. For this blob, we obtained contrastsof 13.31 ± 0.28 mag in Y , 13.34 ± 0.23 mag in J band, 13.29 ± 0.21 mag in H band, 12.94 ± 0.5 mag in K1 band, and 12.35 ± 0.5 mag in K2 band (the last two values being obtained from the IRDIS data set). As expected, the object is beyond the 5σ contrast limit in each individual image. However, we expect a detection with an S/N in the range from 2.3 to 3.7 in the individual images, and at an S/N~ 6 in the combination of the images. It is then not surprising that we detected it only by combining them. While error bars are quite large, the absolute K1 and K2 magnitude of 14.07 ± 0.50 and 13.48 ± 0.50 mag corresponds to a ~ 3M_J object usingdusty isochrones (Allard et al. 2001) with an age of 5 Myr, which is at the lower edge of the age range according to Blondel & Djie (2006) and Manoj et al. (2007). This model has an effective temperature of about 1260 K (see Fig. 9). Of course, the mass estimated from photometry depends on the model, the age used to derive it and the possible extinction, and it assumes that the object is in hydrostatic equilibrium, which may be incorrect for a very young planet. This result is then highly uncertain.

	Fig. 8 Contrast of blobs as a function of wavelength. Blob B: asterisks and dotted line. Blob C: triangles and solid line. Blob D: diamonds and dashed line.
Open with DEXTER

	Fig. 9 Absolute magnitude of blob D in various bands (diamonds). The solid line is the prediction for a 3 MJ , 5 Myr old planet using dusty isochrones by Allard et al. (2001).
Open with DEXTER

3.5 Comparison with previous detection claims

We note that none of these blobs coincides with either the sub-stellar companions proposed by Biller et al. (2014) and Reggianiet al. (2014), nor with the structure observed by Osorio et al. (2014). More in detail, after taking into account their motion (see Sect. 3.3), the expected position angles for blobs A, B, and C at the observation epochs of Biller et al. (2014) and Reggiani et al. (2014; both acquired at an epoch about 2013.5), are 276°, 34°, and 324°, respectively (blob D is much farther away from the star). For comparison, the object of Biller et al. (2014) is at PA = 0 ± 14° (separation of 110 ± 30 mas) and the object of Reggiani et al. (2014) is at PA = 7.4 ± 11.3 degree (separation of 156 ± 32 mas). In addition, the objects proposed by Biller et al. (2014) and Reggiani et al. (2014), with a contrast of ΔL ~ 6.5, are brighter than our blobs B and C, even after the finite-size correction is taken into account, see Sect. 3.4. However, the object proposed by Reggianiet al. (2014) might be the combination of blobs B and C, within the errors of their astrometry; the combination of their luminosity is also not that far from the value of Reggiani et al. (2014). We note that the resolution of their observation is lower than ours because they observed at much longer wavelength, and their object appears elongated (in the E-W direction, i.e., the direction expected at the epoch of their observation) in their published image beyond the diffraction limit.

On the other hand, the inner and brighter object detected by Biller et al. (2014) is too close to the star to coincide with any of the objects we observed, while a fainter object they found might be blob B, as discussed by Ligi et al. (2018). However, when we examine the image published by Biller et al. (2014), it seems that the two brightest sources have a relative separation and orientation that coincides with those of blobs B and C. In this case, the fainter object should be blob B (as discussed in Ligi et al. 2018) and the brighter object might coincide with our blob C. Of course, this would require that the stellar position in their images does not correspond with the position assumed in their paper.

4 Spiral arms

Most bright features in the coincidence images (including blobs B and D) can be reproduced by a three-arm spiral design (see below), while blob C differs slightly. A similar three-arm structure is predicted by models for not very massive companions (Fung & Dong 2015) and it has been observed around other stars (see, e.g., the case of MWC758 recently published by Reggiani et al. 2018). While it is not at all obvious that spiral arms indicate a planet (see, e.g., Dong et al. 2018), we may interpret it as due to a planet in the location of blob D. The structure of this blob appears to resemble the structure expected for an accreting object with a leading and a trailing arm. If this hypothesis were correct, the radial separation between leading (Sep = 310 mas) and trailing (Sep = 343 mas) arms should be about twice the Hill radius (see Machida et al. 2010); the Hill radius would then be 16.5 ± 5 mas and the planet mass in the range 0.25–1.6 M_J, which is lower than the mass estimated with DUSTY isochrones. Because the Hill radius is not accurately estimated and the dependence of the mass on the Hill radius is strong, the error on the planetary mass is quite large. The photosphereof such an object would be too faint for detection in YJH, while it might be compatible with detection in K1 and K2 bands.

4.1 Separation of spiral arms

We may also estimate the mass of the object exciting the spiral design observed in Ring 1 by different criteria, using the calibration by Fung & Dong (2015; their Eq. (9)). After transforming into a polar coordinates system (see Fig. 10), we could identify the three spiral arms, which we may call primary, secondary, and tertiary, following the approach of Fung & Dong (2015). The view in Cartesian coordinates is given in Fig. 11. Position angle and separation of the arms in some reference positions are given in Table 6.

These arms may be density waves excited by a planet at the location of blob D, which is indeed along the primary arm of the spiral design: the predicted PA at the separation of blob D, 334 mas, is 36 ± 6°, in very good agreement with the observed value of ~ 35° (as measured in the K2 band).

The phase difference between the primary and secondary arm (127.2 ± 3.1°) can be used to estimate the mass of the planet exciting the spiral design, using the calibration by Fung & Dong (2015). We obtain a mass ratio of q = 0.0030 ± 0.0004, which translates into a mass of M_p = 5.1 ± 1.1M_J, adopting the stellar mass derived above. The phase difference between secondary and tertiary arms (69.9 ± 3.9°) agrees with the expectations by Fung & Dong (2015) given the pitch angle and the expected ratio for resonances 1:2 and 1:3.

	Fig. 10 Median over time of the individual S/N maps in polar coordinates. Each image has been rotated to the last image for the rotation angle of blob D before the median was made. Arrows mark the location of the primary (white), secondary (cyan), and tertiary arms (yellow). The location of the blobs is marked.
Open with DEXTER

	Fig. 11 Upper panel: median over time of the individual S/N maps; each image has been rotated to the last image for the rotation angle of blob D before the median was made. Lower panel: same as the upper panel, with a three-arm spiral design overplotted. The putative planet is at the position of blob D. In both panels, N is up and E to the left, and a white cross marks the stellar position.
Open with DEXTER

4.2 Pitch angle

The pitch angle is the angle between a spiral arm and the tangent to a circle at the same distance from the star. Zhu et al. (2015) showed that the pitch angle can be used to estimate the mass of the planet exciting the spiral design. We measured the pitch angle at a separation of 183 mas to be 17.5 ± 1.7°. This separation is about r∕r_p = 0.55. This value agrees with the results they obtained from their simulations for a mass ratio of q = 0.006, supporting the mass determination obtained from the separation of primary and secondary spiral arm; moreover, the larger pitch for the tertiary arm agrees with expectations from models.

4.3 Disk gap

Using the relation by Kanagawa et al. (2016), we expect that there is a planet at ~ 0.36 arcsec from the star with a mass ratio with respect to the star of , where h_p ∕R_p is the disk thickness and α is the disk viscosity. For R_p = 0.36 arcsec, W = 0.2 arcsec, h_p ∕R_p = 0.05, and α = 1E − 3, a value of q = 0.00044 is obtained, which means a planet of 0.75 M_J.

Dong & Fung (2017) considered the case of HD 169142 and concluded for a value of q² ∕α = 1.1E − 4 for R_p = 0.37 arcsec, W = 0.17 arcsec, and h_p ∕R_p = 0.079. For α = 1E − 3, their formula implies q = 0.00033, which suggests a 0.56 M_J planet. We note that the formula by Dong & Fung (2017) produces planets that are smaller by a factor of 2.6 with respect to that by Kanagawa et al. (2016); however, the value they suggest for h_p ∕R_p is higher than considered above. The value considered by Dong & Fung (2017) is similar to the value obtained by Fedele et al. (2017) by modeling the ALMA observations (h_p∕R_p = 0.07).

There are considerable uncertainties in these formulas that are due to the exact values to be adopted for R_p , W, h_p ∕R_p, α, and the difference of a factor of 2.5 in the constant factor. While a mass around 1 M_J seems favored, we cannot exclude values different by as much as an order of magnitude. We conclude that a planet with about one Jupiter mass likely causes the gap seen in HD 169142, but its mass is not yet well defined from the gap alone.

4.4 Summary of mass determination

A summary of the mass determinations is given in Table 7. All these estimates are quite uncertain. The higher values are given by thespiral arm parameters. If we make an harmonic mean of the various estimates, we would conclude for a planet with a mass of M_J. This mass seems lower than what we can detect with our SPHERE images (about 3 M_J from photometry), but is within the error bar. This value is also within the range 1–10 M_J suggested by Fedele et al. (2017) to justify the dust cavity observed with ALMA between rings 1 and 2, and it is on the same order as the missing mass in the disk within the gap, as given by their disk model (4.3 M_J). For comparison, we note that if we were to try to interpret the spiral arms of MWC 758 (Reggiani et al. 2018) using the same approach, we would conclude for a more massive faint companion because in that case the separation between the primary and secondary arm is much closer to 180°.

5 Conclusion

We performed an analysis of faint structures around HD 169142 that are persistent among several data sets obtained with SPHERE and analyzed them using differential image techniques. We found a number of blobs that rotate around the star as well as spiral arms. These structures represent small fluctuations of the overall disk structure around this star. The blobs are found to consistently rotate around the star with Keplerian circular motion.

Although we cannot exclude other hypotheses, blob D might correspond to a low-mass (~ 1–4M_J, best guess of 2.2 M_J), 5 Myr old, and still-accreting planet at about 335 mas (38 au) from the star, causing the gap between Rings 1 and 2 and exciting the spiral arm design observed within Ring 1. The separation between the outer edge of Ring 1 and blob D is 55 mas, which is about twice the proposed value for the Hill radius of the planet. There is a clear excess of flux at short wavelengths with respect to the flux expected for a planetary photosphere (see Fig. 9). In our proposed scenario, the planetary photosphere is not detected in YJH band, where we only see the accreting material fueling through the spiral arm and reflecting star light (consistent with its detection in the Q_Φ image), while it might have been detected in the K1 and K2 bands. A planet of 2.2 M_J at 335 mas (38 au) from HD 169142 would have a Hill radius of about 25 mas (3.2 au). A disk around such an object would have an FWHM slightly larger than the resolution limit of SPHERE and may well reflect some 10⁻⁵ of the stellar light, which is required to justify the flux observed in the YJH bands. On the other hand, it is also possible that no other planet exists, and we merely observe a dust cloud. Detection of a planet could be confirmed by observations in the L′ band. According to the AMES-dusty isochrones (Allard et al. 2001), a 5 Myr old planet of 2.2 M_J should have an absolute L′ magnitude of ~ 11 mag. The contrast in the L′ band should then be of 10.1 mag, which is 3.7 mag fainter than the objects proposed by Biller et al. (2014) and Reggiani et al. (2014) and likely too faint for a detection in their data set. However, a future deeper data set can solve this issue.

The location of blob B (and C to a lesser extent) suggests at first sight that the blobs might be related to the secondary and tertiary spiral arms (see, e.g., Crida et al. 2017). If this were the case, they would follow the same angular speed as the perturbing object, that is, the putative planet. However, we showed (along with Ligi et al. 2018) that those blobs follow a Keplerian motion appropriate for their separation from the star.

Finally, we note that Ligi et al. (2018) proposed that blobs B and C could be vortices (Meheut et al. 2012). This explanation might very well be true. Another scenario might be suggested by the possibility that they are in 1:2 resonance with a putative planet related to blob D. It concerns planetesimals or asteroid giant impacts that generate dust clouds. This might be a manifestation of the general phenomenon of planetesimal erosion that is expected to follow the formation of giant planets (see, e.g., Turrini et al. 2012, 2018). However, the probability of observing such clouds is low in a gas-rich disk such as that of HD 169142 because large planetesimals are required to generate clouds as large as blobs B and C, unless the impact occurs far from the disk plane. The debris cloud from an impact roughly expands until the debris sweeps a gas mass that is no more than an order of magnitude higher than the mass of the debris itself.Because the volume of the clouds is at least one hundredth of the total volume of ring 1, this requires that the mass of the interacting bodies is higher than 1/1000 of the mass of the disk when we assume a disk gas-to-dust ratio of unity and that the impact occurs close to the disk plane. The impacting bodies should then have a masson the order of that of Mars or at least the Moon. Since it is not likely that many such objects are present in the disk of HD 169142, the probability of observing one or even more similar debris clouds is likely very low.

Acknowledgements

The authors thank A. Pohl for allowing them to use the original reduction of the Q_Φ data set and the ESO Paranal Staff for support for conducting the observations. E.S., R.G., D.M., S.D. and R.U.C. acknowledge support from the “Progetti Premiali” funding scheme of the Italian Ministry of Education, University, and Research. E.R. and R.L. are supported by the European Union’s Horizon 2020 research and innovation programme under the Marie Skłodowska-Curie grant agreement No 664931. This work has been supported by the project PRININAF 2016 The Cradle of Life – GENESIS-SKA (General Conditions in Early Planetary Systems for the rise of life with SKA). The authors acknowledge financial support from the Programme National de Planétologie (PNP) and the Programme National de PhysiqueStellaire (PNPS) of CNRS-INSU. This work has also been supported by a grant from the French Labex OSUG020 (Investissements d’avenir - ANR10 LABX56). The project is supported by CNRS, by the Agence Nationale de la Recherche (ANR-14-CE33-0018). This work is partly based on data products produced at the SPHERE Data Centre hosted at OSUG/IPAG, Grenoble. We thank P. Delorme and E. Lagadec (SPHERE Data Centre) for their efficient help during the data reduction process. 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 Astronomico di Padova (Italy), Observatoire de Genéve (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).

Appendix A Radial and rotational velocities for HD 169142

We measured relative radial velocities and projected rotational velocities of HD 169142 from a series of eight high-resolution archive spectra of HD 169142 acquired with the HARPS spectrograph at the ESO 3.6 m telescope in La Silla in 2008 (Program 080.C-0712, PI: Desort). The spectra were reduced using the ESO pipeline. Radial velocities were obtained by cross-correlating them with the average of the last two spectra that have the highest S/N. The radial velocities are then relative. Relevant data are given in Table A.1.

References

All Tables

All Figures

	Fig. 1 View of the surroundings of HD 169142 obtained from polarimetric observations. Left panel: QΦ image in the J band acquired with SPHERE (Pohl et al. 2017) on a linear scale. The two rings are clearly visible. Right panel: pseudo-ADI image of the inner regions obtained by differentiating the QΦ image (see Ligi et al. 2018, for more details). The white cross marks the position of the star, and the cyan circle shows the size of the coronagraphic mask. The other labels refer to the blobs we discuss in this paper that are also visible in these images. Thecolor scale of the differential image is five times less extended to show the faint structures better. In both panels, N is up and E to the left; a segment represents 1 and 0.5 arcsec in the left and right panel, respectively.
Open with DEXTER
In the text

	Fig. 2 Signal-to-noise ratio (S/N) maps for the individual epochs for IFS. The individual S/N maps are obtained from the ASDI PCA algorithm using 50 components and making a median over the wavelength (see Mesa et al. 2015). Upper row, left panel: JD = 57180.17, middle panel: JD = 57201.12, and right panel: JD = 57499.34. Lower row, left panel: JD = 57566.15, middle panel: JD = 57873.30, and right panel: JD = 58288.19. In all panels, the central 0.1 arcsec is masked, the solid white line at the bottom represents 1 arcsec, the white cross represents the position of the star. N is up and E to the left.
Open with DEXTER
In the text

	Fig. 3 Upperleft panel: median over time of the wavelength-collapsed images of HD 169142 obtained with an ADI PCA algorithm, one mode per wavelength. Upper right panel: image obtained by averaging the S/N maps for the individual epochs for IFS. The individual S/N maps are obtained from the ASDI PCA algorithm using 50 components and making a median of the wavelength (see Mesa et al. 2015), and they were rotated for Keplerian motion to the last image before making the median. Lower panel: coincidence image obtained from the same data set. In all panels, the solid line at the bottom of each panel represents 1 arcsec, and a white cross shows the position of the star. N is up and E to the left.
Open with DEXTER
In the text

	Fig. 4 Same as the lower panel of Fig 3, but showing the edges of the two disk rings (Ring 1 in green, Ring 2 in magenta). The ring edges are drawn according to Fedele et al. (2017). The blobs are labeled. The white solid line represents 1 arcsec, and the white cross shows the position of the star. N is up and E to the left.
Open with DEXTER
In the text

	Fig. 5 Zoom of a two-color image of the region around blob D. This image was constructed using the K2 observation of JD = 2458288.19 (red) and the weighted sum of all the IFS images (collapsed against wavelength) and rotated for a Keplerian motion assuming that the star has a mass of 1.7 M⊙ (blue). This last image is for the same epoch as the K2 observation. For clarity, the region within 0.28 arcsec from the star (i.e., within the outer edge of Ring 1) was masked in the K2-band image. The green circle is centered on the position of the blob measured in the K2 image. We note the different aspect and small offset between the position of the blob in the K2 image with respect to that at shorter wavelength. The white solid line represents 1 arcsec. N is up and E to the left. The white cross marks the position of the star.
Open with DEXTER
In the text

	Fig. 6 Variation in PA of the blobs with time. Upper left panel: blob A. Upper right panel: blob B. Lower left panel: blob C. Lower right panel: blob D. The dashed lines are best-fit lines through the points. Dash-dotted lines are predictions for circular orbits assuming a mass of 1.85 M⊙ for the star.
Open with DEXTER
In the text

	Fig. 7 Run of the angular speed as a function of separation for the four blobs around HD 169142; these values are on the disk plane. Overimposed we show predictions for circular Keplerian motion for three different values for the stellar mass (1.5, 2.0, and 2.5 M⊙ ).
Open with DEXTER
In the text

	Fig. 8 Contrast of blobs as a function of wavelength. Blob B: asterisks and dotted line. Blob C: triangles and solid line. Blob D: diamonds and dashed line.
Open with DEXTER
In the text

	Fig. 9 Absolute magnitude of blob D in various bands (diamonds). The solid line is the prediction for a 3 MJ , 5 Myr old planet using dusty isochrones by Allard et al. (2001).
Open with DEXTER
In the text

	Fig. 10 Median over time of the individual S/N maps in polar coordinates. Each image has been rotated to the last image for the rotation angle of blob D before the median was made. Arrows mark the location of the primary (white), secondary (cyan), and tertiary arms (yellow). The location of the blobs is marked.
Open with DEXTER
In the text

	Fig. 11 Upper panel: median over time of the individual S/N maps; each image has been rotated to the last image for the rotation angle of blob D before the median was made. Lower panel: same as the upper panel, with a three-arm spiral design overplotted. The putative planet is at the position of blob D. In both panels, N is up and E to the left, and a white cross marks the stellar position.
Open with DEXTER
In the text