626

Lecture 2: number counts and functions

How much of the extragalactic background light can we identify, and how much is unidentified (unresolved)?

The next step is to count all the galaxies we can find and see how much light they contain

So we now want to form the galaxy number counts at all wavelengths B, R, z/

15/ x 15/

B (1.7hrs) R (5.2hrs) z’ (3.9hrs)

B (27) R (26.4) z’ (25.4)

Optical Image

Capak et al. 2003 Measuring Galaxy

Galaxies, unlike , are not point sources

The Hubble Space Telescope can resolve (i.e. detect the extended nature of) essentially all galaxies

Even from the ground, most galaxies can easily be distinguished from stars morphologically and on the basis of their colors

Measuring Galaxy Luminosities

Define the surface brightness of a galaxy I as the amount of light from the galaxy per square arcsecond on the sky

Consider a small square patch, of side D, in a galaxy at distance d:

d D

Angle patch subtends on sky α = D/d Finding Galaxies in an Image

Finding galaxies in an image and constructing the number counts is the subject of the first student project

Essentially we look for all objects above a limiting surface brightness and covering more than a specified area

SExtractor is a commonly used package 20 cm Radio Image from the VLA Number Counts

We then need to measure the fluxes of all the objects that we find

The number counts are simply the number of objects we find in a given flux or bin per unit area of the sky Measuring Galaxy Luminosities

Consider again the small square patch, of side D, in a galaxy at distance d

d D

Angle patch subtends on sky α = D/d If the luminosity of all the stars within the patch is L, then the total flux is L F = 4"d 2 Define surface brightness as:

L /4"d 2 L ! I = = D2 /d 2 4"D2

Units of I are mag arcsec-2; i.e., if a galaxy has a surface -2 brightness! of 20 mag arcsec , then we receive as many photons from one square arcsecond of the galaxy’s image as from a of 20th magnitude

-2 Centers of galaxies have IB ~ 18 mag arcsec To measure the total amount of light coming from a galaxy, need to integrate the surface brightness across the galaxy image. This leads to a related problem: galaxies do not have sharp well-defined edges

Typically integrate out to some limiting isophote; e.g., sum up all the light coming from regions with surface -2 brightness IB < 25 mag arcsec . Alternatively, we can measure to a surface brightness level relative to the central surface brightness (say 1 percent)

This measure of the total galaxy brightness is called an isophotal magnitude---numerous variations are possible An alternative and commonly used procedure is to measure the magnitudes in an aperture that contains most of the galaxy light and to apply an average correction for the light outside the aperture

These are usually called corrected aperture magnitudes The UBL – Unidentified (Unresolved) Background Light Key Issues: • We need to know the EBL to high accuracy (removal of spurious foregrounds, calibration issues) • We need to determine the extragalactic galaxy counts to very faint levels and properly measure their total fluxes (calibration and getting the fainter surface brightness parts of the galaxies) -2 -1 -1 -1 N(m) deg 0.5mag iν cgs Hz sr

γ=0.4

Counts from U through K extend well into γ<0.4 regime

Madau & Pozzetti (2000) MNRAS 312, L9 The UBL (above & beyond known counts)

# UBL= EBL - I" = $ i" (m)dm -#

where #0.4(m+const) i" (m) =10 N(m) ! Multiplying flux of the sources by N(m), the! d ifferential counts per unit magnitude per unit area

Units 1 Jansky = 10-26 W m-2 Hz-1 1 nW m-2 sr-1 = 3000/λ(µm) MJy sr-1

Contributions to the Optical EBL

Capak et al. 2004 Pozzetti et al. 2000 Ground-based data Ground + HST data Contributions from integrated counts

nW m-2 sr-1

• Extrapolated counts enable contribution to integrated light from known populations to be evaluated as function of wavelength

-2 -1 • Find EBL = ∫ ν Iν dν = 55 nW m sr ; largest contribution in NIR 4 • Observational issues: isophotal losses and (1+z) redshift dimming CCDDF-FN-N 503 sources 12 MMss

RED: 0.5-2 keV

GREEN: 2-4 keV

BLUE: 4-8 keV

3D70. sAoulercxesander 460 sq arcmin CDF-N + HAWAII

COMBINED

CDF-S COMBINED CDF-N CDF-S

SSA13 Moretti et al. 2003 Measurements of the 2-8 keV background How much unidentified light do we have?

Resolved background Hickox & Markevitch 2006 Moretti et al. 2003 Measured backgrounds

Contribution from counts The Galaxy Luminosity Function

The luminosities of galaxies span a very wide range---most luminous ellipticals are 107 more luminous than the faintest dwarfs

Luminosity function, Φ(L), describes the relative number of galaxies of different luminosities

Definition: If we count galaxies in a representative volume of the Universe, Φ(L)dL is the number of galaxies with luminosities between L and L + dL

Identical to the definition of the stellar luminosity function

Luminosity functions are easiest to measure in clusters of galaxies, where all the galaxies have the same distance Coma Cluster The Schecter Luminosity Function A convenient approximation to the luminosity function was suggested by Paul Schecter in 1976:

* $ L ' $ L ' dL "(L)dL = n#& ) exp& + ) % L# ( % L # ( L#

In this expression:

• n* is a normalization factor which defines the overall ! density of galaxies (number per cubic Mpc) • L* is a characteristic galaxy luminosity. An L* galaxy is a bright galaxy, roughly comparable in luminosity

to the . A galaxy with L < 0.1 L* is a dwarf. • α defines the `faint-end slope’ of the luminosity function. α is typically negative, implying large numbers of galaxies with low luminosities. The Schecter Luminosity Function

The Schecter function is: • a fitting formula that does not distinguish between galaxy types • as with the stellar mass function, parameters must be determined observationally: Illustrative -3 -3 n* = 8 x 10 Mpc - related to mean galaxy numbers density 10 L* = 1.4 x 10 Lsun - luminosity of galaxies that dominate light output α = −0.7 - lots of faint galaxies

33 -1 …where Lsun = 3.9 x 10 erg s is the The Schecter function plotted for different faint-end slopes, α = -0.5 (red), α = -0.75 (green), α = -1 (blue): Properties of the Schecter Luminosity Function The total number of galaxies per unit volume with luminosity greater than L is:

# n = (L)dL This integral is not expressible in $ " terms of elementary functions L If α < -1, n diverges as L tends toward zero. Obviously unphysical, so Schecter law must fail for very low luminosity! galaxies (or have a larger alpha).

The luminosity density (units Solar luminosities per cubic Megaparsec) is given by: # $ "(L)LdL L

Dominated by galaxies with L ~ L* for typical value of α.

! Significance of the Schecter Function

Does the Schecter form of the luminosity function have any deeper significance, or is it just a good fitting function?

• Luminous galaxies become exponentially rarer at high enough luminosities. Similar to a general result in cosmology (confusingly, Press-Schecter theory) - the number density of massive objects (e.g. clusters) drops off exponentially at high masses.

• But remember that Schecter function is a composite which includes all types of galaxies. Luminosity function of any individual class of galaxies looks very different - irregular galaxies typically have much lower luminosities than ellipticals for example. Mass Function of Galaxies

For stars, measurements of the luminosity function can be used to derive the Initial Mass Function (IMF).

For galaxies, this is more difficult:

• Mass to light ratio (M/L) of the stellar population depends upon the star formation history of the galaxy.

• Image of the galaxy tells us nothing about the amount and distribution of the dark matter.

More difficult measurements are needed to try and get at the mass function of galaxies. Local Inventory of Stars

Relevant papers: Fukugita et al. 1998, ApJ, 503, 518 Fukugita & Peebles 2004, ApJ, 616, 643 Survey data: Cole et al. 2001, MNRAS, 326, 255 (2dF+2MASS) Kauffmann et al. 2003, MNRAS, 341, 33 (SDSS)

Stellar density: measure local infrared LF, Φ (LK), and scale by a mean mass/light ratio (M/LK), which depends on the initial mass function (IMF) Most of the mass is in lower-mass stars, and most of the light is from solar mass stars. One empirically tries to determine the M/LK ratio, but the only information below 1 solar mass is from the disk of the Milky Way One calculates all of the galaxy properties one would expect for a particular IMF to see if it gives a self-consistent description Some Initial Mass Functions where the contributions are coming from

Salpeter > 1 M reproduces colors and Hα properties of spirals

IMF < 1 M makes minor contribution to light but is very important for mass inventory (Salpeter diverges at low- Mass (M) mass end) Stellar Mass Function at z=0 17,000 K-band selected galaxies with K<13.0 2dFGRS redshifts & 2MASS photometry

Stellar K mass Stellar Mass Density at Present-Day

Cole et al find Ωstars h = 0.0014 ± 0.00013; M/LK = 0.73 (Miller-Scalo)

Ωstars h = 0.0027 ± 0.00027; M/LK = 1.32 (Salpeter)

Fukugita & Peebles: Ωstars h = 0.0027 ± 0.0005 (5% in brown dwarfs)

Ωgas = 0.00078 ± 0.00016 (H I, He I, H2) NB: a comparison with WMAP result for the baryon density, only 6% of baryons are in stars! Baryon density is dominated by gas (WHIM), not stars Would expect the integral of the past activity to agree with the locally-determined stellar density (Fukugita & Peebles 2004) If both the local value and the cumulative star formation are calibrated with the same IMF, then at least the description is consistent However, the IMF may not be invariant; one could be averaging over a lot of IMFs over time, so the IMF uncertainties are significant Student Project 2

The evolution of the galaxy luminosity functions with redshift will be the topic of the second project Summary

• We have made astonishing progress in understanding the extragalactic populations across the spectrum in the last 10- 20 years with EBL measurements and number counts across the spectrum • The EBL observations remain challenging and controversial because of instrumental effects and dominant foregrounds and determining how much light remains unidentified (the UBL) is still a subtle problem • However, while there may be small (10%-ish) residual excesses at many wavelengths there is not a large amount of unidentified light • Local galaxy luminosity functions are also now well established as a consequence of wide field surveys such as 2MASS, SDSS and the 2DF and we can use these to estimate the local baryon mass function of galaxies