An 80-kpc Lyα halo around a high-redshift type-2 quasi-stellar object

We announce the discovery of an extended emission-line region associated with a high-redshift type-2 quasi-stellar object (QSO). The halo, which was discovered in our new wide-field narrow-band survey, resides at   z = 2.85  in the Spitzer First Look Survey region and is extended over ∼80 kpc. Deep very long baseline interferometry (VLBI) observations imply that approximately 50 per cent of the radio emission is extended on scales >200 pc. The inferred active galactic nuclei (AGN) luminosity is sufficient to ionize the extended halo, and the optical emission is consistent with being triggered coevally with the radio source. The Lyα halo is as luminous as those found around high-redshift radio galaxies; however, the active nucleus is several orders of magnitude less luminous at radio wavelengths than those Fanarof–Riley type II (FRIIs) more commonly associated with extended emission-line regions. AMS05 appears to be a high-redshift analogue to the radio-quiet quasar E1821+643 which is core dominated, but which also exhibits extended Fanarof–Riley type I (FRI)-like structure and contains an optically powerful AGN. We also find evidence for more quiescent kinematics in the Lyα emission line in the outer regions of the halo, reminiscent of the haloes around the more powerful FRIIs. The optical to mid-infrared spectral energy distribution is well described by a combination of an obscured QSO  ( L _bol∼ 3.4 ± 0.2 × 10¹³ L_⊙)  and a 1.4 Gyr old simple stellar population with mass  ∼3.9 ± 0.3 × 10¹¹ M_⊙  .     

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.     

Downsizing of supermassive black holes from the SDSS quasar survey

Starting from the ∼50 000 quasars of the Sloan Digital Sky Survey for which Mg  ii line width and 3000 Å monochromatic flux are available, we aim to study the dependence of the mass of active black holes on redshift. We focus on the observed distribution in the full width at half-maximum–nuclear luminosity plane, which can be reproduced at all redshifts assuming a limiting M _BH, a maximum Eddington ratio and a minimum luminosity (due to the survey flux limit). We study the z -dependence of the best-fitting parameters of assumed distributions at increasing redshift and find that the maximum mass of the quasar population evolves as  log ( M _BH(max)/M_⊙) ∼ 0.3 z + 9  , while the maximum Eddington ratio (∼0.45) is practically independent of cosmic time. These results are unaffected by the Malmquist bias.     

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.     

A constant dark matter halo surface density in galaxies

We confirm and extend the recent finding that the central surface density  μ_0D≡ r ₀ρ₀  of galaxy dark matter haloes, where r ₀ and  ρ₀  are the halo core radius and central density, is nearly constant and independent of galaxy luminosity. Based on the co-added rotation curves (RCs) of ∼1000 spiral galaxies, the mass models of individual dwarf irregular and spiral galaxies of late and early types with high-quality RCs, and the galaxy–galaxy weak-lensing signals from a sample of spiral and elliptical galaxies, we find that  log μ_0D= 2.15 ± 0.2  in units of  log(M_⊙ pc⁻²)  . We also show that the observed kinematics of Local Group dwarf spheroidal galaxies are consistent with this value. Our results are obtained for galactic systems spanning over 14 mag, belonging to different Hubble types and whose mass profiles have been determined by several independent methods. In the same objects, the approximate constancy of  μ_0D  is in sharp contrast to the systematical variations, by several orders of magnitude, of galaxy properties, including  ρ₀  and central stellar surface density.     

Baryonic conversion tree: the global assembly of stars and dark matter in galaxies from the Sloan Digital Sky Survey

Using the spectroscopic sample of the Sloan Digital Sky Survey Data Release 1 (SDSS DR1), we measure how gas was transformed into stars as a function of time and stellar mass: the baryonic conversion tree (BCT). There is a clear correlation between early star formation activity and present-day stellar mass: the more massive galaxies have formed approximately 80 per cent of their stars at   z > 1  , while for the less massive ones the value is only approximately 20 per cent. By comparing the BCT with the dark matter merger tree, we find indications that star formation efficiency at   z > 1  had to be approximately a factor of two higher than today (∼10 per cent) in galaxies with present-day stellar mass larger than  2 × 10¹¹ M_⊙  , if this early star formation occurred in the main progenitor. Therefore, the λ cold dark matter (LCDM) paradigm can accommodate a large number of red objects. On the other hand, in galaxies with present-day stellar mass less than  10¹¹ M_⊙  , efficient star formation seems to have been triggered at   z ∼ 0.2  . We show that there is a characteristic mass  ( M _*∼ 10¹⁰ M_⊙)  for feedback efficiency (or lack of star formation). For galaxies with masses lower than this, feedback (or star formation suppression) is very efficient while for higher masses it is not. The BCT, determined here for the first time, should be an important observable with which to confront theoretical models of galaxy formation.     

Semi-analytic simulations of galactic winds: volume filling factor, ejection of metals and parameter study

We present a semi-analytic treatment of galactic winds within high-resolution, large-scale cosmological N -body simulations of a Λ cold dark matter (ΛCDM) universe. The evolution of winds is investigated by following the expansion of supernova-driven superbubbles around the several hundred thousand galaxies that form in an approximately spherical region of space with diameter 52  h ⁻¹ Mpc and mean density close to the mean density of the universe. We focus our attention on the impact of winds on the diffuse intergalactic medium. Initial conditions for mass loss at the base of winds are taken from Shu, Mo & Mao. Results are presented for the volume filling factor and the mass fraction of the intergalactic medium (IGM) affected by winds, and their dependence on the model parameters is carefully investigated. The mass-loading efficiency of bubbles is a key factor to determine the evolution of winds and their global impact on the IGM: the higher the mass loading, the later the IGM is enriched with metals. Galaxies with 10⁹ < M _★ < 10¹⁰ M_⊙ are responsible for most of the metals ejected into the IGM at   z = 3  , while galaxies with   M _★ < 10⁹ M_⊙   give a non-negligible contribution only at higher redshifts, when larger galaxies have not yet assembled. We find a higher mean IGM metallicity than Lyα forest observations suggest, and we argue that the discrepancy may be explained by the high temperatures of a large fraction of the metals in winds, which may not leave detectable imprints in absorption in the Lyα forest.     

An iterative filter to reconstruct planetary transit signals in the presence of stellar variability

We report the identification, from a photometric, astrometric and spectroscopic study, of a massive white dwarf member of the nearby, approximately solar metallicity, Coma Berenices open star cluster (Melotte 111). We find the optical to near-infrared energy distribution of WD 1216+260 to be entirely consistent with that of an isolated DA and determine the effective temperature and surface gravity of this object to be   T _eff= 15 739⁺¹⁹⁷₋₁₉₆ K  and  log  g = 8.46^+0.03_−0.02  . We set tight limits on the mass of a putative cool companion,   M ≳ 0.036 M_⊙  (spatially unresolved) and   M ≳ 0.034 M_⊙  (spatially resolved and   a ≲ 2500 au  ). Based on the predictions of CO core, thick H layer evolutionary models we determine the mass and cooling time of WD 1216+260 to be   M _WD= 0.90 ± 0.04 M_⊙  and  τ_cool= 363⁺⁴⁶₋₄₁ Myr  , respectively. For an adopted cluster age of  τ= 500 ± 100 Myr  we infer the mass of its progenitor star to be   M _init= 4.77^+5.37_−0.97 M_⊙  . We briefly discuss this result in the context of the form of the stellar initial mass–final mass relation.     

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.     

Microlensing in dark matter haloes

Using eight dark matter haloes extracted from fully self-consistent cosmological N -body simulations, we perform microlensing experiments. A hypothetical observer is placed at a distance of 8.5 kpc from the centre of the halo measuring optical depths, event durations and event rates towards the direction of the Large Magellanic Cloud. We simulate 1600 microlensing experiments for each halo. Assuming that the whole halo consists of massive astronomical compact halo objects (MACHOs),   f = 1.0  , and a single MACHO mass is   m _M= 1.0 M_⊙  , the simulations yield mean values of  τ= 4.7^+5.0_−2.2× 10⁻⁷  and  Γ= 1.6^+1.3_−0.6× 10⁻⁶  events star⁻¹ yr⁻¹. We find that triaxiality and substructure can have major effects on the measured values so that τ and Γ values of up to three times the mean can be found. If we fit our values of τ and Γ to the MACHO collaboration observations, we find   f = 0.23^+0.15_−0.13  and   m _M= 0.44^+0.24_−0.16  . Five out of the eight haloes under investigation produce f and m _M values mainly concentrated within these bounds.     

The distribution of ejected subhaloes and its implication for halo assembly bias

Using a high-resolution cosmological N -body simulation, we identify the ejected population of subhaloes, which are haloes at redshift   z = 0  but were once contained in more massive 'host' haloes at high redshifts. The fraction of the ejected subhaloes in the total halo population of the same mass ranges from 9 to 4 per cent for halo masses from  ∼10¹¹  to  ∼10¹²  h ⁻¹ M_⊙  . Most of the ejected subhaloes are distributed within four times the virial radius of their hosts. These ejected subhaloes have distinct velocity distribution around their hosts in comparison to normal haloes. The number of subhaloes ejected from a host of given mass increases with the assembly redshift of the host. Ejected subhaloes in general reside in high-density regions, and have a much higher bias parameter than normal haloes of the same mass. They also have earlier assembly times, so that they contribute to the assembly bias of dark matter haloes seen in cosmological simulations. However, the assembly bias is not dominated by the ejected population, indicating that large-scale environmental effects on normal haloes are the main source for the assembly bias.     

New Praesepe white dwarfs and the initial mass–final mass relation

We report the spectroscopic confirmation of four further white dwarf members of Praesepe. This brings the total number of confirmed white dwarf members to 11, making this the second largest collection of these objects in an open cluster identified to date. This number is consistent with the high-mass end of the initial mass function of Praesepe being Salpeter in form. Furthermore, it suggests that the bulk of Praesepe white dwarfs did not gain a substantial recoil kick velocity from possible asymmetries in their loss of mass during the asymptotic giant branch phase of evolution. By comparing our estimates of the effective temperatures and the surface gravities of WD0833+194, WD0840+190, WD0840+205 and WD0843+184 to modern theoretical evolutionary tracks, we have derived their masses to be in the range  0.72–0.76 M_⊙  and their cooling ages ∼300 Myr. For an assumed cluster age of 625 ± 50 Myr, the inferred progenitor masses are between 3.3 and  3.5 M_⊙  . Examining these new data in the context of the initial mass–final mass relation, we find that it can be adequately represented by a linear function  ( a ₀= 0.289 ± 0.051,  a ₁= 0.133 ± 0.015)  over the initial mass range 2.7–6  M_⊙  . Assuming an extrapolation of this relation to larger initial masses is valid and adopting a maximum white dwarf mass of  1.3 M_⊙  , our results support a minimum mass for core-collapse supernovae progenitors in the range  ∼6.8–8.6 M_⊙  .     

Quantifying the coexistence of massive black holes and dense nuclear star clusters

In large spheroidal stellar systems, such as elliptical galaxies, one invariably finds a  10⁶–10⁹ M_⊙  supermassive black hole at their centre. In contrast, within dwarf elliptical galaxies one predominantly observes a  10⁵–10⁷ M_⊙  nuclear star cluster. To date, few galaxies have been found with both types of nuclei coexisting and even less have had the masses determined for both central components. Here, we identify one dozen galaxies housing nuclear star clusters and supermassive black holes whose masses have been measured. This doubles the known number of such hermaphrodite nuclei – which are expected to be fruitful sources of gravitational radiation. Over the host spheroid (stellar) mass range  10⁸–10¹¹ M_⊙  , we find that a galaxy's nucleus-to-spheroid (baryon) mass ratio is not a constant value but decreases from a few per cent to ∼0.3 per cent such that  log[( M _BH+ M _NC)/ M _sph]=−(0.39 ± 0.07) log[ M _sph/10¹⁰ M_⊙]− (2.18 ± 0.07)  . Once dry merging commences and the nuclear star clusters disappear, this ratio is expected to become a constant value.
As a byproduct of our investigation, we have found that the projected flux from resolved nuclear star clusters is well approximated with Sérsic functions having a range of indices from ∼0.5 to ∼3, the latter index describing the Milky Way's nuclear star cluster.     

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.     

Dark matter subhaloes in numerical simulations

We use cosmological Λ cold dark matter (CDM) numerical simulations to model the evolution of the substructure population in 16 dark matter haloes with resolutions of up to seven million particles within the virial radius. The combined substructure circular velocity distribution function (VDF) for hosts of 10¹¹ to  10¹⁴ M_⊙  at redshifts from zero to two or higher has a self-similar shape, is independent of host halo mass and redshift, and follows the relation  d n /d v = (1/8)( v _cmax/ v _cmax,host)⁻⁴  . Halo to halo variance in the VDF is a factor of roughly 2 to 4. At high redshifts, we find preliminary evidence for fewer large substructure haloes (subhaloes). Specific angular momenta are significantly lower for subhaloes nearer the host halo centre where tidal stripping is more effective. The radial distribution of subhaloes is marginally consistent with the mass profile for   r ≳ 0.3 r _vir  , where the possibility of artificial numerical disruption of subhaloes can be most reliably excluded by our convergence study, although a subhalo distribution that is shallower than the mass profile is favoured. Subhalo masses but not circular velocities decrease towards the host centre. Subhalo velocity dispersions hint at a positive velocity bias at small radii. There is a weak bias towards more circular orbits at lower redshift, especially at small radii. We additionally model a cluster in several power-law cosmologies of   P ∝ kⁿ   , and demonstrate that a steeper spectral index, n , results in significantly less substructure.     

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.     

Dry mergers: a crucial test for galaxy formation

We investigate the role that dry mergers play in the build-up of massive galaxies within the cold dark matter paradigm. Implementing an empirical shut-off mass scale for star formation, we find a nearly constant dry merger rate of  ∼6 × 10⁻⁵ Mpc⁻³ Gyr⁻¹  at   z ≤ 1  and a steep decline at larger z . Less than half of these mergers are between two galaxies that are morphologically classified as early-types, and the other half is mostly between an early- and late-type galaxy. Latter are prime candidates for the origin of tidal features around red elliptical galaxies. The introduction of a transition mass scale for star formation has a strong impact on the evolution of galaxies, allowing them to grow above a characteristic mass scale of   M _{*, c} ∼ 6.3 × 10¹⁰ M_⊙  by mergers only. As a consequence of this transition, we find that around   M _{*, c}   , the fraction of 1:1 mergers is enhanced with respect to unequal mass major mergers. This suggests that it is possible to detect the existence of a transition mass scale by measuring the relative contribution of equal mass mergers to unequal mass mergers as a function of galaxy mass. The evolution of the high-mass end of the luminosity function is mainly driven by dry mergers at low z . We however find that only 10–20 per cent of galaxies more massive than   M _{*, c}   experience dry major mergers within their last Gyr at any given redshift   z ≤ 1  .     

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.     

H i and dark matter in the windy starburst dwarf galaxy NGC 1705

We present 21-cm H  i line observations of the blue compact dwarf galaxy NGC 1705. Previous optical observations show a strong outflow powered by an ongoing starburst dominating the H  ii morphology and kinematics. In contrast, most of the H  i lies in a rotating disc. An extraplanar H  i spur accounts for ∼8 per cent of the total H  i mass, and is possibly associated with the H  ii outflow. The inferred mass loss rate out of the core of the galaxy is significant, ∼0.2 − 2 M_⊙ yr⁻¹, but does not dominate the H  i dynamics. Mass model fits to the rotation curve show that the dark matter (DM) halo is dominant at nearly all radii and has a central density ρ₀ ≈ 0.1 M_⊙ pc⁻³: ten times higher than typically found in dwarf irregular galaxies, but similar to the only other mass-modelled blue compact dwarf, NGC 2915. This large difference strongly indicates that there is little evolution between dwarf irregular and blue compact dwarf types. Instead, dominant DM haloes may regulate the morphology of dwarf galaxies by setting the critical surface density for disc star formation. Neither our data nor catalogue searches reveal any likely external trigger to the starburst in NGC 1705.     

Ly+ and ultraviolet emission from high-redshift gamma-ray burst hosts: to what extent do gamma-ray bursts trace star formation?

We report the result of a search for Lyα emission from the host galaxies of the gamma-ray bursts  (GRBs) 030226 ( z = 1.986), 021004 ( z = 2.335)  and  020124 ( z = 3.198)  . We find that the host galaxy of GRB 021004 is an extended (around 8 kpc) strong Lyα emitter with a rest-frame equivalent width (EW) of 68⁺¹²₋₁₁Å, and a star formation rate of  10.6 ± 2.0 M_⊙ yr⁻¹  . We do not detect the hosts of GRB 030226 and GRB 020124, but the upper limits on their Lyα fluxes do not rule out large rest-frame EWs. In the fields of GRB 021004 and GRB 030226 we find seven and five other galaxies, respectively, with excess emission in the narrow-band filter. These galaxies are candidate Lyα-emitting galaxies in the environment of the host galaxies. We have also compiled a list of all   z ≳ 2  GRB hosts, and demonstrate that a scenario where they trace star formation in an unbiased way is compatible with current observational constraints. Fitting the   z = 3  luminosity function (LF) under this assumption results in a characteristic luminosity of   R *= 24.6  and a faint-end slope of  α=−1.55  , consistent with the LF measured for Lyman-break galaxies.