Sanjit Mitra

Probing the anisotropies of a stochastic gravitational-wave background using a network of ground-based laser interferometers

We present a maximum-likelihood analysis for estimating the angular distribution of power in an anisotropic stochastic gravitational-wave background using ground-based laser interferometers. The standard isotropic and gravitational-wave radiometer searches (optimal for point sources) are recovered as special limiting cases. The angular distribution can be decomposed with respect to any set of basis functions on the sky, and the single-baseline, cross-correlation analysis is easily extended to a network of three or more detectors--that is, to multiple baselines. A spherical-harmonic decomposition, which provides maximum-likelihood estimates of the multipole moments of the gravitational-wave sky, is described in detail. We also discuss (i) the covariance matrix of the estimators and its relationship to the detector response of a network of interferometers, (ii) a singular-value decomposition method for regularizing the deconvolution of the detector response from the measured sky map, (iii) the expected increase in sensitivity obtained by including multiple baselines, and (iv) the numerical results of this method when applied to simulated data consisting of both pointlike and diffuse sources. Comparisons between this general method and the standard isotropic and radiometer searches are given throughout, to make contact with the existing literature on stochastic background searches.

Gravitational wave radiometry : Mapping a stochastic gravitational wave background

The problem of the detection and mapping of a stochastic gravitational wave background (SGWB), either cosmological or astrophysical, bears a strong semblance to the analysis of the cosmic microwave background (CMB) anisotropy and polarization, which too is a stochastic field, statistically described in terms of its correlation properties. An astrophysical gravitational wave background (AGWB) will likely arise from an incoherent superposition of unmodelled and/or unresolved sources and cosmological gravitational wave backgrounds (CGWB) are also predicted in certain scenarios. The basic statistic we use is the cross correlation between the data from a pair of detectors. In order to “point” the pair of detectors at different locations one must suitably delay the signal by the amount it takes for the gravitational waves (GW) to travel to both detectors corresponding to a source direction. Then the raw (observed) sky map of the SGWB is the signal convolved with a beam response function that varies with location in the sky. We first present a thorough analytic understanding of the structure of the beam response function using an analytic approach employing the stationary phase approximation. The true sky map is obtained by numerically deconvolving the beam function in the integral (convolution) equation. We adopt the maximum likelihood framework to estimate the true sky map using the conjugate gradient method that has been successfully used in the broadly similar, well-studied CMB map-making problem. We numerically implement and demonstrate the method on signal generated by simulated (unpolarized) SGWB for the GW radiometer consisting of the LIGO pair of detectors at Hanford and Livingston. We include “realistic” additive Gaussian noise in each data stream based on the LIGO-I noise power spectral density. The extension of the method to multiple baselines and polarized GWB is outlined. In the near future the network of GW detectors, including the Advanced LIGO and Virgo detectors that will be sensitive to sources within a thousand times larger spatial volume, could provide promising data sets for GW radiometry.

Non-circular beam correction to the CMB power spectrum

Abstract In the era of high precision CMB measurements, systematic effects are beginning to limit the ability to extract subtler cosmological information. The non-circularity of the experimental beam has become progressively important as CMB experiments strive to attain higher angular resolution and sensitivity. The effect of non-circular beam on the power spectrum is important at multipoles larger than the beam-width. For recent experiments with high angular resolution, optimal methods of power spectrum estimation are computationally prohibitive and sub-optimal approaches, such as the Pseudo- C l method are used. We provide an analytic framework for correcting the power spectrum for the effect of beam non-circularity and non-uniform sky coverage (including incomplete/masked sky maps). The approach is perturbative in the distortion of the beam from non-circularity allowing for rapid computations when the beam is mildly non-circular. We advocate that when the non-circular beams are important, it is computationally advantageous to employ ‘soft’ azimuthally apodized masks whose spherical harmonic transforms die down fast with m .

Towards a first design of a Newtonian-noise cancellation system for Advanced LIGO

Newtonian gravitational noise from seismic fields is predicted to be a limiting noise source at low frequency for second generation gravitational-wave detectors. Mitigation of this noise will be achieved by Wiener filtering using arrays of seismometers deployed in the vicinity of all test masses. In this work, we present optimized configurations of seismometer arrays using a variety of simplified models of the seismic field based on seismic observations at LIGO Hanford. The model that best fits the seismic measurements leads to noise reduction limited predominantly by seismometer self-noise. A first simplified design of seismic arrays for Newtonian-noise cancellation at the LIGO sites is presented, which suggests that it will be sufficient to monitor surface displacement inside the buildings.

Orthogonal bipolar spherical harmonics measures: Scrutinizing sources of isotropy violation

The two-point correlation function of the cosmic microwave background temperature anisotropies is generally assumed to be statistically isotropic (SI). Deviations from this assumption could be traced to physical or observational artifacts and systematic effects. Measurement of nonvanishing power in the bipolar spherical harmonic spectra is a standard statistical technique to search for isotropy violations. Although this is a neat tool allowing a blind search for SI violations in the cosmic microwave background sky, it is not easy to discern the cause of isotropy violation by using this measure. In this article, we propose a novel technique of constructing orthogonal bipolar spherical harmonic estimators, which can be used to discern between models of isotropy violation.

Effect of noncircularity of experimental beam on CMB parameter estimation

Measurement of Cosmic Microwave Background (CMB) anisotropies has been playing a lead role in precision cosmology by providing some of the tightest constrains on cosmological models and parameters. However, precision can only be meaningful when all major systematic effects are taken into account. Non-circular beams in CMB experiments can cause large systematic deviation in the angular power spectrum, not only by modifying the measurement at a given multipole, but also introducing coupling between different multipoles through a deterministic bias matrix. Here we add a mechanism for emulating the effect of a full bias matrix to the Planck likelihood code through the parameter estimation code SCoPE. We show that if the angular power spectrum was measured with a non-circular beam, the assumption of circular Gaussian beam or considering only the diagonal part of the bias matrix can lead to huge error in parameter estimation. We demonstrate that, at least for elliptical Gaussian beams, use of scalar beam window functions obtained via Monte Carlo simulations starting from a fiducial spectrum, as implemented in Planck analyses for example, leads to em only few percent of sigma deviation of the best-fit parameters. However, we notice more significant differences in the posterior distributions for some of the parameters, which would in turn lead to incorrect errorbars. These differences can be reduced, so that the errorbars match within few percent, by adding an iterative reanalysis step, where the beam window function would be recomputed using the best-fit spectrum estimated in the first step.

Statistical isotropy violation in WMAP CMB maps resulting from non-circular beams

Statistical isotropy (SI) of cosmic microwave background (CMB) fluctuations is a key observational test to validate the cosmological principle underlying the standard model of cosmology. While a detection of SI violation would have immense cosmological ramification, it is important to recognise their possible origin in systematic effects of observations. The WMAP seven year (WMAP-7) release claimed significant deviation from SI in the bipolar spherical harmonic (BipoSH) coefficients A ll 20 and A l -2 l 20 >. Here we present the first explicit reproduction of the measurements reported in WMAP-7, confirming that beam systematics alone can completely account for the measured SI violation. The possibility of such a systematic origin was alluded to in WMAP-7 paper itself and other authors but not as explicitly so as to account for it accurately. We simulate CMB maps using the actual WMAP non-circular beams and scanning strategy. Our estimated BipoSH spectra from these maps match the WMAP-7 results very well. It is also evident that only a very careful and adequately detailed modelling, as carried out here, can conclusively establish that the entire signal arises from non-circular beam effect. This is important since cosmic SI violation signals are expected to be subtle and dismissing a large SI violation signal as observational artefact based on simplistic plausibility arguments run the serious risk of “throwing the baby out with the bathwater”.

All-sky, narrowband, gravitational-wave radiometry with folded data

Gravitational-wave radiometry is a powerful tool by which weak signals with unknown signal morphologies are recovered through a process of cross correlation. Radiometry has been used, e.g., to search for persistent signals from known neutron stars such as Scorpius X-1. In this paper, we demonstrate how a more ambitious search--for persistent signals from unknown neutron stars--can be efficiently carried out using folded data, in which an entire ~year-long observing run is represented as a single sidereal day. The all-sky, narrowband radiometer search described here will provide a computationally tractable means to uncover gravitational-wave signals from unknown, nearby neutron stars in binary systems, which can have modulation depths of ~0.1-2 Hz. It will simultaneously provide a sensitive search algorithm for other persistent, narrowband signals from unexpected sources.

Fast gravitational wave radiometry using data folding

Gravitational Waves (GWs) from the early universe and unresolved astrophysical sources are expected to create a stochastic GW background (SGWB). The GW radiometer algorithm is well suited to probe such a background using data from ground based laser interferometric detectors. Radiometer analysis can be performed in different bases, e.g., isotropic, pixel or spherical harmonic. Each of these analyses possesses a common temporal symmetry which we exploit here to fold the whole dataset for every detector pair, typically a few hundred to a thousand days of data, to only one sidereal day, without any compromise in precision. We develop the algebra and a software pipeline needed to fold data, accounting for the effect of overlapping windows and non-stationary noise. We implement this on LIGOs fifth science run data and validate it by performing a standard anisotropic SGWB search on both folded and unfolded data. Folded data not only leads to orders of magnitude reduction in computation cost, but it results in a conveniently small data volume of few gigabytes, making it possible to perform an actual analysis on a personal computer, as well as easy movement of data. A few important analyses, yet unaccomplished due to computational limitations, will now become feasible. Folded data, being independent of the radiometer basis, will also be useful in reducing processing redundancies in multiple searches and provide a common ground for mutual consistency checks. Most importantly, folded data will allow vast amount of experimentation with existing searches and provide substantial help in developing new strategies to find unknown sources.

Estimating statistical isotropy violation in CMB due to non-circular beam and complex scan in minutes

Mild, unavoidable deviations from circular-symmetry of instrumental beams along with scan strategy can give rise to measurable Statistical Isotropy (SI) violation in Cosmic Microwave Background (CMB) experiments. If not accounted properly, this spurious signal can complicate the extraction of other SI violation signals (if any) in the data. However, estimation of this effect through exact numerical simulation is computationally intensive and time consuming. A generalized analytical formalism not only provides a quick way of estimating this signal, but also gives a detailed understanding connecting the leading beam anisotropy components to a measurable BipoSH characterisation of SI violation. In this paper, we provide an approximate generic analytical method for estimating the SI violation generated due to a non-circular (NC) beam and arbitrary scan strategy, in terms of the Bipolar Spherical Harmonic (BipoSH) spectra. Our analytical method can predict almost all the features introduced by a NC beam in a complex scan and thus reduces the need for extensive numerical simulation worth tens of thousands of CPU hours into minutes long calculations. As an illustrative example, we use WMAP beams and scanning strategy to demonstrate the easability, usability and efficiency of our method. We test all our analytical results against that from exact numerical simulations.