DEMOS. II. What the full Lupus demographics truly tell us about protoplanetary disc evolution

Assessing whether disc evolution is driven by viscosity or MHD winds is crucial to determine the conditions in which planets form. The rich statistical information encoded in nearby star-forming regions provides an excellent opportunity to probe disc evolution. In this paper, we apply the Bayesian inference tool DEMOS, designed to exploit the statistical properties of disc populations and infer initial conditions and model parameters, to the Lupus disc population. We consider the observed distributions of stellar mass M, disc mass M (from continuum flux), mass accretion rate and CO disc size R 90%, CO. We constrain all parameters for both the viscous model, with internal photoevaporation, and the MHD winds model, with unprecedented precision. The inferred initial conditions suggest that positive correlations between initial disc properties and M must be present. Both models yield a typical R 0 of 3-5 au for a disc around a 0.3 M star, with positive stellar-mass dependence. We infer α (or α DW) of 3-4 times 10-4 for discs around 0.3 M stars, with an approximately quadratic dependence on M. For the viscous-photoevaporative model, the typical wind mass-loss rate is w about 10-9 M yr-1, scaling with M with an exponent of about 2.5. The inferred mean age is 1.6 Myr for the MHD winds model and 0.9 Myr for the viscous-photoevaporative model. Despite a slight preference for the latter, Bayesian evidence shows that both models reproduce the data comparably well, and neither can be ruled out. However, their best-fit synthetic populations predict substantially different evolutionary pathways. Applying DEMOS to multiple populations, especially evolved ones like Upper Scorpius, will allow to place stronger constraints on disc evolution mechanisms.