• ### Aurigaia: mock Gaia DR2 stellar catalogues from the Auriga cosmological simulations(1804.08549)

April 27, 2018 astro-ph.GA
We present and analyse mock stellar catalogues that match the selection criteria and observables (including uncertainties) of the Gaia satellite data release 2 (DR2). The source are six cosmological high-resolution magneto-hydrodynamic $\Lambda$CDM zoom simulations of the formation of Milky Way analogues from the AURIGA project. Mock data are provided for stars with $V < 16$ mag, and $V < 20$ mag at $|b|>20$ degrees. The mock catalogues are made using two different methods: the public SNAPDRAGONS code, and a method based on that of Lowing et al. that preserves the phase-space distribution of the model stars. The catalogues contain 5-parameter astrometry, radial velocities, multi-band photometry, stellar parameters, dust extinction values, and uncertainties in all these quantities. In addition, we provide theoretical quantities from the simulations, including the gravitational potential and information on the origin of each star. By way of demonstration we apply the mock catalogues to the young stellar disc and the stellar halo. We show that: i)~the young ($< 300$ Myr) outer stellar disc exhibits a flared distribution that is reflected in the height and vertical velocity distribution of $\rm A$- and $\rm B$-dwarf stars up to radii comparable to $\sim 15$ kpc, and ii)~the spin of the stellar halo out to 100~kpc can be accurately measured with Gaia DR2 RR Lyrae stars. We discuss the limitations of both methods and conclude that these mock stellar catalogues, which we make publicly available, are well suited for comparisons with observations and should help to: i)~develop and test analysis methods for the Gaia DR2 data, ii)~gauge the limitations and biases of the data and iii)~ interpret the data in the light of theoretical predictions from realistic ab initio simulations of galaxy formation in the $\Lambda$CDM cosmological model.
• ### Comparing galaxy formation in semi-analytic models and hydrodynamical simulations(1709.08647)

Sept. 27, 2017 astro-ph.GA
It is now possible for hydrodynamical simulations to reproduce a representative galaxy population. Accordingly, it is timely to assess critically some of the assumptions of traditional semi-analytic galaxy formation models. We use the Eagle simulations to assess assumptions built into the Galform semi-analytic model, focussing on those relating to baryon cycling, angular momentum and feedback. We show that the assumption in Galform that newly formed stars have the same specific angular momentum as the total disc leads to a significant overestimate of the total stellar specific angular momentum of disc galaxies. In Eagle, stars form preferentially out of low specific angular momentum gas in the interstellar medium (ISM) due to the assumed gas density threshold for stars to form, leading to more realistic galaxy sizes. We find that stellar mass assembly is similar between Galform and Eagle but that the evolution of gas properties is different, with various indications that the rate of baryon cycling in Eagle is slower than is assumed in Galform. Finally, by matching individual galaxies between Eagle and Galform, we find that an artificial dependence of AGN feedback and gas infall rates on halo mass doubling events in Galform drives most of the scatter in stellar mass between individual objects. Put together our results suggest that the Galform semi-analytic model can be significantly improved in light of recent advances.
• ### Addressing the too big to fail problem with baryon physics and sterile neutrino dark matter(1611.00005)

Aug. 22, 2017 astro-ph.CO, astro-ph.GA
N-body dark matter simulations of structure formation in the $\Lambda$CDM model predict a population of subhalos within Galactic halos that have higher central densities than inferred for satellites of the Milky Way, a tension known as the too big to fail' problem. Proposed solutions include baryonic effects, a smaller mass for the Milky Way halo, and warm dark matter. We test these three possibilities using a semi-analytic model of galaxy formation to generate luminosity functions for Milky Way halo-analogue satellite populations, the results of which are then coupled to the Jiang & van den Bosch model of subhalo stripping to predict the subhalo $V_\mathrm{max}$ functions for the 10 brightest satellites. We find that selecting the brightest satellites (as opposed to the most massive) and modelling the expulsion of gas by supernovae at early times increases the likelihood of generating the observed Milky Way satellite $V_\mathrm{max}$ function. The preferred halo mass is $6\times10^{11}M_{\odot}$, which has a 14 percent probability to host a $V_\mathrm{max}$ function like that of the Milky Way satellites. This probability is reduced to 8 percent for a $1.0\times10^{12}M_{\odot}$ halo and to 3 percent for a $1.4\times10^{12}M_{\odot}$ halo. We conclude that the Milky Way satellite $V_\mathrm{max}$ function is compatible with a CDM cosmology, as previously found by Sawala et al. using hydrodynamic simulations. Sterile neutrino-warm dark matter models achieve a higher degree of agreement with the observations, with a maximum 35 percent chance of generating the observed Milky Way satellite $V_\mathrm{max}$ function. However, more work is required to check that the semi-analytic stripping model is calibrated correctly in the sterile neutrino cosmology, and to check if our sterile neutrino models produce sufficient numbers of faint satellites.
• ### The Pan-STARRS1 Medium-deep Survey: Star Formation Quenching in Group and Cluster Environments(1708.02850)

Aug. 9, 2017 astro-ph.GA
We make use of a catalog of 1600 Pan-STARRS1 groups produced by the probability friends-of-friends algorithm to explore how the galaxy properties, i.e. the specific star formation rate (SSFR) and quiescent fraction, depend on stellar mass and group-centric radius. The work is the extension of Lin et al. (2014). In this work, powered by a stacking technique plus a background subtraction for contamination removal, a finer correction and more precise results are obtained than in our previous work. We find that while the quiescent fraction increases with decreasing group-centric radius the median SSFRs of star-forming galaxies in groups at fixed stellar mass drop slightly from the field toward the group center. This suggests that the major quenching process in groups is likely a fast mechanism. On the other hand, a reduction in SSFRs by ~0.2 dex is seen inside clusters as opposed to the field galaxies. If the reduction is attributed to the slow quenching effect, the slow quenching process acts dominantly in clusters. In addition, we also examine the density-color relation, where the density is defined by using a sixth-nearest neighbor approach. Comparing the quiescent fractions contributed from the density and radial effect, we find that the density effect dominates over the massive group or cluster galaxies, and the radial effect becomes more effective in less massive galaxies. The results support mergers and/or starvation as the main quenching mechanisms in the group environment, while harassment and/or starvation dominate in clusters.
• ### The core-cusp problem: a matter of perspective(1707.06303)

July 19, 2017 astro-ph.GA
The existence of two kinematically and chemically distinct stellar subpopulations in the Sculptor and Fornax dwarf galaxies offers the opportunity to constrain the density profile of their matter haloes by measuring the mass contained within the well-separated half-light radii of the two metallicity subpopulations. Walker and Pe\~{n}arrubia have used this approach to argue that data for these galaxies are consistent with constant-density cores' in their inner regions and rule out `cuspy' Navarro-Frenk-White (NFW) profiles with high statistical significance, particularly in the case of Sculptor. We test the validity of these claims using dwarf galaxies in the APOSTLE $\Lambda$CDM cosmological hydrodynamics simulations of analogues of the Local Group. These galaxies all have NFW dark matter density profiles and a subset of them develop two distinct metallicity subpopulations reminiscent of Sculptor and Fornax. We apply a method analogous to that of Walker and Penarrubia to a sample of 53 simulated dwarfs and find that this procedure often leads to a statistically significant detection of a core in the profile when in reality there is a cusp. Although multiple factors contribute to these failures, the main cause is a violation of the assumption of spherical symmetry upon which the mass estimators are based. The stellar populations of the simulated dwarfs tend to be significantly elongated and, in several cases, the two metallicity populations have different asphericity and are misaligned. As a result, a wide range of slopes of the density profile are inferred depending on the angle from which the galaxy is viewed.
• ### Globular clusters vs dark matter haloes in strong lensing observations(1707.01849)

July 6, 2017 astro-ph.CO, astro-ph.GA
Small distortions in the images of Einstein rings or giant arcs offer the exciting prospect of detecting dark matter haloes or subhaloes of mass below $10^9$M$_{\odot}$, most of which are too small to have made a visible galaxy. A very large number of such haloes are predicted to exist in the cold dark matter model of cosmogony; in contrast other models, such as warm dark matter, predict no haloes below a mass of this order which depends on the properties of the warm dark matter particle. Attempting to detect these small perturbers could therefore discriminate between different kinds of dark matter particles, and even rule out the cold dark matter model altogether. Globular clusters in the lens galaxy also induce distortions in the image which could, in principle, contaminate the test. Here, we investigate the population of globular clusters in six early type galaxies in the Virgo cluster. We find that the number density of globular clusters of mass $\sim10^6$M$_{\odot}$ is comparable to that of the dark matter perturbers (including subhaloes in the lens and haloes along the line-of-sight). We show that the very different degrees of mass concentration in globular clusters and dark matter haloes result in different lensing distortions. These are detectable with milli-arcsecond resolution imaging which can distinguish between globular cluster and dark matter halo signals.
• ### Revisiting HOD model assumptions: the impact of AGN feedback and assembly bias(1705.01988)

May 1, 2017 astro-ph.GA
The standard Halo Occupation Distribution (HOD) models were originally developed based on results from semi-analytic and hydrodynamical galaxy formation models. Those models have since progressed, in particular to include AGN feedback to match the galaxy luminosity function in a universe with the observed baryon fraction. AGN feedback affects the relationship between galaxy stellar mass and luminosity, in particular making the relationship non-monotonic. For matched number density samples, galaxies in luminosity-threshold samples occupy a different range of halo masses from those in stellar-mass-threshold samples. We find that the shapes of the HODs of luminosity-threshold samples are slightly more complicated in semi-analytic galaxy formation models that include AGN feedback than are assumed by standard HOD models. We also find that subhalo abundance matching (SHAM) does not preserve these non-standard shapes. We show that catalogues created using SHAM and the semi-analytic model Galform that have the same large-scale 2-point clustering by construction have different void probability functions (VPFs) in both real and redshift space. We find that these differences arise from the different HOD shapes, as opposed to assembly bias, which indicates that the VPF could be used to test the suitability of an HOD model with real data.
• ### Projection effects in the strong lensing study of subhaloes(1612.06227)

Dec. 19, 2016 astro-ph.CO
The defining characteristic of the cold dark matter (CDM) hypothesis is the presence of a very large number of low-mass haloes, too small to have made a visible galaxy. Other hypotheses for the nature of the dark matter, such as warm dark matter (WDM), predict a much smaller number of such low-mass haloes. Strong lensing systems offer the possibility of detecting small-mass haloes through the distortions they induce in the lensed image. Here we show that the main contribution to the image distortions comes from haloes along the line of sight rather than subhaloes in the lens as has normally been assumed so far. These interlopers enhance the differences between the predictions of CDM and WDM models. We derive the total perturber mass function, including both subhaloes and interlopers, and show that measurements of approximately 20 strong lens systems with a detection limit of $M_{\rm low}=10^7 h^{-1} M_{\odot}$ would distinguish (at 3 sigma) between CDM and a WDM model consisting of 7 keV sterile neutrinos such as those required to explain the recently detected 3.5 keV X-ray emission line from the centres of galaxies and clusters.
• DESI (Dark Energy Spectroscopic Instrument) is a Stage IV ground-based dark energy experiment that will study baryon acoustic oscillations (BAO) and the growth of structure through redshift-space distortions with a wide-area galaxy and quasar redshift survey. To trace the underlying dark matter distribution, spectroscopic targets will be selected in four classes from imaging data. We will measure luminous red galaxies up to $z=1.0$. To probe the Universe out to even higher redshift, DESI will target bright [O II] emission line galaxies up to $z=1.7$. Quasars will be targeted both as direct tracers of the underlying dark matter distribution and, at higher redshifts ($2.1 < z < 3.5$), for the Ly-$\alpha$ forest absorption features in their spectra, which will be used to trace the distribution of neutral hydrogen. When moonlight prevents efficient observations of the faint targets of the baseline survey, DESI will conduct a magnitude-limited Bright Galaxy Survey comprising approximately 10 million galaxies with a median $z\approx 0.2$. In total, more than 30 million galaxy and quasar redshifts will be obtained to measure the BAO feature and determine the matter power spectrum, including redshift space distortions.
• DESI (Dark Energy Spectropic Instrument) is a Stage IV ground-based dark energy experiment that will study baryon acoustic oscillations and the growth of structure through redshift-space distortions with a wide-area galaxy and quasar redshift survey. The DESI instrument is a robotically-actuated, fiber-fed spectrograph capable of taking up to 5,000 simultaneous spectra over a wavelength range from 360 nm to 980 nm. The fibers feed ten three-arm spectrographs with resolution $R= \lambda/\Delta\lambda$ between 2000 and 5500, depending on wavelength. The DESI instrument will be used to conduct a five-year survey designed to cover 14,000 deg$^2$. This powerful instrument will be installed at prime focus on the 4-m Mayall telescope in Kitt Peak, Arizona, along with a new optical corrector, which will provide a three-degree diameter field of view. The DESI collaboration will also deliver a spectroscopic pipeline and data management system to reduce and archive all data for eventual public use.
• ### Satellite galaxies in semi-analytic models of galaxy formation with sterile neutrino dark matter(1511.04078)

Nov. 7, 2016 astro-ph.CO, astro-ph.GA
The sterile neutrino is a viable dark matter candidate that can be produced in the early Universe via non-equilibrium processes, and would therefore possess a highly non-thermal spectrum of primordial velocities. In this paper we analyse the process of structure formation with this class of dark matter particles. To this end we construct primordial dark matter power spectra as a function of the lepton asymmetry, $L_6$, that is present in the primordial plasma and leads to resonant sterile neutrino production. We compare these power spectra with those of thermally produced dark matter particles and show that resonantly produced sterile neutrinos are much colder than their thermal relic counterparts. We also demonstrate that the shape of these power spectra is not determined by the free-streaming scale alone. We then use the power spectra as an input for semi-analytic models of galaxy formation in order to predict the number of luminous satellite galaxies in a Milky Way-like halo. By assuming that the mass of the Milky Way halo must be no more than $2\times10^{12}M_{\odot}$ (the adopted upper bound based on current astronomical observations) we are able to constrain the value of $L_6$ for $M_s\le 8$~keV. We also show that the range of $L_6$ that is in best agreement with the 3.5~keV line (if produced by decays of 7~keV sterile neutrino) requires that the Milky Way halo has a mass no smaller than $1.5\times10^{12}M_{\odot}$. Finally, we compare the power spectra obtained by direct integration of the Boltzmann equations for a non-resonantly produced sterile neutrino with the fitting formula of Viel~et~al. and find that the latter significantly underestimates the power amplitude on scales relevant to satellite galaxies.
• ### Modelling Galaxy Merger Timescales and Tidal Destruction(1609.09520)

Sept. 29, 2016 astro-ph.CO, astro-ph.GA
We present a model for the dynamical evolution of subhaloes based on an approach combining numerical and analytical methods. Our method is based on tracking subhaloes in an N-body simulation up to the last point that it can be resolved, and applying an analytic prescription for its merger timescale that takes dynamical friction and tidal disruption into account. When applied to cosmological N-body simulations with mass resolutions that differ by two orders of magnitude, the technique produces halo occupation distributions that agree to within 3%.
• ### Galaxies in the EAGLE hydrodynamical simulation and in the Durham and Munich semi-analytical models(1512.00015)

May 27, 2016 astro-ph.CO, astro-ph.GA
We compare global predictions from the EAGLE hydrodynamical simulation, and two semi-analytic (SA) models of galaxy formation, L-GALAXIES and GALFORM. All three models include the key physical processes for the formation and evolution of galaxies and their parameters are calibrated against a small number of observables at $z\approx 0$. The two SA models have been applied to merger trees constructed from the EAGLE dark matter only simulation. We find that at $z\leq 2$, both the galaxy stellar mass functions for stellar masses $M_*<10^{10.5} M_{\odot}$ and the median specific star formation rates (sSFRs) in the three models agree to better than $0.4$~dex. The evolution of the sSFR predicted by the three models closely follows the mass assembly history of dark matter haloes. In both EAGLE and L-GALAXIES there are more central passive galaxies with $M_*<10^{9.5} M_{\odot}$ than in L-GALAXIES. This difference is related to galaxies that have entered and then left a larger halo and which are treated as satellites in GALFORM. In the range $0<z<1$, the slope of the evolution of the star formation rate density in EAGLE is a factor of $\approx 1.5$ steeper than for the two SA models. The median sizes for galaxies with $M_*>10^{9.5} M_{\odot}$ differ in some instances by an order of magnitude, while the stellar mass-size relation in EAGLE is a factor of $\approx 2$ tighter than for the two SA models. Our results suggest the need for a revision of how SA models treat the effect of baryonic self-gravity on the underlying dark matter. The treatment of gas flows in the models needs to be revised based on detailed comparison with observations to understand in particular the evolution of the stellar mass-metallicity relation.
• ### The Clustering Evolution of Dusty Star-Forming Galaxies(1504.04516)

May 5, 2016 astro-ph.CO
We present predictions for the clustering of galaxies selected by their emission at far infra-red (FIR) and sub-millimetre wavelengths. This includes the first predictions for the effect of clustering biases induced by the coarse angular resolution of single-dish telescopes at these wavelengths. We combine a new version of the GALFORM model of galaxy formation with a self-consistent model for calculating the absorption and re-emission of radiation by interstellar dust. Model galaxies selected at $850$ $\mu$m reside in dark matter halos of mass $M_{\rm halo}\sim10^{11.5}-10^{12}$ $h^{-1}$ M$_{\odot}$, independent of redshift (for $0.2\lesssim z\lesssim4$) or flux (for $0.25\lesssim S_{850\mu\rm m}\lesssim4$ mJy). At $z\sim2.5$, the brightest galaxies ($S_{850\mu\rm m}>4$ mJy) exhibit a correlation length of $r_{0}=5.5_{-0.5}^{+0.3}$ $h^{-1}$ Mpc, consistent with observations. We show that these galaxies have descendants with stellar masses $M_{\star}\sim10^{11}$ $h^{-1}$ M$_{\odot}$ occupying halos spanning a broad range in mass $M_{\rm halo}\sim10^{12}-10^{14}$ $h^{-1}$ M$_{\odot}$. The FIR emissivity at shorter wavelengths ($250$, $350$ and $500$ $\mu$m) is also dominated by galaxies in the halo mass range $M_{\rm halo}\sim10^{11.5}-10^{12}$ $h^{-1}$ M$_{\odot}$, again independent of redshift (for $0.5\lesssim z\lesssim5$). We compare our predictions for the angular power spectrum of cosmic infra-red background anisotropies at these wavelengths with observations, finding agreement to within a factor of $\sim2$ over all scales and wavelengths, an improvement over earlier versions of the model. Simulating images at $850$ $\mu$m, we show that confusion effects boost the measured angular correlation function on all scales by a factor of $\sim4$. This has important consequences, potentially leading to inferred halo masses being overestimated by an order of magnitude.
• ### Constraints on the identity of the dark matter from strong gravitational lenses(1512.06507)

April 5, 2016 astro-ph.CO
The cold dark matter (CDM) cosmological model unambigously predicts that a large number of haloes should survive as subhaloes when they are accreted into a larger halo. The CDM model would be ruled out if such substructures were shown not to exist. By contrast, if the dark matter consists of warm particles (WDM), then below a threshold mass that depends on the particle mass far fewer substructures would be present. Finding subhaloes below a certain mass would then rule out warm particle masses below some value. Strong gravitational lensing provides a clean method to measure the subhalo mass function through distortions in the structure of Einstein rings and giant arcs.Using mock lensing observations constructed from high-resolution N-body simulations, we show that measurements of approximately 100 strong lens systems with a detection limit of $10^7 h^{-1} M_{\odot}$ would clearly distinguish CDM from WDM in the case where this consists of 7 keV sterile neutrinos such as those that might be responsible for the 3.5 keV X-ray emission line recently detected in galaxies and clusters.
• ### The orbital PDF: general inference of the gravitational potential from steady-state tracers(1507.00769)

Dec. 4, 2015 astro-ph.CO, astro-ph.GA
We develop two general methods to infer the gravitational potential of a system using steady-state tracers, i.e., tracers with a time-independent phase-space distribution. Combined with the phase-space continuity equation, the time independence implies a universal Orbital Probability Density Function (oPDF) $\mathrm{d} P(\lambda|{\rm orbit})\propto \mathrm{d} t$, where $\lambda$ is the coordinate of the particle along the orbit. The oPDF is equivalent to Jeans theorem, and is the key physical ingredient behind most dynamical modelling of steady-state tracers. In the case of a spherical potential, we develop a likelihood estimator that fits analytical potentials to the system, and a non-parametric method ("phase-mark") that reconstructs the potential profile, both assuming only the oPDF. The methods involve no extra assumptions about the tracer distribution function and can be applied to tracers with any arbitrary distribution of orbits, with possible extension to non-spherical potentials. The methods are tested on Monte Carlo samples of steady-state tracers in dark matter haloes to show that they are unbiased as well as efficient. A fully documented \textsc{C/Python} code implementing our method is freely available at a GitHub repository linked from \url{http://icc.dur.ac.uk/data/#oPDF}.
• ### The orbital PDF: the dynamical state of Milky Way sized haloes and the intrinsic uncertainty in the determination of their masses(1507.00771)

Nov. 29, 2015 astro-ph.CO, astro-ph.GA
Using realistic cosmological simulations of Milky Way sized haloes, we study their dynamical state and the accuracy of inferring their mass profiles with steady-state models of dynamical tracers. We use a new method that describes the phase-space distribution of a steady-state tracer population in a spherical potential without any assumption regarding the distribution of their orbits. Applying the method to five haloes from the Aquarius $\Lambda$CDM N-body simulation, we find that dark matter particles are an accurate tracer that enables the halo mass and concentration parameters to be recovered with an accuracy of $5\%$. Assuming a potential profile of the NFW form does not significantly affect the fits in most cases, except for halo A whose density profile differs significantly from the NFW form, leading to a $30\%$ bias in the dynamically fitted parameters. The existence of substructures in the dark matter tracers only affects the fits by $\sim 1\%$. Applying the method to mock stellar haloes generated by a particle-tagging technique, we find the stars are farther from equilibrium than dark matter particles, yielding a systematic bias of $\sim 20\%$ in the inferred mass and concentration parameter. The level of systematic biases obtained from a conventional distribution function fit to stars is comparable to ours, while similar fits to DM tracers are significantly biased in contrast to our fits. In line with previous studies, the mass bias is much reduced near the tracer half-mass radius.
• ### Recovering dark-matter clustering from galaxies with Gaussianization(1511.02034)

Nov. 6, 2015 astro-ph.CO
The Gaussianization transform has been proposed as a method to remove the issues of scale-dependent galaxy bias and nonlinearity from galaxy clustering statistics, but these benefits have yet to be thoroughly tested for realistic galaxy samples. In this paper, we test the effectiveness of the Gaussianization transform for different galaxy types by applying it to realistic simulated blue and red galaxy samples. We show that in real space, the shapes of the Gaussianized power spectra of both red and blue galaxies agree with that of the underlying dark matter, with the initial power spectrum, and with each other to smaller scales than do the statistics of the usual (untransformed) density field. However, we find that the agreement in the Gaussianized statistics breaks down in redshift space. We attribute this to the fact that red and blue galaxies exhibit very different fingers of god in redshift space. After applying a finger-of-god compression, the agreement on small scales between the Gaussianized power spectra is restored. We also compare the Gaussianization transform to the clipped galaxy density field and find that while both methods are effective in real space, they have more complicated behaviour in redshift space. Overall, we find that Gaussianization can be useful in recovering the shape of the underlying dark matter power spectrum to k ~ 0.5 h/Mpc and of the initial power spectrum to k ~ 0.4 h/Mpc in certain cases at z = 0.
• ### The evolution of the stellar mass versus halo mass relationship(1510.08463)

Nov. 4, 2015 astro-ph.GA
We present an analysis of the predictions made by the Galform semi-analytic galaxy formation model for the evolution of the relationship between stellar mass and halo mass. We show that for the standard implementations of supernova feedback and gas reincorporation used in semi-analytic models, this relationship is predicted to evolve weakly over the redshift range 0<z<4. Modest evolution in the median stellar mass versus halo mass (SHM) relationship implicitly requires that, at fixed halo mass, the efficiency of stellar mass assembly must be almost constant with cosmic time. We show that in our model, this behaviour can be understood in simple terms as a result of a constant efficiency of gas reincorporation, and an efficiency of SNe feedback that is, on average, constant at fixed halo mass. We present a simple explanation of how feedback from active galactic nuclei (AGN) acts in our model to introduce a break in the SHM relation whose location is predicted to evolve only modestly. Finally, we show that if modifications are introduced into the model such that, for example, the gas reincorporation efficiency is no longer constant, the median SHM relation is predicted to evolve significantly over 0<z<4. Specifically, we consider modifications that allow the model to better reproduce either the evolution of the stellar mass function or the evolution of average star formation rates inferred from observations.
• ### Formation of In Situ Stellar Haloes in Milky Way-Mass Galaxies(1501.04630)

Oct. 4, 2015 astro-ph.GA
We study the formation of stellar haloes in three Milky Way-mass galaxies using cosmological smoothed particle hydrodynamics simulations, focusing on the subset of halo stars that form in situ, as opposed to those accreted from satellites. In situ stars in our simulations dominate the stellar halo out to 20 kpc and account for 30-40 per cent of its total mass. We separate in situ halo stars into three straightforward, physically distinct categories according to their origin: stars scattered from the disc of the main galaxy ('heated disc'), stars formed from gas smoothly accreted on to the halo ('smooth' gas) and stars formed in streams of gas stripped from infalling satellites ('stripped' gas). We find that most belong to the stripped gas category. Those originating in smooth gas outside the disc tend to form at the same time and place as the stripped-gas population, suggesting that their formation is associated with the same gas-rich accretion events. The scattered disc star contribution is negligible overall but significant in the Solar neighbourhood, where ~90 per cent of stars on eccentric orbits once belonged to the disc. However, the distinction between halo and thick disc in this region is highly ambiguous. The chemical and kinematic properties of the different components are very similar at the present day, but the global properties of the in situ halo differ substantially between the three galaxies in our study. In our simulations, the hierarchical buildup of structure is the driving force behind not only the accreted stellar halo, but also those halo stars formed in situ.
• ### Detection of a Supervoid Aligned with the Cold Spot of the Cosmic Microwave Background(1405.1566)

March 5, 2015 astro-ph.CO
We use the WISE-2MASS infrared galaxy catalog matched with Pan-STARRS1 (PS1) galaxies to search for a supervoid in the direction of the Cosmic Microwave Background Cold Spot. Our imaging catalog has median redshift $z\simeq 0.14$, and we obtain photometric redshifts from PS1 optical colours to create a tomographic map of the galaxy distribution. The radial profile centred on the Cold Spot shows a large low density region, extending over 10's of degrees. Motivated by previous Cosmic Microwave Background results, we test for underdensities within two angular radii, $5^\circ$, and $15^\circ$. The counts in photometric redshift bins show significantly low densities at high detection significance, $\gtrsim 5 \sigma$ and $\gtrsim 6 \sigma$, respectively, for the two fiducial radii. The line-of-sight position of the deepest region of the void is $z\simeq 0.15-0.25$. Our data, combined with an earlier measurement by Granett et al. 2010, are consistent with a large $R_{\rm void}=(220 \pm 50) h^{-1}Mpc$ supervoid with $\delta_{m} \simeq -0.14 \pm 0.04$ centered at $z=0.22\pm0.03$. Such a supervoid, constituting at least a $\simeq 3.3\sigma$ fluctuation in a Gaussian distribution of the $\Lambda CDM$ model, is a plausible cause for the Cold Spot.
• ### Orbital parameters of infalling satellite haloes in the hierarchical $\Lambda$CDM model(1409.1179)

Jan. 27, 2015 astro-ph.CO
We present distributions of orbital parameters of infalling satellites of $\Lambda$CDM haloes in the mass range $10^{12}-10^{14}$M$_\odot$, which represent the initial conditions for the subsequent evolution of substructures within the host halo. We use merger trees constructed in a high resolution cosmological N-body simulation to trace satellite haloes, and identify the time of infall. We find signficant trends in the distribution of orbital parameters with both the host halo mass and the ratio of satellite-to-host halo masses. For all host halo masses, satellites whose infall mass is a larger fraction of the host halo mass have more eccentric, radially biased orbits. At fixed satellite-to-host halo mass ratio, high mass haloes are biased towards accreting satellites on slightly more radial orbits. To charactise the orbital distributions fully requires fitting the correlated bivariate distribution of two chosen orbital parameters (e.g. radial and tangential velocity or energy and angular momentum). We provide simple fits to one choice of the bivariate distributions, which when transformed faithfully, captures the behaviour of any of the projected one-dimensional distributions.
• ### Simulated observations of sub-millimetre galaxies: the impact of single-dish resolution and field variance(1406.0855)

Oct. 24, 2014 astro-ph.CO, astro-ph.GA
Recent observational evidence suggests that the coarse angular resolution ($\sim20"$ FWHM) of single-dish telescopes at sub-mm wavelengths has biased the observed galaxy number counts by blending together the sub-mm emission from multiple sub-mm galaxies (SMGs). We use lightcones computed from an updated implementation of the \galform semi-analytic model to generate $50$ mock sub-mm surveys of $0.5$ deg$^2$ at $850$ $\mu$m, taking into account the effects of the finite single-dish beam in a more accurate way than has been done previously. We find that blending of SMGs does lead to an enhancement of source extracted number counts at bright fluxes ($S_{\mathrm{850}\mu\mathrm{m}}\gtrsim1$ mJy). Typically, $\sim3{-}6$ galaxies contribute $90\%$ of the flux of an $S_{850\mu\mathrm{m}}=5$ mJy source and these blended galaxies are physically unassociated. We find that field-to-field variations are comparable to Poisson fluctuations for our $S_{850\mu\mathrm{m}}>5$ mJy SMG population, which has a median redshift $z_{50}=2.0$, but are greater than Poisson for the $S_{850\mu\mathrm{m}}>1$ mJy population ($z_{50}=2.8$). In a detailed comparison to a recent interferometric survey targeted at single-dish detected sources, we reproduce the difference between single-dish and interferometer number counts and find a median redshift ($z_{50}=2.5$) in excellent agreement with the observed value ($z_{50}=2.5\pm 0.2$). We also present predictions for single-dish survey number counts at $450$ and $1100$ $\mu$m, which show good agreement with observational data.
• ### Supervoid Origin of the Cold Spot in the Cosmic Microwave Background(1407.1470)

July 6, 2014 astro-ph.CO
We use a WISE-2MASS-Pan-STARRS1 galaxy catalog to search for a supervoid in the direction of the Cosmic Microwave Background Cold Spot. We obtain photometric redshifts using our multicolor data set to create a tomographic map of the galaxy distribution. The radial density profile centred on the Cold Spot shows a large low density region, extending over 10's of degrees. Motivated by previous Cosmic Microwave Background results, we test for underdensities within two angular radii, $5^\circ$, and $15^\circ$. Our data, combined with an earlier measurement by Granett et al 2010, are consistent with a large $R_{\rm void}=(192 \pm 15)h^{-1} Mpc$ $(2\sigma)$ supervoid with $\delta \simeq -0.13 \pm 0.03$ centered at $z=0.22\pm0.01$. Such a supervoid, constituting a $\sim3.5 \sigma$ fluctuation in the $\Lambda CDM$ model, is a plausible cause for the Cold Spot.
• ### Herschel-ATLAS/GAMA:How does the far-IR luminosity function depend on galaxy group properties?(1401.0986)

July 4, 2014 astro-ph.CO
We use the Herschel ATLAS (H-ATLAS) Phase I data to study the conditional luminosity function of far-IR (250 um) selected galaxies in optically-selected galaxy groups from the GAMA spectroscopic survey, as well as environmental effects on the far-IR-to-optical colour. We applied two methods, which gave consistent results for the far-IR conditional luminosity functions. The direct matching method matches H-ATLAS sources to GAMA/SDSS galaxies, then links the optical counterparts to GAMA groups. The stacking method counts the number of far-IR sources within the projected radii of GAMA groups, subtracting the local background. We investigated the dependence of the far-IR (250 um) luminosity function on group mass in the range 10d12 < M_h < 10d14Msun/h and on redshift in the range 0 < z < 0.4, using a sample of 3000 groups containing H-ATLAS sources with GAMA redshifts over an area of 126 sq. deg. We find that the characteristic 250 um luminosity, L^*(250), increases with group mass up to M_h ~ 10d13Msun/h, but is roughly constant above this, while it increases with redshift at high group masses, but less so at low masses. We also find that the group far-IR luminosity-to-mass ratio L(250)/M_h increases with redshift and is higher in low-mass groups. We estimate that around 70% of the 250 um luminosity density in the local universe is contributed by groups with M_h > 10d12Msun/h. We also find that the far-IR-to-optical colours of H-ATLAS galaxies are independent of group mass over the range 10d12 < M_h <10d14Msun/h in the local universe. We also compare our observational results with recent semi-analytical models, and find that none of these galaxy formation model can reproduce the conditional far-IR luminosity functions of galaxy groups.