Analytic solutions for coupled linear perturbations

Analytic solutions for the evolution of cosmological linear density perturbations in the baryonic gas and collisionless dark matter are derived. The solutions are expressed in a closed form in terms of elementary functions, for arbitrary baryonic mass fraction. They are obtained assuming =1 and a time-independent comoving Jeans wavenumber, k _J. By working with a time variable ln( t ^2/3), the evolution of the perturbations is described by linear differential equations with constant coefficients. The new equations are then solved by means of Laplace transformation, assuming that the gas and dark matter trace the same density field before a sudden heating epoch. In a dark matter-dominated Universe, the ratio of baryonic to dark matter density perturbation decays with time roughly as exp(5 /4) t ^5/6 to the limiting value 1/[1+( k k _J)²]. For wavenumbers the decay is accompanied by oscillations with a period in . In comparison, as increases in a baryonic matter-dominated Universe, the ratio approaches 1( k k _J)² for k k _J, and zero otherwise.     

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.     

Microwave polarization in the direction of galaxy clusters induced by the CMB quadrupole anisotropy

Electron scattering induces a polarization in the cosmic microwave background (CMB) signal measured in the direction of a galaxy cluster owing to the presence of a quadrupole component in the CMB temperature distribution. Measuring the polarization towards distant clusters provides the unique opportunity to observe the evolution of the CMB quadrupole at moderate redshifts, z ∼0.5–3. We demonstrate that for the local cluster population the polarization degree will depend on the cluster celestial position. There are two extended regions in the sky, which are opposite to each other, where the polarization is maximal, ∼0.1( τ /0.02) μK in the Rayleigh–Jeans part of the CMB spectrum ( τ being the Thomson optical depth across the cluster). This value exceeds the polarization introduced by the cluster transverse peculiar motion if v _t<1300 km s⁻¹. One can hope to detect this small signal by measuring a large number of clusters, thereby effectively removing the systematic contribution from other polarization components produced in clusters. These polarization effects, which are of the order of ( v _t c )² τ , ( v _t c ) τ ² and ( kT _e m _e c ²) τ ², as well as the polarization owing to the CMB quadrupole, were previously given by Sunyaev and Zel'dovich for the Rayleigh–Jeans part of the spectrum. We fully confirm their earlier results and present exact frequency dependences for all these effects. The polarization degree is considerably higher in the Wien region.     

Electron density diagnostic potential of the Si x ion and its application in Procyon

Theoretical electron density sensitive line ratios   R ₁– R ₆  of Si  x soft X-ray emission lines are presented. We found that these line ratios are sensitive to electron density n _e, and the ratio R ₁ is insensitive to electron temperature T _e. For reliable determination of the electron density of laboratory and astrophysical plasmas, atomic data, such as electron impact excitation rates, are very important. Our results reveal that the discrepancy of the line ratios from different atomic data calculated with the distorted wave (DW) approximation and the R-matrix method is up to 19 per cent at   n _e= 2 × 10⁸ cm⁻³  . We applied the theoretical intensity ratio R ₁ to the Low Energy Transmission Grating Spectrometer (LETGS) spectrum of the solar-like star Procyon. By comparing the observed value (1.29) with the theoretical calculation, the derived electron density n _e is  2.6 × 10⁸ cm⁻³  , which is consistent with that derived from  (C  v < 8.3 × 10⁸ cm⁻³)  . When the temperature structure of the Procyon corona is taken into account, the derived electron density increases from   n _e= 2.6 × 10⁸  to  2.8 × 10⁸ cm⁻³  .     

Non-linear evolution of suppressed dark matter primordial power spectra

We address the degree and rapidity of generation of small-scale power over the course of structure formation in cosmologies where the primordial power spectrum is strongly suppressed beyond a given wavenumber. We first summarize the situations where one expects such suppressed power spectra and point out their diversity. We then employ an exponential cut-off, which characterizes warm dark matter (WDM) models, as a template for the shape of the cut-off and focus on damping scales ranging from 10⁶ to  10⁹  h ⁻¹ M_⊙  . Using high-resolution simulations, we show that the suppressed part of the power spectrum is quickly (re)generated and catches up with both the linear and the non-linear evolution of the unsuppressed power spectrum. From   z = 2  onwards, a power spectrum with a primordial cut-off at  10⁹  h ⁻¹ M_⊙  becomes virtually indistinguishable from an evolved cold dark matter (CDM) power spectrum. An attractor such as that described in Zaldarriaga, Scoccimarro & Hui for power spectra with different spectral indices also emerges in the case of truncated power spectra. Measurements of   z ∼ 0  non-linear power spectra at  ∼100  h ⁻¹ kpc  cannot rule out the possibility of linear power spectra damped below  ∼10⁹  h ⁻¹ M_⊙  . Therefore, WDM or scenarios with similar features should be difficult to exclude in this way.     

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.     

Constraints on perfect fluid and scalar field dark energy models from future redshift surveys

We discuss the constraints that future photometric and spectroscopic redshift surveys can put on dark energy through the baryon oscillations of the power spectrum. We model the dark energy either with a perfect fluid or a scalar field and take into account the information contained in the linear growth function. We show that the growth function helps to break the degeneracy in the dark energy parameters and reduce the errors on   w ₀, w ₁  roughly by 30 per cent, making more appealing multicolour surveys based on photometric redshifts. We find that a 200-deg² spectroscopic survey reaching   z ≈ 3  can constrain   w ₀, w ₁  to within  Δ w ₀= 0.21, Δ w ₁= 0.26  , to  Δ w ₀= 0.39, Δ w ₁= 0.54  using photometric redshifts with an absolute uncertainty of 0.02, and to  Δ w ₀= 0.43, Δ w ₁= 0.66  with an uncertainty of 0.04. In the scalar field case, we show that the slope n of the inverse power-law potential for dark energy can be constrained to  Δ n = 0.26  (spectroscopic redshifts) or  Δ n = 0.40  (photometric redshifts), i.e. better than with future ground-based supernovae surveys or cosmic microwave background data.     

A Bayesian Kepler periodogram detects a second planet in HD 208487

An automatic Bayesian Kepler periodogram has been developed for identifying and characterizing multiple planetary orbits in precision radial velocity data. The periodogram is powered by a parallel tempering Markov chain Monte Carlo (MCMC) algorithm which is capable of efficiently exploring a multiplanet model parameter space. The periodogram employs an alternative method for converting the time of an observation to true anomaly that enables it to handle much larger data sets without a significant increase in computation time. Improvements in the periodogram and further tests using data from HD 208487 have resulted in the detection of a second planet with a period of 909⁸²₋₉₂ d, an eccentricity of 0.37^0.26_−0.20, a semimajor axis of 1.87^0.13_−0.14 au and an M sin  i = 0.45^+0.11_−0.13 M _J. The revised parameters of the first planet are period = 129.8 ± 0.4 d, eccentricity = 0.20 ± 0.09, semimajor axis = 0.51 ± 0.02 au and M sin  i = 0.41 ± 0.05  M _J. Particular attention is paid to several methods for calculating the model marginal likelihood which is used to compare the probabilities of models with different numbers of planets.     

Principal-component analysis of the cosmic microwave background anisotropies: revealing the tensor degeneracy

A principal-component analysis of cosmic microwave background (CMB) anisotropy measurements is used to investigate degeneracies among cosmological parameters. The results show that a degeneracy with tensor modes – the 'tensor degeneracy'– dominates uncertainties in estimates of the baryon and cold dark matter densities,   ω _b=Ω_b  h ²  ,   ω _c=Ω_c  h ²  , ¹ from an analysis of CMB anisotropies alone. The principal-component analysis agrees well with a maximum-likelihood analysis of the observations, identifying the main degeneracy directions and providing an impression of the effective dimensionality of the parameter space.     

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.     

Consequences for some dark energy candidates from the type Ia supernova SN 1997ff

We examine the status of various dark energy models in light of the recently observed SN 1997ff at   z ≈1.7  . The modified data still fit a pure cosmological constant Λ or a quintessence with an equation of state similar to that of Λ. The kinematical Λ models,  Λ∼ S ^-2  and  Λ∼ H ²  , also fit the data reasonably well and require less dark energy density (hence more matter energy density) than is required by the constant Λ model. However, the model  Λ∼ S ^-2  with low energy density becomes unphysical as it cannot accommodate higher redshift objects.
We also examine an alternative explanation of the data, namely the absorption by the intervening whisker-like dust, and find that the quasi-steady state (QSS) model and the Friedmann–Robertson–Walker (FRW) model  Ω_m0=0.33  without any dark energy also fit the data reasonably well.
We notice that the addition of SN 1997ff to the old data has worsened the fit to most of the models, except a closed FRW model with a constant Λ and a closed quintessence model with   ω _φ =-0.82  , and the models have started departing from each other as we go above   z =1  . However, to make a clear discrimination possible, a few more supernovae with   z >1  are required.
We have also calculated the age of the Universe in these models and find that, in the models with a constant Λ, the expansion age is uncomfortably close to the age of the globular clusters. Quintessence models show even lower age. The kinematical Λ models are, however, interesting in this connection (especially the model  Λ∼ H ²)  , as they give a remarkably large age of the Universe.     

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

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

Maximum-likelihood method for estimating the mass and period distributions of extrasolar planets

We investigate the distribution of mass M and orbital period P of extrasolar planets, taking account of selection effects caused by the limited velocity precision and duration of existing surveys. We fit the data on 72 planets to a power-law distribution of the form  d n = CM ^−α P ^−β(d M / M )(d P / P )  , and find  α= 0.11 ± 0.10  ,  β=−0.27 ± 0.06  for   M ≲ 10 M _J  , where   M _J  is the mass of Jupiter. The correlation coefficient between these two exponents is −0.31, indicating that uncertainties in the two distributions are coupled. We estimate that 4 per cent of solar-type stars have companions in the range  1 M _J < M < 10 M _J  ,  2_d < P < 10 yr  .     

Constraints on a non-Gaussian cold dark matter model

We consider constraints on the structure formation model based on non-Gaussian fluctuations generated during inflation, which have     distributions. Using three data sets, the abundance of the clusters at z =0, moderate z and the correlation length, we show that constraints on the non-Gaussianity and the amplitude of fluctuations and the density parameter can be obtained. We obtain an upper bound for _m, and a lower bound for the non-Gaussianity and the amplitude of the fluctuations. Using the abundance of clusters at z 0.6, for the spectrum parametrized by cold dark matter (CDM) shape parameter =0.23, we obtain an upper bound for the density parameter of _m0.5 and lower bounds for the amplitude of ₈0.7 and for the non-Gaussianity of fluctuations of G 2 ( m 200), where G =1 for Gaussian.     

Neutralino annihilations and the gas temperature in the dark ages

Assuming that the dark matter is entirely made up of neutralinos, we re-visit the role of their annihilation on the temperature of diffuse gas in the high-redshift universe  ( z > 10)  , before the formation of luminous structures. We consider neutralinos of particle mass 36 and 100 GeV. The former is able to produce  ∼7  e ⁻ e ⁺  particles per annihilation through the fremionic channel, and the latter ∼53 particles assuming a purely bosonic channel. High-energy   e ⁻ e ⁺  particles up-scatter the cosmic microwave background (CMB) photons into higher energies via the inverse-Compton scattering. The process produces a power-law   e ⁻ e ⁺  energy spectrum of index −1 in the energy range of interest, independent of the initial energy distribution. The corresponding energy spectrum of the up-scattered photons is a power law of index −1/2, if absorption by the gas is not included. The scattered photons photoheat the gas by releasing electrons which deposit a fraction (14 per cent) of their energy as heat into the ambient medium. For uniformly distributed neutralinos, the heating is insignificant. The effect is greatly enhanced by the clumping of neutralinos into dense haloes. We use a time-dependent clumping model which takes into account the damping of density fluctuations on mass-scales smaller than  ∼10⁻⁶ M_⊙  . With this clumping model, the heating mechanism boosts the gas temperature above that of the CMB after a redshift of   z ∼ 30  . By   z ≈ 10  , the gas temperature is nearly 100 times its temperature when no heating is invoked. Similar increase is obtained for the two neutralino masses considered.     

Clustering of galaxy clusters in cold dark matter universes

We use very large cosmological N -body simulations to obtain accurate predictions for the two-point correlations and power spectra of mass-limited samples of galaxy clusters. We consider two currently popular cold dark matter (CDM) cosmogonies, a critical density model ( τ CDM) and a flat low density model with a cosmological constant (ΛCDM). Our simulations each use 10⁹ particles to follow the mass distribution within cubes of side 2  h ⁻¹ Gpc ( τ CDM) and 3  h ⁻¹ Gpc (ΛCDM) with a force resolution better than 10⁻⁴ of the cube side. We investigate how the predicted cluster correlations increase for samples of increasing mass and decreasing abundance. Very similar behaviour is found in the two cases. The correlation length increases from     for samples with mean separation     to     for samples with     The lower value here corresponds to τ CDM and the upper to ΛCDM. The power spectra of these cluster samples are accurately parallel to those of the mass over more than a decade in scale. Both correlation lengths and power spectrum biases can be predicted to better than 10 per cent using the simple model of Sheth, Mo & Tormen. This prediction requires only the linear mass power spectrum and has no adjustable parameters. We compare our predictions with published results for the automated plate measurement (APM) cluster sample. The observed variation of correlation length with richness agrees well with the models, particularly for ΛCDM. The observed power spectrum (for a cluster sample of mean separation     ) lies significantly above the predictions of both models.     

The Millennium Galaxy Catalogue: 16 ≤BMGC < 24 galaxy counts and the calibration of the local galaxy luminosity function

The Millennium Galaxy Catalogue (MGC) is a 37.5 deg², medium-deep, B -band imaging survey along the celestial equator, taken with the Wide Field Camera on the Isaac Newton Telescope. The survey region is contained within the regions of both the Two Degree Field Galaxy Redshift Survey (2dFGRS) and the Sloan Digital Sky Survey Early Data Release (SDSS-EDR). The survey has a uniform isophotal detection limit of 26 mag arcsec⁻² and it provides a robust, well-defined catalogue of stars and galaxies in the range  16 ≤ B _MGC < 24 mag  .
Here we describe the survey strategy, the photometric and astrometric calibration, source detection and analysis, and present the galaxy number counts that connect the bright and faint galaxy populations within a single survey. We argue that these counts represent the state of the art and use them to constrain the normalizations (φ*) of a number of recent estimates of the local galaxy luminosity function. We find that the 2dFGRS, SDSS Commissioning Data (CD), ESO Slice Project, Century Survey, Durham/UKST, Mt Stromlo/APM, SSRS2 and NOG luminosity functions require a revision of their published φ* values by factors of  1.05 ± 0.05, 0.76 ± 0.10, 1.02 ± 0.22, 1.02 ± 0.16, 1.16 ± 0.28, 1.75 ± 0.37, 1.40 ± 0.26  and  1.01 ± 0.39  , respectively. After renormalizing the galaxy luminosity functions we find a mean local b _J luminosity density of     . ¹     

A Chandra detection of the radio hotspot of 3C 123

Chandra X-ray Observatory observations of the powerful, peculiar radio galaxy 3C 123 have resulted in an X-ray detection of the bright eastern hotspot, with a 1-keV flux density of ∼5 nJy. The X-ray flux and spectrum of the hotspot are consistent with the X-rays being inverse-Compton scattering of radio synchrotron photons by the population of electrons responsible for the radio emission ('synchrotron self-Compton emission') if the magnetic fields in the hotspot are close to their equipartition values. 3C 123 is thus the third radio galaxy to show X-ray emission from a hotspot which is consistent with being in equipartition. Chandra also detects emission from a moderately rich cluster surrounding 3C 123, with L _X(2–10 keV)=2×10⁴⁴ erg s⁻¹ and kT ∼5 keV, and absorbed emission from the active nucleus, with an inferred intrinsic column density of 1.7×10²² cm⁻² and an intrinsic 2–10 keV luminosity of 10⁴⁴ erg s⁻¹.     

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  .