H 0 from an orientation-unbiased sample of Sunyaev–Zel'dovich and X-ray clusters

We have observed the Sunyaev–Zel'dovich (SZ) effect in a sample of five moderate-redshift clusters with the Ryle Telescope, and used them in conjunction with X-ray imaging and spectral data from ROSAT and ASCA to measure the Hubble constant. This sample was chosen with a strict X-ray flux limit using both the Bright Cluster Sample and the Northern ROSAT All-Sky Survey (RASS) cluster catalogues to be well above the surface brightness limit of the RASS, and hence to be unbiased with respect to the orientation of the cluster. This controls a major potential systematic effect in the SZ/X-ray method of measuring H ₀. Taking the weighted geometric mean of the results and including the main sources of error, namely the noise in the SZ measurement, the uncertainty in the X-ray temperatures and the unknown ellipticity and substructure of the clusters, we find   H ₀= 59⁺¹⁰₋₉ (random)⁺⁸₋₇(systematic) km s⁻¹ Mpc⁻¹  assuming a standard cold dark matter model with  Ω_M= 1.0, Ω_Λ= 0.0  or   H ₀= 66⁺¹¹₋₁₀ ⁺⁹₋₈ km  s⁻¹ Mpc⁻¹  if  Ω_M= 0.3, Ω_Λ= 0.7  .     

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.     

RR Pic (1925): a Chandra X-ray view

We present the Chandra ACIS-S3 data of the old classical nova RR Pic (1925). The source has a count rate of 0.067 ± 0.002 count s⁻¹ in the 0.3–5.0 keV energy range. We detect the orbital period of the underlying binary system in the X-ray wavelengths. We also find that the neutral hydrogen column density differs for orbital minimum and orbital maximum spectra with values  0.25^+0.23_−0.18× 10²²  and  0.64^+0.13_−0.14× 10²² cm⁻²  at 3σ confidence level. The X-ray spectrum of RR Pic can be represented by a composite model of bremsstrahlung with a photoelectric absorption, two absorption lines centered around 1.1–1.4 keV and five Gaussian lines centered at emission lines around 0.3–1.1 keV corresponding to various transitions of S, N, O, C, Ne and Fe. The bremsstrahlung temperature derived from the fits ranges from 0.99 to 1.60 keV and the unabsorbed X-ray flux is found to be  2.5^+0.4_−1.2× 10⁻¹³ erg  cm⁻² s⁻¹  in the 0.3–5.0 keV range with a luminosity of 1.1 ± 0.2  10³¹ erg s⁻¹  at 600 pc. We also detect excess emission in the spectrum possibly originating from the reverse shock in the ejecta. A fit with a cooling flow plasma emission model shows enhanced abundances of He, C, N, O and Ne in the X-ray emitting region indicating existence of diffusive mixing.     

ROSAT ASCA observations of X-ray luminous starburst galaxies: NGC 3310 and 3690

We present ROSAT [High Resolution Imager (HRI) and Position Sensitive Proportional Counter (PSPC)] and ASCA observations of the two luminous ( L _x ∼ 10⁴¹⁻⁴² erg s⁻¹) star-forming galaxies NGC 3310 and 3690. The HRI shows clearly that the sources are extended with the X-ray emission in NGC 3690 coming from at least three regions. The combined 0.1–10 keV spectrum of NGC 3310 can be described by two components, a Raymond–Smith plasma with temperature kT  = 0.81^+0.09_−0.12 keV and a hard power law, Γ = 1.44^−0.20_−0.11 (or alternatively a harder Raymond–Smith plasma with kT  ∼ 15 keV), while there is no substantial excess absorption above the Galactic column value. The soft component emission is probably a super wind while the nature of the hard emission is more uncertain with the likely origins being X-ray binaries, inverse Compton scattering of infrared photons, an active galactic nucleus or a very hot gas component (∼10⁸ K). The spectrum of NGC 3690 is similar, with kT  = 0.83^+0.02_−0.04 keV and Γ = 1.56^+0.11_−0.11. We also employ more complicated models such as a multi-temperature thermal plasma, a non-equilibrium ionization code or the addition of a third softer component, which improve the fit but not at a statistically significant level (2σ). These results are similar to recent results on the archetypal star-forming galaxies M82 and NGC 253.     

Deep galaxy counts, extragalactic background light and the stellar baryon budget

We assess the constraints imposed by the observed extragalactic background light (EBL) on the cosmic history of star formation and the stellar-mass density today. The logarithmic slope of the galaxy number–magnitude relation from the Southern Hubble Deep Field imaging survey is flatter than 0.4 in all seven UBVIJHK optical bandpasses, i.e. the light from resolved galaxies has converged from the UV to the near-IR. We find a lower limit to the surface brightness of the optical extragalactic sky of about 15 nW m⁻² sr⁻¹, comparable to the intensity of the far-IR background from COBE data. Assuming a Salpeter initial mass function with a lower cut-off consistent with observations of M subdwarf disc stars, we set a lower limit of Ω_g+s h ²>0.0013  I ₅₀ to the visible (processed gas + stars) mass density required to generate an EBL at a level of 50  I ₅₀ nW m⁻² sr⁻¹; our 'best-guess' value is Ω_g+s h ²≈0.0031  I ₅₀. Motivated by the recent microlensing results of the MACHO collaboration, we consider the possibility that massive dark haloes around spiral galaxies are composed of faint white dwarfs, and show that only a small fraction (≲5 per cent) of the nucleosynthetic baryons can be locked in the remnants of intermediate-mass stars forming at z _F≲5, as the bright early phases of such haloes would otherwise overproduce the observed EBL.     

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.     

Bulk flow and shear moments of the SFI++ survey

We find the nine bulk flow and shear moments from the SFI++ survey, as well as for subsamples of group and field galaxies. We constrain the velocity power spectrum shape parameter Γ in linear theory using these moments. A likelihood function for Γ was found after marginalizing over the power spectrum amplitude  σ₈Ω^0.6_m  using constraints obtained from comparisons between redshift surveys and peculiar velocity data. We have estimated the velocity noise  σ_*  from the data since without it our results may be biased. We also performed a statistical analysis of the difference between the field and group catalogues and found that the results from each reflect the same underlying large-scale flows. We found that we can constrain the power spectrum shape parameter to be  Γ= 0.15^+0.18_−0.08  for the groups catalogue and  Γ= 0.09^+0.04_−0.04  for the field galaxy catalogue in fair agreement with the value from Wilkinson Microwave Anisotropy Probe .     

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.     

Sakurai's Object: characterizing the near-infrared CO ejecta between 2003 and 2007

We present observations of Sakurai's Object obtained at 1–5 μm between 2003 and 2007. By fitting a radiative transfer model to an echelle spectrum of CO fundamental absorption features around  4.7 μm  , we determine the excitation conditions in the line-forming region. We find  ¹²C/¹³C = 3.5^+2.0_−1.5  , consistent with CO originating in ejecta processed by the very late thermal pulse, rather than in the pre-existing planetary nebula. We demonstrate the existence of  2.2 × 10⁻⁶≤ M _CO≤ 2.7 × 10⁻⁶ M_⊙  of CO ejecta outside the dust, forming a high-velocity wind of  500 ± 80 km s⁻¹  . We find evidence for significant weakening of the CO band and cooling of the dust around the central star between 2003 and 2005. The gas and dust temperatures are implausibly high for stellar radiation to be the sole contributor.     

Understanding the halo-mass and galaxy-mass cross-correlation functions

We use the Millennium Simulation (MS) to measure the cross-correlation between halo centres and mass (or equivalently the average density profiles of dark haloes) in a Lambda cold dark matter (ΛCDM) cosmology. We present results for radii in the range  10  h ⁻¹ kpc < r < 30  h ⁻¹ Mpc  and for halo masses in the range  4 × 10¹⁰ < M ₂₀₀ < 4 × 10¹⁴  h ⁻¹ M_⊙  . Both at   z = 0  and at   z = 0.76  these cross-correlations are surprisingly well fitted if the inner region is approximated by a density profile of NFW or Einasto form, the outer region by a biased version of the linear mass autocorrelation function, and the maximum of the two is adopted where they are comparable. We use a simulation of galaxy formation within the MS to explore how these results are reflected in cross-correlations between galaxies and mass. These are directly observable through galaxy–galaxy lensing. Here also we find that simple models can represent the simulation results remarkably well, typically to ≲10 per cent. Such models can be used to extend our results to other redshifts, to cosmologies with other parameters, and to other assumptions about how galaxies populate dark haloes. Our galaxy formation simulation already reproduces current galaxy–galaxy lensing data quite well. The characteristic features predicted in the galaxy–galaxy lensing signal should provide a strong test of the ΛCDM cosmology as well as a route to understanding how galaxies form within it.     

The distribution function of dark matter in massive haloes

We study the distribution function (DF) of dark matter particles in haloes of mass range  10¹⁴–10¹⁵ M_⊙  . In the numerical part of this work we measure the DF for a sample of relaxed haloes formed in the simulation of a standard Λ cold dark matter (ΛCDM) model. The DF is expressed as a function of energy E and the absolute value of the angular momentum L , a form suitable for comparison with theoretical models. By proper scaling we obtain the results that do not depend on the virial mass of the haloes. We demonstrate that the DF can be separated into energy and angular momentum components and propose a phenomenological model of the DF in the form     . This formulation involves three parameters describing the anisotropy profile in terms of its asymptotic values (β₀ and  β_∞  ) and the scale of transition between them ( L ₀). The energy part   f _E ( E )  is obtained via inversion of the integral for spatial density. We provide a straightforward numerical scheme for this procedure as well as a simple analytical approximation for a typical halo formed in the simulation. The DF model is extensively compared with the simulations: using the model parameters obtained from fitting the anisotropy profile, we recover the DF from the simulation as well as the profiles of the dispersion and kurtosis of radial and tangential velocities. Finally, we show that our DF model reproduces the power-law behaviour of phase-space density   Q =ρ( r )/σ³( r )  .     

Resolving cosmic structure formation with the Millennium-II Simulation

We present the Millennium-II Simulation (MS-II), a very large N -body simulation of dark matter evolution in the concordance Λ cold dark matter (ΛCDM) cosmology. The MS-II assumes the same cosmological parameters and uses the same particle number and output data structure as the original Millennium Simulation (MS), but was carried out in a periodic cube one-fifth the size  (100  h ⁻¹ Mpc)  with five times better spatial resolution (a Plummer equivalent softening of  1.0  h ⁻¹ kpc  ) and with 125 times better mass resolution (a particle mass of  6.9 × 10⁶  h ⁻¹ M_⊙  ). By comparing results at MS and MS-II resolution, we demonstrate excellent convergence in dark matter statistics such as the halo mass function, the subhalo abundance distribution, the mass dependence of halo formation times, the linear and non-linear autocorrelations and power spectra, and halo assembly bias. Together, the two simulations provide precise results for such statistics over an unprecedented range of scales, from haloes similar to those hosting Local Group dwarf spheroidal galaxies to haloes corresponding to the richest galaxy clusters. The 'Milky Way' haloes of the Aquarius Project were selected from a lower resolution version of the MS-II and were then resimulated at much higher resolution. As a result, they are present in the MS-II along with thousands of other similar mass haloes. A comparison of their assembly histories in the MS-II and in resimulations of 1000 times better resolution shows detailed agreement over a factor of 100 in mass growth. We publicly release halo catalogues and assembly trees for the MS-II in the same format within the same archive as those already released for the MS.     

Towards a precision cosmology from starburst galaxies at z > 2

This work investigates the use of a well-known empirical correlation between the velocity dispersion, metallicity and luminosity in Hβ of nearby H  ii galaxies to measure the distances to H  ii -like starburst galaxies at high redshifts. This correlation is applied to a sample of 15 starburst galaxies with redshifts between   z = 2.17  and   z = 3.39  to constrain  Ω_m  , using data available from the literature. A best-fitting value of  Ω_m= 0.21^+0.30_−0.12  in a Λ-dominated universe and of  Ω_m= 0.11^+0.37_−0.19  in an open universe is obtained. A detailed analysis of systematic errors, their causes and their effects on the values derived for the distance moduli and  Ω_m  is carried out. A discussion of how future work will improve constraints on  Ω_m  by reducing the errors is also presented.     

High-energy emission from the starburst galaxy NGC 253

Measurement sensitivity in the energetic γ-ray region has improved considerably and is about to increase further in the near future, motivating a detailed calculation of high-energy (HE; ≥100 MeV) and very high-energy (VHE; ≥100 GeV) γ-ray emission from the nearby starburst galaxy NGC 253. Adopting the convection–diffusion model for energetic electron and proton propagation, and accounting for all the relevant hadronic and leptonic processes, we determine the steady-state energy distributions of these particles by a detailed numerical treatment. The electron distribution is directly normalized by the measured synchrotron radio emission from the central starburst region; a commonly expected theoretical relation is then used to normalize the proton spectrum in this region. Doing so fully specifies the electron spectrum throughout the galactic disc and, with an assumed spatial profile of the magnetic field, the predicted radio emission from the full disc matches well the observed spectrum, confirming the validity of our treatment. The resulting radiative yields of both particles are calculated; the integrated HE and VHE fluxes from the entire disc are predicted to be   f (≥100 MeV) ≃ (1.8^+1.5_−0.8) × 10⁻⁸ cm⁻² s⁻¹  and   f (≥100 GeV) ≃ (3.6^+3.4_−1.7) × 10⁻¹² cm⁻² s⁻¹  , with a central magnetic field value   B ₀≃ 190 ± 10 μ  G. We discuss the feasibility of measuring emission at these levels with the space-borne Fermi and ground-based Cherenkov telescopes.     

Figure rotation of dark haloes in cold dark matter simulations

We investigate the figure rotation of dark matter haloes identified in Λ cold dark matter (CDM) simulations. We find that when strict criteria are used to select suitable haloes for study, five of the 222 haloes identified in our   z = 0  simulation output undergo coherent figure rotation over a  5 h ⁻¹ Gyr  period. We discuss the effects of varying the selection criteria and find that pattern speeds for a much larger fraction of the haloes can be measured when the criteria are relaxed. Pattern speeds measured over a  1 h ⁻¹ Gyr  period follow a lognormal distribution, centred at  Ω_p= 0.2 h rad Gyr⁻¹  with a maximum value of 0.94 h rad Gyr⁻¹. Over a  5 h ⁻¹ Gyr  period, the average pattern speed of a halo is about  0.1 h rad Gyr⁻¹  and the largest pattern speed found is  0.24 h rad Gyr⁻¹  . Less than half of the selected haloes showed alignment between their figure rotation axis and minor axis, the exact fraction being somewhat dependent on how one defines a halo. While the pattern speeds observed are lower than those generally thought capable of causing spiral structure, we note that coherent figure rotation is found over very long periods and argue that further simulations would be required before strong conclusions about spiral structure in all galaxies could be drawn. We find no correlation between halo properties such as total mass and the pattern speed.     

The coarse geometry of merger trees in Λ cold dark matter

We introduce the contour process to describe the geometrical properties of merger trees. The contour process translates a tree structure into a one-dimensional object: the contour walk. We characterize the contour walk measuring its length and action. The length is proportional to the number of progenitors in the tree, and the action is a proxy for the mean length of a branch in the tree.
We construct the contour walk for merger trees extracted from the public data base of the Millennium Run and also for merger trees constructed with a public Monte Carlo code which implements a Markovian algorithm. The trees correspond to haloes of final masses between 10¹¹ and  10¹⁴  h ⁻¹ M_⊙  . We study how the length and action scale with the mass of the final halo.
In all cases, except for the action measured from Markovian trees, we find a transitional scale around  1–6 × 10¹²  h ⁻¹ M_⊙  , which can be interpreted as the point where the halo growth starts to be dominated by mergers and not accretion. The transitional scale depends weakly on the mass resolution and time discretization of the merging trees.     

Power spectrum shape from peculiar velocity data

We constrain the velocity power spectrum shape parameter Γ in linear theory using the nine bulk flow and shear moments estimated from four recent peculiar velocity surveys. For each survey, a likelihood function for Γ was found after marginalizing over the power spectrum amplitude  σ₈Ω^0.6_m  using constraints obtained from comparisons between redshift surveys and peculiar velocity data. In order to maximize the accuracy of our analyses, the velocity noise σ_* was estimated directly for each survey. A statistical analysis of the differences between the values of the moments estimated from different surveys showed consistency with theoretical predictions, suggesting that all the surveys investigated reflect the same large-scale flows. The peculiar velocity surveys were combined into a composite survey yielding the constraint  Γ= 0.13^+0.09_−0.05  . This value is lower than, but consistent with, values obtained using redshift surveys and cosmic microwave background data.     

Halo masses for optically selected and for radio-loud AGN from clustering and galaxy–galaxy lensing

We compute two-point correlation functions and measure the shear signal due to galaxy–galaxy lensing for 80 000 optically identified and 5700 radio-loud active galactic nuclei (AGN) from Data Release 4 of the Sloan Digital Sky Survey. Halo occupation models are used to estimate halo masses and satellite fractions for these two types of AGN. The large sample size allows us to separate AGN according to the stellar mass of their host galaxies. We study how the halo masses of optical and radio AGN differ from those of the parent population at fixed   M _*  . Halo masses deduced from clustering and from lensing agree satisfactorily. Radio AGN are found in more massive haloes than optical AGN: in our samples, their mean halo masses are  1.6 × 10¹³  and  8 × 10¹¹  h ⁻¹ M_⊙  , respectively. Optical AGN follow the same relation between stellar mass and halo mass as galaxies selected without regard to nuclear properties, but radio-loud AGN deviate significantly from this relation. The dark matter haloes of radio-loud AGN are about twice as massive as those of control galaxies of the same stellar mass. This boost is independent of radio luminosity, and persists even when our analysis is restricted to field galaxies. The large-scale gaseous environment of the galaxy clearly plays a crucial role in producing observable radio emission. The dark matter halo masses that we derive for the AGN in our two samples are in good agreement with recent models in which feedback from radio AGN becomes dominant in haloes where gas cools quasi-statically.     

Revisiting the baryon fractions of galaxy clusters: a comparison with WMAP 3-yr results

The universal baryonic mass fraction  (Ω_b/Ω_m)  can be sensitively constrained using X-ray observations of galaxy clusters. In this paper, we compare the baryonic mass fraction inferred from measurements of the cosmic microwave background with the gas mass fractions ( f _gas) of a large sample of clusters taken from the recent literature. In systems cooler than 4 keV, f _gas declines as the system temperature decreases. However, in higher temperature systems, f _gas( r ₅₀₀) converges to  ≈(0.12 ± 0.02)( h /0.72)^−1.5  , where the uncertainty reflects the systematic variations between clusters at r ₅₀₀. This is significantly lower than the maximum-likelihood value of the baryon fraction from the recently released Wilkinson Microwave Anisotropy Probe ( WMAP ) 3-yr results. We investigate possible reasons for this discrepancy, including the effects of radiative cooling and non-gravitational heating, and conclude that the most likely solution is that Ω_m is higher than the best-fitting WMAP value (we find  Ω_m= 0.36^+0.11_−0.08  ), but consistent at the 2σ level. Degeneracies within the WMAP data require that σ₈ must also be greater than the maximum likelihood value for consistency between the data sets.