The Anomalous 21-cm Absorption at High Redshifts

The EDGES collaboration has reported the detection of a global 21-cm signal with a plateau centered at 76 MHz (i.e., redshift 17.2), with an amplitude of 500^(+200)_(-500) mK. This anomalous measurement does not comport with standard cosmology, which can only accommodate an amplitude<230 mK. Nevertheless, the line profile's redshift range (15<z<20) suggests a possible link to Pop III star formation and an implied evolution out of the `dark ages.' Given this tension with the standard model, we here examine whether the observed 21-cm signal is instead consistent with the results of recent modeling based on the alternative Friedmann-Lemaitre-Robertson-Walker cosmology known as the R_h=ct universe, showing that--in this model--the CMB radiation might have been rethermalized by dust ejected into the IGM by the first-generation stars at redshift z<16. We find that the requirements for this process to have occurred would have self-consistently established an equilibrium spin temperature T_s~3.4 K in the neutral hydrogen, via the irradiation of the IGM by deep penetrating X-rays emitted at the termination shocks of Pop III supernova remnants. Such a dust scenario has been strongly ruled out for the standard model, so the spin temperature (~3.3 K) inferred from the 21-cm absorption feature appears to be much more consistent with the R_h=ct profile than that implied by LCDM, for which adiabatic cooling would have established a spin temperature T_s(z=17.2)~6 K.


Introduction
The EDGES collaboration has recently reported the detection of a 21-cm signal in absorption between redshifts 20 and 15, with an amplitude of 500 mK, roughly twice the strength expected from cosmological simulations [1]. The most likely origin of this signal is the absorption of CMB photons by neutral atomic hydrogen and, given that it arises from the redshift interval in which the first generation of (Pop III) stars is believed to have formed, may have important implications for the early phase of cosmic structure formation.
A complication with the analysis of this signal, however, is the challenge of making this measurement, given the very large foreground of Galactic diffuse synchrotron emission. Though this foreground is spectrally smooth above 50 MHz, it unfortunately requires several components for a proper modeling [2]. The reality of the anomalous 21-cm absorption feature associated with the cosmic microwave background (CMB) has therefore been called into question. Nevertheless, Bowman et al. (2018) performed numerous hardware and processing tests to validate their detection, demonstrating that no other astronomical or atmospheric mechanism could have produced the observed profile. Until future observations, e.g., with the Square Kilometre Array (SKA; ref. [3]), can confirm or reject this finding, the implications of such a global a John Woodruff Simpson Fellow.
high-redshift 21-cm absorption line on cosmological models remain somewhat ambiguous and warrant further investigation. In this paper, we shall advocate for the reality of this anomalous signal and adopt the following essential characteristics: the 21-cm line extends over the redshift range 15 z 20, with an amplitude of 500 +200 −500 mK at the 99% confidence level. Its plateau is centered at a frequency of 78 MHz, corresponding to a redshift of 17.2. The depth of this feature is effectively 3.8σ away from the prediction of standard cosmology (ΛCDM) which, as we shall discuss shortly, cannot accommodate an amplitude greater than ≈ 230 mK.
Given the disparity between the measured and predicted profile of the global 21-cm line, our principal focus in this paper will be to gauge whether this anomalous measurement might instead offer some observational support for the alternative Friedmann-Lemaître-Robertson-Walker (FLRW) cosmology known as the R h = ct universe [4,5], which has much in common with ΛCDM, differing with it principally via an essential constraint from general relativity referred to as the "zero active mass condition," i.e., an equation-of-state ρ + 3p = 0 in the cosmic fluid, in terms of the total energy density, ρ, and pressure, p [6].
As we shall demonstrate in this paper, the gas temperature at z = 17.2 required to produce the observed 21-cm signal is consistent with the physical conditions of the intergalactic medium (IGM) expected during this epoch, if the expansion history of the Universe reflects the dynamics implied by this equation-of-state. In § 2, we shall summarize the status of this recent work and then, in § 3, describe the stellar and IGM environment expected in this model. We end with a discussion of our results and a conclusion in § 4.

The R h = ct Universe
This FLRW cosmology has been under development for over a decade, and its predictions have been tested using a broad range of data, incorporating both integrated and differential signatures at all redshifts (for some recent developmental papers, see refs. [7,8,6]; for a summary of the comparative tests, see Table 2 in ref. [9]; a complete description of this model is presented in ref. [5]). The zero active mass condition implies an equation-of-state p = wρ, with w = (ρ r /3 + w de ρ de )/ρ = −1/3 at all epochs, under the assumption that the pressure of matter is always negligible. Throughout this paper, ρ r , ρ m and ρ de denote, respectively, the radiation, matter and dark-energy densities, and p de = w de ρ de . We shall write ρ Λ in the special case where dark-energy is assumed to be a cosmological constant.
In this picture, the luminosity distance d L (z) and Hubble rate of expansion H(z) have very simple forms, dependent only on the parameter H 0 : and The R h = ct cosmology is flat (i.e., k = 0), for otherwise the gravitational radius would not equal ct [10]. Therefore, from the first Friedmann equation [4], one easily finds that where is the critical energy density. For simplicity, we scale all relevant quantities to the Planck value of H 0 though, in reality, the actual H 0 in R h = ct is somewhat different. Even so, it would differ from this value by no more than ∼ 10 − 15% (see, e.g., ref. [11]).
Our imposition of the zero active mass constraint on the cosmic fluid implies that the fractional representation of the various species is modified from the evolution seen in ΛCDM. As we shall discuss further below, this is partly the reason why the rate of structure formation in this model, notably the halo growth rate and stellar formation rate, differ considerably from their counterparts in the standard model. Let us define the dark energy density as ρ de = ̟(z) ρ. The equation-of-state in R h = ct does not tell us the value of ̟, but we know from observations of the local Universe that ∼ 30% of ρ at low redshifts is comprised of (luminous and dark) matter. Thus, if ̟ = 2/3 at z → 0, the expression for p implies that w de = −1/2, precluding any possibility of dark energy being a cosmological constant. We do not know whether this equation-of-state changes with redshift, but if we assume for simpicity that w de is constant, then the early Universe must have both radiation and dark energy, in proportions such that ̟ = 4/5 at z >> 1. Dark energy is therefore always present in R h = ct, though its fractional contribution to ρ evolves slightly from 0.8 at early times to 2/3 locally. In concert with this change, radiation is gradually replaced by matter, with and If we further reasonably assume that the radiation is always a blackbody, the temperature evolves according to the standard free-streaming relation after decoupling (at z cmb ), but follows the modified form at high redshifts.
In recent work [12], we discussed the reasons (both theoretical and observational) why z cmb in this model cannot be as large (i.e., ∼ 1080) as in the standard model. Indeed, we may see right away from Equations (7) and (8) that this transition must occur, crudely, at the redshift where these two expressions overlap, i.e., at z cmb 42. More quantitatively, Planck [13] has identified an angular scale θ s = (0.596724 ± 0.00038) • on the last scattering surface (LSS), interpreted as the acoustic horizon-the maximum distance traveled by sound waves in the early universe (see, e.g., refs. [14,15,16,17]). The comoving distance r s corresponding to this scale is thought to have remained constant thereafter.
A characteristic scale has also been seen in the correlation function of galaxies and the Ly-α forest [18,17,19,20,21,22,23,24,25,26]. Interpreting the peak seen in large galaxy surveys as the subsequent manifestation of the acoustic horizon, the baryon acoustic oscillation (BAO) comoving scale in the context of R h = ct has been determined to have the value r BAO = 131.1 ± 4.3 Mpc [27,9], and it is straightforward to show from Equation (1) that where R h (t 0 ) = c/H 0 is the gravitational (or Hubble) radius today. Taking the measurement errors into account, we therefore infer that (in the R h = ct cosmology), corresponding to a cosmic time t cmb ≈ 849 Myr. Matter could therefore not have been ionized at this redshift, so the recombination picture for the origin of the CMB could not work in this model. A recent analysis [12] of the dust injected into the IGM by Pop III stars instead suggests that the primordial radiation field would have been rethermalized by grains at z cmb ∼ 16 to produce the CMB we see today. In this picture, the CMB anisotropies therefore reflect the large-scale structure seen just prior to the onset of reionization at z ∼ 15. Additional evidence for the value of z cmb in Equation (10) is provided by the CMB angular correlation function [28]. The lack of large-angle correlations in the CMB temperature fluctuations conflicts with the predictions of inflationary ΛCDM at a high level of confidence ( 3σ; ref. [29]). Following a careful analysis of the latest Planck release, we now understand that the weakness of the large-angle correlations is probably due to a non-zero minimum wavenumber k min in the primordial power spectrum P (k) [28]. Its relevance to the subject of this paper is that, if real, the inferred value of k min , in the context of R h = ct, corresponds to a redshift z cmb = 17.05 +8 −5 -a compelling confirmation of the estimate shown in Equation (10), based on an entirely different argument.
But is it really possible to even contemplate a dusty origin for the CMB, given that all of the evidence today appears to overwhelmingly favour recombination at z ∼ 1080? In the next section, we shall summarize the extensive analysis carried out to address this far-fetched idea. As it turns out, such a model was actually seriously considered in the late 1900's before its evident flaws made it redundant. But we now know that, while a dust model for the CMB is untenable in ΛCDM, it may still be viable in other cosmologies. As it turns out, dust could still play a critical role in forming the microwave background if R h = ct is the correct model. And as we shall see, this scenario also self-consistently accounts for the 21-cm signal measured by the EDGES collaboration, which (subject to confirmation by future measurements) is otherwise not consistent with the standard model.

Reassessing Dust's Role in Forming the CMB
Placing the last scattering surface (LSS) at z cmb ∼ 16 may seem to be at odds with many kinds of observation, but this is only true in the case of ΛCDM. As shown in ref. [12], there are at least three principal observational signatures one may use to distinguish a CMB originating via recombination at z ∼ 1080 from one due to reprocessing by dust at z ∼ 16. These are: the spectrum itself, which appears to be a true blackbody [30]; the presence or absence of recombination lines [31,32], and a test of whether the angular power spectrum varies with frequency [13].
The physical attributes required of an LSS at z ∼ 16 in the context of R h = ct ironically echo some of the theoretical ideas explored for a dusty origin of the CMB several decades ago [33,34,35]. Before attempting to rescue this now defunct model in order to explain the anomalous 21-cm signal, however, it is essential to consider whether such a proposal even makes sense based on what we know today.
The first issue is quite obvious: recombination lines ought to be present at some level in the CMB spectrum if the standard picture is correct, whereas dust rethermalization at z ∼ 16 would have wiped all of them out. Extensive simulations have already been carried out to predict the level at which we should see such lines in ΛCDM [31,32], but unfortunately the effect of recombination lines on the angular power spectrum is expected to be quite small. Though it may eventually be separated from other effects with the improved sensitivities of future experiments, there is no evidence today of recombination lines in the CMB, so this test must await future developments.
Arguably the most convincing observation made to support the recombination scenario was COBE's discovery [30] that the CMB's spectrum is a near perfect blackbody. We shall not repeat the argument made in ref. [12], but it is not difficult to show that optically thick dust in thermal equilibrium with the radiation it rethermalizes near the photosphere (at the LSS) also produces a near perfect blackbody field, as was already suggested by some of the earlier work [34]. Thus, the key issue is not that dust opacity is frequency dependent but, rather, that dust would need to reach local thermal equilibrium with the radiation. This question hinges on how much dust was produced by Pop III stars in the redshift range 15 z 20, and we shall discuss this in § 3.1 below.
The third issue is whether the angular power spectrum of the CMB is frequency-dependent, as one might expect if the anisotropies in the temperature distribution across the sky vary among surveys conducted at different wavelengths. One would not expect photospheric optical depths to affect the observed distribution of fluctuations with Thomson scattering because the optical depth it produces is independent of frequency. Maps made at different wavelengths should therefore be identical in ΛCDM. But this would not be the case if the opacity were frequency dependent, as would happen in the case of dust. Though photospheric depth effects might not significantly change the larger fluctuations from one map to another, they could still alter the observed anisotropy distribution on smaller scales, which would in turn produce variations in the inferred CMB power spectrum constructed at different wavelengths. A careful analysis of this frequency dependence [36] shows that indeed the location of the LSS changes somewhat with frequency. Variations can be as high as ∼ 2% at ℓ ∼ 400, increasing to ∼ 5% for ℓ 800. If real, these changes would argue against a Thomson scattering opacity, but would be consistent with the dust model in R h = ct (see ref. [12]) as long as the location of the LSS were restricted to the redshift range ∼ 14 − 16. This period would have lasted ∼ 100 Myr, after which the Pop III and the early Pop II supernovae would have completely destroyed the dust, making the IGM transparent again and initiating the epoch of reionization at z ∼ 15. Under these constraints, a frequency-dependent opacity would have impacted the CMB spectrum by no more than a few percent, consistent with the current Planck observational measurements.

Pop III Stars in the R h = ct Universe
Assessing the CMB-imposed constraints on the Pop III star formation rate and evolution is critical to understanding whether the IGM at 15 z 20 could also have produced the observed global 21-cm signal. Again, we won't repeat the detailed analysis carried out in ref. [12], but we shall here merely summarize the key findings from that work. In order for every photon in the CMB to have been absorbed by dust prior to z ∼ 16, we would require the dust number density to have been n s (z = 16) ∼ 5 × 10 −15 f Z cm −3 for a bulk density of ∼ 2 g cm −3 of silicate grains, and a grain radius r s ∼ 0.1 micron. This restricts the IGM metallicity f Z near the end of the Pop III star formation and evolution era to ∼ 20% of the solar value, comfortably small to avoid any obvious inconsistency with the general perception that the bulk of today's abundances grew across subsequent phases of star formation.
In order for the dust temperature to have remained in equilibrium with the CMB during that period, two important factors had to be satisfied. The first has to do with the average heating and cooling rates for a given dust particle, while the second is based on the fact that each absorption of a photon produces a quantum change in the dust particle's temperature, and is thus strongly dependent on its size [37,38]. Under the conditions expected in the IGM at that redshift, it would have taken a ∼ 0.1µm sized particle roughly 50 seconds to reach equilibrium at a dust temperature T d ∼ 46 K, so the first condition would have been satisfied trivially. The second issue is more constraining. The assumption of a smooth evolution in T d starts to break down for grains smaller than r s ∼ 0.003 µm [38], at which point the heating starts to produce temperature spikes. Putting these estimates together, we see that the dust model required for consistency with the observed CMB spectrum would therefore be based on silicates with size ∼ 0.003 − 0.3 µm, perhaps even larger, though sizes much larger than ∼ 0.3 µm would violate our earlier estimate of n s (z = 16) and the reasonable value of the metallicity f Z ∼ 0.2.
Let us now see what these constraints have to say about the Pop III star formation rate and evolution. Believed to have reached masses 500 M ⊙ M 21 M ⊙ [39,40], these stars emitted copious high-energy radiation that ionized the halos within which they were born. At the end of their brief (∼ 10 6 −10 7 yr) lives, most of them exploded as SNe [41], ejecting heavy elements into the IGM [42], also producing X-rays via shock acceleration that penetrated deeply into the IGM (more on this below). Given the above estimates for the dust size and number density, one can infer that roughly 9 × 10 44 g Mpc −3 (co-moving volume) of dust material had to be injected into the IGM during this epoch (20 z 15). Thus, adopting a typical mass M ∼ 100 M ⊙ , with a typical ejection fraction of 30% [43], we see that ∼ 1.5 × 10 8 Mpc −3 Pop III stars must have exploded as SNe between z = 20 and 15 to provide the dust required to completely thermalize with the CMB radiation. As we shall see below, this star-formation rate (SFR) is critical to understanding how the global 21cm signal could have been produced within the context of R h = ct.
Those familiar with current models of star formation during this epoch (see, e.g., refs. [44,45,46,47,48,49] would realize that such an aggressive SFR would be orders of magnitude greater than what is expected in the context of ΛCDM. Comparable simulations have not yet been completed using R h = ct as the background cosmology, but already we have firm indications that large-scale structure (LSS) formation proceeded at a very different rate in this model, due to several influential factors, including: (i) the seeding of a primordial fluctuation spectrum [50]; (ii) the alternative redshift dependence of the relative abundances (Eqns. 5, 6); (iii) the different timeline versus redshift [51]; (iv) the modified growth equation, which has now been solved semi-analytically [52,53,54]. What emerges from a comparison of LSS formation in these two cosmologies is a significant difference in the growth rate of stellar mass at z 4, reaching values exceeding ∼ 10 4 at z 10. The observational evidence one may use to distinguish between these two scenarios is still controversial and strongly debated, but a quick inspection of figures 2-7 in ref. [53] would reveal that the data do not yet rule out the higher rates predicted by R h = ct. Indeed, the currently inferred halo mass function appears to differ from what is expected at z ∼ 10 in ΛCDM by over 4 orders of magnitude [55]. Some additional support for the timeline in Rh = ct is also provided by its consistency with the early appearance of supermassive black holes [56], which would otherwise pose quite a challenge to the time versus redshift relation in the standard model.
In addition to the dust they subsequently ejected into the IGM, Pop III stars also emitted copious amounts of UV light during their lives, and X-rays via their supernova remnants. A typical Pop III star with mass M ∼ 100 M ⊙ was a blackbody emitter with radius R * = 3.9 R ⊙ and surface effective temperature T * = 10 5 K. It is straightforward to estimate from this that, during their evolution, they would have bathed the IGM with an energy density U UV ∼ 4 × 10 63 erg Mpc −3 . Compared to the energy density in the CMB at that time, this would have amounted to no more than ∼ 0.5% of the total radiation field, but these UV photons would have had important consequences concerning the HI gas and spin temperature we shall be discussing shortly.

The Gas and Spin Temperature at
The principal effort in modeling the global 21-cm signal is directed towards understanding the gas kinetic (T g ) and HI spin temperature (T s ) during the epoch (z ∼ 17.2) when the absorption line formed. The 21-cm line is the triplet-to-singlet hyperfine transition of the atomic hydrogen ground state due to the coupling of the magnetic moments of the proton and electron (for a review, see, e.g., ref. [57]). For a system in thermal equilibrium, the relative populations of the two spin levels are given by the ratio where ∆E = 5.9×10 −6 eV is the energy difference between the two distinct spin states, and g 1 /g 0 = 3 is the ratio of the statistical degeneracy factors of the two levels.
The spin temperature coupled to either the gas or CMB temperature, depending on which physical processes dominated the excitation of the line at 15 z 20. There appear to be three dominant effects: (1) the absorption of CMB photons which, in the absence of other interactions, would imply that T s = T γ ; (2) nucleon-nucleon collisions, which would couple T s to the gas temperature T g ; and (3) excitation due to the so-called Wouthuysen-Field effect [58]. In this process, Ly-α photons produced, say, by the Pop III stars, excited the ground state of atomic hydrogen to the 2P level, which then re-emited Ly-α photons and entered either of the two spin states, creating an asymmetric redistribution of the electrons between the hyperfine levels. In equilibrium [59], one has In this expression, x c is the collisional coupling coefficient for H-H and H-e −1 scatterings [60]. The middle term accounts for the Wouthuysen-Field effect, with a coupling coefficient [59] x α = 1.81 where S α is a factor of order unity describing the detailed atomic physics of the scattering mechanism, and J α is the Ly-α background (number) intensity in units of cm −2 s −1 Hz −1 sr −1 . This process couples T s to an effective color temperature T eff which, however, never deviates from T g by more than ∼ 20% (see fig. 4 in ref. [59]), due to the recoils induced by repeated scatterings. It is straightforward to estimate J α from the total Pop III emissivity discussed in the previous subsection. One finds for the redshift range 15 z 20 in R h = ct that This value does not include the possible contribution from nascent quasars, which might have started forming during this epoch as well, so it should be viewed as a lower limit. Thus, the coefficient x α in Equation (12) would have been 2, 400, and consequently the spin temperature T s would have been tightly coupled to the gas temperature throughout the period where the global 21-cm absorption line was formed.
The 21-cm brightness temperature (relative to the CMB) averaged over the whole sky is given by ref. [61] δT b = 27 mK x HI (1 + δ) where x HI is the fraction of neutral hydrogen (≈ 1 in this regime), δ is the fractional overdensity (which is smaller than 1, so that 1+δ ≈ 1), and Ω m and Ω b are, respectively, the matter and baryonic densities today in units of the critical density. The quantity ∂ r v r is the velocity gradient along the line of sight and is approximately equal to (1 + z)H(z) within this redshift range. The application of Equation (15) to the observed global 21-cm signal thus implies that T s (and therefore T g ) must have been approximately 3.3 K during this epoch (remembering that T b is actually negative), compared to the CMB temperature T γ = 2.725(1 + 17.2) ≈ 49.6 K. So the key question is now "How did the gas reach this equilibrium at 15 z 20?" Let us first consider the situation in the standard model. The conventional assumption in ΛCDM is that, while the radiation effectively decoupled from the baryons at z ∼ 1080 (LSS), the particles were much less numerous and therefore remained thermalized with the radiation for much longer [62]. The rate of energy transfer per unit volume between the CMB photons and free electrons may be written [63] where σ T is the Thomson cross section, a r is the radiation constant, k B is the Boltzmann constant, m e is the electron mass, T γ is the photon temperature, T e is the electron temperature, n H is the nucleon number density and x e is the free electron fraction. One further assumes that the electrons preserved thermal contact with the nucleons, so that T e was effectively the gas temperature T g . Using this expression, one may find the redshift z ad at which thermal contact between the matter and the CMB was broken, which occured when the ratio of the heat content of the matter to the heat-transfer rate in Equation (16) exceeded the characteristic expansion time of the Universe. Beyond this time, matter continued to cool adiabatically, T g ∝ (1+ z) 2 , at a rate faster than radiation. In ΛCDM, one finds that z ad ≈ 150, and so one expects the gas temperature at z = 17.2 to have been But as we have just seen, in order to account for the observed global 21-cm signal, the gas temperature could not have been higher than ≈ 3.3 K. Thus, to correctly account for the observed profile of the 21-cm absorption line, adiabatic cooling in the standard model would have to have started much earlier, at z ad ∼ 250. This is the reason the 21-cm signal measured by EDGES is considered to be an anomaly.
In the R h = ct universe, the temperature at high redshifts (Eq. 8) is lower than that in the standard model, so adiabatic cooling of the gas would have started even earlier than in ΛCDM, due to the much smaller free electron fraction x e . Applying the same matter-radiation decoupling algorithm described above, one finds that z ad ≫ 150 in this model, producing a much colder gas than implied by Equation (17). This too is inconsistent with the 21-cm line profile, unless some other mechanism intervened to reheat the gas at 15 z 20. This is where the Pop III stars we discussed in the previous subsection enter the picture. We pointed out that these early stellar sources not only injected large quantities of dust and UV photons (which were critical in coupling the spin temperature to the gas temperature via Eqn. 12) into the IGM, but also irradiated the cosmic background with an X-ray intensity emitted at the termination shocks of their supernova remnants. As we pointed out, however, an important caveat with this analysis is that the Pop III star formation rate required for this to work would have greatly exceeded the predictions in ΛCDM. It would have been much more in line with the expansion scenario expected in R h = ct, though the question still remains open regarding which of them is favoured by the data (see, e.g., ref. [55]).
To produce the amount of dust required to rethermalize the CMB radiation at z ∼ 16 (see Eqn. 10), we found in § 3.1 that n III ≈ 1.5 × 10 8 Mpc −3 (18) Pop III stars of average mass M III = 100 M ⊙ must have terminated their lives as supernova explosions. Such stars reach this terminal state as pair-instability explosions [43], releasing E III 10 51 erg (19) of energy per event. Within the ejecta, a fraction ǫ e of the total supernova energy is converted into accelerated electrons. The efficiency of this process depends on the density of the medium surrounding the supernova. For low densities, the shock is weak and the efficiency is low. Supernova remnants impacting high-density clouds in the local Universe can have much larger efficiencies. We still do not have a precise understanding of how the shocks would have evolved in the early IGM, though simulations suggest that the efficiency might have been near the bottom of the range [64], i.e., ǫ e ∼ 0.01 .
The relativistic electrons produced in this fashion radiated a fraction f X of their energy as X-rays, the dominant radiation field responsible for heating the IGM. This fraction is also subject to uncertainties, e.g., the strength of the magnetic field in the radiation zone but, based on the modeling of supernova remnants in the local Universe, one would expect f X ∼ 0.1 .
Some attention has already been paid to the fate of this X-radiation in the high-redshift IGM, once the first generation of Pop III stars brought an end to the dark ages. Prior to the reionization of HI gas, UV photons would have been mostly trapped near their sources, creating localized HII bubbles, but X-rays would have penetrated to much larger distances. They therefore most likely constituted the high-energy radiative background throughout the IGM, initiating the ionization of the neutral gas far from the stars (see, e.g., ref. [65]. Most of their energy would have been deposited indirectly, starting with photoionization and subsequently dispersed by the fast electrons released in this process. It is believed that the electrons energized by X-rays from the Pop III supernovae [66,67], were responsible for heating the IGM before reionization began in earnest [68,69]. Some have even speculated-well before the EDGES measurement was known-that this heating and ionization could have important observational consequences on the highly redshifted 21-cm signal produced by the early IGM [69,70,71].
In their detailed simulations of the energy deposited into the IGM from the photoionization of neutral atoms by X-rays, ref. [65] found that ejected electrons with energy < 10.2 eV could not interact with any atoms or ions and thus dispersed all of their energy as heat. As the photon energy increased, however, approaching the X-ray region of the spectrum, more and more excitation and ionization processes became available and the fraction of energy deposited as heat decreased. But this trend came to an end once the number of such process had been saturated, and the fraction of incoming energy deposited as heat eventually approached a reasonably constant value at ∼ 1 − 10 keV. These authors found that in highly neutral gas, with n HII /n HI < 10 −3 , approximately 20% of the electron's energy was deposited as heat, with the remainder split roughly equally between ionization and excitation. They also found that these results are insensitive to n HII /n HI at high energies, because most of the heating comes from secondary electrons with energies below 10 eV. Since, in addition, only a small fraction of the X-ray photon's energy is used to photoionize hydrogen in this regime, we may simply write for the fraction of X-ray energy deposited as heat throughout the IGM during the time (i.e., 15 z 20) that Pop III stars ended their lives as supernova explosions.
Thus, all told, the Pop III X-ray energy converted into IGM heat during the 15 z 20 epoch may be estimated as follows: As we have seen, adiabatic cooling of the gas in R h = ct would have resulted in an IGM at z ∼ 20 much colder than T γ , so it is safe to assume that T g at z ∼ 17.2 would have been entirely due to the X-ray re-heating implied by Equation (23). And if we were to put (given that both the electron and n I densities were negligible compared to n H ), we would find that One would like to think that the remarkable consistency between the spin temperature (T s ≈ 3.3 K) measured by the EDGES collaboration and the gas temperature we have just estimated (T g ≈ 3.4 K) is strong evidence in favour of this model. Of course, this is almost certainly partially a coincidence, given that every parameter in Equation (25) contributing to this outcome was assigned its fiducial value. On the other hand, the fact that this analysis produced this level of consistency without any optimization of the physical inputs does suggest that the scenario we have explored here is at least viable. None of the estimates had to be 'stretched' to unlikely values, and all were internally consistent with the basic premise that the CMB we see was produced by dust rethermalization of the primordial radiation field during the Pop III star-formation era. The principal caveat with this analysis concerns the Pop III star formation rate, which differs by several orders of magnitude between the ΛCDM and R h = ct cosmologies. The picture we have painted in this paper holds together self-consistently as long as the R h = ct rate is correct, but would break down if future observations reveal that the large-scale structure formation rate at z 10 was more in line with the predictions of the standard model.
Once the global 21-cm signal is confirmed by upcoming observations, a more detailed simulation than that attempted here ought to be carried out, highlighting the evolution of the physical state of the IGM from z ∼ 20 to z ∼ 15. The EDGES absorption line was seen throughout this region, so a natural question to address is whether the brightness temperature calculated with Equation (15) is consistent with the radiation (T γ ) and spin (T s ) temperatures we have estimated thus far. The answer appears to be yes, given that T γ changes slowly at z 16, evolving from the free-streaming relation in Equation (7) to the expression in Equation (8). So, for example, while T γ ≈ 49. the IGM is truly consistent with the observed 21-cm line profile.
Nevertheless, even at this level of approximation, the brightness temperature of the 21-cm signal predicted by the R h = ct universe appears to be more consistent with the EDGES observation than that calculated using the standard model (see fig. 1). As noted earlier, the limitation with our simplified approach in this paper is that the physical conditions within the redshift range 15 z 20 are estimated more or less as averages, so we are not yet in a position to show a true time-evolution across this transition region. The brightness temperature (Eqn. 15), however, does trace the evolution in CMB temperature according to Equation (7) in ΛCDM, and according to a transition from Equation (8) to (7) in R h = ct, so we may at least gauge how well the 21-cm line profile matches the CMB temperature in these two models within this redshift range. As one may see in figure 1, the correspondence is much better for R h = ct than for ΛCDM, largely due to two principal effects: (i) the difference in spin temperature predicted by these two cosmologies (T s ≈ 3.4 K in the former, versus ≈ 6 K in the latter); and (ii) the somewhat different redshift dependence of the CMB temperature at z 16, given by Equation (7) in ΛCDM and Equation (8) in R h = ct.

Discussion and Conclusion
There is no question that the debate between a recombination versus dust origin for the CMB has shifted strongly in favour of the former following the dramatic discoveries by COBE, WMAP and Planck over the past two decades. Yet this issue has slowly resurfaced in the face of growing tension between the predictions of the standard model and the ever improving precision of cosmological measurements. A good example is the Hubble constant H 0 , which characterizes the current expansion rate of the Universe and determines its absolute distance scale. The accuracy of measuring H 0 has been significantly improved recently, but its value (67.4 ± 0.5 km s −1 Mpc −1 ; ref. [13] inferred from the CMB observations in the context of flat ΛCDM disagrees with that based on local Type Ia SNe calibrated by the Cepheid distance ladder (74.03 ± 1.42 km s −1 Mpc −1 ; ref. [72]) at a 4.4σ level of significance.
A previously completed thorough analysis of the evidence in favour of recombination [12] demonstrated that a dust scenario for the origin of the CMB is still viable, though not in the context of ΛCDM. Instead, the observed characteristics of the microwave background, such as its distribution of anisotropies and the implied acoustic scale, could be consistent with the alternative FLRW cosmology known as the R h = ct universe, which is effectively ΛCDM, though with an additional constraint on its overall equation-of-state derived from the zero active mass condition in general relativity [5]. The work we have reported in this paper has a significant impact on this discussion because it represents one of the earliest tests of the dust model based on a very different kind of observation.
The observed profile of the global 21-cm signal associated with the CMB is largely consistent with the expectations of ΛCDM, except that its amplitude is more than a factor of two greater than the largest predictions. Like the disparity with H 0 , this anomaly represents a statistically significant indication that the standard model may not be quite right. It may need an infusion of new physics or, at worst, a complete replacement. The R h = ct universe is a compromise, because it actually does not replace any of the essential features of the standard model (other than the inflationary phase, which itself appears to be problematic; ref. [73]), but only adds the requirement that ρ + 3p = 0 at all times. From a theoretical standpoint, the Local Flatness Theorem in general relativity [6] leaves no doubt that the use of FLRW is valid only with the inclusion of this equation-of-state. And its observational support continues to grow with the completion of each new comparative test (see, e.g., Table 2 in ref. [9]). The fact that R h = ct presents a viable mechanism by which the measured 21-cm signal could have formed at 15 z 20 constitutes additional observational support for this model.
If confirmed, the EDGES detection [1] is critically important to astrophysics in general, perhaps even to particle physics, because it opens up a window on the early phase of cosmic structure formation, and provides us with the physical state of the cosmic environment at the start of the epoch of reionization. As we have seen in this paper, the requirements for the rethermalization of the CMB by dust at z ∼ 16 in the context of R h = ct seemlessly weave together a narrative in which the onset of Pop III star formation at 15 z 20 self-consistently initiated a reheating of the IGM, an increase in metallicity and the concomitant injection of dust into the background medium. A single population of stars produced all of the physical attributes needed to account for the CMB spectrum, the apparent frequency dependence of its distribution of anisotropies, and an equilibrium IGM gas (and spin) temperature reflected in the global 21-cm profile produced during this epoch. This represents a comprehensive set of developments fulfilling the cosmic transition away from the dark ages and into the subsequent era of reionization and large-scale formation of structure.
Fortunately, it appears that we shall not have to wait very long to see a confirmation of the EDGES measurement. Several other similar experiments are underway, including the Large-Aperture Experiment to Detect the Dark Ages (LEDA [74]; the Sonda Cosmológica de las Islas para la Detección de Hidrógeno Neutro (SCI-HI [75]); and the Shaped Antenna measurement of the background Radio Spectrum 2 (SARAS2 [76]). Farther afield, the observation of the 21-cm line should be significantly enhanced by the use of interferometric arrays, such as the Hydrogen Epoch of Reionization Array (HERA [77]); and the Square Kilometre Array (SKA; https://www.skatelescope.org), among others. When constructed, the SKA Low-Frequency Aperture Array will detect the power spectrum associated with the EDGES absorption profile, and should also be able to image the 21-cm signal, providing more fertile ground for testing the scenario we have explored in this paper, in which the physical conditions producing the 21-cm absorption line are inextricably linked to the requirements for rethermalizing the CMB at z ∼ 16, just prior to the epoch of reionization.
I am very grateful to the anonymous referee for their exceptional review of this paper and very helpful suggestions that have led to an improved presentation of the results. I am also grateful to Amherst College for its support through a John Woodruff Simpson Lectureship, and Purple Mountain Observatory in Nanjing, China, for its hospitality while part of this work was being carried out.