GALICS: capturing the panchromaticity of galaxies

This contribution describes results obtained with the GALICS model (for Galaxies In Cosmological Simulations), which is a hybrid model for hierarchical galaxy formation studies, combining the outputs of large cosmological N-body simulations with simple, semi-analytic recipes to describe the fate of the baryons within dark matter halos. Designed to predict the overall statistical properties of galaxies, with special emphasis on the panchromatic spectral energy distribution emitted by galaxies in the UV/optical and IR/submm wavelength ranges, such an approach can be used to predict the galaxy luminosity function evolution from the ultraviolet to far infrared, along with individual galaxies star formation histories. This revised version was published online in September 2006 with corrections to the Cover Date.     

Young Stellar Populations at High z: Where are They?

Our view on the deep universe is strongly biassed towards the starlight that directly escapes from high-redshift galaxies, since we know very little on the fraction of luminosity absorbed by dust. Attempts to correct for this effect directly from the slope of the UV spectra seem to suggest that a significant fraction of the UV flux is extinguished. New constraints are now set on the fraction of dust-enshrouded young stellar populations by the detection of the Cosmic Infrared Background, and faint galaxy counts at IR and submm wavelengths. We briefly review the observations and use a semi-analytic model of galaxy formation and evolution to predict number counts consistent with the background. This revised version was published online in July 2006 with corrections to the Cover Date.     

The clustering evolution of the galaxy distribution

We follow the evolution of the galaxy population in a ΛCDM cosmology by means of high-resolution N -body simulations in which the formation of galaxies and their observable properties are calculated using a semi-analytic model. We display images of the spatial distribution of galaxies in the simulations that illustrate its evolution and provide a qualitative understanding of the processes responsible for the various biases that develop. We consider three specific statistical measures of clustering at     and     : the correlation length (in both real and redshift space) of galaxies of different luminosity, the morphology–density relation and the genus curve of the topology of galaxy isodensity surfaces. For galaxies with luminosity below L ∗, the     correlation length depends very little on the luminosity of the sample, but for brighter galaxies it increases very rapidly, reaching values in excess of 10  h ⁻¹ Mpc. The 'accelerated' dynamical evolution experienced by galaxies in rich clusters, which is partly responsible for this effect, also results in a strong morphology–density relation. Remarkably, this relation is already well-established at     . The genus curves of the galaxies are significantly different from the genus curves of the dark matter, however this is not a result of genuine topological differences but rather of the sparse sampling of the density field provided by galaxies. The predictions of our model at     will be tested by forthcoming data from the 2dF and Sloan galaxy surveys, and those at     by the DEEP and VIRMOS surveys.     

Galaxy–dark matter correlations applied to galaxy–galaxy lensing: predictions from the semi-analytic galaxy formation models

We use semi-analytic models of galaxy formation combined with high-resolution N -body simulations to make predictions for galaxy–dark matter correlations and apply them to galaxy–galaxy lensing. We analyse cross-power spectra between the dark matter and different galaxy samples selected by luminosity, colour or star formation rate. We compare the predictions with the recent detection by the Sloan Digital Sky Survey (SDSS). We show that the correlation amplitude and the mean tangential shear depend strongly on the luminosity of the sample on scales below 1  h ⁻¹ Mpc, reflecting the correlation between the galaxy luminosity and the halo mass. The cross-correlation cannot, however, be used to infer the halo profile directly because different halo masses dominate on different scales and because not all galaxies are at the centres of the corresponding haloes. We compute the redshift evolution of the cross-correlation amplitude and compare it with those of galaxies and dark matter. We also compute the galaxy–dark matter correlation coefficient and show that it is close to unity on scales above 1  h ⁻¹ Mpc for all considered galaxy types. This would allow one to extract the bias and the dark matter power spectrum on large scales from the galaxy and galaxy–dark matter correlations.     

Constraints on σ8 from galaxy clustering in N-body simulations and semi-analytic models

We generate mock galaxy catalogues for a grid of different cosmologies, using rescaled N -body simulations in tandem with a semi-analytic model run using consistent parameters. Because we predict the galaxy bias, rather than fitting it as a nuisance parameter, we obtain an almost pure constraint on σ₈ by comparing the projected two-point correlation function we obtain to that from the Sloan Digital Sky Survey (SDSS). A systematic error arises because different semi-analytic modelling assumptions allow us to fit the r -band luminosity function equally well. Combining our estimate of the error from this source with the statistical error, we find  σ₈= 0.97 ± 0.06  . We obtain consistent results if we use galaxy samples with a different magnitude threshold, or if we select galaxies by b _J-band rather than r -band luminosity and compare to data from the 2dF Galaxy Redshift Survey (2dFGRS). Our estimate for σ₈ is higher than that obtained for other analyses of galaxy data alone, and we attempt to find the source of this difference. We note that in any case, galaxy clustering data provide a very stringent constraint on galaxy formation models.     

Simulating wide-field quasar surveys from the optical to near-infrared

A number of deep, wide-field, near-infrared (NIR) surveys employing new infrared cameras on 4-m class telescopes are about to commence. These surveys have the potential to determine the fraction of luminous dust-obscured quasars that may have eluded surveys undertaken at optical wavelengths. In order to understand the new observations it is essential to make accurate predictions of surface densities and number–redshift relations for unobscured quasars in the NIR based on information from surveys at shorter wavelengths. The accuracy of the predictions depends critically on a number of key components. The commonly used single power-law representation for quasar spectral energy distributions (SEDs) is inadequate and the use of an SED incorporating the upturn in continuum flux at  λ∼ 12 000 Å  is essential. The presence of quasar host galaxies is particularly important over the rest-frame wavelength interval  8000 < λ < 16 000 Å  and we provide an empirical determination of the magnitude distribution of host galaxies using a low-redshift sample of quasars from the Sloan Digital Sky Survey Data Release 3 quasar catalogue. A range of models for the dependence of host galaxy luminosity on quasar luminosity is investigated, along with the implications for the NIR surveys. Even adopting a conservative model for the behaviour of host galaxy luminosity the number counts for shallow surveys in the K band increase by a factor of 2. The degree of morphological selection applied to define candidate quasar samples in the NIR is found to be an important factor in determining the fraction of the quasar population included in such samples.     

The effect of dwarf galaxy disruption in semi-analytic models

We present results for a galaxy formation model that includes a simple treatment for the disruption of dwarf galaxies by gravitational forces and galaxy encounters within galaxy clusters. This is implemented a posteriori in a semi-analytic model by considering the stability of cluster dark matter subhaloes at   z = 0  . We assume that a galaxy whose dark matter substructure has been disrupted will itself disperse, while its stars become part of the population of intracluster stars responsible for the observed intracluster light. Despite the simplicity of this assumption, our results show a substantial improvement over previous models and indicate that the inclusion of galaxy disruption is indeed a necessary ingredient of galaxy formation models. We find that galaxy disruption suppresses the number density of dwarf galaxies by about a factor of 2. This makes the slope of the faint end of the galaxy luminosity function shallower, in agreement with observations. In particular, the abundance of faint, red galaxies is strongly suppressed. As a result, the luminosity function of red galaxies and the distinction between the red and the blue galaxy populations in colour–magnitude relationships are correctly predicted. Finally, we estimate a fraction of intracluster light comparable to that found in clusters of galaxies.     

Signatures of recent star formation in ring S0 galaxies

We present a study of the stellar populations of ring and/or arm-like structures in a sample of S0 galaxies using GALEX far- and near-ultraviolet imaging and SDSS optical data. Such structures are prominent in the UV and reveal recent star formation. We quantitatively characterize these rejuvenation events, estimating the average age and stellar mass of the ring structures, as well as of the entire galaxy. The mass fraction of the UV-bright rings is a few percent of the total galaxy mass, although the UV ring luminosity reaches 70% of the galaxy luminosity. The integrated colors of these S0s locates them in the red sequence (NGC 2962) and in the so-called green valley. We suggest that the star formation episodes may be induced by different triggering mechanisms, such as the inner secular evolution driven by bars, and interaction episodes.     

A photometric study of the field around the candidate recoiling/binary black hole SDSS J092712.65+294344.0

We present a photometric far-ultraviolet (FUV) to K _s-band study of the field around quasar SDSS J092712.65+294344.0. The SDSS spectrum of this object shows various emission lines with two distinct redshifts, at   z = 0.699  and 0.712. Because of this peculiar spectroscopic feature, this source has been proposed as a candidate recoiling or binary black hole. A third alternative model involves two galaxies moving in the centre of a rich galaxy cluster. Here, we present a study addressing the possible presence of such a rich cluster of galaxies in the SDSS J092712.65+294344.0 field. We observed the  3.6 × 2.6  arcmin² field in the K _s band and matched the near-infrared data with the FUV and near-ultraviolet images in the Galaxy Evolution Explorer archive and the ugriz observations in the SDSS. From various colour–colour diagrams, we were able to classify the nature of 32 sources, only 6–11 of which have colours consistent with galaxies at   z ≈ 0.7  . We compare these numbers with the surface density of galaxies, stars and quasars and the expectations for typical galaxy clusters both at low and high redshift. Our study shows that the galaxy cluster scenario is in clear disagreement with the new observations.     

The K-band Hubble diagram for the brightest cluster galaxies: a test of hierarchical galaxy formation models

We analyse the K -band Hubble diagram for a sample of brightest cluster galaxies (BCGs) in the redshift range 0< z <1. In good agreement with earlier studies, we confirm that the scatter in the absolute magnitudes of the galaxies is small (0.3 mag). The BCGs exhibit very little luminosity evolution in this redshift range: if q ₀=0.0, we detect no luminosity evolution; for q ₀=0.5, we measure a small negative evolution (i.e., BCGs were about 0.5 mag fainter at z =1 than today). If the mass in stars of these galaxies had remained constant over this period of time, substantial positive luminosity evolution would be expected: BCGs should have been brighter in the past, since their stars were younger. A likely explanation for the observed zero or negative evolution is that the stellar mass of the BCGs has been assembled over time through merging and accretion, as expected in hierarchical models of galaxy formation. The colour evolution of the BCGs is consistent with that of an old stellar population ( z _for>2) that is evolving passively. We can thus use evolutionary population synthesis models to estimate the rate of growth in stellar mass for these systems. We find that the stellar mass in a typical BCG has grown by a factor ≃2 since z ≃1 if q ₀=0.0, or by factor ≃4 if q ₀=0.5. These results are in good agreement with the predictions of semi-analytic models of galaxy formation and evolution set in the context of a hierarchical scenario for structure formation. The models predict a scatter in the luminosities of the BCGs that is somewhat larger than the observed one, but that depends on the criterion used to select the model clusters.     

Hα and UV luminosities and star formation rates in a large sample of luminous compact galaxies

We present the results of a statistical study of the star formation rates (SFR) derived from the Galaxy Evolution Explorer (GALEX) observations in the ultraviolet continuum and in the Hα emission line for a sample of about 800 luminous compact galaxies (LCGs). Galaxies in this sample have a compact structure and include one or several regions of active star formation. Global galaxy characteristics (metallicity, luminosity, stellar mass) are intermediate between ones of the nearby blue compact dwarf (BCD) galaxies and Lyman-break galaxies (LBGs) at high redshifts z>2–3. SFRs were corrected for interstellar extinction which was derived from the optical Sloan Digital Sky Survey (SDSS) spectra. We find that SFRs derived from the galaxy luminosities in the far ultraviolet (FUV) and near ultraviolet (NUV) ranges vary in a wide range from 0.18 M _⊙?yr^?1 to 113 M _⊙?yr^?1 with median values of 3.8 M _⊙?yr^?1 and 5.2 M _⊙?yr^?1, respectively. Simple regression relations are found for luminosities L(Hα) and L(UV) as functions of the mass of the young stellar population, the starburst age, and the galaxy metallicity. We consider the evolution of L(Hα), L(FUV) and L(NUV) with a starburst age and introduce new characteristics of star formation, namely the initial Hα, FUV and NUV luminosities at zero starburst age.     

Optical and infrared signatures of ultra-luminous X-ray sources

We have constructed a model to describe the optical emission from ultra-luminous X-ray sources (ULXs). We assume a binary model with a black hole accreting matter from a Roche lobe filling companion star. We consider the effects of radiative transport and radiative equilibrium in the irradiated surfaces of both the star and a thin accretion disc. We have developed this model as a tool with which to positively identify the optical counterparts of ULXs, and subsequently derive parameters such as the black hole mass and the luminosity class and spectral type of the counterpart. We examine the dependence of the optical emission on these and other variables. We extend our model to examine the magnitude variation at infrared wavelengths, and we find that observations at these wavelengths may have more diagnostic power than in the optical. We apply our model to existing HST observations of the candidates for the optical counterpart of ULX X-7 in NGC 4559. All candidates could be consistent with an irradiated star alone, but we find that a number of them are too faint to fit with an irradiated star and disc together. Were one of these the optical counterpart to X-7, it would display a significant temporal variation.     

Evolutionary synthesis models for galaxy transformation in clusters

The galaxy population in rich local galaxy clusters shows a ratio of one quarter elliptical galaxies, two quarters S0 galaxies, and one quarter spiral galaxies. Observations of clusters at redshift 0.5 show a perspicuously different ratio, the dominant galaxy type are spiral galaxies with a fraction of two quarters while the number of S0 galaxies decreases to a fraction of one quarter (Dressler et al. 1997). This shows an evolution of the galaxy population in clusters since redshift 0.5 and it has been suspected that galaxy transformation processes during the infall into a cluster are responsible for this change. These could be merging, starburst or ram-pressure stripping. We use our evolutionary synthesis models to describe various possible effects of those interactions on the star formation of spiral galaxies infalling into clusters. We study the effects of starbursts of various strengths as well as of the truncation of star formation at various epochs on the color and luminosity evolution of model galaxies of various spectral types. As a first application we present the comparison of our models with observed properties of the local S0 galaxy population to constrain possible S0 formation mechanisms in clusters. Application to other types of galaxies is planned for the future. This revised version was published online in August 2006 with corrections to the Cover Date.     

Chemo-spectrophotometric evolution of spiral galaxies – III. Abundance and colour gradients in discs

We study the relations between luminosity and chemical-abundance profiles of spiral galaxies, using detailed models for the chemical and spectrophotometric evolution of galactic discs. The models are 'calibrated' on the Milky Way disc and are successfully extended to other discs with the help of simple 'scaling' relations, obtained in the framework of semi-analytic models of galaxy formation. We find that our models exhibit oxygen abundance gradients that increase in absolute value with decreasing disc luminosity (when expressed in dex kpc⁻¹) and are independent of disc luminosity (when expressed in dex scalelength⁻¹), both in agreement with observations. We notice an important strong correlation between abundance gradient and disc scalelength. These results support the idea of 'homologous evolution' of galactic discs.     

Clustering of galaxies in a hierarchical universe – III. Mock redshift surveys

This is the third paper in a series which combines N -body simulations and semi-analytic modelling to provide a fully spatially resolved simulation of the galaxy formation and clustering processes. Here we extract mock redshift surveys from our simulations: a cold dark matter model with either Ω₀=1 ( τ CDM) or Ω₀=0.3 and Λ=0.7 (ΛCDM). We compare the mock catalogues with the northern region (CfA2N) of the Center for Astrophysics (CfA) Redshift Surveys. We study the properties of galaxy groups and clusters identified using standard observational techniques, and also the relation of these groups to real virialized systems. Most features of CfA2N groups are reproduced quite well by both models with no obvious dependence on Ω₀. Redshift‐space correlations and pairwise velocities are also similar in the two cosmologies. The luminosity functions predicted by our galaxy formation models depend sensitively on the treatment of star formation and feedback. For the particular choices of Paper I they agree poorly with the CfA survey. To isolate the effect of this discrepancy on our mock redshift surveys, we modify galaxy luminosities in our simulations to reproduce the CfA luminosity function exactly. This adjustment improves agreement with the observed abundance of groups, which depends primarily on the galaxy luminosity density, but other statistics, connected more closely with the underlying mass distribution, remain unaffected. Regardless of the luminosity function adopted, modest differences with observation remain. These can be attributed to the presence of the 'Great Wall' in the CfA2N. It is unclear whether the greater coherence of the real structure is a result of cosmic variance, given the relatively small region studied, or reflects a physical deficiency of the models.