Search

We investigate the performance of group finding algorithms that reconstruct galaxy groups from the positional information of tracer galaxies that are observed in redshift surveys carried out with multiplexed spectrographs. We use mock light-cones produced by the L-Galaxies semi-analytic model of galaxy evolution in which the underlying reality is known. We particularly focus on the performance at high redshift, and how this is affected by choices of the mass of the tracer galaxies (largely equivalent to their co-moving number density) and the (assumed random) sampling rate of these tracers. We first however compare two different approaches to group finding as applied at low redshift, and conclude that these are broadly comparable. For simplicity we adopt just one of these, "Friends-of-Friends" (FoF) as the basis for our study at high redshift. We introduce 12 science metrics that are designed to quantify the performance of the group-finder as relevant for a wide range of science investigations with a group catalogue. These metrics examine the quality of the recovered group catalogue, the median halo masses of different richness structures, the scatter in dark matter halo mass and how successful the group-finder classifies singletons, centrals and satellites. We analyze how these metrics vary with the limiting stellar mass and random sampling rate of the tracer galaxies, allowing quantification of the various trade-offs between different possible survey designs. Finally, we look at the impact of these same design parameters on the relative "costs" in observation time of the survey using as an example the potential MOONRISE survey using the MOONS instrument., ERC STFC

We compare predictions for galaxy-galaxy lensing profiles and clustering from the Henriques et al. (2015) public version of the Munich semi-analytical model of galaxy formation (SAM) and the IllustrisTNG suite, primarily TNG300, with observations from KiDS+GAMA and SDSS-DR7 using four different selection functions for the lenses (stellar mass, stellar mass and group membership, stellar mass and isolation criteria, stellar mass and colour). We find that this version of the SAM does not agree well with the current data for stellar mass-only lenses with $M_\ast > 10^{11}\,M_\odot$. By decreasing the merger time for satellite galaxies as well as reducing the radio-mode AGN accretion efficiency in the SAM, we obtain better agreement, both for the lensing and the clustering, at the high mass end. We show that the new model is consistent with the signals for central galaxies presented in Velliscig et al. (2017). Turning to the hydrodynamical simulation, TNG300 produces good lensing predictions, both for stellar mass-only ($��^2 = 1.81$ compared to $��^2 = 7.79$ for the SAM), and locally brightest galaxies samples ($��^2 = 3.80$ compared to $��^2 = 5.01$). With added dust corrections to the colours it matches the SDSS clustering signal well for red low mass galaxies. We find that both the SAMs and TNG300 predict $\sim 50\,\%$ excessive lensing signals for intermediate mass red galaxies with $10.2 < \log_{10} M_\ast [ M_\odot ] < 11.2$ at $r \approx 0.6\,h^{-1}\,\mathrm{Mpc}$, which require further theoretical development., 30 pages, 20 figures, 10 tables, minor revisions, matches version accepted by MNRAS

We develop a conceptual framework and methodology to study the drivers of the quenching of galaxies, including the drivers of galactic conformity. The framework is centred on the statistic $\Delta$, which is defined as the difference between the observed star-formation state of a galaxy, and a prediction of its state based on an empirical model of quenching. In particular, this work uses the average quenching effects of stellar mass and local density to construct an empirical model of quenching. $\Delta$ is therefore a residual which reflects the effects of drivers of quenching not captured by stellar mass and local density, or so-called 'hidden variables'. Through a toy model, we explore how the statistical properties of $\Delta$ can be used to learn about the internal and external hidden variables which control the quenching of a sample of galaxies. We then apply this analysis to a sample of local galaxies and find that, after accounting for the average quenching effects of stellar mass and local density, $\Delta$ remains correlated out to separations of 3 Mpc. Furthermore, we find that external hidden variables remain important for driving the residual quenching of low-mass galaxies, while the residual quenching of high-mass galaxies is driven mostly by internal properties. These results, along with a similar analysis of a semi-analytical mock catalogue, suggest that it is necessary to consider halo-related properties as candidates for hidden variables. A preliminary halo-based analysis indicates that much of the correlation of $\Delta$ can be attributed to the physics associated with individual haloes., Comment: 19 pages, 11 figures, submitted to MNRAS

We present a modified version of the L-GALAXIES 2020 semi-analytic model of galaxy evolution, which includes significantly increased direct metal enrichment of the circumgalactic medium (CGM) by supernovae (SNe). These more metal-rich outflows do not require increased mass-loading factors, in contrast to some other galaxy evolution models. This modified L-GALAXIES 2020 model is able to simultaneously reproduce the gas-phase metallicity $(Z_{\rm g})$ and stellar metallicity $(Z_{*})$ radial profiles observed in nearby disc galaxies by MaNGA and MUSE, as well as the observed mass - metallicity relations for gas and stars at $z=0$ and their evolution back to $z\sim{}2-3$. A direct CGM enrichment fraction of $\sim{}90\%$ for SNe-II is preferred. We find that massive disc galaxies have slightly flatter $Z_{\rm g}$ profiles than their lower-mass counterparts in L-GALAXIES 2020, due to more efficient enrichment of their outskirts via inside-out growth and metal-rich accretion. Such a weak, positive correlation between stellar mass and $Z_{\rm g}$ profile slope is also seen in our MaNGA-DR15 sample of 571 star-forming disc galaxies. Although, below ${\rm log}(M_{*}/{\rm M}_{\odot})\sim{}10.0$ this observational result is strongly dependent on the metallicity diagnostic and morphological selection chosen. In addition, a lowered maximum SN-II progenitor mass of $25{\rm M}_{\odot}$, reflecting recent theoretical and observational estimates, can also provide a good match to observed metallicity profiles at $z=0$ in L-GALAXIES 2020. However, this model version fails to reproduce an evolution in $Z_{\rm g}$ at fixed mass over cosmic time, or the magnesium abundances observed in the intracluster medium (ICM)., Comment: 22 pages, 19 figures. Accepted by MNRAS. (Missing citation hyperlinks corrected)

We study the quenching of star formation as a function of redshift, environment and stellar mass in the galaxy formation simulations of Henriques et al. (2015), which implement an updated version of the Munich semi-analytic model (L-GALAXIES) on the two Millennium Simulations after scaling to a Planck cosmology. In this model, massive galaxies are quenched by active galactic nucleus (AGN) feedback depending on both black hole and hot gas mass, and hence indirectly on stellar mass. In addition, satellite galaxies of any mass can be quenched by ram-pressure or tidal stripping of gas and through the suppression of gaseous infall. This combination of processes produces quenching efficiencies which depend on stellar mass, host halo mass, environment density, distance to group centre and group central galaxy properties in ways which agree qualitatively with observation. Some discrepancies remain in dense regions and close to group centres, where quenching still seems too efficient. In addition, although the mean stellar age of massive galaxies agrees with observation, the assumed AGN feedback model allows too much ongoing star formation at late times. The fact that both AGN feedback and environmental effects are stronger in higher density environments leads to a correlation between the quenching of central and satellite galaxies which roughly reproduces observed conformity trends inside haloes., Monthly Notices of the Royal Astronomical Society, 469 (3), ISSN:0035-8711, ISSN:1365-2966, ISSN:1365-8711

We have updated the Munich galaxy formation model, L-Galaxies, to follow the radial distributions of stars and atomic and molecular gas in galaxy discs. We include an H2-based star-formation law, as well as a detailed chemical-enrichment model with explicit mass-dependent delay times for SN-II, SN-Ia and AGB stars. Information about the star formation, feedback and chemical-enrichment histories of discs is stored in 12 concentric rings. The new model retains the success of its predecessor in reproducing the observed evolution of the galaxy population, in particular, stellar mass functions and passive fractions over the redshift range 0, published in MNRAS. A full description of the model is available as MNRAS online supplementary material and at https://lgalaxiespublicrelease.github.io/. Data from the model are publicly released in SQL-queryable form and in flat files at http://www.mpa-garching.mpg.de/millennium

In this work we update the L-Galaxies semi-analytic model (SAM) to better follow the physical processes responsible for the growth of bulges via disc instabilities (leading to pseudo-bulges) and mergers (leading to classical bulges). We address the former by considering the contribution of both stellar and gaseous discs in the stability of the galaxy, and we update the latter by including dissipation of energy in gas-rich mergers. Furthermore, we introduce angular momentum losses during cooling and find that an accurate match to the observed correlation between stellar disc scale length and mass at z ∼ 0.0 requires that the gas loses 20% of its initial specific angular momentum to the corresponding dark matter halo during the formation of the cold gas disc. We reproduce the observed trends between the stellar mass and specific angular momentum for both disc- and bulge-dominated galaxies, with the former rotating faster than the latter of the same mass. We conclude that a two-component instability recipe provides a morphologically diverse galaxy sample which matches the observed fractional breakdown of galaxies into different morphological types. This recipe also enables us to obtain an excellent fit to the morphology-mass relation and stellar mass function of different galactic types. Finally, we find that energy dissipation during mergers reduces the merger remnant sizes and allows us to match the observed mass-size relation for bulge-dominated systems.

We implement a detailed dust model into the L-Galaxies semi-analytical model which includes: injection of dust by type II and type Ia supernovae (SNe) and AGB stars; grain growth in molecular clouds; and destruction due to supernova-induced shocks, star formation, and reheating. Our grain growth model follows the dust content in molecular clouds and the inter-cloud medium separately, and allows growth only on pre-existing dust grains. At early times, this can make a significant difference to the dust growth rate. Above $z\sim8$, type II SNe are the primary source of dust, whereas below $z\sim8$, grain growth in molecular clouds dominates, with the total dust content being dominated by the latter below $z\sim6$. However, the detailed history of galaxy formation is important for determining the dust content of any individual galaxy. We introduce a fit to the dust-to-metal (DTM) ratio as a function of metallicity and age, which can be used to deduce the DTM ratio of galaxies at any redshift. At $z\lesssim3$, we find a fairly flat mean relation between metallicity and the DTM, and a positive correlation between metallicity and the dust-to-gas (DTG) ratio, in good agreement with the shape and normalisation of the observed relations. We also match the normalisation of the observed stellar mass -- dust mass relation over the redshift range of $0-4$, and to the dust mass function at $z=0$. Our results are important in interpreting observations on the dust content of galaxies across cosmic time, particularly so at high redshift., 19 pages, 16 figures. Edits: Corrected equation 15 and its parameters. Accepted in MNRAS

We present an exploration of the expected detection of the earliest Active Galactic Nuclei (AGN) in the Universe from state-of-art galaxy formation and evolution semi-analytic models and hydro-dynamical simulations. We estimate the number and radiative characteristics of Super Massive Black Holes (SMBHs) at $z\geq 6$, a redshift range that will be intensively explored by the next generation of telescopes, in particular in the radio through the Square Kilometre Array (SKA) and at high energies with ESA's Athena X-ray Observatory. We find that Athena will be able to observe over 5000 AGN/$\rm deg^2$ at the Epoch of Re-ionization (EoR), $6\leq z \leq 10$. Similarly, for the same redshift range the models/simulations suggest that SKA will detect at least 400 AGN/$\rm deg^2$. Additionally, we stress the importance of the volume of the simulation box as well as the initial physical conditions of the models/simulations on their effect on the luminosity functions (LFs) and the creation of the most massive SMBHs that we currently observe at the EoR. Furthermore, following the evolution of the accretion mode of the SMBHs in each model/simulation, we show that, while the quasar dominates over the radio mode at the EoR, detection at radio wavelengths still reaches significant numbers even at the highest redshifts. Finally, we present the effect that the radiative efficiency has on the LFs by comparing results produced with a constant value for the radiative efficiency and more complex calculations based on the spin of each SMBH., Comment: Accepted for publication in MNRAS, 17 pages, 8 figures, 2 tables

We present the results of a search for galaxy clusters and groups in the $\sim2$ square degree of the COSMOS field using all available X-ray observations from the XMM-Newton and Chandra observatories. We reach an X-ray flux limit of $3\times10^{-16}\;ergs\;cm^{-2}\;s^{-1}$ in 0.5--2 keV range, and identify 247 X-ray groups with $M_{200c}=8\times10^{12}-3\times10^{14}\;M_{\odot}$ at a redshift range of $0.08\leq z, Accepted 2018 November 16 (MNRAS)

We study the cosmological build-up of pseudobulges using the LGalaxies semi-analytical model for galaxy formation with a new approach for following separately the assembly of classical bulges and pseudobulges. Classical bulges are assumed to be the result of violent processes (i.e., mergers and starbursts), while the formation of pseudobulges is connected to the secular growth of disks. We apply the model to both the Millennium and the Millennium II simulations, in order to study our results across a wide range of stellar masses ($10^{7} - 10^{11.5} M_{\odot}$). We find that $z=0$ pseudobulges mainly reside in galaxies of $ \rm M_{stellar} \sim 10^{10} - 10^{10.5} M_{\odot}$ ($\rm M_{halo} \sim 10^{11.5}-10^{12} M_{\odot}$) and we recover structural properties of these objects (e.g., sizes and bulge-to-total ratios) that are in good agreement with observational results. Tracing their formation history, we find that pseudobulges assembled in galaxies with a very quiet merger history, as opposed to the host galaxies of classical bulges. Regarding the bulge structure, we find that $\sim$ 30\% of the galaxies with a predominant pseudobulge feature a composite structure, hosting both a pseudo and a classical bulge component. The classical component typically constitutes $\sim$10\% of the total bulge galaxy mass. When looking at the properties of the host galaxies, we find that $z = 0$ pseudobulges are hosted by main sequence galaxies, characterized by a stellar population which is generally younger compared to the one of the hosts of classical bulges., Comment: Submitted to MNRAS

We investigate the dependence of galaxy structure on a variety of galactic and environmental parameters for ~500,000 galaxies at z, Accepted to MNRAS. 31 pages, 15 figures

We use the Henriques et al. (2015) version of the Munich galaxy formation model (L-GALAXIES) to investigate why the halo and stellar mass scales above which galaxies are quenched are constant with redshift and coincide with the scale where baryons are most efficiently converted into stars. This model assumes that central galaxies are quenched by AGN feedback when hot halo gas accretes onto a supermassive black hole. Nevertheless, we find that supernova (SN) feedback sets both mass scales. As haloes grow above a threshold mass, SNe can no longer eject material so their hot gas content increases, enhancing the cooling rate onto the central galaxy, its cold gas content, its star formation rate and the growth rate of its central black hole. Strong AGN feedback terminates this short-lived phase by suppressing the fuel supply for star formation. Despite strong evolution of the halo mass - temperature relation, quenching occurs at a redshift-independent halo and stellar mass which coincides with the mass where baryons have been converted into stars with maximal efficiency. These regularities and coincidences are a result of the specific parameters selected by MCMC tuning of the model to fit the observed abundance and passive fraction of galaxies over the redshift range 0, Comment: 12 pages, 8 figures, submitted to MNRAS

Using the L-Galaxies semi-analytic model we simultaneously fit the H I mass function, stellar mass function and the fraction of red galaxies. We find good fits to all three observations at z = 0 and to the stellar mass function and red fraction at z = 2. Using Markov Chain Monte Carlo (MCMC) techniques we adjust the L-Galaxies parameters to best fit the constraining data. In order to fit the H I mass function we must greatly reduce the gas surface density threshold for star formation, thus lowering the number of low H I mass galaxies. A simultaneous reduction in the star formation efficiency prevents the overproduction of stellar content. A simplified model in which the surface density threshold is eliminated altogether also provides a good fit to the data. Unfortunately, these changes weaken the fit to the Kennicutt–Schmidt relation and raise the star formation rate density at recent times, suggesting that a change to the model is required to prevent accumulation of gas on to dwarf galaxies in the local Universe., Monthly Notices of the Royal Astronomical Society, 472 (2), ISSN:0035-8711, ISSN:1365-2966, ISSN:1365-8711

We present a relationship between the black hole mass, stellar mass, and star formation rate of a diverse group of 91 galaxies with dynamically-measured black hole masses. For our sample of galaxies with a variety of morphologies and other galactic properties, we find that the specific star formation rate is a smoothly decreasing function of the ratio between black hole mass and stellar mass, or what we call the specific black hole mass. In order to explain this relation, we propose a physical framework where the gradual suppression of a galaxy's star formation activity results from the adjustment to an increase in specific black hole mass and, accordingly, an increase in the amount of heating. From this framework, it follows that at least some galaxies with intermediate specific black hole masses are in a steady state of partial quiescence with intermediate specific star formation rates, implying that both transitioning and steady-state galaxies live within this region known as the "green valley." With respect to galaxy formation models, our results present an important diagnostic with which to test various prescriptions of black hole feedback and its effects on star formation activity., 15 pages, 4 figures, 2 tables. Accepted for publication in The Astrophysical Journal

We present a clustering comparison of 12 galaxy formation models [including semi-analytic models (SAMs) and halo occupation distribution (HOD) models] all run on halo catalogues and merger trees extracted from a single Λ cold dark matter N-body simulation. We compare the results of the measurements of the mean halo occupation numbers, the radial distribution of galaxies in haloes and the two-point correlation functions (2PCF). We also study the implications of the different treatments of orphan (galaxies not assigned to any dark matter subhalo) and non-orphan galaxies in these measurements. Our main result is that the galaxy formation models generally agree in their clustering predictions but they disagree significantly between HOD and SAMs for the orphan satellites. Although there is a very good agreement between the models on the 2PCF of central galaxies, the scatter between the models when orphan satellites are included can be larger than a factor of 2 for scales smaller than 1 h(−1) Mpc. We also show that galaxy formation models that do not include orphan satellite galaxies have a significantly lower 2PCF on small scales, consistent with previous studies. Finally, we show that the 2PCF of orphan satellites is remarkably different between SAMs and HOD models. Orphan satellites in SAMs present a higher clustering than in HOD models because they tend to occupy more massive haloes. We conclude that orphan satellites have an important role on galaxy clustering and they are the main cause of the differences in the clustering between HOD models and SAMs., Monthly Notices of the Royal Astronomical Society, 469 (1), ISSN:0035-8711, ISSN:1365-2966, ISSN:1365-8711

We investigate the effects of galaxy formation on the baryonic acoustic oscillations (BAO) peak by applying semi-analytic modelling techniques to the Millennium-XXL, a 3×10 11 particle N-body simulation of similar volume to the future EUCLID survey. Our approach explicitly incorporates the effects of tidal fields and stochasticity on halo formation, as well as the presence of velocity bias, spatially correlated merger histories, and the connection of all these with the observable and physical properties of galaxies. We measure significant deviations in the shape of the BAO peak from the expectations of a linear bias model built on top of the nonlinear dark matter distribution. We find that the galaxy correlation function shows an excess close to the maximum of the BAO peak (r � 110h 1 Mpc) and a deficit at r � 90h 1 Mpc. Depending on the redshift, selection criteria and number density of the galaxy samples, these bias distortions can be up to 5% in amplitude. They are, however, largely absorbed by marginalization over nuisance parameters in current analytical modelling of the BAO peak in configuration space, in particular into the parameter that controls the broadening due to nonlinear evolution. As a result, the galaxy formation effects detected here are unlikely to bias the high-precision measurements planned by the upcoming generation of wide-field galaxy surveys.

We use the SDSS/DR8 galaxy sample to study the radial distribution of satellite galaxies around isolated primaries, comparing to semi-analytic models of galaxy formation based on the Millennium and Millennium-II simulations. SDSS satellites behave differently around high- and low-mass primaries: those orbiting objects with $M_*>10^{11}M_\odot$ are mostly red and are less concentrated towards their host than the inferred dark matter halo, an effect that is very pronounced for the few blue satellites. On the other hand, less massive primaries have steeper satellite profiles that agree quite well with the expected dark matter distribution and are dominated by blue satellites, even in the inner regions where strong environmental effects are expected. In fact, such effects appear to be strong only for primaries with $M_* > 10^{11}M_\odot$. This behaviour is not reproduced by current semi-analytic simulations, where satellite profiles always parallel those of the dark matter and satellite populations are predominantly red for primaries of all masses. The disagreement with SDSS suggests that environmental effects are too efficient in the models. Modifying the treatment of environmental and star formation processes can substantially increase the fraction of blue satellites, but their radial distribution remains significantly shallower than observed. It seems that most satellites of low-mass primaries can continue to form stars even after orbiting within their joint halo for 5 Gyr or more., 19 pages, 11 figures, accepted by MNRAS

We calculate stellar masses for massive luminous galaxies at redshift 0.2-0.7 using the first two years of data from the Baryon Oscillation Spectroscopic Survey (BOSS). Stellar masses are obtained by fitting model spectral energy distributions to u,g,r,i,z magnitudes, and simulations with mock galaxies are used to understand how well the templates recover the stellar mass. Accurate BOSS spectroscopic redshifts are used to constrain the fits. We find that the distribution of stellar masses in BOSS is narrow (Delta log M~0.5 dex) and peaks at about logM ~ 11.3 (for a Kroupa initial stellar mass function), and that the mass sampling is uniform over the redshift range 0.2 to 0.6, in agreement with the intended BOSS target selection. The galaxy masses probed by BOSS extend over ~10^{12} M, providing unprecedented measurements of the high-mass end of the galaxy mass function. We find that the galaxy number density above ~ 2.5 10^{11} M agrees with previous determinations. We perform a comparison with semi-analytic galaxy formation models tailored to the BOSS target selection and volume, in order to contain incompleteness. The abundance of massive galaxies in the models compare fairly well with the BOSS data, but the models lack galaxies at the massive end. Moreover, no evolution with redshift is detected from ~0.6 to 0.4 in the data, whereas the abundance of massive galaxies in the models increases to redshift zero. Additionally, BOSS data display colour-magnitude (mass) relations similar to those found in the local Universe, where the most massive galaxies are the reddest. On the other hand, the model colours do not display a dependence on stellar mass, span a narrower range and are typically bluer than the observations. We argue that the lack of a colour-mass relation for massive galaxies in the models is mostly due to metallicity, which is too low in the models., Comment: 31 pages, 32 figures, accepted for publication to the Monthly Notices of the Royal Astronomical Society

We re-examine the observational evidence for large-scale (4 Mpc) galactic conformity in the local Universe, as presented in Kauffmann et al. We show that a number of methodological features of their analysis act to produce a misleadingly high amplitude of the conformity signal. These include a weighting in favour of central galaxies in very high density regions, the likely misclassification of satellite galaxies as centrals in the same high-density regions and the use of medians to characterize bimodal distributions. We show that the large-scale conformity signal in Kauffmann et al. clearly originates from a very small number of central galaxies in the vicinity of just a few very massive clusters, whose effect is strongly amplified by the methodological issues that we have identified. Some of these ‘centrals’ are likely misclassified satellites, but some may be genuine centrals showing a real conformity effect. Regardless, this analysis suggests that conformity on 4 Mpc scales is best viewed as a relatively short-range effect (at the virial radius) associated with these very large neighbouring haloes, rather than a very long-range effect (at tens of virial radii) associated with the relatively low-mass haloes that host the nominal central galaxies in the analysis. A mock catalogue constructed from a recent semi-analytic model shows very similar conformity effects to the data when analysed in the same way, suggesting that there is no need to introduce new physical processes to explain galactic conformity on 4 Mpc scales., Monthly Notices of the Royal Astronomical Society, 471 (1), ISSN:0035-8711, ISSN:1365-2966, ISSN:1365-8711

Roughly half of all stars reside in galaxies without significant ongoing star formation. However, galaxy formation models indicate that it is energetically challenging to suppress the cooling of gas and the formation of stars in galaxies that lie at the centers of their dark matter halos. In this Letter, we show that the dependence of quiescence on black hole and stellar mass is a powerful discriminant between differing models for the mechanisms that suppress star formation. Using observations of 91 star-forming and quiescent central galaxies with directly measured black hole masses, we find that quiescent galaxies host more massive black holes than star-forming galaxies with similar stellar masses. This observational result is in qualitative agreement with models that assume that effective, more-or-less continuous active galactic nucleus feedback suppresses star formation, strongly suggesting the importance of the black hole in producing quiescence in central galaxies.

We quantify the impact that a variety of galactic and environmental properties have on the quenching of star formation. We collate a sample of $\sim$ 400,000 central and $\sim$ 100,000 satellite galaxies from the Sloan Digital Sky Survey Data Release 7 (SDSS DR7). Specifically, we consider central velocity dispersion ($\sigma_{c}$), stellar, halo, bulge and disk mass, local density, bulge-to-total ratio, group-centric distance and galaxy-halo mass ratio. We develop and apply a new statistical technique to quantify the impact on the quenched fraction ($f_{\rm Quench}$) of varying one parameter, while keeping the remaining parameters fixed. For centrals, we find that the $f_{\rm Quench} - \sigma_{c}$ relationship is tighter and steeper than for any other variable considered. We compare to the Illustris hydrodynamical simulation and the Munich semi-analytic model (L-Galaxies), finding that our results for centrals are qualitatively consistent with their predictions for quenching via radio-mode AGN feedback, hinting at the viability of this process in explaining our observational trends. However, we also find evidence that quenching in L-Galaxies is too efficient and quenching in Illustris is not efficient enough, compared to observations. For satellites, we find strong evidence that environment affects their quenched fraction at fixed central velocity dispersion, particularly at lower masses. At higher masses, satellites behave identically to centrals in their quenching. Of the environmental parameters considered, local density affects the quenched fraction of satellites the most at fixed central velocity dispersion., Comment: Accepted for publication in MNRAS. 28 pages, 18 figures, 2 tables. This is the final accepted version of an earlier publication: arXiv:1412.3862

We address the fundamental question of matching the rest-frame K-band luminosity function (LF) of galaxies over the Hubble time using semi-analytic models, after modification of the stellar population modelling. We include the Maraston evolutionary synthesis models, that feature a higher contribution by the Thermally Pulsating - Asymptotic Giant Branch (TP-AGB) stellar phase, into three different semi-analytic models, namely the De Lucia and Blaizot version of the Munich model, MORGANA and the Menci model. We leave all other input physics and parameters unchanged. We find that the modification of the stellar population emission can solve the mismatch between models and the observed rest-frame K-band luminosity from the brightest galaxies derived from UKIDSS data at high redshift. For all explored semi-analytic models this holds at the redshifts - between 2 and 3 - where the discrepancy was recently pointed out. The reason for the success is that at these cosmic epochs the model galaxies have the right age (~1 Gyr) to contain a well-developed TP-AGB phase which makes them redder without the need of changing their mass or age. At the same time, the known overestimation of the faint end is enhanced in the K-band when including the TP-AGB contribution. At lower redshifts (z

We introduce a new physical recipe into the De Lucia and Blaizot version of the Munich semi-analytic model built upon the Millennium dark matter simulation: the tidal stripping of stellar material from satellite galaxies during mergers. To test the significance of the new physical process we apply a Monte Carlo Markov Chain parameter estimation technique constraining the model with the $K$-band luminosity function, $B-V$ colours and the black hole-bulge mass relation. The differences in parameter correlations, and in the allowed regions in likelihood space, reveal the impact of the new physics on the basic ingredients of the model, such as the star-formation laws, feedback recipes and the black hole growth model. With satellite disruption in place, we get a model likelihood four times higher than in the original model, indicating that the new process seems to be favoured by observations. This is achieved mainly due to a reduction in black hole growth that produces a better agreement between the properties of central black holes and host galaxies. Compared to the best-fit model without disruption, the new model removes the excess of dwarf galaxies in the original recipe with a more modest supernova heating. The new model is now consistent with the three observational data sets used to constrain it, while significantly improving the agreement with observations for the distribution of metals in stars. Moreover, the model now follows the build up of intra-cluster light.

We present an analysis of the iron abundance in the hot gas surrounding galaxy groups and clusters. To do this, we first compile and homogenise a large dataset of 79 low-redshift (|z| = 0.03) systems (159 individual measurements) from the literature. Our analysis accounts for differences in aperture size, solar abundance, and cosmology, and scales all measurements using customised radial profiles for the temperature (T), gas density, and iron abundance (Z). We then compare this dataset to groups and clusters in the L-Galaxies galaxy evolution model. Our homogenised dataset reveals a tight T-Z relation for clusters, with a scatter in Z of only 0.10 dex and a slight negative gradient. After examining potential measurement biases, we conclude that at least some of this negative gradient has a physical origin. Our model suggests greater accretion of hydrogen in the hottest systems, via stripping of gas from infalling satellites, as a cause. At lower temperatures, L-Galaxies over-estimates Z in groups, indicating that metal-rich gas removal (via e.g. AGN feedback) is required. L-Galaxies provides a reasonable match to the observed Z in the intracluster medium (ICM) of the hottest clusters from at least z ~ 1.3 to 0.3. This is achieved without needing to modify any of the galactic chemical evolution (GCE) model parameters. However, the Z in intermediate-temperature clusters appears to be under-estimated in our model at z = 0. The merits and problems with modifying the GCE modelling to correct this are discussed., Comment: 19 pages, 14 figures (excluding appendices), submitted to MNRAS

We use the semi-analytic model developed by Henriques et al. (2015) to explore the origin of star formation history diversity for galaxies that lie at the centre of their dark matter haloes and have present-day stellar masses in the range 5-8 $\times$ 10$^{10}$ M$_{\odot}$, similar to that of the Milky Way. In this model, quenching is the dominant physical mechanism for introducing scatter in the growth histories of these galaxies. We find that present-day quiescent galaxies have a larger variety of growth histories than star-formers since they underwent 'staggered quenching' - a term describing the correlation between the time of quenching and present-day halo mass. While halo mass correlates broadly with quiescence, we find that quiescence is primarily a function of black hole mass, where galaxies quench when heating from their active galactic nuclei becomes sufficient to offset the redshift-dependent cooling rate. In this model, the emergence of a prominent quiescent population is the main process that flattens the stellar mass-halo mass relation at mass scales at or above that of the Milky Way., Comment: 19 pages, 9 figures, 1 table. Accepted for publication in MNRAS

We use galaxy-galaxy lensing to study the dark matter halos surrounding a sample of Locally Brightest Galaxies (LBGs) selected from the Sloan Digital Sky Survey. We measure mean halo mass as a function of the stellar mass and colour of the central galaxy. Mock catalogues constructed from semi-analytic galaxy formation simulations demonstrate that most LBGs are the central objects of their halos, greatly reducing interpretation uncertainties due to satellite contributions to the lensing signal. Over the full stellar mass range, $10.3 < \log [M_*/M_\odot] < 11.6$, we find that passive central galaxies have halos that are at least twice as massive as those of star-forming objects of the same stellar mass. The significance of this effect exceeds $3��$ for $\log [M_*/M_\odot] > 10.7$. Tests using the mock catalogues and on the data themselves clarify the effects of LBG selection and show that it cannot artificially induce a systematic dependence of halo mass on LBG colour. The bimodality in halo mass at fixed stellar mass is reproduced by the astrophysical model underlying our mock catalogue, but the sign of the effect is inconsistent with recent, nearly parameter-free age-matching models. The sign and magnitude of the effect can, however, be reproduced by halo occupation distribution models with a simple (few-parameter) prescription for type-dependence., 20 pages, 14 figures, submitted to MNRAS; v2 has minor changes in presentation in response to referee suggestions, but no changes in results

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 sub-haloes 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 two. 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.

We introduce methods which allow observed galaxy clustering to be used together with observed luminosity or stellar mass functions to constrain the physics of galaxy formation. We show how the projected two-point correlation function of galaxies in a large semi-analytic simulation can be estimated to better than ~10% using only a very small subsample of the subhalo merger trees. This allows measured correlations to be used as constraints in a Monte Carlo Markov Chain exploration of the astrophysical and cosmological parameter space. An important part of our scheme is an analytic profile which captures the simulated satellite distribution extremely well out to several halo virial radii. This is essential to reproduce the correlation properties of the full simulation at intermediate separations. As a first application, we use low-redshift clustering and abundance measurements to constrain a recent version of the Munich semi-analytic model. The preferred values of most parameters are consistent with those found previously, with significantly improved constraints and somewhat shifted "best" values for parameters that primarily affect spatial distributions. Our methods allow multi-epoch data on galaxy clustering and abundance to be used as joint constraints on galaxy formation. This may lead to significant constraints on cosmological parameters even after marginalising over galaxy formation physics., 17 pages, 11 figures. Replaced to match the version accepted by MNRAS

We use weak gravitational lensing to measure mean mass profiles around Locally Brightest Galaxies (LBGs). These are selected from the SDSS/DR7 spectroscopic and photometric catalogues to be brighter than any neighbour projected within 1.0 Mpc and differing in redshift by $ 83\%$) are expected to be the central galaxies of their dark matter halos. Previous stacking analyses have used this LBG sample to measure mean Sunyaev-Zeldovich flux and mean X-ray luminosity as a function of LBG stellar mass. In both cases, a simulation of the formation of the galaxy population was used to estimate effective halo mass for LBGs of given stellar mass, allowing the derivation of scaling relations between the gas properties of halos and their mass. By comparing results from a variety of simulations to our lensing data, we show that this procedure has significant model dependence reflecting: (i) the failure of any given simulation to reproduce observed galaxy abundances exactly; (ii) a dependence on the cosmology underlying the simulation; and (iii) a dependence on the details of how galaxies populate halos. We use our lensing results to recalibrate the scaling relations, eliminating most of this model dependence and explicitly accounting both for residual modelling uncertainties and for observational uncertainties in the lensing results. The resulting scaling relations link the mean gas properties of dark halos to their mass over an unprecedentedly wide range, $10^{12.5}, Comment: 22 pages, 11 figures, accepted by MNRAS

We adapt the L-Galaxies semi-analytic model to follow the star-formation histories (SFH) of galaxies -- by which we mean a record of the formation time and metallicities of the stars that are present in each galaxy at a given time. We use these to construct stellar spectra in post-processing, which offers large efficiency savings and allows user-defined spectral bands and dust models to be applied to data stored in the Millennium data repository. We contrast model SFHs from the Millennium Simulation with observed ones from the VESPA algorithm as applied to the SDSS-7 catalogue. The overall agreement is good, with both simulated and SDSS galaxies showing a steeper SFH with increased stellar mass. The SFHs of blue and red galaxies, however, show poor agreement between data and simulations, which may indicate that the termination of star formation is too abrupt in the models. The mean star-formation rate (SFR) of model galaxies is well-defined and is accurately modelled by a double power law at all redshifts: SFR proportional to $1/(x^{-1.39}+x^{1.33})$, where $x=(t_a-t)/3.0\,$Gyr, $t$ is the age of the stars and $t_a$ is the loopback time to the onset of galaxy formation; above a redshift of unity, this is well approximated by a gamma function: SFR proportional to $x^{1.5}e^{-x}$, where $x=(t_a-t)/2.0\,$Gyr. Individual galaxies, however, show a wide dispersion about this mean. When split by mass, the SFR peaks earlier for high-mass galaxies than for lower-mass ones, and we interpret this downsizing as a mass-dependence in the evolution of the quenched fraction: the SFHs of star-forming galaxies show only a weak mass dependence., Accepted version of the paper, to appear in MNRAS. Compared to the original version, contains more detail on the post-processing of magnitudes, including a table of rms magnitude errors. SFHs available on Millennium database http://gavo.mpa-garching.mpg.de/MyMillennium/

We present high-redshift predictions of the star-formation-rate distribution function (SFR DF), UV luminosity function (UV LF), galactic stellar mass function (GSMF), and specific star-formation rates (sSFRs) of galaxies from the latest version of the Munich semi-analytic model L-Galaxies. We find a good fit to both the shape and normalisation of the SFR DF at $z=4-7$, apart from a slight under-prediction at the low SFR end at $z=4$. Likewise, we find a good fit to the faint number counts for the observed UV LF; at brighter magnitudes our predictions lie below the observations, increasingly so at higher redshifts. At all redshifts and magnitudes, the raw (unattenuated) number counts for the UV LF lie above the observations. Because of the good agreement with the SFR we interpret our under-prediction as an over-estimate of the amount of dust in the model for the brightest galaxies, especially at high-redshift. While the shape of our GSMF matches that of the observations, we lie between (conflicting) observations at $z=4-5$, and under-predict at $z=6-7$. The sSFRs of our model galaxies show the observed trend of increasing normalisation with redshift, but do not reproduce the observed mass dependence. Overall, we conclude that the latest version of L-Galaxies, which is tuned to match observations at $z\leq3$, does a fair job of reproducing the observed properties of galaxies at $z\geq4$. More work needs to be done on understanding observational bias at high-redshift, and upon the dust model, before strong conclusions can be drawn on how to interpret remaining discrepancies between the model and observations., Comment: Updated to reflect the published version in MNRAS (12 pages, 7 figures). Note that this was previously paper IV, but was renamed to III by the journal. Some references have been fixed, and the Bowler2014b observational errors in Figure 3 have been corrected to their proper values

We have updated the Munich galaxy formation model to the Planck first-year cosmology, while modifying the treatment of baryonic processes to reproduce recent data on the abundance and passive fractions of galaxies from z= 3 down to z=0. Matching these more extensive and more precise observational results requires us to delay the reincorporation of wind ejecta, to lower the surface density threshold for turning cold gas into stars, to eliminate ram-pressure stripping in haloes less massive than ~10^14 Msun, and to modify our model for radio mode feedback. These changes cure the most obvious failings of our previous models, namely the overly early formation of low-mass galaxies and the overly large fraction of them that are passive at late times. The new model is calibrated to reproduce the observed evolution both of the stellar mass function and of the distribution of star formation rate at each stellar mass. Massive galaxies (M>10^11 [Msun]) assemble most of their mass before z=1 and are predominantly old and passive at z=0, while lower mass galaxies assemble later and, for M, Revised to match published MNRAS version. New databases containing the full galaxy data at all redshifts and for both the Millennium and Millennium-II simulations are publicly released at http://www.mpa-garching.mpg.de/millennium/. Submission includes 15 pages of Supplementary Material with a full description of the treatment of astrophysical processes

Mergers of gas-poor galaxies, so-called dry mergers, may play a fundamental role in the assembly of the most massive galaxies, and therefore, in galaxy formation theories. Using the SDSS, we have serendipitously discovered a rare system in the observational and theoretical context, possibly a quintuple dry merger at low redshift. As a follow-up, we have obtained NOT long-slit spectra of the group, in order to measure the individual redshifts and gain insight into its merger fate. Our results show an isolated, low-redshift galaxy group consisting of massive, quiescent, early-type galaxies, composed of two clumps (possibly themselves in the process of merging), which we estimate will hypothetically merge in roughly less than a Gyr. With the possible exception of the high line-of-sight velocity dispersion, the overall properties of the system may be comparable to a compact Shakhbazyan group. However, when the small projected separations and relative mass ratios of the galaxies are taken into account in cosmological simulations, we find that this system is rather unique. We hypothesize that this group is a dry merger, whose fate will result in the assembly of an isolated, massive elliptical galaxy at low redshift., 12 pages

We present a catalogue of bulgeless galaxies, which includes 19225 objects selected in four of the deepest, largest multi-wavelength datasets available -- COSMOS, AEGIS, GEMS and GOODS -- at intermediate redshift ($0.4 \leq z \leq 1.0$). The morphological classification was provided by the Advanced Camera for Surveys General Catalogue (ACS-GC), which used publicly available data obtained with the ACS instrument on the Hubble Space Telescope. Rest-frame photometric quantities were derived using kcorrect. We analyse the properties of the sample and the evolution of pure-disc systems with redshift. Very massive [$\log (M_\star/M_{\odot}) > 10.5$] bulgeless galaxies contribute to ~30% of the total galaxy population number density at $z \geq 0.7$, but their number density drops substantially with decreasing redshift. We show that only a negligible fraction of pure discs appear to be quiescent systems, and red sequence bulgeless galaxies show indications of dust-obscured star formation. X-ray catalogues were used to search for X-ray emission within our sample. After visual inspection and detailed parametric morphological fitting we identify 30 AGN that reside in galaxies without a classical bulge. The finding of such peculiar objects at intermediate redshift shows that while AGN growth in merger-free systems is a rare event (0.2% AGN hosts in this sample of bulgeless galaxies), it can indeed happen relatively early in the Universe history.

We update the treatment of chemical evolution in the Munich semi-analytic model, L-GALAXIES. Our new implementation includes delayed enrichment from stellar winds, supernovae type II (SNe-II) and supernovae type Ia (SNe-Ia), as well as metallicity-dependent yields and a reformulation of the associated supernova feedback. Two different sets of SN-II yields and three different SN-Ia delay-time distributions (DTDs) are considered, and eleven heavy elements (including O, Mg and Fe) are self-consistently tracked. We compare the results of this new implementation with data on a) local, star-forming galaxies, b) Milky Way disc G dwarfs, and c) local, elliptical galaxies. We find that the z=0 gas-phase mass-metallicity relation is very well reproduced for all forms of DTD considered, as is the [Fe/H] distribution in the Milky Way disc. The [O/Fe] distribution in the Milky Way disc is best reproduced when using a DTD with less than or equal to 50 per cent of SNe-Ia exploding within ~400 Myrs. Positive slopes in the mass-[alpha/Fe] relations of local ellipticals are also obtained when using a DTD with such a minor `prompt' component. Alternatively, metal-rich winds that drive light alpha elements directly out into the circumgalactic medium also produce positive slopes for all forms of DTD and SN-II yields considered. Overall, we find that the best model for matching the wide range of observational data considered here should include a power-law SN-Ia DTD, SN-II yields that take account of prior mass loss through stellar winds, and some direct ejection of light alpha elements out of galaxies., 24 pages, 22 figures. Submitted to MNRAS

We apply Monte Carlo Markov Chain (MCMC) methods to large-scale simulations of galaxy formation in a LambdaCDM cosmology in order to explore how star formation and feedback are constrained by the observed luminosity and stellar mass functions of galaxies. We build models jointly on the Millennium and Millennium-II simulations, applying fast sampling techniques which allow observed galaxy abundances over the ranges 7, Comment: 25 pages, 14 figures, submitted to MNRAS

We have updated our radially-resolved SAMs of galaxy formation, which track both the atomic and molecular gas phases of the ISM. The models are adapted from those of Guo et al. using similar methodology as in Fu et al. and are run on halo merger trees from the MS and MS II with the following main changes: (1) We adopt a simple star formation law where \Sigma_SFR \propto \Sigma_H2 (2) We inject the heavy elements produced by supernovae directly into the halo hot gas, instead of first mixing them with the cold gas in the disk. (3) We include radial gas inflows in disks using a model of the form v_inflow = \alpha r. The models are used to study the radial profiles of star formation rate and gas-phase metallicity in present-day galaxies. The \Sigma_H2 profiles in L* galaxies place strong constraints on inflow velocities, favouring models where v_inflow~7km/s at a galactocentric radius of 10kpc. Radial gas inflow has little influence on gas-phase and stellar metallicity gradients, which are affected much more strongly by the fraction of metals that are directly injected into the halo gas, rather than mixed with the cold gas. Metals ejected out of the galaxy in early epochs result in late infall of pre-enriched gas and flatter present-day gas-phase metallicity gradients. A prescription in which 80% of the metals are injected into the halo gas results in good fits to the flat observed metallicity gradients in galaxies with stellar masses greater than 10^10 M_sun, as well as the relations between gas-phase metallicity and sSFR in the outer parts of galactic disks. We examine the correlation between gas-phase metallicity gradient and global galaxy properties, finding that it is most strongly correlated with the bulge-to-total ratio of the galaxy. This is because gas is consumed when the bulge forms during galaxy mergers, and the gas-phase metallicity gradient is then set by newly-accreted gas., Comment: 20 pages, 13 figures, accepted for publication in MNRAS. Newer version rephrases the descriptions on the physical recipes and model results according to the refree's report

Simulations of galaxy evolution aim to capture our current understanding as well as to make predictions for testing by future experiments. Simulations and observations are often compared in an indirect fashion: physical quantities are estimated from the data and compared to models. However, many applications can benefit from a more direct approach, where the observing process is also simulated and the models are seen fully from the observer's perspective. To facilitate this, we have developed the Millennium Run Observatory (MRObs), a theoretical virtual observatory which uses virtual telescopes to `observe' semi-analytic galaxy formation models based on the suite of Millennium Run dark matter simulations. The MRObs produces data that can be processed and analyzed using the standard software packages developed for real observations. At present, we produce images in forty filters from the rest-frame UV to IR for two stellar population synthesis models, three different models of IGM absorption, and two cosmologies (WMAP1/7). Galaxy distributions for a large number of mock lightcones can be `observed' using models of major ground- and space-based telescopes. The data include lightcone catalogues linked to structural properties of galaxies, pre-observation model images, mock telescope images, and Source Extractor products that can all be traced back to the higher level dark matter, semi-analytic galaxy, and lightcone catalogues available in the Millennium database. Here, we describe our methods and announce a first public release of simulated surveys (e.g., SDSS, CFHT-LS, GOODS, GOODS/ERS, CANDELS, and HUDF). The MRObs browser, an online tool, further facilitates exploration of the simulated data. We demonstrate the benefits of a direct approach through a number of example applications (galaxy number counts in CANDELS, clusters, morphologies, and dropout selections)., MNRAS, in press. Millennium Run Observatory data products, online tools, and more available through http://galformod.mpa-garching.mpg.de/mrobs/

We present the Spitzer Extragalactic Representative Volume Survey (SERVS), an 18 square degrees medium-deep survey at 3.6 and 4.5 microns with the post-cryogenic Spitzer Space Telescope to ~2 microJy (AB=23.1) depth of five highly observed astronomical fields (ELAIS-N1, ELAIS-S1, Lockman Hole, Chandra Deep Field South and XMM-LSS). SERVS is designed to enable the study of galaxy evolution as a function of environment from z~5 to the present day, and is the first extragalactic survey both large enough and deep enough to put rare objects such as luminous quasars and galaxy clusters at z>1 into their cosmological context. SERVS is designed to overlap with several key surveys at optical, near- through far-infrared, submillimeter and radio wavelengths to provide an unprecedented view of the formation and evolution of massive galaxies. In this paper, we discuss the SERVS survey design, the data processing flow from image reduction and mosaicing to catalogs, as well as coverage of ancillary data from other surveys in the SERVS fields. We also highlight a variety of early science results from the survey., 21 pages, 23 figures, published in PASP. This version differs slightly from PASP, mainly due to formatting issues. Figure 2 was also corrected from the previous version

Using the technique of Angulo & White (2010) we scale the Millennium and Millennium-II simulations of structure growth in a LCDM universe from the cosmological parameters with which they were carried out (based on first-year results from the Wilkinson Microwave Anisotropy Probe, WMAP1) to parameters consistent with the seven-year WMAP data (WMAP7). We implement semi-analytic galaxy formation modelling on both simulations in both cosmologies to investigate how the formation, evolution and clustering of galaxies are predicted to vary with cosmological parameters. The increased matter density Omega_m and decreased linear fluctuation amplitude sigma8 in WMAP7 have compensating effects, so that the abundance and clustering of dark halos are predicted to be very similar to those in WMAP1 for z 1 galaxies are predicted to be more strongly clustered for WMAP7. Differences in galaxy properties, including, clustering, in these two cosmologies are rather small up to redshift 3. Given that there are still considerable residual uncertainties in galaxy formation models, it is very difficult to distinguish WMAP1 from WMAP7 through observations of galaxy properties or their evolution., 16 pages, accepted version (MNRAS)

We use the Spitzer Wide-area InfraRed Extragalactic Legacy Survey (SWIRE) to explore the specific star-formation activity of galaxies and their evolution near the peak of the cosmic far-infrared (FIR) background at 70 and 160um. We use a stacking analysis to determine the mean FIR properties of well defined subsets of galaxies at flux levels well below the FIR catalogue detection limits of SWIRE and other Spitzer surveys. We tabulate the contribution of different subsets of galaxies to the FIR background at 70um and 160um. These long wavelengths provide a good constraint on the bolometric, obscured emission. The large area provides good constraints at low z and in finer redshift bins than previous work. At all redshifts we find that the specific FIR Luminosity (sLFIR) decreases with increasing mass, following a trend L_FIR/M* propto M_* ^beta with beta =-0.38\pm0.14. This is a more continuous change than expected from the {Delucia2007} semi-analytic model suggesting modifications to the feedback prescriptions. We see an increase in the sLFIR by about a factor of ~100 from 0, Comment: v2 Update doesn't change the content of the paper, but now includes data files for the plots Fig 5-13 (all.plotdat, spi.plotdat and ell.plotdat on arXiv package)

[abridged] We present a statistical exploration of the parameter space of the De Lucia and Blaizot version of the Munich semi-analytic model built upon the millennium dark matter simulation. This is achieved by applying a Monte Carlo Markov Chain method to constrain the 6 free parameters that define the stellar and black-hole mass functions at redshift zero. The model is tested against three different observational data sets, including the galaxy K-band luminosity function, B-V colours, and the black hole-bulge mass relation, separately and combined, to obtain mean values, confidence limits and likelihood contours for the best fit model. Using each observational data set independently, we discuss how the SA model parameters affect each galaxy property and to what extent the correlations between them can lead to improved understandings of the physics of galaxy formation. When all the observations are combined, we find reasonable agreement between the majority of the previously published parameter values and our confidence limits. However, the need to suppress dwarf galaxy formation requires the strength of the supernova feedback to be significantly higher in our best-fit solution than in previous work. To balance this, we require the feedback to become ineffective in halos of lower circular velocity than before, so as to permit the formation of sufficient high-luminosity galaxies: unfortunately, this leads to an excess of galaxies around L*. Although the best-fit is formally consistent with the data, there is no region of parameter space that reproduces the shape of galaxy luminosity function across the whole magnitude-range. We discuss modifications to the semi-analytic model that might simultaneously improve the fit to the observed luminosity function and reduce the reliance on excessive supernova feedback in small halos., accepted by MNRAS