Source author record

Guinevere Kauffmann

Guinevere Kauffmann appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

83works
4topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

83 published item(s)

preprint2022arXiv

A panchromatic view of star cluster formation in a simulated dwarf galaxy starburst

We present a photometric analysis of star and star cluster (SC) formation in a high-resolution simulation of a dwarf galaxy starburst that allows the formation of individual stars to be followed. Previous work demonstrated that the properties of the SCs formed in the simulation are in good agreement with observations. In this paper, we create mock spectral energy distributions and broad-band photometric images using the radiative transfer code SKIRT 9. We test several observational star formation rate (SFR) tracers and find that $24$ $μ$m, total infrared and H$α$ trace the underlying SFR during the (post)starburst phase, while UV tracers yield a more accurate picture of star formation during quiescent phases prior to and after the merger. We then place the simulated galaxy at distances of $10$ and $50$ Mpc and use aperture photometry at Hubble Space Telescope resolution to analyse the simulated SC population. During the starburst phase, a hierarchically forming set of SCs leads inaccurate source separation because of crowding. This results in estimated SC mass function slopes that are up to $\sim0.3$ shallower than the true slope of $\sim-1.9$ to $-2$ found for the bound clusters identified from the particle data in the simulation. The masses of the largest clusters are overestimated by a factor of up to $2.9$ due to unresolved clusters within the apertures. The aperture-based analysis also produces a relation between cluster formation efficiency and SFR surface density that is slightly flatter than that recovered from bound clusters. The differences are strongest in quiescent SF environments.

preprint2022arXiv

A study of 1000 galaxies with unusually young and massive stars in the SDSS: a search for hidden black holes

We select 1076 galaxies with extinction-corrected Halpha equivalent widths too large to be explained with a Kroupa (2001) IMF, and compare these with a control sample of galaxies that is matched in stellar mass, redshift and 4000 AA break strength, but with normal Halpha equivalent widths. Our goal is to study how processes such as black hole growth and energetic feedback processes from massive stars differ between galaxies with extreme central Halpha emission and galaxies with normal young central stellar populations. The stellar mass distribution of Halpha excess galaxies is peaked at 3 \times 10^10 Msun and almost all fall well within the star-forming locus in the [OIII]/Hbeta versus [NII]/Halpha BPT disgram. Halpha excess galaxies are twice as likely to exhibit Halpha line asymmetries and 1.55 times more likely to be detected at 1 GHz in the VLA FIRST survey compared to control sample galaxies. The radio luminosity per unit stellar mass decreases with the stellar age of the system. Using stacked spectra, we demonstrate that [NeV] emission is not present in the very youngest of the radio-quiet Halpha excess galaxies with detectable Wolf-Rayet features, suggesting that black hole growth has not yet commenced in such systems. [NeV] emission is detected in Halpha excess galaxies with radio detections and the strength of the line correlates with the radio luminosity. This is the clearest indication for a possible population of black holes that may be forming in a subset of the Halpha excess population.

preprint2022arXiv

Cool circumgalactic gas in galaxy clusters: connecting the DESI legacy imaging survey and SDSS DR16 MgII absorbers

We investigate the cool gas absorption in galaxy clusters by cross-correlating MgII absorbers detected in quasar spectra from Data Release 16 of the Sloan Digital Sky Survey (SDSS) with galaxy clusters identified in the Dark Energy Spectroscopic Instrument (DESI) survey. We find significant covering fractions ($1-5\, \%$ within $r_{500}$, depending on the chosen redshift interval), $\sim 4-5$ times higher than around random sightlines. While the covering fraction of cool gas in clusters decreases with increasing mass of the central galaxy, the total MgII mass within $r_{\rm 500}$ is nonetheless $\sim 10$ times higher than for SDSS luminous red galaxies (LRGs). The MgII covering fraction versus impact parameter is well described by a power law in the inner regions and a exponential function at larger distances. The characteristic scale of the transition between these two regimes is smaller for large equivalent width absorbers. Cross-correlating MgII absorption with photo-z selected cluster member galaxies from DESI reveals a statistically significant connection. The median projected distance between MgII absorbers and the nearest cluster member is $\sim200$ kpc, compared to $\sim500$ kpc in random mocks with the same galaxy density profiles. We do not find a correlation between MgII strength and the star formation rate of the closest cluster neighbour. This suggests that cool gas in clusters, as traced by MgII absorption, is: (i) associated with satellite galaxies, (ii) dominated by cold gas clouds in the intracluster medium, rather than by the interstellar medium of galaxies, and (iii) may originate in part from gas stripped from these cluster satellites in the past.

preprint2022arXiv

SDSS-IV MaNGA: Exploring the local scaling relations for N/O

We present, for the first time, the relationship between local stellar mass surface density, $\mathrm{Σ_{*}}$, and N/O derived from SDSS-IV MaNGA data, using a sample of $792765$ high signal-to-noise ratio star-forming spaxels. Using a combination of phenomenological modelling and partial correlation analysis, we find that $\mathrm{Σ_{*}}$ alone is insufficient to predict the N/O in MaNGA spaxels, and that there is an additional dependence on the local star formation rate surface density, $\mathrm{Σ_{SFR}}$. This effect is a factor of $3$ stronger than the dependence of 12+log(O/H) on $\mathrm{Σ_{SFR}}$. Surprisingly, we find that the local N/O scaling relations also depend on the total galaxy stellar mass at fixed $Σ_{*}$ as well as the galaxy size at fixed stellar mass. We find that more compact galaxies are more nitrogen rich, even when $\mathrm{Σ_{*}}$ and $\mathrm{Σ_{SFR}}$ are controlled for. We show that $\sim50\%$ of the variance of N/O is explained by the total stellar mass and size. Thus, the evolution of nitrogen in galaxies is set by more than just local effects and does not simply track the build up of oxygen in galaxies. The precise form of the N/O-O/H relation is therefore sensitive to the sample of galaxies from which it is derived. This result casts doubt on the universal applicability of nitrogen-based strong-line metallicity indicators derived in the local universe.

preprint2020arXiv

Comparing galaxy formation in the L-GALAXIES semi-analytical model and the IllustrisTNG simulations

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.

preprint2020arXiv

L-GALAXIES 2020: Spatially resolved cold gas phases, star formation and chemical enrichment in galactic discs

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.

preprint2020arXiv

Lyman-αabsorption beyond the disk of simulated spiral galaxies

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-$α$ 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-$α$ 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.

preprint2016arXiv

On the stellar halo metallicity profile of Milky Way-like galaxies in the Auriga simulations

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.

preprint2016arXiv

Resolution-independent modeling of environmental effects in semi-analytic models of galaxy formation that include ram-pressure stripping of both hot and cold gas

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.

preprint2016arXiv

The Data Reduction Pipeline for the SDSS-IV MaNGA IFU Galaxy Survey

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.

preprint2016arXiv

The morphology and kinematics of neutral hydrogen in the vicinity of z=0 galaxies with Milky Way masses -- a study with the Illustris simulation

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.

preprint2016arXiv

The Properties of the Circumgalactic Medium in Red and Blue Galaxies: Results from the COS-GASS+COS-Halos Surveys

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~$α$ 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~$α$ 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~$α$ lines. The ratio of the velocity offset to the velocity dispersion for the Lyman~$α$ 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~$α$ 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.

preprint2015arXiv

A systematic study of the inner rotation curves of galaxies observed as part of the GASS and COLD GASS surveys

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.

preprint2015arXiv

Connection Between the Circumgalactic Medium and the Interstellar Medium of Galaxies: Results from the COS-GASS Survey

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$α$ absorption-lines) in the CGM of 92% of the galaxies. We find the radial profile of the CGM as traced by the Ly$α$ 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$α$ 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$α$ equivalent width. This is stronger than the analogous correlation between corrected Ly$α$ 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.

preprint2015arXiv

HI scaling relations of galaxies in the environment of HI-rich and control galaxies observed by the Bluedisk project

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.

preprint2015arXiv

Physical origin of the large-scale conformity in the specific star formation rates of galaxies

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.

preprint2015arXiv

The distribution of atomic hydrogen in EAGLE galaxies: morphologies, profiles, and HI holes

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 Σ_HI > 1 M_Sun pc^-2 (by a factor of <~ 2 compared to data from the Bluedisk survey). In the outer region (Σ_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)

preprint2015arXiv

The Outer Stellar Populations and Environments of Unusually HI-rich Galaxies

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.

preprint2015arXiv

The spatially resolved Kennicutt-Schmidt relation in the HI dominated regions of spiral and dwarf irregular galaxies

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 $Σ_{SFR}$ at a given value of $Σ_{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.

preprint2015arXiv

The variation in molecular gas depletion time among nearby galaxies: II the impact of galaxy internal structures

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.

preprint2015arXiv

The Variation in Molecular Gas Depletion Time among Nearby Galaxies: What are the Main Parameter Dependencies?

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.

preprint2014arXiv

An observational and theoretical view of the radial distribution of HI gas in 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 $Λ$ 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.

preprint2014arXiv

Overview of the SDSS-IV MaNGA Survey: Mapping Nearby Galaxies at Apache Point Observatory

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.

preprint2014arXiv

Quantitative constraints on starburst cycles in galaxies with stellar masses in the range 10**8-10**10 Msol

We have used 4000 Å break and H$δ_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.

preprint2014arXiv

The nature of obscuration in AGN: II. insights from clustering properties

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.

preprint2013arXiv

Dilution in elliptical galaxies: Implications for the relation between metallicity, stellar mass and star formation rate

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.

preprint2013arXiv

Estimating gas masses and dust-to-gas ratios from optical spectroscopy

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.

preprint2013arXiv

Galactic accretion and the outer structure of galaxies in the CDM 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

preprint2013arXiv

Modelling Element Abundances in Semi-analytic Models of Galaxy Formation

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.

preprint2013arXiv

On the Evolution of the Bi-Modal Distribution of Galaxies in the Plane of Specific Star Formation Rate versus Stellar Mass

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.

preprint2013arXiv

Spatially resolved star formation histories of nearby galaxies: evidence for episodic star formation in discs

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.

preprint2013arXiv

The Bluedisks project, a study of unusually HI-rich galaxies: I. HI Sizes and Morphology

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 $Δ$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.

preprint2013arXiv

The Nature of Obscuration in AGN: I. Insights from Host Galaxies

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.

preprint2012arXiv

A Re-examination of Galactic Conformity and a Comparison with Semi-analytic Models of Galaxy Formation

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.

preprint2012arXiv

ALFALFA HI Data Stacking III. Comparison of environmental trends in HI gas mass fraction and specific star formation rate

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.

preprint2012arXiv

COLD GASS, an IRAM Legacy Survey of Molecular Gas in Massive Galaxies: III. Comparison with semi-analytic models of galaxy formation

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.

preprint2012arXiv

Evolution of the Most Massive Galaxies to z ~ 0.6: II. The link between radio AGN activity and star formation

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.

preprint2012arXiv

Gas depletion in cluster galaxies depends strongly on their internal structure

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 $μ_*$ 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.

preprint2012arXiv

Quantifying the role of bars in the build-up of central mass concentrations in disk 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.

preprint2012arXiv

The Baryon Oscillation Spectroscopic Survey of SDSS-III

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.

preprint2012arXiv

The clustering of galaxies as a function of their photometrically-estimated atomic gas content

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.

preprint2012arXiv

The Effect of Star Formation on the Redshift Evolution of Interstellar Metals, Atomic and Molecular Gas in Galaxies

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.

preprint2012arXiv

The GALEX Arecibo SDSS Survey. VI. Second Data Release and Updated Gas Fraction Scaling Relations

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.

preprint2012arXiv

The Ninth Data Release of the Sloan Digital Sky Survey: First Spectroscopic Data from the SDSS-III Baryon Oscillation Spectroscopic Survey

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.

preprint2011arXiv

Autocorrelations of stellar light and mass at z~0 and ~1: From SDSS to DEEP2

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.

preprint2011arXiv

COLD GASS, an IRAM legacy survey of molecular gas in massive galaxies: I. Relations between H2, HI, stellar content and structural properties

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.

preprint2011arXiv

COLD GASS, an IRAM Legacy Survey of Molecular Gas in Massive Galaxies: II. The non-universality of the Molecular Gas Depletion Timescale

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.

preprint2011arXiv

Evolution of the Most Massive Galaxies to z=0.6: I. A New Method for Physical Parameter Estimation

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.

preprint2011arXiv

Extreme Feedback and the Epoch of Reionization: Clues in the Local Universe

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.

preprint2011arXiv

From dwarf spheroidals to cDs: Simulating the galaxy population in a LCDM cosmology

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.

preprint2011arXiv

SDSS-III: Massive Spectroscopic Surveys of the Distant Universe, the Milky Way Galaxy, and Extra-Solar Planetary Systems

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)

preprint2011arXiv

The GALEX Arecibo SDSS Survey V: The Relation between the HI Content of Galaxies and Metal Enrichment at their Outskirts

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.

preprint2011arXiv

The GALEX Arecibo SDSS Survey. IV. Baryonic Mass-Velocity-Size Relations of Massive Galaxies

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]

preprint2011arXiv

The relation between metallicity, stellar mass and star formation in galaxies: an analysis of observational and model data

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.

preprint2010arXiv

Absorption-line probes of the prevalence and properties of outflows in present-day star-forming galaxies

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 ($Σ_{\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 $Σ_{\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$Σ_{\rm SFR}$ and 1 dex in log M*).

preprint2010arXiv

ALFALFA HI Data Stacking I. Does the Bulge Quench Ongoing Star Formation in Early-Type Galaxies?

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).

preprint2010arXiv

Cluster Galaxies Die Hard

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.

preprint2010arXiv

Semi-Analytic Models on the Molecular and Neutral Gas in Galaxies

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.

preprint2010arXiv

The Accretion of Gas onto Galaxies as Traced by their Satellites

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.

preprint2010arXiv

The Atomic to Molecular Transition and its Relation to the Scaling Properties of Galaxy Disks in the Local Universe

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 $λ$. 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 $Σ_{SFR}\proptoΣ_{H_2}$ in the regime where the molecular gas dominates the total gas surface density, and $Σ_{SFR}\propto Σ_{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.

preprint2010arXiv

The GALEX Arecibo SDSS Survey II: The Star Formation Efficiency of Massive Galaxies

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]

preprint2010arXiv

The GALEX Arecibo SDSS survey: III. Evidence for the Inside-Out Formation of Galactic Disks

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.

preprint2010arXiv

UGC8802: A Massive Disk Galaxy in Formation

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.

preprint2009arXiv

Clustering of Radio Galaxies and Quasars

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.

preprint2009arXiv

Estimating the H I gas fractions of galaxies in the local Universe

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.

preprint2009arXiv

Feast and Famine: Regulation of Black Hole Growth in Low Redshift Galaxies

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.

preprint2009arXiv

Star formation and AGN activity in SDSS cluster galaxies

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.

preprint2009arXiv

The clustering of barred galaxies in the local Universe

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.

preprint2009arXiv

The GALEX Arecibo SDSS Survey. I. Gas Fraction Scaling Relations of Massive Galaxies and First Data Release

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]

preprint2009arXiv

The growth of supermassive black holes in pseudo-bulges, classical bulges and elliptical galaxies

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.

preprint2009arXiv

The UV-optical colours of brightest cluster galaxies in optically and X-ray selected clusters

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.

preprint2008arXiv

The Lopsidedness of Present-Day Galaxies: Connections to the Formation of Stars, the Chemical Evolution of Galaxies, and the Growth of Black Holes

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.

preprint2007arXiv

Bursty stellar populations and obscured AGN in galaxy bulges

[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.

preprint2005arXiv

Active Galactic Nuclei in the Sloan Digital Sky Survey: I. Sample Selection

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$α$ 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.

preprint2005arXiv

Active Galactic Nuclei in the Sloan Digital Sky Survey: II. Emission-Line Luminosity Function

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]$λ5007$ luminosity functions for Seyferts cover luminosity range of $10^{5-9}$$L_\odot$ in H$α$ 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.

preprint2005arXiv

The Colours of Elliptical galaxies

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 $α$-enhancement, while the optical-infrared colours are sensitive to metallicity and to $α$-enhancement, but are somewhat less sensitive to age.

preprint1996arXiv

Detection of Strong Evolution in the Population of Early-Type Galaxies

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)^{-γ}, 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.

preprint1995arXiv

Disk galaxies at Z=0 and at High Redshift: an explanation of the observed evolution of damped Ly-alpha systems

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.

preprint1995arXiv

Galaxy Formation and Large-Scale Bias

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.

preprint1995arXiv

The Ages of Elliptical Galaxies in a Merger Model

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.

preprint1994arXiv

Constraints on Models of Galaxy Formation from the Evolution of Damped Ly$α$ Absorption Systems

There is accumulating observational evidence suggesting that damped Ly$α$ 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 $Ω_{HI}$ contributed by damped Ly$α$ 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 $Ω_ν=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.

preprint1994arXiv

Hierarchical Clustering and the Butcher-Oemler Effect

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 $Ω$ models, predict rather little recent star formation and merging activity in clusters at redshifts of around 0.4.

preprint1994arXiv

The Observed Properties of High-Redshift Cluster Galaxies

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.