The hyperfine transition of $^3$He$^+$ at 3.5cm has been thought as a probe of the high-z IGM since it offers a unique insight into the evolution of the helium component of the gas, as well as potentially give an independent constraint on the 21cm signal from neutral hydrogen. In this paper, we use radiative transfer simulations of reionization driven by sources such as stars, X-ray binaries, accreting black holes and shock heated interstellar medium, and simulations of a high-z quasar to characterize the signal and analyze its prospects of detection. We find that the peak of the signal lies in the range 1-50 $\mu$K for both environments, but while around the quasar it is always in emission, in the case of cosmic reionization a brief period of absorption is expected. As the evolution of HeII is determined by stars, we find that it is not possible to distinguish reionization histories driven by more energetic sources. On the other hand, while a bright QSO produces a signal in 21cm that is very similar to the one from a large collection of galaxies, its signature in 3.5cm is very peculiar and could be a powerful probe to identify the presence of the QSO. We analyze the prospects of the signal's detectability using SKA1-mid as our reference telescope. We find that the noise power spectrum dominates over the power spectrum of the signal, although a modest S/N ratio can be obtained when the wavenumber bin width and the survey volume are sufficiently large.

We investigate whether the galaxy and star formation model used for the Auriga simulations can produce a realistic globular cluster (GC) population. We compare statistics of GC candidate star particles in the Auriga haloes with catalogues of the Milky Way (MW) and Andromeda (M31) GC populations. We find that the Auriga simulations do produce sufficient stellar mass for GC candidates at radii and metallicities that are typical for the MW GC system (GCS). We also find varying mass-ratios of the simulated GC candidates relative to the observed mass in the MW and M31 GC systems for different bins of galactocentric radius-metallicity (r$_{\text{gal}}$ -[Fe/H]). Overall, the Auriga simulations produce GC candidates with higher metallicities than the MW and M31 GCS and they are found at larger radii than observed. The Auriga simulations would require bound cluster formation efficiencies higher than ten percent for the metal-poor GC candidates, and those within the Solar radius should experience negligible destruction rates to be consistent with observations. GC candidates in the outer halo, on the other hand, should either have low formation efficiencies, or experience high mass loss for the Auriga simulations to produce a GCS that is consistent with that of the MW or M31. Finally, the scatter in the metallicity as well as in the radial distribution between different Auriga runs is considerably smaller than the differences between that of the MW and M31 GCSs. The Auriga model is unlikely to give rise to a GCS that can be consistent with both galaxies.

We use the Millennium and Millennium-II simulations to illustrate the Tessellation-Level-Tree (TLT), a hierarchical tree structure linking density peaks in a field constructed by voronoi tessellation of the particles in a cosmological N-body simulation. The TLT uniquely partitions the simulation particles into disjoint subsets, each associated with a local density peak. Each peak is a subpeak of a unique higher peak. The TLT can be persistence filtered to suppress peaks produced by discreteness noise. Thresholding a peak's particle list at $\sim 80\left<\rho\right>$ results in a structure similar to a standard friend-of-friends halo and its subhaloes. For thresholds below $\sim 7\left<\rho\right>$, the largest structure percolates and is much more massive than other objects. It may be considered as defining the cosmic web. For a threshold of $5\left<\rho\right>$, it contains about half of all cosmic mass and occupies $\sim 1\%$ of all cosmic volume; a typical external point is then $\sim 7h^{-1}\mathrm{Mpc}$ from the web. We investigate the internal structure and clustering of TLT peaks. Defining the saddle point density $\rho_{\mathrm{lim}}$ as the density at which a peak joins its parent peak, we show the median value of $\rho_{\mathrm{lim}}$ for FoF-like peaks to be similar to the density threshold at percolation. Assembly bias as a function of $\rho_{\mathrm{lim}}$ is stronger than for any known internal halo property. For peaks of group mass and below, the lowest quintile in $\rho_{\mathrm{lim}}$ has $b\approx 0$, and is thus uncorrelated with the mass distribution.

We use hydrodynamics and radiative transfer simulations to study the 21~cm signal around a bright QSO at $z \sim 10$. Due to its powerful UV and X-ray radiation, the QSO quickly increases the extent of the fully ionized bubble produced by the pre-existing stellar type sources, in addition to partially ionize and heat the surrounding gas. As expected, a longer QSO lifetime, $t_{\rm QSO}$, results in a 21~cm signal in emission located at increasingly larger angular radii, $\theta$, and covering a wider range of $\theta$. Similar features can be obtained with a higher galactic emissivity efficiency, $f_{\rm UV}$, so that determining the origin of a large ionized bubble (i.e. QSO vs stars) is not straightforward. Such degeneracy could be reduced by taking advantage of the finite light traveltime effect, which is expected to affect an HII region produced by a QSO differently from one created by stellar type sources. From an observational point of view, we find that the 21 cm signal around a QSO at various $t_{\rm QSO}$ could be detected by SKA1-low with a high signal-noise ratio (S/N). As a reference, for $t_{\rm QSO} = 10\,\rm Myr$, a S/N $\sim 8$ is expected assuming that no pre-heating of the IGM has taken place due to high-$z$ energetic sources, while it can reach value above 10 in case of pre-heating. Observations of the 21~cm signal from the environment of a high-$z$ bright QSO could then be used to set constraints on its lifetime, as well as to reduce the degeneracy between $f_{\rm UV}$ and $t_{\rm QSO}$.

We use publicly available data for the Millennium Simulation to explore the implications of the recent detection of assembly bias and splashback signatures in a large sample of galaxy clusters. These were identified in the SDSS/DR8 photometric data by the redMaPPer algorithm and split into high- and low-concentration subsamples based on the projected positions of cluster members. We use simplified versions of these procedures to build cluster samples of similar size from the simulation data. These match the observed samples quite well and show similar assembly bias and splashback signals. Previous theoretical work has found the logarithmic slope of halo density profiles to have a well-defined minimum whose depth decreases and whose radius increases with halo concentration. Projected profiles for the observed and simulated cluster samples show trends with concentration which are opposite to these predictions. In addition, for high-concentration clusters the minimum slope occurs at significantly smaller radius than predicted. We show that these discrepancies all reflect confusion between splashback features and features imposed on the profiles by the cluster identification and concentration estimation procedures. The strong apparent assembly bias is not reflected in the three-dimensional distribution of matter around clusters. Rather it is a consequence of the preferential contamination of low-concentration clusters by foreground or background groups.

Despite the fact that the mean matter density of the universe has been measured to an accuracy of a few percent within the standard $\Lambda$CDM paradigm, its median density is not known even to order of magnitude. Typical points lie in low-density regions and are not part of a collapsed structure of any scale. Locally, the dark matter distribution is then simply a stretched version of that in the early universe. In this single-stream regime, the distribution of unsmoothed density is sensitive to the initial power spectrum on all scales, in particular on very small scales, and hence to the nature of the dark matter. It cannot be estimated reliably using conventional cosmological simulations because of the enormous dynamic range involved, but a suitable excursion set procedure can be used instead. For the Planck cosmological parameters, a 100 GeV WIMP, corresponding to a free-streaming mass $\sim 10^{-6}M_\odot$, results in a median density of $\sim 4\times 10^{-3}$ in units of the mean density, whereas a 10 $\mu$eV axion with free-streaming mass $\sim 10^{-12}M_\odot$ gives $\sim 3\times 10^{-3}$, and Warm Dark Matter with a (thermal relic) mass of 1 keV gives $\sim 8\times 10^{-2}$. In CDM (but not in WDM) universes, single-stream regions are predicted to be topologically isolated by the excursion set formalism. A test by direct N-Body simulations seems to confirm this prediction, although it is still subject to finite size and resolution effects. Unfortunately, it is unlikely that any of these properties is observable and so suitable for constraining the properties of dark matter.

We introduce a novel technique, called "granulometry", to characterize and recover the mean size and the size distribution of HII regions from 21-cm tomography. The technique is easy to implement, but places the previously not very well defined concept of morphology on a firm mathematical foundation. The size distribution of the cold spots in 21-cm tomography can be used as a direct tracer of the underlying probability distribution of HII region sizes. We explore the capability of the method using large-scale reionization simulations and mock observational data cubes while considering capabilities of SKA1-low and a future extension to SKA2. We show that the technique allows the recovery of the HII region size distribution with a moderate signal-to-noise ratio from wide-field imaging ($\rm SNR\lesssim3$), for which the statistical uncertainty is sample variance dominated. We address the observational requirements on the angular resolution, the field-of-view, and the thermal noise limit for a successful measurement. To achieve a full scientific return from 21-cm tomography and to exploit a synergy with 21-cm power spectra, we suggest an observing strategy using wide-field imaging (several tens of square degrees) by an interferometric mosaicking/multi-beam observation with additional intermediate baselines (~2-4 km).