Source author record

Marcelo A. Alvarez

Marcelo A. Alvarez 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

18works
7topics
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

18 published item(s)

preprint2022arXiv

AMBER: A Semi-numerical Abundance Matching Box for the Epoch of Reionization

The Abundance Matching Box for the Epoch of Reionization (AMBER) is a semi-numerical code for modeling the cosmic dawn. The new algorithm is not based on the excursion set formalism for reionization, but takes the novel approach of calculating the reionization-redshift field $z_\mathrm{re}(\boldsymbol{x})$ assuming that hydrogen gas encountering higher radiation intensity are photoionized earlier. Redshift values are assigned while matching the abundance of ionized mass according to a given mass-weighted ionization fraction $\bar{x}_\mathrm{i}(z)$. The code has the unique advantage of allowing users to directly specify the reionization history through the redshift midpoint $z_\mathrm{mid}$, duration $Δ_\mathrm{z}$, and asymmetry $A_\mathrm{z}$ input parameters. The reionization process is further controlled through the minimum halo mass $M_\mathrm{min}$ for galaxy formation and the radiation mean free path $l_\mathrm{mfp}$ for radiative transfer. We implement improved methods for constructing density, velocity, halo, and radiation fields, which are essential components for modeling reionization observables. We compare AMBER with two other semi-numerical methods and find that our code more accurately reproduces the results from radiation-hydrodynamic simulations. The parallelized code is over four orders of magnitude faster than radiative transfer simulations and will efficiently enable large-volume models, full-sky mock observations, and parameter-space studies. AMBER will be made publicly available to facilitate and transform studies of the EoR.

preprint2022arXiv

Snowmass2021 Computational Frontier White Paper: Cosmological Simulations and Modeling

Powerful new observational facilities will come online over the next decade, enabling a number of discovery opportunities in the "Cosmic Frontier", which targets understanding of the physics of the early universe, dark matter and dark energy, and cosmological probes of fundamental physics, such as neutrino masses and modifications of Einstein gravity. Synergies between different experiments will be leveraged to present new classes of cosmic probes as well as to minimize systematic biases present in individual surveys. Success of this observational program requires actively pairing it with a well-matched state-of-the-art simulation and modeling effort. Next-generation cosmological modeling will increasingly focus on physically rich simulations able to model outputs of sky surveys spanning multiple wavebands. These simulations will have unprecedented resolution, volume coverage, and must deliver guaranteed high-fidelity results for individual surveys as well as for the cross-correlations across different surveys. The needed advances are as follows: (1) Development of scientifically rich and broadly-scoped simulations, which capture the relevant physics and correlations between probes (2) Accurate translation of simulation results into realistic image or spectral data to be directly compared with observations (3) Improved emulators and/or data-driven methods serving as surrogates for expensive simulations, constructed from a finite set of full-physics simulations (4) Detailed and transparent verification and validation programs for both simulations and analysis tools. (Abridged)

preprint2022arXiv

Snowmass2021: Opportunities from Cross-survey Analyses of Static Probes

Cosmological data in the next decade will be characterized by high-precision, multi-wavelength measurements of thousands of square degrees of the same patches of sky. By performing multi-survey analyses that harness the correlated nature of these datasets, we will gain access to new science, and increase the precision and robustness of science being pursued by each individual survey. However, effective application of such analyses requires a qualitatively new level of investment in cross-survey infrastructure, including simulations, associated modeling, coordination of data sharing, and survey strategy. The scientific gains from this new level of investment are multiplicative, as the benefits can be reaped by even present-day instruments, and can be applied to new instruments as they come online.

preprint2021arXiv

Constraining Cluster Virialization Mechanism and Cosmology using Thermal-SZ-selected clusters from Future CMB Surveys

We forecast the number of galaxy clusters that can be detected via the thermal Sunyaev-Zeldovich (tSZ) signals by future cosmic microwave background (CMB) experiments, primarily the wide area survey of the CMB-S4 experiment but also CMB-S4's smaller delensing survey and the proposed CMB-HD experiment. We predict that CMB-S4 will detect 75,000 clusters with its wide survey of $f_{\rm sky}$ = 50% and 14,000 clusters with its deep survey of $f_{\rm sky}$ = 3%. Of these, approximately 1350 clusters will be at $z \ge 2$, a regime that is difficult to probe by optical or X-ray surveys. We assume CMB-HD will survey the same sky as the S4-Wide{}, and find that CMB-HD will detect $\times3$ more overall and an order of magnitude more $z \ge 2$ clusters than CMB-S4. These results include galactic and extragalactic foregrounds along with atmospheric and instrumental noise. Using CMB-cluster lensing to calibrate cluster tSZ-mass scaling relation, we combine cluster counts with primary CMB to obtain cosmological constraints for a two parameter extension of the standard model ($ΛCDM+\sum m_ν+w_{0}$). Besides constraining $σ(w_{0})$ to $\lesssim 1\%$, we find that both surveys can enable a $\sim 2.5-4.5σ$ detection of $\sum m_ν$, substantially strengthening CMB-only constraints. We also study the evolution of intracluster medium by modelling the cluster virialization ${\rm v}(z)$ and find tight constraints from CMB-S4, with further factors of 3-4 improvement for CMB-HD. The binned cluster counts, Fisher matrices, and other associated products can be downloaded from https://github.com/sriniraghunathan/tSZ_cluster_forecasts.

preprint2021arXiv

Simulated catalogs and maps of radio galaxies at millimeter wavelengths in Websky

We present simulated millimeter-wavelength maps and catalogs of radio galaxies across the full sky that trace the nonlinear clustering and evolution of dark matter halos from the Websky simulation at $z<4.6$ and $M_{\rm halo}>10^{12} M_{\odot}/h$, and the accompanying framework for generating a new sample of radio galaxies from any halo catalog of positions, redshifts, and masses. Object fluxes are generated using a hybrid approach that combines (1) existing astrophysical halo models of radio galaxies from the literature to determine the positions and rank-ordering of the observed fluxes with (2) empirical models from the literature based on fits to the observed distribution of flux densities and (3) spectral indices drawn from an empirically-calibrated frequency-dependent distribution. The resulting population of radio galaxies is in excellent agreement with the number counts, polarization fractions, and distribution of spectral slopes from the data from observations at millimeter wavelengths from 20-200~GHz, including \emph{Planck}, ALMA, SPT, and ACT. Since the radio galaxies are correlated with the existing cosmic infrared background (CIB), Compton-$y$ (tSZ), and CMB lensing maps from Websky, our model makes new predictions for the cross-correlation power spectra and stacked profiles of radio galaxies and these other components. These simulations will be important for unbiased analysis of a wide variety of observables that are correlated with large-scale structure, such as gravitational lensing and SZ clusters.

preprint2015arXiv

Evolution of Low Mass Galactic Subhalos and Dependence on Concentration

We carry out a detailed study of the orbital dynamics and structural evolution of over 6000 subhalos in the Via Lactea II simulation, from infall to present. By analyzing subhalos with masses down to m = 4e5 Msun, we find that lower mass subhalos, which are not strongly affected by dynamical friction, exhibit behaviors qualitatively different from those found previously for more massive ones. Furthermore, there is a clear trend of subhalos that fell into the host earlier being less concentrated. We show that the concentration at infall characterizes various aspects of subhalo evolution. In particular, tidal effects truncate the growth of less concentrated subhalos at larger distances from the host; subhalos with smaller concentrations have larger infall radii. The concentration at infall is further shown to be a determining factor for the subsequent mass loss of subhalos within the host, and also for the evolution of their internal structure in the v_max-r_max plane. Our findings raise the prospects of using the concentration to predict the tidal evolution of subhalos, which will be useful for obtaining analytic models of galaxy formation, as well as for near field cosmology.

preprint2015arXiv

Probing First Galaxies and Their Impact on the Intergalactic Medium through the 21-cm Observation of the Cosmic Dawn with the SKA

We present an overview of the theory of high-redshift star and X-ray source formation, and how they affect the 21-cm background. Primary focus is given to Lyman alpha pumping and X-ray heating mechanisms at cosmic dawn, opening a new observational window for high-redshift astrophysics by generating sizable fluctuations in the 21-cm background. We describe observational prospects for power spectrum analysis and 3D tomography (imaging) of the signature of these early astrophysical sources by SKA1-LOW and SKA2.

preprint2015arXiv

Spatially Extended 21 cm Signal from Strongly Clustered UV and X-Ray Sources in the Early Universe

We present our prediction for the local 21 cm differential brightness temperature ($δT_{b}$) from a set of strongly clustered sources of Population III (Pop III) and II (Pop II) objects in the early Universe, by a numerical simulation of their formation and radiative feedback. These objects are located inside a highly biased environment, which is a rare, high-density peak ("Rarepeak") extending to $\sim7$ comoving Mpc. We study the impact of ultraviolet (UV) and X-ray photons on the intergalactic medium (IGM) and the resulting $δT_{b}$, when Pop III stars are assumed to emit X-ray photons by forming X-ray binaries very efficiently. We parameterize the rest-frame spectral energy distribution (SED) of X-ray photons, which regulates X-ray photon-trapping, IGM-heating, secondary Lyman-alpha pumping and the resulting morphology of $δT_{b}$. A combination of emission ($δT_{b}>0$) and absorption ($δT_{b}<0$) regions appears in varying amplitudes and angular scales. The boost of the signal by the high-density environment ($δ\sim0.64$) and on a relatively large scale combine to make Rarepeak a discernible, spatially-extended ($θ\sim10'$) object for 21 cm observation at $13\lesssim z\lesssim17$, which is found to be detectable as a single object by SKA with integration time of $\sim1000$ hours. Power spectrum analysis by some of the SKA precursors (LOFAR, MWA, PAPER) of such rare peaks is found difficult due to the rarity of these peaks, and the contribution only by these rare peaks to the total power spectrum remains subdominant compared to that by all astrophysical sources.

preprint2015arXiv

The Kinetic Sunyaev-Zel'dovich Effect from Reionization: Simulated Full Sky Maps at Arcminute Resolution

The kinetic Sunyaev-Zel'dovich (kSZ) effect results from Thomson scattering by coherent flows in the reionized intergalactic medium. We present new results based on ray-tracing a 10 Gpc scale simulation at 2-3 Mpc scale resolution to create a full sky kSZ map. The simulation includes, self-consistently, the effects of reionization on scales corresponding to multipoles $10\lesssim\ell\lesssim{5000}$. We separate the kSZ map into Doppler ($\mathbf{v}$), Ostriker-Vishniac ($δ\mathbf{v}$), patchy ($x\mathbf{v}$), and third-order ($xδ\mathbf{v}$) components, and compute explicitly all the auto and cross correlations (e.g., $\langle\mathbf{v}\mathbf{v}\rangle$, $\langleδ\mathbf{v}{x}\mathbf{v}\rangle$, etc.) that contribute to the total power. We find a complex dependence on the duration of reionization at $\ell\sim{300}$ and evidence for a non-negligible (10-30 per cent) contribution from connected four point ionization-velocity correlations, $\langle{x}\mathbf{v}x\mathbf{v}\rangle_c$, that are usually neglected in analytical models. We also investigate the Doppler-large scale structure (LSS) correlation, focusing on two different probes: (1) cross power spectra with linearly biased tracers of LSS and (2) cold spots from infall onto large, rare \ion{H}{2} regions centered on peaks in the matter distribution at redshifts $z>10$. Finally, the reionization history can be reconstructed at 5-10$σ$ significance by correlating 21-cm maps with existing CMB temperature maps at $\ell\!<\!500$ --sophisticated velocity reconstruction methods can probe the distribution of electrons in the IGM by using combined CMB and LSS measurements well into the epoch of reionization. The resulting kSZ maps have been made publicly available at http://www.cita.utoronto.ca/~malvarez/research/ksz-data.

preprint2015arXiv

UV driven evaporation of close-in planets: energy-limited; recombination-limited and photon-limited flows

We have investigated the evaporation of close-in exoplanets irradiated by ionizing photons. We find that the properties of the flow are controlled by the ratio of the recombination time to the flow time-scale. When the recombination time-scale is short compared to the flow time-scale the the flow is in approximate local ionization equilibrium with a thin ionization front, where the photon mean free path is short compared to flow scale. In this "recombination limited" flow the mass-loss scales roughly with the square root of the incident flux. When the recombination time is long compared to the flow time-scale the ionization front becomes thick and encompasses the entire flow, with the mass-loss rate scaling linearly with flux. If the planet's potential is deep the flow is approximately "energy-limited"; however, if the planet's potential is shallow we identify a new limiting mass-loss regime, which we term "photon-limited". In this scenario the mass-loss rate is purely limited by the incoming flux of ionizing photons. We have developed a new numerical approach that takes into account the frequency dependence of the incoming ionizing spectrum and performed a large suite of 1D simulations to characterise UV driven mass-loss around low mass planets. We find the flow is "recombination-limited" at high fluxes but becomes "energy-limited" at low fluxes; however, the transition is broad occurring over several order of magnitude in flux. Finally, we point out the transitions between the different flow types does not occur at a single flux value, but depends on the planet's properties, with higher mass planets becoming "energy-limited" at lower fluxes.

preprint2013arXiv

Reionization Histories of Milky Way Mass Halos

We investigate the connection between the epoch of reionization and the present day universe, by examining the extended mass reionization histories of dark matter halos identified at z=0. We combine an N-body dark matter simulation of a 600 Mpc volume with a three-dimensional, seminumerical reionization model. This provides reionization redshifts for each particle, which can then be connected with the properties of their halos at the present time. We find that the vast majority of present-day halos with masses larger than ~ few x 10^11 Msun reionize earlier than the rest of the universe. We also find significant halo-to-halo diversity in mass reionization histories, and find that in realistic inhomogenous models, the material within a given halo is not expected to reionize at the same time. In particular, the scatter in reionization times within individual halos is typically larger than the scatter among halos. From our fiducial reionization model, we find that the typical 68% scatter in reionization times within halos is ~ 115 Myr for 10^(12 \pm 0.25) Msun halos, decreasing slightly to ~ 95 Myr for 10^(15 \pm 0.25) Msun halos. We find a mild correlation between reionization history and environment: halos with shorter reionization histories are typically in more clustered environments, with the strongest trend on a scale of ~ 20 Mpc. Material in Milky Way mass halos with short reionization histories is preferentially reionized in relatively large HII regions, implying reionization mostly by sources external to the progenitors of the present-day halo. We investigate the impact on our results of varying the reionization model parameters, which span a range of reionization scenarios with varying timing and morphology.

preprint2013arXiv

The Opacity of the Intergalactic Medium During Reionization: Resolving Small-Scale Structure

Early in the reionization process, the intergalactic medium (IGM) would have been quite inhomogeneous on small scales, due to the low Jeans mass in the neutral IGM and the hierarchical growth of structure in a cold dark matter Universe. This small-scale structure acted as an important sink during the epoch of reionization, impeding the progress of the ionization fronts that swept out from the first sources of ionizing radiation. Here we present results of high-resolution cosmological hydrodynamics simulations that resolve the cosmological Jeans mass of the neutral IGM in representative volumes several Mpc across. The adiabatic hydrodynamics we follow are appropriate in an unheated IGM, before the gas has had a chance to respond to the photoionization heating. Our focus is determination of the resolution required in cosmological simulations in order to sufficiently sample and resolve small-scale structure regulating the opacity of an unheated IGM. We find that a dark matter particle mass of m_dm < 50 Msun and box size of L > 1 Mpc are required. With our converged results we show how the mean free path of ionizing radiation and clumping factor of ionized hydrogen depends upon the ultraviolet background (UVB) flux and redshift. We find, for example at z = 10, clumping factors typically of 10 to 20 for an ionization rate of Gamma ~ 0.3 - 3 x 1e-12 s^-1, with corresponding mean free paths of ~ 3 - 15 Mpc, extending previous work on the evolving mean free path to considerably smaller scales and earlier times.

preprint2012arXiv

Constraints on the Ionizing Efficiency of the First Galaxies

Observations of the Lyman-alpha forest and of high-redshift galaxies at z~5-10 imply that there were just enough photons to maintain the universe in an ionized state at z~5-6, indicating a "photon-starved" end to reionization. The ionizing emissivity must have been larger at earlier times in order to yield the extended reionization history implied by the electron scattering optical depth constraint from WMAP. Here we address the possibility that a faint population of galaxies with host halo masses of ~1e8-1e9 Msun dominated the ionizing photon budget at redshifts of about z>9, due to their much higher escape fractions. Such faint, early galaxies, would not have formed in ionized regions due to suppression by heating from the UV background (UVB), and would therefore not contribute to the ionizing background at z<6, after reionization is complete. Our model matches: (1) the low escape fractions observed for high-redshift galaxies, (2) the WMAP constraint of tau_es~0.09, (3) the low values for the UVB at z<6, and (4) the observed star formation rate density inferred from Lyman-break galaxies. A top heavy IMF from Pop III stars is not required in this scenario. We compare our model to recent ones in the literature that were forced to introduce an escape fraction that increases strongly towards high redshift, and show that a similar evolution occurs naturally if low mass galaxies possess high escape fractions.

preprint2012arXiv

The Effect of Absorption Systems on Cosmic Reionization

We use large-scale simulations to investigate the morphology of reionization during the final, overlap phase. Our method uses an efficient three-dimensional smoothing technique which takes into account the finite mean free path due to absorption systems, lambda, by only smoothing over scales R_s<lambda. The large dynamic range of our calculations is necessary to resolve the neutral patches left at the end of reionization within a representative volume; we find that simulation volumes exceeding several hundred Mpc on a side are necessary in order to properly model reionization when the neutral fraction is ~0.01-0.3. Our results indicate a strong dependence of percolation morphology on a large and uncertain region of model parameter space. The single most important parameter is the mean free path to absorption systems, which serve as opaque barriers to ionizing radiation. If these absorption systems were as abundant as some realistic estimates indicate, the spatial structure of the overlap phase is considerably more complex than previously predicted. In view of the lack of constraints on the mean free path at the highest redshifts, current theories that do not include absorption by Lyman-limit systems, and in particular three-dimensional simulations, may underestimate the abundance of neutral clouds at the end of reionization. This affects predictions for the 21 cm signal associated with reionization, interpretation of absorption features in quasar spectra at z ~5-6, the connection between reionization and the local universe, and constraints on the patchiness and duration of reionization from temperature fluctuations measured in the cosmic microwave background arising from the kinetic Sunyaev-Zel'dovich effect.

preprint2012arXiv

Topology and Sizes of HII Regions during Cosmic Reionization

We use the results of large-scale simulations of reionization to explore methods for characterizing the topology and sizes of HII regions during reionization. We use four independent methods for characterizing the sizes of ionized regions. Three of them give us a full size distribution: the friends-of-friends (FOF) method, the spherical average method (SPA) and the power spectrum (PS) of the ionized fraction. These latter three methods are complementary: While the FOF method captures the size distribution of the small scale H II regions, which contribute only a small amount to the total ionization fraction, the spherical average method provides a smoothed measure for the average size of the H II regions constituting the main contribution to the ionized fraction, and the power spectrum does the same while retaining more details on the size distribution. Our fourth method for characterizing the sizes of the H II regions is the average size which results if we divide the total volume of the H II regions by their total surface area, (i.e. 3V/A), computed in terms of the ratio of the corresponding Minkowski functionals of the ionized fraction field. To characterize the topology of the ionized regions, we calculate the evolution of the Euler Characteristic. We find that the evolution of the topology during the first half of reionization is consistent with inside-out reionization of a Gaussian density field. We use these techniques to investigate the dependence of size and topology on some basic source properties, such as the halo mass-to-light ratio, susceptibility of haloes to negative feedback from reionization, and the minimum halo mass for sources to form. We find that suppression of ionizing sources within ionized regions slows the growth of H II regions, and also changes their size distribution. Additionally, the topology of simulations including suppression is more complex. (abridged)

preprint2010arXiv

Enhanced Detectability of Pre-reionization 21-cm Structure

Before the universe was reionized, it was likely that the spin temperature of intergalactic hydrogen was decoupled from the CMB by UV radiation from the first stars through the Wouthuysen-Field effect. If the IGM had not yet been heated above the CMB temperature by that time, then the gas would appear in absorption relative to the CMB. Large, rare sources of X-rays could inject sufficient heat into the neutral IGM, so that the differential brightness temperature was greater than zero at comoving distances of tens to hundreds of Mpc, resulting in large 21-cm fluctuations with amplitudes of about 250 mK on arcminute to degree angular scales, an order of magnitude larger in amplitude than that caused by ionized bubbles during reionization, about 25 mK. This signal could therefore be easier to detect and probe higher redshifts than that due to patchy reionization. For the case in which the first objects to heat the IGM are QSOs hosting 10^7-solar mass black holes with an abundance exceeding about 1 per Gpc^3 at z~15, observations with either the Arecibo Observatory or the Five Hundred Meter Aperture Spherical Telescope (FAST) could detect and image their fluctuations at greater than 5-sigma significance in about a month of dedicated survey time. Additionally, existing facilities such as MWA and LOFAR could detect the statistical fluctuations arising from a population of 10^5-solar mass black holes with an abundance of about 10^4 per Gpc^3 at z~10-12.

preprint2009arXiv

The Impact of Inhomogeneous Reionization on the Satellite Galaxy Population of the Milky Way

We use the publicly available subhalo catalogs from the Via Lactea simulation along with a Gpc-scale N-body simulation to understand the impact of inhomogeneous reionization on the satellite galaxy population of the Milky Way. The large-volume simulation is combined with a model for reionization that allows us to predict the distribution of reionization times for Milky Way mass halos. Motivated by this distribution, we identify candidate satellite galaxies in the simulation by requiring that any subhalo must grow above a specified mass threshold before it is reionized; after this time the photoionizing background will suppress both the formation of stars and the accretion of gas. We show that varying the reionization time over the range expected for Milky Way mass halos can change the number of satellite galaxies by roughly two orders of magnitude. This conclusion is in contradiction with a number of studies in the literature, and we conclude that this is a result of inconsistent application of the results of Gnedin (2000). We compare our satellite galaxies to observations using both abundance matching and stellar population synthesis methods to assign luminosities to our subhalos and account for observational completeness effects. Additionally, if we assume that the mass threshold is set by the virial temperature Tvir = 8e3K we find that our model accurately matches the vmax distribution, radial distribution, and luminosity function of observed Milky Way satellites for a reionization time zreion = 9.6^{1.0}_{-2.1}, assuming that the Via Lacteasubhalo distribution is representative of the Milky Way. This results in the presence of 119^{+202}_{-50} satellite galaxies.

preprint2006arXiv

Recognizing the First Radiation Sources Through Their 21-cm Signature

At the beginning of the reionization epoch, radiation sources produce fluctuations in the redshifted 21-cm background. We show that different types of sources (such as miniquasars, Pop II and III stars, supernovae, etc.) produce distinct signatures in the 21-cm signal radial profiles and statistical fluctuations, through which they can be identified. Further, we show that the 21-cm signal from X-ray emitting sources is much easier to observe than was expected, due to a previously neglected pumping mechanism.