Search

The Sagittarius Dwarf Spheroidal galaxy (Sgr) is investigated as a target for DM annihilation searches utilizing 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}\ \text{M}^2_\odot\ \text{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 (Crocker, Macias et al. 2022; arXiv:2204.12054) to DM annihilation. We also calculate a J-factor value using an 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: 13 pages, 20 figures. Submitted to MNRAS. Comments welcome

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.

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., 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

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., 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

We present the Dark-ages Reionization and Galaxy-formation Observables from Numerical Simulations (DRAGONS) program and Tiamat, the collisionless N-body simulation program upon which DRAGONS is built. The primary trait distinguishing Tiamat from other large simulation programs is its density of outputs at high redshift (100 from z=35 to z=5; roughly one every 10 Myr) enabling the construction of very accurate merger trees at an epoch when galaxy formation is rapid and mergers extremely frequent. We find that the friends-of-friends halo mass function agrees well with the prediction of Watson et al. at high masses, but deviates at low masses, perhaps due to our use of a different halo finder or perhaps indicating a break from "universal" behaviour. We then analyse the dynamical evolution of galaxies during the Epoch of Reionization finding that only a small fraction (~20%) of galactic halos are relaxed. We illustrate this using standard relaxation metrics to establish two dynamical recovery time-scales: i) halos need ~1.5 dynamical times following formation, and ii) ~2 dynamical times following a major (3:1) or minor (10:1) merger to be relaxed. This is remarkably consistent across a wide mass range. Lastly, we use a phase-space halo finder to illustrate that major mergers drive long-lived massive phase-space structures which take many dynamical times to dissipate. This can yield significant differences in the inferred mass build-up of galactic halos and we suggest that care must be taken to ensure a physically meaningful match between the galaxy-formation physics of semi-analytic models and the halo finders supplying their input., 15 pages, 10 figures. Accepted for publication by MNRAS

We use the Dark-ages, Reionization And Galaxy-formation Observables from Numerical Simulations (DRAGONS) framework to investigate the effect of galaxy-formation physics on the morphology and statistics of ionized hydrogen (HII) regions during the Epoch of Reioinization (EoR). DRAGONS self-consistently couples a semi-analytic galaxy-formation model with the inhomogeneous ionizing UV background, and can therefore be used to study the dependence of morphology and statistics of reionization on feedback phenomena of the ionizing source galaxy population. Changes in galaxy-formation physics modify the sizes of HII regions and the amplitude and shape of 21-cm power spectra. Of the galaxy physics investigated, we find that supernova feedback plays the most important role in reionization, with HII regions up to $\approx 20$ per cent smaller and a fractional difference in the amplitude of power spectra of up to $\approx 17$ per cent at fixed ionized fraction in the absence of this feedback. We compare our galaxy-formation-based reionization models with past calculations that assume constant stellar-to-halo mass ratios and find that with the correct choice of minimum halo mass, such models can mimic the predicted reionization morphology. Reionization morphology at fixed neutral fraction is therefore not uniquely determined by the details of galaxy formation, but is sensitive to the mass of the haloes hosting the bulk of the ionizing sources. Simple EoR parametrizations are therefore accurate predictors of reionization statistics. However, a complete understanding of reionization using future 21-cm observations will require interpretation with realistic galaxy-formation models, in combination with other observations., 14 pages, 12 figures, published in MNRAS

In this paper we present calculations of the UV luminosity function from the Dark-ages Reionization And Galaxy-formation Observables from Numerical Simulations (DRAGONS) project, which combines N-body, semi-analytic and semi-numerical modelling designed to study galaxy formation during the Epoch of Reionization. Using galaxy formation physics including supernova feedback, the model naturally reproduces the UV LFs for high-redshift star-forming galaxies from $z{\sim}5$ through to $z{\sim}10$. We investigate the luminosity--star formation rate (SFR) relation, finding that variable SFR histories of galaxies result in a scatter around the median relation of $0.1$--$0.3$ dex depending on UV luminosity. We find close agreement between the model and observationally derived SFR functions. We use our calculated luminosities to investigate the luminosity function below current detection limits, and the ionizing photon budget for reionization. We predict that the slope of the UV LF remains steep below current detection limits and becomes flat at $M_\mathrm{UV}{\gtrsim}{-14}$. We find that $48$ ($17$) per cent of the total UV flux at $z{\sim}6$ ($10$) has been detected above an observational limit of $M_\mathrm{UV}{\sim}{-17}$, and that galaxies fainter than $M_\mathrm{UV}{\sim}{-17}$ are the main source of ionizing photons for reionization. We investigate the luminosity--stellar mass relation, and find a correlation for galaxies with $M_\mathrm{UV}{, 17 pages, 14 figures, Accepted for publication in MNRAS

We use high resolution N-Body simulations to study the concentration and spin parameters of dark matter haloes in the mass range $10^8\, {\rm M}_{\odot}\, h^{-1} < {\rm M} < 10^{11}\, {\rm M}_{\odot}\, h^{-1}$ and redshifts $5{}5$ concentration$-$mass ($c(M)$) relation that is almost flat and well described by a simple power-law for both NFW and Einasto fits. The intrinsic scatter around the mean relation is $��c_{\rm{vir}}{\sim1}$ (or 20 per cent) at $z=5$. We also find that the analytic model proposed by Ludlow et al. reproduces the mass and redshift-dependence of halo concentrations. Our best-fit Einasto shape parameter, $��$, depends on peak height, $��$, in a manner that is accurately described by $��{=}0.0070��^2{+}0.1839$. The distribution of the spin parameter, $��$, has a weak dependence on equilibrium state; $��$ peaks at roughly ${\sim}0.033$ for our relaxed sample, and at ${\sim}0.04$ for the full population. The spin--virial mass relation has a mild negative correlation at high redshift.

Galaxy redshift surveys using optical telescopes have, in combination with other cosmological probes, enabled precision measurements of the nature of dark energy. We show that radio telescopes are rapidly becoming competitive with optical facilities in spectroscopic surveys of large numbers of galaxies. Two breakthroughs are driving this change. Firstly, individual radio telescopes are more efficient at mapping the sky thanks to the large field-of-view of new phased-array feeds. Secondly, ever more dishes can be correlated in a cost-effective manner with rapid increases in computing power. The next decade will see the coming of age of the 21cm radio wavelength as a cosmological probe as first the Pathfinders then, ultimately, the Square Kilometre Array is constructed. The latter will determine precise 3D positions for a billion galaxies, mapping the distribution of matter in the Universe over the last 12 billion years. This radio telescope will be able to constrain the equation of state of dark energy, and its potential evolution, to a precision rivalling that of future optical facilities such as DESI and Euclid., Comment: 17 pages, 2 figures, 1 table. Accepted to Annalen der Physik for the Special Issue "The Accelerating Universe"

The Australian Square Kilometer Array Pathfinder (ASKAP) will revolutionise our knowledge of gas-rich galaxies in the Universe. Here we present predictions for two proposed extragalactic ASKAP neutral hydrogen (HI) emission-line surveys, based on semi-analytic models applied to cosmological N-body simulations. The ASKAP HI All-Sky Survey, known as WALLABY, is a shallow 3 Pi survey (z = 0 - 0.26) which will probe the mass and dynamics of over 600,000 galaxies. A much deeper small-area HI survey, called DINGO, aims to trace the evolution of HI from z = 0 - 0.43, a cosmological volume of 40 million Mpc^3, detecting potentially 100,000 galaxies. The high-sensitivity 30 antenna ASKAP core (diameter ~2 km) will provide an angular resolution of 30 arcsec (at z=0). Our simulations show that the majority of galaxies detected in WALLABY (87.5%) will be resolved. About 5000 galaxies will be well resolved, i.e. more than five beams (2.5 arcmin) across the major axis, enabling kinematic studies of their gaseous disks. This number would rise to 160,000 galaxies if all 36 ASKAP antennas could be used; the additional six antennas provide baselines up to 6 km, resulting in an angular resolution of 10 arcsec. For DINGO this increased resolution is highly desirable to minimise source confusion; reducing confusion rates from a maximum of 10% of sources at the survey edge to 3%. We estimate that the sources detected by WALLABY and DINGO will span four orders of magnitude in total halo mass (from 10^{11} to 10^{15} Msol) and nearly seven orders of magnitude in stellar mass (from 10^{5} to 10^{12} Msol), allowing us to investigate the process of galaxy formation across the last four billion years., Comment: 21 pages, accepted for publication in MNRAS, minor updates to published version and fixed links. Movies and images available at http://ict.icrar.org/store/Movies/Duffy12c/

The characterisation of the atomic and molecular hydrogen content of high-redshift galaxies is a major observational challenge that will be addressed over the coming years with a new generation of radio telescopes. We investigate this important issue by considering the states of hydrogen across a range of structures within high-resolution cosmological hydrodynamical simulations. Additionally, our simulations allow us to investigate the sensitivity of our results to numerical resolution and to sub-grid baryonic physics (especially feedback from supernovae and active galactic nuclei). We find that the most significant uncertainty in modelling the neutral hydrogen distribution arises from our need to model a self-shielding correction in moderate density regions. Future simulations incorporating radiative transfer schemes will be vital to improve on our empirical self-shielding threshold. Irrespective of the exact nature of the threshold we find that while the atomic hydrogen mass function evolves only mildly from redshift two to zero, the molecular hydrogen mass function increases with increasing redshift, especially at the high-mass end. Interestingly, the weak evolution of the neutral hydrogen mass function is insensitive to the feedback scheme utilised, but the opposite is true for the molecular gas, which is more closely associated with the star formation in the simulations., Comment: 22 pages, 12 figures (4 figures in appendix). Accepted for publication in MNRAS, minor text changes

The back-reaction of baryons on the dark matter halo density profile is of great interest, not least because it is an important systematic uncertainty when attempting to detect the dark matter. Here, we draw on a large suite of high resolution cosmological hydrodynamical simulations, to systematically investigate this process and its dependence on the baryonic physics associated with galaxy formation. The inclusion of baryons results in significantly more concentrated density profiles if radiative cooling is efficient and feedback is weak. The dark matter halo concentration can in that case increase by as much as 30 (10) per cent on galaxy (cluster) scales. The most significant effects occur in galaxies at high redshift, where there is a strong anti-correlation between the baryon fraction in the halo centre and the inner slope of both the total and the dark matter density profiles. If feedback is weak, isothermal inner profiles form, in agreement with observations of massive, early-type galaxies. However, we find that AGN feedback, or extremely efficient feedback from massive stars, is necessary to match observed stellar fractions in groups and clusters, as well as to keep the maximum circular velocity similar to the virial velocity as observed for disk galaxies. These strong feedback models reduce the baryon fraction in galaxies by a factor of 3 relative to the case with no feedback. The AGN is even capable of reducing the baryon fraction by a factor of 2 in the inner region of group and cluster haloes. This in turn results in inner density profiles which are typically shallower than isothermal and the halo concentrations tend to be lower than in the absence of baryons., Comment: 20 pages, 14 figures, 1 table. MNRAS in press. Version 2: added a few references.

We discuss the possibility of performing a substantial spectroscopic galaxy redshift survey selected via the 21cm emission from neutral hydrogen using the Five-hundred metre Aperture Spherical Telescope (FAST) to be built in China. We consider issues related to the estimation of the source counts and optimizations of the survey, and discuss the constraints on cosmological models that such a survey could provide. We find that a survey taking around two years could detect ~10^7 galaxies with an average redshift of ~0.15 making the survey complementary to those already carried out at optical wavelengths. These conservative estimates have used the z=0 HI mass function and have ignored the possibility of evolution. The results could be used to constrain Gamma = (Omega_m h) to 5 per cent and the spectral index, n_s, to 7 per cent independent of cosmic microwave background data. If we also use simulated power spectra from the Planck satellite, we can constrain w to be within 5 per cent of -1., Comment: 13 pages, 15 figures, Accepted by MNRAS, minor corrections