Draft version April 27, 2018 Preprint typeset using LATEXstyleemulateapjv.12/16/11
⇤CDM COSMOLOGY FOR ASTRONOMERS J. J. Condon1 National Radio Astronomy Observatory, 520 Edgemont Road, Charlottesville, VA 22903, USA and A. M. Matthews Department of Astronomy, University of Virginia, Charlottesville, VA 22904, USA Draft version April 27, 2018
ABSTRACT The homogeneous, isotropic, and flat ⇤CDM universe favored by observations of the cosmic microwave background can be described using only Euclidean geometry, locally correct Newtonian mechanics, and the basic postulates of special and general relativity. We present simple derivations of the most useful equations connecting astronomical observables (redshift, flux density, angular diameter, brightness, local space density, . . . ) with the corresponding intrinsic properties of distant sources (lookback time, distance, spectral luminosity, linear size, specific intensity, source counts, . . . ). We also present an analytic equation for lookback time that is accurate within 0.1% for all redshifts z. The exact equation for comoving distance is an elliptic integral that must be evaluated numerically, but we found a simple approximation with errors < 0.2% for all redshifts up to z 50. ⇡ Keywords: cosmology: distance scale — cosmology: theory — galaxies: distances and redshifts — galaxies: evolution
1. INTRODUCTION “Innocent, light-minded men, who think that astron- omy can be learnt by looking at the stars without knowl- edge of mathematics will, in the next life, be birds.”— Plato, Timaeus According to the cosmological principle and confirmed by observation (e.g., Figure 1), the universe is isotropic and homogeneous on large scales. It began with a very dense “big bang” and has been expanding uniformly ever since. General relativity can describe both the ge- ometry and expansion dynamics of the universe. How- ever, general relativity permits spatially curved universes which are mathematically complicated. Extragalactic as- tronomers were once faced with the choice of learning general relativity or applying published relativistic re- sults without really understanding them, at the risk of being birds in the next life. In today’s ⇤CDM (⇤for dark energy with constant energy density and CDM for cold dark matter) con- cordance model, the universe is spatially “flat,” so its geometry is Euclidian, its expansion is not a↵ected by curvature, and locally correct Newtonian calculations arXiv:1804.10047v1 [astro-ph.CO] 26 Apr 2018 can be extended to cosmological scales. Fortunately for non-mathematicians, flatness allows simple (no tensors) Figure 1. Positions of the N 4 104 radio sources stronger derivations of accurate equations for the kinematics and than S =2.5mJyat1.4GHzareindicatedbypointsonthis⇠ ⇥ dynamics of cosmic expansion that an undergraduate equal-area plot covering the sky within 15 of the north celestial physics major can understand. Such derivations are pre- pole. Nearly all of these sources are extragalactic and so distant (median redshift z 1) that their distribution is quite isotropic. sented in Sections 2 and 3, and the main results used by h i⇠ observational astronomers are developed in Section 4. See David Hogg’s useful “cheat sheet” (Hogg 1999) vature, and the books by Peebles (1993) and Wein- listing results from relativistic models with nonzero cur- berg (1972) for their derivations. The astropy.cosmology Python package at http://docs.astropy.org/ en/ [email protected] [email protected] stable/cosmology/ contains utilities for calculating 1 The National Radio Astronomy Observatory is a facility many of the quantities discussed in this paper. of the National Science Foundation operated under cooperative We use the subscript 0 to distinguish present (red- agreement by Associated Universities, Inc. shift zero) values of evolving quantities. Unless oth- 2 Condon & Matthews erwise noted, all numerical results are based on a Equation 1 defines the recession velocity of the nearby 1 1 Hubble constant H0 = 70 km s Mpc so h observer: 1 1 ⌘ H /(100 km s Mpc )=0.7, plus the following nor- 0 ˙ a˙ malized densities at redshift zero: total ⌦ = 1, bary- vr D =˙aD0 = D = HD . (3) 0 ⌘ a onic and cold dark matter ⌦0,m =0.3, radiation and ✓ ◆ 2 5 relic neutrinos ⌦ = h 4.2 10 , and dark energy Successive wave crests of light emitted with frequency 0,r · ⇥ ⌦0,⇤ =1 (⌦0,m +⌦0,r) 0.7. Note that many authors ⌫e and wavelength e = c/⌫e are separated in time by ⇡ 1 just write ⌦,⌦m,⌦r, and ⌦⇤ without the subscript 0 dt = ⌫e in the source frame. If the source is receding to indicate the present values of these densities. Also, from the observer with velocity vr c,successivewaves 19 7 ⌧ 1Mpc 3.0857 10 km and 1 yr 3.1557 10 s. must travel an extra distance vr dt = vr/⌫e, so their ob- ⇡ ⇥ ⇡ ⇥ served wavelength is 2. EXPANSION KINEMATICS c v v According to the equivalence principle, all fundamen- = + r = + r (v c) (4) o ⌫ ⌫ e e c r ⌧ tal observers locally at rest relative to their surround- e e ✓ ◆ ings anywhere in the isotropically expanding (and hence and vr is measurable via the first-order Doppler shift: homogeneous) universe are in inertial frames, and their clocks all agree on the time t elapsed since the big bang. vr o e = (vr c) . (5) This universal time t is sometimes called world time or c e ⌧ cosmic time, and it equals the proper time of all funda- mental observers. Fundamental observers didn’t have to The redshift z of a source is defined by be present at the creation to synchronize their clocks; the o e temperature of the cosmic microwave background (CMB) z , (6) ⌘ radiation is a suitable proxy for time. Likewise, the e CMB appears isotropic only to fundamental observers, and the domain of this definition extends to all z. Note and others can use the CMB dipole anisotropy (Kogut that most recession “velocities” reported by astronomers 2 2 et al. 1993) to deduce their usually small (v c ,where are actually vr = cz and may be much larger than the 1 ⌧ c 299792 km s is the vacuum speed of light) peculiar vacuum speed of light. velocities⇡ and correct for them if necessary. Combining Equations 3 and 5 for a nearby source at Homogeneity and isotropy are preserved if and only if distance D = c t gives the small proper distance D(t) between any close pair of vr D a˙ a˙ a observers expands as = = = t = (vr c) . (7) c c a a a ⌧ D(t)=a(t) D0 , (1) Integrating the local Equation 7 over time: 1 where D0 is the proper distance now and a(t)istheuni- o d da versal (meaning, it is the same at every position in the = (8) a universe) dimensionless scale factor that grew with time Z e Za from a 0 just after the big bang to a(t ) 1today ⇡ 0 ⌘ gives the global result that o/ e =1/a and ⌫o/⌫e = a. (Figure 2). Equation 1 applies to any expansion that Thus the observable redshift z of any distant source can preserves homogeneity and isotropy—the separations of be used to calculate the scale factor a of the universe dots on a photo being enlarged, for example. The cos- when it emitted the light seen today: mological expansion a↵ects only to the separations of 1 non-interacting objects, and the dots representing rigid a =(1+z) . (9) rulers, gravitationally bound galaxies, etc. do not ex- pand with the universe. Because 1 (1 + z) < is the reciprocal of the scale factor 0
16 3 4 /c 7.566 10 Jm is the radiation con- 4 Gyr ago. In the distant future completely dominated ⇡ ⇥ stant. The T0 2.73 K CMB has energy density by dark energy, the Hubble parameter will asymptoti- 14 ⇡ 3 1/2 U 4.20 10 Jm and gravitational mass density 1 1 0 cally approach H = H0 ⌦0,⇤ 84 h km s Mpc ⇡ ⇥2 31 3 ⇢ = U /c 4.67 10 kg m . Massless relic neu- 1 1 ⇡ ⇡ 0,r 0 59 km s Mpc and the scale factor will grow ex- ⇡ ⇥ 14 3 1 trinos have energy density U0 2.86 10 Jm and ponentially [a exp(Ht)] with a time scale H ⇡ ⇥ 31 3 10 1 / ⇡ gravitational mass density ⇢0,⌫ 3.18 10 kg m 1.17 10 h yr 17 Gyr. (Peebles 1993). The total from photons⇡ ⇥ plus neutrinos is ⇥ ⇡ 5 2 5 ⌦0,r 4.2 10 h 8.6 10 . 4. RESULTS Mass⇡ conservation⇥ of⇡ non-relativistic⇥ matter implies 3 3 4.1. Cosmic Times ⇢m a =(1+z) . Inthe⇤CDM model, dark en- ergy/ is assumed to behave like a cosmological constant: The lookback time tL(z) to a source at redshift z can be 0 0 calculated by inserting Equations 25 and 26 into Equa- ⇢⇤ a =(1+z) . The density of radiation (and mass- / 4 4 tion 18: less neutrinos) scales as ⇢r a =(1+z) because / 3 3 z the number density of photons is a =(1+z) tL dz0 2 2 / = . (27) and the mass E/c = h⌫/c of each photon scales as tH0 0 (1 + z0)E(z0) 1 1 1 Z E (1 + z) a .Insertingtheseresultsinto Equations/ / 21, 22, and/ 23 leads to This integral cannot be expressed in terms of elementary functions, but the integrand is smooth enough that the 2 ⇢ H ⌦0,m ⌦0,⇤ ⌦0,r integral can be evaluated numerically by Simpson’s rule, = = + + . (24) 2 3 0 4 at least for finite z. ⇢0 H0 a a a The proper age of the universe t(z) at redshift z is Using Equation 9 to replace the scale factor a by the 1 1 better evaluated in terms of a =(1+z) : observable (1 + z) yields the dynamical equation spec- ifying the evolution of H: 1 t (1+z) ada = 4 1/2 . (28) H 3 4 1/2 tH0 0 (⌦0,ma +⌦0,⇤a +⌦0,r) =[⌦0,m(1 + z) +⌦0,⇤ +⌦0,r(1 + z) ] . (25) Z H0 The ratios tL/tH0 and t/tH0 are plotted as functions of The symbol redshift in Figure 5. The present age of the ⇤CDM universe starting with the big bang (z = )ist0 3 4 1/2 9 1 1 ⇡ E(z) [⌦ (1 + z) +⌦ +⌦ (1 + z) ] (26) 0.964 tH 9.42 10 h yr 13.5 Gyr for ⌦0,m =0.3. ⌘ 0,m 0,⇤ 0,r 0 ⇡ ⇥ ⇡ is a convenient shorthand for subsequent calculations. Figure 4 shows H/H0 = E(z) as a function of z for ⌦0,m =0.3.
Figure 5. The normalized lookback time (tL/tH0 )andthenor- malized age (t/tH0 )inaflat⇤CDMuniversewith⌦0,m =0.3.
Figure 4. The normalized Hubble parameter H/H0 is nearly 1/3 Figure 5 shows that redshifts z 1 contribute little to constant after (1 + z) . (⌦0,⇤/⌦0,m) 1.33 and ⇢⇤ >⇢m. the age of the universe, so extremely good analytic ap- 3/2 ⇡ H/H0 (1 + z) at higher redshifts when ⇢m >⇢⇤,and / 2 proximations to tL/tH0 , t/tH0 , and t0/tH0 can be made by H/H0 (1 + z) at the highest redshifts z>zeq 3500 when / ⇠ 4 ⇢r >⇢m. ignoring the radiation term ⌦0,r(1 + z) that dominates E(z) only during the brief period when z>zeq 3500 6 4 1 ⇠4 The densities of radiation and matter were equal when (t/tH0 3.4 10 , or t 5 10 h yr 7 10 yr). 4 3 ⇠ ⇥ ⇠ ⇥ ⇠ ⇥ ⌦0,r(1 + zeq) =⌦0,m(1 + zeq) at zeq =(⌦0,m/⌦0,r) This simplifies Equation 27 to 1 3500. The density of matter fell below that of z t dz dark⇡ energy at z (⌦ /⌦ )1/3 1 0.33 about L 0 , (29) ⇡ 0,⇤ 0,m ⇡ t ⇡ (1 + z )[⌦ (1 + z )3 +⌦ ]1/2 H0 Z0 0 0,m 0 0,⇤ ⇤CDM Cosmology for Astronomers 5 which can be integrated analytically. See Appendix A for analytic approximations to tL/tH0 (Equation A5) and t/tH0 (Equation A10). The current age of the universe normalized by the Hubble time is very nearly
1/2 t0 2 1+⌦0,⇤ 1/2 ln 1/2 . (30) tH ⇡ (1 ⌦0,⇤) 0 3⌦0,⇤ " #
Figure 6 plots t0/tH0 from Equation 30 as a function of the matter density parameter ⌦0,m =1 ⌦0,⇤.Inthe limit ⌦ =1 ⌦ 0, the universe would expand 0,m 0,⇤ ! exponentially and t0/tH0 would diverge.
Figure 6. The normalized age of the universe t0/tH0 as a function of ⌦0,m =1 ⌦ 0,⇤
The observable redshift z is the traditional proxy for lookback time tL in models of cosmological evolution. However, it can be misleading because the lookback time is an extremely nonlinear function of redshift when z & 1. For example, the top panel in Figure 7 shows the Madau & Dickinson (2014) best fit to the star formation rate density (SFRD) in solar masses per year of per (co- moving) Mpc3 as a linear of function of lookback time back to z = 8. This plot accurately displays the fact that only 10% of today’s stellar mass was assembled be- fore z 2.9. The very nonlinear upper abscissa indi- cates the⇡ redshifts z matching the lookback times on the lower abscissa. The time between z = 0 and z =1is 7.8Gyr,butthetimebetweenz = 2 and z = 3 is only ⇡ 1.1 Gyr. The middle panel plots the same function ⇡ Figure 7. The three panels in this figure show the Madau & as a linear function of redshift, with lookback time now Dickinson (2014) fit to the evolving star formation rate density as on the very nonlinear upper abscissa. This plot makes linear functions of lookback time tL (top panel), redshift z (middle it look like 10% of today’s stellar mass was assem- panel), and scale factor a (bottom panel). Clearly a is a much bled before z 2.9. The nonlinearity at high redshifts better proxy for tL than z is. is primarily caused⇡ by the fact that (1 + z)istherecip- The vacuum speed of light c is invariant, so the light rocal of the scale factor a.When is plotted as a linear travel distance DT corresponding to lookback time tL is function of a (bottom panel), the upper abscissa showing z lookback time is nearly linear and the plot of (a) looks dz0 DT = ctL = c . (31) much more like the plot of (tL). (1 + z )H Z0 0 4.2. Light Travel Distance The light travel distance in meters can be interpreted physically as the number of meter sticks laid end-to-end 6 Condon & Matthews that the photon must pass on its journey from the source Inserting Equations 17 and 26 into Equation 34 gives to the observer. z The light travel distance is of limited use in cosmog- dz0 DC = DH0 . (35) raphy because it is the distance between two events oc- 0 E(z0) curring at two di↵erent proper times, t and t0.This Z limitation can be illustrated by a nonrelativistic terres- Unfortunately, this indefinite integral for DC cannot be trial example: two trains moving in opposite directions expressed in terms of elementary functions, only ellip- with speed v passed each other at time t = 0 (Figure 8). tic integrals that must be evaluated numerically. It is The horn on one train sounded at time te when the dis- smooth enough to be evaluated by Simpson’s rule (Ap- tance between the trains was de =2vte. The sound pendix B). Alternatively, the simple empirical fit reached the other train at a later time to when the dis- DC(fit) DH0 / [a/(1 a)+0.2278+ tance between the trains was do =2vto. The sound ⇡ travel distance is d = vt + vt = c (t t ), where 0.2070 (1 a)/(0.785 + a) T e o s o e cs >vis the speed of sound. Some algebra reveals that 0.0158 (1 a)/(0.312 + a)2] , (36) dT = de [cs/(cs v)] = do [cs/(cs + v)] equals neither de 1 nor do. where a =(1+z) , can be used to avoid the nu- merical integration for most astronomical applications. Equation 36 is accurate to within 0.2% for z . 50 and ⌦0,m =0.3 (Figure 9).
Figure 8. Two trains moving with speed v passed each other at time t =0,andthetrainmovingleftemittedasoundattimete when the trains were separated by de =2vte.Thesoundreached train moving right at time to when the trains were separated by Figure 9. Equation 36 fits the comoving distance up to z 50 d =2vt .Thesoundtraveldistanced = vt +vt equals neither ⇡ o o T e o within 0.2% for ⌦m =0.3. de nor do. The comoving distance in meters equals the number 4.3. Comoving Coordinates of meter sticks laid end-to-end that would connect the source to the observer today. It is the proper distance Comoving coordinates expand with the universe. The between us and the source at time t0,sothecomov- comoving distance D between any close (D D ) C ⌧ H0 ing distance is the most fundamental distance for use in pair of fundamental observers defined by ⇤CDM cosmology. D(t) The “observable” universe refers to the sphere in which DC = D0 (32) signals emitted at any time t>0 and traveling at the ⌘ a(t) vacuum speed of light could have reached the observer (see Equation 1) is independent of t and equals the today. Its light-travel radius is therefore ct0, and the present proper distance D0. Thus comoving rulers are current comoving radius of the observable universe is like rubber bands connecting neighboring fundamental 1 dz0 observers, and their markings agree with rigid rulers to- R0,C = DH0 . (37) day. At redshift z, the vacuum speed of light in comoving 0 E(z0) 1 Z coordinates is a c =(1+z)c and If h =0.7 and ⌦0,m =0.3, then R0,C 3.24 DH0 1 ⇡ ⇡ dDC =(c/a) dt =(1+z)cdt. (33) 3.24 2998 h Mpc 13.9 Gpc. This number is not par- ticularly· significant,⇡ but it is often used to make amus- In the homogeneous universe, summing over the “lo- ing calculations like the following: The comoving vol- cal” comoving distances dD yields the global “line of 3 C ume of the observable universe is V0 =4⇡R / 3 sight” comoving distance to a distant source at any red- 0,C ⇡ 1.12 104 Gpc3 3.3 1080 m3 and the present mean shift z: ⇥ ⇡ ⇥27 3 t density is ⇢0 9.2 10 kg m , so the total mass of 0 the observable⇡ universe⇥ is M = V ⇢ 3 1054 kg. D = c (1 + z) dt0 0 0 0 C Figure 10 shows how the observable⇡ radius⇥ R varies Zt 0,C 0 with the normalized matter density ⌦0,m. dt0 = c (1 + z0) dz0 . (34) The comoving volume VC measured in comoving co- dz Zz ✓ 0 ◆ ordinates is valuable for tracking the cosmic evolution ⇤CDM Cosmology for Astronomers 7
for a rigid source. Equations 43 and 44 imply
D D = C (45) A (1 + z)
in a flat universe. The angular diameter in arcsec (1 arcsec = ⇡/648000 rad) of a source with fixed proper diameter l = 1 kpc is shown as a function of redshift in Figure 11.?
Figure 10. The normalized comoving radius of the observable universe today R0,C/DH0 as a function of ⌦0,m. of source populations because the number of permanent objects (e.g., immortal fundamental observers, baryons, or galaxies if mergers are ignored) in a comoving vol- ume element is constant. (See Appendix C for a sample application of VC to calculate source counts from local luminosity functions.) The Euclidean geometry of a flat universe implies that the total comoving volume out to redshift z is Figure 11. The angular diameter ✓ of a source with a fixed proper 3 diameter l =1kpchasaminimum✓ 0 00. 118 at z 1.6ina 4⇡DC ? ⇡ ⇡ V = , (38) ⇤CDM universe with h =0.7and⌦m =0.3. C 3 and the comoving volume in a shell covering solid angle The proper motion µ of a source is its observed angular ! sr between z and z + dz is speed across the sky: 2 d✓ dVC = !DC dDC . (39) µ . (46) ⌘ dt Di↵erentiating Equation 35 yields For example, the radio source in the quasar 3C 279 at dz dD = D (40) z =0.5362 appears to consist of a stationary “core” C H0 E(z) plus moving components whose proper motions µ 1 ⇠ so 0 00. 0005 yr were measured by very long baseline inter- 2 ferometry (Piner et al. 2003). !DC DH0 dVC = dz . (41) The proper-motion distance DM of a source with E(z) proper transverse velocity v = dl /dt(z)definedby ? ? v 4.4. Angular-diameter and Proper-motion Distances D ? (47) M ⌘ µ The flat ⇤CDM universe is Euclidean, so the angular distance ✓ (rad) between two fundamental observers at is also called the “transverse” comoving distance. In the same redshift with comoving transverse separation l0 terms of the angular diameter distance, at comoving distance D is simply C dl dt l DM = ? = DA(1 + z) , (48) ✓ = 0 (✓ 1) . (42) d✓ dt(z) D ⌧ C where dt(z) is the proper time measured by a fundamen- A rigid source (e.g., a transverse meter stick or a gravi- tal observer at the source redshift z. The proper-motion tationally bound galaxy) has a fixed proper transverse distance is just the angular diameter distance multiplied length l , so its comoving transverse length is l0 = by the (1 + z) time dilation factor. In a flat universe, (1 + z) l ? and ? (1 + z)l DM = DC , (49) ✓ = ? . (43) DC so the “line of sight” comoving distance DC equals the The angular diameter distance DA is a “convenience” “transverse” comoving distance DM and the two can be distance defined to satisfy the static Euclidean equation treated simply as a “the” comoving distance. In the l non-Euclidean geometry of a curved universe DM = DC D ? (44) (Hogg 1999). 6 A ⌘ ✓ 8 Condon & Matthews
4.5. Luminosities and Fluxes Let L be the total (or bolometric) luminosity (emitted power, SI units W) of an isotropic source measured in the source frame and F be the total flux (power received 2 per unit area, SI units W m ) in the observer’s frame. (If the source is not isotropic, L should be replaced by 4⇡ times the power per steradian beamed in the direction of the observer.) In a static Euclidean universe, the inverse-square law accounts for the transverse spatial dilution of photons spread over the area A0 of the spherical surface contain- ing the observer and centered on the source: F = L/A0. In the Euclidean but expanding ⇤CDM universe, the present area of the spherical shell centered on a source 2 at redshift z and containing the observer is A0 =4⇡DC. It is not 4⇡ times the square of the distance DT cov- ered by the photons in an expanding universe; the name “inverse-square (of the distance) law” is misleading and “inverse area law” would be better. In an expanding but Euclidean flat universe, the ob- Figure 12. The luminosity distance DL,thecomovingdistance served flux is lower than in a static universe because (1) DC,thelighttraveldistanceDT,andtheangulardiameterdis- tance DA,allnormalizedbythecurrentHubbledistanceDH0 ,are the observed rate at which photons cross the t = t0 sur- compared for ⌦ =0.3. face centered on the source and containing the observer 0,m is a factor (1 + z) lower than the rate at which they were emitted and (2) the observed energy Eo = hc/ o of each redshifted photon is a factor (1+z) lower than its energy Ee = hc/ e in the source frame. Consequently L 1 2 F = . (50) 4⇡D2 1+z ✓ C ◆✓ ◆ The luminosity distance DL is another “convenience” distance defined by the form of the static Euclidean inverse-square law: L F 2 . (51) ⌘ 4⇡DL Figure 13. The comoving distances DC are plotted for ⌦0,m = Equation 50 implies 0.28 (top curve) through ⌦0,m =0.34 (bottom curve) in steps of ⌦0,m =0.01, all normalized by DC(⌦0,m =0.3) (heavy line). D =(1+z)D . (52) L C in the source frame and doesn’t refer to any particular frequency. The ⌫ in parentheses is the actual frequency 4.6. Comparison of Distance Types in the source frame, so L⌫ (1.4 GHz) is the power per unit Figure 12 compares the luminosity distance DL,the frequency emitted at ⌫ =1.4 GHz in the source frame. comoving distance DC, the light travel distance DT, and The spectral flux density F⌫ (⌫) (or S) of a source is the the angular-diameter distance DA in a ⇤CDM universe observed flux per unit frequency in the observer’s frame 2 1 with current matter density parameter ⌦0,m =0.3. All of (SI units W m Hz , or the astronomically practical the plotted distances are normalized by the current Hub- 26 2 1 1 Jy 10 Wm Hz ). ble distance D = c/H 2998 h Mpc 4280 Mpc. ⌘ H0 0 ⇡ ⇡ All of the photons received in a narrow logarithmic These distances also vary slowly and smoothly with frequency range centered on frequency ⌫ were emitted in ⌦0,m, but the e↵ect of changing ⌦0,m cannot be repre- the equally narrow logarithmic frequency range centered sented by a scale factor like h. Figure 13 shows the ratios on [(1 + z)⌫], so the bolometric Equation 50 implies of DC(⌦0,m)toDC(⌦0,m =0.3) for matter densities from ⌦0,m =0.28 (top curve) through 0.34 (bottom curve) L⌫ [(1 + z)⌫] ⌫F⌫ (⌫)=[(1+z)⌫] (53) that include the best measured ⌦0,m =0.315 0.013 4⇡D2 (Planck Collaboration XLVI 2016) and its quoted± uncer- L tainty. Near ⌦ =0.3, dD /d⌦ 0.006 at z =1 and 0,m C 0,m ⇡ and dDC/d⌦0,m 0.013 when z 1. L [(1 + z)⌫] ⇡ F (⌫)=(1+z) ⌫ . (54) ⌫ 4⇡D2 4.7. Spectral Luminosities and Flux Densities L
The spectral luminosity L⌫ (⌫) of a source is its lumi- The leading factor of (1 + z) in Equation 54 comes from 1 nosity per unit frequency (SI units W Hz ). In this no- bandwidth compression: photons emitted over the fre- tation, the subscript ⌫ just means “per unit frequency” quency range [(1+z) ⌫] are squeezed into the frequency ⇤CDM Cosmology for Astronomers 9 range ⌫ in the observer’s frame. The wavelength counterparts of Equations 53 and 54 are The spectral index ↵ between frequencies ⌫1 and ⌫2 is defined by [ /(1 + z)] L [ /(1 + z)] F ( )= 2 . (61) 4⇡DL ln[L(⌫1)/L(⌫2)] ↵(⌫1,⌫2) + . (55) and ⌘ ln(⌫1/⌫2) 1 L [ /(1 + z)] [Beware that some authors define ↵ with the opposite F ( )=(1+z) . (62) 4⇡D2 sign.] In terms of ↵[⌫, (1 + z)⌫], Equation 54 becomes L
↵+1 L⌫ (⌫) 4.8. Magnitudes and K Corrections F⌫ (⌫)=(1+z) 2 . (56) 4⇡DL The apparent magnitude m and absolute magnitude M of a source are related by If the spectral luminosity distance DL⌫ is defined by analogy with Equation 51: DL L m M = 5 log10 + K , (63) F ⌫ , (57) 10 pc ⌫ 4⇡D2 ✓ ◆ ⌘ L⌫ where then Equation 54 implies DL DM 5 log10 (64) (↵+1)/2 ⌘ 10 pc DL⌫ = DL (1 + z) . (58) ✓ ◆
Figure 14 shows DL (z)/DH for a range of spectral in- is the bolometric distance modulus (Figure 15) in magni- ⌫ 0 0.4 dices. Steep-spectrum (↵ 1) synchrotron sources tudes (1 mag 10 ). The K correction converts the ⇡ apparent magnitude⌘ measured through a filter covering have DL⌫ DL and consequently are quite faint at high redshifts.⇡ Flat-spectrum self-absorbed synchrotron a fixed wavelength range (Oke & Sandage 1968; Hogg sources and optically thin free-free emitters (↵ 0) et al. 2002) in the observer’s frame to yield the absolute are only slightly stronger. For a source with ↵ =⇡ +1, magnitude over the same wavelength range in the source frame. The bolometric K correction is zero. There are DL⌫ = DC. Blackbody emission in the long wavelength Rayleigh-Jeans limit has ↵ +2. Optically thin dusty many magnitude systems covering di↵erent bandpasses ⇡ galaxies have ↵ +3 at wavelengths 0.1 . e . 1mm, and having di↵erent zero points (Blanton & Roweis 2007) ⇠ relating m =0tofluxdensity. so their DL⌫ DA and their submillimeter flux densi- ties are nearly⇠ independent of redshift over a broad range centered on z 1.6. ⇠
Figure 15. The bolometric distance modulus DM as a function of redshift z for h =0.7, ⌦0,m =0.3. Figure 14. Spectral luminosity distances normalized by the cur- rent Hubble distance are plotted for spectral indices ↵ = 1, 0, +1, +2, and +3. ⌦0,m =0.3inallcases. Equation 54 implies that L [(1 + z)⌫] K = 2.5 log (1 + z) ⌫ + The spectral flux density F measured per unit wave- 10 2 3 4⇡DL length (SI units W m ) is related to F by ⇢ ⌫ L (⌫) 2.5 log ⌫ F d = F⌫ d⌫ (59) 10 4⇡D2 | | | | ⇢ L so c c L⌫ [(1 + z)⌫] F = F and L = L . (60) K = 2.5 log10 (1 + z) . (65) 2 ⌫ 2 ⌫ L (⌫) ⇢ ⌫ 10 Condon & Matthews
3 In terms of the spectral index ↵, B⌫0(⌫)=(1+z) B⌫ [(1 + z)⌫] . (74) K = 2.5 log (1 + z)↵+1 10 For a source with spectral index ↵, = 2.5(↵ +⇥ 1) log10(1 +⇤z) . (66) ↵ 3 B⌫0(⌫)=(1+z) B⌫ (⌫) . (75) In terms of wavelengths, The Planck brightness spectrum of a blackbody source at temperature T is 1 L [ /(1 + z)] K = 2.5 log (1 + z) . (67) 10 L ( ) 1 ⇢ 2h⌫3 h⌫ B (⌫ T )= exp 1 . (76) ⌫ | c2 kT 4.9. Spectral Lines " ✓ ◆ # Spectral lines are narrow (line width ⌫ ⌫) emission If the source is at redshift z, the observed spectrum or absorption features in the spectra of gaseous⌧ or ionized 3 sources. The total line luminosity L is related to the total B (⌫)=(1+z) B [(1 + z)⌫ T ] ⌫0 ⌫ | line flux F and the line flux density F⌫ by Equation 51, 3 1 2h (1 + z)⌫ h⌫ (1 + z) so = exp 1 L =4⇡D2 F =4⇡D2 F ⌫. (68) c2 (1 + z) k T L L ⌫ ⇢ 2 (77) In Equation 68, the SI units of F are W m = 26 10 Jy Hz. However, line fluxes are often reported in is just the Planck spectrum B⌫ (⌫ T0) of a blackbody 1 | the dimensionally confusing units Jy km s based on with temperature T0 = T/(1 + z). The source of the the nonrelativistic Doppler equation T0 2.73 K CMB seen today is the T 3000 K black- body⇡ surface of last scattering when the⇠ universe became ⌫ v 1 . (69) transparent at z 1100. ⌫ ⇡ c ⌧ ⇠ Solving Equation 69 for v yields the factor needed to 4.11. Nonrelativistic Approximation Errors 1 convert from Jy km s to Jy Hz: In the low-redshift limit, it is tempting to calculate intrinsic source parameters from the observables using 1 ⌫ 1kms Hz . (70) the nonrelativistic distance approximation ⇡ 299792 cz See Carilli & Walter (2013) for a detailed discussion of DN . (78) this equation and its uses. ⌘ H0 Figure 16 displays ratios of the relativistically correct lu- 4.10. Total Intensity and Specific Intensity minosity distance DL, comoving distance DC, light travel The total intensity or bolometric brightness B of a distance DT, and angular diameter distance DA in a source is the power it emits per unit area per unit solid ⇤CDM universe with ⌦0,m =0.3toDN. 2 1 angle ! (SI units W m sr ). In a static Euclidian universe, intensity is conserved along any ray passing through empty space and the brightness B0 seen by an observer at rest with respect to the source equals B. For a source at redshift z in the Euclidean but expanding ⇤CDM universe, the observed brightness B0 can be cal- culated with the aid of Equations 45 and 52. For a source at any comoving distance DC, expansion multiplies its luminosity distance by (1 + z) and divides its angular diameter distance by (1 + z), so
B0 F0 ! 2 2 = =(1+z) (1 + z) (71) B F !0 and B0 4 =(1+z) . (72) B
The specific intensity or spectral brightness B⌫ (⌫) of a source is its power per unit frequency per unit solid 2 1 1 angle (SI units W m Hz sr ) at frequency ⌫ in the source frame. Figure 16. The ratios of the relativistically correct distances DL, All of the photons received in a narrow logarithmic DC, DT,andDA to the nonrelativistic distance approximation D cz/H0 at low redshifts z indicate the errors that can result N ⌘ frequency range centered on ⌫ were emitted in the equally from using DN. narrow logarithmic frequency range centered on [(1 + z)⌫], so the bolometric brightness Equation 72 implies For example, the bolometric luminosity L of a source 2 [(1 + z)⌫] B [(1 + z)⌫] at redshift z is proportional to DL, so the luminosity ⌫B = ⌫ . (73) calculated using D instead will be too low by a factor of ⌫0 (1 + z)4 N ⇤CDM Cosmology for Astronomers 11
2 (DL/DN) . It will be 5% too low for a source at redshift At 4.85 GHz 4C +39.05 has flux density S 200 mJy, 4 1 ⇡ z 0.032 (cz 10 km s ), and the maximum redshift so its spectral index (Equation 55) is at⇡ which the⇠ luminosity error is < 10% is z 0.065 4 1 ⇡ ln(635 mJy/200 mJy) (cz 2 10 km s ). ↵ =+ 0.9 . (79) The⇠ calculated⇥ linear size of a source is proportional ln(1.4 GHz/4.85 GHz) ⇡ to DA Figure 17. The bent radio source 4C +39.05 (1.4 GHz contours at 12 1 0.4 25.448 1 L⌫ 7.98 10 W Hz 10 · 1, 2, 4,..., 128 mJy beam in a 4500 FWHM Gaussian beam) ⇡ ⇥ · originated± in the elliptical galaxy 2MASX 02005301+3935003 at 23 1 the center of this 6 6 optical finding chart (gray scale). 1.2 10 W Hz (84) 0 ⇥ 0 ⇡ ⇥ 14 at ⌫e 1.390 10 Hz in the source frame. The light-travel time from 4C +39.05 can be found by From⇡ Equation⇥ 75, the observed 1.4 GHz spectral inserting z =0.0718 and ⌦0,⇤ 1 ⌦0,m =1 0.3=0.7 brightness of the ↵ 0.9 radio source at z 0.0718 ⇡ ⇡ ⇡ into Equation A5; it is tL/tH0 0.0682. For h =0.7, is lower than its 1.4 GHz specific intensity in the source ⇡ 8 tH0 14 Gyr (Equation 12) and tL 9.5 10 yr. frame by the factor ⇡ ⇡ ⇥ For h =0.7 and ⌦0,m =0.3, its comoving distance from either Equation 35 or Equation 36 is DC 303 Mpc. B⌫0 ↵ 3 3.9 ⇡ =(1+z) 1.0718 0.76 . (85) The comoving volume within the sphere (! =4⇡ sr) of B⌫ ⇡ ⇡ this radius is (Equation 38) V =4⇡D3 /3 4 3.14 C C ⇡ · · (303Mpc)3/3 1.17 108 Mpc3. 4.13. Example Calculation: A Source Population The angular-size⇡ distance⇥ (Equation 45) to 4C +39.05 The 1.4 GHz spectral luminosity of a star-forming is DA = DC/(1 + z) 283 Mpc, so its projected linear galaxy is a linear and dust-unbiased tracer of the recent ⇡ size is l = ✓DA 200 arcsec ⇡/648000 rad/arcsec star formation rate (SFR) (Murphy et al. 2011): 283 Mpc? 0.27 Mpc.⇡ Notice that· the angular diameter· of this source⇡ would remain ✓>30 if it were moved to 00 SFR 21 L1.4GHz any redshift (Figure 11). 1 =1.0 0.1 10 1 . (86) M yr ± ⇥ W Hz ! ! 12 Condon & Matthews If the comoving space density of 1.4 GHz sources in star- forming galaxies is ⇢(L1.4GHz), then the local luminosity- weighted spectral power density function (Equation C8) is U (L z) = ln(10)L2 ⇢(L z) . (87) dex ⌫ | ⌫ ⌫ | Figure 19. The brightness-weighted count S2n(S)ofall1.4GHz radio sources is shown by the heavy curve. The lower light curve is the count of radio sources that would be produced by a non- evolving population of star-forming galaxies, and the upper light curve is the count that would result from pure luminosity evolution consistent with the Madau & Dickinson (2014) formula for the evolving SFRD (z). Udex(L1.4GHz z)=[ (z)/ (0)] Udex(L1.4GHz z 0). Inserting this |U (L z) into Equation C10| ⇡ yields dex 1.4GHz| Figure 18. The local 1.4 GHz spectral energy density function the higher thin curve in Figure 19. Both thin curves agree Udex of galaxies whose radio emission is powered primarily by star for the low-redshift sources at high flux densities, but the formation, not by AGNs. evolving SFRD produces a peak at log[S(Jy)] 5 that is much closer to the observed faint-source count.⇠ The observed Udex(L1.4GHz z 0) (Condon & Matthews 2018) is shown by the| data⇡ points in Figure 18, and it can be approximated by the function (dotted curve We thank Ken Kellermann, Eric Murphy, Kristina Ny- in Figure 18) land, and Mark Whittle for their valuable comments and suggestions. 1 U (L z 0) L (W Hz ) dex 1.4GHz| ⇡ ⇡ 1.4GHz · 3 3 1 L1.4GHz 4.0 10 Mpc dex ⇥ L · ✓ ⌫⇤ ◆ 1 2 L1.4GHz exp 2 log 1+ , (88) " 2 L⌫⇤ !# 21 1 where L1⇤.4GHz =1.7 10 W Hz , = 0.24, and =0.585. ⇥ The evolution of the SFRD (z) can be constrained by comparing the observed brightness-weighted 1.4 GHz source count S2n(S) (Condon et al. 2012) shown as the heavy curve in Figure 19 with counts predicted by Equa- tion C10 for various evolving Udex(L1.4GHz z). The total count has two peaks, the peak at log[S(Jy)]| 1 pro- duced by AGN-powered radio sources and the⇠ peak at log[S(Jy)] 4.5 attributed to star-forming galaxies. For the case⇠ of no evolution in the SFRD ,inserting Udex(L1.4GHz z)=Udex(L1.4GHz z 0) and the me- dian spectral index| ↵ 0.7 yields| the⇡ lower light curve in Figure 19. The slope⇡ of this curve is 0.5 at high flux densities because the stronger star-forming⇡ galaxies are at such low redshifts that their counts approach the static Euclidean limit S5/2n(S) = constant. The Madau & Dickinson (2014) model for the evolution of the SFRD is shown in Figure 7. Their result might be interpreted in terms of pure luminosity evolution: the comoving density of star-forming galaxies is constant and the luminosity of each galaxy is proportional to ,so ⇤CDM Cosmology for Astronomers 13 APPENDIX ANALYTIC APPROXIMATIONS FOR LOOKBACK TIME AND AGE Equation 18 for the lookback time tL: D z dz t = L = t 0 (A1) L c H0 (1 + z )E(z ) Z0 0 0 can be integrated analytically if the currently small radiation term ⌦0,r is ignored in E(z0). Then t z dz L 0 . (A2) t ⇡ (1 + z )[⌦ (1 + z )3 +⌦ ]1/2 H0 Z0 0 0,m 0 0,⇤ 3/2 Substituting ⌦ =1 ⌦ and x =(1+z0) reduces Equation A2 to 0,⇤ 0,m t 2 1 dx L . (A3) 1/2 3/2 2 tH0 ⇡ 3⌦ (1+z) (1 ⌦0,⇤)/⌦0,⇤ + x 0,⇤ Z The indefinite integral p dx =ln x + a2 + x2 + C. (A4) pa2 + x2 Z p can be found in integral tables or evaluated via the trigonometric substitution x = a tan(u). Thus 1/2 t 2 1+⌦ L ln 0,⇤ (A5) 1/2 3/2 3 tH0 ⇡ 3⌦ "(1 + z) + (1 + z) +(1 ⌦0,⇤)/⌦0,⇤ # 0,⇤ In the limit z , t becomes the present age of the universep t : !1 L 0 1/2 t0 2 1+⌦0,⇤ 1/2 ln 1/2 . (A6) tH ⇡ (1 ⌦0,⇤) 0 3⌦0,⇤ " # 3 The fractional errors in Equations A5 and A6 are < 10 for all z and all ⌦0,m > 0.1 only because the high redshifts at 4 which the omitted ⌦0,r(1 + z0) term is significant contribute little to tL and t0. Equation A5 can be used to calculate accurate lookback time di↵erences tL = tL(z + z) tL(z) only if z 3500. The fractional errors in di↵erential 3 2 ⌧1 lookback times are < 10 if z<10 and rise to 10 at z 60 and to 10 at z 500. The age of the universe at redshift z was ⇠ ⇠ t 1 dz 0 (A7) t ⇡ (1 + z )[⌦ (1 + z )3 +⌦ ]1/2 H0 Zz 0 0,m 0 0,⇤ 4 8 without the radiation term ⌦0,r(1 + z0) . This approximation is safe at redshifts z . 25 (t & 0.045 tH0 6 10 yr) 6 4 1 4 ⇠ ⇥ because radiation dominated E(z) for only the first t 3.4 10 tH 5 10 h yr 7 10 yr after the big bang (Section 4.1). ⇠ ⇥ ⇠ ⇥ ⇠ ⇥ 1 t (1+z) a1/2 da 3 1/2 . (A8) tH ⇡ [(1 ⌦ )+⌦ a ] 0 Z0 0,⇤ 0,⇤ Substituting x a3/2⌦1/2 yields the elementary form ⌘ 0,⇤ 3/2 1/2 (1+z) ⌦ t 2 ⇤ dx 1/2 2 1/2 . (A9) tH ⇡ [(1 ⌦0,⇤)+x ] 0 3⌦0,⇤ Z0 This is similar to Equation A4, so 1/2 1/2 t 2 ⌦0,⇤ 3/2 ⌦0,⇤ 1/2 ln (1 + z) + 3 +1 . (A10) tH ⇡ 1 ⌦0,⇤ (1 + z) (1 ⌦0,⇤) 0 3⌦⇤ (✓ ◆ ) 2 Equation A10 directly gives the same t0/tH0 as Equation A6, t/tH0 with fractional errors < 10 for all z<25, and smaller errors in time di↵erences t = t at high z than Equation A5. L 14 Condon & Matthews NUMERICAL CALCULATION OF COMOVING DISTANCE Equation 35 for comoving distance: z dz D = D 0 , (B1) C H0 E(z ) Z0 0 where E(z) [⌦ (1 + z)3 +⌦ +⌦ (1 + z)4]1/2 , (B2) ⌘ 0,m 0,⇤ 0,r cannot be expressed in terms of elementary functions. However, DC varies smoothly with both redshift z (Figure 12) and normalized matter density ⌦0,m (Figure 13), so the integral can be approximated by Simpson’s rule. For very 1 large redshifts z 1, it is more e cient to integrate over a =(1+z) instead. Integrating Equation 33: dDC =(c/a) dt (B3) gives t0 1 1 1 1 dt 1 dt 1 a0 da0 da0 DC = c = c da0 = c da0 = c = DH (B4) a a da a 2 a˙ a 2H 0 a 2E(a ) Zt Za 0 0 Za 0 0 Zz 0 Za 0 0 and finally 1 da0 DC = DH . (B5) 0 (⌦ a +⌦ a 4 +⌦ )1/2 Za 0,m 0 0,⇤ 0 0,r The FORTRAN function dcmpc below evaluates Equation B5 to return an accurate DC (in Mpc) for redshifts z even into the photon-dominated era z>zeq 3500, given h 0.7 and ⌦0,m 0.3. It can be copied and pasted directly into a text editor such as Emacs.⇡ The corresponding⇠ Python function⇠ dc.py is available at https://github.com/allison-matthews/astro-cosmo C ********************************************************************** function dcmpc (z, h, Omega0m) C ********************************************************************** real*8 aofz, da, simp aofz = z aofz = 1. / (1. + aofz) Omega0r = 4.2e-05 / h**2 Omega0l = 1. - Omega0m - Omega0r dh0mpc = 2997.92458 / h C Evaluate Condon & Matthews 2018, PASP, Eq. B5 using Simpson’s rule: C simp = (da / 3.) * (y0 + 4.*y1 + 2.*y2 + 4.*y3 + 2.*y4 + ... + yn) C where n must be even n = 10 if (z > 0.1) n = 100 if (z > 3.0) n = 1000 if (z > 90.) n = 10000 if (z > 2700.) n = 100000 da = n da = (1. - aofz) / da simp = 0.0 do 100 isub = 0, n a = isub a = a * da + aofz y = 1. / sqrt (Omega0m * a + Omega0l * a**4 + Omega0r) C test for even or odd y subscript and assign weighting factor itest = (isub / 2) * 2 if (itest == isub) factor = 2. if (itest /= isub) factor = 4. if (isub == 0) factor = 1. if (isub == n) factor = 1. simp = simp + factor * y 100 continue dcmpc = dh0mpc * (da / 3.) * simp return end ⇤CDM Cosmology for Astronomers 15 LUMINOSITY FUNCTIONS, SOURCE COUNTS, AND SKY BRIGHTNESS Let N(>S) be the number of sources per steradian stronger than flux density S F⌫ , n(S) dN/dS be the di↵erential source count, and ⌘(S, z) dS dz be the number of sources per steradian with⌘ flux densities⌘ S to S + dS in the redshift range z to z + dz. The spectral luminosity function ⇢(L z) dL is the comoving number density of sources ⌫ | at a given redshift z having spectral luminosities L⌫ to L⌫ + dL⌫ , and dVC is the comoving volume element covering ! = 1 sr of sky between z and z + dz. The number of sources equals the comoving density times the comoving volume: ⌘(S, z) dS dz = ⇢(L z) dL dV . (C1) ⌫ | C For sources with spectral indices ↵, 2 1 ↵ 2 1 ↵ L⌫ =4⇡DL (1 + z) F⌫ =4⇡DC (1 + z) S (C2) (Equation 56) and D2 D dV = C H0 dz (C3) C E(z) (Equation 41). Thus 2 2 1 ↵ DCDH0 ⌘(S, z) dS dz = ⇢(L z)4⇡D (1 + z) dS dz . (C4) ⌫ | C E(z) Multiplying both sides by (1 + z)↵ 1 L 2 S2 = ⌫ (C5) 4⇡D2 C gives (1 + z)↵ 1D S2⌘(S, z)=L2 ⇢(L z) H0 . (C6) ⌫ ⌫ | 4⇡E(z) Frequently the spectral luminosity function is specified as the density of sources per decade of luminosity ⇢ (L z) = ln(10)L ⇢(L z) . (C7) dex ⌫ | ⌫ ⌫ | The luminosity-weighted spectral luminosity function U (L z) L ⇢ (L z) = ln(10)L2 ⇢(L z) (C8) dex ⌫ | ⌘ ⌫ dex ⌫ | ⌫ ⌫ | 1 3 1 3 1 (SI units W Hz m dex =Jm dex , the same as energy density) emphasizes the luminosity ranges con- tributing the most to the spectral luminosity density. In terms of these quantities, (1 + z)↵ 1D S2⌘(S, z)=U (L z) H0 (C9) dex ⌫ | 4⇡ ln(10)E(z) and the total brightness-weighted source count is obtained by integrating over all redshifts: ↵ 1 D 1 (1 + z) S2n(S)= H0 U (L z) dz . (C10) 4⇡ ln(10) dex ⌫ | E(z) Z0 2 2 1 1 Note that S n(S) has dimensions of spectral brightness (SI units W m Hz sr or astronomically practical units 1 Jy sr ). REFERENCES Kochanek, K., Pahre, M. A., Falco, E. E., et al. 2001, ApJ, 560, 566 Birkho↵, G. D. 1923, Relativity and Modern Physics (Cambridge, Madau, P., & Dickinson, M. 2014, ARA&A, 52, 415 MA: Harvard University Press) Murphy, E. J. et al. 2011, ApJ, 737:67 Blanton, M. R., & Roweis, S. 2007, AJ, 133, 734 Oke, J. B., & Sandage, A. 1968, ApJ, 154, 210 Carilli, C. L., & Walter, F. 2013, ARA&A, 51:105 Peebles, P. J. E. 1993, Principles of Physical Cosmology, Condon, J. J., Cotton, W. D., Fomalont, E. B., Kellermann, K. I., (Princeton, NJ: Princeton University Press) Miller, N., Perley, R. A., Scott, D., Vernstrom, T., & Wall, Piner, B. G., et al. 2003, ApJ, 588, 716 J. V. 2012, ApJ, 758:23 Planck Collaboration XLVI 2016, A&A, 596, A107 Condon, J. J., & Matthews, A. M. 2018, in prep. Riess, A., et al. 2016, ApJ, 826:56 Freedman, W. L. 2017, arXiv:astro-ph/1706.02739v1 Weinberg, S. 1972, Gravitation and Cosmology: Principles and Hogg, D. W. 1999, arXiv:astro-ph/9905116v4 Applications of the General Theory of Relativity,(NewYork: Hogg, D. W., Baldry, I. K., Blanton, M. R., & Eisenstein, D. J. John Wiley & Sons) 2002, arXiv:astro-ph/0210394 Kogut, A., et al. 1993, ApJ, 419:1