Source author record

Jose Oñorbe

Jose Oñorbe 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

8works
4topics
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

8 published item(s)

preprint2022arXiv

Forecasts for WEAVE-QSO: 3D clustering and connectivity of critical points with Lyman-$α$ tomography

The upcoming WEAVE-QSO survey will target a high density of quasars over a large area, enabling the reconstruction of the 3D density field through Lyman-$α$ tomography over unprecedented volumes smoothed on intermediate scales ($\approx$ 16 Mpc/$h$). We produce mocks of the Lyman-$α$ forest using LyMAS, and reconstruct the 3D density field between sightlines through Wiener filtering in a configuration compatible with the future WEAVE-QSO observations. The fidelity of the reconstruction is assessed by measuring one- and two-point statistics from the distribution of critical points in the cosmic web. In addition, initial Lagrangian statistics are predicted from first principles, and measurements of the connectivity of the cosmic web are performed. The reconstruction captures well the expected features in the auto- and cross-correlations of the critical points. This remains true after a realistic noise is added to the synthetic spectra, even though sparsity of sightlines introduces systematics, especially in the cross-correlations of points with mixed signature. Specifically, for walls and filaments, the most striking clustering features could be measured with up to 4 sigma of significance with a WEAVE-QSO-like survey. Moreover, the connectivity of each peak identified in the reconstructed field is globally consistent with its counterpart in the original field, indicating that the reconstruction preserves the geometry of the density field not only statistically, but also locally. Hence the critical points relative positions within the tomographic reconstruction could be used as standard rulers for dark energy by WEAVE-QSO and similar surveys.

preprint2022arXiv

Measuring the thermal and ionization state of the low-$z$ IGM using likelihood free inference

We present a new approach to measure the power-law temperature density relationship $T=T_0 (ρ/ \barρ)^{γ-1}$ and the UV background photoionization rate $Γ_{\rm HI}$ of the IGM based on the Voigt profile decomposition of the Ly$α$ forest into a set of discrete absorption lines with Doppler parameter $b$ and the neutral hydrogen column density $N_{\rm HI}$. Previous work demonstrated that the shape of the $b$-$N_{\rm HI}$ distribution is sensitive to the IGM thermal parameters $T_0$ and $γ$, whereas our new inference algorithm also takes into account the normalization of the distribution, i.e. the line-density d$N$/d$z$, and we demonstrate that precise constraints can also be obtained on $Γ_{\rm HI}$. We use density-estimation likelihood-free inference (DELFI) to emulate the dependence of the $b$-$N_{\rm HI}$ distribution on IGM parameters trained on an ensemble of 624 Nyx hydrodynamical simulations at $z = 0.1$, which we combine with a Gaussian process emulator of the normalization. To demonstrate the efficacy of this approach, we generate hundreds of realizations of realistic mock HST/COS datasets, each comprising 34 quasar sightlines, and forward model the noise and resolution to match the real data. We use this large ensemble of mocks to extensively test our inference and empirically demonstrate that our posterior distributions are robust. Our analysis shows that by applying our new approach to existing Ly$α$ forest spectra at $z\simeq 0.1$, one can measure the thermal and ionization state of the IGM with very high precision ($σ_{\log T_0} \sim 0.08$ dex, $σ_γ\sim 0.06$, and $σ_{\log Γ_{\rm HI}} \sim 0.07$ dex).

preprint2015arXiv

Characterizing the Pressure Smoothing Scale of the Intergalactic Medium

The thermal state of the intergalactic medium (IGM) at z < 6 constrains the nature and timing of cosmic reionization events, but its inference from the Ly-alpha forest is degenerate with the 3-D structure of the IGM on ~100 kpc scales, where, analogous to the classical Jeans argument, the pressure of the T~$10^4$ K gas supports it against gravity. We simulate the IGM using smoothed particle hydrodynamics, and find that, at z < 6, the gas density power spectrum does not exhibit the expected Jeans filtering cutoff, because dense gas in collapsed halos dominates the small-scale power masking pressure smoothing effects. We introduce a new statistic, the real-space Ly-alpha flux, $F_\mathrm{real}$, which naturally suppresses dense gas, and is thus robust against the poorly understood physics of galaxy formation, revealing pressure smoothing in the diffuse IGM. The $F_\mathrm{real}$ power spectrum is accurately described by a simple fitting function with cutoff at $λ_F$, allowing us to rigorously quantify the filtering scale for the first time: we find $λ_F$ = 79 kpc (comoving) at z=3 for our fiducial thermal model. This statistic has the added advantage that it directly relates to observations of correlated Ly-alpha forest absorption in close quasar pairs, recently proposed as a method to measure the filtering scale. Our results enable one to quantify the filtering scale in simulations, and ask meaningful questions about its dependence on reionization and thermal history. Accordingly, the standard description of the IGM in terms of the amplitude $T_0$ and slope $γ$ of the temperature-density relation $T = T_0Δ^{γ-1}$ should be augmented with a third filtering scale parameter $λ_F$.

preprint2015arXiv

Effects of Coupled Dark Energy on the Milky Way and its Satellites

We present the first numerical simulations in coupled dark energy cosmologies with high enough resolution to investigate the effects of the coupling on galactic and sub-galactic scales. We choose two constant couplings and a time-varying coupling function and we run simulations of three Milky-Way-size halos ($\sim$10$^{12}$M$_{\odot}$), a lower mass halo (6$\times$10$^{11}$M$_{\odot}$) and a dwarf galaxy halo (5$\times$10$^{9}$M$_{\odot}$). We resolve each halo with several millions dark matter particles. On all scales the coupling causes lower halo concentrations and a reduced number of substructures with respect to LCDM. We show that the reduced concentrations are not due to different formation times, but they are related to the extra terms that appear in the equations describing the gravitational dynamics. On the scale of the Milky Way satellites, we show that the lower concentrations can help in reconciling observed and simulated rotation curves, but the coupling values necessary to have a significant difference from LCDM are outside the current observational constraints. On the other hand, if other modifications to the standard model allowing a higher coupling (e.g. massive neutrinos) are considered, coupled dark energy can become an interesting scenario to alleviate the small-scale issues of the LCDM model.

preprint2015arXiv

Forged in FIRE: cusps, cores, and baryons in low-mass dwarf galaxies

We present ultra-high resolution cosmological hydrodynamic simulations of $M_*\simeq10^{4-6.3}M_{\odot}$ dwarf galaxies that form within $M_{v}=10^{9.5-10}M_{\odot}$ dark matter halos. Our simulations rely on the FIRE implementation of star formation feedback and were run with high enough force and mass resolution to directly resolve stellar and dark matter structure on the ~200 pc scales of interest for classical and ultra-faint dwarfs in the Local Group. The resultant galaxies sit on the $M_*$ vs. $M_{v}$ relation required to match the Local Group stellar mass function. They have bursty star formation histories and also form with half-light radii and metallicities that broadly match those observed for local dwarfs at the same stellar mass. We demonstrate that it is possible to create a large (~1 kpc) dark matter core in a cosmological simulation of an $M_*\simeq10^{6.5}M_{\odot}$ dwarf galaxy that resides within an $M_{v}=10^{10}M_{\odot}$ halo -- precisely the scale of interest for resolving the Too Big to Fail problem. However, these large cores are not ubiquitous and appear to correlate closely with the star formation histories of the dwarfs: dark matter cores are largest in systems that form their stars late ($z\lesssim2$), after the early epoch of cusp building mergers has ended. Our $M_*\simeq10^4M_{\odot}$ dwarf retains a cuspy dark matter halo density profile that matches almost identically that of a dark-matter only run of the same system. Despite forming in a field environment, this very low mass dwarf has observable properties that match closely to those of ultra-faint satellite galaxies of the Milky Way, including a uniformly old stellar population (>10 Gyr). Though ancient, most of the stars in our ultra-faint form after reionization; the UV field acts mainly to suppress fresh gas accretion, not to boil away gas that is already present in the proto-dwarf.

preprint2014arXiv

Core Formation in Dwarf Halos with Self Interacting Dark Matter: No Fine-Tuning Necessary

We investigate the effect of self-interacting dark matter (SIDM) on the density profiles of $V_{\rm max} \simeq 40~km~s^{-1}$ isolated dwarf dark matter halos -- the scale of relevance for the too big to fail problem (TBTF) -- using very high-resolution cosmological zoom simulations. Each halo has millions of particles within its virial radius. We find that SIDM models with cross sections per unit mass spanning the range σ/m = $0.5 - 50$ $cm^2~g^{-1}$ alleviate TBTF and produce constant density cores of size 300-1000 pc, comparable to the half-light radii of $M_\star$ ~ $10^{5-7}$ $M_\odot$ dwarfs. The largest, lowest density cores develop for cross sections in the middle of this range, σ/m ~ $5-10~cm^2~g^{-1}$. Our largest SIDM cross section run (σ/m = $50~cm^2~g^{-1}$) develops a slightly denser core owing to mild core-collapse behavior, but it remains less dense than the CDM case and retains a constant density core profile. Our work suggests that SIDM cross sections as large or larger than $50~cm^2~g^{-1}$ remain viable on velocity scales of dwarf galaxies ($v_{\rm rms}$ ~ $40~km~s^{-1}$). The range of SIDM cross sections that alleviate TBTF and the cusp/core problem spans at least two orders of magnitude and therefore need not be particularly fine-tuned.

preprint2014arXiv

The High-$z$ Universe Confronts Warm Dark Matter: Galaxy Counts, Reionization and the Nature of Dark Matter

We use $N$-body simulations to show that high-redshift galaxy counts provide an interesting constraint on the nature of dark matter, specifically Warm Dark Matter (WDM), owing to the lack of early structure formation these models. Our simulations include three WDM models with thermal-production masses of 0.8 keV, 1.3 keV, and 2.6 keV, as well as CDM. Assuming a relationship between dark halo mass and galaxy luminosity that is set by the observed luminosity function at bright magnitudes, we find that 0.8 keV WDM is disfavored by direct galaxy counts in the Hubble Ultra Deep Field at $>\!\!10σ$. Similarly, 1.3 keV WDM is statistically inconsistent at $2.2σ$. Future observations with JWST (and possibly HST via the Frontier Fields) could rule out $1.3$ keV WDM at high significance, and may be sensitive to WDM masses greater than 2.6 keV. We also examine the ability of galaxies in these WDM models to reionize the universe, and find that 0.8 keV and 1.3 keV WDM produce optical depths to the Cosmic Microwave Background (CMB) that are inconsistent at 68% C.L. with current Planck results, even with extremely high ionizing radiation escape fractions, and 2.6 keV WDM requires an optimistic escape fraction to yield an optical depth consistent with Planck data. Although CMB optical depth calculations are model dependent, we find a strong challenge for stellar processes alone to reionize the universe in a 0.8 keV and 1.3 keV WDM cosmology.

preprint2013arXiv

How to Zoom: Bias, Contamination, and Lagrange Volumes in Multimass Cosmological Simulations

We perform a suite of multimass cosmological zoom simulations of individual dark matter halos and explore how to best select Lagrangian regions for resimulation without contaminating the halo of interest with low-resolution particles. Such contamination can lead to significant errors in the gas distribution of hydrodynamical simulations, as we show. For a fixed Lagrange volume, we find that the chance of contamination increases systematically with the level of zoom. In order to avoid contamination, the Lagrangian volume selected for resimulation must increase monotonically with the resolution difference between parent box and the zoom region. We provide a simple formula for selecting Lagrangian regions (in units of the halo virial volume) as a function of the level of zoom required. We also explore the degree to which a halo's Lagrangian volume correlates with other halo properties (concentration, spin, formation time, shape, etc.) and find no significant correlation. There is a mild correlation between Lagrange volume and environment, such that halos living in the most clustered regions have larger Lagrangian volumes. Nevertheless, selecting halos to be isolated is not the best way to ensure inexpensive zoom simulations. We explain how one can safely choose halos with the smallest Lagrangian volumes, which are the least expensive to resimulate, without biasing one's sample.