The 21-cm bispectrum as a probe of non-Gaussianities due to X-ray heating

We present analysis of the normalized 21-cm bispectrum from fully-numerical simulations of intergalactic-medium heating by stellar sources and high-mass X-ray binaries (HMXBs) during the cosmic dawn. Ly-alpha coupling is assumed to be saturated, we therefore probe the nature of non-Gaussianities produced by X-ray heating processes. We find the evolution of the normalized bispectrum to be very different from that of the power spectrum. It exhibits a turnover whose peak moves from large to small scales with decreasing redshift, and corresponds to the typical separation of emission regions. This characteristic scale reduces as more and more regions move into emission with time. Ultimately, small-scale fluctuations within heated regions come to dominate the normalized bispectrum, which at the end of the simulation is almost entirely driven by fluctuations in the density field. To establish how generic the qualitative evolution of the normalized bispectrum we see in the stellar + HMXB simulation is, we examine several other simulations - two fully numerical simulations that include quasi-stellar object (QSO) sources, and two with contrasting source properties produced with the semi-numerical simulation 21CMFAST. We find the qualitative evolution of the normalized bispectrum during X-ray heating to be generic, unless the sources of X-rays are, as with QSOs, less numerous and so exhibit more distinct isolated heated profiles. Assuming mitigation of foreground and instrumental effects are ultimately effective, we find that we should be sensitive to the normalized bispectrum during the epoch of heating, so long as the spin temperature has not saturated by z approximate to 19.