Toward Early-Warning Detection of Gravitational Waves from Compact Binary Coalescence

Rapid detection of compact binary coalescence (CBC) with a network of advanced gravitational-wave detectors will offer a unique opportunity for multi-messenger astronomy. Prompt detection alerts for the astronomical community might make it possible to observe the onset of electromagnetic emission from (CBC). We demonstrate a computationally practical filtering strategy that could produce early-warning triggers before gravitational radiation from the final merger has arrived at the detectors.Comment: 16 pages, 7 figures, published in ApJ. Reformatted preprint with emulateap

Efficient and Accurate Optimal Linear Phase FIR Filter Design Using Opposition-Based Harmony Search Algorithm

In this paper, opposition-based harmony search has been applied for the optimal design of linear phase FIR filters. RGA, PSO, and DE have also been adopted for the sake of comparison. The original harmony search algorithm is chosen as the parent one, and opposition-based approach is applied. During the initialization, randomly generated population of solutions is chosen, opposite solutions are also considered, and the fitter one is selected as a priori guess. In harmony memory, each such solution passes through memory consideration rule, pitch adjustment rule, and then opposition-based reinitialization generation jumping, which gives the optimum result corresponding to the least error fitness in multidimensional search space of FIR filter design. Incorporation of different control parameters in the basic HS algorithm results in the balancing of exploration and exploitation of search space. Low pass, high pass, band pass, and band stop FIR filters are designed with the proposed OHS and other aforementioned algorithms individually for comparative optimization performance. A comparison of simulation results reveals the optimization efficacy of the OHS over the other optimization techniques for the solution of the multimodal, nondifferentiable, nonlinear, and constrained FIR filter design problems

Time-delay interferometric ranging for LISA: Statistical analysis of bias-free ranging using laser noise minimization

Die Laser Interferometer Space Antenna (LISA) ist eine Mission der europäischen Weltraumagentur (ESA) zur Detektion von Gravitationswellen im Frequenzbereich zwischen 10^-4 Hz und 1 Hz. Gravitationswellen induzieren relative Abstandsänderungen, die LISA mithilfe von Laserinterferometrie mit Picometerpräzision misst. Ein großes Problem hierbei ist das Frequenzrauschen der Laser. Um dieses zu unterdrücken, ist es notwendig, mithilfe eines Algorithmus namens TDI (engl. time-delay interferometry), virtuelle Interferometer mit gleichlangen Armen zu konstruieren, wie z.B. das klassische Michelson-Interferometer. In dieser Arbeit untersuchen wir die Performanz von TDI unter realistischen Bedingungen und identifizieren verschiedene Kopplungsmechanismen des Laserfrequenzrauschens. Als erstes betrachten wir die Datenverarbeitung an Bord der Satelliten, die benötigt wird, um die Abtastrate der interferometrischen Messungen zu reduzieren. Hierfür sind Anti-Alias-Filter vorgesehen, die der Faltung von Laserrauschleistung in das Beobachtungsband vorbeugen. Außerdem wirkt sich die Ebenheit der Filter auf die Effektivität von TDI aus (engl. flexing-filtering-effect). Dieser Effekt ist bereits in der Literatur beschrieben und wir demonstrieren in dieser Arbeit die Möglichkeit, ihn mithilfe von Kompensationsfiltern effektiv zu reduzieren. Als zweites betrachten wir Kopplungsmechanismen von Laserfrequenzrauschen im TDI-Algorithmus selbst. Fehler in der Interpolation der interferometrischen Messungen und Ungenauigkeiten in den absoluten Abstandsmessungen zwischen den Satelliten führen ebenfalls zu einer unzureichenden Reduzierung des Laserfrequenzrauschens. Wir beschreiben die oben genannten Kopplungsmechanismen analytisch und validieren die zugrundeliegenden Modelle mithilfe von numerischen Simulationen. Das tiefere Verständnis dieser Residuen ermöglicht es uns, geeignete instrumentelle Parameter zu wählen, die von hoher Relevanz für das Missionsdesign von LISA sind. Des Weiteren beschäftigen wir uns in dieser Arbeit mit der möglichst genauen Bestimmung der absoluten Abständen zwischen den Satelliten, die für den TDI Algorithmus erforderlich sind. Hierfür werden die Abstandsinformationen aus den Seitenbändern und der PRN-Modulation (engl. pseudo-random noise) kombiniert. Wir zeigen, dass die PRN-Messung von systematischen Verzerrungen betroffen ist, die zu Laserrauschresiduen in den TDI-Variablen führen. Um diesen Fehler zu korrigieren, schlagen wir als zusätzliche Abstandsmessung TDI-Ranging (TDI-R) vor. TDI-R ist zwar ungenauer, aber frei von systematischen Verzerrungen und kann daher zur Kalibrierung der PRN-Messungen herangezogen werden. Wir präsentieren in dieser Arbeit eine ausführliche statistische Studie, um die Performanz von TDI-R zu charakterisieren. Dafür formulieren wir die Likelihood-Funktion der interferometrischen Messungen und berechnen die Fisher-Informationsmatrix, um die theoretisch mögliche untere Grenze der Schätzvarianz zu finden. Diese verhält sich invers proportional zur Integrationszeit und dem Verhältnis von Sekundärrauschleistung, die die interferometrische Messung fundamental limitiert, und Laserrauschleistung. Zusätzlich validieren wir die analytische untere Grenze der Schätzvarianz mithilfe von numerischen Simulationen und zeigen damit, dass unsere Implementierung von TDI-R optimal ist. Der entwickelte TDI-R-Algorithmus wird Teil der Datenverarbeitungspipeline sein und Konsistenzprüfungen und Kalibrierung der primären Abstandsmessmethoden ermöglichen.The Laser Interferometer Space Antenna (LISA) is a future ESA-led space-based observatory to explore the gravitational universe in the frequency band between 10^-4 Hz and 1 Hz. LISA implements picometer-precise inter-satellite ranging to measure tiny ripples in spacetime induced by gravitational waves (GWs). However, the single-link measurements are dominated by laser frequency noise, which is about nine orders of magnitude larger than the GW signals. Therefore, in post-processing, the time-delay interferometry (TDI) algorithm is used to synthesize virtual equal-arm interferometers to suppress laser frequency noise. In this work we identify several laser frequency noise coupling channels that limit the performance of TDI. First, the on-board processing, which is used to decimate the sampling rate from tens of megahertz down to the telemetry rate of a few hertz, requires careful design. Appropriate anti-aliasing filters must be implemented to mitigate folding of laser noise power into the observation band. Furthermore, the flatness of these filters is important to limit the impact of the flexing-filtering effect. We demonstrate that this effect can be effectively reduced by using compensation filters on ground. Second, the post-processing delays applied in TDI are subject to interpolation and ranging errors. We study these laser and timing noise residuals analytically and perform simulations to validate the models numerically. Our findings have direct implications for the design of the LISA instrument as we identify the instrumental parameters that are essential for successful laser noise suppression and provide methods for designing appropriate filters for the on-board processing. In addition, we discuss a dedicated ranging processing pipeline that produces high-precision range estimates that are the input for TDI by combining the sideband and pseudo-random noise (PRN) ranges. We show in this thesis that biases in the PRN measurements limit the laser noise suppression performance. Therefore, we propose time-delay interferometric ranging (TDI-R) as a third ranging sensor to estimate bias-free ranges that can be used to calibrate the biases in the PRN measurements. We present a thorough statistical study of TDI-R to evaluate its performance. Therefore, we formulate the likelihood function of the interferometric data and use the Fisher information formalism to find a lower bound on the estimation variance of the inter-satellite ranges. We find that the ranging uncertainty is proportional to the inverse of the integration time and the ratio of secondary noise power, that limits the interferometric readout, to the laser noise power. To validate our findings we implement prototype TDI-R pipelines and perform numerical simulations. We show that we are able to formulate optimal estimators of the unbiased range that reach the Cramér-Rao lower bound previously expressed analytically. The developed TDI-R pipeline will be integrated into the ranging processing pipeline to perform consistency checks and ensure well-calibrated inter-satellite ranges

BICEP2 II: Experiment and Three-Year Data Set

We report on the design and performance of the BICEP2 instrument and on its three-year data set. BICEP2 was designed to measure the polarization of the cosmic microwave background (CMB) on angular scales of 1 to 5 degrees ($\ell$=40-200), near the expected peak of the B-mode polarization signature of primordial gravitational waves from cosmic inflation. Measuring B-modes requires dramatic improvements in sensitivity combined with exquisite control of systematics. The BICEP2 telescope observed from the South Pole with a 26~cm aperture and cold, on-axis, refractive optics. BICEP2 also adopted a new detector design in which beam-defining slot antenna arrays couple to transition-edge sensor (TES) bolometers, all fabricated on a common substrate. The antenna-coupled TES detectors supported scalable fabrication and multiplexed readout that allowed BICEP2 to achieve a high detector count of 500 bolometers at 150 GHz, giving unprecedented sensitivity to B-modes at degree angular scales. After optimization of detector and readout parameters, BICEP2 achieved an instrument noise-equivalent temperature of 15.8 $\mu$K sqrt(s). The full data set reached Stokes Q and U map depths of 87.2 nK in square-degree pixels (5.2 $\mu$K arcmin) over an effective area of 384 square degrees within a 1000 square degree field. These are the deepest CMB polarization maps at degree angular scales to date. The power spectrum analysis presented in a companion paper has resulted in a significant detection of B-mode polarization at degree scales.Comment: 30 pages, 24 figure

Design of optimal band-stop FIR filter using L1-norm based RCGA

AbstractIn this paper, an optimal design of linear phase digital finite impulse response (FIR) band stop (BS) filter using the L1-norm based real-coded genetic algorithm (L1-RCGA) is presented. Although RCGA has proved its ability to overcome the drawbacks associated with conventional gradient-based optimization methods of filter design, it is applied here with a novel fitness function based on the L1-norm. This leads to a global optimal solution along with the improvement in filter design with same specifications. The designed filter pursues a better response in terms of flat passband, high stopband attenuation and fast convergence. The simulation results justify that the proposed FIR BS filter using L1-RCGA outperforms the existing optimization techniques, the L1-method and particle swarm optimization (PSO) and the conventional methods such as least-squares (LS) approach, Kaiser window method and the Parks McClellan (PM) algorithm. A detailed analysis is performed to evaluate the performance of the designed filters

Optimizing Advanced Ligo\u27s Scientific Output with Fast, Accurate, Clean Calibration

Since 2015, the direct observation of gravitational waves has opened a new window to observe the universe and made strong-field tests of Einstein\u27s general theory of relativity possible for the first time. During the first two observing runs of the Advanced gravitational-wave detector network, the Laser Interferometer Gravitational-wave Observatory (LIGO) and the Virgo detector have made 10 detections of binary black hole mergers and one detection of a binary neutron star merger with a coincident gamma-ray burst. This dissertation discusses methods used in low and high latency to produce Advanced LIGO\u27s calibrated strain data, highlighting improvements to accuracy, latency, and noise reduction that have been made since the beginning of the second observing run (O2). Systematic errors in the calibration during O2 varied by frequency, but were generally no greater that 5% in amplitude and 3 deg in phase from 20 Hz to 1 kHz. Due in part to this work, it is now possible to achieve calibration accuracy at the level of ~1% in amplitude and ~1 deg in phase, offering improvements to downstream astrophysical analyses. Since the beginning of O2, latency intrinsic to the calibration procedure has decreased from ~12 s to ~3 s. As latency in data distribution and the sending of automated alerts to astronomers is minimized, reduction in calibration latency will become important for follow-up of events like the binary neutron star merger GW170817. A method of removing spectral lines and broadband noise in the calibration procedure has been developed since O2, offering increases in total detectable volume during future observing runs. High-latency subtraction of lines and broadband noise had a large impact on astrophysical analyses during O2. A similar data product can now be made available in low latency for the first time

The balloon-borne large-aperture submillimeter telescope for polarimetry: BLAST-Pol

The Balloon-borne Large Aperture Submillimeter Telescope for Polarimetry (BLAST-Pol) is a suborbital mapping experiment designed to study the role played by magnetic fields in the star formation process. BLAST-Pol is the reconstructed BLAST telescope, with the addition of linear polarization capability. Using a 1.8 m Cassegrain telescope, BLAST-Pol images the sky onto a focal plane that consists of 280 bolometric detectors in three arrays, observing simultaneously at 250, 350, and 500 um. The diffraction-limited optical system provides a resolution of 30'' at 250 um. The polarimeter consists of photolithographic polarizing grids mounted in front of each bolometer/detector array. A rotating 4 K achromatic half-wave plate provides additional polarization modulation. With its unprecedented mapping speed and resolution, BLAST-Pol will produce three-color polarization maps for a large number of molecular clouds. The instrument provides a much needed bridge in spatial coverage between larger-scale, coarse resolution surveys and narrow field of view, and high resolution observations of substructure within molecular cloud cores. The first science flight will be from McMurdo Station, Antarctica in December 2010.Comment: 14 pages, 9 figures Submitted to SPIE Astronomical Telescopes and Instrumentation Conference 201

A survey on tidal analysis and forecasting methods for Tsunami detection

Accurate analysis and forecasting of tidal level are very important tasks for human activities in oceanic and coastal areas. They can be crucial in catastrophic situations like occurrences of Tsunamis in order to provide a rapid alerting to the human population involved and to save lives. Conventional tidal forecasting methods are based on harmonic analysis using the least squares method to determine harmonic parameters. However, a large number of parameters and long-term measured data are required for precise tidal level predictions with harmonic analysis. Furthermore, traditional harmonic methods rely on models based on the analysis of astronomical components and they can be inadequate when the contribution of non-astronomical components, such as the weather, is significant. Other alternative approaches have been developed in the literature in order to deal with these situations and provide predictions with the desired accuracy, with respect also to the length of the available tidal record. These methods include standard high or band pass filtering techniques, although the relatively deterministic character and large amplitude of tidal signals make special techniques, like artificial neural networks and wavelets transform analysis methods, more effective. This paper is intended to provide the communities of both researchers and practitioners with a broadly applicable, up to date coverage of tidal analysis and forecasting methodologies that have proven to be successful in a variety of circumstances, and that hold particular promise for success in the future. Classical and novel methods are reviewed in a systematic and consistent way, outlining their main concepts and components, similarities and differences, advantages and disadvantages

BICEP2. II. Experiment and three-year Data Set

We report on the design and performance of the BICEP2 instrument and on its three-year data set. BICEP2 was designed to measure the polarization of the cosmic microwave background (CMB) on angular scales of 1°-5°(ℓ = 40-200), near the expected peak of the B-mode polarization signature of primordial gravitational waves from cosmic inflation. Measuring B-modes requires dramatic improvements in sensitivity combined with exquisite control of systematics. The BICEP2 telescope observed from the South Pole with a 26 cm aperture and cold, on-axis, refractive optics. BICEP2 also adopted a new detector design in which beam-defining slot antenna arrays couple to transition-edge sensor (TES) bolometers, all fabricated on a common substrate. The antenna-coupled TES detectors supported scalable fabrication and multiplexed readout that allowed BICEP2 to achieve a high detector count of 500 bolometers at 150 GHz, giving unprecedented sensitivity to B-modes at degree angular scales. After optimization of detector and readout parameters, BICEP2 achieved an instrument noise-equivalent temperature of 15.8 µK√s. The full data set reached Stokes Q and U map depths of 87.2 nK in square-degree pixels (5.'2 μK) over an effective area of 384 deg^2 within a 1000 deg^2 field. These are the deepest CMB polarization maps at degree angular scales to date. The power spectrum analysis presented in a companion paper has resulted in a significant detection of B-mode polarization at degree scales