Search

In this paper, we use hydrodynamic zoom-in simulations of Milky Way-type haloes to explore using dust as an observational tracer to discriminate between cold and warm dark matter (WDM) universes. Comparing a cold and 3.5 keV WDM particle model, we tune the efficiency of galaxy formation in our simulations using a variable supernova rate to create Milky Way systems with similar satellite galaxy populations while keeping all other simulation parameters the same. Cold dark matter (CDM), having more substructure, requires a higher supernova efficiency than WDM to achieve the same satellite galaxy number. These different supernova efficiencies create different dust distributions around their host galaxies, which we generate by post-processing the simulation output with the \powderday{} codebase. Analysing the resulting dust in each simulation, we find $\sim$4.5 times more dust in our CDM Milky Way haloes compared with WDM. The distribution of dust out to $R_{200\text{c}}$ is then explored, revealing that the WDM simulations are noticeably less concentrated than their CDM counterparts, although differences in substructure complicate the comparison. Our results indicate that dust is a possible unique probe to test theories of dark matter., Comment: Accepted for publication to MNRAS, 12 pages, 5 figures

The Sagittarius Dwarf Spheroidal galaxy (Sgr) is investigated as a target for DM annihilation searches utilising J-factor distributions calculated directly from a high-resolution hydrodynamic simulation of the infall and tidal disruption of Sgr around the Milky Way. In contrast to past studies, the simulation incorporates DM, stellar and gaseous components for both the Milky Way and the Sgr progenitor galaxy. The simulated distributions account for significant tidal disruption affecting the DM density profile. Our estimate of the J-factor value for Sgr, $J_{\text{Sgr}}=1.48\times 10^{10}$ M$_\odot^2$ kpc$^{-5}$ ($6.46\times10^{16}\ \text{GeV}\ \text{cm}^{-5}$), is significantly lower than found in prior studies. This value, while formally a lower limit, is likely close to the true J-factor value for Sgr. It implies a DM cross-section incompatibly large in comparison with existing constraints would be required to attribute recently observed $\gamma$-ray emission from Sgr to DM annihilation. We also calculate a J-factor value using a NFW profile fitted to the simulated DM density distribution to facilitate comparison with past studies. This NFW J-factor value supports the conclusion that most past studies have overestimated the dark matter density of Sgr on small scales. This, together with the fact that the Sgr has recently been shown to emit $\gamma$-rays of astrophysical origin, complicate the use of Sgr in indirect DM detection searches., Comment: 15 pages, 15 figures. Accepted for publication in MNRAS. Accepted manuscript uploaded

We use high-resolution, hydrodynamic, galaxy simulations from the Latte suite of FIRE-2 simulations to investigate the inherent variation of dark matter in sub-sampled regions around the Solar Circle of a Milky Way-type analogue galaxy and its impact on direct dark matter detection. These simulations show that the baryonic backreaction, as well as the assembly history of substructures, has lasting impacts on the dark matter's spatial and velocity distributions. These are experienced as 'gusts' of dark matter wind around the Solar Circle, potentially complicating interpretations of direct detection experiments on Earth. We find that the velocity distribution function in the galactocentric frame shows strong deviations from the Maxwell Boltzmann form typically assumed in the fiducial Standard Halo Model, indicating the presence of high-velocity substructures. By introducing a new numerical integration technique which removes any dependencies on the Standard Halo Model, we generate event-rate predictions for both single-element Germanium and compound Sodium Iodide detectors, and explore how the variability of dark matter around the Solar Circle influences annual modulation signal predictions. We find that these velocity substructures contribute additional astrophysical uncertainty to the interpretation of event rates, although their impact on summary statistics such as the peak day of annual modulation is generally low., Comment: 22 pages, 27 figures, 4 tables. Submitted to MNRAS. Comments welcome

The Fermi Bubbles are giant, gamma-ray emitting lobes emanating from the nucleus of the Milky Way discovered in ~1-100 GeV data collected by the Large Area Telescope on board the Fermi Gamma-Ray Space Telescope. Previous work has revealed substructure within the Fermi Bubbles that has been interpreted as a signature of collimated outflows from the Galaxy's super-massive black hole. Here we show via a spatial template analysis that much of the gamma-ray emission associated to the brightest region of substructure -- the so-called cocoon -- is likely due to the Sagittarius dwarf spheroidal (Sgr dSph) galaxy. This large Milky Way satellite is viewed through the Fermi Bubbles from the position of the Solar System. As a tidally and ram-pressure stripped remnant, the Sgr dSph has no on-going star formation, but we nevertheless demonstrate that the dwarf's millisecond pulsar (MSP) population can plausibly supply the gamma-ray signal that our analysis associates to its stellar template. The measured spectrum is naturally explained by inverse Compton scattering of cosmic microwave background photons by high-energy electron-positron pairs injected by MSPs belonging to the Sgr dSph, combined with these objects' magnetospheric emission. This finding plausibly suggests that MSPs produce significant gamma-ray emission amongst old stellar populations, potentially confounding indirect dark matter searches in regions such as the Galactic Centre, the Andromeda galaxy, and other massive Milky Way dwarf spheroidals., Comment: Updated to match version accepted for publication in Nature Astronomy (2022). 14 pages main text, 3 main figures, 7 extended data figures. For published version of the paper, see https://rdcu.be/cUZ6X

Fast Radio Bursts (FRBs) are a promising tool for studying the low-density universe as their dispersion measures (DM) are extremely sensitive probes of electron column density. Active Galactic Nuclei (AGN) inject energy into the intergalactic medium, affecting the DM and its scatter. To determine the effectiveness of FRBs as a probe of AGN feedback, we analysed three different AGN models from the EAGLE simulation series. We measured the mean DM-redshift relation, and the scatter around it, using $2.56 \times 10^8$ sightlines at 131 redshift ($z$) bins between $0 \leq z \leq 3$. While the DM-redshift relation itself is highly robust against different AGN feedback models, significant differences are detected in the scatter around the mean: weaker feedback leads to more scatter. We find $\sim 10^4$ localised FRBs are needed to discriminate between the scatter in standard feedback and stronger, more intermittent feedback models. The number of FRBs required is dependent on the redshift distribution of the detected population. A log-normal redshift distribution at $z=0.5$ requires approximately 50% fewer localised FRBs than a distribution centred at $z=1$. With the Square Kilometre Array expected to detect $>10^3$ FRBs per day, in the future, FRBs will be able to provide constraints on AGN feedback., Comment: 5 pages, 4 figures; Submitted to MNRAS Letters. Comments welcome

The lunar surface is extremely harsh and current mineral processing and metal extraction technologies are not adequately equipped to address this environment. In this paper we review the metals available for extraction and conditions at the lunar surface, and analyse the challenges associated with comminution, beneficiation, and metal extraction operations. The potential beneficial effects of the natural lunar conditions are also evaluated. This investigation concludes that process plant design on the lunar surface will favour lightweight, schematically simple flow sheets that enable automation, and that utilise the local environment wherever possible. The elimination of traditional comminution and beneficiation stages and their replacement with basic classification could be economically favourable. The most promising metal reduction pathways are identified as molten regolith electrolysis, and vacuum thermal dissociation, other processes with merit are hydrogen reduction, carbothermal reduction, and solid electrolysis. Finally, it is identified that a significant research effort in all areas of astrometallurgy will be required before industrial-sized extra-terrestrial mineral processing and metal extraction operations will be viable., Comment: 31 Pages, 5 Figures

Modern deep Super-Resolution (SR) networks have established themselves as valuable techniques in image reconstruction and enhancement. However, these networks are normally trained and tested on benchmark image data that lacks the typical image degrading noise present in real images. In this paper, we test the feasibility of using deep SR in real remote sensing payloads by assessing SR performance in reconstructing realistically degraded satellite images. We demonstrate that a state-of-the-art SR technique called Enhanced Deep Super-Resolution Network (EDSR), without domain specific pre-training, can recover encoded pixel data on images with poor ground sampling distance, provided the ground resolved distance is sufficient. However, this recovery varies amongst selected geographical types. Our results indicate that custom training has potential to further improve reconstruction of overhead imagery, and that new satellite hardware should prioritise optical performance over minimising pixel size as deep SR can overcome a lack of the latter but not the former., Comment: 6 pages, 6 figures, for supplementary results, see https://smpetrie.github.io/superres/

The dispersion measure (DM) of fast radio bursts (FRBs) provides a unique way to probe ionised baryons in the intergalactic medium (IGM). Cosmological models with different parameters lead to different DM-redshift ($\mathrm{DM}-z$) relations. Additionally, the over/under-dense regions in the IGM and the circumgalactic medium of intervening galaxies lead to scatter around the mean $\mathrm{DM}-z$ relations. We have used the Evolution and Assembly of GaLaxies and their Environments (EAGLE) simulations to measure the mean $\mathrm{DM}-z$ relation and the scatter around it using over one billion lines-of-sight between redshifts $0

We demonstrate the potential of an emulator-based approach to analyzing galaxy formation models in the domain where constraining data is limited. We have applied the open-source Python package PRISM to the galaxy formation model Meraxes. Meraxes is a semi-analytic model, purposefully built to study the growth of galaxies during the Epoch of Reionization (EoR). Constraining such models is however complicated by the scarcity of observational data in the EoR. PRISM's ability to rapidly construct accurate approximations of complex scientific models using minimal data is therefore key to performing this analysis well. This paper provides an overview of our analysis of Meraxes using measurements of galaxy stellar mass densities; luminosity functions; and color-magnitude relations. We demonstrate the power of using PRISM instead of a full Bayesian analysis when dealing with highly correlated model parameters and a scarce set of observational data. Our results show that the various observational data sets constrain Meraxes differently and do not necessarily agree with each other, signifying the importance of using multiple observational data types when constraining such models. Furthermore, we show that PRISM can detect when model parameters are too correlated or cannot be constrained effectively. We conclude that a mixture of different observational data types, even when they are scarce or inaccurate, is a priority for understanding galaxy formation and that emulation frameworks like PRISM can guide the selection of such data., Comment: 19 pages, 120 figures, accepted for publication in ApJS

The Cosmic Ray Extremely Distributed Observatory (CREDO) is a newly formed, global collaboration dedicated to observing and studying cosmic rays (CR) and cosmic ray ensembles (CRE): groups of a minimum of two CR with a common primary interaction vertex or the same parent particle. The CREDO program embraces testing known CR and CRE scenarios, and preparing to observe unexpected physics, it is also suitable for multi-messenger and multi-mission applications. Perfectly matched to CREDO capabilities, CRE could be formed both within classical models (e.g. as products of photon-photon interactions), and exotic scenarios (e.g. as results of decay of Super Heavy Dark Matter particles). Their fronts might be significantly extended in space and time, and they might include cosmic rays of energies spanning the whole cosmic ray energy spectrum, with a footprint composed of at least two extensive air showers with correlated arrival directions and arrival times. Since CRE are mostly expected to be spread over large areas and, because of the expected wide energy range of the contributing particles, CRE detection might only be feasible when using available cosmic ray infrastructure collectively, i.e. as a globally extended network of detectors. Thus, with this review article, the CREDO Collaboration invites the astroparticle physics community to actively join or to contribute to the research dedicated to CRE, and in particular to share any cosmic ray data useful for the specific CRE detection strategies., Comment: 59 pages, 23 figures

The Cosmic-Ray Extremely Distributed Observatory (CREDO) is a project created a few years ago in the Institute of Nuclear Physics PAS in Krak\'ow and dedicated is to global studies of extremely extended cosmic-ray phenomena. The main reason for creating such a project was that the cosmic-ray ensembles (CRE) are beyond the capabilities of existing detectors and observatories. Until now, cosmic ray studies, even in major observatories, have been limited to the recording and analysis of individual air showers therefore ensembles of cosmic-rays, which may spread over a significant fraction of the Earth were neither recorded nor analyzed. In this paper the status and perspectives of the CREDO project are presented., Comment: 8 pages, 5 figures

We study the feasibility of detecting preshower initiated by ultra-high energy photons using Monte-Carlo simulations of nearly horizontal air showers for the example of the La Palma site of the Cherenkov Telescope Array. We investigate the efficiency of multivariate analysis in correctly identifying preshower events initiated by 40 EeV photons and cosmic-ray dominated background simulated in the energy range 10 TeV -- 10 EeV. The effective areas for such kind of events are also investigated and event rate predictions related to different ultra-high energy photons production models are presented. While the expected number of preshowers from diffuse emission of UHE photon for 30 hours of observation is estimated around $3.3\times10^{-5}$ based on the upper limits put by the Pierre Auger Observatory, this value is at the level of $2.7\times10^{-4}$ ($5.7\times 10^{-5}$) when considering the upper limits of the Pierre Auger Observatory (Telescope Array) on UHE photon point sources. However, UHE photon emission may undergo possible "boosting" due to gamma-ray burst, increasing the expected number of preshower events up to 0.17 and yielding a minimum required flux of $\sim 0.2$ $\mathrm{km^{-2}yr^{-1}}$ to obtain one preshower event, which is about a factor 10 higher than upper limits put by the Pierre Auger Observatory and Telescope Array (0.034 and 0.019 $\mathrm{km^{-2}yr^{-1}}$, respectively)., Comment: 16 pages, 16 figures, 4 tables

We use KiDS weak lensing data to measure variations in mean halo mass as a function of several key galaxy properties (namely: stellar colour, specific star formation rate, Sersic index, and effective radius) for a volume-limited sample of GAMA galaxies in a narrow stellar mass range ($M_* \sim 2$--$5 \times 10^{10}$ Msol). This mass range is particularly interesting, inasmuch as it is where bimodalities in galaxy properties are most pronounced, and near to the break in both the galaxy stellar mass function and the stellar-to-halo mass relation (SHMR). In this narrow mass range, we find that both size and Sersic index are better predictors of halo mass than either colour or SSFR, with the data showing a slight preference for Sersic index. In other words, we find that mean halo mass is more tightly correlated with galaxy structure than either past star formation history or current star formation rate. Our results lead to an approximate lower bound on the dispersion in halo masses among $\log M_* \approx {10.5}$ galaxies: we find that the dispersion is $\gtrsim 0.3$ dex. This would imply either that offsets from the mean SHMR are closely coupled to size/structure, or that the dispersion in the SHMR is larger than past results have suggested. Our results thus provide new empirical constraints on the relationship between stellar and halo mass assembly at this particularly interesting mass range., Comment: 16 pages, 7 figures; submitted to MNRAS; updated to reflect accepted version

The measurement of an annual modulation in the event rate of direct dark matter detection experiments is a powerful tool for dark matter discovery. Indeed, several experiments have already claimed such a discovery in the past decade. While most of them have later revoked their conclusions, and others have found potentially contradictory results, one still stands today. This paper explains the potential as well as the challenges of annual modulation measurements, and gives an overview of past, present and future direct detection experiments., Comment: Accepted review article in Journal of Physics G: Nuclear and Particle Physics, 49 pages, 3 figures, v3 updated to accepted version

Muons and other ionizing radiation produced by cosmic rays and radiative decays affect CMOS/CCD sensor. When particles colliding with sensors atoms cause specific kind of noise on images recorded by cameras. We present a concept and preliminary implementation of method for recognizing those events and algorithms for image processing and their classification by machine learning. Our method consists of analyzing the shape of traces present in images recorded by a camera sensor and metadata related to an image like camera model, GPS location of camera, vertical and horizontal orientation of a camera sensor, timestamp of image acquisition, and other events recognized near-by sensors. The so created feature vectors are classified as either a muon-like event, an electron-like event or the other event, possibly noise. For muon-like events our method estimates azimuth of a muon track. Source of the data is database of CREDO (Cosmic-Ray Extremely Distributed Observatory) project and ESO (European Southern Observatory) archives. The telescope dark frames from ESO are analysed. CREDO project collected so far over 2 millions images of events from many kinds of cameralike: smartphones camera, laptop webcams and Internet of Things cameras localised around the globe., Comment: 36th International Cosmic Ray Conference (ICRC2019)

The search for Cosmic-Ray Ensembles (CRE), groups of correlated cosmic rays that might be distributed over very large areas, even of the size of the planet, requires a globally spread and dense network of detectors, as proposed by the Cosmic-Ray Extremely Distributed Observatory (CREDO) Collaboration. This proposal motivates an effort towards exploring the potential of using even very much diversified detection technologies within one system, with detection units located even in hard-to-reach places, where, nevertheless, the sensors could work independently - without human intervention. For these reasons we have developed a dedicated communication solution enabling the connection of many different types of detectors, in a range of environments. The proposed data transmission system uses radio waves as an information carrier on the 169MHz frequency band in a contrast to the typical commercially used frequencies in a IoT systems (868MHz). The connectivity within the system is based on the star topology, which ensures the least energy consumption. The solution is now being prepared to being implemented using the prototype detection system based on the CosmicWatch open hardware design: a portable, pocket size, and economy particle detector using the scintillation technique. Our prototype detector is equipped with a dedicated software that integrates it with the already operational CREDO server system., Comment: ICRC2019

The Cosmic-Ray Extremely Distributed Observatory (CREDO) uses the hunt for particle cascades from deep space as a vehicle for a unique "bottom-up" approach to scientific research. By engaging the non-specialist public of all ages as "citizen scientists" we create opportunities for lifelong learning for individuals as well as for cooperation and the sharing of common educational tools amongst institutions. The discoveries of these citizen scientists will feed directly into a pioneering new area of scientific research oriented on Cosmic Ray Ensembles (CRE). The detection (or non-detection) of such particle groups promises to open up a new method for exploring our universe, and a new channel on the multi-messenger stage, oriented on both astro- and geo-investigations. The opportunities this would create for cross-disciplinary research are significant and beneficial for individuals, networks of institutions and the global communities of both professional scientists and science enthusiasts., Comment: 4 pages, 3 figures, Proceedings The New Era of Multi-Messenger Astrophysics - Asterics2019, 25 - 29 March, 2019 Groningen, The Netherlands

Propagation of ultra-high energy photons in the galactic and intergalactic space gives rise to cascades comprising thousands of photons. Using Monte Carlo simulations, we investigate the development of such cascades in the solar magnetosphere, and find that the photons in the cascades are distributed over hundreds of kilometers as they arrive at the top of the Earth's atmosphere. We also perform similar study for cascades starting as far as 10 Mpc away from us using relevant magnetic field models. A few photons correlated in time are expected to arrive at the Earth from the latter type of cascade. We present our simulation results and discuss the prospects for detection of these cascades with the Cosmic-Ray Extremely Distributed Observatory.

The Cosmic-Ray Extremely Distributed Observatory (CREDO) is a project dedicated to global studies of extremely extended cosmic-ray phenomena, the cosmic-ray ensembles (CRE), beyond the capabilities of existing detectors and observatories. Up to date cosmic-ray research has been focused on detecting single air showers, while the search for ensembles of cosmic-rays, which may spread over a significant fraction of the Earth, is a scientific terra incognita. The key idea of CREDO is to combine existing cosmic-ray detectors (large professional arrays, educational instruments, individual detectors, such as smartphones, etc.) into a worldwide network, thus enabling a global analysis. The second goal of CREDO involves a large number of participants (citizen science!), assuring the geographical spread of the detectors and managing manpower necessary to deal with vast amount of data to search for evidence for cosmic-ray ensembles. In this paper the status and perspectives of the project are presented., Comment: Proceedings of 36th International Cosmic Ray Conference - ICRC2019 - July 24th - August 1st, 2019 Madison, WI, U.S.A

We have built PRISM, a "Probabilistic Regression Instrument for Simulating Models". PRISM uses the Bayes linear approach and history matching to construct an approximation ('emulator') of any given model, by combining limited model evaluations with advanced regression techniques, covariances and probability calculations. It is designed to easily facilitate and enhance existing Markov chain Monte Carlo (MCMC) methods by restricting plausible regions and exploring parameter space efficiently. However, PRISM can additionally be used as a standalone alternative to MCMC for model analysis, providing insight into the behavior of complex scientific models. With PRISM, the time spent on evaluating a model is minimized, providing developers with an advanced model analysis for a fraction of the time required by more traditional methods. This paper provides an overview of the different techniques and algorithms that are used within PRISM. We demonstrate the advantage of using the Bayes linear approach over a full Bayesian analysis when analyzing complex models. Our results show how much information can be captured by PRISM and how one can combine it with MCMC methods to significantly speed up calibration processes (>15 times faster). PRISM is an open-source Python package that is available under the BSD 3-Clause License (BSD-3) at https://github.com/1313e/PRISM and hosted at https://prism-tool.readthedocs.io. PRISM has also been reviewed by "The Journal of Open Source Software" (https://doi.org/10.21105/joss.01229)., Comment: 28 pages, 13 figures, 1 table. Updated to reflect some changes made to published version

One of the main objectives of cosmic-ray studies are precise measurements of energy and chemical composition of particles with extreme energies. Large and sophisticated detectors are used to find events seen as showers starting in the Earth's atmosphere with recorded energies larger than 100 EeV. However, a Cosmic-Ray Ensemble (CRE) developing before reaching the Earth as a bunch of correlated particles may spread over larger areas and requires an extended set of detectors to be discovered. The Cosmic-Ray Extremely Distributed Observatory (CREDO) is a solution to find such phenomena. Even simple detectors measuring the particle arrival time only are useful in this approach, as they are sufficient both to provide candidate CRE events and to determine the direction from which they are arriving., Comment: 20th International Symposium on Very High Energy Cosmic Ray Interactions, May 21-15, 2018, Nagoya, Japan

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 1e910$), 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., Comment: 18 pages, 9 figures; submitted to MNRAS, comments are welcome

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., Comment: 17 pages, 9 figures. Accepted for publication in MNRAS. Comments welcome

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., Comment: 20 pages, 10 figures; Updated to match the published version. All results and conclusions remain unchanged

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., Comment: 24 pages, 15 figures, accepted for publication in MNRAS. Comments welcome

Using a series of high-resolution hydrodynamical simulations we show that during the rapid growth of high-redshift (z > 5) galaxies, reserves of molecular gas are consumed over a time-scale of 300Myr, almost independent of feedback scheme. We find that there exists no such simple relation for the total gas fractions of these galaxies, with little correlation between gas fractions and specific star formation rates. The bottleneck or limiting factor in the growth of early galaxies is in converting infalling gas to cold star-forming gas. Thus, we find that the majority of high redshift dwarf galaxies are effectively in recession, with demand (of star formation) never rising to meet supply (of gas), irrespective of the baryonic feedback physics modelled. We conclude that the basic assumption of self-regulation in galaxies - that they can adjust total gas consumption within a Hubble time - does not apply for the dwarf galaxies thought to be responsible for providing most UV photons to reionize the high redshift Universe. We demonstrate how this rapid molecular time-scale improves agreement between semi-analytic model predictions of the early Universe and observed stellar mass functions., Comment: 17 pages, 27 figures, accepted for publication in MNRAS, minor updates to align with final published version

The SKA and its pathfinders will enable studies of HI emission at higher redshifts than ever before. In moving beyond the local Universe, this will require the use of cosmologically appropriate formulae that have traditionally been simplified to their low-redshift approximations. In this paper, we summarise some of the most important relations for tracing HI emission in the SKA era, and present an online calculator to assist in the planning and analysis of observations (hifi.icrar.org)., Comment: submitted to PASA

Direct detection of regions of ionized hydrogen (HII) has been suggested as a promising probe of cosmic reionization. Observing the redshifted 21-cm signal of hydrogen from the epoch of reionization (EoR) is a key scientific driver behind new-generation, low-frequency radio interferometers. We investigate the feasibility of combining low-frequency observations with the Square Kilometre Array and near infra-red survey data of the Wide-Field Infrared Survey Telescope to detect cosmic reionization by imaging HII bubbles surrounding massive galaxies during the cosmic dawn. While individual bubbles will be too small to be detected, we find that by stacking redshifted 21-cm spectra centred on known galaxies, it will be possible to directly detect the EoR at $z \sim 9-12$, and to place qualitative constraints on the evolution of the spin temperature of the intergalactic medium (IGM) at $z \geq 9$. In particular, given a detection of ionized bubbles using this technique, it is possible to determine if the IGM surrounding them is typically in absorption or emission. Determining the globally-averaged neutral fraction of the IGM using this method will prove more difficult due to degeneracy with the average size of HII regions., Comment: 14 pages, 11 figures, 2 tables, published in MNRAS. Updated to match published version. Additional results and comments added from previous version. All other results and conclusions remain unchanged

Massive quiescent galaxies (MQGs) are thought to have formed stars rapidly at early times followed by a long period of quiescence. The recent discovery of a MQG, ZF-COSMOS-20115 at $z\sim4$, only 1.5 Gyr after the big bang, places new constraints on galaxy growth and the role of feedback in early star formation. Spectroscopic follow-up confirmed ZF-COSMOS-20115 as a MQG at $z=3.717$ with an estimated stellar mass of ${\sim}10^{11}\mathrm{M}_\odot$, showing no evidence of recent star formation. We use the Meraxes semi-analytic model to investigate how ZF-COSMOS-20115 analogues build stellar mass, and why they become quiescent. We identify three analogue galaxies with similar properties to ZF-COSMOS-20115. We find that ZF-COSMOS-20115 is likely hosted by a massive halo with virial mass of ${\sim}10^{13}\mathrm{M}_\odot$, having been through significant mergers at early times. These merger events drove intense growth of the nucleus, which later prevented cooling and quenched star formation. ZF-COSMOS-20115 likely remained quiescent at $z<3.7$. We find that the analogues host the most massive black holes in our simulation and were luminous quasars at $z\sim5$, indicating that ZF-COSMOS-20115 and other MQGs may be the descendants of high-redshift quasars. In addition, the model suggests that ZF-COSMOS-20115 formed in a region of intergalactic medium that was reionized early., Comment: 11 pages, 7 figures; Updated to match the published version. All results and conclusions remain unchanged

We investigate the clustering properties of Lyman-break galaxies (LBGs) at $z\sim6$ - $8$. Using the semi-analytical model {\scshape Meraxes} constructed as part of the Dark-ages Reionization And Galaxy-formation Observables from Numerical Simulation (DRAGONS) project, we predict the angular correlation function (ACF) of LBGs at $z\sim6$ - $8$. Overall, we find that the predicted ACFs are in good agreement with recent measurements at $z\sim 6$ and $z\sim 7.2$ from observations consisting of the Hubble eXtreme Deep Field (XDF), the Hubble Ultra-Deep Field (HUDF) and Cosmic Assembly Near-infrared Deep Extragalactic Legacy Survey (CANDELS) field. We confirm the dependence of clustering on luminosity, with more massive dark matter haloes hosting brighter galaxies, remains valid at high redshift. The predicted galaxy bias at fixed luminosity is found to increase with redshift, in agreement with observations. We find that LBGs of magnitude $M_{{\rm AB(1600)}} < -19.4$ at $6\lesssim z \lesssim 8$ reside in dark matter haloes of mean mass $\sim 10^{11.0}$- $10^{11.5} M_{\rm \odot}$, and this dark matter halo mass does not evolve significantly during reionisation., Comment: 14 pages, 10 figures, published in MNRAS

Motivated by recent measurements of the number density of faint AGN at high redshift, we investigate the contribution of quasars to reionization by tracking the growth of central supermassive black holes in an update of the Meraxes semi-analytic model. The model is calibrated against the observed stellar mass function at $z\sim0.6-7$, the black hole mass function at $z\lesssim0.5$, the global ionizing emissivity at $z\sim2-5$ and the Thomson scattering optical depth. The model reproduces a Magorrian relation in agreement with observations at $z<0.5$ and predicts a decreasing black hole mass towards higher redshifts at fixed total stellar mass. With the implementation of an opening angle of 80 deg for quasar radiation, corresponding to an observable fraction of ${\sim}23.4$ per cent due to obscuration by dust, the model is able to reproduce the observed quasar luminosity function at $z\sim0.6-6$. The stellar light from galaxies hosting faint AGN contributes a significant or dominant fraction of the UV flux. At high redshift, the model is consistent with the bright end quasar luminosity function and suggests that the recent faint $z\sim4$ AGN sample compiled by Giallongo et al. (2015) includes a significant fraction of stellar light. Direct application of this luminosity function to the calculation of AGN ionizing emissivity consequently overestimates the number of ionizing photons produced by quasars by a factor of 3 at $z\sim6$. We conclude that quasars are unlikely to make a significant contribution to reionization., Comment: 21 pages, 12 figures; Updated to match the published version. All results and conclusions remain unchanged

We investigate how the hydrostatic suppression of baryonic accretion affects the growth rate of dark matter halos during the Epoch of Reionization. By comparing halo properties in a simplistic hydrodynamic simulation in which gas only cools adiabatically, with its collisionless equivalent, we find that halo growth is slowed as hydrostatic forces prevent gas from collapsing. In our simulations, at the high redshifts relevant for reionization (between ${\sim}6$ and ${\sim}11$), halos that host dwarf galaxies ($\lesssim 10^{9} \mathrm{M_\odot}$) can be reduced by up to a factor of 2 in mass due to the hydrostatic pressure of baryons. Consequently, the inclusion of baryonic effects reduces the amplitude of the low mass tail of the halo mass function by factors of 2 to 4. In addition, we find that the fraction of baryons in dark matter halos hosting dwarf galaxies at high redshift never exceeds ${\sim}90\%$ of the cosmic baryon fraction. When implementing baryonic processes, including cooling, star formation, supernova feedback and reionization, the suppression effects become more significant with further reductions of ${\sim}30\%$ to 60\%. Although convergence tests suggest that the suppression may become weaker in higher resolution simulations, this suppressed growth will be important for semi-analytic models of galaxy formation, in which the halo mass inherited from an underlying N-body simulation directly determines galaxy properties. Based on the adiabatic simulation, we provide tables to account for these effects in N-body simulations, and present a modification of the halo mass function along with explanatory analytic calculations., Comment: 17 pages, 11 figures; Updated to match the published version. Two changes in Figures 1 and 3 in order to 1) correct bin sizes of the 10^8 and 10^8.5 Msol bins for NOSN_NOZCOOL_NoRe (was 0.5, should be 0.25); 2) include stellar mass in baryon fraction (was missed in Fig. 3). Quantitative description of Fig. 3 changed slightly in Section 2.2. All other results and conclusions remain unchanged