We interpret observations of intergalactic low ionisation metal absorption systems at redshifts z $\gtrsim$5 in terms of pressure-confined clouds. We find clouds confined by the expected pressure of galactic haloes with masses $11<\log M_h/h^{-1}M_\odot<12$ provide a good description of the column density ratios between low ionisation metal absorbers. Some of the ratios, however, require extending conventional radiative transfer models of irradiated slabs to spherical (or cylindrical) clouds to allow for lines of sight passing outside the cores of the clouds. Moderate depletion of silicon onto dust grains is also indicated in some systems. The chemical abundances inferred span the range between solar and massive-star dominated stellar populations as may arise in starburst galaxies. The typical HI column densities matching the data correspond to Damped Lyman-$\alpha$ Absorbers (DLAs) or sub-DLAs, with sizes of 40 pc to 3 kpc, gas masses $3.5<\log M_c/M_\odot<8$ and metallicites $0.001-0.01Z_\odot$. Such systems continue to pose a challenge for galaxy-scale numerical simulations to reproduce.

Using the Ly$\alpha$ mass assignment scheme (LyMAS), we make theoretical predictions for the 3-dimensional 3-point correlation function (3PCF) of the Ly$\alpha$ forest at redshift $z=2.3$. We bootstrap results from the (100 $h^{-1} \mbox{ Mpc}$)$^3$ Horizon hydrodynamic simulation to a (1 $h^{-1}$ Gpc)$^3$ $N$-body simulation, considering both a uniform UV background (UVB) and a fluctuating UVB sourced by quasars with a comoving $n_q \approx 10^{-5}$ $h^3$ Mpc$^{-3}$ placed either in massive halos or randomly. On scales of $10-30$ $h^{-1} \mbox{ Mpc}$, the flux 3PCF displays hierarchical scaling with the square of the 2PCF, but with an unusual value of $Q \equiv \zeta_{123}/(\xi_{12} \xi_{13} + \xi_{12} \xi_{23} + \xi_{13} \xi_{23}) \approx -4.5$ that reflects the low bias of the Ly$\alpha$ forest and the anti-correlation between mass density and transmitted flux. For halo-based quasars and an ionizing photon mean free path of $\lambda = 300$ $h^{-1} \mbox{ Mpc}$ comoving, UVB fluctuations moderately depress the 2PCF and 3PCF, with cancelling effects on $Q$. For $\lambda = 100$ $h^{-1} \mbox{ Mpc}$ or 50 $h^{-1} \mbox{ Mpc}$, UVB fluctuations substantially boost the 2PCF and 3PCF on large scales, shifting the hierarchical ratio to $Q \approx -3$. We scale our simulation results to derive rough estimate of the 3PCF detectability in observational data sets for the redshift range $z=2.1 - 2.6$. At $r = 10$ $h^{-1} \mbox{ Mpc}$ and 20 $h^{-1} \mbox{ Mpc}$, we predict a signal-to-noise (SNR) of $\sim$ 9 and $\sim$ 7, respectively, for both BOSS and eBOSS, and $\sim$ 37 and $\sim$ 25 for DESI. At $r = 40$ $h^{-1} \mbox{ Mpc}$ the predicted SNR is lower by $\sim$ 3$-$5 times. Measuring the flux 3PCF would be a novel test of the conventional paradigm of the Ly$\alpha$ forest and help separate the contributions of UVB fluctuations and density fluctuations to Ly$\alpha$ forest clustering.

When modelling the ionisation of gas in the intergalactic medium after reionisation, it is standard practice to assume a uniform radiation background. This assumption is not always appropriate; models with radiative transfer show that large-scale ionisation rate fluctuations can have an observable impact on statistics of the Lyman-alpha forest. We extend such calculations to include beaming of sources, which has previously been neglected but which is expected to be important if quasars dominate the ionising photon budget. Beaming has two effects: first, the physical number density of ionising sources is enhanced relative to that directly observed; and second, the radiative transfer itself is altered. We calculate both effects in a hard-edged beaming model where each source has a random orientation, using an equilibrium Boltzmann hierarchy in terms of spherical harmonics. By studying the statistical properties of the resulting ionisation rate and HI density fields at redshift $z\sim 2.3$, we find that the two effects partially cancel each other; combined, they constitute a maximum $5\%$ correction to the power spectrum $P_{\mathrm{HI}}(k)$ at $k=0.04 \, h/\mathrm{Mpc}$. On very large scales ($k<0.01\, h/\mathrm{Mpc}$) the source density renormalisation dominates; it can reduce, by an order of magnitude, the contribution of ionising shot-noise to the intergalactic HI power spectrum. The effects of beaming should be considered when interpreting future observational datasets.

We study how outflows of gas launched from a central galaxy undergoing repeated starbursts propagate through the circumgalactic medium (CGM), using the simulation code RAMSES. We assume that the outflow from the disk can be modelled as a rapidly moving bubble of hot gas at $\mathrm{\sim1\;kpc}$ above disk, then ask what happens as it moves out further into the halo around the galaxy on $\mathrm{\sim 100\;kpc}$ scales. To do this we run 60 two-dimensional simulations scanning over parameters of the outflow. Each of these is repeated with and without radiative cooling, assuming a primordial gas composition to give a lower bound on the importance of cooling. In a large fraction of radiative-cooling cases we are able to form rapidly outflowing cool gas from in situ cooling of the flow. We show that the amount of cool gas formed depends strongly on the 'burstiness' of energy injection; sharper, stronger bursts typically lead to a larger fraction of cool gas forming in the outflow. The abundance ratio of ions in the CGM may therefore change in response to the detailed historical pattern of star formation. For instance, outflows generated by star formation with short, intense bursts contain up to 60 per cent of their gas mass at temperatures $<5 \times 10^4\,\mathrm{K}$; for near-continuous star formation the figure is $\lesssim$ 5 per cent. Further study of cosmological simulations, and of idealised simulations with e.g., metal-cooling, magnetic fields and/or thermal conduction, will help to understand the precise signature of bursty outflows on observed ion abundances.

We use a specific form of the interpolation function in the MOND formalism, which optimally accounts for the internal structure of dwarf spheroidal (dSph) galaxies, to explore the consequences it has on the scaling relations seen in these systems.}} {{The particular form of the interpolation function we used leads to a law of gravity that does not degrade the good fit of the MOND proposal on galactic scales, and in fact, slightly improves the accordance with observations at dSph scales. This formalism yields a good description of gravitational phenomena without the need of invoking any still undetected and hypothetically dominant dark matter, in the weak field regime probed by local dSph galaxies.}} {{Isothermal equilibrium density profiles then yield projected surface density profiles for the local dSph galaxies in very good agreement with observational determinations, for values of the relevant parameters as inferred from recent observations of these Galactic satellites. The observed scaling relations for these systems are also naturally accounted for within the proposed scheme, including a previously unrecognised correlation of the inferred mass-to-light ratios of local dSph's with the ages of their stellar populations, which is natural in modified gravity schemes in the absence of dark matter.}} {The results shed some light on the form that the MOND interpolating function may have in the most challenging regime, which occurs at moderate accelerations and intermediate mass--weighted lengths.