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

A. N. Lasenby

A. N. Lasenby 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.GAgr-qcastro-ph.IMInformation Theorymath.ITastro-ph.HEphysics.comp-phmath.NAnucl-exNumerical Analysisphysics.data-anphysics.ins-det

Catalog footprint

What is connected

35works
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

35 published item(s)

preprint2023arXiv

QUIJOTE scientific results -- IV. A northern sky survey in intensity and polarization at 10-20GHz with the Multi-Frequency Instrument

We present QUIJOTE intensity and polarization maps in four frequency bands centred around 11, 13, 17 and 19GHz, and covering approximately 29000 deg$^2$, including most of the Northern sky region. These maps result from 9000 h of observations taken between May 2013 and June 2018 with the first QUIJOTE instrument (MFI), and have angular resolutions of around $1^\circ$, and sensitivities in polarization within the range 35-40 $μ$K per 1-degree beam, being a factor $\sim 2$-$4$ worse in intensity. We discuss the data processing pipeline employed, and the basic characteristics of the maps in terms of real space statistics and angular power spectra. A number of validation tests have been applied to characterise the accuracy of the calibration and the residual level of systematic effects, finding a conservative overall calibration uncertainty of 5%. We also discuss flux densities for four bright celestial sources (Tau A, Cas A, Cyg A and 3C274) which are often used as calibrators at microwave frequencies. The polarization signal in our maps is dominated by synchrotron emission. The distribution of spectral index values between the 11GHz and WMAP 23GHz map peaks at $β=-3.09$ with a standard deviation of 0.14. The measured BB/EE ratio at scales of $\ell=80$ is $0.26\pm 0.07$ for a Galactic cut $|b|>10^\circ$. We find a positive TE correlation for 11GHz at large angular scales ($\ell \lesssim 50$), while the EB and TB signals are consistent with zero in the multipole range $30 \lesssim \ell \lesssim 150$. The maps discussed in this paper are publicly available.

0 citations0 reviewsNo code linkedOpen access
preprint2023arXiv

QUIJOTE scientific results -- IX. Radio sources in the QUIJOTE-MFI wide survey maps

We present the catalogue of Q-U-I JOint TEnerife (QUIJOTE) Wide Survey radio sources extracted from the maps of the Multi-Frequency Instrument compiled between 2012 and 2018. The catalogue contains 786 sources observed in intensity and polarization, and is divided into two separate sub-catalogues: one containing 47 bright sources previously studied by the \emph{Planck} collaboration and an extended catalogue of 739 sources either selected from the \emph{Planck} Second Catalogue of Compact Sources or found through a blind search carried out with a Mexican Hat 2 wavelet. A significant fraction of the sources in our catalogue (38.7 per cent) are within the $|b| \leq 20^\circ$ region of the Galactic plane. We determine statistical properties for those sources that are likely to be extragalactic. We find that these statistical properties are compatible with currently available models, with a $\sim$1.8 Jy completeness limit at 11 GHz. We provide the polarimetric properties of (38, 33, 31, 23) sources with P detected above the $99.99\%$ significance level at (11, 13, 17, 19) GHz, respectively. Median polarization fractions are in the $2.8$-$4.7$\% range in the 11-19 GHz frequency interval. We do not distinguish between Galactic and extragalactic sources here. The results presented here are consistent with those reported in the literature for flat- and steep-spectrum radio sources.

0 citations0 reviewsNo code linkedOpen access
preprint2023arXiv

QUIJOTE scientific results -- V. The microwave intensity and polarisation spectra of the Galactic regions W49, W51 and IC443

We present new intensity and polarisation maps obtained with the QUIJOTE experiment towards the Galactic regions W49, W51 and IC443, covering the frequency range from 10 to 20 GHz at $\sim$ 1 deg angular resolution, with a sensitivity in the range 35-79 $μ$K/beam for total intensity and 13-23 $μ$K/beam for polarisation. For each region, we combine QUIJOTE maps with ancillary data at frequencies ranging from 0.4 to 3000 GHz, reconstruct the spectral energy distribution and model it with a combination of known foregrounds. We detect anomalous microwave emission (AME) in total intensity towards W49 at 4.7$σ$ and W51 at 4.0$σ$ with peak frequencies $ν_{AME}$ = (20.0 $\pm$ 1.4) GHz and $ν_{AME}$ = (17.7 $\pm$ 3.6) GHz respectively; this is the first detection of AME towards W51. The contamination from ultra-compact HII regions to the residual AME flux density is estimated at 10% in W49 and 5% in W51, and does not rule out the AME detection. The polarised SEDs reveal a synchrotron contribution with spectral indices $α_s$ = -0.67 $\pm$ 0.10 in W49 and $α_s$ = -0.51 $\pm$ 0.07 in W51, ascribed to the diffuse Galactic emission and to the local supernova remnant respectively. Towards IC443 in total intensity we measure a broken power-law synchrotron spectrum with cut-off frequency $ν_{0,s}$ = (114 $\pm$ 73) GHz, in agreement with previous studies; our analysis, however, rules out any AME contribution which had been previously claimed towards IC443. No evidence of polarised AME emission is detected in this study.

0 citations0 reviewsNo code linkedOpen access
preprint2023arXiv

QUIJOTE scientific results -- VI. The Haze as seen by QUIJOTE

The Haze is an excess of microwave intensity emission surrounding the Galactic centre. It is spatially correlated with the $γ$-ray Fermi bubbles, and with the S-PASS radio polarization plumes, suggesting a possible common provenance. The models proposed to explain the origin of the Haze, including energetic events at the Galactic centre and dark matter decay in the Galactic halo, do not yet provide a clear physical interpretation. In this paper we present a re-analysis of the Haze including new observations from the Multi-Frequency Instrument (MFI) of the Q-U-I JOint TEnerife (QUIJOTE) experiment, at 11 and 13 GHz. We analyze the Haze in intensity and polarization, characterizing its spectrum. We detect an excess of diffuse intensity signal ascribed to the Haze. The spectrum at frequencies 11$\,\leqν\leq\,$70 GHz is a power-law with spectral index $β^{\rm H}=-2.79\pm0.08$, which is flatter than the Galactic synchrotron in the same region ($β^{\rm S}=-2.98\pm0.04$), but steeper than that obtained from previous works ($β^{\rm H}\sim-2.5$ at 23$\,\leq\,ν\leq\,$70 GHz). We also observe an excess of polarized signal in the QUIJOTE-MFI maps in the Haze area. This is a first hint detection of polarized Haze, or a consequence of curvature of the synchrotron spectrum in that area. Finally, we show that the spectrum of polarized structures associated with Galactic centre activity is steep at low frequencies ($β\sim -3.2$ at 2.3 $\leqν\leq$ 23 GHz), and becomes flatter above 11 GHz.

0 citations0 reviewsNo code linkedOpen access
preprint2023arXiv

QUIJOTE Scientific Results -- VII. Galactic AME sources in the QUIJOTE-MFI Northern Hemisphere Wide-Survey

The QUIJOTE-MFI Northern Hemisphere Wide-Survey has provided maps of the sky above declinations $-30^\circ$ at 11, 13, 17 and 19$\,$GHz. These data are combined with ancillary data to produce Spectral Energy Distributions in intensity in the frequency range 0.4--3\,000$\,$GHz on a sample of 52 candidate compact sources harbouring anomalous microwave emission (AME). We apply a component separation analysis at 1$^\circ$ scale on the full sample from which we identify 44 sources with high AME significance. We explore correlations between different fitted parameters on this last sample. QUIJOTE-MFI data contribute to notably improve the characterisation of the AME spectrum, and its separation from the other components. In particular, ignoring the 10--20\,GHz data produces on average an underestimation of the AME amplitude, and an overestimation of the free-free component. We find an average AME peak frequency of 23.6 $\pm$ 3.6$\,$GHz, about 4$\,$GHz lower than the value reported in previous studies. The strongest correlation is found between the peak flux density of the thermal dust and of the AME component. A mild correlation is found between the AME emissivity ($A_{\rm AME}/τ_{250}$) and the interstellar radiation field. On the other hand no correlation is found between the AME emissivity and the free-free radiation Emission Measure. Our statistical results suggest that the interstellar radiation field could still be the main driver of the intensity of the AME as regards spinning dust excitation mechanisms. On the other hand, it is not clear whether spinning dust would be most likely associated with cold phases of the interstellar medium rather than with hot phases dominated by free-free radiation.

0 citations0 reviewsNo code linkedOpen access
preprint2023arXiv

QUIJOTE scientific results -- VIII. Diffuse polarized foregrounds from component separation with QUIJOTE-MFI

We derive linearly polarized astrophysical component maps in the Northern Sky from the QUIJOTE-MFI data at 11 and 13 GHz in combination with the WMAP K and Ka bands (23 and 33 GHz) and all Planck polarized channels (30-353 GHz), using the parametric component separation method B-SeCRET. The addition of QUIJOTE-MFI data significantly improves the parameter estimation of the low-frequency foregrounds, especially the estimation of the synchrotron spectral index, $β_s$. We present the first detailed $β_s$ map of the Northern Celestial Hemisphere at a smoothing scale of $2^{\circ}$. We find statistically significant spatial variability across the sky. We obtain an average value of $-3.08$ and a dispersion of $0.13$, considering only pixels with reliable goodness-of-fit. The power law model of the synchrotron emission provides a good fit to the data outside the Galactic plane but fails to track the complexity within this region. Moreover, when we assume a synchrotron model with uniform curvature, $c_s$, we find a value of $c_s = -0.0797 \pm 0.0012$. However, there is insufficient statistical significance to determine which model is favoured, either the power law or the power law with uniform curvature. Furthermore, we estimate the thermal dust spectral parameters in polarization. Our CMB, synchrotron, and thermal dust maps are highly correlated with the corresponding products of the PR4 Planck release, although some large-scale differences are observed in the synchrotron emission. Finally, we find that the $β_s$ estimation in the high signal-to-noise synchrotron emission areas is prior-independent while, outside these regions, the prior governs the $β_s$ estimation.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Improved cosmological fits with quantized primordial power spectra

We observationally examine cosmological models based on primordial power spectra with quantized wavevectors. Introducing a linearly quantized power spectrum with $k_0=3.225\times10^{-4}\mathrm{Mpc}^{-1}$ and spacing $Δk = 2.257 \times 10^{-4} \mathrm{Mpc}^{-1}$ provides a better fit to the Planck 2018 observations than the concordance baseline, with $Δχ^2 = -8.55$. Extending the results of Lasenby et al [1], we show that the requirement for perturbations to remain finite beyond the future conformal boundary in a universe containing dark matter and a cosmological constant results in a linearly quantized primordial power spectrum. It is found that the infrared cutoffs for this future conformal boundary quantized cosmology do not provide cosmic microwave background power spectra compatible with observations, but future theories may predict more observationally consistent quantized spectra.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Perturbations and the Future Conformal Boundary

The concordance model of cosmology predicts a universe which finishes in a finite amount of conformal time at a future conformal boundary. We show that for particular cases we study, the background variables and perturbations may be analytically continued beyond this boundary and that the "end of the universe" is not necessarily the end of their physical development. Remarkably, these theoretical considerations of the end of the universe might have observable consequences today: perturbation modes consistent with these boundary conditions have a quantised power spectrum which may be relevant to features seen in the large scale cosmic microwave background. Mathematically these cosmological models may either be interpreted as a palindromic universe mirrored in time, a reflecting boundary condition, or a double cover, but are identical with respect to their observational predictions and stand in contrast to the predictions of conformal cyclic cosmologies.

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

Blueshifted absorption lines from X-ray reflection in IRAS 13224-3809

We explore a disc origin for the highly-blueshifted, variable absorption lines seen in the X-ray spectrum of the Narrow Line Seyfert 1 galaxy IRAS13224-3809. The blueshift corresponds to a velocity of about 0.25c. Such features in other Active Galactic Nuclei are often interpreted as UltraFast Outflows (UFOs). The velocity is of course present in the orbital motions of the inner disk. The absorption lines in IRAS13224-3809 are best seen when the flux is low and the reflection component of the disk is strong relative to the power-law continuum. The spectra are consistent with a model in which the reflection component passes through a thin, highly-ionized absorbing layer at the surface of the inner disc, the blue-shifted side of which dominates the flux due to relativistic aberration (the disc inclination is about 70 deg). No fast outflow need occur beyond the disc.

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

Optical validation and characterisation of Planck PSZ1 sources at the Canary Islands observatories. II. Second year of ITP13 observations

We report new galaxy clusters previously unknown included in the PSZ1 catalogue. The results here presented have been achieved during the second year of a 2-year observational programme, the ITP13, developed at the Roque de los Muchachos Observatory (La Palma, Spain). We characterise 75 SZ sources with low SZ significance, SZ S/N$<5.32$ by performing deep optical imaging and spectroscopy in order to associate actual galaxy clusters to the SZ Planck source. We adopt robust criteria, based on the 2D-spatial distribution, richness and velocity dispersions to confirm actual optical counterparts up to $z<0.85$. At the end of the ITP13 observational programme, we study 256 SZ sources with $Dec \geq -15^{\circ}$ (212 of them completely unknown), finding optical counterparts for 152 SZ sources. The ITP13 validation programme has allowed us to update the PSZ1 purity, which is now more refined, increasing from 72\% to 83\% in the low SZ S/N regime. Our results are consistent with the predicted purity curve for the full PSZ1 catalogue and with the expected fraction of false detections caused by the non-Gaussian noise of foreground signals. Indeed, we find a strong correlation between the number of unconfirmed sources and the thermal emission of diffuse galactic dust at 857 GHz, thus increasing the fraction of false Planck SZ detections at low galactic latitudes.

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

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
preprint2013arXiv

Bayesian analysis of anisotropic cosmologies: Bianchi VII_h and WMAP

We perform a definitive analysis of Bianchi VII_h cosmologies with WMAP observations of the cosmic microwave background (CMB) temperature anisotropies. Bayesian analysis techniques are developed to study anisotropic cosmologies using full-sky and partial-sky, masked CMB temperature data. We apply these techniques to analyse the full-sky internal linear combination (ILC) map and a partial-sky, masked W-band map of WMAP 9-year observations. In addition to the physically motivated Bianchi VII_h model, we examine phenomenological models considered in previous studies, in which the Bianchi VII_h parameters are decoupled from the standard cosmological parameters. In the two phenomenological models considered, Bayes factors of 1.7 and 1.1 units of log-evidence favouring a Bianchi component are found in full-sky ILC data. The corresponding best-fit Bianchi maps recovered are similar for both phenomenological models and are very close to those found in previous studies using earlier WMAP data releases. However, no evidence for a phenomenological Bianchi component is found in the partial-sky W-band data. In the physical Bianchi VII_h model we find no evidence for a Bianchi component: WMAP data thus do not favour Bianchi VII_h cosmologies over the standard Lambda Cold Dark Matter (LCDM) cosmology. It is not possible to discount Bianchi VII_h cosmologies in favour of LCDM completely, but we are able to constrain the vorticity of physical Bianchi VII_h cosmologies at $(ω/H)_0 < 8.6 \times 10^{-10}$ with 95% confidence.

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 Comparison of Algorithms for the Construction of SZ Cluster Catalogues

We evaluate the construction methodology of an all-sky catalogue of galaxy clusters detected through the Sunyaev-Zel'dovich (SZ) effect. We perform an extensive comparison of twelve algorithms applied to the same detailed simulations of the millimeter and submillimeter sky based on a Planck-like case. We present the results of this "SZ Challenge" in terms of catalogue completeness, purity, astrometric and photometric reconstruction. Our results provide a comparison of a representative sample of SZ detection algorithms and highlight important issues in their application. In our study case, we show that the exact expected number of clusters remains uncertain (about a thousand cluster candidates at |b|> 20 deg with 90% purity) and that it depends on the SZ model and on the detailed sky simulations, and on algorithmic implementation of the detection methods. We also estimate the astrometric precision of the cluster candidates which is found of the order of ~2 arcmins on average, and the photometric uncertainty of order ~30%, depending on flux.

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
preprint2010arXiv

Power Spectrum Analyses of Nuclear Decay Rates

We provide the results from a spectral analysis of nuclear decay data displaying annually varying periodic fluctuations. The analyzed data were obtained from three distinct data sets: 32Si and 36Cl decays reported by an experiment performed at the Brookhaven National Laboratory (BNL), 56Mn decay reported by the Children's Nutrition Research Center (CNRC), but also performed at BNL, and 226Ra decay reported by an experiment performed at the Physikalisch-Technische Bundesanstalt (PTB) in Germany. All three data sets exhibit the same primary frequency mode consisting of an annual period. Additional spectral comparisons of the data to local ambient temperature, atmospheric pressure, relative humidity, Earth-Sun distance, and their reciprocals were performed. No common phases were found between the factors investigated and those exhibited by the nuclear decay data. This suggests that either a combination of factors was responsible, or that, if it was a single factor, its effects on the decay rate experiments are not a direct synchronous modulation. We conclude that the annual periodicity in these data sets is a real effect, but that further study involving additional carefully controlled experiments will be needed to establish its origin.

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 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

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
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
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
preprint1999arXiv

Microwave background anisotropies due to non-linear structures in open and Lambda universes

A new method arising from a gauge-theoretic approach to general relativity is applied to the formation of clusters in an expanding universe. The three cosmological models (Omega_0=1, Omega_Lambda=0), (Omega_0=0.3, Omega_Lambda=0.7), (Omega_0=0.3, Omega_Lambda=0) are considered, which extends our previous application (Lasenby et al. 1999 -"http://xxx.soton.ac.uk/abs/astro-ph/9810123"-, Dabrowski et al. 1999 -"http://xxx.soton.ac.uk/abs/astro-ph/9810149"-). A simple initial velocity and density perturbation of finite extent is imposed at the epoch z=1000 and we investigate the subsequent evolution of the density and velocity fields for clusters observed at redshifts z=1, z=2 and z=3. Photon geodesics and redshifts are also calculated so that the Cosmic Microwave Background (CMB) anisotropies due to collapsing clusters can be estimated. We find that the central CMB temperature decrement is slightly stronger and extends to larger angular scales in the non-zero Omega_Lambda case. This effect is strongly enhanced in the open case. Gravitational lensing effects are also considered and we apply our model to the reported microwave decrement observed towards the quasar pair PC 1643+4631 A & B.

0 citations0 reviewsNo code linkedOpen access
preprint1997arXiv

Studies of CMB structure at Dec=40. II: Analysis and cosmological interpretation

We present a detailed analysis of the cosmic microwave background structure in the Tenerife Dec=+40 degrees data. The effect of local atmospheric contributions on the derived fluctuation amplitude is considered, resulting in an improved separation of the intrinsic CMB signal from noise. Our analysis demonstrates the existence of common structure in independent data scans at 15 and 33 GHz. For the case of fluctuations described by a Gaussian auto-correlation function, a likelihood analysis of our combined results at 15 and 33 GHz implies an intrinsic rms fluctuation level of 48^{+21}_{-15} uK on a coherence scale of 4 degrees; the equivalent analysis for a Harrison-Zel'dovitch model gives a power spectrum normalisation of Q_{rms-ps} = 22^{+10}_{-6} uK. The fluctuation amplitude is seen to be consistent at the 68% confidence level with that reported for the COBE two-year data for primordial fluctuations described by a power law model with a spectral index in the range 1.0 \le n \le 1.6. This limit favours the large scale CMB anisotropy being dominated by scalar fluctuations rather than tensor modes from a gravitational wave background. The large scale Tenerife and COBE results are considered in conjunction with observational results from medium scale experiments in order to place improved limits on the fluctuation spectral index; we find n=1.10 +/- 0.10 assuming standard CDM with H_{0}=50 kms^{-1}Mpc^{-1}.

0 citations0 reviewsNo code linkedOpen access