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.     

High-resolution observations and mass modelling of the CLASS gravitational lens B1152+199

We present a series of high-resolution radio and optical observations of the CLASS gravitational lens system B1152+199 obtained with the Multi-Element Radio-Linked Interferometer Network, Very Long Baseline Array and Hubble Space Telescope . Based on the milliarcsecond-scale substructure of the lensed radio components and precise optical astrometry for the lensing galaxy, we construct models for the system and place constraints on the galaxy mass profile. For a single galaxy model with surface mass density  Σ(r)∝r^{− β}   , we find that  0.95 β 1.21  at 2 σ confidence. Including a second deflector to represent a possible satellite galaxy of the primary lens leads to slightly steeper mass profiles.     

Feast and Famine: regulation of black hole growth in low-redshift galaxies

We analyse the observed distribution of Eddington ratios  ( L / L _Edd)  as a function of supermassive black hole mass for a large sample of nearby galaxies drawn from the Sloan Digital Sky Survey. We demonstrate that there are two distinct regimes of black hole growth in nearby galaxies. The first is associated with galaxies with significant star formation [   M _*/star formation rate (SFR) ∼  a Hubble time] in their central kiloparsec regions, and is characterized by a broad lognormal distribution of accretion rates peaked at a few per cent of the Eddington limit. In this regime, the Eddington ratio distribution is independent of the mass of the black hole and shows little dependence on the central stellar population of the galaxy. The second regime is associated with galaxies with old central stellar populations (   M _*/SFR ≫  a Hubble time), and is characterized by a power-law distribution function of Eddington ratios. In this regime, the time-averaged mass accretion rate on to black holes is proportional to the mass of stars in the galaxy bulge, with a constant of proportionality that depends on the mean stellar age of the stars. This result is once again independent of black hole mass. We show that both the slope of the power law and the decrease in the accretion rate on to black holes in old galaxies are consistent with population synthesis model predictions of the decline in stellar mass loss rates as a function of mean stellar age. Our results lead to a very simple picture of black hole growth in the local Universe. If the supply of cold gas in a galaxy bulge is plentiful, the black hole regulates its own growth at a rate that does not further depend on the properties of the interstellar medium. Once the gas runs out, black hole growth is regulated by the rate at which evolved stars lose their mass.     

The lens and source of the optical Einstein ring gravitational lens ER 0047–2808

We perform a detailed analysis of the optical gravitational lens ER 0047–2808 imaged with the Wide Field Planetary Camera 2 on the Hubble Space Telescope . Using software specifically designed for the analysis of resolved gravitational lens systems, we focus on how the image alone can constrain the mass distribution in the lens galaxy. We find that the data are of sufficient quality to strongly constrain the lens model with no a priori assumptions about the source. Using a variety of mass models, we find statistically acceptable results for elliptical isothermal-like models with an Einstein radius of 1.17 arcsec. An elliptical power-law model  (Σ∝ R ^−β)  for the surface mass density favours a slope slightly steeper than isothermal with  β= 1.08 ± 0.03  . Other models including a constant mass-to-light ratio (M/L), pure Navarro, Frenk & White halo and (surprisingly) an isothermal sphere with external shear are ruled out by the data. We find the galaxy light profile can only be fit with a Sérsic plus point-source model. The resulting total  M/L _B   contained within the images is  4.7  h ₆₅± 0.3  . In addition, we find the luminous matter is aligned with the total mass distribution within a few degrees. This is the first time a resolved optical gravitational lens image has been quantitatively reproduced using a non-parametric source.
The source, reconstructed by the software, is revealed to have two bright regions, with an unresolved component inside the caustic and a resolved component straddling a fold caustic. The angular size of the entire source is ∼0.1 arcsec and its (unlensed) Lyα flux is  3 × 10⁻¹⁷ erg s⁻¹ cm⁻²  .     

Faraday ghosts: depolarization canals in the Galactic radio emission

An analysis is presented of the power spectrum of X-ray variability of the bright Seyfert 1 galaxy Mrk 766 as observed by XMM–Newton . Over the 0.2–10 keV energy range the power spectral density (PSD) is well-represented by a power-law with a slope of  α_low≈ 1  at low frequencies, breaking to a slope of  α_hi= 2.8^+0.2_−0.4  at a frequency   f _br≈ 5 × 10⁻⁴ Hz  . As has been noted before, this broken power-law PSD shape is similar to that observed in the Galactic black hole candidate Cygnus X-1. If it is assumed that Mrk 766 shows a power spectrum similar in form to that of Cyg X-1, and that the break time-scale scales linearly with black hole mass, then the mass of the black hole in Mrk 766 is inferred to be  ≲ 5 × 10⁵ M_⊙  . This rather low mass would mean Mrk 766 radiates above the Eddington limit. The coherence between different energy bands is significantly below unity implying that variations in the different energy bands are rather poorly correlated. The low coherence can be explained in the framework of standard Comptonization models if the properties of the Comptonizing medium are rapidly variable or if there are several distinct emission sites.     

On the galaxy stellar mass function, the mass–metallicity relation and the implied baryonic mass function

A comparison between published field galaxy stellar mass functions (GSMFs) shows that the cosmic stellar mass density is in the range 4–8 per cent of the baryon density (assuming  Ω_b= 0.045  ). There remain significant sources of uncertainty for the dust correction and underlying stellar mass-to-light ratio even assuming a reasonable universal stellar initial mass function. We determine the   z < 0.05  GSMF using the New York University Value-Added Galaxy Catalog sample of 49 968 galaxies derived from the Sloan Digital Sky Survey and various estimates of stellar mass. The GSMF shows clear evidence for a low-mass upturn and is fitted with a double Schechter function that has  α₂≃−1.6  . At masses below  ∼10^8.5 M_⊙  , the GSMF may be significantly incomplete because of missing low-surface-brightness galaxies. One interpretation of the stellar mass–metallicity relation is that it is primarily caused by a lower fraction of available baryons converted to stars in low-mass galaxies. Using this principle, we determine a simple relationship between baryonic mass and stellar mass and present an 'implied baryonic mass function'. This function has a faint-end slope,  α₂≃−1.9  . Thus, we find evidence that the slope of the low-mass end of the galaxy mass function could plausibly be as steep as the halo mass function. We illustrate the relationship between halo baryonic mass function → galaxy baryonic mass function → GSMF. This demonstrates the requirement for peak galaxy formation efficiency at baryonic masses  ∼10¹¹ M_⊙  corresponding to a minimum in feedback effects. The baryonic-infall efficiency may have levelled off at lower masses.     

The effect of supernova heating on cluster properties and constraints on galaxy formation models

We report the detection of a 5.8 Å– 10⁴ s periodicity in the 0.5–10 keV X-ray light curve of the Seyfert galaxy IRAS 18325–5926, obtained from a 5-d ASCA observation. Nearly nine cycles of the periodic variation are seen; it shows no strong energy dependence and has an amplitude of about 15 per cent. Unlike most other well-studied Seyfert galaxies, there is no evidence for strong power-law red noise in the X-ray power spectrum of IRAS 18325–5926. Scaling from the QPOs found in Galactic black hole candidates suggests that the mass of the black hole in IRAS 18325–5926 is ∼ 6 Å– 10⁶–4 Å– 10⁷ M_⊙.     

The tidal disruption rate in dense galactic cusps containing a supermassive binary black hole

We consider the problem of tidal disruption of stars in the centre of a galaxy containing a supermassive binary black hole with unequal masses. We assume that over the separation distance between the black holes, the gravitational potential is dominated by the more massive primary black hole. Also, we assume that the number density of stars is concentric with the primary black hole and has a power-law cusp. We show that the bulk of stars with a small angular-momentum component normal to the black hole binary orbit can reach a small value of total angular momentum through secular evolution in the gravitational field of the binary, and hence they can be tidally disrupted by the larger black hole. This effect is analogous to the so-called Kozai effect well known in celestial mechanics. We develop an analytical theory for the secular evolution of the stellar orbits and calculate the rate of tidal disruption. We compare our analytical theory with a simple numerical model and find very good agreement.
Our results show that for a primary black hole mass of  ∼10⁶–10⁷ M_⊙  , the black hole mass-ratio   q > 10⁻²  , cusp size ∼1 pc, the tidal disruption rate can be as large as  ∼10⁻²–1 M_⊙ yr⁻¹  . This is at least 10²–10⁴ times larger than estimated for the case of a single supermassive black hole. The duration of the phase of enhanced tidal disruption is determined by the dynamical-friction time-scale, and it is rather short: ∼10⁵ yr. The dependence of the tidal disruption rate on the mass ratio, and on the size of the cusp, is also discussed.     

Two new large-separation gravitational lenses from SDSS

We present discovery images, together with follow-up imaging and spectroscopy, of two large-separation gravitational lenses found by our survey for wide arcs [the CAmbridge Sloan Survey Of Wide ARcs in the skY (CASSOWARY)]. The survey exploits the multicolour photometry of the Sloan Digital Sky Survey to find multiple blue components around red galaxies. CASSOWARY 2 (or 'the Cheshire Cat') is composed of two massive early-type galaxies at   z = 0.426  and 0.432, respectively, lensing two background sources, the first a star-forming galaxy at   z = 0.97  and the second a high -redshift galaxy  ( z > 1.4)  . There are at least three images of the former source and probably four or more of the latter, arranged in two giant arcs. The mass enclosed within the larger arc of radius ∼11 arcsec is  ∼33 × 10¹² M_⊙  . CASSOWARY 3 comprises an arc of three bright images of a   z = 0.725  source, lensed by a foreground elliptical at   z = 0.274  . The radius of the arc is ∼4 arcsec and the enclosed mass is  ∼2.5 × 10¹² M_⊙  . Together with earlier discoveries like the Cosmic Horseshoe and the 8 o'clock Arc, these new systems, with separations intermediate between the arcsecond-separation lenses of typical strong galaxy lensing and arcminute-separation cluster lenses, probe the very high end of the galaxy mass function.     

A possible energy mechanism for cosmological γ-ray bursts

We suggest that an extreme Kerr black hole with a mass ∼10⁶ M_⊙, a dimensionless angular momentum     and a marginally stable orbital radius     located in a normal galaxy may produce a γ -ray burst (GRB) by capturing and disrupting a star. During the capture period, a transient accretion disc is formed and a strong transient magnetic field ∼     lasting for     may be produced at the inner boundary of the accretion disc. A large amount of rotational energy of the black hole is extracted and released in an ultrarelativistic jet with a bulk Lorentz factor Γ larger than 10³ via the Blandford–Znajek process. The relativistic jet energy can be converted into γ -radiation via an internal shock mechanism. The GRB duration should be the same as the lifetime of the strong transient magnetic field. The maximum number of sub-bursts is estimated to be     because the disc material is likely to break into pieces with a size about the thickness of the disc h at the cusp     The shortest risetime of the burst estimated from this model is ∼     The model GRB density rate is also estimated.     

Lens models with density cusps

Lens models appropriate for representing cusped galaxies and clusters are developed. The analogue of the odd-number theorem for cusped density distributions is given. Density cusps are classified into strong, isothermal or weak, according to their lensing properties. Strong cusps cause multiple imaging for any source position, whereas isothermal and weak cusps give rise to only one image for distant sources. Isothermal cusps always possess a pseudo-caustic. When the source crosses the pseudo-caustic, the number of images changes by unity.   Two families of cusped galaxy and cluster models are examined in detail. The double power-law family has an inner cusp, followed by a transition region and an outer envelope. One member of this family — the isothermal double power-law model — possesses an exceedingly scarce property, namely the lens equation is exactly solvable for any source position. This means that the magnifications, the time delay and the lensing cross-sections are all readily available. The model has a three-dimensional density that is cusped like r ⁻² at small radii and falls off like r ⁻⁴ asymptotically. Thus, it provides a very useful representation of the lensing properties of a galaxy or cluster of finite total mass with a flat rotation curve. The second set of models studied is the single power-law family. These are single density cusps of infinite extent. The properties of the critical curves and caustics and the behaviour of the lenses in the presence of external shear are all discussed in some detail.     

X-ray variability of the Seyfert 1 Markarian 335: power spectrum and time lags

To investigate further the comparison between active galactic nuclei (AGN) and black hole X-ray binaries, we have studied the main X-ray variability properties of the Seyfert 1 galaxy Mrk 335. We put particular emphasis on the X-ray time lags, which is a potentially important diagnostic of physical models. From a 100 ks observation by XMM–Newton , we show that the power spectrum of this source is well fitted by a bending power-law model, and the bend time-scale T _b is precisely at the value predicted by the T _b versus Hβ linewidth relation of McHardy et al. Variations in different energy bands show time-scale-dependent time lags, where higher energy bands lag lower ones. The lag, τ, varies as a function of the Fourier frequency, f , of the variability component in the light curves as  τ∝ f ⁻¹  at low frequencies, but there is a sharp cut-off in the lags at a frequency close to the bend frequency in the power spectrum. Similar behaviour is seen in black hole X-ray binary systems. The length of the time lags increases continuously with energy separation, in an almost loglinear relation. We show that the lag spectra can be produced by fluctuations propagating through the accretion flow as long as the energy spectrum of the X-ray emitting region hardens towards the centre.     

Supermassive Black Holes in Early-Type Galaxies: Relationship with Radio Emission and Constraints on the Black Hole Mass Function

Using recently published estimates — based on high spatial resolution spectroscopy — of the mass M _BH of nuclear black holes for a sample of nearby galaxies, we explore the dependence of galaxy nucleus emissivity at various wavelengths on M _BH. We confirm an almost linear scaling of the black hole mass with the baryonic mass of the host spheroidal galaxy. A remarkably tight relationship is also found with both nuclear and total radio centimetric flux, with a very steep dependence of the radio flux on M _BH ( P  ∝  M ^2.5_BH). The high-frequency radio power is thus a very good tracer of a supermassive black hole, and a good estimator of its mass. This, together with the lack of significant correlations with the low-energy X-ray and far-IR flux, supports the view that advection-dominated accretion is ruling the energy output in the low accretion rate regime. Using the tight dependence of total radio power on M _BH and the rich statistics of radio emission of galaxies, we derive an estimate of the mass function of remnants in the nearby Universe. This is compared with current models of quasar and active galactic nucleus (AGN) activity and of the origin of the hard X-ray background (HXRB). As for the former, continuous long-lived AGN activity is excluded by the present data with high significance, whereas the assumption of a short-lived, possibly recurrent, activity pattern gives remarkable agreement. The presently estimated black hole mass function also implies that the HXRB has been produced by a numerous population (∼ 10⁻² Mpc⁻³) of moderately massive ( M _BH ∼ 10⁷ M⊙) black holes.     

Models of the Cosmic Horseshoe gravitational lens J1004+4112

We model the extremely massive and luminous lens galaxy in the Cosmic Horseshoe Einstein ring system J1004+4112, recently discovered in the Sloan Digital Sky Survey. We use the semilinear method of Warren & Dye, which pixelizes the source surface brightness distribution, to invert the Einstein ring for sets of parametrized lens models. Here, the method is refined by exploiting Bayesian inference to optimise adaptive pixelization of the source plane and to choose between three differently parametrized models: a singular isothermal ellipsoid, a power-law model and a Navarro, Frenk & White (NFW) profile. The most probable lens model is the power law with a volume mass density  ρ∝ r ^{−1.96±0.02}  and an axis ratio of ∼0.8. The mass within the Einstein ring (i.e. within a cylinder with projected distance of ∼30 kpc from the centre of the lens galaxy) is  (5.02 ± 0.09) × 10¹² M _⊙  , and the mass-to-light ratio is ∼30. Even though the lens lies in a group of galaxies, the preferred value of the external shear is almost zero. This makes the Cosmic Horseshoe unique amongst large separation lenses, as almost all the deflection comes from a single, very massive galaxy with little boost from the environment.     

Galaxy density profiles and shapes – II. Selection biases in strong lensing surveys

Many current and future astronomical surveys will rely on samples of strong gravitational lens systems to draw conclusions about galaxy mass distributions. We use a new strong lensing pipeline (presented in Paper I of this series) to explore selection biases that may cause the population of strong lensing systems to differ from the general galaxy population. Our focus is on point-source lensing by early-type galaxies with two mass components (stellar and dark matter) that have a variety of density profiles and shapes motivated by observational and theoretical studies of galaxy properties. We seek not only to quantify but also to understand the physics behind selection biases related to: galaxy mass, orientation and shape; dark matter profile parameters such as inner slope and concentration; and adiabatic contraction. We study how all of these properties affect the lensing Einstein radius, total cross-section, quad/double ratio and image separation distribution, with a flexible treatment of magnification bias to mimic different survey strategies. We present our results for two families of density profiles: cusped and deprojected Sérsic models. While we use fixed lens and source redshifts for most of the analysis, we show that the results are applicable to other redshift combinations, and we also explore the physics of how our results change for very different redshifts. We find significant (factors of several) selection biases with mass; orientation, for a given galaxy shape at fixed mass; cusped dark matter profile inner slope and concentration; concentration of the stellar and dark matter deprojected Sérsic models. Interestingly, the intrinsic shape of a galaxy does not strongly influence its lensing cross-section when we average over viewing angles. Our results are an important first step towards understanding how strong lens systems relate to the general galaxy population.     

Measuring the black hole masses of high-redshift quasars

A new technique is presented for determining the black hole masses of high-redshift quasars from optical spectroscopy. The new method utilizes the full-width at half-maximum (FWHM) of the low-ionization Mg  ii emission line and the correlation between the broad-line region (BLR) radius and the continuum luminosity at 3000 Å. Using archival ultraviolet (UV) spectra it is found that the correlation between BLR radius and 3000-Å luminosity is tighter than the established correlation with 5100-Å luminosity. Furthermore, it is found that the correlation between BLR radius and 3000-Å continuum luminosity is consistent with a relation of the form   R _BLR∝λ L ^1/2_λ  , as expected for a constant ionization parameter. Using a sample of objects with broad-line radii determined from reverberation mapping it is shown that the FWHM of Mg  ii and Hβ are consistent with following an exact one-to-one relation, as expected if both Hβ and Mg  ii are emitted at the same radius from the central ionizing source. The resulting virial black hole mass estimator based on rest-frame UV observables is shown to reproduce black hole mass measurements based on reverberation mapping to within a factor of 2.5 (1σ). Finally, the new UV black hole mass estimator is shown to produce identical results to the established optical (Hβ) estimator when applied to 128 intermediate-redshift  (0.3 < z < 0.9)  quasars drawn from the Large Bright Quasar Survey and the radio-selected Molonglo quasar sample. We therefore conclude that the new UV virial black hole mass estimator can be reliably used to estimate the black hole masses of quasars from   z ∼ 0.25  through to the peak epoch of quasar activity at   z ∼ 2.5  via optical spectroscopy alone.     

Black hole demographics from the relation

We analyse a sample of 32 galaxies for which a dynamical estimate of the mass of the hot stellar component, M _bulge, is available. For each of these galaxies, we calculate the mass of the central black hole, M _•, using the tight empirical correlation between M _• and bulge stellar velocity dispersion. The frequency function     is reasonably well described as a Gaussian with     and standard deviation ∼0.45; the implied mean ratio of black hole mass to bulge mass is a factor of ∼5 smaller than generally quoted in the literature. We present marginal evidence for a lower, average black hole mass fraction in more massive galaxies. The total mass density in black holes in the local Universe is estimated to be ∼     consistent with that inferred from high-redshift     active galactic nuclei.     

Mass and gas profiles in A1689: joint X-ray and lensing analysis

We carry out a comprehensive joint analysis of high-quality HST /ACS and Chandra measurements of A1689, from which we derive mass, temperature, X-ray emission and abundance profiles. The X-ray emission is smooth and symmetric, and the lensing mass is centrally concentrated indicating a relaxed cluster. Assuming hydrostatic equilibrium we deduce a 3D mass profile that agrees simultaneously with both the lensing and X-ray measurements. However, the projected temperature profile predicted with this 3D mass profile exceeds the observed temperature by ∼30 per cent at all radii, a level of discrepancy comparable to the level found for other relaxed clusters. This result may support recent suggestions from hydrodynamical simulations that denser, more X-ray luminous small-scale structure can bias observed temperature measurements downward at about the same (∼30 per cent) level. We determine the gas entropy at  0.1 r _vir  (where r _vir is the virial radius) to be ∼800 keV cm², as expected for a high-temperature cluster, but its profile at  >0.1 r _vir  has a power-law form with index ∼0.8, considerably shallower than the ∼1.1 index advocated by theoretical studies and simulations. Moreover, if a constant entropy 'floor' exists at all, then it is within a small region in the inner core,   r < 0.02 r _vir  , in accord with previous theoretical studies of massive clusters.     

J0041+3224: a new double–double radio galaxy

We report the discovery of a double–double radio galaxy (DDRG), J0041+3224, with the Giant Metrewave Radio Telescope (GMRT) and subsequent high-frequency observations with the Very Large Array (VLA). The inner and outer doubles are aligned within ∼4° and are reasonably collinear with the parent optical galaxy. The outer double has a steeper radio spectrum compared to the inner one. Using an estimated redshift of 0.45, the projected linear sizes of the outer and inner doubles are 969 and 171 kpc, respectively. The time-scale of interruption of jet activity has been estimated to be ∼20 Myr, similar to other known DDRGs. We have compiled a sample of known DDRGs, and have re-examined the inverse correlation between the ratio of the luminosities of the outer to the inner double and the size of the inner double, l _in. Unlike the other DDRGs with   l _in≳ 50 kpc  , the inner double of J0041+3224 is marginally more luminous than the outer one. The two DDRGs with   l _in≲  few kpc have a more luminous inner double than the outer one, possibly due to a higher efficiency of conversion of beam energy as the jets propagate through the dense interstellar medium. We have examined the symmetry parameters and found that the inner doubles appear to be more asymmetric in both its armlength and its flux density ratios compared to the outer doubles, although they appear marginally more collinear with the core than the outer double. We discuss briefly the possible implications of these trends.     

Chandra measurements of the distribution of mass in the luminous lensing cluster Abell 2390

We present spatially resolved X-ray spectroscopy of the luminous lensing cluster Abell 2390, using observations made with the Chandra observatory. The temperature of the X-ray gas rises with increasing radius within the central ∼ 200 kpc of the cluster, and then remains approximately isothermal, with kT =11.5_−1.6^+1.5 keV , out to the limits of the observations at r ∼1.0 Mpc . The total mass profile determined from the Chandra data has a form in good agreement with the predictions from numerical simulations. Using the parametrization of Navarro, Frenk and White, we measure a scale radius r _s∼0.8 Mpc and a concentration parameter c ∼3 . The best-fitting X-ray mass model is in good agreement with independent gravitational lensing results and optical measurements of the galaxy velocity dispersion in the cluster. The X-ray gas to total mass ratio rises with increasing radius with f _gas∼21 per cent at r =0.9 Mpc . The azimuthally averaged 0.3–7.0 keV surface brightness profile exhibits a small core radius and a clear 'break' at r ∼500 kpc , where the slope changes from S _X ∼^∝  r ^−1.5 to S _X ∼^∝  r ^−3.6 . The data for the central region of the cluster indicate the presence of a cooling flow with a mass deposition rate of 200–300 M_⊙ yr⁻¹ and an effective age of 2–3 Gyr .