-
scida: scalable analysis for scientific big data
Authors:
Chris Byrohl,
Dylan Nelson
Abstract:
scida is a Python package for reading and analyzing large scientific data sets with support for various cosmological and galaxy formation simulations out-of-the-box. Data access is provided through a hierarchical dictionary-like data structure after a simple load() function. Using the dask library for scalable, parallel and out-of-core computation, all computation requests from a user session are…
▽ More
scida is a Python package for reading and analyzing large scientific data sets with support for various cosmological and galaxy formation simulations out-of-the-box. Data access is provided through a hierarchical dictionary-like data structure after a simple load() function. Using the dask library for scalable, parallel and out-of-core computation, all computation requests from a user session are first collected in a task graph. Arbitrary custom analysis, as well as all available dask (array) operations, can be performed. The subsequent computation is executed only upon request, on a target resource (e.g. a HPC cluster).
△ Less
Submitted 27 February, 2024;
originally announced February 2024.
-
Absorption Troughs of Lyman Alpha Emitters in HETDEX
Authors:
Laurel H. Weiss,
Dustin Davis,
Karl Gebhardt,
Simon Gazagnes,
Mahan Mirza Khanlari,
Erin Mentuch Cooper,
John Chisholm,
Danielle Berg,
William P. Bowman,
Chris Byrohl,
Robin Ciardullo,
Maximilian Fabricius,
Daniel Farrow,
Caryl Gronwall,
Gary J. Hill,
Lindsay R. House,
Donghui Jeong,
Hasti Khoraminezhad,
Wolfram Kollatschny,
Eiichiro Komatsu,
Maja Lujan Niemeyer,
Shun Saito,
Donald P. Schneider,
Gregory R. Zeimann
Abstract:
The Hobby-Eberly Telescope Dark Energy Experiment (HETDEX) is designed to detect and measure the redshifts of more than one million Ly$α$ emitting galaxies (LAEs) between $1.88 < z < 3.52$. In addition to its cosmological measurements, these data enable studies of Ly$α$ spectral profiles and the underlying radiative transfer. Using the roughly half a million LAEs in the HETDEX Data Release 3, we s…
▽ More
The Hobby-Eberly Telescope Dark Energy Experiment (HETDEX) is designed to detect and measure the redshifts of more than one million Ly$α$ emitting galaxies (LAEs) between $1.88 < z < 3.52$. In addition to its cosmological measurements, these data enable studies of Ly$α$ spectral profiles and the underlying radiative transfer. Using the roughly half a million LAEs in the HETDEX Data Release 3, we stack various subsets to obtain the typical Ly$α$ profile for the $z \sim 2-3$ epoch and to understand their physical properties. We find clear absorption wings around Ly$α$ emission, which extend $\sim 2000$ km $\mathrm{s}^{-1}$ both redward and blueward of the central line. Using far-UV spectra of nearby ($0.002 < z < 0.182$) LAEs in the CLASSY treasury and optical/near-IR spectra of $2.8 < z < 6.7$ LAEs in the MUSE-Wide survey, we observe absorption profiles in both redshift regimes. Dividing the sample by volume density shows that the troughs increase in higher density regions. This trend suggests that the depth of the absorption is dependent on the local density of objects near the LAE, a geometry that is similar to damped Lyman-$α$ systems. Simple simulations of Ly$α$ radiative transfer can produce similar troughs due to absorption of light from background sources by HI gas surrounding the LAEs.
△ Less
Submitted 4 January, 2024;
originally announced January 2024.
-
Resonant scattering of the OVII X-ray emission line in the circumgalactic medium of TNG50 galaxies
Authors:
Dylan Nelson,
Chris Byrohl,
Anna Ogorzalek,
Maxim Markevitch,
Ildar Khabibullin,
Eugene Churazov,
Irina Zhuravleva,
Akos Bogdan,
Priyanka Chakraborty,
Caroline Kilbourne,
Ralph Kraft,
Annalisa Pillepich,
Arnab Sarkar,
Gerrit Schellenberger,
Yuanyuan Su,
Nhut Truong,
Stephan Vladutescu-Zopp,
Nastasha Wijers
Abstract:
We study the impact of resonantly scattered X-ray line emission on the observability of the hot circumgalactic medium (CGM) of galaxies. We apply a Monte Carlo radiative transfer post-processing analysis to the high-resolution TNG50 cosmological magnetohydrodynamical galaxy formation simulation. This allows us to model the resonant scattering of OVII(r) X-ray photons within the complex, multi-phas…
▽ More
We study the impact of resonantly scattered X-ray line emission on the observability of the hot circumgalactic medium (CGM) of galaxies. We apply a Monte Carlo radiative transfer post-processing analysis to the high-resolution TNG50 cosmological magnetohydrodynamical galaxy formation simulation. This allows us to model the resonant scattering of OVII(r) X-ray photons within the complex, multi-phase, multi-scale CGM. The resonant transition of the OVII He-like triplet is one of the brightest, and most promising, X-ray emission lines for detecting the hot CGM and measuring its physical properties. We focus on galaxies with stellar masses 10 < log(M*/Msun) < 11 at z ~ 0. After constructing a model for OVII(r) emission from the central galaxy as well as from CGM gas, we forward model these intrinsic photons to derive observable surface brightness maps. We find that scattering significantly boosts the observable OVII(r) surface brightness of the extended and diffuse CGM. This enhancement can be large -- an order of magnitude on average at a distance of 200 projected kpc for high-mass M* = 10^10.7 Msun galaxies. The enhancement is larger for lower mass galaxies, and can even reach a factor of 100, across the extended CGM. Galaxies with higher star formation rates, AGN luminosities, and central OVII(r) luminosities all have larger scattering enhancements, at fixed stellar mass. Our results suggest that next-generation X-ray spectroscopic missions including XRISM, LEM, ATHENA, and HUBS -- which aim to detect the hot CGM in emission -- could specifically target halos with significant enhancements due to resonant scattering.
△ Less
Submitted 8 June, 2023;
originally announced June 2023.
-
The cosmic web in Lyman-alpha emission
Authors:
Chris Byrohl,
Dylan Nelson
Abstract:
We develop a comprehensive theoretical model for Lyman-alpha emission, from the scale of individual Lyman-alpha emitters (LAEs) to Lyman-alpha halos (LAHs), Lyman-alpha blobs (LABs), and Lyman-alpha filaments (LAFs) of the diffuse cosmic web itself. To do so, we post-process the high-resolution TNG50 cosmological magnetohydrodynamical simulation with a Monte Carlo radiative transfer method to capt…
▽ More
We develop a comprehensive theoretical model for Lyman-alpha emission, from the scale of individual Lyman-alpha emitters (LAEs) to Lyman-alpha halos (LAHs), Lyman-alpha blobs (LABs), and Lyman-alpha filaments (LAFs) of the diffuse cosmic web itself. To do so, we post-process the high-resolution TNG50 cosmological magnetohydrodynamical simulation with a Monte Carlo radiative transfer method to capture the resonant scattering process of Lyman-alpha photons. We build an emission model incorporating recombinations and collisions in diffuse gas, including radiative effects from nearby AGN, as well as emission sourced by stellar populations. Our treatment includes a physically motivated dust model, which we empirically calibrate to the observed LAE luminosity function. We then focus on the observability, and physical origin, of the $z=2$ Lyman-alpha cosmic web, studying the dominant emission mechanisms and spatial origins. We find that diffuse Lyman-alpha filaments are, in fact, illuminated by photons which originate, not from the intergalactic medium itself, but from within galaxies and their gaseous halos. In our model, this emission is primarily sourced by intermediate mass halos ($10^{10} - 10^{11}\,$M$_{\odot}$), principally due to collisional excitations in their circumgalactic media as well as central, young stellar populations. Observationally, we make predictions for the abundance, area, linear size, and embedded halo/emitter populations within filaments. Adopting an isophotal surface brightness threshold of $10^{-20}\,$erg$\,$s$^{-1}\,$cm$^{-2}\,$arcsec$^{-2}$, we predict a volume abundance of Lyman-alpha filaments of ${\sim}10^{-3}$ cMpc$^{-3}\,$ for lengths above $400\,$pkpc. Given sufficiently large survey footprints, detection of the Lyman-alpha cosmic web is within reach of modern integral field spectrographs, including MUSE, VIRUS, and KCWI.
△ Less
Submitted 16 December, 2022;
originally announced December 2022.
-
Surface Brightness Profile of Lyman-$α$ Halos out to 320 kpc in HETDEX
Authors:
Maja Lujan Niemeyer,
Eiichiro Komatsu,
Chris Byrohl,
Dustin Davis,
Maximilian Fabricius,
Karl Gebhardt,
Gary J. Hill,
Lutz Wisotzki,
William P. Bowman,
Robin Ciardullo,
Daniel J. Farrow,
Steven L. Finkelstein,
Eric Gawiser,
Caryl Gronwall,
Donghui Jeong,
Martin Landriau,
Chenxu Liu,
Erin Mentuch Cooper,
Masami Ouchi,
Donald P. Schneider,
Gregory R. Zeimann
Abstract:
We present the median-stacked Lyman-$α$ surface brightness profile of 968 spectroscopically selected Lyman-$α$ emitting galaxies (LAEs) at redshifts $1.9<z<3.5$ in the early data of the Hobby-Eberly Telescope Dark Energy Experiment (HETDEX). The selected LAEs are high-confidence Lyman-$α$ detections with large signal-to-noise ratios observed with good seeing conditions (point-spread-function full-…
▽ More
We present the median-stacked Lyman-$α$ surface brightness profile of 968 spectroscopically selected Lyman-$α$ emitting galaxies (LAEs) at redshifts $1.9<z<3.5$ in the early data of the Hobby-Eberly Telescope Dark Energy Experiment (HETDEX). The selected LAEs are high-confidence Lyman-$α$ detections with large signal-to-noise ratios observed with good seeing conditions (point-spread-function full-width-at-half-maximum $<1.4"$), excluding active galactic nuclei (AGN). The Lyman-$α$ luminosities of the LAEs are $10^{42.4}-10^{43}\, \mathrm{erg}\, \mathrm{s}^{-1}$. We detect faint emission in the median-stacked radial profiles at the level of $(3.6\pm 1.3)\times 10^{-20}\,\mathrm{erg}\,\mathrm{s}^{-1}\,\mathrm{cm}^{-2}\,\mathrm{arcsec}^{-2}$ from the surrounding Lyman-$α$ halos out to $r\simeq 160$ kpc (physical). The shape of the median-stacked radial profile is consistent at $r<80\,\mathrm{kpc}$ with that of much fainter LAEs at $3<z<4$ observed with the Multi Unit Spectroscopic Explorer (MUSE), indicating that the median-stacked Lyman-$α$ profiles have similar shapes at redshifts $2<z<4$ and across a factor of $10$ in Lyman-$α$ luminosity. While we agree with the results from the MUSE sample at $r<80\,\mathrm{kpc}$, we extend the profile over a factor of two in radius. At $r>80\,\mathrm{kpc}$, our profile is flatter than the MUSE model. The measured profile agrees at most radii with that of galaxies in the Byrohl et al. (2021) cosmological radiative transfer simulation at $z=3$. This suggests that the surface brightness of a Lyman-$α$ halo at $r\lesssim 100$ kpc is dominated by resonant scattering of Lyman-$α$ photons from star-forming regions in the central galaxy, whereas at $r > 100$ kpc it is dominated by photons from galaxies in surrounding dark matter halos.
△ Less
Submitted 9 March, 2022;
originally announced March 2022.
-
The cold circumgalactic medium in emission: MgII halos in TNG50
Authors:
Dylan Nelson,
Chris Byrohl,
Celine Peroux,
Kate H. R. Rubin,
Joseph N. Burchett
Abstract:
We outline theoretical predictions for extended emission from MgII, tracing cool ~10^4 K gas in the circumgalactic medium (CGM) of star-forming galaxies in the high-resolution TNG50 cosmological magnetohydrodynamical simulation. We synthesize surface brightness maps of this strong rest-frame ultraviolet metal emission doublet (2796, 2803), adopting the assumption that the resonant scattering of Mg…
▽ More
We outline theoretical predictions for extended emission from MgII, tracing cool ~10^4 K gas in the circumgalactic medium (CGM) of star-forming galaxies in the high-resolution TNG50 cosmological magnetohydrodynamical simulation. We synthesize surface brightness maps of this strong rest-frame ultraviolet metal emission doublet (2796, 2803), adopting the assumption that the resonant scattering of MgII can be neglected and connecting to recent and upcoming observations with the Keck/KCWI, VLT/MUSE, and BlueMUSE optical integral field unit spectrographs. Studying galaxies with stellar masses 7.5 < log(M*/M_sun) < 11 at redshifts z=0.3, 0.7, 1 and 2 we find that extended MgII halos in emission, similar to their Lyman-alpha counterparts, are ubiquitous across the galaxy population. Median surface brightness profiles exceed 10^-19 erg/s/cm^2/arcsec^2 in the central ~10s of kpc, and total halo MgII luminosity increases with mass for star-forming galaxies, reaching 10^40 erg/s for M* ~ 10^9.5 Msun. MgII halo sizes increase from a few kpc to > 20 kpc at the highest masses, and sizes are larger for halos in denser environments. MgII halos are highly structured, clumpy, and asymmetric, with isophotal axis ratio increasing with galaxy mass. Similarly, the amount and distribution of MgII emission depends on the star formation activity of the central galaxy. Kinematically, inflowing versus outflowing gas dominates the MgII luminosity at high and low galaxy masses, respectively, although the majority of MgII halo emission at z~0.7 traces near-equilibrium fountain flows and gas with non-negligible rotational support, rather than rapidly outflowing galactic winds.
△ Less
Submitted 26 July, 2021; v1 submitted 16 June, 2021;
originally announced June 2021.
-
Cosmological 3D HI Gas Map with HETDEX Ly$α$ Emitters and eBOSS QSOs at $z=2$: IGM-Galaxy/QSO Connection and a $\sim$ 40-Mpc Scale Giant HII Bubble Candidate
Authors:
Shiro Mukae,
Masami Ouchi,
Gary J. Hill,
Karl Gebhardt,
Erin Mentuch Cooper,
Donghui Jeong,
Shun Saito,
Maximilian Fabricius,
Eric Gawiser,
Robin Ciardullo,
Daniel Farrow,
Dustin Davis,
Greg Zeimann,
Steven L. Finkelstein,
Caryl Gronwall,
Chenxu Liu,
Yechi Zhang,
Chris Byrohl,
Yoshiaki Ono,
Donald P. Schneider,
Matt J. Jarvis,
Caitlin M. Casey,
Ken Mawatari
Abstract:
We present cosmological ($30-400$ Mpc) distributions of neutral hydrogen (HI) in the inter-galactic medium (IGM) traced by Ly$α$ Emitters (LAEs) and QSOs at $z=2.1-2.5$, selected with the data of the on-going Hobby-Eberly Telescope Dark Energy Experiment (HETDEX) and the eBOSS survey. We investigate spatial correlations of LAEs and QSOs with HI tomography maps reconstructed from HI Ly$α$ forest ab…
▽ More
We present cosmological ($30-400$ Mpc) distributions of neutral hydrogen (HI) in the inter-galactic medium (IGM) traced by Ly$α$ Emitters (LAEs) and QSOs at $z=2.1-2.5$, selected with the data of the on-going Hobby-Eberly Telescope Dark Energy Experiment (HETDEX) and the eBOSS survey. We investigate spatial correlations of LAEs and QSOs with HI tomography maps reconstructed from HI Ly$α$ forest absorption in the spectra of background galaxies and QSOs obtained by the CLAMATO survey and this study, respectively. In the cosmological volume far from QSOs, we find that LAEs reside in regions of strong HI absorption, i.e. HI rich, which is consistent with results of previous galaxy-background QSO pair studies. Moreover, there is an anisotropy in the HI-distribution plot of transverse and line-of-sight distances; on average the HI absorption peak is blueshifted by $\sim 200$ km s$^{-1}$ from the LAE Ly$α$ redshift, reproducing the known average velocity offset between the Ly$α$ emission redshift and the galaxy systemic redshift. We have identified a $\sim$ 40-Mpc scale volume of HI underdensity that is a candidate for a giant HII bubble, where six QSOs and an LAE overdensity exist at $\left < z \right > =2.16$. The coincidence of the QSO and LAE overdensities with the HI underdensity indicates that the ionizing photon radiation of the QSOs has created a highly ionized volume of multiple proximity zones in a matter overdensity. Our results suggest an evolutionary picture where HI gas in an overdensity of galaxies becomes highly photoionized when QSOs emerge in the galaxies.
△ Less
Submitted 15 September, 2020;
originally announced September 2020.
-
The physical origins and dominant emission mechanisms of Lyman-alpha halos: results from the TNG50 simulation in comparison to MUSE observations
Authors:
Chris Byrohl,
Dylan Nelson,
Christoph Behrens,
Ivan Kostyuk,
Martin Glatzle,
Annalisa Pillepich,
Lars Hernquist,
Federico Marinacci,
Mark Vogelsberger
Abstract:
Extended Lyman-alpha emission is now commonly detected around high redshift galaxies through stacking and even on individual basis. Despite recent observational advances, the physical origin of these Lyman-alpha halos (LAHs), as well as their relationships to galaxies, quasars, circumgalactic gas, and other environmental factors remains unclear. We present results from our new Lyman-alpha full rad…
▽ More
Extended Lyman-alpha emission is now commonly detected around high redshift galaxies through stacking and even on individual basis. Despite recent observational advances, the physical origin of these Lyman-alpha halos (LAHs), as well as their relationships to galaxies, quasars, circumgalactic gas, and other environmental factors remains unclear. We present results from our new Lyman-alpha full radiative transfer code voroILTIS which runs directly on the unstructured Voronoi tessellation of cosmological hydrodynamical simulations. We make use of the TNG50 simulation and simulate LAHs from redshift $z=2$ to $z=5$, focusing on star-forming galaxies with $8.0 < \log_{10}{(M_\star/\rm{M}_\odot)} < 10.5$. While TNG50 does not directly follow ionizing radiation, it includes an on-the-fly treatment for active galactic nuclei and ultraviolet background radiation with self-shielding, which are important processes impacting the cooling and ionization of the gas. Based on this model, we present the predictions for the stacked radial surface brightness profiles of Ly$α$ as a function of galaxy mass and redshift. Comparison with data from the MUSE UDF at $z>3$ reveals a promising level of agreement. We measure the correlations of LAH size and central brightness with galaxy properties, finding that at the masses of $8.5 \leq \log_{10} \left(M_\star/\rm{M}_\odot\right) \leq 9.5$, physical LAH sizes roughly double from $z=2$ to $z=5$. Finally, we decompose the profiles into contributions from diffuse emission and scattered photons from star-forming regions. In our simulations, we find rescattered photons from star-forming regions to be the major source in observed LAHs. Unexpectedly, we find that the flattening of LAH profiles at large radii becomes dominated by photons originating from other nearby halos rather than diffuse emission itself.
△ Less
Submitted 5 October, 2021; v1 submitted 15 September, 2020;
originally announced September 2020.
-
Variations of Observed Lyman-α Spectra Shapes due to the Intergalactic Absorption
Authors:
Chris Byrohl,
Max Gronke
Abstract:
Lyman-$α$ (Ly$α$) spectra provide insights into the small-scale structure and kinematics of neutral hydrogen (HI) within galaxies as well as the ionization state of the intergalactic medium (IGM). The former defines the intrinsic spectrum of a galaxy, which is modified by the latter. These two effects are degenerate. Using the IllustrisTNG100 simulation, we study the impact of the IGM on Ly$α$ spe…
▽ More
Lyman-$α$ (Ly$α$) spectra provide insights into the small-scale structure and kinematics of neutral hydrogen (HI) within galaxies as well as the ionization state of the intergalactic medium (IGM). The former defines the intrinsic spectrum of a galaxy, which is modified by the latter. These two effects are degenerate. Using the IllustrisTNG100 simulation, we study the impact of the IGM on Ly$α$ spectral shapes between $z\sim 0$ and $5$. We compute the distribution of the expected Ly$α$ peaks and of the peak asymmetry for different intrinsic spectra, redshifts and large-scale environments. We find that the commonly used averaged transmission curves give an incorrect perception of the observed spectral properties. We show that the distributions of peak counts and asymmetry can lift the degeneracy between the intrinsic spectrum and the IGM absorption. For example, we expect a significant number of triple peaked Ly$α$ spectra (up to 30% at $z\sim 3$) if the galaxies' HI distribution become more porous at higher redshift as predicted by cosmological simulations. We provide a public catalog of transmission curves to be used in future simulations and observations to allow a more realistic IGM treatment.
△ Less
Submitted 17 June, 2020;
originally announced June 2020.
-
Radiative Transfer Distortions of Lyman-$α$ Emitters: a New Fingers-of-God Damping in the Clustering in Redshift Space
Authors:
Chris Byrohl,
Shun Saito,
Christoph Behrens
Abstract:
Complex radiative transfer (RT) of the Lyman-$α$ photons poses a theoretical challenge to galaxy surveys which infer the large-scale structure with Lyman-$α$ emitters (LAEs). Guided by RT simulations, prior studies investigated the impact of RT on the large-scale LAE clustering, and claimed that RT induces a selection effect which results in an anisotropic distortion even in real space but in an o…
▽ More
Complex radiative transfer (RT) of the Lyman-$α$ photons poses a theoretical challenge to galaxy surveys which infer the large-scale structure with Lyman-$α$ emitters (LAEs). Guided by RT simulations, prior studies investigated the impact of RT on the large-scale LAE clustering, and claimed that RT induces a selection effect which results in an anisotropic distortion even in real space but in an otherwise negligible effect in redshift space. However, our previous study, which relies on a full radiative transfer code run on the Illustris simulations, shows that the anisotropic selection effect was drastically reduced with higher spatial resolution. Adopting the same simulation framework, we further study the impact of RT on the LAE clustering in redshift space. Since we measure LAE's radial position through a spectral peak of a Lyman-$α$ emission, the frequency shift due to RT contaminates the redshift measurement and hence the inferred radial position in redshift space. We demonstrate that this additional RT offset suppresses the LAE clustering along the line of sight, which can be interpreted as a novel Fingers-of-God (FoG) effect. To assess the FoG effect, we develop a theoretical framework with particular emphasis on its connection with the underlying one-point velocity distribution which simultaneously takes into account the RT offset as well as the peculiar velocity that is commonly studied in the context of the Redshift Space Distortion (RSD). Although our findings strongly encourage a more careful RSD modeling in LAE surveys, we also seek a method to mitigate the additional FoG effect due to RT by making use of other information in a Lyman-$α$ spectrum.
△ Less
Submitted 6 June, 2019; v1 submitted 5 June, 2019;
originally announced June 2019.
-
The Intrinsic Stochasticity of the $^{56}$Ni Distribution of Single-Degenerate Type Ia Supernovae
Authors:
Chris Byrohl,
Robert T. Fisher,
Dean M. Townsley
Abstract:
Binary Chandrasekhar-mass white dwarfs accreting mass from non-degenerate stellar companions through the single-degenerate channel have reigned for decades as the leading explanation of Type Ia supernovae. Yet, a comprehensive theoretical explanation has not yet emerged to explain the expected properties of the canonical near-Chandrasekhar-mass white dwarf model. A simmering phase within the conve…
▽ More
Binary Chandrasekhar-mass white dwarfs accreting mass from non-degenerate stellar companions through the single-degenerate channel have reigned for decades as the leading explanation of Type Ia supernovae. Yet, a comprehensive theoretical explanation has not yet emerged to explain the expected properties of the canonical near-Chandrasekhar-mass white dwarf model. A simmering phase within the convective core of the white dwarf leads to the ignition of one or more flame bubbles scattered across the core. Consequently, near-Chandrasekhar-mass single-degenerate SNe Ia are inherently stochastic, and are expected to lead to a range of outcomes, from subluminous SN 2002cx-like events, to overluminous SN 1991T-like events. However, all prior simulations of the single-degenerate channel carried through the detonation phase have set the ignition points as free parameters. In this work, for the first time, we place ignition points as predicted by {\it ab initio} models of the convective phase leading up to ignition, and follow through the detonation phase in fully three-dimensional simulations. Single-degenerates in this framework are characteristically overluminous. Using a statistical approach, we determine the $^{56}$Ni mass distribution arising from stochastic ignition. While there is a total spread of $\gtrsim 0.2 M_{\odot}$ for detonating models, the distribution is strongly left-skewed, and with a narrow standard deviation of $\simeq 0.03 M_{\odot}$. Conversely, if single-degenerates are not overluminous but primarily yield normal or failed events, then the models require fine-tuning of the ignition parameters, or otherwise require revised physics or progenitor models. We discuss implications of our findings for the modeling of single-degenerate SNe Ia.
△ Less
Submitted 18 October, 2018;
originally announced October 2018.
-
The impact of Lyman-$α$ radiative transfer on large-scale clustering in the Illustris simulation
Authors:
Christoph Behrens,
Chris Byrohl,
Shun Saito,
Jens C. Niemeyer
Abstract:
Lyman-$α$ emitters (LAEs) are a promising probe of the large-scale structure at high redshift, $z\gtrsim 2$. In particular, the Hobby-Eberly Telescope Dark Energy Experiment aims at observing LAEs at 1.9 $<z<$ 3.5 to measure the Baryon Acoustic Oscillation (BAO) scale and the Redshift-Space Distortion (RSD). However, Zheng et al. (2011) pointed out that the complicated radiative transfer (RT) of t…
▽ More
Lyman-$α$ emitters (LAEs) are a promising probe of the large-scale structure at high redshift, $z\gtrsim 2$. In particular, the Hobby-Eberly Telescope Dark Energy Experiment aims at observing LAEs at 1.9 $<z<$ 3.5 to measure the Baryon Acoustic Oscillation (BAO) scale and the Redshift-Space Distortion (RSD). However, Zheng et al. (2011) pointed out that the complicated radiative transfer (RT) of the resonant Lyman-$α$ emission line generates an anisotropic selection bias in the LAE clustering on large scales, $s\gtrsim 10$ Mpc. This effect could potentially induce a systematic error in the BAO and RSD measurements. Also, Croft et al. (2016) claims an observational evidence of the effect in the Lyman-$α$ intensity map, albeit statistically insignificant. We aim at quantifying the impact of the Lyman-$α$ RT on the large-scale galaxy clustering in detail. For this purpose, we study the correlations between the large-scale environment and the ratio of an apparent Lyman-$α$ luminosity to an intrinsic one, which we call the `observed fraction', at $2<z<6$. We apply our Lyman-$α$ RT code by post-processing the full Illustris simulations. We simply assume that the intrinsic luminosity of the Lyman-$α$ emission is proportional to the star formation rate of galaxies in Illustris, yielding a sufficiently large sample of LAEs to measure the anisotropic selection bias. We find little correlations between large-scale environment and the observed fraction induced by the RT, and hence a smaller anisotropic selection bias than what was claimed by Zheng et al. (2011). We argue that the anisotropy was overestimated in the previous work due to the insufficient spatial resolution: it is important to keep the resolution such that it resolves the high density region down to the scale of the interstellar medium, $\sim1$ physical kpc. (abridged)
△ Less
Submitted 22 June, 2018; v1 submitted 17 October, 2017;
originally announced October 2017.
-
Viscosity, pressure, and support of the gas in simulations of merging cool-core clusters
Authors:
W. Schmidt,
C. Byrohl,
J. F. Engels,
C. Behrens,
J. C. Niemeyer
Abstract:
Major mergers are considered to be a significant source of turbulence in clusters. We performed a numerical simulation of a major merger event using nested-grid initial conditions, adaptive mesh refinement, radiative cooling of primordial gas, and a homogeneous ultraviolet background. By calculating the microscopic viscosity on the basis of various theoretical assumptions and estimating the Kolmog…
▽ More
Major mergers are considered to be a significant source of turbulence in clusters. We performed a numerical simulation of a major merger event using nested-grid initial conditions, adaptive mesh refinement, radiative cooling of primordial gas, and a homogeneous ultraviolet background. By calculating the microscopic viscosity on the basis of various theoretical assumptions and estimating the Kolmogorov length from the turbulent dissipation rate computed with a subgrid-scale model, we are able to demonstrate that most of the warm-hot intergalactic medium can sustain a fully turbulent state only if the magnetic suppression of the viscosity is considerable. Accepting this as premise, it turns out that ratios of turbulent and thermal quantities change only little in the course of the merger. This confirms the tight correlations between the mean thermal and non-thermal energy content for large samples of clusters in earlier studies, which can be interpreted as second self-similarity on top of the self-similarity for different halo masses. Another long-standing question is how and to which extent turbulence contributes to the support of the gas against gravity. From a global perspective, the ratio of turbulent and thermal pressures is significant for the clusters in our simulation. On the other hand, a local measure is provided by the compression rate, i.e. the growth rate of the divergence of the flow. Particularly for the intracluster medium, we find that the dominant contribution against gravity comes from thermal pressure, while compressible turbulence effectively counteracts the support. For this reason it appears to be too simplistic to consider turbulence merely as an effective enhancement of thermal energy.
△ Less
Submitted 19 May, 2017;
originally announced May 2017.
-
Constraining The Single-Degenerate Channel of Type Ia Supernovae With Stable Iron-Group Elements in SNR 3C 397
Authors:
Pranav Dave,
Rahul Kashyap,
Robert Fisher,
Frank Timmes,
Dean Townsley,
Chris Byrohl
Abstract:
Recent Suzaku X-ray spectra of SNR 3C 397 indicate enhanced stable iron-group element abundances of Ni, Mn, Cr, and Fe. Seeking to address key questions about the progenitor and explosion mechanism of 3C 397, we compute nucleosynthetic yields from a suite of multidimensional hydrodynamics models in the near-Chandrasekhar mass, single-degenerate paradigm for supernova Type Ia. Varying the progenito…
▽ More
Recent Suzaku X-ray spectra of SNR 3C 397 indicate enhanced stable iron-group element abundances of Ni, Mn, Cr, and Fe. Seeking to address key questions about the progenitor and explosion mechanism of 3C 397, we compute nucleosynthetic yields from a suite of multidimensional hydrodynamics models in the near-Chandrasekhar mass, single-degenerate paradigm for supernova Type Ia. Varying the progenitor white dwarf internal structure, composition, ignition, and explosion mechanism, we find the best match to the observed iron-peak elements of 3C 397 are dense (central density $\ge$ 6$\times$10$^{9}$ g cm$^{-3}$), low-carbon white dwarfs that undergo a weak, centrally-ignited deflagration, followed by a subsequent detonation. The amount of $^{56}$Ni produced is consistent with a normal or bright normal supernova Type Ia. A pure deflagration of a centrally-ignited, low central density ($\simeq$ 2$\times$10$^{9}$ g cm$^{-3}$) progenitor white dwarf, frequently considered in the literature, is also found to produce good agreement with 3C 397 nucleosynthetic yields, but leads to a subluminous SN Ia event, in conflict with X-ray linewidth data. Additionally, in contrast to prior work which suggested a large super-solar metallicity for the white dwarf progenitor for SNR 3C 397, we find satisfactory agreement for solar and sub-solar metallicity progenitors. We discuss a range of implications our results have for the single-degenerate channel.
△ Less
Submitted 5 May, 2017;
originally announced May 2017.