HOD modelling of high redshift galaxies using the BLUETIDES simulation
  We construct halo occupation distribution (HOD) models of high redshift ($z
\gtrsim 7.5$) galaxies with $M_{*}>10^8~M_{\odot}/h$ using the BlueTides
hydrodynamic simulation suite, with a particular emphasis on modelling the
small scale / 1-halo clustering ($0.01\lesssim r \lesssim 1~ h^{-1}\rm{Mpc}$).
Similar to low redshift studies, we find that the central and satellite mean
HODs ($\left<N_{\mathrm{cen}}\right>$ and $\left<N_{\mathrm{sat}}\right>$) can
be modeled by a smoothed step function and a power law respectively. The number
density of satellite galaxies is however significantly suppressed compared to
low redshift (satellite fractions drop from $\sim 50 \%$ at $z=0$ to $\lesssim
10 \%$ at $z=7.5$). The mean number of satellites,
$\left<N_{\mathrm{sat}}\right> < 1$ for halo masses below $3 \times 10^{11}
M_{\odot}/h$ (a rare halo at these redshifts). For the radial number density
profiles, satellites with $10^8 \lesssim M^* \lesssim 10^{9} M_{\odot}/h$ in
halos with $M_H \gtrsim 3 \times10^{11} M_{\odot}/h$ are consistent with NFW
(with concentrations $c_{\mathrm{sat}} \sim 10-40$). Within halos of mass
$M_H\lesssim 3 \times 10^{11} M_{\odot}/h$ satellites exhibit a power law
profile with slope -3. Because these halos dominate the small scale clustering,
the resulting 1-halo term is steeper than predicted using standard NFW
profiles. Using this power-law profile for satellites, we can successfully
reproduce the small-scale clustering exhibited by BlueTides galaxies using HOD
modelling. We predict the highest probability of detecting satellites at
$z>7.5$ is around centrals of $M^*\sim 3 \times 10^{10} M_{\odot}/h$ (with
$M^{*}\gtrsim$ a few $10^{7} M_{\odot}/h$ ). This should be achievable with the
James Webb Space Telescope (JWST).
