Low-density Intergalactic Medium Research Articles

Emergence of the temperature–density relation in the low-density intergalactic medium

ABSTRACT We examine the evolution of the phase diagram of the low-density intergalactic medium during the Epoch of Reionization in simulation boxes with varying reionization histories from the Cosmic Reionization on Computers project. The probability density function (PDF) of gas temperature at fixed density exhibits two clear modes: a warm and a cold temperature mode, corresponding to the gas inside and outside of ionized bubbles. We find that the transition between the two modes is ‘universal’ in the sense that its timing is accurately parametrized by the value of the volume-weighted neutral fraction for any reionization history. This ‘universality’ is more complex than just a reflection of the fact that ionized gas is warm and neutral gas is cold: it holds for the transition at a fixed value of gas density, and gas at different densities transitions from the cold to the warm mode at different values of the neutral fraction, reflecting a non-trivial relationship between the ionization history and the evolving gas density PDF. Furthermore, the ‘emergence’ of the tight temperature–density relation in the warm mode is also approximately ‘universally’ controlled by the volume-weighted neutral fraction for any reionization history. In particular, the ‘emergence’ of the temperature–density relation (as quantified by the rapid decrease in its width) occurs when the neutral fraction is 10−4 ≲ XH i ≲ 10−3 for any reionization history. Our results indicate that the neutral fraction is a primary quantity controlling the various properties of the temperature–density relation, regardless of reionization history.

Constraining the reionization and thermal history of the Universe using a seminumerical photon-conserving code script

ABSTRACT Given that the reionization history of cosmic hydrogen is yet to be stringently constrained, it is worth checking the prospects of doing so using physically motivated models and available observational data. For this purpose, we use an extended version of the explicitly photon-conserving seminumerical model of reionization, script, which also includes thermal evolution of the intergalactic medium (IGM). The model incorporates the effects of inhomogeneous recombination and radiative feedback self-consistently and is characterized by five free parameters (two for the redshift-dependent ionization efficiency, two for the ionizing escape fraction, and another for reionization temperature increment). We constrain these free parameters by simultaneously matching with various observational probes, e.g. estimates of the ionized hydrogen fraction, the CMB scattering optical depth and the galaxy UV luminosity function. In addition, we include the low-density IGM temperature measurements obtained from Lyman-α absorption spectra at z∼5.5, a probe not commonly used for Bayesian analysis of reionization parameters. We find that the interplay of the various data sets, particularly inclusion of the temperature data, leads to tightening of the parameter constraints. Our default models prefer a late end of reionization (at z ≲ 6), in agreement with other recent studies. We can also derive constraints on the duration of reionization, $\Delta z=1.81^{+0.51}_{-0.67}$ and the mid-point of reionization, $z_{\mathrm{mid}}=7.0^{+0.30}_{-0.40}$. The constraints can be further tightened by including other available and upcoming data sets.

The Hosts of X-Ray Absorption Lines Toward AGNs

Most baryonic matter in the universe exists in gaseous form and can be found in structures such as galactic halos and the low-density intergalactic medium. Proposed X-ray spectroscopy missions such as Athena, Arcus, and Lynx will have the capability to identify absorption lines in spectra toward bright active galactic nuclei (AGNs), which can be used as tools to probe this missing matter. In this study, we examine the optical fields surrounding 15 primary observing targets and identify the foreground galaxies and galaxy groups that are potential hosts of absorption. We record the basic properties of the potential host and their angular and physical separation from the AGN line of sight. This process is done by marking the location of various galaxies and groups in optical images of the field surrounding the target and plotting their angular separation versus redshift to gauge the physical proximity to the background source. We identify the surrounding objects according to those with measured redshifts and those that require accurate redshift measurements.

Constraining blazar heating with the 2 ≲z≲ 3 Lyman-α forest

ABSTRACTThe intergalactic medium (IGM) acts like a calorimeter recording energy injection by cosmic structure formation, shocks and photoheating from stars and active galactic nuclei. It was recently proposed that spatially inhomogeneous TeV-blazars could significantly heat up the underdense IGM, resulting in patches of both cold and warm IGM around z ≃ 2−3. The goal of this study is to compare predictions of different blazar heating models with recent observations of the IGM. We perform a set of cosmological simulations and carefully compute mock observables of the Lyman-α (Lyα) forest. We perform a detailed assessment of different systematic uncertainties which typically impact this type of observables and find that they are smaller than the differences between our models. We find that our inhomogeneous blazar heating model is in good agreement with the Lyα line properties and the rescaled flux probability distribution function at high redshift (2.5 < z < 3) but that our blazar heating models are challenged by lower redshift data (2 < z < 2.5). Our results could be explained by HeII reionization although state-of-the-art models fall short on providing enough heating to the low-density IGM, thus motivating further radiative transfer studies of inhomogeneous HeII reionization. If blazars are indeed hosted by group-mass haloes of $2\times 10^{13}\, \rm {M}_\odot$, a later onset of blazar heating in comparison with previous models would be favoured, which could bring our findings here in agreement with the evidence of blazar heating from local gamma-ray observations.

A Hydro-particle-mesh Code for Efficient and Rapid Simulations of the Intracluster Medium

Abstract We introduce the cosmological HYPER code based on an innovative hydro-particle-mesh (HPM) algorithm for efficient and rapid simulations of gas and dark matter. For the HPM algorithm, we update the approach of Gnedin & Hui to expand the scope of its application from the lower-density intergalactic medium (IGM) to the higher-density intracluster medium (ICM). While the original algorithm tracks only one effective particle species, the updated version separately tracks the gas and dark matter particles, as they do not exactly trace each other on small scales. For the approximate hydrodynamics solver, the pressure term in the gas equations of motion is calculated using robust physical models. In particular, we use a dark matter halo model, ICM pressure profile, and IGM temperature–density relation, all of which can be systematically varied for parameter-space studies. We show that the HYPER simulation results are in good agreement with the halo model expectations for the density, temperature, and pressure radial profiles. Simulated galaxy cluster scaling relations for Sunyaev–Zel’dovich (SZ) and X-ray observables are also in good agreement with mean predictions, with scatter comparable to that found in hydrodynamic simulations. HYPER also produces lightcone catalogs of dark matter halos and full-sky tomographic maps of the lensing convergence, SZ effect, and X-ray emission. These simulation products are useful for testing data analysis pipelines, generating training data for machine learning, understanding selection and systematic effects, and for interpreting astrophysical and cosmological constraints.

Probing the thermal history during reionization using a seminumerical photon-conserving code script

ABSTRACT The ionization and thermal state of the intergalactic medium (IGM) during the epoch of reionization has been of interest in recent times because of their close connection to the first stars. We present in this paper a seminumerical code which computes the large-scale temperature and ionized hydrogen fields in a cosmologically representative volume accounting for the patchiness in these quantities arising from reionization. The code is an extension to a previously developed version for studying the growth of ionized regions, namely, Semi Numerical Code for ReionIzation with PhoTon Conservation (script). The main additions in the present version are the inhomogeneous recombinations which are essential for temperature calculations. This extended version of script also implements physical consequences of photoheating during reionization, e.g. radiative feedback. These enhancements allow us to predict observables which were not viable with the earlier version. These include the faint-end of the ultra-violet luminosity function of galaxies (which can get affected by the radiative feedback) and the temperature-density relation of the low-density IGM at z ∼ 6. We study the effect of varying the free parameters and prescriptions of our model on a variety of observables. The conclusion of our analysis is that it should be possible to put constraints on the evolution of thermal and ionization state of the IGM using available observations accounting for all possible variations in the free parameters. A detailed exploration of the parameter space will be taken up in the future.

Recovering Density Fields inside Quasar Proximity Zones at z ∼ 6

Abstract The matter density field at z ∼ 6 is very challenging to probe. One of the traditional probes of the low density intergalactic medium that works successfully at lower redshifts is the Lyα forest in quasar spectra. However, at the end of reionization, the residual neutral hydrogen usually creates saturated absorption, thus much of the information about the gas density is lost. Luckily, in a quasar proximity zone, the ionizing radiation is exceptionally intense, thus creating a large region with nonzero transmitted flux. In this study we use the synthetic spectra from simulations to investigate how to recover the density fluctuations inside the quasar proximity zones. We show that, under ideal conditions, the density can be recovered accurately with a small scatter. We also discuss how systematics such as the quasar continuum fitting and reionization models affect the results. This study shows that by analyzing the absorption features inside quasar proximity zones we can potentially constrain quasar properties and the environments they reside in.

Effect of separate initial conditions on the lyman-α forest in simulations

ABSTRACT Using a set of high resolution simulations, we quantify the effect of species-specific initial transfer functions on probes of the intergalactic medium (IGM) via the Lyman-α forest. We focus on redshifts 2–6, after H i reionization. We explore the effect of these initial conditions on measures of the thermal state of the low density IGM: the curvature, Doppler width cutoff, and Doppler width distribution. We also examine the matter and flux power spectrum, and potential consequences for constraints on warm dark matter models. We find that the curvature statistic is at most affected at the $\approx 2{{\ \rm per\ cent}}$ level at z = 6. The Doppler width cutoff parameters are affected by $\approx 5{{\ \rm per\ cent}}$ for the intercept, and $\approx 8{{\ \rm per\ cent}}$ for the fit slope, though this is subdominant to sample variation. The Doppler width distribution shows a $\approx 30{{\ \rm per\ cent}}$ effect at z = 3, however the distribution is not fully converged with simulation box size and resolution. The flux power spectrum is at most affected by $\approx 5{{\ \rm per\ cent}}$ at high redshift and small scales. We discuss numerical convergence with simulation parameters.

New constraints on the free-streaming of warm dark matter from intermediate and small scale Lyman- α forest data

We present new measurements of the free-streaming of warm dark matter (WDM) from Lyman-$\alpha$ flux-power spectra. We use data from the medium resolution, intermediate redshift XQ-100 sample observed with the X-shooter spectrograph ($z=3 - 4.2$) and the high-resolution, high-redshift sample used in Viel et al. (2013) obtained with the HIRES/MIKE spectrographs ($z=4.2 - 5.4$). Based on further improved modelling of the dependence of the Lyman-$\alpha$ flux-power spectrum on the free-streaming of dark matter, cosmological parameters, as well as the thermal history of the intergalactic medium (IGM) with hydrodynamical simulations, we obtain the following limits, expressed as the equivalent mass of thermal relic WDM particles. The XQ-100 flux power spectrum alone gives a lower limit of 1.4 keV, the re-analysis of the HIRES/MIKE sample gives 4.1 keV while the combined analysis gives our best and significantly strengthened lower limit of 5.3 keV (all 2$\sigma$ C.L.). The further improvement in the joint analysis is partly due to the fact that the two data sets have different degeneracies between astrophysical and cosmological parameters that are broken when the data sets are combined, and more importantly on chosen priors on the thermal evolution. These results all assume that the temperature evolution of the IGM can be modelled as a power law in redshift. Allowing for a non-smooth evolution of the temperature of the IGM with sudden temperature changes of up to 5000K reduces the lower limit for the combined analysis to 3.5 keV. A WDM with smaller thermal relic masses would require, however, a sudden temperature jump of $5000\,K$ or more in the narrow redshift interval $z=4.6-4.8$, in disagreement with observations of the thermal history based on high-resolution resolution Lyman-$\alpha$ forest data and expectations for photo-heating and cooling in the low density IGM at these redshifts.

Simulated metal and H i absorption lines at the conclusion of reionization

We present a theoretical study of intergalactic metal absorption lines imprinted in the spectra of distant quasars during and after the Epoch of Reionization (EoR). We use high resolution hydrodynamical simulations at high redshift ($4 <z<8$), assuming a uniform UV background Haardt--Madau 12, post-processing with CLOUDY photoionization models and Voigt profile fitting to accurately calculate column densities of the ions CII, CIV, SiII, SiIV and OI in the intergalactic medium (IGM). In addition, we generate mock observations of neutral Hydrogen (HI) at $z<6$. Our simulations successfully reproduce the evolution of the cosmological mass density ($\Omega$) of CII and CIV, with $\Omega_{CII}$ exceeding $\Omega_{CIV}$ at $z >6$, consistent with the current picture of the tail of the EoR. The simulated CII exhibits a bimodal distribution with large absorptions in and around galaxies, and some traces in the lower density IGM. We find some discrepancies between the observed and simulated column density relationships among different ionic species at $z=6$, probably due to uncertainties in the assumed UV background. Finally, our simulations are in good agreement with observations of the HI column density distribution function at $z = 4$ and the HI cosmological mass density $\Omega_{HI}$ at $4 < z < 6$.

Exploring the thermal state of the low-density intergalactic medium atz = 3 with an ultrahigh signal-to-noise QSO spectrum

At low densities the standard ionisation history of the intergalactic medium (IGM) predicts a decreasing temperature of the IGM with decreasing density once hydrogen (and helium) reionisation is complete. Heating the high-redshift, low-density IGM above the temperature expected from photo-heating is difficult, and previous claims of high/rising temperatures in low density regions of the Universe based on the probability density function (PDF) of the opacity in Lyman-$\alpha$ forest data at $2<z<4$ have been met with considerable scepticism, particularly since they appear to be in tension with other constraints on the temperature-density relation (TDR). We utilize here an ultra-high signal-to-noise spectrum of the QSO HE0940-1050 and a novel technique to study the low opacity part of the PDF. We show that there is indeed evidence (at 90% confidence level) that a significant volume fraction of the under-dense regions at $z \sim 3$ has temperatures as high or higher than those at densities comparable to the mean and above. We further demonstrate that this conclusion is nevertheless consistent with measurements of a slope of the TDR in over-dense regions that imply a decreasing temperature with decreasing density, as expected if photo-heating of ionised hydrogen is the dominant heating process. We briefly discuss implications of our findings for the need to invoke either spatial temperature fluctuations, as expected during helium reionization, or additional processes that heat a significant volume fraction of the low-density IGM.

Models of the thermal evolution of the intergalactic medium after reionization

Recent years have brought more precise temperature measurements of the low-density intergalactic medium (IGM). These new measurements constrain the processes that heated the IGM, such as the reionization of H I and of He II. We present a semi-analytical model for the thermal history of the IGM that follows the photoheating history of primordial gas. Our model adopts a multizone approach that, compared to previous models, more accurately captures the inhomogeneous heating and cooling of the IGM during patchy reionization. We compare our model with recent temperature measurements spanning z = 1.6-4.8, finding that these measurements are consistent with scenarios in which the He II was reionized at z = 3-4 by quasars. Significantly longer duration or higher redshift He II reionization scenarios are ruled out by the measurements. For hydrogen reionization, we find that only low-redshift and high-temperature scenarios are excluded. For example, a model in which the IGM was heated to 30 000 K when an ionization front passed, and with hydrogen reionization occurring over 6 < z < 9, is ruled out. Finally, we place constraints on how much heating could owe to TeV blazars, cosmic rays, and other non-standard mechanisms. We find that by z = 2, a maximum of 1 eV of additional heat could be injected per baryon over standard photoheating-only models, with this limit becoming <~0.5 eV at z > 3.

On the intergalactic temperature–density relation

Cosmological simulations of the low-density intergalactic medium exhibit a strikingly tight power-law relation between temperature and density that holds over two decades in density. It is found that this relation should roughly apply Delta z ~ 1-2 after a reionization event, and this limiting behavior has motivated the power-law parameterizations used in most analyses of the Ly-alpha forest. This relation has been explained by using equations linearized in the baryonic overdensity (which does not address why a tight power-law relation holds over two decades in density) or by equating the photoheating rate with the cooling rate from cosmological expansion (which we show is incorrect). Previous explanations also did not address why recombination cooling and Compton cooling off of the cosmic microwave background, which are never negligible, do not alter the character of this relation. We provide an understanding for why a tight power-law relation arises for unshocked gas at all densities for which collisional cooling is unimportant. We also use our results to comment on (1) how quickly fluctuations in temperature redshift away after reionization processes, (2) how much shock heating occurs in the low-density intergalactic medium, and (3) how the temperatures of collapsing gas parcels evolve.

Mapping metals at high redshift with far-infrared lines

Cosmic metal enrichment is one of the key physical processes regulating galaxy formation and the evolution of the intergalactic medium (IGM). However, determining the metal content of the most distant galaxies has proven so far almost impossible; also, absorption line experiments at $z\sim6$ become increasingly difficult because of instrumental limitations and the paucity of background quasars. With the advent of ALMA, far-infrared emission lines provide a novel tool to study early metal enrichment. Among these, the [CII] line at 157.74 $\mu$m is the most luminous line emitted by the interstellar medium of galaxies. It can also resonant scatter CMB photons inducing characteristic intensity fluctuations ($\Delta I/I_{CMB}$) near the peak of the CMB spectrum, thus allowing to probe the low-density IGM. We compute both [CII] galaxy emission and metal-induced CMB fluctuations at $z\sim 6$ by using Adaptive Mesh Refinement cosmological hydrodynamical simulations and produce mock observations to be directly compared with ALMA BAND6 data ($\nu_{obs}\sim 272$ GHz). The [CII] line flux is correlated with $M_{UV}$ as $\log(F_{peak}/\mu{\rm Jy})=-27.205-2.253\,M_{UV}-0.038\,M_{UV}^2$. Such relation is in very good agreement with recent ALMA observations (e.g. Maiolino et al. 2015; Capak et al. 2015) of $M_{UV}<-20$ galaxies. We predict that a $M_{UV}=-19$ ($M_{UV}=-18$) galaxy can be detected at $4\sigma$ in $\simeq40$ (2000) hours, respectively. CMB resonant scattering can produce $\simeq\pm 0.1\,\mu$Jy/beam emission/absorptions features that are very challenging to be detected with current facilities. The best strategy to detect these signals consists in the stacking of deep ALMA observations pointing fields with known $M_{UV}\simeq-19$ galaxies. This would allow to simultaneously detect both [CII] emission from galactic reionization sources and CMB fluctuations produced by $z\sim6$ metals.

ON MODELING AND MEASURING THE TEMPERATURE OF THEz∼ 5 INTERGALACTIC MEDIUM

The temperature of the low-density intergalactic medium (IGM) at high redshift is sensitive to the timing and nature of hydrogen and HeII reionization, and can be measured from Lyman-alpha forest absorption spectra. Since the memory of intergalactic gas to heating during reionization gradually fades, measurements as close as possible to reionization are desirable. In addition, measuring the IGM temperature at sufficiently high redshifts should help to isolate the effects of hydrogen reionization since HeII reionization starts later, at lower redshift. Motivated by this, we model the IGM temperature at z>5 using semi-numeric models of patchy reionization. We construct mock Lyman-alpha forest spectra from these models and consider their observable implications. We find that the small-scale structure in the Lyman-alpha forest is sensitive to the temperature of the IGM even at redshifts where the average absorption in the forest is as high as 90%. We forecast the accuracy at which the z~5 IGM temperature can be measured using existing samples of high resolution quasar spectra, and find that interesting constraints are possible. For example, an early reionization model in which reionization ends at z~10 should be distinguishable -- at high statistical significance -- from a lower redshift model where reionization completes at z~6. We discuss improvements to our modeling that may be required to robustly interpret future measurements.

THE BIGGEST EXPLOSIONS IN THE UNIVERSE

Supermassive primordial stars are expected to form in a small fraction of massive protogalaxies in the early universe, and are generally conceived of as the progenitors of the seeds of supermassive black holes (BHs). Supermassive stars with masses of ~55,000 M_Sun, however, have been found to explode and completely disrupt in a supernova (SN) with an energy of up to ~10^55 erg instead of collapsing to a BH. Such events, ~10,000 times more energetic than typical SNe today, would be among the biggest explosions in the history of the universe. Here we present a simulation of such a SN in two stages. Using the RAGE radiation hydrodynamics code we first evolve the explosion from an early stage through the breakout of the shock from the surface of the star until the blast wave has propagated out to several parsecs from the explosion site, which lies deep within an atomic cooling dark matter (DM) halo at z ~ 15. Then, using the GADGET cosmological hydrodynamics code we evolve the explosion out to several kiloparsecs from the explosion site, far into the low-density intergalactic medium. The host DM halo, with a total mass of 4 x 10^7 M_Sun, much more massive than typical primordial star-forming halos, is completely evacuated of high density gas after < 10 Myr, although dense metal-enriched gas recollapses into the halo, where it will likely form second-generation stars with metallicities of ~ 0.05 Z_Sun after > 70 Myr. The chemical signature of supermassive star explosions may be found in such long-lived second-generation stars today.

Generation of Primordial Magnetic Fields on Linear Overdensity Scales

Magnetic fields appear to be present in all galaxies and galaxy clusters. Recent measurements indicate that a weak magnetic field may be present even in the smooth low density intergalactic medium. One explanation for these observations is that a seed magnetic field was generated by some unknown mechanism early in the life of the Universe, and was later amplified by various dynamos in nonlinear objects like galaxies and clusters. We show that a primordial magnetic field is expected to be generated in the early Universe on purely linear scales through vorticity induced by scale-dependent temperature fluctuations, or equivalently, a spatially varying speed of sound of the gas. Residual free electrons left over after recombination tap into this vorticity to generate magnetic field via the Biermann battery process. Although the battery operates even in the absence of any relative velocity between dark matter and gas at the time of recombination, the presence of such a relative velocity modifies the predicted spatial power spectrum of the magnetic field. At redshifts of order a few tens, we estimate a root mean square field strength of order 10(-25)-10(-24) G on comoving scales ~10 kpc. This field, which is generated purely from linear perturbations, is expected to be amplified significantly after reionization, and to be further boosted by dynamo processes during nonlinear structure formation.

Planckintermediate results

About half of the baryons of the Universe are expected to be in the form of filaments of hot and low density intergalactic medium. Most of these baryons remain undetected even by the most advanced X-ray observatories which are limited in sensitivity to the diffuse low density medium. The Planck satellite has provided hundreds of detections of the hot gas in clusters of galaxies via the thermal Sunyaev-Zel'dovich (tSZ) effect and is an ideal instrument for studying extended low density media through the tSZ effect. In this paper we use the Planck data to search for signatures of a fraction of these missing baryons between pairs of galaxy clusters. Cluster pairs are good candidates for searching for the hotter and denser phase of the intergalactic medium (which is more easily observed through the SZ effect). Using an X-ray catalogue of clusters and the Planck data, we select physical pairs of clusters as candidates. Using the Planck data we construct a local map of the tSZ effect centered on each pair of galaxy clusters. ROSAT data is used to construct X-ray maps of these pairs. After having modelled and subtracted the tSZ effect and X-ray emission for each cluster in the pair we study the residuals on both the SZ and X-ray maps. For the merging cluster pair A399-A401 we observe a significant tSZ effect signal in the intercluster region beyond the virial radii of the clusters. A joint X-ray SZ analysis allows us to constrain the temperature and density of this intercluster medium. We obtain a temperature of kT = 7.1 +- 0.9, keV (consistent with previous estimates) and a baryon density of (3.7 +- 0.2)x10^-4, cm^-3. The Planck satellite mission has provided the first SZ detection of the hot and diffuse intercluster gas.

Synthetic X-ray and radio maps for two different models of Stephan's Quintet

We present simulations of the compact galaxy group Stephan's Quintet (SQ) including magnetic fields, performed with the N-body/smoothed particle hydrodynamics (SPH) code \textsc{Gadget}. The simulations include radiative cooling, star formation and supernova feedback. Magnetohydrodynamics (MHD) is implemented using the standard smoothed particle magnetohydrodynamics (SPMHD) method. We adapt two different initial models for SQ based on Renaud et al. and Hwang et al., both including four galaxies (NGC 7319, NGC 7320c, NGC 7318a and NGC 7318b). Additionally, the galaxies are embedded in a magnetized, low density intergalactic medium (IGM). The ambient IGM has an initial magnetic field of $10^{-9}$ G and the four progenitor discs have initial magnetic fields of $10^{-9} - 10^{-7}$ G. We investigate the morphology, regions of star formation, temperature, X-ray emission, magnetic field structure and radio emission within the two different SQ models. In general, the enhancement and propagation of the studied gaseous properties (temperature, X-ray emission, magnetic field strength and synchrotron intensity) is more efficient for the SQ model based on Renaud et al., whose galaxies are more massive, whereas the less massive SQ model based on Hwang et al. shows generally similar effects but with smaller efficiency. We show that the large shock found in observations of SQ is most likely the result of a collision of the galaxy NGC 7318b with the IGM. This large group-wide shock is clearly visible in the X-ray emission and synchrotron intensity within the simulations of both SQ models. The order of magnitude of the observed synchrotron emission within the shock front is slightly better reproduced by the SQ model based on Renaud et al., whereas the distribution and structure of the synchrotron emission is better reproduced by the SQ model based on Hwang et al..

THE COSMOLOGICAL IMPACT OF LUMINOUS TeV BLAZARS. II. REWRITING THE THERMAL HISTORY OF THE INTERGALACTIC MEDIUM

The Universe is opaque to extragalactic very high-energy gamma rays (VHEGRs, E>100 GeV) because they annihilate and pair produce on the extragalactic background light. The resulting ultra-relativistic pairs are assumed to lose energy through inverse Compton scattering of CMB photons. In Broderick et al. (2011, Paper I of this three paper series), we argued that instead powerful plasma instabilities in the ultra-relativistic pair beam dissipate the kinetic energy of the TeV-generated pairs locally, heating the intergalactic medium (IGM). Here, we explore the effect of this heating upon the thermal history of the IGM. We collate the observed extragalactic VHEGR sources to determine a local VHEGR heating rate and correct for the pointed nature of VHEGR observations using Fermi observations of high and intermediate peaked BL Lacs. Because the local extragalactic VHEGR flux is dominated by TeV blazars, we tie the TeV blazar luminosity density to the quasar luminosity density, and produce a VHEGR heating rate as a function of redshift. This heating is relatively homogeneous for z<~4 with increasing spatial variation at higher redshift (order unity at z~6). This new heating process dominates photoheating at low redshift and the inclusion of TeV blazar heating qualitatively and quantitatively changes the structure and history of the IGM. TeV blazars produce a uniform volumetric heating rate that is sufficient to increase the temperature of the mean density IGM by nearly an order of magnitude, and at low densities by substantially more, naturally producing an inverted equation of state inferred by observations of the Ly-alpha forest, a feature that is difficult to reconcile with standard reionization models. Finally, we close with a discussion on the possibility of detecting this hot low-density IGM, but find that such measurements are currently not feasible. (abridged)