The evolution of star formation in quasar host galaxies

We have used far-infrared data from IRAS , Infrared Space Observatory ( ISO ), Spitzer Wide-Area Infrared Extragalactic (SWIRE), Submillimetre Common User Bolometer Array (SCUBA) and Max-Planck Millimetre Bolometer (MAMBO) to constrain statistically the mean far-infrared luminosities of quasars. Our quasar compilation at redshifts  0 < z < 6.5  and I -band luminosities  −20 < I _AB < −32  is the first to distinguish evolution from quasar luminosity dependence in such a study. We carefully cross-calibrate IRAS against Spitzer and ISO , finding evidence that IRAS 100-μm fluxes at <1 Jy are overestimated by ∼30 per cent. We find evidence for a correlation between star formation in quasar hosts and the quasar optical luminosities, varying as star formation rate (SFR)  ∝ L ^0.44±0.07_opt  at any fixed redshift below   z = 2  . We also find evidence for evolution of the mean SFR in quasar host galaxies, scaling as  (1 + z )^1.6±0.3  at   z < 2  for any fixed quasar I -band absolute magnitude fainter than −28. We find no evidence for any correlation between SFR and black hole mass at  0.5 < z < 4  . Our data are consistent with feedback from black hole accretion regulating stellar mass assembly at all redshifts.     

ASCA spectroscopy of the luminous infrared galaxy NGC 6240: X-ray emission from a starburst and a buried active nucleus

We present an X-ray spectroscopic study of the prototype far-infrared galaxy NGC 6240 from ASCA . The soft X-ray spectrum (below 2 keV) shows clear signatures of thermal emission well described by a multitemperature optically thin plasma, which probably originates in a powerful starburst. Strong hard X-ray emission is also detected with ASCA and its spectrum above 3 keV is extremely flat with a prominent iron K line complex, very similar to that seen in the Seyfert 2 galaxy NGC 1068 but about an order of magnitude more luminous ( L _3−10keV ≈ 1.4 × 10⁴² erg s⁻¹). The hard X-ray spectrum indicates that only reflected X-rays of an active galactic nucleus (AGN) buried in a heavy obscuration ( N _H > 2 × 10²⁴ cm⁻²) are visible. This is evidence for an AGN in NGC 6240, emitting possibly at a quasar luminosity (∼ 10⁴⁵ erg s⁻¹), and suggests its significant contribution to the far-infrared luminosity.     

The far-infrared–submillimetre spectral energy distribution of high-redshift quasars

We combine photometric observations of high-redshift     quasars, obtained at submillimetre to millimetre wavelengths, to obtain a mean far-infrared (rest-frame) spectral energy distribution (SED) of the thermal emission from dust, parametrized by a single temperature ( T ) and power-law emissivity index ( β ). The best-fitting values are     and     . Our method exploits the redshift spread of this set of quasars, which allows us to sample the SED at a larger number of rest wavelengths than is possible for a single object: the wavelength range extends down to ∼60 μm, and therefore samples the turnover in the greybody curve for these temperatures. This parametrization is of use to any studies that extrapolate from a flux at a single wavelength, for example to infer dust masses and far-infrared luminosities.
We interpret the cool, submillimetre component as arising from dust heated by star formation in the host galaxy of the quasar, although we do not exclude the presence of dust heated directly by the active galactic nucleus (AGN). Applying the mean SED to the data, we derive consistent star formation rates ∼1000 M_⊙ yr⁻¹ and dust masses ∼10⁹ M_⊙, and investigate a simple scheme of AGN and host galaxy co-evolution to account for these quantities. The time-scale for formation of the host galaxy is     , and the luminous quasar phase occurs towards the end of this period, just before the reservoir of cold gas is depleted. Given the youth of the Universe at     (1.6 Gyr), the coexistence of a massive black hole and a luminous starburst at high redshifts is a powerful constraint on models of quasar host galaxy formation.     

Hyperluminous infrared galaxies

39 galaxies are now known, from follow-up of faint IRAS sources and from submillimetre observations of high-redshift AGN, with far-infrared luminosities >10¹³ L_⊙. 13 of these, which have been found in 60- or 850-μm surveys, form an important unbiased subsample. 12 have been found by comparison of 60-μm surveys with quasar or radio galaxy catalogues, or from infrared surveys with colour selection biased towards AGN, while a further 14 have been found through submillimetre observations of known high-redshift AGN. In this paper I argue, on the basis of detailed modelling of the spectral energy distributions of hyperluminous galaxies with accurate radiative transfer models, and from evidence of high gas mass in several cases, that the bulk of the emission from these galaxies at rest frame wavelengths ≥50 μm is caused by star formation. Even after correction for the effects of lensing, hyperluminous galaxies with emission peaking at rest frame wavelengths ≥50 μm are therefore undergoing star formation at rates >10³ M_⊙ yr⁻¹ and are strong candidates for being primeval galaxies, in the process of a major episode of star formation.     

Far-infrared-loud quasars – I. Disturbed and quiescent quasars in the PG survey

We use host galaxy imaging studies of the PG quasar survey to compare the far-infrared (FIR) properties of quasars with disturbed and undisturbed host galaxies. By using survival analysis, we show that the quasars with disturbed host galaxies, with morphologies classified from a homogeneous data set, have a 60-μm luminosity distribution that is different from that of those with undisturbed hosts with >97 per cent confidence. For morphological classifications using an inhomogeneous data set, including HST data for some objects, this confidence rises to >99 per cent confidence. The mean 60-μm luminosity of the disturbed-host quasars is several times greater than that of the undisturbed-host quasars. However, possible biases in the PG survey might affect these conclusions. Our results are interpreted as supporting the idea that quasars are related to at least some ultraluminous infrared galaxies. We discuss the implications of this result for studies of quasar and galaxy evolution.     

An irradiated accretion disc in the narrow-line Seyfert 1 RE J1034+396?

We model the optical to X-ray continuum spectrum of the narrow-line Seyfert 1 galaxy RE J1034+396. We show that the flat optical spectrum is consistent with emission from an irradiated accretion disc. The X-ray emission can be modelled with a disc blackbody and a Comptonized component. The temperature at the inner edge of the disc     Using this constraint, we show that the flat optical spectrum is consistent with emission from the irradiatively heated outer part of the accretion disc. We constrain the outer radius of the optically thick disc     and the inner radius of the irradiation-dominated region     . Our optical and X-ray spectral fits indicate a mass     , and do not rule out a low (i.e. face-on) inclination angle for the system.     

The radio properties of radio-quiet quasars

Although radio-quiet quasars (RQQs) constitute ≳ 90 per cent of optically identified quasar samples, their radio properties are only poorly understood. In this paper we present the results of a multi-frequency VLA study of 27 low-redshift RQQs. We detect radio emission from 20 objects, half of which are unresolved (≤ 0.24 arcsec). In cases where significant structure can be resolved, double, triple and linear radio sources on scales of a few kpc are found. The radio emission (typically) has a steep spectrum (α ∼ 0.7, where S  ∝ ν^−α), and high brightness temperatures ( T _B ≥ 10⁵ K) are measured in some of the radio components. The RQQs form a natural extension to the radio luminosity–absolute magnitude distribution of nearby Seyfert 1s. We conclude that a significant fraction of the radio emission in RQQs originates in a compact nuclear source directly associated with the quasar. There are no significant differences between the radio properties of RQQs with elliptical hosts and those in disc galaxies within the current sample.     

Red companions to a z=2.15 radio-loud quasar

We have conducted observations of the environment around the z =2.15 radio-loud quasar 1550–269 in search of distant galaxies associated either with it or the z =2.09 C  iv absorber along its line of sight. Such objects will be distinguished by their red colours, R − K >4.5. We find five such objects in a 1.5 arcmin² field around the quasar, with typical K ' magnitudes of ∼20.4 and no detected R -band emission. We also find a sixth object with K =19.6±0.3, and undetected at R , just two arcsec from the quasar. The nature of all these objects is currently unclear, and will remain so until we have determined their redshifts. We suggest that it is likely that they are associated with either the quasar or the C  iv absorber, in which case their properties might be similar to those of the z =2.38 red Ly α emitting galaxies discovered by Francis et al. The small separation between the quasar and the brightest of our objects suggests that it may be the galaxy responsible for the C  iv metal line absorption system. The closeness to the quasar and the red colour might have precluded similar objects from being uncovered in previous searches for emission from C  iv and damped absorbers.     

The unified model and the Seyfert 2 infrared dichotomy

We present X-ray imaging spectroscopy of the extremely luminous infrared galaxy IRAS 09104+4109     obtained with the Chandra X-ray Observatory. With the arcsec resolution of Chandra , an unresolved source at the nucleus is separated from the surrounding cluster emission. A strong iron K line at 6.4 keV on a very hard continuum is detected from the nuclear source, rendering IRAS 09104+4109 the most distant reflection-dominated X-ray source known. Combined with the BeppoSAX detection of the excess hard X-ray emission, it provides further strong support for the presence of a hidden X-ray source of quasar luminosity in this infrared galaxy. Also seen is a faint linear structure to the north, which coincides with the main radio jet. An X-ray deficit in the corresponding region suggests an interaction between the cluster medium and the jet driven by the active nucleus.     

Emission-line outflows in PKS1549−79: the effects of the early stages of radio-source evolution?

We present new spectroscopic observations of the southern radio galaxy         . Despite the flat-spectrum character of the radio emission from this source, our optical spectra show no sign of the broad permitted lines and non-stellar continuum characteristic of quasar nuclei and broad-line radio galaxies. However, the high-ionization forbidden lines, including [O  iii ] λλ 5007, 4959, are unusually broad for a narrow-line radio galaxy     , and are blueshifted by 600 km s⁻¹ relative to the low-ionization lines such as [O  ii ] λλ 3726,3729. The [O  ii ] lines are also considerably narrower     than the [O  iii ] lines, and have a redshift consistent with that of the recently detected H  i 21-cm absorption-line system. Whereas the kinematics of the [O  iii ] emission lines are consistent with outflow in an inner narrow-line region, the properties of the [O  ii ] emission lines suggest that they are emitted by a more extended and quiescent gaseous component. We argue that, given the radio properties of the source, our line of sight is likely to be lying close to the direction of bulk outflow of the radio jets. In this case it is probable that the quasar nucleus is entirely obscured at optical wavelengths by the material responsible for the H  i absorption-line system. The unusually broad [O  iii ] emission lines suggest that the radio source is intrinsically compact. Overall, our data are consistent the idea that     is a radio source in an early stage of evolution.     

Modelling the cosmological co-evolution of supermassive black holes and galaxies – II. The clustering of quasars and their dark environment

We use semi-analytic modelling on top of the Millennium simulation to study the joint formation of galaxies and their embedded supermassive black holes. Our goal is to test scenarios in which black hole accretion and quasar activity are triggered by galaxy mergers, and to constrain different models for the light curves associated with individual quasar events. In the present work, we focus on studying the spatial distribution of simulated quasars. At all luminosities, we find that the simulated quasar two-point correlation function is fit well by a single power law in the range  0.5 ≲ r ≲ 20  h ⁻¹ Mpc  , but its normalization is a strong function of redshift. When we select only quasars with luminosities within the range typically accessible by today's quasar surveys, their clustering strength depends only weakly on luminosity, in agreement with observations. This holds independently of the assumed light-curve model, since bright quasars are black holes accreting close to the Eddington limit, and are hosted by dark matter haloes with a narrow mass range of a few  10¹²  h ⁻¹ M_⊙  . Therefore, the clustering of bright quasars cannot be used to disentangle light-curve models, but such a discrimination would become possible if the observational samples can be pushed to significantly fainter limits. Overall, our clustering results for the simulated quasar population agree rather well with observations, lending support to the conjecture that galaxy mergers could be the main physical process responsible for triggering black hole accretion and quasar activity.     

The CO content of the Local Group dwarf irregular galaxies IC 5152, UGCA 438 and the Phoenix dwarf

We present a search for  CO(1 → 0)  emission in three Local Group dwarf irregular galaxies: IC 5152, the Phoenix dwarf and UGCA 438, using the ATNF Mopra radio telescope. Our scans largely cover the optical extent of the galaxies and the stripped H  i cloud west of the Phoenix dwarf. Apart from a tentative but non-significant emission peak at one position in the Phoenix dwarf, no significant emission was detected in the CO spectra of these galaxies. For a velocity width of 6 km s⁻¹, we derive 4σ upper limits of 0.03, 0.04 and 0.06 K km s⁻¹ for IC 5152, the Phoenix dwarf and UGCA 438, respectively. This is an improvement of over a factor of 10 compared with previous observations of IC 5152; the other two galaxies had not yet been observed at millimetre wavelengths. Assuming a Galactic CO-to-H₂ conversion factor, we derive upper limits on the molecular gas mass of  6.2 × 10⁴, 3.7 × 10³  and  1.4 × 10⁵ M_⊙  for IC 5152, the Phoenix dwarf and UGCA 438, respectively. We investigate two possible causes for the lack of CO emission in these galaxies. On the one hand, there may be a genuine lack of molecular gas in these systems, in spite of the presence of large amounts of neutral gas. However, in the case of IC 5152 which is actively forming stars, molecular gas is at least expected to be present in the star-forming regions. On the other hand, there may be a large increase in the CO-to-H₂ conversion factor in very low-metallicity dwarfs  (−2 ≤[Fe/H]≤−1)  , making CO a poor tracer of the molecular gas content in dwarf galaxies.     

XMM observations of the high-redshift quasar RX J1028.6 – 0844 at z= 4.276 : soft X-ray spectral flattening

We present results from a new XMM–Newton observation of the high-redshift quasar RX J1028.6 – 0844 at a redshift of 4.276. The soft X-ray spectral flattening, as reported by a previous study with ASCA , is confirmed to be present, with, however, a reduced column density when modelled by absorption. The inferred column density for absorption intrinsic to the quasar is  2.1(^+0.4_−0.3) × 10²²  cm⁻²  for cold matter, and higher for ionized gas. The spectral flattening shows remarkable similarity with that of two similar object, namely GB 1428 + 4217 and PMN J0525 − 3343. The results improve upon those obtained from a previous short-exposure observation for RX J1028.6 – 0844 with XMM–Newton . A comparative study of the two XMM–Newton observations reveals a change in the power-law photon index from  Γ≃ 1.3  to 1.5 on time-scales of about one year. A tentative excess emission feature in the rest-frame 5–10 keV band is suggested, which is similar to that marginally suggested for GB 1428 + 4217.     

The quasar fraction in low-frequency-selected complete samples and implications for unified schemes

Low-frequency radio surveys are ideal for selecting orientation-independent samples of extragalactic sources because the sample members are selected by virtue of their isotropic steep-spectrum extended emission. We use the new 7C Redshift Survey along with the brighter 3CRR and 6C samples to investigate the fraction of objects with observed broad emission lines – the 'quasar fraction'– as a function of redshift and of radio and narrow-emission-line luminosity. We find that the quasar fraction is more strongly dependent upon luminosity (both narrow-line and radio) than it is on redshift. Above a narrow [O  ii ] emission-line luminosity of log₁₀( L _{[O  ii ]}/W)≳35 [or radio luminosity log₁₀( L ₁₅₁/W Hz⁻¹ sr⁻¹)≳ 26.5], the quasar fraction is virtually independent of redshift and luminosity; this is consistent with a simple unified scheme with an obscuring torus with a half-opening angle θ _trans≈53°. For objects with less luminous narrow lines, the quasar fraction is lower. We show that this is not due to the difficulty of detecting lower luminosity broad emission lines in a less luminous, but otherwise similar, quasar population. We discuss evidence which supports at least two probable physical causes for the drop in quasar fraction at low luminosity: (i) a gradual decrease in θ _trans and/or a gradual increase in the fraction of lightly reddened (0≲ A _V ≲5) lines of sight with decreasing quasar luminosity; and (ii) the emergence of a distinct second population of low-luminosity radio sources which, like M87, lack a well-fed quasar nucleus and may well lack a thick obscuring torus.     

ASCA and ROSAT observations of the Seyfert 1 galaxy RX J0437.4−4711

Results of ASCA and ROSAT observations of the Seyfert 1 galaxy RX J0437.4−4711 are presented. The X-ray continuum spectrum can be described by the sum of a power law with photon index 2.15 ± 0.04 and a soft emission component characterized by a blackbody with temperature 29 ± 2 eV. The total luminosity of the soft component is larger than that of the power-law component if the power law is cut off around a few hundred keV. A weak absorption edge with τ = 0.26 ± 0.13 at the rest-frame energy of E  = 0.83 ± 0.05 keV and an Fe Kα line with EW = 430 ± 220 eV at an energy E  = 6.47 ± 0.15 keV are also detected. The X-ray flux showed a 47 per cent increase between two ASCA observations 4 months apart, but no spectral variability was seen. We argue that reprocessing of the hard X-ray emission cannot produce all the soft X-ray emission, since the total luminosity of the soft component is larger than that of the integrated power-law component. Similarities with some stellar black hole candidates are briefly discussed.     

A high-velocity ionized outflow and XUV photosphere in the narrow emission line quasar PG1211+143

We report on the analysis of a ∼60-ks XMM–Newton observation of the bright, narrow emission line quasar PG1211+143. Absorption lines are seen in both European Photon Imaging Camera and Reflection Grating Spectrometer spectra corresponding to H- and He-like ions of Fe, S, Mg, Ne, O, N and C. The observed line energies indicate an ionized outflow velocity of ∼24 000 km s⁻¹. The highest energy lines require a column density of   N _H∼ 5 × 10²³ cm⁻²  , at an ionization parameter of  log ξ∼ 3.4  . If the origin of this high-velocity outflow lies in matter being driven from the inner disc, then the flow is likely to be optically thick within a radius of ∼130 Schwarzschild radii, providing a natural explanation for the big blue bump (and strong soft X-ray) emission in PG1211+143.     

The infrared jet in 3C 66B

We present images of infrared (IR) emission from the radio jet in 3C 66B. Data at three wavelengths (4.5, 6.75 and 14.5 μm) were obtained using the Infrared Space Observatory . The 6.75-μm image clearly shows an extension aligned with the radio structure. The jet was also detected in the 14.5-μm image, but not at 4.5 μm. The radio–infrared–optical spectrum of the jet can be interpreted as synchrotron emission from a population of electrons with a high-energy break of 4×10¹¹ eV. We place upper limits on the IR flux from the radio counter-jet. A symmetrical, relativistically beamed twin-jet structure is consistent with our results if the jets consist of multiple components.     

Contribution of galaxies to the background hydrogen-ionizing flux

We estimate the evolution of the contribution of galaxies to the cosmic background flux at 912 Å by means of a semi-analytic model of galaxy formation and evolution. Such modelling has been quite successful in reproducing the optical properties of galaxies. We assume that high-redshift damped Lyman α systems are the progenitors of present-day galaxies, and we design a series of models that are consistent with the evolution of cosmic comoving emissivities in the available near-infrared, optical, ultraviolet and far-infrared bands along with the evolution of the neutral hydrogen content and average metallicity of damped Lyman α systems. We use these models to compute the galactic contribution to the Lyman-limit emissivity and background flux for 0 ≃  z  ≤ 4. We take into account the absorption of Lyman-limit photons by H  I and dust in the interstellar medium of the galaxies. We find that the background Lyman-limit flux due to galaxies might dominate (or be comparable to) the contribution from quasars at almost all redshifts if the absorption by H  I in the interstellar medium is neglected. Such H  I absorption would result in a severe diminishing of this flux — by almost three orders of magnitude at high redshifts and by one to two orders at z  ≃ 0. Though the resulting galaxy flux is completely negligible at high redshifts, it is comparable to the quasar flux at z  ≃ 0.     

Bayesian inference of jet bulk-flow speeds in Fanaroff–Riley type II radio sources

Radio jet and core data for a complete sample of 98 Fanaroff–Riley type II sources with   z < 1  are analysed with a Markov Chain Monte Carlo model fitting method to obtain constraints on bulk-flow speeds in the beam. The Bayesian parameter-inference method is described and demonstrated to be capable of providing meaningful constraints on the Lorentz factor at both kpc and parsec scales. For both jets and cores, we show that models in which some intrinsic dispersion is present in the features' intrinsic prominence, bulk-flow speeds or both provide the best fit to the data. The constraints on the Lorentz factor on parsec scales are found to be consistent with the expected values given very long baseline interferometry observations and other evidence, with     . On kpc scales, the Lorentz factor is found to be ≈1.18–1.49, in agreement with the results of previous analyses of radio jet data. These values are clearly not consistent with the  γ≈ 10  speeds required by beamed inverse-Compton models of X-ray emission from quasar jets; our results therefore support models that require velocity structure in powerful jets.