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

Michael P. Hobson

Michael 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-ph.GAgr-qcastro-ph.IMastro-ph.SRastro-phastro-ph.HEhep-thMachine Learningphysics.data-anNeural and Evolutionary Computing

Catalog footprint

What is connected

46works
11topics
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

46 published item(s)

preprint2020arXiv

Comparison of physical and observational galaxy cluster modelling

We present a comparison between three cluster models applied to data obtained by the Arcminute Microkelvin Imager radio interferometer system. The physical model (PM) parameterises a cluster in terms of its physical quantities to model the dark matter and baryonic components of the cluster using NFW and GNFW profiles respectively. The observational models (OM I and OM II) model only the gas content of the cluster. The two OMs vary only in the priors they use in Bayesian inference: OM I has a joint prior on angular radius $θ$ and integrated Comptonisation $Y$, derived from simulations, while OM II uses separable priors on $θ$ and $Y$ which are based on calculations of the physical model. For the comparison we consider a sample of $54$ clusters which are a subsample of the second Planck catalogue of Sunyaev-Zel'dovich sources. We first compare the $Y$ estimates of the three models, and find that the PM generally yields lower estimates relative to the OMs. We then compute the Earth Mover's Distance between the $θ$ - $Y$ posterior distributions obtained from each model for each cluster, and find that the two models which are most discrepant are PM and OM I. Finally, we compare the Bayesian evidence values obtained from each model for each cluster. OM I generally provides the best fit to the data but not at a statistically significant level, according to the Jeffreys scale. The highest evidence ratio obtained is actually in favour of the PM over OM I.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Physical modelling of galaxy clusters detected by Planck

We present a comparison of mass estimates for $54$ galaxy cluster candidates from the second Planck catalogue (PSZ2) of Sunyaev-Zel'dovich sources. We compare the mass values obtained with data taken from the Arcminute Microkelvin Imager (AMI) radio interferometer system and from the Planck satellite. The former of these uses a Bayesian analysis pipeline that parameterises a cluster in terms of its physical quantities, and models the dark matter & baryonic components of a cluster using NFW and GNFW profiles respectively. Our mass estimates derived from Planck data are obtained from the results of the Bayesian detection algorithm PowellSnakes (PwS), are based on the methodology detailed in the PSZ2 paper, and produce two sets of mass estimates; one estimate is calculated directly from the angular radius $θ$ - integrated Comptonisation parameter $Y$ posterior distributions, and the other uses a `slicing function' to provide information on $θ$ based on X-ray measurements and previous Planck mission samples. We find that for $37$ of the clusters, the AMI mass estimates are lower than both values obtained from Planck data. However the AMI and slicing function estimates are within one combined standard deviation of each other for $31$ clusters. We also generate cluster simulations based on the slicing-function mass estimates, and analyse them in the same way as we did the real AMI data. We find that inclusion in the simulations of radio-source confusion & CMB noise and measurable radio-sources causes AMI mass estimates to be systematically low.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Power-counting renormalizable, ghost-and-tachyon-free Poincaré gauge theories

We present 48 further examples, in addition to the 10 identified in [1], of ghost-and-tachyon-free critical cases of parity-conserving Poincaré gauge theories of gravity (PGT$^+$) that are also power-counting renormalizable (PCR). This is achieved by extending the range of critical cases considered. Of the new PCR theories, seven have 2 massless degrees of freedom (d.o.f.) in propagating modes and a massive $0^-$ or $2^-$ mode, eight have only 2 massless d.o.f., and 33 have only massive mode(s). We also clarify the treatment of nonpropagating modes in determining whether a theory is PCR.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Sunyaev--Zel'dovich profile fitting with joint AMI-Planck analysis

We develop a Bayesian method of analysing Sunyaev-Zel'dovich measurements of galaxy clusters obtained from the Arcminute Microkelvin Imager (AMI) radio interferometer system and from the Planck satellite, using a joint likelihood function for the data from both instruments. Our method is applicable to any combination of Planck data with interferometric data from one or more arrays. We apply the analysis to simulated clusters and find that when the cluster pressure profile is known a-priori, the joint dataset provides precise and accurate constraints on the cluster parameters, removing the need for external information to reduce the parameter degeneracy. When the pressure profile deviates from that assumed for the fit, the constraints become biased. Allowing the pressure profile shape parameters to vary in the analysis allows an unbiased recovery of the integrated cluster signal and produces constraints on some shape parameters, depending on the angular size of the cluster. When applied to real data from Planck-detected cluster PSZ2 G063.80+11.42, our method resolves the discrepancy between the AMI and Planck $Y$-estimates and usefully constrains the gas pressure profile shape parameters at intermediate and large radii.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

AMI observations of ten CLASH galaxy clusters: SZ and X-ray data used together to determine cluster dynamical states

Using Arcminute Microkelvin Imager (AMI) SZ observations towards ten CLASH clusters we investigate the influence of cluster mergers on observational galaxy cluster studies. Although selected to be largely relaxed, there is disagreement in the literature on the dynamical states of CLASH sample members. We analyse our AMI data in a fully Bayesian way to produce estimated cluster parameters and consider the intrinsic correlations in our NFW/GNFW-based model. Varying pressure profile shape parameters, illustrating an influence of mergers on scaling relations, induces small deviations from the canonical self-similar predictions -- in agreement with simulations of Poole et al. 2007 who found that merger activity causes only small scatter perpendicular to the relations. We demonstrate this effect observationally using the different dependencies of SZ and X-ray signals to $n_{\rm e}$ that cause different sensitivities to the shocking and/or fractionation produced by mergers. Plotting $Y_{\rm X}$--$M_{\rm gas}$ relations (where $Y_{\rm X}=M_{\rm gas}T$) derived from AMI SZ and from $Chandra$ X-ray gives ratios of AMI and $Chandra$ $Y_{\rm X}$ and $M_{\rm gas}$ estimates that indicate movement of clusters \textit{along} the scaling relation, as predicted by Poole et al. 2007. Clusters that have moved most along the relation have the most discrepant $T_{\rm SZ}$ and $T_{\rm X}$ estimates: all the other clusters (apart from one) have SZ and X-ray estimates of $M_{\rm gas}$, $T$ and $Y_{\rm X}$ that agree within $r_{500}$. We use SZ vs X-ray discrepancies in conjunction with $Chandra$ maps and $T_{\rm X}$ profiles, making comparisons with simulated cluster merger maps in Poole et al. 2006, to identify disturbed members of our sample and estimate merger stages.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Friedmann-Robertson-Walker models do not require zero active mass

The $R_h = ct$ cosmological model has received considerable attention in recent years owing to claims that it is favoured over the standard $Λ$ cold dark matter ($Λ$CDM) model by most observational data. A key feature of the $R_h = ct$ model is that the zero active mass condition $ρ+ 3p = 0$ holds at all epochs. Most recently, Melia has claimed that this condition is a requirement of the symmetries of the Friedmann-Robertson-Walker (FRW) spacetime. We demonstrate that this claim is false and results from a flaw in the logic of Melia's argument.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Spherically-symmetric solutions in general relativity

We present a tetrad-based method for solving the Einstein field equations for spherically-symmetric systems and compare it with the widely-used Lemaître-Tolman-Bondi (LTB) model. In particular, we focus on the issues of gauge ambiguity and the use of comoving versus 'physical' coordinate systems. We also clarify the correspondences between the two approaches, and illustrate their differences by applying them to the classic examples of the Schwarzschild and Friedmann-Robertson-Walker spacetimes. We demonstrate that the tetrad-based method does not suffer from the gauge freedoms inherent to the LTB model, naturally accommodates non-zero pressure and has a more transparent physical interpretation. We further apply our tetrad-based method to a generalised form of 'Swiss cheese' model, which consists of an interior spherical region surrounded by a spherical shell of vacuum that is embedded in an exterior background universe. In general, we allow the fluid in the interior and exterior regions to support pressure, and do not demand that the interior region be compensated. We pay particular attention to the form of the solution in the intervening vacuum region and verify the validity of Birkhoff's theorem at both the metric and tetrad level. We then reconsider critically the original theoretical arguments underlying the so-called $R_h = ct$ cosmological model, which has recently received considerable attention. These considerations in turn illustrate the interesting behaviour of a number of 'horizons' in general cosmological models.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Bayesian estimation of non-Gaussianity in pulsar timing analysis

We introduce a method for performing a robust Bayesian analysis of non-Gaussianity present in pulsar timing data, simultaneously with the pulsar timing model, and additional stochastic parameters such as those describing red spin noise and dispersion measure variations. The parameters used to define the presence of non-Gaussianity are zero for Gaussian processes, giving a simple method of defining the strength of non-Gaussian behaviour. We use simulations to show that assuming Gaussian statistics when the noise in the data is drawn from a non-Gaussian distribution can significantly increase the uncertainties associated with the pulsar timing model parameters. We then apply the method to the publicly available 15 year Parkes Pulsar Timing Array data release 1 dataset for the binary pulsar J0437$-$4715. In this analysis we present a significant detection of non-Gaussianity in the uncorrelated non-thermal noise, but we find that it does not yet impact the timing model or stochastic parameter estimates significantly compared to analysis performed assuming Gaussian statistics. The methods presented are, however, shown to be of immediate practical use for current European Pulsar Timing Array (EPTA) and International Pulsar Timing Array (IPTA) datasets.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

SKYNET: an efficient and robust neural network training tool for machine learning in astronomy

We present the first public release of our generic neural network training algorithm, called SkyNet. This efficient and robust machine learning tool is able to train large and deep feed-forward neural networks, including autoencoders, for use in a wide range of supervised and unsupervised learning applications, such as regression, classification, density estimation, clustering and dimensionality reduction. SkyNet uses a `pre-training' method to obtain a set of network parameters that has empirically been shown to be close to a good solution, followed by further optimisation using a regularised variant of Newton's method, where the level of regularisation is determined and adjusted automatically; the latter uses second-order derivative information to improve convergence, but without the need to evaluate or store the full Hessian matrix, by using a fast approximate method to calculate Hessian-vector products. This combination of methods allows for the training of complicated networks that are difficult to optimise using standard backpropagation techniques. SkyNet employs convergence criteria that naturally prevent overfitting, and also includes a fast algorithm for estimating the accuracy of network outputs. The utility and flexibility of SkyNet are demonstrated by application to a number of toy problems, and to astronomical problems focusing on the recovery of structure from blurred and noisy images, the identification of gamma-ray bursters, and the compression and denoising of galaxy images. The SkyNet software, which is implemented in standard ANSI C and fully parallelised using MPI, is available at http://www.mrao.cam.ac.uk/software/skynet/.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

A joint analysis of AMI and CARMA observations of the recently discovered SZ galaxy cluster system AMI-CL J0300+2613

We present CARMA observations of a massive galaxy cluster discovered in the AMI blind SZ survey. Without knowledge of the cluster redshift a Bayesian analysis of the AMI, CARMA and joint AMI & CARMA uv-data is used to quantify the detection significance and parameterise both the physical and observational properties of the cluster whilst accounting for the statistics of primary CMB anisotropies, receiver noise and radio sources. The joint analysis of the AMI & CARMA uv-data was performed with two parametric physical cluster models: the β-model; and the model described in Olamaie et al. 2012 with the pressure profile fixed according to Arnaud et al. 2010. The cluster mass derived from these different models is comparable but our Bayesian evidences indicate a preference for the β-profile which we, therefore, use throughout our analysis. From the CARMA data alone we obtain a Bayesian probability of detection ratio of 12.8:1 when assuming that a cluster exists within our search area; alternatively assuming that Jenkins et al. 2001 accurately predicts the number of clusters as a function of mass and redshift, the Bayesian probability of detection is 0.29:1. From the analysis of the AMI or AMI & CARMA data the probability of detection ratio exceeds 4.5x10^3:1. Performing a joint analysis of the AMI & CARMA data with a physical cluster model we derive the total mass internal to r200 as MT,200 = 4.1x10^14Msun. Using a phenomenological β-model to quantify the temperature decrement as a function of angular distance we find a central SZ temperature decrement of 170μK in the AMI & CARMA data. The SZ decrement in the CARMA data is weaker than expected and we speculate that this is a consequence of the cluster morphology. In a forthcoming study we will assess the impact of cluster morphology on the SZ decrements that are observed with interferometers such as AMI and CARMA.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

AMI SZ observations and Bayesian analysis of a sample of six redshift-one clusters of galaxies

We present 16-GHz Sunyaev-Zel'dovich observations using the Arcminute Microkelvin Imager (AMI) and subsequent Bayesian analysis of six galaxy clusters at redshift ($z \approx 1$) chosen from an X-ray and Infrared selected sample from Culverhouse et al. (2010). In the subsequent analysis we use two cluster models, an isothermal β-model and a Dark Matter GNFW (DM-GNFW) model in order to derive a formal detection probability and the cluster parameters. We detect two clusters (CLJ1415+3612 & XMJ0830+5241) and measure their total masses out to a radius of 200 $\times$ the critical density at the respective cluster's redshift. For CLJ1415+3612 and XMJ0830+5241, we find M_{\mathrm{T},200} for each model, which agree with each other for each cluster. We also present maps before and after source subtraction of the entire sample and provide 1D and 2D posterior marginalised probability distributions for each fitted cluster profile parameter of the detected clusters. Using simulations which take into account the measured source environment from the AMI Large Array (LA), source confusion noise, CMB primordials, instrument noise, we estimate from low-radius X-ray data from Culverhouse et al. (2010), the detectability of each cluster in the sample and compare it with the result from the Small Array (SA) data. Furthermore, we discuss the validity of the assumptions of isothermality and constant gas mass fraction. We comment on the bias that these small-radius estimates introduce to large-radius SZ predictions. In addition, we follow-up the two detections with deep, single-pointed LA observations. We find a 3 sigma tentative decrement toward CLJ1415+3612 at high-resolution and a 5 sigma high-resolution decrement towards XMJ0830+5241.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Hyper-efficient model-independent Bayesian method for the analysis of pulsar timing data

A new model independent method is presented for the analysis of pulsar timing data and the estimation of the spectral properties of an isotropic gravitational wave background (GWB). We show that by rephrasing the likelihood we are able to eliminate the most costly aspects of computation normally associated with this type of data analysis. When applied to the International Pulsar Timing Array Mock Data Challenge data sets this results in speedups of approximately 2 to 3 orders of magnitude compared to established methods. We present three applications of the new likelihood. In the low signal to noise regime we sample directly from the power spectrum coefficients of the GWB signal realization. In the high signal to noise regime, where the data can support a large number of coefficients, we sample from the joint probability density of the power spectrum coefficients for the individual pulsars and the GWB signal realization. Critically in both these cases we need make no assumptions about the form of the power spectrum of the GWB, or the individual pulsars. Finally we present a method for characterizing the spatial correlation between pulsars on the sky, making no assumptions about the form of that correlation, and therefore providing the only truly general Bayesian method of confirming a GWB detection from pulsar timing data.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

SZ observations with AMI of the hottest galaxy clusters detected in the XMM-Newton Cluster Survey

We have obtained deep SZ observations towards 15 of the apparently hottest XMM Cluster Survey (XCS) clusters that can be observed with the Arcminute Microkelvin Imager (AMI). We use a Bayesian analysis to quantify the significance of our SZ detections. We detect the SZ effect at high significance towards three of the clusters and at lower significance for a further two clusters. Towards the remaining ten clusters, no clear SZ signal was measured. We derive cluster parameters using the XCS mass estimates as a prior in our Bayesian analysis. For all AMI-detected clusters, we calculate large-scale mass and temperature estimates while for all undetected clusters we determine upper limits on these parameters. We find that the large- scale mean temperatures derived from our AMI SZ measurements (and the upper limits from null detections) are substantially lower than the XCS-based core-temperature estimates. For clusters detected in the SZ, the mean temperature is, on average, a factor of 1.4 lower than temperatures from the XCS. For clusters undetected in SZ, the average 68% upper limit on the mean temperature is a factor of 1.9 below the XCS temperature.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

TempoNest: A Bayesian approach to pulsar timing analysis

A new Bayesian software package for the analysis of pulsar timing data is presented in the form of TempoNest which allows for the robust determination of the non-linear pulsar timing solution simultaneously with a range of additional stochastic parameters. This includes both red spin noise and dispersion measure variations using either power law descriptions of the noise, or through a model-independent method that parameterises the power at individual frequencies in the signal. We use TempoNest to show that at noise levels representative of current datasets in the European Pulsar Timing Array (EPTA) and International Pulsar Timing Array (IPTA) the linear timing model can underestimate the uncertainties of the timing solution by up to an order of magnitude. We also show how to perform Bayesian model selection between different sets of timing model and stochastic parameters, for example, by demonstrating that in the pulsar B1937+21 both the dispersion measure variations and spin noise in the data are optimally modelled by simple power laws. Finally we show that not including the stochastic parameters simultaneously with the timing model can lead to unpredictable variation in the estimated uncertainties, compromising the robustness of the scientific results extracted from such analysis.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

AMI Galactic Plane Survey at 16 GHz: I -- Observing, mapping and source extraction

The AMI Galactic Plane Survey (AMIGPS) is a large area survey of the outer Galactic plane to provide arcminute resolution images at milli-Jansky sensitivity in the centimetre-wave band. Here we present the first data release of the survey, consisting of 868 deg^2 of the Galactic plane, covering the area 76 deg \lessapprox l \lessapprox 170 deg between latitudes of |b| \lessapprox 5 deg, at a central frequency of 15.75 GHz (1.9 cm). We describe in detail the drift scan observations which have been used to construct the maps, including the techniques used for observing, mapping and source extraction, and summarise the properties of the finalized datasets. These observations constitute the most sensitive Galactic plane survey of large extent at centimetre-wave frequencies greater than 1.4 GHz.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

AMI radio continuum observations of young stellar objects with known outflows

We present 16 GHz (1.9 cm) deep radio continuum observations made with the Arcminute Microkelvin Imager (AMI) of a sample of low-mass young stars driving jets. We combine these new data with archival information from an extensive literature search to examine spectral energy distributions (SEDs) for each source and calculate both the radio and sub-mm spectral indices in two different scenarios: (1) fixing the dust temperature (Td) according to evolutionary class; (2) allowing Td to vary. We use the results of this analysis to place constraints on the physical mechanisms responsible for the radio emission. From AMI data alone, as well as from model fitting to the full SED in both scenarios, we find that 80 per cent of the objects in this sample have spectral indices consistent with free-free emission. We find an average spectral index in both Td scenarios consistent with free-free emission. We examine correlations of the radio luminosity with bolometric luminosity, envelope mass, and outflow force and find that these data are consistent with the strong correlation with envelope mass seen in lower luminosity samples. We examine the errors associated with determining the radio luminosity and find that the dominant source of error is the uncertainty on the opacity index, beta. We examine the SEDs for variability in these young objects, and find evidence for possible radio flare events in the histories of L1551 IRS 5 and Serpens SMM 1.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

BAMBI: blind accelerated multimodal Bayesian inference

In this paper we present an algorithm for rapid Bayesian analysis that combines the benefits of nested sampling and artificial neural networks. The blind accelerated multimodal Bayesian inference (BAMBI) algorithm implements the MultiNest package for nested sampling as well as the training of an artificial neural network (NN) to learn the likelihood function. In the case of computationally expensive likelihoods, this allows the substitution of a much more rapid approximation in order to increase significantly the speed of the analysis. We begin by demonstrating, with a few toy examples, the ability of a NN to learn complicated likelihood surfaces. BAMBI's ability to decrease running time for Bayesian inference is then demonstrated in the context of estimating cosmological parameters from Wilkinson Microwave Anisotropy Probe and other observations. We show that valuable speed increases are achieved in addition to obtaining NNs trained on the likelihood functions for the different model and data combinations. These NNs can then be used for an even faster follow-up analysis using the same likelihood and different priors. This is a fully general algorithm that can be applied, without any pre-processing, to other problems with computationally expensive likelihood functions.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Bayesian analysis of weak gravitational lensing and Sunyaev-Zel'dovich data for six galaxy clusters

We present an analysis of observations made with the Arcminute Microkelvin Imager (AMI) and the Canada-France-Hawaii Telescope (CFHT) of six galaxy clusters in a redshift range of 0.16--0.41. The cluster gas is modelled using the Sunyaev--Zel'dovich (SZ) data provided by AMI, while the total mass is modelled using the lensing data from the CFHT. In this paper, we: i) find very good agreement between SZ measurements (assuming large-scale virialisation and a gas-fraction prior) and lensing measurements of the total cluster masses out to r_200; ii) perform the first multiple-component weak-lensing analysis of A115; iii) confirm the unusual separation between the gas and mass components in A1914; iv) jointly analyse the SZ and lensing data for the relaxed cluster A611, confirming our use of a simulation-derived mass-temperature relation for parameterizing measurements of the SZ effect.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Detailed SZ study of 19 LoCuSS galaxy clusters: masses and temperatures out to the virial radius

We present 16-GHz AMI SZ observations of 19 clusters with L_X >7x10^37 W (h50=1) selected from the LoCuS survey (0.142<z<0.295) and of A1758b, in the FoV of A1758a. We detect 17 clusters with 5-23sigma peak surface brightnesses. Cluster parameters are obtained using a Bayesian cluster analysis. We fit isothermal beta-models to our data and assume the clusters are virialized (with all the kinetic energy in gas internal energy). Our gas temperature, T_AMI, is derived from AMI SZ data, not from X-ray spectroscopy. Cluster parameters internal to r500 are derived assuming HSE. We find: (i) Different gNFW parameterizations yield significantly different parameter degeneracies. (ii) For h70 = 1, we find the virial radius r200 to be typically 1.6+/-0.1 Mpc and the total mass M_T(r200) typically to be 2.0-2.5xM_T(r500).(iii) Where we have found M_T X-ray (X) and weak-lensing (WL) values in the literature, there is good agreement between WL and AMI estimates (with M_{T,AMI}/M_{T,WL} =1.2^{+0.2}_{-0.3} and =1.0+/-0.1 for r500 and r200, respectively). In comparison, most Suzaku/Chandra estimates are higher than for AMI (with M_{T,X}/M_{T,AMI}=1.7+/-0.2 within r500), particularly for the stronger mergers.(iv) Comparison of T_AMI to T_X sheds light on high X-ray masses: even at large r, T_X can substantially exceed T_AMI in mergers. The use of these higher T_X values will give higher X-ray masses. We stress that large-r T_SZ and T_X data are scarce and must be increased. (v) Despite the paucity of data, there is an indication of a relation between merger activity and SZ ellipticity. (vi) At small radius (but away from any cooling flow) the SZ signal (and T_AMI) is less sensitive to ICM disturbance than the X-ray signal (and T_X) and, even at high r, mergers affect n^2-weighted X-ray data more than n-weighted SZ, implying significant shocking or clumping or both occur even in the outer parts of mergers.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Parameterization Effects in the analysis of AMI Sunyaev-Zel'dovich Observations

Most Sunyaev--Zel'dovich (SZ) and X-ray analyses of galaxy clusters try to constrain the cluster total mass and/or gas mass using parameterised models and assumptions of spherical symmetry and hydrostatic equilibrium. By numerically exploring the probability distributions of the cluster parameters given the simulated interferometric SZ data in the context of Bayesian methods, and assuming a beta-model for the electron number density we investigate the capability of this model and analysis to return the simulated cluster input quantities via three rameterisations. In parameterisation I we assume that the T is an input parameter. We find that parameterisation I can hardly constrain the cluster parameters. We then investigate parameterisations II and III in which fg(r200) replaces temperature as a main variable. In parameterisation II we relate M_T(r200) and T assuming hydrostatic equilibrium. We find that parameterisation II can constrain the cluster physical parameters but the temperature estimate is biased low. In parameterisation III, the virial theorem replaces the hydrostatic equilibrium assumption. We find that parameterisation III results in unbiased estimates of the cluster properties. We generate a second simulated cluster using a generalised NFW (GNFW) pressure profile and analyse it with an entropy based model to take into account the temperature gradient in our analysis and improve the cluster gas density distribution. This model also constrains the cluster physical parameters and the results show a radial decline in the gas temperature as expected. The mean cluster total mass estimates are also within 1 sigma from the simulated cluster true values. However, we find that for at least interferometric SZ analysis in practice at the present time, there is no differences in the AMI visibilities between the two models. This may of course change as the instruments improve.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Sunyaev-Zel'dovich observations of LoCuSS clusters with the Arcminute Microkelvin Imager: moderate X-ray luminosity sample

We present 16-GHz observations using the Arcminute Microkelvin Imager (AMI) of 11 clusters with 7 x 10^{37}W < L_X < 11 x 10^{37}W (h_{50}=1.0) selected from the Local Cluster Substructure Survey (LoCuSS) and compare them to X-ray data. We use a fast, Bayesian cluster analysis to explore the high-dimensional parameter space of the cluster-plus-sources model and obtain robust cluster parameter estimates in the presence of radio point sources, receiver noise and primordial CMB anisotropy. Our analysis fits a spherical, isothermal beta-model to our data and assumes the cluster follows the theoretical mass-temperature relation. Large-scale cluster parameters internal to r_{500} are derived under the assumption of hydrostatic equilibrium. Posterior distributions for the large-scale parameters of 8 of our clusters are given; SZ effects towards Abell 1704 and Zw0857.9+2107 were not detected and our spherical beta-profile was found to be an inadequate fit to the decrement on our map for Abell 2409.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

The effect of a massive object on an expanding universe

A tetrad-based procedure is presented for solving Einstein's field equations for spherically-symmetric systems; this approach was first discussed by Lasenby et al. in the language of geometric algebra. The method is used to derive metrics describing a point mass in a spatially-flat, open and closed expanding universe respectively. In the spatially-flat case, a simple coordinate transformation relates the metric to the corresponding one derived by McVittie. Nonetheless, our use of non-comoving (`physical') coordinates greatly facilitates physical interpretation. For the open and closed universes, our metrics describe different spacetimes to the corresponding McVittie metrics and we believe the latter to be incorrect. In the closed case, our metric possesses an image mass at the antipodal point of the universe. We calculate the geodesic equations for the spatially-flat metric and interpret them. For radial motion in the Newtonian limit, the force acting on a test particle consists of the usual $1/r^2$ inwards component due to the central mass and a cosmological component proportional to $r$ that is directed outwards (inwards) when the expansion of the universe is accelerating (decelerating). For the standard $Λ$CDM concordance cosmology, the cosmological force reverses direction at about $z\approx 0.67$. We also derive an invariant fully general-relativistic expression, valid for arbitrary spherically-symmetric systems, for the force required to hold a test particle at rest relative to the central point mass.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

The effect of an expanding universe on massive objects

We present some astrophysical consequences of the metric for a point mass in an expanding universe derived in Nandra, Lasenby & Hobson, and of the associated invariant expression for the force required to keep a test particle at rest relative to the central mass. We focus on the effect of an expanding universe on massive objects on the scale of galaxies and clusters. Using Newtonian and general-relativistic approaches, we identify two important time-dependent physical radii for such objects when the cosmological expansion is accelerating. The first radius, $r_F$, is that at which the total radial force on a test particle is zero, which is also the radius of the largest possible circular orbit about the central mass $m$ and where the gas pressure and its gradient vanish. The second radius, $r_S$, which is \approx r_F/1.6$, is that of the largest possible stable circular orbit, which we interpret as the theoretical maximum size for an object of mass $m$. In contrast, for a decelerating cosmological expansion, no such finite radii exist. Assuming a cosmological expansion consistent with a $Λ$CDM concordance model, at the present epoch we find that these radii put a sensible constraint on the typical sizes of both galaxies and clusters at low redshift. For galaxies, we also find that these radii agree closely with zeroes in the radial velocity field in the neighbourhood of nearby galaxies, as inferred by Peirani & Pacheco from recent observations of stellar velocities. We then consider the future effect on massive objects of an accelerating cosmological expansion driven by phantom energy, for which the universe is predicted to end in a `Big Rip' at a finite time in the future at which the scale factor becomes singular. In particular, we present a novel calculation of the time prior to the Big Rip that an object of a given mass and size will become gravitationally unbound.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

10C Survey of Radio Sources at 15.7 GHz: I - Observing, mapping and source extraction

We have observed an area of approximatley 27 deg^2 to an rms noise level of less than 0.2 mJy at 15.7 GHz, using the Arcminute Microkelvin Imager Large Array. These observations constitute the most sensitive radio-source survey of any extent (greater than approximately 0.2 deg^2) above 1.4 GHz. This paper presents the techniques employed for observing, mapping and source extraction. We have used a systematic procedure for extracting information and producing source catalogues, from maps with varying noise and uv-coverage. We have performed simulations to test our mapping and source-extraction procedures, and developed methods for identifying extended, overlapping and spurious sources in noisy images. In an accompanying paper, AMI Consortium: Davies et al. 2010, the first results from the 10C survey, including the deep 15.7-GHz source count, are presented.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

10C Survey of Radio Sources at 15.7 GHz: II - First Results

The first results from the Tenth Cambridge (10C) Survey of Radio Sources, carried out using the AMI Large Array (LA) at an observing frequency of 15.7 GHz, are presented. The survey fields cover an area of approximately 27 sq. degrees to a flux-density completeness of 1 mJy. Results for some deeper areas, covering approximately 12 sq. degrees, wholly contained within the total areas and complete to 0.5 mJy, are also presented. The completeness for both areas is estimated to be at least 93 per cent. The source catalogue contains 1897 entries and is available at www.mrao.cam.ac.uk/surveys/10C. It has been combined with that of the 9C Survey to calculate the 15.7-GHz source counts. A broken power law is found to provide a good parameterisation of the differential count between 0.5 mJy and 1 Jy. The measured count has been compared to that predicted by de Zotti et al. (2005). The model displays good agreement with the data at the highest flux densities but under-predicts the integrated count between 0.5 mJy and 1 Jy by about 30 per cent. Entries from the source catalogue have been matched to those contained in the catalogues of NVSS and FIRST (both of which have observing frequencies of 1.4 GHz). This matching provides evidence for a shift in the typical 1.4-to-15.7-GHz spectral index of the 15.7-GHz-selected source population with decreasing flux density towards sub-mJy levels - the spectra tend to become less steep. Automated methods for detecting extended sources have been applied to the data; approximately 5 per cent of the sources are found to be extended relative to the LA synthesised beam of approximately 30 arcsec. Investigations using higher-resolution data showed that most of the genuinely extended sources at 16 GHz are classical doubles, although some nearby galaxies and twin-jet sources were also identified.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

AMI observations of unmatched Planck ERCSC LFI sources at 15.75 GHz

The Planck Early Release Compact Source Catalogue includes 26 sources with no obvious matches in other radio catalogues (of primarily extragalactic sources). Here we present observations made with the Arcminute Microkelvin Imager Small Array (AMI SA) at 15.75 GHz of the eight of the unmatched sources at declination > +10 degrees. Of the eight, four are detected and are associated with known objects. The other four are not detected with the AMI SA, and are thought to be spurious.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

AMI-LA radio continuum observations of Spitzer c2d small clouds and cores: Perseus region

We present deep radio continuum observations of the cores identified as deeply embedded young stellar objects in the Perseus molecular cloud by the Spitzer c2d programme at a wavelength of 1.8 cm with the Arcminute Microkelvin Imager Large Array (AMI-LA). We detect 72% of Class 0 objects from this sample and 31% of Class I objects. No starless cores are detected. We use the flux densities measured from these data to improve constraints on the correlations between radio luminosity and bolometric luminosity, infrared luminosity and, where measured, outflow force. We discuss the differing behaviour of these objects as a function of protostellar class and investigate the differences in radio emission as a function of core mass. Two of four possible very low luminosity objects (VeLLOs) are detected at 1.8 cm.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

AMI-LA radio continuum observations of Spitzer c2d small clouds and cores: Serpens region

We present deep radio continuum observations of the cores identified as deeply embedded young stellar objects in the Serpens molecular cloud by the Spitzer c2d programme at a wavelength of 1.8cm with the Arcminute Microkelvin Imager Large Array (AMI-LA). These observations have a resolution of ~30arcsec and an average sensitivity of 19microJy/beam. The targets are predominantly Class I sources, and we find the detection rate for Class I objects in this sample to be low (18%) compared to that of Class 0 objects (67%), consistent with previous works. For detected objects we examine correlations of radio luminosity with bolometric luminosity and envelope mass and find that these data support correlations found by previous samples, but do not show any indiction of the evolutionary divide hinted at by similar data from the Perseus molecular cloud when comparing radio luminosity with envelope mass. We conclude that envelope mass provides a better indicator for radio luminosity than bolometric luminosity, based on the distribution of deviations from the two correlations. Combining these new data with archival 3.6cm flux densities we also examine the spectral indices of these objects and find an average spectral index of 0.53+/-1.14, consistent with the canonical value for a partially optically thick spherical or collimated stellar wind. However, we caution that possible inter-epoch variability limits the usefulness of this value, and such variability is supported by our identification of a possible flare in the radio history of Serpens SMM 1.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Further Sunyaev-Zel'dovich observations of two Planck ERCSC clusters with the Arcminute Microkelvin Imager

We present follow-up observations of two galaxy clusters detected blindly via the Sunyaev-Zel'dovich (SZ) effect and released in the Planck Early Release Compact Source Catalogue. We use the Arcminute Microkelvin Imager, a dual-array 14-18 GHz radio interferometer. After radio source subtraction, we find a SZ decrement of integrated flux density -1.08+/-0.10 mJy toward PLCKESZ G121.11+57.01, and improve the position measurement of the cluster, finding the centre to be RA 12 59 36.4, Dec +60 04 46.8, to an accuracy of 20 arcseconds. The region of PLCKESZ G115.71+17.52 contains strong extended emission, so we are unable to confirm the presence of this cluster via the SZ effect.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Radio continuum observations of Class I protostellar disks in Taurus: constraining the greybody tail at centimetre wavelengths

We present deep 1.8 cm (16 GHz) radio continuum imaging of seven young stellar objects in the Taurus molecular cloud. These objects have previously been extensively studied in the sub-mm to NIR range and their SEDs modelled to provide reliable physical and geometrical parametres.We use this new data to constrain the properties of the long-wavelength tail of the greybody spectrum, which is expected to be dominated by emission from large dust grains in the protostellar disk. We find spectra consistent with the opacity indices expected for such a population, with an average opacity index of beta = 0.26+/-0.22 indicating grain growth within the disks. We use spectra fitted jointly to radio and sub-mm data to separate the contributions from thermal dust and radio emission at 1.8 cm and derive disk masses directly from the cm-wave dust contribution. We find that disk masses derived from these flux densities under assumptions consistent with the literature are systematically higher than those calculated from sub-mm data, and meet the criteria for giant planet formation in a number of cases.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Sunyaev-Zel'dovich observation of the Bullet-like cluster Abell 2146 with the Arcminute Microkelvin Imager

We present 13.9-18.2 GHz observations of the Sunyaev-Zel'dovich (SZ) effect towards Abell 2146 using the Arcminute Microkelvin Imager (AMI). The cluster is detected with a peak SNR ratio of 13 sigma in the radio source subtracted map. Comparison of the SZ and X-ray images suggests that they both have extended regions which lie approximately perpendicular to one another, with their emission peaks significantly displaced. These features indicate non-uniformities in the distributions of the gas temperature and pressure, indicative of a cluster merger. We use a Bayesian cluster analysis to explore the high-dimensional parameter space of the cluster-plus-sources model to obtain cluster parameter estimates in the presence of radio point sources, receiver noise and primordial CMB anisotropy; the probability of SZ + CMB primordial structure + radio sources + receiver noise to CMB + radio sources + receiver noise is 3 x 10^{6}:1. We compare the results from three different cluster models. Our preferred model exploits the observation that the gas fractions do not appear to vary greatly between clusters. Given the relative masses of the two merging systems in Abell 2146, the mean gas temperature can be deduced from the virial theorem (assuming all of the kinetic energy is in the form of internal gas energy) without being affected significantly by the merger event, provided the primary cluster was virialized before the merger. In this model we fit a simple spherical isothermal beta-model, despite the inadequacy of this model for a merging system like Abell 2146, and assume the cluster follows the mass-temperature relation of a virialized, singular, isothermal sphere. We note that this model avoids inferring large-scale cluster parameters internal to r_200 under the widely used assumption of hydrostatic equilibrium. We find that at r_200 M_T= 4.1 \pm 0.5 x 10^{14} h^{-1}M_sun and T=4.5 \pm 0.5 keV.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Sunyaev-Zel'dovich observations of LoCuSS clusters with the Arcminute Microkelvin Imager: high X-ray luminosity sample

We present observations from the Small Array of the Arcminute Microkelvin Imager (AMI) of eight high X-ray luminosity galaxy cluster systems selected from the Local Cluster Substructure Survey (LoCuSS) sample.We detect the Sunyaev-Zel'dovich (SZ) effect in seven of these clusters. With the assumptions that galaxy clusters are isothermal, have a density profile described by a spherical b -model and obey the theoretical M-T relation, we are able to derive cluster parameters at r200 from our SZ data. With the additional assumption of hydrostatic equilibrium we are able to derive parameters at r500. We present posterior probability distributions for cluster parameters such as mass, radius and temperature (TSZ, MT). Combining our sample with that of AMI Consortium: Rodr'\iguez-Gonzálvez et al. (2011) and using large-radius X-ray temperature estimates (TX) from Chandra and Suzaku observations, we find that there is reasonable correspondence between TX and TSZ,MT values at low TX, but that for clusters with TX above around 6keV the correspondence breaks down with TX exceeding TSZ, MT; we stress that this finding is based on just ten clusters.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

AMI Large Array radio continuum observations of Spitzer c2d small clouds and cores

We perform deep 1.8 cm radio continuum imaging towards thirteen protostellar regions selected from the Spitzer c2d small clouds and cores programme at high resolution (25") in order to detect and quantify the cm-wave emission from deeply embedded young protostars. Within these regions we detect fifteen compact radio sources which we identify as radio protostars including two probable new detections. The sample is in general of low bolometric luminosity and contains several of the newly detected VeLLO sources. We determine the 1.8 cm radio luminosity to bolometric luminosity correlation, L_rad -L_bol, for the sample and discuss the nature of the radio emission in terms of the available sources of ionized gas. We also investigate the L_rad-L_IR correlation and suggest that radio flux density may be used as a proxy for the internal luminosity of low luminosity protostars.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Cracking the Taub-NUT

We present further analysis of an anisotropic, non-singular early universe model that leads to the viable cosmology presented in Dechant et al (arXiv:0809.4335). Although this model (the DLH model) contains scalar field matter, it is reminiscent of the Taub-NUT vacuum solution in that it has biaxial Bianchi IX geometry and its evolution exhibits a dimensionality reduction at a quasi-regular singularity that one can identify with the big-bang. We show that the DLH and Taub-NUT metrics are related by a coordinate transformation, in which the DLH time coordinate plays the role of conformal time for Taub-NUT. Since both models continue through the big-bang, the coordinate transformation can become multivalued. In particular, in mapping from DLH to Taub-NUT, the Taub-NUT time can take only positive values. We present explicit maps between the DLH and Taub-NUT models, with and without a scalar field. In the vacuum DLH model, we find a periodic solution expressible in terms of elliptic integrals. Mapping the vacuum solution over to Taub-NUT coordinates, recovers the standard (non-periodic) Taub-NUT solution in the Taub region, where Taub-NUT time takes positive values, but does not exhibit the two NUT regions known in the standard Taub-NUT solution. Conversely, mapping the complete Taub-NUT solution to the DLH case reveals that the NUT regions correspond to imaginary time and space in DLH coordinates. We show that many of the well-known `pathologies' of the Taub-NUT solution arise because the traditional coordinates are connected by a multivalued transformation to the physically more meaningful DLH coordinates. In particular, the `open-to-closed-to-open' transition and the Taub and NUT regions of the (Lorentzian) Taub-NUT model are replaced by a closed pancaking universe with spacelike homogeneous sections at all times.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

High resolution AMI Large Array imaging of spinning dust sources: spatially correlated 8 micron emission and evidence of a stellar wind in L675

We present 25 arcsecond resolution radio images of five Lynds Dark Nebulae (L675, L944, L1103, L1111 & L1246) at 16 GHz made with the Arcminute Microkelvin Imager (AMI) Large Array. These objects were previously observed with the AMI Small Array to have an excess of emission at microwave frequencies relative to lower frequency radio data. In L675 we find a flat spectrum compact radio counterpart to the 850 micron emission seen with SCUBA and suggest that it is cm-wave emission from a previously unknown deeply embedded young protostar. In the case of L1246 the cm-wave emission is spatially correlated with 8 micron emission seen with Spitzer. Since the MIR emission is present only in Spitzer band 4 we suggest that it arises from a population of PAH molecules, which also give rise to the cm-wave emission through spinning dust emission.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Microwave observations of spinning dust emission in NGC6946

We report new cm-wave measurements at five frequencies between 15 and 18GHz of the continuum emission from the reportedly anomalous "region 4" of the nearby galaxy NGC6946. We find that the emission in this frequency range is significantly in excess of that measured at 8.5GHz, but has a spectrum from 15-18GHz consistent with optically thin free-free emission from a compact HII region. In combination with previously published data we fit four emission models containing different continuum components using the Bayesian spectrum analysis package radiospec. These fits show that, in combination with data at other frequencies, a model with a spinning dust component is slightly preferred to those that possess better-established emission mechanisms.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Sunyaev-Zel'dovich observations of galaxy clusters out to the virial radius with the Arcminute Microkelvin Imager

We present observations using the Small Array of the Arcminute Microkelvin Imager (AMI; 14-18 GHz) of four Abell and three MACS clusters spanning 0.171-0.686 in redshift. We detect Sunyaev-Zel'dovich (SZ) signals in five of these without any attempt at source subtraction, although strong source contamination is present. With radio-source measurements from high-resolution observations, and under the assumptions of spherical $β$-model, isothermality and hydrostatic equilibrium, a Bayesian analysis of the data in the visibility plane detects extended SZ decrements in all seven clusters over and above receiver noise, radio sources and primary CMB imprints. Bayesian evidence ratios range from 10^{11}:1 to 10^{43}:1 for six of the clusters and 3000:1 for one with substantially less data than the others. We present posterior probability distributions for, e.g., total mass and gas fraction averaged over radii internal to which the mean overdensity is 1000, 500 and 200, r_200 being the virial radius. Reaching r_200 involves some extrapolation for the nearer clusters but not for the more-distant ones. We find that our estimates of gas fraction are low (compared with most in the literature) and decrease with increasing radius. These results appear to be consistent with the notion that gas temperature in fact falls with distance (away from near the cluster centre) out to the virial radius.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

AMI observations of Lynds Dark Nebulae: further evidence for anomalous cm-wave emission

Observations at 14.2 to 17.9 GHz made with the AMI Small Array towards fourteen Lynds Dark Nebulae with a resolution of 2' are reported. These sources are selected from the SCUBA observations of Visser et al. (2001) as small angular diameter clouds well matched to the synthesized beam of the AMI Small Array. Comparison of the AMI observations with radio observations at lower frequencies with matched uv-plane coverage is made, in order to search for any anomalous excess emission which can be attributed to spinning dust. Possible emission from spinning dust is identified as a source within a 2' radius of the Scuba position of the Lynds dark nebula, exhibiting an excess with respect to lower frequency radio emission. We find five sources which show a possible spinning dust component in their spectra. These sources have rising spectral indices in the frequency range 14.2--17.9 GHz. Of these five one has already been reported, L1111, we report one new definite detection, L675, and three new probable detections (L944, L1103 and L1246). The relative certainty of these detections is assessed on the basis of three criteria: the extent of the emission, the coincidence of the emission with the Scuba position and the likelihood of alternative explanations for the excess. Extended microwave emission makes the likelihood of the anomalous emission arising as a consequence of a radio counterpart to a protostar or a proto-planetary disk unlikely. We use a 2' radius in order to be consistent with the IRAS identifications of dark nebulae (Parker 1988), and our third criterion is used in the case of L1103 where a high flux density at 850 microns relative to the FIR data suggests a more complicated emission spectrum.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

AMI observations of northern supernova remnants at 14-18 GHz

We present observations between 14.2 and 17.9 GHz of 12 reported supernova remnants (SNRs) made with the Arcminute Microkelvin Imager Small Array (AMI SA). In conjunction with data from the literature at lower radio frequencies, we determine spectra of these objects. For well-studied SNRs (Cas A, Tycho's SNR, 3C58 and the Crab Nebula), the results are in good agreement with spectra based on previous results. For the less well-studied remnants the AMI SA observations provide higher-frequency radio observations than previously available, and better constrain their radio spectra. The AMI SA results confirm a spectral turnover at ~11 GHz for the filled-centre remnant G74.9+1.2. We also see a possible steepening of the spectrum of the filled-centre remnant G54.1+0.3 within the AMI SA frequency band compared with lower frequencies. We confirm that G84.9+0.5, which had previously been identified as a SNR, is rather an HII region and has a flat radio spectrum.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Follow-up observations at 16 and 33 GHz of extragalactic sources from WMAP 3-year data: I - Spectral properties

We present follow-up observations of 97 point sources from the Wilkinson Microwave Anisotropy Probe (WMAP) 3-year data, contained within the New Extragalactic WMAP Point Source (NEWPS) catalogue between declinations of -4 and +60 degrees; the sources form a flux-density-limited sample complete to 1.1 Jy (approximately 5 sigma) at 33 GHz. Our observations were made at 16 GHz using the Arcminute Microkelvin Imager (AMI) and at 33 GHz with the Very Small Array (VSA). 94 of the sources have reliable, simultaneous -- typically a few minutes apart -- observations with both telescopes. The spectra between 13.9 and 33.75 GHz are very different from those of bright sources at low frequency: 44 per cent have rising spectra (alpha < 0.0), where flux density is proportional to frequency^-alpha, and 93 per cent have spectra with alpha < 0.5; the median spectral index is 0.04. For the brighter sources, the agreement between VSA and WMAP 33-GHz flux densities averaged over sources is very good. However, for the fainter sources, the VSA tends to measure lower values for the flux densities than WMAP. We suggest that the main cause of this effect is Eddington bias arising from variability.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Follow-up observations at 16 and 33 GHz of extragalactic sources from WMAP 3-year data: II - Flux-density variability

Using the Arcminute Microkelvin Imager (AMI) at 16 GHz and the Very Small Array (VSA) at 33 GHz to make follow-up observations of sources in the New Extragalactic WMAP Point Source Catalogue, we have investigated the flux-density variability in a complete sample of 97 sources over timescales of a few months to approximately 1.5 years. We find that 53 per cent of the 93 sources, for which we have multiple observations, are variable, at the 99 per cent confidence level, above the flux density calibration uncertainties of approximately 4 per cent at 16 GHz; the fraction of sources having varied by more than 20 per cent is 15 per cent at 16 GHz and 20 per cent at 33 GHz. Not only is this common occurrence of variability at high frequency of interest for source physics, but strategies for coping with source contamination in CMB work must take this variability into account. There is no strong evidence of a correlation between variability and flux density for the sample as a whole. Using a maximum-likelihood method, we calculate the correlation in the variability at the two frequencies in a subset of sources classified as variable from both the AMI and VSA data and find the Pearson product-moment correlation coefficient to be very high (0.955 +/- 0.034). We also find the degree of variability at 16 GHz (0.202 +/- 0.028) to be very similar to that at 33 GHz (0.224 +/- 0.039). Finally, we have investigated the relationship between variability and spectral index, alpha_13.9^33.75 (where S proportional to nu^-alpha), and find a significant difference in the spectral indices of the variable sources (-0.06 +/- 0.05) and non-variable sources (0.13 +/- 0.04).

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

G64.5+0.9, a new shell supernova remnant with unusual central emission

We present observations between 1.4 and 18 GHz confirming that G64.5+0.9 is new Galactic shell supernova remnant, using the Very Large Array and the Arcminute Microkelvin Imager. The remnant is a shell ~8 arcmin in diameter with a spectral index of alpha = 0.47 +/- 0.03 (with alpha defined such that flux density S varies with frequency nu as S proportional to nu to the power of -alpha). There is also emission near the centre of the shell, ~1 arcmin in extent, with a spectral index of alpha = 0.81 +/- 0.02. We do not find any evidence for spectral breaks for either source within our frequency range. The nature of the central object is unclear and requires further investigation, but we argue that is most unlikely to be extragalactic. It is difficult to avoid the conclusion that it is associated with the shell, although its spectrum is very unlike that of known pulsar wind nebulae.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Use of the MultiNest algorithm for gravitational wave data analysis

We describe an application of the MultiNest algorithm to gravitational wave data analysis. MultiNest is a multimodal nested sampling algorithm designed to efficiently evaluate the Bayesian evidence and return posterior probability densities for likelihood surfaces containing multiple secondary modes. The algorithm employs a set of live points which are updated by partitioning the set into multiple overlapping ellipsoids and sampling uniformly from within them. This set of live points climbs up the likelihood surface through nested iso-likelihood contours and the evidence and posterior distributions can be recovered from the point set evolution. The algorithm is model-independent in the sense that the specific problem being tackled enters only through the likelihood computation, and does not change how the live point set is updated. In this paper, we consider the use of the algorithm for gravitational wave data analysis by searching a simulated LISA data set containing two non-spinning supermassive black hole binary signals. The algorithm is able to rapidly identify all the modes of the solution and recover the true parameters of the sources to high precision.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

VSA Observations of the Anomalous Microwave Emission in the Perseus Region

The dust feature G159.6--18.5 in the Perseus region has previously been observed with the COSMOSOMAS experiment \citep{Watson:05} on angular scales of $\approx$ 1$^{\circ}$, and was found to exhibit anomalous microwave emission. We present new observations of this dust feature, performed with the Very Small Array (VSA) at 33 GHz, to help increase the understanding of the nature of this anomalous emission. On the angular scales observed with the VSA ($\approx$ 10 -- 40$^{\prime}$), G159.6--18.5 consists of five distinct components, each of which have been individually analysed. All five of these components are found to exhibit an excess of emission at 33 GHz, and are found to be highly correlated with far-infrared emission. We provide evidence that each of these compact components have anomalous emission that is consistent with electric dipole emission from very small, rapidly rotating dust grains. These components contribute $\approx$ 10 % to the flux density of the diffuse extended emission detected by COSMOSOMAS, and are found to have a similar radio emissivity.

0 citations0 reviewsNo code linkedOpen access
preprint2008arXiv

Observations of the Corona Borealis supercluster with the superextended Very Small Array: further constraints on the nature of the non-Gaussian CMB cold spot

We present interferometric imaging at 33 GHz, with the new superextended configuration of the Very Small Array (VSA), of a very deep decrement in the cosmic microwave background (CMB) temperature. This decrement is located in the direction of the Corona Borealis supercluster, at a position with no known galaxy clusters, and was discovered by a previous VSA survey (Genova-Santos et al.). A total area of 3 sq.deg. has now been imaged, with an angular resolution of 7 arcmin and a flux sensitivity of 5 mJy/beam. These observations confirm the presence of this strong and resolved negative spot at -41+/-5 mJy/beam (-258+/-29 muK), with a signal to noise level of 8. This structure is also present in the WMAP 5-year data. The temperature of the W-band (94 GHz) data at the position of the decrement agrees within 1.2-sigma with that observed by the VSA at 33 GHz, and within 0.2-sigma with the Sunyaev-Zel'dovich (SZ) spectrum. Our analyses show that it is a non-Gaussian feature in the CMB at a level of 4.8-sigma. The probability of finding such a deviation or larger in CMB Gaussian simulations is only 0.19 per cent. Therefore, an explanation other than primordial CMB is required. We have considered the possibility of an SZ effect generated in a diffuse, extended warm/hot gas distribution. This hypothesis is especially relevant, as the presence of such structures, if confirmed, could provide the location for a significant fraction of the missing baryons in the Local Universe.

0 citations0 reviewsNo code linkedOpen access
preprint2008arXiv

Radio source calibration for the VSA and other CMB instruments at around 30 GHz

Accurate calibration of data is essential for the current generation of CMB experiments. Using data from the Very Small Array (VSA), we describe procedures which will lead to an accuracy of 1 percent or better for experiments such as the VSA and CBI. Particular attention is paid to the stability of the receiver systems, the quality of the site and frequent observations of reference sources. At 30 GHz the careful correction for atmospheric emission and absorption is shown to be essential for achieving 1 percent precision. The sources for which a 1 percent relative flux density calibration was achieved included Cas A, Cyg A, Tau A and NGC7027 and the planets Venus, Jupiter and Saturn. A flux density, or brightness temperature in the case of the planets, was derived at 33 GHz relative to Jupiter which was adopted as the fundamental calibrator. A spectral index at ~30 GHz is given for each. Cas A,Tau A, NGC7027 and Venus were examined for variability. Cas A was found to be decreasing at $0.394 \pm 0.019$ percent per year over the period March 2001 to August 2004. In the same period Tau A was decreasing at $0.22\pm 0.07$ percent per year. A survey of the published data showed that the planetary nebula NGC7027 decreased at $0.16\pm 0.04$ percent per year over the period 1967 to 2003. Venus showed an insignificant ($1.5 \pm 1.3$ percent) variation with Venusian illumination. The integrated polarization of Tau A at 33 GHz was found to be $7.8\pm 0.6$ percent at pa $ = 148^\circ \pm 3^\circ$.}

0 citations0 reviewsNo code linkedOpen access