Free Access
Issue
A&A
Volume 638, June 2020
Article Number L6
Number of page(s) 9
Section Letters to the Editor
DOI https://doi.org/10.1051/0004-6361/202038131
Published online 12 June 2020

© ESO 2018

1. Introduction

As young stellar objects accrete gas as part of their formation (e.g. Hartmann et al. 2016), the heated gas produces emission lines that constitute the main observational probe of this event (e.g. Muzerolle et al. 1998). Hydrogen, which is the main constituent of the accreted gas, produces the most prominent features. While hydrogen emission lines can also result from chromospheric activity (e.g. Manara et al. 2017), sufficiently strong and wide lines are thought to result exclusively from accretion (Jayawardhana et al. 2003; Natta et al. 2004; Rigliaco et al. 2012). In analogy to this process in young stars, young brown dwarfs (BDs) and planetary mass objects have been increasingly examined for accretion signatures (e.g. Bowler et al. 2011; Joergens et al. 2013; Zhou et al. 2014; Manara et al. 2015; Petrus et al. 2020). With the advent of increasingly sophisticated adaptive optics (AO) systems, such studies can now be performed even for BDs and planetary companions at close separations from their parent stars (e.g. Sallum et al. 2015; Cugno et al. 2019). One prominent example is the PDS 70 system, where two accreting planet candidates have been observed both photospherically in the near-infrared (Keppler et al. 2018; Mesa et al. 2019) and through hydrogen emission (Haffert et al. 2019). In both cases there is also additional evidence for the circumplanetary discs from which the accreting gas is thought to originate (Christiaens et al. 2019; Isella et al. 2019).

In 2013, Delorme et al. (2013) published the discovery of a 12−14 MJup circumbinary companion in the system 2MASS J01033563-5515561. Based on the detection of the system architecture, the system is registered in the Washington Double Star catalog (WDS), Mason et al. (2001), as “** DLR 1”. Following the example of common nomenclature for the BD binary “Luhman 16” (Mamajek 2013) based on its WDS identifier, and to avoid the cumbersome identifier in 2MASS and similar catalogues, we therefore refer to the system henceforth as “Delorme 1”. The stellar components are referred to as “A” and “B”, respectively, and the circumstellar companion is referred to as “(AB)b”. In this Letter, we present the discovery of hydrogen and helium emission lines from Delorme 1 (AB)b, using the Integral Field Spectrograph (IFS) of the Multi Unit Spectroscopic Explorer (MUSE; Bacon et al. 2010) in its newly commissioned narrow-field mode (NFM). The lines probably indicate ongoing accretion.

2. System characteristics and observations

Here we present the known characteristics of the Delorme 1 system, and introduce the MUSE NFM alongside the fully resolved system. The observing log and a description of the data reduction process are available in Appendix A.1.

2.1. Characteristics of the Delorme 1 system

The Delorme 1 system was resolved in 2012 by Delorme et al. (2013) in the L′ band, revealing three components. They identified Delorme 1 AB as a ∼0.25″ (12 au) low-mass binary of spectral type (SpT) M5.5, with isochronal masses of 0.19 M and 0.17 M for A and B, respectively, at a distance of 47.2 ± 3.1 pc (Riedel et al. 2014). Delorme 1 (AB)b was discovered at a projected separation of ∼1.77″ (84 au) as a 12−14 MJup substellar companion, assuming a system age of 30 Myr. This young age is based on a high probability of membership in the Tucana-Horologium (THA) association (Gagné et al. 2015). Riedel et al. (2014) noted that the central binary pair is over-luminous, possibly indicating an even younger age for the system (more similar to β Pic), or further multiplicity in both components. The very red JHKS colours of Delorme 1 (AB)b are also compatible with a young planetary-mass object of an early-L SpT, providing an independent indication of youth. Using archival NACO data, Delorme et al. (2013) confirmed it to be a bound companion, and it has subsequently also been imaged in z′ using AstraLux by Janson et al. (2017). Based on the location of Delorme 1(AB)b in a J–K colour magnitude diagram, Liu et al. (2016) found that the companion likely has very low surface gravity (VL-G; another indication of youth). It is very similar in colour to several VL-G field BDs and young companions, for example the 13−14 MJup AB Pic b (SpT L0; Chauvin et al. 2005) and the ∼31MJup CD-35 2722 b (SpT L4; Wahhaj et al. 2011). Finally, the system has been studied by Binks & Jeffries (2017) using WISE data. The authors searched for an infrared excess of the unresolved stellar system that might indicate a disc. They found no clear detection of excess, but ultimately rejected the analysis because of Delorme 1 (AB)b.

2.2. MUSE in narrow-field mode

The observations were carried out with the MUSE IFS, located at UT4 of the European Southern Observatory (ESO) Very Large Telescope (VLT), in the new NFM, using science verification (SV) time (ID: 60.A-9485(A)). A total of 12 data cubes were obtained over 2.5 h. The IFS of MUSE covers a spectral range of 4650−9350 Å, with a resolving power R of 2000−4000 (Bacon et al. 2010). The NFM uses a new laser tomography adaptive optics (LTAO) system to provide near diffraction limited observations over a 7.5″ × 7.5″ field of view (FOV) and spatial sampling of 0.025″ × 0.025″. This increases the spatial resolution by a factor 8 over the wide-field mode. Figure 1 shows the system as resolved by MUSE.

thumbnail Fig. 1.

Composite-colour image of the Delorme 1 system, created from the MUSE NFM observations. The slightly brighter primary component A is located to the south-west in the image, with B to the north-east and the significantly redder companion Delorme 1 (AB)b to the north. The central binary is scaled down in flux by a factor of 200.

3. Analysis and results

We detect very strong Hα and Hβ emission from Delorme 1 (AB)b, accompanied by emission in three He I lines (Fig. 2). This likely indicates that accretion onto this planetary-mass object is ongoing. For key lines (Table 1) we provide line fluxes and line profile widths based on which we estimate the accretion rate using three different models (Table 2). Supplementary figures and tables can be found in Appendix B.

thumbnail Fig. 2.

Full averaged (over 11 cubes) MUSE spectrum of the companion Delorme 1 (AB)b (black line) showing very strong Hα emission, and highlighting other atomic lines (black text and red lines) and molecular absorption features (dashed red lines). The two narrow gaps in the spectrum are caused by sections of unusable data at 5800−6051 Å and 7592−7681 Å.

Table 1.

Individual line characteristics (extinction-corrected flux, EW, and the 10% and 50% line profile widths) of Delorme 1 (AB)b for detected emission lines, obtained from the averaged spectrum.

Table 2.

Accretion luminosities and mass accretion rates for Delorme 1 (AB)b estimated from Table 1, and upper limits inferred from the non-detections of the Ca IRT.

3.1. Detection of very strong line emission in Delorme 1 (AB)b

The detected Hα emission in the spectrum of Delorme 1 (AB)b (Fig. 2) clearly stands out, with an integrated flux of 1.28 ± 0.70 × 10−15 erg cm−2 s−1. Hβ does not have a well-defined local continuum for comparison, and peaks at ∼10% of Hα. Three weak but clearly visible helium emission lines are also present at (in-air wavelengths) 6678.1517 Å, 7065.1771 Å, and 7281.35 Å, in decreasing line strength. We detect no lines in the Ca infrared triplet (IRT), and while these are also possible indicators of accretion (e.g. White & Basri 2003), they are not always observed in accreting objects (e.g. Mohanty et al. 2005). The average Hα equivalent width (EW) of Delorme 1 (AB)b is −135 Å, which is comparable to that of other BDs that are possible accretors (e.g. Jayawardhana et al. 2003; Natta et al. 2004; Herczeg et al. 2009; Joergens et al. 2013). Extinction-corrected line fluxes, line luminosities, EWs, and line profile widths are available in Table 1.

Based on the spectrum of Delorme 1 (AB)b, and using qualitative criteria from Kirkpatrick et al. (1999), we tentatively classify Delorme 1 (AB)b as very late-M or early-L, most likely SpT L0 due to its similarity to the L0 standard BD 2MASP J0345+2540. The Rb I and Cs I doublets are weak but present, and the depth of TiO (8432 Å) is similar to that of CrH (8611 Å) and FeH (8692 Å). Stronger-than-expected absorption depth in the VO molecular bands, particularly at 7800−8000 Å, might indicate a lower-gravity atmosphere (e.g. McGovern et al. 2004), which would be in line with the expectations of youth in this object.

It is clear that the line emission is intrinsic to Delorme 1 (AB)b. We therefore turn our attention to the question of the mechanism behind it.

3.2. Mass accretion rates

For the purposes of the continued analysis of Delorme 1 (AB)b, we make use of the DUSTY00 atmosphere isochronal models for very-low-mass stars and BDs from Chabrier et al. (2000) and Baraffe et al. (2002) in lieu of a more in-depth spectral modelling. Based on the photometry of Delorme 1 (AB)b we obtain an absolute Gunn-i magnitude of 15.7 ± 0.3. Combining this with the JHKSL′ photometry from Delorme et al. (2013) and the z′ from Janson et al. (2017), we find very good agreement with the expected absolute IzJHKSL′ magnitudes for the 0.012 M and 0.015 M isochrones at 10 Myr and 50 Myr, respectively (Table 3). This indicates a probable hot-start mass for Delorme 1 (AB)b in the range of 12−15 MJup, spanning the deuterium-burning limit. This is consistent with Delorme et al. (2013). The use of cold- or warm-start models would imply a higher mass (e.g. Marleau & Cumming 2014), but we restrict ourselves here to hot starts, especially because recent modelling suggests that these are more likely (Berardo et al. 2017; Marleau et al. 2019). Thus, to calculate the mass accretion rate, we adopt a companion mass of MP = 0.012 M with a bolometric luminosity log(Lbol/L) = − 3.58, and take the mean radius of the two models, RP = 0.163R (1.59 RJup).

Table 3.

Absolute magnitudes of Delorme 1 (AB)b compared with the DUSTY00 isochrones (Chabrier et al. 2000; Baraffe et al. 2002) that agree well with the photometry.

We estimated the mass accretion rate acc onto Delorme 1 (AB)b with several different methods. We took two that are in common use for stars and BDs (Natta et al. 2004 and Alcalá et al. 2017), another from Thanathibodee et al. (2019), and we also applied the new gas giant accretion model of Aoyama et al. (2018) and Aoyama & Ikoma (2019).

Line luminosities can be used to estimate acc (e.g. Mohanty et al. 2005; Herczeg & Hillenbrand 2008; Rigliaco et al. 2012; Alcalá et al. 2017). Rigliaco et al. (2012) characterised the empirical correlation between the line luminosity (Lline) and the accretion luminosity (Lacc), and from this, estimated acc. This method was further developed by Alcalá et al. (2017) and expanded to other lines. Using the relations from Alcalá et al. (2017), we calculated the accretion luminosity and mass accretion rates for Hβ, Hα, HeIλ6678, and HeI7065. Table 2 shows that these all provide consistent estimates, with a Hα derived acc of 3.4 × 10−10 ± 0.4MJup yr−1. Furthermore, we estimated upper limits of the mass accretion rate from the non-detection of the Ca IRT using the same relations. We assumed an upper limit of the peak flux at 5σ of the noise level at each line centre. From this we then obtained an integrated flux and an accretion luminosity, and the acc upper limits for the Ca IRT shown in Table 2 are consistent with a non-detection of the Ca IRT in an accretion scenario.

The line profile width can also be used to estimate acc, and a width of >200 km s−1 is commonly used to separate accreting BDs from chromospherically active BDs (e.g. Jayawardhana et al. 2003). Natta et al. (2004) used the 10% width of the Hα line in the following empirical relation:

log ( M ˙ acc ) = 12.89 ( ± 0.3 ) + 9.7 ( ± 0.7 ) × 10 3 H α 10 % . $$ \begin{aligned} \mathrm{log} (\dot{M}_{\rm acc}) = -12.89(\pm 0.3) + 9.7(\pm 0.7) \times 10^{-3} \mathrm{H}\alpha 10\% .\end{aligned} $$(1)

We measure an Hα 10% width for Delorme 1 (AB)b of 105−241 km s−1. The wide range is due to the line being only marginally resolved by MUSE (RHα = 2516, FWHM ∼ 119 km s−1). We discuss this, and the implication for other studies of accreting planets using MUSE, in more detail in Sect. 4. From this width, we tentatively estimate M ˙ acc = 3 × 10 8 1.4 + 0.4 M Jup $ \dot{M}_{\mathrm{acc}} = 3\times10^{-8^{+0.4}_{-1.4}}\,M_{\mathrm{Jup}} $ yr−1 using Eq. (1), which is a factor ∼102 greater than what we obtain from the line luminosity above.

A similar estimate can be obtained from the planet-surface shock model of Aoyama & Ikoma (2019), which combines the line luminosity with the Hα line profile widths to estimate acc, as well as the planetary mass MP (independently of evolutionary models). Using a radius RP = 0.163 R, we tentatively estimate an accretion rate of acc = 1 × 10−8 MJup yr−1, and a planetary mass MP = 11 MJup. If the true line profile widths are closer to the lower limits of the width, this would indicate a similar acc but a substantially lower mass, down to MP ∼ 6 MJup. We note that this lower estimate is discrepant from the estimate provided by the evolutionary models (Table 3). The accretion rates derived from the methods of Natta et al. (2004) and Aoyama & Ikoma (2019) are consistent with the acc = 0.8 × 10−8 MJup yr−1 obtained for Delorme 1 (AB)b using the Hα line luminosity based model of Thanathibodee et al. (2019). It is important to note that the applicability of this method to the case of Delorme 1 (AB)b depends on it having a comparable mass and radius to the PDS 70 planets (Thanathibodee et al. 2019). While both parameters are still highly uncertain for these planets, the mass estimates for Delorme 1 (AB)b thus far are comparable to PDS 70 b, and to the radius and mass assumed for the model by Thanathibodee et al. (2019) (MP = 6 MJup, RP = 1.3 RJup).

4. Discussion

We have presented photometry and spectroscopy from MUSE NFM observations of the Delorme 1 system. From the analysis of the detected emission lines (Fig. 2), we found several indicators that accretion is ongoing in Delorme 1 (AB)b. Using four different methods, we estimated mass accretion rates acc of 0.8−3.0 × 10−8MJup yr−1 (from Natta et al. 2004; Aoyama & Ikoma 2019; Thanathibodee et al. 2019) and 3.4 × 10−10MJup yr−1 (from Alcalá et al. 2017). In Appendix A.2 we also investigate the potential emission contribution by chromospheric activity, which is likely negligible. Based on these results, there are strong indications that Delorme 1 (AB)b is actively accreting. We now aim to discuss this in more detail, and conclude with a brief discussion on the age of the Delorme 1 system. We also determine whether it excludes the hypothesis of ongoing accretion in Delorme 1 (AB)b.

The limited resolution of MUSE affects the precision of acc estimates that make use of line profile shape measurements, such as Eq. (1) (Natta et al. 2004) and in part the model by Aoyama & Ikoma (2019) (which is still constrained by the Hα line luminosity). The Hβ line FWHM is essentially equal to the line-broadening FWHM of MUSE at 4861 Å (177 km s−1), so it is unresolved and cannot be used to distinguish the intrinsic line shape. The Hα line, on the other hand, has a somewhat larger FWHM (133 km s−1) than the line-broadening profile at 6563 Å (119 km s−1), therefore it is marginally but not fully resolved. The reported Hα FWHMs for PDS 70 b and c from Haffert et al. (2019) are smaller than Delorme 1 (AB)b and are effectively unresolved, with widths similar to the FWHM of the line-broadening profile (see also Thanathibodee et al. 2019; Hashimoto et al. 2020), which means that they cannot be meaningfully compared with Delorme 1 (AB)b in this context. While there is clear evidence that the PDS 70 (proto)planets are accreting, this suggests that as MUSE becomes an ever more important instrument in the hunt for accreting planets (Girard et al. 2020), its limitations should be kept in mind as well.

As mentioned in Sect. 3.2, there is a factor ∼102 difference in the mass accretion rates estimated for Delorme 1 (AB)b by the results of Natta et al. (2004), Aoyama & Ikoma (2019) and Thanathibodee et al. (2019) on the one hand, and Alcalá et al. (2017) on the other. This discrepancy could be related to the uncertainty on the line profiles and line ratios, or to that on the origin of the emission (and potential absorption thereof). For example, in the stellar case, hydrogen is excited in the pre-shock gas and is ionised post-shock, while for planetary mass objects, the pre-shock temperatures are too low for the excitation of hydrogen (Aoyama & Ikoma 2019). Following the accretion shock, however, temperatures are high enough to produce hydrogen emission.

In Sect. 2.1 we noted that the age of the Delorme 1 system is still not well constrained. Several works assign Delorme 1 to the THA (e.g. Kraus et al. 2014; Malo et al. 2014; Gagné et al. 2015), while Riedel et al. (2014) noted that the over-luminosity of both Delorme 1 A and B indicates an even younger age. Ujjwal et al. (2020) used data from Gaia DR2 and found a much greater age span for THA, 3−60.5 Myr, with an estimated average age of ∼20 Myr, and several other moving groups. This might add uncertainty to the true age of THA. It is worth noting that an age for Delorme 1 of 30−45 Myr does not automatically exclude it as a possible site of accretion. Murphy et al. (2018) identified the 45 Myr old 0.1 M M5 WISE J0808-6443 as the host of a primordial disc, and these authors were recently joined by Lee et al. (2020), who reported a likely transitional disc around another 0.1 M, M5 star, the 55 Myr old 2MASS J15460752-6258042. Silverberg et al. (2020) identified another four, coining the term “Peter Pan discs” for these types of long-lived accretion discs that continue to be found around very low-mass stars. Perhaps even more interestingly, Boucher et al. (2016) reported the detection of a disc around the L0, ∼14 MJup BD and fellow THA member 2MASS J02265658-5327032. Taken together with our results, we find that an in-depth investigation of the Delorme 1 system is well motivated.

To conclude, the implied masses of the three components indicate a high ratio of companion to host star/system mass of ∼0.03−0.04. This makes it difficult to constrain possible formation scenarios for Delorme 1 (AB)b (Delorme et al. 2013). In light of this, the possibility that the companion Delorme 1 (AB)b is still actively accreting could be an important piece to the puzzle. For example, higher resolution observations of the emission lines could provide more precise accretion estimates and allow for comparison with different classes of very low-mass objects, such as the populations studied in Zhou et al. (2014) and Manara et al. (2015), respectively.

Acknowledgments

We thank the anonymous referee for the very useful comments that improved this Letter. R. Asensio-Torres acknowledges support from the European Research Council under the Horizon 2020 Framework Program via the ERC Advanced Grant Origins 83 24 28. MJ gratefully acknowledges support from the Knut and Alice Wallenberg Foundation. G-DM acknowledges the support of the DFG priority program SPP 1992 “Exploring the Diversity of Extrasolar Planets” (KU 2849/7-1) and the support from the Swiss National Science Foundation under grant BSSGI0_155816 “PlanetsInTime”. Parts of this work have been carried out within the framework of the NCCR PlanetS supported by the Swiss National Science Foundation. This research has made use of the VizieR catalogue access tool, CDS, Strasbourg, France (Ochsenbein et al. 2000), the Washington Double Star Catalog maintained at the US Naval Observatory, and Photutils, an Astropy package for detection and photometry of astronomical sources (Bradley et al. 2019).

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. Allard, F., Hauschildt, P. H., Alexander, D. R., Tamanai, A., & Schweitzer, A. 2001, ApJ, 556, 357 [NASA ADS] [CrossRef] [Google Scholar]
  3. Aoyama, Y., & Ikoma, M. 2019, ApJ, 885, L29 [NASA ADS] [CrossRef] [Google Scholar]
  4. Aoyama, Y., Ikoma, M., & Tanigawa, T. 2018, ApJ, 866, 84 [NASA ADS] [CrossRef] [Google Scholar]
  5. Bacon, R., Accardo, M., & Adjali, L. 2010, Ground-based and Airborne Instrumentation for Astronomy III eds. I. S. McLean, S. K. Ramsay, & H. Takami, Int. Soc. Opt. Photonics (SPIE), 7735, 131 [Google Scholar]
  6. Baraffe, I., Chabrier, G., Allard, F., & Hauschildt, P. H. 2002, A&A, 382, 563 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  7. Berardo, D., Cumming, A., & Marleau, G.-D. 2017, ApJ, 834, 149 [NASA ADS] [CrossRef] [Google Scholar]
  8. Binks, A. S., & Jeffries, R. D. 2017, MNRAS, 469, 579 [NASA ADS] [CrossRef] [Google Scholar]
  9. Boucher, A., Lafrenière, D., Gagné, J., et al. 2016, ApJ, 832, 50 [CrossRef] [Google Scholar]
  10. Bowler, B. P., Liu, M. C., Kraus, A. L., Mann, A. W., & Ireland, M. J. 2011, ApJ, 743, 148 [NASA ADS] [CrossRef] [Google Scholar]
  11. Bradley, L., Sipocz, B., & Robitaille, T. 2019, astropy/photutils: v0.6 [Google Scholar]
  12. Cardelli, J. A., Clayton, G. C., & Mathis, J. S. 1989, ApJ, 345, 245 [NASA ADS] [CrossRef] [Google Scholar]
  13. Chabrier, G., Baraffe, I., Allard, F., & Hauschildt, P. 2000, ApJ, 542, 464 [NASA ADS] [CrossRef] [Google Scholar]
  14. Chauvin, G., Lagrange, A. M., Zuckerman, B., et al. 2005, A&A, 438, L29 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  15. Christiaens, V., Casassus, S., Absil, O., et al. 2019, MNRAS, 486, 5819 [NASA ADS] [CrossRef] [Google Scholar]
  16. Cugno, G., Quanz, S. P., Hunziker, S., et al. 2019, A&A, 622, A156 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  17. Delorme, P., Gagné, J., Girard, J. H., et al. 2013, A&A, 553, L5 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  18. DENIS Consortium 2005, VizieR Online Data Catalog: B/denis [Google Scholar]
  19. Gagné, J., Lafrenière, D., Doyon, R., Malo, L., & Artigau, É. 2015, ApJ, 798, 73 [NASA ADS] [CrossRef] [Google Scholar]
  20. Girard, J. H., de Boer, J., Haffert, S., et al. 2020, O4ELT6 Conference Proceedings (2019) [arXiv:2003.02145] [Google Scholar]
  21. Haffert, S. Y., Bohn, A. J., de Boer, J., et al. 2019, Nat. Astron., 3, 749 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  22. Hartmann, L., Herczeg, G., & Calvet, N. 2016, ARA&A, 54, 135 [NASA ADS] [CrossRef] [Google Scholar]
  23. Hashimoto, J., Aoyama, Y., Konishi, M., et al. 2020, AJ, 159, 222 [NASA ADS] [CrossRef] [Google Scholar]
  24. Herczeg, G. J., & Hillenbrand, L. A. 2008, ApJ, 681, 594 [NASA ADS] [CrossRef] [Google Scholar]
  25. Herczeg, G. J., Cruz, K. L., & Hillenbrand, L. A. 2009, ApJ, 696, 1589 [Google Scholar]
  26. Isella, A., Benisty, M., Teague, R., et al. 2019, ApJ, 879, L25 [NASA ADS] [CrossRef] [Google Scholar]
  27. Janson, M., Durkan, S., Hippler, S., et al. 2017, A&A, 599, A70 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  28. Jayawardhana, R., Mohanty, S., & Basri, G. 2003, ApJ, 592, 282 [NASA ADS] [CrossRef] [Google Scholar]
  29. Joergens, V., Bonnefoy, M., Liu, Y., et al. 2013, A&A, 558, L7 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  30. Keppler, M., Benisty, M., Müller, A., et al. 2018, A&A, 617, A44 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  31. Kirkpatrick, J. D., Reid, I. N., Liebert, J., et al. 1999, ApJ, 519, 802 [NASA ADS] [CrossRef] [Google Scholar]
  32. Kraus, A. L., Shkolnik, E. L., Allers, K. N., & Liu, M. C. 2014, AJ, 147, 146 [NASA ADS] [CrossRef] [Google Scholar]
  33. Lee, J., Song, I., & Murphy, S. 2020, MNRAS, 494, 62 [CrossRef] [Google Scholar]
  34. Liu, M. C., Dupuy, T. J., & Allers, K. N. 2016, ApJ, 833, 96 [NASA ADS] [CrossRef] [Google Scholar]
  35. Malo, L., Artigau, É., Doyon, R., et al. 2014, ApJ, 788, 81 [NASA ADS] [CrossRef] [Google Scholar]
  36. Mamajek, E. E. 2013, ArXiv e-prints [arXiv:1303.5345] [Google Scholar]
  37. Manara, C. F., Testi, L., Rigliaco, E., et al. 2013, A&A, 551, A107 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  38. Manara, C. F., Testi, L., Natta, A., & Alcalá, J. M. 2015, A&A, 579, A66 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  39. Manara, C. F., Frasca, A., Alcalá, J. M., et al. 2017, A&A, 605, A86 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  40. Marleau, G.-D., & Cumming, A. 2014, MNRAS, 437, 1378 [NASA ADS] [CrossRef] [Google Scholar]
  41. Marleau, G.-D., Mordasini, C., & Kuiper, R. 2019, ApJ, 881, 144 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  42. Mason, B. D., Wycoff, G. L., Hartkopf, W. I., Douglass, G. G., & Worley, C. E. 2001, AJ, 122, 3466 [Google Scholar]
  43. McGovern, M. R., Kirkpatrick, J. D., McLean, I. S., et al. 2004, ApJ, 600, 1020 [NASA ADS] [CrossRef] [Google Scholar]
  44. Mesa, D., Keppler, M., Cantalloube, F., et al. 2019, A&A, 632, A25 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  45. Mohanty, S., Basri, G., & Jayawardhana, R. 2005, Astron. Nachr., 326, 891 [NASA ADS] [CrossRef] [Google Scholar]
  46. Murphy, S. J., Mamajek, E. E., & Bell, C. P. M. 2018, MNRAS, 476, 3290 [NASA ADS] [CrossRef] [Google Scholar]
  47. Musso Barcucci, A., Cugno, G., Launhardt, R., et al. 2019, A&A, 631, A84 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  48. Muzerolle, J., Calvet, N., & Hartmann, L. 1998, ApJ, 492, 743 [NASA ADS] [CrossRef] [Google Scholar]
  49. Natta, A., Testi, L., Muzerolle, J., et al. 2004, A&A, 424, 603 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  50. Ochsenbein, F., Bauer, P., & Marcout, J. 2000, A&AS, 143, 23 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  51. Oke, J. B., & Gunn, J. E. 1983, ApJ, 266, 713 [NASA ADS] [CrossRef] [Google Scholar]
  52. Petrus, S., Bonnefoy, M., Chauvin, G., et al. 2020, A&A, 633, A124 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  53. Riedel, A. R., Finch, C. T., Henry, T. J., et al. 2014, AJ, 147, 85 [NASA ADS] [CrossRef] [Google Scholar]
  54. Rigliaco, E., Natta, A., Testi, L., et al. 2012, A&A, 548, A56 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  55. Sallum, S., Follette, K. B., Eisner, J. A., et al. 2015, Nature, 527, 342 [NASA ADS] [CrossRef] [Google Scholar]
  56. Schlafly, E. F., & Finkbeiner, D. P. 2011, ApJ, 737, 103 [NASA ADS] [CrossRef] [Google Scholar]
  57. Silverberg, S. M., Wisniewski, J. P., Kuchner, M. J., et al. 2020, ApJ, 890, 106 [CrossRef] [Google Scholar]
  58. Soderblom, D. R. 2010, ARA&A, 48, 581 [NASA ADS] [CrossRef] [Google Scholar]
  59. Thanathibodee, T., Calvet, N., Bae, J., Muzerolle, J., & Hernández, R. F. 2019, ApJ, 885, 94 [NASA ADS] [CrossRef] [Google Scholar]
  60. Ujjwal, K., Kartha, S. S., Mathew, B., Manoj, P., & Narang, M. 2020, AJ, 159, 166 [CrossRef] [Google Scholar]
  61. Wahhaj, Z., Liu, M. C., Biller, B. A., et al. 2011, ApJ, 729, 139 [NASA ADS] [CrossRef] [Google Scholar]
  62. Weilbacher, P. M., Streicher, O., Urrutia, T., et al. 2012, in Software and Cyberinfrastructure for Astronomy II, eds. N. M. Radziwill, & G. Chiozzi, Int. Soc. Opt. Photonics (SPIE), 8451, 102 [Google Scholar]
  63. White, R. J., & Basri, G. 2003, ApJ, 582, 1109 [NASA ADS] [CrossRef] [Google Scholar]
  64. Zhou, Y., Herczeg, G. J., Kraus, A. L., Metchev, S., & Cruz, K. L. 2014, ApJ, 783, L17 [NASA ADS] [CrossRef] [Google Scholar]
  65. Zuckerman, B., & Song, I. 2004, ARA&A, 42, 685 [NASA ADS] [CrossRef] [MathSciNet] [Google Scholar]

Appendix A: MUSE data reduction and chromospheric activity

A.1. MUSE data reduction

Table A.1 lists the 12 integrations obtained over 2.5 h during overall favourable conditions. These were executed using the NFM generic offset observing template, which includes a small-step dither of 0.05″ (2 pixels) and a field rotation of 90° between each exposure in a four-point pattern.

Table A.1.

Observing log of the 2018 September 18 observation of Delorme 1 under ESO SV programme ID 60.A-9485(A).

We used the MUSE data reduction pipeline (Weilbacher et al. 2012) to process the observations. Each observation resulted in a data cube consisting of 3861 individual frames. The wide angular separation between the central binary and the companion, combined with the high-performance AO, results in the companion being easily resolved, as shown in Fig. 1. It is located well outside the halo of the stellar point spread function (PSF), which means that stellar subtraction techniques are not required.

Aperture photometry using the Python routine photutils was then used to extract the spectrum for Delorme 1 (AB)b. Circular aperture radii in the range of 2.0−6.0 pixels were used in increments of 0.5 pixels. Background was subtracted by placing two circular apertures of the same radius on either side of the target, equidistant from the photometric centre, providing an average value of the background at the target position. An aperture radius of 4.0 pixels was found to provide a good balance between a high signal-to-noise ratio (S/N) and good coverage of the companion PSF. A similar technique was used for Delorme 1 A and B, with an aperture radius of 2.5 pixels.

Flux was calibrated automatically in the pipeline, and was later verified by integrating the extinction-corrected (AV = 0.0582; Schlafly & Finkbeiner 2011) flux from the central pair over a suitable band pass. For individual line fluxes, the line-specific extinction correction was calculated from Cardelli et al. (1989), and this value was used for AV. Delorme 1 was observed in the Gunn-i (Oke & Gunn 1983) band during the DEep Near-Infrared Survey (DENIS; DENIS Consortium 2005) at 11.38 mag. Using a wide aperture for the central pair and a four-pixel aperture for the companion, we obtain apparent Gunn-i magnitudes of 11.31 ± 0.06 mag for Delorme 1 AB and 19.08 ± 0.16 mag for Delorme 1 (AB)b. This corresponds to an absolute magnitude of 15.7 ± 0.3, assuming a distance of 47.2 ± 3.1 pc. This brightness is comparable to BDs of similar SpT and distance as Delorme 1 (AB)b.

A.2. Chromospheric activity

An alternative to the hypothesis that Delorme 1 (AB)b is actively accreting is that the emission is instead due to chromospheric activity. Hα emission thought to be originating from chromospheric activity has been observed in a number of BDs. Recent examples include the more massive BD companion PZ TEL B (Musso Barcucci et al. 2019) and a survey of Sco-Cen BDs by Petrus et al. (2020), but typically with lower levels of emission than we find for Delorme 1 (AB)b. Manara et al. (2013) obtained an analytical relation from observations of a number of young (< 10 Myr) K7 – M9.5 stars and BDs for the purpose of estimating the chromospheric activity contribution as an accretion-luminosity noise. This work was extended by Manara et al. (2017) by including more stars in the lowest mass-range, resulting in the relation

log ( L acc , noise / L bol ) = ( 6.2 ± 0.5 ) · log T eff ( 24.5 ± 1.9 ) . $$ \begin{aligned} \mathrm{log}(L_{\mathrm{acc,noise}}/L_{\mathrm{bol}}) = (6.2 \pm 0.5) \cdot \mathrm{log}T_{\mathrm{eff}} - (24.5\pm 1.9) .\end{aligned} $$(A.1)

While this relation includes substantial uncertainties and Delorme 1 (AB)b may fall outside it because it is a very low-mass object, it is worth investigating the potential contribution from chromospheric activity. The 10 Myr, 0.012 M isochrone (Sect. 3.2) yields an effective temperature Teff = 1801 K and log(Lbol/L) = −3.58, resulting in log(Lacc, noise/L) = − 8.0 ± 2.6. Based on the relation from Alcalá et al. (2017), this implies a mass accretion rate noise acc,noise = 5.2×10−12±2.6 MJup yr−1, which is a factor 10−2–10−4 lower than the previously estimated accretion rates in this work. Activity is generally expected to decrease with age (e.g. Zuckerman & Song 2004; Soderblom 2010; Malo et al. 2014). Consequently, from the point of view of activity, it should be even less likely to measure such high Hα emission in Delorme 1 (AB)b at ages later than 10 Myr.

Appendix B: Supplementary material

This appendix includes a number of supplementary figures and tables. Figure B.1 shows the spectra of the central binary pair, with both components showing signs of strong hydrogen emission. The nature of this emission will be explored further in Eriksson et al. (in prep.). By comparing the spectrum of Delorme 1 (AB)b with the spectra of the central pair (Fig. B.1), we can immediately note a number of differences that indicate that detected spectral features and line emission are intrinsic to each object. (i) The primary component, A, is undergoing a flare event at the time of observation. This is evident from comparing the temporal evolution of the EW during the observation (Table B.1 and Fig. B.2), which for Delorme 1 (AB)b shows an evolution that does not correlate with either A or B. (ii) He I emission can only be seen in the companion, and not in the primary or secondary component. (iii) The Hα/Hβ line flux ratios of A and B are ∼1/3 of (AB)b, and the 10% and 50% line profile widths of Hα and Hβ are narrower than the stellar counterparts.

thumbnail Fig. B.1.

Upper panel: spectrum of Delorme 1 (AB)b as shown in Fig. 2. Bottom left: full averaged (over 12 cubes) MUSE IFS spectrum of the two components of the central binary pair, Delorme 1 AB, indicating nearly identical spectral types (∼M5 − M6). Bottom right: secondary component (B) shifted down in flux by 10−14 for clarity. Prominent features include very strong Hα and Na I doublet, and weaker Hβ lines in both components. The primary component (A) shows signs of Ca IRT emission, with the first two Ca II lines visible at 8498/8542 Å. This is likely associated with the recent flare because it is not detectable in B. We also note the lack of any He I features in the 6678−7281 Å range, in contrast to the companion spectrum. We detect a blueshift of ∼ − 15 km s−1 in the Hα line of Delorme 1 A during the flare event, which is not detected in B. We detect no significant blue- or redshift in the lines of Delorme 1 (AB)b.

thumbnail Fig. B.2.

Hα EWs during the observation (from Table B.1) of the central binary components Delorme 1 A and Delorme 1 B (top), and companion Delorme 1 (AB)b (bottom). Negative values indicate emission. From the evolution of the EW over 2.5 h, it is clear that A is undergoing a rapid decrease in width (weaker emission) until it plateaus at around −40 Å, while B stays at an average of −32 Å, and (AB)b fluctuates around an average of −135 Å. There is no clear correlation between the changes in EW in the central binary pair and the companion, indicating that emission is intrinsic to each object. We also note that the rapid change in the EW of Delorme 1 A, as the excess emission due to the flare dissipates, is absent in B and (AB)b.

Table B.1.

Hα EWs of the three components, A, B, and (AB)b, in the Delorme 1 system during the observation.

Table B.2 lists the line characteristics for Delorme 1 (AB)b (Hβ and Hα) for the individual datacubes, with accompanying mass accretion rates listed in Table B.3. Figure B.3 shows the resolution curve for MUSE over its wavelength coverage. The interpolated R values for the lines are discussed in Sect. 3, along with the corresponding line profile Δv. Finally, Fig. B.4 illustrates to which extent the Hβ and Hα line profiles are solved by MUSE in Delorme 1 (AB)b, and Fig. B.5 shows that the line profiles are better resolved for Delorme 1 A and B than for Delorme 1 (AB)b. To estimate how well Hα is resolved, we convolved the MUSE line-broadening profile with continuously narrower widths and locate the smallest width that would produce a convolution matching the observed profile. To illustrate how the uncertainty is reduced at greater widths, we repeated this estimation for Delorme 1 A and find a 10% width of 217−307 km s−1 (FWHM 169 km s−1).

thumbnail Fig. B.3.

Interpolated values for R at the four emission lines used in the data analysis (coloured points). The blue line is obtained from the values listed in 500 Å increments, in the MUSE User Manual (version 11.3, Fig. 18). The Δv values given for each point represent the FWHM of the MUSE line profile at that resolution.

thumbnail Fig. B.4.

Comparison between observed spectral line profiles with errors (solid lines) and the line-broadening profile of MUSE at each line except for He Iλ 7065 (horizontal bars at FWHM and dashed lines). The normalised (pseudo-)continuum level for He Iλ 6678 is fit by eye to fmin = 0.5. Hα is marginally resolved (wider than the dashed line-broadening profile), while Hβ is unresolved, and we therefore consider the measured Hβ line profile widths to be upper limits.

thumbnail Fig. B.5.

Similar to Fig. B.4, but for the central pair Delorme 1 A and Delorme 1 B. The fmin values denote the continuum level. Compared with Delorme 1 (AB)b, the Hα profiles are much better resolved here, while Hβ is marginally resolved. We estimate the Delorme 1 A 10% widths of Hβ and Hα to be 195−353 km s−1 and 217−307 km s−1, respectively.

Table B.2.

Apparent and extinction-corrected (AV = 0.0582, with Aline calculated from Cardelli et al. 1989) Hβ and Hα emission line fluxes, line luminosities, accretion luminosities, and the 10% Hα width range during the observation on 2018 September 18.

Table B.3.

Resulting mass accretion rates from the accretion luminosities and 10% Hα widths in Table B.2.

All Tables

Table 1.

Individual line characteristics (extinction-corrected flux, EW, and the 10% and 50% line profile widths) of Delorme 1 (AB)b for detected emission lines, obtained from the averaged spectrum.

Table 2.

Accretion luminosities and mass accretion rates for Delorme 1 (AB)b estimated from Table 1, and upper limits inferred from the non-detections of the Ca IRT.

Table 3.

Absolute magnitudes of Delorme 1 (AB)b compared with the DUSTY00 isochrones (Chabrier et al. 2000; Baraffe et al. 2002) that agree well with the photometry.

Table A.1.

Observing log of the 2018 September 18 observation of Delorme 1 under ESO SV programme ID 60.A-9485(A).

Table B.1.

Hα EWs of the three components, A, B, and (AB)b, in the Delorme 1 system during the observation.

Table B.2.

Apparent and extinction-corrected (AV = 0.0582, with Aline calculated from Cardelli et al. 1989) Hβ and Hα emission line fluxes, line luminosities, accretion luminosities, and the 10% Hα width range during the observation on 2018 September 18.

Table B.3.

Resulting mass accretion rates from the accretion luminosities and 10% Hα widths in Table B.2.

All Figures

thumbnail Fig. 1.

Composite-colour image of the Delorme 1 system, created from the MUSE NFM observations. The slightly brighter primary component A is located to the south-west in the image, with B to the north-east and the significantly redder companion Delorme 1 (AB)b to the north. The central binary is scaled down in flux by a factor of 200.

In the text
thumbnail Fig. 2.

Full averaged (over 11 cubes) MUSE spectrum of the companion Delorme 1 (AB)b (black line) showing very strong Hα emission, and highlighting other atomic lines (black text and red lines) and molecular absorption features (dashed red lines). The two narrow gaps in the spectrum are caused by sections of unusable data at 5800−6051 Å and 7592−7681 Å.

In the text
thumbnail Fig. B.1.

Upper panel: spectrum of Delorme 1 (AB)b as shown in Fig. 2. Bottom left: full averaged (over 12 cubes) MUSE IFS spectrum of the two components of the central binary pair, Delorme 1 AB, indicating nearly identical spectral types (∼M5 − M6). Bottom right: secondary component (B) shifted down in flux by 10−14 for clarity. Prominent features include very strong Hα and Na I doublet, and weaker Hβ lines in both components. The primary component (A) shows signs of Ca IRT emission, with the first two Ca II lines visible at 8498/8542 Å. This is likely associated with the recent flare because it is not detectable in B. We also note the lack of any He I features in the 6678−7281 Å range, in contrast to the companion spectrum. We detect a blueshift of ∼ − 15 km s−1 in the Hα line of Delorme 1 A during the flare event, which is not detected in B. We detect no significant blue- or redshift in the lines of Delorme 1 (AB)b.

In the text
thumbnail Fig. B.2.

Hα EWs during the observation (from Table B.1) of the central binary components Delorme 1 A and Delorme 1 B (top), and companion Delorme 1 (AB)b (bottom). Negative values indicate emission. From the evolution of the EW over 2.5 h, it is clear that A is undergoing a rapid decrease in width (weaker emission) until it plateaus at around −40 Å, while B stays at an average of −32 Å, and (AB)b fluctuates around an average of −135 Å. There is no clear correlation between the changes in EW in the central binary pair and the companion, indicating that emission is intrinsic to each object. We also note that the rapid change in the EW of Delorme 1 A, as the excess emission due to the flare dissipates, is absent in B and (AB)b.

In the text
thumbnail Fig. B.3.

Interpolated values for R at the four emission lines used in the data analysis (coloured points). The blue line is obtained from the values listed in 500 Å increments, in the MUSE User Manual (version 11.3, Fig. 18). The Δv values given for each point represent the FWHM of the MUSE line profile at that resolution.

In the text
thumbnail Fig. B.4.

Comparison between observed spectral line profiles with errors (solid lines) and the line-broadening profile of MUSE at each line except for He Iλ 7065 (horizontal bars at FWHM and dashed lines). The normalised (pseudo-)continuum level for He Iλ 6678 is fit by eye to fmin = 0.5. Hα is marginally resolved (wider than the dashed line-broadening profile), while Hβ is unresolved, and we therefore consider the measured Hβ line profile widths to be upper limits.

In the text
thumbnail Fig. B.5.

Similar to Fig. B.4, but for the central pair Delorme 1 A and Delorme 1 B. The fmin values denote the continuum level. Compared with Delorme 1 (AB)b, the Hα profiles are much better resolved here, while Hβ is marginally resolved. We estimate the Delorme 1 A 10% widths of Hβ and Hα to be 195−353 km s−1 and 217−307 km s−1, respectively.

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.