We introduce the First Light And Reionisation Epoch Simulations (FLARES), a suite of zoom simulations using the EAGLE model. We resimulate a range of overdensities during the Epoch of Reionisation (EoR) in order to build composite distribution functions, as well as explore the environmental dependence of galaxy formation and evolution during this critical period of galaxy assembly. The regions are selected from a large $(3.2 \;\mathrm{cGpc})^{3}$ parent volume, based on their overdensity within a sphere of radius $14\,h^{-1}\;\mathrm{cMpc}$. We then resimulate with full hydrodynamics, and employ a novel weighting scheme that allows the construction of composite distribution functions that are representative of the full parent volume. This significantly extends the dynamic range compared to smaller volume periodic simulations. We present an analysis of the galaxy stellar mass function (GSMF), the star formation rate distribution function (SFRF) and the star forming sequence (SFS) predicted by \flares, and compare to a number of observational and model constraints. We also analyse the environmental dependence over an unprecedented range of overdensity. Both the GSMF and the SFRF exhibit a clear double-Schechter form, up to the highest redshifts ($z = 10$). We also find no environmental dependence of the SFS normalisation. The increased dynamic range probed by FLARES will allow us to make predictions for a number of large area surveys that will probe the EoR in coming years, such as WFIRST and Euclid.

Nebular emission lines associated with galactic HII regions carry information about both physical properties of the ionised gas and the source of ionising photons as well as providing the opportunity of measuring accurate redshifts and thus distances once a cosmological model is assumed. While nebular line emission has been extensively studied at lower redshift there are currently only few constraints within the epoch of reionisation (EoR, $z>6$), chiefly due to the lack of sensitive near-IR spectrographs. However, this will soon change with the arrival of the Webb Telescope providing sensitive near-IR spectroscopy covering the rest-frame UV and optical emission of galaxies in the EoR. In anticipation of Webb we combine the large cosmological hydrodynamical simulation Bluetides with photoionisation modelling to predict the nebular emission line properties of galaxies at $z=8\to 13$. We find good agreement with the, albeit limited, existing direct and indirect observational constraints on equivalent widths though poorer agreement with luminosity function constraints.

We study the synchrotron radio emission from extra-planar regions of star forming galaxies. We use ideal magneto-hydrodynamical (MHD) simulations of a rotating Milky Way-type disk galaxy with distributed star formation sites for three star formation rates (SFRs) (0.3, 3, 30 M$_{\odot}$ yr$^{-1}$). From our simulations, we see emergence of galactic-scale magnetised outflows, carrying gas from the disk. We compare the morphology of the outflowing gas with hydrodynamic (HD) simulations. We look at the spatial distribution of magnetic field in the outflows. Assuming that a certain fraction of gas energy density is converted into cosmic ray energy density, and using information about the magnetic field, we obtain synchrotron emissivity throughout the simulation domain. We generate the surface brightness maps at a frequency of 1.4 GHz. The outflows are more extended in the vertical direction than radial and hence have an oblate shape. We further find that the matter right behind the outer shock, shines brighter in these maps than that above or below. To understand whether this feature can be observed, we produce vertical intensity profiles. We convolve the vertical intensity profile with the typical beam sizes of radio telescopes, for a galaxy located at 10 Mpc (similar to NGC 891) in order to estimate the radio scale height to compare with observations. We find that for our SFRs this feature will lie below the RMS noise limit of instruments. The radio scale height is found to be $\sim 300-1200$ pc , depending on the resolution of the telescope. We relate the advection speed of the outer shock with the surface density of star formation as $\rm{v}_{\rm adv} \propto \Sigma_{\rm SFR}^{0.3}$ which is consistent with earlier observations and analytical estimates.

Galactic outflows produced by stellar feedback are known to be multiphase in nature. Both observations and simulations indicate that the material within several kpc of galactic disk mid-planes consists of warm clouds embedded within a hot wind. A theoretical understanding of the outflow phenomenon, including both winds and fountain flows, requires study of the interactions among thermal phases. We develop a method to quantify these interactions via measurements of mass, momentum, and energy flux exchanges using temporally and spatially averaged quantities and conservation laws. We apply this method to a star-forming ISM MHD simulation based on the TIGRESS framework, for Solar neighbourhood conditions. To evaluate the extent of interactions among the phases, we first examine the validity of the ``ballistic model,'' which predicts trajectories of the warm phase ($5050\,\rm{K}<T<2\times10^4\,\rm{K}$) treated as non-interacting clouds. This model is successful at intermediate vertical velocities ($ 50$ km s$^{-1}$ $\lesssim |v_z| \lesssim 100 $ km s$^{-1}$), but at higher velocities we observe an excess in simulated warm outflow compared to the ballistic model. This discrepancy cannot be fully accounted for by cooling of high-velocity intermediate-temperature ($2\times10^4\,\rm{K}<T<5\times10^5\,\rm{K}$) gas. By examining the fluxes of mass, momentum and energy, we conclude that warm phase gains mass via cooling of the intermediate phase, while momentum transfer occurs from the hot ($T>5\times10^5\,\rm{K}$) to the warm phase. The large energy flux from the hot outflow that is transferred to the warm and intermediate phases is quickly radiated away. A simple interaction model implies an effective warm cloud size in the fountain flow of a few 100~pc, showing that warm-hot flux exchange mainly involves a few large clouds rather than many small ones.

We implement a detailed dust model into the L-Galaxies semi-analytical model which includes: injection of dust by type II and type Ia supernovae (SNe) and AGB stars; grain growth in molecular clouds; and destruction due to supernova-induced shocks, star formation, and reheating. Our grain growth model follows the dust content in molecular clouds and the inter-cloud medium separately, and allows growth only on pre-existing dust grains. At early times, this can make a significant difference to the dust growth rate. Above $z\sim8$, type II SNe are the primary source of dust, whereas below $z\sim8$, grain growth in molecular clouds dominates, with the total dust content being dominated by the latter below $z\sim6$. However, the detailed history of galaxy formation is important for determining the dust content of any individual galaxy. We introduce a fit to the dust-to-metal (DTM) ratio as a function of metallicity and age, which can be used to deduce the DTM ratio of galaxies at any redshift. At $z\lesssim3$, we find a fairly flat mean relation between metallicity and the DTM, and a positive correlation between metallicity and the dust-to-gas (DTG) ratio, in good agreement with the shape and normalisation of the observed relations. We also match the normalisation of the observed stellar mass -- dust mass relation over the redshift range of $0-4$, and to the dust mass function at $z=0$. Our results are important in interpreting observations on the dust content of galaxies across cosmic time, particularly so at high redshift.

We study the effects of mass and energy injection due to OB associations spread across the rotating disc of a Milky Way-type galaxy, with the help of 3D hydrodynamic simulations. We compare the resulting X-ray emission with that produced from the injection of mass and energy from a central region. We find that the predicted X-ray image shows a filamentary structure that arises even in the absence of disc gas inhomogeneity. This structure stems from warm clumps made of disc material being lifted by the injected gas. We show that as much as half of the total X-ray emission comes from regions surrounding warm clumps that are made of a mix of disk and injected gas. This scenario has the potential to explain the origin of the observed extra-planar X-ray emission around star forming galaxies and can be used to understand the observed sublinear relation between the $L_X$ and SFR. We quantify the mass contained in these `bow-shock' regions. We also show that the top-most region of the outer shock above the central area emits harder X-rays than the rest. Further, we find that the mass distribution in different temperature ranges is bimodal, peaking at $10^4\hbox{-}10^5$ K (in warm clumps) and $10^6\hbox{-}10^7$ K (X-ray emitting gas). The mass loading factor is found to decrease with increasing SFR, consistent with previous theoretical estimates and simulations.