The topology of the IRAS Point Source Catalogue Redshift Survey

We investigate the topology of the new Point Source Catalogue Redshift Survey (PSCz) of IRAS galaxies by means of the genus statistic. The survey maps the local Universe with approximately 15 000 galaxies over 84.1 per cent of the sky, and provides an unprecedented number of resolution elements for the topological analysis. For comparison with the PSCz data we also examine the genus of large N -body simulations of four variants of the cold dark matter (CDM) cosmogony. The simulations are part of the Virgo project to simulate the formation of structure in the Universe. We assume that the statistical properties of the galaxy distribution can be identified with those of the dark matter particles in the simulations. We extend the standard genus analysis by examining the influence of sampling noise on the genus curve and introducing a statistic able to quantify the amount of phase correlation present in the density field, the amplitude drop of the genus compared to a Gaussian field with identical power spectrum. The results for PSCz are consistent with the hypothesis of random-phase initial conditions. In particular, no strong phase correlation is detected on scales ranging from 10 to 32 h ⁻¹ Mpc, whereas there is a positive detection of phase correlation at smaller scales. Among the simulations, phase correlations are detected in all models at small scales, albeit with different strengths. When scaled to a common normalization, the amplitude drop depends primarily on the shape of the power spectrum. We find that the constant-bias standard CDM model can be ruled out at high significance, because the shape of its power spectrum is not consistent with PSCz. The other CDM models with more large-scale power all fit the PSCz data almost equally well, with a slight preference for a high-density τCDM model.     

The size and shape of local voids

We study the size and shape of low-density regions in the local Universe, which we identify in the smoothed density field of the PSCz flux-limited IRAS galaxy catalogue. After quantifying the systematic biases that enter the detection of voids using our data set and method, we identify, using a smoothing length of 5  h ⁻¹ Mpc, 14 voids within 80  h ⁻¹ Mpc (having volumes 10³  h ⁻³ Mpc³) and, using a smoothing length of 10  h ⁻¹ Mpc, eight voids within 130  h ⁻¹ Mpc (having volumes  8×10³ h⁻³ Mpc³)  . We study the void size distribution and morphologies and find that there is roughly an equal number of prolate and oblate-like spheroidal voids. We compare the measured PSCz void shape and size distributions with those expected in six different cold dark matter (CDM) models and find that only the size distribution can discriminate between models. The models preferred by the PSCz data are those with intermediate values of   σ ₈(≃0.83)  , independent of cosmology.     

The one-point PDF of the initial conditions of our local Universe from the IRAS PSC redshift catalogue

The algorithm ztrace of Monaco & Efstathiou is applied to the IRAS PSCz catalogue to reconstruct the initial conditions of our local Universe with a resolution down to ~5  h ¹ Mpc. The one-point probability distribution function (PDF) of the reconstructed initial conditions is consistent with the assumptions that: (i) IRAS galaxies trace mass on scales of ~5  h ¹ Mpc and (ii) the statistics of the primordial density fluctuations are Gaussian. We use simulated PSCz catalogues, constructed from N -body simulations with Gaussian initial conditions, to show that local non-linear bias can cause the recovered initial PDF (assuming no bias) to be non-Gaussian. However, for plausible bias models, the distortions of the recovered PDF would be difficult to detect using the volume finely sampled by the PSCz catalogue. So, for Gaussian initial conditions, a range of bias models remain compatible with our PSCz reconstruction results.     

Evolution of peaks in weakly non-linear density field and dark halo profiles

Using the two-point Edgeworth series up to second order in the linear rms density fluctuation we construct the weakly non-linear conditional probability distribution function for the density field around an overdense region. This requires calculating the two-point analogues of the skewness parameter S ₃. We test the dependence of the two-point skewness on distance from the peak for scale-free power spectra and Gaussian smoothing. The statistical features of such a conditional distribution are given as the values obtained within linear theory corrected by the terms that arise as a result of weakly non-linear evolution. The expected density around the peak is found to be always below the linear prediction while its dispersion is always larger than in the linear case. For large enough overdensities the weakly non-linear corrections can be more significant than the peak constraint introduced by Bardeen et al. We apply these results to the spherical model of collapse as developed by Hoffman & Shaham and find that in general the effect of weakly non-linear interactions is to decrease the scale from which a peak gathers mass and therefore also the mass itself. In the case of an open universe this results in steepening of the final profile of the virialized proto-object.     

Wavelet analysis and the detection of non-Gaussianity in the cosmic microwave background

We investigate the use of wavelet transforms in detecting and characterizing non-Gaussian structure in maps of the cosmic microwave background (CMB). We apply the method to simulated maps of the KaiserStebbins effect resulting from cosmic strings, on to which Gaussian signals of varying amplitudes are superposed. We find that the method significantly outperforms standard techniques based on measuring the moments of the pixel temperature distribution. We also compare the results with those obtained using techniques based on Minkowski functionals, and we again find the wavelet method to be superior. In particular, using the wavelet technique, we find that it is possible to detect non-Gaussianity even in the presence of a superposed Gaussian signal with 3 times the rms amplitude of the original cosmic string map. We also find that the wavelet technique is useful in characterizing the angular scales at which the non-Gaussian signal occurs.     

Reconstruction of cosmological initial conditions from galaxy redshift catalogues

We present and test a new method for the reconstruction of cosmological initial conditions from a full-sky galaxy catalogue. This method, called ZTRACE, is based on a self-consistent solution of the growing mode of gravitational instabilities according to the Zel'dovich approximation and higher order in Lagrangian perturbation theory. Given the evolved redshift-space density field, smoothed on some scale, ZTRACE finds, via an iterative procedure, an approximation to the initial density field for any given set of cosmological parameters; real-space densities and peculiar velocities are also reconstructed. The method is tested by applying it to N -body simulations of an Einstein–de Sitter and an open cold dark matter universe. It is shown that errors in the estimate of the density contrast dominate the noise of the reconstruction. As a consequence, the reconstruction of real-space density and peculiar velocity fields using non-linear algorithms is little improved over those based on linear theory. The use of a mass-preserving adaptive smoothing, equivalent to a smoothing in Lagrangian space, allows an unbiased (although noisy) reconstruction of initial conditions, as long as the (linearly extrapolated) density contrast does not exceed unity. The probability distribution function of the initial conditions is recovered to high precision, even for Gaussian smoothing scales of ∼5  h ⁻¹ Mpc, except for the tail at δ ≥1. This result is insensitive to the assumptions of the background cosmology.     

Cosmological implications of the PSCz PDF and its moments

We compare the probability density function (PDF) and its low-order moments (variance and skewness) of the smoothed IRAS Point Source Catalogue Redshift Survey (PSC z ) galaxy density field and of the corresponding simulated PSC z look-alikes, generated from N -body simulations of six different dark matter models: four structure-normalized with     and     , one COBE -normalized, and the old standard cold dark matter model. The galaxy distributions are smoothed with a Gaussian window at three different smoothing scales,     , 10 and 15  h ⁻¹ Mpc. We find that the simulation PSC z look-alike PDFs are sensitive only to the normalization of the power spectrum, probably owing to the shape similarity of the simulated galaxy power spectrum on the relevant scales. We find that the only models that are consistent, at a high significance level, with the observed PSC z PDF are models with a relatively low power spectrum normalization     . From the phenomenologically derived σ ₈–moments relation, fitted from the simulation data, we find that the PSC z moments suggest     .     

On the mass of dense star clusters in starburst galaxies from spectrophotometry

The mass of unresolved young star clusters derived from spectrophotometric data may well be off by a factor of 2 or more once the migration of massive stars driven by mass segregation is accounted for. We quantify this effect for a large set of cluster parameters, including variations in the stellar initial mass function (IMF), the intrinsic cluster mass, and mean mass density. Gas-dynamical models coupled with the Cambridge stellar evolution tracks allow us to derive a scheme to recover the real cluster mass given measured half-light radius, one-dimensional velocity dispersion and age. We monitor the evolution with time of the ratio of real to apparent mass through the parameter η. When we compute η for rich star clusters, we find non-monotonic evolution in time when the IMF stretches beyond a critical cut-off mass of  25.5 M_⊙  . We also monitor the rise of colour gradients between the inner and outer volume of clusters: we find trends in time of the stellar IMF power indices overlapping well with those derived for the Large Magellanic Cloud cluster NGC 1818 at an age of 30 Myr. We argue that the core region of massive Antennae clusters should have suffered from much segregation despite their low ages. We apply these results to a cluster mass function, and find that the peak of the mass distribution would appear to observers shifted to lower masses by as much as 0.2 dex. The star formation rate derived for the cluster population is then underestimated by from 20 to 50 per cent.     

An excursion set model for the distribution of dark matter and dark matter haloes

A model of the gravitationally evolved dark matter distribution, in the Eulerian space, is developed. It is a simple extension of the excursion set model that is commonly used to estimate the mass function of collapsed dark matter haloes. In addition to describing the evolution of the Eulerian space distribution of the haloes, the model allows one to describe the evolution of the dark matter itself. It can also be used to describe density profiles, on scales larger than the virial radius of these haloes, and to quantify the way in which matter flows in and out of Eulerian cells. When the initial Lagrangian space distribution is white noise Gaussian, the model suggests that the Inverse Gaussian distribution should provide a reasonably good approximation to the evolved Eulerian density field, in agreement with numerical simulations. Application of this model to clustering from more general Gaussian initial conditions is discussed at the end.     

Low-redshift Galaxies in the Infrared and Submillimetre

In recent years the evolution of dust-enshrouded galaxies has become an important issue in cosmology. We present the results of a study linking the high- and low-redshift populations of such galaxies. Previously, a simple hierarchical clustering model was found to be successful in describing the high-redshift populations of dusty galaxies. This study showed that strong evolution could not be avoided in order to satisfy all the high-redshift infrared and submillimetre constraints. We now apply the same model to the low-redshift data available, and find that the strong evolution inferred is compatible with the count and redshift data we have from IRAS and ISO. This revised version was published online in September 2006 with corrections to the Cover Date.     

Smoothing supernova data to reconstruct the expansion history of the Universe and its age

We propose a non-parametric method of smoothing supernova data over redshift using a Gaussian kernel in order to reconstruct important cosmological quantities including   H ( z )  and   w ( z )  in a model-independent manner. This method is shown to be successful in discriminating between different models of dark energy when the quality of data is commensurate with that expected from the future Supernova Acceleration Probe ( SNAP ). We find that the Hubble parameter is especially well determined and useful for this purpose. The look-back time of the Universe may also be determined to a very high degree of accuracy (≲0.2 per cent) using this method. By refining the method, it is also possible to obtain reasonable bounds on the equation of state of dark energy. We explore a new diagnostic of dark energy – the ' w -probe'– which can be calculated from the first derivative of the data. We find that this diagnostic is reconstructed extremely accurately for different reconstruction methods even if Ω_{0 m} is marginalized over. The w -probe can be used to successfully distinguish between Λ cold dark matter and other models of dark energy to a high degree of accuracy.     

The evolution of the core and surface magnetic fields in isolated neutron stars

We apply the model of flux expulsion from the superfluid and superconductive core of a neutron star, developed by Konenkov & Geppert, both to neutron star models based on different equations of state and to different initial magnetic field structures. Initially, when the core and the surface magnetic field are of the same order of magnitude, the rate of flux expulsion from the core is almost independent of the equation of state, and the evolution of the surface field decouples from the core field evolution with increasing stiffness. When the surface field is initially much stronger than the core field, the magnetic and rotational evolution resembles that of a neutron star with a purely crustal field configuration; the only difference is the occurrence of a residual field. In the case of an initially submerged field, significant differences from the standard evolution only occur during the early period of the life of a neutron star, until the field has been re-diffused to the surface. The reminder of the episode of submergence is a correlation of the residual field strength with the submergence depth of the initial field. We discuss the effect of the re-diffusion of the magnetic field on the difference between the real and the active age of young pulsars and on their braking indices. Finally, we estimate the shear stresses built up by the moving fluxoids at the crust–core interface and show that these stresses may cause crust cracking, preferentially in neutron stars with a soft equation of state.     

On the errors on Ω0: Monte Carlo simulations of the EMSS cluster sample

We perform Monte Carlo simulations of synthetic EMSS cluster samples, to quantify the systematic errors and the statistical uncertainties on the estimate of Ω₀ derived from fits to the cluster number density evolution and to the X-ray temperature distribution up to z =0.83 . We identify the scatter around the relation between cluster X-ray luminosity and temperature to be a source of systematic error, of the order of Δ_systΩ₀=0.09 , if not properly taken into account in the modelling. After correcting for this bias, our best Ω₀ is 0.66. The uncertainties on the shape and normalization of the power spectrum of matter fluctuations imply relatively large uncertainties on this estimate of Ω₀, of the order of Δ_statΩ₀=0.1 at the 1 σ level. On the other hand, the statistical uncertainties due to the finite size of the high-redshift sample are twice as small. Therefore, what is needed in order to improve the accuracy of Ω₀ estimates based on cluster number density evolution is a more reliable measure of the local temperature function and a better understanding of the cluster observed properties both in the local Universe and at high redshift, that is the relation between cluster mass, temperature and luminosity. This requires detailed observations of X-ray selected cluster samples, in comparison with hydrodynamic simulations including refined physics.     

asmooth: a simple and efficient algorithm for adaptive kernel smoothing of two-dimensional imaging data

An efficient algorithm for adaptive kernel smoothing (AKS) of two-dimensional imaging data has been developed and implemented using the Interactive Data Language ( idl ). The functional form of the kernel can be varied (top-hat, Gaussian, etc.) to allow different weighting of the event counts registered within the smoothing region. For each individual pixel, the algorithm increases the smoothing scale until the signal-to-noise ratio (S/N) within the kernel reaches a pre-set value. Thus, noise is suppressed very efficiently, while at the same time real structure, that is, signal that is locally significant at the selected S/N level, is preserved on all scales. In particular, extended features in noise-dominated regions are visually enhanced. The asmooth algorithm differs from other AKS routines in that it allows a quantitative assessment of the goodness of the local signal estimation by producing adaptively smoothed images in which all pixel values share the same S/N above the background .
We apply asmooth to both real observational data (an X-ray image of clusters of galaxies obtained with the Chandra X-ray Observatory) and to a simulated data set. We find the asmooth ed images to be fair representations of the input data in the sense that the residuals are consistent with pure noise, that is, they possess Poissonian variance and a near-Gaussian distribution around a mean of zero, and are spatially uncorrelated.     

A DWT Power Spectrum Analysis of PSCz Galaxies

The power spectrum estimator based on the Discrete Wavelet Transformation (DWT) is applied to detect the clustering power in the IRAS Point Source Catalog Redshift Survey(PSCz). Comparison with mock samples extracted from N-body simulation shows that the DWT power spectrum estimator could provide a robust measurement of banded fluctuation power over a range of wavenumbers 0.1-2.0h Mpc^-1. We have fitted three typical CDM models (SCDM, TCDM and ACDM) using the Peacock-Dodds formula including non-linear evolution and redshift distortion. We find that, our results are in good agreement with other statistical measurements of the PSCz.     

Unravelling the mystery of the M31 bar

The inclination of M31 is too close to edge-on for a bar component to be easily recognized and is not sufficiently edge-on for a boxy/peanut bulge to protrude clearly out of the equatorial plane. Nevertheless, a sufficient number of clues allow us to argue that this galaxy is barred. We use fully self-consistent N -body simulations of barred galaxies and compare them with both photometric and kinematic observational data for M31. In particular, we rely on the near-infrared photometry presented in a companion paper. We compare isodensity contours to isophotal contours and the light profile along cuts parallel to the galaxy major axis and offset towards the north, or the south, to mass profiles along similar cuts on the model. All these comparisons, as well as position–velocity diagrams for the gaseous component, give us strong arguments that M31 is barred. We compare four fiducial N -body models to the data and thus set constraints on the parameters of the M31 bar, as its strength, length and orientation. Our 'best' models, although not meant to be exact models of M31, reproduce in a very satisfactory way the main relevant observations. We present arguments that M31 has both a classical and a boxy/peanut bulge. Its pseudo-ring-like structure at roughly 50 arcmin is near the outer Lindblad resonance of the bar and could thus be an outer ring, as often observed in barred galaxies. The shape of the isophotes also argues that the vertically thin part of the M31 bar extends considerably further out than its boxy bulge, that is, that the boxy bulge is only part of the bar, thus confirming predictions from orbital structure studies and from previous N -body simulations. It seems very likely that the backbone of M31's boxy bulge is families of periodic orbits, members of the x₁-tree and bifurcating from the x₁ family at its higher order vertical resonances, such as the x1v3 or x1v4 families.     

PSCz superclusters: detection, shapes and cosmological implications

We study the possibility of correctly identifying, from the smooth galaxy density field of the PSC z flux-limited catalogue, high-density regions (superclusters) and recovering their true shapes in the presence of a bias introduced by the coupling between the selection function and the constant radius smoothing. We quantify such systematic biases in the smoothed PSC z density field and after applying the necessary corrections we study supercluster multiplicity and morphologies using a differential geometry definition of shape. Our results strongly suggest that filamentary morphology is the dominant feature of PSC z superclusters. Finally, we compare our results with those expected in three different cosmological models and find that the Λ cold dark matter (CDM) model (Ω_Λ=1−Ω_m=0.7) performs better than Ω_m=1 CDM models.     

Large-scale magnetic field of the G8 dwarf ξ Bootis A

We investigate the magnetic geometry of the active G8 dwarf ξ Bootis A (ξ Boo A), from spectropolarimetric observations obtained in 2003 with the MuSiCoS échelle spectropolarimeter at the Télescope Bernard Lyot (Observatoire du Pic du Midi, France). We repeatedly detect a photospheric magnetic field, with periodic variations consistent with rotational modulation. Circularly polarized (Stokes V) line profiles present a systematic asymmetry, showing up as an excess in amplitude and area of the blue lobe of the profiles. Direct modelling of Stokes V profiles suggests that the global magnetic field is composed of two main components, with an inclined dipole and a large-scale toroidal field. We derive a dipole intensity of about 40 G, with an inclination of 35° of the dipole with respect to the rotation axis. The toroidal field strength is of the order of 120 G. A noticeable evolution of the field geometry is observed over the 40 nights of our observation window and results in an increase in field strength and dipole inclination.
This study is the first step of a long-term monitoring of ξ Boo A and other active solar-type stars, with the aim of investigating secular fluctuations of stellar magnetic geometries induced by activity cycles.