Probing dark energy with steerable wavelets through correlation of WMAP and NVSS local morphological measures

Using local morphological measures on the sphere defined through a steerable wavelet analysis, we examine the three-year WMAP and the NVSS data for correlation induced by the integrated Sachs-Wolfe (ISW) effect. The steerable wavelet constructed from the second derivative of a Gaussian allows one to define three local morphological measures, namely the signed-intensity, orientation and elongation of local features. Detections of correlation between the WMAP and NVSS data are made with each of these morphological measures. The most significant detection is obtained in the correlation of the signed-intensity of local features at a significance of 99.9%. By inspecting signed-intensity sky maps, it is possible for the first time to see the correlation between the WMAP and NVSS data by eye. Foreground contamination and instrumental systematics in the WMAP data are ruled out as the source of all significant detections of correlation. Our results provide new insight on the ISW effect by probing the morphological nature of the correlation induced between the cosmic microwave background and large scale structure of the Universe. Given the current constraints on the flatness of the Universe, our detection of the ISW effect again provides direct and independent evidence for dark energy. Moreover, this new morphological analysis may be used in future to help us to better understand the nature of dark energy.Comment: 12 pages, 10 figures, replaced to match version accepted by MNRA

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

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

Classifying LISA gravitational wave burst signals using Bayesian evidence

We consider the problem of characterisation of burst sources detected with the Laser Interferometer Space Antenna (LISA) using the multi-modal nested sampling algorithm, MultiNest. We use MultiNest as a tool to search for modelled bursts from cosmic string cusps, and compute the Bayesian evidence associated with the cosmic string model. As an alternative burst model, we consider sine-Gaussian burst signals, and show how the evidence ratio can be used to choose between these two alternatives. We present results from an application of MultiNest to the last round of the Mock LISA Data Challenge, in which we were able to successfully detect and characterise all three of the cosmic string burst sources present in the release data set. We also present results of independent trials and show that MultiNest can detect cosmic string signals with signal-to-noise ratio (SNR) as low as ~7 and sine-Gaussian signals with SNR as low as ~8. In both cases, we show that the threshold at which the sources become detectable coincides with the SNR at which the evidence ratio begins to favour the correct model over the alternative.Comment: 21 pages, 11 figures, accepted by CQG; v2 has minor changes for consistency with accepted versio

Computational science and re-discovery: open-source implementations of ellipsoidal harmonics for problems in potential theory

We present two open-source (BSD) implementations of ellipsoidal harmonic expansions for solving problems of potential theory using separation of variables. Ellipsoidal harmonics are used surprisingly infrequently, considering their substantial value for problems ranging in scale from molecules to the entire solar system. In this article, we suggest two possible reasons for the paucity relative to spherical harmonics. The first is essentially historical---ellipsoidal harmonics developed during the late 19th century and early 20th, when it was found that only the lowest-order harmonics are expressible in closed form. Each higher-order term requires the solution of an eigenvalue problem, and tedious manual computation seems to have discouraged applications and theoretical studies. The second explanation is practical: even with modern computers and accurate eigenvalue algorithms, expansions in ellipsoidal harmonics are significantly more challenging to compute than those in Cartesian or spherical coordinates. The present implementations reduce the "barrier to entry" by providing an easy and free way for the community to begin using ellipsoidal harmonics in actual research. We demonstrate our implementation using the specific and physiologically crucial problem of how charged proteins interact with their environment, and ask: what other analytical tools await re-discovery in an era of inexpensive computation?Comment: 25 pages, 3 figure

Dust heating by the interstellar radiation field in models of turbulent molecular clouds

We have calculated the radiation field, dust grain temperatures, and far infrared emissivity of numerical models of turbulent molecular clouds. When compared to a uniform cloud of the same mean optical depth, most of the volume inside the turbulent cloud is brighter, but most of the mass is darker. There is little mean attenuation from center to edge, and clumping causes the radiation field to be somewhat bluer. There is also a large dispersion, typically by a few orders of magnitude, of all quantities relative to their means. However, despite the scatter, the 850 micron emission maps are well correlated with surface density. The fraction of mass as a function of intensity can be reproduced by a simple hierarchical model of density structure.Comment: 32 pages, 14 figures, submitted to Ap

Characterizing Planetary Orbits and the Trajectories of Light

Exact analytic expressions for planetary orbits and light trajectories in the Schwarzschild geometry are presented. A new parameter space is used to characterize all possible planetary orbits. Different regions in this parameter space can be associated with different characteristics of the orbits. The boundaries for these regions are clearly defined. Observational data can be directly associated with points in the regions. A possible extension of these considerations with an additional parameter for the case of Kerr geometry is briefly discussed.Comment: 49 pages total with 11 tables and 10 figure

Non-invasive, near-field terahertz imaging of hidden objects using a single pixel detector

Terahertz (THz) imaging has the ability to see through otherwise opaque materials. However, due to the long wavelengths of THz radiation ({\lambda}=300{\mu}m at 1THz), far-field THz imaging techniques are heavily outperformed by optical imaging in regards to the obtained resolution. In this work we demonstrate near-field THz imaging with a single-pixel detector. We project a time-varying optical mask onto a silicon wafer which is used to spatially modulate a pulse of THz radiation. The far-field transmission corresponding to each mask is recorded by a single element detector and this data is used to reconstruct the image of an object placed on the far side of the silicon wafer. We demonstrate a proof of principal application where we image a printed circuit board on the underside of a 115{\mu}m thick silicon wafer with ~100{\mu}m ({\lambda}/4) resolution. With subwavelength resolution and the inherent sensitivity to local conductivity provided by the THz probe frequencies, we show that it is possible to detect fissures in the circuitry wiring of a few microns in size. Imaging systems of this type could have other uses where non-invasive measurement or imaging of concealed structures with high resolution is necessary, such as in semiconductor manufacturing or in bio-imaging

Faecal incontinence persisting after childbirth : a 12 year longitudinal study

© 2012 The Authors BJOG An International Journal of Obstetrics and Gynaecology © 2012 RCOG.Peer reviewedPostprin