Observations of the early Universe suggest that reionization was complete by $z\sim6$, however, the exact history of this process is still unknown. One method for measuring the evolution of the neutral fraction throughout this epoch is via observing the Ly$\alpha$ damping wings of high-redshift quasars. In order to constrain the neutral fraction from quasar observations, one needs an accurate model of the quasar spectrum around Ly$\alpha$, after the spectrum has been processed by its host galaxy but before it is altered by absorption and damping in the intervening IGM. In this paper, we present a novel machine learning approach, using artificial neural networks, to reconstruct quasar continua around Ly$\alpha$. Our QSANNdRA algorithm improves the error in this reconstruction compared to the state-of-the-art PCA-based model in the literature by 14.2% on average, and provides an improvement of 6.1% on average when compared to an extension thereof. In comparison with the extended PCA model, QSANNdRA further achieves an improvement of 22.1% and 16.8% when evaluated on low-redshift quasars most similar to the two high-redshift quasars under consideration, ULAS J1120+0641 at $z=7.0851$ and ULAS J1342+0928 at $z=7.5413$, respectively. Using our more accurate reconstructions of these two $z>7$ quasars, we estimate the neutral fraction of the IGM using a homogeneous reionization model and find $\bar{x}_\mathrm{HI} = 0.25^{+0.05}_{-0.05}$ at $z=7.0851$ and $\bar{x}_\mathrm{HI} = 0.60^{+0.11}_{-0.11}$ at $z=7.5413$. Our results are consistent with the literature and favour a rapid end to reionization.

Identifying low-redshift galaxies that emit Lyman Continuum radiation (LyC leakers) is one of the primary, indirect methods of studying galaxy formation in the epoch of reionization. However, not only has it proved challenging to identify such systems, but it also remains uncertain whether the low-redshift LyC leakers are true "analogues" of the sources that reionized the Universe. Here, we use high-resolution cosmological radiation hydrodynamics simulations to examine whether simulated galaxies in the epoch of reionization share similar emission-line properties to observed LyC leakers at $z\sim3$ and $z\sim0$. We find that the simulated galaxies with high LyC escape fractions ($f_{\rm esc}$) often exhibit high O32 and populate the same regions of the R23-O32 plane as $z\sim3$ LyC leakers. However, we show that viewing angle, metallicity, and ionisation parameter can all impact where a galaxy resides on the O32-$f_{\rm esc}$ plane. Based on emission line diagnostics and how they correlate with $f_{\rm esc}$, lower-metallicity LyC leakers at $z\sim3$ appear to be good analogues of reionization-era galaxies. In contrast, we find that identifying low-redshift galaxies based on [SII] deficiencies does not seem to produce true analogues. We use our simulated galaxies to develop multiple new diagnostics to identify LyC leakers using IR and nebular emission lines. We show that our model using only [CII]$_{\rm 158\mu m}$ and [OIII]$_{\rm 88\mu m}$ can identify potential leakers from non-leakers from the local Dwarf Galaxy Survey. Finally, we apply this diagnostic to known high-redshift galaxies and find that MACS1149\_JD1 at $z=9.1$ is the most likely galaxy to be actively contributing to the reionization of the Universe.

We use the SPHINX suite of high-resolution cosmological radiation hydrodynamics simulations to study how spatially and temporally inhomogeneous reionization impacts the baryonic content of dwarf galaxies and cosmic filaments. The SPHINX simulations simultaneously model an inhomogeneous reionization, follow the escape of ionising radiation from thousands of galaxies, and resolve haloes well below the atomic cooling threshold. This makes them an ideal tool for examining how reionization impacts star formation and the gas content of dwarf galaxies. We compare simulations with and without stellar radiation to isolate the effects of radiation feedback from that of supernova, cosmic expansion, and numerical resolution. We find that the gas content of cosmic filaments can be reduced by more than 80% following reionization. The gas inflow rates into haloes with $M_{vir}<10^8M_{\odot}$ are strongly affected and are reduced by more than an order of magnitude compared to the simulation without reionization. A significant increase in gas outflow rates is found for halo masses $M_{vir}<7\times10^7M_{\odot}$. Our simulations show that inflow suppression, rather than photoevaporation, is the dominant mechanism by which the baryonic content of high-redshift dwarf galaxies is regulated. At fixed redshift and halo mass, there is a large scatter in the halo baryon fractions that is entirely dictated by the timing of reionization in the local region surrounding a halo which can change by $\Delta z>3$ at fixed mass. Finally, although the gas content of high-redshift dwarf galaxies is significantly impacted by reionization, we find that most haloes with $M_{vir}<10^8M_{\odot}$ can remain self-shielded and form stars long after reionization, until their local gas reservoir is depleted, suggesting that local group dwarf galaxies do not necessarily exhibit star formation histories that peak prior to $z=6$...

We present the Obelisk project, a cosmological radiation-hydrodynamics simulation following the assembly and reionization of a protocluster progenitor during the first two billions of years from the big bang, down to $z = 3.5$. The simulation resolves haloes down to the atomic cooling limit, and tracks the contribution of different sources of ionization: stars, active galactic nuclei, and collisions. The Obelisk project is designed specifically to study the coevolution of high redshift galaxies and quasars in an environment favouring black hole growth. In this paper, we establish the relative contribution of these two sources of radiation to reionization and their respective role in establishing and maintaining the high redshift ionizing background. Our volume is typical of an overdense region of the Universe and displays star formation rate and black hole accretion rate densities similar to high redshift protoclusters. We find that hydrogen reionization happens inside-out and is completed by $z \sim 6$ in our overdensity, and is predominantly driven by galaxies, while accreting black holes only play a role at $z \sim 4$.

Infrared and nebular lines provide some of our best probes of the physics regulating the properties of the interstellar medium (ISM) at high-redshift. However, interpreting the physical conditions of high-redshift galaxies directly from emission lines remains complicated due to inhomogeneities in temperature, density, metallicity, ionisation parameter, and spectral hardness. We present a new suite of cosmological, radiation-hydrodynamics simulations, each centred on a massive Lyman-break galaxy that resolves such properties in an inhomogeneous ISM. Many of the simulated systems exhibit transient but well defined gaseous disks that appear as velocity gradients in [CII]~158.6$\mu$m emission. Spatial and spectral offsets between [CII]~158.6$\mu$m and [OIII]~88.33$\mu$m are common, but not ubiquitous, as each line probes a different phase of the ISM. These systems fall on the local [CII]-SFR relation, consistent with newer observations that question previously observed [CII]~158.6$\mu$m deficits. Our galaxies are consistent with the nebular line properties of observed $z\sim2-3$ galaxies and reproduce offsets on the BPT and mass-excitation diagrams compared to local galaxies due to higher star formation rate (SFR), excitation, and specific-SFR, as well as harder spectra from young, metal-poor binaries. We predict that local calibrations between H$\alpha$ and [OII]~3727$\AA$ luminosity and galaxy SFR apply up to $z>10$, as do the local relations between certain strong line diagnostics (R23 and [OIII]~5007$\AA$/H$\beta$) and galaxy metallicity. Our new simulations are well suited to interpret the observations of line emission from current (ALMA and HST) and upcoming facilities (JWST and ngVLA).

The scatter in the relationship between the strength of [CII] 158$\mu$m emission and the star formation rate at high-redshift has been the source of much recent interest. Although the relationship is well-established locally, several intensely star-forming galaxies have been found whose [CII] 158$\mu$m emission is either weak, absent or spatially offset from the young stars. Here we present new ALMA data for the two most distant, gravitationally-lensed and spectroscopically-confirmed galaxies, A2744\_YD4 at $z=$8.38 and MACS1149\_JD1 at $z=$9.11, both of which reveal intense [OIII] 88$\mu$m emission. In both cases we provide stringent upper limits on the presence of [CII] 158$\mu$m with respect to [OIII] 88$\mu$m. We review possible explanations for this apparent redshift-dependent [CII] deficit in the context of our recent hydrodynamical simulations. Our results highlight the importance of using several emission line diagnostics with ALMA to investigate the nature of the interstellar medium in early galaxies.

Understanding the escape of Lyman continuum (LyC) and Lyman alpha (Lya) photons from molecular clouds is one of the keys to constraining the reionization history of the Universe and the evolution of galaxies at high redshift. Using a set of radiation-hydrodynamic simulations with adaptive mesh refinement, we investigate how photons propagate and escape from turbulent clouds with different masses, star formation efficiencies (SFEs), and metallicities, as well as with different models of stellar spectra and supernova feedback. We find that the escape fractions in both LyC and Lya are generally increasing with time if the cloud is efficiently dispersed by radiation and supernova feedback. When the total SFE is low (1% of the cloud mass), 0.1-5% of LyC photons leave the metal-poor cloud, whereas the fractions increase to 20-70% in clouds with a 10% SFE. LyC photons escape more efficiently if gas metallicity is lower, if the upper mass limit in the stellar initial mass function is higher, if binary interactions are allowed in the evolution of stars, or if additional strong radiation pressure, such as Lya pressure, is present. As a result, the number of escaping LyC photons can easily vary by a factor of $\sim4$ on cloud scales. The escape fractions of Lya photons are systemically higher (60-80%) than those of LyC photons despite large optical depths at line centre ($\tau_0\sim10^6-10^9$). Scattering of Lya photons is already significant on cloud scales, leading to double-peaked profiles with peak separations of $v_{\rm sep}\sim400\,{\rm km\,s^{-1}}$ during the initial stage of the cloud evolution, while it becomes narrower than $v_{\rm sep} \le 150 \, {\rm km\,s^{-1}}$ in the LyC bright phase. Comparisons with observations of low-redshift galaxies suggest that Lya photons require further interactions with neutral hydrogen to reproduce their velocity offset for a given LyC escape fraction.

The presence of spectroscopically confirmed Balmer breaks in galaxy spectral energy distributions (SEDs) at $z>9$ provides one of the best probes of the assembly history of the first generations of stars in our Universe. Recent observations of the gravitationally lensed source, MACS 1149_JD1 (JD1), indicate that significant amounts of star formation likely occurred at redshifts as high as $z\simeq15$. The inferred stellar mass, dust mass, and assembly history of JD1, or any other galaxy at these redshifts that exhibits a strong Balmer break, can provide a strong test of our best theoretical models from high-resolution cosmological simulations. In this work, we present the results from a cosmological radiation-hydrodynamics simulation of the region surrounding a massive Lyman-break galaxy. For two of our most massive systems, we show that dust preferentially resides in the vicinity of the young stars thereby increasing the strength of the measured Balmer break such that the simulated SEDs are consistent with the photometry of JD1 and two other $z>9$ systems (GN-z10-3 and GN-z9-1) that have proposed Balmer breaks at high redshift. We find strong variations in the shape and luminosity of the SEDs of galaxies with nearly identical stellar and halo masses, indicating the importance of morphology, assembly history, and dust distribution in making inferences on the properties of individual galaxies at high redshifts. Our results stress the importance that dust may play in modulating the observable properties of galaxies, even at the extreme redshifts of $z>9$.

We study the baryonic Tully-Fisher relation (BTFR) at z=0 using 153 galaxies from the SPARC sample. We consider different definitions of the characteristic velocity from HI and H-alpha rotation curves, as well as HI line-widths from single-dish observations. We reach the following results: (1) The tightest BTFR is given by the mean velocity along the flat part of the rotation curve. The orthogonal intrinsic scatter is extremely small (6%) and the best-fit slope is 3.85+/-0.09, but systematic uncertainties may drive the slope from 3.5 to 4.0. Other velocity definitions lead to BTFRs with systematically higher scatters and shallower slopes. (2) We provide statistical relations to infer the flat rotation velocity from HI line-widths or less extended rotation curves (like H-alpha and CO data). These can be useful to study the BTFR from large HI surveys or the BTFR at high redshifts. (3) The BTFR is more fundamental than the relation between angular momentum and galaxy mass (the Fall relation). The Fall relation has about 7 times more scatter than the BTFR, which is merely driven by the scatter in the mass-size relation of galaxies. The BTFR is already the "fundamental plane" of galaxy discs: no value is added with a radial variable as a third parameter.

The mass--velocity--size relation of late-type galaxies decouples into independent correlations between mass and velocity (the Tully-Fisher relation), and between mass and size. This behaviour is different to early-type galaxies which lie on a Fundamental Plane. We study the coupling of the Tully-Fisher and mass-size relations in observations (the SPARC sample) and in empirical galaxy formation models based on halo abundance matching, and rotation curve fits with a hydrodynamically motivated halo profile. We systematically investigate the correlation coefficient between the Tully-Fisher residuals $\Delta V_r$ and mass-size residuals $\Delta R$ as a function of the radius $r$ at which the velocity is measured, and thus present the $\Delta V_r-\Delta R$ relation across rotation curves. We find no significant correlation in either the data or models for any $r$, aside from $r \ll R_\text{eff}$ where baryonic mass dominates. We show that this implies an anticorrelation between galaxy size and halo concentration (or halo mass) at fixed baryonic mass, and provides evidence against the hypothesis that galaxy and halo specific angular momentum are proportional. Finally, we study the $\Delta V_r-\Delta R$ relations produced by the baryons and dark matter separately by fitting halo profiles to the rotation curves. The balance between these components illustrates the "disk-halo conspiracy" required for no overall correlation.

Accurate models of radiative cooling are a fundamental ingredient of modern cosmological simulations. Without cooling, accreted baryons will not efficiently dissipate their energy and collapse to the centres of haloes to form stars. It is well established that local variations in the amplitude and shape of the spectral energy distribution of the radiation field can drastically alter the cooling rate. Here we introduce deepCool, deepHeat, and deepMetal: methods for accurately modelling the total cooling rates, total heating rates, and metal-line only cooling rates of irradiated gas using artificial neural networks. We train our algorithm on a high-resolution cosmological radiation hydrodynamics simulation and demonstrate that we can predict the cooling rate, as measured with the photoionisation code CLOUDY, under the influence of a local radiation field, to an accuracy of ~5%. Our method is computationally and memory efficient, making it suitable for deployment in state-of-the-art radiation hydrodynamics simulations. We show that the circumgalactic medium and diffuse gas surrounding the central regions of a galaxy are most affected by the interplay of radiation and gas, and that standard cooling functions that ignore the local radiation field can incorrectly predict the cooling rate by more than an order of magnitude, indicating that the baryon cycle in galaxies is affected by the influence of a local radiation field on the cooling rate.

Despite their ubiquity, the origin of cosmic magnetic fields remains unknown. Various mechanisms have been proposed for their existence including primordial fields generated by inflation, or amplification and injection by compact astrophysical objects. Separating the potential impact of each magnetogenesis scenario on the magnitude and orientation of the magnetic field and their impact on gas dynamics may give insight into the physics that magnetised our Universe. In this work, we demonstrate that because the induction equation and solenoidal constraint are linear with $B$, the contribution from different sources of magnetic field can be separated in cosmological magnetohydrodynamics simulations and their evolution and influence on the gas dynamics can be tracked. Exploiting this property, we develop a magnetic field tracer algorithm for cosmological simulations that can track the origin and evolution of different components of the magnetic field. We present a suite of cosmological magnetohydrodynamical RAMSES simulations that employ this algorithm where the primordial field strength is varied to determine the contributions of the primordial and supernovae-injected magnetic fields to the total magnetic energy as a function of time and spatial location. We find that, for our specific model, the supernova-injected fields rarely penetrate far from haloes, despite often dominating the total magnetic energy in the simulations. The magnetic energy density from the supernova-injected field scales with density with a power-law slope steeper than 4/3 and often dominates the total magnetic energy inside of haloes. However, the star formation rates in our simulations are not affected by the presence of magnetic fields, for the ranges of primordial field strengths examined. These simulations represent a first demonstration of the magnetic field tracer algorithm (abridged).

We use the photon flux and absorption tracer algorithm presented in Katz et al.~2018, to characterise the contribution of haloes of different mass and stars of different age and metallicity to the reionization of the Universe. We employ a suite of cosmological multifrequency radiation hydrodynamics AMR simulations that are carefully calibrated to reproduce a realistic reionization history and galaxy properties at $z \geq 6$. In our simulations, haloes with mass $10^9{\rm M_{\odot}}h^{-1}<M<10^{10}{\rm M_{\odot}}h^{-1}$, stars with metallicity $10^{-3}Z_{\odot}<Z<10^{-1.5}Z_{\odot}$, and stars with age $3\,{\rm Myr} < t < 10 \, {\rm Myr}$ dominate reionization by both mass and volume. We show that the sources that reionize most of the volume of the Universe by $z=6$ are not necessarily the same sources that dominate the meta-galactic UV background at the same redshift. We further show that in our simulations, the contribution of each type of source to reionization is not uniform across different gas phases. The IGM, CGM, filaments, ISM, and rarefied supernova heated gas have all been photoionized by different classes of sources. Collisional ionisation contributes at both the lowest and highest densities. In the early stages of the formation of individual HII bubbles, reionization proceeds with the formation of concentric shells of gas ionised by different classes of sources, leading to large temperature variations as a function of galacto-centric radius. The temperature structure of individual HII bubbles may thus give insight into the star formation history of the galaxies acting as the first ionising sources. Our explorative simulations highlight how the complex nature of reionization can be better understood by using our photon tracer algorithm.

We present a new empirical relation between galaxy dark matter halo mass (${\rm M_{halo}}$) and the velocity along the flat portion of the rotation curve (${\rm V_{flat}}$), derived from 120 late-type galaxies from the SPARC database. The orthogonal scatter in this relation is comparable to the observed scatter in the baryonic Tully-Fisher relation (BTFR), indicating a tight coupling between total halo mass and galaxy kinematics at $r\ll R_{\rm vir}$. The small vertical scatter in the relation makes it an extremely competitive estimator of total halo mass. We demonstrate that this conclusion holds true for different priors on $M_*/L_{[3.6\mu]}$ that give a tight BTFR, but requires that the halo density profile follows DC14 rather than NFW. We provide additional relations between ${\rm M_{halo}}$ and other velocity definitions at smaller galactic radii (i.e. ${\rm V_{2.2}}$, ${\rm V_{eff}}$, and ${\rm V_{max}}$) which can be useful for estimating halo masses from kinematic surveys, providing an alternative to abundance matching. Furthermore, we constrain the dark matter analog of the Radial Acceleration Relation and also find its scatter to be small, demonstrating the fine balance between baryons and dark matter in their contribution to galaxy kinematics.

In a $\Lambda$CDM cosmology, galaxy formation is a globally inefficient process: it is often the case that far fewer baryons are observed in galaxy disks than expected from the cosmic baryon fraction. The location of these "missing baryons" is unclear. By fitting halo profiles to the rotation curves of galaxies in the SPARC data set, we measure the "missing baryon" mass for individual late-type systems. Assuming that haloes initially accrete the cosmological baryon fraction, we show that the maximum energy available from supernovae is typically not enough to completely eject these "missing baryons" from a halo, but it is often sufficient to heat them to the virial temperature. The energy available from supernovae has the same scaling with galaxy mass as the energy needed to heat or eject the "missing baryons", indicating that the coupling efficiency of the feedback to the ISM may be constant with galaxy virial mass. We further find that the energy available from supernova feedback is always enough to convert a primordial cusp into a core and has magnitude consistent with what is required to heat the "missing baryons" to the virial temperature. Taking a census of the baryon content of galaxies with ${\rm 10^9<M_{vir}/M_{\odot}<10^{12}}$ reveals that $\sim86\%$ of baryons are likely to be in a hot phase surrounding the galaxies and possibly observable in the X-ray, $\sim7\%$ are in the form of cold gas, and $\sim7\%$ are in stars.

The early Universe was composed almost entirely of hydrogen and helium, with only trace amounts of heavy elements. It was only after the first generation of star formation that the Universe became sufficiently polluted to produce a second generation (Population II) of stars which are similar to those in our local Universe. Evidence of massive star cluster formation is nearly ubiquitous among the observed galaxy population and if this mode of star formation occurred at early enough epochs, the higher densities in the early Universe may have caused many of the stars in the cluster to strongly interact. In this scenario, it may be possible to form a very massive star by repeated stellar collisions that may directly collapse into a black hole and form a supermassive black hole seed. In this chapter, we will explore this scenario in detail to understand the dynamics which allow for this process to ensue and measure the probability for this type of seed to represent the supermassive black hole population observed at z > 6.

We present the SPHINX suite of cosmological adaptive mesh refinement simulations, the first radiation-hydrodynamical simulations to simultaneously capture large-scale reionization and the escape of ionizing radiation from thousands of resolved galaxies. Our $5$ and $10$ co-moving Mpc volumes resolve haloes down to the atomic cooling limit and model the inter-stellar medium with better than $\approx10$ pc resolution. The project has numerous goals in improving our understanding of reionization and making predictions for future observations. In this first paper we study how the inclusion of binary stars in computing stellar luminosities impacts reionization, compared to a model that includes only single stars. Owing to the suppression of galaxy growth via strong feedback, our galaxies are in good agreement with observational estimates of the galaxy luminosity function. We find that binaries have a significant impact on the timing of reionization: with binaries, our boxes are $99.9$ percent ionized by volume at $z\approx 7$, while without them our volumes fail to reionize by $z=6$. These results are robust to changes in volume size, resolution, and feedback efficiency. The escape of ionizing radiation from individual galaxies varies strongly and frequently. On average, binaries lead to escape fractions of $\approx 7-10$ percent, about $3.5$ times higher than with single stars only. The higher escape fraction is a result of a shallower decline in ionizing luminosity with age, and is the primary reason for earlier reionization, although the higher integrated luminosity with binaries also plays a sub-dominant role.

We present a new, on-the-fly photon flux and absorption tracer algorithm designed to directly measure the contribution of different source populations to the metagalactic UV background and to the ionisation fraction of gas in the Universe. We use a suite of multifrequency radiation hydrodynamics simulations that are carefully calibrated to reproduce a realistic reionization history and galaxy properties at $z \ge 6$, to disentangle the contribution of photons emitted by different mass haloes and by stars with different metallicities and ages to the UV background during reionization. While at very early cosmic times low mass, metal poor haloes provide most of the LyC photons, their contribution decreases steadily with time. At $z = 6$ it is the photons emitted by massive systems (${\rm M_{halo}}/{\rm M_\odot} > 10^{10} \, {\rm h ^{-1}}$) and by the metal enriched stars ($10^{-3} < Z/Z_{\rm \odot} < 10^{-1.5}$) that provide the largest contribution to the ionising UV background. We demonstrate that there are large variations in the escape fraction depending on the source, with the escape fraction being highest ($\sim 45-60\%$) for photons emitted by the oldest stars that penetrate into the IGM via low opacity channels carved by the ionising photons and supernova from younger stars. Before HII regions begin to overlap, the photoionisation rate strongly fluctuates between different, isolated HII bubbles, depending on the embedded ionising source, which we suggest may result in spatial variations in the properties of dwarf galaxies.

Understanding the origin of strong galactic outflows and the suppression of star formation in dwarf galaxies is a key problem in galaxy formation. Using a set of radiation-hydrodynamic simulations of an isolated dwarf galaxy embedded in a $10^{10}\,M_\odot$ halo, we show that the momentum transferred from resonantly scattered Lyman-$\alpha$ (Lya) photons is an important source of stellar feedback which can shape the evolution of galaxies. We find that Lya feedback suppresses star formation by a factor of two in metal-poor galaxies by regulating the dynamics of star-forming clouds before the onset of supernova explosions (SNe). This is possible because each Lya photon resonantly scatters and imparts 10-300 times greater momentum than in the single scattering limit. Consequently, the number of star clusters predicted in the simulations is reduced by a factor of $\sim 5$, compared to the model without the early feedback. More importantly, we find that galactic outflows become weaker in the presence of strong Lya radiation feedback, as star formation and associated SNe become less bursty. We also examine a model in which radiation field is arbitrarily enhanced by a factor of up to 10, and reach the same conclusion. The typical mass loading factors in our metal-poor dwarf system are estimated to be $\sim5-10$ near the mid plane, while it is reduced to $\sim1$ at larger radii. Finally, we find that the escape of ionizing radiation and hence the reionization history of the Universe is unlikely to be strongly affected by Lya feedback.

We present cosmological, radiation-hydrodynamics simulations of galaxy formation during the epoch of reionization in an effort towards modelling the interstellar medium (ISM) and interpreting ALMA observations. Simulations with and without stellar radiation are compared at large (Mpc), intermediate (tens of kpc), and small (sub kpc) scales. At large scales, the dense regions around galaxies reionize first before UV photons penetrate the voids; however, considerable amounts of neutral gas remain present within the haloes. The spatial distribution of neutral gas is highly dynamic and is anti-correlated with the presence of stars older than a few Myrs. For our specific feedback implementation, most of the metals remain inside the virial radii of haloes and they are proportionally distributed over the ionized and neutral medium by mass. For our most massive galaxy with ${\rm M_h}\sim10^{11}$M$_{\odot}$, the majority of the CII and OI mass are associated with cold neutral clumps. NII is more diffuse and arises in warmer gas while OIII arises in hotter gas with a higher ionization parameter, produced by photo-heating and supernovae. If smaller pockets of high metallicity gas exist in the ISM, the emission from these ions may be observable by ALMA while the low metallicity of the galaxy may cause these systems to fall below the local [CII]-SFR relation. The presence of dust can cause spatial offsets between UV/Ly$\alpha$ and [CII] emission as suggested by the recent observations of Maiolino et al. [OIII] may be spatially offset from both of these components since it arises from a different part of density-temperature phase-space.

Cosmological $N$-body simulations predict dark matter (DM) haloes with steep central cusps (e.g. NFW, Navarro et al. 1996). This contradicts observations of gas kinematics in low-mass galaxies that imply the existence of shallow DM cores. Baryonic processes such as adiabatic contraction and gas outflows can, in principle, alter the initial DM density profile, yet their relative contributions to the halo transformation remain uncertain. Recent high resolution, cosmological hydrodynamic simulations (Di Cintio et al. 2014, DC14) predict that inner density profiles depend systematically on the ratio of stellar to DM mass (M$_*$/M$_{\text{halo}}$). Using a Markov Chain Monte Carlo approach, we test the NFW and the M$_*$/M$_{\text{halo}}$-dependent DC14 halo models against a sample of 147 galaxy rotation curves from the new {\it Spitzer} Photometry and Accurate Rotation Curves (SPARC) data set. These galaxies all have extended H{\small I} rotation curves from radio interferometry as well as accurate stellar mass density profiles from near-infrared photometry. The DC14 halo profile provides markedly better fits to the data compared to the NFW profile. Unlike NFW, the DC14 halo parameters found in our rotation curve fits naturally fall within two standard deviations of the mass-concentration relation predicted by $\Lambda$CDM and the stellar mass-halo mass relation inferred from abundance matching with few outliers. Halo profiles modified by baryonic processes are therefore more consistent with expectations from $\Lambda$ cold dark matter ($\Lambda$CDM) cosmology and provide better fits to galaxy rotation curves across a wide range of galaxy properties than do halo models that neglect baryonic physics. Our results offer a solution to the decade long cusp-core discrepancy.

Reionisation in the early Universe is likely driven by dwarf galaxies. Using cosmological radiation-hydrodynamic simulations, we study star formation and the escape of Lyman continuum (LyC) photons from mini-haloes with $M_{\rm halo} \le 10^8\,M_\odot$. Our simulations include a new thermo-turbulent star formation model, non-equilibrium chemistry, and relevant stellar feedback processes (photoionisation by young massive stars, radiation pressure, and mechanical supernova explosions). We find that feedback reduces star formation very efficiently in mini-haloes, resulting in the stellar mass consistent with the slope and normalisation reported in Kimm \& Cen and the empirical stellar mass-to-halo mass relation derived in the local Universe. Because star formation is stochastic and dominated by a few gas clumps, the escape fraction in mini-haloes is generally determined by radiation feedback (heating due to photo-ionisation), rather than supernova explosions. We also find that the photon number-weighted mean escape fraction in mini-haloes is higher ($\sim20$-$40\%$) than that in atomic-cooling haloes, although the instantaneous fraction in individual haloes varies significantly. The escape fraction from Pop III stars is found to be significant ($\ge10\%$) only when the mass is greater than $\sim$100\,\msun. Based on simple analytic calculations, we show that LyC photons from mini-haloes are, despite their high escape fractions, of minor importance for reionisation due to inefficient star formation. We confirm previous claims that stars in atomic-cooling haloes with masses $10^8\,M_\odot\le M_{\rm halo} \le 10^{11}\,M_\odot$ are likely to be the most important source of reionisation.

We study how runaway stellar collisions in high-redshift, metal-poor star clusters form very massive stars (VMSs) that can directly collapse to intermediate-mass black holes (IMBHs). We follow the evolution of a pair of neighbouring high-redshift mini-haloes with high-resolution, cosmological hydrodynamical zoom-in simulations using the adaptive mesh refinement code RAMSES combined with the non-equilibrium chemistry package KROME. The first collapsing mini-halo is assumed to enrich the central nuclear star cluster (NSC) of the other to a critical metallicity, sufficient for Population II (Pop. II) star formation at redshift $z\approx27$. Using the spatial configuration of the flattened, asymmetrical gas cloud forming in the core of the metal enriched halo, we set the initial conditions for simulations of an initially non-spherical star cluster with the direct summation code NBODY6 which are compared to about 2000 NBODY6 simulations of spherical star clusters for a wide range of star cluster parameters. The final mass of the VMS that forms depends strongly on the initial mass and initial central density of the NSC. For the initial central densities suggested by our RAMSES simulations, VMSs with mass > 400 M$_{\odot}$ can form in clusters with stellar masses of $\approx10^4$ M$_{\odot}$, and this can increase to well over 1000 M$_{\odot}$ for more massive and denser clusters. The high probability we find for forming a VMS in these mini-haloes at such an early cosmic time makes collisional runaway of Pop. II star clusters a promising channel for producing large numbers of high-redshift IMBHs that may act as the seeds of supermassive black holes.

Globular clusters are unique tracers of ancient star formation. We determine the formation efficiencies of globular clusters across cosmic time by modeling the formation and dynamical evolution of the globular cluster population of a Milky Way type galaxy in hierarchical cosmology, using the merger tree from the Via Lactea II simulation. All of the models are constrained to reproduce the observed specific frequency and initial mass function of globular clusters in isolated dwarfs. Globular cluster orbits are then computed in a time varying gravitational potential after they are either accreted from a satellite halo or formed in situ, within the Milky Way halo. We find that the Galactocentric distances and metallicity distribution of globular clusters are very sensitive to the formation efficiencies of globular clusters as a function of redshift and halo mass. Our most accurate models reveal two distinct peaks in the globular cluster formation efficiency at z~2 and z~7-12 and prefer a formation efficiency that is mildly increasing with decreasing halo mass, the opposite of what expected for feedback-regulated star formation. This model accurately reproduces the positions, velocities, mass function, metallicity distribution, and age distribution of globular clusters in the Milky Way and predicts that ~ 40% formed in situ, within the Milky Way halo, while the other ~ 60% were accreted from about 20 satellite dwarf galaxies with Vc > 30 km/s, and about 29% or all globular clusters formed at redshifts z > 7. These results further strengthen the notion that globular cluster formation was an important mode of star formation in high-redshift galaxies and likely played a significant role in the reionization of the intergalactic medium

We utilize Young's algorithm to model the adiabatic compression of the dark matter haloes of galaxies in the THINGS survey to determine the relationship between the halo fit to the rotation curve and the corresponding primordial halo prior to compression. Young's algorithm conserves radial action and angular momentum, resulting in less halo compression than more widely utilized approximations. We find that estimates of the parameters of the NFW haloes fit to the current dark matter distribution systematically overestimate the concentration and underestimate the virial velocity of the corresponding primordial halo. It is the latter that is predicted by dark matter simulations; so accounting for compression is a necessary step for evaluating whether massive galaxies are consistent with dark matter-only simulations. The inferred primordial haloes broadly follow the c-V200 relation expected in a LCDM cosmogony, but often scatter to lower concentrations. We are unable to obtain fits at all for those galaxies whose current dark matter haloes are poorly described by the NFW form. We thus find a mixed bag: some galaxies are reasonably well described by adiabatic compression within a primordial NFW halo, while others require additional mechanisms that reduce the density of dark matter below the primordial initial condition.

We examine the formation of clusters of galaxies in numerical simulations of a QUMOND cosmogony with massive sterile neutrinos. Clusters formed in these exploratory simulations develop higher velocities than those found in {\Lambda}CDM simulations. The bulk motions of clusters attain about 1000 km/s by low redshift, comparable to observations whereas {\Lambda}CDM simulated clusters tend to fall short. Similarly, high pairwise velocities are common in cluster-cluster collisions like the Bullet cluster. There is also a propensity for the most massive clusters to be larger in QUMOND and to appear earlier than in {\Lambda}CDM, potentially providing an explanation for 'pink elephants' like El Gordo. However, it is not obvious that the cluster mass function can be recovered.

The ages of globular clusters in our own Milky Way are known with precision of about \pm 1 Gyr, hence their formation at redshifts z>~3 and their role in hierarchical cosmology and the reionization of the intergalactic medium remain relatively undetermined. Here we analyze the effect of globular cluster formation on the observed rest-frame UV luminosity functions (LFs) and UV continuum slopes of high redshift galaxies in the Hubble Ultra Deep Fields. We find that the majority of present day globular clusters have formed during two distinct epochs: at redshifts z ~ 2-3 and at redshifts z>~6. The birth of proto-GC systems produce the steep, faint-end slopes of the galaxy LFs, and because the brightness of proto-GCs fades 5 Myrs after their formation, their blue colors are in excellent agreement with observations. Our results suggest that: i) the bulk of the old globular cluster population with estimated ages >~12 Gyr (about 50% of the total population), formed in the relatively massive dwarf galaxies at redshifts z>~6; ii) proto-GC formation was an important mode of star formation in those dwarf galaxies, and likely dominated the reionization process. Another consequence of this scenario is that some of the most massive Milky Way satellites may be faint and yet undiscovered because tidal stripping of a dominant GC population precedes significant stripping of the dark matter halos of these satellites. This scenario may alleviate some remaining tensions between CDM simulations and observations.