”Observations and models of high redshift Radio Galaxies and Quasars from the 3rd Cambridge catalog”
Dissertation
zur Erlangung des Grades ”Doktor der Naturwissenschaften” an der Fakultat¨ fur¨ Physik und Astronomie der Ruhr-Universitat¨ Bochum
von Frank Heymann
aus Leipzig
Bochum 2010
b.w. 2
1. Gutachter Prof. Dr. Rolf Chini 2. Gutachter Priv.-Doz. Dr. Dominik Bomans Datum der Disputation 05.07.2010 3
Abstract
This thesis provides new observations of the most powerful high redshift Active Galactic Nuclei (AGN), namely the complete z > 1 3CR sample, and new dust radiative transfer modeling of their measured spectral energy distribution in the infrared. This work is separated into three main parts, two observational sections and one section containing the modeling.
The first part shows observational results in the near (NIR) and mid (MIR) infrared obtained with the Spitzer Space Telescope, to extend the knowledge on high redshift sources. The main aspect of these observations is to study orien- tation dependence of the NIR and MIR emission and to confirm the unification scheme for the most powerful high redshift AGNs.
The second part reports on a pilot study to detect galaxy clustering around high redshift radio sources using the Spitzer data. Because the radio AGN reside in massive host galaxies, they are expected to serve as signposts for cosmic mass peaks. These galaxy clusters are among the most distant known structures and therefore of particular cosmological interest.
The third part explains a newly developed method to solve the radiative trans- fer equation in three dimensional configurations. This method makes use of the parallelization capabilities of modern vector computing units, like the graphics cards. The speed improvement is about a factor of 100. This enables us to model the close environment of AGN in so far unprecedented detail within reasonable computing time. 4 Contents
1 Introduction 7 1.1 Aimofthisthesis...... 7 1.2 ActiveGalacticNuclei ...... 7 1.2.1 Seyfert Galaxies ...... 9 1.2.2 Quasars & Radio Galaxies ...... 10 1.2.3 Unification...... 10 1.3 GalaxyClustering...... 11 1.4 Dustradiation...... 12 1.4.1 Theemissivityofdust ...... 13 1.4.2 Thetemperatureofbiggrains ...... 13 1.4.3 Verysmallgrains ...... 14 1.5 BasicRadiativetransport ...... 15 1.5.1 Definitions...... 15 1.5.2 The general transfer equation ...... 16 1.5.3 Analytical solutions ...... 17
2 Near- and mid infrared photometry 19 2.1 Motivation...... 20 2.2 ObservationsandData ...... 21 2.3 ResultsandDiscussion ...... 24 2.3.1 Radio galaxies as obscured quasars ...... 24 2.3.2 Evolutionary effects and non-thermal contributions . . . . 27 2.4 Conclusions ...... 28
3 The cluster search 35 3.1 Motivation...... 35
5 6 CONTENTS
3.2 Clusteringaround3C270.1 ...... 37 3.2.1 Observationsanddata ...... 37 3.2.2 Results...... 38 3.3 No evidence for clustering around 3C 437 ...... 51 3.3.1 Observationsanddata ...... 51 3.3.2 Results...... 51 3.3.3 Preliminary Conclusion ...... 53
4 Parallel 3D radiative transport 55 4.1 Theory...... 55 4.1.1 MonteCarlomethod ...... 55 4.1.2 Parallelization ...... 58 4.1.3 Pseudo random number generator in parallel ...... 60 TM 4.1.4 CUDA parallelization on graphic cards ...... 60 4.2 Imaging ...... 61 4.2.1 Solarsystem...... 65 4.3 Benchmarktest ...... 65 4.3.1 Sphericalsymmetry(1D)...... 65 4.3.2 Diskgeometry(2D)...... 67 4.3.3 Dustproperties ...... 69 4.3.4 Spiral expansion of disk structure (3D) ...... 69 4.3.5 Clumpy Torus geometry (3D) ...... 70 4.4 Modelingaveragespectra&SEDs...... 75
5 Summary and Outlook 81
Bibliography 87
List of Figures 93
List of Acronyms 101 Chapter 1 Introduction
1.1 Aim of this thesis
At the begin of my thesis, new unprecedented infrared observations of the com- plete high-redshift 3CR sample have been obtained with the Spitzer Space Tele- scope. Therefore, the aim of this thesis is
• to explore the near- and mid-infrared spectral energy distributions of this sample, comprising the most powerful radio-loud AGN,
• to test, how far it is possible with these data to detect galaxy clustering around these mass peaks of the early universe,
• and to develop a new proper 3D Monte Carlo radiative transfer code to model the spectral energy distributions.
These three tasks represent new challenges and thus lead to new results.
1.2 Active Galactic Nuclei
Active Galactic Nuclei (AGN) belong to the the most luminous objects in the universe. The luminosity of an AGN is provided by accretion of matter onto the central supermassive black hole:
˙ ˙ 2 12 ǫ M LAGN = ǫ M c ≈ 1.2 10 L⊙ , (1.1) 0.1 M⊙/yr where ǫ is the efficiency of the mass to radiative energy transfer, L⊙ solar lu- minosity, M⊙ solar mass and M˙ the mass accretion. This leads to a theoretical
7 8 1.2. ACTIVE GALACTIC NUCLEI
upper limit for the central luminosity, named after Arthur Stanley Eddington:
4πGMBH mpc 11 MBH Lmax where mp is the mass of the proton, σT the Thomson cross section for interaction between electrons and protons. The wavelength spectral energy distribution of an AGN exhibits three characteristic features (Elvis et al. 1994) as shown in Fig 1.2: Figure 1.1: Sketch of an AGN continuum spectrum of the nuclear region, without stellar contribution. Three different bumps can be seen (Big Blue Bump in the middle, Infrared Bump on the left and X-ray ’Bump’ on the right). Figure from Manners (2002) • Infrared Bump This feature consists of several components. Dominant is the emission from the hot and warm dust torus (red dashed line) and cooler dust from the host galaxy. The starburst activity in the host galaxy contributes to the far infrared (purple dotted line) followed by a steep decrease of the Infrared bump to the submillimeter (Chini et al. 1989). The local minimum at around 1 m is given by the sublimation temperature of the dust around ∼ 1500 K. CHAPTER 1. INTRODUCTION 9 • Big Blue Bump At shorter waverlengths the minimum turns into the Big Blue Bump (blue dashed-tripple-dotted line). This bump comes from the thermal emission of hot gas (5 000K − 100 000 K) heated by viscous processes, in the accretion disk. The gap in the bump results from absorption of neutral hydrogen and therefore missing data. This optical/UV radiation is efficiently transfered into infrared emission and therefore powers the Infrared Bump. (Miley et al. 1985) • X-ray ’Bump’ The final feature in the AGN continuum is the high energy X-ray ’Bump’. The radiation in this bump is produced by the hot corona above the accre- tion disk (green dashed line) and reflection of the disk (blue dashed-dotted line). The observational data of my thesis provide new constraints on the Infrared Bump (chapter 2), and the model part makes use of all three bumps (chapter 4). 1.2.1 Seyfert Galaxies These galaxies, discovered by Seyfert (1943), contain a bright nucleus with strong emission lines from highly ionised gas (hydrogen, helium, nitrogen, oxygen). The Seyfert galaxies can be divided into two subclasses depending on the existence of broad and narrow emission lines (type 1) or only narrow lines (type 2) (Khachikian & Weedman 1974). The broad lines have velocities of 1 500−10 000 km s−1, may vary on short timescales and can be explained by Doppler broadening. These high velocities can be explained by gas clouds, orbiting the black hole at small distances. It is also possible that these lines are emitted from the accretion disc itself. However due to the extremly high resolution which is neccessary to resolve the innermost part, it is difficult to observe the exact geometry of these objects. The narrow emission lines may by emitted by gas clouds further out. This is strengthened by the fact that the narrow lines are detected in all types of Seyfert galaxies, which implies that the emitting region is large. Breakthrough spectropolarimetric observations revealed, that some Seyfert 2 galaxies contain a hidden broad line region, leading to the AGN unification scheme (section 1.2.3). 10 1.2. ACTIVE GALACTIC NUCLEI 1.2.2 Quasars & Radio Galaxies These two object classes are the powerful radio-loud cousins of the Seyfert galax- ies. On early optical images quasars appeared starlike, which gave these objects the name quasi stellar radio source. In quasars and radio galaxies large struc- tures, the radio lobes, are prominent. Depending on the morphology and power of the radio lobes, the radio sources are subdivided into the FR I and FR II classes (Fanaroff & Riley 1974). We here consider only the edge-brightened pow- erful FR II sources. The radio emission is powered by synchrotron radiation of outflowing material reaching 95% of the speed of light. Assuming an accretion disc and perpendicular to this disc an outflowing jet producing the synchrotron radiation, it is possible to explain these two objects with the orientation of the jet to the line of sight. The quasars, where the jet points to the observer, repre- sent the type 1 AGN and the Radio Galaxies, where the jet is perpendicular to the line of sight, can be classified as type 2 (Barthel 1989; Orr & Browne 1982). The isotropic radiation at meter wavelengths enables to select objects with the same radio power and to study orientation dependent effects. While for the local universe (z < 1) orientation dependencies are fairly well studied, the knowledge is poor at high redshift (z > 1). 1.2.3 Unification Figure 1.2 shows the actual picture of the unification model. It is believed that Type 1 and Type 2 galaxies are in essence the same, and they differ only by the angle at which they are observed. This unification scheme was proposed by Antonucci (1993). In Type 2 Seyferts the broad line component, which is produced by high velocity clouds in the vicinity of the black hole, is obscured by dust in a torus like structure. The obscuration then depends on the observing angle onto the dusty torus. The narrow emission lines are produced in clouds with lower rotation velocities and therefore larger distances to the black hole presumably in a bicone far outside the obscuring torus scale hight. Perpendicular to the accretion disk exists a jet of outflowing material, to remove the angular moment from the accretion disk. This outflow reaches speeds of up to 95% of the speed of light, and produces synchrotron radiation in the radio wavelength. In some Type 2 Seyfert galaxies the broad component can be observed in polarized light. In these objects the broad-line region is scattered by material in the bicone (dust or electrons), which allows us to see it indirectly. This effect was first discovered by Antonucci (1983) in the Type 2 Seyfert NGC 1068. CHAPTER 1. INTRODUCTION 11 Narrow Line Region Broad Line Region Jet Black Accretion Hole Disk Obscuring Torus Figure 1.2: This figure shows the unification scheme from Urry & Padovani (1995). In the center the black hole with its accreating disk is shown. Per- pendicular to this disk is a jet of outflowing material. Close to the black hole are clouds which produce the broad emission line features with speeds of up to 10 000 kms−1. The torus of dust shields the broad line region depending on the observing orientation. Further out is the narrow line region. 1.3 Galaxy Clustering Clusters of nebulae (Zwicky 1937) have been recognised later as cluster of galaxies. When observed visually, clusters appear to be collections of galaxies held together by mutual gravitational attraction. However, their velocities are too large for 12 1.4. DUST RADIATION them to remain gravitationally bound by their attractions, implying the presence of either an additional invisible mass component, or an additional force besides gravity. X-ray studies have revealed the presence of large amounts of intergalactic gas known as the intracluster medium (ICM). This gas is very hot, between 107K and 108K, and therefore emits X-rays in the form of bremsstrahlung and atomic line emission. The total mass of the gas is greater than that of the galaxies by roughly a factor of two. However, this is still not enough mass to keep the galaxies in the cluster. Since this gas is in approximate hydrostatic equilibrium with the overall cluster gravitational field, the total mass distribution can be determined. It turns out that the total mass deduced from this measurement is approximately six times larger than the mass of the galaxies or the hot gas. The missing component is called dark matter and its nature is unknown. In a typical cluster perhaps only 5% of the total mass is in the form of galaxies, maybe 10% in the form of hot X-ray emitting gas and the remainder is dark matter. Noteworthy, Brownstein & Moffat (2006) use a theory of modified gravity to explain X-ray cluster masses without dark matter. Clusters typically have the following properties: • They contain 50 to 1,000 galaxies, hot X-ray emitting gas and large amounts of dark matter 14 15 • They have total masses of 10 to 10 M⊙ and sizes of 2 to 10 Mpc. • The velocity distribution for the individual galaxies is about 800−1000 kms−1. Well known galaxy clusters in the relatively nearby universe include the Virgo cluster, Hercules Cluster, and the Coma Cluster. A very large aggregation of galaxies known as the Great Attractor, dominated by the Norma cluster, is mas- sive enough to affect the local expansion of the universe (Hubble flow). However, at large distances (redshift z > 1) cluster knowledge is very sparse, because of observational difficulties to discern cluster member galaxies from foreground galaxies. In chapter 3 a new pilot study to find high redshift (z ∼ 1.5) clusters is presented. 1.4 Dust radiation This section describes the physical background of the dust emission. A single grain emits according to Kirchhoff’s law, which can be derived from the assump- tion of harmonic oscillators in thermal equilibrium. In section 1.4.1 we discuss CHAPTER 1. INTRODUCTION 13 the thermal emission. The big grain radiation is discussed in section 1.4.2 and the very small grains in section 1.4.3. 1.4.1 The emissivity of dust A dust particle in a radiation field has an equilibrium temperature derived from the condition that it absorbs as much energy as it emits. The ratio of the emission abs coefficient ǫν to the absorption coefficient Kν depends only on the frequency of the radiation field and temperature of the dust grain. Planck discovered that this dependence can be described by a function Bν(T ), named the Planck function. Astrophysical grains are usually not in local thermal equilibrium (LTE). But whenever it is heated by a photon the energy is distributed among the energy states. The distribution of the energy states derived by the Boltzmann function depends only on the dust temperature T . Therefore the grain emission only depends on T and is given, as in LTE, by abs ǫν = Kν Bν(T ) . (1.3) Then the emission into all directions equals 4πǫν. The value ǫν is called emissivity and can reference to a single particle, a unit volume or a unit mass. The dimension for the emissivity per unit volume is expressed in erg s−1 cm−3 Hz−1 ster−1. 1.4.2 The temperature of big grains The dust temperature can be calculated with Equ. 1.3 and the dust temperature T from the balance between heating and cooling by Kabs J dν = Kabs B (T) dν (1.4) ν ν ν ν where the left side of this equation describes the heating and the right side the cooling. Jν is the average of the radiation intensity over all directions. Equ. 1.4 neglects other heating or cooling mechanisms like collisions with gas particles. To account for these additional heating (cooling) one has to add them to Equ. abs 1.4. It is obvious that not only the absolute value of Kν determines the dust abs temperature, which means the grain is at high temperature if Kν is high at the absorption wavelength and low at the emission wavelength. For the simple case of abs a black-body the absorption coefficient is Kν constant. This simplifies equation 14 1.4. DUST RADIATION 1.4 to σ J dν = T 4 , (1.5) ν π with σ the radiation constant. The temperature of a dust black-body at distance r from a star with bolometric luminosity L is given by L = 4σT 4 . (1.6) 4πr2 For the astrophysical example of our solar system, where L = L⊙, one finds: r − 1 T = 279 ( ) 2 K , (1.7) AU which is close to the surface temperature of the earth. This holds even if our planet emits like a perfect black body and the atmosphere reflects nearly 30% of the radiation from the sun. 1.4.3 Very small grains When a very small dust particle absorbs an energetic photon its temperature fluc- tuates. To compute these time variable temperatures their optical and thermal properties are needed. The optical behavior depends on the particle shape and two dielectric functions ǫ1(ω) and ǫ2(ω) and the thermal behavior is determined from the specific heat (see Kr¨ugel 2008). For a large number of very small iden- tical interstellar dust grains the chance for one of them to have a temperature in the interval T...T + dT is P (T )dT. Thus the normalised probability density is: ∞ P (T )dT = 1 . (1.8) 0 For big interstellar grains the temperature fluctuations are small, so that the probability density P (T ) approaches the δ-function δ(T ). The balance between emission and absorption is therefore given in equation 1.4. For very small particles we assume that the radiation fulfills Kirchhoff’s law at any time. In the case of spherical particles we can compute the average monochromatic emission per solid angle by ǫ = π a2 Qabs B (T )P (T ) dT (1.9) ν ν ν The emission of one of these grains is not constant over time, but the whole ensemble radiates at any frequency at a steady rate. The solution of 1.9 is reduced to find P (T ). A more detailed view of the emission of very small particles is CHAPTER 1. INTRODUCTION 15 Figure 1.3: Definition of the radiative intensity Iν described by Guhathakurta & Draine (1989); Draine & Li (2001); Kr¨ugel (2008). 1.5 Basic Radiative transport 1.5.1 Definitions The intensity Iν describes the radiation field and is defined as shown in Fig. 1.3 A surface element dσ at location r receives the power dW = Iν( r, e) cos (θ)dσdΩdν , (1.10) where e is the direction, n is the unit vector of the surface element, θ is the angle between e and n, dΩ is the solid angle and dν is the width of the frequency interval. The dimension of the intensity is erg s−1 cm−2 Hz−1 ster−1. For extended astrophysical sources of solid angle ΩS the observed total flux received at the observer is given by F obs = I (θ, φ)dΩ , (1.11) ν ν ΩS where dΩ= sin(θ)dθdφ (1.12) is the element of solid angle. For constant Iν(θ, φ) over the source Equ. 1.11 simplifies to obs Fν = IνΩS . (1.13) 16 1.5. BASIC RADIATIVE TRANSPORT The mean intensity J is the average over the total solid angle, 1 J = I( e)dΩ . (1.14) 4π 4π The net flux F through a unit surface is: 2π π F = I(θ, φ) cos (θ) sin (φ)dθdφ . (1.15) 0 0 1.5.2 The general transfer equation The light changes its intensity by traveling through a dusty medium. The inten- sity is reduced by absorption and increased by the thermal emission of the dust and the scattering of photons out of and into the line of sight. The intensity change at position r into direction e can be written as dI ν = −KextI ( r, e)+ ǫ , (1.16) ds ν ν ν ext abs sca where Kν = Kν + Kν is the extinction coefficient per unit volume includ- ing the absorption and scattering. The emission and scattering of the dust at temperature T can be written as sca abs Kν ′ ′ ′ ǫν = Kν Bν(T )+ pν(e , e)Iν(e )dΩ . (1.17) 4π 4π ′ ′ The factor pν(e , e) is the phase function for scattering from e into e. A funda- mental quantity in the light propagation is the optical depth, τν. When the light travels from a point P1 to P2 the optical depth is given by P2 τ = Kextds , (1.18) ν ν P1 or ext dτν = Kν ds . (1.19) For dusty media the extinction coefficient is independent of temperature, which simplifies the integral in Equ. 1.18. With the column density Nd, the number of grains in a cylinder with a base of 1 cm2, the optical depth is given by ext τν = NdCν , (1.20) CHAPTER 1. INTRODUCTION 17 ext where Cν is the extinction coefficient for one grain. 1.5.3 Analytical solutions It is possible to solve Equ. 1.16 for pure extinction without emission. This can be the light from an object, which is only decreased by the intervening dust. The solution is then −τ I(τν)= I(0)e , (1.21) where I(0) is the intensity of the source. The intensity I(τ) decreases exponential with the optical depth. The solution of Equ. 1.16 for pure emission without extinction (Kext = 0) is s I(s)= I(0) + ǫ(s′) ds′ , (1.22) 0 where s is the path length and ǫ the emission coefficient (Equ. 1.17). This situation is a good approximation for the sub-millimeter radiation of interstellar dust. The intensity for pure emission without a background source scales with the path length s. The more general case of background source emission plus dust emission and absorption is τ −τ ′ τ−τ ′ ′ Iν(τ)= I(0)e + Sν(τ )e dτ , (1.23) 0 where S(τ) is the source function. For the simple case of a constant source function it is possible to solve the integral in equ. 1.23. We assume a dust cloud of uniform temperature T , the source function is given by the Planck function Sν(τ)= Bν(T ) . (1.24) This gives the solution for the intensity and flux of the absorption plus emission of dust on a background source by τ −τ Iν(τ)= I(0) e + Bν(T ) (1 − e ) (1.25) and τ −τ Fν(τ)= F (0) e + Bν(T ) (1 − e )Ω . (1.26) If the source emission is known or neglectable this equation gives the possibility to get the dust temperature and optical depth and therefore the dust mass. The 18 1.5. BASIC RADIATIVE TRANSPORT limiting solutions for Equ. 1.25 are B(T )τ + I(0) : for τ << 1 Iν(τ)= (1.27) B(T ) : for τ >> 1 . For low optical depth and no background source the radiation intensity received at the observer is directly proportional to the optical depth. The intensity of the received dust emission can be used to investigate dust properties like the dust mass or the grain frequency dependence. For very high optical depth the intensity is given by the black-body radiation. This means all information about the dust, with the exception of the dust temperature, is lost. Chapter 2 Near- and mid infrared photometry This part represents one of the developed papers during my PhD thesis. My contribution in this large collaboration of many coauthors was mainly: • Photometry of the entire high redshift 3C sample (Tables 5.1 and 5.2 in appendix 3C sources). • Developing a pipeline to extract the photometry from the Spitzer mosaic images. • Interactive inspection, verification and improvement of the spectral energy distribution of these faint sources. Abstract: The 178 MHz-selected sample is unbiased with respect to orientation and there- fore suited to study orientation-dependent effects in the most powerful active galactic nuclei (AGN). Quasar and radio galaxy subsamples matched in isotropic radio luminosity are compared. The quasars all have similar spectral energy dis- tributions (SEDs), nearly constant in ν Fν through the rest 1.6–10 m range, con- sistent with a centrally heated dust distribution which outshines the host galaxy contribution. The radio galaxy SEDs show larger dispersion, but the mean radio galaxy SED declines from rest 1.6 to 3 m and then rises from 3 to 8 m. The radio galaxies are on average a factor 3–10 less luminous in this spectral range than the quasars. These characteristics are consistent with composite emission from a heavily reddened AGN plus starlight from the host galaxy. The mid- infrared colors and radio to mid-infrared spectral slopes of individual galaxies are also consistent with this picture. Individual galaxies show different amounts of extinction and host galaxy starlight, consistent with the orientation-dependent unified scheme. 19 20 2.1. MOTIVATION 2.1 Motivation When exploring the general evolution of galaxy populations across cosmic times, a particular challenge is to distinguish between black hole and star-forming ac- tivity. Star formation and obscuring dust go hand in hand, and black-hole-driven active galactic nuclei (AGN) are also surrounded by dust mainly distributed in a disk/torus-like geometry (Antonucci 1993). There is evidence that AGN mainly power the near- and mid-IR emission (NIR, ∼2 m; MIR, ∼10 m) from hot nu- clear dust, while starbursts contribute mainly to the far-infrared (FIR, ∼60 m) luminosity (Rowan-Robinson 1995; Vernet et al. 2001; Schweitzer et al. 2006). Using the MIR/FIR luminosity ratio as an indicator for the relative AGN and starburst contributions, numerous studies have found an increase of AGN/star- burst activity with total luminosity and redshift, but the validity of this trend is still under discussion because of selection effects on the various samples. More se- riously, an unfavorable AGN orientation could cause MIR obscuration (e.g. Pier & Krolik 1993), leading to a fundamental observational degeneracy: a low MIR/- FIR luminosity ratio can be due to either a high star-forming contribution or to an AGN in which the hot dust is obscured. The spectral energy distribution (SED) of an obscured AGN may thus mimic that of a starburst-powered source. While this degeneracy has now been widely examined at low/intermediate lumi- nosity and redshift (z < 1), it has still to be explored for the most luminous sources at high redshift (z > 1). In order to assess galaxy and AGN evolution in the universe, we therefore need to understand this AGN/starburst degeneracy for a population of luminous high-redshift sources. A crucial step towards this is to study the orientation dependence of the NIR and MIR emission of high-redshift AGN. Orientation-dependent effects can only be tested and quantified with AGN samples having type 1 (unobscured) and type 2 (obscured) subsamples matched in isotropic emission. The clean AGN tracers — optical, [O III] λ5007 A,˚ NIR, and X-ray (<10 keV) — all fail to fulfill this requirement. The [O II] λ3727A˚ emission, while isotropic (Hes et al. 1993), is probably dominated by extended starbursts and shocks (Best 2000) rather than by the AGN. Therefore, the only feasible way is low-frequency (meter-wavelength) radio selection because the integrated emis- sion from the radio lobes is optically thin and essentially isotropic. This makes radio-loud AGN particularly attractive for studying orientation-dependent prop- erties at other wavelengths and, after sorting out the influence of radio jets/lobes on the emission, for generalizing conclusions about orientation-dependent effects CHAPTER 2. NEAR- AND MID INFRARED PHOTOMETRY 21 to the much larger population of radio-quiet AGN. The brightest low-frequency-selected AGN sample is the 3CR compilation Spinrad et al. (1985). The powerful double-lobed radio galaxies (henceforth sim- ply called radio galaxies) are supposed to be misaligned quasars Barthel (1989). Based on IRAS co-added scans and a few individual detections, Heckman et al. (1992, 1994) already noted an average MIR/FIR difference between 3CR quasars and radio galaxies. More comprehensive MIR and FIR spectrophotometry from ISO is in hand (as compiled by Siebenmorgen et al. (2004) and by Haas et al. (2004)) as well as from Spitzer (e.g., Shi et al. (2005); Haas et al. (2005); Ogle et al. (2006); Cleary et al. (2007), providing a basis to study the z < 1 3CR ob- jects. These sources are, however, a factor of five less radio-luminous on average than the most powerful radio sources seen at higher redshift, and the lower indi- cated accretion power may reflect different source physics. The higher-luminosity population can be sampled by the 3CR sources at 1 H0 = 71 km/s/Mpc, Ωm = 0.27 and ΩΛ = 0.73. 2.2 Observations and Data With the Spitzer Space Telescope (Werner et al. 2004), we have obtained the entire sample of 64 high-z 3CR sources using the instruments IRAC (3.6–8.0 m, Fazio et al. 2004), the IRS blue peak-up array (16 m, Houck et al. 2004) and MIPS (24 m, Rieke et al. 2004). Most observations are performed in our guaranteed time program (PID 40072; PI G. Fazio) with on-source exposure times 4×30 s (each IRAC band), 4×14 s (IRS), and 10×10 s (MIPS). A few sources have been observed in other programs, and we use the published photometry if available (e.g., PID 3329; PI D. Stern, Seymour et al. 2007). For IRAC, we used the basic calibrated data products (BCD, version S16) and co-added them to 0′′.869 pixels using the latest version of IRACProc (Schuster et al. 2006). This optimally handles the slightly under-sampled IRAC PSF in order to assure accurate point-source photometry. For IRS, we used the post- 22 2.2. OBSERVATIONS AND DATA 1047 3CR at z > 1 298 191 9 257 46 220.2 119 m, erg/s ) 469.1 10 68.2 µ 305.1 208.1 45 ( rest 8 10 173 ν L quasars galaxies ν Si absorption 1044 1045 ν Lν ( rest 178 MHz, erg/s ) Figure 2.1: Infrared versus radio luminosity of the 3CR sample at z > 1 prior to normalization. ’x’ symbols denote quasars; circles and squares denote radio galaxies. Superposed crosses indicate radio galaxies that show evidence of sili- cate absorption (§2.3.1). The vertical long-dashed lines mark the range of our luminosity-matched quasar and radio galaxy subsamples. The dotted lines indi- cate L8 m/L178 MHz ratios of 1, 10, and 100. The radio galaxies are grouped into several SED classes in Fig. 2.3 and §2.3.1. The color-coding and symbols are: green circle (A), red circle (B), red square (C), blue square (D), blue circle (E). The two low-excitation radio galaxies 3C 68.1 and 3C 469.1 are labeled with their 3C numbers, as are sources outside the luminosity range of our analysis. BCD pipeline product, version S16. For MIPS, we used custom routines to modify the version S17 BCD files to remove instrumental artifacts (e.g., residual images) before shifting and co-adding to create the final mosaics. All sources are well seen on the images in all filters. The sources were extracted and matched using the SExtractor tool Bertin & Arnouts (1996). We used sufficiently large apertures so that aperture corrections are small (<5%). The photometric errors are typically smaller than 10% but increase for faint sources; exceptions are 3C 225A and 3C 294, where nearby bright stars make the photometry uncertain in the shorter IRAC bands. As of 2008 April, 24 quasars and 38 radio galaxies have been observed, covering CHAPTER 2. NEAR- AND MID INFRARED PHOTOMETRY 23 the complete high-z 3CR sample with the exception of the quasar 3C 245 and the radio galaxy 3C 325. All 62 sources have IRAC measurements and are observed in at least one of the 16 and 24 m bands (54 sources at 16 m and 60 sources at 24 m). For the analysis it is desirable to compare rest frame SEDs with the same wavelength sampling. Depending on the redshift of our sources (1