Time Evolution of Galaxy Formation and Bias in Cosmological Simulations

The clustering of galaxies relative to the mass distribution declines with time because: first, nonlinear peaks become less rare events; second, the densest regions stop forming new galaxies because gas there becomes too hot to cool and collapse; third, after galaxies form, they are gravitationally ``debiased'' because their velocity field is the same as the dark matter. To show these effects, we perform a hydrodynamic cosmological simulation and examine the density field of recently formed galaxies as a function of redshift. We find the bias b_* of recently formed galaxies (the ratio of the rms fluctuations of these galaxies and mass), evolves from 4.5 at z=3 to around 1 at z=0, on 8 h^{-1} Mpc comoving scales. The correlation coefficient r_* between recently formed galaxies and mass evolves from 0.9 at z=3 to 0.25 at z=0. As gas in the universe heats up and prevents star formation, star-forming galaxies become poorer tracers of the mass density field. After galaxies form, the linear continuity equation is a good approximation to the gravitational debiasing, even on nonlinear scales. The most interesting observational consequence of the simulations is that the linear regression of the star-formation density field on the galaxy density field evolves from about 0.9 at z=1 to 0.35 at z=0. These effects also provide a possible explanation for the Butcher-Oemler effect, the excess of blue galaxies in clusters at redshift z ~ 0.5. Finally, we examine cluster mass-to-light ratio estimates of Omega, finding that while Omega(z) increases with z, one's estimate Omega_est(z) decreases. (Abridged)Comment: 31 pages of text and figures; submitted to Ap

The Faint End of the Luminosity Function and Low Surface Brightness Galaxies

SHELS (Smithsonian Hectospec Lensing Survey) is a dense redshift survey covering a 4 square degree region to a limiting R = 20.6. In the construction of the galaxy catalog and in the acquisition of spectroscopic targets, we paid careful attention to the survey completeness for lower surface brightness dwarf galaxies. Thus, although the survey covers a small area, it is a robust basis for computation of the slope of the faint end of the galaxy luminosity function to a limiting M_R = -13.3 + 5logh. We calculate the faint end slope in the R-band for the subset of SHELS galaxies with redshif ts in the range 0.02 <= z < 0.1, SHELS_{0.1}. This sample contains 532 galaxies with R< 20.6 and with a median surface brightness within the half light radius of SB_{50,R} = 21.82 mag arcsec^{-2}. We used this sample to make one of the few direct measurements of the dependence of the faint end of the galaxy luminosity function on surface brightness. For the sample as a whole the faint end slope, alpha = -1.31 +/- 0.04, is consistent with both the Blanton et al. (2005b) analysis of the SDSS and the Liu et al. (2008) analysis of the COSMOS field. This consistency is impressive given the very different approaches of th ese three surveys. A magnitude limited sample of 135 galaxies with optical spectroscopic reds hifts with mean half-light surface brightness, SB_{50,R} >= 22.5 mag arcsec^{-2} is unique to SHELS_{0.1}. The faint end slope is alpha_{22.5} = -1.52+/- 0.16. SHELS_{0.1} shows that lower surface brightness objects dominate the faint end slope of the l uminosity function in the field, underscoring the importance of surface brightness limits in evaluating measurements of the faint end slope and its evolution.Comment: 34 pages, 13 figures, 3 tables, Astronomical Journal, in press (updated based on review

The Millennium Galaxy Catalogue: the space density and surface brightness distribution(s) of galaxies

We recover the joint and individual space density and surface brightness distribution(s) of galaxies from the Millennium Galaxy Catalogue. The MGC is a local survey spanning 30.9 sq deg and probing approximately one--two mag/sq arcsec deeper than either the Two-Degree Field Galaxy Redshift Survey (2dFGRS) or the Sloan Digital Sky Survey (SDSS). The MGC contains 10,095 galaxies to B_mgc < 20 mag with 96 per cent spectroscopic completeness. We implement a joint luminosity-surface brightness step-wise maximum likelihood method to recover the bivariate brightness distribution (BBD) inclusive of most selection effects. Integrating the BBD over surface brightness we recover the following Schechter function parameters: phi* = (0.0177 +/- 0.0015) h^3 Mpc^{-3}, M_{B}* - 5 log h = (-19.60 +/- 0.04) mag and alpha =-1.13 +/- 0.02. Compared to the 2dFGRS (Norberg et al 2002) we find a consistent M* value but a slightly flatter faint-end slope and a higher normalisation, resulting in a final luminosity density j_{b_J} = (1.99 +/- 0.17) x 10^8 h L_{odot} Mpc^{-3}. The MGC surface brightness distribution is a well bounded Gaussian at the M* point with phi* = (3.5 +/- 0.1) x 10^{-2} h^3 Mpc^{-3}, mu^{e*} = (21.90 +/- 0.01) mag/sq arcsec and sigma_{ln R_e} = 0.35 +/- 0.01. The characteristic surface brightness for luminous systems is invariant to M_{B} - 5 log h ~ -19 mag faintwards of which it moves to lower surface brightness. Higher resolution (FWHM 26 mag/sq arcsec in the B-band) observations of the local universe are now essential to probe to lower luminosity and lower surface brightness levels. [abridged]Comment: Accepted for publication in MNRAS, 26 pages with 21 figures (some degraded). A full pdf version, along with MGC data release, is available from the MGC website at, http://www.eso.org/~jliske/mg

Influences of physical oceanographic processes on chlorophyll distributions in coastal and estuarine waters of the South Atlantic Bight

Coastal and estuarine waters of the South Atlantic Bight are highly productive, with primary production of 600-700 gC/m2/y. While controls and fate of this production are conceptually well understood, the importance of meteorology and physical circulation processes on phytoplankton has not received equivalent attention. Here, we describe the effects of wind stress and tidal currents on temporal and spatial distributions of phytoplankton biomass represented as chlorophyll a (chl a). Moored instruments were deployed and shipboard sampling was conducted in the North Edisto estuary (South Carolina) and adjacent inner shelf waters during four, two-week field studies in May and August 1993, and June and September 1994. Local wind regimes induced upwelling- and downwelling-favorable conditions which strengthened or reduced vertical density stratification in the coastal frontal zone, respectively, and shifted the location of the front. Chl a in shelf waters was more or less homogenous independent of the wind regime, while chl a on the estuary delta was generally vertically stratified. Within the estuary, chl a concentrations were positively correlated with the alongshore component of wind stress; chl a was not correlated with the weaker cross-shelf component of wind stress. Highest chl a occurred during strong downwelling-favorable events. The quick response time to wind forcing (6-12 hrs) implied a direct effect on chl a distributions and not a stimulation of growth processes. The source of the elevated chl a in response to wind forcing was apparently resuspension of settled and epibenthic algal cells. Tidal currents also influenced the vertical distribution and concentration of chl a. Time series sampling on the estuary delta showed that, with increasing velocity of ebb and flood tide currents, the relative contributions of pennate and centric diatoms with attached detritus and sand grains also increased, indicating that tidal resuspension of settled and epibenthic microalgae also occurred. Vertical stratification of chl a (highest concentrations near the bottom) began to degrade upon mixing by tidal currents with velocities as low as 10 cm/sec. Homogenization of 5-7 m water columns was fully achieved at velocities of 20-30 cm/sec. The data document the direct and comparatively immediate (timescales of minuteshours) impact of tidal and wind energy on concentrations and distribution patterns of phytoplankton in coastal and estuarine waters of the South Atlantic Bight

The UV-Optical Color Dependence of Galaxy Clustering in the Local Universe

We measure the UV-optical color dependence of galaxy clustering in the local universe. Using the clean separation of the red and blue sequences made possible by the NUV - r color-magnitude diagram, we segregate the galaxies into red, blue and intermediate "green" classes. We explore the clustering as a function of this segregation by removing the dependence on luminosity and by excluding edge-on galaxies as a means of a non-model dependent veto of highly extincted galaxies. We find that \xi (r_p, \pi) for both red and green galaxies shows strong redshift space distortion on small scales -- the "finger-of-God" effect, with green galaxies having a lower amplitude than is seen for the red sequence, and the blue sequence showing almost no distortion. On large scales, \xi (r_p, \pi) for all three samples show the effect of large-scale streaming from coherent infall. On scales 1 Mpc/h < r_p < 10 Mpc/h, the projected auto-correlation function w_p(r_p) for red and green galaxies fits a power-law with slope \gamma ~ 1.93 and amplitude r_0 ~ 7.5 and 5.3, compared with \gamma ~ 1.75 and r_0 ~ 3.9 Mpc/h for blue sequence galaxies. Compared to the clustering of a fiducial L* galaxy, the red, green, and blue have a relative bias of 1.5, 1.1, and 0.9 respectively. The w_p(r_p) for blue galaxies display an increase in convexity at ~ 1 Mpc/h, with an excess of large scale clustering. Our results suggest that the majority of blue galaxies are likely central galaxies in less massive halos, while red and green galaxies have larger satellite fractions, and preferentially reside in virialized structures. If blue sequence galaxies migrate to the red sequence via processes like mergers or quenching that take them through the green valley, such a transformation may be accompanied by a change in environment in addition to any change in luminosity and color.Comment: accepted by MNRA

The overdensities of galaxy environments as a function of luminosity and color

We study the mean environments of galaxies in the Sloan Digital Sky Survey as a function of rest-frame luminosity and color. Overdensities in galaxy number are estimated in $8 h^{-1} \mathrm{Mpc}$ and $1 h^{-1} \mathrm{Mpc}$ spheres centered on $125,000$ galaxies taken from the SDSS spectroscopic sample. We find that, at constant color, overdensity is independent of luminosity for galaxies with the blue colors of spirals. This suggests that, at fixed star-formation history, spiral-galaxy mass is a very weak function of environment. Overdensity does depend on luminosity for galaxies with the red colors of early types; both low-luminosity and high-luminosity red galaxies are found to be in highly overdense regions.Comment: submitted to ApJ

The Sloan Digital Sky Survey: The Cosmic Spectrum and Star-Formation History

We present a determination of the `Cosmic Optical Spectrum' of the Universe, i.e. the ensemble emission from galaxies, as determined from the red-selected Sloan Digital Sky Survey main galaxy sample and compare with previous results of the blue-selected 2dF Galaxy Redshift Survey. Broadly we find good agreement in both the spectrum and the derived star-formation histories. If we use a power-law star-formation history model where star-formation rate $\propto (1+z)^\beta$ out to z=1, then we find that $\beta$ of 2 to 3 is still the most likely model and there is no evidence for current surveys missing large amounts of star formation at high redshift. In particular `Fossil Cosmology' of the local universe gives measures of star-formation history which are consistent with direct observations at high redshift. Using the photometry of SDSS we are able to derive the cosmic spectrum in absolute units (i.e.$W \AA$^{-1}$Mpc$^{-3}$) at 2--5\AA resolution and find good agreement with published broad-band luminosity densities. For a Salpeter IMF the best fit stellar mass/light ratio is 3.7--7.5$\Msun/\Lsun$in the r-band (corresponding to$\omstars h = 0.0025$--0.0055) and from both the stellar emission history and the H$\alpha$luminosity density independently we find a cosmological star-formation rate of 0.03--0.04 h$\Msun$yr$^{-1}$Mpc$^{-3}$ today.Comment: 17 pages, 11 figures, ApJ in press (April 10th 2003

A List of Galaxies for Gravitational Wave Searches

We present a list of galaxies within 100 Mpc, which we call the Gravitational Wave Galaxy Catalogue (GWGC), that is currently being used in follow-up searches of electromagnetic counterparts from gravitational wave searches. Due to the time constraints of rapid follow-up, a locally available catalogue of reduced, homogenized data is required. To achieve this we used four existing catalogues: an updated version of the Tully Nearby Galaxy Catalog, the Catalog of Neighboring Galaxies, the V8k catalogue and HyperLEDA. The GWGC contains information on sky position, distance, blue magnitude, major and minor diameters, position angle, and galaxy type for 53,255 galaxies. Errors on these quantities are either taken directly from the literature or estimated based on our understanding of the uncertainties associated with the measurement method. By using the PGC numbering system developed for HyperLEDA, the catalogue has a reduced level of degeneracies compared to catalogues with a similar purpose and is easily updated. We also include 150 Milky Way globular clusters. Finally, we compare the GWGC to previously used catalogues, and find the GWGC to be more complete within 100 Mpc due to our use of more up-to-date input catalogues and the fact that we have not made a blue luminosity cut.Comment: Accepted for publication in Classical and Quantum Gravity, 13 pages, 7 figure

The Optical, Infrared and Radio Properties of Extragalactic Sources Observed by SDSS, 2MASS and FIRST Surveys

We positionally match sources observed by the Sloan Digital Sky Survey (SDSS), the Two Micron All Sky Survey (2MASS), and the Faint Images of the Radio Sky at Twenty-cm (FIRST) survey. Practically all 2MASS sources are matched to an SDSS source within 2 arcsec; ~11% of them are optically resolved galaxies and the rest are dominated by stars. About 1/3 of FIRST sources are matched to an SDSS source within 2 arcsec; ~80% of these are galaxies and the rest are dominated by quasars. Based on these results, we project that by the completion of these surveys the matched samples will include about 10^7 stars and 10^6 galaxies observed by both SDSS and 2MASS, and about 250,000 galaxies and 50,000 quasars observed by both SDSS and FIRST. Here we present a preliminary analysis of the optical, infrared and radio properties for the extragalactic sources from the matched samples. In particular, we find that the fraction of quasars with stellar colors missed by the SDSS spectroscopic survey is probably not larger than ~10%, and that the optical colors of radio-loud quasars are ~0.05 mag. redder (with 4-sigma significance) than the colors of radio-quiet quasars.Comment: 10 pages, 6 color figures, presented at IAU Colloquium 184. AGN Survey

The velocity function in the local environment from LCDM and LWDM constrained simulations

Using constrained simulations of the local Universe for generic cold dark matter and for 1keV warm dark matter, we investigate the difference in the abundance of dark matter halos in the local environment. We find that the mass function within 20 Mpc/h of the Local Group is ~2 times larger than the universal mass function in the 10^9-10^13 M_odot/h mass range. Imposing the field of view of the on-going HI blind survey ALFALFA in our simulations, we predict that the velocity function in the Virgo-direction region exceeds the universal velocity function by a factor of 3. Furthermore, employing a scheme to translate the halo velocity function into a galaxy velocity function, we compare the simulation results with a sample of galaxies from the early catalog release of ALFALFA. We find that our simulations are able to reproduce the velocity function in the 80-300 km/s velocity range, having a value ~10 times larger than the universal velocity function in the Virgo-direction region. In the low velocity regime, 35-80 km/s, the warm dark matter simulation reproduces the observed flattening of the velocity function. On the contrary, the simulation with cold dark matter predicts a steep rise in the velocity function towards lower velocities; for V_max=35 km/s, it forecasts ~10 times more sources than the ones observed. If confirmed by the complete ALFALFA survey, our results indicate a potential problem for the cold dark matter paradigm or for the conventional assumptions about energetic feedback in dwarf galaxies.Comment: 24 pages, 14 figures, 1 table, accepted for publication in Ap