Search

This work considers the impact of resolution in the construction of mock observations of simulated galaxies. In particular, when building mock integral field spectroscopic observations from galaxy formation models in cosmological simulations, we investigate the possible systematics that may arise given the assumption that all galaxies above some stellar mass limit will provide unbiased and meaningful observable stellar kinematics. We build a catalogue of N-body simulations to sample the range of stellar particle resolutions within the EAGLE Ref0050N0752 simulation box and examine how their observable kinematics vary relative to a higher-resolution N-body control. We use these models to compile a table of the minimum number of particles-per-pixel to reach a given uncertainty in the fitted line-of-sight velocity distribution parameters. Further, we introduce a Voronoi-binning module to the mock observation code, SimSpin, in order to meet these minimum numbers. Using EAGLE, we show the impact of this shot noise on the observed spin-ellipticity plane and the recovery of this space when observations are binned with increasing numbers of particles. In conclusion, we advise binning mock images to meet at least 200 particles-per-pixel to avoid systematically under-estimating the velocity dispersion along a given line-of-sight. We demonstrate that this is important for comparing galaxies extracted from the same simulation, as well as between simulations of varying mass resolution and observations of real galaxies., Comment: 19 pages, 11 figures (including 3 pages of appendix with 2 figures). Accepted 2024 November 3 for publication in MNRAS. Received 2024 October 28; in original form 2024 August 9

In this work, we present a methodology and a corresponding code-base for constructing mock integral field spectrograph (IFS) observations of simulated galaxies in a consistent and reproducible way. Such methods are necessary to improve the collaboration and comparison of observation and theory results, and accelerate our understanding of how the kinematics of galaxies evolve over time. This code, SimSpin, is an open-source package written in R, but also with an API interface such that the code can be interacted with in any coding language. Documentation and individual examples can be found at the open-source website connected to the online repository. SimSpin is already being utilised by international IFS collaborations, including SAMI and MAGPI, for generating comparable data sets from a diverse suite of cosmological hydrodynamical simulations., Comment: 20 pages, 15 figures, 2 tables. Accepted for publication in PASA. 30/08/23

We describe the scientific goals and survey design of the First Large Absorption Survey in HI (FLASH), a wide field survey for 21-cm line absorption in neutral atomic hydrogen (HI) at intermediate cosmological redshifts. FLASH will be carried out with the Australian Square Kilometre Array Pathfinder (ASKAP) radio telescope and is planned to cover the sky south of $\delta \approx +40$deg at frequencies between 711.5 and 999.5MHz. At redshifts between $z = 0.4$ and $1.0$ (look back times of 4 - 8Gyr), the HI content of the Universe has been poorly explored due to the difficulty of carrying out radio surveys for faint 21-cm line emission and, at ultra-violet wavelengths, space-borne searches for Damped Lyman-$\alpha$ absorption in quasar spectra. The ASKAP wide field of view and large spectral bandwidth, in combination with a radio-quiet site, will enable a search for absorption lines in the radio spectra of bright continuum sources over 80% of the sky. This survey is expected to detect at least several hundred intervening 21-cm absorbers, and will produce an HI-absorption-selected catalogue of galaxies rich in cool, star-forming gas, some of which may be concealed from optical surveys. Likewise, at least several hundred associated 21-cm absorbers are expected to be detected within the host galaxies of radio sources at $0.4 < z < 1.0$, providing valuable kinematical information for models of gas accretion and jet-driven feedback in radio-loud active galactic nuclei. FLASH will also detect OH 18-cm absorbers in diffuse molecular gas, megamaser OH emission, radio recombination lines, and stacked HI emission., Comment: 31 pages, 10 figures, 6 tables, accepted for publication in PASA

Observers experience a series of limitations when measuring galaxy kinematics, such as variable seeing conditions and aperture size. These effects can be reduced using empirical corrections, but these equations are usually applicable within a restrictive set of boundary conditions (e.g. S\'ersic indices within a given range) which can lead to biases when trying to compare measurements made across a full kinematic survey. In this work, we present new corrections for two widely used kinematic parameters, $\lambda_R$ and $V/\sigma$, that are applicable across a broad range of galaxy shapes, measurement radii and ellipticities. We take a series of mock observations of N-body galaxy models and use these to quantify the relationship between the observed kinematic parameters, structural properties and different seeing conditions. Derived corrections are then tested using the full catalogue of galaxies, including hydro-dynamic models from the EAGLE simulation. Our correction is most effective for regularly-rotating systems, yet the kinematic parameters of all galaxies -- fast, slow and irregularly rotating systems -- are recovered successfully. We find that $\lambda_R$ is more easily corrected than $V/\sigma$, with relative deviations of 0.02 and 0.06 dex respectively. The relationship between $\lambda_R$ and $V/\sigma$, as described by the parameter $\kappa$, also has a minor dependence on seeing conditions. These corrections will be particularly useful for stellar kinematic measurements in current and future integral field spectroscopic (IFS) surveys of galaxies., Comment: Accepted for publication in MNRAS 2020 June 20. 29 pages, 26 figures, including further supplementary material included online in Appendix D

We present SimSpin, a new, public, software framework for generating integral field spectroscopy (IFS) data cubes from N-body/hydrodynamical simulations of galaxies, which can be compared directly with observational datasets. SimSpin provides a consistent method for studying a galaxy's stellar component. It can be used to explore how observationally inferred measurements of kinematics, such as the spin parameter $\lambda_R$, are impacted by the effects of, for example, inclination, seeing conditions, distance, etc. SimSpin is written in R and has been designed to be highly modular, flexible, and extensible. It is already being used by the astrophysics community to generate IFS-like cubes and FITS files for direct comparison of simulations to observations. In this paper, we explain the conceptual framework of SimSpin; how it is implemented in R; and we demonstrate SimSpin's current capabilities, providing as an example a brief investigation of how numerical resolution affects how reliably we can recover the intrinsic stellar kinematics of a simulated galaxy., Comment: 14 pages, 10 figures. Accepted for publication in PASA

We investigate the impact of time-resolved `gradual' stellar feedback processes in high redshift dwarf spheroidal galaxies. Here `gradual' feedback refers to individual stellar feedback events which deposit energy over a period of time. We conduct high resolution hydrodynamical simulations of dwarf spheroidal galaxies with halo masses of 10$^7$ M$_{\odot}$ - 10$^8$ M$_{\odot}$, based on z = 6 progenitors of the Milky Way's dwarf spheroidal galaxies. We also include a novel feedback prescription for individual massive stars, which includes stellar winds and a HMXB (High Mass X-ray Binary) phase, on top of supernovae. We find the mass of gas unbound across a 1 Gyr starburst is uniformly lowered if gradual feedback mechanisms are included across the range of metallicities, halo concentration parameters and galaxy masses studied here. Furthermore, we find including gradual feedback in the smallest galaxies delays the unbinding of the majority of the gas and facilitates the production of `chimneys' in the dense shell surrounding the feedback generated hot, pressurised `superbubble'. These `chimneys' vent hot gas from the galaxy interior, lowering the temperature of the central 10 kpc of the gaseous halo. Additionally, we find radiative cooling has little effect on the energetics of simulations that include a short, violent starburst compared with those that have a longer, less concentrated starburst. Finally, we investigate the relative impact of HMXB feedback and stellar winds on our results, finding the ubiquity of stellar winds throughout each starburst makes them a defining factor in the final state of the interstellar medium., Comment: accepted for publication in MNRAS, 18 pages, 35 figures in main body of text (bibliography and appendices extra)

Accurately predicting the shape of the HI velocity function of galaxies is regarded widely as a fundamental test of any viable dark matter model. Straightforward analyses of cosmological $N$-body simulations imply that the $\Lambda$CDM model predicts an overabundance of low circular velocity galaxies when compared to observed HI velocity functions. More nuanced analyses that account for the relationship between galaxies and their host haloes suggest that how we model the influence of baryonic processes has a significant impact on HI velocity function predictions. We explore this in detail by modelling HI emission lines of galaxies in the SHARK semi-analytic galaxy formation model, built on the SURFS suite of $\Lambda$CDM $N$-body simulations. We create a simulated ALFALFA survey, in which we apply the survey selection function and account for effects such as beam confusion, and compare simulated and observed HI velocity width distributions, finding differences of $\lesssim 50$%, orders of magnitude smaller than the discrepancies reported in the past. This is a direct consequence of our careful treatment of survey selection effects and, importantly, how we model the relationship between galaxy and halo circular velocity - the HI mass-maximum circular velocity relation of galaxies is characterised by a large scatter. These biases are complex enough that building a velocity function from the observed HI line widths cannot be done reliably., Comment: accepted to MNRAS 19 pages, 14 figures

WAVES is designed to study the growth of structure, mass and energy on scales of ~1 kpc to ~10 Mpc over a 7 Gyr timeline. On the largest length scales (1-10 Mpc) WAVES will measure the structures defined by groups, filaments and voids, and their emergence over recent times. Comparisons with bespoke numerical simulations will be used to confirm, refine or refute the Cold Dark Matter paradigm. At intermediate length scales (10 kpc-1 Mpc) WAVES will probe the size and mass distribution of galaxy groups, as well as the galaxy merger rates, in order to directly measure the assembly of dark matter halos and stellar mass. On the smallest length scales (1-10 kpc) WAVES will provide accurate distance and environmental measurements to complement high-resolution space-based imaging to study the mass and size evolution of galaxy bulges, discs and bars. In total, WAVES will provide a panchromatic legacy dataset of ~1.6 million galaxies, firmly linking the very low ($z < 0.1$) and intermediate ($z \sim 0.8$) redshift Universe., Comment: Part of the 4MOST issue of The Messenger, published in preparation of the 4MOST Community Workshop, see http://www.eso.org/sci/meetings/2019/4MOST.html

We introduce the 4-metre Multi-Object Spectroscopic Telescope (4MOST), a new high-multiplex, wide-field spectroscopic survey facility under development for the four-metre-class Visible and Infrared Survey Telescope for Astronomy (VISTA) at Paranal. Its key specifications are: a large field of view (FoV) of 4.2 square degrees and a high multiplex capability, with 1624 fibres feeding two low-resolution spectrographs ($R = \lambda/\Delta\lambda \sim 6500$), and 812 fibres transferring light to the high-resolution spectrograph ($R \sim 20\,000$). After a description of the instrument and its expected performance, a short overview is given of its operational scheme and planned 4MOST Consortium science; these aspects are covered in more detail in other articles in this edition of The Messenger. Finally, the processes, schedules, and policies concerning the selection of ESO Community Surveys are presented, commencing with a singular opportunity to submit Letters of Intent for Public Surveys during the first five years of 4MOST operations., Comment: Part of the 4MOST issue of The Messenger, published in preparation of 4MOST Community Workshop, see http://www.eso.org/sci/meetings/2019/4MOST.html

A primary goal of integral field spectroscopic (IFS) surveys is to provide a statistical census of galaxies classified by their internal kinematics. As a result, the observational spin parameter, $\lambda_R$, has become one of the most popular methods of quantifying the relative importance of velocity dispersion and rotation in supporting a galaxy's inner structure. The goal of this paper is to examine the relationship between the observationally deduced $\lambda_R$ and one of the most commonly used theoretical spin parameters in the literature, the Bullock et al. (2001) $\lambda'$. Using a set of $N$-body realisations of galaxies from which we construct mock IFS observations, we measure $\lambda_R$ as an observer would, incorporating the effects of beam smearing and seeing conditions. Assuming parameters typical of current IFS surveys, we confirm that there are strong positive correlations between $\lambda_R$ and measurement radius, and strong negative correlations between $\lambda_R$ and size of the PSF, for late-type galaxies; these biases can be reduced using a recently proposed empirical correction. Once observational biases are corrected for, we find that $\lambda_R$ provides a good approximation to $\sim \sqrt{3}/2 \; \lambda'(\rm R_{\rm eff})$, where $\lambda'$ is evaluated for the galactic stellar component within 1 R$_{\rm eff}$., Comment: Accepted for publication in MNRAS. 15 pages, 14 figures

Galaxy cluster outskirts mark the transition region from the mildly non-linear cosmic web to the highly non-linear, virialised, cluster interior. It is in this transition region that the intra-cluster medium (ICM) begins to influence the properties of accreting galaxies and groups, as ram pressure impacts a galaxy's cold gas content and subsequent star formation rate. Conversely, the thermodynamical properties of the ICM in this transition region should also feel the influence of accreting substructure (i.e. galaxies and groups), whose passage can drive shocks. In this paper, we use a suite of cosmological hydrodynamical zoom simulations of a single galaxy cluster, drawn from the nIFTy comparison project, to study how the dynamics of substructure accreted from the cosmic web influences the thermodynamical properties of the ICM in the cluster's outskirts. We demonstrate how features evident in radial profiles of the ICM (e.g. gas density and temperature) can be linked to strong shocks, transient and short-lived in nature, driven by the passage of substructure. The range of astrophysical codes and galaxy formation models in our comparison are broadly consistent in their predictions (e.g. agreeing when and where shocks occur, but differing in how strong shocks will be); this is as we would expect of a process driven by large-scale gravitational dynamics and strong, inefficiently radiating, shocks. This suggests that mapping such shock structures in the ICM in a cluster's outskirts (via e.g. radio synchrotron emission) could provide a complementary measure of its recent merger and accretion history., Comment: 15 pages, 13 figures, version accepted for publication in MNRAS

We investigate the effects of gradual heating on the evolution of turbulent molecular clouds of mass $2\times 10^6$ M$_\odot$ and virial parameters ranging between $0.7-1.2$. This gradual heating represents the energy output from processes such as winds from massive stars or feedback from High Mass X-ray binaries (HMXBs), contrasting the impulsive energy injection from supernovae (SNe). For stars with a mass high enough that their lifetime is shorter than the life of the cloud, we include a SN feedback prescription. Including both effects, we investigate the interplay between slow and fast forms of feedback and their effectiveness at triggering/suppressing star formation. We find that SN feedback can carve low density chimneys in the gas, offering a path of least resistance for the energy to escape. Once this occurs the more stable, but less energetic, gradual feedback is able to keep the chimneys open. By funneling the hot destructive gas away from the centre of the cloud, chimneys can have a positive effect on both the efficiency and duration of star formation. Moreover, the critical factor is the number of high mass stars and SNe (and any subsequent HMXBs) active within the free-fall time of each cloud. This can vary from cloud to cloud due to the stochasticity of SN delay times and in HMXB formation. However, the defining factor in our simulations is the efficiency of the cooling, which can alter the Jeans mass required for sink particle formation, along with the number of massive stars in the cloud., Comment: 35 pages, 46 figures, accepted for publication in MNRAS

We introduce a method for modeling disk galaxies designed to take full advantage of data from integral field spectroscopy (IFS). The method fits equilibrium models to simultaneously reproduce the surface brightness, rotation and velocity dispersion profiles of a galaxy. The models are fully self-consistent 6D distribution functions for a galaxy with a Sersic-profile stellar bulge, exponential disk and parametric dark matter halo, generated by an updated version of GalactICS. By creating realistic flux-weighted maps of the kinematic moments (flux, mean velocity and dispersion), we simultaneously fit photometric and spectroscopic data using both maximum-likelihood and Bayesian (MCMC) techniques. We apply the method to a GAMA spiral galaxy (G79635) with kinematics from the SAMI Galaxy Survey and deep $g$- and $r$-band photometry from the VST-KiDS survey, comparing parameter constraints with those from traditional 2D bulge-disk decomposition. Our method returns broadly consistent results for shared parameters, while constraining the mass-to-light ratios of stellar components and reproducing the HI-inferred circular velocity well beyond the limits of the SAMI data. While the method is tailored for fitting integral field kinematic data, it can use other dynamical constraints like central fibre dispersions and \HI circular velocities, and is well-suited for modelling galaxies with a combination of deep imaging and HI and/or optical spectra (resolved or otherwise). Our implementation (MagRite) is computationally efficient and can generate well-resolved models and kinematic maps in under a minute on modern processors., Comment: ApJ accepted; 21 pages, 12 figures

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 ($\sigma$) maps are created and compared to observations of disc galaxies obtained with the Sydney-AAO Multi-object Integral field spectrograph (SAMI), where $\sigma$-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 ($\Sigma_{\rm SFR}$), with low-M$_{\star}/$low-$\Sigma_{\rm SFR}$ galaxies showing a narrow peak at low $\sigma$ ($\sim30$ km s$^{-1}$) and more active, high-M$_{\star}/$high-$\Sigma_{\rm SFR}$ galaxies reaching $\sigma>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 $\sigma$-distribution. The following equivalence relations hold in EAGLE: $i)$ low-$\sigma$ peak $\,\Leftrightarrow\,$ disc of the galaxy $\,\Leftrightarrow\,$ gas with T $<10^5$ K; $ii)$ high-$\sigma$ tail $\,\Leftrightarrow\,$ galactic winds $\,\Leftrightarrow\,$ gas with T $\ge 10^5$ K., Comment: 18 pages, 12 figures, accepted for publication in MNRAS

The formation and evolution of galaxies with low neutral atomic hydrogen (HI) masses, M$_{\rm HI}$$<$10$^{8}h^{-2}$M$_{\odot}$, are affected by host dark matter halo mass and photoionisation feedback from the UV background after the end of reionization. We study how the physical processes governing the formation of galaxies with low HI mass are imprinted on the distribution of neutral hydrogen in the Universe using the hierarchical galaxy formation model, GALFORM. We calculate the effect on the correlation function of changing the HI mass detection threshold at redshifts $0 \le z \le 0.5$. We parameterize the clustering as $\xi(r)=(r/r_{0})^{-\gamma}$ and we find that including galaxies with M$_{\rm HI}$$<$10$^{8}h^{-2}$M$_{\odot}$ increases the clustering amplitude $r_{0}$ and slope $\gamma$ compared to samples of higher HI masses. This is due to these galaxies with low HI masses typically being hosted by haloes with masses greater than 10$^{12}{h}^{-1}$M$_{\odot}$, and is in contrast to optically selected surveys for which the inclusion of faint, blue galaxies lowers the clustering amplitude. We show the HI mass function for different host dark matter halo masses and galaxy types (central or satellite) to interpret the values of $r_{0}$ and $\gamma$ of the clustering of HI-selected galaxies. We also predict the contribution of low HI mass galaxies to the 21cm intensity mapping signal. We calculate that a dark matter halo mass resolution better than $\sim$10$^{10}{h}^{-1}$M$_{\odot}$ at redshifts higher than 0.5 is required in order to predict converged 21cm brightness temperature fluctuations., Comment: 14 pages, 10 figures, accepted for pubication in MNRAS

We explore the galaxy formation physics governing the low mass end of the HI mass function in the local Universe. Specifically, we predict the effects on the HI mass function of varying i) the strength of photoionisation feedback and the redshift of the end of the epoch of reionization, ii) the cosmology, iii) the supernovae feedback prescription, and iv) the efficiency of star formation. We find that the shape of the low-mass end of the HI mass function is most affected by the critical halo mass below which galaxy formation is suppressed by photoionisation heating of the intergalactic medium. We model the redshift dependence of this critical dark matter halo mass by requiring a match to the low-mass end of the HI mass function. The best fitting critical dark matter halo mass decreases as redshift increases in this model, corresponding to a circular velocity of $\sim 50 \, {\rm km \,s}^{-1}$ at $z=0$, $\sim 30 \, {\rm km\, s}^{-1}$ at $z \sim 1$ and $\sim 12 \, {\rm km \, s}^{-1}$ at $z=6$. We find that an evolving critical halo mass is required to explain both the shape and abundance of galaxies in the HI mass function below $M_{\rm HI} \sim 10^{8} h^{-2} {\rm M_{\odot}}$. The model makes specific predictions for the clustering strength of HI-selected galaxies with HI masses > $10^{6} h^{-2} {\rm M_{\odot}}$ and $> 10^{7} h^{-2} {\rm M_{\odot}}$ and for the relation between the HI and stellar mass contents of galaxies which will be testable with upcoming surveys with the Square Kilometre Array and its pathfinders. We conclude that measurements of the HI mass function at $z \ge 0$ will lead to an improvement in our understanding of the net effect of photoionisation feedback on galaxy formation and evolution., Comment: 13 pages, 9 figures, Minor corrections added, revised to match version accepted by MNRAS

Theoretical galaxy formation models are an established and powerful tool for interpreting the astrophysical significance of observational data, particularly galaxy surveys. Such models have been utilised with great success by optical surveys such as 2dFGRS and SDSS, but their application to radio surveys of cold gas in galaxies has been limited. In this chapter we describe recent developments in the modelling of the cold gas properties in the models, and how these developments are essential if they are to be applied to cold gas surveys of the kind that will be carried out with the SKA. By linking explicitly a galaxy's star formation rate to the abundance of molecular hydrogen in the galaxy rather than cold gas abundance, as was assumed previously, the latest models reproduce naturally many of the global atomic and molecular hydrogen properties of observed galaxies. We review some of the key results of the latest models and highlight areas where further developments are necessary. We discuss also how model predictions can be most accurately compared with observational data, what challenges we expect when creating synthetic galaxy surveys in the SKA era, and how the SKA can be used to test models of dark matter., Comment: 20 pages, 9 figures. Contribution to the conference 'Advancing Astrophysics with the Square Kilometre Array', June 8-13, 2014, Giardini Naxos, Italy

One of the key science drivers for the development of the SKA is to observe the neutral hydrogen, HI, in galaxies as a means to probe galaxy evolution across a range of environments over cosmic time. Over the past decade, much progress has been made in theoretical simulations and observations of HI in galaxies. However, recent HI surveys on both single dish radio telescopes and interferometers, while providing detailed information on global HI properties, the dark matter distribution in galaxies, as well as insight into the relationship between star formation and the interstellar medium, have been limited to the local universe. Ongoing and upcoming HI surveys on SKA pathfinder instruments will extend these measurements beyond the local universe to intermediate redshifts with long observing programmes. We present here an overview of the HI science which will be possible with the increased capabilities of the SKA and which will build upon the expected increase in knowledge of HI in and around galaxies obtained with the SKA pathfinder surveys. With the SKA1 the greatest improvement over our current measurements is the capability to image galaxies at reasonable linear resolution and good column density sensitivity to much higher redshifts (0.2 < z < 1.7). So one will not only be able to increase the number of detections to study the evolution of the HI mass function, but also have the sensitivity and resolution to study inflows and outflows to and from galaxies and the kinematics of the gas within and around galaxies as a function of environment and cosmic time out to previously unexplored depths. The increased sensitivity of SKA2 will allow us to image Milky Way-size galaxies out to redshifts of z=1 and will provide the data required for a comprehensive picture of the HI content of galaxies back to z~2 when the cosmic star formation rate density was at its peak., Comment: 25 pages, 5 figures, 3 tables. Contribution to the conference 'Advancing Astrophysics with the Square Kilometre Array', June 8-13, 2014, Giardini Naxos, Italy

Angular momentum is one of the most fundamental physical quantities governing galactic evolution. Differences in the colours, morphologies, star formation rates and gas fractions amongst galaxies of equal stellar/baryon mass M are potentially widely explained by variations in their specific stellar/baryon angular momentum j. The enormous potential of angular momentum science is only just being realised, thanks to the emergence of the first simulations of galaxies with converged spins, paralleled by a dramatic increase in kinematic observations. Such observations are still challenged by the fact that most of the stellar/baryon angular momentum resides at large radii. In fact, the radius that maximally contributes to the angular momentum of an exponential disk (3Re-4Re) is twice as large as the radius that maximally contributes to the disk mass; thus converged measurements of angular momentum require either extremely deep IFS data or, alternatively, kinematic measurements of neutral atomic hydrogen (HI), which naturally resides at the large disk radii that dominate the angular momentum. The SKA has a unique opportunity to become the world-leading facility for angular momentum studies due to its ability to measure the resolved and/or global HI kinematics in very large and well-characterised galaxy samples. These measurements will allow, for example, (1) a very robust determination of the two-dimensional distribution of galaxies in the (M,j)-plane, (2) the largest, systematic measurement of the relationship between M, j, and tertiary galaxy properties, and (3) the most accurate measurement of the large-scale distribution and environmental dependence of angular momentum vectors, both in terms of norm and orientation. All these measurements will represent exquisite tools to build a next generation of galaxy evolution models., Comment: 13 pages, 5 figures, 3 tables

We present the Early Data Release of the Sydney-AAO Multi-object Integral field spectrograph (SAMI) Galaxy Survey. The SAMI Galaxy Survey is an ongoing integral field spectroscopic survey of ~3400 low-redshift (z<0.12) galaxies, covering galaxies in the field and in groups within the Galaxy And Mass Assembly (GAMA) survey regions, and a sample of galaxies in clusters. In the Early Data Release, we publicly release the fully calibrated datacubes for a representative selection of 107 galaxies drawn from the GAMA regions, along with information about these galaxies from the GAMA catalogues. All datacubes for the Early Data Release galaxies can be downloaded individually or as a set from the SAMI Galaxy Survey website. In this paper we also assess the quality of the pipeline used to reduce the SAMI data, giving metrics that quantify its performance at all stages in processing the raw data into calibrated datacubes. The pipeline gives excellent results throughout, with typical sky subtraction residuals in the continuum of 0.9-1.2 per cent, a relative flux calibration uncertainty of 4.1 per cent (systematic) plus 4.3 per cent (statistical), and atmospheric dispersion removed with an accuracy of 0."09, less than a fifth of a spaxel., Comment: MNRAS accepted. 16 pages, 20 figures. Galaxy datacubes and related data available from http://sami-survey.org/edr . v2: Minor edits to match accepted version

We study the contribution of galaxies with different properties to the global densities of star formation rate (SFR), atomic (HI) and molecular hydrogen (H2) as a function of redshift. We use the GALFORM model of galaxy formation, which is set in the LCDM framework. This model includes a self-consistent calculation of the SFR, which depends on the H2 content of galaxies. The predicted SFR density and how much of this is contributed by galaxies with different stellar masses and infrared luminosities are in agreement with observations. The model predicts a modest evolution of the HI density at z<3, which is also in agreement with the observations. The HI density is predicted to be always dominated by galaxies with SFR<1Msun/yr. This contrasts with the H2 density, which is predicted to be dominated by galaxies with SFR>10Msun/yr. Current high-redshift galaxy surveys are limited to detect carbon monoxide in galaxies with SFR>30Msun/yr, which in our model make up, at most, 20% of the H2 in the universe. In terms of stellar mass, the predicted H2 density is dominated by massive galaxies, Mstellar>10^10Msun, while the HI density is dominated by low mass galaxies, Mstellar<10^9Msun. In the context of upcoming neutral gas surveys, we suggest that the faint nature of the galaxies dominating the HI content of the Universe will hamper the identification of optical counterparts, while for H2, we expect follow up observations of molecular emission lines of already existing galaxy catalogues to be able to uncover the H2 density of the Universe., Comment: 23 pages (18 of main body and 5 of Appendices), 16 figures. Accepted by MNRAS

Abridged: We simulate a massive galaxy cluster in a LCDM Universe using three different approaches to solving the equations of non-radiative hydrodynamics: `classic' Smoothed Particle Hydrodynamics (SPH); a novel SPH with a higher order dissipation switch (SPHS); and adaptive mesh refinement (AMR). We find that SPHS and AMR are in excellent agreement, with both forming a well-defined entropy core that rapidly converges with increasing mass and force resolution. By contrast, SPH exhibits rather different behaviour. At low redshift, entropy decreases systematically with decreasing cluster-centric radius, converging on ever lower central values with increasing resolution. At higher redshift, SPH is in better agreement with SPHS and AMR but shows much poorer numerical convergence. We trace these discrepancies to artificial surface tension in SPH at phase boundaries. At early times, the passage of massive substructures close to the cluster centre stirs and shocks gas to build an entropy core. At later times, artificial surface tension causes low entropy gas to sink artificially to the centre of the cluster. We use SPHS to study the contribution of numerical versus physical dissipation on the entropy core, and argue that numerical dissipation is required to ensure single-valued fluid quantities in converging flows. However, provided this dissipation occurs only at the resolution limit, and provided that it does not propagate errors to larger scales, its effect is benign. There is no requirement to build `sub-grid' models of unresolved turbulence for galaxy cluster simulations. We conclude that entropy cores in non-radiative simulations of galaxy clusters are physical, resulting from entropy generation in shocked gas during cluster assembly, putting to rest the long-standing puzzle of cluster entropy cores in AMR simulations versus their apparent absence in classic SPH simulations., Comment: 14 pages, 10 figures; submitted to MNRAS

Analysing all Galaxy and Mass Assembly (GAMA) galaxies within a factor two (+/- 0.3 dex) of the stellar mass of the Milky Way (MW), there is a 11.9% chance that one of these galaxies will have a close companion (within a projected separation of 70 kpc and radial separation of 400 km/s) that is at least as massive as the Large Magellanic Cloud (LMC). Two close companions at least as massive as the Small Magellanic Cloud (SMC) are rare at the 3.4% level. Two full analogues to the MW- LMC-SMC system were found in GAMA (all galaxies late-type and star forming), suggesting such a combination of close together, late-type, star-forming galaxies is rare: only 0.4% of MW mass galaxies (in the range where we could observe both the LMC and SMC) have such a system. In summary, the MW-LMC-SMC system is a 2.7? event (when recast into Gaussian statistics). Using cross-correlation comparisons we find that there is a preference for SMC- LMC binary pair analogues to be located within 2 Mpc of a range of different lumi- nosity groups. There is a particular preference is for such binaries to be located near LG luminosity systems. When these groups are subdivided into small magnitude gap and large magnitude gap subsets, the binaries prefer to be spatially associated with the small magnitude gap systems. These systems will be dynamically less evolved, but still offer the same amount of gravitational dark matter. This suggests that binaries such as the SMC-LMC might be transient systems, usually destroyed during vigorous merger events. Details of a particularly striking analogue to the MW-SMC-LMC and M31 complex are included., Comment: 6 pages 2 figures, Published in MNRAS, presented at IAU SpS2 Thu 23rd Aug 11:45am

We investigate the impact of feedback - from supernovae (SNe), active galactic nuclei (AGN) and a photo-ionizing background at high redshifts - on the neutral atomic hydrogen (HI) mass function, the $b_{\rm J}$ band luminosity function, and the spatial clustering of these galaxies at $z$=0. We use a version of the semi-analytical galaxy formation model GALFORM that calculates self-consistently the amount of HI in a galaxy as a function of cosmic time and links its star formation rate to its mass of molecular hydrogen (H$_2$). We find that a systematic increase or decrease in the strength of SNe feedback leads to a systematic decrease or increase in the amplitudes of the luminosity and HI mass functions, but has little influence on their overall shapes. Varying the strength of AGN feedback influences only the numbers of the brightest or most HI massive galaxies, while the impact of varying the strength of photo-ionization feedback is restricted to changing the numbers of the faintest or least HI massive galaxies.Our results suggest that the HI mass function is a more sensitive probe of the consequences of cosmological reionization for galaxy formation than the luminosity function. We find that increasing the strength of any of the modes of feedback acts to weaken the clustering strength of galaxies, regardless of their HI-richness. In contrast, weaker AGN feedback has little effect on the clustering strength whereas weaker SNe feedback increases the clustering strength of HI-poor galaxies more strongly than HI-rich galaxies. These results indicate that forthcoming HI surveys on next generation radio telescopes such as the Square Kilometre Array and its pathfinders will be exploited most fruitfully as part of multiwavelength survey campaigns., Comment: 10 pages, 3 figures, Accepted version by MNRAS

We discuss massive outflows in galaxy bulges, particularly ones driven by accretion episodes where the central supermassive black hole reaches the Eddington limit. We show that the quasar radiation field Compton--cools the wind shock until this reaches distances $\sim 1$ kpc from the black hole, but becomes too dilute to do this at larger radii. Radiative processes cannot cool the shocked gas within the flow time at any radius. Outflows are therefore momentum-driven at small radii (as required to explain the $M - \sigma$ relation). At large radii they are energy-driven, contrary to recent claims. We solve analytically the motion of an energy--driven shell after the central source has turned off. This shows that the thermal energy in the shocked wind can drive further expansion for a time $\sim 10$ times longer than the active time of the central source. Outflows observed at large radii with no active central source probably result from an earlier short (few Myr) active phase of this source., Comment: 5 pages, 2 figures, accepted for publication in MNRAS

Due to their large dynamical mass-to-light ratios, dwarf spheroidal galaxies (dSphs) are promising targets for the indirect detection of dark matter (DM) in gamma-rays. We examine their detectability by present and future gamma-ray observatories. The key innovative features of our analysis are: (i) We take into account the angular size of the dSphs; while nearby objects have higher gamma ray flux, their larger angular extent can make them less attractive targets for background-dominated instruments. (ii) We derive DM profiles and the astrophysical J-factor (which parameterises the expected gamma-ray flux, independently of the choice of DM particle model) for the classical dSphs directly from photometric and kinematic data. We assume very little about the DM profile, modelling this as a smooth split-power law distribution, with and without sub-clumps. (iii) We use a Markov Chain Monte Carlo (MCMC) technique to marginalise over unknown parameters and determine the sensitivity of our derived J-factors to both model and measurement uncertainties. (iv) We use simulated DM profiles to demonstrate that our J-factor determinations recover the correct solution within our quoted uncertainties. Our key findings are: (i) Sub-clumps in the dSphs do not usefully boost the signal; (ii) The sensitivity of atmospheric Cherenkov telescopes to dSphs within 20 kpc with cored halos can be up to ~50 times worse than when estimated assuming them to be point-like. Even for the satellite-borne Fermi-LAT the sensitivity is significantly degraded on the relevant angular scales for long exposures, hence it is vital to consider the angular extent of the dSphs when selecting targets; (iii) No DM profile has been ruled out by current data, but using a prior on the inner dark matter cusp slope 0<=gamma<=1 provides J-factor estimates accurate to a factor of a few if an appropriate angular scale [abridged], Comment: 31 pages, 49 figures (2 new figures, section 5.1.5 extended, new App. H2, clarifications), accepted for publication in MNRAS

We show that star formation in galaxy bulges is self-regulating through momentum feedback, limiting the stellar bulge mass to M_b ~ sigma^4. Together with a black hole mass M_BH ~ sigma^4 set by AGN momentum feedback, this produces a linear M_BH - M_b relation. At low redshift this gives M_BH/M_b ~ 0.001, close to the observed ratio. We show that AGN feedback can remove any remaining gas from the bulge and terminate star formation once the central black hole reaches the M_BH - sigma value, contrary to earlier claims. We find a mild upward deviation from the sigma^4 law at higher redshift and at higher sigma., Comment: 6 pages, no figures, accepted for publication in MNRAS

The distribution of cold gas in dark matter haloes is driven by key processes in galaxy formation: gas cooling, galaxy mergers, star formation and reheating of gas by supernovae. We compare the predictions of four different galaxy formation models for the spatial distribution of cold gas. We find that satellite galaxies make little contribution to the abundance or clustering strength of cold gas selected samples, and are far less important than they are in optically selected samples. The halo occupation distribution function of present-day central galaxies with cold gas mass > 10^9 h^-1 Msun is peaked around a halo mass of ~ 10^11 h^-1 Msun, a scale that is set by the AGN suppression of gas cooling. The model predictions for the projected correlation function are in good agreement with measurements from the HI Parkes All-Sky Survey. We compare the effective volume of possible surveys with the Square Kilometre Array with those expected for a redshift survey in the near-infrared. Future redshift surveys using neutral hydrogen emission will be competitive with the most ambitious spectroscopic surveys planned in the near-infrared., Comment: 19 pages, 15 figures, To appear in MNRAS

(Abridged) Accurately predicting how the cosmic abundance of neutral hydrogen evolves with redshift is a challenging problem facing modellers of galaxy formation. We investigate the predictions of four currently favoured semi-analytical galaxy formation models applied to the Millennium simulation for the mass function of cold neutral gas (atomic and molecular) in galaxies as a function of redshift, and we use these predictions to construct number counts for the next generation of all-sky neutral atomic hydrogen (HI) surveys. Despite the different implementations of the physical ingredients of galaxy formation, we find that the model predictions are broadly consistent with one another; the key differences reflect how the models treat AGN feedback and how the timescale for star formation evolves with redshift. The models produce mass functions of cold gas in galaxies that are generally in good agreement with HI surveys at $z$=0. Interestingly we find that these mass functions do not evolve significantly with redshift. Adopting a simple conversion factor for cold gas mass to HI mass that we apply to all galaxies at all redshifts, we derive mass functions of HI in galaxies from the predicted mass functions of cold gas, which we use to predict the number counts of sources likely to be detected by HI surveys on next generation radio telescopes such as the Square Kilometre Array and its pathfinders. We show also how adopting a conversion factor for cold gas mass to HI mass that varies from galaxy to galaxy impacts on number counts. In addition, we examine how the typical angular sizes of galaxies vary with redshift. Taken together, these results make clear that forthcoming HI surveys will provide important and powerful tests of theoretical galaxy formation models., Comment: 19 pages, 10 figures, accepted for publication in MNRAS. Version contains a number of substantial changes, including a detailed investigation of how the conversion from cold gas mass to HI mass impacts on the number counts of HI sources

We investigate the influence of the external tidal field of a dark matter halo on the dynamical evolution of star clusters using direct N-body simulations, where we assume that the halo is described by a Navarro, Frenk & White mass profile which has an inner density cusp. We assess how varying the mass and concentration of the halo affects the rate at which the star cluster loses mass and we find that increasing halo mass and concentration drives enhanced mass loss rates and in principle shorter cluster disruption timescales. In addition, we examine disruption timescales in a three-component model of a galaxy (bulge, disk and dark matter halo) and find good agreement with results based on an empirical model of the Galactic potential if we assume a halo mass of ~1e12 solar masses. In general, dark matter halos are expected to contribute significantly to the masses of galaxies and should not be ignored when modelling the evolution of star clusters. We extend our results to discuss how this can have a potentially profound effect on the disruption timescales of globular clusters, suggesting that we may underestimate the rate at which primordial globular clusters are disrupted., Comment: 8 pages, 4 figures, accepted for publication in New Astronomy

Globular clusters are dense stellar systems that have typical ages of ~13 billion years, implying that they formed at redshifts of z>~6. Massive stars in newly formed or primordial globular clusters could have played an important role during the epoch of cosmological reionisation (z>~6) as sources of energetic, neutral hydrogen ionising UV photons. We investigate whether or not these stars could have been as important in death as sources of energetic X-ray photons as they were during their main sequence lives. Most massive stars are expected to form in binaries, and an appreciable fraction of these (as much as ~30%) will evolve into X-ray luminous (L_X~10^38 erg/s) high-mass X-ray binaries (HMXBs). These sources would have made a contribution to the X-ray background at z>~6. Using Monte Carlo models of a globular cluster, we estimate the total X-ray luminosity of a population of HMXBs. We compare and contrast this with the total UV luminosity of the massive stars during their main sequence lives. For reasonable estimates, we find that the bolometric luminosity of the cluster peaks at ~10^42 erg/s during the first few million years, but declines to ~10^41 erg/s after ~5 million years as the most massive stars evolve off the main sequence. From this time onwards, the total bolometric luminosity is dominated by HMXBs and falls gradually to ~10^40 erg/s after ~50 million years. Assuming a power-law spectral energy distribution for the HMXBs, we calculate the effective number of neutral hydrogen ionisations per HMXB and show that HMXBs can be as important as sources of ionising radiation as massive stars. Finally we discuss the implications of our results for modelling galaxy formation at high redshift and the prospects of using globular clusters as probes of reionisation., Comment: 7 pages, 3 figures, accepted for publication in MNRAS

[ABRIDGED VERSION] The future of cm and m-wave astronomy lies with the Square Kilometre Array (SKA), a telescope under development by a consortium of 17 countries. The SKA will be 50 times more sensitive than any existing radio facility. A majority of the key science for the SKA will be addressed through large-area imaging of the Universe at frequencies from 300 MHz to a few GHz. The Australian SKA Pathfinder (ASKAP) is aimed squarely in this frequency range, and achieves instantaneous wide-area imaging through the development and deployment of phase-array feed systems on parabolic reflectors. This large field-of-view makes ASKAP an unprecedented synoptic telescope poised to achieve substantial advances in SKA key science. The central core of ASKAP will be located at the Murchison Radio Observatory in inland Western Australia, one of the most radio-quiet locations on the Earth and one of the sites selected by the international community as a potential location for the SKA. Following an introductory description of ASKAP, this document contains 7 chapters describing specific science programmes for ASKAP. The combination of location, technological innovation and scientific program will ensure that ASKAP will be a world-leading radio astronomy facility, closely aligned with the scientific and technical direction of the SKA. A brief summary chapter emphasizes the point, and considers discovery space., Comment: Accepted for publication in Experimental Astronomy. This submission contains only a summary of the full 120 page article which can be found at http://wwwatnf.atnf.csiro.au/projects/askap/newdocs/askap_expast08.pdf . Jasper Wall (UBC, Overall Editor)