KICC papers

To improve Type Ia supernova (SN Ia) standardisability, the consistency of distance estimates to siblings -- SNe in the same host galaxy -- should be investigated. We present Young Supernova Experiment Pan-STARRS-1 $grizy$ photometry of SN 2021hpr, the third spectroscopically confirmed SN Ia in the high-stellar-mass Cepheid-calibrator galaxy NGC 3147. We analyse NGC 3147's trio of SN Ia siblings: SNe 1997bq, 2008fv and 2021hpr, using a new version of the BayeSN model of SN Ia spectral-energy distributions, retrained simultaneously using optical-NIR $BgVrizYJH$ (0.35--1.8 $\mu$m) data. The distance estimates to each sibling are consistent, with a sample standard deviation $\lesssim$0.01 mag, much smaller than the total intrinsic scatter in the training sample: $\sigma_0\approx0.09$ mag. Fitting normal SN Ia siblings in three additional galaxies, we estimate a $\approx$90% probability that the siblings' intrinsic scatter is smaller than $\sigma_0$. We build a new hierarchical model that fits light curves of siblings in a single galaxy simultaneously; this yields more precise estimates of the common distance and the dust parameters. Fitting the trio for a common dust law shape yields $R_V=2.69\pm0.52$. Our work motivates future hierarchical modelling of more siblings, to tightly constrain their intrinsic scatter, and better understand SN-host correlations. Finally, we estimate the Hubble constant, using a Cepheid distance to NGC 3147, the siblings trio, and 109 Hubble flow ($0.01 < z_{\rm{CMB}} < 0.08$) SNe Ia; marginalising over the siblings' and population's intrinsic scatters, and the peculiar velocity dispersion, yields $H_0=77.9\pm6.5 \text{ km s}^{-1}\text{Mpc}^{-1}$.

We use the BayeSN hierarchical probabilistic SED model to analyse the optical-NIR ($BVriYJH$) light curves of 86 Type Ia supernovae (SNe Ia) from the Carnegie Supernova Project to investigate the SN Ia host galaxy dust law distribution and correlations between SN Ia Hubble residuals and host mass. Our Bayesian analysis simultaneously constrains the mass step and dust $R_V$ population distribution by leveraging optical-NIR colour information. We demonstrate how a simplistic analysis where individual $R_V$ values are first estimated for each SN separately, and then the sample variance of these point estimates is computed, overestimates the $R_V$ population variance $\sigma_R^2$. This bias is exacerbated when neglecting residual intrinsic colour variation beyond that due to light curve shape. Instead, Bayesian shrinkage estimates of $\sigma_R$ are more accurate, with fully hierarchical analysis of the light curves being ideal. For the 75 SNe with low-to-moderate reddening (peak apparent $B-V\leq0.3$), we estimate an $R_V$ distribution with population mean $\mu_R=2.59\pm0.14$, and standard deviation $\sigma_R=0.62\pm0.16$. Splitting this subsample at the median host galaxy mass ($10^{10.57}~\mathrm{M}_\odot$) yields consistent estimated $R_V$ distributions between low- and high-mass galaxies, with $\mu_R=2.79\pm0.18$, $\sigma_R=0.42\pm0.24$, and $\mu_R=2.35\pm0.27$, $\sigma_R=0.74\pm0.36$, respectively. When estimating distances from the full optical-NIR light curves while marginalising over various forms of the dust $R_V$ distribution, a mass step of $\gtrsim0.06$ mag persists in the Hubble residuals at the median host mass.

We study the effect of magnification in the Dark Energy Survey Year 3 analysis of galaxy clustering and galaxy-galaxy lensing, using two different lens samples: a sample of Luminous red galaxies, redMaGiC, and a sample with a redshift-dependent magnitude limit, MagLim. We account for the effect of magnification on both the flux and size selection of galaxies, accounting for systematic effects using the Balrog image simulations. We estimate the impact of magnification on the galaxy clustering and galaxy-galaxy lensing cosmology analysis, finding it to be a significant systematic for the MagLim sample. We show cosmological constraints from the galaxy clustering auto-correlation and galaxy-galaxy lensing signal with different magnifications priors, finding broad consistency in cosmological parameters in $\Lambda$CDM and $w$CDM. However, when magnification bias amplitude is allowed to be free, we find the two-point correlations functions prefer a different amplitude to the fiducial input derived from the image simulations. We validate the magnification analysis by comparing the cross-clustering between lens bins with the prediction from the baseline analysis, which uses only the auto-correlation of the lens bins, indicating systematics other than magnification may be the cause of the discrepancy. We show adding the cross-clustering between lens redshift bins to the fit significantly improves the constraints on lens magnification parameters and allows uninformative priors to be used on magnification coefficients, without any loss of constraining power or prior volume concerns.

Cosmological experiments often employ Bayesian workflows to derive constraints on cosmological and astrophysical parameters from their data. It has been shown that these constraints can be combined across different probes such as Planck and the Dark Energy Survey and that this can be a valuable exercise to improve our understanding of the universe and quantify tension between multiple experiments. However, these experiments are typically plagued by differing systematics, instrumental effects and contaminating signals, which we collectively refer to as `nuisance' components, that have to be modelled alongside target signals of interest. This leads to high dimensional parameter spaces, especially when combining data sets, with > 20 dimensions of which only around 5 correspond to key physical quantities. We present a means by which to combine constraints from different data sets in a computationally efficient manner by generating rapid, reusable and reliable marginal probability density estimators, giving us access to nuisance-free likelihoods. This is possible through the unique combination of nested sampling, which gives us access to Bayesian evidences, and the marginal Bayesian statistics code MARGARINE. Our method is lossless in the signal parameters, resulting in the same posterior distributions as would be found from a full nested sampling run over all nuisance parameters, and typically quicker than evaluating full likelihoods. We demonstrate our approach by applying it to the combination of posteriors from the Dark Energy Survey and Planck.

Reconstruction of gravitational lensing effects in the CMB from current and upcoming surveys is still dominated by temperature anisotropies. Extragalactic foregrounds in temperature maps can induce significant biases in the lensing power spectrum obtained with the standard quadratic estimators. Techniques such as masking cannot remove these foregrounds fully, and the residuals can still lead to large biases if unaccounted for. In this paper, we study the "shear-only" estimator, an example of a class of geometric methods that suppress extragalactic foreground contamination while making only minimal assumptions about foreground properties. The shear-only estimator has only been formulated in the flat-sky limit and so is not easily applied to wide surveys. Here, we derive the full-sky version of the shear-only estimator and its generalisation to an $m=2$ multipole estimator that has improved performance for lensing reconstruction on smaller scales. The multipole estimator is generally not separable, and so is expensive to compute. We explore separable approximations based on a singular-value decomposition, which allow efficient evaluation of the estimator with real-space methods. Finally, we apply these estimators to simulations that include extragalactic foregrounds and verify their efficacy in suppressing foreground biases.

Our view of the last-scattering surface in the cosmic microwave background (CMB) is obscured by secondary anisotropies, sourced by scattering, extragalactic emission and gravitational processes between recombination and observation. Whilst it is established that non-Gaussianity from the correlation between the integrated-Sachs-Wolfe (ISW) effect and gravitational lensing can significantly bias primordial non-Gaussianity (PNG) searches, recent work by Hill (2018) has suggested that other combinations of secondary anisotropies can also produce significant biases. Building on that work, we use the WebSky and Sehgal et al.(2010) simulations to perform an extensive examination of possible biases to PNG measurements for the local, equilateral and orthogonal shapes. For a Planck-like CMB experiment, without foreground cleaning, we find significant biases from cosmic infrared background (CIB)-lensing and thermal Sunyaev-Zel'dovich (tSZ)-lensing bispectra for the local and orthogonal templates, and from CIB and tSZ bispectra for the equilateral template. For future experiments, such as the Simons Observatory, biases from correlations between the ISW effect and the tSZ and CIB will also become important. Finally we investigate whether foreground-cleaning techniques are able to suppress these biases sufficiently. We find that the majority of these biases are effectively suppressed by the internal-linear-combination method and the total bias to Planck-like and SO-like experiments is less than the $1\,\sigma$ statistical error. However, the small total bias arises from the cancellation of several $1\,\sigma$ biases for Planck-like experiments and $2\,\sigma$ biases for SO-like. As this cancellation is likely sensitive to the precise modelling, to ensure robustness against these biases explicit removal methods should be used, likely at the cost of decreased constraining power.

Relying on the dramatic increase in the number of stars with full 6D phase-space information provided by the Gaia Data Release 3, we discover unambiguous signatures of phase-mixing in the stellar halo around the Sun. We show that for the stars likely belonging to the last massive merger, the (v_r,r) distribution contains a series of long and thin chevron-like overdensities. These phase-space sub-structures are predicted to emerge following the dissolution of a satellite, when its tidal debris is given time to wind up, thin out and fold. Additionally, the observed energy and angular momentum (E, L_z) distribution appears more prograde at high energies, possibly revealing the original orbital angular momentum of the in-falling galaxy. The energy distribution of the debris is strongly asymmetric with a peak at low E -- which, we surmise, may be evidence of the dwarf's rapid sinking -- and riddled with wrinkles and bumps. If these small-scale energy inhomogeneities have been seeded during or immediately after the interaction with the Milky Way, and are not due to the spatial restriction of our study, then making use of the (v_r,r) chevrons to constrain the time of the merger becomes cumbersome. Nonetheless, we demonstrate that similar phase-space and (E,L_z) sub-structures are present in numerical simulations of galaxy interactions, both in bespoke N-body runs and in cosmological hydrodynamical zoom-in suites. The remnant traces of the progenitor's disruption and the signatures of the on-going phase-mixing discovered here will not only help to constrain the properties of our Galaxy's most important interaction, but also can be used as a novel tool to map out the Milky Way's current gravitational potential and its perturbations.

The distribution of dark matter halo masses can be accurately predicted in the $\Lambda$CDM cosmology. The presence of a single massive halo or galaxy at a particular redshift, assuming some baryon and stellar fraction for the latter, can therefore be used to test the underlying cosmological model. A number of recent measurements of very large galaxy stellar masses at high redshift ($z > 8$) motivate an investigation into whether any of these objects are in tension with $\Lambda$CDM. We use extreme value statistics to generate confidence regions in the mass-redshift plane for the most extreme mass haloes and galaxies. Tests against numerical models show no tension, neither in their dark matter halo masses nor their galaxy stellar masses. However, we find tentative $> 3\sigma$ tension with recent observational determinations of galaxy masses at high redshift from both HST & JWST, despite using conservative estimates for the stellar fraction ($f_{\star} \sim 1$). Either these galaxies are in tension with $\Lambda$CDM, or there are unaccounted for uncertainties in their stellar mass or redshift estimates.

The fuzzy dark matter (FDM) scenario has received increased attention in recent years due to the small-scale challenges of the vanilla Lambda cold dark matter ($\Lambda$CDM) cosmological model and the lack of robust experimental evidence for any constituent particle. In this study, we use cosmological $N$-body simulations to investigate the high-redshift cosmic web and its responsiveness to an FDM-like power spectrum cutoff in the primordial density perturbations by looking at three distinct properties of virialised FDM dark matter halos as a function of the particle mass $m$. First, compared to $\Lambda$CDM the concentrations of their mass density profiles are lower, peaking at an $m$-dependent halo mass and thus breaking the approximate universality of density profiles in $\Lambda$CDM even further. The halo profiles of the intermediate-to-major and minor-to-major shape parameters are monotonically increasing with ellipsoidal radius in $N$-body simulations of $\Lambda$CDM, yet become non-monotonic owing to baryonic physics at lower redshifts and an FDM-like power spectrum cutoff at higher redshifts. Finally, intrinsic alignment correlations, stemming from the deformation of initially spherically collapsing halos in an ambient gravitational tidal field, become stronger with decreasing FDM particle mass. At $z\sim 4$, we find a $6.4 \sigma$-significance in the fractional differences between the inferred isotropised linear alignment magnitudes $D_{\text{iso}}$ in $\Lambda$CDM and the rather extreme $m=10^{-22}$ eV FDM model. Such FDM-like imprints on the internal properties of virialised halos are strikingly visible in the pristine high-$z$ cosmic web whose evolution is governed largely by linear structure formation physics.

The phase transition between galaxies and quasars is often identified with the rare population of hyper-luminous, hot dust-obscured galaxies. Galaxy formation models predict these systems to grow via mergers, that can deliver large amounts of gas toward their centers, induce intense bursts of star formation and feed their supermassive black holes. Here we report the detection of 24 galaxies emitting Lyman-alpha emission on projected physical scales of about 400 kpc around the hyper-luminous hot dust-obscured galaxy W0410-0913, at redshift z = 3.631, using Very Large Telescope observations. While this indicates that W0410-0913 evolves in a very dense environment, we do not find clear signs of mergers that could sustain its growth. Data suggest that if mergers occurred, as models expect, these would involve less massive satellites, with only a moderate impact on the internal interstellar medium of W0410-0913, which is sustained by a rotationally-supported fast-rotating molecular disk, as Atacama Large Millimeter Array observations suggest.

We present the first detailed study comparing the populations of stellar streams in cosmological simulations to observed Milky Way dwarf galaxy streams. In particular, we compare streams identified around Milky Way analogs in the FIRE-2 simulations to stellar streams observed by the Southern Stellar Stream Spectroscopic Survey (S5). For an accurate comparison between the stream populations, we produce mock Dark Energy Survey (DES) observations of the FIRE streams and estimate the detectability of their tidal tails and progenitors. The number and stellar mass distributions of detectable stellar streams is consistent between observations and simulations. However, there are discrepancies in the distributions of pericenters and apocenters, with the detectable FIRE streams, on average, forming at larger pericenters (out to > 110 kpc) and surviving only at larger apocenters (> 40 kpc) than those observed in the Milky Way. We find that the population of high-stellar mass dwarf galaxy streams in the Milky Way is incomplete. Interestingly, a large fraction of the FIRE streams would only be detected as satellites in DES-like observations, since their tidal tails are too low-surface brightness to be detectable. We thus predict a population of yet-undetected tidal tails around Milky Way satellites, as well as a population of fully undetected low surface brightness stellar streams, and estimate their detectability with the Rubin Observatory. Finally, we discuss the causes and implications of the discrepancies between the stream populations in FIRE and the Milky Way, and explore future avenues for tests of satellite disruption in cosmological simulations.

We analyse the chemical properties of three z~8 galaxies behind the galaxy cluster SMACS J0723.3-7327, observed as part of the Early Release Observations programme of the James Webb Space Telescope (JWST). Exploiting [O III]4363 auroral line detections in NIRSpec spectra, we robustly apply the direct Te method for the very first time at such high redshift, measuring metallicities ranging from extremely metal poor (12+log(O/H)~7) to about one-third solar. We also discuss the excitation properties of these sources, and compare them with local strong-line metallicity calibrations. We find that none of the considered diagnostics match simultaneously the observed relations between metallicity and strong-line ratios for the three sources, implying that a proper re-assessment of the calibrations may be needed at these redshifts. On the mass-metallicity plane, the two galaxies at z~7.6 (log(M*/M_sun) = 8.1, 8.7) have metallicities that are consistent with the extrapolation of the mass-metallicity relation at z~2-3, while the least massive galaxy at z~8.5 (log(M*/M_sun) = 7.8) shows instead a significantly lower metallicity . The three galaxies show different level of offset relative to the Fundamental Metallicity Relation, with two of them (at z~7.6) being marginally consistent, while the z~8.5 source deviating significantly, being probably far from the smooth equilibrium between gas flows, star formation and metal enrichment in place at later epochs.

We present a detailed analysis of the rest-frame optical emission line ratios for three spectroscopically confirmed galaxies at $z>7.5$. The galaxies were identified in the \emph{James Webb Space Telescope} (\emph{JWST}) Early Release Observations field SMACS J0723.3$-$7327. By quantitatively comparing Balmer and oxygen line ratios of these galaxies with various low-redshift "analogue" populations (e.g. Green Peas, Blueberries, etc.), we show that no single analogue population captures the diversity of line ratios of all three galaxies observed at $z>7.5$. We find that S06355 at $z=7.67$ and S10612 at $z=7.66$ are similar to local Green Peas and Blueberries. In contrast, S04590 at $z=8.50$ appears to be significantly different from the other two galaxies, most resembling extremely low-metallicity systems in the local Universe. Perhaps the most striking spectral feature in S04590 is the curiously high [O {\small III}]\,$\lambda4363$/[O {\small III}]\,$\lambda5007$ ratio (RO3) of $0.047$ (or $0.059$ when dust-corrected), implying either extremely high electron temperatures, $>3\times10^4$~K, or gas densities $>10^4\ {\rm cm^{-3}}$. Observed line ratios indicate that this galaxy is unlikely to host an AGN. Using photoionization modelling, we show that the inclusion of high-mass X-ray binaries or a high cosmic ray background in addition to a young, low-metallicity stellar population can provide the additional heating necessary to explain the observed high RO3 while remaining consistent with other observed line ratios. Our models represent a first step at accurately characterising the dominant sources of photoionization and heating at very high redshifts, demonstrating that non-thermal processes may become important as we probe deeper into the Epoch of Reionization.

We present the Sherwood-Relics simulations, a new suite of large cosmological hydrodynamical simulations aimed at modelling the intergalactic medium (IGM) during and after the cosmic reionization of hydrogen. The suite consists of over 200 simulations that cover a wide range of astrophysical and cosmological parameters. It also includes simulations that use a new lightweight hybrid scheme for treating radiative transfer effects. This scheme follows the spatial variations in the ionizing radiation field, as well as the associated fluctuations in IGM temperature and pressure smoothing. It is computationally much cheaper than full radiation hydrodynamics simulations and circumvents the difficult task of calibrating a galaxy formation model to observational constraints on cosmic reionization. Using this hybrid technique, we study the spatial fluctuations in IGM properties that are seeded by patchy cosmic reionization. We investigate the relevant physical processes and assess their impact on the z > 4 Lyman-alpha forest. Our main findings are: (i) Consistent with previous studies patchy reionization causes large scale temperature fluctuations that persist well after the end of reionization, (ii) these increase the Lyman-alpha forest flux power spectrum on large scales, and (iii) result in a spatially varying pressure smoothing that correlates well with the local reionization redshift. (iv) Structures evaporated or puffed up by photoheating cause notable features in the Lyman-alpha forest, such as flat-bottom or double-dip absorption profiles.

The Hubble Frontier Fields represent the opportunity to probe the high-redshift evolution of the main sequence of star-forming galaxies to lower masses than possible in blank fields thanks to foreground lensing of massive galaxy clusters. We use the BEAGLE SED-fitting code to derive stellar masses, $\mathrm{M_{\star}}=\log(M/\mathrm{M_{\odot}})$, SFRs, $\Psi=\log(\psi/\mathrm{M_{\odot}}\,\mathrm{yr}^{-1})$ and redshifts from galaxies within the ASTRODEEP catalogue. We fit a fully Bayesian hierarchical model of the main sequence over $1.25<z<6$ of the form $\Psi = \alpha_\mathrm{9.7}(z) + \beta(\mathrm{M_{\star}}-9.7) + \mathcal{N}(0,\sigma^2)$ while explicitly modelling the outlier distribution. The redshift-dependent intercept at $\mathrm{M_{\star}}=9.7$ is parametrized as $\alpha_\mathrm{9.7}(z) = \log[N (1+z)^{\gamma}] + 0.7$. Our results agree with an increase in normalization of the main sequence to high redshifts that follows the redshift-dependent rate of accretion of gas onto dark matter halos with $\gamma=2.40^{+0.18}_{-0.18}$. We measure a slope and intrinsic scatter of $\beta=0.79^{+0.03}_{-0.04}$ and $\sigma=0.26^{+0.02}_{-0.02}$. We find that the sampling of the SED provided by the combination of filters (Hubble + ground-based Ks-band + Spitzer 3.6 and 4.5 $\mathrm{\mu m}$) is insufficient to constrain $\mathrm{M_{\star}}$ and $\Psi$ over the full dynamic range of the observed main sequence, even at the lowest redshifts studied. While this filter set represents the best current sampling of high-redshift galaxy SEDs out to $z>3$, measurements of the main sequence to low masses and high redshifts still strongly depend on priors employed in SED fitting (as well as other fitting assumptions). Future data-sets with JWST should improve this.

We present constraints on primordial B modes from large angular scale cosmic microwave background polarisation anisotropies measured with the Planck satellite. To remove Galactic polarised foregrounds, we use a Bayesian parametric component separation method, modelling synchrotron radiation as a power law and thermal dust emission as a modified blackbody. This method propagates uncertainties from the foreground cleaning into the noise covariance matrices of the maps. We construct two likelihoods: (i) a semi-analytical cross-spectrum-based likelihood-approximation scheme (momento) and (ii) an exact polarisation-only pixel-based likelihood (pixlike). Since momento is based on cross-spectra it is statistically less powerful than pixlike, but is less sensitive to systematic errors correlated across frequencies. Both likelihoods give a tensor-to-scalar ratio, r, that is consistent with zero from low multipole (2 <= ell < 30) Planck polarisation data. From full-mission maps we obtain r_0.05<0.274, at 95 per cent confidence, at a pivot scale of k = 0.05 Mpc^-1, using pixlike. momento gives a qualitatively similar but weaker 95 per cent confidence limit of r_0.05<0.408.

We present new [OIII] 88 {\mu}m observations of five bright $z \sim 7$ Lyman-break galaxies spectroscopically confirmed by ALMA through the [CII] 158 {\mu}m line, unlike recent [OIII] detections where Lyman-{\alpha} was used. This nearly doubles the sample of Epoch of Reionisation galaxies with robust ($5 \sigma$) detections of [CII] and [OIII]. We perform a multi-wavelength comparison with new deep HST images of the rest-frame UV, whose compact morphology aligns well with [OIII] tracing ionised gas. By contrast, we find more spatially extended [CII] emission likely produced in neutral gas, as indicated by a [NII] 205 {\mu}m non-detection in one source. We find a positive correlation between the equivalent width of the optical [OIII] and H{\beta} lines and the [OIII]/[CII] ratio, as seen in local metal-poor dwarf galaxies. Cloudy models of a nebula of typical density harbouring a young stellar population with a high ionisation parameter appear to adequately reproduce the far-infrared lines. Surprisingly, however, our models fail to reproduce the strength of [OIII] 88 {\mu}m, unless we assume an {\alpha}/Fe enhancement and a near-solar nebular oxygen abundance. On spatially resolved scales, we find [OIII]/[CII] shows a tentative anti-correlation with infrared excess, L_IR/L_UV, also seen on global scales in the local Universe. Finally, we introduce the far-infrared spectral energy distribution fitting code \textsc{mercurius} to show that dust-continuum measurements of one source appear to favour a low dust temperature coupled with a high dust mass. This implies a high stellar metallicity yield and may point towards the need of dust production or grain-growth mechanisms beyond supernovae.

Weak galaxy lensing surveys have consistently reported a lower amplitude for the matter fluctuation spectrum, as measured by the $S_8$ parameter, than expected in the $\Lambda{\rm CDM}$ cosmology favoured by $Planck$. However, the expansion history follows the predictions of the $Planck$ $\Lambda{\rm CDM}$ cosmology to high accuracy, as do measurements of lensing of the cosmic microwave background anisotropies. Redshift space distortion measurements also appear to be consistent with $Planck$ $\Lambda{\rm CDM}$. In this paper, we argue that these observations can be reconciled with the $Planck$ $\Lambda{\rm CDM}$ cosmology if the matter power spectrum is suppressed more strongly on non-linear scales than assumed in analyses of weak galaxy lensing. We demonstrate this point by fitting a one-parameter model, characterising a suppression of the non-linear power spectrum, to the KiDS-1000 weak lensing measurements. Such a suppression could be attributed to new properties of the dark matter that affect non-linear scales, or to a response of the matter fluctuations to baryonic feedback processes that are stronger than expected from recent cosmological simulations. Our proposed explanation can be tested using measurements of the amplitude of the matter fluctuation spectrum on linear scales, in particular via high precision redshift space distortion measurements from forthcoming galaxy and quasar redshift surveys.

There is growing observational evidence for dwarf galaxies hosting active galactic nuclei (AGN), including hints of AGN-driven outflows in dwarfs. However, in the common theoretical model of galaxy formation, efficient supernova (SN) feedback is the tool of choice for regulating star formation in the low-mass regime. In this paper, we present a suite of high-resolution cosmological dwarf zoom-in simulations relaxing the assumption of strong SN feedback, with the goal to determine whether more moderate SN feedback in combination with an efficient AGN could be a suitable alternative. Importantly, we find that there are sufficient amounts of gas to power brief Eddington-limited accretion episodes in dwarfs. This leads to a variety of outcomes depending on the AGN accretion model: from no additional suppression to moderate regulation of star formation to catastrophic quenching. Efficient AGN can drive powerful outflows, depleting the gas reservoir of their hosts via ejective feedback and then maintaining a quiescent state through heating the circumgalactic medium. Moderate AGN outflows can be as efficient as the strong SN feedback commonly employed, leading to star formation regulation and HI gas masses in agreement with observations of field dwarfs. All efficient AGN set-ups are associated with overmassive black holes (BHs) compared to the (heavily extrapolated) observed BH mass - stellar mass scaling relations, with future direct observational constraints in this mass regime being crucially needed. Efficient AGN activity is mostly restricted to high redshifts, with hot, accelerated outflows and high X-ray luminosities being the clearest tell-tale signs for future observational campaigns.

Cosmic Dawn, the onset of star formation in the early universe, can in principle be studied via the 21cm transition of neutral hydrogen, for which a sky-averaged absorption signal, redshifted to MHz frequencies, is predicted to be {\it O}(10-100)\,mK. Detection requires separation of the 21cm signal from bright chromatic foreground emission due to Galactic structure, and the characterisation of how it couples to instrumental response. In this work, we present characterisation of antenna gain patterns for the Large-aperture Experiment to detect the Dark Ages (LEDA) via simulations, assessing the effects of the antenna ground-plane geometries used, and measured soil properties. We then investigate the impact of beam pattern uncertainties on the reconstruction of a Gaussian absorption feature. Assuming the pattern is known and correcting for the chromaticity of the instrument, the foregrounds can be modelled with a log-polynomial, and the 21cm signal identified with high accuracy. However, uncertainties on the soil properties lead to \textperthousand\ changes in the chromaticity that can bias the signal recovery. The bias can be up to a factor of two in amplitude and up to few \% in the frequency location. These effects do not appear to be mitigated by larger ground planes, conversely gain patterns with larger ground planes exhibit more complex frequency structure, significantly compromising the parameter reconstruction. Our results, consistent with findings from other antenna design studies, emphasise the importance of chromatic response and suggest caution in assuming log-polynomial foreground models in global signal experiments.