Open Access
Issue
A&A
Volume 691, November 2024
Article Number A169
Number of page(s) 10
Section Interstellar and circumstellar matter
DOI https://doi.org/10.1051/0004-6361/202452120
Published online 08 November 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

The timescale for planet formation corresponds to the first few million years of stars’ lifetimes, while they still host a ‘protoplanetary disc’ of dust and gas (e.g. Haisch et al. 2001). This timescale is comparable to the typical timescale over which stars form from overdensities in the galactic interstellar medium (ISM, e.g. Jeffreson et al. 2024). As a result, planet formation in the protoplanetary disc occurs simultaneously with processes governing star formation, in regions of enhanced stellar and ISM density. Theoretical and observational evidence indicates that several processes, including dynamical encounters (Cuello et al. 2023) and external photoevaporation (Winter & Haworth 2022) can feedback on the planet formation process. As a result, the external star formation environment may contribute to the diversity of exoplanets. However, how exactly the environment may sculpt planetary systems remains poorly understood.

One process that has been suggested as an important driver of disc evolution is late stage infall of gas from the ISM via the Bondi-Hoyle-Lyttleton (BHL) accretion mechanism, which would explain the observed scaling between accretion rates and stellar mass (Padoan et al. 2005; Throop & Bally 2008; Klessen & Hennebelle 2010). However, a lack of apparent correlation between stellar accretion rates and the environment (based on substantial accretion rates in Trumpler 37), as well as some theoretical problems in redistributing material to the inner disc, led to the widely held conclusion that BHL accretion is not a primary driver of stellar accretion (Hartmann et al. 2006, see also discussion in Section 4.1 of Winter et al. 2024). Recent years have resulted in a wealth of new observational evidence that motivates a second look at this conclusion. In particular, a recent surge in discoveries of examples of extended structures in molecular line emission (Huang et al. 2021) and infrared scattered light around class II discs (Benisty et al. 2023; Garufi et al. 2024) have provoked renewed interest in late stage infall. Theoretical studies (Dullemond et al. 2019; Kuffmeier et al. 2020, 2023; Hanawa et al. 2024) have suggest infall is a viable driver of observed structures such as spirals and streamers (Huang et al. 2023; Boccaletti et al. 2020), which are often associated with reflection nebulosity (Gupta et al. 2023), accretion outbursts (Pérez et al. 2020; Hales et al. 2024; Dong et al. 2022) and misaligned inner discs (Ginski et al. 2021). The latter misalignments are apparent among a high fraction of discs (Villenave et al. 2024). Ongoing infall may also help explain an apparent mass budget problem for protoplanetary discs compared with the observed exoplanet population (Manara et al. 2018), and the non-monotonic evolution of disc dust masses with time (Testi et al. 2016, 2022; Cieza et al. 2019; Cazzoletti et al. 2019; Williams et al. 2019). Infall was also suggested as a possible solution for the presence of old but still accreting discs in star-forming regions (Scicluna et al. 2014), and even a partial explanation for chemical inhomogeneities in globular clusters (Bastian et al. 2013; Winter & Clarke 2023). From a theoretical perspective, recent simulations have suggested that BHL accretion is a substantial source of mass and angular momentum for discs up to ~1 Myr after the formation of the host star (Kuffmeier et al. 2023; Pelkonen et al. 2024; Padoan et al. 2024). Warps instigated by late-stage infall may propagate throughout the disc (Kimmig & Dullemond 2024; Tanious et al. 2024), possibly driving instabilities and short-lived high levels of turbulence in the outer disc (Deng et al. 2020; Fairbairn & Ogilvie 2023). Shadows resulting from misalignment of the inner disc may also drive spirals and turbulence in the outer disc (Montesinos et al. 2016; Benisty et al. 2017; Su & Bai 2024; Zhang & Zhu 2024). Even if inner disc material loses angular momentum by some other mechanism, such as magnetic winds (e.g. Bai & Stone 2013; Tabone et al. 2022), infall appears a plausible mechanism to replenish this gas, possibly regulating stellar accretion. This replenishment seems particularly plausible given that the near infrared excess, probing inner disc material, is correlated with outer disc substructures such as spirals and shadows (Garufi et al. 2018), suggesting that processes that influence one also influence the other.

The question of the timescale over which late-stage accretion occurs is of critical importance. If the infall onto the disc is only substantial in the first <1 Myr of disc evolution, then it may be reasonable to assume that this simply sets initial conditions, as conventionally assumed in theoretical studies of planet formation. On the other hand, if substantial infall occurs over the entire disc lifetime, disc replenishment may be a critical ingredient for planet formation models. Observationally, examples of infall have not only been uncovered for very young stars that are embedded in their natal environment (Jørgensen et al. 2009; Maury et al. 2019; Pineda et al. 2020; Codella et al. 2024), but also stars older than 1 Myr such as AB Aur (Nakajima & Golimowski 1995; Grady et al. 1999; Fukagawa et al. 2004), DR Tau (Mesa et al. 2022), S CrA (Gupta et al. 2024) and SU Aur (Ginski et al. 2021). Garufi et al. (2024) find evidence of interaction between the star-disc system and the ISM around 16% of their sample of discs in Taurus observed in scattered light. From a theoretical perspective, Winter et al. (2024) recently estimated that ~20–70% of discs in the age range 1-3 Myr should be mostly composed of recently accreted material, based on the timescale for turbulent fluctuations of the ISM (although this finding may be contingent on stellar feedback). Given the potential critical importance to planet formation models, these studies motivate urgent empirical tests as to whether protoplanetary discs are undergoing significant environmental replenishment.

One approach for testing the importance of replenishment is to correlate observed stellar accretion rates with the external environment. In a medium where the gas relative velocity ∆υgascs, the sound speed, the BHL accretion rate BhL is proportional to the cross-section carved out by the radius RBHL=2Gm/Δvgas2${R_{{\rm{BHL}}}}{\rm{ = }}2G{m_ * }/{\rm{\Delta }}\v _{{\rm{gas}}}^{\rm{2}}$ within which the gas is captured, where m*, is the stellar mass. Then M˙BHLm2ρgas/Δvgas3${\dot M_{{\rm{BHL}}}} \propto m_ * ^2{\rho _{{\rm{gas}}}}/{\rm{\Delta }}\v _{{\rm{gas}}}^{\rm{3}}$, where ρgas is the local gas density. Thus if BHL accretion is substantially altering disc evolution, then the stellar accretion rate acc corrected for the stellar mass dependence − i.e. M˙acc/m2${\dot M_{{\rm{acc}}}}/m_ * ^2 - $ may correlate with local density and anti-correlate with relative velocity of gas. However, finding evidence of such correlations is not necessarily a straight forward task for several reasons. Most importantly, we lack complete (3D) position and velocity information for both the stars and gas. For example, we cannot directly map ISM overdensities along the line-of-sight. The effective noise that this introduces into any statistical signal, alongside numerous other potential complicating factors, is problematic when we consider the modest sample sizes of homogeneously determined accretion rates in most star-forming regions.

In this work, we investigated whether spatial correlations in stellar accretion are evident in existing observational samples. There remain some approaches for searching for such correlations, even in the absence of very large, homogeneous samples of measured accretion rates (and stellar masses) in individual star-forming regions. One approach is to identify overdensities in tracers for ISM dust or gas that coincide with stellar over-densities. Barring improbable projection effects, such alignment would indicate that the local ISM from which the local stars form has not yet dispersed, and therefore stars and gas share a common physical location, not just on the plane of the sky. A second approach for a more dispersed stellar population is to correlate the spatial proximity of neighbouring stars with relative accretion rates. In a sub-structured star-forming region, stars that are closer together in projected separation are also more likely to be close together in (6D) position-velocity space. If the turbulent scale that dominates the fluctuations in local density and velocity of gas on the scale at which they accrete is greater than the interstellar spacing, then stars that are closer together spatially should also have more similar BHL accretion rates.

In the remainder of this work, we present our investigation into possible correlations between stellar accretion and spatial location in the Lupus star-forming region by adopting these approaches. In Section 2 we discuss our approach, including the reasons for choosing Lupus and the datasets we employ. In Section 3 we report on our findings and assess whether there is evidence of spatially dependent stellar accretion within Lupus. We summarise our findings and steps for the future in Section 4.

2 Approach and data

2.1 Normalised accretion rate

As discussed in Section 1, the BHL accretion rate is dependent both on the stellar mass and the properties of the local ISM. In order to search for correlations between stellar accretion acc and the environment it is necessary to isolate the environmental contribution. We therefore define the normalised stellar accretion rate: ˙M˙accm2.$\dot {\cal M} \equiv {\dot M_{{\rm{acc}}}}m_ * ^{ - 2}.$(1)

This factors out the stellar mass dependence of the accretion rate that is expected theoretically (if BHL accretion regulates stellar accretion and until feedback cuts off the accretion flow – see Appendix G of Winter et al. 2024) and found empirically (e.g. Hartmann et al. 2016; Manara et al. 2017; Testi et al. 2022). We are interested in determining whether ˙$\dot {\cal M}$ is randomly distributed or is correlated with position. We therefore require a measurement of both stellar mass and accretion rate for a sufficient sample in a local star-forming region with similar ages.

2.2 Data

In this work, we use the database of homogeneously compiled star and disc properties for the PP7 chapter by Manara et al. (2023). This sample is described in detail in Section 2.5 of that chapter, to which we refer the interested reader. In particular, for this work we required the stellar masses and accretion rates for stars in the Lupus region, which are compiled from the data of Alcalá et al. (2014) and Alcalá et al. (2017). We note that these properties do not have associated individual uncertainties, which are only provided as typical values related to the method used to derive the properties. However, in the context of this work we are interested in inferring correlations, which statistically should only be weakened by the scatter introduced by ignored uncertainties. We limited our sample to stars that have measurements or upper limits for the standardised PP7 accretion rates, stellar mass estimates. This leaves a sample of 74 stars across the whole of Lupus. We further cut this sample by cross-matching with the Luhman (2020) catalogue of 191 kinematically identified Lupus member candidates. Cross-matching with the PPVII catalogue left a sample of 61. We further reduced this to 57 stars that were described as ‘on cloud’ by Luhman (2020), which are those that are spatially co-located with the Lupus 1–4 clouds.

thumbnail Fig. 1

Distribution of candidate Lupus members with Gaia kinematics as identified by Luhman (2020), which we show as hollow shapes (fuchsia circles the ‘on cloud’ sample that coincide with Lupus 1–4, and red squares are the ‘off cloud’ sample). The locations of stars with discs included in the PP7 catalogue of Manara et al. (2023) but without measurements of stellar mass, accretion rate or disc dust mass are shown as lime triangles. Those stars that have stellar masses, dust masses and accretion rates recorded by Manara et al. (2023) are shown as filled coloured circles (or downward triangles for upper limits in accretion rate) that represent the normalised accretion rate according to the colour bar. The greyscale overlay is the IRIS 100 µm map of the region (Miville-Deschênes & Lagache 2005). We mark with cyan stars and labels the locations of systems that exhibit some evidence of interaction with the ISM or high levels of outer disc turbulence, discussed in Section 3.5. We also show a zoom-in on a region centred on Lupus 3, which we define as the ‘clustered’ region (inside the solid red ellipse). Within this circle, we show the ‘on cloud’ sample with fuchsia circles and inner points to make the sample clear in the crowded central region.

2.3 Lupus star-forming region

The Lupus star-forming region is one of the closest low-mass star-forming regions, at a distance of ~160 pc (Lombardi et al. 2008; Galli et al. 2020). It is an ideal laboratory to explore whether accretion rates are correlated with local stellar gas, both because its young stellar population is well-studied, and because the discs have intermediate ages (~1-3 Myr according to Galli et al. 2020, although see also Luhman 2020, arguing for an older typical age of ~6 Myr – we discuss this further in Section 3.3). If these discs are still substantially influenced by BHL accretion, this would strongly support the importance of this process throughout the disc lifetime.

Lupus is composed of several molecular clouds, among which the physical conditions of star formation vary. In many of the clouds, star formation is low density and the stellar population is dispersed. However, Lupus 3 is a much denser region, appearing to be a young, low-mass cluster (Schwartz 1977; Hughes et al. 1994; Nakajima et al. 2000). Galli et al. (2020) find that the ages of stars in Lupus 3 are typically in the range 2.5–3 Myr, commensurate with the stars throughout the complex. However, Rygl et al. (2013) argue that star formation is also more advanced in Lupus 3 with respect to the other clouds, based on the high star formation efficiency and slowing star formation rate inferred from the comparative dearth of proto- and prestellar sources.

We show the distribution of stars surrounding the Lupus clouds in Figure 1. Luhman (2020) demonstrate that this field contains stars associated with not only the Lupus clouds, but also stars associated with V1062 Sco and Upper Centaurus-Lupus. However, the PP7 sample already contains disc-hosting stars that have been previously identified as being associated with one of the Lupus clouds. Indeed, the vast majority of the sample are defined as ‘on cloud’ (i.e. spatially coinciding with the Lupus clouds), with only four stars identified as ‘off cloud’ by Luhman (2020). In our fiducial sample we retain only the ‘on cloud’ stars, but this does not affect our results.

For this work, we performed a separate analysis of the clustered star-forming region centred on Lupus 3. We show a zoom-in panel on the central region of Lupus 3 in the bottom right of Figure 1. The stars within the solid red circle (radius 0.5°) we label the ‘clustered’ sample in this work. In the clustered sample, there are 32 stars with stellar mass and accretion rate constraints in the PP7 catalogue with an associated membership inferred from Gaia kinematics by Luhman (2020). This region is not only densely populated with young stars with respect to the other parts of the cloud complex, but those young stars are centred on a peak in the 100 µm IRIS map (Miville-Deschênes & Lagache 2005). This would suggest that the stellar population shares the same spatial location with residual gas from the star formation process. We may therefore expect enhanced BHL accretion in the dense core, with respect to the lower density outer regions – i.e. a radial gradient in ˙$\dot {\cal M}$.

The stars with accretion rate and stellar mass measurements in the less dense parts of the Lupus complex are distributed among the clouds Lupus 1–4. We label these the ‘distributed’ stars, of which 25 kinematically identified ‘on cloud’ Lupus members have measured masses and accretion rates in the PP7 sample. Most of the distributed stars appear to be similar ages ~2–3 Myr (Rygl et al. 2013), but none are clearly associated with a stellar overdensity that is coincident with a peak in 100 µm emission. Even though we cannot determine the local gas density on a star-by-star basis among the distributed stars, we can assume that the local properties of the turbulent cloud are spatially correlated. In this case, if BHL is contributing to accretion rates, we should expect close neighbouring stars in the distributed population to have more similar ˙$\dot {\cal M}$.

Apart from Lupus, we could also have chosen to study correlations in Taurus or Chameleon I (Cham I), which are also nearby regions that host intermediate age stellar populations. Taurus does not presently contain a sufficient sample of homogeneously derived stellar masses and accretion rates, ruling it out. Cham I does contain such a sample in the PP7 database, however interpreting spatial correlations of accretion rates in this region is not straight forward. Stars in Cham I are distributed along a filament, but there is no clear stellar overdensity (cluster) associated with a monolithic ISM overdensity. Instead, the geometry is complex, and the cloud may be experiencing stellar feedback and dispersal. It is therefore not possible to define ‘clustered’ and ‘distributed’ populations as we have done for Lupus. Such a distinction is necessary since strong local gas density gradients in a bound region may confound efforts to spatially correlate accretion rates. We therefore restricted our attention to the Lupus region. We focus the remainder of this work on the hypothesised correlations within the ‘clustered’ and ‘distributed’ samples in the Lupus star-forming region.

3 Results and discussion

3.1 Clustered stars

We first considered whether there is a radial gradient in ˙$\dot {\cal M}$ among the clustered population. To refine the nominal centre of our population, we adopted the median RA and Dec from the Lupus sample within our chosen spatial region (red circle, Figure 1). We then calculated the angular distance from the centre for each of the stars inside this circle. We show the distribution of normalised accretion rates as a function of separation in Figure 2a. We exclude from our analysis an outlier, Sz 106, which we discuss in Section 3.4. We justified our choice to normalise the accretion rate by m2$m_ * ^2$ based on the theoretical expectation. However, if this scaling is not the ‘true’ physical scaling, this choice may also introduce spurious correlations if the distribution of stellar masses with measured disc and star properties in the PPVII sample are not spatially homogeneous. To ensure the robustness of our findings, we constructed a fiducial sample by restricting the stellar mass range 0.05 M < m*, < 0.4 M, as shown in Figure 2b. These limits are somewhat arbitrary, but motivated by an attempt to maximise the sample while minimising the range in stellar masses. We show in Figure 2b that stellar masses in this fiducial sample are not spatially correlated.

Among our fidudicial sample, from a Spearman’s rank correlation test we find a significant anti-correlation between the separation from the centre and ˙$\dot {\cal M}$, with significance pS = 0.014 (pS = 0.023 if we remove the stellar mass restriction). We also found a similar result by fitting a power-law using LINMIX1 (Kelly 2007), with a form: ˙=A(d1)b,$\dot {\cal M} = A \cdot {\left( {{d \over {1'}}} \right)^b},$(2)

where d is the angular separation from the centre of the cluster. We find normalisation A=3.5×108M1  yr1$A = 3.5 \times {10^{ - 8}}M_ \odot ^{ - 1}\,\,{\rm{y}}{{\rm{r}}^{ - 1}}$ with standard deviation 1.8×108M1  yr1$1.8 \times {10^{ - 8}}M_ \odot ^{ - 1}\,\,{\rm{y}}{{\rm{r}}^{ - 1}}$ and power-law index b = −0.56 ± 0.24 (pink line in Figure 2a). We also find a similar result when including stars of all masses (cyan line in Figure 2a) with b = −0.71 ± 0.23.

We may also ask whether the ISM conditions in Lupus 3 can feasibly drive BHL accretion onto the disc at a rate comparable to the observed stellar accretion rates. While we cannot directly observe the local density and relative gas velocity for individual stars, we can appeal to existing observational constraints to make order of magnitude estimates. The H2 column density in the central regions of Lupus 3 is N0 ~ 1022 cm−2, over a region of radius ~ 1 pc (see Figure A2 of Rygl et al. 2013). Although we do not know the true density profile, from these constraints and by adopting a Plummer density profile we can analytically estimate a reasonable local density and virialised velocity dispersion throughout the cloud. We estimated the normalised BHL accretion rate: ˙BHLπG2ρgas/σv3,${\dot {\cal M}_{{\rm{BHL}}}} \approx \pi {G^2}{\rho _{{\rm{gas}}}}/\sigma _\v ^3,$(3)

where the velocity dispersion συ and ISM density ρgas are both a function of separation from the centre of the cluster. The result of this exercise is shown by red lines in Figure 2a, where we fix N0 ~ 1022 cm−2 and adopt a = 0.5 pc and a = 1 pc. From this simple estimate, we infer that the expected BHL accretion rate is comparable to or slightly exceeds the typical observed stellar accretion rates. This suggests that, even if material accreted onto the disc does not accrete onto the star with 100% efficiency, environmental infall is a plausible driver of disc evolution and stellar accretion.

Our findings suggest that stars in the higher density central regions of Lupus 3 accrete at a higher rate than stars in the outer regions. The BHL accretion rate implied by the approximate density and velocity dispersion in the central regions is comparable to the observed accretion rate. Taken together, these findings support the hypothesis that BHL accretion plays a role in regulating stellar accretion, although the substantial scatter may indicate that internal angular momentum transport processes in the inner disc also play an important role. We discuss caveats to this conclusion in Section 3.4.

thumbnail Fig. 2

Distribution of normalised accretion rates (Figure 2a) and stellar masses (Figure 2b) as a function of the projected distance from the centre of our clustered sample, which is cut at 0.5°. Stars used to calculate correlation coefficients are marked as circles, while excluded stars are shown as squares. The latter lie outside the mass limits shown as horizontal dashed red lines in Figure 2b. All points are coloured by the local 100 µm intensity. Figure 2a also shows the outcome of fitting a power-law to the full sample (cyan) and the stellar mass restricted sample (pink). We also show the estimated rate of BHL accretion based on two Plummer models (red lines, see Section 3.1 for details). We highlight with a red circle Sz 106, which is an outlier and sub-luminous star (see Section 3.4).

3.2 Distributed stars

Outside of the clustered region the stellar population does not have a regular geometry. However, we would still expect stars that are close together to experience more similar BHL accretion rates, since the ISM properties should also be spatially correlated. If BHL accretion influences stellar accretion rates, then such a spatial correlation may also be evident across the measured ˙$\dot {\cal M}$ distribution. We therefore explored the ratio between the normalised accretion rates between stars i and j: log˙r,ij=| log˙ilog˙j |,$\log {\dot {\cal M}_{{\rm{r}},ij}} = \left| {\log {{\dot {\cal M}}_i} - \log {{\dot {\cal M}}_j}} \right|,$(4)

such that ˙r,ij=˙r,ji>1${\dot {\cal M}_{{\rm{r}},ij}} = {\dot {\cal M}_{{\rm{r}},ji}} > 1$ for all i, j. A large (small) value for ˙r,ij${\dot {\cal M}_{{\rm{r}},ij}}$ means that the normalised accretion rates are very different (similar). We calculated this metric for all the pairs (counting each only once). For this calculation, we exclude upper limits leaving 26 accretion rate measurements, although we confirmed that our conclusions are not affected if we take upper limits as measurements. We then divided up the sample based on whether the pairs are closer than ∆θth = 1°. This choice is arbitrary (we subsequently varied it), but must be larger than typical interstellar spacing while substantially smaller than the scale of the whole region. The outcome of this experiment is shown in Figure 3a. We find that close pairs have much more similar ˙r${\dot {\cal M}_{\rm{r}}}$, with significance from a Kologorov-Smirnov (KS) test pKS = 5.7 × 10−3.

Strictly, since each measurement of separation is not independent, the KS test may somewhat overestimate significance. However, we performed a direct bootstrapping experiment by shuffling accretion rates between all individual star locations, measuring the probability of obtaining the observed KS test statistic among these randomly distributed synthetic samples. From this exercise we obtain a similarly significant result pboot = 8.0 × 10−3, with the bootstrap experiments shown as faint dots in Figure 3a. We note that the dispersion among the close stars (blue dots) is larger than the distant stars (red dots) because the subsample size is smaller.

We can further ask if there is a bias in terms of the stellar mass distribution within the observational sample. However, we found no significant correlation between stellar mass and position (Figure 3b). The (insignificant) differences observed between the stellar mass ratios actually go in the wrong direction, with slightly larger differences in stellar masses among nearest neighbours.

Finally, we confirmed that our (arbitrary) choice of ∆θth is not special. In Figure 3c we show KS test results as a function of ∆θth obtain a significant result over a wide range of ∆θth.

We conclude that the normalised accretion rates among the distributed population are correlated, and that this correlation is not related to an inhomogeneous distribution of stellar masses. We discuss the degree to which this finding supports BHL accretion as a driver of stellar accretion in Section 3.4.

thumbnail Fig. 3

Outcome of the correlation experiment for the distributed stellar population across the Lupus cloud complex. Figure 3a shows the cumulative distribution function (CDF) of the ratio of normalised accretion rates between pairs separated by less than (blue) and more than (red) the threshold separation ∆θth = 1°. Faint blue and red points are bootstrapping experiments where we shuffle the normalised accretion rates between star locations. Figure 3b a similar CDF for stellar masses, which are not correlated. Figure 3c shows statistical significance as a function of the threshold ∆θth.

3.3 Age versus environment

3.3.1 Age estimates

If angular momentum transport is mediated by internal processes, either viscous or via MHD winds, we expect stellar accretion to decrease with stellar age (e.g. Manara et al. 2023). Stellar age gradients have been inferred in other star-forming regions, such as the Orion Nebula cluster (over few parsec scales Hillenbrand 1997; Getman et al. 2014; Beccari et al. 2017). It is therefore plausible that stars on the outskirts of Lupus 3 are somewhat older than those in the core. Age differences may be even more substantial among the distributed population, which may form from separate gravitationally unstable regions of the ISM at different times. Age gradients generally may therefore partially or fully explain spatial correlations in accretion rates.

We estimate what kind of age gradients across Lupus might be needed to explain our results as follows. We first consider Upper Sco, in which stars have typical ages inferred by Luhman (2020) ~ 10–12 Myr (similar to Ratzenböck et al. 2023). The typical normalised accretion rate for a star of log m*, ~ −0.5 M in Upper Sco is ˙1.3×109M1  yr1$\dot {\cal M} \approx 1.3 \times {10^{ - 9}}M_ \odot ^{ - 1}\,\,{\rm{y}}{{\rm{r}}^{ - 1}}$, compared to ˙7.9×109M1  yr1$\dot {\cal M} \approx 7.9 \times {10^{ - 9}}M_ \odot ^{ - 1}\,\,{\rm{y}}{{\rm{r}}^{ - 1}}$ in Lupus, based on the fits in Figure 7 of that Almendros-Abad et al. (2024). If disc evolution proceeds somewhat self-similarly in this age range, then we would need a substantial age gradient (a factor ≳2) between the inner and outer regions to reproduce the order of magnitude difference in ˙$\dot {\cal M}$ across our clustered sample (Figure 2a). The factor few systematic differences in ˙$\dot {\cal M}$ across the distributed population (Figure 3a) are more comparable to the differences between Upper Sco and Lupus, but would still require age gradients comparable to the age differences between the two regions. Overall, it would appear that the systematic variation in stellar age across Lupus would have to be of order the age itself in order to explain the observed correlations without appealing to an environmental origin.

Stellar ages inferred from isochrone fitting, for example, are notoriously uncertain for individual stars. It is therefore challenging to test whether age gradients may drive our results by directly estimating these ages. Based on the relative number of stars at different evolutionary stages (defined by the spectral energy distribution – SED – classification), Rygl et al. (2013) suggested that Lupus 3 is more evolved than Lupus 1, which is in turn older than Lupus 4. However, Galli et al. (2020) used photometry from Gaia DR2, 2MASS, AllWISE and Spitzer c2d to fit SEDs to different stellar evolution models, and found that for Lupus 3 and Lupus 4 stellar ages were similar, approximately 2.5–3.5 Myr, while sample sizes were prohibitively small across the other clouds. Luhman (2020) and Ratzenböck et al. (2023) used Gaia optical photometry to estimate ages averaged over all the clouds, both inferring an older age ~6 Myr across clouds 1–4. The typical (model-dependent) uncertainty quoted based on the larger sample by Ratzenböck et al. (2023) is approximately ± 0.5 Myr, although this is a statistical uncertainty and not a physical spread. The large scale systematic spatial age correlations across the Sco-Cen region hint that typically age differences of several Myr between stellar populations correspond to spatial distances of several tens of parsec2. It would therefore seem unlikely for several Myr age differences across the Lupus clouds (on scales ~10–20 pc). To yield the correlations we observe, this would also suggest that the intra-cloud age dispersion for the distributed population must be comparatively small, while in Lupus 3 the dispersion is large. These conditions, while not impossible, appear contrived.

thumbnail Fig. 4

As in Figure 2a, but for normalised disc dust masses rather than accretion rates.

3.3.2 Disc mass as a proxy

We also performed perhaps a more direct test of age differences by considering the properties of our sample. To do so, we reason that if a stellar age gradient were responsible for our findings, we would expect a spatial correlation of the disc masses. Observationally disc (dust) masses decrease on few Myr timescales (e.g. Figure 5 of Manara et al. 2023). We may therefore assume that normalised dust masses are a reasonable proxy for stellar age. Based on this premise, we tested for spatial dust mass correlations across Lupus as follows.

Analogously to accretion rates, we define the normalised dust mass: dust =Mdust /m*2,${{\cal M}_{{\rm{dust }}}} = {M_{{\rm{dust }}}}/m_*^2,$(5)

where Mdust is the dust mass in the disc as in the PPVII dataset. We have used the same normalisation factor m*2$m_*^2$ in this case, despite there being no clear theoretical motivation for such a scaling with stellar mass. However, observationally the accretion timescale (or Mdust/acc) is not strongly dependent on stellar mass (Almendros-Abad et al. 2024).

To test for correlations among the clustered sample, we restricted the range of stellar masses in our fiducial sample, as for stellar accretion rates. We show the outcome of this exercise in Figure 4. We found that, regardless of whether we consider the stellar mass-restricted or entire sample, there is no significant correlation between normalised dust mass and distance from the centre of Lupus 3. Fitting power-laws yielded indices b = −0.14 ± 0.42 and b = −0.25 ± 0.33 respectively. We conclude that disc masses do not appear correlated with separation from cluster centre.

For the distributed population, as for the clustered population, if age were responsible for the spatial correlations of the normalised accretion rates, we would expect this correlation to be evident also among the normalised disc masses. Performing a spatial correlation test for the normalised disc dust masses among the distributed population, we find no such evidence (Figure 5).

We conclude that in Lupus normalised disc (dust) masses do not correlate strongly with spatial location, in contrast to normalised stellar accretion rates. We interpret this as evidence that infall rather than age differences drive the latter correlation. This is because:

  1. Among the broader sample protoplanetary discs, there is no evidence that accretion rates evolve substantially faster. In fact, based on Figure 8 of Almendros-Abad et al. (2024), ℳdust for log m* = −0.5 M decreases from 1.6 × 10−2 at the age of Lupus to 1.8 × 10−3 at the age of Upper Sco (a factor 9 compared to a factor 6 in .$\dot {\cal M}$). This hints that accretion rates may even evolve marginally slower than dust masses (see also Figure 8 of Manara et al. 2023). Therefore if age gradients were resulting in differences in stellar accretion, they should also result in differences in disc (dust) masses.

  2. If the origin of our findings is environmental infall, we may not expect such a strong correlation between spatial location and disc mass. This is because disc mass is a time-integrated quantity, which is dependent on the whole disc evolution history. The present-day BHL accretion rate may not be the same as the typical historic rate. By contrast, stellar accretion may be enhanced by turbulence driven by infall onto the outer disc, which is short-lived (e.g. Winter et al. 2024). The accretion rate is therefore a more direct probe of the local environment than disc mass.

Despite the above arguments, we do not categorically rule out age as a factor in our findings. Interpreting our findings as having an environmental origin rests on the assumption that normalised disc dust masses (continuum fluxes) are at least as closely correlated with stellar age as stellar accretion rates are. While present data would suggest this is a safe assumption, large samples with homogeneous age determinations, or correlation experiments with different probes of the total disc mass may help to strengthen or rebuke the infall scenario. If age is responsible for our findings, this suggests that dust masses must evolve substantially more slowly than accretion rates. In either case, benchmarking models without environmental effects at a specific age to individual star-forming regions may result in incorrect conclusions.

thumbnail Fig. 5

As in Figure 3a, but for normalised disc dust masses rather than accretion rates.

thumbnail Fig. 6

As in Figure 2a, but for a clustered sample contained within a larger radius 0.7° from the cluster centre.

3.4 Sample choice and outliers

Although our results are statistically significant, the sample size remains small. Our approach has also made it necessary to make some choices, some of which can alter our results. In particular, we have chosen a cut-off radius of 0.5° in which to define our clustered sample. This choice is arbitrary, and is made to simply contain the majority of the stars that appear to be approximately spherically distributed around the 100 µm peak.

We therefore tested our conclusions regarding the correlation with separation from the cluster centre against this cut-off radius. First, we find that if we choose instead to limit the sample to those stars within ≲0.4°, the significance drops below 2σ due to the reduced sample size. The power-law fit to the between the separation and .$\dot {\cal M}$ remains similar regardless of the cut-off radius down to ~0.2°.

If we instead increase the cut-off radius that defines the clustered sample to 0.7°, as shown by the dashed red circle in Figure 1, then we include the curious case of J16092697- 3836269, with stellar mass 0.15 M. This is an interesting outlier, which has the largest .$\dot {\cal M}$ across the entire region. It clearly does not fit the trend of the majority of the rest of the clustered sample as shown in Figure 6. It also has an estimated dust mass of Md = 1.34 M and an accretion rate of 7.9 × 10−9 M yr−1. Assuming a gas-to-dust ratio of 100, this implies an accretion timescale of τacc ~ 5 × 104 years, which does not seem compatible with the ~3 Myr age of the majority of the cluster. Given the anomalously short τacc, we are justified in excluding it from the clustered sample. In this case, we retained a significant correlation for our broader cluster sample (Spearman rank test pS = 0.010). If we instead included the outlier, this erases the significance according to a Spearman test. This highlights that, while we have justified our choices, we are limited by a relatively small sample size that means a small number of outliers can alter our results.

There is a second outlier in the clustered region, which we have excluded in our statistical test due to our choice of mass cut. This is Sz 106, which accretes at a very low level despite being positioned close to the centre of the clustered region. Interestingly, this object has been designated ‘sub-luminous’ by Comerón et al. (2003) and Alcalá et al. (2014), who pointed out that the star is significantly less luminous than expected for its spectral type. Alcalá et al. (2014) hypothesise that this is an extinction effect, either due to circumstellar material or an edge on disc. An alternative explanation is an episodic accretion event that could produce a smaller radius, higher temperature, and lower luminosity (Baraffe & Chabrier 2010), but this would also imply enhanced surface gravity that is not observed. Whether or not extinction is the origin for the luminosity dip, the accretion rate may be underestimated for this star and we are therefore justified in considering this anomalous. In fact, the known exceptional state of the star in a sense further supports the apparent trend inferred from Figure 2b, which similarly identifies Sz 106 as an outlier.

3.5 Systems with evidence of infall

In Figure 1, we have highlighted a number of systems which there is some evidence of infall in the literature. While this list may not be complete, they represent interesting case studies that are plausible evidence of ongoing infall:

  • IM Lup has a very extended (~400 au) gas disc with high degree of outer disc turbulence with a turbulent velocity ~0.2–0.6 cs (Paneque-Carreño et al. 2024; Flaherty et al. 2024). IM Lup also has a relatively high normalised accretion rate .=2.7×108  M  yr1$\dot {\cal M} = 2.7 \times {10^{ - 8}}$. While IM Lup has not previously been suggested as a candidate for infall, it does exhibit an extended and asymmetric halo in 12CO out to ~1000 au (Cleeves et al. 2016). This has been previously suggested to be evidence of an externally driven weak pho- toevaporative wind (Haworth et al. 2017), although such asymmetric extended structures could also be interpreted as infalling material. Such infall may also explain the large scale m = 2 spiral structures observed in the mm-continuum (Huang et al. 2018). Modelling is required to examine if infall is a plausible driver of the observed structure.

  • RU Lup has a disc that exhibits large-scale spiral-like structures, extending out to ~1500 au, far beyond the ~120 au Keplerian disc (Huang et al. 2020), as well as reflection nebulosity (Gupta et al. 2023). The spiral structures may be signatures of infalling material (Hennebelle et al. 2017; Kuffmeier et al. 2017). Interestingly, this large scale is comparable to the extent of the halo around IM Lup. Interpreted in the context of BHL accretion, this would suggest a similar accretion radius, corresponding to similar local gas relative velocity. RU Lup also has a very high .=3.3×107M1  yr1$\dot {\cal M} = 3.3 \times {10^{ - 7}}M_ \odot ^{ - 1}$.

  • HD 142527 also exhibits both reflection nebulosity ~700 au spiral structures (Garg et al. 2021), exceeding the disc radius of about 200–300 au (Christiaens et al. 2014; Garg et al. 2021). HD 142527 also shows a misalignment between the inner and outer disks (Bohn et al. 2022), which may also be a consequence of late infall (Thies et al. 2011; Kuffmeier et al. 2021). HD 142527 is not included in the PP7 catalogue of star and disc properties, and we therefore did not consider it in the statistics we present in this work. However, based on the accretion rate inferred by Mendigutía et al. (2014), the normalised accretion rate (.~6×108M1  yr1)$\dot {\cal M}\~6 \times {10^{ - 8}}M_ \odot ^{ - 1}$ is high.

  • HT Lup exhibits reflection nebulosity, as well as a crescentshaped emission extended beyond observed in near-IR polarimetric observations (Garufi et al. 2020) and Herschel maps (Cieza et al. 2013), which may also be a signpost of infall (e.g., Ginski et al. 2021). However, HT Lup is in fact a triple system (A, B and C; Correia et al. 2006), where the overall luminosity is dominated by HT Lup A (Anthonioz et al. 2015), with a mass 1.3 ± 0.2 M (Rich et al. 2022) or 1.32 M in the PP7 database. Interestingly, HT Lup B hosts a disc that is apparently counter-rotating (or possibly perpendicular) with respect to the disc around HT Lup A (Kurtovic et al. 2018). In the PP7 database, HT Lup also only has an upper limit constraint on the stellar accretion rate, such that .<4.4×109M1  yr1$\dot {\cal M} < 4.4 \times {10^{ - 9}}M_ \odot ^{ - 1}$. However, recent observations of HT Lup with VLT/MUSE by Jorquera et al. (2024) have uncovered evidence of variability in the accretion rate onto HT Lup B (which has a dynamical mass ~0.09 M) by a factor several, corresponding to acc ~ 0.7 − 2.9 × 10−9 M yr−1, or a very large .~0.93.7×107M1  yr1$\dot {\cal M}\~0.9 - 3.7 \times {10^{ - 7}}M_ \odot ^{ - 1}$. Accretion onto HT Lup A is only somewhat variable (~30%), with .~3×109M1  yr1$\dot {\cal M}\~3 \times {10^{ - 9}}M_ \odot ^{ - 1}$. The variable accretion rates in this system may be due to mutual gravitational perturbations between stellar components. However, this would suggest a small physical (line-of-sight) separation that would make explaining counter-rotating discs challenging. HT Lup is also very close to 2MASS J15450887- 3417333, which has one of the largest normalised accretion rates (.=8.2×107M1  yr1$\dot {\cal M} = 8.2 \times {10^{ - 7}}M_ \odot ^{ - 1}$) across the PP7 sample.

  • HN Lup exhibits the least compelling evidence of infall due to the absence of necessary ALMA observations. However, it exhibits reflection nebulosity, suggesting that it has substantial local clouds from which it may accrete material (Gupta et al. 2023). Like RU Lup, HN Lup also has relatively large .=5.1×108M1  yr1$\dot {\cal M} = 5.1 \times {10^{ - 8}}M_ \odot ^{ - 1}$.

To summarise, two of the three cases where we have constraints on .$\dot {\cal M}$ from the PP7 database and evidence of infall have particularly large normalised accretion rates, and one component of the third system (HT Lup B) exhibits a high normalised accretion rate and variability. Added to these three systems, HD 142527 has a high .$\dot {\cal M}$, and IM Lup has a more moderate but still substantial .$\dot {\cal M}$. All of these cases are in the distributed population. This is probably a selection bias. For example, reflection nebulae catalogues used in Gupta et al. (2023) are expected to be highly incomplete, particularly in high-density regions (such as Lupus 3) that suffer from higher extinction. A more complete sample of infall candidates is required to statistically correlate the environment with stellar accretion. Nonetheless, the evidence of ongoing infall onto discs around several young stars in Lupus is compelling. Given the short free-fall timescale material at ~1000 au scales (≲104 yr), this would suggest that this process is either relatively common or sustained by the broader environment even up to the ~1–3 Myr age of Lupus.

4 Conclusions

In this work, we have investigated the possibility that accretion rates in the Lupus star-forming region are spatially correlated, which may be evidence of environmentally regulated stellar accretion. To do so, we divided stars into the clustered population in Lupus 3 centred on a substantial remaining gaseous reservoir, and the distributed population across the rest of the Lupus cloud complex, for which stars are not clearly associated with a particular gaseous overdensity. We found statistically significant correlations both within the clustered population and among the distributed population. In the clustered population, the stellar accretion rates are commensurate with the BHL accretion rate we estimate based on the local gas density. Our results cannot be explained by spatial variations of the stellar mass distribution among the observational sample. Age gradients are also disfavoured, since unlike the normalised accretion rates, the normalised disc masses do not exhibit any similar spatial gradient. However, we cannot categorically rule out the role of age gradients, which may be present particularly among the distributed population.

Our results suggest that late stage infall from the ISM onto the protoplanetary disc plays a role in regulating stellar accretion, probably alongside internal processes that drive angular momentum transport in the inner part of the disc. The apparent dependence of star-disc evolution on external environment would not be expected if these latter internal processes were solely responsible for disc evolution. Future studies focused on accurate and homogeneous age determinations or correlating infall occurrence with disc properties may in future strengthen or rebuke this interpretation.

Regardless of their origin, our findings are strong evidence that disc populations should not be considered the result of isolated disc evolution from a well-defined initial condition. Even for star-disc systems at the typical age of Lupus (~1–3 Myr), stellar accretion rates systematically vary across the region.

Acknowledgements

We thank the referee for a considered report, which substantially improved the discussion presented in this paper. AJW has received funding from the European Union’s Horizon 2020 research and innovation programme under the Marie Skłodowska-Curie grant agreement No 101104656. This project has received funding from the European Research Council (ERC) under the European Union’s Horizon research and innovation programme for the ‘PROTOPLANETS’ project, grant No. 101002188, and ‘WANDA’, grant No. 101039452. Views and opinions expressed are however those of the author(s) only and do not necessarily reflect those of the European Union or the European Research Council Executive Agency. Neither the European Union nor the granting authority can be held responsible for them. This work has made use of data from the European Space Agency (ESA) mission Gaia (https://www.cosmos.esa.int/gaia), processed by the Gaia Data Processing and Analysis Consortium (DPAC, https://www.cosmos.esa.int/web/gaia/dpac/consortium). Funding for the DPAC has been provided by national institutions, in particular the institutions participating in the Gaia Multilateral Agreement.

References

  1. Alcalá, J. M., Natta, A., Manara, C. F., et al. 2014, A&A, 561, A2 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  2. Alcalá, J. M., Manara, C. F., Natta, A., et al. 2017, A&A, 600, A20 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  3. Almendros-Abad, V., Manara, C. F., Testi, L., et al. 2024, Evolution of the relation between the mass accretion rate and the stellar and disk mass from brown dwarfs to stars [Google Scholar]
  4. Anthonioz, F., Ménard, F., Pinte, C., et al. 2015, A&A, 574, A41 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  5. Bai, X.-N., & Stone, J. M. 2013, ApJ, 769, 76 [Google Scholar]
  6. Baraffe, I., & Chabrier, G. 2010, A&A, 521, A44 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  7. Bastian, N., Lamers, H. J. G. L. M., de Mink, S. E., et al. 2013, MNRAS, 436, 2398 [CrossRef] [Google Scholar]
  8. Beccari, G., Petr-Gotzens, M. G., Boffin, H. M. J., et al. 2017, A&A, 604, A22 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  9. Benisty, M., Stolker, T., Pohl, A., et al. 2017, A&A, 597, A42 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  10. Benisty, M., Dominik, C., Follette, K., et al. 2023, in Astronomical Society of the Pacific Conference Series, 534, Protostars and Planets VII, eds. S. Inutsuka, Y. Aikawa, T. Muto, K. Tomida, & M. Tamura, 605 [NASA ADS] [Google Scholar]
  11. Boccaletti, A., Di Folco, E., Pantin, E., et al. 2020, A&A, 637, L5 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  12. Bohn, A. J., Benisty, M., Perraut, K., et al. 2022, A&A, 658, A183 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  13. Cazzoletti, P., Manara, C. F., Liu, H. B., et al. 2019, A&A, 626, A11 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  14. Christiaens, V., Casassus, S., Perez, S., van der Plas, G., & Ménard, F. 2014, ApJ, 785, L12 [NASA ADS] [CrossRef] [Google Scholar]
  15. Cieza, L. A., Olofsson, J., Harvey, P. M., et al. 2013, ApJ, 762, 100 [NASA ADS] [CrossRef] [Google Scholar]
  16. Cieza, L. A., Ruíz-Rodríguez, D., Hales, A., et al. 2019, MNRAS, 482, 698 [Google Scholar]
  17. Cleeves, L. I., Öberg, K. I., Wilner, D. J., et al. 2016, ApJ, 832, 110 [Google Scholar]
  18. Codella, C., Podio, L., Simone, M. D., et al. 2024, FAUST XII. Accretion streamers and jets in the VLA 1623-2417 protocluster [Google Scholar]
  19. Comerón, F., Fernández, M., Baraffe, I., Neuhäuser, R., & Kaas, A. A. 2003, A&A, 406, 1001 [CrossRef] [EDP Sciences] [Google Scholar]
  20. Correia, S., Zinnecker, H., Ratzka, T., & Sterzik, M. F. 2006, A&A, 459, 909 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  21. Cuello, N., Ménard, F., & Price, D. J. 2023, Eur. Phys. J. Plus, 138, 11 [NASA ADS] [CrossRef] [Google Scholar]
  22. Deng, H., Ogilvie, G. I., & Mayer, L. 2020, MNRAS, 500, 4248 [CrossRef] [Google Scholar]
  23. Dong, R., Liu, H. B., Cuello, N., et al. 2022, Nat. Astron., 6, 331 [NASA ADS] [CrossRef] [Google Scholar]
  24. Dullemond, C. P., Küffmeier, M., Goicovic, F., et al. 2019, A&A, 628, A20 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  25. Fairbairn, C. W., & Ogilvie, G. I. 2023, MNRAS, 520, 1022 [NASA ADS] [CrossRef] [Google Scholar]
  26. Flaherty, K., Hughes, A. M., Simon, J. B., et al. 2024, MNRAS, 532, 363 [NASA ADS] [CrossRef] [Google Scholar]
  27. Fukagawa, M., Hayashi, M., Tamura, M., et al. 2004, ApJ, 605, L53 [NASA ADS] [CrossRef] [Google Scholar]
  28. Galli, P. A. B., Bouy, H., Olivares, J., et al. 2020, A&A, 643, A148 [EDP Sciences] [Google Scholar]
  29. Garg, H., Pinte, C., Christiaens, V., et al. 2021, MNRAS, 504, 782 [Google Scholar]
  30. Garufi, A., Benisty, M., Pinilla, P., et al. 2018, A&A, 620, A94 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  31. Garufi, A., Avenhaus, H., Pérez, S., et al. 2020, A&A, 633, A82 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  32. Garufi, A., Ginski, C., van Holstein, R. G., et al. 2024, A&A, 685, A53 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  33. Getman, K. V., Feigelson, E. D., & Kuhn, M. A. 2014, ApJ, 787, 109 [NASA ADS] [CrossRef] [Google Scholar]
  34. Ginski, C., Facchini, S., Huang, J., et al. 2021, ApJ, 908, L25 [NASA ADS] [CrossRef] [Google Scholar]
  35. Grady, C. A., Woodgate, B., Bruhweiler, F. C., et al. 1999, ApJ, 523, L151 [NASA ADS] [CrossRef] [Google Scholar]
  36. Gupta, A., Miotello, A., Manara, C. F., et al. 2023, A&A, 670, L8 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  37. Gupta, A., Miotello, A., Williams, J. P., et al. 2024, A&A, 683, A133 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  38. Haisch, Jr., K. E., Lada, E. A., & Lada, C. J. 2001, ApJ, 553, L153 [NASA ADS] [CrossRef] [Google Scholar]
  39. Hales, A. S., Gupta, A., Ruíz-Rodríguez, D., et al. 2024, ApJ, 966, 96 [NASA ADS] [CrossRef] [Google Scholar]
  40. Hanawa, T., Garufi, A., Podio, L., Codella, C., & Segura-Cox, D. 2024, MNRAS, 528, 6581 [NASA ADS] [CrossRef] [Google Scholar]
  41. Hartmann, L., D’Alessio, P., Calvet, N., & Muzerolle, J. 2006, ApJ, 648, 484 [Google Scholar]
  42. Hartmann, L., Herczeg, G., & Calvet, N. 2016, ARA&A, 54, 135 [Google Scholar]
  43. Haworth, T. J., Facchini, S., Clarke, C. J., & Cleeves, L. I. 2017, MNRAS, 468, L108 [CrossRef] [Google Scholar]
  44. Hennebelle, P., Lesur, G., & Fromang, S. 2017, A&A, 599, A86 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  45. Hillenbrand, L. A. 1997, AJ, 113, 1733 [Google Scholar]
  46. Huang, J., Andrews, S. M., Pérez, L. M., et al. 2018, ApJ, 869, L43 [Google Scholar]
  47. Huang, J., Andrews, S. M., Öberg, K. I., et al. 2020, ApJ, 898, 140 [NASA ADS] [CrossRef] [Google Scholar]
  48. Huang, J., Bergin, E. A., Öberg, K. I., et al. 2021, ApJS, 257, 19 [NASA ADS] [CrossRef] [Google Scholar]
  49. Huang, J., Bergin, E. A., Bae, J., Benisty, M., & Andrews, S. M. 2023, ApJ, 943, 107 [NASA ADS] [CrossRef] [Google Scholar]
  50. Hughes, J., Hartigan, P., Krautter, J., & Kelemen, J. 1994, AJ, 108, 1071 [NASA ADS] [CrossRef] [Google Scholar]
  51. Jeffreson, S. M. R., Semenov, V. A., & Krumholz, M. R. 2024, MNRAS, 527, 7093 [Google Scholar]
  52. Jørgensen, J. K., van Dishoeck, E. F., Visser, R., et al. 2009, A&A, 507, 861 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  53. Jorquera, S., Bonnefoy, M., Pérez, L. M., et al. 2024, arXiv e-prints [arXiv:2408.15976] [Google Scholar]
  54. Kelly, B. C. 2007, ApJ, 665, 1489 [Google Scholar]
  55. Kimmig, C. N., & Dullemond, C. P. 2024, A&A, 689, A45 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  56. Klessen, R. S., & Hennebelle, P. 2010, A&A, 520, A17 [NASA ADS] [CrossRef] [Google Scholar]
  57. Kuffmeier, M., Haugbølle, T., & Nordlund, Å. 2017, ApJ, 846, 7 [NASA ADS] [CrossRef] [Google Scholar]
  58. Kuffmeier, M., Goicovic, F. G., & Dullemond, C. P. 2020, A&A, 633, A3 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  59. Kuffmeier, M., Dullemond, C. P., Reissl, S., & Goicovic, F. G. 2021, A&A, 656, A161 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  60. Kuffmeier, M., Jensen, S. S., & Haugbølle, T. 2023, Eur. Phys. J. Plus, 138, 272 [NASA ADS] [CrossRef] [Google Scholar]
  61. Kurtovic, N. T., Pérez, L. M., Benisty, M., et al. 2018, ApJ, 869, L44 [Google Scholar]
  62. Lombardi, M., Lada, C. J., & Alves, J. 2008, A&A, 489, 143 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  63. Luhman, K. L. 2020, AJ, 160, 186 [NASA ADS] [CrossRef] [Google Scholar]
  64. Manara, C. F., Testi, L., Herczeg, G. J., et al. 2017, A&A, 604, A127 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  65. Manara, C. F., Morbidelli, A., & Guillot, T. 2018, A&A, 618, L3 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  66. Manara, C. F., Ansdell, M., Rosotti, G. P., et al. 2023, in Astronomical Society of the Pacific Conference Series, 534, eds. S. Inutsuka, Y. Aikawa, T. Muto, K. Tomida, & M. Tamura, 539 [NASA ADS] [Google Scholar]
  67. Maury, A. J., André, P., Testi, L., et al. 2019, A&A, 621, A76 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  68. Mendigutía, I., Fairlamb, J., Montesinos, B., et al. 2014, ApJ, 790, 21 [CrossRef] [Google Scholar]
  69. Mesa, D., Ginski, C., Gratton, R., et al. 2022, A&A, 658, A63 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  70. Miville-Deschênes, M.-A., & Lagache, G. 2005, ApJS, 157, 302 [Google Scholar]
  71. Montesinos, M., Perez, S., Casassus, S., et al. 2016, ApJ, 823, L8 [Google Scholar]
  72. Nakajima, T., & Golimowski, D. A. 1995, AJ, 109, 1181 [NASA ADS] [CrossRef] [Google Scholar]
  73. Nakajima, Y., Tamura, M., Oasa, Y., & Nakajima, T. 2000, AJ, 119, 873 [NASA ADS] [CrossRef] [Google Scholar]
  74. Padoan, P., Kritsuk, A., Norman, M. L., & Nordlund, Å. 2005, ApJ, 622, L61 [NASA ADS] [CrossRef] [Google Scholar]
  75. Padoan, P., Pan, L., Pelkonen, V.-M., Haugboelle, T., & Nordlund, A. 2024, arXiv e-prints [arXiv:2405.07334] [Google Scholar]
  76. Paneque-Carreño, T., Izquierdo, A. F., Teague, R., et al. 2024, A&A, 684, A174 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  77. Pelkonen, V.-M., Padoan, P., Juvela, M., Haugbølle, T., & Nordlund, Å. 2024, A&A, submitted [arXiv:2405.06520] [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  78. Pérez, S., Hales, A., Liu, H. B., et al. 2020, ApJ, 889, 59 [CrossRef] [Google Scholar]
  79. Pineda, J. E., Segura-Cox, D., Caselli, P., et al. 2020, Nat. Astron., 4, 1158 [NASA ADS] [CrossRef] [Google Scholar]
  80. Ratzenböck, S., Großschedl, J. E., Alves, J., et al. 2023, A&A, 678, A71 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  81. Rich, E. A., Monnier, J. D., Aarnio, A., et al. 2022, AJ, 164, 109 [NASA ADS] [CrossRef] [Google Scholar]
  82. Rygl, K. L. J., Benedettini, M., Schisano, E., et al. 2013, A&A, 549, L1 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  83. Schwartz, R. D. 1977, ApJS, 35, 161 [Google Scholar]
  84. Scicluna, P., Rosotti, G., Dale, J. E., & Testi, L. 2014, A&A, 566, L3 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  85. Su, Z., & Bai, X.-N. 2024, arXiv e-prints [arXiv:2407.12659] [Google Scholar]
  86. Tabone, B., Rosotti, G. P., Cridland, A. J., Armitage, P. J., & Lodato, G. 2022, MNRAS, 512, 2290 [NASA ADS] [CrossRef] [Google Scholar]
  87. Tanious, M., Le Gal, R., Neri, R., et al. 2024, A&A, 687, A92 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  88. Testi, L., Natta, A., Scholz, A., et al. 2016, A&A, 593, A111 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  89. Testi, L., Natta, A., Manara, C. F., et al. 2022, A&A, 663, A98 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  90. Thies, I., Kroupa, P., Goodwin, S. P., Stamatellos, D., & Whitworth, A. P. 2011, MNRAS, 417, 1817 [NASA ADS] [CrossRef] [Google Scholar]
  91. Throop, H. B., & Bally, J. 2008, AJ, 135, 2380 [NASA ADS] [CrossRef] [Google Scholar]
  92. Villenave, M., Stapelfeldt, K. R., Duchêne, G., et al. 2024, ApJ, 961, 95 [NASA ADS] [CrossRef] [Google Scholar]
  93. Williams, J. P., Cieza, L., Hales, A., et al. 2019, ApJ, 875, L9 [Google Scholar]
  94. Winter, A. J., & Clarke, C. J. 2023, MNRAS, 521, 1646 [NASA ADS] [CrossRef] [Google Scholar]
  95. Winter, A. J., & Haworth, T. J. 2022, Eur. Phys. J. Plus, 137, 1132 [NASA ADS] [CrossRef] [Google Scholar]
  96. Winter, A. J., Benisty, M., & Andrews, S. M. 2024, ApJ, 972, L9 [NASA ADS] [CrossRef] [Google Scholar]
  97. Zhang, S., & Zhu, Z. 2024, arXiv e-prints [arXiv:2409.08373] [Google Scholar]

All Figures

thumbnail Fig. 1

Distribution of candidate Lupus members with Gaia kinematics as identified by Luhman (2020), which we show as hollow shapes (fuchsia circles the ‘on cloud’ sample that coincide with Lupus 1–4, and red squares are the ‘off cloud’ sample). The locations of stars with discs included in the PP7 catalogue of Manara et al. (2023) but without measurements of stellar mass, accretion rate or disc dust mass are shown as lime triangles. Those stars that have stellar masses, dust masses and accretion rates recorded by Manara et al. (2023) are shown as filled coloured circles (or downward triangles for upper limits in accretion rate) that represent the normalised accretion rate according to the colour bar. The greyscale overlay is the IRIS 100 µm map of the region (Miville-Deschênes & Lagache 2005). We mark with cyan stars and labels the locations of systems that exhibit some evidence of interaction with the ISM or high levels of outer disc turbulence, discussed in Section 3.5. We also show a zoom-in on a region centred on Lupus 3, which we define as the ‘clustered’ region (inside the solid red ellipse). Within this circle, we show the ‘on cloud’ sample with fuchsia circles and inner points to make the sample clear in the crowded central region.

In the text
thumbnail Fig. 2

Distribution of normalised accretion rates (Figure 2a) and stellar masses (Figure 2b) as a function of the projected distance from the centre of our clustered sample, which is cut at 0.5°. Stars used to calculate correlation coefficients are marked as circles, while excluded stars are shown as squares. The latter lie outside the mass limits shown as horizontal dashed red lines in Figure 2b. All points are coloured by the local 100 µm intensity. Figure 2a also shows the outcome of fitting a power-law to the full sample (cyan) and the stellar mass restricted sample (pink). We also show the estimated rate of BHL accretion based on two Plummer models (red lines, see Section 3.1 for details). We highlight with a red circle Sz 106, which is an outlier and sub-luminous star (see Section 3.4).

In the text
thumbnail Fig. 3

Outcome of the correlation experiment for the distributed stellar population across the Lupus cloud complex. Figure 3a shows the cumulative distribution function (CDF) of the ratio of normalised accretion rates between pairs separated by less than (blue) and more than (red) the threshold separation ∆θth = 1°. Faint blue and red points are bootstrapping experiments where we shuffle the normalised accretion rates between star locations. Figure 3b a similar CDF for stellar masses, which are not correlated. Figure 3c shows statistical significance as a function of the threshold ∆θth.

In the text
thumbnail Fig. 4

As in Figure 2a, but for normalised disc dust masses rather than accretion rates.

In the text
thumbnail Fig. 5

As in Figure 3a, but for normalised disc dust masses rather than accretion rates.

In the text
thumbnail Fig. 6

As in Figure 2a, but for a clustered sample contained within a larger radius 0.7° from the cluster centre.

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.