We present an analysis of the origin and properties of the circum-galactic medium (CGM) in a suite of 11 cosmological zoom simulations resembling present day spiral galaxies. On average the galaxies retain about 50\% of the cosmic fraction in baryons, almost equally divided into disc (interstellar medium) gas, cool CGM gas and warm-hot CGM gas. At radii smaller than 50 kpc the CGM is dominated by recycled warm-hot gas injected from the central galaxy, while at larger radii it is dominated by cool gas accreted onto the halo. The recycled gas typically accounts for one-third of the CGM mass. We introduce the novel publicly available analysis tool \textsc{pygad} to compute ion abundances and mock absorption spectra. For Lyman-${\alpha}$ absorption we find good agreement of the simulated equivalent width (EW) distribution and observations out to large radii. Disc galaxies with quiescent assembly histories show significantly more absorption along the disc major axis. By comparing the EW and HI column densities we find that CGM Lyman-${\alpha}$ absorbers are best represented by an effective line-width $b\approx 50 - 70$ km s$^{-1}$ that increases mildly with halo mass, larger than typically assumed.

We perform a comparison, object-by-object and statistically, between the Munich semi-analytical model, L-Galaxies, and the IllustrisTNG hydrodynamical simulations. By running L-Galaxies on the IllustrisTNG dark matter-only merger trees, we identify the same galaxies in the two models. This allows us to compare the stellar mass, star formation rate and gas content of galaxies, as well as the baryonic content of subhaloes and haloes in the two models. We find that both the stellar mass functions and the stellar masses of individual galaxies agree to better than $\sim0.2\,$dex. On the other hand, specific star formation rates and gas contents can differ more substantially. At $z=0$ the transition between low-mass star-forming galaxies and high-mass, quenched galaxies occurs at a stellar mass scale $\sim0.5\,$dex lower in IllustrisTNG than in L-Galaxies. IllustrisTNG also produces substantially more quenched galaxies at higher redshifts. Both models predict a halo baryon fraction close to the cosmic value for clusters, but IllustrisTNG predicts lower baryon fractions in group environments. These differences are due primarily to differences in modelling feedback from stars and supermassive black holes. The gas content and star formation rates of galaxies in and around clusters and groups differ substantially, with IllustrisTNG satellites less star-forming and less gas-rich. We show that environmental processes such as ram-pressure stripping are stronger and operate to larger distances and for a broader host mass range in IllustrisTNG. We suggest that the treatment of galaxy evolution in the semi-analytic model needs to be improved by prescriptions which capture local environmental effects more accurately.

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<=z<=3 and mass range 8<=log(M_*/Msun)<=12, the black hole-bulge mass relation at z=0, galaxy morphology as a function of stellar mass and the mass-metallicity relations of both stellar and gas components. In addition, its detailed modelling of the radial structure of discs allows qualitatively new comparisons with observation, most notably with the relative sizes and masses of the stellar, atomic and molecular components in discs. Good agreement is found with recent data. Comparison of results obtained for simulations differing in mass resolution by more than two orders of magnitude shows that all important distributions are numerically well converged even for this more detailed model. An examination of metallicity and surface-density gradients in the stars and gas indicates that our new model, with star formation, chemical enrichment and feedback calculated self-consistently on local disc scales, reproduces some but not all of the trends seen in recent many-galaxy IFU surveys.

Using a set of 15 high-resolution magnetohydrodynamic cosmological simulations of Milky Way formation, we investigate the origin of the baryonic material found in stars at redshift zero. We find that roughly half of this material originates from subhalo/satellite systems and half is smoothly accreted from the Inter-Galactic Medium (IGM). About $90 \%$ of all material has been ejected and re-accreted in galactic winds at least once. The vast majority of smoothly accreted gas enters into a galactic fountain that extends to a median galactocentric distance of $\sim 20$ kpc with a median recycling timescale of $\sim 500$ Myr. We demonstrate that, in most cases, galactic fountains acquire angular momentum via mixing of low-angular momentum, wind-recycled gas with high-angular momentum gas in the Circum-Galactic Medium (CGM). Prograde mergers boost this activity by helping to align the disc and CGM rotation axes, whereas retrograde mergers cause the fountain to lose angular momentum. Fountain flows that promote angular momentum growth are conducive to smooth evolution on tracks quasi-parallel to the disc sequence of the stellar mass-specific angular momentum plane, whereas retrograde minor mergers, major mergers and bar-driven secular evolution move galaxies towards the bulge-sequence. Finally, we demonstrate that fountain flows act to flatten and narrow the radial metallicity gradient and metallicity dispersion of disc stars, respectively. Thus, the evolution of galactic fountains depends strongly on the cosmological merger history and is crucial for the chemo-dynamical evolution of Milky Way-sized disc galaxies.

We have analyzed the emission line properties of 6019 Type II AGN at redshifts between 0.4-0.8 with [OIII] luminosities greater than 3x10^8 Lsun, characteristic of the Type II quasars first identified in population studies by Zakamska et al (2003). The AGN are drawn from the CMASS sample of galaxies with stellar masses greater than 10^11 Msun that were studied as part of the Baryon Oscillation Spectroscopic Survey (BOSS) and comprise 0.5% of the total population of these galaxies. Individual spectra have low S/N, so the analysis is carried out on stacked spectra in bins of [OIII] luminosity and estimated stellar age. The emission line ratios of the stacks are well fit with simple uniform-density photoionization models with metallicities between solar and twice solar. In the stacks, a number of emission lines are found to have distinct broad components requiring a double Gaussian rather than a single Gaussian fit, indicative of outflowing ionized gas. These are: [OIII]4959, [OIII]5007, [OII]3727,3729 and Halpha$. Higher ionization lines such as [NeIII]3869 and [NeV]3345 are detected in the stacks, but are well fit by single Gaussians. The broad components typically contain a third of the total line flux and have widths of 600 km/s for the oxygen lines and 900 km/s for Halpha. The fraction of the flux in the broad component and its width are independent of [OIII] luminosity, stellar age, radio and mid-IR luminosity. The stellar mass of the galaxy is the only parameter we could identify that influences the width of the broad line component.

We have carried out a controlled comparison of the structural and kinematic properties of the circumgalactic medium (CGM) around Milky Way mass galaxies in the Illustris and IllustrisTNG simulations. Very striking differences are found. At z=0, gas column density and temperature profiles at large radii ($\sim 100$ kpc) correlate strongly with disk gas mass fraction in Illustris, but not in TNG. The neutral gas at large radii is preferentially aligned in the plane of the disk in TNG, whereas it is much more isotropic in Illustris. The vertical coherence scale of the rotationally supported gas in the CGM is linked to the gas mass fraction of the galaxy in Illustris, but not in TNG. A tracer particle analysis allows us to show how these differences can be understood as a consequence of the different sub-grid models of feedback in the two simulations. A study of spatially matched galaxies in the two simulations shows that in TNG, feedback by supernovae and AGN helps to create an extended smooth reservoir of hot gas at high redshifts, that then cools to form a thin, rotationally-supported disk at later times. In Illustris, AGN dump heat in the form of hot gas bubbles that push diffuse material at large radii out of the halo. The disk is formed by accretion of colder, recycled material, and this results in more vertically extended gas distributions above and below the Galactic plane. We conclude that variations in the structure of gas around Milky Way mass galaxies are a sensitive probe of feedback physics in simulations and are worthy of more observational consideration in future.

We introduce a Local Background Environment (LBE) estimator that can be measured in and around every galaxy or its dark matter subhalo in high-resolution cosmological simulations. The LBE is designed to capture the influence of environmental effects such as ram-pressure stripping on the formation and evolution of galaxies in semi-analytical models. We define the LBE directly from the particle data within an adaptive spherical shell, and devise a Gaussian mixture estimator to separate background particles from previously unidentified subhalo particles. Analyzing the LBE properties, we find that the LBE of satellite galaxies is not at rest with respect to their host halo, in contrast to typical assumptions. The orientations of the velocities of a subhalo and its LBE are well aligned in the outer infall regions of haloes, but decorrelated near halo center. Significantly, there is no abrupt change in LBE velocity or density at the halo virial radius. This suggests that stripping should also happen beyond this radius. Therefore, we use the time-evolving LBE of galaxies to develop a method to better account for ram-pressure stripping within the Munich semi-analytical model, L-Galaxies. Overall, our new approach results in a significant increase in gas stripping across cosmic time. Central galaxies, as well as satellites beyond the virial radius, can lose a significant fraction of their hot halo gas. As a result, the gas fractions and star formation rates of satellite galaxies are suppressed relative to the fiducial model, although the stellar masses and global stellar mass functions are largely unchanged.

We study an idealized 1D model for the evolution of hot gas in dark matter halos for redshifts $z=[0,6]$. We introduce a numerical setup incorporating cosmological accretion of gas, along with the growth of the halo, based on the Van den Bosch model for the average growth of halos as a function of cosmic time. We evolve one-dimensional Lagrangian shells with radiative cooling of the gas and heating due to feedback from the gas cooling and moving in toward the center. A simple Bondi accretion model on to a central black hole is used to include feedback heating. The setup captures some of the key characteristics of spherically symmetric accretion onto the halos: formation of virial shocks slightly outside $r_{200}$ and long-term thermal balance in the form of cooling and heating cycles. The gas density outside our initial halos at $z=6$ is constrained by requiring that the baryon fraction within the virial radius for non-radiative evolution be equal to the universal value at almost all times. The total mass in the cold phase (taken to be $\sim 10^4$ K) within $40$ kpc is tracked as a function of the halo mass and redshift. We compare the evolution of the cold gas mass to the observed stellar-mass versus halo mass relations, following which, we can constrain the feedback energy required for different halo masses and redshifts. We also compare and match the hot gas density and temperature profiles for our most massive halo to those of clusters observed upto redshift $2$. Our model is thus an improvement over the semi-analytic models in which isothermal condition and $\rho \propto r^{-2}$ are assumed.

We examine the stellar haloes of the Auriga simulations, a suite of thirty cosmological magneto-hydrodynamical high-resolution simulations of Milky Way-mass galaxies performed with the moving-mesh code AREPO. We study halo global properties and radial profiles out to $\sim 150$ kpc for each individual galaxy. The Auriga haloes are diverse in their masses and density profiles; mean metallicity and metallicity gradients; ages; and shapes, reflecting the stochasticity inherent in their accretion and merger histories. A comparison with observations of nearby late-type galaxies shows very good agreement between most observed and simulated halo properties. However, Auriga haloes are typically too massive. We find a connection between population gradients and mass assembly history: galaxies with few significant progenitors have more massive haloes, possess large negative halo metallicity gradients and steeper density profiles. The number of accreted galaxies, either disrupted or under disruption, that contribute 90% of the accreted halo mass ranges from 1 to 14, with a median of 6.5, and their stellar masses span over three orders of magnitude. The observed halo mass--metallicity relation is well reproduced by Auriga and is set by the stellar mass and metallicity of the dominant satellite contributors. This relationship is found not only for the accreted component but also for the total (accreted + in-situ) stellar halo. Our results highlight the potential of observable halo properties to infer the assembly history of galaxies.

We study the narrow emission line properties and stellar populations of a sample of 1385 AGN selected to have strong excess emission at mid-infrared wavelengths based on comparing Wide-field Infrared Survey Explorer W1-W2 band colours with optical stellar absorption line indicators. Our goal is to understand whether the physical conditions in the interstellar medium of these objects differ from those of local AGN selected by their optical emission line ratios. To enable this comparison, we construct a control sample of 50,000 optically-selected AGN with the same redshifts that do not have strong mid-IR excess emission. The mid-IR excess and control samples differ strongly in [OIII] line luminosity, ionized gas excitation mechanism, ionization state and electron density. We show that the radio-detected, mid-IR excess AGN constitute the most luminous and highly ionized AGN in the local Universe and they contribute primarily to the growth of black holes in the most massive galaxies. At least half of this black hole growth is occurring in galaxies with recent starbursts. The morphologies of these systems indicate that the starbursts have probably been triggered by galaxy-galaxy mergers and interactions. The most luminous AGN in our mid-IR excess sample have properties that are similar to the Type II quasars studied at higher redshifts. In contrast, the control sample constitute a class of lower ionization, less luminous AGN in more quiescent galaxies that contribute primarily to the growth of low mass black holes.

We explore the abundance, spatial distribution, and physical properties of the OVI, OVII, and OVIII ions of oxygen in circumgalactic and intergalactic media (the CGM, IGM, and WHIM). We use the TNG100 and TNG300 large volume cosmological magneto-hydrodynamical simulations. Modeling the ionization states of simulated oxygen, we find good agreement with observations of the low-redshift OVI column density distribution function (CDDF), and present its evolution for all three ions from z=0 to z=4. Producing mock quasar absorption line spectral surveys, we show that the IllustrisTNG simulations are fully consistent with constraints on the OVI content of the CGM from COS-Halos and other low redshift observations, producing columns as high as observed. We measure the total amount of mass and average column densities of each ion using hundreds of thousands of simulated galaxies spanning 10^11 < Mhalo/Msun < 10^15 corresponding to 10^9 < M*/Msun < 10^12 in stellar mass. Stacked radial profiles of OVI are computed in 3D number density and 2D projected column, decomposing into the 1-halo and 2-halo terms. Relating halo OVI to properties of the central galaxy, we find a correlation between the (g-r) color of a galaxy and the total amount of OVI in its CGM. In comparison to the COS-Halos finding, this leads to a dichotomy of columns around star-forming versus passive galaxies at fixed stellar (or halo) mass. We demonstrate that this correlation is a direct result of blackhole feedback associated with quenching and represents a causal consequence of galactic-scale baryonic feedback impacting the physical state of the circumgalactic medium.

We introduce the first two simulations of the IllustrisTNG project, a next generation of cosmological magnetohydrodynamical simulations, focusing on the optical colors of galaxies. We explore TNG100, a rerun of the original Illustris box, and TNG300, which includes 2x2500^3 resolution elements in a volume twenty times larger. Here we present first results on the galaxy color bimodality at low redshift. Accounting for the attenuation of stellar light by dust, we compare the simulated (g-r) colors of 10^9 < M*/Msun < 10^12.5 galaxies to the observed distribution from the Sloan Digital Sky Survey (SDSS). We find a striking improvement with respect to the original Illustris simulation, as well as excellent quantitative agreement in comparison to the observations, with a sharp transition in median color from blue to red at a characteristic M* ~ 10^10.5 Msun. Investigating the build-up of the color-mass plane and the formation of the red sequence, we demonstrate that the primary driver of galaxy color transition in the TNG model is supermassive blackhole feedback in its low-accretion state. Across the entire population we measure a median color transition timescale dt_green of ~1.6 Gyr, a value which drops for increasingly massive galaxies. We find signatures of the physical process of quenching: at fixed stellar mass, redder galaxies have lower SFRs, gas fractions, and gas metallicities; their stellar populations are also older and their large-scale interstellar magnetic fields weaker than in bluer galaxies. Finally, we measure the amount of stellar mass growth on the red sequence. Galaxies with M* > 10^11 Msun which redden at z<1 accumulate on average ~25% of their final z=0 mass post-reddening; at the same time, ~18% of such massive galaxies acquire half or more of their final stellar mass while on the red sequence.

The physical origin of the large-scale conformity in the colours and specific star formation rates of isolated low mass central galaxies and their neighbours on scales in excess of 1 Mpc is still under debate. One possible scenario is that gas is heated over large scales by feedback from active galactic nuclei (AGN), leading to coherent modulation of cooling and star formation between well-separated galaxies. In this Letter, the metal line absorption catalogue of Zhu & Menard (2013) is used to probe gas out to large projected radii around a sample of a million galaxies with stellar masses ~10^{10} M_{sun} and photometric redshifts in the range 0.4<z<0.8 selected from Sloan Digital Sky Survey imaging data. This galaxy sample covers an effective volume of 2.2 Gpc^3. A statistically significant excess of MgII absorbers is present around the red low mass galaxies compared to their blue counterparts out to projected radii of 10 Mpc. In addition, the equivalent width distribution function of MgII absorbers around low mass galaxies is shown to be strongly affected by the presence of a nearby (R_p<2 Mpc) radio-loud AGN out to projected radii of 5 Mpc.

The goal of this paper is to investigate the physical nature of galaxies in the redshift range $0.02<z<0.15$ that have strong excess emission at mid-IR wavelengths and to determine whether they host a population of accreting black holes that cannot be identified using optical emission lines. We show that at fixed stellar mass $M_*$ and $D_n(4000)$, the distribution of [3.4]-[4.6] $\mu$m (WISE W1-W2 band) colours is sharply peaked, with a long tail to much redder W1-W2 colours. We introduce a procedure to pull out the red outlier population based on a combination of three stellar population diagnostics. When compared with optically-selected AGN, red outliers are more likely to be found in massive galaxies, and they tend to have lower stellar mass densities, younger stellar ages and higher dust content than optically-selected AGN hosts. They are twice as likely to be detected at radio wavelengths. We examine W1-W2 colour profiles for a subset of the nearest, reddest outliers and find that most are not centrally peaked, indicating that the hot dust emission is spread throughout the galaxy. We find that radio luminosity is the quantity that is most predictive of a redder central W1-W2 colour. Radio-loud galaxies with centrally concentrated hot dust emission are almost always morphologically disturbed, with compact, unresolved emission at 1.4 Ghz. Eighty percent of such systems are identifiable as AGN using optical emission line diagnostics.

We introduce xCOLD GASS, a legacy survey providing a census of molecular gas in the local Universe. Building upon the original COLD GASS survey, we present here the full sample of 532 galaxies with CO(1-0) measurements from the IRAM-30m telescope. The sample is mass-selected in the redshift interval $0.01<z<0.05$ from SDSS, and therefore representative of the local galaxy population with M$_{\ast}>10^9$M$_{\odot}$. The CO(1-0) flux measurements are complemented by observations of the CO(2-1) line with both the IRAM-30m and APEX telescopes, HI observations from Arecibo, and photometry from SDSS, WISE and GALEX. Combining the IRAM and APEX data, we find that the CO(2-1) to CO(1-0) luminosity ratio for integrated measurements is $r_{21}=0.79\pm0.03$, with no systematic variations across the sample. The CO(1-0) luminosity function is constructed and best fit with a Schechter function with parameters {$L_{\mathrm{CO}}^* = (7.77\pm2.11) \times 10^9\,\mathrm{K\,km\,s^{-1}\, pc^{2}}$, $\phi^{*} = (9.84\pm5.41) \times 10^{-4} \, \mathrm{Mpc^{-3}}$ and $\alpha = -1.19\pm0.05$}. With the sample now complete down to stellar masses of $10^9$M$_{\odot}$, we are able to extend our study of gas scaling relations and confirm that both molecular gas fraction and depletion timescale vary with specific star formation rate (or offset from the star-formation main sequence) much more strongly than they depend on stellar mass. Comparing the xCOLD GASS results with outputs from hydrodynamic and semi-analytic models, we highlight the constraining power of cold gas scaling relations on models of galaxy formation.

We describe the sample design for the SDSS-IV MaNGA survey and present the final properties of the main samples along with important considerations for using these samples for science. Our target selection criteria were developed while simultaneously optimizing the size distribution of the MaNGA integral field units (IFUs), the IFU allocation strategy, and the target density to produce a survey defined in terms of maximizing S/N, spatial resolution, and sample size. Our selection strategy makes use of redshift limits that only depend on i-band absolute magnitude ($M_i$), or, for a small subset of our sample, $M_i$ and color (NUV-i). Such a strategy ensures that all galaxies span the same range in angular size irrespective of luminosity and are therefore covered evenly by the adopted range of IFU sizes. We define three samples: the Primary and Secondary samples are selected to have a flat number density with respect to $M_i$ and are targeted to have spectroscopic coverage to 1.5 and 2.5 effective radii (Re), respectively. The Color-Enhanced supplement increases the number of galaxies in the low-density regions of color-magnitude space by extending the redshift limits of the Primary sample in the appropriate color bins. The samples cover the stellar mass range $5\times10^8 \leq M_* \leq 3\times10^{11} M_{\odot}$ and are sampled at median physical resolutions of 1.37 kpc and 2.5 kpc for the Primary and Secondary samples respectively. We provide weights that will statistically correct for our luminosity and color-dependent selection function and IFU allocation strategy, thus correcting the observed sample to a volume limited sample.

We use the latest version of the metal line absorption catalogue of Zhu & M\'enard (2013) to study the clustering of MgII absorbers around massive galaxies (~10^11.5 M_sun), quasars and radio-loud AGN with redshifts between 0.4 and 0.75. Clustering is evaluated in two dimensions, by binning absorbers both in projected radius and in velocity separation. Excess MgII is detected around massive galaxies out to R_p=20 Mpc. At projected radii less than 800 kpc, the excess extends out to velocity separations of 10,000 km/s. The extent of the high velocity tail within this radius is independent of the mean stellar age of the galaxy and whether or not it harbours an active galactic nucleus. We interpret our results using the publicly available Illustris and Millennium simulations. Models where the MgII absorbers trace the dark matter particle or subhalo distributions do not fit the data. They overpredict the clustering on small scales and do not reproduce the excess high velocity separation MgII absorbers seen within the virial radius of the halo. The Illustris simulations which include thermal, but not mechanical feedback from AGN, also do not provide an adequate fit to the properties of the cool halo gas within the virial radius. We propose that the large velocity separation MgII absorbers trace gas that has been pushed out of the dark matter halos, possibly by multiple episodes of AGN-driven mechanical feedback acting over long timescales.

We have conducted a study of extra-planar diffuse ionized gas using the first year data from the MaNGA IFU survey. We have stacked spectra from 49 edge-on, late-type galaxies as a function of distance from the midplane of the galaxy. With this technique we can detect the bright emission lines Halpha, Hbeta, [OII]3726, 3729, [OIII]5007, [NII]6549, 6584, and [SII]6717, 6731 out to about 4 kpc above the midplane. With 16 galaxies we can extend this analysis out to about 9 kpc, i.e. a distance of ~2R_e, vertically from the midplane. In the halo, the surface brightnesses of the [OII] and Halpha emission lines are comparable, unlike in the disk where Halpha dominates. When we split the sample by specific star formation rate, concentration index, and stellar mass, each subsample's emission line surface brightness profiles and ratios differ, indicating that extra-planar gas properties can vary. The emission line surface brightnesses of the gas around high specific star formation rate galaxies are higher at all distances, and the line ratios are closer to ratios characteristic of HII regions compared with low specific star formation rate galaxies. The less concentrated and lower stellar mass samples exhibit line ratios that are more like HII regions at larger distances than their more concentrated and higher stellar mass counterparts. The largest difference between different subsamples occurs when the galaxies are split by stellar mass. We additionally infer that gas far from the midplane in more massive galaxies has the highest temperatures and steepest radial temperature gradients based on their [NII]/Halpha and [OII]/Halpha ratios between the disk and the halo.

We use the combined data from the COS-GASS and COS-Halos surveys to characterize the Circum-Galactic Medium (CGM) surrounding typical low-redshift galaxies in the mass range $\rm~M_*\sim~10^{9.5-11.5}~M_{\odot} $, and over a range of impact parameters extending to just beyond the halo virial radius ($\rm~R_{vir}$). We find the radial scale length of the distributions of the equivalent widths of the Lyman~$\alpha$ and Si III absorbers to be 0.9 and 0.4 $\rm~R_{vir}$, respectively. The radial distribution of equivalent widths is relatively uniform for the blue galaxies, but highly patchy (low covering fraction) for the red galaxies. We also find that the Lyman~$\alpha$ and Si III equivalent widths show significant positive correlations with the specific star-formation rate (sSFR) of the galaxy. We find a surprising lack of correlations between the halo mass (virial velocity) and either the velocity dispersions or velocity offsets of the Lyman~$\alpha$ lines. The ratio of the velocity offset to the velocity dispersion for the Lyman~$\alpha$ absorbers has a mean value of $\sim$ 4, suggesting that a given the line-of-sight is intersecting a dynamically coherent structure in the CGM rather than a sea of orbiting clouds. The kinematic properties of the CGM are similar in the blue and red galaxies, although we find that a significantly larger fraction of the blue galaxies have large Lyman~$\alpha$ velocity offsets (>200 km s$^{-1}$). We show that - if the CGM clouds represent future fuel for star-formation - our new results could imply a large drop in the specific star-formation rate across the galaxy mass-range we probe.

We analyze the properties of the circumgalactic gas (CGM) around 120 galaxies with stellar and dark matter halo masses similar to that of the Milky Way. We focus on the morphology and kinematics of the neutral hydrogen and how this depends on f_g, the ratio of gas-to-stellar mass within the optical radius. In gas-rich galaxies with f_g > 0.1, gas temperatures rise slowly from center of the halo out to the virial radius and average neutral gas column densities remain above 10^19 atoms cm^-2 out to radii of 50-70 kpc. In gas-poor galaxies with f_g < 0.1, gas temperatures rise quickly outside the edge of the disk to 10^6 K, and then remain fixed out to radii of 100 kpc. The column density of neutral gas quickly drops below 10^19 atoms cm^-2 at radii of 10 kpc. Neutral gas distributions are also more asymmetric in gas-poor galaxies. Most of the differences between gas-poor and gas-rich galaxies in the Illustris simulation can be attributed to the effects of "radio-mode" AGN feedback. In the Illustris simulation, the circumgalactic gas is found to rotate coherently about the center of the galaxy with a maximum rotational velocity of around 200 km/s. In gas-rich galaxies, the average coherence length of the rotating gas is 40 kpc, compared to 10 kpc in gas-poor galaxies. In the very most gas-rich systems, the CGM can rotate coherently over scales of 70-100 kpc. We discuss our results in the context of recent observations of the CGM in low mass galaxies via UV absorption-line spectroscopy and deep 21cm observations of edge-on spiral galaxies.

Mapping Nearby Galaxies at Apache Point Observatory (MaNGA) is an optical fiber-bundle integral-field unit (IFU) spectroscopic survey that is one of three core programs in the fourth-generation Sloan Digital Sky Survey (SDSS-IV). With a spectral coverage of 3622 - 10,354 Angstroms and an average footprint of ~ 500 arcsec^2 per IFU the scientific data products derived from MaNGA will permit exploration of the internal structure of a statistically large sample of 10,000 low redshift galaxies in unprecedented detail. Comprising 174 individually pluggable science and calibration IFUs with a near-constant data stream, MaNGA is expected to obtain ~ 100 million raw-frame spectra and ~ 10 million reduced galaxy spectra over the six-year lifetime of the survey. In this contribution, we describe the MaNGA Data Reduction Pipeline (DRP) algorithms and centralized metadata framework that produces sky-subtracted, spectrophotometrically calibrated spectra and rectified 3-D data cubes that combine individual dithered observations. For the 1390 galaxy data cubes released in Summer 2016 as part of SDSS-IV Data Release 13 (DR13), we demonstrate that the MaNGA data have nearly Poisson-limited sky subtraction shortward of ~ 8500 Angstroms and reach a typical 10-sigma limiting continuum surface brightness mu = 23.5 AB/arcsec^2 in a five arcsec diameter aperture in the g band. The wavelength calibration of the MaNGA data is accurate to 5 km/s rms, with a median spatial resolution of 2.54 arcsec FWHM (1.8 kpc at the median redshift of 0.037) and a median spectral resolution of sigma = 72 km/s.

A recent observational study of haloes of nearby Milky Way-like galaxies shows that only half (four out of eight) of the current sample exhibits strong negative metallicity ([Fe/H]) gradients. This is at odds with predictions from hydrodynamical simulations where such gradients are ubiquitous. In this Letter, we use high resolution cosmological hydrodynamical simulations to study the [Fe/H] distribution of galactic haloes. We find that kinematically selected stellar haloes, including both in-situ and accreted particles, have an oblate [Fe/H] distribution. Spherical [Fe/H] radial profiles show strong negative gradients within 100 kpc, in agreement with previous numerical results. However, the projected median [Fe/H] profiles along the galactic disc minor axis, typically obtained in observations, are significantly flatter. The median [Fe/H] values at a given radius are larger for the spherical profiles than for the minor axis profiles by as much as 0.4 dex within the inner 50 kpc. Similar results are obtained if only the accreted stellar component is considered indicating that the differences between spherical and minor axis profiles are not purely driven by heated disc star particles formed in situ. Our study highlights the importance of performing careful comparisons between models and observations of halo [Fe/H] distributions.

The quenching of star formation in satellite galaxies is observed over a wide range of dark matter halo masses and galaxy environments. In the recent Guo et al (2011) and Fu et al (2013) semi-analytic + N-body models, the gaseous environment of the satellite galaxy is governed by the properties of the dark matter subhalo in which it resides. This quantity depends of the resolution of the N-body simulation, leading to a divergent fraction of quenched satellites in high- and low-resolution simulations. Here, we incorporate an analytic model to trace the subhaloes below the resolution limit. We demonstrate that we then obtain better converged results between the Millennium I and II simulations, especially for the satellites in the massive haloes ($\rm log M_{halo}=[14,15]$). We also include a new physical model for the ram-pressure stripping of cold gas in satellite galaxies. However, we find very clear discrepancies with observed trends in quenched satellite galaxy fractions as a function of stellar mass at fixed halo mass. At fixed halo mass, the quenched fraction of satellites does not depend on stellar mass in the models, but increases strongly with mass in the data. In addition to the over-prediction of low-mass passive satellites, the models also predict too few quenched central galaxies with low stellar masses, so the problems in reproducing quenched fractions are not purely of environmental origin. Further improvements to the treatment of the gas-physical processes regulating the star formation histories of galaxies are clearly necessary to resolve these problems.

We compare the mass and internal distribution of atomic hydrogen (HI) in 2200 present-day central galaxies with M_star > 10^10 M_Sun from the 100 Mpc EAGLE Reference simulation to observational data. Atomic hydrogen fractions are corrected for self-shielding using a fitting formula from radiative transfer simulations and for the presence of molecular hydrogen using an empirical or a theoretical prescription from the literature. The resulting neutral hydrogen fractions, M_(HI+H2)/M_star, agree with observations to better than 0.1 dex for galaxies with M_star between 10^10 and 10^11 M_Sun. Our fiducial, empirical H2 model based on gas pressure results in galactic HI mass fractions, M_HI/M_star, that agree with observations from the GASS survey to better than 0.3 dex, but the alternative theoretical H2 formula leads to a negative offset in M_HI/M_star of up to 0.5 dex. Visual inspection reveals that most HI disks in simulated HI-rich galaxies are vertically disturbed, plausibly due to recent accretion events. Many galaxies (up to 80 per cent) contain spuriously large HI holes, which are likely formed as a consequence of the feedback implementation in EAGLE. The HI mass-size relation of all simulated galaxies is close to (but 16 per cent steeper than) observed, and when only galaxies without large holes in the HI disc are considered, the agreement becomes excellent (better than 0.1 dex). The presence of large HI holes also makes the radial HI surface density profiles somewhat too low in the centre, at \Sigma_HI > 1 M_Sun pc^-2 (by a factor of <~ 2 compared to data from the Bluedisk survey). In the outer region (\Sigma_HI < 1 M_Sun pc^-2), the simulated profiles agree quantitatively with observations. Scaled by HI size, the simulated profiles of HI-rich (M_HI > 10^9.8 M_Sun) and control galaxies (10^9.1 M_Sun > M_HI > 10^9.8 M_Sun) follow each other closely, as observed. (Abridged)

We derive average flux corrections to the \texttt{Model} magnitudes of the Sloan Digital Sky Survey (SDSS) galaxies by stacking together mosaics of similar galaxies in bins of stellar mass and concentration. Extra flux is detected in the outer low surface brightness part of the galaxies, leading to corrections ranging from 0.05 to 0.32 mag for the highest stellar mass galaxies. We apply these corrections to the MPA-JHU (Max-Planck Institute for Astrophysics - John Hopkins University) stellar masses for a complete sample of half a million galaxies from the SDSS survey to derive a corrected galaxy stellar mass function at $z=0.1$ in the stellar mass range $9.5<\log(M_\ast/M_\odot)<12.0$. We find that the flux corrections and the use of the MPA-JHU stellar masses have a significant impact on the massive end of the stellar mass function, making the slope significantly shallower than that estimated by Li \& White (2009), but steeper than derived by Bernardi et al. (2013). This corresponds to a mean comoving stellar mass density of galaxies with stellar masses $\log(M_\ast/M_\odot) \ge 11.0$ that is a factor of 3.36 larger than the estimate by Li \& White (2009), but is 43\% smaller than reported by Bernardi et al. (2013).

We present a study exploring the nature and properties of the Circum-Galactic Medium (CGM) and its connection to the atomic gas content in the interstellar medium (ISM) of galaxies as traced by the HI 21cm line. Our sample includes 45 low-z (0.026-0.049) galaxies from the GALEX Arecibo SDSS Survey. Their CGM was probed via absorption in the spectra of background Quasi-Stellar Objects at impact parameters of 63 to 231kpc. The spectra were obtained with the Cosmic Origins Spectrograph aboard the Hubble Space Telescope. We detected neutral hydrogen (Ly$\alpha$ absorption-lines) in the CGM of 92% of the galaxies. We find the radial profile of the CGM as traced by the Ly$\alpha$ equivalent width can be fit as an exponential with a scale length of roughly the virial radius of the dark matter halo. We found no correlation between the orientation of sightline relative to the galaxy major axis and the Ly$\alpha$ equivalent width. The velocity spread of the circumgalactic gas is consistent with that seen in the atomic gas in the interstellar medium. We find a strong correlation (99.8% confidence) between the gas fraction (M(HI)/M*) and the impact-parameter-corrected Ly$\alpha$ equivalent width. This is stronger than the analogous correlation between corrected Ly$\alpha$ equivalent width and SFR/M* (97.5% confidence). These results imply a physical connection between the HI disk and the CGM, which is on scales an order-of-magnitude larger. This is consistent with the picture in which the HI disk is nourished by accretion of gas from the CGM.

Two explanations have been put forward to explain the observed conformity between the colours and specific star formation rates (SFR/$M_*$) of galaxies on large scales: 1) the formation times of their surrounding dark matter halos are correlated (commonly referred to as "assembly bias"), 2) gas is heated over large scales at early times, leading to coherent modulation of cooling and star formation between well-separated galaxies (commonly referred to as "pre-heating") . To distinguish between the pre-heating and assembly bias scenarios, we search for relics of energetic feedback events in the neighbourhood of central galaxies with different specific star formation rates. We find a significant excess of very high mass ($\log M_* > 11.3$) galaxies out to a distance of 2.5 Mpc around low SFR/$M_*$ central galaxies compared to control samples of higher SFR/$M_*$ central galaxies with the same stellar mass and redshift. We also find that very massive galaxies in the neighbourhood of low SFR/$M_*$ galaxies have much higher probability of hosting radio loud active galactic nuclei. The radio-loud AGN fraction in neighbours with $\log M_* > 11.3$ is four times higher around passive, non star-forming centrals at projected distances of 1 Mpc and two times higher at projected distances of 4 Mpc. Finally, we carry out an investigation of conformity effects in the recently publicly-released Illustris cosmological hydrodynamical simulation, which includes energetic input both from quasars and from radio mode accretion onto black holes. We do not find conformity effects of comparable amplitude on large scales in the simulations and we propose that gas needs to be pushed out of dark matter halos more efficiently at high redshifts.

We use a data set of nearby galaxies drawn from the HERACLES, ATLAS3D, and COLD GASS surveys to study variations in molecular gas depletion time (Tdep) in galaxy structures such as bulges, grand-design spiral arms, bars and rings. Molecular gas is traced by CO line emission and star formation rate (SFR) is derived using the combination of far-ultraviolet and mid-infrared (MIR) data. The contribution of old stars to MIR emission for the ATLAS3D sample is corrected using 2MASS K-band images. We apply a two-dimensional image decomposition algorithm to decompose galaxies into bulges and discs. Spiral arms, bars and rings are identified in the residual maps, and molecular gas depletion times are derived on a square grid of 1 kpc^2 size. In previous work, we showed that Tdep correlates strongly with specific star formation rate (sSFR). We now find that at a given sSFR, the bulge has shorter Tdep than the disc. The shift to shorter depletion times is most pronounced in the inner bulge (R < 0.1Re). Grids from galaxies with bars and rings are similar to those from galactic bulges in that they have reduced Tdep at a given sSFR. In contrast, the Tdep versus sSFR relation in the discs of galaxies with spiral arms is displaced to longer Tdep at fixed sSFR. We then show that the differences in the Tdep-sSFR relation for bulges, discs, arms, bars and rings can be linked to variations in "stellar", rather than gas surface density between different structures. Our best current predictor for Tdep, both globally and for 1 kpc grids, is given by Tdep= -0.36log(Sigma_SFR)-0.5log(Sigma_*)+5.87.

We present a systematic analysis of the rotation curves of 187 galaxies with masses greater than 10^10 M_sol, with atomic gas masses from the GALEX Arecibo Sloan Survey (GASS), and with follow-up long-slit spectroscopy from the MMT. Our analysis focuses on stellar rotation curves derived by fitting stellar template spectra to the galaxy spectra binned along the slit. In this way, we are able to obtain accurate rotation velocity measurements for a factor of 2 more galaxies than possible with the Halpha line. Galaxies with high atomic gas mass fractions are the most dark-matter dominated galaxies in our sample and have dark matter halo density profiles that are well fit by Navarro, Frenk & White profiles with an average concentration parameter of 10. The inner slopes and of the rotation curves correlate more strongly with stellar population age than with galaxy mass or structural parameters. At fixed stellar mass, the rotation curves of more actively star-forming galaxies have steeper inner slopes than less actively star-forming galaxies. The ratio between the galaxy specific angular momentum and the total specific angular momentum of its dark matter halo, R_j, correlates strongly with galaxy mass, structure and gas content. Low mass, disk-dominated galaxies with atomic gas mass fractions greater than 20% have median values of R_j of around 1, but massive, bulge-dominated galaxies have R_j=0.2-0.3. We argue that these trends can be understood in a picture where gas inflows triggered by disk instabilities lead to the formation of passive, bulge-dominated galaxies with low specific angular momentum.

Our work is based on the "Bluedisk" project, a program to map the neutral gas in a sample of 25 HI-rich spirals and a similar number of control galaxies with the Westerbork Synthesis Radio Telescope (WSRT). In this paper we focus on the HI properties of the galaxies in the environment of our targeted galaxies. In total, we extract 65 galaxies from the WSRT cubes with stellar masses between $10^8M_{\odot}$ and $10^{11}M_{\odot}$. Most of these galaxies are located on the same HI mass-size relation and "HI-plane" as normal spiral galaxies. We find that companions around HI-rich galaxies tend to be HI-rich as well and to have larger R90,HI/R50,HI. This suggests a scenario of "HI conformity", similar to the colour conformity found by Weinmann et al. (2006): galaxies tend to adopt the HI properties of their neighbours. We visually inspect the outliers from the HI mass-size relation and galaxies which are offset from the HI plane and find that they show morphological and kinematical signatures of recent interactions with their environment. We speculate that these outliers have been disturbed by tidal or ram-pressure stripping processes, or in a few cases, by accretion events.

We study the Kennicutt-Schmidt relation between average star formation rate and average cold gas surface density in the Hi dominated ISM of nearby spiral and dwarf irregular galaxies. We divide the galaxies into grid cells varying from sub-kpc to tens of kpc in size. Grid-cell measurements of low SFRs using H-alpha emission can be biased and scatter may be introduced because of non-uniform sampling of the IMF or because of stochastically varying star formation. In order to alleviate these issues, we use far-ultraviolet emission to trace SFR, and we sum up the fluxes from different bins with the same gas surface density to calculate the average $\Sigma_{SFR}$ at a given value of $\Sigma_{gas}$. We study the resulting Kennicutt-Schmidt relation in 400 pc, 1 kpc and 10 kpc scale grids in nearby massive spirals and in 400 pc scale grids in nearby faint dwarf irregulars. We find a relation with a power law slope of 1.5 in the HI-dominated regions for both kinds of galaxies. The relation is offset towards longer gas consumption timescales compared to the molecular hydrogen dominated centres of spirals, but the offset is an order-of-magnitude less than that quoted by earlier studies. Our results lead to the surprising conclusion that conversion of gas to stars is independent of metallicity in the HI dominated regions of star-forming galaxies. Our observed relations are better fit by a model of star formation based on thermal and hydrostatic equilibrium in the ISM, in which feedback driven turbulence sets the thermal pressure.

We re-analyze correlations between global molecular gas depletion time (Tdep) and galaxy parameters including stellar mass, specific star formation rate (sSFR), stellar mass surface density and concentration index. The analysis is based on the COLD GASS survey, which includes galaxies with stellar mass in the range 10^10 - 10^11.5 Msun with molecular gas mass estimates derived from CO(1-0) line measurements. We improve on previous work by Saintonge et al. (2011b) by estimating SFRs using the combination of GALEX FUV and WISE 22 micron data and by deriving Tdep within a fixed aperture set by the IRAM beam size. In our new study we find correlations with much smaller scatter. Dependences of the Tdep on galaxy structural parameters such as stellar mass surface density and concentration index are now weak or absent. Differences with previous work arise because dust extinction correlates strongly with galaxy structural parameters. We further demonstrate that the 'primary' global parameter correlation is between Tdep and sSFR; all other remaining correlations can be shown to be induced by this primary dependence. This implies that galaxies with high current-to-past-averaged star formation activity, will drain their molecular gas reservoir sooner. We then analyze molecular gas depletion times on 1-kpc scales in galactic disks using data from the HERACLES survey. There is remarkably good agreement between the global Tdep versus sSFR relation for the COLD GASS galaxies and that derived for 1 kpc scale grid regions in disks. The strong correlation between Tdep and sSFR extends continuously over a factor of 10 in Tdep from log(SFR/M*) = -11.5 to -9, i.e. from nearly quiescent patches of the disk to disk regions with very strong star formation. This leads to the conclusion that the local molecular gas depletion time in galactic disks is dependent on the local fraction of young-to-old stars.

We investigate the nature of HI-rich galaxies from the ALFALFA and GASS surveys, which are defined as galaxies in the top 10th percentile in atomic gas fraction at a given stellar mass. We analyze outer (R>1.5 Re) stellar populations for a subset of face-on systems using optical g-r versus r-z colour/colour diagrams. The results are compared with those from control samples that are defined without regard to atomic gas content, but are matched in redshift, stellar mass and structural parameters. HI-rich early-type (C>2.6) and late-type (C<2.6) galaxies are studied separately. When compared to the control sample, the outer stellar populations of the majority of HI-rich early-type galaxies are shifted in the colour/colour plane along a locus consistent with younger stellar ages, but similar metallicities. The outer colours of HI-rich late-type galaxies are much bluer in r-z than the HI-rich early types, and we infer that they have outer disks which are both younger and more metal-poor. We then proceed to analyze the galaxy environments of HI-rich galaxies on scales of 500 kpc. HI-rich early-type galaxies with low (log M* < 10.5) stellar masses differ significantly from the control sample in that they are more likely to be central rather than satellite systems. Their satellites are also less massive and have younger stellar populations. Similar, but weaker effects are found for low mass HI-rich late-type galaxies. In addition, we find that the satellites of HI-rich late-types exhibit a greater tendency to align along the major axis of the primary. No environmental differences are found for massive (log M* > 10.5) HI-rich galaxies, regardless of type.

Based on large optical and mid-infrared (IR) surveys, we investigate the relation between nuclear activity in local Seyfert 2 galaxies and galaxy interactions using a statistical neighbour counting technique. At the same level of host galaxy star formation (SF), we find that active galactic nuclei (AGNs) with stronger [OIII] emission lines do not show an excess of near neighbours, while AGNs with stronger mid-IR emission do have more near neighbours within a projected distance of 100 kpc. The excess neighbour count increases with decreasing projected radius. These results suggest a phase of torus formation during galaxy interactions.

We present an overview of a new integral field spectroscopic survey called MaNGA (Mapping Nearby Galaxies at Apache Point Observatory), one of three core programs in the fourth-generation Sloan Digital Sky Survey (SDSS-IV) that began on 2014 July 1. MaNGA will investigate the internal kinematic structure and composition of gas and stars in an unprecedented sample of 10,000 nearby galaxies. We summarize essential characteristics of the instrument and survey design in the context of MaNGA's key science goals and present prototype observations to demonstrate MaNGA's scientific potential. MaNGA employs dithered observations with 17 fiber-bundle integral field units that vary in diameter from 12" (19 fibers) to 32" (127 fibers). Two dual-channel spectrographs provide simultaneous wavelength coverage over 3600-10300 A at R~2000. With a typical integration time of 3 hr, MaNGA reaches a target r-band signal-to-noise ratio of 4-8 (per A, per 2" fiber) at 23 AB mag per sq. arcsec, which is typical for the outskirts of MaNGA galaxies. Targets are selected with stellar mass greater than 1e9 Msun using SDSS-I redshifts and i-band luminosity to achieve uniform radial coverage in terms of the effective radius, an approximately flat distribution in stellar mass, and a sample spanning a wide range of environments. Analysis of our prototype observations demonstrates MaNGA's ability to probe gas ionization, shed light on recent star formation and quenching, enable dynamical modeling, decompose constituent components, and map the composition of stellar populations. MaNGA's spatially resolved spectra will enable an unprecedented study of the astrophysics of nearby galaxies in the coming 6 yr.

We study the stellar haloes of galaxies out to 70-100 kpc as a function of stellar mass and galaxy type by stacking aligned $r$ and $g$ band images from a sample of 45508 galaxies from SDSS DR9 in the redshift range $0.06\,\le\,z\,\le\,0.1$ and in the mass range $10^{10.0} M_{\odot} < M_{*} < 10^{11.4} M_{\odot}$r. We derive surface brightness profiles to a depth of almost $\mu_r \sim 32 \,\mathrm{mag\,arcsec}^{-2}$. We find that the ellipticity of the stellar halo is a function of galaxy stellar mass and that the haloes of high concentration ($C > 2.6$) galaxies are more elliptical than those of low concentration ($C < 2.6$) galaxies. The $g$-$r$ colour profile of high concentration galaxies reveals that the $g$-$r$ colour of the stellar population in the stellar halo is bluer than in the main galaxy, and the colour of the stellar halo is redder for higher mass galaxies. We further demonstrate that the full two-dimensional surface intensity distribution of our galaxy stacks can only be fit through multi-component S\'{e}rsic models. Using the fraction of light in the outer component of the models as a proxy for the fraction of accreted stellar light, we show that this fraction is a function of stellar mass and galaxy type. For high concentration galaxies, the fraction of accreted stellar light rises from $30\%$ to $70\%$ for galaxies in the stellar mass range from $10^{10.0} M_{\odot}$ to $10^{11.4} M_{\odot}$. The fraction of accreted light is much smaller in low concentration systems, increasing from $2\%$ to $25\%$ over the same mass range. This work provides important constraints for the theoretical understanding of the formation of stellar haloes of galaxies.

We analyze the radial distribution of HI gas for 23 disk galaxies with unusually high HI content from the Bluedisk sample, along with a similar-sized sample of "normal" galaxies. We propose an empirical model to fit the radial profile of the HI surface density, an exponential function with a depression near the center. The radial HI surface density profiles are very homogeneous in the outer regions of the galaxy; the exponentially declining part of the profile has a scale-length of $\sim 0.18$ R1, where R1 is the radius where the column density of the HI is 1 M$_{\odot}$ pc$^{-2}$. This holds for all galaxies, independent of their stellar or HI mass. The homogenous outer profiles, combined with the limited range in HI surface density in the non-exponential inner disk, results in the well-known tight relation between HI size and HI mass. By comparing the radial profiles of the HI-rich galaxies with those of the control systems, we deduce that in about half the galaxies, most of the excess gas lies outside the stellar disk, in the exponentially declining outer regions of the HI disk. In the other half, the excess is more centrally peaked. We compare our results with existing smoothed-particle hydrodynamical simulations and semi-analytic models of disk galaxy formation in a $\Lambda$ Cold Dark Matter universe. Both the hydro simulations and the semi-analytic models reproduce the HI surface density profiles and the HI size-mass relation without further tuning of the simulation and model inputs. In the semi-analytic models, the universal shape of the outer HI radial profiles is a consequence of the {\em assumption} that infalling gas is always distributed exponentially. The conversion of atomic gas to molecular form explains the limited range of HI surface densities in the inner disk. These two factors produce the tight HI mass-size relation.

We have used 4000 \AA\ break and H$\delta_A$ indices in combination with SFR/M* derived from emission line flux measurements, to constrain the recent star formation histories of galaxies with stellar masses in the range 10**8-10**10 Msol. The fraction of the total SFR density in galaxies with ongoing bursts is a strong function of stellar mass, declining from 0.85 at a stellar mass of 10**8 Msol to 0.25 for galaxies with M* = 10**10 Msol. Low mass galaxies are not all young. The distribution of half mass formation times for galaxies with stellar masses less than 10**9 Msol is broad, spanning the range 1-10 Gyr. The peak-to-trough variation in star formation rate among the bursting population ranges lies in the range 10-25. In low mass galaxies, the average duration of the burst is comparable to the dynamical time of the galaxy. Galaxy structure is correlated with estimated burst mass fraction, but in different ways in low and in high mass galaxies. High mass galaxies with large burst mass fractions are more centrally concentrated, indicating that bulge formation is at work. In low mass galaxies, stellar surface mass densities decrease as a function of F_{burst}. These results are in good agreement with the observational predictions of Teyssier et al (2013) and lend further credence to the idea that the cuspy halo problem can be solved by energy input from multiple starbursts over the lifetime of the galaxy. We note that there is no compelling evidence for IMF variations in the population of star-forming galaxies in the local Universe.

We investigate whether gradual dilution of the gas in some elliptical galaxies is the cause of a positive correlation between star formation rate (SFR) and gas-phase metallicity (Zg) at high stellar mass (M*) in the local Universe. To do this, two classes of massive (M* >= 10^10.5 Msun) galaxy are selected from the Sloan Digital Sky Survey (SDSS) and the Munich semi-analytic model of galaxy formation, L-Galaxies. The first class is selected by high specific star formation rates (sSFR) and high Zg, and the second class by low sSFR and low Zg. These criteria roughly distinguish disc-dominant galaxies from metal-poor, elliptical galaxies. In the semi-analytic model, the second class of galaxies obtain low sSFR and low Zg due to gradual dilution of the interstellar medium by accretion of metal-poor gas via infalling clumps and low-mass satellites. This occurs after a merger-induced starburst and the associated supernova feedback have quenched most of the original gas reservoir. A number of signatures of this evolution are present in these model galaxies at z=0, including low gas fractions, large central black holes, elliptical morphologies, old ages, and importantly, low (Zg-Z*) indicating dilution after star formation. Remarkably, all of these properties are also found in low-sSFR, low-Zg, massive galaxies in the SDSS-DR7. This provides strong, indirect evidence that some elliptical galaxies are undergoing gradual dilution after a gas-rich merger in the local Universe. This dilution scenario also explains the positive correlation between SFR and Zg measured in high-M* galaxies, and therefore has consequences for the local fundamental metallicity relation (FMR), which assumes a weak anti-correlation between SFR and Zg above 10^10.5 Msun.

We analyze a sample of 30,000 nearby obscured AGNs with optical spectra from SDSS and mid-IR photometry from WISE. Our aim is to investigate the AGN host galaxy properties with mid-IR luminosities as AGN activity indicator, and to compare with previous studies based on [OIII] emission lines. First we find that the [3.4] - [4.6] colour has weak dependence on host stellar age, but strong dependence on AGN activity. We then use a "pair-matching" technique to subtract the host 4.6 micron contribution. By combining Seyferts with a sample of SDSS quasars at z < 0.7, we show that the [OIII] and the intrinsic AGN 4.6 micron luminosities correlate roughly linearly over 4 orders of magnitude, but with substantial scatter. We also compare the partition functions of the total integrated 4.6 micron and [OIII] line luminosities from Seyferts and a sub-population of LINERs with significant nuclear 4.6 micron emission, as function of a variety of host galaxy properties, finding that they are identical. We conclude, therefore, that [OIII] as an AGN indicator shows no particular biases as compared to the 4.6 micron luminosity. Our results also demonstrate that some LINERs do fit in with the expectations of the simple Unified Model.

We have combined the semi-analytic galaxy formation model of Guo et al. (2011) with the particle-tagging technique of Cooper et al. (2010) to predict galaxy surface brightness profiles in a representative sample of ~1900 massive dark matter haloes (10^12--10^14 M_sol) from the Millennium II Lambda-CDM N-body simulation. Here we present our method and basic results focusing on the outer regions of galaxies, consisting of stars accreted in mergers. These simulations cover scales from the stellar haloes of Milky Way-like galaxies to the 'cD envelopes' of groups and clusters, and resolve low surface brightness substructure such as tidal streams. We find that the surface density of accreted stellar mass around the central galaxies of dark matter haloes is well described by a Sersic profile, the radial scale and amplitude of which vary systematically with halo mass (M_200). The total stellar mass surface density profile breaks at the radius where accreted stars start to dominate over stars formed in the galaxy itself. This break disappears with increasing M_200 because accreted stars contribute more of the total mass of galaxies, and is less distinct when the same galaxies are averaged in bins of stellar mass, because of scatter in the relation between M_star and M_200. To test our model we have derived average stellar mass surface density profiles for massive galaxies at z~0.08 by stacking SDSS images. Our model agrees well with these stacked profiles and with other data from the literature, and makes predictions that can be more rigorously tested by future surveys that extend the analysis of the outer structure of galaxies to fainter isophotes. We conclude that it is likely that the outer structure of the spheroidal components of galaxies is largely determined by collisionless merging during their hierarchical assembly

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.

We use long-slit spectroscopy from Moran et al. to study the radial dependence of the recent star formation histories of nearby galaxies with stellar masses greater than 10^10M_sun. We fit stellar population models to the combination of SSFR, D4000 and Hdelta_A and show that many galaxies have Balmer absorption line equivalent widths that require recent short-lived episodes or bursts of star formation. The fraction of galaxies that have experienced episodic rather than continuous star formation is highest for late-type galaxies with low stellar masses. In these systems, bursts occur both in the inner and outer regions of the galaxy. The fraction of stars formed in a single burst episode is typically around 15% of the total stellar mass in the inner regions of the galaxy and around 5% of the mass in the outer regions. When we average over the population, we find that such bursts contribute around a half of the total mass in stars formed in the last 2 Gyr. In massive galaxies, bursts occur predominantly in the outer disk. Around a third of all massive, bulge-dominated galaxies have experienced recent star formation episodes that are fully confined to their outer (R > 0.7R_90) regions. The fraction of stars formed in a single episode is only 2 - 3 % of the underlying stellar mass, but such bursts contribute nearly all the stellar mass formed in the last 2 Gyr. Recent star formation in outer disks is strongly correlated with the global atomic gas fraction of the galaxy, but not its global molecular gas fraction. We suggest that outer episodic star formation is triggered by gas accretion events.

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.

We present a method to estimate the total gas column density, dust-to-gas and dust-to-metal ratios of distant galaxies from rest-frame optical spectra. The technique exploits the sensitivity of certain optical lines to changes in depletion of metals onto dust grains and uses photo-ionization models to constrain these physical ratios along with the metallicity (Z) and dust column density. We compare our gas column density estimates with HI and CO gas mass estimates in nearby galaxies to show that we recover their total gas mass surface density to within a factor of 2 up to a total surface gas mass density of ~75 Mo/pc^2. Our technique is independent of the conversion factor of CO to H2 and we show that a Z-dependent XCO is required to achieve good agreement between our measurements and that provided by CO and HI. However we also show that our method can not be reliably aperture corrected to total gas mass. We calculate dust-to-gas ratios for all star-forming galaxies in the SDSS DR7 and show the resulting dependence on Z agrees well with the trend inferred from modelling of the dust emission of nearby galaxies using FIR data. We also present estimates of the variation of the dust-to-metal ratio with Z and show that this is poorly constrained at Z<~0.5Z_sun. We end with a study of the inventory of gas in the central regions, defined both in terms of a fixed physical radius and as a fixed fraction of the half-light radius, of ~70,000 star-forming galaxies from the Sloan Digital Sky Survey. We show that their central gas content and gas depletion times are not accurately predicted by a single parameter, but in agreement with recent studies we find that a combination of the M* and some measure of central concentration provides a good predictor of gas content in galaxies. We also identify a population of galaxies with low surface densities of stars and very long gas depletion times.

We introduce the "Bluedisk" project, a large program at the Westerbork Synthesis Radio Telescope (WSRT) that has mapped the HI in a sample of 23 nearby galaxies with unusually high HI mass fractions, along with a similar-sized sample of control galaxies matched in stellar mass, size, inclination and redshift. This paper presents the sample selection, observational set-up, data reduction strategy, and a first analysis of the sizes and structural properties of the HI disks. We find that the HI-rich galaxies lie on the same HI mass versus HI size relation as normal spiral galaxies, extending it to total HI masses of $2 \times 10^{10} M_{\odot}$ and radii R1 of $\sim 100$ kpc (where R1 is defined as the radius where the HI column density reaches 1 $M_{\odot}$ pc$^{-2}$). HI-rich galaxies have significantly larger values of HI-to-optical size ratio at fixed stellar mass, concentration index, stellar and star formation rate surface density compared to the control sample. The disks of HI-rich galaxies are also significantly more clumpy (i.e. have higher HI Gini and $\Delta$Area coefficient) than those of normal spirals. There is no evidence that the disks of HI-rich galaxies are more disturbed: HI-rich galaxies exhibit no difference with respect to control samples in their distributions of HI asymmetry indices or optical/HI disk position angle differences. In fact, the center of the HI distribution corresponds more closely with the center of the optical light in the HI-rich galaxies than in the controls. All these results argue against a scenario in which new gas has been brought in by mergers. It is possible that they may be more consistent with cooling from a surrounding quasi-static halo of warm/hot gas.

We have compared the observed distribution of galaxies in the plane of specific star formation rate versus stellar mass with the predictions of the Garching semi-analytic model at redshifts 0, 1 and 2. The goal is to test whether the implementation of radio mode AGN feedback, which is responsible for terminating the formation of stars in high mass galaxies, provides an adequate match to current high-redshift observations. The fraction of quenched galaxies as a function of stellar mass in the models is in good agreement with data at z=0 and z=1. By z=2, there are too few quenched galaxies with low stellar masses in the models. At z=2, the population of galaxies with no ongoing star formation is clearly separated from the `main sequence' of star-forming galaxies in the data. This is not found in the models, because z=2 galaxies with stellar masses less than 10^11 solar masses are predicted to host black holes with relatively low masses (less than 10^8 solar masses). The current implementation of radio mode feedback from such black holes reduces the cooling rates from the surrounding halo, but does not generate sufficient energy to stop star formation entirely. We suggest that the models may be brought into better agreement with the data if black hole growth is triggered by disc instabilities in addition to major mergers, and if feedback mechanisms associated with the formation of galactic bulges act to quench star formation in galaxies.

The Square Kilometre Array (SKA) is the most ambitious radio telescope ever planned. With a collecting area of about a square kilometre, the SKA will be far superior in sensitivity and observing speed to all current radio facilities. The scientific capability promised by the SKA and its technological challenges provide an ideal base for interdisciplinary research, technology transfer, and collaboration between universities, research centres and industry. The SKA in the radio regime and the European Extreme Large Telescope (E-ELT) in the optical band are on the roadmap of the European Strategy Forum for Research Infrastructures (ESFRI) and have been recognised as the essential facilities for European research in astronomy. This "White Paper" outlines the German science and R&D interests in the SKA project and will provide the basis for future funding applications to secure German involvement in the Square Kilometre Array.

We analyze galaxies in 300 nearby groups and clusters identified in the Sloan Digital Sky Survey using a photometric gas mass indicator that is useful for estimating the degree to which the interstellar medium of a cluster galaxy has been depleted. We study the radial dependence of inferred gas mass fractions for galaxies of different stellar masses and stellar surface densities. At fixed clustercentric distance and at fixed stellar mass, lower density galaxies are more strongly depleted of their gas than higher density galaxies. An analysis of depletion trends in the two-dimensional plane of stellar mass $M_*$ and stellar mass surface density $\mu_*$ reveals that gas depletion at fixed clustercentric radius is much more sensitive to the density of a galaxy than to its mass. We suggest that low density galaxies are more easily depleted of their gas, because they are more easily affected by ram-pressure and/or tidal forces. We also look at the dependence of our gas fraction/radius relations on the velocity dispersion of the cluster, finding no clear systematic trend.

The Baryon Oscillation Spectroscopic Survey (BOSS) is designed to measure the scale of baryon acoustic oscillations (BAO) in the clustering of matter over a larger volume than the combined efforts of all previous spectroscopic surveys of large scale structure. BOSS uses 1.5 million luminous galaxies as faint as i=19.9 over 10,000 square degrees to measure BAO to redshifts z<0.7. Observations of neutral hydrogen in the Lyman alpha forest in more than 150,000 quasar spectra (g<22) will constrain BAO over the redshift range 2.15<z<3.5. Early results from BOSS include the first detection of the large-scale three-dimensional clustering of the Lyman alpha forest and a strong detection from the Data Release 9 data set of the BAO in the clustering of massive galaxies at an effective redshift z = 0.57. We project that BOSS will yield measurements of the angular diameter distance D_A to an accuracy of 1.0% at redshifts z=0.3 and z=0.57 and measurements of H(z) to 1.8% and 1.7% at the same redshifts. Forecasts for Lyman alpha forest constraints predict a measurement of an overall dilation factor that scales the highly degenerate D_A(z) and H^{-1}(z) parameters to an accuracy of 1.9% at z~2.5 when the survey is complete. Here, we provide an overview of the selection of spectroscopic targets, planning of observations, and analysis of data and data quality of BOSS.

The observed correlation between star-formation in central galaxies and in their neighbours (a phenomenon dubbed galactic conformity) is in need of a convincing physical explanation. We use a volume-limited sample of galaxies with redshifts less than 0.03 drawn from the SDSS DR7 to investigate the scale dependence of the effect and how it changes as a function of the mass of the central galaxy. Conformity extends over a central galaxy stellar mass range spanning two orders of magnitude. In central galaxies with masses less than 10^10 M_sun, conformity extends out to scales in excess of 4 Mpc, well beyond the virial radii of their dark matter halos. For low mass central galaxies, large-scale conformity with neighbours is only seen when the centrals have low star formation rate or gas content. In contrast, at high stellar masses, conformity with neighbours applies in the gas-rich regime and is clearly confined to scales comparable to the virial radius of the dark matter halo of the central galaxy. Our analysis of a mock catalogue from the Guo et al (2011) models shows that conformity-like effects are much weaker than observed, and apply only to the low SFR/M* tail of neighbouring galaxies. In the models, the median and the upper percentiles of the SFR/M* distribution remain almost unchanged, which is in contradiction with the data. Conformity between low-mass, gas-poor central galaxies and their distant neighbours cannot be explained within the framework of halo occupation distribution (HOD) models. It is likely a signature of pre-heating of the intergalactic gas at an earlier epoch. The smaller-scale conformity between high-mass, gas-rich central galaxies and their close neighbours may be a signature of ongoing gas accretion onto central galaxies in a minority of massive dark matter halos.

It is well known that both the star formation rate and the cold gas content of a galaxy depend on the local density out to distances of a few Megaparsecs. In this paper, we compare the environmental density dependence of the atomic gas mass fractions of nearby galaxies with the density dependence of their central and global specific star formation rates. We stack HI line spectra extracted from the Arecibo Legacy Fast ALFA survey centered on galaxies with UV imaging from GALEX and optical imaging/spectroscopy from SDSS. We use these stacked spectra to evaluate the mean atomic gas mass fraction of galaxies in bins of stellar mass and local density. For galaxies with stellar masses less than 10^10.5 M_sun, the decline in mean atomic gas mass fraction with density is stronger than the decline in mean global and central specific star formation rate. The same conclusion does not hold for more massive galaxies. We interpret our results as evidence for ram-pressure stripping of atomic gas from the outer disks of low mass satellite galaxies. We compare our results with the semi-analytic recipes of Guo et al. (2011) implemented on the Millennium II simulation. These models assume that only the diffuse gas surrounding satellite galaxies is stripped, a process that is often termed "strangulation". We show that these models predict relative trends in atomic gas and star formation that are in disagreement with observations. We use mock catalogues generated from the simulation to predict the halo masses of the HI-deficient galaxies in our sample. We conclude that ram-pressure stripping is likely to become effective in dark matter halos with masses greater than 10^13 M_sun.

We analyze the optical spectra of massive (log M*/Msun > 11.4) radio-loud galaxies at z~0.2 and z~0.6. By comparing stellar population parameters of these radio-loud samples with radio-quiet control samples, we investigate how the presence of a radio-emitting jet relates to the recent star formation history of the host galaxy. We also investigate how the emission-line properties of the radio galaxies evolve with redshift by stacking their spectra. Our main results are the following. (1) Both at low and at high redshift, half as many radio-loud as radio-quiet galaxies have experienced significant star formation in the past Gyr. (2) The Balmer absorption line properties of massive galaxies that have experienced recent star formation show that star formation occurred as a burst in many of these systems. (3) Both the radio and the emission-line luminosity of radio AGN evolve significantly with redshift. However, radio galaxies with similar stellar population parameters, have similar emission-line properties both at high- and at low-redshift. These results suggest that massive galaxies experience cyclical episodes of gas accretion, star formation and black hole growth, followed by the production of a radio jet that shuts down further activity. The behaviour of galaxies with log M*/Msun > 11.4 is the same at z = 0.6 as it is at z = 0.2, except that higher redshift galaxies experience more star formation and black hole growth and produce more luminous radio jets during each accretion cycle.

The Sloan Digital Sky Survey III (SDSS-III) presents the first spectroscopic data from the Baryon Oscillation Spectroscopic Survey (BOSS). This ninth data release (DR9) of the SDSS project includes 535,995 new galaxy spectra (median z=0.52), 102,100 new quasar spectra (median z=2.32), and 90,897 new stellar spectra, along with the data presented in previous data releases. These spectra were obtained with the new BOSS spectrograph and were taken between 2009 December and 2011 July. In addition, the stellar parameters pipeline, which determines radial velocities, surface temperatures, surface gravities, and metallicities of stars, has been updated and refined with improvements in temperature estimates for stars with T_eff<5000 K and in metallicity estimates for stars with [Fe/H]>-0.5. DR9 includes new stellar parameters for all stars presented in DR8, including stars from SDSS-I and II, as well as those observed as part of the SDSS-III Sloan Extension for Galactic Understanding and Exploration-2 (SEGUE-2). The astrometry error introduced in the DR8 imaging catalogs has been corrected in the DR9 data products. The next data release for SDSS-III will be in Summer 2013, which will present the first data from the Apache Point Observatory Galactic Evolution Experiment (APOGEE) along with another year of data from BOSS, followed by the final SDSS-III data release in December 2014.

We present the second data release from the GALEX Arecibo SDSS Survey (GASS), an ongoing large Arecibo program to measure the HI properties for an unbiased sample of ~1000 galaxies with stellar masses greater than 10^10 Msun and redshifts 0.025<z<0.05. GASS targets are selected from the Sloan Digital Sky Survey (SDSS) spectroscopic and Galaxy Evolution Explorer (GALEX) imaging surveys, and are observed until detected or until a gas mass fraction limit of a few per cent is reached. This second data installment includes new Arecibo observations of 240 galaxies, and marks the 50% of the complete survey. We present catalogs of the HI, optical and ultraviolet parameters for these galaxies, and their HI-line profiles. Having more than doubled the size of the sample since the first data release, we also revisit the main scaling relations of the HI mass fraction with galaxy stellar mass, stellar mass surface density, concentration index, and NUV-r color, as well as the gas fraction plane introduced in our earlier work.

We examine how the atomic and molecular gas components of galaxies evolve to higher redshifts using the semi-analytic galaxy formation models of Fu et al. (2010) in which we track the surface density profiles of gas in disks. We adopt two different prescriptions based either on gas surface density and metallicity, or on interstellar pressure, to compute the molecular fraction as a function of radius in each disk. We demonstrate that the adopted star formation law determines how the {\em balance} between gas, stars and metals changes with time in the star-forming galaxy population, but does not influence the total mass in stars formed into galaxies at redshifts below $z\sim 2.5$. The redshift evolution of the mass-metallicity relation places strong constraints on the timescale over which cold gas is converted into stars in high redshift galaxies, and favours models where this remains constant at values around 1-2 Gyr. Future observations of the evolution of the average molecular-to-atomic gas ratio in galaxies as a function of stellar mass and redshift will constrain models of the atomic-to-molecular transition.

We introduce a new photometric estimator of the HI mass fraction (M_HI/M_*) in local galaxies, which is a linear combination of four parameters: stellar mass, stellar surface mass density, NUV-r colour, and g-i colour gradient. It is calibrated using samples of nearby galaxies (0.025<z<0.05) with HI line detections from the GASS and ALFALFA surveys, and it is demonstrated to provide unbiased M_HI/M_* estimates even for HI-rich galaxies. We apply this estimator to a sample of ~24,000 galaxies from the SDSS/DR7 in the same redshift range. We then bin these galaxies by stellar mass and HI mass fraction and compute projected two point cross-correlation functions with respect to a reference galaxy sample. Results are compared with predictions from current semi-analytic models of galaxy formation. The agreement is good for galaxies with stellar masses larger than 10^10 M_sun, but not for lower mass systems. We then extend the analysis by studying the bias in the clustering of HI-poor or HI-rich galaxies with respect to galaxies with normal HI content on scales between 100 kpc and ~5 Mpc. For the HI-deficient population, the strongest bias effects arise when the HI-deficiency is defined in comparison to galaxies of the same stellar mass and size. This is not reproduced by the semi-analytic models, where the quenching of star formation in satellites occurs by "starvation" and does not depend on their internal structure. HI-rich galaxies with masses greater than 10^10 M_sun are found to be anti-biased compared to galaxies with "normal" HI content. Interestingly, no such effect is found for lower mass galaxies.

We analyze the role of bars in the build-up of central mass concentrations in massive, disk galaxies. Our parent sample consists of 3757 face-on disk galaxies with redshifts between 0.01 and 0.05, selected from the seventh Data Release of the Sloan Digital Sky Survey. 1555 galaxies with bars are identified using position angle and ellipticity profiles of the $i$-band light. We compare the ratio of the specific star formation rate measured in the 1-3 kpc central region of the galaxy to that measured for the whole galaxy. Galaxies with strong bars have centrally enhanced star formation; the degree of enhancement depends primarily on the ellipticity of the bar, and not on the size of the bar or on the mass or structure of the host galaxy. The fraction of galaxies with strong bars is highest at stellar masses greater than $3 \times 10^{10} M_{\odot}$, stellar surface densities less than $3 \times 10^8 M_{\odot}$ and concentration indices less than 2.5. In this region of parameter space, galaxies with strong bars either have enhanced central star formation rates, or star formation that is {\em suppressed} compared to the mean. This suggests that bars may play a role in the eventual quenching of star formation in galaxies. Only 50% of galaxies with strongly concentrated star formation have strong bars, indicating that other processes such as galaxy interactions also induce central star-bursts. We also find that the ratio of the size of the bar to that of the disk depends mainly on the colour of the galaxy, suggesting that the growth and destruction of bars are regulated by gas accretion, as suggested by simulations.

We have used the European VLBI Network (EVN) to observe a sample of Lyman Break Analogs (LBAs), nearby (z < 0.3) galaxies with properties similar to the more distant Lyman Break Galaxies (LBGs). The study of LBGs may help define the feed-back relationship between black holes (BHs) and their host galaxies. Previous VLA observations have shown that the kpc-scale radio emission from LBAs is dominated by starbursts. The main targets of this VLBI experiment were selected because they possessed emission-line properties between starbursts and Type 2 (obscured) AGN. Eight targets (three star-forming LBAs, four composite LBAs, and one Type 1 AGN) were observed at 5 GHz, four of which were also observed at 1.7 GHz. One star-forming LBA and one composite LBA were detected above 5 \sigma at 1.7 GHz (only), while the AGN was detected at 5 GHz. In both LBAs, the radio luminosity (LR) exceeds that expected from supernovae (remnants) based on a comparison with Arp220, Arp229A and Mrk273, by factors of 2 - 8. The composite LBA exhibits a compact core emitting around 10% of the VLA flux density. The high Tb of 3.5E7 K and excess core L_R with respect to the L_R/L_X relation of radio-quiet AGN indicate that this LBA possesses an obscured AGN (MBH ~ E5-E7 M_sun). While weak AGN may co-exist with the starbursts as shown in at least one of the LBAs, their contribution to the total radio flux is fairly minimal. Our results show that the detection of such weak AGN presents a challenge at radio, X-ray and optical emission-line wavelengths at z ~ 0.2, indicating the great difficulties that need to be overcome in order to study similar processes at high redshift when these types of galaxies were common.

We compare the semi-analytic models of galaxy formation of Fu et al. (2010), which track the evolution of the radial profiles of atomic and molecular gas in galaxies, with gas fraction scaling relations derived from the COLD GASS survey (Saintonge et al 2011). The models provide a good description of how condensed baryons in galaxies with gas are partitioned into stars, atomic and molecular gas as a function of galaxy stellar mass and surface density. The models do not reproduce the tight observed relation between stellar surface density and bulge-to-disk ratio for this population. We then turn to an analysis of the"quenched" population of galaxies without detectable cold gas. The current implementation of radio-mode feedback in the models disagrees strongly with the data. In the models, gas cooling shuts down in nearly all galaxies in dark matter halos above a mass of 10**12 M_sun. As a result, stellar mass is the observable that best predicts whether a galaxy has little or no neutral gas. In contrast, our data show that quenching is largely independent of stellar mass. Instead, there are clear thresholds in bulge-to-disk ratio and in stellar surface density that demarcate the location of quenched galaxies. We propose that processes associated with bulge formation play a key role in depleting the neutral gas in galaxies and that further gas accretion is suppressed following the formation of the bulge, even in dark matter halos of low mass.

We have obtained long-slit spectra of 174 star-forming galaxies with stellar masses greater than 10^10 M_\odot from the GALEX Arecibo SDSS (GASS) survey. These galaxies have both HI and H_2 mass measurements. The average metallicity profile is strikingly flat out to R_90, the radius enclosing 90% of the r-band light. Metallicity profiles which decline steadily with radius are found primarily for galaxies in our sample with low stellar mass (Log(M_*)<10.2), concentration, and/or mean stellar mass density. Beyond ~R_90, however, around 10 percent of the galaxies in our sample exhibit a sharp downturn in metallicity. Remarkably, we find that the magnitude of the outer metallicity drop is well correlated with the total HI content of the galaxy (measured as f_HI=M_HI/M_*). We examine the radial profiles of stellar population ages and star formation rate densities, and conclude that the galaxies with largest outer metallicity drops are actively growing their stellar disks, with mass doubling times across the whole disk only one third as long as a typical GASS galaxy. We also describe a correlation between local stellar mass density and metallicity, which is valid across all galaxies in our sample. We argue that much of the recent stellar mass growth at the edges of these galaxies can be linked to the accretion or radial transport of relatively pristine gas from beyond the galaxies' stellar disks.

We use principal component analysis (PCA) to estimate stellar masses, mean stellar ages, star formation histories (SFHs), dust extinctions and stellar velocity dispersions for ~290,000 galaxies with stellar masses greater than $10^{11}Msun and redshifts in the range 0.4<z<0.7 from the Baryon Oscillation Spectroscopic Survey (BOSS). We find the fraction of galaxies with active star formation first declines with increasing stellar mass, but then flattens above a stellar mass of 10^{11.5}Msun at z~0.6. This is in striking contrast to z~0.1, where the fraction of galaxies with active star formation declines monotonically with stellar mass. At stellar masses of 10^{12}Msun, therefore, the evolution in the fraction of star-forming galaxies from z~0.6 to the present-day reaches a factor of ~10. When we stack the spectra of the most massive, star-forming galaxies at z~0.6, we find that half of their [OIII] emission is produced by AGNs. The black holes in these galaxies are accreting on average at ~0.01 the Eddington rate. To obtain these results, we use the stellar population synthesis models of Bruzual & Charlot (2003) to generate a library of model spectra with a broad range of SFHs, metallicities, dust extinctions and stellar velocity dispersions. The PCA is run on this library to identify its principal components over the rest-frame wavelength range 3700-5500A. We demonstrate that linear combinations of these components can recover information equivalent to traditional spectral indices such as the 4000A break strength and HdA, with greatly improved S/N. This method is able to recover physical parameters such as stellar mass-to-light ratio, mean stellar age, velocity dispersion and dust extinction from the relatively low S/N BOSS spectra. We examine the sensitivity of our stellar mass estimates to the input parameters in our model library and the different stellar population synthesis models.

We present dynamical scaling relations for a homogeneous and representative sample of ~500 massive galaxies, selected only by stellar mass (>10^10 Msun) and redshift (0.025<z<0.05) as part of the ongoing GALEX Arecibo SDSS Survey. We compare baryonic Tully-Fisher (BTF) and Faber-Jackson (BFJ) relations for this sample, and investigate how galaxies scatter around the best fits obtained for pruned subsets of disk-dominated and bulge-dominated systems. The BFJ relation is significantly less scattered than the BTF when the relations are applied to their maximum samples, and is not affected by the inclination problems that plague the BTF. Disk-dominated, gas-rich galaxies systematically deviate from the BFJ relation defined by the spheroids. We demonstrate that by applying a simple correction to the stellar velocity dispersions that depends only on the concentration index of the galaxy, we are able to bring disks and spheroids onto the same dynamical relation -- in other words, we obtain a generalized BFJ relation that holds for all the galaxies in our sample, regardless of morphology, inclination or gas content, and has a scatter smaller than 0.1 dex. We find that disks and spheroids are offset in the stellar dispersion-size relation, and that the offset is removed when corrected dispersions are used instead. The generalized BFJ relation represents a fundamental correlation between the global dark matter and baryonic content of galaxies, which is obeyed by all (massive) systems regardless of morphology. [abridged]

We present measurements of projected autocorrelation functions w_p(r_p) for the stellar mass of galaxies and for their light in the U, B and V bands, using data from the third data release of the DEEP2 Galaxy Redshift Survey and the final data release of the Sloan Digital Sky Survey (SDSS). We investigate the clustering bias of stellar mass and light by comparing these to projected autocorrelations of dark matter estimated from the Millennium Simulations (MS) at z=1 and 0.07, the median redshifts of our galaxy samples. All of the autocorrelation and bias functions show systematic trends with spatial scale and waveband which are impressively similar at the two redshifts. This shows that the well-established environmental dependence of stellar populations in the local Universe is already in place at z=1. The recent MS-based galaxy formation simulation of Guo et al. (2011) reproduces the scale-dependent clustering of luminosity to an accuracy better than 30% in all bands and at both redshifts, but substantially overpredicts mass autocorrelations at separations below about 2 Mpc. Further comparison of the shapes of our stellar mass bias functions with those predicted by the model suggests that both the SDSS and DEEP2 data prefer a fluctuation amplitude of sigma_8 0.8 rather than the sigma_8=0.9 assumed by the MS.

Building on the legacy of the Sloan Digital Sky Survey (SDSS-I and II), SDSS-III is a program of four spectroscopic surveys on three scientific themes: dark energy and cosmological parameters, the history and structure of the Milky Way, and the population of giant planets around other stars. In keeping with SDSS tradition, SDSS-III will provide regular public releases of all its data, beginning with SDSS DR8 (which occurred in Jan 2011). This paper presents an overview of the four SDSS-III surveys. BOSS will measure redshifts of 1.5 million massive galaxies and Lya forest spectra of 150,000 quasars, using the BAO feature of large scale structure to obtain percent-level determinations of the distance scale and Hubble expansion rate at z<0.7 and at z~2.5. SEGUE-2, which is now completed, measured medium-resolution (R=1800) optical spectra of 118,000 stars in a variety of target categories, probing chemical evolution, stellar kinematics and substructure, and the mass profile of the dark matter halo from the solar neighborhood to distances of 100 kpc. APOGEE will obtain high-resolution (R~30,000), high signal-to-noise (S/N>100 per resolution element), H-band (1.51-1.70 micron) spectra of 10^5 evolved, late-type stars, measuring separate abundances for ~15 elements per star and creating the first high-precision spectroscopic survey of all Galactic stellar populations (bulge, bar, disks, halo) with a uniform set of stellar tracers and spectral diagnostics. MARVELS will monitor radial velocities of more than 8000 FGK stars with the sensitivity and cadence (10-40 m/s, ~24 visits per star) needed to detect giant planets with periods up to two years, providing an unprecedented data set for understanding the formation and dynamical evolution of giant planet systems. (Abridged)

We study relations between stellar mass, star formation and gas-phase metallicity in a sample of 177,071 unique emission line galaxies from the SDSS-DR7, as well as in a sample of 43,767 star forming galaxies at z=0 from the cosmological semi-analytic model L-GALAXIES. We demonstrate that metallicity is dependent on star formation rate at fixed mass, but that the trend is opposite for low and for high mass galaxies. Low-mass galaxies that are actively forming stars are more metal-poor than quiescent low-mass galaxies. High-mass galaxies, on the other hand, have lower gas-phase metallicities if their star formation rates are small. Remarkably, the same trends are found for our sample of model galaxies. We find that massive model galaxies with low gas-phase metallicities have undergone a gas-rich merger in the past, inducing a starburst which exhausted their cold gas reservoirs and shut down star formation. This led to a gradual dilution in the gas-phase metallicities of these systems via accretion of gas. These model galaxies have lower-than-average gas-to-stellar mass ratios and higher-than-average central black hole masses. We confirm that massive galaxies with low gas-phase metallicities in our observational sample also have very massive black holes. We propose that accretion may therefore play a significant role in regulating the gas-phase metallicities of present-day massive galaxies.

We present the results of the analysis of a sample of 17 low-luminosity (L_x < 1e42 erg/s), radio loud AGNs in massive galaxies. The sample is extracted from the SDSS database and it spans uniformly a wide range in optical [OIII] emission line and radio luminosity, but within a narrow redshift range (0.05 < z < 0.11) and a narrow super massive black hole mass range (~ 1e8 M_sun). For these sources we measured core X-ray emission with the Chandra X-ray telescope and radio emission with the VLA. Our main goal is to establish which emission component, if any, can be regarded as the most reliable accretion/jet-power estimator at these regimes. In order to do so, we studied the correlation between emission line properties, radio luminosity, radio spectral slopes and X-ray luminosity, as well as more complex multi-variate relations involving black hole mass, such as the fundamental plane of black hole activity. We find that 15 out of 17 sources of our sample can be classified as Low-Excitation Galaxies (LEG) and their observed properties suggest X-ray and radio emission to originate from the jet basis. We also find that X-ray emission does not appear to be affected by nuclear obscuration and can be used as a reliable jet-power estimator. More generally, X-ray, radio and optical emission appear to be related, although no tight correlation is found. In accordance with a number of recent studies of this class of objects these findings may be explained by a lack of cold (molecular) gaseous structures in the innermost region of these massive galaxies.

We use a stacking technique to measure the average HI content of a volume-limited sample of 1871 AGN host galaxies from a parent sample of galaxies selected from the SDSS and GALEX imaging surveys with stellar masses greater than 10^10 M_sun and redshifts in the range 0.025<z<0.05. HI data are available from the Arecibo Legacy Fast ALFA (ALFALFA) survey. In previous work, we found that the HI gas fraction in galaxies correlates most strongly with the combination of optical/UV colour and stellar surface mass density. We therefore build a control sample of non-AGN matched to the AGN hosts in these two properties. We study trends in HI gas mass fraction (M(HI)/M_*), where M_* is the stellar mass) as a function of black hole accretion rate indicator L[OIII]/M(BH). We find no significant difference in HI content between AGN and control samples at all values of black hole accretion rate probed by the galaxies in our sample. This indicates that AGN do not influence the large-scale gaseous properties of galaxies in the local Universe. We have studied the variation in HI mass fraction with black hole accretion rate in the blue and red galaxy populations. In the blue population, the HI gas fraction is independent of accretion rate, indicating that accretion is not sensitive to the properties of the interstellar medium of the galaxy on large scales. However, in the red population accretion rate and gas fraction do correlate. The measured gas fractions in this population are not too different from the ones expected from a stellar mass loss origin, implying that the fuel supply in the red AGN population could be a mixture of mass loss from stars and gas present in disks.

We study the relation between molecular gas and star formation in a volume-limited sample of 222 galaxies from the COLD GASS survey, with measurements of the CO(1-0) line from the IRAM 30m telescope. The galaxies are at redshifts 0.025<z<0.05 and have stellar masses in the range 10.0<log(M*/Msun)<11.5. The IRAM measurements are complemented by deep Arecibo HI observations and homogeneous SDSS and GALEX photometry. A reference sample that includes both UV and far-IR data is used to calibrate our estimates of star formation rates from the seven optical/UV bands. The mean molecular gas depletion timescale, tdep(H2), for all the galaxies in our sample is 1 Gyr, however tdep(H2) increases by a factor of 6 from a value of ~0.5 Gyr for galaxies with stellar masses of 10^10 Msun to ~3 Gyr for galaxies with masses of a few times 10^11 Msun. In contrast, the atomic gas depletion timescale remains contant at a value of around 3 Gyr. This implies that in high mass galaxies, molecular and atomic gas depletion timescales are comparable, but in low mass galaxies, molecular gas is being consumed much more quickly than atomic gas. The strongest dependences of tdep(H2) are on the stellar mass of the galaxy (parameterized as log tdep(H2)= (0.36+/-0.07)(log M* - 10.70)+(9.03+/-0.99)), and on the specific star formation rate. A single tdep(H2) versus sSFR relation is able to fit both "normal" star-forming galaxies in our COLD GASS sample, as well as more extreme starburst galaxies (LIRGs and ULIRGs), which have tdep(H2) < 10^8 yr. Normal galaxies at z=1-2 are displaced with respect to the local galaxy population in the tdep(H2) versus sSFR plane and have molecular gas depletion times that are a factor of 3-5 times longer at a given value of sSFR due to their significantly larger gas fractions.

We are conducting COLD GASS, a legacy survey for molecular gas in nearby galaxies. Using the IRAM 30m telescope, we measure the CO(1-0) line in a sample of ~350 nearby (D=100-200 Mpc), massive galaxies (log(M*/Msun)>10.0). The sample is selected purely according to stellar mass, and therefore provides an unbiased view of molecular gas in these systems. By combining the IRAM data with SDSS photometry and spectroscopy, GALEX imaging and high-quality Arecibo HI data, we investigate the partition of condensed baryons between stars, atomic gas and molecular gas in 0.1-10L* galaxies. In this paper, we present CO luminosities and molecular hydrogen masses for the first 222 galaxies. The overall CO detection rate is 54%, but our survey also uncovers the existence of sharp thresholds in galaxy structural parameters such as stellar mass surface density and concentration index, below which all galaxies have a measurable cold gas component but above which the detection rate of the CO line drops suddenly. The mean molecular gas fraction MH2/M* of the CO detections is 0.066+/-0.039, and this fraction does not depend on stellar mass, but is a strong function of NUV-r colour. Through stacking, we set a firm upper limit of MH2/M*=0.0016+/-0.0005 for red galaxies with NUV-r>5.0. The average molecular-to-atomic hydrogen ratio in present-day galaxies is 0.3, with significant scatter from one galaxy to the next. The existence of strong detection thresholds in both the HI and CO lines suggests that "quenching" processes have occurred in these systems. Intriguingly, atomic gas strongly dominates in the minority of galaxies with significant cold gas that lie above these thresholds. This suggests that some re-accretion of gas may still be possible following the quenching event.

The source responsible for reionizing the universe at z > 6 remains uncertain. While an energetically adequate population of star-forming galaxies may be in place, it is unknown whether a large enough fraction of their ionizing radiation can escape into the intergalactic medium. Attempts to measure this escape-fraction in intensely star-forming galaxies at lower redshifts have largely yielded upper limits. In this paper we present new HST COS and archival FUSE far-UV spectroscopy of a sample of eleven Lyman Break Analogs (LBAs), a rare population of local galaxies that strongly resemble the high-z Lyman Break galaxies. We combine these data with SDSS optical spectra and Spitzer photometry. We also analyze archival FUSE observations of fifteen typical UV-bright local starbursts. We find evidence of small covering factors for optically-thick neutral gas in 3 cases. This is based on two independent pieces of evidence: a significant residual intensity in the cores of the strongest interstellar absorption-lines tracing neutral gas and a small ratio of extinction-corrected H-alpha to UV plus far-IR luminosities. These objects represent three of the four LBAs that contain a young, very compact (~100pc), and highly massive (~10^9 Mo) dominant central object (DCO). These three objects also differ from the other galaxies in showing a significant amount of blueshifted Ly-alpha emission, which may be related to the low covering factor of neutral gas. All four LBAs with DCOs in our sample show extremely high velocity outflows of interstellar gas, with line centroids blueshifted by about 700km/s and maximum outflow velocities reaching at least 1500km/s. We show that these properties are consistent with an outflow driven by a powerful starburst that is exceptionally compact. We speculate that such extreme feedback may be required to enable the escape of ionizing radiation from star forming galaxies.

We apply updated semi-analytic galaxy formation models simultaneously to the stored halo/subhalo merger trees of the Millennium and Millennium-II simulations. These differ by a factor of 125 in mass resolution, allowing explicit testing of resolution effects on predicted galaxy properties. We have revised the treatments of the transition between the rapid infall and cooling flow regimes of gas accretion, of the sizes of bulges and of gaseous and stellar disks, of supernova feedback, of the transition between central and satellite status as galaxies fall into larger systems, and of gas and star stripping once they become satellites. Plausible values of efficiency and scaling parameters yield an excellent fit not only to the observed abundance of low-redshift galaxies over 5 orders of magnitude in stellar mass and 9 magnitudes in luminosity, but also to the observed abundance of Milky Way satellites. This suggests that reionisation effects may not be needed to solve the "missing satellite" problem except, perhaps, for the faintest objects. The same model matches the observed large-scale clustering of galaxies as a function of stellar mass and colour. The fit remains excellent down to ~30kpc for massive galaxies. For M* < 6 x 10^10Msun, however, the model overpredicts clustering at scales below 1 Mpc, suggesting that the sigma_8 adopted in the simulations (0.9) is too high. Galaxy distributions within rich clusters agree between the simulations and match those observed, but only if galaxies without dark matter subhalos (so-called orphans) are included. Our model predicts a larger passive fraction among low-mass galaxies than is observed, as well as an overabundance of ~10^10Msun galaxies beyond z~0.6, reflecting deficiencies in the way star-formation rates are modelled.

The evolutionary path followed by a galaxy shapes its internal structure, and, in particular, its internal colour variation. We present a study of the internal colour variation within galaxies from the Seventh Data Release of the Sloan Digital Sky Survey (SDSS DR7). We statistically study the connection between the internal colour variation and global galactic properties, looking for hints of the recent galactic evolution. Considering only galaxies with good photometry and spectral measurements, we define four luminosity-threshold samples within the redshift range 0.01<z<0.17, each containing more than 48000 galaxies. Colour gradients are calculated for these galaxies from the surface brightness measurements provided by the SDSS DR7. Possible systematic effects in their determination have been analysed. We find that, on average, galaxies have redder cores than their external parts. We also find that it is more likely to find steep colour gradients among late-type galaxies. This result holds for a range of classifications based on both morphological and spectral characteristics. In fact, our results relate, on average, steep colour gradients to a higher presence of young stars within a galaxy. Our results also suggest that nuclear activity is a marginal driver for creating steep colour gradients in massive galaxies. We have selected pairs of interacting galaxies, with a separation of 5 arcsec, in projected radius, and a difference in redshift of 100 km/s, finding that they present steeper gradients than the average population, skewed towards bluer cores. Our analysis implies that colour gradients can be useful for selecting galaxies that have suffered a recent (minor) burst of star formation.

We analyze a sample of galaxies with stellar masses greater than $10^{10} M_{\odot}$ and with redshifts in the range $0.025<z<0.05$ for which HI mass measurements are available from the GALEX Arecibo SDSS Survey (GASS) or from the Arecibo Legacy Fast ALFA survey (ALFALFA). At a given value of $M_*$, our sample consists primarily of galaxies that are more HI-rich than average. We constructed a series of three control samples for comparison with these HI-rich galaxies. As expected, HI-rich galaxies differ strongly from galaxies of same stellar mass that are selected without regard to HI content. The majority of these differences are attributable to the fact that galaxies with more gas are bluer and more actively star-forming. In order to identify those galaxy properties that are causally connected with HI content, we compare results derived for the HI sample with those derived for galaxies matched in stellar mass, size and NUV-$r$ colour. The only photometric property that is clearly attributable to increasing HI content, is the colour gradient of the galaxy. Galaxies with larger HI fractions have bluer, more actively star-forming outer disks compared to the inner part of the galaxy. HI-rich galaxies also have larger $g$-band radii compared to $i$-band radii. Our results are consistent with the "inside-out" picture of disk galaxy formation, which has commonly served as a basis for semi-analytic models of the formation of disks in the context of Cold Dark Matter cosmologies. The lack of any intrinsic connection between HI fraction and galaxy asymmetry suggests that gas is accreted smoothly onto the outer disk.

We have carried out an HI stacking analysis of a volume-limited sample of ~5000 galaxies with imaging and spectroscopic data from GALEX and the Sloan Digital Sky Survey, which lie within the current footprint of the Arecibo Legacy Fast ALFA (ALFALFA) Survey. Our galaxies are selected to have stellar masses greater than 10^10 Msun and redshifts in the range 0.025<z<0.05. We extract a sub-sample of 1833 "early-type" galaxies with inclinations less than 70deg, with concentration indices C>2.6 and with light profiles that are well fit by a De Vaucouleurs model. We then stack HI line spectra extracted from the ALFALFA data cubes at the 3-D positions of the galaxies from these two samples in bins of stellar mass, stellar mass surface density, central velocity dispersion, and NUV-r colour. We use the stacked spectra to estimate the average HI gas fractions M_HI/M_* of the galaxies in each bin. Our main result is that the HI content of a galaxy is not influenced by its bulge. The average HI gas fractions of galaxies in both our samples correlate most strongly with NUV-r colour and with stellar surface density. The relation between average HI fraction and these two parameters is independent of concentration index C. We have tested whether the average HI gas content of bulge-dominated galaxies on the red sequence, differs from that of late-type galaxies on the red sequence. We find no evidence that galaxies with a significant bulge component are less efficient at turning their available gas reservoirs into stars. This result is in contradiction with the "morphological quenching" scenario proposed by Martig et al. (2009).

We report new observations of the galaxy UGC8802 obtained through GASS, the GALEX Arecibo SDSS Survey, which show this galaxy to be in a remarkable evolutionary state. UGC8802 (GASS35981) is a disk galaxy with stellar mass M*=2x10^10 Msolar which appears to contain an additional 2.1x10^10 Msolar of HI gas. New millimeter observations with the IRAM 30m telescope indicate a molecular gas mass only a tenth this large. Using deep long-slit spectroscopy, we examine the spatially resolved star formation rate and metallicity profiles of GASS35981 for clues to its history. We find that the star formation surface density in this galaxy is low (Sigma_SFR=0.003 Msolar/yr/kpc^2) and that the star formation is spread remarkably evenly across the galaxy. The low molecular gas masses measured in our three IRAM pointings are largely consistent with the total star formation measured within the same apertures. Our MMT long-slit spectrum reveals a sharp drop in metallicity in the outer disk of GASS35981. The ratio of current star formation rate to existing stellar mass surface density in the outer disk is extremely high, implying that all the stars must have formed within the past ~1Gyr. At current star formation rates, however, GASS35981 will not consume its HI reservoir for another 5-7 Gyr. Despite its exceptionally large gas fraction for a galaxy this massive, GASS35981 has a regular rotation curve and exhibits no sign of a recent interaction or merger. We speculate that GASS35981 may have acquired its gas directly from the inter-galactic medium, and that it and other similar galaxies identified in the GASS survey may provide rare local glimpses of gas accretion processes that were more common during the prime epoch of disk galaxy formation at z~1.

We extend existing semi-analytic models of galaxy formation to track atomic and molecular gas in disk galaxies. Simple recipes for processes such as cooling, star formation, supernova feedback, and chemical enrichment of the stars and gas are grafted on to dark matter halo merger trees derived from the Millennium Simulation. Each galactic disk is represented by a series of concentric rings. We assume that surface density profile of infalling gas in a dark matter halo is exponential, with scale radius r_d that is proportional to the virial radius of the halo times its spin parameter $\lambda$. As the dark matter haloes grow through mergers and accretion, disk galaxies assemble from the inside out. We include two simple prescriptions for molecular gas formation processes in our models: one is based on the analytic calculations by Krumholz, McKee & Tumlinson (2008), and the other is a prescription where the H_2 fraction is determined by the kinematic pressure of the ISM. Motivated by the observational results of Leroy et al. (2008), we adopt a star formation law in which $\Sigma_{SFR}\propto\Sigma_{H_2}$ in the regime where the molecular gas dominates the total gas surface density, and $\Sigma_{SFR}\propto \Sigma_{gas}^2$ where atomic hydrogen dominates. We then fit these models to the radial surface density profiles of stars, HI and H_2 drawn from recent high resolution surveys of stars and gas in nearby galaxies. We explore how the ratios of atomic gas, molecular gas and stellar mass vary as a function of global galaxy scale parameters, including stellar mass, stellar surface density, and gas surface density. We elucidate how the trends can be understood in terms of three variables that determine the partition of baryons in disks: the mass of the dark matter halo, the spin parameter of the halo, and the amount of gas recently accreted from the external environment.

We use measurements of the HI content, stellar mass and star formation rates in ~190 massive galaxies with stellar masses greater than 10^10 Msun, obtained from the Galex Arecibo SDSS Survey (GASS) described in Paper I (Catinella et al. 2010) to explore the global scaling relations associated with the bin-averaged ratio of the star formation rate over the HI mass, which we call the HI-based star formation efficiency (SFE). Unlike the mean specific star formation rate, which decreases with stellar mass and stellar mass surface density, the star formation efficiency remains relatively constant across the sample with a value close to SFE = 10^-9.5 yr^-1 (or an equivalent gas consumption timescale of ~3 Gyr). Specifically, we find little variation in SFE with stellar mass, stellar mass surface density, NUV-r color and concentration. We interpret these results as an indication that external processes or feedback mechanisms that control the gas supply are important for regulating star formation in massive galaxies. An investigation into the detailed distribution of SFEs reveals that approximately 5% of the sample shows high efficiencies with SFE > 10^-9 yr^-1, and we suggest that this is very likely due to a deficiency of cold gas rather than an excess star formation rate. Conversely, we also find a similar fraction of galaxies that appear to be gas-rich for their given specific star-formation rate, although these galaxies show both a higher than average gas fraction and lower than average specific star formation rate. Both of these populations are plausible candidates for "transition" galaxies, showing potential for a change (either decrease or increase) in their specific star formation rate in the near future. We also find that 36+/-5% of the total HI mass density and 47+/-5% of the total SFR density is found in galaxies with stellar mass greater than 10^10 Msun. [abridged]

We analyze star forming galaxies drawn from SDSS DR7 to show how the interstellar medium (ISM) Na I 5890, 5896 (Na D) absorption lines depend on galaxy physical properties, and to look for evidence of galactic winds. We combine the spectra of galaxies with similar geometry/physical parameters to create composite spectra with signal-to-noise ~300. The stellar continuum is modeled using stellar population synthesis models, and the continuum-normalized spectrum is fit with two Na I absorption components. We find that: (1) ISM Na D absorption lines with equivalent widths EW > 0.8A are only prevalent in disk galaxies with specific properties -- large extinction (Av), high star formation rates (SFR), high star formation rate per unit area ($\Sigma_{\rm SFR}$), or high stellar mass (M*). (2) the ISM Na D absorption lines can be separated into two components: a quiescent disk-like component at the galaxy systemic velocity and an outflow component; (3) the disk-like component is much stronger in the edge-on systems, and the outflow component covers a wide angle but is stronger within 60deg of the disk rotation axis; (4) the EW and covering factor of the disk component correlate strongly with dust attenuation, highlighting the importance that dust shielding may play the survival of Na I. (5) The EW of the outflow component depends primarily on $\Sigma_{\rm SFR}$ and secondarily on Av; (6) the outflow velocity varies from ~120 to 160km/s but shows little hint of a correlation with galaxy physical properties over the modest dynamic range that our sample probes (1.2 dex in log$\Sigma_{\rm SFR}$ and 1 dex in log M*).

We have compiled a large sample of isolated central galaxies from the Sloan Digital Sky Survey, which do not have a neighbour of comparable brightness within a projected distance of 1 Mpc. We use the colours, luminosities and surface brightnesses of satellite galaxies in the vicinity of these objects to estimate their atomic gas content and to derive the average total mass of HI gas contained in satellites as a function of projected radius from the primary. Recent calibrations of merging timescales from N-body simulations are used to estimate the rate at which this gas will accrete onto the primaries. Our estimated accretion rates fall short of those needed to maintain the observed level of star formation in these systems by nearly two orders of magnitude. Nevertheless, there are strong correlations between the total mass of gas in satellites and the colours and specific star formation rates of central galaxies of all stellar masses. The correlations are much weaker if we consider the total stellar mass in the satellites, rather than their total gas mass. We ask why star formation in the central galaxies should be correlated with gas contained in satellites at projected separations of a Mpc or more, well outside the virial radius of the dark matter halos of these systems. We suggest that gas-rich satellites trace an underlying reservoir of ionized gas that is accreted continuously, and that provides the fuel for ongoing star formation in galaxies in the local Universe.

We investigate how the specific star formation rates of galaxies of different masses depend on cluster-centric radius and on the central/satellite dichotomy in both field and cluster environments. Recent data from a variety of sources, including the cluster catalogue of von der Linden et al. are compared to the semi-analytic models of De Lucia & Blaizot. We find that these models predict too many passive satellite galaxies in clusters, too few passive central galaxies with low stellar masses, and too many passive central galaxies with high masses. We then outline a series of modifications to the model necessary to solve these problems: a) Instead of instantaneous stripping of the external gas reservoir after a galaxy becomes a satellite, the gas supply is assumed to decrease at the same rate that the surrounding halo loses mass due to tidal stripping, b) The AGN feedback efficiency is lowered to bring the fraction of massive passive centrals in better agreement with the data. We also allow for radio mode AGN feedback in satellite galaxies. c) We assume that satellite galaxies residing in host haloes with masses below 10^12 M_sun do not undergo any stripping. We highlight the fact that in low mass galaxies, the external reservoir is composed primarily of gas that has been expelled from the galactic disk by supernovae driven winds. This gas must remain available as a future reservoir for star formation, even in satellite galaxies. Finally, we present a simple recipe for the stripping of gas and dark matter in satellites that can be used in models where subhalo evolution is not followed in detail.

We introduced a semi-analytic model to study the detail of radial related processes on galaxy disks. Based on this recipe, we can calculate neutral and molecular gas in interstellar medium, and include the molecular phase into star formation recipes. From our model, we reproduced the radial surface density profiles and mass functions shown in the observational results. We can also give some predictions and explanations to the scale relations of neutral and molecular gas in disk galaxies.

The GALEX Arecibo SDSS Survey (GASS) is an ambitious program designed to investigate the cold gas properties of massive galaxies, a challenging population for HI studies. Using the Arecibo radio telescope, GASS is gathering high-quality HI-line spectra for an unbiased sample of ~1000 galaxies with stellar masses greater than 10^10 Msun and redshifts 0.025 < z < 0.05, uniformly selected from the SDSS spectroscopic and GALEX imaging surveys. The galaxies are observed until detected or until a low gas mass fraction limit (1.5-5%) is reached. We present initial results based on the first Data Release, which consists of ~20% of the final GASS sample. We use this data set to explore the main scaling relations of HI gas fraction with galaxy structure and NUV-r colour, and show our best fit plane describing the relation between gas fraction, stellar mass surface density and NUV-r colour. Interesting outliers from this plane include gas-rich red sequence galaxies that may be in the process of regrowing their disks, as well as blue, but gas-poor spirals.

We introduce the GALEX Arecibo SDSS Survey (GASS), an on-going large program that is gathering high quality HI-line spectra using the Arecibo radio telescope for an unbiased sample of ~1000 galaxies with stellar masses greater than 10^10 Msun and redshifts 0.025<z<0.05, selected from the SDSS spectroscopic and GALEX imaging surveys. The galaxies are observed until detected or until a low gas mass fraction limit (1.5-5%) is reached. This paper presents the first Data Release, consisting of ~20% of the final GASS sample. We use this data set to explore the main scaling relations of HI gas fraction with galaxy structure and NUV-r colour. A large fraction (~60%) of the galaxies in our sample are detected in HI. We find that the atomic gas fraction decreases strongly with stellar mass, stellar surface mass density and NUV-r colour, but is only weakly correlated with galaxy bulge-to-disk ratio (as measured by the concentration index of the r-band light). We also find that the fraction of galaxies with significant (more than a few percent) HI decreases sharply above a characteristic stellar surface mass density of 10^8.5 Msun kpc^-2. The fraction of gas-rich galaxies decreases much more smoothly with stellar mass. One of the key goals of GASS is to identify and quantify the incidence of galaxies that are transitioning between the blue, star-forming cloud and the red sequence of passively-evolving galaxies. Likely transition candidates can be identified as outliers from the mean scaling relations between gas fraction and other galaxy properties. [abridged]

We compute the cross-correlation between a sample of 14,000 radio-loud AGN (RLAGN) with redshifts between 0.4 and 0.8 selected from the Sloan Digital Sky Survey and a reference sample of 1.2 million luminous red galaxies in the same redshift range. We quantify how the clustering of radio-loud AGN depends on host galaxy mass and on radio luminosity. Radio-loud AGN are clustered more strongly on all scales than control samples of radio-quiet galaxies with the same stellar masses and redshifts, but the differences are largest on scales less than 1 Mpc. In addition, the clustering amplitude of the RLAGN varies significantly with radio luminosity on scales less than 1 Mpc. This proves that the gaseous environment of a galaxy on the scale of its dark matter halo, plays a key role in determining not only the probability that a galaxy is radio-loud AGN, but also the total luminosity of the radio jet. Next, we compare the clustering of radio galaxies with that of radio-loud quasars in the same redshift range. Unified models predict that both types of active nuclei should cluster in the same way. Our data show that most RLAGN are clustered more strongly than radio-loud QSOs, even when the AGN and QSO samples are matched in both black hole mass and radio luminosity. Only the most extreme RLAGN and RLQSOs in our sample, with radio luminosities in excess of 10^26 W/Hz, have similar clustering properties. The majority of the strongly evolving RLAGN population at z~0.5 are found in different environments to the quasars, and hence must be triggered by a different physical mechanism.

We investigate the recent and current star formation activity of galaxies as function of distance from the cluster center in a sample of 521 SDSS clusters at z<0.1. We show that when the BCGs are excluded from the galaxy sample, there is no evidence for mass segregation in the clusters, so that differences in cluster and field populations cannot simply be attributed to different mass functions. We find a marked star formation-radius relation in that almost all galaxies in the cluster core are quiescent, i.e. have terminated star formation a few Gyr ago. This star formation-radius relation is most pronounced for low-mass galaxies and is very weak or absent beyond the virial radius. The typical star formation rate of non-quiescent galaxies declines by approximately a factor of two towards the cluster center. However, the fraction of galaxies with young stellar populations indicating a recently completed starburst or a truncation of star formation does not vary significantly with radius. The fraction of star-forming galaxies which host a powerful optical AGN is also independent of clustercentric radius, indicating that the link between star formation and AGN in these galaxies operates independent of environment. The fraction of red galaxies which host a weak optical AGN decreases, however, towards the cluster center, with a similar timescale as the decline of star forming galaxies. Our results favor a scenario in which star formation is quenched slowly, on timescales similar to the cluster crossing time, and rule out significant contributions from more violent processes at least beyond cluster radii >0.1 R_200.

Many brightest cluster galaxies (BCGs) at the centers of X-ray selected clusters exhibit clear evidence for recent star formation. However, studies of BCGs in optically-selected clusters show that star formation is not enhanced when compared to control samples of non-BCGs of similar stellar mass. Here we analyze a sample of 113 BCGs in low redshift (z<0.1), optically-selected clusters, a matched control sample of non-BCGs, and a smaller sample of BCGs in X-ray selected clusters. We convolve the SDSS images of the BCGs to match the resolution of the GALEX data and we measure UV-optical colours in their inner and outer regions. We find that optically-selected BCGs exhibit smaller scatter in optical colours and redder inner NUV-r colours than the control galaxies, indicating that they are a homogenous population with very little ongoing star formation. The BCGs in the X-ray selected cluster sample span a similar range in optical colours, but have bluer NUV-r colours. Among X-ray selected BCGs, those located in clusters with central cooling times of less than 1 Gyr are significantly bluer than those located in clusters where the central gas cooling times are long. Our main conclusion is that the location of a galaxy at the centre of its halo is not sufficient to determine whether or not it is currently forming stars. One must also have information about the thermodynamic state of the gas in the core of the halo.

Using results from structural analysis of a sample of nearly 1000 local galaxies from the Sloan Digital Sky Survey, we estimate how the mass in central black holes is distributed amongst elliptical galaxies, classical bulges and pseudo-bulges, and investigate the relation between their stellar masses and central stellar velocity dispersion sigma. Assuming a single relation between elliptical galaxy/bulge mass, M_Bulge, and central black hole mass, M_BH, we find that 55^{+8}_{-4} per cent of the mass in black holes in the local universe is in the centres of elliptical galaxies, 41^{+4}_{-2} per cent in classical bulges and 4^{+0.9}_{-0.4} per cent in pseudo-bulges. We find that ellipticals, classical bulges and pseudo-bulges follow different relations between their stellar masses and sigma, and the most significant offset occurs for pseudo-bulges in barred galaxies. This structural dissimilarity leads to discrepant black hole masses if single M_BH-M_Bulge and M_BH-sigma relations are used. Adopting relations from the literature, we find that the M_BH-sigma relation yields an estimate of the total mass density in black holes that is roughly 55 per cent larger than if the M_BH-M_Bulge relation is used.

We use a sample of 800 galaxies with H I mass measurements from the HyperLeda catalogue and optical photometry from the fourth data release of the Sloan Digital Sky Survey to calibrate a new photometric estimator of the H I to-stellar mass ratio for nearby galaxies. Our estimator, which is motivated by the Kennicutt-Schmidt star formation law, is log(G_{HI}/S) = -1.73238(g-r) + 0.215182mu_i - 4.08451, where mu_i is the i-band surface brighteness and g-r is the optical colour estimated from the g- and r-band Petrosian apparent agnitudes. This estimator has a scatter of sigma = 0.31 dex in log(G_{HI}/S), compared to sigma ~ 0.4 dex for previous estimators that were based on colour alone. We investigate whether the residuals in our estimate of log(G_{HI}/S) depend in a systematic way on a variety of different galaxy properties. We find no effect as a function of stellar mass or 4000A break strength, but there is a systematic effect as a function of the concentration index of the light. We then apply our estimator to a sample of 10^5 emission-line galaxies in the SDSS DR4 and derive an estimate of the H I mass function, which is in excellent agreement with recent results from H I blind surveys. Finally, we re-examine the well-known relation between gas-phase metallicity and stellar mass and ask whether there is a dependence on H I-to-stellar mass ratio, as predicted by chemical evolution models. We do find that gas-poor galaxies are more metal rich at fixed stellar mass. We compare our results with the semi-analytic models of De Lucia & Blaizot, which include supernova feedback, as well as the cosmological infall of gas.

We study the clustering properties of barred galaxies using data from the Sloan Digital Sky Survey (SDSS). We compute projected redshift-space two-point cross-correlation functions w_p(r_p) for a sample of nearly 1000 galaxies for which we have performed detailed structural decompositions using the methods described in Gadotti (2009). The sample includes 286 barred galaxies. The clustering of barred and unbarred galaxies of similar stellar mass is indistinguishable over all the scales probed (~20kpc - 30Mpc). This result also holds even if the sample is restricted to bars with bluer g-i colours (and hence younger ages). Our result also does not change if we split our sample of barred galaxies according to bar-to-total luminosity ratio, bar boxyness, effective surface brightness, length, or the shape of the surface density profile within the bar. There is a hint that red, elliptical bars are more strongly clustered than red and less elliptical bars, on scales >~ 1 Mpc, although the statistical significance is not high. We conclude that there is no significant evidence that bars are a product of mergers or interactions. We tentatively interpret the stronger clustering of the more elliptical bars as evidence that they are located in older galaxies, which reside in more massive haloes.

We analyze the observed distribution of Eddington ratios as a function of supermassive black hole mass for a large sample of nearby galaxies drawn from the Sloan Digital Sky Survey. We demonstrate that there are two distinct regimes of black hole growth in nearby galaxies. The first is associated with galaxies with significant star formation in their central kiloparsec regions, and is characterized by a broad log-normal distribution of accretion rates peaked at about one percent of the Eddington limit. In this regime, the Eddington ratio distribution is independent of the mass of the black hole and shows no further dependence on the central stellar population of the galaxy. The second regime is associated with galaxies with old central stellar populations, and is characterized by a power-law distribution function of Eddington ratios. In this regime, the time-averaged mass accretion rate onto black holes is proportional to the mass of stars in the galaxy bulge, with a constant of proportionality that depends on the mean stellar age of the stars. This result is once again independent of black hole mass. We show that both the slope of the power-law and the decrease in the accretion rate onto black holes in old galaxies are consistent with population synthesis model predictions of the decline in stellar mass loss rates as a function of mean stellar age. Our results lead to a very simple picture of black hole growth in the local Universe. If the supply of cold gas in a galaxy bulge is plentiful, the black hole regulates its own growth at a rate that does not further depend on the properties of the interstellar medium. Once the gas runs out, black hole growth is regulated by the rate at which evolved stars lose their mass.

We present a detailed study of the Galaxy Evolution Explorer's photometric catalogs with special focus on the statistical properties of the All-sky and Medium Imaging Surveys. We introduce the concept of primaries to resolve the issue of multiple detections and follow a geometric approach to define clean catalogs with well-understood selection functions. We cross-identify the GALEX sources (GR2+3) with Sloan Digital Sky Survey (DR6) observations, which indirectly provides an invaluable insight about the astrometric model of the UV sources and allows us to revise the band merging strategy. We derive the formal description of the GALEX footprints as well as their intersections with the SDSS coverage along with analytic calculations of their areal coverage. The crossmatch catalogs are made available for the public. We conclude by illustrating the implementation of typical selection criteria in SQL for catalog subsets geared toward statistical analyses, e.g., correlation and luminosity function studies.

Over the last decade optical spectroscopic surveys have characterized the low redshift galaxy population and uncovered populations of star-forming galaxies back to z ~ 7. This work has shown that the primary epoch of galaxy building and black hole growth occurs at redshifts of 2 to 3. The establishment of the concordance LCDM cosmology shifted the focus of galaxy population studies from constraining cosmological parameters to characterizing the processes which regulate the formation and evolution of galaxies.In the next decade, high redshift observers will attempt to formulate a coherent evolutionary picture connecting galaxies in the high redshift Universe to galaxies today. In order to link galaxy populations at different redshifts, we must not only characterize their evolution in a systematic way, we must establish which physical processes are responsible for it. Considerable progress has already been made in understanding how galaxies evolved from z ~ 1 to the present day. Large spectroscopic surveys in the near infrared are required to push these studies back towards the main epoch of galaxy building. Only then will we understand the full story of the formation of L* galaxies like our own Milky Way. A large near-IR spectroscopic survey will also provide the calibration needed to avoid systematics in the large photometric programs proposed to study the nature of dark matter and dark energy. We provide an outline design for a multi-object 0.4 to 1.8 micron spectrograph, which could be placed on an existing telescope, and which would allow a full characterization of the galaxy population out to z ~ 2. We strongly recommend a serious further study to design a real instrument, which will be required for galaxy formation studies to advance to the next frontier.

One of the main themes in extragalactic astronomy for the next decade will be the evolution of galaxies over cosmic time. Many future observatories, including JWST, ALMA, GMT, TMT and E-ELT will intensively observe starlight over a broad redshift range, out to the dawn of the modern Universe when the first galaxies formed. It has, however, become clear that the properties and evolution of galaxies are intimately linked to the growth of their central black holes. Understanding the formation of galaxies, and their subsequent evolution, will therefore be incomplete without similarly intensive observations of the accretion light from supermassive black holes (SMBH) in galactic nuclei. To make further progress, we need to chart the formation of typical SMBH at z>6, and their subsequent growth over cosmic time, which is most effectively achieved with X-ray observations. Recent technological developments in X-ray optics and instrumentation now bring this within our grasp, enabling capabilities fully matched to those expected from flagship observatories at longer wavelengths.

We analyse a sample of 52,000 Milky Way (MW) type galaxies drawn from the publicly available galaxy catalogue of the Millennium Simulation with the aim of studying statistically the differences and similarities of their properties in comparison to our Galaxy. Model galaxies are chosen to lie in haloes with maximum circular velocities in the range 200-250 km/seg and to have bulge-to-disk ratios similar to that of the Milky Way. We find that model MW galaxies formed quietly through the accretion of cold gas and small satellite systems. Only 12 per cent of our model galaxies experienced a major merger during their lifetime. Most of the stars formed in situ, with only about 15 per cent of the final mass gathered through accretion. Supernovae and AGN feedback play an important role in the evolution of these systems. At high redshifts, when the potential wells of the MW progenitors are shallower, winds driven by supernovae explosions blow out a large fraction of the gas and metals. As the systems grow in mass, SN feedback effects decrease and AGN feedback takes over, playing a more important role in the regulation of the star formation activity at lower redshifts. Although model Milky Way galaxies have been selected to lie in a narrow range of maximum circular velocities, they nevertheless exhibit a significant dispersion in the final stellar masses and metallicities. Our analysis suggests that this dispersion results from the different accretion histories of the parent dark matter haloes. Statically, we also find evidences to support the Milky Way as a typical Sb/Sc galaxy in the same mass range, providing a suitable benchmark to constrain numerical models of galaxy formation

Using the SDSS DR4 group catalogue of Yang et al. (2007), we investigate sizes, concentrations, colour gradients and surface brightness profiles of central and satellite galaxies. We compare central and satellite galaxies at fixed stellar mass, in order to disentangle environmental from stellar mass dependencies. Early and late type galaxies are defined according to concentration. We find that at fixed stellar mass, late type satellite galaxies have smaller radii and larger concentrations than late type central galaxies. No such differences are found for early-type galaxies. We have also constructed surface brightness and colour profiles for the central and satellite galaxies in our sample. We find that late-type satellite galaxies have a lower surface brightness and redder colours than late-type central galaxies. We show that all observed differences between satellite and central galaxies can be explained by a simple fading model, in which the star formation in the disk decreases over timescales of 2-3 Gyr after a galaxy becomes a satellite. Processes that induce strong morphological changes (e.g. harassment) and processes that strip the galaxy of its entire ISM need not to be invoked in order to explain the environmental dependencies we find.

Motivated by recent observational studies of the environment of z~6 QSOs, we have used the Millennium Run (MR) simulations to construct a very large (~20 deg^2) mock redshift survey of star-forming galaxies at z~6. We use this simulated survey to study the relation between density enhancements in the distribution of i-dropouts and Lya-emitters, and their relation to the most massive halos and protocluster regions at z~6. Our simulation predicts significant variations in surface density across the sky with some voids and filaments extending over scales of 1 degree, much larger than probed by current surveys. Approximately one third of all z~6 halos hosting i-dropouts brighter than z=26.5 mag (~M* at z=6) become part of z=0 galaxy clusters. i-dropouts associated with protocluster regions are found in regions where the surface density is enhanced on scales ranging from a few to several tens of arcminutes on the sky. We analyze two structures of i-dropouts and Lya-emitters observed with the Subaru Telescope and show that these structures must be the seeds of massive clusters-in-formation. In striking contrast, six z~6 QSO fields observed with HST show no significant enhancements in their i-dropout number counts. With the present data, we cannot rule out the QSOs being hosted by the most massive halos. However, neither can we confirm this widely used assumption. We conclude by giving detailed recommendations for the interpretation and planning of observations by current and future ground- and space based instruments that will shed new light on questions related to the large-scale structure at z~6.

We present a new method to estimate the average star formation rate per unit stellar mass (SSFR) of a stacked population of galaxies. We combine the spectra of 600-1000 galaxies with similar stellar masses and parameterise the star formation history of this stacked population using a set of exponentially declining functions. The strength of the Hydrogen Balmer absorption line series in the rest-frame wavelength range 3750-4150\AA is used to constrain the SSFR by comparing with a library of models generated using the BC03 stellar population code. Our method, based on a principal component analysis (PCA), can be applied in a consistent way to spectra drawn from local galaxy surveys and from surveys at $z \sim 1$, and is only weakly influenced by attenuation due to dust. We apply our method to galaxy samples drawn from SDSS and DEEP2 to study mass-dependent growth of galaxies from $z \sim 1$ to $z \sim 0$. We find that, (1) high mass galaxies have lower SSFRs than low mass galaxies; (2) the average SSFR has decreased from $z=1$ to $z=0$ by a factor of $\sim 3-4$, independent of galaxy mass. Additionally, at $z \sim 1$ our average SSFRs are a factor of $2-2.5$ lower than those derived from multi-wavelength photometry using similar datasets. We then compute the average time (in units of the Hubble time, $t_{\rm H}(z)$) needed by galaxies of a given mass to form their stars at their current rate. At both $z=0$ and at $z=1$, this timescale decreases strongly with stellar mass from values close to unity for galaxies with masses $\sim 10^{10} M_{\odot}$, to more than ten for galaxies more massive than $ 10^{11} M_{\odot}$. Our results are in good agreement with models in which AGN feedback is more efficient at preventing gas from cooling and forming stars in high mass galaxies.

We have used the Sloan Digital Sky Survey (SDSS) to undertake an investigation of lopsidedness in a sample of ~25,000 nearby galaxies (z < 0.06). We use the m=1 azimuthal Fourier mode between the 50% and 90% light radii as our measure of lopsidedness. The SDSS spectra are used to measure the properties of the stars, gas, and black hole in the central-most few-kpc-scale region. We show that there is a strong link between lopsidedness in the outer parts of the galactic disk and the youth of the stellar population in the central region. This link is independent of the other structural properties of the galaxy. These results provide a robust statistical characterization of the connections between accretion/interactions/mergers and the resulting star formation. We also show that residuals in the galaxy mass-metallicity relation correlate with lopsidedness (at fixed mass, the more metal-poor galaxies are more lopsided). This suggests that the events causing lopsidedness and enhanced star formation deliver lower metallicity gas into the galaxy's central region. Finally, we find that there is a trend for the more powerful active galactic nuclei to be hosted by more lopsided galaxies (at fixed galaxy mass, density, or concentration). However if we compare samples matched to have both the same structures and central stellar populations, we then find no difference in lopsidedness between active and non-active galaxies. This leads to the following picture. The presence of cold gas in the central region of a galaxy (irrespective of its origin) is essential for both star-formation and black hole growth. The delivery of cold gas is aided by the processes that produce lopsidedness. Other processes on scales smaller than we can probe with our data are required to transport the gas to the black hole.

We compute two-point correlation functions and measure the shear signal due to galaxy-galaxy lensing for 80,000 optically identified and 5,700 radio-loud AGN from Data Release 4 (DR4) of the Sloan Digital Sky Survey. Halo occupation models are used to estimate halo masses and satellite fractions for these two types of AGN. The large sample size allows us to separate AGN according to the stellar mass of their host galaxies. We study how the halo masses of optical and radio AGN differ from those of the parent population at fixed M*. Halo masses deduced from clustering and from lensing agree satisfactorily. Radio AGN are found in more massive halos than optical AGN: in our samples their mean halo masses are 1.6 x 10^{13} and 8 x 10^{11} M_{sun}/h, respectively. Optical AGN follow the same relation between stellar mass and halo mass as galaxies selected without regard to nuclear properties, but radio-loud AGN deviate significantly from this relation. The dark matter halos of radio-loud AGN are about twice as massive as those of control galaxies of the same stellar mass. This boost is independent of radio luminosity, and persists even when our analysis is restricted to field galaxies. The large-scale gaseous environment of the galaxy clearly plays a crucial role in producing observable radio emission. The dark matter halo masses that we derive for the AGN in our two samples are in good agreement with recent models in which feedback from radio AGN becomes dominant in halos where gas cools quasi-statically.

A catalogue of 14453 radio-loud AGN with 1.4 GHz fluxes above 3.5 mJy in the redshift range 0.4<z<0.8, has been constructed from the cross-correlation of the NVSS and FIRST radio surveys with the MegaZ-LRG catalogue of luminous red galaxies derived from Sloan Digital Sky Survey imaging data. New techniques were developed for extending the cross-correlation algorithm to FIRST detections that are below the nominal 1 mJy S/N limit of the catalogued sources. We estimate a reliability of ~98.3%, and completeness level (for LRGS) of about 95% for our new catalogue. We present a new determination of the luminosity function of radio AGN at z~0.55 and compare this to the luminosity function of nearby (z~0.1) radio sources from the SDSS main survey. The comoving number density of radio AGN with luminosities less than 10^{25} W Hz^{-1} increases by a factor ~1.5 between z=0.1 and z=0.55. At higher lumiosities, this factor increases sharply, reaching values of more than 10 at radio luminosities larger than 10^{26} W Hz^{-1}. We then study how the relation between radio AGN and their host galaxies evolves with redshift. Our main conclusion is that the fraction of radio-loud AGN increases towards higher redshift in all massive galaxies, but the evolution is particularly strong for the lower mass galaxies in our sample. These trends may be understood if there are two classes of radio galaxies (likely associated with the "radio" and "quasar mode" dichotomy) that have different fuelling/triggering mechanisms and hence evolve in different ways.

We use the ROSAT All-Sky Survey to study the X-ray properties of a sample of 625 groups and clusters of galaxies selected from the Sloan Digital Sky Survey. We stack clusters with similar velocity dispersions and investigate whether their average X-ray luminosities and surface brightness profiles vary with the radio activity level of their central galaxies. We find that at a given value of $\sigma$, clusters with a central radio AGN have more concentrated X-ray surface brightness profiles, larger central galaxy masses, and higher X-ray luminosities than clusters with radio-quiet central galaxies. The enhancement in X-ray luminosity is more than a factor of two, is detected with better than 6$\sigma$ significance, and cannot be explained by X-ray emission from the radio AGN itself. This difference is largely due to a subpopulation of radio-quiet, high velocity dispersion clusters with low mass central galaxies. These clusters are underluminous at X-ray wavelengths when compared to otherwise similar clusters where the central galaxy is radio-loud, more massive, or both.

We describe the scientific motivations, the mission concept and the instrumentation of SPACE, a class-M mission proposed for concept study at the first call of the ESA Cosmic-Vision 2015-2025 planning cycle. SPACE aims to produce the largest three-dimensional evolutionary map of the Universe over the past 10 billion years by taking near-IR spectra and measuring redshifts for more than half a billion galaxies at 0<z<2 down to AB~23 over 3\pi sr of the sky. In addition, SPACE will also target a smaller sky field, performing a deep spectroscopic survey of millions of galaxies to AB~26 and at 2<z<10+. These goals are unreachable with ground-based observations due to the ~500 times higher sky background. To achieve the main science objectives, SPACE will use a 1.5m diameter Ritchey-Chretien telescope equipped with a set of arrays of Digital Micro-mirror Devices (DMDs) covering a total field of view of 0.4 deg2, and will perform large-multiplexing multi-object spectroscopy (e.g. ~6000 targets per pointing) at a spectral resolution of R~400 as well as diffraction-limited imaging with continuous coverage from 0.8mum to 1.8mum.

The GALEX Arecibo SDSS Survey (GASS) is a large targeted survey that started at Arecibo in March 2008. GASS is designed to measure the neutral hydrogen content of ~1000 massive galaxies (with stellar mass Mstar > 10^10 Msun) at redshift 0.025<z<0.05, uniformly selected from the SDSS spectroscopic and GALEX imaging surveys. Our selected mass range straddles the recently identified "transition mass" (Mstar ~3x10^10 Msun) above which galaxies show a marked decrease in their present to past-averaged star formation rates. GASS will produce the first statistically significant sample of massive "transition" galaxies with homogeneously measured stellar masses, star formation rates and gas properties. The analysis of this sample will allow us to investigate if and how the cold gas responds to a variety of different physical conditions in the galaxy, thus yielding insights on the physical processes responsible for the transition between blue, star-forming and red, passively evolving galaxies. GASS will be of considerably legacy value not only in isolation but also by complementing ongoing HI-selected surveys.

Using data from the Sloan Digital Sky Survey data release 3 (SDSS DR3) we investigate how narrow (<700km/s) CIV and MgII quasar absorption line systems are distributed around quasars. The CIV absorbers lie in the redshift range 1.6 < z < 4 and the MgII absorbers in the range 0.4<z<2.2. By correlating absorbers with quasars on different but neighbouring lines-of-sight, we measure the clustering of absorbers around quasars on comoving scales between 4 and 30Mpc. The observed comoving correlation lengths are r_o~5h^-1Mpc, similar to those observed for bright galaxies at these redshifts. Comparing with correlations between absorbers and the quasars in whose spectra they are identified then implies: (i) that quasars destroy absorbers to comoving distances of ~300kpc (CIV) and ~800kpc (MgII) along their lines-of-sight; (ii) that >40% of CIV absorbers within 3,000km/s of the QSO are not a result of large-scale clustering but rather are directly associated with the quasar itself; (iii) that this intrinsic absorber population extends to outflow velocities of order 12,000km/s; (iv) that this outflow component is present in both radio-loud and radio-quiet quasars; and (v) that a small high-velocity outflow component is observed in the MgII population, but any further intrinsic absorber component is undetectable in our clustering analysis. We also find an indication that absorption systems within 3,000km/s are more abundant for radio-loud than for radio-quiet quasars. This suggests either that radio-loud objects live in more massive halos, or that their radio activity generates an additional low-velocity outflow, or both.

It has long been known that galaxy interactions are associated with enhanced star formation. In a companion paper, we explored this connection by applying a variety of statistics to SDSS data. In particular, we showed that specific star formation rates of galaxies are higher if they have close neighbours. Here we apply exactly the same techniques to AGN in the survey, showing that close neighbours are not associated with any similar enhancement of nuclear activity. Star formation is enhanced in AGN with close neighbours in exactly the same way as in inactive galaxies, but the accretion rate onto the black hole, as estimated from the extinction-corrected [O III] luminosity, is not influenced by the presence or absence of companions. Previous work has shown that galaxies with more strongly accreting black holes contain more young stars in their inner regions. This leads us to conclude that star formation induced by a close companion and star formation associated with black hole accretion are distinct events. These events may be part of the same physical process, for example a merger, provided they are separated in time. In this case, accretion onto the black hole and its associated star formation would occur only after the two interacting galaxies have merged. The major caveat in this work is our assumption that the extinction-corrected [O III] luminosity is a robust indicator of the bolometric luminosity of the central black hole. It is thus important to check our results using indicators of AGN activity at other wavelengths.

We use semi-analytic models implemented in the Millennium Simulation to analyze the merging histories of dark matter haloes and of the galaxies that reside in them. We assume that supermassive black holes only exist in galaxies that have experienced at least one major merger. Only a few percent of galaxies with stellar masses less than $M_* < 10^{10} M_{\odot}$ are predicted to have experienced a major merger and to contain a black hole. The fraction of galaxies with black holes increases very steeply at larger stellar masses. This agrees well with the observed strong mass dependence of the fraction of nearby galaxies that contain either low-luminosity (LINER-type) or higher-luminosity (Seyfert or composite-type) AGN. We then investigate when the major mergers that first create the black holes are predicted to occur. High mass galaxies are predicted to have formed their black holes at very early epochs. The majority of low mass galaxies never experience a major merger and hence do not contain a black hole, but a significant fraction of the supermassive black holes that do exist in low mass galaxies are predicted to have formed recently.

We investigate the clustering properties of a complete sample of 10^5 star-forming galaxies drawn from the SDSS DR4. On scales less than 100 kpc, the amplitude of the correlation function exhibits a strong dependence on the specific star formation rate of the galaxy. We interpret this as the signature of enhanced star formation induced by tidal interactions. We then explore how the average star formation rate in a galaxy is enhanced as the projected separation r_p between the galaxy and its companions decreases. We find that the enhancement depends strongly on r_p, but very weakly on the relative luminosity of the companions. The enhancement is also stronger in low mass galaxies than in high mass galaxies. In order to explore whether a tidal interaction is not only sufficient, but also necessary to trigger enhanced star formation in a galaxy, we compute background subtracted neighbour counts for the galaxies in our sample. The average number of close neighbours around galaxies with low to average values of SFR/M* is close to zero. At the highest specific star formation rates, however, more than 40% of the galaxies in our sample have a companion within a projected radius of 100 kpc. Visual inspection of the highest SFR/M* galaxies without companions reveals that more than 50% of these are clear interacting or merging systems. We conclude that tidal interactions are the dominant trigger of enhanced star formation in the most strongly star-forming systems. Finally, we find clear evidence that tidal interactions not only lead to enhanced star formation in galaxies, but also cause structural changes such as an increase in concentration.

Heckman et al. (2005) used the Galaxy Evolution Explorer (GALEX) UV imaging survey to show that there exists a rare population of nearby compact UV-luminous galaxies (UVLGs) that closely resembles high redshift Lyman break galaxies (LBGs). We present HST images in the UV, optical, and Ha, and resimulate them at the depth and resolution of the GOODS/UDF fields to show that the morphologies of UVLGs are also similar to those of LBGs. Our sample of 8 LBG analogs thus provides detailed insight into the connection between star formation and LBG morphology. Faint tidal features or companions can be seen in all of the rest-frame optical images, suggesting that the starbursts are the result of a merger or interaction. The UV/optical light is dominated by unresolved (~100-300 pc) super starburst regions (SSBs). A detailed comparison with the galaxies Haro 11 and VV 114 at z=0.02 indicates that the SSBs themselves consist of diffuse stars and (super) star clusters. The structural features revealed by the new HST images occur on very small physical scales and are thus not detectable in images of high redshift LBGs, except in a few cases where they are magnified by gravitational lensing. We propose, therefore, that LBGs are mergers of gas-rich, relatively low-mass (~10^10 Msun) systems, and that the mergers trigger the formation of SSBs. If galaxies at high redshifts are dominated by SSBs, then the faint end slope of the luminosity function is predicted to have slope alpha~2. Our results are the most direct confirmation to date of models that predict that the main mode of star formation in the early universe was highly collisional.

We report on the properties of a sample of ultraviolet luminous galaxies (UVLGs) selected by matching the Galaxy Evolution Explorer (GALEX) Surveys with the Sloan Digital Sky Survey Third Data Release. Out of 25362 galaxies between 0.0<z<0.3 detected by GALEX, there are 215 galaxies with L>2x10^10 L_solar at 1530 Angstroms (observed wavelength). The properties of this population are well correlated with ultraviolet surface brightness. We find that the galaxies with low UV surface brightness are primarily large spiral systems with a mixture of old and young stellar populations, while the high surface brightness galaxies consist primarily of compact starburst systems. In terms of the behavior of surface brightness with luminosity, size with luminosity, the mass-metallicity relation, and other parameters, the compact UVLGs clearly depart from the trends established by the full sample of galaxies. The subset of compact UVLGs with the highest surface brightness (``supercompact UVLGs'') have characteristics that are remarkably similar to Lyman Break Galaxies at higher redshift. They are much more luminous than typical local ultraviolet-bright starburst galaxies and blue compact dwarf galaxies. They have metallicities that are systematically lower than normal galaxies of the same stellar mass, indicating that they are less chemically evolved. In all these respects, they are the best local analogs for Lyman Break Galaxies.

Large-scale asymmetries in the stellar mass distribution in galaxies are believed to trace non-equilibrium situations in the luminous and/or dark matter component. These may arise in the aftermath of events like mergers, accretion, and tidal interactions. These events are key in the evolution of galaxies. In this paper we quantify the large-scale lopsidedness of light distributions in 25155 galaxies at z < 0.06 from the Sloan Digital Sky Survey Data Release 4 using the m = 1 azimuthal Fourier mode. We show that the lopsided distribution of light is primarily due to a corresponding lopsidedness in the stellar mass distribution. Observational effects, such as seeing, Poisson noise, and inclination, introduce only small errors in lopsidedness for the majority of this sample. We find that lopsidedness correlates strongly with other basic galaxy structural parameters: galaxies with low concentration, stellar mass, and stellar surface mass density tend to be lopsided, while galaxies with high concentration, mass, and density are not. We find that the strongest and most fundamental relationship between lopsidedness and the other structural parameters is with the surface mass density. We also find, in agreement with previous studies, that lopsidedness tends to increase with radius. Both these results may be understood as a consequence of several factors. The outer regions of galaxies and low-density galaxies are more susceptible to tidal perturbations, and they also have longer dynamical times (so lopsidedness will last longer). They are also more likely to be affected by any underlying asymmetries in the dark matter halo.

The majority of nearby radio-loud AGN are found in massive, old elliptical galaxies with weak emission lines. At high redshifts,however, most known radio AGN have strong emission lines. In this paper, we examine a subset of radio AGN with emission lines selected from the Sloan Digital Sky Survey. The probability for a nearby radio AGN to have emission lines is a strongly decreasing function of galaxy mass and an increasing function of radio luminosity above 10^25 W/Hz. Emission line and radio luminosities are correlated, but with large dispersion. At a given radio power, AGN with small black holes have higher [OIII] luminosities (which we interpret as higher accretion rates) than AGN with big black holes. However, if we scale both radio and emission line luminosities by the black hole mass, we find a correlation between normalized radio power and accretion rate in Eddington units that is independent of black hole mass. There is also a clear correlation between normalized radio power and the age of the stellar population in the galaxy. Present-day AGN with the highest normalized radio powers are confined to galaxies with small black holes. High-redshift, high radio-luminosity AGN could be explained if big black holes were similarly active at earlier cosmic epochs. To investigate why only a small fraction of emission line AGN become radio loud, we create matched samples of radio-loud and radio-quiet AGN and compare their host galaxy properties and environments. The main difference lies in their environments; our local density estimates are a factor 2 larger around the radio-loud AGN. We propose a scenario in which radio-loud AGN with emission lines are located in galaxies where accretion of both cold and hot gas can occur simultaneously. (Abridged)

We use empirical techniques to interpret the near-infrared colours of a sample of 5800 galaxies drawn from Sloan Digital Sky Survey (SDSS) main spectroscopic sample with YJHK photometry from the UK Infrared Deep Sky Survey (UKIDSS) data release one. We study correlations between near-IR colours measured within SDSS fibre and physical parameters derived from the spectra. These parameters include specific star formation rate, stellar age, metallicity and dust attenuation. All correlations are analyzed for samples of galaxies that are closely matched in redshift, in stellar mass and in concentration index. Whereas more strongly star-forming galaxies have bluer optical colours, the opposite is true at near-IR wavelengths -- galaxies with higher specific star formation rate have redder near-IR colours. This result agrees qualitatively with the predictions of models in which Thermally Pulsing Asymptotic Giant Branch (TP-AGB) stars dominate the H and K-band light of a galaxy following a burst of star formation. We also find a surprisingly strong correlation between the near-IR colours of star-forming galaxies and their dust attenuation as measured from the Balmer decrement. Unlike optical colours, however, near-IR colours exhibit very little dependence on galaxy inclination. This suggests that the correlation of near-IR colours with dust attenuation arises because TP-AGB stars are the main source of dust in the galaxy. Finally, we compare the near-IR colours of the galaxies in our sample to the predictions of three different stellar population models: the Bruzual & Charlot 2003 model, a preliminary version of a new model under development by Charlot & Bruzual, which includes a new prescription for AGB star evolution, and the Maraston 2005 model.

[Abridged] We investigate trends between the recent star formation history and black hole growth in galaxy bulges in the Sloan Digital Sky Survey (SDSS). The galaxies lie at 0.01<z<0.07 where the fibre aperture covers only the central 0.6-4.0kpc diameter of the galaxy. We find strong trends between black hole growth, as measured by dust-attenuation-corrected OIII luminosity, and the recent star formation history of the bulges. We conclude that our results support the popular hypothesis for black hole growth occurring through gas inflow into the central regions of galaxies, followed by a starburst and triggering of the AGN. However, while this is a significant pathway for the growth of black holes, it is not the dominant one in the present-day Universe. More unspectacular processes are apparently responsible for the majority of this growth. In order to arrive at these conclusions we have developed a set of new high signal-to-noise ratio (SNR) optical spectral indicators, designed to allow a detailed study of stellar populations which have undergone recent enhanced star formation. Working in the rest-frame wavelength range 3750-4150AA, ideally suited to many recent and ongoing spectroscopic surveys at low and high redshift, the first two indices are equivalent to the previously well studied 4000AA break strength and Hdelta equivalent width. The primary advantage of this new method is a greatly improved SNR for the latter index, allowing the present study to use spectra with SNR-per-pixel as low as 8.

We use the Sloan Digital Sky Survey to construct a sample of 625 brightest group and cluster galaxies (BCGs) together with control samples of non-BCGs matched in stellar mass, redshift, and color. We investigate how the systematic properties of BCGs depend on stellar mass and on their privileged location near the cluster center. The groups and clusters that we study are drawn from the C4 catalogue of Miller et al. (2005) but we have developed improved algorithms for identifying the BCG and for measuring the cluster velocity dispersion. Since the SDSS photometric pipeline tends to underestimate the luminosities of large galaxies in dense environments, we have developed a correction for this effect which can be readily applied to the published catalog data. We find that BCGs are larger and have higher velocity dispersions than non-BCGs of the same stellar mass, which implies that BCGs contain a larger fraction of dark matter. In contrast to non-BCGs, the dynamical mass-to-light ratio of BCGs does not vary as a function of galaxy luminosity. Hence BCGs lie on a different fundamental plane than ordinary elliptical galaxies. BCGs also follow a steeper Faber-Jackson relation than non-BCGs, as suggested by models in which BCGs assemble via dissipationless mergers along preferentially radial orbits. We find tentative evidence that this steepening is stronger in more massive clusters. BCGs have similar mean stellar ages and metallicities to non-BCGs of the same mass, but they have somewhat higher alpha/Fe ratios, indicating that star formation may have occurred over a shorter timescale in the BCGs. Finally, we find that BCGs are more likely to host radio-loud active galactic nuclei than other galaxies of the same mass, but are less likely to host an optical AGN. The differences we find are more pronounced for the less massive BCGs.

(Abridged) The prevalence of radio-loud AGN activity in present-day massive halos is determined using a sample of 625 nearby groups and clusters from the SDSS. Brightest group and cluster galaxies (BCGs) are more likely to host a radio-loud AGN than other galaxies of the same stellar mass (by less than a factor of two at a stellar mass of 5e11 M_sun, but over an order of magnitude below 1e11 M_sun). The distribution of radio luminosities for BCGs does not depend on stellar mass, however, and is similar to that of field galaxies of the same mass. Neither the radio-loud fraction nor the radio luminosity distribution of BCGs depends strongly on the cluster velocity dispersion. The radio-AGN fraction is also studied as a function of clustercentric distance: only within 0.2 r_200 do cluster galaxies exhibit an enhanced likelihood of radio-loud AGN activity, approaching that of the BCGs. The radio-loud AGN properties of both BCGs and non-BCGs can naturally be explained if this activity is fuelled by cooling from hot gas surrounding the galaxy. Using observational estimates of the mechanical output of the radio jets, the time-averaged energy output associated with recurrent radio source activity is estimated for all group/cluster galaxies. Within the cooling radius of the cluster, the radio-mode heating from the BCG dominates over that of all other galaxies combined. The scaling between total radio-AGN energy output and cluster velocity dispersion is considerably shallower than the sigma^4 scaling of the radiative cooling rate. Thus, unless either the mechanical-to-radio luminosity ratio or the efficiency of converting AGN mechanical energy into heating increases by 2-3 orders of magnitude between groups and rich clusters, radio-mode heating will not balance cooling in systems of all masses.

The structural parameters of bulges, disks and bars of a sample of nearly 1000 nearby galaxies are being determined through sophisticated image decomposition in the g, r and i bands. The sample is carefully drawn from the Sloan Digital Sky Survey Data Release 2 (SDSS DR2), contains 963 galaxies, of which 407 host AGN, and we show that it is representative of the galaxy and AGN populations in the local universe. A large number of other physical properties have also been determined for these galaxies within the SDSS collaboration. With these data, we reinforce several recent studies and find a number of results leading to new insights into how the different galaxy components relate to each other and the extent to which the hosts galaxies of AGN differ from quiescent galaxies.

We extend our previous physically-based halo occupation distribution models to include the dependence of clustering on the spectral energy distributions of galaxies. The high resolution Millennium Simulation is used to specify the positions and the velocities of the model galaxies. The stellar mass of a galaxy is assumed to depend only on M_{infall}, the halo mass when the galaxy was last the central dominant object of its halo. Star formation histories are parametrized using two additional quantities that are measured from the simulation for each galaxy: its formation time (t_{form}), and the time when it first becomes a satellite (t_{infall}). Central galaxies begin forming stars at time t_{form} with an exponential time scale tau_c. If the galaxy becomes a satellite, its star formation declines thereafter with a new time scale tau_s. We compute 4000 \AA break strengths for our model galaxies using stellar population synthesis models. By fitting these models to the observed abundances and projected correlations of galaxies as a function of break strength in the Sloan Digital Sky Survey, we constrain tau_c and tau_s as functions of galaxy stellar mass. We find that central galaxies with large stellar masses have ceased forming stars. At low stellar masses, central galaxies display a wide range of different star formation histories, with a significant fraction experiencing recent starbursts. Satellite galaxies of all masses have declining star formation rates, with similar e-folding times, tau_s ~ 2.5Gyr. One consequence of this long e-folding time is that the colour-density relation is predicted to flatten at redshifts > 1.5, because star formation in the majority of satellites has not yet declined by a significant factor. This is consistent with recent observational results from the DEEP and VVDS surveys.

We use the Millennium Simulation, a 10 billion particle simulation of the growth of cosmic structure, to construct a new model of galaxy clustering. We adopt a methodology that falls midway between the traditional semi-analytic approach and the halo occupation distribution(HOD) approach. In our model, we adopt the positions and velocities of the galaxies that are predicted by following the orbits and merging histories of the substructures in the simulation. Rather than using star formation and feedback `recipes' to specify the physical properties of the galaxies, we adopt parametrized functions to relate these properties to the quantity M_{infall}, defined as the mass of the halo at the epoch when the galaxy was last the central dominant object in its own halo. We test whether these parametrized relations allow us to recover the basic statistical properties of galaxies in the semi-analytic catalogues, including the luminosity function, the stellar mass function and the shape and amplitude of the two-point correlation function evaluated in different stellar mass and luminosity ranges. We then use our model to interpret recent measurements of these quantities from Sloan Digital Sky Survey data. We derive relations between the luminosities and the stellar masses of galaxies in the local Universe and their host halo masses. Our results are in excellent agreement with recent determinations of these relations by Mandelbaum et al using galaxy-galaxy weak lensing measurements from the SDSS.

By comparing semi-analytic galaxy catalogues with data from the Sloan Digital Sky Survey (SDSS), we show that current galaxy formation models reproduce qualitatively the dependence of galaxy clustering and pairwise peculiar velocities on luminosity, but some subtle discrepancies with the data still remain. The comparisons are carried out by constructing a large set of mock galaxy redshift surveys that have the same selection function as the SDSS Data Release Four (DR4). The mock surveys are based on two sets of semi-analytic catalogues presented by Croton et al. and Kang et al. From the mock catalogues, we measure the redshift-space projected two-point correlation function, the power spectrum, and the pairwise velocity dispersion (PVD) in Fourier space and in configuration space, for galaxies in different luminosity intervals. We then compare these theoretical predictions with the measurements derived from the SDSS DR4. On large scales and for galaxies brighter than L*, both sets of mock catalogues agree well with the data. For fainter galaxies, however, both models predict stronger clustering and higher pairwise velocities than observed. We demonstrate that this problem can be resolved if the fraction of faint satellite galaxies in massive haloes is reduced by ~30% compared to the model predictions. A direct look into the model galaxy catalogues reveals that a signifcant fraction (15%) of faint galaxies ($-18<M_{^{0.1}r}<-17$) reside in haloes with $M_{vir}>10^{13}\msun$, and this population is predominantly red in colour. These faint red galaxies are responsible for the high PVD values of low-luminosity galaxies on small scales.

We have analyzed the clustering of ~ 90,000 narrow-line AGN drawn from the Data Release 4 (DR4) of the SDSS. We compute the cross-correlation between AGN and a reference sample of galaxies, and compare this to results for control samples of inactive galaxies matched simultaneously in redshift,stellar mass,concentration, velocity dispersion and the 4000A break strength. We also compare near-neighbour counts around AGN and around the control galaxies. On scales larger than a few Mpc, AGN have almost the same clustering amplitude as the control sample. This demonstrates that AGN host galaxies and inactive galaxies populate dark matter halos of similar mass.On scales between 100kpc and 1Mpc,AGN are clustered more weakly than the control galaxies. We use mock catalogues constructed from high-resolution N-body simulations to interpret this anti-bias, showing that the observed effect is easily understood if AGN are preferentially located at the centres of their dark matter halos. On scales less than 70 kpc, AGN cluster marginally more strongly than the control sample, but the effect is weak. When compared to the control sample, we find that only one in a hundred AGN has an extra neighbour within a radius of 70 kpc. This excess increases as a function of the accretion rate onto the black hole, but it does not rise above the few percent level. Although interactions between galaxies may be responsible for triggering nuclear activity in a minority of nearby AGN, some other mechanism is required to explain the activity seen in the majority of the objects in our sample. (abridged)

We analyze a volume-limited sample of massive bulge-dominated galaxies with data from both the Sloan Digital Sky Survey and the Galaxy Evolution Explorer (GALEX) satellite. The galaxies have central velocity dispersions greater than 100 km/s and stellar surface mass densities that lie above the value where galaxies transition from actively star forming to passive systems. The sample is limited to redshifts 0.03<z<0.07. At these distances, the SDSS spectra sample the light from the bulge-dominated central regions of the galaxies. The GALEX NUV data provide high sensitivity to low rates of global star formation in these systems. Our sample of bulge-dominated galaxies exhibits a much larger dispersion in NUV-r colour than in optical g-r colour. Nearly all of the galaxies with bluer NUV-r colours are AGN. Both GALEX images and SDSS colour profiles demonstrate that the excess UV light is associated with an extended disk. We find that galaxies with red outer regions almost never have a young bulge or a strong AGN. Galaxies with blue outer regions have bulges and black holes that span a wide range in age and accretion rate. Galaxies with young bulges and strongly accreting black holes almost always have blue outer disks. Our suggested scenario is one in which the source of gas that builds the bulge and black hole is a low mass reservoir of cold gas in the disk.The presence of this gas is a necessary, but not sufficient condition for bulge and black hole growth. Some mechanism must transport this gas inwards in a time variable way. As the gas in the disk is converted into stars, the galaxies will turn red, but further gas infall can bring them back into the blue NUV-r sequence.(Abridged)

Current emission-line based estimates of the metallicity of active galactic nuclei (AGN) at both high and low redshifts indicate that AGN have predominantly solar to supersolar metallicities. This leads to the question: do low metallicity AGN exist? In this paper we use photoionization models to examine the effects of metallicity variations on the narrow emission lines from an AGN. We explore a variety of emission-line diagnostics that are useful for identifying AGN with low metallicity gas. We find that line ratios involving [NII] are the most robust metallicity indicators in galaxies where the primary source of ionization is from the active nucleus. Ratios involving [SII] and [OI] are strongly affected by uncertainties in modelling the density structure of the narrow line clouds. To test our diagnostics, we turn to an analysis of AGN in the Sloan Digital Sky Survey (SDSS). We find a clear trend in the relative strength of [NII] with the mass of the AGN host galaxy. The metallicity of the ISM is known to be correlated with stellar mass in star-forming galaxies; our results indicate that a similar trend exists for AGN. We also find that the best-fit models for typical Seyfert narrow line regions have supersolar abundances. Although there is a mass-dependent range of a factor of 2-3 in the NLR metallicities of the AGN in our sample, AGN with sub-solar metallicities are very rare in the SDSS. Out of a sample of ~23000 Seyfert 2 galaxies we find only ~40 clear candidates for AGN with NLR abundances that are below solar.

We simulate the growth of galaxies and their central supermassive black holes by implementing a suite of semi-analytic models on the output of the Millennium Run, a very large simulation of the concordance LCDM cosmogony. Our procedures follow the detailed assembly history of each object and are able to track the evolution of all galaxies more massive than the Small Magellanic Cloud throughout a volume comparable to that of large modern redshift surveys. In this first paper we supplement previous treatments of the growth and activity of central black holes with a new model for `radio' feedback from those AGN that lie at the centre of a quasistatic X-ray emitting atmosphere in a galaxy group or cluster. We show that for energetically and observationally plausible parameters such a model can simultaneously explain: (i) the low observed mass drop-out rate in cooling flows; (ii) the exponential cut-off at the bright end of the galaxy luminosity function; and (iii) the fact that the most massive galaxies tend to be bulge-dominated systems in clusters and to contain systematically older stars than lower mass galaxies. This success occurs because static hot atmospheres form only in the most massive structures, and radio feedback (in contrast, for example, to supernova or starburst feedback) can suppress further cooling and star formation without itself requiring star formation. We discuss possible physical models which might explain the accretion rate scalings required for our phenomenological `radio mode' model to be successful.

We synthesize the rest-frame Stroemgren colors using SDSS spectra for 99,088 galaxies selected from Data Release 1. This narrow-band ~200 AA photometric system (uz, vz, bz, yz), first designed for the determination of effective temperature, metallicity and gravity of stars, measures the continuum spectral slope of galaxies in the rest-frame 3200-5800 AA wavelength range. Galaxies form a remarkably narrow locus (~0.03 mag) in the resulting color-color diagram. The Bruzual & Charlot population synthesis models suggest that the position of a galaxy along this locus is controlled by a degenerate combination of metallicity and age of the dominant stellar population. Galaxy distribution along the locus is bimodal, with the local minimum corresponding to an ~1 Gyr old single stellar population. The position perpendicular to the locus is independent of metallicity and age, and reflects the galaxy's dust content, as implied by both the models and the statistics of IRAS detections. A comparison of this locus with the galaxy locus in the H_delta-D_n(4000) diagram, utilized by Kauffmann et al. (2003) to estimate stellar masses, reveals a tight correlation, although the two analyzed spectral ranges barely overlap. Overall, the galaxy spectral energy distribution in the entire UV to near-IR range can be described as a single-parameter family with an accuracy of 0.1 mag, or better. This nearly one-dimensional distribution of galaxies in the multi-dimensional space of measured parameters strongly supports the conclusion of Yip et al. (2004), based on a principal component analysis, that SDSS galaxy spectra can be described by a small number of eigenspectra. Apparently, the contributions of stellar populations that dominate the optical emission from galaxies are combined in a simple and well-defined way.

We present an analysis of the host properties of 85224 emission-line galaxies from the Sloan Digital Sky Survey. We show that Seyferts and LINERs form clearly separated branches on the standard optical diagnostic diagrams. We derive a new empirical classification scheme which cleanly separates star-forming galaxies, composite AGN-HII galaxies, Seyferts and LINERs and we study the host galaxy properties of these different classes of objects. LINERs are older, more massive, less dusty and more concentrated, and have higher velocity dispersions and lower [OIII] luminosities than Seyfert galaxies. We consider the quantity L[OIII]/sigma^4, which is an indicator of the black hole accretion rate relative to the Eddington rate. Remarkably, we find that at fixed L[OIII]/sigma^4, all differences between Seyfert and LINER host properties disappear. LINERs and Seyferts form a continuous sequence, with LINERs dominant at low L/L(EDD) and Seyferts dominant at high L/L(EDD). These results suggest that the majority of LINERs are AGN and that the Seyfert/LINER dichotomy is analogous to the high/low-state transition for X-ray binary systems. We apply theoretical photoionization models and show that pure LINERs require a harder ionizing radiation field with lower ionization parameter than Seyfert galaxies, consistent with the low and high X-ray binary states.

We combine data from SDSS and the GALEX and Spitzer observatories to create a sample of galaxies observed homogeneously from the UV to the Far-IR. This sample, consisting of ~460 galaxies observed spectroscopically by SDSS provides a multiwavelength (0.15-24 micron) view of obscured and unobscured star formation in nearby (z<0.3) galaxies with SFRs ranging from 0.01 to 100 M_solar/yr. We calculate a robust dust measure from the infrared to UV ratio (IRX) and explore the influence of star formation history (SFH) on the dust-UV color relation (i.e. the IRX-beta relation). We find that the UV colors of galaxies are only weakly dependent on their SFH as measured by the 4000A break. However, we find that the contributions of dust and SFH are distinguishable when colors at widely separated wavelengths (e.g. 0.23-3.6 micron) are introduced. We show this explicitly by recasting the IRX-beta relation as a more general IRX-SFH-color relation, which we examine in different projections. We also determine simple fits to this relation.

The relationship between galaxies and dark matter can be characterized by the halo mass of the central galaxy and the fraction of galaxies that are satellites. Here we present observational constraints from the SDSS on these quantities as a function of r-band luminosity and stellar mass using galaxy-galaxy weak lensing, with a total of 351,507 lenses. We use stellar masses derived from spectroscopy and virial halo masses derived from weak gravitational lensing to determine the efficiency with which baryons in the halo of the central galaxy have been converted into stars. We find that an L* galaxy with a stellar mass of 6x10^{10} M_{sun} is hosted by a halo with mass of 1.4x10^{12} M_{sun}/h, independent of morphology, yielding baryon conversion efficiencies of 17_{-5}^{+10} (early types) and 16_{-6}^{+15} (late types) per cent at the 95 per cent CL (statistical, not including systematic uncertainty due to assumption of a universal initial mass function, or IMF). We find that for a given stellar mass, the halo mass is independent of morphology below M_{stellar}=10^{11} M_{sun}, in contrast to typically a factor of two difference in halo mass between ellipticals and spirals at a fixed luminosity. This suggests that stellar mass is a good proxy for halo mass in this range and should be used preferentially whenever a halo mass selected sample is needed. For higher stellar masses, the conversion efficiency is a declining function of stellar mass, and the differences in halo mass between early and late types become larger, reflecting the fact that most group and cluster halos with masses above 10^{13} M_{sun} host ellipticals at the center, while even the brightest central spirals are hosted by halos of mass below 10^{13} M_{sun}. (Abridged)

(abridged) We present measurements of the pairwise velocity dispersion (PVD) for different classes of galaxies in the Sloan Digital Sky Survey. For a sample of about 200,000 galaxies, we study the dependence of the PVD on galaxy properties such as luminosity, stellar mass (M_*), colour (g-r), 4000A break strength (D4000), concentration index (C), and stellar surface mass density (\mu_*). The luminosity dependence of the PVD is in good agreement with the results of Jing & B\"orner (2004) for the 2dFGRS catalog. The value of \sigma_{12} measured at k=1 h/Mpc decreases as a function of increasing galaxy luminosity for galaxies fainter than L*, before increasing again for the most luminous galaxies in our sample. Each of the galaxy subsamples selected according to luminosity or stellar mass is divided into two further subsamples according to colour, D4000, C and \mu_*. We find that galaxies with redder colours and higher D4000, C, and \mu_* values have larger PVDs on all scales and at all luminosities/stellar masses. The dependence of the PVD on parameters related to recent star formation(colour, D4000) is stronger than on parameters related to galaxy structure (C, \mu_*), especially on small scales and for faint galaxies. The reddest galaxies and galaxies with high surface mass densities and intermediate concentrations have the highest pairwise peculiar velocities, i.e. these move in the strongest gravitational fields. We conclude that the faint red population located in rich clusters is responsible for the high PVD values that are measured for low-luminosity galaxies on small scales.

(abridged)We use a sample of ~200,000 galaxies drawn from the Sloan Digital Sky Survey to study how clustering depends on properties such as stellar mass (M*), colour (g-r), 4000A break strength (D4000), concentration index (C), and stellar surface mass density (\mu_*). We find that more massive galaxies cluster more strongly than less massive galaxies, with the difference increasing above the characteristic stellar mass of the Schechter mass function. When divided by physical quantities, galaxies with redder colours, larger D4000, higher C and larger \mu_* cluster more strongly. The clustering differences are largest on small scales and for low mass galaxies. At fixed stellar mass,the dependences of clustering on colour and 4000A break strength are similar. Different results are obtained when galaxies are split by concentration or surface density. The dependence of w(r_p) on g-r and D4000 extends out to physical scales that are significantly larger than those of individual dark matter haloes (> 5 Mpc/h). This large-scale clustering dependence is not seen for the parameters C or \mu_*. On small scales (< 1 Mpc/h), the amplitude of the correlation function is constant for ``young'' galaxies with 1.1 < D4000< 1.5 and a steeply rising function of age for ``older'' galaxies with D4000>1.5. In contrast, the dependence of the amplitude of w(r_p) on concentration on scales less than 1 Mpc/h is strongest for disk-dominated galaxies with C<2.6. This demonstrates that different processes are required to explain environmental trends in the structure and in star formation history of galaxies.

We use the ROSAT All Sky Survey (RASS) to study the soft X-ray properties of a homogeneous sample of 46,420 quasars selected from the third data release of the Sloan Digital Sky Survey (SDSS). Optical luminosities, both at rest-frame 2500\AA ($L_{2500}$) and in [OIII] ($L_{[\rm{OIII}]}$) span more than three orders of magnitude, while redshifts range over $0.1<z<5.4$. We detect 3366 quasars directly in the observed 0.1--2.4 keV band. Sub-samples of radio-loud and radio-quiet objects (RLQs and RQQs) are obtained by cross-matching with the FIRST catalogue. We study the distribution of X-ray luminosity as a function of optical luminosity, redshift and radio power using both individual detections and stacks of complete sets of similar quasars. At every optical luminosity and redshift $\log L_{2\kev}$ is, to a good approximation, normally distributed with dispersion $\sim 0.40$, at least brightwards of the median X-ray luminosity. This median X-ray luminosity of quasars is a power law of optical luminosity with index $\sim 0.53$ for $L_{2500}$ and $\sim 0.30$ for $L_{[\rm{OIII}]}$. RLQs are systematically brighter than RQQs by about a factor of 2 at given optical luminosity. The zero-points of these relations increase systematically with redshift, possibly in different ways for RLQs and RQQs. Evolution is particularly strong at low redshift and if the optical luminosity is characterised by $L_{[\rm{OIII}]}$. At low redshift and at given $L_{[\rm{OIII}]}$ the soft X-ray emission from type II AGN is more than 100 times weaker than that from type I AGN.

A long-standing problem for models of galaxy formation has been the mismatch between the predicted shape of the mass function of dark matter halos and the observed shape of the luminosity function of galaxies. The number of massive halos is predicted to decrease as a power law (N proportional to M^-2) out to very large masses, while the galaxy luminosity function cuts off exponentially at luminosities above L*. This implies that the efficiency with which gas cools onto massive systems is lower than expected. This letter investigates the role of radio-loud active galactic nuclei (AGN) in continually re-heating the cooling gas. By combining two observational results, the time-averaged energy output associated with recurrent radio source activity is determined, as a function of the black hole mass of the host galaxy: H = 10^21.4 (M_BH / M_sun)^1.6 W. It is shown that for massive elliptical galaxies this radio-source heating balances the radiative energy losses from the hot gas surrounding the galaxy. The recurrent radio-loud AGN activity may therefore provide a self-regulating feedback mechanism capable of controlling the rate of growth of galaxies.

In Figure 6 we inadvertently labeled the proxy circular velocity as the virial velocity of the dark matter halo instead of what is actually plotted, the maximum circular velocity of the dark matter halo. The maximum halo circular velocity is a much better estimate of the disk V_c than is V_vir. This confusion influenced the discussion of the Tully-Fisher relation in our paper. In fact, Figure 6 demonstrates that it is possible to simultaneously reproduce both the local Tully-Fisher relation and luminosity function using semi-analytic techniques applied to the standard LCDM cosmology, thus contradicting previous studies of this issue and our own discussion in Section 3.6.

We take advantage of the largest high-resolution simulation of cosmic structure growth ever carried out -- the Millennium Simulation of the concordance LambdaCDM cosmogony -- to study how the star formation histories, ages and metallicities of elliptical galaxies depend on environment and on stellar mass. We concentrate on a galaxy formation model which is tuned to fit the joint luminosity/colour/morphology distribution of low redshift galaxies. Massive ellipticals in this model have higher metal abundances, older luminosity-weighted ages, shorter star formation timescales, but lower assembly redshifts than less massive systems. Within clusters the typical masses, ages and metal abundances of ellipticals are predicted to decrease, on average, with increasing distance from the cluster centre. We also quantify the effective number of progenitors of ellipticals as a function of present stellar mass, finding typical numbers below 2 for M* < 10^{11} Msun, rising to about 5 for the most massive systems. These findings are consistent with recent observational results that suggest ``down-sizing'' or ``anti-hierarchical'' behaviour for the star formation history of the elliptical galaxy population, despite the fact that our model includes all the standard elements of hierarchical galaxy formation and is implemented on the standard, LambdaCDM cosmogony.

We have compiled a sample of 2728 nearby ($z<0.08$) elliptical galaxies with photometry in the $g,r,i,z$ bands from the Sloan Digital Sky Survey (SDSS) and $J,H,K$ photometry from the Two Micron All Sky Survey (2MASS). Stellar masses, stellar velocity dispersions and structural parameters such as sizes and surface mass densities are also available for these objects. In order to correct the aperture mismatch between SDSS and 2MASS, we correct the SDSS magnitudes to the isophotal circular radius where the 2MASS magnitudes are measured. We compare the correlations between optical, optical-infrared and infrared colours and galaxy luminosity, stellar mass, velocity dispersion and surface mass density. We find that all galaxy colours correlate more strongly with stellar mass and velocity dispersion than with any other structural parameter. The dispersion about these two relations is also smaller. We also study the correlations between a variety of stellar absorption line indices and the same set of galaxy parameters and we reach very similar conclusions. Finally, we analyze correlations between absorption line indices and colour. Our results suggest that the optical colours of elliptical galaxies are sensitive to a combination of age, metallicity and $\alpha$-enhancement, while the optical-infrared colours are sensitive to metallicity and to $\alpha$-enhancement, but are somewhat less sensitive to age.

We study the recent star formation histories of local galaxies by analyzing the scatter in their colours and spectral properties. We present evidence that the distribution of star formation histories changes qualitatively above a characteristic stellar surface mass density of 3x10^8 M_sol kpc^-2, corresponding to the transition between disk-dominated and bulge-dominated galaxies. When we average over subpopulations of galaxies with densities below this value, we find that subpopulations of all masses form their stars at the same average rate per unit stellar mass. However, the scatter in galaxy colours and emission line strengths is larger for more compact galaxies of a given mass. This suggests that star formation occurs in shorter, higher amplitude events in galaxies with smaller sizes. Above the characteristic density, galaxy growth through star formation shuts down. We propose that star formation events are triggered when cold gas is accreted onto a galaxy. We have used the new high resolution Millennium Simulation of structure formation in a concordance LCDM Universe to quantify the incidence of these accretion events and we show that the observational data are well fit by a model in which the consumption time of accreted gas decreases with the surface density of the galaxy as t_cons \propto mu*^-1. The dark matter halos hosting massive galaxies with high densities are also expected to grow through accretion, but the observations indicate that star formation in bulge- dominated galaxies is no longer coupled to the hierarchical build-up of their surrounding dark matter halos.

Observations have shown that the Eddington ratios (the ratio of the bolometric luminosity to the Eddington luminosity) in QSOs/active galactic nuclei (AGNs) cover a wide range. In this paper we connect the demography of AGNs obtained by the Sloan Digital Sky Survey with the accretion physics around massive black holes and propose that the diversity in the Eddington ratios is a natural result of the long-term evolution of accretion disks in AGNs. The observed accretion rate distribution of AGNs (with host galaxy velocity dispersion sigma~70-200 km/s) in the nearby universe (z<0.3) is consistent with the predictions of simple theoretical models in which the accretion rates evolve in a self-similar way. We also discuss the implications of the results for the issues related to self-gravitating disks, coevolution of galaxies and QSOs/AGNs, and the unification picture of AGNs.

We present the ESO Distant Cluster Survey (EDisCS) a survey of 20 fields containing distant galaxy clusters with redshifts ranging from 0.4 to almost 1.0. Candidate clusters were chosen from among the brightest objects identified in the Las Campanas Distant Cluster Survey, half with estimated redshift z_est~0.5 and half with z_est~0.8. They were confirmed by identifying red sequences in moderately deep two colour data from VLT/FORS2. For confirmed candidates we have assembled deep three-band optical photometry using VLT/FORS2, deep near-infrared photometry in one or two bands using NTT/SOFI, deep optical spectroscopy using VLT/FORS2, wide field imaging in two or three bands using the ESO Wide Field Imager, and HST/ACS mosaic images for 10 of the most distant clusters. This first paper presents our sample and the VLT photometry we have obtained. We present images, colour-magnitude diagrams and richness estimates for our clusters, as well as giving redshifts and positions for the brightest cluster members. Subsequent papers will present our infrared photometry, spectroscopy, HST and wide-field imaging, as well as a wealth of further analysis and science results. Our reduced data become publicly available as these papers are accepted.

In this paper we assess the relationship of the population of Active Galactic Nuclei (AGN) selected by hard X-rays to the traditional population of AGN with strong optical emission lines. First, we study the emission-line properties of a new hard X-ray selected sample of 47 local AGN (classified optically as both Type 1 and Type 2 AGN). We find that the hard X- ray (3-20 keV) and [OIII]$\lambda$5007 optical emission-line luminosities are well-correlated over a range of about four orders-of-magnitude in luminosity (mean luminosity ratio 2.15 dex with a standard deviation of $\sigma$ = 0.51 dex). Second, we study the hard X-ray properties of a sample of 55 local AGN selected from the literature on the basis of the flux in the [OIII] line. The correlation between the hard X-ray (2-10 keV) and [OIII] luminosity for the Type 1 AGN is consistent with what is seen in the hard X-ray selected sample. However, the Type 2 AGN have a much larger range in the luminosity ratio, and many are very weak in hard X-rays (as expected for heavily absorbed AGN). We then compare the hard X-ray (3-20 keV) and [OIII] luminosity functions of AGN in the local universe. These have similar faint-end slopes with a luminosity ratio of 1.60 dex (0.55 dex smaller than the mean value for individual hard X-ray selected AGN). We conclude that at low redshift, selection by narrow optical emission- lines will recover most AGN selected by hard X-rays (with the exception of BL Lac objects). However, selection by hard X-rays misses a significant fraction of the local AGN population with strong emission lines.

We have cross-correlated the Sloan Digital Sky Survey (SDSS) second data release spectroscopic galaxy sample with the IRAS faint-source catalogue (FSC). Optical emission line ratios are used to classify the galaxies with reliable IRAS 60 and 100 microns detections into AGN and normal star-forming galaxies. We then create subsamples of normal galaxies and AGN that are very closely matched in terms of key physical properties such as stellar mass, redshift, size, concentration and mean stellar age (as measured by absorption line indicators in the SDSS spectra). We then quantify whether there are systematic differences between the IR luminosities of the galaxies and the AGN in the matched subsamples. We find that the AGN exhibit a significant excess in far-IR emission relative to the star-forming galaxies in our sample. The excesses at 60 and 100 microns are 0.21 +/- 0.03 dex and 0.12 +/- 0.035 dex in log[L(60)/M*] and log[L(100)/M*], respectively. We then discuss whether the far-IR excess is produced by radiation from the active nucleus that is absorbed by dust or alternatively, by an extra population of young stars that is not detectable at optical wavelengths.

The properties of the host galaxies of a well-defined sample of 2215 radio-loud AGN with redshifts 0.03 < z < 0.3, defined from the SDSS, are investigated. These are predominantly low radio luminosity sources, with 1.4GHz luminosities of 10^23 to 10^25 W/Hz. The fraction of galaxies that host radio-loud AGN with L(1.4GHz) > 10^23 W/Hz is a strong function of stellar mass, rising from nearly zero below a stellar mass of 10^10 Msun to more than 30% at 5x10^11 Msun. The integral radio luminosity function is derived in six ranges of stellar and black hole mass. Its shape is very similar in all of these ranges and can be well fitted by a broken power-law. Its normalisation varies strongly with mass, as M_*^2.5 or M_BH^1.6; this scaling only begins to break down when the predicted radio-loud fraction exceeds 20-30%. There is no correlation between radio and emission line luminosities for the radio-loud AGN in the sample and the probability that a galaxy of given mass is radio-loud is independent of whether it is optically classified as an AGN. The host galaxies of the radio-loud AGN have properties similar to those of ordinary galaxies of the same mass. All of these findings support the conclusion that the optical AGN and low radio luminosity AGN phenomena are independent and are triggered by different physical mechanisms. Intriguingly, the dependence on black hole mass of the radio-loud AGN fraction mirrors that of the rate at which gas cools from the hot atmospheres of elliptical galaxies. It is speculated that gas cooling provides a natural explanation for the origin of the radio-loud AGN activity, and it is argued that AGN heating could plausibly balance the cooling of the gas over time. [Abridged]

A sample of 2712 radio-luminous galaxies is defined from the second data release of the Sloan Digital Sky Survey (SDSS) by cross-comparing the main spectroscopic galaxy sample with two radio surveys: the NRAO VLA Sky Survey (NVSS) and the Faint Images of the Radio Sky at Twenty centimetres (FIRST) survey. The comparison is carried out in a multi-stage process and makes optimal use of both radio surveys by exploiting the sensitivity of the NVSS to extended and multi-component radio sources in addition to the high angular resolution of the FIRST images. A radio source sample with 95% completeness and 98.9% reliability is achieved, far better than would be possible for this sample if only one of the surveys was used. The radio source sample is then divided into two classes: radio-loud AGN and galaxies in which the radio emission is dominated by star formation. The division is based on the location of a galaxy in the plane of 4000Ang break strength versus radio luminosity per unit stellar mass and provides a sample of 2215 radio-loud AGN and 497 star forming galaxies brighter than 5mJy at 1.4GHz. A full catalogue of positions and radio properties is provided for these sources. The local radio luminosity function is then derived both for radio-loud AGN and for star-forming galaxies and is found to be in agreement with previous studies. By using the radio to far-IR correlation, the radio luminosity function of star forming galaxies is also compared to the luminosity function derived in the far-infrared. It is found to agree well at high luminosities but less so at lower luminosities, confirming that the linearity of the radio to far-IR correlation breaks down below about 10^22 W/Hz at 1.4GHz.

The emission line luminosity function of active galactic nuclei (AGN) is measured from about 3000 AGN included in the main galaxy sample of the Sloan Digital Sky Survey within a redshift range of $0<z<0.15$. The $\Ha$ and [OIII]$\lambda 5007$ luminosity functions for Seyferts cover luminosity range of $10^{5-9}$$L_\odot$ in H$\alpha$ and the shapes are well fit by broken power laws, without a turnover at fainter nuclear luminosities. Assuming a universal conversion from emission line strength to continuum luminosity, the inferred B band magnitude luminosity function is comparable both to the AGN luminosity function of previous studies and to the low redshift quasar luminosity function derived from the 2dF redshift survey. The inferred AGN number density is approximately 1/5 of all galaxies and about $6\times 10^{-3}$ of the total light of galaxies in the $r$-band comes from the nuclear activity. The numbers of Seyfert 1s and Seyfert 2s are comparable at low luminosity, while at high luminosity, Seyfert 1s outnumber Seyfert 2s by a factor of 2-4. In making the luminosity function measurements, we assumed that the nuclear luminosity is independent of the host galaxy luminosity, an assumption we test {\it a posteriori}, and show to be consistent with the data. Given the relationship between black hole mass and host galaxy bulge luminosity, the lack of correlation between nuclear and host luminosity suggests that the main variable that determines the AGN luminosity is the Eddington ratio, not the black hole mass. This appears to be different from luminous quasars, which are most likely to be shining near the Eddington limit.

We have compiled a large sample of low-redshift active galactic nuclei (AGN) identified via their emission line characteristics from the spectroscopic data of the Sloan Digital Sky Survey. Since emission lines are often contaminated by stellar absorption lines, we developed an objective and efficient method of subtracting the stellar continuum from every galaxy spectrum before making emission line measurements. The distribution of the measured H$\alpha$ Full Width at Half Maxima values of emission line galaxies is strongly bimodal, with two populations separated at about 1,200km s$^{-1}$. This feature provides a natural separation between narrow-line and broad-line AGN. The narrow-line AGN are identified using standard emission line ratio diagnostic diagrams. 1,317 broad-line and 3,074 narrow-line AGN are identified from about 100,000 galaxy spectra selected over 1151 square degrees. This sample is used in a companion paper to determine the emission-line luminosity function of AGN.

We have used the first matched set of GALEX and SDSS data to investigate the properties of a sample of 74 nearby galaxies with far-ultraviolet luminosities chosen to overlap the luminosity range of typical high-z Lyman Break Galaxies (LBGs). GALEX deep surveys have shown that ultraviolet-luminous galaxies (UVLGs) similar to these are the fastest evolving component of the UV galaxy population. Model fits to the combined GALEX and SDSS photometry yield typical FUV extinctions similar to LBGs. The implied star formation rates are SFR ~ 3 to 30 solar mass per year. This overlaps the range of SFRs for LBGs. We find a strong inverse correlation between galaxy mass and far-ultraviolet surface brightness, and on this basis divide the sample into ``large'' and ``compact'' UVLGs. The compact UVLGs have half-light radii of a few kpc or less (similar to LBGs). They are relatively low mass galaxies (~10 billion solar masses) with typical velocity dispersions of 60 to 150 km/s. They span a range in metallicity from 0.3 to 1 times solar, have blue optical-UV colors, and are forming stars at a rate sufficient to build the present galaxy in ~a Gigayear. In all these respects they appear similar to the LBG population. These ``living fossils'' may therefore provide an opportunity for detailed investigation of the physical processes occurring in typical star forming galaxies in the early universe.

We present a sample of over 50 luminous field bulges (including ellipticals) found in the Groth Strip Survey (GSS), with 0.73< z < 1.04 and with bulge magnitudes I <= 23. The exponential disk light is removed via decomposition of HST images using GIM2D. We find that 85% of these bulges are nearly as red as local E/S0's and have a shallow slope and a small color dispersion in the color-luminosity relation, suggesting roughly coeval formation. The surface brightnesses of these bulges are about 1 mag higher than local bulges. These results are explained adopting a "drizzling" scenario where a metal-rich early formation is later polluted by small amounts of additional star formation. Almost all disks have the same or bluer colors than their accompanying bulges, regardless of the bulge-disk ratio and bulge luminosity, as expected from semi-analytic hierarchical galaxy formation models. We present evidence that the few blue bulge candidates are not likely to be genuine blue ellipticals or bulges. Our deeper, more extensive, and less disk-contaminated observations challenge prior claims that 30% to 50% of field bulges or ellipticals are in a blue, star-forming phase at z < 1. We conclude that field bulges and ellipticals at z ~ 1, like luminous early- type cluster galaxies at the same redshift, are already dominated by metal-rich, old stellar populations that have been fading from a formation epoch earlier than z ~ 1.5. (abridged)

We present ultraviolet photometry for a sample of morphologically early-type galaxies selected by matching the Sloan Digital Sky Survey Data Release 1 with the GALEX Medium and All-sky Imaging Surveys. We obtain a working sample of 1032 early-type galaxies with GALEX FUV detections, SDSS spectroscopy, and $z<0.2$. Using the SDSS spectra to identify galaxies with even weak star formation or evidence of AGN, and further removing galaxies with any evidence of non early-type morphology, we derive a final sample of 172 red quiescent early-type galaxies. We find that the $FUV-r$ color has a full range of 5 mag. Plotting against the $FUV-r$ color the metallicity sensitive Lick $\rm Mg_2$ and D4000 indices, and the stellar velocity dispersion, we find no correlation between our measurement of UV rising flux, and any parameter sensitive to metallicity.

We derive a variety of physical parameters including star formation rates (SFRs), dust attenuation and burst mass fractions for 6472 galaxies observed by the Galaxy Evolution Explorer (GALEX) and present in the SDSS DR1 main spectroscopic sample. Parameters are estimated in a statistical way by comparing each observed broad-band SED (two GALEX and five SDSS bands) with an extensive library of model galaxy SEDs, which cover a wide range of star formation histories and include stochastic starbursts. We compare the constraints derived using SDSS bands only with those derived using the combination of SDSS and GALEX photometry. We find that the addition of the GALEX bands leads to significant improvement in the estimation of both the dust optical depth and the star formation rate over timescales of 100 Myr to 1 Gyr in a galaxy. We are sensitive to SFRs as low as 10^{-3} M_sun/yr, and we find that low levels of star formation (SF) are mostly associated with early-type, red galaxies. The least massive galaxies have ratios of current to past-averaged SF rates (b-parameter) consistent with constant SF over a Hubble time. For late-type galaxies, this ratio on average decreases with mass. We find that b correlates tightly with NUV-r color, implying that the SF history of a galaxy can be constrained on the basis of the NUV-r color alone. The fraction of galaxies that have undergone a significant starburst episode within the last 1 Gyr steeply declines with mass-from ~20% for galaxies with ~10^8 M_sun to ~5% for ~10^11 M_sun galaxies.

We examine the relationship between galaxies, supermassive black holes and AGN using a sample of 22,000 narrow-emission-line AGN drawn from a a sample of 122,000 galaxies from the Sloan Digital Sky Survey. We have studied how AGN host properties compare with those of normal galaxies and how they depend on the luminosity of the active nucleus. We find that AGN reside in massive galaxies and have distributions of sizes and concentrations that are similar to those of the early-type galaxies in our sample. The host galaxies of low-luminosity AGN have stellar populations similar to normal early-types. The hosts of high-luminosity AGN have much younger mean stellar ages and a significant fraction have experienced recent starbursts. High-luminosity AGN are also found in lower density environments. We use the stellar velocity dispersions of the AGN hosts to estimate black hole masses and their [OIII]$\lambda$5007 emission line luminosities to estimate black hole accretion rates. We find that the volume averaged ratio of star formation to black hole accretion is ~1000, in remarkable agreement with the observed ratio of stellar mass to black hole mass in nearby bulges. Our estimated accretion rates imply that low mass black holes are growing on a timescale that is comparable to the age of the Universe. The growth timescale increases by more than an order of magnitude for the most massive black holes in our sample. We conclude that the evolution of the AGN luminosity function is driven by a decrease in the characteristic mass scale of actively accreting black holes.

By providing homogeneous photometric and spectroscopic data of high quality for very large and objectively selected samples of galaxies, the Sloan Digital Sky Survey allows statistical studies of the physical properties of galaxies and AGN to be carried out at an unprecedented level of precision and detail. Here we publicly release catalogues of derived physical properties for 211,894 galaxies, including 33,589 narrow-line AGN. These are complete samples with well understood selection criteria drawn from the normal galaxy spectroscopic sample in the second SDSS public data release (DR2). We list properties obtained from the SDSS spectroscopy and photometry using modelling techniques presented in papers already published by our group. These properties include: stellar masses; stellar mass-to-light ratios; attenuation of the starlight by dust; indicators of recent major starbursts; current total and specific star-formation rates, both for the regions with spectroscopy and for the galaxies as a whole; gas-phase metallicities; AGN classifications based on the standard emission line ratio diagnostic diagrams and AGN [OIII] emission line luminosities. We also list our own measurements of absorption line indices and emission line fluxes from which these quantities were derived, together with a few quantities from the standard SDSS pipelines which play an important role in our work. Many other observed properties of these galaxies can be obtained from the SDSS DR2 catalogues themselves. We will add further physical properties to this release site as the relevant papers are accepted for publication. Catalogies containing these parameters may be accessed at http://www.mpa-garching.mpg.de/SDSS/.

We investigate the accretion-driven growth of supermassive black holes in the low-redshift Universe using 23,000 ``Type 2'' AGN and the complete sample of 123,000 galaxies in the SDSS from which they were drawn. We use the stellar velocity dispersions of the early type galaxies and AGN hosts to estimate their black hole masses and we use the AGN [OIII]5007 emission line luminosities to estimate black hole accretion rates. We find that most present-day accretion occurs onto black holes with masses < 10E8 Msun that reside in moderately massive galaxies (M ~ 10E10 to 10E11.5 Msun) with high stellar densities (~10E8.5 to 10E9.5 Msun/kpc^2) and young stellar populations. The volume-averaged accretion rates of low mass black holes (<10E7.5 Msun) imply that this population is growing on a timescale that is comparable to the age of the Universe. Around half this growth takes place in AGN that are radiating within a factor of five of the Eddington luminosity. Such systems are rare, making up only 0.2% of the low mass black hole population at the present day. The rest of the growth occurs in lower luminosity AGN. The growth timescale is more than two orders of magnitude longer for the population of the most massive black holes in our sample. The volume averaged ratio of star formation to black hole accretion in bulge-dominated galaxies is ~1000, in remarkable agreement with the observed ratio of stellar mass to black hole mass in nearby galaxy bulges. We conclude: a) that bulge formation and black hole formation are tightly coupled, even in present-day galaxies; and b) that the evolution of the AGN luminosity function documented in recent optical and x-ray surveys is driven by a decrease in the characteristic mass scale of accreting black holes.

We utilize Sloan Digital Sky Survey imaging and spectroscopy of ~53,000 star-forming galaxies at z~0.1 to study the relation between stellar mass and gas-phase metallicity. We derive gas-phase oxygen abundances and stellar masses using new techniques which make use of the latest stellar evolutionary synthesis and photoionization models. We find a tight (+/-0.1 dex) correlation between stellar mass and metallicity spanning over 3 orders of magnitude in stellar mass and a factor of 10 in metallicity. The relation is relatively steep from 10^{8.5} - 10^{10.5} M_sun, in good accord with known trends between luminosity and metallicity, but flattens above 10^{10.5} M_sun. We use indirect estimates of the gas mass based on the H-alpha luminosity to compare our data to predictions from simple closed box chemical evolution models. We show that metal loss is strongly anti-correlated with baryonic mass, with low mass dwarf galaxies being 5 times more metal-depleted than L* galaxies at z~0.1. Evidence for metal depletion is not confined to dwarf galaxies, but is found in galaxies with masses as high as 10^{10} M_sun. We interpret this as strong evidence both of the ubiquity of galactic winds and of their effectiveness in removing metals from galaxy potential wells.

We summarize the detection rates at wavelengths other than optical for \~99,000 galaxies from the Sloan Digital Sky Survey (SDSS) Data Release 1 ``main'' spectroscopic sample. The analysis is based on positional cross-correlation with source catalogs from ROSAT, 2MASS, IRAS, GB6, FIRST, NVSS and WENSS surveys. We find that the rest-frame UV-IR broad-band galaxy SEDs form a remarkably uniform, nearly one parameter, family. As an example, the SDSS u and r band data, supplemented with redshift, can be used to predict K band magnitudes measured by 2MASS with an rms scatter of only 0.2 mag; when measurement uncertainties are taken into account, the astrophysical scatter appears not larger than ~0.1 mag.

(modified) We present a comprehensive study of the physical properties of \~10^5 galaxies with measurable star formation in the SDSS. By comparing physical information extracted from the emission lines with continuum properties, we build up a picture of the nature of star-forming galaxies at z<0.2. We take out essentially all aperture bias using resolved imaging, allowing an accurate estimate of the total SFRs in galaxies. We determine the SFR density to be 1.915^{+0.02}_{-0.01}(rand.)^{+0.14}_{-0.42} (sys.) h70 10^{-2} Msun/yr/Mpc^3 at z=0.1 (for a Kroupa IMF) and we study the distribution of star formation as a function of various physical parameters. The majority of the star formation in the low redshift universe takes place in moderately massive galaxies (10^10-10^11 Msun), typically in HSB disk galaxies. Roughly 15% of all star formation takes place in galaxies that show some sign of an active nucleus. About 20% occurs in starburst galaxies. We show that the present to past-average star formation rate, the Scalo b-parameter; is almost constant over almost three orders of magnitude in mass, declining only at M*>10^10 Msun. The volume averaged b parameter is 0.408^{+0.005}_{-0.002} (rand).^{+0.029}_{-0.090} (sys.) h70^{-1}. We use this value constrain the star formation history of the universe. In agreement with other work we find a correlation between $b$ and morphological type, as well as a tight correlation between the 4000AA break (D4000) and b. We discuss how D4000 can be used to estimate b parameters for high redshift galaxies.

We use a complete sample of galaxies drawn from the SDSS to study how structure, star formation and nuclear activity depend on local density and on stellar mass. Local density is estimated by counting galaxies above a fixed absolute magnitude limit within cylinders 2 Mpc in projected radius and +-500 km/s in depth. The stellar mass distribution of galaxies shifts by nearly a factor of two towards higher masses between low and high density regions. At fixed stellar mass, both star formation and nuclear activity depend strongly on local density, while structural parameters such as size and concentration are almost independent of it. The galaxy property most sensitive to environment is specific star formation rate. For galaxies with stellar masses less than 3 x 10^10 M_sun, the median SFR/M* decreases by more than a factor of 10 from low to high densities. This decrease is less marked for massive galaxies. At fixed stellar mass, twice as many galaxies host AGN with strong [OIII] emission in low density regions as in high. Massive galaxies in low-density environments also contain more dust. We have analyzed correlations between spectroscopic indicators that probe SFH on different timescales (D4000, Hdelta_A and SFR/M*). The correlations do not depend on environment, suggesting that the decrease in star formation has ocurred over long (>1 Gyr) timescales. Since structure does not depend on environment for more massive galaxies, trends in recent SFH, dust and AGN for these systems cannot be driven by processes that alter structure, e.g. mergers. The SFH-density correlation is strongest for small scale (< 1 Mpc) estimates of local density. Finally, we highlight a striking similarity between changes in the galaxy population as a function of density and as a function of redshift and we interpret this using N-body simulations.

(Abridged) We use a combination of high resolution N-body simulations and semi-analytic techniques to follow the formation, the evolution and the chemical enrichment of cluster galaxies in a Lambda-CDM Universe. We model the transport of metals between the stars, the cold gas in galaxies, the hot gas in dark matter haloes, and the intergalactic gas outside virialized haloes. We have compared three different feedback schemes. The `retention' model assumes that material reheated by supernova explosions is able to leave the galaxy, but not the dark matter halo. The `ejection' model assumes that this material leaves the halo and is then re-incorporated when structure collapses on larger scales. The `wind' model uses prescriptions that are motivated by observations of local starburst galaxies. We require that our models reproduce the cluster LF from the 2dF survey, the relations between stellar mass, gas mass and metallicity inferred from new SDSS data, and the observed amount of metals in the ICM. With suitable adjustment of the free parameters in the model, a reasonable fit to the observational results at redshift zero can be obtained for all three feedback schemes. All three predict that the chemical enrichment of the ICM occurs at high redshift with 60-80 per cent of the metals currently in the ICM ejected at redshifts larger than 1. Massive galaxies are important contributors to the chemical pollution. The observed decline in baryon fraction from rich clusters to galaxy groups is reproduced only in an `extreme' ejection scheme, where material ejected from dark matter haloes is re-incorporated on a timescale comparable to the age of the Universe. We explore how the metal abundance in the intergalactic medium as a function of redshift can constraint how and when galaxies ejected their metals.

We analyse the properties of substructures within dark matter halos (subhalos) using a set of high-resolution numerical simulations of the formation of structure in a Lambda-CDM Universe. Our simulation set includes 11 high-resolution simulations of massive clusters as well as a region of mean density, allowing us to study the spatial and mass distribution of substructures down to a mass resolution limit of 10^9 h^(-1)Mo. We also investigate how the properties of substructures vary as a function of the mass of the `parent' halo in which they are located. We find that the substructure mass function depends at most weakly on the mass of the parent halo and is well described by a power-law. The radial number density profiles of substructures are steeper in low mass halos than in high mass halos. More massive substructures tend to avoid the centres of halos and are preferentially located in the external regions of their parent halos. We also study the mass accretion and merging histories of substructures, which we find to be largely independent of environment. We find that a significant fraction of the substructures residing in clusters at the present day were accreted at redshifts z < 1. This implies that a significant fraction of present-day `passive' cluster galaxies should have been still outside the cluster progenitor and more active at z~1.

We present two different projects. The first one is related to the development of a semi-analytic model to follow the formation, the evolution and chemical enrichment of cluster galaxies in a hierarchical dark matter model. The second project concerns the comparison between high resolution N-body simulations of clusters of galaxies with the EDisCS observational dataset.

We study the variation of the dark matter mass fraction of elliptical galaxies as a function of their luminosity, stellar mass, and size using a sample of 29,469 elliptical galaxies culled from the Sloan Digital Sky Survey. We model ellipticals as a stellar Hernquist profile embedded in an adiabatically compressed dark matter halo. This model allows us to estimate a dynamical mass ($M_{dynm}$) at the half-light radius from the velocity dispersion of the spectra, and to compare these to the stellar mass estimates ($M_{*}$) from Kauffmann et al (2003). We find that $M_{*}/L$ is independent of luminosity, while $M_{dynm}/L$ increases with luminosity, implying that the dark matter fraction increases with luminosity. We also observe that at a fixed luminosity or stellar mass, the dark matter fraction increases with increasing galaxy size or, equivalently, increases with decreasing surface brightness: high surface brightness galaxies show almost no evidence for dark matter, while in low surface brightness galaxies, the dark matter exceeds the stellar mass at the half light radius. We relate this to the fundamental plane of elliptical galaxies, suggesting that the tilt of this plane from simple virial predictions is due to the dark matter in galaxies. We find that a simple model where galaxies are embedded in dark matter halos and have a star formation efficiency independent of their surface brightness explains these trends. We estimate the virial mass of ellipticals as being approximately 7-30 times their stellar mass, with the lower limit suggesting almost all of the gas within the virial radius is converted into stars.

abridged: We use a complete sample of about 140,000 galaxies from the Sloan Digital Sky Survey (SDSS) to study the size distribution of galaxies and its dependence on their luminosity, stellar mass, and morphological type. The large SDSS database provides statistics of unprecedented accuracy. For each type of galaxy, the size distribution at given luminosity (or stellar mass) is well described by a log-normal function, characterized by its median $\bar{R}$ and dispersion $\sigma_{\ln R}$. For late-type galaxies, there is a characteristic luminosity at $M_{r,0}\sim -20.5$ (assuming $h=0.7$) corresponding to a stellar mass $M_0\sim 10^{10.6}\Msun$. Galaxies more massive than $M_0$ have $\bar{R}\propto M^{0.4}$ and $\sigma_{\ln R}\sim 0.3$, while less massive galaxies have $\bar{R}\propto M^{0.15}$ and $\sigma_{\ln R}\sim 0.5$. For early-type galaxies, the $\bar{R}$ - $M$ relation is significantly steeper, $\bar{R}\propto M^{0.55}$, but the $\sigma_{\ln R}$ - $M$ relation is similar to that of late-type galaxies. Faint red galaxies have sizes quite independent of their luminosities.

Previous modelling has demonstrated that it is difficult to reproduce the SCUBA source counts within the framework of standard hierarchical structure formation models if the sources are assumed to be the high-redshift counterparts of local ultra-luminous infrared galaxies with dust temperatures in the range 40-60 K. Here, we show that the counts are more easily reproduced in a model in which the bulk of the sub-millimetre emission comes from extended, cool (20-25 K) dust in objects with star formation rates of 50-100 solar masses per year. The low temperatures imply typical sizes of about 1 (S_{850}/1mJy})^{1/2} arcsec, a factor two to three larger than those predicted using starburst-like spectral energy distributions. Low dust temperatures also imply a ratio of optical/UV to 850-micron flux which is 30-100 times smaller, for the same optical depth, than expected for objects with a hot, starburst-like SED. This may help explain the small overlap between SCUBA sources and Lyman-break galaxies.

We examine the properties of the host galaxies of 22,623 narrow-line AGN with 0.02<z<0.3 selected from a complete sample of 122,808 galaxies from the Sloan Digital Sky Survey. We focus on the luminosity of the [OIII]$\lambda$5007 emission line as a tracer of the strength of activity in the nucleus. We study how AGN host properties compare to those of normal galaxies and how they depend on L[OIII]. We find that AGN of all luminosities reside almost exclusively in massive galaxies and have distributions of sizes, stellar surface mass densities and concentrations that are similar to those of ordinary early-type galaxies in our sample. The host galaxies of low-luminosity AGN have stellar populations similar to normal early-types. The hosts of high- luminosity AGN have much younger mean stellar ages. The young stars are not preferentially located near the nucleus of the galaxy, but are spread out over scales of at least several kiloparsecs. A significant fraction of high- luminosity AGN have strong H$\delta$ absorption-line equivalent widths, indicating that they experienced a burst of star formation in the recent past. We have also examined the stellar populations of the host galaxies of a sample of broad-line AGN. We conclude that there is no significant difference in stellar content between type 2 Seyfert hosts and QSOs with the same [OIII] luminosity and redshift. This establishes that a young stellar population is a general property of AGN with high [OIII] luminosities.

We develop a new method to constrain the star formation histories, dust attenuation and stellar masses of galaxies. It is based on two stellar absorption line indices, the 4000 Angstrom break strength and the Balmer absorption line index Hdelta_A. Together, these indices allow us to constrain the mean stellar ages of galaxies and the fractional stellar mass formed in bursts over the past few Gyr. A comparison with broad band photometry then yields estimates of dust attenuation and of stellar mass. We generate a large library of Monte Carlo realizations of different star formation histories, including starbursts of varying strength and a range of metallicities. We use this library to generate median likelihood estimates of burst mass fractions, dust attenuation strengths, stellar masses and stellar mass-to-light ratios for a sample of 122,208 galaxies drawn from the Sloan Digital Sky Survey. The typical 95% confidence range in our estimated stellar masses is +-40%. We study how the stellar mass-to-light ratios of galaxies vary as a function of absolute magnitude, concentration index and photometric pass-band and how dust attenuation varies as a function of absolute magnitude and 4000 Angstrom break strength. We also calculate how the total stellar mass of the present Universe is distributed over galaxies as a function of their mass, size, concentration, colour, burst mass fraction and surface mass density. We find that most of the stellar mass in the local Universe resides in galaxies that have stellar masses \~5\times 10^10 M_sol, half light radii ~3 kpc, and half- light surface mass densities ~10^9 M_sol/kpc^2. The distribution of D(4000) is strongly bimodal, showing a clear division between galaxies dominated by old stellar populations and galaxies with more recent star formation.

We study the relations between stellar mass, star formation history, size and internal structure for a complete sample of 122,808 galaxies drawn from the Sloan Digital Sky Survey. We show that low-redshift galaxies divide into two distinct families at a stellar mass of 3 \times 10^10 M_sol. Lower mass galaxies have young stellar populations, low surface mass densities and the low concentrations typical of disks. A significant fraction of the lowest mass galaxies in our sample have experienced recent starbursts. At given stellar mass, the sizes of low mass galaxies are log- normally distributed with dispersion sigma(ln R_50) \sim 0.5, in excellent agreement with the idea that they form with little angular momentum loss through cooling and condensation in a gravitationally dominant dark matter halo. Their median stellar surface mass density scales with stellar mass as mu* propto M_*^0.54, suggesting that the stellar mass of a disk galaxy is proprtional to the three halves power of its halo mass. This suggests that the efficiency of the conversion of baryons into stars in low mass galaxies increases in propor- tion to halo mass, perhaps as a result of supernova feedback processes. At stellar masses above 3 \times 10^10 M_sol, there is a rapidly increasing frac- tion of galaxies with old stellar populations, high surface mass densities and high concentrations typical of bulges. In this regime, the size distribution is log-normal, but its dispersion decreases rapidly with increasing stellar mass and the median mass surface density is approximately constant. This suggests that the star formation efficiency decreases in the highest mass halos, and that little star formation occurs in massive galaxies once they have assembled.

We study the number and interaction rates of supermassive black holes in galactic bulges as predicted by hierarchical models of galaxy formation in which the spheroidal components of galaxies are formed by mergers. In bright ellipticals, the number of events that can eject a central supermassive binary black hole is large. Central binaries must therefore merge in less than a Hubble time -- otherwise there will be too much scatter in the M-sigma relation and too many off-center galactic nuclei. We propose that binary black holes are able to merge during the major gas accretion events that trigger QSO activity in galaxies. If this is the case, less than 10 percent of faint ellipticals and 40 percent of bright ellipticals are predicted to harbour binary black holes with near equal masses at their centres. This binary may be ejected away from the centre of the galaxy or even into intergalactic space in up to 20 percent of the most luminous ellipticals. The number of low mass black holes that can interact with the central object is predicted to be a strong function of galaxy luminosity. In most faint ellipticals, no black holes fall into the centre of the galaxy after the last major gas accretion event, but in the most luminous ellipticals, an average of 10 low mass black holes interact with the central supermassive object after this time. It is expected that stars will be ejected from galaxy cores as these low mass-ratio binaries harden. Multiple black holes in galactic bulges thus provide a natural explanation for the strong systematic trends in the observed central density profiles of ellipticals as a function of luminosity.

McKay et al. (2002) have recently used measurements of weak galaxy-galaxy lensing in the Sloan Digital Sky Survey to estimate the cross correlation between galaxies and the projected dark matter density field. They derive a relation between aperture mass within a radius of $260\kpch$, $M_{260}$, and lens galaxy luminosity, that does not depend on galaxy luminosity, type or environment. In this paper, we study the cross-correlation between galaxies and dark matter using galaxy catalogs constructed from a high-resolution N-body simulation of a $\Lambda$CDM Universe. We show that our simulations reproduce the McKay et al. results reasonably well. In the simulation, $M_{260}$ is approximately equal to the halo virial mass for $L_*$ galaxies. $M_{260}$ overestimates the virial mass for fainter galaxies and underestimates it for brighter galaxies. We use the simulations to show that under certain circumstances the halo virial mass may be recovered by fitting an NFW profile to the projected galaxy-mass correlation function. If we apply our method to the observations, we find that $L_*$ galaxies typically reside in halos of $\sim 2\times 10^{12} h^{-1} \msun$, consistent with halo masses estimated from the observed Tully-Fisher relation. In the simulations, the halo virial mass scales with galaxy luminosity as $L^{1.5}$ for central galaxies in halos and for galaxies in low-density regions. For satellite galaxies, and for galaxies in high-density regions, there is no simple relation between galaxy luminosity and halo mass and care must be exercised when interpreting the lensing results.

We have derived star formation rates (SFRs), gas-phase oxygen abundances and effective dust absorption optical depths for a sample of galaxies drawn from the Stromlo-APM redshift survey using the new Charlot and Longhetti (2001; CL01) models, which provide a physically consistent description of the effects of stars, gas and dust on the integrated spectra of galaxies. Our sample consists of 705 galaxies with measurements of the fluxes and equivalent widths of Halpha, [OII], and one or both of [NII] and [SII]. For a subset of the galaxies, 60 and 100 micron IRAS fluxes are available. We compare the star formation rates derived using the models with those derived using standard estimators based on the Halpha, the [OII] and the far-infrared luminosities of the galaxies. The CL01 SFR estimates agree well with those derived from the IRAS fluxes, but are typically a factor of ~3 higher than those derived from the Halpha or the [OII] fluxes, even after the usual mean attenuation correction of A_Halpha=1 mag is applied to the data. We show that the reason for this discrepancy is that the standard Halpha estimator neglects the absorption of ionizing photons by dust in HII regions and the contamination of Halpha emission by stellar absorption. We also use our sample to study variations in star formation and metallicity as a function of galaxy absolute bJ magnitude. For this sample, the star formation rate per unit bJ luminosity is independent of magnitude. The gas-phase oxygen abundance does increase with bJ luminosity, although the scatter in metallicity at fixed magnitude is large.

We simulate the formation and evolution of the local galaxy population starting from initial conditions with a smoothed linear density field which matches that derived from the IRAS 1.2 Jy galaxy survey. Our simulations track the formation and evolution of all dark matter haloes more massive than 10e+11 solar masses out to a distance of 8000 km/s from the Milky Way. We implement prescriptions similar to those of Kauffmann et al. (1999) to follow the assembly and evolution of the galaxies within these haloes. We focus on two variants of the CDM cosmology: an LCDM and a tCDM model. Galaxy formation in each is adjusted to reproduce the I-band Tully-Fisher relation of Giovanelli et al. (1997). We compare the present-day luminosity functions, colours, morphology and spatial distribution of our simulated galaxies with those of the real local population, in particular with the Updated Zwicky Catalog, with the IRAS PSCz redshift survey, and with individual local clusters such as Coma, Virgo and Perseus. We also use the simulations to study the clustering bias between the dark matter and galaxies of differing type. Although some significant discrepancies remain, our simulations recover the observed intrinsic properties and the observed spatial distribution of local galaxies reasonably well. They can thus be used to calibrate methods which use the observed local galaxy population to estimate the cosmic density parameter or to draw conclusions about the mechanisms of galaxy formation. To facilitate such work, we publically release our z=0 galaxy catalogues, together with the underlying mass distribution.

We study the cross-correlation between quasars and galaxies by embedding models for the formation and evolution of the two populations in cosmological N-body simulations. We adopt the quasar evolution model of Kauffmann & Haehnelt (2000), in which supermassive black holes are formed and fuelled during major mergers. We define the `bias' parameter b_QG as the ratio of the cross-correlation function xi_QG to the galaxy auto-correlation function xi_GG. On scales larger than 1 h^{-1} Mpc, the values of b_QG predicted by our models at low redshift depend very little on galaxy selection. They measure the characteristic mass of the dark matter halos that host quasars and can be used to estimate the typical quasar lifetime. In current redshift surveys, such measurements will constrain the lifetimes of low z quasars more accurately than measurements of the quasar auto-correlation function, because galaxies have much higher space densities than quasars. On scales smaller than 1 h^{-1} Mpc, the main contribution to xi_QG comes from quasar/galaxy pairs in the same dark matter halo. The amplitude of xi_QG depends both on the location of the host galaxy and on the density profile of other galaxies within the halo. As a result, measurements on these scales yield information about the processes responsible for fuelling super- massive black holes. At high redshifts, our models predict that quasars of fixed luminosity are located in less massive halos than at low redshift. They are therefore less biased relative to galaxies of given luminosity or stellar mass. We have used the simulations to calculate the evolution of the quasar auto- correlation function. We find that models with quasar lifetimes in the range 10^6-10^7 years provide a good match to the results of the 2dF QSO survey.

We explore a model in which efficient star formation in galaxies is triggered by merging satellites. We show that the merger/interaction rates predicted by hierarchical galaxy formation models depend strongly on galaxy mass. If a merger of a satellite larger than 1% the mass of the primary triggers efficient star formation, low-mass dwarf galaxies experience strong bursts separated by quiescent periods lasting several Gyr. Massive galaxies, such as our own Milky Way, are perturbed by a 1% satellite every few hundred million years, and thus have fluctuating, but relatively continuous star formation histories. We study the spectral signatures of a population of galaxies undergoing intermittent star formation by combining the models with the latest version of the Bruzual & Charlot spectral synthesis code. We concentrate on spectral indicators that are sensitive primarily to stellar age, rather than to metallicity and extinction. We show that if a population undergoes intermittent rather than continuous star formation, the signatures of the bursts should be evident from the observed dispersion in the star formation rates, H-delta equivalent widths and gas mass fractions of galaxies with 4000 Angstrom break strengths indicative of recent or ongoing star formation.

We simulate the assembly of a massive rich cluster and the formation of its constituent galaxies in a flat, low-density universe. Our most accurate model follows the collapse, the star-formation history and the orbital motion of all galaxies more luminous than the Fornax dwarf spheroidal, while dark halo structure is tracked consistently throughout the cluster for all galaxies more luminous than the SMC. Within its virial radius this model contains about 2.0e7 dark matter particles and almost 5000 distinct dynamically resolved galaxies. Simulations of this same cluster at a variety of resolutions allow us to check explicitly for numerical convergence both of the dark matter structures produced by our new parallel N-body and substructure identification codes, and of the galaxy populations produced by the phenomenological models we use to follow cooling, star formation, feedback and stellar aging. This baryonic modelling is tuned so that our simulations reproduce the observed properties of isolated spirals outside clusters. Without further parameter adjustment our simulations then produce a luminosity function, a mass-to-light ratio, luminosity, number and velocity dispersion profiles, and a morphology-radius relation which are similar to those observed in real clusters. In particular, since our simulations follow galaxy merging explicitly, we can demonstrate that it accounts quantitatively for the observed cluster population of bulges and elliptical galaxies.

Recent work has demonstrated that there is a tight correlation between the mass of a black hole and the velocity dispersion of the bulge of its host galaxy. We show that the model of Kauffmann & Haehnelt, in which bulges and supermassive black holes both form during major mergers, produces a correlation between M_bh and sigma with slope and scatter comparable to the observed relation. In the model, the M_bh - sigma relation is significantly tighter than the correlation between black hole mass and bulge luminosity or the correlation between bulge luminosity and velocity dispersion. There are two reasons for this: i) the gas masses of bulge progenitors depend on the velocity dispersion but not on the formation epoch of the bulge, whereas the stellar masses of the progenitors depend on both; ii) mergers between galaxies move black holes along the observed M_bh - sigma relation, even at late times when the galaxies are gas-poor and black holes grow mainly by merging of pre-existing black holes. We conclude that the small scatter in the observed M_bh - sigma relation is consistent with a picture in which bulges and black holes form over a wide range in redshift.

We combine dissipationless N-body simulations and semi-analytic models of galaxy formation to study the spatial and kinematic distributions of cluster galaxies in a LCDM cosmology. We investigate how the star formation rates, colours and morphologies of galaxies vary as a function of distance from the cluster centre and compare our results with the CNOC1 survey of galaxies from 15 X-ray luminous clusters in the redshift range 0.18 to 0.55. In our model, gas no longer cools onto galaxies after they fall into the cluster and their star formation rates decline on timescales of 1-2 Gyr. Galaxies in cluster cores have lower star formation rates and redder colours than galaxies in the outer regions because they were accreted earlier. Our colour and star formation gradients agree with those those derived from the data. The difference in velocity dispersions between red and blue galaxies observed in the CNOC1 clusters is also well reproduced by the model. We assume that the morphologies of cluster galaxies are determined solely by their merging histories. Morphology gradients in clusters arise naturally, with the fraction of bulge- dominated galaxies highest in cluster cores. We compare these gradients with the CNOC1 data and find excellent agreement for bulge-dominated galaxies. The simulated clusters contain too few galaxies of intermediate bulge-to-disk ratio, suggesting that additional processes may influence the morphological evolution of disk-dominated galaxies in clusters. Although the properties of the cluster galaxies in our model agree extremely well with the data, the same is not true of field galaxies. Both the star formation rates and the colours of bright field galaxies appear to evolve much more strongly from redshift 0.2 to 0.4 in the CNOC1 field sample than in our simulations.

We present new near-infrared J and K imaging data for 67 galaxies from the Universidad Complutense de Madrid survey used to determine the SFR density of the local universe by Gallego et al. (1995). This is a sample of local star-forming galaxies with redshift lower than 0.045, and they constitute a representative subsample of the galaxies in the complete UCM survey. From the new data, complemented with our own Gunn-r images and long-slit optical spectroscopy, we have measured integrated K-band luminosities, r-J and J-K colours, and Halpha luminosities and equivalent widths. Using a maximum likelihood estimator and a complete set of evolutionary synthesis models, these observations have allowed us to estimate the strength of the current (or most recent) burst of star formation, its age, the star-formation rate and the total stellar mass of the galaxies. An average galaxy in the sample has a stellar mass of 5E10 Msun and is undergoing (or recently completed) a burst of star formation involving about 2 per cent of its total stellar mass. We have identified two separate classes of star-forming galaxies in the UCM sample: low luminosity, high excitation galaxies (HII-like) and relatively luminous spirals galaxies (starburst disk-like). The former show higher specific star formation rates (SFR per unit mass) and burst strengths, and lower stellar masses than the latter. With regard to their specific star formation rates, the UCM galaxies are intermediate objects between normal quiescent spirals and the most extreme HII galaxies.

We study the biasing relation between dark-matter halos or galaxies and the underlying mass distribution, using cosmological $N$-body simulations in which galaxies are modelled via semi-analytic recipes. The nonlinear, stochastic biasing is quantified in terms of the mean biasing function and the scatter about it as a function of time, scale and object properties. The biasing of galaxies and halos shows a general similarity and a characteristic shape, with no galaxies in deep voids and a steep slope in moderately underdense regions. At $\sim 8\hmpc$, the nonlinearity is typically $\lsim 10$ percent and the stochasticity is a few tens of percent, corresponding to $\sim 30$ percent variations in the cosmological parameter $\beta=\Omega^{0.6}/b$. Biasing depends weakly on halo mass, galaxy luminosity, and scale. The time evolution is rapid, with the mean biasing larger by a factor of a few at $z\sim 3$ compared to $z=0$, and with a minimum for the nonlinearity and stochasticity at an intermediate redshift. Biasing today is a weak function of the cosmological model, reflecting the weak dependence on the power-spectrum shape, but the time evolution is more cosmology-dependent, relecting the effect of the growth rate. We provide predictions for the relative biasing of galaxies of different type and color, to be compared with upcoming large redshift surveys. Analytic models in which the number of objects is conserved underestimate the evolution of biasing, while models that explicitly account for merging provide a good description of the biasing of halos and its evolution, suggesting that merging is a crucial element in the evolution of biasing.

We discuss constraints on the assembly history of supermassive black holes from the observed remnant black holes in nearby galaxies and from the emission caused by accretion onto these black holes. We also summarize the results of a specific model for the evolution of galaxies and their central black holes which traces their hierachical build-up in CDM-like cosmogonies. The model assumes (i) that black holes, ellipticals and starburts form during major mergers of galaxies (ii) that the gas fraction in galaxies decreases with decreasing redshift (iii) that the optical bright phase of a QSO lasts for about 10^7 years. The model succesfully reproduces the evolution of cold gas as traced by damped damped Lyman alpha systems, the evolution of optically bright QSOs, the remnant black hole mass distribution and the host galaxy luminosities of QSOs.

We incorporate a simple scheme for the growth of supermassive black holes into semi-analytic models that follow the formation and evolution of galaxies in a cold dark matter dominated Universe. We assume that supermassive black holes are formed and fuelled during major mergers. If two galaxies of comparable mass merge, their central black holes coalesce and a few percent of the gas in the merger remnant is accreted by the new black hole over a timescale of a few times 10^7 years. With these simple assumptions, our model not only fits many aspects of the observed evolution of galaxies, but also reproduces quantitatively the observed relation between bulge luminosity and black hole mass in nearby galaxies, the strong evolution of the quasar population with redshift and the relation between the luminosities of nearby quasars and those of their host galaxies. The strong decline in the number density of quasars from z=2 to z=0 is due to the combination of three effects: i) a decrease in the merging rate, ii) a decrease in the amount of cold gas available to fuel black holes, and iii) an increase in the timescale for gas accretion. In a LCDM cosmology the predicted decline in the total content of cold gas in galaxies is consistent with that inferred from observations of damped Lyman-alpha systems. Our results strongly suggest that the evolution of supermassive black holes, quasars and starbursts is inextricably linked to the hierarchical build-up of galaxies.

In models of galaxy formation in a hierarchical Universe, elliptical galaxies form through the merging of smaller disk systems. These models yield a number of testable predictions if reliable techniques for determining the relative ages and compositions of the stellar populations of different galaxies can be found: 1) ellipticals in low-density environments form later than ellipticals in clusters, 2) more massive ellipticals form later, 3) more massive ellipticals form in dissipationless mergers from disk galaxies with low gas content. While colours and the Balmer line strengths of galaxies can be used to infer the average ages of the stellar populations of ellipticals, alpha/Fe element ratios carry information about the timescale over which star formation took place. Here we present preliminary results from semi-analytic models for the distribution of Mg/Fe ratios in galaxies as a function of morphological type, luminosity and environment.

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 Omega_0=1 (tauCDM) or Omega_0=0.3 and Lambda=0.7 (LambdaCDM). 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 we study the relation of these groups to real virialised systems. Most features of CfA2N groups are reproduced quite well by both models with no obvious dependence on Omega_0. 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.

In hierarchical models of galaxy formation, ellipticals form from the merging of disk galaxies drawn together by gravity as their surrounding dark halos coalesce. Using semi-analytic techniques, we are able to follow the merging, star formation and chemical enrichment histories of galaxies for a given choice of cosmological parameters. We calculate the ages and metallicities of elliptical galaxies in a standard CDM cosmology and show that we can reproduce the observed correlation between magnitude and colour, and the scatter about this relation. We also study the evolution of ellipticals to high redshift. In rich clusters, the models predict a substantial population of "evolved" ellipticals at z>1. However, the average number density of bright ellipticals must decrease at high redshift. The rate of decrease is sensitive to cosmological parameters. In a high-density model, the abundance of bright ellipticals declines by factor of 3 by z=1. In a low-density model (Omega=0.3, Lambda=0.7), the abundance drops by 30% at z=1, but is down by a factor 3 at z=2

In hierarchical cosmologies the evolution of galaxy clustering depends both on cosmological quantities such as Omega and Lambda, which determine how dark matter halos form and evolve, and on the physical processes - cooling, star formation and feedback - which drive the formation of galaxies within these merging halos. In this paper, we combine dissipationless cosmological N-body simulations and semi-analytic models of galaxy formation in order to study how these two aspects interact. We focus on the differences in clustering predicted for galaxies of differing luminosity, colour, morphology and star formation rate and on what these differences can teach us about the galaxy formation process. We show that a "dip" in the amplitude of galaxy correlations between z=0 and z=1 can be an important diagnostic. Such a dip occurs in low-density CDM models because structure forms early and dark matter halos of 10**12 solar masses, containing galaxies with luminosities around L*, are unbiased tracers of the dark matter over this redshift range; their clustering amplitude then evolves similarly to that of the dark matter. At higher redshifts bright galaxies become strongly biased and the clustering amplitude increases again. In high density models, structure forms late and bias evolves more rapidly. As a result, the clustering amplitude of L* galaxies remains constant from z=0 to 1. The strength of these effects is sensitive to sample selection. The dip is weaker for galaxies with lower star formation rates, redder colours, higher luminosities and earlier morphological types. We explain why this is the case and how it is related to the variation with redshift of the abundance and and environment of the observed galaxies. Studies of clustering evolution as a function of galaxy properties should place strong constraints on models.

We introduce a new technique for following the formation and evolution of galaxies in cosmological N-body simulations. Dissipationless simulations are used to track the formation and merging of dark matter halos as a function of redshift. Simple prescriptions, taken directly from semi-analytic models of galaxy formation, are adopted for cooling, star formation, supernova feedback and the merging of galaxies within the halos. This scheme enables us to study the clustering properties of galaxies and to investigate how selection by type, colour or luminosity influences the results. In this paper, we study properties of the galaxy distribution at z=0. These include luminosity functions, colours, correlation functions, pairwise peculiar velocities, cluster M/L ratios and star formation rates. We focus on two variants of a CDM cosmology: a high- density model with Gamma=0.21 (TCDM) and a low-density model with Omega=0.3 and Lambda=0.7 (LCDM). Both are normalized to reproduce the I-band Tully-Fisher relation near a circular velocity of 220 km/s. Our results depend strongly both on this normalization and on the prescriptions for star formation and feedback. Very different assumptions are required to obtain an acceptable model in the two cases. For TCDM, efficient feedback is required to suppress the growth of galaxies low-mass field halos. Without it, there are too many galaxies and the correlation function turns over below 1 Mpc. For LCDM, feedback must be weak, otherwise too few L* galaxies are produced and the correlation function is too steep. Given the uncertainties in modelling some of the key physical processes, we conclude that it is not yet possible to draw conclusions about the values of cosmological parameters from studies of this kind. Further work on global star formation and feedback effects is required to narrow the range of possibilities

There are two major approaches to modelling galaxy evolution. The traditional view is that the most massive galaxies were assembled early and have evolved with steeply declining star formation rates since a redshift of 2 or higher. According to hierarchical theories, massive galaxies were assembled much more recently from mergers of smaller subunits. Here we present a simple observational test designed to differentiate between the two. The observed K-band flux from a galaxy is a good measure of its stellar mass even at high redshift. It is probably only weakly affected by dust extinction. We compute the evolution of the observed K-band luminosity function for traditional, pure luminosity evolution (PLE) models and for hierarchical models. At z=0, both models can fit the observed local K-band luminosity function. By redshift 1, they differ greatly in the predicted abundance of bright galaxies. We calculate the redshift distributions of K-band selected galaxies and compare these with available data. We show that the number of K<19 galaxies with redshifts greater than 1 is well below the numbers predicted by the PLE models. In the Songaila et al (1994) redshift sample of 118 galaxies with 16<K<18, 33 galaxies are predicted to lie at z>1. Only 2 are observed. In the Cowie et al. (1996) redshift sample of 52 galaxies with 18<K<19, 28 galaxies are predicted to lie at at z>1. Only 5 are observed. Both these samples are more than 90% complete. We conclude that there is already strong evidence that the abundance of massive galaxies at z=1 is well below the local value. This is inconsistent with the traditional model, but similar to the expectations of hierarchical models.

The K-band Hubble diagram for a sample of brightest cluster galaxies (BCGs) in the redshift range 0<z<1 shows a very small scatter (0.3 magnitudes r.m.s). The BCGs exhibit very little luminosity evolution in this redshift range: if q_0=0.0 we detect no luminosity evolution; for q_0=0.5 we measure a small negative evolution. 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. This suggests that the stellar mass of the BCGs has been assembled over time through merging and accretion. We estimate that the stellar mass in a typical BCG has grown by a factor ~ 2 since z=1 if q_0=0.0 or by factor ~ 4 if q_0=0.5. These results are in remarkably 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.

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 magnitudes). The BCGs exhibit very little luminosity evolution in this redshift range: if q_0=0.0 we detect no luminosity evolution; for q_0=0.5 we measure a small negative evolution (i.e., BCGs were about 0.5 magnitudes 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_form >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.0 or by factor ~ 4 if q_0=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.

We use large N-body simulations of dissipationless gravitational clustering in cold dark matter (CDM) cosmologies to study whether the properties of dark matter halos are affected by their environment. We look for correlations between the masses, formation redshifts, concentrations, shapes and spins of halos and the overdensity of their local environment. We also look for correlations of these quantities with the local tidal field. Our conclusion is extremely simple. Only the mass distribution varies as a function of environment. This variation is well described by a simple analytic formula based on the conditional Press-Schechter theory. We find no significant dependence of any other halo property on environment. Our results do not depend on our choice of cosmology. According to current hierarchical models, the structure and evolutionary history of a galaxy is fully determined by the structure and evolutionary history of the dark halo in which it is embedded. If these models are correct, clustering variations between galaxies of differing morphological types, luminosities, colours and surface brightnesses, must all arise because the halo mass function is skewed towards high mass objects in overdense regions of the Universe and towards low mass objects in underdense regions.

In this paper, we present a model of the formation and chemical enrichment of elliptical galaxies that differs from the conventional picture in two ways: 1)Ellipticals do not form in a single monolithic collapse and burst of star formation at high redshift. Instead, most of their stars form at modest rates in disk galaxies, which then merge to form the ellipticals. 2)Galaxies do not undergo closed-box chemical evolution. Instead, metals can be transferred between the stars, cold gas and the hot gas halos of the galaxies. It is assumed that metals are ejected out of disk galaxies during supernova explosions and these metals enter the hot gas component. The fact that metals are more easily ejected from small galaxies leads to the establishment of a mass-metallicity relation for the disk systems. Big ellipticals are more metal rich because they are formed from the mergers of bigger disks. We use semi-analytic techniques to follow the formation, evolution, and chemical enrichment of cluster ellipticals in a merging hierarchy of dark matter halos. The inclusion of the new metallicity-dependent spectral synthesis models of Bruzual & Charlot allow us to compute the colours, line indices and mass-to- light ratios of these galaxies. If feedback is assumed to be efficient, we are able to reproduce the slope and scatter of the colour-magnitude and the Mg2-sigma relations. We also study the evolution of these relations to high redshift. We show that the luminosity-metallicity relation remains fixed, but the mean stellar ages of the galaxies scale with the age of the Universe. Finally, we study the enrichment history of the intracluster gas. We find that 80% of the metals were ejected by low-mass galaxies at redshifts greater than 1. The metallicity of the ICM thus evolves very little out to z > 1.

The standard picture holds that giant elliptical galaxies formed in a single burst at high redshift. Aging of their stellar populations subsequently caused them to fade and become redder. The Canada-France Redshift Survey provides a sample of about 125 galaxies with the luminosities and colours of passively evolving giant ellipticals and with 0.1 < z < 1. This sample is inconsistent with the standard evolutionary picture with better than 99.9% confidence. The standard Schmidt test gives <V/V_{max}> = 0.398 when restricted to objects with no detected star formation, and <V/V_{max}> = 0.410 when objects with emission lines are also included. A smaller sample of early-type galaxies selected from tha Hawaii Deep Survey gives equally significant results. With increasing redshift a larger and larger fraction of the nearby elliptical and S0 population must drop out of of the sample, either because the galaxies are no longer single units or because star formation alters their colours. If the remaining fraction is modelled as F=(1+z)^{-\gamma}, the data imply gamma = 1.5 +- 0.3. At z=1 only about one third of bright E and S0 galaxies were already assembled and had the colours of old passively evolving systems. We discuss the sensitivity of these results to the incompleteness corrections and stellar population models we have adopted. We conclude that neither is uncertain enough to reconcile the observations with the standard picture. Hierarchical galaxy formation models suggest that both merging and recent star formation play a role in the strong evolution we have detected.

The analysis of disk formation is based on the White & Rees (1978) picture, in which disk galaxies form by continuous cooling and accretion of gas within a merging hierarchy of dark matter halos. A simple Kennicutt law of star formation for disks, based on a single- fluid gravitational instability model is introduced. Since the gas supply in the disk is regulated by infall from the surrounding halo, the gas is always maintained at a critical threshold surface density. Chemical enrichment of the disks occurs when the surrounding hot halo gas is enriched with heavy elements ejected during surpernova explosions. This gas then cools onto the disk producing a new generation of metal-rich stars. I first show that models of this type can reproduce many of the observed properties of a typical spiral galaxy like the Milky Way, including its gas and stellar surface density profiles and the observed relationship between the ages and metallicities of solar neighbourhood stars. I then use the models to make inferences about the properties of disk galaxies at high redshift. The total neutral hydrogen density Omega(HI) increases at higher z. The predicted increase is mild, but is roughly consistent with the latest derivation of Omega(HI) as a function of z by Storrie-Lombardi & MacMahon (1995). The models are also able to account for some of the other trends other trends seen in the high-redshift data, including the increase in the number of high column-density systems at high redshift, as well as the metallicity distribution of damped Lyman-alpha systems at redshifts 2-3.

We outline a simple approach to understanding the physical origin of bias in the distribution of galaxies relative to that of dark matter. The first step is to specify how collapsed, virialized halos of dark matter trace the overall matter distribution. The next step is to make a connection between halos and the luminous galaxies we observe. We appeal to the results of semi-analytic models of galaxy formation that are tuned to fit the observed luminosity functions of local groups and clusters. We have also used a high-resolution N-body simulation of a cold dark matter (CDM) universe to study the bias relation in more detail. The differences between the galaxy and dark matter distributions are quantified using a number of different clustering statistics. We arrive at the following general conclusions: 1) A comparison of the galaxy and dark matter density fields shows that linear biasing is a good description on large scales. 2) The bias factor b depends on the shape and normalization of the power spectrum. The lower the normalization, the larger the bias. More bias is obtained for spectra with more power on large scales. For "realistic" models, b ranges from 1 to 2.5. 3) Galaxies of different luminosity or morphology have different bias factors. 4) The scale dependence of the bias factor is weak.

The tightness of the observed colour-magnitude and Mg$_{2}$- velocity dispersion relations for elliptical galaxies has often been cited as an argument against a picture in which ellipticals form by the merging of spiral disks. A common view is that merging would mix together stars of disparate ages and produce a large scatter in these relations. Here I use semi-analytic models of galaxy formation to derive the distribution of the mean ages, colours and metallicities of the stars in elliptical galaxies formed by mergers in a flat CDM universe. It is seen that most of the stars in ellipticals form at relatively high redshift (z > 1.9) and that the predicted scatter in the colour-magnitude and Mg_2 - sigma relations falls within observational bounds. I conclude that the apparent homogeneity in the properties of the stellar populations of ellipticals is not inconsistent with a merger scenario for the origin of these systems.

We use the semi-analytic models of galaxy formation developed by Kauffmann, White \& Guiderdoni to generate predictions for the observed properties of cluster and group galaxies at redshifts between 0 and 0.6. We examine four sets of cosmological initial conditions: low-density CDM models with and without cosmological constant, a flat CDM model and a mixed dark matter model. These models were selected because they span a wide range in cluster formation epoch. The semi-analytic models that we employ are able to follow both the evolution of the dark matter component of clusters and the formation and evolution of the stellar populations of the cluster galaxies. We are thus able to generate model predictions that can be compared directly with the observational data. In the low-density CDM models, clusters form at high red- shift and accrete very little mass at recent times. Our models predict that essentially no evolution in the observed properties of clusters will have occurred by a redshift of 0.6, in direct contradiction with the data. In contrast, in the MDM model, both galaxies and clusters form extremely late. This model predicts evolution which appears to be too extreme to be in agreement with the observations. The flat CDM model, which is intermediate in structure formation epoch, is most successful. This model is able to account for the evolution of the blue fraction of rich clusters with redshift, the relationship between blue fraction and cluster richness at different epochs, and the changes in the distribution of the morphologies of cluster galaxies by a redshift of 0.4.

We show that the rapid evolution in the fraction of blue, star-forming galaxies seen in clusters as a function of redshift (the Butcher-Oemler effect) can be explained very simply if structure formation in the universe proceeeds hierarchically. We show that a rich cluster observed at high redshift has had a significantly different evolutionary history to a cluster of the same richness observed today. High redshift clusters take longer to assemble and thus under- go more merging at small lookback times. We have investigated two models of star formation in cluster galaxies: 1) a model in which star formation is induced by galaxy-galaxy mergers and interactions and 2) a model in which star formation is regulated by the infall of galaxies onto larger systems such as groups and clusters. Both models produce trends consistent with the Butcher- Oemler effect. Our models of cluster formation and evolution allow us to make predictions about trends in the observed properties of clusters with redshift. We find that there should be a correlation between the mass of the cluster or group and the strength of the observed Butcher-Oemler effect, with more massive systems exhibiting more evolution than less massive systems. We also predict that both the blue galaxy fraction and the incidence of interacting galaxies in rich clusters should rise continuously with redshift. Finally, we have explored the influence of cosmological parameters on our predictions for cluster evolution. We find that models in which structure formation occurs at very early epochs, such as low $\Omega$ models, predict rather little recent star formation and merging activity in clusters at redshifts of around 0.4.

There is accumulating observational evidence suggesting that damped Ly$\alpha$ absorption systems systems are the progenitors of present-day spiral galaxies. We use the observed properties of these systems to place constraints on the history of star formation in galactic disks, and on cosmological theories of structure formation in the universe. We show that the observed increase in $\Omega_{HI}$ contributed by damped Ly$\alpha$ systems at high redshift implies that star formation must have been considerably less efficient in the past. We also show that the data can constrain cosmological models in which structure forms at late epochs. A mixed dark matter (MDM) model with $\Omega_{\nu}=0.3$ is unable to reproduce the mass densities of cold gas seen at high redshift, even in the absence of any star formation. We show that at redshifts greater than 3, this model predicts that the total baryonic mass contained in dark matter halos with circular velocities $V_c > 35$ km s$^{-1}$ is less than the observed mass of HI in damped systems. At these redshifts, the photo-ionizing background would prevent gas from dissipating and collapsing to form high column density systems in halos smaller than 35 km s$^{-1}$. MDM models are thus ruled out by the observations.