Open Access
Issue
A&A
Volume 688, August 2024
Article Number A149
Number of page(s) 23
Section Interstellar and circumstellar matter
DOI https://doi.org/10.1051/0004-6361/202348555
Published online 13 August 2024

© The Authors 2024

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

This article is published in open access under the Subscribe to Open model. Subscribe to A&A to support open access publication.

1 Introduction

Disk evolution seems inherently connected to the formation of planets in circumstellar disks around low-mass pre-main-sequence (PMS) stars. Some of the big questions are when, where, and how planets, in all their observed varieties, exactly form. A thorough understanding of disk evolution is key to understand the onset of planet formation in these environments.

An important aspect of the evolution of a disk is how matter is distributed, for instance how dust grains with different sizes are distributed vertically, radially and within substructures. In recent years, major progress has been made using telescopes and instruments that are able to probe the mid-plane structure and surface of the disk (e.g., ALMA and VLT/SPHERE). Extensive observations show great variety in disk shapes (e.g., Avenhaus et al. 2018; Law et al. 2023) and substructures in the form of cavities, rings, spiral arms, and arcs (e.g., Long et al. 2018; Andrews et al. 2018; Garufi et al. 2020; Andrews 2020; Benisty et al. 2023; Bae et al. 2023).

These substructures are thought to trace active sites in the disk, where planetesimals may grow and act as an indication of planet formation (Andrews 2020). This makes the study of these substructures important to uncover disk evolution and diskplanet interactions. For example, Pinilla et al. (2018) studied transition disks with cavities and identified giant planet formation and dead zones as potential explanations for the observed properties of these disks. Another large study of disk substructures from Garufi et al. (2020) with ALMA and SPHERE shows that rings and cavities detected with ALMA do not always have a similar SPHERE counterpart, demonstrating the importance of combined and complementary observations.

Disks around PMS stars have limited time to form planets since they are thought to dissipate rapidly. Large surveys of T Tauri stars have shown that only a few percent of the disks survive beyond an age of 10Myr (Padgett et al. 2006; Carpenter et al. 2009; Wahhaj et al. 2010). However, recently so-called “Peter Pan Disks” have been discovered around M-type stars. These systems have accretion disks older than 20Myr that “refuse to grow up”. These disks are thought to be long-lived due to a high initial disk mass with a low turbulence parameter and a lack of external photo-evaporation (Sicilia-Aguilar et al. 2011; Silverberg et al. 2020; Coleman & Haworth 2020).

In this paper, we combine SPHERE and ALMA observations with medium-resolution spectroscopy of the star and the hot inner disk in order to infer a detailed picture of the system PDS 111. This TTauri star, often classified as G-type (Rojas et al. 2008), is known to be in the foreground of the Orion star-forming region at a distance of 158 pc (Gaia Collaboration 2022) with a relatively high age (15–42 Myr, Rojas et al. 2008; Bell et al. 2015; Launhardt et al. 2022). However, this system seems to be one of these rare cases where the circumstellar disk has survived up until this significant age. The aim of this work is to reassess the age determination of the system and to measure its disk properties in order to understand how PDS 111 fits in the current understanding of disk evolution, and how this relates to the long-lived disks in the “Peter Pan Disk” scenario.

In the following, we describe the observations in Sect. 2 and the data analysis of the photospheric light, scattered light, dust continuum and 12CO emission in Sect. 3. The age of the system and the disk properties are presented in Sect. 4. We discuss the results in Sect. 5 and summarize our findings in Sect. 6.

2 Observations

PDS 111 was observed in the framework of the Disk Evolution Study Through Imaging of Nearby Stars (DESTINYS) ESO large program (Ginski et al. 2020, 2021). Resolved high-resolution imaging observations were carried out with the extreme adaptive optics instrument VLT/SPHERE (Beuzit et al. 2019). Additionally, we obtained sub-mm observations with ALMA, medium-resolution spectroscopic observations of the central star with VLT/X-shooter (Vernet et al. 2011), high-resolution spectroscopy with Mercator/HERMES (Raskin et al. 2011) and we used photometric observations from the Transiting Exoplanet Survey Satellite mission (TESS; Ricker et al. 2015).

2.1 Spectroscopy

2.1.1 X-shooter medium resolution spectroscopy

Spectroscopic observations in the optical to near infrared (NIR) (300–2500 nm) have been taken with X-shooter, mounted on the UT2 of the VLT, with program ID 108.2216.001. PDS 111 was observed with X-shooter during two separate epochs on November 28, 2021 and December 7, 2021. The first observation has been taken with a poor seeing of ~2.6″ and the second under good seeing conditions of ~0.7″. Both observations can be used in the study. The spectra are taken in nodding mode. The used slit widths are 0.5″ in the UVB arm leading to a resolving power of 9700, and 0.4″ in the VIS and NIR arms, which corresponds to a resolving power of 18 400 and 11600, respectively.

The spectroscopic data were reduced with the ESO X-shooter pipeline 3.3.5 provided by ESO (Modigliani et al. 2010). The pipeline consists of a flat-field correction, bias subtraction, and wavelength calibration. Spectrophotomeric standards from the ESO database were used for flux calibration. The reduction of the nodding spectra was done with the xsh_scired_slit_nod recipe provided in the software package from ESO. The Molecfit 1.5.9 tool was used to correct for the telluric sky lines in the spectra (Smette et al. 2015; Kausch et al. 2015).

The arms are normalized separately using second-order spline functions. The continuum points for the spline fit are selected at the same wavelength intervals for each epoch.

2.1.2 HERMES high resolution spectroscopy

High-resolution follow-up observations have been taken with the HERMES spectrograph (380–900 nm) mounted on the Mercator 1.2 m telescope on La Palma (Raskin et al. 2011) to obtain a longer temporal baseline for line variability. HERMES has a spectral resolving power of R ~ 85,000. Observations were taken on two nights, October 20, 2022 and October 28, 2022. The data were reduced with the standard HERMES pipeline, consisting of a flat-field correction, bias subtraction, wavelength calibration, and cosmic removal. The humidity was relatively high (~50%) for the second observation compared to the first observation, introducing stronger telluric features. The spectra were rebinned on the same resolution as X-shooter to increase the signal-to-noise ratio and similarly normalized in order to compare the spectra of HERMES and X-shooter.

2.2 SPHERE/IRDIS near-infrared imaging

PDS 111 was observed with the InfraRed Dual-band Imager and Spectrograph (IRDIS, Dohlen et al. 2008) subsystem of SPHERE on December 30, 2020. The observations were carried out in the BB_H filter with a central wavelength of 1.65 µm. We used the dual-beam polarimetric imaging mode of the instrument (de Boer et al. 2020; van Holstein et al. 2020) in pupil-tracking mode (van Holstein et al. 2017). The central star was placed behind an apodized Lyot coronagraph with an effective radius of 92.5 mas to avoid saturation and to limit photon noise. We recorded a total of 104 science frames split into 26 polarimetric cycles. The individual frame exposure time was 32 s resulting in a total exposure time of 55.5 min. The atmospheric conditions during the observation sequence were excellent with an average seeing of 0.72″ and a coherence time of 9.8 ms. Bracketing the science observations we acquired dedicated astrometric star-center and flux calibration frames, as well as sky calibration frames.

We reduced the data using the IRDIS Data reduction for Accurate Polarimetry1 pipeline (IRDAP, van Holstein et al. 2020), which performs polarimetric differential imaging (PDI, Kuhn et al. 2001) to remove the light from the bright central star. We used the default setup parameters. IRDAP first performs image pre-processing, which includes sky subtraction, flat-fielding, and bad-pixel masking. A static master flat-field image, as well as a static badpixel mask were used, which are included in IRDAP. Then the frames were separated in the left and right detector side and individually centered, using the astro-metric star-center calibration frames. After these pre-processing steps IRDAP combines images from both detector sides, as well as images at different half-wave plate position to compute the Stokes Q and U images. To remove residual instrumental polarization a full Mueller matrix model of the telescope and instrument is implemented that describes the instrumental polarization effects. Additionally, residual (inter)stellar polarization is removed by measuring the signal at the adaptive optics correction radius within the Q and U images and subtracting a scaled total intensity image. The final Q and U images are then converted into the azimuthal Stokes parameters Qϕ and Uϕ, following the definition in de Boer et al. (2020). The Qϕ image contains all azimuthally polarized light as positive signal. Since this is the expected orientation of the polarization direction for single scattering events it gives a good representation of the disk polarized light image. We show the Qϕ image in Fig. 1 and the full set of Stokes images in Fig. A.1.

In addition to the PDI reduction, we performed angular differential imaging (ADI, Marois et al. 2006) post-processing with the IRDAP pipeline. We used the principal component analysis approach to fit the reference catalog (consisting of all science images themselves) to each science frame. For the processing the default settings of IRDAP were utilized. The ADI processed total intensity image is likewise shown in Fig. 1. For this image we specifically used 4 principal components and optimization annuli between 10 and 30 pixels (0.12” to 0.37”), 30 and 100 pixels (0.37” to 1.23”), and 100 and 512 pixels (1.23” to 6.27”).

thumbnail Fig. 1

SPHERE/IRDIS H-band (1.65 µm) observations of PDS 111. In all images the star is at the center and the images are aligned with the North-axis up and East-axis to the left. Left: polarized light observation. We show the azimuthal Stokes parameter Qϕ which contains all azimuthally polarized scattered light as a positive signal. Middle: same data set as in the left panel, but reduced with angular differential imaging, showing total intensity. Due to the reduction approach, the disk signal is strongly self-subtracted removing particularly low spatial frequencies. The shape of the detected high spatial frequency sub-structure is likely distorted. For the same reason, the given surface brightness is a lower limit only, which depends in a complex way on the local self subtraction in the image. Right: same image as in the left panel, overlayed with contours (5 σ, 14 σ, 68 σ, 140 σ, 205 σ; relative to the image background noise) from the center panel to highlight the structures that are traced in both images. In all images, the gray, hashed circle in the center marks the area that was covered by the coronagraphic mask and where the star is located. All images are shown on a log scale to highlight the disk morphology. Appendix A shows a flux calibrated version of the polarimetric image.

2.3 ALMA observations

The 1.25 mm emission of PDS 111 was observed with ALMA Band 6 as part of the ALMA project 2021.1.01705.S (PI: C. Ginski), with 1.31 min of exposure on August 28, 2022 and 1.31 min in September 2, 2022. The correlator was configured to observe four spectral windows: one included the 12CO molecular line at 230.538 GHz in the J = 2 − 1 transition with a frequency spacing of 488.3 kHz (~0.65 km s−1 per channel), while the remaining three were targeted at dust continuum emission centered at 233.519 GHz, 246.020 GHz, and 248.021 GHz, with a channel frequency spacing of 15.625 MHz and a total bandwidth of 2 GHz. The antenna baselines span 15.1 m to 1210.6 m.

Executing the scriptforPI provided by ALMA, we obtain the pipeline calibrated data, which is self-calibrated using dust continuum emission observations, so that the signal-to-noise ratio can be improved. Using CASA 5.6.2 we flagged the channels located at ±20 km s−1 from the line center to obtain the dust continuum emission. A pseudo-continuum measurement set is created by combining the remaining channels with the continuum spectral windows. To reduce data volume, the set was averaged into 125 MHz channels and ‘6s’ bins. The self-calibration was applied on the continuum by imaging with the CLEAN algorithm, with a Briggs robust parameter of 0.5. The robust parameter is used to control the trade-off between the sensitivity and angular resolution of the data by applying a natural and/or uniform weighting scheme to the data. The value of 0.5 slightly favors angular resolution over sensitivity (Briggs 1995; Boone 2013). We applied four phase calibrations, using solution intervals of length “inf”, “40s”, “18s”, and “6s”, while the amplitude calibration was done with a solution interval of “6s”, combining spectral windows and scans to find the solution. The correction tables were later applied to the 12CO spectral window.

After self-calibration, the continuum emission was imaged with a robust parameter of −0.5 to maximize the balance between angular resolution and sensitivity. The beam size of the 12CO is 0.429″. The gas emission was imaged with a robust parameter of 0.0 and no visibility tapering to focus on best possible angular resolution. Both images were cleaned until the 4σ noise level to apply the JvM correction (JvM “effect” from Jorsater & van Moorsel 1995), which takes in account the volume ratio (ϵ) values of the point spread function of the images and the restored Gaussian from the CLEAN algorithm following the method as described in Czekala et al. (2021). We found e values of 0.909 for the dust continuum image and 0.694 for the 12CO, which means that the dirty beam shape for the dust continuum is more closely described by an elliptical Gaussian, while that for the 12CO imaging deviates more significantly.

2.4 TESS photometry

TESS observed PDS 111 during its runs starting in November and December 2018 (target name 24697724). The target was observed with exposures of 120 s during 42 days (spread between November, December, and January) with the TESS detector bandpass that spans from 600 to 1000 nm.

thumbnail Fig. 2

Stacked spectrum of PDS 111 shown together with a subset of Kurucz models with different temperatures. The models have a log 𝑔 of 4.2 and a υ sini of 21 km s−1.

3 Data analysis

3.1 Photospheric spectrum PDS 111

Based on the strength of the G-band, Fe I lines and Hy around 432 nm, the star is classified as early G-type (Gray 2000). A comparison of the spectrum with stars in the X-shooter library indicates a spectral class between G1 and G4. The luminosity class V is based on the relative strength of MnI 403.0nm and Fe I 404.6 nm and on the ratio between Ti II, Fe II 417.2−8 nm, and Ti II, Fe II 444.4 nm (see Fig. 2).

The temperature, surface gravity, and rotational velocity were constrained through the fitting of Kurucz models (Kurucz 1993) in the spectral range of 402–430 nm with varying grid parameters based on the spectral classification; Teff (5500–6500 K in steps of 100 K), log 𝑔 (4.0–4.5 in steps of 0.1), and υ sin/ (20–50 km s−1 in steps of 1 km s−1). The best fit model is with a Teff of 5900 K, log 𝑔 of 4.2, and v sin/ of 21 km s−1; see Appendix B for the X minimization results and Fig. 2 for the best fit and sensitivity of the spectrum to changes in effective temperature. The Teff value is typical for G2-type pre-MS stars (Pecaut & Mamajek 2013) and the gravity is in agreement with luminosity class V.

Similar to Woitke et al. (2016) and Valegård et al. (2021), the Kurucz-model spectrum was fitted to optical photometry to determine the star’s luminosity, radius, extinction parameters, and log 𝑔 from the spectral energy distribution (SED) fit. The results are reported in Table 1. The optical photometry is known to vary (see e.g., Torres et al. 1995; Rojas et al. 2008; Zacharias et al. 2009; da Silva et al. 2009; Kiraga 2012); therefore, we selected the mean V-band measurement of the ASAS-SN survey with 279 measurements over a period of five years in order to account for the variations (∆V ~ 0.43 mag, Jayasinghe et al. 2018). A similar B-band study for PDS 111 is lacking. We used Gaia magnitudes to inspect the behavior of B–V, applying Gaia photometry conversion tables from Busso et al. (2022). The mean B–V is 0.729 mag and varies with 0.075 mag. Although these variations are significant compared to the average errors on the B–V (<1%), we neglected color variability for our use because the absolute variations of the V-band were substantially higher. Therefore, we used the mean B–V from Gaia to obtain a B-band photometric point for the SED fit. For the fitting procedure, we fixed the distance to the Gaia distance from DR32 (Gaia Collaboration 2022) and effective temperature. The SED fit can be found in Appendix C, which shows an infrared excess compared to the stellar SED. We might observe excess flux compared to the bluest photometric point from SDSS in the u-band (λeff = 360.804 nm), but more photometry blueward from the u-band is needed to confirm this. We used the variability of the V-band magnitude in the ASAS-SN survey, extinction, effective temperature, distance, and bolometric correction from Pecaut & Mamajek (2013) to infer errors on the luminosity and radius.

The luminosity and Teff of the star result in a radius of 1.52 R. Alternatively, the radius was derived from the υ sini, rotation period, and inclination of the star. The TESS photometry displays a dominant period with the use of a Lomb-Scargle peri-odogram (Lomb 1976; Scargle 1982) of 3.137 days. This period is often attributed to the rotation period of the star. Using this period, the υ sini of 21 km s−1 from Kurucz model fitting and an inclination of 58° (see Sect. 3.3), we calculated a radius of 1.53 R, which is within errors similar to the derived radius from the SED fit.

Table 1

Stellar parameters of PDS 111.

3.2 Scattered light observations

Our SPHERE H-band observations reveal an extended disk seen in polarized scattered light. We traced disk signal out to ~1″ from the central star along the major axis of the disk, corresponding to ~158 au. We show the SPHERE polarized and total intensity reductions in Fig. 1. The polarized light image shows a classical flaring disk with an illuminated top side, a dark mid-plane and an illuminated outer edge of the disk bottom side (going from North to South in the image). The disk is seen at an intermediate inclination. Tracing of the outer edge of the top side of the disk by eye yielded an inclination of ~60° at a position angle (PA) of the major axis of ~68°.

The total intensity image in Fig. 1 traces sub-structures on the disk top side. The employed angular differential imaging processing removed most of the signal with low spatial frequencies due to image self-subtraction (see e.g., Milli et al. 2012; Ginski et al. 2016; Stapper & Ginski 2022). We identified four pairs of arc-like structures, which indicate the presence of rings on the disk surface. As discussed in Ginski et al. (2016) the shape of these arcs is likely distorted due to the same image self-subtraction effects that filtered out the low spatial frequency data. In the rightmost panel of Fig. 1 we overlayed the contours of the total intensity image on the polarized light image which does not suffer from distortions of the disk morphology. The brightest pair of arcs in total intensity is tracing the outer edge of the disk top-side, while the two pairs of arcs further North trace rings in the inner part of the disk, also visible in the polarized light data. The faintest and southernmost pair of arcs in total intensity is tracing the edge of the disk mid-plane.

thumbnail Fig. 3

ALMA observations of PDS 111. Left panel: dust continuum emission at 1.25 mm of PDS 111. The resolution of the observations is 0.357” × 0.25” with a position angle of −71.9° and is shown in the lower left part of the figure. Right panel: overlap of ALMA (contours at 10,…,100% of the peak of emission) and SPHERE observations (in color).

3.3 Dust and line observations

3.3.1 Dust continuum emission

Figure 3 shows the dust continuum emission observed with ALMA with a final beam of 0.357” × 0.25” with a PA of −71.9° at 1.25 mm after self-calibration as described in Sect. 2.3. To check if the emission is optically thin, we calculated the optical depth of the peak of the continuum by assuming (e.g., Dullemond et al. 2018): Iν(r)=Bν(Td(r))(1exp[ τν(r) ]) thusτ=ln(1I(rpeak )B(Td(rpeak ))),$\matrix{ {{I_\nu }(r) = {B_\nu }({T_{\rm{d}}}(r))(1 - \exp \left[ { - {\tau _\nu }(r)} \right])\quad {\rm{thus}}} \hfill \cr {\tau = - \ln \left( {1 - {{I\left( {{r_{{\rm{peak}}}}} \right)} \over {B\left( {{T_d}\left( {{r_{{\rm{peak}}}}} \right)} \right)}}} \right),} \hfill \cr } $(1)

with Td being the dust temperature at the peak location (rpeak). Assuming a dust temperature of 20 K, the optical depth at the peak is τpeak = 0.27, similar to the values found in the DSHARP sample (Huang et al. 2018; Dullemond et al. 2018). This value of the optical depth may be underestimated because dust scattering could be neglected (Zhu et al. 2019), which would result from unresolved substructures in the ring, and/or a lower dust temperature or the large beam size relative to the disk.

If the emission is optically thin, the dust disk mass can be calculated as Mdust d2FvκvBv(T),${M_{{\rm{dust}}}} \simeq {{{d^2}{F_v}} \over {{\kappa _v}{B_v}(T)}},$(2)

where d is the distance to the source, Fv is the total flux at 1.3 mm, and Bv is the black-body surface brightness at a given temperature (Hildebrand 1983). Taking for the mass absorption coefficient κv = 2.3 cm2 g−1 × (v/230GHz)0.4 (Beckwith et al. 1990; Andrews et al. 2013), we calculated a dust disk mass of 45.8 M or ~ 0.14 MJup, similar to typical disks around T Tauri stars in different star-forming regions (e.g., Pascucci et al. 2016).

3.3.2 12CO emission

Figure G.1 presents the channel maps of 12CO of PDS 111 observed with ALMA, showing emission from ~0.4 km s−1 to ~ 14.7 km s−1. For comparison, the 5σ level of the continuum emission is overplotted, and it shows that 12CO is more extended than the dust continuum emission and the scattered light detection.

Figure 4 shows the moment 0, 1, and 8 (peak value of the spectrum) of the 12CO of PDS 111 generated with a Keplerian mask and clipping at 3σ of the JvM corrected image as described in Sect. 2.3 (the beam size is 0.429″ or 0.342″ with a PA of −76°.114). To extract the stellar mass from the 12CO Keplerian rotation pattern, we used eddy3 (Teague 2019), which gives 1.290.18+0.82$1.29_{ - 0.18}^{ + 0.82}$ M when assuming a simple disk model and fixing the inclination to be that of the continuum. The large uncertainty in the stellar mass may come from the unknown vertical structure and disk orientation. Higher angular resolution is needed to better resolve the surface height with ALMA.

thumbnail Fig. 4

Moment 0, 1, and 8 of the 12CO of PDS 111. Moment 0 (left) is the integrated value of the spectrum, moment 1 (middle) the intensity weighted coordinate to get the velocity dispersion and moment 8 (right) shows the peak value of the spectrum. The beam size is 0.429″ × 0.342″ with a PA of −76°.114, shown as the oval in the lower left corner of the panels.

4 Age and disk properties

By combining the observations at different wavelengths, we are able to sketch a more complete picture of the system. Figure 5 shows what different parts of the disk are traced by the different observing techniques.

4.1 Age

Using the luminosity and effective temperature, we placed the star in the Hertzsprung Russell Diagram (HRD) on evolutionary tracks from Siess et al. (2000), and infer an age of 15.9 Myr and a evolutionary mass of 1.20 M, see Fig. 6. The error on the effective temperature and conservative error on the luminosity lead to an uncertainty range of 12.2−17.6 Myr for the age and 1.15−1.30 M for the mass. The mass is similar to the best fit mass value obtained in Sect. 3.3. Additionally, using the radius and surface gravity of the star in Table 1, we measured the spectroscopic mass. Using the log 𝑔 from the SED fit resulted in a mass of 1.14 M and using the log 𝑔 from the spectrum fit resulted in a mass of 1.34 M. Both are within the errors consistent with the evolutionary mass.

Manara et al. (2013) show that the age of T Tauri stars could be over-estimated when the extinction, accretion, and PMS stellar models are not fitted simultaneously. However, their method cannot be applied to PDS 111 since their stellar models go up to spectral type G4. The accretion rate is measured from line luminosities in Sect. 4.6 demonstrating that the accretion rate for PDS 111 is low for a star of this mass (Manara et al. 2023).

The estimated age of the system of 15.9 Myr could be wrong if the distance to PDS 111 obtained from the Gaia DR3 parallax (Gaia Collaboration 2022) is not correct. This may be the case considering the large RUWE value of 2.6, despite the small error on the parallax (0.035 marcsec). The RUWE statistic shows whether a single-star light profile fits the observations. The circumstellar material of PDS 111 might inflate its value. Fitton et al. (2022) determine that for disk-bearing sources a RUWE cutoff value of 2.5 is more appropriate. The RUWE value of PDS 111 is slightly above this cutoff4. We checked the Gaia-database to see whether the visibility_period_used is more than 10 and ipd_gof_harmonic_amplitude is less than 0.1. Both is the case: the former is 22 and the latter is 0.029. Moreover, the Gaia parallax cross-matches that obtained with HIPPARCOS (Gaia Collaboration 2022; ESA 1997). Therefore, we do not find any obvious indications that the Gaia measurements are inaccurate.

Throughout the SED fitting procedure, the total-to-selective visual extinction RV is fixed to 3.1 (Cardelli et al. 1989). However, it has been noticed that RV is sometimes higher toward star-forming regions (e.g., Ramírez-Tannus et al. 2017). If this value would for example be 4.1, the extinction toward the source would be higher and the luminosity would increase to a value of 3.6 L and the corresponding age in the HRD would be 11.3 Myr for the same temperature, lowering the age of the system. More simultaneous photometric measurements over a period of time could help to improve the SED fit and fit for a range of different RV values in order to find the best fitting parameter in the sightline of PDS 111. We conclude that the uncertainty in distance or extinction does not reject our hypothesis that PDS 111 is relatively old, that is at least ~11 Myr. This relatively old age is similar to what has been determined for the stars that are surrounded by a “Peter Pan Disk”; we further discuss this scenario in Sect. 5.6.

4.2 Disk surface

As discussed in de Boer et al. (2016), Ginski et al. (2016), and Avenhaus et al. (2018), we used the identified disk sub-structures to extract the height profile of the disk surface layer at which the disk becomes optically thick for NIR radiation. To facilitate the tracing of the structures we employed high pass filtering as well as deconvolution techniques using the Adaptive Image Deconvolution Algorithm (AIDA) code by Hom et al. (2007). We show the high-pass filtered and the de-convolved data in Fig. E.1. Because the disk structure is complex, we traced the individual rings by eye as was done by Avenhaus et al. (2018) for the case of IM Lup. Tracing was done for all structures in the original polarized light images, in the high-pass filtered images and in the de-convolved images to ensure that there was no mismatch in structures. To reduce the number of free parameters, the inclination of the individual rings as well as their PA were fixed to the values extracted from the ALMA data discussed in Sect. 3.3. Thus, for each ring we traced an ellipse by re-scaling to match the semi-major axis and then applying an offset along the disk minor axis to match the height profile of the disk. Our results are shown in Fig. 7 (top). Additionally we show a saturated image of the disk for the overlay of the ellipse on the outer disk edge in Fig. E.2.

We overlay the resulting ellipses on the polarized light image and calculate the disk aspect ratio from the ellipse center-offset along the minor axis and the ellipse semi-major axis following de Boer et al. (2016). We find a steep increase in the disk aspect ratio from the innermost ring at 37 au to the second ring at 52 au from 0.12 to 0.26. While the aspect ratio of the inner ring is compatible with results presented by Avenhaus et al. (2018) for a small sample of T Tauri stars, the subsequent structures all show a higher aspect ratio, indicative of a strong flaring geometry of the disk, see Fig. 7 (bottom). In addition, to the two rings at 37 au and 52 au we identify a third more diffuse ring-structure in the deconvolved and high-pass filtered images (see Fig. E.1) at a separation of 74 au. Since this structure is of lower spatial frequency than the inner two rings it is entirely filtered out by the ADI post-processing and thus not visible in the total intensity image in Fig. 1. We tentatively overlay an ellipse on this ring as well and find an aspect ratio of 0.23, roughly consistent with the inner ring at 52 au, however, with larger uncertainties.

The increase of the aspect ratio between the innermost ring at 38 au and the outer disk edge at 140 au, for which we found a value of 0.26, is in principle compatible with the extreme flaring profile found for the disk around the Herbig star HD 97048 by Ginski et al. (2016). However, since the aspect ratio at intermediate separations from the central star appears to be significantly higher, we could not fit the entire surface height profile with a single power law. As briefly discussed in Ginski et al. (2016), this may be due to an optical depth effect, disproportionately affecting the outer disk region. As the dust density in the outer disk decreases, the optically thick layer sinks closer to the mid-plane.

thumbnail Fig. 5

Overview of PDS 111 indicating the parts of the system that are traced by the different observational techniques.

thumbnail Fig. 6

Hertzsprung-Russell diagram with PMS tracks and isochrones from Siess et al. (2000). The location of PDS 111 coincides with an age of 15.93.7+1.7$15.9_{ - 3.7}^{ + 1.7}$ Myr.

thumbnail Fig. 7

Ellipses fitted to the radial sub-structures traced in scattered light (top) and resulting height of the scattering surface of the disk (bottom). We compare the scattering surface height with the average surface height profile of several T Tauri stars found by Avenhaus et al. (2018) and the extreme flaring case of the Herbig star HD 97048, discussed by Ginski et al. (2016).

thumbnail Fig. 8

Scattered light disk flipped around the minor axis and divided. The minor axis is indicated by the white, dashed line. The scattered light disk is asymmetric relative to the minor axis. The colorbar gives the flux ratio between the two sides of the disk, showing that the flux deviates up to a factor of two between both sides. Notably, the flux ratio turns around for the region inside of ~0.3″ along the major axis. This may indicate a local change in scale height symmetric around the disk midplane rather than a warp or misalignment since the bottom side displays a similar flux ratio.

4.3 Disk asymmetry

For a disk with an azimuthally homogeneous dust and gas distribution, we would expect that the polarized light image is symmetric with respect to the disk’s minor axis. However, this does not appear to be the case for PDS 111. To highlight the disk asymmetry we flipped the image around the minor axis and then divided the original image by the flipped image. This highlights the areas in which both sides of the disk deviate. We show the result in Fig. 8. We see two main effects; first of all the disk is indeed asymmetric with flux deviations of up to a factor two between both sides. Second, we find a flip in the asymmetric brightness distribution (the white line between the red and blue areas) around 0.3″ (47 au) along the major axis. Inside of this border the North-East (NE) side of the disk is fainter than the South-West (SW) side, while it is the opposite way in the outer disk.

Since the SPHERE observations are tracing scattered stellar light, the disk brightness is intricately linked to the disk height profile and geometry. Broad brightness asymmetries as observed here can be caused by warps or misalignment in inner disk regions and subsequent shadowing of the outer disk (e.g., Debes et al. 2017; Benisty et al. 2018; Muro-Arena et al. 2020). In such a scenario the breaking or warping point would be at the location where we observe the flip in the asymmetric signal of the disk. In that case, the scale height inside of the flipping point would be increased in the SW on the disk top and decreased in the NE. This could then explain simultaneously why the disk is brighter inside the flipping point (it is more exposed to stellar light) and fainter outside the flipping point (the outer disk is shadowed by the inner disk) in the SW.

However, we note that in a warping scenario we would expect that the disk bottom side shows the opposite behavior from the disk top side, that is, we would expect the disk bottom side to be bright in the SW. While the S/N is lower on the disk bottom, it nevertheless appears that it instead exhibits the identical behavior as the disk top side, which means it is faint in the SW and bright in the NE. This may then indicate that the disk is not radially warped or misaligned as a whole but rather shows a local change in scale height, symmetric around the disk mid-plane. This may be compatible with the fact that we are measuring an increased aspect ratio for the disk ring at 52 au. However, since we traced the ring at all azimuthal angles simultaneously a localized effect may be averaged out.

thumbnail Fig. 9

Observations versus model images from our best fit with galario after convolving with the same beam as the observations and the residuals. The contours of the residual map are [−4σ, −2σ, 2σ, 4σ], where negative values are dashed contours and positive values solid contours.

4.4 Disk structure

The dust continuum emission morphology was analyzed with a parametric brightness modeling using the packages galario (Tazzari et al. 2018) to calculate the Fourier Transform of the models and zeus (Karamanis & Beutler 2021; Karamanis et al. 2021) to run an MCMC and sample the posterior likelihood of the parameter space. We considered four components to describe the PDS 111 emission: A centrally peaked Gaussian (𝑔1) motivated by the inner emission found in several single-ringed disks (e.g., Pinilla et al. 2021; Pyerin et al. 2021), two Gaussian profiles for the ring (𝑔2, 𝑔ı) to allow for radially asymmetric ring morphology (rather than one broken Gaussian profile), and a Gaussian azimuthal asymmetry (𝑔4) to describe the over-emission observed in the South-western side of the disk. Figure F.1 illustrates each of the components of the model used to fit the morphology of PDS 111. We also include as free parameters the central position of the disk (x0, y0) and the geometry, given by inclination and PA. Every parameter has a flat prior and boundaries representative of physical limitations, such as only allowing positive distances and brightness. The model images have a pixel size of 4 mas and an image size of 2048 pixels. We ran the MCMC with eight times the number of parameters as the number of walkers (nwalk = 136). For each walker, the burn stage was 5000, and 6000 for the parameter space sampling. The summary of the modeling results is given in Table 2.

Figure 9 shows the comparison between the dust continuum observations and model images from our best fit with galario after convolving with the same beam as the observations and the residuals. The residual map demonstrates that our model recovers most of the morphology of the dust continuum emission. The main characteristic of the dust morphology is a cavity with a ring-like emission peaking at ~30 au (see Fig. F.1 for the model components). The radius that encloses 90% of the total flux from the visibility fitting is ~55 au.

Table 2

Dust continuum visibility modeling results.

4.5 Radial extent of the disk

The radial distribution and extension of PDS 111 in dust continuum, scattered light, and 12CO is shown in Fig. 10. For the ALMA profiles, the images have been deprojected and then azimuthally averaged, while for the SPHERE data we show radial cuts at the North-East (NE) and South-West (SW) of the major axis of the disk. The vertical lines show the radius that encloses 90% of the emission in each of the maps. The 90% radius of the 12CO map (~295 au) is much more extended than that tracing the dust, both in scattered light and continuum emission. The 90% radius of the dust continuum emission is at 95 au, three times smaller than observed in12 CO (both profiles being affected by convolution of a similar beam). The latter is very similar to the average ratio observed in several protoplanetary disks that are on average younger (Long et al. 2022).

The radius that enclose 90% of the flux measured directly from the dust continuum emission (95 au) and that measured from the visibility fitting of the dust continuum emission (55 au, Sect. 4.4) are different. This is a natural result of the large ALMA beam and the high inclination of the source. The ~55 au estimate is probably more representative of the actual size, therefore of the distribution of mm-grains that dominate the emission at ALMA wavelengths. Figure 10 shows that this mm-grain thermal emission disk is smaller than the disk size from scattered light (~ 105 au for the NE and ~96 au for the SW cuts), which would indicate that the small grains responsible for the scattering are more radially extended than the large grains, as a product of radial drift (e.g., Guilloteau et al. 2011; Tazzari et al. 2016; Trapman et al. 2019, 2020).

When the radial profiles from scattered light are multiplied by r2 to compensate for the stellar illumination, the 90% radius shifts to ~ 155 au for the NE and ~ 136 au for the SW emission (in line with the 158 au mentioned in Sect. 3.2 at which the disk signal drops below the noise floor along the major axis). The size of the disk in scattered light is thus ~2.4−2.7 times larger than in the millimeter emission. Such a difference is similar to that found in the ALMA images of dust vs. gas emission. Trapman et al. (2019, 2020) propose that optical depth effects cannot explain such a big difference, and this is likely evidence for small and large grain separation through radial drift.

thumbnail Fig. 10

Azimuthally averaged radial intensity proflies of the deprojected images from ALMA dust continuum and the moment 0 maps of the 12CO from ALMA Band-6 observations of PDS 111. Each profile is normalized to the peak. The ALMA profiles include the standard deviation of each elliptical bin divided by the square root of the number of beams spanning the full azimuthal angle at each radial bin. For the SPHERE observations, we take two radial cuts at the North-East (NE) and South-West (SW) of the major axis of the disk. The vertical lines show the radius that encloses 90% of the emission in each of the maps. The gray region shows the size of the coronagraph in the SPHERE observations.

Table 3

Derived mass-accretion rates of Ha and Hβ at each epoch in logṀacc [M yr−1].

4.6 Inner gaseous disk

The presence of hot gaseous circumstellar material is revealed by Hα and Hβ emission and strong variability in these and other H I, He I, Ca II and Na I lines. We do not detect any forbidden emission lines. The equivalent width (EW) is determined after subtracting the best fit Kurucz model from the observations. This ensures that the measured EW is an indication of the emission or absorption from the circumstellar material. The EW measurements are listed in Table D.1. The EW (Hα) emission is < 10 Å, which classifies PDS 111 as a weak-line T Tauri star.

Following Alcalá et al. (2017), we determined line luminosities by measuring the EW of the lines, the absolute flux from the SED in Ha and Hβ, and converting this to luminosity using the distance to the star. Their empirical relations were used to convert the line luminosities to accretion luminosities and the resulting mass-accretion rates assuming a typical inner radius of 5 R, are reported in Table 3. This yields a rate that varies between 1−5 × 10−10 M yr−1. Due to the lack of photometric observations in the UV, the mass-accretion rate could be constrained only indirectly.

The temporal variance method (TVS; Fullerton et al. 1996; de Jong et al. 2001) is employed to study spectroscopic variability. The observed standard deviation σobs is compared to the expected standard deviation σexp at a certain wavelength bin j. If the ratio of the two standard deviations exceeds unity, the variation in the spectral line between the epochs is larger than what is expected from the noise in that wavelength bin. Likewise, when no variation is detected, the ratio (i.e., the value of the TVS) is randomly distributed around unity, reflecting the noise in the data.

Similar to de Jong et al. (2001) and Derkink et al. (2024), σobs, j is determined from the average of the standard deviations of the bins in the continuum region around the spectral line in all epochs. This yields σcon, which is scaled by the ratio of the mean normalized flux of all epochs of the evaluated wavelength bin j and the continuum. The σobs,j is calculated from observed fluxes in wavelength bin j. The TVS value is given by the ratio: TVSj=σobs,jσexp,j=σjσconF¯conF¯j.${\rm{TV}}{{\rm{S}}_j} = {{{\sigma _{{\rm{obs}},j}}} \over {{\sigma _{{\rm{exp}},j}}}} = {{{\sigma _j}} \over {{\sigma _{{\rm{con}}}}}}\sqrt {{{{{\bar F}_{{\rm{con}}}}} \over {{{\bar F}_j}}}} .$(3)

Here _F¯j${\bar F_j}$ is the average flux in wavelength bin j for all epochs and F¯con ${\bar F_{{\rm{con}}}}$ the average continuum flux. The significance of the variability is determined by comparing a χ2-distribution with the degrees of freedom set to the number of observations, to the TVS value. A 95% level is used for a significant detection.

The resulting TVS for all significantly varying lines are shown in Fig. 11 and were used to measure the maximum velocity at which significant variability is detected. The lines that vary in a similar manner are shown together. The highest velocities and strongest variations are seen in Hα and He I, up to 320 km s−1 on the red side. The blue side shows variability up to 220 km s−1 which is similar to the highest variable velocities on both sides of the other H I lines. Na I lines show variations up to 120 km s−1 on the red side and up to 90 km s−1 on the blue side. The latter is similar to that seen in the red side of the components of the Ca II triplet.

The circumstellar Hα emission, after star model subtraction and radial-velocity correction, has a strongly varying line profile (Fig. 12). The first and last epochs are characterized by a strong and sharp absorption component around ~55 km s−1 on the blue side of the feature. The second epoch shows an inverse P Cygni profile and the third is characterized by two emission peaks, one around the rest wavelength and the other blue shifted. These profiles are discussed in Sect. 5.5.

thumbnail Fig. 11

Variable lines in PDS 111. From left to right: Hα and He I λ 1083.0 nm lines in the spectra of PDS 111 are shown on velocity scale in the upper panel. They show significant variations. This is quantified in the lower panel by the use of the temporal variance spectrum, showing the significant variable velocity regime and similarity in the amplitude of the variation. The second panel shows three Balmer (Hβ, Hγ and Hδ>) lines and Paβ, all displaying variations. The variations are quantified in the lower plot of that panel, showing variations in a similar velocity range as in the lower left panel. The third panel shows the Ca II triplet, and the last panel the Na I doublet. The He I λ 1083.0 nm and Paβ line only include two spectra due to the limited wavelength range of HERMES compared to X-shooter. The HERMES observation on 28-10-2022 shows strong telluric contamination.

4.7 Binaríty

We calculated the contrast curves for the angular differential imaging final imaging. Those have been calculated using the Principal Component Analysis method. Then we assumed an age for the system of 16 Myr and a distance of 158 pc. The contrast curves were converted in mass limits of putative companions around the object. The mass limits are shown in Fig. 13.

We observe no obvious detection of an inner binary component in the SPHERE data. Specifically, we can rule out objects down to ~0.25 M and down to the coronagraphic edge at 0.1”, that is 16 au.

With these limits, we calculated the maximum expected radial velocity shift in the spectrum of the primary for a secondary mass of 0.25 M. At 16 au the radial velocity of the primary would be ~ 1 km s−1, which would increase to ~2 km s−1 when the distance between the primary and secondary is 8 au and to ~5 km s−1 for a separation of 1 au).

The mid-resolution X-shooter spectra and limiting signal-to-noise ratio of the HERMES spectra would be sufficient to detect radial velocity shifts down to at least 5 km s−1 if the epochs would be taken at opportune moments. This would mean that we should detect a secondary mass of 0.25 M if this star would have a separation < 1 au from the primary. The corresponding period of a few months at this separation is shorter than the time between the X-shooter and HERMES observations. Therefore, given our imaging and spectroscopic detection limits, the disk could potentially be a circumbinary disk, hiding a low mass secondary companion of a few tenths of M at a separation between ~ 1−16 au.

thumbnail Fig. 12

Hα line of PDS 111. The upper panel shows the Hα lines of the four epochs and a star model. This star model is used to subtract the stellar contribution to Hα. The result is shown in the lower panel. The blue and gray dotted lines are the intervals used for normalization and measuring the equivalent width.

thumbnail Fig. 13

Detectable minimum mass of co-eval (sub)stellar companions to PDS 111 based on our SPHERE data processed with angular differential imaging.

5 Discussion

5.1 A very flared, but old disk

One of the most intriguing properties of PDS 111 is the amount of flaring of the dust in the disk, which is comparable to systems that are much younger than PDS 111 (Fig. 7). If it is really this old, then such strong flaring suggests that vertical stirring of small grains is still efficient after ~ 16 Myr of evolution, implying that vertical turbulence remains at play. The origin and evolution of turbulence as well as the origin of small grains in protoplanetary disks is a matter of debate. Regarding turbulence, the traditional picture is that the magneto-rotational instabilities (MRI) would be the main source of disk turbulence, but the efficiency of MRI is a complex process that depends on dust properties and abundance, the magnetic field strength and configuration, the gas and dust chemistry, and various ionization sources (see recent review in PPVIJ by Lesur et al. 2022); most of these are not constrained by observations. Other hydro-dynamical instabilities can contribute to vertical mixing, such as the vertical shearing instability (e.g., Flores-Rivera et al. 2020; Barraza-Alfaro et al. 2021). In PDS 111, whatever is dominating the dust vertical mixing still needs to be efficient after ~ 16 Myr of evolution, and be comparable to much younger disks, such as the flared disk of IM Lup (Avenhaus et al. 2018). Franceschi et al. (2023) model this star and suggest that the flaring could be caused by small particles that do not settle, because they do not grow after an initial growth time scale; some remain at high altitude because of long collisional time scales.

An alternative explanation is that the disk around PDS 111 lacks turbulence and that the scattered light observations trace very small fluffy particles, which would not settle in the absence of vertical turbulence (Dullemond & Dominik 2004; Kataoka et al. 2013). Such particles can remain in the surface layers of the disk for several Myr if they do not grow and do not become compact.

Another possible explanation is that the flaring is caused by small dust particles that are added to the disk surface by continuing infall on the disk surface or by blowing dust from the inner system to the disk surface. The possibility of infall from the surroundings to the disk (Dominik & Dullemond 2008) is explored in Sect. 5.4, where we discuss the near environment of PDS 111.

We exclude a disk wind from causing the flared surface layer, as seen in RYTau (Valegård et al. 2022). A wind would make the surface of the disk more smooth without observable substructures in SPHERE, which is not the case in PDS 111. We do observe hints of sporadic outflows (see Sect. 5.5), but we do not observe any forbidden emission lines in the spectrum often related to a disk wind (Nisini et al. 2018).

5.2 Origin of the observed millimeter-cavity

PDS 111 displays a shallow cavity in the dust continuum emission that peaks at 30 au, and which is not detected in scattered light or in the 12CO emission because the emission is hot and optically thick, similar to what has been seen in larger samples of transition disks (e.g., Villenave et al. 2019). This cavity can originate from planet-disk interaction or (magneto-) hydrodynamical instabilities. In the case of planet-disk interaction, the planet mass should be ≲MJup because otherwise a cavity would have been opened in small grains and that would leave a signature in scattered light (de Juan Ovelar et al. 2013), and potentially in the 12CO (Facchini et al. 2018). The mass of such a planet would allow gas to flow through the gap maintaining mass-accretion onto the star (Zhu et al. 2011; Zhang et al. 2018), which is moderate in PDS 111 (Table 3). The detection of one or more planets is challenging with current telescopes (e.g., with JWST, Wagner et al. 2024) because of their potential location (≲30 au) and low mass.

Alternatively, stellar binarity could also explain the presence of a cavity. Its shape of depends on the mass ratio, eccentricity, and disk viscosity (e.g., Artymowicz & Lubow 1994; Miranda et al. 2017). The location of the ring in the binary case can be 3.5–4 times the binary separation. For PDS 111 that would be a binary separation around 8 au. The binary scenario is further discussed in Sect. 5.3.

Other origins of the cavity such as dead zones (e.g., Flock et al. 2015; Delage et al. 2022) cannot be excluded. More observational constraints, such as disk kinematics, may help to understand its origin.

5.3 Binarity

A circumbinary disk might explain several of the observed disk features such as the scattered light asymmetry, as well as the large disk cavity in the mm-continuum, as discussed in Sect. 5.2, for which the predicted separation is around ~8 au. Binarity could also explain the high RUWE value. Juhász & Facchini (2017) performed 3D SPH hydrodynamic simulations of a disk warp induced by a central binary. The synthetic scattered light images show that a disk warp can create strong illumination asymmetries, due to disk shadowing, especially when the warp “axis” is aligned with the disk minor axis. However, we note that their Fig. 5 (panel p) shows that for a continuous warp we would expect the illumination pattern to be reversed between disk top and bottom side. As we discuss in Sect. 4.3, this does not appear to be the case in PDS 111, which may indicate a local increase in scale height rather than a continuous warp.

Recent SPH hydrodynamic models by Ragusa et al. (2020), Muñoz & Lithwick (2020) and Penzlin et al. (2022) suggest that a circumbinary disk that is co-planar with the orbital plane will always develop an eccentricity even if the initial binary orbit was close to circular. The main condition that they find is that the mass ratio between the binary components is q > 0.05. This means in the case of PDS 111 a binary component more massive than 55 MJup would be required. This is not ruled out by our detection limits. Following the relation between cavity size and binary orbit recovered by Ragusa et al. (2020), we would expect the binary component to be located at a semi-major axis of ~8.5 au. This is again consistent with our detection limits. However, we note that the presence of such a massive perturber at this relatively large separation might be in tension with our non-detection of the disk cavity in small grains, that is scattered light.

Ragusa et al. (2020) find that their disks become increasingly eccentric because of resonant excitation or periodic “pumping” via oblique spiral shocks (Ogilvie 2008; Shi et al. 2012). Due to the eccentricity, a “traffic jam” makes the density of the gas and dust higher at the apocenter than at pericenter (Ataiee et al. 2013; Ragusa et al. 2020), potentially explaining the observed asymmetries in PDS 111, as well as the possible local puffed up disk in scattered light.

The ALMA data puts an upper limit on a potential cavity gap size caused by a stellar companion. If there is a gas cavity, we would have detected this in the gas emission if the cavity was larger than the beam size of the 12CO measurement, like in CS Cha (Kurtovic et al. 2022). Therefore, the cavity must be smaller than 0.429″.

Literature only shows a few confirmed cases of prominent young, circumbinary disks seen with SPHERE, for instance HD 142527 (Fukagawa et al. 2006; Biller et al. 2012; Close et al. 2014), HD 34700 (Monnier et al. 2019), GGTau (Roddier et al. 1996; Keppler et al. 2020), and CS Cha (Ginski et al. 2018). In most of these cases, the resolved scattered light disk shows strong azimuthal asymmetries like spiral arms, which are not in our observations of PDS 111. However, this may strongly depend on the binary separation, eccentricity, mass ratio, and orbit inclination relative to the disk. We note that in the case of the CS Cha circumbinary disk the outer scattered light disk appears completely smooth and featureless, highlighting the diversity in disk morphology.

5.4 Environment

PDS 111 is known to be a foreground star of the Orion cloud (Torres et al. 1995), recently also confirmed by its Gaia DR3 distance (ϖ = 6.3136 ± 0.0352 mas; Gaia Collaboration 2022). No strong arguments are found that suggest the star is a member of a young cluster. PDS 111 has been associated with the Columba association (Torres et al. 2008; da Silva et al. 2009; Launhardt et al. 2022). However this predates the availability of an accurate trigonometric parallax measurement (Tycho parallax was ϖ = 28.70 ± 21.60 mas; ESA 1997). Columba is, however, older (424+6$42_{ - 4}^{ + 6}$ Myr; Bell et al. 2015) and an early G-type memberwould already be on the main sequence. The amount of circumstellar matter and our evolutionary age imply PDS 111 is considerably younger. Another membership assignment from the literature is with Theia 160 (Kounkel & Covey 2019), a clumpy string of co-distant and co-moving stars associated with PDS 111. This string shows young stars, but also includes stars that have been attributed to a much older (62 ± 7 Myr) group, the µ Tau association (Gagné et al. 2020). More investigation into the young stars associated with Theia 160 could lead to a PMS track age, which would give additional constraints on the age of the star. So, its genetic connection to any characterized young stellar groups is uncertain at present.

Furthermore, its environment could be important for a rejuvenation scenario, where streamers feed material to the disk (Kuffmeier et al. 2020; Ginski et al. 2021). Arc-like structures, spiral arms, disk misalignment, or other disk structures in later stages of star formation, might be related to late infall of material to the disk. Therefore, late infall could have influenced the flaring and asymmetry in the disk of PDS 111. The STILISM 3D maps (Lallement et al. 2018) predict reddening toward PDS 111 at distance of 158 pc to be very modest E(B–V) = 0.008 ± 0.016 mag (AV = 0.025 ± 0.050 mag), and the reddening is not expected to exceed E(B–V) ~ 0.05 mag until ~475 pc (i.e., near Orion complex), so there does not seem to be any dense, high extinction interstellar dust clouds in its vicinity. Given its disk, any significant reddening or changes in reddening must be due to its own circumstellar disk. Therefore, it seems unlikely that the disk is rejuvenated by its local environment.

5.5 Hα and He I variability

The hot gas from which Ha originates is highly variable in line morphology and shows velocities up to 320 km s−1. These velocities, assuming a Keplerian rotating disk, are formed close to the stellar surface. The inverse PCygni profile, seen in the epoch of 2021-12-07, is widely interpreted as a signature of accretion (e.g., Cauley & Johns-Krull 2014) and has a red-shifted absorption component in the feature around ~100 km s−1. The deep narrow blue-shifted absorption component in the first and last epoch is strongly variable and resides at relatively low velocities (<100 km s−1, Nisini et al. 2018; Bouvier et al. 2023). This He I feature shows a similar blue-shifted absorption feature in the one epoch that also has an X-shooter observation. Since the blue-shifted absorption is deep, narrow, and variable, it probably arises when some cooler optically thick material moves in front of the hot gas from the inner disk where the Hα and He I lines are formed. These blue-shifted features are generally associated with outflows (Kwan et al. 2007; Nisini et al. 2018; Bouvier et al. 2023).

Cauley & Johns-Krull (2014) observed similar profiles in a large sample of He I emission lines from PMS stars and imply an origin of the line morphology from a misalignment between the inner and outer disk. Bouvier et al. (2023) also showed a similar line morphology of the He I and Balmer lines in the classical T Tauri star GMAur with a mass accretion rate of ~10−8 M yr−1 and an inclination of ~68°. By combining optical and NIR photometry, and high-resolution spectroscopic time series, they concluded that these line features are caused by sporadic mass loss from the disk surface. The blue-shifted narrow absorption components in GM Aur start at high velocities, move to lower velocity over a period of days and disappear. Whether PDS 111 shows similar mass loss or an inner-disk misalignment is to be confirmed by follow-up observations with a higher time cadence (e.g., daily).

5.6 Long-lived accretion disks

Recently, Silverberg et al. (2020) found a class of long-lived accretion disks around M-dwarfs with ages above 20 Myr, so-called “Peter Pan disks” for their refusal to “grow up” and disperse. They concluded that these disks were most likely gas-poor hybrid disks, exhibiting signs of young and debris disks, but could not explain why such old systems still harbored these disks. Coleman & Haworth (2020) used disk evolutionary models to study possible conditions for a gas-rich disk to survive around these old stars. They concluded that a low turbulence parameter a ~ 10−4, a high initial disk mass (>0.25 M*), and very low external photo-evaporation (<10−9 M yr−1 ) can reproduce the “Peter Pan disks” around M-dwarfs. Importantly, they also studied how these conditions would change for higher central star masses, as would be the case for the PDS 111 system. They conclude that for a solar mass star, a disk can be maintained up to 50 Myr under similar conditions, if external photo-evaporation is extremely low (<10−10 M yr−1).

5.7 A coherent picture of the disk of PDS 111

We speculate that both radial and azimuthal sub-structures could be caused by an internal perturber. Either a low-mass planet close to the cavity edge at ~30 au, which may explain the continuum cavity as well as the absence of such a cavity in the gas or scattered light, or a close-in stellar binary located between 1.5 and 16 au, with a preference for 8.5 au to explain the mm cavity size. Such a binary may then also explain the azimuthal asymmetry by exciting disk eccentricity. The presence of a perturber in the disk cavity might also go some way to explain the presence of a gas-rich disk at the old age of the system by trapping material in the outer disk. This may be consistent with the low accretion rate of 1−5×10−10 M yr−1 that we inferred from our Hα and Hβ measurements. Still, planet formation and thus cavity opening is presumably a common occurrence in young disks. Related sub-structures are indeed detected in numerous systems with high-resolution ALMA or scattered light observations (e.g., Andrews et al. 2018; Benisty et al. 2023). Yet studies that measure the disk occurrence rate in young clusters typically find that most disks have dispersed after 10 Myr (e.g., Mamajek 2009; Fedele et al. 2010; Ribas et al. 2014).

6 Summary and conclusions

We present the results of an extensive observational campaign of the hitherto relatively unexplored PDS 111 system in the framework of the DESTINYS program. Our spectroscopic, imaging, and mm-interferometry observations explore the system from the stellar photosphere out to the edge of the gaseous disk at ~295 au. They reveal a system that appears archetypal in many ways for a young gas-rich circumstellar disk. Yet our spec-troscopic analysis of the central star as well as the possible membership in the Columba association, place the system age at an age of ~16 Myr (for spectroscopy) or even significantly older (for association membership). Equally puzzling is the relative isolation of PDS 111 in the foreground of the Orion star-forming region.

Our NIR scattered light observations, as well as ALMA continuum observations, resolve the outer disk and reveal radial sub-structures and azimuthal asymmetries. The azimuthal asymmetry might trace a warp or local overdensity in the outer disk. Interestingly, the complex and variable profiles of the Hα line might also indicate a misalignment of the inner disk region.

The “Peter Pan disk” scenario for higher mass central stars fits well with the relative isolation of PDS 111. Thus, it is possible that the combination of low external photo-evaporation due to lack of close association neighbors, as well as the cavity opening and subsequent low viscous transport of material within the disk, conspired in this system to maintain a young-looking disk around an older star. Thus, PDS 111 might well be the first high mass equivalent of the “Peter Pan” population among M-dwarfs. Our extensive data set, and in particular our spatially resolved observations of the disk, make PDS 111 the ideal laboratory to study this class of long-lasting disks with dedicated simulations.

Acknowledgements

SPHERE was designed and built by a consortium made 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 CNRS (France), MPIA (Germany), INAF (Italy), FINES (Switzerland) and NOVA (Netherlands). Additional funding from EC’s 6th and 7th Framework Programmes as part of OPTICON was received (grant number RII3-Ct-2004-001566 for FP6 (2004-2008); 226604 for FP7 (2009-2012); 312430 for FP7 (2013-2016)). We acknowledge the Programme National de Planétologie (PNP) and the Programme National de Physique Stellaire (PNPS) of CNRS-INSU, France, the French Labex OSUG@2020 (Investissements d’avenir – ANR10 LABX56) and LIO (Lyon Institute of Origins, ANR-10-LABX-0066 within the programme Investissements d’Avenir, ANR-11-IDEX-0007), and the Agence Nationale de la Recherche (ANR-14-CE33-0018) for support. A.Z. acknowledges support from ANID - Millennium Science Initiative Program – Center Code NCN2021_080. This paper makes use of the following ALMA data: ADS/JAO.ALMA#2021.1.01705.S. ALMA is a partnership of ESO (representing its member states), NSF (USA) and NINS (Japan), together with NRC (Canada), MOST and ASIAA (Taiwan), and KASI (Republic of Korea), in cooperation with the Republic of Chile. The Joint ALMA Observatory is operated by ESO, AUI/NRAO and NAOJ. Based on observations obtained with the HERMES spec-trograph, which is supported by the Research Foundation - Flanders (FWO), Belgium, the Research Council of KU Leuven, Belgium, the Fonds National de la Recherche Scientifique (F.R.S.-FNRS), Belgium, the Royal Observatory of Belgium, the Observatoire de Genève, Switzerland and the Thüringer Landessternwarte Tautenburg, Germany. This paper includes data collected by the TESS mission. Funding for the TESS mission is provided by the NASA’s Science Mission Directorate. P.P. acknowledge funding from the UK Research and Innovation (UKRI) under the UK government’s Horizon Europe funding guarantee from ERC (grant agreement No 101076489) This project has received funding from the European Research Council (ERC) under the European Union’s Horizon Europe research and innovation program (DUST2PLANETS, grant agreement No. 101053020). This project has received funding from the European Research Council (ERC) under the European Union’s Horizon 2020 research and innovation programme (PROTOPLANETS, grant agreement No. 101002188). Part of this research was carried out at the Jet Propulsion Laboratory, California Institute of Technology, under a contract with the National Aeronautics and Space Administration (80NM0018D0004).

Appendix A Flux calibrated stacked total intensity and polarimetric images

In Fig. A.2 we show the de-rotated and stacked total intensity image, which contains the disk and stellar signal. In Fig. A.1 we show the flux calibrated Stokes Q, U, Qϕ, Uϕ images. Flux calibration was done with the central star as the reference source. For this purpose we used the flux calibration frames which were taken before and after the science images and in which the star is removed from the coronagraph and a neutral density filter inserted to prevent saturation. We measured the central star flux using aperture photometry. We then translate detector counts into units of surface brightness via the known H-band magnitude of the primary star.

thumbnail Fig. A.1

Flux calibrated polarimetric images. The directly measured Stokes Q and U parameters are shown in the top row and the derived azimuthal Stokes Qϕ and Uϕ parameters in the bottom row. All images are shown on the same (logarithmic) scale, symmetric around 0. In all images, the gray, hashed circle in the center marks the area that was covered by a coronagraphic mask.

thumbnail Fig. A.2

Flux calibrated stacked total intensity image. We note that in this image the stellar light was not removed. The bright ring visible on the outer edge of the image shows the adaptive optics correction radius which shows bright stellar speckle signal (smoothed out by the field rotation during observations). The signal of the disk is faintly visible close to the coronagraph. The image is shown on a (logarithmic) scale, symmetric around 0. The grey, hashed circle in the center marks the area that was covered by a coronagraphie mask.

Appendix B Stellar atmosphere model fit

The χ2 values from atmospheric model fitting are shown in Fig. B.1 and a subset of the models, including the best fit, is shown in Fig. 2. The fitting range is from 402 nm to 430 nm based on Gray (2000).

thumbnail Fig. B.1

Resulting reduced χ2 values from atmospheric model fitting. The Kurucz models have different temperatures and surface gravities and a constant υ sini. The lowest χ2 value is indicated with a red dot.

Appendix C SED

In this section we show the B- and V-band from Gaia photometry and the SED discussed in Sect. 4.1. Figs. C.1 and C.2 shows the B- and V-band measurements converted from Gaia measurements (Gaia Collaboration 2022; Jayasinghe et al. 2018). Fig. C.3 shows the SED and Table C.1.

thumbnail Fig. C.1

B- and V-band photometry converted from Gaia measurements. The B- and V-band measurements show variations.

thumbnail Fig. C.2

B-V color and AV shown for V-band measurements determined from Gaia. B–V shows a slight increase with increasing V-band magnitude. The reported AV follows from AV = RV * E(B–V) where E(B–V) is determined by comparing the observed B–V to the intrinsic (B–V)0 of 0.597 mag (interpolated from table 6 in Pecaut & Mamajek (2013) for a star of 5900 K) and also increases with an increasing V-band measurement. The AV values are slightly higher than determined by the SED fit.

thumbnail Fig. C.3

The SED following from Kurucz model fitting to optical photometry, discussed in Sect. 4.1. The model and photometry are dereddend using Cardelli’s extinction law. Infrared excess is observed from the J-band on toward redder wavelengths.

Table C.1

Photometry used for SED fitting.

Appendix D EW measurements

The EW measurements of several lines are reported in Tab. D.1.

Table D.1

EW measurements in Å after stellar model subtraction.

Appendix E Ellipse fitting

Fig. E.1 and E.2 highlight the disk structures and the overlayed ellipses.

thumbnail Fig. E.1

SPHERE scattered light data after high pass filtering (top row) and deconvolution (bottom row). The color map was adjusted in each panel to highlight the most prominent sub-structure in the disk. We overplot ellipses tracing the two inner rings and the following gap from left to right.

thumbnail Fig. E.2

SPHERE Qϕ image shown on a saturated color map to highlight the outer edge of the disk. The red dashed line marks the ellipse also shown in Fig. 7, which traces the outer edge. Note that due to the scattering phase function the disk far side is faint, and thus appears only at low S/N.

Appendix F Components of the visibility fitting

Fig. F.1 shows a sketch of the components of the intensity profile that are assumed for the visibility fitting of the dust continuum emission of the disk around PDS 111 observed with ALMA.

thumbnail Fig. F.1

Schematic of PDS 111 components of the model used to fit the visibilities from the ALMA observations. The model components do not reflect the relative flux of contribution.

Appendix G Channel map of the 12CO

Fig. G.1 shows the channel maps of 12CO of PDS 111 from ALMA observations

thumbnail Fig. G.1

Channel maps of 12CO of PDS 111. The gray contours are at the 5 × σ level of the continuum emission. The scale bar in the left panel represents a scale of 50 au. The dashed line is the PA of the continuum emission in the disk.

References

  1. Alcalá, J. M., Manara, C. F., Natta, A., et al. 2017, A&A, 600, A20 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  2. Andrews, S. M. 2020, ARA&A, 58, 483 [Google Scholar]
  3. Andrews, S. M., Rosenfeld, K. A., Kraus, A. L., & Wilner, D. J. 2013, ApJ, 771, 129 [Google Scholar]
  4. Andrews, S. M., Huang, J., Pérez, L. M., et al. 2018, ApJ, 869, L41 [NASA ADS] [CrossRef] [Google Scholar]
  5. Artymowicz, P., & Lubow, S. H. 1994, ApJ, 421, 651 [Google Scholar]
  6. Ataiee, S., Pinilla, P., Zsom, A., et al. 2013, A&A, 553, L3 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  7. Avenhaus, H., Quanz, S. P., Garufi, A., et al. 2018, ApJ, 863, 44 [NASA ADS] [CrossRef] [Google Scholar]
  8. Bae, J., Isella, A., Zhu, Z., et al. 2023, in ASP Conf. Ser., 534, eds. S. Inutsuka, Y. Aikawa, T. Muto, K. Tomida, & M. Tamura, 423 [NASA ADS] [Google Scholar]
  9. Bailer-Jones, C. A. L., Rybizki, J., Fouesneau, M., Demleitner, M., & Andrae, R. 2021, AJ, 161, 147 [Google Scholar]
  10. Barraza-Alfaro, M., Flock, M., Marino, S., & Pérez, S. 2021, A&A, 653, A113 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  11. Beckwith, S. V. W., Sargent, A. I., Chini, R. S., & Guesten, R. 1990, AJ, 99, 924 [Google Scholar]
  12. Bell, C. P. M., Mamajek, E. E., & Naylor, T. 2015, MNRAS, 454, 593 [Google Scholar]
  13. Benisty, M., Juhász, A., Facchini, S., et al. 2018, A&A, 619, A171 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  14. Benisty, M., Dominik, C., Follette, K., et al. 2023, in ASP Conf. Ser., 534, Protostars and Planets VII, eds. S. Inutsuka, Y. Aikawa, T. Muto, K. Tomida, & M. Tamura, 605 [NASA ADS] [Google Scholar]
  15. Beuzit, J. L., Vigan, A., Mouillet, D., et al. 2019, A&A, 631, A155 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  16. Biller, B., Lacour, S., Juhász, A., et al. 2012, ApJ, 753, L38 [Google Scholar]
  17. Boone, F. 2013, Exp. Astron., 36, 77 [NASA ADS] [CrossRef] [Google Scholar]
  18. Bouvier, J., Sousa, A., Pouilly, K., et al. 2023, A&A, 672, A5 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  19. Briggs, D. S. 1995, in Amer. Astron. Soc. Meeting Abstracts, 187, 112.02 [Google Scholar]
  20. Busso, G., Cacciari, C., Bellazzini, M., et al. 2022, Gaia DR3 documentation Chapter 5: Photometric data, https://gea.esac.esa.int/archive/documentation/GDR3/index.html, 5 [Google Scholar]
  21. Cardelli, J. A., Clayton, G. C., & Mathis, J. S. 1989, ApJ, 345, 245 [Google Scholar]
  22. Carpenter, J. M., Bouwman, J., Mamajek, E. E., et al. 2009, ApJS, 181, 197 [Google Scholar]
  23. Cauley, P. W., & Johns-Krull, C. M. 2014, ApJ, 797, 112 [NASA ADS] [CrossRef] [Google Scholar]
  24. Close, L. M., Follette, K. B., Males, J. R., et al. 2014, ApJ, 781, L30 [NASA ADS] [CrossRef] [Google Scholar]
  25. Coleman, G. A., & Haworth, T. J. 2020, MNRAS, 496, L111 [NASA ADS] [CrossRef] [Google Scholar]
  26. Czekala, I., Loomis, R. A., Teague, R., et al. 2021, ApJS, 257, 2 [NASA ADS] [CrossRef] [Google Scholar]
  27. da Silva, L., Torres, C. A. O., de La Reza, R., et al. 2009, A&A, 508, 833 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  28. Debes, J. H., Poteet, C. A., Jang-Condell, H., et al. 2017, ApJ, 835, 205 [NASA ADS] [CrossRef] [Google Scholar]
  29. de Boer, J., Salter, G., Benisty, M., et al. 2016, A&A, 595, A114 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  30. de Boer, J., Langlois, M., van Holstein, R. G., et al. 2020, A&A, 633, A63 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  31. de Jong, J. A., Henrichs, H. F., Kaper, L., et al. 2001, A&A, 368, 601 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  32. de Juan Ovelar, M., Min, M., Dominik, C., et al. 2013, A&A, 560, A111 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  33. Delage, T. N., Okuzumi, S., Flock, M., Pinilla, P., & Dzyurkevich, N. 2022, A&A, 658, A97 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  34. Derkink, A. R., Ramírez-Tannus, M. C., Kaper, L., et al. 2024, A&A, 681, A112 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  35. Dohlen, K., Langlois, M., Saisse, M., et al. 2008, in SPIE Conf. Ser., 7014, Ground-based and Airborne Instrumentation for Astronomy II, ed. I. S. McLean, & M. M. Casali, 70143L [NASA ADS] [Google Scholar]
  36. Dominik, C., & Dullemond, C. P. 2008, A&A, 491, 663 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  37. Dullemond, C. P., & Dominik, C. 2004, A&A, 421, 1075 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  38. Dullemond, C. P., Birnstiel, T., Huang, J., et al. 2018, ApJ, 869, L46 [NASA ADS] [CrossRef] [Google Scholar]
  39. ESA 1997, VizieR Online Data Catalog: I/239 [Google Scholar]
  40. Facchini, S., Pinilla, P., van Dishoeck, E. F., & de Juan Ovelar, M. 2018, A&A, 612, A104 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  41. Fedele, D., van den Ancker, M. E., Henning, T., Jayawardhana, R., & Oliveira, J. M. 2010, A&A, 510, A72 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  42. Fitton, S., Tofflemire, B. M., & Kraus, A. L. 2022, Res. Notes Am. Astron. Soc., 6, 18 [Google Scholar]
  43. Flock, M., Ruge, J. P., Dzyurkevich, N., et al. 2015, A&A, 574, A68 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  44. Flores-Rivera, L., Flock, M., & Nakatani, R. 2020, A&A, 644, A50 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  45. Franceschi, R., Birnstiel, T., Henning, T., & Sharma, A. 2023, A&A, 671, A125 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  46. Fukagawa, M., Tamura, M., Itoh, Y., et al. 2006, ApJ, 636, L153 [Google Scholar]
  47. Fullerton, A. W., Gies, D. R., & Bolton, C. T. 1996, ApJS, 103, 475 [NASA ADS] [CrossRef] [Google Scholar]
  48. Gagné, J., David, T. J., Mamajek, E. E., et al. 2020, ApJ, 903, 96 [CrossRef] [Google Scholar]
  49. Gaia Collaboration 2022, VizieR Online Data Catalog: I/355 [Google Scholar]
  50. Garufi, A., Avenhaus, H., Pérez, S., et al. 2020, A&A, 633, A82 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  51. Ginski, C., Stolker, T., Pinilla, P., et al. 2016, A&A, 595, A112 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  52. Ginski, C., Benisty, M., van Holstein, R. G., et al. 2018, A&A, 616, A79 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  53. Ginski, C., Ménard, F., Rab, C., et al. 2020, A&A, 642, A119 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  54. Ginski, C., Facchini, S., Huang, J., et al. 2021, ApJ, 908, L25 [NASA ADS] [CrossRef] [Google Scholar]
  55. Gray, R. 2000, A Digital Spectral Classification Atlas (Appalachian State University), https://ned.ipac.caltech.edu/level5/Gray/frames.html [Google Scholar]
  56. Guilloteau, S., Dutrey, A., Piétu, V., & Boehler, Y. 2011, A&A, 529, A105 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  57. Hildebrand, R. H. 1983, QJRAS, 24, 267 [NASA ADS] [Google Scholar]
  58. Hom, E. F. Y., Marchis, F., Lee, T. K., et al. 2007, J. Opt. Soc. Am. A, 24, 1580 [NASA ADS] [CrossRef] [Google Scholar]
  59. Huang, J., Andrews, S. M., Pérez, L. M., et al. 2018, ApJ, 869, L43 [Google Scholar]
  60. Jayasinghe, T., Kochanek, C. S., Stanek, K. Z., et al. 2018, MNRAS, 477, 3145 [Google Scholar]
  61. Jorsater, S., & van Moorsel, G. A. 1995, AJ, 110, 2037 [Google Scholar]
  62. Juhász, A., & Facchini, S. 2017, MNRAS, 466, 4053 [NASA ADS] [Google Scholar]
  63. Karamanis, M., & Beutler, F. 2021, Stat. Comput., 31, 61 [Google Scholar]
  64. Karamanis, M., Beutler, F., & Peacock, J. A. 2021, MNRAS, 508, 3589 [NASA ADS] [CrossRef] [Google Scholar]
  65. Kataoka, A., Tanaka, H., Okuzumi, S., & Wada, K. 2013, A&A, 557, L4 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  66. Kausch, W., Noll, S., Smette, A., et al. 2015, A&A, 576, A78 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  67. Keppler, M., Penzlin, A., Benisty, M., et al. 2020, A&A, 639, A62 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  68. Kiraga, M. 2012, Acta Astron., 62, 67 [NASA ADS] [Google Scholar]
  69. Kounkel, M., & Covey, K. 2019, AJ, 158, 122 [Google Scholar]
  70. Kuffmeier, M., Goicovic, F. G., & Dullemond, C. P. 2020, A&A, 633, A3 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  71. Kuhn, J. R., Potter, D., & Parise, B. 2001, ApJ, 553, L189 [NASA ADS] [CrossRef] [Google Scholar]
  72. Kurtovic, N. T., Pinilla, P., Penzlin, A. B. T., et al. 2022, A&A, 664, A151 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  73. Kurucz, R. L. 1993, VizieR Online Data Catalog: VI/39 [Google Scholar]
  74. Kwan, J., Edwards, S., & Fischer, W. 2007, ApJ, 657, 897 [Google Scholar]
  75. Lallement, R., Capitanio, L., Ruiz-Dern, L., et al. 2018, A&A, 616, A132 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  76. Launhardt, R., Loinard, L., Dzib, S. A., et al. 2022, ApJ, 931, 43 [NASA ADS] [CrossRef] [Google Scholar]
  77. Law, C. J., Teague, R., Oberg, K. I., et al. 2023, ApJ, 948, 60 [NASA ADS] [CrossRef] [Google Scholar]
  78. Lesur, G., Ercolano, B., Flock, M., et al. 2022, arXiv e-prints [arXiv:2203.09821] [Google Scholar]
  79. Lomb, N. R. 1976, Ap&SS, 39, 447 [Google Scholar]
  80. Long, F., Pinilla, P., Herczeg, G. J., et al. 2018, ApJ, 869, 17 [Google Scholar]
  81. Long, F., Andrews, S. M., Rosotti, G., et al. 2022, ApJ, 931, 6 [NASA ADS] [CrossRef] [Google Scholar]
  82. Mamajek, E. E. 2009, in AIP Conf. Ser., 1158, Exoplanets and Disks: Their Formation and Diversity, eds. T. Usuda, M. Tamura, & M. Ishii, 3 [NASA ADS] [CrossRef] [Google Scholar]
  83. Manara, C. F., Beccari, G., Da Rio, N., et al. 2013, A&A, 558, A114 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  84. Manara, C. F., Ansdell, M., Rosotti, G. P., et al. 2023, in ASP Conf. Ser., 534, Protostars and Planets VII, eds. S. Inutsuka, Y. Aikawa, T. Muto, K. Tomida, & M. Tamura, 539 [Google Scholar]
  85. Marois, C., Lafrenière, D., Doyon, R., Macintosh, B., & Nadeau, D. 2006, ApJ, 641, 556 [Google Scholar]
  86. Milli, J., Mouillet, D., Lagrange, A. M., et al. 2012, A&A, 545, A111 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  87. Miranda, R., Muñoz, D. J., & Lai, D. 2017, MNRAS, 466, 1170 [NASA ADS] [CrossRef] [Google Scholar]
  88. Modigliani, A., Goldoni, P., Royer, F., et al. 2010, in SPIE Conf. Ser., 7737, Observatory Operations: Strategies, Processes, and Systems III, 773728 [Google Scholar]
  89. Monnier, J. D., Harries, T. J., Bae, J., et al. 2019, ApJ, 872, 122 [Google Scholar]
  90. Muñoz, D. J., & Lithwick, Y. 2020, ApJ, 905, 106 [CrossRef] [Google Scholar]
  91. Muro-Arena, G. A., Benisty, M., Ginski, C., et al. 2020, A&A, 635, A121 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  92. Nisini, B., Antoniucci, S., Alcalá, J. M., et al. 2018, A&A, 609, A87 [EDP Sciences] [Google Scholar]
  93. Ogilvie, G. I. 2008, MNRAS, 388, 1372 [NASA ADS] [Google Scholar]
  94. Padgett, D., Koerner, D., Wahhaj, Z., et al. 2006, in ASP Conf. Ser., 357, The Spitzer Space Telescope: New Views of the Cosmos, eds. L. Armus, & W. T. Reach, 97 [Google Scholar]
  95. Pascucci, I., Testi, L., Herczeg, G. J., et al. 2016, ApJ, 831, 125 [Google Scholar]
  96. Pecaut, M. J., & Mamajek, E. E. 2013, ApJS, 208, 9 [Google Scholar]
  97. Penzlin, A. B. T., Kley, W., Audiffren, H., & Schäfer, C. M. 2022, A&A, 660, A101 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  98. Pinilla, P., Tazzari, M., Pascucci, I., et al. 2018, ApJ, 859, 32 [Google Scholar]
  99. Pinilla, P., Kurtovic, N. T., Benisty, M., et al. 2021, A&A, 649, A122 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  100. Pyerin, M. A., Delage, T. N., Kurtovic, N. T., et al. 2021, A&A, 656, A150 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  101. Ragusa, E., Alexander, R., Calcino, J., Hirsh, K., & Price, D. J. 2020, MNRAS, 499, 3362 [NASA ADS] [CrossRef] [Google Scholar]
  102. Ramírez-Tannus, M. C., Kaper, L., De Koter, A., et al. 2017, A&A, 604 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  103. Raskin, G., van Winckel, H., Hensberge, H., et al. 2011, A&A, 526, A69 [CrossRef] [EDP Sciences] [Google Scholar]
  104. Ribas, Á., Merín, B., Bouy, H., & Maud, L. T. 2014, A&A, 561, A54 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  105. Ricker, G. R., Winn, J. N., Vanderspek, R., et al. 2015, J. Astron. Telesc. Instrum. Syst., 1, 014003 [Google Scholar]
  106. Roddier, C., Roddier, F., Northcott, M. J., Graves, J. E., & Jim, K. 1996, ApJ, 463, 326 [NASA ADS] [CrossRef] [Google Scholar]
  107. Rojas, G., Gregorio-Hetem, J., & Hetem, A. 2008, MNRAS, 387, 1335 [Google Scholar]
  108. Scargle, J. D. 1982, ApJ, 263, 835 [Google Scholar]
  109. Shi, J.-M., Krolik, J. H., Lubow, S. H., & Hawley, J. F. 2012, ApJ, 749, 118 [Google Scholar]
  110. Sicilia-Aguilar, A., Henning, T., Dullemond, C. P., et al. 2011, ApJ, 742, 39 [NASA ADS] [CrossRef] [Google Scholar]
  111. Siess, L., Dufour, E., & Forestini, M. 2000, A&A, 358, 593 [Google Scholar]
  112. Silverberg, S. M., Wisniewski, J. P., Kuchner, M. J., et al. 2020, ApJ, 890, 106 [NASA ADS] [CrossRef] [Google Scholar]
  113. Skrutskie, M. F., Cutri, R. M., Stiening, R., et al. 2006, AJ, 131, 1163 [NASA ADS] [CrossRef] [Google Scholar]
  114. Smette, A., Sana, H., Noll, S., et al. 2015, A&A, 576, A77 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  115. Stapper, L. M., & Ginski, C. 2022, A&A, 668, A50 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  116. Tazzari, M., Testi, L., Ercolano, B., et al. 2016, A&A, 588, A53 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  117. Tazzari, M., Beaujean, F., & Testi, L. 2018, MNRAS, 476, 4527 [Google Scholar]
  118. Teague, R. 2019, J. Open Source Softw., 4, 1220 [Google Scholar]
  119. Torres, C. A. O., Quast, G., de La Reza, R., Gregorio-Hetem, J., & Lepine, J. R. D. 1995, AJ, 109, 2146 [NASA ADS] [CrossRef] [Google Scholar]
  120. Torres, C. A. O., Quast, G. R., Melo, C. H. F., & Sterzik, M. F. 2008, in Handbook of Star Forming Regions, Volume II, ed. B. Reipurth, 5, 757 [NASA ADS] [Google Scholar]
  121. Trapman, L., Facchini, S., Hogerheijde, M. R., van Dishoeck, E. F., & Bruderer, S. 2019, A&A, 629, A79 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  122. Trapman, L., Ansdell, M., Hogerheijde, M. R., et al. 2020, A&A, 638, A38 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  123. Valegård, P. G., Waters, L. B. F. M., & Dominik, C. 2021, A&A, 652, A133 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  124. Valegård, P. G., Ginski, C., Dominik, C., et al. 2022, A&A, 668, A25 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  125. van Holstein, R. G., Snik, F., Girard, J. H., et al. 2017, in SPIE Conf. Ser., 10400, ed. S. Shaklan, 1040015 [NASA ADS] [Google Scholar]
  126. van Holstein, R. G., Girard, J. H., de Boer, J., et al. 2020, A&A, 633, A64 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  127. Vernet, J., Dekker, H., D’Odorico, S., et al. 2011, A&A, 536, A105 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  128. Villenave, M., Benisty, M., Dent, W. R. F., et al. 2019, A&A, 624, A7 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  129. Wagner, K., Leisenring, J., Cugno, G., et al. 2024, AJ, 167, 181 [NASA ADS] [CrossRef] [Google Scholar]
  130. Wahhaj, Z., Cieza, L., Koerner, D. W., et al. 2010, ApJ, 724, 835 [Google Scholar]
  131. Woitke, P., Min, M., Pinte, C., et al. 2016, A&A, 586, A103 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  132. Wolf, C., Onken, C. A., Luvaul, L. C., et al. 2018, PASA, 35, e010 [Google Scholar]
  133. Wright, E. L., Eisenhardt, P. R. M., Mainzer, A. K., et al. 2010, AJ, 140, 1868 [Google Scholar]
  134. Zacharias, N., Finch, C., Girard, T., et al. 2009, VizieR Online Data Catalog: I/315 [Google Scholar]
  135. Zhang, S., Zhu, Z., Huang, J., et al. 2018, ApJ, 869, L47 [Google Scholar]
  136. Zhu, Z., Nelson, R. P., Hartmann, L., Espaillat, C., & Calvet, N. 2011, ApJ, 729, 47 [NASA ADS] [CrossRef] [Google Scholar]
  137. Zhu, Z., Zhang, S., Jiang, Y.-F., et al. 2019, ApJ, 877, L18 [NASA ADS] [CrossRef] [Google Scholar]

2

The Renormalized Unit Weight Error (RUWE) is 2.6 which is discussed in Sect. 4.1.

3

Tool to extract kinematic information from spatially and spectrally resolved line data, https://eddy.readthedocs.io/en/latest/

4

The circumstellar disk may inflate the RUWE value.

All Tables

Table 1

Stellar parameters of PDS 111.

Table 2

Dust continuum visibility modeling results.

Table 3

Derived mass-accretion rates of Ha and Hβ at each epoch in logṀacc [M yr−1].

Table C.1

Photometry used for SED fitting.

Table D.1

EW measurements in Å after stellar model subtraction.

All Figures

thumbnail Fig. 1

SPHERE/IRDIS H-band (1.65 µm) observations of PDS 111. In all images the star is at the center and the images are aligned with the North-axis up and East-axis to the left. Left: polarized light observation. We show the azimuthal Stokes parameter Qϕ which contains all azimuthally polarized scattered light as a positive signal. Middle: same data set as in the left panel, but reduced with angular differential imaging, showing total intensity. Due to the reduction approach, the disk signal is strongly self-subtracted removing particularly low spatial frequencies. The shape of the detected high spatial frequency sub-structure is likely distorted. For the same reason, the given surface brightness is a lower limit only, which depends in a complex way on the local self subtraction in the image. Right: same image as in the left panel, overlayed with contours (5 σ, 14 σ, 68 σ, 140 σ, 205 σ; relative to the image background noise) from the center panel to highlight the structures that are traced in both images. In all images, the gray, hashed circle in the center marks the area that was covered by the coronagraphic mask and where the star is located. All images are shown on a log scale to highlight the disk morphology. Appendix A shows a flux calibrated version of the polarimetric image.

In the text
thumbnail Fig. 2

Stacked spectrum of PDS 111 shown together with a subset of Kurucz models with different temperatures. The models have a log 𝑔 of 4.2 and a υ sini of 21 km s−1.

In the text
thumbnail Fig. 3

ALMA observations of PDS 111. Left panel: dust continuum emission at 1.25 mm of PDS 111. The resolution of the observations is 0.357” × 0.25” with a position angle of −71.9° and is shown in the lower left part of the figure. Right panel: overlap of ALMA (contours at 10,…,100% of the peak of emission) and SPHERE observations (in color).

In the text
thumbnail Fig. 4

Moment 0, 1, and 8 of the 12CO of PDS 111. Moment 0 (left) is the integrated value of the spectrum, moment 1 (middle) the intensity weighted coordinate to get the velocity dispersion and moment 8 (right) shows the peak value of the spectrum. The beam size is 0.429″ × 0.342″ with a PA of −76°.114, shown as the oval in the lower left corner of the panels.

In the text
thumbnail Fig. 5

Overview of PDS 111 indicating the parts of the system that are traced by the different observational techniques.

In the text
thumbnail Fig. 6

Hertzsprung-Russell diagram with PMS tracks and isochrones from Siess et al. (2000). The location of PDS 111 coincides with an age of 15.93.7+1.7$15.9_{ - 3.7}^{ + 1.7}$ Myr.

In the text
thumbnail Fig. 7

Ellipses fitted to the radial sub-structures traced in scattered light (top) and resulting height of the scattering surface of the disk (bottom). We compare the scattering surface height with the average surface height profile of several T Tauri stars found by Avenhaus et al. (2018) and the extreme flaring case of the Herbig star HD 97048, discussed by Ginski et al. (2016).

In the text
thumbnail Fig. 8

Scattered light disk flipped around the minor axis and divided. The minor axis is indicated by the white, dashed line. The scattered light disk is asymmetric relative to the minor axis. The colorbar gives the flux ratio between the two sides of the disk, showing that the flux deviates up to a factor of two between both sides. Notably, the flux ratio turns around for the region inside of ~0.3″ along the major axis. This may indicate a local change in scale height symmetric around the disk midplane rather than a warp or misalignment since the bottom side displays a similar flux ratio.

In the text
thumbnail Fig. 9

Observations versus model images from our best fit with galario after convolving with the same beam as the observations and the residuals. The contours of the residual map are [−4σ, −2σ, 2σ, 4σ], where negative values are dashed contours and positive values solid contours.

In the text
thumbnail Fig. 10

Azimuthally averaged radial intensity proflies of the deprojected images from ALMA dust continuum and the moment 0 maps of the 12CO from ALMA Band-6 observations of PDS 111. Each profile is normalized to the peak. The ALMA profiles include the standard deviation of each elliptical bin divided by the square root of the number of beams spanning the full azimuthal angle at each radial bin. For the SPHERE observations, we take two radial cuts at the North-East (NE) and South-West (SW) of the major axis of the disk. The vertical lines show the radius that encloses 90% of the emission in each of the maps. The gray region shows the size of the coronagraph in the SPHERE observations.

In the text
thumbnail Fig. 11

Variable lines in PDS 111. From left to right: Hα and He I λ 1083.0 nm lines in the spectra of PDS 111 are shown on velocity scale in the upper panel. They show significant variations. This is quantified in the lower panel by the use of the temporal variance spectrum, showing the significant variable velocity regime and similarity in the amplitude of the variation. The second panel shows three Balmer (Hβ, Hγ and Hδ>) lines and Paβ, all displaying variations. The variations are quantified in the lower plot of that panel, showing variations in a similar velocity range as in the lower left panel. The third panel shows the Ca II triplet, and the last panel the Na I doublet. The He I λ 1083.0 nm and Paβ line only include two spectra due to the limited wavelength range of HERMES compared to X-shooter. The HERMES observation on 28-10-2022 shows strong telluric contamination.

In the text
thumbnail Fig. 12

Hα line of PDS 111. The upper panel shows the Hα lines of the four epochs and a star model. This star model is used to subtract the stellar contribution to Hα. The result is shown in the lower panel. The blue and gray dotted lines are the intervals used for normalization and measuring the equivalent width.

In the text
thumbnail Fig. 13

Detectable minimum mass of co-eval (sub)stellar companions to PDS 111 based on our SPHERE data processed with angular differential imaging.

In the text
thumbnail Fig. A.1

Flux calibrated polarimetric images. The directly measured Stokes Q and U parameters are shown in the top row and the derived azimuthal Stokes Qϕ and Uϕ parameters in the bottom row. All images are shown on the same (logarithmic) scale, symmetric around 0. In all images, the gray, hashed circle in the center marks the area that was covered by a coronagraphic mask.

In the text
thumbnail Fig. A.2

Flux calibrated stacked total intensity image. We note that in this image the stellar light was not removed. The bright ring visible on the outer edge of the image shows the adaptive optics correction radius which shows bright stellar speckle signal (smoothed out by the field rotation during observations). The signal of the disk is faintly visible close to the coronagraph. The image is shown on a (logarithmic) scale, symmetric around 0. The grey, hashed circle in the center marks the area that was covered by a coronagraphie mask.

In the text
thumbnail Fig. B.1

Resulting reduced χ2 values from atmospheric model fitting. The Kurucz models have different temperatures and surface gravities and a constant υ sini. The lowest χ2 value is indicated with a red dot.

In the text
thumbnail Fig. C.1

B- and V-band photometry converted from Gaia measurements. The B- and V-band measurements show variations.

In the text
thumbnail Fig. C.2

B-V color and AV shown for V-band measurements determined from Gaia. B–V shows a slight increase with increasing V-band magnitude. The reported AV follows from AV = RV * E(B–V) where E(B–V) is determined by comparing the observed B–V to the intrinsic (B–V)0 of 0.597 mag (interpolated from table 6 in Pecaut & Mamajek (2013) for a star of 5900 K) and also increases with an increasing V-band measurement. The AV values are slightly higher than determined by the SED fit.

In the text
thumbnail Fig. C.3

The SED following from Kurucz model fitting to optical photometry, discussed in Sect. 4.1. The model and photometry are dereddend using Cardelli’s extinction law. Infrared excess is observed from the J-band on toward redder wavelengths.

In the text
thumbnail Fig. E.1

SPHERE scattered light data after high pass filtering (top row) and deconvolution (bottom row). The color map was adjusted in each panel to highlight the most prominent sub-structure in the disk. We overplot ellipses tracing the two inner rings and the following gap from left to right.

In the text
thumbnail Fig. E.2

SPHERE Qϕ image shown on a saturated color map to highlight the outer edge of the disk. The red dashed line marks the ellipse also shown in Fig. 7, which traces the outer edge. Note that due to the scattering phase function the disk far side is faint, and thus appears only at low S/N.

In the text
thumbnail Fig. F.1

Schematic of PDS 111 components of the model used to fit the visibilities from the ALMA observations. The model components do not reflect the relative flux of contribution.

In the text
thumbnail Fig. G.1

Channel maps of 12CO of PDS 111. The gray contours are at the 5 × σ level of the continuum emission. The scale bar in the left panel represents a scale of 50 au. The dashed line is the PA of the continuum emission in the disk.

In the text

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

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

Initial download of the metrics may take a while.