Low radio frequency experiments performed on Earth are contaminated by both ionospheric effects and radio frequency interference (RFI) from Earth-based sources. The lunar farside provides a unique environment above the ionosphere where RFI is heavily attenuated by the presence of the Moon. We present electrodynamics simulations of the propagation of radio waves around and through the Moon in order to characterize the level of attenuation on the farside. The simulations are performed for a range of frequencies up to 100 kHz, assuming a spherical lunar shape with an average, constant density. Additionally, we investigate the role of the topography and density profile of the Moon in the propagation of radio waves and find only small effects on the intensity of RFI. Due to the computational demands of performing simulations at higher frequencies, we propose a model for extrapolating the width of the quiet region above 100 kHz that also takes into account height above the lunar surface as well as the intensity threshold chosen to define the quiet region. This model, which we make publicly available through a Python package, allows the size of the radio quiet region to be easily calculated both in orbit or on the surface, making it directly applicable for lunar satellites as well as surface missions.

The assumptions inherent to global 21-cm signal analyses are rarely delineated. In this paper, we formulate a general list of suppositions underlying a given claimed detection of the global 21-cm signal. Then, we specify the form of these assumptions for two different analyses: 1) the one performed by the EDGES team showing an absorption trough in brightness temperature that they modeled separately from the sky foreground and 2) a new, so-called Minimum Assumption Analysis (MAA), that makes the most conservative assumptions possible for the signal. We show fits using the EDGES analysis on various beam-weighted foreground simulations from the EDGES latitude with no signal added. Depending on the beam used, these simulations produce large false troughs due to the invalidity of the foreground model to describe the combination of beam chromaticity and the shape of the Galactic plane in the sky, the residuals of which are captured by the ad hoc flattened Gaussian signal model. On the other hand, the MAA provides robust fits by including many spectra at different time bins and allowing any possible 21-cm spectrum to be modeled exactly. We present uncertainty levels and example signal reconstructions found with the MAA for different numbers of time bins. With enough time bins, one can determine the true 21-cm signal with the MAA to $<10$ times the noise level.

When using valid foreground and signal models, the uncertainties on extracted signals in global 21-cm signal experiments depend principally on the overlap between signal and foreground models. In this paper, we investigate two strategies for decreasing this overlap: (i) utilizing time dependence by fitting multiple drift-scan spectra simultaneously and (ii) measuring all four Stokes parameters instead of only the total power, Stokes I. Although measuring polarization requires different instruments than are used in most existing experiments, all existing experiments can utilize drift-scan measurements merely by averaging their data differently. In order to evaluate the increase in constraining power from using these two techniques, we define a method for connecting Root-Mean-Square (RMS) uncertainties to probabilistic confidence levels. Employing simulations, we find that fitting only one total power spectrum leads to RMS uncertainties at the few K level, while fitting multiple time-binned, drift-scan spectra yields uncertainties at the $\lesssim 10$ mK level. This significant improvement only appears if the spectra are modeled with one set of basis vectors, instead of using multiple sets of basis vectors that independently model each spectrum. Assuming that they are simulated accurately, measuring all four Stokes parameters also leads to lower uncertainties. These two strategies can be employed simultaneously and fitting multiple time bins of all four Stokes parameters yields the best precision measurements of the 21-cm signal, approaching the noise level in the data.

We present the completion of a data analysis pipeline that self-consistently separates global 21-cm signals from large systematics using a pattern recognition technique. In the first paper of this series, we obtain optimal basis vectors from signal and foreground training sets to linearly fit both components with the minimal number of terms that best extracts the signal given its overlap with the foreground. In this second paper, we utilize the spectral constraints derived in the first paper to calculate the full posterior probability distribution of any signal parameter space of choice. The spectral fit provides the starting point for a Markov Chain Monte Carlo (MCMC) engine that samples the signal without traversing the foreground parameter space. At each MCMC step, we marginalize over the weights of all linear foreground modes and suppress those with unimportant variations by applying priors gleaned from the training set. This method drastically reduces the number of MCMC parameters, augmenting the efficiency of exploration, circumvents the need for selecting a minimal number of foreground modes, and allows the complexity of the foreground model to be greatly increased to simultaneously describe many observed spectra without requiring extra MCMC parameters. Using two nonlinear signal models, one based on EDGES observations and the other on phenomenological frequencies and temperatures of theoretically expected extrema, we demonstrate the success of this methodology by recovering the input parameters from multiple randomly simulated signals at low radio frequencies (10-200 MHz), while rigorously accounting for realistically modeled beam-weighted foregrounds.

Detecting the cosmological sky-averaged (global) 21 cm signal as a function of observed frequency will provide a powerful tool to study the ionization and thermal history of the intergalactic medium (IGM) in the early Universe ($\sim$ 400 million years after the Big Bang). The greatest challenge in conventional total-power global 21 cm experiments is the removal of the foreground synchrotron emission ($\sim 10^3$-$10^4$ K) to uncover the weak cosmological signal (tens to hundreds of mK), especially since the intrinsic smoothness of the foreground spectrum is corrupted by instrumental effects. Although the EDGES team has recently reported an absorption profile at 78 MHz in the sky-averaged spectrum, it is necessary to confirm this detection with an independent approach. The projection effect from observing anisotropic foreground source emission with a wide-view antenna pointing at the North Celestial Pole (NCP) can induce a net polarization, referred as the Projection-Induced Polarization Effect (PIPE). Due to Earth's rotation, observation centered at the circumpolar region will impose a dynamic sky modulation on the net polarization's waveforms which is unique to the foreground component. In this study, we review the implementation practicality and underlying instrumental effects of this new polarimetry-based technique with detailed numerical simulation and a testbed instrument, the Cosmic Twilight Polarimeter (CTP). In addition, we explore an SVD-based analysis approach for separating the foreground and instrumental effects from the background global 21 cm signal using the sky-modulated PIPE.

Most of the current Global 21-cm experiments include ground screens that help moderate effects from the Earth. In this paper, we report on a possible systematic artifact within the ground plane that may produce broad absorption features in the spectra observed by these experiments. Using analytical approximations and numerical modeling, the origin of the artifact and its impact on the sky-averaged spectrum are described. The publicly released EDGES dataset, from which a 78 MHz absorption feature was recently suggested, is used to probe for the potential presence of ground plane resonances. While the lack of a noise level for the EDGES spectrum makes traditional goodness-of-fit statistics unattainable, the rms residual can be used to assess the relative goodness of fits performed under similar circumstances. The fit to the EDGES spectrum using a model with a simple 2-term foreground and three cavity-mode resonances is compared to a fit to the same spectrum with a model used by the EDGES team consisting of a 5-term foreground and a flattened Gaussian signal. The fits with the physically motivated resonance and empirical flattened Gaussian models have rms residuals of 20.8 mK (11 parameters) and 24.5 mK (9 parameters), respectively, allowing us to conclude that ground plane resonances constitute another plausible explanation for the EDGES data.

The Dark Ages, probed by the redshifted 21-cm signal, is the ideal epoch for a new rigorous test of the standard LCDM cosmological model. Divergences from that model would indicate new physics, such as dark matter decay (heating) or baryonic cooling beyond that expected from adiabatic expansion of the Universe. In the early Universe, most of the baryonic matter was in the form of neutral hydrogen (HI), detectable via its ground state's spin-flip transition. A measurement of the redshifted 21-cm spectrum maps the history of the HI gas through the Dark Ages and Cosmic Dawn and up to the Epoch of Reionization (EoR). The Experiment to Detect the Global EoR Signature (EDGES) recently reported an absorption trough at 78 MHz (redshift z of 17), similar in frequency to expectations for Cosmic Dawn, but about 3 times deeper than was thought possible from standard cosmology and adiabatic cooling of HI. Interactions between baryons and slightly-charged dark matter particles with electron-like mass provide a potential explanation of this difference but other cooling mechanisms are also being investigated to explain these results. The Cosmic Dawn trough is affected by cosmology and the complex astrophysical history of the first luminous objects. Another trough is expected during the Dark Ages, prior to the formation of the first stars and thus determined entirely by cosmological phenomena (including dark matter). Observations on or in orbit above the Moon's farside can investigate this pristine epoch (15-40 MHz; z=100-35), which is inaccessible from Earth. A single cross-dipole antenna or compact array can measure the amplitude of the 21-cm spectrum to the level required to distinguish (at >5$\sigma$}) the standard cosmological model from that of additional cooling derived from current EDGES results. This observation constitutes a powerful, clean probe of exotic physics in the Dark Ages.

The reduced chi-squared statistic is a commonly used goodness-of-fit measure, but it cannot easily detect features near the noise level, even when a large amount of data is available. In this paper, we introduce a new goodness-of-fit measure that we name the reduced psi-squared statistic. It probes the two-point correlations in the residuals of a fit, whereas chi-squared accounts for only the absolute values of each residual point, not considering the relationship between these points. The new statistic maintains sensitivity to individual outliers, but is superior to chi-squared in detecting wide, low level features in the presence of a large number of noisy data points. After presenting this new statistic, we show an instance of its use in the context of analyzing radio spectroscopic data for 21-cm cosmology experiments. We perform fits to simulated data with four components: foreground emission, the global 21-cm signal, an idealized instrument systematic, and noise. This example is particularly timely given the ongoing efforts to confirm the first observational result for this signal, where this work found its original motivation. In addition, we release a Python script dubbed $\texttt{psipy}$ which allows for quick, efficient calculation of the reduced psi-squared statistic on arbitrary data arrays, to be applied in any field of study.

The sky-averaged (global) highly redshifted 21-cm spectrum from neutral hydrogen is expected to appear in the VHF range of $\sim20-200$ MHz and its spectral shape and strength are determined by the heating properties of the first stars and black holes, by the nature and duration of reionization, and by the presence or absence of exotic physics. Measurements of the global signal would therefore provide us with a wealth of astrophysical and cosmological knowledge. However, the signal has not yet been detected because it must be seen through strong foregrounds weighted by a large beam, instrumental calibration errors, and ionospheric, ground and radio-frequency-interference effects, which we collectively refer to as "systematics". Here, we present a signal extraction method for global signal experiments which uses Singular Value Decomposition (SVD) of "training sets" to produce systematics basis functions specifically suited to each observation. Instead of requiring precise absolute knowledge of the systematics, our method effectively requires precise knowledge of how the systematics can vary. After calculating eigenmodes for the signal and systematics, we perform a weighted least square fit of the corresponding coefficients and select the number of modes to include by minimizing an information criterion. We compare the performance of the signal extraction when minimizing various information criteria and find that minimizing the Deviance Information Criterion (DIC) most consistently yields unbiased fits. The methods used here are built into our widely applicable, publicly available Python package, $\texttt{pylinex}$, which analytically calculates constraints on signals and systematics from given data, errors, and training sets.

The redshifted 21-cm monopole is expected to be a powerful probe of the epoch of the first stars and galaxies ($10<z<35$). The global 21-cm signal is sensitive to the thermal and ionization state of hydrogen gas and thus provides a tracer of sources of energetic photons -- primarily hot stars and accreting black holes -- which ionize and heat the high redshift intergalactic medium (IGM). This paper presents a strategy for observations of the global spectrum with a realizable instrument placed in a low altitude lunar orbit, performing night-time 40-120 MHz spectral observations, while on the farside to avoid terrestrial radio frequency interference, ionospheric corruption, and solar radio emissions. The frequency structure, uniformity over large scales, and unpolarized state of the redshifted 21-cm spectrum are distinct from the spectrally featureless, spatially-varying, and polarized emission from the bright foregrounds. This allows a clean separation between the primordial signal and foregrounds. For signal extraction, we model the foreground, instrument, and 21-cm spectrum with eigenmodes calculated via Singular Value Decomposition analyses. Using a Markov Chain Monte Carlo algorithm to explore the parameter space defined by the coefficients associated with these modes, we illustrate how the spectrum can be measured and how astrophysical parameters (e.g., IGM properties, first star characteristics) can be constrained in the presence of foregrounds using the Dark Ages Radio Explorer (DARE).