Variable dynamics in the inner disk of HD 135344B revealed with multi-epoch scattered light imaging1footnote 11footnote 1Based on observations collected at the European Organisation for Astronomical Research in the Southern Hemisphere under ESO programmes 087.C-0702(A,B), 087.C-0458(B,C), 087.C-0703(B), 088.C-0670(B), 088.D-0185(A), 088.C-0763(D), 089.C-0211(A), 091.C-0570(A), 095.C-0273(A), 097.C-0885(A), 097.C-0702(A), and 297.C-5023(A).

Variable dynamics in the inner disk of HD 135344B revealed with multi-epoch scattered light imaging111Based on observations collected at the European Organisation for Astronomical Research in the Southern Hemisphere under ESO programmes 087.C-0702(A,B), 087.C-0458(B,C), 087.C-0703(B), 088.C-0670(B), 088.D-0185(A), 088.C-0763(D), 089.C-0211(A), 091.C-0570(A), 095.C-0273(A), 097.C-0885(A), 097.C-0702(A), and 297.C-5023(A).

[ Anton Pannekoek Institute for Astronomy, University of Amsterdam, Science Park 904, 1098 XH Amsterdam, The Netherlands Mike Sitko Department of Physics, University of Cincinnati, Cincinnati OH 45221, USA Bernard Lazareff Université Grenoble Alpes, IPAG, F-38000 Grenoble, France CNRS, IPAG, F-38000 Grenoble, France Myriam Benisty Université Grenoble Alpes, IPAG, F-38000 Grenoble, France CNRS, IPAG, F-38000 Grenoble, France Unidad Mixta Internacional Franco-Chilena de Astronoma (CNRS UMI 3386), Departamento de Astronomía, Universidad de Chile, Camino El Observatorio 1515, Las Condes, Santiago, Chile Carsten Dominik Anton Pannekoek Institute for Astronomy, University of Amsterdam, Science Park 904, 1098 XH Amsterdam, The Netherlands Rens Waters Anton Pannekoek Institute for Astronomy, University of Amsterdam, Science Park 904, 1098 XH Amsterdam, The Netherlands SRON Netherlands Institute for Space Research, Sorbonnelaan 2, 3584 CA Utrecht, The Netherlands Michiel Min Anton Pannekoek Institute for Astronomy, University of Amsterdam, Science Park 904, 1098 XH Amsterdam, The Netherlands SRON Netherlands Institute for Space Research, Sorbonnelaan 2, 3584 CA Utrecht, The Netherlands Sebastian Perez Departamento de Astronomía, Universidad de Chile, Casilla 36-D, Santiago, Chile Millennium Nucleus ”Protoplanetary Disks,” Chile Julien Milli ESO, Alonso de Córdova 3107, Vitacura, Casilla 19001, Santiago de Chile, Chile Antonio Garufi Universidad Autónoma de Madrid, Dpto. Física Teórica, Módulo 15, Facultad de Ciencias, Campus de Cantoblanco, E-28049 Madrid, Spain Jozua de Boer Leiden Observatory, Leiden University, P.O. Box 9513, 2300 RA Leiden, The Netherlands Christian Ginski Leiden Observatory, Leiden University, P.O. Box 9513, 2300 RA Leiden, The Netherlands Stefan Kraus University of Exeter, Astrophysics Group, Stocker Road, Exeter, EX4 4QL, UK Jean-Philippe Berger Université Grenoble Alpes, IPAG, F-38000 Grenoble, France CNRS, IPAG, F-38000 Grenoble, France Henning Avenhaus Departamento de Astronomía, Universidad de Chile, Casilla 36-D, Santiago, Chile Millennium Nucleus ”Protoplanetary Disks”, Chile Institute for Astronomy, ETH Zurich, Wolfgang-Pauli-Strasse 27, 8093 Zurich, Switzerland

We present multi-epoch Very Large Telescope/Spectro-Polarimetric High-contrast Exoplanet REsearch observations of the protoplanetary disk around HD 135344B (SAO 206462). The -band scattered light imagery reveal, with high spatial resolution (41 mas, 6.4 au), the disk surface beyond 20 au. Temporal variations are identified in the azimuthal brightness distributions of all epochs, presumably related to the asymmetrically shading dust distribution in the inner disk. These shadows manifest themselves as narrow lanes, cast by localized density enhancements, and broader features which possibly trace the larger scale dynamics of the inner disk. We acquired visible and near-infrared photometry which shows variations up to 10% in the bands, possibly correlated with the presence of the shadows. Analysis of archival Very Large Telescope Interferometer/PIONIER -band visibilities constrain the orientation of the inner disk to and , consistent with an alignment with the outer disk or a minor disk warp of several degrees. The latter scenario could explain the broad, quasi-stationary shadowing in N-NW direction in case the inclination of the outer disk is slightly larger. The correlation between the shadowing and the near-infrared excess is quantified with a grid of radiative transfer models. The variability of the scattered light contrast requires extended variations in the inner disk atmosphere (). Possible mechanisms that may cause asymmetric variations in the optical depth () through the atmosphere of the inner disk include turbulent fluctuations, planetesimal collisions, or a dusty disk wind, possibly enhanced by a minor disk warp. A fine temporal sampling is required to follow day-to-day changes of the shadow patterns which may be a face-on variant of the UX Orionis phenomenon.

protoplanetary disks — radiative transfer — scattering — stars: individual (HD 135344B) — techniques: high angular resolution — techniques: interferometric

Tomas Stolker

0000-0002-5823-3072]Tomas Stolker

1 Introduction

Spatially resolved observations provide detailed insight into the physical and chemical processes occurring in protoplanetary disks. The disk around the intermediate-mass star HD 135344B (SAO 206462) is a suitable target to be observed with high-resolution due to its proximity ( pc; Gaia Collaboration et al., 2016), spatial extent (115 in the scattered light; Grady et al., 2009), low inclination (; van der Marel et al., 2016), and brightness from visible to millimeter wavelengths (e.g., Carmona et al., 2014). The disk is classified as a transition disk (e.g., Espaillat et al., 2014) with a large dust cavity resolved at continuum (sub)millimeter wavelengths ( au at 156 pc; Andrews et al., 2011). In scattered light, two symmetric spiral arms have been detected (Muto et al., 2012), which might indicate the presence of a massive gas giant in the outer disk (Dong et al., 2015; Fung & Dong, 2015; Dong & Fung, 2017). So far, only upper limits on planet masses have been derived from direct imaging observations (3  at 07, assuming hot-start evolutionary models; Maire et al., 2017). The cavity radius in scattered light ( au at 156 pc; Stolker et al., 2016), tracing micron-sized grains in the disk surface, is located inward with respect to the large grains in the midplane which can be explained by planet-induced dust filtration (Garufi et al., 2013). The scattered light cavity coincides with a region in which the surface density of CO gas is significantly reduced (van der Marel et al., 2016).

Pre-main-sequence stars are commonly variable at optical and near-infrared wavelengths on various timescales, for example due to rotational modulation by stellar spots, variable accretion, dust obscuration, and structural changes in the inner disk (e.g., Eiroa et al., 2002). Variability also occurs at mid-infrared wavelengths, for example, Spitzer/IRS spectra show a typical anti-correlation between the amplitude of the near- and mid-infrared emission, indicating changes in the height of the inner disk at sub-au distance and consequent shadowing of the disk further outward (Espaillat et al., 2011). The spectral energy distribution (SED) of HD 135344B contains a large near-infrared excess (; Garufi et al., 2017) due to the presence of hot dust in the innermost disk region (Brown et al., 2007). It was noted by Garufi et al. (2017) that HD 135344B belongs to a sub-category of group I protoplanetary disks with both a large near-infrared excess and spiral arms, as well as shadows in most cases. The near-infrared continuum emission is variable up to 20–30% while the 10 m flux exhibits fluctuations of 60% (Grady et al., 2009; Sitko et al., 2012). Furthermore, Grady et al. (2009) observed an anti-correlation between the strength of the - and /-band fluxes which was linked to geometrical changes of the inner disk.

Multi-wavelength polarimetric differential imaging observations by Stolker et al. (2016) revealed three shadow lanes in the band and a broader shaded region bound by two of the shadow lanes. A comparison with optical images from a month earlier showed that the southern -band shadow lane was not present in the bands, pointing toward a transient or variable origin. Those shadows are presumably cast by dust in the inner disk which is asymmetrically perturbed and/or misaligned with respect to the outer disk (Stolker et al., 2016). Similarly, Wisniewski et al. (2008) found that the scattered light flux from the protoplanetary disk around HD 163296 showed variations between different imagery epochs obtained with the Hubble Space Telescope (HST). This might indicate a time-variable shadowing of the outer disk by scale height variations of the inner disk wall (Wisniewski et al., 2008), in line with the photometric variability in the near-infrared due to structural disk changes near the dust sublimation zone (Sitko et al., 2008). Furthermore, Ellerbroek et al. (2014) reported enhanced extinction in the optical for HD 163296, lasting from a few days up to a year, which was interpreted as caused by a dusty disk wind.

In this paper, we present multi-epoch, polarized scattered light imagery of the protoplanetary disk around HD 135344B that were obtained with the Spectro-Polarimetric High-contrast Exoplanet REsearch (SPHERE; Beuzit et al., 2008) instrument. We aim to detect and characterize brightness variations caused by shading dust in the inner disk. The shadow patterns and their variability allow us probe to the physical processes occurring in the innermost disk region which are not directly accessible by high-contrast imaging instruments. The scattered light images are complemented with multi-epoch visible and near-infrared photometry, that we aim to link to the scattered light variations, and near-infrared interferometry, allowing us to place constraints on the orientation of the inner disk. Furthermore, we will use radiative transfer simulations to quantify the correlation between the scattered light contrast and near-infrared excess in order to estimate the extent of the inner disk variations.

2 Observations and data reduction

2.1 SPHERE/IRDIS dual-polarimetric imaging (DPI)

UT date IntegrationaaMultiplication of the integration time (DIT), the number of integrations (NDIT), the number of polarimetric cycles (NPOL), and the number of half-wave plate orientations (NHWP=4). Airmass SeeingbbSeeing measured by the differential image motion monitor (DIMM) at 0.5 m Wind speedccWind speed at ground level. Coherence timeddCoherence time measured by the multi-aperture scintillation sensor (MASS), except the first epoch which is estimated from the DIMM. Strehl ratioee-band Strehl ratio estimated by SAXO. PSF FWHMffPoint spread function full width at half maximum fitted to the non-coronagraphic flux frames.
(minutes) (arcsec) (m s) (ms) (%) (mas)
2015 May 03ggArchival data from Stolker et al. (2016). 076.8 1.11–1.37 0.69(0.07) 06.1(0.3) 02.4(0.4) 73(4)
2016 May 04 102.4 1.06–1.38 0.52(0.16) 01.4(0.7) 13.1(3.4) 78(2)
2016 May 12 017.1 1.14–1.18 2.14(0.30) 14.1(0.8) 03.0(0.6) 60(7)
2016 June 22 034.1 1.09–1.15 0.63(0.07) 08.1(0.3) 08.3(1.4) 72(2)
2016 June 30 034.1 1.02–1.03 0.37(0.05) 08.0(0.5) 09.5(1.5) 79(1)

Note. – Values in parenthesis provide the standard deviation of the average measured value.

Table 1: SPHERE/IRDIS observations

Imaging polarimetry data sets were obtained on 2016 May 03, 2016 May 11, 2016 June 22, and 2016 June 29 with the near-infrared imager (IRDIS; Dohlen et al., 2008) of SPHERE at the European Southern Observatory’s Very Large Telescope (VLT). Observations were carried out with the broadband  filter (BB_J, 1.245 m) in DPI (Langlois et al., 2014) mode. The pixel scale of the detector is 12.26 mas pix (Maire et al., 2016). An apodized Lyot coronagraph was employed (N_ALC_YJH_S, 185 mas mask diameter), allowing for an integration time of 32 s. The four standard half-wave plate orientations were cycled with 2 or 4 subsequent integrations per half-wave plate orientation. The extreme adaptive optics system (SAXO; Fusco et al., 2006) provided a typical Strehl ratio of 75% in the band (see Table 1).

Seeing conditions were mostly good (07) except on 2016 May 11 when the observations were executed with an average seeing of 21 and the presence of strong winds. Nonetheless, the AO loop remained closed for a total integration time of 17 minutes. The integration time was 34 minutes or longer for the other observations. The point spread function (PSF) of the sequence on 2016 May 03 was partly affected by low-order aberrations caused by low wind speeds (1.4 m s on average). The PSF quality was evaluated from the non-coronagraphic images recorded by the differential tip-tilt sensor after which 13 polarimetric cycles were removed, leaving a total of 11 cycles (47 minutes).

Standard calibration procedures were applied with the SPHERE data reduction and handling (DRH) pipeline (v0.18.0; Pavlov et al., 2008) which included sky subtraction, flat field correction, and bad pixel interpolation. The frames with the horizontally and vertically polarized flux were separated and subsequently processed by a custom pipeline for DPI data. We obtained coronagraphic images with four symmetric satellite spots, induced by a periodic modulation applied to the deformable mirror, before and after the science sequence. These frames were used to determine the position of the star behind the coronagraph mask. To center the coronagraphic DPI data, we interpolated linearly between the start and end position of the star. Stokes  and images were obtained with the double-difference method (Hinkley et al., 2009) and subsequently collapsed with a mean stacking.

To correct for instrumental polarization, we use the method described by Canovas et al. (2011) which assumes that the central star is unpolarized. The signal, which provides an estimate of the noise level in the single-scattering limit, was minimized by stepwise changing the inner and outer radius of the annulus used to measure the (assumed to be unpolarized) signal close to the star. The azimuthal counterparts of the Stokes and images, and , were calculated with an additional minimization applied on the image by correcting for a minor rotational offset of the half-wave plate (Avenhaus et al., 2014) for which the optimized values ranged from –20 to –13. We note that the procedure of minimizing the signal is not strictly valid because part of the scattered light flux from the disk will be present in the image due to multiple scattering (Canovas et al., 2015). The effect will be small when the disk inclination is low, however, the high signal-to-noise ratio (S/N) of the disk detection around HD 135344B might reveal a real signal in the image (Stolker et al., 2016). Finally, the images are rotated by –18 toward the true north orientation (Maire et al., 2016).

Flux frames were obtained at the start of each sequence, by shifting the star away from the coronagraph, with a shorter integration of 0.87 s and an additional neutral density filter (ND1.0) in the optical path to avoid saturation. The flux frames are used to determine the angular resolution of the images (see Table 1), as well as the scattered light contrast of the disk (see Section 3.1). The data reduction procedure included a dark frame subtraction, flat field correction, and bad pixel interpolation. To remove any residual background and bias from the images, we calculated for each detector column the mean pixel value (with the inner 15 masked) and subtracted that value from all pixels in that column. The PSF of HD 135344B was fitted with a 2D Gaussian profile which yielded a typical FWHM of 41 mas. More details on the observations and conditions are provided in Table 1.

2.2 Rapid Eye Mount (REM) visible and near-infrared photometry

HD 135344B was observed with the REM at La Silla, Chile, in June 2016. The visible camera, ROSS2, is a simultaneous multi-channel imaging camera which delivers the bands onto four quadrants of the same CCD detector. A dichroic enables simultaneous observations with the infrared camera, REMIR, with a similar field of view of . Observations were executed with 3 s exposures in the bands and 1 s exposures in the bands. For the photometry, a standard five-position dither pattern was used and additional sky frames were obtained. Data was acquired during a total of twelve nights but the data of two nights were rejected due to thick clouds.

The photometry of HD 135344B was measured differentially with respect to HD 135344A (SAO 206463), an A0V star with a separation of from HD 135344B, which appears to be photometrically stable (Sitko et al., 2012). Differential photometry allowed us to measure with high precision the absolute fluxes of HD 135344B also in variable conditions or with the presence of thin clouds. The magnitudes of HD 135344A were retrieved from the 2MASS catalog (Skrutskie et al., 2006) while the magnitudes were calculated through a transformation of the magnitudes with the relations from Jordi et al. (2006). The magnitudes of HD 135344A were obtained from Sitko et al. (2012):  mag,  mag,  mag, and  mag.

2.3 PIONIER interferometry

We retrieved all the available archival near-infrared interferometric data of HD 135344B from the PIONIER instrument (Le Bouquin et al., 2011) at the Very Large Telescope Interferometer (VLTI). The data were taken during multiple epochs from 2011 to 2013222UT dates: 2011 April 27 and 29, 2011 May 26 and 27, 2011 June 3, 2011 August 7 and 8, 2012 March 6, 28, 29 and 30, 2012 April 27 and 2013 May 15.. The instrument recombines the four auxiliary telescopes which were positioned in the short (A1-B2-C1-D0), intermediate (D0-G1-H0-I1), and long (A1-G1-I1-K0) baseline configurations. The projected baseline, , ranged from 7 to 135 m, enabling a maximum angular resolution of  mas across the seven spectrally dispersed channels in the band (). Each observation of HD 134453B was preceded and followed by an observation of a calibration star to characterize the instrumental and atmospheric contribution to the visibilities and closure phases (i.e., the transfer function). Calibration stars were identified with the SearchCal tool (Bonneau et al., 2006, 2011). The data were reduced with the pndrs package, described in detail by Le Bouquin et al. (2011). The 27 calibrated OIFITS files (Pauls et al., 2005) are available in the Optical interferometry DataBase (OiDB) (Haubois et al., 2016).

The observing conditions were best during the two epochs in 2011 April with a coherence time of 5 ms and a seeing of 05–07, while the conditions were significantly poorer when the other data sets were taken. We rejected low S/N measurements by selecting only the data points where the error estimates were in the first quartile of the total distribution of errors. The dominating factor in the error estimate is the stability of the transfer function during the night which is determined by the temporal scatter of the calibrator visibilities. For each data set, the quality assessment was done independently for the six visibilities, four closure phases, and seven spectral channels. The selected measurements were retrieved from 11 out of 13 epochs, with 2011 May 26 and 2012 April 27 excluded.

3 Results

3.1 Multi-epoch polarized light imagery

Figure 1: Multi-epoch polarized scattered light images of HD 135344B in the band. The columns show from left to right the unscaled images, unscaled images, -scaled images, and an unsharp-masked version of the -scaled images. The field of view of each image is with north and east in upward and leftward direction, respectively. The surface brightness of the images has been normalized to the integrated flux (see main text for details). The dynamical range of the color stretch is fixed in each column except for the unsharp-masked images. The dynamical range of the images is a factor 10 smaller than the partner images. Orange corresponds to positive values, blue to negative values, and black is the zero point. The extent of the coronagraph has been masked out. The major axis position angle of the outer disk, (purple line; van der Marel et al., 2015), and the inner disk, (yellow line; see Section 3.4), are shown in the top row.

Scattered light images are displayed in chronological order (top to bottom) in Figure 1. Besides the newly obtained data sets from 2016, we also show the 2015 -band imagery from Stolker et al. (2016). The first column shows the unscaled images, defined such that positive values correspond to azimuthally polarized flux. The images are normalized to the disk-integrated flux, measured with an annulus aperture between 01 and 10, and shown with identical dynamical range. The second column shows the corresponding images, containing flux with a rotational offset of the direction of polarization with respect to the flux. The third column contains the images with a stellar irradiation correction applied (i.e., -scaling). The fourth column shows unsharp-masked images that were obtained by smoothing the -scaled images with a Gaussian kernel ( mas) and subtracting the smoothed images from the original -scaled images. This procedure enhances the contrast of small scale features by removing low spatial frequencies. The dynamical range of the unsharp-masked images is limited to positive values and for each image separately normalized to the peak intensity.

The SW direction () of the disk, located around the major axis (see Figure 1), appears relatively bright in all -scaled images. In contrast to the shadowing variations, the origin of that brightness enhancement is presumably intrinsic since the bright wedge in scattered light coincides with the (sub)millimeter emission peak of the crescent-shaped dust continuum (Pérez et al., 2014; Stolker et al., 2016). An enhancement of the surface density and/or midplane temperature will elevate the height of the scattering surface, therefore, a larger geometrical cross section of the disk surface is irradiated which increases the scattered light flux. In this work, we focus on local brightness variations between the five epochs. We refer the reader to Muto et al. (2012), Garufi et al. (2013), Stolker et al. (2016) for a detailed analysis and discussion of the scattered light detection of the spiral arms and cavity edge.

The images in Figure 1 contain contributions from noise, residuals of instrumental polarization, and possibly multiple scattered light from the disk. The residual signal could not be removed with the minimization steps explained in Section 2.1 but a more detailed analysis is required to determine if the remaining signal is an instrumental artifact. We distinguish between two different type of signals in the images which appear to be related to the total integration time and therefore the S/N. First, the relative contribution of noise is particularly well visible at separations  mas from the star in the last three epochs for which the total integration time was shortest while the relative contribution is lower in the first two epochs. Second, the images show an enhanced signal within 300 mas of which the relative strength is larger in the first two epochs. The inner signal reaches only mildly above the background noise level in the remaining epochs. The inner signal appears to be variable between epochs revealing a variety of brightness patterns. For example, the first epoch shows a complex pattern of multiple positive and negative lobes whereas the second epoch displays an anti-symmetric signal which is bisected in a positive and negative side. The relative strength of the flux with respect to the flux is quantified in Section 3.3.

Although we deem it likely that the remaining signal is an uncorrected instrumental artifact, we speculate that the increasing strength of the innermost signal with increasing S/N could be a result of multiple scattered light from the cavity edge. A minor fraction of the scattered light will be non-azimuthally polarized because of the small inclination of the outer disk (; van der Marel et al., 2015, 2016). Furthermore, a small fraction of the stellar light will scatter in the extended inner disk atmosphere (see Section 4.3) before it scatters from the outer disk. The signal from those photons will get modulated by the temporal variations in the inner disk such that each epoch may show a different image. However, more detailed analysis and modeling is required to determine if the signal is real and if the temporal variations could be caused by the inner disk.

Figure 2: Polar projections of the -scaled, unsharp-masked images shown in chronological order (top to bottom). North corresponds to and positive position angles are measured east from north. Localized and broad shadow features are indicated with solid and dashed arrows, respectively.

3.2 Asymmetric Q brightness variations

A comparison of the -scaled images in Figure 1 shows epoch-to-epoch brightness variations. Azimuthal brightness minima are visible in all images with variations in their location, shape, and strength. The locality of the brightness minima (e.g., the shadow lane at on 2015 May 03) points toward a shadowing effect, likely caused by dust in the (unresolved) inner disk (Stolker et al., 2016). Furthermore, the brightness variations occur on a timescale similar to the dynamical timescale of the inner disk (the finest temporal resolution is 8 days) and the variability timescale of the near-infrared photometry. Minor brightness variations in the unscaled images are also visible in the region between the cavity edge and the coronagraph, possibly related to the flow of gas and small dust grains from the outer disk.

Azimuthal brightness variations, and their changes between epochs, are more evidently revealed with polar projections of the scattered light images which are displayed in Figure 2. For clarity, we choose the unsharp-masked images for the identification of the shadow features in the polar projections. We caution that applying an unsharp mask may introduce a bias in the identification of brightness variations. However, the shadow features, that we will discuss below and are marked with arrows in Figure 2, are also visible in the regular -scaled images, but the contrast between shadowed and non-shadowed regions is smaller.

Asymmetric illumination/shadowing variations are visible in the scattered light imagery of all five epochs. Here we list the main characteristics of the shadow features, in consonance with the locations that are pointed out in Figure 2:

  1. Epoch 1, 2015 May 03 - Three narrow shadow lanes are present at position angles of , , and , and an azimuthally broader dimming is visible in N-NW direction which is bound by two of the shadow lanes (Stolker et al., 2016).

  2. Epoch 2, 2016 May 04 - The eastern half of the disk is mildly shadowed, approximately in the position angle range of . Deeper shadows are superimposed near the edges of the global shadow. The deepening is particularly well visible in the north, extending at the cavity edge from to . Radially outward, the location of the shadow shows an azimuthal gradient. There is a hint of a localized shadow lane at , approximately colocated with the southern shadow lane in the first epoch.

  3. Epoch 3, 2016 May 12 - The bisection of the brightness distribution from the second epoch seems to have disappeared although the poor observing conditions and the short total integration time (see Table 1) make the identification of the shadows challenging. The broad, northern shadow from the previous epoch is still present. There is also a hint of the broad, southern shadow whereas the narrow southern shadow has disappeared.

  4. Epoch 4, 2016 June 22 - A broad shadow is present between and upon which finer shadow variations are superimposed, including narrow shadow lanes at the boundary with the non-shadowed region, similar to the northern shadow features in the first epoch. The shadow lane at possibly coincides with the location of the NE shadow lane detected in first epoch.

  5. Epoch 5, 2016 June 30 - The cavity edge is shadowed between and while shadowing of the exterior spiral arm only occurs from onwards. The radial extent of the broad shadow increases with increasing position angle, similar to the broad shadows in the first and second epoch. The shadow covers the full radial extent of the disk between and , similar to the shadow feature at the same location in the second and third epoch. The narrow shadow lanes from the fourth epoch seem to have disappeared.

Figure 3: Polarized brightness of the disk normalized to the total Stokes  flux. The plot shows the contrast from four of the epochs (colored dashed lines), the mean contrast (black solid line), and the relative optical depth variation between the minimum and maximum contrast (black dashed line, right -axis). The gray shaded area covers the total variation of the contrast between the epochs. The mean error on the contrast, across all epochs and position angles, is shown on the bottom of the figure (see main text for details). The image from 2016 May 12 has been excluded since it was affected by the poor observing conditions.

The brightness depth of the shadow lanes appears typically larger than the broader shadows (e.g., epoch 1 in Figure 1) indicating larger density enhancements in the inner disk. Indeed, there is no evident radial dependence in the depth of the shadow lanes even though the height of the outer disk increases with radius. The brightness depth of the broad shadows is typically deepest at the cavity edge but weakens toward larger radii (e.g., epoch 4 in Figure 2), presumably an effect of the increasing outer disk height. The broad shadow in N-NW direction (– ) of all epochs shows an azimuthal gradient which is possibly a result of the asymmetric spiral arm perturbation of the disk surface.

The absence of strong azimuthal gradients, caused by a light travel time effect, in the location of the narrow shadows provides a lower limit on the radius from where the shadows are cast. For example, the position angle of a shadow will change by between the inner disk and 80 au (500 mas) if the responsible dust clump is located at 0.15 au (see also Kama et al., 2016). Arguably, some of the shadow lanes show very minor tilts in Figure 2 but the precision and angular resolution of the observations challenge the identification of light travel time effects. The best candidate is the shadow lane at in the 2015 epoch. Stolker et al. (2016) speculate that its azimuthal tilt could be caused by orbital motion in the inner disk from which an orbital radius of 0.06 au (at 140 pc) was estimated. For the other shadow lanes, we may conclude that the dust clumps responsible are presumably located at distances 0.15 au.

A quantification of the azimuthal brightness variations is shown in Figure 3. The mean flux is measured in position angle bins of wide across a radial separation of 01–07 and divided by the angular area of a pixel. The polarized surface brightness (in counts s arcsec) is normalized to the total Stokes  flux (in counts s) which is measured with a circular aperture on the unsaturated, non-coronagraphic flux frames after a correction for the integration time and response of the neutral density filter. The optimal aperture size (15) was determined by measuring the photometric flux with a large range of aperture sizes (up to 30) from which it was established that the total encompassed flux flattened for apertures larger than . The mean error bar is calculated from the standard error on the individual contrast points.

Figure 4: Visible and near-infrared photometry obtained with the Rapid Eye Mount during 10 nights in 2016 June. The mean and standard deviation of the fluxes are provided by the horizontally dashed lines and shaded regions, respectively. Photometric monitoring overlapped with the two most recent SPHERE epochs of which the UT dates are indicated with vertically dotted lines. The arrows point in the direction of an increasing and decreasing near-infrared variation.

The polarized surface brightness contrast in Figure 3 shows typical values in the range of 2–6, except in the SW direction where the contrast goes up to at . The integrated disk brightness consists mainly of signal from the cavity edge and the spiral arms (see unscaled images in Figure 1) which are both asymmetrically shadowed resulting in dimming variations of 20–45%. The contrast variation is minimal (10%) in S-SW direction where the large-scale density and/or scale height enhancement affects the disk surface. The opening angle of the scattering surface is larger in that direction such that shadowing by the inner disk requires dust to be located at higher altitude above the midplane. The contrast variations provide an upper limit on the local changes in optical depth through the inner disk atmosphere. The stellar radiation that is transmitted through the atmosphere will be attenuated by a factor , therefore, the relative change in optical depth, , can be calculated from the minimum and maximum contrast (see Figure 3). However, this only provides an upper limit on the optical depth variations because (i) the total flux does also affect the contrast and is correlated with the shadowing of the outer disk (see Section 4.3), (ii) part of the thermal emission from the inner disk may illuminate the outer disk without being affected by any local obscurations, and (iii) PSF smearing will lower the brightness contrast between shadowed and non-shadowed regions.

3.3 Photometry and scattered light contrast

Multi-epoch photometry in the and bands are displayed in Figure 4, covering ten nights between 2016 June 07 and 2016 July 04 with nearly a daily sampling from June 16 till June 22. The error bars reflect the uncertainty on both the science and calibration star. The early June photometry in the bands starts with a minor increase of , remains approximately constant in mid-June, and increases with by the end of June.

The temporal course of the photometry appears more irregular with variations up to 10% with respect to the mean (horizontally dashed lines in Figure 4) in the bands. During the first half of June 2016, the fluxes increased with approximately 6%, 9%, and 10% in the , , and bands, respectively, following the trend of the fluxes but with a larger fractional increase. In mid-June, also the fluxes remained approximately constant but they increased further from June 22 onwards, in contrast to the fluxes. The final epoch shows a decline in the and bands while the band photometry remained constant. Although the variability in the first part of June seems correlated with the photometry, in the second half no correlation is apparent, that is, the fluxes increased up to 10% with respect to the mean while the photometry remained constant. The total temporal coverage of the photometry is too short to reveal any trends and the sampling is too sparse to resolve possible variations on timescales less than one day.

In addition to the absolute REM photometry, we measured the relative disk photometry of the SPHERE data which is presented in Figure 5. The disk-integrated polarized flux was determined from the and images with an annulus aperture (01–10) centered on the star. The Stokes  flux was measured with a circular aperture (15 radius) from the dedicated flux frames. Absolute pixel values were used for the photometry. The photometric contrast (top panel in Figure 5) is calculated as the ratio of the and Stokes  flux after correcting for the difference in integration time and the response of the neutral density filter. Relative photometry allows for an epoch-to-epoch analysis without requiring an absolute flux calibration, assuming that both the coronagraphic sequence and the total flux data were obtained during similar observing conditions. The and surface brightness errors are computed as the standard deviation within an aperture centered on each pixel with a radius of 62 mas (i.e., 1.5 resolution elements) and propagated accordingly to an integrated error (3%–5%). The uncertainty on the total flux (1%–2%) is computed in a background-limited region as the standard error of the sum, , with an annulus aperture equal in size to the aperture used for the Stokes  photometry.

Figure 5: Top: integrated polarized scattered light contrast (black crosses) of the five -band SPHERE epochs. Bottom: the integrated Stokes  flux (purple squares) and the disk-integrated flux (red circles), shown in arbitrary normalized units. The disk-integrated flux (green circles) is computed from the absolute pixel values and shown relative to the flux with a factor 10 enhancement. The uncertainties are given at a level.

The integrated contrast in Figure 5 varies between and . The photometry is not calibrated so only relative variations of and are meaningful. The second epoch shows a consistent decrease of both the and Stokes  flux with respect to the first epoch. In the third epoch, the contrast decreased by 20% possibly due to the poor observing conditions (see Section 2.1), particularly affecting the photometry. During the last two epochs, the relative increase of the total flux is large compared to the flux, resulting in relatively low contrast. The increase of the total flux during the last epochs seems consistent with the REM photometry (see June 22 and 30 in Figure 4). The relative photometry shows an increase in the second epoch due to the residual within 200 mas (see Section 3.1), possibly caused by multiple scattered light from the inner and/or outer disk, while the relative photometry in the third epoch is larger due to the enhanced noise residual.

3.4 Parametric model fitting of the visibilities

The normalized, squared visibilities, , across the spectrally dispersed -band channels of the multi-epoch PIONIER observations are displayed in the left panel of Figure 6. The visibilities decrease continuously with increasing spatial frequency (i.e., ), indicating that the region from which the -band flux originates is resolved. At the longest baselines, there appears no turnover point to an asymptotic value of the stellar flux so the circumstellar emission is not over-resolved. The closure phases, , are consistent with zero within the error bars (), therefore, the brightness distribution of the inner disk is point symmetric within the uncertainties and at the spatial resolution probed by the observations. There is no significant dispersion visible, for visibility points obtained with baselines of similar lengths but different position angles, which implies that the inclination of the inner disk is small. The diagonal scatter of the visibilities is an effect of chromatic dispersion (Lazareff et al., 2017). Coverage of the -plane is shown in the right panel of Figure 6.

The orientation and characteristic radius of the inner disk -band emission is determined by fitting, in Fourier space, a parametric model to the visibilities, following the procedure described in detail by Lazareff et al. (2017). This allowed us to apply a minimization and to assign formal error bars to the inferred parameter values. We parameterize the -band emission with an elliptical brightness distribution that is radially parameterized by a weighted combination of a Gaussian and pseudo-Lorentzian profile. The inner rim is not fully resolved by the longest baselines which justifies an ellipsoidal distribution instead of a broadened ring.

Figure 6: Left: squared visibilities () of in the VLTI/PIONIER -band channels with (top) and the fitting residuals of the ellipsoidal brightness model (bottom). Right: coverage of the -plane, shown with the same color coding as the visibilities.

The best-fit model () corresponds to an inclination and major axis position angle of and , respectively. The values are, within the uncertainties, very similar to those of the outer disk (see discussion in Section 4.1.1). The best-fit position angle is shown in Figure 1 in comparison with the outer disk value from van der Marel et al. (2015). The half-flux semi-major axis is  mas ( au) which implies that a significant fraction of the -band emission originates from within the silicate sublimation radius ( au; Carmona et al., 2014). A detailed overview of the fitting results is provided in Appendix A where the best-fit values of all parameters are listed with their dependence on the cutoff level of the selection criterion for the points.

We caution that the visibilities were combined from multiple epochs while the inner disk is variable on a timescale of days or less (see Section 3.3). The visibilities in Figure 6 depend on the absolute -band flux and the relative contributions of the star and disk. Therefore, an additional uncertainty has been introduced by combining the visibilities from multiple epochs while the absolute -band flux is variable and not measured at the nights of the observations. Also, we made the assumption that the orientation of the inner disk did not change due to precession between those epochs.

4 Discussion

In Section 4.1 we will discuss the available constraints on the orientation of the inner and outer disk, as well as the presence of the broad, quasi-stationary shadow. In Section 4.2, we will provide an observational perspective on the variability and we will discuss several processes that may affect the inner disk dynamics. In Section 4.3, we will estimate the extent of the inner disk variations by quantifying the correlation between the thermal emission and the scattered light flux with a grid of radiative transfer models.

4.1 Constraints on the inner disk (mis)alignment

4.1.1 Near-infrared and submillimeter interferometry

The relative orientation of the inner and outer disk is determined by their inclination with respect to the plane of the sky and their position angle, with the misalignment defined by the angle between the normal vectors of the two midplane orientations. The outer disk’s orientation of HD 135344B has been determined in several studies (see Carmona et al., 2014, for an overview). Here we list those values that were derived from spatially resolved (sub)millimeter CO observations: and (Lyo et al., 2011), and (van der Marel et al., 2015), and (van der Marel et al., 2016). In Section 3.4, we found that the inner disk inclination and position angle are and , respectively. The values are consistent with the orientation of the outer disk within the uncertainties of the model fitting and the available values for the outer disk. This result is in contrast to several earlier studies which suggested a significant misalignment between the inner and outer disk (Fedele et al., 2008; Grady et al., 2009; Stolker et al., 2016). However, we can not exclude a minor disk warp given the uncertainties on both the inner and outer disk orientation.

Figure 7: Raytraced images of the DIANA radiative transfer model of HD 135344B (Woitke et al., 2016). Left: the inner disk is aligned with the outer disk. Right: the misalignment of the inner disk is , by adopting the best-fit inclination and position angle from Section 3.4. Images have been convolved with a Gaussian kernel ( mas) to match the angular resolution of the SPHERE imagery. Surface brightness values are provided along the major and minor axis direction of the outer disk.

A possible misalignment of the two disk components relies also on the identification of the near and far side of both the inner and outer disk. Spatially resolved observations of CO gas show a redshifted and blueshifted velocity in SW and NE direction, respectively (e.g., Pérez et al., 2014). This implies that the near side of the outer disk is along SE direction of the minor axis if we assume that the spiral arms in scattered light follow a trailing motion. For the inner disk, there is no direct constraint on the near and far side at the angular resolution of the PIONIER observations. However, the misalignment will be if the near side of the inner disk is in NW direction, a scenario that can be excluded from the absence of two stationary shadow lanes similar to HD 142527 (; Marino et al., 2015) and HD 100453 (; Benisty et al., 2017). This means that the near side of the inner disk is located in SE direction and a possible misalignment will be minor.

4.1.2 Quasi-stationary shadowing by a minor disk warp?

The location and width of the shadows provide further constraints on the orientation of the inner disk with respect to the outer disk. Three scenarios remain possible: (i) the inner disk and outer disk are aligned such that shadowing occurs only through uplifting of dust in the inner disk atmosphere, (ii) a minor misalignment () is present which might not cast a stationary shadow but additional variations in the inner disk will cause preferential shadowing in the direction where elevation of the inner disk above the outer disk midplane is largest, (iii) an intermediate misalignment () is present which casts a broad, stationary shadow when the misalignment of the inner disk is similar to the opening angle of the scattering surface of the outer disk. An example of the third scenario is the broad shadow detected with HST on the TW Hya disk (; Rosenfeld et al., 2012).

The inclination of the inner disk could be either smaller or larger than the outer disk given the uncertainties on its orientation (see Section 2.3) which could result in a broad shadow in approximately NW or SE direction, respectively. Interestingly, a broad shadow seems present in all images in N-NW direction (see Figure 2) although with small variations in its precise location, shape, and depth. This may imply that the inclination of the inner disk is slightly smaller (i.e., more face-on) than the outer disk. Variable fine structure is present in the N-NW shadow which requires additional optical depth variations through the atmosphere of the inner disk as will be discussed in more detail in Section 4.2.

A broad dimming was also present in NW direction of the Subaru/HiCIAO -band imagery by Muto et al. (2012) and the VLT/NACO and -bands imagery by Garufi et al. (2013). This was interpreted in both studies as a depolarization effect which can occur when the disk is inclined because the polarization efficiency peaks around the major axis where scattering angles are close to (e.g., Murakawa, 2010; Min et al., 2012). The inclination of the outer disk is relatively small so a strong depolarization effect might not be expected. We speculate that the broad dimming in the HiCIAO and NACO imagery might be the same quasi-stationary shadowing effect that is seen in the SPHERE imagery, possibly related to a minor disk warp.

To illustrate this scenario, we have adopted the DIANA (Disc Analysis; Woitke et al., 2016) radiative transfer model of HD 135344B. The left image in Figure 7 displays the raytraced image for a setup in which the inner disk is aligned with the outer disk (, ; van der Marel et al., 2015). A mild depolarization effect is seen in both directions of the minor axis but the effect is slightly stronger on the near side (SE) due to the flaring geometry of the disk surface (Min et al., 2012). In the right image of Figure 7, we changed the orientation of the inner disk to the best-fit values from Section 3.4. The misalignment here is , that is, comparable to the disk warp seen in the debris disk around  Pic (; Heap et al., 2000). The atmosphere of the inner disk casts a mild shadow on the outer disk in NW direction, reminiscent of the broad shadow on the cavity edge of the SPHERE images in Figures 1 and 2. In opposite direction, the outer disk becomes more strongly irradiated which introduces an azimuthal brightness modulation along the cavity edge (see also Rosenfeld et al., 2012), although intertwined with the modulation by the polarization efficiency of the outer disk.

4.2 Origin of the shadows and their variability

4.2.1 Observational perspective and variability timescale

The most prominent azimuthal brightness variations were identified in Section 3.2 and we noticed that the shadows can be classified into two categories, (i) localized shadow lanes and (ii) broader shadows that are tens of degrees wide. The shadow variations are presumably caused by small density enhancements in the atmosphere of the inner disk with the optical depth variations being approximately smaller than unity (see Figure 3). Although a broad shadow is present in N-NW direction of all epochs, none of the shadows are fully stationary. This indicates that the atmosphere of the inner disk is a dynamical environment in which the gas and the micron-sized dust grains (which are dynamically coupled) are subjected to processes that change their vertical distribution on a fast timescale.

Variability of the shadows might be caused by orbital and/or vertical motion of dust enhancements in the inner disk. The vertical response of the disk occurs approximately on the Keplerian timescale meaning that a local perturbation of the disk will settle to an equilibrium state within approximately one orbit. However, the orbital timescale provides only a lower limit on the response timescale as it depends on the heating and cooling timescales of the gas. The -band flux is emitted from a characteristic radius of 0.11 au (see Section 3.4) which corresponds to a Keplerian timescale of 10 days, similar to the finest temporal sampling of the SPHERE imagery (8 days). However, the shadows are presumably cast further outward as inferred from the absence of significant light travel time effects (see Section 3.2). The temporal sampling of the SPHERE observations is too sparse to determine the timescale by which the shadow features appear and disappear, therefore, disentangling variability due to orbital motion and the (dis)appearance of shadows is not possible.

The fluxes in Figure 4 show variations up to 10% on a timescale of days to weeks (see also Grady et al., 2009; Sitko et al., 2012) while the fluxes varied only by 1–2%, indicating that mainly the thermal emission from the inner disk is affecting the near-infrared variability. For reference, the -band flux of HD 135344B consists of 69% stellar radiation and 31% inner disk emission (see Section 4.3). The fast variability timescale of both the shadows and the near-infrared photometry may point toward a common origin in the inner disk.

Small variations of the visible photometry might be related to episodic accretion events, photospheric or chromospheric activity, stellar pulsations, minor attenuation variations by an (optically thin) dust envelope, or increased scattering from the inner disk. The scattered light flux from the inner disk is in the optical 1% of the total flux so structural changes in the inner disk may cause a change both in thermal emission and scattered light. In that case, the photometric variations in the visible should be correlated with variations in the near-infrared which seems the case in the first half of the REM photometry but not the second half.

4.2.2 Inner disk processes affecting the dust dynamics

Several processes may have an effect on the dynamics and distribution of the gas and dust in the inner disk, possibly causing shadow variations on the outer disk. For example, hydrodynamical fluctuations, such as turbulent eddies and filaments, may produce short-lived obscuration events (Dullemond et al., 2003; Flock et al., 2017). Or, catastrophic collisions between planetesimals, possibly stirred-up by a planet (Kenyon & Bromley, 2004), will locally enhance the dust density, although a gaseous environment will affect the dust dynamics differently than in a debris disk. Dedicated simulations are required to determine if such processes could produce disk perturbations that are localized and strong enough to explain the narrow shadows. Precession of a disk warp, for example driven by a companion (Lai, 2014), will result in a variable location of the casted shadow (Debes et al., 2017). The broad shadow in N-NW direction of HD 135344B appears approximately stationary so a fast precession can be excluded if the shadow is cast by a minor disk warp.

Variations in the inner disk might also be related to star-disk interactions. HD 135344B is an F4V-type star with a weak magnetic field ( G; Hubrig et al., 2009) so a magnetic coupling to and warping of the inner disk seems unlikely. Indeed, with the absence of a breaking mechanism, the star has been able to spin-up to a near break-up rotational velocity (Müller et al., 2011). The fast rotation might drive a viscous decretion disk, similar to classical Be stars (Rivinius et al., 2013), by which gas from the stellar atmosphere spreads outward, possibly creating disturbances in the inner disk by an interaction with the inward accretion flow. Also accretion may play a role in the distribution of material in the inner disk. Fairlamb et al. (2015) measured a rate of  yr and Sitko et al. (2012) determined a factor of two variation during the course of a few months. Therefore, accretion of gas and dust from the outer disk might be an irregular process, possibly mediated by one or multiple companions inside the large dust cavity, such that the inner disk gets asymmetrically replenished.

The near-infrared emission, hydrogen line fluxes, and the He i line profile are all variable on various timescales, related to the processes occurring in the inner disk and near the stellar surface (Grady et al., 2009; Sitko et al., 2012). The variable P Cygni profile of the He i line is likely related to a wind whose orientation changes on timescales of a day or less (Sitko et al., 2012) which is launched in the star-disk interaction region (e.g., Edwards, 2009). Micron-sized dust grains in the inner disk atmosphere could be entrained by a photoevaporative wind that is driven by the UV radiation of the star (Owen et al., 2011; Hutchison et al., 2016) or dust could be uplifted by a centrifugally driven disk wind (Bans & Königl, 2012). An extended low-density atmosphere could be supported by the magnetic field of the inner disk which may explain a large near-infrared excess and possible shadowing (Turner et al., 2014). Alternatively, the central star could drive a wind from the circumplanetary disk of a planetary companion (Tambovtseva et al., 2006) or disk perturbations by a companion on an inclined orbit may also cause an asymmetric illumination of the disk (Demidova et al., 2013). Three-dimensional radiation nonideal magnetohydrodynamical simulations show turbulent velocities in the inner disk up to 10% of the sound speed and a nonaxisymmetric shadow on the outer disk cast by a dead zone-induced vortex (Flock et al., 2017). An inner disk vortex will orbit the star with a Keplerian velocity so it is likely not responsible for the broad, quasi-stationary shadow.

4.2.3 Face-on variant of the UX Orionis phenomenon?

Although the origin of the shadows remains uncertain, we notice a resemblance between the photometric variations of UX Orionis stars (UXORs) and the spatially resolved shadow variations on the disk surface around HD 135344B. UXORs are a subclass of Herbig Ae/Be stars which are characterized by sudden declines in brightness up to several magnitudes in the optical, associated with increased extinction and polarization, which suggests changes of the column density in the line of sight toward the star (Waters & Waelkens, 1998). It has been proposed that such photometric variations could be caused by orbiting dust clouds when the disk is observed almost edge-on (Grinin et al., 1994) although alternative explanations involving disk winds (Grinin & Tambovtseva, 2003) and turbulent filaments exist (Dullemond et al., 2003). Similar processes could be invoked to explain the variability of UXORs and HD 135344B, which may suggest that HD 135344B is a face-on variant of the UXOR phenomenon.

4.3 Radiative transfer models of a shadowed outer disk

Figure 8: Top: radiative transfer simulations of the intensity of the outer disk (colored map) and the -band photometry (white contours). Center: scattered light contrast (colored map) with the mean and 1 of the REM -band flux (black dashed lines) and the contrast of the SPHERE imagery (white dashed lines). The white cross denotes to the DIANA model of HD 135344B (Woitke et al., 2016). Bottom: spectral energy distributions (SEDs) of all models (black solid lines), superimposed by SED of the DIANA model (yellow dashed line) and the photometry (red points; Carmona et al., 2014). Error bars of the photometry have been excluded when they are smaller than the symbol.

Shadow variations on the outer disk depend on changes in the vertical distribution of dust in the inner disk. More specifically, the strength of the shadowing is set by the inner disk radius where the transition region from optically thin to thick occurs highest above the midplane. The height of this transition region is mainly determined by the pressure scale height and the surface density if the dust opacities are constant throughout the inner disk. For a flaring geometry, the thickest part will be close to the outer edge of the inner disk. Alternatively, a puffed-up inner rim may shadow the outer disk in case the exterior of the inner disk is fully shadowed by the rim (Dullemond & Dominik, 2004; Dong, 2015).

4.3.1 Parametric model setup

To understand quantitatively the effect of the inner disk on the near-infrared excess and the scattered light contrast, without making an assumption about the origin of the shadows, we have constructed a grid of radiative transfer models which are meant as a proof of concept rather than an accurate fit to the data. We adopted the DIANA model setup of HD 135344B which provides a multiwavelength fit to the SED and several other gas and dust observables (Woitke et al., 2016). The inner disk is aligned with the outer disk and ranges from 0.16 au to 0.21 au (at 140 pc), slightly beyond the characteristic radius inferred from the -band visibilities (see Section 3.4), with a surface density profile parameterized as and a sharp inner rim. The pressure scale height profile is parameterized as with the normalization provided by a reference aspect ratio, , at  au. The negative flaring index implies a decrease of the scale height with increasing radius which is expected due to direct heating of the inner rim by the star. The radial extent of the inner disk is small so the flaring index has only a minor impact on the disk geometry. The grid covers values of the inner disk dust mass in the range of and the inner disk aspect ratio, , in the range of 0.005–0.25. The flaring index of the outer disk is 1.14 with a reference aspect ratio of 0.16 at 50 au. The radiative transfer and image raytracing were done with MCMax3D, a Monte Carlo continuum radiative transfer code (Min et al., 2009).

4.3.2 Polarized scattered light versus thermal emission

For each model, we computed the disk-integrated intensity of the outer disk and the total -band flux which are displayed in the top panel of Figure 8. Several effects of the inner disk mass and the aspect ratio on the polarized intensity are evident in the colored map. Increasing the aspect ratio from very small values up to 0.05 results in a larger fraction of the stellar light being reprocessed by the inner disk and reemitted in the band, thereby increasing the scattered light flux from the outer disk. Similarly, the increase of the polarized intensity with increasing dust mass is also the result of a larger fraction of thermal radiation from the inner disk scattering from the outer disk. In this regime of the aspect ratio, the opening angle of the inner disk is too small to shadow the outer disk.

Shadowing of the outer disk by the inner disk starts to have an effect for such that the opening angle of the inner disk atmosphere is comparable to or larger than the opening angle of the scattering surface of the outer disk. The polarized intensity decreases with increasing aspect ratio when the inner disk dust mass is smaller than because the optical depth through the inner disk atmosphere toward the outer disk increases. Consequently, a larger fraction of the stellar light is attenuated by extinction in the inner disk. For inner disk masses larger than , the effect of enhanced reprocessing by the inner disk starts to dominate over the shadowing effect. More specifically, even though increasing the aspect ratio beyond 0.05 will increase the amount of shadowing of the outer disk, enhanced reprocessing by a larger amount of dust in the inner disk results in an increased irradiation of the outer disk in the band. Since the exponent is negative for both the surface density and the scale height profile, the inner disk is thickest at the location of the inner rim (i.e., the radius where the vertical surface is located highest above the midplane). As a result, reprocessed stellar light by the inner rim can directly irradiate the outer disk, in contrast to an inner rim which is partly shielded from the outer disk by a flaring scale height profile.

The -band photometry of the radiative transfer models is superimposed on the polarized surface brightness in the top panel of Figure 1. The inner disk dust mass and aspect ratio have a similar effect on the -band flux because both parameters affect the amount of reprocessing by the inner disk. While the stellar -band flux remains constant, the thermal emission from the inner disk is positively correlated with both the inner disk dust mass and the reference aspect ratio. The SEDs of all models are shown in the bottom panel of Figure 1. As expected, the near- and mid-infrared excess varies greatly in the explored parameter space.

4.3.3 Observed and simulated scattered light contrast

An absolute flux calibration of the SPHERE data is not possible with high-precision so we use a relative brightness measurement to compare the results from the grid of models with the observations. The central panel of Figure 8 displays the polarized surface brightness contrast that has been computed with the same procedure described in Section 3.2, here simply given by the ratio of the mean polarized intensity of the outer disk and the total -band flux. The -band flux and shadowing of the outer disk are interlinked in the lower part of the disk mass regime, therefore, the contrast distribution across the grid appears similar to the individual components in the top panel. In the regime of the large inner disk mass and aspect ratio, the contrast balances between the increasing scattered light flux and the increasing -band photometry. The latter dominates over the former such that the contrast decreases toward the top right of the distribution, yet, the polarized intensity has a significant effect on the contrast. This means that using a relative flux measurement instead of an absolute one may lead to a degeneracy in the interpretation because a low contrast could correspond to either a high or low polarized intensity.

The azimuthally and temporally averaged scattered light contrast of the SPHERE observations, as well as the REM -band photometry, are superimposed on the contrast from the radiative transfer models in the central panel of Figure 8. We rejected the third quadrant because the brightness of the SW part of the disk is mainly determined by the enhanced midplane density which would otherwise bias the characteristic contrast variations by the shadowing. Significant changes in the total dust mass of the inner disk are not likely to occur as the local viscous timescale will be multiple orders of magnitude longer than the timescale probed with the SPHERE observations. For a fixed dust mass of from the DIANA model, the observed scattered light contrast corresponds to aspect ratio values in the range of 0.11–0.22, likely beyond what is expected from a hydrostatically supported disk. The observed -band fluxes on the other hand cover a lower aspect ratio regime (0.07–0.1). A less extended atmosphere will be required when a small misalignment is present between the inner and outer disk (see Section 4.1) in which case the estimated value of will be smaller.

4.3.4 Evidence for extended variations in the inner disk

We infer from Figure 8 that extended variations in the inner disk atmosphere () can explain the variations of the scattered light contrast and the related shadowing. However, the predicted -band flux from an extended atmosphere is too large compared to the REM photometry. The discrepancy between the scattered light contrast and -band flux possibly points to the uncertain origin of the near-infrared excess. We want to stress that the goal of the radiative transfer models is to provide a quantitative estimate of the required extent of the inner disk atmosphere without making an assumption about the origin of the near-infrared excess (e.g., super-refractory grains, disk wind, magnetically supported atmosphere). Nonetheless, alternative disk structures or dust properties could be considered to dissolve the discrepancy. For example, an inner disk consisting of an optically thick but geometrically thin component and an optically thin but geometrically extended component. The relative strength of the near-infrared excess and the outer disk shadowing is also sensitive to the abundances of amorphous silicates, amorphous carbon, and metallic iron. Amorphous laboratory silicates have a single scattering albedo close to unity (Dorschner et al., 1995) while the absorption cross section of carbon and iron grains is larger (Zubko et al., 1996). Therefore, increasing the relative amount of silicate grains will lower the near-infrared excess while the extinction through the atmosphere may remain large enough to cast a shadow on the outer disk.

We also tested the dependence of the turbulence mixing strength, set by the dimensionless viscosity parameter , on the scattered light contrast and near-infrared excess. The parameter controls the dust settling in the radiative transfer models by assuming an equilibrium between upward turbulent mixing and downward gravitational settling (Woitke et al., 2016). The turbulence/settling parameter mainly affects the brightness contrast when its value is such that even the micron-sized grains slightly settle. The maximum value affecting the brightness contrast increases with increasing aspect ratio. The dependence on the flaring index of the inner disk was also tested. As expected, for a given reference aspect ratio there appeared no significant dependence on the flaring index because of the narrow radial extent of the inner disk.

5 Summary and conclusions

We have presented a multi-epoch scattered light study of the protoplanetary disk surrounding HD 135344B, an isolated pre-main-sequence F4V-type star in the Scorpius OB2-3 association. Polarimetric differential imaging observations with VLT/SPHERE in the band revealed, with a spatial resolution of 6.4 au, azimuthal shadowing variations on the outer disk, related to the vertical dust distribution of dust in the inner disk. The imagery shows irregularly changing shadow patterns between all epochs although similarities have been identified. Shadows appear both as localized lanes and broader structures, typically colocated, which likely trace small scale perturbations and large scale dynamics, respectively. A broad, quasi-stationary shadow is present in N-NW direction of all scattered light images, in particular well visible at the cavity edge where the outer disk scattering surface is lowest. This might suggest that the inner disk is misaligned by several degrees, requiring additional optical depth enhancements through the inner disk atmosphere to shadow the disk further outward. However, fitting of a parametric brightness model to VLTI/PIONIER -band visibilities provided a best-fit inclination and position angle which is, within the uncertainties, consistent with the outer disk. The photometry showed only minor variations (1–2%) in the bands while the fluxes varied up to 10%, indicating significant changes in the amount of reprocessing of stellar light by the inner disk. Variability of the photometry and shadowing appear to be correlated, both are related to structural changes in the inner disk which we have quantified with a grid of radiative transfer models. The observed variations in scattered light contrast require extended variations in the inner disk atmosphere (), beyond what is inferred from the near-infrared excess alone, highlighting the uncertainty about the origin of the near-infrared excess. The variability of the shadows is likely related to the structure and dynamics of the inner disk, whose -band emission originates from a characteristic radius of 0.11 au, that is, inside the silicate sublimation zone. Asymmetric shadowing variations by the inner disk might be caused by mechanisms such as turbulent fluctuations, planetesimal collisions, a dusty disk wind, or asymmetric and episodic accretion. Simultaneous observations of scattered light, photometry, spectroscopy, and/or near-infrared interferometry will provide more stringent constraints on the driving processes in the inner disk for which an approximate daily sampling is required to trace the fast disk dynamics.

We would like to thank Doug Lin, Mario Flock, and the other attendees of the Ringberg workshop on ‘The Atmospheres of Disks and Planets’ for insightful discussions. We also thank the referee for providing constructive comments and suggestions. This research has made use of the OiDB333Available at and SearchCal444Available at services at the Jean-Marie Mariotti Center. MB acknowledges funding from ANR of France under contract number ANR-16-CE31-0013 (Planet Forming Disks). SP acknowledges support from CONICYT-Gemini grant 32130007. SK acknowledges support from an STFC Rutherford fellowship (ST/J004030/1) and an ERC Starting Grant (Grant Agreement No. 639889). \facilitiesVLT:Melipal (SPHERE), VLTI (PIONIER) \softwareMCMax3D (Min et al., 2009), SPHERE DRH (Pavlov et al., 2008), SearchCal (Bonneau et al., 2006, 2011), Astropy (Astropy Collaboration et al., 2013)

Appendix A Details on the visibility fitting

The PIONIER -band visibilities can be expressed as (Lazareff et al., 2017),


where is the visibility of the circumstellar component, and are the spectral index of the stellar and circumstellar component, respectively, and are the fractional flux of the stellar and circumstellar component at a reference wavelength , and is the wavelength of the spectral channel. The flux fractions of the star and circumstellar disk are normalized such that .

As described in Section 3.4, the visibilities are fitted with an ellipsoidal brightness distribution which is parameterized by a weighted combination of a Gaussian, , and a pseudo-Lorentzian, , radial distribution which allows for some freedom in the steepness of the asymptotic decay. The brightness distributions are derived from,


followed by an anamorphosis along the minor axis, and where is the emission radius and the semi-major axis of the half-light isophote. The Hankel transforms of the brightness distributions take a simple analytical form. The visibility of the circumstellar component is,


where is a weighting factor for the contribution of the pseudo-Lorentzian component, and and are the Hankel transforms of the Gaussian and pseudo-Lorentzian brightness distribution, respectively.

The best fit results in Section 3.4 were obtained by selecting the first quartile of the visibility data, on the basis of individual error estimates. Table 2 provides an overview of all the fitting results of the PIONIER visibilities to show the effect of the cutoff level on the best-fit values. The cutoff level in the cumulative distribution of the error estimates is stepwise loosened from 25% up to 100%, that is, the first quartile and all available data points, respectively. The error estimates for individual data points are derived by the PIONIER reduction pipeline from the internal dispersion on a short timescale. However, experience shows that other errors, not captured by the estimate, become increasingly prevalent as the observing conditions degrade. Our choice of the first quartile is a compromise between a severe selection that discards valid information and a loose selection that includes corrupted data, as is reflected in the increase of the reduced in Table 2.

Cutoff PA
(%) (deg) (deg) (mas)
25 199 1.10
35 278 1.44
45 358 1.67
55 437 1.85
65 517 2.07
75 597 2.28
85 676 2.76 0
95 756 2.99
100 796 3.07

Note. – Table columns (from left to right): cutoff level in the cumulative distribution of the error estimates, number of points, reduced , inclination, position angle of the major axis, half-flux semi-major axis, spectral index of the circumstellar component, fractional flux of the circumstellar component, and weighting factor of the pseudo-Lorentzian profile. The uncertainties are provided at a 1 level. Most parameter values appear stable over a broad range of quality cutoff levels.

Table 2: Visibility fitting results


  • Andrews et al. (2011) Andrews, S. M., Wilner, D. J., Espaillat, C., et al. 2011, ApJ, 732, 42
  • Astropy Collaboration et al. (2013) Astropy Collaboration, Robitaille, T. P., Tollerud, E. J., et al. 2013, A&A, 558, A33
  • Avenhaus et al. (2014) Avenhaus, H., Quanz, S. P., Schmid, H. M., et al. 2014, ApJ, 781, 87
  • Bans & Königl (2012) Bans, A., & Königl, A. 2012, ApJ, 758, 100
  • Benisty et al. (2017) Benisty, M., Stolker, T., Pohl, A., et al. 2017, A&A, 597, A42
  • Beuzit et al. (2008) Beuzit, J.-L., Feldt, M., Dohlen, K., et al. 2008, Proc. SPIE, 7014, 701418
  • Bonneau et al. (2006) Bonneau, D., Clausse, J.-M., Delfosse, X., et al. 2006, A&A, 456, 789
  • Bonneau et al. (2011) Bonneau, D., Delfosse, X., Mourard, D., et al. 2011, A&A, 535, A53
  • Le Bouquin et al. (2011) Le Bouquin, J.-B., Berger, J.-P., Lazareff, B., et al. 2011, A&A, 535, A67
  • Brown et al. (2007) Brown, J. M., Blake, G. A., Dullemond, C. P., et al. 2007, ApJ, 664, L107
  • Canovas et al. (2011) Canovas, H., Rodenhuis, M., Jeffers, S. V., Min, M., & Keller, C. U. 2011, A&A, 531, A102
  • Canovas et al. (2015) Canovas, H., Ménard, F., de Boer, J., et al. 2015, A&A, 582, L7
  • Carmona et al. (2014) Carmona, A., Pinte, C., Thi, W. F., et al. 2014, A&A, 567, A51
  • Debes et al. (2017) Debes, J. H., Poteet, C. A., Jang-Condell, H., et al. 2017, ApJ, 835, 205
  • Demidova et al. (2013) Demidova, T. V., Grinin, V. P., & Sotnikova, N. Y. 2013, Astronomy Letters, 39, 26
  • Dohlen et al. (2008) Dohlen, K., Langlois, M., Saisse, M., et al. 2008, Proc. SPIE, 7014, 70143L
  • Dong et al. (2015) Dong, R., Zhu, Z., Rafikov, R. R., & Stone, J. M. 2015, ApJ, 809, L5
  • Dong (2015) Dong, R. 2015, ApJ, 810, 6
  • Dong & Fung (2017) Dong, R., & Fung, J. 2017, ApJ, 835, 38
  • Dorschner et al. (1995) Dorschner, J., Begemann, B., Henning, T., Jaeger, C., & Mutschke, H. 1995, A&A, 300, 503
  • Dullemond et al. (2003) Dullemond, C. P., van den Ancker, M. E., Acke, B., & van Boekel, R. 2003, ApJ, 594, L47
  • Dullemond & Dominik (2004) Dullemond, C. P., & Dominik, C. 2004, A&A, 417, 159
  • Edwards (2009) Edwards, S. 2009, 15th Cambridge Workshop on Cool Stars, Stellar Systems, and the Sun, 1094, 29
  • Eiroa et al. (2002) Eiroa, C., Oudmaijer, R. D., Davies, J. K., et al. 2002, A&A, 384, 1038
  • Ellerbroek et al. (2014) Ellerbroek, L. E., Podio, L., Dougados, C., et al. 2014, A&A, 563, A87
  • Espaillat et al. (2011) Espaillat, C., Furlan, E., D’Alessio, P., et al. 2011, ApJ, 728, 49
  • Espaillat et al. (2014) Espaillat, C., Muzerolle, J., Najita, J., et al. 2014, Protostars and Planets VI, 497, University of Arizona Press
  • Fairlamb et al. (2015) Fairlamb, J. R., Oudmaijer, R. D., Mendigutía, I., Ilee, J. D., & van den Ancker, M. E. 2015, MNRAS, 453, 976
  • Fedele et al. (2008) Fedele, D., van den Ancker, M. E., Acke, B., et al. 2008, A&A, 491, 809
  • Flock et al. (2017) Flock, M., Fromang, S., Turner, N. J., & Benisty, M. 2017, ApJ, 835, 230
  • Fung & Dong (2015) Fung, J., & Dong, R. 2015, ApJ, 815, L21
  • Fusco et al. (2006) Fusco, T., Rousset, G., Sauvage, J.-F., et al. 2006, Optics Express, 14, 7515
  • Gaia Collaboration et al. (2016) Gaia Collaboration, Brown, A. G. A., Vallenari, A., et al. 2016, A&A, 595, A2
  • Garufi et al. (2013) Garufi, A., Quanz, S. P., Avenhaus, H., et al. 2013, A&A, 560, A105
  • Garufi et al. (2017) Garufi, A., Meeus, G., Benisty, M., et al. 2017, A&A, 603, A21
  • Grady et al. (2009) Grady, C. A., Schneider, G., Sitko, M. L., et al. 2009, ApJ, 699, 1822
  • Grinin et al. (1994) Grinin, V. P., The, P. S., de Winter, D., et al. 1994, A&A, 292, 165
  • Grinin & Tambovtseva (2003) Grinin, V. P., & Tambovtseva, L. V. 2003, Earths: DARWIN/TPF and the Search for Extrasolar Terrestrial Planets, 539, 429
  • Haubois et al. (2016) Haubois, X., Mella, G., Duvert, G., et al. 2016, Proc. SPIE, 9907, 99073V
  • Heap et al. (2000) Heap, S. R., Lindler, D. J., Lanz, T. M., et al. 2000, ApJ, 539, 435
  • Hinkley et al. (2009) Hinkley, S., Oppenheimer, B. R., Soummer, R., et al. 2009, ApJ, 701, 804-810
  • Hubrig et al. (2009) Hubrig, S., Stelzer, B., Schöller, M., et al. 2009, A&A, 502, 283
  • Hutchison et al. (2016) Hutchison, M. A., Laibe, G., & Maddison, S. T. 2016, MNRAS, 463, 2725
  • Jordi et al. (2006) Jordi, K., Grebel, E. K., & Ammon, K. 2006, A&A, 460, 339
  • Kama et al. (2016) Kama, M., Pinilla, P., & Heays, A. N. 2016, A&A, 593, L20
  • Kenyon & Bromley (2004) Kenyon, S. J., & Bromley, B. C. 2004, AJ, 127, 513
  • Lai (2014) Lai, D. 2014, MNRAS, 440, 3532
  • Langlois et al. (2014) Langlois, M., Dohlen, K., Vigan, A., et al. 2014, Proc. SPIE, 9147, 91471R
  • Lazareff et al. (2017) Lazareff, B., Berger, J.-P., Kluska, J., et al. 2017, A&A, 599, A85
  • Lyo et al. (2011) Lyo, A.-R., Ohashi, N., Qi, C., Wilner, D. J., & Su, Y.-N. 2011, AJ, 142, 151
  • Maire et al. (2016) Maire, A.-L., Bonnefoy, M., Ginski, C., et al. 2016, A&A, 587, A56
  • Maire et al. (2017) Maire, A.-L., Stolker, T., Messina, S., et al. 2017, A&A, 601, A134
  • Marino et al. (2015) Marino, S., Perez, S., & Casassus, S. 2015, ApJ, 798, L44
  • Min et al. (2009) Min, M., Dullemond, C. P., Dominik, C., de Koter, A., & Hovenier, J. W. 2009, A&A, 497, 155
  • Min et al. (2012) Min, M., Canovas, H., Mulders, G. D., & Keller, C. U. 2012, A&A, 537, A75
  • Müller et al. (2011) Müller, A., van den Ancker, M. E., Launhardt, R., et al. 2011, A&A, 530, A85
  • Muto et al. (2012) Muto, T., Grady, C. A., Hashimoto, J., et al. 2012, ApJ, 748, L22
  • Murakawa (2010) Murakawa, K. 2010, A&A, 518, A63
  • Owen et al. (2011) Owen, J. E., Ercolano, B., & Clarke, C. J. 2011, MNRAS, 411, 1104
  • Pauls et al. (2005) Pauls, T. A., Young, J. S., Cotton, W. D., & Monnier, J. D. 2005, PASP, 117, 1255
  • Pavlov et al. (2008) Pavlov, A., Möller-Nilsson, O., Feldt, M., et al. 2008, Proc. SPIE, 7019, 701939
  • Pérez et al. (2014) Pérez, L. M., Isella, A., Carpenter, J. M., & Chandler, C. J. 2014, ApJ, 783, L13
  • Rivinius et al. (2013) Rivinius, T., Carciofi, A. C., & Martayan, C. 2013, A&A Rev., 21, 69
  • Rosenfeld et al. (2012) Rosenfeld, K. A., Qi, C., Andrews, S. M., et al. 2012, ApJ, 757, 129
  • Sitko et al. (2008) Sitko, M. L., Carpenter, W. J., Kimes, R. L., et al. 2008, ApJ, 678, 1070-1087
  • Sitko et al. (2012) Sitko, M. L., Day, A. N., Kimes, R. L., et al. 2012, ApJ, 745, 29
  • Skrutskie et al. (2006) Skrutskie, M. F., Cutri, R. M., Stiening, R., et al. 2006, AJ, 131, 1163
  • Stolker et al. (2016) Stolker, T., Dominik, C., Avenhaus, H., et al. 2016, A&A, 595, A113
  • Tambovtseva et al. (2006) Tambovtseva, L. V., Grinin, V. P., & Weigelt, G. 2006, A&A, 448, 633
  • Turner et al. (2014) Turner, N. J., Benisty, M., Dullemond, C. P., & Hirose, S. 2014, ApJ, 780, 42
  • 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
  • van der Marel et al. (2016) van der Marel, N., van Dishoeck, E. F., Bruderer, S., et al. 2016, A&A, 585, A58
  • Waters & Waelkens (1998) Waters, L. B. F. M., & Waelkens, C. 1998, ARA&A, 36, 233
  • Wisniewski et al. (2008) Wisniewski, J. P., Clampin, M., Grady, C. A., et al. 2008, ApJ, 682, 548-558
  • Woitke et al. (2016) Woitke, P., Min, M., Pinte, C., et al. 2016, A&A, 586, A103
  • Zubko et al. (1996) Zubko, V. G., Mennella, V., Colangeli, L., & Bussoletti, E. 1996, MNRAS, 282, 1321
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
Loading ...
This is a comment super asjknd jkasnjk adsnkj
The feedback must be of minumum 40 characters
The feedback must be of minumum 40 characters

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 description