BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

M. P. Hobson

M. P. Hobson appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
astro-ph.COastro-phastro-ph.IMgr-qcInformation Theorymath.ITastro-ph.EPhep-phphysics.comp-phastro-ph.GAhep-exmath.NANumerical Analysisphysics.data-an

Catalog footprint

What is connected

47works
14topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

47 published item(s)

preprint2022arXiv

COSMOPOWER: emulating cosmological power spectra for accelerated Bayesian inference from next-generation surveys

We present $\it{CosmoPower}$, a suite of neural cosmological power spectrum emulators providing orders-of-magnitude acceleration for parameter estimation from two-point statistics analyses of Large-Scale Structure (LSS) and Cosmic Microwave Background (CMB) surveys. The emulators replace the computation of matter and CMB power spectra from Boltzmann codes; thus, they do not need to be re-trained for different choices of astrophysical nuisance parameters or redshift distributions. The matter power spectrum emulation error is less than $0.4\%$ in the wavenumber range $k \in [10^{-5}, 10] \, \mathrm{Mpc}^{-1}$, for redshift $z \in [0, 5]$. $\it{CosmoPower}$ emulates CMB temperature, polarisation and lensing potential power spectra in the $5σ$ region of parameter space around the $\it{Planck}$ best fit values with an error $\lesssim 10\%$ of the expected shot noise for the forthcoming Simons Observatory. $\it{CosmoPower}$ is showcased on a joint cosmic shear and galaxy clustering analysis from the Kilo-Degree Survey, as well as on a Stage IV $\it{Euclid}$-like simulated cosmic shear analysis. For the CMB case, $\it{CosmoPower}$ is tested on a $\it{Planck}$ 2018 CMB temperature and polarisation analysis. The emulators always recover the fiducial cosmological constraints with differences in the posteriors smaller than sampling noise, while providing a speed-up factor up to $O(10^4)$ to the complete inference pipeline. This acceleration allows posterior distributions to be recovered in just a few seconds, as we demonstrate in the $\it{Planck}$ likelihood case. $\it{CosmoPower}$ is written entirely in Python, can be interfaced with all commonly used cosmological samplers and is publicly available at https://github.com/alessiospuriomancini/cosmopower .

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Addressing $H_0$ tension with emergent dark radiation in unitary gravity

We propose a one-parameter extension to $Λ$CDM, expected to strongly affect cosmological tensions. An effective dark radiation component in the early universe redshifts away as hot dark matter, then quintessence, tracking the dominant equation-of-state parameter and leaving a falsifiable torsion field in the current epoch. This picture results from a new Poincaré gauge theory (PGT), one of the most promising among the latest batch of 58 PGTs found to be both power-counting renormalisable and free from ghosts and tachyons. We systematically categorise the cosmologies of 33 of these PGTs, as special cases of the most general parity-preserving, Ostrogradsky-stable PGT with a purely Yang-Mills action. The theory we consider contains two propagating massless gravitons, which may be $J^P=2^+$ (long-range gravitation and gravitational waves). A conspiracy among the coupling constants eliminates the spatial curvature $k\in\{\pm 1,0\}$ from the field equations. We show that this `$k$-screening' is not restricted to conformal gravity theories. The flat Friedmann equations are then emergent, with potentially tension-resolving freedom at the early scale-invariant epoch that reliably gives way to an attractor-like state of modern $Λ$CDM evolution. We compare with related theories and promising special cases, such as $k$-screened theories with negative-definite effective $k$, and more traditional theories with effective $Λ$ and a $J^P=0^-$ massive graviton (dark matter candidate). As a bonus, we analyse similarly constrained actions in the new extended Weyl gauge theory (eWGT). We show that in cosmology, PGT and eWGT span exactly the same classical phenomenology up to a linear map between their coupling constants, hinting at a deeper relationship between the two.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Dense output for highly oscillatory numerical solutions

We present a method to construct a continuous extension (otherwise known as dense output) for a numerical routine in the special case of the numerical solution being a scalar-valued function exhibiting rapid oscillations. Such cases call for numerical routines that make use of the known global behaviour of the solution, one example being methods using asymptotic expansions to forecast the solution at each step of the independent variable. An example is oscode, numerical routine which uses the Wentzel-Kramers-Brillouin (WKB) approximation when the solution oscillates rapidly and otherwise behaves as a Runge-Kutta (RK) solver. Polynomial interpolation is not suitable for producing the solution at an arbitrary point mid-step, since efficient numerical methods based on the WKB approximation will step through multiple oscillations in a single step. Instead we construct the continuous solution by extending the numerical quadrature used in computing a WKB approximation of the solution with no additional evaluations of the differential equation or terms within, and provide an error estimate on this dense output. Finally, we draw attention to previous work on the continuous extension of Runge-Kutta formulae, and construct an extension to a RK method based on Gauss--Lobatto quadrature nodes, thus describing how to generate dense output from each of the methods underlying oscode.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Quantum initial conditions for inflation and canonical invariance

We investigate the transformation of initial conditions for primordial curvature perturbations under two types of transformations of the associated action: simultaneous redefinition of time and the field to be quantised, and the addition of surface terms. The latter encompasses all canonical transformations, whilst the time- and field-redefinition is a distinct, non-canonical transformation since the initial and destination systems use different times. Actions related to each other via such transformations yield identical equations of motion and preserve the commutator structure. They further preserve the time-evolution of expectation values of quantum operators unless the vacuum state also changes under the transformation. These properties suggest that it is of interest to investigate vacuum prescriptions that also remain unchanged under canonical transformations. We find that initial conditions derived via minimising the vacuum expectation value of the Hamiltonian and those obtained using the Danielsson vacuum prescription are not invariant under these transformations, whereas those obtained by minimising the local energy density are. We derive the range of physically distinct initial conditions obtainable by Hamiltonian diagonalisation, and illustrate their effect on the scalar primordial power spectrum and the Cosmic Microwave Background under the just enough inflation model. We also generalise the analogy between the dynamics of a quantum scalar field on a curved, time-dependent spacetime and the gauge-invariant curvature perturbation. We argue that the invariance of the vacuum prescription obtained by minimising the renormalised stress--energy tensor should make it the preferred procedure for setting initial conditions for primordial perturbations. All other procedures reviewed in this work yield ambiguous initial conditions, which is problematic both in theory and practice.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

An efficient method for solving highly oscillatory ordinary differential equations with applications to physical systems

We present a novel numerical routine (oscode) with a C++ and Python interface for the efficient solution of one-dimensional, second-order, ordinary differential equations with rapidly oscillating solutions. The method is based on a Runge-Kutta-like stepping procedure that makes use of the Wentzel-Kramers-Brillouin (WKB) approximation to skip regions of integration where the characteristic frequency varies slowly. In regions where this is not the case, the method is able to switch to a made-to-measure Runge-Kutta integrator that minimises the total number of function evaluations. We demonstrate the effectiveness of the method with example solutions of the Airy equation and an equation exhibiting a burst of oscillations, discussing the error properties of the method in detail. We then show the method applied to physical systems. First, the one-dimensional, time-independent Schrödinger equation is solved as part of a shooting method to search for the energy eigenvalues for a potential with quartic anharmonicity. Then, the method is used to solve the Mukhanov-Sasaki equation describing the evolution of cosmological perturbations, and the primordial power spectrum of the perturbations is computed in different cosmological scenarios. We compare the performance of our solver in calculating a primordial power spectrum of scalar perturbations to that of BINGO, an efficient code specifically designed for such applications.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Novel quantum initial conditions for inflation

We present a novel approach for setting initial conditions on the mode functions of the Mukhanov Sazaki equation. These conditions are motivated by minimisation of the renormalised stress-energy tensor, and are valid for setting a vacuum state even in a context where the spacetime is changing rapidly. Moreover, these alternative conditions are potentially observationally distinguishable. We apply this to the kinetically dominated universe, and compare with the more traditional approach.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Comparison of Sunyaev-Zel'dovich measurements from Planck and from the Arcminute Microkelvin Imager for 99 galaxy clusters

We present observations and analysis of a sample of 123 galaxy clusters from the 2013 Planck catalogue of Sunyaev-Zel'dovich sources with the Arcminute Microkelvin Imager (AMI), a ground-based radio interferometer. AMI provides an independent measurement with higher angular resolution, 3 arcmin compared to the Planck beams of 5-10 arcmin. The AMI observations thus provide validation of the cluster detections, improved positional estimates, and a consistency check on the fitted 'size' ($θ_{s}$) and 'flux' ($Y_{\rm tot}$) parameters in the Generalised Navarro, Frenk and White (GNFW) model. We detect 99 of the clusters. We use the AMI positional estimates to check the positional estimates and error-bars produced by the Planck algorithms PowellSnakes and MMF3. We find that $Y_{\rm tot}$ values as measured by AMI are biased downwards with respect to the Planck constraints, especially for high Planck-SNR clusters. We perform simulations to show that this can be explained by deviation from the 'universal' pressure profile shape used to model the clusters. We show that AMI data can constrain the $α$ and $β$ parameters describing the shape of the profile in the GNFW model for individual clusters provided careful attention is paid to the degeneracies between parameters, but one requires information on a wider range of angular scales than are present in AMI data alone to correctly constrain all parameters simultaneously.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Generative pulsar timing analysis

A new Bayesian method for the analysis of folded pulsar timing data is presented that allows for the simultaneous evaluation of evolution in the pulse profile in either frequency or time, along with the timing model and additional stochastic processes such as red spin noise, or dispersion measure variations. We model the pulse profiles using `shapelets' - a complete ortho-normal set of basis functions that allow us to recreate any physical profile shape. Any evolution in the profiles can then be described as either an arbitrary number of independent profiles, or using some functional form. We perform simulations to compare this approach with established methods for pulsar timing analysis, and to demonstrate model selection between different evolutionary scenarios using the Bayesian evidence. %s The simplicity of our method allows for many possible extensions, such as including models for correlated noise in the pulse profile, or broadening of the pulse profiles due to scattering. As such, while it is a marked departure from standard pulsar timing analysis methods, it has clear applications for both new and current datasets, such as those from the European Pulsar Timing Array (EPTA) and International Pulsar Timing Array (IPTA).

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

PolyChord: nested sampling for cosmology

PolyChord is a novel nested sampling algorithm tailored for high dimensional parameter spaces. In addition, it can fully exploit a hierarchy of parameter speeds such as is found in CosmoMC and CAMB. It utilises slice sampling at each iteration to sample within the hard likelihood constraint of nested sampling. It can identify and evolve separate modes of a posterior semi-independently and is parallelised using openMPI. PolyChord is available for download at: http://ccpforge.cse.rl.ac.uk/gf/project/polychord/

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

PolyChord: next-generation nested sampling

PolyChord is a novel nested sampling algorithm tailored for high-dimensional parameter spaces. This paper coincides with the release of PolyChord v1.3, and provides an extensive account of the algorithm. PolyChord utilises slice sampling at each iteration to sample within the hard likelihood constraint of nested sampling. It can identify and evolve separate modes of a posterior semi-independently, and is parallelised using openMPI. It is capable of exploiting a hierarchy of parameter speeds such as those present in CosmoMC and CAMB, and is now in use in the CosmoChord and ModeChord codes. PolyChord is available for download at: http://ccpforge.cse.rl.ac.uk/gf/project/polychord/

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Bayes-X: a Bayesian inference tool for the analysis of X-ray observations of galaxy clusters

We present the first public release of our Bayesian inference tool, Bayes-X, for the analysis of X-ray observations of galaxy clusters. We illustrate the use of Bayes-X by analysing a set of four simulated clusters at z=0.2-0.9 as they would be observed by a Chandra-like X-ray observatory. In both the simulations and the analysis pipeline we assume that the dark matter density follows a spherically-symmetric Navarro, Frenk and White (NFW) profile and that the gas pressure is described by a generalised NFW (GNFW) profile. We then perform four sets of analyses. By numerically exploring the joint probability distribution of the cluster parameters given simulated Chandra-like data, we show that the model and analysis technique can robustly return the simulated cluster input quantities, constrain the cluster physical parameters and reveal the degeneracies among the model parameters and cluster physical parameters. We then analyse Chandra data on the nearby cluster, A262, and derive the cluster physical profiles. To illustrate the performance of the Bayesian model selection, we also carried out analyses assuming an Einasto profile for the matter density and calculated the Bayes factor. The results of the model selection analyses for the simulated data favour the NFW model as expected. However, we find that the Einasto profile is preferred in the analysis of A262. The Bayes-X software, which is implemented in Fortran 90, is available at http://www.mrao.cam.ac.uk/facilities/software/bayesx/.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Testing the mutual consistency of different supernovae surveys

It is now common practice to constrain cosmological parameters using supernovae (SNe) catalogues constructed from several different surveys. Before performing such a joint analysis, however, one should check that parameter constraints derived from the individual SNe surveys that make up the catalogue are mutually consistent. We describe a statistically-robust mutual consistency test, which we calibrate using simulations, and apply it to each pairwise combination of the surveys making up, respectively, the UNION2 catalogue and the very recent JLA compilation by Betoule et al. We find no inconsistencies in the latter case, but conclusive evidence for inconsistency between some survey pairs in the UNION2 catalogue.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Bayesian constraints on dark matter halo properties using gravitationally-lensed supernovae

A hierarchical Bayesian method is applied to the analysis of Type-Ia supernovae (SNIa) observations to constrain the properties of the dark matter haloes of galaxies along the SNIa lines-of-sight via their gravitational lensing effect. The full joint posterior distribution of the dark matter halo parameters is explored using the nested sampling algorithm {\sc MultiNest}, which also efficiently calculates the Bayesian evidence, thereby facilitating robust model comparison. We first demonstrate the capabilities of the method by applying it to realistic simulated SNIa data, based on the real 3-year data release from the Supernova Legacy Survey (SNLS3). Assuming typical values for the halo parameters in our simulations, we find that a catalogue analogous to the existing SNLS3 data set is incapable of detecting the lensing signal, but a catalogue containing approximately three times as many SNIa does produce robust and accurate parameter constraints and model selection results for two halo models: a truncated singular isothermal sphere (SIS) and a Navarro--Frenk--White (NFW) profile, thereby validating our analysis methodology. In the analysis of the real SNLS3 data, contrary to previous studies, we obtain only a very marginal detection of a lensing signal and weak constraints on the halo parameters for the truncated SIS model, although these constraints are tighter than those obtained from the equivalent simulated SNIa data set. This difference is driven by a preferred value of $η\approx 1$ in the assumed scaling-law $σ\propto L^η$ between velocity dispersion and luminosity, which is somewhat higher than the canonical values of $η= \tfrac{1}{4}$ and $η= \tfrac{1}{3}$ for early and late-type galaxies, respectively, and leads to a stronger lensing effect by the halo. No detection of a lensing signal is made for the NFW model.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Constraints on the Tensor-to-Scalar ratio for non-power-law models

Recent cosmological observations hint at a deviation from the simple power-law form of the primordial spectrum of curvature perturbations. In this paper we show that in the presence of a tensor component, a turn-over in the initial spectrum is preferred by current observations, and hence non-power-law models ought to be considered. For instance, for a power-law parameterisation with both a tensor component and running parameter, current data show a preference for a negative running at more than $2.5σ$ C.L. As a consequence of this deviation from a power-law, constraints on the tensor-to-scalar ratio $r$ are slightly broader. We also present constraints on the inflationary parameters for a model-independent reconstruction and the Lasenby & Doran (LD) model. In particular, the constraints on the tensor-to-scalar ratio from the LD model are: $r_{\rm LD}=0.11\pm{0.024}$. In addition to current data, we show expected constraints from Planck-like and CMB-Pol sensitivity experiments by using Markov-Chain-Monte-Carlo sampling chains. For all the models, we have included the Bayesian Evidence to perform a model selection analysis. The Bayes factor, using current observations, shows a strong preference for the LD model over the standard power-law parameterisation, and provides an insight into the accuracy of differentiating models through future surveys.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Exploring local fNL estimators based on the binned bispectrum

We explore different estimators of the local non-linear coupling parameter, fNL, based on the binned bispectrum presented in Bucher et al. Using simulations of Wilkinson Microwave Anisotropy Probe (WMAP)-7yr data, we compare the performance of a regression neural network with a χ^2-minimization and study the dependence of the results on the presence of the linear term in the analysis and on the use of inpainting for masked regions. Both methods obtain similar results and are robust to the use of inpainting, but the neural network estimator converges considerably faster. We also examine the performance of a simplified χ^2 estimator that assumes a diagonal matrix and has the linear term subtracted, which considerably reduces the computational time; in this case inpainting is found to be crucial. The estimators are also applied to real WMAP-7yr data, yielding constraints at 95% confidence level of -3< fNL <83.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

A Bayesian study of the primordial power spectrum from a novel closed universe model

We constrain the shape of the primordial power spectrum using recent measurements of the cosmic microwave background (CMB) from the Wilkinson Microwave Anisotropy Probe (WMAP) 7-year data and other high-resolution CMB experiments. We also include observations of the matter power spectrum from the luminous red galaxy (LRG) subset DR7 of the Sloan Digital Sky Survey (SDSS). We consider two different models of the primordial power spectrum. The first is the standard nearly scale-invariant spectrum in the form of a generalised power-law parameterised in terms of the spectral amplitude $A_{\rm s}$, the spectral index $n_{\rm s}$ and (possibly) the running parameter $n_{\rm run}$. The second spectrum is derived from the Lasenby and Doran (LD) model. The LD model is based on the restriction of the total conformal time available in a closed Universe and the predicted primordial power spectrum depends upon just two parameters. An important feature of the LD spectrum is that it naturally incorporates an exponential fall-off on large scales, which might provide a possible explanation for the lower-than-expected power observed at low multipoles in the CMB. In addition to parameter estimation, we compare both models using Bayesian model selection. We find there is a significant preference for the LD model over a simple power-law spectrum for a CMB-only dataset, and over models with an equal number of parameters for all the datasets considered.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

A blind detection of a large, complex, Sunyaev--Zel'dovich structure

We present an interesting Sunyaev-Zel'dovich (SZ) detection in the first of the Arcminute Microkelvin Imager (AMI) 'blind', degree-square fields to have been observed down to our target sensitivity of 100μJy/beam. In follow-up deep pointed observations the SZ effect is detected with a maximum peak decrement greater than 8 \times the thermal noise. No corresponding emission is visible in the ROSAT all-sky X-ray survey and no cluster is evident in the Palomar all-sky optical survey. Compared with existing SZ images of distant clusters, the extent is large (\approx 10') and complex; our analysis favours a model containing two clusters rather than a single cluster. Our Bayesian analysis is currently limited to modelling each cluster with an ellipsoidal or spherical beta-model, which do not do justice to this decrement. Fitting an ellipsoid to the deeper candidate we find the following. (a) Assuming that the Evrard et al. (2002) approximation to Press & Schechter (1974) correctly gives the number density of clusters as a function of mass and redshift, then, in the search area, the formal Bayesian probability ratio of the AMI detection of this cluster is 7.9 \times 10^4:1; alternatively assuming Jenkins et al. (2001) as the true prior, the formal Bayesian probability ratio of detection is 2.1 \times 10^5:1. (b) The cluster mass is MT,200 = 5.5+1.2\times 10^14h-1M\odot. (c) Abandoning a physical model with num- -1.3 70 ber density prior and instead simply modelling the SZ decrement using a phenomenological β-model of temperature decrement as a function of angular distance, we find a central SZ temperature decrement of -295+36 μK - this allows for CMB primary anisotropies, receiver -15 noise and radio sources. We are unsure if the cluster system we observe is a merging system or two separate clusters.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

A simple and robust method for automated photometric classification of supernovae using neural networks

A method is presented for automated photometric classification of supernovae (SNe) as Type-Ia or non-Ia. A two-step approach is adopted in which: (i) the SN lightcurve flux measurements in each observing filter are fitted separately; and (ii) the fitted function parameters and their associated uncertainties, along with the number of flux measurements, the maximum-likelihood value of the fit and Bayesian evidence for the model, are used as the input feature vector to a classification neural network (NN) that outputs the probability that the SN under consideration is of Type-Ia. The method is trained and tested using data released following the SuperNova Photometric Classification Challenge (SNPCC). We consider several random divisions of the data into training and testing sets: for instance, for our sample D_1 (D_4), a total of 10% (40%) of the data are involved in training the algorithm and the remainder used for blind testing of the resulting classifier; we make no selection cuts. Assigning a canonical threshold probability of p_th=0.5 on the NN output to classify a SN as Type-Ia, for the sample D_1 (D_4) we obtain a completeness of 0.78 (0.82), purity of 0.77 (0.82), and SNPCC figure-of-merit of 0.41 (0.50). Including the SN host-galaxy redshift and its uncertainty as additional inputs to the NN results in a modest 5-10% increase in these values. We find that the classification quality does not vary significantly with SN redshift. Moreover, our probabilistic classification method allows one to calculate the expected completeness, purity or other measures of classification quality as a function of the p_th, without knowing the true classes of the SNe in the testing sample. The method may thus be improved further by optimising p_th and can easily be extended to divide non-Ia SNe into their different classes.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

A simple parametric model for spherical galaxy clusters

We present an analytic parametric model to describe the baryonic and dark matter distributions in clusters of galaxies with spherical symmetry. It is assumed that the dark matter density follows a Navarro, Frenk and White (NFW) profile and that the gas pressure is described by a generalised NFW (GNFW) profile. By further demanding hydrostatic equilibrium and that the gas fraction is small throughout the cluster, one obtains unique functional forms, dependent on basic cluster parameters, for the radial profiles of all the properties of interest in the cluster. We show these profiles are consistent both with numerical simulations and multi-wavelength observations of clusters. We also use our model to analyse six simulated SZ clusters as well as A611 SZ data from the Arcminute Microkelvin Imager (AMI). In each case, we derive the radial profile of the enclosed total mass and the gas pressure and show that the results are in good agreement with our model prediction.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Comparison of cosmological parameter inference methods applied to supernovae lightcurves fitted with SALT2

We present a comparison of two methods for cosmological parameter inference from supernovae Ia lightcurves fitted with the SALT2 technique. The standard chi-square methodology and the recently proposed Bayesian hierarchical method (BHM) are each applied to identical sets of simulations based on the 3-year data release from the Supernova Legacy Survey (SNLS3), and also data from the Sloan Digital Sky Survey (SDSS), the Low Redshift sample and the Hubble Space Telescope (HST), assuming a concordance LCDM cosmology. For both methods, we find that the recovered values of the cosmological parameters, and the global nuisance parameters controlling the stretch and colour corrections to the supernovae lightcurves, suffer from small biasses. The magnitude of the biasses is similar in both cases, with the BHM yielding slightly more accurate results, in particular for cosmological parameters when applied to just the SNLS3 single survey data sets. Most notably, in this case, the biasses in the recovered matter density $Ω_{\rm m,0}$ are in opposite directions for the two methods. For any given realisation of the SNLS3-type data, this can result in a $\sim 2 σ$ discrepancy in the estimated value of $Ω_{\rm m,0}$ between the two methods, which we find to be the case for real SNLS3 data. As more higher and lower redshift SNIa samples are included, however, the cosmological parameter estimates of the two methods converge.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Constraints on fNL from Wilkinson Microwave Anisotropy Probe 7-year data using a neural network classifier

We present a multi-class neural network (NN) classifier as a method to measure nonGaussianity, characterised by the local non-linear coupling parameter fNL, in maps of the cosmic microwave background (CMB) radiation. The classifier is trained on simulated non-Gaussian CMB maps with a range of known fNL values by providing it with wavelet coefficients of the maps; we consider both the HealPix (HW) wavelet and the spherical Mexican hat wavelet (SMHW). When applied to simulated test maps, the NN classfier produces results in very good agreement with those obtained using standard chi2 minimization. The standard deviations of the fNL estimates for WMAPlike simulations were σ = 22 and σ = 33 for the SMHW and the HW, respectively, which are extremely close to those obtained using classical statistical methods in Curto et al. and Casaponsa et al. Moreover, the NN classifier does not require the inversion of a large covariance matrix, thus avoiding any need to regularise the matrix when it is not directly invertible, and is considerably faster.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Mass and pressure constraints on galaxy clusters from interferometric SZ observations

Following on our previous study of an analytic parametric model to describe the baryonic and dark matter distributions in clusters of galaxies with spherical symmetry, we perform an SZ analysis of a set of simulated clusters and present their mass and pressure profiles. The simulated clusters span a wide range in mass, 2.0 x 10^14 Msun < M200 < 1.0 x 10^15Msun, and observations with the Arcminute Microkelvin Imager (AMI) are simulated through their Sunyaev- Zel'dovich (SZ) effect. We assume that the dark matter density follows a Navarro, Frenk and White (NFW) profile and that the gas pressure is described by a generalised NFW (GNFW) profile. By numerically exploring the probability distributions of the cluster parameters given simulated interferometric SZ data in the context of Bayesian methods, we investigate the capability of this model and analysis technique to return the simulated clusters input quantities. We show that considering the mass and redshift dependency of the cluster halo concentration parameter is crucial in obtaining an unbiased cluster mass estimate and hence deriving the radial profiles of the enclosed total mass and the gas pressure out to r200.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Model selection applied to reconstruction of the Primordial Power Spectrum

The preferred shape for the primordial spectrum of curvature perturbations is determined by performing a Bayesian model selection analysis of cosmological observations. We first reconstruct the spectrum modelled as piecewise linear in \log k between nodes in k-space whose amplitudes and positions are allowed to vary. The number of nodes together with their positions are chosen by the Bayesian evidence, so that we can both determine the complexity supported by the data and locate any features present in the spectrum. In addition to the node-based reconstruction, we consider a set of parameterised models for the primordial spectrum: the standard power-law parameterisation, the spectrum produced from the Lasenby & Doran (LD) model and a simple variant parameterisation. By comparing the Bayesian evidence for different classes of spectra, we find the power-law parameterisation is significantly disfavoured by current cosmological observations, which show a preference for the LD model.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Reconstruction of the Dark Energy equation of state

One of the main challenges of modern cosmology is to investigate the nature of dark energy in our Universe. The properties of such a component are normally summarised as a perfect fluid with a (potentially) time-dependent equation-of-state parameter $w(z)$. We investigate the evolution of this parameter with redshift by performing a Bayesian analysis of current cosmological observations. We model the temporal evolution as piecewise linear in redshift between `nodes', whose $w$-values and redshifts are allowed to vary. The optimal number of nodes is chosen by the Bayesian evidence. In this way, we can both determine the complexity supported by current data and locate any features present in $w(z)$. We compare this node-based reconstruction with some previously well-studied parameterisations: the Chevallier-Polarski-Linder (CPL), the Jassal-Bagla-Padmanabhan (JBP) and the Felice-Nesseris-Tsujikawa (FNT). By comparing the Bayesian evidence for all of these models we find an indication towards possible time-dependence in the dark energy equation-of-state. It is also worth noting that the CPL and JBP models are strongly disfavoured, whilst the FNT is just significantly disfavoured, when compared to a simple cosmological constant $w=-1$. We find that our node-based reconstruction model is slightly disfavoured with respect to the $Λ$CDM model.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Bayesian evidence for two companions orbiting HIP 5158

We present results of a Bayesian analysis of radial velocity (RV) data for the star HIP 5158, confirming the presence of two companions and also constraining their orbital parameters. Assuming Keplerian orbits, the two-companion model is found to be e^{48} times more probable than the one-planet model, although the orbital parameters of the second companion are only weakly constrained. The derived orbital periods are 345.6 +/- 2.0 d and 9017.8 +/- 3180.7 d respectively, and the corresponding eccentricities are 0.54 +/- 0.04 and 0.14 +/- 0.10. The limits on planetary mass (m \sin i) and semimajor axis are (1.44 +/- 0.14 M_{J}, 0.89 +/- 0.01 AU) and (15.04 +/- 10.55 M_{J}, 7.70 +/- 1.88 AU) respectively. Owing to large uncertainty on the mass of the second companion, we are unable to determine whether it is a planet or a brown dwarf. The remaining `noise' (stellar jitter) unaccounted for by the model is 2.28 +/- 0.31 m/s. We also analysed a three-companion model, but found it to be e^{8} times less probable than the two-companion model.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Constraints on general primordial non-Gaussianity using wavelets for the Wilkinson Microwave anisotropy probe 7-year data

We present constraints on the non-linear coupling parameter fnl with the Wilkinson Microwave Anisotropy Probe (WMAP) data. We use the method based on the spherical Mexican hat wavelet (SMHW) to measure the fnl parameter for three of the most interesting shapes of primordial non-Gaussianity: local, equilateral and orthogonal. Our results indicate that this parameter is compatible with a Gaussian distribution within the two sigma confidence level (CL) for the three shapes and the results are consistent with the values presented by the WMAP team. We have included in our analysis the impact on fnl due to contamination by unresolved point sources. The point sources add a positive contribution of Delta(fnl) = 2.5 \pm 3.0, Delta(fnl) = 37 \pm 18 and Delta(fnl) = 25 \pm 14 for the local, equilateral and orthogonal cases respectively. As mentioned by the WMAP team, the contribution of the point sources to the orthogonal and equilateral form is expected to be larger than the local one and thus it cannot be neglected in future constraints on these parameters. Taking into account this contamination, our best estimates for fnl are -16.0 \leq fnl \leq 76.0, -382 \leq fnl \leq 202 and -394 \leq fnl \leq 34 at 95% CL for the local, equilateral and orthogonal cases respectively. The three shapes are compatible with zero at 95% CL (2σ). Our conclusion is that the WMAP 7-year data are consistent with Gaussian primordial fluctuations within ~2σ CL. We stress however the importance of taking into account the unresolved point sources in the measurement of fnl in future works, especially when using more precise data sets such as the forthcoming Planck data.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Detecting extrasolar planets from stellar radial velocities using Bayesian evidence

Stellar radial velocity (RV) measurements have proven to be a very successful method for detecting extrasolar planets. Analysing RV data to determine the parameters of the extrasolar planets is a significant statistical challenge owing to the presence of multiple planets and various degeneracies between orbital parameters. Determining the number of planets favoured by the observed data is an even more difficult task. Bayesian model selection provides a mathematically rigorous solution to this problem by calculating marginal posterior probabilities of models with different number of planets, but the use of this method in extrasolar planetary searches has been hampered by the computational cost of the evaluating Bayesian evidence. Nonetheless, Bayesian model selection has the potential to improve the interpretation of existing observational data and possibly detect yet undiscovered planets. We present a new and efficient Bayesian method for determining the number of extrasolar planets, as well as for inferring their orbital parameters, without having to calculate directly the Bayesian evidence for models containing a large number of planets. Instead, we work iteratively and at each iteration obtain a conservative lower limit on the odds ratio for the inclusion of an additional planet into the model. We apply this method to simulated data-sets containing one and two planets and successfully recover the correct number of planets and reliable constraints on the orbital parameters. We also apply our method to RV measurements of HD 37124, 47 Ursae Majoris and HD 10180. For HD 37124, we confirm that the current data strongly favour a three-planet system. We find strong evidence for the presence of a fourth planet in 47 Ursae Majoris, but its orbital period is suspiciously close to one year, casting doubt on its validity. For HD 10180 we find strong evidence for a six-planet system.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Powellsnakes II: a fast Bayesian approach to discrete object detection in multi-frequency astronomical data sets

Powellsnakes is a Bayesian algorithm for detecting compact objects embedded in a diffuse background, and was selected and successfully employed by the Planck consortium in the production of its first public deliverable: the Early Release Compact Source Catalogue (ERCSC). We present the critical foundations and main directions of further development of PwS, which extend it in terms of formal correctness and the optimal use of all the available information in a consistent unified framework, where no distinction is made between point sources (unresolved objects), SZ clusters, single or multi-channel detection. An emphasis is placed on the necessity of a multi-frequency, multi-model detection algorithm in order to achieve optimality.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Weak lensing by triaxial galaxy clusters

Weak gravitational lensing studies of galaxy clusters often assume a spherical cluster model to simplify the analysis, but some recent studies have suggested this simplifying assumption may result in large biases in estimated cluster masses and concentration values, since clusters are expected to exhibit triaxiality. Several such analyses have, however, quoted expressions for the spatial derivatives of the lensing potential in triaxial models, which are open to misinterpretation. In this paper, we give a clear description of weak lensing by triaxial NFW galaxy clusters and also present an efficient and robust method to model these clusters and obtain parameter estimates. By considering four highly triaxial NFW galaxy clusters, we re-examine the impact of simplifying spherical assumptions and found that while the concentration estimates are largely unbiased except in one of our traixial NFW simulated clusters, for which the concentration is only slightly biased, the masses are significantly biased, by up to 40%, for all the clusters we analysed. Moreover, we find that such assumptions can lead to the erroneous conclusion that some substructure is present in the galaxy clusters or, even worse, that multiple galaxy clusters are present in the field. Our cluster fitting method also allows one to answer the question of whether a given cluster exhibits triaxiality or a simple spherical model is good enough.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

A 33 GHz VSA survey of the Galactic plane from 27 to 46 degrees

The Very Small Array (VSA) has been used to survey the l = 27 to 46 deg, |b|<4 deg region of the Galactic plane at a resolution of 13 arcmin. The survey consists of 44 pointings of the VSA, each with a r.m.s. sensitivity of ~90 mJy/beam. These data are combined in a mosaic to produce a map of the area. The majority of the sources within the map are HII regions. We investigated anomalous radio emission from the warm dust in 9 HII regions of the survey by making spectra extending from GHz frequencies to the FIR IRAS frequencies. Acillary radio data at 1.4, 2.7, 4.85, 8.35, 10.55, 14.35 and 94 GHz in addition to the 100, 60, 25 and 12 micron IRAS bands were used to construct the spectra. From each spectrum the free-free, thermal dust and anomalous dust emission were determined for each HII region. The mean ratio of 33 GHz anomalous flux density to FIR 100 micron flux density for the 9 selected HII regions was 1.10 +/-0.21x10^(-4). When combined with 6 HII regions previously observed with the VSA and the CBI, the anomalous emission from warm dust in HII regions is detected with a 33 GHz emissivity of 4.65 +/- 0.4 micro K/ (MJy/sr) at 11.5σ. The anomalous radio emission in HII regions is on average 41+/-10 per cent of the radio continuum at 33 GHz.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Comment on "Bayesian evidence: can we beat MultiNest using traditional MCMC methods", by Rutger van Haasteren (arXiv:0911.2150)

In arXiv:0911.2150, Rutger van Haasteren seeks to criticize the nested sampling algorithm for Bayesian data analysis in general and its MultiNest implementation in particular. He introduces a new method for evidence evaluation based on the idea of Voronoi tessellation and requiring samples from the posterior distribution obtained through MCMC based methods. He compares its accuracy and efficiency with MultiNest, concluding that it outperforms MultiNest in several cases. This comparison is completely unfair since the proposed method can not perform the complete Bayesian data analysis including posterior exploration and evidence evaluation on its own while MultiNest allows one to perform Bayesian data analysis end to end. Furthermore, their criticism of nested sampling (and in turn MultiNest) is based on a few conceptual misunderstandings of the algorithm. Here we seek to set the record straight.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Selecting a Model of Supersymmetry Breaking Mediation

We study the problem of selecting between different mechanisms of supersymmetry breaking in the Minimal Supersymmetric Standard Model using current data. We evaluate the Bayesian evidence of four supersymmetry breaking scenarios: mSUGRA, mGMSB, mAMSB and moduli mediation. The results show a strong dependence on the dark matter assumption. Using the inferred cosmological relic density as an upper bound, minimal anomaly mediation is at least moderately favoured over the CMSSM. Our fits also indicate that evidence for a positive sign of the mu parameter is moderate at best. We present constraints on the anomaly and gauge mediated parameter spaces and some previously unexplored aspects of the dark matter phenomenology of the moduli mediation scenario. We use sparticle searches, indirect observables and dark matter observables in the global fit and quantify robustness with respect to prior choice. We quantify how much information is contained within each constraint.

0 citations0 reviewsNo code linkedOpen access
preprint2008arXiv

An excess of emission in the dark cloud LDN 1111 with the Arcminute Microkelvin Imager

We present observations of the Lynds' dark nebula LDN 1111 made at microwave frequencies between 14.6 and 17.2 GHz with the Arcminute Microkelvin Imager (AMI). We find emission in this frequency band in excess of a thermal free--free spectrum extrapolated from data at 1.4 GHz with matched uv-coverage. This excess is > 15 sigma above the predicted emission. We fit the measured spectrum using the spinning dust model of Drain & Lazarian (1998a) and find the best fitting model parameters agree well with those derived from Scuba data for this object by Visser et al. (2001).

0 citations0 reviewsNo code linkedOpen access
preprint2008arXiv

Bayesian modelling of clusters of galaxies from multi-frequency pointed Sunyaev--Zel'dovich observations

We present a Bayesian approach to modelling galaxy clusters using multi-frequency pointed observations from telescopes that exploit the Sunyaev--Zel'dovich effect. We use the recently developed MultiNest technique (Feroz, Hobson & Bridges, 2008) to explore the high-dimensional parameter spaces and also to calculate the Bayesian evidence. This permits robust parameter estimation as well as model comparison. Tests on simulated Arcminute Microkelvin Imager observations of a cluster, in the presence of primary CMB signal, radio point sources (detected as well as an unresolved background) and receiver noise, show that our algorithm is able to analyse jointly the data from six frequency channels, sample the posterior space of the model and calculate the Bayesian evidence very efficiently on a single processor. We also illustrate the robustness of our detection process by applying it to a field with radio sources and primordial CMB but no cluster, and show that indeed no cluster is identified. The extension of our methodology to the detection and modelling of multiple clusters in multi-frequency SZ survey data will be described in a future work.

0 citations0 reviewsNo code linkedOpen access
preprint2008arXiv

Bayesian optimal reconstruction of the primordial power spectrum

The form of the primordial power spectrum has the potential to differentiate strongly between competing models of perturbation generation in the early universe and so is of considerable importance. The recent release of five years of WMAP observations have confirmed the general picture of the primordial power spectrum as deviating slightly from scale invariance with a spectral tilt parameter of n_s ~ 0.96. Nonetheless, many attempts have been made to isolate further features such as breaks and cutoffs using a variety of methods, some employing more than ~ 10 varying parameters. In this paper we apply the robust technique of Bayesian model selection to reconstruct the optimal degree of structure in the spectrum. We model the spectrum simply and generically as piecewise linear in ln k between `nodes' in k-space whose amplitudes are allowed to vary. The number of nodes and their k-space positions are chosen by the Bayesian evidence so that we can identify both the complexity and location of any detected features. Our optimal reconstruction contains, perhaps, surprisingly few features, the data preferring just three nodes. This reconstruction allows for a degree of scale dependence of the tilt with the `turn-over' scale occuring around k ~ 0.016 Mpc^{-1}. More structure is penalised by the evidence as over-fitting the data, so there is currently little point in attempting reconstructions that are more complex.

0 citations0 reviewsNo code linkedOpen access
preprint2008arXiv

MultiNest: an efficient and robust Bayesian inference tool for cosmology and particle physics

We present further development and the first public release of our multimodal nested sampling algorithm, called MultiNest. This Bayesian inference tool calculates the evidence, with an associated error estimate, and produces posterior samples from distributions that may contain multiple modes and pronounced (curving) degeneracies in high dimensions. The developments presented here lead to further substantial improvements in sampling efficiency and robustness, as compared to the original algorithm presented in Feroz & Hobson (2008), which itself significantly outperformed existing MCMC techniques in a wide range of astrophysical inference problems. The accuracy and economy of the MultiNest algorithm is demonstrated by application to two toy problems and to a cosmological inference problem focussing on the extension of the vanilla $Λ$CDM model to include spatial curvature and a varying equation of state for dark energy. The MultiNest software, which is fully parallelized using MPI and includes an interface to CosmoMC, is available at http://www.mrao.cam.ac.uk/software/multinest/. It will also be released as part of the SuperBayeS package, for the analysis of supersymmetric theories of particle physics, at http://www.superbayes.org

0 citations0 reviewsNo code linkedOpen access
preprint2008arXiv

Optimal filters on the sphere

We derive optimal filters on the sphere in the context of detecting compact objects embedded in a stochastic background process. The matched filter and the scale adaptive filter are derived on the sphere in the most general setting, allowing for directional template profiles and filters. The performance and relative merits of the two optimal filters are discussed. The application of optimal filter theory on the sphere to the detection of compact objects is demonstrated on simulated mock data. A naive detection strategy is adopted, with an initial aim of illustrating the application of the new optimal filters derived on the sphere. Nevertheless, this simple object detection strategy is demonstrated to perform well, even a low signal-to-noise ratio. Code written to compute optimal filters on the sphere (S2FIL), to perform fast directional filtering on the sphere (FastCSWT) and to construct the simulated mock data (COMB) are all made publicly available from http://www.mrao.cam.ac.uk/~jdm57/

0 citations0 reviewsNo code linkedOpen access
preprint2008arXiv

The impact of priors and observables on parameter inferences in the Constrained MSSM

We use a newly released version of the SuperBayeS code to analyze the impact of the choice of priors and the influence of various constraints on the statistical conclusions for the preferred values of the parameters of the Constrained MSSM. We assess the effect in a Bayesian framework and compare it with an alternative likelihood-based measure of a profile likelihood. We employ a new scanning algorithm (MultiNest) which increases the computational efficiency by a factor ~200 with respect to previously used techniques. We demonstrate that the currently available data are not yet sufficiently constraining to allow one to determine the preferred values of CMSSM parameters in a way that is completely independent of the choice of priors and statistical measures. While b->s gamma generally favors large m_0, this is in some contrast with the preference for low values of m_0 and m_1/2 that is almost entirely a consequence of a combination of prior effects and a single constraint coming from the anomalous magnetic moment of the muon, which remains somewhat controversial. Using an information-theoretical measure, we find that the cosmological dark matter abundance determination provides at least 80% of the total constraining power of all available observables. Despite the remaining uncertainties, prospects for direct detection in the CMSSM remain excellent, with the spin-independent neutralino-proton cross section almost guaranteed above sigma_SI ~ 10^{-10} pb, independently of the choice of priors or statistics. Likewise, gluino and lightest Higgs discovery at the LHC remain highly encouraging. While in this work we have used the CMSSM as particle physics model, our formalism and scanning technique can be readily applied to a wider class of models with several free parameters.

0 citations0 reviewsNo code linkedOpen access
preprint2007arXiv

Multimodal nested sampling: an efficient and robust alternative to MCMC methods for astronomical data analysis

In performing a Bayesian analysis of astronomical data, two difficult problems often emerge. First, in estimating the parameters of some model for the data, the resulting posterior distribution may be multimodal or exhibit pronounced (curving) degeneracies, which can cause problems for traditional MCMC sampling methods. Second, in selecting between a set of competing models, calculation of the Bayesian evidence for each model is computationally expensive. The nested sampling method introduced by Skilling (2004), has greatly reduced the computational expense of calculating evidences and also produces posterior inferences as a by-product. This method has been applied successfully in cosmological applications by Mukherjee et al. (2006), but their implementation was efficient only for unimodal distributions without pronounced degeneracies. Shaw et al. (2007), recently introduced a clustered nested sampling method which is significantly more efficient in sampling from multimodal posteriors and also determines the expectation and variance of the final evidence from a single run of the algorithm, hence providing a further increase in efficiency. In this paper, we build on the work of Shaw et al. and present three new methods for sampling and evidence evaluation from distributions that may contain multiple modes and significant degeneracies; we also present an even more efficient technique for estimating the uncertainty on the evaluated evidence. These methods lead to a further substantial improvement in sampling efficiency and robustness, and are applied to toy problems to demonstrate the accuracy and economy of the evidence calculation and parameter estimation. Finally, we discuss the use of these methods in performing Bayesian object detection in astronomical datasets.

0 citations0 reviewsNo code linkedOpen access
preprint2006arXiv

A directional continuous wavelet transform on the sphere

A new construction of a directional continuous wavelet analysis on the sphere is derived herein. We adopt the harmonic scaling idea for the spherical dilation operator recently proposed by Sanz et al. but extend the analysis to a more general directional framework. Directional wavelets are a powerful extension that allow one to also probe oriented structure in the analysed function. Our spherical wavelet methodology has the advantage that all functions and operators are defined directly on the sphere. The construction of wavelets in our framework is demonstrated with an example.

0 citations0 reviewsNo code linkedOpen access
preprint2006arXiv

A high-significance detection of non-Gaussianity in the WMAP 3-year data using directional spherical wavelets

We repeat the directional spherical real Morlet wavelet analysis used to detect non-Gaussianity in the WMAP 1-year data (McEwen et al. 2005a), on the WMAP 3-year data. The non-Gaussian signal previously detected is indeed present in the 3-year data, although the significance of the detection is reduced. Using our most conservative method for constructing significance measures, we find the significance of the detection of non-Gaussianity drops from 98.3+/-0.4% to 94.9+/-0.7%; the significance drops from 99.3+/-0.3% to 97.2+/-0.5% using a method based on the $χ^2$ statistic. The wavelet analysis allows us to localise most likely sources of non-Gaussianity on the sky. We detect very similar localised regions in the WMAP 1-year and 3-year data, although the regions extracted appear more pronounced in the 3-year data. When all localised regions are excluded from the analysis the 3-year data is consistent with Gaussianity.

0 citations0 reviewsNo code linkedOpen access
preprint2006arXiv

Fast directional continuous spherical wavelet transform algorithms

We describe the construction of a spherical wavelet analysis through the inverse stereographic projection of the Euclidean planar wavelet framework, introduced originally by Antoine and Vandergheynst and developed further by Wiaux et al. Fast algorithms for performing the directional continuous wavelet analysis on the unit sphere are presented. The fast directional algorithm, based on the fast spherical convolution algorithm developed by Wandelt and Gorski, provides a saving of O(sqrt(Npix)) over a direct quadrature implementation for Npix pixels on the sphere, and allows one to perform a directional spherical wavelet analysis of a 10^6 pixel map on a personal computer.

0 citations0 reviewsNo code linkedOpen access
preprint2001arXiv

The discriminating power of wavelets to detect non-Gaussianity in the CMB

We investigate the power of wavelet techniques in detecting non-Gaussianity in the cosmic microwave background (CMB). We use the method to discriminate between an inflationary and a cosmic strings model using small simulated patches of the sky. We show the importance of the choice of a good test statistic in order to optimise the discriminating power of the wavelet technique. In particular, we construct the Fisher discriminant function, which combines all the information available in the different wavelet scales. We also compare the performance of different decomposition schemes and wavelet bases. For our case, we find that the Mallat and {\it à trous} algorithms are superior to the 2D-tensor wavelets. Using this technique, the inflationary and strings models are clearly distinguished even in the presence of a superposed Gaussian component with twice the rms amplitude of the original cosmic string map.

0 citations0 reviewsNo code linkedOpen access
preprint2000arXiv

Testing the Gaussianity of the COBE-DMR data with spherical wavelets

We investigate the Gaussianity of the 4-year COBE-DMR data (in HEALPix pixelisation) using an analysis based on spherical Haar wavelets. We use all the pixels lying outside the Galactic cut and compute the skewness, kurtosis and scale-scale correlation spectra for the wavelet coefficients at each scale. We also take into account the sensitivity of the method to the orientation of the input signal. We find a detection of non-Gaussianity at $> 99$ per cent level in just one of our statistics. Taking into account the total number of statistics computed, we estimate that the probability of obtaining such a detection by chance for an underlying Gaussian field is 0.69. Therefore, we conclude that the spherical wavelet technique shows no strong evidence of non-Gaussianity in the COBE-DMR data.

0 citations0 reviewsNo code linkedOpen access
preprint1996arXiv

Ideal scales for weighing the Universe

We investigate the performance of a large class of cosmic microwave background experiments with respect to their ability to measure various cosmological parameters. We pay special attention to the measurement of the total cosmological density, $Ω$. We consider interferometer experiments, all-sky single-dish experiments, and also single-dish experiments with a deep-patch technique. Power spectrum estimates for these experiments are studied, and their induced errors in cosmological parameter estimates evaluated. Given this motivation we find various promising corners in the experiment parameter space surveyed. Low noise all-sky satellite experiments are the expensive option, but they are best suited for dealing with large sets of cosmological parameters. At intermediate noises we find a useful corner in high-resolution deep patch single-dish experiments. Interferometers are limited by sample variance, but provide the best estimates based on the very small angular scales. For all these experiments we present conservative, but still promising estimates of the accuracy of the measurement of $Ω$. In these estimates we consider a variety of possible signals not necessarily in the vicinity of standard cold dark matter.

0 citations0 reviewsNo code linkedOpen access
preprint1996arXiv

Observability of secondary Doppler peaks in the CMBR power spectrum by experiments with small fields

We investigate the effects of finite sky coverage on the spectral resolution $Δ\ell$ in the estimation of the CMBR angular power spectrum $C^{\ell}$. A method is developed for obtaining quasi-independent estimates of the power spectrum, and the cosmic/sample variance of these estimates is calculated. The effect of instrumental noise is also considered for prototype interferometer and single-dish experiments. By proposing a statistic for the detection of secondary (Doppler) peaks in the CMBR power spectrum, we then compute the significance level at which such peaks may be detected for a large range of model CMBR experiments. In particular, we investigate experimental design features required to distinguish between competing cosmological theories, such as cosmic strings and inflation, by establishing whether or not secondary peaks are present in the CMBR power spectrum.

0 citations0 reviewsNo code linkedOpen access