Clustering dependence on Lyα luminosity from MUSE surveys at 3 < z < 6

We investigate the dependence of Ly alpha emitter (LAE) clustering on Ly alpha luminosity and connect the clustering properties of asymptotic to L* LAEs with those of much fainter ones, namely, asymptotic to 0.04L*. We use 1030 LAEs from the MUSE-Wide survey, 679 LAEs from MUSE Deep, and 367 LAEs from the to-date deepest ever spectroscopic survey, the MUSE Extremely Deep Field. All objects have spectroscopic redshifts of 3 < z < 6 and cover a large dynamic range of Ly alpha luminosities: 40.15 < log(LLy alpha/erg s(-1)) < 43.35. We apply the Adelberger et al. K-estimator as the clustering statistic and fit the measurements with state-of-the-art halo occupation distribution (HOD) models. We find that the large-scale bias factor increases weakly with an increasing line luminosity. For the low-luminosity (loghLLy alpha/[erg s(-1)]i = 41.22) and intermediate-luminosity (loghLLy alpha/[erg s-1]i = 41.64) LAEs, we compute consistent bias factors blow = 2.43(+0.15) (-0.15) and b(interm). = 2.42(-0.09)(+0.10), whereas for the high-luminosity (loghLLy alpha/[erg s(-1)]i = 42.34) LAEs we calculated bhigh = 2.65(-0.11.)(+0.13 ) Consequently, high-luminosity LAEs occupy dark matter halos (DMHs) with typical masses of log(Mh/[h-1M(?)]) = 11.09+0.10 -0.09, while low-luminosity LAEs reside in halos of log(Mh/[h-1M(?)]) = 10.77+0.13 -0.15. The minimum masses to host one central LAE, M-min, and (on average) one satellite LAE, M1, also vary with Ly alpha luminosity, growing from log(M-min/[h-1M(?)]) = 10.3+0.2 -0.3 and log(M1/[h-1M(?)]) = 11.7+0.3-0.2 to log(M-min/[h-1M(?)]) = 10.7+0.2-0.3 and log(M1/[h-1M(?)]) = 12.4+0.4 -0.6 from low-to high-luminosity samples, respectively. The satellite fractions are 10% (20%) at 1 sigma (3 sigma) confidence level, supporting a scenario in which DMHs typically host one single LAE. We next bisected the three main samples into disjoint subsets to thoroughly explore the dependence of the clustering properties on LLy alpha. We report a strong (8 sigma) clustering dependence on Ly alpha luminosity, not accounting for cosmic variance effects, where the highest luminosity LAE subsample (log(LLy alpha/erg s(-1))asymptotic to 42.53) clusters more strongly (bhighest = 3.13+0.08-0.15) and resides in more massive DMHs (log(Mh/[h(-1)M(?)]) = 11.43(-0.10)(+0.04 )) than the lowest luminosity one (log(LLy alpha/erg s-1)asymptotic to 40.97), which presents a bias of blowest = 1.79(-0.06)(+0.08) and occupies log(Mh/[h-1M(?)]) = 10.00(+0.12) (-0.09 )halos. We discuss the implications of these results for evolving Ly alpha luminosity functions, halo mass dependent Ly alpha escape fractions, and incomplete reionization signatures.