Showing 1–50 of 55 results for author: Onorbe, J

The impact of varying inhomogeneous reionization histories on metrics of Ly$α$ opacity

Authors: Caitlin C. Doughty, Joseph F. Hennawi, Jose Oñorbe, Frederick B. Davies, Zarija Lukić

Abstract: The epoch of hydrogen reionization is complete by $z=5$, but its progression at higher redshifts is uncertain. Measurements of Ly$α$ forest opacity show large scatter at $z<6$, suggestive of spatial fluctuations in neutral fraction ($x_\mathrm{HI}$), temperature, or ionizing background, either individually or in combination. However, these effects are degenerate, necessitating modeling these physi… ▽ More The epoch of hydrogen reionization is complete by $z=5$, but its progression at higher redshifts is uncertain. Measurements of Ly$α$ forest opacity show large scatter at $z<6$, suggestive of spatial fluctuations in neutral fraction ($x_\mathrm{HI}$), temperature, or ionizing background, either individually or in combination. However, these effects are degenerate, necessitating modeling these physics in tandem in order to properly interpret the observations. We begin this process by developing a framework for modeling the reionization history and associated temperature fluctuations, with the intention of incorporating ionizing background fluctuations at a later time. To do this, we generate several reionization histories using semi-numerical code AMBER, selecting histories with volume-weighted neutral fractions that adhere to the observed CMB optical depth and dark pixel fractions. Implementing these histories in the \texttt{Nyx} cosmological hydrodynamics code, we examine the evolution of gas within the simulation, and the associated metrics of the Ly$α$ forest opacity. We find that the pressure smoothing scale within the IGM is strongly correlated with the adiabatic index of the temperature-density relation. We find that while models with 20,000 K photoheating at reionization are better able to reproduce the shape of the observed $z=5$ 1D flux power spectrum than those with 10,000 K, they fail to match the highest wavenumbers. The simulated autocorrelation function and optical depth distributions are systematically low and narrow, respectively, compared to the observed values, but are in better agreement when the reionization history is longer in duration, more symmetric in its distribution of reionization redshifts, or if there are remaining neutral regions at $z<6$. The systematically low variance likely requires the addition of a fluctuating UVB. △ Less

Submitted 3 December, 2023; originally announced December 2023.

Comments: To be submitted to MNRAS. Comments welcome

Measurements of the Thermal and Ionization State of the Intergalactic Medium during the Cosmic Afternoon

Authors: Teng Hu, Vikram Khaire, Joseph F. Hennawi, Todd M. Tripp, Jose Oñorbe, Michael Walther, Zarija Lukic

Abstract: We perform the first measurement of the thermal and ionization state of the intergalactic medium (IGM) across 0.9 < z < 1.5 using 301 \lya absorption lines fitted from 12 HST STIS quasar spectra, with a total pathlength of Δz=2.1. We employ the machine-learning-based inference method that uses joint b-N distributions obtained from \lyaf decomposition. Our results show that the HI photoionization r… ▽ More We perform the first measurement of the thermal and ionization state of the intergalactic medium (IGM) across 0.9 < z < 1.5 using 301 \lya absorption lines fitted from 12 HST STIS quasar spectra, with a total pathlength of Δz=2.1. We employ the machine-learning-based inference method that uses joint b-N distributions obtained from \lyaf decomposition. Our results show that the HI photoionization rates, Γ, are in good agreement with the recent UV background synthesis models, with \log (Γ/s^{-1})={-11.79}^{0.18}_{-0.15}, -11.98}^{0.09}_{-0.09}, and {-12.32}^{0.10}_{-0.12} at z=1.4, 1.2, and 1 respectively. We obtain the IGM temperature at the mean density, T_0, and the adiabatic index, γ, as [\log (T_0/K), γ]= [{4.13}^{+0.12}_{-0.10}, {1.34}^{+0.10}_{-0.15}], [{3.79}^{+0.11}_{-0.11}, {1.70}^{+0.09}_{-0.09}] and [{4.12}^{+0.15}_{-0.25}, {1.34}^{+0.21}_{-0.26}] at z=1.4, 1.2 and 1 respectively. Our measurements of T_0 at z=1.4 and 1.2 are consistent with the expected trend from z<3 temperature measurements as well as theoretical expectations that, in the absence of any non-standard heating, the IGM should cool down after HeII reionization. Whereas, our T_0 measurements at z=1 show unexpectedly high IGM temperature. However, because of the relatively large uncertainty in these measurements of the order of ΔT_0~5000 K, mostly emanating from the limited redshift path length of available data in these bins, we can not definitively conclude whether the IGM cools down at z<1.5. Lastly, we generate a mock dataset to test the constraining power of future measurement with larger datasets. The results demonstrate that, with redshift pathlength Δz \sim 2 for each redshift bin, three times the current dataset, we can constrain the T_0 of IGM within 1500K. Such precision would be sufficient to conclusively constrain the history of IGM thermal evolution at z < 1.5. △ Less

Submitted 2 February, 2024; v1 submitted 29 November, 2023; originally announced November 2023.

Forecasting constraints on the high-z IGM thermal state from the Lyman-$α$ forest flux auto-correlation function

Authors: Molly Wolfson, Joseph F. Hennawi, Frederick B. Davies, Zarija Lukić, Jose Oñorbe

Abstract: The auto-correlation function of the Lyman-$α$ (Ly$α$) forest flux from high-z quasars can statistically probe all scales of the intergalactic medium (IGM) just after the epoch of reionization. The thermal state of the IGM, which is determined by the physics of reionization, sets the amount of small-scale power seen in the \lya forest. To study the sensitivity of the auto-correlation function to t… ▽ More The auto-correlation function of the Lyman-$α$ (Ly$α$) forest flux from high-z quasars can statistically probe all scales of the intergalactic medium (IGM) just after the epoch of reionization. The thermal state of the IGM, which is determined by the physics of reionization, sets the amount of small-scale power seen in the \lya forest. To study the sensitivity of the auto-correlation function to the thermal state of the IGM, we compute the auto-correlation function from cosmological hydrodynamical simulations with semi-numerical models of the thermal state of the IGM. We create mock data sets of 20 quasars to forecast constraints on $T_0$ and $γ$, which characterize a tight temperature-density relation in the IGM, at $5.4 \leq z \leq 6$. At $z = 5.4$ we find that an ideal data set constrains $T_0$ to 29\% and $γ$ to 9\%. In addition, we investigate four realistic reionization scenarios that combine temperature and ultra-violet background (UVB) fluctuations at $z = 5.8$. We find that, when using mock data generated from a model that includes temperature and UVB fluctuations, we can rule out a model with no temperature or UVB fluctuations at $>1σ$ level 50.5\% of the time. △ Less

Submitted 11 September, 2023; originally announced September 2023.

Comments: 18 pages, 13 figures, comments welcome. arXiv admin note: text overlap with arXiv:2208.09013

The Impact of the WHIM on the IGM Thermal State Determined from the Low-$z$ Lyman-$α$ Forest

Authors: Teng Hu, Vikram Khaire, Joseph F. Hennawi, Jose Onorbe, Michael Walther, Zarija Lukic, Frederick Davies

Abstract: At $z \lesssim 1$, shock heating caused by large-scale velocity flows and possibly violent feedback from galaxy formation, converts a significant fraction of the cool gas ($T\sim 10^4$ K) in the intergalactic medium (IGM) into warm-hot phase (WHIM) with $T >10^5$K, resulting in a significant deviation from the previously tight power-law IGM temperature-density relationship,… ▽ More At $z \lesssim 1$, shock heating caused by large-scale velocity flows and possibly violent feedback from galaxy formation, converts a significant fraction of the cool gas ($T\sim 10^4$ K) in the intergalactic medium (IGM) into warm-hot phase (WHIM) with $T >10^5$K, resulting in a significant deviation from the previously tight power-law IGM temperature-density relationship, $T=T_0 (ρ/ {\barρ})^{γ-1}$. This study explores the impact of the WHIM on measurements of the low-$z$ IGM thermal state, $[T_0,γ]$, based on the $b$-$N_{H I}$ distribution of the Lyman-$α$ forest. Exploiting a machine learning-enabled simulation-based inference method trained on Nyx hydrodynamical simulations, we demonstrate that [$T_0$, $γ$] can still be reliably measured from the $b$-$N_{H I}$ distribution at $z=0.1$, notwithstanding the substantial WHIM in the IGM. To investigate the effects of different feedback, we apply this inference methodology to mock spectra derived from the IllustrisTNG and Illustris simulations at $z=0.1$. The results suggest that the underlying $[T_0,γ]$ of both simulations can be recovered with biases as low as $|Δ\log(T_0/\text{K})| \lesssim 0.05$ dex, $|Δγ| \lesssim 0.1$, smaller than the precision of a typical measurement. Given the large differences in the volume-weighted WHIM fractions between the three simulations (Illustris 38\%, IllustrisTNG 10\%, Nyx 4\%) we conclude that the $b$-$N_{H I}$ distribution is not sensitive to the WHIM under realistic conditions. Finally, we investigate the physical properties of the detectable Lyman-$α$ absorbers, and discover that although their $T$ and $Δ$ distributions remain mostly unaffected by feedback, they are correlated with the photoionization rate used in the simulation. △ Less

Submitted 28 August, 2023; originally announced August 2023.

Comments: 23 pages, 24 figures

First measurement of the Mg II forest correlation function in the Epoch of Reionization

Authors: Suk Sien Tie, Joseph F. Hennawi, Feige Wang, Silvia Onorato, Jinyi Yang, Eduardo Bañados, Frederick B. Davies, Jose Oñorbe

Abstract: In the process of producing the roughly three ionizing photons per atom required to reionize the IGM, the same massive stars explode and eject metals into their surroundings. While the overly sensitive Lya transition makes Gunn-Peterson absorption of background quasar light an ineffective probe of reionization at z > 6, strong low-ionization transitions like the MgII doublet will give rise to a de… ▽ More In the process of producing the roughly three ionizing photons per atom required to reionize the IGM, the same massive stars explode and eject metals into their surroundings. While the overly sensitive Lya transition makes Gunn-Peterson absorption of background quasar light an ineffective probe of reionization at z > 6, strong low-ionization transitions like the MgII doublet will give rise to a detectable "metal-line forest", if metals pollute the neutral IGM. We measure the auto-correlation of the MgII forest transmission using a sample of ten ground based z >= 6.80 quasar spectra probing the redshift range 5.96 < z_MgII < 7.42 (z_MgII,median = 6.47). The correlation function exhibits strong small-scale clustering and a pronounced peak at the doublet velocity (768 km/s) arising from strong absorbers in the CGM of galaxies. After these strong absorbers are identified and masked the signal is consistent with noise. Our measurements are compared to a suite of models generated by combining a large hydrodynamical simulation with a semi-numerical reionization topology, assuming a simple uniform enrichment model. We obtain a 95% credibility upper limit of [Mg/H] < -3.73 at z_MgII,median = 6.47, assuming uninformative priors on [Mg/H] and the IGM neutral fraction x_HI. Splitting the data into low-z (5.96 < z_MgII < 6.47; z_MgII,median = 6.235) and high-z (6.47 < z_MgII < 7.42; z_MgII,median = 6.72) subsamples again yields null-detections and 95% upper limits of [Mg/H] < -3.75 and [Mg/H] < -3.45, respectively. These first measurements set the stage for an approved JWST Cycle 2 program (GO 3526) targeting a similar number of quasars that will be an order of magnitude more sensitive, making the Mgii forest an emerging powerful tool to deliver precision constraints on the reionization and enrichment history of the Universe. △ Less

Submitted 22 August, 2023; originally announced August 2023.

Comments: Submitted to MNRAS; 23 pages, 21 figures, 1 table

Convergence of small scale Ly$α$ structure at high-$z$ under different reionization scenarios

Authors: Caitlin C. Doughty, Joseph F. Hennawi, Frederick B. Davies, Zarija Lukić, Jose Oñorbe

Abstract: The Ly$α$ forest (LAF) at $z>5$ probes the thermal and reionization history of the intergalactic medium (IGM) and the nature of dark matter, but its interpretation requires comparison to cosmological hydrodynamical simulations. At high-$z$, convergence of these simulations is more exacting since transmission is dominated by underdense voids that are challenging to resolve. With evidence mounting f… ▽ More The Ly$α$ forest (LAF) at $z>5$ probes the thermal and reionization history of the intergalactic medium (IGM) and the nature of dark matter, but its interpretation requires comparison to cosmological hydrodynamical simulations. At high-$z$, convergence of these simulations is more exacting since transmission is dominated by underdense voids that are challenging to resolve. With evidence mounting for a late end to reionization, small structures down to the sub-kpc level may survive to later times than conventionally thought due to the reduced time for pressure smoothing to impact the gas, further tightening simulation resolution requirements. We perform a suite of simulations using the Eulerian cosmological hydrodynamics code Nyx, spanning domain sizes of 1.25-10 $h^{-1}$ Mpc and 5-80 $h^{-1}$ kpc cells, and explore the interaction of these variables with the timing of reionization on the properties of the matter distribution and the simulated LAF at $z=5.5$. In observable Ly$α$ power, convergence within 10% is achieved for $k< 0.1$ s/km, but larger $k$ shows deviation of up to 20 percent. While a later reionization retains more small structure in the density field, because of the greater thermal broadening there is little difference in the convergence of LAF power between early ($z=9$) and later ($z=6$) reionizations. We conclude that at $z\sim5.5$, resolutions of 10 kpc are necessary for convergence of LAF power at $k<0.1$ s/km, while higher-$k$ modes require higher resolution, and that the timing of reionization does not significantly impact convergence given realistic photoheating. △ Less

Submitted 25 May, 2023; originally announced May 2023.

Comments: Submitted to MNRAS. Comments welcome

Modelling the cosmological Lyman-Werner background radiation field in the Early Universe

Authors: Andrea Incatasciato, Sadegh Khochfar, Jose Oñorbe

Abstract: The Lyman-Werner (LW) radiation field is a key ingredient in the chemo-thermal evolution of gas in the Early Universe, as it dissociates H2 molecules, the primary cooling channel in an environment devoid of metals and dust. Despite its important role, it is still not implemented in cosmological simulations on a regular basis, in contrast to the ionising UV background. This is in part due to uncert… ▽ More The Lyman-Werner (LW) radiation field is a key ingredient in the chemo-thermal evolution of gas in the Early Universe, as it dissociates H2 molecules, the primary cooling channel in an environment devoid of metals and dust. Despite its important role, it is still not implemented in cosmological simulations on a regular basis, in contrast to the ionising UV background. This is in part due to uncertainty in the source modelling, their spectra and abundance, as well as the detailed physics involved in the propagation of the photons and their interactions with the molecules. The goal of this work is to produce an accurate model of the LW radiation field at $z\geq6$, by post-processing the physics-rich high-resolution FiBY simulation. Our novelties include updated cross sections for H$_2$, H$^-$ and H$^+_2$ chemical species, IGM absorption by neutral Hydrogen and various spectral models for Population III and Population II stars. With our fiducial set of parameters, we show that the mean LW intensity steadily increases by three orders of magnitude from $z\sim23$ to $z\sim6$, while spatial inhomogeneities originate from massive star-forming galaxies that dominate the photon budget up to a distance of $\sim100$ proper kpc. Our model can be easily applied to other simulations or semi-analytical models as an external radiation field that regulates the formation of stars and massive black hole seeds in high-$z$ low-mass halos. △ Less

Submitted 11 April, 2023; v1 submitted 19 January, 2023; originally announced January 2023.

Comments: 20 pages, 16 figures, plus 4 figures in the appendices. Main result in Figure 16. Published by MNRAS

The wide-field, multiplexed, spectroscopic facility WEAVE: Survey design, overview, and simulated implementation

Authors: Shoko Jin, Scott C. Trager, Gavin B. Dalton, J. Alfonso L. Aguerri, J. E. Drew, Jesús Falcón-Barroso, Boris T. Gänsicke, Vanessa Hill, Angela Iovino, Matthew M. Pieri, Bianca M. Poggianti, D. J. B. Smith, Antonella Vallenari, Don Carlos Abrams, David S. Aguado, Teresa Antoja, Alfonso Aragón-Salamanca, Yago Ascasibar, Carine Babusiaux, Marc Balcells, R. Barrena, Giuseppina Battaglia, Vasily Belokurov, Thomas Bensby, Piercarlo Bonifacio , et al. (190 additional authors not shown)

Abstract: WEAVE, the new wide-field, massively multiplexed spectroscopic survey facility for the William Herschel Telescope, will see first light in late 2022. WEAVE comprises a new 2-degree field-of-view prime-focus corrector system, a nearly 1000-multiplex fibre positioner, 20 individually deployable 'mini' integral field units (IFUs), and a single large IFU. These fibre systems feed a dual-beam spectrogr… ▽ More WEAVE, the new wide-field, massively multiplexed spectroscopic survey facility for the William Herschel Telescope, will see first light in late 2022. WEAVE comprises a new 2-degree field-of-view prime-focus corrector system, a nearly 1000-multiplex fibre positioner, 20 individually deployable 'mini' integral field units (IFUs), and a single large IFU. These fibre systems feed a dual-beam spectrograph covering the wavelength range 366$-$959\,nm at $R\sim5000$, or two shorter ranges at $R\sim20\,000$. After summarising the design and implementation of WEAVE and its data systems, we present the organisation, science drivers and design of a five- to seven-year programme of eight individual surveys to: (i) study our Galaxy's origins by completing Gaia's phase-space information, providing metallicities to its limiting magnitude for $\sim$3 million stars and detailed abundances for $\sim1.5$ million brighter field and open-cluster stars; (ii) survey $\sim0.4$ million Galactic-plane OBA stars, young stellar objects and nearby gas to understand the evolution of young stars and their environments; (iii) perform an extensive spectral survey of white dwarfs; (iv) survey $\sim400$ neutral-hydrogen-selected galaxies with the IFUs; (v) study properties and kinematics of stellar populations and ionised gas in $z<0.5$ cluster galaxies; (vi) survey stellar populations and kinematics in $\sim25\,000$ field galaxies at $0.3\lesssim z \lesssim 0.7$; (vii) study the cosmic evolution of accretion and star formation using $>1$ million spectra of LOFAR-selected radio sources; (viii) trace structures using intergalactic/circumgalactic gas at $z>2$. Finally, we describe the WEAVE Operational Rehearsals using the WEAVE Simulator. △ Less

Submitted 31 October, 2023; v1 submitted 7 December, 2022; originally announced December 2022.

Comments: 41 pages, 27 figures, accepted for publication by MNRAS; updated version including information on individual grants in a revised Acknowledgements section, corrections to the affiliation list, and an updated references list

Forecasting constraints on the mean free path of ionizing photons at $z \geq 5.4$ from the Lyman-$α$ forest flux auto-correlation function

Authors: Molly Wolfson, Joseph F. Hennawi, Frederick B. Davies, Jose Oñorbe

Abstract: Fluctuations in Lyman-$α$ (Ly$α$) forest transmission towards high-$z$ quasars are partially sourced from spatial fluctuations in the ultraviolet background (UVB), the level of which are set by the mean free path of ionizing photons ($λ_{\text{mfp}}$). The auto-correlation function of Ly$α$ forest flux characterizes the strength and scale of transmission fluctuations and, as we show, is thus sensi… ▽ More Fluctuations in Lyman-$α$ (Ly$α$) forest transmission towards high-$z$ quasars are partially sourced from spatial fluctuations in the ultraviolet background (UVB), the level of which are set by the mean free path of ionizing photons ($λ_{\text{mfp}}$). The auto-correlation function of Ly$α$ forest flux characterizes the strength and scale of transmission fluctuations and, as we show, is thus sensitive to $λ_{\text{mfp}}$. Recent measurements at $z \sim 6$ suggest a rapid evolution of $λ_{\text{mfp}}$ at $z>5.0$ which would leave a signature in the evolution of the auto-correlation function. For this forecast, we model mock Ly$α$ forest data with properties similar to the XQR-30 extended data set at $5.4 \leq z \leq 6.0$. At each $z$ we investigate 100 mock data sets and an ideal case where mock data matches model values of the auto-correlation function. For ideal data with $λ_{\text{mfp}}=9.0$ cMpc at $z=6.0$, we recover $λ_{\text{mfp}}=12^{+6}_{-3}$ cMpc. This precision is comparable to direct measurements of $λ_{\text{mfp}}$ from the stacking of quasar spectra beyond the Lyman limit. Hypothetical high-resolution data leads to a $\sim40\%$ reduction in the error bars over all $z$. The distribution of mock values of the auto-correlation function in this work is highly non-Gaussian for high-$z$, which should caution work with other statistics of the high-$z$ Ly$α$ forest against making this assumption. We use a rigorous statistical method to pass an inference test, however future work on non-Gaussian methods will enable higher precision measurements. △ Less

Submitted 3 March, 2023; v1 submitted 18 August, 2022; originally announced August 2022.

Comments: Accepted for publication in MNRAS, 16 pages, 11 figures

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

Authors: Teng Hu, Vikram Khaire, Joseph F. Hennawi, Michael Walther, Hector Hiss, Justin Alsing, Jose Oñorbe, Zarija Lukic, Frederick Davies

Abstract: 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… ▽ More 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). △ Less

Submitted 14 July, 2022; originally announced July 2022.

Comments: 20 pages, 17 figures, accepted for publication in MNRAS

The Cosmic Mach Number as an Environment Measure for the Underlying Dark Matter Density Field

Authors: Romain Meriot, Sadegh Khochfar, Jose Onorbe, Britton Smith

Abstract: Using cosmological dark matter only simulations of a $(1.6$ Gpc$/h)^3$ volume from the Legacy simulation project, we calculate Cosmic Mach Numbers (CMN) and perform a theoretical investigation of their relation with halo properties and features of the density field to gauge their use as an measure of the environment. CMNs calculated on individual spheres show correlations with both the overdensi… ▽ More Using cosmological dark matter only simulations of a $(1.6$ Gpc$/h)^3$ volume from the Legacy simulation project, we calculate Cosmic Mach Numbers (CMN) and perform a theoretical investigation of their relation with halo properties and features of the density field to gauge their use as an measure of the environment. CMNs calculated on individual spheres show correlations with both the overdensity in a region and the density gradient in the direction of the bulk flow around that region. To reduce the scatter around the median of these correlations, we introduce a new measure, the rank ordered Cosmic Mach number ($\hat{\mathcal{M}}_g$), which shows a tight correlations with the overdensity $δ=\frac{ρ-\barρ}{\barρ}$. Measures of the large scale density gradient as well as other average properties of the halo population in a region show tight correlations with $\hat{\mathcal{M}}_g$ as well. Our results in this first empirical study suggest that $\hat{\mathcal{M}}_g$ is an excellent proxy for the underlying density field and hence environment that can circumvent reliance on number density counts in a region. For scales between $10$ and $100 Mpc$/h, Mach numbers calculated using dark matter halos $(> 10^{12}$ M$_{\odot})$ that would typically host massive galaxies are consistent with theoretical predictions of the linear matter power spectrum at a level of $10\%$ due to non-linear effects of gravity. At redshifts $z\geq 3$, these deviations disappear. We also quantify errors due to missing large scale modes in simulations. Simulations of box size $\leq 1 $ Gpc/$h$ typically predict CMNs 10-30\% too small on scales of$\sim 100$ Mpc$/h$. △ Less

Submitted 13 February, 2022; v1 submitted 8 February, 2022; originally announced February 2022.

Comments: 16 pages, 20 figures. Accepted for publication in MNRAS on 08/02/2022. Typos corrected

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

Authors: Katarina Kraljic, Clotilde Laigle, Christophe Pichon, Sebastien Peirani, Sandrine Codis, Junsup Shim, Corentin Cadiou, Dmitri Pogosyan, Stéphane Arnouts, Matthiew Pieri, Vid Iršič, Sean S. Morrison, Jose Oñorbe, Ignasi Pérez-Ràfols, Gavin Dalton

Abstract: 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 configura… ▽ More 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. △ Less

Submitted 7 January, 2022; originally announced January 2022.

Comments: 26 pages, 22 figures, submitted to MNRAS

Improving IGM temperature constraints using wavelet analysis on high-redshift quasars

Authors: Molly Wolfson, Joseph F. Hennawi, Frederick B. Davies, Jose Oñorbe, Hector Hiss, Zarija Lukić

Abstract: The thermal state of the intergalactic medium (IGM) contains vital information about the epoch of reionization, one of the most transformative yet poorly understood periods in the young universe. This thermal state is encoded in the small-scale structure of Lyman-$α$ (Ly$α$) absorption in quasar spectra. The 1D flux power spectrum measures the average small-scale structure along quasar sightlines.… ▽ More The thermal state of the intergalactic medium (IGM) contains vital information about the epoch of reionization, one of the most transformative yet poorly understood periods in the young universe. This thermal state is encoded in the small-scale structure of Lyman-$α$ (Ly$α$) absorption in quasar spectra. The 1D flux power spectrum measures the average small-scale structure along quasar sightlines. At high redshifts, where the opacity is large, averaging mixes high signal-to-noise ratio transmission spikes with noisy absorption troughs. Wavelet amplitudes are an alternate statistic that maintains spatial information while quantifying fluctuations at the same spatial frequencies as the power spectrum, giving them the potential to more sensitively measure the small-scale structure. Previous Ly$α$ forest studies using wavelet amplitude probability density functions (PDFs) used limited spatial frequencies and neglected strong correlations between PDF bins and across wavelets scales, resulting in sub-optimal and unreliable parameter inference. Here we present a novel method for performing statistical inference using wavelet amplitude PDFs that spans the full range of spatial frequencies probed by the power spectrum and that fully accounts for these correlations. We applied this procedure to realistic mock data drawn from a simple thermal model parameterized by the temperature at mean density, $T_0$, and find that wavelets deliver 1$σ$ constraints on $T_0$ that are on average 7% more sensitive at $z=5$ (12% at $z=6$) than those from the power spectrum. We consider the possibility of combing wavelet PDFs with the power, but find that this does not lead to improved sensitivity. △ Less

Submitted 6 October, 2021; originally announced October 2021.

Comments: Accepted for publication in MNRAS, 22 Pages, 16 Figures

Gaseous Dynamical Friction: A Challenge to Modern Hydrodynamical Schemes

Authors: Ben Morton, Sadegh Khochfar, Jose Oñorbe

Abstract: The process of momentum and energy transfer from a massive body moving through a background medium, known as dynamical friction (DF), is key to our understanding of many astrophysical systems. We present a series of high-resolution simulations of gaseous DF using Lagrangian hydrodynamics solvers, in the state-of-the-art multi-physics code, GIZMO. The numerical setup is chosen to allow direct compa… ▽ More The process of momentum and energy transfer from a massive body moving through a background medium, known as dynamical friction (DF), is key to our understanding of many astrophysical systems. We present a series of high-resolution simulations of gaseous DF using Lagrangian hydrodynamics solvers, in the state-of-the-art multi-physics code, GIZMO. The numerical setup is chosen to allow direct comparison to analytic predictions for DF in the range of Mach 0.2<M<3. We investigate, in detail, the DF drag force, the radial structure of the wake, and their time evolution. The subsonic forces are shown to be well resolved, except at Mach numbers close to M=1. The supersonic cases, close to M=1, significantly under-shoot the predicted force. We find that for scenarios with 0.7<M<2, between 10%-50% of the expected DF force is missing. The origin of this deficit is mostly related to the wake structure close to the perturber, where the density profile of the Mach cone face shows significant smoothing, which does not improve with time. The spherical expanding perturbation of the medium is captured well by the hydro scheme, but it is the sharp density structure, at the transition from Mach cone to average density, that introduces the mismatch. However, we find a general improvement of the force deficit with time, though significant differencesremain, in agreement with other studies. This is due to (1) the structure of the far field wake being reproduced well, and (2) the fraction of total drag from the far field wake increasing with time. Dark matter sub-haloes, in typical cosmological simulations, occupy parameters similar to those tested here, suggesting that the DF which these sub-haloes experience is significantly underestimated, and hence their merger rate. Dynamical friction is a relevant benchmark and should be included as one of the standard hydro tests for astrophysical simulations. △ Less

Submitted 29 March, 2021; originally announced March 2021.

Comments: 16 pages, 18 figures, submitted to MNRAS

Introducing SPHINX-MHD: The Impact of Primordial Magnetic Fields on the First Galaxies, Reionization, and the Global 21cm Signal

Authors: Harley Katz, Sergio Martin-Alvarez, Joakim Rosdahl, Taysun Kimm, Jérémy Blaizot, Martin G. Haehnelt, Léo Michel-Dansac, Thibault Garel, Jose Oñorbe, Julien Devriendt, Adrianne Slyz, Omar Attia, Romain Teyssier

Abstract: We present the first results from SPHINX-MHD, a suite of cosmological radiation-magnetohydrodynamics simulations designed to study the impact of primordial magnetic fields (PMFs) on galaxy formation and the evolution of the intergalactic medium during the epoch of reionization. The simulations are among the first to employ multi-frequency, on-the-fly radiation transfer and constrained transport id… ▽ More We present the first results from SPHINX-MHD, a suite of cosmological radiation-magnetohydrodynamics simulations designed to study the impact of primordial magnetic fields (PMFs) on galaxy formation and the evolution of the intergalactic medium during the epoch of reionization. The simulations are among the first to employ multi-frequency, on-the-fly radiation transfer and constrained transport ideal MHD in a cosmological context to simultaneously model the inhomogeneous process of reionization as well as the growth of PMFs. We run a series of $(5\,\text{cMpc})^3$ cosmological volumes, varying both the strength of the seed magnetic field ($B_0$) and its spectral index ($n_B$). We find that PMFs that have $n_B > -0.562\log_{10}\left(\frac{B_0}{1{\rm n}G}\right) - 3.35$ produce electron optical depths ($τ_e$) that are inconsistent with CMB constraints due to the unrealistically early collapse of low-mass dwarf galaxies. For $n_B\geq-2.9$, our constraints are considerably tighter than the $\sim{\rm n}G$ constraints from Planck. PMFs that do not satisfy our constraints have little impact on the reionization history or the shape of the UV luminosity function. Likewise, detecting changes in the Lya forest due to PMFs will be challenging because photoionisation and photoheating efficiently smooth the density field. However, we find that the first absorption feature in the global 21cm signal is a sensitive indicator of the properties of the PMFs, even for those that satisfy our $τ_e$ constraint. Furthermore, strong PMFs can marginally increase the escape of LyC photons by up to 25\% and shrink the effective radii of galaxies by $\sim44\%$ which could increase the completeness fraction of galaxy surveys. Finally, our simulations show that surveys with a magnitude limit of ${\rm M_{UV,1500}=-13}$ can probe the sources that provide the majority of photons for reionization out to $z=12$. △ Less

Submitted 22 July, 2021; v1 submitted 27 January, 2021; originally announced January 2021.

Comments: 29 pages, 27 figures, accepted MNRAS

Probing Reionization and Early Cosmic Enrichment with the MgII Forest

Authors: Joseph F. Hennawi, Frederick B. Davies, Feige Wang, Jose Oñorbe

Abstract: Because the same massive stars that reionized the intergalactic medium (IGM) inevitably exploded as supernovae that polluted the Universe with metals, the history of cosmic reionization and enrichment are intimately intertwined. While the overly sensitive Ly-alpha transition completely saturates in a neutral IGM, strong low-ionization metal lines like the MgII 2796,2804 doublet will give rise to a… ▽ More Because the same massive stars that reionized the intergalactic medium (IGM) inevitably exploded as supernovae that polluted the Universe with metals, the history of cosmic reionization and enrichment are intimately intertwined. While the overly sensitive Ly-alpha transition completely saturates in a neutral IGM, strong low-ionization metal lines like the MgII 2796,2804 doublet will give rise to a detectable `metal-line forest' if the metals produced during reionization (Z ~ 10^{-3}Z_sol) permeate the neutral IGM. We simulate the MgII forest for the first time by combining a large hydrodynamical simulation with a semi-numerical reionization topology, assuming a simple enrichment model where the IGM is uniformly suffused with metals. In contrast to the traditional approach of identifying discrete absorbers, we treat the absorption as a continuous random field and measure its two-point correlation function, leveraging techniques from precision cosmology. We show that a realistic mock dataset of 10 JWST spectra can simultaneously determine the Mg abundance, [Mg/H], with a 1sigma precision of 0.02 dex and measure the global neutral fraction <x_HI> to 5% for a Universe with <x_HI> = 0.74 and [Mg/H] = -3.7. Alternatively, if the IGM is pristine, a null-detection of the MgII forest would set a stringent upper limit on the IGM metallicity of [Mg/H] < -4.4 at 95% credibility, assuming <x_HI> > 0.5 from another probe. Concentrations of metals in the circumgalactic environs of galaxies can significantly contaminate the IGM signal, but we demonstrate how these discrete absorbers can be easily identified and masked such that their impact on the correlation function is negligible. The MgII forest thus has tremendous potential to precisely constrain the reionization and enrichment history of the Universe. △ Less

Submitted 30 July, 2020; originally announced July 2020.

Comments: 22 pages, 13 figures, submitted to MNRAS

Anomaly in the Opacity of the Post-Reionization Intergalactic Medium in the Ly$α$ and Ly$β$ Forest

Authors: Anna-Christina Eilers, Joseph F. Hennawi, Frederick B. Davies, Jose Oñorbe

Abstract: We measure the intergalactic medium (IGM) opacity in the Ly$α$ as well as in the Ly$β$ forest along $19$ quasar sightlines between $5.5\lesssim z_{\rm abs}\lesssim 6.1$, probing the end stages of the reionization epoch. Owing to its lower oscillator strength the Ly$β$ transition is sensitive to different gas temperatures and densities than Ly$α$, providing additional constraints on the ionization… ▽ More We measure the intergalactic medium (IGM) opacity in the Ly$α$ as well as in the Ly$β$ forest along $19$ quasar sightlines between $5.5\lesssim z_{\rm abs}\lesssim 6.1$, probing the end stages of the reionization epoch. Owing to its lower oscillator strength the Ly$β$ transition is sensitive to different gas temperatures and densities than Ly$α$, providing additional constraints on the ionization and thermal state of the IGM. A comparison of our measurements to different inhomogeneous reionization models, derived from post-processing the Nyx cosmological hydrodynamical simulation to include spatial fluctuations in the ultraviolet background (UVB) or the gas temperature field, as well as to a uniform reionization model with varying thermal states of the IGM, leads to two primary conclusions: First, we find that including the effects of spectral noise is key for a proper data to model comparison. Noise effectively reduces the sensitivity to high opacity regions, and thus even stronger spatial inhomogeneities are required to match the observed scatter in the observations than previously inferred. Second, we find that models which come close to reproducing the distribution of Ly$α$ effective optical depths nevertheless underpredict the Ly$β$ opacity at the same spatial locations. The origin of this disagreement is not entirely clear but models with an inversion in the temperature-density relation of the IGM just after reionization is completed match our measurements best, although they still do not fully capture the observations at $z\gtrsim 5.8$. △ Less

Submitted 18 June, 2019; v1 submitted 13 June, 2019; originally announced June 2019.

Comments: 14 pages+appendix. Accepted for publication in ApJ

A Novel Statistical Method for Measuring the Temperature-Density Relation in the IGM Using the $b$-$N_{\text{HI}}$ Distribution of Absorbers in the Ly$α$ Forest

Authors: Hector Hiss, Michael Walther, José Oñorbe, Joseph F. Hennawi

Abstract: We present a new method for determining the thermal state of the intergalactic medium based on Voigt profile decomposition of the Ly$α$ forest. The distribution of Doppler parameter and column density ($b$-$N_{\text{HI}}$ distribution) is sensitive to the temperature density relation $T=T_0 (ρ/ρ_0)^{γ-1}$, and previous work has inferred $T_0$ and $γ$ by fitting its low-$b$ cutoff. This approach di… ▽ More We present a new method for determining the thermal state of the intergalactic medium based on Voigt profile decomposition of the Ly$α$ forest. The distribution of Doppler parameter and column density ($b$-$N_{\text{HI}}$ distribution) is sensitive to the temperature density relation $T=T_0 (ρ/ρ_0)^{γ-1}$, and previous work has inferred $T_0$ and $γ$ by fitting its low-$b$ cutoff. This approach discards the majority of available data, and is susceptible to systematics related to cutoff determination. We present a method that exploits all information encoded in the $b$-$N_{\text{HI}}$ distribution by modeling its entire shape. We apply kernel density estimation to discrete absorption lines to generate model probability density functions, then use principal component decomposition to create an emulator which can be evaluated anywhere in thermal parameter space. We introduce a Bayesian likelihood based on these models enabling parameter inference via Markov chain Monte Carlo. The method's robustness is tested by applying it to a large grid of thermal history simulations. By conducting 160 mock measurements we establish that our approach delivers unbiased estimates and valid uncertainties for a 2D $(T_0, γ)$ measurement. Furthermore, we conduct a pilot study applying this methodology to real observational data at $z=2$. Using 200 absorbers, equivalent in pathlength to a single Ly$α$ forest spectrum, we measure $\log T_0 =4.092^{+0.050}_{-0.055}$ and $γ=1.49^{+0.073}_{-0.074}$ in excellent agreement with cutoff fitting determinations using the same data. Our method is far more sensitive than cutoff fitting, enabling measurements of $\log T_0$ and $γ$ with precision on $\log T_0$ ($γ$) nearly two (three) times higher for current dataset sizes. △ Less

Submitted 28 March, 2019; originally announced March 2019.

Comments: Accepted for publication in ApJ, 23 pages, 15 figures

Inhomogeneous Reionization Models in Cosmological Hydrodynamical Simulations

Authors: Jose Oñorbe, F. B. Davies, Z. Lukić, J. F. Hennawi, D. Sorini

Abstract: In this work we present a new hybrid method to simulate the thermal effects of the reionization in cosmological hydrodynamical simulations. The method improves upon the standard approach used in simulations of the intergalactic medium (IGM) and galaxy formation without a significant increase of the computational cost allowing for efficient exploration of the parameter space. The method uses a smal… ▽ More In this work we present a new hybrid method to simulate the thermal effects of the reionization in cosmological hydrodynamical simulations. The method improves upon the standard approach used in simulations of the intergalactic medium (IGM) and galaxy formation without a significant increase of the computational cost allowing for efficient exploration of the parameter space. The method uses a small set of phenomenological input parameters and combines a semi-numerical reionization model to solve for the topology of reionization and an approximate model of how reionization heats the IGM, with the massively parallel \texttt{Nyx} hydrodynamics code, specifically designed to solve for the structure of diffuse IGM gas. We have produced several large-scale high resolution cosmological hydrodynamical simulations ($2048^3$, $L_{\rm box} = 40$ Mpc/h) with different instantaneous and inhomogeneous HI reionization models that use this new methodology. We study the IGM thermal properties of these models and find that large scale temperature fluctuations extend well beyond the end of reionization. Analyzing the 1D flux power spectrum of these models, we find up to $\sim 50\%$ differences in the large scale properties (low modes, $k\lesssim0.01$ s/km) of the post-reionization power spectrum due to the thermal fluctuations. We show that these differences could allow one to distinguish between different reionization scenarios already with existing Ly$α$ forest measurements. Finally, we explore the differences in the small-scale cutoff of the power spectrum and we find that, for the same heat input, models show very good agreement provided that the reionization redshift of the instantaneous reionization model happens at the midpoint of the inhomogeneous model. △ Less

Submitted 7 June, 2019; v1 submitted 27 October, 2018; originally announced October 2018.

Comments: 24 pages, 16 figures. Accepted by MNRAS. Minor changes to match published version

Mapping quasar light echoes in 3D with Lyα forest tomography

Authors: Tobias M. Schmidt, Joseph F. Hennawi, Khee-Gan Lee, Zarija Lukic, Jose Onorbe, Martin White

Abstract: The intense radiation emitted by luminous quasars dramatically alters the ionization state of their surrounding IGM. This so-called proximity effect extends out to tens of Mpc, and manifests as large coherent regions of enhanced Lyman-$α$ (Ly$α$) forest transmission in absorption spectra of background sightlines. Here we present a novel method based on Ly$α$ forest tomography, which is capable of… ▽ More The intense radiation emitted by luminous quasars dramatically alters the ionization state of their surrounding IGM. This so-called proximity effect extends out to tens of Mpc, and manifests as large coherent regions of enhanced Lyman-$α$ (Ly$α$) forest transmission in absorption spectra of background sightlines. Here we present a novel method based on Ly$α$ forest tomography, which is capable of mapping these quasar `light echoes' in three dimensions. Using a dense grid (10-100) of faint ($m_r\approx24.7\,\mathrm{mag}$) background galaxies as absorption probes, one can measure the ionization state of the IGM in the vicinity of a foreground quasar, yielding detailed information about the quasar's radiative history and emission geometry. An end-to-end analysis - combining cosmological hydrodynamical simulations post-processed with a quasar emission model, realistic estimates of galaxy number densities, and instrument + telescope throughput - is conducted to explore the feasibility of detecting quasar light echoes. We present a new fully Bayesian statistical method that allows one to reconstruct quasar light echoes from thousands of individual low S/N transmission measurements. Armed with this machinery, we undertake an exhaustive parameter study and show that light echoes can be convincingly detected for luminous ($M_{1450} < -27.5\,\mathrm{mag}$ corresponding to $m_{1450} < 18.4\,\mathrm{mag}$ at $z\simeq 3.6$) quasars at redshifts $3<z_\mathrm{QSO}<5$, and that a relative precision better than $20\,\%$ on the quasar age can be achieved for individual objects, for the expected range of ages between 1 Myr and 100 Myr. The observational requirements are relatively modest - moderate resolution ($R\gtrsim750$) multi object spectroscopy at low $\rm{}S/N > 5$ is sufficient, requiring three hour integrations using existing instruments on 8m class telescopes. △ Less

Submitted 11 October, 2018; originally announced October 2018.

Comments: 22 pages, 21 figures

QSO MUSEUM I: A sample of 61 extended Ly$α$-emission nebulae surrounding $z\sim3$ quasars

Authors: F. Arrigoni Battaia, J. F. Hennawi, J. X. Prochaska, Jose Oñorbe, E. P. Farina, S. Cantalupo, E. Lusso

Abstract: Motivated by the recent discovery of rare Enormous Lyman-Alpha Nebulae (ELAN) around z~2 quasars, we have initiated a long-term observational campaign with the MUSE instrument to directly uncover the astrophysics of the gas around quasars. We present here the first 61 targets of our effort under the acronym QSO MUSEUM (Quasar Snapshot Observations with MUse: Search for Extended Ultraviolet eMissio… ▽ More Motivated by the recent discovery of rare Enormous Lyman-Alpha Nebulae (ELAN) around z~2 quasars, we have initiated a long-term observational campaign with the MUSE instrument to directly uncover the astrophysics of the gas around quasars. We present here the first 61 targets of our effort under the acronym QSO MUSEUM (Quasar Snapshot Observations with MUse: Search for Extended Ultraviolet eMission). These quasars have a median redshfit of z=3.17, absolute $i$ magnitude in the range $-29.67\leq M_i(z=2)\leq-27.03$, and different levels of radio-loudness. This sample unveils diverse specimens of Ly$α$ nebulosities extending for tens of kiloparsecs around these quasars (on average out to a maximum projected distance of 80 kpc) above a surface brightness SB$>8.8\times10^{-19}$ erg s$^{-1}$ cm$^{-2}$ arcsec$^{-2}$ ($2σ$). The bulk of the extended Ly$α$ emission is within R<50 kpc, and is characterized by relatively quiescent kinematics, with average velocity dispersions of $\langle σ_{\rm Lyα}\rangle < 400$~km~s$^{-1}$. Therefore, the motions within all these Ly$α$ nebulosities have amplitudes consistent with gravitational motions expected in dark matter halos hosting quasars at these redshifts, possibly reflecting the complexity in propagating a fast wind on large scales. Our current data suggest a combination of photoionization and resonant scattering as powering mechanisms of the Ly$α$ emission. We discover the first $z\sim3$ ELAN, which confirms a very low probability ($\sim1\%$) of occurrence of such extreme systems at these cosmic epochs. Finally, we discuss the redshift evolution currently seen in extended Ly$α$ emission around radio-quiet quasars from $z\sim3$ to $z\sim2$, concluding that it is possibly linked to a decrease of cool gas mass within the quasars' CGM from $z\sim3$ to $z\sim2$, and thus to the balance of cool vs hot media. △ Less

Submitted 31 August, 2018; originally announced August 2018.

Comments: 48 pages, 3 page-long tables, 28 figures, 3 appendices. Submitted to MNRAS. Comments are welcome

The Power Spectrum of the Lyman-$α$ Forest at z < 0.5

Authors: Vikram Khaire, Michael Walther, Joseph F. Hennawi, Jose Oñorbe, Zarija Lukić, J. Xavier Prochaska, Todd M. Tripp, Joseph N. Burchett, Christian Rodriguez

Abstract: We present new measurements of the flux power-spectrum P(k) of the $z<0.5$ HI Lyman-$α$ forest spanning scales k ~ 0.001-0.1 s/km. These results were derived from 65 far ultraviolet quasar spectra (resolution R~18000) observed with the Cosmic Origin Spectrograph (COS) on board the Hubble Space Telescope. The analysis required careful masking of all contaminating, coincident absorption from HI and… ▽ More We present new measurements of the flux power-spectrum P(k) of the $z<0.5$ HI Lyman-$α$ forest spanning scales k ~ 0.001-0.1 s/km. These results were derived from 65 far ultraviolet quasar spectra (resolution R~18000) observed with the Cosmic Origin Spectrograph (COS) on board the Hubble Space Telescope. The analysis required careful masking of all contaminating, coincident absorption from HI and metal-line transitions of the Galactic interstellar medium and intervening absorbers as well as proper treatment of the complex COS line-spread function. From the P(k) measurements, we estimate the HI photoionization rate ($Γ_{\rm HI}$) in the z<0.5 intergalactic medium. Our results confirm most of the previous $Γ_{\rm HI}$ estimates. We conclude that previous concerns of a photon underproduction crisis are now resolved by demonstrating that the measured $Γ_{\rm HI}$ can be accounted for by ultraviolet emission from quasars alone. In a companion paper, we will present constraints on the thermal state of the $z<0.5$ intergalactic medium from the P(k) measurements presented here. △ Less

Submitted 9 April, 2019; v1 submitted 16 August, 2018; originally announced August 2018.

Comments: 15 pages, 8 figures, Accepted for publication in MNRAS

New Constraints on IGM Thermal Evolution from the Lyα Forest Power Spectrum

Authors: Michael Walther, Jose Oñorbe, Joseph F. Hennawi, Zarija Lukić

Abstract: We determine the thermal evolution of the intergalactic medium (IGM) over $3\, \mathrm{Gyr}$ of cosmic time $1.8<z<5.4$ by comparing measurements of the Lyα forest power spectrum to a suite of $\sim70$ hydrodynamical simulations. We conduct Bayesian inference of IGM thermal parameters using an end-to-end forward modeling framework whereby mock spectra generated from our simulation grid are used to… ▽ More We determine the thermal evolution of the intergalactic medium (IGM) over $3\, \mathrm{Gyr}$ of cosmic time $1.8<z<5.4$ by comparing measurements of the Lyα forest power spectrum to a suite of $\sim70$ hydrodynamical simulations. We conduct Bayesian inference of IGM thermal parameters using an end-to-end forward modeling framework whereby mock spectra generated from our simulation grid are used to build a custom emulator which interpolates the power spectrum between thermal grid points. The temperature at mean density $T_0$ rises steadily from $T_0\sim 6000\, \mathrm{K}$ at $z=5.4$, peaks at $14000\, \mathrm{K}$ for $z\sim 3.4$, and decreases at lower redshift reaching $T_0\sim 7000\, \mathrm{K}$ by $z\sim1.8$. This evolution provides conclusive evidence for photoionization heating resulting from the reionization of He II, as well as the subsequent cooling of the IGM due to the expansion of the Universe after all reionization events are complete. Our results are broadly consistent with previous measurements of thermal evolution based on a variety of approaches, but the sensitivity of the power spectrum, the combination of high precision BOSS measurements of large-scale modes ($k\lesssim 0.02\, \mathrm{s/km}$) with our recent determination of the small-scale power, our large grid of models, and our careful statistical analysis allow us to break the well known degeneracy between the temperature at mean density $T_0$ and the slope of the temperature density relation $γ$ that has plagued previous analyses. At the highest redshifts $z\geq5$ we infer lower temperatures than expected from the standard picture of IGM thermal evolution leaving little room for additional smoothing of the Lyα forest by free streaming of warm dark matter. △ Less

Submitted 20 December, 2018; v1 submitted 13 August, 2018; originally announced August 2018.

Comments: 32 pages, 18 figures, accepted for publication in ApJ

Modeling the HeII Transverse Proximity Effect: Constraints on Quasar Lifetime and Obscuration

Authors: Tobias M. Schmidt, Joseph F. Hennawi, Gábor Worseck, Frederick B. Davies, Zarija Lukić, Jose Oñorbe

Abstract: The HeII transverse proximity effect - enhanced HeII Lyα transmission in a background sightline caused by the ionizing radiation of a foreground quasar - offers a unique opportunity to probe the emission properties of quasars, in particular the emission geometry (obscuration, beaming) and the quasar lifetime. Building on the foreground quasar survey published in Schmidt+2017, we present a detailed… ▽ More The HeII transverse proximity effect - enhanced HeII Lyα transmission in a background sightline caused by the ionizing radiation of a foreground quasar - offers a unique opportunity to probe the emission properties of quasars, in particular the emission geometry (obscuration, beaming) and the quasar lifetime. Building on the foreground quasar survey published in Schmidt+2017, we present a detailed model of the HeII transverse proximity effect, specifically designed to include light travel time effects, finite quasar ages, and quasar obscuration. We post-process outputs from a cosmological hydrodynamical simulation with a fluctuating HeII UV background model, plus the added effect of the radiation from a single bright foreground quasar. We vary the age $t_\mathrm{age}$ and obscured sky fractions $Ω_\mathrm{obsc}$ of the foreground quasar, and explore the resulting effect on the HeII transverse proximity effect signal. Fluctuations in IGM density and the UV background, as well as the unknown orientation of the foreground quasar, result in a large variance of the HeII Lyα transmission along the background sightline. We develop a fully Bayesian statistical formalism to compare far UV HeII Lyα transmission spectra of the background quasars to our models, and extract joint constraints on $t_\mathrm{age}$ and $Ω_\mathrm{obsc}$ for the six Schmidt+2017 foreground quasars with the highest implied HeII photoionization rates. Our analysis suggests a bimodal distribution of quasar emission properties, whereby one foreground quasar, associated with a strong HeII transmission spike, is relatively old $(22\,\mathrm{Myr})$ and unobscured $Ω_\mathrm{obsc}<35\%$, whereas three others are either younger than $(10\,\mathrm{Myr})$ or highly obscured $(Ω_\mathrm{obsc}>70\%)$. △ Less

Submitted 27 October, 2017; v1 submitted 12 October, 2017; originally announced October 2017.

Comments: 19 pages, 6 figures, submitted to ApJ

A New Measurement of the Temperature Density Relation of the IGM From Voigt Profile Fitting

Authors: Hector Hiss, Michael Walther, Joseph F. Hennawi, José Oñorbe, John M. O'Meara, Alberto Rorai

Abstract: We decompose the Lyman-α (Lyα) forest of an extensive sample of 74 high signal-to-noise ratio and high-resolution quasar spectra into a collection of Voigt profiles. Absorbers located near caustics in the peculiar velocity field have the smallest Doppler parameters, resulting in a low-$b$ cutoff in the $b$-$N_{\text{HI}}$ set by the thermal state of intergalactic medium (IGM). We fit this cutoff a… ▽ More We decompose the Lyman-α (Lyα) forest of an extensive sample of 74 high signal-to-noise ratio and high-resolution quasar spectra into a collection of Voigt profiles. Absorbers located near caustics in the peculiar velocity field have the smallest Doppler parameters, resulting in a low-$b$ cutoff in the $b$-$N_{\text{HI}}$ set by the thermal state of intergalactic medium (IGM). We fit this cutoff as a function of redshift over the range $2.0\leq z \leq 3.4$, which allows us to measure the evolution of the IGM temperature-density ($T= T_0 (ρ/ ρ_0)^{γ-1}$) relation parameters $T_0$ and $γ$. We calibrate our measurements against Ly$α$ forest simulations, using 21 different thermal models of the IGM at each redshift, also allowing for different values of the IGM pressure smoothing scale. We adopt a forward-modeling approach and self-consistently apply the same algorithms to both data and simulations, propagating both statistical and modeling uncertainties via Monte Carlo. The redshift evolution of $T_0$ shows a suggestive peak at $z=2.8$, while our evolution of $γ$ is consistent with $γ\simeq 1.4$ and disfavors inverted temperature-density relations. Our measured evolution of $T_0$ and $γ$ are generally in good agreement with previous determinations in the literature. Both the peak in the evolution of $T_0$ at $z = 2.8$, as well as the high temperatures $T_0\simeq 15000-20000\,$K that we observe at $2.4 < z < 3.4$, strongly suggest that a significant episode of heating occurred after the end of HI reionization, which was most likely the cosmic reionization of HeII. △ Less

Submitted 16 August, 2018; v1 submitted 2 October, 2017; originally announced October 2017.

Comments: Accepted for publication in ApJ, 23 pages, 26 figures, machine readable tables available online

A New Precision Measurement of the Small-Scale Line-of-Sight Power Spectrum of the Ly$α$ Forest

Authors: Michael Walther, Joseph F. Hennawi, Hector Hiss, Jose Oñorbe, Khee-Gan Lee, Alberto Rorai, John O'Meara

Abstract: We present a new measurement of the Ly$α$ forest power spectrum at $1.8 < z < 3.4$ using 74 Keck/HIRES and VLT/UVES high-resolution, high-S/N quasar spectra. We developed a custom pipeline to measure the power spectrum and its uncertainty, which fully accounts for finite resolution and noise, and corrects for the bias induced by masking missing data, DLAs, and metal absorption lines. Our measureme… ▽ More We present a new measurement of the Ly$α$ forest power spectrum at $1.8 < z < 3.4$ using 74 Keck/HIRES and VLT/UVES high-resolution, high-S/N quasar spectra. We developed a custom pipeline to measure the power spectrum and its uncertainty, which fully accounts for finite resolution and noise, and corrects for the bias induced by masking missing data, DLAs, and metal absorption lines. Our measurement results in unprecedented precision on the small-scale modes $k > 0.02\,\mathrm{s\,km^{-1}}$, unaccessible to previous SDSS/BOSS analyses. It is well known that these high-$k$ modes are highly sensitive to the thermal state of the intergalactic medium, however contamination by narrow metal lines is a significant concern. We quantify the effect of metals on the small-scale power, and find a modest effect on modes with $k < 0.1\,\mathrm{s\,km^{-1}} $. As a result, by masking metals and restricting to $k < 0.1\,\mathrm{s\,km^{-1}}$ their impact is completely mitigated. We present an end-to-end Bayesian forward modeling framework whereby mock spectra with the same noise, resolution, and masking as our data are generated from Ly$α$ forest simulations. These mocks are used to build a custom emulator, enabling us to interpolate between a sparse grid of models and perform MCMC fits. Our results agree well with BOSS on scales $k < 0.02\,\mathrm{s\,km^{-1}}$ where the measurements overlap. The combination of BOSS' percent level low-$k$ precision with our $5-15\%$ high-$k$ measurements, results in a powerful new dataset for precisely constraining the thermal history of the intergalactic medium, cosmological parameters, and the nature of dark matter. The power spectra and their covariance matrices are provided as electronic tables. △ Less

Submitted 30 November, 2017; v1 submitted 21 September, 2017; originally announced September 2017.

Comments: 24 pages, 12 figures, accepted for publication in ApJ, machine readable tables will be made available after publication in the journal

A Fundamental Test for Galaxy Formation Models: Matching the Lyman-$α$ Absorption Profiles of Galactic Halos over Three Decades in Distance

Authors: Daniele Sorini, José Oñorbe, Joseph F. Hennawi, Zarija Lukić

Abstract: Galaxy formation depends critically on the physical state of gas in the circumgalactic medium (CGM) and its interface with the intergalactic medium (IGM), determined by the complex interplay between inflows from the IGM and outflows from supernovae or AGN feedback. The average Lyman-alpha (Ly-a) absorption profile around galactic halos represents a powerful tool to probe their gaseous environments… ▽ More Galaxy formation depends critically on the physical state of gas in the circumgalactic medium (CGM) and its interface with the intergalactic medium (IGM), determined by the complex interplay between inflows from the IGM and outflows from supernovae or AGN feedback. The average Lyman-alpha (Ly-a) absorption profile around galactic halos represents a powerful tool to probe their gaseous environments. We compare predictions from Illustris and Nyx hydrodynamical simulations with the observed absorption around foreground quasars, damped Ly-a systems, and Lyman-break galaxies. We show how large-scale BOSS and small-scale quasar pair measurements can be combined to precisely constrain the absorption profile over three decades in transverse distance 20kpc$\lesssim b\lesssim$20Mpc. Far from galaxies $\gtrsim2$Mpc, the simulations converge to the same profile and provide a reasonable match to the observations. This asymptotic agreement arises because the $Λ$CDM model successfully describes the ambient IGM, and represents a critical advantage of studying the mean absorption profile. However, significant differences between the simulations, and between simulations and observations are present on scales 20kpc$\lesssim b\lesssim$2Mpc, illustrating the challenges of accurately modeling and resolving galaxy formation physics. It is noteworthy that these differences are observed as far out as $\sim2$Mpc, indicating that the `sphere-of-influence' of galaxies could extend to approximately $\sim7$ times the halo virial radius ($\sim100$kpc). Current observations are very precise on these scales and can thus strongly discriminate between different galaxy formation models. We demonstrate that the Ly-a absorption profile is primarily sensitive to the underlying temperature-density relationship of diffuse gas around galaxies, and argue that it thus provides a fundamental test of galaxy formation models. △ Less

Submitted 15 November, 2021; v1 submitted 12 September, 2017; originally announced September 2017.

Comments: Published in ApJ. This is the Accepted Manuscript version of an article accepted for publication in The Astrophysical Journal. IOP Publishing Ltd is not responsible for any errors or omissions in this version of the manuscript or any version derived from it. The Version of Record is available online at https://iopscience.iop.org/article/10.3847/1538-4357/aabb52

Measurement of the small-scale structure of the intergalactic medium using close quasar pairs

Authors: Alberto Rorai, Joseph F. Hennawi, Jose Oñorbe, Martin White, J. Xavier Prochaska, Girish Kulkarni, Michael Walther, Zarija Lukić, Khee-Gan Lee

Abstract: The distribution of diffuse gas in the intergalactic medium (IGM) imprints a series of hydrogen absorption lines on the spectra of distant background quasars known as the Lyman-$α$ forest. Cosmological hydrodynamical simulations predict that IGM density fluctuations are suppressed below a characteristic scale where thermal pressure balances gravity. We measured this pressure-smoothing scale by qua… ▽ More The distribution of diffuse gas in the intergalactic medium (IGM) imprints a series of hydrogen absorption lines on the spectra of distant background quasars known as the Lyman-$α$ forest. Cosmological hydrodynamical simulations predict that IGM density fluctuations are suppressed below a characteristic scale where thermal pressure balances gravity. We measured this pressure-smoothing scale by quantifying absorption correlations in a sample of close quasar pairs. We compared our measurements to hydrodynamical simulations, where pressure smoothing is determined by the integrated thermal history of the IGM. Our findings are consistent with standard models for photoionization heating by the ultraviolet radiation backgrounds that reionized the universe. △ Less

Submitted 26 April, 2017; originally announced April 2017.

Comments: Accepted for publication on Science

Constraining Reionization with the $z \sim 5-6$ Lyman-$α$ Forest Power Spectrum: the Outlook after Planck

Authors: J. Oñorbe, J. F. Hennawi, Z. Lukić, M. Walther

Abstract: The latest measurements of CMB electron scattering optical depth reported by Planck significantly reduces the allowed space of HI reionization models, pointing toward a later ending and/or less extended phase transition than previously believed. Reionization impulsively heats the intergalactic medium (IGM) to $\sim10^4$ K, and owing to long cooling and dynamical times in the diffuse gas, comparabl… ▽ More The latest measurements of CMB electron scattering optical depth reported by Planck significantly reduces the allowed space of HI reionization models, pointing toward a later ending and/or less extended phase transition than previously believed. Reionization impulsively heats the intergalactic medium (IGM) to $\sim10^4$ K, and owing to long cooling and dynamical times in the diffuse gas, comparable to the Hubble time, memory of reionization heating is retained. Therefore, a late ending reionization has significant implications for the structure of the $z\sim5-6$ Lyman-$α$ (ly$α$) forest. Using state-of-the-art hydrodynamical simulations that allow us to vary the timing of reionization and its associated heat injection, we argue that extant thermal signatures from reionization can be detected via the Ly$α$ forest power spectrum at $5< z<6$. This arises because the small-scale cutoff in the power depends not only the the IGM temperature at these epochs, but is also particularly sensitive to the pressure smoothing scale set by the IGM full thermal history. Comparing our different reionization models with existing measurements of the Ly$α$ forest flux power spectrum at $z=5.0-5.4$, we find that models satisfying Planck's $τ_e$, constraint favor a moderate amount of heat injection consistent with galaxies driving reionization, but disfavoring quasar-driven scenarios. We explore the impact of different reionization histories and heating models on the shape of the power spectrum, and find that they can produce similar effects, but argue that this degeneracy can be broken with high enough quality data. We study the feasibility of measuring the flux power spectrum at $z\simeq 6$ using mock quasar spectra and conclude that a sample of $\sim10$ high-resolution spectra with an attainable signal-to-noise ratio will allow distinguishing between different reionization scenarios. △ Less

Submitted 27 April, 2019; v1 submitted 24 March, 2017; originally announced March 2017.

Comments: 14 pages, 7 Figures. Matches version accepted by ApJ. UVB models used in this work can be downloaded from http://thermal.joseonorbe.com/uvbtables.html

FIRE in the Field: Simulating the Threshold of Galaxy Formation

Authors: Alex Fitts, Michael Boylan-Kolchin, Oliver D. Elbert, James S. Bullock, Philip F. Hopkins, Jose Onorbe, Andrew R. Wetzel, Coral Wheeler, Claude-Andre Faucher-Giguere, Dusan Keres, Evan D. Skillman, Daniel R. Weisz

Abstract: We present a suite of 15 cosmological zoom-in simulations of isolated dark matter halos, all with masses of $M_{\rm halo} \approx 10^{10}\,{\rm M}_\odot$ at $z=0$, in order to understand the relationship between halo assembly, galaxy formation, and feedback's effects on the central density structure in dwarf galaxies. These simulations are part of the Feedback in Realistic Environments (FIRE) proj… ▽ More We present a suite of 15 cosmological zoom-in simulations of isolated dark matter halos, all with masses of $M_{\rm halo} \approx 10^{10}\,{\rm M}_\odot$ at $z=0$, in order to understand the relationship between halo assembly, galaxy formation, and feedback's effects on the central density structure in dwarf galaxies. These simulations are part of the Feedback in Realistic Environments (FIRE) project and are performed at extremely high resolution. The resultant galaxies have stellar masses that are consistent with rough abundance matching estimates, coinciding with the faintest galaxies that can be seen beyond the virial radius of the Milky Way ($M_\star/{\rm M}_\odot\approx 10^5-10^7$). This non-negligible spread in stellar mass at $z=0$ in halos within a narrow range of virial masses is strongly correlated with central halo density or maximum circular velocity $V_{\rm max}$. Much of this dependence of $M_\star$ on a second parameter (beyond $M_{\rm halo}$) is a direct consequence of the $M_{\rm halo}\sim10^{10}\,{\rm M}_\odot$ mass scale coinciding with the threshold for strong reionization suppression: the densest, earliest-forming halos remain above the UV-suppression scale throughout their histories while late-forming systems fall below the UV-suppression scale over longer periods and form fewer stars as a result. In fact, the latest-forming, lowest-concentration halo in our suite fails to form any stars. Halos that form galaxies with $M_\star\gtrsim2\times10^{6}\,{\rm M}_\odot$ have reduced central densities relative to dark-matter-only simulations, and the radial extent of the density modifications is well-approximated by the galaxy half-mass radius $r_{1/2}$. This apparent stellar mass threshold of $M_\star \approx 2\times 10^{6} \approx 2\times 10^{-4} \,M_{\rm halo}$ is broadly consistent with previous work and provides a testable prediction of FIRE feedback models in LCDM. △ Less

Submitted 30 September, 2017; v1 submitted 7 November, 2016; originally announced November 2016.

Comments: 16 pages, 16 figures; matches version published in MNRAS

Journal ref: Mon.Not.Roy.Astron.Soc.471:3547-3562,2017

Self-Consistent Modeling of Reionization in Cosmological Hydrodynamical Simulations

Authors: Jose Oñorbe, Joseph F Hennawi, Zarija Lukić

Abstract: The ultraviolet background (UVB) emitted by quasars and galaxies governs the ionization and thermal state of the intergalactic medium (IGM), regulates the formation of high-redshift galaxies, and is thus a key quantity for modeling cosmic reionization. The vast majority of cosmological hydrodynamical simulations implement the UVB via a set of spatially uniform photoionization and photoheating rate… ▽ More The ultraviolet background (UVB) emitted by quasars and galaxies governs the ionization and thermal state of the intergalactic medium (IGM), regulates the formation of high-redshift galaxies, and is thus a key quantity for modeling cosmic reionization. The vast majority of cosmological hydrodynamical simulations implement the UVB via a set of spatially uniform photoionization and photoheating rates derived from UVB synthesis models. We show that simulations using canonical UVB rates reionize and, perhaps more importantly, spuriously heat the IGM, much earlier z ~ 15 than they should. This problem arises because at z > 6, where observational constraints are nonexistent, the UVB amplitude is far too high. We introduce a new methodology to remedy this issue, and we generate self-consistent photoionization and photoheating rates to model any chosen reionization history. Following this approach, we run a suite of hydrodynamical simulations of different reionization scenarios and explore the impact of the timing of reionization and its concomitant heat injection on the the thermal state of the IGM. We present a comprehensive study of the pressure smoothing scale of IGM gas, illustrating its dependence on the details of both hydrogen and helium reionization, and argue that it plays a fundamental role in interpreting Lyman-alpha forest statistics and the thermal evolution of the IGM. The premature IGM heating we have uncovered implies that previous work has likely dramatically overestimated the impact of photoionization feedback on galaxy formation, which sets the minimum halo mass able to form stars at high redshifts. We make our new UVB photoionization and photoheating rates publicly available for use in future simulations. △ Less

Submitted 16 March, 2017; v1 submitted 14 July, 2016; originally announced July 2016.

Comments: 28 pages, 15 figures, 5 tables, matches version accepted in ApJ

High Angular Momentum Halo Gas: a Feedback and Code-Independent Prediction of LCDM

Authors: Kyle Stewart, Ariyeh Maller, Jose Oñorbe, James Bullock, M. Ryan Joung, Julien Devriendt, Daniel Ceverino, Dušan Kereš, Phil Hopkins, Claude-André Faucher-Giguère

Abstract: We investigate angular momentum acquisition in Milky Way-sized galaxies by comparing five high resolution zoom-in simulations, each implementing identical cosmological initial conditions, but utilizing different hydrodynamic codes: Enzo, Art, Ramses, Arepo, and Gizmo-PSPH. Each code implements a distinct set of feedback and star formation prescriptions. We find that while many galaxy and halo prop… ▽ More We investigate angular momentum acquisition in Milky Way-sized galaxies by comparing five high resolution zoom-in simulations, each implementing identical cosmological initial conditions, but utilizing different hydrodynamic codes: Enzo, Art, Ramses, Arepo, and Gizmo-PSPH. Each code implements a distinct set of feedback and star formation prescriptions. We find that while many galaxy and halo properties vary between the different codes (and feedback prescriptions), there is qualitative agreement on the process of angular momentum acquisition in the galaxy's halo. In all simulations, cold filamentary gas accretion to the halo results in ~4 times more specific angular momentum in cold halo gas ($λ_{cold} \gtrsim 0.1$) than in the dark matter halo. At z>1, this inflow takes the form of inspiraling cold streams that are co-directional in the halo of the galaxy and are fueled, aligned, and kinematically connected to filamentary gas infall along the cosmic web. Due to the qualitative agreement among disparate simulations, we conclude that the buildup of high angular momentum halo gas and the presence of these inspiraling cold streams are robust predictions of Lambda Cold Dark Matter galaxy formation, though the detailed morphology of these streams is significantly less certain. A growing body of observational evidence suggests that this process is borne out in the real universe. △ Less

Submitted 23 August, 2017; v1 submitted 27 June, 2016; originally announced June 2016.

Comments: 15 pages, 7 figures, 1 table. Revised to match published version, including significantly expanded and discussion. Primary qualitative results unchanged

Journal ref: ApJ, 843, 47S, 2017

Modeling the Lyman-alpha Forest in Collisionless Simulations

Authors: Daniele Sorini, José Oñorbe, Zarija Lukić, Joseph F. Hennawi

Abstract: Cosmological hydrodynamic simulations can accurately predict the properties of the intergalactic medium (IGM), but only under the condition of retaining high spatial resolution necessary to resolve density fluctuations in the IGM. This resolution constraint prohibits simulating large volumes, such as those probed by BOSS and future surveys, like DESI and 4MOST. To overcome this limitation, we pres… ▽ More Cosmological hydrodynamic simulations can accurately predict the properties of the intergalactic medium (IGM), but only under the condition of retaining high spatial resolution necessary to resolve density fluctuations in the IGM. This resolution constraint prohibits simulating large volumes, such as those probed by BOSS and future surveys, like DESI and 4MOST. To overcome this limitation, we present Iteratively Matched Statistics (IMS), a novel method to accurately model the Lyman-alpha forest with collisionless N-body simulations, where the relevant density fluctuations are unresolved. We use a small-box, high-resolution hydrodynamic simulation to obtain the probability distribution function (PDF) and the power spectrum of the real-space Lyman-alpha forest flux. These two statistics are iteratively mapped onto a pseudo-flux field of an N-body simulation, which we construct from the matter density. We demonstrate that our method can perfectly reproduce line-of-sight observables, such as the PDF and power spectrum, and accurately reproduce the 3D flux power spectrum (5-20%). We quantify the performance of the commonly used Gaussian smoothing technique and show that it has significantly lower accuracy (20-80%), especially for N-body simulations with achievable mean inter-particle separations in large-volume simulations. In addition, we show that IMS produces reasonable and smooth spectra, making it a powerful tool for modeling the IGM in large cosmological volumes and for producing realistic "mock" skies for Lyman-alpha forest surveys. △ Less

Submitted 25 February, 2016; originally announced February 2016.

Comments: 25 pages, 15 figures, submitted to ApJ

Properties of Resonantly Produced Sterile Neutrino Dark Matter Subhalos

Authors: Shunsaku Horiuchi, Brandon Bozek, Kevork N. Abazajian, Michael Boylan-Kolchin, James S. Bullock, Shea Garrison-Kimmel, Jose Onorbe

Abstract: The anomalous 3.55 keV X-ray line recently detected towards a number of massive dark matter objects may be interpreted as the radiative decays of 7.1 keV mass sterile neutrino dark matter. Depending on its parameters, the sterile neutrino can range from cold to warm dark matter with small-scale suppression that differs in form from commonly-adopted thermal warm dark matter. Here, we numerically in… ▽ More The anomalous 3.55 keV X-ray line recently detected towards a number of massive dark matter objects may be interpreted as the radiative decays of 7.1 keV mass sterile neutrino dark matter. Depending on its parameters, the sterile neutrino can range from cold to warm dark matter with small-scale suppression that differs in form from commonly-adopted thermal warm dark matter. Here, we numerically investigate the subhalo properties for 7.1 keV sterile neutrino dark matter produced via the resonant Shi-Fuller mechanism. Using accurate matter power spectra, we run cosmological zoom-in simulations of a Milky Way-sized halo and explore the abundance of massive subhalos, their radial distributions, and their internal structure. We also simulate the halo with thermal 2.0 keV warm dark matter for comparison and discuss quantitative differences. We find that the resonantly produced sterile neutrino model for the 3.55 keV line provides a good description of structures in the Local Group, including the number of satellite dwarf galaxies and their radial distribution, and largely mitigates the too-big-to-fail problem. Future searches for satellite galaxies by deep surveys, such as the Dark Energy Survey, Large Synoptic Survey Telescope, and Wide Field Infrared Survey Telescope, will be a strong direct test of warm dark matter scenarios. △ Less

Submitted 14 December, 2015; originally announced December 2015.

Comments: 9 pages, 6 figures, 1 table, accepted for publication in mnras

Report number: INT-PUB-15-072, UCI-HEP-TR-2015-23

Journal ref: MNRAS 456, 4346-4353 (2016)

The no-spin zone: rotation vs dispersion support in observed and simulated dwarf galaxies

Authors: Coral Wheeler, Andrew B. Pace, James S. Bullock, Michael Boylan-Kolchin, Jose Onorbe, Oliver D. Elbert, Alex Fitts, Philip F. Hopkins, Dusan Keres

Abstract: We perform a systematic Bayesian analysis of rotation vs. dispersion support ($v_{\rm rot} / σ$) in $40$ dwarf galaxies throughout the Local Volume (LV) over a stellar mass range $10^{3.5} M_{\rm \odot} < M_{\star} < 10^8 M_{\rm \odot}$. We find that the stars in $\sim 80\%$ of the LV dwarf galaxies studied -- both satellites and isolated systems -- are dispersion-supported. In particular, we show… ▽ More We perform a systematic Bayesian analysis of rotation vs. dispersion support ($v_{\rm rot} / σ$) in $40$ dwarf galaxies throughout the Local Volume (LV) over a stellar mass range $10^{3.5} M_{\rm \odot} < M_{\star} < 10^8 M_{\rm \odot}$. We find that the stars in $\sim 80\%$ of the LV dwarf galaxies studied -- both satellites and isolated systems -- are dispersion-supported. In particular, we show that $6/10$ *isolated* dwarfs in our sample have $v_{\rm rot} / σ< 1.0$. All have $v_{\rm rot} / σ\lesssim 2.0$. These results challenge the traditional view that the stars in gas-rich dwarf irregulars (dIrrs) are distributed in cold, rotationally-supported stellar disks, while gas-poor dwarf spheroidals (dSphs) are kinematically distinct in having dispersion-supported stars. We see no clear trend between $v_{\rm rot} / σ$ and distance to the closest $\rm L_{\star}$ galaxy, nor between $v_{\rm rot} / σ$ and $M_{\star}$ within our mass range. We apply the same Bayesian analysis to four FIRE hydrodynamic zoom-in simulations of isolated dwarf galaxies ($10^9 M_{\odot} < M_{\rm vir} < 10^{10} M_{\rm \odot}$) and show that the simulated *isolated* dIrr galaxies have stellar ellipticities and stellar $v_{\rm rot} / σ$ ratios that are consistent with the observed population of dIrrs *and* dSphs without the need to subject these dwarfs to any external perturbations or tidal forces. We posit that most dwarf galaxies form as puffy, dispersion-dominated systems, rather than cold, angular momentum-supported disks. If this is the case, then transforming a dIrr into a dSph may require little more than removing its gas. △ Less

Submitted 12 October, 2016; v1 submitted 3 November, 2015; originally announced November 2015.

Comments: 13 pages, 1 table, 4 figures. Accepted in MNRAS

The Impact of Baryonic Physics on the Structure of Dark Matter Halos: the View from the FIRE Cosmological Simulations

Authors: T. K. Chan, D. Kereš, J. Oñorbe, P. F. Hopkins, A. L. Muratov, C. -A. Faucher-Giguère, E. Quataert

Abstract: We study the distribution of cold dark matter (CDM) in cosmological simulations from the FIRE (Feedback In Realistic Environments) project, for $M_{\ast}\sim10^{4-11}\,M_{\odot}$ galaxies in $M_{\rm h}\sim10^{9-12}\,M_{\odot}$ halos. FIRE incorporates explicit stellar feedback in the multi-phase ISM, with energetics from stellar population models. We find that stellar feedback, without "fine-tuned… ▽ More We study the distribution of cold dark matter (CDM) in cosmological simulations from the FIRE (Feedback In Realistic Environments) project, for $M_{\ast}\sim10^{4-11}\,M_{\odot}$ galaxies in $M_{\rm h}\sim10^{9-12}\,M_{\odot}$ halos. FIRE incorporates explicit stellar feedback in the multi-phase ISM, with energetics from stellar population models. We find that stellar feedback, without "fine-tuned" parameters, greatly alleviates small-scale problems in CDM. Feedback causes bursts of star formation and outflows, altering the DM distribution. As a result, the inner slope of the DM halo profile ($α$) shows a strong mass dependence: profiles are shallow at $M_{\rm h}\sim10^{10}-10^{11}\,M_{\odot}$ and steepen at higher/lower masses. The resulting core sizes and slopes are consistent with observations. This is broadly consistent with previous work using simpler feedback schemes, but we find steeper mass dependence of $α$, and relatively late growth of cores. Because the star formation efficiency $M_{\ast}/M_{\rm h}$ is strongly halo mass dependent, a rapid change in $α$ occurs around $M_{\rm h}\sim 10^{10}\,M_{\odot}$ ($M_{\ast}\sim10^{6}-10^{7}\,M_{\odot}$), as sufficient feedback energy becomes available to perturb the DM. Large cores are not established during the period of rapid growth of halos because of ongoing DM mass accumulation. Instead, cores require several bursts of star formation after the rapid buildup has completed. Stellar feedback dramatically reduces circular velocities in the inner kpc of massive dwarfs; this could be sufficient to explain the "Too Big To Fail" problem without invoking non-standard DM. Finally, feedback and baryonic contraction in Milky Way-mass halos produce DM profiles slightly shallower than the Navarro-Frenk-White profile, consistent with the normalization of the observed Tully-Fisher relation. △ Less

Submitted 23 September, 2015; v1 submitted 8 July, 2015; originally announced July 2015.

Comments: 23 pages, 14 figures. MNRAS, in press; minor changes to match the accepted version

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

Authors: Camilla Penzo, Andrea V. Macciò, Marco Baldi, Luciano Casarini, Jose Oñorbe

Abstract: 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… ▽ More 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. △ Less

Submitted 27 April, 2015; originally announced April 2015.

Comments: 12 pages, 16 figures, submitted to MNRAS

Lagrangian Volume Deformations around Simulated Galaxies

Authors: S. Robles, R. Domínguez-Tenreiro, J. Oñorbe, F. J. Martínez-Serrano

Abstract: We present a detailed analysis of the local evolution of 206 Lagrangian Volumes (LVs) selected at high redshift around galaxy seeds, identified in a large-volume $Λ$ cold dark matter ($Λ$CDM) hydrodynamical simulation. The LVs have a mass range of $1 - 1500 \times 10^{10} M_\odot$. We follow the dynamical evolution of the density field inside these initially spherical LVs from $z=10$ up to… ▽ More We present a detailed analysis of the local evolution of 206 Lagrangian Volumes (LVs) selected at high redshift around galaxy seeds, identified in a large-volume $Λ$ cold dark matter ($Λ$CDM) hydrodynamical simulation. The LVs have a mass range of $1 - 1500 \times 10^{10} M_\odot$. We follow the dynamical evolution of the density field inside these initially spherical LVs from $z=10$ up to $z_{\rm low} = 0.05$, witnessing highly non-linear, anisotropic mass rearrangements within them, leading to the emergence of the local cosmic web (CW). These mass arrangements have been analysed in terms of the reduced inertia tensor $I_{ij}^r$, focusing on the evolution of the principal axes of inertia and their corresponding eigendirections, and paying particular attention to the times when the evolution of these two structural elements declines. In addition, mass and component effects along this process have also been investigated. We have found that deformations are led by dark matter dynamics and they transform most of the initially spherical LVs into prolate shapes, i.e. filamentary structures. An analysis of the individual freezing-out time distributions for shapes and eigendirections shows that first most of the LVs fix their three axes of symmetry (like a skeleton) early on, while accretion flows towards them still continue. Very remarkably, we have found that more massive LVs fix their skeleton earlier on than less massive ones. We briefly discuss the astrophysical implications our findings could have, including the galaxy mass-morphology relation and the effects on the galaxy-galaxy merger parameter space, among others. △ Less

Submitted 8 June, 2015; v1 submitted 23 April, 2015; originally announced April 2015.

Comments: 23 pages, 20 figures. Minor editorial improvements

Journal ref: MNRAS 451, 5005-5026 (2015)

Sweating the small stuff: simulating dwarf galaxies, ultra-faint dwarf galaxies, and their own tiny satellites

Authors: Coral Wheeler, Jose Onorbe, James S. Bullock, Michael Boylan-Kolchin, Oliver D. Elbert, Shea Garrison-Kimmel, Philip F. Hopkins, Dusan Keres

Abstract: We present FIRE/Gizmo hydrodynamic zoom-in simulations of isolated dark matter halos, two each at the mass of classical dwarf galaxies ($M_{\rm vir} \simeq 10^{10} M_{\odot}$) and ultra-faint galaxies ($M_{\rm vir} \simeq 10^9 M_{\odot}$), and with two feedback implementations. The resultant central galaxies lie on an extrapolated abundance matching relation from $M_{\star} \simeq 10^6$ to… ▽ More We present FIRE/Gizmo hydrodynamic zoom-in simulations of isolated dark matter halos, two each at the mass of classical dwarf galaxies ($M_{\rm vir} \simeq 10^{10} M_{\odot}$) and ultra-faint galaxies ($M_{\rm vir} \simeq 10^9 M_{\odot}$), and with two feedback implementations. The resultant central galaxies lie on an extrapolated abundance matching relation from $M_{\star} \simeq 10^6$ to $10^4 M_{\odot}$ without a break. Every host is filled with subhalos, many of which form stars. Our dwarfs with $M_{\star} \simeq 10^6 M_{\odot}$ each have 1-2 well-resolved satellites with $M_{\star} = 3-200 \times 10^3 M_{\odot}$. Even our isolated ultra-faint galaxies have star-forming subhalos. If this is representative, dwarf galaxies throughout the universe should commonly host tiny satellite galaxies of their own. We combine our results with the ELVIS simulations to show that targeting $\sim 50~ \rm kpc$ regions around nearby isolated dwarfs could increase the chances of discovering ultra-faint galaxies by $\sim 35\%$ compared to random halo pointings, and specifically identify the region around the Phoenix dwarf galaxy as a good potential target. The well-resolved ultra-faint galaxies in our simulations ($M_{\star} \simeq 3 - 30 \times 10^3 M_{\odot}$) form within $M_{\rm peak} \simeq 0.5 - 3 \times 10^9 M_{\odot}$ halos. Each has a uniformly ancient stellar population ($ > 10~ \rm Gyr$) owing to reionization-related quenching. More massive systems, in contrast, all have late-time star formation. Our results suggest that $M_{\rm halo} \simeq 5 \times 10^9 M_{\odot}$ is a probable dividing line between halos hosting reionization "fossils" and those hosting dwarfs that can continue to form stars in isolation after reionization. △ Less

Submitted 6 May, 2015; v1 submitted 9 April, 2015; originally announced April 2015.

Comments: 12 pages, 6 figures, 1 table, submitted to MNRAS

Characterizing the Pressure Smoothing Scale of the Intergalactic Medium

Authors: Girish Kulkarni, Joseph F. Hennawi, Jose Oñorbe, Alberto Rorai, Volker Springel

Abstract: 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 hydrod… ▽ More 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$. △ Less

Submitted 13 October, 2015; v1 submitted 1 April, 2015; originally announced April 2015.

Comments: Now published. New from previous version: Figure 14 and some nomenclature

Journal ref: ApJ (2015) 812 30

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

Authors: Jose Oñorbe, Michael Boylan-Kolchin, James S. Bullock, Philip F. Hopkins, Dušan Kerěs, Claude-André Faucher-Giguère, Eliot Quataert, Norman Murray

Abstract: 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 inter… ▽ More 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. △ Less

Submitted 17 October, 2015; v1 submitted 6 February, 2015; originally announced February 2015.

Comments: 16 pages, 13 figures, accepted in MNRAS

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

Authors: Oliver D. Elbert, James S. Bullock, Shea Garrison-Kimmel, Miguel Rocha, Jose Oñorbe, Annika H. G. Peter

Abstract: 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… ▽ More 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. △ Less

Submitted 3 December, 2014; originally announced December 2014.

Comments: 9 pages, 7 figures

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

Authors: Christian Schultz, Jose Oñorbe, Kevork N. Abazajian, James S. Bullock

Abstract: 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… ▽ More 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. △ Less

Submitted 14 May, 2014; v1 submitted 15 January, 2014; originally announced January 2014.

Journal ref: Mon. Not. Roy. Astr. Soc. 442, p.1597-1609 (2014)

Galaxies on FIRE (Feedback In Realistic Environments): Stellar Feedback Explains Cosmologically Inefficient Star Formation

Authors: Philip F. Hopkins, Dusan Keres, Jose Onorbe, Claude-Andre Faucher-Giguere, Eliot Quataert, Norm Murray, James S. Bullock

Abstract: We present a series of high-resolution cosmological simulations of galaxy formation to z=0, spanning halo masses ~10^8-10^13 M_sun, and stellar masses ~10^4-10^11. Our simulations include fully explicit treatment of both the multi-phase ISM (molecular through hot) and stellar feedback. The stellar feedback inputs (energy, momentum, mass, and metal fluxes) are taken directly from stellar population… ▽ More We present a series of high-resolution cosmological simulations of galaxy formation to z=0, spanning halo masses ~10^8-10^13 M_sun, and stellar masses ~10^4-10^11. Our simulations include fully explicit treatment of both the multi-phase ISM (molecular through hot) and stellar feedback. The stellar feedback inputs (energy, momentum, mass, and metal fluxes) are taken directly from stellar population models. These sources of stellar feedback, with zero adjusted parameters, reproduce the observed relation between stellar and halo mass up to M_halo~10^12 M_sun (including dwarfs, satellites, MW-mass disks, and small groups). By extension, this leads to reasonable agreement with the stellar mass function for M_star<10^11 M_sun. We predict weak redshift evolution in the M_star-M_halo relation, consistent with current constraints to z>6. We find that the M_star-M_halo relation is insensitive to numerical details, but is sensitive to the feedback physics. Simulations with only supernova feedback fail to reproduce the observed stellar masses, particularly in dwarf and high-redshift galaxies: radiative feedback (photo-heating and radiation pressure) is necessary to disrupt GMCs and enable efficient coupling of later supernovae to the gas. Star formation rates agree well with the observed Kennicutt relation at all redshifts. The galaxy-averaged Kennicutt relation is very different from the numerically imposed law for converting gas into stars in the simulation, and is instead determined by self-regulation via stellar feedback. Feedback reduces star formation rates considerably and produces a reservoir of gas that leads to rising late-time star formation histories significantly different from the halo accretion history. Feedback also produces large short-timescale variability in galactic SFRs, especially in dwarfs. Many of these properties are not captured by common 'sub-grid' galactic wind models. △ Less

Submitted 15 September, 2014; v1 submitted 8 November, 2013; originally announced November 2013.

Comments: 23 pages, 13 figures, accepted to MNRAS. Revised to match published version. For movies of the simulations here, see http://www.tapir.caltech.edu/~phopkins/Site/Movies_cosmo.html or the FIRE project site http://fire.northwestern.edu

Journal ref: MNRAS, 445, 581 (2014)

Sterile neutrino dark matter bounds from galaxies of the Local Group

Authors: Shunsaku Horiuchi, Philip J. Humphrey, Jose Onorbe, Kevork N. Abazajian, Manoj Kaplinghat, Shea Garrison-Kimmel

Abstract: We show that the canonical oscillation-based (non-resonant) production of sterile neutrino dark matter is inconsistent at $>99$% confidence with observations of galaxies in the Local Group. We set lower limits on the non-resonant sterile neutrino mass of $2.5$ keV (equivalent to $0.7$ keV thermal mass) using phase-space densities derived for dwarf satellite galaxies of the Milky Way, as well as li… ▽ More We show that the canonical oscillation-based (non-resonant) production of sterile neutrino dark matter is inconsistent at $>99$% confidence with observations of galaxies in the Local Group. We set lower limits on the non-resonant sterile neutrino mass of $2.5$ keV (equivalent to $0.7$ keV thermal mass) using phase-space densities derived for dwarf satellite galaxies of the Milky Way, as well as limits of $8.8$ keV (equivalent to $1.8$ keV thermal mass) based on subhalo counts of $N$-body simulations of M 31 analogues. Combined with improved upper mass limits derived from significantly deeper X-ray data of M 31 with full consideration for background variations, we show that there remains little room for non-resonant production if sterile neutrinos are to explain $100$% of the dark matter abundance. Resonant and non-oscillation sterile neutrino production remain viable mechanisms for generating sufficient dark matter sterile neutrinos. △ Less

Submitted 1 November, 2013; originally announced November 2013.

Comments: 10 pages, 4 figures, 2 tables. Submitted to PRD

Journal ref: Phys. Rev. D 89, 025017 (2014)

The AGORA High-Resolution Galaxy Simulations Comparison Project

Authors: Ji-hoon Kim, Tom Abel, Oscar Agertz, Greg L. Bryan, Daniel Ceverino, Charlotte Christensen, Charlie Conroy, Avishai Dekel, Nickolay Y. Gnedin, Nathan J. Goldbaum, Javiera Guedes, Oliver Hahn, Alexander Hobbs, Philip F. Hopkins, Cameron B. Hummels, Francesca Iannuzzi, Dusan Keres, Anatoly Klypin, Andrey V. Kravtsov, Mark R. Krumholz, Michael Kuhlen, Samuel N. Leitner, Piero Madau, Lucio Mayer, Christopher E. Moody , et al. (21 additional authors not shown)

Abstract: We introduce the AGORA project, a comprehensive numerical study of well-resolved galaxies within the LCDM cosmology. Cosmological hydrodynamic simulations with force resolutions of ~100 proper pc or better will be run with a variety of code platforms to follow the hierarchical growth, star formation history, morphological transformation, and the cycle of baryons in and out of 8 galaxies with halo… ▽ More We introduce the AGORA project, a comprehensive numerical study of well-resolved galaxies within the LCDM cosmology. Cosmological hydrodynamic simulations with force resolutions of ~100 proper pc or better will be run with a variety of code platforms to follow the hierarchical growth, star formation history, morphological transformation, and the cycle of baryons in and out of 8 galaxies with halo masses M_vir ~= 1e10, 1e11, 1e12, and 1e13 Msun at z=0 and two different ("violent" and "quiescent") assembly histories. The numerical techniques and implementations used in this project include the smoothed particle hydrodynamics codes GADGET and GASOLINE, and the adaptive mesh refinement codes ART, ENZO, and RAMSES. The codes will share common initial conditions and common astrophysics packages including UV background, metal-dependent radiative cooling, metal and energy yields of supernovae, and stellar initial mass function. These are described in detail in the present paper. Subgrid star formation and feedback prescriptions will be tuned to provide a realistic interstellar and circumgalactic medium using a non-cosmological disk galaxy simulation. Cosmological runs will be systematically compared with each other using a common analysis toolkit, and validated against observations to verify that the solutions are robust - i.e., that the astrophysical assumptions are responsible for any success, rather than artifacts of particular implementations. The goals of the AGORA project are, broadly speaking, to raise the realism and predictive power of galaxy simulations and the understanding of the feedback processes that regulate galaxy "metabolism." The proof-of-concept dark matter-only test of the formation of a galactic halo with a z=0 mass of M_vir ~= 1.7e11 Msun by 9 different versions of the participating codes is also presented to validate the infrastructure of the project. △ Less

Submitted 24 December, 2013; v1 submitted 12 August, 2013; originally announced August 2013.

Comments: 21 pages, 6 figures, Accepted for publication in the Astrophysical Journal Supplement, Image resolution greatly reduced, High-resolution version of this article and more information about the AGORA Project including the science goals of the 13 Working Groups are available at http://www.AGORAsimulations.org/ and at http://sites.google.com/site/santacruzcomparisonproject/details/

Journal ref: Astrophysical Journal Supplement Series, 210, 14, 2014

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

Authors: Jose Oñorbe, Shea Garrison-Kimmel, Ariyeh H. Maller, James S. Bullock, Miguel Rocha, Oliver Hahn

Abstract: 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 ch… ▽ More 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. △ Less

Submitted 19 October, 2013; v1 submitted 29 May, 2013; originally announced May 2013.

Comments: 15 pages, 13 Figures, all figures include colors. MNRAS accepted

Cosmological Simulations with Self-Interacting Dark Matter I: Constant Density Cores and Substructure

Authors: Miguel Rocha, Annika H. G. Peter, James S. Bullock, Manoj Kaplinghat, Shea Garrison-Kimmel, Jose Onorbe, Leonidas A. Moustakas

Abstract: We use cosmological simulations to study the effects of self-interacting dark matter (SIDM) on the density profiles and substructure counts of dark matter halos from the scales of spiral galaxies to galaxy clusters, focusing explicitly on models with cross sections over dark matter particle mass σ/m = 1 and 0.1 cm^2/g. Our simulations rely on a new SIDM N-body algorithm that is derived self-consis… ▽ More We use cosmological simulations to study the effects of self-interacting dark matter (SIDM) on the density profiles and substructure counts of dark matter halos from the scales of spiral galaxies to galaxy clusters, focusing explicitly on models with cross sections over dark matter particle mass σ/m = 1 and 0.1 cm^2/g. Our simulations rely on a new SIDM N-body algorithm that is derived self-consistently from the Boltzmann equation and that reproduces analytic expectations in controlled numerical experiments. We find that well-resolved SIDM halos have constant-density cores, with significantly lower central densities than their CDM counterparts. In contrast, the subhalo content of SIDM halos is only modestly reduced compared to CDM, with the suppression greatest for large hosts and small halo-centric distances. Moreover, the large-scale clustering and halo circular velocity functions in SIDM are effectively identical to CDM, meaning that all of the large-scale successes of CDM are equally well matched by SIDM. From our largest cross section runs we are able to extract scaling relations for core sizes and central densities over a range of halo sizes and find a strong correlation between the core radius of an SIDM halo and the NFW scale radius of its CDM counterpart. We construct a simple analytic model, based on CDM scaling relations, that captures all aspects of the scaling relations for SIDM halos. Our results show that halo core densities in σ/m = 1 cm^2/g models are too low to match observations of galaxy clusters, low surface brightness spirals (LSBs), and dwarf spheroidal galaxies. However, SIDM with σ/m ~ 0.1 cm^2/g appears capable of reproducing reported core sizes and central densities of dwarfs, LSBs, and galaxy clusters without the need for velocity dependence. (abridged) △ Less

Submitted 15 August, 2012; originally announced August 2012.

Comments: 26 pages, 16 figures, all figures include colors, submitted for publication in MNRAS

Massive Galaxies at High-z: Assembly Patterns, Structure & Dynamics in the Fast Phase of Galaxy Formation

Authors: J. Oñorbe, F. J. Martínez-Serrano, R. Domínguez-Tenreiro, A. Knebe, A. Serna

Abstract: Relaxed, massive galactic objects have been identified at redshifts z = 4;5; and 6 in hydrodynamical simulations run in a large cosmological volume. This allowed us to analyze the assembly patterns of the high mass end of the galaxy distribution at these high zs, by focusing on their structural and dynamical properties. Our simulations indicate that massive objects at high redshift already follow… ▽ More Relaxed, massive galactic objects have been identified at redshifts z = 4;5; and 6 in hydrodynamical simulations run in a large cosmological volume. This allowed us to analyze the assembly patterns of the high mass end of the galaxy distribution at these high zs, by focusing on their structural and dynamical properties. Our simulations indicate that massive objects at high redshift already follow certain scaling relations. These relations define virial planes at the halo scale, whereas at the galactic scale they define intrinsic dynamical planes that are, however, tilted relative to the virial plane. Therefore, we predict that massive galaxies must lie on fundamental planes from their formation. We briefly discuss the physical origin of the tilt in terms the physical processes underlying massive galaxy formation at high z, in the context of a two-phase galaxy formation scenario. Specifically, we have found that it lies on the different behavior of the gravitationally heated gas as compared with cold gas previously involved in caustic formation, and the mass dependence of the energy available to heat the gas. △ Less

Submitted 22 March, 2011; originally announced March 2011.

Comments: 12 pages, 1 figure. Accepted to Astrophysical Journal Letters

Journal ref: Astrophys.J. 732 (2011) L32

Large-Scale Gas Dynamics in the Adhesion Model: Implications for the Two-Phase Massive Galaxy Formation Scenario

Authors: R. Domínguez-Tenreiro, J. Oñorbe, F. Martínez-Serrano, A. Serna

Abstract: The mass assembly and star formation histories of massive galaxies identified at low redshift z in different cosmological hydrodynamical simulations, have been studied through a detailed follow-up backwards in time of their constituent mass elements (sampled by particles) of different types. Then, the configurations they depict at progressively higher zs have been analysed. The analyses show tha… ▽ More The mass assembly and star formation histories of massive galaxies identified at low redshift z in different cosmological hydrodynamical simulations, have been studied through a detailed follow-up backwards in time of their constituent mass elements (sampled by particles) of different types. Then, the configurations they depict at progressively higher zs have been analysed. The analyses show that these histories share common generic patterns, irrespective of particular circumstances. In any case, the results we have found are different depending on the particle type. The most outstanding differences follow. We have found that by z ~ 3.5 - 6, mass elements identified as stellar particles at z=0 exhibit a gaseous cosmic-web-like morphology with scales of ~ 1 physical Mpc, where the densest mass elements have already turned into stars by z ~ 6. These settings are in fact the densest pieces of the cosmic web, where no hot particles show up, and dynamically organized as a hierarchy of flow convergence regions, that is, attraction basins for mass flows. On the other hand, mass elements identified at the diffuse hot coronae surrounding massive galaxies at z = 0, do not display a clear web-like morphology at any z. Diffuse gas is heated when flow convergence regions go through contractive deformations, and most of it keeps hot and with low density along the evolution. To shed light on the physical foundations of the behaviour our analyses show up, as well as on their possible observational implications, these patterns have been confronted with some generic properties of singular flows as described by the adhesion model. We have found that these common patterns simulations show can be interpreted as a consequence of flow properties, that, moreover, could explain different generic observational results on massive galaxies or their samples. We briefly discuss some of them.[Abridged] △ Less

Submitted 23 March, 2011; v1 submitted 29 October, 2010; originally announced October 2010.

Comments: Aceppted to MNRAS. 28 pages, 8 Figures. Minor changes to match the published version

Journal ref: MNRAS 413 (2011) 3022