• ### A comparison between semi-analytical gas cooling models and hydrodynamical simulations(1803.01923)

March 5, 2018 astro-ph.GA
We compare the mass cooling rates and cumulative cooled-down masses predicted by several semi-analytical (SA) cooling models with hydrodynamical simulations performed using the grid-based code AREPO. The SA cooling models are the new cooling model introduced in Hou et al. (2017), along with two earlier GALFORM cooling models and the cooling models in the L-GALAXIES and MORGANA SA galaxy formation models. We find that the predictions of the new cooling model generally agree the best with the simulations. For halos with $M_{\rm halo}\lesssim 3\times 10^{11}\,{\rm M}_{\odot}$ , the SA models predict that radiative cooling is faster than gravitational infall. Even though SA models assume that gas falls onto galaxies from a spherical gas halo, while the simulations show that the cold gas is accreted through non-spherical filaments, both methods predict similar mass cooling rates, because in both cases the gas accretion onto galaxies is limited by the gravitational infall timescale. For halos with $M_{\rm halo}\gtrsim 10^{12}\,{\rm M}_{\odot}$ , gas in the simulations typically cools from a roughly spherical hot gas halo, as assumed in the SA models, but the halo gas only gradually contracts during cooling, leading to compressional heating. SA models ignore this compressional heating, and so overestimate mass cooling rates by factors of a few. At low redshifts halo major mergers or a sequence of successive smaller mergers are seen in the simulations to strongly heat the hot gas halo and suppress cooling, while mergers at high redshifts do not suppress cooling, because the gas filaments are difficult to heat up. The new SA cooling model best captures these effects of halo mergers.
• ### A new gas cooling model for semi-analytical galaxy formation models(1708.02950)

Dec. 11, 2017 astro-ph.GA
Semi-analytic galaxy formation models are widely used to gain insight into the astrophysics of galaxy formation and in model testing, parameter space searching and mock catalogue building. In this work we present a new model for gas cooling in halos in semi-analytic models, which improves over previous cooling models in several ways. Our new treatment explicitly includes the evolution of the density profile of the hot gas driven by the growth of the dark matter halo and by the dynamical adjustment of the gaseous corona as gas cools down. The effect of the past cooling history on the current mass cooling rate is calculated more accurately, by doing an integral over the past history. The evolution of the hot gas angular momentum profile is explicitly followed, leading to a self-consistent and more detailed calculation of the angular momentum of the cooled down gas. This model predicts higher cooled down masses than the cooling models previously used in GALFORM, closer to the predictions of the cooling models in L-GALAXIES and MORGANA, even though those models are formulated differently. It also predicts cooled down angular momenta that are higher than in previous GALFORM cooling models, but generally lower than the predictions of L-GALAXIES and MORGANA. When used in a full galaxy formation model, this cooling model improves the predictions for early-type galaxy sizes in GALFORM.
• ### 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.
• ### Understanding the non-linear clustering of high redshift galaxies(1702.00853)

June 6, 2017 astro-ph.CO
We incorporate the non-linear clustering of dark matter halos, as modelled by Jose et al. (2016) into the halo model to better understand the clustering of Lyman break galaxies (LBGs) in the redshift range $z=3-5$. We find that, with this change, the predicted LBG clustering increases significantly on quasi-linear scales ($0.1 \leq r\,/\,h^{-1} \,{\rm Mpc} \leq 10$) compared to that in the linear halo bias model. This in turn results in an increase in the clustering of LBGs by an order of magnitude on angular scales $5" \leq \theta \leq 100"$. Remarkably, the predictions of our new model on the whole remove the systematic discrepancy between the linear halo bias predictions and the observations. The correlation length and large scale galaxy bias of LBGs are found to be significantly higher in the non-linear halo bias model than in the linear halo bias model. The resulting two-point correlation function retains an approximate power-law form in contrast with that computed using the linear halo bias theory. We also find that the non-linear clustering of LBGs increases with increasing luminosity and redshift. Our work emphasizes the importance of using non-linear halo bias in order to model the clustering of high-z galaxies to probe the physics of galaxy formation and extract cosmological parameters reliably.
• 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.
• ### Constraining SN feedback: a tug of war between reionization and the Milky Way satellites(1512.04595)

Aug. 11, 2016 astro-ph.GA
Theoretical models of galaxy formation based on the cold dark matter cosmogony typically require strong feedback from supernova (SN) explosions in order to reproduce the Milky Way satellite galaxy luminosity function and the faint end of the field galaxy luminosity function. However, too strong a SN feedback also leads to the universe reionizing too late, and the metallicities of Milky Way satellites being too low. The combination of these four observations therefore places tight constraints on SN feedback. We investigate these constraints using the semi-analytical galaxy formation model galform. We find that these observations favour a SN feedback model in which the feedback strength evolves with redshift. We find that, for our best fit model, half of the ionizing photons are emitted by galaxies with rest-frame far-UV absolute magnitudes $M_{\rm AB}(1500{\rm \AA})<-17.5$, which implies that already observed galaxy populations contribute about half of the photons responsible for reionization. The $z=0$ descendants of these galaxies are mainly galaxies with stellar mass $M_*>10^{10}\,{\rm M}_{\odot}$ and preferentially inhabit halos with mass $M_{\rm halo}>10^{13}\,{\rm M}_{\odot}$.
• ### A unified multi-wavelength model of galaxy formation(1509.08473)

Aug. 6, 2016 astro-ph.GA
We present a new version of the GALFORM semi-analytical model of galaxy formation. This brings together several previous developments of GALFORM into a single unified model, including a different initial mass function (IMF) in quiescent star formation and in starbursts, feedback from active galactic nuclei supressing gas cooling in massive halos, and a new empirical star formation law in galaxy disks based on their molecular gas content. In addition, we have updated the cosmology, introduced a more accurate treatment of dynamical friction acting on satellite galaxies, and updated the stellar population model. The new model is able to simultaneously explain both the observed evolution of the K-band luminosity function and stellar mass function, and the number counts and redshift distribution of sub-mm galaxies selected at 850 mu. This was not previously achieved by a single physical model within the LambdaCDM framework, but requires having an IMF in starbursts that is somewhat top-heavy. The new model is tested against a wide variety of observational data covering wavelengths from the far-UV to sub-mm, and redshifts from z=0 to z=6, and is found to be generally successful. These observations include the optical and near-IR luminosity functions, HI mass function, fraction of early type galaxies, Tully-Fisher, metallicity-luminosity and size-luminosity relations at z=0, as well as far-IR number counts, and far-UV luminosity functions at z ~ 3-6. [abridged]
• ### The clustering of dark matter halos: scale-dependent bias on quasi-linear scales(1509.06715)

July 18, 2016 astro-ph.CO
We investigate the spatial clustering of dark matter halos, collapsing from $1-4 \sigma$ fluctuations, in the redshift range $0 - 5$ using N-body simulations. The halo bias of high redshift halos ($z \geq 2$) is found to be strongly non-linear and scale-dependent on quasi-linear scales that are larger than their virial radii ($0.5-10$ Mpc/h). However, at lower redshifts, the scale-dependence of non-linear bias is weaker and and is of the order of a few percent on quasi-linear scales at $z \sim 0$. We find that the redshift evolution of the scale dependent bias of dark matter halos can be expressed as a function of four physical parameters: the peak height of halos, the non-linear matter correlation function at the scale of interest, an effective power law index of the {\it rms} linear density fluctuations and the matter density of the universe at the given redshift. This suggests that the scale-dependence of halo bias is not a universal function of the dark matter power spectrum, which is commonly assumed. We provide a fitting function for the scale dependent halo bias as a function of these four parameters. Our fit reproduces the simulation results to an accuracy of better than 4 % over the redshift range $0\leq z \leq 5$. We also extend our model by expressing the non-linear bias as a function of the linear matter correlation function. It is important to incorporate our results into the clustering models of dark matter halos at any redshift, including those hosting early generations of stars and galaxies before reionization.
• ### 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.
• ### Reionisation in sterile neutrino cosmologies(1605.03179)

May 10, 2016 astro-ph.CO, astro-ph.GA
We investigate the process of reionisation in a model in which the dark matter is a warm elementary particle such as a sterile neutrino. We focus on models that are consistent with the dark matter decay interpretation of the recently detected line at 3.5 keV in the X-ray spectra of galaxies and clusters. In warm dark matter models the primordial spectrum of density perturbations has a cut-off on the scale of dwarf galaxies. Structure formation therefore begins later than in the standard cold dark matter (CDM) model and very few objects form below the cut-off mass scale. To calculate the number of ionising photons, we use the Durham semi-analytic model of galaxy formation, GALFORM. We find that even the most extreme 7 keV sterile neutrino we consider is able to reionise the Universe early enough to be compatible with the bounds on the epoch of reionisation from Planck. This, perhaps surprising, result arises from the rapid build-up of high redshift galaxies in the sterile neutrino models which is also reflected in a faster evolution of their far-UV luminosity function between $10>z>7$ than in CDM. The dominant sources of ionising photons are systematically more massive in the sterile neutrino models than in CDM. As a consistency check on the models, we calculate the present-day luminosity function of satellites of Milky Way-like galaxies. When the satellites recently discovered in the DES survey are taken into account, strong constraints are placed on viable sterile neutrino models.
• ### 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.
• ### The environments of high redshift radio galaxies and quasars: probes of protoclusters(1509.01254)

Dec. 10, 2015 astro-ph.CO, astro-ph.GA
We use the GALFORM semi-analytical model to study high density regions traced by radio galaxies and quasars at high redshifts. We explore the impact that baryonic physics has upon the properties of galaxies in these environments. Star-forming emission-line galaxies (Ly{\alpha} and H{\alpha} emitters) are used to probe the environments at high redshifts. Radio galaxies are predicted to be hosted by more massive haloes than quasars, and this is imprinted on the amplitude of galaxy overdensities and cross-correlation functions. We find that Ly{\alpha} radiative transfer and AGN feedback indirectly affect the clustering on small scales and also the stellar masses, star- formation rates and gas metallicities of galaxies in dense environments. We also investigate the relation between protoclusters associated with radio galaxies and quasars, and their present- day cluster descendants. The progenitors of massive clusters associated with radio galaxies and quasars allow us to determine an average protocluster size in a simple way. Overdensities within the protoclusters are found to correlate with the halo descendant masses. We present scaling relations that can be applied to observational data. By computing projection effects due to the wavelength resolution of modern spectrographs and narrow-band filters we show that the former have enough spectral resolution to map the structure of protoclusters, whereas the latter can be used to measure the clustering around radio galaxies and quasars over larger scales to determine the mass of dark matter haloes hosting them.
• ### 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.
• ### Clustering of Extremely Red Objects in Elais-N1 from the UKIDSS DXS with optical photometry from Pan-STARRS1 and Subaru(1311.4624)

Nov. 19, 2013 astro-ph.CO
We measure the angular clustering of 33 415 extremely red objects (EROs) in the Elais-N1 field covering 5.33 deg$^{2}$, which cover the redshift range $z=0.8$ to $2$. This sample was made by merging the UKIDSS Deep eXtragalactic Survey (DXS) with the optical Subaru and Pan-STARRS PS1 datasets. We confirm the existence of a clear break in the angular correlation function at $\sim 0.02^{\circ}$ corresponding to $1 h^{-1}$ Mpc at $z\sim1$. We find that redder or brighter EROs are more clustered than bluer or fainter ones. Halo Occupation Distribution (HOD) model fits imply that the average mass of dark matter haloes which host EROs is over $10^{13} h^{-1} M_{\odot}$ and that EROs have a bias ranging from 2.7 to 3.5. Compared to EROs at $z\sim1.1$, at $z\sim1.5$ EROs have a higher bias and fewer are expected to be satellite galaxies. Furthermore, EROs reside in similar dark matter haloes to those that host $10^{11.0} M_{\odot}<M_{*}<10^{11.5} M_{\odot}$ galaxies. We compare our new measurement and HOD fits with the predictions of the GALFORM semi-analytical galaxy formation model. Overall, the clustering predicted by GALFORM gives an encouraging match to our results. However, compared to our deductions from the measurements, GALFORM puts EROs into lower mass haloes and predicts that a larger fraction of EROs are satellite galaxies. This suggests that the treatment of gas cooling may need to be revised in the model. Our analysis illustrates the potential of clustering analyses to provide observational constraints on theoretical models of galaxy formation.
• ### Extending the halo mass resolution of $N$-body simulations(1310.3880)

Oct. 14, 2013 astro-ph.CO
We present a scheme to extend the halo mass resolution of N-body simulations of the hierarchical clustering of dark matter. The method uses the density field of the simulation to predict the number of sub-resolution dark matter haloes expected in different regions. The technique requires as input the abundance of haloes of a given mass and their average clustering, as expressed through the linear and higher order bias factors. These quantities can be computed analytically or, more accurately, derived from a higher resolution simulation as done here. Our method can recover the abundance and clustering in real- and redshift-space of haloes with mass below $\sim 7.5 \times 10^{13}h^{-1}M_{\odot}$ at $z=0$ to better than 10%. We demonstrate the technique by applying it to an ensemble of 50 low resolution, large-volume $N$-body simulations to compute the correlation function and covariance matrix of luminous red galaxies (LRGs). The limited resolution of the original simulations results in them resolving just two thirds of the LRG population. We extend the resolution of the simulations by a factor of 30 in halo mass in order to recover all LRGs. With existing simulations it is possible to generate a halo catalogue equivalent to that which would be obtained from a $N$-body simulation using more than 20 trillion particles; a direct simulation of this size is likely to remain unachievable for many years. Using our method it is now feasible to build the large numbers of high-resolution large volume mock galaxy catalogues required to compute the covariance matrices necessary to analyse upcoming galaxy surveys designed to probe dark energy.
• ### Single Colour Diagnostics of the Mass-to-light Ratio: Predictions from Galaxy Formation Models(1304.4421)

April 16, 2013 astro-ph.CO
Accurate galaxy stellar masses are crucial to better understand the physical mechanisms driving the galaxy formation process. We use synthetic star formation and metal enrichment histories predicted by the {\sc galform} galaxy formation model to investigate the precision with which various colours $(m_{a}-m_{b})$ can alone be used as diagnostics of the stellar mass-to-light ratio. As an example, we find that, at $z=0$, the {\em intrinsic} (B$_{f435w}-$V$_{f606w}$) colour can be used to determine the intrinsic rest-frame $V$-band stellar mass-to-light ratio ($\log_{10}\Gamma_{V}=\log_{10}[(M/M_{\odot})/(L_{V}/L_{V\odot})]$) with a precision of $\sigma_{lg\Gamma}\simeq 0.06$ when the initial mass function and redshift are known beforehand. While the presence of dust, assuming a universal attenuation curve, can have a systematic effect on the inferred mass-to-light ratio using a single-colour relation, this is typically small as it is often possible to choose a colour for which the dust reddening vector is approximately aligned with the $(m_{a}-m_{b})-\log_{10}\Gamma_{V}$ relation. The precision with which the stellar mass-to-light ratio can be recovered using a single colour diagnostic rivals implementations of SED fitting using more information but in which simple parameterisations of the star formation and metal enrichment histories are assumed. To facilitate the wide use of these relations, we provide the optimal observer frame colour to estimate the stellar mass-to-light ratio, along with the associated parameters, as a function of redshift ($0<z<1.5$) for two sets of commonly used filters.
• ### The accuracy of the UV continuum as an indicator of the star formation rate in galaxies(1209.3473)

Sept. 16, 2012 astro-ph.CO, astro-ph.IM
The rest-frame intrinsic UV luminosity is often used as an indicator of the instantaneous star formation rate (SFR) in a galaxy. While it is in general a robust indicator of the ongoing star formation activity, the precise value of the calibration relating the UV luminosity to the SFR ($B_{\nu}$), is sensitive to various physical properties, such as the recent star formation and metal enrichment histories, along with the choice of stellar initial mass function. The distribution of these properties for the star-forming galaxy population then suggests that the adoption of a single calibration is not appropriate unless properly qualified with the uncertainties on the calibration. We investigate, with the aid of the {\sc galform} semi-analytic model of galaxy formation, the distribution of UV-SFR calibrations obtained using realistic star formation and metal enrichment histories. At $z=0$, we find that when the initial mass function is fixed (to the Kennicutt IMF), the median calibration is $B_{\rm fuv}=0.9$ where ${\rm SFR}/[{\rm M_{\odot}\,yr^{-1}}]=B_{\nu}\times 10^{-28}\times L_{\nu}/[{\rm ergs\,s^{-1}\,Hz^{-1}}]$. However, the width of the distribution $B_{\rm fuv}$ suggests that for a single object there is around a 20% {\em intrinsic} uncertainty (at $z=0$, rising to $\simeq 30%$ at $z=6$) on the star formation rate inferred from the FUV luminosity without additional constraints on the star formation history or metallicity. We also find that the median value of the calibration $B_{\rm fuv}$ is correlated with the star formation rate and redshift (at $z>3$) raising implications for the correct determination of the star formation rate from the UV.
• ### Predictions for the CO emission of galaxies from a coupled simulation of galaxy formation and photon dominated regions(1204.0795)

Aug. 9, 2012 astro-ph.CO
We combine the galaxy formation model GALFORM with the Photon Dominated Region code UCL_PDR to study the emission from the rotational transitions of 12CO (CO) in galaxies from z=0 to z=6 in the Lambda CDM framework. GALFORM is used to predict the molecular (H2) and atomic hydrogen (HI) gas contents of galaxies using the pressure-based empirical star formation relation of Blitz & Rosolowsky. From the predicted H2 mass and the conditions in the interstellar medium, we estimate the CO emission in the rotational transitions 1-0 to 10-9 by applying the UCL_PDR model to each galaxy. We find that deviations from the Milky-Way CO-H2 conversion factor come mainly from variations in metallicity, and in the average gas and star formation rate surface densities. In the local universe, the model predicts a CO(1-0) luminosity function (LF), CO-to-total infrared (IR) luminosity ratios for multiple CO lines and a CO spectral line energy distribution (SLED) which are in good agreement with observations of luminous and ultra-luminous IR galaxies. At high redshifts, the predicted CO SLED of the brightest IR galaxies reproduces the shape and normalization of the observed CO SLED. The model predicts little evolution in the CO-to-IR luminosity ratio for different CO transitions, in good agreement with observations up to z~5. We use this new hybrid model to explore the potential of using colour selected samples of high-redshift star-forming galaxies to characterise the evolution of the cold gas mass in galaxies through observations with the Atacama Large Millimeter Array.
• ### Predictions for the intrinsic UV continuum properties of star forming galaxies and the implications for inferring dust extinction(1206.2732)

June 13, 2012 astro-ph.CO
The observed ultraviolet continuum (UVC) slope is potentially a powerful diagnostic of dust obscuration in star forming galaxies. However, the intrinsic slope is also sensitive to the form of the stellar initial mass function (IMF) and to the recent star formation and metal enrichment histories of a galaxy. Using the galform semi-analytical model of galaxy formation, we investigate the intrinsic distribution of UVC slopes. For star-forming galaxies, we find that the intrinsic distribution of UVC slopes at z=0, parameterised by the power law index beta, has a standard deviation of sigma_beta=0.30. This suggests an uncertainty on the inferred UV attenuation of A_fuv=0.7\$ (assuming a Calzetti attenuation curve) for an individual object, even with perfect photometry. Furthermore, we find that the intrinsic UVC slope correlates with star formation rate, intrinsic UV luminosity, stellar mass and redshift. These correlations have implications for the interpretation of trends in the observed UVC slope with these quantities irrespective of the sample size or quality of the photometry. Our results suggest that in some cases the attenuation by dust has been incorrectly estimated.
• ### Numerical overcooling in shocks(1106.0306)

June 1, 2011 astro-ph.CO
We present a study of cooling in radiative shocks simulated with smoothed particle hydrodynamics (SPH) and adaptive mesh refinement codes. We obtain a similarity solution for a shock-tube problem in the presence of radiative cooling, and test how well the solution is reproduced in Gadget and Flash. Shock broadening governed by the details of the numerical scheme (artificial viscosity or Riemann solvers) leads to potentially significant overcooling in both codes. We interpret our findings in terms of a resolution criterion, and apply it to realistic simulations of cosmological accretion shocks onto galaxy haloes, cold accretion and thermal feedback from supernovae or active galactic nuclei. To avoid numerical overcooling of accretion shocks onto haloes that should develop a hot corona requires a particle or cell mass resolution of 10^6 M_sun, which is within reach of current state-of-the-art simulations. At this mass resolution, thermal feedback in the interstellar medium of a galaxy requires temperatures of supernova or AGN driven bubbles to be in excess of 10^7 K at densities of n_H=1.0 cm^-3, in order to avoid spurious suppression of the feedback by numerical overcooling.
• ### Which halos host Herschel-ATLAS galaxies in the local Universe?(1011.3048)

Nov. 12, 2010 astro-ph.CO
We measure the projected cross-correlation between low redshift (z < 0.5) far-IR selected galaxies in the SDP field of the Herschel-ATLAS (H-ATLAS) survey and optically selected galaxies from the Galaxy and Mass Assembly (GAMA) redshift survey. In order to obtain robust correlation functions, we restrict the analysis to a subset of 969 out of 6900 H-ATLAS galaxies, which have reliable optical counterparts with r<19.4 mag and well-determined spectroscopic redshifts. The overlap region between the two surveys is 12.6 sq. deg; the matched sample has a median redshift of z ~ 0.2. The cross-correlation of GAMA and H-ATLAS galaxies within this region can be fitted by a power law, with correlation length r_0 ~ 4.63 +/- 0.51 Mpc. Comparing with the corresponding auto-correlation function of GAMA galaxies within the SDP field yields a relative bias (averaged over 2-8 Mpc) of H-ATLAS and GAMA galaxies of b_H/b_G ~ 0.6. Combined with clustering measurements from previous optical studies, this indicates that most of the low redshift H-ATLAS sources are hosted by halos with masses comparable to that of the Milky Way. The correlation function appears to depend on the 250 um luminosity, L_250, with bright (median luminosity \nu L_250 ~ 1.6 x 10^10 L_sun) objects being somewhat more strongly clustered than faint (\nu L_250 ~ 4.0 x 10^9 L_sun) objects. This implies that galaxies with higher dust-obscured star formation rates are hosted by more massive halos.
• ### Designing a space-based galaxy redshift survey to probe dark energy(1006.3517)

July 8, 2010 hep-ph, gr-qc, astro-ph.CO
A space-based galaxy redshift survey would have enormous power in constraining dark energy and testing general relativity, provided that its parameters are suitably optimized. We study viable space-based galaxy redshift surveys, exploring the dependence of the Dark Energy Task Force (DETF) figure-of-merit (FoM) on redshift accuracy, redshift range, survey area, target selection, and forecast method. Fitting formulae are provided for convenience. We also consider the dependence on the information used: the full galaxy power spectrum P(k, P(k) marginalized over its shape, or just the Baryon Acoustic Oscillations (BAO). We find that the inclusion of growth rate information (extracted using redshift space distortion and galaxy clustering amplitude measurements) leads to a factor of ~ 3 improvement in the FoM, assuming general relativity is not modified. This inclusion partially compensates for the loss of information when only the BAO are used to give geometrical constraints, rather than using the full P(k) as a standard ruler. We find that a space-based galaxy redshift survey covering ~20,000 deg^2 over 0.5< z < 2 with \sigma_z/(1+z)<= 0.001 exploits a redshift range that is only easily accessible from space, extends to sufficiently low redshifts to allow both a vast 3-D map of the universe using a single tracer population, and overlaps with ground-based surveys to enable robust modeling of systematic effects. We argue that these parameters are close to their optimal values given current instrumental and practical constraints.
• ### Ram pressure stripping the hot gaseous halos of galaxies in groups and clusters(0710.0964)

Oct. 11, 2007 astro-ph
We use a large suite of carefully controlled full hydrodynamic simulations to study the ram pressure stripping of the hot gaseous halos of galaxies as they fall into massive groups and clusters. The sensitivity of the results to the orbit, total galaxy mass, and galaxy structural properties is explored. For typical structural and orbital parameters, we find that ~30% of the initial hot galactic halo gas can remain in place after 10 Gyr. We propose a physically simple analytic model that describes the stripping seen in the simulations remarkably well. The model is analogous to the original formulation of Gunn & Gott (1972), except that it is appropriate for the case of a spherical (hot) gas distribution (as opposed to a face-on cold disk) and takes into account that stripping is not instantaneous but occurs on a characteristic timescale. The model reproduces the results of the simulations to within approximately 10% at almost all times for all the orbits, mass ratios, and galaxy structural properties we have explored. The one exception involves unlikely systems where the orbit of the galaxy is highly non-radial and its mass exceeds about 10% of the group or cluster into which it is falling (in which case the model under-predicts the stripping following pericentric passage). The proposed model has several interesting applications, including modelling the ram pressure stripping of both observed and cosmologically-simulated galaxies and as a way to improve current semi-analytic models of galaxy formation. One immediate consequence is that the colours and morphologies of satellite galaxies in groups and clusters will differ significantly from those predicted with the standard assumption of complete stripping of the hot coronae.
• ### The metal enrichment of elliptical galaxies in hierarchical galaxy formation models(astro-ph/0504618)

July 14, 2005 astro-ph
We investigate the metal enrichment of elliptical galaxies in the framework of hierarchical models of galaxy formation. The semi-analytical model we use in this Letter, which has been used to study the metal enrichment of the intracluster medium (ICM) by Nagashima et al., includes the effects of flows of gas and metals both into and out of galaxies and the processes of metal enrichment due to both type Ia and type II supernovae. We adopt a solar neighbourhood IMF for star formation in discs, but consider models in which starbursts have either a solar neighbourhood IMF or a top-heavy IMF. We find that the alpha-element abundance in ellipticals is consistent with observed values only if the top-heavy IMF is used. This result is consistent with our previous study on the metal enrichment of the ICM. We also discuss the abundance ratio of alpha elements to iron as a function of velocity dispersion and metallicity. We find that models with a top-heavy IMF match the alpha/Fe ratios observed in typical L_* ellipticals, but none of the models reproduce the observed increase of alpha/Fe with velocity dispersion.