Predictions for the ability of 21-cm interferometric experiments to discriminate Epoch of Reionization (EoR) signal models are typically limited by the simplicity of data models, whereby foreground signals and characteristics of the instrument are often simplified or neglected.To move towards more realistic scenarios, we explore the effects of applying more realistic foreground and instrument models to the 21cm signal, and the ability to estimate astrophysical parameters with these additional complexities. We use a highly-optimized version of \textsc{21cmFAST}, integrated into \textsc{21cmMC}, to generate lightcones of the brightness temperature fluctuation for Bayesian parameter estimation. We include a statistical point-source foreground model and an instrument model based on the Murchison Widefield Array (MWA) scaled in observation time to have an effective sensitivity similar to the future Square Kilometre Array (SKA). We also extend the current likelihood prescription to account for the presence of beam convolution and foregrounds, the 2-Dimensional Power Spectrum (PS), and the correlation of PS modes. We use frequency bands between 150 and 180 MHz to constrain the ionizing efficiency ($\zeta$), the minimum virial temperature of halos ($T_{\mathrm{vir}}$), the soft X-ray emissivity per unit Star Formation Rate (SFR) ($L_X/SFR$ ), and the X-ray energy threshold ($E_0$). We find that the inclusion of realistic foregrounds and instrumental components biases the parameter constraints due to unaccounted for cross-power between the EoR signal, foregrounds and thermal noise. This causes estimates of $\zeta$ to be biased by up to $5\sigma$ but the estimates of $T_{vir}$, L$_X$/SFR and E$_0$ remain unaffected and are all within $1\sigma$.

Detecting a signal from the Epoch of Reionisation (EoR) requires an exquisite understanding of galactic and extra-galactic foregrounds, low frequency radio instruments, instrumental calibration, and data analysis pipelines. In this work we build upon existing work that aims to understand the impact of calibration errors on 21-cm power spectrum (PS) measurements. It is well established that calibration errors have the potential to inhibit EoR detections by introducing additional spectral features that mimic the structure of EoR signals. We present a straightforward way to estimate the impact of a wide variety of modelling residuals in EoR PS estimation. We apply this framework to the specific case of broken dipoles in Murchison Widefield Array (MWA) to understand its effect and estimate its impact on PS estimation. Combining an estimate of the percentage of MWA tiles that have at least one broken dipole (15%-40%) with an analytic description of beam errors induced by such dipoles, we compute the residuals of the foregrounds after calibration and source subtraction. We find that that incorrect beam modelling introduces bias in the 2D-PS on the order of $\sim 10^3\, \mathrm{mK}^2 \,h^{-3}\, \mathrm{Mpc}^{3}$. Although this is three orders of magnitude lower than current lowest limits, it is two orders of magnitude higher than the expected signal. Determining the accuracy of both current beam models and direction dependent calibration pipelines is therefore crucial in our search for an EoR signal.

A very large dynamic range with simultaneous capture of both large- and small-scales in the simulations of cosmic structures is required for correct modelling of many cosmological phenomena, particularly at high redshift. This is not always available, or when it is, it makes such simulations very expensive. We present a novel sub-grid method for modelling low-mass ($10^5\,M_\odot\leq M_{\rm halo}\leq 10^9\,M_\odot$) haloes, which are otherwise unresolved in large-volume cosmological simulations limited in numerical resolution. In addition to the deterministic halo bias that captures the average property, we model its stochasticity that is correlated in time. We find that the instantaneous binned distribution of the number of haloes is well approximated by a log-normal distribution, with overall amplitude modulated by this "temporal correlation bias". The robustness of our new scheme is tested against various statistical measures, and we find that temporally correlated stochasticity generates mock halo data that is significantly more reliable than that from temporally uncorrelated stochasticity. Our method can be applied for simulating processes that depend on both the small- and large-scale structures, especially for those that are sensitive to the evolution history of structure formation such as the process of cosmic reionization. As a sample application, we generate a mock distribution of medium-mass ($ 10^{8} \leq M/M_{\odot} \leq 10^{9}$) haloes inside a 500 Mpc$\,h^{-1}$, $300^3$ grid simulation box. This mock halo catalogue bears a reasonable statistical agreement with a halo catalogue from numerically-resolved haloes in a smaller box, and therefore will allow a very self-consistent sets of cosmic reionization simulations in a box large enough to generate statistically reliable data.