Source author record

Russell J. Smith

Russell J. Smith 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

29works
5topics
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

29 published item(s)

preprint2022arXiv

Comparing Lensing and Stellar Orbital Models of a Nearby Massive Strong-Lens Galaxy

Exploiting the relative proximity of the nearby strong-lens galaxy SNL-1, we present a critical comparison of the mass estimates derived from independent modelling techniques. We fit triaxial orbit-superposition dynamical models to spatially-resolved stellar kinematics, and compare to the constraints derived from lens modelling of high-resolution photometry. From the dynamical model, we measure the total (dynamical) mass enclosed within a projected aperture of radius the Einstein radius to be $\log_{10} M_{\mathrm{Ein.}} = 11.00 \pm 0.02$, which agrees with previous measurements from lens modelling to within $5\%$. We then explore the intrinsic (de-projected) properties of the best-fitting dynamical model. We find that SNL-1 has approximately-constant, intermediate triaxiality at all radii. It is oblate-like in the inner regions (around the Einstein radius) and tends towards spherical at larger radii. The stellar velocity ellipsoid gradually transforms from isotropic in the very central regions to radially-biased in the outskirts. We find that SNL-1 is dynamically consistent with the broader galaxy population, as measured by the relative fraction of orbit `temperatures' compared to the CALIFA survey. On the mass--size plane, SNL-1 occupies the most-compact edge given its mass, compared to both the MaNGA and SAMI surveys. Finally, we explore how the observed lensing configuration is affected by the orientation of the lens galaxy. We discuss the implications of such detailed models on future combined lensing and dynamical analyses.

preprint2020arXiv

A Triple Rollover: A third multiply-imaged source at z~6 behind the Jackpot gravitational lens

UUsing a five-hour adaptive-optics-assisted observation with MUSE, we have identified a doubly-imaged Ly alpha source at redshift 5.975 behind the z=0.222 lens galaxy J0946+1006 ('the Jackpot'). The source separation implies an Einstein radius of ~2.5 arcsec. Combined with the two previously-known Einstein rings in this lens (radii 1.4 arcsec at z = 0.609 and 2.1 arcsec at z ~ 2.4), this system is now a unique galaxy-scale triple-source-plane lens. We show that existing lensing models for J0946+1006 successfully map the two new observed images to a common point on the z=5.975 source plane. The new source will provide further constraints on the mass distribution in the lens and in the two previously known sources. The third source also probes two new distance scaling factors which are sensitive to the cosmological parameters of the Universe. We show that detection of a new multiply imaged emission-line source is not unexpected in observations of this depth; similar data for other known lenses should reveal a larger sample of multiple-image-plane systems for cosmography and other applications.

preprint2020arXiv

Discovering novae in early-type galaxies with MUSE: A chance find in NGC 1404, and twelve more candidates from an archival search

I report the discovery of a transient broad-H$α$ point source in the outskirts of the giant elliptical galaxy NGC 1404, discovered in archival observations taken with the MUSE integral field spectrograph. The H$α$ line width of 1950 km s$^{-1}$ FWHM, and luminosity of (4.1$\pm$0.1)$\times$10$^{36}$ erg s$^{-1}$, are consistent with a nova outburst, and the source is not visible in MUSE data obtained nine months later. A transient soft X-ray source was detected at the same position (within $<$1 arcsec), 14 years before the H$α$ transient. If the X-ray and H$α$ emission are from the same object, the source may be a short-timescale recurrent nova with a massive white dwarf accretor, and hence a possible Type-Ia supernova progenitor. Selecting broad-H$α$ point sources in MUSE archival observations for a set of nearby early-type galaxies, I discovered twelve more nova candidates with similar properties to the NGC 1404 source, including five in NGC 1380 and four in NGC 4365. Multi-epoch data are available for four of these twelve sources; all four are confirmed to be transient on $\sim$1 year timescales, supporting their identification as novae.

preprint2020arXiv

MNELLS: The MUSE Nearby Early-Type Galaxy Lens Locator Survey

Low-redshift strong-lensing galaxies can provide robust measurements of the stellar mass-to-light ratios in early-type galaxies (ETG), and hence constrain variations in the stellar initial mass function (IMF). At present, only a few such systems are known. Here, we report the first results from a blind search for gravitationally-lensed emission line sources behind 52 massive $z$ $<$ 0.07 ETGs with MUSE integral field spectroscopy. For 16 galaxies, new observations were acquired, whilst the other 36 were analysed from archival data. This project has previously yielded one confirmed galaxy-scale strong lens (J0403-0239) which we report in an earlier paper. J0403-0239 has since received follow-up observations, presented here, which indicate support for our earlier IMF results. Three cluster-scale, and hence dark-matter-dominated, lensing systems were also discovered (central galaxies of A4059, A2052 and AS555). For nine further galaxies, we detect a singly-imaged but closely-projected source within 6 arcsec (including one candidate with sources at three different redshifts); such cases can be exploited to derive upper limits on the IMF mass-excess factor, $α$. Combining the new lens and new upper limits, with the previously-discovered systems, we infer an average $\langle α\rangle$ = 1.06 $\pm$ 0.08 (marginalised over the intrinsic scatter), which is inconsistent with a Salpeter-like IMF ($α$ = 1.55) at the 6$σ$ level. We test the detection threshold in these short-exposure MUSE observations with the injection and recovery of simulated sources, and predict that one in twenty-five observations is expected to yield a new strong-lens system. Our observational results are consistent with this expected yield.

preprint2019arXiv

Subaru FOCAS IFU observations of two z=0.12 strong-lensing elliptical galaxies from SDSS MaNGA

We present new observations of two z=0.12 strong-lensing elliptical galaxies, originally discovered from the SDSS-IV MaNGA survey, using the new FOCAS IFU spectrograph on the Subaru Telescope. For J1436+4943, our observations confirm the identification of this system as a multiple-image lens, in a cusp configuration, with Einstein radius $θ_{Ein}$=2.0 arcsec. For J1701+3722, the improved data confirm earlier hints of a complex source plane, with different configurations evident in different emission lines. The new observations reveal a previously unseen inner counter-image to the [OIII] arc found from MaNGA, leading to a smaller revised Einstein radius of $θ_{Ein}$=1.6 arcsec. The inferred projected masses within the Einstein apertures (3.7-4.7kpc) are consistent with being dominated by stars with an initial mass function (IMF) similar to that of the Milky Way, and a dark matter contribution of ~35 per cent as supported from cosmological simulations. These results are consistent with `pure lensing' analyses of lower-redshift lenses, but contrast with claims for heavier IMFs from combined lensing-and-dynamical studies of more distant early-type galaxies.

preprint2016arXiv

A strong-lensing elliptical galaxy in the MaNGA survey

I report discovery of a new galaxy-scale gravitational lens system, identified using public data from the MaNGA survey, as part of a systematic search for lensed background line-emitters. The lens is SDSS J170124.01+372258.0, a giant elliptical galaxy with velocity dispersion $σ=256$ km/s, at a redshift of $z_l=0.122$. After modelling and subtracting the target galaxy light, the integral-field data-cube reveals [OII], [OIII] and H$β$ emission lines corresponding to a source at $z_s=0.791$, forming an identifiable ring around the galaxy center. The Einstein radius is $R_{Ein} \approx 2.3$ arcsec, projecting to ~5 kpc at the distance of the lens. The total projected lensing mass is $(3.6\pm0.6) \times 10^{11} M_\odot$, and the total J-band mass-to-light ratio is $3.0\pm0.7$ solar units. Plausible estimates of the likely dark matter content could reconcile this with a Milky-Way-like initial mass function (for which M/L~1.5 is expected), but heavier IMFs are by no means excluded with the present data. An alternative interpretation of the system, with a more complex source plane, is also discussed. The discovery of this system bodes well for future lens searches based on MaNGA and other integral-field spectroscopic surveys.

preprint2015arXiv

The IMF-sensitive 1.14-micron Na I doublet in early-type galaxies

We present J-band spectroscopy of passive galaxies focusing on the Na I doublet at 1.14 μm. Like the Na I 0.82 μm doublet, this feature is strong in low-mass stars and hence may provide a useful probe of the initial mass function (IMF). From high signal-to-noise composite spectra, we find that Na I 1.14 μm increases steeply with increasing velocity dispersion, σ, and for the most massive galaxies (σ> 300 km/s) is much stronger than predicted from synthetic spectra with Milky-Way-like IMFs and solar abundances. Reproducing Na I 1.14 μm at high σ likely requires either a very high [Na/H], or a bottom-heavy IMF, or a combination of both. Using the Na D line to break the degeneracy between IMF and abundance, we infer [Na/H] $\approx$ +0.5 and a steep IMF (single-slope-equivalent x $\approx$ 3.2, where x = 2.35 for Salpeter), for the high-σgalaxies. At lower mass (σ = 50-100 km/s), the line strengths are compatible with MW-like IMFs and near-solar [Na/H]. We highlight two galaxies in our sample where strong gravitational lensing masses favour MW-like IMFs. Like the high-σ sample on average, these galaxies have strong Na I 1.14 μm; taken in isolation their sodium indices imply bottom-heavy IMFs which are hard to reconcile with the lensing masses. An alternative full-spectrum-fitting approach, applied to the high-σsample, recovers an IMF less heavy than Salpeter, but under-predicts the Na I 1.14 μm line at the 5σ level. We conclude that current models struggle to reproduce this feature in the most massive galaxies without breaking other constraints, and caution against over-reliance on the sodium lines in spectroscopic IMF studies.

preprint2015arXiv

The SINFONI Nearby Elliptical Lens Locator Survey: Discovery of two new low-redshift strong lenses and implications for the initial mass function in giant early-type galaxies

We present results from a blind survey to identify strong gravitational lenses among the population of low-redshift early-type galaxies. The SINFONI Nearby Elliptical Lens Locator Survey (SNELLS) uses integral-field infrared spectroscopy to search for lensed emission line sources behind massive lens candidates at $z$<0.055. From 27 galaxies observed, we have recovered one previously-known lens (ESO325-G004) at $z$=0.034, and discovered two new systems, at $z$=0.031 and $z$=0.052. All three lens galaxies have high velocity dispersions (σ>300 km/s) and α-element abundances ([Mg/Fe]>0.3). From the lensing configurations we derive total J-band mass-to-light ratios of 1.8$\pm$0.1, 2.1$\pm$0.1 and 1.9$\pm$0.2 within the $\sim$2 kpc Einstein radius. Correcting for estimated dark-matter contributions, and comparing to stellar population models with a Milky Way (Kroupa) initial mass function (IMF), we determine the "mass excess factor", α. Assuming the lens galaxies have "old" stellar populations (10$\pm$1 Gyr), the average IMF mass factor is $\langleα\rangle$=1.10$\pm$0.08$\pm$0.10, where the first error is random and the second is systematic. If we instead fit the stellar populations from 6dF optical survey spectra, all three galaxies are consistent with being old, but the age errors are 3-4 Gyr, due to limited signal-to-noise. The IMF constraints are therefore looser in this case, with $\langleα\rangle$ = $1.23^{+0.16}_{-0.13}\pm{0.10}$. Our results are thus consistent with a Kroupa IMF (α=1.00) on average, and strongly reject very heavy IMFs with α>2. A Salpeter IMF (α=1.55) is inconsistent at the 3.5$σ$ level if the galaxies are old, but cannot be excluded using age constraints derived from the currently-available optical spectra.

preprint2014arXiv

Barred S0 Galaxies in the Coma Cluster

This study uses r-band images from the Eighth Data Release of the Sloan Digital Sky Survey (SDSS DR8) to study bars in lenticular (S0) galaxies in one of the nearest rich cluster environments, the Coma cluster. We develop techniques for bar detection, and assess their success when applied to SDSS image data. To detect and characterise bars we perform 2D bulge+disk+bar light decompositions of galaxy images with GALFIT. Using a sample of artificial galaxy images we determine the faintest magnitude at which bars can be successfully measured at the depth and resolution of SDSS. We perform detailed decompositions of 83 S0 galaxies in Coma, 64 from a central sample, and 19 from a cluster outskirts sample. For the central sample, the S0 bar fraction is 72^{+5}_{-6}%. This value is significantly higher than that obtained using an ellipse fitting method for bar detection, 48^{+6}_{-6}%. At a fixed luminosity, barred S0s are redder in (g-r) colour than unbarred S0s by 0.02 mag. The frequency and strength of bars increase towards fainter luminosities. Neither central metallicity nor stellar age distributions differ significantly between barred and unbarred S0s. There is an increase in the bar fraction towards the cluster core, but this is at a low significance level. Bars have at most a weak correlation with cluster-centric radius.

preprint2014arXiv

Dissecting the Red Sequence: The Bulge and Disc Colours of Early-Type Galaxies in the Coma Cluster

We explore the internal structure of red sequence galaxies in the Coma cluster across a wide range of luminosities ($-17>M_g>-22$) and cluster-centric radii ($0<r_{\rm{cluster}}<1.3 r_{200}$). We present the 2D bulge-disc decomposition of galaxies in deep Canada-France-Hawaii Telescope $u,g,i$ imaging using GALFIT. Rigorous filtering is applied to identify an analysis sample of 200 galaxies which are well described by an `archetypal' S0 structure (central bulge + outer disc). We consider internal bulge and/or disc colour gradients by allowing component sizes to vary between bands. Gradients are required for $30\%$ of analysis sample galaxies. Bulge half-light radii are found to be uncorrelated with galaxy luminosity ($R_e \sim 1$ kpc, $n\sim2$) for all but the brightest galaxies ($M_g<-20.5$). The S0 discs are brighter (at fixed size, or smaller at fixed luminosity) than those of star-forming spirals. A similar colour-magnitude relation is found for both bulges and discs. The global red sequence for S0s in Coma hence results from a combination of both component trends. We measure an average bulge $-$ disc colour difference of $0.09\pm0.01$ mag in $g-i$, and $0.16\pm0.01$ mag in $u-g$. Using simple stellar population models, bulges are either $\sim2$-$3\times$ older, or $\sim2\times$ more metal-rich than discs. The trend towards bluer global S0 colours observed further from Coma's core is driven by a significant correlation in disc colour with cluster-centric radius. An equivalent trend is detected in bulge colours at a marginal significance level. Our results therefore favour environment-mediated mechanisms of disc fading as the dominant factor in S0 formation.

preprint2014arXiv

Variations in the initial mass function in early-type galaxies: A critical comparison between dynamical and spectroscopic results

I present a comparison between published dynamical (ATLAS3D) and spectroscopic (Conroy & van Dokkum) constraints on the stellar initial mass function (IMF) in early-type galaxies, using the 34 galaxies in common between the two works. Both studies infer an average IMF mass factor $α$ (the stellar mass relative to a Kroupa-IMF population of similar age and metallicity) greater than unity, i.e. both methods favour an IMF which is heavier than that of the Milky Way, on average over the sample. However, on a galaxy-by-galaxy basis, there is no correlation between $α$ inferred from the two approaches. I investigate how the two estimates of $α$ are correlated systematically with the galaxy velocity dispersion, $σ$, and with the Mg/Fe abundance ratio. The spectroscopic method, based on the strengths of metal absorption lines, yields a correlation only with metal abundance ratios: at fixed Mg/Fe, there is no residual correlation with $σ$. The dynamical method, applied to exactly the same galaxy sample, yields the opposite result: the IMF variation correlates only with dynamics, with no residual correlation with Mg/Fe after controlling for $σ$. Hence although both methods indicate a heavy IMF on average in ellipticals, they lead to incompatible results for the systematic trends, when applied to the same set of galaxies. The sense of the disagreement could suggest that one (or both) of the methods has not accounted fully for the main confounding factors, i.e. element abundance ratios or dark matter contributions. Alternatively, the poor agreement might indicate additional variation in the detailed shape of the IMF, beyond what can currently be inferred from the spectroscopic features.

preprint2013arXiv

A giant elliptical galaxy with a lightweight initial mass function

We present new VLT observations of the closest-known strong-lensing galaxy, the sigma=330 km/s giant elliptical ESO325-G004. The low redshift of the lens (z_lens=0.035) results in arcs being formed at a small fraction of the effective radius, (R_Ein ~ R_Eff/4). At such small radii, stars dominate the lensing mass, so that lensing provides a direct probe of the stellar mass-to-light ratio, with only small corrections needed for dark matter. The redshift of the galaxy lensed by ESO325-G004 was unknown until now, so the lensing mass was not securely determined. Using X-SHOOTER, we measure a secure source redshift of z_src=2.141. Combined with the lensing configuration, this yields a total mass inside the Einstein radius of (1.50+/-0.06) 10^11 M_sun. We estimate the range of possible contribution of dark matter to the lensing mass, using statistics from cosmological N-body simulations. Subtracting this component yields a stellar M/L ratio for the lens of 3.14^+0.24_-0.42 (M/L)_sun (in F814W band). Using VIMOS, we have also obtained very high-S/N spectroscopy for the lens galaxy. The VIMOS data confirm that ESO325-G004 has a very old stellar population. For a Kroupa IMF, the stellar population fit yields a predicted stellar M/L ratio of 3.01+/-0.25 (M/L)_sun. Hence the mass attributable to stars with this IMF is consistent with the lensing estimate. By contrast, a Salpeter or heavier IMF is disfavoured at the 99.8% confidence level. A "heavyweight" IMF, with a mass twice as large as the Kroupa case, is firmly excluded. Such an IMF has been proposed for more distant elliptical lenses, and also to explain strong dwarf-star spectral features, in particular the NaI 8200-Ang doublet. A FORS2 far-red spectrum shows that this feature is as strong in ESO325-G004 as it is in other high-sigma ellipticals, suggesting tension between dwarf-star indicators and lensing-mass constraints for this galaxy.

preprint2013arXiv

Quenching star formation in cluster galaxies

In order to understand the processes that quench star formation within rich clusters, we construct a library of subhalo orbits drawn from $Λ$CDM cosmological N-body simulations of four rich clusters. The orbits are combined with models of star formation followed by quenching in the cluster environment. These are compared with observed bulge and disc colours and stellar absorption linestrength indices of satellite galaxies. Models in which the bulge stellar populations depend only on the galaxy subhalo mass while the disc quenching depends on the cluster environment are acceptable fits to the data. An exponential disc quenching timescale of 3 - 3.5 Gyr is preferred. Models with short ($\lesssim 1$ Gyr) quenching timescales yield cluster-centric gradients in disc colours and Balmer line indices that are too steep compared to observations. We also examine models in which there is quenching in lower mass groups prior to cluster infall ("pre-processing"), finding that such models are a better fit to the data than models without pre-processing and require similar quenching times. The data slightly prefer models where quenching occurs only for galaxies falling within about 0.5 $r_{200}$. Finally, we have examined models with short quenching timescales of 1 Gyr, but a long delay time of 3 Gyr prior to quenching. All models with short quenching timescales, even such "delayed-then-rapid" quenching models, produce excessively red galaxies near the cluster core and are strongly disfavoured by the data. These results imply that the environments of rich clusters must impact star formation rates of infalling galaxies on relatively long timescales -- several times longer than a typical halo spends within the virial radius of a cluster. This scenario favours gentler quenching mechanisms such as slow "strangulation" over more rapid ram-pressure stripping.

preprint2013arXiv

S0 galaxies in the Coma cluster: Environmental dependence of the S0 offset from the Tully-Fisher relation

We present deep GMOS long-slit spectroscopy of 15 Coma cluster S0 galaxies, and extract kinematic properties along the major axis to several times the disc scale-length. Supplementing our dataset with previously published data, we create a combined sample of 29 Coma S0s, as well as a comparison sample of 38 Coma spirals. Using photometry from SDSS and 2MASS, we construct the Tully-Fisher relation (TFR; luminosity versus maximum rotational velocity) for S0 galaxies. At fixed rotational velocity, the Coma S0 galaxies are on average fainter than Coma spirals by 1.10$\pm$0.18, 0.86$\pm$0.19 and 0.83$\pm$0.19 mag in the g, i and Ks bands respectively. The typical S0 offsets remain unchanged when calculated relative to large field-galaxy spiral samples. The observed offsets are consistent with a simple star formation model in which S0s are identical to spirals until abrupt quenching occurs at some intermediate redshift. The offsets form a continuous distribution tracing the time since the cessation of star formation, and exhibit a strong correlation (>6σ) with residuals from the optical colour-magnitude relation. Typically, S0s which are fainter than average for their rotational velocity are also redder than average for their luminosity. The S0 TFR offset is also correlated with both the projected cluster-centric radius and the Σ (projected) local density parameter. Since current local environment is correlated with time of accretion into the cluster, our results support a scenario in which transformation of spirals to S0s is triggered by cluster infall.

preprint2012arXiv

The stellar initial mass function in red-sequence galaxies: 1-micron spectroscopy of Coma Cluster galaxies with Subaru/FMOS

To investigate possible variations in the stellar initial mass function (IMF) in red-sequence galaxies, we have obtained infrared spectroscopy with Subaru/FMOS for a sample of 92 red-sequence galaxies in the Coma cluster. Velocity dispersions, ages and element abundances for these galaxies have been previously determined from optical data. By stacking the FMOS spectra in the rest frame, removing sky-subtraction residuals and other artefacts fixed in the observed frame, we derive composite spectra in the 9600-10500 Angstrom range for galaxies grouped according to their velocity dispersion or Mg/Fe ratio. We measure the Wing-Ford band of FeH and a new index centred on a CaI line at 10345 Angstrom; these features are strong in cool dwarf stars, and hence reflect the form of the IMF at low mass (<0.5M_sun) The CaI line, unlike the Wing-Ford band and other `classical' IMF indicators (NaI doublet, CaII triplet), is unaffected by the abundance of sodium. We compare the measured indices against predictions from spectral synthesis models matched to the element abundances estimated from the optical data. Binning galaxies by velocity dispersion, we find that both the Wing-Ford and CaI index measurements are best reproduced by models with the Salpeter IMF. There is no clear evidence for an increase in dwarf-star content with velocity dispersion. Binning the observed galaxies instead by Mg/Fe ratio, the behaviour of both indices implies a trend of IMF from Chabrier-like, at abundance ratios close to solar, to Salpeter or heavier for highly alpha-enhanced populations. At face value, this suggests that the IMF depends on the mode of star formation, with intense rapid star-bursts generating a larger population of low-mass stars.

preprint2012arXiv

What drives the ultra-violet colours of passive galaxies?

We present and analyse optical and ultra-violet colours for passive and optically-red Coma cluster galaxies for which we have spectroscopic age and element abundance estimates. Our sample of 150 objects covers a wide range in mass, from giant ellipticals to the bright end of the dwarf-galaxy regime. We focus on the colours FUV-i, NUV-i, FUV-NUV, u*-g and g-i. We find that all of these colours are correlated with both luminosity and velocity dispersion at the >5 sigma level, with FUV-i and FUV-NUV becoming bluer with increasing `mass' while the other colours become redder. We perform an empirical analysis to assess what fraction of the variation in each colour can be accounted for by variations in the average stellar populations, as traced by the optical spectra. For u*-g and g-i, most of the observed scatter (~80% after allowing for measurement errors and for systematic errors in u*-g) is attributable to stellar population variations, with colours becoming redder with increasing age and metallicity (Mg/H). The FUV-i colour becomes bluer with increasing age and with increasing Mg/H, favouring a `metal-rich single-star' origin for the UV upturn. However, correlations with the optically-dominant stellar populations account for only about half of the large observed scatter. We propose that the excess scatter in FUV-i may be due to a varying proportion of ancient stars in galaxies with younger average ages. The NUV-i colour is sensitive to age and Mg/H, but exhibits excess scatter that can be attributed to `leakage' of the FUV upturn. Applying a correction based on the FUV-i colour we account for ~80% of the variance in NUV-i, as for the optical colours. The FUV-NUV colour shows strong correlations with age and Mg/H, and little residual scatter. Interpreting this colour is complicated however, since it mixes the effects from the main-sequence turn-off with those from the hot evolved stars.

preprint2011arXiv

Environmental quenching and hierarchical cluster assembly: Evidence from spectroscopic ages of red-sequence galaxies in Coma

We explore the variation in stellar population ages for Coma cluster galaxies as a function of projected cluster-centric distance, using a sample of 362 red-sequence galaxies with high signal-to-noise spectroscopy. The sample spans a wide range in luminosity (0.02-4 L*) and extends from the cluster core to near the virial radius. We find a clear distinction in the observed trends of the giant and dwarf galaxies. The ages of red-sequence giants are primarily determined by galaxy mass, with only weak modulation by environment, in the sense that galaxies at larger cluster-centric distance are slightly younger. For red-sequence dwarfs (with mass <10^10 Msun), the roles of mass and environment as predictors of age are reversed: there is little dependence on mass, but strong trends with projected cluster-centric radius are observed. The average age of dwarfs at the 2.5 Mpc limit of our sample is approximately half that of dwarfs near the cluster centre. The gradient in dwarf galaxy ages is a global cluster-centric trend, and is not driven by the ongoing merger of the NGC 4839 group to the south west of Coma. We interpret these results using environmental histories extracted from the Millennium Simulation for members of massive clusters. Hierarchical cluster assembly naturally leads to trends in the accretion times of galaxies as a function of projected cluster-centric radius. On average, simulated galaxies now located in cluster cores joined halos above any given mass threshold earlier than those now located in the outskirts of clusters. We test environmental quenching models, in which star formation is halted in galaxies when they enter halos of a given mass, or become satellite galaxies. The models broadly reproduce the gradients observed in Coma, but for dwarf galaxies the efficiency of environmental quenching must be very high to match the strong trends observed.

preprint2011arXiv

The HST/ACS Coma Cluster Survey IV. Intergalactic Globular Clusters and the Massive Globular Cluster System at the Core of the Coma Galaxy Cluster

Intracluster stellar populations are a natural result of tidal interactions in galaxy clusters. Measuring these populations is difficult, but important for understanding the assembly of the most massive galaxies. The Coma cluster is one of the nearest truly massive galaxy clusters, and is host to a correspondingly large system of globular clusters (GCs). We use imaging from the HST/ACS Coma Cluster Survey to present the first definitive detection of a large population of intracluster GCs (IGCs) that fills the Coma cluster core and is not associated with individual galaxies. The GC surface density profile around the central massive elliptical galaxy, NGC 4874, is dominated at large radii by a population of IGCs that extend to the limit of our data (R<520 kpc). We estimate that there are 47000+/-1600 (random) +4000/-5000 (systematic) IGCs out to this radius, and that they make up ~70% of the central GC system, making this the largest GC system in the nearby Universe. Even including the GC systems of other cluster galaxies, IGCs still make up ~30-45% of the GCs in the cluster core. Observational limits from previous studies of the intracluster light (ICL) suggest that the IGC population has a high specific frequency. If the IGC population has a specific frequency similar to high-S_N dwarf galaxies, then the ICL has a total stellar mass of ~10^12 M_sun within the cluster core. The ICL makes up approximately half of the stellar luminosity and one-third of the stellar mass of the central (NGC4874+ICL) system. The color distribution of the IGC population is bimodal, with blue, metal-poor GCs outnumbering red, metal-rich GCs by a ratio of 4:1. The fraction of red IGCs (20%), and the red color of those GCs, implies that IGCs can originate from the halos of relatively massive, L* galaxies, and not solely from the disruption of dwarf galaxies. (Abridged)

preprint2011arXiv

The spatial distribution and origin of the FUV excess in early-type galaxies

We present surface photometry of a sample of 52 galaxies from the GALEX and 2MASS data archives, these include 32 normal elliptical galaxies, 10 ellipticals with weak Liner or other nuclear activity, and 10 star forming ellipticals or early-type spirals. We examine the spatial distribution of the Far Ultra-Violet excess in these galaxies, and its correlation with dynamical and stellar population properties of the galaxies. From aperture photometry we find that all galaxies except for recent major remnants and galaxies with ongoing star formation show a positive gradient in the (FUV-NUV) colour determined from the GALEX images. The logarithmic gradient does not correlate with any stellar population parameter, but it does correlate with the central velocity dispersion. The strength of the excess on the other hand, correlates with both [alpha/Fe] and [Z/H], but more strongly with the former. We derive models of the underlying stellar population from the 2MASS H-band images, and the residual of the image from this model reveals a map of the centrally concentrated FUV excess. We examine a possible hypothesis for generating the FUV excess and the radial gradient in its strength, involving a helium abundance gradient set up early in the formation process of the galaxies. If this hypothesis is correct, the persistence of the gradients to the present day places a strong limit on the importance of dry mergers in the formation of ellipticals.

preprint2010arXiv

ACCESS II: A Complete Census of Star Formation in the Shapley Supercluster - UV and IR Luminosity Functions

We present panoramic Spitzer/MIPS mid- and far-infrared and GALEX ultraviolet imaging of the the most massive and dynamically active system in the local Universe, the Shapley supercluster at z=0.048, covering the 5 clusters which make up the supercluster core. We combine these data with existing spectroscopic data from 814 confirmed supercluster members to produce the first study of a local rich cluster including both ultraviolet and infrared luminosity functions (LFs). This joint analysis allows us to produce a complete census of star-formation (both obscured and unobscured), extending down to SFRs~0.02-0.05Msun/yr, and quantify the level of obscuration of star formation among cluster galaxies, providing a local benchmark for comparison to ongoing and future studies of cluster galaxies at higher redshifts with Spitzer and Herschel. The GALEX NUV and FUV LFs obtained have steeper faint-end slopes than the local field population, due largely to the contribution of massive, quiescent galaxies at M_FUV>-16. The 24um and 70um galaxy LFs for the Shapley supercluster instead have shapes fully consistent with those obtained for the Coma cluster and for the local field galaxy population. This apparent lack of environmental dependence for the shape of the FIR luminosity function suggests that the bulk of the star-forming galaxies that make up the observed cluster infrared LF have been recently accreted from the field and have yet to have their star formation activity significantly affected by the cluster environment. We estimate a global SFR of 327 Msun/yr over the whole supercluster core, of which just ~20% is visible directly in the UV continuum and ~80% is reprocessed by dust and emitted in the infrared. The level of obscuration (L_IR/L_FUV) in star-forming galaxies is seen to increase linearly with L_K over two orders of magnitude in stellar mass.

preprint2010arXiv

ACCESS III: The Nature of Star Formation in the Shapley Supercluster

We present a joint analysis of panoramic Spitzer/MIPS mid-infrared and GALEX ultraviolet imaging of the Shapley supercluster at z=0.048. Combining this with spectra of 814 supercluster members and 1.4GHz radio continuum maps, this represents the largest complete census of star-formation (both obscured and unobscured) in local cluster galaxies to date, reaching SFRs~0.02Msun/yr. We take advantage of this comprehensive panchromatic dataset to perform a detailed analysis of the nature of star formation in cluster galaxies, using several quite independent diagnostics of the quantity and intensity of star formation to develop a coherent view of the types of star formation within cluster galaxies. We observe a robust bimodality in the infrared (f_24/f_K) galaxy colours, which we are able to identify as another manifestation of the broad split into star-forming spiral and passive elliptical galaxy populations seen in UV-optical surveys. This diagnostic also allows the identification of galaxies in the process of having their star formation quenched as the infrared analogue to the UV "green valley" population. The bulk of supercluster galaxies on the star-forming sequence have specific-SFRs consistent with local field specific-SFR-M* relations, and form a tight FIR-radio correlation confirming that their FIR emission is due to star formation. We show that 85% of the global SFR is quiescent star formation within spiral disks, as manifest by the observed sequence in the IRX-beta relation being significantly offset from the starburst relation of Kong et al. (2004), while their FIR-radio colours indicate dust heated by low-intensity star formation. Just 15% of the global SFR is due to nuclear starbursts. The vast majority of star formation seen in cluster galaxies comes from normal infalling spirals who have yet to be affected by the cluster environment.

preprint2010arXiv

The HST/ACS Coma Cluster Survey. II. Data Description and Source Catalogs

The Coma cluster was the target of a HST-ACS Treasury program designed for deep imaging in the F475W and F814W passbands. Although our survey was interrupted by the ACS instrument failure in 2007, the partially completed survey still covers ~50% of the core high-density region in Coma. Observations were performed for 25 fields that extend over a wide range of cluster-centric radii (~1.75 Mpc) with a total coverage area of 274 arcmin^2. The majority of the fields are located near the core region of Coma (19/25 pointings) with six additional fields in the south-west region of the cluster. In this paper we present reprocessed images and SExtractor source catalogs for our survey fields, including a detailed description of the methodology used for object detection and photometry, the subtraction of bright galaxies to measure faint underlying objects, and the use of simulations to assess the photometric accuracy and completeness of our catalogs. We also use simulations to perform aperture corrections for the SExtractor Kron magnitudes based only on the measured source flux and half-light radius. We have performed photometry for ~73,000 unique objects; one-half of our detections are brighter than the 10-sigma point-source detection limit at F814W=25.8 mag (AB). The slight majority of objects (60%) are unresolved or only marginally resolved by ACS. We estimate that Coma members are 5-10% of all source detections, which consist of a large population of unresolved objects (primarily GCs but also UCDs) and a wide variety of extended galaxies from a cD galaxy to dwarf LSB galaxies. The red sequence of Coma member galaxies has a constant slope and dispersion across 9 magnitudes (-21<M_F814W<-13). The initial data release for the HST-ACS Coma Treasury program was made available to the public in 2008 August. The images and catalogs described in this study relate to our second data release.

preprint2010arXiv

Ultraviolet tails and trails in cluster galaxies: A sample of candidate gaseous stripping events in Coma

We have used new deep observations of the Coma cluster from GALEX to identify 13 star-forming galaxies with asymmetric morphologies in the ultraviolet. Aided by optical broad-band and H-alpha imaging, we interpret the asymmetric features as being due to star formation within gas stripped from the galaxies by interaction with the cluster environment. The selected objects display a range of structures from broad fan-shaped systems of filaments and knots (`jellyfish') to narrower and smoother tails extending up to 100 kpc in length. Some of the features have been discussed previously in the literature, while others are newly identified here. As an ensemble, the candidate stripping events are located closer to the cluster centre than other star-forming galaxies; their radial distribution is similar to that of all cluster members, dominated by passive galaxies. The fraction of blue galaxies which are undergoing stripping falls from 40% in the central 500 kpc, to less than 5% beyond 1 Mpc. We find that tails pointing away from (i.e. galaxies moving towards) the cluster centre are strongly favoured (11/13 cases). From the small number of `outgoing' galaxies with stripping signatures we conclude that the stripping events occur primarily on first passage towards the cluster centre, and are short-lived compared to the cluster crossing time. Using infall trajectories from simulations, the observed fraction of blue galaxies undergoing stripping can be reproduced if the events are triggered at a threshold radius of ~1 Mpc and detectable for ~500 Myr. HST images are available for two galaxies from our sample and reveal compact blue knots coincident with UV and H-alpha emission, apparently forming stars within the stripped material. Our results confirm that stripping of gas from infalling galaxies, and associated star formation in the stripped material, is a widespread phenomenon in rich clusters.

preprint2009arXiv

Abundance ratios in red-sequence galaxies over a wide mass range: The "X-planes" for magnesium, calcium, carbon and nitrogen

We analyse the abundance ratios of the light elements Mg, Ca, C and N, relative to Fe, for 147 red-sequence galaxies in the Coma cluster and the Shapley Supercluster. The sample covers a six-magnitude range in luminosity, from giant ellipticals to dwarfs at M^*+4. We exploit the wide mass range to investigate systematic trends in the abundance ratios Mg/Fe, Ca/Fe, C/Fe and N/Fe. We find that each of these ratios can be well modelled using two-parameter relations of the form [X/Fe] = a0 + a1 log sigma + a2 [Fe/H], where sigma is the velocity dispersion. Analysing these "X-planes" reveals new structure in the abundance patterns, beyond the traditional one-parameter (e.g. Mg/Fe-sigma) correlations. The X-planes for the alpha elements, Mg and Ca, indicate a positive correlation with velocity dispersion, and simultaneously an anti-correlation with Fe/H (i.e. a1>0 and a2<0). Taking both effects into account dramatically reduces the scatter, compared to the traditional X/Fe-sigma relations. For C and N, a similar correlation with velocity dispersion is recovered, but there is no additional dependence on Fe/H (i.e. a1>0 and a2~0). The explicit dependence of X/Fe on two parameters is evidence that at least two physical processes are at work in setting the abundance patterns. The Fe/H dependence of Mg/Fe and Ca/Fe, at fixed sigma, may result from different durations of star formation, from galaxy to galaxy. The absence of corresponding Fe/H dependence for C and N is consistent with these elements being generated in lower-mass stars. The increase with sigma, at fixed Fe/H, is similar for elements Mg, C and N, and slightly shallower for Ca. This pattern of trends cannot be explained solely by a systematic variation of star-formation time-scale with sigma.

preprint2009arXiv

Ages and metallicities for quiescent galaxies in the Shapley Supercluster: Driving parameters of the stellar populations

We use high signal-to-noise spectroscopy for a sample of 232 quiescent galaxies in the Shapley Supercluster, to investigate how their stellar populations depend on velocity dispersion, luminosity and stellar mass. The sample spans a large range in velocity dispersion (sigma from 30-300 km/s) and in luminosity (M_R from -18.7 to -23.2). Estimates of age, total metallicity (Z/H) and alpha-element abundance ratio (a/Fe) were derived from absorption-line analysis, using single-burst models. Age, Z/H and a/Fe are all correlated positively with velocity dispersion, but we also find significant residual trends with luminosity: at given sigma, the brighter galaxies are younger, less alpha-enriched, and have higher Z/H. At face value, these results might suggest that the stellar populations depend on stellar mass as well as on velocity dispersion. However, we show that the observed trends can be reproduced by models in which the stellar populations depend systematically only on sigma, and are independent of stellar mass Mstel. For age, the observed luminosity correlation arises because young galaxies are brighter, at fixed Mstel. For metallicity, the observed luminosity dependence arises because metal-rich galaxies, at fixed mass, tend also to be younger, and hence brighter. We find a good match to the observed luminosity correlations with Age ~ sigma^0.40, Z/H ~ sigma^0.35, a/Fe ~ sigma^0.20, where the slopes are close to those found when fitting traditional scaling relations. We conclude that the star formation and enrichment histories of galaxies are determined primarily by the depth of their gravitational potential wells. The observed residual correlations with luminosity do not imply a corresponding dependence on stellar mass.

preprint2009arXiv

Stellar population gradients in early-type cluster galaxies

We present a study of internal stellar population gradients in early-type cluster galaxies. Using the VLT VIMOS integral field unit, we observed 19 galaxies in the core of the Shapley Supercluster (z = 0.048). The radial trends in nine absorption lines (HdF to Fe5406) were measured to the effective radius for 14 galaxies, from which we derived the gradients in age, total metallicity and alpha-element over-abundance. We combine these with results from 11 galaxies studied in our previous VIMOS work (Rawle et al 2008). We observe a mean metallicity gradient of -0.13 +/- 0.04 per dex and, in common with the findings of previous studies, galaxies with log(sigma) > 2.1 have a sizeable intrinsic scatter in metallicity gradient. The mean log(age/Gyr) gradient is -0.02 +/- 0.06 per dex, although several galaxies have significant positive or negative age gradients. The mean gradient in alpha-element enhancement is -0.10 +/- 0.04 per dex. We find that stellar population gradients are primarily related to the central metallicity: early-type galaxies with super-solar centres have steep negative metallicity gradients and positive age gradients; those with solar metallicity centres have negligible [Z/H] gradients and negative age gradients. There is a strong observed anti-correlation between the gradients in age and metallicity. While a part of this trend can be attributed to the correlation of measurement errors, we demonstrate that there is an underlying intrinsic relation. For the Shapley galaxies, B-R colour gradients predicted from spectroscopic age and metallicity generally agree well with those measured directly from photometry.

preprint2005arXiv

NOAO Fundamental Plane Survey -- II. Age and Metallicity along the Red Sequence

We present spectroscopic linestrength data for 4097 red-sequence galaxies in 93 low-redshift galaxy clusters, and use these to investigate variations in average stellar populations as a function of galaxy mass. Our analysis includes an improved treatment of nebular emission contamination, which affects ~10% of the sample galaxies. Using the stellar population models of D. Thomas and collaborators, we simultaneously fit twelve observed linestrength-sigma relations in terms of common underlying trends of age, [Z/H] (total metallicity) and a/Fe (alpha-element enhancement). We find that the observed linestrength-sigma relations can be explained only if higher-mass red-sequence galaxies are, on average, older, more metal rich, and more alpha-enhanced than lower-mass galaxies. Quantitatively, the scaling relations are age=sigma^(0.59+/-0.13), Z/H=sigma^(0.53+/-0.08) and a/Fe=sigma^(0.31+/-0.06), where the errors reflect the range obtained using different subsets of indices. We conclude that although the stars in giant red galaxies in clusters formed early, most of the galaxies at the faint end joined the red sequence only at recent epochs. This "down-sizing" trend is in good qualitative agreement with observations of the red sequence at higher redshifts, but is not predicted by semi-analytic models of galaxy formation.

preprint2005arXiv

Synthesis of H-alpha absorption in old stellar systems: Formation of the cluster red sequence by "Downsizing"

We compute population synthesis models for the variation of H-alpha absorption indices (HaA and HaF), as a function of age and metallicity in old stellar systems. The models are based on the STELIB spectral library, and defined at resolution of 3 Angstrom FWHM. The indices are found to be highly sensitive to age variation, with only moderate response to metallicity. For galaxies uncontaminated by nebular emission, our HaA index is more powerful in breaking the age/metallicity degeneracy than H-beta or H-gamma_F. Using a sample of red cluster galaxies carefully selected to exclude objects with emission, we find a steep decline of HaA with velocity dispersion (slope -0.75 +/- 0.07 Angstrom/dex). The slope can be translated to constraints on age and metallicity scaling relations, incorporating measurement and model errors. If the HaA-sigma slope is due only to age, we obtain Age proportional to sigma^{0.95 +/- 0.12}. Because HaA depends quite weakly on [Fe/H], a metallicity interpretation would require Fe/H proportional to sigma^1.2 or steeper. The HaA-sigma slope is consistent with the combined age and metallicity scaling relations reported by Nelan et al. from classical Lick indices. The relations obtained by Thomas et al. significantly under-predict the observed slope. The discrepancy could arise from differences in the sample compositions.We discuss in detail the effects of emission contamination on the results, and conclude that these are unlikely to yield the observed behaviour in the Ha-sigma relations. The steep age--mass relation supports a "downsizing" formation scenario: fainter red-sequence galaxies became quiescent at lower redshifts, z < 0.5. This picture accords with recent observations of truncated red sequences in clusters at z ~ 0.7.

preprint2000arXiv

On the dependence of spectroscopic indices of early-type galaxies on age, metallicity and velocity dispersion

We investigate the Mg-sigma and <Fe>-sigma relations in a sample of 72 early-type galaxies drawn mostly from cluster and group environments using a homogeneous data-set which is well-calibrated onto the Lick/IDS system. The small intrinsic scatter in Mg at a given sigma gives upper limits on the spread in age and metallicity of 49% and 32% respectively, if the spread is attributed to one quantity only and if the variations in age and metallicity are uncorrelated. The age/metallicity distribution as inferred from the Hbeta vs <Fe> diagnostic diagram reinforces this conclusion, as we find mostly galaxies with large luminosity weighted ages spanning a range in metallicity. In our sample we do not find significant evidence for an anti-correlation of ages and metallicities which would keep the index-sigma relations tight while hiding a large spread in age and metallicity. As a result of correlated errors in the age-metallicity plane, a mild age-metallicity anti-correlation cannot be completely ruled out given the current data. Correcting the line-strengths indices for non-solar abundance ratios following the recent paper by Trager et al., leads to higher mean metallicity and slightly younger age estimates while preserving the metallicity sequence. The [Mg/Fe] ratio is mildly correlated with the central velocity dispersion and ranges from [Mg/Fe]=0.05 to 0.3 for galaxies with sigma > 100 km/s. Under the assumption that there is no age gradient along the index-sigma relations, the abundance-ratio corrected Mg-sigma, <Fe>-sigma and Hbeta-sigma relations give consistent estimates of Delta [M/H]/ Delta log sigma = 0.9 (+- 0.1). The slope of the Hbeta-sigma relation limits a potential age trend as a function of sigma to 2-3 Gyrs along the sequence.(abriged)