Multiple rings in the transition disk and companion candidates around RX J1615.3-3255.

Multiple rings in the transition disk and companion candidates around RX J1615.3-3255.

High contrast imaging with VLT/SPHERE.thanks: Based on observations made with ESO Telescopes at the La Silla Paranal Observatory under programme IDs 095.C-0298(A), 095.C-0298(B), and 095.C-0693(A) during guaranteed and open time observations of the SPHERE consortium, and on NACO observations: program IDs: 085.C-0012(A), 087.C-0111(A), and 089.C-0133(A).
J. de Boer 1 Leiden Observatory, Leiden University, P.O. Box 9513, 2300RA Leiden, The Netherlands.
1deboer@strw.leidenuniv.nl
   G. Salter 2Aix Marseille Université, CNRS, LAM (Laboratoire d’Astrophysique de Marseille) UMR 7326, 13388, Marseille, France 2    M. Benisty 3Université Grenoble Alpes, IPAG, F-38000 Grenoble, France 3    A. Vigan 2Aix Marseille Université, CNRS, LAM (Laboratoire d’Astrophysique de Marseille) UMR 7326, 13388, Marseille, France 2    A. Boccaletti 4LESIA, CNRS, Observatoire de Paris, Université Paris Diderot, UPMC, 5 place J. Janssen, 92190 Meudon, France. 4    P. Pinilla 1 Leiden Observatory, Leiden University, P.O. Box 9513, 2300RA Leiden, The Netherlands.
1deboer@strw.leidenuniv.nl
   C. Ginski 1 Leiden Observatory, Leiden University, P.O. Box 9513, 2300RA Leiden, The Netherlands.
1deboer@strw.leidenuniv.nl
   A. Juhasz 5Institute of Astronomy, Madingley Road, Cambridge CB3 OHA, UK. 5    A.-L. Maire 6Max-Planck-Institut fuer Astronomie, Koenigstuhl 17, 69117 Heidelberg, Germany 6    S. Messina 7INAF – Osservatorio Astronomico di Padova, Vicolo dell’Osservatorio 5, 35122 Padova, Italy. 7    S. Desidera 8INAF Catania Astrophysical Observatory, via S. Sofia 78, 95123 Catania, Italy. 8    A. Cheetham 9Observatoire de Genève, Université de Genève, 51 Chemin des Maillettes, 1290, Versoix, Switzerland. 9    J. H. Girard 10European Southern Observatory, Alonso de Cordova 3107, Casilla 19001 Vitacura, Santiago 19, Chili. 10    Z. Wahhaj 10European Southern Observatory, Alonso de Cordova 3107, Casilla 19001 Vitacura, Santiago 19, Chili. 10    M. Langlois 11Observatoire de Lyon, Centre de Recherche Astrophysique de Lyon, Ecole Normale Supérieure de Lyon, CNRS, Université Lyon 1, UMR 5574, 9 avenue Charles André, Saint-Genis Laval, 69230, France 11 2 2    M. Bonnefoy 3Université Grenoble Alpes, IPAG, F-38000 Grenoble, France 3    J.-L. Beuzit 3Université Grenoble Alpes, IPAG, F-38000 Grenoble, France 3    E. Buenzli 6Max-Planck-Institut fuer Astronomie, Koenigstuhl 17, 69117 Heidelberg, Germany 6    G. Chauvin 3Université Grenoble Alpes, IPAG, F-38000 Grenoble, France 3    C. Dominik 12Sterrenkundig Instituut Anton Pannekoek, Science Park 904, 1098 XH Amsterdam, The Netherlands 12    M. Feldt 6Max-Planck-Institut fuer Astronomie, Koenigstuhl 17, 69117 Heidelberg, Germany 6    R. Gratton 7INAF – Osservatorio Astronomico di Padova, Vicolo dell’Osservatorio 5, 35122 Padova, Italy. 7    J. Hagelberg 9Observatoire de Genève, Université de Genève, 51 Chemin des Maillettes, 1290, Versoix, Switzerland. 9    A. Isella 13Department of Physics & Astronomy, Rice University, 6100 Main Street, Houston, TX 77005, USA 13    M. Janson 14Department of Astronomy, Stockholm University, AlbaNova University Center, 106 91 Stockholm, Sweden 14    C. U. Keller 1 Leiden Observatory, Leiden University, P.O. Box 9513, 2300RA Leiden, The Netherlands.
1deboer@strw.leidenuniv.nl
   A.-M. Lagrange 3Université Grenoble Alpes, IPAG, F-38000 Grenoble, France 3    J. Lannier 3Université Grenoble Alpes, IPAG, F-38000 Grenoble, France 3    F. Menard 3Université Grenoble Alpes, IPAG, F-38000 Grenoble, France 3    D. Mesa 7INAF – Osservatorio Astronomico di Padova, Vicolo dell’Osservatorio 5, 35122 Padova, Italy. 7    D. Mouillet 3Université Grenoble Alpes, IPAG, F-38000 Grenoble, France 3    M. Mugrauer 15Astrophysical Institute and University Observatory Jena, Schillergäßchen 2, 07745 Jena, Germany 15    S. Peretti 9Observatoire de Genève, Université de Genève, 51 Chemin des Maillettes, 1290, Versoix, Switzerland. 9    C. Perrot 4LESIA, CNRS, Observatoire de Paris, Université Paris Diderot, UPMC, 5 place J. Janssen, 92190 Meudon, France. 4    E. Sissa 7INAF – Osservatorio Astronomico di Padova, Vicolo dell’Osservatorio 5, 35122 Padova, Italy. 7    F. Snik 1 Leiden Observatory, Leiden University, P.O. Box 9513, 2300RA Leiden, The Netherlands.
1deboer@strw.leidenuniv.nl
   N. Vogt 16Instituto de Física y Astronomía, Facultad de Ciencias, Universidad de Valparaíso, Av. Gran Bretaña 1111, Playa Ancha, Valparaíso, Chile. 16    A. Zurlo 17Núcleo de Astronomía, Facultad de Ingeniería, Universidad Diego Portales, Av. Ejercito 441, Santiago, Chile 17 2 2    SPHERE Consortium
Received July 8, 2016; accepted September 22, 2016
Key Words.:
protoplanetary disks – planet-disk interactions – circumstellar matter – stars: pre-main sequence – panets and satellites: detection – planets and satellites: formation
Abstract

Context:The effects of a planet sculpting the disk from which it formed are most likely to be found in disks in transition between classical protoplanetary and debris disk. Recent direct imaging of transition disks has revealed structures such as dust rings, gaps, and spiral arms, but an unambiguous link between these structures and sculpting planets is yet to be found.

Aims: We search for signs of ongoing planet-disk interaction and study the distribution of small grains at the surface of the transition disk around RX J1615.3-3255 (RX J1615).

Methods: We observed RX J1615 with VLT/SPHERE: we obtained polarimetric imaging with ZIMPOL (-band) and IRDIS (); and IRDIS () dual-band imaging with simultaneous spatially resolved spectra with the IFS ().

Results: We image the disk for the first time in scattered light and detect two arcs, two rings, a gap and an inner disk with marginal evidence for an inner cavity. The shapes of the arcs suggest that they probably are segments of full rings. Ellipse fitting for the two rings and inner disk yield a disk inclination and find semi-major axes of (278 au), (196 au) and (56 au), respectively. We determine the scattering surface height above the midplane, based on the projected ring center offsets. Nine point sources are detected between 2.1 and 8.0 separation and considered as companion candidates. With NACO data we recover four of the nine point sources, which we determine not to be co-moving, and therefore unbound to the system.

Conclusions: We present the first detection of the transition disk of RX J1615 in scattered light. The height of the rings indicate limited flaring of the disk surface, which enables partial self-shadowing in the disk. The outermost arc either traces the bottom of the disk or it is another ring with semi-major axis (435 au). We explore both scenarios, extrapolating the complete shape of the feature, which will allow to distinguish between the two in future observations. The most interesting scenario, where the arc traces the bottom of the outer ring, requires the disk truncated at  au. The closest companion candidate, if indeed orbiting the disk at  au, would then be the most likely cause for such truncation. This companion candidate, as well as the remaining four, require follow up observations to determine if they are bound to the system.

Date Instrument Mode Coronagraph filter m) FWHM (nm) DIT (s) (min) Seeing () SR(%) FWHM (mas) 12-05-2015 IRDIS IRDIFS ALC_YJH_S 1.59 & 1.67 53 & 55 64 136.5 0.65 65-70 44.5 & 46.6 12-05-2015 IFS IRDIFS ALC_YJH_S 0.96 - 1.34 55.1 64 136.5 0.65 45-65 30.0 - 36.1 15-05-2015 IRDIS IRDIFS ALC_YJH_S 1.59 & 1.67 53 & 55 64 68.3 0.65 65-70 44.5 & 46.6 15-05-2015 IFS IRDIFS ALC_YJH_S 0.96 - 1.34 55.1 64 68.3 0.65 45-65 30.0 - 36.1 06-06-2015 IRDIS DPI ALC_YJH_S 1.26 197 64 76.8 1.3 30-35 42.5 09-06-2015 ZIMPOL P2 - 0.626 149 120 96 1.0 95.5
Table 1: SPHERE observations of RX J1615. Full Width at Half Maximum of the filter transmission. For the IFS, the FWHM is given per spectral channel. Strehl Ratios are measured at the observing wavelength using the unsaturated PSFs (Girard et al. in prep.). For the low level of AO correction obtained for ZIMPOL, the FWHM is the appropriate image quality metric. Full Width at Half Maximum of the PSF.

1 Introduction

The evolution of circumstellar disks around pre-main sequence stars, and planet-disk interactions constitute two of the major components in our understanding of planet formation. Determining the disk structure (such as its geometry, gas and dust distribution) is essential to our understanding of both disk evolution and planet-disk interactions. Direct imaging of protoplanetary disks shows a variety of features, such as inner cavities, gaps, rings and spirals that can be observed both in early (HL Tau;  Myr, ALMA Partnership et al. 2015) and late evolutionary phases (e.g HD100546,  Myr, Grady et al. 2001). This poses the questions on whether the shapes and features of protoplanetary disks are caused by the presence of planets, or if the disk structures are regulating the formation of planets, or both. If the disk structures are not created by nearby planets, we still need to determine what did create these structures.

Strom et al. (1989) classified a subset of protoplanetary disks by a dip in the InfraRed (IR) range of their Spectral Energy Distribution (SED), caused by an inner dust depleted region (cavity). The authors suggested these disks represent an evolutionary stage ‘in transition’ between classical (younger) protoplanetary disks and more evolved debris disks. Multiple explanations have been suggested for the origin of cavities within transition disks, such as photo-evaporation (e.g. Hollenbach & Gorti 2005), and clearing by massive planets (e.g. Strom et al. 1989; Alexander & Armitage 2009; Pinilla et al. 2015). Transition disks also display other potential tracers of planet-disk interaction: the large scale structures in disks such as spiral arms in e.g. MWC 758 (Grady et al. 2013; Benisty et al. 2015); ring structures, as seen in e.g. RX J1604-2130 (Mayama et al. 2012; Pinilla et al. 2015) and TW Hydrae (Andrews et al. 2016, van Boekel et al. submitted;); and dust trapping vortices in e.g. IRS 48 (van der Marel et al. 2013).

RX J1615.3-3255, or 2MASS J16152023-3255051 (hereafter RX J1615) has previously been determined to be a 1.1 M pre-main sequence star with an age of 1.4 Myr (Wahhaj et al. 2010), and a member of the Lupus star forming region at a distance d  pc (Krautter et al. 1997). From its very strong H emission Wahhaj et al. (2010) determine that RX J1615 is a Classical T-Tauri Star (CTTS), and classify this target as having an IR-excess with a ‘Turn-on wavelength () in the range (3.6 - 8 m) of the IR Array Camera’ (TIRAC, where turn-on is the smallest with disk emission distinguishable from stellar emission). The authors describe TIRAC targets as young ( Myr) objects with a small inner clearing in the disks, i.e. transition disks. Andrews et al. (2011) have resolved the disk with the SMA interferometer at 880m, determined a disk Position Angle PA , inclination and found in their best-fit model a low-density cavity in the dust disk within  au. The disk is resolved with ALMA (cycle 0) in CO (6-5) and 690 GHz continuum data by van der Marel et al. (2015). The authors find PA and from their model a dust cavity inside  au. Both Andrews et al. (2011) and van der Marel et al. (2015) have used a gas-to-dust mass ratio of 100:1 in their models for the outer disk, and find a total disk mass of  M and  M, respectively.

While sub-mm observations probe emission of the mm-sized grains in the deeper layers of the disk, the surface of the micron sized dust disk can be detected by its scattering of Near-IR starlight. The large contrast between starlight and disk surface brightness at NIR wavelengths requires us to remove the stellar speckle halo, in order to detect the light scattered by the disk.

Since 2014 the Spectro-Polarimetric High-contrast Exoplanet REsearch (SPHERE, Beuzit et al. 2008) instrument has been commissioned at the Very Large Telescope (VLT). SPHERE is an extreme adaptive optics (AO) assisted instrument designed for high-contrast imaging of young giant exoplanets and circumstellar disks.

We use SPHERE to image the disk of RX J1615 in scattered light at multiple wavelengths with the aim of constraining the 3D disk geometry. In addition, we perform a search for possible companions that could be responsible for sculpting the disk. We report the new observations and the data reduction in Sections 2 and 3, respectively. Section 4 describes the new detections of the disk structures and companion candidates. In Section 5, we give constraints on the disk geometry and discuss possible scenarios for the disk vertical structure. We list our conclusions in Section 6. Based on archival data, we study the stellar properties of RX J1615 in Appendix A.

2 Observations

2.1 VLT/SPHERE IRDIS, IFS and ZIMPOL

We observed RX J1615 at multiple wavelengths using different modes of SPHERE (listed in Table 1). The extreme AO system, SAXO (SPHERE AO for eXoplanet Observation, Fusco et al. 2014) includes a 4141-actuator deformable mirror, pupil stabilization, differential tip tilt control and stress polished toric mirrors (Hugot et al. 2012) for beam transportation to the coronagraphs (Boccaletti et al. 2008; Martinez et al. 2009) and science instruments. The latter comprise a near-InfraRed Dual-band Imager and Spectrograph (IRDIS, Dohlen et al. 2008), a near-infrared Integral Field Spectrometer (IFS, Claudi et al. 2008), and the Zurich IMaging POLarimeter (ZIMPOL, Thalmann et al. 2008).

The SpHere INfrared survey for Exoplanets (SHINE), executed during the Guaranteed Time Observations (GTO) of the SPHERE consortium, includes observations of RX J1615 on May 12 and 15 of 2015. The observations were performed with the IRDIFS mode (Zurlo et al. 2014) in pupil tracking, using IRDIS in dual-band imaging mode (Vigan et al. 2010) in the -bands and IFS in -band. The observations are taken with the Apodized pupil Lyot Coronagraph ALC_YJH_S, which has a diameter of 185 mas. During the observations of May 12 the field rotated by . On May 15, we observed the system with a field rotation of .

On June 6, 2015, the SPHERE-Disk GTO program has observed RX J1615 with IRDIS in Dual-band Polarimetric Imaging mode (DPI, Langlois et al. 2014, de Boer et al., in prep.) in -band with the ALC_YJH_S coronagraph. For the six polarimetric cycles, we observed with the Half Wave Plate (HWP) at the angles , and to modulate the linear Stokes parameters ( and ). Per we have taken three exposures.

During open time observations on June 9, 2015 we have observed RX J1615 with SPHERE’s visible light polarimeter, ZIMPOL in field tracking (P2) non-coronagraphic mode in -band with the dichroic beamsplitter. We have recorded six polarimetric (or HWP) cycles, with ,, , and to modulate and . For each two frames were recorded for each of the two detectors of ZIMPOL.

The faint guide star ( mag, Makarov 2007) poses a challenge for the AO system. Especially in the visible light which is split between ZIMPOL and the SAXO wave front sensor, the moderate seeing ( 1.0) yielded very poor Strehl ratio SR and FWHM mas.

2.2 VLT/NACO and Keck/NIRC2 SAM

Additionally, RX J1615 was observed before with VLT/NACO (Lenzen et al. 2003; Rousset et al. 2003) as part of a survey for substellar companions in the Lupus star forming region during three epochs on May 7th 2010, May 8th 2011 and August 7th 2012 (Mugrauer et al. in prep.). The system was imaged in the -band (m,  nm) in standard jitter imaging (field stabilized) mode. NACO’s cube mode was used to save each individual frame (Girard et al. 2010). The Detector Integration Time (DIT) for a single exposure was 1 s for the three epochs

To investigate the possibility of stellar or high-mass substellar companions at close separations, we include the results from archival Keck NIRC2 Sparse Aperture Masking (SAM) data. RX J1615 was observed on April 4, 2012 with  s, 4 coadds per frame, 12 frames in total; on July 8, 2012 with  s, 4 coadds per frame, 62 frames in total; and on June 10, 2014 with  s, 2 coadds per frame, 40 frames in total. During all three epochs, the -filter and the nine hole mask have been used.

3 Data Reduction

3.1 ZIMPOL P2 observations in -band

Figure 1: a: ZIMPOL R’-band Q image, smoothed with a boxcar of 3 pixels. No radial scaling applied b: IRDIS/DPI J-band Q image, without radial scaling. c: IRDIS/DPI J-band Q image, with inclination () corrected scaling. This panel has the following disk features annotated: The southern ansa of a ring (), a full ring , a clear depression/gap in surface density () and an inner disk structure . Panels d,e, and f display the images corresponding to and with the same dynamic range as the images of panels a, b, and c respectively.

The Polarimetric Differential Imaging (PDI, Kuhn et al. 2001) reduction of the ZIMPOL data is performed according to the description of de Boer et al. (submitted) which is briefly summarized below. After dark subtraction and flat-fielding, we separate the orthogonal polarization states from each frame, and bin the resulting images to a pixel scale of  mas. For the first , we obtain the intensity image by adding the two polarization states per frame, and sum over the two frames (which measure the ‘’ and ‘’ phase, Schmid et al. 2012). We compute by subtracting the two polarization states per frame and subtracting the residual for the -phase (2nd frame) from the residual image of the -phase (1st frame). In the same manner, we measure and with ; and with ; and and with .

We center each Stokes I image by determining the offset between the image and a Moffat function using a cross-correlation. We use this offset to shift the image to the center and apply the same shift to the corresponding and images, and compute for both detectors (d1/d2):

(1)
(2)

Because we used the same () filter for both ZIMPOL detectors, we get the final , and images according to:

(3)
(4)

where the minus sign is necessary because the polarizing beamsplitter yields orthogonal polarization signals for the two detectors. We correct the and images for remaining instrumental and background polarization according to Canovas et al. (2011).

Finally, we compute the azimuthal Stokes components (cf. Schmid et al. 2006, where the azimuth is defined with respect to the star-center):

(5)
(6)

For cases of single scattering or simple symmetries, should contain all the signal, while provides an indication of the measurement error. The and images are smoothed with a boxcar of three pixel width and shown in Figures 1a & 1d, respectively.

Figure 2: TLOCI reductions of SPHERE-IRDIS ADI H2H3 (mean combination of the H2 and H3 filters) data. a: May 12 2015; b: May 15 2015; c: Mean combination of May 12 and May 15, with features annotated. From the outside in, we see a point source, which we consider as a companion candidate (cc1), an arc (), two full rings ( & ), another arc (), the location of the gap seen in Figure 1c (), and an innermost disk structure (I1); d: Same as b, for a larger dynamic range.
Figure 3: TLOCI reductions of SPHERE-IFS ADI Y-J data (0.96-1.07 m) a: Median of 13 channels (0.96-1.07 m). b: Median of 13 channels (1.08-1.21 m). c: Median of 13 channels (1.22-1.33 m). c: The median over the entire Y-J range of the IFS.
Figure 4: Companion candidates detected in NACO and SPHERE data. a: NACO image in K-band of May 7, 2010. We self-subtracted the image after it was rotated by . A small section of is seen west of the star. b: SPHERE-IRDIS cADI in H23-band of May 15, 2015. The green squares in a and b show the 4 point sources first detected with NACO; the white squares highlight the point sources detected with SPHERE-IRDIS only.

3.2 IRDIS/DPI observations in -band

The IRDIS DPI mode splits the beam before two orthogonal polarizers create the two orthogonal images simultaneously on different halves of the detector: the ordinary intensity image and the extraordinary intensity image . After we subtract the dark, and flat-field the images, we center the image of both detector halves on a Moffat function, using a cross-correlation. This first centering does not guarantee that the star is placed exactly at the center yet, but does ensure that all images of the star are at the same place with respect to the center, which suffices to subtract the two orthogonal polarization states on each frame. For the first is computed by subtracting the two simultaneouly measured beams:

(7)

By performing this operation for each of the three frames observed per we obtain three ; ;; and images for and , respectively.

We stack the three difference images obtained per and apply for each of the six HWP cycles the double difference:

(8)
(9)

For all HWP cycles we apply a correction on the and images to remove a detector artefact which creates continuous vertical bands on the IRDIS detector which vary with time. This correction is similar to the correction of Avenhaus et al. (2014) for comparable artefacts on the NACO detector: for each pixel column we take the median over the top 20 and bottom 20 pixels and subtract this signal from the entire column. Next, we compute and according to Equations 5 and 6.

An additional centering per HWP cycle uses a minimization of the signal in the image, which is based on the assumption that no astrophysical signal is measured in the image, ony noise (including reduction artefacts). We shift the and image at subpixel steps in the x and y direction and compute for each step. The shift which has the image with the lowest value over a centered but co-shifted annulus ( pixels) is assumed to place the star closest to the center of the image.

After we stack the centered and images, we correct for instrumental and sky polarization, as we did for the ZIMPOL reduction. From these stacked images, we compute the final and images for the -band, which are shown at the same intensity scale in Figures 1b and 1e, respectively. Figures 1c and 1d show these and images after scaling with the inclination-corrected radius squared, for which we used from van der Marel et al. (2015).

3.3 Irdis () and IFS () pupil tracking observations

The IRDIS and IFS data are reduced using the SPHERE Data Center. We use the SPHERE pipeline (Pavlov et al. 2008) to process cosmetic reductions including sky subtraction, bad pixel removal, flat field and distortion corrections, IFS wavelength calibration, IFU flat correction, instrument anamorphism correction (0.60  0.02%, Maire et al. 2016), and frame registering.

Then, several types of Angular Differential Imaging (ADI, Marois et al. 2006) based algorithms are implemented in a dedicated tool (SpeCal, Galicher et al., in prep.) to perform starlight subtraction independently for each of the two IRDIS filters, and the 39 IFS spectral chanels.

For IRDIS, we present the results obtained with classical ADI (cADI) for the full field of view (Figure 4b) and the Template Locally Optimized Combination of Images algorithm (TLOCI, Marois et al. 2014) for a 4x4 field (Figure 2). The May 12 data shows a strong vertical negative residual in the reduction, for both IRDIS (Figure 2a) and IFS. We ascribe this to an aberration of the PSF, due to a loss of AO performance when coping with high altitude winds. The May 15 data is not plagued by this effect. Figure 3 shows the IFS TLOCI+ADI reductions for May 15, where the first three panels show the median combination of 13 spectral chanels between 0.96-1.07 m, 1.08-1.21 m, and 1.22-1.33 m, respectively. The righthand panel shows the median combination of all 39 spectral bands.

3.4 Naco jitter imaging observations

The NACO data reduction is similar for all epochs. First, we average all frames in each data cube, and use the jitter routine in the ESO Eclipse package (Devillard 2001) to flat-field, shift and combine all averaged images. An initial offset of each frame is determined from the image header and is then refined using jitter’s cross correlation routine. Since the star was moved to different positions on the detector, each pair of consecutive frames can be used to estimate the sky background in -band and subtract it. To remove the stellar halo, we stack the centered image, and subtract the same image after we have rotated it with . We present the resulting differential image in Figure 4a.

The astrometric calibration of the NACO epochs was taken from Ginski et al. (2014). They imaged the core of the globular cluster 47 Tuc for this purpose in the same filter and imaging mode. These calibrations are within two days of the RX J1615 observations for the 2010 and the 2012 epoch. Due to bad conditions during the 2011 observations, no companion candidates were recovered.

3.5 Keck NIRC2 SAM observations

Data were reduced using the aperture masking pipeline developed at the University of Sydney. An in-depth description of the reduction process can be found in Tuthill et al. (2000) and Kraus et al. (2008), but a brief summary follows: data were dark subtracted, flat-fielded, cleaned of bad pixels and cosmic rays, then windowed with a super-Gaussian function. The complex visibilities were extracted from the cleaned cubes and turned into closure phases. The closure phases were then calibrated by subtracting a weighted average of those measured on several point-source calibrator stars observed during the same night.

4 Results

4.1 Disk

The disk of RX J1615 is detected in all the datasets included in this study. For the ADI as well as the PDI images, all features described below are brighter on the northeastern side from the major axis.

- Outer rings: In the -band scaled Q image (Figure 1c) and the (not scaled) TLOCI ADI images of the disk, we see multiple arc-like and continuous ellipses which we consider to be ring-like features in the scattering () surface of the disk, projected with the inclination of the system. For decreasing separation, Figure 2c shows:

  • an arc ();

  • two full rings and ;

  • a second arc ().

Where the two arcs cross the minor axis of the disk, they seem to lie parallel to both rings. We therefore consider it most likely that both arcs are segments of full rings. The feature does not appear to be as clearly separated from in the -band scaled polarization image, and is therefore not annotated in the image. However, in this image the disk signal at radii appears continuous in the azimuthal direction, which confirms that is indeed a full ring. The disk feature is detected (at least as a ring segment) in all datasets except the ZIMPOL -band polarization image. In Figure 3, we see that the segment which lies within the IFS field of view is detected more clearly at longer wavelengths. Only the southern ansa of is detected in the -band Q image. In Figures 1b,c and 2 we can discern that and are clearly not concentric with the inner disk component and the star-center.

- Gap: A gap (feature ) in between features and is detected, most clearly in Figure 1c. Figure 2c shows the gap as well, but not for all azimuth angles.

- Inner disk: Closest to the star, we see an elliptical inner disk component (feature in Figure 2b). The surface brightness is continuous inward in the Q images in both and -band (Figures 1a & 1c). However, for a disk with continuous surface density with a linearly increasing scattering surface we would expect the surface brightness to drop off with the distance to the star squared. Conversely, if we create an image with the surface brightness scaled with (corrected for the inclination), the previous example of a continuous surface density would show a continuous surface brightness. However, in our inclination-corrected -band image (Figure 1c) appears more ring-like, with a cavity inside, which agrees with the outer radius of the inner dust cavity as determined from the 880 m interferometric observations of Andrews et al. (2011,  au). Still, we consider the detection of the cavity as tentative because is bordering the coronagraph. We also see a non-continuous surface brightness of in the ADI reductions of Figure 2. However, observations of a disk with continuous surface brightess is likely to be plagued by self-subtraction, which is often seen in ADI for low and intermediate inclination circumstellar disks (Milli et al. 2012). This kind of self-subtraction does not occur in PDI, which is extremely efficient at isolating the polarized disk signal, hence revealing the disk structure with high fidelity. Note that based on our results we cannot rule out a disk component at a separation  au (inside ). However, Andrews et al. (2011) mention that their model requires a very low density inside  au to adequately fit the SED.

Ellipse fitting

Parameter ADI- PDI- Semi-major axis () 1.50 0.01 1.50 0.01 Semi-minor axis () 1.01 0.01 1.02 0.01 RA offset () -0.15 0.01 -0.15 0.01 Dec offset () -0.10 0.01 -0.10 0.01 Offset angle () 238 1 236 1 PA () 145.7 1.0 144.2 0.8 Inclination angle () 47.3 1.0 47.0 0.8 (au) 44.9 2.2 44.7 1.7 0.162 0.009 0.162 0.007 Semi-major axis () 1.06 0.01 1.06 0.01 Semi-minor axis () 0.70 0.01 0.72 0.01 RA offset () -0.10 0.01 -0.10 0.01 Dec offset () -0.07 0.01 -0.06 0.01 Offset angle () 235 1 236 2 PA () 145.4 1.3 144.3 1.4 Inclination angle () 48.5 1.3 46.8 1.4 (au) 30.9 2.4 29.6 2.2 0.158 0.014 0.152 0.013 Semi-major axis () 0.30 0.01 0.35 0.01 Semi-minor axis () 0.20 0.01 0.24 0.01 RA offset () -0.01 0.01 0.00 0.01 Dec offset () 0.00 0.01 0.00 0.01 Offset angle () 261 209 PA () 145.5 4.2 144.5 4.3 Inclination angle () 49.0 3.9 47.7 4.1 (au) 3.5 1.4
Table 2: Ellipse parameters for the fits to the features listed in Figure 5a and 5b. ) For , the errors on the ellipse offset () are larger than the measured value. We therefore do not consider the value for the ellipse offset angle and to be significant.
Figure 5: a: IRDIS -band PDI image after inclination-corrected scaling. b: IRDIS -band TLOCI + ADI image. The overplotted green ellipses are the fits to , and , as listed in Table 2.
Figure 6: Schematic view of a double ringed thick disk. Observed at inclination , circular rings will appear elliptical. The height of the scattering surface of Ring 1 (()) can be found by determining the semi-minor axis (), and the distance between the star and the center of the ellipse (). Note that this figure has the purpose of explaining how to determine for the different rings, not to portray ideas about the radial density distribution.

For both the -band PDI image and the -band TLOCI ADI image we have fitted ellipses111 using the routine mpfitellipse.pro with the Interactive Data Language (IDL) to the two rings and to the inner disk. The resulting ellipses are overplotted in green in Figures 5a and 5b, and the ellipse parameters are listed in Table 2. From the assumption that the ellipses trace the highest surface above the midplane for circular rings (i.e. not the slope facing the star/wall), we determine the inclination of each ellipse according to semi-minor axissemi-major axis. With a weighted mean we find an inclination of the disk , where the error represents the random error on the fitted ellipses but does not include the systematic errors from our method. A conservative estimate of the systematic errors brings us to a final value of , which is in good agreement with the inclination derived by van der Marel et al. (2015). The weighted mean for the position angles is PA . Including our estimate of the systematic errors gives a final PA , in between the values of Andrews et al. (2011, PA ) and van der Marel et al. (2015, PA ).

We also find that and are not centered around the star (listed as X () and Y () offset in Table 2). The offset of the ellipse-centers with respect to the position of the stars means that the rings are either:

  • eccentric rings;

  • concentric circular disk components with considerable radially increasing thickness, as described by the height of the surface () above the disk midplane;

  • an intermediate of the two extremes: an eccentric and a thick ring.

The directions in which and are offset from the star (listed as ‘Offset angle’ in Table 2) are roughly PA , i.e. along the minor axes of the ellipses. These apparent ellipse displacements along the minor axes is a necessary condition for the offsets to be caused by a projection of the surface on the midplane due to the inclination of the system. Below, we explore the scenario of thick circular rings viewed at an inclination away from face-on. In this scenario, the disk in the direction opposite to the ellipse offsets (i.e. the eastern side, PA ) forms the near side. Without trying to explain the surface density distribution of the disk, in Figure 6 we portray such a disk with thick rings. From the sketch we can derive that for any given ellipse with its center offset from the star-center with distance , we can determine the height of the scattering surface () of this ring according to:

(10)
(11)

In a similar fashion, Lagage et al. (2006) have determined the tickness of HD 97048 based on PAH emission maps. The main difference between the method used by Lagage et al. (2006) and our method is that the former have used isophotes, while we can use the sharp rings in the disk of RX J1615. It is clearly visible in Figure 5 that the surface brightness of the rings is strongly varying with azimuth angle. To accurately determine based on isophotes it is necessary to correct for any azimuthal surface brightness variation, which requires radiative transfer modeling. To the best of our knowledge, this is the first example where can be determined strictly from geometrical constraints (i.e. model independent) in scattered light images.

Polarized intensity profile

Figure 7: intensity profile of scaled -band image (Figure 1c) along PA . The green dashed lines mark the semi-major axes of the ellipses , and . Note that although this profile is measured for a PA similar to those of the ellipses, the profile does include the star center, and therefore does not lie on the major axes of and . This explains why the peaks of the rings lie at closer separations than their semi-major axes.

Figure 7 shows the intensity profile of the scaled image in -band (Figure 1c), after smoothing with 4 pixels. The profile is measured along PA , centered on the star (i.e. parallel to but not on the major axes of and ). The errors are the standard deviation measured over the same aperture in the image (Figure 1f), divided by the square root of the number of pixels. The disk features which can be distinguished in the profile are listed along the top axis of the plot. Due to the ellipse-center offset with respect to the star, the profile does not cross the ansae of the ellipses, which becomes visible when we compare the position of peaks in the intensity profile with the semi-major axis of the ellipses for and , annotated with green dashed lines (values adopted from Table 2). Although the feature is not clearly visible in Figure 1c, we can recognize the two peaks of in the the intensity profile. The gap is visible at  au. van der Marel et al. (2015) find a tentative dust gap between 110 - 130 au (0.6 - 0.7) in the 690 GHz continuum profile. This sub-mm gap lies outside the gap we detect in the -band data.

4.2 Point sources

Figure 8: Astrometry for the 4 point sources detected (solid points) by both NACO (black and green) and IRDIS (red) and the expected values for background objects (open points). Data points are color coded by date, not by instrument. These point sources are marked with the green squares in Figure 4. The numbers of the candidate companions listed in the title of each panel column are the same as shown in Figure 4. All 4 point sources are shown to not be co-moving with RX J1615.

Four point sources have been detected in the NACO -band image, and marked with green boxes in Figure 4a & 4b. In the IRDIS data (Figure 4b), we detect no point sources within the ringed disk structure, out to . We do detect nine point sources outside the outermost disk component , including all four point sources detected with NACO. The five new companion candidates are marked with white boxes in Figure 4b. For reference, we included the white boxes at the same locations in Figure 4a, even though these point sources are not detected by NACO. The astrometry and photometry of the six innermost candidates detected with IRDIS were derived using the LAM-ADI pipeline (Vigan et al. 2012, 2016) using injection of fake negative planets in the pre-processed ADI data cubes. The position and flux of the negative fake companion are adjusted using a Levenberg-Marquardt least-squares minimisation routine where we try to minimise the residual noise after ADI processing in a circular aperture of radius centered on the position of the companion candidate. The error bars for the fitting process are then calculated by varying the position and contrast of the fake companion until the variation of the reduced reaches a level of 1. For the speckle subtraction, we have used a PCA analysis (Soummer et al. 2012) where five modes are subtracted. At the distances of the outer three candidates there is negligible speckle noise negating the need for such post-processing. For these we perform basic relative photometry after subtracting a radial profile and stacking the de-rotated frames. The astrometry and photometry measurements for the companion candidates are reported in Table 3. For IRDIS, we have adopted a true north direction of and a plate scale of  mas/pix, based on astrometric calibrations performed on the theta Ori B field (Maire et al. 2016, and Vigan private communication). For the NACO observations of 2010, we adopted a true-north of and a plate scale of  mas; for 2012 we used a true-north of and a plate scale of  mas. The final error bars for the photometry include the fitting error detailed above; the uncertainties on the star center (1 mas); and the level of noise residuals estimated at the same separation as the detections.

Obj Date Instr. filt. cc1 15/05/2015 IRDIS H2 15/05/2015 IRDIS H3 cc2 05/04/2010 NACO Ks 07/05/2012 NACO Ks 15/05/2015 IRDIS H2 15/05/2015 IRDIS H3 cc3 15/05/2015 IRDIS H2 15/05/2015 IRDIS H3 cc4 05/04/2010 NACO Ks 15/05/2015 IRDIS H2 15/05/2015 IRDIS H3 cc5 15/05/2015 IRDIS H2 15/05/2015 IRDIS H3 cc6 15/05/2015 IRDIS H2 15/05/2015 IRDIS H3 cc7 05/04/2010 NACO Ks 07/05/2012 NACO Ks 15/05/2015 IRDIS H2 15/05/2015 IRDIS H3 cc8 15/05/2015 IRDIS H2 15/05/2015 IRDIS H3 cc9 05/04/2010 NACO Ks 07/05/2012 NACO Ks 15/05/2015 IRDIS H2 15/05/2015 IRDIS H3
Table 3: Astrometry and photometry for nine companion candidates relative to RX J1615

In Figure 8 we compare the astrometry for the four point sources detected with NACO with that of our IRDIS detection. We determine them to not be co-moving, and therefore not associated to the RX J1615 system. The five remaining candidates, seen in Figure 4b, require follow up observations to determine whether they are bound to RX J1615.

No significant point-source signals were found in any of the Keck NIRC2 SAM datasets, and so we place limits on their detectability by drawing 10,000 simulated closure phase datasets consistent with Gaussian random noise using the measured uncertainties. For each combination of separation, contrast and PA on a 3D grid, the simulated datasets were compared to a binary model. The () detection limits were calculated as the point at which the binary model gave a worse fit to 99.9% of the simulated datasets. The datasets probe separation ranges as small as 30 mas, and achieve contrast limits that are approximately flat at larger separations. The 2014-06-10 data allow us to rule out objects with  mag (  M for a system age of 1.4 Myr), while the 2012-07-08 and 2012-04-14 data reach similar contrasts of  mag and  mag, respectively.

5 Discussion

5.1 Disk geometry

The explanation for the apparent ellipse offsets by the projection of the surface at height above the midplane, as given in Section 4.1, implies that the northeast (PA ) is the near side of the disk (i.e. pointing towards earth). Min et al. (2012) and Dong et al. (2016) show that the predominance of forward scattering in the near side of disks gives it in general a larger surface brightness in total intensity than the far side. The ADI intensity image of Figure 2 shows the northeastern sides of , and to be brighter thant their southwestern sides, confirming that the northeast is the near side of the disk.

In Section 4.1, we find the disk inclination , which is in good agreement with both the disk inclination derived by van der Marel et al. (2015) and the inclination of the stellar rotation axis derived in Appendix A from  km s (Wichmann et al. 1999, where the stellar rotational velocity depends on the determination of the stellar radius from its luminosity, which in turn depends on the interstellar extinction  mag, as suggested by Manara et al. (2014)). Higher extinction (e.g.  mag, Andrews et al. 2011) would lead to lower inclinations of the stellar rotation axis, which do not agree with the inclination of the disk. However, the disk inclination would be altered if the observed ellipses do not trace the rings’ highes surfaces above the midplane (called “edge ring 1/2” in Figure 6). Without advanced radiative transfer modeling it is not possible to differentiate in our scattered light images between starlight scattered off the ring edges and light scattered off the slopes/walls of the rings (called “ring slope/wall” in Figure 6). Consequently, the ring edges might truly lie further out than our ellipse fits. This effect would not be symmetric; we are more likely to have a larger contribution of the wall in the backward scattering side (southwest, near the minor axis) of the rings. Therefore, if the apparent ellipses are affected by scattering by the ring walls, the ellipticity of the true ring edges (and their inclination) will be slightly lower and the ring offsets larger.

5.1.1 Nature of the ring structures

It is tempting to interpret the rings and arcs detected in the disk as spatial variations in surface density. However, our NIR scattered light detections of the disk trace the disk surface for the micron-sized dust grains. We cannot unambiguously determine whether the rings are either a manifestation of variations in the scale height caused by spatial variations in temperature (e.g. due to shocks), or variations in the surface density either caused by dead zones or by massive planets carving a gap in the gas surface density. Dust trapping by local peaks in the gas pressure will appear different for small and large dust grains. Pinilla et al. (submitted) show that if we can measure a difference between the density enhancements for the different grain sizes, we can discriminate between dead zones creating a bump in the gas pressure and massive planets carving a gap in the gas disk. Furthermore, the mass of a gap carving planet can be predicted by the amount of the displacement between large grain and small grain peaks in the surface density for a given gas viscosity (de Juan Ovelar et al. 2016). Large baseline sub-mm (ALMA) observations are therefore required to study the origin of the ring structures. A resolution comparable to our SPHERE observations will be needed in order to resolve the ring structure (when present for large grains), and accurately compare their radius with those of the rings in the surface of the small grain dust disk presented in this study.

Another asymmetry is detected in the polarized intensity profile of Figure 7, which roughly traces the major axes of the ellipses. The peaks of and are brighter in the northwest than their southeastern counterparts, while the opposite (brighter in the southeast) is true for and . The fact that these brightness asymmetries along the rings are oscillating between northeast and southwest can possibly be explained with shadowing: the brighter parts of the rings might have a larger scale height than their fainter counterparts. If each ring outside of is just marginally rising out from the shadow of the ring directly inside, this would cause a brighter ring segment in the inner ring to cast a larger shadow on the outer ring, hindering the stellar radiation to heat up the segment in the outer ring with similar azimuth angle. The opposite happens for the faint segments, which have a smaller scale height, casting less of a shadow on the next ring: this next ring receives more stellar radiation, heating it up and allowing it to ‘puff up’ more, making it brighter than the opposite side of the same ring. Shadowing of the outer rings by the inner rings is only possible when the flaring of the disk is very small. When we compare the disk angle (parameter in Figure 6) for both rings, we find that is marginally larger than , which confirms that the disk flaring is minimal.

Although ALMA is detecting an increasing number of protoplanetary disks with multiple rings (e.g. HL Tau), very few have been detected in scattered light. To the best of our knowledge, only TW Hya (Rapson et al. 2015, van Boekel et al. submitted), HD 141569A (e.g.  Weinberger et al. 1999; Perrot et al. 2016) and HD 97048 (Ginski et al. 2016) display multiple rings in scattered light. The inclination of HD 141569A (between and for the different rings, Biller et al. 2015, and for the entire disk, Mazoyer et al. 2016) is comparable to RX J1615. The ringed structure also looks very similar to RX J1615, because the rings are relatively sharp compared to the larger radial extent of its gaps. As we discussed above, in a disk with low flaring (hereafter ‘flat’, i.e.  constant), a small ripple in the disk surface can cast large shadows outward. Indeed, Thi et al. (2014) suggest that the disk of HD 141569A is very flat, while the radial extent of rings in the surface of more flaring disks, such as HD 97048 and TW Hya, is of similar size or larger than the width of their gaps. We therefore suggest that the sharpness of rings in the surface of a primordial disk can be used as a tracer for the degree of flaring of the disk surface. This argument is based on the assumption that the apparent gaps are (mainly) due to shadows cast by ripples in the scattering surface, rather than true gaps in the surface density of the disk. To test this hypothesis we would need better knowledge of the scale height of the disk (e.g. through high angular resolution gas observations with ALMA) in several disks that show multiple ring structures in scattered light.

5.1.2 : additional ring or bottom of

Parameter Red Purple Semi-major axis () 1.66 2.35 (no fit) Semi-minor axis () 1.14 1.61 RA offset () 0.16 -0.23 Dec offset () 0.11 -0.16 Offset angle () 56 236 (au) 49.2 70.7 0.16 0.16
Table 4: Ellipse parameters for feature in Figure 2c. In Figure 5, the red ellipse shows the scenario where is the bottom side of ; the purple ellipse shows as a separate outermost ellipse.
Figure 9: b: IRDIS -band TLOCI + ADI image. The overplotted green ellipses (solid lines) are the fits to , and , as listed in Table 2, while the red (dashed) and purple (dash dot) ellipses depict the two scenarios where is either the backward-facing side of ring 1, or the forward-facing side of a 3rd ring further out (sketched in panels c and d), respectively. c: First explanation for , for the same considerations as Figure 6: is the bottom/backward facing side of Ring 1, for which is the top/forward facing side, i.e. . d: Second explanation for : is the top/forward facing side of a third ring, with . Both the red and the purple ellipses are selected based on the assumption that compared to and , would have roughly the same aspect ratio between major and minor axes as well as the same ratio between ellipse-star offset and the major axis.

For the arc-like structure () in Figure 2, we consider two explanations to be equally plausible: it could either be an additional ring, at a separation from the star () larger than for . However, the similarity of this arc to the shape of at similar PAs () is much stronger than when we compare and at these PAs. Therefore, we suggest an alternative explanation for as it being the backward facing end of , comparable to the bottom part of ring 1 in Figure 6.

To speculate on the shape of the ellipse for either explanation of , we added two ellipses (red and purple) to Figure 9b. The red and purple ellipses (explained in the cartoons of Figures 5c and 5d) are fixed at the intersection between and the minor axis of the disk; have the same PA as , and the same ratio of major/minor axes (i.e. inclination). The final constraints we used for the two ellipses is that

  • the ring has either a similar (or larger) for the scenario where is the backward facing side of (red ellipse in Figure 9b and 9c);

  • of a new ring further out (purple ellipse in Figure 9b and 9d) needs to be large enough to rise from the shadow of . As we discussed in Section 5.1, this means that .

From Equation 10 and the constraints above we derive that the absolute value of the ellipse (x-y) offset () divided by the minor axis () for should be:

(12)

We have used for both ellipses in Figure 9b, which can be considered as a best guess for the red ellipse and a lower limit for the purple ellipse: larger purple ellipses (which will have larger ) are not ruled out. Deeper observations of this structure will reveal which of the two scenarios is true: detecting a larger (azimuth coverage of the) ring segment will enable us to distinguish between the red and the purple ellipse scenarios.

In order to fix the red ellipse on the intersection between the minor axis and the detected feature, both a larger semi-major axis and X and Y offset (in opposite direction) are needed for the red ellipse (see Table 4) than what we found for the green ellipse of in Table 2. This either hints at a structure with or . The scattering-angle for the backward facing part of ring 1 will be different than the forward facing part: moving the beam back into denser regions of the disk instead of away from higher density (as for the forward facing ). This will bring the backward facing surface to lie further from the midplane than for the forward facing surface at the same distance .

5.2 Possible disk sculpting companion

Should the ring turn out to be the backward facing side of ring 1, as we illustrate in Figure 9c, we would expect there to be a fairly massive planetary companion to the system. To make the backward facing side of ring 1 visible, the disk needs to be truncated at most  au beyond ring 1, which lies at  au (where we used the semi-major axis of for , instead of the for , since the former will be dominated by the outer edge, the latter by the inner edge of ring 1). The disk truncation at  au can be a consequence of a massive planet beyond this radius, opening a gap in the disk. Detection of such a companion should be possible with our IRDIS ADI observations as long as the planet signal is not obscured by the disk (e.g. when the planet is on the western side of its orbit, or the truncation is not performed by multiple lower mass planets. Should companion candidate 1 (cc1, see Figure 2) be associated with the system it could provide the required disk truncation.

To determine the mass of cc1 (if it is indeed bound to the system) we first need to find the age of the system. In Appendix A we determine the rotation period ( d), spectral type (K5-K6) and temperature ( K) of RX J1615, which we use to determine the inclination of the stellar rotation axis and constrain the stellar age and mass. Overall, the properties of the star such as accretion, characteristics of the disk, the kinematic, and the limits from lithium and rotation period indicate an age less than 5 Myr. Following the same process as Wahhaj et al. (2010) we use the models of Siess et al. (2000) to determine the age of RX J1615 from isochrones. Adopting the luminosity and effective temperature , as stated before, we find age range of  Myr and a mass of  M.

Using the COND model (Baraffe et al. 2003), our observed -band magnitude, and the assumed age of the system, we estimate the mass of cc1 to be  M. In this estimation, we have assumed the signal to be dominated by thermal emission of the planet, rather than the emission by a circumplanetary accretion disk. If we assume that cc1 is a companion orbiting RX J1615 along the midplane of the disk, its deprojected distance to the star would be  au. Assuming that the maximum gap size created by a companion on a circular orbit is Hill Radii (Dodson-Robinson & Salyk 2011; Pinilla et al. 2012), at this separation we would need a planet with higher mass ( M) to fully truncate the gas disk at  au, or the planet needs to have an elliptical orbit. However, if the gas density is sufficiently reduced in the outer disk, in particular at the planet gap, the m sized grains can become decoupled to start drifting inwards, leading to a smaller disk truncation of the micron-sized particles compared to the gas. Without the presence of a planetary gap, the micron sized particles would remain coupled and the truncation in scattered NIR light would be further out than observed. Whether the potential drift of the m grains in a continuous gas disk suffices to reduce the disk optical depth enough to enable the transmission of light scattered by the backward-facing side of ring 1 requires radiative transfer modelling, which is beyond the scope of this study.

We show in Figure 10 the locations of our candidate companions on a color magnitude diagram, and included previously known brown dwarfs and planetary companions for reference. Upon first inspection of this diagram we might assume that cc1 is more likely to be a background object as it, along with all remaining candidates, has an color typical of a background M dwarf. However, we know that some planetary companions (such as those of HR8799 and 2M1707) do have colors within the range of cc1 so we can claim nothing for certain. Indeed RX J1615 is treading on new territory being so young and low mass. If cc1 is attributed to the RX J1615 system we would normally expect there to be a methane signature (which would yield a negative color). However, to date no observational evidence for methane emission in such young and low mass systems exists, leaving the possibility open that cc1 is bound to the RX J1615 system. Follow-up observations of cc1 are needed to confirm whether the candidate is co-moving with RX J1615. If bound cc1 would be an ideal target for characterization with the IRDIS low resolution long-slit spectroscopy mode (Vigan et al. 2008).

Figure 10: Color magnitude diagram displaying our candidate companions compared to already known planetary companions and brown dwarfs. Candidate companions are displayed, and numbered, in red with corresponding error bars.

Besides the option of the disk being truncated by a single massive planet as an explanation of why could possibly be the backward facing side of , multiple lower mass planets can provide a similar truncation of the disk. Such lower mass planets might remain undetectabel by SPHERE. For this reason we chose to include the contrast plot seen in Figure 11, showing the contrast in two cuts across the center of RX J1615, parallel to the major and along minor axes of the rings. A conventional contrast curve, where the contrast is measured over concentric annuli would contain the disk signal in such a way that one can no longer tell if a peak is caused by a ring or by something else. However, outwards of R1 this is not an issue. We have inserted the measured signal of cc1 at its appropriate radial separation for reference.

Figure 11: The contrast measured in the H23 TLOCI reduction of May 15 2015, along PA , parallel to the ellipse’s major axes (black solid line, negative is southeast); and PA , the ellipse’s minor axes (blue dash-dotted line, negative is northeast). The vertical green dashed lines display the semi-major axes of , , and as annotated above the top axis of the plot. The red diamond shows the measured signal of cc1 (which, we find, corresponds to ). Notice that the actual location of cc1 does not coincide with either the minor or the major axes of the disk.

5.3 Wavelength dependent surface brightness of

is clearly detected in -band ADI image; it is fainter in the -band PDI data; and can no longer be distinguished against the background disk signal in the -band PDI image. We detect a similar wavelength dependence of the surface brightness when we compare the three wavelength regimes of the IFS in Figure 3. If this wavelength dependence is truely astrophysical, it could be indicative of spatial/radial variations in the chemistry or grain size distribution throughout the disk. However, the low Strehl ratio and FWHM of the ZIMPOL data described in Section 2 can also be responsible for washing out a structure which is not resolved in the radial direction, smoothing out the unresolved structure so much that we can no longer distinguish from the surface brightness of its surroundings. The simultaneous measurements at the three IFS wavelength regimes suffer from the same systematic effect: the Strehl ratio increases for longer wavelength.

Since disk signal is present both directly inside and outside of , we cannot determine if the low resolution of the image has washed out the disk signal. Convolving a radiative transfer model with the PSFs of the different observations will be the best way to answer whether we are looking at a true astrophysical change with wavelength of the surface brightness or rather at systematic effects due to the difference in Strehl ratio and FWHM. However, creating a realistic radiative transfer model lies outside the scope of this paper and is left for future work.

6 Conclusion

We have studied the system of RX J1615 with four observing modes of the high contrast imager VLT/SPHERE. We detect the disk of RX J1615 in scattered light for the first time, from the optical -band to the NIR -band, with all high contrast imaging modes used in this study. The and band images show three elliptical disk components surrounding the star and two arc-like features. The two outer ellipses and have their centers offset from the position of the star. The simplest explanation for the elliptical features is that we see a disk at inclination which contains an inner disk surrounded by two circular rings, for which the height of the surface increases with its distance to the star. Ellipse fitting yields the major axes (i.e. radii (r) of the rings) for to be  au, for  au, and for  au. A tentative cavity is detected inside the inner disk. The apparent offsets between star and ring centers allow us to determine the height of the scattering surface above the midplane:  au and  au.

For the outermost disk feature , our detection was not deep enough to determine its origin. Deeper observations are needed to determine whether is the backward facing side of or a new ring at larger separation from the star. ALMA observations might allow us to directly image the thickness of the gas disk, as was done for HD163296 by de de Gregorio-Monsalvo et al. (2013). If we detect the dust disk truncated between  au, and the gas disk with comparable to the thickness we derived, this would be a strong confirmation of our understanding of the disk geometry.

Nine companion candidates are detected between 2.1 - 8.0 with SPHERE/IRDIS. In VLT/NACO data we detect four of these nine point sources, and determine that they are not co-moving, and therefore not bound to the system. Follow-up observations of the remaining five viable companion candidates need to be made to determine if they are co-moving. If cc1 indeed turns out to be orbiting the disk at  au and if is indeed showing the backward facing side of a truncated disk, RX J1615 provides the most unambiguous example of ongoing planet-disk interaction yet detected.

Acknowledgements.
We thank the anonymous referee for his/her very rapid and constructive comments. Many thanks go out to the instrument scientists and operators of the ESO Paranal observatory for their support during the observations. We also thank Michiel Min for the insightful discussion about the geometry of the disk. P. Pinilla is supported by a Royal Netherlands Academy of Arts and Sciences (KNAW) professor prize. AJ is supported by the DISCSIM project, grant agreement 341137 funded by the European Research Council under ERC-2013-ADG. MM and CG acknowledge the German science foundation for support in the programme MU2695/13-1. This research has made use of the SIMBAD database, operated at the CDS, Strasbourg, France. NASA’s Astrophysics Data System Bibliographic Services has been very useful for this research. 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 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 the CNRS (France), MPIA (Germany), INAF (Italy), FINES (Switzerland) and NOVA (Netherlands). SPHERE also received funding from the European Commission Sixth and Seventh Framework Programs 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).

References

  • Alexander & Armitage (2009) Alexander, R. D. & Armitage, P. J. 2009, ApJ, 704, 989
  • Allard et al. (2012) Allard, F., Homeier, D., & Freytag, B. 2012, Philosophical Transactions of the Royal Society of London Series A, 370, 2765
  • ALMA Partnership et al. (2015) ALMA Partnership, Brogan, C. L., Pérez, L. M., et al. 2015, ApJ, 808, L3
  • Andrews et al. (2011) Andrews, S. M., Wilner, D. J., Espaillat, C., et al. 2011, ApJ, 732, 42
  • Andrews et al. (2016) Andrews, S. M., Wilner, D. J., Zhu, Z., et al. 2016, ApJ, 820, L40
  • Avenhaus et al. (2014) Avenhaus, H., Quanz, S. P., Schmid, H. M., et al. 2014, The Astrophysical Journal, 781, 87
  • Baraffe et al. (2003) Baraffe, I., Chabrier, G., Barman, T. S., Allard, F., & Hauschildt, P. H. 2003, A&A, 402, 701
  • Benisty et al. (2015) Benisty, M., Juhasz, A., Boccaletti, A., et al. 2015, A&A, 578, L6
  • Beuzit et al. (2008) Beuzit, J.-L., Feldt, M., Dohlen, K., et al. 2008, in Proc. SPIE, Vol. 7014, Ground-based and Airborne Instrumentation for Astronomy II, 701418
  • Bianchi et al. (2011) Bianchi, L., Herald, J., Efremova, B., et al. 2011, Ap&SS, 335, 161
  • Biller et al. (2015) Biller, B. A., Liu, M. C., Rice, K., et al. 2015, MNRAS, 450, 4446
  • Boccaletti et al. (2008) Boccaletti, A., Abe, L., Baudrand, J., et al. 2008, in Proc. SPIE, Vol. 7015, Adaptive Optics Systems, 70151B
  • Butters et al. (2010) Butters, O. W., West, R. G., Anderson, D. R., et al. 2010, A&A, 520, L10
  • Canovas et al. (2011) Canovas, H., Rodenhuis, M., Jeffers, S. V., Min, M., & Keller, C. U. 2011, Astronomy & Astrophysics, 531, A102
  • Claudi et al. (2008) Claudi, R. U., Turatto, M., Gratton, R. G., et al. 2008, in Proc. SPIE, Vol. 7014, Ground-based and Airborne Instrumentation for Astronomy II, 70143E
  • Cutri & et al. (2013) Cutri, R. M. & et al. 2013, VizieR Online Data Catalog, 2328
  • Cutri et al. (2003) Cutri, R. M., Skrutskie, M. F., van Dyk, S., et al. 2003, 2MASS All Sky Catalog of point sources.
  • de Gregorio-Monsalvo et al. (2013) de Gregorio-Monsalvo, I., Ménard, F., Dent, W., et al. 2013, A&A, 557, A133
  • de Juan Ovelar et al. (2016) de Juan Ovelar, M., Pinilla, P., Min, M., Dominik, C., & Birnstiel, T. 2016, MNRAS, 459, L85
  • Devillard (2001) Devillard, N. 2001, in Astronomical Society of the Pacific Conference Series, Vol. 238, Astronomical Data Analysis Software and Systems X, ed. F. R. Harnden, Jr., F. A. Primini, & H. E. Payne, 525
  • Dodson-Robinson & Salyk (2011) Dodson-Robinson, S. E. & Salyk, C. 2011, ApJ, 738, 131
  • Dohlen et al. (2008) Dohlen, K., Langlois, M., Saisse, M., et al. 2008, in Proc. SPIE, Vol. 7014, Ground-based and Airborne Instrumentation for Astronomy II, 70143L
  • Dong et al. (2016) Dong, R., Fung, J., & Chiang, E. 2016, ArXiv e-prints [\eprint[arXiv]1602.04814]
  • Fusco et al. (2014) Fusco, T., Sauvage, J. F., Petit, C., et al. 2014, in SPIE Astronomical Telescopes + Instrumentation, ed. E. Marchetti, L. M. Close, & J.-P. Véran (SPIE), 91481U
  • Galli et al. (2013) Galli, P. A. B., Bertout, C., Teixeira, R., & Ducourant, C. 2013, A&A, 558, A77
  • Ginski et al. (2014) Ginski, C., Schmidt, T. O. B., Mugrauer, M., et al. 2014, MNRAS, 444, 2280
  • Ginski et al. (2016) Ginski, C., Stolker, T., Pinilla, P., et al. 2016, ArXiv e-prints [\eprint[arXiv]1609.04027]
  • Girard et al. (2010) Girard, J. H. V., Kasper, M., Quanz, S. P., et al. 2010, in Society of Photo-Optical Instrumentation Engineers (SPIE) Conference Series, Vol. 7736, Society of Photo-Optical Instrumentation Engineers (SPIE) Conference Series
  • Grady et al. (2013) Grady, C. A., Muto, T., Hashimoto, J., et al. 2013, ApJ, 762, 48
  • Grady et al. (2001) Grady, C. A., Polomski, E. F., Henning, T., et al. 2001, AJ, 122, 3396
  • Helou & Walker (1988) Helou, G. & Walker, D. W., eds. 1988, Infrared astronomical satellite (IRAS) catalogs and atlases. Volume 7: The small scale structure catalog, Vol. 7, 1–265
  • Hollenbach & Gorti (2005) Hollenbach, D. & Gorti, U. 2005, in Protostars and Planets V Posters, Vol. 1286, 8433
  • Horne & Baliunas (1986) Horne, J. H. & Baliunas, S. L. 1986, ApJ, 302, 757
  • Hugot et al. (2012) Hugot, E., Ferrari, M., El Hadi, K., et al. 2012, A&A, 538, A139
  • Ishihara et al. (2010) Ishihara, D., Onaka, T., Kataza, H., et al. 2010, A&A, 514, A1
  • Kraus et al. (2008) Kraus, A. L., Ireland, M. J., Martinache, F., & Lloyd, J. P. 2008, ApJ, 679, 762
  • Krautter et al. (1997) Krautter, J., Wichmann, R., Schmitt, J. H. M. M., et al. 1997, A&AS, 123
  • Kuhn et al. (2001) Kuhn, J. R., Potter, D., & Parise, B. 2001, ApJ, 553, L189
  • Lagage et al. (2006) Lagage, P.-O., Doucet, C., Pantin, E., et al. 2006, Science, 314, 621
  • Langlois et al. (2014) Langlois, M., Dohlen, K., Vigan, A., et al. 2014, in Proc. SPIE, Vol. 9147, Ground-based and Airborne Instrumentation for Astronomy V, 91471R
  • Lenzen et al. (2003) Lenzen, R., Hartung, M., Brandner, W., et al. 2003, in Astronomical Telescopes and Instrumentation, ed. M. Iye & A. F. M. Moorwood (SPIE), 944–952
  • Maire et al. (2016) Maire, A.-L., Bonnefoy, M., Ginski, C., et al. 2016, A&A, 587, A56
  • Makarov (2007) Makarov, V. V. 2007, ApJ, 658, 480
  • Manara et al. (2014) Manara, C. F., Testi, L., Natta, A., et al. 2014, A&A, 568, A18
  • Marois et al. (2014) Marois, C., Correia, C., Galicher, R., et al. 2014, in Proc. SPIE, Vol. 9148, Adaptive Optics Systems IV, 91480U
  • Marois et al. (2006) Marois, C., Lafrenière, D., Doyon, R., Macintosh, B., & Nadeau, D. 2006, ApJ, 641, 556
  • Martinez et al. (2009) Martinez, P., Dorrer, C., Aller-Carpentier, E., et al. 2009, The Messenger, 137, 18
  • Mayama et al. (2012) Mayama, S., Hashimoto, J., Muto, T., et al. 2012, ApJ, 760, L26
  • Mazoyer et al. (2016) Mazoyer, J., Boccaletti, A., Choquet, É., et al. 2016, ApJ, 818, 150
  • Merín et al. (2010) Merín, B., Brown, J. M., Oliveira, I., et al. 2010, ApJ, 718, 1200
  • Messina et al. (2010) Messina, S., Desidera, S., Turatto, M., Lanzafame, A. C., & Guinan, E. F. 2010, A&A, 520, A15
  • Milli et al. (2012) Milli, J., Mouillet, D., Lagrange, A.-M., et al. 2012, A&A, 545, A111
  • Min et al. (2012) Min, M., Canovas, H., Mulders, G. D., & Keller, C. U. 2012, Astronomy & Astrophysics, 537, A75
  • Ochsenbein et al. (2000) Ochsenbein, F., Bauer, P., & Marcout, J. 2000, A&AS, 143, 23
  • Pavlov et al. (2008) Pavlov, A., Möller-Nilsson, O., Feldt, M., et al. 2008, in Proc. SPIE, Vol. 7019, Advanced Software and Control for Astronomy II, 701939
  • Pecaut & Mamajek (2013) Pecaut, M. J. & Mamajek, E. E. 2013, ApJS, 208, 9
  • Perrot et al. (2016) Perrot, C., Boccaletti, A., Pantin, E., et al. 2016, A&A, 590, L7
  • Pinilla et al. (2012) Pinilla, P., Benisty, M., & Birnstiel, T. 2012, A&A, 545, A81
  • Pinilla et al. (2015) Pinilla, P., de Boer, J., Benisty, M., et al. 2015, A&A, 584, L4
  • Pojmanski (1997) Pojmanski, G. 1997, Acta Astron., 47, 467
  • Rapson et al. (2015) Rapson, V. A., Kastner, J. H., Millar-Blanchaer, M. A., & Dong, R. 2015, ApJ, 815, L26
  • Roberts et al. (1987) Roberts, D. H., Lehar, J., & Dreher, J. W. 1987, AJ, 93, 968
  • Rousset et al. (2003) Rousset, G., Lacombe, F., Puget, P., et al. 2003, in Astronomical Telescopes and Instrumentation, ed. P. L. Wizinowich & D. Bonaccini (SPIE), 140–149
  • Scargle (1982) Scargle, J. D. 1982, ApJ, 263, 835
  • Schmid et al. (2012) Schmid, H. M., Downing, M., Roelfsema, R., et al. 2012, in Ground-based and Airborne Instrumentation for Astronomy IV. Proceedings of the SPIE, ETH Zürich (Switzerland)
  • Schmid et al. (2006) Schmid, H. M., Joos, F., & Tschan, D. 2006, A&A, 452, 657
  • Siess et al. (2000) Siess, L., Dufour, E., & Forestini, M. 2000, A&A, 358, 593
  • Soummer et al. (2012) Soummer, R., Pueyo, L., & Larkin, J. 2012, ApJ, 755, L28
  • Strom et al. (1989) Strom, K. M., Strom, S. E., Edwards, S., Cabrit, S., & Skrutskie, M. F. 1989, AJ, 97, 1451
  • Thalmann et al. (2008) Thalmann, C., Schmid, H. M., Boccaletti, A., et al. 2008, in Proc. SPIE, Vol. 7014, Ground-based and Airborne Instrumentation for Astronomy II, 70143F
  • Thi et al. (2014) Thi, W.-F., Pinte, C., Pantin, E., et al. 2014, A&A, 561, A50
  • Tuthill et al. (2000) Tuthill, P. G., Monnier, J. D., Danchi, W. C., Wishnow, E. H., & Haniff, C. A. 2000, PASP, 112, 555
  • van der Marel et al. (2013) van der Marel, N., van Dishoeck, E. F., Bruderer, S., et al. 2013, Science, 340, 1199
  • van der Marel et al. (2015) van der Marel, N., van Dishoeck, E. F., Bruderer, S., Pérez, L., & Isella, A. 2015, A&A, 579, A106
  • Vigan et al. (2016) Vigan, A., Bonnefoy, M., Ginski, C., et al. 2016, A&A, 587, A55
  • Vigan et al. (2008) Vigan, A., Langlois, M., Moutou, C., & Dohlen, K. 2008, A&A, 489, 1345
  • Vigan et al. (2010) Vigan, A., Moutou, C., Langlois, M., et al. 2010, MNRAS, 407, 71
  • Vigan et al. (2012) Vigan, A., Patience, J., Marois, C., et al. 2012, A&A, 544, A9
  • Wahhaj et al. (2010) Wahhaj, Z., Cieza, L., Koerner, D. W., et al. 2010, ApJ, 724, 835
  • Weinberger et al. (1999) Weinberger, A. J., Becklin, E. E., Schneider, G., et al. 1999, ApJ, 525, L53
  • Wichmann et al. (1999) Wichmann, R., Covino, E., Alcalá, J. M., et al. 1999, MNRAS, 307, 909
  • Zurlo et al. (2014) Zurlo, A., Vigan, A., Mesa, D., et al. 2014, A&A, 572, A85

Appendix A Stellar Properties

Rotation Period

We have taken advantage of publicly available SuperWASP (Butters et al. 2010) photometry time series data, collected during 2004 - 2006, to determine the rotation period. This consists of 12684 -band measurements with an average photometric precision  mag. After the removal of outliers and low-precision measurements from the time series by applying a moving boxcar filter with 3 threshold, we average consecutive data collected within 1 hr. Finally we are left with 545 averaged magnitudes with an associated standard deviation  mag for the subsequent analysis. To search for the rotation period of RX J1615, we apply the Lomb-Scargle method (LS, Scargle 1982), with the prescription of Horne & Baliunas (1986), as well as the Clean (Roberts et al. 1987) periodogram analyses of the data. In all three SuperWASP observing seasons (2004, 2005 and 2006) we have searched in the period range 0.1 - 100 d, and find with both LS and Clean the same rotation period  d, which is the first determination for the rotation period of RX J1615. Although other peaks due to observation timings and beat frequencies are present in the periodogram, we do not detect any other significant peaks of interest.

RX J1615 was also observed by the All Sky Automated Survey (ASAS, Pojmanski 1997) in the years 2001-2009. Despite the lower photometric precision, our LS and Clean analyses allowed us to measure the same  d rotation period and a peak-to-peak light curve amplitude of  mag.

Figure 12: Spectral Energy Distributions (SED, see legend) and our best fit BT-NextGen model (grey solid line,  K) for RX J1615.

Spectral Type and Temperature

In the literature the spectral type is found in the range from K4 (Merín et al. 2010, based on Spitzer/IRS spectra) to K7 (Manara et al. 2014, using spectroscopy with VLT/X-Shooter), with Wichmann et al. (1999, with visible and NIR photometry) and Krautter et al. (1997, based on visible light spectroscopy) favouring a K5 classification. Using the sequence of intrinsic colors and temperatures of pre-main sequence stars by Pecaut & Mamajek (2013) and comparing it to the measured colors by Makarov (2007) we find that the photometric colors are fully compatible with a K5-K6 star suffering a small or even negligible amount of reddening. The colors of a K4 star would instead indicate a reddening of about E(BV) = 0.15 mag while a K7 star is expected to have redder colors than those observed. Previous estimates of extinction span between  mag (Manara et al. 2014) to  mag (Wahhaj et al. 2010) with an intermediate value  mag by Andrews et al. (2011).

Assuming the spectral type K5-K6 and a low amount of reddening (E(BV)  mag), the corresponding effective temperatures in the scale by Pecaut & Mamajek (2013) are  K for K5 and K6 spectral types respectively. We have also used the available optical, near-IR, and IR photometry to build the observed SED (Figure 12). We have obtained the Far Ultra Violet magnitudes from the GALEX catalogs (Bianchi et al. 2011); the UBVRI magnitudes are taken from Makarov (2007); Denis magnitudes from DENIS Catalogue (Ochsenbein et al. 2000); magnitudes from the 2MASS project (Cutri et al. 2003); W1-W4 magnitudes from the WISE project (Cutri & et al. 2013); IRAS magnitudes from the IRAS Catalog of Point Sources (Helou & Walker 1988); AKARI magnitudes from AKARI/IRC mid-IR all-sky Survey (Ishihara et al. 2010). The SED was fit with a grid of theoretical spectra from the BT-NextGen Model (Allard et al. 2012) and the best fit, shown in Figure 12 is obtained with a model of  K, when adopting a distance  pc and E(BV) = 0.00 mag. We note the presence of a significant UV and IR excess which is common when there is an associated disk. Combining our two estimates we conservatively take the temperature to be  K.

Stellar rotation axis inclination

Using the brightest visual magnitude  mag from the ASAS time series, the distance  pc,  mag, bolometric correction from Pecaut & Mamajek (2013) tables, we derive the stellar luminosity . For  K, we derive the stellar radius . With the stellar radius and the rotation period determined above we can compute the rotational velocity at the stellar surface  km s. From the projected rotational velocity  km s (Wichmann et al. 1999), we derive the inclination of the rotation axis , which is compatible within errors with the inclination of the disk derived by van der Marel et al. (2015, & private communication). Therefore, we infer that the disk is coplanar with the stellar equatorial plane. If we instead assume an extinction  mag (Andrews et al. 2011) to determine and , the inferred inclination .

Age, mass and distance

A comparison of the rotation period of RX J1615 with the distribution of rotation periods of stellar associations of known age helps us to constrain the age. We find that the rotation period  d is longer than for members of the TW Hya association ( d, age = 9 Myr, Messina et al. 2010). Considering that RX J1615 is a pre-main sequence star whose radius is still contracting and its rotation still spinning up, we can state that RX J1615 is not older than about 9 Myr. The measured stellar lithium (Li) equivalent width of  mÅ (Wichmann et al. 1999) fits very well into the Li distribution of the 5 Myr NGC 2264 open cluster (Bouvier et al., submitted), leading us to expect an even younger age compatible with that expected for a CTTS with a gas-rich disk.

Makarov (2007) classified the star as a member of the Lupus association, deriving a kinematic distance of 184 pc. They also noted a group of stars that appear younger ( Myr as opposed to  Myr) and at a greater distance than the bulk of the members of the association. However, the star is not considered in the most recent study of Lupus (Galli et al. 2013).

Comments 0
Request Comment
You are adding the first comment!
How to quickly get a good reply:
  • Give credit where it’s due by listing out the positive aspects of a paper before getting into which changes should be made.
  • Be specific in your critique, and provide supporting evidence with appropriate references to substantiate general statements.
  • Your comment should inspire ideas to flow and help the author improves the paper.

The better we are at sharing our knowledge with each other, the faster we move forward.
""
The feedback must be of minimum 40 characters and the title a minimum of 5 characters
   
Add comment
Cancel
Loading ...
171770
This is a comment super asjknd jkasnjk adsnkj
Upvote
Downvote
""
The feedback must be of minumum 40 characters
The feedback must be of minumum 40 characters
Submit
Cancel

You are asking your first question!
How to quickly get a good answer:
  • Keep your question short and to the point
  • Check for grammar or spelling errors.
  • Phrase it like a question
Test
Test description