Source author record

Anson D'Aloisio

Anson D'Aloisio 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

17works
2topics
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

17 published item(s)

preprint2021arXiv

Hydrodynamic Response of the Intergalactic Medium to Reionization

The intergalactic medium is expected to clump on scales down to $10^4-10^8$ M$_{\odot}$ before the onset of reionization. The impact of these small-scale structures on reionization is poorly understood despite the modern understanding that gas clumpiness limits the growth of H II regions. We use a suite of radiation-hydrodynamics simulations that capture the $\sim 10^4$ $M_\odot$ Jeans mass of unheated gas to study density fluctuations during reionization. Our simulations track the complex ionization and hydrodynamical response of gas in the wake of ionization fronts. The clumping factor of ionized gas (proportional to the recombination rate) rises to a peak value of $5-20$ approximately $Δt = 10$ Myr after ionization front passage, depending on the incident intensity, redshift, and degree to which the gas had been pre-heated by the first X-ray sources. The clumping factor reaches its relaxed value of $\approx 3$ by $Δt = 300$ Myr. The mean free path of Lyman-limit photons evolves in unison, being up to several times shorter in un-relaxed, recently reionized regions compared to those that were reionized much earlier. Assessing the impact of this response on the global reionizaton process, we find that un-relaxed gaseous structures boost the total number of recombinations by $\approx 50$ % and lead to spatial fluctuations in the mean free path that persist appreciably for several hundred million years after the completion of reionization.

preprint2021arXiv

Hydrodynamic Response of the Intergalactic Medium to Reionization II: Physical Characteristics and Dynamics of Ionizing Photon Sinks

Becker et al. 2021 measured the mean free path of Lyman limit photons in the IGM at $z=6$. The short value suggests that absorptions may have played a prominent role in reionization. Here we study physical properties of ionizing photon sinks in the wake of ionization fronts (I-fronts) using radiative hydrodynamic simulations. We quantify the contributions of gaseous structures to the Lyman limit opacity by tracking the column density distributions in our simulations. Within $Δt = 10$ Myr of I-front passage, we find that self-shielding systems ($N_{\rm HI} > 10^{17.2}$ cm$^{-2}$) are comprised of two distinct populations: (1) over-density $Δ\sim 50$ structures in photo-ionization equilibrium with the ionizing background; (2) $Δ\gtrsim 100$ density peaks with fully neutral cores. The self-shielding systems contribute more than half of the opacity at these times, but the IGM evolves considerably in $Δt \sim 100$ Myr as structures are flattened by pressure smoothing and photoevaporation. By $Δt = 300$ Myr, they contribute $\lesssim 10 \%$ to the opacity in an average 1 Mpc$^3$ patch of the Universe. The percentage can be a factor of a few larger in over-dense patches, where more self-shielding systems survive. We quantify the characteristic masses and sizes of self-shielding structures. Shortly after I-front passage, we find $M=10^{4} - 10^8$ M$_\odot$ and effective diameters $d_{\rm eff} = 1 - 20$ ckpc$/h$. These scales increase as the gas relaxes. The picture herein presented may be different in dark matter models with suppressed small-scale power.

preprint2020arXiv

A Model-Insensitive Baryon Acoustic Oscillation Feature in the 21 cm Signal from Reionization

We examine the impact of baryon-dark matter relative velocities on intergalactic small-scale structure and the 21 cm signal during reionization. Streaming velocities reduced clumping in the intergalactic medium (IGM) on mass scales of $\sim 10^4 - 10^8$ M$_{\odot}$. This effect produced a distinct baryon acoustic oscillation (BAO) feature in the 21 cm power spectrum at wave numbers $k\sim 0.1$ h/Mpc, near which forthcoming surveys will be most sensitive. In contrast to the highly uncertain impact of streaming velocities on star formation, the effect on clumping is better constrained because it is set mainly by cosmology and straightforward gas dynamics. We quantify the latter using coupled radiation-hydrodynamic simulations that capture the Jeans scale of pre-reionization gas. The clumping factor of ionized gas is reduced by 5-10\% in regions with RMS streaming velocities. The suppression peaks $\approx 5$ Myr after a region is reionized, but disappears within 200 Myr due to pressure smoothing. We model the corresponding impact on the 21 cm signal and find that the BAO feature is most likely to appear at $\approx$ 10 \% ionization. During this phase, the feature may appear at the 1 \% (5 \%) level at $k \sim 0.1 (0.06)$ h/Mpc with an amplitude that varies by a factor of $< 10$ across a range of reionization histories. We also provide a model for the signal originating from streaming velocity's impact on ionizing sources, which can vary by 4 orders of magnitude depending on highly uncertain source properties. We find that the clumping signal probably dominates the source one unless Population III star formation in $10^6 - 10^8$ M$_{\odot}$ halos contributed significantly to the first 10\% of reionization.

preprint2020arXiv

Observing the tail of reionization: neutral islands in the $z=5.5$ Lyman-$α$ Forest

Previous studies have noted difficulties in modeling the highest opacities of the $z>5.5$ Ly$α$ forest, epitomized by the extreme Ly$α$ trough observed towards quasar ULAS J0148+0600. One possibility is that the most opaque regions at these redshifts contain significant amounts of neutral hydrogen. This explanation, which abandons the common assumption that reionization ended before $z=6$, also reconciles evidence from independent observations of a significantly neutral Universe at $z=7.5$. Here we explore a model in which the neutral fraction is still $\approx 10\%$ at $z=5.5$. We confirm that this model can account for the observed scatter in Ly$α$ forest opacities, as well as the observed Ly$β$ transmission in the J0148 trough. We contrast the model with a competing "earlier" reionization scenario characterized by a short mean free path and large fluctuations in the post-reionization ionizing background. We consider Ly$α$ and Ly$β$ effective optical depths, their correlations, trough size distributions, dark pixel fractions, the IGM thermal history, and spatial distributions of Lyman-$α$ emitters around forest sight lines. We find that the models are broadly similar in almost all of these statistics, suggesting that it may be difficult to distinguish between them definitively. We argue that improved constraints on the mean free path and the thermal history at $z>5$ could go a long way towards diagnosing the origin of the $z>5.5$ opacity fluctuations.

preprint2016arXiv

Cosmic Dawn (CoDa): the First Radiation-Hydrodynamics Simulation of Reionization and Galaxy Formation in the Local Universe

Cosmic reionization by starlight from early galaxies affected their evolution, thereby impacting reionization, itself. Star formation suppression, for example, may explain the observed underabundance of Local Group dwarfs relative to N-body predictions for Cold Dark Matter. Reionization modelling requires simulating volumes large enough [~(100Mpc)^3] to sample reionization "patchiness", while resolving millions of galaxy sources above ~10^8 Msun , combining gravitational and gas dynamics with radiative transfer. Modelling the Local Group requires initial cosmological density fluctuations pre-selected to form the well-known structures of the local universe today. Cosmic Dawn ("CoDa") is the first such fully-coupled, radiation-hydrodynamics simulation of reionization of the local universe. Our new hybrid CPU-GPU code, RAMSES-CUDATON, performs hundreds of radiative transfer and ionization rate-solver timesteps on the GPUs for each hydro-gravity timestep on the CPUs. CoDa simulated (91Mpc)^3 with 4096^3 particles and cells, to redshift 4.23, on ORNL supercomputer Titan, utilizing 8192 cores and 8192 GPUs. Global reionization ended slightly later than observed. However, a simple temporal rescaling which brings the evolution of ionized fraction into agreement with observations also reconciles ionizing flux density, cosmic star formation history, CMB electron scattering optical depth and galaxy UV luminosity function with their observed values. Photoionization heating suppressed the star formation of haloes below ~2 x 10^9 Msun , decreasing the abun- dance of faint galaxies around MAB_1600 = [-10,-12]. For most of reionization, star formation was dominated by haloes between 10^10 - 10^11 Msun , so low-mass halo suppression was not reflected by a distinct feature in the global star formation history. (Abridged)

preprint2016arXiv

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.

preprint2015arXiv

Large Opacity Variations in the High-Redshift Lyman-alpha Forest: the Signature of Relic Temperature Fluctuations from Patchy Reionization

Recent observations of the Lyman-alpha forest show large-scale spatial variations in the intergalactic Lyman-alpha opacity that grow rapidly with redshift at z>5, far in excess of expectations from empirically motivated models. Previous studies have attempted to explain this excess with spatial fluctuations in the ionizing background, but found that this required either extremely rare sources or problematically low values for the mean free path of ionizing photons. Here we report that much -- or potentially all -- of the observed excess likely arises from residual spatial variations in temperature that are an inevitable byproduct of a patchy and extended reionization process. The amplitude of opacity fluctuations generated in this way depends on the timing and duration of reionization. If the entire excess is due to temperature variations alone, the observed fluctuation amplitude favors a late-ending but extended reionization process that was roughly half complete by z~9 and that ended at z~6. In this scenario, the highest opacities occur in regions that reionized earliest, since they have had the most time to cool, while the lowest opacities occur in the warmer regions that reionized most recently. This correspondence potentially opens a new observational window into patchy reionization.

preprint2015arXiv

The Linear Perturbation Theory of Reionization in Position-Space: Cosmological Radiative Transfer Along the Light-Cone

The linear perturbation theory of inhomogeneous reionization (LPTR) has been developed as an analytical tool for predicting the global ionized fraction and large-scale power spectrum of ionized density fluctuations during reionization. In the original formulation of the LPTR, the ionization balance and radiative transfer equations are linearized and solved in Fourier space. However, the LPTR's approximation to the full solution of the radiative transfer equation is not straightforward to interpret, since the latter is most intuitively conceptualized in position space. To bridge the gap between the LPTR and the language of numerical radiative transfer, we present a new, equivalent, position-space formulation of the LPTR that clarifies the approximations it makes and facilitates its interpretation. We offer a comparison between the LPTR and the excursion-set model of reionization (ESMR), and demonstrate the built-in capability of the LPTR to explore a wide range of reionization scenarios, and to go beyond the ESMR in exploring scenarios involving X-rays.

preprint2014arXiv

The Effect of Large-Scale Structure on the Magnification of High-Redshift Sources by Cluster-Lenses

Cluster gravitational lensing surveys like the Hubble Space Telescope Frontier Fields survey will detect distant galaxies 10-50 times fainter than any yet discovered. Using these surveys to measure the luminosity function of such faint, distant galaxies, however, requires that magnification maps built from the constraints of strongly-lensed images be accurate. For models that assume the cluster and nearby (correlated) structures are the only significant sources of lensing, a potential source of error in these maps comes from the fact that light rays also suffer weak deflections by uncorrelated large-scale structure along the line-of-sight, i.e. cosmic weak lensing (CWL). To demonstrate the magnitude of this effect, we calculate the magnification change which results when the same cluster-lens is placed along different lines of sight. Using a simple density profile for a cluster-lens at z~0.3-0.5 and the power spectrum of the matter density fluctuations responsible for CWL, we show that the typical magnifications of ~5(10) of sources at z=6-10 can differ by ~10-20(20-30)% from one line-of-sight to another. However, these fluctuations rise to greater than order unity near critical curves, indicating that CWL tends to make its greatest contribution to the most magnified images. We conclude that the neglect of CWL in determining the intrinsic luminosities of highly-magnified galaxies may introduce errors significant enough to warrant further effort to include this contribution in cluster-lens modeling. We suggest that methods of modeling CWL in galaxy-strong-lensing systems should be generalized to cluster-lensing systems.

preprint2013arXiv

Primordial Non-Gaussianity Estimation using 21 cm Tomography from the Epoch of Reionization

Measuring the small primordial nonGaussianity (PNG) predicted by cosmic inflation theories may help diagnose them. The detectability of PNG by its imprint on the 21cm power spectrum from the epoch of reionization is reassessed here in terms of $f_{NL}$, the local nonlinearity parameter. We find that an optimum, multi-frequency observation by SKA can achieve $Δf_{NL} \sim 3$ (comparable to recent Planck CMB limits), while a cosmic-variance-limited array of this size like Omniscope can even detect $Δf_{NL} \sim 0.2$. This substantially revises the methods and results of previous work.

preprint2013arXiv

The scale-dependent signature of primordial non-Gaussianity in the large-scale structure of cosmic reionization

(ABRIDGED)The rise of cosmic structure depends upon the statistical distribution of initial density fluctuations generated by inflation. While the simplest models predict an almost perfectly Gaussian distribution, more-general models predict a level of primordial non-Gaussianity (PNG) that observations might yet be sensitive enough to detect. Recent Planck Collaboration measurements of the CMB temperature anisotropy bispectrum significantly tighten the observational limits, but they are still far from the PNG level predicted by the simplest models of inflation. Probing levels below CMB sensitivities will require other methods, such as searching for the statistical imprint of PNG on galactic halo clustering. During the epoch of reionization (EoR), the first stars and galaxies released radiation into the intergalactic medium (IGM) that created ionized patches whose large-scale geometry and evolution reflected the underlying abundance and large-scale clustering of the star-forming galaxies. This statistical connection between ionized patches in the IGM and galactic halos suggests that observing reionization may be another way to constrain PNG. We employ the linear perturbation theory of reionization and semi-analytic models based on the excursion-set formalism to model the effects of PNG on the EoR. We quantify the effects of PNG on the large-scale structure of reionization by deriving the ionized density bias, i.e. ratio of ionized atomic to total matter overdensities in Fourier space, at small wavenumber. Just as previous studies found that PNG creates a scale-dependent signature in the halo bias, so, too, we find a scale-dependent signature in the ionized density bias. Our results, which differ significantly from previous attempts in the literature to characterize this PNG signature, will be applied elsewhere to predict its observable consequences, e.g. in the cosmic 21cm background.

preprint2012arXiv

Halo statistics in non-Gaussian cosmologies: the collapsed fraction, conditional mass function, and halo bias from the path-integral excursion set method

Characterizing the level of primordial non-Gaussianity (PNG) in the initial conditions for structure formation is one of the most promising ways to test inflation and differentiate among different scenarios. The scale-dependent imprint of PNG on the large-scale clustering of galaxies and quasars has already been used to place significant constraints on the level of PNG in our observed Universe. Such measurements depend upon an accurate and robust theory of how PNG affects the bias of galactic halos relative to the underlying matter density field. We improve upon previous work by employing a more general analytical method - the path-integral extension of the excursion set formalism - which is able to account for the non-Markovianity caused by PNG in the random-walk model used to identify halos in the initial density field. This non-Markovianity encodes information about environmental effects on halo formation which have so far not been taken into account in analytical bias calculations. We compute both scale-dependent and -independent corrections to the halo bias, along the way presenting an expression for the conditional collapsed fraction for the first time, and a new expression for the conditional halo mass function. To leading order in our perturbative calculation, we recover the halo bias results of Desjacques et. al. (2011), including the new scale-dependent correction reported there. However, we show that the non-Markovian dynamics from PNG can lead to marked differences in halo bias when next-to-leading order terms are included. We quantify these differences here. [abridged]

preprint2012arXiv

The Effects of Primordial Non-Gaussianity on Giant-Arc Statistics: A Scale Dependent Example

In a recently published article, we quantified the impact of primordial non-Gaussianity on the probability of giant-arc formation. In that work, we focused on the local form of non-Gaussianity and found that it can have only a modest effect given the most recent constraints from Cosmic Microwave Background (CMB) measurements. Here, we present new calculations using a parameterization of scale-dependent non-Gaussianity in which the primordial bispectrum has the equilateral shape and the effective f_NL parameter depends on scale. We find that non-Gaussianity of this type can yield a larger effect on the giant-arc abundance compared to the local form due to both the scale dependence and the relatively weaker constraints on the equilateral shape from CMB measurements. In contrast to the maximum ~40% effect (within the latest CMB constraints) previously found for the local form, we find that the predicted giant-arc abundance for the scale-dependent equilateral form can differ by a factor of a few with respect to the Gaussian case.

preprint2011arXiv

Cosmography with cluster strong lenses: the influence of substructure and line-of-sight halos

We explore the use of strong lensing by galaxy clusters to constrain the dark energy equation of state and its possible time variation. The cores of massive clusters often contain several multiply imaged systems of background galaxies at different redshifts. The locations of lensed images can be used to constrain cosmological parameters due to their dependence on the ratio of angular diameter distances. We employ Monte-Carlo simulations of cluster lenses, including the contribution from substructures, to assess the feasibility of this potentially powerful technique. At the present, parametric lens models use well motivated scaling relations between mass and light to incorporate cluster member galaxies, and do not explicitly model line-of-sight structure. Here, we quantify modeling errors due to scatter in the cluster galaxy scaling relations and un-modeled line-of-sight halos. These errors are of the order of a few arcseconds on average for clusters located at typical redshifts (z ~ 0.2 - 0.3). Using Bayesian Markov Chain Monte-Carlo techniques, we show that the inclusion of these modeling errors is critical to deriving unbiased constraints on dark energy. However, when the uncertainties are properly quantified, we show that constraints competitive with other methods may be obtained by combining results from a sample of just 10 simulated clusters with 20 families each. Cosmography with a set of well studied cluster lenses may provide a powerful complementary probe of the dark energy equation of state. Our simulations provide a convenient method of quantifying modeling errors and assessing future strong lensing survey strategies.

preprint2011arXiv

The effects of primordial non-Gaussianity on giant-arc statistics

For over a decade, it has been debated whether the concordance LCDM model is consistent with the observed abundance of giant arcs in clusters. While previous theoretical studies have focused on properties of the lens and source populations, as well as cosmological effects such as dark energy, the impact of initial conditions on the giant-arc abundance is relatively unexplored. Here, we quantify the impact of non-Gaussian initial conditions with the local bispectrum shape on the predicted frequency of giant arcs. Using a path-integral formulation of the excursion set formalism, we extend a semi-analytic model for calculating halo concentrations to the case of primordial non-Gaussianity, which may be useful for applications outside of this work. We find that massive halos tend to collapse earlier in models with positive f_NL, relative to the Gaussian case, leading to enhanced concentration parameters. The converse is true for f_NL < 0. In addition to these effects, which change the lensing cross sections, non-Gaussianity also modifies the abundance of supercritical clusters available for lensing. These combined effects work together to either enhance (f_NL > 0) or suppress (f_NL < 0) the probability of giant-arc formation. Using the best value and 95% confidence levels currently available from the Wilkinson Microwave Anisotropy Probe, we find that the giant-arc optical depth for sources at z_s~2 is enhanced by ~20% and ~45% for f_NL = 32 and 74 respectively. In contrast, we calculate a suppression of ~5% for f_NL = -10. These differences translate to similar relative changes in the predicted all-sky number of giant arcs.

preprint2011arXiv

Virialization of high redshift dark matter haloes

We present results of a study of the virial state of high redshift dark matter haloes in an N-body simulation. We find that the majority of collapsed, bound haloes are not virialized at any redshift slice in our study ($z=15-6$) and have excess kinetic energy. At these redshifts, merging is still rampant and the haloes cannot strictly be treated as isolated systems. To assess if this excess kinetic energy arises from the environment, we include the surface pressure term in the virial equation explicitly and relax the assumption that the density at the halo boundary is zero. Upon inclusion of the surface term, we find that the haloes are much closer to virialization, however, they still have some excess kinetic energy. We report trends of the virial ratio including the extra surface term with three key halo properties: spin, environment, and concentration. We find that haloes with closer neighbors depart more from virialization, and that haloes with larger spin parameters do as well. We conclude that except at the lowest masses ($M < 10^6 \Msun$), dark matter haloes at high redshift are not fully virialized. This finding has interesting implications for galaxy formation at these high redshifts, as the excess kinetic energy will impact the subsequent collapse of baryons and the formation of the first disks and/or baryonic structures.

preprint2010arXiv

Cosmological Constraints from Strong Gravitational Lensing in Galaxy Clusters

Current efforts in observational cosmology are focused on characterizing the mass-energy content of the Universe. We present results from a geometric test based on strong lensing in galaxy clusters. Based on Hubble Space Telescope images and extensive ground-based spectroscopic follow-up of the massive galaxy cluster Abell 1689, we used a parametric model to simultaneously constrain the cluster mass distribution and dark energy equation of state. Combining our cosmological constraints with those from X-ray clusters and the Wilkinson Microwave Anisotropy Probe 5-year data gives Ωm = 0.25 +/- 0.05 and wx = -0.97 +/- 0.07 which are consistent with results from other methods. Inclusion of our method with all other techniques available brings down the current 2σ contours on the dark energy equation of state parameter wx by about 30%.