Fast Large Volume Simulations of the 21 cm Signal from the Reionization and pre-Reionization Epochs

While limited to low spatial resolution, the next generation low-frequency radio interferometers that target 21 cm observations during the era of reionization and prior will have instantaneous fields-of-view that are many tens of square degrees on the sky. Predictions related to various statistical measurements of the 21 cm brightness temperature must then be pursued with numerical simulations of reionization with correspondingly large volume box sizes, of order 1000 Mpc on one side. We pursue a semi-numerical scheme to simulate the 21 cm signal during and prior to Reionization by extending a hybrid approach where simulations are performed by first laying down the linear dark matter density field, accounting for the non-linear evolution of the density field based on second-order linear perturbation theory as specified by the Zel'dovich approximation, and then specifying the location and mass of collapsed dark matter halos using the excursion-set formalism. The location of ionizing sources and the time evolving distribution of ionization field is also specified using an excursion-set algorithm. We account for the brightness temperature evolution through the coupling between spin and gas temperature due to collisions, radiative coupling in the presence of Lyman-alpha photons and heating of the intergalactic medium, such as due to a background of X-ray photons. The hybrid simulation method we present is capable of producing the required large volume simulations with adequate resolution in a reasonable time so a large number of realizations can be obtained with variations in assumptions related to astrophysics and background cosmology that govern the 21 cm signal.Comment: 14 pages and 15 figures. New version to match accepted version for MNRAS. Code available in: http://www.SimFast21.or

Semi-numeric simulations of helium reionization and the fluctuating radiation background

Recent He II Lyman-alpha forest observations from 2.0 2.7. These results point to a fluctuating He-ionizing background, which may be due to the end of helium reionization of this era. We present a fast, semi-numeric procedure to approximate detailed cosmological simulations. We compute the distribution of dark matter halos, ionization state of helium, and density field at z = 3 in broad agreement with recent simulations. Given our speed and flexibility, we investigate a range of ionizing source and active quasar prescriptions. Spanning a large area of parameter space, we find order-of-magnitude fluctuations in the He II ionization rate in the post-reionization regime. During reionization, the fluctuations are even stronger and develop a bimodal distribution, in contrast to semi-analytic models and the hydrogen equivalent. These distributions indicate a low-level ionizing background even at significant He II fractions

The Impact of HI in Galaxies on 21-cm Intensity Fluctuations During the Reionisation Epoch

We investigate the impact of neutral hydrogen (HI) in galaxies on the statistics of 21-cm fluctuations using analytic and semi-numerical modelling. Following the reionisation of hydrogen the HI content of the Universe is dominated by damped absorption systems (DLAs), with a cosmic density in HI that is observed to be constant at a level equal to ~2% of the cosmic baryon density from z~1 to z~5. We show that extrapolation of this constant fraction into the reionisation epoch results in a reduction of 10-20% in the amplitude of 21-cm fluctuations over a range of spatial scales. The assumption of a different percentage during the reionisation era results in a proportional change in the 21-cm fluctuation amplitude. We find that consideration of HI in galaxies/DLAs reduces the prominence of the HII region induced shoulder in the 21-cm power spectrum (PS), and hence modifies the scale dependence of 21-cm fluctuations. We also estimate the 21cm-galaxy cross PS, and show that the cross PS changes sign on scales corresponding to the HII regions. From consideration of the sensitivity for forthcoming low-frequency arrays we find that the effects of HI in galaxies/DLAs on the statistics of 21-cm fluctuations will be significant with respect to the precision of a PS or cross PS measurement. In addition, since overdense regions are reionised first we demonstrate that the cross-correlation between galaxies and 21-cm emission changes sign at the end of the reionisation era, providing an alternative avenue to pinpoint the end of reionisation. The sum of our analysis indicates that the HI content of the galaxies that reionise the universe will need to be considered in detailed modelling of the 21-cm intensity PS in order to correctly interpret measurements from forthcoming low-frequency arrays.Comment: 11 pages, 6 figures. Submitted to MNRA

Ultraviolet Line Emission from Metals in the Low-Redshift Intergalactic Medium

We use a high-resolution cosmological simulation that includes hydrodynamics, multiphase star formation, and galactic winds to predict the distribution of metal line emission at z~0 from the intergalactic medium (IGM). We focus on two ultraviolet doublet transitions, OVI 1032,1038 and CIV 1548,1551. Emission from filaments with moderate overdensities is orders of magnitude smaller than the background, but isolated emission from enriched, dense regions with T~10^5-10^5.5 K and characteristic sizes of 50-100 kpc can be detected above the background. We show that the emission from these regions is substantially greater when we use the metallicities predicted by the simulation (which includes enrichment through galactic winds) than when we assume a uniform IGM metallicity. Luminous regions correspond to volumes that have recently been influenced by galactic winds. We also show that the line emission is clustered on scales ~1 h^-1 Mpc. We argue that although these transitions are not effective tracers of the warm-hot intergalactic medium, they do provide a route to study the chemical enrichment of the IGM and the physics of galactic winds.Comment: replaced by version to appear in ApJ (conclusions unchanged, one new figure), 16 pages (emulateapj), 11 figures, version with higher resolution figures available at http://www.tapir.caltech.edu/~sfurlane/metals/coverpage.htm

Constraining the Collisional Nature of the Dark Matter Through Observations of Gravitational Wakes

We propose to use gravitational wakes as a direct observational probe of the collisional nature of the dark matter. We calculate analytically the structure of a wake generated by the motion of a galaxy in the core of an X-ray cluster for dark matter in the highly-collisional and collisionless limits. We show that the difference between these limits can be recovered from detailed X-ray or weak lensing observations. We also discuss the sizes of sub-halos in these limits. Preliminary X-ray data on the motion of NGC 1404 through the Fornax group disfavors fluid-like dark matter but does not exclude scenarios in which the dark matter is weakly collisional.Comment: 29 pages, 3 figures, submitted to Ap

Modification of the 21-cm power spectrum by X-rays during the epoch of reionisation

We incorporate a contribution to reionization from X-rays within analytic and semi-numerical simulations of the 21-cm signal arising from neutral hydrogen during the epoch of reionization. We explore the impact that X-ray ionizations have on the power spectrum (PS) of 21-cm fluctuations by varying both the average X-ray MFP and the fractional contribution of X-rays to reionization. In general, prior to the epoch when the intergalactic medium is dominated by ionized regions (H {\sevensize II} regions), X-ray-induced ionization enhances fluctuations on spatial scales smaller than the X-ray MFP, provided that X-ray heating does not strongly supress galaxy formation. Conversely, at later times when \H2 regions dominate, small-scale fluctuations in the 21-cm signal are suppressed by X-ray ionization. Our modelling also shows that the modification of the 21-cm signal due to the presence of X-rays is sensitive to the relative scales of the X-ray MFP, and the characteristic size of \H2 regions. We therefore find that X-rays imprint an epoch and scale-dependent signature on the 21-cm PS, whose prominence depends on fractional X-ray contribution. The degree of X-ray heating of the IGM also determines the extent to which these features can be discerned. We show that the MWA will have sufficient sensitivity to detect this modification of the PS, so long as the X-ray photon MFP falls within the range of scales over which the array is most sensitive ($\sim0.1$ Mpc$^{-1}$). In cases in which this MFP takes a much smaller value, an array with larger collecting area would be required.Comment: 15 pages, 6 figures, Accepted for publication in MNRAS X-ray heating contribution now adde

Searching for the earliest galaxies in the 21 cm forest

We use a model developed by Xu et al. (2010) to compute the 21 cm line absorption signatures imprinted by star-forming dwarf galaxies (DGs) and starless minihalos (MHs). The method, based on a statistical comparison of the equivalent width (W_\nu) distribution and flux correlation function, allows us to derive a simple selection criteria for candidate DGs at very high (z >= 8) redshift. We find that ~ 18% of the total number of DGs along a line of sight to a target radio source (GRB or quasar) can be identified by the condition W_\nu < 0; these objects correspond to the high-mass tail of the DG distribution at high redshift, and are embedded in large HII regions. The criterion W_\nu > 0.37 kHz instead selects ~ 11% of MHs. Selected candidate DGs could later be re-observed in the near-IR by the JWST with high efficiency, thus providing a direct probe of the most likely reionization sources.Comment: 8 pages, 3 figures. Accepted for publication in Science in China Series

Formation of the seed black holes: a role of quark nuggets?

Strange quark nuggets (SQNs) could be the relics of the cosmological QCD phase transition, and they could very likely be the candidate of cold quark matter if survived the cooling of the later Universe, although the formation and evolution of these SQNs depend on the physical state of the hot QGP (quark-gluon plasma) phase and the state of cold quark matter. We reconsider the possibility of SQNs as cold dark matter, and find that the formation of black holes in primordial halos could be significantly different from the standard scenario. In a primordial halo, the collision between gas and SQNs could be frequent enough, and thus the viscosity acting on each SQN would decrease its angular momentum and make it to sink into the center of the halo, as well as heat the gas. The SQNs with baryon numbers less than $10^{35}$ could assemble in the center of the halo before the formation of primordial stars. A black hole could form by merger of these SQNs, and then its mass could quickly become about $10^3\ M_\odot$ or higher, by accreting the surrounding SQNs or gas. The black holes formed in this way could be the seeds for the supermassive black holes at redshift as high as $z\sim 6$.Comment: 15 page

21-cm cosmology

Imaging the Universe during the first hundreds of millions of years remains one of the exciting challenges facing modern cosmology. Observations of the redshifted 21 cm line of atomic hydrogen offer the potential of opening a new window into this epoch. This would transform our understanding of the formation of the first stars and galaxies and of the thermal history of the Universe. A new generation of radio telescopes is being constructed for this purpose with the first results starting to trickle in. In this review, we detail the physics that governs the 21 cm signal and describe what might be learnt from upcoming observations. We also generalize our discussion to intensity mapping of other atomic and molecular lines.Comment: 64 pages, 20 figures, submitted to Reports on Progress in Physics, comments welcom

Evolution and Distribution of Magnetic Fields from AGNs in Galaxy Clusters. I. The Effect of Injection Energy and Redshift

We present a series of cosmological magnetohydrodynamic (MHD) simulations that simultaneously follow the formation of a galaxy cluster and evolution of magnetic fields ejected by an Active Galactic Nucleus (AGN). Specifically, we investigate the influence of both the epoch of AGN (z $\sim$ 3-0.5) and the AGN energy ($\sim$ 3 $\times$ 10$^{57}$ - 2 $\times$ 10$^{60}$ ergs)on the final magnetic field distribution in a relatively massive cluster (M$_{vir}$ $\sim$10$^{15}$ M$_\odot$). We find that as long as the AGN magnetic fields are ejected before the major mergers in the cluster formation history, magnetic fields can be transported throughout the cluster and can be further amplified by the intra-cluster medium (ICM) turbulence cause by hierarchical mergers during the cluster formation process. The total magnetic energy in the cluster can reach $\sim$ $10^{61}$ ergs, with micro Gauss fields distributed over $\sim$ Mpc scale. The amplification of the total magnetic energy by the ICM turbulence can be significant, up to $\sim$1000 times in some cases. Therefore even weak magnetic fields from AGNs can be used to magnetize the cluster to the observed level. The final magnetic energy in the ICM is determined by the ICM turbulent energy, with a weak dependence on the AGN injection energy. We discuss the properties of magnetic fields throughout the cluster and the synthetic Faraday rotation measure maps they produce. We also show that high spatial resolution over most of the magnetic regions of the cluster is very important to capture the small scale dynamo process and maintain the magnetic field structure in our simulations.Comment: 38 pages, 18 figures, Accepted for publication in Ap