Learning Dictionaries with Bounded Self-Coherence

Sparse coding in learned dictionaries has been established as a successful approach for signal denoising, source separation and solving inverse problems in general. A dictionary learning method adapts an initial dictionary to a particular signal class by iteratively computing an approximate factorization of a training data matrix into a dictionary and a sparse coding matrix. The learned dictionary is characterized by two properties: the coherence of the dictionary to observations of the signal class, and the self-coherence of the dictionary atoms. A high coherence to the signal class enables the sparse coding of signal observations with a small approximation error, while a low self-coherence of the atoms guarantees atom recovery and a more rapid residual error decay rate for the sparse coding algorithm. The two goals of high signal coherence and low self-coherence are typically in conflict, therefore one seeks a trade-off between them, depending on the application. We present a dictionary learning method with an effective control over the self-coherence of the trained dictionary, enabling a trade-off between maximizing the sparsity of codings and approximating an equiangular tight frame.Comment: 4 pages, 2 figures; IEEE Signal Processing Letters, vol. 19, no. 12, 201

Bayesian inference on compact binary inspiral gravitational radiation signals in interferometric data

Presented is a description of a Markov chain Monte Carlo (MCMC) parameter estimation routine for use with interferometric gravitational radiational data in searches for binary neutron star inspiral signals. Five parameters associated with the inspiral can be estimated, and summary statistics are produced. Advanced MCMC methods were implemented, including importance resampling and prior distributions based on detection probability, in order to increase the efficiency of the code. An example is presented from an application using realistic, albeit fictitious, data.Comment: submitted to Classical and Quantum Gravity. 14 pages, 5 figure

Parameter estimation of spinning binary inspirals using Markov-chain Monte Carlo

We present a Markov-chain Monte-Carlo (MCMC) technique to study the source parameters of gravitational-wave signals from the inspirals of stellar-mass compact binaries detected with ground-based gravitational-wave detectors such as LIGO and Virgo, for the case where spin is present in the more massive compact object in the binary. We discuss aspects of the MCMC algorithm that allow us to sample the parameter space in an efficient way. We show sample runs that illustrate the possibilities of our MCMC code and the difficulties that we encounter.Comment: 10 pages, 2 figures, submitted to Classical and Quantum Gravit

Coherent Bayesian analysis of inspiral signals

We present in this paper a Bayesian parameter estimation method for the analysis of interferometric gravitational wave observations of an inspiral of binary compact objects using data recorded simultaneously by a network of several interferometers at different sites. We consider neutron star or black hole inspirals that are modeled to 3.5 post-Newtonian (PN) order in phase and 2.5 PN in amplitude. Inference is facilitated using Markov chain Monte Carlo methods that are adapted in order to efficiently explore the particular parameter space. Examples are shown to illustrate how and what information about the different parameters can be derived from the data. This study uses simulated signals and data with noise characteristics that are assumed to be defined by the LIGO and Virgo detectors operating at their design sensitivities. Nine parameters are estimated, including those associated with the binary system, plus its location on the sky. We explain how this technique will be part of a detection pipeline for binary systems of compact objects with masses up to 20 \sunmass, including cases where the ratio of the individual masses can be extreme.Comment: Accepted for publication in Classical and Quantum Gravity, Special issue for GWDAW-1

Upper limits on the strength of periodic gravitational waves from PSR J1939+2134

The first science run of the LIGO and GEO gravitational wave detectors presented the opportunity to test methods of searching for gravitational waves from known pulsars. Here we present new direct upper limits on the strength of waves from the pulsar PSR J1939+2134 using two independent analysis methods, one in the frequency domain using frequentist statistics and one in the time domain using Bayesian inference. Both methods show that the strain amplitude at Earth from this pulsar is less than a few times $10^{-22}$.Comment: 7 pages, 1 figure, to appear in the Proceedings of the 5th Edoardo Amaldi Conference on Gravitational Waves, Tirrenia, Pisa, Italy, 6-11 July 200

Improving the sensitivity to gravitational-wave sources by modifying the input-output optics of advanced interferometers

We study frequency dependent (FD) input-output schemes for signal-recycling interferometers, the baseline design of Advanced LIGO and the current configuration of GEO 600. Complementary to a recent proposal by Harms et al. to use FD input squeezing and ordinary homodyne detection, we explore a scheme which uses ordinary squeezed vacuum, but FD readout. Both schemes, which are sub-optimal among all possible input-output schemes, provide a global noise suppression by the power squeeze factor, while being realizable by using detuned Fabry-Perot cavities as input/output filters. At high frequencies, the two schemes are shown to be equivalent, while at low frequencies our scheme gives better performance than that of Harms et al., and is nearly fully optimal. We then study the sensitivity improvement achievable by these schemes in Advanced LIGO era (with 30-m filter cavities and current estimates of filter-mirror losses and thermal noise), for neutron star binary inspirals, and for narrowband GW sources such as low-mass X-ray binaries and known radio pulsars. Optical losses are shown to be a major obstacle for the actual implementation of these techniques in Advanced LIGO. On time scales of third-generation interferometers, like EURO/LIGO-III (~2012), with kilometer-scale filter cavities, a signal-recycling interferometer with the FD readout scheme explored in this paper can have performances comparable to existing proposals. [abridged]Comment: Figs. 9 and 12 corrected; Appendix added for narrowband data analysi

Nuclear astrophysics: the unfinished quest for the origin of the elements

Half a century has passed since the foundation of nuclear astrophysics. Since then, this discipline has reached its maturity. Today, nuclear astrophysics constitutes a multidisciplinary crucible of knowledge that combines the achievements in theoretical astrophysics, observational astronomy, cosmochemistry and nuclear physics. New tools and developments have revolutionized our understanding of the origin of the elements: supercomputers have provided astrophysicists with the required computational capabilities to study the evolution of stars in a multidimensional framework; the emergence of high-energy astrophysics with space-borne observatories has opened new windows to observe the Universe, from a novel panchromatic perspective; cosmochemists have isolated tiny pieces of stardust embedded in primitive meteorites, giving clues on the processes operating in stars as well as on the way matter condenses to form solids; and nuclear physicists have measured reactions near stellar energies, through the combined efforts using stable and radioactive ion beam facilities. This review provides comprehensive insight into the nuclear history of the Universe and related topics: starting from the Big Bang, when the ashes from the primordial explosion were transformed to hydrogen, helium, and few trace elements, to the rich variety of nucleosynthesis mechanisms and sites in the Universe. Particular attention is paid to the hydrostatic processes governing the evolution of low-mass stars, red giants and asymptotic giant-branch stars, as well as to the explosive nucleosynthesis occurring in core-collapse and thermonuclear supernovae, gamma-ray bursts, classical novae, X-ray bursts, superbursts, and stellar mergers.Comment: Invited Review. Accepted for publication in "Reports on Progress in Physics" (version with low-resolution figures

Localization and broadband follow-up of the gravitational-wave transient GW 150914

A gravitational-wave (GW) transient was identified in data recorded by the Advanced Laser Interferometer Gravitational-wave Observatory (LIGO) detectors on 2015 September 14. The event, initially designated G184098 and later given the name GW150914, is described in detail elsewhere. By prior arrangement, preliminary estimates of the time, significance, and sky location of the event were shared with 63 teams of observers covering radio, optical, near-infrared, X-ray, and gamma-ray wavelengths with ground- and space-based facilities. In this Letter we describe the low-latency analysis of the GW data and present the sky localization of the first observed compact binary merger. We summarize the follow-up observations reported by 25 teams via private Gamma-ray Coordinates Network circulars, giving an overview of the participating facilities, the GW sky localization coverage, the timeline, and depth of the observations. As this event turned out to be a binary black hole merger, there is little expectation of a detectable electromagnetic (EM) signature. Nevertheless, this first broadband campaign to search for a counterpart of an Advanced LIGO source represents a milestone and highlights the broad capabilities of the transient astronomy community and the observing strategies that have been developed to pursue neutron star binary merger events. Detailed investigations of the EM data and results of the EM follow-up campaign are being disseminated in papers by the individual teams

Search for gravitational-lensing signatures in the full third observing run of the LIGO-Virgo network

Gravitational lensing by massive objects along the line of sight to the source causes distortions of gravitational wave-signals; such distortions may reveal information about fundamental physics, cosmology and astrophysics. In this work, we have extended the search for lensing signatures to all binary black hole events from the third observing run of the LIGO--Virgo network. We search for repeated signals from strong lensing by 1) performing targeted searches for subthreshold signals, 2) calculating the degree of overlap amongst the intrinsic parameters and sky location of pairs of signals, 3) comparing the similarities of the spectrograms amongst pairs of signals, and 4) performing dual-signal Bayesian analysis that takes into account selection effects and astrophysical knowledge. We also search for distortions to the gravitational waveform caused by 1) frequency-independent phase shifts in strongly lensed images, and 2) frequency-dependent modulation of the amplitude and phase due to point masses. None of these searches yields significant evidence for lensing. Finally, we use the non-detection of gravitational-wave lensing to constrain the lensing rate based on the latest merger-rate estimates and the fraction of dark matter composed of compact objects

Search for eccentric black hole coalescences during the third observing run of LIGO and Virgo

Despite the growing number of confident binary black hole coalescences observed through gravitational waves so far, the astrophysical origin of these binaries remains uncertain. Orbital eccentricity is one of the clearest tracers of binary formation channels. Identifying binary eccentricity, however, remains challenging due to the limited availability of gravitational waveforms that include effects of eccentricity. Here, we present observational results for a waveform-independent search sensitive to eccentric black hole coalescences, covering the third observing run (O3) of the LIGO and Virgo detectors. We identified no new high-significance candidates beyond those that were already identified with searches focusing on quasi-circular binaries. We determine the sensitivity of our search to high-mass (total mass M>70 M⊙) binaries covering eccentricities up to 0.3 at 15 Hz orbital frequency, and use this to compare model predictions to search results. Assuming all detections are indeed quasi-circular, for our fiducial population model, we place an upper limit for the merger rate density of high-mass binaries with eccentricities 0<e≤0.3 at 0.33 Gpc−3 yr−1 at 90\% confidence level