Showing 1–29 of 29 results for author: Behrens, C

A survey of high-$z$ galaxies: SERRA simulations

Authors: A. Pallottini, A. Ferrara, S. Gallerani, C. Behrens, M. Kohandel, S. Carniani, L. Vallini, S. Salvadori, V. Gelli, L. Sommovigo, V. D'Odorico, F. Di Mascia, E. Pizzati

Abstract: We introduce SERRA, a suite of zoom-in high-resolution ($\sim 10\,\rm pc$) cosmological simulations including non-equilibrium chemistry and on-the-fly radiative transfer. The outputs are post-processed to derive galaxy UV+FIR continuum and emission line properties. Results are compared with available multi-wavelength data to constrain the physical properties (e.g., star formation rates, stellar/ga… ▽ More We introduce SERRA, a suite of zoom-in high-resolution ($\sim 10\,\rm pc$) cosmological simulations including non-equilibrium chemistry and on-the-fly radiative transfer. The outputs are post-processed to derive galaxy UV+FIR continuum and emission line properties. Results are compared with available multi-wavelength data to constrain the physical properties (e.g., star formation rates, stellar/gas/dust mass, metallicity) of high-redshift $6 \lesssim z \lesssim 15$ galaxies. This flagship paper focuses on the $z=7.7$ sub-sample, including 202 galaxies with stellar mass $10^7 M_\odot \lesssim M_\star \lesssim 5\times 10^{10}M_\odot$, and specific star formation ranging from ${\rm sSFR} \sim 100\,{\rm Gyr}^{-1}$ in young, low-mass galaxies to $\sim 10\,{\rm Gyr}^{-1}$ for older, massive ones. At this redshift, SERRA galaxies are typically bursty, i.e. they are located above the Schmidt-Kennicutt relation by a factor $κ_s = 3.03^{+4.9}_{-1.8}$, consistent with recent findings for [OIII] and [CII] emitters at high-$z$. They also show relatively large ${\rm IRX} = L_{\rm FIR}/L_{\rm UV}$ values as a result of their compact/clumpy morphology effectively blocking the stellar UV luminosity. Note that this conclusion might be affected by insufficient spatial resolution at the molecular cloud level. We confirm that early galaxies lie on the standard $\rm [CII]-SFR$ relation; their observed $L_{\rm [OIII]}/L_{\rm [CII]} \simeq 1-10$ ratios can be reproduced by a part of the SERRA galaxies without the need of a top-heavy IMF and/or anomalous C/O abundances. [OI] line intensities are similar to local ones, making ALMA high-$z$ detections challenging but feasible ($\sim 6\,\rm hr$ for a SFR of $50\,M_\odot\,{\rm yr}^{-1}$). △ Less

Submitted 4 May, 2022; v1 submitted 7 January, 2022; originally announced January 2022.

Comments: 23 pages, 14 figures, accepted by MNRAS; part of the data used for this study is available at the website http://cosmology.sns.it/data_access.html

Effects of Spatial Discretization in Lyman-alpha Line Radiation Transfer Simulations

Authors: Peter Camps, Christoph Behrens, Maarten Baes, Anand Utsav Kapoor, Robert Grand

Abstract: We describe the addition of Lyman-alpha resonant line transfer to our dust continuum radiation transfer code SKIRT, verifying our implementation with published results for spherical problems and using some self-designed three-dimensional setups. We specifically test spatial discretization through various grid types, including hierarchical octree grids and unstructured Voronoi tessellations. We the… ▽ More We describe the addition of Lyman-alpha resonant line transfer to our dust continuum radiation transfer code SKIRT, verifying our implementation with published results for spherical problems and using some self-designed three-dimensional setups. We specifically test spatial discretization through various grid types, including hierarchical octree grids and unstructured Voronoi tessellations. We then use a radiation transfer post-processing model for one of the spiral galaxies produced by the Auriga cosmological zoom simulations to investigate the effect of spatial discretization on the synthetic observations. We find that the calculated Lyman-alpha line profiles exhibit an extraordinarily strong dependence on the type and resolution of the spatial grid, rendering the results untrustworthy at best. We attribute this effect to the large gradients in the hydrogen density distribution over small distances, which remain significantly under-resolved in the input model. We therefore argue that further research is needed to determine the required spatial resolution of a hydrodynamical simulation snapshot to enable meaningful Lyman-alpha line transfer post-processing. △ Less

Submitted 1 June, 2021; originally announced June 2021.

Comments: Accepted for publication in The Astrophysical Journal

Infrared emission of $z \sim 6$ galaxies: AGN imprints

Authors: F. Di Mascia, S. Gallerani, C. Behrens, A. Pallottini, S. Carniani, A. Ferrara, P. Barai, F. Vito, T. Zana

Abstract: We investigate the infrared (IR) emission of high-redshift ($z\sim 6$), highly star-forming (${ {\rm SFR} > 100}$ $M_{\rm \odot} {\rm yr}^{-1}$) galaxies, with/without Active Galactic Nuclei (AGN), using a suite of cosmological simulations featuring dust radiative transfer. Synthetic Spectral Energy Distributions (SEDs) are used to quantify the relative contribution of stars/AGN to dust heating. I… ▽ More We investigate the infrared (IR) emission of high-redshift ($z\sim 6$), highly star-forming (${ {\rm SFR} > 100}$ $M_{\rm \odot} {\rm yr}^{-1}$) galaxies, with/without Active Galactic Nuclei (AGN), using a suite of cosmological simulations featuring dust radiative transfer. Synthetic Spectral Energy Distributions (SEDs) are used to quantify the relative contribution of stars/AGN to dust heating. In dusty (${M_{\rm d}\gtrsim 3\times 10^7 M_{\rm \odot}}$) galaxies, $\gtrsim 50-90 \%$ of the UV radiation is obscured by dust inhomogeneities on scales ${\gtrsim 100}$ pc. In runs with AGN, a clumpy, warm ($\approx 250$ K) dust component co-exists with a colder ($\approx 60$ K) and more diffuse one, heated by stars. Warm dust provides up to ${50 \%}$ of the total IR luminosity, but only $\lesssim 0.1 \%$ of the total mass content. The AGN boosts the MIR flux by ${10-100 \times}$ with respect to star forming galaxies, without significantly affecting the FIR. Our simulations successfully reproduce the observed SED of bright (${M_{\rm UV}\sim -26}$) ${z\sim 6}$ quasars, and show that these objects are part of complex, dust-rich merging systems, containing multiple sources (accreting BHs and/or star forming galaxies) in agreement with recent HST and ALMA observations. Our results show that the proposed ORIGINS missions will be able to investigate the MIR properties of dusty star forming galaxies and to obtain good quality spectra of bright quasars at $z\sim 6$. Finally, the MIR-to-FIR flux ratio of faint (${M_{\rm UV}\sim -24}$) AGN is ${>10\times}$ higher than for normal star forming galaxies. This implies that combined JWST/ORIGINS/ALMA observations will be crucial to identify faint and/or dust-obscured AGN in the distant Universe. △ Less

Submitted 17 February, 2021; originally announced February 2021.

Comments: 22 pages, 13 figures, accepted for publication on MNRAS

A Novel Estimator for the Equation of State of the IGM by Ly$α$ Forest Tomography

Authors: Hendrik Müller, Christoph Behrens, David James Edward Marsh

Abstract: We present a novel procedure to estimate the Equation of State of the intergalactic medium in the quasi-linear regime of structure formation based on Ly$α$ forest tomography and apply it to 21 high quality quasar spectra from the UVES\_SQUAD survey at redshift $z=2.5$. Our estimation is based on a full tomographic inversion of the line of sight. We invert the data with two different inversion algo… ▽ More We present a novel procedure to estimate the Equation of State of the intergalactic medium in the quasi-linear regime of structure formation based on Ly$α$ forest tomography and apply it to 21 high quality quasar spectra from the UVES\_SQUAD survey at redshift $z=2.5$. Our estimation is based on a full tomographic inversion of the line of sight. We invert the data with two different inversion algorithms, the iterative Gauss-Newton method and the regularized probability conservation approach, which depend on different priors and compare the inversion results in flux space and in density space. In this way our method combines fitting of absorption profiles in flux space with an analysis of the recovered density distributions featuring prior knowledge of the matter distribution. Our estimates are more precise than existing estimates, in particular on small redshift bins. In particular, we model the temperature-density relation with a power law and observe for the temperature at mean density $T_0 = 13400^{+1700}_{-1300}\,\mathrm{K}$ and for the slope of the power-law (polytropic index) $γ= 1.42 \pm 0.11$ for the power-law parameters describing the temperature-density relation. Moreover, we measure an photoionization rate $Γ_{-12} = 1.1^{+0.16}_{-0.17}$. An implementation of the inversion techniques used will be made publicly available. △ Less

Submitted 6 April, 2021; v1 submitted 12 October, 2020; originally announced October 2020.

Comments: 20 pages, 25 figures, 4 tables, accepted for publication in MNRAS, a python toolkit containing our implementations is publicly available under this https URL https://github.com/hmuellergoe/reglyman

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.

Comments: accepted in MNRAS

Velocity dispersion in the interstellar medium of early galaxies

Authors: M. Kohandel, A. Pallottini, A. Ferrara, S. Carniani, S. Gallerani, L. Vallini, A. Zanella, C. Behrens

Abstract: We study the structure of spatially resolved, line-of-sight velocity dispersion for galaxies in the Epoch of Reionization (EoR) traced by [CII] $158μ\rm{m}$ line emission. Our laboratory is a simulated prototypical Lyman-break galaxy, "Freesia", part of the SERRA suite. The analysis encompasses the redshift range 6 < z < 8, when Freesia is in a very active assembling phase. We build velocity dispe… ▽ More We study the structure of spatially resolved, line-of-sight velocity dispersion for galaxies in the Epoch of Reionization (EoR) traced by [CII] $158μ\rm{m}$ line emission. Our laboratory is a simulated prototypical Lyman-break galaxy, "Freesia", part of the SERRA suite. The analysis encompasses the redshift range 6 < z < 8, when Freesia is in a very active assembling phase. We build velocity dispersion maps for three dynamically distinct evolutionary stages (Spiral Disk at z=7.4, Merger at z=8.0, and Disturbed Disk at z=6.5) using [CII] hyperspectral data cubes. We find that, at a high spatial resolution of 0.005" ($\simeq 30 pc$), the luminosity-weighted average velocity dispersion is $σ_{\rm{CII}}$~23-38 km/s with the highest value belonging to the highly-structured Disturbed Disk stage. Low resolution observations tend to overestimate $σ_{\rm CII}$ values due to beam smearing effects that depend on the specific galaxy structure. For an angular resolution of 0.02" (0.1"), the average velocity dispersion is 16-34% (52-115%) larger than the actual one. The [CII] emitting gas in Freesia has a Toomre parameter $\mathcal{Q}$~0.2 and a rotational-to-dispersion ratio of $v_{\rm c}/σ$~ 7 similar to that observed in z=2-3 galaxies. The primary energy source for the velocity dispersion is due to gravitational processes, such as merging/accretion events; energy input from stellar feedback is generally subdominant (< 10%). Finally, we find that the resolved $σ_{\rm{CII}} - Σ_{\rm SFR}$ relation is relatively flat for $0.02<Σ_{\rm SFR}/{{\rm M}_{\odot}} \mathrm{yr}^{-1} {\mathrm kpc}^{-2} < 30$, with the majority of data lying on the derived analytical relation $σ\propto Σ_{\rm SFR}^{5/7}$. At high SFR, the increased contribution from stellar feedback steepens the relation, and $σ_{\rm{CII}}$ rises slightly. △ Less

Submitted 10 September, 2020; originally announced September 2020.

Comments: Accepted for publication in MNRAS; 16 pages (plus appendix), 8 figures, 3 tables

Early galaxy growth: mergers or gravitational instability?

Authors: A. Zanella, A. Pallottini, A. Ferrara, S. Gallerani, S. Carniani, M. Kohandel, C. Behrens

Abstract: We investigate the spatially-resolved morphology of galaxies in the early Universe. We consider a typical redshift z = 6 Lyman Break galaxy, "Althaea" from the SERRA hydrodynamical simulations. We create mock rest-frame ultraviolet, optical, and far-infrared observations, and perform a two-dimensional morphological analysis to de-blend the galaxy disk from substructures (merging satellites or star… ▽ More We investigate the spatially-resolved morphology of galaxies in the early Universe. We consider a typical redshift z = 6 Lyman Break galaxy, "Althaea" from the SERRA hydrodynamical simulations. We create mock rest-frame ultraviolet, optical, and far-infrared observations, and perform a two-dimensional morphological analysis to de-blend the galaxy disk from substructures (merging satellites or star-forming regions). We find that the [CII]158um emitting region has an effective radius 1.5 - 2.5 times larger than the optical one, consistent with recent observations. This [CII] halo in our simulated galaxy arises as the joint effect of stellar outflows and carbon photoionization by the galaxy UV field, rather than from the emission of unresolved nearby satellites. At the typical angular resolution of current observations (> 0.15") only merging satellites can be detected; detection of star-forming regions requires resolutions of < 0.05". The [CII]-detected satellite has a 2.5 kpc projected distance from the galaxy disk, whereas the star-forming regions are embedded in the disk itself (distance < 1 kpc). This suggests that multi-component systems reported in the literature, which have separations > 2 kpc, are merging satellites, rather than galactic substructures. Finally, the star-forming regions found in our mock maps follow the local L[CII] - SFR_UV relation of galaxy disks, although sampling the low-luminosity, low-SFR tail of the distribution. We show that future JWST observations, bridging UV and [CII] datasets, will be exceptionally suited to characterize galaxy substructures thanks to their exquisite spatial resolution and sensitivity to both low-metallicity and dust-obscured regions that are bright at infrared wavelengths. △ Less

Submitted 8 September, 2020; originally announced September 2020.

Comments: Accepted for publication in MNRAS; 17 pages (plus appendix), 7 figures, 4 tables

AxioNyx: Simulating Mixed Fuzzy and Cold Dark Matter

Authors: Bodo Schwabe, Mateja Gosenca, Christoph Behrens, Jens C. Niemeyer, Richard Easther

Abstract: The distinctive effects of fuzzy dark matter are most visible at non-linear galactic scales. We present the first simulations of mixed fuzzy and cold dark matter, obtained with an extended version of the Nyx code. Fuzzy (or ultralight, or axion-like) dark matter dynamics are governed by the comoving Schrödinger-Poisson equation. This is evolved with a pseudospectral algorithm on the root grid, and… ▽ More The distinctive effects of fuzzy dark matter are most visible at non-linear galactic scales. We present the first simulations of mixed fuzzy and cold dark matter, obtained with an extended version of the Nyx code. Fuzzy (or ultralight, or axion-like) dark matter dynamics are governed by the comoving Schrödinger-Poisson equation. This is evolved with a pseudospectral algorithm on the root grid, and with finite differencing at up to six levels of adaptive refinement. Cold dark matter is evolved with the existing N-body implementation in Nyx. We present the first investigations of spherical collapse in mixed dark matter models, focusing on radial density profiles, velocity spectra and soliton formation in collapsed halos. We find that the effective granule masses decrease in proportion to the fraction of fuzzy dark matter which quadratically suppresses soliton growth, and that a central soliton only forms if the fuzzy dark matter fraction is greater than 10\%. The Nyx framework supports baryonic physics and key astrophysical processes such as star formation. Consequently, AxioNyx will enable increasingly realistic studies of fuzzy dark matter astrophysics. △ Less

Submitted 16 July, 2020; originally announced July 2020.

Journal ref: Phys. Rev. D 102, 083518 (2020)

Axion Miniclusters Made Easy

Authors: David Ellis, David J. E. Marsh, Christoph Behrens

Abstract: We use a modified version of the Peak Patch excursion set formalism to compute the mass and size distribution of QCD axion miniclusters from a fully non-Gaussian initial density field obtained from numerical simulations of axion string decay. We find strong agreement with N-Body simulations at a significantly lower computational cost. We employ a spherical collapse model and provide fitting functi… ▽ More We use a modified version of the Peak Patch excursion set formalism to compute the mass and size distribution of QCD axion miniclusters from a fully non-Gaussian initial density field obtained from numerical simulations of axion string decay. We find strong agreement with N-Body simulations at a significantly lower computational cost. We employ a spherical collapse model and provide fitting functions for the modified barrier in the radiation era. The halo mass function at $z=629$ has a power-law distribution $M^{-0.6}$ for masses within the range $10^{-15}\lesssim M\lesssim 10^{-10}M_{\odot}$, with all masses scaling as $(m_a/50μ\mathrm{eV})^{-0.5}$. We construct merger trees to estimate the collapse redshift and concentration mass relation, $C(M)$, which is well described using analytical results from the initial power spectrum and linear growth. Using the calibrated analytic results to extrapolate to $z=0$, our method predicts a mean concentration $C\sim \mathcal{O}(\text{few})\times10^4$. The low computational cost of our method makes future investigation of the statistics of rare, dense miniclusters easy to achieve. △ Less

Submitted 11 May, 2021; v1 submitted 15 June, 2020; originally announced June 2020.

Journal ref: Phys. Rev. D 103, 083525 (2021)

An Optimized Ly$α$ Forest Inversion Tool Based on a Quantitative Comparison of Existing Reconstruction Methods

Authors: Hendrik Müller, Christoph Behrens, David James Edward Marsh

Abstract: We present a same-level comparison of the most prominent inversion methods for the reconstruction of the matter density field in the quasi-linear regime from the Ly$α$ forest flux. Moreover, we present a pathway for refining the reconstruction in the framework of numerical optimization. We apply this approach to construct a novel hybrid method. The methods which are used so far for matter reconstr… ▽ More We present a same-level comparison of the most prominent inversion methods for the reconstruction of the matter density field in the quasi-linear regime from the Ly$α$ forest flux. Moreover, we present a pathway for refining the reconstruction in the framework of numerical optimization. We apply this approach to construct a novel hybrid method. The methods which are used so far for matter reconstructions are the Richardson-Lucy algorithm, an iterative Gauss-Newton method and a statistical approach assuming a one-to-one correspondence between matter and flux. We study these methods for high spectral resolutions such that thermal broadening becomes relevant. The inversion methods are compared on synthetic data (generated with the lognormal approach) with respect to their performance, accuracy, their stability against noise, and their robustness against systematic uncertainties. We conclude that the iterative Gauss-Newton method offers the most accurate reconstruction, in particular at small S/N, but has also the largest numerical complexity and requires the strongest assumptions. The other two algorithms are faster, comparably precise at small noise-levels, and, in the case of the statistical approach, more robust against inaccurate assumptions on the thermal history of the intergalactic medium (IGM). We use these results to refine the statistical approach using regularization. Our new approach has low numerical complexity and makes few assumptions about the history of the IGM, and is shown to be the most accurate reconstruction at small S/N, even if the thermal history of the IGM is not known. Our code will be made publicly available under https://github.com/hmuellergoe/reglyman. △ Less

Submitted 12 June, 2020; originally announced June 2020.

Comments: 21 pages, 13 figures, 3 tables, submitted to MNRAS, a python toolkit containing our implementations will be made publicly available under https://github.com/hmuellergoe/reglyman

A random-walk model for dark matter halo spins

Authors: Andrew Benson, Christoph Behrens, Yu Lu

Abstract: We extend the random-walk model of Vitvitska et al. for predicting the spins of dark matter halos from their merger histories. Using updated merger rates, orbital parameter distributions, and N-body constraints we show that this model can accurately reproduce the distribution of spin parameters measured in N-body simulations when we include a weak correlation between the spins of halos and the ang… ▽ More We extend the random-walk model of Vitvitska et al. for predicting the spins of dark matter halos from their merger histories. Using updated merger rates, orbital parameter distributions, and N-body constraints we show that this model can accurately reproduce the distribution of spin parameters measured in N-body simulations when we include a weak correlation between the spins of halos and the angular momenta of infalling subhalos. We further show that this model is in approximate agreement with the correlation of the spin magnitude over time as determined from N-body simulations, while it slightly underpredicts the correlation in the direction of the spin vector measured from the same simulations. This model is useful for predicting spins from merger histories derived from non-N-body sources, thereby circumventing the need for very high resolution simulations to permit accurate measurements of spins. It may be particularly relevant to modeling systems which accumulate angular momentum from halos over time (such as galactic discs) - we show that this model makes small but significant changes in the distribution of galactic disc sizes computed using the Galacticus semi-analytic galaxy formation model. △ Less

Submitted 16 June, 2020; v1 submitted 24 January, 2020; originally announced January 2020.

Comments: 11 pages, MNRAS accepted, comments welcome

A physical model for [CII] line emission from galaxies

Authors: A. Ferrara, L. Vallini, A. Pallottini, S. Gallerani, S. Carniani, M. Kohandel, D. Decataldo, C. Behrens

Abstract: A tight relation between the [CII]158$μ$m line luminosity and star formation rate is measured in local galaxies. At high redshift ($z>5$), though, a much larger scatter is observed, with a considerable (15-20\%) fraction of the outliers being [CII]-deficient. Moreover, the [CII] surface brightness ($Σ_{\rm CII}$) of these sources is systematically lower than expected from the local relation. To cl… ▽ More A tight relation between the [CII]158$μ$m line luminosity and star formation rate is measured in local galaxies. At high redshift ($z>5$), though, a much larger scatter is observed, with a considerable (15-20\%) fraction of the outliers being [CII]-deficient. Moreover, the [CII] surface brightness ($Σ_{\rm CII}$) of these sources is systematically lower than expected from the local relation. To clarify the origin of such [CII]-deficiency we have developed an analytical model that fits local [CII] data, and has been validated against radiative transfer simulations performed with CLOUDY. The model predicts an overall increase of $Σ_{\rm CII}$ with the surface star formation rate ($Σ_*$). However, for $Σ_* > 1 M_\odot~{\rm yr}^{-1}~{\rm kpc}^{-2}$, $Σ_{\rm CII}$ saturates. We conclude that underluminous [CII] systems can result from a combination of three factors: (a) large upward deviations from the Kennicutt-Schmidt relation ($κ_s \gg 1$), parameterized by the "burstiness" parameter $κ_s$; (b) low metallicity; (c) low gas density, at least for the most extreme sources (e.g. CR7). Observations of [CII] emission alone cannot break the degeneracy among the above three parameters; this requires additional information coming from other emission lines (e.g. [OIII]88$μ$m, CIII]1909A, CO lines). Simple formulae are given to interpret available data for low and high-$z$ galaxies. △ Less

Submitted 20 August, 2019; originally announced August 2019.

Comments: Accepted for publication in MNRAS. 12 pages, 6 Figures. Comments welcome

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.

Kinematics of $z\geq 6$ galaxies from [CII] line emission

Authors: M. Kohandel, A. Pallottini, A. Ferrara, A. Zanella, C. Behrens, S. Carniani, S. Gallerani, L. Vallini

Abstract: We study the kinematical properties of galaxies in the Epoch of Reionization via the [CII] 158$μ$m line emission. The line profile provides information on the kinematics as well as structural properties such as the presence of a disk and satellites. To understand how these properties are encoded in the line profile, first we develop analytical models from which we identify disk inclination and gas… ▽ More We study the kinematical properties of galaxies in the Epoch of Reionization via the [CII] 158$μ$m line emission. The line profile provides information on the kinematics as well as structural properties such as the presence of a disk and satellites. To understand how these properties are encoded in the line profile, first we develop analytical models from which we identify disk inclination and gas turbulent motions as the key parameters affecting the line profile. To gain further insights, we use "Althaea", a highly-resolved ($30\, \rm pc$) simulated prototypical Lyman Break Galaxy, in the redshift range $z = 6-7$, when the galaxy is in a very active assembling phase. Based on morphology, we select three main dynamical stages: I) Merger , II) Spiral Disk, and III) Disturbed Disk. We identify spectral signatures of merger events, spiral arms, and extra-planar flows in I), II), and III), respectively. We derive a generalised dynamical mass vs. [CII]-line FWHM relation. If precise information on the galaxy inclination is (not) available, the returned mass estimate is accurate within a factor $2$ ($4$). A Tully-Fisher relation is found for the observed high-$z$ galaxies, i.e. $L_{\rm[CII]}\propto (FWHM)^{1.80\pm 0.35}$ for which we provide a simple, physically-based interpretation. Finally, we perform mock ALMA simulations to check the detectability of [CII]. When seen face-on, Althaea is always detected at $> 5σ$; in the edge-on case it remains undetected because the larger intrinsic FWHM pushes the line peak flux below detection limit. This suggests that some of the reported non-detections might be due to inclination effects. △ Less

Submitted 27 May, 2019; originally announced May 2019.

Comments: 14 pages, 12 figures, accepted for publication in MNRAS

Deep into the structure of the first galaxies: SERRA views

Authors: A. Pallottini, A. Ferrara, D Decataldo, S. Gallerani, L. Vallini, S. Carniani, C. Behrens, M. Kohandel, S. Salvadori

Abstract: We study the formation and evolution of a sample of Lyman Break Galaxies in the Epoch of Reionization by using high-resolution ($\sim 10 \,{\rm pc}$), cosmological zoom-in simulations part of the SERRA suite. In SERRA, we follow the interstellar medium (ISM) thermo-chemical non-equilibrium evolution, and perform on-the-fly radiative transfer of the interstellar radiation field (ISRF). The simulati… ▽ More We study the formation and evolution of a sample of Lyman Break Galaxies in the Epoch of Reionization by using high-resolution ($\sim 10 \,{\rm pc}$), cosmological zoom-in simulations part of the SERRA suite. In SERRA, we follow the interstellar medium (ISM) thermo-chemical non-equilibrium evolution, and perform on-the-fly radiative transfer of the interstellar radiation field (ISRF). The simulation outputs are post-processed to compute the emission of far infrared lines ([CII], [NII], and [OIII]). At $z=8$, the most massive galaxy, `Freesia', has an age $t_\star \simeq 409\,{\rm Myr}$, stellar mass $M_{\star} \simeq 4.2\times 10^9 {\rm M}_{\odot}$, and a star formation rate ${\rm SFR} \simeq 11.5\,{\rm M}_{\odot}{\rm yr}^{-1}$, due to a recent burst. Freesia has two stellar components (A and B) separated by $\simeq 2.5\, {\rm kpc}$; other 11 galaxies are found within $56.9 \pm 21.6 \, {\rm kpc}$. The mean ISRF in the Habing band is $G = 7.9\, G_0$ and is spatially uniform; in contrast, the ionisation parameter is $U = 2^{+20}_{-2} \times 10^{-3}$, and has a patchy distribution peaked at the location of star-forming sites. The resulting ionising escape fraction from Freesia is $f_{\rm esc}\simeq 2\%$. While [CII] emission is extended (radius 1.54 kpc), [OIII] is concentrated in Freesia-A (0.85 kpc), where the ratio $Σ_{\rm [OIII]}/Σ_{\rm [CII]} \simeq 10$. As many high-$z$ galaxies, Freesia lies below the local [CII]-SFR relation. We show that this is the general consequence of a starburst phase (pushing the galaxy above the Kennicutt-Schmidt relation) which disrupts/photodissociates the emitting molecular clouds around star-forming sites. Metallicity has a sub-dominant impact on the amplitude of [CII]-SFR deviations. △ Less

Submitted 20 May, 2019; originally announced May 2019.

Comments: 22 pages, 14 figures, accepted by MNRAS

Lyα emission from galaxies in the Epoch of Reionization

Authors: Christoph Behrens, Andrea Pallottini, Andrea Ferrara, Simona Gallerani, Livia Vallini

Abstract: The intrinsic strength of the Ly$α$ line in young, star-forming systems makes it a special tool for studying high-redshift galaxies. However, interpreting observations remains challenging due to the complex radiative transfer involved. Here, we combine state-of-the-art hydrodynamical simulations of 'Althaea', a prototypical Lyman Break Galaxy (LBG, stellar mass $M_{\star}$$\simeq$… ▽ More The intrinsic strength of the Ly$α$ line in young, star-forming systems makes it a special tool for studying high-redshift galaxies. However, interpreting observations remains challenging due to the complex radiative transfer involved. Here, we combine state-of-the-art hydrodynamical simulations of 'Althaea', a prototypical Lyman Break Galaxy (LBG, stellar mass $M_{\star}$$\simeq$ $10^{10}{\rm M}_{\odot})$ at $z=7.2$, with detailed radiative transfer computations of dust/continuum, [CII] 158 $μ$m, and Ly$α$ to clarify the relation between the galaxy properties and its Ly$α$ emission. Althaea exhibits low ($f_α< 1\%$) Ly$α$ escape fractions and Equivalent Widths, EW $\lesssim 6$ Angstrom for the simulated lines of sight, with a large scatter. The correlation between escape fraction and inclination is weak, as a result of the rather chaotic structure of high-redshift galaxies. Low $f_α$ values persist even if we artificially remove neutral gas around star forming regions to mimick the presence of HII regions. The high attenuation is primarily caused by dust clumps co-located with young stellar clusters. We can turn Althaea into a Lyman Alpha Emitter (LAE) only if we artificially remove dust from the clumps, yielding EWs up to $22$ Angstrom. Our study suggests that the LBG-LAE duty-cycle required by recent clustering measurements poses the challenging problem of a dynamically changing dust attenuation. Finally, we find an anti-correlation between the magnitude of Ly$α$-[CII] line velocity shift and Ly$α$ luminosity. △ Less

Submitted 4 June, 2019; v1 submitted 14 March, 2019; originally announced March 2019.

Comments: published in MNRAS

Journal ref: Monthly Notices of the Royal Astronomical Society, Volume 486, Issue 2, p.2197-2209, 2019

First Identification of 10-kpc Scale [CII] 158um Halos around Star-Forming Galaxies at z=5-7

Authors: Seiji Fujimoto, Masami Ouchi, Andrea Ferrara, Andrea Pallottini, R. J. Ivison, Christoph Behrens, Simona Gallerani, Shohei Arata, Hidenobu Yajima, Ken Nagamine

Abstract: We report the discovery of 10-kpc scale [CII] 158um halos surrounding star-forming galaxies in the early Universe. We choose deep ALMA data of 18 galaxies each with a star-formation rate of ~ 10-70 Msun with no signature of AGN whose [CII] lines are individually detected at z=5.153-7.142, and conduct stacking of the [CII] lines and dust-continuum in the uv-visibility plane. The radial profiles of… ▽ More We report the discovery of 10-kpc scale [CII] 158um halos surrounding star-forming galaxies in the early Universe. We choose deep ALMA data of 18 galaxies each with a star-formation rate of ~ 10-70 Msun with no signature of AGN whose [CII] lines are individually detected at z=5.153-7.142, and conduct stacking of the [CII] lines and dust-continuum in the uv-visibility plane. The radial profiles of the surface brightnesses show a 10-kpc scale [CII] halo at the 9.2sigma level significantly extended more than the HST stellar continuum data by a factor of ~5 on the exponential-profile basis, as well as the dust continuum. We also compare the radial profiles of [CII] and Lya halos universally found in star-forming galaxies at this epoch, and find that the scale lengths agree within the 1sigma level. While two independent hydrodynamical zoom-in simulations match the dust and stellar continuum properties, the simulations cannot reproduce the extended [CII] line emission. The existence of the extended [CII] halo is the evidence of outflow remnants in the early galaxies and suggest that the outflows may be dominated by cold-mode outflows expelling the neutral gas. △ Less

Submitted 28 September, 2019; v1 submitted 18 February, 2019; originally announced February 2019.

Comments: 19 pages, 13 figures. ApJ in press

A semi-analytical perspective on massive galaxies at $z\sim0.55$

Authors: D. Stoppacher, F. Prada, A. D. Montero-Dorta, S. Rodríguez-Torres, A. Knebe, G. Favole, W. Cui, A. J. Benson, C. Behrens, A. A. Klypin

Abstract: The most massive and luminous galaxies in the Universe serve as powerful probes to study the formation of structure, the assembly of mass, and cosmology. However, their detailed formation and evolution is still barely understood. Here we extract a sample of massive mock galaxies from the semi-analytical model of galaxy formation (SAM) GALACTICUS from the MultiDark-Galaxies, by replicating the CMAS… ▽ More The most massive and luminous galaxies in the Universe serve as powerful probes to study the formation of structure, the assembly of mass, and cosmology. However, their detailed formation and evolution is still barely understood. Here we extract a sample of massive mock galaxies from the semi-analytical model of galaxy formation (SAM) GALACTICUS from the MultiDark-Galaxies, by replicating the CMASS photometric selection from the SDSS-III Baryon Oscillation Spectroscopic Survey (BOSS). The comparison of the GALACTICUS CMASS-mock with BOSS-CMASS data allows us to explore different aspects of the massive galaxy population at $0.5<z<0.6$, including the galaxy-halo connection and the galaxy clustering. We find good agreement between our modelled galaxies and observations regarding the galaxy-halo connection, but our CMASS-mock over-estimates the clustering amplitude of the 2-point correlation function, due to a smaller number density compared to BOSS, a lack of blue objects, and a small intrinsic scatter in stellar mass at fixed halo mass of $<0.1$ dex. To alleviate this problem, we construct an alternative mock catalogue mimicking the CMASS colour-magnitude distribution by randomly down-sampling the SAM catalogue. This CMASS-mock reproduces the clustering of CMASS galaxies within 1$σ$ and shows some environmental dependency of star formation properties that could be connected to the quenching of star formation and the assembly bias. △ Less

Submitted 14 February, 2019; originally announced February 2019.

Comments: 15 pages, 10 figures, 2 tables, submitted to MNRAS

Dusty galaxies in the Epoch of Reionization: simulations

Authors: C. Behrens, A. Pallottini, A. Ferrara, S. Gallerani, L. Vallini

Abstract: The recent discovery of dusty galaxies well into the Epoch of Reionization (redshift $z>6$) poses challenging questions about the properties of the interstellar medium in these pristine systems. By combining state-of-the-art hydrodynamic and dust radiative transfer simulations, we address these questions focusing on the recently discovered dusty galaxy A2744_YD4 ($z=8.38$, Laporte et al. 2017}). W… ▽ More The recent discovery of dusty galaxies well into the Epoch of Reionization (redshift $z>6$) poses challenging questions about the properties of the interstellar medium in these pristine systems. By combining state-of-the-art hydrodynamic and dust radiative transfer simulations, we address these questions focusing on the recently discovered dusty galaxy A2744_YD4 ($z=8.38$, Laporte et al. 2017}). We show that we can reproduce the observed spectral energy distribution (SED) only using different physical values with respect to the inferred ones by Laporte et al(2017), i.e. a star formation rate of $\mathrm{SFR} = 78\;\rm M_\odot \rm yr^{-1}$, a factor $\approx 4$ higher than deduced from simple Spectral Energy Distribution fitting. In this case we find: (a) dust attenuation (corresponding to $τ_V=1.4$) is consistent with a Milky Way extinction curve; (b) the dust-to-metal ratio is low, $f_\mathrm{d} \sim 0.08$, implying that early dust formation is rather inefficient; (c) the luminosity-weighted dust temperature is high, $T_d=91\pm 23\, \rm K$, as a result of the intense ($\approx 100\times$ MW) interstellar radiation field; (d) due to the high $T_d$, the ALMA Band 7 detection can be explained by a limited dust mass, $M_d=1.6\times 10^6 $M$_\odot$. Finally, the high dust temperatures might solve the puzzling low infrared excess recently deduced for high-$z$ galaxies from the IRX-$β$ relation. △ Less

Submitted 21 February, 2018; originally announced February 2018.

Comments: 15 pages, accepted for publication in MNRAS

MultiDark-Galaxies: data release and first results

Authors: Alexander Knebe, Doris Stoppacher, Francisco Prada, Christoph Behrens, Andrew Benson, Sofia A. Cora, Darren J. Croton, Nelson D. Padilla, Andrés N. Ruiz, Manodeep Sinha, Adam R. H. Stevens, Cristian A. Vega-Martínez, Peter Behroozi, Violeta Gonzalez-Perez, Stefan Gottlöber, Anatoly A. Klypin, Gustavo Yepes, Harry Enke, Noam I. Libeskind, Kristin Riebe, Matthias Steinmetz

Abstract: We present the public release of the MultiDark-Galaxies: three distinct galaxy catalogues derived from one of the Planck cosmology MultiDark simulations (i.e. MDPL2, with a volume of (1 Gpc/$h$)$^{3}$ and mass resolution of $1.5 \times 10^{9} M_{\odot}/h$) by applying the semi-analytic models GALACTICUS, SAG, and SAGE to it. We compare the three models and their conformity with observational data… ▽ More We present the public release of the MultiDark-Galaxies: three distinct galaxy catalogues derived from one of the Planck cosmology MultiDark simulations (i.e. MDPL2, with a volume of (1 Gpc/$h$)$^{3}$ and mass resolution of $1.5 \times 10^{9} M_{\odot}/h$) by applying the semi-analytic models GALACTICUS, SAG, and SAGE to it. We compare the three models and their conformity with observational data for a selection of fundamental properties of galaxies like stellar mass function, star formation rate, cold gas fractions, and metallicities - noting that they sometimes perform differently reflecting model designs and calibrations. We have further selected galaxy subsamples of the catalogues by number densities in stellar mass, cold gas mass, and star formation rate in order to study the clustering statistics of galaxies. We show that despite different treatment of orphan galaxies, i.e. galaxies that lost their dark-matter host halo due to the finite mass resolution of the N-body simulation or tidal stripping, the clustering signal is comparable, and reproduces the observations in all three models - in particular when selecting samples based upon stellar mass. Our catalogues provide a powerful tool to study galaxy formation within a volume comparable to those probed by on-going and future photometric and redshift surveys. All model data consisting of a range of galaxy properties - including broad-band SDSS magnitudes - are publicly available. △ Less

Submitted 23 October, 2017; originally announced October 2017.

Comments: 29 pages, 16 figures, 8 tables, accepted for publication in MNRAS. All data incl. the complete galaxy catalogues for all models are publicly available from the CosmoSim database (http://www.cosmosim.org); a selected set of galaxy properties is available via the Skies & Universes website (http://www.skiesanduniverses.org)

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.

Comments: 11 pages, published in A&A

Journal ref: A&A 614, A31 (2018)

CO line emission from galaxies in the Epoch of Reionization

Authors: L. Vallini, A. Pallottini, A. Ferrara, S. Gallerani, E. Sobacchi, C. Behrens

Abstract: We study the CO line luminosity ($L_{\rm CO}$), the shape of the CO Spectral Line Energy Distribution (SLED), and the value of the CO-to-$\rm H_2$ conversion factor in galaxies in the Epoch of Reionization (EoR). To this aim, we construct a model that simultaneously takes into account the radiative transfer and the clumpy structure of giant molecular clouds (GMCs) where the CO lines are excited. W… ▽ More We study the CO line luminosity ($L_{\rm CO}$), the shape of the CO Spectral Line Energy Distribution (SLED), and the value of the CO-to-$\rm H_2$ conversion factor in galaxies in the Epoch of Reionization (EoR). To this aim, we construct a model that simultaneously takes into account the radiative transfer and the clumpy structure of giant molecular clouds (GMCs) where the CO lines are excited. We then use it to post-process state-of-the-art zoomed, high resolution ($30\, \rm{pc}$), cosmological simulation of a main-sequence ($M_{*}\approx10^{10}\, \rm{M_{\odot}}$, $SFR\approx 100\,\rm{M_{\odot}\, yr^{-1}}$) galaxy, "Althæa", at $z\approx6$. We find that the CO emission traces the inner molecular disk ($r\approx 0.5 \,\rm{kpc}$) of Althæa with the peak of the CO surface brightness co-located with that of the [CII] 158$\rm μm$ emission. Its $L_{\rm CO(1-0)}=10^{4.85}\, \rm{L_{\odot}}$ is comparable to that observed in local galaxies with similar stellar mass. The high ($Σ_{gas} \approx 220\, \rm M_{\odot}\, pc^{-2}$) gas surface density in Althæa, its large Mach number (\mach$\approx 30$), and the warm kinetic temperature ($T_{k}\approx 45 \, \rm K$) of GMCs yield a CO SLED peaked at the CO(7-6) transition, i.e. at relatively high-$J$, and a CO-to-$\rm H_2$ conversion factor $α_{\rm CO}\approx 1.5 \, \rm M_{\odot} \rm (K\, km\, s^{-1}\, pc^2)^{-1} $ lower than that of the Milky Way. The ALMA observing time required to detect (resolve) at 5$σ$ the CO(7-6) line from galaxies similar to Althæa is $\approx13$ h ($\approx 38$ h). △ Less

Submitted 12 September, 2017; originally announced September 2017.

Comments: 16 pages, 14 figures, accepted for publication in MNRAS

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.

Comments: 16 pages, 16 figures, accepted for publication by MNRAS

Core-halo mass relation of ultralight axion dark matter from merger history

Authors: Xiaolong Du, Christoph Behrens, Jens C. Niemeyer, Bodo Schwabe

Abstract: In the context of structure formation with ultralight axion dark matter, we offer an alternative explanation for the mass relation of solitonic cores and their host halos observed in numerical simulations. Our argument is based entirely on the mass gain that occurs during major mergers of binary cores and largely independent of the initial core-halo mass relation assigned to hosts that have just c… ▽ More In the context of structure formation with ultralight axion dark matter, we offer an alternative explanation for the mass relation of solitonic cores and their host halos observed in numerical simulations. Our argument is based entirely on the mass gain that occurs during major mergers of binary cores and largely independent of the initial core-halo mass relation assigned to hosts that have just collapsed. We find a relation between the halo mass $M_h$ and corresponding core mass $M_c$, $M_c\propto M_h^{2β-1}$, where $(1-β)$ is the core mass loss fraction. Following the evolution of core masses in stochastic merger trees, we find empirical evidence for our model. Our results are useful for statistically modeling the effects of dark matter cores on the properties of galaxies and their substructures in axion dark matter cosmologies. △ Less

Submitted 14 February, 2017; v1 submitted 29 September, 2016; originally announced September 2016.

Comments: 6 pages, 5 figures, accepted for publication in PRD

Journal ref: Phys. Rev. D 95, 043519 (2017)

Substructure of fuzzy dark matter haloes

Authors: Xiaolong Du, Christoph Behrens, Jens C. Niemeyer

Abstract: We derive the halo mass function (HMF) for fuzzy dark matter (FDM) by solving the excursion set problem explicitly with a mass-dependent barrier function, which has not been done before. We find that compared to the naive approach of the Sheth--Tormen HMF for FDM, our approach has a higher cut off mass and the cut off mass changes less strongly with redshifts. Using merger trees constructed with a… ▽ More We derive the halo mass function (HMF) for fuzzy dark matter (FDM) by solving the excursion set problem explicitly with a mass-dependent barrier function, which has not been done before. We find that compared to the naive approach of the Sheth--Tormen HMF for FDM, our approach has a higher cut off mass and the cut off mass changes less strongly with redshifts. Using merger trees constructed with a modified version of the Lacey & Cole formalism that accounts for suppressed small scale power and the scale-dependent growth of FDM haloes and the semi-analytic GALACTICUS code, we study the statistics of halo substructure including the effects from dynamical friction and tidal stripping. We find that if the dark matter is a mixture of cold dark matter (CDM) and FDM, there will be a suppression on the halo substructure on small scales which may be able to solve the Missing Satellites Problem faced by the pure CDM model. The suppression becomes stronger with increasing FDM fraction or decreasing FDM mass. Thus, it may be used to constrain the FDM model. △ Less

Submitted 18 December, 2016; v1 submitted 8 August, 2016; originally announced August 2016.

Comments: 12 pages, 18 figures, accepted for publication in MNRAS

Journal ref: Mon. Not. R. Astron. Soc. 465, 941 (2017)

HST Emission Line Galaxies at z ~ 2: Comparing Physical Properties of Lyman Alpha and Optical Emission Line Selected Galaxies

Authors: Alex Hagen, Gregory R. Zeimann, Christoph Behrens, Robin Ciardullo, Henry S. Grasshorn Gebhardt, Caryl Gronwall, Joanna S. Bridge, Derek B. Fox, Donald P. Schneider, Jonathan R. Trump, Guillermo A. Blanc, Yi-Kuan Chiang, Taylor S. Chonis, Steven L. Finkelstein, Gary J. Hill, Shardha Jogee, Eric Gawiser

Abstract: We compare the physical and morphological properties of z ~ 2 Lyman-alpha emitting galaxies (LAEs) identified in the HETDEX Pilot Survey and narrow band studies with those of z ~ 2 optical emission line galaxies (oELGs) identified via HST WFC3 infrared grism spectroscopy. Both sets of galaxies extend over the same range in stellar mass (7.5 < logM < 10.5), size (0.5 < R < 3.0 kpc), and star-format… ▽ More We compare the physical and morphological properties of z ~ 2 Lyman-alpha emitting galaxies (LAEs) identified in the HETDEX Pilot Survey and narrow band studies with those of z ~ 2 optical emission line galaxies (oELGs) identified via HST WFC3 infrared grism spectroscopy. Both sets of galaxies extend over the same range in stellar mass (7.5 < logM < 10.5), size (0.5 < R < 3.0 kpc), and star-formation rate (~1 < SFR < 100). Remarkably, a comparison of the most commonly used physical and morphological parameters -- stellar mass, half-light radius, UV slope, star formation rate, ellipticity, nearest neighbor distance, star formation surface density, specific star formation rate, [O III] luminosity, and [O III] equivalent width -- reveals no statistically significant differences between the populations. This suggests that the processes and conditions which regulate the escape of Ly-alpha from a z ~ 2 star-forming galaxy do not depend on these quantities. In particular, the lack of dependence on the UV slope suggests that Ly-alpha emission is not being significantly modulated by diffuse dust in the interstellar medium. We develop a simple model of Ly-alpha emission that connects LAEs to all high-redshift star forming galaxies where the escape of Ly-alpha depends on the sightline through the galaxy. Using this model, we find that mean solid angle for Ly-alpha escape is 2.4+/-0.8 steradians; this value is consistent with those calculated from other studies. △ Less

Submitted 9 December, 2015; originally announced December 2015.

Comments: Accepted to the ApJ. 32 pages, 4 figures, and 2 tables

Inclination Dependence of Lyman-Alpha Properties in a Turbulent Disk Galaxy

Authors: Christoph Behrens, Harald Braun

Abstract: We present simulations of Lyman-Alpha radiation transfer in an isolated disk galaxy with a turbulence sub-grid model, multi-phase interstellar medium and detailed star formation modelling. We investigate the influence of inclination on the observed Lyman-Alpha properties for different snapshots. The Lyman-alpha spectrum, equivalent width distribution and escape fractions vary significantly with th… ▽ More We present simulations of Lyman-Alpha radiation transfer in an isolated disk galaxy with a turbulence sub-grid model, multi-phase interstellar medium and detailed star formation modelling. We investigate the influence of inclination on the observed Lyman-Alpha properties for different snapshots. The Lyman-alpha spectrum, equivalent width distribution and escape fractions vary significantly with the detailed morphology of the disk, leading to variations from one snapshot to another. In particular, we find that supernova-driven cavities near star-forming regions in the simulation can dominate the transmitted Lyman-alpha fraction, suggesting a variability of LAEs on the timescales of the star formation activity. △ Less

Submitted 14 October, 2014; originally announced October 2014.

Comments: accepted for publication, 9 pages

Journal ref: A&A 572, A74 (2014)

Beamed Lyman Alpha Emission through Outflow-Driven Cavities

Authors: Christoph Behrens, Mark Dijkstra, Jens Niemeyer

Abstract: We investigate the radiative transfer of Lyman alpha photons through simplified anisotropic gas distributions, which represent physically motivated extensions of the popular 'shell-models'. Our study is motivated by the notion that (i) shell models do not always reproduce observed Lyman alpha spectral line profiles, (ii) (typical) shell models do not allow for the escape of ionizing photons, and (… ▽ More We investigate the radiative transfer of Lyman alpha photons through simplified anisotropic gas distributions, which represent physically motivated extensions of the popular 'shell-models'. Our study is motivated by the notion that (i) shell models do not always reproduce observed Lyman alpha spectral line profiles, (ii) (typical) shell models do not allow for the escape of ionizing photons, and (iii) the observation & expectation that winds are more complex, anisotropic phenomena. We examine the influence of inclination on the Lyman alpha spectra, relative fluxes and escape fractions. We find the flux to be enhanced/suppressed by factors up to a few depending on the parameter range of the models, corresponding to a boost in equivalent width of the same amplitude if we neglect dust. In general, lower mean optical depths tend to reduce the impact of anisotropies as is expected. We find a correlation between an observed peak in the -- occasionally triple-peaked -- spectrum at the systemic velocity and the existence of a low optical depth cavity along the line of sight. This can be of importance in the search for ionizing photons leaking from high- redshift galaxies since these photons will also be able to escape through the cavity. △ Less

Submitted 22 January, 2014; v1 submitted 20 January, 2014; originally announced January 2014.

Comments: 11 pages, accepted for publication in A&A

Effects of Lyman alpha scattering in the IGM on clustering statistics of Lyman alpha emitters

Authors: Christoph Behrens, Jens Niemeyer

Abstract: We investigate the correlations between the observed fraction of Lyman alpha emission from star-forming galaxies and the large scale structure by post-processing snapshots of a large, high resolution hydrodynamical simulation with a Lyman alpha radiative transfer code at redshifts z=4, 3, 2. We find correlations of the observed fraction with density, density gradient along the line of sight, veloc… ▽ More We investigate the correlations between the observed fraction of Lyman alpha emission from star-forming galaxies and the large scale structure by post-processing snapshots of a large, high resolution hydrodynamical simulation with a Lyman alpha radiative transfer code at redshifts z=4, 3, 2. We find correlations of the observed fraction with density, density gradient along the line of sight, velocity and velocity gradient along the line of sight, all within the same order of magnitude (tens of percent). Additionally, a correlation with the angular momentum of the dark matter halo is detected. In contrast to a previous study, we find no significant deformation of the 2-point correlation function due to selection effects from radiative transfer in the IGM within the limited statistics of the simulation volume. △ Less

Submitted 4 June, 2013; originally announced June 2013.

Comments: 13 pages; accepted for publication in A&A

Journal ref: A&A 556 A5 (2013)