Epoch Of Heating Research Articles

The ionizing and heating power of ultraluminous X-ray sources under the geometrical beaming model

While there is now a consensus that X-ray binaries (XRBs) are the dominant X-ray sources in the early Universe and play a significant role during the epoch of heating of the intergalactic medium (IGM), recent studies report contradicting results regarding their contribution in the nebular emission of local Universe galaxies. Ultraluminous X-ray sources (ULXs), which dominate the X-ray budget of normal galaxies, may be important interstellar-medium (ISM) ionizing sources. However, their output in the extreme-ultraviolet (EUV) and soft X-ray part of the spectrum remains observationally unconstrained. In this paper, we predict the ionizing and heating power from ULX populations under the geometrical beaming scenario, and three models describing the emission from super-critical accretion disks. We find that our theoretical spectra for ULX populations cannot (can) explain the He II (Ne V) emission observed in some galaxies, with their contribution being less (more) important than the underlying stellar population. Stochastic fluctuations in the number of ULXs may allow for equal contributions in the He II emission, in a fraction of galaxies. We provide average spectra of ULX populations as an input to local, and early-Universe studies. We find that the soft X-ray emission arising from super-critical accretion is significant for the heating of the IGM, and consistent with recent constraints from the 21-cm cosmic signal. Based on the dependence on the adopted compact-object (CO) mass and accretion model, we encourage efforts in modeling ULX spectra via simulations, and their combination with detailed binary population synthesis models.

The 21-cm signal from the cosmic dawn: metallicity dependence of high-mass X-ray binaries

ABSTRACT X-rays from high-mass X-ray binaries (HMXBs) are likely the main source of heating of the intergalactic medium (IGM) during Cosmic Dawn (CD), before the completion of reionization. This Epoch of Heating (EoH; z ∼10–15) should soon be detected via the redshifted 21-cm line from neutral hydrogen, allowing us to indirectly study the properties of HMXBs in the unseen, first galaxies. Low-redshift observations, as well as theoretical models, imply that the integrated X-ray luminosity to star formation rate of HMXBs (LX/SFR) should increase in metal-poor environments, typical of early galaxies. Here, we study the impact of the metallicity (Z) dependence of LX/SFR during the EoH. For our fiducial models, galaxies with star formation rates of order 10−3–$10^{-1}\, M_\odot$ yr−1 and metallicities of order 10−3–$10^{-2}\, Z_\odot$ are the dominant contributors to the X-ray background (XRB) during this period. Different LX/SFR–Z relations result in factors of ∼ 3 differences in these ranges, as well as in the mean IGM temperature and the large-scale 21-cm power, at a given redshift. We compute mock 21-cm observations adopting as a baseline a 1000 h integration with the upcoming Square Kilometer Array (SKA) for two different LX/SFR–Z relations. We perform inference on these mock observations using the common simplification of a constant LX/SFR, finding that constant LX/SFR models can recover the IGM evolution of the more complicated LX/SFR–Z simulations only during the EoH. At z < 10, where the typical galaxies are more polluted, constant LX/SFR models overpredict the XRB and its relative contribution to the early stages of the reionization.

Robust velocity-induced acoustic oscillations at cosmic dawn

The redshifted 21-cm line of hydrogen holds great potential for the study of cosmology, as it can probe otherwise unobservable cosmic epochs. In particular, measurements of the 21-cm power spectrum during cosmic dawn---the era when the first stars were formed---will provide us with a wealth of information about the astrophysics of stellar formation, as well as the origin of fluctuations in our Universe. In addition to their usually considered density fluctuations, dark matter and baryons possess large relative velocities with respect to each other, due to the baryon acoustic oscillations (BAOs) suffered by the latter, which suppress the formation of stars during cosmic dawn, leaving an imprint on 21-cm observables during this era. Here we present 21cmvFAST, a version of the publicly available code 21cmFAST modified to account for this effect. Previous work has shown that the inclusion of relative velocities produces an acoustic modulation on the large-scale 21-cm power spectrum during cosmic dawn. By comparing analytic calculations with simulations from 21cmvFAST, here we demonstrate that this modulation takes the form of robust velocity-induced acoustic oscillations (VAOs), during both the Lyman-$\alpha$ coupling era and the subsequent epoch of heating. The unique shape of these VAOs, which is determined by the acoustic physics that generated the relative velocities, can be analytically computed and is easily distinguishable from the usual density-sourced fluctuations. We find that, under a broad range of astrophysical assumptions, the VAOs are detectable at high significance by the upcoming HERA interferometer, which could therefore confirm the presence of acoustic oscillations at cosmic dawn.

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

We present analysis of the normalised 21-cm bispectrum from fully-numerical simulations of intergalactic-medium heating by stellar sources and high-mass X-ray binaries (HMXB) during the cosmic dawn. Lyman-$\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 normalised 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 normalised 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 normalised bispectrum we see in the stellar + HMXB simulation is, we examine several other simulations - two fully-numerical simulations that include QSO sources, and two with contrasting source properties produced with the semi-numerical simulation 21cmFAST. We find the qualitative evolution of the normalised 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 normalised bispectrum during the epoch of heating, so long as the spin temperature has not saturated by $z \approx 19$.

The epoch of cosmic heating by early sources of X-rays

Observations of the 21 cm line from neutral hydrogen indicate that an Epoch of Heating (EoH) might have preceded the later Epoch of Reionization (EoR). Here we study the effects on the ionization state and the thermal history of the Intergalactic Medium (IGM) during the EoH induced by different assumptions on ionizing sources in the high redshift Universe: (i) stars, (ii) X-ray binaries (XRBs), (iii) thermal bremsstrahlung of the hot Interstellar Medium (ISM), and (iv) accreting nuclear black holes (BHs). To this aim, we post-process outputs from the ($100 h^{-1}$ cMpc)$^3$ hydrodynamical simulation MassiveBlack-II with the cosmological 3D radiative transfer code CRASH, which follows the propagation of UV and X-ray photons, computing the thermal and ionization state of hydrogen and helium through the EoH. We find that stars determine the fully ionized morphology of the IGM, while the spectrally hard XRBs pave way for efficient subsequent heating and ionization by the spectrally softer ISM. With the seeding prescription in MassiveBlack-II, BHs do not contribute significantly to either ionization or heating. With only stars, most of the IGM remains in a cold state (with a median $T=11$ K at $z=10$), however, the presence of more energetic sources raises the temperature of regions around the brightest and more clustered sources above that of the CMB, opening the possibility to observing the 21 cm signal in emission.

Simultaneously constraining the astrophysics of reionisation and the epoch of heating with 21CMMC

AbstractWe extend our MCMC sampler of 3D EoR simulations, 21CMMC, to perform parameter estimation directly on light-cones of the cosmic 21cm signal. This brings theoretical analysis one step closer to matching the expected 21-cm signal from next generation interferometers like HERA and the SKA. Using the light-cone version of 21CMMC, we quantify biases in the recovered astrophysical parameters obtained from the 21cm power spectrum when using the co-eval approximation to fit a mock 3D light-cone observation. While ignoring the light-cone effect does not bias the parameters under most assumptions, it can still underestimate their uncertainties. However, significant biases (∼few – 10 σ) are possible if all of the following conditions are met: (i) foreground removal is very efficient, allowing large physical scales (k ∼ 0.1 Mpc−1) to be used in the analysis; (ii) theoretical modelling is accurate to ∼10 per cent in the power spectrum amplitude; and (iii) the 21cm signal evolves rapidly (i.e. the epochs of reionisation and heating overlap significantly

Simultaneously constraining the astrophysics of reionization and the epoch of heating with 21CMMC

Simultaneously constraining the astrophysics of reionization and the epoch of heating with 21CMMC

High Mass X-ray Binaries: Beacons in a Stormy Universe.

The discovery of gravity waves from the mergers of black hole binaries has focused the astronomical community on the high mass X-ray binaries (HMXBs) as the potential progenitors of close pairs of compact stars. This symposium gathered experts in observational and theoretical work for a very timely review of our understanding of the processes that drive the X-ray luminosity of the diverse kinds of binaries and what evolutionary stages are revealed in the observed cases. Here I offer a condensed summary of some of the results about massive star properties, the observational categories of HMXBs, their accretion processes, their numbers in the Milky Way and other galaxies, and how they may be related to the compact binaries that merge in a burst of gravity waves.

Simulating the impact of X-ray heating during the cosmic dawn

Upcoming observations of the 21-cm signal from the Epoch of Reionization will soon provide the first direct detection of this era. This signal is influenced by many astrophysical effects, including long range X-ray heating of the intergalactic gas. During the preceding Cosmic Dawn era the impact of this heating on the 21-cm signal is particularly prominent, especially before spin temperature saturation. We present the largest-volume (349\,Mpc comoving=244~$h^{-1}$Mpc) full numerical radiative transfer simulations to date of this epoch that include the effects of helium and multi-frequency heating, both with and without X-ray sources. We show that X-ray sources contribute significantly to early heating of the neutral intergalactic medium and, hence, to the corresponding 21-cm signal. The inclusion of hard, energetic radiation yields an earlier, extended transition from absorption to emission compared to the stellar-only case. The presence of X-ray sources decreases the absolute value of the mean 21-cm differential brightness temperature. These hard sources also significantly increase the 21-cm fluctuations compared the common assumption of temperature saturation. The 21-cm differential brightness temperature power spectrum is initially boosted on large scales, before decreasing on all scales. Compared to the case of the cold, unheated intergalactic medium, the signal has lower rms fluctuations and increased non-Gaussianity, as measured by the skewness and kurtosis of the 21-cm probability distribution functions. Images of the 21-cm signal with resolution around 11~arcmin still show fluctuations well above the expected noise for deep integrations with the SKA1-Low, indicating that direct imaging of the X-ray heating epoch could be feasible.

Constraining high-redshift X-ray sources with next generation 21-cm power spectrum measurements

We use the Fisher matrix formalism and semi-numerical simulations to derive quantitative predictions of the constraints that power spectrum measurements on next-generation interferometers, such as the Hydrogen Epoch of Reionization Array (HERA) and the Square Kilometre Array (SKA), will place on the characteristics of the X-ray sources that heated the high redshift intergalactic medium. Incorporating observations between $z=5$ and $z=25$, we find that the proposed 331 element HERA and SKA phase 1 will be capable of placing $\lesssim 10\%$ constraints on the spectral properties of these first X-ray sources, even if one is unable to perform measurements within the foreground contaminated "wedge" or the FM band. When accounting for the enhancement in power spectrum amplitude from spin temperature fluctuations, we find that the observable signatures of reionization extend well beyond the peak in the power spectrum usually associated with it. We also find that lower redshift degeneracies between the signatures of heating and reionization physics lead to errors on reionization parameters that are significantly greater than previously predicted. Observations over the heating epoch are able to break these degeneracies and improve our constraints considerably. For these two reasons, 21\,cm observations during the heating epoch significantly enhance our understanding of reionization as well.

Cassini Observes the Active South Pole of Enceladus

Cassini has identified a geologically active province at the south pole of Saturn's moon Enceladus. In images acquired by the Imaging Science Subsystem (ISS), this region is circumscribed by a chain of folded ridges and troughs at approximately 55 degrees S latitude. The terrain southward of this boundary is distinguished by its albedo and color contrasts, elevated temperatures, extreme geologic youth, and narrow tectonic rifts that exhibit coarse-grained ice and coincide with the hottest temperatures measured in the region. Jets of fine icy particles that supply Saturn's E ring emanate from this province, carried aloft by water vapor probably venting from subsurface reservoirs of liquid water. The shape of Enceladus suggests a possible intense heating epoch in the past by capture into a 1:4 secondary spin/orbit resonance.

Analytic solutions for coupled linear perturbations

Analytic solutions for the evolution of cosmological linear density perturbations in the baryonic gas and collisionless dark matter are derived. The solutions are expressed in a closed form in terms of elementary functions, for arbitrary baryonic mass fraction. They are obtained assuming $\Omega=1$ and a time independent comoving Jeans wavenumber, $k_J$. By working with a time variable $\tau\equiv \ln(t^{2/3})$, the evolution of the perturbations is described by linear differential equations with constant coefficients. The new equations are then solved by means of Laplace transformation assuming that the gas and dark matter trace the same density field before a sudden heating epoch. In a dark matter dominated Universe, the ratio of baryonic to dark matter density perturbation decays with time roughly like $\exp(-5\tau/4)\propto t^{-5/6}$ to the limiting value $1/[1+(k/k_J)^2]$. For wavenumbers $k>k_J/\sqrt{24}$, the decay is accompanied with oscillations of a period $ 8\pi/\sqrt{24 (k/k_J)^2 -1}$ in $\tau$. In comparison, as $\tau $ increases in a baryonic matter dominated Universe, the ratio approaches $1-(k/k_J)^2$ for $k\le k_J$, and zero otherwise.

Primordial electrical induction heating of asteroids

Recently observed systematic trends of asteroidal composition with respect to size and heliocentric distance are interpreted according to the hypothesis that primordial asteroidal thermal evolution was driven by electromagnetic induction in a dense pre-main-sequence solar wind. It is found that the spatial distribution of induction heating corresponds well with that inferred from the present-day distribution of asteroidal compositional types when reasonable proto-Solar-System characteristics are assumed. From this match is deduced a rough constraint on the induction heating epoch's duration (≈10 5years) and the temperature-dependent electrical conductivity of the protoasteroidal material (about the same as that of the C2 meteorite Murcheson).