Decoding of the light changes in eclipsing Wolf-Rayet binaries I. A non-classical approach to the solution of light curves

We present a technique to determine the orbital and physical parameters of eclipsing eccentric Wolf-Rayet + O-star binaries, where one eclipse is produced by the absorption of the O-star light by the stellar wind of the W-R star. Our method is based on the use of the empirical moments of the light curve that are integral transforms evaluated from the observed light curves. The optical depth along the line of sight and the limb darkening of the W-R star are modelled by simple mathematical functions, and we derive analytical expressions for the moments of the light curve as a function of the orbital parameters and the key parameters of the transparency and limb-darkening functions. These analytical expressions are then inverted in order to derive the values of the orbital inclination, the stellar radii, the fractional luminosities, and the parameters of the wind transparency and limb-darkening laws. The method is applied to the SMC W-R eclipsing binary HD 5980, a remarkable object that underwent an LBV-like event in August 1994. The analysis refers to the pre-outburst observational data. A synthetic light curve based on the elements derived for the system allows a quality assessment of the results obtained.Comment: Accepted for publication in Astronomy & Astrophysic

A Comparison of search templates for gravitational waves from binary inspiral

We compare the performances of the templates defined by three different types of approaches: traditional post-Newtonian templates (Taylor-approximants), ``resummed'' post-Newtonian templates assuming the adiabatic approximation and stopping before the plunge (P-approximants), and further ``resummed'' post-Newtonian templates going beyond the adiabatic approximation and incorporating the plunge with its transition from the inspiral (Effective-one-body approximants). The signal to noise ratio is significantly enhanced (mainly because of the inclusion of the plunge signal) by using these new effective-one-body templates relative to the usual post-Newtonian ones for binary masses greater than $30 M_\odot$, the most likely sources for initial laser interferometers. Independently of the question of the plunge signal, the comparison of the various templates confirms the usefulness of using resummation methods. The paper also summarizes the key elements of the construction of various templates and thus can serve as a resource for those involved in writing inspiral search software.Comment: eta-dependent tail terms corrected after related errata by Blanchet (2005

STROOPWAFEL: Simulating rare outcomes from astrophysical populations, with application to gravitational-wave sources

Gravitational-wave observations of double compact object (DCO) mergers are providing new insights into the physics of massive stars and the evolution of binary systems. Making the most of expected near-future observations for understanding stellar physics will rely on comparisons with binary population synthesis models. However, the vast majority of simulated binaries never produce DCOs, which makes calculating such populations computationally inefficient. We present an importance sampling algorithm, STROOPWAFEL, that improves the computational efficiency of population studies of rare events, by focusing the simulation around regions of the initial parameter space found to produce outputs of interest. We implement the algorithm in the binary population synthesis code COMPAS, and compare the efficiency of our implementation to the standard method of Monte Carlo sampling from the birth probability distributions. STROOPWAFEL finds $\sim$25-200 times more DCO mergers than the standard sampling method with the same simulation size, and so speeds up simulations by up to two orders of magnitude. Finding more DCO mergers automatically maps the parameter space with far higher resolution than when using the traditional sampling. This increase in efficiency also leads to a decrease of a factor $\sim$3-10 in statistical sampling uncertainty for the predictions from the simulations. This is particularly notable for the distribution functions of observable quantities such as the black hole and neutron star chirp mass distribution, including in the tails of the distribution functions where predictions using standard sampling can be dominated by sampling noise.Comment: Accepted. Data and scripts to reproduce main results is publicly available. The code for the STROOPWAFEL algorithm will be made publicly available. Early inquiries can be addressed to the lead autho

Phenomenological model for the gravitational-wave signal from precessing binary black holes with two-spin effects

The properties of compact binaries, such as masses and spins, are imprinted in the gravitational-waves they emit and can be measured using parameterised waveform models. Accurately and efficiently describing the complicated precessional dynamics of the various angular momenta of the system in these waveform models is the object of active investigation. One of the key models extensively used in the analysis of LIGO and Virgo data is the single-precessing-spin waveform model IMRPhenomPv2. In this article we present a new model IMRPhenomPv3 which includes the effects of two independent spins in the precession dynamics. Whereas IMRPhenomPv2 utilizes a single-spin frequency-dependent post-Newtonian rotation to describe precession effects, the improved model, IMRPhenomPv3, employs a double-spin rotation that is based on recent developments in the description of precessional dynamics. Besides double-spin precession, the improved model benefits from a more accurate description of precessional effects. We validate our new model against a large set of precessing numerical-relativity simulations. We find that IMRPhenomPv3 has better agreement with the inspiral portion of precessing binary-black-hole simulations and is more robust across a larger region of the parameter space than IMRPhenomPv2. As a first application we analyse, for the first time, the gravitational-wave event GW151226 with a waveform model that describes two-spin precession. Within statistical uncertainty our results are consistent with published results. IMRPhenomPv3 will allow studies of the measurability of individual spins of binary black holes using GWs and can be used as a foundation upon which to build further improvements, such as modeling precession through merger, extending to higher multipoles, and including tidal effects.Comment: 15 pages, 5 figure

Range Analysis of Binaries with Minimal Effort

COTS components are ubiquitous in military, industrial and governmental systems. However, the bene?fits of reduced development and maintainance costs are compromised by security concerns. Since source code is unavailable, security audits necessarily occur at the binary level. Push-button formal method techniques, such as model checking and abstract interpretation, can support this process by, among other things, inferring ranges of values for registers. Ranges aid the security engineer in checking for vulnerabilities that relate, for example, to integer wrapping, uninitialised variables and bu?er over ows. Yet the lack of structure in binaries limits the e?ffectiveness of classical range analyses based on widening. This paper thus contributes a simple but novel range analysis, formulated in terms of linear programming, which calculates ranges without manual intervention