Source author record

Rennan Barkana

Rennan Barkana appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

40works
3topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

40 published item(s)

preprint2023arXiv

Strong 21-cm fluctuations and anisotropy due to the line-of-sight effect of radio galaxies at cosmic dawn

The reported detection of the global 21-cm signal by the EDGES collaboration is significantly stronger than standard astrophysical predictions. One possible explanation is an early radio excess above the cosmic microwave background. Such a radio background could have been produced by high redshift galaxies, if they were especially efficient in producing low-frequency synchrotron radiation. We have previously studied the effects of such an inhomogeneous radio background on the 21-cm signal; however, we made a simplifying assumption of isotropy of the background seen by each hydrogen cloud. Here we perform a complete calculation that accounts for the fact that the 21-cm absorption occurs along the line of sight, and is therefore sensitive to radio sources lying behind each absorbing cloud. We find that the complete calculation strongly enhances the 21-cm power spectrum during cosmic dawn, by up to two orders of magnitude; on the other hand, the effect on the global 21-cm signal is only at the $5\%$ level. In addition to making the high-redshift 21-cm fluctuations potentially more easily observable, the line of sight radio effect induces a new anisotropy in the 21-cm power spectrum. While these effects are particularly large for the case of an extremely-enhanced radio efficiency, they make it more feasible to detect even a moderately-enhanced radio efficiency in early galaxies. This is especially relevant since the EDGES signal has been contested by the SARAS experiment.

preprint2021arXiv

Effect of the cosmological transition to metal-enriched star-formation on the hydrogen 21-cm signal

Mapping Cosmic Dawn with 21-cm tomography offers an exciting new window into the era of primordial star formation. However, self-consistent implementation of both the process of star formation and the related 21-cm signal is challenging, due to the multi-scale nature of the problem. In this study, we develop a flexible semi-analytical model to follow the formation of the first stars and the process of gradual transition from primordial to metal-enriched star formation. For this transition we use different in scenarios with varying time-delays (or recovery times) between the first supernovae and the formation of the second generation of stars. We use recovery times between 10 and 100\,Myr and find that these delays have a strong impact on the redshift at which the transition to metal-enriched star formation occurs. We then explore the effect of this transition on the 21-cm signal and find that the recovery time has a distinctive imprint in the signal. Together with an improved understanding of how this time-delay relates to the properties of Population~III stars, future 21-cm observations can give independent constraints on the earliest epoch of star formation.

preprint2021arXiv

Shot noise and scatter in the star formation efficiency as a source of 21-cm fluctuations

The 21-cm signal from cosmic dawn and the epoch of reionization (EoR) probes the characteristics of the high redshift galaxy population. Many of the astrophysical properties of galaxies at high redshifts are currently unconstrained due to the lack of observations. This creates a vast space of possible astrophysical scenarios where the 21-cm signal needs to be modeled in order to plan for, and eventually fit, future observations. This is done with fast numerical methods which make simplifying approximations for the underlying physical processes. In this work we quantify the effect of Poisson fluctuations and scatter in the star formation efficiency; while Poisson fluctuations are included in some works and not in others, scatter in the star formation efficiency is usually neglected, and all galaxies of a given mass are assumed to have the same properties. We show that both features can have a significant effect on the 21-cm power spectrum, most importantly in scenarios where the signal is dominated by massive galaxies. Scatter in the star formation efficiency does not simply enhance the effect of Poisson fluctuations; for example we show that the power spectrum shape at cosmic dawn has a feature corresponding to the width of the galaxy brightness distribution. We also discuss some of the consequences for 21-cm imaging, and the signature of reduced correlation between the density and radiation fields.

preprint2020arXiv

Emulating the Global 21-cm Signal from Cosmic Dawn and Reionization

The 21-cm signal of neutral hydrogen is a sensitive probe of the Epoch of Reionization (EoR) and Cosmic Dawn. Currently operating radio telescopes have ushered in a data-driven era of 21-cm cosmology, providing the first constraints on the astrophysical properties of sources that drive this signal. However, extracting astrophysical information from the data is highly non-trivial and requires the rapid generation of theoretical templates over a wide range of astrophysical parameters. To this end emulators are often employed, with previous efforts focused on predicting the power spectrum. In this work we introduce 21cmGEM - the first emulator of the global 21-cm signal from Cosmic Dawn and the EoR. The smoothness of the output signal is guaranteed by design. We train neural networks to predict the cosmological signal using a database of ~30,000 simulated signals which were created by varying seven astrophysical parameters: the star formation efficiency and the minimal mass of star-forming halos; the efficiency of the first X-ray sources and their spectrum parameterized by spectral index and the low energy cutoff; the mean free path of ionizing photons and the CMB optical depth. We test the performance with a set of ~2,000 simulated signals, showing that the relative error in the prediction has an r.m.s. of 0.0159. The algorithm is efficient, with a running time per parameter set of 0.16 sec. Finally, we use the database of models to check the robustness of relations between the features of the global signal and the astrophysical parameters that we previously reported.

preprint2020arXiv

Quantiles as Robust Probes of Non-Gaussianity in 21-cm Images

The early epoch in which the first stars and galaxies formed is among the most exciting unexplored eras of the Universe. A major research effort focuses on probing this era with the 21-cm spectral line of hydrogen. While most research focused on statistics like the 21-cm power spectrum or the sky-averaged global signal, there are other ways to analyze tomographic 21-cm maps, which may lead to novel insights. We suggest statistics based on quantiles as a method to probe non-Gaussianities of the 21-cm signal. We show that they can be used in particular to probe the variance, skewness, and kurtosis of the temperature distribution, but are more flexible and robust than these standard statistics. We test these statistics on a range of possible astrophysical models, including different galactic halo masses, star-formation efficiencies, and spectra of the X-ray heating sources, plus an exotic model with an excess early radio background. Simulating data with angular resolution and thermal noise as expected for the Square Kilometre Array (SKA), we conclude that these statistics can be measured out to redshifts above 20 and offer a promising statistical method for probing early cosmic history.

preprint2019arXiv

Extracting the Global Signal from 21-cm Fluctuations: the Multi-Tracer Approach

The multi-tracer technique employs a ratio of densities of two differently biased galaxy samples that trace the same underlying matter density field, and was proposed to alleviate the cosmic variance problem. Here we propose a novel application of this approach, applying it to two different tracers one of which is the 21-cm signal of neutral hydrogen from the epochs of reionization and comic dawn. The second tracer is assumed to be a sample of high-redshift galaxies, but the approach can be generalized and applied to other high-redshift tracers. We show that the anisotropy of the ratio of the two density fields can be used to measure the sky-averaged 21-cm signal, probe the spectral energy distribution of radiative sources that drive this signal, and extract large-scale properties of the second tracer, e.g., the galaxy bias. Using simulated 21-cm maps and mock galaxy samples, we find that the method works well for an idealized galaxy survey. However, in the case of a more realistic galaxy survey which only probes highly biased luminous galaxies, the inevitable Poisson noise makes the reconstruction far more challenging. This difficulty can be mitigated with the greater sensitivity of future telescopes along with larger survey volumes.

preprint2016arXiv

Constraining the redshifted 21-cm signal with the unresolved soft X-ray background

We use the observed unresolved cosmic X-ray background (CXRB) in the 0.5-2 keV band and existing upper limits on the 21-cm power spectrum to constrain the high-redshift population of X-ray sources, focusing on their effect on the thermal history of the Universe and the cosmic 21-cm signal. Because the properties of these sources are poorly constrained, we consider hot gas, X-ray binaries and mini-quasars (i.e., sources with soft or hard X-ray spectra) as possible candidates. We find that (1) the soft-band CXRB sets an upper limit on the X-ray efficiency of sources that existed before the end of reionization, which is one-to-two orders of magnitude higher than typically assumed efficiencies, (2) hard sources are more effective in generating the CXRB than the soft ones, (3) the commonly-assumed limit of saturated heating is not valid during the first half of reionization in the case of hard sources, with any allowed value of X-ray efficiency, (4) the maximal allowed X-ray efficiency sets a lower limit on the depth of the absorption trough in the global 21-cm signal and an upper limit on the height of the emission peak, while in the 21-cm power spectrum it sets a minimum amplitude and frequency for the high-redshift peaks, and (5) the existing upper limit on the 21-cm power spectrum sets a lower limit on the X-ray efficiency for each model. When combined with the 21-cm global signal, the CXRB will be useful for breaking degeneracies and helping constrain the nature of high-redshift heating sources.

preprint2016arXiv

The 21-cm BAO signature of enriched low-mass galaxies during cosmic reionization

Studies of the formation of the first stars have established that they formed in small halos of $\sim 10^5 - 10^6 M_{\odot}$ via molecular hydrogen cooling. Since a low level of ultraviolet radiation from stars suffices to dissociate molecular hydrogen, under the usually-assumed scenario this primordial mode of star formation ended by redshift $z \sim 15$ and much more massive halos came to dominate star formation. However, metal enrichment from the first stars may have allowed the smaller halos to continue to form stars. In this Letter we explore the possible effect of star formation in metal-rich low-mass halos on the redshifted 21-cm signal of neutral hydrogen from $z = 6-40$. These halos are significantly affected by the supersonic streaming velocity, with its characteristic baryon acoustic oscillation (BAO) signature. Thus, enrichment of low-mass galaxies can produce a strong signature in the 21-cm power spectrum over a wide range of redshifts, especially if star formation in the small halos was more efficient than suggested by current simulations. We show that upcoming radio telescopes can easily distinguish among various possible scenarios.

preprint2016arXiv

The Rise of the First Stars: Supersonic Streaming, Radiative Feedback, and 21-cm Cosmology

Understanding the formation and evolution of the first stars and galaxies represents one of the most exciting frontiers in astronomy. Since the universe was filled with hydrogen atoms at early times, the most promising probe of the epoch of the first stars is the prominent 21-cm spectral line of hydrogen. Current observational efforts are focused on the cosmic reionization era, but observations of the pre-reionization cosmic dawn are also promising. While observationally unexplored, theoretical studies predict a rich variety of observational signatures from cosmic dawn. As the first stars formed, their radiation (plus that from stellar remnants) produced significant cosmic events including Lyman-alpha coupling at z~25, and early X-ray heating. Much focus has gone to studying the angle-averaged power spectrum of 21-cm fluctuations. Additional probes include the global (sky-averaged) 21-cm spectrum, and the line-of-sight anisotropy of the 21-cm power spectrum. A particularly striking signature may result from the recently recognized effect of a supersonic relative velocity between the dark matter and gas. Work in this field, focused on understanding the whole era of reionization and cosmic dawn with analytical models and numerical simulations, is likely to grow in intensity and importance, as the theoretical predictions are finally expected to confront 21-cm observations in the coming years. [Abridged]

preprint2015arXiv

High-Redshift Star Formation in a Time-Dependent Lyman-Werner Background

The first generation of stars produces a background of Lyman-Werner (LW) radiation which can photo-dissociate molecular hydrogen, increasing the mass of dark matter halos required to host star formation. Previous studies have determined the critical mass required for efficient molecular cooling with a constant LW background. However, the true background is expected to increase rapidly at early times. Neglecting this evolution could underestimate star formation in small halos that may have started to cool in the past when the LW intensity was much lower. Background evolution is a large source of uncertainty in pre-reionization predictions of the cosmological 21cm signal, which can be observed with future radio telescopes. To address this, we perform zero-dimentional one-zone calculations that follow the density, chemical abundances, and temperature of gas in the central regions of dark matter halos, including hierarchical growth and an evolving LW background. We begin by studying the physics of halos subjected to a background that increases exponentially with redshift. We find that when the intensity increases more slowly than $J_{\rm LW}(z) \propto 10^{-z/5}$, cooling in the past is a relatively small effect. We then self-consistently compute the cosmological LW background over $z=15-50$ and find that cooling in the past due to an evolving background has a modest impact. Finally, we compare these results to three-dimensional hydrodynamical cosmological simulations with varying LW histories. While only a small number of halos were simulated, the results are consistent with our one-zone calculations.

preprint2015arXiv

Reconstructing the nature of the first cosmic sources from the anisotropic 21-cm signal

The redshifted 21-cm background is expected to be a powerful probe of the early Universe, carrying both cosmological and astrophysical information from a wide range of redshifts. In particular, the power spectrum of fluctuations in the 21-cm brightness temperature is anisotropic due to the line-of-sight velocity gradient, which in principle allows for a simple extraction of this information in the limit of linear fluctuations. However, recent numerical studies suggest that the 21-cm signal is actually rather complex, and its analysis likely depends on detailed model fitting. We present the first realistic simulation of the anisotropic 21-cm power spectrum over a wide period of early cosmic history. We show that on observable scales, the anisotropy is large and thus measurable at most redshifts, and its form tracks the evolution of 21-cm fluctuations as they are produced early on by Lyman-a radiation from stars, then switch to X-ray radiation from early heating sources, and finally to ionizing radiation from stars. In particular, we predict a redshift window during cosmic heating (at z ~ 15), when the anisotropy is small, during which the shape of the 21-cm power spectrum on large scales is determined directly by the average radial distribution of the flux from X-ray sources. This makes possible a model-independent reconstruction of the X-ray spectrum of the earliest sources of cosmic heating.

preprint2014arXiv

Detecting the First Stars and Black Holes with 21-cm Cosmology

Understanding the formation and evolution of the first stars and galaxies represents one of the most exciting frontiers in astronomy. Since the universe was filled with neutral hydrogen at early times, the most promising method for observing the epoch of the first stars is using the prominent 21-cm spectral line of the hydrogen atom. Current observational efforts are focused on the reionization era (cosmic age t~500 Myr), with earlier times considered much more challenging. However, the next frontier of even earlier galaxy formation (t~200 Myr) is emerging as a promising observational target. This is made possible by a recently noticed effect of a significant relative velocity between the baryons and dark matter at early times. The velocity difference significantly suppresses star formation. The spatial variation of this suppression enhances large-scale clustering and produces a prominent cosmic web on 100 comoving Mpc scales in the 21-cm intensity distribution. This structure makes it much more feasible for radio astronomers to detect these early stars, and should drive a new focus on this era, which is rich with little-explored astrophysics.

preprint2014arXiv

The observable signature of late heating of the Universe during cosmic reionization

Models and simulations of the epoch of reionization predict that spectra of the 21-cm transition of atomic hydrogen will show a clear fluctuation peak, at a redshift and scale, respectively, that mark the central stage of reionization and the characteristic size of ionized bubbles. This is based on the assumption that the cosmic gas was heated by stellar remnants - particularly X-ray binaries - to temperatures well above the cosmic microwave background at that time (~ 30 K). Here we show instead that the hard spectra (that is, spectra with more high-energy photons than low-energy photons) of X-ray binaries make such heating ineffective, resulting in a delayed and spatially uniform heating that modifies the 21-cm signature of reionization. Rather than looking for a simple rise and fall of the large-scale fluctuations (peaking at several millikelvin), we must expect a more complex signal also featuring a distinct minimum (at less than a millikelvin) that marks the rise of the cosmic mean gas temperature above the microwave background. Observing this signal, possibly with radio telescopes in operation today, will demonstrate the presence of a cosmic background of hard X-rays at that early time.

preprint2014arXiv

The rich complexity of 21-cm fluctuations produced by the first stars

We explore the complete history of the 21-cm signal in the redshift range z = 7-40. This redshift range includes various epochs of cosmic evolution related to primordial star formation, and should be accessible to existing or planned low-frequency radio telescopes. We use semi-numerical computational methods to explore the fluctuation signal over wavenumbers between 0.03 and 1 Mpc$^{-1}$, accounting for the inhomogeneous backgrounds of Ly-$α$, X-ray, Lyman-Werner and ionizing radiation. We focus on the recently noted expectation of heating dominated by a hard X-ray spectrum from high-mass X-ray binaries. We study the resulting delayed cosmic heating and suppression of gas temperature fluctuations, allowing for large variations in the minimum halo mass that contributes to star formation. We show that the wavenumbers at which the heating peak is detected in observations should tell us about the characteristic mean free path and spectrum of the emitted photons, thus giving key clues as to the character of the sources that heated the primordial Universe. We also consider the line-of-sight anisotropy, which allows additional information to be extracted from the 21-cm signal. For example, the heating transition at which the cosmic gas is heated to the temperature of the cosmic microwave background should be clearly marked by an especially isotropic power spectrum. More generally, an additional cross-power component $P_X$ directly probes which sources dominate 21-cm fluctuations. In particular, during cosmic reionization (and after the just-mentioned heating transition), $P_X$ is negative on scales dominated by ionization fluctuations and positive on those dominated by temperature fluctuations.

preprint2013arXiv

Complete history of the observable 21-cm signal from the first stars during the pre-reionization era

We present the first complete calculation of the history of the inhomogeneous 21-cm signal from neutral hydrogen during the era of the first stars. We use hybrid computational methods to capture the large-scale distribution of the first stars, whose radiation couples to the neutral hydrogen emission, and to evaluate the 21-cm signal from z ~ 15-35. In our realistic picture large-scale fluctuations in the 21-cm signal are sourced by the inhomogeneous density field and by the Ly-alpha and X-ray radiative backgrounds. The star formation is suppressed by two spatially varying effects: negative feedback provided by the Lyman-Werner radiative background, and supersonic relative velocities between the gas and dark matter. Our conclusions are quite promising: we find that the fluctuations imprinted by the inhomogeneous Ly-alpha background in the 21-cm signal at z ~ 25 should be detectable with the Square Kilometer Array.

preprint2013arXiv

Luminosity Bias (II): The Cosmic Web of the First Stars

Understanding the formation and evolution of the first stars and galaxies represents one of the most exciting frontiers in astronomy. Since the universe was filled with neutral hydrogen at early times, the most promising method for observing the epoch of the first stars is using the prominent 21-cm spectral line of the hydrogen atom. Current observational efforts are focused on the reionization era (cosmic age t ~ 500 Myr), with earlier times considered much more challenging. However, the next frontier of even earlier galaxy formation (t ~ 200 Myr) is emerging as a promising observational target. This is made possible by a recently noticed effect of a significant relative velocity between the baryons and dark matter at early times. The velocity difference suppresses star formation, causing a unique form of early luminosity bias. The spatial variation of this suppression enhances large-scale clustering and produces a prominent cosmic web on 100 comoving Mpc scales in the 21-cm intensity distribution. This structure makes it much more feasible for radio astronomers to detect these early stars, and should drive a new focus on this era, which is rich with little-explored astrophysics.

preprint2012arXiv

Measuring the History of Cosmic Reionization using the 21-cm Difference PDF

During cosmic reionization, the 21-cm brightness fluctuations were highly non-Gaussian, and complementary statistics can be extracted from the distribution of pixel brightness temperatures that are not derivable from the 21-cm power spectrum. One such statistic is the 21-cm difference PDF, the probability distribution function of the difference in the 21-cm brightness temperatures between two points, as a function of the distance between the points. Guided by 21-cm difference PDFs extracted from simulations, we perform a maximum likelihood analysis on mock observational data, and analyze the ability of present and future low-frequency radio array experiments to estimate the shape of the 21-cm difference PDF, and measure the history of cosmic reionization. We find that one-year data with an experiment such as the Murchison Wide-field Array should suffice for probing large scales during the mid-to-late stages of reionization, while a second-generation experiment should yield detailed measurements over a wide range of scales during most of the reionization era.

preprint2012arXiv

Profiles of Dark Matter Velocity Anisotropy in Simulated Clusters

We report statistical results for dark matter (DM) velocity anisotropy, β, from a sample of some 6000 cluster-size halos (at redshift zero) identified in a ΛCDM hydrodynamical adaptive mesh refinement simulation performed with the Enzo code. These include profiles of β in clusters with different masses, relaxation states, and at several redshifts, modeled both as spherical and triaxial DM configurations. Specifically, although we find a large scatter in the DM velocity anisotropy profiles of different halos (across elliptical shells extending to at least ~$1.5 r_{vir}$), universal patterns are found when these are averaged over halo mass, redshift, and relaxation stage. These are characterized by a very small velocity anisotropy at the halo center, increasing outward to about 0.27 and leveling off at about $0.2 r_{vir}$. Indirect measurements of the DM velocity anisotropy fall on the upper end of the theoretically expected range. Though measured indirectly, the estimations are derived by using two different surrogate measurements - X-ray and galaxy dynamics. Current estimates of the DM velocity anisotropy are based on very small cluster sample. Increasing this sample will allow testing theoretical predictions, including the speculation that the decay of DM particles results in a large velocity boost. We also find, in accord with previous works, that halos are triaxial and likely to be more prolate when unrelaxed, whereas relaxed halos are more likely to be oblate. Our analysis does not indicate that there is significant correlation (found in some previous studies) between the radial density slope, γ, and β at large radii, $0.3 r_{vir} < r < r_{vir}$.

preprint2012arXiv

Studying cosmic reionization with observations of the global 21-cm signal

We explore the ability of observations of the global brightness temperature of the 21-cm signal to constrain the reionization history and the properties of the ionizing sources. In order to describe the reionization signal, we employ either a commonly-used toy model or a structure formation model that parameterizes the properties of the ionizing sources. If the structure formation model captures the actual evolution of the reionization signal, then detecting the signal is somewhat easier than it would be for the toy model; using the toy model in this case also leads to systematic errors in reconstructing the reionization history, though a sufficiently sensitive experiment should be able to distinguish between the two models. We show that under optimistic assumptions regarding systematic noise and foreground removal, one-year observations of the global 21-cm spectrum should be able to detect a wide range of realistic models and measure the main features of the reionization history while constraining the key properties of the ionizing sources.

preprint2012arXiv

The 21-cm signature of the first stars during the Lyman-Werner feedback era

The formation of the first stars is an exciting frontier area in astronomy. Early redshifts z ~ 20 have become observationally promising as a result of a recently recognized effect of a supersonic relative velocity between the dark matter and gas. This effect produces prominent structure on 100 comoving Mpc scales, which makes it much more feasible to detect 21-cm fluctuations from the epoch of first heating. We use semi-numerical hybrid methods to follow for the first time the joint evolution of the X-ray and Lyman-Werner radiative backgrounds, including the effect of the supersonic streaming velocity on the cosmic distribution of stars. We incorporate self-consistently the negative feedback on star formation induced by the Lyman-Werner radiation, which dissociates molecular hydrogen and thus suppresses gas cooling. We find that the feedback delays the X-ray heating transition by a Delta z ~ 2, but leaves a promisingly large fluctuation signal over a broad redshift range. The large-scale power spectrum is predicted to reach a maximal signal-to-noise ratio of S/N ~ 3-4 at z ~ 18 (for a projected first-generation instrument), with S/N > 1 out to z ~ 22-23. We hope to stimulate additional numerical simulations as well as observational efforts focused on the epoch prior to cosmic reionization.

preprint2012arXiv

The signature of the first stars in atomic hydrogen at redshift 20

Dark and baryonic matter moved at different velocities in the early Universe, which strongly suppressed star formation in some regions. This was estimated to imprint a large-scale fluctuation signal of about 2 mK in the 21-cm spectral line of atomic hydrogen associated with stars at a redshift of 20, although this estimate ignored the critical contribution of gas heating due to X-rays and major enhancements of the suppression. A large velocity difference reduces the abundance of halos and requires the first stars to form in halos of about a million solar masses, substantially greater than previously expected. Here we report a simulation of the distribution of the first stars at z=20 (cosmic age of ~180 Myr), incorporating all these ingredients within a 400 Mpc box. We find that the 21-cm signature of these stars is an enhanced (10 mK) fluctuation signal on the 100-Mpc scale, characterized by a flat power spectrum with prominent baryon acoustic oscillations. The required sensitivity to see this signal is achievable with an integration time of a thousand hours with an instrument like the Murchison Wide-field Array or the Low Frequency Array but designed to operate in the range of 50-100 MHz.

preprint2012arXiv

Unusually Large Fluctuations in the Statistics of Galaxy Formation at High Redshift

We show that various milestones of high-redshift galaxy formation, such as the formation of the first stars or the complete reionization of the intergalactic medium, occurred at different times in different regions of the universe. The predicted spread in redshift, caused by large-scale fluctuations in the number density of galaxies, is at least an order of magnitude larger than previous expectations that argued for a sharp end to reionization. This cosmic scatter in the abundance of galaxies introduces new features that affect the nature of reionization and the expectations for future probes of reionization, and may help explain the present properties of dwarf galaxies in different environments. The predictions can be tested by future numerical simulations and may be verified by upcoming observations. Current simulations, limited to relatively small volumes and periodic boundary conditions, largely omit cosmic scatter and its consequences. In particular, they artificially produce a sudden end to reionization, and they underestimate the number of galaxies by up to an order of magnitude at redshift 20.

preprint2011arXiv

Impact of the Relative Motion between the Dark Matter and Baryons on the First Stars

Recently the initial supersonic relative velocity between the dark matter and baryons was shown to have an important effect on galaxy formation at high redshift. We study the impact of this relative motion on the distribution of the star-forming halos and on the formation redshift of the very first star. We include a new aspect of the relative velocity effect found in recent simulations by fitting their results to obtain the spatially-varying minimum halo mass needed for molecular cooling. Thus, the relative velocities have three separate effects: suppression of the halo abundance, suppression of the gas content within each halo, and boosting of the minimum cooling mass. We show that the two suppressions (of gas content and of halo abundance) are the primary effects on the small minihalos that cannot form stars, while the cooling mass boost combines with the abundance suppression to produce order unity fluctuations in stellar density. We quantify the large-scale inhomogeneity of galaxies, finding that 68% of the star formation (averaged on a 3 Mpc scale) is confined to 35% of the volume at z=20 (and just 18% at z=40). In addition, we estimate the redshift of the first star to be z ~ 65, which includes a delay of Dz ~ 5 due to the relative velocity.

preprint2011arXiv

The nonlinear evolution of baryonic overdensities in the early universe: Initial conditions of numerical simulations

We run very large cosmological N-body hydrodynamical simulations in order to study statistically the baryon fractions in early dark matter halos. We critically examine how differences in the initial conditions affect the gas fraction in the redshift range z = 11--21. We test three different linear power spectra for the initial conditions: (1) A complete heating model, which is our fiducial model; this model follows the evolution of overdensities correctly, according to Naoz & Barkana (2005), in particular including the spatial variation of the speed of sound of the gas due to Compton heating from the CMB. (2) An equal-δ model, which assumes that the initial baryon fluctuations are equal to those of the dark matter, while conserving sigma8 of the total matter. (3) A mean cs model, which assumes a uniform speed of sound of the gas. The latter two models are often used in the literature. We calculate the baryon fractions for a large sample of halos in our simulations. Our fiducial model implies that before reionization and significant stellar heating took place, the minimum mass needed for a minihalo to keep most of its baryons throughout its formation was ~ 3 * 10^4 Msun. However, the alternative models yield a wrong (higher by about 50%) minimum mass, since the system retains a memory of the initial conditions. We also demonstrate this using the "filtering mass" from linear theory, which accurately describes the evolution of the baryon fraction throughout the simulated redshift range.

preprint2011arXiv

Triaxiality and non-thermal gas pressure in Abell 1689

Clusters of galaxies are uniquely important cosmological probes of the evolution of the large scale structure, whose diagnostic power depends quite significantly on the ability to reliably determine their masses. Clusters are typically modeled as spherical systems whose intracluster gas is in strict hydrostatic equilibrium (i.e., the equilibrium gas pressure is provided entirely by thermal pressure), with the gravitational field dominated by dark matter, assumptions that are only rough approximations. In fact, numerical simulations indicate that galaxy clusters are typically triaxial, rather than spherical, and that turbulent gas motions (induced during hierarchical merger events) provide an appreciable pressure component. Extending our previous work, we present results of a joint analysis of X-ray, weak and strong lensing measurements of Abell 1689. The quality of the data allows us to determine both the triaxial shape of the cluster and the level of non-thermal pressure that is required if the intracluster gas is in hydrostatic equilibrium. We find that the dark matter axis ratios are 1.24 +/- 0.13 and 2.02 +/- 0.01 on the plane of the sky and along the line of sight, respectively, and that about 20% of the pressure is non-thermal. Our treatment demonstrates that the dynamical properties of clusters can be determined in a (mostly) bias-free way, enhancing the use of clusters as more precise cosmological probes.

preprint2010arXiv

Quantifying the collisionless nature of dark matter and galaxies in A1689

We use extensive measurements of the cluster A1689 to assess the expected similarity in the dynamics of galaxies and dark matter (DM) in their motion as collisionless `particles' in the cluster gravitational potential. To do so we derive the radial profile of the specific kinetic energy of the cluster galaxies from the Jeans equation and observational data. Assuming that the specific kinetic energies of galaxies and DM are roughly equal, we obtain the mean value of the DM velocity anisotropy parameter, and the DM density profile. Since this deduced profile has a scale radius that is higher than inferred from lensing observations, we tested the validity of the assumption by repeating the analysis using results of simulations for the profile of the DM velocity anisotropy. Results of both analyses indicate a significant difference between the kinematics of galaxies and DM within $r \lesssim 0.3r_{\rm vir}$. This finding is reflected also in the shape of the galaxy number density profile, which flattens markedly with respect to the steadily rising DM profile at small radii. Thus, $r \sim 0.3r_{\rm vir}$ seems to be a transition region interior to which collisional effects significantly modify the dynamical properties of the galaxy population with respect to those of DM in A1689

preprint2010arXiv

Scale-Dependent Bias of Galaxies from Baryonic Acoustic Oscillations

Baryonic acoustic oscillations (BAOs) modulate the density ratio of baryons to dark matter across large regions of the Universe. We show that the associated variation in the mass-to-light ratio of galaxies should generate an oscillatory, scale-dependent bias of galaxies relative to the underlying distribution of dark matter. A measurement of this effect would calibrate the dependence of the characteristic mass-to-light ratio of galaxies on the baryon mass fraction in their large scale environment. This bias, though, is unlikely to significantly affect measurements of BAO peak positions.

preprint2010arXiv

Statistics of 21-cm Fluctuations in Cosmic Reionization Simulations: PDFs and Difference PDFs

In the coming decade, low-frequency radio arrays will begin to probe the epoch of reionization via the redshifted 21-cm hydrogen line. Successful interpretation of these observations will require effective statistical techniques for analyzing the data. Due to the difficulty of these measurements, it is important to develop techniques beyond the standard power spectrum analysis in order to offer independent confirmation of the reionization history, probe different aspects of the topology of reionization, and have different systematic errors. In order to assess the promise of probability distribution functions (PDFs) as statistical analysis tools in 21-cm cosmology, we first measure the 21-cm brightness temperature (one-point) PDFs in six different reionization simulations. We then parametrize their most distinct features by fitting them to a simple model. Using the same simulations, we also present the first measurements of difference PDFs in simulations of reionization. We find that while these statistics probe the properties of the ionizing sources, they are relatively independent of small-scale, sub-grid astrophysics. We discuss the additional information that the difference PDF can provide on top of the power spectrum and the one-point PDF.

preprint2010arXiv

Strong-Lensing Analysis of a Complete Sample of 12 MACS Clusters at z>0.5: Mass Models and Einstein Radii

We present the results of a strong-lensing analysis of a complete sample of 12 very luminous X-ray clusters at $z>0.5$ using HST/ACS images. Our modelling technique has uncovered some of the largest known critical curves outlined by many accurately-predicted sets of multiple images. The distribution of Einstein radii has a median value of $\simeq28\arcsec$ (for a source redshift of $z_{s}\sim2$), twice as large as other lower-$z$ samples, and extends to $55\arcsec$ for MACS J0717.5+3745, with an impressive enclosed Einstein mass of $7.4\times10^{14} M_{\odot}$. We find that 9 clusters cover a very large area ($>2.5 \sq \arcmin$) of high magnification ($μ> \times10$) for a source redshift of $z_{s}\sim8$, providing primary targets for accessing the first stars and galaxies. We compare our results with theoretical predictions of the standard $Λ$CDM model which we show systematically fall short of our measured Einstein radii by a factor of $\simeq1.4$, after accounting for the effect of lensing projection. Nevertheless, a revised analysis once arc redshifts become available, and similar analyses of larger samples, are needed in order to establish more precisely the level of discrepancy with $Λ$CDM predictions.

preprint2010arXiv

Suppression and Spatial Variation of Early Galaxies and Minihalos

We study the effect of the relative velocity of dark matter and baryonic fluids after the epoch of recombination on the evolution of the first bound objects in the early universe. Recent work has shown that, although relative motion of the two fluids is formally a second order effect in density, it has a dramatic impact on the formation and distribution of the first cosmic structures. Focusing on the gas content, we analyze the effect of relative velocity on the properties of halos over a wide range of halo masses and redshifts. We calculate accurately the linear evolution of the baryon and dark matter fluctuations, and quantify the resulting effect on halos based on an analytical formalism that has been carefully checked with simulations in the case with no relative velocity. We estimate the effect on the abundance of and gas fraction in early halos. We find that the relative velocity effect causes several changes: (i) the characteristic mass that divides gas-rich and gas-poor halos is increased by roughly an order of magnitude, from 2 10^4 Msun to about 2 10^5 Msun; (ii) this characteristic mass has a large scatter (full width at half maximum is ~ 1.5 10^5 Msun at z=20); (iii) the fraction of baryons in star-less gas minihalos is suppressed by a factor of 4 at z=20; (iv) the fraction of baryons in halos that can cool and form stars is suppressed by a factor of 1.5 at z=20; and (v) there are enhanced spatial variations of these various fractions.

preprint2009arXiv

Concentrating the Dark Matter in Galaxy Clusters through Tidal Stripping of Baryonically-Compressed Galactic Halos

Gravitational lensing observations of massive X-ray clusters imply a steep characteristic density profile marked by a central concentration of dark matter. The observed mass fraction within a projected radius of 150 kpc is twice that found in state-of-the-art dark matter simulations of the standard Lambda-CDM cosmology. A central baryon enhancement that could explain this discrepancy is not observed, leaving a major puzzle. We propose a solution based on the merger histories of clusters. A significant fraction of the final dark matter content of a cluster halo originates within galaxy-sized halos, in which gas can cool and compress the dark matter core to high densities. The subsequent tidal stripping of this compressed dark matter occurs in denser regions that are closer to the center of the cluster halo. Eventually, the originally cooled gas must be dispersed into the intracluster medium through feedback, for consistency with observations that do not find central baryon enhancements in clusters. Still, the early adiabatic compression of the galactic dark matter leaves a net effect on the cluster. Using a simple model for this process, we show that the central cluster profile is substantially modified, potentially explaining the observed discrepancy.

preprint2009arXiv

Dynamical Study of A1689 from Wide-Field VLT/VIMOS Spectroscopy: Mass Profile, Concentration Parameter, and Velocity Anisotropy

We examine the dynamics structure of the rich cluster A1689, combining VLT/VIMOS spectroscopy with Subaru/Suprime-Cam imaging. The radial velocity distribution of $\sim 500$ cluster members is bounded by a pair of clearly defined velocity caustics, with a maximum amplitude of $\sim|4000|$ km/s at $\simeq$ 300 h$^{-1}$ kpc, beyond which the amplitude steadily declines, approaching zero velocity at a limiting radius of $\sim$ 2 h$^{-1}$ Mpc. We derive the 3D velocity anisotropy and galaxy number density profiles using a model-independent method to solve the Jeans equation, simultaneously incorporating the observed velocity dispersion profile, the galaxy counts from deep Subaru imaging, and our previously derived cluster mass profile from a joint lensing and X-ray analysis. The velocity anisotropy is found to be predominantly radial at large radius, becoming increasingly tangential towards the center, in accord with expectations. We also analyze the galaxy data independently of our previous analysis using two different methods: The first is based on a solution of the Jeans equation assuming an NFW form for the mass distribution, whereas in the second method the caustic amplitude is used to determine the escape velocity. The cluster virial mass derived by both of these dynamical methods is in good agreement with results from our earlier lensing and X-ray analysis. We also confirm the high NFW concentration parameter, with results from both methods combined to yield $c_{\rm vir}>13$ (1$σ$). The inferred virial radius is consistent with the limiting radius where the caustics approach zero velocity and where the counts of cluster members drop off, suggesting that infall onto A1689 is currently not significant.

preprint2009arXiv

Gas in Simulations of High Redshift Galaxies and Minihalos

We study the gas content of halos in the early universe using high resolution hydrodynamical simulations. We extract from the simulations and also predict based on linear theory the halo mass for which the enclosed baryon fraction equals half of the mean cosmic fraction. We find a rough agreement between the simulations and the predictions, which suggests that during the high-redshift era before stellar heating, the minimum mass needed for a minihalo to keep most of its baryons throughout its formation was $\sim 3 \times 10^4$ M$_\odot$. We also carry out a detailed resolution analysis and show that in order to determine a halo's gas fraction even to 20% accuracy the halo must be resolved into at least 500 dark matter particles.

preprint2009arXiv

Measuring the History of Cosmic Reionization using the 21-cm PDF from Simulations

The 21-cm PDF (i.e., distribution of pixel brightness temperatures) is expected to be highly non-Gaussian during reionization and to provide important information on the distribution of density and ionization. We measure the 21-cm PDF as a function of redshift in a large simulation of cosmic reionization and propose a simple empirical fit. Guided by the simulated PDF, we then carry out a maximum likelihood analysis of the ability of upcoming experiments to measure the shape of the 21-cm PDF and derive from it the cosmic reionization history. Under the strongest assumptions, we find that upcoming experiments can measure the reionization history in the mid to late stages of reionization to 1-10% accuracy. Under a more flexible approach that allows for four free parameters at each redshift, a similar accuracy requires the lower noise levels of second-generation 21-cm experiments.

preprint2008arXiv

Studying the Sources of Cosmic Reionization with 21-cm Fluctuations

We explore the ability of measurements of the 21-cm power spectrum during reionization to enable the simultaneous reconstruction of the reionization history and the properties of the ionizing sources. For various sets of simulated 21-cm measurements, we perform maximum likelihood fits in order to constrain the reionization and galaxy formation histories. We employ a flexible six-parameter model that parametrizes the uncertainties in the properties of high-redshift galaxies. The computational speed needed is attained through the use of an analytical model that is in reasonable agreement with numerical simulations of reionization. We find that one-year observations with the MWA array should measure the cosmic ionized fraction to ~ 1% accuracy at the very end of reionization, and a few percent accuracy around the mid-point of reionization. The mean halo mass of the ionizing sources should be measureble to 10% accuracy when reionization is 2/3 of the way through, and to 20% accuracy throughout the central stage of reionization, if this mass is anywhere in the range 1/3 to 100 billion solar masses.

preprint2005arXiv

Growth of Linear Perturbations before the Era of the First Galaxies

We calculate the evolution of linear density and temperature perturbations in a universe with dark matter, baryons, and radiation, from cosmic recombination until the epoch of the first galaxies. In addition to gravity, the perturbations are effected by electron scattering with the radiation, by radiation pressure, and by gas pressure. We include the effect of spatial fluctuations in the baryonic sound speed and show that they induce a >10% change in the baryonic density power spectrum on small scales, and a larger change on all scales in the power spectrum of gas temperature fluctuations. A precise calculation of the growth of linear perturbations is essential since they provide the initial conditions for the formation of galaxies and they can also be probed directly via cosmological 21cm fluctuations. We also show that in general the thermal history of the cosmic gas can be measured from 21cm fluctuations using a small-scale anisotropic cutoff due to the thermal width of the 21cm line.

preprint2003arXiv

Spectral Signature of Cosmological Infall of Gas Around the First Quasars

Recent observations have shown that, only a billion years after the Big Bang, the Universe was already lit up by bright quasars fuelled by the infall of gas onto supermassive black holes. The masses of these early black holes are inferred from their luminosities to be > 10^9 solar masses, which is a difficult theoretical challenge to explain. Like nearby quasars, the early objects could have formed in the central cores of massive host galaxies. The formation of these hosts could be explained if, like local large galaxies, they were assembled gravitationally inside massive (> 10^12 solar mass) halos of dark matter. There has hitherto been no observational evidence for the presence of these massive hosts or their surrounding halos. Here we show that the cosmic gas surrounding each halo must respond to its strong gravitational pull, where absorption by the infalling hydrogen produces a distinct spectral signature. That signature can be seen in recent data.

preprint2002arXiv

The Star Formation Rate Intensity Distribution Function -- Comparison of Observations with Hierarchical Galaxy Formation

Recently, Lanzetta et al. (2002) have measured the distribution of star formation rate intensity in galaxies at various redshifts. This data set has a number of advantages relative to galaxy luminosity functions; the effect of surface-brightness dimming on the selection function is simpler to understand, and this data set also probes the size distribution of galactic disks. We predict this function using semi-analytic models of hierarchical galaxy formation in a LCDM cosmology. We show that the basic trends found in the data follow naturally from the redshift evolution of dark matter halos. The data are consistent with a constant efficiency of turning gas into stars in galaxies, with a best-fit value of 2%, where dust obscuration is neglected; equivalently, the data are consistent with a cosmic star formation rate which is constant to within a factor of two at all redshifts above two. However, the practical ability to use this kind of distribution to measure the total cosmic star formation rate is limited by the predicted shape of an approximate power law with a smoothly varying power, without a sharp break.

preprint2001arXiv

In the Beginning: The First Sources of Light and the Reionization of the Universe

The formation of the first stars and quasars marks the transformation of the universe from its smooth initial state to its clumpy current state. In popular cosmological models, the first sources of light began to form at redshift 30 and reionized most of the hydrogen in the universe by redshift 7. Current observations are at the threshold of probing the hydrogen reionization epoch. The study of high-redshift sources is likely to attract major attention in observational and theoretical cosmology over the next decade.

preprint1997arXiv

Analysis of Time Delays in the Gravitational Lens PG1115+080

We present a new method for determining time delays among the light curves of various images in a gravitational lens. The method is based on constructing a simple model for the source variation and forming a chi squared measure of the agreement of this same variation with all of the lightcurves. While inspired by Press et al. (1992a, b) our approach is different since we do not assume a Gaussian process for the source variation. We apply our method to the light curves of the quadruple gravitational lens PG1115+080 measured by Schechter et al. (1997). Unlike Schechter et al. we include correlated measurement errors in the analysis, as well as the possibility that microlensing may cause different images to vary by different factors in flux. We find a value of 25.0 (+3.3,-3.8) days (95% confidence) for the delay between components B and C (close to the 24 day value of Schechter et al). However, the ratio t_AC/t_BA of the two intermediate delays is poorly determined at 1.13 (+.18,-.17) (68% confidence), close to the value predicted by lens models (1.4) unlike the Schechter et al. value (0.7). The variation ratios of C with respect to A and of A with respect to B are both different from 1, 1.39 (+.16,-.20) and .79 (+.10,-.12) (95% confidence), respectively. This is an indication of a microlensing gradient, and this type of microlensing may allow us to conclude that the size of the quasar optical emission region is about 1000 AU.