Constraints on cosmological parameters from recent measurements of cosmic microwave background anisotropy

A key prediction of cosmological theories for the origin and evolution of structure in the Universe is the existence of a 'Doppler peak' in the angular power spectrum of cosmic microwave background (CMB) fluctuations. We present new results from a study of recent CMB observations which provide the first strong evidence for the existence of a 'Doppler peak' localized in both angular scale and amplitude. This first estimate of the angular position of the peak is used to place a new direct limit on the curvature of the Universe, corresponding to a density of Ω = 0.7^+0.8_−0.5, consistent with a flat universe. Very low-density 'open' universe models are inconsistent with this limit unless there is a significant contribution from a cosmological constant. For a flat standard cold dark matter dominated universe we use our results in conjunction with big bang nucleosynthesis constraints to determine the value of the Hubble constant as H ₀ = 30 − 70 km s⁻¹ Mpc⁻¹ for baryon fractions Ω_b = 0.05 to 0.2. For H ₀ = 50 km s⁻¹ Mpc⁻¹ we find the primordial spectral index of the fluctuations to be n  = 1.1 ± 0.1, in close agreement with the inflationary prediction of n  ≃ 1.0.     

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  .     

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.     

Observations of the Sunyaev–Zel'dovich effect in the z=0.78 cluster MS 1137.5+6625

We have observed the   z =0.78  cluster MS 1137.5+6625 with the Ryle Telescope (RT) at 15 GHz. After subtraction of contaminating radio sources in the field, we find a Sunyaev–Zel'dovich flux decrement of  -421±60 μJy  on the ≈0.65 k λ baseline of the RT, spatially coincident with the optical and X-ray positions for the cluster core.
For a spherical King-profile cluster model, the best fit to our flux measurement has a core radius   θ _C=20 arcsec  , consistent with previous X-ray observations, and a central temperature decrement  Δ T =650±92 μK  .
Using this model, we calculate that the cluster has a gas mass inside a     radius of  2.9×10¹³ M_⊙  for an  Ω _M =1  universe and  1.6×10¹³ M_⊙  for  Ω _M =0.3  ,  Ω_Λ=0.7  . We compare this model with existing measurements of the total mass of the cluster, based on gravitational lensing, and estimate a gas fraction for MS 1137.5+6625 of ≈8 per cent.     

Survival and disruption of subsystems during a cold collapse

Cold collapse of a cluster composed of small identical clumps, each of which is in virial equilibrium, is considered. Since the clumps have no relative motion with respect to each other initially, the cluster collapses under its own gravity. At the first collapse of the cluster, most of the clumps are destroyed, but some survive. In order to find the condition for the clumps to survive, we made a systematic study in two-parameter space: the number of the clumps N _c and the size of the clump r _v . We obtained the condition N _c ≫ 1 and n _k  ≥ 1, where n _k is related to r _v and the initial radius of the cluster R _ini through the relation R _ini/ r _v  = 2 N ^{( n k +5)/6}_c. A simple analytical argument supports the numerical result. This n _k corresponds to the index of the power spectrum of the density fluctuation in the cosmological hierarchical clustering, and thus our result may suggest that in the systems smaller than 2/Ω h ²)Mpc, the first violent collapse is strong enough to sweep away all the substructures that exist before the collapse.     

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.     

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 dipole anisotropy of the 2 Micron All-Sky Redshift Survey

We estimate the acceleration on the Local Group (LG) from the 2 Micron All-Sky Redshift Survey (2MRS). The sample used includes about 23 200 galaxies with extinction-corrected magnitudes brighter than   K _s= 11.25  and it allows us to calculate the flux-weighted dipole. The near-infrared flux-weighted dipoles are very robust because they closely approximate a mass-weighted dipole, bypassing the effects of redshift distortions and require no preferred reference frame. This is combined with the redshift information to determine the change in dipole with distance. The misalignment angle between the LG and the cosmic microwave background (CMB) dipole drops to  12°± 7°  at around  50  h ⁻¹ Mpc  , but then increases at larger distances, reaching  21°± 8°  at around  130  h ⁻¹ Mpc  . Exclusion of the galaxies Maffei 1, Maffei 2, Dwingeloo 1, IC342 and M87 brings the resultant flux dipole to  14°± 7°  away from the CMB velocity dipole. In both cases, the dipole seemingly converges by  60  h ⁻¹ Mpc  . Assuming convergence, the comparison of the 2MRS flux dipole and the CMB dipole provides a value for the combination of the mass density and luminosity bias parameters  Ω^0.6_m/ b _L= 0.40 ± 0.09  .     

The effective acceleration of plasma outflow in the paraboloidal magnetic field

The problem of the efficiency of particle acceleration for a paraboloidal poloidal magnetic field is considered within the approach of steady axisymmetric magnetohydrodynamic (MHD) flow. For the large Michel magnetization parameter σ it is possible to linearize the stream equation near the force-free solution and to solve the problem self-consistently as was done by Beskin, Kuznetsova & Rafikov for a monopole magnetic field. It is shown that, on the fast magnetosonic surface (FMS), the particle Lorentz factor γ does not exceed the standard value  σ^1/3  . On the other hand, in the supersonic region, the Lorentz factor grows with the distance z from the equatorial plane as  γ≈ ( z / R _L)^1/2  up to the distance   z ≈σ² R _L  , where   R _L= c /Ω_F  is the radius of the light cylinder. Thus, the maximal Lorentz factor is  γ_max≈σ  , which corresponds to almost the full conversion of the Poynting energy flux into the particle kinetic one.     

Detection of OH and wide H i absorption toward B0218+357

Large-scale polarization of the cosmic microwave background measured by the WMAP satellite requires a mean optical depth to Thomson scattering,  τ_e∼ 0.17  . The reionization of the Universe must therefore have begun at relatively high redshift. We have studied the reionization process using supercomputer simulations of a large and representative region of a universe which has cosmological parameters consistent with the WMAP results (  Ω_m= 0.3, Ω_Λ= 0.7, h = 0.7, Ω_b= 0.04, n = 1  and  σ₈= 0.9  ). Our simulations follow both the radiative transfer of ionizing photons and the formation and evolution of the galaxy population which produces them. A previously published model with ionizing photon production as expected for zero-metallicity stars distributed according to a standard stellar initial mass function (IMF) (10⁶¹ photons per unit solar mass of formed stars) and with a moderate photon escape fraction from galaxies (5 per cent), produces  τ_e= 0.104  , which is within 1.0 to  1.5σ  of the 'best' WMAP value. Values of up to 0.16 can be produced by taking larger escape fractions or a top-heavy IMF. The data do not require a separate populations of 'miniquasars' or of stars forming in objects with total masses below  10⁹ M_⊙  . Reconciling such early reionization with the observed Gunn–Peterson troughs in   z > 6  quasars may be challenging. Possible resolutions of this problem are discussed.     

Modelling the bar in the centre of the starburst galaxy M82

We present VLA A-array 21-cm atomic hydrogen (H  i ) absorption observed against the central region of the starburst galaxy M82 with an angular resolution of ∼1.3 arcsec (≃20 pc). These observations, together with MERLIN H  i absorption measurements, are compared with the molecular (CO) and ionized ([Ne  ii ]) gas distributions and are used to constrain the dynamics and structure of the ionized, neutral and molecular gas in this starburst.
A position–velocity diagram of the H  i distribution reveals an unusual 'hole' feature which, when previously observed in CO, has been interpreted as an expanding superbubble contained within a ring of gas in solid body rotation. However, we interpret this feature as a signature of a nearly edge-on barred galaxy. In addition, we note that the CO, H  i and [Ne  ii ] position–velocity diagrams reveal two main velocity gradients, and we interpret these as gas moving on x₁- and x₂-orbits within a bar potential. We find the best fit to the data to be produced using a bar potential with a flat rotation curve velocity v _b=140 km s⁻¹ and a total length of 1 kpc, a non-axisymmetry parameter q =0.9, an angular velocity of the bar Ω_b=217 km s⁻¹ arcsec⁻¹, a core radius R _c=25 pc, an inclination angle i =80° and a projected angle between the bar and the major axis of the galaxy φ '=4°. We also discuss the orientation of the disc and bar in M82.     

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.     

Spectroscopic and photometric observations of the eclipsing star UV Leo

High-resolution spectroscopic observations around the Hα line and BVRI photometry of the eclipsing short-period RS CVn star UV Leo are presented. The simultaneous light-curve solution and radial velocity-curve solution led to the following values of the global parameters of the binary: temperatures   T ₁= 6000 ± 100 K  and   T ₂= 5970 ± 20 K  ; masses   M ₁= 0.976 ± 0.067 M_⊙  and   M ₂= 0.931 ± 0.052 M_⊙  ; separation   a = 3.716 ± 0.048 R_⊙  ; orbital inclination     ; radii   R ₁= 1.115 ± 0.052 R_⊙  and   R ₂= 1.078 ± 0.051 R_⊙  ; equatorial velocities   V ₁= 98.8 ± 2.3 km s⁻¹  and   V ₂= 89.6 ± 2.7 km s⁻¹  . These results lead to the conclusion that the two components of UV Leo are slightly oversized for their masses and lie within the main-sequence band on the mass–radius diagram, close to the isochrone 9 × 10¹⁰ yr.     

Masses for the Local Group and the Milky Way

We use the very large Millennium Simulation of the concordance Λ cold dark matter cosmogony to calibrate the bias and error distribution of Timing Argument estimators of the masses of the Local Group and of the Milky Way. From a large number of isolated spiral–spiral pairs similar to the Milky Way/Andromeda system, we find the interquartile range of the ratio of timing mass to true mass to be a factor of 1.8, while the 5 and 95 per cent points of the distribution of this ratio are separated by a factor of 5.7. Here, we define true mass as the sum of the 'virial' masses, M ₂₀₀, of the two dominant galaxies. For present best values of the distance and approach velocity of Andromeda, this leads to a median likelihood estimate of the true mass of the Local Group of  5.27 × 10¹² M_⊙  or  log  M _LG/M_⊙= 12.72  , with an interquartile range of [12.58, 12.83] and a 5–95 per cent range of [12.26, 13.01]. Thus, a 95 per cent lower confidence limit on the true mass of the Local Group is  1.81 × 10¹² M_⊙  . A timing estimate of the Milky Way's mass based on the large recession velocity observed for the distant satellite Leo I works equally well, although with larger systematic uncertainties. It gives an estimated virial mass for the Milky Way of  2.43 × 10¹² M_⊙  with a 95 per cent lower confidence limit of  0.80 × 10¹² M_⊙  .     

The Gunn–Peterson effect and the Lyman alpha forest

We show that spatial correlations in a stochastic large-scale velocity field in an otherwise smooth intergalactic medium (homogeneous comoving density) superposed on the general Hubble flow may cause a 'line-like' structure in QSO spectra similar to the population of unsaturated Lyα forest lines which usually are attributed to individual clouds with 10¹¹ ≲ N _{H i}  5 × 10¹³ cm⁻². Therefore there is no clear observational distinction between a diffuse intergalactic medium and discrete intergalactic clouds. It follows that the H  i density in the diffuse intergalactic medium might be substantially underestimated if it is determined from the observed intensity distribution near the apparent continuum in high-resolution spectra of QSOs. Our tentative estimate implies a diffuse neutral hydrogen opacity τ_GP ∼ 0.3 at z  ∼ 3 and a current baryon density Ω_IGM ≃ 0.08, assuming a Hubble constant H ₀ = 70 km s⁻¹ Mpc⁻¹.     

Optimal filters for the construction of the ensemble pulsar time

An algorithm of the ensemble pulsar time based on the optimal Wiener filtration method has been constructed. This algorithm allows the separation of the contributions to the post-fit pulsar timing residuals of the atomic clock and the pulsar itself. Filters were designed using the cross- and auto-covariance functions of the timing residuals. The method has been applied to the timing data of millisecond pulsars PSR B1855+09 and B1937+21 and allowed the filtering out of the atomic-scale component from the pulsar data. Direct comparison of the terrestrial time TT(BIPM06) and the ensemble pulsar time PT_ens revealed that the fractional instability of TT(BIPM06)−PT_ens is equal to  σ _z = (0.8 ± 1.9) × 10⁻¹⁵  . Based on the  σ _z   statistics of TT(BIPM06)−PT_ens, a new limit of the energy density of the gravitational wave background was calculated to be equal to  Ω_g h ²∼ 3 × 10⁻⁹  .     

Continuous frequency spectrum of the global hydromagnetic oscillations of a magnetically confined mountain on an accreting neutron star

We compute the continuous part of the ideal-magnetohydrodynamic (ideal-MHD) frequency spectrum of a polar mountain produced by magnetic burial on an accreting neutron star. Applying the formalism developed by Hellsten & Spies, extended to include gravity, we solve the singular eigenvalue problem subject to line-tying boundary conditions. This spectrum divides into an Alfvén part and a cusp part. The eigenfunctions are chirped and anharmonic with an exponential envelope, and the eigenfrequencies cover the whole spectrum above a minimum ω_low. For equilibria with accreted mass  1.2 × 10⁻⁶≲ M _a/M_⊙≲ 1.7 × 10⁻⁴  and surface magnetic fields  10¹¹≲ B _*/G ≲ 10¹³, ω_low  is approximately independent of   B _*  , and increases with M _a. The results are consistent with the Alfvén spectrum excited in numerical simulations with the zeus-mp solver. The spectrum is modified substantially by the Coriolis force in neutron stars spinning faster than ∼100 Hz. The implications for gravitational-wave searches for low-mass X-ray binaries are considered briefly.     

The cooling of shock-compressed primordial gas

We find that at redshifts   z ≳ 10, HD  line cooling allows strongly shocked primordial gas to cool to the temperature of the cosmic microwave background (CMB). This temperature is the minimum value attainable via radiative cooling. Provided that the abundance of HD, normalized to the total number density, exceeds a critical level of  ∼10⁻⁸  , the CMB temperature floor is reached in a time which is short in comparison to the Hubble time. We estimate the characteristic masses of stars formed out of shocked primordial gas in the wake of the first supernovae, and resulting from the virialization of dark matter haloes during hierarchical structure formation to be  ∼10 M_⊙  . In addition, we show that cooling by HD enables the primordial gas in relic H  ii regions to cool to temperatures considerably lower than those reached via H₂ cooling alone. We confirm that HD cooling is unimportant in cases where the primordial gas does not go through an ionized phase, as in the formation process of the very first stars in   z ≳ 20  minihaloes of mass  ∼10⁶ M_⊙  .     

On the MHD effects on the force-free monopole outflow

The stationary axisymmetric outflow from a rotating sphere with a (split) monopole magnetic field is considered. The stream equation describing the outflow is linearized in terms of the Michel magnetization parameter σ⁻¹ ≪ 1, which allows a self-consistent analysis of the direct problem. It is shown that for a finite σ the fast magnetosonic surface is located at a finite distance ∼ σ^1/3 R _L ( R _L =  c /Ω_F is the light cylinder). We have also found that the particle energy at the fast surface is just equal to the Michel value γ ∼ 1/3σ. The particle acceleration and magnetic field collimation are shown to become ineffective outside the fast magnetosonic surface.     

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.