Showing 1–50 of 249 results for author: Wyithe, S

A model for galaxy-galaxy strong lensing statistics in surveys

Authors: Giovanni Ferrami, Stuart Wyithe

Abstract: Photometric wide-area observations in the next decade will be capable of detecting a large number of galaxy-scale strong gravitational lenses, increasing the gravitational lens sample size by orders of magnitude. To aid in forecasting and analysis of these surveys, we construct a flexible model based on observed distributions for the lens and source properties and test it on the results of past le… ▽ More Photometric wide-area observations in the next decade will be capable of detecting a large number of galaxy-scale strong gravitational lenses, increasing the gravitational lens sample size by orders of magnitude. To aid in forecasting and analysis of these surveys, we construct a flexible model based on observed distributions for the lens and source properties and test it on the results of past lens searches, including SL2S, SuGOHI and searches on the COSMOS HST and DES fields. We use this model to estimate the expected yields of some current and planned surveys, including Euclid Wide, Vera Rubin LSST, and Roman High Latitude Wide Area. The model proposed includes a set of free parameters to constrain on the identifiability of a lens in an image, allowing construction of prior probability distributions for different lens detection methods. The code used in this work is made publicly available. △ Less

Submitted 29 June, 2024; v1 submitted 3 April, 2024; originally announced April 2024.

Comments: 18 pages, 15 figures. Accepted in MNRAS

The Impact of ionization Morphology and X-ray Heating on the Cosmological 21cm Skew Spectrum

Authors: J. H. Cook, S. Balu, B. Greig, C. M. Trott, J. L. B. Line, Y. Qin, J. S. B. Wyithe

Abstract: The cosmological 21cm signal offers a potential probe of the early Universe and the first ionizing sources. Current experiments probe the spatially-dependent variance (Gaussianity) of the signal through the power spectrum (PS). The signal however is expected to be highly non-Gaussian due to the complex topology of reionization and X-ray heating. We investigate the non-Gaussianities of X-ray heatin… ▽ More The cosmological 21cm signal offers a potential probe of the early Universe and the first ionizing sources. Current experiments probe the spatially-dependent variance (Gaussianity) of the signal through the power spectrum (PS). The signal however is expected to be highly non-Gaussian due to the complex topology of reionization and X-ray heating. We investigate the non-Gaussianities of X-ray heating and reionization, by calculating the skew spectrum (SS) of the 21cm signal using MERAXES, which couples a semi-analytic galaxy population with semi-numerical reionization simulations. The SS is the cross-spectrum of the quadratic temperature brightness field with itself. We generate a set of seven simulations from $z = 30$ to $z = 5$, varying the halo mass threshold for hosting star-formation, the X-ray luminosity per star-formation rate, and the minimum X-ray energy escaping host galaxies. We find the SS is predominantly negative as a function of redshift, transitioning to positive towards the start of reionization, and peaking during the midpoint of reionization. We do not see a negative dip in the SS during reionization, likely due to the specifics of modelling ionization sources. We normalise the SS by the PS during reionization isolating the non-Gaussianities. We find a trough ($k\sim\,0.1\,\textrm{Mpc}^{-1}$) and peak ($k\sim\,0.4-1\,\textrm{Mpc}^{-1}$) in the normalised SS during the mid to late periods of reionization. These correlate to the ionization topology, and neutral islands in the IGM. We calculate the cosmic variance of the normalised SS, and find these features are detectable in the absence of foregrounds with the SKA_LOW. △ Less

Submitted 27 February, 2024; originally announced February 2024.

Comments: 19 pages, 14 figures

Simulated host galaxy analogs of high-z quasars observed with JWST

Authors: Sabrina Berger, Madeline A. Marshall, J. Stuart B. Wyithe, Tiziana di Matteo, Yueying Ni, Stephen M. Wilkins

Abstract: The hosts of two low-luminosity high-z quasars, J2255+0251 and J2236+0032, were recently detected using JWST's NIRCam instrument. These represent the first high-z quasar host galaxy stellar detections and open a new window into studying high-z quasars. We examine the implications of the measured properties of J2255+0251 and J2236+0032 within the context of the hydrodynamic simulation BlueTides at… ▽ More The hosts of two low-luminosity high-z quasars, J2255+0251 and J2236+0032, were recently detected using JWST's NIRCam instrument. These represent the first high-z quasar host galaxy stellar detections and open a new window into studying high-z quasars. We examine the implications of the measured properties of J2255+0251 and J2236+0032 within the context of the hydrodynamic simulation BlueTides at z = 6.5. We find that these observed quasars fall on the BlueTides stellar to black hole mass relation and have similar luminosities to the brightest simulated quasars. We predict their star formation rates, estimating approximately $10^{2-3}$ $M_{\odot}/ \rm yr$ for both quasar hosts. J2255+0251 and J2236+0032's host galaxy radii also fall within estimates of the radii of the simulated host galaxies of similar luminosity quasars. We generate mock JWST NIRCam images of analogs to the observed quasars within BlueTides and perform a point source removal to illustrate both a qualitative and quantitative comparison of the measured and simulated radii and magnitudes. The quasar subtraction works well for similar luminosity quasars, and the recovered host images are consistent with what was observed for J2255+0251 and J2236+0032, further supporting the success of those observations. We also use our mock imaging pipeline to make predictions for the detection of J2255+0251 and J2236+0032's hosts in upcoming JWST observations. We anticipate that the simulation analogs of future high-z quasar host discoveries will allow us to make accurate predictions of their properties beyond the capabilities of JWST. △ Less

Submitted 18 April, 2024; v1 submitted 16 January, 2024; originally announced January 2024.

Comments: Accepted to MNRAS. 15 pages, 11 figures

Semi-analytic modelling of Pop. III star formation and metallicity evolution -- I. Impact on the UV luminosity functions at z = 9-16

Authors: Emanuele M. Ventura, Yuxiang Qin, Sreedhar Balu, J. Stuart B. Wyithe

Abstract: We implemented Population III (Pop. III) star formation in mini-halos within the MERAXES semi-analytic galaxy formation and reionisation model, run on top of a N-body simulation with $L = 10 h^{-1}$ cMpc with 2048$^3$ particles resolving all dark matter halos down to the mini-halos ($\sim 10^5 M_\odot$). Our modelling includes the chemical evolution of the IGM, with metals released through superno… ▽ More We implemented Population III (Pop. III) star formation in mini-halos within the MERAXES semi-analytic galaxy formation and reionisation model, run on top of a N-body simulation with $L = 10 h^{-1}$ cMpc with 2048$^3$ particles resolving all dark matter halos down to the mini-halos ($\sim 10^5 M_\odot$). Our modelling includes the chemical evolution of the IGM, with metals released through supernova-driven bubbles that expand according to the Sedov-Taylor model. We found that SN-driven metal bubbles are generally small, with radii typically of 150 ckpc at z = 6. Hence, the majority of the first galaxies are likely enriched by their own star formation. However, as reionization progresses, the feedback effects from the UV background become more pronounced, leading to a halt in star formation in low-mass galaxies, after which external chemical enrichment becomes more relevant. We explore the sensitivity of the star formation rate density and stellar mass functions on the unknown values of free parameters. We also discuss the observability of Pop. III dominated systems with JWST, finding that the inclusion of Pop. III galaxies can have a significant effect on the total UV luminosity function at z = 12 - 16. Our results support the idea that the excess of bright galaxies detected with JWST might be explained by the presence of bright top-heavy Pop. III dominated galaxies without requiring an increased star formation efficiency. △ Less

Submitted 21 February, 2024; v1 submitted 14 January, 2024; originally announced January 2024.

Comments: 19 Pages, 22 Figures, 2 tables, 2 Appendices - Accepted for publication in MNRAS (14th of Feb. 2024) - Meraxes is now publicly available!

Massive Optically Dark Galaxies Unveiled by JWST Challenge Galaxy Formation Models

Authors: Mengyuan Xiao, Pascal Oesch, David Elbaz, Longji Bing, Erica Nelson, Andrea Weibel, Rohan Naidu, Emanuele Daddi, Rychard Bouwens, Jorryt Matthee, Stijn Wuyts, John Chisholm, Gabriel Brammer, Mark Dickinson, Benjamin Magnelli, Lucas Leroy, Pieter van Dokkum, Daniel Schaerer, Thomas Herard-Demanche, Laia Barrufet, Ryan Endsley, Yoshinobu Fudamoto, Carlos Gómez-Guijarro, Rashmi Gottumukkala, Garth Illingworth , et al. (12 additional authors not shown)

Abstract: Over the past decade, the existence of a substantial population of optically invisible, massive galaxies at $z\gtrsim3$ has been implied from mid-infrared to millimeter observations. With the unprecedented sensitivity of the JWST, such extremely massive galaxy candidates have immediately been identified even at $z>7$, in much larger numbers than expected. These discoveries raised a hot debate. If… ▽ More Over the past decade, the existence of a substantial population of optically invisible, massive galaxies at $z\gtrsim3$ has been implied from mid-infrared to millimeter observations. With the unprecedented sensitivity of the JWST, such extremely massive galaxy candidates have immediately been identified even at $z>7$, in much larger numbers than expected. These discoveries raised a hot debate. If confirmed, early, high-mass galaxies challenge the current models of galaxy formation. However, the lack of spectroscopic confirmations leads to uncertain stellar mass ($M_{\star}$) estimates, and the possible presence of active galactic nuclei (AGN) adds further uncertainty. Here, we present the first sample of 36 dust-obscured galaxies with robust spectroscopic redshifts at $z_{\rm spec}=5-9$ from the JWST FRESCO survey. The three most extreme sources at $z\sim5-6$ ($\sim$1 billion years after the Big Bang) are so massive (log$M_{\star}/M_{\odot}$ $\gtrsim11.0$) that they would require, on average, about 50% of the baryons in their halos to be converted into stars -- two to three times higher than even the most efficient galaxies at later times. The extended emission of these galaxies suggests limited contribution by AGN. This population of ultra-massive galaxies accounts for 20% of the total cosmic star formation rate density at $z\sim5-6$, suggesting a substantial proportion of extremely efficient star formation in the early Universe. △ Less

Submitted 5 September, 2023; originally announced September 2023.

Comments: Submitted to Nature. 22 pages, 4 main figures, 7 supplementary figures, 3 supplementary tables. Comments are welcome

Efficient Reionization in a Large Hydrodynamic Galaxy Formation Simulation

Authors: James E. Davies, Simeon Bird, Simon Mutch, Yueying Ni, Yu Feng, Rupert Croft, Tiziana Di Matteo, J. Stuart B. Wyithe

Abstract: Accuracy in the topology and statistics of a simulated Epoch of Reionization (EoR) are vital to draw connections between observations and physical processes. While full radiative transfer models produce the most accurate reionization models, they are highly computationally expensive, and are infeasible for the largest cosmological simulations. Instead, large simulations often include EoR models th… ▽ More Accuracy in the topology and statistics of a simulated Epoch of Reionization (EoR) are vital to draw connections between observations and physical processes. While full radiative transfer models produce the most accurate reionization models, they are highly computationally expensive, and are infeasible for the largest cosmological simulations. Instead, large simulations often include EoR models that are pre-computed via the initial density field, or post-processed where feedback effects are ignored. We introduce Astrid-ES, a resimulation of the Astrid epoch of reionisation $20 > z > 5.5$ which includes an on-the-fly excursion-set reionization algorithm. Astrid-ES produces more accurate reionization histories without significantly impacting the computational time. This model directly utilises the star particles produced in the simulation to calculate the EoR history and includes a UV background which heats the gas particles after their reionization. We contrast the reionization topology and statistics in Astrid-ES with the previously employed parametric reionisation model, finding that in Astrid-ES, ionised regions are more correlated with galaxies, and the 21cm power-spectrum shows an increase in large scale power. We calculate the relation between the size of HII regions and the UV luminosity of the brightest galaxy within them. Prior to the overlap phase, we find a power-law fit of $\mathrm{log} (R) = -0.314 M_\mathrm{UV} - 2.550 \mathrm{log}(1+z) + 7.408$ with a standard deviation $σ_R < 0.15 \mathrm{dex}$ across all mass bins. We also examine the properties of halos throughout reionization, finding that while the properties of halos in the simulation are correlated with the redshift of reionisation, they are not greatly affected by reionisation itself. △ Less

Submitted 8 August, 2023; v1 submitted 13 June, 2023; originally announced June 2023.

Comments: 12 pages, 11 figures. Accepted for publication in Monthly Notices of the Royal Astronomical Society

Implications of $z>{\sim}12$ JWST galaxies for galaxy formation at high redshift

Authors: Yuxiang Qin, Sreedhar Balu, J. Stuart B. Wyithe

Abstract: Using a semi-analytic galaxy-formation model, we study analogues of 8 recently discovered JWST galaxies at $z>{\sim}12$. We select analogues from a cosmological simulation with a $(311{\rm cMpc})^3$ volume and an effective particle number of $10^{12}$ enabling resolution of every atomic-cooling galaxy at $z{\le}20$. We vary model parameters to reproduce the observed UV luminosity function at… ▽ More Using a semi-analytic galaxy-formation model, we study analogues of 8 recently discovered JWST galaxies at $z>{\sim}12$. We select analogues from a cosmological simulation with a $(311{\rm cMpc})^3$ volume and an effective particle number of $10^{12}$ enabling resolution of every atomic-cooling galaxy at $z{\le}20$. We vary model parameters to reproduce the observed UV luminosity function at $5{<}z{<}13$, aiming for a statistically representative high-redshift galaxy mock catalogue. Using the forward-modelled JWST photometry, we identify analogues from this catalogue and study their properties as well as possible evolutionary paths and local environments. We find faint JWST galaxies ($M_{\rm UV}>{\sim}-19.5$) to remain consistent with standard galaxy-formation model and that our fiducial catalogue includes large samples of their analogues. The properties of these analogues broadly agree with conventional SED fitting results, except for having systematically lower redshifts due to the evolving UV luminosity function, and for having higher specific star formation rates as a result of burstier histories in our model. On the other hand, only a handful of bright galaxy analogues can be identified for the observed $z{\sim}12$ galaxies. Moreover, in order to reproduce the $z>{\sim}16$ JWST galaxy candidates, boosted star-forming efficiencies and reduced feedback regulation are necessary relative to models of lower-redshift populations. This suggests star formation in the first galaxies could differ significantly from their lower-redshift counterparts. We also find that these candidates are subject to low-redshift contamination, which is present in our fiducial results as both the dusty or quiescent galaxies at $z{\sim}5$. △ Less

Submitted 5 June, 2023; v1 submitted 29 May, 2023; originally announced May 2023.

Comments: updated with a couple of references and some typos fixed. 18 pages, 11 figures, submitted to MNRAS, comments welcomed

Fisher matrix forecasts on the astrophysics of galaxies during the epoch of reionisation from the 21-cm power spectra

Authors: Sreedhar Balu, Bradley Greig, J. Stuart B. Wyithe

Abstract: The hyperfine 21-cm transition of neutral hydrogen from the early Universe ($z>5$) is a sensitive probe of the formation and evolution of the first luminous sources. Using the Fisher matrix formalism we explore the complex and degenerate high-dimensional parameter space associated with the high-$z$ sources of this era and forecast quantitative constraints from a future 21-cm power spectrum (21-cm… ▽ More The hyperfine 21-cm transition of neutral hydrogen from the early Universe ($z>5$) is a sensitive probe of the formation and evolution of the first luminous sources. Using the Fisher matrix formalism we explore the complex and degenerate high-dimensional parameter space associated with the high-$z$ sources of this era and forecast quantitative constraints from a future 21-cm power spectrum (21-cm PS) detection. This is achieved using MERAXES, a coupled semi-analytic galaxy formation model and reionisation simulation, applied to an $N$-body halo merger tree with a statistically complete population of all atomically cooled galaxies out to $z\sim20$. Our mock observation assumes a 21-cm detection spanning $z \in [5, 24]$ from a 1000 h mock observation with the forthcoming Square Kilometre Array and is calibrated with respect to ultraviolet luminosity functions (UV LFs) at $z\in[5, 10]$, the optical depth of CMB photons to Thompson scattering from Planck, and various constraints on the IGM neutral fraction at $z > 5$. In this work, we focus on the X-ray luminosity, ionising UV photon escape fraction, star formation and supernova feedback of the first galaxies. We demonstrate that it is possible to recover 5 of the 8 parameters describing these properties with better than $50$ per cent precision using just the 21-cm PS. By combining with UV LFs, we are able to improve our forecast, with 5 of the 8 parameters constrained to better than 10 per cent (and all below 50 per cent). △ Less

Submitted 22 August, 2023; v1 submitted 8 May, 2023; originally announced May 2023.

Comments: 14 pages, 5 figures (+3 appendix), Revised and MNRAS accepted version

The JWST FRESCO Survey: Legacy NIRCam/Grism Spectroscopy and Imaging in the two GOODS Fields

Authors: P. A. Oesch, G. Brammer, R. P. Naidu, R. J. Bouwens, J. Chisholm, G. D. Illingworth, J. Matthee, E. Nelson, Y. Qin, N. Reddy, A. Shapley, I. Shivaei, P. van Dokkum, A. Weibel, K. Whitaker, S. Wuyts, A. Covelo-Paz, R. Endsley, Y. Fudamoto, E. Giovinazzo, T. Herard-Demanche, J. Kerutt, I. Kramarenko, I. Labbe, E. Leonova , et al. (16 additional authors not shown)

Abstract: We present the JWST Cycle 1 53.8hr medium program FRESCO, short for "First Reionization Epoch Spectroscopically Complete Observations". FRESCO covers 62 arcmin$^2$ in each of the two GOODS/CANDELS fields for a total area of 124 arcmin$^2$ exploiting JWST's powerful new grism spectroscopic capabilities at near-infrared wavelengths. By obtaining ~2 hr deep NIRCam/grism observations with the F444W fi… ▽ More We present the JWST Cycle 1 53.8hr medium program FRESCO, short for "First Reionization Epoch Spectroscopically Complete Observations". FRESCO covers 62 arcmin$^2$ in each of the two GOODS/CANDELS fields for a total area of 124 arcmin$^2$ exploiting JWST's powerful new grism spectroscopic capabilities at near-infrared wavelengths. By obtaining ~2 hr deep NIRCam/grism observations with the F444W filter, FRESCO yields unprecedented spectra at R~1600 covering 3.8 to 5.0 $μ$m for most galaxies in the NIRCam field-of-view. This setup enables emission line measurements over most of cosmic history, from strong PAH lines at z~0.2-0.5, to Pa$α$ and Pa$β$ at z~1-3, HeI and [SIII] at z~2.5-4.5, H$α$ and [NII] at z~5-6.5, up to [OIII] and H$β$ for z~7-9 galaxies, and possibly even [OII] at z~10-12. FRESCO's grism observations provide total line fluxes for accurately estimating galaxy stellar masses and calibrating slit-loss corrections of NIRSpec/MSA spectra in the same field. Additionally, FRESCO results in a mosaic of F182M, F210M, and F444W imaging in the same fields to a depth of ~28.2 mag (5 $σ$ in 0.32" diameter apertures). Together with this publication, the v1 imaging mosaics are released as high-level science products via MAST. Here, we describe the overall survey design and the key science goals that can be addressed with FRESCO. We also highlight several, early science results, including: spectroscopic redshifts of Lyman break galaxies that were identified almost 20 years ago, the discovery of broad-line active galactic nuclei at z>4, and resolved Pa$α$ maps of galaxies at z~1.4. These results demonstrate the enormous power for serendipitous discovery of NIRCam/grism observations. △ Less

Submitted 16 August, 2023; v1 submitted 4 April, 2023; originally announced April 2023.

Comments: 12 pages, 11 figures; MNRAS in press; for more information on the survey and data releases, see http://jwst-fresco.astro.unige.ch/ and https://archive.stsci.edu/hlsp/fresco

Dark-ages reionization and galaxy formation simulation -- XXI. Constraining the evolution of the ionizing escape fraction

Authors: Simon J. Mutch, Bradley Greig, Yuxiang Qin, Gregory B. Poole, J. Stuart B. Wyithe

Abstract: The fraction of ionizing photons that escape their host galaxies to ionize hydrogen in the inter-galactic medium (IGM) is a critical parameter in analyses of the reionization era. In this paper we use the Meraxes semi-analytic galaxy formation model to infer the mean ionizing photon escape fraction and its dependence on galaxy properties through joint modelling of the observed high redshift galaxy… ▽ More The fraction of ionizing photons that escape their host galaxies to ionize hydrogen in the inter-galactic medium (IGM) is a critical parameter in analyses of the reionization era. In this paper we use the Meraxes semi-analytic galaxy formation model to infer the mean ionizing photon escape fraction and its dependence on galaxy properties through joint modelling of the observed high redshift galaxy population and existing constraints on the reionization history. Using a Bayesian framework, and under the assumption that escape fraction is primarily related to halo mass, we find that the joint constraints of the UV luminosity function, CMB optical depth, and the Ly$α$ forest require an escape fraction of $(18\pm5)\%$ for galaxies within haloes of $M\lesssim10^{9}$M$_\odot$ and $(5\pm2)\%$ for more massive haloes. In terms of galaxy properties, this transition in escape fraction occurs at stellar masses of $M_\star\sim10^7$M$_\odot$, nearly independent of redshift. As a function of redshift, reionization is dominated by the smaller $M_\star\lesssim10^7$M$_\odot$ galaxies with high escape fractions at $z\gtrsim6$ and by the larger $M_\star\gtrsim10^7$M$_\odot$ galaxies with lower escape fractions at $z\lesssim6$. Galaxies with star formation rates of $10^{-2.5}$M$_\odot$yr$^{-1}$ to $10^{-1.5}$M$_\odot$yr$^{-1}$ provide the dominant source of ionizing photons throughout reionization. Our results are consistent with recent direct measurements of a $\sim5\%$ escape fraction from massive galaxies at the end of reionization and support the picture of low mass galaxies being the dominant sources of ionizing photons during reionization. △ Less

Submitted 7 December, 2023; v1 submitted 13 March, 2023; originally announced March 2023.

Comments: Updated to match accepted version for publication in MNRAS. Includes two new figures showing model predictions and related discussion

Gravitational lensing modification of the high redshift galaxy luminosity function

Authors: Giovanni Ferrami, J. Stuart B. Wyithe

Abstract: The bright end of the rest-frame UV luminosity function (UVLF) of high-redshift galaxies is modified by gravitational lensing magnification bias. Motivated by recent discoveries of very high-z galaxies with JWST, we study the dependence of magnification bias on the finite size of sources at $6<z<14$. We calculate the magnification probability distributions and use these to calculate the magnificat… ▽ More The bright end of the rest-frame UV luminosity function (UVLF) of high-redshift galaxies is modified by gravitational lensing magnification bias. Motivated by recent discoveries of very high-z galaxies with JWST, we study the dependence of magnification bias on the finite size of sources at $6<z<14$. We calculate the magnification probability distributions and use these to calculate the magnification bias assuming a rest-frame Schechter UVLF for galaxies at redshift $6<z<14$. We find that the finite size of bright high-redshift galaxies together with lens ellipticity significantly suppresses magnification bias, producing an observed bright end which declines more sharply than the power-law resulting from assumption of point sources. By assuming a luminosity-size relation for the source population and comparing with the observed $z=6$ galaxy luminosity function from Harikane+(2022), we show that the UVLF can be used to set mild constraints on the galaxies intrinsic size, favoring smaller galaxies compared to the fiducial luminosity-size relation. In the future, wide surveys using Euclid and Roman Space Telescope will place stronger constraints. We also tabulate the maximum magnification possible as a function of source size and lens ellipticity. △ Less

Submitted 8 May, 2023; v1 submitted 1 December, 2022; originally announced December 2022.

Comments: 5 pages, 4 figures, 1 table. Accepted for publication in MNRAS Letters

The role of Pop III stars and early black holes in the 21cm signal from Cosmic Dawn

Authors: Emanuele M. Ventura, Alessandro Trinca, Raffaella Schneider, Luca Graziani, Rosa Valiante, J. Stuart B. Wyithe

Abstract: Modeling the 21cm global signal from the Cosmic Dawn is challenging due to the many poorly constrained physical processes that come into play. We address this problem using the semi-analytical code "Cosmic Archaeology Tool" (CAT). CAT follows the evolution of dark matter halos tracking their merger history and provides an ab initio description of their baryonic evolution, starting from the formati… ▽ More Modeling the 21cm global signal from the Cosmic Dawn is challenging due to the many poorly constrained physical processes that come into play. We address this problem using the semi-analytical code "Cosmic Archaeology Tool" (CAT). CAT follows the evolution of dark matter halos tracking their merger history and provides an ab initio description of their baryonic evolution, starting from the formation of the first (Pop III) stars and black holes (BHs) in mini-halos at z > 20. The model is anchored to observations of galaxies and AGN at z < 6 and predicts a reionization history consistent with constraints. In this work we compute the evolution of the mean global 21cm signal between $4\leq z \leq 40$ based on the rate of formation and emission properties of stars and accreting black holes. We obtain an absorption profile with a maximum depth $δ{\rm T_b} = -95$ mK at $z \sim 26.5$ (54 MHz). This feature is quickly suppressed turning into an emission signal at $z = 20$ due to the contribution of accreting BHs that efficiently heat the IGM at $z < 27$. The high-$z$ absorption feature is caused by the early coupling between the spin and kinetic temperature of the IGM induced by Pop III star formation episodes in mini-halos. Once we account for an additional radio background from early BHs, we are able to reproduce the timing and the depth of the EDGES signal only if we consider a smaller X-ray background from accreting BHs, but not the shape. △ Less

Submitted 18 January, 2023; v1 submitted 18 October, 2022; originally announced October 2022.

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

Thermal and Reionisation History within a Large-Volume Semi-Analytic Galaxy Formation Simulation

Authors: Sreedhar Balu, Bradley Greig, Yisheng Qiu, Chris Power, Yuxiang Qin, Simon Mutch, J. Stuart B. Wyithe

Abstract: We predict the 21-cm global signal and power spectra during the Epoch of Reionisation using the MERAXES semi-analytic galaxy formation and reionisation model, updated to include X-ray heating and thermal evolution of the intergalactic medium. Studying the formation and evolution of galaxies together with the reionisation of cosmic hydrogen using semi-analytic models (such as MERAXES) requires N-bo… ▽ More We predict the 21-cm global signal and power spectra during the Epoch of Reionisation using the MERAXES semi-analytic galaxy formation and reionisation model, updated to include X-ray heating and thermal evolution of the intergalactic medium. Studying the formation and evolution of galaxies together with the reionisation of cosmic hydrogen using semi-analytic models (such as MERAXES) requires N-body simulations within large volumes and high mass resolutions. For this, we use a simulation of side-length $210~h^{-1}$ Mpc with $4320^3$ particles resolving dark matter haloes to masses of $5\times10^8~h^{-1}~M_\odot$. To reach the mass resolution of atomically cooled galaxies, thought to be the dominant population contributing to reionisation, at $z=20$ of $\sim 2\times10^7~h^{-1}~M_\odot$, we augment this simulation using the DARKFOREST Monte-Carlo merger tree algorithm (achieving an effective particle count of $\sim10^{12}$). Using this augmented simulation we explore the impact of mass resolution on the predicted reionisation history as well as the impact of X-ray heating on the 21-cm global signal and the 21-cm power spectra. We also explore the cosmic variance of 21-cm statistics within $70^{3}$ $h^{-3}$ Mpc$^3$ sub-volumes. We find that the midpoint of reionisation varies by $Δz\sim0.8$ and that the cosmic variance on the power spectrum is underestimated by a factor of $2-4$ at $k\sim 0.1-0.4$ Mpc$^{-1}$ due to the non-Gaussian nature of the 21-cm signal. To our knowledge, this work represents the first model of both reionisation and galaxy formation which resolves low-mass atomically cooled galaxies while simultaneously sampling sufficiently large scales necessary for exploring the effects of X-rays in the early Universe. △ Less

Submitted 17 February, 2023; v1 submitted 17 October, 2022; originally announced October 2022.

Comments: 17 pages, 13 figures, revised and MNRAS accepted version

JWST's PEARLS: Prime Extragalactic Areas for Reionization and Lensing Science: Project Overview and First Results

Authors: Rogier A. Windhorst, Seth H. Cohen, Rolf A. Jansen, Jake Summers, Scott Tompkins, Christopher J. Conselice, Simon P. Driver, Haojing Yan, Dan Coe, Brenda Frye, Norman Grogin, Anton Koekemoer, Madeline A. Marshall, Rosalia O'Brien, Nor Pirzkal, Aaron Robotham, Russell E. Ryan, Jr., Christopher N. A. Willmer, Timothy Carleton, Jose M. Diego, William C. Keel, Paolo Porto, Caleb Redshaw, Sydney Scheller, Stephen M. Wilkins , et al. (60 additional authors not shown)

Abstract: We give an overview and describe the rationale, methods, and first results from NIRCam images of the JWST "Prime Extragalactic Areas for Reionization and Lensing Science" ("PEARLS") project. PEARLS uses up to eight NIRCam filters to survey several prime extragalactic survey areas: two fields at the North Ecliptic Pole (NEP); seven gravitationally lensing clusters; two high redshift proto-clusters;… ▽ More We give an overview and describe the rationale, methods, and first results from NIRCam images of the JWST "Prime Extragalactic Areas for Reionization and Lensing Science" ("PEARLS") project. PEARLS uses up to eight NIRCam filters to survey several prime extragalactic survey areas: two fields at the North Ecliptic Pole (NEP); seven gravitationally lensing clusters; two high redshift proto-clusters; and the iconic backlit VV 191 galaxy system to map its dust attenuation. PEARLS also includes NIRISS spectra for one of the NEP fields and NIRSpec spectra of two high-redshift quasars. The main goal of PEARLS is to study the epoch of galaxy assembly, AGN growth, and First Light. Five fields, the JWST NEP Time-Domain Field (TDF), IRAC Dark Field (IDF), and three lensing clusters, will be observed in up to four epochs over a year. The cadence and sensitivity of the imaging data are ideally suited to find faint variable objects such as weak AGN, high-redshift supernovae, and cluster caustic transits. Both NEP fields have sightlines through our Galaxy, providing significant numbers of very faint brown dwarfs whose proper motions can be studied. Observations from the first spoke in the NEP TDF are public. This paper presents our first PEARLS observations, their NIRCam data reduction and analysis, our first object catalogs, the 0.9-4.5 $μ$m galaxy counts and Integrated Galaxy Light. We assess the JWST sky brightness in 13 NIRCam filters, yielding our first constraints to diffuse light at 0.9-4.5 μm. PEARLS is designed to be of lasting benefit to the community. △ Less

Submitted 28 November, 2022; v1 submitted 9 September, 2022; originally announced September 2022.

Comments: Accepted to AJ, comments welcome. We ask anyone who uses our public PEARLS (NEP TDF) data to refer to this overview paper

JWST's PEARLS: dust attenuation and gravitational lensing in the backlit-galaxy system VV 191

Authors: William C. Keel, Rogier A. Windhorst, Rolf A. Jansen, Seth H. Cohen, Jake Summers, Benne Holwerda, Sarah T. Bradford, Clayton D. Robertson, Giovanni Ferrami, Stuart Wyithe, Haojing Yan, Christopher J. Conselice, Simon P. Driver, Aaron Robotham, Norman A. Grogin, Christopher N. A. Willmer, Anton M. Koekemoer, Brenda L. Frye, Nimish P. Hathi, Russell E. Ryan, Jr., Nor Pirzkal, Madeline A. Marshall, Dan Coe, Jose M. Diego, Thomas J. Broadhurst , et al. (6 additional authors not shown)

Abstract: We derive the spatial and wavelength behavior of dust attenuation in the multiple-armed spiral galaxy VV191b using backlighting by the superimposed elliptical system VV191a in a pair with an exceptionally favorable geometry for this measurement. Imaging using JWST and HST spans the wavelength range 0.3-4.5 microns with high angular resolution, tracing the dust in detail from 0.6 to 1.5 microns. Di… ▽ More We derive the spatial and wavelength behavior of dust attenuation in the multiple-armed spiral galaxy VV191b using backlighting by the superimposed elliptical system VV191a in a pair with an exceptionally favorable geometry for this measurement. Imaging using JWST and HST spans the wavelength range 0.3-4.5 microns with high angular resolution, tracing the dust in detail from 0.6 to 1.5 microns. Distinct dust lanes continue well beyond the bright spiral arms, and trace a complex web, with a very sharp radial cutoff near 1.7 Petrosian radii. We present attenuation profiles and coverage statistics in each band at radii 14-21 kpc. We derive the attenuation law with wavelength; the data both within and between the dust lanes clearly favor a stronger reddening behavior (R ~ 2.0 between 0.6 and 0.9 microns, approaching unity by 1.5 microns) than found for starbursts and star-forming regions of galaxies. Power-law extinction behavior lambda^(-beta) gives beta=2.1 from 0.6-0.9 microns. R decreases at increasing wavelengths (R~1.1 between 0.9 and 1.5 microns), while beta steepens to 2.5. Mixing regions of different column density flattens the wavelength behavior, so these results suggest a different grain population than in our vicinity. The NIRCam images reveal a lens arc and counterimage from a background galaxy at z~1, spanning 90 degrees azimuthally at 2.8" from the foreground elliptical galaxy nucleus, and an additional weakly-lensed galaxy. The lens model and imaging data give a mass/light ratio 7.6 in solar units within the Einstein radius 2.0 kpc. △ Less

Submitted 21 February, 2023; v1 submitted 30 August, 2022; originally announced August 2022.

Comments: Accepted by Astron. J. Analysis redone since submission, using updated JWST calibrations. Dust reddening behavior is steeper with wavelength and lensed galaxy redshift lower than we first derived

Generating extremely large-volume reionisation simulations

Authors: Bradley Greig, J. Stuart B. Wyithe, Steven G. Murray, Simon J. Mutch, Cathryn M. Trott

Abstract: Preparing for the first detection of the cosmic 21-cm signal from large-scale interferometer experiments requires rigorous testing of the data analysis and reduction pipelines. To validate that these pipelines do not erroneously remove or add features that can mimic the cosmic signal (e.g. from side-lobes or large-scale power leakage), we require reionisation simulations larger than the experiment… ▽ More Preparing for the first detection of the cosmic 21-cm signal from large-scale interferometer experiments requires rigorous testing of the data analysis and reduction pipelines. To validate that these pipelines do not erroneously remove or add features that can mimic the cosmic signal (e.g. from side-lobes or large-scale power leakage), we require reionisation simulations larger than the experiments primary field of view. For an experiment such as the MWA, with a field of view of $\sim25^{2}$ deg.$^{2}$, this would require a simulation of several Gpcs, which is currently infeasible. To overcome this, we developed a simplified version of the semi-numerical reionisation simulation code 21CMFAST preferencing large volumes over some physical accuracy by assuming linear theory for structure formation. With this, we constructed a 7.5 Gpc comoving volume with voxel resolution of $\sim1.17$ cMpc tailored specifically to the binned spectral resolution of the MWA. This simulation was used for validating the pipelines for the 2020 MWA 21-cm power spectrum (PS) upper limits (Trott et al.). We then use this large-volume simulation to explore: (i) whether smaller volume simulations are biased by the missing large-scale modes, (ii) non-Gaussianity in estimates of the cosmic variance, (iii) biases in the recovered 21-cm PS following foreground wedge removal and (iv) the impact of tiling smaller volume simulations to achieve extremely large volumes. In summary, we find: (i) no biases from missing large-scale power, (ii) significant contribution from non-Gaussianity in the cosmic variance as expected following Mondal et al. (iii) an over-estimate of the 21-cm PS of 10-20 per cent following wedge mode excision for our particular model and (iv) tiling smaller volume simulations under-estimates the large-scale power and also the estimated cosmic variance. △ Less

Submitted 20 May, 2022; originally announced May 2022.

Comments: 13 pages, 9 figures and 1 table. Submitted to MNRAS, comments welcome

The Prevalence of Galaxy Overdensities Around UV-Luminous Lyman $\mathbfα$ Emitters in the Epoch of Reionization

Authors: E. Leonova, P. A. Oesch, Y. Qin, R. P. Naidu, J. S. B. Wyithe, S. de Barros, R. J. Bouwens, R. S. Ellis, R. M. Endsley, A. Hutter, G. D. Illingworth, J. Kerutt, I. Labbe, N. Laporte, D. Magee, S. J. Mutch, G. W. Roberts-Borsani, R. Smit, D. P. Stark, M. Stefanon, S. Tacchella, A. Zitrin

Abstract: Before the end of the epoch of reionization, the Hydrogen in the Universe was predominantly neutral. This leads to a strong attenuation of Ly$α$ lines of $z\gtrsim6$ galaxies in the intergalactic medium. Nevertheless, Ly$α$ has been detected up to very high redshifts ($z\sim9$) for several especially UV luminous galaxies. Here, we test to what extent the galaxy's local environment might impact the… ▽ More Before the end of the epoch of reionization, the Hydrogen in the Universe was predominantly neutral. This leads to a strong attenuation of Ly$α$ lines of $z\gtrsim6$ galaxies in the intergalactic medium. Nevertheless, Ly$α$ has been detected up to very high redshifts ($z\sim9$) for several especially UV luminous galaxies. Here, we test to what extent the galaxy's local environment might impact the Ly$α$ transmission of such sources. We present an analysis of dedicated Hubble Space Telescope (HST) imaging in the CANDELS/EGS field to search for fainter neighbours around three of the most UV luminous and most distant spectroscopically confirmed Ly$α$ emitters: EGS-zs8-1, EGS-zs8-2 and EGSY-z8p7 at $z_\mathrm{spec}=7.73$, 7.48, and 8.68, respectively. We combine the multi-wavelength HST imaging with Spitzer data to reliably select $z\sim7-9$ galaxies around the central, UV-luminous sources. In all cases, we find a clear enhancement of neighbouring galaxies compared to the expected number in a blank field (by a factor $\sim 3-9\times$). Our analysis thus reveals ubiquitous overdensities around luminous Ly$α$ emitting sources in the heart of the cosmic reionization epoch. We show that our results are in excellent agreement with expectations from the Dragons simulation, confirming the theoretical prediction that the first ionized bubbles preferentially formed in overdense regions. JWST follow-up observations of the neighbouring galaxies identified here will be needed to confirm their physical association and to map out the ionized regions produced by these sources. △ Less

Submitted 14 December, 2021; originally announced December 2021.

Comments: 13 pages, 14 figures, Submitted to MNRAS

Epoch of Reionization Power Spectrum Limits from Murchison Widefield Array Data Targeted at EoR1 Field

Authors: M. Rahimi, B. Pindor, J. L. B. Line, N. Barry, C. M. Trott, R. L. Webster, C. H. Jordan, M. Wilensky, S. Yoshiura, A. Beardsley, J. Bowman, R. Byrne, A. Chokshi, B. J. Hazelton, K. Hasegawa, E. Howard, B. Greig, D. Jacobs, R. Joseph, M. Kolopanis, C. Lynch, B. McKinley, D. A. Mitchell, S. Murray, M. F. Morales , et al. (6 additional authors not shown)

Abstract: Current attempts to measure the 21cm Power Spectrum of neutral hydrogen during the Epoch of Reionization are limited by systematics which produce measured upper limits above both the thermal noise and the expected cosmological signal. These systematics arise from a combination of observational, instrumental, and analysis effects. In order to further understand and mitigate these effects, it is ins… ▽ More Current attempts to measure the 21cm Power Spectrum of neutral hydrogen during the Epoch of Reionization are limited by systematics which produce measured upper limits above both the thermal noise and the expected cosmological signal. These systematics arise from a combination of observational, instrumental, and analysis effects. In order to further understand and mitigate these effects, it is instructive to explore different aspects of existing datasets. One such aspect is the choice of observing field. To date, MWA EoR observations have largely focused on the EoR0 field. In this work, we present a new detailed analysis of the EoR1 field. The EoR1 field is one of the coldest regions of the Southern radio sky, but contains the very bright radio galaxy Fornax-A. The presence of this bright extended source in the primary beam of the interferometer makes the calibration and analysis of EoR1 particularly challenging. We demonstrate the effectiveness of a recently developed shapelet model of Fornax-A in improving the results from this field. We also describe and apply a series of data quality metrics which identify and remove systematically contaminated data. With substantially improved source models, upgraded analysis algorithms and enhanced data quality metrics, we determine EoR power spectrum upper limits based on analysis of the best $\sim$14-hours data observed during 2015 and 2014 at redshifts 6.5, 6.8 and 7.1, with the lowest $2σ$ upper limit at z=6.5 of $Δ^2 \leq (73.78 ~\mathrm{mK)^2}$ at $k=0.13~\mathrm{h~ Mpc^{-1}}$, improving on previous EoR1 measurement results. △ Less

Submitted 7 October, 2021; originally announced October 2021.

Comments: 18 pages, 17 figures, Accepted for publication in MNRAS

Dark-ages Reionization and Galaxy Formation Simulation XX. The Ly$α$ IGM transmission properties and environment of bright galaxies during the Epoch of Reionization

Authors: Yuxiang Qin, J. Stuart B. Wyithe, Pascal A. Oesch, Garth D. Illingworth, Ecaterina Leonova, Simon J. Mutch, Rohan P. Naidu

Abstract: The highly neutral inter-galactic medium (IGM) during the Epoch of Reionization (EoR) is expected to suppress Ly$α$ emission with damping-wing absorption, causing nearly no Ly$α$ detection from star-forming galaxies at $z{\sim}8$. However, spectroscopic observations of the 4 brightest galaxies (${\rm H}_{160}{\sim}25$ mag) at these redshifts do reveal prominent Ly$α$ line, suggesting locally ionis… ▽ More The highly neutral inter-galactic medium (IGM) during the Epoch of Reionization (EoR) is expected to suppress Ly$α$ emission with damping-wing absorption, causing nearly no Ly$α$ detection from star-forming galaxies at $z{\sim}8$. However, spectroscopic observations of the 4 brightest galaxies (${\rm H}_{160}{\sim}25$ mag) at these redshifts do reveal prominent Ly$α$ line, suggesting locally ionised IGM. In this paper, we explore the Ly$α$ IGM transmission and environment of bright galaxies during the EoR using the Meraxes semi-analytic model. We find brighter galaxies to be less affected by damping-wing absorption as they are effective at ionizing surrounding neutral hydrogen. Specifically, the brightest sources (${\rm H}_{160}{\lesssim}25.5$ mag) lie in the largest ionized regions in our simulation, and have low attenuation of their Ly$α$ from the IGM (optical depth ${<}1$). Fainter galaxies (25.5 mag${<}{\rm H}_{160}{<}27.5$ mag) have transmission that depends on UV luminosity, leading to a lower incidence of Ly$α$ detection at fainter magnitudes. This luminosity-dependent attenuation explains why Ly$α$ has only been observed in the brightest galaxies at $z{\sim}8$. Follow-up observations have revealed counterparts in the vicinity of these confirmed $z{\sim}8$ Ly$α$ emitters. The environments of our modelled analogues agree with these observations in the number of nearby galaxies, which is a good indicator of whether Ly$α$ can be detected among fainter galaxies. At the current observational limit, galaxies with ${\ge}2$--5 neighbours within $2'{\times}2'$ are ${\sim}2$--3 times more likely to show Ly$α$ emission. JWST will discover an order of magnitude more neighbours, revealing ${\gtrsim}50$ galaxies in the largest ionizing bubbles and facilitating direct study of reionization morphology. △ Less

Submitted 8 August, 2021; originally announced August 2021.

Comments: 10 pages, 3 figures, submitted to MNRAS, comments are welcome

Constraining the 21cm brightness temperature of the IGM at $z$=6.6 around LAEs with the Murchison Widefield Array

Authors: Cathryn M. Trott, C. H. Jordan, J. L. B. Line, C. R. Lynch, S. Yoshiura, B. McKinley, P. Dayal, B. Pindor, A. Hutter, K. Takahashi, R. B. Wayth, N. Barry, A. Beardsley, J. Bowman, R. Byrne, A. Chokshi, B. Greig, K. Hasegawa, B. J. Hazelton, E. Howard, D. Jacobs, M. Kolopanis, D. A. Mitchell, M. F. Morales, S. Murray , et al. (7 additional authors not shown)

Abstract: The locations of Ly-$α$ emitting galaxies (LAEs) at the end of the Epoch of Reionisation (EoR) are expected to correlate with regions of ionised hydrogen, traced by the redshifted 21~cm hyperfine line. Mapping the neutral hydrogen around regions with detected and localised LAEs offers an avenue to constrain the brightness temperature of the Universe within the EoR by providing an expectation for t… ▽ More The locations of Ly-$α$ emitting galaxies (LAEs) at the end of the Epoch of Reionisation (EoR) are expected to correlate with regions of ionised hydrogen, traced by the redshifted 21~cm hyperfine line. Mapping the neutral hydrogen around regions with detected and localised LAEs offers an avenue to constrain the brightness temperature of the Universe within the EoR by providing an expectation for the spatial distribution of the gas, thereby providing prior information unavailable to power spectrum measurements. We use a test set of 12 hours of observations from the Murchison Widefield Array (MWA) in extended array configuration, to constrain the neutral hydrogen signature of 58 LAEs, detected with the Subaru Hypersuprime Cam in the \textit{Silverrush} survey, centred on $z$=6.58. We assume that detectable emitters reside in the centre of ionised HII bubbles during the end of reionization, and predict the redshifted neutral hydrogen signal corresponding to the remaining neutral regions using a set of different ionised bubble radii. A prewhitening matched filter detector is introduced to assess detectability. We demonstrate the ability to detect, or place limits upon, the amplitude of brightness temperature fluctuations, and the characteristic HII bubble size. With our limited data, we constrain the brightness temperature of neutral hydrogen to $Δ{\rm T}_B<$30 mK ($<$200 mK) at 95% (99%) confidence for lognormally-distributed bubbles of radii, $R_B =$ 15$\pm$2$h^{-1}$cMpc. △ Less

Submitted 30 July, 2021; originally announced July 2021.

Comments: 9 pages, 6 figures, accepted for publication in MNRAS

A new MWA limit on the 21 cm Power Spectrum at Redshifts $\sim$ 13 $-$ 17

Authors: S. Yoshiura, B. Pindor, J. L. B. Line, N. Barry, C. M. Trott, A. Beardsley, J. Bowman, R. Byrne, A. Chokshi, B. J. Hazelton, K. Hasegawa, E. Howard, B. Greig, D. Jacobs, C. H. Jordan, R. Joseph, M. Kolopanis, C. Lynch, B. McKinley, D. A. Mitchell, M. F. Morales, S. G. Murray, J. C. Pober, M. Rahimi, K. Takahashi , et al. (7 additional authors not shown)

Abstract: Observations in the lowest MWA band between $75-100$ MHz have the potential to constrain the distribution of neutral hydrogen in the intergalactic medium at redshift $\sim 13-17$. Using 15 hours of MWA data, we analyse systematics in this band such as radio-frequency interference (RFI), ionospheric and wide field effects. By updating the position of point sources, we mitigate the direction indepen… ▽ More Observations in the lowest MWA band between $75-100$ MHz have the potential to constrain the distribution of neutral hydrogen in the intergalactic medium at redshift $\sim 13-17$. Using 15 hours of MWA data, we analyse systematics in this band such as radio-frequency interference (RFI), ionospheric and wide field effects. By updating the position of point sources, we mitigate the direction independent calibration error due to ionospheric offsets. Our calibration strategy is optimized for the lowest frequency bands by reducing the number of direction dependent calibrators and taking into account radio sources within a wider field of view. We remove data polluted by systematics based on the RFI occupancy and ionospheric conditions, finally selecting 5.5 hours of the cleanest data. Using these data, we obtain two sigma upper limits on the 21 cm power spectrum in the range of $0.1\lessapprox k \lessapprox 1 ~\rm ~h~Mpc^{-1}$ and at $z$=14.2, 15.2 and 16.5, with the lowest limit being $6.3\times 10^6 ~\rm mK^2$ at $\rm k=0.14 \rm ~h~Mpc^{-1}$ and at $z=15.2$ with a possibility of a few \% of signal loss due to direction independent calibration. △ Less

Submitted 26 May, 2021; originally announced May 2021.

Comments: 17 pages, 18 figures, accepted to MNRAS

Observing the host galaxies of high-redshift quasars with JWST: predictions from the BlueTides simulation

Authors: Madeline A. Marshall, J. Stuart B. Wyithe, Rogier A. Windhorst, Tiziana Di Matteo, Yueying Ni, Stephen Wilkins, Rupert A. C. Croft, Mira Mechtley

Abstract: The bright emission from high-redshift quasars completely conceals their host galaxies in the rest-frame ultraviolet/optical, with detection of the hosts in these wavelengths eluding even the Hubble Space Telescope (HST) using detailed point spread function (PSF) modelling techniques. In this study we produce mock images of a sample of z=7 quasars extracted from the BlueTides simulation, and apply… ▽ More The bright emission from high-redshift quasars completely conceals their host galaxies in the rest-frame ultraviolet/optical, with detection of the hosts in these wavelengths eluding even the Hubble Space Telescope (HST) using detailed point spread function (PSF) modelling techniques. In this study we produce mock images of a sample of z=7 quasars extracted from the BlueTides simulation, and apply Markov Chain Monte Carlo-based PSF modelling to determine the detectability of their host galaxies with the James Webb Space Telescope (JWST). While no statistically significant detections are made with HST, we predict that at the same wavelengths and exposure times JWST NIRCam imaging will detect ~50% of quasar host galaxies. We investigate various observational strategies, and find that NIRCam wide-band imaging in the long-wavelength filters results in the highest fraction of successful quasar host detections, detecting >80% of the hosts of bright quasars in exposure times of 5 ks. Exposure times of ~5 ks are required to detect the majority of host galaxies in the NIRCam wide-band filters, however even 10 ks exposures with MIRI result in <30% successful host detections. We find no significant trends between galaxy properties and their detectability. The PSF modelling can accurately recover the host magnitudes, radii, and spatial distribution of the larger-scale emission, when accounting for the central core being contaminated by residual quasar flux. Care should be made when interpreting the host properties measured using PSF modelling. △ Less

Submitted 9 June, 2021; v1 submitted 4 January, 2021; originally announced January 2021.

Comments: 21 pages, 14 figures. Resubmitted to MNRAS following referee feedback, with most notable improvements in Section 5

Exploring reionisation and high-z galaxy observables with recent multi-redshift MWA upper limits on the 21-cm signal

Authors: Bradley Greig, Cathryn M. Trott, Nichole Barry, Simon J. Mutch, Bart Pindor, Rachel L. Webster, J. Stuart B. Wyithe

Abstract: We use the latest multi-redshift ($z=6.5-8.7$) upper limits on the 21-cm signal from the Murchison Widefield Array (MWA) to explore astrophysical models which are inconsistent with the data. These upper limits are achieved using 298 h of carefully excised data over four observing seasons. To explore these upper limits in the context of reionisation astrophysics, we use 21CMMC. We then connect the… ▽ More We use the latest multi-redshift ($z=6.5-8.7$) upper limits on the 21-cm signal from the Murchison Widefield Array (MWA) to explore astrophysical models which are inconsistent with the data. These upper limits are achieved using 298 h of carefully excised data over four observing seasons. To explore these upper limits in the context of reionisation astrophysics, we use 21CMMC. We then connect the disfavoured regions of parameter space to existing observational constraints on reionisation such as high-$z$ galaxy ultra-violet (UV) luminosity functions, background UV photoionisation rate, intergalactic medium (IGM) neutral fraction, the electron scattering optical depth and the soft-band X-ray emissivity. We find the vast majority of models disfavoured by the MWA limits are already inconsistent with existing observational constraints. These inconsistent models arise from two classes of models: (i) `cold' reionisation and (ii) pure matter density fluctuations (i.e. no reionisation). However, a small subsample of models are consistent implying the existing MWA limits provide unique information in disfavouring models of reionisation, albeit extremely weakly. We also provide the first limits on the soft-band X-ray emissivity from galaxies at high redshifts, finding $1σ$ lower limits of $ε_{{\rm X},0.5-2~{\rm keV}}\gtrsim10^{34.5}$ erg s$^{-1}$ Mpc$^{-3}$. Finally, we recover 95 per cent disfavoured limits on the IGM spin temperature of $\bar{T}_{\rm S}\lesssim$ 1.3, 1.4, 1.5, 1.8, 2.1, 2.4 K at $z=6.5, 6.8, 7.1, 7.8, 8.2, 8.7$. With this we infer the IGM must have undergone, at the very least, a small amount of X-ray heating. Note, the limits on $ε_{{\rm X},0.5-2~{\rm keV}}$ and $\bar{T}_{\rm S}$ are conditional on the IGM neutral fraction. △ Less

Submitted 5 August, 2020; originally announced August 2020.

Comments: 14 pages, 7 figures, 1 table. Submitted to MNRAS, comments welcome. arXiv admin note: text overlap with arXiv:2006.03203

An efficient hybrid method to produce high resolution large volume dark matter simulations for semi-analytic models of reionisation

Authors: Yisheng Qiu, Simon J. Mutch, Pascal J. Elahi, Rhys J. J. Poulton, Chris Power, J. Stuart B. Wyithe

Abstract: Resolving faint galaxies in large volumes is critical for accurate cosmic reionisation simulations. While less demanding than hydrodynamical simulations, semi-analytic reionisation models still require very large N-body simulations in order to resolve the atomic cooling limit across the whole reionisation history within box sizes $\gtrsim 100 \, h^{-1} {\rm Mpc}$. To facilitate this, we extend the… ▽ More Resolving faint galaxies in large volumes is critical for accurate cosmic reionisation simulations. While less demanding than hydrodynamical simulations, semi-analytic reionisation models still require very large N-body simulations in order to resolve the atomic cooling limit across the whole reionisation history within box sizes $\gtrsim 100 \, h^{-1} {\rm Mpc}$. To facilitate this, we extend the mass resolution of N-body simulations using a Monte Carlo algorithm. We also propose a method to evolve positions of Monte Carlo halos, which can be an input for semi-analytic reionisation models. To illustrate, we present an extended halo catalogue that reaches a mass resolution of $M_\text{halo} = 3.2 \times 10^7 \, h^{-1} \text{M}_\odot$ in a $105 \, h^{-1} {\rm Mpc}$ box, equivalent to an N-body simulation with $\sim 6800^3$ particles. The resulting halo mass function agrees with smaller volume N-body simulations with higher resolution. Our results also produce consistent two-point correlation functions with analytic halo bias predictions. The extended halo catalogues are applied to the \textsc{meraxes} semi-analytic reionisation model, which improves the predictions on stellar mass functions, star formation rate densities and volume-weighted neutral fractions. Comparison of high resolution large volume simulations with both small volume or low resolution simulations confirms that both low resolution and small volume simulations lead to reionisation ending too rapidly. Lingering discrepancies between the star formation rate functions predicted with and without our extensions can be traced to the uncertain contribution of satellite galaxies. △ Less

Submitted 17 October, 2020; v1 submitted 29 July, 2020; originally announced July 2020.

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

Limits to Rest-Frame Ultraviolet Emission From Far-Infrared-Luminous z~6 Quasar Hosts

Authors: Madeline A. Marshall, Mira Mechtley, Rogier A. Windhorst, Seth H. Cohen, Rolf A. Jansen, Linhua Jiang, Victoria R. Jones, J. Stuart B. Wyithe, Xiaohui Fan, Nimish P. Hathi, Knud Jahnke, William C. Keel, Anton M. Koekemoer, Victor Marian, Keven Ren, Jenna Robinson, Huub J. A. Röttgering, Russell E. Ryan Jr., Evan Scannapieco, Donald P. Schneider, Glenn Schneider, Brent M. Smith, Haojing Yan

Abstract: We report on a Hubble Space Telescope search for rest-frame ultraviolet emission from the host galaxies of five far-infrared-luminous $z\simeq{}6$ quasars and the $z=5.85$ hot-dust free quasar SDSS J0005-0006. We perform 2D surface brightness modeling for each quasar using a Markov-Chain Monte-Carlo estimator, to simultaneously fit and subtract the quasar point source in order to constrain the und… ▽ More We report on a Hubble Space Telescope search for rest-frame ultraviolet emission from the host galaxies of five far-infrared-luminous $z\simeq{}6$ quasars and the $z=5.85$ hot-dust free quasar SDSS J0005-0006. We perform 2D surface brightness modeling for each quasar using a Markov-Chain Monte-Carlo estimator, to simultaneously fit and subtract the quasar point source in order to constrain the underlying host galaxy emission. We measure upper limits for the quasar host galaxies of $m_J>22.7$ mag and $m_H>22.4$ mag, corresponding to stellar masses of $M_\ast<2\times10^{11}M_\odot$. These stellar mass limits are consistent with the local $M_{\textrm{BH}}$-$M_\ast$ relation. Our flux limits are consistent with those predicted for the UV stellar populations of $z\simeq6$ host galaxies, but likely in the presence of significant dust ($\langle A_{\mathrm{UV}}\rangle\simeq 2.6$ mag). We also detect a total of up to 9 potential $z\simeq6$ quasar companion galaxies surrounding five of the six quasars, separated from the quasars by 1.4''-3.2'', or 8.4-19.4 kpc, which may be interacting with the quasar hosts. These nearby companion galaxies have UV absolute magnitudes of -22.1 to -19.9 mag, and UV spectral slopes $β$ of -2.0 to -0.2, consistent with luminous star-forming galaxies at $z\simeq6$. These results suggest that the quasars are in dense environments typical of luminous $z\simeq6$ galaxies. However, we cannot rule out the possibility that some of these companions are foreground interlopers. Infrared observations with the James Webb Space Telescope will be needed to detect the $z\simeq6$ quasar host galaxies and better constrain their stellar mass and dust content. △ Less

Submitted 27 July, 2020; originally announced July 2020.

Comments: 22 pages, 13 figures. Accepted for publication in ApJ

Stacking Redshifted 21cm Images of HII Regions Around High Redshift Galaxies as a Probe of Early Reionization

Authors: James E Davies, Rupert A. Croft, Tiziana Di-Matteo, Bradley Greig, Yu Feng, Stuart Wyithe

Abstract: A number of current and future experiments aim to detect the reionization of neutral hydrogen by the first stars and galaxies in the Universe via the redshifted 21cm line. Using the \textsc{BlueTides} simulation, we investigate the measurement of an \textit{average} ionised region towards the beginning of reionization by stacking redshifted 21cm images around optically identified bright galaxies u… ▽ More A number of current and future experiments aim to detect the reionization of neutral hydrogen by the first stars and galaxies in the Universe via the redshifted 21cm line. Using the \textsc{BlueTides} simulation, we investigate the measurement of an \textit{average} ionised region towards the beginning of reionization by stacking redshifted 21cm images around optically identified bright galaxies using mock observations. We find that with an SKA 1000 hour observation, assuming perfect foreground subtraction, a $5σ$ detection of a stacked HII region can be made with 30 images around some of the brightest galaxies in \textsc{bluetides} (brighter than $M_{UV} < -22.75$) at $z=9$ (corresponding to a neutral fraction of 90.1 \% in our model). We present simulated relationships between the UV magnitude of galaxies, the sizes of the ionised regions they reside in, and the shape of the stacked profiles. These mock observations can also distinguish between scenarios where the IGM is in net emission or absorption of 21cm photons. Once 21cm foreground contamination is included, we find that even with up to 200 images around these rare, bright galaxies, only a tentative $> 1σ$ detection will be possible. However, partial foreground subtraction substantially improves signal-to-noise. For example, we predict that reducing the area of Fourier space dominated by foregrounds by 50 (80) percent will allow $> 3σ$ ($> 5σ$) detections of ionised regions at $z=9$. △ Less

Submitted 19 November, 2020; v1 submitted 27 July, 2020; originally announced July 2020.

Comments: 12 pages, 17 figures, Accepted for publication in MNRAS

A Strong-Lensing Model for the WMDF JWST/GTO Very Rich Cluster Abell 1489

Authors: Adi Zitrin, Ana Acebron, Dan Coe, Patrick L. Kelly, Anton M. Koekemoer, Mario Nonino, Rogier A. Windhorst, Brenda Frye, Massimo Pascale, Tom Broadhurst, Seth H. Cohen, Jose M. Diego, Steven L. Finkelstein, Rolf A. Jansen, Rebecca L. Larson, Haojing Yan, Mehmet Alpaslan, Christopher J. Conselice, Alex Griffiths, Louis-Gregory Strolger, J. Stuart B. Wyithe

Abstract: We present a first strong-lensing model for the galaxy cluster RM J121218.5+273255.1 ($z=0.35$; hereafter RMJ1212; also known as Abell 1489). This cluster is amongst the top 0.1\% richest clusters in the redMaPPer catalog; it is significantly detected in X-ray and through the Sunyaev-Zel'dovich effect in ROSAT and \emph{Planck} data, respectively; and its optical luminosity distribution implies a… ▽ More We present a first strong-lensing model for the galaxy cluster RM J121218.5+273255.1 ($z=0.35$; hereafter RMJ1212; also known as Abell 1489). This cluster is amongst the top 0.1\% richest clusters in the redMaPPer catalog; it is significantly detected in X-ray and through the Sunyaev-Zel'dovich effect in ROSAT and \emph{Planck} data, respectively; and its optical luminosity distribution implies a very large lens, following mass-to-light scaling relations. Based on these properties it was chosen for the Webb Medium Deep Fields (WMDF) JWST/GTO program. In preparation for this program, RMJ1212 was recently imaged with GMOS on Gemini North and in seven optical and near-infrared bands with the \emph{Hubble Space Telescope}. We use these data to map the inner mass distribution of the cluster, uncovering various sets of multiple images. We also search for high-redshift candidates in the data, as well as for transient sources. We find over a dozen high-redshift ($z\gtrsim6$) candidates based on both photometric redshift and the dropout technique. No prominent ($\gtrsim5 σ$) transients were found in the data between the two HST visits. Our lensing analysis reveals a relatively large lens with an effective Einstein radius of $θ_{E}\simeq32\pm3''$ ($z_{s}=2$), in broad agreement with the scaling-relation expectations. RMJ1212 demonstrates that powerful lensing clusters can be selected in a robust and automated way following the light-traces-mass assumption. △ Less

Submitted 22 July, 2020; originally announced July 2020.

Comments: 17 pages, 6 figures, 2 tables; To be submitted

Deep multi-redshift limits on Epoch of Reionisation 21cm Power Spectra from Four Seasons of Murchison Widefield Array Observations

Authors: Cathryn M. Trott, C. H. Jordan, S. Midgley, N. Barry, B. Greig, B. Pindor, J. H. Cook, G. Sleap, S. J. Tingay, D. Ung, P. Hancock, A. Williams, J. Bowman, R. Byrne, A. Chokshi, B. J. Hazelton, K. Hasegawa, D. Jacobs, R. C. Joseph, W. Li, J. L. B Line, C. Lynch, B. McKinley, D. A. Mitchell, M. F. Morales , et al. (11 additional authors not shown)

Abstract: We compute the spherically-averaged power spectrum from four seasons of data obtained for the Epoch of Reionisation (EoR) project observed with the Murchison Widefield Array (MWA). We measure the EoR power spectrum over $k= 0.07-3.0~h$Mpc$^{-1}$ at redshifts $z=6.5-8.7$. The largest aggregation of 110 hours on EoR0 high-band (3,340 observations), yields a lowest measurement of (43~mK)$^2$ = 1.8… ▽ More We compute the spherically-averaged power spectrum from four seasons of data obtained for the Epoch of Reionisation (EoR) project observed with the Murchison Widefield Array (MWA). We measure the EoR power spectrum over $k= 0.07-3.0~h$Mpc$^{-1}$ at redshifts $z=6.5-8.7$. The largest aggregation of 110 hours on EoR0 high-band (3,340 observations), yields a lowest measurement of (43~mK)$^2$ = 1.8$\times$10$^3$ mK$^2$ at $k$=0.14~$h$Mpc$^{-1}$ and $z=6.5$ (2$σ$ thermal noise plus sample variance). Using the Real-Time System to calibrate and the CHIPS pipeline to estimate power spectra, we select the best observations from the central five pointings within the 2013--2016 observing seasons, observing three independent fields and in two frequency bands. This yields 13,591 2-minute snapshots (453 hours), based on a quality assurance metric that measures ionospheric activity. We perform another cut to remove poorly-calibrated data, based on power in the foreground-dominated and EoR-dominated regions of the two-dimensional power spectrum, reducing the set to 12,569 observations (419 hours). These data are processed in groups of 20 observations, to retain the capacity to identify poor data, and used to analyse the evolution and structure of the data over field, frequency, and data quality. We subsequently choose the cleanest 8,935 observations (298 hours of data) to form integrated power spectra over the different fields, pointings and redshift ranges. △ Less

Submitted 6 February, 2020; originally announced February 2020.

Comments: 19 pages, 29 figures, accepted for publication in Monthly Notices of the Royal Astronomical Society

The host galaxies of z=7 quasars: predictions from the BlueTides simulation

Authors: Madeline A. Marshall, Yueying Ni, Tiziana Di Matteo, J. Stuart B. Wyithe, Stephen Wilkins, Rupert A. C. Croft, Jussi K. Kuusisto

Abstract: We examine the properties of the host galaxies of $z=7$ quasars using the large volume, cosmological hydrodynamical simulation BlueTides. We find that the 10 most massive black holes and the 191 quasars in the simulation (with $M_{\textrm{UV,AGN}}<M_{\textrm{UV,host}}$) are hosted by massive galaxies with stellar masses $\log(M_\ast/M_\odot)=10.8\pm0.2$, and $10.2\pm0.4$, which have large star for… ▽ More We examine the properties of the host galaxies of $z=7$ quasars using the large volume, cosmological hydrodynamical simulation BlueTides. We find that the 10 most massive black holes and the 191 quasars in the simulation (with $M_{\textrm{UV,AGN}}<M_{\textrm{UV,host}}$) are hosted by massive galaxies with stellar masses $\log(M_\ast/M_\odot)=10.8\pm0.2$, and $10.2\pm0.4$, which have large star formation rates, of $513\substack{+1225 \\ -351}M_\odot/\rm{yr}$ and $191\substack{+288 \\ -120}M_\odot/\rm{yr}$, respectively. The hosts of the most massive black holes and quasars in BlueTides are generally bulge-dominated, with bulge-to-total mass ratio $B/T\simeq0.85\pm0.1$, however their morphologies are not biased relative to the overall $z=7$ galaxy sample. We find that the hosts of the most massive black holes and quasars are significantly more compact, with half-mass radii $R_{0.5}=0.41\substack{+0.18 \\ -0.14}$ kpc and $0.40\substack{+0.11 \\ -0.09}$ kpc respectively; galaxies with similar masses and luminosities have a wider range of sizes with a larger median value, $R_{0.5}=0.71\substack{+0.28 \\ -0.25}$ kpc. We make mock James Webb Space Telescope (JWST) images of these quasars and their host galaxies. We find that distinguishing the host from the quasar emission will be possible but still challenging with JWST, due to the small sizes of quasar hosts. We find that quasar samples are biased tracers of the intrinsic black hole--stellar mass relation, following a relation that is 0.2 dex higher than that of the full galaxy sample. Finally, we find that the most massive black holes and quasars are more likely to be found in denser environments than the typical $M_{\textrm{BH}}>10^{6.5}M_\odot$ black hole, indicating that minor mergers play at least some role in growing black holes in the early Universe. △ Less

Submitted 25 September, 2020; v1 submitted 6 December, 2019; originally announced December 2019.

Comments: 19 pages, 20 figures. Accepted for publication in MNRAS

First Season MWA Phase II EoR Power Spectrum Results at Redshift 7

Authors: W. Li, J. C. Pober, N. Barry, B. J. Hazelton, M. F. Morales, C. M. Trott, A. Lanman, M. Wilensky, I. Sullivan, A. P. Beardsley, T. Booler, J. D. Bowman, R. Byrne, B. Crosse, D. Emrich, T. M. O. Franzen, K. Hasegawa, L. Horsley, M. Johnston-Hollitt, D. C. Jacobs, C. H. Jordan, R. C. Joseph, T. Kaneuji, D. L. Kaplan, D. Kenney , et al. (22 additional authors not shown)

Abstract: The compact configuration of Phase II of the Murchison Widefield Array (MWA) consists of both a redundant subarray and pseudo-random baselines, offering unique opportunities to perform sky-model and redundant interferometric calibration. The highly redundant hexagonal cores give improved power spectrum sensitivity. In this paper, we present the analysis of nearly 40 hours of data targeting one of… ▽ More The compact configuration of Phase II of the Murchison Widefield Array (MWA) consists of both a redundant subarray and pseudo-random baselines, offering unique opportunities to perform sky-model and redundant interferometric calibration. The highly redundant hexagonal cores give improved power spectrum sensitivity. In this paper, we present the analysis of nearly 40 hours of data targeting one of the MWA's EoR fields observed in 2016. We use both improved analysis techniques presented in Barry et al. (2019) as well as several additional techniques developed for this work, including data quality control methods and interferometric calibration approaches. We show the EoR power spectrum limits at redshift 6.5, 6.8 and 7.1 based on our deep analysis on this 40-hour data set. These limits span a range in $k$ space of $0.18$ $h$ $\mathrm{Mpc^{-1}}$ $<k<1.6$ $h$ $\mathrm{Mpc^{-1}}$, with a lowest measurement of $Δ^2\leqslant2.39\times 10^3$ $\mathrm{mK}^2$ at $k=0.59$ $h$ $\mathrm{Mpc^{-1}}$ and $z=6.5$. △ Less

Submitted 20 December, 2019; v1 submitted 22 November, 2019; originally announced November 2019.

Comments: 26 pages, 15 figures, 3 tables. Updated acknowledgements. Published in ApJ

Journal ref: ApJ (2019) 887 141

Dark-ages reionization and galaxy formation simulation -- XVIII. The high-redshift evolution of black holes and their host galaxies

Authors: Madeline A. Marshall, Simon J. Mutch, Yuxiang Qin, Gregory B. Poole, J. Stuart B. Wyithe

Abstract: Correlations between black holes and their host galaxies provide insight into what drives black hole-host co-evolution. We use the Meraxes semi-analytic model to investigate the growth of black holes and their host galaxies from high redshift to the present day. Our modelling finds no significant evolution in the black hole-bulge and black hole-total stellar mass relations out to a redshift of 8.… ▽ More Correlations between black holes and their host galaxies provide insight into what drives black hole-host co-evolution. We use the Meraxes semi-analytic model to investigate the growth of black holes and their host galaxies from high redshift to the present day. Our modelling finds no significant evolution in the black hole-bulge and black hole-total stellar mass relations out to a redshift of 8. The black hole-total stellar mass relation has similar but slightly larger scatter than the black hole-bulge relation, with the scatter in both decreasing with increasing redshift. In our modelling the growth of galaxies, bulges and black holes are all tightly related, even at the highest redshifts. We find that black hole growth is dominated by instability-driven or secular quasar-mode growth and not by merger-driven growth at all redshifts. Our model also predicts that disc-dominated galaxies lie on the black hole-total stellar mass relation, but lie offset from the black hole-bulge mass relation, in agreement with recent observations and hydrodynamical simulations. △ Less

Submitted 6 April, 2020; v1 submitted 17 October, 2019; originally announced October 2019.

Comments: 14 pages, 11 figures. Accepted for publication in MNRAS

Improving the Epoch of Reionization Power Spectrum Results from Murchison Widefield Array Season 1 Observations

Authors: N. Barry, M. Wilensky, C. M. Trott, B. Pindor, A. P. Beardsley, B. J. Hazelton, I. S. Sullivan, M. F. Morales, J. C. Pober, J. Line, B. Greig, R. Byrne, A. Lanman, W. Li, C. H. Jordan, R. C. Joseph, B. McKinley, M. Rahimi, S. Yoshiura, J. D. Bowman, B. M. Gaensler, J. N. Hewitt, D. C. Jacobs, D. A. Mitchell, N. Udaya Shankar , et al. (5 additional authors not shown)

Abstract: Measurements of 21 cm Epoch of Reionization (EoR) structure are subject to systematics originating from both the analysis and the observation conditions. Using 2013 data from the Murchison Widefield Array (MWA), we show the importance of mitigating both sources of contamination. A direct comparison between results from Beardsley et al. 2016 and our updated analysis demonstrates new precision techn… ▽ More Measurements of 21 cm Epoch of Reionization (EoR) structure are subject to systematics originating from both the analysis and the observation conditions. Using 2013 data from the Murchison Widefield Array (MWA), we show the importance of mitigating both sources of contamination. A direct comparison between results from Beardsley et al. 2016 and our updated analysis demonstrates new precision techniques, lowering analysis systematics by a factor of 2.8 in power. We then further lower systematics by excising observations contaminated by ultra-faint RFI, reducing by an additional factor of 3.8 in power for the zenith pointing. With this enhanced analysis precision and newly developed RFI mitigation, we calculate a noise-dominated upper limit on the EoR structure of $Δ^2 \leq 3.9 \times 10^3$ mK$^2$ at $k=0.20$ $\textit{h}$ Mpc$^{-1}$ and $z=7$ using 21 hr of data, improving previous MWA limits by almost an order of magnitude. △ Less

Submitted 8 October, 2019; v1 submitted 2 September, 2019; originally announced September 2019.

Comments: 21 pages, 10 figures, published in ApJ, updated title and text

Gridded and direct Epoch of Reionisation bispectrum estimates using the Murchison Widefield Array

Authors: Cathryn M. Trott, Catherine A. Watkinson, Christopher H. Jordan, Shintaro Yoshiura, Suman Majumdar, N. Barry, R. Byrne, B. J. Hazelton, K. Hasegawa, R. Joseph, T. Kaneuji, K. Kubota, W. Li, J. Line, C. Lynch, B. McKinley, D. A. Mitchell, M. F. Morales, S. Murray, B. Pindor, J. C. Pober, M. Rahimi, J. Riding, K. Takahashi, S. J. Tingay , et al. (20 additional authors not shown)

Abstract: We apply two methods to estimate the 21~cm bispectrum from data taken within the Epoch of Reionisation (EoR) project of the Murchison Widefield Array (MWA). Using data acquired with the Phase II compact array allows a direct bispectrum estimate to be undertaken on the multiple redundantly-spaced triangles of antenna tiles, as well as an estimate based on data gridded to the $uv$-plane. The direct… ▽ More We apply two methods to estimate the 21~cm bispectrum from data taken within the Epoch of Reionisation (EoR) project of the Murchison Widefield Array (MWA). Using data acquired with the Phase II compact array allows a direct bispectrum estimate to be undertaken on the multiple redundantly-spaced triangles of antenna tiles, as well as an estimate based on data gridded to the $uv$-plane. The direct and gridded bispectrum estimators are applied to 21 hours of high-band (167--197~MHz; $z$=6.2--7.5) data from the 2016 and 2017 observing seasons. Analytic predictions for the bispectrum bias and variance for point source foregrounds are derived. We compare the output of these approaches, the foreground contribution to the signal, and future prospects for measuring the bispectra with redundant and non-redundant arrays. We find that some triangle configurations yield bispectrum estimates that are consistent with the expected noise level after 10 hours, while equilateral configurations are strongly foreground-dominated. Careful choice of triangle configurations may be made to reduce foreground bias that hinders power spectrum estimators, and the 21~cm bispectrum may be accessible in less time than the 21~cm power spectrum for some wave modes, with detections in hundreds of hours. △ Less

Submitted 17 May, 2019; originally announced May 2019.

Comments: 19 pages, 10 figures, accepted for publication in PASA

Dark-ages Reionization and Galaxy Formation Simulation -- XIX: Predictions of infrared excess and cosmic star formation rate density from UV observations

Authors: Yisheng Qiu, Simon J. Mutch, Elisabete da Cunha, Gregory B. Poole, J. Stuart B. Wyithe

Abstract: We present a new analysis of high-redshift UV observations using a semi-analytic galaxy formation model, and provide self-consistent predictions of the infrared excess (IRX) -- $β$ relations and cosmic star formation rate density. We combine the Charlot & Fall dust attenuation model with the Meraxes semi-analytic model, and explore three different parametrisations for the dust optical depths, link… ▽ More We present a new analysis of high-redshift UV observations using a semi-analytic galaxy formation model, and provide self-consistent predictions of the infrared excess (IRX) -- $β$ relations and cosmic star formation rate density. We combine the Charlot & Fall dust attenuation model with the Meraxes semi-analytic model, and explore three different parametrisations for the dust optical depths, linked to star formation rate, dust-to-gas ratio and gas column density respectively. A Bayesian approach is employed to statistically calibrate model free parameters including star formation efficiency, mass loading factor, dust optical depths and reddening slope directly against UV luminosity functions and colour-magnitude relations at z ~ 4-7. The best-fit models show excellent agreement with the observations. We calculate IRX using energy balance arguments, and find that the large intrinsic scatter in the IRX -$β$ plane is driven by the specific star formation rate. Additionally, the difference among the three dust models suggests a factor of two systematic uncertainty in the dust-corrected star formation rate when using the Meurer IRX - $β$ relation at z > 4. △ Less

Submitted 7 August, 2019; v1 submitted 7 May, 2019; originally announced May 2019.

Comments: 15 pages, 10 figure. Accepted for publication in MNRAS

Robust statistics toward detection of the 21 cm signal from the Epoch of Reionisation

Authors: Cathryn M. Trott, Shih Ching Fu, Steven Murray, Christopher Jordan, Jack Line, N. Barry, R. Byrne, B. J. Hazelton, K. Hasegawa, R. Joseph, T. Kaneuji, K. Kubota, W. Li, C. Lynch, B. McKinley, D. A. Mitchell, M. F. Morales, B. Pindor, J. C. Pober, M. Rahimi, K. Takahashi, S. J. Tingay, R. B. Wayth, R. L. Webster, M. Wilensky , et al. (4 additional authors not shown)

Abstract: We explore methods for robust estimation of the 21 cm signal from the Epoch of Reionisation (EoR). A Kernel Density Estimator (KDE) is introduced for measuring the spatial temperature fluctuation power spectrum from the EoR. The KDE estimates the underlying probability distribution function of fluctuations as a function of spatial scale, and contains different systematic biases and errors to the t… ▽ More We explore methods for robust estimation of the 21 cm signal from the Epoch of Reionisation (EoR). A Kernel Density Estimator (KDE) is introduced for measuring the spatial temperature fluctuation power spectrum from the EoR. The KDE estimates the underlying probability distribution function of fluctuations as a function of spatial scale, and contains different systematic biases and errors to the typical approach to estimating the fluctuation power spectrum. Extraction of histograms of visibilities allows moments analysis to be used to discriminate foregrounds from 21 cm signal and thermal noise. We use the information available in the histograms, along with the statistical dis-similarity of foregrounds from two independent observing fields, to robustly separate foregrounds from cosmological signal, while making no assumptions about the Gaussianity of the signal. Using two independent observing fields to robustly discriminate signal from foregrounds is crucial for the analysis presented in this paper. We apply the techniques to 13 hours of Murchison Widefield Array (MWA) EoR data over two observing fields. We compare the output to that obtained with a comparative power spectrum estimation method, and demonstrate the reduced foreground contamination using this approach. Using the second moment obtained directly from the KDE distribution functions yields a factor of 2-3 improvement in power for k < 0.3hMpc^{-1} compared with a matched delay space power estimator, while weighting data by additional statistics does not offer significant improvement beyond that available for thermal noise-only weights. △ Less

Submitted 25 April, 2019; originally announced April 2019.

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

Dark-Ages Reionisation & Galaxy Formation Simulation XVI: The Thermal Memory of Reionisation

Authors: James E. Davies, Simon J. Mutch, Yuxiang Qin, Andrei Mesinger, Gregory B. Poole, J. Stuart B. Wyithe

Abstract: Intergalactic medium temperature is a powerful probe of the epoch of reionisation, as information is retained long after reionisation itself. However, mean temperatures are highly degenerate with the timing of reionisation, with the amount heat injected during the epoch, and with the subsequent cooling rates. We post-process a suite of semi-analytic galaxy formation models to characterise how diff… ▽ More Intergalactic medium temperature is a powerful probe of the epoch of reionisation, as information is retained long after reionisation itself. However, mean temperatures are highly degenerate with the timing of reionisation, with the amount heat injected during the epoch, and with the subsequent cooling rates. We post-process a suite of semi-analytic galaxy formation models to characterise how different thermal statistics of the intergalactic medium can be used to constrain reionisation. Temperature is highly correlated with redshift of reionisation for a period of time after the gas is heated. However as the gas cools, thermal memory of reionisation is lost, and a power-law temperature-density relation is formed, $T = T_0(1+δ)^{1-γ}$ with $γ\approx 1.5$. Constraining our model against observations of electron optical depth and temperature at mean density, we find that reionisation likely finished at $z_{\rm{reion}} = 6.8 ^{+ 0.5} _{-0.8}$ with a soft spectral slope of $α= 2.8 ^{+ 1.2} _{-1.0}$. By restricting spectral slope to the range $[0.5,2.5]$ motivated by population II synthesis models, reionisation timing is further constrained to $z_{\rm{reion}} = 6.9 ^{+ 0.4} _{-0.5}$. We find that, in the future, the degeneracies between reionisation timing and background spectrum can be broken using the scatter in temperatures and integrated thermal history. △ Less

Submitted 23 August, 2019; v1 submitted 12 April, 2019; originally announced April 2019.

Comments: 17 pages, 17 figures, Accepted for publication in MNRAS

Dark-ages Reionization and Galaxy Formation Simulation -- XVII. Sizes, angular momenta and morphologies of high redshift galaxies

Authors: Madeline A. Marshall, Simon J. Mutch, Yuxiang Qin, Gregory B. Poole, J. Stuart B. Wyithe

Abstract: We study the sizes, angular momenta and morphologies of high-redshift galaxies using an update of the Meraxes semi-analytic galaxy evolution model. Our model successfully reproduces a range of observations from redshifts $z=0$-$10$. We find that the effective radius of a galaxy disc scales with UV luminosity as $R_e\propto L_{\textrm{UV}}^{0.33}$ at $z=5$-$10$, and with stellar mass as… ▽ More We study the sizes, angular momenta and morphologies of high-redshift galaxies using an update of the Meraxes semi-analytic galaxy evolution model. Our model successfully reproduces a range of observations from redshifts $z=0$-$10$. We find that the effective radius of a galaxy disc scales with UV luminosity as $R_e\propto L_{\textrm{UV}}^{0.33}$ at $z=5$-$10$, and with stellar mass as $R_e\propto M_\ast^{0.24}$ at $z=5$ but with a slope that increases at higher redshifts. Our model predicts that the median galaxy size scales with redshift as $R_e \propto (1+z)^{-m}$, where $m=1.98\pm0.07$ for galaxies with $(0.3$-$1)L^\ast_{z=3}$ and $m=2.15\pm0.05$ for galaxies with $(0.12$-$0.3)L^\ast_{z=3}$. We find that the ratio between stellar and halo specific angular momentum is typically less than one and decreases with halo and stellar mass. This relation shows no redshift dependence, while the relation between specific angular momentum and stellar mass decreases by $\sim0.5$ dex from $z=7$ to $z=2$. Our model reproduces the distribution of local galaxy morphologies, with bulges formed predominantly through galaxy mergers for low-mass galaxies, disc-instabilities for galaxies with $M_\ast\simeq10^{10}$-$10^{11.5}M_\odot$, and major mergers for the most massive galaxies. At high redshifts, we find galaxy morphologies that are predominantly bulge-dominated. △ Less

Submitted 29 June, 2019; v1 submitted 2 April, 2019; originally announced April 2019.

Comments: 21 pages, 12 figures. Accepted for publication in MNRAS

On the observability of individual Population III stars and their stellar-mass black hole accretion disks through cluster caustic transits

Authors: Rogier A. Windhorst, M. Alpaslan, S. Andrews, T. Ashcraft, T. Broadhurst, D. Coe, S. Cohen, C. Conselice, J. Diego, M. Dijkstra, S. Driver, K. Duncan, S. Finkelstein, B. Frye, A. Griffiths, N. Grogin, N. Hathi, A. Hopkins, R. Jansen, B. Joshi, A. Kashlinsky, W. Keel, P. Kelly, D. Kim, A. Koekemoer , et al. (22 additional authors not shown)

Abstract: Recent near-IR power-spectra and panchromatic Extragalactic Background Light measurements provide upper limits on the near-IR surface brightness (SB>31 mag/arcsec^2) that may come from Pop III stars and accretion disks around resulting stellar-mass black holes (BHs) in the epoch of First Light (z=7-17). Physical parameters for zero metallicity Pop III stars at z>7 can be estimated from MESA stella… ▽ More Recent near-IR power-spectra and panchromatic Extragalactic Background Light measurements provide upper limits on the near-IR surface brightness (SB>31 mag/arcsec^2) that may come from Pop III stars and accretion disks around resulting stellar-mass black holes (BHs) in the epoch of First Light (z=7-17). Physical parameters for zero metallicity Pop III stars at z>7 can be estimated from MESA stellar evolution models through helium-depletion, and for BH accretion disks from quasar microlensing results and multicolor accretion models. Second-generation stars can form at higher multiplicity, so that BH accretion disks may be fed by Roche-lobe overflow from lower-mass companions in their AGB stage. The near-IR SB constraints can be used to calculate the number of caustic transits behind lensing clusters that JWST and the 25~39 m ground-based telescopes may detect for both Pop III stars and stellar mass BH accretion disks. Because Pop III stars and stellar mass BH accretion disks have sizes of a few x 10^-11 arcsec at z>7, typical caustic magnifications can be mu=10^4~10^5, with rise times of hours and decline times of < 1 year for cluster transverse velocities of v<1000 km/s. Microlensing by intracluster medium objects can modify transit magnifications, and lengthen visibility times. Depending on BH masses, accretion-disk radii and feeding efficiencies, stellar-mass BH accretion-disk caustic transits could outnumber those from Pop III stars. To observe Pop III caustic transits directly may require monitoring 3~30 lensing clusters to AB< 29 mag over a decade or more. Such a program must be started with JWST in Cycle 1, and -- depending on the role of microlensing in the Intra Cluster Light -- should be continued for decades with the GMT and TMT, where JWST and the ground-based telescopes each will play a unique and strongly complementary role. △ Less

Submitted 13 March, 2019; originally announced March 2019.

Comments: 6 pages + references, 3 figures, science white paper submitted to the Astro2020 decadal survey. arXiv admin note: substantial text overlap with arXiv:1801.03584

Dependence of galaxy clustering on UV-luminosity and stellar mass at $z \sim 4 - 7$

Authors: Yisheng Qiu, J. Stuart B. Wyithe, Pascal A. Oesch, Simon J. Mutch, Yuxiang Qin, Ivo Labbé, Rychard J. Bouwens, Mauro Stefanon, Garth D. Illingworth

Abstract: We investigate the dependence of galaxy clustering at $z \sim 4 - 7$ on UV-luminosity and stellar mass. Our sample consists of $\sim$ 10,000 Lyman-break galaxies (LBGs) in the XDF and CANDELS fields. As part of our analysis, the $M_\star - M_{\rm UV}$ relation is estimated for the sample, which is found to have a nearly linear slope of $d\log_{10} M_\star / d M_{\rm UV} \sim 0.44$. We subsequently… ▽ More We investigate the dependence of galaxy clustering at $z \sim 4 - 7$ on UV-luminosity and stellar mass. Our sample consists of $\sim$ 10,000 Lyman-break galaxies (LBGs) in the XDF and CANDELS fields. As part of our analysis, the $M_\star - M_{\rm UV}$ relation is estimated for the sample, which is found to have a nearly linear slope of $d\log_{10} M_\star / d M_{\rm UV} \sim 0.44$. We subsequently measure the angular correlation function and bias in different stellar mass and luminosity bins. We focus on comparing the clustering dependence on these two properties. While UV-luminosity is only related to recent starbursts of a galaxy, stellar mass reflects the integrated build-up of the whole star formation history, which should make it more tightly correlated with halo mass. Hence, the clustering segregation with stellar mass is expected to be larger than with luminosity. However, our measurements suggest that the segregation with luminosity is larger with $\simeq 90\%$ confidence (neglecting contributions from systematic errors). We compare this unexpected result with predictions from the \textsc{Meraxes} semi-analytic galaxy formation model. Interestingly, the model reproduces the observed angular correlation functions, and also suggests stronger clustering segregation with luminosity. The comparison between our observations and the model provides evidence of multiple halo occupation in the small scale clustering. △ Less

Submitted 26 September, 2018; originally announced September 2018.

Comments: 10 pages, 6 figures, 2 tables, accepted for publication in MNRAS

Assessment of ionospheric activity tolerances for Epoch of Reionisation science with the Murchison Widefield Array

Authors: Cathryn M. Trott, C. H. Jordan, S. G. Murray, B. Pindor, D. A. Mitchell, R. B. Wayth, J. Line, B. McKinley, A. Beardsley, J. Bowman, F. Briggs, B. J. Hazelton, J. Hewitt, D. Jacobs, M. F. Morales, J. C. Pober, S. Sethi, U. Shankar, R. Subrahmanyan, M. Tegmark, S. J. Tingay, R. L. Webster, J. S. B. Wyithe

Abstract: Structure imprinted in foreground extragalactic point sources by ionospheric refraction has the potential to contaminate Epoch of Reionisation (EoR) power spectra of the 21~cm emission line of neutral hydrogen. The alteration of the spatial and spectral structure of foreground measurements due to total electron content (TEC) gradients in the ionosphere create a departure from the expected sky sign… ▽ More Structure imprinted in foreground extragalactic point sources by ionospheric refraction has the potential to contaminate Epoch of Reionisation (EoR) power spectra of the 21~cm emission line of neutral hydrogen. The alteration of the spatial and spectral structure of foreground measurements due to total electron content (TEC) gradients in the ionosphere create a departure from the expected sky signal. We present a general framework for understanding the signatures of ionospheric behaviour in the two-dimensional (2D) neutral hydrogen power spectrum measured by a low-frequency radio interferometer. Two primary classes of ionospheric behaviour are considered, corresponding to dominant modes observed in Murchison Widefield Array (MWA) EoR data; namely, anisotropic structured wave behaviour, and isotropic turbulence. Analytic predictions for power spectrum bias due to this contamination are computed, and compared with simulations. We then apply the ionospheric metric described in Jordan et al. (2017) to study the impact of ionospheric structure on MWA data, by dividing MWA EoR datasets into classes with good and poor ionospheric conditions, using sets of matched 30-minute observations from 2014 September. The results are compared with the analytic and simulated predictions, demonstrating the observed bias in the power spectrum when the ionosphere is active (displays coherent structures or isotropic turbulence). The analysis demonstrates that unless ionospheric activity can be quantified and corrected, active data should not be included in EoR analysis in order to avoid systematic biases in cosmological power spectra. When data are corrected with a model formed from the calibration information, bias reduces below the expected 21~cm signal level. Data are considered `quiet' when the median measured source position offsets are less than 10-15~arcseconds. △ Less

Submitted 17 September, 2018; originally announced September 2018.

Comments: 19 pages, 11 figures, Accepted for publication in the Astrophysical Journal

The Phase II Murchison Widefield Array: Design Overview

Authors: Randall B. Wayth, Steven J. Tingay, Cathryn M. Trott, David Emrich, Melanie Johnston-Hollitt, Ben McKinley, B. M. Gaensler, A. P. Beardsley, T. Booler, B. Crosse, T. M. O. Franzen, L. Horsley, D. L. Kaplan, D. Kenney, M. F. Morales, D. Pallot, G. Sleap, K. Steele, M. Walker, A. Williams, C. Wu, Iver. H. Cairns, M. D. Filipovic, S. Johnston, T. Murphy , et al. (4 additional authors not shown)

Abstract: We describe the motivation and design details of the "Phase II" upgrade of the Murchison Widefield Array (MWA) radio telescope. The expansion doubles to 256 the number of antenna tiles deployed in the array. The new antenna tiles enhance the capabilities of the MWA in several key science areas. Seventy-two of the new tiles are deployed in a regular configuration near the existing MWA core. These n… ▽ More We describe the motivation and design details of the "Phase II" upgrade of the Murchison Widefield Array (MWA) radio telescope. The expansion doubles to 256 the number of antenna tiles deployed in the array. The new antenna tiles enhance the capabilities of the MWA in several key science areas. Seventy-two of the new tiles are deployed in a regular configuration near the existing MWA core. These new tiles enhance the surface brightness sensitivity of the MWA and will improve the ability of the MWA to estimate the slope of the Epoch of Reionisation power spectrum by a factor of ~3.5. The remaining 56 tiles are deployed on long baselines, doubling the maximum baseline of the array and improving the array u,v coverage. The improved imaging capabilities will provide an order of magnitude improvement in the noise floor of MWA continuum images. The upgrade retains all of the features that have underpinned the MWA's success (large field-of-view, snapshot image quality, pointing agility) and boosts the scientific potential with enhanced imaging capabilities and by enabling new calibration strategies. △ Less

Submitted 17 September, 2018; originally announced September 2018.

Comments: 11 pages, 8 figs. Accepted for publication in PASA

Dark-ages Reionization and Galaxy Formation Simulation - XV. Stellar evolution and feedback in dwarf galaxies at high redshift

Authors: Yuxiang Qin, Alan R. Duffy, Simon J. Mutch, Gregory B. Poole, Andrei Mesinger, J. Stuart B. Wyithe

Abstract: We directly compare predictions of dwarf galaxy properties in a semi-analytic model (SAM) with those extracted from a high-resolution hydrodynamic simulation. We focus on galaxies with halo masses of 1e9<Mvir/Msol<1e11 at high redshift ($z\ge5$). We find that, with the modifications previously proposed in Qin et al. (2018), including to suppress the halo mass and baryon fraction, as well as to mod… ▽ More We directly compare predictions of dwarf galaxy properties in a semi-analytic model (SAM) with those extracted from a high-resolution hydrodynamic simulation. We focus on galaxies with halo masses of 1e9<Mvir/Msol<1e11 at high redshift ($z\ge5$). We find that, with the modifications previously proposed in Qin et al. (2018), including to suppress the halo mass and baryon fraction, as well as to modulate gas cooling and star formation efficiencies, the SAM can reproduce the cosmic evolution of galaxy properties predicted by the hydrodynamic simulation. These include the galaxy stellar mass function, total baryonic mass, star-forming gas mass and star formation rate at $z\sim5-11$. However, this agreement is only possible by reducing the star formation threshold relative to that suggested by local observations. Otherwise, too much star-forming gas is trapped in quenched dwarf galaxies. We further find that dwarf galaxies rapidly build up their star-forming reservoirs in the early universe ($z>10$), with the relevant time-scale becoming significantly longer towards lower redshifts. This indicates efficient accretion in cold mode in these low-mass objects at high redshift. Note that the improved SAM, which has been calibrated against hydrodynamic simulations, can provide more accurate predictions of high-redshift dwarf galaxy properties that are essential for reionization study. △ Less

Submitted 10 August, 2018; originally announced August 2018.

Comments: 18 pages, 9 figures; submitted to MNRAS, comments are welcome

Comparing Redundant and Sky Model Based Interferometric Calibration: A First Look with Phase II of the MWA

Authors: W. Li, J. C. Pober, B. J. Hazelton, N. Barry, M. F. Morales, I. Sullivan, A. R. Parsons, Z. S. Ali, J. S. Dillon, A. P. Beardsley, J. D. Bowman, F. Briggs, R. Byrne, P. Carroll, B. Crosse, D. Emrich, A. Ewall-Wice, L. Feng, T. M. O. Franzen, J. N. Hewitt, L. Horsley, D. C. Jacobs, M. Johnston-Hollitt, C. Jordan, R. C. Joseph , et al. (31 additional authors not shown)

Abstract: Interferometric arrays seeking to measure the 21 cm signal from the Epoch of Reionization must contend with overwhelmingly bright emission from foreground sources. Accurate recovery of the 21 cm signal will require precise calibration of the array, and several new avenues for calibration have been pursued in recent years, including methods using redundancy in the antenna configuration. The newly u… ▽ More Interferometric arrays seeking to measure the 21 cm signal from the Epoch of Reionization must contend with overwhelmingly bright emission from foreground sources. Accurate recovery of the 21 cm signal will require precise calibration of the array, and several new avenues for calibration have been pursued in recent years, including methods using redundancy in the antenna configuration. The newly upgraded Phase II of Murchison Widefield Array (MWA) is the first interferometer that has large numbers of redundant baselines while retaining good instantaneous UV-coverage. This array therefore provides a unique opportunity to compare redundant calibration with sky-model based algorithms. In this paper, we present the first results from comparing both calibration approaches with MWA Phase II observations. For redundant calibration, we use the package OMNICAL, and produce sky-based calibration solutions with the analysis package Fast Holographic Deconvolution (FHD). There are three principal results. (1) We report the success of OMNICAL on observations of ORBComm satellites, showing substantial agreement between redundant visibility measurements after calibration. (2) We directly compare OMNICAL calibration solutions with those from FHD, and demonstrate these two different calibration schemes give extremely similar results. (3) We explore improved calibration by combining OMNICAL and FHD. We evaluate these combined methods using power spectrum techniques developed for EoR analysis and find evidence for marginal improvements mitigating artifacts in the power spectrum. These results are likely limited by signal-to-noise in the six hours of data used, but suggest future directions for combining these two calibration schemes. △ Less

Submitted 13 July, 2018; originally announced July 2018.

Comments: 20 pages, 11 figures. Accepted to ApJ

The impact of feedback and the hot halo on the rates of gas accretion onto galaxies

Authors: Camila A. Correa, Joop Schaye, Freeke van de Voort, Alan R. Duffy, J. Stuart B. Wyithe

Abstract: We investigate the physics that drives the gas accretion rates onto galaxies at the centers of dark matter haloes using the EAGLE suite of hydrodynamical cosmological simulations. We find that at redshifts $z{\le}2$ the accretion rate onto the galaxy increases with halo mass in the halo mass range $10^{10}-10^{11.7}M_{\odot}$, flattens between the halo masses $10^{11.7}-10^{12.7}M_{\odot}$, and in… ▽ More We investigate the physics that drives the gas accretion rates onto galaxies at the centers of dark matter haloes using the EAGLE suite of hydrodynamical cosmological simulations. We find that at redshifts $z{\le}2$ the accretion rate onto the galaxy increases with halo mass in the halo mass range $10^{10}-10^{11.7}M_{\odot}$, flattens between the halo masses $10^{11.7}-10^{12.7}M_{\odot}$, and increases again for higher-mass haloes. However, the galaxy gas accretion does not flatten at intermediate halo masses when AGN feedback is switched off. To better understand these trends, we develop a physically motivated semi-analytic model of galaxy gas accretion. We show that the flattening is produced by the rate of gas cooling from the hot halo. The ratio of the cooling radius and the virial radius does not decrease continuously with increasing halo mass as generally thought. While it decreases up to ${\sim}10^{13}M_{\odot}$ haloes, it increases for higher halo masses, causing an upturn in the galaxy gas accretion rate. This may indicate that in high-mass haloes AGN feedback is not sufficiently efficient. When there is no AGN feedback, the density of the hot halo is higher, the ratio of the cooling and virial radii does not decrease as much and the cooling rate is higher. Changes in the efficiency of stellar feedback can also increase or decrease the accretion rates onto galaxies. The trends can plausibly be explained by the re-accretion of gas ejected by progenitor galaxies and by the suppression of black hole growth, and hence AGN feedback, by stellar feedback. △ Less

Submitted 4 April, 2018; originally announced April 2018.

Comments: 17 pages, 9 figures. Accepted for publication in MNRAS. Comments welcome

Intensity mapping cross-correlations II: HI halo models including shot noise

Authors: L. Wolz, S. G. Murray, C. Blake, J. S. Wyithe

Abstract: HI intensity mapping data traces the large-scale structure matter distribution using the integrated emission of neutral hydrogen gas (HI). The cross-correlation of the intensity maps with optical galaxy surveys can mitigate foreground and systematic effects, but has been shown to significantly depend on galaxy evolution parameters of the HI and the optical sample. Previously, we have shown that th… ▽ More HI intensity mapping data traces the large-scale structure matter distribution using the integrated emission of neutral hydrogen gas (HI). The cross-correlation of the intensity maps with optical galaxy surveys can mitigate foreground and systematic effects, but has been shown to significantly depend on galaxy evolution parameters of the HI and the optical sample. Previously, we have shown that the shot noise of the cross-correlation scales with the HI content of the optical samples, such that the shot noise estimation infers the average HI masses of these samples. In this article, we present an adaptive framework for the cross-correlation of HI intensity maps with galaxy samples using our implementation of the halo model formalism (Murray et al 2018, in prep) which utilises the halo occupation distribution of galaxies to predict their power spectra. We compare two HI population models, tracing the spatial halo and the galaxy distribution respectively, and present their auto- and cross-power spectra with an associated galaxy sample. We find that the choice of the HI model and the distribution of the HI within the galaxy sample have minor significance for the shape of the auto- and cross-correlations, but highly impact the measured shot noise amplitude of the estimators, a finding we confirm with simulations. We demonstrate parameter estimation of the HI halo occupation models and advocate this framework for the interpretation of future experimental data, with the prospect of determining the HI masses of optical galaxy samples via the cross-correlation shot noise. △ Less

Submitted 6 March, 2018; originally announced March 2018.

Comments: 15 pages, 8 figures, 3 tables. Comments welcome

Dark-ages Reionization and Galaxy Formation Simulation - XIV. Gas accretion, cooling and star formation in dwarf galaxies at high redshift

Authors: Yuxiang Qin, Alan R. Duffy, Simon J. Mutch, Gregory B. Poole, Paul M. Geil, Andrei Mesinger, J. Stuart B. Wyithe

Abstract: We study dwarf galaxy formation at high redshift ($z\ge5$) using a suite of high- resolution, cosmological hydrodynamic simulations and a semi-analytic model (SAM). We focus on gas accretion, cooling and star formation in this work by isolating the relevant process from reionization and supernova feedback, which will be further discussed in a companion paper. We apply the SAM to halo merger trees… ▽ More We study dwarf galaxy formation at high redshift ($z\ge5$) using a suite of high- resolution, cosmological hydrodynamic simulations and a semi-analytic model (SAM). We focus on gas accretion, cooling and star formation in this work by isolating the relevant process from reionization and supernova feedback, which will be further discussed in a companion paper. We apply the SAM to halo merger trees constructed from a collisionless N-body simulation sharing identical initial conditions to the hydrodynamic suite, and calibrate the free parameters against the stellar mass function predicted by the hydrodynamic simulations at z = 5. By making comparisons of the star formation history and gas components calculated by the two modelling techniques, we find that semi-analytic prescriptions that are commonly adopted in the literature of low-redshift galaxy formation do not accurately represent dwarf galaxy properties in the hydrodynamic simulation at earlier times. We propose 3 modifications to SAMs that will provide more accurate high-redshift simulations. These include 1) the halo mass and baryon fraction which are overestimated by collisionless N-body simulations; 2) the star formation efficiency which follows a different cosmic evolutionary path from the hydrodynamic simulation; and 3) the cooling rate which is not well defined for dwarf galaxies at high redshift. Accurate semi-analytic modelling of dwarf galaxy formation informed by detailed hydrodynamical modelling will facilitate reliable semi-analytic predictions over the large volumes needed for the study of reionization. △ Less

Submitted 6 August, 2018; v1 submitted 11 February, 2018; originally announced February 2018.

Comments: 20 pages, 10 figures; Updated to match the published version. All results and conclusions remain unchanged

Journal ref: Monthly Notices of the Royal Astronomical Society, 2018, Volume 477, Issue 1, p.1318-1335

On the observability of individual Population III stars and their stellar-mass black hole accretion disks through cluster caustic transits

Authors: Rogier A. Windhorst, F. X. Timmes, J. Stuart B. Wyithe, Mehmet Alpaslan, Stephen K. Andrews, Daniel Coe, Jose M. Diego, Mark Dijkstra, Simon P. Driver, Patrick L. Kelly, Duho Kim

Abstract: We summarize panchromatic Extragalactic Background Light data to place upper limits on the integrated near-infrared surface brightness (SB) that may come from Population III stars and possible accretion disks around their stellar-mass black holes (BHs) in the epoch of First Light, broadly taken from z$\simeq$7-17. Theoretical predictions and recent near-infrared power-spectra provide tighter const… ▽ More We summarize panchromatic Extragalactic Background Light data to place upper limits on the integrated near-infrared surface brightness (SB) that may come from Population III stars and possible accretion disks around their stellar-mass black holes (BHs) in the epoch of First Light, broadly taken from z$\simeq$7-17. Theoretical predictions and recent near-infrared power-spectra provide tighter constraints on their sky-signal. We outline the physical properties of zero metallicity Population III stars from MESA stellar evolution models through helium-depletion and of BH accretion disks at z$\gtrsim$7. We assume that second-generation non-zero metallicity stars can form at higher multiplicity, so that BH accretion disks may be fed by Roche-lobe overflow from lower-mass companions. We use these near-infrared SB constraints to calculate the number of caustic transits behind lensing clusters that the James Webb Space Telescope and the next generation ground-based telescopes may observe for both Population III stars and their BH accretion disks. Typical caustic magnifications can be $μ$$\simeq$10$^4$-10$^5$, with rise times of hours and decline times of $\lesssim$1 year for cluster transverse velocities of $v_{T}$$\lesssim$1000 km s$^{-1}$. Microlensing by intracluster medium objects can modify transit magnifications, but lengthen visibility times. Depending on BH masses, accretion-disk radii and feeding efficiencies, stellar-mass BH accretion-disk caustic transits could outnumber those from Population III stars. To observe Population III caustic transits directly may require to monitor 3-30 lensing clusters to AB$\lesssim$29 mag over a decade. △ Less

Submitted 10 January, 2018; originally announced January 2018.

Comments: 53 pages, 5 figures, Accepted for publication in ApJS

MUSE spectroscopy and deep observations of a unique compact JWST target, lensing cluster CLIO

Authors: Alex Griffiths, Christopher J. Conselice, Mehmet Alpaslan, Brenda L. Frye, Jose M. Diego, Adi Zitrin, Haojing Yan, Zhiyuan Ma, Robert Barone-Nugent, Rachana Bhatawdekar, Simon. P. Driver, Aaron S. G. Robotham, Rogier A. Windhorst, J. Stuart B. Wyithe

Abstract: We present the results of a VLT MUSE/FORS2 and Spitzer survey of a unique compact lensing cluster CLIO at z = 0.42, discovered through the GAMA survey using spectroscopic redshifts. Compact and massive clusters such as this are understudied, but provide a unique prospective on dark matter distributions and for finding background lensed high-z galaxies. The CLIO cluster was identified for follow up… ▽ More We present the results of a VLT MUSE/FORS2 and Spitzer survey of a unique compact lensing cluster CLIO at z = 0.42, discovered through the GAMA survey using spectroscopic redshifts. Compact and massive clusters such as this are understudied, but provide a unique prospective on dark matter distributions and for finding background lensed high-z galaxies. The CLIO cluster was identified for follow up observations due to its almost unique combination of high mass and dark matter halo concentration, as well as having observed lensing arcs from ground based images. Using dual band optical and infra-red imaging from FORS2 and Spitzer, in combination with MUSE optical spectroscopy we identify 89 cluster members and find background sources out to z = 6.49. We describe the physical state of this cluster, finding a strong correlation between environment and galaxy spectral type. Under the assumption of a NFW profile, we measure the total mass of CLIO to be M$_{200} = (4.49 \pm 0.25) \times 10^{14}$ M$_\odot$. We build and present an initial strong-lensing model for this cluster, and measure a relatively low intracluster light (ICL) fraction of 7.21 $\pm$ 1.53% through galaxy profile fitting. Due to its strong potential for lensing background galaxies and its low ICL, the CLIO cluster will be a target for our 110 hour JWST 'Webb Medium-Deep Field' (WMDF) GTO program. △ Less

Submitted 3 January, 2018; originally announced January 2018.

Comments: Accepted at MNRAS, 19 pages, 16 figures

Journal ref: Monthly Notices of the Royal Astronomical Society, Volume 475, Issue 3, p.2853-2869. 04/2018

The SAMI Galaxy Survey: understanding observations of large-scale outflows at low redshift with EAGLE simulations

Authors: E. Tescari, L. Cortese, C. Power, J. S. B. Wyithe, I. -T. Ho, R. A. Crain, J. Bland-Hawthorn, S. M. Croom, L. J. Kewley, J. Schaye, R. G. Bower, T. Theuns, M. Schaller, L. Barnes, S. Brough, J. J. Bryant, M. Goodwin, M. L. P. Gunawardhana, J. S. Lawrence, S. K. Leslie, Á. R. López-Sánchez, N. P. F. Lorente, A. M. Medling, S. N. Richards, S. M. Sweet , et al. (1 additional authors not shown)

Abstract: This work presents a study of galactic outflows driven by stellar feedback. We extract main sequence disc galaxies with stellar mass $10^9\le$ M$_{\star}/$M$_{\odot} \le 5.7\times10^{10}$ at redshift $z=0$ from the highest resolution cosmological simulation of the Evolution and Assembly of GaLaxies and their Environments (EAGLE) set. Synthetic gas rotation velocity and velocity dispersion ($σ$) ma… ▽ More This work presents a study of galactic outflows driven by stellar feedback. We extract main sequence disc galaxies with stellar mass $10^9\le$ M$_{\star}/$M$_{\odot} \le 5.7\times10^{10}$ at redshift $z=0$ from the highest resolution cosmological simulation of the Evolution and Assembly of GaLaxies and their Environments (EAGLE) set. Synthetic gas rotation velocity and velocity dispersion ($σ$) maps are created and compared to observations of disc galaxies obtained with the Sydney-AAO Multi-object Integral field spectrograph (SAMI), where $σ$-values greater than $150$ km s$^{-1}$ are most naturally explained by bipolar outflows powered by starburst activity. We find that the extension of the simulated edge-on (pixelated) velocity dispersion probability distribution depends on stellar mass and star formation rate surface density ($Σ_{\rm SFR}$), with low-M$_{\star}/$low-$Σ_{\rm SFR}$ galaxies showing a narrow peak at low $σ$ ($\sim30$ km s$^{-1}$) and more active, high-M$_{\star}/$high-$Σ_{\rm SFR}$ galaxies reaching $σ>150$ km s$^{-1}$. Although supernova-driven galactic winds in the EAGLE simulations may not entrain enough gas with T $<10^5$ K compared to observed galaxies, we find that gas temperature is a good proxy for the presence of outflows. There is a direct correlation between the thermal state of the gas and its state of motion as described by the $σ$-distribution. The following equivalence relations hold in EAGLE: $i)$ low-$σ$ peak $\,\Leftrightarrow\,$ disc of the galaxy $\,\Leftrightarrow\,$ gas with T $<10^5$ K; $ii)$ high-$σ$ tail $\,\Leftrightarrow\,$ galactic winds $\,\Leftrightarrow\,$ gas with T $\ge 10^5$ K. △ Less

Submitted 6 September, 2017; originally announced September 2017.

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

The formation of hot gaseous haloes around galaxies

Authors: Camila A. Correa, Joop Schaye, J. Stuart B. Wyithe, Alan R. Duffy, Tom Theuns, Robert A. Crain, Richard Bower

Abstract: We use a suite of hydrodynamical cosmological simulations from the Evolution and Assembly of GaLaxies and their Environments (EAGLE) project to investigate the formation of hot hydrostatic haloes and their dependence on feedback mechanisms. We find that the appearance of a strong bimodality in the probability density function (PDF) of the ratio of the radiative cooling and dynamical times for halo… ▽ More We use a suite of hydrodynamical cosmological simulations from the Evolution and Assembly of GaLaxies and their Environments (EAGLE) project to investigate the formation of hot hydrostatic haloes and their dependence on feedback mechanisms. We find that the appearance of a strong bimodality in the probability density function (PDF) of the ratio of the radiative cooling and dynamical times for halo gas provides a clear signature of the formation of a hot corona. Haloes of total mass $10^{11.5}-10^{12}{\rm{M}}_{\odot}$ develop a hot corona independent of redshift, at least in the interval $z=0-4$ where the simulation has sufficiently good statistics. We analyse the build up of the hot gas mass in the halo, $M_{\rm{hot}}$, as a function of halo mass and redshift and find that while more energetic galactic winds powered by SNe increases $M_{\rm{hot}}$, AGN feedback reduces it by ejecting gas from the halo. We also study the thermal properties of gas accreting onto haloes and measure the fraction of shock-heated gas as a function of redshift and halo mass. We develop analytic and semianalytic approaches to estimate a `critical halo mass', $M_{\rm{crit}}$, for hot halo formation. We find that the mass for which the heating rate produced by accretion shocks equals the radiative cooling rate, reproduces the mass above which haloes develop a significant hot atmosphere. This yields a mass estimate of $M_{\rm{crit}} \approx 10^{11.7}{\rm{M}}_{\odot}$ at $z=0$, which agrees with the simulation results. The value of $M_{\rm{crit}}$ depends more strongly on the cooling rate than on any of the feedback parameters. △ Less

Submitted 6 September, 2017; originally announced September 2017.

Comments: 24 pages, 15 figures, accepted for publication in MNRAS. Comments welcome