Search

Environment has long been known to impact the evolution of galaxies, but disentangling its effects from mass evolution requires careful analysis of statistically significant samples. By implementing advanced visualisation methods to test group-finding algorithms, we utilise a mass-complete sample of galaxies to z < 0.1, comprising spectroscopic redshifts from prominent surveys such as the 2dFGRS and GAMA. Our group-finding methods identify 1,413 galaxy groups made up of 8,990 galaxies, corresponding to 36% of galaxies associated with group environments. We also search for close pairs, with separations of $r_{sep}$ < 50 h$^{-1}$ kpc and $v_{sep}$ < 500 km/s, and classify them into major ($M_{sec}/M_{prim} \leq$ 0.25) and minor ($M_{sec}/M_{prim}$ > 0.25) pairs. To examine the impact of environmental factors, we employ bespoke WISE photometry to derive a star-forming main sequence relation that shows star-formation (SF) within galaxies is pre-processed as a function of group membership. Our analysis reveals that SF in galaxies is pre-processed as a function of group membership. We observe an increase in the fraction of quiescent galaxies relative to the field as group membership rises, quantified using the environmental quenching efficiency metric ($\epsilon_{env}$). Within the star-forming population, we detect pre-processing with the relative difference in specific SF rates (${\Delta} sSFR$), showing a net decrease in SF as group membership increases, particularly at larger stellar masses. Our sample of close pairs at low stellar masses shows enhanced SF efficiencies compared to the field, while at larger masses, deficiencies are evident. Our results indicate that the small-scale environments of galaxies influence SF properties, demonstrating that galaxies do not evolve in isolation over cosmic time but are shaped by complex interactions between internal dynamics and external influences., Comment: Accepted by PASA 17/10/2024. 24 pages, 22 figures, 5 tables. Abridged abstract

Nebular emission lines are powerful diagnostics for the physical processes at play in galaxy formation and evolution. Moreover, emission-line galaxies (ELGs) are one of the main targets of current and forthcoming spectroscopic cosmological surveys. We investigate the contributions to the line luminosity functions (LFs) of different galaxy populations in the local Universe, providing a benchmark for future surveys of earlier cosmic epochs. The large statistics of the observations from the SDSS DR7 Main galaxy sample and the MPA-JHU spectral catalogue enabled us to precisely measure the H$\alpha$, H$\beta$, [OII], [OIII], and, for the first time, the [NII], and [SII] emission-line LFs over ~2.4 Gyrs in the low-z Universe, 0.0210$^{-11}$/yr. When classifying ELGs using the BPT+WHAN diagrams, we find that 63.3 per cent are star forming, only 0.03 are passively evolving, and 1.3 have nuclear activity (Seyfert). The rest are low-ionization narrow emission-line regions (LINERs) and composite ELGs. We found that a Saunders function is the most appropriate to describe all of the emission-line LFs, both observed and dust-extinction-corrected. They are dominated by star-forming regions, except for the bright end of the [OIII] and [NII] LFs (i.e. L[NII]>10$^{42}$ erg/s, L[OIII]>10$^{43}$ erg/s), where the contribution of Seyfert galaxies is not negligible.

In the submm regime, spectral line scans and line intensity mapping (LIM) are new promising probes for the cold gas content and star formation rate of galaxies across cosmic time. However, both of these two measurements suffer from field-to-field variance. We study the effect of field-to-field variance on the predicted CO and [CII] power spectra from future LIM experiments such as CONCERTO, as well as on the line luminosity functions (LFs) and the cosmic molecular gas mass density that are currently derived from spectral line scans. We combined a 117 $\rm deg^2$ dark matter lightcone from the Uchuu cosmological simulation with the simulated infrared dusty extragalactic sky (SIDES) approach. We find that in order to constrain the CO LF with an uncertainty below 20%, we need survey sizes of at least 0.1 $\rm deg^2$. Furthermore, accounting for the field-to-field variance using only the Poisson variance can underestimate the total variance by up to 80%. The lower the luminosity is and the larger the survey size is, the higher the level of underestimate. At $z$<3, the impact of field-to-field variance on the cosmic molecular gas density can be as high as 40% for the 4.6 arcmin$^2$ field, but drops below 10% for areas larger than 0.2 deg$^2$. However, at $z>3$ the variance decreases more slowly with survey size and for example drops below 10% for 1 deg$^2$ fields. Finally, we find that the CO and [CII] LIM power spectra can vary by up to 50% in $\rm 1 deg^2$ fields. This limits the accuracy of the constraints provided by the first 1 deg$^2$ surveys. The level of the shot noise power is always dominated by the sources that are just below the detection thresholds. We provide an analytical formula to estimate the field-to-field variance of current or future LIM experiments. The code and the full SIDES-Uchuu products (catalogs, cubes, and maps) are publicly available.

We use three semi-analytic models (SAMs) of galaxy formation and evolution, run on the same 1$h^{-1}$Gpc MultiDark Planck2 cosmological simulation, to investigate the properties of [OII] emission line galaxies in the redshift range $0.6

We present a simple semi-numerical model designed to explore black hole growth and galaxy evolution. This method builds on a previous model for black hole accretion that uses a semi-numerical galaxy formation model and universal Eddington ratio distribution to describe the full AGN population by independently connecting galaxy and AGN growth to the evolution of the host dark matter halos. We fit observed X-ray luminosity functions up to a redshift of z ~ 4, as well as investigate the evolution of the Eddington ratio distributions. We find that the Eddington ratio distribution evolves with redshift such that the slope of the low-Eddington accretion rate distribution increases with cosmic time, consistent with the behavior predicted in hydrodynamical simulations for galaxies with different gas fractions. We also find that the evolution of our average Eddington ratio is correlated with observed star formation histories, supporting a picture in which black holes and galaxies evolve together in a global sense. We further confirm the impact of luminosity limits on observed galaxy and halo properties by applying selection criteria to our fiducial model and comparing to surveys across a wide range of redshifts., Comment: 9 pages, 8 figures, Accepted to ApJ

We introduce the 4-metre Multi-Object Spectroscopic Telescope (4MOST), a new high-multiplex, wide-field spectroscopic survey facility under development for the four-metre-class Visible and Infrared Survey Telescope for Astronomy (VISTA) at Paranal. Its key specifications are: a large field of view (FoV) of 4.2 square degrees and a high multiplex capability, with 1624 fibres feeding two low-resolution spectrographs ($R = \lambda/\Delta\lambda \sim 6500$), and 812 fibres transferring light to the high-resolution spectrograph ($R \sim 20\,000$). After a description of the instrument and its expected performance, a short overview is given of its operational scheme and planned 4MOST Consortium science; these aspects are covered in more detail in other articles in this edition of The Messenger. Finally, the processes, schedules, and policies concerning the selection of ESO Community Surveys are presented, commencing with a singular opportunity to submit Letters of Intent for Public Surveys during the first five years of 4MOST operations., Comment: Part of the 4MOST issue of The Messenger, published in preparation of 4MOST Community Workshop, see http://www.eso.org/sci/meetings/2019/4MOST.html

The primary difficulty in measuring dynamical masses of galaxy clusters from galaxy data lies in the separation between true cluster members from interloping galaxies along the line of sight. We study the impact of membership contamination and incompleteness on cluster mass estimates obtained with 25 commonly used techniques applied to nearly 1000 mock clusters. We show that all methods overestimate or underestimate cluster masses when applied to contaminated or incomplete galaxy samples respectively. This appears to be the main source of the intrinsic scatter in the mass scaling relation. Applying corrections based on a prior knowledge of contamination and incompleteness can reduce the scatter to the level of shot noise expected for poorly sampled clusters. We establish an empirical model quantifying the effect of imperfect membership on cluster mass estimation and discuss its universal and method-dependent features. We find that both imperfect membership and the response of the mass estimators depend on cluster mass, effectively causing a flattening of the estimated - true mass relation. Imperfect membership thus alters cluster counts determined from spectroscopic surveys, hence the cosmological parameters that depend on such counts., Comment: 18 pages, 11 figures, 3tables; accepted for publication in MNRAS

While major mergers have long been proposed as a driver of both AGN activity and the M-sigma relation, studies of moderate to high redshift Seyfert-luminosity AGN hosts have found little evidence for enhanced rates of interactions. However, both theory and observation suggest that while these AGN may be fueled by stochastic accretion and secular processes, high-luminosity, high-redshift, and heavily obscured AGN are the AGN most likely to be merger-driven. To better sample this population of AGN, we turn to infrared selection in the CANDELS/COSMOS field. Compared to their lower-luminosity and less obscured X-ray-only counterparts, IR-only AGN (luminous, heavily obscured AGN) are more likely to be classified as either irregular (50$^{+12}_{-12}$% vs. 9$^{+5}_{-2}$%) or asymmetric (69$^{+9}_{-13}$% vs. 17$^{+6}_{-4}$%) and are less likely to have a spheroidal component (31$^{+13}_{-9}$% vs. 77$^{+4}_{-6}$%). Furthermore, IR-only AGN are also significantly more likely than X-ray-only AGN (75$^{+8}_{-13}$% vs. 31$^{+6}_{-6}$%) to be classified either as interacting or merging in a way that significantly disturbs the host galaxy or disturbed though not clearly interacting or merging, which potentially represents the late stages of a major merger. This suggests that while major mergers may not contribute significantly to the fueling of Seyfert luminosity AGN, interactions appear to play a more dominant role in the triggering and fueling of high-luminosity heavily obscured AGN., Comment: 13 pages, 7 figures, 1 table, Accepted for publication in the Astrophysical Journal

JWST was conceived and built to answer one of the most fundamental questions that humans can address empirically: "How did the Universe make its first stars?". Our First Lights At REionization (FLARE) project transforms the quest for the epoch of reionization from the static to the time domain. It targets the complementary question: "What happened to those first stars?". It will be answered by observations of the most luminous events: supernovae and accretion on to black holes formed by direct collapse from the primordial gas clouds. These transients provide direct constraints on star-formation rates (SFRs) and the truly initial Initial Mass Function (IMF), and they may identify possible stellar seeds of supermassive black holes (SMBHs). Furthermore, our knowledge of the physics of these events at ultra-low metallicity will be much expanded. JWST's unique capabilities will detect these most luminous and earliest cosmic messengers easily in fairly shallow observations. However, these events are very rare at the dawn of cosmic structure formation and so require large area coverage. Time domain astronomy can be advanced to an unprecedented depth by means of a shallow field of JWST reaching 27 mag AB in 2 and 4.4 microns over a field as large as 0.1 square degree visited multiple times each year. Such a survey may set strong constraints or detect massive Pop III SNe at redshifts beyond 10, pinpointing the redshift of the first stars, or at least their death. Based on our current knowledge of superluminous supernovae (SLSNe), such a survey will find one or more SLSNe at redshifts above 6 in five years and possibly several direct collapse black holes. Although JWST is not designed as a wide field survey telescope, we show that such a wide field survey is possible with JWST and is critical in addressing several of its key scientific goals.

With the advent of wide-field cosmological surveys, we are approaching samples of hundreds of thousands of galaxy clusters. While such large numbers will help reduce statistical uncertainties, the control of systematics in cluster masses becomes ever more crucial. Here we examine the effects of an important source of systematic uncertainty in galaxy-based cluster mass estimation techniques: the presence of significant dynamical substructure. Dynamical substructure manifests as dynamically distinct subgroups in phase-space, indicating an 'unrelaxed' state. This issue affects around a quarter of clusters in a generally selected sample. We employ a set of mock clusters whose masses have been measured homogeneously with commonly-used galaxy-based mass estimation techniques (kinematic, richness, caustic, radial methods). We use these to study how the relation between observationally estimated and true cluster mass depends on the presence of substructure, as identified by various popular diagnostics. We find that the scatter for an ensemble of clusters does not increase dramatically for clusters with dynamical substructure. However, we find a systematic bias for all methods, such that clusters with significant substructure have higher measured masses than their relaxed counterparts. This bias depends on cluster mass: the most massive clusters are largely unaffected by the presence of significant substructure, but masses are significantly overestimated for lower mass clusters, by $\sim10\%$ at $10^{14}$ and $\geq20\%$ for $\leq10^{13.5}$. The use of cluster samples with different levels of substructure can, therefore, bias certain cosmological parameters up to a level comparable to the typical uncertainties in current cosmological studies., Comment: 14 pages, 8 figures, 3 tables. Submitted to MNRAS

The velocity function is a fundamental observable statistic of the galaxy population, similarly impor- tant as the luminosity function, but much more difficult to measure. In this work we present the first directly measured circular velocity function that is representative between 60 < v_circ < 320 km/s for galaxies of all morphological types at a given rotation velocity. For the low mass galaxy population (60 < v_circ < 170 km/s), we use the HIPASS velocity function. For the massive galaxy population (170 < v_circ < 320 km/s), we use stellar circular velocities from the Calar Alto Legacy Integral Field Area Survey (CALIFA). In earlier work we obtained the measurements of circular velocity at the 80% light radius for 226 galaxies and demonstrated that the CALIFA sample can produce volume- corrected galaxy distribution functions. The CALIFA velocity function includes homogeneous velocity measurements of both late and early-type rotation-supported galaxies and has the crucial advantage of not missing gas-poor massive ellipticals that HI surveys are blind to. We show that both velocity functions can be combined in a seamless manner, as their ranges of validity overlap. The resulting observed velocity function is compared to velocity functions derived from cosmological simulations of the z = 0 galaxy population. We find that dark matter-only simulations show a strong mismatch with the observed VF. Hydrodynamic simulations fare better, but still do not fully reproduce observations., Comment: Accepted to ApJL. Keywords: Galaxies: kinematics and dynamics - galaxies: statistics - galaxies: evolution 7 pages, 3 figures

We examine the properties of galaxies in the Galaxies and Mass Assembly (GAMA) survey located in voids with radii $>10~h^{-1}$ Mpc. Utilising the GAMA equatorial survey, 592 void galaxies are identified out to z~0.1 brighter than $M_{r} = -18.4$, our magnitude completeness limit. Using the $W_{\rm{H\alpha}}$ vs. [NII]/H$\alpha$ (WHAN) line strength diagnostic diagram, we classify their spectra as star forming, AGN, or dominated by old stellar populations. For objects more massive than $5\times10^{9}$ M$_{\odot}$, we identify a sample of 26 void galaxies with old stellar populations classed as passive and retired galaxies in the WHAN diagnostic diagram, else they lack any emission lines in their spectra. When matched to WISE mid-IR photometry, these passive and retired galaxies exhibit a range of mid-IR colour, with a number of void galaxies exhibiting [4.6]-[12] colours inconsistent with completely quenched stellar populations, with a similar spread in colour seen for a randomly drawn non-void comparison sample. We hypothesise that a number of these galaxies host obscured star formation, else they are star forming outside of their central regions targeted for single fibre spectroscopy. When matched to a randomly drawn sample of non-void galaxies, the void and non-void galaxies exhibit similar properties in terms of optical and mid-IR colour, morphology, and star formation activity, suggesting comparable mass assembly and quenching histories. A trend in mid-IR [4.6]-[12] colour is seen, such that both void and non-void galaxies with quenched/passive colours <1.5 typically have masses higher than $10^{10}$ M$_{\odot}$, where internally driven processes play an increasingly important role in galaxy evolution., Comment: 22 pages, 12 figures. Accepted for publication in MNRAS

This article is the second in a series in which we perform an extensive comparison of various galaxy-based cluster mass estimation techniques that utilise the positions, velocities and colours of galaxies. Our aim is to quantify the scatter, systematic bias and completeness of cluster masses derived from a diverse set of 25 galaxy-based methods using two contrasting mock galaxy catalogues based on a sophisticated halo occupation model and a semi-analytic model. Analysing 968 clusters, we find a wide range in the RMS errors in log M200c delivered by the different methods (0.18 to 1.08 dex, i.e., a factor of ~1.5 to 12), with abundance matching and richness methods providing the best results, irrespective of the input model assumptions. In addition, certain methods produce a significant number of catastrophic cases where the mass is under- or over-estimated by a factor greater than 10. Given the steeply falling high-mass end of the cluster mass function, we recommend that richness or abundance matching-based methods are used in conjunction with these methods as a sanity check for studies selecting high mass clusters. We see a stronger correlation of the recovered to input number of galaxies for both catalogues in comparison with the group/cluster mass, however, this does not guarantee that the correct member galaxies are being selected. We do not observe significantly higher scatter for either mock galaxy catalogues. Our results have implications for cosmological analyses that utilise the masses, richnesses, or abundances of clusters, which have different uncertainties when different methods are used., Comment: 25 pages, 19 figures, 7 tables. Accepted for publication in MNRAS

We study the relationship between the structure and star-formation rate (SFR) of X-ray selected low and moderate luminosity active galactic nuclei (AGNs) in the two Chandra Deep Fields, using Hubble Space Telescope imaging from the Cosmic Assembly Near Infrared Extragalactic Legacy Survey (CANDELS) and deep far-infrared maps from the PEP+GOODS-Herschel survey. We derive detailed distributions of structural parameters and FIR luminosities from carefully constructed control samples of galaxies, which we then compare to those of the AGNs. At z~1, AGNs show slightly diskier light profiles than massive inactive (non-AGN) galaxies, as well as modestly higher levels of gross galaxy disturbance (as measured by visual signatures of interactions and clumpy structure). In contrast, at z~2, AGNs show similar levels of galaxy disturbance as inactive galaxies, but display a red central light enhancement, which may arise due to a more pronounced bulge in AGN hosts or due to extinguished nuclear light. We undertake a number of tests of these alternatives, but our results do not strongly favour one interpretation over the other. The mean SFR and its distribution among AGNs and inactive galaxies are similar at z>1.5. At z<1, however, clear and significant enhancements are seen in the SFRs of AGNs with bulge-dominated light profiles. These trends suggest an evolution in the relation between nuclear activity and host properties with redshift, towards a minor role for mergers and interactions at z>1.5., Comment: 25 pages (double-column), 17 Figures, under review at A&A

The SAMI Galaxy Survey will observe 3400 galaxies with the Sydney-AAO Multi-object Integral-field spectrograph (SAMI) on the Anglo-Australian Telescope (AAT) in a 3-year survey which began in 2013. We present the throughput of the SAMI system, the science basis and specifications for the target selection, the survey observation plan and the combined properties of the selected galaxies. The survey includes four volume limited galaxy samples based on cuts in a proxy for stellar mass, along with low-stellar mass dwarf galaxies all selected from the Galaxy And Mass Assembly (GAMA) survey. The GAMA regions were selected because of the vast array of ancillary data available, including ultraviolet through to radio bands. These fields are on the celestial equator at 9, 12, and 14.5 hours, and cover a total of 144 square degrees (in GAMA-I). Higher density environments are also included with the addition of eight clusters. The clusters have spectroscopy from 2dFGRS and SDSS and photometry in regions covered by the Sloan Digital Sky Survey (SDSS) and/or VLT Survey Telescope/ATLAS. The aim is to cover a broad range in stellar mass and environment, and therefore the primary survey targets cover redshifts 0.004 < z < 0.095, magnitudes r$_{pet}$ < 19.4, stellar masses $10^{7} - 10^{12}$ M$_{sol}$, and environments from isolated field galaxies through groups to clusters of $10^{15}$ M$_{sol}$., Comment: 19 figures; Published in MNRAS 447, 2857, 2015

Recent observations have probed the formation histories of nearby elliptical galaxies by tracking correlations between the stellar population parameters, age and metallicity, and the structural parameters that enter the Fundamental Plane, radius and velocity dispersion sigma. These studies have found intriguing correlations between these four parameters. In this work, we make use of a semi-analytic model, based on halo merger trees extracted from the Bolshoi cosmological simulation, that predicts the structural properties of spheroid-dominated galaxies based on an analytic model that has been tested and calibrated against an extensive suite of hydrodynamic+N-body binary merger simulations. We predict the radius, sigma, luminosity, age, and metallicity of spheroid-dominated galaxies, enabling us to compare directly to observations. Our model predicts a strong correlation between age and sigma for early-type galaxies, and no significant correlation between age and radius, in agreement with observations. In addition we predict a strong correlation between metallicity and sigma, and a weak correlation between metallicity and radius, in qualitative agreement with observations. We find that the correlations with sigma arise as a result of the strong link between sigma and the galaxy assembly time. Minor mergers produce a large change in radius while leaving sigma nearly the same, which explains the weaker trends with radius., Comment: 15 pages, 9 figures, submitted to MNRAS. arXiv admin note: text overlap with arXiv:1201.5918

This paper is the first in a series in which we perform an extensive comparison of various galaxy-based cluster mass estimation techniques that utilise the positions, velocities and colours of galaxies. Our primary aim is to test the performance of these cluster mass estimation techniques on a diverse set of models that will increase in complexity. We begin by providing participating methods with data from a simple model that delivers idealised clusters, enabling us to quantify the underlying scatter intrinsic to these mass estimation techniques. The mock catalogue is based on a Halo Occupation Distribution (HOD) model that assumes spherical Navarro, Frenk and White (NFW) haloes truncated at R_200, with no substructure nor colour segregation, and with isotropic, isothermal Maxwellian velocities. We find that, above 10^14 M_solar, recovered cluster masses are correlated with the true underlying cluster mass with an intrinsic scatter of typically a factor of two. Below 10^14 M_solar, the scatter rises as the number of member galaxies drops and rapidly approaches an order of magnitude. We find that richness-based methods deliver the lowest scatter, but it is not clear whether such accuracy may simply be the result of using an over-simplistic model to populate the galaxies in their haloes. Even when given the true cluster membership, large scatter is observed for the majority non-richness-based approaches, suggesting that mass reconstruction with a low number of dynamical tracers is inherently problematic., Comment: 25 pages, 15 figures, 5 tables. Accepted for publication in MNRAS

The processes that trigger Active Galactic Nuclei (AGN) remain poorly understood. While lower luminosity AGN may be triggered by minor disturbances to the host galaxy, stronger disturbances are likely required to trigger luminous AGN. Major wet mergers of galaxies are ideal environments for AGN triggering since they provide large gas supplies and galaxy scale torques. There is however little observational evidence for a strong connection between AGN and major mergers. We analyse the morphological properties of AGN host galaxies as a function of AGN and host galaxy luminosity and compare them to a carefully matched sample of control galaxies. AGN are X-ray selected in the redshift range 0.5 < z < 0.8 and have luminosities 41 < log(L_X [erg/s]) < 44.5. 'Fake AGN' are simulated in the control galaxies by adding point sources with the magnitude of the matched AGN. We find that AGN host and control galaxies have comparable assymetries, Sersic indices and ellipticities at restframe ~950nm. AGN host galaxies show neither higher average asymmetries nor higher fractions of very disturbed objects. There is no increase in the prevalence of merger signatures with AGN luminosity. At 95% confidence we find that major mergers are responsible for <6% of all AGN in our sample as well as <40% of the highest luminosity AGN log(L_X [erg/s]) ~ 43.5). Major mergers therefore either play only a very minor role in the triggering of AGN in the luminosity range studied or time delays are too long for merger features to remain visible., Comment: Accepted for publication in MNRAS

We analyze the star-forming and structural properties of 45 massive (log(M/Msun)>10) compact star-forming galaxies (SFGs) at 2

We have used high-resolution, HST WFC3/IR, near-infrared imaging to conduct a detailed bulge-disk decomposition of the morphologies of ~200 of the most massive (M_star > 10^11 M_solar) galaxies at 12 they are mostly disk-dominated. Interestingly, we find that while most of the quiescent galaxies are bulge-dominated, a significant fraction (25-40%) of the most quiescent galaxies, have disk-dominated morphologies. Thus, our results suggest that the physical mechanisms which quench star-formation activity are not simply connected to those responsible for the morphological transformation of massive galaxies., Comment: Contributed article (4 pages, 2 figures), IAU Symposium 295 on "The intriguing life of massive galaxies", to appear in the proceedings of the IAU XXVIII General Assembly, Beijing, 2012, (Eds: Daniel Thomas, Anna Pasquali, Ignacio Ferreras)

We have used deep, HST, near-IR imaging to study the morphological properties of the most massive galaxies at high z, modelling the WFC3/IR H-band images of the ~200 galaxies in the CANDELS-UDS field with 1 < z_phot < 3, and stellar masses M_star > 10^11 M_sun. We have used both single-Sersic and bulge+disk models, have investigated the errors/biases introduced by uncertainties in the background and the PSF, and have obtained formally-acceptable model fits to >90% of the galaxies. Our results indicate that these massive galaxies at 1 < z < 3 lie both on and below the local size-mass relation, with a median R_e~2.6 kpc, a factor of ~2.3 smaller than comparably-massive local galaxies. Moreover, we find that bulge-dominated objects in particular show evidence for a growing bimodality in the size-mass relation with increasing z, and by z > 2 the compact bulges display effective radii a factor ~4 smaller than local ellipticals of comparable mass. These trends appear to extend to the bulge components of disk-dominated galaxies, and vice versa. We also find that, while such massive galaxies at low z are bulge-dominated, at 1 < z < 2 they are predominantly mixed bulge+disk systems, and by z > 2 they are mostly disk-dominated. The majority of the disk-dominated galaxies are actively forming stars, but this is also true for many of the bulge-dominated systems. Interestingly, however, while most of the quiescent galaxies are bulge-dominated, we find that a significant fraction (25-40%) of the most quiescent galaxies have disk-dominated morphologies. Thus, while our results show that the massive galaxy population is undergoing dramatic changes at this crucial epoch, they also suggest that the physical mechanisms which quench star-formation activity are not simply connected to those responsible for the morphological transformation of massive galaxies into present-day giant ellipticals., Comment: 39 pages, 19 figures, 2 tables, submitted to MNRAS

Recent low-redshift observations have attempted to determine the star formation histories of elliptical galaxies by tracking correlations between the stellar population parameters (age and metallicity) and the structural parameters that enter the fundamental plane (size and velocity dispersion). These studies have found that velocity dispersion, rather than effective radius or dynamical mass, is the main predictor of a galaxy's stellar age and metallicity. In this work, we apply an analytic model that predicts the structural properties of remnants formed in major mergers to progenitor disk galaxies with properties taken from two different semi-analytic models. We predict the effective radius, velocity dispersion, luminosity, age, and metallicity of the merger remnants, enabling us to compare directly to observations of early-type galaxies. While we find a tight correlation between age and velocity dispersion, we find a stronger dependence of age and metallicity on effective radius than observations report. The correlations arise as a result of the dependence of gas fraction, age, and metallicity on the stellar mass in the progenitor disk galaxies. These dependences induce a rotation in the radius-velocity plane between the correlations with effective radius and circular velocity in the disk galaxy progenitors, and the correlations with effective radius and velocity dispersion in the elliptical galaxy remnants. The differences between our results and those from observations suggest that major mergers alone cannot produce the observed lack of correlation between effective radius and stellar population parameters. Simulations have suggested that subsequent minor mergers introduce scatter in the effective radius while leaving the velocity dispersion essentially unchanged. Incorporating such minor mergers into the model may, then, bring the simulations into closer agreement with observations., Comment: 20 pages, 17 figures. Submitted to MNRAS

There are strong correlations between the three structural properties of elliptical galaxies -- stellar mass, velocity dispersion and size -- in the form of a tight "fundamental plane" and a "scaling relation" between each pair. Major mergers of disk galaxies are assumed to be a mechanism for producing ellipticals, but semi-analytic galaxy formation models (SAM) have encountered apparent difficulties in reproducing the observed slope and scatter of the size-mass relation. We study the scaling relations of merger remnants using progenitor properties from two SAMs. We apply a simple merger model that includes gas dissipation and star formation based on theoretical considerations and simulations. Combining the SAMs and the merger model allows calculation of the structural properties of the remnants of major mergers that enter the population of elliptical galaxies at a given redshift. Without tuning the merger model parameters for each SAM, the results roughly match the slope and scatter in the observed scaling relations and their evolution in the redshift range $z=0-3$. Within this model, the observed scaling relations, including the tilt of the fundamental plane relative to the virial plane, result primarily from the decrease of gas fraction with increasing progenitor mass. The scatter in the size-mass relation of the remnants is reduced from that of the progenitors because of a correlation between progenitor size and gas fraction at a given mass., Comment: in press

We combine SAURON integral field data of a representative sample of local early-type, red sequence galaxies with Spitzer/IRAC imaging in order to investigate the presence of trace star formation in these systems. With the Spitzer data, we identify galaxies hosting low-level star formation, as traced by PAH emission, with measured star formation rates that compare well to those estimated from other tracers. This star formation proceeds according to established scaling relations with molecular gas content, in surface density regimes characteristic of disk galaxies and circumnuclear starbursts. We find that star formation in early-type galaxies happens exclusively in fast-rotating systems and occurs in two distinct modes. In the first, star formation is a diffuse process, corresponding to widespread young stellar populations and high molecular gas content. The equal presence of co- and counter-rotating components in these systems strongly implies an external origin for the star-forming gas, and we argue that these star formation events may be the final stages of (mostly minor) mergers that build up the bulges of red sequence lenticulars. In the second mode of star formation, the process is concentrated into well-defined disk or ring morphologies, outside of which the host galaxies exhibit uniformly evolved stellar populations. This implies that these star formation events represent rejuvenations within previously quiescent stellar systems. Evidence for earlier star formation events similar to these in all fast rotating early-type galaxies suggests that this mode of star formation may be common to all such galaxies, with a duty cycle of roughly 1/10, and likely contributes to the embedded, co-rotating inner stellar disks ubiquitous in this population., Comment: Accepted for publication in MNRAS. 48 pages, 36 figures. Revision fixes typos, including column headers in Table 2. The revised paper with full resolution figures is available at http://astro.berkeley.edu/~shapiro/Papers/Shapiro_SAURON_PaperXV.pdf

We present a simulation of the cosmic evolution of the atomic and molecular phases of the cold hydrogen gas in about 3e7 galaxies, obtained by post-processing the virtual galaxy catalog produced by (De Lucia et al. 2007) on the Millennium Simulation of cosmic structure (Springel et al. 2005). Our method uses a set of physical prescriptions to assign neutral atomic hydrogen (HI) and molecular hydrogen (H2) to galaxies, based on their total cold gas masses and a few additional galaxy properties. These prescriptions are specially designed for large cosmological simulations, where, given current computational limitations, individual galaxies can only be represented by simplistic model-objects with a few global properties. Our recipes allow us to (i) split total cold gas masses between HI, H2, and Helium, (ii) assign realistic sizes to both the HI- and H2-disks, and (iii) evaluate the corresponding velocity profiles and shapes of the characteristic radio emission lines. The results presented in this paper include the local HI- and H2-mass functions, the CO-luminosity function, the cold gas mass--diameter relation, and the Tully-Fisher relation (TFR), which all match recent observational data from the local Universe. We also present high-redshift predictions of cold gas diameters and the TFR, both of which appear to evolve markedly with redshift., Comment: 19 pages, 15 figures

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

We use a semianalytic galaxy catalogue constructed from the Millennium Simulation to study the satellites of isolated galaxies in the LCDM cosmogony. This sample (~80,000$ bright primaries, surrounded by ~178,000 satellites) allows the characterization, with minimal statistical uncertainty, of the dynamical properties of satellite/primary galaxy systems in a LCDM universe. We find that, overall, the satellite population traces the dark matter rather well: its spatial distribution and kinematics may be approximated by an NFW profile with a mildly anisotropic velocity distribution. Their spatial distribution is also mildly anisotropic, with a well-defined ``anti-Holmberg'' effect that reflects the misalignment between the major axis and angular momentum of the host halo. The isolation criteria for our primaries picks not only galaxies in sparse environments, but also a number of primaries at the centre of ''fossil'' groups. We find that the abundance and luminosity function of these unusual systems are in reasonable agreement with the few available observational constraints. We recover the expected L_{host} \sigma_{sat}^3 relation for LCDM models for truly-isolated primaries. Less strict primary selection, however, leads to substantial modification of the scaling relation. Our analysis also highlights a number of difficulties afflicting studies that rely on blind stacking of satellite systems to constrain the mean halo mass of the primary galaxies., Comment: 18 pages, 14 figures, MNRAS in press. Accepted version with minor changes. Version with high resolution figures available at: http://www.astro.uvic.ca/~lsales/SatPapers/SatPapers.html

We investigate halo mass selection properties of red-sequence cluster finders using galaxy populations of the Millennium Simulation (MS). A clear red sequence exists for MS galaxies in massive halos at redshifts z < 1, and we use this knowledge to inform a cluster-finding algorithm applied to 500 Mpc/h projections of the simulated volume. At low redshift (z=0.4), we find that 90% of the clusters found have galaxy membership dominated by a single, real-space halo, and that 10% are blended systems for which no single halo contributes a majority of a cluster's membership. At z=1, the fraction of blends increases to 22%, as weaker redshift evolution in observed color extends the comoving length probed by a fixed range of color. Other factors contributing to the increased blending at high-z include broadening of the red sequence and confusion from a larger number of intermediate mass halos hosting bright red galaxies of magnitude similar to those in higher mass halos. Our method produces catalogs of cluster candidates whose halo mass selection function, p(M|\Ngal,z), is characterized by a bimodal log-normal model with a dominant component that reproduces well the real-space distribution, and a redshift-dependent tail that is broader and displaced by a factor ~2 lower in mass. We discuss implications for X-ray properties of optically selected clusters and offer ideas for improving both mock catalogs and cluster-finding in future surveys., Comment: final version to appear in MNRAS. Appendix added on purity and completeness, small shift in red sequence due to correcting an error in finding it

We present a 200 ksec Chandra observation of seven spectroscopically selected, high redshift (0.75 < z < 1.03) galaxy groups and clusters discovered by the DEEP2 Galaxy Redshift Survey in the Extended Groth Strip (EGS). X-ray emission at the locations of these systems is consistent with background. The 3-sigma upper limits on the bolometric X-ray luminosities (L_X) of these systems put a strong constraint on the relation between L_X and the velocity dispersion of member galaxies sigma_gal at z~1; the DEEP2 systems have lower luminosity than would be predicted by the local relation. Our result is consistent with recent findings that at high redshift, optically selected clusters tend to be X-ray underluminous. A comparison with mock catalogs indicates that it is unlikely that this effect is entirely caused by a measurement bias between sigma_gal and the dark matter velocity dispersion. Physically, the DEEP2 systems may still be in the process of forming and hence not fully virialized, or they may be deficient in hot gas compared to local systems. We find only one possibly extended source in this Chandra field, which happens to lie outside the DEEP2 coverage., Comment: 5 pages, 3 figures. Accepted for publication in AEGIS ApJ Letters special edition

In this the first of a series of Letters, we present a description of the panchromatic data sets that have been acquired in the Extended Groth Strip region of the sky. Our survey, the All-wavelength Extended Groth strip International Survey (AEGIS), is intended to study the physical properties and evolutionary processes of galaxies at z ~ 1. It includes the following deep, wide-field imaging data sets: Chandra/ACIS X-ray (0.5 - 10 keV), GALEX ultraviolet (1200 - 2500 Angstrom), CFHT/MegaCam Legacy Survey optical (3600 - 9000 Angstroms), CFHT/CFH12K optical (4500 - 9000 Angstroms), Hubble Space Telescope/ACS optical (4400 - 8500 Angstroms), Palomar/WIRC near-infrared (1.2 - 2.2 microns), Spitzer/IRAC mid-infrared (3.6 - 8.0 microns), Spitzer/MIPS far-infrared (24 - 70 microns), and VLA radio continuum (6 - 20 cm). In addition, this region of the sky has been targeted for extensive spectroscopy using the DEIMOS spectrograph on the Keck II 10 m telescope. Our survey is compared to other large multiwavelength surveys in terms of depth and sky coverage., Comment: 21 pages, 1 fig; intro letter for AEGIS ApJL special issue; ps & pdf at http://www.ucolick.org/~raja/aegis.ps & http://www.ucolick.org/~raja/aegis.pdf

We explore the environment of z~1 AGN using a sample of 53 spectroscopically identified X-ray sources in the All-wavelength Extended Groth strip International Survey. We quantify the local density in the vicinity of an X-ray source by measuring the projected surface density of spectroscopically identified optical galaxies within a radius defined by the 3rd nearest neighbour. Our main result is that X-ray selected AGN at z~1 avoid underdense regions at the 99.89% confidence level. Moreover, although we find that the overall population shares the same (rich) environment with optical galaxies of similar U-B and M_B, there is also tentative evidence (96%) that AGN with blue colors (U-B<1) reside in denser environments compared to optical galaxies. We argue that the results above are a consequence of the whereabouts of massive galaxies, capable of hosting supermassive black holes at their centers, with available cold gas reservoirs, the fuel for AGN activity. At z~1 an increasing fraction of such systems are found in dense regions., Comment: accepted by ApJL

We discuss the relationship between rest-frame color and optical luminosity for X-ray sources in the range 0.6

We present new results for the 3-point correlation function, \zeta, measured as a function of scale, luminosity and colour from the final version of the two-degree field galaxy redshift survey (2dFGRS). The reduced three point correlation function, Q_3 is estimated for different triangle shapes and sizes, employing a full covariance analysis. The form of Q_3 is consistent with the expectations for the \Lambda-cold dark matter model, confirming that the primary influence shaping the distribution of galaxies is gravitational instability acting on Gaussian primordial fluctuations. However, we find a clear offset in amplitude between Q_3 for galaxies and the predictions for the dark matter. We are able to rule out the scenario in which galaxies are unbiased tracers of the mass at the 9-sigma level. On weakly non-linear scales, we can interpret our results in terms of galaxy bias parameters. We find a linear bias term that is consistent with unity, b_1 = 0.93^{+0.10}_{-0.08} and a quadratic bias c_2 = b_2 /b_1 = -0.34^{+0.11}_{-0.08}. This is the first significant detection of a non-zero quadratic bias, indicating a small but important non-gravitational contribution to the three point function. Our estimate of the linear bias from the three point function is independent of the normalisation of underlying density fluctuations, so we can combine this with the measurement of the power spectrum of 2dFGRS galaxies to constrain the amplitude of matter fluctuations. We find that the rms linear theory variance in spheres of radius 8Mpc/h is \sigma_8 = 0.88^{+0.12}_{-0.10}, providing an independent confirmation of values derived from other techniques. On non-linear scales, where \xi>1, we find that Q_3 has a strong dependence on scale, colour and luminosity., Comment: 16 pages, 10 figures, minor changes, extended comparison to previous results, accepted for publication in MNRAS

We measure the clustering of galaxy groups in the 2dFGRS Percolation-Inferred Galaxy Group (2PIGG) catalogue. The 2PIGG sample has 29,000 groups with at least two members. The clustering amplitude of the full 2PIGG catalogue is weaker than that of 2dFGRS galaxies, in agreement with theoretical predictions. We have subdivided the 2PIGG catalogue into samples that span a factor of 25 in median total luminosity. Our correlation function measurements span an unprecedented range of clustering strengths, connecting the regimes probed by groups fainter than L* galaxies and rich clusters. There is a steady increase in clustering strength with group luminosity; the most luminous groups are ten times more strongly clustered than the full 2PIGG catalogue. We demonstrate that the 2PIGG results are in very good agreement with the clustering of groups expected in the LCDM model., Comment: 16 pages. Submitted to MNRAS

We measure moments of the galaxy count probability distribution function in the two-degree field galaxy redshift survey (2dFGRS). The survey is divided into volume limited subsamples in order to examine the dependence of the higher order clustering on galaxy luminosity. We demonstrate the hierarchical scaling of the averaged p-point galaxy correlation functions, xibar_p, up to p=6. The hierarchical amplitudes, S_p = xibar_p/xibar_2^{p-1}, are approximately independent of the cell radius used to smooth the galaxy distribution on small to medium scales. On larger scales we find the higher order moments can be strongly affected by the presence of rare, massive superstructures in the galaxy distribution. The skewness S_3 has a weak dependence on luminosity, approximated by a linear dependence on log luminosity. We discuss the implications of our results for simple models of linear and non-linear bias that relate the galaxy distribution to the underlying mass., Comment: 14 pages, 10 figures. Slight revisions to match version published in MNRAS

We use the two-degree field Galaxy Redshift Survey (2dFGRS) to test the hierarchical scaling hypothesis: namely, that the $p$-point galaxy correlation functions can be written in terms of the two point correlation function or variance. This scaling is expected if an initially Gaussian distribution of density fluctuations evolves under the action of gravitational instability. We measure the volume averaged $p$-point correlation functions using a counts in cells technique applied to a volume limited sample of 44,931 $L_*$ galaxies. We demonstrate that $L_{*}$ galaxies display hierarchical clustering up to order $p=6$ in redshift space. The variance measured for $L_{*}$ galaxies is in excellent agreement with the predictions from a $\Lambda$-cold dark matter N-body simulation. This applies to all cell radii considered, $0.3<(R/h^{-1}{\rm Mpc})<30$. However, the higher order correlation functions of $L_*$ galaxies have a significantly smaller amplitude than is predicted for the dark matter for $R<10h^{-1}$Mpc. This disagreement implies that a non-linear bias exists between the dark matter and $L_*$ galaxies on these scales. We also show that the presence of two rare, massive superclusters in the 2dFGRS has an impact on the higher-order clustering moments measured on large scales., Comment: 5 pages, 4 figures. Slight revisions to match version published in MNRAS

We measure the redshift space reduced void probability function (VPF) for 2dFGRS volume limited galaxy samples covering the absolute magnitude range M_bJ-5logh=-18 to -22. Theoretically, the VPF connects the distribution of voids to the moments of galaxy clustering of all orders, and can be used to discriminate clustering models in the weakly non-linear regime. The reduced VPF measured from the 2dFGRS is in excellent agreement with the paradigm of hierarchical scaling of the galaxy clustering moments. The accuracy of our measurement is such that we can rule out, at a very high significance, popular models for galaxy clustering, including the lognormal distribution. We demonstrate that the negative binomial model gives a very good approximation to the 2dFGRS data over a wide range of scales, out to at least 20h-1Mpc. Conversely, the reduced VPF for dark matter in a LambdaCDM universe does appear to be lognormal on small scales but deviates significantly beyond \approx 4h-1Mpc. We find little dependence of the 2dFGRS reduced VPF on galaxy luminosity. Our results hold independently in both the north and south Galactic pole survey regions., Comment: 10 pages, 6 figures. Dark matter comparison added. Accepted MNRAS

The Spitzer Extended Deep Survey (SEDS) is a very deep infrared survey within five well-known extragalactic science fields: the UKIDSS Ultra-Deep Survey, the Extended Chandra Deep Field South, COSMOS, the Hubble Deep Field North, and the Extended Groth Strip. SEDS covers a total area of 1.46 deg(exp 2) to a depth of 26 AB mag (3sigma) in both of the warm Infrared Array Camera (IRAC) bands at 3.6 and 4.5 micron. Because of its uniform depth of coverage in so many widely-separated fields, SEDS is subject to roughly 25% smaller errors due to cosmic variance than a single-field survey of the same size. SEDS was designed to detect and characterize galaxies from intermediate to high redshifts (z = 2-7) with a built-in means of assessing the impact of cosmic variance on the individual fields. Because the full SEDS depth was accumulated in at least three separate visits to each field, typically with six-month intervals between visits, SEDS also furnishes an opportunity to assess the infrared variability of faint objects. This paper describes the SEDS survey design, processing, and publicly-available data products. Deep IRAC counts for the more than 300,000 galaxies detected by SEDS are consistent with models based on known galaxy populations. Discrete IRAC sources contribute 5.6 +/- 1.0 and 4.4 +/- 0.8 nW / square m/sr at 3.6 and 4.5 micron to the diffuse cosmic infrared background (CIB). IRAC sources cannot contribute more than half of the total CIB flux estimated from DIRBE data. Barring an unexpected error in the DIRBE flux estimates, half the CIB flux must therefore come from a diffuse component.