Angular correlation functions of bright Lyman-break galaxies at $\mathbf{3 \lesssim z \lesssim 5}$

We investigate the clustering of Lyman-break galaxies at redshifts of 3 $\lesssim z \lesssim$ 5 within the COSMOS field by measuring the angular two-point correlation function. Our robust sample of $\sim$60,000 bright ($m_{\rm UV}\lesssim 27$) Lyman-break galaxies was selected based on spectral energy distribution fitting across 14 photometric bands spanning optical and near-infrared wavelengths. We constrained both the 1- and 2-halo terms at separations up to 300 arcsec, finding an excess in the correlation function at scales corresponding to $<20$ kpc, consistent with enhancement due to clumps in the same galaxy or interactions on this scale. We then performed Bayesian model fits on the correlation functions to infer the Halo Occupation Distribution parameters, star formation duty cycle, and galaxy bias in three redshift bins. We examined several cases where different combinations of parameters were varied, showing that our data can constrain the slope of the satellite occupation function, which previous studies have fixed. For an $M_{\rm{UV}}$-limited sub-sample, we found galaxy bias values of $b_g=3.18^{+0.14}_{-0.14}$ at $z\simeq3$, $b_g=3.58^{+0.27}_{-0.29}$ at $z\simeq4$, $b_g=4.27^{+0.25}_{-0.26}$ at $z\simeq5$. The duty cycle values are $0.62^{+0.25}_{-0.26}$, $0.40^{+0.34}_{-0.22}$, and $0.39^{+0.31}_{-0.20}$, respectively. These results suggest that, as the redshift increases, there is a slight decrease in the host halo masses and a shorter timescale for star formation in bright galaxies, at a fixed rest-frame UV luminosity threshold.