Source author record

Marc van der Sluys

Marc van der Sluys appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

11works
5topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

11 published item(s)

preprint2015arXiv

Robust parameter estimation for compact binaries with ground-based gravitational-wave observations using the LALInference software library

The Advanced LIGO and Advanced Virgo gravitational wave (GW) detectors will begin operation in the coming years, with compact binary coalescence events a likely source for the first detections. The gravitational waveforms emitted directly encode information about the sources, including the masses and spins of the compact objects. Recovering the physical parameters of the sources from the GW observations is a key analysis task. This work describes the LALInference software library for Bayesian parameter estimation of compact binary signals, which builds on several previous methods to provide a well-tested toolkit which has already been used for several studies. We show that our implementation is able to correctly recover the parameters of compact binary signals from simulated data from the advanced GW detectors. We demonstrate this with a detailed comparison on three compact binary systems: a binary neutron star, a neutron star black hole binary and a binary black hole, where we show a cross-comparison of results obtained using three independent sampling algorithms. These systems were analysed with non-spinning, aligned spin and generic spin configurations respectively, showing that consistent results can be obtained even with the full 15-dimensional parameter space of the generic spin configurations. We also demonstrate statistically that the Bayesian credible intervals we recover correspond to frequentist confidence intervals under correct prior assumptions by analysing a set of 100 signals drawn from the prior. We discuss the computational cost of these algorithms, and describe the general and problem-specific sampling techniques we have used to improve the efficiency of sampling the compact binary coalescence parameter space.

preprint2013arXiv

Reconstructing the sky location of gravitational-wave detected compact binary systems: methodology for testing and comparison

The problem of reconstructing the sky position of compact binary coalescences detected via gravitational waves is a central one for future observations with the ground-based network of gravitational-wave laser interferometers, such as Advanced LIGO and Advanced Virgo. Different techniques for sky localisation have been independently developed. They can be divided in two broad categories: fully coherent Bayesian techniques, which are high-latency and aimed at in-depth studies of all the parameters of a source, including sky position, and "triangulation-based" techniques, which exploit the data products from the search stage of the analysis to provide an almost real-time approximation of the posterior probability density function of the sky location of a detection candidate. These techniques have previously been applied to data collected during the last science runs of gravitational-wave detectors operating in the so-called initial configuration. Here, we develop and analyse methods for assessing the self-consistency of parameter estimation methods and carrying out fair comparisons between different algorithms, addressing issues of efficiency and optimality. These methods are general, and can be applied to parameter estimation problems other than sky localisation. We apply these methods to two existing sky localisation techniques representing the two above-mentioned categories, using a set of simulated inspiral-only signals from compact binary systems with total mass $\le 20\,M_\odot$ and non-spinning components. We compare the relative advantages and costs of the two techniques and show that sky location uncertainties are on average a factor $\approx 20$ smaller for fully coherent techniques than for the specific variant of the "triangulation-based" technique used during the last science runs, at the expense of a factor $\approx 1000$ longer processing time.

preprint2013arXiv

Using electromagnetic observations to aid gravitational-wave parameter estimation of compact binaries observed with LISA II: The effect of knowing the sky position

In this follow-up paper, we continue our study of the effect of using knowledge from electromagnetic observations in the gravitational wave (GW) data analysis of Galactic binaries that are predicted to be observed by the new \textit{Laser Interferometer Space Antenna} in the low-frequency range, $10^{-4} \:\mathrm{Hz}<f<1 \:\mathrm{Hz}$. In the first paper, we have shown that the strong correlation between amplitude and inclination can be used for mildly inclined binaries to improve the uncertainty in amplitude, and that this correlation depends on the inclination of the system. In this paper we investigate the overall effect of the other orientation parameters, namely the sky position and the polarisation angle. We find that after the inclination, the ecliptic latitude of the source has the strongest effect in determining the GW parameter uncertainties. We ascertain that the strong correlation we found previously, only depends on the inclination of the source and not on the other orientation parameters. We find that knowing the sky position of the source from electromagnetic data can reduce the GW parameter uncertainty up to a factor of $\sim 2$, depending on the inclination and the ecliptic latitude of the system. Knowing the sky position and inclination can reduce the uncertainty in amplitude by a factor larger than 40. We also find that unphysical errors in the inclinations, which we found when using the Fisher matrix, can affect the corresponding uncertainties in the amplitudes, which need to be corrected.

preprint2012arXiv

Using electromagnetic observations to aid gravitational-wave parameter estimation of compact binaries observed with LISA

We present a first-stage study of the effect of using knowledge from electromagnetic (EM) observations in the gravitational wave (GW) data analysis of Galactic binaries that are predicted to be observed by the new Laser Interferometer Space Antenna in the low-frequency range, $10^{-4} \mathrm{Hz}<f<1 \mathrm{Hz}$. In particular, we examine the extent to which the accuracy of GW parameter estimation improves if we use available information from EM data. We do this by investigating whether correlations exist between the GW parameters that describe these binaries and whether some of these parameters are also available from EM observations. We used verification binaries, which are known as the guaranteed sources for eLISA and will test the functioning of the instrument. We find that of the seven parameters that characterise such a binary, only a few are correlated. The most useful result is the strong correlation between amplitude and inclination, which can be used to constrain the parameter uncertainty in amplitude by making use of the constraint of inclination from EM measurements. The improvement can be up to a factor of $\sim6.5$, but depends on the signal-to-noise ratio of the source data. Moreover, we find that this strong correlation depends on the inclination. For mildly face-on binaries ($ι\lesssim 45^{\circ}$), EM data on inclination can improve the estimate of the GW amplitude by a significant factor. However, for edge-on binaries ($ι\sim 90^{\circ}$), the inclination can be determined accurately from GW data alone, thus GW data can be used to select systems that will likely be eclipsing binaries for EM follow-up.

preprint2011arXiv

Estimating parameters of coalescing compact binaries with a detector network including LIGO Australia

One of the goals of gravitational-wave astronomy is simultaneous detection of gravitational-wave signals from merging compact-object binaries and the electromagnetic transients from these mergers. With the next generation of advanced ground-based gravitational wave detectors under construction, we examine the benefits of the proposed extension of the detector network to include a fourth site in Australia in addition to the network of Hanford, Livingston and Cascina sites. Using Bayesian parameter-estimation analyses of simulated gravitational-wave signals from a range of coalescing-binary locations and orientations, we study the improvement in parameter estimation. We find that an Australian detector can break degeneracies in several parameters; in particular, the localization of the source on the sky is improved by a factor of ~4, with more modest improvements in distance and binary inclination estimates. This enhanced ability to localize sources on the sky will be crucial in any search for electromagnetic counterparts to detected gravitational-wave signals.

preprint2011arXiv

Gravitational waves from compact binaries

In this review, I give a summary of the history of our understanding of gravitational waves and how compact binaries were used to transform their status from mathematical artefact to physical reality. I also describe the types of compact (stellar) binaries that LISA will observe as soon as it is switched on. Finally, the status and near future of LIGO, Virgo and GEO are discussed, as well as the expected detection rates for the Advanced detectors, and the accuracies with which binary parameters can be determined when BH/NS inspirals are detected.

preprint2010arXiv

Analytical expressions for the envelope binding energy of giants as a function of basic stellar parameters

The common-envelope (CE) phase is an important stage in the evolution of binary stellar populations. The most common way to compute the change in orbital period during a CE is to relate the binding energy of the envelope of the Roche-lobe filling giant to the change in orbital energy. Especially in population-synthesis codes, where the evolution of millions of stars must be computed and detailed evolutionary models are too expensive computationally, simple approximations are made for the envelope binding energy. In this study, we present accurate analytic prescriptions based on detailed stellar-evolution models that provide the envelope binding energy for giants with metallicities between Z = 10-4 and Z = 0.03 and masses between 0.8 Msun and 100 Msun, as a function of the metallicity, mass, radius and evolutionary phase of the star. Our results are also presented in the form of electronic data tables and Fortran routines that use them. We find that the accuracy of our fits is better than 15% for 90% of our model data points in all cases, and better than 10% for 90% of our data points in all cases except the asymptotic giant branches for three of the six metallicities we consider. For very massive stars (M > 50 Msun), when stars lose more than ~20% of their initial mass due to stellar winds, our fits do not describe the models as accurately. Our results are more widely applicable - covering wider ranges of metallicity and mass - and are of higher accuracy than those of previous studies.

preprint2010arXiv

Masses and envelope binding energies of primary stars at the onset of a common envelope

We present basic properties of primary stars that initiate a common envelope (CE) in a binary, while on the giant branch. We use the population-synthesis code described in Politano et al. (2010) and follow the evolution of a population of binary stars up to the point where the primary fills its Roche lobe and initiates a CE. We then collect the properties of each system, in particular the donor mass and the binding energy of the donor's envelope, which are important for the treatment of a CE. We find that for most CEs, the donor mass is sufficiently low to define the core-envelope boundary reasonably well. We compute the envelope-structure parameter {λ_\mathrm{env}} from the binding energy and compare its distribution to typical assumptions that are made in population-synthesis codes. We conclude that {λ_\mathrm{env}} varies appreciably and that the assumption of a constant value for this parameter results in typical errors of 20--50%. In addition, such an assumption may well result in the implicit assumption of unintended and/or unphysical values for the CE parameter {α_\mathrm{CE}}. Finally, we discuss accurate existing analytic fits for the envelope binding energy, which make these oversimplified assumptions for {λ_\mathrm{env}}, and the use of {λ_\mathrm{env}} in general, unnecessary.

preprint2010arXiv

Population Synthesis of Common Envelope Mergers: I. Giant Stars with Stellar or Substellar Companions

Using a population synthesis technique, we have calculated detailed models of the present-day field population of objects that have resulted from the merger of a giant primary and a main-sequence or brown dwarf secondary during common-envelope evolution. We used a grid of 116 stellar and 32 low-mass/brown dwarf models, a crude model of the merger process, and followed the angular momentum evolution of the binary orbit and the primary's rotation prior to merger, as well as the merged object's rotation after the merger. We find that present-day merged objects that are observable as giant stars or core-helium burning stars in our model population constitute between 0.24% and 0.33% of the initial population of ZAMS binaries, depending upon the input parameters chosen. The median projected rotational velocity of these merged objects is ~16 km/sec, an order of magnitude higher than the median projected rotational velocity in a model population of normal single stars calculated using the same stellar models and initial mass function. The masses of the merged objects are typically less than ~2 solar masses, with a median mass of 1.28 solar masses, which is slightly more than, but not significantly different from, their normal single star counterparts. The luminosities in our merged object population range from ~10-100 solar luminosities, with a strong peak in the luminosity distribution at ~60 solar luminosities, since the majority of the merged objects (57%) lie on the horizontal branch at the present epoch. The results of our population synthesis study are discussed in terms of possible observational counterparts either directly involving the high rotational velocity of the merger product or indirectly, via the effect of rotation on envelope abundances and on the amount and distribution of circumstellar matter.

preprint2009arXiv

Status of NINJA: the Numerical INJection Analysis project

The 2008 NRDA conference introduced the Numerical INJection Analysis project (NINJA), a new collaborative effort between the numerical relativity community and the data analysis community. NINJA focuses on modeling and searching for gravitational wave signatures from the coalescence of binary system of compact objects. We review the scope of this collaboration and the components of the first NINJA project, where numerical relativity groups shared waveforms and data analysis teams applied various techniques to detect them when embedded in colored Gaussian noise.

preprint2009arXiv

Testing gravitational-wave searches with numerical relativity waveforms: Results from the first Numerical INJection Analysis (NINJA) project

The Numerical INJection Analysis (NINJA) project is a collaborative effort between members of the numerical relativity and gravitational-wave data analysis communities. The purpose of NINJA is to study the sensitivity of existing gravitational-wave search algorithms using numerically generated waveforms and to foster closer collaboration between the numerical relativity and data analysis communities. We describe the results of the first NINJA analysis which focused on gravitational waveforms from binary black hole coalescence. Ten numerical relativity groups contributed numerical data which were used to generate a set of gravitational-wave signals. These signals were injected into a simulated data set, designed to mimic the response of the Initial LIGO and Virgo gravitational-wave detectors. Nine groups analysed this data using search and parameter-estimation pipelines. Matched filter algorithms, un-modelled-burst searches and Bayesian parameter-estimation and model-selection algorithms were applied to the data. We report the efficiency of these search methods in detecting the numerical waveforms and measuring their parameters. We describe preliminary comparisons between the different search methods and suggest improvements for future NINJA analyses.