Researcher profile

Marcelo A. Alvarez

Marcelo A. Alvarez contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
10works
0followers
5topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

10 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.

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

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&#39;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)

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.