We map the lensing-inferred substructure in the first three clusters observed by the Hubble Space Telescope Frontier Fields Initiative (HSTFF): Abell 2744 (z = 0.308), MACSJ0416, (z = 0.396) and MACSJ1149 (z = 0.543). Statistically resolving dark-matter subhaloes down to ~10^{9.5} solar masses, we compare the derived subhalo mass functions (SHMFs) to theoretical predictions from analytical models and with numerical simulations in a Lambda Cold Dark Matter (LCDM) cosmology. Mimicking our observational cluster member selection criteria in the HSTFF, we report excellent agreement in both amplitude and shape of the SHMF over four decades in subhalo mass (10^{9-13} solar masses). Projection effects do not appear to introduce significant errors in the determination of SHMFs from simulations. We do not find evidence for a substructure crisis, analogous to the missing satellite problem in the Local Group, on cluster scales, but rather excellent agreement of the count-matched HSTFF SHMF down to M_{sub halo}/M_{halo} ~ 10^{-5}. However, we do find discrepancies in the radial distribution of sub haloes inferred from HSTFF cluster lenses compared to determinations from simulated clusters. This suggests that although the selected simulated clusters match the HSTFF sample in mass, they do not adequately capture the dynamical properties and complex merging morphologies of these observed cluster lenses. Therefore, HSTFF clusters are likely observed in a transient evolutionary stage that is presently insufficiently sampled in cosmological simulations. The abundance and mass function of dark matter substructure in cluster lenses continues to offer an important test of the LCDM paradigm, and at present we find no tension between model predictions and observations.

Using adiabatic high-resolution numerical simulations we quantify the effect of the streaming motion of baryons with respect to dark matter at the time of recombination on structure formation and evolution. Formally a second order effect, the baryonic stream velocity has proven to have significant impact on dark matter halo abundance, as well as on the gas content and morphology of small galaxy clusters. In this work, we study the impact of stream velocity on the formation and gas content of haloes with masses up to $10^9 M_{\odot}$, an order of magnitude larger than previous studies. We find that the non-zero stream velocity has a sizable impact on the number density of haloes with masses $\lesssim$ few $\times 10^7 M_{\odot}$ up to $z=10$, the final redshift of our simulations. Furthermore, the gas stream velocity induces a suppression of the gas fraction in haloes, which at z=10 is $\sim 10 \%$ for objects with $M\sim10^7M_{\odot}$, as well as a flattening of the gas density profiles in the inner regions of haloes. We further identify and study the formation, in the context of a non-zero stream velocity, of moderately long lived gas dominated structures at intermediate redshifts $10 < z < 20$, which Naoz and Narayan have recently proposed as potential progenitors of globular clusters.