Weighing a galaxy bar in the lens Q2237 +=0305

In the gravitational lens system Q2237+0305 the cruciform quasar image geometry is twisted by 10c by the lens effect of a bar in the lensing galaxy. This effect can be used to measure the mass of the bar. We construct a new lensing model for this system with a power-law elliptical bulge and a Ferrers bar. The observed ellipticity of the optical isophotes of the galaxy leads to a nearly isothermal elliptical profile for the bulge, with a total quasar magnification of 16⁺⁵₋₄. We measure a bar mass of (7.5 ∼ 1.5) −10⁸ h ⁻¹₇₅ M⊙ in the region inside the quasar images.     

Constraining quasar host halo masses with the strength of nearby Lyα forest absorption

Using cosmological hydrodynamic simulations, we measure the mean transmitted flux in the Lyα forest for quasar sightlines that pass near a foreground quasar. We find that the trend of absorption with pixel quasar separation distance can be fitted using a simple power-law form including the usual correlation function parameters r ₀ and γ, so that     . From the simulations, we find the relation between r ₀ and quasar host mass, and formulate this as a way to estimate quasar host dark matter halo masses, quantifying uncertainties due to cosmological and IGM parameters, and redshift errors. With this method, we examine data for ∼9000 quasars from the Sloan Digital Sky Survey (SDSS) Data Release 5, assuming that the effect of ionizing radiation from quasars (the so-called transverse proximity effect) is unimportant (no evidence for it is seen in the data). We find that the best-fitting host halo mass for SDSS quasars with mean redshift z = 3 and absolute G -band magnitude −27.5 is  log  M /M_⊙= 12.68^+0.81_−0.67  . We also use the Lyman-Break Galaxy (LBG) and Lyα forest data of Adelberger et al. in a similar fashion to constrain the halo mass of LBGs to be  log₁₀  M /M_⊙= 11.41^+0.54_−0.59  , a factor of ∼20 lower than the bright quasars. In addition, we study the redshift distortions of the Lyα forest around quasars, using the simulations. We use the quadrupole to monopole ratio of the quasar Lyα forest correlation function as a measure of the squashing effect. We find its dependence on halo mass difficult to measure, but find that it may be useful for constraining cosmic geometry.     

Non-parametric strong lens inversion of SDSS J1004+4112

In this article, we study the well-known strong lensing system SDSS J1004+4112. Not only does it host a large-separation lensed quasar with measured time-delay information, but several other lensed galaxies have been identified as well. A previously developed strong lens inversion procedure that is designed to handle a wide variety of constraints is applied to this lensing system and compared to results reported in other works. Without the inclusion of a tentative central image of one of the galaxies as a constraint, we find that the model recovered by the other constraints indeed predicts an image at that location. An inversion which includes the central image provides tighter constraints on the shape of the central part of the mass map. The resulting model also predicts a central image of a second galaxy where indeed an object is visible in the available Advanced Camera for Surveys images. We find masses of  2.5 × 10¹³  and  6.1 × 10¹³ M_⊙  within a radius of 60 and 110 kpc, respectively, confirming the results from other authors. The resulting mass map is compatible with an elliptical generalization of a projected NFW profile, with   r _s= 58⁺²¹₋₁₃  arcsec and   c _vir= 3.91 ± 0.74  . The orientation of the elliptical NFW profile closely follows the orientation of the central cluster galaxy and the overall distribution of cluster members.     

Implications for unified schemes from submillimetre and far-infrared follow-up of radio-selected samples

We extend our previous analysis which used generalized luminosity functions (GLFs) to predict the number of quasars and galaxies in low-radio-frequency-selected samples as a function of redshift, radio luminosity, narrow-emission-line luminosity and type of unified scheme. Our extended analysis incorporates the observed submillimetre (850-μm) flux densities of radio sources, employs a new method which allows us to deal with non-detections, and focuses on the high-luminosity population. First, we conclude that the submillimetre luminosity L ₈₅₀ of low-frequency-selected radio sources is correlated with the bolometric luminosity L _bol of their quasar nuclei via an approximate scaling relation   L ₈₅₀∝ L ^0.7±0.2_bol  . Secondly, we conclude that there is quantitative evidence for a receding-torus-like physical process for the high-luminosity population within a two-population unified scheme for radio sources; this evidence comes from the fact that radio quasars are brighter in both narrow emission lines and submillimetre luminosity than radio galaxies matched in radio luminosity and redshift. Thirdly, we note that the combination of a receding-torus-like scheme and the assumption that the observed submillimetre emission is dominated by quasar-heated dust yields a scaling relation   L ₈₅₀∝ L ^1/2_bol  which is within the errors of that determined here for radio-selected quasars, and consistent with that inferred for radio-quiet quasars.     

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.     

Constraints on a quintessence model from gravitational lensing statistics

Constraints on an exact quintessence scalar-field model with an exponential potential are derived from gravitational lens statistics. An exponential potential can account for data from both optical quasar surveys and radio-selected sources. Based on the Cosmic Lens All-Sky Survey (CLASS) sample, lensing statistics provides, for the pressureless matter density parameter, an estimate of  Ω_M0= 0.31^+0.12_−0.14  .     

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.     

Discovery of radio-loud quasars with z=4.72 and z=4.010

We report the discovery of two radio-loud quasars with redshifts greater than 4: GB1428+4217, with z =4.72, and GB1713+2148 with z =4.01. This doubles the number of published radio-selected quasars with z >4, bringing the total to 4. GB1428+4217 is the third most distant quasar known and the highest redshift radio and X-ray source currently known. It has a radio flux density at 5 GHz of 259±31 mJy and an optical magnitude of R ∼20.9. The rest frame absolute UV magnitude, M_v (1450 Å), is −26.7, similar to that of the archetypal radio-selected quasar 3C273 [ z =0.158; M_v (1450 Å)=−26.4]. GB1428+4217 is tentatively detected in ROSAT PSPC observations, which has been confirmed by more recent ROSAT observations described in a companion paper by Fabian et al. Both quasars were discovered during the CCD imaging phase of an investigation into the evolution of the space density of radio-loud quasars at high redshift. Combined with our earlier survey results, these objects give a lower limit on the space density of quasars with radio power P _{5 GHz}>5.8×10²⁶ W Hz⁻¹ sr⁻¹ between z =4 and z =5 of 1.4±0.9×10⁻¹⁰ Mpc⁻³. This can be compared to 2.9±0.2×10⁻¹⁰ Mpc⁻³ at z =2 from Dunlop & Peacock for flat-spectrum sources of the same luminosity.     

On the streaming motions of haloes and galaxies

Little is known about the statistics of gravitationally lensed quasars at large (7–30 arcsec) image separations, which probe masses on the scale of galaxy clusters. We have carried out a survey for gravitationally lensed objects, among sources in the FIRST 20-cm radio survey that have unresolved optical counterparts in the digitizations of the Palomar Observatory Sky Survey. From the statistics of ongoing surveys that search for quasars among FIRST sources, we estimate that there are about 9100 quasars in this source sample, making this one of the largest lensing surveys to date. Using broad-band imaging, we have isolated all objects with double radio components separated by 5–30 arcsec that have unresolved optical counterparts with similar BVI colours. Our criteria for similar colours conservatively allow for observational error and for colour variations due to time delays between lensed images. Spectroscopy of these candidates shows that none of the pairs are lensed quasars. This sets an upper limit (95 per cent confidence) on the lensing fraction in this survey of 3.3×10⁻⁴, assuming 9100 quasars. Although the source redshift distribution is poorly known, a rough calculation of the expected lensing frequency and the detection efficiencies and biases suggests that simple theoretical expectations are of the same order of magnitude as our observational upper limit. Our procedure is novel in that our exhaustive search for lensed objects does not require prior identification of the quasars in the sample as such. Characterization of the FIRST-selected quasar population will enable use of our result to constrain quantitatively the mass properties of clusters.     

Separation of the visible and dark matter in the Einstein ring LBG J213512.73−010143

We model the mass distribution in the recently discovered Einstein ring LBG J213512.73−010143 (the 'Cosmic Eye') using archival Hubble Space Telescope imaging. We reconstruct the mass density profile of the z = 0.73 lens and the surface brightness distribution of the z = 3.07 source and find that the observed ring is best fitted with a dual-component lens model consisting of a baryonic Sersic component nested within a dark matter halo. The dark matter halo has an inner slope of 1.42^+0.24_−0.22, consistent with cold dark matter simulations after allowing for baryon contraction. The baryonic component has a mass-to-light ratio of  1.71^+0.28_−0.38 M_⊙/L _{B ⊙}  which when evolved to the present day is in agreement with local ellipticals. Within the Einstein radius of 0.77 arcsec (5.6 kpc), the baryons account for 46 ± 11 per cent of the projected lens mass. External shear from a nearby foreground cluster is accurately predicted by the model. The reconstructed surface brightness distribution in the source plane clearly shows two peaks. Through a generalization of our lens inversion method, we conclude that the redshifts of both peaks are consistent with each other, suggesting that we are seeing structure within a single galaxy.     

Clustering of galaxies around radio quasars at 0.5≤z≤0.8

We have observed the galaxy environments around a sample of 21 radio-loud, steep-spectrum quasars at 0.5≤ z ≤0.82, spanning several orders of magnitude in radio luminosity. The observations also include background control fields used to obtain the excess number of galaxies in each quasar field. The galaxy excess was quantified using the spatial galaxy–quasar correlation amplitude, B _gq, and an Abell-type measurement, N _0.5. A few quasars are found in relatively rich clusters, but on average, they seem to prefer galaxy groups or clusters of approximately Abell class 0. We have combined our sample with literature samples extending down to z ≈0.2 and covering the same range in radio luminosity. By using the Spearman statistic to disentangle redshift and luminosity dependences, we detect a weak, but significant, positive correlation between the richness of the quasar environment and the radio luminosity of the quasar. However, we do not find any epoch dependence in B _gq, as has previously been reported for radio quasars and galaxies. We discuss the radio luminosity–cluster richness link and possible explanations for the weak correlation that is seen.     

Interpreting the transmission windows of distant quasars

We propose the apparent shrinking criterion (ASC) to interpret the spatial extent, R _w, of transmitted flux windows in the absorption spectra of high- z quasars. The ASC can discriminate between the two regimes in which R _w corresponds either to the physical size,   R _{H  ii}   , of the quasar H  ii region or to the distance,   R ^max_w  , at which the transmitted flux drops to  =0.1  and a Gunn–Peterson (GP) trough appears. In the first case [H  ii region (HR) regime], one can determine the intergalactic medium mean H  i fraction,   x _H I  ; in the second [proximity region (PR) regime], the value of R _w allows one to measure the local photoionization rate and the local enhancement of the photoionization rate,  Γ_G  , due to nearby/intervening galaxies. The ASC has been tested against radiative transfer+smoothed particle hydrodynamics numerical simulations, and applied to a sample of 15 high-   z ( z > 5.8  ) quasar spectra. All sample quasars are found to be in the PR regime; hence, their observed spectral properties (inner flux profile, extent of transmission window) cannot reliably constrain the value of   x _{H  i}   . Four sample quasars show evidence for a local enhancement (up to 50 per cent) in the local photoionization rate possibly produced by a galaxy overdensity. We discuss the possible interpretations and uncertainties of this result.     

The quasar-frame velocity distribution of narrow C iv absorbers

We report on a survey for narrow (full widths at half-minimum <600 km s⁻¹) C  iv absorption lines in a sample of bright quasars at redshifts  1.8 ≤ z < 2.25  in the Sloan Digital Sky Survey. Our main goal is to understand the relationship of narrow C  iv absorbers to quasar outflows and, more generally, to quasar environments. We determine velocity zero-points using the broad Mg  ii emission line, and then measure the absorbers' quasar-frame velocity distribution. We examine the distribution of lines arising in quasar outflows by subtracting model fits to the contributions from cosmologically intervening absorbers and absorption due to the quasar host galaxy or cluster environment. We find that a substantial number (  ≥43 ± 6  per cent) of absorbers with   W ^λ1548₀ > 0.3  Å in the velocity range  +750 ≲ v ≲+ 12 000  km s⁻¹ are intrinsic to the active galactic nucleus outflow. This 'outflow fraction' peaks near   v =+2000  km s⁻¹ with a value of   f _outflow≃ 0.81 ± 0.13  . At velocities below   v ≈+ 2000  km s⁻¹, the incidence of outflowing systems drops, possibly due to geometric effects or to the over-ionization of gas that is nearer the accretion disc. Furthermore, we find that outflow absorbers are on average broader and stronger than cosmologically intervening systems. Finally, we find that ∼14 per cent of the quasars in our sample exhibit narrow, outflowing C  iv absorption with   W ^λ1548₀ > 0.3  Å, slightly larger than that for broad absorption line systems.     

A simple model for the evolution of supermassive black holes and the quasar population

An empirically motivated model is presented for accretion-dominated growth of supermassive black holes (SMBH) in galaxies, and the implications are studied for the evolution of the quasar population in the Universe. We investigate the core aspects of the quasar population, including space density evolution, evolution of the characteristic luminosity, plausible minimum masses of quasars, the mass function of SMBH and their formation epoch distribution. Our model suggests that the characteristic luminosity in the quasar luminosity function arises primarily as a consequence of a characteristic mass scale above which there is a systematic separation between the black hole and the halo merging rates. At lower mass scales, black hole merging closely tracks the merging of dark haloes. When combined with a declining efficiency of black hole formation with redshift, the model can reproduce the quasar luminosity function over a wide range of redshifts. The observed space density evolution of quasars is well described by formation rates of SMBH above  ∼10⁸  M_⊙  . The inferred mass density of SMBH agrees with that found independently from estimates of the SMBH mass function derived empirically from the quasar luminosity function.     

Multiwavelength observations of serendipitous Chandra X-ray sources in the field of A 2390

We present optical spectra and near-infrared imaging of a sample of 31 serendipitous X-ray sources detected in the field of Chandra observations of the A 2390 cluster of galaxies. The sources have  0.5–7 keV  fluxes of  (0.6–8)×10^-14 erg cm^-2 s^-1  and lie around the break in the  2–10 keV  source counts. They are therefore typical of sources dominating the X-ray Background in that band. 12 of the 15 targets for which we have optical spectra show emission lines at a range of line luminosities, and half of these show broad lines. These active galaxies and quasars have soft X-ray spectra. Including photometric redshifts and published spectra, we have redshifts for 17 of the sources, ranging from   z ∼0.2  up to   z ∼3  , with a peak between   z =1–2  . 10 of our sources have hard X-ray spectra indicating a spectral slope flatter than that of a typical unabsorbed quasar. Two hard sources that are gravitationally lensed by the foreground cluster are obscured quasars, with intrinsic  2–10 keV  luminosities of  (0.2–3)×10⁴⁵ erg s^-1  , and absorbing columns of   N _H>10²³ cm^-2  . Both of these sources were detected in the mid-infrared by ISOCAM on the Infrared Space Observatory , which when combined with radiative transfer modelling leads to the prediction that the bulk of the reprocessed flux emerges at ∼100 μm.     

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.     

HST Planetary Camera images of quasar host galaxies

We present Hubble Space Telescope ( HST ) images of seven low-redshift quasars (six taken with the Planetary Camera, one with the Wide Field Camera). These complete the sample of 14 quasars observed by the Faint Object Camera Investigation Definition Team (FOC IDT). Following subtraction of the quasar nuclear light, host galaxies can be seen in all seven cases. A combination of the optical morphology and luminosity profiles of the residual host galaxies and the results of 2D cross-correlation model fitting implies that five of the objects have elliptical host galaxies and two have disc host galaxies. The luminosities vary from slightly fainter than L ^* to about 1.3 mag brighter than L ^*.   We discuss the properties of the complete sample of 14 quasars. Nine of the objects appear to have elliptical host galaxies (all six of the radio-loud quasars in the sample as well as three radio-quiet quasars). Two further radio-quiet quasars appear to lie in disc galaxies. The other three objects (radio-quiet, ultraluminous infrared quasars) all lie in violently interacting systems. The sample as a whole has an average luminosity about 0.8 mag brighter than L ^*, although the radio-loud objects have hosts on average 0.7 mag brighter than the radio-quiet objects.   We compare our results with those from HST imaging of quasars by other authors. Taken together, our observations are in broad agreement with those of Bahcall et al. Radio-loud quasars appear to lie in luminous elliptical galaxies whereas radio-quiet quasars are found to lie in either elliptical or spiral hosts. Host galaxy luminosities (of radio-quiet and radio-loud quasars) are much brighter than would be expected if they followed a Schechter luminosity function.     

PDS 456: an extreme accretion rate quasar?

We present quasi-simultaneous ASCA and RXTE observations of the most luminous known active galactic nucleus in the local ( z <0.3) Universe, the recently discovered quasar PDS 456. Multiwavelength observations have been conducted that show that PDS 456 has a bolometric luminosity of ∼10⁴⁷ erg s⁻¹ peaking in the ultraviolet part of the spectrum. In the X-ray band the 2–10 keV (rest-frame) luminosity is 10⁴⁵ erg s⁻¹. The broad-band X-ray spectrum obtained with ASCA and RXTE contains considerable complexity. The most striking feature observed is a very deep, ionized iron K edge, observed at 8.7 keV in the quasar rest-frame. We find that these features are consistent with reprocessing from highly ionized matter, probably the inner accretion disc. PDS 456 appeared to show a strong (factor of ∼2.1) outburst in just ∼17 ks, although non-intrinsic sources cannot be completely ruled out. If confirmed, this would be an unusual event for such a high-luminosity source, with a light-crossing-time corresponding to ∼2 R _S . The implication would be that flaring occurs within the very central regions, or else that PDS 456 is a 'super-Eddington' or relativistically beamed system. Overall we conclude on the basis of the extreme blue/UV luminosity, the rapid X-ray variability and from the imprint of highly ionized material on the X-ray spectrum, that PDS 456 is a quasar with an unusually high accretion rate.     

Effects of quasar feedback in galaxy groups

We study the effect of quasar feedback on distributions of baryons in galaxy groups using high-resolution numerical simulations. We use the entropy-conserving gadget code that includes gas cooling and star formation, modified to include a physically based model of quasar feedback. For a sample of 10 galaxy group-sized dark matter haloes with masses in the range of  1–5 × 10¹³ M_⊙  h ⁻¹  , star formation is suppressed by more than 50 per cent in the inner regions due to the additional pressure support by quasar feedback, while gas is driven from the inner region towards the outer region of the haloes. As a result, the average gas density is 50 per cent lower in the inner region and 10 per cent higher in the outer region in the simulation, compared to a similar simulation with no quasar feedback. Gas pressure is also higher in the outer region, while temperature and entropy are enhanced in the inner region. The total group gas fraction in the two simulations generally differs by less than 10 per cent. We also find a small change of the total thermal Sunyaev–Zeldovich distortion, leading to 10 per cent changes in the microwave angular power spectrum at angular scales below 2 arcmin.     

Infrared photometry of z ∼ 1 3C quasars

We present JHKL ' photometry of a complete sample of steep-spectrum radio-loud quasars from the revised 3CR catalogue in the redshift range 0.65 z <1.20. After correcting for contributions from emission lines and the host galaxies, we investigate their spectral energy distributions (SEDs) around 1 μm. About 75 per cent of the quasars are tightly grouped in the plane of optical spectral index, α _opt, versus near-infrared spectral index, α _IR, with the median value of α _opt close to the canonical value, and the median α _IR slightly flatter. We conclude that the fraction of moderately obscured, red quasars decreases with increasing radio power, in accordance with the 'receding torus' model which can also explain the relatively flat median near-infrared spectra of the 3CR quasars. Two of the red quasars have inverted infrared spectral indices, and we suggest that their unusual SEDs might result from a combination of dust-scattered and transmitted quasar light.