Subtraction-noise projection in gravitational-wave detector networks

In this paper, we present a successful implementation of a subtraction-noise projection method into a simple, simulated data analysis pipeline of a gravitational-wave search. We investigate the problem to reveal a weak stochastic background signal which is covered by a strong foreground of compact-binary coalescences. The foreground which is estimated by matched filters, has to be subtracted from the data. Even an optimal analysis of foreground signals will leave subtraction noise due to estimation errors of template parameters which may corrupt the measurement of the background signal. The subtraction noise can be removed by a noise projection. We apply our analysis pipeline to the proposed future-generation space-borne Big Bang Observer (BBO) mission which seeks for a stochastic background of primordial GWs in the frequency range $\sim 0.1-1$Hz covered by a foreground of black-hole and neutron-star binaries. Our analysis is based on a simulation code which provides a dynamical model of a time-delay interferometer (TDI) network. It generates the data as time series and incorporates the analysis pipeline together with the noise projection. Our results confirm previous ad hoc predictions which say that BBO will be sensitive to backgrounds with fractional energy densities below $\Omega=10^{-16}$Comment: 54 pages, 15 figure

BAT - The Bayesian Analysis Toolkit

We describe the development of a new toolkit for data analysis. The analysis package is based on Bayes' Theorem, and is realized with the use of Markov Chain Monte Carlo. This gives access to the full posterior probability distribution. Parameter estimation, limit setting and uncertainty propagation are implemented in a straightforward manner. A goodness-of-fit criterion is presented which is intuitive and of great practical use.Comment: 31 pages, 10 figure

Characterization of the Crab Pulsar's Timing Noise

We present a power spectral analysis of the Crab pulsar's timing noise, mainly using radio measurements from Jodrell Bank taken over the period 1982-1989. The power spectral analysis is complicated by nonuniform data sampling and the presence of a steep red power spectrum that can distort power spectra measurement by causing severe power ``leakage''. We develop a simple windowing method for computing red noise power spectra of uniformly sampled data sets and test it on Monte Carlo generated sample realizations of red power-law noise. We generalize time-domain methods of generating power-law red noise with even integer spectral indices to the case of noninteger spectral indices. The Jodrell Bank pulse phase residuals are dense and smooth enough that an interpolation onto a uniform time series is possible. A windowed power spectrum is computed revealing a periodic or nearly periodic component with a period of about 568 days and a 1/f^3 power-law noise component with a noise strength of 1.24 +/- 0.067 10^{-16} cycles^2/sec^2 over the analysis frequency range 0.003 - 0.1 cycles/day. This result deviates from past analyses which characterized the pulse phase timing residuals as either 1/f^4 power-law noise or a quasiperiodic process. The analysis was checked using the Deeter polynomial method of power spectrum estimation that was developed for the case of nonuniform sampling, but has lower spectral resolution. The timing noise is consistent with a torque noise spectrum rising with analysis frequency as f implying blue torque noise, a result not predicted by current models of pulsar timing noise. If the periodic or nearly periodic component is due to a binary companion, we find a companion mass > 3.2 Earth masses.Comment: 53 pages, 9 figures, submitted to MNRAS, abstract condense

Bayesian power-spectrum inference for Large Scale Structure data

We describe an exact, flexible, and computationally efficient algorithm for a joint estimation of the large-scale structure and its power-spectrum, building on a Gibbs sampling framework and present its implementation ARES (Algorithm for REconstruction and Sampling). ARES is designed to reconstruct the 3D power-spectrum together with the underlying dark matter density field in a Bayesian framework, under the reasonable assumption that the long wavelength Fourier components are Gaussian distributed. As a result ARES does not only provide a single estimate but samples from the joint posterior of the power-spectrum and density field conditional on a set of observations. This enables us to calculate any desired statistical summary, in particular we are able to provide joint uncertainty estimates. We apply our method to mock catalogs, with highly structured observational masks and selection functions, in order to demonstrate its ability to reconstruct the power-spectrum from real data sets, while fully accounting for any mask induced mode coupling.Comment: 25 pages, 15 figure

Calculation of the Performance of Communication Systems from Measured Oscillator Phase Noise

Oscillator phase noise (PN) is one of the major problems that affect the performance of communication systems. In this paper, a direct connection between oscillator measurements, in terms of measured single-side band PN spectrum, and the optimal communication system performance, in terms of the resulting error vector magnitude (EVM) due to PN, is mathematically derived and analyzed. First, a statistical model of the PN, considering the effect of white and colored noise sources, is derived. Then, we utilize this model to derive the modified Bayesian Cramer-Rao bound on PN estimation, and use it to find an EVM bound for the system performance. Based on our analysis, it is found that the influence from different noise regions strongly depends on the communication bandwidth, i.e., the symbol rate. For high symbol rate communication systems, cumulative PN that appears near carrier is of relatively low importance compared to the white PN far from carrier. Our results also show that 1/f^3 noise is more predictable compared to 1/f^2 noise and in a fair comparison it affects the performance less.Comment: Accepted in IEEE Transactions on Circuits and Systems-I: Regular Paper

EChOSim: The Exoplanet Characterisation Observatory software simulator

EChOSim is the end-to-end time-domain simulator of the Exoplanet Characterisation Observatory (EChO) space mission. EChOSim has been developed to assess the capability EChO has to detect and characterize the atmospheres of transiting exoplanets, and through this revolutionize the knowledge we have of the Milky Way and of our place in the Galaxy. Here we discuss the details of the EChOSim implementation and describe the models used to represent the instrument and to simulate the detection. Software simulators have assumed a central role in the design of new instrumentation and in assessing the level of systematics affecting the measurements of existing experiments. Thanks to its high modularity, EChOSim can simulate basic aspects of several existing and proposed spectrometers for exoplanet transits, including instruments on the Hubble Space Telescope and Spitzer, or ground-based and balloon borne experiments. A discussion of different uses of EChOSim is given, including examples of simulations performed to assess the EChO mission

Source-intrinsic near-infrared properties of Sgr A*: Total intensity measurements

We present a comprehensive data description for Ks-band measurements of Sgr A*. We characterize the statistical properties of the variability of Sgr A* in the near-infrared, which we find to be consistent with a single-state process forming a power-law distribution of the flux density. We discover a linear rms-flux relation for the flux-density range up to 12 mJy on a timescale of 24 minutes. This and the power-law flux density distribution implies a phenomenological, formally non-linear statistical variability model with which we can simulate the observed variability and extrapolate its behavior to higher flux levels and longer timescales. We present reasons why data with our cadence cannot be used to decide on the question whether the power spectral density of the underlying random process shows more structure at timescales between 25 min and 100 min compared to what is expected from a red noise random process.Comment: Accepted to ApJS, August 27, 201

Report by the ESA-ESO Working Group on Fundamental Cosmology

ESO and ESA agreed to establish a number of Working Groups to explore possible synergies between these two major European astronomical institutions. This Working Group's mandate was to concentrate on fundamental questions in cosmology, and the scope for tackling these in Europe over the next ~15 years. One major resulting recommendation concerns the provision of new generations of imaging survey, where the image quality and near-IR sensitivity that can be attained only in space are naturally matched by ground-based imaging and spectroscopy to yield massive datasets with well-understood photometric redshifts (photo-z's). Such information is essential for a range of new cosmological tests using gravitational lensing, large-scale structure, clusters of galaxies, and supernovae. Great scope in future cosmology also exists for ELT studies of the intergalactic medium and space-based studies of the CMB and gravitational waves; here the synergy is less direct, but these areas will remain of the highest mutual interest to the agencies. All these recommended facilities will produce vast datasets of general applicability, which will have a tremendous impact on broad areas of astronomy.Comment: ESA-ESO Working Groups Report No. 3, 125 pages, 28 figures. A PDF version including the cover is available from http://www.stecf.org/coordination/esa_eso/cosmology/report_cover.pdf and a printed version (A5 booklet) is available in limited numbers from the Space Telescope-European Coordinating Facility (ST-ECF): [email protected]

High precision simulations of weak lensing effect on Cosmic Microwave Background polarization

We study accuracy, robustness and self-consistency of pixel-domain simulations of the gravitational lensing effect on the primordial CMB anisotropies due to the large-scale structure of the Universe. In particular, we investigate dependence of the results precision on some crucial parameters of such techniques and propose a semi-analytic framework to determine their values so the required precision is a priori assured and the numerical workload simultaneously optimized. Our focus is on the B-mode signal but we discuss also other CMB observables, such as total intensity, T, and E-mode polarization, emphasizing differences and similarities between all these cases. Our semi-analytic considerations are backed up by extensive numerical results. Those are obtained using a code, nicknamed lenS2HAT -- for Lensing using Scalable Spherical Harmonic Transforms (S2HAT) -- which we have developed in the course of this work. The code implements a version of the pixel-domain approach of Lewis (2005) and permits performing the simulations at very high resolutions and data volumes, thanks to its efficient parallelization provided by the S2HAT library -- a parallel library for a calculation of the spherical harmonic transforms. The code is made publicly available.Comment: 20 pages, 14 figures, submitted to A&A, matches version accepted for publication in A&